特厚煤层综放开采底板变形破坏深度研究

侯俊华

（兖煤菏泽能化有限公司 赵楼煤矿，山东 郓城 272400）

0 引言

随着华北型煤矿开采年限增加，延深带压开采带来的底板水害威胁问题逐渐引起重视。在特厚煤层综放开采采动影响下，工作面底板岩层受矿山压力及卸载影响而发生变形破坏，产生的裂隙很有可能沟通底部含水层。准确获得采动底板变形破坏深度对于特厚煤层综放开采底板水害防治具有重要意义[1-2]。

许多学者对特厚煤层综放开采底板变形破坏深度进行了研究。刘杰[3]采用现场微震监测技术对特厚煤层综放工作面顶底板岩层破坏情况进行了监测分析。朱术云等[4]应用应变感应法和钻孔成像技术对特厚煤层综放采动底板岩层变形破坏特征进行了现场实测。王一栋等[5]和李进军等[6]采用数值模拟方法获得了煤层底板变形破坏深度。杨本水等[7]在对特厚煤层底板地层微观及水理性质分析的基础上，采用预埋分布式光纤感应法和孔间电阻率CT法确定了采动底板变形破坏最大深度。范红伟等[8]通过理论分析、数值模拟和现场实测，对近距离特厚煤层采动影响下底板变形破坏深度进行了研究。张郑伟[9]和陈洋[10]采用应变感应法获得了特厚煤层综放开采条件下底板变形破坏深度。高银贵等[11]采用光纤传感技术和高密度并行电法勘探技术对特厚煤层工作面底板破坏特征进行了研究。王升阳等[12]采用数值模拟和现场实测的方法对煤层底板采动破坏及渗流特征进行了分析。郭国强[13]采用钻孔应力-应变技术开展现场监测，获得了三位一体“空间-时间”全方位底板扰动破坏数据并分析规律。本文在上述研究的基础上，以兖矿能源集团股份有限公司东滩煤矿特厚煤层综放开采工作面为背景，采用现场实测、数值模拟和理论分析方法对该工作面底板变形破坏深度进行研究，以期为矿井特厚煤层综放开采底板水害防治提供量化依据。

1 工程背景

东滩煤矿1305工作面位于矿井一采区，工作面东起开切眼（运输巷侧切眼外帮距FS18断层10 m，轨道巷侧较运输巷侧向东调斜30 m）；西至设计终采线（运输巷侧位于F2断层以东10 m处，轨道巷侧较运输巷侧向西调斜30 m）；北邻1306综放工作面（未采）；南邻1304综放工作面采空区。工作面标高为-510～-590 m，平均标高为-550 m；其上方地面标高为+46.38～+50.34 m，平均标高为+48.36 m。工作面开采山西组3号煤层，厚度为8.1～9.5 m，平均厚度为8.8 m，采用综放开采方式。

工作面范围内煤层沉积稳定，煤岩层走向NE-SW，倾向SE；煤层倾角为0～12°，平均倾角为6°。基本顶由灰至灰白色中细砂岩组成，岩性致密、坚硬，厚度为16.3～23.2 m，平均厚度为20.2 m；直接顶由深灰色粉砂岩组成，厚度为1.3～5.2 m，平均厚度为3.2 m，呈东薄西厚变化特征；直接底由深灰色粉砂岩组成，厚度为3.6～4.5 m，平均厚度为4.0 m；基本底为强度较高的灰白色中细砂岩，东厚西薄，厚度为7.2～18.4 m，平均厚度为12.8 m。工作面顶底板地层柱状图如图1所示。

2 现场实测

2.1 测试钻孔及传感器布设

根据应变感应法及钻孔成像技术特征[4]，结合东滩煤矿开采实际条件，在不影响正常生产的情况下选择在1305工作面下部的东翼第一回风巷中布设2个测试钻孔（图2）。1号测试钻孔中布设应变传感器对底板变形破坏深度进行实测，2号测试钻孔用于成像观测，2个测试钻孔水平间距约10 m。1号测试钻孔孔口在底板与巷道交接处，垂直于巷道方向向工作面底板施工；2号测试钻孔在距离巷道底板约2 m处垂直巷道方向向上施工。测试钻孔剖面如图3所示，其中1号测试钻孔布置的4个应变传感器自上至下依次为测点1-4，1-3，1-2，1-1，分别位于3号煤层底板下垂深16，20，24，29 m处。测试钻孔主要参数见表1。

表 1 测试钻孔主要参数Table 1 Main parameters of test boreholes

表 2 工作面采动底板变形感应距Table 2 Deformation induction distances of mining floor of working face

图 1 1305工作面顶底板地层柱状图Fig. 1 Strata column of 1305 working face roof and floor

图 2 测试钻孔布置平面Fig. 2 Plane of test boreholes layout

图 3 测试钻孔剖面Fig. 3 Profile of test boreholes

2.2 实测结果分析

整个测试过程历时46 d，1号测试钻孔共进行30次应变测试，2号测试钻孔进行了4次成像观测。

1号测试钻孔各测点应变增量随工作面距传感器距离变化曲线如图4所示（每个传感器有10个通道，其中8个工作通道，1个地线通道，1个公共线通道；8个工作通道中1和5、2和6、3和7、4和8是对称的）。测点1-4在距工作面115 m出现了较明显的应变感应，直至距工作面44.9 m时，该测点的应变感应信号全部中断，认为是由于受底板破坏影响剧烈，导致传感器受损[7，9-10]；与测点1-4相比，测点1-3、测点1-2和测点1-1的应变感应程度较弱，应变增量明显变小，超前距分别为110，102，96 m。随着工作面逐渐推进到钻孔位置和推过钻孔，4个测点的应变增量差异性较明显：测点1-4应变增量最大，且在距工作面44.9 m位置出现信号中断现象；其他3个测点应变增量较小，且在推过一定距离后各个通道均有可靠的数据。

图 4 1号测试钻孔各测点应变增量随工作面推进距离变化曲线Fig. 4 Strain increment curves of measuring points in No.1 test borehole during working face advance

图 5 工作面推进不同距离时底板下垂深10 m处2号测试钻孔孔壁图像Fig. 5 Images of No.2 test borehole wall at 10 m vertical depth of floor with different working face advance distances

根据图4可得工作面采动底板变形感应距（表2）。可见采动矿压在底板的传播具有较好的周期性，不同深度应变传感器所反映的水平方向超前距及滞后距基本一致。

工作面推进不同距离时底板下垂深10 m处2号测试钻孔孔壁图像如图5所示。2号测试钻孔距工作面50 m时，底板下垂深10 m处发生塑性破坏，出现裂隙；随着工作面不断推进，该处破坏程度越来越大，裂隙逐渐扩展直至贯通。这与1号测试钻孔的应变感应具有很好的一致性：测点1-4距工作面44.9 m时出现了传感器损坏，表明该处底板存在较强的塑性变形，底板变形破坏深度已达16 m；底板下垂深20 m以下各测点的应变增量较小且同步，表明底板岩层以弹性变形为主。

将应变测试和钻孔成像结果对比分析，可获得1305工作面综放开采底板变形破坏深度不超过20 m。

3 数值模拟

3.1 模型建立和边界条件

结合1305工作面开采实际和地层赋存条件，进行3号煤层顶底板地层岩组划分，建立工作面顶底板工程地质模型，尺寸为400 m×300 m×154 m（长×宽×高），模型纵剖面如图6所示。工作面顶底板岩层物理力学参数见表3。模型边界条件：侧向边界水平方向固定，垂直方向自由；底面边界垂直和水平方向均固定；基于埋深考虑，模型顶部施加11.0 MPa等效荷载。采用FLAC3D软件，模型共划分91 362个单元、96 596个节点。

图 6 工作面顶底板工程地质模型纵剖面Fig. 6 Profile of engineering geological model for working face roof and floor

表 3 工作面顶底板岩层物理力学参数Table 3 Physical and mechanical parameters of working face roof and floor strata

3.2 模拟结果分析

3号煤层开采完毕后塑性区分布如图7所示。由图7（a）可得底板变形破坏深度为20 m；由图7（b）-图7（e）可知，底板下4 m内呈面状破坏，底板下12 m呈环状破坏，且采空区四周较大范围呈剪切破坏，底板下18 m破坏范围明显变小，底板下20 m基本没发生破坏。表明随着距工作面底板距离越远，受采动矿压影响越小，塑性区范围越小，工作面综放开采底板变形破坏深度不超过20 m。

图 7 3号煤层开采完毕后塑性区分布Fig. 7 Distribution of plastic zone after No.3 coal seam mining

4 理论计算

根据断裂力学和弹塑性理论[14]，可得底板变形破坏深度计算公式：

式中：h为底板变形破坏深度，m；x0为塑性区宽度，即煤层屈服区长度，m；φ0为底板岩体平均内摩擦角，（°）。

根据采动矿压理论[15]，塑性区宽度计算公式为

式中：M为采高，m；φ为煤层与顶底板间的内摩擦角，（°）；K为应力集系数；ρ为顶板岩层平均密度，kg/m3；g为重力加速度，m/s2；H为平均采深，m；C为煤层与顶底板间的黏聚力，Pa。

综合1305工作面地质及开采条件，基于式（2）得到的塑性区宽度计算结果，选取经验公式x0=0.015H，并代入式（1），可得

根据1305工作面实际条件，平均采深H=600 m，底板岩体平均内摩擦角φ0=38°，代入式（3）可得1305工作面底板变形破坏深度为19.2 m。

5 结论

（1） 采用应变感应法和钻孔成像技术对1305工作面进行现场实测，发现工作面推进过程中，底板受采动矿压作用在水平和垂直方向上存在明显的显现特征，水平方向上超前测点50 m附近、底板浅部10 m深度位置开始受采动矿压影响，工作面推过一定距离后底板变形破坏剧烈；底板不同深度水平方向上超前距和滞后距变化范围分别为96～115 m和48～52 m，工作面综放开采底板变形破坏深度为16～20 m，垂深20 m以下底板岩层以弹性变形为主。

（2） 通过数值模拟分析了采动底板不同深度塑性区分布特征，发现底板下18 m破坏范围明显变小，底板下20 m基本没发生破坏。

（3） 通过理论计算获得该工作面底板变形破坏深度为19.2 m。

（4） 通过现场实测、数值模拟和理论计算3种方法综合分析，可知1305工作面综放开采底板变形破坏深度不超过20 m。