巷道深孔掘进不同药径比爆破动力研究.pdf
第34卷 第3期 2017年9月 爆 破 BLASTING Vol. 34 No. 3 Sep. 2017 doi10. 3963/ j. issn. 1001 -487X. 2017. 03. 012 巷道深孔掘进不同药径比爆破动力研究* 柴修伟 1, 阎要锋1, 倪小山2, 习本军2, 徐 亮 2, 陈国祥3, 许 奎 1 (1.武汉工程大学资源与土木工程学院,武汉430074;2.湖北兴发化工集团股份有限公司, 兴山443711; 3.河南工程学院安全工程学院, 郑州451191) 摘 要 近年来, 随着矿山机械化程度的不断提高, 大型先进凿岩台车的普遍推广应用, 某些矿山在深孔掘 进爆破过程中出现了炮孔药径比不匹配的问题。为了确定合理的炮孔药径比, 运用有限元数值模拟软件 ANSYS/ LS-DYNA构建巷道深孔爆破中装药结构的计算模型,4种炮孔药径比分别取1∶ 1、1∶ 1. 2、1∶ 1. 225、 1∶ 1. 25, 对典型单元的应力云图和压力时程曲线进行对比分析。计算结果表明 当炸药与孔壁之间有空气介 质存在时, 应力波压力峰值会快速衰减, 能够减少炮孔壁岩石的过粉碎或造成“再生岩石”现象。此外当炮 孔药径比为1∶ 1时, 单元压力峰值随着远离炮孔装药轴线衰减最快, 而炮孔药径比1∶ 1. 2时, 压力峰值则衰 减最慢。对于钻头直径φ 45mm的凿岩台车, 合适的炸药药卷直径为φ 38 ~40mm, 可依此指导矿山企业合 理定购爆破器材。 关键词 地下矿山;巷道掘进;深孔掏槽;药径比;动力分析 中图分类号 TU235. 1 文献标识码 A 文章编号 1001 -487X(2017)03 -0068 -06 Study on Explosion Dynamics of Different Ratio of Charge to Hole Diameter in Deep Hole Tunnel Blasting CHAI Xiu -wei1,YAN Yao-feng1,NI Xiao-shan2,XI Ben-jun2,XU Liang2,CHEN Guo-xiang3,XU Kui1 (1. School of Environment and Civil Engineering,Wuhan Institute of Technology, Wuhan 430074,China;2. Hubei Xingfa Chemicals Group Co Ltd,Xingshan 443711,China; 3. Department of Safety Engineering,Henan Institute of Engineering,Zhengzhou 451191,China) Abstract In recent years,following the continuously improvement of the mine mechanization,there exist the problem of bad ratio of charge to hole diameter during the deep hole tunnel blasting in some mine with the wide- spread application of the advanced drill jumbo. To decide the reasonable ratio of charge to hole diameter,the software ANSYS/ LS-DYNA was used to establish the numeration model of charge structure of deep hole tunnel blasting. There are 4 types ratios with 1∶ 1、1∶ 1. 2、1∶ 1. 225、1∶ 1. 25,and the stress nephogram and the pressure time interval curve of typical unit were compared and analyzed. The results show that,when there is air material between charge and hole wall,the peak of stress pressure will attenuate rapidly,which could reduce the over broken of hole wall or“renew rock”. Moreover,when the ratio is 1∶ 1,the peak of unit pressure reduces fastest far away to the charge,but when the ratio is 1∶ 1. 2,the peak of unit pressure reduces slowest relatively. For the bit diameter of drill jumbo is φ 45 mm, the befitting charge diameter is φ 38 ~40 mm,which is of benefit to order the blasting materials for mine company. Key words underground mine;tunnel advance;deep hole cutting;ratio of charge to hole diameter;dynamical analysis 万方数据 收稿日期2017 -08 -01 作者简介柴修伟(1980 -) , 男, 博士、 副教授, 从事爆破工程及安全 监测技术研究, (E-mail)342650018@ qq. com。 通讯作者徐 亮(1976 -) , 男, 本科、 工程师, 从事地下矿山生产系 统管理及矿山机械化研究, (E-mail)xuliang62333@ 163. com。 基金项目湖北省教育厅科学技术研究项目(D20171506) ; 湖北省安 全生产专项资金项目( 鄂安监发(2017)35号) ; 武汉工程 大学科学研究基金项目(13125042) 目前, 地下矿山巷道工程主要还是采用炮孔爆 破法进行掘进, 为了提高巷道掘进的效率, 国内外很 多学者从优化巷道掘进的掏槽方式和爆破参数方面 对巷道掘进爆破进行了系统、 详细的理论和实验研 究, 如郑祥滨[ 1]、 李萍对螺旋空孔直眼掏糟成腔过 程进行了数值模拟研究[ 2], 单仁亮[3]、 郭东明对岩 巷掘进直眼掏槽爆破技术进行了实验研究[ 4], 凌同 华等对菱形楔形混合掏槽爆破方式进行数值模拟和 应用研究[ 5]。除此之外, 炮孔中炸药和炮孔的空间 耦合状态, 对炸药的爆炸性能、 炮孔壁初始压力峰值 大小、 炸药能量利用率以及破岩效果有直接的影响。 钟芮英对炮孔和装药直径的匹配关系进行了分析说 明[ 6]。Yilmaz O[7], Jae Kwang Ahn从应力波的角度 研究了爆破作用下岩体的破坏[ 8]。杨浩[9]、罗 强[ 10]、 宗琦等对巷道掘进爆破采用耦合装药和不耦 合装药对岩体中的动态应力场、 爆破效果的影响进 行了系统的研究[ 11]。 在矿山巷道掘进爆破中多采用水平钻孔, 装填 卷装炸药, 为了方便炮孔装药, 一般装药直径要比钻 孔直径小, 合理的炮孔药径比在方便装药和确保充 足破岩能量的前提下, 能够降低冲击波作用在炮孔 壁上的压力, 增加爆生气体的作用时间, 改进爆破效 果。近年来, 随着矿山机械化程度的不断提高, 大型 先进凿岩台车的普遍推广应用, 某些矿山掘进爆破 过程中出现了炮孔药径比不匹配的问题, 凿岩台车 的钻头直径采用φ 45 mm, 而炸药药卷还是多采用 φ 32mm标准药卷, 炮孔药径比达到1∶ 1. 5左右, 而 合理的炮孔药径比一般不宜大于1∶ 1. 25。针对目 前地下矿山采用凿岩台车进行钻爆作业, 为了确定 合理的炮孔药径比, 根据钻孔设备来选定合适的炸 药药卷规格, 运用有限元数值模拟软件ANSYS/ LS - DYNA构建巷道深孔爆破中装药结构的计算模 型, 对炮孔不同药径比1∶ 1( 装药机或混装车装药) 、 1∶ 1. 2、1∶ 1. 225、1∶ 1. 25等药径比进行数值模拟, 对 比分析不同药径比的岩石破碎效果。 1 理论分析 根据岩石爆破破碎基本理论观点, 岩石的破坏 是应力波和爆生气体共同作用的结果, 爆炸冲击波 传递到岩石中, 以应力波方式使岩石产生大量小裂 隙, 随后爆生气体膨胀作用使这些裂隙继续扩展、 贯 穿。有实验研究表明[ 12], 爆炸冲击波能量占炸药总 能量的比例较小, 约为10% ~ 20%; 而爆生气体膨 胀作用能量占炸药总能量比重较大, 超过一半以上。 耦合装药时, 爆炸冲击波直接作用在岩石上, 将爆炸 能量直接传给岩石; 当炮孔中炸药与孔壁之间存在 空气介质时, 爆炸冲击波要先压缩空气介质后, 再以 应力波方式作用于岩石。压缩空气需要消耗能量, 虽然直接作用在岩石上的能量受到一定程度的影 响, 但由于空气介质的存在, 产生冲击波的反射和折 射, 也间接提高了爆生气体的作用能量。 2 数值模拟计算 2. 1 本构模型建立 结合某磷矿地下巷道深孔掘进爆破工程, 采用 凿岩台车进行钻孔,孔深一般为3. 5m,孔径为 45 mm, 为了简化模型和便于后续计算, 建立炸药、 空气、 岩石三种材料的准三维计算模型, 由于岩体纵 向尺寸较大不考虑端部效应的影响, 将此类问题归 结为平面应变问题进行处理, 建立1/2模型, 采取对 称面施加对称约束和设置无反射边界,采用 cm/ g/ μs进行建模, 模型尺寸长为(Y方向)420 cm, 宽(X方向) 为200 cm,Z方向1个单位长度。 利用相同的建模思路, 分别建立炮孔药径比为 1∶ 1、1∶ 1. 2、1∶ 1. 225、1∶ 1. 25时4个计算模型, 其一 如图1所示。 图1 计算模型和典型单元选择的示意图 Fig. 1 Schematic diagram of model and units 96第34卷 第3期 柴修伟, 阎要锋, 倪小山, 等 巷道深孔掘进不同药径比爆破动力研究 万方数据 2. 2 材料本构方程及参数 ANSYS/ LS-DYNA软件是大型的有限元显式求 解程序, 其在岩土爆破中已经得到广泛应用, 本模型 包括炸药、 岩石、 空气三种材料, 采用流固耦合算法, 以Lagrangian算法为主, 兼有ALE和Euler算法, 通 过*CONSTRAINED-LAGRANGE-IN-SOLID进行耦 合, 避免计算过程中单元严重畸变引起的计算困难。 炸药采用MAT HIGH EXPLOSIVE_ BURN模型 以及JWL状态方程, 其表达式为 P = A 1 - ω R1 J V e-R1V+ B 1 - ω R2 J V e-R2V+ ωE 0 V (1) 式中P为压力;V为体积变化;R1、R2、ω、B和 A均为材料常数;E0为初始比内能。 岩石单元选用的LS-DYNA材料模型库内的各 向异性弹塑性模型。 炸药、 岩石、 空气三种材料的计算参数如表1 ~ 表3所示。 表1 炸药参数 Table 1 Explosive parameters 密度 ρ/ (gcm -3) 爆速D/ (cm μs -1) PCJ/ MPaK/ MPaG/ MPa 1. 00. 363. 24E -0200 表2 岩石材料参数 Table 2 Rock material parameters 密度 ρ/ (gcm -3) 弹性模量 E/ GPa 泊松比 μ 应变率 参数SRC 硬化参数 BETA 202. 850. 350. 0421. 0 表3 空气参数 Table 3 Air material parameters 密度 ρ/ (gcm -3) 截断 压力pc 动态黏度 系数MU 相对 体积 弹性模量 E/ GPa 泊松比 μ 1. 2929E -030. 00. 00. 00. 00. 0 2. 3 计算过程应力云图 将炮孔药径比为1∶ 1、1∶ 1. 2、1∶ 1. 225、1∶ 1. 25 时4个计算模型输入到LS - DYNA后处理程序中 分别计算, 得到4种不同炮孔药径比的计算过程。 图2所示为孔药径比为1∶ 1. 2时的计算应力云图。 从炮孔药径比为1∶ 1. 2时的计算应力云图中可 以看出, 炸药在孔底开始起爆, 爆炸冲击波在炮孔周 围岩体的有效应力场呈现椭圆状, 随着应力波在岩 体中继续传播, 应力场会以炸药起爆点为中心沿着 炮孔朝炮孔自由面方向传播, 并在自由界面处形成 反射拉伸波, 并与正应力场形成叠加。 图2 炮孔药径比为1∶ 1. 2时的计算应力云图 Fig. 2 The stress nephogram when ratio of charge to hole diameter is 1∶ 1. 2 3 结果分析 3. 1 单元压力峰值对比分析 为了对比分析不同炮孔药径比条件下, 爆炸应 力波在炮孔周围岩体的分布和衰减规律, 对每个计 算模型选取9个典型单元, 沿着炮孔装药竖向方向 选取编号21773、84353、64853、45158、36968的单 元, 垂直炮孔方向选取编号64853、64803、64753、 64703、64668的单元。不同炮孔药径比时典型单元 的压力峰值如表4所列。 从表4中可以看出, 在典型单元位置处, 随着炮 孔药径比的减小, 初始应力波峰值压力是逐渐降低。 从应力分布规律来看, 在沿炮孔走向的竖向单元中, 64853单元的应力波压力峰值最大,分别为 3917 MPa、2980 MPa、2867 MPa、2734 MPa, 而36968 单元的压力峰值最小,为760 MPa、496 MPa、 462 MPa、417 MPa, 符合柱状药包爆炸应力波分布 规律; 在沿炮孔垂向的横向单元中, 随着单元远离炮 孔轴线,单元的应力波压力峰值迅速衰减,从 3917 MPa、2980 MPa、2867 MPa、2734 MPa迅速降低 到542 MPa、624 MPa、449 MPa、438 MPa。将竖向单 元和横向单元随着炮孔药径比变化时的变化趋势拟 合成变化曲线, 如图3所示。 从图3中变化曲线可以看出, 在沿炮孔走向方 07爆 破 2017年9月 万方数据 向, 不同药径比条件下, 单元的压力峰值变化趋势是 一致的, 药径比1∶ 1时各单元的压力比其他三种药 径比条件下的压力峰值高很多, 说明在炮孔和孔壁 之间有空气存在时, 爆炸冲击波经空气介质后再作 用到孔壁岩石, 应力波压力峰值会快速衰减, 能够减 少炮孔壁岩石的过粉碎或造成“再生岩石”现象; 从 在垂直炮孔轴线的横向单元变化曲线看, 炮孔药径 比为1∶ 1时, 单元压力峰值随着远离炮孔装药轴线 衰减最快, 从3917 MPa降低至542 MPa; 炮孔药径 比为1 ∶ 1. 225(从2867 MPa降低至449MPa)和 1∶ 1. 25( 从2734 MPa降低至438 MPa)的衰减趋势 次之, 炮孔药径比1∶ 1. 2时,从2980 MPa降低至 624 MPa, 单元压力峰值则衰减最慢, 应力作用时间 最长, 在方便装药的前提下, 炮孔药径比1∶ 1. 2比其 他药径比更合适。 表4 不同的药径比下在的竖向、 横向单元压力 Table 4 Vertical and horizontal pressure under different ratio of charge to hole diameter 竖向单元 横向单元 单元编号炮孔药径比单元压力/ MPa单元编号炮孔药径比单元压力/ MPa 21773 1∶ 12280 64853 1∶ 13917 1∶ 1. 222491∶ 1. 22980 1∶ 1. 22521491∶ 1. 2252867 1∶ 1. 2520301∶ 1. 252734 84353 1∶ 12498 64803 1∶ 11485 1∶ 1. 220761∶ 1. 22979 1∶ 1. 22520201∶ 1. 2251233 1∶ 1. 2519731∶ 1. 251218 64853 1∶ 13917 64753 1∶ 1920 1∶ 1. 229801∶ 1. 21237 1∶ 1. 22528671∶ 1. 225750 1∶ 1. 2527341∶ 1. 25737 45158 1∶ 13854 64703 1∶ 1720 1∶ 1. 228541∶ 1. 2772 1∶ 1. 22527441∶ 1. 225612 1∶ 1. 2526571∶ 1. 25598 36968 1∶ 1760 64668 1∶ 1542 1∶ 1. 24961∶ 1. 2624 1∶ 1. 2254621∶ 1. 225449 1∶ 1. 254171∶ 1. 25438 图3 不同的药径比下的竖向、 横向典型单元压力变化曲线 Fig. 3 Curve of pressure change of typical vertical and horizontal units under different ratio of charge to hole diameter 3. 2 单元压力时程曲线对比分析 从上面可以看出, 炮孔药径比的变化会会影响 岩石中的应力波压力峰值。此外, 为了对比分析炮 孔药径比对爆炸冲击波在岩石中的传播和衰减规律 17第34卷 第3期 柴修伟, 阎要锋, 倪小山, 等 巷道深孔掘进不同药径比爆破动力研究 万方数据 的影响, 将选择典型单元64853和64668, 将这两个 单元在不同炮孔药径比条件下的单元压力时程曲线 进行对比, 各单元的压力时程曲线如图4、5所示。 图4 不同的药径比下的典型单元(64853) 压力曲线 Fig. 4 Pressure curve of horizontal special units(64853)in different ratio of charge and diameter 图5 不同的药径比下的典型单元(64668) 压力曲线 Fig. 5 Pressure curve of horizontal special units(64668)in different ratio of charge and diameter 从图4中看出, 单元64853在不同炮孔药径比 下, 初始应力波都是快速达到压力峰值, 并很快衰减 下来, 由于距离炮孔装药很近, 岩石在初始应力波作 用下已受压破碎;此外,由于在药径比为1∶ 1. 2、 1∶ 1. 225、1∶ 1. 25时, 装药与炮孔之间存在空气介 质, 使得初始应力波衰减曲线出现更多的上下起伏 变化。 从图5中看出, 单元64668的初始应力波在快 27爆 破 2017年9月 万方数据 速达到压力峰值后, 衰减的相对较缓, 并出现了反射 波。虽然在不同炮孔药径比条件下, 初始应力波的 压力峰值变化较大, 但随后出现的反射波压力峰值 基本相近, 可能由于爆生气体的膨胀静压作用引 起的。 4 结论 目前, 在矿山安全生产“机械化换人, 自动化减 人” 的指导思想下, 大型先进的凿岩台车会逐步成 为矿山钻孔作业的主力设备, 为直径φ 45mm的凿 岩台车钻头所钻凿的炮孔确定合适的炸药药卷, 已 成为矿山爆破作业的迫切任务。 从对所选的炮孔不同药径比1∶ 1(装药机或混 装车装药) 、1∶ 1. 2、1∶ 1. 225、1∶ 1. 25的模拟结果来 看, 当炮孔与孔壁之间有空气介质存在时, 爆炸冲击 波经空气介质后再作用到孔壁岩石, 应力波压力峰 值会快速衰减, 能够减少炮孔壁岩石的过粉碎或造 成“ 再生岩石”现象。另外从岩石中应力波的传播 和衰减规律来看, 炮孔药径比为1∶ 1时, 单元压力峰 值随着远离炮孔装药轴线衰减最快, 而炮孔药径比 1∶ 1. 2时, 单元压力峰值则衰减最慢, 应力作用时间 最长, 在方便装药的前提下, 炮孔药径比1∶ 1. 2比其 他药径比更合适。对于钻头直径φ 45mm的凿岩台 车, 在成孔的条件下, 根据炮孔药径比1∶ 1. 2的原 则, 合适的炸药药卷直径为φ 38 ~40 mm, 指导矿山 企业合理定购爆破器材。 参考文献(References) [1] 郑祥滨, 璩世杰, 范利华, 等.单螺旋空孔直眼掏糟成 腔过程数值模拟研究[J].岩土力学,2008,29(9) 2589-2594. [1] ZHENG Xiang-bin,QIU Shi-jie,FAN Li-hua,et al. Nu- merical simulation of cavity formation process by single helix parallel hole cut blasting[J]. Rock and Soil Me- chanics,2008,29(9) 2589-2594.(in Chinese) [2] 李 萍.渐进式大直径空孔螺旋掏槽爆破参数的研究 [D].武汉 武汉理工大学硕士学位论文,2010. [2] LI Ping. Study on blasting parameter of progressive large hole spiral cutting[D]. WuhanWuhan University of Technology,2010.(in Chinese) [3] 单仁亮, 黄宝龙, 高文蛟, 等.岩巷掘进准直眼掏槽爆 破新技术应用实例分析[J].岩石力学与工程学报, 2011,30(2) 224-232. [3] SHAN Ren-liang,HUANU Bao-long,GAO Wen-jiao,et al. Case studies of new technology application of quasi- parallel cut blasting in rock road way drivagc[J]. Chinese Journal of RockMechanicsandEngineering,2011, 30(2) 224-232.(in Chinese) [4] 郭东明, 王 豪, 李孝林, 等.益新矿中心大空孔掏槽 爆破扩腔过程数值分析[J].爆破,2017,34(2) 9-14. [4] GUO Dong-ming,WANG Hao,LI Xiao-lin,et al. Numeri- cal analysis on expansion chamber process of large hollow hole cut blasting in Yixin coal mine[J]. Blasting,2016, 34(2) 9-14.(in Chinese) [5] 凌同华, 曹 峰, 李 洁.菱形楔形混合掏槽爆破的数 值模拟研究与应用[J].爆破,2016,33(4) 9-14,33. [5] LING Tong-hua,CAO Feng,LI Jie. Numerical simulation research and application of diamond and wedge-shaped cutting blasting[J]. Blasting,2016,33(4) 9-14,33.(in Chinese) [6] 钟芮英, 李公照.炮孔与药卷直径之最佳匹配[J].世 界采矿快报,1994(3) 10-11 [6] ZHONG Rui-ying,LI Gong-zhao. The best match between the hole diameter and the diameter of the cartridge[J]. The World Mining Letters,1994(3) 10-11.(in Chinese) [7] YILMAZ O,UNLU T. Three dimensional numerical rock damage analysis under blasting load[J]. Tunn Undergr Space Technol,2013,38266-278. [8] JAE Kwang Ahn,DUHEE Park. Prediction of near-field wave attenuation due to a spherical blast source[J]. Rock Mech Rock Eng,20171274-1276. [9] 杨 浩.巷道掘进爆破耦合装药对爆破效果的影响研 究[D].淮南 安徽理工大学硕士学位论文,2016. [9] YANG Hao. Tunneling blasting charge coupled to study the influence of blasting effect[D]. HuainanAnhui Uni- versity of Science and Technology,2016.(in Chinese) [10] 罗 强.装药结构对岩体内爆炸应力场的影响研究 [D].淮南 安徽理工大学硕士学位论文,2006. [10] LUO Qiang. Study on the influence of blast stress field in rock mass due to modification of charge structure[D]. HuainanAnhui University of Science and Technology, 2006.(in Chinese) [11] 宗 琦, 马亚东, 汪海波.煤矿硬岩巷道掘进大直径 炮孔爆破试验研究[J].爆破,2017,34(1) 47-51. [11] ZONG Qi,MA Ya-dong,WANG Hai-bo. Large diameter hole blasting experimental in hard rock roadway of coalmine[J]. Blasting,2017,34(1) 47-51.(in Chi- nese) [12] J R布里克曼.分离冲击波和气体膨胀作用的破碎机 理[C]∥第二届爆破破岩国际会议论文集.长江科学 研究院水利水电爆破咨询服务部, 译. 1990(12) 7- 15. 37第34卷 第3期 柴修伟, 阎要锋, 倪小山, 等 巷道深孔掘进不同药径比爆破动力研究 万方数据