近年来,随着地质找矿工作的持续推进,地表及浅部矿产资源不断枯竭,“攻深探盲”和“探边摸底”已成为我国今后找矿的主攻方向(
赵鹏大,2007),寻找深部矿产资源的需求日趋迫切。矿床原生晕是将岩(矿)石作为主要研究对象,因此,原生晕找矿方法是目前深部找矿的最直观且高效的方法之一,特别是在金矿地质找矿和危机矿山深边部找矿预测中取得了良好的找矿效果(
Guo et al.,2022;
刘鹤群等,2024)。国内外众多学者利用原生晕找矿方法在金矿深部预测中取得了丰硕的成果(
贺昌坤等,2020;
Helba et al.,2020;
Aliyari et al.,2021;
Wu et al.,2022;
李彦强等,2022),充分说明原生晕找矿方法在深部矿产资源预测和隐伏矿体勘探过程中发挥着重要的作用,对于深边部找矿工作具有重要的理论和实践指导意义。
青海省都兰县达热尔金矿床位于东昆仑东段昆中断裂北侧的沟里矿集区(
图1),该区分布有果洛龙洼(金资源量超过40 t)、坑得弄舍(金资源量超过39 t)、瓦勒尕(金资源量超过12 t)、按纳格(金资源量超过8 t)、阿斯哈(金资源量超过6 t)和德龙(金资源量超过5 t)等多处大中型金矿床(
陈加杰,2018),金矿找矿潜力巨大。许多学者对该区域金矿床地质特征、岩石地球化学特征、矿床成因及找矿预测等进行了大量研究(
岳维好等,2013;
陈加杰等,2016;
李治华,2023;
王凤林,2023;
王晓云等,2024;
李兴辉等,2025),认为沟里地区金矿床属于与印支晚期造山运动有关的造山型金矿,具有明显的构造控矿特征,指出果洛龙洼、瓦勒尕等矿区的主矿脉深部以及沿走向方向均具有较大的找矿潜力,达热尔、迈龙和色日等矿区也具有良好的找矿前景。达热尔金矿是青海省有色第二地质勘查院新发现的一处小型金矿床(金资源量约为2.29 t,金品位为13.96×10
-9),对该矿床1/2.5万水系沉积物综合异常HS4和HS5进行查证过程中共发现6条矿体,实现了地表找矿突破。但由于该矿区目前尚处于预查阶段,对成矿地质背景和矿床特征的研究较为薄弱,尚未开展岩石地球化学、矿床原生晕和成矿预测等方面的系统研究工作。
为了摸清深部赋矿规律,根据金矿原生叠加晕理论,对达热尔金矿床钻孔样品进行原生晕元素测试,并运用聚类分析、轴向分带序列和矿体剥蚀系数等数学统计方法,构建原生晕元素分带模型,深入分析元素组合特征及变化规律,并利用可控源音频大地电磁法(CSAMT)对深部找矿前景进行预测,从而为沟里地区深部找矿提供更精准的理论依据。
1 成矿地质背景
东昆仑成矿带呈近EW向展布于东昆仑北坡—鄂拉山一带,北部以昆北断裂为界,与柴达木地块分开,南部以昆南断裂为界,与北巴颜喀拉成矿带为邻[
图1(a)]。该造山带主要经历了加里东—印支期复杂的演化过程。带内显生宙地层发育齐全,岩浆活动频繁而强烈,形成巨大的花岗岩带。晚奥陶世、晚石炭—早二叠世和晚三叠世3个期次的构造活动强烈,断裂构造复杂,主干构造由昆北、昆中和昆南3条断裂带构成(
岳维好等,2013),呈近EW向展布,构成带内的构造格架。
沟里金矿田位于东昆仑造山带东段,区内出露地层由老到新主要有元古宇、下古生界、上古生界、中生界和新生界。元古宇包括金水口岩群、小庙岩组和万保沟群,岩性主要为中高级变质岩,共同组成区内基底岩系,也是区内部分金银矿床的赋矿围岩。下古生界包括纳赤台混杂变质岩和祁漫塔格群浅变质岩,该地层是区内金矿床的重要赋矿围岩。上古生界包括哈拉郭勒组、大干沟组、缔敖苏组、浩特洛洼组、树维门科组、切吉组、马尔争组和格曲组,岩性以沉积碎屑岩为主,其中哈拉郭勒组和浩特洛洼组为区内重要的金矿赋矿围岩(
陈加杰,2018)。中生界包括隆务河组、洪水川组、闹仓坚沟组、希里可特组、八宝山组、鄂拉山组和羊曲组,其中洪水川组和鄂拉山组为区内重要的银多金属矿赋矿地层。新生界包括沱沱河组、干柴沟组、油砂山组、临夏组、羌塘组和第四系。
区内断裂十分发育,尤以昆中断裂带规模最大,贯穿全区[
图1(b)],是区内主干断裂,断裂性质主要表现为压性或压扭性,控制着区内地层和岩浆岩的展布,具有多期活动的特点(
李兴辉等,2025),并发育有一系列NW、NWW-近EW向的次级断裂。近EW向构造规模大且延伸长,具有多期活动的特点,是区内主要的导矿和控矿构造,NW向构造为区内主要的容矿构造。
区域内岩浆活动强烈,尤其是古生代—中生代与原特提斯洋和古特提斯洋扩张、俯冲、闭合碰撞和后碰撞事件相关的岩浆岩最为发育,在全区均有分布。区内岩浆活动强度呈波动变化特点,其中以奥陶—泥盆纪(加里东期)和二叠—三叠纪(华力西期—印支期)岩浆活动强度最大,岩浆岩早期主要发育花岗闪长岩和闪长岩类,晚期主要发育二长花岗岩和钾长花岗岩等,代表了2期不同的造山旋回。
2 矿床地质特征
达热尔金矿床出露地层简单,主要为古元古界金水口岩群中深成变质岩系,分布于矿区中东部(
图2),其他地段为小规模零散分布,地层走向总体呈NE向,往东逐渐转为NW向,倾角为50°~70°,主要岩性有黑云斜长片麻岩、黑云石英片岩、斜长角闪片岩和大理岩。区内断裂发育,主要有近EW-NW向、NNE-近SN向和NE向3组,其中,NW向(F
2)和NNE-近SN向(F
4、F
5、F
11、F
14、F
16和F
17等)断裂是主要含银和含金构造。矿区岩浆活动强烈,火山岩不发育,以华力西期和加里东期侵入岩为主,岩性以灰白色中—粗粒花岗闪长岩、石英二长闪长岩和英云闪长岩为主,与金水口岩群呈断层或侵入接触,局部具混合岩化现象。
青海省有色第二地质勘查院于2016—2020年在达热尔金矿开展了预查工作,在主矿区圈定出9条含金矿带(AuⅠ、AuⅡ、AuⅢ、AuⅣ、AuⅤ、AuⅥ、AuⅦ、AuⅧ和AuⅨ),其中AuⅡ和AuⅢ为主矿带,1条含银矿带(AgⅠ),21条金矿体、金银矿体和银矿体。矿体均产于构造蚀变带内,除了AgⅠ矿体外,其余矿体走向均为NNE向,产状近直立,呈条带状、细脉状和透镜状产出,形态复杂、大小不一(
图3)。达热尔金矿床的矿体特征与南部的迈龙、色日金矿床矿体特征相似,具有尖灭再现的特征,随着矿体尖灭品位也趋于矿化,而在断层交会处或拐弯处矿体膨大,金矿石品位较高。
矿体围岩主要有花岗闪长岩、二长花岗岩和黑云斜长片麻岩,近矿围岩以蚀变闪长岩为主,地表岩石普遍具碎裂—片理化。围岩蚀变主要为硅化、高岭土化、绿泥石化和绢云母化等,靠近矿体的围岩蚀变较强,围岩与矿体界限较清晰,蚀变现象在空间上从蚀变带向两侧逐渐变弱,呈渐变过渡关系。
矿石类型主要有石英脉型和构造角砾岩型,矿石矿物主要有黄铁矿、黄铜矿、方铅矿、毒砂和次生褐铁矿,脉石矿物主要有石英、长石和黑云母等。矿石构造主要有胶状、块状、细脉—网脉状和浸染状等,矿石结构有半自形—他形粒(柱)状、交代、粒状变晶和碎裂结构等。
3 样品采集与分析
本次研究选择矿区内控制AuⅡ和AuⅢ主矿带的63号勘探线,共采集钻孔(ZK6301和ZK6302)岩芯地球化学测试样品209件。采样方法为全孔连续捡块法,一般点距为3~5 m,在矿体或蚀变较强的地段采样间距为3 m,要求采集新鲜岩(矿)石样品,并以多点采样组合成一个样品,避免随意采样,单个样品重量一般不低于500 g。本研究对Au、As、Sb、Cu、Pb、Zn、Ag、W、Sn、Mo、Hg、F、Mn、Bi和Co共15种元素进行分析。
为确保样品测试分析的准确性,针对不同元素采用不同的分析方法。其中,Au、W和Mo元素采用电感耦合等离子体质谱仪进行分析;Cu、Pb、Zn、Co、F和Mn元素采用X射线荧光光谱仪进行分析;As、Sb、Bi和Hg元素采用原子荧光仪进行分析;Ag和Sn元素采用一米平面光栅摄谱仪和测微光度计进行分析。样品测试由青海省有色第一地质矿产勘查院测试中心完成,所有元素的平均报出率均大于90%,样品内检合格率均大于95%,按样品总数5%的比例抽取样品作为密码样,交由分析人员单独进行重复分析,为提高测试的准确性,质量监控采用国家一级标准物质,仪器检出限、样品测试报出率以及分析方法的准确度和精密度等各项指标均符合相关技术规范要求。
4 结果与讨论
4.1 原生晕地球化学特征
(1)异常下限的确定
叶天竺等(2014)研究指出,矿体品位的高低与矿床原生晕异常含量之间具有一定的函数关联,即:当某个成矿元素的原生晕异常含量超过其矿石边界品位的1/10时,有可能在异常区域内找到相关矿化,因此通常情况下将这一浓度定义为各元素异常的内带浓度/内带含量(异常内带),而中带为内带含量的1/2,外带为中带含量的1/2。依此原则,将矿床的相关元素含量确定为其矿种边界品位的1/10,中带为内带的1/2,外带为中带的1/2。将矿床相关伴生成分的基础信息经过对数变换分析,并剔除离群信息,直至相关数据达到正态分布特征,而后按背景值加2倍的标准离差选取其下限值,再按背景值1倍、2倍和4倍分别圈定异常的外带、中带和内带。另外,由于每个钻孔的岩芯数据含量可能会出现高、低不稳定,所以对圈定异常的基本数据先进行三点均匀化处理,即根据浓度分级指标最大限度地圈出完整异常(
刘玉军等,2018)。按此原则圈定出达热尔金矿床各元素浓度分带参数,具体见
表1。
根据分析结果确定的异常含量分带值,圈定出达热尔金矿床AuⅡ号矿体63号勘探线原生晕剖面等值线图(
图4)。由
图4可知,大部分元素异常分布在4 400 m标高以上,在4 400 m标高以下出现As、Pb和F元素异常;主成矿元素Au、Ag、Bi、Sb、As、Cu、Mn、Co和Sn等元素异常均与矿体基本吻合,深部未封闭,且往NW向深部呈延伸趋势,个别元素异常在深部逐渐变窄或尖灭,在矿体部位具有明显的异常内、中、外带,其中W、Mn、Co和Hg元素异常出现在矿体附近,F、Pb和As元素异常主要分布在矿体下部。
(2)元素组合特征
聚类分析是通过对变量进行分类分析,从而了解各个变量及变量组合之间的亲疏程度(
An et al.,2021),本文采用R型聚类分析。
通过对达热尔金矿63号勘探线2个钻孔所收集的209件岩芯样原始数据进行R型聚类分析(
图4),以15个元素的相关水平来确定其成矿阶段及矿物共生组合特征。首先,以相关系数0.5为界,Co、Mn、Cu、Ag、Bi、Au、Sb、As、Pb和Sn构成一个大类,基本代表了成矿元素组合,F、Zn、Mo、W和Hg等元素较为分散,各自单独为一类。然后,以相关系数0.7为界进一步细分,第一大类中又划分为Co、Mn、Cu和Ag、Bi、Au、Sb共2个亚类,其余元素则单独为一类。第二亚类中Ag和Au为主要成矿元素,其中,Ag-Bi和Au-Sb相关性分析结果较一致,说明该类元素关系密切;同时(Ag-Bi)(Au-Sb)As与(Co-Mn)Cu相关性分析结果亦相近,与Pb-Sn-F相关性较弱,反映出金矿化热液活动除了与含Pb-Sn中—高温热液具有一定的联系外,还与(Ag-Bi)(Au-Sb)As中—低温含矿热液活动有关。由此可见,Au的成矿热液活动是多期次的。(Zn-Mo)之间相关较紧密,但与其他元素不相关,反映Zn和Mo元素可能与矿区高背景地层及后期热液再富集有关;W和Hg元素与其他微量元素无聚类关系,且相关性较弱,说明这2个元素与成矿无紧密联系。
李惠等(2006)综合研究提出,我国金矿床原生晕综合轴向分带序列从上至下为:B-I-As-Hg-F-Sb-Ba→Pb-Ag-Au-Zn-Cu→W-Bi-Mo-Mn-Ni-Cd-Co-V-Ti。结合R型聚类分析结果,认为达热尔金矿床63号勘探线AuⅡ矿体的近矿晕元素为Au、Cu、Pb、Zn和Ag,前缘晕元素为F、Sb、As和(Hg),尾晕元素为Mo、Co、Bi、Mn、Sn和(W)。
(3)原生晕轴向分带特征
研究原生晕的轴向分带序列特征,对于矿体在地表的剥蚀程度研究和深部延伸情况判断具有指导作用,从而对深部盲矿体预测起到重要指示意义(
魏江等,2024)。本文采用C.B.格里戈良法原生晕轴向分带序列,通过运算得出原生晕的分带指数,在此基础上计算出元素变异性指数(
表2),从而确定原生晕元素轴向分带序列。经计算可得达热尔金矿床63号勘探线的原生晕轴向分带为(由浅至深)Zn-W-Hg-Ag-Sn-Cu-Mo-Co-Au-Sb-Mn-F-Pb-As-Bi。
前人研究认为,矿体附近聚集的往往是与成矿元素相关性较好的元素或主成矿元素,因此越容易富集在矿体前部的微量元素,其活动性越强(如 Hg、As和Ba等),而在矿体尾部出现的微量元素(如W、Bi和Mo等)活动性较差(
李惠等,2006)。但是,达热尔金矿床的原生晕轴向分带序列却有所不同,显示出鲜明的“反向分带”特征,前缘晕元素F、As和Sb等反常出现在分带序列下部,而尾晕元素Mo、Co、W、Bi和Sn分布混杂,且出现在接近矿体的部位,认为可能是下一个成矿期次的开始,预示着深部可能存在有盲矿体或有新阶段矿体的叠加,表明矿体向深部仍有较大的延伸。
(4)矿体剥蚀系数评价
矿体剥蚀系数通常用前缘晕元素组合与尾晕元素组合的累加比值或累乘比值来表示,根据原生晕轴向分带序列并结合研究区地球化学特征,对元素分布情况进行研究,以此判断是否存在隐伏矿体。矿体剥蚀系数是确定矿体在地表的剥蚀程度,以及预测深部成矿能力的关键指标,该指标可以定量反映矿体在空间上的变化规律(
段启超等,2020)。选取并对比研究具有一定组合关系的成矿元素,包括前缘晕和尾晕特征元素,从而得到原生晕轴向地球化学参数的变化特征,以此评价矿体在深部的含矿性。
通过对达热尔金矿区AuⅡ矿体前缘晕元素与尾晕元素间比值、累加比值及乘积比值等进行计算(
叶庆森,2014),其中3种比值的线性关系较明显,分别为
w(F)/
w(Co)、
w(As+Sb+F)/
w(Bi+Co)、
w(As×Sb×F)/
w(Bi×Co),得出如
图5(a)所示的变化曲线。由
图5(a)可以看出:达热尔金矿地球化学参数的变化曲线由浅至深具有明显的“降—升”特征。
4.2 地球物理特征
鉴于可控源音频大地电磁测深(CSAMT)法对深部信息的采集能力,结合区域相邻矿区找矿经验,达热尔矿区无论岩性或物性均适宜开展相关电磁法测深工作。因此,选择矿区内工程控制程度较高的AuⅡ和AuⅢ矿带,开展可控源音频大地电磁剖面(CS63)测量工作,设计剖面长度为1 km,点距为20 m。获得的测量成果反演解译图[
图5(b)]显示:该剖面自西向东分别在420 m、540 m和780 m处呈现多条低阻带。其中,420 m处低阻带对应地面已知的F
4断层构造(AuⅢ矿带),540 m处低阻带对应地面已知的F
5断层构造(AuⅡ矿带),这2处异常规模大且向深部延伸长,与之对应的含矿构造蚀变带规模大,有望扩大找矿规模。结合勘查资料,综合矿区地质背景和岩(矿)石电物性特征等信息,推断该剖面780 m处低阻带也为隐伏断裂,应引起重视。
此外,青海省有色第二地质勘查院(
袁克全等,2021)在矿区实施了1/5 000激电剖面测量,并圈定出5处激电异常(编号为
ηs-2、
ηs-3、
ηs-5、
ηs-6和
ηs-7),其中
ηs-2和
ηs-3异常分别对应AuⅡ和AuⅢ矿带。2处异常均呈NNE向条带状展布,
ηs-2异常长度约为2.1 km,宽度为30~70 m,极化率为3.4%~6.0%,最高值为6.89%,电阻率为974~2 086 Ω,最低值为455.59 Ω,整体表现为相对中低阻、高极化特征;
ηs-3异常长度约为2.1 km,宽度为30~100 m,极化率为2.33%~5.69%,电阻率为692~1 915 Ω,整体表现为相对中低阻、高极化特征。在对应异常内开展了1/5 000土壤剖面测量工作,圈定出Au和Ag元素高含量段。从物化探异常展布和矿带分布来看,物探异常比化探异常宽,物化异常基本套合,物化探异常展布方向与矿带方向一致,后期通过探矿工程发现并圈定出金矿体。
4.3 深部找矿预测
一般情况下,前缘晕元素的含量随着矿体向深部延伸而逐渐降低,相应的剥蚀系数也呈有规律的下降趋势。若剥蚀系数从下降向上升产生转折变化,则预示着叠加有多期次成矿过程,表明在深部存在盲矿体或矿体向下继续延伸。达热尔金矿床AuⅡ矿体的原生晕轴向分带序列具有“反向分带”特征:地表至深部标高4 400 m,矿体前缘晕元素F、As和Sb在矿尾晕出现,而尾晕元素W、Mo和Mn出现在前缘,存在前缘晕、近矿晕和尾晕元素叠加;在标高4 400 m以下,出现F和As前缘晕元素,Co、Mo、Bi、Mn和Sn尾晕未出现,说明矿体深部沿轴向仍有延伸。同时,深部出现了较强的Pb元素异常晕和弱的Au异常晕,说明矿体向深部仍有延伸;矿体剥蚀系数由浅至深表现出“降—升”的波状变化特征,说明矿体经历了多期次成矿过程,进而表明深部存在盲矿体或矿体向深部仍有较大延伸。可控源音频大地电磁测深剖面(CS63)探测成果显示,AuⅡ矿带所对应的含矿构造蚀变带异常规模大且向深部延伸长,具有较大的找矿潜力。达热尔金矿区ZK6302钻孔深部仅见矿化,未圈定出矿体,但该区矿体具有尖灭再现的特征,结合原生晕特征认为深部存在盲矿体的可能性比较大。
结合物化探异常,在63号勘探线南北两侧的55号和71号勘探线开展了钻孔验证,均显示有矿体,金品位最高达27×10-9,进一步验证了达热尔矿区深部具有良好的成矿远景。
5 结论
(1)对达热尔金矿床控制的AuⅡ矿体63号勘探线剖面209件原生晕样品的分析测试结果进行R型聚类分析,确定近矿晕元素为Au、Cu、Pb、Zn和Ag,前缘晕元素为F、Sb、As和(Hg),尾晕元素为Mo、Co、Bi、Mn、Sn和(W)。
(2)根据格里戈良分带指数法,确定达热尔AuⅡ矿体原生晕自上而下的轴向分带次序为Zn-W-Hg-Ag-Sn-Cu-Mo-Co-Au-Sb-Mn-F-Pb-As-Bi,显示出鲜明的“反向分带”特征,预示着深部可能存在盲矿体或有新阶段矿体的叠加。
(3)矿体原生晕轴向分带序列具有反向分带现象,矿体前缘晕元素F、As和Sb在尾晕出现,尾晕元素W、Mo和Mn出现在前缘,存在前缘晕、近矿晕和尾晕元素叠加现象;矿体剥蚀系数由浅至深表现出“降—升”的波状变化特征,表明矿体向深部具有一定延伸或深部存在盲矿体。
(4)结合可控源音频大地电磁测深等方法,推测深部可能存在矿(化)体,通过钻探验证发现了矿体,说明原生晕与物探方法相结合开展综合找矿,取得了良好效果。
青海省地质勘查基金项目“青海省柴周缘主要矿集区矿产资源深部勘查方法技术示范研究”(2021074008ky008)