基于震波CT探测的宽煤柱冲击地压防控技术

孙刘伟，鞠文君，潘俊锋，赵忠显，夏永学，冯美华，付玉凯

(1.煤炭科学研究总院 开采研究分院,北京 100013; 2.天地科技股份有限公司 开采设计事业部,北京 100013; 3.山东龙郓煤业有限公司,山东 菏泽 274700)

近年来，伴随着矿井开采深度与开采强度的增加，我国冲击地压灾害日益严重，冲击发生强度和频度显著增加，冲击地压矿井数量明显上升，已成为制约深部矿井安全开采的重大动力灾害之一[1-2]。深井区段煤柱是冲击地压易发、多发、难防治区域，其冲击危险性与周围应力场[3]分布密切相关。姜福兴等[4]建立了采场“载荷三带”覆岩结构模型，可用于区段煤柱冲击危险性评价;王书文[5]通过应力监测系统实测得到采空区侧向宽煤柱分阶段、分区域受力特征，指出临空煤柱帮内弹性区不断增长的垂直应力是主导冲击载荷源;潘立友等[6]提出软弱冲击层在上部力源层和下部稳定层夹持下产生应力集中和能量积聚，外界采掘活动诱发积聚的能量释放，是两软煤层区段小煤柱沿空巷道发生冲击地压的主要原因;杨伟利等[7]指出当煤柱冲击力大于阻抗力时，具有发生冲击的危险性;成云海等[8]研究了“应力场特征-减冲”关系，提出了巷道外错合理位态减冲设计方法;李学华等[9]基于应力控制和支护系统设计2个方面，提出了沿空掘巷的防冲机理。上述研究多以理论分析、数值模拟为主，具有重要的理论指导意义;而实测研究更贴近现场实际，更具有工程实用价值。

本文以某矿80 m区段宽煤柱冲击地压为例，分析了冲击地压致灾机制及类型，鉴于静载荷在冲击地压启动过程中的主导作用，采用震波CT原位探测技术反演冲击区域静载荷分布特征，评估其冲击危险性，并在此基础上提出针对性防治方法。

1 工程背景

1301工作面为一采区南翼西侧第2个工作面，埋深约865～970 m，平均煤厚6.8 m，综放开采，走向长度2 265.8 m，倾向长度223.4 m，工作面西侧前部与八里庄支四断层(H=5～110 m)相邻、后部隔80 m煤柱与1300工作面采空区相邻，东、南两侧为实体煤，北邻采区轨道大巷。直接顶为厚1.2 m的砂质泥岩，基本顶为厚12.9 m的中砂岩，底板为厚1.8 m的泥质砂岩。采用钻孔应力解除法对一采区进行地应力测量，地应力场以水平应力为主，最大水平主应力为垂直应力的1.40～2.20倍。1301工作面运输巷位于一采区背斜构造翼部，初始作为1300工作面泄水巷使用，先于1300工作面回采已施工掘进。

当1300工作面回采496 m时，80 m宽煤柱内发生2.1×105J大能量事件，经微震定位，震源距1301运输巷为13 m，距1300回采工作面煤壁为310 m。经现场排查，2 号联络巷以北125～280 m范围内的1301工作面运输巷具有不同程度的冲击破坏现象，破坏严重段10 m范围顶板冒落达2.5～3.0 m，两帮内挤，巷道断面由4.8 m×4.0 m缩减为1.1 m×0.8 m;靠近冒顶区30 m范围巷道顶板下沉0.5 m;巷道顶板隔爆棚被冲击脱落，局部两帮双抗网开裂，管路、电缆被冲到巷道中部，工作面及冲击显现位置如图1所示，巷道破坏状况如图2所示。

图1 工作面冲击显现位置Fig.1 Position of rockburst

图2 巷道冲击破坏情况Fig.2 Impact failure of roadway

2 宽煤柱冲击地压致灾机制分析

2.1 宽煤柱静载荷分布特征

采用FLAC3D数值软件，根据工作面开采条件及实验室测定煤岩体参数，基于摩尔-库伦本构模型，模拟1300工作面开挖前、后宽煤柱静载荷分布特征，并提取宽煤柱煤体垂直应力，如图3所示。

图3 宽煤柱垂直应力分布Fig.3 Vertical stress distribution in wide coal pillar

两巷开挖形成后，宽煤柱内静载荷变化不明显，此时高自重应力与强构造应力叠加作用，构成宽煤柱冲击地压发生的基础静载荷。1300工作面开采后，侧向支承压力影响范围约为120 m，宽煤柱震源位置煤体应力σ=34 MPa>2σc，已具备冲击地压发生的冲击临界支承压力条件[10]。当1301工作面回采后，宽煤柱受到两侧支承压力叠加影响，静载荷水平将明显提高，应力峰值向煤柱深处转移，冲击危险性进一步增大。采空区侧向支承压力为宽煤柱冲击地压发生的增量静载荷。

2.2 宽煤柱静载荷演化规律

当应力集中超过极限后煤岩体突然破坏释放能量产生微震事件，可通过其时空分布特征反映宽煤柱静载荷演化规律[11]。基于“最佳D值”设计准则，除1300回采工作面前方布置ARAMIS M/E微震系统探头外，1301工作面运输巷补充布置了拾震器，使得宽煤柱区域处在探头-拾震器包络的区域内，保证了定位精度，实时记录了本次冲击区域微震活动全过程。微震事件(≥102J)平面投影如图4所示。

图4 微震演化平面投影Fig.4 Plane of microseismic events

当1300工作面回采至冲击区域时(图4(a))，微震事件分布靠近宽煤柱内侧(采空区侧44 m范围)，3次方微震事件较多，主要为1300工作面采空区低位顶板回转下沉导致应力升高;伴随着工作面逐渐远离(图4(b))，微震事件分布向宽煤柱外侧(1301工作面侧)转移，表明采空区侧向支承压力向外侧转移;此后，该区域仍有微震活动(图4(c))，表明1300工作面采空区侧向高位覆岩仍在回转下沉加载煤岩体[5]，造成宽煤柱内煤岩体应力持续升高，弹性应变能不断积聚;直至宽煤柱内弹性变性能超过煤岩体动力破坏的最小能量时，发生冲击破坏(图4(d))。

2.3 宽煤柱冲击地压致灾机制

经鉴定，3 号煤层具有弱冲击倾向性，顶板岩层具有强冲击倾向性，已具备发生冲击地压的内在条件;高自重应力、强构造应力提供基础静载荷，采空区后方持续增长的侧向支承压力提供了增量静载荷，当两者叠加导致宽煤柱垂直应力超过冲击临界支承压力时，为冲击地压的发生提供了力源条件。并且，冲击震源距1300工作面煤壁310 m，已不再受到1300工作面采动影响，因此该冲击地压属于静载荷型，防治应以调控静载荷为主。

3 主导冲击启动的宽煤柱静载荷震波CT原位探测

基于静载荷在宽煤柱冲击地压启动过程中的主导作用[12]，并且空间赋存状态较为明确，因此采用震波CT原位探测技术实际揭示宽煤柱静载荷分布是评估其冲击危险性、制定针对性防冲方案的基础。

3.1 震波CT原位探测理论基础

震波CT(震波层析成像)原位探测技术采用波兰PASAT-M型便携式微震系统，通过检波器接收人为爆破激发的穿越工作面煤岩体的纵波走时信号。纵波在走时成像的情况下以射线形式在探测区域介质内传播，将“激发点-检波器”包络的探测区域划分为一系列小矩形网格，通过高频近似反演，走时成像公式可表示为

(1)

式中，ti为纵波走时,ms;dij为第i条射线在第j个网格中的射线路径长度;M为射线总数;N为网格数量。

表示为矩阵方程形式为

T=D·S

(2)

式中，T为纵波走时列向量;S为慢度向量;D为射线长度矩阵。

通过联合迭代重建算法(SIRT)进行速度场图像重建，可获得煤岩体内部波速分布情况。

研究表明，冲击地压多发生在高应力区及应力异常区，对应速度场图像中纵波波速区及波速梯度区，因此，建立以波速异常系数CA和波速梯度系数CG为主要因子的冲击地压危险性评估模型[13]:

(3)

利用PASAT-SSA震波CT后处理软件进行反演计算，最终形成冲击危险性指数C值分布云图，C值与冲击危险等级对应标准[13]见表1。当C<0.25时，无需处理，当C≥0.25时，需根据不同危险等级进行相应的卸压及解危处理。震波CT原位探测技术以穿透煤岩体的实际震动波射线进行波速反演，可有效反映实际条件下煤岩体内静载荷分布特征及结构特性，实现大范围静载荷探测，较常规解析和数值方法更真实、全面，但该技术观测系统的布置受到井下采掘空间、炸药激发震动波能量等客观条件限制[12]。

3.2 震波CT原位探测方案

为掌握80 m宽煤柱的静载荷实际分布状态，考虑1301工作面运输巷冲击地压显现区域分布情况、工作面条件及设备探测能力，采用近完全观测系统观测方式(两边布置激发震源、一边布置接收探头)可提高数据覆盖率，能有效提高探测精度。因此，冲击发生后，在煤柱内钻取65 m深孔，孔内按炮间距5 m布置激发震源炮至孔口，并在1301运输巷按炮间距10 m、孔深2 m补充设计激发炮，每炮装药量为200 g，一炮一放，短断触发;采集端检波器布置在80 m宽煤柱内原有2号联络巷，道间距为7 m，观测系统布置如图5所示。设定检波器采样频率为2 000 Hz，增益20 dB，采样长度0.5 s，每次激发有11通道检波器同时接收信号。

表1冲击危险等级分类
Table1Classificationofrockbursthazardlevel

C值<0.250.25～0.500.50～0.750.75～1.00等级无弱中强

图5 宽煤柱冲击危险指数分布Fig.5 Distribution of burst hazard indices in wide coal pillar

3.3 宽煤柱冲击危险性评估

图5为80 m宽煤柱探测区域冲击危险性指数C分布图，图中以蓝色到红色从小到大来代表探测区域内冲击危险性指数，区域内C最大值为0.7，最小值为-0.5。可以看出:

(1)宽煤柱测区内约2/3区域冲击危险指数C=0.5～0.7，对照表1，表明80 m宽煤柱仍然存在静载荷集中区域，整体具有中等冲击危险。

(2)宽煤柱运输巷侧冲击危险指数较采空区侧高，表明采空区侧煤柱已塑性破坏，承载能力低，导致应力向运输巷侧转移，静载荷集中程度高。

(3)在冲击震源区域，-0.50

4 基于静载荷疏导的多层次防控技术

基于震波CT探测反演结果，为阻止宽煤柱静载荷型冲击地压发生启动，关键在于避免畸高静载荷的产生，削弱宽煤柱弹性能积聚水平，实现的主要技术途径有以下3种:

(1)大直径钻孔预卸压。

在1301运输巷两帮超前工作面350 m范围内实施大直径钻孔卸压，钻孔直径125 mm，间距1 m，工作面侧钻孔深度25 m，煤柱侧钻孔深度为50 m，使得应力集中区向深部转移，耗散弹性核区积聚的弹性变性能[14]。

(2)降低开采扰动。

统计对比分析1301工作面微震活动规律与推进度的关系如图6所示，日进尺为2.4 m时微震频次和能量较为平稳，日进尺为3.2 m时微震活动较为剧烈，因此，确定合理日进尺为2.4 m，均匀推进。

图6 微震活动与推进度关系Fig.6 Relationship of microseismic activities and advancement

(3)巷道全断面补强支护。

对1301运输巷冲击地压显现及巷道破碎、变形严重的区域进行返修处理，巷道断面扩刷至设计尺寸，支护布置如图7所示，顶板采用φ22 mm×2 500 mm高强左旋螺纹钢锚杆，间排距900 mm×1 000 mm;采用φ21.6 mm×7 300 mm锚索配合L4.0 m T型钢带加强支护，锚索间排距1 800 mm×1 300 mm，钢带顺巷道方向布置，间距1.3 m。两帮采用φ20×2 500 mm等强右旋螺纹钢锚杆，间排距700 mm×2 000 mm;采用φ21.6 mm×4 300 mm锚索加强支护，间排距1 200 mm×2 000 mm;锚杆索均配合L1.6 m W型钢带支护，改良围岩应力状态，提高围岩抗冲击能力[15-17]，降低冲击地压显现程度。

图7 巷道补强支护设计Fig.7 Support reinforcement design

由于大直径钻孔预卸压造成宽煤柱深部煤体松散破碎，施工的65 m爆破激发孔塌孔严重无法布置激发震源，未能实现卸压后的震波CT探测检验。在宽煤柱布置了应力及微震监测系统，回采过程中监测应力未发生突增，微震能量及频次变化平缓，目前，1301工作面已安全回采宽煤柱区。

需要指出的是，上述防冲方案只是针对已经出现冲击危险的区段煤柱采取的补救性措施。要从根本上避免或消除区段煤柱冲击地压，应以静载荷区域疏导为防治理念[18]，避免煤柱内形成高应力集中为原则，在工作面开采设计阶段超前防范，采用小煤柱或无煤柱护巷技术是解决煤柱冲击地压的有效路径。

5 结 论

(1)高自重应力、强构造应力为冲击启动提供基础静载荷，采空区后方持续增长的侧向支承压力提供了增量静载荷，当两者叠加导致宽煤柱垂直应力超过冲击临界支承压力时，为静载荷型冲击地压的发生提供了力源条件，因此，防治应以调控静载荷为主。

(2)基于静载荷是冲击地压启动的基础，采用震波CT原位探测技术反演得到宽煤柱冲击危险指数C=0.5～0.7，表明冲击发生后，宽煤柱仍然存在静载荷集中区域，具有中等冲击危险。

(3)为阻止宽煤柱静载荷型冲击地压发生启动，制定了“大直径钻孔预卸压-降低开采扰动-巷道全断面补强支护”的多层次防冲技术，现场监测应力未发生突增，微震能量及频次变化平缓。

(4)以静载荷区域疏导为理念，探索采用小煤柱或无煤柱护巷技术是防治煤柱冲击地压的有效路径。