氯气泄漏扩散计算模拟

氯气泄漏扩散计算模拟

（1）氯气泄漏扩散过程

氯气管道、阀门、压力表或液氯钢瓶的材质、焊缝以及腐蚀造成的强度下降等设备缺陷原因都可能引起容器破裂从而引发氯气泄漏，使氯气迅速扩展。

氯气属于重气，它泄漏时的扩散机理与一般烟囱热污染性烟羽的扩散完全不同，它在地表的释放可分为三个过程。

初始阶段：气云刚形成的阶段，主导其运动的作用力为释放的惯性及外界的平均风速。

重力扩展阶段：初始的动量消失后，主导作用为重力及外界的湍流扰动，由于重力使气团下降到地表拓展范围，而且稀释作用主要靠大气湍流以及气云下降引起的湍流卷增作用引入外界空气。

被动扩散阶段：当气云经一段时间混合稀释后，其密度和温度或浓度逐渐接近外界空气，主要靠大气的湍流。

（2）氯气危害浓度

根据有关资料，氯气的半数致死浓度，即引起实验染毒动物半数死亡的毒物浓度为850mg/m3。氯的最高允许浓度为1mg/m3，即在一个工作日内任何时间都不应超过的浓度。选取浓度850mg/m3（283ppm）、300mg/m3（100ppm）、90mg/m3（30ppm）、30－90mg/m3（10－30ppm）、15mg/m3（5ppm）、5mg/m3（1.7ppm）和1mg/m3（0.3ppm）7种浓度分别计算氯气泄漏的影响范围。

浓度

伤害分区危害程度Mg/m3ppm

850283深吸入少量可能危及生命

300100致死区（A区）可能造成致命性损害

9030重伤区（B区）引起剧咳

（3）液氯泄漏量计算

流体泄漏速度可用流体力学的柏努利方程计算。其泄漏速度为：

()

gh p p A Cd Q 220+-••=ρ

ρ

（7-1）

式中Q 0——液体泄漏速度，kg/s ； Cd ——液体泄漏系数； A ——裂口面积，m 2; ρ——泄漏液体密度，kg/m 3; P ——容器内介质压力，Pa ； P 0——环境压力，Pa; g ——重力加速度，9.8m/s 2; h ——裂口之上液位高度，m 。

对于非常压下的液体泄漏速度，主要取决于内介质压力与环境压力之差和液位高低。

（4）液氯容器破裂形成毒害区半径估算

液化介质液氯在容器破裂时会发生蒸气扩散，从而造成大面积的毒害区域，其毒害区按如下方法估算。

设液氯质量为W （单位：kg ）,容器破裂前器内温度为t （单位：℃），液体介质比热为C[单位：KJ/（kg •℃）]当容器破裂时，器内压力降至大气压，处于过热状态的液化气温度迅速降至标准沸点t 0（单位：℃），此时全部液体所吸收的热量为：

()0t t C W Q -••= （7-2）

设这些热量全部用于器内液体的蒸发，如它的汽化热为q （单位：kJ/kg ），则其蒸发量：

()q

t t C W q Q W 0'-•==

（7-3）

如介质的分子量为M ，则在沸点下蒸发蒸气的体积Vg （单位：m 3）为：

()273

2734.222732734.220

00t q M t t C W t M W Vg +••-••=+•=

（7-4）

假设这些带有氯气的空气以半球形向地面扩散,则可求出该氯气扩散半径为:

30

3

00994.2/3

421/C Vg C Vg R =

⨯=

π （7-5） 式中R ——有毒气体的半径，m ； Vg ——有毒介质的蒸气体积，m 3； C O ——有毒介质在空气中的危险浓度，%; （5）气体连续扩散高斯分布

采用世界银行提供的中等密度云连续扩散模型（高斯扩散模型），分析时，根据分析源周边的实际情况，考虑地面有效粗糙度并进行修正。连续排放：

()⎥⎥⎦

⎤

⎢⎢⎣⎡⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+-=222221exp ,,z y z y Z y u Q

z y x c σσσπσ （7-6）

式中c （ x,y,z ）——连续排放时，给定地点的污染物浓度，mg/m 3；

Q ——连续排放的物料流量，mg/s ； u ——平均风速，m/s;

x ——下风向距离，m; y ——横风向距离，m ； z ——离地面的距离，m;

σx ,σy ,σz ——x,y,z 方向扩散参数，按下表选择：

表2 不同大气稳定度下的扩散参数

有效粗糙度Z 0＞0.1m 的粗糙地形扩散参数按如下计算（地面有效粗糙度长度按下表选取）：

y yo y f σσ=

z z z f 0σσ=

(

)0001Z Z f y α+=

()()()x g f z Z x e d x c b Z x f ln 01

0000000ln ln ,--+-=

上式中，大气稳定度为D ，各参数按下表选择：

表4 不同大气稳定度下的系数值

（6）气体瞬时扩散高斯分布

采用世界银行提供的中等密度云瞬时扩散模型（高斯扩散模型）。

()()()⎥⎥⎦

⎤

⎢⎢⎣⎡⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧++--=22222

22/321exp 22,,,z y x z y x Z y ut x Q

t z y x c σσσμσσσπ （6-7） 式中c （ x,y,z,t ）——瞬时排放时，给定地点（ x,y,z ）和时间t 的污染物浓度，mg/m 3；

Q ——瞬时排放的物料流量，mg/s ； t ——瞬时排放时，污染物运行的时间，s ； u ——平均风速，m/s ； x ——下风向距离，m ； y ——横风向距离，m ； z ——离地面的距离，m ；

σx ,σy ,σz ——x,y,z 方向扩散参数，地面有效粗糙度修正等因素的修正方法同上。瞬时源时，给定取样时间t （ s ），则按下式修正：

()()2.0min 10600/t y y σσ=