东昆仑成矿带是我国著名的金腰带,具有多阶段成矿活动、成矿物质多期叠加的特征(
王成勇等,2018)。其中,东昆仑东段先后发现了瓦勒尕、五龙沟、巴隆、阿斯哈、开荒北和果洛龙洼等一批大中型金矿床(
闫臻等,2000;
肖晔等,2014;
国显正等,2018;
张爱奎等,2021;
梁改忠等,2022;
张勇等,2023)。东昆仑西段已发现的金矿床数量少且规模小,研究人员从矿床地质特征和找矿潜力等方面对东昆仑西段少数金矿床进行了研究,初步总结了该区域金矿成矿规律(
潘彤等,2003,
2008;
李金超,2017;
张爱奎等,2021;
刘永乐等,2022)。截至目前,东昆仑成矿带西段尚未有大型金矿床的报道,达到中型规模的金矿床仅有它温查汉西和肯德可克,其余基本为小型或矿点规模,这与该地区优越的成矿地质背景和物化探异常不相匹配。
自20世纪90年代初开始,一直到2010年,研究人员在茫崖河东及其周边地区主要围绕矽卡岩型铁矿开展工作,先后发现了全红山、那东和茫崖河东铁矿(
李海宾等,2018),很少开展金多金属矿的勘查评价工作。2010—2012年,前人在那陵格勒河中游地区开展的1/5万水系沉积物测量工作(
张伟等,2023),圈定了多处以Au为主的综合异常,但异常查证效果不明显,针对金矿找矿获得的有用找矿信息不足,地表找矿效果不理想。2015年,通过开展青藏专项大调查项目——“青海省格尔木市那东—黑石沟地区铜多金属矿调查评价”,在研究区发现了一些金矿化线索,但是由于受找矿方法和技术手段的限制,金矿找矿一直未能取得重大突破。
鉴于研究区具备良好的金矿成矿条件,为实现浅覆盖区金矿找矿突破,亟需建立一套高效的地球化学勘查方法体系。1∶2.5万大比例尺地球化学测量作为一种成熟有效的地球化学勘查手段,在青海省东昆仑成矿带乌兰乌珠尔、沟里和五龙沟等地区取得了显著找矿成效,尤其在金矿勘查工作中发挥着重要作用(
邱瑜等,2019)。综合研判研究区地球化学景观特征,本次采用1∶2.5万地球化学测量与1∶1万土壤测量相结合的大比例尺化探方法,通过分析元素地球化学和共生组合特征,并结合成矿地质条件,圈定出具有找矿指示意义的综合异常。同时,择优开展异常详细查证,新发现5条含金构造蚀变带,圈定金矿体18条,初步估算矿床规模达中型,实现了该地区金矿找矿的首次重大突破。通过进一步开展构造叠加晕分析,证实研究区金矿主矿体深部具有较大的延伸空间。该套找矿方法组合有效揭示了不同尺度地球化学测量对金矿勘查的指示作用,取得了良好效果,为区域找矿工作提供了示范。
1 地质背景
1.1 区域地质背景
茫崖河东金矿大地构造位置处于秦祁昆造山系—东昆仑造山带—昆北复合岩浆弧内(
祁生胜等,2019)。Ⅳ级成矿区带处于秦祁昆成矿域、昆仑(造山带)成矿省、东昆仑成矿带(青海段)(Ⅲ-26)祁漫塔格—都兰成矿亚带(Ⅳ-26-1)(
潘彤等,2022;
何书跃等,2025;
李少南等,2025)[
图1(a)]。区内出露的主体地层主要为奥陶—志留系祁漫塔格群碎屑岩、火山岩、碳酸盐岩及古元古代金水口群片岩和片麻岩等,局部出露早泥盆世陆相火山岩和少量早石炭世石拐子组(
安勇胜等,2008)。区内侵入岩分布广泛,多以中—晚三叠世和晚泥盆世中酸性岩为主。由于受区域性断裂的影响,区内主要构造线多呈NW-SE向展布(
张伟等,2023)。
1.2 矿区地质
区内出露的主体地层为奥陶—志留系祁漫塔格群,为一套弧后盆地深海热水沉积岩、火山岩、碳酸盐岩和浊积岩建造(
李少南等,2025)。其中,与金矿关系较为密切的祁漫塔格群火山岩组所在的巴音郭勒—苏海图火山岩带是一个相对独立的次级火山盆地,火山岩出露规模相对较大,在区域上有较大的延展性,玄武岩地层体总体呈NWW向展布,火山岩厚度由西向东具有减薄的趋势(
王秉璋,2012)。
区内金矿带南、北两侧侵入岩大面积发育,岩石类型较为复杂(
张爱奎等,2015),成带性好,岩性主要为晚三叠世二长花岗岩和正长花岗岩(
马文君等,2024)。金矿带北东部出露少量印支期花岗闪长岩[
图1(b)],出露面积约为1.03 km
2,呈NWW向带状展布,与围岩祁漫塔格群碎屑岩组呈侵入接触关系。该岩体与地层接触面可见黄铁矿,晶形完整。
由于受区域深大断裂的影响,研究区内断裂发育,通过开展1∶1万地质草测、高精度磁法解译及后期地表系统追索工作,划分出66条断裂,主要为NW向深大断裂和NW向次级断裂,其次为NE和近EW向断裂,NW向断裂是区内主要导矿和储矿构造(
李少南等,2025)[
图1(b)]。断裂附近岩石较破碎,形成构造破碎蚀变带,带内发育构造角砾岩、碎裂岩及黄褐色断层泥。蚀变带局部区段可见由断裂活动挤压及牵引形成的扭曲构造。该地区已发现茫崖河东金矿、茫崖河东铅锌矿、茫崖河东铁矿和茫崖河东钴镍矿,以及较好的钨矿化线索。
通过开展1∶20万地球化学测量工作(
安勇胜等,2008),在茫崖河东金矿及其周边圈定出AS
CuV-CoNiAu(那陵格勒河)和AS
CuSbVCoNiAu(开木棋河西)2处综合异常。异常内NW-SE向断裂发育,出露较大规模的中、晚三叠世、晚泥盆世中酸性侵入岩体,出露地层为奥陶—志留系祁漫塔格群(
马文君等,2024)。Au元素峰值为5.5×10
-9~14.5×10
-9,面积为32~80 km
2,Au、As和Sb元素套合较好。
通过开展1∶5万化探工作,在茫崖河东金矿及其南东部圈定了中央山HSCu(V-Ni-Au-Co-As)和红卫山HSAu(W-V-Ni-Sb-Cu-Co-As)2处以Au为主的综合异常,金异常具有面积大(14.40~44.54 km2)和峰值高(15.1×10-9~257.9×10-9)的特点,总体呈NW-SE走向带状展布,与奥陶—志留系祁漫塔格群在区位上吻合,异常重现性较好。茫崖河东金矿主矿带位于中央山HSCu(V-Ni-Au-Co-As)综合异常Au20子异常浓集中心。
2 地球化学特征
2.1 1∶2.5万地球化学测量样品采集及分析测试
水系沉积物测量作为一种高效的化探找矿方法,被广泛应用于区域地球化学勘探及矿产普查工作中(
胡兆国等,2018;
徐云峰等,2021;
梁东等,2023)。本次采样严格按照《青海省1∶25 000地球化学测量规范》执行。遵循最有效控制汇水域、以有限的采样点取得最佳效果以及大致均匀的原则,并充分考虑研究区景观特点及自然地理条件后扣除不可采区。在样品采集过程中,尽可能避开风成沙和有机质的干扰,使样品最大限度地代表采样单元基岩的化学成分。对采集样品进行系统编号,并做好野外定点、标记及影像资料的采集。加工后的分析样重量均大于150 g。本次工作共采集水系及岩屑样品4 321件(不含50件重复样),采样粒级选择-10~+60目,完成采样面积为210.0 km
2,采样平均密度为20.8个/km
2。样品分析测试工作由青海省地质矿产测试应用中心完成,共分析15种元素,各项指标符合规范要求,样品分析数据真实、质量可靠(
王振等,2023)。
2.2 1∶2.5万地球化学异常特征
为了查明茫崖河东地区地球化学特征,圈定地球化学异常,进一步缩小金矿勘查靶区,在长山—红卫山地区开展了1∶2.5万地球化学测量工作,全区共圈定综合异常39处,其中本次研究区内圈定Au单元素异常63处,以Au为主元素的综合异常5处。
(1)元素组合特征
元素组合是反映研究对象地球化学特征的一组特殊指示元素,能够为找矿勘查提供有效的指示元素(
张嘉升等,2023),为后续地球化学异常的准确解释、评价及确定矿(化)体成因类型提供理论依据。R型聚类分析(
席明杰等,2013)可以揭示各因子元素之间的相关关系,从而去除一些冗余的地球化学信息,将能够反映研究区主要成矿特征的地球化学元素组合类型信息提炼出来,定量研究元素之间的聚集和分离(
于俊博等,2014)。
为了解主要地质环境中元素聚集趋势和内在成因关联,掌握元素之间的组合关系,揭示区内经历的主要地质活动及作用(
赵娟等,2021),本次利用全区Au、As和Sb等15种元素的原始数据集进行R型聚类分析和因子分析,发现2种方法的处理结果具有共同之处,均揭示了某种特定的地质环境和地质作用。进一步利用GeoExpl软件(多元地学空间数据管理与分析系统)对研究区内15种元素分析结果进行R型聚类分析,选取累计方差贡献达80%以上的8个主要因子,各因子变量载荷以0.6的绝对值为界,按大小排序形成结构式(
表1)。区内地球化学元素分布特征与区内地质构造展布关系密切,通过应用旋转因子成果表达测区地球化学变化特征(
邱炜等,2015),结合R型聚类分析、区内地质背景和特征元素变化系数等参数特征进行分析,R型聚类分析图谱如
图2所示。
由
表1可知,F4因子特征根为1.2113,特征根百分比为8.0756%,是本区第4个因子,元素组合为Au和As,载荷值均大于0.7,反映了与区内低温构造活动有关元素的因子。结合因子特征及其他地球化学特征可知,异常与区内奥陶—志留纪祁漫塔格群火山岩组及其内部发育的NW-SE向构造关系密切,属于奥陶—志留纪祁漫塔格群和NW向断裂构造控制的矿致异常。茫崖河东研究区已发现5条含金矿化蚀变带,金矿床规模达中型以上,金矿成矿事实清楚。
(2)Au元素分布及富集特征
元素变化系数反映元素含量的相对离散程度(
陈巧妹等,2019)。成矿作用是一个长期、多期且复杂的地质过程,导致不同元素在各期次地质体中的富集或贫化现象不尽相同。这种富集或贫化现象是基于元素背景值基础上的高低变化,是地质—地球化学共同作用的结果(
马文君等,2024)。通常运用富集系数和变异系数等参数体现成矿过程中每个阶段地球化学元素的成矿能力。研究区内Au、As和Sb元素地球化学变量统计结果见
表2。
全区4 371件样品中Au元素含量平均丰度值为2.02×10
-9,大于青海省全省Au元素含量平均丰度值(1.35×10
-9)(青海省三轮区划水系沉积物统计样本容量为73 778件)。根据变异系数(
CV<1.0,1.0≤
CV≤2.0,
CV>2.0)和富集系数(
EF<1.25,1.25≤
EF≤2.0,
EF>2.0)双重指标对元素进行分类,将区内15种元素划分为均匀元素、不均匀元素和很不均匀元素。由
表2和
图3可知,Au、As和Sb元素变异系数均大于2.0,呈强分异性(其中,奥陶纪祁漫塔格群火山岩组中Au元素变异系数达12.12);富集系数均大于2.0,呈富集状态。Au、As和Sb元素具有富集程度高,分布极不均匀的特性,成矿可能性极大,这与目前在该地区发现的矿种十分吻合。
进一步分析茫崖河东矿区各元素含量概率发现,Au、As和Sb元素呈现出一定的偏态分布,但元素峰值较高,证实了与深大断裂有关的金属成矿元素富集的事实,通过多期地质—地球化学作用的叠加,成矿可能性较大。
(3)Au元素时空分布特征及变化趋势
为了解茫崖河东地区Au、As和Sb元素在不同地质单元的分布情况,通过水系沉积物数据对各地质单元中元素分配情况进行统计,结果见
表3。Au元素在金水口群(Pt
1J)、晚泥盆世牦牛山组(D
3m)及中酸性岩体中的含量较低,普遍低于1.0×10
-9,祁漫塔格群碎屑岩组(OSQ
1)和碳酸盐组(OSQ
3)内Au元素含量平均值为1.9×10
-9~2.2×10
-9,与全区平均值接近,而火山岩组(OSQ
2)内Au元素含量平均值为5.4×10
-9,接近全区平均值的3倍(
表3和
图4)。同时,在茫崖河东金矿采集的光谱样品测试结果显示,祁漫塔格群火山岩组内Au元素丰度值为6.14×10
-9,普遍高于祁漫塔格群碎屑岩组内Au元素丰度值(1.29×10
-9~2.64×10
-9),与构造作用有关的片理化蚀变安山质玄武岩(9.83×10
-9)、碎裂状蚀变安山质玄武岩(100.36×10
-9)和构造角砾岩(45.25×10
-9)相较完整的蚀变安山质玄武岩具有更明显的Au元素富集趋势,区内发育的各类侵入岩体和脉体Au元素丰度值为0.57×10
-9~3.35×10
-9,普遍偏低,说明区内祁漫塔格群火山岩组为金的初始矿源层(
李少南等,2025)。
全区Au元素平均值为2.02×10
-9,最大值为1 710×10
-9;As元素平均值为21.74×10
-6,最大值为4 343×10
-6;Sb元素平均值为1.39×10
-6,最大值为152×10
-6。全区共圈定Au异常65处,其中发育外中带或外中内带的异常有33处(
表4和
图5),主要分布在茫崖河东中央山、金牛滩和黑石沟一带出露的祁漫塔格群中,异常分布形态及区位与该套地层套合性较好。异常强度较高,多具有三级浓度分带,主要伴生元素为Cu、Pb、Zn、Ni和Co。茫崖河东金矿一带异常面积大,带状分布明显,沿断裂两侧分布,展布方向呈NW-SE向,与断裂走向一致。通过衬值累加地球化学分析发现,Au、As和Sb地球化学场变化规律基本一致,在全区呈带状分布(
图6)。高背景浓集区主要沿NW向断裂浓集分布,表明在该范围Au、As和Sb低温热液矿化活动强烈,该地段是寻找构造蚀变岩型金矿的前景区。
(4)典型异常特征
通过对研究区圈定的5个以Au为主的综合异常进行了查证,发现奥陶—志留系祁漫塔格群火山岩组(OQ2)及该套地层中发育的断裂是Au异常产生的原因。就全区而言,GAAu(Cu-Ni-As-Co-Pb)综合异常评序第一。通过地表追索及实施槽探、钻探工程后,发现了5条含金构造蚀变带,蚀变带内见矿效果好,充分说明研究区具有较好的成矿地球化学和地质条件,有必要作为典型异常进行研究。
GA
Au异常区分布在研究区东北部中央山西段,呈带状NW-SE向展布。综合面积约为6.94 km
2。异常区北侧主要出露奥陶—志留系祁漫塔格群碎屑岩组,岩性主要为硅质黏土质板岩、硅质板岩和硅质岩;异常区中南部出露火山岩组,岩性为灰绿色蚀变安山质玄武岩。异常内侵入岩很少发育,零星可见长英质岩脉露头。断裂十分发育,以NW-SE向为主,亦发育有少量NE向断裂(
图7)。
异常主元素为Au,特征组合元素为Cu-Ni-As-Co-Pb。Au异常由Au12、Au13、Au16、Au17、Au20和Au22共6个子异常组成。其中,Au13子异常,异常点数有8个,异常面积为0.39 km
2,Au元素峰值为1 710×10
-9,平均值为332.73×10
-9;Au22子异常,异常点数有38个,异常面积为1.80 km
2,元素峰值为38.1×10
-9,平均值为9.24×10
-9。As异常由As9、As12、As13、As15、As18和As21共6个子异常组成。其中,As9子异常,异常点数有3个,异常面积为1.07 km
2,元素峰值为659×10
-6,平均值为95.07×10
-6;As13子异常,异常点数有3个,异常面积为0.16 km
2,元素峰值为1 563×10
-6,平均值为862.3×10
-6。其余各元素也具有多个子异常,异常规模大小均有所不同,各元素之间套合性较好。异常特征见
表5。
综上所述,1∶2.5万地球化学测量工作基本确定了Au异常的范围,显现了金矿找矿的地球化学空间和潜力。当前,异常查证中发现的金矿(化)体矿化线索,在空间分布上与水系沉积物地球化学测量所圈定的Au单元素异常呈现显著耦合性。异常晕的展布范围受控于地质体边界,其延伸方位与区域地层走向呈现良好的一致性。同时,其空间配置与NW-SE向断裂构造系统存在成因联系,综合显示典型层控—构造复合型矿床的成矿特征。
2.3 1∶1万地球化学特征
针对区内1∶2.5万GA
Au(CuNiAsCoPb)等多处综合异常,进一步开展了1∶1万土壤剖面测量进行解剖,并结合以往研究资料(
杜承文等,2016),在茫崖河东金矿区累计圈定了以金、铜为主的1∶1万土壤综合异常12处。异常具有规模大和强度高的特点,其中区内发现的含金矿化带及主矿体均穿越AP3综合异常Au异常浓集中心。
AP3异常呈不规则状NW-SE向展布,主元素为Au,组合元素为As、Sb、Cu、Pb、Zn和Ag,各元素分带明显,特征值见
表6。异常区出露地层为奥陶—志留系祁漫塔格群火山岩组(OSQ
2 )。Au异常面积为0.36 km
2,峰值高达1 890×10
-9。Au与As、Sb元素套合好(
图8)。该异常与C3-1磁异常南东对应性较好。经33条槽探工程检查,新发现的Ⅰ、Ⅱ和Ⅳ号含金构造蚀变带主体位于该异常内,异常范围内3条金矿带中目前圈定金矿体(Ⅰ-1、Ⅰ-2、Ⅰ-3、Ⅱ-1、Ⅱ-2、Ⅳ-1、Ⅳ-2和Ⅳ-3)8条,锌矿体1条(Ⅱ-3),该异常属矿致异常,其赋矿岩性为构造角砾岩、构造片理化蚀变安山质玄武岩和糜棱岩,赋矿围岩为蚀变安山质玄武岩,矿体属构造蚀变岩型。目前3条含金构造蚀变带及带内矿体控制程度普遍不高,因此该异常内具有较好的金矿找矿潜力。
1∶1万土壤测量金异常,基本定位了含金构造蚀变带的展布空间,并为地表金矿体的发现奠定了基础。
2.4 构造叠加晕特征
金矿体沿倾斜方向上延伸稳定,斜深最大为340 m的主矿体往深部均未圈闭[
图9(a)],优选主矿体开展了构造叠加晕研究工作,自地表槽探工程至深部钻探工程共采集10件样品进行分析测试(
表7)。对比全国及东昆仑五龙沟金矿床构造叠加晕模式(
魏占浩等,2015;
刘正军,2019;
李惠等,2021),厘定出研究区前缘晕地球化学标志组合为As-Sb-Hg;近矿晕特征指示元素组合为Au-Ag-Cu-Pb-Zn;尾晕特征指示元素组合为 Co-W-Mo-Ni。其中,Ⅰ号含金构造蚀变带南东Ⅰ-4和Ⅰ-6矿体自3 500 m至3 100 m高程段,深部近矿晕呈内带或强带异常,前缘晕呈外带,尾晕呈中内带,指示矿体(Ⅰ-6)往深部仍有很大延伸;Ⅰ号含金构造蚀变带中部Ⅰ-3矿体自3 350 m至3 100 m高程段,深部近矿晕呈中内带异常,前缘晕与尾晕呈中内带异常共存,指示深部有盲矿体;Ⅱ号含金构造蚀变带中部Ⅱ-3矿体自3 410 m至3 175 m高程段,深部近矿晕呈中内带异常,前缘晕与尾晕呈中内带异常共存,指示深部有盲矿体。
3 金矿找矿效果及矿体特征
3.1 找矿效果
以往研究表明,在茫崖河东地区开展的1∶20万和1∶5万多尺度地球化学测量工作,揭示出具有重要找矿指示意义的Au元素地球化学组合异常,其异常强度与规模均达到Ⅱ级勘查异常标准。近20年来,该地区以化探异常为依托的金矿找矿工作持续进行,但异常源追溯一直未取得突破,找矿成效不显著。基于茫崖河东地区优越的成矿地质条件和物化探综合异常,近年来以金矿等战略矿产为重点,陆续开展了1∶2.5万地球化学测量及异常查证工作,旨在解剖1∶20万和1∶5万化探异常,逐步缩小找矿靶区。在茫崖河东地区圈定了多处以Au为主的综合异常,通过进一步查证,确定了异常源,累计发现含金构造蚀变带5条,蚀变带长度为0.5~5.2 km,宽度为0.8~9.9 m,走向NW,倾向SW。带内共圈定金矿体18条,金矿体长度为140~1 300 m,真厚度为0.80~9.34 m,矿体金品位为0.8×10-6~26.4×10-6,金矿床平均品位为2.27×10-6,找矿效果显著。
3.2 主要矿化带特征
(1) Ⅰ号含金构造蚀变带。Ⅰ号含金构造蚀变带沿区内F27断裂形成,呈NW向展布,走向延伸约为4.6 km,宽度为0.8~9.9 m,倾向SW,倾角为40°~80°,断层性质为逆断层。该构造带分布于祁漫塔格群中,其中北西侧位于火山岩组安山质玄武岩与碎屑岩组硅质黏土质板岩接触部位,围岩主体位于火山岩组安山质玄武岩内。带内普遍发育褐铁矿化、高岭土化、硅化和黄钾铁矾化,局部具黄铁矿化、碳酸盐化和毒砂矿化,且带内岩石褪色蚀变现象明显。
硅化地表以石英宽脉(1~20 cm)或石英透镜体形式分布,往浅深部硅化以石英细脉和网脉状形式分布。经槽探和钻探工程揭露验证,在该带北西部圈定出Ⅰ-1和Ⅰ-2金矿体;中部圈定出Ⅰ-3金矿体;南东部圈定出Ⅰ-4和Ⅰ-5金矿体,南东深部圈定出Ⅰ-6金矿体(
李少南等,2025)。Ⅰ号含金构造蚀变带矿体特征见
表8。
(2)Ⅱ号含金构造蚀变带。该蚀变带沿F
30断裂形成,与Ⅰ号含矿带平行,呈NW向展布,倾向SE,倾角为55°~78°,带宽为1.44~12.5 m,延伸长度大于5.2 km。断层性质为逆断层。构造带分布于祁漫塔格群火山岩组中,构造带围岩为蚀变安山质玄武岩。带内岩石普遍具较强的褐铁矿化、高岭土化和硅化(
图10)。硅化地表以石英宽脉(1~10 cm)或石英透镜体形式分布,往浅深部硅化以石英细脉和网脉状形式分布。在15、28、64、272勘探线附近各有1条NE向断裂穿切,12勘探线南东发育多条分支带,76勘探线附近Ⅰ、Ⅱ号含金构造蚀变带会合形成一膨大部位,通过捡块化学样和刻槽化学样对该带及其分支带进行厘定。224线南东含矿带倾角逐渐变缓。经槽探和钻探工程揭露控制,该带北西部圈定出Ⅱ-1金矿体,中部圈定出Ⅱ-2金矿体,南东部圈定出Ⅱ-4和Ⅱ-5金矿体(
表9)。含矿岩性为构造角砾岩和糜棱岩,推测原岩为蚀变安山质玄武岩(
李少南等,2025)。
4 结论
(1)采用1∶2.5万地球化学测量方法开展地质找矿工作,能够快速缩小找矿靶区,提供较直接且可靠的找矿信息。在茫崖河东地区圈定出65处Au异常,5处以Au为主元素的1∶2.5万综合异常,Au异常内发现了规模不等的含矿构造带,带内圈定出金矿(化)体,且矿床规模已达中型,取得了良好的找矿效果。
(2)不同尺度地球化学测量为该金矿勘查提供了不同阶段的找矿信息。1∶2.5万地球化学异常反映出区内Au元素地球化学空间和找矿潜力。1∶1万土壤异常基本精确定位了含金构造蚀变带的展布空间,为圈定金矿体奠定了基础。构造叠加晕对深部找矿具有较好的指示作用,可作为深部找矿工程部署的依据之一。
(3)茫崖河东地区Au、As和Sb异常发育于茫崖河东—金牛滩—黑石沟一线,异常呈NW-SE向带状展布,延伸长度超过22 km。茫崖河东金矿位于该异常带的北西部,异常带中部和南东部仅开展了简单的异常查证工作,找矿空间较大。通过实施系统的地质—地球化学—地球物理综合勘查工作,有望在茫崖河东地区实现金矿找矿重大突破。
青海省地质矿产勘查开发局基金项目“东昆仑奥陶纪钴矿成矿环境及找矿潜力”(青地矿科[2024]67号)
“青海省格尔木市长山—红卫山地区1∶2.5万地球化学测量”(青地矿科[2022]6号)