Study on evolution law of overburden fracture in fully-mechanized top-coal caving face with large mining height in thick coal seam
-
摘要:
大采高综放开采易形成强烈的矿压显现，上覆岩层的断裂演化规律对工作面安全生产至关重要。以羊场湾煤矿160206工作面为工程背景，运用相似模拟试验、数值模拟和理论分析的综合研究方法，对大采高综放工作面覆岩的断裂过程与覆岩运移规律进行系统研究。研究表明：导水裂隙带内岩层随工作面推进表现为“台阶下沉”，同层岩层下沉趋势沿走向表现为“急剧下降—稳定（最大值）—快速上升—稳定（最小值）”。工作面覆岩运动场由两区分布（加速下沉区、缓慢下沉区）演化为三区分布（加速下沉区、缓慢下沉区、稳定区）。对离层演化与地表下沉规律进行了定量描述，运用理论计算表达式深入地剖析了地表的动态下沉机理及其相关因素。结合相似模拟与数值模拟的试验结果，提出了覆岩断裂演化的形态变化特征：覆岩断裂形态由“单等腰梯形”演化为“双等腰梯形”，表土层影响区由“矩形”演化为“倒梯形”。分析了覆岩中垮落区、离层区、压实区、裂隙富集区的动态演化过程：垮落区逐渐增大至一定程度，高度小幅度降低并趋于稳定，离层区由下至上逐渐发育并随工作面向前移动，离层区逐渐闭合形成压实区，且压实区逐渐增大并最终保持稳定，裂隙富集区位于采空区前后端部并随工作面向前移动。
Abstract:Fully mechanized top-coal caving mining with large mining height is easy to form strong strata behavior, and the fracture evolution law of overlying strata is very important for safe production of working face. Taking the 160206 working face of Yangchangwan Coal Mine as the engineering background, the comprehensive research methods of similar simulation test, numerical simulation and theoretical analysis are used to systematically study the overburden fracture process and overburden migration law of the fully mechanized top coal caving working face with large mining height. The research shows that the strata in the water flowing fracture zone show “bench sinking” with the advance of the working face, and the subsidence trend of the same strata along the strike shows “sharp decline - stability (maximum) - rapid rise - stability (minimum)”. The overburden movement field of the working face has evolved from two areas (accelerated subsidence area and slow subsidence area) to three areas (accelerated subsidence area, slow subsidence area and stable area). The evolution of abscission layer and the law of surface subsidence are quantitatively described, and the dynamic subsidence mechanism of the surface and its related factors are deeply analyzed by using theoretical calculation expressions. Combined with the experimental results of similar simulation and numerical simulation, the morphological change characteristics of overburden fracture evolution are proposed: the overburden fracture morphology evolves from “single isosceles trapezoid” to “double isosceles trapezoid”, and the influence area of topsoil layer evolves from “rectangle” to “inverted trapezoid”. The dynamic evolution process of collapse area, separation area, compaction area and fracture enrichment area in overburden is analyzed: The caving zone gradually increases to a certain extent, the height decreases slightly and tends to be stable, the separation zone gradually develops from bottom to top and moves forward with the working face, the separation zone gradually closes to form a compaction zone, and the compaction zone gradually increases and finally remains stable, and the fracture enrichment zone is located at the front and back ends of goaf and moves forward with the working face.
-
0. 引 言
我国现有煤炭储量中，厚煤层约占45%，是我国实现高产高效开采的主力煤层[1-2]。厚煤层一般采用分层开采、放顶煤开采、大采高综采3种开采方法，特殊情况下也可同时采用2种或多种方法[3-5]。大采高综放开采是厚煤层的主要开采方法，能够实现煤炭的快速高效开采[6-10]。
随着割煤和放煤高度的增加，弯曲下沉空间显著增大，当顶板断裂时，将释放出大量的能量，从而形成强烈的矿压显现，引起工作面煤壁片帮严重、顶煤采出率低、瓦斯涌出量大等突出问题。因此，要实现厚煤层大采高综放工作面的正常安全生产，正确认识工作面覆岩的断裂演化规律尤为重要。近年来，我国诸多学者运用数值模拟、理论计算、相似模拟试验与现场实测等方法，王云广等[11]给出了高强度开采覆岩岩“两带”高度计算公式，提出了三维立体空间条件下的高强度开采覆岩运移的“弹性薄板”+“平行压力拱”复合机理模型。蒋金泉等[12]提出了离层空间演化特征的确定方法，建立了主关键层底部最大离层空间的分析模型与预测方法。许满贵等[13]研究了工作面受采动影响后采场覆岩运移规律。殷伟等[14]研究了混合综采工作面覆岩运移规律和空间结构特征，得到了不同充实率状态下混合综采工作面覆岩空间结构特征。赵毅鑫等[15]分析了8.8 m 超大采高工作面采动覆岩裂隙演化及能量耗散规律。张广超等[16]提出了基于极限拉伸应变的软弱覆岩破坏高度计算方法。杨科等[17]获得了厚硬顶板条件下厚煤层开采覆岩破断运移规律和“三带”动态演化特征。采场覆岩中对岩体活动起主要控制作用的坚硬岩层称为关键层，根据控制不同岩层可分为主关键层和亚关键层，关键层理论为覆岩破断运移规律的研究提供了新的思路[18-21]。闫少宏等[22]提出了大采高综采顶板“短悬臂梁－铰接岩梁”结构。付玉平等[23]对浅埋厚煤层大采高工作面顶板岩层断裂演化规律进行了模拟研究。柴敬等[24]对采场上覆关键层的运移进行了模拟试验并结合现场观测检测。任艳芳[25]提出了浅埋深工作面覆岩“悬臂梁－铰接岩梁”结构并进行验证。
在此基础上，根据羊场湾煤矿160206厚煤层综放工作面的来压特征，采用理论分析、相似模拟、数值模拟三者结合的方法进一步研究覆岩断裂演化规律，指导现场生产，具有重要的理论意义和工程应用价值。
1. 工程概况
160206工作面距地面平均垂深592 m，目前主开采二层煤。160206工作面的煤岩层综合柱状图如图1所示。二层煤层理发育较破碎，顶板坚硬，回采过程受采动矿压影响容易造成工作面大面积悬顶、煤壁片帮，煤层顶板富水性强，开采过程中一层煤顶板粗砂岩含水层水涌入工作面下端头。需要对160206工作面覆岩断裂演化规律进行研究，保证现场安全生产。工作面支架测点布置如图2所示。
![]()
图 1 工作面煤岩层综合柱状Figure 1. Comprehensive histogram of coal and rock strata in working face160206工作面倾斜长度为227 m，共有131台液压支架，编号分别为1、2、···，131号，从5号支架开始，每隔10台液压支架布设一组液压支架阻力监测，共计13个监测测点。
收集的液压支架数据见表1，由表1可知，160206综放工作面平均工作阻力为31.45 MPa，来压时的判据为36.31 MPa。
表 1 液压支架支护工作阻力Table 1. Working resistance of hydraulic support分站编号 P/MPa (P+σ)/MPa 动载系数 5号 30.78 35.72 1.16 15号 31.44 36.25 1.15 25号 31.90 36.45 1.14 35号 31.67 36.39 1.15 45号 31.49 36.33 1.15 55号 32.02 37.56 1.17 65号 32.37 37.82 1.17 75号 32.11 36.92 1.15 85号 31.44 36.20 1.15 95号 31.39 36.22 1.15 105号 31.36 35.94 1.15 115号 30.87 35.34 1.14 125号 30.05 34.93 1.16 平均值 31.45 36.31 1.15 注：P为工作阻力；σ为均方差。 根据支架阻力数据绘制成热图如图3所示。
由图3可知，工作面支架阻力在走向方向表现为中间高、两端低。前100 m支架阻力变化不大，之后支架阻力迅速增加，工作面推进至130.5 m，整个工作面支架阻力超过来压判据，工作面初次来压，顶板大面积垮落，初次来压步距为130.5 m。
以现场数据为基础，展开相似模拟与数值模拟研究，揭示大采高综放工作面的覆岩断裂演化规律。
2. 相似模拟试验
2.1 试验方案
2.1.1 模型设计
试验采用长×宽×高=1 800 mm×160 mm×1 300 mm的试验台。根据相似定律和试验装置，确定几何相似比$ {C_1} = 1/200 $；容重相似比$ {C_{\text{γ}} } = 0.57 $；应力相似比$ {C_{\mathrm{a}}} = {C_1}{C_{\text{γ}} } = 0.002\;9 $；时间相似比$ {C_{\mathrm{t}}} = \sqrt {{C_1}} = 0.070\;7 $。模型左右两边各留设100 mm煤柱以避免边界效应产生的影响。
2.1.2 试验过程
模型开挖高度为46.5 mm，代表井下实际煤厚9.3 m；开挖长度为1 600 mm，代表井下实际推进长度约320 m；每21 min回采一次，代表井下实际时间3 d；每次开挖进尺50 mm，代表井下实际进尺10 m。
岩层相似材料配比与关键层位置见表2。
表 2 岩层相似材料配比与关键层位置Table 2. Proportion of similar materials in rock strata and location of key layers层号 岩性 厚度/m 容重/（kN·m−3） 模拟层厚/mm 分层数 配比号 关键层位置 18 （配重）上覆岩层 402.00 — — — — — 17 粉砂岩 8.27 24.00 41.4 2 755 16 细砂岩 15.88 25.30 79.4 4 655 15 泥岩 12.56 24.20 62.8 3 855 14 中粒砂岩 20.26 24.50 101.3 5 655 13 粗砂岩 5.88 23.50 29.4 1 655 12 粉砂岩 18.46 24.00 92.3 3 755 11 中粒砂岩 7.61 24.50 38.1 2 655 10 细砂岩 11.10 25.30 55.5 2 655 9 泥岩 8.00 24.20 40.0 2 855 8 粗砂岩 55.70 23.50 278.5 9 655 主关键层 7 一层煤 1.60 14.50 8.0 1 855 6 粉砂岩 3.80 24.00 19.0 1 755 5 中粒砂岩 18.40 24.50 92.0 4 655 亚关键层 4 粉砂岩 2.40 24.00 12.0 1 755 3 二层煤 9.30 15.00 46.5 2 855 2 泥岩 3.00 24.20 15.0 1 855 1 粉砂岩 7.90 24.00 39.5 2 755 注：配比号755 表示细砂、石灰、石膏质量比为7∶0.5∶0.5；岩层总厚度210.12 m，模拟层总厚1 050.6 mm，总计分层数45层。 2.1.3 测点布置
为了精确监测煤层推进过程中，覆岩移动和位移变化情况，在模型正面不同层位布置了位移基点，位移基点布置9行，18列。对每个位移基点黏贴测标并进行编号，垂直测线由下至上编号为1～9，水平测线从左到右编号为1～18。为尽可能展示整个覆岩的运移规律，选择5条测线均布于上覆岩层。沿工作面推进方向布设3条垂直测线，分别为1号测线、2号测线、3号测线，相对开切眼距离分别为70、170、270 m。在覆岩水平方向布设2条水平测线，分别为4号测线、5号测线，4号测线布置于亚关键层中粒砂岩，5号测线布置于主关键层上方45 m处的粗砂岩。位移基点实际布置如图4所示。
2.2 覆岩断裂规律
2.2.1 基本顶初次垮落
基本顶初次垮落如图5所示。工作面推进100 m（图5a），亚关键层第1分层达到极限跨距，开始垮落。但由于采空区自由空间很大导致回转量较大，基本顶无法在采空区形成铰接稳定结构，呈现出顶板的台阶下沉。
推进至120 m（图5b），亚关键层全部垮落，断裂岩块相互铰接形成砌体梁结构，上覆岩层亦随之垮落。工作面初次来压，步距为120 m，与现场初次来压步距相对误差7.7%，误差较小，能较好地研究覆岩的断裂演化规律。
2.2.2 基本顶周期垮落
基本顶周期垮落如图6所示。推进130 m（图6a），工作面发生第1次周期来压。覆岩断裂形态呈“梯形”，断裂岩层与主关键层之间形成一条最宽为8.95 m的离层裂隙，主关键层尚无明显变化。推进270 m（图6b），基本顶发生第8次周期来压，导水裂隙带达到最大高度，采空区中部覆岩被压实，且压实区向推进方向和垂直方向扩张。推进320 m（图6c），基本顶发生第10次周期来压，导水裂隙带高度不变，压实区达到最大高度，继续向推进方向扩张。采空区两端纵向与水平裂隙发育较高，属于裂隙富集区，且跟随工作面向前移动，处于动态平衡。
2.3 覆岩运移规律
煤层开采引起覆岩的断裂失稳，导致顶板发生水平错动与竖直下沉，相应地产生水平位移与垂直位移，而水平位移相对较小，故主要对垂直位移进行监测分析。
1）垂直测线在不同推进距离下的位移（向下为负）变化趋势如图7所示。综合1、2、3号测线垂直位移变化趋势可知：
同一测点的位移变化随工作面推进表现为“台阶”下沉，不同测点的台阶数量、台阶高度、台阶宽度不同。下沉量最大值与测点编号呈负相关，即覆岩下沉量由下至上逐层递减。
亚关键层和主关键层的断裂共同决定覆岩的下沉，但不同时期主导因素不同。前期亚关键层的断裂为主导因素，中后期主关键层的断裂为主导因素。
初次来压后，由于工作面上方悬臂梁的存在，基本顶的大幅下沉不同程度地滞后于推进距离（10～20 m）。
工作面推进一定距离后，工作面后方50 m以内的覆岩尚处于下沉阶段，而距离较远的下层覆岩率先达到最大下沉量，进入压实阶段。覆岩由下至上沿推进方向逐渐压实。
2）水平测线在不同推进距离下的位移（向下为负）变化趋势如图8所示（图8a为3D曲面，图8b为水平投影）。
4号测线（图8a）位于亚关键层中粒砂岩，推进120 m时，4-2、4-3、4-4、4-5、4-6测点急剧下沉8.25～9.30 m，说明亚关键层整体断裂，工作面初次来压。工作面继续推进，测点位移逐渐向右发育，中部测点下沉量达到最大值9.3 m，处于稳定状态，采空区从左到右逐渐被压实。
5号测线（图8b）位于主关键层上方45 m处的粗砂岩，反映着主关键层所控制的岩层变化趋势。左右两端测点位移为0，测点两端到中部下沉量逐渐增大。随着工作面推进，中部测点下沉量逐渐增加，中部测点首先被压实处于稳定状态，且逐渐向右压实。
水平测线的位移变化趋势沿走向表现为急剧下降→稳定（最大值）→快速上升→稳定（最小值）。位移变化趋势在三维图中表现为一个“深渊+盆地”形状。
3. 覆岩断裂数值模拟试验
3.1 数值模型建立
基于羊场湾煤矿160206工作面实际地质条件，建立数值计算模型，如图9所示。
尺寸为520 m×1 m×289.6 m（长×宽×高），根据岩层的力学性质对网格进行大小不同的划分，模型底部为固定边界，四周为位移约束边界，模型上部距地表330 m，施加等效载荷8.25 MPa。为去除边界效应，模型两端各留设100 m煤柱。模拟工作面采高9.3 m，走向推进320 m，采用分步开挖，每步进尺10 m，共计32步。块体采用摩尔−库伦模型，节理采用库仑滑移模型。
工作面覆岩力学参数见表3。
表 3 工作面覆岩力学参数Table 3. Mechanical parameters of working face overburden岩性 体积模量/GPa 剪切模量/GPa 黏聚力/MPa 内摩擦角/(°) 抗拉强度/MPa 砂质泥岩 2.0 1.6 0.8 20 1.5 粉砂岩 13 3.33 8.0 33 3.46 细砂岩 14.4 5.85 4.83 55 3.86 泥岩 2 1.5 1.6 20 1.5 中砂岩 12.5 5.12 4.27 39 3.57 粗砂岩 7.8 2.91 6.89 35 1.5 煤 1.86 0.74 0.75 20 1.3 3.2 覆岩运移规律
3.2.1 覆岩运动场
工作面推进过程中，覆岩的运移情况如图10所示。
根据覆岩运移速度的大小分为3种区域，分别是加速下沉区（v≥0.2 m/s）、缓慢下沉区（0<v<0.2 m/s）、稳定区（v=0）。随着工作面推进，覆岩运动场由两区分布（加速下沉区、缓慢下沉区）逐渐演化为三区分布（加速下沉区、缓慢下沉区、稳定区），且稳定区逐渐向工作面推进方向与采空区上方扩张。
工作面推进120 m，工作面初次来压。亚关键层及其控制岩层垮落，垮落区上部覆岩下沉速度由采空区中部向两端递减。采空区中部覆岩最大下沉速度为1.01 m/s，而主关键层下沉速度为0.03 m/s，仅为垮落覆岩的3%。推进140 m，覆岩运动范围在水平与垂直方向均增大。此时采空区两端下沉速度较大，而采空区中部覆岩由于达到最大下沉量而不再下沉，下沉速度为0，成为稳定区。推进200 m，覆岩向上运动至地表，地表开始下沉，下沉速度为0.03 m/s。采空区左端与中部发展成稳定区，采空区右端和覆岩中部为加速下沉区，缓慢下沉区已发育至地表。推进290 m，地表下沉形成盆地。稳定区继续向垂直和水平方向扩展，而加速下沉区集中在靠近工作面的覆岩中下部。
为研究推进过程中覆岩下沉速度的变化规律，分别以工作面上方7 m、相距开切眼56 m作为水平与垂直测线，监测不同推进距离下覆岩的下沉速度。取4种典型情况下的测线数据，绘制成折线图，如图11所示。
![]()
图 11 不同推进距离的覆岩下沉速度Figure 11. Overburden subsidence speed with different advancing distance由图11a可知，初期（推进120 m）覆岩下沉速度较小，中期（推进140 m）覆岩下沉速度急剧增加，最高可达612.9 mm/s，后期（推进200～290 m）下沉速度大幅降低，在190 mm/s上下波动。说明同一层岩层下沉速度会经历急剧加速→大幅减速→稳定阶段。
由图11b可知，随着工作面推进，覆岩逐渐向上运动，且最大下沉速度逐渐减小，初期（推进120 m）最大下沉速度可达990.8 mm/s，后期（推进290 m）仅为31.3 mm/s，减小至初始速度的3.2%。
煤层开采后，工作面覆岩由下至上成组运动，发展到地面引起地表下沉，各岩层在时空耦合下非连续破断变形产生离层裂隙，接下来从离层变化与地表下沉的角度对覆岩运移规律进行阐述。
3.2.2 离层与地表下沉动态演化规律
1）离层演化规律。离层量的时空动态分布呈现为截然不同的两阶段，如图12所示。第一阶段：工作面初次来压（亚关键层初次断裂）→主关键层断裂前夕。此阶段内，离层量在煤层走向分布呈现为“高帽状”，工作面向前推进，高帽宽度逐渐增加，即离层范围随之扩大。一般情况下，最大离层量与推进距离正相关，初期（推进120 m）最大离层量为4 468 mm，末期最大离层量高达8 508 mm，为初期的1.9倍。第二阶段：主关键层断裂之后。此阶段内，离层量在煤层走向分布呈现为“双驼峰状”，且前驼峰（近工作面处驼峰）低于后驼峰（近开切眼处驼峰）。前后驼峰初始跨度相近，随着推进距离增加，后驼峰跨度基本不变，前驼峰跨度逐渐增加，离层范围向前扩张。前后驼峰分界处离层量逐渐降低，直至为零，即采空区中部覆岩首先被压实。
![]()
图 12 离层量时空动态分布Figure 12. Spatiotemporal dynamic distribution of stratification separation2）地表下沉规律。推进320 m时，记录模型顶部的垂直位移与水平位移，绘制成折线图，如图13所示（水平位移向右为正，垂直位移向上为正），代表地表的位移分布情况。
由图13可知，垂直位移类似1/2正弦图像，水平位移类似一个正弦图像。地表下沉量由采空区中部向两端递减，水平移动量在左右两端都经历增加至峰值→减小至0的过程，地表最大下沉处即为水平位移为0处。地表位移变化表现为中间深、边缘浅的盆地，主要下沉区与主要移近区均位于采空区中部，且主要移近区略大于主要下沉区。
洞悉地表位移分布情况后，分别以相距开切眼75、126 m两处典型位置作为地表水平位移与垂直位移测点，记录相关位移增加量（或减小量）以及下沉与移近速度（速度按推进1 m的地表下沉量或移近量计算，单位为mm/m），如图14所示，研究推进过程中地表的位移变化情况。
由图14可知，地表移近量与下沉量的变化趋势接近一致，但地表下沉量比移近量高出3～12倍，地表最大下沉量为7 909 mm，最大移近量为1 939 mm，即地表位移主要表现为下沉现象。初期地表下沉速度较大，而移近速度几乎从0开始增长，中后期地表下沉速度与移近速度变化趋势大致接近。地表初次下沉后，移近速度与下沉速度经历加速→减速→平稳3个阶段。
随着工作面向前推进，地表下沉值不断增大。地表最大下沉值及其增长速度（速度按推进1 m的地表最大下沉量计算，单位为mm/m）、最大下沉位置的动态演化过程如图15所示。
由图15可知，地表初次下沉最大值仅为962 mm，而后期下沉最大值高达7 909 mm，增长了7.2倍，且后期最大下沉值几乎不变。初次下沉后，地表最大下沉值经历匀加速增长→加速度减小的增长。最大下沉位置随着工作面推进不断向前运动，且一直处于采空区中部，后期最大下沉位置为一段距离，而非一点，即盆地中央。
3.3 地表下沉机理
地表下沉示意图如图16所示。亚关键层断裂后，随着工作面向前推进，覆岩由下至上逐层运动，主关键层与亚关键层之间岩层的不协调变形产生离层裂隙。主关键层断裂后，直至地表的上覆岩层共同下沉，地表开始下沉，离层裂隙在煤层走向方向呈“双驼峰状”。
工作面继续推进，在岩层自重与上覆岩层载荷作用下，采空区中部岩层由下至上逐渐压实，离层裂隙闭合，地表下沉范围与下沉值增大。离层裂隙的动态演化、破碎岩体的压实过程与地表下沉息息相关，有助于研究地表下沉规律。地表开始下沉时，任一点处的地表下沉值理论计算式可表达为
$$ u = M - d - \sum\limits_{i = 1}^n {({K_i} - 1)} {h_i} $$ $$ \Delta h = \sum\limits_{i = 1}^n {({K_i} - 1)} {h_i} $$ 式中：u为任一点处地表下沉值，m；M为煤层厚度，m；d为该点下方离层量之和，m；n为导水裂隙带间岩层总数，层；Ki、hi分别为各岩层的碎胀系数、厚度，m；$ \Delta h $为岩层破断后体积膨胀导致的空间上升高度，m。
直接顶破碎较为严重，碎胀系数较大，基本顶破断主要形成层间结构，破断块体较大，碎胀系数较小，所以$ \Delta h $以直接顶为主。地表下沉初期，离层量d与$ \Delta h $较大，故而地表下沉值较小。工作面向前推进，采空区中部岩层逐渐被压实，地表下沉值理论计算式可表达为：
$$ u = M - d - \left[ {\sum\limits_{i = 1}^m {({K_{{\mathrm{P}}i}} - 1)} {h_{1i}} + \sum\limits_{i = 1}^x {({K_i} - 1)} {h_{2i}}} \right] $$ $$ \Delta h = \left[ {\sum\limits_{i = 1}^m {({K_{{\mathrm{P}}i}} - 1)} {h_{1i}} + \sum\limits_{i = 1}^x {({K_i} - 1)} {h_{2i}}} \right] $$ $$ h = {h_{1i}} + {h_{2i}} $$ 式中：m为压实的岩层总数，层；KPi、h1i分别为压实各岩层的残余碎胀系数与厚度，m；x为未被压实的岩层总数，层；h2i为未压实的各岩层厚度，m；h为导水裂隙带的岩层总厚度，m。
离层量d由采空区中部向两端递增，且采空区中部覆岩被压实导致其$ \Delta h $值低于两端，所以地表越靠近采空区中部，下沉值越大，两端下沉值较小，呈现为盆地。岩层被压实导致裂隙闭合，相应的d、$ \Delta h $减小，地表下沉值逐渐增大，直至岩层被完全压实，离层量d减小为0，$ \Delta h $降为最低，地表下沉达到最大值，之后进入充分采动阶段。地表最大下沉值理论计算式可表达为
$$ {u_{\max }} = M - \sum\limits_{i = 1}^n {({K_{{\mathrm{P}}i}} - 1){h_i}} $$ 由此可知，煤层厚度、离层量、导水裂隙带内岩层的厚度、覆岩岩性及物理性质共同决定地表下沉量。而煤层厚度与离层量较大，起决定性作用，一般情况下，地表最大下沉量与煤层厚度正相关。在一定开采条件下，通过离层区注浆充填，使得充填体代替离层量d，且充填体高度不会随工作面推进而减小，便能有效控制地表下沉，实现煤矿绿色开采。
4. 覆岩断裂演化过程
工作面推进过程中，覆岩的断裂演化过程如图17所示。工作面推进过程中，覆岩断裂形态会经历3个阶段。第1阶段：亚关键层断裂，主关键层处于稳定阶段，工作面初次来压。采空区下部为矩形垮落区，上部为砌体梁结构，覆岩断裂形态为等腰梯形。工作面继续推进，垮落区高度增加，覆岩离层位置逐渐向上发育，且离层区朝水平与垂直方向扩张。第2阶段：主关键层断裂，覆岩断裂形态演化为双等腰梯形，表土层影响区为矩形。采空区两端纵向与水平裂隙发育较高，为裂隙富集区。此时，离层区纵向高度最大，矩形垮落区高度达到最大，之后由于覆岩压实作用导致矩形垮落区高度小幅度降低并趋于稳定。过渡期：覆岩断裂位置发育至最高，即导水裂隙带达到最大高度，双等腰梯形高度达到最大。表土层影响区由矩形演化为倒梯形，且地表影响范围逐渐扩大，地表影响角$ \alpha $逐渐降低。采空区中部由下至上为压实区和离层区，压实区逐渐向上发育而增大，离层逐渐闭合导致离层区逐渐减小；两端为裂隙富集区，且裂隙富集区达到最大。第3阶段：地表下沉达到最大值，进入充分采动阶段。此时，地表影响角降至最低，采空区中部覆岩全部被压实，压实区达到最大，离层区消失。工作面继续推进，覆岩断裂位置跟随工作面前进，裂隙富集区处于动态平衡状态，覆岩断裂双等腰梯形宽度增加，但高度不变。
5. 结 论
1）亚关键层和主关键层的断裂共同决定覆岩的下沉，覆岩下沉量由下至上逐层递减。导水裂隙带内岩层随工作面推进表现为台阶下沉，同层岩层下沉趋势沿走向表现为急剧下降→稳定（最大值）→快速上升→稳定（最小值）。
2）覆岩运动场由初次来压时的两区分布（加速下沉区、缓慢下沉区）逐渐演化为3区分布（加速下沉区、缓慢下沉区、稳定区），同一层岩层下沉速度会经历急剧加速→大幅减速→稳定阶段，覆岩逐渐向上运动，且最大下沉速度逐渐减小。
3）通过理论分析，得出了地表开始下沉时任意点下沉值的理论计算式与地表最大下沉值的理论计算式，且煤层厚度与离层量对地表最大下沉值起决定性作用。
4）通过相似模拟与数值模拟的试验结果，提出了覆岩断裂的动态演化过程。其演化过程分为3阶段和1个过渡期，覆岩断裂形态由单等腰梯形演化为双等腰梯形，表土层影响区由矩形演化为倒梯形。分析了覆岩中垮落区、离层区、压实区、裂隙富集区的动态演化过程。
-
![]()
图 1 工作面煤岩层综合柱状
Figure 1. Comprehensive histogram of coal and rock strata in working face
![]()
图 11 不同推进距离的覆岩下沉速度
Figure 11. Overburden subsidence speed with different advancing distance
![]()
图 12 离层量时空动态分布
Figure 12. Spatiotemporal dynamic distribution of stratification separation
表 1 液压支架支护工作阻力
Table 1 Working resistance of hydraulic support
分站编号 P/MPa (P+σ)/MPa 动载系数 5号 30.78 35.72 1.16 15号 31.44 36.25 1.15 25号 31.90 36.45 1.14 35号 31.67 36.39 1.15 45号 31.49 36.33 1.15 55号 32.02 37.56 1.17 65号 32.37 37.82 1.17 75号 32.11 36.92 1.15 85号 31.44 36.20 1.15 95号 31.39 36.22 1.15 105号 31.36 35.94 1.15 115号 30.87 35.34 1.14 125号 30.05 34.93 1.16 平均值 31.45 36.31 1.15 注：P为工作阻力；σ为均方差。 
下载: 导出CSV
表 2 岩层相似材料配比与关键层位置
Table 2 Proportion of similar materials in rock strata and location of key layers
层号 岩性 厚度/m 容重/（kN·m−3） 模拟层厚/mm 分层数 配比号 关键层位置 18 （配重）上覆岩层 402.00 — — — — — 17 粉砂岩 8.27 24.00 41.4 2 755 16 细砂岩 15.88 25.30 79.4 4 655 15 泥岩 12.56 24.20 62.8 3 855 14 中粒砂岩 20.26 24.50 101.3 5 655 13 粗砂岩 5.88 23.50 29.4 1 655 12 粉砂岩 18.46 24.00 92.3 3 755 11 中粒砂岩 7.61 24.50 38.1 2 655 10 细砂岩 11.10 25.30 55.5 2 655 9 泥岩 8.00 24.20 40.0 2 855 8 粗砂岩 55.70 23.50 278.5 9 655 主关键层 7 一层煤 1.60 14.50 8.0 1 855 6 粉砂岩 3.80 24.00 19.0 1 755 5 中粒砂岩 18.40 24.50 92.0 4 655 亚关键层 4 粉砂岩 2.40 24.00 12.0 1 755 3 二层煤 9.30 15.00 46.5 2 855 2 泥岩 3.00 24.20 15.0 1 855 1 粉砂岩 7.90 24.00 39.5 2 755 注：配比号755 表示细砂、石灰、石膏质量比为7∶0.5∶0.5；岩层总厚度210.12 m，模拟层总厚1 050.6 mm，总计分层数45层。
下载: 导出CSV
表 3 工作面覆岩力学参数
Table 3 Mechanical parameters of working face overburden
岩性 体积模量/GPa 剪切模量/GPa 黏聚力/MPa 内摩擦角/(°) 抗拉强度/MPa 砂质泥岩 2.0 1.6 0.8 20 1.5 粉砂岩 13 3.33 8.0 33 3.46 细砂岩 14.4 5.85 4.83 55 3.86 泥岩 2 1.5 1.6 20 1.5 中砂岩 12.5 5.12 4.27 39 3.57 粗砂岩 7.8 2.91 6.89 35 1.5 煤 1.86 0.74 0.75 20 1.3
下载: 导出CSV
-
[1] 王再岚,智颖飙,张东海,等. 我国煤炭资源禀赋与国际储量格局分析[J]. 中国人口·资源与环境,2010,20(S1):318−320. WANG Z L,ZHI Y B,ZHANG D H,et al. China’s coal resources endowments with the pattern of international reserves[J]. China Population,Resources and Environment,2010,20(S1):318−320.
[2] 李 霞,崔 涛. 我国煤炭资源可持续发展的保障分析[J]. 中国煤炭,2019,45(1):33−37. LI X,CUI T. Supportability analysis of sustainable development for coal resources in China[J]. China Coal,2019,45(1):33−37.
[3] 王志强,武 超,罗健侨,等. 特厚煤层巨厚顶板分层综采工作面区段煤柱失稳机理及控制[J]. 煤炭学报,2021,46(12):3756−3770. WANG Z Q,WU C,LUO J Q,et al. Instability mechanism and control of section coal pillar in fully mechanized mining face with super thick roof and extra thick seam[J]. Journal of China Coal Society,2021,46(12):3756−3770.
[4] 于斌,匡铁军,杨敬轩,等. 特厚煤层开采坚硬顶板覆岩结构及其演化特征分析[J]. 煤炭科学技术,2023,51(1):95−104. YU Bin,KUANG Tiejun,YANG Jingxuan,et al. Analysis of overburden structure and evolution characteristics of hard roof mining in extremely thick coal seam[J]. Coal Science and Technology,2023,51(1):95−104.
[5] 王国法,庞义辉. 特厚煤层大采高综采综放适应性评价和技术原理[J]. 煤炭学报,2018,43(1):33−42. WANG G F,PANG Y H. Full-mechanized coal mining and caving mining method evaluation and key technology for thick coal seam[J]. Journal of China Coal Society,2018,43(1):33−42.
[6] 于 斌. 大同矿区特厚煤层综放开采强矿压显现机理及顶板控制研究[D]. 徐州:中国矿业大学,2014:60−82. YU B. Study on strong pressure behavior mechanism and roof control of fully mechanized top coal caving in extra thickness seam in Datong Coal Mine [D]. Xuzhou:China University of Mining and Technology,2014:60−82.
[7] 薛吉胜,赵铁林,潘黎明. “高位-低位”厚硬岩层综放面特厚煤层矿压显现特征研究[J]. 煤炭技术,2021,4(7):55−59. XUE J S,ZHAO T L,PAN L M. Research on characteristics of underground pressure behavior of extra-thick coal seams in “high-low” thick and hard rock layers in fully mechanized caving working face[J]. Coal Technology,2021,4(7):55−59.
[8] 刘垚鑫,高明仕,贺永亮,等. 倾斜特厚煤层综放沿空掘巷围岩稳定性研究[J]. 中国矿业大学学报,2021,50(6):1051−1059. LIU Y X,GAO M S,HE Y L,et al. Study of control technology about gob-side entry driving with top-coal caving in inclined extra-thick coal seam[J]. Journal of China University of Mining and Technology,2021,50(6):1051−1059.
[9] 孔令海. 特厚煤层大空间综放采场覆岩运动及其来压规律研究[J]. 采矿与安全工程学报,2020,37(5):943−950. KONG L H. Overlying strata movement law and its strata pressure mechanism in fully mechanized top-coal caving workface with large space[J]. Journal of Mining and Safety Engineering,2020,37(5):943−950.
[10] 窦桂东,林建成,杜 鑫,等. 特厚煤层综采工作面矿压显现规律研究[J]. 煤炭工程,2019,51(8):84−88. DOU G D,LIN J C,DU X,et al. Study on mining pressure behavior of fully mechanized top-coal caving face in extra-thick coal seam[J]. Coal Engineering,2019,51(8):84−88.
[11] 王云广,郭文兵,白二虎,等. 高强度开采覆岩运移特征与机理研究[J]. 煤炭学报,2018,43(S1):28−35. WANG Y G,GUO W B,BAI E H,et al. Characteristics and mechanism of overlying strata movement due to high-intensity mining[J]. Journal of China Coal Society,2018,43(S1):28−35.
[12] 蒋金泉,王 普,武泉林,等. 上覆高位岩浆岩下离层空间的演化规律及其预测[J]. 岩土工程学报,2015,37(10):1769−1779. JIANG J Q,WANG P,WU Q L,et al. Evolution laws and prediction of separated stratum space under overlying high-position magmatic rocks[J]. Chinese Journal of Geotechnical Engineering,2015,37(10):1769−1779.
[13] 许满贵,魏 攀,李树刚,等. “三软”煤层综采工作面覆岩运移和裂隙演化规律实验研究[J]. 煤炭学报,2017,42(S1):122−127. XU M G,WEI P,LI S G,et al. Experimental study on overburden migration and fracture evolution law of “three soft” coal seam fully mechanized working-face[J]. Journal of China Coal Society,2017,42(S1):122−127.
[14] 殷 伟,张 强,韩晓乐,等. 混合综采工作面覆岩运移规律及空间结构特征分析[J]. 煤炭学报,2017,42(2):388−396. YIN W,ZHANG Q,HAN X L,et al. Overlying strata movement law and spatial structure analysis of fully mechanized mixed mining of backfilling and caving[J]. Journal of China Coal Society,2017,42(2):388−396.
[15] 赵毅鑫,令春伟,刘 斌,等. 浅埋超大采高工作面覆岩裂隙演化及能量耗散规律研究[J]. 采矿与安全工程学报,2021,38(1):9−18,30. ZHAO Y X,LING C W,LIU B,et al. Fracture evolution and energy dissipation of overlying strata in shallow-buried underground mining with ultra-high working face[J]. Journal of Mining and Safety Engineering,2021,38(1):9−18,30.
[16] 张广超,陶广哲,孟祥军,等. 巨厚松散层下软弱覆岩破坏规律[J]. 煤炭学报,2022,47(11):3998−4010. doi: 10.13225/j.cnki.jccs.2021.1994 ZHANG G C,TAO G Z,MENG X J,et al. Failure law of weak overburden stratum underlying extra-thick alluvium[J]. Journal of China Coal Society,2022,47(11):3998−4010. doi: 10.13225/j.cnki.jccs.2021.1994
[17] 杨 科,刘文杰,焦 彪,等. 深部厚硬顶板综放开采覆岩运移三维物理模拟试验研究[J]. 岩土工程学报,2021,43(1):85−93. YANG K,LIU W J,JIAO B,et al. Three-dimensional physical simulation of overburden migration in deep thick hard roof fully-mechanized caving mining[J]. Chinese Journal of Geotechnical Engineering,2021,43(1):85−93.
[18] 钱鸣高,缪协兴,许家林. 岩层控制中的关键层理论研究[J]. 煤炭学报,1996,21(3):2−7. QIAN M G,MIAO X X,XU J L. Theoretical study of key stratum in ground control[J]. Journal of China Coal Society,1996,21(3):2−7.
[19] 许家林,钱鸣高. 覆岩关键层位置的判别方法[J]. 中国矿业大学学报,2000,29(5):21−25. XU J L,QIAN M G. Method to distinguish key strata in overburden[J]. Journal of China University of Mining and Technology,2000,29(5):21−25.
[20] 缪协兴,陈荣华,浦 海,等. 采场覆岩厚关键层破断与冒落规律分析[J]. 岩石力学与工程学报,2005,24(8):1289−1295. MIAO X X,CHEN R H,PU H,et al. Analysis of breakage and collapse of thick key strata around coal face[J]. Chinese Journal of Rock Mechanics and Engineering,2005,24(8):1289−1295.
[21] 许家林,鞠金峰. 特大采高综采面关键层结构形态及其对矿压显现的影响[J]. 岩石力学与工程学报,2011,30(8):1547−1556. XU J L,JU J F. Structural morphology of key stratum and its influence on strata behaviors in fully-mechanized face with super-large mining height[J]. Chinese Journal of Rock Mechanics and Engineering,2011,30(8):1547−1556.
[22] 闫少宏,尹希文,许红杰,等. 大采高综采顶板短悬臂梁-铰接岩梁结构与支架工作阻力的确定[J]. 煤炭学报,2011,36(11):1816−1820. YAN S H,YIN X W,XU H J,et al. Roof structure of short cantilever-articulated rock beam and calculation of support resistance in full-mechanized face with large mining height[J]. Journal of China Coal Society,2011,36(11):1816−1820.
[23] 付玉平,宋选民,邢平伟,等. 浅埋厚煤层大采高工作面顶板岩层断裂演化规律的模拟研究[J]. 煤炭学报,2012,37(3):366−371. FU Y P,SONG X M,XING P W,et al. Study on simulation of caving and evolution law of roof strata of large mining height workface in shallow thick coal seam[J]. Journal of China Coal Society,2012,37(3):366−371.
[24] 柴 敬,汪志力,刘文岗,等. 采场上覆关键层运移的模拟实验检测[J]. 煤炭学报,2015,40(1):35−41. CHAI J,WANG Z L,LIU W G,et al. Monitoring movement laws of overlying key strata for coal mining in similar model[J]. Journal of China Coal Society,2015,40(1):35−41.
[25] 任艳芳. 浅埋深工作面覆岩“悬臂梁-铰接岩梁”结构的提出与验证[J]. 煤炭学报,2019,44(S1):1−8. REN Y F. Presentation and verification of “cantilever beam-articulated rock beam” composite structure in shallow buried working face[J]. Journal of China Coal Society,2019,44(S1):1−8.
-
期刊类型引用(2)
1. 刘旺海. 大采高综采工作面采动覆岩时空演化特征. 采矿技术. 2025(01): 23-28 . 
百度学术
2. 张乐强，王伟，刘晓彤，冯昌如，朱昊. 深部开采超宽工作面覆岩运移破断规律研究. 煤炭科技. 2024(04): 26-32 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 205
- HTML全文浏览量: 27
- PDF下载量: 120
- 被引次数: 7
