胶东地区是我国最重要的金资源基地,累计探获金资源储量超过5 000 t(
宋明春,2015;
杨立强等,2014,
2024;
黄鑫,2021;
于晓卫等,2023;
张泽涛等,2025)。胶西北地区位于胶东半岛西北部,区域内分布有3条近平行分布的NNE-NE向断裂带,自西向东依次为三山岛、焦家和招平断裂(
肖风利等,2018),随着深部找矿工作的持续推进,该地区已成为找矿突破的重要区域,金资源储量占胶东地区已探明金矿床的90%(
宋明春等,2011;
朱平平等,2023;
冯雅杰等,2025)。
谢家沟金矿位于招远市西南部,产于招平断裂带以西5 km的下盘次级谢家沟断裂带上,严格受控于NNE向韧性剪切带,属于同韧性剪切带蚀变岩型矿化类型(
辛洪波等,2007;
姚晓峰等,2020),区域内大面积出露花岗岩,脉岩多呈近SN向分布。前人在谢家沟金矿开展了控矿构造、矿床地质特征和成矿时代等大量工作(
田锋,2005;
张群喜等,2007;
辛洪波等,2006,
2007;
姚晓峰等,2020;
王文钰等,2021;
鲍兴隆等,2024),但对成矿前后岩脉的地球化学研究相对较少,导致对末期岩浆活动的性质认识不足。为此,通过对谢家沟金矿的两套岩脉开展系统的地球化学研究,并结合现有研究成果,探讨了岩脉与金矿化的潜在关系,进一步明确了谢家沟金矿的成矿地质背景,为认识胶西北金成矿作用和下一步找矿勘查提供理论依据。
1 矿床地质特征
矿区内地层仅见第四系,主要为残积物、残坡积物和冲积物,其中残积物和残坡积物主要分布在低山丘陵地段,冲积物主要分布在山间河流地带(
丁东胜,2019;
王文钰等,2021)。
矿区内构造较为发育,主要为NNE和NNW向断裂(
图1)。其中,NNE向断裂为矿区内主要断裂,发育规模较大,具有多期活动特征,早期韧性剪切带内可见S-C组构和定向排列现象,与矿床空间展布形态一致,晚期韧性变形带破坏矿体,可见断层泥和碎裂岩等;NNW向断裂规模较小,带内可见断层泥和角砾岩(
姚晓峰等,2020;
鲍兴隆等,2024)。
矿区内发育玲珑花岗岩,主要岩性为片麻状黑云母花岗岩,主要矿物成分为钾长石、斜长石、石英和黑云母等(
胡文萱等,2024),还发育有大量的燕山期岩脉,主要有煌斑岩、辉石闪长岩、闪长玢岩、石英闪长玢岩和花岗闪长岩等,大部分岩脉呈NW和NNW向展布。围岩蚀变主要有钾长石化、绢云母化、硅化和黄铁矿化等。矿体沿NE向展布,倾向SE,倾角为70°~80°,严格受控于韧性剪切带,具有尖灭再现特征。矿石类型主要为黄铁矿(黄铜矿)绢英岩化和石英硫化物脉型(
辛洪波等,2006),其中前者品位较高但空间连续性较差,后者叠加于绢英岩化矿体之上,矿化较为均匀但品位相对较低。矿石矿物主要有黄铁矿、黄铜矿、方铅矿、闪锌矿、自然金和银金矿,还有少量的碲银矿和辉银矿,脉石矿物主要有石英、绢云母、钾长石、碳酸盐矿物和萤石(
杜泽忠等,2020;
姚晓峰等,2020;
Du et al,2023;
鲍兴隆等,2024)。结合矿物共生组合和穿插关系,建立谢家沟金矿的矿物生成顺序如
图2所示。
2 侵入岩特征
矿区主要脉岩类型有闪长玢岩、花岗闪长斑岩、闪长岩、石英闪长玢岩、辉石闪长岩、辉绿岩和煌斑岩等(
图3)(
张群喜等,2007)。目前对于成矿前脉岩缺乏直接同位素测年研究,但
韦延光(2005)通过单矿物及全岩K-Ar法测年,获得谢家沟金矿的主成矿年龄为115~125 Ma,据此推断成矿前脉岩形成时代应早于125 Ma,成矿后脉岩形成时代约为118 Ma(
姚晓峰等,2020)。本研究的成矿前蚀变岩(原岩为闪长岩)具有半自形粒状结构,主要矿物为斜长石(50%~55%)、角闪石(20%~25%)、黑云母(5%~10%)和少量石英(5%)。成矿前细粒花岗岩具有细晶结构,块状构造,主要矿物为钾长石(40%~45%)、石英(30%~35%)和斜长石(15%~20%)。成矿后石英闪长玢岩具有斑状结构,斑晶含量为20%~25%,主要为斜长石(10%~15%)、角闪石(5%~8%)和石英(3%~5%),基质具有显微晶质结构,由长石和石英微晶组成。成矿后花岗闪长岩具有中—细粒花岗结构,主要矿物为斜长石(35%~40%)、钾长石(15%~20%)、石英(20%~25%)、角闪石(10%~15%)和黑云母(3%~5%)。
3 采样及分析测试
3.1 采样位置
选取坑道内两类岩脉进行取样,分别为成矿前细粒花岗岩、蚀变岩和成矿后石英闪长玢岩、花岗闪长岩,采集无蚀变及脉体穿插的新鲜样品,用纯化水冲洗干净,烘干并粉碎至200目以备测试使用。
3.2 测试方法
样品分析测试工作由北京燕都中实测试技术有限公司完成,主量元素采用XRF测试,微量元素采用ICP-MS测试。
主量元素测试:首先将样品研磨至200目以下,称量后加入Li2B4O7(1∶8)助熔剂混合,放置于融样机中加热至1 150 ℃,使其在金铂坩埚中熔融成均一玻璃片体,后使用XRF(Zetium,PANalytical)测试。测试结果保证数据误差小于1%。
微量元素测试:首先将粉末样品称量后放入聚四氟乙烯溶样罐内,加入HF+HNO3,将高压消解罐放置于干燥箱中,温度控制在190 ℃,72 h后取出,经过赶酸并将溶液稀释后上机测试。测试使用ICP-MS(M90,analytikjena)完成,所测数据根据监控标样GSR-2显示误差小于5%,部分挥发性元素和极低含量元素的分析误差小于10%。
4 测试数据
4.1 主量元素特征
主量元素测试分析结果见
表1。测试结果显示,样品烧失量变化范围为2.14%~4.97%,差异性较小。烧失量校正后样品氧化物含量变化如下:
成矿前,细粒花岗岩和蚀变岩主要氧化物含量变化较大。细粒花岗岩SiO2含量为73.06%~76.07%,平均值为74.68%;Al2O3含量为11.93%~14.18%,平均值为13.28%;MgO含量为0.16%~0.32%,平均值为0.25%;CaO含量为1.37%~2.37%,平均值为1.80%,全碱含量为6.89%~7.86%,平均值为7.47%。蚀变岩SiO2含量为56.70%~62.70%,平均值为59.06%;Al2O3含量为11.27%~15.26%,平均值为13.65%;MgO含量为4.60%~6.45%,平均值为5.65%;CaO含量为5.21%~6.41%,平均值为5.81%;全碱含量为6.19%~7.11%,平均值为6.76%。
成矿后,石英闪长玢岩和花岗闪长岩中主要氧化物含量变化较为稳定。石英闪长玢岩SiO2含量为63.12%~66.25%,平均值为64.60%;Al2O3含量为15.95%~16.70%,平均值为16.35%;MgO含量为1.99%~2.63%,平均值为2.31%;CaO含量为3.59%~4.76%,平均值为4.13%;全碱含量为7.00%~7.23%,平均值为7.13%。花岗闪长岩SiO2含量为66.24%~66.98%,平均值为66.60%;Al2O3含量为15.79%~15.96%,平均值为15.85%;MgO含量为1.89%~2.16%,平均值为2.01%;CaO含量为3.00%~3.72%,平均值为3.37%;全碱含量为7.04%~7.68%,平均值为7.43%。
由
图4可知,成矿前,细粒花岗岩位于钙碱性和高钾钙碱性系列,蚀变岩全部位于钾玄岩系列,谢家沟金矿蚀变地球化学研究表明,钾长石化是该矿床早期重要的蚀变类型,K元素发生明显迁入(
丁东胜等,2017;
杜泽忠等,2020);成矿后,石英闪长玢岩和花岗闪长岩全部位于高钾钙碱性系列。在A/NK-A/CNK图解(
图5)中,样品基本全部落入偏铝质岩石区域,仅部分成矿前的细粒花岗岩位于过铝质区域,属于准铝质钙碱性花岗岩类,为I型花岗岩(
图6)。
4.2 稀土、微量元素分析特征
由
图7和
表2可知,成矿前后岩脉微量元素配分模式表现为高度一致性变化,所有样品均表现出明显的Nb和Ta亏损,指示在俯冲背景下,富集Nb、Ta的矿物可能在源区发生残留或在地幔楔受到俯冲流体交代,同时,成矿后岩脉Zr/Hf比值介于成矿前两岩脉之间,且均低于球粒陨石值(约为36.3),可能指示其发生岩浆混合作用,并形成于高压环境。
由
表3和
图8可知,成矿前后岩脉稀土元素均呈右倾模式,富集轻稀土(LREE)元素,亏损重稀土(HREE)元素,稀土元素配分曲线较为平滑,
δEu值为0.77~1.38,多数样品无显著Eu异常,说明斜长石分离结晶不强烈或源区压力较高,
δCe值为0.88~0.91,所以样品均显示无Ce异常或异常极弱,表明岩浆体系处于还原环境。成矿前岩脉稀土元素含量变化较大,总量为53.79×10
-6~227.16×10
-6,LREE含量为41.44×10
-6~190.80×10
-6,HREE含量为10.58×10
-6~36.36×10
-6,LREE/HREE比值为3.36~6.52,除细粒花岗岩的一个样品La
N/Yb
N值较低之外,其余样品比值范围为13.30~17.95。其中,蚀变岩稀土元素总量较高,为175.47×10
-6~276.95×10
-6,LREE含量为145.00×10
-6~190.80×10
-6,HREE含量为30.47×10
-6~36.36×10
-6;细粒花岗岩含量较低,为53.79×10
-6~97.74×10
-6,LREE含量为41.44×10
-6~84.73×10
-6,HREE含量为10.58×10
-6~13.39×10
-6。成矿后岩脉稀土元素含量变化较为稳定,总量为183.18×10
-6~205.46×10
-6,LREE含量为161.84×10
-6~177.85×10
-6,HREE含量为21.34×10
-6~27.73×10
-6,LREE/HREE比值较高(6.41~7.87),La
N/Yb
N比值范围为21.21~27.64。
5 讨论
5.1 谢家沟金矿成矿地质背景
太平洋板块俯冲是胶东地区构造作用及岩浆活动的主要源动力,构造体制转换为成矿提供了有利的地质条件(
韦延光,2005)。早白垩纪早期太平洋板块发生近EW向俯冲,至123~119 Ma,区域应力场向NW-SE向转换,含金流体从深部地幔经拆离断层快速上升;到早白垩纪后期,EW-SE向伸展持续,热液系统进入脆性主导阶段,促使金在浅部地壳大规模沉淀(
Ni et al,2025)。
通过典型金矿对比,有助于理解不同类型金矿对同一区域构造背景的响应差异。胶西北金矿集区位于华北克拉通东南缘,是我国最大的金矿集区,包含三山岛、焦家和玲珑等大型—超大型金矿,成矿时代主要集中在(120±10)Ma(
尹业长,2020),与胶东地区大规模金成矿事件高峰期高度一致(
韦延光,2005;
辛洪波,2005)。谢家沟金矿与焦家金矿、玲珑金矿的成矿时代相近,与古太平洋板块俯冲方向转换时期吻合,均受到区域断裂系统控制,与中生代岩浆活动密切相关,成矿流体均显示壳幔混合特征,反映了幔源物质的参与,但谢家沟金矿严格受控于韧性剪切带,其矿石与围岩来源差异明显,成矿物质可能具有复杂的来源,成矿压力约为2 200×10
5 Pa,远高于焦家金矿的成矿压力,反映其形成深度更大(
辛洪波,2005)。
成矿后岩脉具有高Sr(658×10
-6~937×10
-6)、低Y(11.7×10
-6~15.5×10
-6)和低Yb(1.09×10
-6~1.50×10
-6)特征,(La/Yb)
N为21.21~27.64(20)、Eu无明显负异常,表现出极为明显的埃达克岩特征(
图9)。俯冲洋壳部分熔融形成的埃达克岩需要高压条件(深度大于50 km),下地壳部分熔融形成的埃达克岩则要求源区残留相中有斜长石存在,因此形成深度相对较浅(
Sun et al,2012)。其高Sr低Y、(La/Yb)
N明显高于典型岛弧岩浆岩以及强烈的轻重稀土分馏特征指示源区残留相中有石榴子石存在,石榴子石是高压条件下的特征残留矿物(
Atherton et al,1993;
Rapp et al,1995;张旗等,2001a,2001b)。在SiO
2-MgO图解(
图10)和Sr/Y-(La/Yb)
N图解(
图11)中,成矿后岩脉样品投点均落入下地壳与大洋板块重叠部位,但趋势偏向大洋板块部分。这一分布特征可能是岩浆混合作用和后期流体微弱改造作用的结果,成矿流体中CO
2与金的正相关关系间接指示流体对微量元素存在改造能力(
韦延光,2005),岩石中Sr可能发生迁移或再分配,导致Sr/Y比值变化。谢家沟成矿后岩脉具有相对较高的MgO含量(1.89%~2.63%)和Mg
#值(42~49),实验岩石学研究表明,单纯的下地壳部分熔融形成的熔体Mg
#值通常较低(小于40),较高的Mg
#值指示熔体在上升过程中可能与地幔橄榄岩发生相互作用(
Xie et al,2012),这种相互作用也是俯冲洋壳熔融型埃达克岩的典型特征。
5.2 岩脉与岩浆演化的关系
区域研究表明,壳幔岩浆的混合反应和演化是胶东地区巨量金质来源的主要形成机制(
田瑞聪等,2022)。成矿前蚀变岩的SiO
2含量(56.70%~62.70%)较低,MgO含量(4.60%~6.45%)很高,全铁含量为7.32%~8.45%,指示其来源于未经充分分异演化的岩浆。成矿前细粒花岗岩SiO
2含量为73.06%~76.07%,MgO含量(0.16%~0.32%)很低,全铁含量仅为1.61%~2.85%,与下地壳部分熔融产生的酸性岩浆成分相符。成矿后的石英闪长玢岩和花岗闪长岩成分位于2个端元中间,在哈克图解(
图12)上表现出良好的线性关系,反映了2种不同成分岩浆的混合过程,如果是单一岩浆通过分离结晶作用演化,其元素变化趋势通常呈曲线而非直线(
韦延光,2005),指示谢家沟地区在白垩纪发生岩浆混合作用。轻重稀土分馏程度(La
N/Yb
N比值)在成矿后岩石中更高(21.21~27.64),高于蚀变岩(13.30~17.95)和细粒花岗岩(6.92~17.95),这表明混合后的岩浆可能又经历了进一步的演化,或者混合过程中源区的残留矿物组合发生了变化。
起源于交代岩石圈地幔的岩浆中富含水、硫和碳等挥发分,硫在热液迁移中能够与金形成络合物,因此基性岩浆上侵的过程中,不断地演化和分异使得岩浆中Au元素含量不断升高,在与酸性岩浆发生混合作用后,迅速改变混合体系的物理化学条件,运移至地壳浅部有利构造部位后沉淀成矿。谢家沟金矿主成矿阶段成矿流体的均一温度主要集中在260~340 ℃。
韦延光(2005)提出更详细的阶段划分数据,第一阶段流体包裹体均一温度为320~350 ℃,第二阶段为270~330 ℃。第三阶段为250~290 ℃。
流体包裹体的气相成分分析揭示了成矿流体的挥发分组成。谢家沟金矿流体包裹体的气相成分主要为H
2O和CO
2,从主成矿阶段至成矿晚期,流体中H
2O含量逐渐增加,而CO
2含量逐渐减少。值得注意的是,流体中CO
2的含量与矿石的金品位之间存在一定的正相关关系,表明CO₂在Au元素的运移和沉淀过程中可能扮演了重要角色(
韦延光,2005;
孙莉,2006)。
丁东胜等(2017)研究进一步明确了流体体系的类型,成矿早期流体属于H
2O-CO
2-NaCl体系,主成矿期流体属于H
2O-CO
2-NaCl±CH₄体系。H、O同位素组成是示踪成矿流体来源的有效工具。
孙莉(2006)对谢家沟金矿的石英样品进行了H、O同位素分析,研究结果显示,在Ⅰ~Ⅲ主成矿阶段,石英的
δ18O平均值为10.98‰,与含矿黑云母花岗岩和郭家岭花岗岩的
δ18O值相当,相应的成矿流体的
δ18平均值为4.15‰,
δD
平均值为-86.68‰,成矿晚阶段形成的萤石流体包裹体的
δD
值为-122.71‰,这一数值接近中生代大气降水的值。
5.3 区域脉岩对比及其对成矿背景的指示
关于脉岩与金矿的成因联系,区域研究提出了胶东群的两阶段部分熔融模式(
杨冀,2010)。该模式认为,早期(第一阶段)部分熔融产生中酸性岩浆,此时温度较低,金不易活化迁移,从而在胶东群残留体中富集。后期(第二阶段)在燕山运动强烈的构造—热事件影响下,胶东群再次部分熔融产生中基性岩浆,此时金达到活化迁移的条件,并在幔源流体的参与下进入岩浆系统,含矿岩浆在上升过程中,成矿流体与中基性岩浆发生分离,最终在构造有利部位充填或交代形成金矿床,中基性岩浆则在矿体附近的构造薄弱带侵位,形成各种脉岩。
招平断裂带是胶西北地区重要的控矿构造带。该断裂带及其次级断裂(如谢家沟断裂)在中生代发育大量脉岩。区域研究表明,招平断裂带内的脉岩类型复杂,主要划分为基性、中性和酸性岩类。在时间上,这些脉岩呈现“脉动式”活动特征,在金成矿前、成矿期和成矿后均有发育(
杨冀,2010);在空间上,除部分脉岩本身构成矿体之外,大多数脉岩与金矿体相伴生,并且均受到断裂带的严格控制(
王建国等,2008)。谢家沟金矿成矿后石英闪长玢岩和花岗闪长岩的U-Pb年龄分别为(118.0±1.1)Ma和(118.4±0.9)Ma(
姚晓峰等,2020),与区域上招平断裂带中基性脉岩的主要活动时代(130~110 Ma)和胶西北地区金成矿高峰期(120 Ma左右)高度吻合(
王建国等,2008;
杨冀,2010)。这种时间上的相近性表明,脉岩活动与金成矿作用处于同一构造—岩浆事件框架内。
6 结论
(1)成矿前后岩脉均显示为高钾钙碱性、准铝质的I型花岗岩,具有富集LREE和LILE、亏损HREE和HFSE的特征,弱的负Eu异常和负Ce异常,表明来源于壳幔混合源区;成矿前岩浆活动复杂,有显著的幔源物质注入和强烈的壳幔相互作用,成矿后岩脉具有高度演化、成分均一的深部岩浆房末期特征。
(2)成矿后岩脉具有高Sr、低Y和无负Eu异常等埃达克岩特征,暗示其形成于洋壳俯冲背景下的高压环境,与古太平洋板块俯冲引发的壳幔作用密切相关。
(3)古太平洋板块俯冲导致地幔楔部分熔融,产生基性岩浆并向上运移,底侵古老下地壳产生酸性岩浆,两类成分差距较大的岩浆以不同比例发生岩浆混合作用,在成矿前产出中性蚀变岩和酸性细粒花岗岩。在主成矿期,充分混合的岩浆分异出含矿热液,在有利部位成矿,最终充分混合后形成成分均一的中酸性岩浆,在主成矿期后形成石英闪长玢岩和花岗闪长岩。
(4)起源于交代岩石圈地幔的基性岩浆富含水、硫和碳等挥发分,伴随着岩浆上侵演化过程,流体中Au元素逐渐富集,经历了与酸性岩浆的混合作用,物理化学条件发生改变而导致金在地壳浅部析出成矿。
新一轮找矿突破战略行动科技支撑项目“胶东牟乳金矿带构造控矿规律及找矿预测”(ZKKJ202405)
深地国家科技重大专项“重点金矿集区大数据集与三维勘探模型构建”(2024ZD1001902)
部省合作项目“胶西北金矿集区深部成矿机理研究与数智化找矿增储示范”(2024ZRBSHZ123)
泰山产业领军人才工程专项经费(tscx202507001)