对Surpac软件与地质统计学的矿床资源量估算结合应用的探究

■陈媛娟光有明

（1西宁特殊钢股份有限公司有限公司青海西宁 810005；2青海请齐鑫地质矿产勘查股份有限公司青海西宁 810016）

对Surpac软件与地质统计学的矿床资源量估算结合应用的探究

■陈媛娟1光有明2

（1西宁特殊钢股份有限公司有限公司青海西宁 810005；2青海请齐鑫地质矿产勘查股份有限公司青海西宁 810016）

利用Surpac软件将某铜矿床的钻探工程数据进行处理，建立其地质数据库，铜矿体的三维模型及品位变异函数模型，结合块体估值技术选择合理的地质统计学方法估算其资源量，并与距离幂次反比法估算结果进行对比研究，分析铜矿体不同方向上的品位变化趋势。结果表明，基于Surpac软件与地质统计学方法对该矿床进行资源量估算，结果可靠，操作简便，工作效率高，有利于推动地质统计学在国内矿产资源评价中的应用。

Surpac软件 地质统计学 矿床资源量 估算 结合

1　地质数据库

地质数据库是管理地质数据信息的数据库，其记录勘探工程的空间位置、样品化验数据、岩性等信息。强大的数据库管理功能是Surpac软件的一大特色，用户可根据需要对数据进行一系列处理，如查看、编辑、基础统计分析等。本文依据Surpac软件格式要求，建立了所研究矿床的地质数据库，结构见表1。

表1　地质数据库结构

2　样品组合与统计分析

区域化变量具有相同的承载是得到参数的无偏估计量的前提。样品组合即将空间不等长的样长和品位，量化到一些离散点。在地质数据库中提取全部样品长度，基础统计分析窗口。绝大多数样品长度为2.00m左右，故选择2.00m为组合样长度，采用长度加权法组合样品。

确定组合样长度后，在Surpac中以2.00m对Cu品位进行组合，并对组合后的Cu品位进行基础统计分析，得到其直方分布图。Cu品位近似服从正太分布，故选择普通克立格法对矿体品位进行插值。该矿床属于斑岩型铜矿床，Cu品位分布均匀，变化系数为0.58，特高品位的样品较少，无需处理。

3　块体模型与品位插值

3.1块体模型

Surpac软件将块体模型与地质统计学相结合，根据构建的变异函数模型对其他块体进行品位插值。对地质数据库中的Cu品位进行解译，以铜矿床工业指标一般要求圈定矿体（见图1）。以矿体模型为约束条件，在Surpac中创建块体模型，块尺寸为10m×10m× 5m，次级分块尺寸为5m×5m×2.5m。块体模型的功能与数据库相似，记载矿体的品位、岩矿类型、比重等诸多属性，可根据需要进行相关操作。

3.2品位插值

根据地质统计学原理，距离待估块体越近的已知点对待估块体的影响越大，即插值时所分配的权系数越高。在插值过程中，Surpac软件根据变异函数模型与各向异性椭球体参数分配已知点对待估块体的权系数。通过分析，本文选择普通克立格法对每个块体进行品位插值，并将估值结果与距离幂次反比法估值结果进行对比分析。通过分析可得，两种估算方法的Cu品位几乎相等，金属量的差异也非常微小，相对误差在1%内。这表明了普通克立格法估值的准确性，也验证了变异函数模型的合理性。

图1　矿体模型

4　校验及评估

任何一种估值都不可能与实际情况完全相符，运用普通克立格法对矿体品位进行估值亦然。因此，对估值结果进行校验及评估异常重要。其目的是再次检验所构建的变异函数模型是否合理，同时验证选择的估值方法是否得当。本次工作选择品位-吨位曲线和品位趋势分析对估值结果进行评估。品位-吨位曲线从Surpac软件中，可直接导出矿体的资源量报告，根据报告作出矿体的品位-吨位曲线（见图2）。从图2可以看出，普通克立格法与距离幂次反比法估值结果异常相近，仅在高品位区间（Cu品位大于0.60%）存在微弱差异。两种估值方法Cu品位在0～0.20%和0.60%～1.40%的金属量很少，大部分金属量品位落在Cu品位0.20%～0.60%之间。本例是一个斑岩型铜矿床，品位分布较均匀。因而，这也与矿体的地质意义相吻合。

图2　品味—吨位曲线

5　结语

本文分析并处理了某铜矿床的钻探工程数据，利用Surpac软件建立其三维地质模型（地质数据库、矿体模型等），对Cu品位进行了统计学分析，构建了其变异函数模型，选择普通克立格法对该矿床进行品位插值，估算其资源量，将估算结果与距离幂次反比法估算结果进行对比研究，从侧面证实了资源量估算工作的可靠性。同时，本文通过品位-吨位曲线揭示了两种资源量估算方法得出的结果相吻合，且符合矿床本身的地质意义；通过品位变化趋势分析揭示了两种方法估值后，矿体Cu品位在北、东两个方向的变化趋势保持高度一致且与组合样品位变化趋势保持较好的协调性，进一步校验了估算结果的可靠性。综上所述，基于Surpac软件与地质统计学方法对该矿床进行资源量估算，结果可靠，操作简便，工作效率高，有利于推动地质统计学在国内矿产资源评价中的应用。

P208[文献码]B

1000-405X（2016）-7-269-1