藏北尼玛地区白垩纪岩浆岩对班公湖-怒江缝合带演化的制约*

王伟 翟庆国 胡培远 唐跃 朱志才 吴昊

自然资源部深部动力学重点实验室, 中国地质科学院地质研究所, 北京 100037

洋盆的消减闭合与陆-陆碰撞是古缝合带形成的重要阶段。这个过程往往伴随着多样性的岩浆活动，这些活动记录了洋陆转换的关键信息。通过对岩浆岩进行岩相学、年代学、地球化学和同位素研究，我们可以深入了解洋盆闭合和陆-陆碰撞过程。

1 区域地质概况

青藏高原主要位于中国西藏自治区内，地处特提斯构造域的东段。高原上分布多个近东西向的缝合带和板块。班公湖-怒江缝合带作为班公湖-怒江特提斯洋的残余，南北分割羌塘和拉萨板块(图1a)。它西起班公湖，向东经过改则、丁青和怒江等地，并沿东南方向延伸，全长超过2000km。该缝合带主要由蛇绿岩和复理石沉积组成。蛇绿岩主体形成于侏罗纪，多数表现出与俯冲相关的SSZ型特征(史仁灯，2007；Liuetal., 2016；Wangetal., 2016)。木嘎岗日群普遍出露在缝合带内，主要由砂岩、页岩和玄武岩组成，代表了班公湖-怒江特提斯洋在中生代时期半深海-深海沉积(Fanetal., 2015；Huangetal., 2015；Zengetal., 2016a)。前人研究表明，班公湖-怒江特提斯洋可能在二叠-侏罗纪打开，并在侏罗-白垩纪发生双向俯冲(Yin and Harrison, 2000；Zhuetal., 2011, 2016；Panetal., 2012；Metcalfe, 2013；Caoetal., 2016；Lietal., 2016；Liuetal., 2016；Wangetal., 2016)。缝合带及其两侧广泛发育侏罗-白垩纪岩浆岩，这些岩浆活动记录了班公湖-怒江特提斯洋消减闭合以及板块碰撞的过程，是重塑古缝合带演化的重要窗口(Suietal., 2013；Chenetal., 2014, 2015；Lietal., 2014；Wangetal., 2014, 2018, 2020；Wuetal., 2015, 2019a, b; Caoetal., 2016；Fanetal., 2016；Haoetal., 2016；Zengetal., 2016b；Zhuetal., 2016；Huetal., 2017；Liuetal., 2017; Yietal., 2018)。

研究区地处西藏尼玛县，大地构造位于羌塘板块南缘和拉萨板块北部(图1b，c)。区域内未见蛇绿岩露头。北部研究区主要出露地层为中下侏罗统木嘎岗日群(J1-2M)和上侏罗统索瓦组(J3s；图1b)。木嘎岗日群含有砂岩、灰岩以及枕状玄武岩。索瓦组主要由中薄层状灰岩组成，代表了近滨浅海环境。索瓦组与木嘎岗日群呈断层接触关系。南部研究区出露的地层有志留-石炭纪地层、下白垩统朗山组(K1l)和上白垩统竟柱山组(K2j；图1c)。志留-石炭纪地层包括中上志留统扎弄俄玛组(S2-3z)、下泥盆统达尔东组(D1d)、中上泥盆统查果拉玛组(D2-3c)和下石炭统永珠组(C1y)。朗山组地层中发育大量生物碎屑灰岩，灰岩中可见圆笠虫化石。竟柱山组与朗山组断层接触，代表一套磨拉石沉积，由砾岩、砂岩和火山岩组成。夹层英安岩锆石U-Pb年龄为91Ma，表明竟柱山组沉积时代为晚白垩世早期(Wangetal., 2014)。同时，北部和南部研究区广泛分布中酸性侵入体，这些岩石后期侵入到古老地层中，形成时代为早白垩世到晚白垩世(Kappetal., 2007；Wangetal., 2014；Yangetal., 2018, 2019)。

图1 青藏高原大地构造简图(a，据Zhai et al., 2016)、虾别错地区地质简图(b)和张乃错地区地质简图(c)引用年龄数据来自Kapp et al.(2007)，Yang et al.(2018)和Wang et al.(2014)Fig.1 Tectonic framework of the Tibet Plateau(a, modified after Zhai et al., 2016)，geological sketch maps of the Xiabie Co area(b)and the Zhangnai Co area(c)

图2 西藏虾别错和张乃错地区花岗岩与包体野外和显微照片Amp-角闪石；Bi-黑云母；Or-正长石；Pl-斜长石；Q-石英Fig.2 Field photographs and microphotographs of the granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet

本次研究的白垩纪岩浆岩以侵入岩为主，分别出露在尼玛县北部虾别错和南部张乃错附近。虾别错岩体以花岗岩为主，侵入到上侏罗统索瓦组地层(图1b)。岩石遭受风化剥蚀，多呈圆柱体(图2a)。花岗岩里发育暗色包体，两者呈截然型接触(图2c)。寄主岩石岩性为含黑云母二长花岗岩，矿物组成为黑云母(10%)、斜长石(25%)、正长石(30%)和石英(35%)。矿物粒径一般在0.1～0.5mm(图2b)。包体野外直径10～30cm，形态各异，多呈椭圆状；岩性为石英闪长岩，含有黑云母(5%)、石英(10%)、角闪石(15%)和斜长石(70%)。矿物粒径在0.5～3mm之间(图2d)。张乃错花岗岩侵入到志留-石炭纪地层中(图1c)，岩性为含黑云母二长花岗岩(图2e)，造岩矿物为黑云母(10%)、斜长石(25%～30%)、正长石(25%～30%)和石英(35%)。矿物粒径为0.5～3mm。镜下结果显示，所有岩石经历了不同程度的蚀变作用(绿泥石化和碳酸岩化)。

2 分析方法

2.1 锆石U-Pb分析

锆石分选在河北省区域地质调查研究所完成。锆石透、反射光照片和阴极发光图像采集在中国地质科学院地质研究所完成。锆石U-Pb定年在北京科荟测试技术有限公司完成。利用LA-ICP-MS系统对锆石开展U-Pb同位素分析，详细测试过程和仪器运行条件可参考侯可军等(2009)。激光剥蚀直径为25μm，能量为6.25J/cm2，脉冲频率为8Hz。锆石标准GJ-1作为外标进行同位素分馏校正。采用软件ICPMSDataCal 8.0进行U-Th-Pb同位素比值漂移、元素含量校正和U-Th-Pb同位素比值及年龄计算(Liuetal., 2010)。锆石U-Pb 年龄谐和图绘制和加权平均年龄计算均采用Isoplot 3.0完成(Ludwig, 2003)。

2.2 锆石Lu-Hf同位素分析

锆石Lu-Hf同位素测试在北京科荟测试技术有限公司完成。测试仪器为Neptune Plus多接收电感耦合等离子体质谱仪与213nm激光剥蚀系统构成的LA-MC-ICP-MS。Lu-Hf同位素分析点位与U-Pb测年点位重合或相邻。具体分析过程见Wuetal.(2006)。分析过程中采用的束斑直径为45μm，剥蚀能量为10～11J/cm2，脉冲频率为10Hz。标准锆石GJ-1作为外标，监测仪器运行稳定性。本次实验过程中获得的GJ-1标样测试值为0.282006±4(2SD,n=90)，这与推荐值在误差范围内一致(Moreletal., 2008)。

2.3 全岩地球化学分析

全岩主量和微量元素分析在北京科荟测试技术有限公司完成。主量元素分析采用熔片X荧光光谱法(XRF)，仪器为SHIMADZU公司生产的XRF-1800型X射线荧光光谱仪。样品烧失量(LOI)的测定采用马弗炉加热烧失法。微量元素分析利用等离子体质谱仪(ICP-MS)完成。测试过程中将样品粉末烘干，称取约40mg加入HF和HNO3溶解，再密封到高压釜中，放至烘箱内190℃加热48小时。蒸干后的样品加入HNO3再次密封到高压釜中，转移至烘箱内加热烘干。冷却后，用浓HNO3充分溶解样品，并送至ICP-MS测定微量元素。分析结果误差一般小于5%。

3 测试结果

3.1 锆石U-Pb年龄

本次研究对3件侵入岩样品(虾别错花岗岩(18T329)及包体(18T335)和张乃错花岗岩(18T357))进行了锆石U-Pb定年，测试分析结果见图3和表1。在阴极发光图像里，锆石多呈长柱状、半自形-自形晶体，发育岩浆振荡环带。锆石长度在100～200μm之间，长宽比1:1～3:1(图3)。分析结果显示，锆石中Th含量为214×10-6～4373×10-6，U含量为353×10-6～3498×10-6，Th/U比值为0.26～1.87。上述这些特征表明锆石均为典型的岩浆成因锆石(Hoskin and Schaltegger, 2003)。

虾别错花岗岩和包体具有相似的206Pb/238U加权平均年龄，分别为122±1Ma(18T329；n=21，MSWD=0.016)和121±1Ma(18T335;n=18, MSWD=0.045)(图3a，b)。张乃错花岗岩样品(18T357)获得锆石206Pb/238U加权平均年龄为97±1Ma(18T357；n=19, MSWD=0.035)(图3c)。上述结果表明这些岩体形成于早白垩世到晚白垩世(122～97Ma)。

本课程的总体设计思路是，通过对高职铁道工程技术专业培养目标，与企业合作，对企业中桥隧施工工作要求及桥隧检查和养护工作任务进行调研，并对其岗位职业能力进行分析，确定课程三大目标，安排课程内容，以项目引领，以任务驱动，构建课程体系。结合工作任务及过程，安排知识的讲授，突出课程讲授内容的实用性和针对性；以学生/学徒职业能力培养为核心，工学结合，“教学做”一体，精心组织教学；在考核设计时，建立与职业岗位能力培养相适应的考评体系，辅助以网络手段，与职业能力相结合，实行过程的考核评价。

3.2 锆石Lu-Hf同位素

对3件锆石样品开展Lu-Hf同位素分析，结果见表2。虾别错花岗岩和包体具有相似的锆石Hf同位素组成。其中花岗岩锆石176Hf/177Hf比值为0.282767～0.282902，εHf(t)值为+2.4～+7.0，对应二阶段tDM2模式年龄为731～1023Ma(图4)。包体锆石176Hf/177Hf值范围在0.282783～0.282843之间，εHf(t)值+3.0～+5.1，对应tDM2模式年龄为854～989Ma(图4)。张乃错花岗岩锆石176Hf/177Hf比值为0.282779～ 0.282884，εHf(t)值为+2.2～+6.0，二阶段tDM2模式年龄为779～1016Ma(图4)。

表1 西藏虾别错和张乃错花岗岩与包体锆石U-Pb定年结果

续表1

图3 西藏虾别错和张乃错花岗岩与包体锆石U-Pb年龄谐和图、代表性锆石阴极发光(CL)图像和加权平均年龄图锆石CL图像里实线和虚线圆圈分别代表锆石U-Pb测年位置和Hf同位素测试位置Fig.3 Zircon U-Pb concordia diagrams, representative CL images, and weighted average age diagrams of the granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet

表2 西藏虾别错和张乃错花岗岩与包体锆石Hf同位素分析结果

图4 西藏虾别错和张乃错花岗岩与包体锆石U-Pb年龄-εHf(t)值图解Fig.4 Plot of U-Pb ages vs. εHf(t) values for granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet

3.3 主量和微量元素

所有样品全岩地球化学测试结果见表3。测试样品具有较低的烧失量(0.29%～ 0.91%)，表明采集的样品新鲜。在下文讨论过程中，所有样品主量元素重新标准化为无烧失量百分比。

虾别错花岗岩具有较高SiO2(71.32%～72.96%)、Al2O3(14.02%～15.09%)，K2O(4.39%～4.81%)和Na2O(3.77%～4.32%)以及较低CaO(1.51%～1.77%)、MgO(0.46%～0.54%)和P2O5(0.07%～0.09%)含量(图5)。在TAS和K2O-SiO2图解里，所有样品点均落入花岗岩区域，并表现出高钾钙碱性特征(图6a，b)。样品A/CNK值为0.99～1.01，属于准铝质到弱过铝质系列(图6c)。原始地幔标准化蛛网图上，所有样品富集Rb、Th、U和Pb，亏损Ba、Nb、Sr和Ti等元素(图7a)。该花岗岩相对富集轻稀土((La/Yb)N=4.59～8.26)，并表现出明显的Eu负异常(Eu/Eu*=0.35～0.46)。

表3 西藏虾别错和张乃错花岗岩与包体全岩主量(wt%)和微量元素(×10-6)组成

续表3

续表3

图5 西藏虾别错和张乃错花岗岩与包体Harker图解黄色区域代表了研究区白垩纪岩浆岩，数据引自Wang et al. (2014)和Yang et al. (2018, 2019)，图6、图7和图8同Fig.5 The Harker diagrams of the granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet

图6 西藏虾别错和张乃错花岗岩与包体(Na2O+K2O)-SiO2(a，Le Bas et al., 1986)、K2O-SiO2(b，LeMaitre, 2002)和A/NK-A/CNK(c，Le Maitre, 1989)图解Fig.6 (Na2O+K2O) vs. SiO2(a, Le Bas et al., 1986), K2O vs. SiO2(b, LeMaitre, 2002), and A/NK vs. A/CNK(c, Le Maitre, 1989)diagrams of the granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet

图7 西藏虾别错和张乃错花岗岩与包体原始地幔标准化微量元素蛛网图(a、c)和球粒陨石标准化稀土元素配分图解(b、d)(标准化值来自Sun and McDonough，1989)Fig.7 Primitive-normalized multi-element spider diagrams (a, c) and chondrite-normalized REE diagrams(b, d)of the granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet(normalizing values from Sun and McDonough, 1989)

与虾别错寄主花岗岩相比，包体样品具有较低SiO2(62.80%～63.88%)和K2O(1.75%～2.71%)含量，较高Al2O3(15.80%～16.11%)、Na2O(4.66%～5.07%)、CaO(4.16%～4.69%)、MgO(2.55%～2.85%)和P2O5(0.17%～0.18%)含量(图5)。这些样品点落入闪长岩至花岗闪长岩区域中，表现出中钾钙碱性特征(图6a，b)。A/CNK值(0.84～ 0.90)偏低，属于准铝质系列(图6c)。样品的微量元素分布与寄主花岗岩相似，富集Rb、Th、U和Pb，亏损Ba、Nb、Sr和Ti等元素(图7a)。在球粒陨石标准化稀土配分曲线图里，这些包体同样表现出轻、重稀土分异((La/Yb)N= 4.45～10.23)和Eu负异常(Eu/Eu*= 0.29～0.64；图7b)。但值得注意的是，部分样品具有微弱Ce负异常(Ce/Ce*=0.70～0.94)。

张乃错花岗岩具有高含量SiO2(76.30%～78.18%)、Al2O3(12.15%～12.81%)、Na2O(3.73%～3.99%)、K2O(4.51%～4.95%)和CaO(4.16%～4.69%)。MgO(0.07% ～0.16%)和P2O5(0.01%～0.03%)含量普遍偏低(图5)。张乃错花岗岩表现出高钾钙碱性特征，属于弱过铝质系列(A/CNK=1.01～1.04；图6b，c)。在原始地幔标准化蛛网图中，样品富集Rh、Th、U和Pb，亏损Ba、Nb、Sr和Ti等元素(图7c)。稀土配分曲线图上，Eu具有明显负异常(Eu/Eu*=0.01～0.07)。同时，一些样品也表现出Ce异常(Ce/Ce*=0.74～1.41；图7d)。

4 讨论

4.1 区域白垩纪岩浆岩

班公湖-怒江缝合带及其两侧广泛分布早-晚白垩世岩浆岩，这些岩浆岩成为探索班公湖-怒江特提斯洋演化的重要线索。本次工作在尼玛县周边获得了花岗岩形成年龄为122～97Ma，指示岩浆活动发生在早白垩世至晚白垩世时期。同期岩浆作用在研究区也有报道，如虾别错和张乃错附近出露的花岗岩体，其锆石U-Pb年龄为118～100Ma(Kappetal., 2007；Yangetal., 2018, 2019)。这些花岗岩均属于高分异花岗岩。此外，Wangetal. (2014)在张乃错附近竟柱山组里报道了～91Ma安山岩和英安岩，这些火山岩具有高MgO (2.8%～5.9%)，并表现出埃达克岩亲缘性。在研究区西侧阿索地区也出露白垩纪岩浆岩，其中包括早白垩世高镁闪长岩(128～124Ma)、A2型花岗岩(117～115Ma)和晚白垩世埃达克质花岗闪长斑岩(89～88Ma；Liuetal., 2019；Luoetal., 2019；Wangetal., 2020)。总的来说，研究区及其周缘在白垩世发育成分多样的中酸性岩浆岩，这些岩浆活动为重塑班公湖-怒江特提斯洋俯冲至闭合过程提供重要资料。

4.2 岩石成因

4.2.1 虾别错花岗岩和包体

地球化学特征表明，虾别错花岗岩属于I型花岗岩：(1)SiO2和P2O5呈负相关关系(图5i；Wuetal., 2003)；(2)较低的A/CNK值(0.99～1.01；图6c；Chappell and White 1974)；(3)低Zr和Zr+Nb+Ce+Y含量以及FeOT/MgO比值(图8a；Whalenetal., 1987)。I型花岗岩一般可由幔源基性岩浆结晶分异、壳源物质部分熔融或壳幔岩浆混合形成(Turneretal., 1992；Chappell, 1999；Collins and Richards, 2008；Pankhurstetal., 2013)。区域内未广泛出露基性岩浆岩，因此结晶分异成岩的可能性不大。虾别错花岗岩锆石具有正的εHf(t)值(+2.4～+7.0；图4)，暗示其可能由新生地壳部分熔融形成。同时，我们也观察到εHf(t)值变化范围较大，并且花岗岩中发育有暗色包体，这些特征表明虾别错花岗岩可能是壳幔岩浆混合作用的产物。对包体进行深入分析，诸多证据表明岩浆混合作用在花岗岩和包体形成过程中发挥了重要作用。

已有研究表明花岗岩中暗色包体可来源于：(1)寄主花岗岩早期堆晶产物(Shellnuttetal., 2010)；(2)花岗岩源区难以熔融的残留物(Chappelletal., 1987)；(3)花岗质岩浆对围岩的捕掳体(Chappelletal., 2012)；(4)幔源镁铁质岩浆加入到花岗质岩浆后混合的产物(Vernon, 1984；Peruginietal., 2003；Feeleyetal., 2008)。在第一个模型里，包体通常具有堆晶结构，其矿物粒径与寄主岩相近(Didier and Barbarin, 1991)。然而虾别错石英闪长质包体不发育堆晶结构，包体矿物粒径(0.5～3mm)明显大于寄主花岗岩矿物粒径(0.1～0.5mm)。第二个模型中，包体具有变质岩或沉积岩的结构构造特征(Chappell and White, 1992)。但是虾别错包体形态一般为椭圆形，镜下矿物具有典型的岩浆岩结构，表明该包体并非源区残留物(图2c，d)。同时，在寄主花岗岩和包体接触边界未发现烘烤现象，说明包体并不是围岩的捕掳体。所以，虾别错包体最有可能是幔源镁铁质岩浆加入到花岗质岩浆后混合的产物。包体具有较高的Mg#值(51～55)、Cr(33×10-6～43×10-6)和Ni(18×10-6～25×10-6)含量，暗示其来源于地幔。寄主花岗岩和包体造岩矿物中均含有黑云母、斜长石和石英，并且包体中出现平衡矿物石英，很可能是岩浆混合的结果。在Harker图解里两者CaO、MgO、TiO2、Al2O3、MgO和P2O5等氧化物与SiO2呈线性关系(图5)，表明寄主花岗岩和包体之间发生一定程度的混合作用(杨高学等，2009)。相似的微量元素分布以及Hf同位素特征(寄主花岗岩εHf(t)值：+2.4～+7.0；包体εHf(t)值：+3.0～+5.1)也表明寄主花岗岩和包体之间可能发生过成分交换(图4、图7a，b)。在FeOT-MgO图解中，两者成分显现出混合趋势(图8b；Zorpietal., 1991)。综合以上特征，虾别错包体可能是幔源岩浆与花岗质熔体混和的产物。因此，岩浆混合作用在虾别错寄主花岗岩和包体形成过程中扮演主要角色。

图8 西藏虾别错和张乃错花岗岩与包体(FeOT/MgO)-(Zr+Nb+Ce+Y)(a，Whalen et al., 1987)、FeOT-MgO(b，Zorpi et al., 1991)、(La/Yb)N-La(c，Wu et al., 2003)、Nb/Ta-Nb(d，Ballouard et al., 2016)和Y/Ho-Zr/Hf(e，Bau, 1996)图解Fig.8 FeOT/MgO vs. Zr+Nb+Ce+Y(a，Whalen et al., 1987)，FeOT vs. MgO(b，Zorpi et al., 1991)，(La/Yb)N vs. La(c，Wu et al., 2003)，Nb/Ta vs. Nb(d，Ballouard et al., 2016)and Y/Ho vs. Zr/Hf(e，Bau, 1996)diagrams of the granites and enclaves in the Xiabie Co and Zhangnai Co areas, Tibet

4.2.2 张乃错花岗岩

4.2.3 稀土异常

在稀土配分曲线图里，虾别错包体和张乃错花岗岩均表现出不规则的稀土元素分布特征，如Ce正负异常(Ce/Ce*=0.70～1.41；图7b，d)。样品TE1,3值在0.92～1.05之间，低于稀土四分组效应的标准(1.1；Irber, 1999)。对于岩浆岩里出现稀土不规则分布特征，目前有以下多个解释：(1)地表风化作用(Takahashietal., 2002；Maetal., 2007)；(2)流体与熔体相互作用(Bau, 1996；Irber, 1999；Zhaoetal., 2002；Veksleretal., 2005)；(3)副矿物结晶分异(Zhao and Cooper, 1993；Pan and Breaks, 1997)；和(4)继承原岩的稀土异常(Neal and Taylor, 1989；Shimizuetal., 1992；Shaoetal., 2015；Bellotetal., 2018)。研究的样品均具有较低的烧失量(0.29%～0.91%)，未经历明显的风化作用。一般高度演化的岩石会受热液蚀变而出现稀土异常，这些异常归因于岩浆熔体与热液相互作用(Zhaoetal., 2002)。此时，高演化的岩浆熔体Nb/Ta值一般小于5(Ballouardetal., 2016)，并具有明显偏离球粒陨石的Y/Ho和Zr/Hf比值(Bau, 1996)。虾别错包体和张乃错花岗岩样品Nb/Ta值在7～11之间(图8d)。Y/Ho和Zr/Hf比值分别为24～26和21～40，样品均落入球粒陨石区域及附近(图8e)。以上特征说明虾别错包体和张乃错花岗岩形成于纯熔体系统，不存在热液和熔体相互作用的过程。高分异花岗岩稀土元素分布由副矿物主导(吴福元等，2007)，副矿物结晶分异可能会导致稀土异常。然而，Rayleigh分异模型表明副矿物的结晶分异过程会导致出现不合理的矿物组分含量，并且模型结果并未获得与之相对应的稀土异常(Bau, 1996；Irber, 1999)。此外，虾别错包体不属于高分异花岗岩(图8a)，但其也表现出稀土异常特征。因此，我们认为研究样品的稀土异常现象可能与副矿物结晶分异过程无关。已有研究显示一些物源(如风化古土壤和海相沉积物)具有不规则的REE分布特征，这些物质加入到源区参与岩浆作用会使形成的岩浆岩具有类似的REE分布((Holeetal., 1984；Neal and Taylor, 1989；Bellotetal., 2018；Yangetal., 2018)。古土壤和海相沉积物均具有明显的Ce正负异常，该特征在虾别错包体和张乃错花岗岩中也同样存在。研究样品虽具有正的εHf(t)值(+2.4～+7.0)，但其部分二阶段模式年龄相对较老(可至1016Ma)。因此，我们认为研究样品源区可能混入这些古老循环物质并继承了它们的稀土特征。

4.3 构造背景

班公湖-怒江缝合带及其两侧发育白垩纪岩浆岩，该期岩浆作用与班公湖-怒江特提斯洋南-北双向俯冲和之后拉萨-羌塘板块碰撞有关(Zhuetal., 2016)。然而，班公湖-怒江特提斯洋在早白垩世时期是否完全闭合依旧处于争论中。部分学者认为班公湖-怒江特提斯洋在早白垩世时期尚未完全闭合，现有资料证据如早白垩世蛇绿岩(～132Ma洞错蛇绿岩；鲍佩声等，2007)、洋岛(～116Ma仲岗洋岛和～108Ma塔仁本洋岛；朱弟成等，2006；Fanetal., 2014)、半深海-深海沉积(～118Ma 扎嘎组；Fanetal., 2015)和俯冲相关的岩浆作用(Wangetal., 2020)。另外学者根据区域地层接触关系、海相到陆相的沉积转变以及碰撞背景深熔的花岗岩等认为班公湖-怒江特提斯洋在晚侏罗-早白垩世时期发生闭合，羌塘和拉萨板块发生碰撞(汪明州和董得源，1984；Deweyetal., 1988；Yin and Harrison, 2000；Kappetal., 2007；Maetal., 2017)。

虾别错早白垩世花岗岩(～122Ma)大地构造位于羌塘板块南缘，可能形成于班公湖-怒江特提斯洋北向俯冲或闭合背景。花岗岩中出现暗色包体，表明幔源物质在成岩过程中作出重要贡献。已有研究成果显示，班公湖-怒江缝合带及其附近早白垩世岩浆岩普遍具有正的εHf(t)值，并且部分岩石中也发育暗色包体(Zhuetal., 2009, 2011；Haoetal., 2016)。这些岩浆活动被解释为是由幔源岩浆底侵作用导致的。同时，尼玛地区沉积地层由海相过渡到非海相发生在125～118Ma，表明此时班公湖-怒江特提斯洋经历闭合的过程(Kappetal., 2007)。在邻区阿索，俯冲背景下形成的高镁闪长岩(128～124Ma)和拉萨-羌塘板块碰撞背景下形成的A2型花岗岩(117～115Ma)也限定了班公湖-怒江特提斯洋的闭合时间为124 ～117Ma(Wangetal., 2020)。因此，我们认为虾别错早白垩世花岗岩(～122Ma)应形成于大洋闭合的过程(图9)。由于洋壳密度较大而发生下沉，最终引发地幔流和幔源基性岩浆的底侵(Zhuetal., 2016)。

张乃错晚白垩世花岗岩(～97Ma)地处拉萨地块北部，其形成应与班公湖-怒江特提斯洋南向俯冲或拉萨-羌塘板块碰撞相关。如前文所述，在早白垩世晚期，班公湖-怒江特提斯洋已经发生闭合，拉萨和羌塘板块发生碰撞。研究区分布的晚白垩世陆相磨拉石沉积(竟柱山组)也表明该时期大洋已经闭合，区域进入陆内演化阶段(Panetal., 2012；Wangetal., 2014；Wuetal., 2019b)。持续的碰撞挤压可导致地壳增厚并发生局部熔融(Chenetal., 2015；Sunetal., 2015)。该认识得到了研究区及邻区晚白垩世岩浆活动的响应，如新生地壳发生部分熔融形成了～97Ma张乃错花岗岩和增厚地壳来源的～90Ma埃达克岩(图9；Wangetal., 2014；Liuetal., 2019；Luoetal., 2019)。

图9 西藏尼玛地区白垩纪岩浆岩(122～97Ma)地球动力学背景示意图Fig.9 Schematic illustrations showing the geodynamic settings of the Cretaceous magmatic rocks (122～97Ma) in the Nyima area, Tibet

本文研究的花岗质岩浆岩约束了白垩纪时期班公湖-怒江缝合带的演化过程。在早白垩世中晚期，南-北向俯冲的班公湖-怒江特提斯洋发生闭合，此时洋壳岩石圈下沉引发幔源基性岩浆底侵。底侵作用不仅提供了充足的热量，也让幔源物质加入到岩浆房中参与成岩过程，导致虾别错花岗岩(～122Ma)发育暗色包体并具有较亏损的锆石Hf同位素组成。在晚白垩世早期，随着洋盆闭合以及拉萨-羌塘板块的持续碰撞，地壳逐渐增厚。岩浆活动以地壳部分熔融为主，形成张乃错花岗岩(～97Ma)和伴生的埃达克岩。因此，区域上122～97Ma岩浆活动是班公湖-怒江缝合带从洋盆闭合到南北两侧板块挤压碰撞的响应。

5 结论

(1)尼玛地区虾别错花岗岩和包体形成于早白垩世(122～121Ma)，张乃错花岗岩形成于晚白垩世(97Ma)。

(2)虾别错花岗岩和包体是壳幔物质混合作用的产物，形成于班公湖-怒江特提斯洋闭合过程；张乃错花岗岩来源于新生地壳，后期经历高度结晶分异作用，形成于洋盆闭合后拉萨与羌塘板块碰撞过程。

(3)班公湖-怒江特提斯洋在早白垩世时期发生闭合。尼玛白垩纪岩浆作用记录了缝合带从洋盆闭合到拉萨-羌塘板块挤压碰撞的演变过程。

致谢感谢项目组成员在野外和室内实验工作中给予了无私帮助；两位审稿人详细审阅了本文并提出宝贵的修改意见，在此表示感谢。