游明亮;伍岳;朱朋
【摘 要】Digital mine is an informational and intelligent product which appears in the 21st century un-der the new technology.Building mine 3D geological modeling as the foundation and core is to realize dig-ital mine.In order to achieve the complex mining character of 3D geological modeling-economics,authen-ticity and effectiveness, how to pre-process the discrete DEM spatial data is the difficult point and pre-condition for building realistic mine 3D geological modeling, the article describe that according to the dif-ferent mine and geomorphic environment conditions,several basic spatial interpolation methods are used to interpolate the discrete DEM spatial data and then the interpolated data is conducted by the quality evalu-ation and cross-validation to determine which the optimal spatial interpolation methods is beast.Analysis results indicate that different l spatial interpolation methods adapt to different geological mining fields, and the DEM spatial data of discrete optimization of spatial interpolation become most effective and central method to accomplish the project
of mining 3D geological modeling realistically.%数字矿山是21世纪新技术条件下信息化与智能化的产物，建立矿山三维地质体建模是实现数字矿山的基础与核心。为了实现复杂矿山三维地质建模的经济性、真实性和有效性，如何对离散的DEM空间数据进行前期处理，是建立逼真的矿山三维地质模型难点与前提。针对不同矿山与地貌环境条件下，采用几种基本空间插值法对离散的DEM空间数据进行插值，然后对其插值后的数据进行质量评定与交叉验证，确定最优化的空间插值方法。分析结果表明不同的空间插值方法适应不同的地矿环境领域，对离散的DEM空间数据进行最优化的空间插值，是实现逼真的矿山三维地质体建模这一目最有效最核心的方法。
【期刊名称】《湖北民族学院学报（自然科学版）》
【年(卷),期】2014(000)004
【总页数】3页(P474-476)
【关键词】空间插值法;三维地质体建模;交叉验证;质量评价
【作 者】游明亮;伍岳;朱朋
【作者单位】三峡大学土木与建筑学院，湖北宜昌443002;三峡大学土木与建筑学院，湖北宜昌443002;三峡大学土木与建筑学院，湖北宜昌443002
【正文语种】中 文
【中图分类】P208
数字矿山(digtal mine，简称DM)是运用计算机技术、信息技术与可视化技术等高新技术对整体矿山的数字化［1－2］．然而矿山的三维地质建模是实现数字矿山的核心技术，也是实现数字矿山的基础．三维地质建模(3D Geoscience Modeling)最早是有加拿大学者Simon W．Houlding 于1993 年提出［3］．所谓三维地质建模，就是运用计算机技术与可视化技术，在三维环境下对地质体进行模拟与表达．然而复杂的三维地质建模依赖于原始数据，矿山的三维地质体建模的数据主要来自于钻孔数据和勘测的DEM 数据，但是这些数据是有限的、无规则的且离散的，然而矿山地质空间数据往往又是复杂多变且不连续的数据［3－4］．稀疏随机且不足的数据对建立矿山三维地质模型是十分困难的，花费高昂的代价来获得充足采样数据是不现实的，所以对原始采样数据进行插值是解决三维地质建模的关键技术之一［5］．因此，为了最大限度的获得最佳的矿山三维地质模型，准确的模拟复杂的地
质现象与地质体内的地质构造，并灵活表达三维地质实体，关键在于如何选取合适的插值方法．
本文通过研究几种常用的不同插值方法:克里金插值法(Kriging)［4］，反距离加权插值法(inverse distance to a power)［6－7］，最小曲率插值法(minimum curvature)［8］，线性插值三角网法(triangulation with linear interpolation)［9］，在矿山三维地质建模的应用，并通过交叉验证不同插值方法的插值精度与质量，合理地选择最佳的插值模型［5］，如图1．
1 插值方法质量评价及应用实例
图1 插值方法的过程Fig．1 Interpolation method process
插值方法质量评价是在地质建模时评价某一插值方法的质量，基本过程是通过不同插值方法对原始数据进行插值，然后对插值后的数据随机选取相同的区域通过不同的指标进行质量评价，以求得最精确的插值方法．评价不同插值方法的质量指标有均方根预测误差(RMSPE，root－mean－square prediction Error)、计算残差(residuals)、相对均方差(RM
SE，relative mean－square error)平均估计误差百分比(PAEE，percent average estimation error)［6］．
本文选取湖北某甲，乙两个不同矿山，甲矿山处于荆当盆地西北部边缘，为鄂西山区与江汉平原接壤地带，属于丘陵低山地貌，地形多具缓坡，地形坡度角5～35°．乙矿山处于荆当向斜远安台褶束西部通城河大断裂西南部，矿区地貌类型属构造侵蚀丘陵地貌，地形多具陡坡、冲沟特征，地形坡度角20～55°．通过几种基本常用的空间差值方法对甲、乙矿山地质建模过程中对不同区域的离散数据进行插值，然后选取相同的某区域插值后的数据进行质量评价与交叉验证，进而确定最优的插值方法，以更加逼真地和准确地建立相应的地质体模型．
图2 甲矿部分区域DEM 模型(线性插值三角网法建立)Fig．2 The first mine some area of the DEM model(established by triangulation with linear interpolation method)
图3 乙矿部分区域DEM 模型(克里金法建立)Fig．3 The second mine some area of the DEM model(established by Kriging method)
图4 甲矿部分区域DEM 线框模型(线性插值三角网法建立)Fig．4 The first mine some area of the DEM wireframe model(established by Triangulation with linear interpolation method)
图5 乙矿部分区域DEM 线框模型(克里金法建立)Fig．5 The second mine some area of the DEM wireframe model(established by kriging method)
三峡大学是几本通过质量评价及交叉验证的方差和绝对平均偏离说明插值后的数据离散性和总体偏差，均方根预测误差和平均估计误差百分比预测插值后的数据预插值与真实值得接近程度，均方根预测误差值越小，插值数值就越接近他们的真实值，同时平均估计误差百分比趋于零，则表明待插值的估计值是无偏的．综合上表交叉验证结果表明，甲、乙矿山空间插值方法的效性比，对甲矿山线性三角网法是最优的，其次是径向基函数法与克里金法，最差的是最小值曲率法．乙矿山克里金法是最优的，其次是线性三角网法和反距离加权法，最差的是径向基函数法．
虽然甲，乙矿山都属于丘陵地形，但是由于甲矿山的地形特征相对比较简单，且地形坡度比较缓，而乙矿山的地形特征比较复杂，地形坡度比较陡且多冲沟．结合对甲、乙矿山插
值法的交叉验证表明，以上几种不同的基本空间插值法都有一定的局限性，对于复杂的地形特征各种插值误差明显增大，即插值效果明显下降．通过误差分析理论说明不同的插值法适应不同的地形特征，但是对于不同复杂程度的地形特征克里金空间插值法均能实现比较有效的空间插值，且插值的精度比较高以及插值效果比较稳定，即说明在实际工程应用中克里金法是一种比较实用的空间插值方法．
表1 几种基本插值法在某甲矿局部区域插值结果的质量评价及交叉验证Tab．1 Several basic interpolation method in a mine local area interpolation results of quality evaluation and cross－validation插值法 均方根预测误差RMSPE/m 平均估计误差百分比PAEE/% 方差/m2 绝对平均偏离/m反距离加权 0．607 0．199 0．368 0．277克里金 0．612 0．124 0．037 0．288最小曲率 1．396 0．634 1．949 0．732线性三角网 0．179 0．011 0．032 0．055径向基函数 0．244 0．020 0．059 0．104局部多项式 0．752 0．170 0．566 0．440
发布评论