液氨储罐事故性泄漏扩散过程模拟术技术总结.docxVIP

液氨储罐事故性泄漏扩散过程模拟分析术 液氨是化工企业常用的原料，用途广泛，而每年因为液氨的泄漏造成的事故也十分频繁，由于其毒性很大，吸入毒性指数(Index of Potential Inhalation Toxicity，Prr)300，危险等级2，属于高度危险物质，一旦泄漏极可能造成严重的事故后果。决定液氨泄漏状况的因素多而复杂，与其理化性质、闪蒸系数、泄漏源的压力和几何形状、泄漏地的地貌情况和气象条件、储存运输的操作程序等都有密切关系。因此，综合考虑各种因素，建立液氨泄漏和扩散膜性，运用数学方法进行模拟，分析其泄漏和扩散的规律，对于救灾、重大危险源编制应急事故预案以及对新建项目进行危险性预评价都具有一定 程度的指导意义。 1 数学模型 通常情况下，液氨在常温下加压压缩，液化储存，一旦泄漏到空气中会在常压下迅速膨胀，大量气化，并扩散到大的空间范围。 1.1泄漏模型 对于灾难性破坏引起的液氨泄漏，可保守地认为容器内所有的贮存物质瞬间全部泄漏，全部泄漏时一般有爆炸发生，对其发生爆炸后的状况再运用数值模拟进行预测意义不大。因此，文中所研究的是液氨储罐连续性泄漏的数值模拟。 通过对建国50年以来我国化工系统所发生的重（特）大、典型事故性泄漏的统计分析表明[1]，阀门或法兰处的密封失效及阀门或管道断裂是造成事故 性泄漏的主要原因，因而可以确定液氨储罐下方的液氨出口接管、储罐上方的气氨出口接管以及安全阀为主要泄漏源。 1.1.1液氨泄漏模型[2] · 液氨通过其出口接管泄漏可等效为液体通过受压储罐上的孔洞泄漏。虽然氨在常温常压下为气体，但是由于泄漏发生在液相空间，流动阻力较大，故系统内压下降缓慢，不会发生因大量液氨闪蒸而造成的蒸气爆炸。另外，由于泄漏路径较短，来不及形成汽化核心而使部分液氨在池漏管道中汽化而形 成闪蒸两相流。因此，其泄漏速率可采用式(1)计 算[3]: Qm=PACo[2 (P0/p+ghr)]1／2 (1) 式中：Qm为质量泄漏速率，kg/s；Co为泄漏系数；A为裂口面积，㎡；PO为储罐内压，Pa;hr是泄漏处与液面之间的距离，m。根据式(1)，随着储罐渐渐变空，液体高度减少，流速和质量流速也随之??少。 泄漏的液氨会发生闪蒸，其闪蒸率的大小可由式（2）计算： =（H1-H2）/Hv 式中：Mv，为闪蒸蒸气的质量，kg; Mo为泄漏液体的总质量，kg;H1为液体储存温度时的焓，J/kg; H2为常压下液体沸点Tb时的焓，J/kg:Hv为泄漏液体的蒸发热，J/kg;泄漏所造成的气氨在空气中的浓度分布只是空间的函数。 通常情况下，当闪蒸率 Mv/Mo≥0.2时，则泄漏的液氨全部闪蒸，当Mv/Moo．2时，则泄漏的液氨会在地面形成液池。 1.1.2 气氨泄漏模型 对于气氨通过其出口接管的泄漏，情况较复杂。由于出口接管处于气相空间，其泄漏形式主要与泄漏面积的大小有关。在泄漏面积较大的情况下，高 压蒸气通过裂缝或孔洞喷出，储罐内压急剧下降，直到环境压力（常温）。由于内压急剧下降，气液平衡遭到破坏，储罐内液氨处于过热状态，过热状态的液 氨为了再次恢复平衡，内部会均匀地产生沸腾核，同时产生大量气泡，液氨体积急剧膨胀，最终导致爆炸；当泄漏面积不大时，即使有蒸气喷出，但由于储罐内压下降不急剧，液氨不会达到过热状态，因此不会发生蒸气爆炸，其泄漏速率可采用下式计算： Qm= CoAPo｛2rM/RgTo(r-1)[(P/Po)2/r-(P/Po)r+1/r]｝1/2 式中：Q为质量泄漏速率，kg/s; Co为泄漏系数；A为裂口面积，m2; PO为储罐内压，Pa;M为气体或蒸气的摩尔质量，kg/mol; Rg为理想气体常数； To为泄漏源温度，K；P为泄漏处压力，Pa;r为绝热指数。 1.2扩散模型 1.2.1 液氨扩散模型 液氨泄漏后发生闪蒸的氨由于夹带有氨的液滴，因此其密度比周围空气的密度大，在一段时间内其扩散行为受到重力沉降的作用，表现为重气的扩散。重气扩散模式可采用式(4)计算： b=bo｛1+1.5[gho(?p-?pa)/pa]1/21/vbo} INCLUDEPICTURE \d C:\\Users\\Administrator\\AppData\\Roaming\\Tencent\\Users\\1615113736\\QQ\\WinTemp\\RichOle\\@{STK5N{QK45GY)H[ADA_D2.png \* MERGEFORMATINET  式中：b为重气云羽的横风向半宽，m；bo为泄漏源点重气云羽的横风向半宽，m；ho为泄漏源点重气云羽高度，m； pa为空气密度，kg/m3；v为重气云羽的轴向蔓延