• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

Y型通风采煤工作面瓦斯爆炸传播规律模拟研究

刘佳佳 张扬 张翔 聂子硕

张之凡, 李海龙, 张桂勇, 宗智, 姜宜辰. 聚能装药水下爆炸冲击波和侵彻体载荷作用时序研究[J]. 爆炸与冲击, 2023, 43(10): 102201. doi: 10.11883/bzycj-2022-0397
引用本文: 刘佳佳, 张扬, 张翔, 聂子硕. Y型通风采煤工作面瓦斯爆炸传播规律模拟研究[J]. 爆炸与冲击, 2023, 43(8): 085401. doi: 10.11883/bzycj-2023-0018
ZHANG Zhifan, LI Hailong, ZHANG Guiyong, ZONG Zhi, JIANG Yichen. Action time sequence of underwater explosion shock waves and shaped charge projectiles[J]. Explosion And Shock Waves, 2023, 43(10): 102201. doi: 10.11883/bzycj-2022-0397
Citation: LIU Jiajia, ZHANG Yang, ZHANG Xiang, NIE Zishuo. Simulation study on propagation characteristics of gas explosion in Y-shaped ventilated coal face[J]. Explosion And Shock Waves, 2023, 43(8): 085401. doi: 10.11883/bzycj-2023-0018

Y型通风采煤工作面瓦斯爆炸传播规律模拟研究

doi: 10.11883/bzycj-2023-0018
基金项目: 国家自然科学基金(52074106);河南省优秀青年科学基金(232300421061);河南理工大学创新型科研团队项目(T2023-3);深部煤矿采动响应与灾害防控国家重点实验室开放基金(SKLMRDPC22KF11)
详细信息
    作者简介:

    刘佳佳(1985- ),男,博士,副教授,博士生导师,liujiajia@hpu.edu.cn

  • 中图分类号: O382

Simulation study on propagation characteristics of gas explosion in Y-shaped ventilated coal face

  • 摘要: 针对Y型通风采煤工作面的瓦斯爆炸传播规律,利用Fluent模拟软件,结合余吾煤矿N2105工作面实际情况开展了数值模拟研究。结果表明:模拟结果与前人的实验结果之间的最大相对误差为11.3%,最小相对误差仅为1.7%,验证了数学模型的可靠性;确定了瓦斯爆炸数值模拟最合理的关键参数网格尺寸、迭代步长和点火温度分别为0.4 m、0.10 ms和1 800 K;进风顺槽、胶带顺槽、回风巷道和工作面的瓦斯爆炸超压峰值与爆源之间的距离符合指数函数关系,到达超压峰值所需时间与爆源之间的距离符合线性函数关系;距巷道分叉口7.5 m处,工作面超压衰减率为41.03%,胶带顺槽超压衰减率为25.99%,发生爆炸时胶带顺槽内更危险;工作面分叉处湍流区由右侧逐渐向左侧移动,且巷道分叉处超压峰值会增大;回风巷道火焰消散时间最短,胶带顺槽火焰消散时间次之,工作面火焰消散时间最长;胶带顺槽和回风巷道火焰消散方向与瓦斯爆炸初期火焰传播方向相反,工作面火焰消散方向与瓦斯爆炸初期火焰传播方向一致。
  • 对于工业生产的厂房安全来说,气体泄爆是十分重要的问题[1-2]。为了阻止泄爆口产生的火焰和冲击波对周围人员和环境造成危害,通常使用泄爆导管将爆炸产物泄放到安全位置[3]。然而,与简单泄爆容器(无泄爆导管)相比,导管的存在会严重增加容器内产生的爆炸超压[4-5]。已有的研究表明容器内超压剧增主要有以下几个原因:(1)导管内发生二次爆炸;(2)导管内摩擦阻力和气柱惯性;(3)声波阻塞和亥姆霍兹震荡;(4)R-T不稳定性造成的火焰前锋扰动[6-7]。目前,研究者们普遍认为容器内压力快速上升的主要原因是由于导管内发生二次爆炸[8-12]

    为了降低容器爆炸超压,已开展许多针对减缓导管内二次爆炸的研究工作。Henneton等[13]在泄爆导管入口处放置金属网来延迟火焰进入导管的时间,结果表明,火焰在导管入口处熄灭,凭此消除了导管内二次爆炸。然而,如果金属网结构使用不当,会极大增加容器内最大爆炸超压。Molkov等[14]通过向导管内喷水来抑制导管内二次爆炸,得出在此配置下容器最大爆炸超压与无导管泄爆时超压几乎相当。这些研究结果虽然可以控制导管内二次爆炸强度,但并没有从源头上减弱二次爆炸。大量证据表明导管内二次爆炸严重依赖于容器内火焰动力学。Ponizy等[15]证实了容器/导管尺寸与容器内火焰行为密切相关,并且得出在特殊的实验配置下,容器内形成的郁金香型火焰与导管内气柱惯性相结合可以大幅降低导管内二次爆炸强度。Ferrara等[16]研究了导管泄爆过程中内外爆炸相互作用,他们认为容器内火焰传播行为直接影响二次爆炸程度,之后,外部爆炸反过来影响容器内最终爆炸超压。意识到容器内火焰动力学与导管内二次爆炸以及容器内最大超压的关系,我们考虑通过改变容器内火焰动力学从源头上降低二次爆炸严重程度以及容器内最大爆炸超压。

    为了保护容器(包围体),本文选择在容器内直接添加粉体抑制剂进行抑爆。由于粉体抑制剂具有高抑制效率,易于长期储存等特点,已经被广泛应用于实际工业生产中[17-18]。事实上,碳酸氢钠(NaHCO3)干粉作为一种典型粉体抑制剂可以有效降低可燃气体和粉尘爆炸严重程度[19]。实验和理论研究表明,由于粉体颗粒和火焰间热传递,在可燃气体爆炸火焰区域添加NaHCO3粉体可以大幅降低爆炸强度。例如:Jiang等[20]通过实验比较了NaHCO3和NH4H2PO4对生物质粉尘爆炸的抑制效果,得出抑制剂会消耗生物质火焰中的关键自由基,导管火焰速度降低,火焰温度降低,最终得出NaHCO3对生物质粉尘爆炸的抑制性能优于NH4H2PO4

    本文中在连接不同长度泄爆导管(Kv=4.87,L=250, 500, 750 mm)的5 L垂直容器内,研究当量比为1的甲烷/空气预混气爆炸特性。考虑的泄爆导管长度符合国内外泄爆导管的设计标准。实验所用NaHCO3粉体质量浓度C依次为0、40、80、120、160、200和240 g/m3。主要研究容器和泄爆导管内火焰传播特性(火焰结构,火焰传播速度)和压力特性(压力历史和最大超压),以期望探讨出粉体对导管泄爆过程的影响机制,为预防可燃气体爆炸和减轻爆炸危害提供理论依据。

    实验系统如图1所示,由爆炸管道系统、点火系统、配气系统、喷粉系统、数据采集系统、同步控制系统六部分组成。爆炸管道系统包括100 mm×100 mm×500 mm的有机透明玻璃容器和顶部通过法兰连接不同长度的泄爆导管。泄爆导管长度分别为250、500和750 mm,导管泄爆系数为4.87,(泄爆系数由公式Kv=V2/3 /Av计算,其中V是容器体积,Av是泄爆导管截面积)。容器下端用不锈钢法兰盘密闭,上端用PVC薄膜密封,PVC膜置于容器与泄爆导管之间。2个压力传感器分别安装在容器底端和距离泄爆导管底端20 mm位置,压力和光信号的采集频率为15 kHz。使用德国Lavision 4G高速相机以2 000 s−1的频率拍摄火焰瞬变图像,捕捉火焰前锋位置和火焰结构。配气系统采用两个质量流量计控制通入体积分数为9.5%的甲烷/空气预混气,持续10 min,以保证3~5倍容器体积预混气流经容器,确保容器内混合物均匀[21]。配气完成后,同时关闭进气阀和排气阀,静置30 s。在此期间将制备的体积分数为9.5%的甲烷/空气预混气通入高压储气瓶中至0.5 MPa,关闭气瓶通气阀。启动同步控制器,电磁阀开启,喷粉400 ms后,点火器放电,引燃容器内预混气体,数据采集系统同时动作。根据大量预实验(冷态实验)获得的喷粉压力和点火延迟时间,可以确保粉体在容器内尽可能均匀分布。实验过程中,为了确保数据准确性,每组工况至少进行3次以上重复性实验。

    图  1  实验系统示意图
    Figure  1.  Schematic diagram of experimental system

    众多研究发现:NaHCO3粉体作为一种化学型抑爆剂可以有效降低火灾和爆炸强度。因此,本文选择使用NaHCO3粉体作为抑爆剂。采用经过标准筛制备的200~220目NaHCO3粉体,利用Mastersizer 3 000超高速智能粒度分析仪测试粒度分布结果如图2所示,粉体平均粒径为106 μm。利用公式δ =D90D10/D50计算NaHCO3粉体粒度分散性为0.895(D10D50D90指的是体积百分比分别为10%,50%,90%的直径;D[3,2]为表面平均粒径;D[4,3]为体积平均粒径)。图3为粉体颗粒放大50倍的扫描电镜图。结合粒度分布图和扫描电镜图来看,可见NaHCO3粉体颗粒大小均匀性较好。

    图  2  NaHCO3粉体的粒度分布
    Figure  2.  Particle size distribution of NaHCO3 powder
    图  3  NaHCO3粉体的扫描电镜图
    Figure  3.  Scanning electron microscope of the NaHCO3 powder

    图4L=750 mm的导管为例分析了不同质量浓度NaHCO3粉体对甲烷/空气爆炸火焰锋面结构的影响。图中红色虚线表示容器和泄爆导管间法兰连接处。每个工况挑选出代表性锋面结构图,每组图挑选的第1幅为点火后初期火焰锋面情况,第2幅表示火焰传播至容器中部,第3幅为火焰到达容器末端,接下来是火焰在导管中传播的锋面结构演化图。火焰在导管-容器初期传播过程与简单泄爆容器(不加泄爆导管)几乎相同,都经历了球型和指型阶段,之后火焰发展则完全不同[22]。在导管泄爆配置中容器末端的收缩区域,火焰在加速流中被严重拉伸且收缩,之后,火焰进入泄爆导管,受容器-导管截面突然变化引起的强烈扰动作用,在湍流混合区(约等于导管截面宽度距离)火焰前锋变形及褶皱,此时,火焰极可能由于热量不断向管壁散失发生淬火现象,当导管内热气体含量持续增加能够被再次点燃时,发生二次爆炸[15]。然而,淬火现象并不是发生二次爆炸的必须条件。当火焰进入导管能够直接点燃之前泄入导管的未燃气时,则不会发生淬火,火焰连续传播至导管出口。

    图  4  火焰锋面结构演化图
    Figure  4.  Evolution diagram of flame front structure

    图4所示,以3种典型NaHCO3粉体质量浓度(C=0,160,200 g/m3)为例分析了质量浓度对甲烷爆炸锋面结构演化影响。由图可知,随着NaHCO3粉体质量浓度增加,容器内火焰前锋逐渐模糊不均匀且撕裂情况加剧。对于低质量浓度NaHCO3粉体(C=0,40,80,120 g/m3),容器内火焰破碎程度较小,未燃气与已燃气体边界清晰,火焰传播至容器末端所需时间较短,依次为24、25、26和29 ms。然而,由于较低质量浓度NaHCO3粉体(0~120 g/m3)导致泄入导管内粉体质量浓度较小,火焰尽管在导管初始部分受到强烈扰动,但并未发生淬熄,发生二次爆炸后,火焰可以连续不断传至导管末端。当粉体质量浓度达到160 g/m3时,容器内火焰前锋出现局部熄灭,火焰形态破碎程度加大,在容器末端未燃区内形成一个小孤立燃烧区,此时火焰传播至容器末端时间为36 ms。火焰进入导管后,由于较高质量浓度NaHCO3粉体导致泄入导管内粉体质量浓度较大,在导管初始部分发生淬火。大量对可燃气爆炸的研究表明,火焰锋面结构对火焰的发展起到至关重要的作用[23]

    淬火现象发生标志火焰速度降低,反应进行缓慢,这将会导致一个较弱的二次爆炸。而后,被扰动的火焰在导管下游能够引燃持续聚集的混合物时,火焰在导管中继续传播至导管出口。随着NaHCO3粉体质量浓度进一步增加(C=200,240 g/m3),容器内火焰破碎程度进一步加大,甚至在容器末端熄灭,火焰不能传至泄爆导管,导管内也不会发生二次爆炸,火焰传播至距离容器底端最远处所需时间分别为48和51 ms。由此判断,高效率抑制的火焰前锋不连续且破碎程度严重。较高质量浓度NaHCO3粉体导致容器内火焰破碎程度更大的原因如下:一方面,NaHCO3粉体吸热分解释放CO2和H2O,导致火焰锋面热量散失以及混合物稀释。另一方面,由于典型自由基清除剂(NaOH)质量浓度增加,加速了爆炸链式反应中自由基与气相分解产物之间重组反应。同时,由于抑制剂的混入并不是绝对均匀,因此,火焰锋面沿轴向分布不规则。

    对于不同长度泄爆导管,火焰传播速度随火焰前锋位置变化趋势几乎一致,因此,图5以500 mm导管为例分析了不同质量浓度NaHCO3粉体作用下各锋面位置点处所对应的火焰传播速度变化趋势。图中蓝绿色虚线表示火焰传播至容器尾部(距离容器底端500 mm处)。如图5所示,不同质量浓度NaHCO3粉体所呈现火焰传播速度趋势不同,火焰传播速度随粉体质量浓度增加主要呈现以下3种典型趋势:无粉体作用和粉体质量浓度较低时(0和40 g/m3),火焰传播速度在容器和导管内演化均呈现单调增加趋势。随着粉体质量浓度增加(80和120 g/m3),火焰传播速度随火焰前锋位置变化趋势发生改变,火焰前锋速度先加速后减速再加速,且后期加速(导管内)明显快于前期(容器内)。从图中观察到火焰前锋速度减速阶段发生在导管入口处,与无粉体或较低粉体质量浓度时相比,较高质量浓度NaHCO3粉体导致容器内火焰速度降低,火焰进入导管时,较低火焰传播速度更易受到容器/导管截面突然变化产生的扰动作用,因此,火焰速度呈现下降趋势,之后沿着导管火焰传播速度继续增加。此现象对于160 g/m3的粉体更加明显,火焰在导管内入口处淬火(图4),这极大地降低了导管内的火焰传播速度。对于200和240 g/m3的粉体,火焰只在容器内传播,火焰传播速度几乎恒定。此外,随着NaHCO3粉体质量浓度的增加,容器与导管内的火焰传播速度均逐渐下降,且与导管长度无关。火焰前锋传播至容器末端(距离容器底端500 mm处),火焰速度降低幅度最大(与无粉体时相比),随着粉体质量浓度的增加,容器末端火焰传播速度依次为40、34、24、14、12和5 m/s,最大速度下降率分别达20%、32%、52%、72%、76%和90%。

    图  5  火焰传播速度随位置变化关系
    Figure  5.  Flame propagation velocity changed with flame front position

    图6所示为泄爆导管内最大火焰传播速度随NaHCO3粉体质量浓度变化趋势。由图6可知,随着NaHCO3粉体质量浓度增加,导管内最大火焰传播速度逐渐下降,这与容器末端火焰传播速度变化趋势一致,表明NaHCO3粉体通过控制容器内火焰传播速度进而降低导管内最大火焰传播速度,这意味着火焰从容器末端进入导管的速度越快,导管内最大火焰传播速度越大。此外,如图6所示,导管长度对最大火焰传播速度也有较大影响,长导管内最大火焰传播速度更大,且这一规律与NaHCO3粉体质量浓度无关,长导管内最大火焰传播速度更大的原因同样是由于容器末端较高的火焰传播速度。由于火焰传播速度主要来源于火焰前沿未燃气速度,与流场湍流强度正相关,因而较高火焰传播速度将产生较强的流场湍流度[24]。Ponizy[15]通过实验研究与理论模拟,认为二次爆炸的强弱与导管内入口处的湍流强度成密切相关,湍流程度越剧烈,二次爆炸强度越大。由此可知,NaHCO3粉体质量浓度越高,容器和导管内火焰传播速度越小,在导管入口处形成一个弱湍流区,产生较弱二次爆炸。因此,NaHCO3粉体可以通过降低容器内火焰传播速度有效抑制导管内二次爆炸。同时,由于长导管较大的最大火焰传播速度,将在导管内产生一个更剧烈的二次爆炸。

    图  6  导管内最大火焰传播速度随NaHCO3质量浓度变化趋势
    Figure  6.  Maximum flame propagation velocity in the duct changed with NaHCO3 powder mass concentration

    图7为不同质量浓度的NaHCO3粉体作用下火焰前锋到达容器末端时间t以及容器内平均火焰速度uf的关系曲线(L=750 mm)。分析这2个参量有利于更好地理解NaHCO3粉体质量浓度与二次爆炸的关系。由图7可知,随着NaHCO3粉体质量浓度增加,t逐渐延长,从24 ms(0 g/m3)延长到51 ms(240 g/m3),理论研究证明,粉体抑爆效果与颗粒在反应区的停留时间有关[25],容器内较高质量浓度粉体延长火焰达到容器末端时间,使得NaHCO3粉体有更充分的时间分解发挥抑制作用,致使容器内抑爆作用增强。同时,火焰进入泄爆导管时间越晚,导管内二次爆炸越弱。这是由于二次爆炸强度与导管入口未燃气相对浓度密切相关,火焰进入泄爆导管越晚,之前由于容器内爆炸泄入导管的未燃气有更多时间从导管末端泄放,最终导致一个较弱的二次爆炸。如图所示,uf随着NaHCO3粉体质量浓度增加不断下降,对于0、40、80、120、160、200、240 g/m3的NaHCO3粉体,uf依次为20.8、20.0、19.2、17.2、13.9、10.4和9.8 m/s。uf与导管内火焰速度正相关,而导管内火焰速度与导管内最大压力pd,max正相关。此外,对于同一工况而言,较高质量浓度的NaHCO3粉体会导致容器内粉体泄入导管内的质量浓度升高,这会加重导管内火焰的淬熄,使得导管内的火焰传播速度更低。根据这一系列分析结果,均证实较高NaHCO3粉体质量浓度会产生一个较弱二次爆炸。

    图  7  火焰前锋到达容器末端时间与容器内平均火焰速度
    Figure  7.  Arrival time of flame front and average flame velocity in the vessel

    不同质量浓度NaHCO3粉体对压力发展的影响明显不同。图8L=750 mm为例比较了40和200 g/m3 NaHCO3粉体的超压波形。如图所示,无论NaHCO3粉体质量浓度如何,压力波形曲线都有2个压力峰值p1p2。分析这2个压力峰值产生的原因如下:p1是由于火焰在容器内传播时触壁引起的,火焰表面积和燃烧速率从这一时刻起下降,导致容器压力降低。对于无粉体或粉体质量浓度较低时(C在0~160 g/m3范围内),p2是由于泄爆导管内发生二次爆炸,容器内爆炸产生的压力波携带未燃气在火焰之前进入泄爆导管,之后,进入导管的火焰点燃未燃气,发生二次爆炸。导管内二次爆炸引起管内压力突升,致使容器泄爆效率下降,使得容器内产生p2;而对于高质量浓度粉体(C=200,240 g/m3),p2主要来源于容器内的残余燃烧。由图8可知,当NaHCO3粉体质量浓度为40 g/m3时,p2(43.4 kPa)大于p1(9.0 kPa),然而,NaHCO3粉体质量浓度为200 g/m3时与40 g/m3情况相反,p2(1.5 kPa)小于p1(4.2 kPa),p2下降96.5%。由于p2主要是导管内发生二次爆炸所致,因此,可以判断高质量浓度NaHCO3粉体对二次爆炸有显著影响。Kasmani等[26]认为泄爆口的高气体流速和随后的火焰速度是产生剧烈二次爆炸的主要原因。不同质量浓度NaHCO3粉体严重影响容器内火焰动力学(火焰结构,火焰传播速度)。高质量浓度NaHCO3粉体极大降低容器内平均火焰速度,延长火焰进入导管时间(见图7),最终产生一个较弱的二次爆炸。与低质量浓度NaHCO3时相反,高质量浓度下容器内的最大压力依赖p1。因此,合适质量浓度NaHCO3粉体可以改变容器压力上升机理。

    图  8  40和200 g/m3 NaHCO3作用下容器内超压变化历程
    Figure  8.  Pressure history in the vessel with 40 and 200 g/m3 NaHCO3

    图9给出了连接不同长度导管的容器内两特征压力峰值(p1p2)随NaHCO3粉体质量浓度变化趋势,如图所示,对于不同长度泄爆导管,p1p2均随着NaHCO3粉体质量浓度增加不断下降,而且,NaHCO3p2的抑制效率高于p1,尤其对于较长导管(500 mm和750 mm)。随着NaHCO3粉体质量浓度增加,p2p1的差值Δp不断下降,甚至当粉体质量浓度达到160 g/m3时,Δp变为负值,且这一规律与导管长度无关。分析Δp负值的原因是,较高NaHCO3质量浓度(160 g/m3)极大削弱导管内二次爆炸,由二次爆炸产生的降低容器内泄爆速率这一作用也被减弱,甚至更高的NaHCO3质量浓度(200 g/m3和240 g/m3)致使火焰在容器内熄灭(图4),消除导管内二次爆炸,增加容器最大超压的相关效应(或机理)也被消除。此时p2主要来源于容器内微弱的残余燃烧。因此,随着NaHCO3粉体质量浓度增加,二次爆炸及其相关作用不断减弱,且当NaHCO3粉体质量浓度达到某一临界值时(200 g/m3),二次爆炸被完全抑制。在这种情况下,容器内最大爆炸超压pv,max由火焰在容器内触壁时产生的第一压力峰值p1决定。

    图  9  容器内两压力峰值随粉体质量浓度的变化
    Figure  9.  Variation of two pressure peaks in the vessel with NaHCO3 mass concentration

    爆炸事故中,爆炸超压是评定爆炸危险程度的重要参数。图10是容器内最大超压峰值(pv,max)及其下降率随NaHCO3粉体质量浓度的变化关系。随着粉体质量浓度增加(0~240 g/m3),不同长度导管的pv,max不断下降,且粉体质量浓度达到一个临界值(160 g/m3)时,pv,max逐渐趋于稳定,这表明NaHCO3粉体的质量浓度效应逐渐减弱。此外,长导管(750 mm)会产生更大pv,max,这是由于连接长导管的容器和导管内火焰传播速度更大(图6),导致一个更剧烈的二次爆炸,由二次爆炸造成的容器泄爆效率的下降程度远大于短导管,致使容器内净体积增加率增加,因此产生更大pv,max。本文定义超压峰值下降率为有无粉体时容器最大超压差值与无粉体时最大超压的比值为θθ可以有效体现导管泄爆过程中不同质量浓度粉体的抑制效率。如图10所示,随着导管长度增加,NaHCO3粉体抑爆效率不断增加。与无导管泄爆时相比,导管泄爆时θ明显得到提升。例如:郑立刚等[27]研究数据表明,在无泄爆导管的容器中加入200 g/m3的NaHCO3θ(与无粉体时相比)只能达到48%,而在本文的结果中,当NaHCO3质量浓度为200 g/m3时,对于不同长度导管(L=250,500,750 mm),θ依次为68%、90%和93%。一方面,泄爆导管的存在降低了容器泄爆速率,使得容器内参与抑制的粉体质量浓度更高,提升粉体抑制效率。另一方面,粉体对容器内火焰动力学的影响削弱导管内二次爆炸强度,进而大大降低容器内超压。此外,随着粉体质量浓度增加,θ不断增加,当粉体质量浓度达到临界值(160 g/m3)时,θ趋于稳定,这一规律与导管长度无关。同时,NaHCO3粉体由40~160 g/m3超压峰值下降率的增量要大于160~240 g/m3。这是因为NaHCO3粉体质量浓度较低时,低热阻导致反应区热量更易传至颗粒内部,在高质量浓度粉体下,热阻增加,限制反应区热量向颗粒内部扩散,导致只有颗粒外围部分分解,而颗粒内部由于热量传导受阻不能有效分解,且随着粉体质量浓度增加,热量向颗粒内部传导时间增加,这不利于粉体分解发挥抑制作用[28]。因此,对于相同导管配置,低质量浓度粉体分解效果比高质量浓度的分解效果好,随着NaHCO3粉体质量浓度增加,质量浓度效应逐渐减弱。

    图  10  容器内最大超压峰值及下降率曲线
    Figure  10.  Dependence of maximum overpressure and its drop rate in the vessel on NaHCO3 powder mass concentration

    图11列举了连接不同长度导管的容器内平均火焰传播速度uf与导管内最大爆炸超压pd,max的关系。由图11可知,对于一个给定长度的导管,uf越大,pd,max越大,也就是说容器内火焰传播越快,导管内二次爆炸越剧烈。从这一结果可以得出,火焰在容器内传播过程直接影响导管内二次爆炸大小。对于相同的uf,长导管(750 mm)内产生更大pd,max的原因是,连接长导管的容器内火焰传播速度更大[29],容器内较快的火焰速度导致导管内火焰传播速度更大(图6),进而产生一个较强的湍流场。二次爆炸的剧烈程度与导管内湍流强度正相关。因此,容器内火焰传播速度越快,pd,max越大。

    图  11  导管内最大超压与容器内平均火焰传播速度的关系
    Figure  11.  Maximum pressure in the duct pd,max vs average flame velocity in the vessel uf

    图12所示为容器最大压力pv,max与导管最大压力pd,max的关系。由图12可知,pv,maxpd,max存在线性关系,随着pd,max增加,pv,max不断增加,且这一线性关系不依赖于NaHCO3粉体质量浓度与导管长度。实际上,导管内发生二次爆炸导致容器内泄爆效率大幅降低,使得容器产生最大爆炸超压pv,max,且二次爆炸越剧烈,泄爆效率下降的幅度越大,容器内净体积增加率增加越快,容器最大超压越大,也就是说导管内二次爆炸越剧烈,容器内产生的最终爆炸超压越大。不同质量浓度粉体虽然通过降低uf来削弱pd,max,这只能导致不同的二次爆炸强度,但pd,maxpv,max的线性关系并不会因此发生改变。根据pd,maxpv,max的线性关系以及ufpd,max的特殊关系,可以判断出ufpv,max并不是严格的正相关(即非单值函数),而是对于某一固定长度导管,pv,maxuf的增加而增加。本文实验在小尺度条件下进行,由于大尺度实验和小尺度实验还存在一些差异,在后续的研究中,将在大尺度实验配置中进一步研究干粉对导管泄爆过程的影响,为现实应用提供更加精细准确的理论参考。

    图  12  容器内最大爆炸超压与导管内最大爆炸超压的关系
    Figure  12.  Maximum pressure in the vessel pv,max vs maximum pressure in the duct pd,max

    本文中研究了NaHCO3粉体对导管泄爆过程的影响,从火焰锋面结构、火焰传播速度,压力波形和超压峰值等多个方面进行分析比较,得到以下结论。

    (1)抑爆剂质量浓度严重影响容器和导管内火焰传播过程。随着粉体质量浓度增加,容器内火焰破碎程度加大,导管内火焰逐渐弱化甚至熄灭;不同质量浓度粉体展现3种火焰速度模式。这是由于不同质量浓度粉体导致容器与导管内粉体分布情况有所差异,进而影响抑制效率。综合分析可得:高效率抑制的火焰前锋越来越弱,容器与导管内火焰传播速度均逐渐降低。

    (2)抑爆剂质量浓度能改变容器最大爆炸超压的控制机理,最终影响容器最大超压。粉体质量浓度较低时,容器最大爆炸超压取决于二次爆炸产生的第二压力峰值,反之取决于火焰在容器触壁时产生的第一压力峰值。

    (3)存在一个临界质量浓度(160 g/m3)使得粉体抑制效率逐渐趋于稳定。随着粉体质量浓度增加,容器内最大爆炸超压和超压峰值下降率的变化量逐渐减少,表明质量浓度效应越来越弱,粉体抑制效率受到限制。

    (4)容器最大爆炸超压与容器内平均火焰传播速度和导管内最大超压密切相关。容器最大爆炸超压与导管内最大超压呈线性关系,且与粉体质量浓度和导管长度均无关。对于给定导管尺度,容器最大爆炸超压和导管内最大超压均随着容器内平均火焰传播速度的增加而增加,即容器内火焰传播速度越快,最终爆炸超压越高,但这种正相关关系不是单值的,与导管长度有关。

  • 图  1  余吾煤矿N2105工作面的三维物理模型

    Figure  1.  A 3D physical model of N2105 working face in Yuwu Coal Mine

    图  2  模拟管道示意图

    Figure  2.  Schematic diagram of simulated pipeline

    图  3  参数优化监测点示意图

    Figure  3.  Parameter optimization monitoring points

    图  4  超压峰值随网格尺寸的变化

    Figure  4.  Variations of overpressure peak with mesh size

    图  5  超压峰值到达时间随网格尺寸的变化

    Figure  5.  Variations of overpressure peak arrival time with mesh size

    图  6  超压峰值随迭代步长的变化

    Figure  6.  Variation of overpressure peak with time step

    图  7  超压峰值到达时间随迭代步长的变化

    Figure  7.  Variation of overpressure peak arrival time with time step

    图  8  超压峰值随点火温度的变化

    Figure  8.  Variation of overpressure peak with ignition temperature

    图  9  超压峰值到达时间随点火温度的变化

    Figure  9.  Variation of overpressure peak arrival time with ignition temperature

    图  10  瓦斯爆炸超压传播云图

    Figure  10.  Contours of gas explosion overpressure propagation

    图  11  回风巷道瓦斯爆炸传播规律

    Figure  11.  Gas explosion propagation law in return airway

    图  12  胶带顺槽瓦斯爆炸传播规律

    Figure  12.  Gas explosion propagation law along belt fluting

    图  13  工作面瓦斯爆炸传播规律

    Figure  13.  Gas explosion propagation law of working face

    图  14  进风顺槽瓦斯爆炸传播规律

    Figure  14.  Gas explosion propagation law in the inlet channel

    图  15  爆炸超压峰值对比

    Figure  15.  Comparison of gas explosion overpressure peaks

    图  16  爆炸超压峰值到达时间对比

    Figure  16.  Comparison of overpressure peak arrival time of gas explosion

    图  17  瓦斯爆炸超压矢量图

    Figure  17.  Vector diagrams of gas explosion overpressure

    图  18  火焰的温度云图

    Figure  18.  Contours of flame temperature variation

    表  1  爆炸超压模拟结果与实验结果的对比

    Table  1.   Comparison between simulation and experimental results of explosion overpressure

    爆心距/m爆炸超压/kPa相对误差/%
    实验模拟
    2.0181.850193.915 6.6
    2.5178.860175.819 –1.7
    3.0167.663156.085 –6.9
    3.5159.170141.142–11.3
    4.0139.672129.495 –7.3
    4.5120.682117.704 –2.4
    下载: 导出CSV

    表  2  网格分布

    Table  2.   Grid distribution

    网格尺寸/m网格节点数网格数平均网格质量
    0.824 40217 5350.985
    0.648 18836 8340.996
    0.577 52861 6001.000
    0.4154 280128 8170.999
    0.21 087 294990 0000.999
    下载: 导出CSV
  • [1] 中国能源中长期发展战略研究项目组. 中国能源中长期2030—2050发展战略研究 [M]. 北京: 科学出版社, 2011.
    [2] EMAMI S D, RAJABI M, HASSAN C R C, et al. Experimental study on premixed hydrogen/air and hydrogen-methane/air mixtures explosion in 90 degree bend pipeline [J]. International Journal of Hydrogen Energy, 2013, 38(32): 14115–14120. DOI: 10.1016/j.ijhydene.2013.08.056.
    [3] BLANCHARD R, ARNDT D, GRATZ R, et al. Explosions in closed pipes containing baffles and 90 degree bends [J]. Journal of Loss Prevention in the Process Industries, 2010, 23(2): 253–259. DOI: 10.1016/j.jlp.2009.09.004.
    [4] 孟显华, 谢岩森. 甲烷体积分数及巷道结构对甲烷爆炸特性影响研究 [J]. 中国安全科学学报, 2021, 31(S1): 136–142. DOI: 10.16265/j.cnki.issn1003-3033.2021.S1.024.

    MENG X H, XIE Y S. Study on influence of methane’s volume fraction and roadway structure on methane explosion characteristics [J]. China Safety Science Journal, 2021, 31(S1): 136–142. DOI: 10.16265/j.cnki.issn1003-3033.2021.S1.024.
    [5] YUAN X Q, ZHOU J, LIN Z Y, et al. Numerical study of detonation diffraction through 90-degree curved channels to expansion area [J]. International Journal of Hydrogen Energy, 2017, 42(10): 7045–7059. DOI: 10.1016/j.ijhydene.2017.01.206.
    [6] CAO X Y, FAN L T, CUI Y Y, et al. Effect of connected vessels structure on methane explosion characteristics [J]. Journal of Loss Prevention in the Process Industries, 2022, 77: 104789. DOI: 10.1016/j.jlp.2022.104789.
    [7] ZHANG L L, YANG Q Y, SHI B M, et al. Influences of a pipeline’s bending angle on the propagation law of coal dust explosion induced by gas explosion [J]. Combustion Science and Technology, 2021, 193(5): 798–811. DOI: 10.1080/00102202.2019.1673381.
    [8] QIU J W, JIANG B Y, TANG M Y, et al. Effect of different bend pipes on the propagation characteristics of premixed methane-air explosion in confined spaces [J]. Geofluids, 2021, 2021: 6635156. DOI: 10.1155/2021/6635156.
    [9] LIU J J, ZHANG Y, CHEN S Q, et al. Simulation study of gas explosion propagation law in coal mining face with different ventilation modes [J]. Frontiers in Energy Research, 2022, 10: 846500. DOI: 10.3389/fenrg.2022.846500.
    [10] 马恒, 陈晓军, 荆德吉. H型通风巷道瓦斯爆炸及泄爆过程模拟研究 [J]. 中国安全科学学报, 2021, 31(1): 45–51. DOI: 10.16265/j.cnki.issn1003-3033.2021.01.007.

    MA H, CHEN X J, JING D J. Simulation study on gas explosion and discharge process in H-type ventilation roadway [J]. China Safety Science Journal, 2021, 31(1): 45–51. DOI: 10.16265/j.cnki.issn1003-3033.2021.01.007.
    [11] 高建良, 吴泽琳, 王文祺, 等. 瓦斯爆炸冲击波在角、并联巷道内传播规律对比研究 [J]. 安全与环境学报, 2021, 21(6): 2494–2499. DOI: 10.13637/j.issn.1009-6094.2020.1012.

    GAO J L, WU Z L, WANG W Q, et al. Comparative study on the propagation law of gas explosion shock wave in the diagonal and parallel roadway [J]. Journal of Safety and Environment, 2021, 21(6): 2494–2499. DOI: 10.13637/j.issn.1009-6094.2020.1012.
    [12] 高智慧, 李雨成, 张欢, 等. 瓦斯爆炸在角联通风管网中的传播特性研究 [J]. 中国安全生产科学技术, 2022, 18(8): 72–78. DOI: 10.11731/j.issn.1673-173x.2022.07.001.

    GAO Z H, LI Y C, ZHANG H, et al. Study on propagation characteristics of gas explosion in diagonal ventilation pipe network [J]. Journal of Safety Science and Technology, 2022, 18(8): 72–78. DOI: 10.11731/j.issn.1673-173x.2022.07.001.
    [13] 景国勋, 朱斯佳, 贺祥. 不同封闭情况下T型中瓦斯爆炸传播规律实验研究 [J]. 中国安全生产科学技术, 2021, 17(7): 84–90. DOI: 10.11731/j.issn.1673-193x.2021.07.014.

    JING G X, ZHU S J, HE X. Experimental study on propagation laws of gas explosion in T-shaped pipe under different enclosed conditions [J]. Journal of Safety Science and Technology, 2021, 17(7): 84–90. DOI: 10.11731/j.issn.1673-193x.2021.07.014.
    [14] 丁浩, 杜玉晶, 解北京. 负压腔体抑制T型瓦斯爆炸数值模拟 [J]. 煤矿安全, 2021, 52(7): 1–8. DOI: 10.13347/j.cnki.mkaq.2021.07.001.

    DING H, DU Y J, XIE B J. Numerical simulation of negative pressure cavity suppression of T-pipe gas explosion [J]. Safety in Coal Mines, 2021, 52(7): 1–8. DOI: 10.13347/j.cnki.mkaq.2021.07.001.
    [15] YAN Z, YUAN S, LI Z, et al. Study on inhibitory effect of cavity on gas explosion propagation [J]. Geofluids, 2021, 2021: 6674425. DOI: 10.1155/2021/6674425.
    [16] 周辉, 穆朝民, 刘伟, 等. 空腔体长度对抑制瓦斯爆炸性能的影响研究 [J]. 中国安全科学学报, 2019, 29(9): 107–112. DOI: 10.16265/j.cnki.issn1003-3033.2019.09.017.

    ZHOU H, MU C M, LIU W, et al. Influence of cavity length on gas explosion suppression [J]. China Safety Science Journal, 2019, 29(9): 107–112. DOI: 10.16265/j.cnki.issn1003-3033.2019.09.017.
    [17] 时本军, 穆朝民, 马海峰, 等. 腔体影响全巷道甲烷爆炸冲击波传播的特性 [J]. 煤炭学报, 2020, 45(S2): 841–849. DOI: 10.13225/j.cnki.jccs.2020.0791.

    SHI B J, MU C M, MA H F, et al. Cavity effect on the characteristics of methane blast wave propagation in the whole roadway [J]. Journal of China Coal Society, 2020, 45(S2): 841–849. DOI: 10.13225/j.cnki.jccs.2020.0791.
    [18] WANG J, WU Y, ZHENG L G, et al, Study on the propagation characteristics of hydrogen/methane/air premixed flames in variable cross-section ducts [J]. Process Safety and Environmental Protection, 2020, 135: 135–143. DOI: 10.1016/j.psep.2019.12.029.
    [19] LI Z, CHEN L, YAN H C, et al. Gas explosions of methane-air mixtures in a large-scale tube [J]. Fuel, 2021, 285: 119239. DOI: 10.1016/j.fuel.2020.119239.
    [20] GAO K, LIU Z Y, WU C Q, et al. Effect of low gas concentration in underground return tunnels on characteristics of gas explosions [J]. Process Safety and Environmental Protection, 2021, 152: 679–691. DOI: 10.1016/j.psep.2021.06.045.
    [21] 罗振敏, 刘利涛, 王涛, 等. C2H6、C2H4、CO与H2对甲烷爆炸压力及动力学特性影响 [J]. 工程科学学报, 2022, 44(3): 339–347. DOI: 10.13374/j.issn2095-9389.2020.10.22.002.

    LUO Z M, LIU L T, WANG T, et al. Effect of C2H6, C2H4, CO and H2 on the explosion pressure and kinetic characteristics of methane [J]. Chinese Journal of Engineering, 2022, 44(3): 339–347. DOI: 10.13374/j.issn2095-9389.2020.10.22.002.
    [22] 梁国栋. 2种气体状态下瓦斯爆炸极限范围的研究 [J]. 现代矿业, 2019, 35(9): 211–213. DOI: 10.3969/j.issn1674-6082.2019.09.064.

    LIANG G D. Study on the limit range of gas explosion in two kinds of gas conditions [J]. Modern Mining, 2019, 35(9): 211–213. DOI: 10.3969/j.issn1674-6082.2019.09.064.
    [23] 贾泉升, 司荣军, 李润之, 等. 定容条件下瓦斯爆炸超压及爆后气体成分试验研究 [J]. 煤矿安全, 2019, 50(12): 1–5. DOI: 10.13347/j.cnki.mkaq.2019.12.001.

    JIA Q S, SI R J, LI R Z, et al. Experimental study on overpressure and gas composition after gases explosion under constant volume condition [J]. Safety in Coal Mines, 2019, 50(12): 1–5. DOI: 10.13347/j.cnki.mkaq.2019.12.001.
    [24] 任搴华, 赵耀江, 周家川. Y型通风系统采空区漏风规律相似模拟实验研究 [J]. 中国煤炭, 2012, 38(3): 98–101. DOI: 10.19880/j.cnki.ccm.2012.03.026.

    REN Q H, ZHAO Y J, ZHOU J C. Similar simulation of air leakage law in goaf with Y-shape ventilation system [J]. China Coal, 2012, 38(3): 98–101. DOI: 10.19880/j.cnki.ccm.2012.03.026.
    [25] 徐景德. 矿井瓦斯爆炸冲击波传播规律及影响因素的研究 [D]. 北京: 中国矿业大学, 2003: 31–33.

    XU J D. The study on gas explosion wave propagation law and affected factors in mine [D]. Beijing, China: China University of Mining and Technology, 2013: 31–33.
    [26] 陈鹏, 李艳超, 黄福军, 等. 方孔障碍物对瓦斯火焰传播影响的实验与大涡模拟 [J]. 爆炸与冲击, 2017, 37(1): 21–26. DOI: 10.11883/1001-1455(2017)01-0021-06.

    CHEN P, LI Y C, HUANG F J, et al. LES approach to premixed methane/air flame propagating in the closed duct with a square-hole obstacle [J]. Explosion and Shock Waves, 2017, 37(1): 21–26. DOI: 10.11883/1001-1455(2017)01-0021-06.
    [27] 李蒙, 杜扬, 李国庆, 等. 含90°直角弯管结构受限空间油气泄压爆炸实验与大涡模拟研究 [J]. 化工学报, 2018, 69(12): 5370–5378. DOI: 10.11949/j.issn.0438-1157.20180519.

    LI M, DU Y, LI G Q, et al. Experimental and large eddy simulation study on gasoline-air mixture explosions in semi-confined pipe with 90º right-angle bend [J]. CIESC Journal, 2018, 69(12): 5370–5378. DOI: 10.11949/j.issn.0438-1157.20180519.
    [28] 温小萍, 余明高, 邓浩鑫, 等. 小尺度受限空间内瓦斯湍流爆燃大涡模拟 [J]. 化工学报, 2016, 67(5): 1837–1843. DOI: 10.11949/j.issn.0438-1157.20151219.

    WEN X P, YU M G, DENG H X, et al. Large eddy simulation of gas turbulent deflagration in small-scale confined space [J]. CIESC Journal, 2016, 67(5): 1837–1843. DOI: 10.11949/j.issn.0438-1157.20151219.
    [29] 洪溢都, 林柏泉, 朱传杰. 开口型内瓦斯爆炸冲击波动压的数值模拟 [J]. 爆炸与冲击, 2016, 36(2): 198–209. DOI: 10.11883/1001-1455(2016)02-0198-12.

    HONG Y D, LIN B Q, ZHU C J. Simulation on dynamic pressure of premixed methane/air explosion in open-end pipes [J]. Explosion and Shock Waves, 2016, 36(2): 198–209. DOI: 10.11883/1001-1455(2016)02-0198-12.
  • 期刊类型引用(2)

    1. 杨克,李雪瑞,纪虹,郑凯,邢志祥,蒋军成. 改性煤矸石-海藻酸钠粉体对管道内甲烷/空气爆炸的抑爆实验. 爆炸与冲击. 2024(07): 174-187 . 本站查看
    2. 余明高,贺涛,李海涛,郑凯. 改性高岭土抑爆剂对瓦斯煤尘复合爆炸压力的影响. 煤炭学报. 2022(01): 348-359 . 百度学术

    其他类型引用(7)

  • 加载中
图(18) / 表(2)
计量
  • 文章访问数:  401
  • HTML全文浏览量:  111
  • PDF下载量:  57
  • 被引次数: 9
出版历程
  • 收稿日期:  2023-01-16
  • 修回日期:  2023-04-06
  • 网络出版日期:  2023-04-25
  • 刊出日期:  2023-08-31

目录

/

返回文章
返回