Evolution of broken structure of hard thick roof and analysis of impact hazard area
-
摘要:
针对硬厚顶板破断及结构演化对工作面冲击显现的影响,以东滩煤矿63上05工作面硬厚顶板条件为背景,运用理论分析和数值模拟分析了硬厚顶板破断失稳的临界条件及其结构演化过程,并划分工作面冲击危险区域。结果表明:顶板厚度和强度与块体破断时端部垂直载荷成正比,直接顶与基本顶的结构形态由“V”形演化成“U”形;顶板初次和周期跨落步距相较于常规情况明显增大,且基本顶结构演化较直接顶具有迟滞性;顶板垂直应力以应力拱形式演化,开采过程中在采空区中形成多处应力集中。基于研究结果,划分63上05工作面冲击危险区域并利用微震监测技术验证准确性。
Abstract:Aiming at the influence of hard and thick roof fracture and structural evolution on the impact of working face, taking the hard and thick roof condition of 63 upper 05 working face in Dongtan Coal Mine as the engineering background, the critical condition and structural evolution process of hard and thick roof fracture and instability are studied by theoretical analysis and numerical simulation, and the dangerous area of working face impact is divided. It is concluded that the thickness and strength of the roof are proportional to the vertical load at the end of the block breaking. The structure of immediate roof and main roof gradually evolved from “V” type to “U” type. The first and periodic steps of the double hard thick roof increase obviously compared with the conventional condition, and the evolution of the basic roof structure is more retarded than that of the direct roof. The vertical stress of the roof evolves in the form of stress arch, and many stress concentration phenomena are formed in the mined-out area during mining. Based on the above research results, the impact danger area of 63 upper 05 working face was divided and the accuracy was verified by microseismic monitoring technology.
-
水资源的保护与利用是我国煤矿绿色开采亟待解决的难题之一,特别是在我国西部地区,煤炭资源约占全国总储量的80%,但是气候干旱,蒸发量是降雨量的6倍左右,再加上西部矿区大规模、高强度开采对其生态水环境造成严重的影响[1-2]。因此针对煤炭资源开采与水资源保护相协调的难题,顾大钊等[3-5]提出了煤矿地下水库的概念及矿井水地下储存利用的理论与技术,有利于西部矿区煤炭开采与水资源保护。
针对煤矿地下水库的稳定性,诸多学者进行了大量的研究,取得了丰富的成果。姚强岭等[6-7]研发了无损浸水实验装置,分析了不同含水状态的煤柱坝体损伤特征,并结合理论分析研究了煤矿地下水库煤柱坝体的合理宽度计算方法及稳定性;张村等[8]、HAN等[9]开发了采动水浸作用下煤矿地下水库流固耦合数值模型,并通过定义煤柱损伤度研究了工程尺度下煤柱坝体的失稳特征;师维刚等[10]通过对防水隔离煤柱进行合理分区提出了计算煤柱宽度理论方法;白东尧等[11]通过数值模拟及理论分析得出了地下水库人工坝体的极限水头,并对地下水库人工坝体的稳定性展开研究;刘少伟等[12]通过理论计算及数值模拟研究了影响缓倾斜煤层沿含水采空区掘巷时影响护巷煤柱稳定性的因素敏感性;王文才等[16]通过实验室测试及数值模拟研究分析了作业条件下矿井地下水库坝体衔接处受力,变形情况及其失效机制;李刚等[14]通过理论分析及数值模拟相结合研究了采动过程中底板应力场分布特征及岩层破坏损伤演化特征;刘祥龙等[15]针对圆柱形试件进行不同侧向渗透压力下的煤样压缩试验,研究煤样的破坏特征与渗透特性,将煤样在侧向渗透-轴向承载下的变形分为5个阶段;刘晓飞等[16]通过理论分析及实验室测试等分析了工作面末采阶段煤柱坝体所受应力环境以及研究了周期应力及水浸作用对媒体的稳定性影响。
综上,对于地下水库煤柱坝体的稳定性研究多集中在实验室尺度下对浸水煤岩样强度及损伤特性研究以及通过理论分析对影响煤柱坝体稳定的临界水头研究,而在工程尺度下,部分学者通过数值模拟提出了力学计算与渗流计算交替进行的间接流固耦合方法,研究采动水浸作用下煤柱坝体的稳定性[8-9]。然而,工程尺度下煤柱坝体在完全流固耦合作用下对其承载能力的动态劣化特征及损伤破坏规律的影响仍值得深入研究。同时基于煤柱的多孔介质特性,FLAC3D数值模拟软件在模拟地下水库煤柱坝体在复杂载荷作用下处于弹性未破坏阶段时的稳定性存在一定的局限性。为此,通过数值模拟及二次开发,构建了基于煤矿地下水库煤柱坝体在流固耦合中多场耦合作用下力学参数动态弱化的Fish语言并嵌入FLAC3D中,实现煤柱坝体在流固耦合作用下力学参数的动态调整以减弱数值模拟进行流固耦合分析时的不足;从而研究了地下水库煤柱坝体在采动及水浸叠加作用下的稳定性,得到了地下水库煤柱坝体在复杂条件下应力分布响应特征及塑性区扩展变化规律,为地下水库煤柱坝体稳定性研究提供借鉴。
1. 工程背景和流固耦合数值模型的构建
1.1 工程背景
以陕西北部某矿主要开采煤层2-2中的2个相邻工作面开采为例,煤层倾角为0°~5°,煤层厚度为8.89~11.73 m,平均厚度为10.41 m。双面开采及采空区蓄水过程如图1所示。
P122106工作面为首采工作面,走向长度为6 200 m,倾斜长度为350 m,开采高度为6 m;P122106长壁工作面南部为P122108长壁工作面,在2个工作面之间留设安全煤柱,煤柱拟留设宽度为40 m,并在工作面开挖稳定后,开始在工作面采空区进行水库建设并进行蓄水。因此,对地下水库煤柱坝体在采动作用及上覆载荷-水浸流固耦合依次作用下的稳定性研究对地下水库长期安全运行非常重要。
1.2 数值模型的构建
按照前述工程概况,建立FLAC3D流固耦合数值仿真计算模型。模型以x方向为倾向,y方向为工作面走向,z方向为垂直于xy平面;模型尺寸为x=1 000 m,y=400 m,z=145 m;在模型的四周及底部边界采用位移约束,上部为自由边界,共模拟2个相邻工作面的开采,为了使煤柱的监测区域更为准确,将其网格细化为0.5 m。数值模拟中材料参数取值见表1,长壁开采工作面数值模型如图2所示。
表 1 数值模拟中材料参数取值Table 1. Values of material parameters in numerical simulation岩性 密度/
(kg·m−3)体积模量/
GPa剪切模量/
GPa内摩擦角/
(°)黏聚力/
MPa抗拉强度/
MPa渗透系数/
(nm·s−1)孔隙率 中粒砂岩 2 500 11.61 9.09 28 1.00 2.30 0.060 0.45 砂质泥岩 2 400 2.57 1.45 26 0.55 0.80 0.006 0.35 粉砂岩 2 500 9.63 8.09 30 2.50 0.84 50.000 0.40 泥岩 2 141 2.21 1.24 24 0.52 0.83 65.000 0.40 细粒砂岩 2 400 2.13 0.93 26 0.45 0.35 0.006 0.30 煤层 1 400 2.13 0.93 25 0.40 0.25 80.000 0.50 模型的中上覆各岩层及底板均采用Mohr-coulomb模型,煤层采用基于Mohr-coulomb的应变软化模型;同时,由于采空区逐渐压实过程影响工作面周围应力分布,开发了基于salamon应力应变曲线的模拟采空区垮落破碎矸石的承载特性的Fish语言嵌入FLAC3D中。
1.3 基于煤柱坝体的水浸软化模型
对地下水库煤柱坝体进行流固耦合分析时,应力和水压是影响煤柱稳定的2个关键物理场,两者的耦合作用是导致煤柱坝体失稳破坏的重要原因;而FLAC3D流固耦合中通过监测模型内部存在孔隙水压力和处于压缩状态时对其力学参数进行软化。实际上,煤柱坝体在弹性状态下,内部处于压密状态下的孔隙闭合会形成弹性压密隔水区,其渗透率显著降低,为煤柱的主要承载区域,此时,这部分煤体应表现为应变硬化特性。
为了更加准确的表现地下水库煤柱坝体在流固耦合作用下的破坏失效形式,对基于煤柱坝体的应变软化模型进行了二次开发,得出煤柱坝体内部应力场、塑性区、能量场多场耦合作用下力学参数动态更新的Fish语言,并嵌入FLAC3D中。基于煤柱的水浸软化模型计算流程图如图3所示。
Fish语言能够捕捉煤柱坝体内部各单元的弹塑性状态:①当目标单元的处于塑性状态时,根据现有研究[17]确定煤柱在饱和状态下的弹性模量弱化系数kE=0.429 2, 黏聚力弱化系数kc=0.572 5, 内摩擦角弱化系数kf=0.960 1,从而模拟煤柱内部塑性破坏区承载能力的劣化及失效特征;②当目标单元处于弹性状态且体积应变小于0,确定煤柱各力学参数的弱化系数为1,从而减弱FLAC3D中流固耦合效应对其力学参数弱化的影响,最终实现对煤柱内部弹性压密隔水区的模拟。
2. 煤柱坝体稳定性
2.1 多工作面采动影响
由于矿方拟留设煤柱宽度为40 m左右,为了研究采动及水浸叠加作用下煤柱坝体的承载特征及渐进失效规律,以煤柱坝体宽40 m作为稳定性的研究对象,同时,为了对比研究煤柱坝体在水浸及采动叠加作用下的承载特征及失稳规律,还研究了煤柱留设宽度小于40 m(20 m和30 m)时的稳定性。煤柱坝体塑性区和垂直应力分布如图4和图5所示,不同宽度煤柱坝体垂直应力分布特征与弹性核区占比如图6所示。
由图6所可以看出,3种宽度的煤柱坝体垂直应力分布曲线均呈现双峰形分布特征,且煤柱右侧的垂直应力峰值均大于左侧,煤柱宽度为20、30、40 m时其垂直应力峰值分别为44.00、35.15、29.13 MPa,呈现出随着煤柱宽度的增加应力峰值随之减小。煤柱坝体2个垂直应力峰值之间为弹性核区,煤柱宽度为20、30、40 m时,弹性核区占比分别为39.73%、73.68%、83.54%。为了保证煤柱坝体的稳定性,煤柱内部需至少维持一定的弹性核区,现场实践表明隔水煤柱内部弹性核区占比应不小于40%,以保证在复杂载荷作用下保持稳定[18],并有效防止动力灾害的发生。因此,在开采过后,煤柱坝体留设宽度为20 m时,承载能力劣化并会发生失稳。表明随着煤柱宽度的增加其承载能力随之增强,煤柱的稳定性越好,且煤柱坝体留设宽度为30 m和40 m时能够保持稳定。
2.2 水浸作用下煤柱坝体稳定性
在工作面开挖稳定后,开始向采空区进行蓄水。此时,地下水库煤柱坝体会在水浸及上覆岩层载荷的流固耦合作用下强度进一步降低,且不同于传统的承载式煤柱,地下水库煤柱坝体在长期水浸作用下更容易失效;因此,在工作面开挖稳定后在采空区通过“watertable”命令建立水平面及各渗流参数,并根据上述煤柱流固耦合数值模拟方法,对不同宽度煤柱坝体在水浸作用下破坏失稳特征展开研究。渗流作用下不同宽度煤柱坝体塑性区如图7所示,渗流作用下不同宽度煤柱坝体垂直应力分布云图如图8所示,渗流作用下不同宽度煤柱坝体垂直应力及塑性区分布特征如图9所示。
在水浸作用下,地下水库煤柱坝体宽度为30 m时,垂直应力曲线呈类似“拱形”分布,“拱顶”部分为煤柱坝体的主要承载区域,也是主要隔水区域。应力峰值为26.75 MPa,位于煤柱左侧,相较于采空区蓄水之前,降幅为23.9%;随着煤柱宽度的增加,垂直应力分布曲线中部主要承载区域的宽度也相应增加,应力峰值为24.11 MPa,位于煤柱左侧,相较于煤柱宽度为30 m时有所下降,采空区蓄水之前降幅为17.2%,结果表明,煤柱坝体宽度越小,流固耦合作用对其承载能力的影响越大,增加煤柱宽度有利于保持其稳定性。
在水浸作用下,不同宽度煤柱坝体的塑性区向煤柱内部进一步扩展,塑性区边界煤柱坝体内部为高垂直应力区域构成煤柱坝体的有效阻水区域。当煤柱宽度为30 m时,弹性核区占比为38.33%,相较于采空区蓄水之前降幅为48%;随着煤柱坝体宽度的增加其弹性核区占比显著增加;当煤柱宽度为40 m时,弹性核区占比为59.53%,相较于采空区蓄水之前降幅为28.7%。结果表明,煤柱宽度越小,流固耦合作用对其强度的影响越大,增加煤柱宽度有利于保持其稳定性。然而,当煤柱留设宽度过大时,一方面会浪费大量的煤炭资源,另一方面会产生一些煤柱型动力灾害。因此,为了保证地下水库蓄水后的系统安全,当煤柱宽度为40 m时,内部弹性核区占比将大于40%,能够在采动及水浸依次作用下保持稳定。
采空区积水水浸作用下煤柱内部力学参数大幅弱化,导致塑性区进一步扩展及承载能力显著降低。运行过程中不同宽度煤柱塑性区扩展演化规律如图10所示,运行过程中煤柱宽度为40 m时黏聚力弱化过程如图11所示。
由图11可知:煤柱坝体右侧的损伤范围逐渐由边缘向内部扩展,并呈现出煤柱下方区域先于上方区域,边界煤体先于内部煤体的渐进失效特征,同时煤柱边界区域受开采扰动作用的影响较大,初始及最终损伤破坏程度均大于煤柱内部区域,且在采动及水浸作用下煤柱坝体中部区域始终未受到影响,进一步表明煤柱坝体宽度为40 m时能够保持稳定。
综上所述,采空区地下水库蓄水后,积水将沿着煤柱坝体在采动作用下形成的塑性区高导水通道侵入煤体,并向煤柱内部不断渗透,使煤体内部原有结构的黏结强度减弱以及裂隙面间摩擦系数减弱,从而增加了煤体破坏的可能性,导致煤体承载能力降低,从而进一步导致煤柱塑性区扩展并引起应力重分布向煤柱内部转移,并促使积水向煤柱内部继续渗透,直至煤柱内部高应力区形成稳定阻水区域或煤柱内部积水区域贯通造成煤柱坝体破坏。为此,构建的数值模型能够很好地体现地下水库煤柱坝体在采动及水浸依次作用下的渐进破坏过程。在采空区积水水浸作用下,煤柱坝体塑性区占比增量呈现“稳定-增大-减小-稳定”的变化规律,能够很好地模拟水浸作用下煤柱坝体的“渗流-弱化-损伤-渗流”的渐进破坏特征。
3. 结 语
1)基于FLAC3D自带的流固耦合分析在模拟煤柱弹性孔隙压密阶段存在的不足,开发了基于煤柱多场耦合作用下参数动态调整的Fish语言并嵌入FLAC3D中,实现煤柱坝体在流固耦合作用下力学参数的动态调整以实现煤柱在弹性压密阶段的模拟。
2)研究了多工作面开采对地下水库煤柱坝体稳定性的影响规律,结果表明:在相邻工作面开采后,不同宽度的隔水煤柱应力分布曲线均呈现双峰形分布,并随着煤柱宽度的增加其峰值应力逐渐降低,同时煤柱内部弹性核区占比也随煤柱留设宽度的增加而增大,当煤柱宽度为20 m时其弹性核区占比为39.73%,将发生破坏失稳,而煤柱宽度为30 m和40 m时煤柱能够保持稳定。
3)运用流固耦合方法对煤柱蓄水水浸作用下的稳定性进行了分析。结果表明:采空区蓄水后,水浸作用下煤柱坝体塑性区占比增量呈现“稳定-增大-减小-稳定”的变化规律,能够很好地模拟水浸作用下煤柱坝体的“渗流-弱化-损伤-渗流”的渐进破坏特征,水浸作用对煤柱坝体稳定性的影响更为显著,是导致其破坏失效的主要因素;同时煤柱内部应力峰值降低,煤柱宽度越小应力峰值降低幅度越大,且煤柱宽度为30 m时其垂直应力分布曲线从双峰形演化为拱形分布,同时其弹性核区占比降为38.33%,降幅为48%,将发生失稳;而煤柱宽度为40 m时,能够在采动及水浸依次作用下保持稳定。
-
表 1 煤岩层力学参数
Table 1 Mechanical parameters of coal and rock
岩 石
类 型密度/
(kg·m−3)体积模
量/GPa剪切模
量/GPa黏聚力/
MPa抗拉强
度/MPa内摩擦
角/(º)泥岩 2 500 5.13 8.75 2.38 3.24 26 中砂岩 2 700 8.10 20.56 8.15 5.78 30 泥岩 2 550 5.13 8.75 2.38 3.24 26 粉砂岩 2 700 8.12 12.87 4.30 5.25 28 煤 1 450 2.45 5.98 2.76 2.69 35 粉砂岩 2 800 8.12 12.87 4.30 5.25 28 细砂岩 2 900 9.23 21.34 7.56 6.89 32 -
[1] 侯玮,霍海鹰,郭晓朋. 深部综放工作面动压微震监测及控制技术研究[J]. 煤炭科学技术,2014,42(3):36−38. HOU Wei, HUO Haiying, GUO Xiaopeng. Research on microseismic monitoring and control technology of dynamic pressure in deep fully-mechanized top coal caving face[J]. Coal Science and Technology, 2014, 42(3): 36−38.
[2] 陈殿赋. 采空区下坚硬顶板动压显现特征及控制技术[J]. 煤炭科学技术,2014,42(10):125−128. CHEN Dianfu. Reserch on strata behaviors features and control technology of hard roof underground goaf[J]. Coal Science and Technology, 2014, 42(10): 125−128.
[3] 谭云亮,张明,徐强,等. 坚硬顶板型冲击地压发生机理及监测预警研究[J]. 煤炭科学技术,2019,47(1):166−172. TAN Yunliang, ZHANG Ming, XU Qiang, et al. Study on occurrence mechanism and monitoring and early warning of rock burst caused by hard roof[J]. Coal Science and Technology, 2019, 47(1): 166−172.
[4] 刘海洋,孟凡林. 厚煤层工作面顶板破断与能量场动态演化规律[J]. 煤矿安全,2023,54(3):153−160. LIU Haiyang, MENG Fanlin. Roof breaking and dynamic evolution of energy field in thick coal seam working face[J]. Safety in Coal Mines, 2023, 54(3): 153−160.
[5] 姜福兴,曲效成,倪兴华,等. 鲍店煤矿硬岩断裂型矿震的预测[J]. 煤炭学报,2013,38(S2):319−324. JIANG Fuxing, QU Xiaocheng, NI Xinghua, et al. Case study on the mine earthquake caused by hard rock fracture in Baodian Coal Mine[J]. Journal of China Coal Society, 2013, 38(S2): 319−324.
[6] 吕鹏飞,陈学华,周年韬. 高位硬厚岩层影响下矿震发生规律及预测[J]. 安全与环境学报,2018,18(1):95−100. LYU Pengfei, CHEN Xuehua, ZHOU Niantao. Mining earthquake incidence features and forecast under the impact of the top-level hard and thick rock and coal seams[J]. Journal of Safety and Environment, 2018, 18(1): 95−100.
[7] 齐庆新,潘一山,李海涛. 煤矿深部开采煤岩动力灾害防控理论基础与关键技术[J]. 煤炭学报,2020,45(5):1567−1584. QI Qingxin, PAN Yishan, LI Haitao. Theoretical basis and key technology of prevention and control of coal-rock dynamic disasters in deep coal mining[J]. Journal of China Coal Society, 2020, 45(5): 1567−1584.
[8] 王兆会,程占博. “两硬”条件下孤岛型短煤柱工作面顶板破断形态及灾害防治分析[J]. 岩石力学与工程学报,2016,35(S2):4018−4028. WANG Zhaohui, CHENG Zhanbo. Hard roof fracturing form and dynamic disaster control in short island mining face[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S2): 4018−4028.
[9] 刘瑞瑞,刘洋,方刚,等. 袁大滩煤矿覆岩破坏规律及顶板水害防治对策[J]. 煤矿安全,2022,53(7):82−91. LIU Ruirui, LIU Yang, FANG Gang, et al. Law of overburden failure and roof water damage in YuandatanCoal Mine[J]. Safety in Coal Mines, 2022, 53(7): 82−91.
[10] 张玉军,张志巍. 煤层采动覆岩破坏规律与控制技术研究进展[J]. 煤炭科学技术,2020,48(11):85−97. ZHANG Yujun, ZHANG Zhiwei. Research progress of mining overlying stratas failure law and control technology[J]. Coal Science and Technology, 2020, 48(11): 85−97.
[11] 窦林名,贺虎. 煤矿覆岩空间结构OX-F-T演化规律研究[J]. 岩石力学与工程学报,2012,31(3):453−460. doi: 10.3969/j.issn.1000-6915.2012.03.003 DOU Linming, HE Hu. Study of OX-F-T spatial structure evolution of overlying strata in coal mines[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(3): 453−460. doi: 10.3969/j.issn.1000-6915.2012.03.003
[12] 姜福兴,张兴民,杨淑华,等. 长壁采场覆岩空间结构探讨[J]. 岩石力学与工程学报,2006,25(5):979−984. doi: 10.3321/j.issn:1000-6915.2006.05.018 JIANG Fuxing, ZHANG Xingmin, YANG Shuhua, et al. Discussion on overlying strata spatial structures of longwall in coal mine[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(5): 979−984. doi: 10.3321/j.issn:1000-6915.2006.05.018
[13] 张宏伟,付兴,于斌,等. 特厚煤层坚硬覆岩柱壳结构特征模型及应用[J]. 中国矿业大学学报,2017,46(6):1226−1230. ZHANG Hongwei, FU Xing, YU Bin, et al. Study of the cylindrical shell structural characteristic of the hard overlying strata of extra thick coal seam and its application[J]. Journal of China University of Mining & Technology, 2017, 46(6): 1226−1230.
[14] 朱卫兵,于斌. 大空间采场远场关键层破断形式及其对矿压显现的影响[J]. 煤炭科学术,2018,46(1):99−104. ZHU Weibing, YU Bin. Breakage form and its effect on strata behavior of far field key stratum in large space stope[J]. Coal Science and Technology, 2018, 46(1): 99−104.
[15] 朱卫兵,于斌,鞠金峰,等. 采场顶板关键层“横U-Y”型周期破断特征的试验研究[J]. 煤炭科学技术,2020,48(2):36−43. ZHU Weibing, YU Bin, JU Jinfeng, et al. Experimental study on “horizontal U-Y” periodical breakage characteristics of key strata in stope roof[J]. Coal Science and Technology, 2020, 48(2): 36−43.
[16] 于斌,朱卫兵,李竹,等. 特厚煤层开采远场覆岩结构失稳机理[J]. 煤炭学报,2018,43(9):2398−2407. YU Bin, ZHU Weibing, LI Zhu, et al. Mechanism of the instability of strata structure in far field for super-thick coal seam mining[J]. Journal of China Coal Society, 2018, 43(9): 2398−2407.
[17] 赵科,张开智,王树立. 巨厚覆岩破断运动与矿震活动规律研究[J]. 煤炭科学技术,2016,44(2):118−122. ZHAO Ke, ZHANG Kaizhi, WANG Shuli. Study on the law of fracture movement and mine earthquake activity of super thick overburden rock[J]. Coal Science and Technology, 2016, 44(2): 118−122.
[18] 梁沙平,陆银龙,郭鹏,等. 特厚煤层坚硬顶板初次破断特征的力学分析[J]. 煤矿安全,2020,51(8):245−250. LIANG Shaping, LU Yinlong, GUO Peng, et al. Mechanical analysis of the first fracture characteristics of hard roof of extra-thick coal seam[J]. Safety in Coal Mines, 2020, 51(8): 245−250.