上覆岩层性质对采空区地表变形影响研究-防灾科技学院学报.PDF

第 12 卷第 4 期 防 灾 科 技 学 院 学 报 Vol.12,No.4 2010 年 12 月 J. of Institute of Disaster Prevention Dec.2010 上覆岩层性质对采空区地表变形影响研究 1 2 2 孙 超 ，薄景山 ，孙有为 （1 吉林建筑工程学院，吉林 长春 130021； 2 防灾科技学院，河北 三河 065201) 摘 要：采空区引发地表变形问题是环境岩土工程领域的重要研究课题。本文以有限元分析软件ANSYS 为基 础，通过数值模拟方法研究了采空区上方不同性质岩层以及不同上覆岩层组合形式对采空区地表沉降变形的 影响。揭示了上覆岩层性质对采空区地表变形的影响规律。 关键词：采空区；上覆岩层；地表变形 中图分类号：TU456 文献标识码：A 文章编号: 1673-8047(2010) 04-0062-05 0 引言[1-4] 矿体开采后，开采区原有的应力平衡状态被破 30 个单元，采空区宽 40 m ，采厚 5m。开采深度为 坏，采空区上方依次形成塌落带、裂隙带及弯曲带， 45 m 。对其单元的划分采用了局部细化，划分为200 甚至发展到地表产生连续或非连续的地表变形，引 个单元。上覆岩层共分为三层。第一层厚 15 m ，第 发一系列的环境岩土工程问题，给矿区人民的生产 二层厚 15 m ，第三层厚 30 m ，通过改变其力学参 生活及工程建设带来重大影响。采空区地表变形破 数，模拟不同的上覆岩层性质对采空区上地表稳定 坏机理、影响因素、治理措施等方面研究工作日益 的影响。采用的三种岩层物理力学参数见表 1。 受到广大科研人员的重视。 本文采用大型有限元分析软件 ANSYS 进行数 值模拟分析，将采空区地表沉陷研究简化为平面应 变问题进行弹塑性有限元分析。在分析过程中每次 计算分为两个步骤，第一步先计算矿体被开挖前只 有自重存在时岩层的应力状态，从而求得计算模拟 研究区域初始的应力状态，以及由重力应力引起的 初始变形；第二步通过 ANSYS 软件中特有的杀死 图1 几何模型 单元命令，“挖除”采空区单元，从而模拟矿体被 采出的过程，并计算此状态下的应力及变形。最后 表1 上覆岩层物理力学参数 通过 ANSYS 软件后处理中的“载荷组”减操作求 容重 得采空区地表处由于开采引起的真实变形量。通过 岩性 弹性模 内摩擦 粘聚力 角(度) 泊松比 N 数值模拟采空区上覆坚硬程度不同的岩层及不同 （Kpa） 3 序号 量 (Gpa) ( mm ) 的岩层结构组合形式，分析研究了上覆岩层性质的 变化对采空区地表变形的影响。 1 2.62 30.1 1043 0.2 2.4 1 计算模型及模型参数[5-6] 2 3.77 36.1 2500 0.2 2.4 计算采用的几何模型如图 1 所示。该模型水平 5100 0.2 2.4