黑龙江小多宝山矽卡岩型Fe-Cu矿床成矿时代及构造背景——来自岩石地球化学、锆石年代学及Hf同位素证据*

杨贺 马万里 蔡文艳 王可勇, 2 孙丰月 王建华 周红英 寇柏燕

1. 吉林大学地球科学学院，长春 1300612. 自然资源部东北亚矿产资源评价重点实验室，长春 1300613. 黑龙江多宝山铜业有限公司，嫩江 1614164. 中国地质调查局天津地质调查中心，天津 3001705. 内蒙古自治区第六地质矿产勘查开发院，海拉尔 0210081.

中国东北地区地处中亚造山带东段(Wuetal., 2003)，自古生代以来，该区先后经历了古亚洲洋构造域、蒙古-鄂霍茨克构造域及滨太平洋构造域发展演化及叠加改造作用过程(Sengor and Natal’in, 1996; Chenetal., 2017; Wangetal., 2017; 唐杰，2016；李宇，2018)，致使区内构造-岩浆活动强烈，不同时期、不同构造背景环境热液成矿作用广泛，形成了区内类型多样的众多重要贵、有色金属矿床(武广，2006)。通过区内典型矿床的深入研究，可以深刻揭示区域不同构造演化阶段热液成矿作用特征，对更好的总结区域成矿作用规律有重要的理论意义

多宝山-三矿沟矿集区是中国东北地区著名的Cu(Mo)多金属成矿带(Chenetal., 2017;葛文春等, 2007a; 郝宇杰, 2015)。该成矿带经历了加里东期、海西-印支期和燕山期强烈多期次构造-岩浆活动及与之相关的成矿作用，造就了区内丰富的矿产资源(褚少雄等, 2012)，主要包括多宝山、铜山斑岩型Cu-Mo-Au多金属矿床，争光热液型金矿床以及三矿沟、小多宝山矽卡岩型Fe-Cu矿床(谭成印等, 2010)。小多宝山Fe-Cu矿床位于是三矿沟-多宝山矿集区内一处典型的矽卡岩型矿床。矿床发现于20世纪50年代，规模虽小，但具有品位高、开采成本低的特点。前人已从成矿流体(白令安, 2013; 白令安等, 2016)、矿化蚀变(白令安等, 2015)、成矿岩体地球化学特征(白令安, 2013; 郝宇杰, 2015)等方面进行了一定的研究，但对于小多宝山矿床形成的时代及构造背景认识尚存争议。

目前对小多宝山矿床成矿时代认识主要存在以下四种观点：早侏罗世(白令安等, 2012;白令安, 2013)、三叠纪(Haoetal., 2017)、二叠纪(赵元艺等, 2011)以及奥陶纪(谭成印等, 2010; 赵忠海等, 2012)。而区域复杂的构造演化致使小多宝山矿床的构造背景亦暂无定论，部分学者认为多宝山-三矿沟多金属成矿带中侏罗世的成矿事件主要是受蒙古鄂霍次克洋俯冲作用影响(Richards, 1999; Meng, 2003; 武广等, 2008; 杨祖龙等, 2009; 佘宏全等, 2012; Chenetal., 2017)，但是另一部分学者认为古太平洋构造域的作用也不可忽视(隋振民等, 2007; 葛文春等, 2007a; Zhangetal., 2008; 褚少雄等, 2012; 郝宇杰等, 2013; Shuetal., 2016)。

本文拟通过对小多宝山花岗闪长岩进行岩石地球化学、锆石U-Pb及Hf同位素等研究，查明小多宝山成岩成矿时代及其岩浆源区、岩石成因和成岩构造背景，对提高该矿床成矿理论认识水平，并为三矿沟-多宝山成矿带构造-岩浆演化与多金属成矿作用研究提供重要依据。

1 地质背景

1.1 区域地质

小多宝山矽卡岩型Fe-Cu矿床所属的多宝山-三矿沟Cu-Mo多金属成矿带地处兴安地块东南部，是我国大兴安岭北段东缘重要的矿产矿集区。兴安地块是中国东北重要的构造单元之一，其北缘被塔源-喜桂图缝合带截断，南部为贺根山-黑河缝合带(Wuetal., 2011; 徐备等, 2014; Liuetal., 2017)(图1a)，由前寒武纪变质基底组成，被古生代和中生代的花岗质岩石侵入(武广, 2006)。

多宝山-三矿沟成矿带区内主要出露地层为古生代奥陶系、志留系、泥盆系，中生代白垩系及新生代第四系。奥陶系地层自下而上分别为铜山组、多宝山组、裸河组、以及爱辉组，主要岩性为安山岩、英安岩、板岩以及火山碎屑岩夹大理岩。志留系地层以卧都河组、八十里小河组和黄花沟组为主，岩性包括砂岩、粉砂岩和板岩。泥盆系地层自下而上依次为泥鳅河组、腰桑南组和根里河组，以凝灰岩、板岩、砂岩和粉砂岩为该地层主要岩性。二叠系主要包括林西组砂岩、板岩以及花朵山组火山岩。矿区东南角内亦有少量白垩系龙江组火山岩和九峰山组陆相含煤地层出露(邓轲等, 2018)。晚古生代造山事件促使区内形成NE向及NW向断裂为主的构造格局， 其中NW向断裂被认为是主要的控矿构造； 大量燕山期花岗质岩石在区内西北部、东部及西南部发育，加里东期花岗质岩石仅出露于多宝山地区，而华力西期花岗质岩石在区内中部及西南部皆有出露(图1b)。

图1 中国东北大地构造略图(a,据Wuetal., 2011)、三矿沟-多宝山矿集区地质图(b,据邓轲等, 2018)及小多宝山Fe-Cu矿床地质图(c,据赵玉明, 1979(1)赵玉明. 1979. 黑龙江省嫩江县小多宝山铜矿点地质普查报告. 哈尔滨: 黑龙江省地质第四队五分队)

Fig.1 Tectonic map of Northeast China (a, after Wuetal., 2011), geological map of Sankuanggou-Duobaoshan metallogenic belt (b, after Dengetal., 2018) and geological map of the Xiaoduobaoshan Fe-Cu deposit (c)

1.2 矿床地质

矿区内出露的地层主要为早古生代海相碎屑岩、火山熔岩和碳酸盐岩等，包括中奥陶统多宝山组和早志留统黄花沟组。多宝山组为一套中性凝灰岩、中性熔岩、安山玢岩夹粉砂岩及条带状大理岩等浅海相火山熔岩和火山碎屑岩(Haoetal., 2015, 2017)，倾向220°～225°，倾角介于60°～80°之间。黄花沟组则属于砂泥岩互层的正常浅海相地层(白令安等, 2016)。从空间上看，多宝山组与小多宝山矿床的成矿作用密切相关。Zhaoetal. (2018)在充分总结前人研究成果基础上，确定多宝山组火山岩成岩年龄介于506～447Ma之间。研究区内岩浆岩出露较少，按照侵位先后顺序划分为闪长岩、安山玢岩和花岗闪长岩(白令安等, 2015)。燕山期花岗闪长岩主要出露于矿区中部，呈NW向展布。矿体产出于花岗闪长岩和多宝山组中性凝灰岩及大理岩接触带内，受接触带构造控制，亦呈NW向展布(图1c)；此外，花岗闪长岩岩体内不发育石英硫化物矿脉，但在接触带附近发育有广泛的内矽卡岩，在岩体深部则发育广泛的硅化、绢云母化以及绿泥石化，因此，可以认为，燕山期花岗闪长岩为该矿床成矿岩体(白令安等, 2016)。小多宝山矿床共圈定7个铜矿体，多呈扁豆状或不规则透镜状，厚度几到几十米，延伸100～400m，铜矿石平均品位约为0.46%，最高可达4.58%(白令安等, 2015)，铁矿石的全铁品位约35%～40%。

矿石矿物主要有镜铁矿(图2a)、黄铜矿(图2b-e)、黄铁矿(图2c、e)、方铅矿、闪锌矿(图2d)、磁铁矿(图2b、e)，另有少量的赤铁矿及铜蓝等表生氧化物。脉石矿物包括石榴子石(图2f)、透辉石、角闪石、绿帘石、绿泥石、石英、方解石等。根据矿物共生组合及矿脉穿切关系，共划分为两个成矿期：矽卡岩期和石英-硫化物期。矽卡岩期可分为三个成矿阶段：石榴子石-辉石干矽卡岩阶段；绿帘石-阳起石等含水硅酸盐湿矽卡岩阶段；磁铁矿-镜铁矿-石英氧化物阶段。石英硫化物期则分为两个成矿阶段：黄铁矿-黄铜矿-石英早期硫化物阶段和方铅矿-闪锌矿-方解石晚期硫化物阶段。

图2 小多宝山Fe-Cu矿床矿物共生组合及矿脉穿切关系照片

图3 小多宝山花岗闪长岩显微照片

2 样品采集与测试方法

本文研究样品采自与小多宝山Fe-Cu矿床成矿密切相关的花岗闪长岩岩体，取样位置见图1c。岩石表面呈灰白色，中细粒结构，块状构造。主要的造岩矿物为石英(10%～15%)、斜长石(～60%)、角闪石(～10%)和黑云母(～15%)(图3)。其中，石英呈他形粒状，粒径为0.3～2.2mm；斜长石呈板状、不规则粒状，粒径0.5～1.8mm；角闪石呈短柱状，粒径0.2～1.2mm；黑云母呈叶片状，粒径0.4～2.5mm。此外，还存在少量榍石、锆石、磁铁矿等副矿物。

锆石U-Pb定年测试工作在吉林大学测试实验中心完成。实验中尽量选择无包裹体无裂纹的锆石，以高纯He作为剥蚀物质载气，采用美国国家标准技术研究院研制的人工合成硅酸盐玻璃标准参考物质NIST SRM610进行仪器最佳化，锆石年龄以国际标样91500作为外标确保标准和样品的仪器条件一致。运用软件ICP-MS DataCal(Liuetal., 2008, 2010)对实验数据进行处理，用Andersen (2002)方法进行普通Pb校正，年龄的计算与协和图绘制均使用国际标准程序Isoplot(ver3.0)(Ludwig, 2003)得出，各点分析得出的同位素比值及年龄误差为1σ。本次测试锆石LA-ICP-MS U-Pb分析结果见表1。

主量、微量及稀土元素的测试选用研究区新鲜花岗闪长岩样品，粉碎研磨至200目以下，在吉林大学测试中心完成。主量元素分析采用X荧光光谱分析技术，精度优于5%；微量及稀土元素采用美国安捷伦科技公司Agilent 7500A 型电感耦合等离子体质谱分析测试，国际标样BHVO-2、BCR-2和国家标样GBW07103、GBW07104为参考对象，当含量>10×10-6时，分析精度优于5%，<10×10-6时则优于10%。本次全岩主、微量元素分析结果见表2。

在LA-ICP-MS锆石U-Pb定年的基础上，参照锆石阴极发光(CL)图像， 选择具有代表性的锆石进行微区Hf同位素测定工作。锆石Lu-Hf同位素原位分析测试在天津地质调查中心实验测试室完成，采用美国Thermo Fisher公司生产的NEPTUNE多接收器电感耦合等离子质谱仪和ESI公司生产的NEW WAVE 193nm FX ArF准分子激光器。实验选取国际标样GJ-1作为参考矿物。对仪器的操作条件和分析程序详见文献(耿建珍等, 2011)。GJ-1的平均176Hf/177Hf值为0.282001±15(2σ, n=29)，与报道数据误差范围内一致(0.282003±18, Gerdes and Zeh, 2006; 0.282013±4, Yuanetal., 2008; 0.282006±24, 耿建珍等, 2011)。Lu-Hf同位素分析结果见表3。

表1小多宝山花岗闪长岩锆石LA-ICP-MSU-Pb分析结果

Table 1 Zircon LA-ICP-MS U-Pb data for granodiorite in the Xiaoduobaoshan area

测点号ThU(×10-6)Th/U207Pb/206Pb207Pb/235U206Pb/238U207Pb/206Pb207Pb/235U206Pb/238U比值1σ比值1σ比值1σ年龄(Ma)1σ年龄(Ma)1σ年龄(Ma)1σGXDBS-011142940.390.050400.001340.191180.004340.027790.000622142417841774 GXDBS-022195070.430.047850.000870.183840.003790.027790.00034922717131772 GXDBS-031593800.420.050450.001180.192470.005000.027680.000462163117941763 GXDBS-043786100.620.049900.000900.190710.003730.027700.000331902417731762 GXDBS-052105020.420.049860.001000.189140.003670.027620.000331882417631762 GXDBS-061904240.450.049260.001280.188370.004850.027730.000391603517541762 GXDBS-071364090.330.049760.001070.188460.004600.027590.000551842617541753 GXDBS-082284460.510.049610.001070.189850.004360.027780.000421772717741773 GXDBS-092054710.440.050160.001570.192200.006330.027800.000482034517951773 GXDBS-102634910.530.049960.001020.190890.003740.027820.000401932217731773 GXDBS-112685470.490.049220.000860.188000.003440.027750.000341582217531762 GXDBS-121894340.440.050210.001290.190610.005950.027510.000542053817751753 GXDBS-13973260.300.050330.001490.191010.005740.027650.000522103717751763 GXDBS-142405030.480.049190.000850.188240.003460.027760.000311572317531772 GXDBS-151875120.360.049520.000930.188520.003760.027620.000301722717531762 GXDBS-162765300.520.050240.000930.191620.003940.027630.000322062717831762 GXDBS-172004680.430.049710.001470.191900.005980.028020.000471824217851783 GXDBS-181353660.370.048790.001170.187980.006220.027760.000521384417551773 GXDBS-192535460.460.049630.001270.188680.004270.027610.000361783017641762 GXDBS-203576540.550.049860.000770.191280.003580.027740.000301882417831762 GXDBS-211183750.310.048870.001250.188440.005160.027880.000311414317541772 GXDBS-221694330.390.049490.001520.188240.005840.027540.000351714917551752 GXDBS-23992800.350.049980.001080.190390.004030.027770.000401942517731772

表2小多宝山花岗闪长岩主量元素(wt%)、稀土和微量元素(×10-6)相关参数

Table 2 The relevant parameters of major elements (wt%) and trace elements (×10-6) of granodiorite in the Xiaoduobaoshan area

样品号GXDBS-1GXDBS-2GXDBS-3样品号GXDBS-1GXDBS-2GXDBS-3SiO266.2065.9465.74Th9.527.507.78TiO20.500.510.51Ta0.10.10.6Al2O315.8615.6615.66Hf3.804.104.30Fe2O35.295.325.27Nb6.506.106.30FeOT7.417.457.37Cs3.814.122.12MnO0.070.060.08Ga18.9518.9018.80MgO1.651.681.66La25.322.423.5CaO3.353.323.68Ce47.944.646.7Na2O4.354.224.44Pr5.255.095.34K2O2.792.912.72Nd18.919.019.6P2O50.160.160.17Sm3.303.613.61LOI0.650.710.65Eu0.950.931.01Total100.2299.7899.93Gd2.892.863.04Na2O+K2O7.147.137.16Tb0.420.420.45K2O/Na2O1.561.451.63Dy2.712.572.61A/CNK0.970.970.92Ho0.490.510.50Mg#43.443.743.7Er1.381.371.38Ba550600580Tm0.210.210.21Cr182122Yb1.411.441.48Sr473485494Lu0.230.240.23V737978Y13.913.814.2Zr135136152∑REE111.3105.3110.0Rb70.477.463.1(La/Yb)N12.8711.1611.39U3.703.002.80δEu0.920.850.91

注：A/CNK=Al/(Ca+Na+K)(原子数)；FeOT=FeO+0.8998×Fe2O3；(La/Yb)N=(La/0.310)/(Yb/0.209)，N为球粒陨石标准化后的值；LOI为烧失量

图4 小多宝山花岗闪长岩代表性锆石阴极发光图像

3 分析结果

3.1 锆石U-Pb定年

本次工作对矿区内与小多宝山矿床成矿关系密切的花岗闪长岩样品(GXDBS)进行锆石LA-ICP-MS U-Pb定年，选取的锆石均为自形-半自形粒状或他形粒状，核边结构明显且发育有良好的震荡生长环带(图4)，显示出岩浆锆石的基本特点。锆石的Th和U含量分别为97×10-6～378×10-6和280×10-6～654×10-6，Th/U比值介于0.30～0.62，与岩浆锆石属性一致。23个分析点曲线位置一致(图5a)，206Pb/238U年龄区间为175～178Ma， 加权平均年龄为176±1Ma(MSWD=0.10, n=23)(图5b)，由此表明花岗闪长岩形成于早侏罗世。

表3小多宝山花岗闪长岩锆石Lu-Hf同位素分析结果

Table 3 Lu-Hf isotopic compositions of zircons from granodiorite in the Xiaoduobaoshan area

测点号年龄(Ma)176Yb177Hf2σ176Lu177Hf2σ176Hf177Hf2σεHf(0)εHf(t)tDM1(Ma)tDM2(Ma)fLu/HfXDBS-01XDBS-02XDBS-03XDBS-04XDBS-05XDBS-06XDBS-07XDBS-08XDBS-09XDBS-101760.0395160.0001840.0007620.0000010.2829050.0000224.78.5488678-0.980.0644680.0002040.0014070.0000080.2828930.0000254.38.0515711-0.960.0557980.0002540.0011990.0000030.2829820.0000257.411.2384506-0.960.0772050.0002230.0016150.0000030.2829900.0000277.711.4378492-0.950.0772970.0001820.0015710.0000040.2828950.0000244.38.0514707-0.950.0618020.0010090.0013190.0000200.2829050.0000234.78.4496683-0.960.0670470.0005340.0014380.0000080.2829290.0000235.59.2464630-0.960.0716770.0005610.0014370.0000100.2829770.0000287.210.9395521-0.960.0505820.0007360.0010510.0000130.2828930.0000244.38.0509707-0.970.0543610.0003440.0012050.0000060.2828830.0000293.97.6526732-0.96

图5 小多宝山花岗闪长岩锆石年龄协和图(a)和加权平均年龄图(b)

图6 小多宝山、三矿沟花岗闪长岩TAS图解(底图据Irvine and Baragar, 1971)

3.2 全岩地球化学特征

3.2.1 主量元素地球化学

研究区花岗闪长岩样品SiO2含量65.74%～66.20%，属酸性岩。Na2O含量为4.22%～4.44%，K2O含量为2.72%～2.91%，全碱含量(Na2O+K2O)范围为7.13%～7.14%，Na2O/K2O值介于1.45%～1.63%，表现出岩石相对富钠的特征。Al2O3含量为15.66%～15.86%，CaO含量为3.32%～3.68%，MgO含量为1.65%～1.68%，显示出岩石相对富铝贫钙、镁的特征。借助火成岩TAS分类图解(图6)，样品全部落入花岗闪长岩区域。A/CNK变化范围在0.92～0.97之间，铝饱和指数图解显示其为准铝质系列(图7a)。里特曼指数σ范围为2.20～2.25，属于钙碱性系列。在SiO2-K2O图解上，岩石落入钙碱性系列与高钾-钙碱性系列的过渡区域，且多数偏向高钾-钙碱性岩石系列，因此可进一步确定小多宝山花岗闪长岩为准铝质-高钾钙碱性系列(图7b)。

3.2.2 稀土和微量元素地球化学

花岗闪长岩稀土元素总含量相对较低(∑REE=105.3×10-6～111.3×10-6)，相对富集轻稀土元素(LREE/HREE=9.94～10.43)，轻、重稀土分馏程度中等[(La/Yb)N=11.16～12.87]，表现出弱的负Eu异常(δEu=0.85～0.92)。从球粒陨石标准化稀土元素配分曲线(图8a)可明显看出，花岗闪长岩样品整体呈明显右倾式，表现出轻稀土元素明显富集，重稀土元素相对亏损的形势。原始地幔标准化蛛网图显示(图8b)，岩石相对富集Rb、Ba、Sr、K等大离子亲石元素，而相对亏损Ta、Nb、Ti等高场强元素。

3.3 锆石Lu-Hf同位素

选取花岗闪长岩样品(GXDBS)中10个能代表小多宝山成矿年龄的锆石，圈定其相同或近似区域进行Lu-Hf同位素测试，结果显示，176Lu/177Hf平均比值为0.0013，表明锆石形成后由Lu衰变为Hf的数量极少，因而所得的176Hf/177Hf值可以反映岩石结晶演化过程中Hf同位素的组成情况(吴福元等, 2007)。另外，锆石样品的fLu/Hf普遍较低，平均值为-0.961，说明实验锆石Lu-Hf同位素二阶段模式年龄(tDM2)指示源区物质在地壳中存留年限或从地幔中抽离的年限是合理的。10个点的Hf同位素分析结果表明，176Hf/177Hf值为0.282883～0.282990，通过其结晶年龄176Ma计算可得，εHf(t)为7.6～11.4，单阶段模式年龄(tDM1)为378～526Ma，二阶段模式年龄(tDM2)为492～732Ma。

4 讨论

4.1 成岩成矿时代

区域构造活动频繁，先后经历了多期构造运动，进而导致区内岩浆活动强烈，衍生了一系列内生热液金属矿床。综合本次研究及前人研究成果，多宝山-三矿沟矿集区主要矿床成岩成矿时代可分为以下三期：奥陶纪斑岩型Cu(Mo)矿化(Zhaoetal., 2018; 佘宏全等, 2012; 向安平等, 2012; 崔根等, 2008; 葛文春等, 2007b; Zengetal., 2014; 白令安, 2013; 赵焕利等, 2012; Wuetal., 2015; 郝宇杰, 2015; 赵一鸣等, 1997; Liuetal., 2012)；三叠纪斑岩型Cu(Mo)-Au矿化(Haoetal., 2017; 郝宇杰, 2015; Zengetal., 2014; 赵元艺等, 2011; 杜琦等, 1988)；早侏罗世矽卡岩型Fe-Cu矿化(Chuetal., 2019; Dengetal., 2018; Haoetal., 2015; 郝宇杰, 2015; 吕鹏瑞等, 2012; 褚少雄等, 2012; 李德荣，2011；李德荣等, 2010, 2011; 葛文春等, 2007a)。

小多宝山Fe-Cu矿体赋存于花岗闪长岩和多宝山组中性凝灰岩及大理岩接触带内，且在花岗闪长岩内广泛出露的由石榴子石、辉石、绿帘石组成的内矽卡岩，同时并未在花岗闪长岩体内发现石英硫化物矿脉，因此可以推断花岗闪长岩为小多宝山Fe-Cu矿的成矿岩体，其成岩年龄可以代表矿床的成矿年龄。本次工作得出小多宝山花岗闪长岩锆石U-Pb年龄为176±1Ma，白令安(2013)报道过小多宝山花岗闪长岩锆石LA-ICP-MS加权平均年龄为171.9±1.7Ma(MSWD=11.4, N=17)，表明该矿床成矿与燕山期早侏罗世岩浆作用密切相关。三矿沟Fe-Cu矿床作为三矿沟-多宝山成矿带上矽卡岩型矿床的典型代表，在空间位置上与小多宝山矿床十分接近，另外，褚少雄等(2012)和Dengetal. (2018) 分别对三矿沟矿床花岗闪长岩的年龄进行测定，加权平均年龄分别为175.9±1.1Ma(MSWD=0.21, N=24)和177±1Ma(MSWD=1.01, N=26)。相较于171.9±1.7Ma，本文对小多宝山年代学的测定(176±1Ma)更接近矿床实际成矿时代。

4.2 岩浆源区特征与岩石成因

小多宝山花岗闪长岩主要矿物组成为石英、斜长石、角闪石和黑云母，符合I型花岗岩矿物组合的基本特征。在主量元素上表现出富硅(SiO2=65.74%～66.20%)、准铝(A/CNK=0.92～0.97)、中等镁值(Mg#=43.4～43.7)及较高的K2O含量的特征；在微量元素上表现出富集轻稀土(LREE)元素、大离子亲石元素(LILE: Rb、Ba、K)和地球化学性质活泼的不相容元素(U、Th、Pb)，相对亏损高场强元素(HFSE: Nb、Ta、Ti)的特征。此外，Zr+Nd+Ce+Y的值介于200.5×10-6～219.2×10-6，低于A型花岗岩的下限值(350×10-6)(Whalenetal., 1987)，在 (Na2O+K2O)/CaO-(Zr+Nb+Ce+Y) 和Ce-SiO2判别图解中，数据点均落入I型花岗岩区域(图9a, b)。综合判断，小多宝山花岗闪长岩应属于准铝质高钾钙碱性I型花岗岩。

图7 小多宝山、三矿沟花岗闪长岩A/CNK-A/NK图解(a,底图据Maniar and Piccoli, 1989)和SiO2-K2O图解(b,底图据Peccerillo and Taylor, 1976)

图8 小多宝山、三矿沟花岗闪长岩球粒陨石标准化稀土元素配分模式图(a)和原始地幔标准化微量元素蛛网图(b)(标准化值据Sun and McDonough, 1989)

图9 小多宝山、三矿沟花岗闪长岩(K2O+Na2O)/CaO-(Zr+Nb+Ce+Y)判别图解(a)和Ce-SiO2判别图解(b)(底图据Whalen et al., 1987)

图10 小多宝山、三矿沟花岗闪长岩锆石Hf同位素特征(底图据吴福元等, 2007; Wu et al., 2015)

普遍认为，花岗岩类岩石主要起源于：(1)地壳物质的部分熔融作用(葛文春等, 2007b; Hofmann, 1988)；(2)玄武质岩浆或安山质岩浆的分离结晶作用(Hanetal., 1997)；(3)酸性岩浆同玄武质岩浆的混合作用(Jahnetal., 2000; Yangetal., 2015)。Nb、Ta和Ti的负异常、弱的负Eu异常，Rb、Ba、K和Sr的正异常(图8b)及LREEs和LILEs的相对富集均表明岩浆演化过程中发生了显著的结晶分异作用，并存在地壳物质的混染。小多宝山花岗闪长岩较高的SiO2和Al2O3(15.66%～15.86%)含量及较低的Ni(6×10-6～6.6×10-6)和Cr(18×10-6～22×10-6)也暗示了地壳物质的参与。此外，花岗闪长岩中可见闪长岩包体，且闪长岩包体的地球化学特征同花岗闪长岩表现出良好的谐和关系(白令安, 2013)亦暗示两者为同源岩浆分异作用的产物。而白令安(2013)提供的Sr-Nd同位素特征也表明亲地幔端元和少量地壳物质的混入。本次Hf同位素测试结果表明，εHf(t)=7.6～11.4，一阶段模式年龄(tDM1)为378～526Ma，二阶段模式年龄(tDM2)为492～732Ma，Hf同位素图解中(图10a, b)，所有数据点均落于球粒陨石与亏损地幔之间，进一步证明小多宝山花岗闪长岩岩浆源区为古生代新生地壳物质，即洋壳的熔融。

4.3 成矿构造背景

包括研究区在内的中国东北地区先后受到古亚洲洋、蒙古鄂霍茨克洋和古太平洋的俯冲作用影响，围绕其构造单元划分及演化历史，前人已开展过大量研究工作，综合起来可将其演化过程概括为：(1)早古生代期间，古亚洲洋向松辽地块及华北地块北缘发生双向俯冲作用(石玉若等, 2005; Shietal., 2010; Liuetal., 2003)；(2)至石炭纪，兴安地块与松辽地块拼合，并伴随大规模的岩浆活动(Chenetal., 2009; Liuetal., 2009; Zhangetal., 2007, 2009; 包志伟等, 1994; Zhouetal., 2015; 施光海等, 2004)；(3)晚二叠世-早三叠世，松辽地块与华北地块北缘拼贴，古亚洲洋完成最终闭合，与此同时蒙古鄂霍茨克洋构造域开始向额尔古纳地块及西伯利亚板块俯冲(唐杰, 2016)；(4)至早-中侏罗世开始，东部的古太平洋构造域持续俯冲作用在欧亚板块之下(Xuetal., 2013)，至晚侏罗世-早白垩世，蒙古鄂霍茨克洋自西向东呈剪式闭合(Tomurtogooetal., 2005)。

本次研究结果显示与Fe-Cu矿化有关的花岗闪长岩锆石U-Pb年龄为176±1Ma，这与古太平洋俯冲作用时间相吻合，同时蒙古鄂霍茨克洋仍处于俯冲环境。地球化学特征表明小多宝山花岗闪长岩属准铝质高钾钙碱性系列I型花岗岩，富集LREEs和LILEs(如Rb、Ba、Sr、K)，亏损HREEs和HFSEs(如Nb、Ta、Ti)。Rb-Y+Nb、Nb-Y、Ta-Yb和Rb-Ta+Yb图解中(图11a-d)，花岗闪长岩样品点均落在火山弧花岗岩(VAG)区域内。以上研究结果与同期成矿事件的三矿沟花岗闪长岩(175.9±1.1Ma, 褚少雄等, 2012；177±1Ma, Dengetal., 2018)相比，两者在地球化学特征上表现出高度相似性(图6-图11)。

图11 小多宝山、三矿沟花岗闪长岩Rb-(Y+Nb) (a)、Nb-Y (b)、Ta-Yb (c)和Rb-(Ta+Yb) (d)图解(底图据Batchelor and Bowden, 1985；Pearce et al., 1984)

本次研究认为古太平洋和蒙古鄂霍茨克洋在中侏罗纪时期都处于俯冲构造环境。但Wuetal. (2011)和Xuetal. (2013)认为与太平洋俯冲缝合带相平行的小兴安岭-张广才地区-吉林中部-延边地区的SN(190～169Ma)向侵入岩带(Zhangetal., 2004; Wuetal., 2007, 2011; 孙德有等, 2001, 2005; 苗来成等, 2003; 隋振民等, 2007)受古太平洋板块俯冲的影响，而同期的额尔古纳地块区域沿蒙古鄂霍茨克缝合带NE-SW向分布的侵入岩则往往受蒙古鄂霍茨克洋俯冲的影响。兴安地块的侵入岩分布形态不能支撑研究区受到了蒙古鄂霍茨克洋俯冲的显著影响的结论，此外在小兴安岭-张广才地区-吉林中部-延边地区侵入岩带的西侧存在一与之平行的145～120Ma的岩浆岩带，自西向东侵入岩呈逐渐年轻的变化趋势，也暗示了古太平洋板块的持续俯冲(Chuetal., 2019)。值得注意的是，Zhouetal. (2009)对黑龙江群中超高压变质带中的蓝片岩(165～180Ma)的变质作用年龄的研究也显示佳木斯地块与松嫩地块的拼合动力主要源于古太平洋西向俯冲增生作用；另外，Guoetal. (2015)认为图们地区与俯冲相关的早侏罗世镁铁质侵入杂岩也是源于古太平洋板块俯冲作用，这些都暗示了古太平洋板块对东亚大陆边缘的持续俯冲作用。Shuetal. (2016)总结了中国东北地区140～200Ma期间形成的45个矿床的时空分布，发现了一个向西北方向年龄逐步降低的趋势，由此认为侏罗纪时期岩浆热液作用是由古太平洋板块的平板俯冲引起的。因此，形成于火山弧环境的小多宝山矿床及其他同期矿床主要受早中侏罗世古太平洋俯冲作用的影响。

5 结论

(1)小多宝山矽卡岩型Fe-Cu矿床成矿岩体锆石U-Pb加权平均年龄为176±1Ma，成矿时代属早侏罗世；

(2)成矿花岗闪长岩属准铝质高钾-钙碱性I型花岗岩，起源于幔源岩浆结晶分异作用，并伴有地壳物质的混染；其εHf(t)为7.6～11.4，暗示岩浆源区为古生代新生地壳的熔融；

(3)矿床形成于早侏罗世古太平洋板块俯冲形成的火山弧构造环境。

致谢感谢吉林大学测试实验中心对本次锆石U-Pb测年过程中提供的帮助。感谢黑龙江多宝山铜业有限公司褚向辉、王宇晨对野外工作的大力支持。感谢项目组成员对室内工作的辅助支持。