基于二分迭代的倾斜条带状工作面开采沉陷预计研究_朱忠华_应用力学.pdfVIP

第 32 卷 第4 期 应 用 力 学 学 报 Vol.32 No.4 2015 年 8 月 CHINESE JOURNAL OF APPLIED MECHANICS Aug. 2015 文章编号：1000- 4939(2015) 04-0664-06 基于二分迭代的倾斜条带状工作面开采沉陷 预计研究 朱忠华 1,2 王李管 1,2 毕林 1,2 肖厚藻1,2 (1 中南大学 资源与安全工程学院 410083 长沙；2 中南大学数字矿山研究中心410083 长沙) 摘要：为方便准确地进行任意倾斜煤层的开采沉陷预计，借鉴变步长数值积分的思想，提出基于 二分迭代的倾斜条带状工作面开采沉陷预计方法，并设计了新的计算模型。具体做法是：将工作 面在埋深方向用二分法划分为多个子工作面，从而将倾斜工作面的开采沉陷预计转化为多个水平 子工作面的开采沉陷预计；对工作面进行二分时细分级数为f ，则要进行f 轮计算；在第 i(i=0 ，1，…， i i f -1)轮计算时，工作面被等分为 2 个子工作面，进行 2 次迭代计算；在每次迭代计算时，均使用 概率积分法，输入地表为上次迭代计算的输出地表，输入工作面为本次计算的子工作面；第 i 轮 迭代计算后，再将每个子工作面在深度方向进行二分，继而进行下轮迭代计算，重复进行，直到 前后两轮计算结果差值的均值和方差分别小于首两轮结果差值的十分之一。在三维可视化环境下 使用 C++编程语言设计并实现了该算法，并从数学上证明了算法的收敛性。算法分析结果表明， 该方法消除了传统概率积分法中倾角和平均深度对倾斜煤层开采沉陷预计的误差，计算模型更加 4 科学。应用实例表明：当f =6 时，前后两轮计算结果差值数据均值为 0.0039m ，方差为2.45×10 ， 具有较好的收敛性；与实测值相比，下沉量和倾斜变形的误差率分别为 0.9%和 0.38% ，误差率较 常规概率积分法小，计算结果更可信。这说明该方法适用于任意倾斜条带状煤层开采沉陷预计， 也可用于非煤地下矿山较规整采场的开采沉陷预计。 关键词：开采沉陷；二分法；变形预计；概率积分法；数值积分 中图分类号：TD672 ；TP391 文献标识码：A DOI ：10.11776/cjam.32.04.B171 民群众的生产、生活。因此，合理准确地进行开采 1 引 言 沉陷预计具有重要的意义。 由于开采沉陷工作的重要性，科研人员和矿 地下矿山开采破坏了采区周围岩体的完整性 山工程技术人员投入了越来越多的时间进行此项 和应力平衡，使上覆岩层断裂、移动，随着开采的 研究[1] 。文献[2]研究了褶皱( 向斜和背斜)煤层的开 进行，这种移动逐渐向上延伸，不可避免地使地表 采沉陷预计，利用已知分布曲线方程来拟合开采煤 产生移动、变形、裂缝、沉降、塌陷，导致不同程 层，采用数值方法求解方程来解决开采曲面分布煤 度的破坏，对生态环境构成严重影响，甚至威胁人 层的沉陷预计问题。文献[2]与本文沉陷预计的“切 基金项目：国家高技术研究发展计划(863 计划)项目(2011AA060407) 收稿日期：2014-12-29 修回日期：2015-07-01 第一作者简介：朱忠华，男，1981 年生，中南大学，博士生；研究方向——数字矿山及数值模拟技术。E-mail ：zzhni