不同表层厚度下含断层煤体开采覆岩下沉规律.pdf

第 11 期 山 西 焦 煤 科 技 No.112023 年 11 月Shanxi Coking Coal Science&TechnologyNov.2023试验研究收稿日期:2023-08-20作者简介:曹呆军(1975),男,山西晋城人,2014 年毕业于太原工业大学,工程师,主要从事煤矿生产管理工作(E-mail)15966757654 不同表层厚度下含断层煤体开采覆岩下沉规律曹呆军(山西兰花同宝煤业有限公司,山西晋城048400)摘要煤层开采与断层滑移将会导致覆岩下沉变形,为探究断层影响下不同表层厚度采场的覆岩变形特征,采用理论分析与数值模拟方法分析了全基岩、薄松散层、厚松散层条件下的地表沉陷规律,建立了厚松散层条件下的基岩稳定性判据,揭示不同表层厚度下的断层带覆岩变形运移规律。结果表明:地表下沉曲线的发育范围与光滑程度会随表层厚度的增加而提升。全基岩条件下地表产生台阶式下沉结构,薄松散层条件下地表产生间断性地表裂缝,厚松散层条件下地表产生明显的非对称式的连续性下沉变形。关键词断层构造;表层厚度;地表变形;覆岩下沉中图分类号:TD325.2文献标识码:B文章编号:1672-0652-(2023)11-0030-04煤层开采后采空区发生应力重新分布,导致上覆岩层 运 移 与 变 形,最 终 影 响 到 地 表 形 成 下 沉 盆地1-4.表层厚度以及地质构造会对覆岩的运移规律与应力分布状态产生影响。针对上述开采沉陷问题,张玉卓等5探究了断层倾角与岩层移动角等对覆岩运移规律的影响,表明以上因素对地表沉陷范围具有控制作用。尹立明、于秋鸽等6-7分析了断层滑移诱发覆岩破坏的内在机理与演化过程。孔宪森等8运用数值模拟方法与理论解析方法研究了断层带活化对表层建筑物沉降的影响。以上研究为断层条件下的地表沉陷规律研究起到了促进作用。为深入探究不同表层厚度对覆岩下沉规律与地表变形的影响,拟采用理论分析与数值模拟方法进行研究。1覆岩下沉规律与地表变形1.1全基岩覆岩下沉规律与地表变形在全基岩地质条件下,地表上方不存在覆盖的松散层,因此,地表受断层的影响较大。断层直接破坏了覆岩的完整性,同时阻隔了覆岩应力的传播。当下方煤层进行开采时,上覆岩层发生连续的垮落变形,而这种连续变形在断层位置被切断。这种断层切断作用控制了地表变形的范围,并在断层带附近形成明显的台阶下沉结构。覆岩下沉规律与地表变形示意图见图 1.图 1全基岩覆岩下沉规律与地表变形示意图当断层结构受到开采的影响产生拉剪滑移破坏变形,破断的采空区上方基岩会向采空区的内部移动,下沉盆地的断层边界区域会出现抬高现象,同时伴随着大量张拉裂缝的衍生。以高度为 h 的基岩衍生裂缝为例,当基岩厚度为 H,基岩的破坏距离 lL时(L 为基岩的极限破坏距离),则基岩的悬臂结构不发生连续性破坏,断层上盘的变形量较小。相应的,当 lL 时,断层的上盘将会向断层方向垮落,并导致大量的地表移动变形。1.2薄松散层覆岩下沉规律与地表变形在薄松散层的地质条件下,下盘煤层开采后导致上覆岩层下沉9-10.断层起到了限制上盘岩层移动和变形的作用。然而,由于差异性沉降,上下盘之间产生了位移不连续。图 2 直观地展示了所涉及的变形状态。松散层在覆盖基岩断层露头位置承受剪切力,这些剪切力是由松散层自身的重力和上盘岩层的支撑力所引起的。当支撑力和重力共同作用产生的剪切力超过松散层的抗剪强度时,松散层发生剪切破坏。这种剪切破坏表现为间断性的地表裂缝。图 2薄松散层覆岩下沉规律与地表变形示意图薄松散层发生剪切破坏后,松散层厚度与覆岩下沉量的比值不同,地表变形形式存在相应差异。随着松散层厚度的占比逐渐增加,地表的台阶式下沉结构将逐渐转变为间断性地表裂缝。1.3厚松散层覆岩下沉规律与地表变形厚松散层受到扰动影响时,松散层下部在采空区边缘将会发生滑移,而中部容易发生剪切破坏11-12.然而,由于厚松散层具有吸收非连续变形的能力,当断层存在于基岩中且上覆厚松散层时,地表表现出连续性的下沉,而非台阶式下沉或间断性地表裂缝。如图 3 所示在无断层情况下,地表呈现对称的下凹形状。然而,当断层存在时,随着煤层开采逐渐靠近断层,断层活化,以断层上盘露头位置为临界点,诱发上覆厚松散层发生移动变形,使得松散层滑移的位置靠近采空区。相比无断层一侧,上盘地表盆地的变形范围较小。随着上覆松散层厚度的增加,松散层以荷载的形式作用于基岩岩层,导致断层上盘基岩岩层以下盘开采扰动区域形成底部临界点,呈现三角形破坏区域。在这种情况下,松散层滑移破坏的位置远离采空区,导致地表移动变形范围和变形值增加。将图 3 所示的断层破坏区域进行力学建模,断层带附近岩石抵抗滑移变形的能力来源于基岩岩石的内聚力与自身重力,而促使其发生滑移变形的能力来图 3厚松散层覆岩下沉规律与地表变形及力学模型示意图源于松散层的重力以及破坏区域的岩石重力。根据二者的比值可推导破坏区域稳定性判据:Fs=Fk/F(1)式中:Fs为稳定性判据;Fk为岩体抵抗滑移变形的阻力,kN;F 为促使岩体滑移变形的驱动力,kN.破坏区域长度与基岩厚度的几何关系:L=h2/tan+h2/tan(2)式中:L 为破坏区域长度,m;h2为岩体破坏区域厚度,m;为岩体滑动面与水平面夹角,();为断层倾角,().F 与 Fk的力学关系:F=h1L+h2L/2tansin(3)Fk=ch2/sin+h1L+h2L/2tancostan(4)式中:h1为松散层厚度,m;为容重,N/m3;为岩体内摩擦角,();c 为岩体内聚力,MPa.将上式进行联立,得到断层带滑移破坏判据:Fs=2ctan(tan+tan)(2h1+h2)cos2+tantan(5)根据式(1)与式(5)进行分析,当 Fs1 时,岩体抗滑能力高于驱动力,岩体将保持稳定;当 0Fs1时,岩体则发生剪切滑移破坏,并且比值越小,发生滑移破坏的范围越大。2数值模拟模型的建立为再现不同表层厚度下含断层煤体开采过程中的覆岩下沉规律,对覆岩的下沉运动进行可视化分132023 年第 11 期曹呆军:不同表层厚度下含断层煤体开采覆岩下沉规律析,采用 FLAC3D数值软件进行研究。2.1数值模拟模型数值模拟的岩层共设置 5 层,分别为松散层、中砂岩、粉砂岩、煤层与泥岩,并在其中布置断层,断层的宽度为 2 m,断层落差为 10 m,倾角为 70.为模拟松散层与断层的节理弱化特点,采用 VORONOI 随机节理划分方式,数值模拟的岩性参数见表 1,表 2.模型的总高度为 325 m,长度为 1000 m,松散层与基岩的总厚度为 300 m,煤层开采厚度为 5 m,煤层底板厚度为 5 m.数值模拟模型见图 4,固定模型的左右边界与底部边界,对各组岩层施加重力,模型采用摩尔-库伦准则进行计算。表 1数值模拟岩性参数表岩层体积模量/GPa剪切模量/GPa内聚力/MPa密度/(kg m-3)内摩擦角/()抗拉强度/MPa松散层0.980.700.15110080.10中砂岩13.8011.265.332600352.55粉砂岩8.896.553.212300332.21煤层0.520.340.351300211.50泥岩4.222.861.561800241.90断层0.500.400.20800100.25表 2数值模拟松散层与断层节理参数表节理参数内摩擦角/()内聚力/MPa法向刚度/GPa切向刚度/GPa抗拉强度/MPa松散层20.041.50.40.01断层20.020.80.20.01图 4数值模拟设计图2.2数值模拟方案实验共设计 3 种方案,分别为全基岩条件、薄松散层条件与厚松散层条件。通过调节松散层与基岩的比例对不同的表层厚度进行模拟,全基岩条件的松散层 占 比 为 0%、薄 松 散 层 条 件 的 松 散 层 占 比 为17%、厚松散层条件的松散层占比为 67%,具体方案参数见表 3.3地表下沉规律比较分析不同方案下的地表下沉曲线见图 5,当覆岩无松散层时,地表的最大下沉量为 3.2 m;当覆岩存在薄松散层时,地表最大下沉量为 3.6 m;当覆岩存在厚松散层时,地表的最大下沉量为 4.2 m.根据图 5 分析可知,在全基岩条件下,模型在走向的 320 m 位置出现了台阶式的跌落点,原因是由于断层切断作用控制了地表变形的范围,这与前文所述的论点一致。随着松散层厚度逐渐增加,地表下沉曲线也逐渐变得光滑,同时由于松散层的岩性较弱,导致了地表下沉的总量逐渐增大。表 3数值模拟方案参数表岩层全基岩薄松散层厚松散层上覆松散层/m050200上覆基岩/m300250100煤层/m555煤层底板/m202020松散层比例/%01767图 5不同表层厚度条件下地表下沉规律图4不同表层厚度下开采覆岩下沉规律不同表层厚度下采场覆岩数值模拟变形云图见图 6.如图 6(a)所示,受到断层的阻隔作用导致断层下盘与上盘出现了明显的滑移错位,在地表产生了明显的下沉台阶,上盘与下盘的错位高度约为 1.06 m.数值模拟结果验证了全基岩覆岩下沉规律与地表变形规律。根据图 6(b)所示,断层在受到采动影响后发生23山 西 焦 煤 科 技2023 年第 11 期剪切滑移变形,其尖端产生的断层裂缝向地表延伸。在延伸线方向,断层上盘与 下盘 的 错 位 高 度 约 为0.11 m,并伴随产生地表裂缝,并且向断层下盘发生偏折。由于松散层吸收了部分断层的剪切滑移变形释放的能量,产生了间断性的地表裂缝。数值模拟结果验证了薄松散层覆岩下沉规律与地表变形。根据图 6(c)中所示的 0.1 m 下沉等值线所示,断层影响下的覆岩变形呈现出明显的非对称形式,断层一侧的地表下沉量大于无断层一侧的地表下沉量。随着松散层厚度的增加,断层带剪切滑移变形释放的能量大部分被松散层所吸收,并未在地表产生明显的台阶式下沉或者间断性地表裂缝。数值模拟结果验证了厚松散层覆岩下沉规律与地表变形。5结语1)分析了全基岩条件下、薄松散层条件下、厚松散层条件下的地表变形与覆岩下沉规律。根据断层影响作用的不同,地表会相应的表现出台阶式下沉变形、间断性地表裂缝与连续性大变形特征,并建立了厚松散层的基岩稳定性判据。2)全基岩条件下、薄松散层条件下、厚松散层条件下的地表最大下沉量分别为 3.2 m、3.6 m 与 4.2 m,同时地表下沉曲线的影响范围和光滑程度也逐渐增加。图 6不同表层厚度下开采覆岩下沉规律图 3)通过对开采条件下的覆岩变形云图进行分析,对 3 种条件下的地表变形规律进行了有效验证。同时表明,断层影响下厚松散层覆岩在开采时将会产生明显的非对称式变形特征。参考文献1杨树才.薄煤层综采地表移动变形规律研究J.煤,2022,31(7):64-66.2赵燕兵.古交矿区重复采动沉陷规律及沉陷参数分析J.山西焦煤科技,2017,41(10):48-52,56.3康新亮.官地矿重复采动造成地表破坏程度分析J.山西焦煤科技,2017,41(1):24-27.4施先盛.大采高工作面两带高度分析研究J.山西焦煤科技,2015(S1):3-5,8.5张玉卓,仲惟林,姚建国.断层影响下地表移动规律的统计和数值模拟研究J.煤炭学报,1989(1):23-31.6尹立明,郭惟嘉,尹增德.断层影响下覆岩破坏规律研究J.山东科技大学学报(自然科学版),2009,28(4):59-63.7于秋鸽.上下盘开采断层滑移失稳诱发地表异常沉陷机理研究J.采矿与安全工程学报,2021,38(1):41-50.8孔宪森,吴侃,陈朋.煤矿开采诱发断层滑移对村庄建筑物影响研究J.煤炭科学技术,2015,43(12):46-50.9陈绍杰,江宁,常西坤,等.采煤塌陷地建设利用关键技术与实践M.北京:科学出版社,2019.10陈绍杰,夏治国,郭惟嘉,等.断层影响下岩体采动灾变响应研究现状与展望J.煤炭科学技术,2018,46(1):20-27.11王诗海,时鹏程,冯宜语.基于离散元方法的煤层覆岩导水裂缝带影响因素研究J.山东煤炭科技,2018(6):149-151,154.12曹慧敏.煤峪口矿地表移动变形量计算及规律分析J.山东煤炭科技,2023,41(3):157-159,163.332023 年第 11 期曹呆军:不同表层厚度下含断层煤体开采覆岩下沉规律