伊旗昊达煤矿围岩活动规律数值模拟探究.docVIP

伊旗昊达煤矿围岩活动规律数值模拟探究 　　[摘要]伊旗昊达煤矿位于内蒙古自治区伊金霍洛旗境内，属浅埋煤层。研究有助于掌握伊旗昊达煤矿综采工作面围岩活动规律、矿山压力显现规律，为煤矿巷道支护设计、工作面开采方案设计提供依据。根据昊达煤矿的地质条件，应用FLAC软件对工作面及回采巷道进行了数值模拟。分析工作面回采过程中，围岩应力与变形的分布特征。保证工作面安全高效生产，使伊旗昊达煤矿在安全生产基础上进步提高煤炭产量，创造更大的经济效益。 [关键词]数值模拟 围岩活动规律 矿压显现规律 支护范围 [中图分类号] U452.1+2[文献码] B[文章编号] 1000-405X（2013）-10-274-2 0引言 20世纪70年代以来，“砌体梁”假说、“关键层”理论和“托板”理论等，极大地丰富了煤层开采过程中岩层运动的研究成果。围岩活动规律研究的主要内容包括：工作面直接顶、老顶的垮落规律、工作面前方支承压力分布规律、工作面开采过程动态围岩应力分布规律、工作面上覆岩层运动规律等。对围岩的活动规律的主要研究方法有：相似材料模拟、数值模拟和现场观测等方法。论文应用数值模拟对围岩的活动规律进行分析研究。迄今为止在研究应力和变形方面有限元方法是最为成熟的方法之一，但是连续性和均匀性假设这两个弱点一直限制着有限元在岩层移动方面的真正运用。岩石内部存在着大量的微孔隙、微裂纹等内部缺陷，并在不同程度上受到非连续面的切割。这些都给岩石变形过程的数值模拟带来了极大的困难。 1昊达煤矿浅埋煤层综采工作面地质概况 昊达煤矿，位于内蒙古自治区鄂尔多斯市伊金霍洛旗新庙镇境内。具体位置在边家壕村以南，丁家圪堵以西一带。伊金霍洛旗昊达煤炭有限责任公司煤矿位于东胜煤田东南部，基本构造形态与东胜煤田一致，为一向南西倾斜的单斜构造，倾角1°～3°，区内无断层及较大褶皱，亦无岩浆岩侵入，地质构造简单，区内无岩浆岩出露。 2 3021工作面围岩活动规律数值模拟 2.1FLAC软件介绍 FLAC（Fast Lagrangian Analysis of Continua）是连续介质快速拉格朗日分析程序。利用显式有限差分方法求解，其基本原理与算法与离散元法相似，它运用节点位移量连续条件，可对连续介质进行大变形分析，基于显式差分法求解运动方程和动力方程。计算中采用莫尔-库仑模型和空单元模型（巷道开挖和工作面回采）。莫尔-库仑（Mohr-Coulomb）屈服准则为： 式中，σ1、 σ3分别是最大主应力和最小主应力，c， φ分别是粘结力和摩擦力。当 fs0时，材料将发生剪切破坏。在通常应力状态下，岩体的抗拉强度很低，因此可根据抗拉强度准则（σ3≥ σ1）判断岩体是否产生拉破坏。 2.2模型建立 根据工作面地质条件和煤岩条件等分别建立以下模型进行计算： 建立3201工作面走向模型图1，模型范围为900 m×83.5 m，共建立39634个单元网格，40100个节点。通过对工作面不断向前回采的数值模拟，研究工作面围岩的应力分布与变形特征。 建立3201工作面倾向模型图2，模型范围为900 m×83.5 m，共建立37234个单元网格，39600个节点。通过对工作面不断向前回采的数值模拟，研究工作面围岩的应力分布与变形特征。 数值计算包括以下内容：（1）工作面矿压显现规律，确定初次垮落步距和周期垮落步距；（2）工作面围岩应力分布规律；（3）工作面上覆岩层活动规律。 2.3模型参数及边界条件 煤岩体物理力学参数的选取，见表1。计算模型边界条件确定如下：（1）模型左右边界施加水平约束，即边界水平位移为零；（2）模型底部边界固定，即底部边界水平、垂直位移均为零；（3）模型顶部为自由边界。 模型上部已达地表，故无需再施加载荷，水平方向施加3 MPa的等效应力。 2.4数值模拟结果分析 （1）3201综采工作面两巷围岩数值计算 对工作面回采巷道进行了数值计算后，运输顺槽与回风顺槽的垂直应力分布与水平应力分布见图3～6。 从图3～6可以看出，3201工作面两巷在锚杆支护下，巷道得到很好的控制，锚杆群形成了组合拱共同支承围岩； 在浅埋深的条件下，由于压力小，巷道支护比较容易，围岩易于控制。可以确定现有的锚杆支护方案是合理的。 （2）3201工作面回采时围岩数值模拟结果分析 随着工作面的不断推进，工作面围岩应力发生着变化，围岩应力分布见图7～11。 （3）工作面走向模型模拟结果分析 工作面向前推进100m时，垂直应力分布情况见图12，煤壁前方受采空区影响出现应力集中现象，在工作面前方0～3m范围内为应力降低区，应力值为2～3MPa，在工作面前方3～2