胶东西北部北截岩体岩石成因：锆石U-Pb年龄、岩石地球化学与Sr-Nd-Pb同位素制约

宋英昕 于学峰 李大鹏,3 耿科 尉鹏飞,3 左晓敏 王秀凤

1. 自然资源部金矿成矿过程与资源利用重点实验室，山东省金属矿产成矿地质过程与资源利用重点实验室，山东省地质科学研究院，济南2500132. 山东离子探针中心，济南 2500133. 中国地质大学地质过程与矿产资源国家重点实验室，北京 1000834. 山东省地质博物馆，济南 2500141.

胶东半岛西北部是我国重要的金矿床密集区，三山岛、焦家、新城、玲珑、台上等超大型金矿床分布于该区(图1)。区内的晚侏罗世玲珑型花岗岩和早白垩世郭家岭型花岗岩赋存了胶东95%以上的金资源储量，是胶东金矿最主要的赋矿围岩(Yangetal., 2006, 2016a; Zhai and Santosh, 2013; Goldfarb and Santosh, 2014; Yanetal., 2014; Dengetal., 2015a, b, 2018)。赋矿的郭家岭型花岗岩一直是国内外学者研究的热点(Dengetal., 2003a, 2011; Chenetal., 2005; Lietal., 2006; Fanetal., 2007; Yangetal., 2008, 2009, 2014; Tanetal., 2012; Yang and Badal, 2013; 杨立强等，2014, 2020)。然而，目前对于郭家岭型花岗岩的成因或其形成的大地构造背景及其与金成矿的关系，尚存在较多的不同意见。在花岗岩类成因方面，一些学者将其归为埃达克岩石，认为俯冲增厚的大陆地壳部分熔融形成郭家岭花岗岩(Houetal., 2007; Zhangetal., 2010; Deng and Wang, 2016; Dengetal., 2017a, b)；许多研究者认为其属于高Ba-Sr花岗岩(刘跃等，2014；王中亮等，2014)，是壳幔相互作用的结果(Yangetal., 2012; Yang and Badal, 2013; 刘跃等，2014；王中亮等，2014)；此外，杨进辉等(2003)则认为郭家岭型花岗岩地球化学特征类似于年轻的TTG(<3.0Ga)和Na质花岗岩，是由下地壳镁铁质岩石脱水部分熔融作用形成的。在花岗岩类与金成矿的关系方面，许多研究者认为金矿床为郭家岭型花岗岩岩浆期后热液矿床(Wangetal., 1998；李士先等，2007；罗贤冬等，2014)；也有研究者提出，胶东金矿成矿主要与早白垩世中期的伟德山型花岗岩及其同期的脉岩、火山岩有关(刘辅臣等，1984；孙景贵等，2000；罗振宽等，2001；宋明春，2014)，早白垩世中期大规模岩浆作用引发了强烈的流体活动和成矿物质活化(Songetal., 2014, 2015; 宋明春等, 2014)。可见，进一步研究郭家岭型花岗岩岩的成因，对深化理解胶东早白垩世大规模金成矿的大地构造背景及岩浆活动与金矿化的关系具有十分重要的意义。

前人对归属于郭家岭型花岗岩的郭家岭岩体、丛家岩体、上庄岩体、新城岩体等的岩石地球化学、同位素地球化学和同位素年代学做过详细系统的研究(杨进辉等，2003；Houetal., 2007; Zhangetal., 2010；Yangetal., 2012, 2006; Yang and Badal, 2013；Wangetal., 2014; 罗贤冬等，2014；刘跃等，2014；王中亮等，2014)，而对同属于郭家岭型花岗岩的北截岩体仅进行过同位素年龄测试工作(Wangetal., 1998；耿科等，2015)。本文通过对北截岩体进行详细的野外地质调查和显微镜下观察，区分出细粒花岗闪长岩、含斑细粒花岗闪长岩、斑状中粒花岗闪长岩、斑状粗粒花岗闪长岩及花岗伟晶岩脉5种侵入岩，并对其进行了系统的矿物化学、主微量元素地球化学、锆石U-Pb年龄和Sr-Nd-Pb同位素研究，确定了北截岩体的地球化学类型，进一步约束了郭家岭型花岗岩的形成时代及岩石成因，揭示了其成岩构造环境，探讨了岩浆活动与金成矿的关系。

1 区域地质概况

胶西北位于华北克拉通东南缘胶东半岛西北部(图1；Dengetal., 2003a, b; Tanetal., 2008)，苏鲁超高压变质带北部，郯庐断裂带东侧，区内广泛发育NNE向展布的郯庐断裂系的各级构造。其结晶基底岩系主要由中-新太古代TTG质花岗片麻岩、新太古代胶东岩群、古元古代粉子山群和荆山群，及新元古代蓬莱群变质岩组成(Tametal., 2011；Lietal., 2012)。区域中生代岩浆活动十分发育，主要为晚侏罗世玲珑型花岗岩、早白垩世早期郭家岭型花岗岩和早白垩世中期伟德山型花岗岩(艾山花岗岩体)(图1；Sunetal., 2007; Maetal., 2014; Wangetal., 2014)。玲珑型花岗岩岩性主要为黑云母花岗岩，其锆石U-Pb年龄为166～149Ma(Jiangetal., 2012; Yangetal., 2012)，被认为是加厚的下地壳(主要是TTG质花岗片麻岩和胶东岩群)部分熔融的产物(Zhangetal., 2010; Jiangetal., 2012)。郭家岭型花岗岩侵入于玲珑型花岗岩体中，由郭家岭、丛家、北截、上庄、新城和三山岛等岩体组成(图1)，总面积约514km2，各岩体呈北东东向不连续带状分布，自东向西岩体面积由大变小。岩性主要为二长闪长岩、石英二长岩、花岗闪长岩和二长花岗岩，其锆石U-Pb年龄为132～126Ma(Houetal., 2007; Yangetal., 2012)。

图2 北截岩体地质图(据耿科等，2015修改)Fig.2 Geological map of Beijie pluton (modified after Geng et al., 2015)

北截岩体位于丛家岩体西南侧，形态为近东西向的椭圆状，出露面积约15km2，主要由花岗闪长岩组成，呈岩株状侵入玲珑岩体中(图2，宋明春和王沛成，2003)。伟德山型花岗岩的艾山岩体侵入郭家岭型花岗岩和玲珑型花岗岩中，呈近南北走向将郭家岭岩体切割为东西两部分，主要由斑状中粗粒二长花岗岩组成，含少量花岗闪长岩，形成时间为116±1Ma(Gossetal., 2010)，被认为是来源于幔源基性岩浆与壳源酸性岩浆相互作用产生的基性岩浆底侵引起的地壳底部岩石部分熔融形成的(Gossetal., 2010；杨宽等，2012)。在研究区东部的藏格庄盆地及研究区以南的胶莱盆地中，发育有与早白垩世花岗岩同时代(130～110Ma)的青山群火山岩(Fanetal., 2001; Liuetal., 2009)，由玄武岩、玄武粗安岩、粗面岩、安山岩、流纹岩等组成(宋明春和王沛成，2003)，其成因是交代的富集岩石圈地幔部分熔融形成(Fanetal., 2001)。

2 样品采集与岩相学特征

以往区域地质调查将北截岩体划归为郭家岭序列大草屋单元，认为这一岩体是由斑状粗中粒含黑云花岗闪长岩组成(宋明春和王沛成，2003)的单一岩性侵入体。我们通过详细的野外观察发现，北截岩体岩性组成并非单一，是由多种岩性侵入体组成的复式岩体。野外宏观上，岩石呈灰白色至浅肉红色，根据粒度大小和有无斑晶，分为细粒花岗闪长岩、含斑细粒花岗闪长岩、斑状中粒花岗闪长岩、斑状粗粒花岗闪长岩和花岗伟晶岩脉。本文在认真观察识别各种岩性侵入体的基础上，于招远市大刘家村北400m处自然露头和招远市北朱家村西北700m处人工采石坑，采集了新鲜的花岗闪长岩和花岗伟晶岩脉样品(图2)。

花岗闪长岩主要造岩矿物为斜长石(45%～50%)、钾长石(20%～22%)、石英(20%～25%)和少量黑云母(5%～7%)、角闪石(1%～5%)，副矿物有锆石、磷灰石、磁铁矿等。具似斑状结构的花岗闪长岩，斑晶主要为钾长石，粒径多1～3cm，偶见石英斑晶；基质由长英质矿物和暗色矿物组成。斜长石呈自形-半自形板状，聚片双晶较为发育，部分大颗粒可见环带构造。钾长石呈半自形-他形粒状，少数大颗粒可见卡式双晶。石英呈他形粒状充填在长石颗粒间，粒度大小不一，常呈集合体状分布，与长石一起形成粗-细粒花岗结构(图3)。黑云母呈细小碎片、长条状，角闪石呈柱状分布于石英与长石的边缘。其中，斑状粗粒花岗闪长岩(BJ1，图3a，b)，呈似斑状结构，斑晶含量约20%，粒径10～30mm，基质呈粗粒半自形粒状结构，粒径0.5～1mm；斑状中粒花岗闪长岩(BJ3，图3e，f)，呈似斑状结构，斑晶含量约20%，粒径10～30mm，基质呈中粒半自形粒状结构，粒径0.25～0.5mm；细粒花岗闪长岩(BJ4，图3g，h)，呈细粒半自形粒状结构，主要矿物粒径0.05～0.25mm；含斑细粒花岗闪长岩(BJ5，图3i，j)，呈似斑状结构，斑晶含量约5%，粒径10mm左右，基质呈细粒半自形粒状结构，粒径0.05～0.25mm。

图3 北截岩体手标本及显微镜照片(a)浅肉红色斑状粗粒花岗闪长岩；(b)钾长石斑晶；(c)浅肉红色花岗伟晶岩；(d)石英呈集合体状分布；(e)浅肉红色斑状中粒花岗闪长岩；(f)中粒半自形粒状结构；(g)浅灰白色细粒斑状花岗闪长岩；(h)细粒半自形粒状结构；(i)浅肉红色-灰白色含斑细粒花岗闪长岩；(j)环带状斜长石斑晶. Q-石英；Kfs-钾长石；Pl-斜长石；Bt-黑云母Fig.3 Petrological characteristics of Beijie pluton(a) light flesh-pink porphyritic megagrained granodiorite; (b) K-feldspar phenocysts; (c) light flesh-pink granite pegmatite; (d) quartz concentrated in distribution; (e) light flesh-pink porphyritic medium-grained granodiorite; (f) medium-grained hypidiomorphic texture; (g) light off-white porphyritic fine-grained granodiorite; (h) fine-grained hypidiomorphic texture; (i) light flesh-pink to off-white porphyritic fine-grained granodiorite; (j) plagioclase phenocryst with zone. Q-quartz; Kfs- K-feldspar; Pl-plagioclase; Bt-biotite

花岗伟晶岩呈脉状侵入花岗闪长岩中，具花岗伟晶结构(BJ2，图3c，d)，主要由斜长石(45%～50%)、钾长石(20%～22%)和石英(20%～25%)组成，含有少量角闪石(1%～5%)和副矿物磁铁矿等。钾长石呈半自形-他形粒状，晶体粒度10～30mm，斜长石呈自形-半自形板状，聚片双晶较为发育。石英呈他形粒状充填在长石颗粒间，粒度大小不一，常呈集合体状分布。

3 分析方法

3.1 电子探针分析

电子探针测试工作在山东省地质科学研究院自然资源部金矿成矿过程与资源利用重点实验室完成，仪器型号为日本电子(JEOL)JXA-8230，波谱分析所用加速15kV，电流1×10-8A，束斑直径1～10μm。所用标准样品均为加拿大Astimex标样。

3.2 主量和微量元素分析

主量元素、微量元素的分析测试工作分别是在自然资源部济南矿产资源监督检测中心和核工业北京地质研究院分析测试研究中心完成的。主量元素分析使用Philips PW2404型X荧光光谱仪(XRF)完成，分析精度优于1%，其中FeO采用化学容量法测定， Fe2O3T为荧光光谱仪(XRF)测试所得的全铁的量，FeO的含量是通过湿化学分析方法得到的；微量元素分析使用Finnigan MAT Element I型电感耦合等离子体质谱仪(ICP-MS)完成，分析精度多小于3%，温度20℃，湿度30%。测试过程中分别使用国际标准进行监控。

3.3 锆石LA-ICP-MS年龄测定

锆石的分选在廊坊市峰泽源岩矿检测技术实验室完成的，先将岩样粉碎至300μm，再用浮选和电磁选方法分离，最后在双目镜下进行单颗粒锆石挑选。靶样、锆石透射光、反射光显微照相以及锆石阴极发光(CL)显微照相均由中国地质大学(北京)科学研究院制作和完成的。锆石LA-ICP-MS U-Pb测试工作在北京锆年领航科技有限公司的激光等离子体质谱实验室完成，分析仪器为美国New Wave Research Inc公司生产的193nm激光剥蚀进样系统(UP193SS)和美国Agilent科技有限公司生产的Agilent7500a型四级杆等离子体质谱仪联合构成的激光等离子体质谱仪(LA-ICP-MS)。工作条件为：测试点束斑直径30μm，激光频率10Hz，预剥蚀时间5s，剥蚀采样时间45s，ICP-MS辅助气为Ar，流速1.13L/min，数据处理软件为Glitter4.4.1。年龄计算以标准锆石91500为外标进行同位素比值校正，标准锆石TEMORA和Qinghu为监控盲样；元素含量以国际标样NIST610为外标，Si为内标计算NIST612和NIST614为监控盲样。204Pb标准值和校正参考Anderson(1984)。实验详细方法参考Songetal.(2010)和Wangetal.(2012)。

3.4 Sr-Nd-Pb同位素分析

Sr-Nd-Pb同位素测试在北京科荟测试技术有限公司完成。使用仪器为Finnigan MAT262多通道质谱仪(LA-MC-ICPMS)。Sr和Nd同位素比值测定分别采用87Sr/88Sr=0.1194、146Nd/144Nd=0.7219进行标准化，详细的实验流程见Zhangetal.(2002)。铅同位素分析采用HF酸(2% HNO3)在高温高压条件下将粉末样品完全溶解，在装有交换柱上分离纯化Pb样品，质量分馏校正系数为每质量单位1‰，该校正系数源于对铅标准物质NBS981的大量测试所获得的数据。实验中国际标样分析结果分别为NBS 98787Sr/86Sr=0.710248±13(2SD，n=35)；GSB Nd143Nd/144Nd=0.512185±8(2SD，n=68)；CAGS Pb206Pb/204Pb=17.9708±0.0009(2SD，n=22)，207Pb/204Pb=15.5626±0.0010(2SD，n=22)，208Pb/204Pb=38.4059±0.0024(2SD，n=22)。

4 分析结果

4.1 矿物化学

对北截岩体斑状粗粒花岗岩(BJ1-1)和细粒花岗闪长岩(BJ4-1)中的黑云母进行了电子探针分析(表1)，结果表明这两组岩石中黑云母成分没有显著差别，均具有富MgO、FeOT、高Al2O3，低TiO2的特点。斑状粗粒花岗岩中黑云母MF值为0.44～0.47，Fe/(Fe+Mg)值为0.48～0.52；细粒花岗闪长岩中黑云母MF值为0.37～0.39，Fe/(Fe+Mg)值为0.56～0.60。

4.2 主量和微量元素

选择北截岩体中五组共20件样品进行岩石地球化学分析，五组样品分别为斑状粗粒花岗闪长岩、斑状中粒花岗闪长岩、细粒花岗闪长岩、含斑细粒花岗闪长岩和花岗伟晶岩，其主量、微量和稀土元素组成见表2。

4.2.1 主量元素

所有花岗闪长岩样品SiO2含量为70.82%～74.40%，K2O含量为1.39%～4.91%，Na2O含量为4.40%～5.74%，总碱值(Na2O+K2O)为6.92%～9.31%。在TAS图解中，样品投点于接近花岗闪长岩的花岗岩范围(图4a)；在SiO2-K2O图解(图4b)上，主要投点于高钾钙碱性-钙碱性系列；岩石的镁指数(Mg#)为21.93～30.11，铝指数(A/CNK)为1.43～1.60，属于过铝质花岗岩(图4c)。

图4 北截岩体TAS图解(a，据Middlemost, 1994)、SiO2-K2O图解(b，据Morrison, 1980)和A/CNK-A/NK图解(c，据Middlemost, 1994)图5、图8、图10、图11、图13图例同此图Fig.4 TAS diagram (a, after Middlemost, 1994), SiO2 vs. K2O diagram (b, after Morrison, 1980) and A/CNK vs. A/NK diagram (c, after Middlemost, 1994) for the Beijie plutonThe legends in Fig.5, Fig.8, Fig.10, Fig.11, Fig.13 are the same as in this figure

图5 北截岩体花岗闪长岩和伟晶岩主量元素与SiO2协变图Fig.5 Diagrams of SiO2 vs. predominant oxides of the granodiorite and pegmatite from Beijie pluton

花岗伟晶岩SiO2含量相对较高，为73.27%～76.80%，K2O含量为4.19%～7.20%，钾含量较高，Na2O含量为3.27%～4.42%，Na2O+K2O值为8.14%～10.47%。在TAS图解中，样品落入花岗岩范围(图4a)；在SiO2-K2O图解上属于高钾钙碱性-钾玄岩系列(图4b)；镁指数(Mg#)为6.69～14.09，铝指数(A/CNK)为1.29～1.39，属于过铝质花岗岩(图4c)。

在Harker图解上(图5)，花岗闪长岩和伟晶岩的SiO2与其余的氧化物之间具有良好的相关性，SiO2含量与TiO2、Fe2O3T、MgO、CaO、P2O5和MnO含量呈负相关，与K2O含量呈正相关，而且由含斑细粒花岗闪长岩-斑状中粒花岗闪长岩-斑状粗粒花岗闪长岩-细粒花岗闪长岩-花岗伟晶岩呈连续演化关系。

图6 北截岩体球粒陨石标准化稀土元素配分图(a)和原始地幔标准化微量元素蛛网图(b)(标准化值据Sun and McDonough, 1989)埃达克岩石数据引自Martin(1999)Fig.6 REE chondrite-normalized patterns (a) and primitive mantle-normalized spider diagrams (b) of the granodiorite and pegmatite from Beijie pluton (normalization values after Sun and McDonough, 1989)The data for adakitic granitoids from Martin(1999)

图7 北截岩体花岗闪长岩和伟晶岩锆石阴极发光图像(a、b)和锆石U-Pb年龄谐和图(c、d)Fig.7 CL images (a, b) and U-Pb concordia diagrams (c, d) of zircons for the granodiorite and pegmatite from Beijie pluton

4.2.2 微量元素与稀土元素

所有花岗闪长岩的样品的微量元素和稀土元素特征都高度一致。其∑REE为25.57×10-6～85.41×10-6，(La/Yb)N=18.39～45.76，LREE/HREE比值为13.13～21.40，指示轻、重稀土元素发生了强烈分异；δEu为1.16～2.61，具有较为明显的正铕异常；δCe为0.79～1.33，具有较弱的铈负异常到较弱的铈正异常；稀土配分曲线为右倾型式(图6a)，LREE高度富集，HREE相对亏损且曲线较平坦。在原始地幔标准化微量元素蛛网图上(图6b)，花岗闪长岩的微量元素显示富集Rb、Sr、Ba、Pb等大离子亲石元素(LILE)，亏损Nd、Ta、Zr、Hf等高场强元素(HFSE)，与郭家岭岩体花岗闪长岩和新城岩体二长花岗岩的特征相似(杨进辉等，2003；王中亮等，2014)。

花岗伟晶岩与花岗闪长岩的稀土特征略有不同，∑REE为7.60×10-6～26.98×10-6，(La/Yb)N=1.58～6.09，指示轻、重稀土元素分馏较弱；LREE/HREE 比值为2.61～7.21，呈LREE一般富集的右倾分布型式(图6a)；δEu为0.94～7.54，具有弱的负铕异常到强的正铕异常，其中一个样品的正铕异常可达7.54；δCe为1.01～1.31，具有弱的铈正异常。花岗伟晶岩的其他微量元素特征与花岗闪长岩相似，在原始地幔标准化微量元素蛛网图上(图6b)，花岗伟晶岩的微量元素显示富集Rb、Sr、Ba、Pb等大离子亲石元素(LILE)，其中高度富集Pb，亏损Nd、Ta、Zr、Hf等高场强元素(HFSE)。

4.3 锆石 U-Pb年龄

对北截岩体斑状粗粒花岗闪长岩(BJ1-1)和花岗伟晶岩(BJ2-1)2件样品进行锆石U-Pb年龄测定。锆石多为浅灰色，自形柱状、短柱状，晶形完整，表面光滑。阴极发光照片显示(图7a，b)，绝大多数锆石具有振荡环带，属于岩浆结晶期锆石。大多数锆石都具有核幔结构，反映其继承性生长的特点。样品中锆石颗粒较大，内部结构清晰，环带规整，长宽比1:1.5～1:5，符合岩浆锆石的特征。内核多为浑圆状，浅白色，表明较老。部分锆石结晶后又经受过流体作用，具有一层暗色边缘。U含量较低，为42×10-6～1388×10-6，Th含量为26×10-6～903×10-6，Th/U比值为0.04～1.17，变化不大且接近于1 (表3)，也反映这些锆石均属岩浆成因。

表5 北截岩体花岗闪长岩和花岗伟晶岩Pb同位素组成

图8 北截岩体La/Sm-La(a)和Rb/Nb-Rb/Zr(b)相关图解 (底图据谭俊等，2006)Fig.8 La/Sm vs. La (a) and Rb/Nb vs. Rb/Zr (b) diagrams for the granodiorite and pegmatite of the Beijie pluton (base map after Tan et al., 2006)

BJ1-1样品的16个测点取得了较为一致的206Pb/238U 年龄，波动范围为128.9～133.8Ma，加权平均值为131.53±0.86Ma(MSWD=1.2)(表3、图7c)。

BJ2-1样品的10个测点取得了较为一致的206Pb/238U年龄，波动范围为125.5～130.4Ma，加权平均值为127.4±1.3Ma(MSWD=1.7)(表3、图7d)。

4.4 Sr-Nd-Pb同位素

对北截岩体中的斑状粗粒花岗闪长岩(BJ1)、花岗伟晶岩脉(BJ2)、斑状中粒花岗闪长岩(BJ3)、细粒花岗闪长岩(BJ4)、含斑细粒花岗闪长岩(BJ5)进行Sr、Nd、Pb同位素测试分析，实验结果列于表4。

结果显示，北截岩体5种样品的Sr-Nd同位素特征比较相似，显示了同源岩浆的特性。花岗闪长岩的87Sr/86Sr值(0.710883～0.712100)略低于花岗伟晶岩(0.713065～0.713373)，εNd(t)值(-16.57～-11.17)略高于花岗伟晶岩(-18.53～-18.06)；二者的143Nd/144Nd值极为相似，为0.511711～0.512002，初始87Sr/86Sr(ISr)在同一范围内，为0.710658～0.711568。计算的亏损地幔模式年龄(tDM2)分别为2107～2153Ma(似斑状粗粒花岗闪长岩)、2387～2425Ma(花岗伟晶岩脉)、2145～2202Ma(似斑状中粒花岗闪长岩)、2220～2267Ma(细粒花岗闪长岩)和1831～1976Ma(含斑细粒花岗闪长岩)。

北截岩体样品的Pb同位素组成列于表5。初始Pb同位素组成中，斑状粗粒花岗闪长岩、斑状中粒花岗闪长岩和含斑细粒花岗闪长岩的(206Pb/204Pb)i极为相近(17.455～17.945)，略高于伟晶岩脉和细粒花岗闪长岩(17.047～17.262)；前三者的(207Pb/204Pb)i也极为相近(15.526～15.577)，略高于伟晶岩脉和细粒花岗闪长岩(15.490～15.510)；三者的(208Pb/204Pb)i(38.061～38.389)也略高于伟晶岩脉和细粒花岗闪长岩(37.744～37.965)。

5 讨论

5.1 北截岩体的形成时代

北截岩体斑状粗粒花岗闪长岩的锆石LA-ICP-MS U-Pb加权平均年龄为131.53±0.86Ma，16个测点的年龄数据集中于5Ma的较小变化范围内，说明这些锆石是同一次岩浆事件的产物，其加权平均年龄指示了花岗闪长岩的侵位时间。这一年龄值与前人获得的该岩体的形成年龄(128±6Ma)比较接近(Wangetal., 1998)。

花岗伟晶岩脉的锆石LA-ICP-MS U-Pb加权平均年龄为127.4±1.3Ma，其年龄值低于花岗闪长岩年龄，与野外观察的伟晶岩脉侵入花岗闪长岩的事实吻合，说明花岗闪长岩的测年结果是可靠的。花岗伟晶岩10个锆石测点的最大年龄值为130.4Ma，位于花岗闪长岩16个测点年龄范围之间，说明二者形成时间具有连续性。花岗伟晶岩与花岗闪长岩在时间和空间上的紧密伴生，以及其SiO2与其余的氧化物之间呈现的良好线性关系(图5)，指示二者具有成因联系，伟晶岩是由形成花岗闪长岩的岩浆进一步演化形成的。

以往对北截岩体周边的其它郭家岭序列花岗闪长岩进行了较多同位素年龄测试，其中郭家岭岩体花岗闪长岩锆石U-Pb年龄为125.4±2.2Ma(罗贤冬等，2014)，丛家岩体的锆石U-Pb年龄为126.0±0.6和123.0±0.5Ma(Yang and Badal, 2013)，上庄岩体花岗闪长岩的锆石U-Pb年龄为128.8±2.0Ma(罗贤冬等，2014)和129±1Ma(Yang and Badal, 2013)，新城岩体石英二长岩锆石U-Pb年龄为128±1Ma～132±1Ma(Wangetal., 2014)、二长花岗岩的锆石U-Pb年龄为123±1Ma(刘跃等，2014)，三山岛岩体的锆石U-Pb年龄为128±2Ma(Wangetal., 1998)。北截岩体的年龄值与新城石英二长岩、上庄花岗闪长岩等一致，而略高于其它岩体，说明北截岩体是郭家岭型花岗岩早阶段的侵入体。

图9 北截岩体花岗闪长岩黑云母源区投点图 (底图据楼亚儿和杜杨松, 2006)Fig.9 Diagram of the source of biotites for the granodiorite of the Beijie pluton (base map after Lou and Du, 2006)

5.2 地球化学性质与岩石成因

上述主量和微量元素测试结果表明，北截岩体的岩石样品总体显示花岗岩地球化学特征，为高钾钙碱性-钙碱性系列、过铝质花岗岩。除花岗伟晶岩外，绝大部分样品的Na2O/K2O>1，属钠质花岗岩。稀土配分型式为轻稀土富集型，具正铕异常，微量元素的Ba、Sr含量高。岩石的这种地球化学特征与太古宙TTG岩系、埃达克岩和科迪勒拉安第斯等地的Na质花岗岩及高Ba-Sr花岗岩相似，而与典型岛弧环境的花岗岩类明显不同(杨进辉等，2003)。在矿物组成上，岩石中斜长石含量为钾长石含量的2倍以上，暗色矿物由少量角闪石和黑云母组成，按实际矿物含量分类为典型的花岗闪长岩类。具有岩浆演化早期阶段的岩相学特征，与根据同位素年龄测试结果推断的岩石形成顺序吻合。岩相学观察表明，斜长石和钾长石均具有反环带结构，而且斜长石有被熔蚀现象，可能与中-基性岩浆的注入有关(Anderson，1984)，指示北截花岗闪长岩具有岩浆混合成因(Wyllieetal., 1962；Xuetal., 2004)。

北截岩体花岗闪长岩和伟晶岩的主量元素具有良好的线性关系，在Harker图解中(图5)，由含斑细粒花岗闪长岩、斑状粗粒花岗闪长岩、细粒花岗闪长岩至花岗伟晶岩，随SiO2含量增加TiO2、Al2O3、Fe2O3、MgO、CaO、P2O5和MnO呈线性降低，K2O则与SiO2呈正相关关系， 呈现典型的同源岩

图10 北截岩体Sr/Y-Y图解(a，底图据Martin et al., 2005)和(Dy/Yb)N-(La/Yb)N图解(b)Fig.10 Diagrams of Sr/Y vs. Y (a, base map after Martin et al., 2005) and (Dy/Yb)N vs. (La/Yb)N (b) for the granodiorite and pegmatite from the Beijie pluton

图11 北截岩体Pb同位素物质源区判别图解(底图据Doe and Zartman, 1979)Fig.11 Source divisional diagram of Pb isotope for the granodiorite and pegmatite from the Beijie pluton (base map after Doe and Zartman, 1979)

图12 北截岩体的εNd(t)-ISr图解(底图据Defant and Drummond, 1990)数据范围来源：胶东岩群据杨进辉等(2003)；基性脉岩据Yang et al. (2012)；玲珑型花岗岩据杨进辉等(2003)，Zhang et al. (2010)，Yang et al. (2012)，Jiang et al. (2012)和Ma et al. (2013)；郭家岭型花岗岩据杨进辉等(2003)，Yang et al. (2012)，王中亮等(2014)；扬子下地壳、华北上地壳和下地壳据Jahn et al. (1999)Fig.12 εNd(t) vs. ISr diagram for Late-Mesozoic magmaic rocks from the Beijie pluton, Jiaodong Peninsula (base map after Defant and Drummond, 1990)Data range sources: Jiaodong Group (Yang et al., 2003); Linglong granitoids (Yang et al., 2003, 2012; Zhang et al., 2010; Jiang et al., 2012; Ma et al., 2013); Guojialing granitoids (Yang et al., 2003, 2012; Wang et al., 2014); Yangze Craton, upper and lower crusts of the North China Craton (Jahn et al., 1999)

浆演化特征，早期岩浆富Fe、Mg、Ca，晚期岩浆富Si、K。在La/Sm-La和Rb/Nb-Rb/Zr相关图解中(图8)，样品数值投点大致呈斜线排列，显示了部分熔融演化趋势。

花岗闪长岩富SiO2，贫MgO、MnO和CaO，A/CNK接近或大于1，富集大离子亲石元素和轻稀土元素，亏损高场强元素，显示了壳源地球化学特征。在黑云母源区投点图上(图9)，样品投点全部落入壳源区，总体显示壳源特征；但黑云母MgO的含量(9.59%～10.49%)高于壳源的玲珑型花岗岩(平均9.5%)(赛盛勋等，2016)，指示有幔源物质混入的特点。岩石具有高Sr/Y比值(41～479)，在Sr/Y-Y图解(图10a)上，所有样品均落入埃达克岩区域(Defant and Drummond, 1990; Martinetal., 2005)，由斑状粗粒花岗闪长岩、斑状中粒花岗闪长岩、含斑细粒花岗闪长岩、细粒花岗闪长岩至花岗伟晶岩Sr/Y逐渐降低，呈现由埃达克岩向岛弧岩石区演化的趋势。岩石在偏高的全碱含量(K2O+Na2O=7.30%～10.47%)、低的Al2O3含量(13.04%～15.82%)和MgO含量(0.06%～0.77%)以及相对平坦的重稀土配分模式(图6a)方面，则与典型埃达克岩有明显差异。岩石具有高Ba(258×10-6～3330×10-6)、Sr(283×10-6～1371×10-6)和低Rb(39.3×10-6～172×10-6)、Y(1.2×10-6～7.28×10-6)、Yb(0.08×10-6～0.77×10-6)含量特征，与高Ba-Sr花岗岩的地球化学特征相似(Fowleretal., 2001; Yeetal., 2008; Pengetal., 2013)。

岩石富集放射性成因铅((206Pb/204Pb)i=17.047～17.945)，在Pb同位素物源判别图解(图11)上，所有数据均落入下地壳区域附近；样品的Sr/Y比值(41～479)特征，与由加厚下地壳部分熔融形成的花岗岩(Heetal., 2011)相似；Nd同位素模式年龄(tDM2)为2107～2425Ma，与古元古代荆山群和粉子山群的同位素年龄一致(宋明春和王沛成，2003；Lietal., 2012)；另外，岩石具有的高Bb、Sr和低Y特征，与来源于加厚下地壳部分熔融的玲珑花岗岩(Zhangetal., 2010; Jiangetal., 2012)类似。因此认为，由胶北早前寒武纪变质基底岩石组成的加厚下地壳部分熔融是北截岩体的重要物质来源。岩石ISr值(0.710658～0.711568)位于华北下地壳范围；εNd(t)值(-18.53～-11.17)高于壳源的玲珑型花岗岩(-25.1，Yangetal., 2012)，与胶东地区的基性脉岩接近(Yangetal., 2004)，在Sr-Nd相关图上落入以往测试的郭家岭型花岗岩范围(图12)，接近于EM2型富集地幔源区，不同于玲珑花岗岩和胶东岩群样品的投点范围，这表明岩体中有幔源物质混入。在同属郭家岭型花岗岩的丛家和郭家岭岩体中见有岩浆成因的微粒闪长质包体，一般认为微粒闪长质包体是岩浆混合的结果(Gossetal., 2010)。综合分析认为，北截岩体既具有地壳物质来源特征，又有地幔物质来源的信息，是下地壳酸性岩浆与幔源基性岩浆混合作用的结果。

北截岩体较高的Sr、Ba丰度、低的Rb/Sr、中等K/Rb比值，暗示岩浆作用过程中不会有大量长石和云母的分离结晶作用。而高的La/Yb、Sr/Y比值和低Y、Yb含量，以及(Dy/Yb)N与(La/Yb)N呈现正相关关系(图10b)，反映出酸性端元的岩浆源区有大量石榴子石残留(Liuetal., 2010；Heetal., 2011)。岩石具正铕异常，重稀土较强烈亏损，暗示源区由石榴石+辉石组成，没有斜长石残留，残留相为榴辉岩。

北截岩体中伟晶岩的元素地球化学特征表明，其与花岗闪长岩有亲缘关系，二者具有同源岩浆演化特点，伟晶岩由富含挥发份的花岗质母岩浆的残余岩浆结晶分异产生(Fuertes-Fuenteetal., 2000)。伟晶岩的稀土元素总量明显偏低，很可能是岩浆演化过程中副矿物(如磷灰石)结晶分离作用所造成的，个别样品稀土元素含量变化较大则可能与磷灰石等矿物的分异有关(禹丽等，2015)。

5.3 成岩构造环境

元素Rb、Y和Nb能有效区分不同构造环境下形成的花岗岩类，包括火山弧花岗岩、同碰撞花岗岩、板内花岗岩和洋脊花岗岩等。在花岗岩Nb-Y和Ta-Yb构造环境判别图解上(图13)，北截岩体样品全部落入火山弧花岗岩范围内；在Hf-Rb-Ta三角图解上(图14)，所有样品也均为火山弧花岗岩。而研究表明，高Ba-Sr花岗岩形成于张性的或者非挤压的构造背景中，包括岩石圈拉张环境或碰撞造山后与重力垮塌作用有关的构造环境等(Fowleretal., 2001; Yeetal., 2008; Pengetal., 2013)。可见，以北截岩体为代表的郭家岭型花岗岩的构造环境与形成于侏罗纪后碰撞挤压环境的玲珑花岗岩(Yangetal., 2012)完全不同，而与早白垩世中期的伟德山型花岗岩构造环境(Songetal., 2019)相似。玲珑型花岗岩形成于三叠纪华北板块与扬子板块强烈碰撞之后，具后碰撞花岗岩类特征，显示有受华北和扬子克拉通共同影响的印记，物质来源包括造山带俯冲杂岩和加厚的华北克拉通下地壳(Yangetal., 2012)；伟德山型花岗岩具有岛弧花岗岩特征，显示了受太平洋板块俯冲和华北克拉通破坏影响的特征，物质来源于华北太古宙地壳和富集岩石圈地幔熔融(Songetal., 2019)。郭家岭型花岗岩则具有二者之间的过渡性特征，其Sr、Y含量具有与玲珑型花岗岩相似的埃达克岩特征(Yangetal., 2012)，Sr、Nd同位素特征则与伟德山型花岗岩、基性脉岩接近(Yangetal., 2004；Songetal., 2019)，而其构造环境也显示了与伟德山型花岗岩及白垩纪基性脉岩的共同性，可能是太平洋板块俯冲过程中在欧亚板块边缘产生的火山弧的组成部分。

图13 北截岩体Nb-Y(a)和Ta-Yb(b)大地构造环境判别图解(据Pearce et al., 1984)虚线区域为洋脊花岗岩异常区Fig.13 Nb vs. Y (a) and Ta vs. Yb (b) discrimination diagrams for the granodiorite and pegmatite from Beijie pluton (after Pearce et al., 1984)Dashed line arises from the boundary line abnormal ridges ORG

图14 北截岩体Hf-Rb/10-3Ta (a)和Hf-Rb/30-3Ta (b)大地构造环境判别图解(据Harris et al., 1986)Fig.14 Hf-Rb/10-3Ta (a) and Hf-Rb/30-3Ta (b) discrimination diagrams for the granodiorite and pegmatite from Beijie pluton (after Harris et al., 1986)

研究表明，胶东地区在中三叠世受华北克拉通与扬子克拉通强烈碰撞的影响(宋明春等，2020a)，陆壳强烈加厚。侏罗纪(约150Ma)时地幔开始上隆，岩石圈由加厚向减薄转化，造山带根部垮塌，由造山带物质或扬子克拉通和华北克拉通基底物质混合组成的下地壳活化，大范围陆壳重熔，岩浆上侵形成了玲珑型花岗岩(图15a)。白垩纪，华北克拉通东缘经历了克拉通破坏、岩石圈减薄并伴随软流圈物质上涌(Menzieetal., 1993; Gaoetal., 2002, 2009; Xuetal., 2009; Yangetal., 2012)，发生了强烈的构造活动、岩浆作用及盆地裂陷，在胶东地区产生了诸如沂沭裂谷系、胶莱盆地、花岗岩类、火山岩系、中-基性脉岩等广泛的构造岩浆活动，北截岩体即形成于这一构造岩浆背景中，其同位素年龄显示，岩体形成于早白垩世构造岩浆活动阶段的早期，是这一构造阶段的初幕。华北克拉通岩石圈减薄及其伴生的软流圈物质上涌与古太平洋板块的俯冲有关(Gossetal., 2010; Santosh, 2010)，早白垩世早期(约130Ma)，处于太平洋板块俯冲的较早阶段，壳幔作用尚较弱，引发的岩浆活动规模不大，在欧亚板块边缘形成少量火山弧环境的郭家岭型花岗岩，并有幔源的基性脉岩产生(图15b)，郯庐断裂开始形成。早白垩世中期(约120Ma)，俯冲板块的前端后撤，岩石圈地幔拆沉，软流圈由深部向浅部上涌，引起减压熔融，富集的岩石圈地幔部分熔融产生基性岩浆。热的基性岩浆上升到地壳底部发生底侵，引起了地壳底部岩石的部分熔融，产生花岗质岩浆。基性岩浆上侵分异出煌斑岩、高镁闪长岩(宋明春等，2020b)等，幔源和壳源岩浆混合及结晶分异形成胶东地区大范围分布的伟德山型花岗岩、崂山型花岗岩和青山群火山岩等。同时，郯庐断裂强烈活动，并产生“热隆-伸展”构造和大规模金成矿作用(宋明春等，2018)(图15c)。

图15 胶东晚中生代岩浆活动的地球动力学机制(据Yang et al., 2012；Ma et al.,2014；宋明春等，2015修改)Fig.15 Geodynamic mechanism of Late Mesozoic magmatic activity in Jiaodong (modified after Yang et al., 2012；Ma et al., 2014；Song et al., 2015)

5.4 岩浆活动与金成矿的关系

关于胶东金矿与中生代岩浆活动的关系及成矿动力学背景和成矿模式前人已进行了较多研究(Yangetal., 2016b； Songetal., 2017, 2019；宋明春等, 2018, 2019；Dengetal., 2020a, b；Grovesetal., 2020；Zhangetal., 2020)。鉴于包括北截岩体在内的郭家岭型花岗岩内部和周边分布有很多金矿床，许多地质工作者认为郭家岭型花岗岩岩浆活动与金成矿密切相关，将金矿床归类为郭家岭型花岗岩岩浆期后热液矿床(Wangetal., 1998；罗贤冬等，2014)。然而，本文测试的北截岩体的形成年龄，明显早于金矿成矿时代，即使岩浆活动末期的花岗伟晶岩的同位素年龄(127.4±1.3Ma)也高于大部分研究者公认的金矿成矿年龄(120Ma左右)(翟明国等，2001；宋明春等，2014)。总体分析，胶东金矿的成矿时间晚于郭家岭型花岗岩5Myr以上(翟明国等，2001)。前人对胶东西北部金矿床成矿深度的估算结果表明，成矿深度上限为2～3km，下限为5～8km(郭春影等，2011)；测算的郭家岭岩体侵位深度约为13±1.6km(张华锋等，2006)。可见，金矿的成矿深度浅于郭家岭花岗岩侵位深度5km以上(宋英昕等，2017)。综合考虑成矿时代与岩浆活动的时差及岩浆期后热液可能存在的时间范围(<0.8Myr)，大量流体参与金成矿的事实(沈昆等，2000；Fanetal.，2007；Songetal.，2014；杨立强等，2014；Deng and Wang, 2016；Dengetal.，2020a，b)，本文和部分研究者(姜晓辉等，2011；宋明春等，2015；Songetal.，2019)认为金成矿流体可能不是来自于这种花岗岩岩浆期后热液。

灵北断裂是胶西北焦家金断裂成矿带下盘的一条重要控矿断裂，断裂全长约50km，走向 40°～60°，断裂西南段穿切玲珑型花岗岩，其北东段至北截一带切割北截岩体(图2)，沿该断裂分布有北截、蚕庄、灵山沟、黄埠岭等中小型金矿床和多个金矿点。在北截地区，灵北断裂明显切割郭家岭型花岗闪长岩，断裂宽4～30m，构造岩主要为花岗质碎裂岩和角砾岩，绢云母化、硅化和团块状黄铁矿化发育，北截金矿赋存于该断裂破碎带中。按照地质体的相互关系，确定北截地区地质事件由早至晚的形成顺序为：玲珑型花岗岩→郭家岭型花岗岩(北截岩体)→灵北断裂→金矿化，北截金矿显然是在郭家岭型花岗岩形成固结并被灵北脆性断裂破坏后再形成的。这种花岗岩显然已不能为金成矿提供流体或热动力条件，但由于岩体中断裂裂隙发育而为金成矿提供了赋矿空间，其在金成矿中的作用同赋矿的玲珑型花岗岩和早前寒武纪基底变质岩系一样。

野外地质调查显示，郭家岭型花岗岩被伟德山型花岗岩侵入，鉴于伟德山型花岗岩的艾山、海阳、牙山、三佛山、伟德山等岩体侵位深度普遍小于3.5km(张华锋等, 2006)，说明伟德山型花岗岩侵位时，郭家岭型花岗闪长岩发生了强烈隆升。大量同位素年龄测试表明，伟德山型花岗岩侵位的峰值年龄为114Ma(Songetal., 2019)，从郭家岭型花岗岩侵位至伟德山型花岗岩形成，在约10Myr内，郭家岭型花岗岩岩体的隆升量达10km左右(豆敬兆等, 2015)。可见，胶东地区早白垩世早中期，发生了大规模岩浆活动和强烈的地壳隆升事件，早白垩世也是华北克拉通破坏的峰期和伸展构造大量发育期，岩浆活动、地壳隆升和伸展构造共同构成了“热隆-伸展”构造系统(宋明春等, 2018)。包括北截岩体在内的郭家岭型花岗岩于早白垩世早期侵位后，快速隆升降温，在区域伸展构造作用下产生一系列断裂裂隙，为金成矿提供了有利空间。随着早白垩世中期伟德山型花岗岩的强烈活动，出现了较广泛分布的含金流体，由于地壳快速隆升，流体显著降温，在郭家岭型花岗岩、玲珑型花岗岩和前寒武变质岩中的先成断裂构造中聚集的丰富含金流体，在适当的温压条件下流体沸腾、金质沉淀成矿。

6 结论

(1)北截岩体斑状粗粒花岗闪长岩和花岗伟晶岩脉的锆石U-Pb年龄分别为131.53±0.86Ma和127.4±1.3Ma，二者为同一期岩浆活动的产物，为郭家岭型花岗岩早阶段的岩浆侵入体。

(2)北截岩体岩石具有高的全碱、Ba、Sr含量，低的Al2O3、MgO、Rb、Y、Yb含量特征，及相对平坦的重稀土配分模式，与高Ba-Sr花岗岩的地球化学特征相似。

(3)北截岩体的ISr值为0.710658～0.711568，εNd(t)为-18.53～-11.17，二阶段模式年龄为2107～2425Ma；Pb同位素数据均落于下地壳区域。结合岩石地球化学认为，北截岩体既具有地壳物质来源特征，又有地幔物质来源的信息，是下地壳酸性岩浆与幔源基性岩浆混合作用的结果，郭家岭型花岗岩是太平洋板块向欧亚板块俯冲产生的欧亚板块边缘火山弧的组成部分。

(4)郭家岭型花岗岩为胶东大规模金成矿提供了有利的赋存空间，而不是成矿流来源或热动力条件。

致谢山东省地质调查院刘汉栋研究员在野外工作期间提供了大力支持；论文写作过程中得到了中国地质大学(北京)邓军教授团队的热情帮助；两位审稿人为论文提出宝贵意见；在此一并表示最诚挚的感谢！