火灾爆炸事故后果模拟计算.docx
- 文档编号:9948149
- 上传时间:2023-05-22
- 格式:DOCX
- 页数:13
- 大小:281.71KB
火灾爆炸事故后果模拟计算.docx
《火灾爆炸事故后果模拟计算.docx》由会员分享,可在线阅读,更多相关《火灾爆炸事故后果模拟计算.docx(13页珍藏版)》请在冰点文库上搜索。
火灾爆炸事故后果模拟计算
火灾、爆炸事故后果模拟计算
在化工生产中,火灾、爆炸和中毒事故不但影响生产的正常运行,而且对人员有较大的身体危害,导致人员的伤亡。
本文运用地面火灾、蒸气云爆炸和中毒的三种数学模型,对年产2万吨顺酐装置的原料库来进行分析,分析各种事故对人员可能造成的危害,借以帮助企业在生产中采取相应的措施。
事故后果分析是危险源危险性分析的一个主要组成部分,其目的在于定量描述一个可能发生的重大事故对工厂、对厂内人员、厂外居民甚至对环境造成危害的严重程度。
一、苯储罐泄漏池火灾后果分析
苯系易燃液体,在苯贮罐区苯泄漏后遇到点火源就会被点燃而着火燃烧。
由于贮罐区设有防火堤,苯泄漏后积聚在防火堤之内,它被点燃后的燃烧方式为池火。
模拟有关数据参数如下。
苯储罐区有两台800m3、两台500m3的苯储罐,苯储罐单罐直径10.5m,每两台罐为一组,贮罐区防火堤尺寸为33×16m,模拟液池半径为18.3m;
苯储罐单台最大贮存量600t,泄漏量为15%时,足以在防火堤内形成液池;
周围环境温度设为25℃;
(1)燃烧速度
当液池中的可燃液体的沸点高于周围环境温度时,液体表面上单位面积的燃烧速度dm/dt为:
0.001Hc
dm/dt=
CP(Tb-T0)+H
式中dm/dt~单位表面积的燃烧速度,kg/m2.s
Hc~液体燃烧热,J/kg。
苯Hc=41792344J/kg。
CP~液体的定压比热容,J/kg.K。
苯CP=1729J/kg.K。
Tb~液体的沸点,Kb=353.1K。
T0~环境温度,环境温度为25℃,K。
=298K。
H~液体的气化热,J/kg。
苯H=428325J/kg。
(25℃)
计算:
dm/dt=0.001×41792344/﹝1729(353.1-298)+428325﹞
=0.0798kg/m2.s
(2)火焰高度
模拟液池为园池,半径为18.3m,其火焰高度可按下式计算:
dm/dt
h=84r﹝﹞0.61
ρ0(2gr)1/2
式中h~火焰高度,m;
r~液池半径,m;取r=18.3m
ρ0~周围空气密度,kg/m3;取ρ0=1.185kg/m3(25℃)
g~重力加速度,9.8m/s2;
dm/dt~单位表面积的燃烧速度,己知0.0798kg/m2.s
计算:
h=84×18.3×{0.0798/[1.185×(2×9.8×18.3)1/2]}0.61
=49.3m
(3)热辐射通量
当液池燃烧时放出的总热辐射通量为:
Q=(兀r2+2兀rh)dm/dt·η·Hc/﹝72(dm/dt)0.6+1﹞…(公式F5-3)
Q~总热辐射通量,W;
η~效率因子,可取0.13~0.35。
取0.15。
其余符号意义同前。
Q=(3.14×18.32+2×3.14×18.3×49.3)×0.0798×0.15×41792344÷(72×0.07980.6+1)
=200068980.5W
(4)目标入射热辐射强度
假设全部辐射热量由液池中心点的小球面辐射出来,则在距离池中心某一距离(X)处的入射热辐射强度为:
Qtc
I=
4兀X2
式中I~热辐射强度,W/m2;
Q~总热辐射通量,W;
tc~热传导系数,在无相对理想的数据时,可取值为1;
X~目标点到液池中心距离,m;
例如,距离液池中心20m处的入射热辐射强度为:
Qtc200068980.5
I===39822.6W/m2
4兀X24×3.14×202
(5)火灾损失
火灾通过辐射热的方式可影响周围环境。
当火灾产生的热辐射强度足够大时,可使周围的物体燃烧或变形,强烈的热辐射可能烧毁设备甚至造成人员伤亡。
火灾损失估算建立在辐射通量与损失等级的相应关系的基础上。
有资料提供不同入射通量造成伤害或损失的情况,结合本篇所要分析的对象,当目标接受的入射通量为某一定值时,目标与液池中心的距离可用下式计算:
X=(Qtc/4兀I)0.5
式中符号意义与以上相同。
按该式计算:
当I=37500,X=(200068980.5/4×3.14×37500)0.5=20.6m
当I=25000,X=25.2m
当I=12500,X=35.7m
当I=4000,X=63.1m
当I=1600,X=99.8m
不同距离、不同入射通量造成的伤害及损失情况表
目标至池火中心的水平距离(m)
入射通量
(W/m2)
对设备的损害
对人的伤害
20.6
37500
操作设备全部损坏
1%死亡/10s
100%死亡/1min
25.2
25000
在无火焰、长时间辐射下,木材燃烧的最小能量
重大烧伤/10s
100%死亡./1min
35.7
12500
有火焰时,木材燃烧,塑料熔化的最低能量
1度烧伤/10s
1%死亡/1min
63.1
4000
20s以上感觉疼痛,未必起泡
99.8
1600
长期辐射无不舒服感
二、苯储罐泄漏爆炸后果分析
采用蒸气云爆炸(VCE)伤害模型计算苯储罐泄漏爆炸的伤害/破坏半径,在伤害/破坏半径所包围的封闭面积内人员多少、财产价值多少将影响事故严重度大小。
(1)有关参数
苯储罐区有两台800m3、两台600m3的苯储罐,储存物质:
苯,单台最大储存量600t。
苯燃烧热41792kJ/kg。
(2)TNT当量计算通常以TNT当量法来预测蒸气云爆炸的严重程度。
蒸气云的TNT当量WTNT计算式如下:
WTNT=1.8·a·Wf·Qf/QTNT…………………(公式F5-5)
式中:
1.8为地面爆炸系数
WTNT——蒸气云的TNT当量,kg;
a——蒸气云的TNT当量系数,取4%;
Wf——蒸气云爆炸中燃烧的总质量,kg;
Qf——为计算对象的燃烧热,kJ/kg;
QTNT——TNT的爆热,KJ/kg,取平均爆破能量值4520kJ/kg。
设:
泄漏出来的苯体积容量为5%,
Wf=600000×5%=30000kg,
WTNT=1.8×0.04×Wf×Qf/4520
=1.8×0.04×30000×41792/4520
=19971kgTNT
死亡区内的人员如缺少防护,则被认为无例外地蒙受严重伤亡或死亡,其内径为零,外径为R1;重伤区内的人员如缺少防护,则被认为无例外地蒙受严重伤害,极少数人员可能死亡或受轻伤。
其内径为R1,外径为R2;轻伤区内的人员如缺少防护,绝大多数人员将遭受轻微伤害,少数人员将受重伤或无事,其内径为R2,外径为R3;安全区内的人员不会受到伤害,该区的内径为R3o各区的计算式如下:
R1=13.6(WTNT/1000)0.37
ΔPS=0.137Z-3+0.119Z-2+0.269Z-1-0.019
Z=R2(E/Po)1/3或Z=R3(E/Po)1/3
ΔPS=44000/Po
E=WTNT·QTNT
式中:
ΔP——引起人员伤害的冲击波峰值:
重伤为44000Pa、轻伤为17000Pa;
E——爆源总能量,J;
Po——为环境压力,Pa;
QTNT——TNT的爆热,KJ/kg。
(3)死亡半径R1
R1=13.6×(WTNT/1000)0.37………………(公式F5-6)
=13.6×(10000/1000)0.37
=41.2m
(4)重伤半径R2
ΔPS=0.137Z-3+0.119Z-2+0.269Z-1-0.019
Z=R2/(E/Po)1/3=0.01039R2
ΔPS=44000/Po=0.4344
E=WTNT·QTNT=90270720000J
式中:
P0为环境大气压,同上。
101300Pa
E为总爆炸能量,同上。
计算得R2=104.8m
(5)轻伤半径R3
ΔPS=0.137Z-3+0.119Z-2+0.269Z-1-0.019
Z=R3/(E/Po)1/3=0.01039R3
ΔPS=17000/Po=0.1678
式中:
P0为环境大气压,取101300Pa
E为总爆炸能量,
计算得R3=188.3m
(6)财产损失半径R财
R财=KWTNT1/3/﹝1+(3175/WTNT)2﹞1/6…………(公式F5-9)
式中K为二级破坏系数取K=5.6
R财=KWTNT1/3/﹝1+(3175/WTNT)2﹞1/6
=5.6×199711/3/﹝1+(3175/19971)2﹞1/6
=151.3m
以上通过苯罐区的火灾爆炸事故案例,对易燃易爆液体储罐区火灾爆炸事故的评价方法选择和评价过程的介绍可以看出,对于易燃易爆液体储罐区火灾爆炸事故的评价,首先应采用美国道化公司的火灾爆炸危险指数法,评价出火灾爆炸事故的影响范围;其次利用池火灾伤害数学模型分析法,评价出事故对周围设施和人员损伤程度,同时还可以对前述道化学方法的评价结论进行验证;最后建立火灾爆炸事故树,分析该事故发生的概率(必须有基本原因事件发生的记录或经验概率值)和找出导致火灾爆炸事故发生的基本原因事件,并求出各原因事件的主次作用,为安全对策措施的轻重缓急排序提供科学依据。
重大事故后果分析方法:
中毒
有毒物质泄漏后生成有毒蒸气云,在空气中飘移、扩散,直接影响现场人员并可能波及居民区。
大量剧毒物质泄漏可能带来严重的人员伤亡和环境污染。
毒物对人员的危害程度取决于毒物的性质、毒物的浓度和人员与毒物接触时间等因素。
有毒物质泄漏初期,其毒气形成气团密集在泄漏源周围,随后由于环境温度、地形、风力和湍流等影响气团飘移、扩散,扩散范围变大,浓度减小。
在后果分析中,往往不考虑毒物泄漏的初期情况,即工厂范围内的现场情况,主要计算毒气气团在空气中飘移、扩散的范围、浓度、接触毒物的人数等。
1描述毒物泄漏后果的概率函数法
概率函数法是通过人们在一定时间接触,一定浓度所造成影响的概率来描述毒物泄漏后果的一种表示法。
概率与中毒死亡百分率有直接关系,二者可以互相换算,如表1所示。
概率值在0~9之间。
表1概率与死亡百分率的换算
概率值Y与接触毒物浓度及接触时间的关系如下:
Y=A+B[n(cnt)
(1)
式中A,B,n——取决于毒物性质的常数,表6—15列出了——些常见有毒物质的有关参数;
c——接触毒物的浓度,ppm;
t——接触毒物的时间,min。
表2一些毒性物质的常数
物质名称
A
B
n
参考资料
氯
-5.3
0.5
2.75
DCMR1984
氨
-9.82
0.71
2.0
DCMR1984
丙稀醛
-9.93
2.05
1.0
USCG1977
四氯化碳
0.54
1.01
0.5
USCG1977
氯化氢
-21.76
2.65
1.0
USCG1977
甲基溴
-19.92
5.16
1.0
USCG1977
光气
-19.27
3.69
1.0
USCG1977
氟氢酸(单体)
-26.4
3.35
1.0
USCG1977
使用概率函数表达式时,必须计算评价点的毒性负荷(cnt),因为在一个已知点,其毒性浓度随着气团的通过和稀释而不断变化,瞬时泄漏就是这种情况。
确定毒物泄漏范围内某点的毒性负荷,可把气团经过该点的时间划分为若干区段,计算每个区段内该点的毒物浓度,得到各时间区段的毒性负荷,然后再求出总毒性负荷:
总毒性负荷=∑时间区段内毒性负荷
一般说来,接触毒物的时间不会超过30min,在这段时间里人员可以逃离现场或采取保护措施。
当毒物连续泄漏时,某点的毒物浓度在整个云团扩散期间没有变化。
当设定某死亡百分率时,由表1查出相应的概率Y值,根据式1有:
(2)
计算出c值,按扩散公式可以算出中毒范围。
如果毒物泄漏是瞬时的,则有毒气团在某点通过时该点处毒物浓度是变化的。
这种情况下,考虑浓度的变化情况,计算气团通过该点的毒性负荷,算出该点的概率值Y,然后查表1就可得出相应的死亡百分率。
2有毒液化气体容器破裂时的毒害区估算
液化介质在容器破裂时会发生蒸气爆炸。
当液化介质为有毒物质,如液氯、液氨、二氧化硫、氢氰酸等,爆炸后如果不燃烧,会造成大面积的毒害区域。
设有毒液化气体质量为W(kg),容器破裂前器内介质温度为t(℃),液体介质比热为c(kJ/(kg·℃)),当容器破裂时,器内压力降至latm(0.1MPa),处于过热状态的液体温度迅速降至标准沸点t0(℃),此时全部液体所放出的热量为:
Q=Wc(t—t0) (3)
设这些热量全部用于器内液体的蒸发,如它的气化热为g(kJ/kg),则其蒸发量为:
(4)
如介质的相对分子质量为M,则在沸点下蒸发蒸气的体积Vg(m3)为:
(5)
为便于计算,现将压力容器最常用的液氨、液氯、氢氰酸等的有关物理化学性能列于表3中。
关于一些有毒气体的危险浓度如表4所示。
表3一些有毒物质的有关物化性能
表4有毒气体的危险浓度
如已知某种有毒物质的危险浓度,则可求出其危险浓度下的有毒空气体积。
如二氧化硫在空气中的浓度达到0.05%时,人吸入5~10min即致死,则Vg(m3)的二氧化硫气可以令人致死的有毒空气体积为:
V=Vg×100/0.05=2000Vg
假设这些有毒空气以半球形向地面扩散,则可求出该有毒气体扩散半径为:
(6)
式中R——有毒气体的半径,m;
Vg——有毒介质的蒸气体积,m3;
c——有毒介质在空气中危险浓度值,%。
3有毒物质喷射泄漏时的毒害区估算
关于有毒物质喷射泄漏的毒害区估算可参考式6—17和式6—18进行。
式中b1、b2——分布函数,其表达式如下:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 火灾 爆炸 事故 后果 模拟 计算