爆炸评价模型及伤害半径计算文档格式.docx
《爆炸评价模型及伤害半径计算文档格式.docx》由会员分享,可在线阅读,更多相关《爆炸评价模型及伤害半径计算文档格式.docx(19页珍藏版)》请在冰点文库上搜索。
SP——环境压力(101300Pa);
0E——爆炸总能量(J),E=W×
Q。
TNTTNT将以上数据代入方程式,解得:
△P=0.43442Z=1.07
21/3)×
1000/101300=1.07R×
(27739×
45002=1.07×
107=115(m)
轻伤半径R,由下列方程式求解:
3-3-2-1-0.019+0.269Z=0.137Z△P+0.119Z33331/3)/(E/P=RZ0332/21
△P=△P/P0S3式中:
△P——引起人员轻伤冲击波峰值,取17000Pa。
S将以上数据代入方程式,解得:
△P=0.168,Z=1.9533轻伤半径R=209(m)32、沸腾液体扩展蒸气爆炸(BLEVE)模型分析计算
(1)沸腾液体扩展蒸气爆炸(BLEVE)模型
液态存贮的易燃液化气体突然瞬间泄漏时,立即遇到火源就会发生剧烈的燃烧,产生巨大的火球,形成强烈的热辐射,此种现象称为沸腾液体扩展蒸气爆炸,简称BLEVE。
沸腾液体扩展蒸气爆炸的主要危险是强烈的热辐射,近场以外的压力效应不重要。
其火球的特征可用国际劳工组织(ILO)建议的蒸气爆炸模型来估算。
火球半径的计算公式为:
1/3R=2.9W式中R——火球半径,m;
W——火球中消耗的可燃物质量,kg。
对单罐储存,W取罐容量的50%;
双罐储存;
W取罐容量的70%;
多罐储存,取W为罐容量的90%。
(2)液氨储罐沸腾液体扩展蒸气爆炸(BLEVE)模型分析计算
由于生产装置液氨贮罐区的液氨罐为多罐贮存,(共六只贮罐,其中三只333罐,则0.6,取50M100M,三只100M)最大库存量为250T。
氨比重约由W=100×
0.6×
1000×
90%=54000(kg)
代入式中,得到:
1/3=109(m)R=2.9(54000)火球半径火球持续时间按下式计算:
1/3t=0.45W式中:
火球持续时间,单位为S.
将数据代入式中,得到:
3/21
1/3=17(s)
(54000)t=0.45×
目标接收到热辐射通量的计算,按下式计算:
(r)2223/2=qR+rr(1-0.058Inr)/(Rq)0式中:
r——目标到火球中心的水平距离,m;
22,球形罐取。
对柱形罐取270kW/mq——火球表面的辐射通量,W/m02。
200kW/m
R——火球半径,m。
R=109m。
有了热辐射q(r),即可求不同伤害、破坏时的热通量及其半径。
下面求不同伤害时的热通量:
死亡可根据下式计算:
4/3)P=-36.38+2.56In(tq1r式中:
P=5
rt——火球持续时间,取t=17s。
2。
解得q=21985W/m1重伤可根据下式计算:
4/3)=-43.143+3.0188In(tqP2r2。
q=18693W/m解得2轻伤可根据下式计算:
4/3)P=-39.83+3.0188In(tq3r2。
=8207W/m解得q3通过q、q、q可以求得对应的死亡半径R、重伤半径R及轻伤半径R。
(由313122于此方程式难以手算解出,故省略)。
(3)小结
通过计算,如果贮存区液氨储罐发生扩展蒸气爆炸,火球半径为109m。
将可能造成其他贮罐的连锁火灾和爆炸,造成灾难性的破坏。
3、液氨泄漏中毒事故的模拟计算
液氨贮存区最大贮存量为250T,假设有1T泄漏量,对蒸发成蒸气扩散造成的危害进行模拟计算。
4/21
)液态气体蒸气体积膨胀计算(1升体积。
根据液态气摩尔气体占有22.41013Mpa),1℃,在标准状态下(0体的相对密度,由下式可计算出它们气化后膨胀的体积:
DV?
00422.?
1000?
V?
M
V——膨胀后的体积(升)V——液态气体的体积(升)0)D——液态气体的相对密度(水=10——液态气体的的分子量M得到=0.597,M=17.03将液氨有关数据代入上式,由D0D?
V00422?
.V?
?
1000M5970.V?
0V785.224?
003.17
倍。
即液态氨若发生泄漏迅速气化,其膨胀体积为原液态体积的7852()液态气体扩散半径模拟计算且液态气液态气体泄漏后在高温下迅速气化并扩散,在一定泄漏量范围内,其短轴与沿地面能扩散到相当远的地方,可模拟为半椭圆形,体比重大于空气,长轴之比将随着扩散半径的增大而减少,可由下式计算:
R23?
3
式中:
V——液态气体膨胀后体积;
5/21
ρ——液态气体在空气中的浓度;
K=h/R。
κ——椭圆形短轴与长轴之比,即)查得:
液氨在工(GBZ2-2002根据我国《工作场所有害因素职业接触限值》33其20mg/m,;
短时间接触容许浓度30mg/m作场所空气中时间加权平均容许浓度-6-6。
ρ=39.5×
在空气中体积浓度换算为:
ρ=26.3×
1010和,其可能发生中毒事故的浓度区域半径计算如1000kg假设泄漏液氨的量为下:
3=1/0.597=1.68m取液氨体积V0-6K=0.10ρ=26.3×
10计算:
3785?
1.68?
23?
6?
26.3?
10?
0.10?
3?
621m
从计算结果可知:
当泄漏1000kg液氨气化成蒸气时可能发生中毒浓度的区域半径为621m,即0.621公里,因此,其扩散时的可能发生中毒浓度的区域面积:
22=1.21(平方公里)=3.14×
S=π×
R0.6214、水煤气泄漏事故的模拟计算
根据我国《工作场所有害因素职业接触限值》(GBZ2-2002)查得:
水煤气(即33。
30mg/m20mg/m一氧化碳)时间加权平均容许浓度;
短时间接触容许浓度经换算,分别为ρ=16×
10-6和ρ=24×
10-6。
水煤气贮罐总容积5000M3,设若泄漏量为100M3,取ρ=16×
10-6,K=0.1则计算如下:
6/21
3100?
16?
0.103=310(M)
22310=301754M=3.14S=π×
R×
5、天然气泄漏形成喷射火模型分析
该项目设计有容积为2m3、工作压力为25MPa的天然气储气瓶组。
当设备损坏、法兰垫片撕裂或接管破裂等条件下,天然气就会在破裂处形成射流,在高速气流摩擦形成的静电火花或其他点火源存在的条件下,就会在裂口处引燃形成喷射火。
火灾通过热辐射方式影响周围环境,当热辐射强度足够大时,可使周围物体燃烧变形,强烈的热辐射可能烧毁设备并造成人员伤亡。
喷射火模拟事故模拟是通过定量的计算,估算出热辐射的不同入射通量所造成的损失程度。
(1)气体泄漏量的计算
假定天然气储气瓶组或管道泄漏,裂口为直径10mm的圆口:
P0=0.1013×
106Pa;
P=25×
κ=1.314(天然气的绝热指数)
则:
P0.10130?
0.004P25?
1.31422?
?
1?
1.314?
0.54?
1314?
.11?
P2?
?
0?
1P?
故气体流动属于音速流动。
其泄漏量为:
7/21
2M?
)A(Q?
C?
d0?
1RT?
式中:
=;
ρT=293K(20℃)的分子量,=1.00;
M=CNG16;
R=8.314J/mol?
K;
CdD为储气瓶裂口直径。
193.5kg/m3,A为裂口面积,A=πD2/4=0.0000785(m2)
计算得:
10-4kg/s8.2×
Q0=)喷射火热辐射通量(2这里所用的喷射火辐射热计算方法是一种包括气流效应在内的喷射火扩散模式每个点热源的把整个喷射火看成是沿喷射中心线上的几个热源点组成,的扩展。
热辐射通量相等。
点热源的热辐射通量按下式计算?
Hq?
QC0式中:
q―点热源热辐射量,W
η―效率因子,可取0.35
Q0―泄漏速度,kg/s
Hc―燃烧热,J/kg
天然气高热值:
Hc=55800kJ/kg,则:
q=0.35×
8.2×
10-4×
55800=16kW
射流轴线上某点热源I到距离该处一点的热源辐射强度为:
qRI?
i2?
x4
Ii―点热源i到目标点x处的热辐射强度,W/m2;
q―点热源的辐射通量,W;
R―辐射率,可取0.2;
x―点热源到目标点的距离,m
某一点处的入射流强度等于喷射火的全部点热源对目标的热辐射强度的总和:
N?
I?
Ii1I?
8/21
。
=5N―计算时选取的点热源系数,一般取n查热辐射的某点总入射热辐射通量的大小,根据喷射火全部点热源在距火焰xm根据可以计算出有代表意义的入射热辐射通量造的不同入射能量所造成的损失。
成的危害范围。
即:
qR?
x?
I4
R―辐射率,取0.2;
表5.6-1不同入射通量所造成的危害范围
目标伤害距入射通I对设备的损对人的损kW/m1死/10s37.5
0.082
操作设备全部损/1min100死亡在无火焰、长时间辐射/10s重大损伤0.125下,木材燃烧的最小能/1min100%死亡量
有火焰时,木材燃烧,1度烧伤/10s12.5
0.14塑料融化的最低能量死亡/1min
1%以上感觉疼痛,未必20s0.254.0
起泡,无不舒服长时间辐射0.401.6感从上表可以看出,如不采取措施,在计算条件下的喷射火如果发生,距射流轴线热源0.082m处的人员在热辐射下10s内1%死亡,1min内100%死亡,另此范围内的操作设备全部损坏。
距射流轴线热源0.25m,热辐射不会造成人员伤亡;
距射流轴线热源0.40m以外为安全距离。
为了安全生产,企业应积极采取应对措施,防止喷射事故发生。
6、天然气泄漏形成蒸气云爆炸模型分析
9/21
假设1个2m3、操作压力25MPa的天然气储气瓶组全部发生泄漏,则参与爆炸的天然气体积为50m3,并与空气混合后形成蒸气云,蒸汽云爆炸冲击波的破坏半径为:
R=CS(NE)1/3
E-爆炸能量,kJ,E=V·
HC
V-参与反应的可燃气体的体积,m3;
HC-可燃气体的燃烧热值;
39.86MJ/Nm3;
N-效率因子,一般取10%;
CS-经验常数,取决于损害等级。
其具体取值情况见下表:
表5.6-2损害等级表
Cs/m?
J-损害等级设备损害人员伤害(1/3)
1%死亡于肺部伤害;
>
50%耳10.03重创建筑物的加工设备膜破裂;
50%被碎片击伤
损害建筑物外表可修复性20.061%耳膜破裂;
1%被碎片击伤破坏
30.15玻璃破碎被碎玻璃击伤
0.4
10%4
玻璃破碎蒸气云爆炸冲击波的破坏半径与损害等级关系,计算结果见下表:
表5.6-3天然气蒸气云爆炸冲击波损害范围表
损害(NE)损害半径CsNE设备损坏人员伤害等级1/3R
m?
J-单位JJ1/3m(1/3)
1%死亡人员肺部重创建筑物和伤害10.035.99
加工设备>
50%耳膜破裂
50%被碎片击伤7.972×
1损坏建筑物外199.771%耳膜破裂0620.0611.98表,可修复性破1%被碎片击伤坏
30.1529.97玻璃破碎被碎玻璃击伤
4
0.479.910%玻璃破碎
爆炸冲击波伤害半径的模拟分析表明,1只2m3、工作压力为25MPa的天然气储气瓶组全部发生泄漏后,天然气与空气形成爆炸性混合气,发生爆炸后在半径5.99m范围内,1%人员死亡及肺部伤害、>
50%耳膜破裂、>
50%被碎片击伤。
7、压缩气体与水蒸气容器爆破能量
当压力容器中介质为压缩气体,即以气态形式存在而发生物理爆10/21
炸时,其释放的爆破能量为:
(1)
式中Eg——气体的爆破能量,kJ;
p——容器内气体的绝对压力,MPa;
3;
mV——容器的容积,
κ——气体的绝热指数,即气体的定压比热与定容比热之比。
常用气体的绝热指数数值如表1所示。
表1常用气体的绝热指数
从表1可看出,空气、氮、氧、氢及一氧化氮、一氧化碳等气体的绝热指数均为1.4或近似1.4,如用κ=1.4代入式
(1)中,得到气体的爆破能量为:
(2)
Eg=CgV(3)
3。
/Cg——常用压缩气体爆破能量系数,kJm式中
压缩气体爆破能量系数Cg是压力p的函数,各种常用压力下的11/21
气体爆破能量系数如表2所示。
表2常用压力下的气体容器爆破能量系数(κ=1.4时)
如将κ=1.135代入式
(1),可得干饱和蒸汽容器爆破能量为:
(4)
用式4计算有较大的误差,因为没有考虑蒸汽干度的变化和其他一些影响,但可以不用查明蒸汽热力性质而直接计算,对危险性评价可提供参考。
对于常用压力下的干饱和蒸汽容器的爆破能量可按下式计算:
Es=CsV(5)
式中Es——水蒸气的爆破能量,kJ;
mV——水蒸气的体积,
/mkJCs——干饱和水蒸气爆破能量系数,各种常用压力下的干饱和水蒸气容器爆破能量系数如表3所示。
表3常用压力下干饱和水蒸气容器爆破能量系数
12/21
1.2介质全部为液体时的爆破能量
通常用液体加压时所做的功作为常温液体压力容器爆炸时释放的能量,计算公式如下:
(6)
式中E——常温液体压力容器爆炸时释放的能量,kJ;
Lp——液体的压力(绝),Pa;
——容器的体积,mV1—。
Pat——液体在压力卢和温度T下的压缩系数,β
1.3液化气体与高温饱和水的爆破能量
液化气体和高温饱和水一般在容器内以气液两态存在,当容器破裂发生爆炸时,除了气体的急剧膨胀做功外,还有过热液体激烈的蒸发过程。
在大多数情况下,这类容器内的饱和液体占有容器介质重量的绝大部分,它的爆破能量比饱和气体大得多,一般计算时不考虑气体膨胀做的功。
过热状态下液体在容器破裂时释放出爆破能量可按下式计算:
—H—S)T]WE=[(H)—(S(7)
11212式中E——过热状态液体的爆破能量,kJ;
H——爆炸前液化液体的焓,kJ/kg;
1H——在大气压力下饱和液体的焓,kJ/kg;
2S1——爆炸前饱和液体的,熵,kJ/(kg·
℃);
S——在大气压力下饱和液体的熵,kJ/(kg·
2T——介质在大气压力下的沸点,℃;
1W——饱和液体时质量,kg。
13/21
饱和水容器的爆破能量按下式计算:
Ew=CwV(8)
式中Ew——饱和水容器的爆破能量,kJ;
——容器内饱和水所占的容积,mV3,其值如表4所示。
m——饱和水爆破能量系数,kJ/Cw
表4常用压力下饱和水爆破能量系数
8、冲击波超压的伤害-破坏作用
压力容器爆破时,爆破能量在向外释放时以冲击波能量、碎片能量和容器残余变形能量3种形式表现出来。
根据介绍,后二者所消耗的能量只占总爆破能量的3%~15%,也就是说大部分能量是产生空气冲击波。
冲击波是由压缩波迭加形成的,是波阵面以突进形式在介质中传播的压缩波。
容器破裂时,容器内的高压气体大量冲出,使它周围的空气受到冲击而发生扰动,使其状态(压力、密度、温度等)发生突跃变化,其传播速度大于扰动介质的声速,这种扰动在空气中传播就成为冲击波。
在离爆破中心一定距离的地方,空气压力会随时间迅速发生而悬殊的变化。
开始时,压力突然升高,产生一个很大的正压力,接着又迅速衰减,在很短时间内正压降至负压。
如此反复循环数次,压力渐次衰减下去。
开始时产生的最大正压力即是冲击波波阵面上的超压△p。
多数情况下,冲击波的伤害一破坏作用是由超压引起的。
14/21
超压△p可以达到数个甚至数十个大气压。
冲击波伤害一破坏作用准则有:
超压准则、冲量准则、超压一冲量准则等。
为了便于操作,下面仅介绍超压准则。
超压准则认为,只要冲击波超压达到一定值时,便会对目标造成一定的伤害或破坏。
超压波对人体的伤害和对建筑物的破坏作用如表5和表6所示。
表5冲击波超压对人体的伤害作用
表6冲击波超压对建筑物的破坏作用
2.2冲击波的超压
冲击波波阵面上的超压与产生冲击波的能量有关,同时也与距离爆炸中心的远近有关。
冲击波的超压与爆炸中心距离的关系:
n—∝△pR(9)
式中△p——冲击波波阵面上的超压,MPa;
R——距爆炸中心的距离,m;
n——衰减系数。
衰减系数在空气中随着超压的大小而变化,在爆炸中心附近内为15/21
2.5~3;
当超压在数个大气压以内时,n=2;
小于1atm(0.1MPa)时,n=1.5。
实验数据表明,不同数量的同类炸药发生爆炸时,如果距离爆炸中心的距离R之比与炸药量q三次方根之比相等,则所产生的冲击波超压相同,用公式表示如下:
如
(10)
式中R——目标与爆炸中心距离,mRm;
——目标与基准爆炸中心的相当距离,0;
TNT,kgq——基准爆炸能量,0q——爆炸时产生冲击波所消耗的能量,TNT,kg;
△p——目标处的超压,MPa;
△p——基准目标处的超压,MPa;
0α——炸药爆炸试验的模拟比。
式10也可写成为:
△p(R)=△p(R/α)(11)
0利用式11就可以根据某些已知药量的试验所测得的超压来确定在各种相应距离下任意药量爆炸时的超压。
表7是1000kgTNT炸药在空气中爆炸时所产生的冲击波超压。
表71000k8TNT爆炸时冲击波超压
16/21
综上所述,计算压力容器爆破时对目标的伤害/破坏作用,可按下列程序进行。
(1)首先根据容器内所装介质的特性,分别选用式2~式8计算出其爆破能量正。
(2)将爆破能量q换算成TNT当量q,因为1kgTNT爆炸所放出0的爆破能量为4230kJ/kg~4836kJ/kg,一般取平均爆破为4500kJ/kg,故其关系为:
q=E/q=E/4500(12)INT(3)按式10求出爆炸的模拟比α,即:
1/31/31/3=(q/q)=(q/1000)=0.1q(13)α0(4)求出在1000kgTNT爆炸试验中的相当距离R,即R=R/α。
00(5)根据R值在表7中找出距离为R0处的超压△p(中间值用插00入法),此即所求距离为R处的超压。
(6)根据超压△p值,从表5和表6中找出对人员和建筑物的伤害一破坏作用。
2.3蒸气云爆炸的冲击波伤害一破坏半径
爆炸性气体以液态储存,如果瞬态泄漏后遇到延迟点火或气态储存时泄漏到空气中遇到火源,则可能发生蒸气云爆炸。
导致蒸气云形成的力来自容器内含有的能量或可燃物含有的内能,或两者兼而有17/21
之。
“能”主要形式是压缩能、化学能或热能。
一般说来,只有压缩能和热能才能单独导致形成蒸气云。
根据荷兰应用科研院(TNO(1979))建议,可按下式预测蒸气云爆炸的冲击波损害半径:
1/3R=Cs(NE)(14)
式中R——损害半径,m;
E——爆炸能量,kJ,可按下式取:
E=VHc(15)
mV——参与反应的可燃气体的体积
Hc——可燃气体的高燃烧热值,取值情况如表8所示;
N——效率因子,其值与燃料浓度持续展开所造成损耗的比例和燃料燃烧所得机械能的数量有关,一般取N=10%;
Cs——经验常数,取决于损害等级,其取值情况如表9所示。
3)
m/8表某些气体的高燃烧热值(kJ
表9损害等级表
18/21
冲击波计算、9氧气钢瓶发生物理爆炸,具体分析如下:
TNT当量计算
当氧气钢瓶发生爆炸时,气体膨胀所释放的能量(即爆破能量)不仅与气体压力和储罐的容积有关而且与介质在容器内的物性相态相关。
氧气为非热力气体,无焓值、熵值;
承压状态下称压缩气体,承压钢瓶破裂时属物理性爆炸;
其能量计算,与瓶内压力、瓶体容积、气体绝热指数有关。
本项目中运用压缩气体爆破能量计算模型计算,其释放的爆破能量为:
k-1/k3[1-(0.1013/p)]×
10Eg=2.5PV/(k-1)式中,Eg-气体的爆破能量,kJ;
P-容器内气体的绝对压力,MPa;
m容器的容积,V-k-气体的绝热指数,即气体的定压比热与定容比热之比,此处取1.4;
0.28573令:
Cg=2.5P[1-(0.1013/P)]×
1019/21
则:
Eg=Cg·
V