个人风险、社会风险分析11.4.1个体风险计算步骤如下：R选择一个泄漏场景(LOC),确定每个LOC的年失效频率fSo选择一种天气等级M和该天气等级下的一种风向，给出天气等级M和风向同时出现的联合概率PM PPM P天气等级M和风向同时出现的联合概率表计算在特定的L0C、天气等级M、风向及点火事件(可燃物)i条件下网格单元 上的死亡概率P个体风险毒性暴露下死亡概率值可按下式计算：Pr abln(Cn t)(7)式中：a, b,n——描述物质毒性的常量，见附录K; C——浓度，单位为mg/m3 ;t暴露于毒物中的时间，单位为min,最大值为30mino其中:氨的常数 A二 - 9. 82;B=0.71;n=2.0;假设：（1）北风、小口泄漏液氨储罐以裂口面积A（单位为m2）连续向外界泄漏液态氨，泄漏出来的液氨瞬 间气化成氨气烟羽，泄漏时储罐内液体压力P（单位为Pa）,外界环境压力P0（单位 为PR,全年主导风向北风风速为u（单位m/s）,液氨密度P （单位为Kg/m3）,泄 漏时间t （单位为s）。以液氨储罐泄漏点为坐标原点0,下风向X轴，以方向东为Y 轴，以垂直地面向上为Z轴，建立坐标系（XY平面垂直地面高度1米）。

XZ0Y液氨经储罐上的孔流出，瞬时质量流率为Q,单位为Kg./s:叩0 Q AC02..（1）式中:CO——液体泄漏系数；连续排放时，形成稳定的流场后，位于地面2m高处的连续稳态源的烟羽，在xyO Cx?y,0的污染物的浓度，单位为kg/m3,给定地点222zHzHQlyll rr式中：y z侧风向和垂直风向的扩散系数，单位为ID；d）计算（LOC、M、、i）条件下对网格单元个体风险的贡献：1RS?M? ?i fS PM P P个体风险 （ii）e）对所有的天气等级和风向，重复b）-d）步的计算;对所有的LOC,重复a）-d） 步的计算，贝9网格点处的个人风险由下式计算：rR. S FRMi SAL ,i....（12）个人风险pP个体风险个体风险式中:P——死亡概率，为年失效频率F年失效频率与死亡口分率P死亡口分率 的乘积；死亡百分率P死亡百分率：在一定时间接触一定浓度毒物所造成影响的概率值Y ABln Cn t0其中:氨的常数 A二-9. 82;B=0.71;n=2.0;C为接触毒物的浓度（单位10-6）,假设:液氨储罐以裂口面积A（单位为m2）连续向外界泄露液态氨，泄露开始时 裂口上方液体高度hL,（单位为m）,泄露出来的液氨瞬间气化成氨气烟羽，泄漏时 储罐内液体压力P（单位为Pa）,外界环境压力P0 （单位为Pa）,全年主导风向北风 风速为U（单位m/s）,液氨密度P （单位为Kg/m3）,泄漏时间t （单位为s）。

以液氨 储罐泄露点为坐标原点0,下风向X轴，以方向东为Y轴，以垂直地面向上为Z轴，建立坐标系（XY平面垂直地面高度1米）。XZ0Y液氨经储罐上的孔流出，瞬时质量流率为Q,单位为Kg./s:叩0 Q AC02 ghT,.....（1）式中:CO——液体泄漏系数;g——重力加速度，9.8 m/s2;连续排放时，形成稳定的流场后，位于地面lm高处的连续稳态源的烟羽，在给定地点的污染物的浓度，单位为kg/m3,X$1 Cx,y,l?22QiwCi!呻即2 y 2ii2 2 yz……⑵y z式中：，一一侧风向和垂直风向的扩散系数，单位为m;吸入氨5~10min致死的浓度为0. 05%,则C取500 X 10-6 ;t取最短时间omin;计算 Y二-9. 82+. 071 In （5002*5）=0. 148对应死亡白分率P死亡口分率二99%——个体风险计算时的死亡概率；P个体风险个体风险一一个体风险计算时的死亡概率修正因子;取值为1;按泄漏量为3吨时，距泄漏点150米处老年房为研究对象310P 个体风险二2. 1X10X40%二8. 4X 10, -5-6-7本项目个人风险不可容许。4. 3. 2社会风险计算PP社会风险社会风险P社会风险一一社会风险讣算时的死亡概率；社会风险一一社会风险计算时的死亡概率修正因子;室内风险修正因子为0.1, 室外为1。

考虑夜间人员均处于室内惜况取0.1;查可容许社会风险标准(F-N)曲线，处于不可容许区。则本项U重大危险源社会风险不可容许。L.2有毒物个体风险和社会风险的死亡概率讣算在确定个体风险和社会风险时需计算网格单元点的死亡概率。假设有毒气云的 宽度有限，仅仅覆盖一个风面。在网格单元位置的有毒气云用具有均衡死亡概率的等效气云代替，网格单元的 死亡概率Pd为等效气云内的死亡概率乘以网格单元位于等效气云内部的概率。计算步骤如下：计算网格单元与释放源的距离R二150米计算全部时间t内，高度h,l in的中心线的浓度C(R, t)计算距离R,参考高度h,l m处烟羽中心的死亡概率Pci,个体风险(R)Pci,社会风险(R)的计算公式为：Pci, ?in fpop,inPcl,?out fpop,out(L. 1) Pel,社会风险社会风险社会风险计算离烟羽中心线距离为y处的死亡概率P(R, y)o离烟羽中心线越远， 其死亡概率越小。计算距离R处的气云有效宽度(ECW)。在距离R的有毒气云用一个等效气云 代替，假设等效气云有一常数的死亡概率，其值等于有毒气云中心线的死亡概率。 通过有毒气云内的死亡数等于等效气云的死亡数，来确定等效气云的宽度。

等效气 云宽度可以使用概率积分来讣算。概率积分PI(R)定义为沿着垂直于烟羽中心线的 轴，对死亡概率进行积分：PI(R)(Rj)dy……(L. 2)P式中：y——离烟羽中心线的距离，单位为m;P(R, y)——位置(R, y)的死亡概率。在数值计算时，积分边界可取1,致死率的距离。等效气云宽度ECW定义为：CW(R) PI(R)EP)……(L. 3)cl(R计算网格单元被等效气云覆盖的概率Pci,如果网格单元在风面内，可近似利用下式求得：PECW(R)ci(R) nws2 R(L. 4)Pci(R) 0如果网格单元不在风面内，。参数nws等于风面数。计算网格单元死亡概率P个体风险Pci PciP个体风险，个体风险(L. 5)计算网格单元死亡概率P社会风险Pel PciP社会风险，社会风险(L. 6)液体经储罐上的孔流出瞬时质量流率为：ppOQmAC02 一 ghL(E. 2)式中：Qm——质量流率，单位为Kg/s;P——储罐泄漏物质在大气中的扩散E. 3.1大气稳定度确定大气稳定度确定通常采用Pasquil 1分类方法，大气稳定度分为A、B、C、D、E 和F六类，A类表示气象条件极不稳定，B类表示气象条件中等程度不稳定，C类 表示气象条件弱不稳定，D类表示气象条件的稳定性在稳定和不稳定之间，E类表 示气象条件弱稳定，F类表示气象条件中等程度稳定。

大气稳定度的具体分类见表 E. 4 和表 E. 5。表E. 4 Pasquill大气稳定度的确定E. 3. 2 Pasquill-Gifford 模型扩散方程辻ford模型扩散方程l)Pasquill-GR位于地面Hr高处的连续稳态源的烟羽1 yQ?Cx?y?z exp2 y zu 2y式中:1 zHrexpz21 zHr.exp2 z(E. 31)Cx,y?z x,y,z——连续排放时，形成稳定的流场后，给定地点的污染物的浓度,单位为kg/m3;Q——连续排放的物料质量流量，单位为kg/s;U风速，单位为m/s;y z,——侧风向和垂直风向的扩散系数，单位为m;X——下风向距离，单位为ID；y——侧风向距离，单位为ID； z——垂直风向距离，单位为nu 2）Pasquill- Gifford模型扩散系数扩散系数确定见下表E. 4和E. 5：表E. 6推荐的烟羽扩散Pasquill-Gifford模型扩散系数方程假设:液氨储罐以裂口面积A（单位为m2）连续向外界泄露液态氨，泄露开始时 裂口上方液体高度hL,（单位为m）,泄露出来的液氨瞬间气化成氨气烟羽，