采空区垮落顶板形态及其演化特征

李 杨,任玉琦,王 楠,金向阳,欧祥吉,罗军波,梅灿旗

(中国矿业大学(北京) 能源与矿业学院,北京 100083)

摘 要:揭示采空区垮落顶板形态特征对于近距离煤层上行开采、含水层下开采以及废弃采空区瓦斯分布等具有十分重要的意义。从开滦矿区12-1号煤层开采后垮落顶板的现场地质雷达探测结果出发,综合运用材料力学、矿压理论与UDEC数值模拟的方法,对采空区垮落顶板形态的典型特征、形成机理进行了充分研究,描述了采空区“正-倒置三角形”垮落顶板形态动态演化特征,给出了采空区“正-倒置三角形”垮落顶板形态的概念。研究结果:① 通过地质雷达探测反演结果可知,煤层顶板分层明显,采空区垮落顶板呈现正置三角形区域与倒置三角形区域交替出现的形态,其中正置三角形区域内下位岩块堆砌整齐、上位岩块堆砌杂乱,而倒置三角形区域内下位岩块堆砌杂乱、上位岩块完整性好;② 采空区垮落顶板形态的出现是基本顶及其控制的上覆岩层周期破断与已破断的关键岩块再次破断共同作用的结果,当各关键岩块的极限抗弯强度小于受到的最大弯矩时即发生再次破断,且基本顶岩块与上覆岩层发生再次破断的位置不同,2者再次破断位置的差异导致采空区垮落顶板被分为正置三角形与倒置三角形区域;③ 提出了采空区垮落顶板“正-倒置三角形”形态的概念,当基本顶及其控制的上覆岩层分层明显时,在煤层开采工作面周期来压后,基本顶及其控制的上覆岩层关键岩块在最大弯矩的作用下发生再次破断,周期破断与再次破断位置将采空区垮落顶板呈现正置三角形与倒置三角形区域,形成“正-倒置三角形”的垮落顶板形态;④ 运用UDCE模拟获得了采空区垮落顶板的形态特征,一定程度上验证了采空区“正-倒置三角形”垮落顶板形态的存在性与合理性。

关键词:采空区;垮落顶板;演化特征;上覆岩层;正-倒置三角形

中图分类号:TD32

文献标志码:A

文章编号:0253-9993(2021)12-3771-10

收稿日期:20210208

修回日期:20210903

责任编辑:王晓珍

DOI:10.13225/j.cnki.jccs.2021.0272

移动阅读

基金项目:国家自然科学基金面上资助项目(52074293);河北省自然科学基金资助项目(E2020402041);国家自然科学基金重点资助项目(51934008)

作者简介:李 杨(1982—),男,河北唐山人,副教授,博士生导师,博士。E-mail:liyangcumtb@163.com

通讯作者:任玉琦(1993—),男,山西大同人,博士研究生。E-mail:ryq2019cumtb@126.com

引用格式:李杨,任玉琦,王楠,等. 采空区垮落顶板形态及其演化特征[J]. 煤炭学报,2021,46(12):3771-3780.

LI Yang,REN Yuqi,WANG Nan,et al. Structure form and evolution characteristics of collapsed roof in goaf[J]. Journal of China Coal Society,2021,46(12):3771-3780.

Structure form and evolution characteristics of collapsed roof in goaf

LI Yang,REN Yuqi,WANG Nan,JIN Xiangyang,OU Xiangji,LUO Junbo,MEI Canqi

(School of Energy and Mining Engineering, China University of Mining and Technology-Beijing, Beijing 100083, China)

Abstract:It is of great significance to reveal the structure form of collapsed roof in gob for the upward mining of multiple coal seams with close distance, mining under aquifer and gas distribution in abandoned gob. Based on the detection results of the collapsed roof structure form in goaf after No.12-1 coal seam mining in Kailuan coal group, the characteristics and the formation mechanism of collapsed roof structure form in goaf were studied with theoretical analysis of material mechanics ground pressure and UDEC numerical simulation. A concept of “upright and inverted triangle” collapsed roof structure form was proposed and the dynamic evolution characteristics of collapsed roof structure form was described. The results were as follows:① The detection results of ground penetrating radar showed that the roof stratification was obvious. The upright triangular zone and inverted triangular zone of the collapsed roof appeared alternately in the goaf with the characteristics of orderly stacking of lower collapsing rock blocks and disorderly stacking of upper collapsing rock blocks in the upright triangular area, and disorderly stacking of lower collapsing blocks and orderly stacking of upper collapsing blocks in inverted triangular area. ② The reason for the collapsed roof structure form in goaf was the periodic weighting of roof strata and re-breaking of key block of roof strata in goaf. When the ultimate bending strength of the key block of roof strata in goaf was less than the maximum bending moment, the re-breaking of the key block of roof strata occurred at different positions. The re-breaking position of key rock blocks of main roof and overlying strata was different. Because of the difference of the re-breaking position of different key rock blocks of roof strata, the collapsed roof was divided into two parts, which was the upright triangular area and the inverted triangular area. ③ The concept of “upright and inverted triangle” collapsed roof structure form was proposed. When the stratification of roof strata is obvious, the key rock blocks of main roof and overlying strata re-break under the action of maximum bending moment after the periodic weighting of panel face. The periodic weighting and re-breaking positions present the collapsed roof in the gob as an upright triangular zone and an inverted triangular zone, forming the “upright and inverted triangle” collapsed roof structure form. ④ The “upright and inverted triangle” collapsed roof structure form was also obviously founded by the UDEC software, which is approximately consistent with field detection and theoretical analysis.

Key words:goaf;collapsed roof;evolution characteristics;overlying strata;upright and inverted triangle

目前,由于原始地层赋存的复杂性与特殊性,我国一些主要生产矿井存在近距离煤层上行开采[1-4]和含水层下开采[5-7]的特殊情况。同时,在系统研究废弃矿井煤层气资源开发利用的大背景下,废弃采空区内煤层气资源空间分布研究也成为重点[8-13]。因此煤层开采后采空区垮落顶板的形态特征研究成为当下采空区研究的重点内容。

近年来,国内外学者对于工作面顶板矿压控制理论等方面进行了深入研究,并取得了丰富的成果。钱鸣高院士等[14-15]建立了采场上覆岩层的“砌体梁”结构模型,并提出采场上覆岩层控制的关键层理论,形成一套统一的采场矿山压力与岩层控制理论。许家林等[16-17]基于岩层控制的“砌体梁”结构模型和关键层理论分别建立了不同地质赋存、开采条件下的关键层破断块体的“砌体梁”结构模型,进一步丰富了采场覆岩结构理论。冯国瑞等[18-20]建立了遗煤资源上行开采层间岩层关键层的面接触块体结构,并基于此提出层间岩层的“块体梁-半拱”控制结构。朱涛等[21]建立了近距离煤层开采下煤层的“散体-块体”顶板结构。黄庆享等[22]建立了浅埋煤层采场基本顶周期来压的“短砌体梁”和“台阶岩梁”结构模型。闫少宏等[23]提出大采高采场顶板的“短悬臂梁-铰接岩梁”的结构。

而在煤层开采覆岩空间结构的研究上,刘天泉院士[24]提出了工作面开采覆岩的“横三区、竖三带”的分区,即采空区纵向上分别为垮落带、断裂带和弯曲下沉带,而在横向上分别为煤壁支承影响区、离层区和重新压实区。进一步地,又将垮落带分为下部不规则垮落带与上部规则垮落带。郭广礼等[25]探究了老矿井采空区破裂岩体结构,分为较完整层状结构、块列层状结构、碎裂结构及散体结构,破裂岩体结构对老采空区的变形及空隙分布存在影响。FENG等[26]探究了不同推进距离与采空区积气空间之间的关系。左建平等[27-29]提出了厚松散层覆岩移动的“类双曲线”模型,将采场岩层控制与地表沉降有机相联,更加精确地反映了基岩于松散层的整体移动与变形。SHI等[30]提出一种煤层开采上覆岩层不对称“W”型裂隙拱概念模型。

总体而言,国内外在对采空区垮落顶板的形态特征研究并不多,且现有研究多采用在相似模拟或数值模拟的方法揭示采空区垮落顶板形态[31-37],这与现场实际情况存在差异。针对此,笔者从现场实际探测结果出发,详细分析了采空区垮落顶板形态的特征,结合理论分析的方法,探究了采空区“正-倒置三角形”垮落顶板形态的形成机理,描述了“正-倒置三角形”垮落顶板形态的动态演化特征,给出了“正-倒置三角形”垮落顶板形态概念,并运用UDEC模拟一定程度上验证了“正-倒置三角形”垮落顶板形态的存在性与合理性,进一步丰富了煤层开采覆岩空间结构。

1 垮落顶板形态现场探测

1.1 探测区概况

本次探测区域为开滦集团钱家营煤矿12-1号煤层采空区。开滦集团钱家营矿井田煤系地层约80 m范围内,由上至下依次分布着5,7,8,9和12-1号5层可采煤层,各煤层的平均厚度依次为1.2,4.0,1.4,1.9和3.4 m,且煤层分布具有层间距小、变质程度高和赋存条件复杂的特点。由于产能、煤质及各煤层之间开采的相互影响等问题,逐步形成了“多层采空区下煤层上行开采”结合“先下行后顶层”的高效协调开采模式,即开采顺序依次为7号煤层→8号煤层→12-1号煤层→9号煤层→5号煤层[38],其中5号煤层与9号煤层为采空区上行开采方式。井田范围内煤层赋存情况如图1所示。

图1 煤层赋存情况

Fig.1 Occurrence of coal seam

钱家营矿12-1号煤层的平均埋深为535.9 m,平均倾角约为10°,与上覆9号煤层之间距离约为27 m,煤系地层范围内煤岩岩性与厚度如图2所示。

图2 煤岩层柱状图

Fig.2 Column condition of coal and rock mass

1.2 探测方法

ZTR12系列地质雷达通过屏蔽天线发射1 MHz~2.5 GHz的高频电磁波,由接收天线采集相应的信号并由计算机存储和显示。ZTR12系列地质雷达的天线中心频率为100 MHz、单个天线的长度小于1 m、步进小于2×10-12 s、输出信号为10×10-9 s、电压90 V,另外当产生100 MHz的发射脉冲时,脉冲与纹波的幅值比明显增加到30 dB以上,使得地质雷达的有效填图深度保持在9号煤层以下30 m岩层范围内,因此可以完全适用于矿井井下环境[39]

采用ZTR12系列地质雷达设备在9号煤层1692回风巷对下伏12-1号煤层1622采空区顶板岩层进行探测。地质雷达探测的深度为30 m,探测距离为70 m,探测区域如图3所示,地质雷达系统测试参数见表1。

图3 地质雷达探测位置示意

Fig.3 Location of GPR detection

表1 地质雷达系统测试参数

Table 1 Parameters of GPR detection system

位置天线频率/MHz采样时窗/s采样点数量1692回风巷100300×10-91 024

1.3 探测结果

通过在对采集到的地质雷达探测数据进行零点标定、去噪、滤波、增益等手段,处理得到地质雷达探测图像,可以清晰地反映12-1号煤层开采后顶板岩层的垮落情况,如图4所示。由图4可知,在距离探测起点30~70 m内出现4处分布规律的强、弱反射区交替出现的现象,即在60~70 m位置出现一处弱反射区域,且其相位与周围介质反射图像基本一致,因此黑色虚线区域内为一轻微破坏的区域,该区域顶板的破碎程度较小,完整性较高;而在50~60 m位置出现一处强反射区域,且其相位与周围介质反射图像也存在差异,因此白色虚线区域为一严重破坏的区域,可认为该区域顶板的破碎程度较大,完整性较差,这一规律现象随工作面推进方向循环发生,与现场窥视结果也基本吻合[40]。基于上述探测结果,反演得到12-1号煤层采空区顶板垮落形态,如图5所示。

图4 地质雷达探测数据结果处理示意

Fig.4 Processed detection results of GPR data

图5 煤层开采顶板垮落形态地质雷达反演

Fig.5 Inversion of breaking structure form with geological radar of coal seammining

由图5可知,12-1号煤层开采后直接顶(炭质泥岩)与基本顶(粉砂岩)裂隙较为发育,且裂隙发育已经波及泥岩岩层,即9号煤层基本底岩层,但9号煤层本身受到12-1号煤层开采影响较小。同时,在12-1号煤层采空区内垮落顶板的形态呈现出一定的规律性特征,即采空区垮落顶板分为4个区域,依次为区域Ⅰ、区域Ⅱ、区域Ⅲ和区域Ⅳ。其中,区域Ⅰ和区域Ⅲ对应轻微破坏区域,区域Ⅱ和区域Ⅳ对应严重破坏区域。在区域Ⅰ和区域Ⅲ内,下位岩块破断程度较为剧烈,破碎程度较大,而上位岩块破断程度缓和,破碎程度较小,形成倒置三角形的垮落顶板区域;而在区域Ⅱ和区域Ⅳ内,下位岩块破碎程度较为缓和,破碎程度较小,而上位岩块破碎程度较为剧烈,破碎程度较大,形成正置三角形的垮落顶板区域。

2 垮落顶板形态形成机理

为进一步探究采空区内正置三角形区域与倒置三角形区域垮落顶板形态的形成机理,本节采用材料力学与矿山压力相关知识进行进一步的分析。

在12-1号煤层工作面基本顶发生初次来压之后,随着工作面的不断推进,工作面基本顶岩层会发生周期来压。依据材料力学[41]相关知识,并结合矿山压力中对顶板来压的分析,1622工作面顶板的周期来压步距L计算公式为

(1)

式中,h为基本顶的厚度,m;RT为基本顶岩层的极限抗拉强度,MPa;q为基本顶受到的载荷,可由顶板载荷公式计算可得。

计算可知,上覆岩层中厚度为7 m的粉砂岩层作为关键层对12-1号煤层的开采起控制作用(表2)。因此,关键层下岩层受到的载荷及其极限破断长度,见表3。通过对12-1号煤层开采液压支架载荷的现场实测结果可知,基本顶周期来压步距为13~17 m,如图6所示,也就是说,第3层粉砂岩基本顶岩层破断长度的计算结果与现场实测结果基本吻合。

表2 12-1号煤层上覆岩层参数

Table 2 Parameters of overlying strata above No. 12-1 coal seam

序号(由下至上)名称岩性容重γi/(MN·m-3)层厚hi/m弹性模量Ei/GPa抗拉强度RTi/MPa12直接顶炭质泥岩0.0230.6152.00炭质泥岩0.0232.6152.003基本顶粉砂岩0.0257.3253.504567上覆岩层粉砂岩0.0253.3253.5011号煤0.0140.95 0.03泥岩0.0264.8111.80粉砂岩0.0257.0253.50

图6 工作面液压支架载荷实测结果

Fig.6 Field measurement results of support load of shield in panel

表3 基本顶及其控制的上覆岩层所受载荷及破断长度

Table 3 Load and breaking length of main roof and overlying strata controlled by main roof

序号岩性所受载荷/kPa破断长度/m3粉砂岩330.413.44粉砂岩94.711.3511号煤12.60.76泥岩124.88.6

2.1 采空区基本顶关键岩块的破断

对采空区基本顶关键岩块的受力进行分析,探究基本顶关键岩块的再次破断的发生位置,如图7所示。

图7 基本顶关键岩块受力分析

Fig.7 Stress analysis of key block of main roof

(2)

其中,Fy为岩块受到的垂直作用力,kN;FBB点受到的垂直作用力,kN;FAA点受到的垂直作用力,kN;l1为基本顶关键岩块的长度,取13 m;α为基本顶关键岩块的回转角度,(°),根据钱家营矿12-1号煤层地质条件及现场实测数据,得α=23°;MBB点受到的弯矩,kN·m;TAA点受到的水平作用力,kN;TBB点受到的水平作用力,kN;f为摩擦因数。

由此可得

(3)

故基本顶关键岩块上任意一点的弯矩为


TAxsin α (0≤xl1)

(4)

(5)

其中,x为基本顶关键岩块上任意一点与A点的距离,m。由此得到基本顶关键岩块上任意一点的弯矩,如图8所示。由图8可知,岩块在x=7.5 m处会产生最大弯矩,即基本顶关键岩块会在0.58 l1处发生再次破断。

图8 基本顶关键岩块上任意一点弯矩

Fig.8 Bending moment at any point on key block of main roof

2.2 采空区上覆岩层关键岩块的破断

对采空区基本顶控制的上覆岩层关键岩块的受力进行分析,探究上覆岩层关键岩块的再次破断的发生位置,如图9所示。

图9 上覆岩层关键岩块受力分析

Fig.9 Stress analysis of key block of overlying strata

(6)

式中,FCC点受到的垂直作用力,kN;FDD点受到的垂直作用力,kN;MDD点受到的弯矩,kN·m;l2n为第n层上覆岩层关键岩块长度,m;qn为上覆岩层第n层关键岩块均布载荷,kPa。

故上覆岩层关键岩块上任意一点的弯矩为

(7)

其中,xn为上覆岩层关键岩块上任意一点与C点的距离,m。由此得到基本顶控制的上覆岩层关键岩块上任意一点的弯矩,如图10所示。

图10 上覆岩层关键岩块上任意一点弯矩

Fig.10 Bending moment at any point on key block of overlying strata

由图10可知,岩块分别在x4=5.5 m,x5=0.35 m,x6=4.5 m处会产生最大弯矩,即上覆岩层关键岩块会在0.5l2n处发生再次破断。

通过上述分析可知,在工作面基本顶经历“稳定—失稳—再稳定”的周期破断过程后,采空区内已经垮落的基本顶与其控制的上覆岩层关键岩块在最大弯矩的作用下均会发生岩块的再次破断,但基本顶与其控制的上覆岩层发生再次破断的位置有所不同,且煤岩层赋存情况与破断模式具有一定的层理性,因此将各岩层周期破断与再次破断位置规律连接,则在采空区呈现出正置三角形区域与倒置三角形区域的垮落顶板形态,如图11所示,图中蓝色实线为周期破断位置,图中红色实线为采空区垮落基本顶与其控制的上覆岩层关键岩块再次破断所在位置。因此可认为,基本顶与其控制的上覆岩层的周期破断与采空区已垮落关键岩块的再次破断共同主导了采空区垮落顶板呈现出正置三角形区域与倒置三角形区域的形态。

图11 采空区垮落顶板结构形态示意

Fig.11 Sketch map of structure form of collapsing roof in goaf

3 垮落顶板形态演化特征

在煤层开采过程中,随着工作面的不断推进,基本顶及其上覆岩层在发生“稳定—失稳—再稳定”的周期来压之后,采空区内已经断裂的基本顶与其控制的上覆岩层关键岩块会再次发生如图12所示的演化过程。

随着工作面推进,在基本顶发生周期破断形成新岩块A,其相邻岩块B即成为基本顶关键岩块,此时,由于岩块A与岩块B的回转与反向回转,导致岩块A的前咬合点产生向上运动趋势,使得岩块A的前端点破碎而发生失稳,由此岩块A与岩块B回转形成一体,如图12(a),(b)所示。之后,采空区内的基本顶关键岩块B在最大弯矩的作用下发生再次破断,而其控制的上覆岩层同样在最大弯矩的作用下发生再次破断,如图12(b)所示。但由于基本顶与其上覆岩层的破断位置存在差异,遂出现偏向工作面方向的破断线,如图12(b)中红色实线所示。煤系地层的形成往往是由于地质沉积作用导致的,因此,采空区内垮落顶板发生再次破断时,煤层顶板各岩层会出现明显的分层现象,如图12(c)所示。综上在考虑采空区垮落顶板周期破断与再次破断位置,即在采空区内基本顶与其控制的上覆岩层破断垮落稳定后,会形成正置三角形区域与倒置三角形区域的垮落顶板形态,如图12(d)所示。重复上述过程,采空区内即出现“正-倒置三角形”垮落顶板的形态,如图11所示。

图12 采空区垮落顶板形态演化过程

Fig.12 Evolution process of the structure form of collapsing roof in goaf

基于此,采空区“正-倒置三角形”垮落顶板的形态概念:当煤层顶板岩层分层明显时,在煤层开采工作面周期来压后,采空区内基本顶与其控制的上覆岩层关键岩块在最大弯矩的作用下发生再次破断,周期破断与再次破断位置将采空区垮落顶板规律地分为正置三角形区域与倒置三角形区域,形成“正-倒置三角形”的垮落顶板形态。

4 垮落顶板形态模拟验证

为更直观验证采空区内“正-倒置三角形”垮落顶板的形态的存在性,本节以钱家营矿12-1号煤层开采的实际情况为背景,建立采空区垮落顶板形态研究的UDEC数值模型(图13),模型长×宽为200 m×160 m。采用Mohr-Coulomb本构模型,模型中各岩层力学参数见表4。模型两侧边界为水平位移约束、下部为边界位移约束,模型上部边界施加10.78 MPa的竖向载荷,以模拟未建立的实际地层载荷,模型如图13所示。为避免模拟过程中的边界效应,模型左右两侧各留宽65.5 m煤柱,并采用全部垮落法处理采空区。

表4 数值模型参数

Table 4 Parameters of numerical simulation

序号颜色岩性密度/(kg·m-3)体积模量/GPa剪切模量/GPa内摩擦角/(°)黏聚力/MPa抗拉强度/MPa1上覆岩层2 7009.604.7037.04.005.212中粒砂岩2 5706.253.5731.49.208.503粉砂岩2 5804.502.8037.03.503.504细砂岩2 2007.356.6340.03.044.005炭质泥岩2 6003.501.8033.03.002.006泥岩2 4202.501.7229.52.111.807煤层1 4002.080.5420.01.200.038下伏岩层2 7009.604.7037.04.005.21

图13 数值模型示意

Fig.13 Numerical model

在12-1号煤层工作面采空区内,会出现正置三角形区域与倒置三角形区域交替出现的垮落顶板形态如图14所示。在顶板的周期性垮落后,顶板岩层由下至上依次发生破断,当工作面继续向前推进时,采空区内基本顶关键岩块在最大弯矩作用下会发生再次破断,而基本顶上控制的上覆岩层同样在最大弯矩的作用下发生再次破断,破断岩块会分别发生回转与反向回转,将周期破断与再次破断位置规律连接,垮落顶板岩块在采空区内形成正置三角形区域与倒置三角区域的垮落顶板形态。当工作面继续推进至下一次周期性垮落时,同样垮落顶板岩块在采空区内又形成正置三角形区域与倒置三角形区域的形态,循环往复,则形成采空区垮落顶板“正-倒置三角形”交替出现的形态。

图14 12-1号煤层采空区垮落顶板形态

Fig.14 Structure form of collapsing roof in goaf of No.12-1 coal seam

其中,正置三角形区域具有下位垮落岩块堆砌整齐、上位垮落岩块堆砌杂乱的特点,而倒置三角形区域具有下位下沉岩块完整性差、上位下沉岩块完整性好的特点。因此,随着12-1号煤层工作面的不断推进,采空区内会依次出现正置三角形区域和倒置三角形区域交替出现的垮落顶板形态。

5 结 论

(1)采用地质雷达对12-1号煤层采空区进行探测,结果显示12-1号煤层顶板分层明显,且12-1号煤层采空区呈现正置三角形区域与倒置三角形区域规律出现的垮落顶板形态,且经地质雷达反演图可知正置三角形垮落顶板区域内下位垮落岩块堆砌整齐、上位垮落岩块堆砌杂乱,而倒置三角形垮落顶板区域内下位下沉岩块堆砌杂乱、上位下沉岩块完整性好。

(2)12-1号煤层采空区呈现的正置三角形与倒置三角形垮落顶板形态的规律性现象是由于顶板岩层的周期破断与采空区已垮落关键岩块的再次破断引起的。基本顶与其控制的上覆岩层关键块体发生再次破断的原因是关键岩块受到的最大弯矩大于其极限抗弯强度。同时,基本顶岩块与上覆岩层发生再次破断的位置不同。基本顶关键岩块发生再次破断的位置为0.58l1,基本顶控制的上覆岩层关键岩块发生再次破断的位置为0.5l2n

(3)探究了采空区垮落顶板形态的演化过程,提出了采空区垮落顶板“正-倒置三角形”形态的概念,当基本顶及其控制的上覆岩层分层明显时,在煤层开采工作面周期来压后,基本顶及其控制的上覆岩层关键岩块在最大弯矩的作用下发生再次破断,周期破断与再次破断位置将采空区垮落顶板呈现正置三角形与倒置三角形区域,形成“正-倒置三角形”的垮落顶板形态。

(4)运用UDEC模拟一定程度上验证了“正-倒置三角形”垮落顶板形态的存在性与合理性。

参考文献(References):

[1] 冯国瑞. 煤矿残采区上行开采基础理论与实践[M]. 北京:煤炭工业出版社,2010.

[2] 张金海. 钱家营矿采煤工作面上行开采可行性分析[J]. 煤炭科学技术,2009,37(1):62-64,29.

ZHANG Jinhai. Feasibility analysis on upward mining for coal mining face of Qianjiaying mine [J]. Coal Science and Technology, 2009, 37(1):62-64,29.

[3] 汪理全,李中颃. 煤层群上行开采技术[M]. 北京:煤炭工业出版社,1995.

[4] 韩军,张宏伟,张普田,等. 开滦矿区近距离煤层群上行开采可行性研究[J]. 煤炭科学技术,2011,39(10):14-17.

HAN Jun,ZHANG Hongwei,ZHANG Putian,et al. Feasibility study on upward mining in seam group with short distance to each other in kailuan mining area [J]. Coal Science and Technology,2011,39(10):14-17.

[5] 许家林,朱卫兵,王晓振. 松散承压含水层下采煤突水机理与防治研究[J]. 采矿与安全工程学报,2011,28(3):333-339.

XU Jialin,ZHU Weibin,WANG Xiaozhen. Study on water-inrush mechanism and prevention during coal mining under unconsolidated confined aquifer [J]. Journal of Mining and Safety Engineering,2011,28(3):333-339.

[6] 许家林,蔡东,傅昆岚. 邻近松散承压含水层开采工作面压架机理与防治[J]. 煤炭学报,2007,32(12):1239-1243.

XU Jialin,CAI Dong,FU Kunlan. Mechanism of supports crushing accident and its prevention measures during coal mining near unconsolidated confined aquifer [J]. Journal of China Coal Society,2017,32(12):1239-1243.

[7] 郭小铭,董书宁. 深埋煤层开采顶板基岩含水层渗流规律及保水技术[J]. 煤炭学报,2019,44(3):805-812.

GUO Xiaoming,DONG Shuning. Seepage law of bedrock aquifer and water-preserved mining technology in deep coal seam mining [J]. Journal of China Coal Society,2019,44(3):805-812.

[8] 钱鸣高,许家林. 覆岩采动裂隙分布的“O”形圈特征研究[J]. 煤炭学报,1998,23(5):3-5.

QIAN Minggao,XU Jialin. Study on the “O-shape” circle distribution characteristics of mining-induced fractures in the overlying strata [J]. Journal of China Coal Society,1998,23(5):3-5.

[9] 邓喀中,周鸣,谭志祥. 采动岩体破裂规律的试验研究[J]. 中国矿业大学学报,1998,27(3):261-264.

DENG Kazhong,ZHOU Ming,TAN Zhixiang. Study on laws of rock mass breaking induced by mining [J]. Journal of China University of Mining and Technology,1998,27(3):261-264.

[10] 王家臣,李杨译. 做得好——德国应对后采矿的挑战[M]. 北京:科学出版社,2020.

[11] 王家臣. 基于采动岩层控制的煤炭科学开采[J]. 采矿与岩层控制工程学报,2019,1(2):40-47.

WANG Jiachen. Sustainable coal mining based on mining ground control [J]. Journal of Mining and Strata Control Engineering,2019,1(2):40-47.

[12] 王家臣,杨胜利,杨宝贵,等. 深井超长工作面基本顶分区破断模型与支架阻力分布特征[J]. 煤炭学报,2019,44(1):54-63.

WANG Jiachen,YANG Shengli,YANG Baogui,et al. Roof sub-regional fracturing and support resistance distribution in deep longwall face with ultra-large length [J]. Journal of China Coal Society,2019,44(1):54-63.

[13] 杨胜利,王家臣,李良晖. 基于中厚板理论的关键岩层变形及破断特征研究[J]. 煤炭学报,2020,45(8):2718-2727.

YANG Shengli,WANG Jiachen,LI Lianghui. Deformation and fracture characteristics of key strata based on the medium thick plate theory [J]. Journal of China Coal Society,2020,45(8):2718-2727.

[14] 钱鸣高,缪协兴,何富连. 采场“砌体梁”结构的关键块分析[J]. 煤炭学报,1994,19(6):557-563.

QIAN Minggao,MIAO Xiexing,HE Fulian. Analysis of key block in the structure of voussoir beam in longwall mining [J]. Journal of China Coal Society,1994,19(6):557-563.

[15] 钱鸣高,缪协兴,许家林. 岩层控制中的关键层理论研究[J]. 煤炭学报,1996,21(3):2-7.

QIAN Minggao,MIAO Xiexing,XU Jialin. Theoretical study of key stratum in ground control [J]. Journal of China Coal Society,1996,21(3):2-7.

[16] 许家林,钱鸣高. 覆岩关键层位置的判别方法[J]. 中国矿业大学学报,2000,29(5):21-25.

XU Jialin,QIAN Minggao. Method to distinguish key strata in overburden [J]. Journal of China University of Mining and Technology,2000,29(5):21-25.

[17] 许家林,钱鸣高. 关键层运动对覆岩及地表移动影响的研究[J]. 煤炭学报,2000,25(22):122-126.

XU Jialin,QIAN Minggao. Study on the influence of key strata movement on subsidence [J]. Journal of China Coal Society,2000,25(2):122-126.

[18] 冯国瑞,王鲜霞,康立勋. 采场上覆岩层面接触块体结构的力学机理分析[J]. 煤炭学报,2008,33(1):33-37.

FENG Guorui,WANG Xianxia,KANG Lixun. Analysis on the mechanism of the face-contacted blocks structure in overlying strata above the longwall face [J]. Journal of China Coal Society,2008,33(1):33-37.

[19] 冯国瑞,闫旭,王鲜霞,等. 上行开采层间岩层控制的关键位置判定[J]. 岩石力学与工程学报,2009,28(S2):3721-3726.

FENG Guorui,YAN Xu,WANG Xianxia,et al. Determination of key position of strata controlling in rocks between coal seams for upward mining [J]. Chinese Journal of Rock Mechanics and Engineering,2009,28(S2):3721-3726.

[20] 冯国瑞,任亚峰,王鲜霞,等. 白家庄煤矿垮落法残采区上行开采相似模拟实验研究[J]. 煤炭学报,2011,36(4):544-550.

FENG Guorui,REN Yafeng,WANG Xianxia,et al. Experimental study on the upward mining of the left-over coal above gob area mined with caving method in Baijiazhuang Coal Mine [J]. Journal of China Coal Society,2011,36(4):544-550.

[21] 朱涛,张百胜,冯国瑞,等. 极近距离煤层下层煤采场顶板结构与控制[J]. 煤炭学报,2010,35(2):190-193.

ZHU Tao,ZHANG Baisheng,FENG Guorui,et al. Roof structure and control in the lower seam mining field in the ultra-close multiple seams [J]. Journal of China Coal Society,2010,35(2):190-193.

[22] 黄庆享,钱鸣高,石平五. 浅埋煤层采场老顶周期来压的结构分析[J]. 煤炭学报,1999,24(6):581-585.

HUANG Qingxiang,QIAN Minggao,SHI Pingwu. Structural analysis of main roof stability during periodic weighting in longwall face [J]. Journal of China Coal Society,1999,24(6):581-585.

[23] 闫少宏,尹希文,许红杰,等. 大采高综采顶板短悬臂梁-铰接岩梁结构与支架工作阻力的确定[J]. 煤炭学报,2011,36(11):1816-1820.

YAN Shaohong,YIN Xiwen,XU Hongjie,et al. Roof structure of short cantilever-articulatedrock 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.

[24] 刘天泉. 矿山岩体采动影响与控制工程学及其应用[J]. 煤炭学报,1995,20(1):1-5.

LIU Tianquan. Influence of mining activities on minerockmass and control engineer [J]. Journal of China Coal Society,1995,20(1):1-5.

[25] 郭广礼,缪协兴,张振南. 老采空区破裂岩体变形性质研究[J]. 科学技术与工程,2002,2(5):44-47.

GUO Guangli,MIAO Xiexing,ZHANG Zhennan. Research on rupturedrock mass deformation characteristics of longwall goafs [J]. Science Technology and Engineering,2002,2(5):44-47.

[26] FENG Guorui,LI Zhen,HU Shengyong,et al. Distribution of gob empty space for methane drainage during the longwall mining:A case study [J]. Journal of Natural Gas Science and Engineering,2018,60:112-124.

[27] 左建平,孙运江,钱鸣高. 厚松散层覆岩移动机理及“类双曲线”模型[J]. 煤炭学报,2017,42(6):1372-1379.

ZUO Jianping,SUN Yunjiang,QIAN Minggao. Movement mechanism and analogous hyperbolamodel of overlying strata with thick alluvium [J]. Journal of China Coal Society,2017,42(6):1372-1379.

[28] 左建平,孙运江,王金涛,等. 充分采动覆岩“类双曲线”破坏移动机理及模拟分析[J]. 采矿与安全工程学报,2018,35(1):71-77.

ZUO Jianping,SUN Yunjiang,WANG Jintao,et al. Mechanical and numerical analysis of “analogous hyperbola” movement of overlying strata after full mining extraction [J]. Journal of Mining and Safety Engineering,2018,35(1):71-77.

[29] 左建平,吴根水,孙运江,等. 岩层移动内外“类双曲线”整体模型研究[J]. 煤炭学报,2021,46(2):333-343.

ZUO Jianping,WU Genshui,SUN Yunjiang,et al. Investigation on the inner and outer analogous hyperbolamodel (AHM) of strata movement [J]. Journal of China Coal Society,2021,46(2):333-343.

[30] SHI Longqing,XU Dongjing,WANG Ying,et al. A novel conceptualmodel of fracture evolution patterns in the overlying strata during horizontal coal seam mining [J]. Arabian Journal of Geosciences,2019,12(10):326.

[31] 张勇,刘传安,张西斌,等. 煤层群上行开采对上覆煤层运移的影响[J]. 煤炭学报,2011,36(12):1990-1995.

ZHANG Yong,LIU Chuanan,ZHANG Xibin,et al. The influence of ascending mining on the movement character of overlying coal seam in coal seams group [J]. Journal of China Coal Society,2011,36(12):1990-1995.

[32] 何尚森,谢生荣,宋宝华,等. 近距离下煤层损伤基本顶破断规律及稳定性分析[J]. 煤炭学报. 2016. 41(10):2596-2605.

HE Shangsen,XIE Shengrong,SONG Baohua,et al. Breaking laws and stability analysis of damage main roof in close distance hypogynous seams [J]. Journal of China Coal Society,2016,41(10):2596-2605.

[33] 侯运炳,何尚森,谢生荣. 近距离煤层层间基本顶损伤及破断规律研究[J]. 岩土力学,2017,38(10):2989-2999,3008.

HOU Yunbing,HE Shangsen,XIE Shengrong. Damage and rupture laws of main roof between coal seams with a close distance [J]. Rock and Soil Mechanics,2017,38(10):2989-2999,3008.

[34] 蒋金泉,孙春江,尹增德,等. 深井高应力难采煤层上行卸压开采的研究与实践[J]. 煤炭学报,2004,29(1):1-6.

JIANG Jinquan,SUN Chunjiang,YIN Zengde,et al. Study and practice of bottom slicing relieving mining in condition of trouble coal seams with high ground stress under deep coal mine [J]. Journal of China Coal Society,2004,29(1):1-6.

[35] 张向阳,任启寒,涂敏,等. 潘一东矿近距离煤层上行开采围岩裂隙演化规律模拟研究[J]. 采矿与安全工程学报,2016,33(2):191-198.

ZHANG Xiangyang,REN Qihan,TU Min,et al. Fracture evolution characteristics of surroundingrock during ascending mining of closer coal seam in Panyidong coal mine[J]. Journal of Mining & Safety Engineering,2016,33(2):191-198.

[36] 李树刚,丁洋,安朝峰,等. 近距离煤层重复采动覆岩裂隙形态及其演化规律实验研究[J]. 采矿与安全工程学报,2016,33(5):904-910.

LI Shugang,DING Yang,AN Zhaofeng,et al. Experimental research on the shape and dynamic evolution of repeated mining-induced fractures in short-distance coal seams[J]. Journal of Mining & Safety Engineering,2016,33(5):904-910.

[37] 王朋飞,赵景礼,王志强,等. 非充分采动采空区与煤岩柱(体)耦合作用机制及应用[J]. 岩石力学与工程学报,2017,36(5):1185-1200.

WANG Pengfei,ZHAO Jingli,WANG Zhiqiang,et al. Mechanism of gob-pillar interaction for subcritical panels and its application [J]. Chinese Journal of Rock Mechanics and Engineering,2017,36(5):1185-1200.

[38] 李杨,雷明星,郑庆学,等. 近距离“薄-中-厚”交错分布煤层群上行协调开采定量判别研究[J]. 煤炭学报,2019,44(S2):410-418.

LI Yang,LEI Mingxing,ZHENG Qingxue,et al. Quantitative criterion on coordinated ascending mining in close multiple “thin-medium-thick” coal seams[J]. Journal of China Coal Society,2019,44(S2):410-418.

[39] LI Yang,WANG Jiachen,CHEN Yiding,et al. Overlying strata movement with ground penetrating radar detection in close-multiple coal seams mining[J]. International Journal of Distributed Sensor Networks,2019,15(8):1550147719869852.

[40] LI Yang,REN Yuqi,SYD S Peng,et al. Measurement of over-burden failure zones in close-multiple coal seams mining [J]. International Journal of Mining Science and Technology,2021,31(1):43-50.

[41] 钱鸣高,石平五,许家林. 矿山压力与岩层控制[M]. 徐州:中国矿业大学出版社,2010.

[42] 孙训方,方孝淑,关来泰. 材料力学(Ⅰ)[M]. 北京:高等教育出版社,2009.