东昆仑五龙沟地区早泥盆世双峰式侵入岩年代学、地球化学及其地质意义*

2021-09-02张亮李碧乐刘磊王盘喜李良

岩石学报 2021年7期

张亮李碧乐刘磊王盘喜李良

1.中国地质科学院郑州矿产综合利用研究所，郑州 450006 2.国家非金属矿资源综合利用工程技术研究中心，郑州 450006 3.吉林大学地球科学学院，长春 130061 4.云南大学地球科学学院，昆明 650500

东昆仑造山带位于青藏高原的东北部，是中央造山带的重要组成部分(许志琴等，2006；杨经绥等，2010)。自中元古代以来，该区经历了多期洋-陆转换、弧-弧和弧-陆碰撞过程(施彬和刘力，2014)，岩浆岩广泛发育，尤其以加里东期和海西-印支期为主要岩浆作用时期，形成了显著的构造-岩浆岩带(丰成友等，2010)。关于加里东期构造-岩浆作用研究，前人已开展了大量研究工作，并取得了丰硕的研究成果。莫宣学等(2007)通过与北祁连造山带进行对比，划分了东昆仑造山带加里东期造山旋回的四个阶段，分别为洋盆打开与扩张(579～518Ma)、俯冲造山(508～450Ma)、碰撞和后碰撞(413～380Ma)、D3后造山崩塌隆升，建立起了东昆仑加里东期运动的构造框架。在早寒武世时原特提斯洋已处于打开和扩张阶段(Yangetal.，1996；陆松年，2002)，东昆仑地区加里东期原特提斯洋壳的俯冲消减多限定在早寒武世末-早志留世初期(李怀坤等，2006；陈能松等，2008；任军虎等，2010；张亚峰等，2010；崔美慧等2011；刘战庆等，2011；刘彬等，2013a, b；Liuetal.，2017)。对于东昆仑地区原特提斯洋的洋盆打开、扩张、消减，前人研究资料相对较多，但对其最终闭合、碰撞及后碰撞等事件的年代学和岩石地球化学方面的证据相对较少，前人通过对该区A型花岗岩研究，认为其多发生于晚志留世-中泥盆世(赵振明等，2008；刘彬等，2012；王冠等，2014；王艺龙等，2018)，而加里东期双峰式侵入岩并没有引起关注。笔者在东昆仑五龙沟地区工作期间发现了晚志留世双峰式侵入岩体，并对其开展相应研究。通过对该区双峰式侵入岩体的年代学、地球化学、构造背景等研究，可为东昆仑加里东晚期-海西早期的构造岩浆演化提供重要资料。同时，东昆仑地区新近发现的夏日哈木、石头坑德、冰沟南等铜镍硫化物矿床的镁铁质-超镁铁质成矿岩体均形成于晚志留世-早泥盆世(王冠等，2014；张照伟等，2018；段雪鹏等，2019；潘彤和张勇，2020)，本研究亦可为该区岩浆铜镍硫化物矿床成矿作用和构造背景研究提供理论依据。

1 地质背景及岩石学特征

东昆仑五龙沟地区位于青海省柴达木盆地南缘，东昆仑中段北缘，大地构造位置处于昆中断裂和昆北断裂夹持的构造岩浆岩带上(图1a, b)，该区地质构造演化历史复杂，经历了极其复杂的洋-陆俯冲碰撞、逆冲推覆造山演化过程，区内构造线总体呈NWW向展布。区内地层从元古代-中生代均有发育，其中以古元古代金水口群角闪岩相-麻粒岩相深变质岩系(黑云斜长片麻岩、黑云石英片岩、石英片岩等)最为发育。区域上岩浆活动强烈，分布广泛，具有期次多、类型全等特点，岩体分布受区域构造演化控制明显(陆露，2011)，主要展布方向呈NW向展布，与区内构造线方向一致，加里东期、华力西期、印支期-燕山期为区内的主要岩浆活动期(郝娜娜等，2014)。

图1 东昆仑大地构造位置图(a)和东昆仑地区岩浆构造图(b，据Dong et al.，2018; Xin et al.，2018修编)

东昆仑五龙沟地区碱长花岗岩岩体出露面积大于100km2，形态近似圆状，长轴方向近NW向，长约13km，宽约8km(图2)。岩体中心与边部的矿物组成、矿物粒度变化不大，岩性比较均一，野外观察岩石新鲜面为肉红色，岩体侵位到金水口群变质岩中，与围岩界限明显。辉长岩和橄榄辉长岩产于同一岩体，整个基性岩体出露面积约5km2，侵位到金水口群白沙河岩组及万保沟组大理岩凝灰岩中，长轴方向近EW向，长约3.8～4.6km，宽约0.6～1.2km(图2)。辉长岩、橄榄辉长岩与该区超基性岩组成辉长岩-橄榄辉长岩-辉石岩-橄辉岩-橄榄岩杂岩体，岩体整体往南-南东倾覆。平面上，岩株中部为中-粗粒、堆晶结构的橄榄岩相，在橄榄岩相外围为中-细粒的橄辉岩相，再外围为中-细粒辉石岩相，最外围为中-细粒辉长岩相，其岩相整体变化特征由中部向外，岩性由超基性向基性过渡，粒度由中-粗粒到中-细粒变化。杂岩体由北向南出露主要岩性依次为浅色辉长岩-辉石岩-橄辉岩-橄榄岩，岩体南端基性程度较高，辉石岩、橄辉岩、橄榄岩等蛇纹石化明显，并具有堆晶结构，局部可见少量闪长岩、辉长岩及片麻岩捕虏体(四川省地质矿产勘查开发局一〇八地质队，2015(1)四川省地质矿产勘查开发局一〇八地质队.2015.青海省都兰县石头坑德铜镍矿普查报告)。碱长花岗岩岩体与辉长岩、橄榄辉长岩杂岩体直线距离最近约8km，野外未发现碱长花岗岩和辉长岩、橄榄辉长岩露头的接触界线。

本次研究选取了东昆仑五龙沟地区的碱长花岗岩(W1-W5取样位置：36°05′N、95°54′E；W6-W7取样位置：36°03′N、95°55′E)和石头坑德地区辉长岩(取样位置N35°56′、E96°09′)、橄榄辉长岩(取样位置35°56′N、96°11′E)(图2)为研究对象进行采样，样品均采自新鲜的岩石露头(图3a-c)。

图2 五龙沟-石头坑德地区地质简图，示取样位置(据青海省地质调查局，2014(2)青海省地质调查局.2014.1:250000地质图修编)

图3 东昆仑五龙沟地区双峰式侵入岩野外(a-c)及显微照片(d-f)

辉长岩呈浅灰色，中细粒结构，块状构造。主要由基性斜长石(d=0.35～2.65mm，含量38%～42%)、单斜辉石(d=0.22～1.85mm，含量15%～25%)、斜方辉石(d=0.25～1.28mm，含量20%～32%)组成，此外含少量黑云母(d=0.15～0.86mm，含量约5%)、石英等(0.06～0.55mm，含量<5%)，副矿物偶见磷灰石和不透明铁质矿物。镜下观察岩石具辉长结构，斜长石呈自形-半自形板条状，聚片双晶和卡钠复合双晶发育，主要为拉长石；单斜辉石斜方辉石呈半自形-他形粒状，单斜辉石可见文象结构，斜方辉石多被包裹或镶嵌在自形的斜长石中。

橄榄辉长岩呈暗灰色，自形中细粒结构，块状构造。主要由基性斜长石(d=0.55～2.15mm，含量20%～30%)、单斜辉石(d=0.37～1.85 mm，含量20%～30%)、斜方辉石(d=0.45～2.05mm，含量20%～35%)、橄榄石(d=0.25～1.28mm，含量15%～20%)组成，副矿物偶见磷灰石和不透明铁质矿物。斜长石呈自形-半自形板条状，聚片双晶和卡钠复合双晶发育，主要为拉长石；单斜辉石斜方辉石呈半自形-他形粒状，单斜辉石部分发生蛇纹石化，斜方辉石多被包裹或镶嵌在自形的斜长石中，部分发生蛇纹石化；橄榄石呈他形粒状，边界见蛇纹石化的辉石组成的反应边结构。

碱长花岗岩呈肉红色，中细粒结构，块状构造。主要由石英(d=0.35～2.85mm，含量40%～45%)、碱性长石(d=0.58～3.86mm，含量45%～52%)组成，此外含少量斜长石(d=0.55～1.52mm，含约5%)、黑云母(d=0.15～0.75mm，含量<3%)，副矿物主要为微量磷灰石、锆石、榍石及不透明金属矿物。石英呈他形粒状，可见波状消光现象；碱性长石他形板状，条纹结构发育，主客晶分别为碱性长石和钠长石。

2 分析方法

本次研究对东昆仑五龙沟地区新鲜的、有代表性的辉长岩、橄榄辉长岩、碱长花岗岩分别进行了全岩主微量元素测试、锆石U-Pb定年、全岩Sr-Nd 同位素测试；对碱长花岗岩和辉长岩进行了锆石Hf同位素测试。测试分析在北京燕都中实测试技术有限公司完成。

锆石U-Pb定年选择有代表性的样品按照标准程序进行锆石分选，将选出的锆石放在双目镜下挑纯，选出晶形和透明度较好的锆石置于环氧树脂中制靶，并进行透反射和阴极发光照相。测试利用LA-ICP-MS分析完成，激光剥蚀系统为New Wave UP213，ICP-MS为布鲁克M90，束斑直径根据实际情况选择25μm。同位素定年采用锆石标准91500和Plesovice作为外标进行同位素分馏校正；锆石微量元素含量利用SRM610作为多外标、Si作内标的方法进行定量计算。

全岩主微量元素测试选取无蚀变及脉体穿插的新鲜样品粉碎至200目以备测试。主量元素测试使用XRF(Zetium，PANalytical)进行，测试数据误差小于1%；微量元素测试使用ICP-MS(M90，analytikjena)进行，测试数据误差小于5%，部分挥发性元素及极低含量元素的分析误差小于10%。

全岩Sr-Nd 同位素测试使用Thermo Fisher Scientific公司多接收电感耦合等离子体质谱仪Neptune Plus MC-ICP-MS分别测定87Sr/86Sr值和143Nd/144Nd值，用88Sr/86Sr值(8.373209)和143Nd/144Nd值(0.7218)对测定的87Sr/86Sr值和143Nd/144Nd值进行校正，87Sr/86Sr值和143Nd/144Nd值的不确定度为2σ。

锆石原位Lu-Hf同位素测试采用美国热电Neptune-plus MC-ICP-MS与NewWave UP213激光烧蚀进样系统完成测试。锆石剥蚀使用频率为8Hz，能量为16 J/cm2的激光剥蚀31s，剥蚀出直径约30μm的剥蚀坑。每个测试点的173Yb/172Yb平均值用于计算Yb的分馏系数，然后再扣除176Yb对176Hf的同质异位素干扰。173Yb/172Yb的同位素比值为1.35274。

3 测试结果

3.1 锆石U-Pb定年

本次研究分别对辉长岩、橄榄辉长岩、碱长花岗岩进行了锆石U-Pb定年以厘定其成岩时代，测试结果见表1。

表1 青海东昆仑五龙沟地区双峰式侵入岩锆石LA-ICP-MS定年分析结果

锆石阴极发光(CL)图像显示，辉长岩的锆石大体可以分为两类，一类为具有较均匀的内部结构、可见条带状吸收的锆石(图4a；点1至点18)，形态以长条状为主，显示了基性岩浆锆石的特征；粒度多为100～120μm，Th/U值范围介于0.08～0.54之间，绝大多数大于0.1，显示其为岩浆成因特征。18粒锆石的206Pb/238U年龄为406±8～413±6Ma，所有数据均落在谐和线上(图5a)，对应的加权平均年龄为409.3±2.7Ma(MSWD=0.1，图5b)，属于早泥盆世，代表了辉长岩的结晶年龄。另一类为具有明显的环带、形态不规则、内部结构多样的锆石(图4a；点19至点27)，206Pb/238U年龄为428±6～435±6 Ma(图5a)，应为辉长岩侵位过程中捕获的锆石。结合区内加里东期岩浆侵入活动频繁，已存在形成于425Ma辉长岩(张照伟等，2018)以及本次研究岩体中存在闪长岩、辉长岩捕掳体的事实，认为该捕获锆石可能来自形成晚奥陶世至中-晚志留世的辉长岩。

图4 东昆仑五龙沟地区双峰式侵入岩锆石阴极发光照片

图5 东昆仑五龙沟地区双峰式侵入岩锆石U-Pb年龄谐和图与加权平均年龄图

与辉长岩类似，橄榄辉长岩的锆石CL图像也显示其存在两种类型的锆石。一类为浑圆状、具有较完整的晶形、较均匀的内部结构的锆石(图4b；点1至点13)，表现出了基性岩岩浆锆石的特征；该类锆石的粒度为80～100μm，Th/U值介于0.14～1.22，为岩浆成因锆石。13粒锆石的206Pb/238U年龄为408±5～411±5Ma，所有数据均落在谐和线上(图5c)，对应的加权平均年龄为410.3±2.7Ma(MSWD=0.1，图5d)，属于早泥盆世，代表了橄榄辉长岩的结晶年龄。另一类锆石形态不规则、部分具有明显环带(图4b；点14至点20)，206Pb/238U年龄为429±4Ma～435±6Ma(中-晚志留世；与辉长岩的捕获锆石年龄一致)以及445±5Ma(晚奥陶世)，也应为橄榄辉长岩侵位过程中捕获的锆石(图5c)。

锆石阴极发光(CL)图像显示，碱长花岗岩的锆石粒度为120～180μm，锆石内部可见清晰的比较窄的振荡环带，部分为断头晶(图4c)，表现出碱性花岗岩锆石的特征；Th/U值范围介于0.47～0.96，为岩浆成因的锆石。23粒锆石的206Pb/238U年龄为406±5Ma～415±5Ma，所有数据均落在谐和线上，对应的加权平均年龄为410.5±1.8Ma(MSWD=0.3，图5e)，属于早泥盆世，代表了碱长花岗岩的侵位年龄。

综上可知，东昆仑五龙沟地区在早泥盆世(约409～411Ma)发生了大规模的岩浆侵入事件。

3.2 岩石地球化学

3.2.1 主量元素特征

碱长花岗岩、辉长岩及橄榄辉长岩的元素地球化学分析结果见表2。其中，碱长花岗岩具有高SiO2(74.74%～77.14%)、高K2O(4.32%～4.97%)、贫MgO(0.08%～0.11%)、CaO(0.18%～0.40%)的特征，Na2O=3.57%～4.19%，Na2O+K2O=8.32%～8.76%，FeOT/MgO=13.82～27.13，在SiO2-K2O图解位于高钾钙碱性区域(图6a)。辉长岩SiO2=52.19%～52.55%，Al2O3=8.8%～9.21%，MgO=19.60%～19.81%，Na2O+K2O=1.49%～1.60%，属于钙碱性系列，TAS图解中位于闪长岩和辉长岩交界处，镜下观察其主要由基性斜长石、单斜辉石、斜方辉石组成，并具有辉长结构(图3e)，显示出辉长岩的特点，因此将其定名为辉长岩。橄榄辉长岩SiO2=42.26%～47.43%，Al2O3=11.65%～12.18%，MgO=25.16%～27.91%，Na2O+K2O=1.11%～1.68%，SiO2-K2O图解位于低钾拉斑质系列区域，并且辉长岩和橄榄辉长岩均落入亚碱性区域(图6b)。在AFM图解中(图7)，基性岩样品均落入镁铁质-超镁铁质堆晶岩区域，反映了该杂岩体形成过程中经历了堆晶作用。TAS图解中(图6b)，该组岩石样品分别落在花岗岩、辉长岩(闪长岩)以及橄榄辉长岩区内，SiO2分布区间主要集中在74.74%～77.14%和42.26%～52.55%两个区间内，在55%～70%之间出现明显的间断，且该组岩石为同一时期同一地区产出的，具有明显的双峰式侵入岩特征。

表2 东昆仑五龙沟地区双峰式侵入岩主量(wt%)和微量(×10-6)元素分析结果

图6 东昆仑五龙沟地区双峰式侵入岩SiO2-K2O图解(a，据Peccerillo and Taylor，1976)和TAS图解(b，据Irvine and Baragar，1971)

图7 东昆仑五龙沟地区基性岩AFM图解(据Irvine and Baragar，1971)

3.2.2 稀土、微量元素特征

在球粒陨石标准化稀土元素配分曲线中(图8a)，碱长花岗岩、辉长岩和橄榄辉长岩均显示为轻稀土元素相对富集、重稀土元素相对亏损的配分模式。碱长花岗岩具有A型花岗岩典型的略右倾的“海鸥”型稀土配分模式，稀土元素总量相对较高，∑REE=146.0×10-6～543.4×10-6，平均306.8×10-6，轻重稀土元素分异较明显(La/Yb)N=1.14～7.82，具有明显的Eu负异常(δEu=0.002～0.005)，Ba、Sr，Eu亏损可能与源区中斜长石的残留有关；辉长岩稀土元素总量∑REE=15.75×106～25.02×106，平均60.83×106，轻重稀土元素分异较明显(La/Yb)N=2.79～8.92，无明显的Eu异常；橄榄辉长岩稀土元素总量∑REE=15.68×106～33.99×106，平均23.37×106，(La/Yb)N=3.37～4.75，无明显的Eu异常。在原始地幔标准化微量元素蛛网图中(图8b)，碱长花岗岩相对富集Rb、Th、U、K、Pb、Nd等元素，相对亏损Ba、Nb、Ta、Sr、P、Eu、Ti等元素；辉长岩和橄榄辉长岩相对富集Rb、U、K、Pb、Nd等，相对亏损Ba、Nb、Ta、P等(图8b)，呈现出Th>Ta，La>Ta的特征。

图8 东昆仑五龙沟地区双峰式侵入岩球粒陨石标准化稀土元素配分图解(a，标准化值据Boynton，1984)及原始地幔标准化微量元素配分图解(b，标准化值据Sun and McDonough，1989)

3.3 全岩Sr-Nd同位素

本次研究对上述3类岩石进行了全岩Sr-Nd同位素分析，测试结果列于表3。结果表明，碱长花岗岩样品εNd(t)值介于-3.33～-3.65之间，对应的两阶段Nd模式年龄为tDM2=1277～1303Ma；辉长岩和橄榄辉长岩87Sr/86Sr初始值介于0.708928～0.710040之间，εNd(t)值介于-3.38～-6.21之间，在ISr-εNd(t)判别图解中落入EMII富集地幔区域(图9)，与同时期的跃进山辉长岩、夏日哈木辉长岩、夏日哈木I号岩体橄榄岩-辉石岩-二辉岩-辉长岩等的演化趋势一致(刘彬等，2012；姜常义等，2015；张照伟等，2018)。同时，辉长岩和橄榄辉长岩的两阶段Nd模式年龄tDM2=1282～1511Ma，均大于其结晶年龄，表明岩浆演化过程中受到过古老地壳物质的混染(吴福元等，2007)。

表3 东昆仑五龙沟地区双峰式侵入岩Sr-Nd同位素分析结果

3.4 Hf同位素

本次研究对辉长岩和碱长花岗岩进行了锆石Lu-Hf同位素分析，结果见表4。结果表明，辉长岩176Yb/177Hf比值为0.007276～0.047558，176Lu/177Hf比值为0.000211～0.001332，176Hf/177Hf比值为0.282406～0.282583，εHf(t)=-4.04～+2.83(大部分为负值)，对应的一阶段模式年龄为tDM(Hf)=927～1352Ma，两阶段模式年龄tDM2(Hf)=1236～1914Ma，均远大于其结晶年龄，表明其来自富集地幔或受过地壳物质的混染；碱长花岗岩176Yb/177Hf比值为0.047410～0.123908，176Lu/177Hf比值为0.001159～0.002945，176Hf/177Hf比值为0.282530～0.282608，εHf(t)=-0.03～2.45，对应的两阶段Hf模式年龄tDM2(Hf)=1245～1396Ma。

表4 东昆仑五龙沟地区双峰式侵入岩锆石Hf同位素分析结果

4 讨论

4.1 岩石成因

4.1.1 基性岩成因

辉长岩与橄榄辉长岩样品的Mg#=82～84，相容元素Cr 和Ni 的含量较高，说明岩浆经历了橄榄石、辉石堆晶作用，这与野外岩相学观察一致。辉长岩和橄榄辉长岩具有较高的ISr值(0.708928～0.710040)和Rb/Sr值(0.07～0.41)，相对富集Rb、Nd、U、K等元素，同时锆石εHf(t)值主要集中在-4.04～+2.83(大部分为负值)，均显示出富集地幔特征，在ISr-εNd(t)判别图解(图9)中样品点落于EMII型富集地幔区域。辉长岩与橄榄辉长岩具有低的Ce/Y比值(0.98～1.58，平均值1.20)，指示其可能来自尖晶石-石榴子石相稳定区域(Mckenzie and Bickle，1988)，岩石低的(Tb/Yb)N值(1.03～1.43，平均1.21)，则可以进一步证实岩石来自尖晶石二辉橄榄岩部分熔融(Wangetal.，2002)，在图10a中，辉长岩与橄榄辉长岩数据大部分位于尖晶石+石榴子石二辉橄榄岩相区域(更靠近尖晶石二辉橄榄岩相)，部分熔融程度大概在3%～20%，说明基性岩浆来自尖晶石二辉橄榄岩相，岩浆源区深度大于80km(Ellam，1992)。

图9 东昆仑五龙沟地区基性岩ISr-εNd(t)判别图解(据Zindler and Hart，1986)

图10 东昆仑五龙沟地区基性岩Sm/Yb-La/Sm(a，据Dilek and Furnes，2011)、(La/Nb)PM-(Th/Ta)PM(b，据Neal et al.，2002)、Th/Yb-Ta/Yb(c，据Pearce，2008)和Nb/Zr-Th/Zr(d，据Woodhead et al.，2001)图解

幔源岩浆在侵入到地壳的过程中，通常会与围岩发生不同程度的同化混染作用，在野外观察中发现基性岩杂岩体中橄榄岩、橄辉岩、辉长岩等基性岩体中局部有少量闪长岩、灰白色辉长岩及片麻岩捕虏体，锆石U-Pb测年结果显示，辉长岩与橄榄辉长岩中也均捕获了晚奥陶世至中-晚志留世(428～445Ma)形成的锆石，这都是该区存在同化混染的直接证据，此外，岩石微量元素呈现出Th>Ta，La>Ta和Nb、Ta 亏损的特征，也暗示岩浆演化过程中存在一定程度的陆壳混染的影响(钱青和王焰，1999)。(La/Nb)PM-(Th/Ta)PM图解(图10b)显示参与混染的地壳物质主要来自于上地壳，这也与辉长岩、橄榄辉长岩Nb、Ta 亏损并捕获片麻岩捕虏体特征相一致。通常认为高La/Sm 值(>4.5)表明地壳物质的混染较强，La/Sm<2则极少受到地壳物质混染(Lassiter and DePaolo，1997)，五龙沟地区基性岩La/Sm 为3.29～4.53(平均值3.77)，指示遭受混染程度并不是特别强烈。需要指出的是该区基性岩石强烈亏损Nb、Ta，不排除这种亏损特征与同化混染作用有关，但由于同化混染作用存在不均匀性，与Nb、Ta亏损在辉长岩和橄榄辉长岩中的普遍性并不协调，所以这种亏损也反映了岩浆固有的特征。

研究表明，岩浆演化和侵位过程中，同化混染和结晶分异往往同时出现，即AFC过程(DePaolo，1981；Halamaetal.，2004；Miretal.，2011)，在主量元素哈克图解上(图11)，MgO 与SiO2、Al2O3、CaO、Na2O + K2O呈一定的负相关性，与FeOT呈正相关性，表明岩浆运移过程中经历了橄榄石、辉石的分离结晶作用，MgO与Ni呈现良好的正相关性，反映了岩浆发生了橄榄石的堆晶作用；未显示明显的Eu负异常，说明斜长石的分离结晶作用不明显。

图11 东昆仑五龙沟地区基性岩主量元素哈克图解

辉长岩和橄榄辉长岩富集轻稀土元素，轻、重稀土元素之间分馏明显，富集Rb、Th、Pb、K等大离子亲石元素，相对亏损Nb、Ta、Zr等高场强元素，这些特征显示了岛弧环境岩浆岩的地球化学特征(Henderson，1984)。此外，Th、Nb、Zr、Yb等元素具有相近的分配系数，分离结晶作用不会影响岩浆中相关元素的分配系数，因此，这些元素对的比值也可以指示源区的特征：辉长岩与橄榄辉长岩样品具有较高的Th/Yb比(0.95～3.41)，在Nb/Yb-Th/Yb图解中(图10c)样品点均明显偏离MORB-OIB演化线，显示了与岛弧环境有关的信息，同时，在Nb/Zr-Th/Zr图解中(图10d)样品大部分落入俯冲交代流体区域，说明岩浆源区受到俯冲组分的影响(Pearce and Peate，1995)。结合其Sr-Nd同位素属富集型地幔范畴特点，判断岩浆源区应该是位于消减板片之上的地幔楔，并且在岩浆生成之前该地幔楔受到了俯冲带流体作用的影响，流体来源可能为加里东期原特提斯洋俯冲洋壳的脱出流体。

综上，笔者认为基性岩石起源于富集地幔部分熔融，并且源区继承了早期大洋俯冲阶段流体交代积累影响，在上升侵位的过程中受到了地壳物质的同化混染并发生了分离结晶最终固结成岩。

4.1.2 碱长花岗岩成因

五龙沟地区碱长花岗岩具有高Si、K，贫Al、Mg、Ca，贫Sr、Eu、Ba、Ti、P等特点，FeOT/MgO=13.82～27.13，显示出典型的A型花岗岩特征(张旗等，2012)。A型花岗岩的成因前人总结主要有以下3种观点(Frost and Frost，2011；刘彬等，2013a)：①长英质地壳的部分熔融；②玄武质岩浆分异；③玄武质岩浆同化地壳物质并发生分异。通常，由陆壳沉积物熔融形成的花岗岩K2O/Na2O>1，而由幔源岩浆分异或地壳中以火成岩为源区重熔形成的花岗岩K2O/Na2O<1(路凤香和桑隆康，2002)。五龙沟地区碱长花岗岩具有高SiO2(74.74%～77.14%)，K2O/Na2O>1的特征，在野外观察中也未发现岩体发育基性岩包体，表明其不是幔源岩浆分异或幔源岩浆与地壳熔体混合成因。Kingetal.(1997)对澳大利亚Lachlen褶皱带的A 型花岗岩进行研究发现其Zr含量为301×10-6，并认为一些高分异铝质A型花岗岩Zr含量甚至低于200×10-6，五龙沟地区碱长花岗岩样品Zr含量为121×10-6～437×10-6，平均225×10-6，并且岩石具有高度的负Eu异常(δEu=0.002～0.005)，强烈亏损Ba、Sr、P和Ti等特点，显示岩浆经历了强烈的分异演化。碱长花岗岩铝指数(A/CNK)为0.97～1.04(平均为1.01)，显示准铝-弱过铝质花岗岩特征，基本可以排除富铝的泥质岩石重熔的可能。花岗岩Rb/Sr比值在15.58～55.84之间，显示出明显的壳源特征(壳源岩浆比相对幔源岩浆更高，幔源岩浆Rb/Sr比通常小于0.05，幔壳混合源介于0.05～0.5之间，壳源岩浆大于0.5；Taylor and McLennan，1985)，此外，元素Ba相对于Rb和Th明显亏损，Nd/Th值在0.73～2.54之间，接近壳源岩石3(Beaetal.，2001)，Nb/Ta值10.04～14，平均11.96，明显不同于幔源岩浆值17.5，而更接近陆壳岩石值11(Taylor and McLennan，1985；Hofmann，1988；Green，1995)，Sr-Nd同位素测试结果岩石具有负的εNd(t)值(-3.33～-3.65，图12a)，这些特征均表明五龙沟地区花岗岩岩浆来源为壳源岩浆，这也与区域上同时期冰沟A型花岗岩(刘彬等，2013a)具有相似的Sr-Nd同位素特征和来源。

图12 东昆仑五龙沟地区双峰式侵入岩εNd(t)图解(a)和锆石Hf同位素组成图解(b)

综上所述，笔者认为五龙沟碱长花岗岩岩浆起源于低压、高温条件下的中元古代长英质地壳部分熔融，热源为基性岩浆底侵供热。

按照岩浆源区，通常认为双峰式火成岩成因有两种(范超鹏，2016):一种为酸性岩和基性岩来自异源母岩浆，二者在空间上的叠加可能仅与一次热事件有关(Sigurdsson，1977)，两个岩石单元一般具有不同的岩浆来源和岩石地球化学特征，微量元素、Sr、Nd和Hf同位素特征差异很大(Davies and MacDonald，l987)，宏观上酸性岩体积要大于基性岩体积；另一种为酸性岩和基性岩来自同源母岩浆，基性岩浆经过后期的分离结晶作用形成SiO2间断的两种岩浆，但二者具有相似的微量元素特征，母岩浆在后期的岩浆演化过程中，基本上没有或者很少受到地壳的混染作用，宏观上基性单元的岩石量要远大于酸性岩。东昆仑五龙沟地区碱长花岗岩显示出壳源特征，基性岩端元则来自于富集地幔并受到地壳的混染，表明该套双峰式侵入岩基性端元与酸性端元来自不同源区，这也与野外观察到酸性岩体出露大于基性岩体地质特征相一致。

4.2 构造环境探讨

东昆仑地区为中央造山带的重要组成部分(许志琴等，2006；杨经绥等，2010)，研究显示，新元古代-早古生代原特提斯洋俯冲消减和石炭纪晚期-三叠纪古特提斯洋俯冲消减是该区最重要的两次造山运动，这两次运动奠定了东昆仑地区现今构造格局(王艺龙等，2018)。因此研究该区古生代构造史对于研究原特提斯洋在我国西部的构造演化都具有重要意义。

研究表明，寒武纪时期，东昆仑地区沿中元古代闭合洋盆再次发生裂解，原特提斯洋开始形成，清水泉蛇绿岩残片中辉长岩的锆石TIMS U-Pb年龄为518～522Ma(Yangetal., 1996；陆松年，2002)，证明了早寒武纪时期昆中洋盆的存在，研究区内主要发育沙松乌拉组变质砂岩、大理岩等(施彬和刘力，2014)。中寒武世时原特提斯洋开始进入俯冲阶段(莫宣学等，2007；张亚峰等，2010)，随着洋壳持续的俯冲，区域上形成了大量与洋壳俯冲有关的弧岩浆岩体(主要以闪长岩为主，少量辉长岩)，如都兰可可沙地区石英闪长岩(515±4Ma；张亚峰等，2010)和镁铁质-超镁铁质杂岩(509±7Ma；冯建赟等，2010)、旺尕秀辉长杂岩(468±2Ma；朱小辉等，2010)、鸭子泉闪长岩(480±3Ma；崔美慧等，2011)、乌兰乌珠尔黑云母花岗(457.5±2.3Ma；韩志辉等，2021)、阿确墩石英闪长岩(448.8±3.9Ma；李婷等，2018)、清水泉变质火山岩(448±4Ma；Chenetal.，2002)、胡晓钦镁铁质岩(438±2Ma；刘彬等，2013b)等。此外，该区一些变质事件的发生也指明洋壳俯冲过程的存在，如昆中缝合带清水泉麻粒岩(507Ma)为大洋岩石圈俯冲过程中发生高温中压麻粒岩相变质作用形成(李怀坤等，2006)；诺木洪南部金水口群夕线黑云二长片麻岩则为昆北带在460Ma时发生岛弧型低压角闪岩相-麻粒岩相变质作用形成(张建新等，2003)。昆中缝合带内具有岛弧玄武岩特征的胡晓钦镁铁质岩石(438±2Ma，刘彬等，2013b)和清水泉基性岩脉(436.4±1.2Ma，任军虎等，2009)基本代表该区洋壳俯冲最晚的岩浆记录。436～431Ma 陆壳型榴辉岩的出现，则代表昆中洋盆已关闭(孟繁聪等，2015；潘彤和张勇，2020)，张照伟等(2018)获得夏日哈木Ⅲ号榴辉岩436Ma递进变质年龄和409Ma 退变质年龄，前者反映碰撞挤压环境，后者反映伸展环境，这也从侧面限定了昆中洋盆的关闭时间。

祁漫塔格地区出现的与陆陆碰撞有关的环斑花岗岩和二长花岗岩年龄分别为428.5±2.2Ma和430.8±1.7Ma(曹世泰等，2011)，而430Ma±该区开始出现A型花岗岩，A型花岗岩出现时间与昆中428Ma的高压榴辉岩相(Mengetal.，2013)和427Ma的中压角闪岩变质岩相(陈能松等，2008)基本一致(刘彬等，2013a)，这说明428Ma开始该区碰撞强烈的陆陆时期结束，东昆仑地区开始进入后碰撞伸展，并延续至晚泥盆世。在这一时期，该区岩浆岩主要以A型花岗岩、过铝质-强过铝质花岗岩、富镁闪长岩和辉长岩等为主，并且区域上双峰式侵入岩广泛发育，如白干湖A型花岗岩(430～429Ma；高永宝等，2014)、乌兰乌珠尔二长花岗岩(422.5±1.6Ma；韩志辉等，2021)、夏日哈木辉长苏长岩(423±1Ma；王冠等，2014)、红旗沟A型正长花岗岩(420±3Ma；王艺龙等，2018)、鸭子沟二长闪长岩(415.5±2.6Ma；王盘喜等，2020)、冰沟正长花岗(391±3Ma；刘彬等，2012)、祁漫塔格喀雅克登塔格区花岗闪长岩(380.52±0.92Ma；郝娜娜等，2014)等，形成时代上与牦牛山组伸展型磨拉石建造的形成时限423～406Ma基本一致(张耀玲等，2010)，上述特点均代表了典型的后碰撞伸展构造背景。

五龙沟地区辉长岩和橄榄辉长岩岩石地球化学显示为钙碱性系列-低钾拉斑质系列，亏损Nb、Ta，上述特点均与碰撞后伸展阶段环境下的基性岩特点相一致(钱青和王焰，1999)。已有的研究证明，在后碰撞伸展环境下所形成的岩浆岩具有岛弧或活动陆缘岩浆岩的地球化学特征，辉长岩和橄榄辉长岩微量元素方面富集Rb、U、Pb、Nd、K等，亏损 Nb、Ta，显示出岛弧岩浆岩的特点，在Hf/3-Th-Ta图解中(图13a)，落入陆缘弧钙碱性玄武岩区域，显示出俯冲流体的交代作用(Pearce and Peate，1995)，并显示出火山弧玄武岩的特点(图13b)，同时在图解区域中明显区别于N型大洋中脊玄武岩，反映出后碰撞伸展阶段刚刚开始。碱长花岗岩在构造判别图解中落入后碰撞区域(图13c, d)，岩石类型显示为A型花岗岩特征；Nb-Y-3×Ga图解进一步确定岩体位于A2型花岗岩区域(图14b)，酸性岩与基性岩构成EMII型富集地幔+A2型花岗岩组合，表明在410Ma左右，该区已经进入碰撞后伸展的阶段。

图13 东昆仑五龙沟地区双峰式侵入岩基性岩Hf/3-Th-Ta(a，据Wood, 1980)、Nb×2-Zr/4-Y(b，据Meschede，1986)、Rb-Y+Nb(c，据Pearce et al.，1984)和R1-R2(d，据Batchelor and Bowden，1985)构造判别图解

图14 东昆仑五龙沟地区碱长花岗岩Zr-104×Ga/Al(a，据Whalen et al.，1987)图解和Nb-Y-3×Ga(b，据Eby, 1992)图解

综合前人研究结果和本次研究，笔者认为五龙沟地区的双峰式侵入岩的形成机制为：加里东期，原特提斯洋在俯冲过程中发生脱水并交代上覆地幔楔，形成了局部具有俯冲组分特征的富集地幔端元；加里东晚期-海西早期该区俯冲作用结束，在后碰撞伸展作用下发生拆沉作用，引发具有俯冲组分特征的富集型岩石圈地幔部分熔融上涌形成基性岩浆，基性岩浆上涌过程中一方面部分基性岩浆发生结晶分异和同化混染作用侵位至地壳形成辉长岩与橄榄辉长岩，另一方面提供热源引起长英质地壳部分熔融，产生A型花岗岩岩浆，形成碱长花岗岩，最终形成了宏观上东昆仑五龙沟地区双峰式侵入岩。