一种基于块体模型和克里格法的矿体储量计算.pdfVIP

SerialNo．508 现 代 矿 业 总第508期 August．2011 M0DERN MINING 2011#-8月第8期 · 技术交流 · 一 种基于块体模型和克里格法的矿体储量计算 王 斌 李发本 (1．洛阳栾川钼业集团有限责任公司；2．西安建筑科技大学) 摘 要 在三维块体模型的基础上，结合普通克里格法研究了储量计算方法，实现了实验变异 函数求取，变异函数球状模型的拟合、结构分析和椭球体的确定，并将三维图形引擎库Hoops与软 件开发方法相结合，设计 出三维可视化储量估算系统。将其应用于国内某大型钼矿 中，实现了该矿 体的储量计算。 关键词 克里格方法 储量计算 变异函数 矿产储量是矿业项 目评价、矿山企业生产中的 系数，最后用加权平均法来对待估块体的未知量进 主要指标之一。受矿山技术水平和矿体复杂性的影 行估值 。 响，传统的手工地质块段法估算矿产储量工作量极 当区域化变量满足内蕴假设并且采用线性估计 大，而且误差较大。本文在三维块体模型的基础上， 量，属于线性平稳地质统计学范畴使用克里格法。 基于空间地质统计学的克里格估值法进行矿块品位 设Z(u)是二阶平稳的，任一待估块体 ，其真 估计，进而对矿产储量进行统计估算，根据矿产地质 值z 的估计值是估计临域内n个信息值的线性组 或矿床勘探所获的矿床 (或矿体)资料、数据，确定 合，其估计值一般表达式为 矿产储量、品位以及空间分布。常用的计算方法有 n Z (“)= ∑． AZ(Z／)． (1) 断面法、开采块段法、算术平均法、地质块段法、等值 线法、三角形法等 1J。20世纪80年代以来，克里格 求出各个加权系数 A (O／=1，2，…，n)，使得 估计技术的理论和应用得到了前所未有的发展 J。 z (z／)是z 的无偏估计量且估计方差最小。 国际矿业界已经把空间地质统计学作为矿山地质储 1．1 无偏条件 量计算的标准方法，三维地质学可视化系统 已广泛 无偏性就是参数估计z ()与真值Z 之间偏 应用于各个工业领域的地质建模 中，如英 国的Dat— 差的数学期望为零，即E{z}=E {z (n)}成立。 amine、澳大利亚Maptek公司Vulcan、ESRI公司的 在二阶平稳条件下，E{Z}=E{z ()}=m，而 ArcGIS等都包括三维矿体可视化模型和简单储量 17, n n E{z (u)}=E{∑AZ }=∑AE{Z }= A m ， 计算l3 。我国矿产资源储量计算方法经历了三个 ． ． ， 阶段的发展历程：传统储量计算方法；空间地质统计 若E{z}=E{z (M)}，则无偏条件就是： 学储量计算方法；多种矿产资源储量计算方法并存 ∑A =1． (2) 与共同发展 。 1．2 最佳性条件 l 克里格估值方法 最佳性是指估计值与真值之间偏差的平方的数 克里格估值法是在考虑了信息样品的形状、大