厚松散含水层失水沉降相似模拟实验研究

陈芳， 张劲满， 徐良骥， 李杰卫， 徐瑞瑞， 张坤

(1.东莞职业技术学院 建筑学院， 广东 东莞 523808； 2.中国矿业大学(北京) 地球科学与测绘工程学院， 北京 100083； 3.安徽理工大学 空间信息与测绘工程学院， 安徽 淮南 232001； 4.安徽理工大学 深部煤矿采动响应与灾害防控国家重点实验室， 安徽 淮南 232001； 5.浙江省地质勘查局 浙江省地矿建设有限公司， 浙江 杭州 310052； 6.安徽省煤田地质局， 安徽 合肥 230088)

0 引言

我国厚松散层矿区分布在东北、华中及华东地区[1-2]，且华东“两淮”矿区的厚松散层分布有丰富的含水层[3]，厚松散含水层地质条件下煤层开采引起的地表沉陷表现为下沉系数偏大、土体受损区域范围广，研究厚松散含水层地质条件下覆岩破断及变形规律具有重要意义[4-5]。

一些学者对厚松散层下开采进行了研究。文献[6-7]通过实测数据研究发现采动后地表含水层失水固结，造成附加采动下沉，导致下沉系数偏大。文献[8-9]研究了含水层失水固结沉降机理，对含水层失水引起的地表下沉进行预计，可提高沉陷预计精度。文献[10]通过FLAC3D流-固耦合计算，得出含水层失水固结条件下最大下沉值较不失水时增加7%。文献[11]针对华东矿区厚松散层失水固结沉降引起的井壁破坏问题，研究了不同水位条件下地表失水固结沉降特征。文献[12]采用现场实测和相似材料模拟实验等手段，针对淮南矿区下沉值偏大现象进行了探究。文献[13]采用理论分析、现场实测及数值模拟等综合手段，探究厚松散层薄基岩综放开采条件下覆岩破断机理及该地层状态下的采动裂隙演化规律。文献[14]指出含水层失水沉降是导致下沉系数偏大的原因，以概率积分法推导出黏土体失水沉降的计算公式。文献[15]指出覆岩损伤是含水层失水的关键，通过数值模拟对采矿引起的覆岩损伤和失水过程进行了研究。文献[16]利用设计的相似材料实验装置对煤矿开采引起的厚松散层条件下承压含水层的失水规律和地表沉降规律进行了探究。文献[17]将基岩和松散层作为2种不同的介质，探究厚松散层矿区岩土体协同作用机理。上述研究对厚松散含水层地质采矿条件下覆岩破断及变形规律未作深入研究。针对上述问题，本文以淮南矿业(集团)有限责任公司潘四东煤矿11111工作面为工程背景，研究了厚松散层条件下含水层失水时覆岩破断及变形规律。根据研究区地层结构构建相似材料模型，采用数字摄影测量提取位移法记录覆岩观测线观测点的破断过程及移动变形情况，在W型及O型剪切应力拱的基础上，根据有效应力揭示厚松散含水层下开采地表下沉值偏大的内在机理，为类似地质采矿条件下安全生产提供理论依据。

1 工程概况

潘四东煤矿11111工作面倾向长度为145 m，实际回采走向长度为411 m，松散层厚度为336 m，属于典型的厚松散层。11111工作面上覆新生界下部含水层的成分是砂土和黏土的混合沉积物，富水性弱。顶部为粉砂岩、细砂岩、砂泥岩互层、黏土层，中部以砂岩为主，夹粉砂岩及泥岩，上部以泥岩、砂质泥岩为主。采煤方法为综采放顶煤一次全采高，顶板管理方式为全部垮落法。研究区地层结构见表1。

表1 研究区地层结构Table 1 Stratigraphic structure in the study area

2 实验设计

2.1 实验模型构建

以潘四东煤矿11111工作面开采情况为模拟对象，选取研究区走向一剖面为原型进行相似材料模拟实验，根据模拟实验原理及实际情况选择1∶100作为几何相似系数，建立相似材料模拟实验模型，如图1所示。11111工作面上覆岩层布设1、2号2条观测线，含水层上方布设3、4号2条观测线，每条观测线上布设23个观测点，共计92个观测点，模型开挖采用2 h模拟实际1 d。

(a) 实验方案

(b) 实验模型

模型观测线具体设置如图2所示。开切眼和终采线均距边界55 cm，回采长度为190 cm，每2 h开采1次，每次推进10 cm，每次推进完成后记录覆岩变化和移动变形情况。

图2 相似材料模型观测线设置Fig.2 Setting of observation lines for similar material models

2.2 覆岩运动观测

在相似材料模型正前方2～5 m处安放工业相机三脚架，通过移动和伸缩三脚架调整工业相机的最佳位置，调整光圈、距离和灯光等直至工业相机记录的覆岩破断及移动变形效果最佳，然后固定三脚架，具体方法为在工业相机下方的悬挂钩装置上配置一定质量的物体，以保持实验过程中三脚架的稳定性(图3)。

图3 工业相机固定Fig.3 Industrial camera fixing

数字摄影测量提取位移法可快速记录覆岩观测线观测点的破断过程及移动变形，如图4所示。将工业相机获取的图像导入GetData数字化软件中，通过模型架上的4个控制点对图像进行校正，建立统一的坐标系统，获取每一幅图像中监测点的坐标。

图4 数字摄影测量提取位移法Fig.4 Digital photogrammetry extraction displacement method

按照开挖时间记录覆岩变化，在实验模型开挖过程中获取了大量的覆岩破断及移动变形图像，并用数字化仪软件对获取的图像进行处理。首先将像平面坐标与物平面坐标进行转换，然后进行数码相机镜头畸变校正，最后根据坐标转换参数和畸变转换参数来计算观测点的下沉值。

3 实验及结果分析

3.1 含水层失水沉降理论分析

随着工作面推进，煤层直接顶开始逐渐垮落，覆岩中主关键层失稳导致上覆岩层分层出现破断。

沿着垮落带两侧底部作切线延伸至含水层，由相似材料模拟实验结果可知垮落角分别为57°和62°。根据覆岩变化过程可知覆岩两端并未发生破断，而在开切眼到采空区中部再到终采线的这部分空间中呈现拉伸-挤压-拉伸的受力状态，对应的剪应力变化为减小→增加→减小→增加，将该剪切应力变化称为W型剪切应力拱。在W型剪切应力拱作用下形成2条纵向的发育至含水层的主导水裂隙带。导水裂隙带的进一步发育引起含水层失水固结，在厚松散层重力作用下进一步压实，增加地表沉降的程度。含水层失水固结及断裂带移动挤压导致含水层与覆岩之间形成一个薄层空间。随着覆岩破断运动的加剧，在弯曲带和覆岩共同挤压下形成O型剪切应力拱，压缩薄层空间，导致地表下沉量增大。覆岩运动及地表沉降过程如图5所示。

(a) W型剪切应力供

(b) O型剪切应力供

工作面开采产生的导水裂隙带发育至含水层时，打破了原有的孔隙水压力平衡状态，加剧了含水层间的流体在水力梯度作用下向采空区的垂向流动，在上覆岩层自重作用下，含水层孔隙体积减小，土体产生固结压缩。含水层失水沉降前覆岩的自重由孔隙水应力和土体颗粒有效应力2部分承受；失水后，含水层中孔隙水承担应力减小，土体颗粒承担应力增大，而工作面开采会打破原岩应力的平衡状态，致使覆岩与含水层产生水力联系，覆岩承担应力增大。在厚松散层地质采矿条件下，松散层厚度及采动程度是含水层失水固结的主要影响因素，松散层厚度影响含水层的水压力，采动程度影响导水裂隙带的发育高度，在孔隙水压力减小且导水裂隙带高度增加作用下，土体颗粒承受的有效应力增大。

工作面开采引起覆岩应力平衡状态破坏，厚松散含水层富含孔隙水，具有成岩晚、孔隙度大、泥质弱胶结、强度低等特点。地下开采使得工作面上方覆岩产生断裂破坏及移动变形，当导水裂隙带发育至厚松散含水层时，含水层中的水就会沿采动裂隙流向采空区，造成岩土体中水位下降。在含水层所受总应力不变的情况下，孔隙水压力降低，导致含水层岩土体颗粒上的有效应力增加。

σ=σ′+(1-β)μ

(1)

式中：σ为总应力；σ′为有效应力；β为孔隙水压力；μ为岩石压缩系数比。

在失水后，含水层固结压密含水层土体颗粒有效应力增大的情况下，含水层岩土体会产生附加下沉量，致使该条件下采矿活动引起的下沉量增大。

(2)

式中：εz为岩土体失水下沉量；αv为岩土体压密系数；Δσ为应力增量；e0为岩土体初始孔隙比。

3.2 失水状态下覆岩损伤分析

在模型达到稳态后，进行11111工作面的开采工作。在实验模型开采工作面的推进过程中，覆岩应力状态处于动态变化之中，11111工作面覆岩演化如图6所示。当工作面推进至距开切眼50 m位置时，顶板初次来压显现，开切眼处的上覆岩层发生破断并垮落至采空区，垮落带上方覆岩出现厚度为2 m的离层裂隙带，垮落带高度为24 m，垮落角为65°；当工作面推进至距开切眼90 m位置时，含水层失水固结及断裂带移动挤压导致含水层和覆岩之间形成一个薄层空间，较厚的松散层增加了含水层的孔隙水压力，导水裂隙带的进一步发育导致基岩上方含水层失水固结，在厚松散层的自重及附加应力作用下进一步压实薄层空间，增加了地表沉降程度；当工作面推进至距开切眼165 m位置时，在终采线外侧10.5 m含水层上方13 m处覆岩发生破断，覆岩产生长20 m的纵向离层裂隙，在距开切眼186 m位置处覆岩发生破断，该破断延伸至含水层下部；在实验模型工作面开采工作完成且覆岩达到稳态后，前垮落角为57°，后垮落角为62°，导水裂隙带高度为63 m，开切眼及终采线上方覆岩在应力集中作用下发生破断，产生纵向裂隙，开切眼及终采线上方垮落带区域内覆岩产生横向离层裂隙，纵向裂隙和横向离层裂隙加剧了覆岩与含水层间的水力联系。

(a) 工作面推进至距开切眼50 m

(b) 工作面推进至距开切眼90 m

(c) 工作面推进至距开切眼165 m

(d) 开采工作完成

3.3 失水状态下覆岩动态运动规律分析

观测实验模型工作面开采过程中观测点的移动变形量，并绘制成图，得到含水层失水状态下各观测线下沉量变化，如图7所示。

(a) 1号观测线动态沉降曲线

(b) 2号观测线动态沉降曲线

(c) 3号观测线动态沉降曲线

(d) 4号观测线动态沉降曲线

由图7可知，不同观测线观测点的下沉量中，随着开采工作面的推进，各观测线覆岩下沉量逐渐增大，接近工作面的1号观测线覆岩下沉量最大，最大值为6 640 mm，大于工作面开采厚度6 630 mm，原因是：原采空区受采动作用致使原先未完全垮落的岩层在其自重的影响下进一步破断压实；覆岩离层裂隙带发育至含水层，致使含水层失水固结及压实，且含水层中下渗的水量致使下部岩层的抗压强度减小，增加了岩层破断弯曲概率。1号及2号观测线在开切眼和终采线附近，其下沉量曲线走势基本类似，3号和4号观测线下沉量曲线走势基本吻合。工作面上方的1号及2号观测线下沉量曲线跳变一致，含水层上方的3号及4号观测线下沉量跳变同步，2号与3号观测线下沉量跳变异步，表明含水层对覆岩移动变形具有重要作用。

4 结论

(1) 根据相似材料模拟实验原理构建了以潘四东煤矿11111工作面为原型的厚松散含水层失水覆岩运动模型，根据模拟实验结果及含水层失水沉降理论分析，提出W型剪切应力拱及O型剪切应力拱的概念，解释了下沉量偏大的现象。

(2) 通过相似材料模拟实验分析研究了该地质条件下覆岩损伤特征。在开采工作面推进过程中，覆岩应力状态处于动态变化之中，在工作面开采工作结束并进入稳态后，前垮落角为57°，后垮落角为62°，导水裂隙带高度为63 m，开切眼及终采线上方覆岩在应力集中作用下断裂，产生纵向裂隙，开切眼及终采线上方断裂带区域内覆岩产生横向离层裂隙，纵向裂隙和横向离层裂隙加剧了覆岩与含水层间的水力联系。

(3) 通过相似材料模拟实验研究了该地质采矿条件下覆岩动态运动规律。随着开采工作面的推进，各观测线覆岩下沉量逐渐增大，接近开采工作面的1号观测线覆岩下沉量最大，大于工作面的开采厚度；1号及2号观测线在开切眼和终采线附近，其下沉量曲线走势基本类似且下沉量跳变一致；3号及4号观测线下沉量曲线走势基本吻合且下沉量跳变同步；2号与3号观测线下沉量跳变异步，表明含水层对覆岩移动变形具有重要作用。