地下开采引起上覆岩层移动的数值模拟.docVIP

地下开采引起上覆岩层移动数值模拟研究 摘要：本文在充分考虑岩石非均匀性及节理存在的基础上，应用RFPA2D-Strata数值统软件分析采动影响下覆岩移动的动态发展过程。模拟结果显示：开采过程中覆岩发生应力重分布，覆岩中破坏裂隙首先发生在节理等薄弱面，并沿结构弱面向四周发展。每一层岩层的冒落都会经历变形、离层、断裂和冒落这样一个过程。模拟研究对煤矿的安全开采提供了重要的保证。 关键词：地下开采；覆岩移动；数值模拟 中图分类号：TD322 文献标识码：B 文章编号： Numerical simulation study on the movement of overlying strata caused by underground mining Abstract:We applied RFPA2D-Strata numerical software to analyze the dynamic development process of cover rock under mining influence which taking full account of the existences of rock joints.Simulation results show that stress redistribution in the process of mining,In the overlying rock, the failure crack occurs in the joints first and then to the surrounding .Every terrane should undergo the process of deformation, separation, fracture and falling.the existence of cave changed the destruction way of cover rock, Which give an important guarantee for the safety of coal mining. Keywords: underground mining; strata movement; numerical simulation  1引言 我国煤炭资源丰富，而地下开采诱发的岩层移动及地表沉降是一个极其复杂的物理、力学过程，“三下”煤层开采引起的上覆岩层移动给国民经济和生态环境带来巨大的破坏的同时还威胁到矿工的生命安全[1]。长期以来，人们通过现场观测，物理模拟以及一些理论[2-5]分析来认识开采过程中的岩层变形、离层、冒落过程。由于数值模拟方法具有周期短、成本低的特点，该方法已经成为研究岩层移动的有力手段。一般连续介质力学理论和数值分析方都只能做一些简单的定性分析，而采用岩石破裂过程分析软件系统RFPA能反演采动过程中覆岩变形-离层-断裂全过程，可以系统的得出岩层移动破坏规律[6-12]。其基本思想是[13]:1、通过考虑细观单元力学参数的非均质性，模拟岩石宏观变形破坏的非线性行为；2、通过考虑材料破坏后单元的参数的弱化 包括刚度退化 ，模拟材料破坏的非连续和不可逆行为。RFPA2D-Strata主要是面向岩土工程应用分析,可针对地下工程诱发的地表沉陷、岩层移动、巷道及隧道施工引起围岩变形破坏进行数值模拟分[7]。 2数值模拟方案设计 2.1 模拟条件 根据某煤矿煤层的赋存条件建立基本数值模型,设定有关参数进行数值模拟。模型尺寸取水平方向100 m，垂直方向取50 m，埋深150 m ，网个划分为100×200 20000个单元，采用平面应变问题。水平设置位移固定，上部自由边界，考虑岩层自重应力。岩层的物理力学参数如表1所示。 表1 岩层模型参数 Table 1 Parameters of strata model 层序岩性厚度m 弹性模量MPa 抗压强度MPa 自重/Nmm-3 摩擦角顶板3砂岩30 8000 80 2.65 30 顶板2砂岩7.5 5000 50 2.5 35 顶板1砂页岩5 3000 30 2.5 37 煤层2.5 1000 25 1.8 38 底版砂岩5 10000 100 2.65 30 位置均质度弹模抗压强度自重M S -T M S -C 节理50 1000 10 2.0 1.0 20 2.2数值模拟计算 模拟过程中，先对模型施加铅垂方向的自重应力，自重应力根据下公式求得， σz γH 式中:γ—覆岩石的重力密度； H—单元立方体所在的深度。 当自重应力达到3.75 MPa时，利用软件的开挖功能，通过开挖步控制，向右开挖30步，每步向右开挖1 m，共开挖30 m开挖宽度（