刘一龙 ，杨天鸿 ，马 凯 ，叶 强 ，赵 永 ，赵乾百

（1.东北大学 资源与土木工程学院, 辽宁 沈阳 110819；2.东北大学 岩石破裂与失稳研究所, 辽宁 沈阳 110819；3.西山煤电(集团)有限责任公司 西曲矿, 山西 太原 030200）

0 引 言

矿井水害是制约矿山安全高效生产的主要灾害之一。其中，地表水害作为矿井水害的一种主要形式，是我国西北与华北地区浅埋藏煤层与地表沟谷发育矿区面临的突出问题[1-3]。沟谷地形浅埋煤层顶板岩层多为单一关键层存在或关键层缺失，采动后地表沉陷剧烈。覆岩裂隙发育贯穿基岩直达地表，不仅会产生因采空区漏风引发的煤矸自燃风险，还会破坏浅、表层水资源系统，并成为汛期地表积水的下渗通道[4-8]。此外，沟谷区两侧坡体在煤层采动扰动下易发生失稳破坏，滑坡堆积体堵塞河道还会形成堰塞湖[9-11]。加之近年极端天气频发，如2021 年山西省现“百年一遇”大范围强降水，受汛情影响，山西全省60 座煤矿被迫停产，采动地质灾害防治不当将严重影响生产与安全。因此，研究沟谷地形下煤层采动覆岩破坏特征及其对地表山体滑坡、径流涌水灾害的影响具有重要意义。

目前国内外学者对于浅埋深与沟谷地形下煤层覆岩移动、破断特征以及裂隙扩展贯通开展了诸多研究，但多仅限于单一煤层开采，且未提及其对地表水害产生的影响与防治。张杰等[12]通过多元素拟合公式及有限元差分软件，研究了工作面过沟谷开采过程中隔水岩层的隔水性、周期性“破坏—恢复”变化特点，并提出了沟谷下留设煤柱的防突水开采措施；李建伟等[13]分析了沟谷区域浅埋煤层开采井上下地质灾害发生机理及其相互作用，认为地表沟谷对覆岩应力的影响作用随着埋深增加而减弱，并提出了针对地表滑坡预防性控制措施；赵杰等[14]以现场煤层地质条件为基础，分析了浅埋煤层过沟谷开采时覆岩关键层破断形式及失稳结构特征，建立了沟谷坡角与关键层破坏类型的关系；刘辉等[15]通过对西部矿区浅埋煤层开采造成的地裂缝的现场观测和分析，提出了采动过程中的临时性地裂缝治理标准和依据，并提出了初步的控制方法；孙魁等[16]研究了“对滑型”黄土滑坡的形成过程，提出了“对滑型”滑坡的形成是由黄土坡体、沟谷地貌、采矿扰动、降雨催化以及相互之间的耦合作用导致。

目前我国煤炭资源开发已进入深部开采阶段，煤层重复采动作用下沟谷区坡体结构损伤加剧，开展沟谷地形与煤层重复采动耦合作用下的覆岩破坏特征研究更符合工程实际需求。同时，针对汛期降雨引发径流水害机理以及沟谷区水害防治方法也需深入开展研究。因此，笔者以西曲煤矿典型的沟谷区多煤层开采为工程背景，综合应用集地表勘察、In-SAR 动态观测、降雨−径流分析和数值模拟为一体的“覆岩移动−地表变形−径流积水”分析方法，探究沟谷地形下的“裂隙发育—山体滑坡—河道堵塞—涌水加剧”灾害链形成过程，最终提出基于裂隙发育与地表淹没范围的沟道防渗方法，以期为沟谷区煤层开采造成的地表地质灾害防治与水资源保护提供借鉴和指导。

1 地质概况与开采条件

研究区位于山西省古交市西曲煤矿，矿区范围内沟谷纵横，地形切割剧烈（图1）。目前西曲煤矿上组煤2+3 号煤层大部分资源已于2010 年前开采完毕，主采下部8 号煤层。本文即以8 号煤层的28306 号工作面与其上覆的2+3 号煤层形成的多煤层开采为研究工况。28306 号工作面位于西曲矿矾石沟沟谷区，矾石沟为西曲矿东部主要沟谷，贯穿井田南北，属季节性沟谷，沟谷平时水量不大，雨季有水。雨季多集中在7、8、9、10 四个月，年平均降水量511.5 mm，降水最终汇入汾河。28306(1)号工作面因遇到风氧化带仅推进了130 m，风氧化带的煤层、岩层均较为破碎，裂隙发育，对井下安全生产带来威胁。为确保安全，此处留设了一段近50 m 煤柱，之后继续开采28306(2)号工作面。在28306(2)号工作面推进15 m 时，地表为一条矾石沟支沟低洼处，支沟倾向南东和东，水流易汇集沿支沟倾斜方向流入东部矾石沟内。

图1 研究区地质与开采条件Fig.1 Geological and mining context of the study area

28306 号工作面于2021 年完成开采，开采方向由北向南，该工作面地表覆盖层厚度90～188 m，西侧为2017 年已开采结束的28307 采空区，东侧为未开采的28305 采区，其上覆靠近沟谷区的2+3 号煤层由小煤窑巷柱式开采，其余2+3 号煤层由西曲矿机械化开采。8 号煤层厚度稳定，煤厚3.90～4.20 m，平均4.10 m，煤层整体倾向西南，倾角1°～8°，一般为4°，皆采用倾斜长壁后退式、一次采全高、全部垮落法综合机械化采煤。

图2 为研究区钻孔综合柱状图。8 号煤上方143 m 范围内赋存有5 层较厚顶板（厚度不小于5 m），其中有3 层顶板岩性为泥岩（厚度小于10 m），1 层近地表的顶板岩性为细粒砂岩（厚度为5 m），顶板中19.64 m 厚的粉砂岩呈深灰−灰黑色，以石英为主，较硬，全层普遍呈致密状，含有较少的孔隙和裂缝，再由关键层的定义与变形特征[17-18]，对于多层顶板组合结构的承载满足如下公式：

图2 综合地层柱状图Fig.2 Synthetical stratum histogram

式中：(qn)m为第n层顶板对第m层顶板的载荷；Em为第m层岩层的弹性模量；γi为第i层岩层的容重；hi为第i层岩层的厚度。

对于满足式（2）的第n层岩层，可判定其为关键层：

综上，判别19.64 m 厚的粉砂岩地层为煤层间关键层。8 号煤顶板地层总体属硬质岩，2+3 号煤上部地层主要为砂泥岩互层。

2 沟谷区岩层移动与地表变形特征

2.1 地表破坏类型及特征

图3 为研究区现场勘察和航拍影像，由图3 总结出了28306 号工作面开采引发的地表破坏类型主要为无明显落差的地裂缝、台阶状裂缝、滑坡与堰塞湖。在开切眼附近采空区边界上方，地表裂隙呈现明显发育的现象（图3a）。这些单条裂缝的宽度在400～500 mm，且宽窄不一，顶端尖灭。裂缝群则呈现弧形或椭圆形的展布特征。在同一地区，大部分地裂缝的展布方向近乎平行。但在矾石沟支沟区域，由于地形的陡峭山坡，地势西高东低，因此地裂缝的展布方向比较杂乱。台阶状裂缝出现在井下采煤工作面切眼位置及采区边界附近的地表，是由浅埋煤层开采顶板初次切落所造成的（图3b）。这些台阶状地裂缝的深度不尽相同，一般不超过3 m。地裂缝在垂直方向上呈现“V”字形，上部较宽，下部较窄，两壁陡峭，呈现出明显的台阶状，其裂缝长度一般在几米到几十米不等。在采空区上方地形切割强烈的沟谷边坡地带地表出现黄土层及砂层崩塌、崩滑，崩塌块呈破碎状，完全丧失原有的结构（图3c）。至2021年7 月31 日工作面回采结束，地裂缝已扩展至全地表。沟谷区地下煤层开采导致地下垮塌引发覆岩渐进式破坏，进一步引发地表沟谷坡体的地质结构改变，使得地表坡体的稳定性受到影响。降雨等自然力的作用加剧了地表坡体的滑动和崩塌（图3d），堵塞河道还引发了土石流和堰塞湖等次生地质灾害（图3e、图3f）。

图3 研究区现场勘察和航拍影像Fig.3 Field survey and UAV image in the study area

2.2 地表沉降的空间扩展过程

InSAR（Interferometric Synthetic Aperture Radar）技术是一种利用合成孔径雷达（SAR）成像的干涉测量方法，可以实现毫米级的地表形变监测，具有范围广、时间跨度大、受气候干扰小等优势，已成功应用于观测滑坡的动态变形过程[19-20]和采煤引起的地表沉降[21-22]。地表沉降是深部岩层变形、破断影响至地表的体现，具有明显的滞后性，且滞后周期较长。28306(2)号工作面开采时间为2020−10−31—2021−07−31，因此分析了2021 年1 月至2022 年9 月期间的时间序列InSAR 形变监测结果，如图4 所示，28306(2)号工作面地表累积沉降自北向南扩展，与开采推进方向一致。图4 为基于InSAR 数据的研究区累计沉降，图4a 展示了此次研究的InSAR 形变范围，至2021 年1 月9 日，工作面累计推进约176.23 m，开采引发的地表沉降已开始显现（图4b），此阶段岩层移动为非充分采动阶段；至2021 年2 月26 日，工作面累计推进约348.18 m，此阶段地表沉降值增大，时间序列InSAR 形变监测得到的地表累计沉降值已达到−69.4 mm（图4c）；至2021 年4 月15 日，工作面累计推进约494.5 m，岩层移动为充分采动阶段，沉降范围扩大（图4d）；至2021 年7 月31 日工作面完成回采，此时累计推进约680.52 m。由图4e 可知，28307 号采空区在28306 号工作面采动影响下地表发生残余沉降，地下开采造成的塌陷面积明显大于实际开采面积，但与平缓地表形成的椭圆形塌陷区不同。本次塌陷区形状与沟谷地形走势大体相同，特别是28306 号工作面东侧由北向南延伸的矾石沟沟谷处的InSAR 结果位移云图，具有明显分界特征。尽管28306 号工作面的煤炭开采活动已终止，但地表下沉量仍在增加。至2022 年9 月，28306 号工作面的最大地面沉降量达到–98.66 mm（图4f）。

图4 基于InSAR 数据的研究区累计沉降Fig.4 Cumulative subsidence in the study area based on In-SAR data

为分析地下煤层开采后造成地表沉降的空间扩展过程，选取了AA’主断面上最大下沉监测点，并绘制了随时间变化的沉降曲线（图5），研究区地表沉降总体可划分为3 个阶段：第1 阶段为初始沉降阶段（2021−01−09—2021−04−15），此阶段煤层开采过后，致使煤层顶板受到破坏，覆岩渐进式坍塌引起地表沉降，这一阶段累计沉降指标为10 mm；第2 阶段为加速沉降阶段（2021−04−15—2021−10−09）。在煤层开采达到充分采动后，地表移动变形活跃，采空区周围的地层受到压缩产生变形，地表沉降持续时间较长；第3 阶段为长期沉降阶段，时间尺度拉长，地下岩层移动会趋于稳定。虽然地下空洞的存在会使地表沉降持续进行，但沉降速度逐渐减缓，沉降量也会逐渐减小，6 个月的沉降值不超过30 mm。

图5 地表最大沉降点累计沉降值Fig.5 Cumulative settlement value of the maximum settlement point on the ground

自2022 年2 月8 日后，地表沉降达到相对稳定状态，估算终采后地表达到充分沉降的滞后周期为5～7 个月，较该矿区记录的中深埋煤层终采后地表达到充分沉降的滞后周期短了1～2 个月。受限于现有SAR 卫星的极轨飞行和侧视成像模式，In-SAR 形变观测值对南北向形变不敏感，导致该方法在对沟谷区地表的应用中，存在局部入射角过大/小的问题，进一步导致视线向位移与实际位移相差较大。但InSAR 技术可展示煤层地下开采造成地表沉降的空间扩展过程的持续性与阶段性，进一步得出煤层终采后地表达到充分沉降的滞后周期，指导未达到充分沉降范围内的地表河道优先采用“柔性治理”技术。

2.3 覆岩垮落和变形特征

在分析岩层渐进式破坏过程中，三维离散单元法程序3DEC 应用广泛，此软件较好还原出了材料介质破坏过程中的不连续特性。基于试验采场的地表等高线和钻孔柱状图，沿着图1 中的I—I’和II—II’剖面分别建立了用于3DEC 数值计算的三维地质模型。图6 为数值计算模型及开采范围，图6a 中的数值方案1 主要用于研究28306 号工作面推进过程中覆岩垮落特征及矾石沟支沟处坡体破坏特征，图6b中的数值方案2 主要用于研究向沟开采对河床的稳定性影响，同时考虑了上部2+3 号煤层巷柱式开采工况。

图6 数值计算模型及开采范围Fig.6 Numerical calculation model and mining range

考虑到块体数量将影响计算的效率，对厚度太薄岩层按照四舍五入的原则进行合并和简化，并将数值计算模型简化为0°。数值模型中涉及到的岩性主要有5 种（砂岩、砂泥岩互层、砂质泥岩、煤和石灰岩），这5 种岩体力学参数见表1。所用参数的可靠性已在此矿山的已有文献[23-24]中得到验证。模型块体采用莫尔−库伦模型，节理模型采用节理面接触库伦滑移模型。模型四周均为法向位移约束，模型底界面限制水平方向和垂直方向位移，地表曲面为自由边界。煤层开采方向从左至右，详细的开采参数如图6 所示。

表1 煤岩物理力学参数Table 1 Physical and mechanical parameters of coal and rock

图7 为岩层移动与地表变形特征数值方案1 计算过程中，首先对2+3 号煤层进行开采，然后对8 号煤层进行开采。结合矿山现场勘察报告，浅部2+3号煤层开采后形成了无明显落差的地裂缝，塌陷坑较少，偶发坡体滑移。下部煤层28306(1)号工作面推进至20 m 时顶板初次垮落（图7a），这与现场记录的工作面顶板初次垮落步距（19 m）接近，且小于西曲矿8 号煤平均顶板初次垮落步距35～40 m，主要原因在于顶板存在风化现象。28306(1)号工作面推进至100 m 时关键层发生破断（图7b），这与CHENG 等[23]通过微震监测得到的西曲矿8 号煤顶板关键层发生破断结果接近。关键层破断后，下部煤层开采扰动已影响至上覆地表，地表发生二次变形、损伤，并产生明显地表沉降（图7c）。28306(2)号工作面推进至110 m 时，关键层以下形成大面积离层，继续开采至120～130 m 时，矾石沟支沟北侧坡体在采动与地形拉伸作用下向临空面一侧发生剧烈变形（图7d），重复采动后地表具有明显落差的台阶裂缝数量增加，塌陷坑范围扩大，地表下沉加剧，地表岩层移动范围扩大，沟谷坡体产生大量裂隙。

图7 岩层移动与地表变形特征Fig.7 Characteristics of strata movement and surface deformation

浅埋沟谷地形下煤层开采后，围岩应力重新分布，岩层不断变形并周期性破断，煤层间关键层以下岩层垮落后岩层移动范围无明显外扩特征，此范围的岩层张拉破裂形成的裂隙居多，剪切破裂次之，以竖向裂隙为主，主要集中在顶板上方23 m 范围内（图7e）；相反，关键层以上岩层剪切破裂形成的裂隙占比增加，关键层破断后覆岩裂隙不断向上发育，如图8 所示。岩层移动破坏了沟谷两侧坡体稳定性，图8 中沟谷地形对浅部岩层移动与裂隙发育特征的影响作用加强，在岩土体重力作用下坡体易发生剪切滑移。同时，坡体地表形成的大量裂隙向深部扩展并与覆岩上行裂隙贯通，由此沟通地表形成导水通道，这些裂隙对地表降水的截流作用会造成井下涌水量的突增。

图8 覆岩裂隙与坡体裂隙发育特征Fig.8 Development characteristics of overlying strata and slope fractures

图9 为重复采动作用下覆岩二次破坏特征，如图9 所示，数值方案2（图1 中II—II’剖面）中首先对上部2+3 号煤巷柱式开采，每开采20 m 后留设20 m煤柱，然后开采其下部的28306 号工作面，开采方向皆朝向矾石沟方向。因浅部煤层“巷柱式”开采，留设的大量煤柱较大程度地减少了应力集中的扰动，覆岩垮落并不充分。在28306 号工作面开采结束后，由于层间关键层破断，浅部2+3 号煤采空区在重复采动作用下发生“活化”，地裂缝剧烈发育，成为了地表水主要下渗通道。

图9 重复采动作用下覆岩二次破坏特征Fig.9 Secondary failure characteristics of overlying strata under repeated mining

28305 号工作面（待采）位于28306 号工作面东侧，且处于矾石沟沟谷正下方。通过离散元程序模拟28305 号工作面开挖，结果如图10 所示，28305 号工作面过沟开采至100 m 时关键层下方产生大范围离层裂隙，关键层弯曲的变形引起了地表沉降与坡体变形。当工作面推进至120 m 时，地表河床发生大范围沉降，采动裂隙贯通至地表以及两侧山坡。28305 号工作面的开采将影响地表河床稳定性，引发沟谷水下渗，因此，开展河床加固以及沟谷区地表防渗治理是必要的。

图10 沟谷下开采引发的地表沉降与裂隙发育特征Fig.10 Characteristics of surface subsidence and fractures development caused by mining under the gully

3 矾石沟降雨−径流特征及径流水害防治

3.1 矾石沟小流域降雨−径流模拟

采用河道水力计算程序HEC-RAS 开展矾石沟小流域降雨−径流模拟。软件二维水动力学模型的数值计算混合了有限体积法和有限差分法，可有效计算二维圣维南方程组[25]。软件中水动力研究模式，可对研究区流域实现降雨影响分析，进行稳定、快速的洪水淹没模拟。非恒定流模拟的计算公式采用连续性方程（式（3））和动量方程（式（4））。

式中：H为水面高程，m；h为水深，m；u，v分别为x、y方向上的平均流速，m/s；q为项源，代表降雨等外部信息。

式中：V为流速，m/s；g为重力加速度，m/s2；vt为水平方向运动黏度，m/s2；cf为河床底部糙率；f为科里奥利系数；k为垂直方向单位矢量。

高精度的DSM 可满足HEC-RAS 模拟小流域的洪水淹没的精度，研究采用无人机倾斜摄影测量获取了研究区DSM 数据。设置二维网格的空间步长为5 m，网格数量为24 428 个，根据研究区植被覆盖和土地利用类型，选取河道糙率值为0.035，草地糙率值为0.04，耕地糙率值为0.03。本案例模型选取全局降水为入流条件。下边界为河道出口，并假设此处因采动滑坡引发堰塞湖。全局降水数据选取2021 年10 月3 日00:00 至2021 年10 月6 日23:00共4 天95 h 的降雨数据（表2），模拟计算间隔为1 h，输出水力要素间隔为1 h。

根据二维非恒定流模拟结果，2021 年10 月4 日至10 月6 日不同时刻降水淹没深度对比如图11 所示（图11a 为4 日0 点降雨淹没范围；图11b 为5 日0点降雨淹没范围；图11c 为6 日0 点降雨淹没范围），可以看出矾石沟汇水由高流向低洼河道，下游堰塞湖积水深度最大处近8 m，成为井下涌水的潜在危险源。同时得到了降雨期间全局淹没时间百分比图（图11d），矾石沟及其支流的淹没时间百分比超过50%，这些区域也在煤矿开采引起的最大沉降区范围内，满足这2 个特征的区域被称为淹没及最大沉降区，这需要保证该区内沟谷河床稳定性。尽管冲沟边坡两侧的淹没百分比仅为20%，但采动滑坡作用下坡体裂隙发育，需要及时填埋沟谷两侧边坡裂缝，符合这2 个特征的区域被称为淹没及裂缝发育区，这是防渗的重点区域，地裂缝采用分层填埋，而其余区域可以被称为径流区，该区域内的地裂缝可直接单层填埋。

图11 降雨引发的研究区地表淹没过程Fig.11 Ground inundation process caused by rainfall in the study area

3.2 矿山地表径流水害防治方法

开展矿山地表径流水害防治的研究不仅会减小地表水引发矿山涌突水事故风险，通过地表裂缝填埋和河床防渗加固，也极大减少了煤层开采对地表水、土壤等环境造成影响，可以有效减少开采过程中对地表水的污染，保证地表水的质量与水资源的可持续利用。通过前述研究总结得到了矾石沟致灾因素主要为多煤层开采形成的地表张拉裂缝、采动沉降引发的边坡失稳以及汛期沟谷积水下渗。本小节针对上述问题提出了以下治理措施，并取得了良好的应用效果。

针对沟谷区不再发育且裂缝宽度小于0.1 m 的地表裂缝可直接填埋，当地表裂缝宽度大于0.1 m时，采用由浅至深的黄土表土层+三合土防渗层+碎石块衬垫层的分层填充，各分层厚度可控制在0.3～0.5 m。沟谷两侧坡体上的地裂缝在蠕变与滞后沉降作用下仍有加剧破坏趋势，提前对潜在滑坡区实施削坡处理，削坡后的坡角为30°～45°。对沟谷区河床防渗治理为主，对地表积水疏放为辅。在工作面开采之前，对地表低洼区域进行填平处理，尽量保持河床坡度一致。实施土方工程延河道加高构筑堤坝，防止河床大面积塌陷后沟流外溢。在工作面开采过程中，及时清理沟谷两侧滑移堆积体；最后，在构筑好的堤坝一旁备足填堵裂缝用的沙土、三合土和碎石。对受山体滑坡影响，易堵塞、积水的河道采用涵管引流的方式进行治理，保证正常行洪，减小井下入渗量。

对沉降未稳定的河床（停采后5～7 个月内）采用红黏土、尼龙网柔性材料进行治理，受重复采动影响河床不易开裂，防渗效果好，治理方法如下：①清理河道碎石，整平河床形成河道基槽并扩宽至4～6 m以上，在河道弯曲系数较大区段进行河道改直；②河道铺设200 mm 厚度黏土垫层后夯实，在垫层上方铺设一层钢塑土工栅格，上方再铺设一层复合土工膜加强防渗水效果，覆盖500 mm 厚红胶土后夯实；③上方再交替施工两层由钢塑土工栅格和红胶土组成的700 mm 厚加筋土层并压实，形成可靠的防渗层；④在表层红胶土上方铺设200 mm 厚碎石并压实嵌入红胶土，提高河床表面强度，防止流水冲刷及车辆行驶破坏河道；⑤河道两侧采用抗紫外线生态袋垒筑800 mm 高度的护坝，形成护坝后进行充填，与山脚形成一定坡度，确保山体流水能正常汇入河道；⑥河道与各支沟交汇处留设汇水口，汇水口两侧采用片石混凝土砌坝20 m，确保各支沟汇水畅通，汇水口两侧不被流水冲刷破坏。实施“柔性”技术治理后的河床断面及治理效果如图12 所示。

图12 河床断面及加固效果Fig.12 River bed section and reinforcement effect

4 结 论

1）通过对矾石沟地表破坏类型的多手段综合分析及运用离散元法数值模拟计算，揭示了多煤层开采地表沉降与沟谷采动滑坡耦合作用下的覆岩破坏及裂隙发育特征。煤层间关键层以下岩层垮落后岩层移动范围无明显外扩特征，以张拉破裂形成的竖向裂隙发育为主。关键层上方岩层剪切破裂形成的裂隙占比增加，在采动沉降效应与重力作用下沟谷坡体易发生剪切滑移，导致坡体下行裂隙与覆岩上行裂隙贯通，贯通裂隙对地表降水的截流作用会导致井下涌水量增大。

2）通过降雨−径流模拟与矾石沟高精度的DSM再现了山西“百年一遇”暴雨期间矾石沟流域不同历时雨水淹没范围与淹没时间百分比。坡体两侧的径流区属于重点防渗区，需要及时采取裂缝填埋、防渗等防治措施。

3）沟谷区地下多煤层开采引发的覆岩沉降效应对地表地貌造成不可逆影响，会引发“裂隙发育—山体滑坡—河道堵塞—涌水加剧”复杂灾害链。提出了沟谷地形下地表径流水害综合防治方法，即径流区、淹没及裂隙发育区地表裂缝填埋、潜在滑坡区削坡治理、淹没及最大沉降区河床“柔性”加固。不仅减小了地表水引发矿山涌突水事故风险，也保证了地表水的质量与水资源的可持续利用。

4）后续将重点探究汛期地表水入渗采场机理及井下涌水量变化规律，指导矿山企业探放水。