MICROMINE的普通克立格法进行资源评价的例子课件.ppt

在澳大利西亚联合储量委员会规范和联合国标准下利用MICROMINE的普通克立格法进行资源评价的例子,JORC规范和联合国标准,联合储量委员会规范（简称JORC规范）在澳大利西亚进行资源和储量评价的标准与规范联合国标准 中国所采用的标准，与JORC规范相似,对数据的说明,在建立模型之前对给定文件中出现的错误进行检查识别并加以改正，包括坐标和品位数据在建立模型之前，将所有文件(孔口文件，测斜和样品分析结果)中的样长与相对应的深度核准，并将诸如钻孔重复、钻孔缺失、长度区间不连续和其它错误一并修正,20%background cutoff,2 g/t 背景值,对样品品位进行直方图分析,拐点显示为2g/t,概率曲线,背景值与矿化的界线,概率曲线的拐点指示了背景与矿化的界线值是2g/t此品位值用于矿化边界的圈定,地表面与地质工程,参数说明,以2g/t为边界值，钻孔被赋与颜色，以此值圈定矿体边界可采厚度为1.1米，最低工业品位3g/t,钻孔平面视图及剖面线,钻孔剖面,钻孔解释,通过鼠标拉出一条剖面线并右击鼠标键后可以平面或剖面的形式显示钻孔分布，并可以定义剖面的窗口“厚度”剖面中的钻孔被赋予颜色，矿体的边界线被分成两个域，1号矿体和号矿体圈定的矿体轮廓线均在三维状态下捕捉在钻孔上,捕捉钻孔点到轮廓线上,地质模型,在两个矿脉中定义了四个域，仅因为稀疏的钻孔在东边，集中的钻孔在西边，区域明显是由于断层而发生了位移。第四个域是分散于大域的那些小块体,创建地质模型,地质模型,地质模型,地质模型,创建3D钻孔坐标,3D坐标是每个样段的中心点坐标，它要引用孔口控制文件和测斜文件,创建3D钻孔坐标,将矿脉赋值到样品分析文件中,三维实体或地质模型被赋值到样品分析文件中，这样用过滤器程序就可以在解释中只利用指定的品位数据观察被赋值字段,将矿脉赋值到样品分析品位文件中,将样品长度组合成等长度,组合样长度,直方图说明超过80%的样品长度为1米，所以选择了1米作为组合样长度组合区间头尾以矿脉边界为限，对品位进行样长加权平均样长相等是为了避免插值偏差,组合样长,左侧为没有组合时的状态，右侧是以1米为组合样长的结果,建模规则,为四个矿化域独立插值，只有存在于矿脉之内的品位值才用于插值计算普通克立格法的最小估值长度为0.5米，块段克立格法计算时北方向和东方向各5个点，垂向2个点,地质模型内的品位呈对数正态分布,搜索椭圆,搜索椭圆,搜索椭圆,搜索椭球,平行于矿脉定义了三个搜索椭球，走向300度，倾伏角30度，尺寸为：1 主矿脉：沿走向250米，垂直走向100米，25米厚2 小矿脉，走向150米，垂直走向90米，20米厚；3 号矿脉，走向200米，垂直走向120米，25米厚搜索数据时使用8个扇区，每个扇区最大6个点，每个椭圆最少2个点 搜索椭球用来包含相关数据和排除冗余数据,体重与块尺寸,在所有的块中缺省的体重值为2.8t/m3块的尺寸为东方向5米、北方向5米和垂高方向5米。次块为东2.5米、北2.5米和高1米。次块用于矿脉地质模型的几何边界与块模型相交的地段,变异函数,变异模型在变异函数的两个结构中，变异值与球形模型拟合很好沿走向和垂直走向的方向用不同的基台值体现了带状各向异性。在第二个结构中使用了一个较长的范围,块模型,块模型用来对品位的分布进行评价和对矿化点进行评估，只有建立了块模型后才可利用经济边界值进行经济评价，通常情况下，如果矿石的品位超过了这一边界值，矿石的价值才会大于开采成本，开采才是经济的。也可能用到另外一个边界值对入选原矿进行品位优化。,模型决策,创建模型时最重要的影响因素：1 地质因素是首要的；2 矿化域：首先要基于地质情况的前提下将不同的样本分开。3 搜索椭球的空间特征对于在预测中使用正确的数据是很重要的；4 插值方法的选择是基于以上因素确定的。,数据搜索,MICROMINE将搜索椭圆体分成若干扇区，且指定进行估值时用到的搜索点的数量，在每个扇中设定允许使用的点的最大数目。每个扇中可以规定最大点数和最小点数。,采样支撑,离散度 在普通克立格方法中，利用块内次级块品位值的平均值对在此块内的点进行估值。,各向异性,各向异性是指矿化在一个方向上的连续程度比另一个方向要强。用克立格法进行估值时，在变异函数的基础上为不同方向的样品数据施以不同的权重。变异函数显示了反应品位相关程度的变程。在三维模型中，由于各向异性体的存在，在三个方向的变异函数有不同的变程。,各向异性体类型,有两种各向异性体：1 带状异向性 基台值和变程在三个方向均不相同。在这种情况下，应该使用两个结构（模型），其中一个具有较长变程值用以将两个基台值变成同一水平，例如此长变程值可以为10000米。2 几何异向性 在三个方向上，基台值相同，但变程不同。,地质统计学,地质统计学是应用数学方法对品位分布进行建模分析，品位分布是在资源中受地质因素控制而明显存在的。,变异函数参数,钻孔走向方向和垂直于走向方向的变异函数进行对比，若它们具有明显不同的基台值时则代表了矿体的带状异向性。变异函数 可以将两个结构（模型）结合为一个球状模型对伽玛值进行拟合。,半变异函数,半变异函数的计算公式：N(h)(h)=1(xi yi)2 2N(h)i=1,（半变异函数模型拟合）结果,半变异函数以点形式显示伽玛值，X轴表示（样品之间的）距离，Y轴表示伽玛值。模型拟合是利用一个数学公式将这些点拟合成一条曲线。,变异函数模型,模型公式：球状模型a在此表示变程指数模型,变异函数规则,三维模型要求有三个互相垂直的变异函数，块金值和基台值应该是相同的，而变程可以不同。,模型拟合,对变异函数进行拟合是一个互动的过程，Noel Cressie 的统计学可以用最小二乘回归值检验拟合的质量，此值越小则说明拟合的质量越好。,Noel Cressie 统计,主要方向半变异函数,第二方向半变异函数,第三方向半变异函数,变异函数检验,交叉验证 交叉验证是在一批原始数据中去掉一个数据，通过运用变异函数加权计算，利用其周围的数据值为缺失数据的样品点估值。如果变异函数是稳定的，则所估计的值应和样品的数据十分接近。,普通克立格规则,普通克立格要求数据服从正态分布。如果在直方图和概率曲线图中有多个众数，则要用到多重指示克立格法。变异函数的模型参数和产状（如变异函数的方位角和倾伏角）可以保存在一个“（格式）Form”（在Micromine存储一系列设置的对象）中。在估值时可以根据存储的信息自动进行加权计算。,将（三个方向的）变异函数结合在一起进行估值,三个变异函数,如果块金值和基台值不同，Micromine不会将此“Form”保存，这是一项很重要的检验。,普通克立格,普通克立格:n(u)Z*OK(u)=OK(u)Z(u)=1 n(u)with OK(u)=1=1,空白模型,在进行估值之前应该在实体框架中创建空白块模型。为了对空白模型中所有的块进行估值，要将搜索椭球增大半径并运行多次。对块模型要进行检验，以确保模型中每个块都被赋与了估计的值,块模型水平视图,块模型倾斜视图,块模型水平视图,模型报告,矿石量=_x _y_z 次级块因子 X 体重品位=被估计的品位,资源总量 4,193,312 t 2.49 g/t总体模型报告,品位/吨位曲线,tonnes,g/t,校验,模型通过了全局校验和局部校验全局校验是指把全部原始数据与每个模型及整体模型的结果对比局部校验是指将块模型加载到钻孔上并且将原始数据与插值所得数据比较两种校验的结果令人满意,(全局),全局校验吨位,地质几何模型吨位11249295吨块模型吨位11208207吨差值0.3%,局部校验，剖面北4155050,模型剖面和品位,模型分级,应用“文件字段创建”命令，依据钻孔的数量、点数和方差等对模型的估值结果分级。在这个案例分析中使用克立格方差作为分级的依据。,克立格方差,没有单独存在的克立格方差，块模型中的对每个块进行估值时都会计算克立格方差和克立格标准差。此误差可以作为块估值精确度的评价依据，而且可以为此矿床进行分级将矿床分成具有不同可靠程度的部分,克立格方差结果,资源量与储量的区别,普通克立格模型得到的结果是资源量，只有对限定为可采的地区进行了开采设计和经济评价后，资源量才能变成储量。,总结,文件和文件校验、数据的三维可视化、交叉验证、搜索椭球的三维可视化、将不同方向的变异函数结合在一起并存储为“格式集”（formset）、检查是否对每个块 都进行了估值、模型的全局校验和局部校验。克立格可以平滑品位值，它考虑了矿化的各向异向性。而距离反比法则没有考虑地质因素和对矿体圈定所采用的方法，不能创建局部模型，而且往往会过高地估计品位值。,总结,这个演示反映了MICROMINE软件可以提供经过数学校正的结果,