考虑水平荷载的采空区群系统灾变失稳模型
Catastrophic Instability Model of Goaf Group System Considering Horizontal Loads
通讯作者:
收稿日期: 2018-06-10 修回日期: 2018-11-07 网络出版日期: 2019-07-08
Received: 2018-06-10 Revised: 2018-11-07 Online: 2019-07-08
作者简介 About authors
谢学斌(1968-),男,湖南祁东人,教授,从事岩土工程、采矿工程和安全工程方面的科研和教学工作
关键词:
Keywords:
本文引用格式
谢学斌, 熊胡晨, 谢和荣, 李建坤, 田听雨.
XIE Xuebin, XIONG Huchen, XIE Herong, LI Jiankun, TIAN Tingyu.
矿山地下开采形成了大量的采空区群,当这些采空区群达到一定规模后,很可能会产生大面积的灾害性地压活动,进而造成巨大损失。因此,开展采空区群系统的灾变失稳研究对于矿山安全生产具有重要意义。目前,国内外学者采用理论分析[1,2,3]、数值模拟[4]、室内实验[5]和现场监测[6,7]等方法对采空区群进行了大量研究。由于采空区群系统的稳定性主要由顶板和矿柱这2个基本要素决定[8],所以许多研究中将采空区群简化为顶板—矿柱模型。高明仕等[9]将顶板和矿柱均视为弹性体,建立了煤柱失稳冲击破坏的尖点突变模型;张钦礼等[10]将顶板视为坚硬固支梁,将矿柱视为弹性体,建立了采场破坏失稳的尖点突变模型;贺广零等[11]将顶板视为弹性薄板,将矿柱视为温克尔弹性地基,建立了采空区煤柱—顶板系统失稳的尖点突变模型;Ma等[12]将顶板视为弹性薄板,将矿柱视为开尔文体,建立了大型采空区塌陷的力学模型;王金安等[13]将顶板视为弹性薄板,将矿柱视为伯格斯体,建立了采空区顶板破坏的流变力学模型。
综上所述,前人对采空区群顶板—矿柱模型稳定性的研究尚存在以下不足:(1)虽然对采空区群顶板—矿柱模型的研究较多,但大多数只考虑竖直方向的荷载,很少考虑水平方向的荷载,所以研究结果与工程实际存在一定的差异;(2)众多研究中将矿柱视为流变模型,研究矿柱的流变特性,而很少考虑矿柱在长期荷载作用下产生的蠕变损伤;(3)多数研究中使用的突变模型是尖点突变模型,很少使用其他突变模型。
本文所研究采空区群顶板的厚度小于长度和宽度的1/5,因此将其视为弹性薄板,同时将矿柱视为在长期荷载作用下会产生蠕变损伤的Poynting-Thomson体,考虑水平荷载和竖直荷载的共同作用,结合燕尾突变和流变理论建立采空区群顶板—矿柱系统灾变失稳模型,并以广西盘龙铅锌矿为工程实例进行验证。
1 考虑矿柱蠕变损伤的采空区群顶板—矿柱系统的流变模型
1.1 模型构建
图1
矿柱的本构关系[14]为
式中:Dʹ(t)为t时刻的损伤变量,其表达式[15]为
式中:E0为矿柱的初始弹性模量;
1.2 考虑水平荷载的顶板扰度计算
当薄板同时受水平方向和竖直方向的荷载时,其控制方程[16]为
式中:D为顶板的抗弯刚度;ω为顶板的扰度;λ为矿柱有效支撑面积与顶板面积之比;σ为作用在矿柱上的应力;q为顶板自重与上覆岩层荷载共同作用下的均布荷载;h为顶板厚度;∂2F/∂x2为y方向上的主应力;∂2F/∂y2为x方向上的主应力;∂2F/∂x∂y为剪应力,其值为0,假设顶板在水平方向和竖直方向上所受的荷载均为均布荷载,且水平方向的荷载在x方向与y方向的主应力相等,那么其值等于ξq,其中ξ为侧压力系数。
由此式(3)可简化为
式(4)左右两边均对时间求导可得:
由式(5)乘以η/E2再加上式(4)并联立式(1)可得:
设顶板扰度ω的表达式为
式中:ω0为顶板的最大扰度,即顶板中心的扰度。
将式(7)代入式(6)并联立
式中:H为矿柱的初始高度。
对式(8)采用伽辽金法求解可得:
令:
那么式(9)可简化为
令:
那么式(11)可简化为
解得:
1.3 求解积分常数及顶板破坏时间
顶板边缘在破坏之前可以看作是四边固支,此时的边界条件为
近似假设:
根据弹性力学中的薄板理论[17]和式(4)可得顶板的控制方程为
式中:k为矿柱的等效弹性系数,其表达式为
联立式(7)、式(15)和式(16),并采用伽辽金法求解可得顶板初始最大扰度:
当时间t=0时,联立式(10)、式(12)、式(14)、式(15)和式(18)可以求出积分常数C=C1。
矿柱的蠕变损伤及风化导致矿柱有效支撑面积逐渐减小,顶板的中心扰度ω0逐渐增大,此时当ω0满足式(19)时顶板由四边固支变为四边简支[18],将顶板为四边固支时的阶段命名为阶段Ⅰ。联立式(10)、式(12)、式(14)、式(15)和式(19)可以求出此阶段的持续时间t=t1。
顶板为四边简支时的边界条件为
近似假设
此时当ω0满足式(21)时顶板由四边简支变为自由边[18],将顶板为四边简支时的阶段命名为阶段Ⅱ。假设阶段Ⅱ开始的时间t=0,联立式(10)、式(12)、式(14)、式(19)和式(20)可以求出积分常数C=C2,联立式(10)、式(12)、式(14)、式(20)和式(21)可以求出此阶段的持续时间t=t2。
当顶板变为自由边之后,其应变能完全释放,并全部压在矿柱之上,此时矿柱还有一定的支撑能力,顶板并未完全塌陷,将顶板为自由边的阶段命名为阶段Ⅲ。
近似假设
假设阶段Ⅲ开始的时间t=0,联立式(10)、式(12)、式(14)、式(21)和式(22)可以求出积分常数C=C3。
2 采空区群顶板—矿柱系统的燕尾突变模型
矿柱的应力—应变关系[20]可表示为
式中:ε1为应力—应变关系曲线中峰值点的应变。为简化计算,将式(26)进行泰勒展开并截取至四次项,可得:
所以模型的势函数为
式中:
式(28)可简化为
式中:
式(30)就是采空区顶板—矿柱灾变失稳的燕尾突变模型的势函数,其平衡曲面为
燕尾突变的分叉集除了需要满足式(32)外,还需满足
联立式(32)和式(33)可得:
由于通过式(32)和式(33)消去u求出分叉集的解析式非常困难,所以把p1看作常数,把u看作参数,根据式(34)画出一条分叉集曲线,当p1依次取不同的数值时,就能得到完整的燕尾突变分叉集曲面,如图2所示。
图2
当突变处于阶段Ⅰ时,联立式(15)、式(23)、式(29)、式(31)和式(34)可以求出突变时λ=λ1、ω0=ω01,联立式(10)、式(12)、式(14)和式(15)可以求出从阶段Ⅰ开始到发生突变的时间t=t1'。当突变处于阶段Ⅱ时,联立式(20)、式(23)、式(29)、式(31)和式(34)可以求出突变时λ=λ2、ω0=ω02,联立式(10)、式(12)、式(14)和式(20)可以求出从阶段Ⅱ开始到发生突变的时间t=t2'。当突变处于阶段Ⅲ时,联立式(22)、式(23)、式(29)、式(31)和式(34)可以求出突变时λ=λ3、ω0=ω03,联立式(10)、式(12)、式(14)和式(22)可以求出从阶段Ⅲ开始到发生突变的时间t=t3'。
3 采空区群顶板—矿柱系统灾变失稳模型的算法流程
假设模型最开始处于阶段Ⅰ,运用燕尾突变理论计算出模型突变时的顶板下沉量ω0,结合流变理论判断出模型突变的阶段,并计算出突变时矿柱有效支撑面积与顶板面积之比λ,最后运用流变理论计算出模型从开始到发生突变持续的时间t。具体求解过程如图3所示。
图3
4 工程实例
广西盘龙铅锌矿位于广西来宾市武宣县城南方向12 km处,矿区出露地层包括碎屑岩、白云岩、泥灰岩和泥岩。矿体及其顶、底板围岩均为白云岩。矿石矿物以闪锌矿和方铅矿为主,属于薄至极厚急倾斜矿体。矿山采用主、副竖井与多斜井联合开拓运输系统,采矿方法为浅孔留矿法,矿山生产能力为3 000 t/d。在-220 m中段16~19线间存在采空区群,该区形状近似矩形,如图4所示。
图4
表1 盘龙铅锌矿相关参数
Table 1
参数 | 数值 |
---|---|
顶板长度2a/m | 140 |
顶板宽度2b/m | 48 |
顶板抗弯刚度D/(MPa·m-3) | 1.35 |
顶板容重γ1/(kN·m-3) | 28.32 |
顶板厚度h/m | 8.5 |
顶板泊松比ν | 0.23 |
顶板抗拉强度σT/ MPa | 5.25 |
顶板上覆岩层容重γ2/(kN·m-3) | 27 |
顶板上覆岩层厚度hʹ/m | 285 |
矿柱初始高度H/m | 32 |
矿柱初始弹性模量E0/GPa | 27 |
侧压力系数ξ | 0.3 |
矿柱初始有效支撑面积与顶板面积之比λ | 0.5 |
应力应变曲线中峰值点的应变ε1 | 0.0036 |
将数据代入式(34)中,求得ω01无实数解,说明阶段Ⅰ不发生突变;
将数据代入式(14)中,求得C1=0.024,t1=3.6个月;C2=0.029,t2=1.7个月;C3=0.032,t
5 各影响因素变化对系统稳定性的影响分析
图5
图5
各影响因素变化时的ω0-λ曲线
Fig.5
Curves of ω0-λ when various influencing factors change
图6
从图5中可以看出,当a、b值增大时,突变点的ω0逐渐减小,λ逐渐增大,且随a、b值的增大λ的变化越来越慢;当D值增大时,突变点的ω0逐渐增大,λ逐渐减小,且随D值的增大λ的变化越来越快;当h值增大时,突变点的ω0逐渐增大,λ逐渐减小;当q值增大时,突变点的ω0逐渐减小,λ逐渐增大,且随q值的增大λ的呈近似线性变化;当E0值增大时,突变点的ω0不变,λ逐渐减小,且随E0值的增大λ变化越来越慢。通过流变理论分析可知,D、h、E0值越大,采空区群系统稳定的时间越长;a、b、q值越大,采空区群系统稳定的时间越短。
图6为侧压力系数ξ值不同时的ω0-λ曲线,从图中可以看出,当ξ值增大时,突变点的ω0逐渐增大,λ逐渐减小。因此,随着ξ的增大,矿柱的有效支撑面积减小到一个更小值时,系统才会发生灾变失稳,系统稳定的时间将会增加。当采空区群系统的矿柱有效支撑面积大于突变点的矿柱有效支撑面积时,系统处于稳定状态,所以随着ξ的增大,回采更多矿柱才会导致矿柱的有效支撑面积减小到突变值,从而导致系统发生灾变失稳。因此,考虑水平荷载时计算出的系统稳定时间更长,系统突变前能回采的矿柱更多。
6 结论
(1)将顶板简化为弹性薄板,矿柱简化为在长期荷载作用下会产生蠕变损伤的Poynting-Thomson体,结合燕尾突变和流变理论建立了在水平荷载和竖直荷载共同作用下的采空区群顶板—矿柱系统灾变失稳模型。
(2)研究了模型灾变失稳的判据,计算出模型突变时顶板中心的下沉量,矿柱有效支撑面积与顶板面积的比值,以及模型从开始到发生突变的时间,并通过工程实例验证了该模型的有效性和实用性。
(3)研究了水平荷载、顶板尺寸、顶板抗弯刚度、顶板厚度、均布荷载q和矿柱初始弹性模量对采空区群系统稳定性的影响,为采空区群系统稳定性分析提供了新思路和新方法。
参考文献
基于组合预测与变精度粗糙模糊集的采空区稳定性评价
[J].
Stability evaluation of goafs based on combined forecasting and variable precision rough fuzzy set
[J].
基于AHP及模糊综合评判法的采空区稳定性评价
[J].
Stability evaluation for mine goaf areas based on AHP and fuzzy synthetic judgement theory
[J].
Buckling failures of reserved thin pillars under the combined action of in-plane and lateral hydrostatic compressive forces
[J].
复杂采空区稳定性数值模拟分析
[J].
Numerical simulation on stability of complicated goaf
[J].
上覆水平煤层采空区衬砌受荷模型试验研究
[J].
Model tests on loading characteristics of linings in overlying horizontal coal mined-out area
[J].
特大采空区上覆岩层地压与地表塌陷灾害监测研究
[J].
Research on disaster monitoring of overburden ground pressure and surface subsidence in extra-large mined-out area
[J].
基于精密探测的金属矿山采空区群稳定性分析
[J].
Analysis of stability of goaf group in metal mines based on precision detection
[J].
地基开采沉陷及其地表建筑
[M].
Foundation Subsidence Induced by Mining and Surface Building
[M].
煤(矿)柱失稳冲击破坏的突变模型及其应用
[J].
Cusp catastrophic model for instability of coal pillar burst damage and analysis of its application
[J].
基于尖点突变模型的采场顶板—矿柱稳定性分析
[J].
Stability analysis of stope roof-pillar based on cusp catastrophe model
[J].
采空区煤柱—顶板系统失稳的力学分析
[J].
Analysis of instability of coal pillar and stiff roof system
[J].
Study on mechanics and domino effect of large-scale goaf cave-in
[J].
采空区矿柱—顶板体系流变力学模型研究
[J].
Study of rheological mechanical model of pillar-roof system in mined-out area
[J].
岩石力学的弹塑粘性理论基础
[M].
The Elastoplastic Viscous Theories of Rock Mechanics
[M].
采空区顶板—矿柱体系蠕变损伤力学模型研究
[J].
A mechanical model for the roof-jamb creep damage system in the goaf
[J].
弹性板理论
[M].
Theory of Finite Plates
[M].
弹性力学
[M].Xu Zhilun
.Theory of Plasticity
[M].
采空区坚硬顶板破断机理与灾变塌陷研究
[J].
Study of fracture mechanism and catastrophic collapse of strong roof strata above the mined-out area
[J].
基于突变和流变理论的采空区群系统稳定性
[J].
Stability of goaf group system based on catastrophe theory and rheological theory
[J].
坚硬顶板下煤柱岩爆的尖点突变理论分析
[J].
Theoretical analysis of a cusp catastrophe bump of coal pillar under hard rocks
[J].
斜坡演化的燕尾突变模型研究
[J].
Coattail catastrophe model of sliding instability
[J].
灾变理论入门
[M].
An Introduction to Catastrophe Theory
[M].
岩土工程流变力学
[M].
Geotechnical Engineering and Rheological Mechanics
[M].
破碎带软岩流变力学试验与参数辨识研究
[J].
Study of rheological tests and its parameters identification of soft rock in fractured belt
[J].33(Supp
.
盘龙铅锌矿地压控制智能化综合监控技术研究(第1阶段)
[R].
Research on ground pressure control technology and intelligent monitoring system for Panlong lead-zinc mine(the first phase)
[R].
/
〈 |
|
〉 |
