煤炭开采引起的地表沉陷是严重的环境地质灾害和安全生产隐患,会给煤炭企业造成巨大的经济损失,已成为影响各大矿区经济发展和社会安定的重要因素(范立民,1995; 郭文兵等,2013)。其中,地表沉陷诱发的地裂缝灾害是地质灾害中一种地面形变灾害,广泛发育于煤矿区内(颜景生,2014; 苏凯峰,2014; 白海波,2002; 张玉宝,刘军波,2013)。在众多煤矿区地裂缝灾害中,尤其以厚硬火成岩覆岩断裂后形成的地裂缝危害性最大。辽宁铁法矿区有5个矿覆岩中侵入了火成岩(于成龙等,2013; 梁越,2006; 刘心广,2010; 朱建建,李娜,2007; 轩大洋等,2012),火成岩厚度在24～400 m。2010年4月,晓南矿西二采区老张庄村附近突然出现了一条大的地裂缝,长1 000 m,宽5.8 m,深约10 m,后经调查,在查采动火成岩覆岩区发现8条大的地裂缝,这些地裂缝长50～1 000 m,宽0.5～8 m,深1～10 m。在煤矿开采形成大量开采地裂缝的同时,还常常会引发滑坡、坍塌等地质灾害,导致矿区生态环境产生不可逆的破坏,并且给井下生产造成重大安全隐患,尤其是当地裂缝与采空区贯通时,常发生漏风、溃水、溃沙等安全事故(康建荣,2008; 邢大韦等,1994; 刘辉等,2013; 许延春等,2009)。因此,加强对火成岩断裂形成地裂缝机理的研究具有十分重要的理论价值和应用价值。

依据铁法矿区物探结果,将火成岩覆岩岩体形状分为楔形、矩形、椭圆形,如图1所示。所有岩体形状均可由这3种组合而成,如有弱面的火成岩(图2a),以及没有弱面火成岩(图2b)。有弱面火成岩,受到均匀荷载时,断裂处就位于弱面处; 没有弱面火成岩,断裂处的形成则与荷载、采空区有关(khachai,Shupletsov,1998; Shapletsov,1990)。

火成岩体在地层中的分布主要有:(1)分布在地表,即火成岩体出露于地表或者火成岩体浅埋于地表下;(2)位于地层中间部位,即火成岩体位于可采煤层和地表的中间部位;(3)位于顶板,即火成岩在煤层上方不远处部位;(4)位于底板,即火成岩体位于开采煤层底板或者距离底板不远处。

采空区与火成岩体的关系可分为:(1)全覆盖火成岩体,即火成岩体位于采空区上方,且火成岩体覆盖整个采空区,如图3a所示;(2)半覆

图1 火成岩形状图

Fig.1 Diagram of igneous shape

(a)wedge;(b)rectangular;(c)oval

图2 火成岩形状组合图

Fig.2 Diagram of igneous shape combination

(a)igneous rock with a weak side;(b)igneous

rocks with no weak side

盖火成岩体,即火成岩体位于采空区上方,但是没有全部覆盖采空区,如图3b所示。平衡拱与采空区之间的分布可分为3类:(1)与平衡拱有两个交点,如图4a所示,由此可形成简支梁模型,根据简支梁模型可准确算出有两个交点时火成岩体的破断距;(2)与平衡拱有一个交点,如图4b所示,形成一端固支的简支梁模型,根据力学知识也能准确得出其火成岩的断裂部位;(3)与平衡拱没有交点,如图4a、b所示,此条件下火成岩与采空区之间关系相对而言较简单,此种情况下的火成岩不论是分布在地表、底层中间还是底板,其与地表的移动关系都较简单,且建模较为方便。

综上,在建立模型时,以均匀层状火成岩为主要研究对象,其他形状(楔形、椭圆形)暂不作为主要研究对象。在此条件下,主要以一端固支且与平衡拱最少有一个交点的情况下为主要研究对象,如图4a、b所示,使得建模更为简便。

图3 火成岩体与采空区关系图

Fig.3 Diagram of the relationship between igneous rock and goaf

(a)full cover igneous rock body;(b)half covered igneous rock body

图4 有两个交点(a)、一个交点(b)和

Fig.4 Diagram of the relationship between igneous rock

and balance arch of two intersection point(a),

a intersection point(b),and

no intersection point(c)

火成岩覆岩断裂处拉应力公式大致分为两种:

{L=((10σ_t+2q-5γh)/(30(q-γh)))^1/2,

L_j=h((4σ_t+q(3μ-2))/q)^1/2.(1)

σ_max=(12μqa⁴b²)/(h²π²(3a⁴+2a²b²+3b⁴)).(2)

式(1)中,L为简支梁模式下发生火成岩覆岩断裂的极限跨距,L_j为悬臂梁模式下火成岩覆岩断裂极限跨距; σ_t为覆岩抗拉强度; q为均布荷载; μ为泊松比; γ为覆岩自身重量; h为覆岩岩层的平均厚度。

式(2)中,σ_max为覆岩的抗拉强度; q为均布荷载,其中q=pg(H-L); μ为泊松比; a为覆岩的悬空长度; b为覆岩的悬空宽度; h为覆岩岩层的平均厚度。

巨厚火成岩具有整体结构,其极限破断距计算公式有3类,其中公式(3)为弯拉破坏模型,式(4)为两端固支梁模型,式(5)为正方形岩板模型:

L=H₀((2σ_max)/q)^1/2.(3)

L=H₀((2σ_max)/(5q₁))^1/2.(4)

b_m=H₀((σ_maxH₀)/(6kq₁))^1/2.(5)

式中,L为岩层极限破断步距; H₀为火成岩厚度; q为岩梁(板)上覆荷载; q₁为硬岩层上的软岩的重量; k为薄板的形状系数; b_m为岩板的极限尺寸。

模型以晓南煤矿82-82'号勘探剖面线地质采矿条件为原型,建立平面应变数学模型,模拟煤层为7号煤层,为缓倾斜煤层,平均倾角为6°,数学模型按水平煤层考虑,采厚0.8～7.3 m,平均采厚为2.8 m,模拟工作面平均采深510 m,工作面长度按7号勘探剖面上开采长度400 m,建立200 m×500 m×623 m的平面模型,数值计算模型如图5所示。模型中均采用Mohr-Coulomb判断岩体的屈服准则破坏,并且均不考虑塑性流动(不考虑剪胀)。由图5的模型可以得到开采厚度为6 m的煤层覆岩的应力分布,如图6～7所示。

图5 数值计算模型

图6 煤层开采采厚6 m水平方向应力云图

Fig.6 Horizontal stress nephogram of the thick of coal seam mining with 6 m

(a)working face advancing 200 m;(b)working face advancing 300 m;(c)working face advancing 400 m

图7 煤层开采采厚6 m垂直应力云图

Fig.7 Vertical stress nephogram of the thick of coal seam mining with 6 m

(a)working face advancing 200 m;(b)working face advancing 300 m;(c)working face advancing 400 m

当采面推进200 m时,切眼处和掘进处前方均出现应力集中区,其垂直应力最大-3.0×10⁷ Pa(-表示压应力),工作面上、下方为减压区,各岩层应力状态由下向上逐渐减小。靠近工作面处存在垂直应力正向分布,说明岩层移动、变形、弯曲过程中存在受拉应力或不受力的情况。工作面推进300 m,采动覆岩应力集中区域范围增大,其垂直应力最大为-3.2×10⁷ Pa,减压区范围也逐渐增加,工作面上方缓慢降低。在工作面不断推进的过程中,在采空区上盘将会出现最大拉应力集中区。顶板拉应力大部分位于顶板中间附近。在采空区顶板中央形成等值应力迹线拱,在未充填情况下靠近顶板中央为拉应力,往上等值迹线拱径逐渐变大,拉应力减小,最终变为压应力,越往上拱径越大,压应力越大。在采空区左右两侧各存在一个高压应力集中区,往远处则压应力又逐渐变小,但变化梯度远远小于矿体附近应力变化梯度。

由图8可见,当工作面推进到火成岩下方时,火成岩所受应力明显集中,由于火成岩厚度较大、强度较高,对上覆岩体的移动和变形起到了控制作用。随着工作面的推进,采空区暴露面积不断增大,火成岩受到的拉应力也不断增大,与岩体抗拉强度逐步接近,当水平应力达到1.243×10⁷ Pa时,火成岩断裂破坏,这是导致地表产生裂缝及塌陷的最主要力源。

图9为采区塑性区分布,从图中可以看出,煤层开采后,塑性区主要分布在煤层顶板岩层上和地表土层下,地表土层主要受剪切破坏; 而煤层顶底板岩层则会出现垮落现象,在采空区顶部出现拉张破坏,这主要是由于采空区外部岩层向采空区垮落,产生拉张应力的作用。

图8 煤层开采水平方向应力云图

图9 煤层开采采厚6 m塑性区分布图

Fig.9 Plastic zone distribution map of the thick of coal seam mining with 6 m

(a)working face advancing 200 m;(b)working face advancing 300 m;(c)working face advancing 400 m

本文针对巨厚火成岩覆岩下开采对地裂缝的形成机理进行探索性研究,采取数值模拟、现场实测资料整理、理论分析相结合的科研手段,对采动地裂缝进行分类,并分析不同类型地裂缝的主要成因。对巨厚火成岩下开采上覆岩层的应力分布情况及火成岩体断裂的部位进行了分析研究,得出以下结论:

(1)火成岩作为关键层控制着地裂缝的发育,火成岩的断裂是形成地裂缝的主要原因。应用平衡拱-梁理论建立火成岩覆岩的力学模型,可以揭示地裂缝发育的力学机理。

(2)采用有限差分软件FLAC3D对晓南矿进行了数值模拟,探讨了巨厚火成岩下煤层开挖后地表的沉陷规律。通过采场上覆岩体的塑性云图的扩展情况分析了覆岩塑性破坏范围,通过主应力云图分析了覆岩内部的应力场分布规律和采区上方地表裂缝及塌陷的形成原因。