地质统计学原到理与地质建模方法.pptVIP

* * * * * * * * * * * * * * * * * 三、地质统计学基本理论 什么是插值？ 用已知测量值对未知位置计算一个值。例如图中位置u周围的六个测量值对位置u的值进行计算。 克立格插值（Kriging interpolation)是根据变异函数模型而发展起来的一系列地统计的空间插值方法，包括： 简单立格法（Simple Kriging）; 普通克立格法（Ordinary Kriging）; 指示克立格法（Indicator Kriging）; 析取克立格法（Disjunctive Kriging）; 协同克立格法（Cokriging）等。 下面仅对普通克立格法作一些简单介绍。 三、地质统计学基本理论 在本征假设下： 估计值： 普通克立格方程组： 三、地质统计学基本理论 四、地质品位模型计算过程 孔口 距离 方位角 倾角 轨迹 建模一般从钻孔数据开始 四、地质品位模型计算过程 地质数据库建立 建立孔口表，测斜表，样品表，岩性表 钻孔数据错误检验，包括钻孔信息缺失，重复记录，样品段缺失、重叠等 建立钻孔数据库 原始样基本统计 确定分布形式 确定分布参数 判定并处理特异值 负值看作无样品 特高值根据直方图尾部状态确定，大小以占分布区域97.5%时的值选取，也可取均值的6~8倍 特高值直接剔出掉，较多时以均值替代 地质数据库建立 四、地质品位模型计算过程 组合样长：最好为平均样长 组合方法： 沿钻孔 沿高程（台阶） 根据岩性 最小样长：取75% 样品组合 地质数据库建立 原始样基本统计 四、地质品位模型计算过程 分布形式、统计参数应与原始样相同和接近 组合样统计分析 地质数据库建立 原始样基本统计 样品组合 四、地质品位模型计算过程 参数取值 容差角：15~30度，据实际情况 容差限：取滞后距（可不设） 滞后距、滞后次数 根据样本空间大小，试算，若曲线离散，可调大容差角、容差限、滞后距； 根据试算曲线趋于平缓时的总距离以及滞后距，重设滞后次数； 变异函数计算 地质数据库建立 原始样基本统计 样品组合 组合样统计分析 四、地质品位模型计算过程 球状模型 微调块金、基台、变程 理论曲线尽量与前几个点趋势逼近 不同方向块金、基台理论上应相等，实践中常不同，一是按相同设置不考虑方向计算并选定参数，一是取各方向参数的平均值 变异参数确定 地质数据库建立 原始样基本统计 样品组合 组合样统计分析 变异函数计算 四、地质品位模型计算过程 以确定的变异函数，根据已知样品进行试估算，并计算估值与真实值的误差，克里格方差等； 误差应属正态分布； 误差均值应为0； 误差方差接近克里格方差； 误差均值为零的正态分布，且峰度较高 交叉验证 地质数据库建立 原始样基本统计 样品组合 组合样统计分析 变异函数计算 变异参数确定 四、地质品位模型计算过程 块段模型用于模拟非均匀矿体的空间分布，由大小不等且非重叠的长方体块构成。每个块包含该块所在位置矿体的各种元素品位信息和其它一些属性。 例如，块段模型中的一个块可以包含铜、铁、金的品位值，这些属性则可以用于各种统计分析。 块段模型是使用传统或者地质统计学储量计算方法的前提。 块段模型建立步骤如下： 地质数据库建立 原始样基本统计 样品组合 组合样统计分析 变异参数确定 变异函数计算 交叉验证 建立块段模型 四、地质品位模型计算过程 矿体剖面轮廓线的圈定 矿体剖面轮廓线 冬瓜山矿体实体模型图 矿体块段模型 四、地质品位模型计算过程 块金、基台选变异函数结果； 理论上搜索半径也应取变异函数变程； 但当模型较大时（如整个矿床），不能全部估值的情况。此时搜索半径可从变程开始，逐渐增大依次估值，可得到不同准确度的结果； 实际应用中最好结合生产建较小的能满足月、季度生产计划需要的模型，并推估。 选一个有准确矿量统计结果的采场或区域，计算变异函数并估值比较，为以后变异函数参数计算提供参考。 品位估值 地质数据库建立 原始样基本统计 样品组合 组合样统计分析 变异参数确定 变异函数计算 交叉验证 建立块段模型 四、地质品位模型计算过程 谢谢！ * * * * * * * * * * * * * * * * * 地质统计学原理与估值方法 中南大学数字矿山研究中心 长沙迪迈信息科技有限公司 冯兴隆 博士 2009年9月 算术平均法，将全区总面积乘以各钻孔矿体厚度的平均值来计算储量，其实质就是把整个复杂形状的矿体当做一个理想的、具有相同厚度的板状体； 地质块段法，实质上是把复杂形态的矿体当做一组密集的、大小不等的柱体，计算各个柱体储量之和； 多角形法，在储量计算平面图上，以各个钻孔见矿点为中