韶关红岭钨矿云英岩型钨矿储量估算方法探讨

付裕

摘要：我队于2014年至2019年在红岭矿区先后开展红岭钨矿接替资源勘查、红岭矿区614线-626线钨矿详查工作，探获的三氧化钨（WO3）达大型矿床规模。准确估算资源储量非常重要，是今后企业对矿山进行矿权评估、开发利用的基础数据。矿产资源储量估算方法主要有断面法、水平投影法、平行剖面法、SD法等，本次根据红岭矿区情况选择水平投影法和平行剖面法。

关键词：红岭矿区；资源储量估算；地质块段法

1.引言

我国矿产资源储量估算采用断面法、地质块段法等传统方法，开发勘探矿产资源过程中，根据不同情况采取相应的储量估算。任务要求不同，储量估算的程度也不相同。矿产资源估算工作占据着矿产勘查任务重要的一个环节。对于当下多金属矿床的开发建设，及时顺应市场，研究更为高效的估算方法，成为了经济建设不可缺少的基础元素。以传统几何算法进行的资源量估算方法（水平投影法、平行剖面法、等值线法、算术平均法等）以及数学地质方法（主要有地质统计学法、SD法和距离平方反比法等）（阳正熙，严冰，2009）。随着矿区勘查程度的提升，选择资源储量估算的方法导致资源储量的结果不同，选择合适的储量估算方法是重要的。

2.地质背景

2.1区域地质背景

区域上处于华南加里东隆起粤北海西-湘南-印支凹陷区内，贵东—大东山—九连山东西向构造岩浆带中部，热水SN向断裂构造与NE向断裂构造复合部位。区域地层分布简单，主要为寒武系八村群、泥盆系中统桂头群、杨溪组、老虎头组、春湾组和东岗岭组，上统天子岭组及帽子峰组。

区内分布水源山背斜，轴线呈NNE向，长约40km，由下古生界和泥盆纪地层组成，岩层倾角40°～70°，东翼较西翼陡；轴部于官渡镇一带向南西方向倾伏，形成一个向南西倾伏的背斜构造。背斜西侧为一轴线呈NNW向的对称正常的翁城向斜；东侧则为轴向NE、轴面倾向NW的翁源倒转向斜。区域上的断裂构造发育，表现为一系列近NNE-NNW向压扭性断裂穿切NNE-NE向压扭性断裂呈直交反接，在交叉部位控制红岭石英脉型钨矿成矿（与云英岩型钨矿关系不明显），区内广泛分布燕山期花岗岩，北部属贵东岩体，南部统称热水岩体，呈岩基状产出。据区测资料贵东岩体侵入的最新地层是中下泥盆统桂头群，覆于其上的最老地层是下第三系丹霞群。热水岩体侵入于下古生界和中下泥盆统桂头群中。热水岩体的锆石U-Pb同位素年龄约160Ma。所在区域经历了燕山岩浆旋回第三、第四期较大的岩浆活动，主要有燕山第三期花岗岩类、燕山第四期花岗岩。

2.2矿区地质概况

（1）地层。矿区地层仅出露第四系（Q）：腐植层及砂粘土层，面积约0.1km2，主要分布在矿区的茶树坳、梅子坑一带，层厚3m～20m。

（2）岩浆岩。矿区位于热水岩体中部，燕山期花岗岩类沿区域性SN向断裂构造多次侵入，形成多阶段、多岩相的复式花岗岩体。按岩性可分为黑云母花岗岩和白云母花岗岩。前者分布面积最广，约占矿区面积80%；后者分布于矿区的中南部，分布面积约占矿区面积20%。岩体划分为三个阶段，属燕山三期产物。红岭钨矿各阶段的岩体进行了锆石U-Pb定年，不同地点、不同岩性的花岗岩锆石U-Pb年龄数值在误差范围内基本一致，都在160Ma左右。

（3）构造。红岭钨矿区位于区域性NE向压扭性断裂与近SN向压扭性断裂交叉部位。区内主要的导矿构造为一组近SN（NNE—NNW）走向的压扭性冲断裂，倾向W，倾角70°～85°。

2.3矿体特征

矿区主要有4个云英岩型矿体，其中主要以Ⅰ、Ⅱ号为主矿体，Ⅰ、Ⅱ号矿体水平投影面积占矿区的70%以上。根据矿体赋存特征，矿体整体上具有向北侧伏的特征，矿体标高逐渐下降，根据矿化特征、矿体形态起伏规律及空间分布特征等圈定4个矿体，基本保障矿体厚度的相对稳定性，以及矿体形态的相似性。东西方向呈似层状或透鏡状赋存于蚀变花岗岩体中，空间上平行重叠，总体形态呈斗笠状分布于近SN向的波状起伏的岩脊中且向四周倾斜，倾角5°～30°。矿体多集中在616～624线的矿化蚀变中心。矿体分布于细粒白云母花岗岩顶面突起部位，单个矿体长865m～1577m，宽746m～779m，平均厚17.63m～19.52m，最大厚度达151.29m。矿体厚度沿走向和倾向均体现为中间大两端小纺锤状，矿石中WO3的含量由上至下逐渐降低，在标高150m左右，钨矿化显著减弱。矿体分布与特征和东西展布类似总体形态呈斗笠状，向南北侧伏。

3.矿区资源量估算

3.1工业指标的选定

按《钨、锡、汞、锑矿产地质勘查规范》（DZT0201-2002），结合接替资源勘查报告确定的工业指标，确定本矿区云英岩型钨矿床地质勘查一般工业指标为：矿石最低工业质量分数（WO3）为0.12%；边界质量分数（WO3）为0.064%；最小可采厚度为1m；夹石剔除厚度为5m，考虑本区云英岩型白钨矿体延伸长、规模大、面状矿化等特点，结合规范要求，夹石剔除厚度取规范值（2m～5m）的上限值。

3.2资源储量估算的方法、原则

（1）水平投影法。根据矿区内进行资源量估算的矿体形态多呈似层状或透镜状、走向、倾向延伸大、倾角平缓，后期破坏矿体的断层不发育，采用水平投影地质块段法计算，按规范划分块段，进行资源量估算。计算公式如下：P（金属量）=S×H×D×C矿体投影面积（S）：块段矿体水平投影图面积采用Mapgis软件计算，单位为m2，小数点后舍去取整数。

矿体铅直厚度（H）：单工程矿体厚度计算，单工程矿体厚度由单工程内对应矿体厚度之和求得；块段平均铅直厚度计算由于云英岩型钨矿主矿体（Ⅰ号、Ⅱ号）的厚度变化系数接近80，相对较为稳定，且工程分布较均匀，因此采用单工程矿体厚度算术平均法求得块段平均厚度（H）。厚度是矿体的铅直厚度。

矿石小体重（D）矿体平均质量分数（C）

（2）平行剖面法。水平投影地质块段法计算，计算公式如下：P=S×H×D矿体勘探线剖面面积（S）：对于采用平行剖面法进行资源量估算的矿块，采用矿体勘探线剖面面积。利用矿山及本次工作的实测资料，用mapgis绘图软件，绘制成1：2000勘探线剖面，再量取并换算出保有块段的剖面面积。

勘探工程间距（H）：由于矿体形状较简单，呈脉状产出，勘探线间距为100m，如两勘探线或工程之間一个见矿，另一个未见矿，取用见矿工程向外推50m，作为矿体边界。

D—矿石小体重（t/m3）

3.3资源储量结果

为了验证水平投影地质块段法计算资源储量的可靠性，现采用水平投影法与平行剖面法进行对比论证，现以矿区332（工业Ⅱ-3、工业Ⅱ-5、）块段为例，进行论证，计算其差异率，不同计算方法计算资源储量差异较小说明该方法准确可行，其差异率均小于10%。计算结果如下：

4.讨论

通过对不同方法对矿区资源储量进行估算，并进行差异化分析对比，对矿区资源储量圈定的逐一筛选，由于矿体形态较为简单，最后发现了平行剖面法与水平投影法估算资源储量相当。不同的资源储量估算方法各自的特点，应结合矿床的地质特征进行选择。

参考文献：

[1]叶天竺.矿床模型综合地质信息预测技术方法理论框架[J].吉林大学学报（地球科学版）， 2013（04）： 1053-1072.

[2]赵鹏大主编.矿产勘查理论与方法[M]武汉：中国地质大学出版社， 2001.

[3]叶天竺，吕志成，庞振山等编著.勘查区找矿预测理论与方法（总论）[M]北京：地质出版社， 2014.

[4]高山，张本仁，韩吟文.广东翁源红岭钨矿床成矿作用与花岗岩浆演化之间的地球化学关系探讨[J].地球科学， 1986（01）： 40.

[5]陈国旭，等编著.传统资源储量估算信息化研究现状及发展方向，地质与测量，社， 2014.

[6]武亚峰，等编著.青海省深部金属矿产资源储量分类与计算方法研究，世界有色金属， 2019.