Research on mining stress distribution law of coal body during open-off cut butt-joining stage of knife handle working face
-
摘要:
以云盖山二矿23304工作面为工程背景,对工作面基本顶的破断形式直接影响刀把式工作面由短工作面过渡到长工作面时切眼对接阶段衔接面周围煤体的应力分布和矿压规律进行了研究。采用理论分析结合3DEC数值模拟的方法分析了云盖山二矿23304刀把式工作面1号切眼工作面与2号切眼对接阶段基本顶不同破断形式和覆岩结构特征,揭示了不同岩梁属性对内应力场范围的影响规律。结果表明,侧向关键块B断裂位置对煤体影响程度由小到大依次为:采空区侧、煤壁前方、煤壁正上方,侧向基本顶断裂位置深入煤体6.21 m,内应力场范围约为12.91~14.43 m,煤体弹性模量对内应力场范围表现为负相关且影响最显著。
Abstract:Taking the 23304 knife handle working face of Yungaishan No.2 Coal Mine as the engineering background, the pressure distribution and ground pressure law of the coal around the joint face during the transition from short working face to long working face are studied when the fracture types of the main roof of the working face directly affects the knife handle working face. Theoretical analysis combined with 3DEC numerical simulation was used to analyze the different fracture types and overlying rock structural characteristics of the main roof at the docking stage between No.1 open-off cut and No.2 open-off cut of Yungaishan No.2 Coal Mine, revealing the influence law of different rock beam properties on the range of internal stress field. The results show that the degree of influence of lateral key block B fracture on the coal body is from small to large: the side of the mining area, in front of the coal wall, and directly above the coal wall, the location of the lateral basic top fracture is 6.21 m deep into the coal body, and the internal stress field range is about 12.91-14.43 m. The elastic modulus of the coal body shows a negative correlation and the most significant influence on the internal stress field range.
-
地下资源开采后,破坏其原有平衡状态,经过覆岩传递到地表,是1个复杂的动态移动过程,涉及地质、采矿、空间等多方面因素,许多学者在开采沉陷动态变形原理、模型、精度等方面进行了研究,并取得了大量成果。开采沉陷动态预测模型主要有动态力学模型和时间函数模型2类,目前主要采用时间函数构建开采沉陷的动态变形过程;常用的时间函数有Knothe时间函数模型、双曲线时间函数模型、Gompertz时间函数模型、logistic 曲线时间函数模型及 Weibull 曲线时间函数模型等。
王悦汉等[1]将动态移动过程分为顶板初次垮落前、顶板垮落未充满采空区、顶板垮落已充满采空区、表土移动4个阶段,采动岩体分为连续介质、拟连续介质、非连续介质3种介质,分阶段建立了动态力学预测模型;吴侃等[2]提出了基于时间序列分析的动态预测模型;崔希民等[3]建立基于Knothe时间函数的概率积分预计方法,根据临界开采尺寸确定时间影响系数;张兵等[4-5]建立了1种新的优化分段Knothe时间函数,并提出了基于实测数据的反算时间函数对比求参法、基于概率积分参数的直接计算法的2种优化分段Knothe时间函数参数求取方法;刘玉成等[6-7]分析了下沉动态过程的Knothe时间函数模型的不足,提出了幂指数的Knothe时间函数,分析了常见地表沉陷时间函数模型的下沉、下沉速度、下沉加速度与时间关系,认为Weibull 曲线时间函数模型能完整地描述地表沉陷的动态过程;孙闯等[8]建立了考虑松散层下沉特点的双因素Knothe时间函数;王军保等[9]借鉴岩石流变力学中非定常流变模型,构建了变时间影响系数的Knothe的地表动态沉陷预测模型;许国胜等[10]回归分析得到地表动态移动变形参数与地质及开采技术参数之间的关系;李春意等[11-12]进行了基于Logistic时间函数、正态分布时间函数的动态沉陷预测;张永胜等[13]提出了以地表点最大下沉速度时刻为分界点,结合偏差改正、生长函数模型的1种分段Weibull时间函数;刘东海[14]根据常村观测资料,指出Weibull时间函数参数不固定,两参数变化规律为距边界变化类偏态分布;JAROSZ等[15]考虑开采产生的空间收敛,建立了双参数Sroka-Schober时间函数;KOWALSKI[16]在Knothe和Sroka-Schober双参数时间函数基础上,提出了广义时间函数;González Nicieza等[17]在Knothe时间函数基础上,引入正态分布函数,建立了正态分布时间函数。
综上,总体认为地表动态下沉量与时间的关系曲线呈“S”形曲线,下沉速度为0→Vmax→0过程,下沉加速度为0→+amax→-amax→0过程。
1. 地表动态移动过程覆岩变化特征
工作面上方的覆岩及地表在井下资源开采后,经历弯曲、垮落、再次压缩等过程,但在工作面不同位置的岩层完整厚度、位移、受力等经历过程不同,根据工作面回采、覆岩受力与垮落情况,将工作面回采结束后的地表分为3个部分,工作面动态分析分区图如图1。
1)竖向整体受压区(Ⅰ区)。处于工作面外侧煤柱区域,工作面开采边界至开采影响边界之间。此部分区域在整个回采期间,覆岩经历了弯曲,但没有垮落、再次压缩等过程,在竖向上是完整的。在整个开采期间竖向方向整体是受压状态。
2)中间区(Ⅱ区)。处于工作面内测靠近煤柱区域,工作面开采边界至覆岩动态完全区边界(覆岩垮落角确定边界与工作面开采边界之间)。此部分区域在整个回采期间,覆岩经历了弯曲、垮落、非完全再次压缩过程等部分或全部过程,主要是由于岩梁的存在导致垮落非充分及垮落岩石受力非充分,在竖向上部分破裂,且离工作面边界(煤柱)越远,破裂高度越大,即完整覆岩厚度减小。在整个回采期间主要经历如下几个阶段:①受到影响到工作面推到下方前,在竖直方向上整体是受到压缩变形,竖向受力增大;②推过后至覆岩垮落前,在竖直方向上整体是受到拉伸变形、竖向受力减小;③垮落后至移动稳定,在竖直方向上整体是受到压缩变形,竖向受力增大至稳定。垮落岩石部分再次受到压缩变形。
3)动态完全区(Ⅲ区)。处于工作面内测远离煤柱区域,覆岩动态完全区(覆岩垮落角确定边界至采空区中心)。此部分区域在整个回采期间,覆岩经历了弯曲、垮落、完整再次压缩过程,在竖向上部分破裂,破裂高度一致。在整个回采期间主要经历如下几个阶段:①受到影响到工作面推到下方(垮落前),在竖直方向上整体是受到压缩变形,竖向受力增大;②到再次受到压缩前,在竖直方向上整体是受到拉伸变形,竖向受力减小;③到移动稳定,在竖直方向上整体是受到压缩变形,竖向受力增大至稳定。
根据开采沉陷岩梁相关理论,地表下沉可简化为各种类型岩梁弯曲,其弯曲下沉值与刚度(弹性模量与惯性矩之积)成反比,其中惯性矩与岩梁高度的三次方成正比。结合图1,在同种情况下,上覆岩层刚度、惯性矩变化规律为:Ⅰ区(整体一致,最大)>Ⅱ区(与离开切眼距离成反比)>Ⅲ区(整体一致,最小)。
目前以时间函数为基础的各种地表动态沉陷预测模型参数主要包括:最大变形值、时间函数(时间函数不同,其参数不同),各类参数与覆岩、开采等因素有关,体现其差异,但目前时间函数都采用相同参数,即在工作面不同位置(煤柱内外)的时间函数参数都一样,这与处在工作面不同位置的受力过程、变形过程、岩层破坏情况及岩石力学性质等明显不符。
2. Weibull时间函数参数变化特点
Weibull时间函数f(t)数学表达式为:
$$ f\left(t\right)=\left\{\begin{array}{cc}km{t}^{m-1}{{\mathrm{e}}}^{-k{t}^{m}}& t\geqslant 0\\ 0& t < 0\end{array}\right. $$ (1) 式中:k、m为参数,正数;t为时间。
据此建立开采沉陷过程中地表下沉时间函数、下沉速度时间函数、下沉加速度时间函数,开采沉陷动态下沉过程曲线如图2。
下沉时间函数W(t):
$$ W\left(t\right)={W}_{0}\left(1-{{\mathrm{e}}}^{-k{t}^{m}}\right) $$ (2) 式中:W0为该点下沉的最终值。
下沉速度时间函数v(t):
$$ v\left(t\right)={W}_{0}km{t}^{m-1}{{\mathrm{e}}}^{-k{t}^{m}} $$ (3) 下沉加速度时间函数a(t):
$$ a\left(t\right)={W}_{0}km\left(m-1-km{t}^{m}\right){t}^{m-2}{{\mathrm{e}}}^{-k{t}^{m}} $$ (4) 根据理论分析及Weibull时间函数上述公式可知:① 当地面达到最大下沉速度时,下沉加速度为0,此时时间tvmax满足${t}_{vmax}=\sqrt[\leftroot{-1}\uproot{8}{{m}}]{{(m-1)}/{km}} $;② 由下沉加速度变化规律可知$m > 1 $;③ k反映了加速时间、加速度大小,k越小,加速度越小,到达最大下沉速度的时间越大;④ m越小,到达最大下沉速度的时间越小,最大下沉速度越小。
根据前面分析可知:工作面不同位置在回采期间,其覆岩移动、破坏、受力、岩梁刚度等不同,则其动态时间函数参数应不同,m、k决定了Weibull时间函数的形态变化过程,因此m、k在工作面不同位置其存在差异,以垮落区边界为原点、指向煤柱方向为正,根据前述岩梁弯曲特征简述及相关岩石力学理论,建立其与地面点距离垮落区域边界距离S的函数,工作面不同位置m、k变化情况图如图3。
在煤柱侧区域(Ⅰ区)m、k为:
$$ m={a}_1 {{\mathrm{e}}}^{-b_1\frac{s}{H}}+c_1 $$ (5) $$ k={a}_2 \frac{{{S}}_{0}}{H}+b_2 $$ (6) 在垮落区外侧至工作面边界(Ⅱ区)m、k为:
$$ {m}={a}_1 {{\mathrm{e}}}^{-b_1\frac{S}{H}}+c_1 $$ (7) $$ {k}={a}_2 \frac{{s}}{{H}}+{b}_2 $$ (8) 在垮落区内侧(Ⅲ区)m、k为:
$$ {m}={a}_1+c_1 $$ (9) $$ {k}={b}_2 $$ (10) 式中: a1、b1、c1、a2、b2为系数;S0为煤柱边界距离;H为工作面采高。
从m、k计算公式及图3中可以看出:
1)在垮落区外侧为负指数关系,以工作面边界为界,此区域分为2部分;Ⅰ区(煤柱范围)m变化较小,基本为常数;Ⅱ区(中间区)m变化较大,在垮落区内侧,m为一常数。
2)以工作面边界、垮落区边界为界,此区域分为3部分;Ⅰ区(煤柱范围)k为常数;Ⅱ区(中间区)m为线性变化;在垮落区内侧,k为一常数。
由于不同位置上覆岩层经历的受力、破坏、变形过程不同,根据开采沉陷相关研究,可采用如下方法确定垮落范围d。① 方法1:采用覆岩力学理论,根据岩梁的相关破裂理论确定垮落边界;② 方法2:采用概率积分法拐点偏移距;③ 方法3:采用顶板覆岩垮落角φ确定,$d={H}/{\mathrm{tan}\;\;\varphi } $。
3. 案例分析
某工作面采用长壁一次性采全高采煤法回采。走向为774 m,倾斜宽为220.7 m。工作面开采深度为470~515 m。根据走向观测数据,采用上面的计算模型,获得的反演模型系数见表1。
表 1 反演模型系数Table 1. Coefficients of inversion model点号 k m R2 Z17 0.271 3 0.294 3 0.469 7 Z18 0.163 1 0.386 6 0.629 3 Z19 0.621 7 0.430 6 0.882 6 Z20 0.131 8 0.401 2 0.945 1 Z21 0.293 1 0.322 1 0.899 5 Z22 0.271 7 0.365 3 0.880 4 Z23 0.135 8 0.530 0 0.925 4 Z24 0.009 9 1.133 0 0.924 2 Z25 0.002 2 1.494 0 0.928 5 Z26 0.001 4 1.545 0 0.952 1 Z27 0.001 7 1.449 0 0.941 8 Z28 0.002 5 1.301 0 0.927 0 Z29 0.002 7 1.261 0 0.939 5 Z30 0.001 0 1.414 0 0.941 2 根据工作面开采情况,以工作面内侧90 m为垮落边界,以此确定各点距离边界的S值,根据表1反演动态时间函数$\varphi \left(t\right) $如下:
$$ \varphi \left(t\right)=1-{{\mathrm{e}}}^{-k{t}^{m}} $$ (11) 在采空区垮落内侧m、k为:
$$ m=1.419\;4, k=0.001\;908 $$ (12) 在采空区垮落外侧m、k为:
$$ m\left(S\right)=1.099\cdot {{\mathrm{e}}}^{-12.97\times \frac{S}{495}}+0.320\;4 \;\;\;(R^{2}=0.921\;6) $$ (13) $$ k=1.517\times \frac{{S}}{495}+0.001\;908 \;\;\;(R^{2}=0.919\text{,}S<0.2H) $$ $$ k=0.305 \qquad\qquad (S>0.2H)$$ (14) 根据工作面开采情况及参数确定情况,分别以开采结束后实测结果作为最终值、概率积分法预测值作为最终值,进行开采第523 d的动态预测,实测结果、概率积分法预测值作为最终值的动态预测与实测对比如图4(图中距离以工作面边界为原点,采空区方向为正)。
以实测结果作为最终值:动态预测中误差为110 mm,为最大下沉值的3.5%;以概率积分法预测结果作为最终值:动态预测中误差为130 mm,为最大下沉值的4.2%。
4. 结 语
1)分析了工作面上方不同位置的覆岩在开采过程中在竖向上的受力、变形、破坏等变化过程,根据煤柱、垮落情况将覆岩与地表划分为竖向整体受压区、中间区、动态完全区3个区域。
2)以Weibull时间函数动态预测模型为例,分析了m、k 2个参数的空间分布规律及特地,构建了2个参数的分区数学模型及选取准则。
3)分别以实测结果、概率积分法预测结果作为工作面开采结束后的最终值,依据本文动态预测模型及参数确定方法,动态预测中误差分别为110、130 mm,为最大下沉值的3.5%、4.2%,证明本模型的可靠性。
-
表 1 岩层物理力学参数
Table 1 Physical and mechanical parameters of rock
岩层 抗拉强度/
MPa抗压强度/
MPa弹性模量/
GPa黏聚力/
MPa内摩擦角/
(°)泊松比 粉砂岩 9.58 126.0 28.65 25.0 41 0.29 细粒砂岩 10.75 119.0 40.12 24.0 36 0.24 中粒砂岩 9.96 109.0 33.00 16.0 29 0.25 二1煤 0.12 1.1 8.00 6.3 26 0.45 炭质泥岩 4.59 38.5 26.21 5.6 32 0.25 粉砂岩 8.75 113.0 32.60 22.0 17 0.20 表 2 岩梁属性对内应力场范围影响计算方案
Table 2 Calculation solution of the influence of rock beam properties on the range of internal stress field
影响因素 S/m $ {h}_{\mathrm{b}} $/m $ {E}_{0} $/GPa $ {h}_{\mathrm{c}} $/m $ {L}_{\mathrm{B}} $/m 1 80 5 4 4 22.5 2 90 7 5 5 24.5 3 100 9 6 6 26.5 4 110 11 7 7 28.5 5 120 13 8 8 30.5 -
[1] 潘锐,蔡毅,黄厚旭,等. 三软煤层回采支承压力分布及支护构件受力规律研究[J]. 采矿与安全工程学报,2021,38(6):1091−1099. PAN Rui, CAI Yi, HUANG Houxu, et al. Distribution of advance abutment pressure and stress law of supporting components in three soft coal seam mining[J]. Journal of Mining & Safety Engineering, 2021, 38(6): 1091−1099.
[2] 高明中,王新丰. 不等长工作面覆岩破坏规律的数值模拟分析[J]. 安徽理工大学学报(自然科学版),2013,33(4):32−36. GAO Mingzhong, WANG Xinfeng. Numerical simulation analysis of overburden failure in working face of unequal length[J]. Journal of Anhui University of Science and Technology (Natural Science), 2013, 33(4): 32−36.
[3] 武东峰. 柏沟矿浅埋变长工作面覆岩运移及压显现规律研究[D]. 淮南: 安徽理工大学, 2019. [4] 王新丰,高明中. 变长工作面采场顶板破断机理的力学模型分析[J]. 中国矿业大学学报,2015,44(1):36−45. WANG Xinfeng, GAO Mingzhong. Mechanical model of fracture mechanism of stope roof for working face with variable length[J]. Journal of China University of Mining & Technology, 2015, 44(1): 36−45.
[5] 刘畅,刘正和,张俊文,等. 工作面长度对覆岩空间结构演化及大采高采场矿压规律的影响[J]. 岩土力学,2018,39(2):691−698. LIU Chang, LIU Zhenghe, ZHANG Junwen, et al. Effect of mining face length on the evolution of spatial structure of overlying strata and the law of underground pressure in large mining height face[J]. Rock and Soil Mechanics, 2018, 39(2): 691−698.
[6] 杨永康,李建胜,康天合,等. 浅埋厚基岩松软顶板综放采场矿压特征工作面长度效应[J]. 岩土工程学报,2012,34(4):709−716. YANG Yongkang, LI Jiansheng, KANG Tianhe, et al. Effect of working face length on underground pressure characteristics by fully-mechanized top-coal caving mining under shallow-buried thick bedrock loose roof[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(4): 709−716.
[7] 刘一扬,史光亮,张光磊,等. 工作面倾向覆岩垮落形态的面长效应及力链拱特征研究[J]. 煤矿安全,2021,52(5):242−249. LIU Yiyang, SHI Guangliang, ZHANG Guanglei, et al. Study on effect of working face length on overburden caving form and force chain arch characteristics[J]. Safety in Coal Mines, 2021, 52(5): 242−249.
[8] 卢邦稳,刘长武,谢辉,等. 不等长工作面覆岩活动及回采巷道采动影响变形规律[J]. 金属矿山,2016(1):34−38. LU Bangwen, LIU Changwu, XIE Hui, et al. Law of overburden strata movement and mining roadway deformation under mining influence in the unequal length of working face[J]. Metal Mine, 2016(1): 34−38.
[9] 高林,刘鹏泽,张盼栋,等. 基本顶断裂形式对倾斜煤层沿空半煤岩巷围岩稳定性影响研究及工程应用[J]. 煤田地质与勘探,2022,50(6):73−80. GAO Lin, LIU Pengze, ZHANG Pandong, et al. Influence of fracture types of main roof on the stability of surrounding rock of the gob-side coal-rock roadway in inclined coal seams and its engineering application[J]. Coal Geoligy and Exploration, 2022, 50(6): 73−80.
[10] 宁静,徐刚,张春会,等. 综放工作面多区支撑顶板的力学模型及破断特征[J]. 煤炭学报,2020,45(10):3418−3426. NING Jing, XU Gang, ZHANG Chunhui, et al. Mechanical model and fracturing characteristics of multi-area supporting roof in fully mechanized mining working face[J]. Journal of China Coal Society, 2020, 45(10): 3418−3426.
[11] 查文华,李雪,华心祝,等. 基本顶断裂位置对窄煤柱护巷的影响及应用[J]. 煤炭学报,2014,39(S2):332−338. ZHA Wenhua, LI Xue, HUA Xinzhu, et al. Impact and application on narrow coal pillar for roadway pro-tecting from fracture position of upper roof[J]. Journal of China Coal Society, 2014, 39(S2): 332−338.
[12] 何文瑞,何富连,陈冬冬,等. 坚硬厚基本顶特厚煤层综放沿空掘巷煤柱宽度与围岩控制[J]. 采矿与安全工程学报,2020,37(2):349−358. HE Wenrui, HE Fulian, CHEN Dongdong, et al. Pillar width and surrounding rock control of gob-side roadway with mechanical caved mining in extra-thick coal seams under hard-thick main roof[J]. Journal of Mining & Safety Engineering, 2020, 37(2): 349−358.
[13] 王中州,李娇娇,秦宾宾,等. 薄基岩破碎顶板工作面回撤通道稳定性及控制研究[J]. 煤矿安全,2021,52(6):230−236. WANG Zhongzhou, LI Jiaojiao, QIN Binbin, et al. Study on stability and control of retracement channel with thin bedrock and broken roof[J]. Safety in Coal Mines, 2021, 52(6): 230−236.
[14] 王红胜. 综空巷道窄帮蠕变特性及稳定性控制技术研究[D]. 徐州: 中国矿业大学, 2011. [15] 邰阳. 坚硬顶板采场定向造缝覆岩三维破断特征及应力场演化规律[D]. 徐州: 中国矿业大学, 2021. [16] 石永奎,宋振骐,王崇革. 软煤层综放工作面沿空掘巷支护设计[J]. 岩土力学,2001,22(4):509−512. SHI Yongkui, SONG Zhenqi, WANG Chongge. The supporting design of road driving along next goaf in soft seam with fully mechanized top-caving[J]. Rock and Soil Mechanics, 2001, 22(4): 509−512.
[17] 王新丰,陆明远. 深部变面长采场顶板破断演化的力学机制分析[J]. 煤田地质与勘探,2022(10):1−15. WANG Xinfeng, LU Mingyuan. Study on mechanical mechanism of roof fracture evolution in deep stope with variable face length[J]. Coal Geology & Exploration, 2022(10): 1−15.
[18] 杨军,王宏宇,王亚军,等. 切顶卸压无煤柱自成巷顶板断裂特征研究[J]. 采矿与安全工程学报,2019,36(6):1137−1144. YANG Jun, WANG Hongyu, WANG Yajun, et al. Fracture characteristics of the roof in gob-side entry retaining with roof cutting and pressure release[J]. Journal of Mining & Safety Engineering, 2019, 36(6): 1137−1144.
[19] 郭金刚,王伟光,何富连,等. 大断面综放沿空巷道基本顶破断结构与围岩稳定性分析[J]. 采矿与安全工程学报,2019,36(3):446−454. GUO Jingang, WANG Weiguang, HE Fulian, et al. Main roof break structure and surrounding stability analysis in gob-side entry with fully-mechanized caving mining[J]. Journal of Mining & Safety Engineering, 2019, 36(3): 446−454.
[20] 宋振骐. 实用矿山压力与控制[M]. 徐州: 中国矿业大学出版社, 1986. [21] 蒋威,鞠文君,汪占领,等. 厚硬基本顶综放沿空巷道覆岩应力分布特征及合理煤柱宽度确定[J]. 采矿与安全工程学报,2020,37(6):1142−1151. JIANG Wei, JU Wenjun, WANG Zhanling, et al. Characteristics of overburden stress distribution and rational pillar width determination of gob-side roadway with thick and hard basic roof in fully m-echanized top coal caving workface[J]. Journal of Mining & Safety Engineering, 2020, 37(6): 1142−1151.
-
期刊类型引用(5)
1. 李鹏. 基于FLAC3D模拟的14103辅运顺槽支护参数设计. 江西煤炭科技. 2025(01): 90-92+96 . 百度学术
2. 马清水,郭瑞,韩伟,宋永明,梁燕翔,刘耀,王佳明. 基于长短期记忆网络的煤巷支护设计参数预测. 能源与环保. 2025(02): 235-240 . 百度学术
3. 李迎富,宋武,杨子仪,赵江梅,曾凡刚. 掘进巷道锚杆支护优化研究. 煤矿机械. 2024(05): 116-118 . 百度学术
4. 陈万辉,郭瑞,韩伟,宋永明,梁燕翔,刘耀,王佳明,许娜,孟波. 煤矿巷道支护方案智能设计研究. 工矿自动化. 2024(08): 76-83+90 . 百度学术
5. 彭启友,缪易辰,潘文,尚江淮. 吸能锚杆数值模拟研究现状与展望. 建筑技术. 2024(20): 2546-2551 . 百度学术
其他类型引用(1)