概率积分参数反演中的开采影响传播角问题研究

1、概率积分参数反演中的开采影响传播角问题研究    摘要：为了根据实测数据反演出可用的开采沉陷概率积分参数，研究概率积分参数的反演，采用VB 语言编制了模矢法求参程序，但是在实际应用中求取开采影响传播角不符合其物理意义。针对此问题，本文提出了6 参数法求取开采沉陷预计参数，并有效解决了概率积分参数反演中的开采影响传播角问题，特别是在数据质量不高时，可以用这种改进的方法获得比较合理的沉陷预计参数。关键词：大地测量技术；概率积分参数；开采影响传播角；模矢法0 引言煤层大面积开采引起覆岩移动，最终导致地表下沉、耕地破坏、岩溶塌陷、潜水出露、沼泽化、山体开裂、滑坡

2、与泥石流、水土流失等一系列问题1。全国各类煤矿开采造成的地表沉陷面积约为60 万hm2，我国每采万吨煤，平均塌陷土地0.2hm2，在村庄稠密的平原矿区，每采出1000 万t 煤需迁移约2000 人。村庄下压煤开采成为矿区面临的主要问题2。要在尽可能多的采出煤炭资源的同时对地表进行保护，就要在开采前对地面损坏进行预计分析。目前开采沉陷的预计方法主要有概率积分法，典型曲线法，剖面函数法，其中我国较多使用概率积分法。要预计结果与实际的地表破坏相吻合，就要求有较为准确的概率积分参数，由于岩层结构的不同各个地区的概率积分参数也不相同3，所以对开采的地表变形进行预计时，要有符合当地地质条件的预计参数，预计

3、参数是由实测数据反演得来的。1 根据实测数据求取概率积分参数的方法概率积分法是因其所用的移动和变形预计公式中含有概率积分（或其导数）而得名。由于这种方法的基础是随机介质理论，所以又叫随机介质理论法。应用概率积分法的相关参数可以进行开采沉陷预计，这一模型在矿山开采沉陷中应用最广，处理数据时用该模型作为基础1。概率积分法关于下沉盆地的共包含8 个参数，即下沉系数q，主要影响角正切tan，水平移动系数b，左、右、上、下拐点偏移距s1、s2、s3、s4，开采影响传播角。前7 个参数是关于下沉盆地的。概率积分的基本原理如下：在倾斜煤层中开采某单元i，单元开采引起地表任意点(x，y)的下沉( 最终值 )为

4、：Weoi (x,y)=(1/r2) ×exp(-(x-xi)2/r2) ×exp(-( y-yi+li)2/r2) ；li=Hi/tan；式中(xi,yi)为 i 单元中心点的平面坐标；(x,y)表示地表任意一点的坐标。设工作面范围为：0p，0a 组成的矩形。则地表任一点的下沉为：W(x,y)=W0Weoi(x,y)dxdy式中：W0 为该地质采矿条件下的最大下沉值（单位mm），W0mqcos；p 为工作面走向长（单位m）；a 为工作面沿倾斜方向的水平距离（单位m）。此公式也可以写为：W(x，y)W°(x)W°(y) /W0式中 W°(x)为

5、倾向方向达到充分采动时走向主断面上横坐标为x 的点的下沉值，W°(y)为走向方向达到充分采动时倾向主断面上横坐标为y 的点的下沉值。公式如下：倾向充分采动，走向半无限开采时走向主断面下沉值：W(x)=W0/2×（erf(sqrt×x/r)+1倾向充分采动，走向有限开采时走向主断面下沉值：W°(x) =W(x)-W(x-l)走向充分采动，倾向有限开采时走向主断面下沉值：W°(y) =W(y)-W(y-L)l=D3-s3-s4；L=(D1-s1-s2)×sin(+)/sin式中 D1 为工作面倾向长；D3 为工作面走向长；l 为走向计算边

6、界；L 为倾向计算边界；为开采影响传播角（近似于最大下沉角，用最大下沉角来计算）；为煤层倾角。2 求参中的开采影响传播角问题本人采用 VB 编制了概率积分求参的程序，在求参的过程中发现求得的参数不符合参数的物理意义，尤其是开采影响传播角的结果与实际不符。针对这个问题，本文认为下沉盆地曲面只需要6 个参数。开采影响传播角是倾向主断面预计的特有参数，一般认为与煤层倾角有关，即=90-k，式中k 为小于1 的常数，一般取值在0.5-0.8 之间2。开采影响传播角和最大下沉角类似，最大下沉角是在倾斜主断面上，由采空区的中点和地表移动盆地的最大下沉点（在地表水平上的投影点）的连线与水平线之间在煤层下山方

7、向一侧的夹角4。分析公式，不难发现:倾向充分采动，走向有限开采时走向主断面下沉值：W°(x) =W(x)-W(x-l) 5走向充分采动，倾向有限开采时走向主断面下沉值：W°(y) =W(y)-W(y-L) 5l=D3-s3-s4；L=(D1-s1-s2)×sin(+)/sin其中s2 和这两个参数用于确定L，并且只用来确定L，所以认为可以用一个参数来代替。本文用s2 代替。而在求参模型中，右拐点偏移距s2 和最大下沉角都只参与了倾向计算边界L 的运算，所以只能用其中一个值参与模型拟合，即先将最大下沉角略去，采用L=(D1-s1-s2)。再根据拟合出的概率积分模型直

8、接求取开采影响传播角，由于开采影响传播角和最大下沉角类似，所以可以使用最大下沉角的定义来求取6。随后对拐点偏移距s2 进行修正。即(D1-s1-s2)= (D1-s1-s2')×sin(+)/sin修正后的右拐点偏移距s2'=D1-s1-(D1-s1-s2)×sin/ sin(+)此时的拐点偏移距是拐点沿竖直方向投影到煤层所取得的值，而拐点偏移距的意义是拐点投影到煤层上所取得的值，也就是沿煤层法向投影，因此，倾向上的左右拐点偏移距均应进行修正S1=s1+H1/tan；S2=s2'-H2/tanS1,S2 即为最终求得的拐点偏移距。3 实验分析该实验对

9、某工作面进行了动态参数求取，使用三维激光扫描数据进行求参。该工作面煤层厚度5.84m，煤层倾角=3.4°，工作面推进速度约为6m/d。分别采用两种方法求参，一是7 个参数同时参与计算，二是先计算6 个参数，然后对反求开采影响传播角，并修正左右拐点偏移距。由实验结果不难看出，两次的拟合中误差相差不大，而且都很小，说明都达到了最优。但是第一组求参结果中的下沉系数与开采影响传播角明显与实际不符，数值已经脱离了其物理意义，这说明其结果不符合实际。一旦所求的参数失去了物理意义，只具有拟合这个工作面下沉曲面的数学意义，那么这一套参数就只适用于这一个工作面，不能够对于相似地质条件的工作面同样适用。

10、由图1 和图2 对比可以看出，两组参数拟合的下沉曲面相差不大，也就证明了该改良的求参方法同样适用于该工作面。所以在数据质量不高的情况下，可以采用该方法获取比较合理的预计参数。4 结论（1）下沉值的轻微变动及其他参数的变化对开采影响传播角影响很大，确定一个下沉盆地用6 个参数控制可获得较好的收敛结果，在数据质量较高的情况下两种方法都可以。（2）下沉系数q，主要影响角正切tan，左右拐点偏移距，共4 个参数即可确定倾向方向的下沉曲线。（3）概率积分参数不是纯数学参数，它有一定的物理意义，所以要符合实际的情况。（4）在数据质量不高的时候，可以用这种改进的方法获得比较合理的沉陷预计参数。本论文网专业从事会计硕士论文业务！如需转载请保留一个链接：代写硕士论文参考文献 (References)1 何国清,杨伦,凌赓娣.矿山开采沉陷学M.江苏徐州：中国矿业大学出版社,19942 连达军,汪云甲.”三下”开采综合评价体系研究J中国矿业大学学报,2005 ,34(1):98