Discussion on key problems in prevention and control system of coal mine rock burst
-
摘要:
针对我国煤矿冲击地压的防治现状，对煤矿冲击地压防治体系中的几个关键问题进行了探讨，得到以下主要结论：从灾害控制的角度界定了冲击地压与矿震及其相互关系，通过主导力源与煤层-围岩结构的分类组合，提出了32类冲击地压和5类矿震的分类，明确了冲击地压以消除灾害为治理目标，而矿震则需要根据致灾程度分级治理；提出并量化了基础应力影响指数、附加应力影响指数、冲击倾向性影响指数和支护强度影响指数4项冲击危险性评价指标，设计了包含应力条件、冲击倾向性条件和支护条件的冲击危险性多参量耦合评价方法；探讨了震动场、应力场和位移场3级冲击危险监测预警机理，提出了单参量特殊预警与多参量常规预警相结合的预警机制，设计了掘进影响区、回采影响区和不受采掘影响区的分区联合监测方法；根据多元评价指标中降低各因素冲击危险性的要求，提出了“低应力”“低扰动”“低倾向”“强支护”的冲击地压防治路线，设计了基于应力和围岩结构分类的局部卸压和解危措施的组合技术；针对矿震诱发井下冲击地压和地面震动损害的动力灾害，将此类灾害的治理等级划分为预防治理、对症治理和源头治理3级，设计了源头减震、分类治灾、卸压抗震的多级防控技术体系。
Abstract:According to the present situation of prevention and control of coal mine rock burst in China, several key problems in the prevention and control system of coal mine rock burst are discussed, and the following main conclusions are obtained: from the perspective of disaster control, rock burst and mine earthquake and their relationship are defined, and 32 types of rock burst and 5 types of mine earthquake are put forward through the classification and combination of main force source and coal seam-surrounding rock structure, and it is clear that the goal of rock burst is to eliminate disasters, while mine earthquake needs to be treated according to the degree of disaster; Four impact risk evaluation indexes, namely, foundation stress impact index, additional stress impact index, impact tendency impact index and support strength impact index, are quantified, and a multi-parameter coupling evaluation method of impact risk including stress condition, impact tendency condition and support condition is designed. This paper probes into the monitoring and early warning mechanism of three-level impact hazards in vibration field, stress field and displacement field, puts forward an early warning mechanism combining single-parameter special warning with multi-parameter conventional warning, and designs a joint monitoring method for the excavation-affected area, mining-affected area and non-excavation-affected area. According to the requirements of reducing the impact risk of various factors in multiple evaluation indexes, the prevention and control routes of rock burst with low stress, low disturbance, low inclination and strong support are put forward, and the combination technology of local pressure relief and risk relief measures based on stress and surrounding rock structure classification is designed. In view of the dynamic disasters caused by underground rock burst and ground vibration caused by mine earthquake, the governance levels of such disasters are divided into three levels: prevention and control, symptomatic control and source control, and a multi-level prevention and control technology system of source shock absorption, classified disaster treatment, pressure relief and earthquake resistance is designed.
-
0. 引 言
我国中东部煤炭资源逐渐濒临枯竭，煤炭开采战略重心逐渐向西部地区转移[1]，随着采煤技术与装备水平的发展，矿井开采强度逐渐增大，开采扰动强度过大时产生的动力灾害更具冲击性和突然性[2]，致使冲击地压已成为目前威胁煤矿安全最为突出的动力灾害之一[3]。而影响冲击地压的因素众多，矿井开采强度作为众多开采技术影响因素之一，近些年来随着对冲击地压研究工作的深入与细化，逐渐受到重视。开采强度的不同主要体现在煤层采高、工作面宽度、推进速度等参数的差异。在采高、工作面宽度等确定的情况下，主要受推进速度的影响，特别是近距离煤层开采存在实体煤下单次采动与采空区下重复采动的复杂空间影响，煤岩赋存条件特殊，不同采动强度下动静载叠加诱发强矿压显现的机制复杂，因此亟需开展冲击地压矿井近距离特厚煤层综放工作面推进速度研究。近距离煤层指井田开采范围内相邻两层煤的层间距离较小，开采时相互影响的煤层[4]，一般采用下行开采，在上层煤开采结束后，其上部覆岩垮落，下层煤二次开采造成重复采动，对覆岩结构的整体性造成严重影响[5]。在工作面推进速度研究方面，现场实践发现，工作面高速推进易诱发冲击地压等动力灾害，诸多学者根据煤层资源赋存条件等取得不同的研究成果，其中张宏伟等[6]认为随推进速度提高，煤体应力升高速率加快，煤体脆性增强，进而导致煤岩冲击倾向性增强。王家臣等[7]认为工作面推进速度提高，基本顶断前积聚的弹性应变能及初始动能增加，发生破断失稳的可能性越大。谭云亮等[8]分析不同开采因素与所产生动能之间关系，得到开采进尺与所产生的动能之间呈正相关关系。冯龙飞等[9]研究发现推进速度越快围岩单位时间内积聚的能量越多，结合微震监测分析得到大采高工作面安全临界推进速度为4 m/d。笔者[10-11]通过对采动影响的微震监测分析，得到宽沟煤矿B4-1煤层推进速度为4 m/d时较为合宜。谢广祥等[12]认为推进速度越快，超前支承压力峰值越大，且峰值位置距离煤壁越近。闫宪磊[13]指出推进速度越快，煤岩发生冲击地压危险性越高。刘金海等[14]揭示了工作面高速推采、非匀速开采时冲击危险性较大，易诱发冲击地压。赵同彬等[15]通过微震活动表现出的低能高频现象认为推进速度为4.0 m/d左右时，更有利于覆岩能量释放。马海峰等[16]发现随着推进速度的增加，煤岩的峰值强度呈不断增大趋势。崔峰等[17]基于加卸载响应比研究认为冲击地压矿井急倾斜巨厚煤层工作面推进速度为4.8 m/d时能较好实现安全高效开采。
上述学者通过理论分析、模拟试验、现场监测等不同方法，对工作面推进速度影响下的采场围岩矿压及冲击危险性进行了有益的探索，为后续相关研究奠定了坚实的基础，但针对冲击地压矿井近距离煤层综放工作面推进速度的研究相对较少。因此，采用数值模拟、理论分析等方法，对冲击地压矿井近距离煤层综放工作面合理推进速度展开综合研究，为类似工程条件下冲击地压灾害防治提供理论基础。
1. 工程概况及数值模型建立
1.1 工程概况
宽沟煤矿隶属国能新疆宽沟矿业有限责任公司，为典型的冲击地压矿井，行政区隶属昌吉州呼图壁县管辖，井田东西长约9.7 km，南北宽约3.15 km，井田面积约20.132 5 km2，矿井主要赋存B4-1和B2两层煤，其中B4-1煤层具有强冲击倾向性，顶板具有弱冲击倾向性，煤层平均厚度3.0 m，平均倾角14°，其主采的W1145工作面倾向宽度约175 m，走向长度约746 m，目前已回采结束。在B4-1煤层下方约43.8 m的B2煤层具有弱冲击倾向性，顶板具有强冲击倾向性，煤层平均埋深约392 m，煤层平均厚度9.5 m，煤层倾角为12°～14°，顶板多为粉砂质泥岩、粉砂岩；底板多为粉砂岩、中砂岩，结构简单，不含夹矸，属稳定特厚煤层。其主采的W1123工作面倾向宽度约192 m，走向长度约1468 m，采用综放开采，采高3.2 m，放煤高度6.3 m，采放比为1∶1.97，B4-1煤层与B2煤层工作面分布情况如图1所示。
由图1可知，W1123工作面回采面临着在上覆B4-1煤层实体煤下单次采动与W1145工作面采空区下重复采动的复杂环境，上覆B4-1煤层实体煤与采空区交界作为边界支承点，承受着采空区覆岩破断后传递的应力，在实体煤一定区域内形成支承压力升高区，在其下方岩层一定范围内形成应力集中区，随着B2煤层W1123工作面的开采，工作面的采动应力和W1145采空区交界实体煤引起的应力集中区叠加，形成一个新的应力集中区，该区域的形成在一定程度上会加剧下部W1123工作面矿压显现强度，煤层交界覆岩结构特征如图2所示。
![]()
图 2 采动覆岩结构及应力分布Figure 2. Overlying rock structure and stress distribution during mining1.2 数值分析模型的构建及方案设计
鉴于3DEC数值模拟分析能够很好反映工作面覆岩运移垮落、能量积聚特征等，根据宽沟煤矿地质条件，借助3DEC数值模拟软件，建立并基于现场实测数据校核了数值模型，采用理想弹塑性本构模型莫尔−库仑准则判断岩体的破坏。
通过数值模拟试验，对W1123工作面不同推进速度下的采场围岩矿压规律及冲击危险性进行研究，其W1123工作面区域钻孔局部综合柱状（ZK201）如图3所示。根据现场地质调查与岩石力学试验结果综合确定煤岩物理力学参数，数值模拟实验中煤层顶底板附近区域主要煤岩物理力学参数见表1。
![]()
图 3 W1123工作面局部钻孔综合柱状（ZK201）Figure 3. Comprehensive histogram of local drilling in W1123 working face (ZK201)表 1 部分岩层煤岩力学参数Table 1. Mechanical parameters of rock strata and coal rocks岩性 密度/(kg·m−3) 抗拉强度/MPa 弹性模量/MPa 抗压强度/MPa 泊松比 体积模量/GPa 剪切模量/GPa 黏聚力/MPa 内摩擦角/(°) 砂砾岩 2 532 2.33 25 400 34.32 0.27 14.17 8.50 8.92 28 B4-1煤 1 303 1.97 7 485 24.80 0.35 6.25 2.88 3.81 21 粗砂岩 2 541 4.48 25 820 48.23 0.19 13.45 10.94 9.73 29 泥岩 2 533 2.28 15 420 37.65 0.20 7.42 3.83 4.39 30 细砂岩 2 592 3.17 33 420 57.60 0.18 19.57 14.07 6.58 27 B2煤 1 640 2.15 4 923 26.76 0.36 4.92 2.74 4.30 24 建立宽沟煤矿走向数值模拟模型的尺寸为：长×宽×高为1 600 m×60 m×378 m，采用摩尔−库伦准则判断煤岩体的破坏，将模型顶部未模拟岩层换算成5.4 MPa的等效载荷施加于模型顶部，四周和底边限制变形，得到W1123工作面走向数值计算模型如图4所示，为去除边界效应，B1煤层及B2煤层边界煤柱均留设65 m。
设计模拟方案采用一次采全高的开采方法，对于冲击地压矿井的推进速度研究，业内学者[7,9-10,15]大多采用4 m/d作为安全或中等推进速度。随着矿井机械化水平的提高及智能监测设备的应用，一定程度上实现了对冲击地压等事故的精准防控，在保证安全生产的前提下，部分矿井将提高工作面推进速度作为提升矿井生产能力的首选措施。笔者对此加以充分考虑，在综合考虑上述因素的基础上，此次将推进速度暂定为4.0 m/d以下、4.0～8.0 m/d，8.0 m/d以上3个区域范围，将推进速度划分为低速、中速、高速3种推进速度。同时为对比工作面推进至同一位置时推进速度对采动应力、能量积聚、覆岩破断特征及煤岩冲击危险性的影响规律，结合宽沟煤矿采煤机滚筒0.8 m的单次截深，设计3.2、6.4、9.6 m/d作为低速、中速、高速推进速度，根据设计的3种推进速度，将模型W1123工作面煤层走向开挖区域网格划分采用3.2、6.4、9.6 m/d推进速度下的最大公因子3.2 m，以此来保证同一数值模型进行不同推进速度下试验效果分析，进一步提高了数值模拟的准确性，在靠近煤层附近区域的顶板岩层中，为了使数值模拟效果更接近于实际，更好地模拟出覆岩运移破断等效果，块体尺寸相对较小，随着逐渐远离煤层，块体尺寸略微增大。并据此展开近距离煤层综放工作面不同推进速度对采动应力、能量积聚、覆岩破断特征及冲击危险性分析。
2. 不同推进速度对采动应力变化与覆岩能量积聚特征的影响
地下煤体开采对煤体来说是卸载过程，相反对周围岩体是加载过程，推进速度作为开采强度的主要影响因素之一，其对煤岩体内部应力、能量演化影响显著，为此本次采用数值模拟试验分析不同推进速度下的超前支承压力分布、覆岩的弹性能积聚等规律。
2.1 推进速度对采动应力的影响
由于W1123工作面煤层赋存空间复杂，地下开采活动的进行造成采场围岩应力的转移及重新分布，破坏了原有的应力平衡状态，推进速度越快，周围岩体的加载速率越大[18]。推进速度提高导致周围岩体加卸载活动剧烈，加剧了围岩局部应力集中现象，易诱发覆岩发生破断失稳。
为研究工作面在实体煤下回采和采空区下回采的超前支承压力演化规律，分析了工作面以不同推进速度推进至相同位置时工作面前方的顶板应力演化规律，并选取工作面推进至280、600、800、1 000、1 300 m五个位置，绘制出超前支承压力演化趋势如图5所示。
![]()
图 5 不同推进速度超前支承压力演化规律Figure 5. Evolution law of advance abutment pressure at different advancing speeds由图5可知，工作面分别在实体煤与采空区下回采时，超前支承压力峰值均随工作面向前推进不断增大，其影响范围随工作面的推进不断向前推移，峰值应力随之向前推移，且随推进速度的提高，超前支承压力峰值也随之不断增大。实体煤下超前支承压力明显大于采空区下，在接近上覆B1煤层W1145工作面开切眼后方100 m左右范围内出现了较为明显的应力集中现象，超前支承分布近似呈现“双峰值”现象。主要是由于上覆W1145工作面回采后，其开切眼左侧实体煤承载了部分未垮落覆岩结构以及边界未压实的铰接结构的应力，W1123回采形成的采动应力与W1145工作面采空区的残余支承压力叠加，致使应力水平升高，当工作面推进至上覆W1145开切眼附近时，超前支承压力峰值快速上升，峰值应力达到最大。
工作面在实体煤下回采时，随着工作面推进速度的不断提高，超前支承压力峰值不断增大。从图5a可以看出，当工作面以3.2 m/d回采时，超前支承压力逐渐由40.85 MPa增加至55.56 MPa，增加了14.71 MPa，增幅约36%；从图5b推进速度提高至6.4 m/d时可以看出，超前支承压力逐渐由46.94 MPa增加至74.48 MPa，增加了27.54 MPa，增幅约58.67%；从图5c推进速度提高至9.6 m/d时可以看出，超前支承压力逐渐由50.35 MPa增加至111.36 MPa，增加了61.01 MPa，增幅约121.17%。结合图5a、图5b、图5c可以看出工作面推进至同一位置时的超前支承压力随推进速度的提高明显升高，当工作面推进至280、600、800 m时，超前支承压力峰值随推进速度的提高分别从40.85 MPa增加至50.35，44.39 MPa加至55.50，55.56 MPa增加至111.36 MPa，3个不同位置的超前支承压力峰值随推进速度的提升其幅值分别增加了23.26%、25.03%、100.43%。
而在采空区下回采时，由于W1145采空区的存在，对采空区与W1123工作面中间的岩层起到一定程度的卸荷作用，使得W1123工作面上覆岩层承担的应力相对比较小，当工作面以3.2 m/d回采时，超前支承压力逐渐由24.86 MPa增加至32.90 MPa，增加了8.04 MPa，增幅约32.34%；当推进速度提高至6.4 m/d时，超前支承压力逐渐由37.49 MPa增加至41.18 MPa，增加了3.69 MPa，增幅约9.84%；当推进速度提高至9.6 m/d时，超前支承压力逐渐由42.33 MPa增加至43.26 MPa，增加了0.93 MPa，增幅约2.2%。
工作面推进至不同位置时，超前支承压力峰值与不同推进速度的关系如图6所示。由图6可知，工作面以3.2、6.4、9.6 m/d的推进速度推进至同一位置时的平均峰值应力分别为39.75、50.35、60.56 MPa，中速（6.4 m/d）、高速（9.6 m/d）推进速度下的峰值应力比低速分别增加约26.67%、52.35%。峰值应力随着推进速度的增加呈明显的非线性递增趋势，且其增幅逐渐增加，其非线性表达式为：
$ y=335.101\;57- 306.337\;71 {\rm{e}}^{(-0.011\;42 v)}$ 。![]()
图 6 峰值应力与不同推进速度的关系Figure 6. Relationship between peak stress and different advancing speed从上述分析可知，推进速度越快，超前支承压力峰值愈大，峰值应力随推进速度增加呈明显的非线性递增趋势，且其增幅逐渐增加。在实体煤下回采时，超前支承压力峰值随推进速度的提高不断增大，其增幅随之不断增加，推进速度较低时，采场围岩内部应力水平相对较为缓和。而在采空区下回采时，超前支承压力较实体煤下明显降低，随推进速度的提高，峰值应力水平呈现略微升高的趋势，但与在实体煤下回采时的应力水平相比还相差甚远，其增幅随推进速度的提高呈现逐渐降低的趋势。
2.2 推进速度对覆岩能量积聚的影响
推进速度的不同造成上覆顶板岩层的应力有所不同，进而导致覆岩内能量变化不一，煤岩体内的能量聚集和突然释放容易诱发动力灾害事故。相关研究表明覆岩破断一般超前于工作面发生[19]。因此，绘制如图7所示工作面回采过程中的覆岩结构演化特征。
![]()
图 7 工作面覆岩结构特征Figure 7. Structural characteristics of overlying strata in working face同时建立如图8所示的顶板岩层简支梁模型。
其弯矩M(x)为
$$ M(x)=\frac{q b x}{2}(L-x) $$ (1) 式中：q为上覆岩层均布荷载，MPa；x为截面到支点的距离，m；L为极限长度，m；b为岩梁宽度，取1 m。
通过积分计算得到该简支梁积聚的弹性能W为
$$ W=\int_0^L \frac{M^2(X)}{2 E I} {\rm{d}} x=\frac{q^2 L^5 b^2}{240 E I} $$ (2) 式中：E为岩层的弹性模量，GPa；I为转动惯量kg·m2。
由式（2）可知，覆岩积聚的能量与覆岩悬空长度呈正相关性，而覆岩的极限悬空长度与推进速度密切相关，推进速度越快，覆岩悬空长度越长，随着覆岩悬空长度的增加，覆岩积聚的弹性能显著增加。结合表1所示岩层物理力学参数，根据3DEC数值模拟试验研究得出不同推进速度下顶板的初次垮落破断距分别为54.4、64、76.8 m，取其上方均布载荷q=8 MPa，转动惯量I=h3/12，并计算覆岩初次破断时的弹性能见表2。
表 2 不同推进速度覆岩初次破断距及弹性能特征Table 2. Initial breaking distance and elastic properties of overlying strata at different advancing speeds推进速度/（m•d−1） 岩层 初次破断距/m 弹性能/kJ 3.2 厚硬直接顶 54.4 11.79 6.4 厚硬直接顶 64.0 26.57 9.6 厚硬直接顶 76.8 66.11 从表2中得出，随着推进速度的不断提高，覆岩顶板初次垮落步距不断增加，覆岩内积聚的弹性能随之不断增加，工作面以3.2、6.4、9.6 m/d推进时覆岩初次破断时覆岩内积聚的弹性能分别为11.79、26.57、66.11 kJ。中速（6.4 m/d）、高速（9.6 m/d）推进速度下的覆岩弹性能分别是低速（3.2 m/d）时的2.25倍、5.61倍。
同时结合三向应力状态下岩体应变能密度公式[20]，计算煤岩体单位积内积聚的弹性能：
$$ U=\frac{1}{2 E}\left[\sigma_1{ }^2+\sigma_2{ }^2+\sigma_3{ }^2-2 \mu\left(\sigma_1 \sigma_2+\sigma_2 \sigma_3+\sigma_1 \sigma_3\right)\right] $$ (3) 式中：U为岩体弹性势能，J/m3；
$ \sigma_1 $ ，$ \sigma_{2} $ ，$ \sigma_{3} $ 为岩体三向应力，Pa；μ为岩体的泊松比。将W1123工作面不同推进速度推进至覆岩初次破断时覆岩应力及煤岩物理力学参数代入式（3），计算并绘制出如图9所示不同推进速度覆岩初次破断时覆岩应变能密度演化趋势。
![]()
图 9 不同推进速度覆岩初次破断时覆岩应变能密度演化趋势Figure 9. Evolution trend of strain energy density of overlying strata during initial fracture of overlying strata at different advancing speeds从图9中可以看出，不同推进速度推进至覆岩初次破断时，应变能密度演化趋势呈相似性，且实体煤下应变能密度远高于上覆W1145采空区下。采空区内由于应力释放导致应变能密度极低，当W1123工作面以3.2 m/d推进时在工作面煤壁前方应变能密度出现激增现象，应变能密度为43.82 kJ/m3，而后出现缓慢下降趋势，在上覆W1145工作面开切眼后方约100 m位置时应变能不断增加，在距离开切眼约50 m时应变能再次出现激增现象，应变能密度为30.27 kJ/m3，与实体煤下超前支承压力分布规律的“双峰值”现象相对应，在实体煤下应变能密度呈现出先减后增的演化趋势；当推进速度提高至6.4 m/d时，煤壁前方应变能密度激增至53.10 kJ/m3，在距离W1145开切眼约50 m时应变能密度激增至36.10 kJ/m3；而当推进速度提高至9.6 m/d时，煤壁前方应变能密度激增至64.71 kJ/m3，在距离W1145开切眼约50 m时应变能密度激增至37.17 kJ/m3。在上覆W1145开切眼附近，应变能密度均出现骤降现象，而后缓慢增长，在停采线附近出现缓慢下降趋势，在W1145采空区下应变能密度均呈现出先增→平稳过渡→后减的演化趋势。
综上可知，随着工作面推进速度的提高，覆岩初次破断距越大，顶板悬空距离越长，覆岩内积聚的能量不断增加，初次破断时应变能密度随之不断增加，同时工作面不同推进速度推进至覆岩初次破断时的应变能密度演化趋势可划分为采空区应变能密度降低阶段、工作面前方应变能密度激增、实体煤下先降低后缓慢增长阶段、W1145采空区下先增后减的演化趋势。
3. 不同推进速度的覆岩运移破断特征及冲击危险性分析
不同推进速度对采动应力与覆岩能量积聚影响较为显著，随着煤岩体采出后，采场应力演化与能量积聚造成覆岩集中破断，引起工作面的矿压显现，其显现的强度与覆岩运移特征密切相关，覆岩运移破断是引起动力失稳、造成冲击的主要原因，因此分析不同推进速度下覆岩破断特征及煤岩冲击危险性对于矿井安全、高效生产尤为重要。
3.1 推进速度对覆岩运移破断特征的影响
随着煤层采出后，其覆岩空间演化特征普遍发生变化，主要表现在围岩开始活动，引起覆岩结构上的变形、破坏和相对位置上的改变，随着开采的不断扰动与破坏，进一步造成覆岩运移破断，而工作面的矿压显现主要是由采动覆岩集中破断运动引起的，为研究推进速度对覆岩破断的影响，选取3种速度推进至相同位置（以推进至141.8 m时）覆岩演化情况进行分析。笔者以位移云图分析覆岩运移与贯通情况，进而得到覆岩变形破坏情况。
由图10工作面不同推进速度的覆岩破坏特征可以看出，当工作面以3.2 m/d的速度推进至141.8 m时，采空区覆岩中部向上产生近似三角形状的发育趋势，最大位移约2.75 m，距离顶板愈远，位移变化量越小；当推进速度提高至6.4 m/d时，采空区左侧覆岩出现明显下沉，覆岩下沉进一步向上扩展，区域进一步加大，采空区中部下沉较为明显的区域不断向上扩展，最大位移为4.68 m；而当速度提高至9.6 m/d时，采空区左侧覆岩下沉愈加强烈，除了采空区覆岩顶板位移变化强烈外，上覆岩层位移变化不断增大，并与周围岩体相互贯通，最大变化量为6.86 m；从位移云图可以看出，采空区右侧边界与上覆W1145工作面开切眼形成近似半椭球形包络线，当推进速度为3.2、6.4 m/d时，采空区覆岩右侧端部与W1145工作面开切眼形成相连通的2个层次分明的半椭球形包络线，当推进速度提高至9.6 m/d时，采空区覆岩上部仅剩一个与W1145工作面开切眼相连通的包络线，且覆岩变化量大、运移破断剧烈，呈现出明显的非连续变形发育特征，采空区下部覆岩运移演化快速向上传递并逐渐贯通。
![]()
图 10 工作面不同推进速度的覆岩破坏特征Figure 10. Failure characteristics of overlying strata at different advancing speeds in working face综上分析，随着推进速度的增加，覆岩位移变化量不断增大，即覆岩破坏高度及范围进一步扩大，呈现出明显的非连续变形发育特征。由于W1123工作面从开切眼开始，需经历实体煤下−上覆W1145工作面采空区下2个阶段回采过程，因此，选取不同推进速度实体煤下的覆岩破断形态进行对比分析，不同推进速度下W1123工作面在实体煤下回采时覆岩破断形态如图11所示。
![]()
图 11 工作面在实体煤下回采覆岩破坏特征Figure 11. Failure characteristics of overlying strata in mining face under solid coal从图11可知，工作面在实体煤下回采结束时，覆岩破断下沉均呈现出梯形结构，在采空区覆岩两端中下部位移变化量交替区域（红、绿色变化交替区）裂隙较为发育，随着推进速度的增加，覆岩破断范围进一步加大，其主要体现在采空区右上方覆岩颜色交替变化区域内，推进速度由3.2 m/d逐渐增加到9.6 m/d时，采空区右侧区域内（方框及椭圆圈内）的位移进一步加大，范围进一步扩展，说明随着推进速度提高，覆岩破断范围进一步加大，导致采空区覆岩变化量进一步加大，覆岩运移愈加强烈，直至W1123工作面回采结束，覆岩破断形态如图12所示。
![]()
图 12 工作面回采结束时覆岩结构破断形态Figure 12. Fracture morphology of overlying rock structure at the end of mining in working face从图12中可以看出，覆岩破断下沉呈现出梯形结构，同时由于上覆B1煤层W1145工作面的存在，使得模型右侧区域覆岩的变形明显高于实体煤下，在采空区两侧中下部位移变化量交替区域，裂隙离层较为发育。随着推进速度的提高，覆岩右侧端部及W1145开切眼附近区域上覆岩层位移波动较大的范围及变化量进一步加大，当推进速度由3.2 m/d增加至9.6 m/d时，右侧端部覆岩位移变化量逐渐加大，W1145开切眼附近上覆岩层位移变化量也随推进速度的增加不断增大，且由于开切眼的存在，导致变化量较大的区域不断向左边实体煤上方覆岩扩展，进一步说明随着推进速度的增加，覆岩位移变化量大、破断范围广，覆岩运移剧烈的特点。
3.2 不同推进速度冲击危险性分析
煤岩冲击倾向性是评估及预测冲击地压危险程度的主要依据之一[21]。由于宽沟煤矿B2煤层顶板具有强冲击倾向性，通常煤岩体冲击倾向性越强，煤层开采引起采场应力重新分布，易在局部形成高度应力集中现象，越易诱发冲击地压。
业界学者公认的冲击危险性工程判据可概括为“引起冲击的应力至少应大于煤岩单轴抗压强度且煤岩具有冲击倾向性”，通常将围岩外应力与围岩抗压强度的比值记为Im，引起冲击地压发生的必要条件是围岩外应力超过煤岩单轴抗压强度，即Im>1。将应力比Im在1～1.5，1.5～3、大于3时相对应的冲击危险性划分为弱、中等、强冲击3种危险性[22]。为方便表述，将应力比值Im称为冲击危险性指数。计算公式如下：
$$ I_{\mathrm{m}}=\frac{\sigma_{\mathrm{p}}}{\sigma_{\rm{c}}}$$ (4) 式中：
$ \sigma_{\mathrm{p}} $ 为前方煤体支承压力，MPa；$\sigma_{{\rm{c}}}$ 为煤样抗压强度，MPa。通过数值模型监测工作面推进至不同区域的超前支承压力峰值，结合B2煤体的单轴抗压强度（26.76 MPa），根据式（4）计算并绘制如图13所示冲击危险性指数分布。
由图13可知，在工作面不同回采区域，由于工作面推进过程中扰动程度不一，造成不同区域冲击危险性指数呈现波浪式分布，在初始采动及采空区下回采时，冲击危险性较为缓和。当工作面以3.2 m/d的速度推进时，冲击危险性逐渐从无冲击−弱冲击−中等冲击区域过渡，在冲击高危区域，最大冲击危险性指数约为2.0；当推进速度提高至6.4 m/d时，冲击危险性从弱冲击−中等冲击区域过渡，回采全过程冲击危险性主要以中等冲击为主，在高危险冲击区域，最大冲击危险性指数接近3，即将触碰到强冲击的临界值；当推进速度提高至9.6 m时，冲击危险性明显升高，逐步从弱冲击−中等冲击−强冲击区域过渡，在接近上覆W1145开切眼附近应力集中区域时，最大冲击危险性指数约为4.2，冲击危险性异常升高。工作面在实体煤下回采时，随着推进速度的提高，冲击危险性不断增加，在距离W1145工作面开切眼约100 m时冲击危险性开始稳步上升，变化幅度较大，上覆W1145采空区与实体煤的交界区域为明显的高危险区域，工作面距离上覆W1145开切眼愈近冲击危险性越高。工作面过渡到W1145采空区下回采时，冲击危险性呈现缓慢下降的趋势，相较于实体煤下回采，采空区下回采时冲击危险性明显较低，冲击危险性随推进速度的变化幅度也明显小于实体煤下。
由上分析可知，随推进速度的提高，工作面冲击危险性明显增加，尤其表现在距W1145开切眼约100 m区域，冲击危险性稳步上升，在开切眼附近冲击危险性快速升高，随推进速度由3.2 m/d逐渐递增至9.6 m/d的过程中，冲击危险性指数随之由2.0递增至4.2，由中等冲击危险性跃升至强冲击危险性。
4. W1123工作面推进速度合理确定
通过上述推进速度对超前支承压力、覆岩能量积聚、覆岩运移破断特征及冲击危险性分析的研究发现：推进速度越快，超前支承压力峰值愈大，采动影响范围越广，峰值应力随推进速度增加呈明显的非线性递增趋势，且其增幅逐渐增加。在实体煤下回采超前支承压力峰值随推进速度的提高不断增大，其增幅随之不断增加，而在采空区下回采超前支承压力随推进速度提高呈现略微升高的趋势，其增幅随推进速度提高呈现逐渐降低的趋势；同时推进速度提高覆岩初次破断距越大，破断时释放的能量随之增加，将覆岩初次破断时应变能密度演化趋势划分为采空区应变能降低阶段、应变能激增阶段、实体煤下先降后增过渡阶段、W1145采空区下先增后降的演化趋势。同时推进速度越快，覆岩位移变化量愈大、破断范围广，运移演化剧烈，呈现出非连续变形发育特征，采空区下部覆岩的运移演化很快向上传递并相互贯通。
工作面冲击危险性随推进速度提高显著增加，实体煤下冲击危险性指数显著高于采空区下，冲击危险性随推进速度的变化也明显大于采空区下，尤其表现在距W1145开切眼约100 m区域，冲击危险性稳步上升，在开切眼附近冲击危险性快速升高，当推进速度大于6.4 m/d时，冲击危险性指数大于3，处于强冲击危险性区域，相较于中等冲击危险性，当工作面煤岩体处于强冲击危险性范畴时，煤岩体内部应力高度集中，通常造成采场围岩难以支护等情况，此时需要加大围岩卸压力度，加强支护强度，以防出现冲击地压等动力灾害。通过上文分析，汇总出近距离煤层综放工作面推进速度综合影响统计表见表3。
表 3 近距离煤层推进速度综合影响统计Table 3. Statistical of comprehensive influence of advancing speed of contugous coal seam推进速度/(m•d−1) 3.2 6.4 9.6 采动应力/MPa 39.75 50.35 60.56 初次破断能量积聚/kJ 11.79 26.57 66.11 初次破断应变能密度/(kJ·m−3) 43.82 53.10 64.71 位移变形/m
（推进至141.8 m）2.75 4.68 6.86 危险性指数 2.0 2.8 4.2 针对宽沟煤矿近距离特厚煤层综放工作面推进速度效应展开综合分析，通过研究B2煤层W1123工作面不同推进速度下超前支承压力、能量积聚、覆岩破断、煤岩冲击危险性等演化特征，综合研究认为W1123工作面在保证安全高效生产前提下的相对适宜的推进速度应不大于6.4 m/d。
5. 结 论
1）随着推进速度的提高，峰值应力不断增大，呈明显的非线性递增趋势，其增幅也随之逐渐增加。随推进速度的提高，覆岩初次破断距不断增大，覆岩能量积聚快速增加，应变能密度也随之显著增加，在煤壁前方呈现明显激增现象。
2）工作面推进至同一位置时覆岩位移变化量随推进速度增加急剧增加，破坏范围广，呈现出明显的非连续变形特征，推进速度愈快，覆岩破断范围越广，运移演化愈剧烈。
3）随推进速度的提高工作面冲击危险性随之显著增加，在实体煤下回采时的冲击危险性指数远大于采空区下，变化幅度随推进速度提高明显增大，上覆W1145采空区与实体煤的交界区域为明显的高危险区域，工作面距离上覆W1145开切眼愈近冲击危险性越高。数值模拟研究结果综合认为W1123工作面推进速度应以不大于6.4 m/d。
-
![]()
图 1 冲击地压与矿震综合分类方法
Figure 1. Comprehensive classification method of rock burst and mine earthquake
![]()
图 2 采掘空间煤层-围岩结构类型
Figure 2. Structure type of coal seam - surrounding rock in mining space
![]()
图 4 冲击地压三级监测预警原理
Figure 4. Three-level monitoring and early warning principle of rock burst
![]()
图 5 冲击危险分类分区监测预警流程
Figure 5. Monitoring and warning process of impact hazard classification and zoning
![]()
图 8 矿震诱发型动力灾害危险等级划分
Figure 8. Risk classification of dynamic disaster induced by mine earthquake
![]()
图 9 矿震诱发型动力灾害防治体系
Figure 9. Prevention and control system of mine earthquake induced dynamic disasters
表 1 Ic与冲击危险性关系
Table 1 Relationship between Ic and impact hazard
Ic <1 1≤Ic<1.5 1.5≤Ic<2 ≥2 冲击危险性 无 弱 中等 强 
下载: 导出CSV
表 2 WET与冲击危险性关系
Table 2 Relationship between WET and impact hazard
WET <2 2≤WET<3.5 3.5≤WET<5 ≥5 冲击危险性 无 弱 中等 强
下载: 导出CSV
表 3 S与冲击危险性关系
Table 3 Relationship between S and impact hazard
S/能级 <2 2≤S<6 6≤S<9 ≥9 冲击危险性 强 中等 弱 无
下载: 导出CSV
表 4 冲击危险性工程判据指标
Table 4 Engineering criterion index of impact hazard
U <0.6 0.6≤U<0.75 0.75≤U<0.9 ≥0.9 冲击危险性 无 弱 中等 强
下载: 导出CSV
表 5 基于工程分类的局部卸压、解危措施组合参考
Table 5 Reference of combination of local pressure relief and risk relief measures based on engineering classification
围岩类型 不同主导力源局部卸压、解危措施 自重应力 煤柱应力 构造应力 开采残余应力 采掘扰动应力 震动附加应力 厚煤层 沿顶结构 B/D C/D B/D B/D A/B/D B/D 全煤结构 B/D/G C/D/G B/D/G B/D/G A/B/D/G A/B/D/G 沿底结构 B/G C/G B/G B/G A/B/G A/B/G 破底半煤岩结构 B/G C/G B/G B/G A/B/G A/B/G 中厚煤层 沿顶沿底结构 B C B B A/B B 软底结构 E/F E/F E/F E/F A/E/F E/F 软顶结构 F/G F/G F/G F/G A/F/G A/F/G 薄煤层 破顶破底结构 C C C C C C 倾斜煤层 厚底煤结构 B/D/G C/D/G B/D/G A/B/D/G A/B/D/G A/B/D/G 薄底煤结构 B/G C/G B/G A/B/G A/B/G A/B/G 注：A为压裂或爆破断顶；B/C为帮部钻孔/爆破；D/E为钻孔/爆破断底；F为煤层注水；G为补强支护。
下载: 导出CSV
表 6 矿震三级治理等级划分
Table 6 Classification of three-level governance of mine earthquake
灾害等级 建筑损害 明显震感 微弱震感 无震感 强冲击 6 5 4 3 中等冲击 5 4 3 2 弱冲击 4 3 2 1 无冲击 3 2 1 0
下载: 导出CSV
-
[1] 齐庆新,王守光,李海涛,等. 冲击地压应力流理论及其数值实现[J]. 煤炭学报,2022,47(1):172−179. doi: 10.13225/j.cnki.jccs.YG21.2003 QI Qingxin,WANG Shouguang,LI Haitao,et al. Stress flow theory for coal bump and its numerical implementation[J]. Journal of China Coal Society,2022,47(1):172−179. doi: 10.13225/j.cnki.jccs.YG21.2003
[2] 窦林名,田鑫元,曹安业,等. 我国煤矿冲击地压防治现状与难题[J]. 煤炭学报,2022,47(1):152−171. doi: 10.13225/j.cnki.jccs.yg21.1873 DOU Linming,TIAN Xinyuan,CAO Anye,et al. Present situation and problem of coal mine rock burst prevention and control in China[J]. Journal of China Coal Society,2022,47(1):152−171. doi: 10.13225/j.cnki.jccs.yg21.1873
[3] 潘一山. 煤矿冲击地压扰动响应失稳理论及应用[J]. 煤炭学报,2018,43(8):2091−2098. doi: 10.13225/j.cnki.jccs.2018.0604 PAN Yishan. Disturbance response instability theory of rock burst in coal mine[J]. Journal of China Coal Society,2018,43(8):2091−2098. doi: 10.13225/j.cnki.jccs.2018.0604
[4] 潘俊锋,闫耀东,马小辉,等. 考虑时变特性的煤层大巷群冲击地压机理及防治[J]. 煤炭学报,2022,47(9):3384−3395. PAN Junfeng,YAN Yaodong,MA Xiaohui,et al. Mechanism and prevention of rock burst in coal seam roadway group considering time-varying characteristics[J]. Journal of China Coal Society,2022,47(9):3384−3395.
[5] 刘金海,姜福兴,朱斯陶,等. 典型深厚表土煤层冲击地压模式研究[J]. 煤炭学报,2020,45(5):1753−1763. doi: 10.13225/j.cnki.jccs.2020.0367 LIU Jinhai,JIANG Fuxing,ZHU Sitao,et al. Rock burst model of typical coal seam under thick alluvium[J]. Journal of China Coal Society,2020,45(5):1753−1763. doi: 10.13225/j.cnki.jccs.2020.0367
[6] 王恩元,冯俊军,张奇明,等. 冲击地压应力波作用机理[J]. 煤炭学报,2020,45(1):100−110. WANG Enyuan,FENG Junjun,ZHANG Qiming,et al. Mechanism of rockbust stress wave[J]. Journal of China Coal Society,2020,45(1):100−110.
[7] 谭云亮,张 明,徐 强,等. 坚硬顶板型冲击地压发生机理及监测预警研究[J]. 煤炭科学技术,2019,47(1):166−172. TAN Yunliang,ZHANG Ming,XU Qiang,et al. Study on occurrence mechanism and monitoring and early warning of rock burst caused by hard roof[J]. Coal Science and Technology,2019,47(1):166−172.
[8] 姜福兴,曲效成,王颜亮,等. 基于云计算的煤矿冲击地压监控预警技术研究[J]. 煤炭科学技术,2018,46(1):199−206,244. JIANG Fuxing,QU Xiaocheng,WANG Yanliang,et al. Study on monitoring & control and early warning technology of mine pressure bump based on cloud computing[J]. Coal Science and Technology,2018,46(1):199−206,244.
[9] 何学秋,陈建强,宋大钊,等. 典型近直立煤层群冲击地压机理及监测预警研究[J]. 煤炭科学技术,2021,49(6):13−22. HE Xuqiu,CHEN Jianqiang,SONG Dazhao,et al. Study on mechanism of rock burst and early warning of typical steeply inclined coal seams[J]. Coal Science and Technology,2021,49(6):13−22.
[10] 潘俊锋,王书文,刘少虹,等. 浅部矿井静载荷主导型冲击地压监测方法与实践[J]. 煤炭科学技术,2016,44(6):64−70,98. PAN Junfeng,WANG Shuwen,LIU Shaohong,et al. Practices and monitoring method of static loading dominant type mine strata pressure bump in shallow depth mine[J]. Coal Science and Technology,2016,44(6):64−70,98.
[11] 朱斯陶,姜福兴,刘金海,等. 复合厚煤层巷道掘进冲击地压机制及监测预警技术[J]. 煤炭学报,2020,45(5):1659−1670. ZHU Sitao,JIANG Fuxing,LIU Jinhai,et al. Mechanism and monitoring and early warning technology of rock burst in the heading face of compound thick coal seam[J]. Journal of China Coal Society,2020,45(5):1659−1670.
[12] 马斌文,邓志刚,赵善坤,等. 钻孔卸压防治冲击地压机理及影响因素分析[J]. 煤炭科学技术,2020,48(5):35−40. MA Binwen,DENG Zhigang,ZHAO Shankun,et al. Analysis on mechanism and influencing factors of drilling pressure relief to prevent rock burst[J]. Coal Science and Technology,2020,48(5):35−40.
[13] 姜福兴,刘 烨,刘 军,等. 冲击地压煤层局部保护层开采的减压机理研究[J]. 岩土工程学报,2019,41(2):368−375. JIANG Fuxing,LIU Ye,LIU Jun,et al. Pressure-releasing mechanism of local protective layer in coal seam with rock burst[J]. Chinese Journal of Geotechnical Engineering,2019,41(2):368−375.
[14] 曹安业,朱亮亮,杜中雨,等. 巷道底板冲击控制原理与解危技术研究[J]. 采矿与安全工程学报,2021,49(11):93−98. CAO Anye,ZHU Liangliang,DU Zhongyu,et al. Control principle and pressure-relief technique of rock burst occurred in roadway floor[J]. Journal of Mining & Safety Engineering,2021,49(11):93−98.
[15] 孔令海,邓志刚,梁开山,等. 深部煤巷顶帮控制防治冲击地压研究[J]. 煤炭科学技术,2018,46(10):83−89. KONG Linghai,DENG Zhigang,LIANG Kaishan,et al. Study on mine rock burst prevention and control with roof and sidewalls control in mine deep seam gateway[J]. Coal Science and Technology,2018,46(10):83−89.
[16] 翟明华,姜福兴,齐庆新,等. 冲击地压分类防治体系研究与应用[J]. 煤炭学报,2017,42(12):3116−3124. QU Minghua,JIANG Fuxing,QI Qingxin,et al. Research and practice of rock burst classified control system[J]. Journal of China Coal Society,2017,42(12):3116−3124.
[17] 姜耀东,潘一山,姜福兴,等. 我国煤炭开采中的冲击地压机理和防治[J]. 煤炭学报,2014,39(2):205−213. JIANG Yaodong,PAN Yishan,JIANG Fuxing,et al. State of the art review on mechanism and prevention of coal bumps in China[J]. Journal of China Coal Society,2014,39(2):205−213.
[18] 潘一山,李忠华,章梦涛. 我国冲击地压分布、类型、机理及防治研究[J]. 岩石力学与工程学报,2003,22(11):1844−1851. doi: 10.3321/j.issn:1000-6915.2003.11.019 PAN Yishan,LI Zhonghua,ZHANG Mengtao,et al. Distribution, type, mechanism and prevention of rockbrust in China[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(11):1844−1851. doi: 10.3321/j.issn:1000-6915.2003.11.019
[19] 齐庆新,陈尚本,王怀新,等. 冲击地压、岩爆、矿震的关系及其数值模拟研究[J]. 岩石力学与工程学报,2003,22(11):1852−1858. doi: 10.3321/j.issn:1000-6915.2003.11.020 QI Qingxin,CHEN Shangben,WANG Huaixin,et al. Study on the relations among coal bump, rock burst and mining tremor with numerical simulation[J]. Chinese Journal of Rock Mechanics and Engineering,2003,22(11):1852−1858. doi: 10.3321/j.issn:1000-6915.2003.11.020
[20] 窦林名,曹晋荣,曹安业,等. 煤矿矿震类型及震动波传播规律研究[J]. 煤炭科学技术,2021,49(6):23−31. DOU Linming,CAO Jinrong,CAO Anye,et al. Research on types of coal mine tremor and propagation law of shock waves[J]. Coal Science and Technology,2021,49(6):23−31.
[21] 杨光宇,姜福兴,李 琳,等. 煤矿冲击地压危险性的工程判据研究[J]. 采矿与安全工程学报,2018,35(6):1200−1207,1216. YANG Guangyu,JIANG Fuxing,CAO Anye,et al. Engineering criterion study on coal mining rock burst hazard[J]. Journal of Mining & Safety Engineering,2018,35(6):1200−1207,1216.
[22] 姜福兴,刘 懿,翟明华,等. 基于应力与围岩分类的冲击地压危险性评价研究[J]. 岩石力学与工程学报,2017,36(5):1041−1052. JIANG Fuxing,LIU Yi,QU Minghua,et al. Evaluation of rock burst hazard on the classification of stress and surrounding rock[J]. Chinese Journal of Rock Mechanics and Engineering,2017,36(5):1041−1052.
[23] 朱斯陶,刘金海,姜福兴,等. 我国煤矿顶板运动型矿震及诱发灾害分类、预测与防控[J]. 煤炭学报,2022,47(2):807−816. doi: 10.13225/j.cnki.jccs.XR21.1800 ZHU Sitao,LIU Jinhai,JIANG Fuxing,et al. Classification, prediction, prevention and control of roof movement-type mine earthquakes and induced disasters in China coal mines[J]. Journal of China Coal Society,2022,47(2):807−816. doi: 10.13225/j.cnki.jccs.XR21.1800
[24] 刘金海,孙 浩,田昭军,等. 煤矿冲击地压的推采速度效应及其动态调控[J]. 煤炭学报,2018,43(7):1858−1865. LIU Jinhai,SUN Hao,TIAN Zhaojun,et al. Effect of advance speed on rock burst in coal mines and its dynamic contral method[J]. Journal of China Coal Society,2018,43(7):1858−1865.
[25] 杨光宇,温经林,李 琳,等. 特厚煤层巷道冲击特征及冲击危险性评价方法研究[J]. 中国安全生产科学技术,2019,15(5):92−98. doi: 10.11731/j.issn.1673-193x.2019.05.015 YANG Guangyu,WEN Jinglin,LI Lin,et al. Study on rock burst characteristics of roadway in extra-thick coal seam and risk assessment method of rock burst[J]. Journal of Safety Science and Technology,2019,15(5):92−98. doi: 10.11731/j.issn.1673-193x.2019.05.015
[26] 姜福兴,杨光宇,魏全德,等. 煤矿复合动力灾害危险性实时预警平台研究与展望[J]. 煤炭学报,2018,43(2):333−339. JIANG Fuxing,YANG Guangyu,WEI Quande,et al. Study and prospect on coal mine composite dynamic disaster real-time prewarning platform[J]. Journal of China Coal Society,2018,43(2):333−339.
[27] GB 6722—2014, 爆破安全规程[S]. 北京: 中国标准出版社, 2014. [28] 高明仕,赵一超,温颖远,等. 震源扰动型巷道冲击矿压破坏力能准则及实践[J]. 煤炭学报,2016,41(4):808−814. GAO Mingshi,ZHAO Yichao,WEN Yingyuan,et al. Stress and energy criterion of the roadway destruction subjected to disturbance type rock burst and its practice[J]. Journal of China Coal Society,2016,41(4):808−814.
-
期刊类型引用(5)
1. 赵兵朝，赵阳，冯杰，王京滨，冯欣怡. 采动地表与输电线塔大板基础协同变形规律研究. 采矿与安全工程学报. 2024(03): 472-480 . 
百度学术
2. 张建，亓佳利，李长青. 厚煤层开采上覆岩层移动与应力分布特征研究. 山东煤炭科技. 2024(06): 139-143+154 . 百度学术
3. 李文福，马小辉，刘咸富，王向阳，何伟，王岩，白晓敏，刘万成. 冲击性煤层综放面采空区漏风区域识别及控制方法研究. 煤炭技术. 2024(09): 133-137 . 百度学术
4. 彭宝山，王永乐，杨学孟. 特厚弱胶结顶板煤层综放开采覆岩破坏特征与强矿压机理. 煤炭技术. 2024(10): 75-80 . 百度学术
5. 李志，杨栋，贾毅超. 深厚煤层开采过程中覆岩运移演变机理及矿压管控研究. 矿业研究与开发. 2024(12): 68-75 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 698
- HTML全文浏览量: 14
- PDF下载量: 753
- 被引次数: 6
