开采沉陷数据可视化研究.docVIP

开采沉陷数据可视化研究 　　摘 要: 利用C# 语言实现了煤矿开采地表移动与变形等值线(下沉、倾斜、曲率、水平移动 、水平变形)的绘制,在此基础上进行了影响区域任意两点间的剖面线的绘制,基本实现了 开采沉陷数据的可视化,可对采区上方建筑物进行了破坏等级划分,从而判别建筑物是否受 到地下开采的影响以及受到开采影响的程度。 　　关键词:开采沉陷;等值线;剖面线 　　中图分类号:TD327 文献标识码:A 文章编号:1672-1098(2009)03-0011-05 　　 　　我国是一个矿产资源丰富的国家,其中煤炭资源在一次能源的生产或消费结构中均占到 75%。当地下煤层被开采出来以后,开采区域周围岩体的原始应力平衡状态受到破坏,应力 重新分布,达到新的平衡。在此过程中,开采煤层的上覆岩层将产生移动、变形与破坏,当 开采面积达到一定范围后,移动与变形将波及到地表,使地表产生沉陷[1],随之 ,沉陷区 范围内的农田、房屋及其他建筑物遭到破坏。因此,开采沉陷数据的可视化可以实现对地表 移动、地表破坏的更直观和更深入的理解[2]。 　　 　　1 等值线的生成 　　 　　等值线是一种形和数的统一,由于等值线图看起来非常直观、形象,因此在水利、土木 、地质等工程和技术领域内得到广泛的应用。 　　目前常用的等值线生成算法主要有不规则三角网格法和规则矩形网格法,由于煤矿开采 沉陷数据采集的规律性,故采用了规则矩形网格法。在实际应用中的数据格式规定方面,拟 采纳如下格式:?? 　　 　　 　　 　　其中,第一行各个字母分别表示行数、列数、x方向上的最小值、x方向上的最大值、y 方向上的最小值、y方向上的最大值、x方向上的间距、y方向上的间距。从第二行开始至最 后一行,每一行上的各个字母分别表示p点的x值、y值、下沉值、倾向倾斜、走向倾斜、倾 向曲率、走向曲率、倾向水平移动、走向水平移动、倾向水平变形、走向水平变形。 　　1.1 等值点的确定 　　绘制等值线,首先要逐一计算每一个网格单元棱边与等值线的交点,然后把这些交点连 接起来。而要判断网格边上是否存在特定的等值点,可以使用如下公式 　　(Avalue-value)×(Bvalue-value) 　　其中Avalue、Bvalue分别是相邻网格点(A点和B点)的属性值,value是所指定的 属 性值。 若上式小于零, 则说明在A点和B点之间存在所需要的等值点。 若上式大于零, 则 说明 在A点和B点之间不存在所需要的等值点。 若上式等于零, 则应该先加上一个微小量shift (一般为区域范围的万分之一),然后再进行判断(见图1)。 　　 　　在图1中,m为行数,n为列数,ContourIJ[i,j]为第i行第j列的属性值,MyP oint是一个自定义的类,其声明如下: 　　class MyPoint 　　{ 　　//存放坐标点 　　public double x; 　　public double y; 　　…… 　　//下沉(毫米) 　　public double w; 　　//倾向倾斜、走向倾斜、倾向曲率、走向曲率、倾向水平移动、走向水平移动、倾向水平 变形、走向水平变形 　　public double ix, iy, kx, ky, ux, uy, ex, ey; 　　} 　　等值点的坐标可通过下列算式得到 　　y=y min+i×dy 　　x=x min+ 　　 　　由上述过程即可得到,等值线与格网横边的交点。类似地,也可很方便地获得等值线 与格网纵边的交点。 　　1.2 等值线的追踪 　　在计算出全部等值点以后,必须有规则、有次序的将它们逐点连成等值线,这就是等值 点的追踪问题。此外,同一属性值下可能有若干条等值线分支,需要加以区分,并完整的绘 制出来。 　　1.2.1 等值点的追踪 对格网中的任何一个单元格而言,若此单元格与等值线相交,则交点数必定为偶数个, 因为对一个封闭的单元格来说,相交的等值线一旦有一个进入点,则必有一个出去点。由此 ,对等值点的追踪,可以先取一个单元格,根据进入点,遍历它的四条边,求出出去点,然 后再根据出去点得到下一个需要遍历的单元格。 　　在求出去点的时候,可能会遇到出去点并不唯一的情况。此时,就需要确定等值线在单 元格中的走向,一般按照最短距离的原则进行判定。值得注意的是,在确定等值点的时候就 已经排除了等值点处于格网点的情况,故此时就不用再作考虑了。 　　1.2.2 等值线的搜索 搜索等值线的关键是如何找到线头,再从线头搜索到线尾。为搜索的便利可设置一个临时点 ,用以记录第一个访问的点。 　　由于定义域是有界的,等值线即可能是闭合的,也可能是不封闭的。所以需要分别从开 曲