Research on surface subsidence law in high-intensity mining of shallow buried with thick coal seam
-
摘要:
以西部神东矿区大柳塔煤矿为工程背景，选取典型的浅埋深高强度开采工作面52307工作面为研究对象，设立地表岩移观测站，并结合RTK观测技术与三维激光扫描技术进行实时监测，基于关键层理论和CISPM综合地表沉陷预计模型软件，综合分析在浅埋深高强度开采条件影响下，地表移动变形特征、地表岩移角值参数变化以及地裂缝发育等问题。研究结果表明：在工作面从开切眼向前推进过程中，地表下沉值在开切眼位置较小，之后突然迅速增大，下沉曲线急剧变陡，当达到最大下沉值时，下沉速度开始变缓；地表移动变形主要集中在工作面中部，工作面四周下沉变形较小，地表移动影响范围较小；大柳塔矿区地表岩移角值参数偏大，基岩移动角和边界角分别达到了87.7°和84.1°。开采引起的地裂缝整体呈“C”字型主要分布在工作面中部；随着工作面的不断推进，地裂缝从开切眼位置开始继续向回风巷一侧缓慢延伸、发育并逐渐闭合，最终形成连续包围“漏斗”型形状；地裂缝总是滞后于工作面出现，地裂缝滞后距与工作面推进速度呈线性增大的关系。根据地表沉陷规律及地裂缝发育分析结果得出，在浅埋深厚煤层高强度开采条件下，由于工作面开采强度大、推进速度快、关键层结构单一且赋存层位低，引起顶板活动剧烈、易滑落失稳，从而造成地表下沉盆地四周收敛快、中部破坏严重，地裂缝发育密集。
Abstract:Taking Daliuta Coal Mine in Western Shendong mining area as the engineering background, this paper selected typical shallow buried high-intensity mining panel 52307 to establish the observation station for surface movement. Combined RTK technology and 3D laser scanning technology for real-time monitoring, based on key layer theory and CISPM comprehensive surface subsidence prediction model software, the characteristics of surface movement and deformation, surface movement angular parameters and surface cracks in high-intensity mining of shallow buried was studied. The results showed that with the advance of the panel from the set-up room, the surface subsidence was small at the set-up room position and then increased suddenly and rapidly. Simultaneously, the subsidence curve became steep sharply. When the maximum subsidence value was reached, the subsidence velocity tended to slow down. The surface movement deformation was mainly concentrated in the middle of the panel, while the subsidence deformation around the panel and the influence range of surface movement was small. The surface movement angular parameters were quite larger in the Daliuta coal mine area, while the bedrock movement angle and boundary angle reached 87.7 ° and 84.1 ° respectively. The surface cracks caused by mining were in an overall “C” shape, mainly distributed in the middle of the panel. Then, As the working face continued to advance, the ground fissures continue to slowly extend, develop, and gradually close from the set-up room to the tailgate side, eventually forming a continuous surrounding "funnel" shape. Ground fissures always lagged behind the working face position and the lag distance of surface cracks increases linearly with the mining speed. According to the analysis of surface subsidence law and surface cracks development, it was concluded that under the condition of high-intensity mining of shallow buried with thick coal seam, due to the high mining intensity of the working face, fast advancing speed, single key layer structure and low occurrence horizon, the roof activity was intense and easy to slide and lose stability. As a result, it was appeared resulting in rapid convergence around the surface movement basin, serious damage in the middle, and intensive development of ground fissures.
-
0. 引 言
煤炭是我国的主体能源，是保障经济社会实现长足高效发展的主干驱动力。近年来，由于对煤炭资源的高强度开采，导致东部地区可采储量面临枯竭，煤炭资源开发逐渐向西部转移[1-3]。西部某些矿区上覆巨厚高承压含水层、覆岩原生裂隙发育等采矿地质条件较为突出，使覆岩导水裂隙网络发育及构造特征更加复杂、涌水形态愈加难辨，进一步加大了采场涌水量，严重影响工作面正常回采进程，且加剧了地下水资源的损耗与流失，使该地区地下水生态环境更加脆弱不堪[4-6]。
为深入了解采动诱发覆岩运移发育导水通道网络特征、突水致灾机理等，相关领域内的专家学者对此展开了研究。在覆岩运移发育导水通道网络特征方面，白贤栖等[7]发现采动诱发下高位巨厚岩层发生“O-X”型破断是诱发矿震的主要原因之一；王炯等[8]的研究表明，红外热辐射演化特征反映了层间相互作用关系的变化，切顶卸压自成巷覆岩运移呈现规律性渐进演变；关守安等[9]利用数值模拟软件分析得出，粗榆金矿矿体岩层沉降显著，但最终倾斜量小于临界值；吕兆海等[10]通过开展相似材料试验发现近距离煤层开采时，下位煤层开采会对上位煤层顶板产生二次扰动；殷伟等[11]研究得出，混合综采工作面顶板覆岩空间运移结构呈现非对称性。在煤矿顶板突水致灾机理方面，ZENG等[12]通过建立固−流耦合模型得出了顶板水砂混合物运移特征及其与覆岩水压和渗透特性之间的动态联动关系；侯恩科等[13]采用RF算法和GIS软件预测了4-2煤掘进巷道顶板突水危险区的位置；杨涛等[14]研究表明，采动使隔水土层孔隙水压聚集直至突破隔水土层，诱使工作面发生突水；许延春等[15]发现，基岩风化带隔水性能较好，有利于导水裂隙的闭合，适合作为防水煤柱保护层；李杨杨等[16]分析得出隧道开挖过断层时，当断层被完全揭露后断层围岩完全破裂发生涌水溃砂。
现有研究成果从理论、试验等角度对煤矿覆岩运移破断规律和顶板涌水溃砂致灾机理进行了大量研究，并取得了较为丰硕的成果。但是，当前研究多基于建立完整地层构造模型展开，鲜有考虑原生裂隙等导水结构对覆岩运移破断发育导水通道网络机制及采场涌（淋）水特征影响的研究报道。
在工程实际中，覆岩内均存在几何尺度及形态结构不一的原生裂隙等导水结构[17-18]，这些导水结构会对覆岩运动及采场涌（淋）水形态产生重要影响。因此，笔者结合煤矿实际开采条件开展固−流耦合相似模拟试验，对原生裂隙诱导下覆岩运移破坏诱发导水裂隙网络形成机制、采场涌（淋）水形态特征等进行研究。研究成果可为具备类似地质条件的工程现场开展防治顶板水事故提供理论参考。
1. 研究区域概况
陕西彬长矿区某矿行政区划隶属于陕西省咸阳市，可采煤层为4上、4−1、4煤，埋深为600～800 m。矿区构造位于鄂尔多斯盆地西南边缘，总体呈现北西向较缓倾斜的单斜构造，短轴背斜与宽缓向斜组成褶皱构造，轴向近NEE-SWW，地层倾角一般小于10°。该矿全区上覆巨厚、高承压洛河组含水层，厚7.35～464.84 m，首采区均厚400.35 m，最大水头压力达5.27 MPa，单位涌水量0.218～1.130 L/(s·m)，渗透系数0.05471～0.32465 m/d，属中等～强富水性含水层。同时，含水层底板与可采煤层顶板间距较小，平均仅84.02 m。隔水层为安定组泥岩，主要由泥岩和砂质泥岩组成，厚0～94.31 m，平均47.15 m，全区大部分地区厚度小于20 m，且有零厚度地段，整体隔水性能不足。同时，该矿全区覆岩中原生裂隙发育，发育类型以高角度构造裂隙为主，成岩裂隙次之。在空间结构特征方面，原生裂隙在平面上延展长度较小，相互叠置，呈雁行式排列，在纵向上呈短−中短蠕虫状垂向延伸，整体展布特征复杂，难以进行精细化探查。长观孔水位监测数据显示，该矿最大涌水量达到5365 m3/h，最小亦达到2886 m3/h，涌水量受工作面推进距离、采放高度、采空区面积等采矿参数影响显著。采空区及工作面顶板存在大范围淋水，且工作面煤壁侧存在超前出水现象，严重影响该矿井正常生产作业。
2. 考虑原生裂隙的覆岩固−流耦合相似模拟试验
2.1 试验系统
试验在山东科技大学自研覆岩涌水溃砂相似模拟物理试验系统上开展[19]，采用DM-YB1820型数据采集器和DMKY型水压计。试验系统（图1）主要由试验平台、供压系统、控制系统3个子系统组成，其中试验平台包括试验舱、压头、承压水仓3个子系统。试验舱几何尺寸为1200 mm×700 mm×400 mm，前置亚克力板，其余各面采用高强度钢材制作；压头和承压水仓为应力和水压的输出端。供压系统通过液压油为试验舱提供载荷，是应力和水压的输入端。控制系统包括位移应力和水压水量伺服控制2个单元，可实现对试验过程中应力、位移、水压、水量等的伺服控制。
![]()
图 1 覆岩涌水溃砂相似模拟物理实验系统Figure 1. Similar simulation physical experiment system of water and sand inrush in overlying strata2.2 相似材料配比优化
结合研究区域内覆岩岩性展布特征及本次试验设备几何尺寸，并基于相似比有关理论公式[20]，综合确定试验几何、时间、容重、强度、渗透系数相似比分别为1∶200、1∶14.14、1∶1.5、1∶300、1∶11.5。矿区覆岩物理力学性质见表1。
表 1 矿区覆岩物理力学性质Table 1. Physical and mechanical properties of overlying strata in mining area岩性 密度
/(kg·m−3）抗压强度
/MPa抗拉强度
/MPa弹性模量
/MPa泊松比 孔隙率
/%渗透系数
/（m·d−1）普氏系数 泥岩 2470 30.20 0.86 8.38 0.25 11.43 0.022 2.88 砂质泥岩 2598 25.90 1.98 12.27 0.20 13.05 2.04×10−4 3.09 粉砂岩 2526 65.60 2.95 13.16 0.18 5.82 0.053 6.69 细粒砂岩 2291 48.40 2.97 13.95 0.20 13.24 0.079 4.66 中粒砂岩 2080 31.80 1.48 9.36 0.24 21.54 0.806 3.42 粗粒砂岩 2276 41.30 1.21 10.21 0.20 15.43 1.064 4.38 试验采用河砂作为骨料，胶结材料选取凡士林、石蜡、碳酸钙，层间铺设云母粉进行分层。为保证相似材料配比的准确性，配制50 mm×100 mm和50 mm×25 mm的圆柱形标准相似材料，试件分别开展单轴压缩和巴西劈裂试验对其力学参数进行检测，进而对配比进行优化使之更符合需求。每种岩性、尺寸各配制3个试件。
配制相似材料试件采用山东科技大学自研恒定加压数显圆柱形相似材料试件制备装置，整体结构如图2所示。该装置具有结构简单，易于操作，脱模便捷等优点，模具直径为50 mm，总长120 mm，有20 mm的可控变长范围。装置底部设置有高精度压感器，可通过数显器对试件的加载力进行实时展示，实现对压力的动态监测调节。压感器最大量程为500 kN，数显器读数精度为0.01 kN。
![]()
图 2 恒定加压数显圆柱形相似材料试件制备装置Figure 2. Constant pressure digital display cylindrical similar material specimen preparation device通过恒定加压数显圆柱形相似材料试件制备装置配制的部分相似材料标准圆柱形试件如图3所示。
利用山东科技大学AG-X250岛津电子万能试验机对相似材料试件进行强度测试试验。根据强度相似比计算各岩层单轴抗压强度、单轴抗拉强度及弹性模量相似值，与相似材料试件强度测试结果进行对比分析，如图4所示。
![]()
图 4 相似材料试件强度试验值和相似值对比Figure 4. Comparison of strength test values and similar values of similar material specimens由图4可见，各岩层强度试验值曲线与理论值曲线大致重合，表明该配比基本符合本次固−流耦合相似材料模拟试验需求。为深度优化配比，基于强度测试结果，调整相似材料试件配比，再开展力学性能测试试验，直至试件强度的试验测试值基本契合相似值。同时，采用相关设备对所配制相似材料试件的孔隙率、渗透系数等进行检测，将最终测定结果与相似值对比，发现基本符合对实际水文地质工况的模拟需求。限于篇幅，此处对孔隙率和渗透系数等参数的测定过程不再赘述。采用最终配比配制的相似材料试件的物理力学性质见表2。
表 2 配比优化后相似材料试件物理力学性质Table 2. Physical and mechanical properties of similar material specimens after ratio optimization岩层 单轴抗压
强度/kN单轴抗
拉
强度/kN弹性模量/
GPa孔隙率/% 渗透系数/
（10−3m·d−1）泥岩 0.196 0.005 0.027 0.903 1.86 砂质泥岩 0.171 0.012 0.040 1.211 0.0172 粉砂岩 0.427 0.020 0.042 0.443 5.06 细粒砂岩 0.315 0.021 0.048 1.071 6.71 中粒砂岩 0.210 0.009 0.029 1.706 67.3 粗粒砂岩 0.271 0.009 0.031 1.286 93.6 2.3 模型的建立及铺设开挖
2.3.1 构建地层模型
试验设计模型高度为70 cm，工作面倾向长度为40 cm，走向长度为100 cm（不含保护煤柱宽度）。在距离煤层底板上14 cm的砂质泥岩和粗粒砂岩之间、26.7 cm的细粒砂岩和粗粒砂岩之间各布设4个水压计（F1—F8），2个水压计水平间距为24 cm。在模型顶部进行预注水处理，模拟地层初始含水状态。
开挖前，在模型顶部施加补偿载荷和水压，补偿载荷计算式如下：
$$ F = {C_\sigma } {\gamma _{\rm{p}}} \left( {H - {H_{\rm{m}}}} \right) $$ (1) 式中：
$F$ 为单位补偿载荷，${\rm{N}}/{{\rm{m}}^3}$ ；${C_\sigma }$ 为应力相似比，1/300；${\gamma _{\rm{p}}}$ 为地层平均容重，${\rm{g}}/{\rm{m}}^3$ ，取25 kN/m3；$ H $ 为煤层顶板至地表的距离，${\rm{m}}$ ；$H_{{\rm{m}}}$ 为相似材料模型中煤层顶板至模型顶部边界对应的实际地层距离，${\rm{m}}$ 。将相关数据代入式（1）计算得施加的补偿载荷为0.149 MPa。根据强度相似比计算得施加水压为0.0175 MPa。建立的地层模型如图5所示。
根据该矿工作面主测线剖面图原生裂隙发育特征及试验实际需求，综合确定预制裂隙参数。设计单位裂隙长度为3 cm，宽0.15 cm。以裂隙顶部尖端为原点，水平方向为基准线，顺时针旋转为正，确定裂隙倾角。为便于表述，按照原生裂隙位置对其进行编号，自上往下划分为顶、中、底3个带，自左至右对裂隙进行顺序编号。具体编号及参数见表3。
表 3 原生裂隙编号及参数Table 3. Primary fracture number and parameters裂隙
分带裂隙
编号单位
数量长度/cm 倾角/（°） 宽度/cm 顶分带 1-1 1 3 90 0.15 1-2 1 3 60 0.15 1-3 1 3 40 0.15 中分带 1-1 1 3 140 0.15 1-2 1 3 120 0.15 底分带 1-1 2 3+3 90+60 0.15 1-2 2 3+3 90+20 0.15 1-3 2 3+3 90+140 0.15 注：“单位数量”指长度为3 cm的预制裂隙数量。 赋参后各组原生裂隙几何形态及在地层模型中的具体位置如图6所示。
2.3.2 模型铺设及开挖
以优化后的相似材料配比配制岩层，依次铺设，层间铺设厚约0.15 cm的云母粉，在既定位置处布设水压计，并按照既定尺寸在对应位置处放置厚0.15 cm的中空板模拟原生裂隙。模型铺设完毕后（图7），在顶部施加设定的地应力补偿载荷和水压，而后开始开挖。
模型两边界各留设10 cm保护煤柱，自模型左侧开始开挖，每次开挖10 cm，共计开挖10次。模型底部均匀铺设尺寸为40 cm×10 cm×3 cm、外侧表面有螺孔的特制钢板，用以模拟煤层开挖。开挖时，将抽采装置前端的螺纹钉拧入需抽出钢板的螺孔内，旋转控制盘将螺纹钉向试验人员一侧拉近，进而带动钢板运动，直至将其抽出，实现对煤层回采的模拟。开挖期间，定时记录覆岩运移特征及工作面涌（淋）水状态。
3. 原生裂隙对覆岩导水通道网络发育影响
结合固−流耦合相似材料模拟试验结果，分阶段分析覆岩原生裂隙对导水通道网络发育规律的诱导效应。
3.1 覆岩导水通道网络发育扩展规律分析
综合工作面不同推进距离内覆岩运移共性特征，将其划分为3个阶段：直接顶常规性沉降冒落阶段、“Z”形离层空间发育阶段、裂隙空间压密闭合阶段。
1）直接顶常规性沉降冒落阶段。该阶段包括工作面推进至100 m前，覆岩发育特征如图8所示。工作面初次开挖后，覆岩未出现显著运移，仅存在轻微的弯曲下沉（图8a）；工作面开挖至60 m时，直接顶受采动影响发生整体性破断冒落，堆矸破碎程度低（图8b）。当工作面回采至100 m时，直接顶底部岩层发生破断冒落，顶板冒落岩层之上层间发生轻度分离（图8c）。
2）“Z”形离层空间发育阶段。图9为“Z”形离层空间发育过程。当工作面推进至120 m时，导水裂隙带高度发育至首层关键层下，关键层与其下位岩层层间产生不同步沉降，形成了显著的离层空间，且采动裂隙已贯通离层空间（图9a）。
随工作面继续推进，离层下位岩层在断裂线左侧发生整体性破断（图9b），断裂后的两段岩层发生不同步运移：左侧段完全冒落，右侧段在岩体支撑下形成悬臂梁，在悬臂梁的支撑作用下，其下部形成一个三角形空间。分析原因，受原生裂隙影响，离层下位岩层自身横向上结构整体性出现突变，使岩层承载能力下降，并最终受采动影响发生超前破断。原生裂隙的存在使受扰动岩层破断位置向采空区中心侧发生偏移。
图10为覆岩导水裂隙网络展布特征。工作面推进至160 m时，导水裂隙带高度已发育至“底分带”原生裂隙所处层位之上。此时，在更高层位关键层支撑作用下，新的离层开始发育，低层位离层被压实闭合。在原生裂隙影响下，该阶段发育的离层下位岩层再次发生超前破断（图9c），断裂点右侧段形成悬臂梁，其下方形成三角形空间，断裂点左侧段发生常规冒落，形成残缺离层空间。由于该岩层在超前断裂点处发生错动破断，使左侧残缺离层空间和右侧悬臂梁下的三角形空间连通，形成一个“Z”形离层空间2号（为便于表述，按照“Z”形离层空间发育先后顺序进行编号，此处为2号）。图9b中破断点左侧离层空间的上位岩层受原生裂隙影响发生破断，破断点如图10中A点所示，进而发育形成“Z”形离层空间1号。
![]()
图 10 覆岩导水裂隙网络展布特征Figure 10. Distribution characteristics of water-conducting fracture network in overlying strata同时在该阶段内，覆岩导水通道逐渐发育完善，2个“Z”形离层空间与原生裂隙、采动裂隙之间相互贯通，发育成同时具备储水和导水功能的裂隙网络，如图10中红色箭头所示区域。
3）裂隙空间压密闭合阶段。图11为覆岩裂隙空间压密闭合过程。当工作面推进至180 m时，低层位离层闭合，更高层位离层形成，且在第二阶段发育的“Z”形离层空间2号及由其形成的导水通道亦开始闭合，如图11a所示。同时，工作面的开挖波及至原生裂隙底1—3处，因此，受该原生裂隙影响，在工作面上方发育形成“Z”形离层空间3号。
![]()
图 11 覆岩裂隙空间压密闭合过程Figure 11. Sealing process of overlying rock fracture space pressure工作面回采完毕后，覆岩中仅工作面处和“Z”形离层空间1号仍存在部分可见空间未被压实，其他离层、裂隙等基本被压实闭合，新生导水通道发育基本停止。
对比1—3号“Z”形离层空间可见：不同倾角原生裂隙使岩层形成不同倾向的软弱结构面，进而使岩层破断面倾向产生差异。但断面倾向对“Z”形离层空间最终形态、体积等无显著影响。因此，原生裂隙倾角仅对岩层破断倾向产生影响，但对“Z”形离层空间及导水裂隙网络构造特征影响不显著。
由图11b可知，当导水裂隙带高度波及至中分带原生裂隙所处层位时，受原生裂隙影响岩层发生破裂，形成宏观可见的裂隙。但由于岩石具有碎胀性，低层位已破断碎裂岩层体积增大，使高层位岩层可运移空间逐渐减小，此时岩层基本呈现同步运移状态。因此，虽然中分带原生裂隙处发育有宏观可见裂隙，但由于岩层无足够的空间发生运移破断，所以中分带原生裂隙附近无“Z”形离层空间发育。回采结束后，顶分带原生裂隙处未见有宏观裂隙发育。由此可见，在纵向上，随原生裂隙层位升高，其对覆岩运移破断特征的影响程度减小。
可见，原生裂隙的存在使其影响范围内的岩层形成软弱结构面，整体结构性受到破坏。当受采动影响时，岩层会沿软弱结构面处发生部分或整体破断，从而发育形成结构不稳定的“Z”形离层空间。在覆岩运移作用下，原生裂隙与采动裂隙及“Z”形离层空间等互馈—贯通形成完整的导水通道网络，加大了工作面及采空区涌（淋）水量及面积。
3.2 “Z”形离层空间结构形态分析
为明确工作面回采过程中受原生裂隙影响发育的“Z”形离层空间结构形态及其与常规离层空间的差异，对二者进行对比分析（图12）。
从空间形态、体积及稳定性3方面对常规型离层空间和“Z”形离层空间进行对比分析：
1）在空间形态方面，常规型离层空间呈现倒三角形（实际情况下，下位岩层近似圆弧形，上位岩层亦存在一定弧度），而“Z”形离层空间呈现字母“Z”的形状。
2）在空间体积方面，“Z”形离层空间是常规型离层下位岩层存在部分未破断冒落而形成，因此理论上而言，“Z”形离层空间体积与对应条件下的常规型离层空间体积近似相同。
3）在空间稳定性方面，常规型离层空间最多由2种岩层构成，而“Z”形离层空间则最多由3种不同岩性的岩层组成，且由于悬臂梁的发育，导致下位岩层（隔水层）厚度相对较薄。因此，“Z”形离层空间整体稳定性更差，赋存时间可能相较于常规型离层空间更短。
综上分析可知，“Z”形离层空间结构形态相对更加复杂，稳定性差，更易受采动影响发生失稳破坏，加大采场涌（淋）水量及面积。
4. 覆岩涌（淋）水时空信息特征
4.1 覆岩涌（淋）水形态规律
根据开挖过程中工作面及采空区涌（淋）水形态将其划分为4个阶段：滴落态出水初期阶段、滴流态出水增加阶段、流动态出水稳定阶段、多元态出水回落阶段。图13为工作面推进过程覆岩涌（淋）水形态。
![]()
图 13 工作面推进过程覆岩涌（淋）水形态Figure 13. Form of water gushing (leaching) in overlying strata during the advancing process of working face1）滴落态出水初期阶段。工作面首次开挖后，顶板出现显著出水点，出水呈滴落状态，滴速均匀，出水点数量较少。与顶板出水状态相反，煤壁上存在大面积、高密度挂水现象，出水速度快，但出水点水量较小，如图13a所示。
2）滴流态出水增加阶段。当工作面开挖至60 m时，工作面及采空区顶板出现局部聚集型出水集中区（图13b），出水量增大，出水速度显著提升，出水状态为滴流。煤壁上出水点水量增大，同时，工作面及采空区堆矸上形成积水区，且积水已在出水口（图13b中红色圆圈处）形成稳定的水流向试验台水槽内流动。
3）流动态出水稳定阶段。如图13c所示，当工作面推进至100 m时，采空区堆矸内形成出水口，存在稳定流量的水体持续向外流出，说明覆岩涌水量随导水通道发育程度的提升而进一步加大。当工作面推进至160 m时，采空区涌水存在间歇性增减特征，但间隔时间较短。该现象表明，在本阶段内，覆岩中离层等空间已具备一定储水能力，且覆岩含水层、离层空间及采空区（或工作面）之间已形成贯通的涌水路径。
4）多元态出水回落阶段。工作面开挖结束后（图13d），采空区及工作面涌水量出现减小趋势，采空区侧呈现多点滴流状态，原生裂隙位置处存在宏观可见水流。工作面侧顶板、煤（岩）壁上及采空区堆矸面向工作面侧均存在出水点，滴速显著大于第一开挖阶段。分析原因，覆岩经过长时间破断运移，已接近充分采动状态，高位岩层的弯曲沉降使低位岩层内的裂隙、离层空间等发生闭合，使其导（储）水性能降低。因此，采空区及工作面涌水量出现下降态势。
综上所述，工作面及采空区涌（淋）水量与覆岩运移过程形成联动响应特征，涌（淋）水形态、水量等的演化均显著受覆岩运移影响。工作面推采前期，覆岩运动破断程度越大，采场涌（淋）水量越大，而原生裂隙的存在加速了覆岩破断冒落进程，使导水通道网络发育更快，分布更加密集。因此，原生裂隙是诱发采场涌（淋）水量增大及出水形态发生变化的关键因素之一。
4.2 覆岩水压演化规律
图14为模型中F1—F4水压计监测数据演化规律，据此分析工作面推进过程中覆岩水压变化特征。
当工作面初次推采过后，F1处水压值出现下降，降幅约0.0035 MPa，表明初次开挖后采场顶板覆岩弯曲沉降范围波及至F1所处层位岩层。随工作面继续推进，F2—F4处水压递次出现拐点，产生下降趋势，其中F4处直至工作面推进到100 m左右时才开始下降。这一现象说明，随工作面持续推进，采动影响范围逐渐向前推移，水压监测点处的覆岩陆续受采场扰动发生弯曲沉降甚至冒落，从而导致水压降低。
工作面推进至60 m时，F1处水压值降至最低，随后工作面推进至80 m左右时，F2处出现最低值区间，F3、F4处水压亦分别在160 m和180 m处降至最低。可见，当监测点所处层位岩层发生破断冒落时，监测水压会降至最小值，且各监测点处降低趋势拐点的出现位置随监测点与开切眼距离的增大而增大。
F1—F4处水压增大趋势拐点的出现位置随监测点与开切眼距离的增大而延后，且F4处最终水压值最小，F1次之，F2和F3处相对最大。该现象表明，随工作面推进，与开切眼较近处的监测点较早受采动影响发生运移，使水压降低，同时，该监测点处亦较早受上覆岩层破断冒落后被压实，导致水压出现回升趋势。通过各监测点最终水压值的区别可知，F2、F3处最终被上覆岩体压实闭合程度高，而两侧的F1、F4处上覆岩体在开切眼侧和工作面煤壁侧未开挖煤（岩）体的支承作用下未发生充分沉降，因此受压实程度较低。
结合前述分析可知，在开挖全过程中，覆岩破断运移程度、采场涌（淋）水量和覆岩水压之间存在联动响应关系。随覆岩运移程度增大，导水裂隙网络进一步发育扩展，进而诱发采场涌（淋）水量增大，监测水压降低。且在“Z”形离层空间发育阶段，覆岩导水裂隙网络分布最为密集，采场涌（淋）水量亦达到最大值，部分点位水压监测值出现最低点。因此可以说明，原生裂隙的存在促进了导水裂隙网络的发育和扩展，进而使采场涌（淋）水形态、水量等发生改变。
5. 结 论
1）将覆岩运移过程划分为直接顶常规性沉降冒落阶段、“Z”形离层空间发育阶段、裂隙空间压密闭合阶段3个阶段，受原生裂隙诱导作用使岩层发生超前破断是“Z”形离层空间形成的主要原因。
2）相较于常规型离层空间，“Z”型离层空间在空间结构形态上呈现区别于常规倒三角形的“Z”形；在空间体积上与常规型无显著差别；在稳定性上，组成“Z”型离层空间的岩性种类多，结构复杂，隔水层薄，稳定性相对更差，更易受采动影响发生失稳诱发采场大量涌水。
3）将采场涌（淋）水过程划分为滴落态出水初期阶段、滴流态出水增加阶段、流动态出水稳定阶段、多元态出水回落阶段4个阶段。整体而言，采场涌（淋）水量呈现先增大后小幅下降的趋势；原生裂隙的存在显著促进了覆岩中水体运移进程及采场涌（淋）水量，对水运移路径的迁移亦产生扰动。
4）随工作面推进，采场顶板区域内水压呈现先减小后增大的趋势，监测水压曲线拐点（包括下降和上升）出现位置均随监测点与开切眼间距离的增大而增延后；各监测点处水压最小值受覆岩采动影响自左至右递次出现；回采结束后，采空区中心区域附近受覆岩压实程度高，水压回升幅度较大，采空区边界处覆岩沉降程度较低，水压回升幅度较小。
-
![]()
图 4 走向方向上A测线地表下沉曲线
Figure 4. Surface subsidence curve of survey line A in strike direction
![]()
图 5 走向方向上主断面地表岩移角量参数
1—裂缝角；2—移动角；3—边界角；4—松散层移动角；5—充分采动角
Figure 5. Surface movement angular parameters on major section along the mining direction
![]()
图 6 三维激光扫描地表最终下沉全盆地云图
1—地裂缝密集区域; 2、3—沟谷地形区域
Figure 6. Cloud image of the final surface subsidence basin with 3D laser scanning technology
![]()
图 7 大柳塔煤矿52307工作面预计参数回归值
Figure 7. Regression of predicted parameters of panel 52307 in Daliuta Coal Mine
![]()
图 9 大柳塔52307工作面地裂缝空间分布
Figure 9. Spatial distribution of ground fissures of panel 52307 in Daliuta Coal Mine
![]()
图 10 地裂缝与工作面推进位置相对位置
Figure 10. Relative location of surface cracks and the working face advancing position
![]()
图 11 不同类型地裂缝滞后距与工作面推进速度之间的关系
Figure 11. Relationship between lag distance of different kind of surface cracks and mining speed
表 1 地表移动观测站设计参数
Table 1 Design parameters of observation station for surface movement
参数 走向观测A线 倾向观测B线 合计 测点号 A1～A18 B1～B20 测线实际长度/m 358 380 738 测点间距/m 20 20 测点数量/个 18 20 38 控制点数量/个 2 0 2 
下载: 导出CSV
表 2 大柳塔52307工作面地表不同岩移角值参数
Table 2 Different angular value parameters of surface movement of panel No.52307 in Daliuta
方向 边界角/(°) 移动角/(°) 充分采动角/(°) 走向 81 85 60 倾向 71 79 52
下载: 导出CSV
表 3 地表移动变形预计参数
Table 3 Predicted parameters of surface movement and deformation
下沉系数q 主要影响角tan β 拐点偏距s/m 主要影响角θ/(°) 水平移动系数b 主要影响半径r/m 0.53 3.69 61.2 74.8° 0.21 51.4
下载: 导出CSV
表 4 地表移动变形最大值对比
Table 4 Results comparison of the maximum values of surface movement and deformation
方法 最大下沉值
${W_{\max } }/ {{\rm{mm}}}$最大倾斜值
${i_{\max } }$/(mm·m−1)最大水平移动值
${U_{\max } }/ {{\rm{mm}}}$最大水平变形值
${\varepsilon _{\max }}$/(mm·m−1)地表沉陷规律 3 560.0 77.0 750.0 32.4 三维激光扫描 3 592.1 78.3 776.1 33.6 地表沉陷预计 3 556.0 76.5 729.4 31.5
下载: 导出CSV
表 5 不同类型地裂缝滞后距与开采速度统计结果
Table 5 Statistics of different surface cracks lag distance and mining speed
裂缝
类型监测日期 裂缝滞后距dL/m 工作面推进速度v/(m·d−1) 拉伸型 2016-03-11 18.4 8.9 2016-03-12 7.2 4.8 2016-03-14 11.6 4.6 2016-03-18 18.3 7.5 2016-03-19 7.6 4.4 2016-03-21 18.9 8.1 2016-03-23 21.1 10 2016-03-24 17 6.7 2016-03-29 16.1 7 2016-03-30 15 6.1 2016-05-12 13.8 7.6 2016-05-13 25.1 11.6 2016-05-14 20 9.9 2016-05-15 22.6 11.2 2016-05-16 15 6.65 2016-05-17 17 8.85 2016-05-18 22 9.8 2016-05-19 23 8.85 2016-05-20 33.6 13.5 2016-05-21 29 12.7 2016-05-22 22 11.5 台阶型 2016-03-11 27.8 8.9 2016-03-14 10.2 4.6 2016-03-23 38 10 2016-03-24 13.3 6.7 2016-03-26 15.6 7
下载: 导出CSV
-
[1] 2020 煤炭行业发展年度报告[R]. 北京: 中国煤炭工业协会, 2021. [2] 王伟东,李少杰,韩九曦. 世界主要煤炭资源国煤炭供需形势分析及行业发展展望[J]. 中国矿业,2015,24(2):5−9. doi: 10.3969/j.issn.1004-4051.2015.02.004 WANG Weidong,LI Shaojie,HAN Jiuxi. Analysis of the main global coal resource countries’ supply-demand structural trend and coal industry outlook[J]. China Mining Magazine,2015,24(2):5−9. doi: 10.3969/j.issn.1004-4051.2015.02.004
[3] 杨正凯,张金虎,张 宁. 特厚软煤大采高综放工作面成套装备技术研究[J]. 煤炭工程,2020,52(8):123−126. YANG Zhengkai,ZHANG Jinhu,ZAHNG Ning. Key technology in complete set of equipment for high cutting fully mechanized caving face in extra-thick soft coal[J]. Coal Engineering,2020,52(8):123−126.
[4] 黄庆享,赵萌烨,黄克军. 浅埋煤层群开采顶板双关键层结构及支护阻力研究[J]. 中国矿业大学学报,2019,48(1):71−77,86. doi: 10.13247/j.cnki.jcumt.000968 HUANG Qingxiang,ZHAO Mengye,HUANG Kejun. Study of roof double key strata structure and support resistance of shallow coal seams group mining[J]. Journal of China University of Mining and Technology,2019,48(1):71−77,86. doi: 10.13247/j.cnki.jcumt.000968
[5] 杨俊哲,刘前进. 8.8 m超大采高工作面矿压显现规律实测及机理分析[J]. 煤炭科学技术,2020,48(1):69−74. YANG Junzhe,LIU Qianjin. Analysis and measured of strata behavior law and mechanism of 8.8 m ultra- high mining height working face[J]. Coal Science and Technology,2020,48(1):69−74.
[6] 高 超,徐乃忠,何标庆,等. 关键层对特厚煤层综放开采地表沉陷规律的影响研究[J]. 煤炭科学技术,2019,47(9):229−234. doi: 10.13199/j.cnki.cst.2019.09.029 GAO Chao,XU Naizhong,HE Biaoqing,et al. Study on influence of key strata on surface subsidence law of fully-mechanized caving mining in extra-thick coal seam[J]. Coal Science and Technology,2019,47(9):229−234. doi: 10.13199/j.cnki.cst.2019.09.029
[7] AN Yanpei, ZHANG Nong, ZHAO Yiming, et al. Field and numerical investigation on roof failure and fracture control of thick coal seam roadway[J]. Engineering Failure Analysis, 2021(1):105594.
[8] 王 伟. 浅埋深厚煤层开采地表岩移观测分析[J]. 江西煤炭科技,2019(3):57−59,62. doi: 10.3969/j.issn.1006-2572.2019.03.020 WANG Wei. Observation and analysis of surface rock movement in mining of shallow depth coal seam[J]. Jiangxi Coal Science and Technology,2019(3):57−59,62. doi: 10.3969/j.issn.1006-2572.2019.03.020
[9] 李 磊. 榆林矿区浅埋深厚土层薄基岩煤层开采覆岩破坏规律研究[J]. 煤矿开采,2017,22(3):62−64,103. doi: 10.13532/j.cnki.cn11-3677/td.2017.03.019 LI Lei. Study of overburden broken law of coal seam mining with deep and thick soil layer with shallow depth in Yulin Mine Area[J]. Coal Mining Technology,2017,22(3):62−64,103. doi: 10.13532/j.cnki.cn11-3677/td.2017.03.019
[10] 郭文兵,白二虎,杨达明. 煤矿厚煤层高强度开采技术特征及指标研究[J]. 煤炭学报,2018,43(8):2117−2125. doi: 10.13225/j.cnki.jccs.2017.1573 GUO Wenbing,BAI Erhu,YANG Daming. Study on the technical characteristics and index of thick coal seam highintensity mining in coalmine[J]. Journal of China Coal Society,2018,43(8):2117−2125. doi: 10.13225/j.cnki.jccs.2017.1573
[11] 杜善周. 神东矿区大规模开采的地表移动及环境修复技术研究[D]. 北京: 中国矿业大学(北京), 2010. DU Shanzhou. The rule of the overlying strata movement and their phytoremediation technology after large-scale underground mining operation in Shengdong coal mine area[D]. Beijing : China University of Mining & Technology-Beijing, 2010.
[12] 王志山. 综放高强度开采地表沉陷变形规律实测研究[J]. 矿山测量,2018,46(4):69−72. doi: 10.3969/j.issn.1001-358X.2018.04.018 WANG Zhishan. Study on law of surface subsidence of fully mechanized caving mining with in - situ measurement[J]. Mine Surveying,2018,46(4):69−72. doi: 10.3969/j.issn.1001-358X.2018.04.018
[13] 谭志祥,王宗胜,李运江,等. 高强度综放开采地表沉陷规律实测研究[J]. 采矿与安全工程学报,2008,84(1):59−62. doi: 10.3969/j.issn.1673-3363.2008.01.012 TAN Zhixiang,WANG Zongsheng,Li Yunjang,et al. Field research on ground subsidence rules of intensive fully-m echanized m ining by sublevel caving[J]. Journal of Mining & Safety Engineering,2008,84(1):59−62. doi: 10.3969/j.issn.1673-3363.2008.01.012
[14] 徐乃忠,高 超,吴太平. 浅埋深高强度综采地表沉陷规律实测研究[J]. 煤炭科学技术,2017,45(10):150−154,202. XU Naizhong,GAO Chao,WU Taiping. Study on actual measurement of surface subsidence law occurred by shallow mining depth and high intensity fully mechanized coal mining[J]. Coal Science and Technology,2017,45(10):150−154,202.
[15] 孙庆先. 浅埋煤层综采(综放)条件下地表移动变形特征分析[J]. 煤炭工程,2019,51(3):97−102. SUN Qingxian. Study on characteristics of shallow coal seam surface subsidence under fully-mechanized mining (top-coal caving)[J]. Coal Engineering,2019,51(3):97−102.
[16] 郭文兵,王金帅,李圣军. 浅埋厚煤层高强度开采地表移动规律实测研究[J]. 河南理工大学学报(自然科学版),2016,35(4):470−475. GUO Wenbing,WANG Jinshuai,LI Shengjun. Surveying study on surface movement characteristics of the high intensity mining under shallow depth thick seam[J]. Journal of Henan Polytechnic Univeirsity (Natural Science),2016,35(4):470−475.
[17] 陈俊杰,南 华,闫伟涛,等. 浅埋深高强度开采地表动态移动变形特征[J]. 煤炭科学技术,2016,44(3):158−162. CHEN Junjie,NAN Hua,YAN Weitao,et al. Features of surface dynamic movement and deformation caused by high intensity mining with shallow depth[J]. Coal Science and Technology,2016,44(3):158−162.
[18] 李德海,许国胜,余华中. 厚松散层煤层开采地表动态移动变形特征研究[J]. 煤炭科学技术,2014,42(7):103−106. LI Dehai,XU Guosheng,YU Huahong. Study on features of surface dynamic movement and deformation caused by coal mining under thick alluvium[J]. Coal Science and Technology,2014,42(7):103−106.
[19] 余学义,穆 驰,张冬冬. 厚松散层大采高开采地表移动变形规律研究[J]. 煤矿安全,2020,51(4):235−239,243. YU Xueyi,MU Chi,ZHANG Dongdong. Study on law of surface movement and deformation in thick loose layer with large mining height[J]. Safety in Coal Mines,2020,51(4):235−239,243.
[20] 任永强,时 代. 万利矿区浅部煤层开采地表移动规律研究[J]. 矿山测量,2013(1):59−61. doi: 10.3969/j.issn.1001-358X.2013.01.020 REN Yongqiang,SHI Dai. Study on surface movement law of shallow coal seam mining in Wanli mining area[J]. Mine Surveying,2013(1):59−61. doi: 10.3969/j.issn.1001-358X.2013.01.020
[21] 张安兵,高井祥,张兆江,等. 老采空区地表沉陷混沌特征及时变规律研究[J]. 中国矿业大学学报,2009,38(2):170−174. doi: 10.3321/j.issn:1000-1964.2009.02.005 ZHANG Anbing,GAO Jingxiang,ZHANG Zhaojiang,et al. Chaotic Characteristics and Time-Variable Law of Surface Subsidence of Goaf[J]. Journal of China University of Mining & Technology,2009,38(2):170−174. doi: 10.3321/j.issn:1000-1964.2009.02.005
[22] 张安兵,刘新侠,高井祥,等. 基于经验模式分解和相空间重构的采空区地表动态变形混沌性态及时变规律研究[J]. 岩土力学,2010,31(10):3191−3196. doi: 10.3969/j.issn.1000-7598.2010.10.026 ZHANG Anbing,LIU Xinxia,GAO Jingxiang,et al. Law of time-dependent deformation and chaos characteristics of goaf surface based on empirical mode decomposition and phase space reconstruction technique[J]. Rock and Soil Mechanics,2010,31(10):3191−3196. doi: 10.3969/j.issn.1000-7598.2010.10.026
[23] 胡青峰,崔希民,袁德宝,等. 厚煤层开采地表裂缝形成机理与危害性分析[J]. 采矿与工程安全学报,2013,29(6):864−869. HU Qingfeng,CUI Ximin,YUAN Debao,et al. Formation mechanism of surface cracks caused by thick seam mining and hazard analysis[J]. Journal of Mining & Safety Engineering,2013,29(6):864−869.
[24] 黄庆享,张 沛,董爱菊. 浅埋煤层地表厚砂土层“拱梁”结构模型研究[J]. 岩土力学,2009,30(9):2722−2726. doi: 10.3969/j.issn.1000-7598.2009.09.030 HUANG Qingxiang,ZHANG Pei,DONG Aiju. Mathematical model of “arch beam” of thick sandy soil layer movement in shallow seam[J]. Rock and Soil Mechanics,2009,30(9):2722−2726. doi: 10.3969/j.issn.1000-7598.2009.09.030
[25] 黄庆享,杜君武,侯恩科,等. 浅埋煤层群覆岩与地表裂隙发育规律和形成机理研究[J]. 采矿与安全工程学报,2019,36(1):7−15. HUANG Qingxiang,DU Junwu,HOU Enke,et al. Research on overburden and ground surface cracks distribution and formation mechanism in shallow coal seams group mining[J]. Journal of Mining & Safety Engineering,2019,36(1):7−15.
[26] 钱鸣高, 石平五, 许家林. 矿山压力与岩层控制[M]. 北京: 中国矿业大学出版社, 2010. QIAN Minggao, SHI Pingwu, XU Jialin. Ground pressure and strata control[M]. Beijing: China University of Mining and Technology Press, 2010.
[27] 张连贵. 兖州矿区非充分开采覆岩破坏机理与地表沉陷规律研究[D]. 徐州: 中国矿业大学, 2009. ZHANG Liangui. Overburden failure mechanism and surface movement law due to sub-critical extraction in Yanzhou Mining Area[D]. Xuzhou: China University of Mining & Technology, 2009.
[28] 国家安全监管总局, 国家煤矿安监局, 国家能源局, 国家铁路局. 建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范[M]. 北京: 煤炭工业出版社, 2017. State administration for safety control and management, state administration of coal mine safety, state energy administration, state railway administration. standard of coal pillar design and mining under buildings, water bodies, railvays and main wells and roadvays[M]. Beijing: Coal Industry Press, 2017.
[29] PENG Syd S, XU Feiya. Application research of 3D laser scanning technology in monitoring subsidence area of coal mining[C]. Proceedings of 37th International Conference on Ground Control in Mining, USA, Morgantown, 2018: 266−273.
[30] GU Yuanyuan,ZHOU Dawei,ZHANG Demin,et al. Study on subsidence monitoring technology using terrestrial 3D laser scanning without a target in a mining area: an example of Wangjiata coal mine, China[J]. Bulletin of Engineering Geology and the Environment,2020:3575−3583.
[31] PENG Syd S, Yi LUO. Comprehensive and Integrated Subsidence Prediction Model -CISPM(V2.0)[C]//Proceedings of 3rd Workshop on Surface Due to Underground Mining, 1992: 22−31.
-
期刊类型引用(1)
1. 任广磊，孙华超，杨京华，张凌筱. 基于流固热三场耦合的致密砂岩破裂压力模拟. 力学季刊. 2024(04): 1097-1106 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 150
- HTML全文浏览量: 37
- PDF下载量: 75
- 被引次数: 1
