地质统计学在地质储量计算中的应用.docxVIP

地质统计学在地质储量计算中的应用 地质统计学是1960年和1970年发展起来的一门新兴的边缘学科，主要结合数学、地质和采矿。它最早由南非地质学家D·G·克立格提出,后来经法国学者G·马特隆进行了系统地研究,使地质统计学成为一门具有坚实理论基础、合理工作方法的独立学科。地质统计学在我国20多年的生产实践表明,它有着广阔的发展前景和强大的生命力,不但可以用于矿产储量计算、矿山设计、矿山生产、矿山地质科研工作,而且还是用作估计各种自然资源的有效方法和手段。 1 现代矿产资源估计的基础 地质统计学是以矿体参数的空间相关性和随机性为前提,以区域化变量为研究对象,以半变异函数为基础,以计算机为运算工具,在估计方差最小和无偏条件下,建立一组待估单元块周围样品最佳权系数方程组,以精确估计矿石储量和质量的现代储量计算方法。在储量计算中,充分考虑到样品形状、大小和空间相互位置等对待估块段的影响,使计算精度大大提高,并且可得到一个相应的估计误差。 1.1 空间变换特征 区域化变量是地质统计学的研究对象,它是具有数值的空间位置的函数,即由一点移往下一点时,函数值是变化的,并具有明显的不同程度的连续性。因为这种变量既有随机性又有结构性,在地质统计学中,常用变异函数和变异曲线来表征区域化变量的空间变化特征与变化程度。如矿石的品位、矿体的厚度、矿石内有害组分的含量等都是区域化变量。 1.2 x型积分式 变异函数是区域化变量增量平方的数学期望,也就是区域化变量增量的方差。但是通常为了数学运算的方便,而使用半变异函数,半变异函数也简称变异函数。半变异函数可分为理论半变异函数,局部半变异函数和实验半变异函数。 理论半变异函数,在内蕴假设条件下定义为: γ(h)=E{[Z(x+h)-Z(x)]2}= 12V∫v[Z(x+h)?Z(x)]2dx12V∫v[Ζ(x+h)-Ζ(x)]2dx 式中积分式是在一个无限大的V几何区域内,无限多对样品增量的平方和。实际上除非矿床全部开采完毕,否则,是不可能得到的。 局部半变异函数,是在一个有限的几何域L中,无限多对点样品增量平方的均值,可写成:γ′(h)=12(L?h)∫L?h0[Z(x+h)?Z(x)]2dxγ′(h)=12(L-h)∫0L-h[Ζ(x+h)-Ζ(x)]2dx,式中h为滞后距,L为几何域。实际上这样的变异函数也是不易得到的。 实验半变异函数,是指在有限的几何区域内,有限多对样品增量平方的均值,可写成:γ?(h)=12NΣi=1N[Z(xi+h)?Z(xi)]2γ*(h)=12ΝΣi=1Ν[Ζ(xi+h)-Ζ(xi)]2。在实际工作中,只有实验半变异函数才可以得到。实验半变异函数是局部半变异函数的估计值。 1.3 空间分布特征 变异曲线可以半变异函数γ(h)为纵轴,以滞后距h为横轴,用通常的作图方法得到(见图1)。 图1是一个理想化的变异曲线图,它清晰地反映了区域化变量的空间分布特征。图中的C0称为块金效应值,它表示h很小时两点间品位的变化;a称为变程,当h≤a时,任意两点间的观测值有相关性,这个相关性随h的变大而减小,当ha时就不再具有相关性,a的大小反映了研究对象(如矿体)中某一区域化变量(如品位)的变化程度;C称为总基台值,它反映某区域化变量在研究范围内变异的强度。 1.4 变异函数理论模型 根据区域化变量特点所绘制的变异函数曲线是一种根据样本实测数据作出的实验变异函数,由于样品数据较少,实际上只是一种锯齿状的非光滑曲线。因此这种变异函数曲线作出后,还必须用一个适当的圆滑曲线或直线对它进行拟合,并用一个特定的函数来描述它,这就是变异函数理论模型,用以反映区域化变量的空间变化特征。理论模型拟合的结果将直接参与克立格储量计算或其他地质统计学研究。最常用的理论模型是球状模型、高斯模型及指数模型。 1.5 地质储量的估计方法 根据不同研究对象和目的,克立格方法有不同的形式,这里只介绍最常用的普通克立格法。 克立格法是用样品的加权平均求估计值的。加权平均就是对于任一待估块段V的真实值ZV的估计值Z*v是通过该待估块段影响范围内n个有效样品值Zα(α=1,2,…,)的线性组合得到的,即:Zv=Σα=1nλαZαΖv=Σα=1nλαΖα,式中λα是加权系数,是各样品在估计Z*v时的影响大小的反映。 估计方法的好坏取决于如何选取和计算权系数λα。克立格法是以最优(方差最小)、无偏为条件来选取和计算权系数λα。它通过列出克立格方程组,解出方程组中的各权系数,然后赋与各样品值,最后进行加权平均来求得待估块段的品位。累加得到每一块段的储量,进而得到整个矿体的储量。 运用克立格法计算地质储量有两个重要的工作要做:第一,列出克立格方程组,解出克立格权系数λα;第二,求出这种估计的最小估计方差——克立格方差。 普通克立格方程组和克立