实现矿山开采沉陷预计的Visual Basic编程思想探讨.pdf

《实现矿山开采沉陷预计的Visual Basic编程思想探讨.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第2期矿山测量No．220l4年4月MINESURVEYINGApr．2014doi：10．3969／j．issn．1001—358X．2014．02．28实现矿山开采沉陷预计的VisualBasic编程思想探讨赵阳，李建伟，谢瑞斌(枣庄矿业集团高庄煤业有限公司，山东枣庄277605)摘要：文中主要探讨了在VB环境中，用概率积分法和特殊函数编程的方法，实现矿山开采沉陷预计的技术和思想，并结合部分模块代码和工程实例进行了说明。关键词：概率积分；沉陷；预计；编程中图分类号：TD327文献标识码-B文章编号：1001—358X(2014)02—0087—03在采

2、矿活动中，往往利用概率积分法进行地表PublicFunctionnerrf(X)下沉及变形预计，本文主要探讨在VB环境下，如何If(X>=01)Then编制程序，进而实现用概率积分法自动进行开采沉nerrf=ngam2(0．5，X$X)陷预计，下面就矩形工作面和不规则工作面两种情Else况分别进行阐述。nerrf=一ngam2(0．5，X：IcX)EndIf1主要原理公式EndFunction根据《开采沉陷学》中地表移动盆地任意点的下子程序中的ngam2()函数属不完全伽马函数，沉计算，在三维条件下，某采区开采后引起的地表点可以直接在

3、程序集》中查找。A的下沉值为：2矩形工作面沉陷预计=×÷}[erf(÷)川一2．1坐标转换[erf(~-)川}×{[e矿(÷)在应用概率积分法进行沉陷预计时，需要建立如图1所示的坐标系，称之为采区坐标系。+1]．[。)+1]}=W⋯CC式中，C、C分别为A(，Y)点沿走向和倾向的下沉分布系数，它可以利用沿走向和倾向两个半无限开采的叠加而求得。c=÷{[e÷)+1]一[erf()+11}c：{[)+1]一[erf()+1]}图1矩形采区计算坐标系示惹图Zrr在进行沉陷预计前，必须将实际的开采边界点上式中，erf()称为误差积分(概率积分)，由矿区坐标系转换为

4、以O点为原点的采区坐标系。在vb环境中，其转换过程如下：其数学表达式为』edA。q读取开采边界点一存入动态数组pt()一对动态对于概率积分部分的编程，可以参考《QBASIC数组进行遍历找出X为最小的点一记录点号，_『，一+常用算法程序集》中的特殊函数编程方法，误差函数计算点至点i的方位一计算0点坐标一将所有要预计的点转换为采区坐标系。erf()=edt函数子程序如下：求取0点坐标，可采用工程测量中支导线坐标87第2期矿山测量2014年4月计算方法，代码如下：Y1=Sqr(3．14159261Y／rDimazimuAsDouble定义坐标方位角Y2=Sqr(

5、3．1415926)(Y—width)／razimu=azimuth(pt(j)，pt(i))+1．571求取isl1=nerf(Y1)至0点的坐标方位角s22：nerrf(Y2)IfaT．imu>23．1415926Thenazimu=azimuaybir=(sl1+width)／2一2}3．14l5926aylbir=(s22+width)／2pt(i)．X=pt(i)．X+hd：IcCos(1tzimu)hd代cyaybir—aylbir表点i至点0的距离resu=wmax术cx’Iccypt(i)．Y=pt(i)．Y+hd}Sin(azimu)En

6、dFunction其中，azimuth()为求取两点间坐标方位角的函上述求取任意点下沉值的函数用到的参数有：数，可自己根据测量学原理在VB环境中定义函数实采区长度length，采区宽度width，最大下沉值wmax现。，影响半径r，及任意点坐标、l，。坐标转换核心代码如下：3不规则工作面沉陷预计DimroteangleAsDouble定义旋转角roteangle=3．1415926／2一azimuth(or(i)，pt对于任意形状的工作面开采沉陷预计，可以直(j))接在矿区坐标系中进行，沿或Y方向将其划分为Ifroteangle<0Thenroteangl

7、e=roteangle+2许多小矩形工作面，然后由概率积分法预计出每个3．14l5926小矩形工作面对任意点引起的下沉和变形值，最后ReDimPreserved(k)‘定义点类型动态数组叠加各小矩形工作面开采的沉陷预计值即可以得到d(k)．X=pt(i)．X+XCos(roteangle)+任意形状工作面的地表下沉变形值。YSin(roteangle)d(k)．Y：pt(i)．Y+YCos(roteangle)一yXSin(roteangle)其中，x、Y坐标为覆盖下沉变形区域的任意点坐标。2．2沉陷预计任意点下沉值的预计，可以采用循环结构语句，对覆盖变形

8、区域的每个点都调用求取下沉值的函数，直到所有点计算完毕。求取下沉值