基于三维克立格的矿产储量估算系统.pdfVIP

986 矿物学报 2011年 基于三维克立格的矿产储量估算系统 余先川，袁瑜琴，邓维科 (北京师范大学信息科学与技术学院空间信息处理与分析校级重点实验室。北京100875) 本文基于国土资源部危机矿山接替资源找 矿储量约70万t(中国地质科学院矿产资源研究 矿项目“危机矿山三维信息评价系统”之“三维克 所)。 立格储量估计系统的算法和实现”，对三维克立 根据矿床倾角，我们选择主平面方向角O。， 格算法进行了深入的研究，并在此基础上建立了 深度角300为主平面方向。主平面变异性分析 一个应用于地质矿产储量计算的三维克立格系 后，选择方向角0。，深度角-60。为主轴方向，主轴 统(Visual 变程11．8m，块金值和基台值分别为0．07m， C++2005面向对象技术4-OpenGL 0．62 技术)。使系统多个模块具备三维立体显示，具 m，副轴方向角15。，深度角29。，变程6．8m， 备三维交互功能。 块金值和基台值分别为0．08m，0．38m，短轴方 向角900，深度角-25。，变程2．4m，块金值和基台 1基于三维克立格的矿产储量估算 m，0．77 值分别为0．13 ITI，网格化数据采用10× 系统的框架结构 10×10分割，矿物比重设为2。 系统总共分为用户接口及总控界面模块、数 普通克立格的计算得到的储量估算约为 据预处理和统计特征分析模块、变异性分析模 69．6万吨，普通克立格估算方差、交叉验证和品 块、克立格和储量估算模块。所涉及的数据可视 位一吨位图结果显示，系统具有良好的可视化和交 化界面包括“原始数据位置图”、“概率分布图”、 互功能，克立格计算所得的储量估计具有较高的 “直方图”、“变差函数图”、“云状图”、“套合结构 估计精度，且能够给出估计值的估计方差、交叉 图”及“品位一吨位曲线”等视图。 验证和品位吨位图。 数据预处理和统计特征分析模块主要用以 3结论 分析原始空间数据的统计分布特征，以及针对数 Visual 据特点做适当必要的调整和处理，以方便后续的 基于MicrosoftC++2005，开发实现 分析计算。变异性分析或者又称为结构分析，主 三维克立格储量估算系统，并以实际矿产采样数 要工作是计算各个方向实验变差函数，然后拟合 据为例。进行克立格估算，实验结果表明，系统具 一个适合与三维空间的理论变差函数模型，使得 有较高的克立格估值精度。 区域化变量的随机性和结构性都能很好的表现 (1)系统三维克立格运算包含了数据特征统 出来。克立格计算中，软件主要涉及到局部计 计分析与预处理，变异性分析(实验变差函数计 算，根据所得的最终拟合后的变异函数模型求得 算、云状图分析、理论变差函数拟合、结构套合)、 有限区域内区域化变量平均值的无偏最优估计。 克立格计算和储量估算等模块，系统克立格估算 储量计算则主要是根据克立格计算的结果在网