On the distribution of explosion strain field and fracture field in segment charge
-
摘要: 为探究分段装药爆炸应变场与裂隙场分布规律,采用数字图像相关分析方法与电子计算机断层扫描实验方法,分析了孔内分段装药爆炸全场应变传播规律,建立了爆后“岩石—爆炸裂隙”的三维重构模型,描述了爆炸裂纹位置与形态的空间分布情况,得到岩石材料爆炸裂隙的分形维数与损伤度。研究结果表明:分段装药改变了连续装药对介质的全场应变形态,由一次应变改变为两次应变,在满足第一段炸药对介质的破坏作用下,同时加大了第二段炸药对介质的作用效应;上分段装药占比0.4时,下分段介质受爆炸作用应变峰值更大,更好满足工程实践中下半段岩体对爆炸能量的需求;相同装药系数下,连续装药结构爆炸裂纹没有贯穿试件整体,炮孔封堵段的爆炸裂纹较少,分段装药结构下,由于提高了炸药的位置,使得上部分岩体能够更好地利用炸药爆炸的能量破碎岩石;分段装药岩石整体损伤度较连续装药提高了23.5%,其中上分段岩石损伤差异较大,分段装药上分段损伤度比连续装药提高46.4%。Abstract: In order to explore the distribution of the explosion strain field and fracture field of segmented charge explosion, used digital image correlation analysis (DIC) and computerized tomography (CT) scanning experiment analyzed the distribution of the explosion strain field and fracture field of the segment charge, established three-dimensional reconstruction model of “rock-explosion crack”, described the spatial distribution of the location and shape of the explosion crack, and obtained the fractal dimension and damage degree of the explosion crack. The research results show that: the segment charge changes the full-field strain morphology of the medium caused by continuous charge. Under the condition of satisfying the damage effect of upper sublevel explosive on medium, the effect of lower sublevel explosive on medium is increased, at the same time, the time of blast stress wave is prolonged. It can be seen from the three-dimensional cracks of segmented charge and continuous charge explosion, the explosion crack mainly expands along the radial direction, the annular crack formed by axial stress and strain is not obvious, and the radial direction is the main direction of rock failure. When the charge ratio of the upper segment is 0.4, the strain peak value of the lower segment is larger, which better meets the demand of the rock mass for explosion energy in the lower segment in engineering practice. Under the same charging coefficient, the explosive cracks in the continuous charging structure do not run through the whole specimen, and the explosive cracks in the plugging section are less, under the segment charge structure, the upper sublevel of rock mass can better use the energy of explosive explosion to break rock because the position of explosive is improve. The overall damage degree of rock in segment charge is 23.5% higher than that in continuous charge, and the damage degree of rock in upper sublevel is 46.4% higher than that in continuous charge.
-
Key words:
- explosion /
- segment charge /
- full-field strain /
- fracture field /
- three dimensional reconstruction
-
油气是生产、生活中必不可少的能源物质,也是典型的易燃易爆物。如何在油气储存、运输、使用的过程中预防或控制火灾、爆炸事故的发生,是被长期关注的问题。油气火灾爆炸事故发生的初始阶段,通常体现为局部油蒸汽-空气混合气着火,发生以亚音速传播、并伴随有压力波的燃烧反应(即爆燃)[1-2],而这个阶段,正是采取有效抑制措施、防止事故扩大的关键阶段[3]。
火焰的形成是可燃混合气爆燃过程的典型外部特征之一,火焰形态即是灾害预警、燃烧模式识别的重要依据,也是探索爆燃发展支配机理、变化规律与分析模型的重要资料。针对燃烧爆炸过程的火焰形态,有了广泛的研究。A.Yoshida等[4]通过火焰图像观察了未燃气体中的湍流尺度对预混火焰褶皱尺度的影响,K.H.Oh等[5]基于火焰形态分析了障碍物对液化石油气爆炸火焰传播的影响。P.Zhang等[6]对坑道内的油气爆燃进行了实验研究,采集了爆燃火焰传播过程中局部的火焰图像,但未对火焰的形态结构进行探讨。M.Zhang等[7]研究了在开敞空间CO/H2/O2/CO2混合气湍流预混火焰的火焰锋面结构,测量了火焰厚度、平均火焰体积等参数,并指出具有褶皱的火焰锋面是湍流预混火焰的普遍特征。S.Jerzembeck等[8]指出,在高于当量浓度时,爆燃火焰会产生褶皱和胞状结构,但他们专注于火焰速度的测量,并未针对火焰特征进行详细阐述。
本文中,基于油气爆燃可视化实验,提出受限空间油气爆燃的3种火焰形态,定量给出区分不同火焰形态的临界判据,分析不同火焰形态的产生机理,并探讨不同火焰形态下油气爆燃特征参数的差异。
1. 实验系统
图 1为实验系统示意图,实验台架主体为四周密闭、水平放置的长方体,内部尺寸为500 mm×200 mm×200 mm,正面设有可视化窗口,上加装钢化玻璃。实验台架顶板设置传感器接口,与压力传感器、火焰强度传感器连接,通过数据采集系统实时记录油气爆燃过程中的瞬时压力与火焰亮度变化;底板中央布置点火装置,点火头伸入实验台架内部50 mm,点火电平为1 500 V,点火能为3 J。压力与火焰数据的采集采用成都泰斯特电子信息有限责任公司的TST6300多用动态信号测试仪。采用自制的雾化装置(包括变径管和油瓶)将油气混入空气中形成可燃混合气。混合气组分的测定采用北京均方理化科技研究所的GXH-1050型红外线分析仪。图像采集采用日本Photron公司的FAST CAM Ultima 512型高速摄影仪。实验所用油品为中石化93号汽油。
2. 不同油气浓度下油气爆燃火焰形态
图 2为油气浓度(混合气中碳氢化合物的体积分数)YCH=1.50%时油气爆燃火焰的时序图像,拍摄频率为250 s-1,图中可视化窗口的实际尺寸为300 mm×200 mm,中部下方黑色柱状体为点火头,点火位置即位于其尖端。当t=0时,点火器启动,电火花释放出的能量在它附近形成高能点火气云,并引燃周围可燃气,形成光滑的淡蓝色椭球形火焰锋面(8~12 ms)。之后,反应区迅速扩大(12~21 ms),此时的火焰形态类似于层流预混火焰,但通过火焰传播至容器边缘所用的时间可以估算,其火焰速度高于层流预混火焰[1]一个数量级。火焰到达容器壁面(28 ms)后,反应区内发出亮白色光(34~42 ms),直至火焰完全熄灭。
图 3为YCH=2.33%时的油气爆燃过程。点火初期(8 ms),点火点附近可见与图 2类似的点火气云和球形火焰锋面,但随着反应的传播,火焰不断产生褶皱(18~32 ms),火焰锋面表面可见凹凸不平的鱼鳞状波纹(37 ms)。在火焰传播过程中,反应区中下方可观察到细小的黄色次生火焰点(26~37 ms)。火焰锋面到达容器边界后,除部分凝结在壁面上的组分仍在燃烧,整个反应区呈现暗黄色,直至反应完全熄灭。
图 4为YCH=2.90%时的油气爆燃过程。与图 2~3比较,图 4中观察不到鲜明的火焰包络面,反应区为卷曲的絮状黄白色火焰。由于此时油气浓度较高,氧气相对不足,爆燃反应相对缓慢,火焰传播速度取决于传质速度,因而更多的体现出扩散燃烧的特征。混合气被点燃(8 ms)后,火焰首先向上方蔓延(32 ms),并形成一条条火焰束,顶部弯向四周未燃区域(53~68 ms)。随着反应的进行,火焰束在容器中上方交织,形成亮白色反应区,并最终扩散到整个容器。之后,火焰自上而下逐渐熄灭。
3. 3种典型火焰形态的临界条件与形成机理
对比图 2~4可以看出,不同浓度下的油气爆燃火焰形态具有显著的差异,然而,在一定的浓度变化范围内,火焰形态又具有一定的相似性。图 5给出了6个不同浓度下爆燃反应的典型火焰图像,可以看出,每列上下两幅图中火焰形态比较相似,而列与列之间的火焰形态则有明显区别。具体而言,从左第1列火焰锋面的形状为光滑的球面(图片中呈弧形),整个流场被火焰锋面清晰地分为已燃区和未燃区,这类火焰称为光滑球形火焰。第2列中,流场仍然被火焰锋面区分为已燃区和未燃区,但火焰锋面并不光滑,而是呈带有明显褶皱的球面,存在大量火焰褶皱形成的鱼鳞状波纹,这类称为褶皱球形火焰。第3列中,已观察不到鲜明的火焰包络面,火焰以絮状传播并向四周弯曲,这类称为卷曲絮状火焰。
进一步的实验从低于爆炸下限的浓度开始,逐步提高油气浓度,并观察其火焰特征、归纳火焰形态,直至浓度高于爆炸上限。结果为:当YCH=1.20%, 1.48%, 1.50%, 1.65%, 1.70%时,为光滑球形火焰;当YCH=2.00%, 2.15%, 2.33%时,为褶皱球形火焰;当YCH=2.75%, 2.90%, 3.25%时,为卷曲絮状火焰;当YCH=0.95%, 3.70%时,未点燃。可见,在爆炸极限范围内,不同浓度的油气爆燃火焰可概括为上述3种模式:油气浓度低于1.70%时,为光滑球形火焰;油气浓度在2.00%~2.33%之间时,为褶皱球形火焰;在油气浓度大于2.75%时,为卷曲絮状火焰。
3种火焰形态对应着3个油气浓度区间,而划分3个区间的临界浓度,可以从油气燃烧的化学当量进行理论分析。由于汽油是由近百种烃类组成的复杂混合物,难以确定具体的化学组分,通常的方法是采用异辛烷、PRF(primary reference fuel)[8]等与汽油燃烧特性近似、且具有确定化学组分的化合物或简单混合物作为替代燃料。以异辛烷为例,根据化学反应方程式可以计算:它在空气中恰好完全燃烧,生成产物完全为CO2和H2O时,对应的体积浓度为1.68%(取空气中氧浓度为21%);生成产物完全为CO和H2O时,对应的体积浓度为2.47%。这两个浓度与上述火焰形态发生变化的临界油气体积分数相近。通常定义燃料的实际浓度与其当量浓度的比为燃料当量比φ,根据上文的讨论,对于油气,取临界当量比φCH, 1=1、φCH, 2=2.47/1.68=1.47。
当φ < 1时,爆燃在富氧条件下进行,混合气被点燃后,由于混合气中氧含量充足,火焰不受氧化剂浓度的制约,向各个方向的传播速度相等;同时,油气浓度较低,反应较弱,对流场的扰动相对较小,因此形成了光滑的球形火焰。当1<φ < 1.47时,爆燃在贫氧条件下进行,由于局部耗/含氧量的不均匀,反应进度在空间上出现差异,火焰锋面不再光滑;同时,油气浓度的增加使反应更剧烈,湍流作用进一步加剧了火焰面的形变,因此形成了褶皱球形火焰。S.Jerzembeck等[8]针对合成气的研究指出:在低于当量浓度时,湍流预混火焰表面是光滑的,而在等于和略大于当量浓度时,合成气湍流预混火焰在点燃后将迅速出现褶皱和胞状结构。本文中所得结果与此一致。当φ>1.47时:爆燃在严重贫氧条件下进行,倾向于发生析碳反应,即CmH2n→mC+nH2,点火头最初的点火能激发了周边的可燃气,释放能量并形成高能活性基团;由于氧气严重不足,反应速度较缓慢,反应区气体升温后,向混合气上方对流,处于其中的剩余可燃气则与上方的氧气混合,继续进行反应;当正上方的氧气不足时,还原性基团向侧上方“寻氧”。这样,就形成了卷曲絮状火焰。
4. 火焰形态与爆燃关键参数的关系
图 6为燃烧产物与油气当量比的关系。可以看出,在不同的火焰形态下,油气爆燃的产物组成显著不同,具体如下。
(1) 当φ≈0.7时,油气相对含量较少,燃烧几乎不产生CO,O2大量剩余。之后,当φ < 1时,随着油气浓度增加,可燃混合气中的碳元素总量增加,O2仍然富余,生成越来越多的CO2;同时,燃烧的不完全程度增加,燃烧产物中的CO含量持续上升。
(2) 当1 < φ < 1.47时,随着可燃物增加,反应速度加剧,加之湍流混合作用变得显著,导致了局部混合气的反应并不充分,O2剩余量增加;随着φ的增大,产物中CO含量不断上升,CO2含量不断下降,这是因为空气含氧量是有限的,油气的增加使得更多的碳不能被充分氧化。应当考虑到,由于O2的大量剩余,在此阶段油气组分的析碳反应已经出现,并相应的随着φ的增大,将生成更多的碳颗粒。
(3) 当φ>1.47时,反应产物中CO与CO2几乎不再变化,这是因为此时氧含量已经严重不足,更多的可燃物只能发生析碳反应,生成更多的碳颗粒。而由于化学反应链竞争的存在,总有部分CO基团能够成功捕获O,生成一定的CO2;产生的碳颗粒也可能继续与O2发生反应,延长了反应的时间,也使得O2进一步被消耗。
图 7为爆燃升压速率ṗ、最大火焰强度I与油气当量比的关系,并通过两个临界当量比(图中虚线)区分了不同火焰形态。当φ < 1时,由于氧含量充足,油气浓度提高增加了可燃混合气中可燃物含量,使单位时间内参与燃烧反应的活性物质数量增加,火焰强度与升压速率因而增大。当1<φ < 1.47时,在油气适度过量的情况下,少部分油气与氧气反应生成CO,加速了反应进程,但在相同耗氧量下降低了释放出的总能量,因此升压速率增加,而火焰强度减弱。W.K.Kim等[9]的研究表明,富燃气条件下湍流混合作用将引起火焰褶皱加速,反应速度将大于富氧条件。本文中所得结果与此一致。当φ>1.47时,随着油气的进一步增加,氧含量相对不足,析碳反应产生的碳颗粒O2的反应延长了总反应时间,大幅降低了升压速率,却促进了氧气的充分燃烧,增加了空间内的热量积累,从而使得火焰强度增强。
图 8为爆燃压力峰值pm、火焰持续时间τ与油气当量比的关系。可以看出,压力峰值随油气当量比的增大呈先增后减的变化规律,这与升压速率的变化是一致的。同时,火焰持续时间随油气当量比的增大呈先减后增的变化规律,即:在油气浓度略高于当量浓度的褶皱球形火焰形态下,爆燃反应速率最快,破坏性也最强。C.Tang等[10]针对天然气的研究指出,当可燃气当量比为1.1时,火焰持续时间最短,压力峰值最大。本文中数据与此十分接近。
然而,从图 5中可以看出,当φ>1.6时,火焰持续时间(约600~900 ms)远大于φ < 1时(约250~400 ms)。这点本文中结果与文献[10]的研究结果并不一致。这主要原因是:相对于文献[10]中的天然气,油气中重质成分含量更多,燃烧反应需要经过更多的步骤,反应速度受贫氧条件的制约更显著。
图 6~8说明,火焰形态与油气爆燃的反应产物、升压速率、火焰强度、压力峰值、火焰持续时间等诸多关键参数密切相关,在火焰形态发生变化的临界油气当量比前后,上述参数随当量比的变化关系往往发生改变。对3种典型火焰形态下的油气爆燃特征进行了汇总,见表 1。
表 1 3种典型火焰形态下的油气爆燃特征变化规律Table 1. Gasoline-air mixture deflagration characteristicsof three different flame patterns爆燃参数 光滑球形火焰 褶皱球形火焰 卷曲絮状火焰 φ < 1 1~1.47 >1.47 YCO 增加 增加 最大 YCO2 增加后减小 减小 最小 YO2 减小 增加后减小 最小 I 增加 减小 增加 ṗ 增加 最大 减小 τ 减小 最小 增加 pm 增加 最大 减小 注:增加或减小均指随初始油气体积分数增加而发生的变化。 5. 结论
通过受限空间油气爆燃可视化实验研究,提出了受限空间油气爆燃的3种火焰形态,通过实验数据与理论计算给出了区分不同火焰形态的临界判据,分析了3种火焰形态产生的机理,并探讨了3种火焰形态下油气爆燃特征参数的差异。主要结论如下。
(1) 初始油气体积分数的不同,使油气爆燃体现出3种典型的火焰形态,即光滑球形火焰、褶皱球形火焰和卷曲絮状火焰,不同火焰形态的形成源自于油气与空气不同配比下燃烧反应机理的差异。
(2) 临界油气当量比φCH, 1=1、φCH, 2=1.47,可作为3种火焰类型的划分依据。
(3) 火焰形态与油气爆燃的关键参数密切相关:对光滑球形火焰,最大火焰强度、升压速率和最大压力随当量比的增加而增加,火焰持续时间随当量比的增加而减小;对褶皱球形火焰,最大火焰强度随当量比的增加而减小,升压速率和最大压力达到其最大值,而火焰持续时间达到最小值;对卷曲絮状火焰,最大火焰强度和火焰持续时间随当量比的增加而增加,升压速率和最大压力随当量比的增加而减小。
-
-
[1] MA X M, WANG Y, LIN T S, et al. Intelligent decision system for vertical shaft blasting scheme based on knowledge base and its application [J]. IEEE Access, 2021, 9: 163831–163842. DOI: 10.1109/ACCESS.2021.3128550. [2] MA Q Y, YUAN P, ZHANG J S, et al. Blast-induced damage on millisecond blasting model test with multicircle vertical blastholes [J]. Shock and Vibration, 2015, 2015: 504043. DOI: 10.1155/2015/504043. [3] 颜事龙, 冯叔瑜, 金孝刚. 有机玻璃中条形药包爆炸破碎区发展过程的高速摄影研究 [J]. 爆炸与冲击, 1998, 18(3): 262–267.YAN S L, FENG S Y, JIN X G. Study on fracture region of cylindric exploding in PMMA by high speed photograph [J]. Explosion and Shock Waves, 1998, 18(3): 262–267. [4] 龚敏, 黎剑华. 延长药包不同位置起爆时的应力场 [J]. 北京科技大学学报, 2002, 24(3): 248–253. DOI: 10.13374/j.issn1001-053x.2002.03.047.GONG M, LI J H. A research on stress field of column and strip-shaped charge in different detonated points [J]. Journal of University of Science and Technology Beijing, 2002, 24(3): 248–253. DOI: 10.13374/j.issn1001-053x.2002.03.047. [5] 龚敏, 王德胜, 黎剑华. 全息干涉法在条形药包离面位移场研究中的应用 [J]. 爆炸与冲击, 2005, 25(3): 227–231. DOI: 10.11883/1001-1455(2005)03-0227-05.GONG M, WANG D S, LI J H. Application of holographic interferometry to study vertical displacement field in linear charge [J]. Explosion and Shock Waves, 2005, 25(3): 227–231. DOI: 10.11883/1001-1455(2005)03-0227-05. [6] 郭洋, 李清, 杨仁树, 等. 三维模型柱状药包爆生裂纹扩展规律研究 [J]. 振动与冲击, 2020, 39(10): 133–140, 184. DOI: 10.13465/j.cnki.jvs.2020.10.018.GUO Y, LI Q, YANG R S, et al. Study on crack propagation law of cylindrical charges in three-dimensional models [J]. Journal of Vibration and Shock, 2020, 39(10): 133–140, 184. DOI: 10.13465/j.cnki.jvs.2020.10.018. [7] 杨仁树, 郭洋, 李清, 等. 中间起爆柱状药包爆炸应力应变场演化规律 [J]. 煤炭学报, 2019, 44(11): 3423–3431. DOI: 10.13225/j.cnki.jccs.2018.1673.YANG R S, GUO Y, LI Q, et al. Evolution law on explosive stress and strain field of column charges at middle detonation position [J]. Journal of China Coal Society, 2019, 44(11): 3423–3431. DOI: 10.13225/j.cnki.jccs.2018.1673. [8] 陈士海, 胡帅伟, 初少凤. 微差时间及柱状装药特征对爆破振动效应影响研究 [J]. 岩石力学与工程学报, 2017, 36(S2): 3974–3983. DOI: 10.13722/j.cnki.jrme.2017.1163.CHEN S H, HU S W, CHU S F. Study on the blasting vibration effect influenced by millisecond time and cylindrical charging characteristics [J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(S2): 3974–3983. DOI: 10.13722/j.cnki.jrme.2017.1163. [9] 李启月, 李夕兵, 范作鹏, 等. 深孔爆破一次成井技术与应用实例分析 [J]. 岩石力学与工程学报, 2013, 32(4): 664–670. DOI: 10.3969/j.issn.1000-6915.2013.04.003.LI Q Y, LI X B, FAN Z P, et al. One time deep hole raise blasting technology and case study [J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(4): 664–670. DOI: 10.3969/j.issn.1000-6915.2013.04.003. [10] 杨国梁, 杨仁树, 姜琳琳. 轴向间隔装药爆破沿炮孔的压力分布 [J]. 爆炸与冲击, 2012, 32(6): 653–657. DOI: 10.11883/1001-1455(2012)06-0653-05.YANG G L, YANG R S, JIANG L L. Pressure distribution along borehole with axial air-deck charge blasting [J]. Explosion and Shock Waves, 2012, 32(6): 653–657. DOI: 10.11883/1001-1455(2012)06-0653-05. [11] 胡涛, 李祥龙, 关思, 等. 分段装药结构爆破效果的数值模拟研究 [J]. 中国工程科学, 2014, 16(11): 36–41. DOI: 10.3969/j.issn.1009-1742.2014.11.005.HU T, LI X L, GUAN S, et al. Numerical simulation of structure performance under blasting with piecewise charge [J]. Strategic Study of CAE, 2014, 16(11): 36–41. DOI: 10.3969/j.issn.1009-1742.2014.11.005. [12] 杨仁树, 鲍舟琦, 王雁冰, 等. 数码电子雷管在高瓦斯矿井岩巷爆破掘进中的应用 [J]. 金属矿山, 2022(7): 27–34. DOI: 10.19614/j.cnki.jsks.202207004.YANG R S, BAO Z Q, WANG Y B, et al. Application of digital electronic detonator in blasting excavation of rock roadway in high gas mine [J]. Metal Mine, 2022(7): 27–34. DOI: 10.19614/j.cnki.jsks.202207004. [13] 杨仁树, 王雁冰, 张召冉, 等. 井巷工程掏槽爆破新技术及应用 [J]. 中国科学基金, 2022, 36(1): 120–127. DOI: 10.16262/j.cnki.1000-8217.2022.01.028.YANG R S, WANG Y B, ZHANG Z R, et al. New technology and application of cutting blasting in shaft and roadway engineering [J]. Bulletin of National Natural Science Foundation of China, 2022, 36(1): 120–127. DOI: 10.16262/j.cnki.1000-8217.2022.01.028. [14] GONZÁLES G L G, GONZÁLEZ J A O, DE CASTRO J T P, et al. Measuring elastoplastic strain loops in the near crack-tip region using a stereo microscope DIC system [J]. International Journal of Fatigue, 2020, 133: 105427. DOI: 10.1016/j.ijfatigue.2019.105427. [15] 左进京, 杨仁树, 汪文良, 等. 条形药包爆炸全场应变以及裂纹动态断裂特性研究 [J]. 工程科学学报, 2022, 44(8): 1306–1314. DOI: 10.13374/j.issn2095-9389.2021.10.19.003.ZUO J J, YANG R S, WANG W L, et al. Explosive full-field strain and crack dynamic fracture characteristics of a linear shaped charge [J]. Chinese Journal of Engineering, 2022, 44(8): 1306–1314. DOI: 10.13374/j.issn2095-9389.2021.10.19.003. [16] JU Y, ZHENG J T, EPSTEIN M, et al. 3D numerical reconstruction of well-connected porous structure of rock using fractal algorithms [J]. Computer Methods in Applied Mechanics and Engineering, 2014, 279: 212–226. DOI: 10.1016/j.cma.2014.06.035. [17] CHEN H D, CHENG Y P, ZHOU H X, et al. Damage and permeability development in coal during unloading [J]. Rock Mechanics and Rock Engineering, 2013, 46(6): 1377–1390. DOI: 10.1007/s00603-013-0370-2. 期刊类型引用(2)
1. 陈程,李全明. 地下矿山巷道掘进凿岩爆破技术参数优化研究进展与展望. 中国矿业. 2023(10): 103-110 . 百度学术
2. 费鸿禄,纪海楠,山杰. 露天台阶水介质间隔装药结构优选及对比试验研究. 黄金科学技术. 2023(06): 930-943 . 百度学术
其他类型引用(1)
-