地质统计学原理及其在矿床建模与储量估算中的应用.pptx

《地质统计学原理及其在矿床建模与储量估算中的应用.pptx》由会员分享，可在线阅读，更多相关《地质统计学原理及其在矿床建模与储量估算中的应用.pptx（97页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、地质统计学原理及其在矿床地质统计学原理及其在矿床建模与储量估算中的应用建模与储量估算中的应用矿床品位建模及其应用需求矿床品位建模及其应用需求矿体表面模型（矿化边界）矿体表面模型（矿化边界）矿体表面模型（矿化边界）矿体表面模型（矿化边界）矿床品位模型矿床品位模型矿床品位模型矿床品位模型勘探线剖面品位分析勘探线剖面品位分析品位吨位曲线分析品位吨位曲线分析矿床品位建模及储量估算流程矿床品位建模及储量估算流程组合样品组合样品分析样品分析样品确定矿床块确定矿床块体模型参数体模型参数选择选择插值类型插值类型设置插值参数等设置插值参数等确定搜索邻域确定搜索邻域精度验证精度验证满意满意估值估值矿床品位模型矿床

2、品位模型否是回顾：地理学第一定律及应用回顾：地理学第一定律及应用地理学第一定律地理学第一定律：距离越近，两点的地理现象相似性越大距离越近，两点的地理现象相似性越大逐点移面内插逐点移面内插：以待插点为中心，确定一个邻域范围，用该邻域内的采样点计算内插点的高程值。反距离加权平均法反距离加权平均法内容介绍内容介绍地质统计学简介地质统计学简介区域化变量变差函数建模克里格插值算法矿体储量估算应用历史背景与产生历史背景与产生为解决矿床从普查勘探、矿山设计到矿山开发整为解决矿床从普查勘探、矿山设计到矿山开发整个过程中各种储量计算和误差估计问题发展起来个过程中各种储量计算和误差估计问题发展起来的。的。地质统计

3、学是数学地质的重要分支，它首先由地质统计学是数学地质的重要分支，它首先由 DGDG克立格克立格（KrigeKrige）工程师在南非的金属矿产）工程师在南非的金属矿产储量计算中使用，后由法国储量计算中使用，后由法国马特隆马特隆（GMathreonGMathreon）教授领导的小组对此作了深入的）教授领导的小组对此作了深入的研究并系统地总结出地质统计学的理论和方法。研究并系统地总结出地质统计学的理论和方法。地质统计学定义地质统计学地质统计学(Geostatistics)是以是以区域化变量区域化变量理论理论作为理论基础，以作为理论基础，以变差函数变差函数作为主要工具，对既作为主要工具，对既具有随机性

4、又具有结构性的变量（如品位值）进具有随机性又具有结构性的变量（如品位值）进行研究的科学。其核心即行研究的科学。其核心即“克里格法克里格法”，它是一，它是一种无偏的最小误差的储量计算方法。种无偏的最小误差的储量计算方法。区域化变量区域化变量变差函数变差函数克里格估值克里格估值与传统储量估算方法相比从传统方法把部分钻孔品位当作一个块段的品位，从传统方法把部分钻孔品位当作一个块段的品位，从而使高品位估计偏高，低品位估计偏低，而且从而使高品位估计偏高，低品位估计偏低，而且没有考虑矿石品位的空间变异性没有考虑矿石品位的空间变异性，在计算块段平，在计算块段平均品位时，每个样品的贡献仅仅是若干个几何因均品位

5、时，每个样品的贡献仅仅是若干个几何因素。素。地质统计学方法避免了传统方法的两个缺陷。其地质统计学方法避免了传统方法的两个缺陷。其加权因子是以加权因子是以矿床的各个方向变差函数的参数矿床的各个方向变差函数的参数为为基础计算出来的基础计算出来的,这种加权方法充分考虑了矿体这种加权方法充分考虑了矿体形态的空间变化及其品位空间变化特征形态的空间变化及其品位空间变化特征,并且采并且采用用了无偏的了无偏的、误差最小误差最小的数理统计方法计算样品的数理统计方法计算样品的加权因子和块段的品位。的加权因子和块段的品位。地质统计学的发展完善的理论基础完善的理论基础 基本概念区域化变量 基本工具变差函数 基本假设本

6、征假设 基本方法克里格法方法与技巧不断涌出方法与技巧不断涌出 析取克里格、多元高斯克里格和各种条件模拟技术的应用和发展地质统计学的软件包及应用软件不断推出地质统计学的软件包及应用软件不断推出 美国斯坦福大学的GSLIB软件包 挪威ODEN公司的STORM随机建模软件 加拿大的Geostat地质统计学软件 澳大利亚的Surpac VisionMicromine矿山工程软件内容介绍地质统计学简介区域化变量区域化变量变差函数建模克里格品位估值矿体储量估算应用区域化变量G.G.马特隆定义区域化变量是：一种在空间上具有马特隆定义区域化变量是：一种在空间上具有数值的实函数，它在空间的每一个点取一个确定数值

7、的实函数，它在空间的每一个点取一个确定的数值，即当由一个点移到下一个点时，函数值的数值，即当由一个点移到下一个点时，函数值是变化的是变化的.特征：特征：随机性随机性和和结构性结构性 随机性随机性 结构性结构性区域化变量从地质及矿业角度来看，区域化变量具有如下性质：（1 1）空间局限性）空间局限性：即它被限制在一个特定的空间（如一个矿体内）；该空间称为区域化的几何域；区域化变量是按几何支撑定义的。（2 2）连续性）连续性：不同的区域化变量具有不同的连续性，这种连续性是通过相邻样品之间的变差函数来描述的。（3 3）异向性）异向性：当区域化变量在各个方向上具有相同的性质时称各向同性，否则称各向异性。

8、（4 4）相关性）相关性：一定范围内、一定程度上的空间相关性，当超出这一范围后相关性减弱以至消失。（5）对于任一区域化变量而言，特殊的变异性是叠加在一般特殊的变异性是叠加在一般规律之上规律之上。内容介绍地质统计学简介区域化变量变差函数建模变差函数建模克里格插值算法矿体储量估算应用变差函数建模 为表征一个矿床金属品位等特征量的变化，经典统计学通常采用均值、方差均值、方差等一类参数，这些统计量只能概括该矿床中金属品位等特征量的全貌，却无法反映局部范围和特定方向局部范围和特定方向上地质特征的变化。地质统计学引入变差函数这一工具，它能够反映区域化变量的空间变化特征相关性和随机性，特别是透过随机性反映区

9、域化变量的结构性，故变差函数又称结构函数结构函数。变差函数定义我们可以把一个矿床看成是空间中的一个域，如图中 为沿 方向被矢量 分割的两个点，其观测值分别为 及 ，该两者的差值差值 就是一个有明确物理意义的结构信息，因而可以看成是一个变量。区域化变量 在空间相距 的任意两点 和 处的值 与 差的方差之半定义为区域化变量区域化变量 的变差函数的变差函数，记为变差函数定义定义：在任一方向，相距的两个区域化变量和的增量的方差的一半。公式：变差函数值与区域化变量位置无关二阶平稳假设二阶平稳假设和本征假设本征假设二阶平稳假设当区域化变量满足下列两个条件时，称该区域化变量满足二阶平稳：（）在整个研究区内，

10、区域化变量 的期望存在且等于常数：（常数）（）在整个研究区内，区域化变量的空间协方差函数存在且平稳：当时，上式变成：即它有有限先验方差。本征假设当区域化变量的增量满足下列两个条件时，称该区域化变量满足本征假设：（）在整个研究区内，区域化变量 的增量的期望为：（）对于所有区域化变量的增量的方差函数存在且平稳：即要求的变差函数存在且平稳实验变差函数计算其中:=两个样本点间的距离 =样本点属性值(位置 )=样本点属性值(位置 )=样本点数变差函数计算公式：变差函数计算实例某地区规则采样数据，数据为属性值，样本间距为100米。实验变差函数计算实例图中表示的是东西方向，相距为100米的样本点对。实验变差

11、函数计算实例通过变差函数计算公式得到东西方向上，滞后距为100米的变差函数值。实验变差函数计算实例变差函数图：滞后距100米的变差函数点024681012141618200100200300400500滞后距变差函数实验变差函数计算实例相距为200米的样本点对。实验变差函数计算实例滞后距为200米的变差函数值。变差函数计算实例变差函数图：滞后距200米的变差函数点024681012141618200100200300400500滞后距变差函数变差函数计算实例变差函数图：滞后距300米、400米的变差函数点024681012141618200100200300400500滞后距变差函数变差函数计

12、算实例计算南北方向滞后距为100米、200米和300米的变差函数。实验变差函数计算实例南北方向400m点数过少，不参与计算。滞后距东西方向南北方向1001.465.352003.39.873004.3118.884006.7变差函数值实验变差函数计算实例变差函数图：东西方向和南北方向024681012141618200100200300400500滞后距变差函数东西方向南北方向实验变差函数计算-距离和角度容差对于不规则采样点不规则采样点：沿某一特定方向和特定滞后距上并没有足够的样本点采用距离和角度容差解决该问题实验变差函数计算步长：步长：4m4m步长容差：步长容差：2m2m方位角：方位角：60

13、60倾角：倾角：0 0方位容差：方位容差：22.522.5倾角容差：倾角容差：22.522.5水平带宽：水平带宽：5m5m垂直带宽：垂直带宽：5m5m实验变差函数计算(3D)变差函数的计算过程是由系统自行完变差函数的计算过程是由系统自行完成的，而合适的参数大小将直接影响计算成的，而合适的参数大小将直接影响计算结果的好坏。结果的好坏。关于参数的选取关于参数的选取实验变差函数参数选择步长大小的选择步长大小的选择：步长间距太小步长间距太小步长间距较合适步长间距较合适实验变差函数参数选择步长个数的选择步长个数的选择：原则:步长大小*步长个数=研究区域长度的一半步长总间距步长总间距理论变差函数实验变差函

14、数并不能定量的反映数据空间相关性，需要对实验变差函数进行拟合得到理论变差函数。理论变差函数三参数：块金值块金值/基台值基台值/变程变程(基台值=先验方差)Samples not spatially correlatedSamplesSpatially Correlated基台值基台值变程变程块金值块金值0.(h)gh样本空间相关样本空间不相关理论变差函数模型Samples not spatially correlatedSamplesSpatially Correlated球状模型线性模型指数模型高斯模型球状模型球状模型公式：接近原点处，变差函数呈线性形状线性形状，在变程处达到基台值。原点处变

15、差函数的切线在变程的2/3处与基台值相交。实验变差函数在大多数情况下可以拟合成球状模型。因此，球状模型是应用最广的一种变差函数模型球状模型是应用最广的一种变差函数模型。指数模型指数模型公式：变差函数渐近地逼近基台值，在实际变程 处，变差函数为0.95，模型在原点处为直线。在原点处连续性最好，是一种较稳定的模型。高斯模型高斯模型公式：变差函数渐近地逼近基台值，在实际变程 处，变差函数为0.95，模型在原点处为抛物线。为一种连续性好但稳定性较差的模型。变差函数拟合用球状模型、指示模型或高斯模型对实验变差函数进行拟合。得到块金值、基台值和变程三个参数。变差函数拟合球状模型变程为4141m，指数模型变

16、程为5823m，高斯模型变程为2884m观察图形：高斯模型拟合最好，其次是球状模型。根据实际情况确定变差函数类型，结果因人而异。变差函数拟合过程几何各向异性 基台值相同变程不同 在不同的方向具有相同的变异程在不同的方向具有相同的变异程度度（基台值相同）（基台值相同）但具有不同的但具有不同的连续程度连续程度（变程不同）（变程不同）为几何各为几何各向异性。向异性。带状各向异性 基台值不同变程可同可不同 在一些不同的方向上具有不同的变异在一些不同的方向上具有不同的变异程度程度（基台值不同）（基台值不同）连续程度连续程度（变程）（变程）可以相同也可不同为带状各向异性可以相同也可不同为带状各向异性。变差

17、函数结构套合不同方向结构套合不同方向结构套合几何各向异性基本思路为通过线性变换将各向异性的坐标向量 转化为各向同性的新坐标向量 设这个线性变换为 ，其中对于各向同性模型，其中 对于几何各向异性变差函数 ，变化为矩阵形式变差函数结构套合不同方向结构套合不同方向结构套合带状各向异性对于带状各向异性，采用分块处理的方法。具体的变差函数模型公式为 ，其中对于 做和几何各向异性相同的处理，对于 做如下处理 ，对于 做如下处理：总的来说，对于带状各向异性的处理方法是将其看作是几何各向异性进行坐标变换后，再分别对次轴和垂直轴方向上多出的基台值进行叠加处理。各向异性椭球各向异性椭球：各向异性椭球：主轴变程次轴

18、变程垂直轴变程方位角倾角旋转角度几何各向异性结构套合变差函数表面图All points that fall in the block are paired with the point at(x,y)to create the variogram maps.The size of the block is the lag size.变差函数建模变差函数是区域化变量空间变异性的一种度量反映了空间变异程度随距离而变化的特征可定量的描述区域化变量的空间相关性地理学第一定律：地理学第一定律：距离越近的点相似性越大距离越近的点相似性越大内容介绍地质统计学简介区域化变量变差函数建模克里格插值算法克里格插值

19、算法矿体储量估算应用 如果要估算 的值，一般情况下 的值应该是 的平均值，并且 随着距离 的增大而减小。克里格插值算法克里格插值算法建立在变差函数及结构分析理论克里格插值算法建立在变差函数及结构分析理论之上之上适用条件是变差函数及相关分析的结果表明样品适用条件是变差函数及相关分析的结果表明样品间存在空间相关性间存在空间相关性其实质是利用其实质是利用区域化变量的原始数据区域化变量的原始数据和和变差函数变差函数的结构特点的结构特点，对未采样点的区域化变量的取值进，对未采样点的区域化变量的取值进行线性、无偏、最优估计。行线性、无偏、最优估计。克里格插值过程组合样品分析样品确定块体模型参数选择克里格类

20、型计算实验变差函数并拟合确定搜索邻域交叉验证满意估值品位模型否是块体模型定义 将整个研究区域划分为多个规则小块，分别对小块属性进行估值。起始点坐标起始点坐标 块大小块大小 块个数块个数搜索邻域确定搜索椭圆直接定义点数交叉验证交叉验证的原理为将原始的样品点去除，然后采用原始样品点周围的点来进行克里格估值得到原始样品的估计值，最后做出原始样品和估值样品的散点图，并对估值误差进行统计。克里格插值算法从矿业上的术语具体来说，它是根据一个从矿业上的术语具体来说，它是根据一个待估块待估块段邻域内的若干信息样品的品位数据段邻域内的若干信息样品的品位数据，在考虑了，在考虑了这些这些样品的形状样品的形状、大小大

21、小及及相互位置关系相互位置关系，它们与，它们与待估块段相互之间的空间位置等几何特征，以及待估块段相互之间的空间位置等几何特征，以及品位的变差函数模型所提供的结构信息之后，为品位的变差函数模型所提供的结构信息之后，为了对该块段品位作出一种了对该块段品位作出一种线性、无偏、最小估计线性、无偏、最小估计方差的估计方差的估计而对每个样品值分别赋予一定的权系而对每个样品值分别赋予一定的权系数，最后进行数，最后进行加权平均加权平均来估计该块段品位的方法。来估计该块段品位的方法。克里格插值算法克里格插值算法克里格插值算法:B.L.U.E-B.L.U.E-best,linear,unbiased estima

22、tor best=最小估计误差 linear=线性估值方式（同距离反比估值）unbiased=无偏估计，估计误差之和为0 estimator=估值方法克里格插值算法定义其中 =待估点位置和其中一个邻接点 位置 =估算未知点 用到的邻接点个数 =和 对应的预测平均值 =对应的克里格权重 克里格插值公式克里格插值算法基础关键在于确定邻接权重最小方差限制条件无偏估计限制条件克里格插值公式克里格插值类型最常用克里格的三种类型最常用克里格的三种类型 简单克里格简单克里格 普通克里格普通克里格 泛克里格泛克里格 其区别在于 的确定方式不同 非线性克里格非线性克里格 指示克里格（Indicator Krig

23、ing）多元高斯克里格（Multi-Gauss kriging)协克里格（Cokriging）块克里格块克里格(Block kriging)(Block kriging)简单克里格插值应用实例六个样本点数据，给出样本点间距、样本点属性值和变差函数变差函数模型：球状模型，块金值0，基台值0.78，变程4141m Pnt1Pnt2Pnt3Pnt4Pnt5Pnt6Pnt1018973130244114001265Pnt2189701281145619702280Pnt3313012810152328003206Pnt4244114561523015231970Pnt51440197028001523

24、0447Pnt612652280320619704470简单克里格插值应用实例由简单克里格插值公式得 简单克里格插值应用实例由简单克里格插值公式得 简单克里格插值应用实例已知该采样数据平均值为14.70六个采样点数据的属性值分别为 13.84,12.15,12.87,12.68,14.41,14.59 由 得 简单克里格插值应用实例普通克里格插值普通克里格插值：未知普通克里格估值公式为由无偏最优估计限制条件无偏最优估计限制条件，构建拉格朗日函数其中拉格朗日乘数法用“拉格朗日乘数法”求函数f(x,y,z)在条件(x,y,z)=0下的极值方法（步骤）是：1.做拉格朗日函数L=f(x,y,z)+(x

25、,y,z),称拉格朗日乘数 2.求L分别对x,y,z,求偏导,得方程组,求出驻点P(x,y,z)如果这个实际问题的最大或最小值存在,一般说来驻点唯一,于是最值可求.普通克里格插值求偏导分别得到下列公式得到求得到求 的方程组的方程组估值方差计算公式估值方差计算公式指示克里格估值在地质、物化探数据处理及矿产储量计算中影响计算精度的因素有很多，但主要有以下几个问题：(1)特异值的出现，所谓特异值是指那些比全部数值的平均值或中位数高得多的数值，它既非分析误差所致，也非采样方法等人为误差引起。而是实际存在于所研究的母体之中。这些特异值只占全部数据的极少部分，但却控制了总金属资源量的很大比例。(2)在一个

26、研究区域或一个矿床中存在几个不同类型的矿化作用，这也影响了品位和储量的精确估计。为解决上述问题，指示克里格法应运而生，它是在不必去掉重要而实际存在的高值数据的条件下来处理不同的现象，而且给出在一定风险概率条件下未知量的估计值及空间分布。指示克里格估值指示克里格是一种非参数地质统计学方法。它是指示克里格是一种非参数地质统计学方法。它是根据一系列的临界值，例如边界品位，先对原根据一系列的临界值，例如边界品位，先对原始数据如下进行转换始数据如下进行转换然后对转换后的数值求变差函数、进行克里格估然后对转换后的数值求变差函数、进行克里格估值。值。总结简单克里格（整体平稳）简单克里格（整体平稳）普通克里格

27、（局部平稳）普通克里格（局部平稳）泛克里格（整体存在趋势）泛克里格（整体存在趋势）指示克里格（数据不平稳存在极值）指示克里格（数据不平稳存在极值）协克里格（存在辅助变量）协克里格（存在辅助变量）已知常量但未知其中 常量但未知几点注意内容几点注意内容变差函数参数变差函数参数 块金值块金值：块金值越小，距离越近的点越重要，这样会导致权值的变化范围变大（从负值到大于1的值变化），使数据出现异常。块金值越大，估值结果越平滑。变程：变程：变程小于任意两点之间的距离，变差函数为纯块金模型，随着变程增大，已知点的位置以及丛聚性变的重要。如果变程很大存在基台值，则变差函数相当于纯块金模型，如果变程很大但不存在

28、基台值，则变差函数为线性模型。比例比例：比例大小只与克里格估值方差有关，与估值结果无关。形状形状：球状模型与指数模型在原点位置接近线性关系，指数模型在原点处更加陡峭一些，与变程较小的球状模型相似。高斯模型在原点处为抛物线形式，这种模型要求原始变量要有很高的连续性，否则会出现大量负的权值情况，使估值结果极不稳定。各向异性各向异性：以各个方向变程不同反映出来，可以理解为对原始数据坐标位置的变化，通过几何校正转化为各向同性情况。几点注意内容屏蔽效应屏蔽效应第一个点在任何情况下获得的权值更大一些。第二个点的权值小于第第一个点在任何情况下获得的权值更大一些。第二个点的权值小于第一个点，在克里格估值过程中

29、，第二个点的权值可能变为负值，这样一个点，在克里格估值过程中，第二个点的权值可能变为负值，这样第二个点的信息变得重要。例如：当一个负值权重可能隐含某种趋势第二个点的信息变得重要。例如：当一个负值权重可能隐含某种趋势（如果第一个点比第二个点要小，那么估值点可能小与这两个点。（如果第一个点比第二个点要小，那么估值点可能小与这两个点。随着块金值的增大屏蔽效应减弱。所用样品点对估值结果作用基本相随着块金值的增大屏蔽效应减弱。所用样品点对估值结果作用基本相同。同。0.22 0.2 0.18几点注意内容几点注意内容搜索策略搜索策略 最小值：最小值：作用是限制克里格估值过程中数据个数，如果在搜索椭球范围内，

30、数据个数过少则不进行估值，否则估值结果可能存在很大误差。最大值最大值：估值过程中用到的数据越多，估值结果越平滑。由于数据的局部平稳性，数据个数不能太多，否则会使不存在相关性的数据仍然参与到估值过程中。搜索半径搜索半径：搜索半径越大估值结果越平稳，搜索半径与最大值相互影响。搜索半径设置与变差函数变程相关。几点注意内容几点注意内容克里格估值类型克里格估值类型 简单克里格简单克里格：简单克里格需要一个固定的平均值。在估值过程中平均值被赋予一个相当大的权值，从而导致估值结果过渡平滑。普通克里格：普通克里格：与简单克里格类似，但平滑性减弱。普通克里格估值权值之和为1，可以合理的减少丛聚效应。泛克里格：泛

31、克里格：与普通克里格类似，但估值结果边缘会存在一个奇异值，这是由趋势面构造的不合理造成的。块克里格：块克里格：由于块克里格估值的平滑效应会使克里格估值结果方差减小。随着块的增大，估值结果越来越平滑。内容介绍地质统计学简介区域化变量变差函数建模克里格品位估值矿体储量估算应用矿体储量估算应用矿床品位建模及储量估算流程矿床品位建模及储量估算流程组合样品组合样品分析样品分析样品确定矿床块确定矿床块体模型参数体模型参数选择选择克里格类型克里格类型计算变差计算变差函数并拟合函数并拟合确定搜索邻域确定搜索邻域精度验证精度验证满意满意估值估值矿床品位模型矿床品位模型否是克里格法资源储量估算流程克里格法资源储量

32、估算流程关键步骤关键步骤否否否否否否否否是是是是是是是是数据分析数据分析数据分析数据分析数据变换数据变换数据变换数据变换选择泛克里格估值选择泛克里格估值选择泛克里格估值选择泛克里格估值选择其它克里格估值选择其它克里格估值选择其它克里格估值选择其它克里格估值剔出趋势剔出趋势剔出趋势剔出趋势导入组合样品数据导入组合样品数据导入组合样品数据导入组合样品数据是否服从是否服从正态分布正态分布是否存是否存在趋势在趋势数数数数据据据据分分分分析析析析与与与与变变变变换换换换克里格法资源储量估算流程克里格法资源储量估算流程关键步骤关键步骤计算实验变差函数并计算实验变差函数并计算实验变差函数并计算实验变差函数并

33、绘制变差函数表面图绘制变差函数表面图绘制变差函数表面图绘制变差函数表面图是否服是否服是否服是否服从各向从各向从各向从各向异性异性异性异性单方向变差单方向变差单方向变差单方向变差函数拟合函数拟合函数拟合函数拟合确定最大确定最大确定最大确定最大变程方向变程方向变程方向变程方向多方向变差多方向变差多方向变差多方向变差函数拟合函数拟合函数拟合函数拟合是是是是否否否否球球状状模模型型拟拟合合变差函数拟合变差函数拟合单方向变差函数拟合单方向变差函数拟合多方向变差函数拟合多方向变差函数拟合指指数数模模型型拟拟合合高高斯斯模模型型拟拟合合几几何何各各向向异异性性带带状状各各向向异异性性变变变变差差差差函函函函

34、数数数数计计计计算算算算及及及及拟拟拟拟合合合合数据分析与变换结果数据分析与变换结果数据分析与变换结果数据分析与变换结果克里格法资源储量估算流程克里格法资源储量估算流程关键步骤关键步骤设置克里格参数设置克里格参数交叉验证交叉验证克里格估值计算克里格估值计算否否是是设置搜索椭球设置搜索椭球结果结果是否满意是否满意数据分析与变换数据分析与变换变差函数拟合变差函数拟合参参数数设设置置及及交交叉叉验验证证矿体储量估算应用矿体储量估算应用以某矿区数据为例介绍地质统计学法在储量估算以某矿区数据为例介绍地质统计学法在储量估算中的应用中的应用在进行矿体储量估算前需要首先对钻孔坐标和品在进行矿体储量估算前需要首

35、先对钻孔坐标和品位进行检查保证数据的正确性。位进行检查保证数据的正确性。该矿区矿体边界品位为该矿区矿体边界品位为0.1g/t0.1g/t确定矿体边界以0.1作为边界品位进行单工程矿体圈定三维矿体连接划分组合样原始样品长度分析原始样品长度平均原始样品长度平均值及中值接近值及中值接近2 2米，米，故选择故选择2 2米作为组合米作为组合样长度样长度样品组合采用长度样品组合采用长度加权平均方式加权平均方式样长相等可避免插样长相等可避免插值偏差值偏差数据分析对组合样数据进行直方图分析变差函数计算及拟合变差函数反映样品品位的空间相关性变差函数反映样品品位的空间相关性计算三个互相垂直方向的变差函数，对三个变

36、差计算三个互相垂直方向的变差函数，对三个变差函数进行拟合，根据各向异性进行结构套合。函数进行拟合，根据各向异性进行结构套合。根据矿区矿体走向和倾向确定主轴、半轴和次轴根据矿区矿体走向和倾向确定主轴、半轴和次轴方向方向变差函数计算及拟合主轴方向主轴方向变差函数及拟合结果变差函数拟合参数需要手动确定变差函数计算及拟合半轴方向半轴方向变差函数及拟合结果变差函数拟合参数需要手动确定变差函数计算及拟合次轴方向次轴方向变差函数及拟合结果变差函数拟合参数需要手动确定结构分析与理论模型套合 利用球状模型对变差函数进行理论模型的拟合利用球状模型对变差函数进行理论模型的拟合利用球状模型对变差函数进行理论模型的拟合

37、利用球状模型对变差函数进行理论模型的拟合获得的三个方向的理论模型，对三个方向拟合的理获得的三个方向的理论模型，对三个方向拟合的理获得的三个方向的理论模型，对三个方向拟合的理获得的三个方向的理论模型，对三个方向拟合的理论模型分析，进行结构套合。论模型分析，进行结构套合。论模型分析，进行结构套合。论模型分析，进行结构套合。主轴：基台值为主轴：基台值为主轴：基台值为主轴：基台值为 1 1，变程为，变程为，变程为，变程为7070半轴：基台值为半轴：基台值为半轴：基台值为半轴：基台值为 1 1，变程为，变程为，变程为，变程为6 60 0次轴：基台值为次轴：基台值为次轴：基台值为次轴：基台值为 1 1，变

38、程为，变程为，变程为，变程为6 60 0看出看出看出看出:矿体在各方向上呈现变异程度相同而变异连续性不同的矿体在各方向上呈现变异程度相同而变异连续性不同的矿体在各方向上呈现变异程度相同而变异连续性不同的矿体在各方向上呈现变异程度相同而变异连续性不同的几何异向性特征，即三方向基台值相同、变程不同。几何异向性特征，即三方向基台值相同、变程不同。几何异向性特征，即三方向基台值相同、变程不同。几何异向性特征，即三方向基台值相同、变程不同。矿床品位模型块体划分将空间矿体按一定的间距划分为一些连续的小立方块，然后用克里格法对小方块进行赋值。块体间距根据矿体实际情况确定如果有矿体边界约束，可对矿体边界进行细分，使块体模型与实体模型更接近搜索邻域设置用搜索邻域内的样本点进行克里格品位插值用搜索邻域内的样本点进行克里格品位插值定义一个搜索椭球，定义一个搜索椭球，一般情况下参数与变差函数一般情况下参数与变差函数变程相同变程相同由于样品沿钻孔方向密集，为保证估值结果准确，由于样品沿钻孔方向密集，为保证估值结果准确，一般采用一般采用八分圆搜索方式八分圆搜索方式。普通克里格品位估值结果估值结果保存到access数据库中普通克里格品位估值结果切面图切面图三维图形显示三维图形显示