金属矿床露天开采品位与储量计算

第一篇金属矿床露天开采第一章品位与储量计算第一节概述投资一个矿床开采项目，首先必须估算其品位和储量。一个矿床的矿量、品位及其空间分布是对矿床进行技术经济评价、可行性研究、矿山规划设计以及开采计划优化的基础，是矿山投资决策的重要依据。因此，品位估算、矿体圈定和储量计算是一项影响深远的工作，其质量直接影响到投资决策的正确性和矿山规划及开采计划的优劣。 从一个市场经济条件下的矿业投资者的角度看，这一工作做不好可能导致两种对投资者不利的决策：（1）矿体圈定与品位、矿量估算结果比实际情况乐观，估计的矿床开采价值在较大程度上高于实际可能实现的最高价值， 致使投资者投资于利润远低于期望值， 甚至带来严重亏损的项目。（2）与第一种情况相反，矿床的矿量与品位的估算值在较大程度上低于实际值，使投资者错误地认为在现有技术经济条件下，矿床的开采 不能带来可以接受的最低利润，从而放弃了一个好的投资机会。然而，准确地估算出一个矿床的矿量、品位绝非易事。大部分矿体被深深地埋于地下，即使有露头，也只能提供靠近地表的局部信息。进行矿体圈定和矿量、品位估算的已知数据主要来源于极其有限的钻孔岩心取样。已知数据量相对于被估算的量往往是一比几十万乃至几百万的关系，即对一吨岩心进行取样化验的结果，可能要用来推算几十万乃至几百万吨的矿量及其品位。可以不过分地说，矿量、品位的估算是世界上最大胆的外推。因此，矿体圈定与矿量、品位估算不仅是一项十分重要的工作，而且是一项极具挑战性的工作。做好这一工作要求掌握现代理论知识与手段，并应用它们对有限的已知数据进行各种详细、深入的定量、定性分析；同时也要求从事这一工作的地质与采矿工程师具有科学的态度和求实精神。本章将较详细地介绍当今世界上常用的矿量、品位估算方法，包括探矿数据的分析、处理和用于品位估值的剖面法、平面法及矿床模型法等。地质统计学作为品位估值的一种方法，从其诞生起就显示了强大的生命力，得到了越来越广泛的应用，本章对此给予较大的篇幅。本章的主要目的不是教会读者如何一步一步地应用所介绍的方法，对一个矿床进行矿量、品位估算，而是使读者了解这些方法的内涵，为读者提供在不同条件下应用最合理的分析、评价方法所需的知识基础。第二节探矿数据及其预处理一、钻孔取样用于矿体圈定与矿量、品位估算的数据主要来源于探矿钻孔的岩心图1-1钻孔与勘探线示意图表1—1钻孔岩芯信息记录钻孔号：zk10； 孔口坐标：6086.21E,6821.68N,205.01；设计深度：135M；实际深度：143.26M；开孔方位角： 开孔倾角：900；开孔日期:_1994年10月12日；终孔日期:_1994年10月23口 换层深度每层提取岩芯长度（M）每层岩芯采取率（%）岩石矿石描述白（M）至(M)共计（M）0.0013.9313.93表十层13.9330.6916.761.69.5云母石英岩：黄绿色，片状结构，主要组成矿物为石英（约25〜30%），云母（约40%）和角闪石（约25%）,其次有些磁铁矿。30.6943.0312.349.778.61阳起磁铁石英岩：钢灰色〜灰白色，细粒结构，主要组成矿物为石英（约40〜45%），磁铁石（约30〜35%），阳起石（约15〜20%）。::::::::::::取样。钻孔一般按照一定的网度布置在一些叫做勘探线的直线上（图1-1）。在钻孔过程中，每钻一定深度（一般在 3米左右）将岩心取出，做好标记后按顺序放在箱中供搬运、贮存和化验。地质人员对取出的岩心进行定性观察和简单的测试，以确定每一段岩心的主要物理特性，如岩心长度、岩性、颜色、图1-2钻孔柱状图硬度等，并记录下来，形成对钻孔穿过地段的地质特性的定性描述。表1-1是一个钻孔的岩心观测结果的部分记录表。为直观起见，常常把表中的数据和文字描述绘成钻孔柱状图（图1-2）。为了确定岩心的化学成分和品位，将岩心的一半送往化验室进行化验，另一半保存下来备用。样品的化验结果记录在如表 1-2岩心长度、岩性、颜色、图1-2钻孔柱状图的表中，或输入计算机的数据库中。手工记录时常将表1-1和表1-2合并为一个表，称为钻孔地质资料记录表。对所有钻孔的定性描述和取样化验结果构成了勘探区域的基本地质数据，这些取样化验数据是进行矿体圈定和矿量、品位估算的依据。在矿量和品位计算前，一般需要对取样数据进行预处理，包括样品组合处理和“极值”样品的处理。

表1-2钻孔岩芯取样化验结果记录钻孔号：zk10； 孔口坐标：6086.21E，6821.68N，205.01；设计深度：135M；实际深度：143.26M；开孔方位角： 开孔开孔日期：1994年10月12日；终孔日期：1994年10月23，倾角：900；日试样号 采样间隔 _化学分析结果1(%)备注白（M）至(M)共计（M）TFeFeSFe108330.6933.693.0029.8016.6022.50108433.6936.693.0032.2015.6025.10108536.6939.693.0032.9516.0028.00108639.6943.033.3426.4014.0021.00 : : : : : : : : : : : : : : 二、样品组合处理样品组合处理就是将几个相邻样品组合成为一个 组合样品，并求出组合样品的品位。当矿岩界限分明，且在矿石段内垂直方向上品位变化不大时，常常将矿石段内（即上下矿岩界限之间）的样品组合成一个组合样品（图13—3），这种组合称为矿段组合。组合样品的品位g是组合段内各样品品位的加权平均值，即(1-1)g=Yigl(1-1)TOC\o"1-5"\h\zii ii=1 i=1式中，z.为第i个样品的长度；g.为第i个样品的品位；n为矿石段内样品个数。’ ’式1-1中用的是长度加权，是最常用的方法。如果不同样品的比重相差较大，可以采用重量加权法。图1-3矿段样品组合示意图图图1-3矿段样品组合示意图图1-4台阶样品组合示意图对于拟用露天开采的矿床，更具实际意义的样品组合处理是 台阶样品组合，即把一个台阶高度内的样品组合成一个组合样品(图 1-4)。组合样品的品位为：g=^l.g^'H (1-2)i=1式中，H为台阶高度。当一个样品跨越台阶分界线时(如图1-4中第一和第五个样品)，在计算中样品的长度取落于本台阶的那部分长度(即图 1-4中的lt’’和15’’)，样品的品位不变。对钻孔取样进行台阶样品组合处理的意义在于：对取样数据进行统计学、地质统计学分析，以及利用取样值进行品位估值时，只有当每个样品具有相同的支持体，即每个样品的体积相同时，分析计算结果才有意义。露天开采在垂直方向上是以台阶为开采单元的， 一旦台阶的参考标高和台阶高度被确定，沿台阶高度无论品位如何变化，也无法进行选别开采。因此，在一个台阶高度内采用不同的取样品 位是毫无意义的。组合样品的品位较原样品品位变化小，在一定程度上减轻了“极值”品位对分析计算的影响，也使样品的统计分布曲线和半变异函数曲线(这些概念将在以后几节讲述)趋于规则。样品组合处理减少了样品总数，节省计算机内存和计算时间。三、极值样品(Outlier)处理极值样品是指那些品位值比绝大多数样品的品位(或样品平均品位)高出许多的样品，它们在贵重金属矿床较为常见。例如，在一金矿床取样1000个，经化验，这些样品的平均品位为 10克/吨，其中有十个样品的品位在100克/吨以上，这十个样品就可以被看成是极值样品。究竟品位比均值高出多少的样品算是极值样品，没有统一的、现成的标准，需视具体情况而定。极值样品虽然数量少，但对金属量影响大，为使品位的分析计算结果不致过分乐观，人们常常在实践中采用以下处理方法：⑴限值处理：即将极值样品的品位降至某一上限值。 比如在上述例子中，将所有高于100克/吨的样品的品位降至100克/吨。删除处理：即将极值样品从样本空间中删去，不参与分析计算。使用上述处理方法时应非常谨慎。虽然极值样品在数量上占样品总数的比例很小，但由于其品位很高，对矿石的总体品位和金属量的贡献值都很大。因此，不加分析地进行降值或删除处理会严重歪曲矿床的实际品位和金属含量，人为地降低矿床的开采价值。这一点可用下面的例子说明。假设对一金矿床进行钻探取样后得知，品位值服从对数正态分布(图1-5)。所有样品的平均品位为g=10克/吨，中值为m=3克/吨(即高于3克/吨和低于3克/吨的样品各占50%)；有1%的样品品位高于100克/吨。若将这1%的极值样品取出，单独计算其平均值，得 190克/吨。那么这1%的样品对矿床总金属量的贡献为(190X1%)/g=1.9/10=19%。也就是说，百分之一的数据量代表的是百分之十九的金属量！假如取边界品位为3克/吨(高于3克/吨为矿石，否则为废石)，矿石的平均品位(即高于3克/吨的那部分样品的平均品位)经计算为 16克/吨。如果把极值样品从样品空间删除，矿石的平均品位变为( 16X50%-190X1%)/(50%-1%)=12.45克/吨，也就是说，矿石品位被低估了22%。如果将极值样品进行限值处理，将其品位值降到 100克/吨，矿石的平均品位变为(12.45X49%+100X1%)/50%=14.2克/吨，也就是说，将矿石品位低估了11%。图1-5金矿取样品位对数正态分布示意图在正常、稳定的经济环境中，采矿的利润率也就是15%左右。因此，不加分析地将极值样品进行删除或限值处理，很可能将本来能够获取正常利润的矿床人为地变为没有开采价值，从而导致错误的投资决策。这对于一个在市场经济条件下，以盈利为主要目的的矿业投资者来说，无疑是一个重大的决策失误。这里必须澄清的是，极值样品是实实在在存在的有效样品，并不是指那些由于化验或数据录入错误造成的、具有“错误品位值”的样品。如果有根据认为某些样品的品位是错误的，将这些样品从样本空间中删除不仅是合理的而且是必要的。对极值样品的最理想的处理方法是，经过对探矿区域的地质构造和成矿机理进行深入分析，将这些样品的发生区域(或构造)划分出来，在进行品位与矿量的分析计算时，这些样品只参与其发生区域的品位与矿量计算，而不把它们外推到发生区域之外。但是在大多数情况下，由于钻孔网度大，已知的地质信息满足不了这种区域划分的要求。这时，可以将矿床看成是由两种不同的矿化作用形成的：样品中占绝大多数的“正常样品”可以看作是由主体矿化作用产生的样本空间；极值样品是由次矿化作用产生的样本空间。然后利用统计学方法计算出空间任一点属于每一类矿化作用的概率，再根据这些概率计算矿床的品位与矿量。这一方法超出了本书的范畴，有兴趣的读者可参阅Journel(1988)和Parker等人(1979)的论文。

第三节取样数据的统计学分析对取样数据进行上述的预处理以后，做一些统计学分析可以提供不少有关矿床的有用信息。因此统计学分析常常是取样数据分析的第一步。对数据进行统计学分析的主要目的是确定：品位的统计分布规律及其特征值；品位变化程度；样品是否属于不同的样本空间；根据样品的分布特征，初步估计矿床的平均品位以及对于给定边界品位的矿量和矿石平均品位。一、取样品位的统计分布规律为了确定取样品位的统计分布规律，首先将取样品位值绘成如图1-6所示的直方图。图中横轴为品位，竖轴为落入每一品位段的样 品数占样品总数的百分比。从直方图的轮廓线形状可以看出品位大体上属于何种分布；从直方图在横轴方向的分散程度可看出取样品位的变化程度。图1-6给出的是几种常见的品位分布情况。图(a)是一品位变化程度中等的正态分布，这样的分布在矿体厚大的层状或块状的硫化类矿床(如铜矿)中最为常见；图(b)是一品位变化小的正态分布，常见于铁、镁等矿床；图(c)是一对数正态分布(即品位的对数值服从正态分布)，品位变化大，此类分布常见于铝、锡、钨以及贵重金属(如金、铂)矿床；图(d)是一“双态”分布，即分布曲线是由两个不同分布组成的，说明样品来源于不同的样本空间。双态分布表明在矿床中很可能存在不同类型的矿石，或在不同区域呈现不同的成矿特征。如果图 (d)所示的情况出现，就需要对矿床地质和成矿机理进行深入分析，尽可能找出对应于不同分布的区域，然后对矿床进行区域划分，把来源于每一区域的样品进行分离，并做单独分析计算。(a)

(a)(d)(d)图1-6常见取样品位分布规律的直方图不同类型的矿床，其取样品位服从不同的统计学分布规律，但大多数矿床的品位服从正态分布或对数正态分布。下面对这两种分布的特征值及置信区间计算作简要介绍。二、正态分布TOC\o"1-5"\h\z检验样品值是否服从正态分布的一个简单方法是将样品的累加发生频率(即小于某一品位的样品数占样品总数的百分比)与品位绘在正态概率纸上(图1-7)。图中横坐标为累加概率，纵坐标为品位。如果 数据点基本落在一条直线上，那么就可以将样品的分布看成是正态分布。正态分布的特征值有均值P和方差。2。P和。2的真值是未知的。当我们获得n个样品，每个样品的值为x.(i=1,2, ,n)时，"和O2的估计值0和cf2可分别用下面的式子计算。 10=x=-龙X. (1-3)i=1Z(xiZ(xii=1-X)2S2=十(Zxi-nX2)(1-4)(1-5)S2的平方根S是样本空间均方差。的估计值。从统计学理论可知，一个正态样本空间的均值p的估计量0也服从正态分布，其均值为p，方差为S2=至 (1-6)0n

J J 1_I__L_1_I；IIIIII Il1 !1 110 1520253040 50 607075808590 9597 99 995累积发生概率（概率坐标,％）图1-7正态分布概率图设均值p小于pp的概率为p，大于pi的概率也为p，那么p落在p卢p侦之间的概率为1-2p（图1-8）。我们称［pp,pq］为均值在置信度为1-2#时的置信区间。当样品数n>25时，均值的68%和95%（即p为16%和2.5%）置信区间可用下面的式子计算图1-8置信区间示意图TOC\o"1-5"\h\z68%置信区间：W，口]=[X-^,X+二] （1-7）p1-p n7n95%置信区间：[口，口]=&-2二，x+2二] （1-8）p 1一p vn vn

当n<25时，计算任意置信度的置信区间的一般公式如下：(1-9)[口，口]=[x-t :—,x+t :—](1-9)p1-p 1-ptn i-p气「n式中、是学生分布(也称为t分布)表中自由度为n-1时，t<ti的概率为1-p的t值。t分布表见附表1-1。例如，当n=10，p=5%，即71-p=95%时，从表中可查得tip=1.833。如果样品的平均品位为X=20%,均方差S=10。那么置信度为"1-2p=90%的置信区间为[20—1.833^^,20+1.833^^]=[14.20,25.80]<10 V10也就是说，平均品位的真值p有90%的可能性是在14.2%和25.8%之间。三、对数正态分布当一个随机变量X的对数loge(X)服从正态分布时，X就服从对数正态分布。检验样品是否服从对数正态分布的方法与检验正态分布的方法相似。将图1-7中纵坐标由算术坐标变为对数坐标可得图 1-9。如果绘于图1-9的数据点基本落在一条直线上，就可认为样品是服从对数正态分布的。对数正态分布有二参数与三参数之分。当 log(X)是正态分布时，X服从二参数对数正态分布。在某些情况下， log(X)不是正态分布，而当X加上一常数6时，ioge(X+6)是正态分布，这时我们说X服从三参数对数正态分布。三参数对数正态分布有三个特征值：即加数6，(X+6)的对数均值和(X+6)的对数方差。当我们有n个样品时，就可以对这三个参数进行估值。如果样品数目足够大，6可用下式估计:(1-10)m也被称为几何均值或m(1-10)m也被称为几何均值或°=ff—2m式中，m为对应于50%累加概率的取样值,中直。f和f2分别为对应于累加概率p和1—p的取样值。理论上讲，p可以取任意值，但p取5%与20%之间时得到的结果最佳。

10累积发生频率(概率座标)I-4 L]J1II1 IIlli| |||5 10152025304050607075808590 9597累积发生概率(概率坐标,％)图1-9对数正态分布概率图令y.为(气+/3)的自然对数，即片.TogeE+S)(1-11)那么(x+6)的对数均值y用下式估计-1寸y=乙yn,i(x+6)的几何均值m的估计值为(1-12)m=eyr.、以.l /\、Iz-t-、r(1-13)对数万差a2的估计值为s2=1Y(y-7)2en ii=1或 S2=-Yy2—y2=1三参数对数正态分布的均值P，几何均值以下关系(1-14)(1-15)m与对数方差a2之间存在口=mea22—p(1-16)

当利用上面的公式从样品值计算出对数正态分布的特征值的估计值ni和S2以后，就可以获得均值p的估计值0:e0=my-p (1-17)式中，Y为从附表1-2中根据n和S；查得的系数，例如当n=10，S2=1.4时、Y为1.936。当n>1000时，Y可用下式计算：y=eSg'2 (1-18)置信度为1-2p的均值置信区间计算公式为：区间上限：01=(0+P)W1-P (1-19)区间下限：0p=(0+P)Wp-P (1-20)式中，当为从附表1-3a中根据n和S2查出的系数。附表1-3a列出的是当p=0.95时的W值，附表1-3b列出的是当p=0.05时的W值。更为完整的表可以在有关概率统计书中找到。当n很大(>1000)时，Wp可用下式计算：Wp=e(。t2+tp。t) (1-21)式中，a：=了(1+甘)，tp是从学生分布表(附表1-1)中查得的数[例1-1]:设从一金矿床取样10个，取样品位服从二参数正态分布，即6=0。应用式(1-12)至(1-14)计算得：对数均值：y=0.600,几何均值：m=1.822，对数方差：S2=0.050。试估计矿体的平均品位和90%置信区间。解：从附表1-2中查得：当S；=0.04和n=10时Y=1.020，当S；=0.06和n=10时，Y=1.030。因此，对于S：=0.05和n=10,线性插值得Y=1.025。应用公式(1-17)，算得平均品位的估计值为：0=1.822X1.025=1.868置信度为90%(即0.9)时，p=0.05,1-p=0.95。从表3a和3b中分别查得：0.950.050.950.05=1.194=0.897因此置信区间为:上限：0 =0W=1.868x1.194=2.2300.95 0.95下限：0 =0W=1.868x0.897=1.676也就是说：有90%的可能性，平均品位的真值p是在1.676和2.230之间。

第四节品位一矿量曲线边界品位是用于区分矿石与废石的临界品位值，矿床中高于边界品位的部分是矿石，低于边界品位的是废石。很显然，边界品位定的越高，矿石量也就越小。因此，边界品位是一个重要的参数，它的取值将通过矿石量及其空间分布影响矿山的生产规模、开采寿命和矿山开采规划。在一定的技术经济条件下就一给定矿床而言，存在着一个使整个矿山的总经济效益达到最大的最佳边界品位。边界品位的优化是当今世界矿业界的重大科研课题之一，但由于超出本节范围，这里不加详述。将一系列边界品位和与之相对应的矿石量绘成曲线就形成所谓的 品位一矿量曲线(图1-10)，由上面对边界品位的定义可知，品位一