地矿工程中空间插值方法的应用.docxVIP

地矿工程中空间插值方法的应用 1 规则网络体系 空间数据分析方法在地矿工程的计算机可视化技术中起着重要作用。众所周知,由于地质勘测工程及地质体的大尺度特征,使得地矿工程的勘测数据一般只能是对三维连续数据场的有限采样。如此取得的原始数据,具有明显的离散性、稀疏性。而为了实现地矿数据的计算机体视化及地矿工程的真三维显示。需要将这样的原始数据处理成规则的体数据,即通过数据的空间插值,形成由均匀网络或规则网络组成的结构化数据,其中的每一个网格是结构化数据的一个元素,称为体元或体素(Voxel),而数据场的函数值则均匀地分布在体元的8个网格顶点上,如图1所示。 可见,数据的空间插值是地矿工程计算机体视化过程中必不可少的一项工作,插值算法的精度和效率将直接影响地矿工程的三维显示效果。 2 构造等价面线的插值方法 空间插值一般可以定义为由已知的空间数据来估计或预测未知的空间数据值。空间数据插值方法有多种,但对于各种插值方法的研究多着眼于数学的角度,其主要目的在于建立各种空间特征曲面(线)。 本文根据地矿工程体视化系统研究与开发的需要,对比分析有关地矿体属性的主要空间插值方法,主要目的是:①补值:地矿工程中经常会由于各种原因造成某点采样(或观测)值的丢失,要利用适当的插值方法对该点进行补充插值;②构造等值面(线):即利用插值方法在二维空间或三维空间构造等值面(线)来表现矿体的品位等属性数据的空间分布特征,以建立所研究对象的表面模型;③构造空间规则体数据:即把离散、稀疏分布的地矿体空间数据通过插值形成规则的体数据,以建立所研究对象的体模型、实现空间数据的体视化。 从这三个目的来看,一般的插值方法多用于第二种目的——构造等值面(线),通过对有限的采样(观测)数据的空间插值,拟合出光滑的曲线或曲面,以表现数据的空间变化,例如,趋势面插值、多元回归插值、傅立叶级数插值、径向基函数插值、有限元插值以及B样条插值方法等,几乎均着眼于构造等值面(线)这一目的,建立研究对象的表面模型。 另一方面,用于地矿体补值的插值与用于构造地矿体空间规则体数据的插值的实质基本相同,其目的均属于获得空间点的信息。两者的差异仅在于补值是要得到地矿体的某一指定点上的信息,属于局部范畴;而构造空间规则体数据则是要得到遍布整个地矿体的一系列点的信息,以反映对象空间属性的整体变化,属于全局范畴。由此可见,用于这两种目的的插值与用于构造等值面(线)的插值的主要区别在于:前者是为了获取对象即地矿体的空间点的信息,以利用这些信息来表现对象的空间分布特征,进而由空间离散分布的点数据形成对象的规则体数据,以揭示对象的内部属性及结构;而后者是间离散分布的点数据形成对象的规则体数据,以揭示对象的内部属性及结构;而后者是为了得到用于反映对象的线或面的信息,利用表面或线来表示对象的空间形态或属性的变化,其插值函数往往需要满足数学连续性要求,以保证空间线、面的数字光滑。 3 规则体数据插值方法的特征 为进行地矿工程体视化系统及技术的研究与开发,我们进一步讨论用来形成规则体数据的插值方法的特征。可以用来形成地矿体空间体数据的插值方法主要有Shepard法、四面体剖分法和克立格法。 3.1 关于权重的关系 Shepard法,通常称为距离幂次反比加权法,是地矿工程领域中经常采用的插值方法之一。该方法将待求点的函数F(x,y,z)定义为各已知点函数值fk对待求点与各已知点之间距离幂次反比的加权平均值,即 F(x,y,z)=∑k=1nfk?Wk(x,y,z)(1)F(x,y,z)=∑k=1nfk?Wk(x,y,z)(1) Wk(x,y,z)=1/dk(x,y,z)μ∑k=1n/dk(x,y,z)μ(2)Wk(x,y,z)=1/dk(x,y,z)μ∑k=1n/dk(x,y,z)μ(2) 在权重Wk(x,y,z)的式2中,dk(x,y,z)=(x?xk)2+(y?yk)2+(z?zk)2??????????????????????????√dk(x,y,z)=(x-xk)2+(y-yk)2+(z-zk)2,表示由点(x,y,z)到点(xk,yk,zk)的距离,μ为幂次,一般取2。 根据这一定义,可见此方法具有如下基本特征:①易于理解,便于计算。②考虑了待求点与已知点之间的距离对插值结果的影响,即距待求点近的已知点数据对插值结果的贡献大,距离远的对结果贡献小。但当距离远到一定程度时,基本上不再对结果有影响,这样的已知点便成为“局外点”。③没有考虑已知点属性之间的相关性对权函数Wk(x,y,z)的影响。例如,假设已知点如图2中的a1,a2,a3和a4所示,它们到待求点a的距离相等,因此,由Shepard方法所算得的各已知点的权值相等。但a1与a2间的距离小于a3与a4间的距离,所以,相距较近的a1与a2两点的属性(例如