甘肃马庄山地区花岗岩类的元素地球化学、U-Pb年代学及Nd-Hf同位素特征

王琦崧 张静 曹双双 于立栋 孙海微

1. 贵州民族大学建筑工程学院，贵阳 5500252. 中国地质大学地质过程与矿产资源国家重点实验室，北京 1000833. 新疆维吾尔族自治区有色地质勘查局七〇四队，哈密 8390001.

中亚造山带(CAOB)，夹持于西伯利亚和塔里木-华北板块之间(图1a)，是世界最大的显生宙增生型造山带之一，经历了漫长而又复杂的增生造山过程(engöretal., 1993; Windleyetal., 2007; Wangetal., 2012; Goldfarbetal., 2014; Zhangetal., 2017; Chenetal., 2020)。北山造山带地处中亚造山带南缘，西邻天山造山带，东接索伦构造带，是剖析中亚造山带南缘构造演化过程的关键区域之一(图1a)。北山造山带南带由花牛山和石板山岛弧组成，记录了古老微陆块的俯冲-碰撞过程(Liuetal., 2011; Zongetal., 2017; Yuanetal., 2018)。北山北带由马鬃山岛弧、旱山岛弧和雀儿山岛弧组成，其性质和构造演化尚存在争议；一种观点认为，晚古生代火成岩与俯冲-增生过程有关(Xiaoetal., 2010; Aoetal., 2016; Han and Zhao, 2018)；另一种观点认为，晚古生代的侵入岩形成于伸展环境(Zhengetal., 2013; Zhangetal., 2017)。基于此，本文对北山北带旱山弧内马庄山地区广泛分布的花岗岩类进行岩石学、元素地球化学、锆石U-Pb年代学及Nd-Hf同位素等研究，结合前人的研究资料，分析研究区花岗岩类的岩石成因和岩浆来源，进一步探讨北山北带的构造演化过程。

1 区域地质背景

北山造山带是由前寒武系的变质基底、古生代岛弧、蛇绿岩带、微陆块和增生杂岩带等拼合而成(毛启贵, 2010; Tianetal., 2014; 郑荣国等, 2016)，经历了古生代-中生代长期的多阶段、复杂的俯冲拼贴过程，是研究中亚造山带的地球动力学演化过程和古亚洲洋的最终闭合时限的关键区域(Xiaoetal., 2010; Songetal., 2013a, b; Tianetal., 2014; Clevenetal., 2015)。受南北板块的构造挤压，造山带、断裂、盆地呈近东西走向，形成了造山带、盆地相间的构造格局(陈衍景, 2000)。由北向南，北山造山带依次可划分为雀儿山弧、黑鹰山-旱山弧、马鬃山弧、双鹰山-花牛山弧、石板山弧和等5个单元(图1b, Xiaoetal., 2010)。雀儿山弧为奥陶纪-二叠纪的活动大陆边缘弧，发育奥陶纪-二叠纪活动陆缘火山岩及碎屑沉积岩地层，并出露古生代侵入岩(聂凤军等, 2002; 刘明强等, 2006; 丁嘉鑫等, 2015)。黑鹰山弧为前寒武纪岛弧，出露石炭系基性-中性火山岩及碎屑岩、古生代及中生代花岗岩(丁嘉鑫等, 2015)，带内发育浅成低温热液金矿床(图1b，如霍勒扎德盖金矿；疏孙平等, 2017)。旱山弧为一古生代的复合型岛弧，主要由高级片麻岩和花岗岩组成(Xiaoetal., 2010)，并发育有多个以浅成低温热液型金矿床为主体的金矿床(如：马庄山、南金山和修翁哈拉金矿床；Chenetal., 2012; Deng and Wang, 2016; Wangetal., 2019; 王琦崧等, 2019)。马鬃山弧为奥陶纪-二叠纪的复合型岛弧，主要由绿片岩-角闪岩相变质岩、火山岩和侵入岩组成(Songetal., 2013c)，该弧内分布金窝子、210等金矿床(Chenetal., 2020)。双鹰山-花牛山岛弧为寒武纪-二叠纪的复合型岛弧，主要发育新元古代-晚古生代地层和侵入岩(王鑫玉, 2017)，已发现国庆、玉山等钨矿床(丁嘉鑫等, 2015)。而最南端的石板山岛弧为向南俯冲的古生代-二叠纪活动大陆边缘弧，主要出露有泥盆系、石炭系和二叠系的碎屑沉积岩、火山岩和火山碎屑岩等(Xiaoetal., 2010)，弧内发育有与韧性剪切带有关的热液金矿床(如：小西弓金矿；江思宏, 2004)。

2 岩体产出及岩相学特征

研究区位于新疆-甘肃交界的泉东山-马庄山一带，大地构造位置属于旱山弧西段，其南北两侧分别受星星峡走滑断层和路井两条大断裂的控制(图1b)。区内出露地层除第四系冲洪积外，主要为蓟县系大理岩、灰岩，长城系星星峡组片麻岩，石炭系灰岩、火山碎屑岩、安山岩、英安岩和碳酸盐岩等(图1c)。岩浆活动以酸性为主，出露二长花岗岩、花岗闪长岩、花岗岩、石英斑岩等，其中以花岗岩类岩石分布最为广泛，多呈岩基、岩株或岩枝状侵入于石炭系火山岩中(图1c、图2a)，局部可见石英脉侵入英云闪长岩(图2g)。野外地质调查和资料显示，各岩相侵入体从早到晚的时间顺序大致为二长花岗岩、花岗闪长岩、花岗岩、石英斑岩，部分侵入体与围岩的外接触带可见烘烤边，表明其接触关系为热侵位。

花岗闪长岩主要出露在马庄山以西约3km处，呈岩基和岩枝产出，灰白色-肉红色(图2b)，中粗粒结构，块状构造。主要由石英(20%～25%)、斜长石(40%～45%)、钾长石(～20%)、少量黑云母(～10%)组成。其中，斜长石表面浑浊，呈自形-半自形柱状，粒径约为0.5～3mm，局部可见环带结构，部分已绢云母化(图2c)；钾长石表面呈肉红色，半自形板状，部分绢云母化，大小为0.3～1.5mm；石英多呈他形粒状，粒径约为0.1～0.3mm，局部可见波状消光；黑云母多呈片状，粒径约0.1～0.5mm。

图1 中亚造山带构造示意图(a, 据Sengor and NataPin, 1996)、北山构造简图(b, 据Xiao et al., 2010修改)和马庄山区域地质图(c, 据新疆维吾尔自治区有色地质勘查局704队, 2018(1)新疆维吾尔自治区有色地质勘查局704队. 2018. 新疆哈密市玉峰金矿预查设计报告修改)

图2 马庄山地区花岗岩类野外及镜下照片(a-c)花岗闪长岩的野外产状、手标本和显微照片(+)，(d-f)钾长花岗岩的野外产状、手标本和显微照片(+)，(g-i)英云闪长岩的野外产状、手标本和显微照片(+).Qtz-石英；Pl-斜长石；Kfs-钾长石；Bt-黑云母；Amp-角闪石；Ser-绢云母Fig.2 Field and petrographic photos of granitoids in the Mazhuangshan area(a-c) occurrence, hand specimen and micrograph of granodiorite (+); (d-f) occurrence, hand specimen and micrograph of K-feldspar granite (+); (g-i) occurrence, hand specimen and micrograph of tonalite (+). Qtz-quartz; Pl-plagioclase; Kfs-K-feldspar; Bt-biotite; Amp-amphibole; Ser-sericite

钾长花岗岩位于泉东山以北1km和玉峰以东约4km处，呈岩珠产出，在空间上总体呈NE向展布，明显受区内断裂控制，出露面积较小，约20km2。岩石呈灰黄色-肉红色、中粗粒结构(图2e)，主要矿物有石英(～40%)、钾长石(～50%)和斜长石(～10%)。其中，钾长石以条纹长石为主，发育格子双晶，粒径约为0.8～3mm；斜长石多呈半自形，斜长石表面浑浊，粒径约为0.5～3mm；石英为灰白色，呈他形粒状，局部可见波状消光，粒径为0.1～1mm(图2f)。

英云闪长岩位于位于玉峰以东约4km处，呈岩枝产出，灰黑色-灰绿色，粒度中等，块状构造(图2h)。主要由黑云母(20%～30%)，石英(～25%)，斜长石(30%～40%)以及角闪石(5%)组成。其中，斜长石多呈自形-半自形，宽板状、聚片双晶发育(图2i)，粒径约为0.5～3.5mm；石英，灰白色，多呈他形粒状，粒径为0.1～0.3mm；角闪石呈他形，含量较少，粒径约为0.2～0.5mm。

二长花岗岩主要出露在研究区东北部的明水附近，呈岩基产出，浅肉红色和灰白色，中-粗粒，似斑状结构，块状构造。主要由微斜长石(20%～30%)、斜长石(25%～35%)、石英(35%～45%)和黑云母(8%～15%)组成。其中微斜长石主要为条纹长石和钾长石，半自形板柱状，卡氏双晶、格子双晶发育；斜长石多呈自形-半自形板状，聚片双晶和环带结构发育，部分已蚀变为绢云母；石英他形粒状；部分黑云母蚀变为绿泥石、绿帘石和黝帘石(Zhangetal., 2017)。

3 样品采集和测试方法

3.1 样品采集

本次研究的花岗质类岩石采自马庄山金矿附近及其以北的10km范围内，岩体围岩主要是石炭系火山碎屑岩、碳酸盐岩和石英斑岩。在野外地质调查的基础上，对钾长花岗岩、英云闪长岩、花岗闪长岩等进行系统采集，样品采自岩体地表或钻孔岩芯，具体位置见图1c。取样过程中，尽量选取新鲜、未蚀变的岩石；在对所有样品进行显微岩相学观察的基础上，筛选8件进行全岩成分分析，2件进行锆石定年分析，3件进行同位素测试。

图3 马庄山地区花岗岩类的TAS图解(a, 底图据Middlemost, 1994)和K2O-SiO2图解(b,底图据Rickwood, 1989)Fig.3 TAS diagram (a, base map after Middlemost, 1994) and K2O vs. SiO2 diagram (b, base map after Rickwood, 1989) of the granitoids in Mazhuangshan area

3.2 测试方法

全岩主量、微量及稀土元素分析测试在澳实矿物实验室(广州)完成。主量元素测试采用ME-XRF06的分析方法，所用仪器为X-荧光光谱分析仪，型号为PANalytical AXIOS，检测下限为0.01n×10-2，分析精度和准确度优于0.01%；微量元素所用分析方法为ME-ICP61，所用仪器为电感耦合等离子体质谱仪，型号为ELAN 9000，分析精度和准确度：Th、U为0.05×10-6，Cs、Sr、Ta为0.1×10-6，Ba为0.5×10-6，Rb、Hf、Nb为0.2×10-6，Zr为2×10-6，V、Co、Ni、Cr为1×10-6，K、P、Ti为0.01%；稀土元素分析方法为ME-MS81(硼酸锂熔融、等离子质谱定量测试)，所用仪器仍为电感耦合等离子体质谱仪，分析精度和准确度除La、Ce、Y为0.5×10-6外，其余均为0.05×10-6。

Nd同位素组成测试在中国科学院广州地球化学研究所同位素地球化学国家重点实验室完成，样品测试所用仪器为Finnigan Neptune多接收器电感耦合等离子质谱仪(MC-ICP-MS)。同位素分析采用Teflon溶样器，加入HNO3和HF混合溶样，用专用的阳离子交换柱进行分离。详细的分析流程及仪器分析情况见(梁细荣等, 2002)。在本文样品分析过程中，EstonJndi-1标准的143Nd/144Nd测定值0.512087±2(2σ，N=18)。

锆石U-Pb定年在中国地质大学(北京)地质过程和矿产资源重点实验室完成，分析过程中使用的激光剥蚀束斑，直径为32μm，使用He、Ar分别作为其载气和补偿气。年龄外标采用国际标准锆石91500校正，MUD作为同位素监测样品，元素含量的外标和内标分别采用的是NIST610、29Si，通过本次实验测得的标样的结果，在推荐值范围内(Ludwig, 2003)。

锆石Hf同位素分析在武汉上谱分析科技有限责任公司完成，仪器采用多接收质谱仪MC-ICP-MS(Neptune Plus)和相干193nm准分子激光剥蚀系统(GeoLasPro HD)。关键参数：束斑：44μm，能量强度：8mJ/cm2，频率：8Hz，载气：600ml/min，标样推荐值(91500：0.282308，GJ-1：0.282013，TEM: 0.282677)。εHf(t)根据每个测点的锆石U-Pb年龄计算而来，采用的176Lu衰变常数λ=1.867×10-11y(Söderlundetal., 2004)，利用平均大陆壳的176Lu/177Hf=0.015(Griffinetal., 2002)计算锆石Hf同位素地壳模式年龄(tDM2)。

4 测试结果

4.1 主、微量元素

钾长花岗岩SiO2含量为76.35%～78.48%、K2O为2.86%～4.95%、P2O5为0.01%、全碱(Na2O+K2O)含量7.14%～8.75%、K2O/Na2O比值0.67～1.33、Zr/Hf比值介于18～24之间(表1)。在TAS图解上，样品点落在花岗岩区域(图3a)；在SiO2-K2O图解(图3b)中，主要落入高钾钙碱性-钙碱性系列范围。钾长花岗岩的稀土元素总量中等，∑REE=83×10-6～154×10-6、(La/Yb)N值为1.4～6.3、δEu=0.04～0.21，相对富集轻稀土元素、亏损重稀土元素，显示强的负Eu异常(图4a)，表明岩浆曾发生明显的斜长石的分离结晶作用(张宏飞和高山, 2012)。微量元素原始地幔标准化蛛网图显示：钾长花岗岩富集Rb、K、Th和U大离子亲石元素，亏损P、Ta、Nb和Ta高场强元素(图4b)。

英云闪长岩SiO2含量为64.54%～64.75%、K2O为3.63%～4.33%、P2O5为0.20%～0.21%、全碱(Na2O+K2O)含量为6.51%～6.94%、K2O/Na2O比值为0.50～0.91，Zr/Hf比值38～39。在TAS图解上，样品落在花岗闪长岩区域(图3a)；在SiO2-K2O图解上落入高钾钙碱性-钙碱性区域(图3b)。英云闪长岩的稀土元素总量中等，∑REE=153.35×10-6～162.49×10-6，(La/Yb)N值为8.6～10.9，δEu=0.89～0.90(表1)，表明轻重稀土分异中等，相对富集轻稀土元素、亏损重稀土元素，显示弱的负Eu异常(图4a)，表明岩浆没有发生明显的斜长石的分离结晶作用(张宏飞和高山, 2012)。总体富集Rb、K、Th和U大离子亲石元素，亏损P、Ta、Nb和Ti高场强等元素(图4b)。

花岗闪长岩SiO2含量为66.64%～67.86%、K2O为1.53%～4.65%、P2O5为0.02%～0.12%、全碱(Na2O+K2O)为4.72%～6.55%、K2O/Na2O比值0.30～0.43、Zr/Hf比值28～36。在TAS图解上，样品落在花岗闪长岩区域(图3a)；全岩样品中硅高钾中磷，1件为钾玄岩系列，另2件是钙碱性系列(图3b)。花岗闪长岩的稀土元素总量中等，∑REE=76×10-6～193×10-6，(La/Yb)N值为7.7～14.0，δEu=0.47～0.97(表1)，表明轻重稀土分异中等，显示弱的负Eu异常(图4a)。微量元素特征与英云闪长岩相似(图4b)。

表1 马庄山地区花岗岩类主量元素(wt%)、稀土、微量元素(×10-6)测试结果

图4 马庄山地区花岗岩类球粒陨石标准化稀土元素配分型式(a, 标准化值据Boynton, 1984)和原始地幔标准化微量元素蜘蛛图(b, 标准化值据Sun and McDonough, 1989)Fig.4 Chondrite-normalized REE patterns (a, normalization values after Boynton, 1984) and primitive mantle-normalized trace element spider diagrams (b, normalization values after Sun and McDonough, 1989) of the granitoids

图5 花岗岩类岩石中锆石阴极发光图像Fig.5 Cathodoluminescence images of the zircon from granitoids

图6 花岗岩类岩石中锆石的球粒陨石标准化稀土元素配分图(标准化值据Boynton, 1984)Fig.6 Chondrite-normalized REE patterns of the zircon from granitoids (normalized after Boynton, 1984)

4.2 锆石U-Pb年龄

选取2件岩石样品进行了锆石U-Pb定年分析。所测定锆石多呈透明长柱或短柱状，自形程度较好，粒径多为80～120μm，阴极发光图像上具清晰的震荡环带(图5)。锆石具有较高的Th(43×10-6～634×10-6)、U(154×10-6～758×10-6)含量和Th/U比值(0.28～0.86)(表2)，稀土元素配分曲线具有轻稀土亏损、重稀土富集、明显Ce正异常、Eu负异常的特点(图6)；显示岩浆锆石的特征(Belousovaetal., 2002)。

图7 马庄山地区花岗岩类LA-ICP-MS U-Pb年龄一致曲线(a、c)和加权平均年龄谱图(b、d)Fig.7 U-Pb concordant age (a, c) and weighted mean zircon 206Pb/238U age (b, d) of the granitoids at Mazhuangshan

钾长花岗岩(YF15-2)8颗锆石颗粒的LA-ICP-MS U-Pb年龄测定表明，其U-Pb一致年龄为317.7±1.0Ma(MSWD=0.031；图7a)，206Pb/238U加权平均年龄值为318.2±4.0Ma(MSWD=1.2；图7b)。

花岗闪长岩(MZS7-6)22颗锆石颗粒的LA-ICP-MS U-Pb年龄测定显示，其U-Pb一致年龄为320.2±0.8Ma(MSWD=1.9；图7c)，206Pb/238U加权平均年龄值为319.3±2.6Ma(MSWD=0.045；图7d)。

4.3 Nd-Hf同位素

本文测试并统计研究区附近的全岩Nd同位素测试样品共10件，统计结果见表3。钾长花岗岩的εNd(t)在-5.30～-4.24之间，(143Nd/144Nd)i值为0.511958～0.512013。花岗闪长岩的εNd(t)为-5.31，(143Nd/144Nd)i值为0.511968。石英斑岩的εNd(t)在-2.27～1.82之间，(143Nd/144Nd)i值为0.512135～0.512326。二长花岗岩的εNd(t)在-2.0～-1.7之间，(143Nd/144Nd)i值为0.512113～0.512128。

花岗闪长岩锆石Hf同位素测点位于LA-ICP-MS U-Pb年龄测点旁(图5b)，Hf同位素共测试12个点(MZS7-6)。测试结果见表4，176Lu/177Hf比值介于0.000804～0.001816，176Hf/177Hf比值在0.282489～0.282590，εHf(t)介于-3.6～-0.2，平均-1.1，tDM介于947～1086Ma，tDM2为1318～1537Ma。

5 讨论

5.1 岩石成因

岩浆岩类型的识别对于岩浆源区、构造环境识别有重要的意义，其成因类型可划分为I型、S型和A型(Chappell and White, 1974, 1991; Collinsetal., 1982)。I型花岗岩被认为是未风化的火成岩经过部分熔融、结晶形成的，Na2O/K2O值一般小于1，特征矿物为角闪石(Chappelletal., 2012)；S型花岗岩被认为是壳源沉积物经部分熔融、结晶形成的，Na2O/K2O值一般大于1，特征性矿物为堇青石、石榴石和夕线石(Chappell and White, 1974)；A型花岗岩是产于稳定大陆板块和裂谷内部的碱性、无水的花岗岩，特征矿物是碱性暗色矿物(如钠闪石；Chappell and White, 1974; Mille, 1985)。

研究区英云闪长岩在显微镜下明显可见角闪石。英云闪长岩和花岗闪长岩的Na2O/K2O比值为0.3～0.9。在(Na2O+K2O)/CaO-Zr+Nb+Ce+Y图解(图8a)中，英云闪长岩和花岗闪长岩的投影点均落入未分异花岗岩区；在10000Ga/Al-Zr和10000Ga/Al-Nb图解(图8b, c)中， 样品点均落在I或S型花岗岩范围内。岩浆分异演化过程中，I型花岗岩的P2O5含量具有随着SiO2含量的增加而递减的特征，S型花岗岩则具有相反的趋势，因此，P2O5可以用来区分低分异的I型和S型花岗岩(Lietal., 2006, 2007; Wuetal., 2003)。在P2O5-SiO2关系图解(图8d)中，P2O5与SiO2相呈现负相关关系，与S型花岗岩的负相关性有明显的区别；在Rb-Th和Rb-Y的图解(图8e, f)中，Y、Th与Rb均呈现正相关的关系，与I型花岗岩的特征一致(Chappell, 1999; Green and Watson, 1982; Lietal., 2007)。综上，认为研究区英云闪长岩和花岗闪长岩应属I型花岗岩，这也与区内广泛发育的二长花岗岩的类型相一致(Zhangetal., 2017)。

表3 马庄山地区花岗岩类全岩Nd同位素组成

表4 马庄山地区花岗闪长岩中锆石的Hf同位素组成

图8 马庄山地区花岗岩类岩石的成因判别图(底图据Whalen et al., 1987)(a) (Na2O+K2O)/CaO与Zr+Nb+Ce+Y图；(b) Zr与10000Ga/Al花岗岩分类图；(c) Nb与10000Ga/Al花岗岩分类图；(d) P2O5-SiO2图解；(e) Rb-Th图解；(f) Rb-Y图解Fig.8 Discrimination diagrams of the granitoids at Mazhuangshan area (base map after Whalen et al., 1987)(a) (Na2O+K2O)/CaO vs. Zr+Nb+Ce+Y diagram; Zr (b) and Nb (c) vs. 10000 Ga/Al diagrams; (d) P2O5 vs. SiO2 diagram; Th (e) and Y(f) vs. Rb diagrams

图9 马庄山地区花岗岩类岩石t-εNd(t)(a)和t-εHf(t)(b)图解除本文之外的其他数据引自聂凤军等(2002)、洪大卫和谢锡林(2000)、Jahn et al. (2000)、Lei et al. (2013)、Song et al. (2013b)、Zhang et al. (2017)及王琦崧等(2019)Fig.9 Diagrams of the age vs. εNd(t) (a) and age vs. εHf(t) (b) for the granitoids in Mazhuangshan areaThe data of other intrusions from Nie et al. (2002), Hong and Xie (2000), Jahn et al. (2000), Lei et al. (2013), Song et al. (2013b), Zhang et al. (2017) and Wang et al. (2019)

相比较而言，区内钾长花岗岩因其较高的10000Ga/Al值，样品总体落入了A型花岗岩和靠近A型花岗岩的区域内(图8b, c)。考虑到A型花岗岩往往与SiO2>72%的I型或S型高分异花岗岩在地球化学特征上表现出一定的相似性(Whalenetal., 1987; Wuetal., 2003)，而钾长花岗岩与花岗闪长岩是同时代的产物(其为年龄分别为317Ma和320Ma左右)(图8a, b)，具有相同的矿物组合和类似的地球化学特征，因此，推测钾长花岗岩为高分异I型花岗岩。

5.2 岩浆来源

中亚造山带大部分显生宙花岗岩以正的εNd(t)值与年轻的Sm-Nd模式年龄为主(Jahnetal., 2000)。研究区内大部分花岗岩类样品的εNd(t)值小于0，反映了前寒武纪古陆的残余(图9a)，或者说岩浆形成过程中有老地壳的混入，与研究区广泛出露的早中元古界所处的浅海相沉积环境吻合(周济元等, 2002)。与此同时，研究区少部分花岗岩类样品的εNd(t)值>0，反映了岩石可能有一定的地幔物质来源。总体而言，研究区花岗岩类的εNd(t)表现了壳幔混合的特征。

Hf同位素具有较高的封闭温度，不会因后期的部分熔融、分离结晶而发生变化，故Hf同位素可以用来揭示岩石源区特征(Griffinetal., 2000; 吴福元等, 2007)。锆石Hf同位素分析目前被广泛用于反应花岗岩的源岩性质和源区特征。一般来说，正的εHf(t)值被解释为其源岩可能来自新生地壳或亏损地幔(Jahnetal., 2000; Vervoortetal., 2000)，负的εHf(t)值指示其源岩为古老地壳物质(Kinny and Maas, 2003)。研究区花岗闪长岩的εHf(t)介于-3.6～-0.2，tDM2为1318～1537Ma。在εHf(t)-t图解中，大部分样品的投点落在亏损地幔(DM)和球粒陨石(CHUR)演化线之间，少部分样品则落在球粒陨石演化线以下，暗示其源岩以古老地壳为主，并有少量新生地壳或亏损地幔参与(图9b)，这也与研究区Nd同位素特征相一致。锆石Hf二阶段模式年龄(tDM2)指示亏损地幔熔融形成玄武质下地壳时代，即模式年龄代表了壳幔分异时代(吴福元等, 2007)。研究区花岗闪长岩锆石Hf二阶段模式年龄(tDM2)为1321～1537Ma(表4)，与岩体的实际侵位年龄(～320Ma)差别较大，代表了旱山弧微古陆在中元古代早期有一次古老地壳增生。

图10 旱山弧花岗质岩石锆石Hf模式年龄频率直方图Fig.10 Zircon Hf isotopic diagram from the granitoids in the Hanshan arc

为了对比研究，同时收集北山地区旱山弧内、不同时代的花岗质岩石的Hf同位素组成，分析旱山弧地壳演化史。从西到东分别选取红柳井(95°E)、马庄山(96°E)、石板井(98°E)的中元古代至晚古生代花岗质片麻岩和花岗岩类岩石(Heetal., 2014, 2015; Wangetal., 2014; Gaoetal., 2015)，它们的锆石Hf模式年龄介于1.3～2.3Ga之间，显示1.5～1.6Ga和1.9～2.0Ga两个峰值(图10)。前人的研究表明，旱山弧与北山锆石Hf模式年龄都具有1.9～2.0Ga峰值(Heetal., 2018)，这一时间与元古代全球陆壳增长的重大地质事件同步(Kröneretal., 2014, 2017)。这进一步印证了，研究区的古地壳可能是在古元古代(约2.0Ga)形成，之后经历了新一期的中元古代(约1.5Ga)地壳生长。

5.3 岩浆侵位时代与构造环境

图11 马庄山地区花岗岩Hf-Rb/30-Ta×3构造环境判别图解(底图据Rollinson, 1993)Fig.11 The Hf-Rb/30-Ta×3 discriminant diagram for granites (base map after Rollinson, 1993)

图12 北山造山带北带奥陶纪-石炭纪构造演化示意图(据Ao et al., 2016修改)Fig.12 Schematic diagram showing the tectonic evolution of the Northern Beishan Orogen from Ordovician to Carboniferous (modified after Ao et al., 2016)

前人曾对研究区及附近的花岗岩类岩石开展了一系列的同位素年代学工作。Zhangetal. (2017)获得研究区东北侧明水地区的二长花岗岩锆石LA-ICP-MS锆石U-Pb一致年龄为328±2Ma；王琦崧等(2019)获得马庄山地区的石英斑岩的锆石LA-ICP-MS锆石U-Pb一致年龄分别为315.4±0.6Ma；本次工作首次获得马庄山地区的钾长花岗岩和花岗闪长岩的锆石LA-ICP-MS锆石U-Pb一致年龄分别为317.7±1.0Ma 和320.2±0.8Ma。综合上述研究结果，研究区花岗岩类岩石侵位时代在315～328Ma，反应了该区岩浆岩侵位时代于早石炭世-晚石炭世。

前人对北山造山带晚古生代的构造环境的认识存在分歧。一种观点认为，晚古生代火成岩与俯冲-增生过程有关(Xiaoetal., 2010; Aoetal., 2016; Han and Zhao, 2018)；而另一种观点认为，晚古生代的侵入岩形成于伸展环境(Zhengetal., 2013; Zhangetal., 2017)。从本文研究区情况看，北山造山带马庄山地区的岩石组合为英云闪长岩、花岗闪长岩和钾长花岗岩。上述部分花岗岩类岩石中含有富水矿物角闪石、黑云母，说明岩浆含水，显示与洋壳俯冲作用有关的弧火山组合特征(肖庆辉等, 2009)。与此同时，岩石相对富集大离子亲石元素，亏损高场强元素(图4b)，Ti-Nb-Ta的负异常，显示了岩石处于俯冲区弧火山岩特征(Pearceetal., 1984)。在花岗岩Hf-Rb/30-Ta×3构造环境判别图上(图11)，英云闪长岩、花岗闪长岩和钾长花岗岩样品点均落入火山弧区，也反映了研究区晚石炭世为火山弧环境。

结合整个区域构造演化情况以及研究区花岗岩类构造环境背景，将北山造山带内奥陶纪-石炭纪期间旱山弧的构造演化总结如下(图12)：奥陶纪-泥盆纪时期，双鹰山弧与马鬃山弧之间的月牙山洋闭合，月牙山蛇绿岩形成，双鹰山弧与马鬃山弧拼贴在一起(图12a)；早石炭世时期，古亚洲洋向东天山和北山下俯冲，俯冲板片后撤使得东天山-北山复合造山带北部进入弧后拉分盆地时期(Xiaoetal., 2018)，红石山洋和星星峡-石板井洋盆可能就是古亚洲洋板片后撤拉张的产物(张元元和郭召杰, 2008; Zhengetal., 2013; 王国强等, 2014; Xiaoetal., 2018)。与此同时，旱山弧南侧的芨芨台子-小黄山洋向南侧俯冲(图12b)；晚石炭世时期，以芨芨台子-小黄山洋代表的洋向南俯冲角度逐渐减小，直至闭合。与此同时，古亚洲洋向南俯冲，在约317～320Ma，研究区处于火山弧环境，下地壳物质的部分熔融形成了钾长花岗岩和花岗闪长岩的初始岩浆，侵位至地壳浅部产出(图12c)。晚石炭世末期，古大洋逐渐闭合，北山研究区开始进入后碰撞造山阶段。

6 结论

(1)马庄山地区的花岗岩类岩石包括花岗闪长岩、钾长花岗岩、英云闪长岩，总体属于富硅、钾，低铁、磷、钛的高钾钙碱性-钙碱性系列，英云闪长岩、花岗闪长岩显示低分异I型花岗岩的特征，而钾长花岗岩显示高分异I型花岗岩的特征。

(2)马庄山地区钾长花岗岩与花岗闪长岩成岩年龄分别为317.7±1.0Ma、320.2±0.8Ma，形成于古亚洲洋俯冲背景下的弧环境。

(3)研究区花岗岩类初始εNd(t)在-5.31～1.82。花岗闪长岩锆石176Hf/177Hf在0.282489～0.282590之间，εHf(t)介于-3.6～-0.2之间，tDM2为1318～1537Ma，认为其可能来源于下地壳物质的部分熔融，原始岩浆在上侵过程中有部分地幔物质的加入。

致谢野外工作得到了西部矿业哈密金矿、新疆有色704队的大力支持；实验工作得到了中国科学院广州地球化学研究所肖兵博士、中国地质大学(北京)地质过程与矿产资源国家重点实验室相鹏老师的帮助和支持；两位审稿人为本文提出了建设性修改建议；在此一并表示感谢。