氨气气体泄漏后果分析

-.z摘要本次课程设计是对液氨储罐事故后果进展定量分析。根据液氨的理化性质和主要危险特性,采用事件树分析法对液氨储罐的事故类型进展了危险源辨识和分析。并通过对泄漏、火灾、爆炸等典型事故影响模型进展了定量计算，分析了液氨泄漏的危害*围及其严重程度,提出了科学、合理、可行的平安对策措施和建议，从而为预防液氨泄漏事故发生和进展平安管理提供了依据。关键词：液氨储罐；危险源；事故影响模型1、储罐区的情况简介1.1储罐区的根本情况本次课程设计以***发电厂为研究对象。该发电厂采用选择性催化复原法进展烟气脱硝,所用复原剂为液氨,共有2个液氨储罐,每个液氨储罐最大储存量为90m3。液氨储存温度为30℃,储存压力为1.1MPa,密度为750kg/m3,液氨体积占储罐容积的最大值为70%(其充装系数为0.70)。则每个贮槽内液氨的总质量为W=90m3×750kg/m3×0.7=47.25t。重大危险源，是指长期地或者临时地生产、搬运、使用或者储存危险物品、且危险物品的数量等于或者超过临界量的单元〔包括产所和设施〕。?危险化学品重大危险源辨识?【2】(GB18218—2021〕规定氨的临界量为10t，该企业布置有两个液氨储罐，每个储罐装存的液氨为47.25t，因此构成了该液氨储罐区构成了重大危险源。1.2液氨的主要危险特性液氨又称为无水氨，是一种无色液体。在温度变化时，液氨体积变化的系数很大。溶于水、乙醇和乙醚，与空气混合能够形成爆炸混合物，火灾危险类别为乙类2项。液氨作为一种重要的化工原料应用广泛，普遍存在于化工生产过程中。为了运输及储存的便利，通常将气态的氨气通过加压或冷却得到液态氨。液氨在工业上应用广泛，由于具有腐蚀性，且易挥发，所以其化学事故发生率相当高，是该储罐区的主要危险物料。液氨物料的危险特性主要表达在燃烧和爆炸、活性反响和安康危害三方面【3】，具体危险特性及理化性质如表1.1所示：表1.1液氨的危险特性及理化性质表危险特性燃烧和爆炸危险性极易燃，能与空气形成爆炸性混合物，遇明火、高热引起燃烧爆炸。活性反响与氟、氯等接触会发生剧烈的化学反响。安康危害对眼、呼吸道粘膜有强烈刺激和腐蚀作用。急性氨中毒引起眼和呼吸道刺激病症，支气管炎或支气管周围炎，肺炎，重度中毒者可发生中毒性肺水肿。高浓度氨可引起反射性呼吸和心搏停顿。可致眼和皮肤灼伤。PC-TWA(时间加权平均容许浓度)(mg/m3):20;PC-STEL(短时间接触容许浓度)(mg/m3):30。理化性质分子量为17.03，熔点-77.7℃，沸点-33.5℃，液氨密度为0.7253g/cm3,临界压力11.40MPa，临界温度132.5℃，饱和蒸气压1013kPa(26℃)，爆炸极限15%～30.2%〔体积比〕，自燃温度630℃，燃点为651℃,最大爆炸压力0.580Mpa2、储罐的事故类型分析2.1液氨泄漏事故模式及统计分析通常情况下,液氨在常温下加压压缩,液化储存,一旦泄漏到空气中会在常压下迅速膨胀,大量气化,并扩散到大的空间*围。液氨事故泄漏后通常有几种事故模式：液氨泄漏后在泄漏出口处立即点火形成喷射火；泄漏处于开放空间且经过一定时间点火形成闪火；泄漏处于局限空间条件且经过一定时间点火形成蒸气云爆炸；假设泄漏过程中没有点火源存在则形成单纯的大气中扩散；储存液氨的储罐有可能发生BLEVE爆炸。根据事故案例,泄漏事故绝大多数可视为连续点源泄漏模式。因此,本文着重分析和模拟氨气连续点源泄漏后果。2.2气氨泄漏事故树模型对于液氨储罐液体泄漏应考虑发生池火+喷射火/BLEVE、蒸汽爆炸、闪火+晚期池火等火灾爆炸场景。具体场景与泄漏类型、点火类型、空间局限情况等有关，可采用事件树方法确定各种可燃物质释放后，各种事件发生的类型及概率【5】。可燃液体泄露后的事件树见图1.1所示：IE完全闪蒸立即点火延迟点火空间局限后果是喷射火是是爆炸是否闪火液氨泄漏否否大气扩散是池火〔+喷射火/BLEVE〕是否爆炸是否否闪火〔+晚期池火〕否大气扩散图1.1可燃液体泄露事件树3氨气泄漏后果定量计算3.1气氨泄漏模型对于灾难性破坏引起的液氨泄漏,可保守地认为容器内所有的贮存物质瞬间全部泄漏,全部泄漏时一般有爆炸发生,对其发生爆炸后的状况再运用数值模拟进展预测意义不大。因此,文中所研究的是液氨储罐连续性泄漏的数值模拟。通过对建国50年以来我国化工系统所发生的重(特)大、典型事故性泄漏的统计分析说明,阀门或法兰处的密封失效及阀门或管道断裂是造成事故性泄漏的主要原因,因而可以确定液氨储罐下方的液氨出口接收、储罐上方的气氨出口接收以及平安阀为主要泄漏源。气氨泄漏模型对于气氨通过其出口接收的泄漏,情况较复杂。由于出口接收处于气相空间,其泄漏形式主要与泄漏面积的大小有关。在泄漏面积较大的情况下,高压蒸气通过裂缝或孔洞喷出,储罐内压急剧下降,直到环境压力(常温)。由于内压急剧下降,气液平衡遭到破坏,储罐内液氨处于过热状态,过热状态的液氨为了再次恢复平衡,内部会均匀地产生沸腾核,同时产生大量气泡,液氨体积急剧膨胀,最终导致爆炸;当泄漏面积不大时,即使有蒸气喷出,但由于储罐内压下降不急剧,液氨不会到达过热状态,因此不会发生蒸气爆炸3.1.1液体泄漏计算由于储罐发生破裂泄漏的概率很小，因此泄漏场景考虑接收泄漏。根据?钢制立式圆筒固定顶储罐系列?【8】HG21502.1-1992，确定接收管径为50mm。液体经管道上的孔流出的质量流率为：〔3.1〕式中:Qm——质量流率，单位为kg/s；A——泄漏孔面积，单位为m2；C0——液体泄漏系数；P——管道内液体压力，单位为Pa；ρ——泄漏液体密度，单位为kg/m3；P0——环境压力，单位为Pa。液氨的密度为0.75g/cm³，即ρ=750kg/m³。储罐内液体压力为1.1×106Pa〔表压力〕，环境压力P0=1.013×105Pa，即储罐内液体压力〔绝对压力〕P=1.1×106+1.013×105=1.2×10^6Pa，当液体泄漏系数不知道或不能确定时，取1.0，使所计算的流量最大。1〕泄漏场景为小孔泄漏时，泄漏孔径为5mm，即半径为R=2.5mm。A=πR²=1.96×10-5m2=〔1.96×10-6〕×1×[2×750×〔1.2×10^6-1.013×105)]1/2=1.106kg/s=0.7957kg/s2〕泄漏场景为中孔泄漏时，泄漏孔径为10mm，即半径为R=12.5mm。A=πR²=4.9×10-4m2=4.9×10-4×1×[2×750×〔1.2×10^6-1.013×105)]1/2=19.89kg/s3〕泄漏场景为大孔破裂时，泄漏孔径为50mm，即半径为R=25mm。A=πR²=1.96×10-3m2=1.96×10-3×1×[2×750×〔1.2×10^6-1.013×105)]1/2=79.57kg/s3.1.2闪蒸计算〔1〕液体中闪蒸局部：〔3.2〕过热液体闪蒸蒸发速率可按下式计算：〔3.3〕式中：Fv——泄漏液体的闪蒸比例；TT——储存温度，单位为K；Tb——泄漏液体的沸点，单位为K；Hv——泄漏液体的蒸发热，单位为kJ/kg；Cp——泄漏液体的定压热容，单位为kJ/(kg·K)；Q1——过热液体闪蒸蒸发速率，单位为kg/s；QL——物质泄漏速率，单位为kg/s。储存温度TT=30+273.15=303.15K，沸点Tb=239.65K，泄漏液体的蒸发热Hv=1366.70kJ/kg，定压热容Cp=4.609kJ/(kg·K)=0.21414由于Fv>0.2时，其泄露的液氨完全闪蒸：/3.2气氨扩散模型根据泄漏气体的密度和泄漏源类型，气体的扩散模式可分为烟团扩散和烟羽扩散两种模式。泄漏量较大且密度比空气的密度大得多的气体扩散呈现烟团式扩散，假设泄漏气体密度与空气接近或经很短时间的空气稀释后密度与空气接近时，其泄漏气体的扩散属于烟羽式扩散，大局部较小流量的连续源易形成烟羽扩散。因此，氨气的扩散属于烟羽扩散。通过查询**地区的气象资料，**一年内的平均风速为2.84m/s，根据?化工企业定量风险评价导则?AQ/T3046-2021表E.5可知，大气稳定度为C级，根据表E.7可知：σy=0.22*〔1+0.0004*〕^(-1/2)σz=0.20*根据Pasquill-Gifford模型扩散方程，位于地面高处的连续稳态源的烟羽〔3.6〕式中：——连续排放时，形成稳定的流场后，给定地点的污染物的浓度，单位为g/m3；——连续排放的物料质量流量，单位为g/s；u——风速，单位为m/s；，——侧风向和垂直风向的扩散系数，单位为m；*——下风向距离，单位为m；y——侧风向距离，单位为m；z——垂直风向距离，单位为m。由液体泄漏中的计算可知，当泄漏为小孔泄漏时，连续排放的物料质量流量Q=795.7g/s。因为液氨储罐水平布置在地面上，所以可以在*轴方向取10m、50m、100m、200m、400m，y轴取0，z轴取0，具体地点污染物浓度见表3.1所示：表3.1不同地点液氨扩散的浓度*〔m〕1050100200400500Y〔m〕000000Z〔m〕000000〔m〕2.195610.0718.9835.7667.3882.62〔m〕210.8921.5742.33981.705100.416C〔g/m3〕20.3090.81320.21780.05890.01620.010753.3火灾模型喷射火计算加压气体泄漏时形成射流，如果在裂口处被点燃，形成喷射火。根据射流的方向，可以分成垂直喷射火和水平喷射火。根据?化工企业定量风险评价导则?【5】AQT3046-2021第条泄漏方向应根据设备安装的实际情况确定。如果没有准确的信息，泄漏方向宜设为水平方向，与风向一样。所以该企业液氨储罐的泄漏方向宜设为水平方向，仅水平方向喷射火计算。加压的可燃物泄漏时形成射流，如果在泄漏裂口处被点燃，则形成喷射火。假定火焰为圆锥形，并用从泄漏处到火焰长度4/5处的点源模型来表示。〔1〕火焰长度计算喷射火的火焰长度可用如下方程得到：〔3.12〕式中：L——火焰长度，单位为m；HC——燃烧热，液氨的燃烧热为11053.315kJ/kg；m——质量流速，单位为kg/s。1〕泄漏场景为小孔泄漏时，泄漏孔径为5mm，质量流速795.7g/s。=0.3489m2〕泄漏场景为中孔泄漏时，泄漏孔径为25mm，质量流速m=19.89kg/s=1.456m3〕泄漏场景为大孔泄露时，泄漏孔径为50mm，质量流速m=79.57kg/s=2.696m(2)热辐射的通量计算距离火焰点源为*(m)处接收到的热辐射通量可用下式表示：(3.13)式中：q——距离*处接收的热辐射的通量，单位为kw/m2；f——热辐射率；τ——大气传输率。大气传输率τ按下式计算：(3.14)选择泄漏场景为完全破裂时，质量流速m=79.57kg/s，热辐射率f取0.25，*分别取10m、50m、100m、150m。不同距离热辐射通量见表3.3所示：表3.3不同距离热辐射通量*〔m〕1050100150τ0.870.780.740.72q〔kw/m²〕0.155.459×10-31.295×10-35.60×10-43.4爆炸模型计算3.4.1沸腾液体扩展为蒸气云爆炸〔BLEVE〕计算沸腾液体扩展为蒸气云爆炸〔BLEVE〕计算:采用国际劳工组织建议的沸腾液体扩展为蒸气云爆炸热辐射模型进展计算，步骤如下：(1)火球直径的计算火球直径计算公式为：(3.15)式中：R——火球直径，单位为m；W——火球中消耗的可燃物质量，单位为kg，对于单罐储存，取罐容量的50%，对于双罐储存，取罐容量的70%；对于多罐储存，取罐容量的90%。该企业属于双液氨储罐，容积都为90m3,每个储罐的实际储存量为47.25T。即三个储罐的实际总储存量为2×47.25=94.5T所以W=70%M=70%×94.5×10^3=6.62×10^4kg即=2.9×〔6.62×10^4〕1/3=117.31m(2)火球持续时间的计算火球持续时间按下式计算：〔3.16〕式中：t——火球持续时间，单位为s；W——火球消耗的可燃物质量，单位为kg。即=18.2s(3)目标接收到热辐射通量的计算〔3.17〕式中：q0——火球外表的辐射通量，单位为W/m2；对于柱形罐取270W/m2，对于球形罐取200W/m2；r——目标到火球中心的平均距离，单位为m。q0取270W/m2，分别取r=10m、15m、20m、50m、100m，来进展目标接收到热辐射通量的计算，具体数据见下表3.1所示：当r=10m时，=19.73W/m2当r=15m时，=29.19W/m2当r=20m时，=38.21W/m2当r=50m时，=77.62W/m2当r=100m时，=87.89W/m2各点目标承受的热辐射通量见表3.4所示：表3.4目标接收的热辐射通量r(m)10152050100qr(W/m2)19.7329.1938.2177.6287.893.4.2蒸气云爆炸〔TNT当量法〕计算当泄漏到空气中的可燃气体与空气的混合物的浓度处于爆炸*围内时遇到点火源发生的爆炸现象称为蒸气云爆炸。其主要危害为爆炸产生的冲击波，能导致人员伤亡及设备、设施、建筑、厂房的破坏。它是一类经常发生、且后果十分严重的爆炸事故。通常采用TNT当量法估计蒸气云爆炸的严重度【8】。WTNT=αWQ/QTNT〔3.18〕式中WTNT—可燃气体的TNT当量，kg；α—可燃气体蒸气云当量系数(一般取值0.04)；W—蒸气云中可燃气体质量，kg；Q—可燃气体的燃烧热，J/kg；QTNT—TNT的燃烧热，J/kg(一般为4.52×106)；E=1.8αWQ〔3.19〕式中E—可燃气体的爆炸总能量，J；1.8—地面爆炸系数。爆炸伤害半径：〔3.20〕式中：C为爆炸实验常数，取决于损害等级，取值为0.03~0.4；N为有限空间内爆炸发生系数，取10%液氨30℃时，单罐实际储存量为47.25T【1】，即两个罐的总储存量为94.5T,氨的分子量为17，燃烧热18.59MJ/kg，沸点-33℃，液体比热容4.6kJ/kg·℃，气化热q为1.37×103kJ/kg。(1)氨气的TNT当量计算WTNT的计算WTNT=αWQ/QTNT=0.04×94.5×10³×18.59×10^6/4.52×10^6=1.5547*10^4kg(2)氨气爆炸总能量E的计算根据公式3.15氨气爆炸总能量：E=1.8αWQ＝1.8×0.04×94.5×10³×18.59=1.26×105MJ(3)氨气蒸气云爆炸冲击波的损害半径根据公式3.16计算氨气蒸气云爆炸冲击波的损害半径为：R=C(NE)1/3=0.2×〔10%×1.26×10^11〕1/3=465.39m〔4〕死亡区伤害半径死亡区内的人员如缺少防护则被认为将无例外地受到严重伤害或死亡，其内径为零，外径为R1,，其计算公式为:式中，WTNT为蒸气云的TNT当量，kg.R1=35.54〔5〕重伤区、轻伤区伤害半径重伤区的人员如缺少防护，则绝大多数人员将遭受严重伤害，极少数人可能死亡或受轻伤。其内径就是死亡半径R1，外径为R2，其计算公式为:式中，E为爆炸产生的总能量，J;P0为环境压力，Pa，取Po=101300Pa,△P为引起重伤冲击波超压峰值，Pa，取△P=44000Pa.轻伤区内的人员如缺少防护，则绝大多数人员将受轻微伤害，少数人将受重伤或平安无事，死亡的可能性极小。内径为R2，其外径为R3，其计算公式与重伤区伤害半径计算公式一样，其中△P为引起轻伤冲击波超压峰值，Pa，取△P=17000Pa.计算得出重重伤区z=1.089，轻伤区z=1.957R2=117.12R3=210.46-.z3.5事故伤害的计算3.5.1火灾热辐射伤害计算火灾通关辐射热的方式影响周边环境，当火灾产生的热辐射足够大时，可使周边的物体燃烧或者变形，强烈的热辐射可能烧毁设备甚至造成人员伤亡。人身伤害半径计算对于喷射火这类稳态火灾，其人身伤害半径的计算是：根据暴露时间10s，死亡时概率为50%的情况下求得导致死亡的热通量q1、重伤热通量q2、轻伤热通量q3分别为81830、69522、305