氯气泄漏扩散计算模拟.docx

1、氯气泄漏扩散计算模拟氯气泄漏扩散计算模拟(1)氯气泄漏扩散过程氯气管道、阀门、压力表或液氯钢瓶的材质、焊缝以及腐蚀造成 的强度下降等设备缺陷原因都可能引起容器破裂从而引发氯气泄漏， 使氯气迅速扩展。氯气属于重气，它泄漏时的扩散机理与一般烟囱热污染性烟羽的 扩散完全不同，它在地表的释放可分为三个过程。初始阶段：气云刚形成的阶段，主导其运动的作用力为释放的 惯性及外界的平均风速。重力扩展阶段：初始的动量消失后，主导作用为重力及外界的 湍流扰动，由于重力使气团下降到地表拓展范围，而且稀释作用 主要靠大气湍流以及气云下降引起的湍流卷增作用引入外界空 气。被动扩散阶段：当气云经一段时间混合稀释后，其密度

2、和温度 或浓度逐渐接近外界空气，主要靠大气的湍流。(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/m(0.3ppm)7种浓度分别计算氯气泄漏的影响范围 本模拟分析的计算依据如下：浓度伤害分区危害程度Mg/m 3ppm850283深吸入少量

3、可能危及生命300100致死区(A区)可能造成致命性损害9030重伤区(B区)引起剧咳30 901030轻伤区（C区）刺激咽喉，引起剧咳155吸入反应区（D区）刺激咽喉62有明显气味，刺激眼、鼻10. 03嗅觉浓度（3）液氯泄漏量计算流体泄漏速度可用流体力学的柏努利方程计算。其泄漏速度为:Cd ?A?2p P0 2gh(7-1)式中Q0 液体泄漏速度，kg/s ；Cd 液体泄漏系数；A 裂口面积，m2;P 泄漏液体密度，kg/m 3;P容器内介质压力，Pa;Po环境压力，Pa;g 重力加速度，9.8m/s2;h裂口之上液位高度，m。对于非常压下的液体泄漏速度，主要取决于内介质压力与环境压 力之

4、差和液位高低。（4）液氯容器破裂形成毒害区半径估算液化介质液氯在容器破裂时会发生蒸气扩散，从而造成大面积的 毒害区域，其毒害区按如下方法估算。设液氯质量为 W（单位：kg）,容器破裂前器内温度为t（单位：C）, 液体介质比热为 C单位：KJ/ （ kg?C）当容器破裂时，器内压力降 至大气压，处于过热状态的液化气温度迅速降至标准沸点 t （单位：C）,此时全部液体所吸收的热量为：如它的汽化热为q （单位:设这些热量全部用于器内液体的蒸发, kJ/kg），则其蒸发量：W Q W?Ct t q q(7-3)如介质的分子量为M，则在沸点下蒸发蒸气的体积 Vg （单位: m3）为：(7-4)x/ 22

5、.4Wc273 t0 22.4W ?C ? t t0 o 273 t0M ?qVg ? ?M 273 M ?q 273假设这些带有氯气的空气以半球形向地面扩散，则可求出该氯气 扩散半径为：Vg/C。14Vg/C。2.0994(7-5)式中R有毒气体的半径，m;Vg有毒介质的蒸气体积，m3 ；Co有毒介质在空气中的危险浓度， %；（5）气体连续扩散高斯分布2 2y z2 2y z(7-6)采用世界银行提供的中等密度云连续扩散模型 （高斯扩散模型）， 分析时，根据分析源周边的实际情况，考虑地面有效粗糙度并进行修 正。连续排放：Qc x,y,z expy zU式中c （ x,y,z） 连续排放时，给

6、定地点的污染物浓度, mg/m3；Q 连续排放的物料流量，mg/s ；u 平均风速，m/s;x 下风向距离，m;y 横风向距离，m；z离地面的距离，m;皈（5 x,y,z方向扩散参数，按下表选择:表2不同大气稳定度下的扩散参数大气稳定度ggA0.22x ( 1+0.0001x) -1/20.20xB0.16x ( 1+0.0001x) -1/20.123xC0.11x ( 1+0.0001x) -1/20.08x ( 1+0.0002x) -1/2D0.08x ( 1+0.0001x) -1/20.06x ( 1+0.0015x) -1/2E0.06x ( 1+0.0001x) -1/20.0

7、3x ( 1+0.0003x) -1/2F0.04x ( 1+0.0001x) -1/20.016x ( 1+0.0003x) -1/2有效粗糙度Z 0.1m的粗糙地形扩散参数按如下计算（地面有效粗 糙度长度按下表选取）：y yo f yz z0 fzfy Zo 1 0Z0fz x,Zo b clnx d e In x g0lnx表3地面有效粗糙度长度地面类型Zdm地面类型Zdm草原、平坦开阔地农作物地区村落、分散的树林 0.10.1-0.30.3-1分散的高矮建筑物（城市）密集的高矮建筑物（大城市）1-44上式中，大气稳定度为D，各参数按下表选择:表4不同大气稳定度下的系数值稳定度ABCDE

8、F0.0420.1150.150.380.30.57b01.101.51.492.532.42.913a0.03640.0450.01820.130.110.0944d00.43640.8530.870.550.860.753e00.050.012380.010460.0420.016820.0228f0.2730.1560.0890.350.270.29g00.0240.01360.00710.030.0220.023（6）气体瞬时扩散高斯分布采用世界银行提供的中等密度云瞬时扩散模型（高斯扩散模型)。2 2 2(6-7)2Q 1 x ut y2 Z2cx,y, z,t 372 exp - 2

9、 2 22 x y z 2 x y z式中c（ x,y,z,t） 瞬时排放时，给定地点（x,y,z）和时间t的污染物浓度，mg/m3;Q 瞬时排放的物料流量，mg/s ；t瞬时排放时，污染物运行的时间，s；u 平均风速，m/s ；x 下风向距离，m;y 横风向距离，m；z离地面的距离，m;笳 x,y,z方向扩散参数，地面有效粗糙度修正等因素的修正方法同上。瞬时源时，给定取样时间 t （ s）,则按下式修正：氯气管道连续泄漏模拟计算(1)事故模式设定 选取氯气从液氯钢瓶连接管到真空加氯机之间的管道、阀门和 压力表泄漏事故进行分析，按管道内径 DN8 计算 : 压力： 0.1MPa 泄漏裂口设定：

10、设裂口直径为 DN=8mm ， 则裂口面积为 5.027X10-5m2(2)泄漏模拟计算结果选取裂口面积为 A = 5.027X10Pm2,压力：O.IMpa, Cd= 0.65；p=1.468kg/m3;P=202660Pa；P0=101330Pa; H=0m 。输入相关参数代入 “有毒有害物质泄漏扩散模拟评价与风险分析 系统”软件计算 ，泄漏模型选择：气体泄漏；扩散模型选择：连续 排放；得计算结果如下：泄漏速度 =0.0376 kg/s； 气体流动情况分析：音速流动。按泄漏时间 15分钟计，得泄漏量为 33.84 kg。(3)风速为1.6m/s时，中毒危害距离取xx市年平均风速：U = 1

11、.6m/s;代入软件计算，得扩 散模拟计算结果详见下表。(4)风速为1.0m/s时，中毒危害距离取风速：U = 1.0m/s;代入软件计算，得扩散模拟计算结果详见 下表。表 5 在 1.6m/s 风速下，管道泄漏氯气扩散的中毒危害距离距离(A =0.00005026吊)氯气浓度对应的中毒危害距离（ m）50g/m3800 m g/m330 m g/m39q0 mg/m351 g/m31m g/m351 mmg/m3下风中毒危害距离618235 97411642 693横风中毒危害距离341840283 99表6在1.0m/s风速下，管道泄漏氯气扩散的中毒危害距离距离(A =0.00005026

12、m)氯气浓度对应的中毒危害距离（ m）50g/m3800mg/m330 m g/m39uJ0 mg/m351 g/m31m g/m351 mmg/m3下风中毒危害距离126396 125841473 934横风中毒危害距离351986205 133液氯钢瓶阀门或易熔塞破损时连续泄漏模拟分析1）事故模式设定选取液氯瓶在起吊、装卸、使用过程中，更换瓶头接管、易熔塞 针阀损坏等无法及时抢修，导致液氯连续泄漏事故模式进行分析， 根 据XX水厂液氯仓库的实际操作情况，选取单个液氯实瓶作为事故模 式进行分析。型号规格：800L储存温度：计算时取293K储存压力：0.6MPa液氯密度:1470kg/m3单瓶

13、液氯重量：1吨2）分析模型简介按连续扩散浓度分布模型进行分析。详见7.1所述。3）单个氯实瓶连续泄漏单个氯实瓶液氯最大存量为1000kg以准静风和平均风速分别为 1.0和1.6m/s，大气稳定度为D级进行计算分析，模拟出氯气在空间 中的浓度分布。将有关数据代入“有毒有害物质泄漏扩散模拟评价与风险分析系 统”软件进行计算，7种空间浓度计算结果如下表所示。单个氯实瓶发生单个阀连续泄漏时的计算结果见下表。表7在1.6m/s风速下，氯气扩散的中毒危害距离距离(A = 0.00005026m2)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/

14、m 35mg/m 31mg/m3下风中毒危害距离31551051942895551570横风中毒危害距离4815284179224表8 在1.0m/s风速下，氯气扩散的中毒危害距离距离(A = 0.000050260?)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/m 35mg/m 31mg/m3下风中毒危害距离40701362543807442178横风中毒危害距离610203654106311单个氯实瓶2个针阀和1个易熔塞同时发生连续泄漏时的计算 结果见下表。表9 在1.6m/s风速下，氯气扩散的中毒危害距离距离(A = 0.

15、00015078m2)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/m 35mg/m 31mg/m3下风中毒危害距离569919436655511143447横风中毒危害距离814285279159493表10 在1.0m/s风速下，氯气扩散的中毒危害距离距离(A 二0.00015078m)氯气浓度对应的中毒危害距离（m）850mg/m3300mg/m390mg/m330mg/m315mg/m35mg/m31mg/m3下风中毒危害距离7312925448574415264949横风中毒危害距离101836691062187074）

16、6个氯实瓶连续泄漏以准静风和平均风速分别为1.0和1.6m/s,大气稳定度为D级进 行计算分析，模拟出氯气在空间中的浓度分布。将有关数据代入“有毒有害物质泄漏扩散模拟评价与风险分析系 统”软件进行计算，7种空间浓度计算结果如下表所示。6个氯实瓶的6个针阀同时发生连续泄漏时的计算结果见下表。表11在1.6m/s风速下，氯气扩散的中毒危害距离距离(A = 0.000301560?)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/m 35mg/m 31mg/m3下风中毒危害距离8214628955585717805912横风中毒危害距离1

17、2214179122254845表12 在1.0m/s风速下，氯气扩散的中毒危害距离距离(A = 0.00030156m2)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/m 35mg/m 31mg/m3下风中毒危害距离106190380744116324818728横风中毒危害距离15275410616635512476个氯实瓶的12个针阀和6个易熔塞同时发生连续泄漏时的计 算结果见下表。表13在1.6m/s风速下，氯气扩散的中毒危害距离距离(A = 0.00090468*)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/m 35mg/m 31mg/m3下风中毒危害距离15027255511141780395815227横风中毒危害距离2239791592545652175表14 在1.0m/s风速下，氯气扩散的中毒危害距离距离(A = 0.00090468*)氯气浓度对应的中毒危害距离（m)850mg/m 3300mg/m 390mg/m 330mg/m 315mg/m 35mg/m 31mg/m 3下风中毒危害距离19635774415262481572023775横风中毒危害距离28511062183558173396