个人风险、社会风险分析

个体风险计算

步骤如下:

a)选择一个泄漏场景(LOC)，确定每个LOC的年失效频率fS。

b)选择一种天气等级M和该天气等级下的一种风向，给出天气等级M和风向同时出现的联合概率

10-6)，

假设:液氨储罐以裂口面积A(单位为m2)连续向外界泄露液态氨，泄露开始时裂口上方液体高度hL，(单位为m)，泄露出来的液氨瞬间气化成氨气烟羽，泄漏时储罐内液体压力P(单位为Pa)，外界环境压力P0(单位为Pa)，全年主导风向北风风速为u(单位m/s)，液氨密度ρ(单位为Kg/m3)，泄漏时间t(单位为s)。以液氨储罐泄露点为坐标原点O，下风向X轴，以方向东为Y轴，以垂

直地面向上为Z轴，建立坐标系(XY平面垂直地面高度1米)。

液氨经储罐上的孔流出，瞬时质量流率为Q，单位为Kg/s:

(1)

式中:C0——液体泄漏系数;g——重力加速度， m/s2;

连续排放时，形成稳定的流场后，位于地面1m高处的连续稳态源的烟羽，在给定地点的污染物的浓度，单位为kg/m3，

……(2)

式中:，——侧风向和垂直风向的扩散系数，单位为m;

吸入氨5~%，则C取500×10-6;t取最短时间5min;计算Y=-+.071ln(5002*5)=

对应死亡百分率P死亡百分率=99%

——个体风险计算时的死亡概率; P个体风险

个体风险——个体风险计算时的死亡概率修正因子;取值为1;

按泄漏量为3吨时，距泄漏点150米处老年房为研究对象

P个体风险=×10×40%=×10,-5-6-7

本项目个人风险不可容许。

社会风险计算

P社会风险社会风险

P社会风险——社会风险计算时的死亡概率;

社会风险——社会风险计算时的死亡概率修正因子;，

室外为1。;

查可容许社会风险标准(F-N)曲线，处于不可容许区。

则本项目重大危险源社会风险不可容许。

有毒物个体风险和社会风险的死亡概率计算

在确定个体风险和社会风险时需计算网格单元点的死亡概率。假设有毒气云的宽度有限，仅仅覆盖一个风面。

在网格单元位置的有毒气云用具有均衡死亡概率的等效气云代替，网格单元的死亡概率Pd为等效气云内的死亡概率乘以网格单元位于等效气云内部的概率。

计算步骤如下:

(1)计算网格单元与释放源的距离R=150米

(2)计算全部时间t内，高度h,1 m的中心线的浓度C(R，t)

(3)计算距离R，参考高度h,1 m处烟羽中心的死亡概率Pcl,个体风险(R)

Pcl,社会风险(R)的计算公式为:

…… () Pcl,社会风险社会风险社会风险

(4)计算离烟羽中心线距离为y处的死亡概率P(R，y)。离烟羽中心线越远，

其死亡概率越小。

(5)计算距离R处的气云有效宽度(ECW)。在距离R的有毒气云用一个等效气云代替，假设等效气云有一常数的死亡概率，其值等于有毒气云中心线的死亡概率。通过有毒气云内的死亡数等于等效气云的死亡数，来确定等效气云的宽度。等效气云宽度可以使用概率积分来计算。概率积分PI(R)定义为沿着垂直于烟羽中心线的轴，对死亡概率进行积分:

()

式中:

y——离烟羽中心线的距离，单位为m;

P(R，y)——位置(R，y)的死亡概率。

在数值计算时，积分边界可取1,致死率的距离。

等效气云宽度ECW定义为:

P)…… ()

cl(R

6. 计算网格单元被等效气云覆盖的概率Pci，如果网格单元在风面内，可近

似利用下式求得:

PECW(R)

()

如果网格单元不在风面内，。参数nws等于风面数。

7. 计算网格单元死亡概率P个体风险

P个体风险，个体风险()

8. 计算网格单元死亡概率P社会风险

P社会风险，社会风险()

液体经储罐上的孔流出

瞬时质量流率为:

()

式中:

Qm——质量流率，单位为Kg/s;

P——储罐泄漏物质在大气中的扩散 大气稳定度确定

大气稳定度确定通常采用Pasquill分类方法,大气稳定度分为A、B、C、D、E和F六类，A类表示气象条件极不稳定，B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定，D类表示气象条件的稳定性在稳定和不稳定之间，E类表示气象条件弱稳定，F类表示气象条件中