1、火灾事故后果分析火灾事故后果分析概述概述历史上著名的火灾历史上著名的火灾火灾事故后果分析研究目的火灾事故后果分析研究目的 火灾事故是工业重大危险事故里的主要类型，对人的生命、火灾事故是工业重大危险事故里的主要类型，对人的生命、财产和环境有巨大的危害。财产和环境有巨大的危害。火灾火灾事故后果事故后果分析分析作为安全评价的一作为安全评价的一种类型，对事故后果进行分析及预测，为预测火灾、爆炸和中毒种类型，对事故后果进行分析及预测，为预测火灾、爆炸和中毒事故造成的后果提供准确有效的数据，大大缩短了灾害后果评估事故造成的后果提供准确有效的数据，大大缩短了灾害后果评估时间，提高灾害后果评估精度，对预防、减
2、少人员伤亡和财产损时间，提高灾害后果评估精度，对预防、减少人员伤亡和财产损失有积极影响失有积极影响，因此对火灾事故进行事故后果因此对火灾事故进行事故后果分析是十分有必要分析是十分有必要的的。火灾的分类火灾的分类 火灾时或失去控制蔓延的一种灾害性燃烧现象，是各种灾害中火灾时或失去控制蔓延的一种灾害性燃烧现象，是各种灾害中发生最频繁且最具毁灭性的一种。火灾的种类很多，通常包括森林、发生最频繁且最具毁灭性的一种。火灾的种类很多，通常包括森林、建筑、油类以及可燃气体和粉尘爆炸所引起的火灾。这里主要针对建筑、油类以及可燃气体和粉尘爆炸所引起的火灾。这里主要针对易燃易爆的气体或液体泄漏后遇到引火源被点燃而
3、着火燃烧而引起易燃易爆的气体或液体泄漏后遇到引火源被点燃而着火燃烧而引起火灾和固体类物质起火成灾进行火灾事故后果分析。火灾和固体类物质起火成灾进行火灾事故后果分析。由此，我们将火灾分为池火灾、喷射火、火球、突发火和固体由此，我们将火灾分为池火灾、喷射火、火球、突发火和固体火灾五类。火灾五类。第一章第一章 池火灾池火灾 由于油罐破裂或运输意外,造成泄漏事故,当遇到点火源,极可能发生池火灾。灾害一旦发生,人员暴露在火灾的热辐射下,将会受到严重的伤害。池火灾通过热辐射方式影响周围环境，当热辐射强度足够大时，可使周围物体燃烧或变形，强烈的热辐射可能烧毁设备，并造成人员伤亡。池火灾实例 1999年，山西
4、省某化工厂5000m3直径为28.5米的汽油储罐发生火灾，顶部塌陷后发展为池火灾，池液半径为12.25 m。顿时周围一定范围内笼罩在强烈的热辐射中。此事故造成16人死亡，百余人受伤，直接经济损失5千万元，是一个典型的池火灾事故。池火灾事故后果分析方法 目前，对池火灾事故后果分析的方法大体有两类：一类是经验公式计算法 这种计算方法方便实用，简单易行，可以规律性地获取相关数据，便于掌握和预测油罐燃烧情况的总体情况。一类是FDS软件模拟法 FDS(Fire Dynamics Simulator)是美国国家技术标准局（NIST）建筑火灾实验室开发的基于场模拟的火灾模拟软件，该软件采用先进的大涡模拟技术
5、，得到众多实例验证，在火灾安全工程领域应用广泛。经验公式计算法步骤1.1.1 燃烧速度的计算 由于燃烧速度极不稳定且不断变化，因此取燃烧中期平均燃烧速度为准。不同直径油罐液体表面上单位面积的燃烧速度为：式中：Vb直线燃烧速度，mm/s；D油罐直径，m。直线燃烧速度与质量燃烧速度的转换利用公式：式中：V质量燃烧速度，kg/(m2s)；液体体积质量，kg/m3。计算结果为：单位面积燃烧速度为0.0254 kg/(m2s)1.1.2 火焰高度计算 当液池为一个半径r的圆形池时，火焰高度按下述公式计算：式中：h火焰高度，m；r液池半径，m，半径为12.25 m；0周围空气的体积质量，取值1.16 kg
6、/m3；g重力加速度，取值9.8 m/s2。计算结果为：火焰高度可达22.41 m（在大型火灾实验中，火焰形状并不清晰，因为在火焰周围存在大量的烟，所以上述火焰高度计算仅是理想情况下的估算值）。1.1.3 液池燃烧时的总辐射量 液池燃烧时的总辐射量按下述公式计算：式中：Q总辐射量，W；效率因子，介于0.13至0.35之间;(根据试验数据，对于大型油盘火灾，由于烟的屏蔽作用，使得辐射分数随着油盘直径的增加而下降，该计算中取平均值0.24);Hc最大发热量，J/mol，取43 728.8 J/mol。计算结果为：由该式计算得总热辐射量Q=61.9103 kW。1.1.4 目标入射热辐射强度的计算
7、假设全部辐射热量都由液池中心释放出来，在距离液池中心点X（m）的某一点的入射热辐射强度(目标入射热辐射强度)为：式中：I入射热辐射强度，W/m3；tc热传导系数，在无相对理想数据时，可取值1；X目标点到火源（罐壁）的距离，m。上述目标入射热辐射强度公式，反映了热辐射能量与受害目标到火源距离之间的关系。当入射热辐射通量是一个定值时，可以计算出受害目标与火源之间的伤害距离，参照辐射危害程度（表1），计算得到：当I1=37.5 kW/m2时,X1=11.4m；I2=25.0 kW/m2时,X2=14.1m；I3=12.5 kW/m2时,X3=19.9m；I4=4.0kW/m2时,X4=35.1m；I
8、5=1.6 kW/m2时,X5=55.5 m。1.1.5 对比热辐射危害程度表热辐射计算所得热辐射曲线 FDS采用NIST的FDS火灾模拟器通过大涡模拟来进行地面油罐火灾发展过程的可视化仿真。选取油罐为中心的100 m3空间进行模拟，持续燃烧20s，取以油罐中心横截面来观察热辐射以及数据导出，经Smokeview截图可看出火焰不停变动但总体趋势稳定，故取中心截面火焰1/2高度处515 s的平均值为导出数据（选取相同距离热辐射最大值），Smokeview模拟视图见图13：图1、2是FDS模拟油罐燃烧热辐射的三维图，图2是包含烟气情况下的辐射图。图3为热辐射截面图。1.2 FDS模拟热辐射计算图1
9、 模拟油罐燃烧热辐射的三维图图2 包含烟气情况下的模拟油罐燃烧热辐射的三维图图3 热辐射截面图模拟热辐射计算所得热辐射曲线经验公式计算法与FDS软件模拟计算法之比较 两组数据对比可看出，计算数据更具有规律性，模拟数据具有多变性和不对称性，但二者对于最大热辐射强度的数值点的走向上是趋于一致的；通过对几个重要数值点据在大部分时线完全拟合。结论：a)传统计算方法方便实用，可以规律性地获取相关数据，便 于掌握和预测油罐燃烧情况的总体情况。b)FDS模拟对于火焰的动态变化情况更具实际意义，同时在数值上对经验公式起到补充和完善作用。c)通过两种方法的结合，可以更全面准确地掌握油罐燃烧的热辐射强度数据，更准
10、确地预测油罐火灾的燃烧特性，对实战具有重要的指导意义。表1 辐射危害程度一览案例分析案例分析 本例中危险单元为输油管道，且无防护本例中危险单元为输油管道，且无防护堤。假定泄漏的液体无蒸发、已充分蔓延且堤。假定泄漏的液体无蒸发、已充分蔓延且地面无渗透。泄漏量为地面无渗透。泄漏量为1000kg，环境温度为，环境温度为30。原油的性质为：比重。原油的性质为：比重0.850.89；相对；相对密度密度0.7800.970（液）；闪点是（液）；闪点是6.6732.22；自燃点；自燃点350；爆炸极限；爆炸极限1.16.4%（V）；）；沸点沸点300325；火焰温度；火焰温度1100；热值；热值41870k
11、J/kg。第二章 喷射火 喷射火事故后果分析中存在大量的不确定性喷射火事故后果分析中存在大量的不确定性,分析困难较大。本分析困难较大。本文是运用喷射火焰对目标的入射热辐射强度计算模型文是运用喷射火焰对目标的入射热辐射强度计算模型,确定喷射火事确定喷射火事故灾害损失面积故灾害损失面积,再根据城市街区平均人口密度和平均财产密度确定再根据城市街区平均人口密度和平均财产密度确定事故中人员伤亡数量和财产损失数量。在目前技术和管理水平下这事故中人员伤亡数量和财产损失数量。在目前技术和管理水平下这种方法是一种较好的方法种方法是一种较好的方法,并能为其他相关事故后果分析提供思路。并能为其他相关事故后果分析提供
12、思路。城市燃气管道喷射火事故实例城市燃气管道喷射火事故实例 2003年，山东省潍坊市某输送天然气的中压年，山东省潍坊市某输送天然气的中压A级管道发生破裂，级管道发生破裂，遇明火被点燃后发展为喷射火灾事故遇明火被点燃后发展为喷射火灾事故,造成了一定量的人员伤亡和财造成了一定量的人员伤亡和财产损失。经测定，管道内燃气压力为产损失。经测定，管道内燃气压力为0.4MPa,外界大气压为外界大气压为1.01105Pa,温度为常温，裂口为圆形温度为常温，裂口为圆形,且面积约为且面积约为0.4cm2。喷射火事故后果分析步骤2.1 事故分析前假设 为了便于计算，首先我们把把整个喷射火看成是由沿喷射中心为了便于计
13、算，首先我们把把整个喷射火看成是由沿喷射中心线上的所有点热源组成，每个点热源的热辐射通量相等。（一般将线上的所有点热源组成，每个点热源的热辐射通量相等。（一般将喷射火焰看成由沿喷射中心线上的喷射火焰看成由沿喷射中心线上的5个点热源组成）个点热源组成）2.2 燃气泄漏速度的计算 当满足当满足 时气体流动属于声速流动，其泄漏速度为：时气体流动属于声速流动，其泄漏速度为：当满足当满足 时，气体呈亚声速流动，其泄漏速度为：时，气体呈亚声速流动，其泄漏速度为：备注：以上两式中p0为环境压力,Pa;p为管道内燃气压力,Pa;k为气体绝热指数（即定压比热CP与定容比热CV之比;Q0为泄漏速度,kg/s;Cd
14、为气体泄漏系数（当裂口形状为圆形时取1.00,三角形时取0.95,长方形时取0.90;M为气体摩尔质量,kg/Kmol;T为气体温度,K;A为裂口面积,m2;R为气体常数,J/molK。根据上述计算公式,则气体流动属于声速流动,在常温下其泄漏速度为Q0=0.868Kg/s2.3 热辐射强度的计算 在将喷射火焰看成由沿喷射中心线上的在将喷射火焰看成由沿喷射中心线上的5个点热源组成的前提个点热源组成的前提下。某点热源下。某点热源i到距离该点到距离该点x处某一目标点的热辐射强度为：处某一目标点的热辐射强度为：式中：式中：为点热源为点热源 i至目标点至目标点z处的热辐射强度，处的热辐射强度，W/m2；
15、为辐射率，一般取为辐射率，一般取0.2；x为点热源到目标点的距离，为点热源到目标点的距离，m。备注：备注：由于将喷射火假设有由于将喷射火假设有5个点热源，所以某点的入射热辐射强度个点热源，所以某点的入射热辐射强度等于喷射火的全部点热源对目标的热辐射强度的总和。同时又由于等于喷射火的全部点热源对目标的热辐射强度的总和。同时又由于喷射火焰长度、形状受到喷射速度、喷射口形状和大小、风速等多喷射火焰长度、形状受到喷射速度、喷射口形状和大小、风速等多种因素影响。喷射火焰对目标的热辐射强度计算存在大量的不确定种因素影响。喷射火焰对目标的热辐射强度计算存在大量的不确定性，为简化计算性，为简化计算,假设喷射火
16、沿喷射中心线的全部点火源集中在某一假设喷射火沿喷射中心线的全部点火源集中在某一点处点处,如下图所示。如下图所示。经计算，热辐射通量为经计算，热辐射通量为48 275.80J；所以单个热源的热辐射强度与；所以单个热源的热辐射强度与距离的关系式为距离的关系式为 。27 6 8.7Ix254qIx喷射火简化模型喷射火简化模型 所以，距离假设点热源所以，距离假设点热源x处的目标的入射热辐射强度即为：处的目标的入射热辐射强度即为：总辐射强度与距离的关系是为：总辐射强度与距离的关系是为：。25768.7Ix2.4 伤害面积的计算伤害面积的计算 火灾通过热辐射的方式影响周围环境火灾通过热辐射的方式影响周围环
17、境,当火灾产生的热辐射强度当火灾产生的热辐射强度足够大时足够大时,可使周围的物体燃烧或变形可使周围的物体燃烧或变形,强烈的热辐射可能烧毁设备强烈的热辐射可能烧毁设备甚至造成人员伤亡甚至造成人员伤亡,表表1中给出了不同热辐射通量造成不同的损失。中给出了不同热辐射通量造成不同的损失。根据表根据表1中不同损失入射通量的阙值中不同损失入射通量的阙值,结合下述公式可以确定喷射火结合下述公式可以确定喷射火造成不同灾害的半径和面积。造成不同灾害的半径和面积。表1 辐射危害程度一览伤害半径伤害半径：伤害面积：伤害面积：式中：式中：I为热辐射强度，为热辐射强度，W/m2；q为点热源热辐射通为点热源热辐射通量，量
18、，W；为辐射率，一般取为辐射率，一般取0.2。O代表喷射火源；代表喷射火源；区域区域SR代表代表“操作设备全部损坏，人员操作设备全部损坏，人员1死亡死亡l0s、100死亡死亡1min”；区域区域So代表代表“在无火焰、长时间辐射下，木材燃烧；人员重大损伤在无火焰、长时间辐射下，木材燃烧；人员重大损伤l0s、100死亡死亡1min”；区域区域SY代表代表“有火焰时，木材燃烧，塑料熔化；人员有火焰时，木材燃烧，塑料熔化；人员1度烧伤度烧伤10s、l死死亡亡1min”。喷射火火灾危险区域的划分 喷射火事故后果分析结果为：喷射火事故后果分析结果为：以喷射火源为中心，以以喷射火源为中心，以10.12米为
19、半径内的米为半径内的SR=321.84m2的的区域内操作设备全部损坏，人员区域内操作设备全部损坏，人员1死亡死亡l0s、100死亡死亡1min；在以喷射火源为中心，以；在以喷射火源为中心，以10.12米到米到12.4米为半径米为半径SO=160.92m2的环形范围内人员重大损伤的环形范围内人员重大损伤l0s、100死亡死亡1min,在一喷射火源为中心，以在一喷射火源为中心，以12.4米到米到17.54米为半径的米为半径的SY=482.76m2的环形区域内人员的环形区域内人员1度烧伤度烧伤10s、l死亡死亡1min。2.5 损失损失估算估算 城市燃气输送不同于长输燃气输送城市燃气输送不同于长输燃
20、气输送,城市燃气管道遍布城市每个街区城市燃气管道遍布城市每个街区,其人口其人口密度和财产密度很大密度和财产密度很大,事故一旦发生事故一旦发生,极有可能造成重特大人员伤亡和财产损失。极有可能造成重特大人员伤亡和财产损失。因此对城市燃气管道喷射火灾害进行损失估计具有重要的实际意义。由于城市因此对城市燃气管道喷射火灾害进行损失估计具有重要的实际意义。由于城市人口流动很大人口流动很大,交通状况随时间和街区的不同而不同交通状况随时间和街区的不同而不同,例如商业文化街、火车站、例如商业文化街、火车站、汽车站等区域人员密度很大汽车站等区域人员密度很大;时间上时间上,在上下班高峰期街道人员密度很大在上下班高峰
21、期街道人员密度很大,而夜晚而夜晚人员密度则很小。所以城市街区人员密度不易确定人员密度则很小。所以城市街区人员密度不易确定,但可以针对街区不同的实际但可以针对街区不同的实际情况通过成立调查小组进行实地调查确定。调查街区人员密度可分三个时间段情况通过成立调查小组进行实地调查确定。调查街区人员密度可分三个时间段进行进行:上下班时间上下班时间(7 00-9 00、11 00-13 00、17 00-19 00)、白天其他时、白天其他时间间(5 00-7 00、9 00-11 00、13 00-17 00、19 00-21 00)、夜晚时间、夜晚时间(21 00-5 00)。最后可以将调查结果统计处理。
22、最后可以将调查结果统计处理,得到各个时间段街区的人员密得到各个时间段街区的人员密度。而街区管道周边的财产密度也可以调查得到度。而街区管道周边的财产密度也可以调查得到,不同的是财产密度趋于稳定不同的是财产密度趋于稳定,受受不确定性因素影响较小。值得注意的是财产密度应包括路边停放和路上行驶的不确定性因素影响较小。值得注意的是财产密度应包括路边停放和路上行驶的车辆。车辆。假设某市某街区经调查三个时间段的平均人员密度分别为假设某市某街区经调查三个时间段的平均人员密度分别为:1、1、1(人人/m2),财产密度为财产密度为2(元元/m2)。根据下列公式可以计算喷。根据下列公式可以计算喷射火事故财产损失程度
23、和人员伤亡数量射火事故财产损失程度和人员伤亡数量,具体公式如下具体公式如下:111121131RYNSNSoNS111121131RYNSNSoNS1 111 211 31RYNSNS oNS1222ROLSLS N11、N11、N11分别表示三个时间段分别表示三个时间段“人员人员1%死亡死亡/10s、100%死亡死亡/1min”;N12、N12、N12分别表示三个时间段分别表示三个时间段“人员人员重大损伤重大损伤/10s、100%死亡死亡/1min”;N13、N13、N13分别表示三个分别表示三个时间段时间段“人员人员1度烧伤度烧伤/10s、1%死亡死亡/1min”。但在一般情况下。但在一般
24、情况下人员可以及时远离喷射火焰人员可以及时远离喷射火焰,不会造成死亡。不会造成死亡。L1、L2分别表示全部财产、部分财产损失。分别表示全部财产、部分财产损失。经对此事故发生后的调查，此事故发生在上午经对此事故发生后的调查，此事故发生在上午10点，不属于上点，不属于上班高峰期，再加上事故发生地点远离闹市区，没有造成人员的重大班高峰期，再加上事故发生地点远离闹市区，没有造成人员的重大伤亡，仅有部分人员轻度烧伤。事故发生地的财产密度为伤亡，仅有部分人员轻度烧伤。事故发生地的财产密度为10万元万元/m，所以最终得出全部损失的财产为所以最终得出全部损失的财产为3218.4万元，部分损失的财产为万元，部分
25、损失的财产为1609.2万元。万元。喷射火事故后果分析总结喷射火事故后果分析总结 喷射火事故后果分析受到大量不确定性因素的影响喷射火事故后果分析受到大量不确定性因素的影响,例如管道裂例如管道裂口面积、喷射火焰长度和形状、平均人员密度、平均财产密度、气口面积、喷射火焰长度和形状、平均人员密度、平均财产密度、气候条件等。尽管这些不确定因素对后果分析造成偏差候条件等。尽管这些不确定因素对后果分析造成偏差,但在目前技术但在目前技术和管理水平下本文所提出的后果分析计算方法不失是一种较好的方和管理水平下本文所提出的后果分析计算方法不失是一种较好的方法法,可以为城市燃气管道风险评估提供重要依据。可以为城市燃
26、气管道风险评估提供重要依据。案例计算案例计算 利用该模型计算甲烷泄漏所造成的喷射利用该模型计算甲烷泄漏所造成的喷射火灾害后果：火灾害后果：估计泄漏口的孔径为估计泄漏口的孔径为0.02m，管道内绝，管道内绝对压力为对压力为15个大气压，折合个大气压，折合1.52106Pa。环境温度为环境温度为300K，相对湿度，相对湿度50%，火焰温度，火焰温度为为2000K。甲烷分子量为。甲烷分子量为16.04，空气分子量，空气分子量为为29，甲烷在空气中含量约为，甲烷在空气中含量约为9.8%。火球 装有过热可燃液体的储存罐或工艺容器突然失效装有过热可燃液体的储存罐或工艺容器突然失效,将导致液体泄将导致液体泄
27、露。由于储罐或容器中的压力远远大于大气压露。由于储罐或容器中的压力远远大于大气压,内部的液体将快速内部的液体将快速喷到空中喷到空中,并有一部分气体几乎在瞬间变为气体并有一部分气体几乎在瞬间变为气体,并且迅速膨胀并且迅速膨胀,形形成一个由可燃液体蒸气、没有气化的液体的液滴以及空气组成的蒸成一个由可燃液体蒸气、没有气化的液体的液滴以及空气组成的蒸气云。此蒸气云被点燃即会形成火球。气云。此蒸气云被点燃即会形成火球。火球迅速增大到最大火球迅速增大到最大,由于火球释放的热量使其内部的液体蒸发由于火球释放的热量使其内部的液体蒸发而导致火球变轻并离开地面而导致火球变轻并离开地面,直到将其内部燃料烧完。在短短
28、的燃烧直到将其内部燃料烧完。在短短的燃烧过程中过程中,火球可以释放出大量的辐射能火球可以释放出大量的辐射能,而且会在比它本身大好几倍而且会在比它本身大好几倍的范围内引起人员伤亡和财产损失。因此的范围内引起人员伤亡和财产损失。因此,在对工艺容器或过热可燃在对工艺容器或过热可燃液体的储存罐进行风险分析的时候液体的储存罐进行风险分析的时候,建立一个精确的火球热辐射后果建立一个精确的火球热辐射后果模型模型,以对其进行精确的定性定量分析具有重要的现实意义。以对其进行精确的定性定量分析具有重要的现实意义。火球事故后果分析方法火球事故后果分析方法 目前，对火球的分析方法包括两类：目前，对火球的分析方法包括两
29、类：一类是一类是静态模型分析法静态模型分析法 此方法将火球简单地看为一件静态事件此方法将火球简单地看为一件静态事件,一个有固定的大小一个有固定的大小,在一个在一个固定位置的火球以一个固定的速率释放能量，而忽略了火球产生过程中固定位置的火球以一个固定的速率释放能量，而忽略了火球产生过程中的变化，从而使其分析结果产生偏差。的变化，从而使其分析结果产生偏差。一类是一类是动态模型分析法动态模型分析法 此方法考虑到了火球本身所具有的产生、燃烧扩大、抬升以及最后此方法考虑到了火球本身所具有的产生、燃烧扩大、抬升以及最后消失的动态过程，从而能够更加准确的对火球事故进行分析。消失的动态过程，从而能够更加准确的
30、对火球事故进行分析。火球事故实例火球事故实例 某液化石油气码头包括某液化石油气码头包括3000 T 级液化石油气专用码头和级液化石油气专用码头和12 台容台容积为积为100 m3 的卧式液化石油气贮罐。液化石油气是由丙烷、丁烷等的卧式液化石油气贮罐。液化石油气是由丙烷、丁烷等碳氢化合物加压液化的一种高热值碳氢化合物加压液化的一种高热值(45144 45980 kJ/kg)的燃料。经的燃料。经测定本港区的液化石油气密度为测定本港区的液化石油气密度为0.16 kg/L。主要工艺设计参数。主要工艺设计参数:设计设计压力压力11765 MPa,设计温度设计温度50 oC,充满系数充满系数85%,安全泄
31、压阀定值安全泄压阀定值116 MPa。由于某次操作失误，储罐区的。由于某次操作失误，储罐区的1 个个100 m3 储罐发生泄漏（泄漏储罐发生泄漏（泄漏百分比为百分比为10%）并引发火球）并引发火球,对周围人员及财产造成了重大的伤害。对周围人员及财产造成了重大的伤害。火球静态模型分析法火球静态模型分析法3.1.1 火球火球半径的计算公式为半径的计算公式为 式中：式中：R火球半径（火球半径（m）；）；W火球中消耗的可燃物质能量。对于单罐储存，火球中消耗的可燃物质能量。对于单罐储存，W取罐取罐容量的容量的50%；对于双罐储存，；对于双罐储存，W 取罐容量的取罐容量的70%；对于多罐储存，；对于多罐储
32、存，W 取罐容量的取罐容量的90%。132.9RW3.1.2 火球持续时间火球持续时间的计算公式的计算公式 式中：式中：t火球持续时间；火球持续时间；W火球中消耗的可燃物质的量。火球中消耗的可燃物质的量。130.45tW3.1.3 目标目标接受到的热辐射通量的计算公接受到的热辐射通量的计算公式式 根据上式，就可以计算出火球周围地区不同距离处的热辐射通量值，从而确定根据上式，就可以计算出火球周围地区不同距离处的热辐射通量值，从而确定热辐射通量的分布。热辐射通量的分布。式中：式中：q(r)对目标的热辐射通量（对目标的热辐射通量（kW/m2）；）；E火球表面辐射通量（火球表面辐射通量（kW/m2）。
33、对于柱状、卧式和立式储罐，）。对于柱状、卧式和立式储罐，E 取取270 kW/m2，对于球罐，对于球罐，E 取取200 kW/m2；R火球半径；火球半径；r目标与评价对象之间的地面水平距离目标与评价对象之间的地面水平距离。32222()(10.058 ln)/(),q rER rrRrrR3.1.4 目标接受到的热量的计算公式为 式中：Q(r)目标接受到的热量（J/m2）；q(r)目标接受到的热辐射通量。由上式可以看出，目标接受到的热量与其接收到的热辐射通量和暴露时间均成正比。暴露时间越长，则其接收到的热量越多，受到的损伤也越大。()()Qrqrt火球动态模型分析法步骤140.9dtM Mar
34、tinsen 和和Marx对对TNO 的经验公式稍做修正之后的经验公式稍做修正之后,提出提出火球持续时间火球持续时间td 为为 式中式中 M 为火球内燃料质量为火球内燃料质量,kg。带入数据可以计算出火球的持续时间为带入数据可以计算出火球的持续时间为7.61s。3.2.1 火球火球持续时间的计算持续时间的计算火球火球直径的计算直径的计算 火球直径在增大阶段持续变大火球直径在增大阶段持续变大,大概在大概在t=td/3 时直径达到最大时直径达到最大值值,然后开始离地升空然后开始离地升空,一直到消散直径都保持最大值不变一直到消散直径都保持最大值不变,计算公计算公式为：式为：在火球扩大阶段，火球直径与
35、实践的关系为：在火球扩大阶段，火球直径与实践的关系为：在火球的稳定阶段，火球的直径为：在火球的稳定阶段，火球的直径为：99.84m。113413m a x1()8.6 6 4;031()5.8;3dddDtMtttDtDMttt13()7 3.2 2Dtt3.2.2 火球中心到地面的高度火球中心到地面的高度 火球在放大阶段保持在地面上火球在放大阶段保持在地面上,而一旦火球直径达到最大时火球就开而一旦火球直径达到最大时火球就开始上升始上升,即在稳定燃烧阶段内即在稳定燃烧阶段内,火球中心以固定的速率从火球中心以固定的速率从 处上升到离地面处上升到离地面 处。处。火球中心到地面的高度的计算公式为：火
36、球中心到地面的高度的计算公式为：m ax2Dm ax32Dm a x()1();02331();23ddddDtHtttDtHttttt 火球在增大阶段火球中心离地面的高度为：火球在增大阶段火球中心离地面的高度为：火球在稳定燃烧阶段火球中心离地面的高度为：火球在稳定燃烧阶段火球中心离地面的高度为：()19.68;2.547.61H ttsts13()36.61;02.54H ttsts3.2.3 目标接受到的热通量的计算目标接受到的热通量的计算 目标接受到的热通量是目标接受到的热通量是火球表面热辐射能火球表面热辐射能E(t),火球和目标之间火球和目标之间的的几何视角系数几何视角系数F(x,t)
37、和和大气传输率大气传输率S(x,t)的函数，即为的函数，即为:(,)(,)(,)()IxTSxtFxtEt火球表面辐射能火球表面辐射能 火球表面的热辐射随时间变化趋势是在火球的增长阶段保持火球表面的热辐射随时间变化趋势是在火球的增长阶段保持最大值不变最大值不变,而且是在而且是在 内以固定的速率由最大值内以固定的速率由最大值减小到减小到0,即即:式中：式中：HC为燃烧热为燃烧热,kJ/kg;Emax为表面辐射能最大值为表面辐射能最大值,kW/m2；f为为燃烧热辐射系数，其计算公式为燃烧热辐射系数，其计算公式为 。11 21()m ax0.0 1 3 3;0331()m ax(1);23cdddd
38、E tEfHMtttE tEtttt13ddttt0.3 20.2 7BfP 代入数据可以得出燃烧代入数据可以得出燃烧热辐射系数热辐射系数f为为0.32，火球表面最大热辐射能，火球表面最大热辐射能Emax为为398.6kW/m2，在火球的增大阶段，火球的表面热辐射能为，在火球的增大阶段，火球的表面热辐射能为()597.978.57;2.547.61E ttsts几何视角系数几何视角系数 目标可能在火球周围的各个方向目标可能在火球周围的各个方向,即目标不可能接受到来自于辐即目标不可能接受到来自于辐射表面上每个点的辐射能射表面上每个点的辐射能,而只能接受到其中一部分而只能接受到其中一部分,因此因此
39、,要引入要引入几何视角系数这个参数。几何视角系数就是每单位面积上目标接受几何视角系数这个参数。几何视角系数就是每单位面积上目标接受的和火球释放的辐射能之比的和火球释放的辐射能之比,换一种说法就是目标能看到的火球的换一种说法就是目标能看到的火球的范围。它不仅与火球的大小范围。它不仅与火球的大小,火球与目标的距离有关火球与目标的距离有关,还与火球与目还与火球与目标的相对方位有关。标的相对方位有关。火球在增大阶段的几何视角为：火球在增大阶段的几何视角为：火球在稳定燃烧阶段的几何视角为：火球在稳定燃烧阶段的几何视角为：232231340.3(,);02.541340.3tF x tststx22249
40、2(,);2.547.61387.3F x tststx大气传输率大气传输率 大气传输率大气传输率(S)为空气中水蒸气和为空气中水蒸气和CO2 对热辐射的吸收率对热辐射的吸收率,其计其计算公式为：算公式为：式中式中 R 为相对湿度为相对湿度;PV 为一定温度下水的饱和蒸汽压为一定温度下水的饱和蒸汽压,Pa。0.0922()(,)2.02()2VD tS x tRPHtx 火球在增大阶段的大气传输率为：火球在增大阶段的大气传输率为：火球在稳定燃烧阶段的大气传输率为：火球在稳定燃烧阶段的大气传输率为：2120.0933(,)1.01 1340.336.61;02.54S x ttxtsts220.
41、09(,)1.01 387.349.92;2.547.61S x ttxsts 图图2 100 m3 的液化石油气储罐的液化石油气储罐10%泄漏所导致火球的泄漏所导致火球的D(t)、H(t)、E(t)。火球在增大阶段和稳定燃烧阶段的热通量之和为：火球在增大阶段和稳定燃烧阶段的热通量之和为：0()(,)(,)(,)dtdoseIxSx t Fx t Ex t dt23212.547.6120.09220.093322202.54(1340.336.61)(597.978.57)(387.349.92)387.31340.3ttxtttxdtdttxx表表2 CCPS静态模型与动态模型对于静态模型
42、与动态模型对于100 m3 液化石油液化石油气贮罐泄漏产生火球的热剂量随距离变化关系气贮罐泄漏产生火球的热剂量随距离变化关系目标距离x/mIdose(x)KW/m2CCPS静态模型动态模型10033819220010857500179100042表表3 热辐射热辐射剂量反应关系剂量反应关系热剂量/(KW/m2)对暴露目标所造成的伤害40皮肤灼疼100一度烧伤150二度烧伤（起水泡）250三度烧伤（1%死亡）500三度烧伤（50%死亡）1200三度烧伤（99%死亡）CCPS静态模型静态模型与动态模型之比较与动态模型之比较 我们这里对我们这里对CCPS静态模型和动态模型所得到的热辐射后静态模型和动
43、态模型所得到的热辐射后果进行了比较果进行了比较,这里计算了当这里计算了当x=100 m,200 m,500 m,1000 m 时时,动态模型和动态模型和CCPC静态模型所得到的热通量的值。静态模型所得到的热通量的值。从表从表2、表、表3中可以看出中可以看出,在相同的位置在相同的位置,CCPS静态模型得到静态模型得到的热剂量明显大于动态模型。比如，在的热剂量明显大于动态模型。比如，在100 m 处处,静态模型引起静态模型引起的是三度烧伤的是三度烧伤,而动态模型引起的是二度烧伤而动态模型引起的是二度烧伤;在在200 m 处处,静态静态模型引起的是一度烧伤模型引起的是一度烧伤,而动态模型引起的是皮肤
44、灼痛。而动态模型引起的是皮肤灼痛。案例计算案例计算 某运输甲烷的管道发生泄漏并着火燃烧某运输甲烷的管道发生泄漏并着火燃烧，泄漏甲烷的体积为，泄漏甲烷的体积为20立方米。已知甲烷的立方米。已知甲烷的密度为密度为716.8kg/m3，相对湿度为，相对湿度为50%。突发火（闪火）突发火（闪火）液体、气体泄漏后在泄漏源附近扩散，在泄源上方形成气团，液体、气体泄漏后在泄漏源附近扩散，在泄源上方形成气团，气团将在大气中进一步扩散，当遇到火源时发生突然燃烧而没有爆气团将在大气中进一步扩散，当遇到火源时发生突然燃烧而没有爆炸，此种火灾为突发火。在此种情况下，处于气体燃烧范围内的室炸，此种火灾为突发火。在此种情
45、况下，处于气体燃烧范围内的室外人员将会全部烧死，建筑物内的将有部分人员烧死。外人员将会全部烧死，建筑物内的将有部分人员烧死。突发火的主要危害来自热辐射和火焰直接接触。可燃物云团的突发火的主要危害来自热辐射和火焰直接接触。可燃物云团的大小决定了可能造成直接火焰接触危害的面积，而云团的大小则部大小决定了可能造成直接火焰接触危害的面积，而云团的大小则部分取决于扩散和泄漏条件。分取决于扩散和泄漏条件。突发火事故后果分析步骤4.1 突发火火焰高度的计算突发火火焰高度的计算 在突发火模型中，假定突发火是一个火焰以恒定速度传播的过程。火焰可在突发火模型中，假定突发火是一个火焰以恒定速度传播的过程。火焰可视高
46、度可用近似的半经验公式计算：视高度可用近似的半经验公式计算：式中：式中：H为火焰可是高度，为火焰可是高度，m；d为云团厚度，为云团厚度，m；s为燃烧速度，为燃烧速度，m；为空气密度，为空气密度，kg/m3;为燃气混合物的密度，为燃气混合物的密度，kg/m3；r为理想配比下空气与燃料的质量比。为理想配比下空气与燃料的质量比。1223302 0(/)(/)2(/(1)aHdsg dw rwa04.2 火焰火焰高度的计算高度的计算 火焰宽度火焰宽度W随时间变化关系为：随时间变化关系为：式中：式中：R为云层半径，为云层半径，m；S为燃烧速度，为燃烧速度，m/s;t为时间，为时间，s。12222()WR
47、RSt4.3 几何视角系数的计算几何视角系数的计算 假设辐射面和接受面是两个互相平行的平面，则假设辐射面和接受面是两个互相平行的平面，则F可用可用Fmax表示，表示，计算关系如下：计算关系如下：220.5m ax()hvFFF11tan(1/)tan()/2hrrFXAXA11tan+tan/2VrrBFHAABH 式中：式中：H为火焰可是高度，为火焰可是高度，m；X为到目标物的距离，为到目标物的距离，m；W为火焰宽度，为火焰宽度，m。220.5()rrAHX20.5(1)rrBHX/rHHb/rXXb12bW4.4 大气传导率的计算大气传导率的计算 在保守计算场合，若是干燥晴朗的天气一般可取
48、在保守计算场合，若是干燥晴朗的天气一般可取 =1，已知湿度，已知湿度时可用下式计算：时可用下式计算：式中：式中：RH为相对湿度；为相对湿度；X为到目标物的距离，为到目标物的距离，m。a0.1080.13log(14.1)aRHX4.5 目标所受热辐射能的计算目标所受热辐射能的计算 平面物体单位面积上接收的辐射能由下式计算：平面物体单位面积上接收的辐射能由下式计算：式中：式中：E为辐射能，为辐射能，KW/m2；F为几何视角系数；为几何视角系数；为大气传导率。为大气传导率。aqE Fa突发火事故后果分析总结突发火事故后果分析总结 突发火由于其发生的突然性和隐蔽性等特点容易被人们所忽略，突发火由于其
49、发生的突然性和隐蔽性等特点容易被人们所忽略，并且此类火灾一旦发生可能造成的损失普遍偏大。并且易燃气体泄并且此类火灾一旦发生可能造成的损失普遍偏大。并且易燃气体泄露点与火灾发生地点可能距离比较远，从而为事故的预防和控制造露点与火灾发生地点可能距离比较远，从而为事故的预防和控制造成很大困难。成很大困难。经过上述计算，在我们可以得出经过上述计算，在我们可以得出平面物体单位面积上接收的辐平面物体单位面积上接收的辐射能然后再对照射能然后再对照 表表1 中的物体和人员破坏和损伤阈值就可以分析出中的物体和人员破坏和损伤阈值就可以分析出突发火的危害半径，从对突发火的危害半径，从对提高灾害后果评估精度，预防、减
50、少人员提高灾害后果评估精度，预防、减少人员伤亡和财产损失伤亡和财产损失有一定帮助。有一定帮助。表表1 1 辐射危害程度一览辐射危害程度一览案例分析案例分析 在一个空旷的地面短时间泄漏了大量的丙烷，在一个空旷的地面短时间泄漏了大量的丙烷，在传播扩散中形成一个扁平的圆形的云团。云团约在传播扩散中形成一个扁平的圆形的云团。云团约有有1.2m厚，直径为厚，直径为80m，云团内燃料的体积百分比，云团内燃料的体积百分比为为15%，从云团边缘点火燃烧。当时风速，从云团边缘点火燃烧。当时风速1m/s，相，相对湿度为对湿度为40%，已知丙烷分子量为，已知丙烷分子量为44kg/kmol，空，空气气29kg/kmo
