高级检索

典型可燃气体小孔泄漏危险范围预测模型研究

赵开功, 张晓蕾, 李长明, 陈刚, 盖泳伶

赵开功,张晓蕾,李长明,等. 典型可燃气体小孔泄漏危险范围预测模型研究[J]. 煤炭科学技术,2023,51(3):281−290

. DOI: 10.13199/j.cnki.cst.2022-1948
引用本文:

赵开功,张晓蕾,李长明,等. 典型可燃气体小孔泄漏危险范围预测模型研究[J]. 煤炭科学技术,2023,51(3):281−290

. DOI: 10.13199/j.cnki.cst.2022-1948

ZHAO Kaigong,ZHANG Xiaolei,LI Zhangming,et al. Numerical simulation on prediction model of risk range of typical gas release through small holes[J]. Coal Science and Technology,2023,51(3):281−290

. DOI: 10.13199/j.cnki.cst.2022-1948
Citation:

ZHAO Kaigong,ZHANG Xiaolei,LI Zhangming,et al. Numerical simulation on prediction model of risk range of typical gas release through small holes[J]. Coal Science and Technology,2023,51(3):281−290

. DOI: 10.13199/j.cnki.cst.2022-1948

典型可燃气体小孔泄漏危险范围预测模型研究

基金项目: 

国家自然科学基金资助项目(52074156);国家重点研发计划资助项目(2016YFC0801805);国家能源集团基金资助项目(10032998IC190047)

详细信息
    作者简介:

    赵开功: (1981—) ,男,山东济宁人,高级工程师,博士研究生。Tel:010-58131536,E-mail:smilevel@163.com

    通讯作者:

    张晓蕾: (1984—) ,女,山西大同人,高级工程师,博士研究生。Tel:010-64954704,E-mail:zhangxl@chinasafety.ac.cn

  • 中图分类号: TE973

Numerical simulation on prediction model of risk range of typical gas release through small holes

Funds: 

National Natural Science Foundation of China (52074156); National Key Research and Development Program of China (2016YFC0801805); Technology Project of China Energy (10032998IC190047)

  • 摘要:

    为了缓解我国石油资源紧张的现状,发展合成油工业作为我国石油资源的补充是合理的途径,而在化工合成中,大多数反应都需要在加压环境中进行。泄漏是最常见的事故之一,而现阶段对于加压装置、管道气体泄漏的研究还有待完善。因此,通过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.

  • 液压支架支护是巷道冲击地压防治重要环节,随着我国冲击地压发生机理和理论研究的不断深入和完善,巷道液压支架在冲击地压发生启动和破坏过程所起的作用逐步得到重视[1-3]。2020年7月1日,国家煤监局颁布实施的《煤矿安全生产标准化管理体系基本要求及评分方法(试行)》中首次提出“冲击地压矿井使用工作面端头支架、两巷超前支护液压支架和吸能装置”,并作为煤矿采煤标准化评分项[4]。在传统的巷道液压支护装备上安装吸能装置,是近几年我国冲击地压矿井支护的重大进步,巷道防冲支架可以有效提升巷道防冲能力[5-7]。吸能液压支架立柱是巷道防冲支架的关键承载部件,支架的抗冲性能主要通过吸能液压支架立柱实现[8-9],开展冲击载荷下巷道吸能液压支架立柱防冲性能研究对于支架参数及结构设计具有重要的理论和实践意义。

    近年来,学者们在巷道支护方面进行了大量研究。潘一山等[2]提出了冲击地压巷道三级支护理论、冲击地压矿井巷道支护理论[10],设计出的防冲液压支架在义马耿村矿安装应用,经受住了ML3.0级冲击地压,巷道支护完好无人员伤亡[11]。唐治等[12]设计出了自移式吸能防冲巷道超前支架。刘欣科等[13]对液压立柱在冲击载荷作用下动态响应进行了数值模拟研究,赵忠辉等[14]分析了立柱的应力、应变分布,研究了立柱的抗冲击性能。安栋等[15]探究了防冲吸能构件的力学特性及其应用效果,肖永惠等[16]对支架吸能构件屈曲吸能可靠性进行了研究,许海亮等[17]依据塑性变形区域最大化原则提出了一种新型矿用吸能构件。学者们在巷道支架支护和吸能防冲领域取得了许多有价值的研究成果,但现有研究大多集中于吸能构件设计和单一立柱的防冲能力,吸能液压支架立柱整体的防冲性能、吸能装置与液压系统的整体受冲过程等方面还需要进一步研究。

    笔者采用有限元与光滑粒子流体动力学方法建立了吸能液压支架立柱流固耦合、流体大变形模型,实现了立柱在冲击载荷作用下的数值模拟。根据防冲支护设计六项原则,对比了普通液压支架立柱与吸能液压支架立柱的防冲性能,为后续吸能液压支架立柱优化设计、防冲支架的改进提供参考。

    冲击地压多发生于巷道中[18],提高支护应力可显著提高冲击地压发生的临界荷载[19-20],液压支架支护被广泛应用于巷道冲击地压的防治。但由于冲击地压的复杂性,某些情况下冲击地压仍会不可避免的发生。传统液压支架受冲时普遍存在抗不住、让不及、吸不了的问题,支架受冲破坏照片如图1所示,难以在冲后形成有效防护。

    图  1  支架受冲破坏照片
    Figure  1.  Photo of support damage

    潘一山等[11]根据冲击地压扰动响应失稳理论、巷道防冲支护理论,提出了防冲支护设计6项原则,从让位位移、让位阻力、让位速度、让位刚度、让位频率、让位能量6个方面设计了防冲装置。该装置是一种带有预折纹的薄壁结构,由高强钢采用固定模具压制而成,是液压立柱防冲吸能的核心构件[8, 15],装置外观及尺寸如图2所示。

    图  2  吸能装置外观及尺寸
    Figure  2.  Appearance and dimensions of the energy-absorbing device

    某单元式巷道吸能液压支架装载的FCZ1650/25/42型吸能液压立柱结构简化模型如图3所示,包括液压活塞、高压缸、低压缸、吸能套筒4个部分,缸内填充乳化液,吸能套筒装载吸能装置。当围岩发生剧烈冲击时,吸能装置达到启动荷载后将高速吸能让位,使液压缸体在引导下发生位移,保护柱体各部分不受损坏。该立柱的核心参数见表1,设计满足国家标准GB25974.2—2010《液压支架立柱技术条件》[21],本课题组研制的不同结构与防冲性能的巷道防冲液压支架如图4所示,目前已在耿村煤矿、龙家堡煤矿、辽宁老虎台等矿开展应用。

    图  3  吸能防冲液压支架立柱结构
    Figure  3.  Model of energy-absorbing hydraulic bracing column
    表  1  吸能液压支架立柱核心参数
    Table  1.  Energy-absorbing hydraulic bracing column core parameters
    (a)液压立柱参数
    缸壁内径/
    mm
    缸壁壁厚/
    mm
    初撑力 工作阻力
    液压/MPa 反力/kN 液压/MPa 反力/kN
    ø230/ø180 21.5/17.5 31.5 1309 39.5 1650
    (b)吸能装置参数
    让位阻力/
    kN
    材料强度/
    MPa
    壁厚/
    mm
    设计让位位移/
    mm
    2500 890 8 130
    下载: 导出CSV 
    | 显示表格
    图  4  吸能防冲液压支架架型
    Figure  4.  Suction anti-impact hydraulic support frame type

    液压缸内充满高压乳化液,乳化液与液压缸壁接触部位对模拟精度有较高要求,光滑粒子流体动力学方法(SPH)可用于处理具有高应力、大变形特点的乳化液[22-24]。首先对乳化液基本属性进行定义,乳化液密度设置为1.08 t/mm3,选用US-UP状态方程,取c0=1400 m/s,s=0,Gamma=0。使用有限元模拟软件提供的无网格粒子化模拟方法,在单元控制属性界面将其转化为粒子,设置为三次(Cubic)样条曲线差值,取PPD=1,实现乳化液的粒子化。

    模拟中防冲吸能装置及缸体材料与实际完全一致。缸体高强钢的屈服强度取835 MPa,弹性模量206 GPa,泊松比0.3;吸能装置钢材屈服强度取890 MPa,弹性模量206 GPa,泊松比0.26。缸体采用实体单元以C3D8R、单元尺寸30划分网格;吸能装置采用壳单元以S4R、单元尺寸6划分网格。吸能套筒底部施加固定约束,活塞顶部限制水平位移,其余部位不设约束。通用接触,摩擦系数取0.3[2, 6],吸能液压立柱有限元计算模型如图5所示,普通液压立柱模型不设置吸能套筒。

    图  5  立柱有限元模型
    Figure  5.  Finite element model of the column

    对吸能液压支架立柱进行冲击试验。立柱底端安装于试验台下方凹槽内,柱顶施加1000 kN静载,0.9 s后叠加4000 kN冲击荷载。使用传感器记录立柱缸内液体压力、支护阻力等数据,使用2000 Hz高速摄相机捕捉立柱冲击过程中油缸壁标识点的位移量,试验现场环境如图6所示。

    图  6  试验环境照片
    Figure  6.  Experimental environment photos

    立柱阻力在0.9 s叠加冲击后突然增大,并在0.9473 s时增至3662.9 kN,位移由118.36 mm增至282.51 mm,静动转换过程在47.3 ms内完成,液体压力迅速由0.9 s时的26.37 MPa增至91.44 MPa,伴随震动现象,固有频率约为12.5 Hz。在3.24 s时立柱缸体产生屈服现象,出现活塞杆密封失效、液控单向阀损坏等多种失效现象,负载阻力与液体压力均下降,位移持续上升。

    同时进行立柱的模拟加载,模拟采用与试验相同的荷载设置。施加冲击荷载后试验与模拟位移结果对比如图7所示。液压立柱在加荷后约10 ms开始变形,立柱变形迅速增大。在冲击后约38 ms时刻位移曲线逐渐平缓并最终保持稳定,试验立柱位移从3.9 mm增至159.0 mm,增加了155.1 mm;模拟立柱位移从9.2 mm增至171.9 mm,增加了162.7 mm,模拟与试验结果基本一致。

    图  7  现场试验与模拟结果对比
    Figure  7.  Comparison of on-site experimental and simulation results

    在液压活塞上端逐步施加0.6倍工作阻力的集中力,代表液压立柱的初撑力,待系统稳定后施加冲击荷载。防冲支架放置于煤矿巷道内,冲击地压发生后仍需其提供足够的承载力避免二次灾害,故冲击全程预压不卸载。参考国家标准[21],该尺寸液压立柱可使用7 t重锤经自由落体进行冲击试验[25-26]。取系数为1.5,使用4.5、7、10.5 t的重锤以8 m/s的速度对普通液压立柱和吸能液压立柱进行模拟加载。

    以7 t重锤冲击为例,2种液压立柱位移时程如图8所示,图中0时刻为冲击开始时刻。对于普通液压立柱,在冲击后39 ms达到最大位移260 mm;对于吸能液压立柱,吸能装置在冲击后23 ms启动变形、迅速让位,立柱位移从200 mm降低至123 mm,39 ms时吸能装置被完全压溃,装置失效后立柱位移回升,出现第2个位移峰值。不同质量重锤冲击下2种立柱位移最大值见表2。在4.5、7、10.5 t重锤冲击下吸能立柱位移比普通立柱分别降低了16%、23%、30%。

    图  8  立柱位移时程曲线
    Figure  8.  Column displacement time history curve
    表  2  立柱位移最大值
    Table  2.  Maximum column displacement
    立柱类型 不同冲击锤质量下的位移/mm
    4.5 t 7 t 10.5 t
    普通立柱 222 261 291
    吸能立柱 187 199 204
    吸能装置 101 135 136
    下载: 导出CSV 
    | 显示表格

    立柱缸壁应力可反映支护阻力水平。以7 t重锤冲击为例,2种立柱缸壁应力最大时刻应力云图如图9所示。

    图  9  缸壁应力云图
    Figure  9.  Stress cloud map of cylinder wall

    7 t重锤冲击下2种立柱高压缸缸壁应力随时间变化曲线如图10所示。吸能立柱在冲击后23 ms缸壁应力升至600 MPa,吸能装置开始变形,2 ms后缸壁应力达到690 MPa后出现下降,冲击39 ms后吸能装置失效,缸壁应力回升;普通立柱在冲击后26 ms缸壁应力达到材料屈服应力不再进一步上升,在冲击后43 ms开始下降。不同重锤冲击下立柱缸壁应力最大值见表3

    图  10  立柱应力时程曲线
    Figure  10.  Column stress time history curve
    表  3  立柱缸壁应力最大值
    Table  3.  Maximum stress on column cylinder wall
    立柱类型 不同冲击锤质量下的应力/MPa
    4.5 t 7 t 10.5 t
    普通立柱 781 835 835
    吸能立柱 677 690 744
    下载: 导出CSV 
    | 显示表格

    冲击荷载下立柱最大让位速度见表4。吸能液压立柱发生弹性变形时让位速度最快,且具有可恢复性,适应频繁的围岩小震冲击;面对大震时吸能装置塑性变形让位速度稍低,理想情况是粘性小塑性大,当理想塑性时同样具有最快变形速度,实现支架的快速让位,满足防冲支护的速度指标。由表4可见,强冲击下吸能装置压溃速度更快、变形持续时间更短。受立柱回弹与外载双重作用,吸能装置让位速度均大于重锤冲击速度,可以瞬时响应高速围岩冲击。

    表  4  立柱让位速度
    Table  4.  Column yielding velocity
    冲击锤
    质量/t
    液压立柱最大
    让位/(m·s−1)
    吸能装置最大
    让位/(m·s−1)
    吸能装置压溃
    总用时/(m·s−1)
    4.5 t 7.89 11.1 18
    7 t 8.12 13.4 15
    10.5 t 8.33 18.5 13
    下载: 导出CSV 
    | 显示表格

    普通液压立柱仅存在弹性刚度;吸能装置具有可变刚度,可调控吸能液压立柱整体结构刚度。冲击荷载下2种立柱位移−支护阻力曲线如图11所示。吸能液压立柱在吸能装置变形前后表现为弹性刚度,吸能装置压溃期间则表现为塑性刚度,弹性刚度为Kt1= 3.08×107 N/m和Kt2= 3.49×107 N/m,塑性刚度Ks= −4.31×107 N/m,其中塑性刚度为负。弹性让位刚度使支架对围岩具有合理的变形约束能力和弹性缓冲作用,防御小震的频繁冲击;塑性让位刚度能够调控支架过载后的结构刚度,以迅速降低围岩冲击力峰值,满足防冲支护的刚度指标。

    图  11  立柱的位移与支护阻力
    Figure  11.  Column displacement and bracing resistance

    围岩震动频率与支架固有频率接近时,会与支架产生共振,此时较小的震级也会造成支架的严重破坏。理想的吸能装置塑性变形阶段频率为0,能够在立柱过载时瞬时改变立柱固有频率。7 t重锤冲击下对2种立柱冲击荷载下的地基反力进行快速傅里叶变换(FFT),立柱反力和反力频谱如图12所示。吸能液压立柱在受冲变形下频率为0,能够彻底防止立柱及支架在围岩冲击、震动下发生共振。

    图  12  立柱基底反力时程与频谱
    Figure  12.  Base reaction force time history and spectrum of the column

    液压立柱应具有储能吸能能力,保证支架结构不受损坏。液压缸体发生弹性变形储存弹性能,吸能装置则发生塑性变形吸收能量。7 t重锤冲击下2种立柱能量时程如图13所示。普通液压立柱缸体弹性能达到280 kJ后不再上升,并开始产生塑性变形,最终塑性变形能达96 kJ;吸能液压立柱缸体弹性能峰值为204 kJ,比普通立柱降低了27%,吸能装置塑性吸能量为266 kJ。

    图  13  立柱能量时程
    Figure  13.  Column energy time history

    不同重锤冲击下2种立柱让位能量见表5。对于普通液压立柱,在7 t和10.5 t重锤冲击下缸体弹性能达到最大值后产生塑性变形;对于吸能液压立柱,立柱最大弹性能在200 kJ附近,缸体不发生塑性变形。

    表  5  立柱让位能量
    Table  5.  Yielding energy of the column
    指标类型 不同冲击锤质量下的能量/kJ
    4.5 t 7 t 10.5 t
    吸能立柱
    缸体弹性变形能
    198 204 223
    吸能装置
    塑性变形能
    112 266 272
    普通立柱
    缸体弹性变形能
    247 279 280
    普通立柱
    塑性变形能
    0 96 210
    下载: 导出CSV 
    | 显示表格

    在7 t重锤冲击下吸能液压立柱防冲响应如图14,可分为4个阶段。第1阶段为弹性承压阶段,从冲后10 ms时起立柱弹性能与缸壁应力开始稳步攀升。第2阶段为吸能装置作用阶段,在25 ms时刻吸能装置开始被压溃,高速让位吸能,39 ms时刻达到设计位移后停止变形,立柱出现震动现象,在39~44 ms内立柱震荡约3次后吸能装置位移稳定在135 mm;在吸能装置变形后2 ms立柱弹性能达到极大值204 kJ、缸壁应力达到690 MPa,随后开始回落,弹性能经震动保持在70 kJ,缸壁应力300 MPa。第3阶段为阻力回升阶段,在44~61 ms内冲击锤动能逐步被立柱消耗吸收,57 ms时刻立柱弹性能回升至98 kJ并开始保持,缸壁应力499 MPa。第4阶段为泄压阶段,冲击锤动能在61 ms时降为0,立柱弹性能保持约12 ms后逐步降低。

    图  14  立柱防冲响应阶段
    Figure  14.  Column impact response stages

    2种液压支架立柱抗冲指标见表6。围岩冲击能量过大时,吸能液压立柱在阻力回升阶段无法负荷剩余的冲击能,缸体结构将发生塑性变形而失去承载能力,该型号吸能液压立柱可承受最大重锤冲击能为368 kJ,是普通液压立柱的2.3倍。矿用液压支架立柱的泄压阀需要在立柱受冲后尽快启动,保护立柱不被冲坏,工程上一般要求泄压阀在立柱缸壁应力过高而发生塑性应变之前完全启动,对缸体和泄压阀有较高要求。该型号吸能液压立柱泄压阀允许开启时间为39 ms,是普通液压立柱的2倍。

    表  6  立柱的抗冲参数
    Table  6.  Column impact resistance parameters
    立柱类型 可承受最大
    冲击能/kJ
    最大弹性让位
    位移/mm
    泄压阀允许
    开启时间/ms
    普通立柱 160 264 18
    吸能立柱 368 396 39
    下载: 导出CSV 
    | 显示表格

    1)采用有限元和SPH方法建立了吸能液压支架立柱数值模型,使用冲击试验台进行了吸能液压支架立柱受冲试验,试验立柱位移增量与模拟结果差距为5%,模拟与试验结果基本一致。

    2)根据防冲支护设计6项原则探讨了防冲液压支架立柱的防冲性能。让位位移方面,吸能立柱最大位移比普通立柱降低了25%;让位阻力方面,吸能立柱比普通立柱最大让位阻力降低了18%;让位速度方面,吸能装置压溃速度远大于重锤冲击和立柱压溃速度,10.5 t重锤冲击下可在13 ms内完全压溃,让位速度达18.5 m/s;让位刚度方面,吸能立柱具有可变刚度,在吸能装置启动前后表现为弹性刚度,吸能装置压溃时表现为塑形刚度;让位频率方面,吸能立柱让位防冲时的固有频率为0,避免与围岩发生共振;让位能量方面,吸能立柱整体吸能量比普通立柱提高了80%。

    3)总结了吸能液压支架立柱的防冲响应过程可分为4个阶段:弹性承压阶段、吸能装置作用阶段、阻力回升阶段、泄压阶段。发现防冲液压立柱可承受最大冲击能为368 kJ,是普通液压立柱的2.3倍,液压阀开启时间则可达到40 ms,为普通液压立柱的2倍。

  • 图  1   计算域及其边界条件

    Figure  1.   Computational domains and boundary conditions

    图  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

    图  20   CH4泄漏危险范围宽度经验公式

    Figure  20.   Empirical formula of the width 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.

图(20)
计量
  • 文章访问数:  133
  • HTML全文浏览量:  16
  • PDF下载量:  40
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-11-16
  • 网络出版日期:  2023-04-26
  • 刊出日期:  2023-03-14

目录

/

返回文章
返回