基于下沉等影响原理的倾斜煤层开采沉陷预计方法及应用_田国灿.pdf

基于下沉等影响原理的倾斜煤层开采沉陷预计方法及应用 田国灿1， 2，徐乃忠1， 2，高超1， 2 1. 煤炭科学研究总院 开采研究分院，北京 100013; 2. 天地科技股份有限公司 开采设计事业部，北京 100013 ［ 摘 要］通过对以往开采沉陷的预计算法进行分析，采用下沉等影响原理建立沉陷预计的模 型，设计适用于任意形状工作面的开采沉陷预计程序。该程序采用 C语言编写，进一步优化了工作 面各角点坐标和最终计算结果的处理，并通过对计算区域采用三角形划分法的对比得知本模型可以进 一步消除倾斜煤层在开采下沉计算中由采深的变化所带来的误差，为倾斜煤层的预计提供了更为准确 的计算方法。 ［ 关键词］下沉等影响原理; 倾斜煤层; 开采沉陷; C语言; 概率积分法 ［ 中图分类号］ TD327［ 文献标识码］ A［ 文章编号］ 1006- 6225 201702- 0065- 04 Prediction and Application of Incline Coal Seam Mining Subsidence Based on Subsidence Principle TIAN Guo- can1， 2，XU Nai- zhong1， 2，GAO Chao1， 2 1. Mining Institute，China Coal Research Institute，Beijing 100013，China; 2. Coal Mining ＆ Designing Department，Tiandi Science ＆ Technology Co. ，Ltd. ，Beijing 100013，China Abstract The past mining subsidence prediction was analyzed，subsidence prediction model was built based on subsidence principle，mining subsidence prediction program that could be applied any shape working face was designed. The program was designed by C language，every corner point coordinate of working face and the final results were optimized，the error that induced by mining depth variation in incline coal seam mining subsidence prediction could be eliminated in the model by triangle division com- pared in calculated area，and then more accurate prediction could be provided for incline coal seam prediction. Key words subsidence principle; incline coal seam; mining subsidence; C language; probability integrity ［ 收稿日期］ 2016－10－27［ DOI］ 10. 13532/j. cnki. cn11－3677/td. 2017. 02. 017 ［ 基金项目］ 国家自然科学基金 51604151 ［ 作者简介］ 田国灿 1991－ ，男，河北衡水人，在读硕士研究生，研究方向为开采沉陷预计。 ［ 引用格式］ 田国灿，徐乃忠，高超 . 基于下沉等影响原理的倾斜煤层开采沉陷预计方法及应用 ［ J］ ． 煤矿开采，2017，22 2 65－68. 地表沉陷预计是矿山开采沉陷学的重要内容之 一，对开采沉陷的理论研究和生产实践均具有重要 意义 ［1 ］。在理论上，利用预计的结果可以定量地 分析研究受开采影响的地表在时间上和空间上的分 布规律。为了提高开采沉陷预计的准确性，有必要 对预计方法所采用的理论模型、预计参数与地质采 矿条件之间的定量、定性关系进行深入地研究，从 而能更进一步地加深对开采沉陷基本规律的认识。 在生产实践上，利用开采沉陷预计结果可以对 “三下” 建筑物下、铁路下和水体下的开采活 动进行技术指导 ［2 ］。 因此，在借鉴已有理论体系的基础上，应用下 沉等影响原理和概率积分法的基本理论，开发一个 能对任意形状的倾斜工作面进行处理并具有足够高 精度的开采沉陷预计程序十分必要 ［3 ］。 1程序设计技术路线 在开发环境方面，C作为目前主流的编程语 言之一，既继承了 C语言的强大功能，又有 VB 的快速编写方式，对于构建系统界面具有无可比拟 的优势。既可以保证快速、稳定地构建应用程序， 又能在进行开发的时候调用现有的、无论是什么语 言编写的 COM 对象，具有极强的交互性。因此决 定采用 C语言来完成系统框架的构建 ［4 ］。系统整 体构架见图 1。 图 1系统整体构架 如图 1 所示，程序的运行首先需要获取工作面 56 第 22 卷 第 2 期 总第 135 期 2017 年 4 月 煤矿开采 COAL MINING TECHNOLOGY Vol. 22No. 2 Series No. 135 April2017 ChaoXing 的相关信息，工作面各角点的坐标不再需要手工录 入，可通过对 CAD 组件的调用实现直接对任意形 状的多边形进行角点坐标捕捉。 地质采矿条件参数的输入包含煤层采高、工作 面最大采深、工作面煤层倾斜方向、煤层倾角、下 沉系数、主要影响角正切、水平移动系数、开采影 响传播角以及上山、下山、开切眼边界、停采线边 界的拐点偏移距。 初始参数设置完之后对所有角点进行沿重心的 坐标旋转，使走向方向沿 X 轴正方向，倾向沿 Y 轴正方向，以便于随后求得沿走向或倾向的地表移 动变形值。然后通过设计好的预计模型对所要预计 的影响范围、任意剖面或任意点的移动变形进行计 算，并对计算结果的坐标点再次进行坐标旋转，以 达到与原坐标系相匹配的目的。最后对数据进行后 处理，生成剖面图或者等值线图。 2等影响原理的算法实现 目前我国应用广泛且较为成熟的预计方法主要 是概率积分法，但由于工作面多边形的不规则性， 使得二重积分的上下限难以界定，所以大部分预计 模型都是将工作面的不规则多边形切割划分为规则 的矩形来进行叠加运算。即使有国内学者提出从工 作面坐标原点出发 如图 2 中的 1 角点分别连 接工作面的其他各点，将整个区域划分为若干个三 角形并对每个三角形进行二重积分的三角形划分 法，以此来确定积分上下限并消除由矩形叠加算法 引起的计算误差 ［5 ］，但该解决方法依然没有对地 表沉陷主要影响半径 r 值的变化因素进行考虑，只 是采用平均采深值将 r 值设定为一个常量。但在实 际开采过程中，倾斜煤层条件下的采深是随着工作 面位置的移动在不断变化的，工作面的倾角越大， 采用平均采深算法的误差也就越大。因此，有必要 建立一个可以对任意形状工作面以及任意煤层倾角 进行沉陷预计的模型。 图 2三角形划分法示意 2. 1等影响工作面及相关变量的转换 由沉陷预计概率积分法计算公式，地下坐标点 U x，y处的单元开采引起地表坐标 P x0， y0 处的下沉值为 w x0，y0 1 r2 e －π x－x02 y－y02 r2 1 r H/tanβ 2 式中，r 为主要影响半径; H 为采空区坐标点 x， y处的采深; tanβ 为地表沉陷预计参数中的主要 影响角正切。 在开采区域较大，并且煤层倾角不为 0 的情况 下，采深 H 值随着煤层的倾向和走向在不断地变 化。若图 3 中 A 点为最大采深位置，其坐标为 xa，ya ，工作面 ABCD 上任意一点 U x，y的 采深值则需要由煤层倾角、煤层倾向以及最大采深 来计算求得。 图 3实际工作面与投影工作面 若煤层倾角为 α ，煤层倾向为 θ ，最大采深为 H0，则开采区域任意点 x，y的采深 H 为 H H0－ ya－ y cotθ xa － x 1 cot2 槡 θ tanα 3 若将所得的各角点坐标沿逆时针方向旋转 90－θ ，则新生成的工作面倾向沿 y 轴正方向， 此时 θ 值为 90，进而消去式 3中的参数 cotθ ， 同时规范了坐标系，并为之后的预计计算做准备。 旋转之后的任意点 x，y的采深 H 为 H H0－ ya－ y tanα 4 将地表点 P x0，y0与工作面上任意点 U x，y相连接并延长至最大点采深所在的水平面 M 得到点 Q x1，y1 ，并由图 4 可知等比关系 y1 － y 0 y － y0 H0 H0－ ya－ y tanα 5 同理 x1 － x 0 x － x0 H0 H0－ ya－ y tanα 6 如果用 r1来表示投影至水平面的等价工作面 开采主要影响半径，则 r1 H 0/tanβ 7 66 总第 135 期煤矿开采2017 年第 2 期 ChaoXing 图 4工作面点与投影点的坐标关系 将式 4 、 5 、 6 、 7带入 1可得 w x0，y0 1 r2 e －π x1 －x 0 2 y 1 －y 0 2 r12 8 可以看出投影面上的 Q 点 x1，y1与工作 面上的 U 点 x，y对地表点 x0，y0 的影响 是等价的 ［6 ］。然后对公式 1 进行面积分得到整 个工作面的开采对任意地表点 x0，y0产生的 下沉值 w x0，y0 w max D 1 r2 e －π x－x02 y－y02 r2dxdy 9 式中， wmax为该工作面开采后地表最大下沉值。 对式 9进行换元并积分可得 w x0，y0 wmax D1 1 r2 e －π x1 －x 0 2 y 1 －y 0 2 r12dx1 dx dx1 dy1 dy dy1 10 同时由式 5 、 6可得 dy dy1 H0－ ya－ y tanα H0－ y1 － y 0 tanα 11 dx dx1 H0－ ya－ y tanα H0 12 最终将式 11和式 12带入式 10 中 可得 w1 x0，y0 w max D1 1 ［ 1 － y1 － y 0 tanα H0 ］ r1 2 e －π x1 －x 0 2 y 1 －y 0 2 r12dx1dy1 13 该式即为转换之后的等价投影面的下沉预计面 积分公式，D1为地表任意点 x0，y0 与工作面 各角点链接之后投影至最大采深平面的多边形 如图 3 中多边形 ABC1D1 ，在新生成的多边形中 主要影响半径值为常量，同时每一个要计算的地表 点都有与之对应的各不相同的等价投影多边形。 2. 2线积分转化 投影面转换完成 D1之后，D1依然是不规则图 形，对面积分的上下限进行判定并求取二重积分的 算法比较复杂，因此将式 13转化为线积分对 其进行求解。 由数学分析可知，当区域 D 是一个单连通域， 且函数 P X，Y ，Q X，Y在 D 上具有一阶连 续偏导数的条件下，则下面等式成立  D  Q  X － P Y dxdy ∮ S Pdx Qdy 这就是线积分和面积分转换的格林公式，其中 S 为区域 D 的边界曲线 ［7－9 ］。 由式 13可知 r1为常量的情况下对 x 进行 积分相对比较简单，又由高斯误差函数 erf x 2 槡π ∫ x 0e －x 2dx 可知  D1 1 ［ 1 － y1 － y 0 tanα H0 ］ r1 2 e －π x1 －x 0 2 y 1 －y 0 2 r12dx1dy1 ∑ n i 1∮ si 1 2r1［ 1 － y1 － y 0 tanα H0 ］ e －π y1 －y 0 2 r12 erf槡π x1 － x 0 r1 dy1 其中区域 D1是投影在平面 M 上的由直线 S1，S2Sn所围成的闭合区域，高斯误差函数 可以单独进行计算。至此将 r 值不断变化的二重积 分，转换为 r 值为常量的定积分，在计算过程中， 只需将地表点所对应的投影面的各角点坐标计算出 来，然后对每一条边进行上式的线积分计算，即可 得到与原开采区域 D 相同的沉陷预计结果。该方 法使工作面边界得以确定，并考虑了开采影响范围 r 值不断变化对地表下沉产生的影响，同时将二重 积分等价转化为定积分减少了由计算机计算二次积 分所带来的误差 ［10 ］。对该式进行相应的求导变换， 即可得到倾斜、曲率、水平变形等地表移动变形函 数，由于这些函数均为定积分，不存在上下限难以 界定的问题，只需沿每条边逐一积分即可得出计算 结果，水平移动的计算函数也可通过上述的换元方 法最终划为相应的线积分形式，在此不再赘述。 3程序结构 本程序主要由如下几个模块所组成 76 田国灿等 基于下沉等影响原理的倾斜煤层开采沉陷预计方法及应用2017 年第 2 期 ChaoXing 书书书 准备模块录入地质条件参数以及工作面各角 点坐标。 地表影响范围计算模块自动计算开采沉陷影 响范围，以 2 倍的主要影响半径为扩展边界，并确 定要计算的地表点坐标。 坐标转换模块将坐标系统转换为以煤层走向 沿 x 轴正方向，倾向 上山方向沿 y 轴正方向。 变形计算模块结合编写完成的预计模型，将 经过左边转换后的各个参数导入，得到任意地表点 的移动变形值。 输出模块生成影响范围内各地表点的地表移 动变形值表格并绘制等值线图。 程序结构框图见图 5。 图 5程序结构 4应用实例 设定 一 不 规 则 形 状 工 作 面，走 向 长 度 约 1000m，倾向长度约 450m，下山方向采深 500m， 煤层采出厚度 3000mm，煤层倾角 13，倾向沿 y 轴正方向，下沉系数为 0. 75，主要影响角正切为 2. 0，最大下沉角 90，水平移动系数 0. 3，拐点偏 移距均为 0。将其输入到编写好的沉陷预计程序界 面中，并在 CAD 中获取工作面各角点坐标值。然 后用建立好的模型对其进行下沉计算，可得下沉等 值线图形，其中多段线为设定的工作面边界，最大 下沉值为 2146mm。采用将工作面划分为多个三角 形的三角形剖分法获得的下沉等值线图形，最大下 沉值为 2095mm。 将 2 个等值线图进行叠加对比得到图 6。可以 看出当考虑煤层倾角所带来的采深变化因素时，相 对于三角形划分算法，开采损害影响的边界范围变 得较小，但是最大下沉值较三角形划分法的计算结 果增大了 51mm。 当下沉值越大的时候，下沉等影响算法的影响 范围就越为明显地大于三角形剖分法的影响范围。 图 6 2 种算法模型计算的下沉等值线对比 煤层倾角越大，这种变化就会越为明显。这种整体 影响范围变小但是一定范围内下沉值变大的变化特 点符合倾斜煤层所引起的实际下沉情况，同时也满 足主要影响角不变，采深越小影响范围越小的理 论。 5结论 1通过下沉等影响原理，将主要影响半径 化为常量，推导出倾斜煤层开采沉陷预计的精确计 算公式。 2将二重积分形式转换为线积分形式，确 定了积分的上下限，进一步简化了积分的过程。 3通过编程语言建立了开采沉陷预计模型， 并与传统的三角形划分模型进行对比，分析了 2 种 模型得出的不同预计结果。得知本模型更为符合实 际情况，同时也使得概率积分法能够更为灵活地应 用在倾斜煤层开采沉陷的相关预计与计算当中。 ［ 参考文献］ ［ 1］ 何国清，杨伦，凌赓娣，等 . 矿山开采沉陷学 ［M］． 徐州 中国矿业大学出版社，1991. ［ 2］ 黄乐亭 . 村庄压煤开采的技术途径 ［J］． 中国煤炭，1996 8 33－35. ［ 3］ 王玲，吴侃，奚新丽 . 基于 ArcGis Engine 的开采沉陷预 计系统 ［J］． 测绘科学，2008 4 153－154. ［ 4］ 李培现，谭志祥，齐公玉 . 基于 MATLAB 的开采沉陷预计系 统 ［J］． 中国矿业，2008 11 72－76. ［ 5］ 康建荣 . 任意形多边形工作面多段线开采沉陷预计系统 MSPS［J］． 矿山测量，2000 1 24－26. ［ 6］ 国家煤炭局 . 建筑物、水体、铁路及主要井巷煤柱留设与压 煤开采规程 ［S］． 北京 煤炭工业出版社，2000. ［ 7］ 张华兴，焦传武 . 地表移动精确预计法 线积分运算 ［J］． 矿山测量，1998 4 21－29. ［ 8］ 周万茂，张华兴，何瑞华 . 任意形状工作面拐点移动距求取 方法 ［J］． 煤矿开采，2000 4 13－16. ［ 9］ 张华兴 . 地表移动与变形的等价转换线积分计算 ［J］． 煤矿 开采，1991 2 14－19. ［ 10］ 田锦州，徐乃忠，李凤明 . 误差函数 erf x近似计算及其 在开采沉陷预计中的应用 ［J］． 煤矿开采，2009，14 2 33－35.［责任编辑 李青］ 86 总第 135 期煤矿开采2017 年第 2 期 ChaoXing