1、摘要本次课程设计是对液氨储罐事故后果进行定量分析。根据液氨的理化性质和主要危险特性, 采用事件树分析法对液氨储罐的事故类型进行了危险源辨识和分析。并通过对泄漏、火灾、爆炸等典型事故影响模型进行了定量计算，分析了液氨泄漏的危害范围及其严重程度,提出了科学、合理、可行的安全对策措施和建议，从而为预防液氨泄漏事故发生和进行安全管理提供了依据。关键词：液氨储罐；危险源；事故影响模型1、储罐区的情况简介1.1 储罐区的基本情况本次课程设计以沈阳某发电厂为研究对象。该发电厂采用选择性催化还原法进行烟气脱硝,所用还原剂为液氨, 共有 2 个液氨储罐, 每个液氨储罐最大储存量为 90 m3。液氨储存温度为 3
2、0, 储存压力为 1. 1 MPa, 密度为 750 kg/m3,液氨体积占储罐容积的最大值为 70%(其充装系数为 0.70)。则每个贮槽内液氨的总质量为 W=90 m3750 kg/m30. 7=47.25t。重大危险源，是指长期地或者临时地生产、搬运、使用或者储存危险物品、且危险物品的数量等于或者超过临界量的单元（包括产所和设施） 。 危险化学品重大危险源辨识 【2】 (GB182182009）规定氨的临界量为 10t，该企业布置有两个液氨储罐，每个储罐装存的液氨为 47.25t，因此构成了该液氨储罐区构成了重大危险源。1.2 液氨的主要危险特性液氨又称为无水氨，是一种无色液体。在温度变
3、化时，液氨体积变化的系数很大。溶于水、乙醇和乙醚，与空气混合能够形成爆炸混合物，火灾危险类别为乙类 2 项。液氨作为一种重要的化工原料应用广泛，普遍存在于化工生产过程中。为了运输及储存的便利，通常将气态的氨气通过加压或冷却得到液态氨。液氨在工业上应用广泛，由于具有腐蚀性，且易挥发，所以其化学事故发生率相当高，是该储罐区的主要危险物料。液氨物料的危险特性主要体现在燃烧和爆炸、活性反应和健康危害三方面【3】 ，具体危险特性及理化性质如表 1.1 所示：表 1.1 液氨的危险特性及理化性质表燃烧和爆炸危险性极易燃，能与空气形成爆炸性混合物，遇明火、高热引起燃烧爆炸。活性反应 与氟、氯等接触会发生剧烈
4、的化学反应。危险特性健康危害对眼、呼吸道粘膜有强烈刺激和腐蚀作用。急性氨中毒引起眼和呼吸道刺激症状，支气管炎或支气管周围炎，肺炎，重度中毒者可发生中毒性肺水肿。高浓度氨可引起反射性呼吸和心搏停止。可致眼和皮肤灼伤。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，燃点为 6
5、51, 最大爆炸压力 0.580Mpa2、储罐的事故类型分析2.1 液氨泄漏事故模式及统计分析通常情况下,液氨在常温下加压压缩, 液化储存,一旦泄漏到空气中会在常压下迅速膨胀, 大量气化,并扩散到大的空间范围。液氨事故泄漏后通常有几种事故模式：液氨泄漏后在泄漏出口处立即点火形成喷射火；泄漏处于开放空间且经过一定时间点火形成闪火；泄漏处于局限空间条件且经过一定时间点火形成蒸气云爆炸；若泄漏过程中没有点火源存在则形成单纯的大气中扩散；储存液氨的储罐有可能发生 BLEVE 爆炸。根据事故案例 , 泄漏事故绝大多数可视为连续点源泄漏模式。因此,本文着重分析和模拟氨气连续点源泄漏后果。2.2 气氨泄漏事
6、故树模型对于液氨储罐液体泄漏应考虑发生池火+喷射火/BLEVE 、蒸汽爆炸、闪火+晚期池火等火灾爆炸场景。具体场景与泄漏类型、点火类型、空间局限情况等有关，可采用事件树方法确定各种可燃物质释放后，各种事件发生的类型及概率【5】 。可燃液体泄露后的事件树见图 1.1 所示：IE 完全闪蒸 立即点火 延迟点火 空间局限 后果是 喷射火是 是 爆炸 是否 闪火液氨泄漏 否 否 大气扩散是 池火（+喷射火/BLEVE）是否 爆炸是否 否 闪火（+晚期池火）否 大气扩散 图 1.1 可燃液体泄露事件树3 氨气泄漏后果定量计算3.1 气氨泄漏模型对于灾难性破坏引起的液氨泄漏,可保守地认为容器内所有的贮存物
7、质瞬间全部泄漏,全部泄漏时一般有爆炸发生,对其发生爆炸后的状况再运用数值模拟进行预测意义不大。因此,文中所研究的是液氨储罐连续性泄漏的数值模拟。 通过对建国 50 年以来我国化工系统所发生的重(特) 大、典型事故性泄漏的统计分析表明 ,阀门或法兰处的密封失效及阀门或管道断裂是造成事故性泄漏的主要原因,因而可以确定液氨储罐下方的液氨出口接管、储罐上方的气氨出口接管以及安全阀为主要泄漏源。气氨泄漏模型对于气氨通过其出口接管的泄漏,情况较复杂。由于出口接管处于气相空间,其泄漏形式主要与泄漏面积的大小有关。在泄漏面积较大的情况下,高压蒸气通过裂缝或孔洞喷出,储罐内压急剧下降,直到环境压力(常温)。由于
8、内压急剧下降, 气液平衡遭到破坏,储罐内液氨处于过热状态,过热状态的液氨为了再次恢复平衡,内部会均匀地产生沸腾核, 同 时产生大量气泡 ,液氨体积急剧膨胀, 最终导致爆炸;当泄漏面积不大时,即使有蒸气喷出 ,但由于储罐内压下降不急剧 ,液氨不会达到过热状态, 因此不会发生蒸气爆炸3.1.1 液体泄漏计算由于储罐发生破裂泄漏的概率很小，因此泄漏场景考虑接管泄漏。根据钢制立式圆筒固定顶储罐系列 【8】HG 21502.1-1992，确定接管管径为 50mm。 液体经管道上的孔流出的质量流率为：002PACQm （3.1）式中:Qm质量流率，单位为 kg/s；A泄漏孔面积，单位为 m2；C0 液体泄
9、漏系数；P管道内液体压力，单位为 Pa；泄漏液体密度，单位为 kg/m3；P0环境压力，单位为 Pa。液氨的密度为 0.75g/cm，即 =750kg/m。储罐内液体压力为 1.1106Pa（表压力） ，环境压力P0=1.013105Pa，即储罐内液体压力（绝对压力）P=1.1106+1.013105=1.2 106 Pa，当液体泄漏系数不知道或不能确定时，取 1.0，使所计算的流量最大。1）泄漏场景为小孔泄漏时，泄漏孔径为 5mm，即半径为 R=2.5mm。A=R=1.9610-5 m2002PACQm=（ 1.9610-6）12750（1.2 106 -1.013105)1/2=1.106
10、kg/s=0.7957kg/s2）泄漏场景为中孔泄漏时，泄漏孔径为 10mm，即半径为 R=12.5mm。A=R=4.910-4 m2002PACQm=4.910-412750（1.2 106 -1.013105)1/2=19.89kg/s3）泄漏场景为大孔破裂时，泄漏孔径为 50mm，即半径为 R=25mm。A=R= 1.9610-3 m2002PACQm=1.9610-312750（1.2 106 -1.013105)1/2=79.57kg/s3.1.2 闪蒸计算（1 ）液体中闪蒸部分：（3.2 ）vbTpvHCF过热液体闪蒸蒸发速率可按下式计算：（3.3）vLQ1式中：Fv泄漏液体的闪蒸
11、比例；TT储存温度，单位为 K；Tb泄漏液体的沸点，单位为 K；Hv泄漏液体的蒸发热，单位为 kJ/kg；Cp泄漏液体的定压热容，单位为 kJ/(kgK)；Q1过热液体闪蒸蒸发速率，单位为 kg/s；QL物质泄漏 速率，单位为 kg/s。已知储存温度 TT=30+273.15=303.15K，沸点 Tb=239.65K，泄漏液体的蒸发热 Hv=1366.70kJ/kg ，定压热容Cp=4.609 kJ/(kgK)=0.21414vbTpvHCF由于 Fv0.2 时，其泄露的液氨完全闪蒸：/3.2 气氨扩散模型根据泄漏气体的密度和泄漏源类型，气体的扩散模式可分为烟团扩散和烟羽扩散两种模式。泄漏量
12、较大且密度比空气的密度大得多的气体扩散呈现烟团式扩散，若泄漏气体密度与空气接近或经很短时间的空气稀释后密度与空气接近时，其泄漏气体的扩散属于烟羽式扩散，大部分较小流量的连续源易形成烟羽扩散。因此，氨气的扩散属于烟羽扩散。通过查询沈阳地区的气象资料，沈阳一年内的平均风速为 2.84m/s，根据化工企业定量风险评价导则AQ/T3046-2013 表 E.5 可知，大气稳定度为 C 级，根据表 E.7 可知：y=0.22x（1+0.0004x）(-1/2) z= 0.20x根据 Pasquill-Gifford 模型扩散方程，位于地面 rH高处的连续稳态源的烟羽 222 1exp1exp1exp2,
13、 zrzyzy HuQzxC （3.6）式中： zyxC,连续排放时，形成稳定的流场后，给定地点 zyx,的污染物的浓度，单位为 g/m3；Q连续排放的物料质量流量，单位为 g/s；u风速，单位为 m/s；y， z侧风向和垂直风向的扩散系数，单位为 m；x下风向距离，单位为 m；y侧风向距离，单位为 m；z垂直风向距离，单位为 m。由液体泄漏中的计算可知，当泄漏为小孔泄漏时，连续排放的物料质量流量 Q=795.7g/s。因为液氨储罐水平布置在地面上，所以可以在 x 轴方向取 10m、50m、100m、200m、400m，y 轴取 0，z 轴取 0，具体地点污染物浓度见表 3.1 所示：表 3.
14、1 不同地点液氨扩散的浓度X（m） 10 50 100 200 400 500Y（ m） 0 0 0 0 0 0Z（m） 0 0 0 0 0 0y（m）2.1956 10.07 18.98 35.76 67.38 82.62z（m）2 10.89 21.57 42.339 81.705 100.416C（ g/m3） 20.309 0.8132 0.2178 0.0589 0.0162 0.010753.3 火灾模型喷射火计算加压气体泄漏时形成射流，如果在裂口处被点燃，形成喷射火。根据射流的方向，可以分成垂直喷射火和水平喷射火。根据化工企业定量风险评价导则 【5】AQT 3046-2013 第
15、 9.2.3 条 泄漏方向应根据设备安装的实际情况确定。如果没有准确的信息，泄漏方向宜设为水平方向，与风向相同。所以该企业液氨储罐的泄漏方向宜设为水平方向，仅水平方向喷射火计算。加压的可燃物泄漏时形成射流，如果在泄漏裂口处被点燃，则形成喷射火。假定火焰为圆锥形，并用从泄漏处到火焰长度 4/5 处的点源模型来表示。（1）火焰长度计算喷射火的火焰长度可用如下方程得到：6.1)(40mHLC（3.12）式中：L火焰长度，单位为 m；HC燃烧热，液氨的燃烧热为 11053.315 kJ/kg；m质量流速，单位为 kg/s。1）泄漏场景为小孔泄漏时，泄漏孔径为 5mm，质量流速 795.7g/s。6.1
16、)(40mHLC=0.3489m2）泄漏场景为中孔泄漏时，泄漏孔径为 25mm，质量流速 m=19.89kg/s6.1)(40mHLC=1.456m3）泄漏场景为大孔泄露时，泄漏孔径为 50mm，质量流速 m=79.57kg/s6.1)(40mHLC=2.696m(2)热辐射的通量计算距离火焰点源为 X(m)处接收到的热辐射通量可用下式表示：(3. 13 )1042XmfHqC式中：q距离 X 处接收的热辐射的通量，单位为 kw/m2；f热辐射率；大气传输率。大气传输率 按下式计算：(3.14 ) Xln056.1选择泄漏场景为完全破裂时，质量流速 m=79.57kg/s，热辐射率 f 取 0
17、.25，X 分别取10m、50m、100m、150m 。不同距离热辐射通量见表 3.3 所示：表 3.3 不同距离热辐射通量X（m） 10 50 100 150 0.87 0.78 0.74 0.72q （kw/m ） 0.15 5.45910-3 1.29510-3 5.6010-43.4 爆炸模型计算3.4.1 沸腾液体扩展为蒸气云爆炸（BLEVE）计算沸腾液体扩展为蒸气云爆炸（BLEVE）计算:采用国际劳工组织建议的沸腾液体扩展为蒸气云爆炸热辐射模型进行计算，步骤如下：(1)火球直径的计算火球直径计算公式为：319.2WR (3.15)式中：R火球直径，单位为 m；W火球中消耗的可燃物质量，单位为 kg，对于单罐储存， W取罐容量的 50%，对于双罐储存， W取罐容量的 70%；对于多罐储存， W取罐容量的 90%。已知该企业属于双液氨储罐，容积都为 90m3,每个储罐的实际储存量为 47.25T。即三个储罐的实际总储存量为 247.25=94.5T所以 W=70%M=70%94.5103=6.62104kg即 319.WR=2.9（6.62104 ）1/3=117.31m(2)火球持续时间的计算
