Numerical simulation on prediction model of risk range of typical gas release through small holes
-
摘要:
为了缓解我国石油资源紧张的现状，发展合成油工业作为我国石油资源的补充是合理的途径，而在化工合成中，大多数反应都需要在加压环境中进行。泄漏是最常见的事故之一，而现阶段对于加压装置、管道气体泄漏的研究还有待完善。因此，通过CFD(计算流体动力学)数值模拟的方法，以化工合成中常见的气体（H2、CO、CH4）为研究对象，气体泄漏后危险浓度的范围作为气体泄漏危险范围，利用Fluent软件开展不同工况下气体泄漏扩散模拟，以研究不同气体、不同泄漏压力、不同泄漏温度以及不同泄漏口直径情况下气体泄漏扩散规律及泄漏危险范围，并构建泄漏危险范围预测模型。研究表明，在其他因素不变的条件下，随着泄漏压力的增大，气体泄漏危险范围随之增加；随着气体温度的增大，气体泄漏危险范围随之明显降低；随着泄漏口直径的增大，气体泄漏危险范围随之明显增加。在3个变量中，泄漏口尺寸的增大对于危险范围的扩大最为明显，当气体泄漏压力由1 MPa提高到2 MPa时，3种气体泄漏危险范围增加10%左右；当气体温度由300 K提高到400 K时，3种气体泄漏危险范围同样增加10%左右，而当泄漏口直径由2 mm增加到3 mm时，3种气体泄漏危险范围增加50%。通过对于3种气体（H2、CO、CH4）数值模拟结果中的泄漏危险范围的长度、宽度数据的定量分析，并对于参量进行无量纲处理，最终分别得出3种气体（H2、CO、CH4）泄漏危险范围（长度、宽度）的预测模型，其中泄漏危险范围长度和宽度为压力比和温度比的函数，预测结果可为泄漏事故救援和应急处理提供参考依据。
Abstract:In order to relieve the shortage of petroleum resources, it is a reasonable way to develop synthetic oil industry as a supplement of petroleum resources. In chemical synthesis, most reactions need to be carried out in pressurized environment. Leakage is one of the most common accidents, and the research of gas leakage in pressurizing device and pipeline is still to be improved. Therefore, by CFD numerical simulation method, this paper takes common gases (H2, CO, CH4) in chemical synthesis as the research object, and uses the range of dangerous concentration after gas leakage as the danger range of gas leakage. Fluent software is used to simulate gas leakage and diffusion under different working conditions. In order to study the gas leakage diffusion law and leakage danger range under different gas, different leakage pressure, different leakage temperature and different leakage port diameter, and build the prediction model of leakage danger range. The results show that the danger range of gas leakage increases with the increase of leakage pressure when other factors remain unchanged. With the increase of gas temperature, the danger range of gas leakage decreases obviously. With the increase of the diameter of the leakage port, the danger range of gas leakage increases obviously. Among the three variables, the increase of leakage port size has the most obvious effect on the expansion of the danger range. When the gas leakage pressure increases from 1 MPa to 2 MPa, the danger range of three kinds of gas leakage increases by about 10%. When the gas temperature increases from 300 K to 400 K, the danger range of three kinds of gas leakage also increases by about 10%, and when the diameter of the leakage port increases from 2 mm to 3 mm, the danger range of three kinds of gas leakage increases by 50%. By quantitative analysis of the length and width data of the leakage danger range in the numerical simulation results of three kinds of gases (H2, CO, CH4), and dimensionless treatment of the parameters, the prediction models of the leakage danger range (length and width) of three kinds of gases (H2, CO, CH4) are finally obtained. The length and width of the leakage danger range are functions of pressure ratio and temperature ratio. The prediction results can provide reference for the rescue and emergency treatment of leakage accidents.
-
Keywords:
- chemical synthesis /
- gas leakage /
- CFD /
- scope of danger /
- predictive models
-
0. 引 言
能源是国家的重要基础产业，能源安全是国家安全的重要组成部分[1-3]，我国是一个缺油、少气，煤炭资源相对丰富的国家[4]，随着时代的发展和国民经济的高速增长，石油消费增速远大于国内石油资源增产的速度，发展合成油工业作为石油资源的补充是合理的[5]。各种化工合成工艺所用原料气的有效成分主要是H2、CO和CH4，而在各种工艺中，大多数都需要在高压环境中进行，如合成油品所需压力为2~3 MPa，合成二甲醚则需要2~4 MPa压力[6]，因此合成装置、运输管道一旦发生泄漏，造成人员及经济损失不可估量。
数值模拟能够表示参数的变化关系，使用数值计算方法进行求解，能够得到过程中的动态模拟结果，极大地促进了科学研究的发展。刘彦青[7]在同时考虑瓦斯含量与压力的情况下，通过数值模拟方法计算得到了使煤层瓦斯涌出量最低的工作面平均推进速度周期为25~30 d。在高压气体泄漏扩散的研究中，李伟等[8]通过流体力学计算软件，对于高压甲烷储罐的泄漏喷射和扩散过程进行了数值模拟，通过数值模拟结果，建立了甲烷容器内压力与破裂口喷射速度的定量公式。在甲烷高压球罐的泄漏研究中，傅杰[9]利用Fluent软件，建立了较为准确的高压气体泄漏扩散的数值模型，并利用该模型探究了甲烷在不同泄漏条件下泄漏口附近以及气体泄漏远端的气体浓度、压力等物理参数的规律。刘延雷等[10]，通过数值模拟研究了管道内高压氢气、天然气的不同扩散特性，发现高压的氢气泄漏时产生的危险气体云团较大，且高压氢气泄漏时初始的泄漏速度远大于高压天然气泄漏时的初始速度，在地面附近区域氢气泄漏扩散产生的危险要小于天然气。同样是在天然气管道气体泄漏的研究中，周伟国等[11]通过数值模拟，对管道发生泄漏后泄漏压力、环境风速对于甲烷气体的扩散特性影响进行了研究，结果表明，管道内压力越大时，泄漏口气体速度越大，气体扩散范围也越大，当环境风速越大时，甲烷气体射流的偏转角度越大，扩散速率也越大。高嵩等[12]，利用Fluent软件研究了高压天然气管道泄漏及其规律，通过建立的天然气长管道运输模型，对不同条件下管道泄漏的压力、浓度、速度及危险区域进行了分析，根据结果对比分析了不同泄漏孔径、不同泄漏压力对于气体泄漏范围的影响。徐坤等[13]利用Fluent软件模拟了在机械通风条件下不同泄漏口方向、位置、泄漏压力等条件下燃气浓度分布特性，发现泄漏口的方向对于甲烷气体泄漏初期的浓度分布具有重大影响。刘敏鸿等[14]利用数值模拟方法，分别研究了天然气在封闭和半封闭小空间内的泄漏扩散规律，发现天然气在封闭空间泄漏发生爆炸的风险小于半封闭空间。袁裕鹏等[15]利用Fluent软件，对不同泄漏条件下，氢气燃料电池仓内氢气泄漏扩散规律进行了探究，并发现氢气探测报警器适宜安装到空间顶部角落。汪建平等[16]采用数值模拟的方法，研究了天然气在厨房内泄漏扩散过程，并且将实验与数值模拟结果进行对比，验证了数值模拟的准确性。万宇飞等[17]利用Fluent软件，研究了不同风向条件对于高压LPG球形储罐中丙烷泄漏扩散规律的影响，可为发生泄漏时人员自救及消防救援提供参考。朱静等[18]采用数值模拟方法，模拟了天然气在典型住宅内泄漏扩散规律及自然通风对气体扩散的影响，结果表明自然通风可以很好地防止天然气聚集，降低爆炸风险。喻健良等[19]通过构建高压CO2管道泄漏扩散数值模型并根据实验结果进行验证，总结出不同大小泄漏孔下CO2浓度和温度分布规律。周宁等[20]利用Fluent软件，研究了环境温度对于LNG泄漏扩散规律的影响，发现环境温度较高时，甲烷扩散最远距离大于环境温度低时。程方明等[21]利用FLACS软件，对储罐内高压天然气非恒定速率泄漏进行了模拟，研究发现风速越大，气体泄漏扩散产生的气云稳定所需的时间越长。陶婧莹等[22]采用数值模拟的方法，利用高斯烟羽模型对LNG泄漏扩散进行研究，发现在泄漏过程中，随着环境风速的增加，危险区域减小；随着大气稳定度的增加，危险区域增加。周宁等[23]利用Fluent软件，对于不同条件下管道内丁烷气体泄漏扩散规律进行了研究，发现气体初始速度越大，自由扩散作用越强，危险范围也越广。张晓蕾等[24]采用PHAST模拟软件，研究了不同工况下风速大小对周围建筑物及人员的热辐射影响，研究发现风速越大，热辐射影响的范围越大。
在上述对于高压气体泄漏的研究中可以发现，由于高压气体泄漏实验危险性较大，实验成本较高，因此许多研究采取了数值模拟的方法，且模拟结果可靠，在之前的大多数研究中，目标气体大多数集中在甲烷和天然气，对于化工合成中常见的气体，如H2、CO、CH4在不同泄漏条件下，泄漏规律及危险范围的研究并不充分，且通常轻视小尺寸泄漏口的危险性，因此开展常见可燃气体小孔泄漏扩散研究是十分必要的。
1. 不同条件下可燃气体泄漏扩散规律及泄漏危险范围
1.1 CFD数值模拟参数设置
本次数值模拟选择CFD数值模拟方法，CFD（Computational Fluid Dynamics）全称为计算流体动力学，是通过计算机模拟流体的特定情境以得到流体相关参数并进行研究。选用ICEM软件进行建模及网格划分，Fluent软件进行高压气体泄漏的模拟，CFD-POST软件进行数据后处理。在模拟过程中，求解器选择压力基求解器Pressure-Based，湍流模型选择考虑可压性影响的双方程k-omega SST模型，该模型是一种在工程上得到广泛应用的混合模型，即在近壁面保留了原始k-ω模型，在远离壁面的位置应用了k-ε模型[25]，其方程与模型具体如下：
可压缩流体质量、动量、能量和组分守恒方程[26]可表示为
$$ \frac{\partial \overline{\rho }}{\partial t}+\frac{\partial }{\partial {x}_{j}}\left(\overline{\rho }{\tilde{ u}}_{j}\right)=0 $$ (1) $$ \begin{gathered} \frac{\partial \overline{\rho }{\tilde{ u}}_{i}}{\partial t}+\frac{\partial }{\partial {x}_{j}}\left(\overline{\rho }{\tilde{ u}}_{i}{\tilde{ u}}_{j}\right)=\\ -\frac{\partial \overline{p}}{\partial {x}_{i}}+\frac{\partial }{\partial {x}_{j}}\left[{\mu }_{\mathrm{e}\mathrm{f}\mathrm{f}}\left(\frac{\partial {\tilde{ u}}_{i}}{\partial {x}_{j}}+\frac{\partial {\tilde{ u}}_{j}}{\partial {x}_{i}}-\frac{2}{3}\frac{\partial {\tilde{ u}}_{k}}{\partial {x}_{k}}{\delta }_{ij}\right)\right]+\overline{\rho }{\mathrm{g}}_{i} \end{gathered} $$ (2) $$ \begin{gathered} \frac{\partial }{{\partial t}}\left( {\bar \rho \tilde E} \right) + \frac{\partial }{{\partial {x_j}}}\left( {{{\tilde u}_j}\left( {\bar \rho \tilde E + \bar p} \right)} \right) = \\ \frac{\partial }{{\partial {x_j}}}\left[ {\frac{{{\mu _{{\rm{eff}}}}{c_{\rm{p}}}}}{{{P_{r{\rm{eff}}}}}}\frac{{\partial \tilde T}}{{\partial {x_j}}} - \sum\limits_{{m}} {{{\tilde h}_{{m}}}} \left( { - \frac{{{\mu _{{\rm{eff}}}}}}{{{S_{{\rm{c}},{\rm{eff}}}}}}\frac{{\partial {{\tilde Y}_{{m}}}}}{{\partial {x_j}}}} \right)} \right. + \\ \left.{{\tilde u}_i} {\mu_{{\rm{eff}}}\left( {\frac{{\partial {{\tilde u}_i}}}{{\partial {x_j}}} + \frac{{\partial {{\tilde u}_j}}}{{\partial {x_i}}} - \frac{2}{3}\frac{{\partial {{\tilde u}_k}}}{{\partial {x_k}}}{\delta _{ij}}} \right)} \right] + \sum\limits_{{m}} {{R_{{m}}}} {H_{\rm{c}}} \end{gathered}$$ (3) $$ \frac{\partial }{\partial t}\left(\overline{\rho }{\tilde{Y}}_{{m}}\right)+\frac{\partial }{\partial {x}_{j}}\left({\overline{\rho }\tilde{u}}_{j}{\tilde{Y}}_{{m}}\right)=\frac{\partial }{\partial {x}_{j}}\left(\frac{{\mu }_{\mathrm{e}\mathrm{f}\mathrm{f}}}{S_{{\rm{c}},\mathrm{e}\mathrm{f}\mathrm{f}}}\frac{\partial {\tilde{Y}}_{{m}}}{\partial {x}_{j}}\right)+{R}_{{m}} $$ (4) 式中：ρ为密度；x为空间坐标；u为速度分量；g为重力加速度；E为总能量；cp为混合物在恒压下的比热容；Pr为普朗特数；T为温度；h为焓；Sc为进度变量守恒方程中的源项；Y为质量分数；R为火焰前锋平均半径；H为燃烧热。上标和下标中，上标：−为LES过滤后的量，~为LES质量加权过滤后的量；下标：i，j，k为空间坐标索引；m为气体混合的第m个组分；c为源项在进行变量方程；eff为有效值。
k方程和ω方程分别是：
$$ \frac{\partial \rho k}{\partial t}+\frac{\partial \rho {u}_{j}k}{\partial {x}_{j}}-\frac{\partial }{\partial {x}_{j}}\left[\left(\mu +{\sigma }_{k}{\mu }_{{\rm{T}}}\right)\frac{\partial k}{\partial {x}_{j}})\right]={\tau }_{ij}\frac{\partial {u}_{i}}{\partial {x}_{i}}-{\beta }^{*}\rho k\omega $$ (5) $$\begin{split} &\frac{\partial \rho \omega }{\partial t}+\frac{\partial \rho {u}_{j}\omega }{\partial {x}_{j}}-\frac{\partial }{\partial {x}_{j}}\left[\left(\mu +{\sigma }_{\omega }{\mu }_{{\rm{T}}}\right)\frac{\partial k}{\partial {x}_{j}})\right]=\\ &\frac{\gamma }{{v}_{{\rm{T}}}}{\tau }_{ij}\frac{\partial {u}_{i}}{\partial {x}_{i}}-\beta {\omega }^{2}+2\rho \left(1-{F}_{1}\right)\frac{{{\sigma }_{\omega} ^2}}{\omega }\frac{\partial k}{\partial {x}_{j}}\frac{\partial \omega }{\partial {x}_{j}} \end{split}$$ (6) 其中：k为湍流动能；ω为湍流动能的特定消散，它们分别定义为
$$ k=\frac{1}{2}{u}_{i}^{'}{u}_{j}^{'} $$ (7) $$ \omega =\frac{\varepsilon }{k{\beta }^{*}} $$ (8) 其中，ε为湍流动能消散率：
$$ \varepsilon =v\frac{\partial {u}_{i}^{'}}{\partial {x}_{k}}\frac{\partial {u}_{i}^{'}}{\partial {x}_{k}} $$ (9) 其中，τij为雷诺应力：
$$ {\tau }_{ij}={\mu }_{t}\left(\frac{\partial {u}_{i}}{\partial {x}_{j}}+\frac{\partial {u}_{j}}{\partial {x}_{i}}\right)-\frac{2}{3}\rho {\delta }_{ij} $$ (10) 式中：t为时间；p为动力压强；μt为涡动黏性；δij为克罗内克函数；vT=μt/ρ为动力涡黏性；σk，β*，σω，γ，β，σω2为模型系数；F1为混合函数[27]。
组分输运共涉及3种气体，分别为H2、CO和CH4，气体密度均为理想气体，时间属性选择定常流动（稳态）。本次CFD数值模拟计算域的边界条件如图1所示。
本次CFD数值模拟网格独立性检验共选择3种网格尺寸，其网格数量分别为101850（high）、71604（medium）、46350（low），泄漏气体为H2，泄漏口直径为3 mm，泄漏压力为1 MPa，气体温度为300 K。图2表示的是中心轴线H2浓度分布曲线，可见选用中质量网格（medium）得到的结果与高质量网格（high）十分接近。因此在考虑时间成本并保证结果准确的前提下，本次CFD数值模拟采用中质量（medium）网格开展相关研究。
![]()
图 2 3组网格中心轴线处H2浓度分布Figure 2. H2 concentration distribution at the central axis of the three groups1.2 不同条件下可燃气体泄漏扩散规律及泄漏危险范围
如前文所述，CFD数值模拟共涉及3种气体，分别为：H2、CO和CH4。对于3种气体，其危险浓度即为爆炸极限，H2：4%~75%，CO：12.4%~75%，CH4：4.9%~16%。
CFD数值模拟共设置81组工况，泄漏气体分别为H2、CO和CH4，泄漏口直径分别为2、3、4 mm，合成油品中低温F-T合成所需压力2~3 MPa[3]，温度280 ℃以下，因此工况中，泄漏压力选定为1、2、3 MPa，气体泄漏温度为300、400、500 K。
1.2.1 不同泄漏压力下可燃气体泄漏扩散规律及泄漏危险范围
通过对H2、CO和CH4三种气体在不同泄漏压力下气体泄漏浓度分布云图进行对比（图3），可以发现：
![]()
图 3 不同泄漏压力下可燃气体扩散规律及泄漏危险范围Figure 3. Diffusion mechanism and danger range of combustible gas leakage under different leakage pressure1）3种气体在温度、直径等其他因素不变的条件下，随着泄漏压力的增大，气体泄漏危险范围随之明显增加，其危险范围除轴向长度增长之外，宽度也有所增大。
2）对于不同气体，压力同为1 MPa时，H2泄漏的危险范围长度与宽度分别为2.36 m和0.316 m，CO与CH4则仅为0.23、0.028 m和0.67、0.09 m，可见H2的泄漏危险范围最大；当泄漏压力从1 MPa上升到2 MPa时，H2泄漏危险范围长度和宽度分别增加了11%和13%，而CO泄漏危险范围长、宽分别增加了17%和11%，CH4则分别增加了9%和7%，可见H2、CO和CH4三种气体泄漏危险范围均受到气体温度的显著影响。
从图4—图6可以看出，泄漏压力越大，中心轴线浓度越高，如H2在泄漏压力1 MPa，泄漏口直径4 mm，泄漏温度300 K时，距泄漏口0.3 m处，体积分数为33%，而当泄漏压力为2 MPa时，同样位置的H2体积分数则升高到38%。
![]()
图 4 不同泄漏压力下中心轴线H2浓度分布曲线（4 mm-300 K）Figure 4. H2 concentration distribution of center axis under different leakage pressure (4 mm-300 K)![]()
图 5 不同泄漏压力下中心轴线CO浓度分布曲线（4 mm-300 K）Figure 5. CO concentration distribution of center axis under different leakage pressure (4 mm-300 K)![]()
图 6 不同泄漏压力下中心轴线CH4浓度分布曲线（4 mm-300 K）Figure 6. CH4 concentration distribution of center axis under different leakage pressure (4 mm-300 K)1.2.2 不同气体温度下可燃气体泄漏扩散规律及泄漏危险范围
通过对H2、CO和CH4三种气体在不同泄漏温度下气体泄漏浓度分布云图进行对比（图7），可以发现：
![]()
图 7 不同泄漏温度下可燃气体扩散规律及泄漏危险范围Figure 7. Diffusion mechanism and danger range of combustible gas leakage under different leakage temperature1）3种气体在泄漏压力、直径等其他因素不变的条件下，随着气体温度的增大，气体泄漏危险范围随之明显降低，其危险范围除轴向长度减小之外，宽度也有所减小。
2）对于不同气体，温度同为300 K时，H2泄漏的危险范围长度与宽度分别为2.04 m和0.28 m，CO与CH4则仅为0.24，0.024和0.62，0.074 m，可见H2的泄漏危险范围最大；当泄漏温度从300 K上升到400 K时，H2泄漏危险范围长度和宽度分别降低了11%和12%，而CO泄漏危险范围长、宽分别降低了4.2%和12.5%，CH4则分别降低了10%和14%。可见H2、CO和CH4三种气体泄漏危险范围均受到气体温度的显著影响。
从图8—图10可以看出，气体温度越高，中心轴线浓度越低，如H2在泄漏温度为300 K，泄漏口直径3 mm，泄漏压力3 MPa时，距泄漏口0.3 m处，体积分数为37%，而当泄漏温度变为400 K时，同样位置的H2体积分数为34%。从图8—图10可以看出，气体温度越高，中心轴线浓度越低，如H2在泄漏温度为300 K，泄漏口直径3 mm，泄漏压力3 MPa时，距泄漏口0.3 m处，体积分数为37%，而当泄漏温度变为400 K时，同样位置的H2体积分数为34%。
![]()
图 8 不同气体温度下中心轴线H2浓度分布曲线（3 mm-3 MPa）Figure 8. H2 concentration distribution of center axis under different leakage temperature (3 mm-3 MPa)![]()
图 9 不同气体温度下中心轴线CO浓度分布曲线（3 mm-3 MPa）Figure 9. CO concentration distribution of center axis under different leakage temperature (3 mm-3 MPa)![]()
图 10 不同气体温度下中心轴线CH4浓度分布曲线（3 mm-3 MPa）Figure 10. CH4 concentration distribution of center axis under different leakage temperature (3 mm-3 MPa)1.2.3 不同泄漏口直径下可燃气体泄漏扩散规律及泄漏危险范围
通过对H2、CO和CH4三种气体在不同泄漏口直径下气体泄漏浓度分布云图进行对比（图11），可以发现：
![]()
图 11 不同泄漏口直径下可燃气体泄漏扩散规律及泄漏危险范围Figure 11. Diffusion mechanism and danger range of combustible gas leakage under different diameter of leak hole1）3种气体在温度、泄漏压力等其他因素不变的条件下，随着泄漏口直径的增大，气体泄漏危险范围随之明显增加，其危险范围除轴向长度增长之外，宽度也有所增大。
2）对于不同气体，直径为2 mm时，H2泄漏的危险范围长度与宽度分别为1.11 m和0.148 m，CO与CH4则仅为0.14 m，0.013 m和0.62 m，0.037 m，可见H2的泄漏危险范围最大；当泄漏口直径从2 mm上升到3 mm时，H2泄漏危险范围长度和宽度分别增加了49%和54%，而CO泄漏危险范围长、宽分别增加了57%和46%，CH4则分别增加了50%和54%。可见H2、CO和CH4三种气体泄漏危险范围均受到泄漏口直径的显著影响。
从图12—图14可以看出，泄漏口直径越大，中心轴线浓度越高，如H2在泄漏口2 mm，泄漏温度500 K，泄漏压力3 MPa时，距泄漏口0.3 m处，浓度为19%，而当泄漏口直径变为3 mm时，同样位置的H2浓度为33%。
![]()
图 12 不同泄漏直径下中心轴线H2浓度分布曲线（500 K-3 MPa）Figure 12. H2 concentration distribution of center axis under different diameter of leak hole (500 K-3 MPa)![]()
图 13 不同泄漏直径下中心轴线CO浓度分布曲线（500 K-3 MPa）Figure 13. CO concentration distribution of center axis under different diameter of leak hole (500 K-3 MPa)![]()
图 14 不同泄漏直径下中心轴线CH4浓度分布曲线（500 K-3 MPa）Figure 14. CH4 concentration distribution of center axis under different diameter of leak hole (500 K-3 MPa)2. 可燃气体泄漏危险范围预测模型
本次CFD数值模拟共考虑3种气体，分别为H2，CO和CH4，并充分考虑泄漏口大小、泄漏压力和泄漏温度后，根据得到的危险范围，对每种泄漏气体危险范围的长和宽分别进行拟合，得出泄漏危险范围定量预测经验公式（图15—图20）。对参量进行无量纲处理，得到的泄漏危险范围定量预测模型为：
![]()
图 15 H2泄漏危险范围长度经验公式Figure 15. Empirical formula of the length of H2 leakage danger range![]()
图 16 H2泄漏危险范围宽度经验公式Figure 16. Empirical formula of the width of H2 leakage danger range![]()
图 17 CO泄漏危险范围长度经验公式Figure 17. Empirical formula of the length of CO leakage danger range![]()
图 18 CO泄漏危险范围宽度经验公式Figure 18. Empirical formula of the width of CO leakage danger range![]()
图 19 CH4泄漏危险范围长度经验公式Figure 19. Empirical formula of the length of CH4 leakage danger range![]()
图 20 CH4泄漏危险范围宽度经验公式Figure 20. Empirical formula of the width of CH4 leakage danger range对于H2：
$$ 长度: L/d=428\times (P_{0}/P_{\infty })^{0.14}(T_{\infty }/T_{0})^{0.43}-0.38 $$ (11) $$ 宽度: W/d = 60.82 \times (P_{0}/P_{\infty })^{0.13}(T_{\infty }/T_{0})^{0.44} - 0.83 $$ (12) 对于CO：
$$ 长度: L/d=44.85\times (P_{0}/P_{\infty })^{0.16}(T_{\infty }/T_{0})^{0.23}-0.69 $$ (13) $$ 宽度: W/d=5.03\times(P_{0}/P_{\infty })^{0.14}(T_{\infty }/T_{0})^{0.43}+0.10 $$ (14) 对于CH4：
$$ 长度: L/d=133.63\times (P_{0}/P_{\infty })^{0.13}(T_{\infty }/T_{0})^{0.40}-0.35 $$ (15) $$ 宽度: W/d=18.05\times (P_{0}/P_{\infty })^{0.10}(T_{\infty }/T_{0})^{0.53}-0.13 $$ (16) 其中：L为泄漏危险范围长度，m；W为泄漏危险范围长度，m；d为泄漏口直径，m；P0为泄漏压力，Pa；P∞为环境压力（101325 Pa）；T0为泄漏气体温度，K，T∞为环境温度，300 K。
3. 结 论
1）3种气体随着泄漏压力的增大，气体泄漏危险范围随之明显增加，且泄漏压力越大，泄漏口中心轴线气体浓度越高。
2）随着泄漏气体温度的增大，气体泄漏危险范围随之明显降低，危险范围除轴向长度及宽度均有所减小，泄漏气体温度越高，泄漏口中心轴线气体浓度越低。
3）随着泄漏口直径的增大，气体泄漏危险范围随之明显增加，危险范围轴向长度和宽度均有所增大，且泄漏口尺寸对于可燃气体泄漏危险范围的影响大于泄漏压力和泄漏温度。泄漏口直径越大，泄漏口中心轴线气体浓度越高。
根据数值模拟得到的危险范围数据，考虑气体、泄漏口大小、泄漏压力和泄漏温度后，对每种泄漏气体危险范围的长和宽分别进行拟合，得出泄漏危险范围预测经验公式。
4）通过可燃气体泄漏扩散规律和危险范围分析结果，可为区域救援协同响应和应急资源快速调配提供技术支持，亦可为重大工业（危险化学品）事故救援或城市公共安全事件处置应急指挥辅助决策提供参考依据。
-
![]()
图 2 3组网格中心轴线处H2浓度分布
Figure 2. H2 concentration distribution at the central axis of the three groups
![]()
图 3 不同泄漏压力下可燃气体扩散规律及泄漏危险范围
Figure 3. Diffusion mechanism and danger range of combustible gas leakage under different leakage pressure
![]()
图 4 不同泄漏压力下中心轴线H2浓度分布曲线（4 mm-300 K）
Figure 4. H2 concentration distribution of center axis under different leakage pressure (4 mm-300 K)
![]()
图 5 不同泄漏压力下中心轴线CO浓度分布曲线（4 mm-300 K）
Figure 5. CO concentration distribution of center axis under different leakage pressure (4 mm-300 K)
![]()
图 6 不同泄漏压力下中心轴线CH4浓度分布曲线（4 mm-300 K）
Figure 6. CH4 concentration distribution of center axis under different leakage pressure (4 mm-300 K)
![]()
图 7 不同泄漏温度下可燃气体扩散规律及泄漏危险范围
Figure 7. Diffusion mechanism and danger range of combustible gas leakage under different leakage temperature
![]()
图 8 不同气体温度下中心轴线H2浓度分布曲线（3 mm-3 MPa）
Figure 8. H2 concentration distribution of center axis under different leakage temperature (3 mm-3 MPa)
![]()
图 9 不同气体温度下中心轴线CO浓度分布曲线（3 mm-3 MPa）
Figure 9. CO concentration distribution of center axis under different leakage temperature (3 mm-3 MPa)
![]()
图 10 不同气体温度下中心轴线CH4浓度分布曲线（3 mm-3 MPa）
Figure 10. CH4 concentration distribution of center axis under different leakage temperature (3 mm-3 MPa)
![]()
图 11 不同泄漏口直径下可燃气体泄漏扩散规律及泄漏危险范围
Figure 11. Diffusion mechanism and danger range of combustible gas leakage under different diameter of leak hole
![]()
图 12 不同泄漏直径下中心轴线H2浓度分布曲线（500 K-3 MPa）
Figure 12. H2 concentration distribution of center axis under different diameter of leak hole (500 K-3 MPa)
![]()
图 13 不同泄漏直径下中心轴线CO浓度分布曲线（500 K-3 MPa）
Figure 13. CO concentration distribution of center axis under different diameter of leak hole (500 K-3 MPa)
![]()
图 14 不同泄漏直径下中心轴线CH4浓度分布曲线（500 K-3 MPa）
Figure 14. CH4 concentration distribution of center axis under different diameter of leak hole (500 K-3 MPa)
![]()
图 15 H2泄漏危险范围长度经验公式
Figure 15. Empirical formula of the length of H2 leakage danger range
![]()
图 16 H2泄漏危险范围宽度经验公式
Figure 16. Empirical formula of the width of H2 leakage danger range
![]()
图 17 CO泄漏危险范围长度经验公式
Figure 17. Empirical formula of the length of CO leakage danger range
![]()
图 18 CO泄漏危险范围宽度经验公式
Figure 18. Empirical formula of the width of CO leakage danger range
![]()
图 19 CH4泄漏危险范围长度经验公式
Figure 19. Empirical formula of the length of CH4 leakage danger range
-
[1] 赵开功,张晓蕾,李 曼,等. 基于一体化运营综合能源企业智能应急救援技术研究[J]. 中国安全生产科学技术,2022,18(5):235−240. ZHAO Kaigong,ZHANG Xiaolei,LI Man,et al. Research on intelligent emergency rescue technology of comprehensive energy enterprise based on integrated operation[J]. Journal of Safety Science and Technology,2022,18(5):235−240.
[2] ZHAO Kaigong, WANG Haiyan, ZHENG Dengfeng. Research on structural similarity design emergency exercise’s scenario [J]. Geofluids, 2022: 6590957.
[3] 赵开功,方 衡. 一体化运营综合能源集团智能应急指挥平台的研究与实现[J]. 矿业安全与环保,2021,48(6):58−64,71. ZHAO Kaigong,FANG Heng. Research and implementation of integrated operation integrated energy group intelligent emergency command platform[J]. Mining Safety & Environmental Protection,2021,48(6):58−64,71.
[4] 步学朋,俞珠峰,王 强,等. 煤制油化工建设项目技术经济指标与参数研究[J]. 煤炭科学技术,2020,48(1):120−124. BU Xuepeng,YU Zhufeng,WANG Qiang,et al. Research on technical-economic indexes and parameters of coal to liquid and chemical construction project[J]. Coal Science and Technology,2020,48(1):120−124.
[5] 唐宏青. 我国煤制油技术的现状和发展[J]. 化学工程,2010,38(10):1−8. doi: 10.3969/j.issn.1005-9954.2010.10.001 TANG Hongqing. Status and development of coal to oil technique in China[J]. Chemical Engineering,2010,38(10):1−8. doi: 10.3969/j.issn.1005-9954.2010.10.001
[6] 徐振刚. 化工合成用煤气化工艺的选择[C]//中国煤炭学会煤化专业委员会年会暨新型煤化工技术研讨会会议文集, 2004: 13–18. [7] 刘彦青. 动态回采工作面煤壁瓦斯涌出数值计算及现场应用[J]. 煤炭科学技术,2021,49(6):195−204. doi: 10.13199/j.cnki.cst.2021.06.023 LIU Yanqing. Numerical calculation and field application of gas emission from coal wall in dynamic stoping face[J]. Coal Science and Technology,2021,49(6):195−204. doi: 10.13199/j.cnki.cst.2021.06.023
[8] 李 伟, 张 奇. 高压气体泄漏喷射和扩散的数值模拟[C]//中国职业安全健康协会2008年学术年会论文集, 2008: 319–325. [9] 傅 杰. 球罐内气体泄漏扩散特性模拟及可燃区域预测[D]. 大连: 大连理工大学, 2014. FU Jie. Numerical simulation and combustible region prediction of gas leakage and diffusion in spherical tank[D]. Dalian : Dalian University of Technology, 2014.
[10] 刘延雷,徐 平,郑津洋,等. 管道输运高压氢气与天然气的泄漏扩散数值模拟[J]. 太阳能学报,2008,29(10):1252−1255. doi: 10.3321/j.issn:0254-0096.2008.10.013 LIU Yanlei,XU Ping,ZHENG Jinyang,et al. Numerical simulation on the dispersion of hydrogen and natural gas due to high pressured pipeline leakage[J]. Acta Energiae Solaris Sinica,2008,29(10):1252−1255. doi: 10.3321/j.issn:0254-0096.2008.10.013
[11] 周伟国,刘东京,滕卯寅. 高压天然气管道泄漏扩散CFD数值模拟[J]. 管道技术与设备,2015(5):5−8. ZHOU Weiguo,LIU Dongjing,TENG Maoyin. CFD simulation on natural gas leakage and diffusion of high-pressure pipeline[J]. Pipeline Technique and Equipment,2015(5):5−8.
[12] 高 嵩,王茂廷. 高压天然气管道泄漏及事故危害规律研究[J]. 当代化工,2017,46(1):82−85. GAO Song,WANG Maoting. Study on leakage accident laws and hazards of high pressure natural gas pipeline[J]. Contemporary Chemical Industry,2017,46(1):82−85.
[13] 陈 坤, 林 浩, 陈 洁, 等. 通风状态下综合管廊燃气管道小孔泄漏扩散模拟研究[J/OL]. 安全与环境学报: 1-11[2023-04-19]. DOI: 10.13637/j.issn.1009-6094.2022.1480. CHEN Kun, LIN Hao, CHEN Jie, et al. Simulation research on small hole leakage and diffusion of gas pipeline in utility tunnel under ventilation condition [J/OL]. Journal of Safety and Environment: 1-11 [2023-04-19]. DOI: 10.13637/j.issn.1009-6094.2022.1480.
[14] 刘敏鸿,刘成良,翁韦强,等. 居民用户天然气泄漏三维数值模拟分析[J]. 煤气与热力,2022,42(8):24−28. doi: 10.13608/j.cnki.1000-4416.2022.08.021 LIU Minhong,LIU Chengliang,WONG Weiqiang,et al. Three-dimensional numerical simulation analysis of natural gas leakage in residential users[J]. Gas & Heat,2022,42(8):24−28. doi: 10.13608/j.cnki.1000-4416.2022.08.021
[15] 袁裕鹏,崔伟逸,沈辉,等. 燃料电池船舶舱内氢气泄漏扩散的数值模拟[J]. 交通运输工程学报,2022,22(4):196−209. YUAN Yupeng,CUI Weiyi,SHEN Hui,et al. Numerical simulation of leakage and diffusion of hydrogen in cabin of fuel cell ship[J]. Journal of Traffic and Transportation Engineering,2022,22(4):196−209.
[16] 汪建平,段卫东,李春青,等. 厨房内天然气泄漏扩散过程数值模拟[J]. 安全与环境学报,2022,22(1):298−305. doi: 10.13637/j.issn.1009-6094.2020.1710 WANG Jianping,DUAN Weidong,LI Chunqing,et al. Study on numerical simulation of indoor natural gas leakage and diffusion process[J]. Journal of Safety and Environment,2022,22(1):298−305. doi: 10.13637/j.issn.1009-6094.2020.1710
[17] 万宇飞,李立婉,王文光,等. 基于Fluent的风向对LPG储罐泄漏扩散规律研究[J]. 石油与天然气化工,2021,50(5):98−103. doi: 10.3969/j.issn.1007-3426.2021.05.019 WAN Yufei,LI Liwan,WANG Wenguang,et al. Influence of wind direction on LPG tank leakage rule with FLUENT software[J]. Chemical Engineering of Oil & Gas,2021,50(5):98−103. doi: 10.3969/j.issn.1007-3426.2021.05.019
[18] 朱 静,綦远磊,辛培刚,等. 多因素下居民室内天然气泄漏扩散数值模拟[J]. 消防科学与技术,2021,40(7):1013−1016. doi: 10.3969/j.issn.1009-0029.2021.07.014 ZHU Jing,QI Yuanlei,XIN Peigang,et al. Numerical simulation of natural gas leakage and diffusion in residence space under multiple factors[J]. Fire Science and Technology,2021,40(7):1013−1016. doi: 10.3969/j.issn.1009-0029.2021.07.014
[19] 喻健良,刘长远,闫兴清,等. 高压超临界CO2管道泄漏扩散灾害模拟[J]. 安全与环境学报,2022,22(3):1495−1502. YU Jianliang,LIU Changyuan,YAN Xingqing,et al. Simulation research of leakage and diffusion disaster in high-pressure supercritical CO2 pipelines[J]. Journal of Safety and Environment,2022,22(3):1495−1502.
[20] 周 宁,陈 力,吕孝飞,等. 环境温度对LNG泄漏扩散影响的数值模拟[J]. 油气储运,2021,40(3):352−360. ZHOU Ning,CHEN Li,LYV Xiaofei,et al. Numerical simulation on the effect of ambient temperature on leakage and diffusion of LNG[J]. Oil & Gas Storage and Transportation,2021,40(3):352−360.
[21] 程方明,张安邦,王 焘,等. 高压天然气非恒定速率泄漏扩散数值模拟研究[J]. 中国安全生产科学技术,2021,17(1):90−95. CHENG Fangming,ZHANG Anbang,WANG Tao,et al. Numerical simulation study on non-constant rate leakage and diffusion of high-pressure natural gas[J]. Journal of Safety Science and Technology,2021,17(1):90−95.
[22] 陶婧莹,李兆慈,郭志超等. LNG泄漏气体浓度和温度扩散过程计算[J]. 低温与超导,2020,48(6):12−19,25. TAO Jingying,LI Zhaoci,GUO Zhichao,et al. Gas concentration and temperature diffusion numerical simulation for LNG leakage[J]. Cryogenics & Superconductivity,2020,48(6):12−19,25.
[23] 周 宁,王伟天,陈 兵,等. 石化管道丁烷气体泄漏扩散数值模拟[J]. 安全与环境工程,2020,27(2):175−182. ZHOU Ning,WANG Weitian,CHEN Bing,et al. Numerical simulation of butane gas leakage and diffusion in petrochemical pipelines[J]. Safety and Environmental Engineering,2020,27(2):175−182.
[24] 张晓蕾,张 洁. 化工厂天然气管道泄漏火灾研究[J]. 能源化工,2021,42(1):66−71. ZHANG Xiaolei,ZHANG Jie. Study on fire accident caused by the leakage of the natural gas pipeline[J]. Energy Chemical Industry,2021,42(1):66−71.
[25] 钱炜祺, 周 宇, 邓有奇, 等. k-ωSST两方程湍流模型参数辨识初步研究[C]//中国力学学会学术大会, 2009: 253. [26] 弓 亮. 管道内高压氢气泄漏自燃机理实验与数值模拟研究[D]. 合肥: 中国科学技术大学, 2019. GONG Liang. Experimental and numerical study on the mechanism of spontaneous ignition during high-pressure hydrogen release into a tube [D]. Hefei: University of Science and Technology of China, 2019.
[27] 闫晓惠,陈 新,李华煜. 含收缩段明渠流的三维k-omega SST模拟[J]. 水利科技与经济,2015,21(11):12−14. YAN Xiaohui,CHEN Xin,LI Huayu. 3D k-omega SST modeling of open channel flow with contraction[J]. Water Conservancy Science and Technology and Economy,2015,21(11):12−14.
-
期刊类型引用(7)
1. 张晓蕾，高进东，赵开功，李严肃，李长明，谢子彬. 煤矿事故智能应急预案生成方法研究. 矿业安全与环保. 2024(01): 78-85+91 . 
百度学术
2. 张晓蕾，赵开功，李长明，高进东，王睿迪，张洁. 基于GIS的煤制油可燃气泄漏事故分析与应急处置. 中国安全生产科学技术. 2024(05): 138-145 . 百度学术
3. 冯占科，杨俊彩，李长明，张晓蕾. 煤矿采掘动态接续平衡预警系统设计及应用. 矿业安全与环保. 2024(03): 166-171 . 百度学术
4. 马坡. 危险气体泄漏监测在工业生产安全中的实践与优化——以甲烷为例. 化学工程与装备. 2024(08): 160-162 . 百度学术
5. 刘江，赵开功，张晓蕾，闫力维，李长明，王睿迪. 基于知识图谱的选煤法规标准培训平台及应用. 中国安全科学学报. 2024(S1): 226-233 . 百度学术
6. 柳璇，高琳，齐文娇，王锦龙，朱宾飞. 液态二氧化碳压力容器泄漏扩散非稳态特性研究. 山东化工. 2024(24): 203-206 . 百度学术
7. 张晓蕾，赵开功，李长明，陈刚. 城市应急队伍救援能力评估与应用研究. 消防科学与技术. 2023(12): 1724-1728 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 129
- HTML全文浏览量: 16
- PDF下载量: 40
- 被引次数: 10
