基于FLAC3D的煤层开采覆岩破坏特征研究（精编版）

1、基于 flac3d 的煤层开采覆岩破坏特征研究【摘 要】本文分析了祁东煤矿7131 工作面工程地质特征，建立工作面的三维地质模型；根据工作面实测岩体力学参数及实际边界条件，利用flac3d软件对煤层开采引起的上覆岩体运移变形破坏特征进行了数值模拟；通过分析由此引起的塑性破坏、应力场变化、位移场量级变化，讨论了工作面上覆岩体的“两带”发育范围、破坏形式及变形特征。【关键词】覆岩破坏；“两带”发育； flac3d ；数值模拟0 引言煤矿生产中煤层开采引起的上覆岩体变形破坏问题一直是与瓦斯、水害并列为三大主要问题，在我国各生产科研单位都对其进行专题研究并取得了很大的进展1，目前确定“两带”形态特征的

2、方法主要有经验公式、物理模拟、现场测试以及数值模拟方法2，常规方法局限颇多；近些年发展起来的数值模拟技术被很好的应用到这一领域中3。本文运用有限差分软件 flac3d 模拟 7131 工作面煤层开采中上覆岩体变形破坏特征，实测数据表明此法在该例“两带”预测中效果良好。1 工程概况祁东煤矿一水平三采区7131 工作面作为模拟对象， 地面标高为+21m， 工作面标高 -393.8m-471.7m；走向长 1630m，倾斜长 160m；煤层结构较简单平均厚3.0m，产状变化不大；倾角平均为12.5。矿井范围内为一单斜构造，采用顶板全部垮落法机械化综采。煤层顶板为复合型顶板，直接顶板为深灰色泥岩或细砂

3、岩，厚度不稳定。老顶为为浅灰色砂岩平均厚16m。其上泥岩类岩层平均厚约20m。粉砂岩平均厚 12m。煤层底板为细粉砂岩。顶底板局部裂隙发育，富水性弱。2 flac3d 数值模拟本次模拟通过有限差分计算模型建立与运行，命令流导出相关特性值图像， 直观反映 71 煤层开采过程中上覆岩体移动、变形、破坏等特征。2.1 模型的建立地质模型的建立遵循综合全面因素、简化又反映实际，尺寸适合的原则。根据工作面特征，模型设计规格为 200（x）*100（y）*80（z） ，开采煤层厚度3.0m，埋深拟定为 400m，煤层底板预留深度15m，上覆岩层 62m，近水平模拟。模型如图1。2.2 模型初始与边界为尽可

4、能的真实体现地下空间开挖情况，更好地进行受力特性模拟和塑性流动分析，需人为封闭单元体，定义左右边界为单约束边界，取u=0，v0（u为 x 方向位移， v 为 y 方向位移）；底部边界为全约束边界，取 u=v=0；顶部为自由边界。顶端施加等效荷载，即自重应力z=?h 式中，为上覆岩层重力密度取25kn/m3；h 为模型顶界距地表的深度，本例中h 为 400，在水平方向上施加由自重力产生的侧向应力 x=?y 式中，为侧压系数，由 =/（1-）确定， 为泊松比，取值 0.3。2.3 本构模型与力学参数选取本次研究确定煤系岩体的本构模型选取较为实用的莫尔 -库仑塑性模型，判别准则采用莫尔-库仑屈服准则

5、描述。模型中各煤岩层的力学参数实验室测定，见表 1。在未进行煤层开挖之前赋予全部属性，得到初始地应力平衡状态，反映原始地层未受扰动影响条件下原位应力分布情况， 见图 2 所示显现了 xx 、yy 、zz 三向主应力均值分布。2.4 模拟计算结果与分析本次计算开挖设计开挖步距为20m，共开挖 5 步计 100m，沿模型走向开挖两端各预留50m 作为保护煤岩柱，两侧煤岩帮各预设35m。选取开挖 40m 和 100m 步长为典型代表，分别从塑性破坏、应力分布、位移变化三个方面来分析说明覆岩破坏随空间和时间变化的动态过程。2.4.1 塑性破坏区模拟结果分析随着工作面回采，煤层顶板相继出现不同程度的破坏

6、，由下而上依次发育拉伸破坏和拉伸裂隙区可用于判断垮落带、裂隙带的发育高度。图3 为工作面沿走向推进的覆岩塑性破坏场模拟结果。顶板破坏类型以剪切与拉伸破坏为主，在煤壁前方和工作面后方塑性破坏形成典型的 “马鞍”型破坏特征； 开挖到 100m时图中显示出了 5 个步距时段塑性破坏区域的特征，对应着“两带”分布区域特征。（a）推进 40m 塑性破坏（b）推进 100m 塑性破坏2.4.2 应力场模拟分析图 4 为工作面推进时的垂直应力场模拟结果。在煤壁的前方和工作面后方产生一定的应力集中的“凸”型应力带。随着工作面推进采空区被填充压实，前方形成超前移动的支撑压力，在每一步开挖中都形成一个类似于“步距

7、”的应力带，至开挖到100m 终止图中记录了 5 个填充步距的应力分布特征。大致可分为未受采动影响区、受采动影响区和受采动影响剧烈区。（a）推进 40m 垂向应力分布（b）推进 100m 垂向应力分布2.4.3 位移场模拟分析图 5 所示为位移变化量级分布图，从图中清楚地看出，顶板岩体随采位移变化特征，位移最大变化量始终在采空区正上方，两帮形成了对上支撑致使侧帮变形较小；随开挖推进而依次向前发展，采空区前方岩体的下沉梯度较后方大，形成非对称性分布随工作面推进距离而增加，而后方经压实使得下沉量逐渐趋于一致。总体上顶板位移变化呈现出移动拱形的特征。（a）开挖 40m 位移量级图（b）开挖 100m 位移量级图3 结论通过数值模拟与分析，结合塑性破坏图与位移变化图，煤层开采上覆岩体塑性破坏总体上表现为“马鞍”型特征，应力的重分布为两端集中分布中间较小的类似拱形特征，位移场的变化体现出阶段性特征，每一步推进都表现为前方下沉较大于后方的不对称型分布。 确定了 7131工作面 “两带” 高度值为 37.539m，基本与现场实测数值基本相符。在资料真实准确的情况下，数值模拟方法是可靠的。【参考文献】1李增学，等 .矿井地质学 m. 北京：煤炭工业出版社， 2009. 2伟韬，武强，李献忠，