非充分开采地表下沉率变化规律研究.pdfVIP

第七届全国矿山测量学术会议论文集 非充分开采地表下沉率变化规律研究 易四海l，戴华阳’，张连贵2，张广学2，刘吉波‘ 0008 73500 1．中国矿业大学(北京)，北京l3；2．兖州煤业股份有限公司，邹城2 摘要：丈中通过离散元数值计算方法，分别对非充分开采地表下沉率与采动程度、采厚的关系进行了 分析研究，总结了非充分开采地表下沉率的变化规律。该研究加强了矿井深部化非充分开采地表移动规律 的认识，为深部“三下”压煤开采提供了依据。 关键字：非充分开采；地表下沉率；采动程度；采厚；数值计算 1引言 随着矿井的生产发展，开采深度的加大，传统宽度布设的．Il作面开采必然形成非充分开采，大 量的实测资料表明：非充分采动条件F地表移动变形值相对较小，有利于地表建筑物的保护，这为 我们今后解决“三下”压煤问题提供了契机¨1。 地表下沉率受很多因素，如采空区的几何尺寸，开采深度，采厚，覆岩岩性，开采方法。本文 结合某矿四采区特定的地质采矿条件，采用离散元数值计算方法，分别对非充分开采地表下沉率与 采动程度、采厚的关系进行了分析，以研究1F充分开采地表下沉率的变化规律。 2基本地质采矿条件 及砂、砾层等组成，基岩平均厚378m，以岩浆岩、砂岩及粘土岩砂岩互层为主，整体岩性中硬 3模型的建立 3．1离散元数值计算软件2D—Block简介 离散元法足20世纪70年代初兴起的一种数值计算方法，特别适用丁．节理岩体的应力分析，在 采矿工程、隧道工程、边坡工程以及放矿力学等方面都有一定的虑用。它的研究对象主要是岩石等 非连续介质的力学行为，它的基本思想是把不连续分离为刚性元素的集合，使各个刚性元素满足运 动方程，用时步迭代的方法求解各刚性元素的运动方程，继而求得不连续体的整体运动形态。离散 元方法允许单元问的相对运动，不‘定要满足位移连续和变形协调条件，计算速度快，所需存储空 问小，尤其适合求解人位移科t非线形的问题㈠1】。 二维离散元法软件系统2D—Block是我国白行开发的有自主版权的数值计算软件，它具有高度的 集成度，从数据输入到最后形成计算报告或论文都无求于第三方程序的帮助。2D—Block的其它一些 现代化的特性包括可视化和联网等功能¨1。 3．2计算模型的建立 为探求地表下沉率与采动程度、采厚的关系，本文结合某矿四采区的地质采矿条件，设计如图l 所示的离散元模型。在数值模型设计过程中，为便于分析计算，本文将上覆岩层归结为5层，即从上 往下依次为表土层、岩浆岩、砂岩、粘土岩细砂岩互层和I粘土岩砂岩互层：由于煤层倾角较小，模 型中将煤层假定为水平煤层，即煤层倾角00；煤厚取7m，采深568m。计算丁作面划分为18个区段， X 每个区段宽为0．1倍的采深，即56．8m，累计计算工作面K度为56．8 18=1022．4m，根据分析要求， 在计算]：作面两测各延长了600m，则二维离散元计算模型总K2222．4m，累高620m。 3．3模型的块体划分 离散单元块体的划分，以实际岩层的层理、 节理和断层的发育及分布情况为准，兼顾岩层的力 23 哺U}i。伞…∥【I【j削州；j求允汉lC-之蜒 块体的尺寸麻小些，反2，块体的凡可廊大些：(3)庄胃落带和裂隙带内．块体的Jt、J府小些：『n】柞弯 曲I、沉带．块体的尺寸麻人些。假定块体为小变形的．备岩层块体划分尺寸1i理性质如表I所示。 表I模型块体划分戍”埋性质 K． K、 JK． m J‘h jr ?÷层 分割比仲 N／m N／m’ N／Ⅲ’ N巾r_ t。j N／In 表I．鹾