地质的空间三维建模插值技术的研究.docVIP

地质的空间三维建模插值技术的研究 　　【摘　要】露天矿复杂地质体中的各种地质信息，包括地表地形、地下水位、地层断面、节理、风化带分布等，都可以通过野外勘探实测或检测仪器获得，但一般都是散乱数据，而且数量有限。利用这些稀疏、有限、不规则的数据建立三维地质模型难以描述矿山地质体的实际情况，以此模型进行的空间分析难以达到理想的效果。为了利用这些原始数据建立比较理想的地质模型，必须分析原始数据间的联系和规律性，结合地质学原理和相关的地质条件选择合适的插值方法来处理原始数据。 　　【关键词】地质空间；三维建模；插值技术 　　0.引言 　　本文分析多种面向地质建模的空间插值方法后，根据数学地质理论、矿床的赋存条件及原始采样数据的分布特点，分析研究参数优化选择的有效方法，研究关键参数选取的方法，为进一步精确构建地质模型提供了数据保障。 　　1.空间插值的概念与分类 　　空间数据的插值可描述为：给定一组已知空间离散点数据，从这些数据中找到一个函数关系式，使该关系式最好地逼近这些已知的空间数据，并根据该函数关系式推求出区域范围内其他任意点的值。空间插值方法很多，分类标准不一。如按照拟插入点的区间范围可以分为内插和外推；按原始数据空间分布完整性可以分为整体插值和局部插值。整体插值利用所有样本点进行全区特征拟合，而局部插值则仅利用相近的数据点对未知点进行估值，根据内插点的分布范围又可分为分块插值和逐点插值。 　　2.整体插值法 　　整体插值的拟合模型是由研究区域所有采样点的观测值建立的。整体插值主要通过多项式函数来实现的，其特点是不能反映插值区域的局部特性，因此该方法不直接用于空间插值，而是用来检测不同于总趋势的最大偏离部分。从数据中去除一些不符合总体趋势的宏观地物特征后，可用剩余残差来进行局部插值。 　　2.1趋势面法 　　多元回归插值法是一种常用的整体插值方法，一般用于确定数据的大规模的趋势，因此也被称为“趋势面拟合”。其原理是用函数代表的面来拟合现象特征的趋势变化。它的理论假设是地理坐标（X，Y）是独立变量，属性值Z也是独立变量且是正态分布的，同样回归误差也是与位置无关的独立变量。 　　用一个简单的示例来说明，地理或环境调查中特征值Z沿一个断面在X1，X2，…，Xn处采样，若Z值随X值增加而线性增大，则该特征值的长期变化可以用回归方程：Z(x)=b0+b1x+ε进行计算（其中b0，b1为回归系数，ε为独立于 的正态分布残差。）。 　　然而许多情况下，不是以线性函数，而是以更为复杂的方式变化，则需用二次多项式Z(x)=b0+b1x+b2x2+ε或更高次的多项式进行拟合。对于二维的情况，X，Y坐标的多元回归分析得到的曲面多项式，形式如下： 　　一次趋势面的数学模型：Z=b0+b1x+b2y+ε 　　二次趋势面的数学模型：Z=b0+b1x+b2y+b3x2+b4xy+b5y2 　　三次趋势面的数学模型：Z=b0+b1x+b2y+b3x2+b4xy+b5y2+b6x3+b7x2y+b8xy2+b9y3 　　计算是一个标准的多元回归问题，大多数情况下可采用低次多项式进行拟合，一般次数为2或3就可以。次数高的多项式可以很好地逼近观测点，但计算复杂，而且其内插效果使整体分离，降低整体趋势；当多项式次数过高时，非观测点部分的估值会产生大幅震荡。 　　趋势面法反映总体拟合情况，有着总体的趋势，但是并不能很好地反映出局部特征，尤其是在边界上的拟合情况偏差较大。因此，趋势面法主要用于在使用某种局部插值方法之前，模拟大范围内的宏观变化趋势，而不直接用其进行空间插值。 　　2.2加权最小二乘拟合法 　　加权最小二乘拟合法是在趋势面分析基础上引入“距离权”的概念，从而使估值更为合理。具体方法如下。 　　已知样本点(X1,Y1,Z1)，(X2,Y2,Z2)，…，(Xn,Yn,Zn)，估值误差平方和为： 　　Q=(f(xi,yi)-Zi) 　　对每个样本点考虑一个距离权系数，通常为w(di)=，其中di=，ε为一个很小的正数。使总估值偏差和Q最小，则 　　Qmin=(f(xi,yi)-Zi)w(di) 令 　　=0 　　=0 i=1,2,3,...,n 　　通过最小二乘方程建立法方程，求得方程系数矩阵，再将插入点（X,Y）带入趋势面数学模型，可求得待估点高程值Z。 　　加权最小二乘拟合法不仅可以反映总体趋势，还可以反映局部的变化特征，但反映出的局部变化不是十分明显。通常这种方法在矿体连续的情况下（即没有断层的情况），也可以用来加密边界点。如果出现断层情况，在边界处拟合会出现下凹问题。 　　3.分块插值法 　　由于实际的估值区域复杂多变，不可能用一个多项式进行拟合，因此一般不用整体函数法进行内插，而是采用局部分块内插拼接的方法。所谓分块内插是把参考空间分成若干分块，对各分