氯气泄漏扩散计算模拟

1 氯氯气泄漏气泄漏扩扩散散计计算模算模拟拟 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 区 引起剧咳 30 9010 30轻伤区 C 区 刺激咽喉 引起剧咳 2 155吸入反应区 D 区 刺激咽喉 62有明显气味 刺激眼 鼻 10 03嗅觉浓度 3 液氯泄漏量计算 流体泄漏速度可用流体力学的柏努利方程计算 其泄漏速度为 7 1 gh pp ACdQ2 2 0 式中 Q0 液体泄漏速度 kg s Cd 液体泄漏系数 A 裂口面积 m2 泄漏液体密度 kg m3 P 容器内介质压力 Pa P0 环境压力 Pa g 重力加速度 9 8m s2 h 裂口之上液位高度 m 对于非常压下的液体泄漏速度 主要取决于内介质压力与环境压 力之差和液位高低 4 液氯容器破裂形成毒害区半径估算 液化介质液氯在容器破裂时会发生蒸气扩散 从而造成大面积的 毒害区域 其毒害区按如下方法估算 设液氯质量为 W 单位 kg 容器破裂前器内温度为 t 单位 液体介质比热为 C 单位 KJ kg 当容器破裂时 器内压力降至大 气压 处于过热状态的液化气温度迅速降至标准沸点 t0 单位 此 时全部液体所吸收的热量为 7 2 0 ttCWQ 3 设这些热量全部用于器内液体的蒸发 如它的汽化热为 q 单位 kJ kg 则其蒸发量 7 3 q ttCW q Q W 0 如介质的分子量为 M 则在沸点下蒸发蒸气的体积 Vg 单位 m3 为 7 4 273 273 4 22 273 273 4 22 000 t qM ttCWt M W Vg 假设这些带有氯气的空气以半球形向地面扩散 则可求出该氯气 扩散半径为 7 5 3 0 3 0 0994 2 3 4 2 1 CVgCVg R 式中 R 有毒气体的半径 m Vg 有毒介质的蒸气体积 m3 CO 有毒介质在空气中的危险浓度 5 气体连续扩散高斯分布 采用世界银行提供的中等密度云连续扩散模型 高斯扩散模型 分析时 根据分析源周边的实际情况 考虑地面有效粗糙度并进行修 正 连续排放 7 6 2 2 2 2 2 1 exp zyzy Zy u Q zyxc 式中 c x y z 连续排放时 给定地点的污染物浓度 mg m3 Q 连续排放的物料流量 mg s u 平均风速 m s x 下风向距离 m 4 y 横风向距离 m z 离地面的距离 m x y z x y z 方向扩散参数 按下表选择 表 2 不同大气稳定度下的扩散参数 大气稳定度 y z A B C D E F 0 22x 1 0 0001x 1 2 0 16x 1 0 0001x 1 2 0 11x 1 0 0001x 1 2 0 08x 1 0 0001x 1 2 0 06x 1 0 0001x 1 2 0 04x 1 0 0001x 1 2 0 20 x 0 123x 0 08x 1 0 0002x 1 2 0 06x 1 0 0015x 1 2 0 03x 1 0 0003x 1 2 0 016x 1 0 0003x 1 2 有效粗糙度 Z0 0 1m 的粗糙地形扩散参数按如下计算 地面有效粗 糙度长度按下表选取 yyoy f zzz f 0 000 1ZZfy xgf z ZxedxcbZxf ln 0 1 00000 00 lnln 表 3 地面有效粗糙度长度 地面类型Z0 m地面类型Z0 m 草原 平坦开阔地农作物地区 村落 分散的树林 0 1 0 1 0 3 0 3 1 分散的高矮建筑物 城市 密集的高矮建筑物 大城市 1 4 4 上式中 大气稳定度为 D 各参数按下表选择 表 4 不同大气稳定度下的系数值 稳定度ABCDEF a00 0420 1150 150 380 30 57 5 b0 c0 d0 e0 f0 g0 1 10 0 0364 0 4364 0 05 0 273 0 024 1 5 0 045 0 853 0 01238 0 156 0 0136 1 49 0 0182 0 87 0 01046 0 089 0 0071 2 53 0 13 0 55 0 042 0 35 0 03 2 4 0 11 0 86 0 01682 0 27 0 022 2 913 0 0944 0 753 0 0228 0 29 0 023 6 气体瞬时扩散高斯分布 采用世界银行提供的中等密度云瞬时扩散模型 高斯扩散模型 6 7 2 2 2 2 2 2 2 3 2 1 exp 2 2 zyxzyx ZyutxQ tzyxc 式中 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 x y z 方向扩散参数 地面有效粗糙度修正等因素 的修正方法同上 瞬时源时 给定取样时间 t s 则按下式修正 2 0 min10 600 t yy 6 氯氯气管道气管道连续连续泄漏模泄漏模拟计拟计算算 1 事故模式设定 选取氯气从液氯钢瓶连接管到真空加氯机之间的管道 阀门和 压力表泄漏事故进行分析 按管道内径 DN8 计算 压力 0 1MPa 泄漏裂口设定 设裂口直径为 DN 8mm 则裂口面积为 5 027X10 5m2 2 泄漏模拟计算结果 选取裂口面积为 A 5 027X10 5m2 压力 0 1Mpa Cd 0 65 1 468kg m3 P Pa P0 Pa H 0m 输入相关参数代入 有毒有害物质泄漏扩散模拟评价与风险分析 系统 软件计算 泄漏模型选择 气体泄漏 扩散模型选择 连续排放 得计算结果如下 泄漏速度 0 0376 kg s 气体流动情况分析 音速流动 按泄漏时间15分钟计 得泄漏量为33 84 kg 3 风速为 1 6m s 时 中毒危害距离 取 xx 市年平均风速 U 1 6m s 代入软件计算 得扩散 模拟计算结果详见下表 4 风速为 1 0m s 时 中毒危害距离 取风速 U 1 0m s 代入软件计算 得扩散模拟计算结果详见 下表 表 5 在 1 6m s 风速下 管道泄漏氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 8 50 3 00 9 0 3 0 1 5 5 m 1 m 7 m g m3 m g m3 m g m3 m g m3 m g m3 g m3g m3 下风中毒 危害距离 1 6 2 8 5 3 9 7 1 41 2 64 6 93 横风中毒 危害距离 348 1 4 2 0 3 8 9 9 表 6 在 1 0m s 风速下 管道泄漏氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 8 50 m g m3 3 00 m g m3 9 0 m g m3 3 0 m g m3 1 5 m g m3 5 m g m3 1 m g m3 下风中毒 危害距离 2 1 3 6 6 9 1 25 1 84 3 47 9 34 横风中毒 危害距离 359 1 8 2 6 5 0 1 33 8 9 液液氯钢氯钢瓶瓶阀门阀门或易熔塞破或易熔塞破损时连续损时连续泄漏模泄漏模拟拟分析分析 1 事故模式设定 选取液氯瓶在起吊 装卸 使用过程中 更换瓶头接管 易熔塞针 阀损坏等无法及时抢修 导致液氯连续泄漏事故模式进行分析 根据 xx 水厂液氯仓库的实际操作情况 选取单个液氯实瓶作为事故模式 进行分析 型号规格 800L 储存温度 计算时取 293K 储存压力 0 6MPa 液氯密度 1470kg m3 单瓶液氯重量 1 吨 2 分析模型简介 按连续扩散浓度分布模型进行分析 详见 7 1 所述 3 单个氯实瓶连续泄漏 单个氯实瓶液氯最大存量为 1000kg 以准静风和平均风速分别 为 1 0 和 1 6m s 大气稳定度为 D 级进行计算分析 模拟出氯气在空 间中的浓度分布 将有关数据代入 有毒有害物质泄漏扩散模拟评价与风险分析系 统 软件进行计算 7 种空间浓度计算结果如下表所示 单个氯实瓶发生单个阀连续泄漏时的计算结果见下表 表 7 在 1 6m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m3 300 mg m3 90 mg m3 30 mg m3 15 mg m3 5 mg m3 1 mg m3 下风中毒危害距离31551051942895551570 横风中毒危害距离4815284179224 10 表 8 在 1 0m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m3 300 mg m3 90 mg m3 30 mg m3 15 mg m3 5 mg m3 1 mg m3 下风中毒危害距离40701362543807442178 横风中毒危害距离610203654106311 单个氯实瓶 2 个针阀和 1 个易熔塞同时发生连续泄漏时的计 算结果见下表 表 9 在 1 6m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m3 300 mg m3 90 mg m3 30 mg m3 15 mg m3 5 mg m3 1 mg m3 下风中毒危害距离569919436655511143447 横风中毒危害距离814285279159493 表 10 在 1 0m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m 3 300 mg m 3 90 mg m 3 30 mg m3 15 mg m 3 5 mg m3 1 mg m3 下风中毒危害距离7312925448574415264949 横风中毒危害距离10183669106218707 4 6 个氯实瓶连续泄漏 以准静风和平均风速分别为 1 0 和 1 6m s 大气稳定度为 D 级进 行计算分析 模拟出氯气在空间中的浓度分布 将有关数据代入 有毒有害物质泄漏扩散模拟评价与风险分析系 统 软件进行计算 7 种空间浓度计算结果如下表所示 6 个氯实瓶的 6 个针阀同时发生连续泄漏时的计算结果见下 11 表 表 11 在 1 6m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m3 300 mg m3 90 mg m3 30 mg m3 15 mg m3 5 mg m3 1 mg m3 下风中毒危害距离8214628955585717805912 横风中毒危害距离12214179122254845 表 12 在 1 0m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m3 300 mg m3 90 mg m3 30 mg m3 15 mg m3 5 mg m3 1 mg m3 下风中毒危害距离106190380744116324818728 横风中毒危害距离1527541061663551247 6 个氯实瓶的 12 个针阀和 6 个易熔塞同时发生连续泄漏时的 计算结果见下表 表 13 在 1 6m s 风速下 氯气扩散的中毒危害距离 氯气浓度对应的中毒危害距离 m 距离 A 0 m2 850 mg m3 300 mg m3 90 mg