爆炸评价模型伤害半径计算.docx
《爆炸评价模型伤害半径计算.docx》由会员分享，可在线阅读，更多相关《爆炸评价模型伤害半径计算.docx（22页珍藏版）》请在冰豆网上搜索。
爆炸评价模型伤害半径计算
爆炸评价模型及伤害半径计算
1、蒸气云爆炸（VCE）模型分析计算
（1）蒸气云爆炸（VCE）模型
当爆炸性气体储存在贮槽内，一旦泄漏，遇到延迟点火则可能发生蒸气云爆炸，如果遇不到火源，则将扩散并消失掉。
用TNT当量法来预测其爆炸严重度。
其原理是这样的：
假定一定百分比的蒸气云参与了爆炸，对形成冲击波有实际贡献，并以TNT当量来表示蒸气云爆炸的威力。
其公式如下：
WtNt=
BAWfQf
QTNT
式中WNt——蒸气云的TNT当量，kg;
B——地面爆炸系数，取6=1.8;
A——蒸气云的TNT当量系数，取值范围为0.02%〜14.9%;
Wf蒸气云中燃料的总质量：
kg;
Qf——燃料的燃烧热，kJ/kg;
Qtnt——TNT的爆热，QTNT=412O4690kJ/kg。
（2）水煤气储罐蒸气云爆炸（VCE分析计算
由于合成氨生产装置使用的原料水煤气为一氧化碳与氢气混合物，具有低闪点、低沸点、爆炸极限较宽、点火能量低等特点，一旦泄漏，极具蒸气云爆炸概率。
若水煤气储罐因泄漏遇明火发生蒸气云爆炸（VCE,设其贮量为70%寸，则
为2.81吨，则其TNT当量计算为：
取地面爆炸系数：
6=1.8；
蒸气云爆炸TNT当量系数，A=4%
蒸气云爆炸燃烧时燃烧掉的总质量，
Wf=2.81X1000=2810（kg）;
水煤气的爆热，以CO30%代43%+（氢为1427700kJ/kg,一氧化碳为10193
kJ/kg）:
取Q=616970kJ/kg;
TNT的爆热,取QN=4500kJ/kg。
将以上数据代入公式，得
1.8X0.04X2810X616970
45（2（7739（kg）
=13.6X27.74°.37
=13.6X3.42=46.5（m）
重伤半径R,由下列方程式求解：
厂
△P2=0.137Z2-3+0.119Z2—2+0.269Z21-0.019
1____1/3
Z2=R/（E/P0）v.
AP2=APS/Po
式中：
△R——引起人员重伤冲击波峰值，取44000P?
P0——环境压力（101300PO）;
E——爆炸总能量（J）,E=WNtXQtnto将以上数据代入方程式，解得：
△P2=0.4344
Z2=1.071/3R=1.07X（27739X4500X1000/101300）
=1.07X107=115（m）
轻伤半径R,由下列方程式求解：
r△R=0.137Z3-3+0.119Z3-2+0.269Z3-1-0.019
ZZ3=R/（E/P0）1/3
AP3=APS/Po
式中：
△Ps——引起人员轻伤冲击波峰值，取17000Pa
将以上数据代入方程式，解得：
△R=0.168,Z3=1.95
轻伤半径R=209（mj）
2、沸腾液体扩展蒸气爆炸（BLEVE）模型分析计算
（1）沸腾液体扩展蒸气爆炸（BLEVE模型
液态存贮的易燃液化气体突然瞬间泄漏时，立即遇到火源就会发生剧烈的燃烧，产生巨大的火球，形成强烈的热辐射，此种现象称为沸腾液体扩展蒸气爆炸，简称BLEVE
沸腾液体扩展蒸气爆炸的主要危险是强烈的热辐射，近场以外的压力效应不重要。
具火球的特征可用国际劳工组织（ILO）建议的蒸气爆炸模型来估算。
火球半径的计算公式为：
R=2.9W1/3
式中R——火球半径，m
W——火球中消耗白可燃物质量，kg。
对单罐储存，W取罐容量的50%双罐储存；W取罐容量的70%多罐储存,取W为罐容量的90%
（2）液氨储罐沸腾液体扩展蒸气爆炸（BLEVE模型分析计算
由于生产装置液氨贮罐区的液氨罐为多罐贮存，（共六只贮罐，其中三只
50M,三只100M）最大库存量为250T。
氨比重约0.6,取100M罐，则
由W=100<0.6X1000X90%=54000（kg）
代入式中，得到：
火球半径R=2.9（54000）1/3=109（m）
火球持续时间按下式计算：
t=0.45W1/3
式中:
火球持续时间，单位为S.
将数据代入式中，得到：
t=0.45X（54000）1/3=17（s）
目标接收到热辐射通量的计算，按下式计算：
q（r）=qoR2r（1-0.058Inr）/（R2+r2）3/2
式中：
r——目标到火球中心的水平距离，m
qo——火球表面白辐射通量，W/m。
对柱形罐取270kW/m,球形罐取
2
200kW/m
R火球半径，miR=109m
有了热辐射q（r）,即可求不同伤害、破坏时的热通量及其半径。
下面求不同伤害时的热通量：
死亡可根据下式计算：
Pr=-36.38+2.56In（tqi4/3）
式中：
Pr=5
t火球持续时间，取t=17s。
解得qi=21985W/rH
重伤可根据下式计算：
4/3
P「=-43.143+3.0188In（tq2）
2
解得q2=18693W/m
轻伤可根据下式计算：
Pr=-39.83+3.0188In（tq34/3）
2
解得q3=8207W/m
通过q〔、q2、q3可以求得对应的死亡半径R、重伤半径R及轻伤半径R。
（由于此方程式难以手算解出，故省略）。
（3）小结
通过计算，如果贮存区液氨储罐发生扩展蒸气爆炸，火球半径为109m将
可能造成其他贮罐的连锁火灾和爆炸，造成灾难性的破坏。
3、液氨泄漏中毒事故的模拟计算
液氨贮存区最大贮存量为250T,假设有仃泄漏量，对蒸发成蒸气扩散造成的危害进行模拟计算。
（1）液态气体蒸气体积膨胀计算
在标准状态下（0C,1013Mpa,1摩尔气体占有22.4升体积。
根据液态气体的相对密度，由下式可计算出它们气化后膨胀的体积：
V=V0D0100022.4
M
V——膨胀后的体积（升）
V0——液态气体的体积（升）
D0——液态气体的相对密度（水=1）
M——液态气体的的分子量
将液氨有关数据代入上式，由D0=0.597,M=17.03得到
V=V0D0100022.4M
即液态氨若发生泄漏迅速气化，其膨胀体积为原液态体积的785倍。
（2）液态气体扩散半径模拟计算
液态气体泄漏后在高温下迅速气化并扩散，在一定泄漏量范围内，且液态气体比重大于空气，沿地面能扩散到相当远的地方，可模拟为半椭圆形，具短轴与长轴之比将随着扩散半径的增大而减少，可由下式计算：
式中：
V——液态气体膨胀后体积;
液态气体在空气中的浓度；
椭圆形短轴与长轴之比，即K=h/R。
根据我国《工作场所有害因素职业接触限值》（GBZ2-2002查得：
液氨在工作场所空气中时间加权平均容许浓度20mg/m；短时间接触容许浓度30mg/m,其在空气中体积浓度换算为：
p=26.3X10-6和p=39.5X10-6O
假设泄漏液氨的量为1000kg,其可能发生中毒事故的浓度区域半径计算如下：
取液氨体积V0=1/0.597=1.68m3
p=26.3X10-6K=0.10
计算：
R-V
32、
V—xnxPxk
3
=1.68-785
32二26.31030.10
=621m
从计算结果可知：
当泄漏1000kg液氨气化成蒸气时可能发生中毒浓度的区域半径为621m即0.621公里，因此，其扩散时的可能发生中毒浓度的区域面积：
S=-XR2=3.14X0.6212=1.21（平方公里）
4、水煤气泄漏事故的模拟计算
根据我国《工作场所有害因素职业接触限值》（GBZ2-2002查得：
水煤气（即一氧化碳）时间加权平均容许浓度20mg/n3;短时间接触容许浓度30mg/mi。
经换算，分另1J为p=16X10-6和p=24X10-6。
水煤气贮罐总容积5000M3设若泄漏量为100M3取p=16X10-6,K=0.1则计算如下：
R=312V
—xnXPXK
二100
3|2^nM16M10上M0.10
=310（MD
S=兀xR=3.14X310=301754M
5、天然气泄漏形成喷射火模型分析
该项目设计有容积为2m3、工作压力为25MPa的天然气储气瓶组。
当设备损坏、法兰垫片撕裂或接管破裂等条件下，天然气就会在破裂处形成射流，在高速气流摩擦形成的静电火花或其他点火源存在的条件下，就会在裂口处引燃形成
喷射火。
火灾通过热辐射方式影响周围环境，当热辐射强度足够大时，可使周围
物体燃烧变形，强烈的热辐射可能烧毁设备并造成人员伤亡。
喷射火模拟事故模拟是通过定量的计算，估算出热辐射的不同入射通量所造成的损失程度。
（1）气体泄漏量的计算
假定天然气储气瓶组或管道泄漏，裂口为直径10mm的圆口：
P0=0.1013106Pa；P=25X06Pa；k=1.314天然气的绝热指数）
则：
F00.1013
P25
一==0.004
1.314
22'1.314」
।=0.54
11.314+1）
P1
故气体流动属于音速流动。
其泄漏量为：
Qo
M,
式中：
Cd=1.00;M=CNG的分子量,16;R=8.314J/mol?
K;T=293K（20C）;p=
193.5kg/m3,A为裂口面积，D为储气瓶裂口直径。
A"D2/40.0000785（m2）
计算得：
Q0=8.2X0-4kg/s
（2）喷射火热辐射通量
这里所用的喷射火辐射热计算方法是一种包括气流效应在内的喷射火扩散模式
的扩展。
把整个喷射火看成是沿喷射中心线上的几个热源点组成，每个点热源的
热辐射通量相等。
点热源的热辐射通量按下式计算
q=QoHc
式中：
q一点热源热辐射量，W
刀一^我率因子，可取0.35
Q0-泄漏速度，kg/s
Hl燃烧热，J/kg
天然气高热值：
Hc=55800kJ/kg,则：
q=0.358.2X0-4）55800=16kW
射流轴线上某点热源I到距离该处一点的热源辐射强度为：
「更
14二x2
式中：
Ii-点热源i到目标点x处的热辐射强度，W/m2；
q一点热源的辐射通量，W;
R一辐射率，可取0.2;
x一点热源到目标点的距离，m
某一点处的入射流强度等于喷射火的全部点热源对目标的热辐射强度的总和:
式中：
N一计算时选取的点热源系数，一般取n=5。
根据喷射火全部点热源在距火焰xm的某点总入射热辐射通量的大小，查热辐射的不同入射能量所造成的损失。
根据可以计算出有代表意义的入射热辐射通量造成的危害范围。
即：
qR
式中：
R—辐射率，取0.2;
表5.6-1不同入射通量所造成的危害范围
目标伤害跑离
x（m）
入射逋量I
（kW/m2）
对设备的损害
对人的损害
0.082
37.5
操作设备全部损坏
1%死亡/10s
100%死亡/1min
0.1
25
在无火焰、长时间辐射下，木材燃烧的最小能量
重大损伤/10s
100%死亡/1min
0.14
12.5
有火焰时，木材燃烧，
塑料融化的最低能量
1度烧伤/10s
1%死亡/1min
0.25
4.0
20s以上感觉疼痛，未必起泡
0.40
1.6
长时间辐射，无不舒服
感
从上表可以看出，如不采取措施，在计算条件下的喷射火如果发生，距射流
轴线热源0.082m处的人员在热辐射下10s内1%死亡，1min内100%死亡，另此
范围内的操作设备全部损坏。
距射流轴线热源0.25m,热辐射不会造成人员伤亡；
距射流轴线热源0.40m以外为安全距离。
为了安全生产，企业应积极采取应对措
施，防止喷射事故发生。
6、天然气泄漏形成蒸气云爆炸模型分析
假设1个2m3、操作压力25MPa的天然气储气瓶组全部发生泄漏，则参与爆炸的天然气体积为50m3,并与空气混合后形成蒸气云，蒸汽云爆炸冲击波的破坏半径为：
R=CS（NE）1/3
E一爆炸能量，kJ,E=VHC
V一参与反应的可燃气体的体积，m3;
HC—可燃气体的燃烧热值；天然气高热值：
39.86MJ/Nm3;
N—效率因子，一般取10%;
CS-经验常数，取决于损害等级。
其具体取值情况见下表：
表5.6-2损害等级表
损害等级
Cs/m?
J（1/3）
设备损害
人员伤害
1
0.03
重创建筑物的加工设备
1%死亡于肺部伤害；＞50%耳膜破裂；＞50%被碎片击伤
2
0.06
损害建筑物外表可修复性破坏
1%耳膜破裂；1%被碎片击伤
3
0.15
玻璃破碎
被碎玻王离击伤
4
0.4
10%玻璃破碎
蒸气云爆炸冲击波的破坏半径与损害等级关系，计算结果见下表:
表5.6-3天然气蒸气云爆炸冲击波损害范围表
损害等级
Cs
NE
（NE）1/3
损害半径
R
设备损坏
人员伤害
单位
m?
J（1/3）
J
J1/3
m
1
0.03
7.9721
06
199.77
5.99
重创建筑物和加工设备
1%死亡人员肺部伤害
＞50%耳膜破裂
＞50%被碎片击伤
2
0.06
11.98
损坏建筑物外表，可修复性破坏
1%耳膜破裂
1%被碎片击伤
3
0.15
29.97
玻璃破碎
被碎玻王离击伤
4
0.4
79.9
10%玻璃破碎
爆炸冲击波伤害半径的模拟分析表明，1只2m3、工作压力为25MPa的大然气储气瓶组全部发生泄漏后，天然气与空气形成爆炸性混合气，发生爆炸后在半径5.99m范围内，1%人员死亡及肺部伤害、＞50%耳膜破裂、＞50%被碎片击伤。
7、压缩气体与水蒸气容器爆破能量
当压力容器中介质为压缩气体，即以气态形式存在而发生物理爆
炸时，其释放的爆破能量为:
式中Eg——气体的爆破能量，kJ;
p——容器内气体的绝对压力，MPa;
V——容器白容积，m3;
K——气体的绝热指数，即气体的定压比热与定容比热之比常用气体的绝热指数数值如表1所示。
表1常用气体的绝热指数
气体名螂
空气
%
M
V烷
乙他
乙舞i
丙烷
i一■化碳
工值
14
1.4
1.397
厂
1.116
!
.18J
1.22
卜13-
气体名戟
二皿化段
氧化区
二战化气
a/
氯气
过热旅汽
施和蒸汽
r值］
1.295
一一
।E4
K31
E.12।
1.35
卜工3
1135
1.32
从表1可看出，空气、氮、氧、氢及一氧化氮、一氧化碳等气体的绝热指数均为1.4或近似1.4,如用k=1.4代入式
（1）中，得到
气体的爆破能量为:
r/ninntO-is5T*
令G=2.541-（匹詈二）J则式673可简化为:
Eg=CgV⑶
式中Cg——常用压缩气体爆破能量系数，kJ/m3。
压缩气体爆破能量系数Cg是压力p的函数，各种常用压力下的
气体爆破能量系数如表2所示
表2常用压力下的气体容器爆破能量系数（6=1.4时）
衰压Jjp/MP#
0.2I
04
06
OJt
to
L6
爆破能■系效
2xio1
4QXIO1
7.5XID1
i.Ixio*
J,4X10'
2.4M101
39xl0J
表压力p/MPa
---a——
4.0
-——1
5.0
64
L5.0
32
40
爆破能・系数
6,7X10J
86MIO-1
LIxID1
2.7m101
*
It）1
B.2x104
如将6=1.135代入式
（1）,可得干饱和蒸汽容器爆破能量为：
x10,
（（4）
用式4计算有较大的误差，因为没有考虑蒸汽干度的变化和其他一些影响，但可以不用查明蒸汽热力性质而直接计算，对危险性评价
可提供参考。
对于常用压力下的干饱和蒸汽容器的爆破能量可按下式计算：
Es=CsV（5）
式中Es水蒸气的爆破能量，kJ;
V——水蒸气的体积，m3;
Cs——干饱和水蒸气爆破能量系数，kJ/m3。
各种常用压力下的干饱和水蒸气容器爆破能量系数如表3所示。
表3常用压力下干饱和水蒸气容器爆破能量系数
森压力p/MFa
0Q
0.5
0.8
13
&5
3.0
爆里能■系数
4.37x10士
800
[-5x|0J
2.75xW*
6.24x103
7,77xio1
1.2介质全部为液体时的爆破能量
通常用液体加压时所做的功作为常温液体压力容器爆炸时释放的能量，计算公式如下：
p_（-T）2V自
匕L一子一
（6）
式中员——常温液体压力容器爆炸时释放的能量，kJ;
p液体的压力（绝），Pa;
V——容器白^体积，m3;
Bt——液体在压力卢和温度T下的压缩系数，Pai1。
1.3液化气体与高温饱和水的爆破能量
液化气体和高温饱和水一般在容器内以气液两态存在，当容器破
裂发生爆炸时，除了气体的急剧膨胀做功外，还有过热液体激烈的蒸发过程。
在大多数情况下，这类容器内的饱和液体占有容器介质重量的绝大部分，它的爆破能量比饱和气体大得多，一般计算时不考虑气体膨胀做的功。
过热状态下液体在容器破裂时释放出爆破能量可按下式计算：
E=[（Hi-H2）-（Si-S2）Ti]W（7）
式中E——过热状态液体的爆破能量，kJ;
Hi——爆炸前液化液体的始，kJ/kg;
H2——在大气压力下饱和液体的始，kJ/kg;
S1——爆炸前饱和液体的，嫡，kJ/（kg•C）；
S2——在大气压力下饱和液体的嫡，kJ/（kg•C）；
Ti——介质在大气压力下的沸点，C；
W——饱和液体时质量，kg。
饱和水容器的爆破能量按下式计算：
Ew=CwV（8）
式中Ew——饱和水容器的爆破能量，kJ;
V——容器内饱和水所占的容积，m3;
Cw——饱和水爆破能量系数，kJ/m3,其值如表4所示
表4常用压力下饱和水爆破能量系数
队3
IJ.5
08
1.3
25
30
能量系数
2乐Ml型
4.56x10.
635*10*
重56x10*
1一06k10*
8、冲击波超压的伤害-破坏作用
压力容器爆破时，爆破能量在向外释放时以冲击波能量、碎片能量和容器残余变形能量3种形式表现出来。
根据介绍，后二者所消耗的能量只占总爆破能量的3%〜15%,也就是说大部分能量是产生空气冲击波。
冲击波是由压缩波迭加形成的，是波阵面以突进形式在介质中传播的压缩波。
容器破裂时，容器内的高压气体大量冲出，使它周围的空气受到冲击而发生扰动，使其状态（压力、密度、温度等）发生突跃变化，其传播速度大于扰动介质的声速，这种扰动在空气中传播就成为冲击波。
在离爆破中心一定距离的地方，空气压力会随时间迅速发生而悬殊的变化。
开始时，压力突然升高，产生一个很大的正压力，接着又迅速衰减，在很短时间内正压降至负压。
如此反复循环数次，压力渐次衰减下去。
开始时产生的最大正压力即是冲击波波阵面上的超压Ap。
多数情况下，冲击波的伤害一破坏作用是由超压引起的。
超压△p可以达到数个甚至数十个大气压。
冲击波伤害一破坏作用准则有：
超压准则、冲量准则、超压一冲量准则等。
为了便于操作，下面仅介绍超压准则。
超压准则认为，只要冲击波超压达到一定值时，便会对目标造成一定的伤害或破坏。
超压波对人体的伤害和对建筑物的破坏作用如表5和表6所示。
表5冲击波超压对人体的伤害作用
物害作用1
图压△户小曲、
桥客作用
内酰产重搞伤或死亡
唠她制官损佛或骨折
>0.101
大部分人员死亡
表6冲击波超压对建筑物的破坏作用
破塔作用
融压
破坏作用
0.005-0.006
门值■•部分破碎
03~047
本建筑厂房房柱折麻,唐颦松动
0.006—0,015
受压面的门四破■大部分破障
0.07-0.10
玲♦陶埸
0.015-（1.02
■据摄嘉
配10〜0.20
防蕤制筋混最上破坏.小房%倒用
0.02-0.03
『塘裂S!
0,20-0.30
大型器架结构胭坏
0.04-0.05
・大裂缝，屋瓦捶下
2.2冲击波的超压
冲击波波阵面上的超压与产生冲击波的能量有关，同时也与距离
爆炸中心的远近有关。
冲击波的超压与爆炸中心距离的关系：
△pxRn（9）
式中Ap——冲击波波阵面上的超压，MPa;
R——距爆炸中心的距离，m;
n衰减系数。
衰减系数在空气中随着超压的大小而变化，在爆炸中心附近内为
2.5〜3;当超压在数个大气压以内时，n=2;小于1atm（0.1MPa）
时，n=1.5。
实验数据表明，不同数量的同类炸药发生爆炸时，如果距离爆炸中心的距离R之比与炸药量q三次方根之比相等，则所产生的冲击波超压相同，用公式表示如下：
旨=/^=。
网&户=4力0
如1一」
（10）
式中R——目标与爆久^中心距离，m;
R0——目标与基准爆炸中心的相当距离，m;
q0——基准爆炸能量，TNT,kg;
q——爆炸时产生冲击波所消耗的能量，TNT,kg;
△p目标处的超压，MPa;
△po基准目标处的超压，MPa;
0t——炸药爆炸试验的模拟比。
式10也可写成为：
Ap（R）=Apo（R/%）（11）
利用式11就可以根据某些已知药量的试验所测得的超压来确定在各种相应距离下任意药量爆炸时的超压。
表7是1000kgTNT炸药在空气中爆炸时所产生的冲击波超压。
表71000k8TNT爆炸时冲击波超压
距Ka飞
5
6
7
8
9
1。
14
JK氏
2.94
2.06
1.67
1.27
0.95
0.76
0.50
队33
距离Ra/m
16
18
20
25
30
35
4。
45
R压小4P*
0235
O.t7
0.126
0J26
队057
0.043
0.033
0.027
50
50
60
65
70
75
超隔apo/MPd
0.0035
0.0205
0.016
0,0143
O.Q13
综上所述，计算压力容器爆破时对目标的伤害/破坏作用，可按
下列程序进行。
（1）首先根据容器内所装介质的特性
升级会员 