Model tests on landslide dam materials improved by dynamic compaction
-
摘要: 为丰富堰塞坝开发利用理论,指导堰塞坝浅层密实加固处理工程,基于相似定律对堰塞坝料进行了不同能级的室内强夯模型试验,综合采用宏-细观方法分析了强夯能量引起的动应力发展和传播规律、内部变形规律、颗粒破碎和加固效果等。试验结果表明:随着夯击次数增加,由于堰塞坝料密实度提高、颗粒重排列、破碎及填充的叠合效应,夯锤有效加固范围内的动土压力峰值整体呈波动上升趋势。强夯过程中能量逐渐从表层传递到深层,同时强夯能量随深度具有较大的耗散,动应力峰值随深度呈现快速衰减的趋势。由于风化严重,强夯引起了明显的颗粒破碎现象。强夯对松散宽级配堰塞坝料的加固效果明显,不同能级强夯后其锥尖阻力大幅提高,但夯击能超过一定数值时,增大夯击能对加固效果提升有限。基于本次易贡堰塞坝料模型试验,最佳夯击能约为6000 kN·m。Abstract: In order to enrich the development and utilization theories of landslide dams and guide the shallow compaction reinforcement projects, based on the similarity law, the dynamic compaction model tests on landslide dam materials with different energy levels are carried out, and the development and propagation law of dynamic stress caused by dynamic compaction energy, the displacement characteristics as well as the particle breakage and reinforcement effects are analyzed by using the macro-meso-method. The test results show that with the increasing tamping times, the peak value of dynamic pressure within the effective reinforcement range of the rammer exhibits a fluctuating upward trend as a whole due to the increasing compactness of landslide dam materials and the superposition effects of particle breakage, rearrangement and filling. During the process of dynamic compaction, the energy is transferred from the surface one to the deep layer gradually. Meanwhile, the energy dissipates greatly with the depth, and the peak value of dynamic stress decreases rapidly with the depth. Due to severe weathering, obvious particle breakage is caused by the dynamic compaction. The reinforcement effects of the dynamic compaction is obvious for loose- and wide-graded landslide dam materials. The cone tip resistance of the dynamic compaction with different energy levels greatly increases, but the reinforcement effects are limited by increasing the tamping energy when it exceeds a certain value. Based on the model tests on the Yigong landslide dam materials, the best tamping energy is about 6000 kN·m.
-
0. 引言
随着中国煤矿开采深度和生产强度的不断增大,冲击地压成为中国煤矿深部开采面临的主要动力灾害之[1-3]。鄂尔多斯深部矿区面临的冲击地压威胁尤为严峻,截止到2020年,该矿区若干矿井(门克庆、母杜柴登、营盘壕、石拉乌素等)已经进行了初期生产,这些新投建的矿井均具备典型的“厚硬岩层-煤柱”结构,具备大采深(多为600~800 m)、大采高(主采煤层3-1煤厚5~8 m)、宽煤柱(区段煤柱宽20~40 m)和厚硬岩层(主关键层厚度均大于30 m)等特点,该矿区所有矿井的第二临空回采工作面均发生过冲击显现,且工作面临空巷为冲击显现的集中发生地点。
厚硬岩层、宽煤柱是鄂尔多斯深部矿区冲击地压的主要诱发因素,中国学者针对厚硬岩层、宽煤柱的诱冲机理进行了大量研究,朱斯陶等[4]研究了基于冲击地压-矿震协同控制的合理隔离煤柱宽度设计方法;张俊文等[5]认为覆岩转移高应力影响下的沿空巷道,受到低位厚硬顶板运移扰动容易诱发冲击地压;郭文豪等[6]提出多关键层宽煤柱采场的覆岩空间结构具有低位竖“O-X”、高位横“OX”的2种破断形式;王博等[7]研究了双巷掘进条件下的宽区段煤柱应力演化规律;赵毅鑫等[8]认为厚硬顶板条件下的采空区侧向静载、超前静载及扰动动载是沿空巷道的冲击力源;翁明月等[9]研究了多关键层运动对于窄煤柱综采面沿空巷的诱冲机制;解嘉豪等[10]提出了沿空巷坚硬顶板爆破前后的覆岩空间结构特征;白贤栖等[11]以震源群为单位、微震台站为对象构建了反演矩阵进行煤矿矿震震源的机制反演;韩刚等[12]认为近距离坚硬顶板破断产生的强扰动是鄂尔多斯深部矿区沿空巷频发动力显现的主要原因;张俊文等[13]研究了圆形断层煤柱对于四周采场覆岩空间结构的影响;丛利等[14]认为多次采掘扰动将导致宽煤柱承载应力高度集中。
厚硬岩层−煤柱结构失稳诱发型冲击地压是目前鄂尔多斯深部矿区大采高综采面典型的冲击地压模式,目前鄂尔多斯深部矿区的冲击地压研究尚处于起步阶段,导致鄂尔多斯深部矿区面临此类冲击地压时处于束手无策的状态,本文以鄂尔多斯深部矿区某矿30202工作面为工程背景,采用现场实测、理论分析、数值模拟等方法,揭示了深部大采高综采面厚硬岩层-煤柱结构失稳诱冲的发生机制,以期为此类型冲击地压的现场防治提供理论依据。
1. 工程案例
鄂尔多斯深部矿区某矿30202工作面,如图 1,该工作面东部为30201采空区,南部为井田边界,西部为设计30203实体煤工作面,北部为工业广场保护煤柱,区段煤柱宽度为19.5 m。30201工作面生产时间为2016年12月—2018年3月,地表最大沉降量小于1.5 m;30202工作面设计生产时间为2018年4月—2021年1月,临空开采后地表最大沉降量小于2 m。判断工作面均处于非充分采动状态。
30202工作面平均采深640 m,回采3-1煤层平均厚度5.7 m,平均倾角约2°,经鉴定3-1煤层具有强冲击倾向性,顶底板均具有弱冲击倾向性。30201与30202工作面宽度分别为241,268 m,均采用一次采全高综采工艺。根据面内B22钻孔揭露情况(表 1),3-1煤层上覆约40 m范围存在一层厚度大于60 m的厚硬砂岩,具备典型的厚硬岩层-煤柱结构。
表 1 B22钻孔地层参数Table 1 Formation parameters of B22 borehole单位:m 编号 岩性 厚度 编号 岩性 厚度 R15 风积砂 255.1 R7 泥岩 41.1 R14 复合砂岩组 120.5 R6 粗粒砂岩 61.6 R13 泥岩 14.1 R5 泥岩 2.1 R12 复合砂岩组 100.4 R4 2-2中煤 1.3 R11 泥岩 40.1 R3 细粒砂岩 11.8 R10 中粒砂岩 7.2 R2 泥岩 10.3 R9 泥岩 57.9 R1 细粒砂岩 11.9 R8 中粒砂岩 20.2 Coal 3-1煤 5.3 30201工作面回采期间面后留巷(现30202回风巷)共发生9次强矿压显现,造成巷帮位移片帮,部分锚杆锚索脱锚崩断、单体支柱压弯折断,说明仅一侧采空条件下的厚硬岩层-煤柱结构失稳时即可发生冲击显现,30202工作面临空后其回风巷冲击风险水平将进一步增加。为了避免30202工作面临空开采阶段发生厚硬岩层-煤柱结构失稳诱发冲击地压事故,需要提前开展深部大采高综采面厚硬岩层-煤柱结构失稳诱冲机理研究,为此类冲击地压防治设计提供依据,保障厚硬岩层-宽煤柱结构工作面的安全开采。
2. 厚硬岩层-煤柱空间结构演化规律
2.1 采场厚硬岩层-煤柱结构分区特征
厚硬岩层的运动状态控制了采场覆岩空间结构,并对区段煤柱应力分布及变形产生重要影响,本节以30202工作面为例针对厚硬岩层-区段煤柱空间结构特征展开分区分析。
临空工作面前方区段煤柱为单侧采空状态,如图 2(a),煤柱承载覆岩主要由临侧采空区垮落带(F)悬臂岩层结构、裂隙带(S)铰接岩层结构组成,处于悬顶状态的弯曲下沉带(B)岩层传递应力峰值位于本工作面实体煤深处,对区段煤柱影响较小。
临空工作面后方区段煤柱为两侧采空状态,由于宽煤柱的阻隔作用,面内垮落带发育受临侧采空区的影响较小,区段煤柱承载采场垮落带覆岩为两侧对称的悬臂结构。垮落带覆岩运动将产生高位顶板的离层空间,非充分采动条件下,采场覆岩裂隙高度发育上限为采空区短边宽度的1/2[15-16],对应单、双采空区覆岩裂隙发育高度分别为(R6,R7)面内厚硬岩层组、(R8,R9)采场厚硬岩层组,考虑高位厚硬岩层的运动时间效应及区段煤柱的压缩变形程度,工作面后方覆岩结构演化可分为面内裂隙带发育阶段和采场裂隙带扩展阶段。面内裂隙带发育阶段如图 2(b),由于面内厚硬岩层组煤柱端的垂直约束远低于实体煤端,煤柱端为简支受力状态[17],当其破断后裂隙带岩层将形成“W”型破断结构;采场裂隙带扩展阶段如图 2(c),随着面内覆岩裂隙高度持续上移,区段煤柱压缩量不断增加甚至失稳破坏,采场厚硬岩层组离层量持续增加直至破断,形成“V”型破断结构,采场裂隙带的“V”型破断结构将控制采场弯曲下沉带岩层的“V”型运动趋势。
由于宽煤柱的隔离和支撑作用,工作面垮落带覆岩运动无法与临侧采空区联动,裂隙带范围内的厚硬岩层组传递基础应力和破断附加应力是区段煤柱失稳冲击的主要应力来源,考虑高位厚硬岩层组“V”型破断具备滞后性,其破断加载对工作面附近区段煤柱的影响较小,所以本文的研究重点为厚硬岩层组“W”型破断运动对区段煤柱的诱冲效应。
2.2 面内厚硬岩层“W”型破断阶段划分
“W”型结构的形成需要面内厚硬岩层组两次破断运动。如图 3,临空工作面采空初期,面内厚硬岩层发生挠曲下沉,其最大拉应力出现在实体煤侧支撑点且靠近岩层上表面,当该处拉应力超过其强度极限时,厚硬岩层发生初次破断,考虑大采高、垮落带砂岩组碎胀系数小等因素,此阶段面内裂隙带岩层组受采空区矸石的承载力较小,进行理论分析时可忽略不计;初次破断后,厚硬岩层承载的最大弯矩向采空区中部转移且靠近其下表面,达到强度极限后发生二次破断,破碎岩块回转并触矸,最终形成“W”型破断结构。
两次破断前后厚硬岩层空间形态与区段煤柱承载基础应力各不相同,且厚硬岩层破断释放的附加应力增加了区段煤柱的应力集中水平,因此应分别研究面内厚硬岩层组两次破断运动条件下的区段煤柱承载叠加应力及对应的失稳冲击工程判据。
3. 厚硬岩层-煤柱结构失稳诱冲机理
3.1 厚硬岩层-煤柱结构力学模型
(1)初次破断阶段煤柱受力分析
由于工作面实体煤侧范围远大于区段煤柱宽度,对比二者垂直方向的变形和受力,区段煤柱可视为“单轴压缩”而实体煤侧煤体相对稳定,可将区段煤柱简化为等效的弹簧,面内厚硬岩层视为弹性岩梁,得到一端固支、一端受弹簧和临侧采空区铰接岩层加载的厚硬岩层-煤柱结构简化力学模型,如图 4,将区段煤柱对厚硬岩层的承载应力简化为集中载荷,即用集中力F表示,由于面内裂隙带厚硬岩层的破断长度远大于垮落带悬臂长度,所以可视集中力F和相邻采空区铰接岩层加载力F2x均作用于裂隙带铰接点。
a)基础应力分析
初次破断时区段煤柱承载基础应力Pj由垮落带岩层自重R和厚硬岩层挠曲加载应力-F(集中力F的反作用力)叠加组成,
Pj=R−F。 (1) 区段煤柱承载垮落带岩层范围S如图 3所示,其自重R为
R=γSS=H(d+Hcotα)。} (2) 式中:γ为岩层重度;d为区段煤柱宽度;H为垮落带岩层总厚度;α为岩层破裂角。
集中力F计算可从区段煤柱垂直压缩变形角度分析,区段煤柱的垂直压缩变形量Δ满足:
Δ=u+w。 (3) 式中:u为垮落带岩层自重垂直压缩煤柱变形量;w为厚硬岩层挠曲加载垂直压缩煤柱变形量,
u=mR/Ed,w=mF/Ed。} (4) 式中:m为煤层厚度;E为煤层弹性模量。
区段煤柱的垂直压缩变形量Δ同样等于厚硬岩层的最大挠曲下沉量,所以Δ等于厚硬岩层承载载荷q(包括厚硬岩层自重)、临侧采空区铰接岩层加载力F2x和煤柱集中力F单独加载条件下的厚硬岩层挠曲下沉量之和,
Δ=wq+wF2x+wF。 (5) 式中:wq,wF2x,wF分别为厚硬岩层在q,F2x,F单独加载条件下的挠曲下沉量:
wq=q(L+d)4/8E1I,wF2x=F2x(L+d)3/3E1I,wF=−F(L+d)3/3E1I∘} (6) 式中:L为采空区宽度;E1,I为分别为厚硬岩层弹性模量、断面距,I=m13/12,m1为厚硬岩层厚度。由于厚硬岩层破断后形成的铰接块长度远大于其回转下沉量,分析时可认为采空区矸石和铰接点各承载一半载荷,所以铰接岩层加载力F2x可表示为
F2x=ql2/2=q(L−2Hcotα)/2。 (7) 式中:l2为临侧采空区厚硬岩层铰接块长度。
联立式(2)~(7)可得
F=−3E1EId[wq+wF2x−u]Ed(L+d)3+3mE1I。 (8) 设A=ωq+ωF2x-u,B=Ed(L+d)3+3mE1I,联立式(1)、(2)、(8),得初次破断时煤柱承载基础应力Pj的表达函数可表示为
Pj=3E1EIdA+γH(d+Hcotα)BB 。 (9) b)附加应力分析
厚硬岩层初次破断时在均布载荷q(包括厚硬岩层自重)、铰接岩层加载力F2x和煤柱集中力F单独加载条件下的弯曲变形为
yq=qx2[x2−4(L+d)+6(L+d)2]/24E1I,yF2x=F2xx2[3(L+d)−x]/6E1I,yF=−Fx2[3(L+d)−x]/6E1I。} (10) 根据功的互等定理,可得单位宽度厚硬岩层积聚能量U为
U=L+d∫0(yq+yF2x−yF)qdx+F2xwF2x−FwF。 (11) 当厚硬岩层破断时仅有约0.1%~1%的积聚能量转变为震动波能量[18],现场分析厚硬岩层破断释能时U应乘以其破断距a。根据弹性波相关理论[19],厚硬岩层初次破断释能产生的区段煤柱附加应力Pd为
Pd=γζ(2KηaUs−λ/γ)1/122。 (12) 式中:ζ为附加应力传播速度;K为附加应力能量作用指数;η为地震效率;s为附加应力源距离;λ为附加应力波衰减系数。
(2)二次破断阶段煤柱受力分析
厚硬岩层初次破断后演变为一端铰接、一端受煤柱承载和临侧采空区铰接岩层加载的简支梁状态,其结构简化力学模型见图 5。
a)基础应力分析
二次破断时区段煤柱承载基础应力Pj׳由垮落带岩层自重R和厚硬岩层直接加载应力-F׳组成,
P′j=R−F′。 (13) 二次破断时,区段煤柱承载的厚硬岩层范围包括临侧采空区铰接岩层、煤柱承载区厚硬岩层和面内悬顶厚硬岩层宽度的1/2,所以F׳为
F′=q(3L+2d−4Hcotα)/q(3L+2d−4Hcotα)22。 (14) 联立式(2)、(13)、(14)可得二次破断时煤柱承载基础应力P′j的表达函数可为
P′j=γH(d+Hcotα)+q(3L+2d−4Hcotα) /2 。 (15) b)附加应力分析
二次破断时厚硬岩层受均布载荷q(包括厚硬岩层自重)加载时的弯曲变形可表示为
y′q=qx[x3−2(L+d)x2 + (L+d)3]/qx[x3−2(L+d)x2 + (L+d)3]24E1I24E1I。 (16) 二次破断时单位宽度厚硬岩层积聚能量为
U′=L+d∫0y′qqdx。 (17) 最终厚硬岩层二次破断释能产生的区段煤柱附加应力Pd׳可用式(12)进行计算。
3.2 厚硬岩层-煤柱结构失稳诱冲工程判据
(1)工程判据推导
厚硬岩层-煤柱结构的冲击方式与类型由厚硬岩层运动阶段和区段煤柱受力状态共同决定。厚硬岩层为应力加载源,区段煤柱为灾害冲击体,所以应系统性分析厚硬岩层破断与区段煤柱失稳的应力条件。
a)厚硬岩层破断应力条件
厚硬岩层初次破断时最大弯矩M位于实体煤侧固支端,二次破断时最大弯矩M ׳位于岩层中部,可表示为
M=q(L+d)2/2+(F2x−F)(L+d),M′=q(L+d)2/2。} (18) 厚硬岩层初次破断、二次破断时最大弯矩位置的拉应力σt,σ′t为
σt=Mm/2I,σ′t=M′m/2I。} (19) 由于铰接岩层加载力F2x小于集中力F,所以厚硬岩层初次破断M,σt小于二次破断M′,σ′t,初次破断后必然发生二次破断,设厚硬岩层抗拉极限强度为[σt],则厚硬岩层破断的应力条件为
[σt]⩽σt<σ′t。 (20) b)区段煤柱失稳应力条件
从煤柱内部是否存在弹性承载区角度划分:屈服煤柱(无弹性核区)和承载煤柱(有弹性核区)。本文研究的鄂尔多斯深部矿区留设区段煤柱均为承载煤柱,由于孤立煤体的应力均化效应,承载煤柱两侧采空后其承载基础应力逐渐向中部集中,工作面附近的煤柱基础应力分布曲线由马鞍形逐渐演化至平台型,而区段煤柱失稳冲击属于应力问题,工程上通常以煤柱弹性承载区的峰值应力Pa与煤体单轴抗压强度[σc]的比值Ic来判断失稳冲击危险程度,可表示为
Ic=Pa/Pa[σc][σc]。 (21) 考虑承载煤柱的弹塑性分布规律,应力平台型分布的峰值应力Pa为
Pa=max(ˉPj + Pd,ˉP′j + P′d)。 (22) 式中:ˉPj,ˉP′j分别为区段煤柱弹性核区的平均基础应力,ˉPj=Pj/(d−ρ),ˉP′j=P′j/(d−ρ),ρ为煤柱两侧塑性区宽度,近似取巷道半径的(3~5)倍。
煤柱失稳冲击阈值与煤体围压近似呈线性关系,结合承载煤柱内部煤体的围压分布规律,一般情况下煤柱弹性核区边缘的煤体可近似为单向受力状态,承载应力达到[σc]即具备局部失稳冲击风险(Ic≥1)[7];弹性核区中部的煤体处于三向受力状态,承载应力达到3[σc]即具备整体失稳冲击风险(Ic≥3)[15]。
(2)工程案例验算
根据30202工作面及周边地质开采技术条件,取煤柱宽度d=19.5 m,采空区宽度L=268 m;岩层重度γ=25k N/m3;垮落带厚度H=37 m;岩层断裂角α=71.5°;煤层厚度m=5.3 m;煤层弹性模量E=0.57 GPa;厚硬岩层弹性模量E1=27 GPa;厚硬岩层厚度m1=61.6 m;抗拉强度[σt]=6.2 MPa;均布载荷q=2.6 MPa;破断距a=50 m;附加应力传播速度ζ=3500 m/s;附加应力能量作用指数K=1%;地震效率η=1%;附加应力源距离s=265 m(初次)和144 m(二次);附加应力波衰减系数λ=1.5;单轴抗压强度[σc]=35.37 MPa;煤柱塑性区宽度ρ=7.5 m。将以上参量代入式(1)~(22),计算结果如下:
a)厚硬岩层R6稳定性分析
厚硬岩层R6承载的最大拉应力:初次破断阶段σt=9.7 Pa>[σt]、二次破断阶段σ′t=14.6 MPa>[σt],满足厚硬岩层R6破断的应力条件。
b)区段煤柱稳定性分析
不同宽度区段煤柱承载峰值应力计算结果见图 6,可知二次破断条件下的煤柱承载叠加应力均大于初次破断,所以计算强度比Ic时应使用二次破断阶段煤柱承载应力。
常规宽度(≤60 m)的区段煤柱均达到局部失稳冲击的应力条件,同时预防整体失稳型冲击地压的安全煤柱宽度为52 m。案例工作面区段煤柱宽度为19.5 m,承载叠加应力峰值Pa=116.7 MPa,强度比Ic=3.3,满足煤柱整体失稳冲击的应力条件。
综上分析,30202工作面回采时厚硬岩层R6将发生“W”型破断运动,区段煤柱承载叠加应力超过整体冲击危险阈值,具备发生整体失稳冲击的风险。
3.3 厚硬岩层-煤柱结构失稳诱冲发生机制
冲击地压本质是处于高应力状态的冲击倾向性煤岩体达到破坏极限时发生的瞬间失稳并冲击[20],对于厚硬岩层-煤柱结构的冲击地压问题,高应力来源于厚硬岩层传递的基础应力、附加应力,冲击倾向性煤岩体位于区段煤柱的弹性承载区,而破坏极限取决于区段煤柱弹性核的失稳阈值。
深部大采高综采面厚硬岩层-煤柱结构失稳诱冲机理为:区段煤柱承载基础应力来源于面内厚硬岩层组和临侧采空区铰接加载岩层,当面内厚硬岩层达到强度极限时发生“W”型破断并传递附加应力,如图 7,区段煤柱承载的叠加应力超过其破坏阈值后将发生局部甚至整体失稳冲击。
4. 数值模拟及现场实测验证
4.1 空间结构演化数值模拟
根据30202工作面的地质、开采技术条件,采用UDEC数值模拟软件沿采场倾向方向建立模型,数值模型尺寸长×高=1260 m×284 m,对模型的侧向和底部边界加以约束,顶部为自由边界,在模型顶部施加10 MPa均布载荷代替所建模型顶部至地表之间400 m范围的岩层重量。模型岩层的屈服准则采用莫尔-库仑准则,根据矿井地质资料及钻孔实测数据确定岩体、节理的物理及力学参数。采场工作面开采顺序为:30201工作面→30202工作面。
统计分析30202工作面开采后不同计算步数条件下的区段煤柱承载应力变化规律。区段煤柱承载应力峰值随覆岩运动加载整体呈递增状态。如图 8,区段煤柱承载应力相邻两次统计结果(500步)的最大增量值呈现4次突增现象,4次突增峰值的计算步数分别为500,3500,5500,13500步,对应的覆岩裂隙发育如图 9所示,可知4次突增对应的采场覆岩运动阶段分别为Ⅰ—垮落带岩层组破断运动、Ⅱ—裂隙带R6岩层组初次破断、Ⅲ—裂隙带R6岩层组二次破断和Ⅳ-裂隙带R8岩层组破断运动。
计算步数增加对应着覆岩运动滞后工作面的距离增大、对回采工作面冲击危险性的影响程度减小,由于Ⅰ阶段区段煤柱的基础应力较小,Ⅳ阶段滞后回采工作面较远,因此Ⅱ、Ⅲ阶段是厚硬岩层-煤柱结构失稳诱冲的关键阶段,区段煤柱承载基础应力达到了局部冲击的应力阈值,同厚硬岩层R6破断释放的附加应力叠加后将增加区段煤柱的冲击危险等级。数值模拟结果验证了深部大采高综采面厚硬岩层-煤柱空间结构演化规律的合理性。
4.2 采场能量释放微震监测
微震监测结果揭示了深部大采高综采面厚硬岩层-煤柱结构的覆岩破坏特征和能量释放规律,图 10为案例矿井不同边界条件的回采工作面微震监测剖面投影,可知实体煤工作面开采时微震事件能级较低且整体分布为一个“等腰梯形”的断裂结构,如图 10(a),工作面上方0~40 m事件密集分布,判断表 1中岩层R1~R5位于覆岩垮落带,由于采高较大,垮落带范围内的岩层破断后不能形成稳定的铰接结构,将以悬臂状态加载至采空区两侧围岩;工作面上方40~140 m微震事件呈现明显的线性分区分布,判断R6、R7厚硬岩层组位于覆岩裂隙带并且破碎岩体处于铰接状态,加载至采空区矸石和两侧围岩;工作面上方140 m以上区间无微震事件分布,判断R8及以上岩层组位于弯曲下沉带并处于悬顶状态。综上可知,非充分采动条件下,深部大采高综采面覆岩加载模式具有分区性,垮落带基岩为悬臂结构加载,裂隙带基岩为铰接结构加载,弯曲下沉带基岩为悬顶结构加载。
沿空工作面开采时监测微震事件高度发育区间同样为0~140 m,对应岩层高度上限为R6厚硬岩层组,说明R6厚硬岩层组即时加载至沿空工作面,为沿空工作面的诱冲关键层,而R8厚硬岩层组为延时加载,其破坏运动滞后于工作面,对回采作业的冲击危险性影响较小。沿空工作面开采时微震事件能级整体高于实体煤工作面,如图 10(b),高能级微震事件主要集中在工作面中部和区段煤柱附近。由于工作面宽度远大于区段煤柱,所以厚硬岩层R6二次破断的最大弯矩位置可视为采空区中部,即中部区域为厚硬岩层R6的二次破断位置,监测结果表明厚硬岩层R6二次破断释放的能量远高于初次破断;区段煤柱附近的高能级事件集中区说明了该区域具备较高的应力水平,此区域煤岩体运动破坏时将释放积聚的弹性势能,产生高能级微震事件。微震监测结果验证了深部大采高综采面厚硬岩层-煤柱结构力学模型、工程判据的正确性。
5. 厚硬岩层-煤柱结构失稳诱冲防治
5.1 防冲技术体系思路
结合厚硬岩层-煤柱结构失稳诱冲机理可知,该类型冲击地压的防控对象为厚硬岩层破断运动传递的叠加应力(应力源)、区段煤柱内部的弹性核区(冲击体),针对防控对象提出适用于深部大采高综采面厚硬岩层-煤柱结构的冲击地压防控体系,在保证监测预警强度和质量的前提下,防控体系包括承载煤柱强卸载、扰动区域强支护和推进速度强管控。
(1)承载煤柱强卸载
根据前章分析可知,煤柱强卸载即降低厚硬岩层组传递应力(Pa)和弱化区段煤柱弹性核区(d)。降低应力源采用顶板深孔爆破处理临侧采空区铰接岩层组(F2x),弱化冲击体采用煤层大直径钻孔破坏煤柱承载能力([σc])。
(2)扰动区域强支护
被动支护系统与区段煤柱共同承载厚硬岩层组传递的叠加应力,由式(22)可知,加强工作面超前影响区域的被动支护系统,相当于增加厚硬岩层运动时的承载体宽度(d),从而减小传递基础应力(Pj和Pj׳)达到降低冲击危险的目的,同时被动支护可以提高巷道支护体系的完整性和稳定性,减小区段煤柱强卸载导致的围岩大变形问题。
(3)推进速度强管控
控制推进速度保持低速、匀速状态可以减少厚硬岩层组积聚能量的集中释放水平,从而降低厚硬岩层破断释能产生的区段煤柱附加应力(Pd和Pd׳),所以可以根据现场监测结果动态调整工作面临界推进速度。
5.2 现场防治设计与评价
(1)防治措施设计
由于30202工作面现场试验顶板深孔爆破时多次出现5次方高能级微震事件,判断仅单侧采空时19.5 m区段煤柱承载应力处于危险水平,所以仅采用煤层大直径钻孔进行全巷道承载煤柱强卸载方案。针对区段煤柱弹性核区宽度大、煤体硬度高的特点,沿空巷道实施四轮煤层大直径钻孔预卸压措施,见表 2,一至三轮均为全巷道煤层预卸压,四轮针对单、双工作面“见方”来压区域集中预卸压。后续工作面区段煤柱设计时选择屈服宽度方案(≤6 m),从根本上消除煤柱弹性核冲击体。
表 2 煤层卸压钻孔施工方案Table 2. Construction plans for pressure relief drilling of coal seam参数 孔深/间距/(m·m-1) 一轮 二轮 三轮 四轮 煤柱帮 10 ≤2 13 ≤2 13 ≤2 13 ≤1 生产帮 20 ≤2 20 ≤2 22 ≤2 35 ≤1 针对30202回风巷回采扰动区域采取的强支护方案为工作面超前200 m以内采用6架垛式液压支架+160架液压单元支架+锚索梁加强支护,工作面超前200 m以外采用单体液压支柱+锚索梁联合支护。
依据现场反馈确定“支护-卸压”条件下的安全极限推速:地表沉降进入稳沉阶段(第三次双工作面“见方”来压后)之前“见方”来压区间4~6刀/d、非“见方”来压区间6~8刀/d,稳沉阶段之后保持8~10刀/d。
(2)防冲效果评价
将工作面倾向微震事件进行能量-频次统计分析,如图 11,监测统计的能量-频次集中区域均位于实体煤内部,其中频次峰值(214个)位于面内距回风巷100~125 m,能量峰值(1.02×106 J)位于面内距回风巷0~25 m,而煤柱区域的微震能量-频次处于高频低能状态(频次107个、8.45 E×105 J能量),说明区段煤柱附近的应力集中程度低,能量释放均匀,高应力集中区转移至实体煤内部,实现回风巷处于“低应力”、“低扰动”、“低倾向”、“强支护”状态。
6. 结论
(1)探讨了深部大采高综采面厚硬岩层-煤柱结构的空间运动特征,建立了厚硬岩层-煤柱结构力学模型并量化了厚硬岩层“W”型破断条件下的煤柱承载应力演化规律,推导了厚硬岩层-煤柱结构失稳诱冲的工程判据。
(2)揭示了深部大采高综采面厚硬岩层-煤柱结构失稳诱冲机理:区段煤柱承载基础应力来源于面内厚硬岩层组和临侧采空区铰接加载岩层,当面内厚硬岩层达到强度极限时发生“W”型破断并传递附加应力,区段煤柱承载的叠加应力超过其破坏阈值后发生局部甚至整体失稳冲击。
(3)数值模拟、微震监测结果验证了空间结构、力学模型与工程判据的合理性。区段煤柱承载应力随厚硬岩层组的破断运动发生阶段性突增现象,厚硬岩层组“W”型两次破断运动传递的载荷是工作面附近区段煤柱失稳冲击的主要应力来源,二次破断释放的能量远高于初次破断。
(4)设计并实施了深部大采高综采面厚硬岩层-煤柱结构失稳型冲击地压的控制体系,包括承载煤柱强卸载、扰动区域强支护和推进速度强管控,现场评价防冲体系安全有效。
-
表 1 各相似系数
Table 1 Similarity coefficients
相似
系数CG=G原型G模型 CD=D原型D模型 Ch=h原型h模型 CE=E原型E模型 Cγd=γd原型γd模型 数值 1728 12 12 20736 1 表 2 1︰12相似系数下强夯原型与模型参数对照
Table 2 Prototype and model parameters of dynamic compaction under similarity coefficient of 1︰12
序号 参数 锤重/
kg落距/
m夯击能/
(kN·m)1 原型 8640 11.58 2000 模型(半模) 5 0.965 0.0965 2 原型 17280 11.58 4000 模型(半模) 10 0.965 0.192 3 原型 25920 11.58 6000 模型(半模) 15 0.965 0.288 4 原型 34560 11.58 8000 模型(半模) 20 0.965 0.384 表 3 易贡堰塞坝料基本性质
Table 3 Basic properties of landslide dam materials for Yigong landslide dam
天然干密度/
(g·cm-3)天然含水率/
%d50/
mmCu Cs 1.81 6 4.18 20.7 2.0 表 4 不同能级强夯后堰塞坝料的颗粒级配统计
Table 4 Statistics of particle gradation of landslide dam materials after dynamic compaction under different energy levels
粒组 粒径/
mm含量/% 未强夯 2000
kN·m4000
kN·m6000
kN·m8000
kN·mA 10~20 21.7 15.8 15.1 14.9 9.4 B 5~10 23.2 19.9 19.2 18.5 16.6 C 2~5 24.0 24.3 24.8 23.3 24.5 D 1~2 6.0 8.3 8.7 8.8 8.9 E 0.5~1 10.0 12.3 13.5 12.3 15.4 F < 0.5 15.0 19.4 18.5 22.1 25.0 -
[1] SHI Z M, XIONG X, PENG M, et al. Risk assessment and mitigation for the Hongshiyan landslide dam triggered by the 2014 Ludian earthquake in Yunnan, China[J]. Landslides, 2017, 14(1): 269-285. doi: 10.1007/s10346-016-0699-1
[2] WU L Z, DENG H, HUANG R Q, et al. Evolution of lakes created by landslide dams and the role of dam erosion: a case study of the Jiajun landslide on the Dadu River, China[J]. Quaternary International, 2019, 503: 41-50. doi: 10.1016/j.quaint.2018.08.001
[3] 沈光泽, 盛金保, 向衍, 等. 堰塞坝漫顶溃决过程数值模拟及应用[J]. 岩土工程学报, 2018, 40(增刊2): 82-86. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2018S2019.htm SHEN Guangze, SHENG Jinbao, XIANG Yan, et al. Numerical modeling of breach process of landslide dams due to overtopping and its application[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(S2): 82-86. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2018S2019.htm
[4] 赵天龙, 陈生水, 付长静, 等. 堰塞坝泄流槽断面型式离心模型试验研究[J]. 岩土工程学报, 2017, 39(10): 1943-1948. doi: 10.11779/CJGE201710025 ZHAO Tianlong, CHEN Shengshui, FU Changjing, et al. Centrifugal model tests on section form of drainage channel of barrier dams[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(10): 1943-1948. (in Chinese) doi: 10.11779/CJGE201710025
[5] 赵天龙, 陈生水, 王俊杰, 等. 堰塞坝漫顶溃坝离心模型试验研究[J]. 岩土工程学报, 2016, 38(11): 1965-1972. doi: 10.11779/CJGE201611005 ZHAO Tianlong, CHEN Shengshui, WANG Junjie, et al. Centrifugal model tests overtopping failure of barrier dams[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(11): 1965-1972. (in Chinese) doi: 10.11779/CJGE201611005
[6] 何宁, 何斌, 张宗亮, 等. 蓄水初期红石岩堰塞坝混凝土防渗墙变形与受力分析[J]. 岩土工程学报, 2021, 43(6): 1125-1130. doi: 10.11779/CJGE202106016 HE Ning, HE Bin, ZHANG Zongliang, et al. Deformation and stress of concrete cut-off wall of Hongshiyan dammed body at initial stage of water storage[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(6): 1125-1130. (in Chinese) doi: 10.11779/CJGE202106016
[7] 张宗亮, 张天明, 杨再宏, 等. 牛栏江红石岩堰塞湖整治工程[J]. 水力发电, 2016, 42(9): 83-86. https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201609024.htm ZHANG Zongliang, ZHANG Tianming, YANG Zaihong, et al. Remediation Project of hongshiyan dammed lake in niulan river[J]. Water Power, 2016, 42(9): 83-86. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLFD201609024.htm
[8] 何宁, 娄炎, 何斌. 堰塞体的加固与开发利用技术[J]. 中国水利, 2008(16): 26-28. doi: 10.3969/j.issn.1000-1123.2008.16.008 HE Ning, LOU Yan, HE Bin. Technologies of dammed lake strengthen and utilization[J]. China Water Resources, 2008(16): 26-28. (in Chinese) doi: 10.3969/j.issn.1000-1123.2008.16.008
[9] 刘宁, 程尊兰, 崔鹏. 堰塞湖及其风险控制[M]. 北京: 科学出版社, 2013. LIU Ning, CHENG Zunlan, CUI Peng. Dammed Lake and Risk Management[M]. Beijing: Science Press, 2013. (in Chinese)
[10] 何长明, 邹金锋, 李亮. 强夯动应力的量测及现场试验研究[J]. 岩土工程学报, 2007, 29(4): 628-632. doi: 10.3321/j.issn:1000-4548.2007.04.027 HE Changming, ZOU Jinfeng, LI Liang. Field tests on measurement of dynamic stress of dynamic compaction[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(4): 628-632. (in Chinese) doi: 10.3321/j.issn:1000-4548.2007.04.027
[11] 费香泽, 王钊, 周正兵. 黄土强夯的模型试验研究[J]. 岩土力学, 2002, 23(4): 437-441. doi: 10.3969/j.issn.1000-7598.2002.04.009 FEI Xiangze, WANG Zhao, ZHOU Zhengbing. Experimental research of dynamic compaction of loess[J]. Rock and Soil Mechanics, 2002, 23(4): 437-441. (in Chinese) doi: 10.3969/j.issn.1000-7598.2002.04.009
[12] 费香泽, 王钊, 周正兵. 强夯加固深度的试验研究[J]. 四川大学学报(工程科学版), 2002, 34(4): 56-59. doi: 10.3969/j.issn.1009-3087.2002.04.014 FEI Xiangze, WANG Zhao, ZHOU Zhengbing. Model test of improvement depth of dynamic compaction[J]. Journal of Sichuan University (Engineering Science Edition), 2002, 34(4): 56-59. (in Chinese) doi: 10.3969/j.issn.1009-3087.2002.04.014
[13] 贾敏才, 王磊, 周健. 砂性土宏细观强夯加固机制的试验研究[J]. 岩石力学与工程学报, 2009, 28(增刊1): 3282-3290. JIA Mincai, WANG Lei, ZHOU Jian. Experimental research on macro-meso consolidation mechanism of sandy soil with dynamic compaction[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(S1): 3282-3290. (in Chinese)
[14] 熊巨华, 胡斌, 冯世进, 等. 强夯法加固粉土地基室内模型试验研究[J]. 世界地震工程, 2010, 26(增刊1): 241-246. https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC2010S1049.htm XIONG Juhua, HU Bin, FENG Shijin, et al. Model test study on silt foundation improvement by dynamic compaction[J]. World Earthquake Engineering, 2010, 26(S1): 241-246. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SJDC2010S1049.htm
[15] 张清峰, 王东权. 强夯法加固煤矸石地基动应力模型试验研究[J]. 岩土工程学报, 2012, 34(6): 1142-1147. http://www.cgejournal.com/cn/article/id/14617 ZHANG Qingfeng, WANG Dongquan. Model tests on dynamic stress in colliery wastes improved by dynamic compaction[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(6): 1142-1147. (in Chinese) http://www.cgejournal.com/cn/article/id/14617
[16] 张清峰, 王东权. 煤矸石地基在强夯冲击荷载作用下的物理模型试验研究[J]. 岩石力学与工程学报, 2013, 32(5): 1049-1056. ZHANG Qingfeng, WANG Dongquan. Physical model tests of coal gangue foundation under dynamic compaction loads[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(5): 1049-1056. (in Chinese)
[17] 高政国, 杜雨龙, 黄晓波, 等. 碎石填筑场地强夯加固机制及施工工艺[J]. 岩石力学与工程学报, 2013, 32(2): 377-384. GAO Zhengguo, DU Yulong, HUANG Xiaobo, et al. Reinforcement mechanism and construction technology of broken stone fills by dynamic consolidation[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(2): 377-384. (in Chinese)
[18] CHEN L, QIAO L, LI Q. Study on dynamic compaction characteristics of gravelly soils with crushing effect[J]. Soil Dynamics and Earthquake Engineering, 2019, 120: 158-169.
[19] 中华人民共和国建设部. 土的工程分类标准: GB/T 50145—2007[S]. 北京: 中国计划出版社, 2008. Ministry of Construction of the People's Republic of China. Standard for Engineering Classification of Soil: GB/T 50145—2007[S]. Beijing: China Planning Press, 2008. (in Chinese)
[20] 中华人民共和国住房和城乡建设部. 土工试验方法标准: GB/T 50123—2019[S]. 北京: 中国计划出版社, 2019. Ministry of Housing and Urban-Rural Development of the People's Republic of China. Standard for Geotechnical Testing Method: GB/T 50123—2019[S]. Beijing: China Planning Press, 2019. (in Chinese)
[21] 王家辉, 饶锡保, 江洎洧, 等. 振冲碎石桩复合地基抗剪机制的模型试验研究[J]. 岩土力学, 2021, 42(4): 1095-1103. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202104022.htm WANG Jiahui, RAO Xibao, JIANG Jiwei, et al. Model experimental study of the shear mechanism of vibroflotation stone column composite foundation[J]. Rock and Soil Mechanics, 2021, 42(4): 1095-1103. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202104022.htm
[22] 建筑地基处理技术规范: JGJ 79—2012[S]. 北京: 中国建筑工业出版社, 2013. Technical Code for Ground Treatment of Buildings: JGJ 79—2012[S]. Beijing: China Construction Industry Press, 2013. (in Chinese)
[23] 冯世进, 胡斌, 张旭, 等. 强夯参数对夯击效果影响的室内模型试验[J]. 同济大学学报(自然科学版), 2012, 40(8): 1147-1153. FENG Shijin, HU Bin, ZHANG Xu, et al. Model test study on impact parameters' influence on tamping effect[J]. Journal of Tongji University (Natural Science), 2012, 40(8): 1147-1153. (in Chinese)
[24] 占鑫杰, 杨守华, 朱群峰, 等. 堰塞坝改良加固关键技术[R]. 南京: 南京水利科学研究院, 2021: 17-60. ZHAN Xinjie, YANG Shouhua, ZHU Qunfeng, et al Key Technologies for Improvement and Reinforcement of Landslide Dam[R] Nanjing: Nanjing Hydraulic Research Institute, 2021: 17-60. (in Chinese)
-
期刊类型引用(3)
1. 支斌,李树忱,朱颖,田垚,万泽恩. 土压平衡盾构螺旋输送机渣土运移及参数影响. 应用力学学报. 2023(01): 107-115 . 百度学术
2. 莫品强,任志文,林玉祥,褚锋,顾瑞海. 基于孔压静探的滨海相软土抗剪强度解译方法研究. 地下空间与工程学报. 2023(S1): 112-123 . 百度学术
3. 王钰轲,冯爽,钟燕辉,张蓓. 基于集成学习模型的正常固结土抗剪强度指标预测方法. 岩土工程学报. 2023(S2): 183-188 . 本站查看
其他类型引用(7)
-
其他相关附件