多种材质障碍物对甲烷-氢气预混燃气的促爆影响

焦一飞; 熊晓曼; 任昊; 米红甫; 何国钦; 李品; 魏鑫

doi:10.11858/gywlxb.20230682

多种材质障碍物对甲烷-氢气预混燃气的促爆影响

doi: 10.11858/gywlxb.20230682

焦一飞^1,,
熊晓曼¹,
任昊¹,
米红甫²,
何国钦^2,*, ,,
李品¹,
魏鑫¹

1.
国家电网四川省电力公司经济技术研究院, 四川成都　610041
2.
重庆科技学院安全工程学院, 重庆　401331

基金项目: 国家自然科学基金（52274177）

详细信息

作者简介:
焦一飞（1985－），男，硕士，高级工程师，主要从事输变电工程变电土建、消防专项评审及相关研究. E-mail：yifeijiao@hotmail.com

通讯作者:
何国钦（2000－），男，硕士研究生，主要从事油气爆炸动力学研究. E-mail：1923715297@qq.com

中图分类号: O389; X932
计量
- 文章访问数: 173
- HTML全文浏览量: 64
- PDF下载量: 26
出版历程
- 收稿日期: 2023-06-16
- 修回日期: 2023-08-28
- 录用日期: 2023-08-28
- 网络出版日期: 2024-02-05
- 刊出日期: 2024-02-05

Effect of Various Material Obstacles on the Promoting Explosion of Methane-Hydrogen Premixed Gas

JIAO Yifei^1
,,
XIONG Xiaoman¹,
REN Hao¹,
MI Hongfu²,
HE Guoqin^{2,*
, ,},
LI Pin¹,
WEI Xin¹

1.
State Grid Sichuan Electric Power Company Economic and Technological Research Institute, Chengdu 610041, Sichuan, China
2.
School of Safety Engineering, Chongqing University of Science and Technology, Chongqing 401331, China

摘要

摘要: 选取3种不同弯曲强度障碍物进行氢气-甲烷混合气体爆炸影响实验研究，以探究不同环境下的促爆危险性。实验过程中采集爆炸管道内火焰图像以及上下游压力。通过分析火焰图像以及爆炸压力数据后发现，爆炸后膨胀燃气推动流场加速在障碍物后方产生涡流，流场在不同材质障碍物后产生不同的涡流强度，导致燃气后期火焰峰值速度差异以及管道内爆炸超压差异。实验证明了促爆强度与障碍物材质的相关性。在本研究中，促爆强度与障碍物弯曲强度成正比，在氢气加入后，燃气基础反应加快，3种材质障碍物管道内爆炸压力峰值开始产生明显差异。通过上述结果可以得出，环境内障碍物以及粗糙壁面会影响燃气促爆效果，由于材料本身特性影响而产生差异，并且这种差异受燃气自身燃烧速率的影响。
- 燃气爆炸 /
- 障碍物 /
- 促爆实验 /
- 预混燃气
Abstract: Three obstacles with different levels of bending strength was selected for experimental research on the impact of hydrogen-methane mixed gas explosions in order to explore the varying environmental hazards of explosion promotion. During the experiment, the images of the flame in the explosion pipeline and the upstream and downstream pressure were collected. Through the analysis of flame images and explosion pressure data, the flow field accelerated by expanding gas after explosion generates eddy currents behind obstacles, and the flow field produces different eddy currents behind obstacles of different materials, which result in the difference of peak flame velocity of gas in the later stage and the difference of explosion overpressure in pipelines. This proves the correlation between the intensity of promoting explosions and the material of obstacles. In the experiments conducted in this paper, there is a proportional relationship between the intensity of promoting explosions and the bending strength of the obstacle. Moreover, after hydrogen was added, the reaction of the gas base was accelerated and the peak explosion pressure in the obstacle pipes of the three materials began to produce obvious differences. It can be concluded that obstacles and rough walls in the environment will affect the gas explosion effect, and the difference is caused by the characteristics of the material itself, and the difference is affected by the combustion rate of the gas itself.
- gas explosion /
- obstacle /
- promoting explosion experiment /
- premixed gas

HTML全文

1. 引言

水力喷射径向侧钻微小井眼技术是利用连续油管，通过特殊的转向设备在套管和水泥环上钻孔，进而在油藏的某一层位或多个层位沿径向钻出多个水平分支井眼，能够极大地提高单井产量，具有适应性强、钻井效率高、钻井成本低的优点^[1-6]。反喷牵引射流钻头具有向后开孔的反向喷嘴，形成的反向射流可为射流钻头提供牵引力，在钻进过程中有效地降低高压软管的屈曲失稳，减小软管进入地层的阻力，同时反向射流还能够清除岩屑，有利于增大钻孔孔径，增强钻头的稳定性^[7-9]。

准确计算反喷牵引射流钻头的自进力是确定径向水平井软管的预安装长度和实际延伸能力的关键因素^[10]。现有的射流钻头自进力计算一般只考虑流量的影响，其值为反向喷嘴牵引力与正向喷嘴反冲力之差^[11-12]。然而研究发现，反喷牵引射流钻头的实测自进力与计算值存在偏差。因此，有必要针对反喷牵引射流钻头自进力的主要影响因素进行研究，分析实测自进力与计算值产生偏差的内在机制。本研究以新型反喷牵引直旋混合射流钻头的结构为基础，实验研究入口流量、井筒直径、正向喷距等关键因素对射流钻头自进力的影响规律，并通过射流钻头流场结构的数值模拟分析3种因素对射流钻头自进力的影响机制，以期为建立多因素条件下射流钻头自进力的计算与控制方法提供基础。

2. 反喷牵引直旋混合射流钻头的结构与工作原理

根据射流的结构特征，目前应用于现场施工的反喷牵引射流钻头主要有反喷牵引多孔射流钻头和反喷牵引旋转射流钻头两种^[13]，其中：多孔射流钻头依靠多个喷嘴产生直射流进行破岩，能量集中，破岩能力强，但是其喷嘴多且直径小，喷嘴布置和安装角度对破岩能力和孔眼圆整性的影响较大，破岩和扩孔能力难以统一，并且对工作液的过滤要求高，容易发生喷嘴堵塞；单出口的旋转射流钻头则会产生具有较强扩散的旋转射流，虽有较强的扩孔能力，但有效喷距较短，在围压条件下破岩能力不足。

廖华林等人^[14-15]设计了一种兼备直射流和旋转射流特性的直旋混合射流钻头，并通过研究证明，直旋混合射流具有直射流和旋转射流的双重优点，比直射流的扩孔能力强，比旋转射流的钻孔深度大，可以通过调节不同组合的钻头结构参数改变其破岩特性，为现有的径向水平井射流钻头存在的问题提供了一种解决方案。本研究以直旋混合射流钻头为基础，在钻头体上增加反向喷嘴，设计出如图 1所示的反喷牵引直旋混合射流钻头。该射流钻头主要由钻头本体、反向喷嘴、带直孔的加旋叶轮、混合腔、喷嘴出口、扩展腔等部分组成。其基本工作原理为：来流进入钻头体后，一部分流体通过反向喷嘴向后喷射，为射流钻头和高压软管提供向前的牵引力，同时辅助扩孔和清除岩屑，另一部分流体则通过带直孔的加旋叶轮，形成向前喷射的直旋混合射流，高效破碎岩石。

图 1 反喷牵引直旋混合射流钻头的结构示意图

Figure 1. Schematic illustration of the self-propelled straight-swirling integrated jet bit

下载: 全尺寸图片幻灯片

3. 射流钻头自进力实验测试

3.1 实验装置

(1) 模拟井筒及实验台架(见图 2)。选用有机玻璃材质的模拟井筒，长度均为2 m，5种尺寸内径，分别为30、40、50、60和70 mm。实验台架的高度为0.5 m，长度为4 m，可用于固定模拟井筒。

图 2 实验装置示意图

Figure 2. Schematic illustration of experimental facilities

下载: 全尺寸图片幻灯片

(2) 射流钻头。为测试反喷牵引射流钻头反向喷嘴的牵引力和正向喷嘴的反冲力，设计并加工了只带反向喷嘴的牵引射流钻头和只带正向喷嘴的直旋混合射流钻头，如图 3、图 4所示。其中反向喷嘴的直径均为1.25 mm，与射流钻头轴线的夹角为20°，正向喷嘴的直径为2 mm，设计正反流量比为4:9。

图 3 牵引射流钻头示意图

Figure 3. Schematic illustration of backward jet bit

下载: 全尺寸图片幻灯片

图 4 直旋混合射流钻头示意图

Figure 4. Schematic illustration of straight-swirling integrated jet bit

下载: 全尺寸图片幻灯片

(3) 高压软管。选用金属编织式软管作为连接射流钻头的实验管线，软管长10 m，外径为18 mm，内径为10 mm，最高耐压强度为50 MPa。

(4) 高压泵。高压柱塞泵1台，额定工作压力为50 MPa，额定排量为63 L/min。

(5) 测试设备。测试设备包括数显推拉力计、压力传感器及高压流量计。

3.2 实验方法

将模拟井筒水平放置，固定于实验台架上。高压软管的一端连接射流钻头并置于模拟井筒内，另一端与拉力测试线相连。数显推拉力计固定于实验台架的尾部，泵压压力传感器连接在高压软管入口处的转换接头上，高压流量计连接于高压管线的入口处。

通过拉力传感器测得不同射流钻头的牵引力、反冲力和自进力。对于牵引射流钻头，因其只有反向喷嘴，拉力计所测的力为射流钻头牵引力F_q；对于直旋混合射流钻头，因其只有正向喷嘴，拉力计所测的力为射流钻头反冲力F_f；而对于反喷牵引直旋混合射流钻头，所测的力则为实际拖动软管的力，定义为自进力F_z。实验时，控制进入牵引射流钻头和直旋混合射流钻头的流量，使二者的流量之和等于进入反喷牵引直旋混合射流钻头的流量。由于牵引射流钻头和直旋混合射流钻头只有一种射流作用，所测的F_q和F_f无其他射流干扰，而反喷牵引射流钻头的自进力F_z在微小井眼内受到正、反射流的综合作用，因此将F_s(F_s=F_q-F_f)与F_z对比，分析实际自进力与计算值之间的差异。同时，改变射流钻头的入口流量Q、井筒直径D、正向喷距H等参数，分析其对反喷牵引射流钻头自进力的影响规律。实验参数列于表 1，每个参数组合测试2次，取其平均值作为实验结果。

表 1 实验及数值模拟参数

Table 1. Parameters of the experiment and numerical simulation

No.	Q/(L/min)	D/(mm)	H/(mm)
1^#	42.5, 48.5, 50.1, 56.1, 60.0	50	20
2^#	56.1	30, 40, 50, 60, 70	20
3^#	56.1	50	10, 20, 30, 40, 50

下载: 导出CSV

| 显示表格

4. 射流钻头流场的数值分析

4.1 控制方程和离散方法

根据表 1的参数建立物理模型，采用RNG湍流模型进行流场模拟^[16]。井眼内射流钻头的三维数值分析模型如图 5所示。各控制方程中，对流项和扩散项的离散采用混合格式，对源项实施负斜率线性化。对于各变量的离散化方程组，采用TDMA(Tri-Diagonal Matrix Algorithm)逐线扫描低松弛迭代法求解。连续方程和动量方程的联立耦合迭代求解则采用SIMPLEC算法。

图 5 井眼内反喷牵引直旋混合射流钻头的三维数值分析模型图

Figure 5. Three dimensional numerical simulation model of self-propelled straight-swirling integrated jet bit in the well

下载: 全尺寸图片幻灯片

4.2 边界条件

入口边界：入口速度取来流速度，湍流脉动动能k及湍流耗散率ε为

$I = 0.16\mathit{R}{\mathit{e}^{\mathit{ - }{\rm{1/8}}}}$

(1)

$k = 1.5{\left( {uI} \right)^2}$

(2)

$\varepsilon = C_\mu ^{3/4}{k^{1.5}}/l$

(3)

式中：Re为雷诺数；u为入口速度；C_μ为常量；I为湍流强度；l=0.07L，其中L为管道特征尺寸。

出口边界：出口压力为0.1 MPa。

壁面条件：采用无滑移固壁边界。

5. 结果分析

5.1 入口流量对射流钻头自进力的影响

图 6给出了井筒直径D=50 mm、正向喷距H=20 mm时入口流量对自进力的影响。由图 6(a)可知，当入口流量Q由42.5 L/min增加至60.0 L/min时，射流钻头牵引力F_q由43 N增至96 N，反冲力F_f由21 N增至57 N，自进力F_z由6 N增至40 N，并且均与入口流量近似呈线性关系。因而，在条件许可的情况下，可通过增大排量提高射流钻头的自进力。分析认为，射流钻头的自进力随入口流量增加的原因是以下3方面综合作用的结果。(1)反向喷嘴的流量大于正向喷嘴，随着入口流量的增加，牵引力和反冲力均呈近线性化增加，故射流钻头的自进力增大。(2)由数值分析结果可知(见图 6(b)，其中p_a为环空压力)，反向喷嘴喷出的高速流体对周围流体形成卷吸，使反向喷嘴出口的局部压力降低，多股反向射流共同作用形成射流封隔，从而在钻头后部产生低于环境压力的低压区，环境压力与低压区的压差作用提高了射流钻头的自进力。随着入口流量的增加，卷吸效果和射流封隔能力增强，钻头后部的低压区压力逐渐降低，有利于增加射流钻头的自进力^[11]。(3)正向喷出的直旋混合射流具有径向和切向速度，在冲击井底的过程中也对周围流体产生卷吸作用，如果环空间隙小，外界流体无法补充，则会在钻头前段形成环形低压区，负压推力作用促进了钻头的自进。从图 6(c)可以看出，随着入口流量的增加，钻头前端的径向速度(v_r)增加，卷吸效果增强，也有利于增大射流钻头的自进力。

图 6 入口流量对自进力的影响

Figure 6. Influence of inlet flow rate onthe self-propelled force

下载: 全尺寸图片幻灯片

实验测试发现，射流钻头的实测自进力F_z与牵引力和反冲力的差值F_s存在差异。当入口流量为42.5 L/min时，测得F_z为6 N，F_s为22 N，F_z远小于F_s；随着入口流量的增加，F_z逐渐逼近F_s；当入口流量为60.0 L/min时，F_z与F_s基本相同，约40 N。分析认为，反喷牵引直旋混合射流钻头的自进力受入口流量、实际正反流量比以及返流强度等多种因素的综合影响，其中反向喷嘴出口处的降压效应在不同流量条件下存在差异，使射流钻头附近的流动更为复杂，实测自进力F_z与计算值F_s之间的差异不同。因此，在计算反喷牵引直旋混合射流的自进力时，需要综合考虑多种因素的影响。

5.2 井筒直径对射流钻头自进力的影响

图 7给出了入口流量Q=56.1 L/min、正向喷距H=20 mm时井筒直径对自进力的影响。从图 7(a)中可以看出，当井筒直径D从30 mm增加至70 mm时，射流钻头的牵引力F_q由93 N降至85 N，反冲力F_f由47 N降至41 N，而自进力F_z则由27 N先增至33 N后降至19 N并稳定。从数值分析结果(见图 7(b)和图 7(c))可以看出：当井筒直径较小时，正向射流冲击壁面形成的漫流会迅速转化为返流，形成的返流流速较高，对钻头前端产生较大的后推力，阻碍钻头自进；随着井筒直径的增加，井底漫流形成的返流流速变缓，对射流钻头前端面的作用减弱，射流钻头的自进力增大；但是在入口流量不变的情况下，反向射流的封隔能力有限，当井筒直径超过反向射流封隔范围时，外部流体对被射流抽吸掉的部分进行补充，从而导致封隔区内、外基本无压差，封隔失效，射流钻头的自进力降低；当井筒直径继续增大时，射流的封隔作用失效，射流钻头的自进力基本无变化。由此可见，存在最优的井筒直径范围，使射流钻头的自进力最大，本实验中射流钻头的最优井筒直径范围为35~45 mm。

图 7 井筒直径对自进力的影响

Figure 7. Influence of wellbore diameter onthe self-propelled force

下载: 全尺寸图片幻灯片

5.3 正向喷距对射流钻头自进力的影响

图 8给出了入口流量Q=56.1 L/min、井筒直径D=50 mm时正向喷距对射流钻头自进力的影响。从图 8(a)中可以看出，正向喷距H从10 mm增加至50 mm时，射流钻头的牵引力F_q基本稳定在93 N左右，反冲力F_f约为45 N，自进力F_z约为27 N，即随着正向喷距的增大，射流钻头的自进力基本不发生改变。分析认为，正向喷距影响射流钻头自进力的方式主要来自于喷距对返流流速和正向射流卷吸效果所产生的影响。从数值模拟结果(见图 8(b)和图 8(c))可以看出：当正向喷距较小时，形成的返流速度较大，对钻头前端产生较大的向后推力，阻碍射流钻头自进，但是较小的喷距也使直旋混合射流冲击井底壁面前的能量耗散较小，对周围流体的卷吸效果增强，促进了射流钻头的自进；随着正向喷距的增加，返流形成的向后推力降低，同时直旋混合射流在冲击壁面前的能量消耗也增加，卷吸降压效果下降，两种作用相互抵消后，对自进力的影响也降到很低，基本可以忽略不计。

图 8 正向喷距对自进力的影响

Figure 8. Influence of spray distance on the self-propelled force

下载: 全尺寸图片幻灯片

6. 结论

(1) 随着入口流量的增大，反喷牵引直旋混合射流钻头的自进力增大。正向射流和反向射流的卷吸降压效果随着入口流量的增加而增强，有利于增大射流钻头的自进力。

(2) 随着井筒直径的增大，反向射流的封隔效果和返流干扰减弱，致使反喷牵引直旋混合射流钻头的自进力先增加后降低。在入口流量相同的条件下，存在最优的井筒直径范围，当入口流量为56.1 L/min、正向喷距为20 mm时，最优的井筒直径范围为35~45 mm。

(3) 反喷牵引直旋混合射流钻头的自进力受入口流量、井筒直径、实际正反流量比、返流强度等多种因素的综合影响，现有的自进力计算方法仅考虑了流量的影响，而自进力的实测值与计算值在不同条件下存在差异，因此需要建立考虑多因素影响的自进力计算方法。

图 1 实验系统示意图

Figure 1. Schematic diagram of experimental system

下载: 全尺寸图片幻灯片

图 2 3种障碍物：刚性障碍物、柔性障碍物A和柔性障碍物B

Figure 2. Three types of obstacle：rigid obstacle，flexible obstacle A and flexible obstacle B

下载: 全尺寸图片幻灯片

图 3 管道内火焰图像（ $\varphi_{{\mathrm{H}}_2}$ =0）

Figure 3. Images of flame in pipeline ( $\varphi_{{\mathrm{H}}_2}$ =0)

下载: 全尺寸图片幻灯片

图 4 管道内火焰图像（ $\varphi_{{\mathrm{H}}_2}$ =10%）

Figure 4. Images of flame in pipeline ( $\varphi_{{\mathrm{H}}_2}$ =10%)

下载: 全尺寸图片幻灯片

图 5 管道内火焰图像（ $\varphi_{{\mathrm{H}}_2}$ =20%）

Figure 5. Images of flame in pipeline ( $\varphi_{{\mathrm{H}}_2}$ =20%)

下载: 全尺寸图片幻灯片

图 6 管道内火焰图像（ $\varphi_{{\mathrm{H}}_2}$ =30%）

Figure 6. Images of flame in pipeline ( $\varphi_{{\mathrm{H}}_2}$ =30%)

下载: 全尺寸图片幻灯片

图 7 火焰锋面速度-位置曲线

Figure 7. Curves of flame front velocity-position

下载: 全尺寸图片幻灯片

图 8 火焰锋面的速度-位置曲线

Figure 8. Curves of flame front velocity-position

下载: 全尺寸图片幻灯片

图 9 $\varphi_{{\mathrm{H}}_2}$ =0时管道内上游的压力-时间曲线

Figure 9. Pipeline upstream pressure-time curves of $\varphi_{{\mathrm{H}}_2}$ =0

下载: 全尺寸图片幻灯片

图 10 $\varphi_{{\mathrm{H}}_2}$ =0时管道内下游压力-时间曲线

Figure 10. Pipeline downstream pressure-time curves of $\varphi_{{\mathrm{H}}_2}$ =0

下载: 全尺寸图片幻灯片

图 11 管道内上游的压力峰值-氢气体积分数曲线

Figure 11. Upstream pressure peak-hydrogen volume fraction curves in pipeline

下载: 全尺寸图片幻灯片

图 12 管道内下游的压力峰值-氢气体积分数曲线

Figure 12. Downstream pressure peak-hydrogen volume fraction curves in pipeline

下载: 全尺寸图片幻灯片

表 1 预混气体的组分

Table 1. Premixed gas components

${\varphi _{{{\mathrm{H}}_2}}}$ /%	${\varPhi _{{{\mathrm{H}}_2}}}$ /%	${\varPhi _{{{\mathrm{CH}}_4}}}$ /%	$\varPhi _{{\mathrm{Air}}}$ /%
0	0	9.50	90.50
10	1.02	9.18	89.80
20	2.20	8.80	89.00
30	3.58	8.53	88.07

下载: 导出CSV

表 2 障碍物参数

Table 2. Obstacle parameter

Obstacle	Material	Dimensions	Mean density/ (g·cm⁻³)	Average elastic modulus/MPa
Rigid obstacle	Polymethyl methacrylate	100 mm×10 mm × 30 mm	1.190	2800.0
Flexible obstacle A	Foamed silicone gel		1.920	31.3
Flexible obstacle B	Polyurethane foam		0.699	7.4

下载: 导出CSV

参考文献(29)

[1]	STARR A, LEE J, NG H D. Detonation limits in rough walled tubes [J]. Proceedings of the Combustion Institute, 2015, 35(2): 1989–1996. doi: 10.1016/j.proci.2014.06.130
[2]	ShCHELKIN K I, TROSHIN Y K. Non-stationary phenomena in the gaseous detonation front [J]. Combustion and Flame, 1963, 7: 143–151. doi: 10.1016/0010-2180(63)90172-X
[3]	SHCHELKIN K I. Instability of combustion and detonation of gases [J]. Soviet Physics Uspekhi, 1966, 8(5): 780–797. doi: 10.1070/PU1966v008n05ABEH003038
[4]	MOEN I O, DONATO M, KNYSTAUTAS R, et al. Flame acceleration due to turbulence produced by obstacles [J]. Combustion and Flame, 1980, 39(1): 21–32. doi: 10.1016/0010-2180(80)90003-6
[5]	NA’INNA A M, PHYLAKTOU H N, ANDREWS G E. The acceleration of flames in tube explosions with two obstacles as a function of the obstacle separation distance [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6): 1597–1603. doi: 10.1016/j.jlp.2013.08.003
[6]	FAIRWEATHER M, HARGRAVE G K, IBRAHIM S S, et al. Studies of premixed flame propagation in explosion tubes [J]. Combustion and Flame, 1999, 116(4): 504–518. doi: 10.1016/S0010-2180(98)00055-8
[7]	YU M G, ZHENG K, CHU T K. Gas explosion flame propagation over various hollow-square obstacles [J]. Journal of Natural Gas Science and Engineering, 2016, 30: 221–227. doi: 10.1016/j.jngse.2016.02.009
[8]	WEN X P, YU M G, JI W T, et al. Methane-air explosion characteristics with different obstacle configurations [J]. International Journal of Mining Science and Technology, 2015, 25(2): 213–218. doi: 10.1016/j.ijmst.2015.02.008
[9]	WEN X P, YU M G, LIU Z C, et al. Large eddy simulation of methane-air deflagration in an obstructed chamber using different combustion models [J]. Journal of Loss Prevention in the Process Industries, 2012, 25(4): 730–738. doi: 10.1016/j.jlp.2012.04.008
[10]	MASRI A R, IBRAHIM S S, NEHZAT N, et al. Experimental study of premixed flame propagation over various solid obstructions [J]. Experimental Thermal and Fluid Science, 2000, 21(1/2/3): 109–106. doi: 10.1016/S0894-1777(99)00060-6
[11]	JOHANSEN C T, CICCARELLI G. Modeling the initial flame acceleration in an obstructed channel using large eddy simulation [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(4): 571–585. doi: 10.1016/j.jlp.2012.12.005
[12]	JOHANSEN C T, CICCARELLI G. Visualization of the unburned gas flow field ahead of an accelerating flame in an obstructed square channel [J]. Combustion and Flame, 2009, 156(2): 405–416. doi: 10.1016/j.combustflame.2008.07.010
[13]	CICCARELLI G, JOHANSEN R T, PARRAVANI R. The role of shock-flame interactions on flame acceleration in an obstacle laden channel [J]. Combustion and Flame, 2010, 157(11): 2125–2136. doi: 10.1016/j.combustflame.2010.05.003
[14]	SALAMANDRA G D, BAZHENOVA T V, NABOKO I M. Formation of detonation wave during combustion of gas in combustion tube [J]. Symposium (International) on Combustion, 1958, 7(1): 851–855. doi: 10.1016/S0082-0784(58)80128-9
[15]	ZHANG B, LIU H, LI Y C. The effect of instability of detonation on the propagation modes near the limits in typical combustible mixtures [J]. Fuel, 2019, 253: 305–310. doi: 10.1016/j.fuel.2019.05.006
[16]	SULAIMAN S Z, KASMANI R M, MUSTAFA A, et al. Effect of obstacle on deflagration to detonation transition (DDT) in closed pipe or channel-an overview [J]. Jurnal Teknologi, 2013, 66(1): 49–52. doi: 10.11113/jt.v66.1326
[17]	OGAWA T, GAMEZO V N, ORAN E S. Flame acceleration and transition to detonation in an array of square obstacles [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(2): 355–362. doi: 10.1016/j.jlp.2011.12.009
[18]	ZHANG B, LIU H, YAN B J. Investigation on the detonation propagation limit criterion for methane-oxygen mixtures in tubes with different scales [J]. Fuel, 2019, 239: 617–622. doi: 10.1016/j.fuel.2018.11.062
[19]	BANG B H, AHN C S, KIM Y T, et al. Deflagration-to-detonation transition in pipes: the analytical theory [J]. Applied Mathematical Modelling, 2019, 66: 332–343. doi: 10.1016/j.apm.2018.09.023
[20]	KIVERIN A D, YAKOVENKO I S. Estimation of critical conditions for deflagration-to-detonation transition in obstructed channels filled with gaseous mixtures [J]. Mathematical Modelling of Natural Phenomena, 2018, 13(6): 54. doi: 10.1051/mmnp/2018071
[21]	COATES A M, MATHIAS D L, CANTWELL B J. Numerical investigation of the effect of obstacle shape on deflagration to detonation transition in a hydrogen-air mixture [J]. Combustion and Flame, 2019, 209: 278–290. doi: 10.1016/j.combustflame.2019.07.044
[22]	LEAL C A, SANTIAGO G F. Do tree belts increase risk of explosion for LPG spheres? [J]. Journal of Loss Prevention in the Process Industries, 2004, 17(3): 217–224. doi: 10.1016/j.jlp.2004.02.003
[23]	BAKKE J R, WINGERDEN K V, HOORELBEKE P, et al. A study on the effect of trees on gas explosions [J]. Journal of Loss Prevention in the Process Industries, 2010, 23(6): 878–884. doi: 10.1016/j.jlp.2010.08.007
[24]	LI Q, LU S X, XU M J, et al. Comparison of flame propagation in a tube with a flexible/rigid obstacle [J]. Energy & Fuels, 2016, 30(10): 8720–8726. doi: 10.1021/acs.energyfuels.6b01594
[25]	LI Q, CICCARELLI G, SUN X X, et al. Flame propagation across a flexible obstacle in a square cross-section channel [J]. International Journal of Hydrogen Energy, 2018, 43(36): 17480–17491. doi: 10.1016/j.ijhydene.2018.07.077
[26]	LI Q, SUN X X, WANG X, et al. Experimental study of flame propagation across flexible obstacles in a square cross-section channel [J]. International Journal of Hydrogen Energy, 2019, 44(7): 3944–3952. doi: 10.1016/j.ijhydene.2018.12.085
[27]	李权. 管道内障碍物对氢-空气预混火焰传播动力学影响研究 [D]. 合肥: 中国科学技术大学, 2019.
[28]	YU M G, ZHENG K, ZHENG L G, et al. Effects of hydrogen addition on propagation characteristics of premixed methane/air flames [J]. Journal of Loss Prevention in the Process Industries, 2015, 34: 1–9. doi: 10.1016/j.jlp.2015.01.017
[29]	YU S W, DUAN Y L, LONG F Y, et al. The influence of flexible/rigid obstacle on flame propagation and blast injuries risk in gas explosion [J]. Energy Sources, Part A: Recovery, Utilization, and Environmental Effects, 2023, 45(2): 4520–4536.