摘要:目前有关巷道底鼓的研究与实践主要探讨巷道底板的变形机理及控制技术,对沿空留巷切顶卸压前后底板力学分析不全面。针对该问题,基于煤体分区破坏特征构建了切顶前后巷道围岩和底板力学模型,分析实体煤、巷旁支护及采空区对底板的作用,获得切顶前后巷道底鼓解析解,得出巷旁煤帮弹塑性区、巷道支护体及顶板下沉区底板所受载荷共同影响巷道底鼓量大小。采用数值模拟验证切顶卸压前后沿空留巷围岩破坏特征、应力分布及底鼓量变化,结果表明:切顶卸压技术可有效缩小巷道实体煤侧面及顶部的破坏区域,维持巷道围岩结构稳定;巷道底板最大应力、巷旁支护阻力、巷道底鼓量均下降,平均降幅分别为25.78%,56.14%,54.07%。现场应用结果表明,厚硬顶板沿空留巷底鼓量由709.345 1 mm 降至320.965 8 mm,切顶卸压技术可以优化巷道围岩应力结构,抑制巷道底鼓,有效改善底板破坏情况。
关键词:底板变形控制;沿空留巷;切顶卸压;力学分析;围岩破坏特征;应力分布;底鼓
中图分类号:TD322.5 文献标志码:A
0 引言
在厚硬顶板沿空留巷条件下,受掘进和采动影响,巷道底板变形,出现底鼓。随着工作面推进,基本顶在采空区上方形成悬顶[1]。一旦发生顶板旋转下沉,会对煤岩体和巷旁支护造成冲击挤压[2],加剧巷道底板变形[3]。底鼓现象造成巷道断面减小,阻碍设备运输、人员行走,影响工作面正常回采,增加大量维修工作和维护费用。切顶卸压技术可以在巷道采空区一侧形成预裂切缝面[4],切断采空区和巷道顶板之间的联系,改善巷道围岩结构,减少底板变形量。因此,根据切顶卸压前后的围岩结构变化进一步研究巷道底鼓的变形机理具有一定意义。
针对沿空留巷切顶卸压的底板变形问题,学者们进行了大量研究。康志鹏等[5]以厚煤层软底留巷为背景,提出对底板进行合适的让压,以实现巷道柔性支护。于光远等[6]提出通过“切顶卸压+柔性让压+补强锚索控顶+双控锚杆控帮”的方式控制底鼓,取得了较好的效果。何满潮等[7]针对深部高应力复合破碎顶板切顶留巷,设计了一种新型可缩U 型钢挡矸结构,可有效控制巷道底板变形。张幼振[8]提出了巷道底鼓机械化治理的技术思路和机械化施工系统的配置原则,研发了巷道修复机、底板锚固钻机等专用主导机械设备,可有效治理巷道底鼓。申斌学等[9]提出了切顶卸压柔模墙支护沿空留巷技术,在切顶后使用该技术可以对巷道顶板进行强力支护,抑制巷道底鼓。张宇旭等[10]通过试验证明了底板注浆与锚杆锚索联合支护可以改善底板应力状态,控制底板位移。以上有关巷道底鼓研究与实践主要探讨巷道底板的变形机理及控制技术,对沿空留巷切顶卸压前后底板力学分析不全面。针对上述问题,本文建立了巷道围岩和底板力学模型,分析各区域对巷道底板的作用,并推导了沿空留巷切顶前后巷道底鼓变形量计算公式。
1 工程背景
以淮南某矿沿空留巷工作面为研究对象,工作面长210 m,走向长1 784 m,1 号煤厚3.5~8.5 m,平均厚度为6.5 m。煤层倾角为6~8°,平均值为7°。工作面煤层直接顶主要为泥岩,平均厚1 m;基本顶主要为中细砂岩,平均厚13.6 m;直接底主要为泥岩,平均厚0.5 m;基本底主要为砂质泥岩,平均厚3.0 m。工作面布置及岩性特征如图1 所示。
2 巷道底鼓力学分析
2.1 沿空留巷围岩结构
工作面开挖后,因直接顶极薄且强度较低,采空区直接顶垮落,基本顶在实体煤体侧断裂。巷旁支护体在顶板载荷作用下产生变形,压缩下沉。此时顶板处于给定变形状态[11],如图2 所示,巷旁支护体主要承受直接顶及悬露基本顶岩块B 的重力。上覆岩层破断岩块A, B, C 三者相互铰接形成承压结构[12],承载上覆岩层压力并控制下方巷道的围岩稳定。岩块C 下沉后,在上覆岩层载荷作用下与矸石接触发生挤压,并逐渐被压实[13]。
2.2 围岩结构力学分析
在给定变形条件下,侧向实体煤承担围岩应力峰值,巷旁支护体无法有效控制基本顶,岩块B 发生旋转下沉。对沿空留巷顶板力学模型作如下假设:① 基本顶于实体煤侧弹塑性交界处断裂,绕直接顶向采空区侧旋转倾斜。② 采空区矸石对岩块C 产生支撑力,对岩块B 支撑力为0[14]。③ 忽略直接顶、基本顶及更上位岩层之间的剪切力。④ 基本顶上覆岩层载荷均匀施加在基本顶上。⑤ 巷内支护的变形忽略不计[15]。根据假设条件,建立沿空留巷力学模型,如图3 所示。
基本顶沿实体煤弹塑性区交界处上方点断裂并向采空区旋转,煤体极限平衡区宽度x0[16]及塑性区煤体对顶板的支撑力σ 分别为
式中:hc 为煤层厚度;λ 为煤体侧压系数;φ 为煤层界面内摩擦角;k 为煤帮侧应力集中系数;γ 为岩层平均容重;H 为开采深度;c 为煤层界面黏聚力;px 为煤帮侧向支护阻力。
基本顶沿倾向垮落长度L 与垮落步距l[17]分别为
式中:h 为基本顶岩层厚度;RT 为基本顶抗拉强度;q2 为基本顶载荷;S 为基本顶沿工作面悬露长度。
岩块C 水平方向受力取0,可得
式中:TB 为岩块B 在处所受水平推力;TC 为岩块C 在C′处所受水平推力;ΔSB为岩块B 在处的下沉量。
岩块C 垂直方向受力取0,可得
NB + FG = NB +KGL = NC +q2L (5)
式中:NB ,NC分别为岩块B,C 所受剪切力;FG为采空区矸石对基本顶的支撑力;KG为采空区矸石支撑系数。
B′点力矩取0,可得
式中MB 为岩块B 在处的残余弯矩。
A′ 对于岩块B, 点力矩取0,可得
式中:MA 为岩块A 在处的残余弯矩;M0 为直接顶对基本顶的抗弯弯矩;ft 为巷旁支护阻力;a 为巷旁支护体宽度;b 为巷道宽度;x 为模型中任一点的横坐标;q1 为直接顶载荷。
联立式(3)、式(4)、式(8)、式(9)得
式中ΔSC为岩块C 在C'处的下沉量。
2.3 巷道底板受力分析
顶板垮落后,底板应力重新分布。为了简化计算,假设煤帮弹性区、塑性区的垂直应力呈线性分布,沿工作面走向方向取剖面,建立沿空留巷底板力学模型,如图4 所示。fg 为顶板破断岩块对底板的作用力;xm 为顶板旋转下沉区宽度;fm 为冒落矸石的重力;x1 为弹性压缩区宽度;f0,f1 分别为塑性区、弹性压缩区煤体对底板的作用力。
实体煤塑性区作用在底板的力为
f0 = (k-1)H(x0 - x) =x0 - x0≤x≤0 (11)
实体煤的弹性压缩区和塑性区共同承担了上覆岩层载荷和采动引起的应力重新分布, 因此(x1+x0) (k+1)γH/2 = γHS0/2,S0为工作面长度,则
x1 = S0/(k+1)- x0 (12)
实体煤弹性压缩区作用在底板的力为
f1 = (k-1)γH(x- x1) /(x1 - x0) - x0 - x1≤x≤- x0(13)
直接顶垮落矸石对底板的作用力为
fm = h1γ1 (14)
式中:h1 为直接顶岩层厚度;γ1 为直接顶岩层容重。
上方基本顶岩块C 垮落后,沿倾向逐渐将矸石堆压实。顶板旋转下沉区宽度[18]为
xm = hcot δ+ L (15)
式中δ 为基本顶的破断角[19]。
式中φ2 为基本顶内摩擦角。
顶板破断下沉岩块对底板的作用力为
fg = 2FG [x-(L- x0)] /x2m(17)
2.4 巷道底板变形分析
引用等效载荷概念 H [20],以消除原岩应力对底板变形的影响。建立沿空留巷底板等效载荷分布力学模型,如图5 所示。
在半无限平面边界上,点I 所受应力为原岩应力,设垂直位移为0,将该点设为基点,则分布载荷q(x)在边界上点J(横坐标为xJ)处的垂直位移[21-22]为
式中:UJ 为点J 处的底鼓量;Ed 为巷道底板岩层弹性模量;s 为基点I 到q(x)的距离;ρ 为点J 到q(x)的距离。
将不同分布载荷在J 处引起的垂直位移相加,求得J 处的底鼓量:
3 切顶卸压对沿空留巷围岩结构影响
3.1 切顶卸压围岩结构力学分析
通过预爆破切顶,基本顶沿切顶缝断裂,侧向悬臂梁长度大幅减小,在巷道上方形成短臂梁结构[23]。由于直接顶厚度极薄,垮落后形成的矸石堆不足以充填整个采空区并支撑上覆岩层,所以采空区上方破断基本顶沿切缝大幅下沉并压实下方矸石堆。切顶后沿空留巷围岩结构如图6 所示。
对切顶后沿空留巷顶板力学模型作如下假设:① 切顶角度为90°,简化力学模型。② 采空区侧基本顶沿切缝直接断裂垮落,岩块下沉时旋转角度为0。③ 垮落后的矸石堆积在采空区,对下沉基本顶起到支撑作用[24]。④ 忽略直接顶、基本顶及更上位岩层之间的剪切力。⑤ 基本顶上覆岩层载荷均匀施加在基本顶上。⑥ 巷内支护的变形忽略不计。⑦ 实体煤弹塑性区域及煤层底板所受应力在切顶前后的变化较覆岩整体的运移程度小,故忽略不计。
根据假设,构建切顶后沿空留巷力学模型,如图7 所示。T′B为岩块A 在处所受水平推力;T′C岩块B 在处所受水平推力;N′B,N′C分别为岩块A,B 受到的剪切力;M′A为岩块A 在处产生的弯矩;M′B为岩块B 在处的残余弯矩; f ′t为切顶后巷旁支护阻力;F′G为切顶后采空区矸石对基本顶支撑力。
根据静力平衡可得
式中ΔS ′C为岩块B 在处的下沉量。
3.2 切顶卸压巷道底板变形分析
切顶后沿空留巷底板等效载荷分布力学模型如图8 所示。
切顶后顶板破断岩块对底板作用力为
f ′g = F′G = K′GL (24)
式中K' G为切顶后采空区矸石支撑系数。
底板上M 点的底鼓量为
4 切顶卸压底鼓数值模拟研究
4.1 模型的建立
采用Mohr-Coulomb 本构模型,从现场获取巷道围岩岩体试样,按照相关要求制备成相应的标准试件。各岩层物理力学参数见表1。
根据地质资料建立立体模型,模拟在上一工作面回采完毕后,接替回采1615A 工作面时,切顶卸压对巷道底板变形的控制效果。
4.2 沿空留巷巷道演化特征
4.2.1 切顶卸压巷道围岩破坏分布特征及应力分析
沿走向取剖面,切顶卸压前后工作面处破坏分布特征如图9 所示。可看出工作面回采时采空区上方形成弹塑性破坏区域,且随工作面的推进逐渐向高位岩层发育,由破坏形态可知采空区上覆中位岩层逐渐趋于稳定。巷道实体煤侧及顶部围岩在切顶卸压后,破坏区域显著减小。随着工作面的推进,采用切顶卸压技术可有效缩小巷道实体煤侧面及顶部的破坏区域,维持巷道围岩结构稳定。
沿走向取剖面,切顶卸压前后工作面处垂直应力如图10 所示。可看出:切顶卸压前,巷道实体煤侧出现应力集中现象,应力峰值随着工作面推进由19.22 MPa 增加至23.13 MPa,巷旁支护阻力峰值为5.75 MPa,巷道底板中部应力峰值为9.77 MPa,对巷道围岩稳定造成严重影响;切顶卸压后,应力集中向巷道顶板上方及深部转移,远离巷道围岩,应力峰值由19.92 MPa 增加至28.20 MPa,巷旁支护阻力峰值减少至2.58 MPa,巷道底板中部应力峰值为7.25 MPa,巷道顶板区域所受应力显著减小。采用切顶卸压技术可有效优化巷道围岩应力结构,巷旁支护阻力降幅为56.14%,底板最大应力平均降幅为25.78%。
4.2.2 切顶卸压巷道底鼓分析
对整条巷道底板变形量进行监测,结果如图11所示。可知切顶前,工作面超前处巷道底板变形量在503~917 mm 波动,均值为664 mm;切顶后,工作面超前处巷道底板变形量在204~409 mm波动,均值为305 mm,较切顶前下降54.07%。沿空留巷切顶卸压有效缓解了沿空留巷底鼓。
5 工程应用
预裂切缝钻孔深12 m,切缝角度为15°。采用双向聚能爆破预裂技术,切缝距煤帮100 mm,切缝孔间距为500 mm。特质聚能管外径为42 mm,内径为36.50 mm,管长1 500 mm。炮空底部装药量为3~5 卷,炮孔处为1~2 卷,炸药规格为ϕ35×330 mm/卷。现场炮孔布置如图12 所示。
结合现场数据对巷道底鼓解析解进行验证。基本参数: hc=6.5 m, h=13.6 m, px=0.1 MPa, φ=24°,H=600 m, γ=25 kN/m3, c=2.4 MPa, k=2.6, S0=200 m,Ed=9.0 GPa,η=0.9,KP0=1.10,a=1 m,b=5 m。
其他参数设置: 直接顶及基本顶载荷根据q=γh确定;垮落长度L=20 m;岩梁断裂处残余弯矩为0;ΔSB = ΔSC;实体煤塑性区宽度x0=10.3 m;切顶前后采空区矸石对基本顶支撑力FG5.237 MN,F′G=10.532 MN;切顶前后巷旁支护阻力ft=5.564 5 MN,f ′t=2.194 2 MN。
将上述参数代入式(19)、式(25),可得切顶前后巷道底鼓量UM = -709:345 1 mm,U′M = 320:965 8 mm。
切顶前后沿空留巷底鼓情况如图13 所示。可看出沿空留巷切顶后,巷道底鼓量减小,切顶卸压有效改善巷道底板破坏情况。
6 结论
1) 针对工作面现场情况,运用给定变形理论,构建了沿空留巷顶板岩层力学模型,推导出切顶卸压前后巷旁支护阻力的计算方法。
2) 运用弹性力学知识,结合底板各区域所受载荷,推导出切顶卸压前后巷道底鼓计算公式。分析得出巷旁煤帮弹塑性区、巷道支护体及顶板下沉区底板所受载荷共同影响巷道底鼓量大小。
3) 切顶卸压后,应力集中现象向巷道顶板上方及深部转移,巷道底板最大应力平均降幅为25.78%,巷旁支护阻力平均降幅为56.14%;厚硬顶板沿空留巷底鼓量由709.345 1 mm 降至320.965 8 mm。切顶卸压技术可以优化巷道围岩应力结构,抑制巷道底鼓,有效改善底板破坏情况。
参考文献(References):
[ 1 ]潘立友,唐鹏,周脉来,等. 悬顶结构巷道冲击地压防控研究[J]. 煤炭科学技术,2022,50(4):42-48.
PAN Liyou,TANG Peng,ZHOU Mailai,et al. Researchon prevention and control of rock burst in entry withsuspended roof structure[J]. Coal Science andTechnology,2022,50(4):42-48.
[ 2 ]何满潮,陈[ 2 ] 上元,郭志飚,等. 切顶卸压沿空留巷围岩结构控制及其工程应用[J]. 中国矿业大学学报,2017,46(5):959-969.
HE Manchao, CHEN Shangyuan, GUO Zhibiao, et al.Control of surrounding rock structure for gob-side entryretaining by cutting roof to release pressure and itsengineering application[J]. Journal of China Universityof Mining & Technology,2017,46(5):959-969.
[ 3 ]杨鸿智, 鞠文君, 王东攀, 等. 综放工作面沿空留巷 “错动” 变形破坏机制及其控制[J]. 采矿与安全工程学报,2023,40(2):224-231.
YANG Hongzhi, JU Wenjun, WANG Dongpan, et al.The "relative movement" deformation mechanism andcontrol technique of gob-side entry retaining in top-coalcaving mining[J]. Journal of Mining & SafetyEngineering,2023,40(2):224-231.
[ 4 ]马新根,何满潮,李先章,等. 切顶卸压自动成巷覆岩变形机理及控制对策研究[J]. 中国矿业大学学报,2019,48(3):474-483.
MA Xingen, HE Manchao, LI Xianzhang, et al.Deformation mechanism and control measures ofoverlying strata with gob-side entry retaining formed byroof cutting and pressure releasing[J]. Journal of ChinaUniversity of Mining & Technology, 2019, 48(3) :474-483.
[ 5 ]康志鹏,段昌瑞,余国锋,等. 厚煤层软底沿空留巷围岩变形特征分析及顶帮强化支护技术[J]. 工矿自动化,2022,48(11):101-109.
KANG Zhipeng, DUAN Changrui, YU Guofeng, et al.Analysis on deformation characteristics of surroundingrock of gob-side entry retaining with soft bottom in thickcoal seam and strengthening support technology of roofand side[J]. Industry and Mine Automation, 2022,48(11):101-109.
[ 6 ]于光远,王炯,孙晗,等. 砌块充填沿空留巷底鼓大变形机理及综合控制技术[J]. 采矿与安全工程学报,2022,39(2):335-346.
YU Guangyuan, WANG Jiong, SUN Han, et al.Mechanism and comprehensive control techniques forlarge deformation of floor heave in block filling gob-sideentry retaining[J]. Journal of Mining & SafetyEngineering,2022,39(2):335-346.
[ 7 ]何满潮, 马资敏, 郭志飚, 等. 深部中厚煤层切顶留巷关键技术参数研究[J]. 中国矿业大学学报, 2018,47(3):468-477.
HE Mancha6/3a9ZVi3NYZJ7fa+k08fw==o, MA Zimin, GUO Zhibiao, et al. Keyparameters of the gob-side entry retaining formed byroof cutting and pressure release in deep mediumthicknesscoal seams[J]. Journal of China University ofMining & Technology,2018,47(3):468-477.
[ 8 ]张幼振. 煤矿巷道底鼓机械化治理技术应用研究[J].煤炭科学技术,2018,46(1):93-98.
ZHANG Youzhen. Applied study on mechanized controltechnology of floor heave in mine roadway[J]. CoalScience and Technology,2018,46(1):93-98.
[ 9 ]申斌学, 周宏范, 朱磊, 等. 深井复合顶板切顶卸压柔模墙支护沿空留巷技术[J]. 工矿自动化, 2021,47(11):101-106.
SHEN Binxue, ZHOU Hongfan, ZHU Lei, et al. Thetechnology of gob-side entry retaining supported byflexible formwork wall with roof cutting and pressurerelief for composite roof in deep mine[J]. Industry andMine Automation,2021,47(11):101-106.
[10]张宇旭,王科. 煤巷底鼓破坏特征及支护技术研究[J].工矿自动化,2019,45(6):73-79,85.
ZHANG Yuxu, WANG Ke. Research on failurecharacteristics of floor heave and support technology incoal roadway[J]. Industry and Mine Automation,2019,45(6):73-79,85.
[11]黄万朋,高延法,文志杰,等. 钢管混凝土支柱巷旁支护沿空留巷技术研究[J]. 中国矿业大学学报,2015,44(4):604-611.
HUANG Wanpeng, GAO Yanfa, WEN Zhijie, et al.Technology of gob-side entry retaining using concretefilledsteel tubular column as roadside supporting[J].Journal of China University of Mining & Technology,2015,44(4):604-611.
[12]王亚军,何满潮,张科学,等. 无煤柱自成巷开采巷道矿压显现特征及控制对策[J]. 采矿与安全工程学报,2018,35(4):677-685.
WANG Yajun, HE Manchao, ZHANG Kexue, et al.Strata behavior characteristics and controlcountermeasures for the gateroad surroundings ininnovative non-pillar mining method with gateroadformed automatically[J]. Journal of Mining & SafetyEngineering,2018,35(4):677-685.
[13]朱珍,张科学,袁红平. 切顶卸压沿空留巷碎石巷帮控制技术及应用[J]. 煤炭科学技术,2018,46(3):25-32.
ZHU Zhen,ZHANG Kexue,YUAN Hongping. Controltechnology and its application of roadway side wallformed by gangue in gob-side entry retaining formed byroof cutting and pressure releasing[J]. Coal Science andTechnology,2018,46(3):25-32.
[14]冯超,刘少伟,张小东,等. 切顶条件下综放沿空留巷充填体承载特征及合理宽度研究[J]. 采矿与安全工程学报,2023,40(2):232-242.
FENG Chao, LIU Shaowei, ZHANG Xiaodong, et al.Study on the bearing characteristics and reasonablewidth of the backfilling body for gob-side entry in fullymechanizedtop-coal caving face under the roof cuttingcondition[J]. Journal of Mining & Safety Engineering,2023,40(2):232-242.
[15]朱永建,任恒,王平,等. 倾斜厚层坚硬顶板条件下沿空留巷稳定性控制[J]. 中南大学学报(自然科学版),2023,54(3):956-966.
ZHU Yongjian,REN Heng,WANG Ping,et al. Stabilitycontrol of gob-side entry retention under the condition ofinclined thick layer and hard roof[J]. Journal of CentralSouth University (Science and Technology) , 2023,54(3):956-966.
[16]郑西贵,安铁梁,郭玉,等. 原位煤柱沿空留巷围岩控制机理及工程应用[J]. 采矿与安全工程学报,2018,35(6):1091-1098.
ZHENG Xigui, AN Tieliang, GUO Yu, et al.Surrounding rock control mechanism and engineeringapplication of in situ coal pillar in gob-side entryretaining[J]. Journal of Mining & Safety Engineering,2018,35(6):1091-1098.
[17]陈磊. 基于弹性基础梁模型的来压步距确定[J]. 煤矿安全,2017,48(5):49-52.
CHEN Lei. Determination of weighting interval basedon elastic foundation beam model[J]. Safety in CoalMines,2017,48(5):49-52.
[18]华心祝,李迎富. 沿空留巷底板变形力学分析及底臌防控[J]. 煤炭学报,2016,41(7):1624-1631.
HUA Xinzhu, LI Yingfu. Mechanics analysis on floordeformation of gob-side entry retaining and preventionand control of floor heave[J]. Journal of China CoalSociety,2016,41(7):1624-1631.
[19]许斌,蒋金泉,代进,等. 采场上覆关键层破断角的力学推导和实验模拟[J]. 煤炭学报, 2018, 43(3) :599-606.
XU Bin, JIANG Jinquan, DAI Jin, et al. Mechanicalderivation and experimental simulation of breaking angleof key strata in overlying strata[J]. Journal of ChinaCoal Society,2018,43(3):599-606.
[20]马念杰,侯朝炯. 采准巷道矿压理论及应用[M]. 北京:煤炭工业出版社,1995.
MA Nianjie,HOU Chaojiong. The underground pressureof sectional roadways and its control[M]. Beijing:ChinaCoal Industry Publishing House,1995.
[21]徐芝纶. 弹性力学简明教程[M]. 3 版. 北京:高等教育出版社,2002.
XU Zhilun. A concise course in elasticity[M]. 3rd ed.Beijing:Higher Education Press,2002.
[22]周道祥,张伟林. 工程弹性力学基础[M]. 合肥:合肥工业大学出版社,2006.
ZHOU Daoxiang, ZHANG Weilin. Fundamentals ofengineering elasticity[M]. Hefei: Hefei Univerisity ofTechnology Publishing House,2006.
[23]何满潮,高玉兵,杨军,等. 无煤柱自成巷聚能切缝技术及其对围岩应力演化的影响研究[J]. 岩石力学与工程学报,2017,36(6):1314-1325.
HE Manchao, GAO Yubing, YANG Jun, et al. Anenergy-gathered roof cutting technique in no-pillarmining and its impact on stress variation in surroundingrocks[J]. Chinese Journal of Rock Mechanics andEngineering,2017,36(6):1314-1325.
[24]陆银龙,韩磊,吴开智,等. 特厚煤层沿空掘巷力源结构特征与围岩协同控制策略[J]. 中国矿业大学学报,2024,53(2):238-249.
LU Yinlong, HAN Lei, WU Kaizhi, et al.Characteristics of stress sources and comprehensivecontrol strategies for surrounding rocks of gob-sidedriving entry in extra thick coal seam[J]. Journal ofChina University of Mining & Technology, 2024,53(2):238-249.
基金项目:国家自然科学基金资助项目(52074008)。