重大危险源分级标准

1、重大危险源分级标准（征求意见稿）1 适用范围本规范规定了重大危险源评估分级的方法和程序。本规范为重大危险源评估分级技术规范，适用于包括储罐区、库区、生产场所等重大危险源。2 规范性引用文件下列文件中的条款，通过本规范的引用而成为本标准的条款。凡是标注日期的引用文件，其随后所有的修改单（不包括勘误的内容）或修订版均不适用于本规范，然而，鼓励根据本规范达成协议的各方研究是否可使用这些文件的最新版本。凡是不注日期的引用文件，其最新版本适用于本规范。中华人民共和国安全生产法危险化学品安全管理条例安全生产许可证条例重大危险源辨识（GB18218）安全评价通则关于规范重大危险源监督与管理工作的通知（安监总

2、协调字2005125号）3 术语和定义卜列术语和定义适用于本规范重大危险源majorhazardinstallations重大危险源是指长期地或者临时地生产、搬运、使用或者储存危险物品，且危险物品的数量等于或超过临界量的单元（包括场所和设施）。4重大危险源分级判据重大危险源分级判据如表1所示。表1重大危险源分级判据危险源等级分级判据死亡人数一级重大危险源可能造成30人（含30人）以上二级重大危险源可能造成1029人三级重大危险源可能造成39人四级重大危险源可能造成1-2人具体判别的依据如下:一级重大危险源：可能造成死亡30人（含30人）以上的重大危险源;二级重大危险源：可能造成死亡10-29人

3、的重大危险源;三级重大危险源：可能造成死亡3-9人的重大危险源;四级重大危险源:可能造成死亡1-2人的重大危险源。5重大危险源死亡人数及财产损失计算方法可能造成的死亡人数评价程序为：将重大危险源的周边区域划分成等间隔的网格区，用一笛卡尔坐标体系的网格覆盖城市的区域地图（如图1所示），网格间距大小取决于当地人口密度，以不影响计算结果为准。确定每一网格内的人员数量，通过火灾（室内火灾除外）、爆炸、毒物泄漏扩散事故后果模型计算重大危险源事故在每一网格中心处产生的热辐射、超压或毒物浓度的数值，然后通过热辐射、冲击波超压、中毒概率函数将其其转化为造成死亡的概率。将每一网格中心的死亡率与人口数量相乘，即得

4、到死亡的人数。将所有网格的死亡人数求和，即得到总的死亡人数。具体用下式表小：(1)nDiSVii1式中，N为总的死亡人数，D为第i个网格的人口密度，S为网格面积，Vi为第i个网格的个人死亡率，n为网格的数目。FeHRlzen industry: Petroleum refinery iFMflfe 1: 十皿PeSroleunS re-finiryiRfBanffukhadlM扁：Y,二 We| lendb-wDbmt& 后二 dChemical industry图1死亡人数计算原理示意图采用财产损失半径的方法评估事故后果造成的损失，并假定此半径内没有损失的财产与此半径外损失的财产相互

5、抵消，或者说此半径内的财产完全损失，此半径外的财产完全无损失。财产损失半径通过火灾、爆炸事故后果模型确定。财产损失半径按下式计算:，1/3Ri-1KiWTNT21763175WTNT(2)式中，R为i区半径，miK为常量。热辐射对建筑物的影响直接取决于热辐射强度的大小及作用时间的长短，以引燃木材的热通量作为对建筑物破坏财产损失半径，计算公式如下：q6730t4/525400(3)tW/Mc(4)式中，q为引燃木材的热通量(W/m),t为热辐射作用时间，即火灾持续时问(s)。6重大危险源评价分级程序重大危险源的评价分级程序如下图所示。如果一种危险物质具有多种事故形态，按照后果最严重的事故形态考虑

6、，即遵循“最大危险原则”。各类重大危险源具体事故情景选择、后果计算及死亡概率计算过程参见附录A。图2重大危险源评价分级程序附录A:重大危险源事故后果模型储罐区重大事故后果分析A.1.1储罐区的主要事故后果类型A.1.1.1池火灾易燃液体如汽油、苯、甲醇、乙酸乙酯等，一旦从储罐及管路中泄漏到地面后，将向四周流淌、扩展，形成一定厚度的液池，若受到防火堤、隔堤的阻挡，液体将在限定区域（相当于围堰）内得以积聚，形成一定范围的液池。这时，若遇到火源，液池可能被点燃，发生地面池火灾。A.1.1.2蒸气云爆炸易燃易爆气体如H2、天然气等，泄漏后随着风向扩散，与周围空气混合成易燃易爆混合物，在扩散扩过程中如遇

7、到点火源，延迟点火，由于存在某些特殊原因和条件，火焰加速传播，产生爆炸冲击波超压，发生蒸气云爆炸。易燃易爆的液化气体如液化石油气、液化丙烷、液化丁烷等，其沸点远小于环境温度，泄漏后将会由于自身的热量、地面传热、太阳辐射、气流运动等迅速蒸发，在液池上面形成蒸气云，与周围空气混合成易燃易爆混合物，并且随着风向扩散，扩散扩过程中如遇到点火源，也会发生蒸气云爆炸。A.1.1.3喷射火对于易燃易爆气体如H2、天然气，以及易燃易爆的液化气体来说，泄漏后可能因摩擦产生的静电立即点火，产生喷射火。A.1.1.4沸腾液体扩展蒸气云爆炸易燃易爆的液化气体容器在外部火焰的烘烤下可能发生突然破裂，压力平衡被破坏，液体

8、急剧气化，并随即被火焰点燃而发生爆炸，产生巨大的火球。这种事故被称为沸腾液体扩展为蒸气云爆炸。A.1.1.5中毒事故毒性的液化气体如液氯、液氨等，由于沸点小于环境温度，泄漏后会因自身热量、地面传热、太阳辐射、气流运动等迅速蒸发，生成有毒蒸气云，密集在泄漏源周围，随后由于环境温度、地形、风力和湍流等因素影响产生漂移、扩散，范围变大，浓度减小。A.1.2储罐区主要事故后果模型A.1.2.1池火灾事故后果模型池火灾火焰的几何尺寸及辐射参数按如下步骤计算。计算池直径根据泄漏的液体量和地面性质，按下式可计算最大可能的池面积。S亚此访（1）式中，S为液池面积（m）,W为泄漏液体的质量（kg）,为液体的密度

9、（kg/m3）Hin为最小油层厚度（m）0最小物料层厚度与地面性质对应关系见表1。表1不同性质地面物料层厚度表地面性质最小物料层厚度（mD草地粗糙地面平整地面混凝土地面平静的水面确定火焰高度计算池火焰高度的经验公式如下：hLD42mf/（0，gD）0.61（2）式中：L为火焰高度（m）,D为池直径（mi）,m为燃烧速率（kg/m2s）P0为空气密度（kg/m3）,g为引力常数。计算火焰表面热通量假定能量由圆柱形火焰侧面和顶部向周围均匀辐射，用下式计算火焰表面的热通量:qo0.25D2HCmff0.25D2DL(3)式中，qo为火焰表面的热通量（kw/m2）,AH为燃烧热（kJ/kg）,冗为圆周

10、率，f为热辐射系数（可取为）,m为燃烧速率（kg/m2s）,其它符号同前目标接收到的热通量的计算目标接收到的热通量q（r）的计算公式为：q（r）q°（10.058lnr）V（4）式中，q（r）为目标接收到的热通量（kw/m2）,q0为由式（3）计算的火焰表面的热通量（kw/m2）,r为目标到油区中心的水平距离（V为视角系数。视角系数的计算角系数V与目标到火焰垂直轴的距离与火焰半径之比s,火焰高度与直径之比h有关。V,（VH2）VHAB1A（b1/s）tan，、1B（a1/s）tanVVtan1（h/（s2J2a0.5tana10.50.51)0.5)/sh(JKtan0.5(5)(6

11、)/(b21)0.5/(a21)0.5K)/s0.5(7)(8)(9)(11)(11)22-a(hs1)/(2s)(b(1s2)/(2s)(13)其中A、RJ、K、VH、VV是为了描述方便而引入的中间变量，冗为圆周率。A.1.2.2蒸气云爆炸事故后果模型蒸气云爆炸产生的冲击波超压是其主要危害。冲击波超压可通过传统的TNT当量系数法进行计算，将事故爆炸产生的爆炸能量等同于一定当量的TNT；也可根据爆炸能量直接计算。（1）TNT当量法确定闪蒸系数在热力学数据资料的基础上，用下式估算燃料的闪蒸部分。CpTF1expL(14)式中，F为蒸发系数，C为燃料的平均比热（kJ/kgK）,AT为环境压力下容器

12、内温度与沸点的温差（K）,L为汽化热（kJ/kg）o计算云团中燃料的质量:Wf2FW(15)式中，W为云团中燃料的质量（kg）,W为泄漏的燃料的质量（kg）,F为闪蒸系数。计算TNT当量：WTNTeWfHf/Htnt（16）式中，WNt为燃料的TNT当量（kg）,W为云团中燃料的质量（kg）,H为燃料的燃烧热（MJ/kg）,Hnt为io*1O2无因次距离mkg-1/3.WtntD1无因次超压IX图1APsR曲线TNT的爆热(MJ/kg),ae为TNT当量系数，推荐y。将实际距离转化为无因次距离：Rr/WtNt(17)式中，R为离爆炸点的实际距离(m),R为无因次距离(m)。在离爆炸点距离为R处

13、，根据相应的R值，查图1得到超压，进而预测人员受伤害和建筑受破坏的情况。(2)直接计算法在得到云团中燃料的质量的情况下，可按下式直接计算爆炸冲击波超压Ap23ln(ps/pa)0.91261.5058(lnZ)0.1675(lnZ)0.0320(lnZ)(18)(WZ012)zR/(马13(19)PaE1.8WQC(20)式中,Aps为冲击波正相最大超压(Pa),Z为无量纲距离，Pa为环境压力，R为目标到爆源的水平距离(项，E为爆源总能量(J),a为蒸气云当量系数，一般取，W为蒸气云中对爆炸冲击波有实际贡献的燃料质量(Kg),Q为燃料的燃烧热(J/Kg)。A.1.2.3喷射火事故后果模型加压的

14、可燃物泄漏时形成射流，如果在泄漏裂口处被点燃，则形成喷射火。假定火焰为圆锥形，并用从泄漏处到火焰长度4/5处的点源模型来表示。火焰长度计算喷射火的火焰长度可用如下方程得到:(21)0444L(Hem).161.66式中，L为火焰长度（n）,HC为燃烧热（J/kg）,m为质量流速（kg/s）热辐射的通量计算距离火焰点源为X（m）处接收到的热辐射通量可用下式表示:(22)fHCmq24X21000式中，q为距离X处接收的热辐射的通量（KW/市，f为热辐射率，p为大气传输率。(23)(24)大气传输率p按下式计算：10.0565InXA.1.2.4沸腾液体扩展为蒸气云爆炸事故后果模型计算主要包括如下

15、步骤火球直径D2.665W0.327式中，D为火球直径（m）,W为火球中消耗的可燃物质量（Kg）。对单罐储存,W取罐容量的50%对双罐储存，W取罐容量的70%对多罐储存，W取罐容量的90%火球持续时间t1.089W0.327（25）式中，t为火球持续时间（s）,W同式（24）。火球抬升高度火球在燃烧时，将抬升到一定高度。火球中心距离地面的高度H由下式估计:HD（26）火球表面热辐射能量SEP由下式计火球表面热辐射能量假设火球表面热辐射能量是均匀扩散的_2SEPFsmHa/(Dt)(27)式中，Fs为火球表面辐射的能量比，H为火球的有效燃烧热(J/Kg)。Fs与储罐破裂瞬间储存物料的饱和蒸气压力

16、P(MPa有关：Fs0.27P0.32(28)对于因外部火灾引起的BLEVE*故，上式中的P值可取储罐安全阀启动压力Pv(MPa的倍，即：P1.21Pv(29)H由下式求得：HaHcHvCpT(30)式中，H为燃烧热(J/kg),T为常沸点下的蒸发热(J/kg),Cp为恒压比热(J/),T为火球表面火焰温度与环境温度之差(K),一般来说T=1700K视角系数视角系数F的计算公式如下：2，一一F(D/2)/r)(31)式中，r为目标到火球中心的距离(m)0令目标与储罐白水平距离为X(M,则：r(X2H2)0.5(32)大气热传递系数火球表面辐射的热能在大气中传输时，由于空气的吸收及散射作用，一部

17、分能量损失掉了。假定能量损失比为a,则大气热传递系数a1oa和大气中的CO和H2O的含量、热传输距离及辐射光谱的特性等因素有关。a可由以下的经验公式来求取：a2.02(pwr)0.09(33)式中，国为环境温度下空气中的水蒸气压(N/m2),r为目标到火球表面的距离(mj)。PwpW/RH(34)0式中，Pw为环境温度下的饱和水蒸气压(N/m),RH为相对湿度。r'rD/2(35)火球热辐射强度分布函数在不考虑障碍物对火球热辐射产生阻挡作用的条件下，距离储罐X处的热辐射强度q(W/m)可由下式计算：qSEPFa(36)A.1.2.5中毒事故后果模型(1)泄漏模型液体泄漏速率模型液体泄漏

18、可根据流体力学中的柏努力方程计算泄漏量。当裂口不规则时，可采取等效尺寸代替；当泄漏过程中压力变化时，则往往采用经验公式。柏努力方程如下：c2(PPo)QCdA,2gh,(37)式中，Q为液体泄漏速率(kg/s),Cd为无量纲泄漏系数，是液体密度(kg/m3),A是泄漏孔面积(m2),P为罐压(Pa),R为大气压力(Pa),g为引力常数(9.8m/s2),h为液压高度(m)。液体出口速度可按下式计算：QuCdA(38)式中，u为液体出口速度(m/s),其他符号如前。持续时间按下式计算：tsuo/(Cdg)(At/A)(39)式中，Uo为初始流速（m/s）,Ar为罐内液面积（m2）。泄漏系数C的取

19、值通常可从标准化学工程手册中查到。对于管道破裂，Cd的典型取值为。表2为常用的液体泄漏系数数据。表2液体泄漏系数Cd雷诺数Re裂口形状圆形（多边形）三角形长方形>100<100这个方法没有考虑泄漏速率对时间的依赖关系（压力随时间而降低以及液压高度下降）。因此，计算出的泄漏速率是保守的最大可能泄漏速率。气体泄漏模型压力气体泄漏通常以射流的方式发生，泄漏的速度与其流动的状态有关，其特征可用临界流（最大出口速度等于声速）或亚临界流来描述。Perry等人用如下的关系式作为临界流的判断准则：当式（40）成立时，气体流动属音速流动；当式（41）成立时，气体流动属亚音速流动(42)k(40)Po

20、21Pk1k(41)Po2此Pk1式中，Po为环境大气压力（Pa）,P为容器压力（Pa）,k为气体的绝热指数,即定压比热G和定容比热C之比。对于很多气体，临界比值（P/P。）c近似为2,也就是说储压近似等于大气压力的两倍，此时流体泄漏的出口速度近似等于声速。临界流的质量泄漏速率可按下式计算:QCdAPk1Mk2k1RTk1气体呈亚音速流动时，其泄漏量为:k1Mk2k1QYCdAPj（43）RTk1式中，Q是气体泄漏速率（kg/s）,Cd为气体泄漏稀疏，A为裂口面积（m2）,M是气体相对分子质量，R是普适气体常数（Jmol-1K1）,T是气体的储存温度（K）,Y为气体膨胀因子，按式（44）计算2

21、Pk1Pok1PokP(44)上述考虑的为理想气体的不可逆绝热扩散过程。止匕外，没有考虑气体泄漏速率随时间的变化，因此使用初始储存条件必然导致保守的结果。两相流泄漏模型Cude在1975年建议了两相流泄漏关系式。假设源容器和泄漏点之间的管道长度和管道直径之比L/D>12泄漏点压力与泄漏点上流压力之比R/P=。具体计算方法如下：第一步，按下式计算两相流的质量分数：(TTc)CpH(45)式中，M为蒸发的液体占液体总量的比例，Tc是对应于泄漏点压力R的平衡温度（K）,T是对应于泄漏点上流压力P的平衡温度（K）,G是液体的定压比热J/（kgK），T是液体的蒸发热（J/kg）。第二步，按下式计算

22、两相流的平均密度：1MV1MVV1(46)式中，、v和1分别是两相流、蒸气和液体的密度（kg/m3）第三步，按下式计算两相流的质量泄漏速率Q（kg/S）:QACda2(PPc)(47)式中，Cd为泄漏系数，多数情况下，取Cd=,A为裂口面积（m2）,P为两相混合物的压力（Pa）,PC为临界压力（Pa）。如果L/D<12,先按前面介绍的方法计算纯液体泄漏速率和两相流泄漏速(48)43）计算出的两相率，再用内插法加以修正。两相流实际泄漏速率的计算公式为:QQvi(QiQvi)(12L/D)/10式中，Q、Qvi和Qi分别为两相流实际泄漏速率、按式（流泄漏速率和纯液体泄漏速率（kg/S）。如果

23、L/D<2,般认为泄漏为纯液体泄漏。（3）非重气云扩散模型瞬间泄漏扩散模型2QC(x, y,z,t,He)(2 )3/2 -exp 喂 expexpx y z(z H)22 Zx (x Vt)22y2 2乙y(49)连续泄漏扩散模型C(x,y,z,t,He)2Q'exo2 V y zexo(zH)22 2exp(z H)22 Z(50)式中，C为气云中危险物质浓度(kg/mi),He为泄漏源有效高度（源瞬间泄漏量（kg）,Q'为源连续泄漏速率（kg/s）,V为风速（m）/s）,t为泄漏后的时间（s） , xy和z分别为x、y和z方向的扩散系数（m）对于连续泄漏，平均时间取

24、iOmin。其中（Tx,（Ty,（Tz与地面的有效粗糙度有关。地面有效粗糙度长度如下表所示。表3地面有效粗糙度长度表地回类型Zo/m地回类型Zo/m草原、平坦开阔地<农作物地区村落、分散的树林1分散的高矮建筑物(城市)14密集的局矮建筑物(大城市)4有效粗糙度ZoW0.1m地区的扩散参数按下表选取。表4ZoW0.1m地区的扩散参数大气稳定度0-y/m0-z/mA/人-1/2(1+B/人-1/2(1+C/人-1/2(1+/人-1/2(1+D-1/2(1+-1/2(1+E(1+-1/2(1+-1/2F/A-1/2(1+(1+-1/2有效粗糙度乙方0.1m的粗糙地形扩散系数为：yyofyzz0

25、fzfy(Zo)1aoZoInxfz(x,Zo)(boColnx)(doeolnx)1Zof0g0式中，(Tyo、(Tzo按表4中的数值取值。其他系数按表5取值。表5不同大气稳定度下的系数值稳定度ABCDEFaoboCodoeofogo式（49）和式（50中泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状时气云中心的离地高度。在大多数问题中，泄漏源有效高度难以与泄漏源实际高度相一致。事实上，它等于泄漏源实际高度加泄漏源抬升高度。泄漏源抬升高度可以用下面的公式近似计算：_1,、HVsd1.50.268Pa（TsTa）Tsd/V（51）或H2.4Vsd/V（52）式中：AH是泄漏源抬升高度（m

26、,Vs是气云出口速度（用7s）,d是出口直径（m）,V是环境风速（m/s）,pa是环境大气压力（Pa）,Ts是气云出口温度（K）,Ta是环境大气温度（K）。计算出泄漏源抬升高度以后，将泄漏源抬升高度与泄漏源实际高度相加就得到了泄漏源有效高度。（3）重气云扩散模型常用模型有盒子模型和平板模型两类。盒子模型用来描述瞬间泄漏形成的重气云团的运动，平板模型用来描述连续泄漏形成的重气云羽的运动。这两类模型的核心是因空气进入而引起的气云质量增加速率方程。盒子模型盒子模型使用如下假设：I、重气云团为正立的坍塌圆柱体，圆柱体初始高度等于初始半径的一半II、在重气云团内部，温度、密度和危险气体浓度等参数均匀分布

27、。III、重气云团中心的移动速度等于风速。重气扩散的盒子模型示意图如下图所示。一7：订”/777777777不77"7»空气通过界面进入云底>云团向外的役向运动图2重气云团盒子模型坍塌圆柱体的径向蔓延速度由下式确定：1/2，一、Vfdr/dtg(a)/ah(53)式中，Vf为圆柱体的径向蔓延速度(m/s),r为圆柱体半径(m),h为圆柱体高度(m),t为泄漏后时间(s)等式两边同时乘以2r,上式变成下面的形式:2dr2/dt 2 g (54)21/2a)/a"2aa1/2由于假设重气云团和环境之间没有热量交换，重气云团的浮力将守恒，即:g(a)/aVg(a)

28、/aV。(55)将上式代入式(53),积分后得到：r2r。22g(oa)/aV。/1/2t(56)式中，r。为重气云团的初始半径(m),V。为重气云团的初始体积(m3),。为重气云团的初始密度(kg/m3)0由于假设重气云团是圆柱体，初始高度等于初始半径的一半，因此重气云团初始半径的计算公式为：r。(2V。/)1/3(57)随着空气的不断进入，云团的高度和体积也将不断变化。云团体积随时间的变化速率由下式确定：dV/dt(R2)VT(2Rh)VE(58)式中，重气云团体积VR2h,VT和VP分别为空气从顶部和边缘进入重气云团的速率(m/s)。由于重气云团内部危险气体质量守恒，因此，在重气云团扩散

29、过程中，下式存立：22C/COVO/V(horo2)/(hr2)(59)式中，CO和C分别为初始时刻和t时刻重气云团内部危险物质浓度(kg/m3)。任意时刻重气云团的半径按式(56)计算。如果知道任意时刻重气云团高度的计算公式，利用上式就可计算任意时刻重气云团内部危险物质浓度。但这里不准备采用先推导重气云团高度的计算公式，然后计算重气云团体积和危险物质浓度的方法。而是先采用量纲分析法求重气云团的体积和浓度，然后利用上式反推重气云团的高度。无量纲量V/Vo与x/VO/3之间存在如下函数关系：VVo(x/Vo1/3)1.5,xVo1/3(60)式中，x为下风向距离(m)。它与时间、风速之间的关系为

30、：xVt(61)将上式代入式(52)，得到：CCo(x/Vo1/3)1.5,xVo1/3(62)将圆柱形重气云团的体积Vr2h代入式(60),可推导出1/31.521/3hVo(x/Vo)/(r),xVo(63)随着空气的不断进入，重气云团的密度将不断减小，重气坍塌引起的扩散将逐步让位于环境湍流引起的扩散。目前，判断重气坍塌过程终止的准则主要有：I、e准则定义（pa）/a。e准则认为，如果e小于或等于某个临界值（在之间），重气坍塌引起的扩散将让位于环境湍流引起的扩散。下面推导转变发生的位置，令：EgV（pa）/a（64）gV将式（60）代入式（64），得到：EgVo（x/Vo1/3）1.5（6

31、5）从上式求出X,得到：xE2/3Vo1/3（g）2/3（66）由于不考虑云团与环境之间的热交换，云团浮力守恒，故E=E0。代入上式得到转变点对应的下风向距离为：2/31/32/3XfEoVo（gcr）（67）II、Ri准则对于瞬间泄漏，定义Richardson数Ri（g（pa）/a）V1/3/V*2。Ri准则认为，如果Ri小于或等于某个临界值（在110之间），重气坍塌引起的扩散将让位于环境湍流引起的扩散。由于云团内部浮力守恒，因此EoEgV（a）/a（68）对上式进行恒等变换，得到：（a）/aEo/（gV）（69）将式（69）代入Richardson数的定义式，得到：RiEo/（XVo1/3

32、V*2）（70）从上式求出转变点下风向距离X，得到：XfEo/（RicrVo1/3V*2）（71）III、Vf准则定义重气云团径向蔓延速度Vf=dr/dt。Vf准则认为，如果V小于或等于某个临界值，重气坍塌引起的扩散将让位于环境湍流引起的扩散。不同研究人员提出的重气云扩散阶段终止时的临界Vf值相差很大。例如，vanUlden认为，重气云扩散阶段终止的条件是Vf=2V*,Germeles和Drake认为，重气云扩散阶段终止的条件是Vf=V,Cox和Carpenter认为，重气云扩散阶段终止的条件是M=d(ry/dt,h为横风向扩散系数(m),Eidsvik(1980)认为，重气云扩散阶段终止的条

33、件是M=*。这些准则覆盖范围很宽，从V=到V=V(假设摩才S速度2=7/15)。Germeles和Drake提出的准则太严，按照他们提出的准则，即使存在重气云扩散阶段，这个阶段持续的时间也很短。不过，多数研究人员认为，V的临界值具有与M相同的数量别。平板模型平板模型使用了如下假设：I、重气云羽横截面为矩形，横风向半宽为b(m)，垂直方向高度为h(m)。在泄漏源点，云羽半宽为高度的两倍，即：bo=2h。II、重气云羽横截面内，浓度、温度、密度等参数均匀分布。III、重气云羽的轴向蔓延速度等于风速。在重气云羽的扩散过程中，横截面半宽的变化由下式确定：Vdb/dx(gh(a)/a)1/2(72)由于

34、假设重气云羽与环境之间无热量交换，重气云羽的浮力通量在扩散过程中守恒，即：2gVbh(a)/a2gVboho(oa)/a(73)将式(73)代入式(72)，积分后得到：bbo11.5gho(oa)/a1/2x(Vbo)12/3(74)由于重气云羽初始半宽等于初始高度的两倍，重气云羽的初始体积通量为:2Vo'2bohoVbo2V(75)从上式可以求出重气云羽的初始半宽：bo2ho(Vo'/V)1/2(76)随着空气的进入，不仅重气云羽的横风向水平尺寸要增大，重气云羽的高度也要增加。重气云羽高度的变化与下风向距离间的关系由下式确定：dh(we/v)dx(77)式中，We为空气卷吸系

35、数，且假设空气卷吸系数由下式确定：'We3.5V*'/(11.67Ri)(78)式中，Ri为当地Richardson数。上式表明，随着Richardson数的增加，空气卷吸系数减小。Richardson数的定义为：Rig(pa)/ah/V*'2(79)式中，V*'是垂直方向的特征湍流速度(m/s)，由下式确定：V*'1.3(V*/V)(4/9)(db/dt)2V21/2(80)式中，V*是摩擦速度。由于xVt，因此，db/dtVdb/dx，结合式(74)，得到：1/21db/dt(2/3)(Vebo)1.5gho(oa)/a1/2(Vebo)1(81)1

36、1.5gho(oa)/a1/2x(Vebo)11/3求解由式(78)式(81)组成的联立方程组，可以求得任意下风向距离重气云羽的高度。由于重气云羽横截面上危险物质通量守恒，因此2bhVC2bohoVCo(82)上式两边同时除以2bhV,得到重气云羽中危险物质浓度的计算公式：CbohoCo/(bh)(83)式中，C表示重气云羽内危险物质浓度(kg/m3),下标0指初始条件。无论是重气云团扩散，还是重气云羽扩散，一旦满足前面讨论过的转变条件，重力驱动扩散将转变为环境湍流驱动扩散。为了将转变前后两个不同的扩散过程有机衔接起来，需要进行虚源计算。所谓虚源，是指位于转变点上游某处的虚拟泄漏源。虚源计算的

37、目的是确定虚源与转变点之间的距离。进行虚源计算时应遵循下面的原则：在相同的泄漏和扩散条件（相同源强、相同地形、相同气象条件等）下，利用重气云扩散模型对实源泄漏进行扩散分析得到的转变点所在位置危险物质浓度等于利用高斯模型对虚源泄漏进行扩散分析得到的转变点所在位置危险物质浓度。虚源计算时假设转变点的下风向扩散系数x、横风向扩散系数y和垂直方向扩散系数z分别由下面三个公式计算：xy（瞬间泄漏才需要）（84）yb/21/2（85）zh/21/2（86）因此，如果知道扩散系数与下风向距离的关系，就可以计算出虚源与转变点之间的距离。例如，如果y0.1x，那么，虚源与转变点之间的距离xv10b/21/2。库

38、区重大事故后果分析A.2.1库区的主要事故类型根据储存场所的不同以及储存危险品特性的不同，库区主要的事故后果类型如下：（1）仓库中存储TNT等爆炸性物品时，容易发生凝聚相含能材料爆炸，产生非常严重的后果。由于爆炸性物品不得露天堆放，只能存储在仓库中，因此通常发生的是受限空间的爆炸。(2)易燃、易爆的气体(包括液化气体)如液化石油气钢瓶等在仓库中存储时，发生泄漏并在扩散过程中遇到点火源，则很容易发生蒸气云爆炸事故。(3)有毒气体(包括液化气体)如液氯、液氨钢瓶在仓库中存储时，发生泄漏并扩散很容易引起中毒事故。(4)易燃液体如苯、甲醇等瓶装、桶装的化工原料在仓库中存储时，泄漏后很容易引发室内池火灾

39、。(5)易燃固体、自燃物品、遇湿易燃物品等瓶装、桶装、袋装的物品在仓库中存储时，容易发生室内固体火灾事故。(6)易燃液体的桶装的化工原料在堆场中存储时，泄漏后很容易引发大面积的池火灾。(7)易燃固体物品在堆场中存储时，容易发生固体火灾事故。A.2.2库区的主要事故后果模型A.2.2.1凝聚相含能材料爆炸凝聚相含能材料爆炸能产生多种破坏效应，如热辐射、一次破片作用、有毒气体产物的致命效应，但破坏力最强，破坏区域最大的是冲击波的破坏效应，因此，爆炸模型主要考虑冲击波的伤害作用。凝聚相含能材料的爆炸冲击波超压Ap可按下式计算：(87)(88)(89)(90)10.1567Z3(p5)Ps0.137Z

40、30.119Z20.269Z10.019(1p10)PsP/P0zR/()13P0E1.8WQC式中，Z为无量纲距离，Ap为目标处的超压值(Pa),p0为环境压力，R为目标到爆源的水平距离(m,E为爆源总能量(J),W为含能材料的质量(Kg),Q为爆炸物的爆热(J/Kg)。A.2.2.2蒸气云爆炸见A.1.3节蒸气云爆炸事故后果模型。A.2.2.3毒物的泄漏扩散见A.1.3节中毒事故后果模型。A.2.2.4池火灾见A.1.3节池火灾事故后果模型。A.2.2.5固体火灾固体火灾的热辐射参数按点源模型估计。此模型认为火焰射出的能量为燃烧的一部分，并且辐射强度与目标至火源中心距离的平方成反比。q(r

41、)fMcHc4 r2(91)式中，q（r）为目标接收到的辐射强度（W/rri）,f为辐射系数，可取f=,M为燃烧速率（kg/s）,H为燃烧热（J/kg）,r为目标至火源中心的水平距离（rr）0生产场所重大事故后果分析A.3.1生产场所主要事故类型生产场所的事故类型非常复杂，因反应介质、工艺设备与机器、操作条件的不同而不同，常见的危害较大的主要包括以下几类：（1）爆炸物理爆炸化工容器及设备由于设计、制造、腐蚀或低温、材料缺陷、交变载荷的作用，使得器壁的平均应力超过材料的屈服点或强度极限，导致脆性疲劳、疲劳破裂和应力腐蚀破裂发生物理爆炸，也可因安全泄放装置失灵、液化气体充装过量、严重受热膨胀、违章

42、超负荷运行等发生物理爆炸。常见的如水夹汽包、化工容器、液化器气瓶等的爆炸。化学爆炸化工设备和机器内的物质发生极迅速、剧烈的化学反应而产生高温高压可引发瞬间的爆炸现象。一般可分为简单分解爆炸、复杂分解爆炸和爆炸性混合物爆炸。在化工、石油化工生产中发生的化学爆炸绝大部分是爆炸性混合物爆炸。例如由于负压操作、系统用气、水封不严或失效，空气用入到装置中，形成爆炸性混合物，发生化学爆炸；再如硝化反应过程中，由于温度控制不良，很容易引起爆炸。蒸气云爆炸化工设备和机器由于密封装置失效、设备管道腐蚀、磨损或疲劳破裂、断裂以及安装检修不良、操作失误等原因，可燃性气体从化工装置、设备、管道内泄漏或喷射，扩散到周围

43、环境中，达到爆炸极限，若遇到明火或高温很就会发生蒸气云爆炸。(2)燃烧当化工设备和机器由于密封装置失效、设备管道腐蚀、磨损或疲劳破裂、断裂以及安装检修不良、操作失误等原因导致物料泄漏时，对于易燃液体而言，泄漏后形成一定范围的液池，若遇到火源，液池可能被点燃，发生地面池火灾；对于易燃气体而言，泄漏后立即遇到明火或高温，或因高速摩擦产生静电而产生喷射火，也可延时点火产生闪火。含有易燃易爆液化气体的容器在外部火焰的烘烤下可能沸腾液体扩展为蒸气爆炸，产生巨大的火球。(3)中毒在化工、石化生产中，由于设备密封不严、严重腐蚀穿孔、疲劳破裂、磨损、超压引起的设备与管道突然断裂、开错阀门、阀门密封不严、水封失

44、效等原因，很容易造成毒性气体的泄漏，向周围环境扩散，造成人员的中毒事故。主要化工设备和机器的事故类型及损坏尺寸见下表。表6主要化工设备和机器的事故类型及损坏尺寸序号设备类型事故类型损坏尺寸1塔（吸收塔、蒸福塔、萃取塔、干燥塔）物理爆炸、化学爆炸、严重泄漏内部爆炸全部泄漏；孔盖泄漏20%T径；喷嘴断裂1。麻径；管路破裂20M100%管径。2储槽、压缩空气储罐、缓冲罐物理爆炸、化学爆炸、严重泄漏物理爆炸全部泄漏；喷嘴断裂1。娟径；管路破裂20%100%f径。3中间储罐物理爆炸、严重泄漏物理爆炸全部泄漏；接头泄漏20%100%管径。4换热器、冷凝器、冷却器、再沸器严重泄漏并可导致爆炸内部爆炸全部泄漏

45、；管路破裂20%100%管径。5反应釜、合成塔、流化床、物理爆炸、化学爆炸、严重泄漏内部爆炸全部泄漏；孔盖泄漏20%T径；喷嘴断裂100%f径；管路破裂20%100%管径。6管道、法兰、接头挠性连接器、过滤器严重泄漏法兰泄漏20%f径；管道泄漏20M100%管径；接头损坏20M100%T径。7阀门、泄放阀严重泄漏壳体泄漏20M100%f径；盖子、杆损坏20%世漏。8压缩机、泵严重泄漏壳体泄漏20M100%!径；密封盖损坏20%世漏。9风机负压空气进入化学爆炸，泄漏内部爆炸全部泄漏；壳体泄漏20M100%管径。10余热锅炉、加热炉炉管泄漏、炉体爆炸内部爆炸全部泄漏；管路破裂20M100%管径。1

46、1火炬严重泄漏接头泄漏20M100%f径。A.3.2生产场所主要事故后果模型A.3.2.1物理爆炸物理爆炸就是物质状态参数（温度、压力、体积）迅速发生变化，在瞬间放出大量能量并对外做功的现象。物理爆炸的特点是在爆炸现象发生过程中，导致爆炸发生的介质的化学性质不发生变化，发生变化的仅是介质的状态参数。物理爆炸如压力容器破裂时，气体膨胀所释放的能量（即爆破能量）不仅与气体压力和容器的体积有关，而且与介质在容器内的物性相态有关。当压力容器中介质为压缩气体，即以气态形式存在而发生物理爆炸时，其释放的爆破能量为：k10.1013Eg-pV1（-）103（92）k1p式中，Eg为气体的爆破能量（kJ）,p

47、为容器内气体的绝对压力（MPa,V为容器的容积（m3）,k为气体的绝热指数，即气体的定压比热与定容比热之比。当介质全部为液体时，鉴于通常用液体加压时所做的功作为常温液体压力容器爆炸时释放的能量，爆破能量计算模型如下：Ei”VA（93）2式中，Ei为常温液体压力容器爆炸时释放的能量（kJ）,p为液体的绝对压力（Pa）,V为容器的体积（mi）,t为液体在压力p和温度T下的压缩系数（PaI。而液化气体一般在容器内以气液两态存在。当容器破裂发生爆炸时，除了气体的急剧膨胀做功外，还有过热液体激烈的蒸发过程。在大多数情况下，这类容器内的饱和液体占有容器介质重量的绝大部分，它的爆破能量比饱和气体大得多，一般

48、计算时不考虑气体膨胀做的功。过热状态下液体在容器破裂时释放出爆破能量可按下式计算:E (H1H 2) (S1 S2)T1W94)式中,E为过热状态液体的爆破能量（kJ）,H为爆炸前液化液体的始（kJ/kg）,H2为在大气压力下饱和液体的焓（kJ/kg），S1为爆炸前饱和液体的熵（kJ/（kg?C）,&为在大气压力下饱和液体的嫡（kJ/（kg?C）,Ti为介质在大气压力下的沸点（C）,W为饱和液体的质量（kg）。压力容器爆破时，爆破能量向外释放以冲击波能量、碎片能量和容器残余变形能量三种形式表现出来。根据介绍，后二者所消耗的能量只占总爆破能量的3%15%也就是说大部分能量产生空气冲击波。

49、A.3.2.2化学爆炸对于简单分解爆炸、复杂分解爆炸产生的爆炸可用凝聚相含能材料爆炸模型计算见A.2.3节；对于爆炸性混合物气体爆炸，可用蒸气云爆炸事故后果模型，见节。A.3.2.3池火灾见A.1.3节池火灾事故后果模型。A.3.2.4沸腾液体扩展蒸气云爆炸见A.1.3节沸腾液体扩展蒸气云爆炸事故后果模型。A.3.2.5蒸气云爆炸见A.1.3节蒸气云爆炸事故后果模型。A.3.2.6喷射火见A.1.3节喷射火事故后果模型。A.3.2.7闪火模型闪火是可燃性气体或蒸气泄漏到空气中，与之混合后被点燃而发生的一种非爆炸性的燃烧过程。闪火的主要危害来自热辐射和火焰直接接触。可燃物云团的大小决定了可能造成直接火焰接触危害的面积，而云团的大小则部分取决于扩散和泄漏条件。在闪火模型中，假定闪火是一个火焰以恒定速度传播的过程。闪火模型的几何关系如图3所小。图3闪火模型火焰高度火焰可视高度可用近似的半经验公式计算：223.1/3H20d（s/gd）（0/a）2（wr/（1w）（95）式中，H为火焰可视高度（m）,d为云团厚度（m）,s为燃烧速度（m/s）,g为重力加速度（m/s2）,pa为空气密度（kg/m3）,p0为燃一气混合物的密度（kg/m3）,r为理想配比下空气与燃料的质量比，w=小-小st/a（1-Tst）（当小st）