【TIEMS2014论文】液氨储罐泄漏扩散事故模拟分析.docVIP

液氨储罐泄漏扩散事故模拟分析 徐亚博 宋冰雪 葛悦 谢昱姝 北京市劳动保护科学研究所，北京 100054 摘要：液氨泄漏初期由于蒸发以及气云带走液滴，在地面形成雾体，多形成重气云团，以气雾状态时的扩散应属于重气扩散；在扩散过程中随着空气的不断卷吸，密度不断变小，重气云的密度会越来越接近空气，从而转变为非重气云。本文应用重气扩散及高斯公式，建立液氨泄漏扩散数学模型，将重气云扩散模型和非重气云的高斯模型有机结合起来，估测液氨泄漏扩散连续变化过程中氨气浓度的分布规律及危害范围。相应的数学模型可为救灾、重大危险源编制应急事故预案以及对新建项目进行危险性预评价提供科学评价方法，具有一定程度的指导意义。 关键词：液氨；泄漏；扩散；应急 0 引言 液氨是化工企业的常用原料，常用于冷库、化肥、制药、塑料、合成纤维、石油精炼等行业。氨在常温常压下，为具有刺激性气味的无色气体，加压易被液化为无色液体（液氨）。氨是一种有毒物质，当空气中氨浓度达到一定时，会使人产生呼吸道严重中毒症状，甚至可使人发生“闪电式”死亡。在生产、储存和使用过程中 2013年8月31日，上海市宝山区丰翔路1258号翁牌冷藏实业有限公司，造成15人死亡、26人受伤 2013年6月3日，吉林省德惠市吉林宝源丰禽业有限公司121人死亡、76人受伤，直接经济损失1.82亿元。 2013年4月21日，四川省眉山市仁寿县凤陵乡金凤食品厂生猪屠宰场冻库由于液氨泄漏导致急性中毒事故，造成4人死亡，22人中毒。 本文针对液氨储罐泄漏引发的急性中毒事故应急救援进行模拟研究，利用数学模型分析计算液氨扩散浓度的分布规律，以及根据毒性对人体的影响，得出不同扩散浓度对人体的伤害程度，从而划定危险区域的影响范围，作为制定应急救援预案的依据，将事故发生时对人员和财产的损失降到最低。 1 数学模型 通过对典型事故性泄漏的统计分析表明，阀门或法兰处的密封失效及阀门或管道断裂是造成事故性泄漏的主要原因，因而可以确定液氨储罐下方的液氨出口接管、储罐上方的气氨出口接管以及安全阀为主要泄漏源[1,2]。 一旦发生泄漏，液氨到空气中会在常压下迅速膨胀，大量气化，并扩散到大的空间范围。液氨泄漏后扩散定量计算比较复杂，在以往的计算多是采用湍流的高斯扩散模型计算，这是不完全正确的。事实上，液氨泄漏扩散计算应分为二个阶段：初始阶段，液氨泄漏初期由于蒸发以及气云带走液滴，在地面形成雾体，多形成重气云团，以气雾状态时的扩散应属于重气扩散；在扩散过程中随着空气的不断卷吸，密度不断变小，重气云的密度会越来越接近空气，从而转变为非重气云，此时应按高斯扩散模型计算。扩散计算步骤：判定为重气——重气扩散计算——扩散转变点——高斯扩散计算 图1 液氨泄漏计算过程 1.1泄漏计算模型 由于泄漏发生液相空间，流动阻力较大，故系统内压下降缓慢，不会发生因大量液氨闪蒸而造成的蒸汽爆炸。另外，由于泄漏路径较短，来不及形成汽化核心而使部分液氨在泄漏管道中汽化而形成闪蒸两相流[3]。因此，其泄漏速率可用流体力学的伯努利方程计算： （1） 式中，Q0液体泄漏速度，kg/s；Cd液体泄漏系数，；A裂口面积，m2；泄漏液体密度，kg/ m3；P容器内介质压力，Pa；P0环境压力，Pag为重力加速度，9.8 m/ s2；h裂口之上液位高度，m。 泄漏时直接蒸发的液体所占百分比F可按下式计算： （2） 式中，Cp液体的定压比热，J/kg·K；T泄漏前液体的温度，K；T0液体在常压下的沸点，K；H液体的汽化热，J/kg。 根据经验，当F＞0.2时，一般不会形成液池；当F＜0.2时，F与带走液体之比，有线性关系，即当F=0时，没有液体被带走（蒸发）；当F=0.1时，有50%的液体被带走。 1.2.1判断是否重气云 如果成立，扩散模型选用非重气扩散模型，否则，选用重气扩散模型。 （3） 式中：瞬时泄漏形成的云团的初始体积，；为液体常压下沸点，液氨的为239.5为分子量，氨为17为闪蒸蒸汽的质量， 为气云初始密度，； 为闪蒸蒸汽的质量， V为环境风速，； 环境空气密度，，约为1.239。 重气云团内部氨气浓度为： 式中： ；初始浓度等于初始密度；x为下风向距离，m为瞬时泄漏形成的云团的初始体积 1.2.3扩散转变点 运用准则判断重气坍塌过程终止位置 定义。准则认为，如果小于或等于某个临界值，重气坍塌引起的扩散将让位于环境湍流引起的扩散。 将式，代入求出，得到 由于不考虑云团和环境之间的热交换，云团能量守恒，，代入上式得到转变点对应的下风向距离为： （5） 式中 ；E0为云团初始能力，计算公式，其中，，为空气密度，1.239 1.2.4高斯扩散模型 高斯模型用来描述