氨气气体泄漏后果分析

摘要
本次课程设计是对液氨储罐事故后果进行定量分析。

根据液氨的理化性质和主要危险特性, 釆用事件树分析法对液氨储罐的事故类型进行了危险源辨识和分析。

并通过对泄漏、火灾、爆炸等典型事故影响模型进行了定量计算，分析了液氨泄漏的危害范围及其严重程度，提出了科学、合理、可行的安全对策措施和建议，从而为预防液氨泄漏事故发生和进行安全管理提供了依据。

关键词：液氨储罐：危险源：事故影响模型
1、储罐区的情况简介
1.1储罐区的基本情况
本次课程设计以沈阳某发电厂为研究对象。

该发电厂采用选择性催化还原法进行烟气脱硝，所用还原剂为液氨，共有2个液氨储罐,每个液氨储罐最大储存量为90 m3。

液氨储存温度为30°C,储存压力为1. 1 MPa, 密度为750 kg/m3,液氨体积占储罐容积的最大值为70%（其充装系数为0.70）。

则每个贮槽内液氨的总质疑为W=90 m3 X 750 kg/m3 X 0. 7=47.25t。

重大危险源，是指长期地或者临时地生产、搬运、使用或者储存危险物品、且危险物品的数量等于或者超过临界量的单元（包括产所和设施）。

《危险化学品重大危险源辨识》【2】（GB18218-2009 ）规泄氨的临界量为20t,该企业布宜有两个液氨储罐，每个储罐装存的液氨为47.25t,因此构成了该液氨储罐区构成了重大危险源。

1.2液氨的主要危险特性
液氨又称为无水氨，是一种无色液体。

在温度变化时，液氨体积变化的系数很大。

溶于水、乙醇和乙醍，与空气混合能够形成爆炸混合物，火灾危险类别为乙类2项。

液氨作为一种重要的化工原料应用广泛，普遍存在于化工生产过程中。

为了运输及储存的便利，通常将气态的氨气通过加压或冷却得到液态氨。

液氨在工业上应用广泛，由于具有腐蚀性，且易挥发，所以苴化学事故发生率相当高，是该储罐区的主要危险物料。

液氨物料的危险特性主要体现在燃烧和爆炸、活性反应和健康危害三方面【3】，具体危险特性及理化性质如表:L1所示：
表1.1
2、储罐的事故类型分析
2.1 液氨泄漏事故模式及统计分析通常情况下，液氨在常温下加压压缩，液化储存，一目泄漏到空气中会在常压下迅速膨胀，大量气化,并扩散到大的空间范围。

液氨事故泄漏后通常有几种事故模式：液氨泄漏后在泄漏出口处立即点火形成喷射火：泄漏处于开放空间且经过一定时间点火形成闪火：泄漏处于局限空间条件且经过一左时间点火形成蒸气云爆炸；若泄漏过程中没有点火源存在则形成单纯的大气中扩散：储存液氨的储罐有可能发生BLEVE爆炸。

根据事故案例，泄漏事故绝大多数可视为连续点源泄漏模式。

因此，本文着重分析和模拟氨气连续点源泄漏后果。

2.2气氨泄漏事故树模型
对于液氨储罐液体泄漏应考虑发生池火+喷射火/BLEVE、蒸汽爆炸、闪火+晚期池火等火灾爆炸场景。

具体场景与泄漏类型、点火类型、空间局限情况等有关，可采用事件树方法确定务种可燃物质释放后，各种事件发生的类型及概率[5] o
1.1
定爆炸
是
否闪火
液氨泄:漏否________
否大气扩散
池火（+喷射火/BLEVE）是
爆炸
闪火（池
否大气扩散
图1丄可燃液体泄露事件树
3氨气泄漏后果左量计算
3.1气氨泄漏模型
对于灾难性破坏引起的液氨泄漏，可保守地认为容器内所有的贮存物质瞬间全部泄漏，全部泄漏时一般有爆炸发生，对其发生爆炸后的状况再运用数值模拟进行预测意义不大。

因此，文中所研究的是液姐储罐连续性泄漏的数值模拟。

通过对建国50年以来我国化工系统所发生的重(特)大、典型事故性泄漏的统计分析表明，阀门或法兰处的密封失效及阀门或管道断裂是造成事故性泄漏的主要原因，因而可以确定液氨储罐下方的液氨岀口接管、储罐上方的气氨出口接管以及安全阀为主要泄漏源。

气氨泄漏模型对于气氨通过其出口接管的泄漏，情况较复杂。

由于出口接管处于气相空间，英泄漏形式主要与泄
漏而积的大小有关。

在泄漏而积较大的情况下，高压蒸气通过裂缝或孔洞喷出，储罐内压急剧下降，直到环境压力(常温)。

由于内压急剧下降，气液平衡遭到破坏，储罐内液氨处于过热状态，过热状态的液氨为了再次恢复平衡，内部会均匀地产生沸腾核，同时产生大量气泡,液氨体积急剧膨胀，最终导致爆炸;当泄漏而积不大时，即使有蒸气喷出,但由于储罐内压下降不急剧，液氨不会达到过热状态，因此不会发生蒸气爆炸
3.1.1液体泄漏计算
由于储罐发生破裂泄漏的概率很小，因此泄漏场景考虑接管泄漏。

根据《钢制立式圆筒固左顶储罐系列》【8】HG 21502.1-1992,确定接管管径为50mm。

液体经管道上的孔流出的质量流率为：
Q m=AC^2p(P-P.) (31)
式中：
Qm——质量流率，单位为kg/s；
A——泄漏孔而积，单位为m2：
C0——液体泄漏系数；
P—一管道内液体压力，单位为Pa：
P—一泄漏液体密度，单位为kg/m3；
P0——环境压力，单位为Pa。

液氨的密度为0.75g/cm3 ,即P =750kg/m3 »储罐内液体压力为1.1 X 106Pa (表压力)，环境压力P0=1.013X105Pa,即储罐内液体压力(绝对压力)P=l.l X 106+1.013 X 105=1.2 X 10A6 Pa,当液体泄漏系数不知道或不能确定时，取1.0,使所计算的流量最大。

1) 泄漏场景为小孔泄漏时，泄漏孔径为5mm,即半径为R=2.5mm-J
A= nR2 =1.96X10-5 m2
fin =ACj2p(P-P.)
=(1.96X10-6) X1X[2X75OX (1.2X10A6-1.013X 105)]l/2
=1.106kg/s=0.7957kg/s
2) 泄漏场景为中孔泄漏时，泄漏孔径为10mm，即半径为R=12.5mm。

A=nR2 =4.9X10-4 m2
=4.9X 10-4X1X[2X75OX (1.2X 10A6 -1.013X105)] 1/2
=19.89kg/s
3) 泄漏场景为大孔破裂时，泄漏孔径为50mm,即半径为R=25mm。

A= nR2 =1.96X10-3 m2
Q m=AC^2p(P-P^
=1.96X10-3X1X[2X75OX (1.2X 10A6 -1.013X 105))1/2 =79.57kg/s
3.1.2闪蒸计算
(1)液体中闪蒸部分：
F7；)
一H、・(3.2)
过热液体闪蒸蒸发速率可按下式计算：
Q =(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=i366.70kJ/kg，定压热容Cp=4.609
kJ/(kg ・ K)
F _C p(J'r-T b)
' H- =0.21414
由于Fv>0.2时，英泄露的液氨完全闪蒸：/
3.2气氨扩散模型
根据泄漏气体的密度和泄漏源类型，气体的扩散模式可分为烟团扩散和烟羽扩散两种模式。

泄漏量较大且密度比空气的密度大得多的气体扩散呈现烟团式扩散，若泄漏气体密度与空气接近或经很短时间的空气稀释后密度与空气接近时，其泄漏气体的扩散属于烟羽式扩散，大部分较小流量的连续源易形成烟羽扩散。

因此，氨气的扩散属于烟羽扩散。

通过査询沈阳地区的气象资料，沈阳一年内的平均风速为2.84m/s,根据《化工企业左量风险评价导则》AQ/T3046-2013表E.5可知，大气稳定度为C级，根据表E.7可知：
o y=0.22x (l+0.0004x) A(-l/2) o z= 0.20x
根据Pasquill-Gifford模型扩散方程，位于地面比髙处的连续稳态源的烟羽
(3.6)
式中:
（Q（x』,z）——连续排放时，形成稳泄的流场后，给立地点（兀”2）的污染物的浓度.单位为g/m3：
°——连续排放的物料质量流量，单位为g/s；
U ------ 风速，单位为m/s；
by,叭一一侧风向和垂宜风向的扩散系数，单位为m：
x——下风向距离，单位为m；
y一一侧风向距离，单位为m：
z——垂直风向距离，单位为m。

C（工9 % s 0）=―——
由液体泄漏中的汁算可知，当泄漏为小孔泄漏时，连续排放的物料质量流量Q=795.7g/s。

因为液氨储罐水平布這在地而上，所以可以在x轴方向取10m、50m. 100m200m. 400m, y轴取0, z轴取0,具体地点污染物浓度见表3.1所示：
表3.1不同地点液氨扩散的浓度
3.3火灾模型
喷射火计算
加压气体泄漏时形成射流，如果在裂口处被点燃，形成喷射火。

根据射流的方向，可以分成垂直喷射火和水平喷射火。

根据《化工企业怎量风险评价导则》（5] AQT 3046-2013第9.2.3条泄漏方向应根据设备安装的实际情况确圧。

如果没有准确的信息，泄漏方向宜设为水平方向，与风向相同。

所以该企业液氨储罐的泄漏方向宜设为水平方向，仅水平方向喷射火计算。

加压的可燃物泄漏时形成射流，如果在泄漏裂口处被点燃，则形成喷射火。

假泄火焰为圆锥形，并用从泄漏处到火焰长度崛处的点源模型来表示。

（1）火焰长度计算
喷射火的火焰长度可用如下方程得到：
161.66 (3.12)
式中：
L——火焰长度，单位为m：
HC——燃烧热，液氨的燃烧热为11053.315 kJ/kg：
m ------ 质疑流速，单位为kg/s。

1) 泄漏场景为小孔泄漏时，泄漏孔径为5mm,质量流速795.7g/s。

161.66 =0.3489m
2) 泄漏场景为中孔泄漏时，泄漏孔径为25mm,质量流速m=19.89kg/s
厶=些型二
161 66=1.456m
3) 泄漏场景为大孔泄露时，泄漏孔径为50mm,质量流速m=79.57kg/s
L =
161.66 =2.696m
(2)热辐射的通星计算
距离火焰点源为X(m)处接收到的热辐射通量可用下式表示：
_ JH c m T
q~ ^7tX~ xlOOO (3.13)
式中：
q——距离X处接收的热辐射的通量，单位为kw/m2；
f—一热辐射率；
T——大气传输率。

大气传输率T按下式计算：
T = 1 —0.0565In X⑶口)
选择泄漏场景为完全破裂时，质量流速m=79.57kg/s,热辐射率f取0.25, X分別取10m、50m. 100m、150m.不冋距离热辐射通量见表3・3所示：
33
34爆炸模型计算
3.4.2沸腾液体扩展为蒸气云爆炸(BLEVE)计算
161.66 (3.12)
沸腾液体扩展为蒸气云爆炸(BLEVE)计算:采用国际劳工组织建议的沸腾液体扩展为蒸气云爆炸热辐射模型进行计算，步骤如下：
当r=15m 时,
=29.19 W/m2
⑴火球直径的讣算 火球直径il •算公式为：
R = 2.9l 沙
(3.15)
式中：
R —一火球直径，单位为m ：
W ——火球中消耗的可燃物质量，单位为kg,对于单罐储存，W 取罐容量的50%,对于双罐储存，W 取
罐容量的70%；对于多罐储存，“取罐容量的go%.
已知该企业属于双液氨储罐，容积都为90m3,每个储罐的实际储存量为47.25T 。

即三个储罐的实际总储存 量为
2X47.25=94.5T
所以 W=70%M=70% X 94.5 X 10A 3=6.62 X10八4kg 即/? = 2.9W.9X (6.62X10M) V3=117.31m ⑵火球持续时间的计算 火球持续时间按下式计算：
心0・451沙 (3.16)
式中：
t ——火球持续时间，单位为S ：
W ——火球消耗的可燃物质量，单位为kg 。

即 Z45W —8.2s
⑶目标接收到热辐射通量的计算
式中：
q0——火球表而的辐射通量，单位为W/m2：对于柱形罐取270 W/m2 ,对于球形罐取200 W/m2 : r —一目标到火球中心的平均距离，单位为m 。

q0取270 W/m2,分别取r=10m. 15m. 20m 、50m 、100m,来进行目标接收到热辐射通量的讣算，具体 数据见下表34
所示：
t/o /e 2r(l-O.O581nr)
t/o /?2r(l-0.0581nr)
g o /?2r(l-0.0581nr)
(宀尸严
(3.17)
当r=10m 时,
(宀打
=19.73 W/m2
（）_ （[一o.O581n / ） 当r=20m 时， （炉+尸）一
=38.21 W/m2 (.(y o /e 2r(l-O.O58In r) 当「沁时， (2严 =77.62W/m2 qO )_%)/?b(i_0・0581n r) 当 r=100m 时， （炉+尸严
=87.89W/m2 各点目标接受的热辐射通量见表 3.4所示: 表3・4目标接收的热辐射通量
342蒸气云爆炸（TNT 当量法）计算 当泄漏到空气中的可燃气体与空气的混合物的浓度处于爆炸范弗I 内时遇到点火源发生的爆炸现象称为蒸气 云爆炸。

其主要危害为爆炸产生的冲击波，能导致人员伤亡及设备、设施、建筑、厂房的破坏。

它是一类 经常发生、且后果十分严重的爆炸事故。

通常采用TNT 当量法估计蒸气云爆炸的严重度【8】。

WTNT= a WQ/ QTNT (3.18)
式 中 WTNT — 可 燃气 体 的 TNT 当 量 ， kg
a —
可燃 气 体蒸气
云当
系数
（一 般取值 0.04)
w —
气
云 中 可 燃 气 体 质
kg
Q - - 可
燃 气
体 的
燃
烧
热
，
J/kg
QTNT - 一 TNT 的燃
烧 热 ，
J/kg(
一 般为 4.52
X
106)
E=1.8 J WQ
(3.19
) 式 中 E — - 可
燃
气 体 的
爆
炸 总 能 量 , J
• ■ 1.8
—
地
面 爆
炸
系
数
o 爆炸伤 害半径 ■ ■ R = • • C(NE),/3
• • MM • •
• •
(3.20
) 式中: C 为 爆 炸实 验常 数，取 决于 损害 等级， 取值为 0.03P4 ■
N 为有限空间内爆炸发生系数，取10% 液氨30 °C 时，单罐实际储存量为47.25T [1],即两个罐的总储存量为94.5T,氨的分子量为17, 燃烧热 18.59
MJ/kg,沸点-33 *C,液体比热容 4.6 kJ/kg • *C ,气化热 q 为 1.37X103 kJ/ kg. ⑴氨气的TNT 当量计算WTNT 的计算 WTNT= a WQ/ QTNT=0.04 X 94.5 X 103 X 18.59 X 10A ^4.52 X 10A 6=1.5547*10A 4kg ⑵氨气爆炸总能量E 的计算 根据公式3.15氨气爆炸总能量： E=1.8 a WQ =1.8 X 0.04 X 94.5 X IO 3 X 18.59
=1.26X105MJ
(3)氨气蒸气云爆炸冲击波的损害半径
根据公式3.16计算氨气蒸气云爆炸冲击波的损害半径为:
R=C(NE)1^=O.2X (10%X1.26X10A ll) l/3=465.39m
(4)死亡区伤害半径
死亡区内的人员如缺少防护则被认为将无
例外地受到严重伤害或死亡，其内径为零，外径为R1,,其计算公式为：
RR3.6(W 刑丁 / 1000)删
式中，WTNT为蒸气云的TNT当量，kg.
Rl=35.54
(5)重伤区、轻伤区伤害半径
重伤区的人员如缺少防护，则绝大多数人员将遭受严重伤害，极少数人可能死亡或受轻伤。

苴内径就是死亡半径R1，外径为R2,其计算公式为：
AP/P0 = 0.137Z-3+0.119Z-2+O.269Z-' -0.019
Z = R/(E/P0),/3
式中，E为爆炸产生的总能量，J;PO为环境压力，Pa,取Po = 101300 Pa, AP为引起重伤冲击波超压峰值，Pa,取△P=44000Pa.
轻伤区内的人员如缺少防护，则绝大多数人员将受轻微伤害，少数人将受重伤或平安无事，死亡的可能性极小。

内径为R2,其外径为R3,其计•算公式与重伤区伤害半径计算公式相同，其中AP为引起轻伤冲击波超压峰值，Pa,取△P=17000Pa.
计算得出重
重伤区z=1.089,轻伤区z=1.957
R2=117.12
R3=210.46
3.5事故伤害的计算
3.5.1火灾热辐射伤害计算
火灾通关辐射热的方式影响周边环境，当火灾产生的热辐射足够大时，可使周边的物体燃烧或者变形，强烈的热辐射可能烧毁设备甚至造成人员伤亡。

人身伤害半径计算
对于喷射火这类稳态火灾，苴人身伤害半径的计算是：根据暴露时间10s,死亡时概率为50%的情况下求得导致死亡的热通虽ql.重伤热通量q2、轻伤热通量q3分别为81830、69522、30548W/m2.然后利用对应的热通量qi来反算得到人身死亡、重伤、和轻伤半径。

选择泄漏场景为大孔泄漏时，质量流速m=79.57kg/s,热辐射率f取0.25,
q=如"
r= 1-0.0565In X
Rl=
R2=
R3=
3.5.2毒物泄漏伤害计算
有毒物质泄漏后生成有毒蒸气云，它在空气中飘移、扩散，直接影响现场人员并波及居民区。

毒物对人员的危害程度取决于毒物的性质、浓度和人员与毒物接触的时间等因素。

毒性址露下死亡概率值可按下式讣算：
£ I =a + bh\(C,t x/)
(3.19 )
C* =2.24X10A7C/M 式中：
'傑——毒性暴露下的伤亡概率单位值
C* ------ 毒气的体积浓度，单位为mg/m3：
C ------- 浓度，单位为kg/m3 :
M——泄漏气体的相对分子量：
a,b,n——描述物质毒性的常数
t——眾露于毒物环境中的时间，单位为min,最大值为30 min。

不同距离毒性暴壺下的伤亡概率见表3.5所示：
表3.5毒性疑露下的伤亡概率
4、结论本次课程设计对液氨储罐事故后果进行了左量分析，并分析了液氨泄漏的危害范帀及其严重程度，提岀了科学、合理、可行的安全对策措施和建议。

经过这次课程设计认识到了在化工及相关行业中进行易燃、
易爆及毒性物质泄漏事故的理论及实验研究，对促进泄量化风险评价方法的成熟和为企业及政府部门采取安全防护措施和制左事故应急救援预案提供科学依据的重要性。

并且对液氨储罐存在的危险有了相应的认识，掌握了相关危险的控制手段。

此外，通过以上对液氨储罐事故后果的立量分析，可见液氨是一种极英危险的化学品，一旦发生泄漏，影响范围大，危害性强，因此必须采取一左的安全措施防止事故的发生。