邓日键,韩昌良,杨厚强,张有志,刘永强,宋 凯,魏 民
(1.中国矿业大学 矿业工程学院,江苏 徐州 221116;2.中天合创能源有限责任公司 葫芦素煤矿,内蒙古 鄂尔多斯 017000;3.山东兖煤黑豹矿业装备有限公司,山东 济宁 273599)
深井巷道围岩裂隙发育,承载能力下降,变形敏感度显著提高,巷道支护困难,常发生大变形失稳乃至冒顶溃帮事故,极大地制约了矿井的安全高效生产[1]。各种应力扰动是造成这种变形破坏的根源问题[2,3]。当前,有诸多学者针对巷道掘进扰动等进行了一系列的研究,取得了较多的成果。
研究表明,许多学者在深部巷道扰动变形[4-7]、掘进应力分布规律[8-10]等方面开展了大量的研究并取得了新的突破,已经认识到了应力扰动对巷道变形破坏的主控作用[11,12]。然而,现有的研究对于近距离采空区卸压作用下的巷道循环掘进扰动机理研究尚存在欠缺,相应的控制对策未能充分考虑重复性扰动影响,尤其是巷道掘进面邻近顶板在循环掘进扰动下的控制问题还没有得到很好的解决[13,14]。本文结合典型工程案例进行模拟分析,探究在上覆工作面采空的环境下巷道掘进工作面重复性推进对邻近顶板的扰动影响,得出开挖卸载与扰动应力之间的关系以及应力的传递规律,基于扰动应力波传播特征进行区域划分。
内蒙古自治区鄂尔多斯市葫芦素煤矿首采2-1煤层基本开采完毕,相邻工作面间遗留宽为30m的煤柱。2-1煤层下方是2-2煤层,平均层间距为21.74m,属于近距离煤层。在21102、21103工作面回采结束的情况下,于2-2煤层中布置接续22103工作面。拟研究的巷道为2-2煤层中的22103工作面运输巷,埋深670m。运输巷位于上覆煤层21103采空区下,即处于卸压区范围内,下行开采巷道受到上煤层开采卸压与残余应力场的耦合影响与制约[15]。巷道与遗留煤柱之间的水平距离为40m。运输巷空间位置如图1所示。
图1 巷道空间位置
2-2煤层层位稳定,结构简单,厚度变化较大,平均厚度3.2m。煤层一般不含夹矸或局部含有2层夹矸,夹矸岩性以泥岩、砂质泥岩为主。直接顶岩性主要为粉砂岩,局部夹砂质泥岩、细粒砂岩和中粒砂岩。而底板岩性以中粒砂岩为主,其余为粉砂岩、细粒砂岩和砂质泥岩。顶底板强度较低,以软弱-半坚硬岩石为主。
运输巷采用综掘机掘进,每天共推进8m。运输巷断面为矩形,净高3.2m,净宽5.4m,净断面17.3m2。全巷采用锚杆、锚索及锚网联合支护。顶板柔性锚杆直径21.8mm,长度4300mm,间距1400mm,排距由原来的1100mm调整为1400mm,每排4根布置,预紧力不小于200kN。柔性锚杆压编织钢筋网支护,规格为∅6.5mm×5400mm×1400mm,网孔100mm×100mm。煤帮采用左旋无纵筋螺纹钢锚杆与菱形金属网支护,其中锚杆直径22mm,长度2500mm,间距850mm,排距调整方案与柔性锚杆一样,每排各帮4根布置,预紧力矩不小于300N·m。两帮帮肩与帮脚位置的锚杆与水平呈15°进行安装,其他位置锚杆和锚索均垂直于岩面布置。菱形金属网规格为3200×1500mm。巷道支护布置如图2所示。
图2 运输巷支护布置(mm)
本模型采用FLAC3D软件建立,模型尺寸为180m(长)×150m(高)×120m(宽),各个岩层相关参数见表1。巷道为宽5.4m、高3.2m的矩形巷道。模型采用位移固定边界,顶部为应力边界。模型上部边界埋深563.3m,下部边界埋深713.3m。巷道上覆岩体平均密度若按照2500kg/m3计算,模型顶部边界所需施加的垂直应力为14.08MPa,侧压系数λ平均值为1.5。模型如图3所示。
图3 三维数值模型
表1 模型各层的岩石物理力学参数
为了符合实际情况,对21102、21103工作面先后进行开挖,设置计算步数不运算平衡,即采空区垮落但并未压实。22103运输巷推进一定的距离,选取70m,按设计方案支护,运算平衡。随着上部采空区逐渐压实,卸压作用慢慢消失。掘进面与邻近柔性锚杆之间最小水平距离为0.6m。扰动应力监测点布置于顶板同排柔性锚杆中间,即顶板正中央。在距离巷道掘进面51m范围内的同一位置均布点监测,测点共计37处。每次循环开挖步距为200步,通过10次的重复性开挖得到顶板监测点的卸载扰动应力数据,分析得出应力传播特征及规律[16,17]。
上覆煤层21102、21103工作面开采引起采场周围应力重新分布,遗留区段煤柱上的支承压力和采空区垮落的顶板对底板岩层产生开采损伤。底板原岩应力与采动应力经过集中、传递、变化等演化过程,最终形成一定的底板岩层应力分布形式。工作面开采之后的残留采动应力对下行开采的运输巷有着不可忽视的影响,这种应力存在一定的不均衡性[18-20]。下位煤层顶板受上层采动影响,覆岩裂隙较为发育。上位煤层区段煤柱附近高应力集中,而运输巷处于遗留煤柱侧下方、21103采空区正下方,回避了高应力的扰动,这对于运输巷具有一定的卸压作用,变形破坏力源趋减。随着上覆采空区顶板下沉、压实,压力又逐渐恢复,卸压效应慢慢消退[21-23]。采掘环境应力分布如图4所示。
图4 采掘环境应力分布
根据顶板所受扰动应力的曲线变化趋势将顶板支护区分为三个区域:近场扰动区、中场扰动区和远场扰动区。扰动应力与顶板某一点距巷道掘进面之间的初始距离的关系如图5所示。
图5 扰动应力与巷道掘进面之间的最初距离的关系
近场扰动区即最初距掘进面0~7.6m的范围内,曲线呈现急剧下降的形态(图5a)。在距掘进面0~2m内,曲线为平均斜率-1.785的直线,扰动应力随距离的增加而急剧减小。而2~7.6m范围内的曲线为弧形线,扰动应力降低的变化量较0~2m范围内的扰动应力远小得多。在近场扰动区内,掘进面在推进中扰动应力变化量为2.84MPa。相较其他两个扰动区,顶板扰动应力递减变化幅度最大、区域范围最小,且破坏影响能力最大。该区域是扰动应力传播的第一时区。
中场扰动区即最初距巷道掘进面7.6~42.6m的范围内,扰动应力以一种振幅越来越小的波形线递减(图5b)。扰动应力进一步被削弱,其平均变化量为0.025MPa,远远小于近场扰动区的应力变化量。中场扰动区是扰动应力传播的第二区,从近场扰动区传播过来的扰动应力被削弱成相对稳定、破坏较小的应力,随后传递至远场扰动区。
远场扰动区即离巷道掘进面初始距离42.6m之外的区域,曲线以波形线的形式缓缓降低(图5c)。根据曲线趋势走向可判断,曲线最终会相交于横轴,即扰动应力消损为0。应力波的传播趋势仍具有中场扰动区应力波的形态特征,但其振幅小,可视为线性传播,直至损耗消失。故远场扰动应力场范围最广、破坏能力最小,但该区域积累的破坏最多。
距离巷道掘进面最初2m范围之外的卸载扰动应力,在消失前其扰动应力波拟合曲线可被当作是“指数函数”。10次开挖卸载对巷道顶板产生了循环扰动作用,每次的卸载应力与距掘进面初始距离之间存在一定的函数关系,即第一至第十次开挖的函数表达式分别为:σ1=0.9718L1;σ2=0.9717L2;σ3=0.9716L3;σ4=0.9715L4;σ5=0.9714L5;σ6=0.9713L6;σ7=0.9712L7;σ8=0.9711L8;σ9=0.9710L9;σ10=0.9709L10。
由上述10组函数表达式拟合出顶板扰动应力与距掘进面初始距离之间的关系,即指数函数为:
σ=aL
(1)