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

第一篇 金属矿床露天开采 第一章 品位与储量计算 第一节 概述 投资一个矿床开采项目，首先必须估算其品位和储量。一个矿床的矿量、品位及其空间分布是对矿床进行技术经济评价、可行性研究、矿山规划设计以及开采计划优化的基础，是矿山投资决策的重要依据。因此，品位估算、矿体圈定和储量计算是一项影响深远的工作，其质量直接影响到投资决策的正确性和矿山规划及开采计划的优劣。从一个市场经济条件下的矿业投资者的角度看，这一工作做不好可能导致两种对投资者不利的决策：(1)矿体圈定与品位、矿量估算结果比实际情况乐观，估计的矿床开采价值在较大程度上高于实际可能实现的最高价值，致使投资者投资于利润远低于期望值，甚至带来严重亏损的项目。(2)与第一种情况相反，矿床的矿量与品位的估算值在较大程度上低于实际值，使投资者错误地认为在现有技术经济条件下，矿床的开采 不能带来可以接受的最低利润，从而放弃了一个好的投资机会。然而，准确地估算出一个矿床的矿量、品位绝非易事。大部分矿体被深深地埋于地下，即使有露头，也只能提供靠近地表的局部信息。进行矿体圈定和矿量、品位估算的已知数据主要来源于极其有限的钻孔岩心取样。已知数据量相对于被估算的量往往是一比几十万乃至几百万的关系，即对一吨岩心进行取样化验的结果，可能要用来推算几十万乃至几百万吨的矿量及其品位。可以不过分地说，矿量、品位的估算是世界上最大胆的外推。因此，矿体圈定与矿量、品位估算不仅是一项十分重要的工作，而且是一项极具挑战性的工作。做好这一工作要求掌握现代理论知识与手段，并应用它们对有限的已知数据进行各种详细、深入的定量、定性分析；同时也要求从事这一工作的地质与采矿工程师具有科学的态度和求实精神。本章将较详细地介绍当今世界上常用的矿量、品位估算方法，包括探矿数据的分析、处理和用于品位估值的剖面法、平面法及矿床模型法等。地质统计学作为品位估值的一种方法，从其诞生起就显示了强大的生命力，得到了越来越广泛的应用，本章对此给予较大的篇幅。本章的主要目的不是教会读者如何一步一步地应用所介绍的方法，对一个矿床进行矿量、品位估算，而是使读者了解这些方法的内涵，为读者提供在不同条件下应用最合理的分析、评价方法所需的知识基础。第二节 探矿数据及其预处理 一、钻孔取样 用于矿体圈定与矿量、品位估算的数据主要来源于探矿钻孔的岩心 2000N 1900N 1800N 1700N 1600N 1500N 2000E 2100E 2200E 2300E 2400E 2500E 图 1-1 钻孔与勘探线示意图 zk2 zk1 表 11 钻孔岩芯信息记录 钻孔号：zk10；孔口坐标：6086.21E，6821.68N，205.01；设计深度：135M；实际深度：143.26M；开孔方位角：开孔倾角：90o；开孔日期：1994 年 10 月 12 日；终孔日期：1994 年 10 月 23 日 换 层 深 度 每层提取 每层岩芯 岩石矿石描述 自(M)至(M)共计(M)岩芯长度(M)采取率(%)0.00 13.93 13.93 表土层 13.93 30.69 16.76 1.6 9.5 云母石英岩：黄绿色，片状结构，主要组成矿物为石英(约 2530)，云母(约)和角闪石（约），其次有些磁铁矿。30.69 43.03 12.34 9.7 78.61 阳起磁铁石英岩：钢灰色灰白色，细粒结构，主要组成 矿 物 为 石 英(约 )，磁铁石(约)，阳起石(约)。：:：取样。钻孔一般按照一定的网度布置在一些叫做勘探线的直线上（图1-1）。在钻孔过程中，每钻一定深度（一般在 3 米左右）将岩心取出，做好标记后按顺序放在箱中供搬运、贮存和化验。地质人员对取出的岩心进行定性观察和简单的测试，以确定每一段岩心的主要物理特性，如岩心长度、岩性、颜色、硬度等，并记录下来，形成对钻孔穿过地段的地质特性的定性描述。表 1-1 是一个钻孔的岩心观测结果的部分记录表。为直观起见，常常把表中的数据和文字描述绘成钻孔柱状图（图1-2）。为了确定岩心的化学成分和品位，将岩心的一半送往化验室进行化验，另一半保存下来备用。样品的化验结果记录在如表 1-2 所示的表中，或输入计算机的数据库中。手工记录时常将表 1-1 和表 1-2 合并为一个表，称为钻孔地质资料记录表。对所有钻孔的定性描述和取样化验结果构成了勘探区域的基本地质数据，这些取样化验数据是进行矿体圈定和矿量、品位估算的依据。在矿量和品位计算前，一般需要对取样数据进行预处理，包括样品组合处理和“极值”样品的处理。0 4 8 11 风化砂岩 页岩 浅灰色 石灰岩 灰色 图 1-2 钻孔柱状图 表 1-2 钻 孔 岩 芯 取 样 化 验 结 果 记 录 钻孔号：zk10；孔口坐标：6086.21E，6821.68N，205.01；设计深度：135M；实际深度：143.26M；开孔方位角：开孔倾角：90o；开孔日期：1994 年 10 月 12 日；终孔日期：1994 年 10 月 23 日 试样号 采 样 间 隔 化学分析结果(%)备 注 自(M)至(M)共计(M)TFe Fe SFe 1083 30.69 33.69 3.00 29.80 16.60 22.50 1084 33.69 36.69 3.00 32.20 15.60 25.10 1085 36.69 39.69 3.00 32.95 16.00 28.00 1086 39.69 43.03 3.34 26.40 14.00 21.00 ：:：:二、样品组合处理 样品组合处理就是将几个相邻样品组合成为一个组合样品，并求出组合样品的品位。当矿岩界限分明，且在矿石段内垂直方向上品位变化不大时，常常将矿石段内（即上下矿岩界限之间）的样品组合成一个组合样品（图 133），这种组合称为矿段组合。组合样品的品位g是组合段内各样品品位的加权平均值，即 gl gliiiniin11 (1-1)式中，li为第i个样品的长度；gi为第i个样品的品位；n为矿石段内样品个数。式 1-1 中用的是长度加权，是最常用的方法。如果不同样品的比重相差较大，可以采用重量加权法。岩石 岩石 l1 l2 l3 li ln g1 g2 g3 gi gn l5,l1,l5 l4 l3 l2 l1 g5,g1,g5 g4 g3 g2 g1 188m 200m H(台阶高度)矿石 图 1-4 台阶样品组合示意图 图 1-3 矿段样品组合示意图 对于拟用露天开采的矿床，更具实际意义的样品组合处理是台阶样品组合，即把一个台阶高度内的样品组合成一个组合样品（图 1-4）。组合样品的品位为：gl gHiiin1 (1-2)式中，H为台阶高度。当一个样品跨越台阶分界线时（如图 1-4 中第一和第五个样品），在计算中样品的长度取落于本台阶的那部分长度（即图 1-4 中的l1和 l5），样品的品位不变。对钻孔取样进行台阶样品组合处理的意义在于：(1)对取样数据进行统计学、地质统计学分析，以及利用取样值进行品位估值时，只有当每个样品具有相同的支持体，即每个样品的体积相同时，分析计算结果才有意义。(2)露天开采在垂直方向上是以台阶为开采单元的，一旦台阶的参考标高和台阶高度被确定，沿台阶高度无论品位如何变化，也无法进行选别开采。因此，在一个台阶高度内采用不同的取样品位是毫无意义的。(3)组合样品的品位较原样品品位变化小，在一定程度上减轻了“极值”品位对分析计算的影响，也使样品的统计分布曲线和半变异函数曲线（这些概念将在以后几节讲述）趋于规则。(4)样品组合处理减少了样品总数，节省计算机内存和计算时间。三、极值样品(Outlier)处理 极值样品是指那些品位值比绝大多数样品的品位（或样品平均品位）高出许多的样品，它们在贵重金属矿床较为常见。例如，在一金矿床取样 1000 个，经化验，这些样品的平均品位为 10 克/吨，其中有 十个样品的品位在 100 克/吨以上，这十个样品就可以被看成是极值样品。究竟品位比均值高出多少的样品算是极值样品，没有统一的、现成的标准，需视具体情况而定。极值样品虽然数量少，但对金属量影响大，为使品位的分析计算结果不致过分乐观，人们常常在实践中采用以下处理方法：(1)限值处理：即将极值样品的品位降至某一上限值。比如在上述例子中，将所有高于 100 克/吨的样品的品位降至 100 克/吨。(2)删除处理：即将极值样品从样本空间中删去，不参与分析计算。使用上述处理方法时应非常谨慎。虽然极值样品在数量上占样品总数的比例很小，但由于其品位很高，对矿石的总体品位和金属量的贡献值都很大。因此，不加分析地进行降值或删除处理会严重歪曲矿床的实际品位和金属含量，人为地降低矿床的开采价值。这一点可用下面的例子说明。假设对一金矿床进行钻探取样后得知，品位值服从对数正态分布（图1-5）。所有样品的平均品位为g=10 克/吨，中值为m=3 克/吨（即高于 3 克/吨和低于 3 克/吨的样品各占 50）；有 1的样品品位高于 100克/吨。若将这 1的极值样品取出，单独计算其平均值，得 190 克/吨。那 么 这1 的 样 品 对 矿 床 总 金 属 量 的 贡 献 为（190 1 ）/g=1.9/10=19%。也就是说，百分之一的数据量代表的是百分之十九的金属量！假如取边界品位为 3 克/吨（高于 3 克/吨为矿石，否则为废石），矿石的平均品位（即高于 3 克/吨的那部分样品的平均品位）经计算为16 克/吨。如果把极值样品从样品空间删除，矿石的平均品位变为（16501901）/（501）12.45 克/吨，也就是说，矿石品位被低估了 22。如果将极值样品进行限值处理，将其品位值降到 100克/吨，矿石的平均品位变为(12.4549%+1001%)/50%=14.2 克/吨，也就是说，将矿石品位低估了 11。在正常、稳定的经济环境中，采矿的利润率也就是 15左右。因此，不加分析地将极值样品进行删除或限值处理，很可能将本来能够获取正常利润的矿床人为地变为没有开采价值，从而导致错误的投资决策。这对于一个在市场经济条件下，以盈利为主要目的的矿业投资者来说，无疑是一个重大的决策失误。这里必须澄清的是，极值样品是实实在在存在的有效样品，并不是指那些由于化验或数据录入错误造成的、具有“错误品位值”的样品。如果有根据认为某些样品的品位是错误的，将这些样品从样本空间中删除不仅是合理的而且是必要的。对极值样品的最理想的处理方法是，经过对探矿区域的地质构造和成矿机理进行深入分析，将这些样品的发生区域（或构造）划分出来，在进行品位与矿量的分析计算时，这些样品只参与其发生区域的品位与矿量计算，而不把它们外推到发生区域之外。但是在大多数情况下，由于钻孔网度大，已知的地质信息满足不了这种区域划分的要求。这时，可以将矿床看成是由两种不同的矿化作用形成的：样品中占绝大多数的“正常样品”可以看作是由主体矿化作用产生的样本空间；极值样品是由次矿化作用产生的样本空间。然后利用统计学方法计算出空间任一点属于每一类矿化作用的概率，再根据这些概率计算矿床的品位与矿量。这一方法超出了本书的范畴，有兴趣的读者可参阅 Journel(1988)和Parker 等人(1979)的论文。g/t m=3 g=10 100 1.0%图 1-5 金矿取样品位对数正态分布示意图 第三节 取样数据的统计学分析 对取样数据进行上述的预处理以后，做一些统计学分析可以提供不少有关矿床的有用信息。因此统计学分析常常是取样数据分析的第一步。对数据进行统计学分析的主要目的是确定：(1)品位的统计分布规律及其特征值；(2)品位变化程度；(3)样品是否属于不同的样本空间；(4)根据样品的分布特征，初步估计矿床的平均品位以及对于给定边界 品位的矿量和矿石平均品位。一、取样品位的统计分布规律 为了确定取样品位的统计分布规律，首先将取样品位值绘成如图 1-6所示的直方图。图中横轴为品位，竖轴为落入每一品位段的样 品数占样品总数的百分比。从直方图的轮廓线形状可以看出品位大体上属于何种分布；从直方图在横轴方向的分散程度可看出取样品位的变化程度。图 1-6 给出的是几种常见的品位分布情况。图(a)是一品位变化程度中等的正态分布，这样的分布在矿体厚大的层状或块状的硫化类矿床（如铜矿）中最为常见；图(b)是一品位变化小的正态分布，常见于铁、镁等矿床；图(c)是一对数正态分布（即品位的对数值服从正态分布），品位变化大，此类分布常见于钼、锡、钨以及贵重金属（如金、铂）矿床；图(d)是一“双态”分布，即分布曲线是由两个不同分布组成的，说明样品来源于不同的样本空间。双态分布表明在矿床中很可能存在不同类型的矿石，或在不同区域呈现不同的成矿特征。如果图(d)所示的情况出现，就需要对矿床地质和成矿机理进行深入分析，尽可能找出对应于不同分布的区域，然后对矿床进行区域划分，把来源于每一区域的样品进行分离，并做单独分析计算。30%15%0 50%25%0 2.0 1.0%Cu(a)50 25%Fe(b)不同类型的矿床，其取样品位服从不同的统计学分布规律，但大多数矿床的品位服从正态分布或对数正态分布。下面对这两种分布的特征值及置信区间计算作简要介绍。二、正态分布 检验样品值是否服从正态分布的一个简单方法是将样品的累加发生频率（即小于某一品位的样品数占样品总数的百分比）与品位绘在正态概率纸上（图 1-7）。图中横坐标为累加概率，纵坐标为品位。如果 数据点基本落在一条直线上，那么就可以将样品的分布看成是正态分布。正态分布的特征值有均值和方差2。和2的真值是未知的。当我们获得n个样品，每个样品的值为xi(i=1,2,n)时，和2的估计值和2可分别用下面的式子计算。xnxiin11(1-3)()222111Snxxiin(1-4)或 Snxnxi22211()(1-5)S2的平方根S是样本空间均方差的估计值。从统计学理论可知，一个正态样本空间的均值的估计量也服从正态分布，其均值为，方差为 Ssn22(1-6)(c)6 30%30%15%15%0 0 3 g/t Au(d)图 1-6 常见取样品位分布规律的直方图 图 1-7 正态分布概率图 设均值小于p的概率为p，大于1-p的概率也为p，那么落在 p与1-p之间的概率为 1-2p（图 1-8）。我们称p，1-p为均值在置信度为 1-2p时的置信区间。当样品数n25 时，均值的 68和 95（即p为 16和 2.5）置信区间可用下面的式子计算 68置信区间：,ppxsnxsn1 (1-7)95 置信区间：,ppxsnxsn122 (1-8)1-2p p p 1-p p 图 1-8 置信区间示意图 当n25 时，计算任意置信度的置信区间的一般公式如下：,ppppx tsnx tsn111 (1-9)式中t1-p是学生分布（也称为t分布）表中自由度为n-1时，t1000 时，可用下式计算：eSe22(1-18)置信度为 1-2p的均值置信区间计算公式为：区间上限：11pp()(1-19)区间下限：pp()(1-20)式中，p为从附表 1-3a 中根据n和Se2查出的系数。附表 1-3a 列出的是当p=0.95 时的值，附表 1-3b 列出的是当p=0.05 时的值。更为完整的表可以在有关概率统计书中找到。当n很大(1000)时，p可用下式计算：ptetp t()22(1-21)式中，tSnSee22221()，tp是从学生分布表（附表 1-1）中查得的数值。例 1-1：设从一金矿床取样 10 个，取样品位服从二参数正态分布，即=0。应用式(1-12)至(1-14)计算得：对数均值：y=0.600，几何均值：m=1.822，对数方差：Se2=0.050。试估计矿体的平均品位和 90置信区间。解：从附表 1-2 中查得：当Se2=0.04 和n=10 时=1.020，当Se2=0.06和n=10 时，=1.030。因此，对于Se2=0.05 和n=10，线性插值得=1.025。应用公式(1-17)，算得平均品位的估计值为：1.8221.025=1.868 置信度为 90（即 0.9）时，p=0.05，1-p=0.95。从表 3a 和 3b 中分别查得：0.95=1.194 0.05=0.897 因此置信区间为：上限：.0 950 95186811942230 下限：.0 050 05186808971676 也就是说：有 90的可能性，平均品位的真值是在 1.676 和 2.230之间。第四节 品位矿量曲线 边界品位是用于区分矿石与废石的临界品位值，矿床中高于边界品位的部分是矿石，低于边界品位的是废石。很显然，边界品位定的越高，矿石量也就越小。因此，边界品位是一个重要的参数，它的取值将通过矿石量及其空间分布影响矿山的生产规模、开采寿命和矿山开采规划。在一定的技术经济条件下就一给定矿床而言，存在着一个使整个矿山的总经济效益达到最大的最佳边界品位。边界品位的优化是当今世界矿业界的重大科研课题之一，但由于超出本节范围，这里不加详述。将一系列边界品位和与之相对应的矿石量绘成曲线就形成所谓的品位矿量曲线（图 1-10），由上面对边界品位的定义可知，品位矿量曲线是一条递减曲线。由于品位矿量曲线指明了任一给定边界品位条件下的矿石量，它是对矿床进行初步技术经济评价的重要依据。当品位服从均值为和方差为2的正态分布时，品位矿量曲线上的每一点可由下式求得：Tc=T(uc)(1-22)式中，T为总矿岩量，即边界品位为零的矿量，对于一给定矿床或矿床中的一给定区域，T是已知的。(uc)是高斯函数，即标准正态分布从uc到的积分：()uedtctuc1222(1-23)式中，uc是边界品位gc的标准正态变量，即：ugcc(1-24)Tc中含有的金属量Qc可由下式计算：50 40 30 20 10 3 2 1 0 边界品位（g/t）矿 量(Mt)图 1-10 品位-矿量曲线示意图 QTuTeccuc()222 (1-25)对应于边界品位gc的矿石平均品位，即品位高于边界品位gc的那部分物料的平均品位为：GQTccc(1-26)当品位服从三参数对数正态分布时，可用下面的公式计算品位矿量曲线：TT ucc()1(1-27)式中，和T与正态分布条件下的定义相同；uc1为：ugceece112log(1-28)Tc中的金属含量Qc为：QTuuucccc()()()221(1-29)式中，uc2由下式算得：uucce21(1-30)应用品位矿量曲线进行品位、矿量分析时，必须注意以下几点：(1)品位分布是从样品值的分布得出的，分布的特征值，或e是未知的，计算中只能用它们的估计值，S或Se。(2)露天开采时，矿石不是以几公斤大小的样品为单位采出的。对于选定的开采设备（电铲）和台阶高度，存在所谓的最小选别单元（SMU），矿床中的矿石是以 SMU 为单元采出的。SMU 在体积上要比样品大得多，如果把整个矿床分为体积为 SMU 的小块（称为单元体），那么这些小块的品位分布较样品品位分布更为集中（即方差更小）。因此根据样品分布计算的品位矿量曲线并不能用来预报将被采出的品位矿量关系。(3)单元体的真实品位是未知的，单元体是否是矿石，不是根据其真实品位确定的，而是根据对单元体的品位的估计值确定的。由于估计有误差，根据估计值得出的品位矿量曲线与实际采出的品位矿量关系有一定的差别。