深孔爆破技术在高地应力低透气性高瓦斯煤层增透防突中的适用性＊

赵宝友，王海东

（1.辽宁工程技术大学力学与工程学院，辽宁 阜新123000；2.中国煤炭科工集团有限公司沈阳研究院煤矿安全技术国家重点实验室，辽宁 沈阳110016）

我国煤层赋存条件差，地质构造复杂，煤层瓦斯含量较高且多以吸附状态存在，这使得我国高瓦斯矿井和煤与瓦斯突出矿井占大部分，相应的瓦斯灾害也比较严重。李润求等［1］的统计资料表明，我国在2 001～2 010的10年内就发生了2 246起煤矿瓦斯事故，累计造成13 155人遇难；在瓦斯事故中，煤与瓦斯突出是发生最多的较大事故，同时造成人员遇难也最多。据第2次全国煤田预测结果［2］，我国煤炭资源总量为5.57万亿吨，其中埋深小于600m 的预测煤炭资源量占全国煤炭预测资源总量的26.8%，埋深在600～1 000m 的约占20%，埋深超过1 000m 的约占53.2%。煤炭作为国民经济的基础产业，在未来一段时间内仍是我国的主体能源［1－3］。迄今我国开采深度千米以上的煤矿至少有17座［3］，随着开采水平向深部的延伸，地质构造更复杂，煤层瓦斯压力增大，不少原来中浅埋深非突出矿井，转化为突出矿井，冲击地压和煤与瓦斯突出的强度和频度随采深增加明显增大。高瓦斯矿井瓦斯灾害治理的根本途径是煤层瓦斯的抽放。但我国大部分煤矿的煤层为低透气性煤层［2］，即便是高渗透性的晋城等煤气田的渗透率也仅为（1.0～18.0）×10－3mm2，比国外相同埋深煤层的渗透率低2～3个数量级。这一特点也决定了对我国煤矿采用传统预抽瓦斯方法开展瓦斯安全治理难度很大。

目前，采煤生产实践中主要采用水力压裂［4］、深孔松动爆破［5］等技术来改善低透气性煤层的渗透特性。实验研究［6］表明，煤层高压注水致裂后，煤炭含水的饱和度增加，高压裂隙水和孔隙水堵塞了瓦斯的渗流通道，又因煤层孔隙压力升高，大量游离的瓦斯转变为吸附态，抑制了瓦斯的解析效果，不利于煤层瓦斯的抽采。深孔爆破是利用在煤层中引爆炸药，使炮孔周围产生大量径向裂隙，形成瓦斯流向抽采孔的渗流通道，同时降低煤体刚度，使局部煤层的应力得到释放，进而达到煤层防冲、防突的目的［5］。深孔爆破增透技术起初是从浅孔采掘爆破中借鉴过来的，但由于煤岩体所承受的应力状态及煤岩体内部爆生裂隙（包括原生裂隙）与采掘临空面的空间位置不同以及爆破最终目的不同，使得煤层深孔爆破增透机理与常规浅孔采掘破岩机理不同［7］：常规浅孔采掘爆破中，爆生气体主要起到抛掷破碎煤岩块的外部动力作用，而主要依靠爆炸波的冲击震裂作用克服煤岩体的强度来破岩；深孔爆破中，爆炸波主要用来生成后继高压爆生气体流入煤体的初始裂隙，煤体裂隙扩展主要依靠高压爆生气体和煤层原始瓦斯压力的准静态尖劈致裂作用，爆破致裂增透期间除克服煤体强度还要克服煤层地应力，尤其对于深埋高应力煤层的深孔爆破，主要是依靠高压爆生气体克服煤层的地应力来产生大量的裂隙，达到增透目的。尽管深孔爆破增透效果已被中浅埋深煤层井工开采的工程实践证实较好［5，7－8］，但随采深的增加，不可避免会遇到高应力的低透气性煤层，因此，有必要研究该爆破技术对高应力煤层增透的适用性及煤层地应力对爆破增透的影响。

1 数值计算模型

鉴于深孔爆破影响范围内煤层的受力为3向应力状态，炮孔装药长度至少几十米，封孔长度也在10m 左右，因此，可将其简化为二维平面应变问题进行研究。以往有关煤体爆破裂隙扩展的理论研究大都只考虑爆炸波的影响，却很少考虑爆生气体对煤体裂隙扩展的影响，本文中，利用FLAC3D 有限差分软件的动力求解模块，同时考虑爆炸波、爆生气体、煤体原始瓦斯压力和煤层地应力的作用，对深埋煤层深孔爆破增透进行数值模拟，通过对动力计算结果的分析及与模型实验、现场实验和前人研究成果的对比，揭示深孔爆破技术在深埋高应力低透气性高瓦斯煤层增透防突中的适用性，详细探究煤层地应力对其增透效果的影响，以期为深埋高应力低透气性煤层深孔爆破工艺的设计提供参考。

1.1 有限差分模型

选取厚为4m 的煤层作为研究对象，只考虑单一炮孔且将其布置在煤层中间，炮孔直径为90mm，炮孔两侧和底部煤（岩）体厚度均取100m，上覆岩层一直建至地表，模型平面法向方向长度为1 m；煤层竖向应力为上覆岩层的自重，上覆岩层平均密度为2.5t／m3。地应力平衡的静力计算期间模型四周和底部为法向位移固定边界条件。动力计算期间模型四周和底部为黏弹性自由边界场条件，以吸收从炮孔爆炸内源传至模型边界的外行波。模型的几何示意图和有限差分网格局部放大图如图1所示。

图1 煤层深孔爆破示意图Fig.1Schematic figure of deep－hole blasting in coal seam

1.2 煤体的本构模型

煤体中大量层理、节理、割理裂隙等不连续面（体）的存在，使得煤体的力学强度显著降低。爆破作用下煤体内爆生裂隙的生成要克服煤层地应力和煤体的抗拉强度，合理的煤体抗拉强度是防突抽采设计的重要参数之一。以往理论分析和数值模拟研究中大都采用Mohr－Coulomb（MC）准则来描述煤岩体的破坏行为。加之利用FLAC3D 内嵌的MC本构可以描述煤体受拉破坏的力学行为，为此，本文数值模拟计算中对煤体同样选取具有拉剪复合破坏准则的MC 模型。采用分段线性拟合HB 模型的方法［9－10］获得等效MC模型参数，如表1所示。表中S 为煤体的地质强度指标参数，σci和mi分别为完整煤块的单轴抗压强度和材料常数，mb、s、a 均为煤体的材料常数，σcc和σct分别为煤体的单轴抗压强度和单轴抗拉强度，ccm为煤体的粘结力，Ei和Ecm分别为煤块和煤体的杨氏模量，φ 为煤体的内摩擦角。

?

1.3 爆炸荷载

煤岩中炸药引爆后，在爆炸波的动态冲击震裂及爆生气体的准静态应力作用下，沿炮孔径向依次生成压缩粉碎区、裂隙区及受应力波扰动未产生裂隙的震动区。鉴于爆炸过程的瞬时性，很难精确测定爆生气体和爆炸波的压力峰值，加之后续煤岩体中生成裂隙时空分布规律的复杂性，更难精确测定裂隙内爆生气体的变化规律。为此，对煤岩体爆破问题进行适当简化，将爆生气体视为等熵绝热膨胀气体，发展了一些理论方法计算炮孔内爆生气体和炮孔壁爆炸波的峰值压力［11－13］。

柱状耦合装药情况下，爆生气体和爆炸波峰值压力的理论公式分别为：

柱状非耦合装药情况下（忽略药柱与炮孔之间的空气间隙），爆生气体和爆炸波峰值压力的理论公式分别为：

式中：pg和ps分别为炮孔内爆生气体和炮孔壁处爆炸波的峰值压力；ρ0和ρc分别为炸药和煤岩体的密度；De和cp分别为炸药的爆速和煤岩体的压缩波波速；r1和r2分别为药柱和炮孔的半径；l1和l2分别为药柱和炮孔的长度；n 为增大因数，取值范围为8～11。

利用有限元软件ABAQUS的Explicit显示动力模块［14］可以模拟炸药的爆炸过程，即可获取炮孔壁的爆炸波时程曲线。已有学者［8，15－17］通过AUTODYN、ABAQUS等有限元软件获得了炸药爆炸的时程曲线。ABAQUS软件中基于Jones－Wilkens－Lee（JWL）状态方程来模拟炸药的爆炸过程，炸药引爆后爆源内任一点的压力p 可表达为：

式中：peos和F 分别为炸药的状态方程和燃烧分数；A、B、R1、R2和ω 为与炸药相关的材料常数；ρ0和ρ分别为炸药和爆生产物的密度；Em为炸药的质量内能；Bs为控制炸药燃烧的速度常数，可取2.5；le为网格单元的特征长度，td为炸药点燃后生成爆生气体的时刻。

选取煤矿许用型水胶（或乳化）炸药，装药形式为柱状耦合装药，采用上述理论公式计算得到爆生气体和爆炸波的压力峰值，利用有限元软件ABAQUS模拟炸药的爆炸过程，获得炮孔壁的爆炸波压力时程曲线。炸药的特性参数［14］及计算结果分别为：De＝2.7km／s，A＝215GPa，B＝0.19GPa，R1＝4.2，R2＝0.9，ω＝0.15，ρ0＝1t／m3，Em＝3.2kJ／g，pg＝0.92GPa，ps＝1.24GPa，pg，c＝0.55GPa。获 得 的炮孔壁爆炸波归一化的压力时程曲线如图2所示。

R.H.Nilson［18］、卢文波等［19］和李宁等［20］认为，爆生气体在裂纹中的衰减规律近似符合指数形式。本文中选取的爆生气体在裂纹中随时间变化的规律为：

式中：pg，c为粉碎区扩腔形成后爆生气体的压力，p′g为裂隙内t时刻的爆生气体压力；t0为爆生气体进入初始裂隙的起始时间；r2和rc分别为炮孔原始半径和生成粉碎区后的扩腔半径；α 为裂纹相对扩展速度，α＝v／cp，v 是裂纹扩展的绝对速度。

Griffith理论［18－19］指出，裂纹扩展的稳定扩展速度α＝v／cp＝0.38，本文中取α＝0.3。张奇［21］和宗琦［22］研究发现，岩石粉碎后生成的扩腔半径为炮孔半径的1.05～1.30倍，本文中取rc／r2＝1.09。那么，由上述公式（9）可计算得到粉碎区边界处裂隙气体的压力峰值pg，c＝0.55GPa。爆生气体的压力时程曲线如图3所示。

图2 归一化的爆炸波峰值压力时程曲线Fig.2 Normalized peak pressure curve of explosion wave

图3 爆生气体压力时程曲线Fig.3 Pressure curve of explosion－induced gas

2 结果分析

鉴于FLAC3D 的动力模块能很好地反映爆炸波在煤岩体内的传播规律［7］，本文中选取FLAC3D的动力模块对煤层深孔爆破增透进行数值模拟。

由于爆炸波作用下炮孔周围粉碎区及径向初始裂隙的生成仅需数百微秒至数毫秒，而后继爆生气体准静态作用下初始裂纹的二次扩展至少需几十毫秒至数秒，因此，可将爆炸波与爆生气体的作用分开来研究。动力计算时，施加1.0 ms爆炸波使炮孔周围煤体形成压缩屈服的粉碎区及初始的径向裂隙后，再通过自编的fish程序时时追踪和判断单元是否发生屈服，将爆生气体以等效节点力的形式施加在屈服单元节点上来近似模拟爆生气体的致裂作用。本文的动力数值计算中，以局部阻尼的形式近似考虑煤岩介质的阻尼效应，煤岩体的临界阻尼比均取5%。动力求解时间为1.0s；煤层的原始瓦斯压力取2.0 MPa。鉴于本文的数值计算是在基于牛顿第二定律的动力模块下进行，较以往以忽略爆破振动引起的惯性力、忽略爆炸波的传播规律及不考虑爆生气体的动态扩展裂纹效应，该数值模拟结果更符合实际情况。

煤体爆破产生的拉伸裂隙在很大程度上可以反映煤层深孔爆破增透的效果。下面对各工况下煤层爆破拉伸破坏区域的数值模拟结果进行分析，探究煤层地应力对煤层深孔爆破增透的影响及该爆破技术在深埋高应力煤层增透防突中的适用性。

2.1 等压地应力对煤层爆破增透的影响

对静水压力为10～40 MPa的煤层（对应的埋深分别为400、800、1 200和1 600m）深孔爆破进行数值模拟，爆破结束后煤层的拉伸破坏区如图4所示。由图4可知，地应力对煤层的深孔爆破效果影响显著；随着地应力的增大，煤层裂隙区域逐渐减小。

图5中给出了煤层单孔爆破影响半径随地应力的变化规律，从图中可以直观地看出，地应力为10～40 MPa时，裂隙区半径随地应力的增大而近线性地减小，线性拟合曲线的相关因数为0.96。图6给出了实验室模型实验得到的不同地应力下距爆破孔75倍炮孔半径处煤体波速与煤层地应力的关系曲线［7］。从图中可以看出，随煤层地应力的增大，距炮孔一定距离位置处煤体的波速随之近线性地增大，即随地应力的增大，炮孔周围裂隙密度和裂隙区范围随之减小。可见，地应力对煤层深孔爆破裂隙区影响的数值模拟结果与实验室模型实验结果基本一致。

另外，数值模拟结果也与一些地应力相当的工程实践的测试结果相近，进一步说明本文中数值计算结果的可信性，如：埋深450m 的鹤壁二矿36061回采工作面，深孔爆破工业实验单孔增透的有效影响半径为4m 左右［23］；埋深780m 的五龙矿3321回采工作面，煤层3向主应力分别为15、20和30MPa，深孔爆破工业实验单孔增透的有效影响半径为3.5 m 左右［7］；埋深约为800 m 淮南矿业集团潘三矿1741（3）回采工作面深孔爆破增透的有效影响半径为3～4m［24］；埋深1 100m 左右的平煤十二矿己15－31010机、风巷穿层控制卸压爆破的有效影响半径为2.5～3.5m［25］。

图4 不同地应力情况下裂隙区分布图（静水压力状态）Fig.4 Crack zone under different hydrostatic geo－stress conditions

图5 裂隙扩展半径随煤层地应力的变化Fig.5 Variation of crack radius with geo－stress endured by coal seam

图6 压缩波速度随煤层地应力的变化［7］Fig.6 Velocity of compressional wave varied with geo－stress endured by coal seam［7］

2.2 非等压地应力对煤层爆破增透的影响

图7 中给出了30 MPa定值竖向应力、不同水平侧压因数（λx、λy）下煤层单孔爆破裂隙区分布情况。从图中可知，随水平侧压因数的增大，裂隙区总体积和裂隙区的竖向半径随之减小，水平向半径有增大的趋势。这是因为，在竖向地应力保持不变的情况下，随水平地应力的增大，由于煤层水平应力的逐渐增大，导致垂直水平向的裂隙，即竖向裂隙的发展逐渐得到抑制，使得有限的爆能将向有利于裂隙发展的水平向传播，进而使水平向的裂隙范围增大。另外，煤层地应力对煤体爆生裂隙的影响可以通过断裂力学进行解释，煤层地应力的增大使得应力场作用下煤体裂纹尖端的应力强度因子K 降低，不利于煤体裂隙的生成［7，12］：式中：D 为裂隙尖端的损伤值；L（t）为t时刻裂隙的扩展长度；σinf为垂直裂隙面的远场地应力；x 为裂隙扩展方向的坐标。

式（10）表明，在距炮孔中心x≤L＋r2处（L 为裂隙最终的扩展长度），若σinf增大，那么K 随之减小，使得L 减小，如图7所示。因此，受地应力影响显著的深埋煤层深孔爆破工艺的设计，应结合实际煤层分布情况及其地应力的分布规律，适当调整爆破孔与控制孔的位置及间距，使煤层深孔爆破能够获得最佳的爆破增透效果。炮孔位置及孔间距的布置建议方案如下：爆破孔与控制孔的中心连线应尽可能地与地应力的最大主应力方向平行；平行于地应力最小主力方向的孔间距s1与平行于地应力最大主应力方向的孔间距s2的比值lk应等于2主应力量值的比值，即lk＝s1／s2＝σmin／σmax，式中σmin和σmax分别为煤层炮孔横断面平面内的最小和最大主应力。

图7 不同地应力侧压力因数下裂隙区的分布Fig.7 Crack zones at different lateral coefficients of initial geo－stress

3 结 论

通过不同地应力下低透气性煤层深孔爆破增透的数值模拟以及与实验室模型实验和相关现场实验的对比，可得如下结论。

（1）煤层地应力对煤层爆破增透具有显著的抑制作用；深孔爆破的裂隙区半径随地应力的增大而近线性地减小；但深孔爆破技术应用于深埋高应力低透气性煤层同样可获得比较理想的增透防突效果。

（2）煤层地应力的主应力方向在一定程度上影响着煤层深孔爆破裂隙的扩展方向：在竖向地应力保持不变的情况下，随水平地应力的逐渐增大，炮孔周围竖向裂隙的发展逐渐得到抑制，有限的爆能将向有利于裂隙发展的水平向传播，使水平向的裂隙范围增大。

（3）在煤层深孔爆破增透的实际工程中，需结合煤层地应力状况来布置炮孔的空间位置，爆破孔与控制孔的中心连线应尽可能地与地应力的最大主应力方向平行。

［1］ 李润求，施式亮，念其锋，等.近10年我国煤矿瓦斯灾害事故规律研究［J］.中国安全科学学报，2011，21（9）：143－151.Li Run－qiu，Shi Shi－liang，Nian Qi－feng，et al.Research on coalmine gas accident rules in China in recent decade［J］.China Safety Science Journal，2011，21（9）：143－151.

［2］ 袁亮，薛俊华.中国煤矿瓦斯治理理论与技术［C］∥2010中国煤矿瓦斯治理国际研讨会论文集.合肥：中国矿业大学出版社，2010：3－16.

［3］ 胡社荣，彭纪超，黄灿，等.千米以上深矿井开采研究现状与进展［J］.中国矿业，2011，20（7）：105－110.Hu She－rong，Peng Ji－chao，Huang Can，et al.An overview of current status and progress in coal mining of deep over a kilometer［J］.China Mining Magazine，2011，20（7）：105－110.

［4］ 林柏泉，李子文，翟成，等.高压脉动水力压裂卸压增透技术及应用［J］.采矿与安全工程学报，2011，28（3）：452－455.Lin Bo－quan，Li Zi－wen，Zhai Cheng，et al.Pressure relief and permeability increasing technology based on high pressure pulsating hydraulic fracturing and its application［J］.Journal of Mining ＆Safety Engineering，2011，28（3）：452－455.

［5］ 邢昭芳，阎永利，李会良.深孔控制卸压爆破防突机理和效果考察［J］.煤炭学报，1991，16（2）：1－9.Xing Zhao－fang，Yan Yong－li，Li Hui－liang.Mechanism and results of prevention of rockburst by controlled pressure relief deep holes［J］.Journal of China Coal Society，1991，16（2）：1－9.

［6］ 赵东，冯增朝，赵阳升.高压注水对煤体瓦斯解吸特性影响的试验研究［J］.岩石力学与工程学报，2011，30（3）：547－555.Zhao Dong，Feng Zeng－chao，Zhao Yang－sheng.Experimental study of effects of high pressure water injection on desorption characteristic of coal－bed methane（CBM）［J］.Chinese Journal of Rock Mechanics and Engineering，2011，30（3）：547－555.

［7］ 王海东.深部开采低渗透煤层预裂控制爆破增透机理研究［D］.哈尔滨：中国地震局工程力学研究所，2012：47－56.

［8］ 龚敏，王华，文斌.岩石深孔爆破对邻近煤层的动应力作用［J］.爆炸与冲击，2012，32（2）：196－202.Gong Min，Wang Hua，Wen Bin.Dynamic stress in adjacent coal seams induced by deep－hole blasting in rock［J］.Explosion and Shock Waves，2012，32（2）：196－202.

［9］ Hoek E.Pratical rock engineering［M／OL］.［2011－03－01］.http：／／www.rocscience.com／hoek／corner／Practical＿Rock＿Engineering.pdf.Chapter 11：1－28.

［10］ 杨圣奇，温森.不同直径煤样强度参数确定方法的探讨［J］.岩土工程学报，2010，32（6）：881－891.Yang Sheng－qi，Wen Sen.Method for strength parameters of coal samples with different diameters［J］.Chinese Journal of Geotechnical Engineering，2010，32（6）：881－891.

［11］ 郭进平，聂兴信.新编爆破工程实用技术大全［M］.北京：光明日报出版社，2002：18－22.

［12］ 吴亮，卢文波，宗琦.岩石中柱状装药爆炸能量分布［J］.岩土力学，2006，27（5）：735－739.Wu Liang，Lu Wen－bo，Zong Qi.Distribution of explosive energy consumed by column charge in rock［J］.Rock and Soil Mechanics，2006，27（5）：735－739.

［13］ 赵铮，陶钢，杜长星.爆轰产物JWL状态方程应用研究［J］.高压物理学报，2009，23（4）：277－284.Zhao Zheng，Tao Gang，Du Chang－xing.Application research on JWL equation of state of detonation products［J］.Chinese Journal of High Pressure Physics，2009，23（4）：277－284.

［14］ ABAQUS Inc.ABAQUS theory manual［M／OL］.［2011－03－11］.http：／／abaqusdoc.ucalgary.ca／v6.9／books／stm／default.htm.

［15］ Liu Li－qing，Katsabanis P D.A numerical study of the effects of accurate timing on rock fragmentation［J］.International Journal of Rock Mechanics and Mining Sciences，1997，34（5）：817－835.

［16］ Chen S G，Zhao J.A study of UDEC modeling for blast wave propagation in jointed rock masses［J］.International Journal of Rock Mechanics and Mining Sciences，1998，35（1）：93－99.

［17］ Wei X Y，Zhao Z Y，Gu J.Numerical simulations of rock mass damage induced by underground explosion［J］.International Journal of Rock Mechanics and Mining Sciences，2009，46（7）：1206－1213.

［18］ Nilson R H.An integral method for predicting hydraulic fracture propagation driven by gasses or liquids［J］.International Journal for Numerical and Analytical Methods in Geomechanics，1986，10（2）：191－211.

［19］ 卢文波，陶振宇.预裂爆破中炮孔压力变化历程的理论分析［J］.爆炸与冲击，1994，14（2）：140－147.Lu Wen－bo，Tao Zhen－yu.Theoretical analysis of the pressure variation in borehole for pre－splitting explosion［J］，Explosion and Shock Waves，1994，14（2）：140－147.

［20］ 李宁，陈莉静，张平.爆生气体驱动岩石裂纹动态扩展分析［J］.岩土工程学报，2006，28（4）：460－463.Li Ning，Chen Li－jing，Zhang Ping.Dynamic analysis for fracturing progress by detonation gas［J］.Chinese Journal of Geotechnical Engineering，2006，28（4）：460－463.

［21］ 张奇.岩石爆破的粉碎区及其空腔膨胀［J］.爆炸与冲击，1990，10（1）：68－95.Zhang Qi.Smash districts and expanding of cavities in rock blasting［J］.Explosion and Shock Waves，1990，10（1）：68－95.

［22］ 宗琦.岩石爆破的扩腔作用及能量消耗［J］.煤炭学报，1997，22（4）：392－396.Zong Qi.Function of cavity expansion and energy consumption in rock blasting［J］.Journal of China Coal Society，1997，22（4）：392－396.

［23］ 李鸿宽，吴继园，刘思远.深孔预裂爆破强化抽放低透气性特厚煤层瓦斯的实践［J］.煤矿安全，2003，34（2）：12－14.Li Hong－kuan，Wu Ji－yuan，Liu Si－yuan.Practice of using deep－hole presplitting blasting technology to drainage methane from extremely thick coal seam with low permeability［J］.Safety in Coal Mines，2003，34（2）：12－14.

［24］ 张兴华.利用深孔控制预裂爆破强化瓦斯抽放消除回采工作面突出危险性［J］.煤矿安全，2006，37（2）：22－24.Zhang Xing－hua.Strengthening methane drainage using deep－hole controlled pre－split blasting for eliminating risk of working face outburst［J］.Safety in Coal Mines，2006，37（2）：22－24.

［25］ 谢友友，张连军，林柏泉，等.穿层深孔控制爆破有效影响半径的确定［J］.煤矿安全，2008，39（11）：11－14.Xie You－you，Zhang Lian－jun，Lin Bo－quan，et al.Determination of effective influence radius for strata penetration controlled blasting［J］.Safety in Coal Mines，2008，39（11）：11－14.