摘要 急倾斜煤层开采时,地表特别是软弱面露头,常常出现台阶、裂缝等非连续变形现象,对建筑物构成极大的威胁,甚至造成建筑物的严重损坏。本文采用数值计算方法,针对地层弱面条件地表非连续变形的形态和机理进行了研究和讨论。
关键词 层间弱面 非连续变形 机理 有限元
1 有限元计算模型和岩性参数
急倾斜煤层开采时,地层弱面受采动影响常出现台阶状裂缝,对建筑物破坏很大。国内外对此研究不多,更无定量的方法。本文采用有限元方法,就急倾斜地层弱面条件下地表非连续变形的规律、机理等,进行了较深入的研究。
本次计算采用平面弹塑性有限元程序模拟急倾斜倾向主断面开采。计算模型自动剖分,计算结果进行自动处理。专门为本次计算设计的前处理和后处理程序十分方便有效,这也为多模型计算提供了快捷的分析手段。开采煤层倾角α=70°,采厚1.5m,平均采深356m,阶段垂高187.9m。
图1为平面基本计算模型的网格剖分图。

图1 计算网格图
若按弹塑性考虑,所需的岩性参数较多,屈服准则选择与参数的波动会对计算结果产生显著影响,考虑到地表非连续变形的影响因素多,现场又只有少数可供利用的实测岩性参数,本项目按弹性模型进行计算。这样,减少了因参数取值不准所产生的误差,增强了计算方案间的可比性,有利于分别研究各单一因素对地表非连续变形的影响。计算中基本岩性参数取值见附表。
附表 岩性参数
| 基 岩 | 弹性模量E ×103MPa |
泊松比μ | 容量γ ×103N/m3 |
| 基岩 | 3.5 | 0.20 | 25 |
| 煤 | 1.5 | 0.26 | 15 |
| 采空区 | 0.01 | 0.35 | 14 |
| 弱面(层) | 0.1 | 0.35 | 12 |
| 计算按不同弱面位置(弱面至煤层的水平距离Lr)、宽度(d)及弱面介质特性(弹性模量E和泊松比μ)等多种方案进行。 2 地表非连续变形机理分析 2.1 岩层与地表移动形态 通过有限元计算,分析了与弱面有关的多种因素对地表非连续变形的影响程度和方式。下面将综合模型实验和数值分析的成果,不仅从移动形态上而且从力学本质上对地表非连续变形的机理进行深入分析。 图2为三种不同弱面位置或弱面宽度有限元计算模型的岩层与地表移动矢量图,2a:Lr=200m,d=0.2m;2b:Lr=200m,d=10.0m。图2中显示,急倾斜煤层开采倾向主断面上岩层移动,是沿煤层法向自煤层顶板向地面传播和扩展的;主体移动矢量指向采空区并沿煤层法线偏向下山方向;在采空区上山煤柱附近,顶底板移动矢量方向发生偏转形成旋涡区。对照相似模型实验的结果,数值计算与模型实验共同揭示了急倾斜煤层开采岩层的移动方式,即在采动作用下,岩层移动沿法向向上扩展,顶板沿层面向采空区弯曲。基于连续介质理论,有限元法虽未直接揭示顶板的垮落形态,但从移动量上也反映了实验所显示的移动方式。这是两种研究手段方法上的差异,分析时需要分辨结果上的差别。
图2 岩层与地表移动矢量图 2.2 采动应力分布与非连续变形机理
图3 岩体采动应力分布等值线图
图4 无弱面时相应位置两侧岩体的应力分布曲线
图5 有弱面时弱面两侧岩层的应力分布曲线 由于弱面介质岩性松软,容易屈服发生塑性变形甚至失去承载力。弱面水平厚度相对较小,介质释放水平应力的能力较弱,故弱面上、下盘的水平正应力变化甚小;弱面沿竖向的尺寸较大,故介质有足够的空间释放竖向正应力和剪应力,因而弱面上、下盘竖向正应力和剪应力变化很大。开采引起的竖向正应力和剪应力在弱面上释放了74%和90%左右,使得弱面间产生不同步相对移动,因而引起地表的非连续变形。 作者单位:戴华阳 王金庄 (中国矿业大学北京校区) |
煤炭网版权与免责声明:
凡本网注明"来源:煤炭网www.coal.com.cn "的所有文字、图片和音视频稿件,版权均为"煤炭网www.coal.com.cn "独家所有,任何媒体、网站或个人在转载使用时必须注明"来源:煤炭网www.coal.com.cn ",违反者本网将依法追究责任。
本网转载并注明其他来源的稿件,是本着为读者传递更多信息的目的,并不意味着本网赞同其观点或证实其内容的真实性。其他媒体、网站或个人从本网转载使用时,必须保留本网注明的稿件来源,禁止擅自篡改稿件来源,并自负版权等法律责任。违反者本网也将依法追究责任。 如本网转载稿件涉及版权等问题,请作者在两周内尽快来电或来函联系。
网站技术运营:北京真石数字科技股份有限公司、喀什中煤远大供应链管理有限公司、喀什煤网数字科技有限公司
总部地址:北京市丰台区总部基地航丰路中航荣丰1层
京ICP备18023690号-1 京公网安备 11010602010109号
