2.2事故后果模拟分析法火灾、爆炸、中毒是常见的重大事故,经常造成严重的人员伤亡和巨大的财产损失,影响社会安定这里重点介绍有关火灾、爆炸和中毒事故(热辐射、爆炸波、中毒)后果分析,在分析过程中运用了数学模型通常一个复杂的问题或现象用数学模型来描述,往往是在一个系列的假设的前提下按理想的情况建立的,有递增模型经过小型试验的验证,有的则可能与实际情况有较大出入,但对辨识危险性来说是可参考的2.2.1 泄漏由于设备损坏或操作失误引起泄漏,大量易燃、易爆、有毒有害物质的释放,将会导致火灾、爆炸、中毒等重大事故发生,因此,后果分析由泄漏分析开始2.2.1.1 泄漏情况分析2.1.1.1.1泄漏的主要设备根据各种设备泄漏情况分析,可将工厂(特别是化工厂)中易发生泄漏的设备归纳为以下10类:管道、挠性连接器、过滤器、阀门、压力容器或反应器、泵、压缩机、储罐、加压或冷冻气体容器,火炬燃烧装置或放散管等 ⑴管道它包括管道、法兰和接头,其典型情况和裂口尺寸分别取管径的20%~100%、20%和20%~100%⑵挠性连接器它包括软管、波纹管和铰接器,其典型泄漏情况和裂口尺寸为:①连接器本体破裂泄漏,裂口尺寸取管径的20%~100%;②接头处的泄漏,裂口尺寸取管径的20%;③连接装置损坏泄漏,裂口尺寸取管径的100%。
⑶过滤器它由过滤器本体、管道、滤网等组成,其典型泄漏情况和裂口尺寸分别取管径的20%~100%和20%其典型泄漏情况和裂口尺寸为:①阀壳体泄漏,裂口尺寸取管径的20%~100%;②阀盖泄漏,裂口尺寸取管径的20%;③阀杆损坏泄漏,裂口尺寸取管径的20%⑸压力容器或反应器包括化工生产中常用的分离器、气体洗涤器、反应釜、热交换器、各种罐和容器等其常见的此类泄漏情况和裂口尺寸为:①容器破裂而泄漏,裂口尺寸取容器本身尺寸;②容器本体泄漏,裂口尺寸取与其连接的粗管道管径的100%;③孔盖泄漏,裂口尺寸取管径的20%;④喷嘴断裂而泄漏,裂口尺寸取管径的100%;⑤仪表管路破裂泄漏,裂口尺寸取管径的20%~100%;⑥容器内部爆炸,全部破裂其典型泄漏情况和裂口尺寸为:①泵体损坏泄漏,裂口尺寸取与其连接管径的20%~100%;②密封压盖处泄漏,裂口尺寸取管径的20%;⑺压缩机包括离心式、轴流式和往复式压缩机,其典型泄漏情况和裂口尺寸为:①压缩机机壳损坏而泄漏,裂口尺寸取与其连接管道管径的20%~100%;②压缩机密封套泄漏,裂口尺寸取管径的20%露天储存危险物质的容器或压力容器,也包括与其连接的管道和辅助设备,其典型泄漏情况和裂口尺寸为:①罐体损坏而泄漏,裂口尺寸为本体尺寸;②接头泄漏,裂口尺寸为与其连接管道管径的20%~100%;③辅助设备泄漏,酌情确定裂口尺寸。
⑼加压或冷冻气体容器包括露天或埋地放置的储存器、压力容器或运输槽车等,其典型泄漏情况和裂口尺寸为:①露天容器内部气体爆炸使容器完全破裂,裂口尺寸取本体尺寸;②容器破裂而泄漏,裂口尺寸取本体尺寸;③焊接点(接管)断裂泄漏,取管径的20%~100%⑽火炬燃烧器或放散管它们包括燃烧装置、放散管、多通接头、气体洗涤器和分离罐等,泄漏主要发生在筒体和多通接头部位,裂口尺寸取管径的20%~100%表2-1 典型泄漏情况和裂口尺寸取值表序号设备典型泄漏情况裂口尺寸取值1管道包括管道、法兰和接头管道泄漏管径的20%~100%法兰泄漏管径的20%接头泄漏管径的20%~100%2绕性连接器包括软管、波纹管和铰接器连接器本体破裂泄漏管径的20%~100%接头处的泄漏管径的20%连接装置损坏泄漏管径的100%3过滤器由过滤器本体、管道、滤网等组成本体泄漏管径的20%~100%管道泄漏管径的20%~100%滤网泄漏管径的20%4阀阀壳体泄漏管径的20%~100%阀盖泄漏管径的20%阀杆损坏泄漏管径的20%5压力容器或反应器包括化工生产中常用的分离器、气体洗涤器、反应釜、热交换器、各种罐和容器等容器破裂而泄漏容器本身尺寸容器本体泄漏与其连接的粗管道管径的100%孔盖泄漏管径的20%喷嘴断裂而泄漏管径的100%仪表管路破裂泄漏管径的20%~100%容器内部爆炸全部破裂6泵泵体损坏泄漏与其连接管径的20%~100%密封压盖处泄漏管径的20%7压缩机包括离心式、轴流式和往复式压缩机压缩机机壳损坏泄漏与其连接管径的20%~100%压缩机密封套泄漏管径的20%8储罐露天储存危险物质的容器或压力容器罐体损坏而泄漏裂口尺寸为本体尺寸接头泄漏与其连接管道管径的20%~100%辅助设备泄漏酌情确定裂口尺寸9加压或冷冻气体容器包括露天或埋地放置的储存器、压力容器或运输槽车等露天容器内部气体爆炸使容器完全破裂本体尺寸容器破裂而泄漏本体尺寸焊接点(接管)断裂泄漏管径的20%~100%10火炬燃烧器或放散管包括燃烧装置、放散管、多通接头、气体洗涤器和分离罐等筒体和多通接头部位泄漏管径的20%~100%2.1.1.1.2造成泄漏的原因从人-机系统来考虑造成各种泄漏事故的原因主要有四类:⑴设计失误①基础设计错误,如地基下沉,造成容器底部产生裂缝,或设备变形、错位等;②选材不当,如强度不够、耐腐蚀性差、规格不符等;③布置不合理,如压缩机和输出管道没有弹性连接,因振动而使管道破裂;④选用机械不合适,如转速过高,耐温、耐压性能差等;⑤选用计测仪器不合适;⑥储罐、储槽未加液位计,反应器(炉)未加溢流管或放散管等。
⑵设备原因①加工不符合要求,或未经检验擅自采用代用材料;②加工质量差,特别是不具有操作证的焊工焊接质量差;③施工和安装精度不高,如泵和电动机不同轴,机械设备不平衡,管道连接不严密等;④选用的标准定型产品质量不合格;⑤对安装的设备未按《机械设备安装工程及验收规范》进行验收;⑥设备长期使用后未按规定检修期进行检修,或检修质量差造成泄漏;⑦计测仪表未定期校验,造成计量不准;⑧阀门损坏或开关泄漏,又未及时更换;⑨设备附件质量差,或长期使用后材料变质、腐蚀或破裂等⑶管理原因①没有制定完善的安全操作规程;②对安全漠不关心,已发现问题不及时解决;③没有严格执行监督检查制度;④指挥错误,甚至违章指挥;⑤让未经培训的工人上岗,知识不足,不能判断错误;⑥检修制度不严,没有及时检修已出现故障的设备,使设备带病运转⑷人为失误①误操作,违反操作规程;②判断失误,如记错阀门位置或开错阀门;③擅自离岗;④思想不集中;⑤发现异常现象不知如何处理;2.1.1.1.3泄漏后果泄漏一旦出现,其后果不单与物质的数量、易燃性、毒性有关,而且与泄漏物质的相态、压力、温度等状态有关这些状态可有多种不同的结合,在后果分析中,常见的可能结合有4种: ①常压液体;②加压液化气体;③低温液化气体;④加压气体。
泄漏物质的物性不同,其泄漏后果也不同⑴可燃气体泄漏可燃气体泄漏后与空气混合达到爆炸极限时,遇到引火源就会发生燃烧或爆炸泄漏后起火的时间不同,泄漏后果也不相同①立即起火可燃气体从容器中往外泄出时即被点燃,发生扩散燃烧,产生喷射性火焰或形成火球,它能迅速地危及泄漏现场,但很少会影响到厂区的外部②泄后起火可燃气体泄出后与空气混合形成可燃蒸气云团,并随风飘移,遇火源发生爆炸或爆轰,能引起较大范围的破坏⑵有毒气体泄漏有毒气体泄漏形成云团在空气中扩散,有毒气体的浓密云团将笼罩很大的空间,影响范围大⑶液体泄漏一般情况下,泄漏的液体在空气中蒸发而生成气体,泄漏后果与液体的性质和储存条件(温度、压力)有关①常温常压下液体泄漏这种液体泄漏后聚集在防液堤内或地势低洼处形成液池,液体由于持表面风的对流而缓慢蒸发,若遇引火源就会发生池火灾②加压液化气体泄漏一些液体泄漏时将瞬间蒸发,剩下的液体将形成一个液池,吸收周围的热量继续蒸发液体瞬时蒸发的比例决定于物质的性质及环境温度有些泄漏物可能在泄漏过程中全部蒸发③低温液体泄漏这种液体泄漏时将形成液池,吸收周围热量蒸发,蒸发量低于加压液化气体的泄漏量,高于常温常压下液体的泄漏量。
无论是气体泄漏还是液体泄漏,泄漏量的多少都是决定后果严重程度的主要因素,而泄漏量又与泄漏时间长短有关2.2.1.2 泄漏量的计算 当发生泄漏的设备的裂口是规则的,而且裂口尺寸及泄漏物质的有关热力学、物理化学性质及参数已知时,可根据流体力学中的有关方程式计算泄漏量当裂口不规则时,可采取等效尺寸代替;当遇到泄漏过程中压力变化等情况时,往往采用经验公式计算1)液体泄漏计算 液体泄漏速度可用流体力学的伯努利方程计算,其泄漏速度为:Q0=CdAρ(2[p+P0]/ρ+2gh)1/2 (1)式中 Q0——液体泄漏速度,kg/s;Cd——液体泄漏系数,按表2-1选取;A ——裂口面积,m2;ρ——泄漏液体密度,㎏∕m3;p ——容器内介质压力,Pa;p0——环境压力,Pa;g ——重力加速度,9.8m∕s2;h ——裂口之上液位高度,m;表2-2 液体泄漏系数Cd雷诺数(Re)裂口形状圆形(多边形)三角形长方形>1000.650.600.55≤1000.500.450.40对于常压下的液体泄漏速度,取决于裂口之上液位的高低;对于非常压下的液体泄漏速度,主要取决于容器内介质与环境压力之差和液位高低。
当容器内液体是过热液体,即液体的沸点低于周围环境温度,液体流过裂口时由于压力减小而突然蒸发蒸发所需热量取自于液体本身,而容器内剩下液体的温度将降至常压沸点在这种情况下,泄漏时直接蒸发的液体所占百分比F可按下式计算:F=Cp(T-T0)/H (2)式中 Cp——液体的定压比热,J∕kgK;T ——泄漏前液体的温度,K;T0——液体在常压下的沸点,K;H ——液体的气化热,J∕kg;按式(2)计算的结果,几乎总是在0~1之间事实上,泄漏时直接蒸发的液体将以细小烟雾的形式形成云团,与空气相混合而吸收蒸发如果空气传给液体烟雾的热量不足以使其蒸发,有一些液体烟雾将凝结成液滴降落到地面,形成液池根据经验,当F>0.2时,一般不会形成液池;当F<0.2时,F与带走液体之比有线性关系,即当F=0时没有液体带走(蒸发),当F=0.1时有50%的液体被带走2)气体泄漏量计算气体从裂口泄漏的速度与其流动状态有关因此,计算泄漏量时首先要判断泄漏时气体流动属于音速还是亚音速流动,前者称为临界流,后者称为次临界流当式(3)成立时,气体流动属音速流动:p0/p≤[2/(k+1)]k/(k-1) (3)当式(4)成立时,气体流动属亚音速流动:p0/p≥[2/(k+1)]k/(k-1) (4)式中p0 、p——符号意义同前;k——气体的绝热指数,即定压比热Cp与定容比热Cv之比。
气体呈音速流动时,其泄漏量为:Q0=CdAρ[(Mk/RT)(2/(k+1))k+1/(k-1)]1/2 (5)气体呈亚音速流动时,其泄漏量为:Q0=YCdAρ[(Mk/RT)(2/(k+1))k+1/(k-1)]1/2 (6)上两式中 Cd——气体泄漏系数,当裂口形状为圆形时取1.00,三角形时取0.95,长方形时取0.90;Y——气体膨胀因子,它由下式计算:Y={(1/(k-1))((k+1)/2)k+1/k-1(p0/p)2/k[1-(p0/p)(k-1)/k]}1/2(7)M——分子量;ρ——气体密度,kg/m3;R——气体常数,J/mol•K;T——气体温度,K当容器内物质随泄漏而减少或压力降低而影响泄漏速度时,泄漏速度的计算比复杂如果流速小或时间短,在后果计算中可采用最初排放速度,否则应计算等效泄漏速度3)两相流动泄漏量计算在过热液体发生泄漏时,有时会出现气、液两相流动均匀两相流动的泄漏速度可按下式计算:Q0=CdA[2ρ(p-pc)]1/2 (8)式中Q0——两相流动混合物泄漏速度,kg/s;Cd——两相流动混合物泄漏系数,可取0.8;A ——裂口面积,m2 ;p ——两相混合物的压力,Pa;pc ——临界压力,Pa,可取pc=0.55Pa;ρ ——两相混合物的平均密度,kg/m3,它由下式计算:ρ=1/(Fv/ρ1+(1-Fv)/ρ2) (9)ρ1 —— 液体蒸发的蒸气密度,kg/m3 ;ρ2 —— 液体密度,kg/m3 ;Fv ——蒸发的液体占液体总量的比例,它由下式计算:Fv=Cp(T-Tc)/H (10)Cp ——两相混合物的定压比热,J/kg•K;T —— 两相混合物的温度,K;Tc—— 临界温度,K;H —— 液体的气化热,J/g。
当F>1时,表面液体将全部蒸发成气体,这时应按气体泄漏公式计算;如果Fv很小,则可近似按液体泄漏公式计算2.2.1.3 泄漏后的扩散如前所述,泄漏物质的特性多种多样,而且还受原有条件的强烈影响,但大多数物质从容器中泄漏出来后,都可发展成弥散的气团向周围空气扩散对可燃气体若遇到引火源会着火这里仅讨论气团圆形释放的开始形式,即液体泄漏后扩散、喷射扩散和绝热扩散关于气团在大气中扩散属环境保护范畴,在此不予考虑1) 液体的扩散液体泄漏后立即扩散到地面,一直流到低洼处或人工边界,如防火堤、岸墙等,形成液池液体泄漏出来不断蒸发,当液体蒸发速度等于泄漏速度时,液池中的液体量将维持不变如果泄漏的液体时是低挥发度的,则从液池中蒸发量较少,不易形成气团,对厂外人员没有危险;如果着火则形成池火灾;如果渗透进土壤,有可能对环境造成影响如果泄漏的时挥发性液体或低温液体,泄漏后液体蒸发量大,大量蒸发在液池上面后会形成蒸气云并对扩散到厂外,对厂外人员有影响1)液池面积如果泄漏的液体已达到人工边界,则液池面积即为人工边界围成的面积如果泄漏的液体未达到人工边界,则将假设液体的泄漏点为中心呈扁圆柱形在光滑表面上扩散,这时液池半径r用下式计算:瞬时泄漏(泄漏时间不超过30s)时: (11)连续泄漏(持续泄漏10min以上)时:r=(32gmt3/πp)1/4 (12)上述两式中:r——液池半径,m;m——泄漏的液体质量,kg;g——重力加速度,9.8m/s2;p——设备中液体压力,Pa;t——泄漏时间,s;(2)蒸发量液池内液体蒸发按其机理可分为闪蒸、热量蒸发和质量蒸发三种。
下面分别介绍①闪蒸:过热液体泄漏后,由于液体的自身热量而直接蒸发称为闪蒸发生闪蒸时液体蒸发速度Q1可由下式计算:Q1=Fv•m/t (13)式中Fv——直接蒸发的液体与液体总量的比例;m——泄漏的液体总量,kg;t——闪蒸时间,s②热量蒸发:当Fv<1或Qt<m时,则液体闪蒸不完全,有一部分液体在地面形成液池并吸收地面热量而气化称为热量蒸发热量蒸发速度Q1按下式计算: (14)式中 A1——液池面积,㎡;T0——环境温度,K;Tb——液体沸点,K;H——液体蒸发热,J/kg;L——液池长度,m;a——热扩散系数,㎡/s,见表;K——导热系数,J/mK,见表;t——蒸发时间,s;Nu——努舍尔特(Nusselt)数表示对流传热系数的准数Nu=C(Gr*Pr)n=C*λ/l*(CpβgΔtl3ρ2/[λμ])n= C*(Cpβg[T0-Ta]l3ρ2/[Kμ])n*λ/lβ——体积膨胀系数;K-1;Ta——液体温度;K;T0——地面温度;μ——粘度;mPasρ——液体密度;kg/m3表2-3 某些地面的热传递性质地面情况K/(Jm-1K-1)a/(㎡s-1)水泥土地(含水率8%)干涸土地湿地砂砾地1.10.90.30.62.51.2910-74.310-72.310-73.310-711.010-7表2-2-1 C、n取值范围(Gr*Pr)范围CN<104104-109109-10121.360.590.101/51/41/3③质量蒸发:当地面传热停止时,热量蒸发终止,转而由液池表面之上气流运动使液体蒸发,称为质量蒸发。
其蒸发速度Q1为: (15)式中 a——分子扩散系数,㎡/s;Sh——舍伍德(sherwood)数;性质:是反映包含有待定传质系数的无因次数群,类似于传热中的努塞特数,以符号Sh或Nsh表示它是由三个物理量组成,即Sh=k′L/DAB式中:k′为传质系数,m/s;L为特性尺寸,m;DAB为溶质A在溶剂中B中的特性系数,m2/s 时而A ——液池面积,m2;L ——液池长度,m;ρ1—— 液体的密度,㎏/m32)喷射扩散气体泄漏时从裂口喷出,形成气体喷射大多数情况下气体直接喷出后,其压力高于周围环境大气压力,温度低于环境温度在进行气体喷射计算时,应以等价喷射孔口直径计算等价喷射的孔口直径下式计算: (16)式中 D——等价喷射孔径,m;D0——孔口孔径,m;ρ0——泄漏气体的密度,㎏/m3;ρ——周围环境条件下气体的密度,㎏/m3如果气体泄漏能瞬时间达到周围环境的温度、压力状况,即ρ0=ρ,则D=D01)喷射的浓度分布 在喷射轴线上距孔口x处的气体浓度C(x)为: (17)式中 b1、b2——分布函数,其表达式如下:b1=50.5+48.2ρ-9.95ρ2b2=23+41ρ其余符号意义同前。
如果把式(17)改写成x是C(x)的函数形式,则给定某浓度值C(x),就可算出具有该浓度的点至孔口的距离x在过喷射轴线上点x且垂直于喷射轴线的平面内任一点处的气体浓度为: (18)式中 C(x,y)——距裂口距离x且垂直于喷口轴线的平面内Y点的气体浓度,㎏/m3;C(x)——喷射轴线上距裂口x处的气体浓度,㎏/m3;b2——分布参数,同前;y——目标点到喷射轴线的距离,m2)喷射轴线上的速度分布喷射速度随着轴线距离增大而减少,直到轴线上的某一点喷射速度等于风速为止,该点称为临界点临界点以后的气体运动不再符合喷射规律沿喷射轴线上的速度分布由下式得出: (19)式中ρ0——泄漏气体的密度,㎏/m3;ρ——周围环境条件下气体的密度,㎏/m3;D——等价喷射孔径,m;b1——分布参数,同前;x——喷射轴线上距裂口某点的距离,m;V(x)——喷射轴线上距裂口x处一点的速度,m/s;V0——喷射初速度,等于气体泄漏时流出裂口时的速度,m/s,按下式计算: (20)Q0——气体泄漏速度,㎏/s;Cd——气体泄漏系数;D0——裂口直径,m;当临界点处的浓度小于允许浓度(如可燃气体的燃烧下限或者有害气体最高允许浓度)时,只需按喷射来分析;若该点浓度大于允许浓度时,则需要进一步分析泄漏气体在大气中扩散的情况。
3)绝热扩散闪蒸液体或加压气体瞬时泄漏后,有一段快速扩散时间,假定此过程相当快以致在混合气团和周围环境之间来不及热交换,则称此扩散为绝热扩散根据TNO(1997年)提出的绝热扩散模式,泄漏气体(或液体闪蒸形成的蒸汽)的气团呈半球形向外扩散根据浓度分布情况,把半球分成内外两层,内层浓度均匀分布,且具有50%的泄漏量;外层浓度呈高斯分布,具有另外50%的泄漏量绝热扩散过程分为两个阶段,第一阶段气团向外扩散至大气压力,在扩散过程中气团获得动能,称为“扩散能”;第二阶段,扩散能再将气团向外推,使紊流混合空气进入气团,从而使气团范围扩大当内层扩散速度降到一定值时,可以认为扩散过程结束1)气团扩散能在气团扩散的第一阶段,扩散的气体(或蒸气)的内能一部分用来增加动能,对周围大气做功假设该阶段的过程为可逆绝热过程,并且是等熵的① 气体泄漏扩散能根据内能变化得出扩散能计算公式如下:E=CV(T1-T2)-0.98P0(V2-V1) (21)式中 E——气体扩散能,J;Cv——定容比热,J/㎏K;T1——气团初始温度,K;T2——气团压力降至大气压力时的温度,K;P0——环境压力,Pa;V1——气团初始体积,m3;V2——气团压力降至大气压力时的体积,m3。
② 闪蒸液体泄漏扩散能蒸发的蒸气团扩散能可以按下式计算:E=[H1-H2-Tb(S1-S2)]W-0.98(P1-P0)V1 (22)式中 E——闪蒸液体扩散能,J;H1——泄漏液体初始焓,J/㎏;H2——泄漏液体最终焓,J/㎏;Tb——液体的沸点,K;S1——液体蒸发前的熵,J/㎏K;S2——液体蒸发后的熵,J/㎏K;W——液体蒸发量,㎏;P1——初始压力,Pa;P0——周围环境压力,Pa;V1——初始体积,m31) 气团半径与浓度在扩散能的推动下气团向外扩散,并与周围空气发生紊流混合① 内层半径与浓度气团内层半径R1和浓度C是时间函数,表达式如下: (23) (24)式中 t——扩散时间,s;V0——在标准温度、压力下的气体体积,m3;Kd——紊流扩散系数,按下式计算: (25)如上所述,当中心扩散速度(dR/dt)降到一定值时,第二阶段才结束临界速度的选择是随机的而且不稳定的。
设扩散结束时扩散速度为1m/s,则扩散结束时内层半径R1和浓度C可按下式计算:R1=0.08837E0.3V01/3 (26)C=172.95E-0.9 (27)② 外层半径与浓度第二阶段末气团外层的大小可根据试验观察得出,即扩散终结时外层气团半径R2由下式求得:R2=1.456R1 (28)式中 R1、R2——分别为气团内层、外层半径,m外层气团浓度自内层向外层呈高斯分布2.2.2 火灾 易燃、易爆的气体及液体泄漏后遇到引火源就会被点燃而着火燃烧它们被点燃后的燃烧方式有池火、喷射火、火球和突发火四种2.2.2.1 池火可燃液体(如汽油、柴油等)泄漏后流到地面形成液池,或流到水面并覆盖水面,遇到火源燃烧而成池火1)燃烧速度当液池中的可燃液体的沸点高于周围环境温度时,液体表面上单位面积的燃烧速度dm∕dt为: (29)式中 dm/dt——单位表面积燃烧速度,kg∕m2s;HC——液体燃烧热,J/kg;Cp——液体的定压比热,J/kgK;Tb——液体的沸点,K;T0——环境温度,K;H——液体的气化热,J/kg。
当液体的沸点低于环境温度时,如加压液化气或冷冻液化气,其单位面积的燃烧速度dm/dt为: (30)式中符号意义同前燃烧速度也可从手册中直接得到下表列出了一些可燃液体的燃烧速度表2-4 一些可燃液体的燃烧速度物质名称汽油煤油柴油重油苯甲苯乙醚丙酮甲醇燃烧速度/(g﹒m2﹒s-1)92-8155.1149.3378.1165.37138.29125.8466.3657.62)火焰高度设液池为一半径为r的圆池子,其火焰高度可按下式计算: (31)式中 h——火焰高度;m;r——液池半径;m;——周围空气密度,kg∕m3;g——重力加速度,9.8m∕s2;dm∕dt——燃烧速度,kg/m2s3)热辐射通量液池燃烧时放出的总热辐射通量为: (32)式中 Q——总热辐射通量,W;η——效率因子,可取0.13~0.35;其余符号意义同前4)目标入射热辐射强度假设全部辐射热量由液池中心点的小球面辐射出来,则在距离池中心某一距离(X)处的入射热辐射强度为 (33)式中 I——热辐射强度,W/㎡;Q——总热辐射通量;W;tc——热传导系数,在无相对理想的数据时,可取值为1;X——目标点到液池中心距离,m;2.2.2.2 喷射火加压的可燃物质泄漏时形成射流,如果在泄漏裂口处被点燃,则形成喷射火。
这里所用的喷射火辐射热计算方法是一种包括气流效应在内的喷射扩散模式的扩散把整个喷射火看成是由沿喷射中心线上的全部点热源组成,每个点热源的热辐射通量相等点热源的热辐射通量按下式计算:q=ηQ0HC (34)式中 q——点热源辐射通量,W;η——效率因子,可取0.35;Q0——泄漏速度,㎏/s;HC——燃烧热,J/kg 从理论上讲,喷射火的火焰长度等于从泄漏口到可燃烧混合气燃烧下限(LFL)的射流轴线长度对表面火焰热通量,则集中在LFL/1.5处对危险评价分析而言,点热源数N点的划分可以是随意的,一般取n=5射流轴线上某点热源i到距离该点x处一点的热辐射强度为: (35)式中 Ii——点热源i至目标点x处的热辐射强度,W/m2;q——点热源的辐射通量,W;x——点热源到目标点的距离,m某一目标点处的入射热辐射强度等于喷射火的全部点热源对目标的热辐射强度的总和: (36)式中 n——计算时选取点热源数,一般取n=5。
2.2.2.3 火球和爆燃低温可燃液化气由于过热,容器内压增大,使容器爆炸,内容物释放并被点燃,发生剧烈的燃烧,产生强大的火球,形成强烈的热辐射1)火球半径R=2.665M0.327 (37)式中 R——火球半径,m;M——急剧蒸发的可燃物质的质量,kg2) 火球持续时间t=1.089M0.327 (38)式中 t——火球持续时间,s;3)火球燃烧时释放出辐射热通量 (39)式中 Q——火球燃烧时辐射热通量,W;Hc——燃烧热,J/kg;η——效率因子,取决于容器内可燃烧物质的饱和蒸气压p,η=0.27-0.32;其他符号同前4)目标接受到的入射热辐射强度 (40)式中 Tc——传导系数,保守取值为1;x——目标距火球中心的水平距离,m;其他符号同前2.2.2.4 固体火灾固体火灾的热辐射参数按点源模型估计。
此模型认为火焰射出的能量为燃烧的一部分,并且辐射强度与目标至火源中心距离的平方成反比,即: (41)式中 qr——目标接受到的辐射强度,W/m2;f——辐射系数,可取f=0.25;Mc——燃烧速度,kg/s;Hc——燃烧热,J/kg;x——目标至火源中心间的水平距离,m2.2.2.5 突发火泄漏的可燃气体、液体蒸发的蒸气在空气中扩散,遇到火源发生突然燃烧而没有爆炸此种情况下,处于气体燃烧范围内的室外人员将会全部死亡;建筑物内将有部分人被烧死突发火后果分析,主要是确定可燃混合气体的燃烧上、下极限的边界线及其下限随气团扩散到达的范围为此,可按气团扩散模型计算气团大小和可燃混合气体的浓度2.2.2.6火灾损失火灾通过辐射热的方式影响周围环境当火灾产生的热辐射强度足够大时,可使周围物体燃烧或变形,强烈的热辐射可能烧毁设备甚至造成人员伤亡等火灾损失估算建立在辐射通量与损失等级的相应关系的基础上,下表为不同入射通量造成伤害或损失的情况表2-5 热辐射的不同入射通量所造成的损失入射通量(kW﹒m-2)对设备的损害对人的伤害37.5操作设备全部损坏1%死亡/10s100%死亡/1min25在无火焰、长时间辐射下,木材燃烧的最小能量重大烧伤/10s100%死亡/1min12.5有火焰时,木材燃烧,塑料熔化的最低能量1度烧伤/10s1%死亡/1min4.020s以上感觉疼痛,未必起泡1.6长期辐射无不舒服感从表中可看出,在较小辐射等级时,致人重伤需要一定的时间,这时人们可以逃离现场或掩蔽起来。
2.2.3爆炸2.2.3.1 简述爆炸是物质的一种非常急剧的物理、化学变化,也是大量能量在短时间内迅速释放或急剧转化成机械功的现象它通常是借助于气体的膨胀来实现从物质运动的表现形式来看,爆炸就是物质剧烈运动的一种表现物质运动急剧增速,由一种状态迅速转变成另一种状态,并在瞬间内释放出大量的能1)爆炸的特征一般说来,爆炸现象具有以下特征:(1) 爆炸过程进行得很快;(2) 爆炸点附近压力急剧升高,产生冲击波;(3) 发出或大或小的响声;(4) 周围介质发生震动或邻近物质遭受破坏一般将爆炸过程分为两个阶段:第一阶段是物质的能量以一定形式(定容、绝热)转变为强压缩能;第二阶段强压缩能急剧绝热膨胀对外做功,引起作用介质变形,移动和破坏2)爆炸的类型按爆炸性质可分为物理爆炸和化学爆炸物理爆炸就是物质状态参数(温度、压力、体积)迅速发生变化,在瞬间放出大量能量并对外做功的现象其特点是在爆炸现象发生过程中,造成爆炸发生的介质的化学性质不发生变化,发生变化的仅是介质的状态参数例如锅炉、压力容器爆炸和各种气体或液化气体钢瓶的超压爆炸以及高温液体金属遇水爆炸等化学爆炸就是物质由一种化学结构迅速转变为另一种化学结构,在瞬间放出大量能量并对外做功的现象。
如可燃气体、蒸气或粉尘与空气混合形成爆炸性混合物的爆炸化学爆炸的特点是:爆炸发生过程中介质的化学性质发生了变化,形成爆炸的能源来自物质迅速发生化学变化时所释放的能量化学爆炸有三个要素,所反应的放热性、反应的快速性和生成气体产物 雷电是一种自然现象,也是一种爆炸从工厂爆炸事故来看,有以下几种化学爆炸类型:(1) 蒸气云团的可燃混合气体遇火源突然燃烧,是在无限空间中的气体爆炸;(2) 受限空间内可燃混合气体的爆炸;(3) 化学反应失控或工艺异常所造成压力容器爆炸;(4) 不稳定的固体或液体爆炸总之,发生化学爆炸时会释放出大量的化学能,爆炸影响范围较大;而物理爆炸仅释放出机械能,其影响范围较小2.2.3.2物理爆炸的能量物理爆炸如压力容器破裂时,气体膨胀所释放的能量(即爆破能量)不仅与气体压力和容器的容积有关,而且与介质在容器内的物性相态有关因为有的介质以气态存在,如空气、氧气、氢气等;有的以液体存在,如液氨、液氯等液化气体、高温饱和水等容积与压力相同而相态不同的介质,在容器破裂时产生的爆破能量也不同,而且爆炸过程也不完全相同,其能量计算公式也不同1)压缩气体与水蒸气容器爆破能量当压力容器中介质为压缩气体,即以气态形式存在而发生物理爆炸时,其释放的爆破能量为: (a)式中 Eg——气体的爆炸能量,kJ;p——容器内气体的绝对压力,MPa;V——容器的容积,m3;k——气体的绝热指数,即气体的定压比热与定容比热之比。
常用气体的绝热指数数值见下表表2-6 常用气体的绝热指数气体名称空气氮氧氢甲烷乙烷乙烯丙烷一氧化碳k值1.41.41.3971.4121.3161.181.221.331.395气体名称二氧化碳一氧化氮二氧化氮氨气氯气过热蒸气干饱和蒸气氢氰酸k值1.2951.41.311.321.351.31.1351.31从表中可看出,空气、氮、氧及一氧化氮、一氧化碳等气体的绝热指数均为1.4或近似1.4,若用k=1.4带入式(a)中,则: (b)令 则式(b)可简化为:Eg=CgV (c)式中 Cg——常用压缩气体爆破能量系数,kg/m3压缩气体爆破能量Cg是压力p的函数,各种常用压力下的气体爆破能量系数列于下表中表2-7 常用压力下气体容器爆破能量系数(k=1.4时)表压力p/MPa0.20.40.60.81.01.62.5爆破能量系数Cg/(kJ﹒m-3)21024.61027.51021.11031.41032.41033.9103表压力p/MPa4.05.06.415.03240爆破能量系数Cg/(kJ﹒m-3)6.71038.61031.11042.71046.51048.2104若将k=1带入式(a),可得干饱和蒸气容器爆破能量为: (d)用上式计算有较大的误差,因为没有考虑蒸气干度的变化和其他的一些影响,但它可以不用查明蒸气热力性质而直接计算,因此可供危险性评价参考。
对于常用压力下的干饱和蒸气容器的爆破能量可按下式计算:Es=CsV (e)式中 Es——水蒸气的爆破能量,kJ;V——水蒸气的体积,m3;Cs——干饱和水蒸气爆破能量系数,kJ/ m3各种常用压力下的干饱和水蒸气容器能量系数列于下表中表2-8 常用压力下干饱和水蒸气容器爆破能量系数表压力p/MPa0.30.50.81.32.53.4爆破能量系数4.371028.311021.51032.751036.241037.771032)介质全部为液体时爆破能量通常用液体加压时所做的功为常温液体压力容器爆炸时释放的能量,计算公式如下: (f)式中 EL——常温液体压力容器爆炸时释放的能量,kJ;p ——液体的压力(绝),Pa-1;V ——容器的体积,m3;t——液体在压力p和温度T下的压缩系数,Pa-13)液化气体与高温饱和水的爆破能量液化气体和高温饱和水一般在容器内以气液两态存在,当容器破裂发生爆炸时,除了气体的急剧膨胀做功外,还有过热液体激烈的蒸发过程。
在大多数情况下,这类容器内的饱和液体占有容器介质重量的绝大部分,它的爆破能量比饱和气体大得多,一般计算时考虑气体膨胀做的功过热状态下液体在容器破裂时释放出爆破能量可按下式计算: (g)式中 E——过热状态液体的爆破能量,kJ;H1——爆炸前饱和液体的焓,kJ/kg;H2——在大气压力下饱和液体的焓,kJ/kg;S1——爆炸前饱和液体的熵,kJ/kg﹒℃;S2——在大气压力下饱和液体的焓,kJ/kg﹒℃;T1——介质在大气压力下的沸点,kJ/kg﹒℃;W——饱和液体的质量,kg饱和水容器的爆破能量按下式计算:EW=CSV (h)式中 EW——饱和水容器的爆破能量,kJ;V——容器内饱和水所占的体积,m3;Cs——饱和水爆破能量系数,kJ/ m3,其值见下表表2-9 常用压力下饱和水爆破能量系数表压力p/MPa0.30.50.81.32.53.4Cs(kJ/ m3)2.381043.251044.561046.351049.561041.061042.2.3.3 爆炸冲击波及其伤害、破坏作用压力容器爆炸时,爆破能量在向外释放时以冲击波能量、碎片能量和容器残余变形能量三种形式表现出来。
后两者所消耗的能量只占总爆破能量的3%~15%,也就是说大部分能量是产生空气冲击波1)爆炸冲击波冲击波是由压缩波叠加形成的,是波阵面以突进形式在介质中传播的压缩波容器破裂时,容器内的高压气体大量冲出,使它周围的空气受到冲击波而发生扰动,使其状态(压力、密度、温度等)发生突跃变化,其传播速度大于扰动介质的声速,这种扰动在空气中的传播就成为冲击波在离爆破中心一定距离的地方,空气压力会随时间发生迅速而悬殊的变化开始时,压力突然升高,产生一个很大的正压力,接着又迅速衰减,在很短时间内正压降至负压如此反复循环数次,压力渐次衰减下去开始时产生的最大正压力即是冲击波波阵面上的超压△p多数情况下,冲击波的伤害、破坏作用是由超压引起的超压△p可以达到数个甚至数十个大气压冲击波伤害、破坏作用准则有:超压准则、冲量准则、超压-冲量准则为了便于操作,下面仅介绍超压准则超压准则认为,只要冲击波超压达到一定值时,便会对目标造成一定的伤害或破坏超压波对人体的伤害和对建筑物的破坏作用见下表表2-10 冲击波超压时对人体的伤害作用超压△p/MPa伤害作用0.02~0.03轻微损伤0.03~0.05听觉器官损伤或骨折0.05~0.10内脏严重损伤或死亡>0.10大部分人员死亡表2-11 冲击波超压时对建筑物的破坏作用超压△p/MPa破坏作用超压△p/MPa破坏作用0.005~0.006门、窗玻璃部分破碎0.06~0.07木建筑厂房房柱折断,房架松动0.006~0.015受压面的门窗玻璃大部分破碎0.07~0.10砖墙倒塌0.015~0.02窗框损坏0.10~0.20防震钢筋混凝土破坏,小房屋倒塌0.02~0.03墙裂缝0.20~0.30大型钢架结构破坏0.04~0.05墙大裂缝,屋瓦掉下2)冲击波的超压冲击波波阵面上的超压与产生冲击波的能量有关,同时也与距离爆炸中心的远近有关。
冲击波的超压与爆破中心距离的关系为:△p∝R-n (i)式中 △p——冲击波波阵面上的超压,MPa;R ——距爆炸中心的距离,m;n ——衰减系数衰减系数在空气中随着超压的大小而变化,在爆炸中心附近内为2.5~3;当超压在数个大气压以内时,n=2;小于1个大气压时,n=1.5实验数据表明,不同数据的同类炸药发生爆炸时,如果距离爆炸中心的距离R之比与炸药量q三次方根之比相等,则所产生的冲击波超压相同,用公式表示如下:若 (j)则 △p=△p0式中 R——目标与爆炸中心距离,m;R0——目标与基准爆炸中心的相当距离,m;q0——基准爆炸能量,TNT,kg;q ——爆炸时产生冲击波所消耗的能量,TNT,kg;△p——目标处的超压,MPa;△p0——基准目标处的超压,MPa;α ——炸药爆炸试验的模拟比上式也可写成为:△p(R)=△p0(R/α) (k)利用式(k)就可以根据某些已知药量的试验所测得的超压来确定任意药量爆炸时在各种相应距离下的超压。
下表是1000kgTNT炸药在空气中爆炸时所产生的冲击波超压表2-12 1000kgTNT爆炸时所产生的冲击波超压距离R0/m56789101214超压△p0/MPa2.942.061.671.270.950.760.500.33距离R0/m1618202530354045超压△p0/MPa0.2350.170.1260.0790.0570.0430.0330.027距离R0/m505560657075超压△p0/MPa0.02350.02050.0180.0160.01430.013综上所述,计算压力容器爆破时对目标的伤害、破坏作用,可按下列程序进行1)首先根据容器内所装介质的特性,分别选用式(b)至式(h)计算出其爆破能量E2)将爆破能量q换算成TNT当量q因为1kgTNT爆炸所放出的爆破能量为4230~4836kJ/kg,一般取平均爆破能量为4500kJ/kg,故其关系为:q=E/qTNT=E/4500 (l)(3)按式(2-51)求出爆炸的模拟比α,即:α=(q/q0)1/3=(1/1000)1/3=0.1q1/3 (m)(4)求出在1000kgTNT爆炸试验中相当距离R0,即R0=R/α。
5)根据R0值在表中找出距离为R0处的超压△p中间值用插入法),此即所求距离为R处的超压6)根据超压△p的值,从表中找出对人员和建筑物的伤害、破坏作用2.2.3.4蒸气云爆炸的冲击波伤害、破坏半径爆炸性气体以液体储存,如果瞬间泄漏后遇到延迟点火或气态储存时泄漏到空气中,遇到火源,则可能发生蒸气云爆炸导致形成蒸气云的力来自容器内含有的能量或可燃物含有的内能,或两者兼而有之能”主要形式是压缩能、化学能或热能一般来说,只有压缩能和热量才能单独导致形成蒸气云根据荷兰应用研究院(TNO(1997))建议,可按下式预测蒸气云爆炸的冲击波的损害半径:R=Cs(NE)1/3 (Ⅰ)式中 R——损害半径,m;E——爆炸能量,kJ,可按下式计算:E=VHc (Ⅱ)V——参与反应的可燃气体的体积,m3;Hc——可燃气体的高燃烧热值,kJ/m3,取值情况见表;N——效率因子,其值与燃烧浓度持续展开所造成损耗的比例和燃料燃烧所得机械能的数量有关,一般取N=10%Cs——经验常数,取决于损害等级,其取值情况见表。
表2-13 某些气体的高燃烧热值(kJ/m3)气体名称高热值气体名称高热值氢气12770乙烯64019氨气17250乙炔58985苯47843丙烷101828一氧化碳17250丙烯94375硫化氨生成SO225708正丁烷134026生成SO330146异丁烷132016甲烷39860丁烯121883乙烷70425表2-14 损害等级表损害等级Cs(mJ-1/3)设备损坏人员伤害10.03重创建筑物和加工设备1%死亡于肺部伤害>50%耳膜破裂>50%被碎片击伤20.06损坏建筑物外表,可修复性破坏1%耳膜破裂1%被碎片击伤30.15玻璃破碎被碎玻璃击伤40.410%玻璃破碎2.2.4 中毒有毒物质泄漏后生成有毒蒸气云,它在空气中飘移、扩散,直接影响现场人员并可能波及居民区大量剧毒物质泄漏可能带来严重的人员伤亡和环境污染毒物对人员的危害程度取决于毒物的性质、浓度和人员与毒。