西藏改则地区沙木罗组中早白垩世晚期闪长玢岩的发现：班公湖-怒江特提斯洋北向俯冲板片折返事件的响应 *

李发桥 唐菊兴 张静 宋扬 李海峰 林彬 王楠

1. 中国地质大学地球科学与资源学院，北京 100083 2. 中国地质科学院矿产资源研究所，自然资源部成矿作用与资源评价重点实验室，北京 100037

作为板块构造运动最重要的表现形式之一，大洋俯冲是研究古大洋演变和古老造山带形成的关键(Stern, 2002)。通常来说，大洋俯冲事件导致在俯冲带之上的活动大陆边缘形成由安山岩-英安岩-流纹岩系列和少量玄武岩等弧火山岩组成的大陆弧岩浆作用(郑永飞等, 2015)。而且大洋俯冲事件通常会经历长期的演变过程，从俯冲伊始到最终大洋的闭合，俯冲板片会经历初始俯冲、俯冲板片折返和断离、俯冲角度的变化等过程，与之响应的岩浆事件也会发生复杂的改变(DeCellesetal., 2009)。板片俯冲过程中发生脱水作用形成一系列弧岩浆作用(Syracuse and Abers, 2006)，加之岩浆弧在后期的地质作用中容易保留下来，因此其时空分布特点可以指示并复原俯冲过程中的运动学和动力学特征(李世民, 2018)。而对弧岩浆的类型等的判别可以指示其来源于俯冲沉积物(如Tatsumi, 2006)或俯冲洋壳(Defant and Drummond, 1990)等不同的源区特征。因此弧岩浆岩的研究对于还原俯冲带演化过程有着重要的意义。

班公湖-怒江缝合带(BNSZ，以下简称班-怒带)分布于青藏高原中部(图1a)，代表班公湖-怒江特提斯洋(以下简称班-怒洋)的残余遗迹(Yin and Harrison, 2000)。班-怒带西段一系列与中酸性岩浆岩组合有关的斑岩-浅成低温热液型铜多金属矿床的发现，例如形成于班-怒洋北向俯冲背景下的多龙矿集区等(唐菊兴等, 2016, 2017; Linetal., 2019; Songetal., 2018)，使得班-怒带成为近年来的研究热点。前人对于班-怒洋的演化研究已经进行了数十年，以沉积地层及古生物、岛弧火山岩、洋岛玄武岩、蛇绿混杂岩等研究对象开展了详细的讨论(Kappetal., 2007; Shietal., 2008; Yangetal., 2009; Lietal., 2009; Panetal., 2012; Zhuetal., 2013; Gengetal., 2016; 许志琴等, 2016)，但许多问题仍悬而未决。例如，学者们对于班-怒洋的开合时间、俯冲极性等有着不同的看法，认为班-怒洋初始裂解于石炭纪时期(Zhuetal., 2013)或二叠纪至早侏罗世、早白垩纪世(Wangetal., 2016; Fanetal., 2018; Lietal., 2019b)。班-怒洋的最终闭合可能发生于中-晚侏罗世(Maetal., 2017; Lietal., 2017a, 2019b)到晚白垩世(Kappetal., 2007; Zhangetal., 2012; Suietal., 2013; Fanetal., 2018)。传统观点认为班-怒洋的俯冲极性以北向俯冲为主(Yin and Harrison, 2000; Guynnetal., 2006; Liuetal., 2017;唐菊兴等, 2016; Zhangetal., 2017; Linetal., 2019; 宋扬等, 2019)，而近年来越来越多的学者们提出班-怒洋存在双向(向南、向北)的俯冲模式(Zhuetal., 2013, 2016; Lietal., 2018; Maetal., 2020)。以上争议势必造成了班-怒带两侧晚侏罗世-早白垩世岩浆作用成因机制的多解性，众多学者认为研究区北侧(图1)发育的晚侏罗世-早白垩世岩浆作用属于班-怒洋向北俯冲的产物(Kappetal., 2003; 唐菊兴等2016; Lietal., 2013, 2016; Dingetal., 2017; Linetal., 2017a, b)，其中多龙矿集区的成矿构造动力学背景亦存在着正常俯冲(Lietal., 2013)和俯冲板片折返(林彬等, 2019)两种模式。而位于多龙矿集区南侧改则县地区的闪长玢岩成岩时的俯冲阶段和深部动力学过程仍缺乏详细的刻画，这对该地区成矿作用的约束非常重要。本文对班-怒带西段北缘改则县地区上侏罗-下白垩统沙木罗组地层中的闪长玢岩开展了地质年代学、元素和同位素地球化学等方面的研究，对其岩石成因、构造地质背景及深部地球动力学过程等进行了探讨，旨在为班-怒洋北向俯冲深部过程与成矿作用的关系进行详细刻画，并为班-怒洋的演化提供新的证据。

图1 研究区大地构造位置及地质简图(a)青藏高原大地构造格架(据宋扬等, 2014; 唐菊兴等, 2017修改; 年龄数据源自雷传扬等, 2018; 刘文等, 2019及其中参考文献)；(b)藏西北大地构造格架及岩浆岩分布地质简图(据刘文等, 2019修改)；(c)研究区地质简图(底图据Huang et al., 2017修改). YZSZ-雅鲁藏布江结合带；BNSZ-班公湖-怒江结合带；JSSZ-金沙江结合带；KLSZ-昆仑结合带；LSSZ-龙木措-双湖结合带；SNMS-狮泉河-纳木错混杂岩带. 1-第四系；2-新近系康托组；3-古近系纳丁错组；4-上侏罗统-下白垩统沙木罗组；5-上三叠-下侏罗统木嘎岗日岩群；6-闪长玢岩；7-藏西北岩浆分布；8-藏西北沙木罗组分布；9-主缝合带；10-断层；11-地名；12-采样点Fig.1 Sketch tectonic map and geological map of the study area(a) sketch tectonic map of the Tibetan Plateau (after Song et al., 2014; Tang et al., 2017. Data of ages from Liu et al., 2019 and therein references); (b) sketch tectonic map of the North Tibet and the distribution of magmatic rocks (after Liu et al., 2019); (c) sketch geological map of study area (modified after Huang et al., 2017). YZSZ-Yarlung Zangbo Suture Zone; BNSZ-Bangongco-Nujiang Suture Zone; JSSZ-Jinsha Suture Zone; KLSZ-Kunlun Suture Zone; LSSZ-Longmu Co-Shuanghu Suture zone; SNMZ-Shiquan River-Namu Co Suture zone. 1-Quaternary; 2-Neogene Kangtuo Fm.; 3-Paleogene Nadingcuo Fm.; 4-Upper Jurassic to Lower Cretaceous Shamuluo Fm.; 5-Upper Triassic to Lower Jurassic Mugagangri Group; 6-diorite porphyrites; 7-the distribution of magmatism in Northwest Tibet; 8-the distribution of Shamuluo Fm. in Northwest Tibet; 9-main suture zone; 10-faults; 11-county; 12-locations of samples

1 区域及研究区地质概况

班-怒带作为西藏三大板块缝合带之一(图1a)，同样也是西藏四大成矿带之一(其余为冈底斯、北喜马拉雅和三江成矿带)(唐菊兴等, 2019)，属于南侧拉萨地体和北侧羌塘地体的分界线。班-怒带西端以班公湖为起点，向东经改则、丁青和嘉玉桥等到八宿，东西跨度大于2000km，以改则和丁青县为界分为西、中、东段(强巴扎西等, 2016)。该缝合带主要由侏罗纪的蛇绿岩和中生代的木嘎岗日群复理石沉积建造组成，其中蛇绿岩主要表现出与俯冲有关的SSZ型(Shietal., 2008; Wangetal., 2016)，木嘎岗日群主要发育有半深海-深海环境沉积的砂岩、页岩以及少量的玄武岩(Fanetal., 2015; Zengetal., 2016)。带内另发育晚侏罗-早白垩世的岩浆岩及沉积地层等。

图2 西藏改则地区闪长玢岩野外(a、b)、手标本(c、d)及镜下(e、f)照片Hb-角闪石；Chl-绿泥石；Cal-方解石；Pl-斜长石Fig.2 The field (a, b)，hand specimen (c, d) and microscope (e, f) photos of diorite porphyrite in Gaize area, TibetHb-hornblende; Chl-chlorite; Cal-calcite; Pl-plagioclase

研究区位于班-怒带内西段改则县西北约20km，大地构造位于羌塘地体南缘和拉萨地体北缘(图1)。区内发育地层主要包括上三叠统-下侏罗统木嘎岗日岩群(Τ3J1M)、上侏罗统-下白垩统沙木罗组(J3K1s)、古近系纳丁错组(En)、新近系康托组(Nk)以及第四系沉积物和残坡积物(Q)等(图1c)。木嘎岗日岩群主要是一套复理石建造，岩性主要包括粉砂质板岩和变质长石石英砂岩；沙木罗组主要含有一套弱变形的碎屑岩，含灰岩夹层，不整合覆盖于木嘎岗日岩群和蛇绿岩之上；纳丁错组为一套火山岩建造，岩性主要包括玄武岩、火山角砾岩、凝灰岩等，以断层方式逆冲推覆于南侧沙木罗地层之上(图1c)；康托组岩性主要为一套砂岩、砾岩，角度不整合覆盖于纳丁错组之上。研究区中酸性侵入岩发育较少，但研究区北侧以及多龙矿集区区域发育大量中生代岩浆活动(图1b)。区内NWW向断裂发育，为一系列北倾的逆冲断层，是班-怒洋盆北向俯冲留下的构造痕迹。

本次研究的闪长玢岩呈岩脉侵位于沙木罗组碎屑岩中，近东西走向，出露面积较小，与围岩有着明显的接触界线，野外观察无显著的烘烤边或冷凝边结构(图2a, b)。本文研究采集了新鲜的样品，手标本样品未发现明显的蚀变和风化现象，亦无脉体发育，呈灰黑色，斑状结构，块状构造，基质为显微晶质结构(图2c, d)。斑晶含量约占10%，杂乱分布，局部矿物发育一定程度的定向性，主要由角闪石(±8%)和斜长石(±2%)组成，其中角闪石可见发育明显的两阶段，早期呈深绿色，被后期淡绿色角闪石所穿切(图2e)，发育暗化边，呈自形程度较高的长柱状或六边形状，粒径大小不均，最小0.2mm，最大可达5mm以上，部分角闪石发生局部绿泥石化蚀变；斜长石多为椭圆状，表面较脏，部分可见环带发育，少数矿物发育方解石化和泥化蚀变。其中基质占比约90%，以斜长石为主，发育较少的石英矿物和角闪石微晶，并发育以锆石和磷灰石为主的副矿物，可见磷灰石以嵌晶结构发育于角闪石中。岩体边缘相发育围岩捕虏体，肉眼观察无明显反应边，而镜下可观察到捕虏体边缘发育薄层的暗色烘烤边，捕虏体由斜长石、石英等岩屑组成，多发生绿泥石、泥质等蚀变(图2f)。

图3 西藏改则地区闪长玢岩锆石U-Pb年龄谐和图(a)和年龄频率分布直方图(b)及117～126Ma段(c)、399～412Ma段(d)、728～762Ma段(e)年龄谐和图及加权平均年龄图Fig.3 The concordia plots (a) and the frequency distribution histogram (b) of all zircons with the concordia plots and weight mean ages of stages of 117～126Ma (c), 399～412Ma (d) and 728～762Ma (e) of zircons from diorite porphyrites in Gaize area, Tibet

图4 西藏改则地区闪长玢岩锆石的球粒陨石标准化稀土元素配分模式图(标准化值据Sun and McDonough, 1989；I、S、A型花岗岩范围据赵志丹等, 2018)Fig.4 Chondrite-normalized REEs patterns of zircons from diorite porphyrites in Gaize area, Tibet (normalized values from Sun and McDonough, 1989; the distributions of I-, S-, A-type granites are based on Zhao et al., 2018)

2 测试方法

2.1 锆石U-Pb年龄及微量元素

用于锆石测年的样品(18GZD04)采自地表闪长玢岩露头新鲜样品，采样位置如图1c所示。锆石单矿物分选及CL拍照在广州拓岩检测技术有限公司完成，将破碎至80～120目的样品粉末进行淘洗，以去除较轻的矿物，最后利用磁选和浮选两种方法分选出锆石矿物颗粒，挑选晶形和色泽较好、通透程度高的锆石颗粒进行环氧树脂制靶，用阴极发光(CL)显微镜照相观察其内部结构并拍照，选点测试时避开内部矿物包裹体、裂痕进行圈点。锆石激光剥蚀等离子体质谱(LA-ICP-MS)U-Pb同位素测年实验在中国地质科学院矿产资源所自然资源部成矿作用与资源评价重点实验室完成。使用RESOlution S-155激光剥蚀系统对矿物进行剥蚀，利用BRUKER auroa M90电感耦合等离子体质谱仪对剥蚀出的成分开展测试，载气为氦气，补充气体为氩气，在进入ICP-MS之前与氮气混合于T形连接器。使用直径为30μm的激光束斑进行剥蚀，频率为10Hz，测试之前进行0.5s的预剥蚀清洁锆石表面，每次分析采集约20s的背景值后再采集40s作为数据值。锆石年龄计算时的外标为国际标准锆石GJ-1和91500。用ICP-MS-DataCal程序计算测试结果(Liuetal., 2010)，未进行铅校正，利用Excel中的Isoplot程序完成平均年龄计算及谐和图绘制(Ludwig, 2003)。详细的分析步骤和数据处理过程可参考文献(侯可军和袁顺达, 2010)。

2.2 全岩主微量元素

样品采集位置与锆石U-Pb测年位置相同，样品新鲜，无蚀变和风化作用影响。在广州拓岩检测技术有限公司完成样品200目碎样工作，在核工业北京地质研究院分析测试中心完成全岩主微量元素测试。参照GB/T14506.28-93硅酸盐岩石化学分析方法X射线荧光光谱法，在飞利浦PW2404X射线荧光光谱仪上完成了主量元素的测试分析。采用DZ/T0223-2001电感耦合等离子体质谱方法在Finnigan MAT ELEMENT型高分辨率等离子体质谱仪(ICP-MS)上测定微量元素和稀土元素，大致过程为：每个样品称取50mg样品粉末置于PTFE溶样器皿中，分别加入1mL浓度38%的HF和0.5mL浓度68%的HNO3，蒸干溶液以除掉大部分的硅。再次加入与第一次浓度一致的1mL的HF和0.5mL的HNO3，放置在190℃烘箱中加热12小时。冷却后加入1mL浓度为0.5μg/mL的Rh内标溶液，加热至约150℃蒸干溶液，并再重复两次以上冷却步骤。用8mL HNO3(浓度40%)提取最终的残留物，重新密封溶样器，将其放入110℃烘箱中加热3h冷却之后，加入去离子水将溶液稀释至100mL。

2.3 全岩Sr-Nd-Pb同位素

全岩Sr-Nd-Pb同位素测试分析在核工业北京地质研究院分析测试研究中心完成，采用ISOPROBE-T热电离质谱仪开展测试，其中Rb-Sr条件为单带，M+，质量分馏用86Sr/88Sr=0.1194校正，标准测量结果为：NBS987=0.710250±7，实验室流程本底：Rb=2×10-10g，Sr=2×10-10g；Sm-Nd条件为三带，M+，质量分馏用146Nd/144Nd=0.7129校正，标准测量结果为：JMC143Nd/144Nd=0.512109±3。全流程本底Sm、Nd<50pg；Pb同位素测定用磷酸硅胶将样品置于铼带上，用静态接收方式测量Pb同位素比值。NBS981未校正结果：208Pb/206Pb=2.164940±15，207Pb/206Pb=0.914338±7，204Pb/206Pb=0.0591107±2，全流程本底Pb<100pg。

3 测试结果

3.1 锆石U-Pb年龄及微量元素

本次工作对闪长玢岩样品中晶型完好、环带清晰的50颗锆石进行了锆石U-Pb测年，结果如表1和图3所示。CL图像显示，锆石呈灰黑色，少数呈灰白色；大部分发育震荡环带结构，较少发育扇形结构，极少数发育有变质边结构(图3)。锆石形态可分为两类：第一类呈柱状的自形晶形态，粒径50～150μm，长宽比3:1～2:1；第二类为磨圆较好，呈次圆状-圆状发育，长宽比2:1～1:1，粒径为50～100μm。锆石的Th/U=0.01～2.81，均值0.71，远大于0.1，绝大部分具有岩浆成因锆石特征(吴元保和郑永飞, 2004)。

所有测试锆石年龄均未明显偏离谐和线(图3a)，表明了有效的锆石测年结果。放射成因的207Pb在年轻锆石积累较少，因此当锆石测试计算年龄<1000Ma时选用206Pb/238U年龄作为其结晶年龄，当>1000Ma时选用207Pb/235U年龄。测试结果显示，闪长玢岩锆石U-Pb年龄介于117～2488Ma之间，集中分布于5个年龄段(图3b)，分别为(1)：117～126Ma阶段(占比12%)，表现出明显的峰值120Ma，该阶段显示出较为谐和的年龄，加权平均年龄为120.4±3.2Ma(图3c)；(2)：214～302Ma阶段，该阶段无明显的峰值；(3)：399～412Ma阶段(占比16%)，峰值为406Ma，亦发育较为谐和的年龄，加权平均年龄为405.6±4.8Ma(图3d)；(4)：728～762Ma阶段(占比14%，峰值736Ma)，年龄较谐和，加权平均年龄为736±9.4Ma(图3e)；(5)：1736～2488Ma阶段无峰值。最新锆石年龄为117Ma，因此将闪长斑岩成岩年龄限制在早白垩世中期。

表1 西藏改则地区闪长玢岩(样品18GZD04)锆石U-Pb年龄测年结果

图5 西藏改则地区闪长玢岩(Na2O+K2O)-SiO2 (a)和K2O-SiO2 (b)图解(底图据Middlemost, 1994)Fig.5 The diagrams of (Na2O+K2O) vs. SiO2 (a) and K2O vs. SiO2 (b) of diorite porphyrites in Gaize area, Tibet (base maps after Middlemost, 1994)

图6 西藏改则地区闪长玢岩球粒陨石标准化稀土元素配分模式图 (a)和原始地幔标准化微量元素蛛网图 (b)(标准化值据Sun and McDonough, 1989)Fig.6 Chondrite-normalized REEs patterns (a) and primitive-normalized trace-elements spider diagrams (b) of the diorite porphyrites in Gaize area, Tibet (normalizing values from Sun and McDonough, 1989)

在进行锆石U-Pb测年的同时开展了锆石原位微量元素测试实验，分析结果如表2所示，整体上锆石稀土元素总含量变化范围较大，介于138×10-6～1972×10-6之间，平均值为672×10-6，∑LREE介于0.33×10-6～119.1×10-6之间，∑HREE变化范围为132.2×10-6～1965×10-6。锆石的球粒陨石标准化稀土元素配分模式表现出明显的岩浆锆石特征(图4)，即强烈的左倾趋势，显著富集重稀土和亏损轻稀土，且轻重稀土元素分馏明显，并具有明显的Ce正异常和不同程度的Eu负异常。各个年龄阶段锆石稀土元素特征较均一，但相互间表现出一定的差异，尤其闪长玢岩结晶阶段(图4a)锆石的稀土配分模式图与其他较老年龄的锆石(图4a-c)相比表现出明显的差异，即相对于其他年龄段(Eu/Eu*=0.02～0.71，平均值0.27)，117～126Ma年龄段(点序号19、23、42、44、47、48)的锆石Eu元素负异常程度较弱(Eu/Eu*=0.68～0.89，平均值0.78)。原始测试数据利用Geo-fO2_V0.2软件计算(Lietal., 2019a)，其中Zr元素未测试，取固定值49600×10-6(Lietal., 2019a)，原岩数据取7件闪长玢岩样品全岩地球化学数据的平均值。结果显示，锆石微量元素Ti地质温度介于638～908℃之间(均值767℃)，其中岩体结晶年龄段(117～126Ma)锆石Ti地质温度介于761～844℃之间(均值810℃)。Eu/Eu*值介于0.01～1.01之间，117～126Ma年龄段的Eu/Eu*=0.77～1.01；Ce/Ce*值介于1～1463之间，117～126Ma年龄段的Ce/Ce*=32.7～192.7。

表3 西藏改则地区闪长玢岩全岩主量(wt%)、微量(×10-6)分析结果

3.2 全岩主微量元素

全岩主微量元素测试分析结果见表3。改则地区闪长玢岩SiO2含量为58.35%～59.40%(平均值为58.79%)，显示中性岩类特征，以及较高的Al2O3(16.39%～16.67%)、全碱(Na2O+K2O=5.79%～8.12%)、FeOT(6.00%～6.41%)和CaO(4.40%～6.12%)，较低的MgO(2.15%～2.53%)和P2O5含量(0.31～0.32)。其K2O/Na2O<1(0.66～0.91)，呈富钠贫钾特征。Mg#介于0.39～0.41之间，均值0.40。A/CNK<1(0.86～0.90)，属于准铝质系列。(Na2O+K2O)-SiO2及K2O-SiO2图解(图5a, b)显示闪长玢岩样品落入闪长岩-二长岩区域，表现出高钾钙碱性特征。

稀土元素分布特征较为均一(图6)，稀土元素总含量介于189.5×10-6～202.2×10-6之间，平均值为195.4×10-6。LREE/HREE=8.30～9.04(均值8.82)，(La/Yb)N介于9.20～10.62(均值9.99)，轻重稀土分馏现象显著。图6a显示其属于富集轻稀土、亏损重稀土的右倾型稀土配分模式，重稀土趋势较为平坦，Eu负异常不明显或较弱(δEu=0.80～0.95)。微量元素蛛网图显示(图6b)，闪长玢岩富集大离子亲石元素Ba、Sr以及高场强元素Th、U、Ce、Zr等，亏损Rb、K、Pb、Nb、Ta、P、Ti等。

3.3 全岩Sr-Nd-Pb同位素

全岩Sr-Nd同位素和Pb同位素分析结果分别见表4和表5。87Sr/86Sr测试比值介于0.705864～0.706328之间，高于原始地幔现代值(0.7045, Depaolo and Wasserburg, 1979)，143Nd/144Nd比值为0.512336～0.512453，低于原始地幔现代值(0.512638, Jacobson and Wasserburg, 1980)。利用前文获得的闪长玢岩最年轻年龄段加权平均年龄120.4Ma计算Sr、Nd同位素初始比值，结果显示，(87Sr/86Sr)i=0.7053～0.7059，(143Nd/144Nd)i=0.51225～0.51236，fSm/Nd值介于-0.45～-0.40之间，Nd同位素二阶段模式年龄tDM2集中于1.1～1.3Ga，εNd(t)=-4.6～-2.4(均值-3.6)。闪长玢岩测试的206Pb/204Pb比值为18.666～18.830，207Pb/204Pb=15.612～15.644，208Pb/204Pb=38.844～39.047。同样的，经加权平均年龄120.4Ma计算后的(206Pb/204Pb)t、(207Pb/204Pb)t、(208Pb/204Pb)t值分别为18.535～18.568、15.605～15.636、38.628～38.744。

表5 西藏改则地区闪长玢岩Pb同位素组成

图7 西藏改则地区闪长玢岩Ce-SiO2 (a,底图据Collins et al., 1982)和A/MF-C/MF (b)图解(底图据Altherr et al., 2000)A/MF=Al2O3/(MgO+FeOT)，C/MF=CaO/(MgO+FeOT)Fig.7 The Ce vs. SiO2 (a, base map after Collins et al., 1982) and A/MF vs. C/MF (b, base map after Altherr et al., 2000) diagrams of the diorite porphyrites in Gaize area, Tibet

图8 西藏改则地区闪长玢岩Sr/Y-Y图解(底图据Defant and Drummond, 1990)Fig.8 The Sr/Y vs. Y diagram of diorite porphyrite in Gaize area, Tibet (base map after Defant and Drummond, 1990)

4 讨论

4.1 闪长玢岩成岩年龄

本文采集的闪长玢岩呈岩脉出露于上侏罗统-下白垩统沙木罗岩屑长石石英砂岩中，前人对改则地区沙木罗地层中发育的碎屑锆石研究发现，其沉积年龄大于130Ma(Huangetal., 2017; Lietal., 2017b)。而闪长玢岩挑选出的锆石显示出碎屑锆石特征，类型较为复杂(图3)，锆石U-Pb测试结果(图3a, b)显示出较大的年龄变化范围(117～2488Ma)，主要集中于5个阶段，即117～126Ma(峰值120Ma)、214～302Ma(无峰值)、399～412Ma(峰值406Ma)、728～762Ma(峰值736Ma)、1736～2488Ma(无峰值)，其中大于200Ma的年龄分布特征与沙木罗组中碎屑锆石年龄分布范围一致，加之岩体边缘相中发育的沙木罗组地层捕虏体现象(图2)，证实了这些捕虏锆石来源于沙木罗地层。而其所在最年轻阶段(117～126Ma)年龄较为谐和，加权平均年龄为120.4Ma±3.4Ma，明显小于改则地区沙木罗组地层最小沉积年龄130Ma(Huangetal., 2017; Lietal., 2017b)，因此可以代表闪长玢岩的结晶年龄，即改则地区沙木罗组地层中发育的闪长玢岩形成于早白垩世中期(120.4±3.2Ma)，前人首次在该组地层中发现的闪长玢岩成岩于103.5±1.4Ma(雷传扬等, 2018)，因此本文报道的闪长玢岩为首次在沙木罗地层中发现的早白垩世中期岩浆活动。

4.2 岩石成因及岩浆源区

图9 西藏改则地区闪长玢岩Th/Hf-Th (a)和(87Sr/86Sr)i-SiO2 (b)图解Fig.9 The diagrams of Th/Hf vs. Th (a) and (87Sr/86Sr)i vs. SiO2 (b) of diorite porphyrites in Gaize area, Tibet

西藏改则地区沙木罗组发育的闪长玢岩有着较低的A/CNK比值(0.86～0.90)和P2O5含量(0.31%～0.32%)，不含有堇青石、白云母等过铝质标型矿物，因此不属于S型花岗岩类别；其1000×Ga/Al=2.19～2.48(均值2.38)，小于2.6，因此可以排除A型花岗岩类别；CIPW计算所得分异指数(Differentiation Index)介于59.33～69.31之间(均值63.55)，从而排除高分异花岗岩类别。改则地区闪长玢岩含有大量I型花岗岩的标志矿物角闪石，测试样品全部落入Ce-SiO2图解中的I型花岗岩类别范围内，大部分均落入A/MF-C/MF图解中的基性岩部分熔融区域，与上文提到的前人首次在沙木罗组地层中发现的闪长玢岩(雷传扬等, 2018)属于相同的成因类别(图7)，即均为I型花岗质岩石类型。岩浆锆石的球粒陨石标准化稀土元素配分模式图同样佐证了本文闪长玢岩属于I型花岗岩质岩石类型(图3a)。

本文闪长玢岩高Y含量(28.4×10-6～30.8×10-6)，低La/Yb比值(12.83～14.81)，高(87Sr/86Sr)i值(0.7053～0.7059)，非常低的Cr(9.70×10-6～20.80×10-6)、Ni(6.05×10-6～8.70×10-6)含量，均与埃达克质岩地球化学特征有明显的差异。Sr/Y-Y图解显示闪长玢岩样品均落入经典岛弧型岩石范围内(图8)，暗示本文闪长玢岩为典型的岛弧型中性岩浆岩，与研究区北部约80km处多龙矿集区中广泛发育的弧岩浆岩类似(唐菊兴等, 2016)。前人对弧环境形成的中酸性岩浆岩开展了大量研究工作，并总结出三类成因机制：(1)俯冲板片的部分熔融(Martinetal., 2005; Defant and Drummond, 1990; Castillo, 2012)；(2)幔源玄武质岩浆的分离结晶(Annen and Sparks, 2002; Groveetal., 2003)；(3)幔源玄武质岩浆与壳源长英质岩浆的混合，及壳幔混合作用(Richards, 2003; Kempetal., 2007; Chiaradiaetal., 2011)。

图10 西藏改则地区闪长玢岩Sr-Nd-Pb同位素组成(底图据Lin et al., 2018)DMM-亏损型地幔；HIUH-异常高238U/204Pb型地幔；EMⅠ-富集Ⅰ型地幔；EMⅡ-富集Ⅱ型地幔(引自Zindler and Hart, 1986; Norman and Leeman, 1989).数据来源：雅鲁藏布MORB引自Hou et al., 2004；玉龙引自Jiang et al., 2006；冈底斯、藏北引自吴伟中等, 2013多龙矿集区数据引自林彬等, 2019；雄村数据引自Hou et al., 2015；安第斯弧岩浆引自Rabbia et al., 2017；全球平均俯冲沉积物引自Petford and Atherton, 1998Fig.10 The Sr-Nd-Pb isotopic compositions of diorite porphyrites in Gaize area, Tibet (base map after Lin et al., 2018)DMM-Depleted MORB Mantle; HIUH-High 238U/204Pb Mantle; EMⅠ-type Ⅰ enriched Mantle; EMⅡ-type Ⅱ enriched Mantle (after Zindler and Hart, 1986; Norman and Leeman, 1989); Data sources: Yarlung Zangbo MORB (Hou et al., 2004); Yulong (Jiang et al., 2006); Gangdese and North Tibet (Wu et al., 2013); Duolong (Lin et al., 2019); Xiongcun (Hou et al., 2015); Andean arc magma (Rabbia et al., 2017); Global Average subduction sediments (Petford and Atherton, 1998)

图11 西藏改则地区闪长玢岩Th/Ce-Th/Sm图解(a,底图据Martin, 1999)和Th/Yb-Sr/Nd图解(b,底图据Woodhead et al., 1998)Fig.11 The diagrams of Th/Ce vs. Th/Sm (a, base map after Martin, 1999) and Th/Yb vs. Sr/Nd (b, base map after Woodhead et al., 1998) of diorite porphyrites in Gaize area, Tibet

前人研究认为俯冲板片一般不会发生部分熔融，但在含水的洋壳俯冲过程中可以发生部分熔融，即年轻的大洋岩石圈发生低角度缓慢俯冲时，洋壳物质可以发生饱和水或脱水熔融，基性岩部分熔融形成埃达克岩(张泽明等, 2020)。但目前没有证据表明班-怒洋洋壳于早白垩世发生了低角度的俯冲作用(Zhuetal., 2019)，且本文闪长玢岩无埃达克岩特征(图8)，因此排除了俯冲板片部分熔融的成因机制。强亲岩浆元素(如Ta、Th、La、Ce)会在岩浆发生分离结晶作用过程中发生富集，而亲岩浆元素(如Zr、Hf、Sm)的丰度也几乎同步上升(Allègre and Minster, 1978)。因此若发生岩浆的分离结晶作用，强亲岩浆与亲岩浆元素的比值应为常数；而在部分熔融过程中，亲岩浆元素进入熔体的速率要小于强亲岩浆岩元素，因此在岩浆发生部分熔融作用时，强亲岩浆元素/亲岩浆元素应呈逐渐上升趋势。强亲岩浆元素(Th)-亲岩浆元素(Hf)比值关系图(图9a)显示，本文闪长玢岩成岩过程以部分熔融为主，可能涉及轻微的分离结晶作用。斜长石的分离结晶作用发生在岩浆锆石结晶前和结晶过程可能会导致锆石中Eu元素的负异常(Hoskin and Ireland, 2000)，闪长玢岩结晶年龄段(117～126Ma)的锆石δEu(Eu/Eu*)=0.77～1.01(均值0.88)，显示轻微或无负异常，表明岩浆结晶过程斜长石的分离结晶作用不明显。加之研究区附近无同时代大规模基性岩浆岩出露，且闪长玢岩稀土配分模式图未出现“四分组”模式效应(图6a)，分异指数(DI)较低(59.33～69.31)，综合以上判断，可以排除幔源玄武质岩浆的分离结晶成因机制。利用实验岩石学方法，前人证明了形成于玄武质岩石部分熔融作用的岩浆岩Mg#值往往低于0.4，有地幔组份参与形成的岩浆岩Mg#值则大于0.4(Rapp and Watson, 1995)。本文闪长玢岩Mg#=0.39～0.41(均值0.40)，表明玄武质岩浆部分熔融是其主要的成因，源区可能有极少部分地幔组份的加入。加之SiO2与(87Sr/86Sr)i无明显的正比例关系(图9b)，暗示地壳物质的混染作用不明显。Sr-Nd同位素组成特征(图10a)显示闪长玢岩样品落入古老下地壳、富集地幔和亏损地幔形成的三角区域，与多龙矿集区岩浆岩范围一致。Pb同位素组成显示样品落入EMⅡ与MORB演化线之间(图10b)，并位于EMⅡ范围附近。一般认为，EMⅡ型地幔与俯冲及陆壳再循环的壳幔相互作用有关(Sun and McDonough, 1989; Weaver, 1991)，即板片熔体或流体与富集地幔发生了混合作用，表明源区可能为俯冲作用下的新生玄武质下地壳。Th/Ce-Th/Sm图解显示(图11a)，样品均沿俯冲沉积物交代混合曲线分布，表明岩浆源区有少量沉积物熔体组份的参与(<5%)，小于晚期(103.5Ma)侵入于沙木罗组的闪长玢岩沉积物残余比例(6.5%～8%, 雷传扬等, 2018)。而Th/Yb-Sr/Nd图解中(图11b)显示本文闪长玢岩源区流体主要为板片来源，与晚期闪长玢岩的洋底沉积物来源明显不同(雷传扬等, 2018)。前人研究认为，因Nd可溶解在沉积物熔体中，造成εNd(t)值的降低(Hawkesworthetal., 1993; Pearce and Peate, 1995; Kesseletal., 2005)，本文闪长玢岩εNd(t)值介于-4.6～-2.4之间(均值-3.6)，证实了源区沉积物熔体的存在。

源区残留斜长石将会引起中酸性岩浆熔体中Eu元素的负异常(Rollinson, 1993)，本文样品δEu值介于0.80～0.95之间(均值0.88)，显示轻微负Eu异常(图6a)，表明源区无或残留极少的斜长石等富Eu矿物。研究表明，Nb、Ta元素亏损而Y不亏损可能表明岩浆源区残留较多的石榴子石或角闪石(Pearce and Norry, 1979; 侯增谦等, 2003)，本文样品具有Nb、Ta亏损和Y不亏损的特征，表明源区可能残留石榴子石或角闪石。根据HREE在石榴石和角闪石中分配系数的不同可对源区残留相进一步区分，呈现出倾斜的HREE配分模式，Y/Yb比值明显大于10，则残留相以石榴石为主；若具有较平坦的HREE配分模式，且Y/Yb比值接近于10，则以角闪石为主要残留相(Defant and Drummond, 1990; 高永丰等, 2003)。改则地区闪长玢岩有着较为平坦的HREE配分模式(图6a)，且Y/Yb值介于8.55～9.57之间(均值9.14)，表明其源区残留相以角闪石为主。综上所述，本文认为改则地区闪长玢岩形成于新生玄武质下地壳的部分熔融，源区残留相以角闪石为主，并有少量沉积物组份的参与(<5%)，岩浆演化过程中分离结晶和地壳混染作用不明显。

4.3 地球动力学背景及意义

利用地球化学微量元素判别岩浆岩成岩构造环境时会存在多解性，因此我们运用多种微量元素结合同位素、锆石微量元素判别图解综合对改则地区闪长玢岩的成岩构造环境进行解译。高场强微量元素具有稳定、不活泼的特性，受到后期风化、蚀变等地质作用的影响较小，能够一定程度上反映岩浆源区特征，Th/Yb-Ta/Yb和Th/Yb-Nb/Yb图解(图12a, b)均显示闪长玢岩属于大陆边缘岛弧构造背景。如上文所述，锆石矿物由于其封闭温度高、稳定性高、抗后期地质作用改造能力强等特性，可以较准确的反映岩浆岩成岩构造背景(赵志丹等, 2018)。前人通过大数据分析处理，在已知构造背景情况下，总结出利用锆石原位微量元素进行判别成岩构造背景的图解(Carleyetal., 2014; Grimesetal., 2015)，并得到了有效的应用。本文闪长玢岩锆石矿物微量元素U/Yb-Hf和U/Yb-Nb/Yb图解(图12c, d)显示，所测试样品均落入大陆弧范围内，暗示其大陆岛弧性质。相比全岩主微量元素，全岩Pb同位素性质稳定，受到外界环境影响较小(唐攀等, 2016)，Pb同位素判别图解显示所测试样品均落入造山带范围甚至造山带集中区域(图12e, f)，与多龙矿集区Pb同位素分布特征一致，暗示与多龙矿集区有着相同的成岩构造背景，即形成于俯冲造山背景下。综上所诉，本文利用高场强微量元素比值、锆石微量元素比值以及全岩Pb同位素构造环境判别图解综合认为，改则地区闪长玢岩形成于俯冲造山背景下大陆弧环境。

区域上来讲，本文闪长玢岩侵位于晚侏罗-早白垩世的浅海陆架和碳酸盐岩台地环境中沉积的沙木罗组地层中(图1b)(Chenetal., 2004; Xieetal., 2010)。Lietal. (2017b)通过对本文研究区沙木罗组的沉积记录、碎屑锆石开展物源分析等研究认为沙木罗组有两种可能的沉积模式：(1)形成于残余海盆，班-怒洋盆的俯冲已经结束，没有明显的地壳缩短但仍存在残余海盆，此时南羌塘地体高于拉萨地体，拉萨地体的碎屑无法越过未填充的前陆盆地，因此沙木罗组沉积物源主要来自羌塘地体；(2)形成于沟槽斜坡盆地，此时班-怒洋仍在北向俯冲，全部来源于南羌塘地体的沉积物不整合覆盖于俯冲杂岩(木嘎岗日岩群)和蛇绿岩之上。研究区北侧多龙斑岩-浅成低温热液型铜多金属矿集区内发育的早白垩世(128～106Ma)火山岩、中酸性侵入岩均形成于大洋板片折返作用下的陆缘弧环境(唐菊兴等, 2016; 林彬等, 2019; Linetal., 2019)，即表明早白垩世班-怒洋仍在北向俯冲。雷传扬等(2018)对革吉县北阿翁错地区沙木罗组地层中首次发现的闪长玢岩开展了地质年代学和元素地球化学研究，提出该岩体形成于班-怒洋北向俯冲于羌塘地体的构造背景下，即103.5Ma时班-怒洋西段仍未完全闭合。综上我们认为在120.4Ma时侵位于改则地区沙木罗组地层中的闪长玢岩成岩于班-怒洋北向俯冲造山背景下的大陆弧环境，由于俯冲板片发生折返引起。侵位于沙木罗组的120.4Ma与103.5Ma闪长玢岩均成岩于相同的构造环境中，但两者深部动力学背景却有所不同：革吉地区沙木罗组地层中侵位的闪长玢岩(103.5Ma)源区物质中的海洋沉积物组份明显多于改则地区闪长玢岩(120.4Ma)(图11a)，暗示着北向俯冲板片角度变大加速了软流圈地幔物质上涌速度，使得更多的俯冲板片上覆沉积物熔融从而导致了源区沉积物组份占比的增加，加之革吉地区闪长玢岩锆饱和温度(759～775℃，均值763℃，数据自雷传扬等，2018)明显高于改则地区闪长玢岩(726～742℃，均值732℃)的现象可证明以上推论，暗示120.4～103.5Ma时，班-怒洋板片俯冲角度继续增大，板片前端折返继续进行(图13)。

图13 班-怒带西段早白垩世成岩、成矿动力学机制模式图(据林彬等, 2019修改)Fig.13 Dynamic model in Early Cretaceous of diagenesis and mineralization of the western section of BNSZ (modified after Lin et al., 2019)

4.4 成矿地质意义

青藏高原发育的多龙矿集区(韦少港等, 2016)、驱龙矿床(Yangetal., 2015)以及朱诺矿床(Sunetal., 2018, 2021)等多个著名的斑岩型铜矿床的形成均与闪长质岩石有所关联，因此闪长质岩石作为一种与成矿密切相关的岩石类型，其含矿性研究有重要的地质意义。研究区北侧多龙矿集区含矿侵入岩与本文闪长玢岩具有同样的成岩时代、构造背景及相似位置的岩浆源区(林彬等, 2019)，两者宿命却天差地别，随着铁格隆南这一西藏首例千万吨级斑岩-浅成低温热液型矿床以及拿若等矿床的发现(方向等, 2014; 唐菊兴等, 2016)，前者一跃成为世界级斑岩-浅成低温热液型铜多金属矿集区，并被划为全国28个对国民经济具有重要价值的矿区之一(《全国矿产资源规划(2016—2020年)》)，而处于缝合带内北缘的改则县附近地区却尚未发现有较大潜力的成矿靶区。

除了成矿围岩条件、成矿后的保存等因素，本文认为改则地区闪长玢岩的时空分布和岩浆源区氧化还原性特征成为其未成矿的众多因素之二：(1)120.4Ma时，多龙矿集区位于俯冲板片前端，随着俯冲作用的进行，俯冲板片前端持续发生脱水作用和部分熔融，较轻的长英质等物质不断析出并以熔体的形式上侵，使得板片前端比重逐渐增大，加之班-怒洋南侧冷的高密度大洋板片发生裂解和断离，由于“跷跷板”原理，导致北侧大洋板片由于重力原因发生下沉和折返，地幔物质加速上涌，板片流体交代上覆地幔发生部分熔融形成玄武质岩浆；基性岩浆继而发生底侵作用至下地壳，发生MASH作用，其形成的中酸性岩浆岩携带金属矿物上侵至增生带中最终形成了多龙矿集区(图13;林彬等, 2019)。而位于班-怒带北缘的同时代发育的改则地区闪长玢岩近乎位于俯冲板片后端根部(图13)，其俯冲角度相对较小，板片及上覆沉积物熔融程度较弱，且源区幔源物质加入较少，形成的中酸性岩浆规模较小，因此无法携带巨量金属上侵至浅部并成矿。(2)大量研究表明，含矿岩体常具有较高的Ce4+/Ce3+比值特征，反映了含矿岩浆具有更高的氧逸度，而前人研究证实岩浆携带S的能力与氧逸度呈正相关关系，而岩浆中S的含量直接决定着其是否成矿(Paratetal., 2008; Botcharnikovetal., 2011)。Ce元素在锆石中含有Ce3+、Ce4+两种价态，由于后者半径更接近Zr4+，因此锆石表现为显著的正Ce异常，当熔体/流体更氧化时，Ce以4价为主，会替代更多的Zr4+，因此锆石的Ce4+/Ce3+比值可作为反应锆石结晶时熔体(或流体)氧逸度大小的标志(Trailetal., 2012; Ballardetal., 2002)。尽管存在部分锆石微量元素无法有效识别含矿与贫矿岩体的实例(Sunetal., 2021)，但前人对多龙矿集区含矿/不含矿岩体锆石微量元素进行了一系列研究，结果显示多龙矿集区含矿和贫矿岩体的锆石微量有着明显的差别，并且认为多龙矿集区含矿岩体形成于相对较高的氧逸度环境，可将Ce4+/Ce3+>300作为岩体成矿的地球化学标志(李云强等, 2020)，与智利北部斑岩型铜矿的含矿岩体的锆石Ce4+/Ce3+比值下限一致(Ce4+/Ce3+>300, Ballardetal., 2002)。而本文研究区临近多龙矿集区且具有相同的构造背景和成岩时代，因此判断该地区闪长玢岩锆石微量元素可能会成为判断含矿性的有效指标，其结晶阶段(117～126Ma)锆石微量元素Ce4+/Ce3+介于30.7～192.7之间(表2)，远低于300，暗示改则地区闪长玢岩源区氧逸度较低，未达到成矿的条件。

如上文所述，班-怒洋壳北向俯冲板片在120.4～103.5Ma时仍在持续发生折返作用，俯冲角度变大，使得形成于103.5Ma的闪长玢岩源区沉积物组份增加。主成矿期约为120Ma的多龙矿集区成矿动力学机制发育位置同样的向南位移(图13)。因此本文推测改则县地区120Ma至班-怒洋俯冲结束的时间内与俯冲背景有关的成矿事件应发育于多龙矿集区与改则县闪长玢岩发育位置之间。但是否成矿还应考虑如构造、围岩、后期保存等多种条件的影响。

5 结论

(1)改则地区侵位于沙木罗组地层中的闪长玢岩成岩于120.4Ma，为首次于沙木罗组发现该时代的岩浆活动。

(2)通过元素地球化学、锆石微量元素分布特征等判断本文闪长玢岩属于I型花岗岩类，具有明显的岛弧型岩石特征，结合其全岩Sr-Nd-Pb同位素组成特征分析得出其成岩于新生玄武质下地壳的部分熔融，源区以角闪石为主要残留相，并有少量沉积物组份的参与，岩浆演变过程中分离结晶和地壳混染作用不明显。

(3)改则地区闪长玢岩形成于由班-怒洋北向俯冲板片折返作用引起的大陆弧环境，而相较之下，革吉地区103.5Ma的闪长玢岩源区具有更多的沉积物组份和更高的成岩温度，暗示120.4～103.5Ma阶段班-怒洋板片俯冲角度继续增大，板片前端折返作用继续进行。

(4)在现有数据分析下，本文认为闪长玢岩未成矿的原因主要有两点：120.4Ma时该岩体成岩于班-怒洋俯冲板片后端根部，俯冲角度相对较小，板片及上覆沉积物熔融程度较弱，且源区幔源物质加入较少，可形成的中酸性岩浆规模较小，因此无法携带巨量金属上侵至浅部并成矿；闪长玢岩源区具有较低的氧逸度，源区无法为成矿提供足够的S。相较于革吉地区闪长玢岩深部动力学成岩背景，我们认为120Ma至俯冲作用结束的时间内与俯冲背景有关的成矿事件应发育于改则地区与多龙矿集区之间，但是否见矿还应考虑构造、围岩、后期保存等条件的影响。

致谢感谢中铝西藏金龙矿业股份有限公司为笔者的野外工作提供的大力支持；感谢王嘉星、张开江、王超等在野外采样、拍照等过程中的帮助；特别感谢两位审稿专家提出的富有建设性的宝贵意见和期刊编辑老师的修改。