《第五章重大事故后果分析》由会员分享，可在线阅读，更多相关《第五章重大事故后果分析（143页珍藏版）》请在金锄头文库上搜索。
1、第六章第六章 重大事故后果分析重大事故后果分析 重大事故：重大事故：是指重大危险源在运行中突然发生重大泄是指重大危险源在运行中突然发生重大泄漏、火灾或爆炸，其中涉及一种或多种有害物质，并漏、火灾或爆炸，其中涉及一种或多种有害物质，并给现场人员、公众或环境造成即刻的或延迟的严重危给现场人员、公众或环境造成即刻的或延迟的严重危害的事件害的事件 重大事故后果分析重大事故后果分析: :是重大危险源评价和管理的一个重是重大危险源评价和管理的一个重要方面，目的是定量描述一个可能发生的事故将造成要方面，目的是定量描述一个可能发生的事故将造成的人员伤亡、财产损失和环境污染情况。根据分析结的人员伤亡、财产损失和
2、环境污染情况。根据分析结果决策者可以采取适当措施，如设置报警系统、压力果决策者可以采取适当措施，如设置报警系统、压力释放系统、防火系统以及编制应急响应程序等，以减释放系统、防火系统以及编制应急响应程序等，以减少事故发生的可能性或降低事故的危害程度少事故发生的可能性或降低事故的危害程度 数学模型是事故后果定量分析的基础。这些模型通常数学模型是事故后果定量分析的基础。这些模型通常是对假想的事故场景在一系列理想假设的前提下，依据是对假想的事故场景在一系列理想假设的前提下，依据一定的物理化学原理建立的灰箱模型，模型的参数通常一定的物理化学原理建立的灰箱模型，模型的参数通常是由实验得到的。还有一些是纯经
3、验的黑箱模型。依据是由实验得到的。还有一些是纯经验的黑箱模型。依据不同的假设和原理，相同的事故场景可以建立不同的模不同的假设和原理，相同的事故场景可以建立不同的模型描述；同时，由于依据不同的实验数据，有些相同的型描述；同时，由于依据不同的实验数据，有些相同的模型其参数却有所不同甚至相差远。显然，采用不同的模型其参数却有所不同甚至相差远。显然，采用不同的模型对同一事故的后果分析结果会有所不同模型对同一事故的后果分析结果会有所不同 当然，每一模型还有其适用范围，因此在进行后果分当然，每一模型还有其适用范围，因此在进行后果分析时，考虑模型的适用范围以选择合适的模型是非常重析时，考虑模型的适用范围以选
4、择合适的模型是非常重要的。另外，在没有可靠依据选择参数值时，采用保守要的。另外，在没有可靠依据选择参数值时，采用保守的估计或考虑最坏后果也是可以接受的的估计或考虑最坏后果也是可以接受的一、后果分析一般程序一、后果分析一般程序1 1后果分析程序后果分析程序重大事故后果分析主要包括以下步骤重大事故后果分析主要包括以下步骤 (1)(1)划分独立功能单元划分独立功能单元 (2)(2)计算单元中有害物质存量计算单元中有害物质存量 根据工艺流程和设备参数计算单元中有害物质的存量，根据工艺流程和设备参数计算单元中有害物质的存量，并记录物质的种类、相态、温度、压力、体积或质量并记录物质的种类、相态、温度、压力
5、、体积或质量(3)(3)找出设备的典型故障找出设备的典型故障(4)(4)计算泄漏量计算泄漏量(5)(5)计算后果计算后果 分析泄漏后可能造成的火灾、爆炸等后果，选择合分析泄漏后可能造成的火灾、爆炸等后果，选择合适的模型计算事故对生产现场内或现场外的影响适的模型计算事故对生产现场内或现场外的影响 将计算结果整理成表格，并在单元平面图上划出将计算结果整理成表格，并在单元平面图上划出影响范围影响范围 (6)(6)整理结果整理结果2 2后果分析需要的参数后果分析需要的参数 (1)(1)有害物质的参数有害物质的参数 有害物质的相态、最大质量或体积、温度、压有害物质的相态、最大质量或体积、温度、压力、密度
6、，热力学性质如沸点、蒸发热、燃烧热、力、密度，热力学性质如沸点、蒸发热、燃烧热、比热容等，有害与毒性参数等比热容等，有害与毒性参数等(2)(2)设备的参数设备的参数 工艺流程、设备类型、设备的可能故障与泄漏位工艺流程、设备类型、设备的可能故障与泄漏位置、泄漏口形状尺寸等置、泄漏口形状尺寸等 (3)(3)现场情况与气象情况现场情况与气象情况 设备布置、人员分布、资金密度、设备地理位置、设备布置、人员分布、资金密度、设备地理位置、堤坝高度面积、常年主导风向、平均风速、大气稳定情堤坝高度面积、常年主导风向、平均风速、大气稳定情况、日照情况、地形情况、地面粗糙度、建筑、树木高况、日照情况、地形情况、地
7、面粗糙度、建筑、树木高度等度等 3 3后果分析模式选择后果分析模式选择 重大事故后果分析关心的是易燃、易爆或有毒的重大事故后果分析关心的是易燃、易爆或有毒的气体和液体，这些物质的泄漏不仅有害而且难以控制。气体和液体，这些物质的泄漏不仅有害而且难以控制。一种泄漏可能带来不同的后果，进行后果分析就需要一种泄漏可能带来不同的后果，进行后果分析就需要对每一种可能后果进行计算。采用系统分析的方法可对每一种可能后果进行计算。采用系统分析的方法可以避免对可能的后果造成遗漏以避免对可能的后果造成遗漏 例如：易燃气体泄漏着火时才有危险性，如果泄漏例如：易燃气体泄漏着火时才有危险性，如果泄漏时立即被点燃，则不形成
8、大的蒸汽云团。根据泄漏性时立即被点燃，则不形成大的蒸汽云团。根据泄漏性质可形成喷射火或火球，它能迅速危及事故现场，但质可形成喷射火或火球，它能迅速危及事故现场，但很少能影响到厂区以外。如果泄漏后延迟点燃，则气很少能影响到厂区以外。如果泄漏后延迟点燃，则气体形成云团飘向下风向，点燃后可能造成闪火或爆炸，体形成云团飘向下风向，点燃后可能造成闪火或爆炸，能引起大面积损害能引起大面积损害 计算燃烧和爆炸的热量或压力，不仅仅用于评价人计算燃烧和爆炸的热量或压力，不仅仅用于评价人员和设备的损失情况。燃烧和爆炸还会波及相邻的危险员和设备的损失情况。燃烧和爆炸还会波及相邻的危险源，产生多米诺效应，因此也要对相
9、邻危险源进行泄漏源，产生多米诺效应，因此也要对相邻危险源进行泄漏后果分析后果分析 气体泄漏分析的一个重要方面是计算蒸汽云的密气体泄漏分析的一个重要方面是计算蒸汽云的密度，密度高于空气或低于空气，对其扩散有较大的影度，密度高于空气或低于空气，对其扩散有较大的影响，应该采用不同的扩散模式响，应该采用不同的扩散模式 毒性气体的泄漏扩散分析较简单，因为不需要考虑起火。毒性气体的泄漏扩散分析较简单，因为不需要考虑起火。主要的问题是根据蒸汽云密度选择适当的扩散模式主要的问题是根据蒸汽云密度选择适当的扩散模式气气体体和和两两相相泄泄漏漏事事故故框框图图 液体泄漏着火一般影响的面积较小，但挥发性液体液体泄漏着
10、火一般影响的面积较小，但挥发性液体的蒸汽应按照气体事故进一步分析的蒸汽应按照气体事故进一步分析 常压液体泄漏后在地面形成液池，池内液体由于表常压液体泄漏后在地面形成液池，池内液体由于表面风的作用而缓慢蒸发。如果点燃则形成池火，火焰的面风的作用而缓慢蒸发。如果点燃则形成池火，火焰的热辐射会危及现场人员和设备。加压液化气体泄漏时将热辐射会危及现场人员和设备。加压液化气体泄漏时将发生闪蒸，剩下的液体形成液池。闪蒸的气体应按气体发生闪蒸，剩下的液体形成液池。闪蒸的气体应按气体事故进一步分析事故进一步分析 冷冻液体泄漏也形成液池，液体吸收周围热量蒸发，蒸冷冻液体泄漏也形成液池，液体吸收周围热量蒸发，蒸发
11、速度虽然比闪蒸慢，但一般比常压液体快发速度虽然比闪蒸慢，但一般比常压液体快 沸腾液体扩展蒸汽爆炸是一种比较特殊但后果极其沸腾液体扩展蒸汽爆炸是一种比较特殊但后果极其严重的事故。通常是装液化气体的容器受到外界火焰加严重的事故。通常是装液化气体的容器受到外界火焰加热，一方面使容器内压力升高，同时使容器强度下降。热，一方面使容器内压力升高，同时使容器强度下降。一旦容器突然破裂，大量沸腾液体立即被点燃，形成巨一旦容器突然破裂，大量沸腾液体立即被点燃，形成巨大火球，影响非常严重大火球，影响非常严重 液液体体泄泄漏漏事事故故框框图图二、泄漏二、泄漏1.1.常见泄漏源常见泄漏源 泄漏可能是设备损坏、失灵造成
12、的，也可能是错误操泄漏可能是设备损坏、失灵造成的，也可能是错误操作引起有害物质排放，甚至安全阀的正常或不正常动作。作引起有害物质排放，甚至安全阀的正常或不正常动作。分析时首先找出引起泄漏的设备及可能的泄漏规模，然后分析时首先找出引起泄漏的设备及可能的泄漏规模，然后按适当模型计算泄漏量按适当模型计算泄漏量 根据泄漏面积的大小和泄漏持续时间的长短，将泄根据泄漏面积的大小和泄漏持续时间的长短，将泄漏源分为两类漏源分为两类: :(1)(1)小孔泄漏小孔泄漏，通常为物料经较小的孔洞，通常为物料经较小的孔洞长时间持续泄漏，如反应器、储罐、管道上出现小孔，长时间持续泄漏，如反应器、储罐、管道上出现小孔，或者
13、是阀门、法兰、机泵、转动设备等处密封失效或者是阀门、法兰、机泵、转动设备等处密封失效(2)(2)大面积泄漏大面积泄漏，是指经较大孔洞在很短时间内泄漏出，是指经较大孔洞在很短时间内泄漏出大量物料，如大管径管线断裂、爆破片爆裂、反应器因大量物料，如大管径管线断裂、爆破片爆裂、反应器因超压爆炸等瞬间泄漏出大量物料超压爆炸等瞬间泄漏出大量物料2.2.泄漏量计算泄漏量计算2.12.1液体经小孔泄漏的源模式液体经小孔泄漏的源模式 系统与外界无热交换，流体流动的不同能量形式遵系统与外界无热交换，流体流动的不同能量形式遵守如下的机械能守恒方程守如下的机械能守恒方程式中，式中，p p压力，压力，PaPa，习惯上
14、将压强也称为压力，习惯上将压强也称为压力 流体密度，流体密度，kgkgm m-3-3 动能校正因子，无因次动能校正因子，无因次, ,从工程计算角度出发，值从工程计算角度出发，值近似取为近似取为1U U流体平均速度，流体平均速度，m ms s-1-1，简称流速，简称流速（1 1）g重力加速度加重力加速度加ms-2z高度，高度，m，以基准面为起始，以基准面为起始F阻力损失，阻力损失，Jkg-1ws轴功轴功,J；稳定流动过程中开口系统所作的功；稳定流动过程中开口系统所作的功m质量，质量，kg对于不可压缩流体，密度恒为常数，有对于不可压缩流体，密度恒为常数，有泄漏过程暂不考虑轴功，泄漏过程暂不考虑轴功
15、，WsWs0 0，则（，则（1 1）式化简为）式化简为(2)(2)(3)(3) 工艺单元中的液体在稳定的压力作用下，经薄壁小工艺单元中的液体在稳定的压力作用下，经薄壁小孔泄漏，容器内的压力为孔泄漏，容器内的压力为p p1 1、孔直径为、孔直径为d d，面积为，面积为A A，容，容器外为大气压力。容器内液体流速可以忽略，不考虑摩器外为大气压力。容器内液体流速可以忽略，不考虑摩擦损失和液位变化，利用擦损失和液位变化，利用(3)(3)，可得，可得(4)(4)(5)(5)(6)(6)液体在稳定压力液体在稳定压力下经薄壁小孔泄下经薄壁小孔泄漏漏 考虑到因惯性引起的截面收缩以及摩擦引起的速度考虑到因惯性引
16、起的截面收缩以及摩擦引起的速度减低，引入孔流系数减低，引入孔流系数C C0 0 为实际流量与理想流量的比值，为实际流量与理想流量的比值，则经小孔泄漏的实际质量流量为则经小孔泄漏的实际质量流量为(7)(7)2.22.2储罐中液体经小孔泄漏的源模式储罐中液体经小孔泄漏的源模式 液体储罐，距液位高度液体储罐，距液位高度Z Z0 0处有一小孔壁，在静处有一小孔壁，在静压能和势能的作用下，储罐中的液体经小孔向外泄压能和势能的作用下，储罐中的液体经小孔向外泄漏。泄漏过程可由机械能守恒方程描述，储罐内的漏。泄漏过程可由机械能守恒方程描述，储罐内的液体流速可以忽略。储罐内的液体压力为液体流速可以忽略。储罐内的
17、液体压力为P Pg g，外部，外部为大气压力为大气压力( (表压表压P P0 0）(8)(8)将（将（8 8）代入（）代入（3 3）求出泄漏速度）求出泄漏速度（9 9）若小孔截面积为若小孔截面积为A A，则质量流量，则质量流量Q Q为为(10)(10)由（由（9 9）（）（1010）看出，随着泄漏过程的延续，储罐内）看出，随着泄漏过程的延续，储罐内液位高度不断下降，泄漏速度和质量流量也均随之降液位高度不断下降，泄漏速度和质量流量也均随之降低。如果储罐通过呼吸阀或弯管与大气连通，则内外低。如果储罐通过呼吸阀或弯管与大气连通，则内外压力差压力差 p p为为0 0。（。（10)10)简化为简化为（1
18、111）若储罐的横截面积为若储罐的横截面积为A A0 0，则可经小孔泄漏的最大液体总量为，则可经小孔泄漏的最大液体总量为(12)(12)取一微元时间内液体的泄漏量取一微元时间内液体的泄漏量（1313）储罐内液体质量的变化速率即为泄漏质量流量储罐内液体质量的变化速率即为泄漏质量流量（1414）将式将式(11)(11)、式、式(13)(13)代入式代入式(14)(14)，得到，得到(15)(15) 由边界条件由边界条件t t0 0，z zz z0 0；t=t=t,zt,z= =z,z,对上式进行分对上式进行分离变量积分，有离变量积分，有(16)(16) 当液体泄漏至泄漏点液位后，泄漏停止点当液体泄
19、漏至泄漏点液位后，泄漏停止点z=0z=0，根据，根据上式可得到总的泄漏时间上式可得到总的泄漏时间（1717）将式将式(16)(16)代入到式代入到式(11)(11)，可以得到随时间变化的质量流量，可以得到随时间变化的质量流量(18)(18) 如果储罐内盛装的是易燃液体，为防止可燃蒸气如果储罐内盛装的是易燃液体，为防止可燃蒸气大量泄漏至空气中，或空气大量进入储罐内的气相空大量泄漏至空气中，或空气大量进入储罐内的气相空间形成爆炸性混合物，通常情况下会采取通氮气保护间形成爆炸性混合物，通常情况下会采取通氮气保护的措施。液体表压为的措施。液体表压为PgPg，内外压差即为，内外压差即为PgPg根据式根据
20、式(10)(10)、式式(12)(12)、式、式(13)(13)、式、式(14)(14)可同理得到可同理得到(19)(19)(20)(20)将式将式(20)(20)代人式代人式(10)(10)得到任意时刻的质量流量得到任意时刻的质量流量Q Q(21)(21)例：某例：某盛装丙酮液体的储罐，上部装设有呼吸阀与盛装丙酮液体的储罐，上部装设有呼吸阀与大气连通。在其下部有一泄漏孔，直径为大气连通。在其下部有一泄漏孔，直径为4cm4cm。已知。已知丙酮的密度为丙酮的密度为800kg800kgm m3 3。求：求：(1)(1)最大泄漏量最大泄漏量 (2)(2)泄漏质量流量随时间变化的表达式泄漏质量流量随时
21、间变化的表达式 (3)(3)最大泄漏时间最大泄漏时间 (4)(4)泄漏量随时间变化的表达式。泄漏量随时间变化的表达式。 解：解：(1)(1)最大泄漏量即为泄漏点最大泄漏量即为泄漏点液位以上的所有液体量液位以上的所有液体量 (2)(2)泄漏质量流量随时间变化的泄漏质量流量随时间变化的表达式，表达式，CoCo取值为取值为1 1，则，则(3)(3)令泄漏质量流量时间表达式的左侧为令泄漏质量流量时间表达式的左侧为0 0，即得最大泄漏时间，即得最大泄漏时间(4)(4)任一时间内总的泄漏量为泄漏质量流量对时间的积分任一时间内总的泄漏量为泄漏质量流量对时间的积分若给定任意泄漏时间，即可得到已经泄漏的液体总量
22、若给定任意泄漏时间，即可得到已经泄漏的液体总量2.32.3气体或蒸气经小孔泄漏的源模式气体或蒸气经小孔泄漏的源模式 在工程上，通常将气体或蒸气近似为理想气体，其在工程上，通常将气体或蒸气近似为理想气体，其压力、密度、温度等参数遵循理想气体状态方程压力、密度、温度等参数遵循理想气体状态方程(2.3.1)(2.3.1) 气体或蒸气在小孔内绝热流动，压力密度关系可用气体或蒸气在小孔内绝热流动，压力密度关系可用绝热方程描述绝热方程描述（2.3.22.3.2） 式中式中 绝热指数，是等压热容与等容热容的比值，绝热指数，是等压热容与等容热容的比值， =C=Cp p/ /C Cv v 气体或蒸气经小孔泄漏的
23、过程。轴功为气体或蒸气经小孔泄漏的过程。轴功为0 0，忽略，忽略势能变化则势能变化则 机械能守恒方程机械能守恒方程(1)(1)简化为简化为(2.3.3)(2.3.3)根据根据(8)(8)式定义孔流系数式定义孔流系数 （2.3.42.3.4）将式将式(2.3.4(2.3.4代人式代人式(2.3.3)(2.3.3)，忽略气体或蒸气的初，忽略气体或蒸气的初始动能，得到始动能，得到(2.3.5)(2.3.5)（2.3.62.3.6）由（由（2.3.22.3.2）得到）得到将（将（2.3.62.3.6）代入（）代入（2.3.52.3.5）并积分可得到）并积分可得到（2.3.72.3.7）由式由式(2.3
24、.6)(2.3.6)、式、式(2.3.7)(2.3.7)得到泄漏质量流量得到泄漏质量流量(2.3.8)(2.3.8)根据理想气体状态方程，有根据理想气体状态方程，有（2.3.92.3.9）（2.3.102.3.10）将（将（2.3.92.3.9）代入（）代入（2.3.82.3.8）可得）可得 从安全工作的角度考虑，我们关心的是经小孔泄漏从安全工作的角度考虑，我们关心的是经小孔泄漏的气体或蒸气的最大流量。式的气体或蒸气的最大流量。式(2.3.10)(2.3.10)表明泄漏质量流表明泄漏质量流量由前后压力的比值所决定。若以压力比量由前后压力的比值所决定。若以压力比P/PP/P0 0为横坐标，为横坐
25、标，以流量以流量Q Q为纵坐标，根据式为纵坐标，根据式(2.3.10)(2.3.10)可得到如图中的可得到如图中的0bc0bc曲线，当曲线，当P Pp p0 01 1时，小孔前后的压力相等，时，小孔前后的压力相等，Q=0Q=0；当；当P/PP/P0 00 0时，气体或蒸气流向绝对真空，时，气体或蒸气流向绝对真空，P P0 0，故，故Q Q0 0流流量曲线存在最大值，令量曲线存在最大值，令dQdQd(P/Pd(P/P0 0) )0 0，可求得极值条，可求得极值条件件(2.3.11)(2.3.11)p pc c称为临界压力称为临界压力 将此极值条件代入式将此极值条件代入式(2.3.7)(2.3.7
26、)、式、式(2.3.10)(2.3.10)可得到最可得到最大流速和最大流量大流速和最大流量（2.3.122.3.12）（2.3.132.3.13） 由曲线可以看到，当由曲线可以看到，当P PP PC C时，气体或蒸气流速低时，气体或蒸气流速低于音速，如图中于音速，如图中bcbc段曲线所示。当段曲线所示。当P=PP=PC C时，气体或蒸气时，气体或蒸气的泄漏速度刚好可能达到的最大流速如式的泄漏速度刚好可能达到的最大流速如式(2.3.12)(2.3.12)所示，所示，实际上就是气体或蒸气中的音速。当实际上就是气体或蒸气中的音速。当P PP PC C时，气体或时，气体或蒸气似乎可以充分降压、膨胀、加
27、速，但是根据气体流蒸气似乎可以充分降压、膨胀、加速，但是根据气体流动力学的原理，泄漏速度不可能超过音速，这时其泄漏动力学的原理，泄漏速度不可能超过音速，这时其泄漏速度和质量流量与速度和质量流量与P PP PC C时相同，因此在图中以时相同，因此在图中以abab线表线表示。在化工生产中发生的气体或蒸气泄漏，很多属于最示。在化工生产中发生的气体或蒸气泄漏，很多属于最后一种情况后一种情况 例：在某生产厂有一空气柜，因外力撞击，在空气柜例：在某生产厂有一空气柜，因外力撞击，在空气柜一侧出现一小孔。小孔面积为一侧出现一小孔。小孔面积为19.6cm19.6cm2 2，空气柜中的空，空气柜中的空气经此小孔泄
28、漏入大气。已知空气柜中的压力为气经此小孔泄漏入大气。已知空气柜中的压力为2.5x102.5x105 5PaPa温度温度ToTo为为330K330K，大气压力为，大气压力为10105 5PaPa绝热指绝热指数数 1.401.40。求空气泄漏的最大质量流量。求空气泄漏的最大质量流量解：先根据式解：先根据式(2.3.11)(2.3.11)判断空气泄漏的临界压力判断空气泄漏的临界压力 大气压力为大气压力为10105 5PaPa，小于临界压力，则空气泄漏，小于临界压力，则空气泄漏的最大质量流量可按式的最大质量流量可按式(2.3.13)(2.3.13)计算计算C C0 01 1，则，则1.09kg/s1.
29、09kg/s1 1若若C C0 0值取为值取为0.610.61，则空气泄漏的最大质量流量为，则空气泄漏的最大质量流量为2.4 2.4 两相泄漏两相泄漏 如果容器中的过热液体泄漏前通过较长的管道就如果容器中的过热液体泄漏前通过较长的管道就会产生两相泄漏，假设系统中出口和上游临界压力会产生两相泄漏，假设系统中出口和上游临界压力比为比为0.550.55，则，则P PC C=0.55P=0.55P泄漏两相中蒸发液体分数泄漏两相中蒸发液体分数F FV V为：为：（2.4.12.4.1） C CP P 为液体恒压比热容，为液体恒压比热容，T TC C为临界压力下的沸点为临界压力下的沸点,H,HV V常压沸
30、点下的汽化热常压沸点下的汽化热 两相流中气相和液相混合物的平均密度为两相流中气相和液相混合物的平均密度为 （设气（设气相密度为相密度为 g g, ,液相密度为液相密度为 l l）（1 1） 两相流排放泄漏流量为：两相流排放泄漏流量为：(2)(2)C C0 0为两相流泄漏系数，一般取为两相流泄漏系数，一般取0.80.8三、泄漏后的扩散三、泄漏后的扩散 泄漏出的介质立即表现出不同的行为，这与其泄漏出的介质立即表现出不同的行为，这与其储存的状态和泄漏情况有关。沸点以下的液体泄漏，储存的状态和泄漏情况有关。沸点以下的液体泄漏，如果挥发性较低，则蒸汽对现场人员有伤害，但一般如果挥发性较低，则蒸汽对现场人
31、员有伤害，但一般不会影响到厂外。如果挥发性高，则蒸汽会在大气中不会影响到厂外。如果挥发性高，则蒸汽会在大气中扩散扩散 对于过热液体泄漏，介质喷出后存在一个绝热膨对于过热液体泄漏，介质喷出后存在一个绝热膨胀过程。液体的泄漏还可以产生池火，气体泄漏则存胀过程。液体的泄漏还可以产生池火，气体泄漏则存在喷射扩散。如果泄漏初期没有被点燃，最终都将发在喷射扩散。如果泄漏初期没有被点燃，最终都将发展成扩散的蒸汽云。可以用大气中的蒸汽扩散描述，展成扩散的蒸汽云。可以用大气中的蒸汽扩散描述，进一步还可以分析火灾、爆炸以及毒害后果进一步还可以分析火灾、爆炸以及毒害后果1.1.液体的扩散与蒸发液体的扩散与蒸发 (1
32、)(1)液体扩散液体扩散 分析液体扩散的分析液体扩散的关键是找到液体扩散形成液池的面积关键是找到液体扩散形成液池的面积，因为随后的蒸发过程和燃烧过程的计算都直接依赖液池面因为随后的蒸发过程和燃烧过程的计算都直接依赖液池面积。如果储罐区建有防火堤，则泄漏液体只能达到防火堤，积。如果储罐区建有防火堤，则泄漏液体只能达到防火堤，液池面积就是防火堤所围面积。没有防火堤则液体流向低液池面积就是防火堤所围面积。没有防火堤则液体流向低洼处，液池面积也可以估计。对于土地较平整的情况，液洼处，液池面积也可以估计。对于土地较平整的情况，液体将扩散至达到最小液体厚度为止体将扩散至达到最小液体厚度为止 液池面积的确定
33、是事故后果分析中最困难也是液池面积的确定是事故后果分析中最困难也是最容易引起误差的地方。可以简单假定扩散在平整、最容易引起误差的地方。可以简单假定扩散在平整、光滑的平面进行，而且没有渗漏损失，扩散期间也光滑的平面进行，而且没有渗漏损失，扩散期间也不考虑挥发不考虑挥发 对于瞬时泄漏：对于瞬时泄漏： 其中其中对于连续泄漏：对于连续泄漏：r r为液池半径为液池半径 对于瞬时泄漏，如果泄漏的液体已经充分扩展，对于瞬时泄漏，如果泄漏的液体已经充分扩展，假定液体无挥发，地面无渗漏，若已知液层的最小假定液体无挥发，地面无渗漏，若已知液层的最小厚度，可求液池面积厚度，可求液池面积V V泄漏液体体积，泄漏液体体
34、积，S S液池面积，液池面积，H Hminmin液层最小厚度液层最小厚度 池火计算一般以圆池为模型，其他形状液池化为等池火计算一般以圆池为模型，其他形状液池化为等面积圆，直径为面积圆，直径为（2 2）液体蒸发）液体蒸发 低温液体或闪蒸后剩余的液体，主要吸收地面低温液体或闪蒸后剩余的液体，主要吸收地面热量进行蒸发，蒸发速率热量进行蒸发，蒸发速率 m m蒸发速率，蒸发速率，kg/(mkg/(m2 2. .s), s), s s表面热导率，表面热导率，w/(m.Kw/(m.K), a), as s热扩散系数，热扩散系数，m m2 2/s, T/s, Ta a环境温度，环境温度，T Tb b液体沸点，
35、液体沸点， H HV V蒸发热，蒸发热，J/kg, tJ/kg, t蒸发时间蒸发时间（1 1） 根据扩散通量正比于液池表面饱和蒸气浓度与根据扩散通量正比于液池表面饱和蒸气浓度与其在大气中的本底浓度之差，结合理想气体状态方其在大气中的本底浓度之差，结合理想气体状态方程，可得液体蒸发速度公式程，可得液体蒸发速度公式（2 2）式中式中 k k扩散传质系数；扩散传质系数； p ps s液体饱和蒸汽液体饱和蒸汽压，压，PaPa；M M摩尔质量，摩尔质量，kgkgmolmol传质系数可以按下式计算传质系数可以按下式计算式中式中 u ul0ml0m高处风速，高处风速，m/hm/h；scscSchmidtSc
36、hmidt数，数， / / D D； 空气黏度，空气黏度，kgkg( (m mh h) )； 空气密度，空气密度，kgkgm m3 3； D D蒸发液体的扩散系数，蒸发液体的扩散系数，m m2 2/h/h传质系数也可简单按下式确定传质系数也可简单按下式确定2.2.喷射扩散喷射扩散 气体在压力从裂口喷出，一般温度会低于环境气体在压力从裂口喷出，一般温度会低于环境温度，如果气体流动是阻塞的，压力会高于环境压力。温度，如果气体流动是阻塞的，压力会高于环境压力。首先需要把喷射流出转化为具有环境条件的等效流动首先需要把喷射流出转化为具有环境条件的等效流动等效裂口直径与实际裂口直径的关系为：等效裂口直径与
37、实际裂口直径的关系为： D Deqeq等效裂口直径；等效裂口直径；D D计算泄漏流量用裂口直径；计算泄漏流量用裂口直径； 0 0气体刚流出时与环境条件下空气相对密度；气体刚流出时与环境条件下空气相对密度； 气气体在环境条件下与同条件下空气的相对密度体在环境条件下与同条件下空气的相对密度喷射轴线上距喷射孔喷射轴线上距喷射孔x x处的浓度处的浓度b1,b2b1,b2为分布系数为分布系数b b1 1=50.5+48.2=50.5+48.2 -9.95-9.95 2 2b b2 2=23.0+41.0=23.0+41.0 垂直于喷射轴的水平面上的浓度分布为垂直于喷射轴的水平面上的浓度分布为C Cx,y
38、x,y 喷射水平面上喷射水平面上(x(x，y)y)处的浓度处的浓度沿轴的喷射速度分布沿轴的喷射速度分布U Ux x喷射轴上距喷射孔喷射轴上距喷射孔x x处的喷射速度；处的喷射速度；u u0 0实实际泄漏气流速度，由泄漏流量和裂口面积计算际泄漏气流速度，由泄漏流量和裂口面积计算 首先计算出喷射速度等于风速的相应位置首先计算出喷射速度等于风速的相应位置x x，再，再计算与计算与x x相应的泄漏气体浓度相应的泄漏气体浓度3.3.绝热扩散绝热扩散 闪蒸的液体或压缩气体瞬时释放后有一个快速膨胀闪蒸的液体或压缩气体瞬时释放后有一个快速膨胀过程，由于该过程非常快，以至于气团和环境之间没有过程，由于该过程非常
39、快，以至于气团和环境之间没有时间讲行热交换，此膨胀过程可以按绝热过程处理时间讲行热交换，此膨胀过程可以按绝热过程处理 该模型假定气云是呈包含两个区间的半球状，内层该模型假定气云是呈包含两个区间的半球状，内层“核核”具有均匀的浓度，包含具有均匀的浓度，包含5050的泄漏质量，外层浓的泄漏质量，外层浓度呈高斯分布，具有另外度呈高斯分布，具有另外5050的泄漏量的泄漏量 这种双层云团扩散假定分两步：（这种双层云团扩散假定分两步：（1 1）气体或气溶）气体或气溶胶膨胀到压力降至大气压，在膨胀过程中气团获得动胶膨胀到压力降至大气压，在膨胀过程中气团获得动能，称为膨胀能（能，称为膨胀能（2 2）在膨胀能作
40、用下气团进一步扩）在膨胀能作用下气团进一步扩张，推动空气紊流混合进入气团。假设第二阶段持续张，推动空气紊流混合进入气团。假设第二阶段持续到核的扩张速度降到某给定值时结束到核的扩张速度降到某给定值时结束 第一步膨胀到大气压，膨胀期间气体反抗大气压做第一步膨胀到大气压，膨胀期间气体反抗大气压做功，部分内能用于增加物质的动能。如果假定增加的功，部分内能用于增加物质的动能。如果假定增加的动能由动能由(P(P1 1一一Pa)dvPa)dv给出，则初始的膨胀到大气压的过给出，则初始的膨胀到大气压的过程可以按可逆绝热过程处理。膨胀能是始态能量和末程可以按可逆绝热过程处理。膨胀能是始态能量和末态能量的差，减去
41、对大气所做的功。按照这种理想化态能量的差，减去对大气所做的功。按照这种理想化的方法，第一步是等熵的的方法，第一步是等熵的 对于气体释放，能量由可逆绝热膨胀决定，即由对于气体释放，能量由可逆绝热膨胀决定，即由P P1 1、T T1 1( (具有内能具有内能U U1 1、体积、体积v v1 1) )膨胀到膨胀到P P2 2、T T2 2( (具有内能具有内能U U2 2、体、体积积v v2 2) )、(1)(1)扩散能为扩散能为内能的变化为内能的变化为(2)(2)对于液体释放，闪蒸分数按等熵过程计算对于液体释放，闪蒸分数按等熵过程计算(3)(3)(4)(4)膨胀能是始、终态内能的变化减去对大气所做
42、的功膨胀能是始、终态内能的变化减去对大气所做的功(5)(5)式中式中 E E膨胀能，膨胀能， U U内能，内能， T T温度，温度， TbTb沸点，沸点， C CV V气体的比热容，气体的比热容，J/(J/(kgkgK K) )， PaPa一一大气压，一一大气压， S S熵，熵，J J( (kgkgK K) )， H Hv v蒸发热，蒸发热，J Jkgkg， H H焓，焓，J Jkgkg 内能、温度、熵和焓的下标内能、温度、熵和焓的下标1 1和和2 2分别表示状态分别表示状态1 1和和状态状态2 2，L L和和v v分别表示液态和气态分别表示液态和气态 第二步空气紊流混合，膨胀的推动产生大范围
43、的紊第二步空气紊流混合，膨胀的推动产生大范围的紊流，这种紊流是气团与空气进一步混合的决定因素流，这种紊流是气团与空气进一步混合的决定因素紊流扩散系数紊流扩散系数(6)(6)内核半径随时间的变化内核半径随时间的变化内核浓度随时间的变化内核浓度随时间的变化(7)(7)(8)(8) 当内核扩张速度当内核扩张速度( (drcdrcdtdt) )降至给定值时第二阶段降至给定值时第二阶段结束。临界速度的选择是任意的，但通常的推荐值结束。临界速度的选择是任意的，但通常的推荐值是是lmlms s。选定此速度再结合扩散能以及内核半径、。选定此速度再结合扩散能以及内核半径、内核浓度与时间的关系，可以得到第二阶段结
44、束时内核浓度与时间的关系，可以得到第二阶段结束时的内核半径和浓度的内核半径和浓度(9)(9)(10)(10)扩散第二阶段结束时，半球形气团的半径按下式计算扩散第二阶段结束时，半球形气团的半径按下式计算气团密度是绝热膨胀后采用其他模型进一步分气团密度是绝热膨胀后采用其他模型进一步分析的重要参数，对于气体释放析的重要参数，对于气体释放, 气团的体积为：气团的体积为：(11(11) )(12)(12) 对于液体释放，如果膨胀后仍有液体存在，则计对于液体释放，如果膨胀后仍有液体存在，则计算就比较复杂。确定气团中是否有液体，需要首先计算就比较复杂。确定气团中是否有液体，需要首先计算气团的浓度算气团的浓度
45、CceCce，如果大干，如果大干1 1则液体存在则液体存在 含有液体的气团计算混合空气的旦时，需要使用含有液体的气团计算混合空气的旦时，需要使用3 3个方程反复计算个方程反复计算 描述最终温度为描述最终温度为T T3 3、最终气体分数为、最终气体分数为F Fv3v3的各组分的的各组分的热平衡热平衡体积平衡体积平衡(13(13) )(14)(14) 第三个方程表明温度和气体分数必须与所要求第三个方程表明温度和气体分数必须与所要求的分压等于的分压等于T T3 3时的饱和蒸汽压一致时的饱和蒸汽压一致(15)(15)空气的质量必须调整至满足以上空气的质量必须调整至满足以上3 3个方程个方程4.4.重气
46、扩散重气扩散 研究重气云扩散的意义研究重气云扩散的意义(1)(1)重气由于受重力作用，重气由于受重力作用，一般是沿地面扩散的，对人的影响更大；一般是沿地面扩散的，对人的影响更大；(2)(2)常见的常见的工业气体，很多密度都大于空气，如液化石油气、氯工业气体，很多密度都大于空气，如液化石油气、氯气、二氧化硫、硫化氢等。（气、二氧化硫、硫化氢等。（3 3）一些液化气体虽然）一些液化气体虽然分子量不一定大于空气的表观分子量，但是泄漏后一分子量不一定大于空气的表观分子量，但是泄漏后一般温度较低，还可能夹带液滴，表观密度仍然大于空般温度较低，还可能夹带液滴，表观密度仍然大于空气，又称为稠密云团，也应该使
47、用重气气，又称为稠密云团，也应该使用重气( (稠密云团稠密云团) )扩扩散模型散模型 重气扩散模型可分为三类：重气扩散模型可分为三类：（1 1）一维模型（一维模型（2 2）二）二维模型（维模型（3 3）三维模型）三维模型 最简单的是一维模型，对于瞬时泄漏，假定气团按最简单的是一维模型，对于瞬时泄漏，假定气团按风速运动，气团内部具有均匀的浓度和密度；对于连续风速运动，气团内部具有均匀的浓度和密度；对于连续泄漏，假定是云羽按风速运动，垂直于风向的截面上具泄漏，假定是云羽按风速运动，垂直于风向的截面上具有均匀的浓度和密度有均匀的浓度和密度 不论是重气泄漏、还是中性或浮性气体泄漏，在建不论是重气泄漏、
48、还是中性或浮性气体泄漏，在建立气体在大气中的扩散模型时，一般做如下的理想化假立气体在大气中的扩散模型时，一般做如下的理想化假设设 瞬时泄漏假定瞬时完成，连续泄漏假定泄漏速瞬时泄漏假定瞬时完成，连续泄漏假定泄漏速 率恒定率恒定 气云在乎整、无障碍物的地面上空扩散气云在乎整、无障碍物的地面上空扩散 气云中不发生化学反应，地面对气云无吸收气云中不发生化学反应，地面对气云无吸收 风向为水平风向，风速和风向不随时间变化风向为水平风向，风速和风向不随时间变化4.1 4.1 瞬时扩散的盒子模型瞬时扩散的盒子模型(1)(1)起始云团形状起始云团形状 盒子模型假定起始云团为高盒子模型假定起始云团为高H H0 0
49、、半径、半径R R0 0的圆柱体。的圆柱体。用于绝热扩散的半球状气团时，可化为等体积的圆柱用于绝热扩散的半球状气团时，可化为等体积的圆柱体，高与半径比通常选体，高与半径比通常选1 1重气云团盒子模型扩散示意图重气云团盒子模型扩散示意图（2)2)重力影响重力影响 在重力作用下，圆柱体高度变小，半径增加，半在重力作用下，圆柱体高度变小，半径增加，半径的变化速度按下式计算径的变化速度按下式计算R R云团半径，云团半径，t t时间，时间，H H云团高度，云团高度，云团密度，云团密度，a a空气的密度，空气的密度，K K常数，这里常数，这里取取1 , g1 , g重力加速度重力加速度(16) (16)
50、(3)(3)空气混入气团空气混入气团 气团随空气运动时，空气从气团顶部和边缘进入，气团随空气运动时，空气从气团顶部和边缘进入，使其浓度和密度逐渐下降。空气从边缘进入的速度正使其浓度和密度逐渐下降。空气从边缘进入的速度正比于边缘面积和半径变化速度比于边缘面积和半径变化速度dMdMa1a1/dt/dt空气从边缘进入气团的速度，空气从边缘进入气团的速度，1 1边缘空气卷入系数，取边缘空气卷入系数，取0 06 6（1717）空气从气团顶部进入速度可按下式计算空气从气团顶部进入速度可按下式计算(18)(18) dMdMa2a2/dt/dt空气从顶部进入气团的速度，空气从顶部进入气团的速度，2 2一一经验
51、常数，取经验常数，取0.10.1，u ut t水平紊流速度函数，取决于水平紊流速度函数，取决于风速、大气稳定性和地面粗糙度风速、大气稳定性和地面粗糙度, , RiRiRichardsonRichardson数数(19)(19)LsLs紊流长度紊流长度(20)(20)(21)(21)k kVon Von KarmansKarmans常数，取常数，取0.40.4； u u高度高度Z Z处的风处的风速，速，Z Z通常取通常取10m10m； u u* *摩擦速度；摩擦速度； Z Z0 0地面粗地面粗糙度糙度大气稳定度大气稳定度A,B3.03.0C,DC,D2.42.4E,FE,F1.61.6u ut
52、t/u/u* *的取的取值紊流长度紊流长度Ls,Ls,取决于大气稳定度和气团高度取决于大气稳定度和气团高度H,H,取值表如下：取值表如下：u u* */u /u 也可以取典型值也可以取典型值0 01 1，紊流长度，紊流长度L Ls s5.88 H5.88 H0.480.48(4)(4)混合热力学混合热力学 从容器中泄漏的气团密度比空气大，但温度较低，从容器中泄漏的气团密度比空气大，但温度较低，与地面接触时将被加热，如果只考虑自然对流传导，与地面接触时将被加热，如果只考虑自然对流传导，不考虑强制对流和太阳辐射等等，则热方程如下不考虑强制对流和太阳辐射等等，则热方程如下(22)(22) 3 3热传
53、导系数，对液化天然气气团典型值为热传导系数，对液化天然气气团典型值为2.7J2.7Jm m2 2s sK K4/34/3；T T云团温度，云团温度，TgTg地面温度地面温度由空气卷入带进的热量为由空气卷入带进的热量为(23)(23)dMadMa卷入空气的总质量，卷入空气的总质量，dMdMa adMdMa1a1十十dMdMa2a2;Cpa;Cpa一空一空气的热容，气的热容，TaTa空气的温度空气的温度总的热平衡方程为总的热平衡方程为(24)(24) MgMg气团中原泄漏气体的质量，气团中原泄漏气体的质量，CpgCpg泄漏气体的热容泄漏气体的热容若气团中含有液滴，要考虑液体蒸发吸收热量若气团中含有
54、液滴，要考虑液体蒸发吸收热量任意时刻的气团密度按下式计算任意时刻的气团密度按下式计算(25)(25)(5)(5)计算方法计算方法采用分步计算的方法采用分步计算的方法（3 3）根据地面热传导和空气带入的热量，可以计算气）根据地面热传导和空气带入的热量，可以计算气团的温度，从而得到新的气团密度和浓度团的温度，从而得到新的气团密度和浓度（1 1）假定在一小的时间间隔内，例如）假定在一小的时间间隔内，例如0.1s0.1s，扩散速度，扩散速度是不变的，则扩散速度可求是不变的，则扩散速度可求（2 2）假定气团体积不变，则可计算卷入空气的量）假定气团体积不变，则可计算卷入空气的量（4 4）根据新的气团密度和
55、体积，计算新的气团高度。）根据新的气团密度和体积，计算新的气团高度。气团的半径由重力坍塌决定，不受混合和加热的影响气团的半径由重力坍塌决定，不受混合和加热的影响（5 5）进入下一个时间间隔，新的密度、高度和温度）进入下一个时间间隔，新的密度、高度和温度作为初始条件。如此循环直到不符合重气扩散条件作为初始条件。如此循环直到不符合重气扩散条件(6)(6)转换标准转换标准 随着空气进入和温度升高，气团密度逐渐下降，随着空气进入和温度升高，气团密度逐渐下降，密度接近空气密度，紊流扩散开始起主要作用时，就密度接近空气密度，紊流扩散开始起主要作用时，就应该转用中等密度云扩散模型应该转用中等密度云扩散模型(
56、 (高斯扩散模型高斯扩散模型) )。转换。转换标准通常可以用密度标准或标准通常可以用密度标准或RichardsonRichardson数标准，计算数标准，计算过程中要随时检查转换标准过程中要随时检查转换标准密度标准定义为密度标准定义为(26)(26) 如果如果 的的临界值通常取临界值通常取0.0010.001或或0.050.05。小于此值则。小于此值则转换为高斯扩散模型转换为高斯扩散模型4.2 连续扩散的平板模型连续扩散的平板模型 平板模型假定连续扩散的云羽截面呈矩形，中心线平板模型假定连续扩散的云羽截面呈矩形，中心线平行于风向，高为平行于风向，高为H H，宽为，宽为2L;2L;在泄漏源，初始
57、矩形的在泄漏源，初始矩形的高为高为H H0 0，宽为，宽为2L2L0 0 对于液他的连续蒸发，取矩形半宽对于液他的连续蒸发，取矩形半宽L L0 0等于液他半径，等于液他半径，矩形高度就取液池上方用于计算空气进入量的高度。对矩形高度就取液池上方用于计算空气进入量的高度。对于喷射扩散，在喷射扩散结束点处，浓度为中心线于喷射扩散，在喷射扩散结束点处，浓度为中心线1010的点形成一个圆，将圆截面折合为同面积的矩形作为扩的点形成一个圆，将圆截面折合为同面积的矩形作为扩散的初始矩形散的初始矩形重气云羽平板模型扩散示意图重气云羽平板模型扩散示意图 连续泄漏是稳态过程，每一位置的密度和浓度不连续泄漏是稳态过程
58、，每一位置的密度和浓度不随时间变化。也可以假想一矩形薄片，从泄漏源开始随时间变化。也可以假想一矩形薄片，从泄漏源开始沿风向运动，这样就完全可以采用前面处理瞬时泄漏沿风向运动，这样就完全可以采用前面处理瞬时泄漏圆柱体的方法处理矩形薄片，所有算法都是相同的，圆柱体的方法处理矩形薄片，所有算法都是相同的，不同的仅仅是底面积、顶面积和侧面积的计算不同的仅仅是底面积、顶面积和侧面积的计算5.5.中等密度云扩散中等密度云扩散 当泄漏气体或气体与空气混合后的密度接近空气当泄漏气体或气体与空气混合后的密度接近空气密度时，重力下沉与浮力上升作用可以忽略，扩散主密度时，重力下沉与浮力上升作用可以忽略，扩散主要是由
59、空气的湍流决定。在假设均匀湍流场的条件下，要是由空气的湍流决定。在假设均匀湍流场的条件下，有害物质在扩散截面的浓度分布呈高斯分布，所以称有害物质在扩散截面的浓度分布呈高斯分布，所以称为高斯扩散。甚至在重气泄漏场合，考虑到重力作用为高斯扩散。甚至在重气泄漏场合，考虑到重力作用影响时间较短，也可以直接采用高斯模型影响时间较短，也可以直接采用高斯模型5.1 5.1 扩散基本方程和基本条件扩散基本方程和基本条件 在均匀流场中，根据菲克定律和质量守恒，可以在均匀流场中，根据菲克定律和质量守恒，可以建立有害气体的三维扩散基本方程建立有害气体的三维扩散基本方程（2727） C C有害物质浓度，有害物质浓度，
60、E Et,xt,x、E Et,yt,y、E Et,zt,z- -x,y,zx,y,z方向方向上的湍流扩散系数，常数；上的湍流扩散系数，常数； u ux x,u,uy y,u,uz zx x、y zy z方方向上的平均风速，向上的平均风速， K K衰减系数，常数衰减系数，常数 浓度对位置的二阶偏导数项表示有害物质的湍流扩浓度对位置的二阶偏导数项表示有害物质的湍流扩散，一阶偏导数项表示有害物质随风的运动，右边最散，一阶偏导数项表示有害物质随风的运动，右边最后一项表示有害物质的衰减速率后一项表示有害物质的衰减速率 在大气场中，只考虑在大气场中，只考虑x x方向风速、风向与方向风速、风向与x x轴一轴
61、一致，即致，即u uy y0 0、u uz z0 0，同时忽略地面吸收等造成的有，同时忽略地面吸收等造成的有害物质衰减，即害物质衰减，即K K0 0，则得到有风时瞬时泄漏气体，则得到有风时瞬时泄漏气体扩散方程扩散方程(28)(28)无风时，无风时，u ux x=0,=0,得到无风条件下的瞬时泄漏扩散模型得到无风条件下的瞬时泄漏扩散模型（2929） 对于有风条件下连续扩散的情况，认为稳态时对于有风条件下连续扩散的情况，认为稳态时某位置的浓度不随时间变化某位置的浓度不随时间变化扩散方程为扩散方程为忽略在风向上的湍流扩散忽略在风向上的湍流扩散（3030）5.25.2无边界点源模型无边界点源模型 不考
62、虑地面对扩散的影响，认为扩散在三维无限空不考虑地面对扩散的影响，认为扩散在三维无限空间进行，没有边界限制间进行，没有边界限制5.2.15.2.1瞬时点源扩散瞬时点源扩散（1 1）无风条件点源扩散模型）无风条件点源扩散模型沿沿x x方向的一维扩散方向的一维扩散（3131）（3232）初始条件：初始条件：t=0,x=0t=0,x=0处，处，c c； x0,c0x0,c0边界条件：边界条件：t t,c0, -x+,c0, -x+ 解方程可得泄漏的单位源在任意时刻在一维空间解方程可得泄漏的单位源在任意时刻在一维空间的浓度分布的浓度分布（3333）源强为源强为Q Q时的浓度分布时的浓度分布(34)(34
63、) 分别求出在分别求出在x,y,zx,y,z方向上单位泄漏的浓度分布，相乘再方向上单位泄漏的浓度分布，相乘再乘以乘以Q Q，可得无风条件下瞬时泄漏在三维空间任意时刻浓，可得无风条件下瞬时泄漏在三维空间任意时刻浓度分布度分布令令（35）(36)（2 2）有风条件）有风条件 有风条件下，气团中心按风速运动，进行坐标变有风条件下，气团中心按风速运动，进行坐标变换后得浓度方程换后得浓度方程(37) 瞬时扩散模型表明：泄漏的气体为随时间膨大的瞬时扩散模型表明：泄漏的气体为随时间膨大的气团，随着气团增大浓度逐渐降低。有风条件下，气气团，随着气团增大浓度逐渐降低。有风条件下，气团随风运动团随风运动5.2.2
64、 5.2.2 连续点源扩散连续点源扩散有风条件下连续点源扩散的浓度分布有风条件下连续点源扩散的浓度分布（3838）令令（3939） 可见连续泄漏向下风方向形成一条烟羽，随着距离泄可见连续泄漏向下风方向形成一条烟羽，随着距离泄漏口的距离增加，泄漏物质浓度下降。在任意点处，泄漏漏口的距离增加，泄漏物质浓度下降。在任意点处，泄漏物质浓度正比于源强，但不随时间变化物质浓度正比于源强，但不随时间变化5.35.3有界点源扩散有界点源扩散 实际扩散需要考虑地面的影响，一般把地面看做实际扩散需要考虑地面的影响，一般把地面看做镜子，可以反射有害物质。有害物质的实际浓度为由镜子，可以反射有害物质。有害物质的实际浓
65、度为由真实源计算的浓度和由与真实源对称的虚源计算的浓真实源计算的浓度和由与真实源对称的虚源计算的浓度之和度之和 如果泄漏口上空存在逆温层，还要考虑逆温层的如果泄漏口上空存在逆温层，还要考虑逆温层的影响影响5.3.15.3.1高架连续点源高架连续点源 设泄漏源有效高度为设泄漏源有效高度为H H，取其在地面投影为坐标，取其在地面投影为坐标原点，原点，x x轴指向风向。考虑地面反射作用，可得高架连轴指向风向。考虑地面反射作用，可得高架连续点源泄漏的浓度分布（高斯烟羽模型）续点源泄漏的浓度分布（高斯烟羽模型）C C污染物浓度，污染物浓度，Q Q源强，源强，kg/s;kg/s;泄漏高度的平均风泄漏高度的
66、平均风速，速， y y, , z z分别用浓度标准偏差表示分别用浓度标准偏差表示y y轴及轴及z z轴上的扩轴上的扩散参数；散参数；H H泄漏有效高度泄漏有效高度注意：高斯烟羽模型不适用于风速小于注意：高斯烟羽模型不适用于风速小于1m/s1m/s的情况，的情况，低风速扩散时因应修正高斯模型低风速扩散时因应修正高斯模型（4040）根据上面的公式，可计算典型的污染的情况根据上面的公式，可计算典型的污染的情况（1 1）高架连续点源地面浓度，当）高架连续点源地面浓度，当z=0z=0时时（2 2）高架连续点源地面轴向浓度，当）高架连续点源地面轴向浓度，当y=0,z=0y=0,z=0时时(41(41) )
67、（4242）（3 3）高架连续点源地面轴向浓度，当）高架连续点源地面轴向浓度，当y=0,z=0y=0,z=0时，假时，假设设 x x/ / y y=a=a=常数时，对常数时，对 z z求导并令等于求导并令等于0 0可得：可得：（4343）(44)(44)(4)(4)地面连续点源扩散模式，当地面连续点源扩散模式，当H=0H=0时，时，(45(45) )（5 5）地面连续点源轴线的浓度，当）地面连续点源轴线的浓度，当y=0,z=0,H=0y=0,z=0,H=0时时（4646）5.3.25.3.2高架瞬时点源高架瞬时点源 设释放源有效高度为设释放源有效高度为H H，取释放源在地面投影为坐，取释放源在
68、地面投影为坐标原点，进行坐标变换并考虑地面反射作用，可将无标原点，进行坐标变换并考虑地面反射作用，可将无边界瞬时点源扩散模型转换为高架瞬时点源模型边界瞬时点源扩散模型转换为高架瞬时点源模型 对于地面源由于泄漏源在地面上，实源与虚源重合，对于地面源由于泄漏源在地面上，实源与虚源重合，所以，地面瞬时点源泄漏时，有害物质的浓度为无界时所以，地面瞬时点源泄漏时，有害物质的浓度为无界时的的2 2倍倍（1 1）无风时）无风时(47)(47)对于地面点源，对于地面点源，H=0,H=0,则则(48)（2 2）有风时）有风时（4949）对于地面点源，对于地面点源，H=0,H=0,则则(50)5.45.4大气稳定
69、度与扩散参数大气稳定度与扩散参数 在均匀湍流场中，扩散参数与下风向距离的关系是在均匀湍流场中，扩散参数与下风向距离的关系是明确的。但是实测的扩散参数与理论预测差别较大，特别明确的。但是实测的扩散参数与理论预测差别较大，特别是垂直扩散参数差别更大。因此，扩散参数是以实际测定是垂直扩散参数差别更大。因此，扩散参数是以实际测定为基础的为基础的 帕斯奎尔帕斯奎尔( (PasquillPasquill) )根据天空中观测的风速、云量、根据天空中观测的风速、云量、云状和日照等天气资料，将大气的扩散稀释能力分为六云状和日照等天气资料，将大气的扩散稀释能力分为六个稳定度级别，古福德个稳定度级别，古福德(Gif
70、ford)(Gifford)在此基础上建立了扩散在此基础上建立了扩散系数与下风向距离的函数关系，并绘成系数与下风向距离的函数关系，并绘成PGPG曲线图。根曲线图。根据大气稳定级别查图即可知道扩散参数据大气稳定级别查图即可知道扩散参数大气稳定度级别划分大气稳定度级别划分1.A-1.A-极不稳定，极不稳定，B B不稳定，不稳定，C C弱不稳定，弱不稳定，D D中性，中性，E E弱稳定，弱稳定，F F稳定；稳定；2.A2.ABB按按A A、B B数据内插数据内插 3.3.规定日落前定日落前1h1h，日出后，日出后1h1h为夜夜间 4.4.不不论什么天气状况，夜晚前后各什么天气状况，夜晚前后各1h1h
71、算中性算中性 5.5.仲夏晴天中午仲夏晴天中午为强日照，寒冬中午日照，寒冬中午为弱日照；云量：目弱日照；云量：目视估估计云蔽天空的份数云蔽天空的份数 由表可知，不论白天还是夜晚风速越大或云由表可知，不论白天还是夜晚风速越大或云量越多，稳定度越接近于中性；晴天和小风时，白量越多，稳定度越接近于中性；晴天和小风时，白天不稳定，夜晚稳定；过渡时间内都是中性天不稳定，夜晚稳定；过渡时间内都是中性 帕斯奎尔的实验数据来自平坦的草原，因此划分稳帕斯奎尔的实验数据来自平坦的草原，因此划分稳定度的方法对开阔平坦的乡村地区比较可靠，但对城市定度的方法对开阔平坦的乡村地区比较可靠，但对城市地区不大可靠，其原因主要
72、是城市地面粗糙度大及热岛地区不大可靠，其原因主要是城市地面粗糙度大及热岛效应的影响。最大的差别出现在晴天静风的夜间，此时，效应的影响。最大的差别出现在晴天静风的夜间，此时，乡村的大气处于稳定状态，而城市由于建筑物多，且有乡村的大气处于稳定状态，而城市由于建筑物多，且有众多的人工热源，因此在相当于城市建筑物平均高度几众多的人工热源，因此在相当于城市建筑物平均高度几倍之内是中性或微不稳定的，在其上面是稳定层倍之内是中性或微不稳定的，在其上面是稳定层 按照帕斯奎尔稳定度分类法确定大气稳定度时，辐射的按照帕斯奎尔稳定度分类法确定大气稳定度时，辐射的强弱欠缺客观标准，在同一天气状况下，不同的人可能定强弱
73、欠缺客观标准，在同一天气状况下，不同的人可能定出不同的稳定度级别。特纳尔出不同的稳定度级别。特纳尔(Turner)(Turner)提出了一套根据太提出了一套根据太阳高度角和云高、云量确定太阳辐射等级，再由辐射等级阳高度角和云高、云量确定太阳辐射等级，再由辐射等级和和10m10m高处风速确定稳定度级别的方法，这种方法称为帕斯高处风速确定稳定度级别的方法，这种方法称为帕斯奎尔奎尔特纳尔法，简称特纳尔法，简称PTPT法法 由由P-TP-T法确定太阳辐射等级时，云量一项比较复杂，法确定太阳辐射等级时，云量一项比较复杂，且需云高的资料，而中国气象站对云量的观测只分低且需云高的资料，而中国气象站对云量的观
74、测只分低云量和总云量，没有云高的资料，因此针对中国的具云量和总云量，没有云高的资料，因此针对中国的具体情况，对体情况，对PTPT分类法进行了适当的修正。中国有关分类法进行了适当的修正。中国有关环境标准中就是采用修正的环境标准中就是采用修正的PTPT法法(51)(51) 0 0360dn360dn365365； 太阳倾角，太阳倾角，dndn一年中一年中日期序数，日期序数，0 0，1 1，2 2，365365（1）首先算出太阳倾角）首先算出太阳倾角（2 2）根据太阳倾角算出太阳高度角）根据太阳倾角算出太阳高度角h h0 0(52)(52)h h0 0太阳高度角，太阳高度角， 当地纬度，当地纬度，t
75、 t北京时间，北京时间， 当地经度当地经度（3 3）再从表中由太阳高度角）再从表中由太阳高度角h h0 0和云量查出太阳辐射等级和云量查出太阳辐射等级太阳辐射等级太阳辐射等级(4)(4)最后从表中由地面风速和太阳辐射等级查出大气稳定度等级最后从表中由地面风速和太阳辐射等级查出大气稳定度等级大气稳定度等级大气稳定度等级后果分折是一种对假想事故的可能后果描述，并没后果分折是一种对假想事故的可能后果描述，并没有事故发生时间和现场的天气观测资料，因此，在有事故发生时间和现场的天气观测资料，因此，在进行分析时宜考虑当地最常出现天气倍况和对扩散进行分析时宜考虑当地最常出现天气倍况和对扩散最不利的天气情况最
76、不利的天气情况 扩散参数和下风向距离的关系以函数形式表示使扩散参数和下风向距离的关系以函数形式表示使用比较方便，常用的表示方法用比较方便，常用的表示方法(53)(53)5.5 5.5 虚拟点源虚拟点源 假设释放源为点源，即释放源的几何尺寸为零，相应假设释放源为点源，即释放源的几何尺寸为零，相应的浓度为无穷大。释放出的气体随风向下风向运动，同的浓度为无穷大。释放出的气体随风向下风向运动，同时在水平和垂直方向上扩散。一般设备的小尺寸泄漏用时在水平和垂直方向上扩散。一般设备的小尺寸泄漏用点源模型处理是可以接受的。但在某些情况下，直接使点源模型处理是可以接受的。但在某些情况下，直接使用点源模型显然是不
77、行的，一般采用虚拟点源技术处理用点源模型显然是不行的，一般采用虚拟点源技术处理(1)(1)重气团、重气云和液池蒸发重气团、重气云和液池蒸发 重气团或重气云羽通常采用重气模型描述，但是重气团或重气云羽通常采用重气模型描述，但是随着空气进入和温度升高，紊流扩散逐渐占主要地位，随着空气进入和温度升高，紊流扩散逐渐占主要地位，就应该采用高斯扩散模型进一步描述。在转变点，重就应该采用高斯扩散模型进一步描述。在转变点，重气云团是高度为气云团是高度为H H、半径为、半径为R R的圆柱体，重气云羽是截的圆柱体，重气云羽是截面高为面高为H H、半宽为、半宽为L L的矩形，很明显该点不能作为高斯的矩形，很明显该点
78、不能作为高斯模型的原点模型的原点 设想在转变点上风向上某点有一虚拟点源，该点源设想在转变点上风向上某点有一虚拟点源，该点源按高斯模型处理，在转变点处中心线的浓度等于转变按高斯模型处理，在转变点处中心线的浓度等于转变点处云团或云羽的浓度，水平扩散系数为点处云团或云羽的浓度，水平扩散系数为 y0y0，垂直扩，垂直扩散系数为散系数为 z0,z0,因为高斯扩散模式中云团或云羽的浓度呈因为高斯扩散模式中云团或云羽的浓度呈高斯分布，为了与均匀浓度分布衔接，令高斯分布，为了与均匀浓度分布衔接，令 根据扩散系数与下风向距离的关系可以反算出转变根据扩散系数与下风向距离的关系可以反算出转变点到虚拟点源的距离点到虚
79、拟点源的距离x xy y和和x xz z, ,确定虚拟点源以后就可以确定虚拟点源以后就可以按一般的高斯扩散模型进行处理。按一般的高斯扩散模型进行处理。x xy y和和x xz z可以不相等可以不相等(2)(2)闪蒸与绝热扩散气团闪蒸与绝热扩散气团 压缩气体或液化气体的突然释放，气团瞬时绝热膨压缩气体或液化气体的突然释放，气团瞬时绝热膨胀，采用绝热扩散模型处理，可得到扩散结束后半球胀，采用绝热扩散模型处理，可得到扩散结束后半球型气团的半径和密度。将气团半径作为该点的高斯扩型气团的半径和密度。将气团半径作为该点的高斯扩散的水平和垂直扩散系数，根据扩散系数与距离的关散的水平和垂直扩散系数，根据扩散系
80、数与距离的关系计算虚拟点源的位置，然后按高斯点源模型处理系计算虚拟点源的位置，然后按高斯点源模型处理(3)(3)喷射扩散喷射扩散 喷射扩散的转变点处，浓度为中轴线浓度喷射扩散的转变点处，浓度为中轴线浓度1010的圆的圆的直径，作为高斯扩散在此点的扩散系数，用它来计的直径，作为高斯扩散在此点的扩散系数，用它来计算虚拟点源的位置。由于喷射扩散的稀释速度比高斯算虚拟点源的位置。由于喷射扩散的稀释速度比高斯扩散快，所以虚拟点源在泄漏口的上风方向扩散快，所以虚拟点源在泄漏口的上风方向5.6 5.6 计算示例计算示例 用高斯连续地面点源模型计算丙烷泄漏扩散情况，用高斯连续地面点源模型计算丙烷泄漏扩散情况，
81、泄漏流量为泄漏流量为lkglkgs s 将扩散系数将扩散系数( (按书中表按书中表4 41414计算计算) )代人高斯连续地代人高斯连续地面点源模式面点源模式(45)(45)，得到下风向任意点浓度的计算式，得到下风向任意点浓度的计算式，给定浓度可求相应坐标。下图为采用给定浓度可求相应坐标。下图为采用mathematicmathematic绘制绘制的等浓度图，的等浓度图， 图中三条曲线从内到外依次为爆炸下限图中三条曲线从内到外依次为爆炸下限(41g(41gm m3 3) )、中度危害浓度中度危害浓度(18g(18gm m3 3) )和最高允许浓度和最高允许浓度(1g(1gm m3 3) )可以看
82、出，可以看出，风速越大风速越大越有利于越有利于有害气体有害气体扩散；相扩散；相同风速条同风速条件下，大件下，大气越不稳气越不稳定，扩散定，扩散越快，危越快，危害范围超害范围超小小6 6 火灾火灾 火灾大致分为火灾大致分为6 6类，即池火、喷射火、火球、闪火、类，即池火、喷射火、火球、闪火、固体火和普通火灾。作为重大事故后果分析主要考虑固体火和普通火灾。作为重大事故后果分析主要考虑前前4 4种，它们是由于可燃液体或气体泄漏造成的，这种，它们是由于可燃液体或气体泄漏造成的，这类事故在重大事故中占相当局的比例类事故在重大事故中占相当局的比例 火灾后果分析首先需要计算燃烧速度、燃烧时火灾后果分析首先需
83、要计算燃烧速度、燃烧时间、火焰几何尺寸、热辐射通量等，然后计算一定间、火焰几何尺寸、热辐射通量等，然后计算一定距离的人或设备接受的辐射大小，再根据伤害准则距离的人或设备接受的辐射大小，再根据伤害准则评价伤害情况评价伤害情况6.1 6.1 池火池火 池火是一种常见的火灾形式，是可燃液体液面上池火是一种常见的火灾形式，是可燃液体液面上的自然燃烧。泄漏到地面上、堤坝内液体的火灾、敞的自然燃烧。泄漏到地面上、堤坝内液体的火灾、敞开的容器内液体的燃烧等称为池火开的容器内液体的燃烧等称为池火 池火模型一般按圆形液面计算，其他形状的液池火模型一般按圆形液面计算，其他形状的液池应换算为等面积的圆池。对于无边界
84、阻拦的连池应换算为等面积的圆池。对于无边界阻拦的连续泄漏，随着液池面积扩大燃烧速度加快，当燃续泄漏，随着液池面积扩大燃烧速度加快，当燃烧速度等于泄漏速度时，液池直径达到最大烧速度等于泄漏速度时，液池直径达到最大最大直径最大直径D D液池直径，液池直径，Q Q液体泄漏流量，液体泄漏流量，m mf f液体单位面液体单位面积燃烧速率积燃烧速率6.1.16.1.1燃烧速率燃烧速率 当液体沸点高于环境温度时，液体单位面积当液体沸点高于环境温度时，液体单位面积燃烧速率的计算公式为：燃烧速率的计算公式为：(54(54) )当液体沸点低于环境温度时当液体沸点低于环境温度时(55)(55)M Mf f液体单位面
85、积燃烧速率，液体单位面积燃烧速率，c c常数，常数，0.0010.001，H HC C液液体燃烧热，体燃烧热，J/kgJ/kg；H HV V液体在常压沸点下的蒸发热，液体在常压沸点下的蒸发热，CpCp液体的比定压热容，液体的比定压热容，J/(J/(kg.Kkg.K) )；T Tb b液体的沸点，液体的沸点，K K；TaTa环境温度，环境温度，K K 上面的公式忽略了液池大小对燃烧速率的影响，实上面的公式忽略了液池大小对燃烧速率的影响，实际上同样条件下的液体在不同大小的池子中的燃烧速率际上同样条件下的液体在不同大小的池子中的燃烧速率是不同的是不同的考虑液池大小后考虑液池大小后（5656）m mf
86、 f液体单位面积燃烧速率，液体单位面积燃烧速率， m mff液体最大单位面积液体最大单位面积燃烧速率燃烧速率, k, k火焰的吸收衰减系数，火焰的吸收衰减系数， 气体有效厚气体有效厚度校正系数，度校正系数， D D液池直径液池直径 (56)(56)式表明：随着池直径增加，单位面积燃烧速式表明：随着池直径增加，单位面积燃烧速率燃烧速率是增加的，达到最大燃烧速率后保持不变率燃烧速率是增加的，达到最大燃烧速率后保持不变6.1.2 6.1.2 火焰高度火焰高度池火的火焰长度（无风时）池火的火焰长度（无风时）（5757）有风时有风时（5858）L L火焰长度，火焰长度，D D液池直径，液池直径，mfmf
87、液体单位面积燃烧液体单位面积燃烧速率，速率， a a 空气密度，空气密度，g g重力加速度，重力加速度，u u1010米高出米高出风速，风速，ucuc特征风速特征风速（5959）如果如果uuu uc c, , u/uu/uc c取取1 1 （5858）（）（5959）式表明：液池直径越大火焰越长，有）式表明：液池直径越大火焰越长，有风时火焰长度有所减少，但火焰向下风向倾斜，加重风时火焰长度有所减少，但火焰向下风向倾斜，加重了下风向的热辐射危害，可能危及附近高大设备了下风向的热辐射危害，可能危及附近高大设备火焰在风作用下向下风方向扩展，风向上直径为火焰在风作用下向下风方向扩展，风向上直径为（60
88、60）与原液池直径之差称为后拖量与原液池直径之差称为后拖量火焰倾斜角按下式计算火焰倾斜角按下式计算（6161） 空气运动粘度，空气运动粘度， 火焰倾角火焰倾角风向风向D D有风时池火示意图有风时池火示意图在计算火焰的公式中，被认为效果较好的关系式为：在计算火焰的公式中，被认为效果较好的关系式为：（6262）6.1.3 6.1.3 热辐射通量热辐射通量池周围距池中心池周围距池中心x x处的热辐射强度为：处的热辐射强度为：（6363）q q接受点热辐射通量，接受点热辐射通量，E E池火表面热辐射通量，池火表面热辐射通量， F F几何视角因子，几何视角因子， 大气透射率大气透射率 视角因子：视角因子
89、是热辐射传递的重要概视角因子：视角因子是热辐射传递的重要概念，是接受体所能接受的发热体辐射能量的分数，念，是接受体所能接受的发热体辐射能量的分数，大小取决于发射体和接受体的形状、距离和相对角大小取决于发射体和接受体的形状、距离和相对角度度(64)(64)这需要对整个火焰表面积分，这需要对整个火焰表面积分，结果往往非常复杂结果往往非常复杂池火视角因子计算示意图池火视角因子计算示意图(65)(65)其中其中 考虑水蒸气、二氧化碳等的影响，辐射的热量在传考虑水蒸气、二氧化碳等的影响，辐射的热量在传播过程中有损失，用大气透射率表示。常用的一种计算播过程中有损失，用大气透射率表示。常用的一种计算方法是方
90、法是(66)(66)另外一种计算方法是另外一种计算方法是 如果按点源模型处理池火，假设全部热量由池如果按点源模型处理池火，假设全部热量由池中心点发出，则下面的简单模型是被广泛使用的中心点发出，则下面的简单模型是被广泛使用的(67(67) )（6868）通常可假定大气透射率为通常可假定大气透射率为1 1 假设燃料燃烧的能量从圆柱状池火焰的侧面和假设燃料燃烧的能量从圆柱状池火焰的侧面和上面均匀向外辐射，则池火焰表面热辐射能量为上面均匀向外辐射，则池火焰表面热辐射能量为(69)(69)f f热辐射系数，范围为热辐射系数，范围为0.130.13一一0.350.35，保守取值，保守取值为为0.350.3
91、5考虑黑烟以及一氧化碳、水蒸气等，火焰表面热辐考虑黑烟以及一氧化碳、水蒸气等，火焰表面热辐射能量可用射能量可用(70)(70)E Ef f火焰可见部分的最大发射能量，取火焰可见部分的最大发射能量，取140kw140kwm m2 2E Es s火焰黑烟部分的最大发射能量，取火焰黑烟部分的最大发射能量，取20kw20kwm m2 2D D池直径池直径6.1.4 6.1.4 热辐射破坏准则与伤害模型热辐射破坏准则与伤害模型1.1.破坏准则破坏准则（1 1）热通量准则）热通量准则 以热通量作为衡量目标是否被破坏的参数，当以热通量作为衡量目标是否被破坏的参数，当目标接受到的热通量大于或等于引起目标破坏所
92、需目标接受到的热通量大于或等于引起目标破坏所需要的临界热通量时，目标被破坏。若目标被破坏目要的临界热通量时，目标被破坏。若目标被破坏目标不被破坏，热通量准则的适用范围为热通量作用标不被破坏，热通量准则的适用范围为热通量作用的时间比目标达到热平衡所需的时间长的时间比目标达到热平衡所需的时间长(2)(2)热强度准则热强度准则 强度准则以目标接受到的热强度作为目标是否被破强度准则以目标接受到的热强度作为目标是否被破坏的唯一参数，当目标接受到的热强度大于或等于目标坏的唯一参数，当目标接受到的热强度大于或等于目标破坏的临界热强度时，目标被破坏。若目标不被破坏，破坏的临界热强度时，目标被破坏。若目标不被破
93、坏，热强度准则的适用范围为作用于目标的热通量持续时间热强度准则的适用范围为作用于目标的热通量持续时间非常短，以至于目标接收到的热量来不及散失掉非常短，以至于目标接收到的热量来不及散失掉（3 3）热通量热强度准则）热通量热强度准则 当热通量准则或热强度准则的适用条件均不具备时，当热通量准则或热强度准则的适用条件均不具备时，应使用热通量热强度准则：目标能否被破坏不能由热通应使用热通量热强度准则：目标能否被破坏不能由热通量或热强度一个参数决定，必须由它们的组合来决定量或热强度一个参数决定，必须由它们的组合来决定 热通量热通量q q作横坐标，热强度作横坐标，热强度Q Q作纵坐标，则目标破坏作纵坐标，则
94、目标破坏的临界状态对应的临界状态对应q qQ Q平面上的一条临界曲线，曲线的平面上的一条临界曲线，曲线的右上方为破坏区，左下方为无伤害区右上方为破坏区，左下方为无伤害区热辐射伤害区示意图热辐射伤害区示意图由图看出：热通量准则和热强度准则都是热通量热由图看出：热通量准则和热强度准则都是热通量热强度准则的极限情况强度准则的极限情况2.2.热辐射伤害概率模型热辐射伤害概率模型概率与伤害百分率的关系为概率与伤害百分率的关系为 当当p pr r=5=5时，伤害百分率为时，伤害百分率为5050% %（7171）皮肤裸露时的死亡概率皮肤裸露时的死亡概率有衣服保护时（有衣服保护时（2020% %皮肤裸露皮肤裸
95、露）的死亡概率）的死亡概率（7272）(73)(73)有衣服保护时（有衣服保护时（2020% %皮肤裸露皮肤裸露）的二度烧伤概率）的二度烧伤概率(74)(74)有衣服保护时（有衣服保护时（2020% %皮肤裸露皮肤裸露）的一度烧伤概率）的一度烧伤概率(75)(75) 根据人体接收的热辐射通量和暴露时间，按上面的公根据人体接收的热辐射通量和暴露时间，按上面的公式计算伤害概率，从而可以进一步计算伤害百分率。后果式计算伤害概率，从而可以进一步计算伤害百分率。后果分析中通常采用相反的方法，即给定伤害概率，确定暴露分析中通常采用相反的方法，即给定伤害概率，确定暴露时间，根据上面的式子计算热辐射通量，根据
96、热辐射通量时间，根据上面的式子计算热辐射通量，根据热辐射通量和距离的关系算出距火源的距离，此距离即为相应的伤害和距离的关系算出距火源的距离，此距离即为相应的伤害距离距离 分析过程中通常都按分析过程中通常都按5050伤害率计算，例如按伤害率计算，例如按5050死亡率划定死亡范围，该范围表明范围内、外死亡死亡率划定死亡范围，该范围表明范围内、外死亡人数各占一半，也可以认为死亡范围内人员全部死亡，人数各占一半，也可以认为死亡范围内人员全部死亡，范围外无一人死亡，这样可以使问题得以简化范围外无一人死亡，这样可以使问题得以简化对于财产损失，可以按引燃木材所需热通量计算对于财产损失，可以按引燃木材所需热通
97、量计算（7676）6.2.6.2.喷射火喷射火 高压气体从裂口高速喷出后被点燃，形成喷射高压气体从裂口高速喷出后被点燃，形成喷射火。喷射火的长度可以认为等于喷口到燃烧浓度下火。喷射火的长度可以认为等于喷口到燃烧浓度下限的长度，热量认为是从中心轴线上一系列相等的限的长度，热量认为是从中心轴线上一系列相等的辐射源发出，每一点源的热通量为辐射源发出，每一点源的热通量为(77)(77)f f燃烧效率因子，取燃烧效率因子，取0.350.35n n假设的点源数；假设的点源数；Q Q泄漏流量，泄漏流量，kgkgs s距离点源距离点源x xi i处某点接收的热辐射通量为处某点接收的热辐射通量为(78)(78)
98、 x xp p发射因子，取发射因子，取0.20.2总热通量是各点辐射的和总热通量是各点辐射的和 辐射点源的数目可以任意选取，但对于后果分析来辐射点源的数目可以任意选取，但对于后果分析来说，取说，取5 5点就可以了。该模型没有考虑风的影响，因为点就可以了。该模型没有考虑风的影响，因为一般喷射速度比风速大得多。在低压喷射时，风速的影一般喷射速度比风速大得多。在低压喷射时，风速的影响比较明显，在下风向接收热量会更多。如果风使喷射响比较明显，在下风向接收热量会更多。如果风使喷射火焰偏离了轴线，则该模型不适用火焰偏离了轴线，则该模型不适用(79)(79)6.3 6.3 火球火球 火球，也称为沸腾液体扩展
99、蒸汽爆炸当压力容器受外火球，也称为沸腾液体扩展蒸汽爆炸当压力容器受外界热量的作用位槽壁强度下降并突然破坏，储存的过热液界热量的作用位槽壁强度下降并突然破坏，储存的过热液体或液化气体突然释放并被点燃，形成巨大火球。火球的体或液化气体突然释放并被点燃，形成巨大火球。火球的危害主要是热辐射而不是爆炸冲击波，强烈的热辐射可能危害主要是热辐射而不是爆炸冲击波，强烈的热辐射可能造成严重的人员伤亡和财产损失造成严重的人员伤亡和财产损失火球直径和火球持续时间的计算公式火球直径和火球持续时间的计算公式(80)(80)(81(81) )W W火球中消耗的可燃物的质量火球中消耗的可燃物的质量 对于单罐储存，对于单罐
100、储存，W W取罐容量的取罐容量的5050；对于双罐储；对于双罐储存，存，W W取罐容量的取罐容量的7070；对于多罐储存，；对于多罐储存，W W取罐容量的取罐容量的9090距火球在地面投影处距火球在地面投影处x x的热辐射通量为的热辐射通量为(82)(82)火球表面热辐射通量为火球表面热辐射通量为（8383）设设f f是燃烧辐射分数，是容器压力的函数是燃烧辐射分数，是容器压力的函数(84)(84) 常数常数f f1 10.27 0.27 ， f f2 20.320.32；p p为容器内压力，单位是为容器内压力，单位是MPaMPa；在没有可靠数据时；在没有可靠数据时f f可取可取0.30.3 大
101、气透射率的取值与池火讨论的相同。考虑最简大气透射率的取值与池火讨论的相同。考虑最简单、最保守的情况，即辐射接收体的平面垂直于接收单、最保守的情况，即辐射接收体的平面垂直于接收体与火球中心连线，视角因子按下式计算体与火球中心连线，视角因子按下式计算（8585）火球辐射示意图火球辐射示意图由图可知：由图可知：（8686）把（把（8383）、）、（8585）（）（8686）代入）代入（8282）得：）得：(87)(87)忽略火球高度的影响，上式可化为简洁形式忽略火球高度的影响，上式可化为简洁形式(88)(88)6.4 6.4 闪火闪火 闪火是可燃蒸汽云的非爆炸燃烧，燃烧速度虽然闪火是可燃蒸汽云的非爆
102、炸燃烧，燃烧速度虽然很快但比爆炸慢得多，危害主要是热辐射面没有冲击很快但比爆炸慢得多，危害主要是热辐射面没有冲击波波 一般后果分析可不考虑闪火的热辐射效应，只考虑一般后果分析可不考虑闪火的热辐射效应，只考虑闪火范围内的伤害。一般可认为闪火范围内的室外人员闪火范围内的伤害。一般可认为闪火范围内的室外人员将全部烧死，建筑物内特有部分人被烧死。在缺乏资料将全部烧死，建筑物内特有部分人被烧死。在缺乏资料时，可以认为室内的死亡率为零时，可以认为室内的死亡率为零6.5 6.5 计算示例计算示例计算计算l00tl00t丁二烯储罐沸腾液体扩展蒸汽云爆炸丁二烯储罐沸腾液体扩展蒸汽云爆炸( (火球火球) )伤害区
103、伤害区火球直径：火球直径：5.85.8（0.50.5100000100000）1/31/3=213.7m=213.7m火球持续时间：火球持续时间：=0.45=0.45(0.5(0.5100000)100000)1/31/3=16.6s=16.6s伤害半径按伤害半径按5050伤害率计算。根据死亡概率与热辐射计量的关系伤害率计算。根据死亡概率与热辐射计量的关系p pr r5 5时，代人火球持续时间，计算热辐射通量时，代人火球持续时间，计算热辐射通量q q=28710.5W/m=28710.5W/m2 2由由f=0.3f=0.3， 1 10.058lnx,Hc=50409KJ/Kg0.058lnx,
104、Hc=50409KJ/Kg代入上式，可求死亡半径代入上式，可求死亡半径X=291mX=291m同理，二度烧伤的热辐射通量：同理，二度烧伤的热辐射通量：q=19551.3W/mq=19551.3W/m2 2, ,伤伤害半径害半径x=350mx=350m一度烧伤的热辐射通量：一度烧伤的热辐射通量：q=8574.5W/mq=8574.5W/m2 2, ,伤害半径伤害半径x=519mx=519m财产损失的热辐射通量：财产损失的热辐射通量：q=26111.1W/mq=26111.1W/m2 2, ,伤害半径伤害半径x=305mx=305m7.7.爆炸爆炸7.17.1分类分类（1 1）物理爆炸锅炉爆炸。压
105、缩气体和液化气体钢瓶）物理爆炸锅炉爆炸。压缩气体和液化气体钢瓶爆炸是典型的物理爆炸，物理爆炸影响范围较小，一爆炸是典型的物理爆炸，物理爆炸影响范围较小，一般不作为重大危险进行后果分析般不作为重大危险进行后果分析7.2 7.2 爆炸伤害准则爆炸伤害准则(1)(1)超压准则超压准则（2 2）化学爆炸两大类。化学爆炸包括不稳定固体或气）化学爆炸两大类。化学爆炸包括不稳定固体或气体的爆炸、受限空间内可燃混合气体爆炸和开放空间体的爆炸、受限空间内可燃混合气体爆炸和开放空间的可燃气体蒸汽云爆炸等，作为重大事故后果分析，的可燃气体蒸汽云爆炸等，作为重大事故后果分析，最重要的是可燃气体泄漏引起的开放空间蒸汽云
106、爆炸最重要的是可燃气体泄漏引起的开放空间蒸汽云爆炸 爆炸波是否对目标造成伤害是由爆炸波超压惟一决爆炸波是否对目标造成伤害是由爆炸波超压惟一决定的，只有当爆炸波超压大于某一临界值时，才会对目定的，只有当爆炸波超压大于某一临界值时，才会对目标造成一定伤害。很明显，超压准则没有考虑正相持续标造成一定伤害。很明显，超压准则没有考虑正相持续时间。理论和实验都表明，爆炸破坏效应不仅与爆炸超时间。理论和实验都表明，爆炸破坏效应不仅与爆炸超压有关，也与超压持续时间有关，持续时间长则破坏更压有关，也与超压持续时间有关，持续时间长则破坏更大。尽管如此，由于爆炸波超压容易测量和估计，所以大。尽管如此，由于爆炸波超压
107、容易测量和估计，所以超压准则是衡量爆炸超压准则是衡量爆炸 破坏效应最常用的准则破坏效应最常用的准则 (2)(2)冲量准则冲量准则 由于伤害效应不仅取决于爆炸波超压，而且与爆炸波由于伤害效应不仅取决于爆炸波超压，而且与爆炸波持续时间有关，爆炸波冲量就是超压和持续时间的函数，持续时间有关，爆炸波冲量就是超压和持续时间的函数，因此用爆炸冲量衡量伤害后果是合理的。冲量准则是指爆因此用爆炸冲量衡量伤害后果是合理的。冲量准则是指爆炸波能否对目标造成伤害，完全取决于爆炸波冲量大小，炸波能否对目标造成伤害，完全取决于爆炸波冲量大小，如果冲量大于临界值，则目标被破坏。但是，如果冲量大于临界值，则目标被破坏。但是
108、， 有一点是有一点是明显的，对于一个很小的超压，作用时间再长也不会产生明显的，对于一个很小的超压，作用时间再长也不会产生任何伤害。因此，仅考虑冲量也是不完全的任何伤害。因此，仅考虑冲量也是不完全的 (3)(3)超压一冲量准则超压一冲量准则 超压一冲量准则综合考虑了超压和冲量两个方面，超压一冲量准则综合考虑了超压和冲量两个方面，如果超压和冲量的共同作用满足某一如果超压和冲量的共同作用满足某一 临界条件，目临界条件，目标就被破坏。超压一冲量准则可以用下式表示标就被破坏。超压一冲量准则可以用下式表示 P P超压，超压，p pcrcr临界超压临界超压, ,I I冲量，冲量，I Icrcr临界冲量，临界
109、冲量，C C常数，与目标性质和破坏等级有关常数，与目标性质和破坏等级有关 超压准则和冲量准则可以视为超压一冲量准则的超压准则和冲量准则可以视为超压一冲量准则的两个极限情况。当冲量小时，伤害主要由超压决定；两个极限情况。当冲量小时，伤害主要由超压决定；当超压小时，伤害主要由冲量决定当超压小时，伤害主要由冲量决定 (89)(89)7.3 7.3 蒸汽云爆炸伤害模型蒸汽云爆炸伤害模型（1 1）TNTTNT当量法当量法 大量泄漏气体延迟点燃可能发生闪火也可能发生爆炸，大量泄漏气体延迟点燃可能发生闪火也可能发生爆炸，这取决于蒸汽的状态和环境条件。空旷的开放空间倾向这取决于蒸汽的状态和环境条件。空旷的开放
110、空间倾向于发生闪火，封闭的或受限空间倾向于发生爆炸于发生闪火，封闭的或受限空间倾向于发生爆炸 蒸汽云爆炸的能量常用蒸汽云爆炸的能量常用TNTTNT当量描述，即将参与当量描述，即将参与爆炸的可燃气体释放的能量折合为能释放相同能量的爆炸的可燃气体释放的能量折合为能释放相同能量的TNTTNT炸药的量，这样，就可以利用有关炸药的量，这样，就可以利用有关TNTTNT爆炸效应的爆炸效应的实验数据预测蒸汽云爆炸效应实验数据预测蒸汽云爆炸效应TNTTNT当量计算公式当量计算公式(90)(90)W WTNTTNT蒸汽云的蒸汽云的TNTTNT当量，当量，kgkg； W Wf f蒸汽云中燃料蒸汽云中燃料的总质量，的
111、总质量，kgkg； 蒸汽云爆炸的效率因子，表明蒸汽云爆炸的效率因子，表明参与爆炸的可燃气体的分数一般取参与爆炸的可燃气体的分数一般取3 3或或4 4； Q Qf f蒸汽的燃烧热，蒸汽的燃烧热，MJ/KgMJ/Kg；Q QTNTTNTTNTTNT的爆炸热，一般取的爆炸热，一般取4.52MJ4.52MJkgkg爆炸中心与给定超压间的距离可以按下式计算爆炸中心与给定超压间的距离可以按下式计算(91)(91)x x距离，距离，m m； p p超压，超压，psi(1psipsi(1psi6.9kPa)6.9kPa) 爆炸涉及的总能量中只有一小部分真正对爆炸有爆炸涉及的总能量中只有一小部分真正对爆炸有贡献
112、，这一分数称为效率因子。效率因子是爆炸后果贡献，这一分数称为效率因子。效率因子是爆炸后果分析中最重要也是最难准确知道的参数，其范围为分析中最重要也是最难准确知道的参数，其范围为2 2一一2020。对于多数脂肪烃，通常推荐值是。对于多数脂肪烃，通常推荐值是3 3；对于某；对于某些烯烃，观察到大约是些烯烃，观察到大约是6 6。含氧燃料趋向于高的效率。含氧燃料趋向于高的效率因子，可以达到因子，可以达到l6l6一一1818（2 2）爆炸伤害概率模型）爆炸伤害概率模型 用概率模型描述超压造成的轻、重伤以及死亡用概率模型描述超压造成的轻、重伤以及死亡情况。模型如下：情况。模型如下：P Pr r概率，概率，
113、 p p超压超压 根据超压一冲量准则根据超压一冲量准则 和概率模型得到的死亡和概率模型得到的死亡半径公式为：半径公式为：(92)(92)（9393）财产损失半径计算公式财产损失半径计算公式7.4 7.4 计算示例计算示例计算计算100100吨丁二烯储罐蒸汽云爆炸伤害区吨丁二烯储罐蒸汽云爆炸伤害区先计算先计算TNTTNT当量当量（9494） 根据伤害的超压计算伤害半径。死亡半径按超压根据伤害的超压计算伤害半径。死亡半径按超压90kp90kpa a计算，重伤半径按计算，重伤半径按44kp44kpa a计算，轻伤半径按计算，轻伤半径按17kp17kpa a计算，财产损失半径按计算，财产损失半径按13
114、.8kp13.8kpa a计算计算计算公式为：计算公式为： 死亡半径：死亡半径：x=104mx=104m；重伤半径：；重伤半径：x=155mx=155m；轻伤半；轻伤半径：径：x=278mx=278m；财产损失半径：；财产损失半径：x=319mx=319m 根据根据(93)(93)计算的死亡半径是计算的死亡半径是62m62m；根据；根据(94)(94)计算的计算的财产损失半径是财产损失半径是180m 180m 8.8.中毒中毒8.1 8.1 毒性释放的影响毒性释放的影响 人接触有毒物质可发生急性的或慢性的中毒。长期人接触有毒物质可发生急性的或慢性的中毒。长期接触低剂量有毒物质造成慢性中毒，而有
115、毒气体泄漏接触低剂量有毒物质造成慢性中毒，而有毒气体泄漏一般是高浓度而持续时间很短，通常引起急性中毒。一般是高浓度而持续时间很短，通常引起急性中毒。事故后果分析中一般只考虑急性中毒事故后果分析中一般只考虑急性中毒急性中毒的表现急性中毒的表现（1 1）刺激）刺激 呼吸系统、皮肤、眼睛等可以感觉到刺激，呼吸系统、皮肤、眼睛等可以感觉到刺激，一些物质在低浓度时就产生刺激这样可以提醒人们寻一些物质在低浓度时就产生刺激这样可以提醒人们寻求防护求防护 （2 2）麻醉）麻醉 有些物质能影响人的反应能力，使采取有些物质能影响人的反应能力，使采取防护措施或提醒别人时反应迟钝防护措施或提醒别人时反应迟钝（3 3）
116、窒息）窒息 大多数气体由于取代了空气中的氧气而造大多数气体由于取代了空气中的氧气而造成窒息。有些物质，如一氧化碳可以置换血液中的氧从成窒息。有些物质，如一氧化碳可以置换血液中的氧从而阻止氧进入人体组织而阻止氧进入人体组织（4 4）系统损害）系统损害 有些物质能损害人体器官，其损害有些物质能损害人体器官，其损害可能是暂时的，也可能是永久的可能是暂时的，也可能是永久的8.28.2毒物影响的后果分析毒物影响的后果分析分析毒物的损害常用概率函数法分析毒物的损害常用概率函数法 p pr r概率函数；概率函数； c c接触有毒物质的浓度，接触有毒物质的浓度，10106 6；t t接触时间，接触时间，min
117、min，A A、B B、n n与特定毒物有与特定毒物有关的常数；关的常数；c cn nt t是毒物接触剂量或称毒性负荷，取决是毒物接触剂量或称毒性负荷，取决于接触浓度和暴露时间于接触浓度和暴露时间(95)(95) 对于连续释放，可认为某一地点的毒物浓度是不变对于连续释放，可认为某一地点的毒物浓度是不变的，则根据暴露时间的，则根据暴露时间( (一般不超过一般不超过30min30min) )计算接触剂量计算接触剂量 对于瞬时泄漏，某一地点的浓度是随时间变化的，对于瞬时泄漏，某一地点的浓度是随时间变化的，那么剂量就应该是浓度对时间的积分。简化的计算是将那么剂量就应该是浓度对时间的积分。简化的计算是将
118、暴露时间分成小段，每一时间段内认为浓度不变，从而暴露时间分成小段，每一时间段内认为浓度不变，从而得到这一时间段的剂量，总暴露剂量为各时间段之和得到这一时间段的剂量，总暴露剂量为各时间段之和 分析有毒气体泄漏的有害区时，可以按照最低有害分析有毒气体泄漏的有害区时，可以按照最低有害浓度、半致死浓度或者不同死亡率浓度、半致死浓度或者不同死亡率( (如如5050死亡率死亡率) )相相应的毒物剂量应的毒物剂量 计算毒物剂量时，暴露时间可以采用实际毒物泄漏计算毒物剂量时，暴露时间可以采用实际毒物泄漏时间，但是，暴露时间一般不超过时间，但是，暴露时间一般不超过30min30min，因为在此时，因为在此时问内有关人员足以逃离现场问内有关人员足以逃离现场
