东沟兰州湾煤矿地表沉陷预测数值分析.pdf

价值工程 0引言 地下煤层的开采势必引起覆岩的运动，进而影响地 表的变形[1]。我国是煤炭资源型国家， 全国一次能源中煤 炭占比达到60左右，且在未来很长一段时间内这种能 源结构占比不会有太大变化[2]， 同时我国人口密集， 地表 土地资源紧缺，煤矿开采与地表土地塌陷形成了显著的 矛盾。为了解决这方面的矛盾我国从绿色发展方面制定 了诸多法律法规，要求煤矿企业进行煤层开采地表塌陷 的相关预测工作，为后期的土地复垦等环保工作奠定基 础[3]。 然而这方面的工作在开展过程中大多评价机构都是 采用传统的岩层移动角及导水裂隙带发育高度计算作为 地表塌陷范围的估计方法[4]， 这种方法的缺点是只能在一 个剖面线上进行平面范围内的预测，对于多个剖面线间 的地表塌陷只能是通过估算进行预测，这将对测量结果 带来一定的误差。 为了提高煤矿开采地表塌陷的预测精度， 本文拟采用 数值计算结合理论计算的方法进行确定， 为开采引起的地 面塌陷预测工作提供新的技术借鉴。 1计算方法及原理 数值分析已成为土木、 水利、 采矿等工程领域进行 要要要要要要要要要要要要要要要要要要要要要要要 作者简介院李永刚 （1974-） ， 男， 新疆乌鲁木齐人， 本科， 采矿工程 师， 主要从事煤矿安全管理及灾害防治相关工作。 的抗爆性能研究[C].中国力学学会.全国强动载效应及防护学术 会议暨复杂介质/结构的动态力学行为创新研究群体学术研讨会 论文集.北京高等教育出版社， 2013 224-233. [2]卢红标， 刘云， 周布奎， 等.高强混凝土板的抗爆性能试验 研究[J].武汉理工大学学报， 2012， 34 （10） 91-95. [3]戴平仁， 黄正祥， 祖旭东， 等.聚脲弹性体 “三明治” 夹层结 构抗爆性能[J].工程塑料应用， 2017， 45 （12） 70-74. [4]潘金龙， 周甲佳， 罗敏.爆炸荷载下 FRP 加固双向板动力 响应数值模拟[J].解放军理工大学学报 （自然科学版） ， 2011， 12 （6） 643-648. [5]周布奎， 王安宝， 杨秀敏， 等.GFRP 加固 RC 双向板抗爆性 能试验研究[J].爆炸与冲击， 2006， 26 （3） 234-239. [6]杜修力， 廖维张.高强钢绞线网-聚合物砂浆加固钢筋混凝 土板的抗爆性能数值分析[J].防灾减灾工程学报， 2010， 30 （6） 595-600. [7]周竞洋， 夏志成， 张建亮， 等.钢板夹泡沫铝组合板抗弯性 能研究[J].工程力学， 2017， 34 （4） 213-220. [8]陈万祥， 卢红标， 候小伟， 等.高强钢筋加强混凝土板抗爆 性能试验研究[J].振动与冲击， 2015， 34 （10） 135-141. [9]刘晶波， 杜义欣， 闫秋实.地下箱型结构在爆炸冲击荷载作 用下的动力反应分析[J].解放军理工大学学报 （自然科学版） ， 2007， 8 （5） 520-524. [10]张秀华， 段忠东， 张春巍.近地空中爆炸作用下钢框架结 构冲击响应[J].地震工程与工程振动， 2009， 29 （4） 70-76. [11]Holmquist T J, Johnson G R. A computational constitutive model for concrete subjected to large strains, high strain rates, and high pressures [A]. 14th International Symposium on Ballistics[C]. Quebec, Canada [s. n.], 1993 591-600. [12]Zhang Yu, Liu Yuan-xue, ZHAO Ji-chang, et al. A new dynamic protection based on the arch axis optimization for underground structures under buried explosion[J]. International Journal of Earth Sciences and Engineering, 2016, 9 （4） 1548- 1558. [13]白金泽.LS-DYNA3D 理论基础与实例分析[M].北京 科学 出版社， 2005 98-105. [14]张斌， 许金余， 李乐， 等.泡沫混凝土回填层在地下复合结 构中的抗爆特性分析[J].四川建筑科学研究， 2010， 36 （6） 135- 138. [15]亨利奇著.爆炸动力学及其应用[M].熊建国， 译.北京 科 学出版社， 1987 123-128. 东沟兰州湾煤矿地表沉陷预测数值分析 Numerical Analysis of Ground Subsidence Prediction in Donggou Lanzhouwan Coal Mine 李永刚 LI Yong-gang （乌鲁木齐市达坂城区发展与改革委员会， 乌鲁木齐 830000） （Development and Re Commission of Daban District， Urumqi， Urumqi 830000， China） 摘要院为确定东沟兰州湾煤矿开采后地表塌陷范围及深度， 采用数值实验的方法， 对 4 号煤层东西翼采区进行了大范围的开采过 程模拟， 获得了地表沉陷的三维变化趋势， 阐明了现有塌陷区的产生过程， 同时与理论计算结果进行了对比， 验证了数值实验的准确 性， 提高了预测工作的精准度， 可为同类工作提供方法借鉴。 Abstract In order to determine the range and depth of surface subsidence after mining in Donggou Lanzhouwan Coal Mine, a large- scale mining process was simulated in the east-west mining area of No.4 coal seam by using numerical experiments, the three-dimensional change trend of surface subsidence was obtained, and the generation process of the existing subsidence area was clarified. At the same time, it was compared with the theoretical calculation results, which verified the accuracy of the numerical experiment and improved the accuracy of the prediction work, which can provide a for similar work. 关键词院煤矿开采； 地表塌陷； 范围预测； 数值模拟 Key words coal mining； surface collapse； range prediction； numerical simulation 中图分类号院TD327文献标识码院A文章编号院1006-4311 （2020） 17-0114-03 114 Value Engineering 岩土工程问题分析的重要手段， 能够非常准确地模拟材 料的塑性破坏和流动， 对土质、 岩石和其它材料的三维 结构受力特性和塑性流动进行模拟具有更为准确和快 速的特点。 本次模拟对象为沉积岩为主的煤矿开采过程， 其中各 岩层的抗拉、 内聚力及摩擦角等力学参数是影响其破断变 形的重要影响因素， 因此， 采用M-C摩尔库伦破坏准则更 为适合， 该准则认为岩石是均质无裂隙的， 岩石内部在法 向应力作用下产生的摩擦力和颗粒间的内聚力形成了抵 抗破坏的力源， 具体可表达为式 （1） [5-6]。 （1） 式中子为拉力；C为内聚力；滓n为法向应力；专为内 摩擦角。 上式为MC准则的基本原理， 为了利用有限元积分实 现对整个模型的计算，需针对单个有限元单元体进行计 算， 此时单元体的摩尔库伦准则可写为 （2） 式中I1为应力张量的第一不变量；J2为应力偏量的 第二不变量；兹滓为应力罗德角。 2实验方案 本次数值分析是对东西翼采区4号煤层开采现状的 模拟， 重点获取该区域沉陷的过程特征， 并与理论计算形 成互补，对地表复垦修复提供区域时空变化信息奠定基 础， 并结合理论计算进行对比分析， 形成相互验证， 提高预 测的精度。 3数值模型的建立及参数选取 根据勘探剖面线及4号煤层底板等高线可知煤层赋 存倾角平均为7-10毅左右，呈现为南高北低的赋存特征， 结合该矿区的水文地质综合柱状图， 以各岩层的平均厚度 作为参考， 并合理简化了地层的空间分布特征， 确定了4 号煤层的上覆岩层的组成分布， 具体见表1。从岩层控制 理论来分析，覆岩中5号、7号、16号岩层分别为主关键 层、 亚关键层和基本顶， 这些岩层强度大、 厚度大， 对其上 覆岩层具有较强的承载作用， 他们的破断规律影响了覆岩 整体的下沉情况。 表 1 覆岩分布特征 序号岩性厚度/m岩性厚度/m 1 2 3 4 6 8 9 10 黄土 散岩 1 散岩 2 泥岩 粉砂岩 粉砂岩 炭质泥岩 中砂岩 10 10 40 85 20 20 10 20 泥岩 细砂岩 粉砂岩 细砂岩 炭质泥岩 泥岩 煤层 地板 10 10 20 20 15 5 10 15 序号 11 12 13 14 15 17 18 19 本次模型为三维立体模型， 根据表1可确定纵向岩层 的空间分布特征， 水平方向的模型范围， 可根据采掘工程 平面图东西翼13个已采工作面的分布情况进行选取， 如 图1（a） 所示， 模型平面范围为2200m*4000m， 完全囊括了 13个采区， 满足塌陷范围的模拟分析需求。同时本次模拟 还根据上述平面范围进行了矿区地表的三维建模， 进一步 提高了数值模拟的可信度， 构建了如图1（b） 所示的三维 地表， 并最终结合地层形成了如图1（c） 所示的三维计算 模型， 该模型长4000m， 宽2200m， 高度600m， 煤层中根据 13个采区的开采顺序以及采区所处空间位置进行了合理 的简化， 将图1（a） 中的1-2，3-4，5-6，7，8-9，10，11-13进 行合并获得了7个模拟开采区域。 图 1 数值模型平面及三维构成情况 （c） 三维模型内外部构成示意图 （a） 采掘工程平面图模型范围示意 （b） 三维地表建立 4地表沉陷现状演变规律分析 模型采用三角网格进行剖分，共产生计算单元 535060个， 可较为精确的现实各个细节， 底部边界设为固 定支撑、 模型四周为辊支撑条件， 顶面为自由面， 更具模型 总体高度设置了最大初始定地应力为16MPa， 地应力侧压 系数取0.5。 115 价值工程 （a） 数值分析影响范围（b） 理论计算影响范围 图 3 地表沉陷范围影响特征 4.1地表位移场演变过程 东西两翼共分为7个模拟开采区域， 开采顺序根据采 掘工程平面图所示时间进行排序， 具体如图1（a） 所示， 通 过模拟开采获得了如图2（a-g）的地表位移场的演变特 征， 首采区开采后地表开始发生沉陷， 由于开采范围及周 边围岩尚未产生大范围的破坏， 因此对地表的影响范围较 小， 第二模拟采区开采后地表影响范围显著扩大， 随着东 翼各个采区的逐一开采， 破坏范围不断扩大， 沉陷深度也 在不断扩大， 沉陷区域中心随采区的移动而变化。西翼区 域目前仅有两个采区开采完毕， 且两个采区为跳采， 相对 各个采区覆岩较为完整， 因此， 地表的沉陷深度较小， 影响 范围也较小， 具体变化过程如图2所示。 4.2最终地表影响范围 如图3（a） 所示， 为本次数值模拟最终的地表沉陷影 响范围的分布特征， 基本与图3（b） 所示的理论计算方式 获得形态和范围相似， 然而， 数值实验的准确度与地层空 间组成特征及岩石力学参数密切相关， 本次所采用的数据 仅为现有资料所提供的合理推测和简化后得到的结果， 因 此， 可能存在一定的误差。 数值实验快速准确， 且能够直观 查看到地表变化的总体过程， 可为土地复垦、 生态修复等 工程提供基础时空演变数据。 5结论 本次数值分析针对大区域的塌陷现状采用 数值和理论进行计算分析，阐明了塌陷区域的 变化规律得到以下结论 淤东翼采区地表沉陷范围在第二模拟采区 开采后开始快速扩大影响范围，且沉陷中心位 置始终处于当前采空区的中部，全部开采完毕 后， 东翼最大下沉量约为7m， 西翼现状为4m左 右， 与理论计算较为接近。 于塌陷现状的平面分布状态与模拟结果的 范围及形态相似，符合一般缓倾斜煤层开采地 表沉陷变化规律。 参考文献院 [1]张向阳， 刘结高， 涂敏.弱黏结顶板厚煤层综放 开采覆岩移动及矿压控制研究[J].采矿与安全工程学 报， 2019 （5） 873-878. [2]谢和平， 吴立新， 郑德志.2025 年中国能源消费 及煤炭需求预测 [J]. 煤炭学报， 2019， 44（07） 1949- 1960. [3]栗欣.绿色矿山文化建设模式的分析[J].中国矿 业， 2019 （7） . [4]童大志， 汪杰.基于数值模拟的红岭铅锌矿岩层 移动带圈定[J].有色金属 矿山部分， 2019 （5） 10-13. [5]柏东良.基于摩尔库伦准则的球壳自重应力弹 塑性分析[J].价值工程， 2014 （32） 326-328. [6]吴晓辉.含软弱夹层砂岩隧道开挖二维数值模 拟分析[J].价值工程， 2019， 38 （06） 144-146. 图 2 地表沉陷位移演变过程 （g） 第七采区 （a） 首采区（b） 第二采区 （c） 第三采区（d） 第四采区 （e） 第五采区（f） 第六采区 116