西藏中拉萨地块北部早白垩世晚期控错A型花岗岩的成因及构造环境研究 *

2022-02-14刘洪李光明李文昌黄瀚霄李佑国欧阳渊张向飞周清

岩石学报 2022年1期

刘洪李光明李文昌,3 黄瀚霄李佑国欧阳渊张向飞周清

1. 中国地质调查局成都地质调查中心, 成都 610081 2. 成都理工大学地球科学学院, 成都 610059 3. 昆明理工大学国土资源工程学院, 昆明 650093

青藏高原是全球最为年轻的碰撞造山带，它是开展地球动力学研究的最佳场所。它由一系列近东西走向的地块(微地块)和蛇绿岩混杂带所组成(图1a)(潘桂棠等, 2001, 2004, 2020; Lietal., 2016, 2017b, 2020, 2021;李文昌和江小均, 2020; 莫宣学, 2020; 朱同兴等, 2020; 王立全等, 2021; 张克信等, 2021; 郑有业等, 2021)，由北至南包括：松潘-甘孜地块(SG)、羌塘地块(QT)、拉萨地块(LS)和喜马拉雅地块(HM)(图1a)。这些地块和缝合带受特提斯洋的俯冲和陆-陆碰撞的影响而发育大量的岩浆岩和矿产资源(侯增谦等, 2008; 杨志明等, 2008; Lietal., 2011; 李光明等, 2011; 宋扬等, 2013; 赵希良等, 2013; 蔡青龙等, 2015; 耿全如等, 2015; Huangetal., 2017, 2021a; 刘洪等, 2017, 2018; 唐菊兴, 2019; Caoetal., 2020)，因而备受地质工作者们的关注。拉萨地块(LS)又名冈底斯地块、念青唐古拉地块，北临班公湖-怒江缝合带(BNS)，南接印度河-雅鲁藏布缝合带(YZSZ)(Ding and Lai, 2003; 杨志明等, 2011; 解超明等, 2020; 李光明等, 2021; 耿全如等, 2021)，它又被狮泉河-纳木错蛇绿岩混杂带(SNS)和洛巴堆-米拉山断裂带(LMF)分割为北拉萨地块(NL)、中拉萨地块(CL)和南拉萨地块(SL)三个次级单元(图1b)(杨经绥等, 2007; 潘桂棠等, 2009; 彭建华等, 2013; Huangetal., 2020; 刘洪等, 2020)。班公湖-怒江缝合带(BNS)是拉萨地块(LS)北部重要的构造带，它是特提斯大洋的残留遗迹(Guynnetal., 2013; 许志琴等, 2013; Huangetal., 2018; 刘洪等, 2021a)，对其地球动力学背景的研究为我们探索特提斯洋在中生代演化史提供了一个重要的窗口。班公湖-怒江结合带(BNS)从东向西延伸大于1500km，以复理石、混杂岩和分段蛇绿岩碎片为特征(Panetal., 2012)。一些学者认为班公湖-怒江洋盆在晚二叠世打开，洋壳俯冲于早-中侏罗世(Yin and Harrison, 2000; Tapponnieretal., 2001)，但由于该地区构造特征的复杂性及宏大的规模，班公湖-怒江洋的闭合时限至今仍存在较大的争议(Xuetal., 2015; 刘洪等, 2016, 2018; Huetal., 2017; Lietal., 2017a;Wuetal., 2018)。

图1 研究区地质简图(a)中国大地构造简图(据刘洪等, 2015修改)；(b)冈底斯地质简图(据刘洪等, 2019a, 2021b修改)；(c)控错地区地质简图(据Liu et al., 2018a修改). 1-第四系；2-上白垩统竟珠山组砂砾岩；3-下白垩统郎山组生物碎屑石灰岩夹粉砂岩；4-下白垩统多尼组砂岩、石灰岩和火山岩；5-中二叠统下拉组石灰岩；6-石炭系永珠组砂岩、页岩、石灰岩；7-中-上泥盆统长蛇湖组砂岩夹石灰岩；8-下泥盆统达尔东组石灰岩；9-晚白垩世中酸性岩；10-早白垩世正长花岗岩；11-早白垩世碱性长石花岗岩；12-断层；13-采样位置. NL-北拉萨地块； SNS-狮泉河-纳木错蛇绿岩混杂带；CL-中拉萨地块；LMF-洛巴堆-米拉山断裂带；SL-南拉萨地块. 数据来源：i区(解龙等, 2015); ii区(Li et al., 2020); iii区(王欣欣等, 2021); iv区(Qu et al., 2012; Chen et al., 2014); v区(史少飞等，2019)Fig.1 Geological maps of the research area(a) tectonic map of China (modified after Liu et al., 2015); (b) tectonic map of Gangdise (modified after Liu et al., 2019a, 2021b); (c) Geological map of the Kong Co (modified after Liu et al., 2018). 1-Quaternary; 2-Upper Cretaceous Jingzhushan Fm.: sandstone and conglomerate; 3-Lower Cretaceous Langshan Fm.: bioclastic limestone with siltstone; 4-Lower Cretaceous Duoni Fm.: Sandstone, limestone and volcanic rock; 5-Middle Permian Xiala Fm.: limestone; 6-Carboniferous Yongzhu Fm.: sandstone, shale and limestone; 7-Middle Upper Devonian Changshehu Fm. sandstone intercalated with limestone; 8-Lower Devonian dardong Fm. Limestone; 9-Early Cretaceous intermediate acid rock; 10-Early Cretaceous syenogranite; 11-Early Cretaceous alkali feldspar granite; 12-fault; 13-sampling positions. NL-Northern Lhasa Subterrane; SNS-Shiquanhe-Namuco Ophiolite Melange Zone; CL-Central Lhasa Subterrane; LMF-Luobadui-Milashan Fault Zone; SL-Southern Lhasa Subterrane. Data resources: i area (Xie et al., 2015); ii area (Li et al., 2020); iii area (Wang et al., 2021); iv area (Qu et al., 2012; Chen et al., 2014); v area (Shi et al., 2019)

近年来，在拉萨地块中北部相继发现了一系列白垩纪岩浆活动(曲晓明等, 2004; 黄瀚霄等, 2012; Wangetal., 2014; Guoetal., 2015; Zhengetal., 2015; 彭智敏等, 2015; 张志等, 2017; 高顺宝等, 2021)，这一岩浆事件的研究为探索班公湖-怒江洋的演化史提供了详细的素材。然而，关于这一岩浆作用的地球动力学背景存在许多争论：一些学者认为它与新特提斯洋向北俯冲有关(Coulonetal., 1986; Matte, 1996; 秦臻等, 2019)；也有一些学者认为它是班公湖-怒江洋向南俯冲期间发生板片断离引发的构造岩浆事件(Lietal., 2020)；还有学者认为它们是羌塘地块和拉萨地块碰撞阶段岩浆作用的产物(Zhuetal., 2013; 王保弟等, 2013)。

岩浆岩的源区、形成过程和岩石组合的研究为探索岩浆岩形成时的构造环境提供了重要而有效的信息，特别是关于特殊岩石类型，如A型花岗岩、埃达克岩和板内玄武岩，这可以为区域地质演化提供重要证据(Huppert and Sparks, 1988; Wilson, 1989; Martin, 1999)。而近些年在中拉萨地块北缘发现的一系列早白垩世A型花岗岩(Quetal., 2012; Chenetal., 2014; 解龙等, 2015; 史少飞等, 2019; Lietal., 2020; 王欣欣等, 2021)，为我们研究班公湖-怒江洋在早白垩世的演化提供了良好的对象。在中拉萨地块北部尼玛县控错地区，我们新发现了早白垩世A型花岗岩(图1c)，并对其进行了地球化学、LA-ICP-MS锆石U-Pb地质年代学、锆石Lu-Hf同位素地球化学和全岩Rb-Sr-Sm-Nd-Pb同位素地球化学分析，精确限定其时代归属，探讨了其形成机制、母岩浆源区和形成的大地构造背景，进而丰富对拉萨地块中北部早白垩世A型岩浆作用动力学背景的认识。

图2 控错岩体野外和镜下岩相学特征照片(a)野外宏观特征；花岗岩细粒相(b)和中粗粒相(c)野外露头特征；细粒相(d)和中粗粒相(e)单偏光(-)和正交偏光(+)显微镜下照片. D1d-下泥盆统达尔东组石灰岩等；χργK1-白垩世花岗岩体；Qz-石英；Kfs-钾长石；Pl-斜长石；Bt-黑云母；Hb-角闪石Fig.2 Field photographs and photomicrographs illustrating the petrographic characteristics of the Kong Co granite(a) the field outcrop of the Kong Co granite; the field outcrops of the fine-grained facies (b) and the medium-coarse grained facies (c) in granite; microphotographs under plane polarized light (-) and perpendicular polarized light (+) of the fine-grained facies (d) and the medium-coarse grained facies (e) in granite. D1d-Lower Devonian Daerdong Fm. limestone, etc; χργ K1-Cretaceous granite; Qz-quartz; Kfs-K-feldspar; Pl-plagioclase; Bt-biotite; Hb-hornblende

1 地质概况及岩体特征

研究区位于中拉萨地块北缘的尼玛县控错地区(图1c)，北临狮泉河-纳木错蛇绿岩混杂带。研究区地质条件较为复杂，整体构造线呈NW-SE向，此外还存在碰撞后伸展作用形成的S-N向构造。区内沿NE-NW方向又可进一步划分为三个构造单元，依次为北拉萨地块、狮泉河-纳木措湖蛇绿岩混杂带以及中拉萨地块(图1c)。区域出露的地层主要为上白垩统竟珠山组(K2j)砂砾岩，下白垩统郎山组(K1l)生物碎屑石灰岩夹粉砂岩、多尼组(K1d)砂岩、石灰岩和火山岩，中二叠统下拉组(P2x)石灰岩，石炭系永珠组(C1-2y)砂岩、页岩、石灰岩，中-上泥盆统长蛇湖组(D2-3c)砂岩夹石灰岩，下泥盆统达尔东组(D1d)石灰岩等。其中达尔东组为一套以深灰色薄层石灰岩、生物碎屑石灰岩、泥石灰岩为主的地层体，含丰富的竹节石Nowakiaacuaria带及腕足类、珊瑚、苔藓虫、三叶虫及海百合等，为深水盆地相沉积(图1c)。研究区白垩纪岩浆活动强烈，出露的侵入岩主要为中酸性岩体，包括早白垩(～105Ma)和晚白垩世(～90Ma)两期(图1c)，这些岩浆岩以复式岩株、岩脉的形式侵位于竟柱山组和下拉组地层中(Liuetal., 2018a; 王欣欣等, 2021)。此外，研究区还有报道了～90Ma和～102Ma的中酸性岩浆活动(Luoetal., 2019; Zengetal., 2020)，地壳熔融随着时间的推移由深变浅，从而为其伸展背景提供了重要证据。

控错早白垩世A型花岗岩体(χργK1)侵入到下泥盆统达尔东组石灰岩中，出露面积约10km2(图1c)。该岩体的岩性主要为碱性长石花岗岩(图1, 图2a-c)，呈灰白色，具有细粒相和中粗粒相两种岩相(图2b, c)，这两种岩相呈渐变过渡的接触关系。主要造岩矿物有石英(35%±)(图2d, e)，呈无色，他形粒状，部分颗粒较为浑圆；钾长石(55%±)，主要为条纹长石(图2d, e)，半自形的宽板状，顺条纹有次生尘化现象，部分条纹长石与石英构成文象结构，部分条纹长石内包裹有斜长石小矿物，可见净边结构；斜长石(5%±)，呈无色，半自形板状，表面轻微次生尘化，镜下观察计算的牌号为5±，长石类型主要为钠长石(酸性斜长石)。同时含少量黑云母、角闪石等透明矿物，以及磁铁矿等暗色矿物。岩体和围岩的接触带发育明显的矽卡岩化，接触带内还普遍可见孔雀石化、黄铜矿化、绢云母化、黄铁矿化和磁铁矿化等热液蚀变现象，附近有曲桑格勒和控错等矽卡岩铜矿化点。

图3 控错花岗岩石学图解(a) QAP图解；(b) K2O-SiO2 图解(Peccerillo and Taylor, 1976; Middlemost, 1985)；(c) A/NK-A/CNK图解(Maniar and Piccoli, 1989)；(d)球粒陨石标准化稀土元素配分图解(标准化值据 Boynton, 1984; (e)原始地幔标准化微量元素蛛网图(标准化值据Sun and MchDonough, 1989). 数据来源：i区据Qu et al., 2012; Chen et al., 2014; 解龙等, 2015; Li et al., 2020; 王欣欣等, 2021. ii区据Zhu et al., 2011; 余红霞等, 2011; Liu et al., 2018a; Ouyang et al., 2017; Dai et al., 2020.图7、图8数据来源同此图Fig.3 Petrological plots of Kong Co granite(a) QAP plot; (b) K2O vs. SiO2 plot (Peccerillo and Taylor, 1976; Middlemost, 1985); (c) A/NK vs. A/CNK plot (Maniar and Piccoli, 1989); (d) chondrite-normalized REE patterns (normalization values after Boynton, 1984); (e) mantle-normalized multi-element diagrams (normalization values after Sun and MchDonough, 1989). Data sources: i area after Qu et al. (2012), Chen et al. (2014), Xie et al. (2015), Li et al. (2020), Wang et al. (2021); ii area after Zhu et al. (2011), Yu et al. (2011), Liu et al. (2018a), Ouyang et al. (2017), Dai et al. (2020). Data resources also in Fig.7 and Fig.8

2 样品及分析方法

8件样品分散采集自控错花岗岩体，采样位置为31°28′01″N、87°31′15″E附近，其中KC01～KC05采自花岗岩细粒相，KC06～KC08采自花岗岩粗粒相。所有样品均进行了野外及手标本观察、薄片鉴定及全岩主量和微量元素分析，对其中2件样品(KC05、KC08)进行了锆石U-Pb年代学分析和Lu-Hf同位素分析，5件样品进行了Rb-Sr、Sm-Nd和Pb同位素测试。

主量和微量元素分析在中国地质大学(武汉)地质过程与矿产资源国家重点实验室(GPMR)分别利用X射线荧光光谱仪(XRF)和Agilent 7700e ICP-MS分析完成，使用的国际标样为AGV-2、BHVO-2、BCR-2和RGM-2。锆石的分选、背散射电子拍摄、锆石阴极发光图像分析、锆石U-Pb同位素定年和微量元素含量在武汉上谱分析科技有限责任公司利用LA-ICP-MS同时分析完成，GeolasPro激光剥蚀系统由COMPexPro 102 ArF 193nm准分子激光器和MicroLas光学系统组成，ICP-MS型号为Agilent 7500a，内标位国际锆石标样91500(1062±4Ma)，详细的仪器参数和分析流程见Zongetal. (2017)。本次分析的激光束斑为32μm。对分析数据的离线处理采用软件ICPMSDataCal (Liuetal., 2008, 2010)完成。锆石原位微区Lu-Hf同位素测试利用Neptune型LA-MC-ICP-MS和Geolas Pro型激光剥蚀系统联用的方法完成的，详细测试流程可参照Mengetal. (2014)。测试束斑直径为32μm，所有测试位置与U-Pb定年点位相同或靠近。国际锆石标样GJ-1(0.282030±40(2σ))作为Lu-Hf同位素测试的监控。Rb-Sr、Sm-Nd、Pb同位素分析在武汉上谱分析科技有限责任公司利完成。同位素分析采用ISOPROBE-T热电离质谱计，单带，M+，可调多法拉第接收器接收。质量分馏用86Sr/88Sr=0.1194校正，标准测量结果：NBS987，0.710250±7，JMC，143Nd/144Nd=0.512109±3，NBS 981,208Pb/206Pb=2.164940±15,207Pb/206Pb=0.914338±7，204Pb/206Pb=0.0591107±2。

表1 控错岩体主量(wt%)和微量(×10-6)元素分析结果

图4 控错花岗岩锆石的典型CL图像(a、d)、协和年龄(b、e)和稀土元素特征(c、f)Fig.4 Cathodoluminescence (CL) images of zircon grains (a, d), concordia diagrams zircon measuring point (b, e) and chondrite-normalized REE patterns (c, f) from Kong Co granite

图5 εHf(t)-U-Pb 年龄图解(底图据Blichert-Toft and Albarède, 1997; Amelin et al., 1999; Griffin et al., 2000, 2002; 吴福元等, 2007)数据来源：i区据Qu et al., 2012; Chen et al., 2014; 解龙等, 2015; Li et al., 2020; 王欣欣等, 2021; ii 区据Zhu et al., 2011; 余红霞等, 2011; 王保弟等, 2012; Meng et al., 2014; Hou et al., 2015; Cao et al., 2021Fig.5 εHf(t) vs. U-Pb ages plot of zircons from Kong Co granite (Base map after Blichert-Toft and Albarède, 1997; Amelin et al., 1999; Griffin et al., 2000, 2002; Wu et al., 2007)Data resources: i area after Qu et al., 2012; Chen et al., 2014; Xie et al., 2015; Li et al., 2020; Wang et al., 2021; ii area after Zhu et al., 2011; Yu et al., 2011; Wang et al., 2012; Meng et al., 2014; Hou et al., 2015; Cao et al., 2021

图6 控错岩体Rb-Sr、Sm-Nd (a, b,底图据Zindler and Hart, 1986; Wilson, 1989; Miller et al., 1999)和Pb (c，底图据朱炳泉，1998)图解Fig.6 Discrimination plots of Rb vs. Sr and Sm vs. Nd (a, b, base map after Zindler and Hart, 1986; Wilson, 1989; Miller et al., 1999) and Pb (c, base map after Zhu, 1998) of Kong Co granite

3 分析结果

3.1 全岩主量元素及微量元素

岩石主量、微量元素测试分析结果显示(表1)，控错花岗岩具有高硅 (SiO2=76.75%～77.51%，平均77.27%)、高钾(K2O=4.61%～4.85%，平均4.77%)、高碱(K2O+Na2O=8.24%～8.57%，平均8.44%)，低钙(CaO=0.28%～0.48%，平均0.35%)、低镁(MgO=0.11%～0.16%，平均0.13%)和中等的铝含量(Al2O3=11.79%～12.22%，平均12.09%), 里特曼指数(σ)为1.96～2.15，平均2.08，属于高钾钙碱性-钾玄岩系列过渡的碱性长石花岗岩(图3a, b)。A/NK值为1.06～1.09，A/CNK值为1.01～1.04，显示微弱的过铝质岩石特征(图3c)。这种从碱性矿物到铝质矿物的共生也出现在中拉萨地块其他早白垩世A型花岗岩体中(图1b, 图3c; Quetal., 2012; Chenetal., 2014; 解龙等, 2015; Lietal., 2020; 王欣欣等, 2021)。镁铁指数(Mg#)为14.92～23.09(平均20.18)，分异系数(DI)为95.50～96.90(平均96.33)，指示岩石分异程度较高。计算得到的控错花岗岩锆石饱和温度为平均875～910℃，平均891℃，表明控错岩体结晶于较高的温度环境。

与中拉萨地块其它A型花岗岩特征相似。控错花岗岩稀土元素总量偏高(ΣREE=85.58×10-6～146.0×10-6，平均113.0×10-6)，轻重稀土元素比值中等(LREE/HREE=2.38～3.15，平均2.96)，(La/Yb)N=1.09～3.15，平均1.86，表明轻稀土富集，而重稀土轻微的亏损，具有明显的负Eu 异常(δEu=0.04～0.09，平均0.06)，轻微的Ce正异常(δCe=0.04～0.09，平均0.06)，在球粒陨石标准化稀土元素分布型式图中(图3d)中，各个样品变化趋势一致，具有向右缓倾的“V型”曲线特征(图3d)。从原始地幔标准化微量元素蛛网图(图3e)中，可以看出，包括控错岩体在内的中拉萨地块A型花岗岩明显相对富集Rb、Th、U、K、Ta、Ce、Nd、Zr、Hf、Sm、Y、Yb和Lu等元素，明显相对亏损Ba、Nb、Sr、P、Eu和Ti。

3.2 锆石U-Pb年龄及Lu-Hf同位素组成

对控错花岗岩CK05(细粒碱性长石花岗岩)和CK08(中粗粒碱性长石花岗岩)两个样品中挑取的锆石进行了U-Pb年代学、微量元素及Lu-Hf同位素测试，测点都选择在韵律环带结构清晰的部位，分析位置和结果见图4、表2、图5、表3和表4。从控错碱性长石花岗岩体样品中提取的锆石具有相似的特征，大多为长度为100～300μm的短柱状晶体。它们的CL图像颜色较深，振荡带密集(图4a, d)，其均匀的纹理、完整的晶体形状和相对较高的Th/U比(0.34～1.23，平均0.63)等典型岩浆锆石特征(Grimesetal., 2007)，以及明显的岩浆锆石稀土配分区线(Hoskin, 2005)，表明所研究的锆石主要为花岗岩结晶过程中形成岩浆锆石，此外，部分锆石样品可能发生了一定程度的热液蚀变作用，导致少量测点显示热液锆石稀土元素的特征(图4f)。

KC05号样品共测试了16个分析点，其206Pb/238U表面年龄在100.3～110.2Ma之间，其加权平均值为104.9±1.4Ma(MSWD=1.5)(图4、表2)，KC08号样品共测试了19个分析点206Pb/238U表面年龄在100.7～109.4Ma之间，其加权平均值为104.6±1.3Ma(MSWD=1.3)(图4、表2)，两个样品的锆石年龄代表了控错花岗岩的结晶年龄，并与Rb-Sr等时线年龄为104.8±4.3Ma(MSWD=0.6)在误差范围内一致(图6)。因此，控错地区的岩浆活动时间应该是早白垩世(约105Ma)。

在锆石U-Pb年代学测试和微量元素测试的基础上，对KC05和KC08样品的部分年代学测点开展了Lu-Hf同位素测试。控错花岗岩锆石176Yb/177Hf值为0.020226～0.056663，平均0.036282(表4)，176Lu/177Hf 变化于0.000804～0.002239，平均0.001416，176Hf/177Hf值为0.282831～0.282891，平均0.282856。根据对应锆石年龄计算的εHf(0)主要分布于2.08～4.20之间，平均为5.15，εHf(t)值为正值(4.26～6.38，平均5.15，图5、表4)，二阶段模式年龄(tDM2)为757～889Ma，平均833Ma(图5、表4)。

图7 控错岩体岩石成因图解(a, b,底图据 Whalen et al., 1987; c-f,底图据Eby, 1992)Fig.7 Discrimination diagrams of genetic types of Kong Co granite (a, b, base maps after Whalen et al., 1987; c-f, base maps after Eby, 1992)

3.3 Rb-Sr、Sm-Nd、Pb同位素

控错花岗岩类的Rb-Sr、Sm-Nd、Pb同位素含量见表5，87Rb/86Sr和87Sr/86Sr分别为91.1801～147.565(平均值为116.547)和0.8768～0.9538(平均值为0.9051)；计算得到的 (87Sr/86Sr)t为0.7194～0.7407(平均值为0.7313)，获得Rb-Sr等时线年龄为104.8±4.3Ma(MSWD=0.6)(图6)，与锆石U-Pb年龄相在一致。147Sm/144Nd和143Nd/144Nd分别为 0.1778～0.1912(平均值为0.1831)和0.5125～0.5125(平均值为0.5125)；计算得到的(143Nd/144Nd)t、εNd(t)分别为 0.512340～0.512352(平均值为0.512346)和-3.39～-3.00(平均值为-3.24)。206Pb /204Pb、207Pb/204Pb和208Pb/204Pb分别为18.948～19.043(平均值为18.991)、15.715～15.728(平均值为15.721)和39.419～39.449(平均值为39.347)；(206Pb /204Pb)t、(207Pb /204Pb)t和(208Pb /204Pb)t分别为18.792～18.845(平均值为18.823)、15.708～15.718(平均值为15.713)和38.870～38.037(平均值为39.979)。

表4 控错岩体LA-ICP-MS 锆石 Lu-Hf 同位素

4 讨论

4.1 岩石类型厘定

图8 控错岩体构造判别图解(a,b,底图据 Pearce et al., 1984; c-e,底图据 et al., 2007; f,底图据Grimes et al., 2007)Fig.8 Tectonic discriminant diagrams of Kong Co granite (a, b, base maps after Pearce et al., 1984; c-e, base maps after al., 2007; f, base map after Grimes et al., 2007)

A型花岗岩在矿物组成上以出现碱性暗色矿物(如钠闪石、铁黑云母、钠质辉石等)为特征，在化学组成上具有低Al2O3(～12.40%)、低CaO(～0.75%)、低MgO(～0.20%)、低A/NK值(～1.08)、高Ca/Al及高含量的高场强元素(HFSE)(邓晋福等, 2009)。控错花岗岩富含钾长石(55%±)，存在黑云母和角闪石等富钾矿物(图2)，在主量元素组成上具有高钾(K2O=4.61%～4.85%，平均4.77%)，高碱(K2O+Na2O=8.24%～8.57%，平均8.44%)，低铝(Al2O3=11.79%～12.22%，平均0.35%)，低钙(CaO=0.28%～0.48%，平均0.35%)，和低镁(MgO=0.11%～0.16%，平均12.09%)的特征(图3、表1)，在铝指数上表现为弱的过铝质特征(A/NK值为1.06～1.09，A/CNK值为1.01～1.04)(图3、表1)，这些特与典型的A性花岗岩矿物和主量元素组成特征(Whalenetal., 1987; 邓晋福等, 2009)相似。在微量和稀土元素方面(图3、表1)，控错花岗岩具有明显的负Eu 异常(δEu=0.04～0.09，平均0.06)，Ce异常不明显，轻稀土略富集于重稀土(LREE/HREE=2.38～3.15, 平均2.96)，在球粒陨石标准化稀土元素配分图中(图3d)中，具有向右缓倾的“V型”曲线特征(图3d)，同时控错花岗岩富含Rb、Th、U、Pb、Zr和Hf，相对亏损Nb、Ta、Ti、P、Ba和Sr，在微量和稀土元素组成上具有A型花岗岩的特征(Whalenetal., 1987; Eby, 1990)，并明显相近于中拉萨地块上报道的早白垩世A型花岗岩(图3)，而有别于中拉萨晚白垩世I型花岗岩(图3)。在FeOT/MgO(5.61～10.22，平均7.26)、10000Ga/Al(2.78～2.56，平均2.84)和Zr(91×10-6～145×10-6，平均112×10-6)等花岗岩成因类型判别指标上(图7a, b)，控错花岗岩均表现出A型花岗岩的特征(Eby, 1992)，并与中拉萨地块上报道的其它早白垩世A型花岗岩非常相似(图7a, b)。S型花岗岩平均温度为～764℃，I型花岗岩平均温度为～781℃，而A型花岗岩的形成温度比S型和I型的更高(Kingetal., 1997; 张旗等, 2008)。全岩锆饱和温度计(Watson and Harrison, 2005)计算出控错岩体的形成温度(Tzr)为 875～910℃(平均890℃)，锆石Ti温度计(Schiller and Finger, 2019)计算出控错岩体锆石的结晶温度为848～919℃(平均890℃)，符合A型花岗岩的形成温度范围。综上所述，我们认为，控错早白垩世碱性长石花岗岩为A型花岗岩。

A型花岗岩通常被定义为具有非造山期的、碱性的和无水的等特点，然而，随着研究的进一步发展，学者们发现A型花岗岩不仅在地幔柱或热点型环境中形成，而且也可以在后造山构造环境中形成(Eby, 1990)。因此，岩石学家们根据A型花岗岩的岩石学和地球化学特征，以及其物质来源和构造背景的差异，将其分为A1型和A2型两个亚类(Eby, 1992)，其中A1型的微量元素含量类似于大洋洲玄武岩(OIB)，被认为产于非造山环境(大陆裂谷或板内环境)，而A2型花岗岩的微量元素含量类似于大陆地壳物质和岛弧玄武岩(IAB)，被认为产于同碰撞/碰撞后环境。A1型花岗岩具有很低的Y/Nb和Rb/Nb比值(Eby, 1992)，而A2型花岗岩则相反，具有较高Y/Nb和Rb/Nb比值。与中拉萨地块已报道早白垩世晚期(～105Ma)A型花岗岩相似，控错花岗岩具有出较高的Rb(373×10-6～506×10-6，平均429×10-6)和Y(50.2×10-6～55.0×10-6，平均51.0×10-6)元素含量和相对较低的Nb(24.6×10-6～33.9×10-6，平均28.6×10-6)元素含量，较高Y/Nb值(2.29～4.97, 平均3.57)和Rb/Nb(11.6～18.2, 平均15.2)，在三角图解上落于A2型花岗岩的区域(表1、图7c-f)，在构造环境判别图上它们落于碰撞后花岗岩的范围(图8a, b)。因此，我们认为控错岩体为弱过铝质的高钾钙碱性-钾玄岩系列的A2型碱性长石花岗岩。

图9 控错岩体Rb-Sr-Ba-Eu结晶分异定量计算图Kfs-钾长石；Pl-斜长石；Bt-黑云母；Amp-角闪石Fig.9 Quantitatively fractional crystllization diagrams with Rb-Sr-Ba-Eu of Kong Co graniteKfs-potassium feldspar; Pl-plagioclase; Bt-biotite; Amp-amphibole

图10 控错岩体形成环境示意图(a)据Zhu et al., 2016; Ma et al., 2014; Wang et al., 2016修改；(b)据Cao et al., 2019修改. QT-羌塘地块； BNS-班公湖-怒江结合带；NL-北拉萨地块； SNS-狮泉河-纳木错蛇绿岩混杂带；CL-中拉萨地块；LMF-洛巴堆-米拉山断裂带；SL-南拉萨地块Fig.10 Sketch maps showing the formation environments of Kong Co granite(a) modified after Zhu et al., 2016; Ma et al., 2014; Wang et al., 2016; (b) modified after Cao et al., 2019). QT-Qiangtang Terrane; BNS-Bangong-Nujiang Suture Zone; NL-Northern Lhasa Subterrane; SNS-Shiquanhe-Namuco Ophiolite Melange Zone; CL-Central Lhasa Subterrane; LMF-Luobadui-Milashan Fault Zone; SL-Southern Lhasa Subterrane

4.2 岩石成因探讨

通常认为A型花岗岩可由多种来源的岩浆形成，主要包括：①残余富F或Cl下地壳麻粒岩的部分熔融(Collinsetal., 1982; Whalenetal., 1987; Kingetal., 1997; Yangetal., 2006)；②浅地壳内英云闪长岩或花岗闪长岩材料的低压脱水熔融(深度<15km)(Skjerlie and Johnston, 1992)；③玄武质岩石的部分熔融(Wuetal., 2002, 2018)；④地幔源镁铁质岩浆分离结晶-同化混染作用(AFC)直接形成(Sparks and Marshall, 1986; Foland and Allen, 1991; Turneretal., 1992)；⑤地幔源熔体和地壳熔体的相互作用形成(Kerr and Fryer, 1993)。实验岩石学表明，来自残余富氟源的熔体中通常MgO含量大于TiO2，且具有强过铝质(Creaseretal., 1991; Dooley and Patio Douce, 1996)。与中拉萨地块中其他早白垩世A型花岗岩一样，控错岩体具有较高的TiO2/MgO值(0.51～0.68，平均0.66)和弱过铝质特征，因此，中拉萨地块中这些早白垩世A型花岗岩并非源自残余富F或Cl下地壳麻粒岩材料的部分熔融(排除了上述①的可能性)。控错花岗岩和中拉萨地块已报道的其他A型花岗岩的Ce/Pb、Ce、Nb/Th和Nb特征上，均表现为弧火山岩的亲缘性(图8c, d)，同时Rb-Sr、Sm-Nd, Pb同位素上显示出存在幔源组分的特征(图6)，因此控错岩花岗岩不可能由纯地壳来源(包括浅地壳和下地壳)形成，排除了上述②和③的可能性。控错花岗岩具有高硅 (SiO2=76.75%～77.51%，平均77.27%)的特征，幔源镁铁质岩浆衍生的高硅岩石涉及巨大规模的幔源镁铁质岩浆岩分离结晶作用(Wilson, 1993)。在控错地区，并没有大规模的晚白垩世镁铁质岩浆岩活动的报道，因此控错花岗岩由AFC过程直接形成的可能性不大。

对于控错A型花岗岩的成因，我们倾向于用地幔源熔体和地壳熔体的相互作用来解释，主要证据如下：①控错花岗岩为弱过铝质的碱性长石花岗岩，这种从碱性矿物到铝质矿物的共生也出现在中拉萨地块其他早白垩世A型花岗岩体(图1b、图3c; Quetal., 2012; 解龙等, 2015; Lietal., 2020; 王欣欣等, 2021)，以及阿根廷北部的A型花岗岩中(Shellnutt and Zhou, 2007)，其形成可能与壳幔相互作用和沉积物熔融有关。②控错岩体及中拉萨地块其他早白垩世A型花岗岩体在Rb/Y-Nb图解上显示出地壳混染的趋势(图8e)，控错锆石的U、Yb、Hf和Y同位素显示出地壳来源锆石的特征(图8f)，而地壳内富含高场强元素副矿物的熔融可产生含有控错花岗质岩浆这样高浓度Zr、Ce、Y和Ga的岩浆(Kerr and Fryer, 1993)，此外在微量元素蛛网图中Nb出于处于的低谷，NbN=34.5～47.6(平均40.1)，也表明该花岗岩具有地壳组分(Barthetal., 2000)；③在U-Pb年龄-εHf(t)关系图解(图5)上，控错样品点落在球粒陨石均一储库(CHUR)和地幔演化线之间，并与中拉萨地块其他早白垩世花岗岩相具有相似的源区，正的锆石εHf(t)值(2.9～9.9，平均6.25)、相对年轻的锆石Hf地壳模式年龄(tDM2=529～994Ma，平均567Ma)、相对低的全岩Rb/Y和Nb/Y值(Rb/Y=7.24～10.54，Nb/Y=0.48～0.71)(表2)，和较高的全岩Ce/Pb值(0.97～1.75)(表2、图8d)指示控错花岗岩来源于新生地壳(图5、表5)，并且有地幔物质参与(图5)；④(87Sr/86Sr)ⅰ、εNd(t)同位素特征(图6b)，显示为地幔和下地壳的混合，Pb同位素特征(图6c, d)上，控错花岗岩类样品主要落在班公湖-怒江蛇绿岩以及地壳与地幔混合的铅的区域内。

一般认为元素的亏损是由某种富集该元素的矿物结晶分异引起(Chappell and White, 1992; Chappell, 1999; Wuetal., 2003)，控错花岗岩类亏损Ba、Nb、Sr、P、Eu和Ti，指示其母岩浆经历了显著的铁镁矿物、富钛矿物、富磷矿物等的分离结晶作用。同时，主要氧化物(FeOT、Al2O3、CaO、MgO、P2O5和TiO2)和Sr等元素与SiO2负相关性(图9)指示长石类矿物、铁镁矿物(如角闪石和黑云母)、含磷矿物(如磷灰石)和含钛(如钛铁矿、钛矿和金红石)矿物的分离结晶，表明在该花岗岩形成过程中发生了明显的晶体分馏。这一点得到了分离结晶定量建模的支持(图9)，它表明控错花岗岩可以解释为从分异程度最高的样品(KC06)代表的假定母岩浆成分开始，控错花岗岩浆按照70%钾长石+1%斜长石+9%黑云母+20%角闪石的比例分离结晶(图10)。总之，控错早白垩世A型花岗岩最有可能由通过幔源镁铁质岩浆和初生下地壳再熔融产生的熔体之间的相互作用形成，并在侵位之前经历了显著的分离结晶。

4.3 构造环境分析

前文已分析，控错花岗岩与中拉萨地块已报道早白垩世晚期(～105Ma)A型花岗岩一样，均属于与陆陆碰撞相关的A2型花岗岩，他们最有可能由通过幔源镁铁质岩浆和初生下地壳再熔融产生的熔体之间的相互作用形成，而地幔物质在碰撞后过程中参与成岩作用的主要方式被认为包括板片撕裂和下地壳拆离等(Houetal., 2004)。如前文所述，拉萨地体北部和中部广泛发育的这一系列早白垩世岩浆作用的地球动力学背景存在争议，一部分学者认为其与雅鲁藏布新特提斯洋洋壳板片向北俯冲有关，也有一部分学者则认为其形成于班公湖-怒江洋的洋壳板片的向南俯冲或板片断离作用，还有一部分学者认为它们是羌塘地块和拉萨地块碰撞后岩浆作用的产物。

一些学者认为，雅鲁藏布新特提斯洋向北俯冲涉及低角度俯冲和俯冲角度的增加以及板片回撤事件触发了西藏中部白垩纪岩浆喷发(Zhangetal., 2004, 2012)。在低角度俯冲模型中，早期俯冲角较小，相应的岩浆作用较少。然而，近年来在拉萨地体南部发现了广泛的侏罗纪-白垩纪岩浆作用(Wenetal., 2008; Zhuetal., 2009; Huangetal., 2019; 刘洪等, 2019b, c)，这与洋壳板片低角度俯冲模型相反。其他一些学者认为，雅鲁藏布新特提斯洋向北俯冲是正常的深俯冲，板块撕裂事件触发了早白垩世岩浆喷发(Daietal., 2015)。考虑到拉萨地体在白垩纪之后经历了显著的地壳缩短，早白垩世晚期，北部岩浆弧与南部雅鲁藏布新特提斯洋俯冲带之间的距离超过600km(Murphyetal., 1997; Kappetal., 2005)，这很难用正常俯冲模型来解释。因此我们认为，早白垩世控错岩体的形成，与雅鲁藏布新特提斯洋的演化不相关。最近对狮泉河-纳木错蛇绿岩混杂带的研究表明，该混杂带是班怒-洋俯冲期间形成的弧后洋盆(Zhuetal., 2011; Xuetal., 2014; Zengetal., 2018)，其小规模和快速演化意味着该洋盆地中的俯冲作用与拉萨地体中北部的大规模岩浆活动不相关(Zhuetal., 2008)。

目前班公湖-怒江洋盆关闭的确切时间和俯冲方向等问题仍存在争论(Volkmeretal., 2007; 陈奇等, 2007; 高顺宝等, 2011; 李光明等, 2011; Lietal., 2014；Wangetal., 2016; Huangetal., 2021b, 2021c; Liuetal., 2021)：一些学者认为，班公湖-怒江洋盆在早侏罗世向北俯冲于羌塘地块之下，而拉萨和羌塘地块最有可能在140～130Ma发生弧-弧“软”碰撞(宋扬等, 2019)；同时，拉萨地块北部尼玛地区在具有125～118Ma海相向非海相转变的记录(Kappetal., 2007)、大约在110Ma 就已经为陆内环境(Zhuetal., 2009, 2011)，这意味着拉萨和羌塘地块在此之前已经发生了“软”碰撞；也有学者研究认为，萨地块中部申扎周围约113Ma 的大规模岩浆活动，是向南俯冲的班公湖-怒江洋板块断裂的产物(Chenetal., 2014)；同时还有学者提供了早白垩世早期(138～134Ma)，班公湖-怒江洋尚未彻底闭合的重要证据(Zengetal., 2021)；此外，大面积分布的竟柱山组(K2j)磨拉石的报道(潘桂棠等, 2006)，也暗示在早白垩世晚期以来，拉萨地块北缘已经进入了陆陆碰撞的阶段。尽管这些研究存在很大的分歧，但可以认为在早白垩世晚期，拉萨地块北缘出于碰撞环境。

班公湖-怒江洋两侧广泛发育的中生代中酸性岩浆作用在白垩纪表现出宽泛的εHf(t)值范围(图5)，从早白垩世早期到早白垩世晚期(～105Ma)，εHf(t)从负值陡然上升到正值，岩浆的地幔贡献组分增加，从早白垩世晚期(～105Ma)到晚白垩世(～65Ma)，εHf(t)则有所下降，甚至到达负值，岩浆的地幔贡献组分有所减少，表明早白垩世晚期存在地幔物质加入到岩浆源区的事件，这可能碰撞后与板片窗口的打开密切相关。当地壳下方存在板片窗时，将产生一系列相应的特殊伸展岩组合，包括双峰火山岩、A型花岗岩、板内玄武岩和埃达克岩。在中拉萨地块北部，已报到了与伸展相关的早白垩世晚期(～105Ma)碰撞岩浆作用事件(Quetal., 2012; 吴浩等，2014; Wuetal., 2015a, b；解龙等, 2015; Lietal., 2020; 王欣欣等, 2021)，包括块玄武岩、埃达克岩、双峰火山岩以及A型花岗岩。这些花岗岩的高Ce/Pb比值、正的εHf(t)值、相对年轻的锆石Hf地壳模式年龄，地幔和下地壳混合(87Sr/86Sr)ⅰ、εNd(t)、Pb同位素特征提供了源于地幔源熔体和新生地壳熔体证据。这与该地区的同期伸展相关岩浆作用一致，并提供了存在板片窗的进一步证据。

因此，中拉萨地块中北部早白垩世晚期的岩浆作用既不可能是雅鲁藏布新特提斯洋壳板片平板俯冲或洋脊俯冲的产物，也不可能是班公湖-怒江洋片南向俯冲消减直接的产物，而更可能是班公湖-怒江洋板片俯冲消减闭合之后羌塘-拉萨地块碰撞过程中板片断离的产物。早白垩世晚期，羌塘-拉萨地块碰撞作用下，俯冲板片的断离形成板片窗，引起地幔物质通过板片窗上涌和上覆岩石圈内的伸展，幔源熔体与新生地壳熔体相互作用而形成了中拉萨地块中北部一系列的A型花岗岩(图10)。