克立格法在SURPAC软件中的实现及应用.pdf

克立格法在 SURPAC 软件中的实现及应用 冯超东 胡乃联 杨鹏 北京科技大学土木与环境工程学院 北京 100083 摘要： 克里格法是地质统计学中最主要和最基本的一种局部估计方法. 本文通过某金矿山实例，介绍了在 SURPAC软件中如何应用克立格法进行品位估计，储量计算及结果校验过程。它利用了脉矿的空间变化及其 品位的变化特征来计算样品的加权因素和品位，改变了传统复杂的求体积、体重的储量计算方法，列出了 在SURPAC软件中用克里格法计算矿区储量过程。 关键词：克立格法；地质统计学；SURAPC软件；品位估计 中图分类号： 0 引言 1 克立格法简介 克里格法储量计算是南非采矿工程师 克里格法是地质统计学中最主要和最 克里格于1951 年首次提出来的。由于它具 基本的一种局部估计的方法。所谓局部估计 有传统地质学方法所无可比拟的优越性，西 是指求局部区域内区域化变量平均值的估 方发达国家在地质学和矿业界一直大力开 计量。克里格法的基本原理为：克里格法是 展地质统计学的研究和应用，并广泛运用于 把矿体划分成许多小块段(待估块段)，在充 地质研究和储量计算过程中。目前，在国外， 分考虑信息样品的形状、大小及其与待估块 地质统计学已形成了一套完整的理论和方 段相互间的空间分布位置等几何特征以及 法体系。先后涌现出一批较为成熟的地质统 品位的空间结构之后，为了达到线性、无偏 计学储量计算应用系统软件，如美国Mintec 和最小估计方差的估计，而对每一信息样品 公司的MineSight软件，加拿大GEMCOM 公司 值分别赋予一定的权系数，最后进行加权平 的GEMS软件，英国MICL公司的DATAMINE软 均来估计块段品位的方法。这些权系数则是 件，澳大利亚SURPAC公司的SURPAC软件等。 通过建立、求解克里格方程组计算而得。它 加拿大GEMCOM公司并购澳大利亚SURPAC公 以最小的估计方差给出对于块段平均品位 司后，成为全球最大的矿业软件公司，占有 的无偏线性估计量。 40％以上的矿业软件市场。 克里格法包括线性克里格法和非线性 克里格法在我国的应用和研究起步相 克里格法。线性克里格法的求解问题实际上 对较晚,但发展相当迅速。地质统计学是在 就是确定线性组合的权系数，这种权系数是 1977年由美国福禄尔采矿金属有限公司 通过求解克里格方程组而获得的。包括简单 （Flour Mining Meta Incorporation） 克里格、普通克里格、泛克里格、协同克里 H.M.Parker博士随美中贸易全国委员会矿 格和正态克里格等。非线性克里格法包括对 业代表团来华访问，传入我国，国内一些地 数正态克里格、析取克里格、指示克里格、 质、冶金单位及部分大专院校等部门逐渐展 概率克里格等。正态克里格法和对数正态克 开了对地质统计学的研究应用工作。1989年 里格法是连接两者的纽带，其过渡关系为： 11月召开的全国第一届地质统计学学术讨 一般的线性克里格法—正态克里格法—对 论会，地质统计学在我国的发展进入了一个 数正态克里格法—一般的非线性克里格法。 新的阶段，理论研究更加深入，涉及的方法 克里格储量计算方法是以矿石品位和矿床 原理更加广泛。