液氨储罐事故性泄漏扩散过程模拟分析术

1、液氨储罐事故性泄漏扩散过程模拟分析术液氨是化工企业常用的原料，用途广泛，而每年因为液氨的泄漏造成的事故也十分频繁，由于其毒性很大，吸入毒性指数(IndexofPotentialInhalationToxicity,Prr)<300，危险等级2,属于高度危险物质，一旦泄漏极可能造成严重的事故后果。决定液氨泄漏状况的因素多而复杂，与其理化性质、闪蒸系数、泄漏源的压力和几何形状、泄漏地的地貌情况和气象条件、储存运输的操作程序等都有密切关系。因此，综合考虑各种因素，建立液氨泄漏和扩散膜性，运用数学方法进行模拟，分析其泄漏和扩散的规律，对于救灾、重大危险源编制应急事故预案以及对新建项目进行危险性预

2、评价都具有一定程度的指导意义。1数学模型通常情况下，液氨在常温下加压压缩，液化储存，一旦泄漏到空气中会在常压下迅速膨胀，大量气化，并扩散到大的空间范围。1.1 泄漏模型对于灾难性破坏引起的液氨泄漏，可保守地认为容器内所有的贮存物质瞬间全部泄漏，全部泄漏时一般有爆炸发生，对其发生爆炸后的状况再运用数值模拟进行预测意义不大。因此，文中所研究的是液氨储罐连续性泄漏的数值模拟。通过对建国50年以来我国化工系统所发生的重(特)大、典型事故性泄漏的统计分析表明1，阀门或法兰处的密封失效及阀门或管道断裂是造成事故性泄漏的主要原因，因而可以确定液氨储罐下方的液氨出口接管、储罐上方的气氨出口接管以及安全阀为主要

3、泄漏源。1.1.1 液氨泄漏模型2液氨通过其出口接管泄漏可等效为液体通过受压储罐上的孔洞泄漏。虽然氨在常温常压下为气体，但是由于泄漏发生在液相空间，流动阻力较大，故系统内压下降缓慢，不会发生因大量液氨闪蒸而造成的蒸气爆炸。另外，由于泄漏路径较短，来不及形成汽化核心而使部分液氨在池漏管道中汽化而形成闪蒸两相流。因此，其泄漏速率可采用式(1)计算3:Qm=PACo2(P0/p+ghr)12(1)式中：Qm为质量泄漏速率，kg/s；Co为泄漏系数；A为裂口面积，近；PO为储罐内压，Pa;hr是泄漏处与液面之间的距离，m。根据式(1),随着储罐渐渐变空，液体高度减少，流速和质量流速也随之减少。泄漏的液

4、氨会发生闪蒸，其闪蒸率的大小可由式(2)计算：=(H1-H2)/Hv式中:Mv，为闪蒸蒸气的质量，kg;Mo为泄漏液体的总质量，kg;H1为液体储存温度时的焓，J/kg;H2为常压下液体沸点Tb时的焓，J/kg:Hv为泄漏液体的蒸发热，J/kg;泄漏所造成的气氨在空气中的浓度分布只是空间的函数。通常情况下，当闪蒸率Mv/Mo三0.2时，则泄漏的液氨全部闪蒸，当Mv/Mo<o2时，则泄漏的液氨会在地面形成液池。1.1.2 气氨泄漏模型对于气氨通过其出口接管的泄漏，情况较复杂。由于出口接管处于气相空间，其泄漏形式主要与泄漏面积的大小有关。在泄漏面积较大的情况下，高压蒸气通过裂缝或孔洞喷出，储

5、罐内压急剧下降，直到环境压力(常温)。由于内压急剧下降，气液平衡遭到破坏，储罐内液氨处于过热状态，过热状态的液氨为了再次恢复平衡，内部会均匀地产生沸腾核，同时产生大量气泡，液氨体积急剧膨胀，最终导致爆炸；当泄漏面积不大时，即使有蒸气喷出，但由于储罐内压下降不急剧，液氨不会达到过热状态，因此不会发生蒸气爆炸，其泄漏速率可采用下式计算：Qm=CoAPo2rM/RgTo(rT)(P/Po)2/r-(P/Po)r+1/r1/2式中：Q为质量泄漏速率，kg/s;Co为泄漏系数；A为裂口面积，m2;PO为储罐内压，Pa;M为气体或蒸气的摩尔质量，kg/mol;Rg为理想气体常数；To为泄漏源温度，K；P为

6、泄漏处压力，Pa;r为绝热指数。1.2扩散模型1.2.1液氨扩散模型液氨泄漏后发生闪蒸的氨由于夹带有氨的液滴，因此其密度比周围空气的密度大，在一段时间内其扩散行为受到重力沉降的作用，表现为重气的扩散。重气扩散模式可采用式(4)计算：b=bo1+15gho(p-pa)/pai/1/vbo(itv2bh式中：b为重气云羽的横风向半宽，m；bo为泄漏源点重气云羽的横风向半宽,m；ho为泄漏源点重气云羽高度，m；pa为空气密度，kg/m3；v为重气云羽的轴向蔓延速度,m/s；h为重气云羽高度,m；we为空气卷吸系数。泄漏后的液氨在空气中稀释后表现为非重气云的扩散，此时其扩散模式为气氨扩散模式。1.2

7、2气氨扩散模型气态氨在大气中的扩散，可以用高斯模型来计算其危害范围。高斯模型包括高斯烟羽模型和高斯烟团模型。其中，高斯烟羽模型适用于连续源的气体扩散，而高斯烟团模型适用于瞬时源的气体扩散4。在泄漏发生的初始一段时间内，其浓度分布是不稳定的，是空间和时间的函数，此时采用高斯烟团迭加模型（高斯烟羽模型与高斯烟团模型的结合）进行描述6，其数学表达式如下：式中：C'为泄漏介质在大气中的摩尔百分比浓度Qm。为质量泄漏速率，kg/s；u为环境平均风速，m/s；t为泄漏时间，s；Hr为有效源高，m；6x、6y、6z为预测点坐标，m；Ox，q，o分别为下风向、横风向和竖直风向的扩散系数，m。1-

8、9;1aL孔.-引(c尤,V盅)二112玛如exp巴任_比广(z+)2-exp2+esp珂-(6)式中：Qm为源的泄放速率，kg/s°6x、6y、6z是高斯模型的重要参数，6x、6y、oz呸是由排放源到计算点的下风向距离和大气稳定度的函数，也与烟羽的排放高度及地面粗糙度有关。1.3 毒负荷氨气为毒性物质。人类对毒物的反应不仅与所接触的毒物的浓度有关，而且还与暴露在该浓度下的时间有关，因此，为了同时反映毒物浓度和暴露时间对中毒反应的影响，人们提出了毒负荷的概念，其定义如下：式中：TL为毒物的毒负荷，它决定中毒程度；tl、t2分别为开始暴露时刻和结束暴露时刻，min；C为t时刻暴露环境中

9、毒物的浓度；n为幂指数，对于氨气，n取2。2模拟分析针对某市化学园区某化工公司液氨储罐工程建设项目应用以上数学模型进行计算机模拟，假设液氨储罐底部由于法兰接头处垫片腐蚀破裂而发生连续性泄漏，泄漏当量直径为2cm，模拟计算的基础数据如表1、袁2、表3所示，得到计算结果如图1图10所示。表1液氨理化特性液氨理化特性相魁分子虽沸点相对密段村对密厦（空气二1）饱和裁汽压空气中允许蔽度(水=n（kPti他了）(%)17.03-315(1S2a6S0&6215.7-27.斗30表裂存储条件装数；充MW鑼型储类306卧罐15H)SO表3气象条件平均气温平均凤速辐太阳辐(E)m/皆)风向射强度IS.3

10、27东南风中等13液氨连续稳定泄漏影响区域图图I图2液氨连续稳定泄漏下凤向100m处地面浓度随时间的变化曲线圉3液氨连续穗定泄漏下凤向WOin处地匱浓度毒员荷随时间的变化曲线圉4液氨连绩稳定泄漏下风向2001齐V沁摊丙济障屢时闸朋变代曲縁5液氨连续稳定泄漏下风向200m处地面浓度毒奂荷随时间的变化曲线tocominutes图$液氨连续稳定泄漏下凤向处地面浓度随时间的变化曲线S«H<I7图7液氨连续稳定泄漏下风向500ni处地面浓度莓员荷随时间的变化曲线."-ymlnuiet團8洒氨连续泄漏下冈向1<MM»m处地面滾度随时间的变化曲线300山12OQIC

11、O'Q“-r0W20304050nuAutes图10液氨连续稳定泄漏速率随时间的变化3模拟分析结果图1为液氨连续泄漏影响区域示意图，图中阴影部分表示氨气浓度大于空气中允许最高浓度，即30mg/m3。由图1可见，液氨泄漏后其最远影响范围已经达到下风向3200m处，横风向的最大影响距离也达到250m。由于模拟过程中假设泄漏工厂周围地势平坦，因此由模拟结果得出的氨气影响范围偏大：一但这对于指导今后的紧急泄漏事故应急处理方案是有益的。图1中阴影部分外围的虚线表示氨气云团有可能影响到的区域，主要是风向改变的影响。液氨连续泄漏时，下风向100m、200m、500m及1000m处室内、外氨气浓度随时间的变化分别见图2、图4、图6和图8。图中LOC线代表空气中的最高允许浓度（3mg/m3），虚线代表室内浓度随时间的变化，实线代表室外浓度随时间的变化。从图中可以看出，室外浓度在某一段时间内很高，远远超过LOC浓度，而室内浓度在泄漏一段时间后才高于LOC浓度，且这一时间随下风向距离的增大而增大。下风向100m、200m、500m及1000m处室内、外氨气毒负荷随