胶东地区东部晚侏罗世花岗岩锆石U-Pb定年、Hf同位素特征及其对金成矿构造背景的限定 *

王栋 李红艳 王天齐,3 于洪 宿建恒

1. 中国地质大学(北京)地球科学与资源学院，北京 100083 2. 中国地质科学院矿产资源研究所，自然资源部成矿作用与资源评价重点实验室，北京 100037 3. 中国科学院地质与地球物理研究所，岩石圈演化国家重点实验室，北京 100029

胶东地区已累积探明黄金储量接近5000t(宋明春等，2019)，是我国目前最大的金矿产区，也是环太平洋成矿系统的重要组成部分之一(Goldfarbetal., 2014; 范宏瑞等，2021)。该地区在中生代被卷入了复杂的构造演化过程，包括以下4期构造事件：(1)扬子板块与华北板块在三叠纪时期发生陆-陆碰撞及后续的大陆深俯冲作用，形成苏鲁超高压变质带(Lietal., 1993; Yeetal., 2000; Zheng, 2008)；(2)郯庐断裂带自三叠纪以来发生了多次左旋走滑运动，引发构造-岩浆(流体)活动及成矿过程(邓军等，2000；朱光等，2009)；(3)晚三叠世发生岩石圈的减薄(徐义刚等，2009；杨进辉和吴福元，2009；朱日祥等，2011；Wuetal., 2019)；(4)至少在侏罗纪时期，中国东部就开始受到太平洋板块俯冲、回撤作用的影响(孙卫东等，2008; Jiangetal., 2010)，胶东地区一系列白垩纪岩浆活动都与此有着密切的关系(Dengetal., 2017)。这些重大地质事件在胶东地区先后叠加，导致该区深部构造体制在中生代发生转换(邓军等，2006；宋明春，2014)、岩石圈地幔遭受破坏和改造(徐义刚等，2009；朱日祥等，2011；Wuetal., 2019)，而同时期强烈的岩浆活动则是对整个构造演化的响应(霍腾飞等，2016)。

复杂的地质演化过程给讨论该地区金矿的物质来源及成矿过程提出了挑战，目前大多数学者认为该地区大规模成矿作用主要发生在110～130Ma(Yang and Zhou, 2001; 陈衍景等，2004；Dengetal., 2015; 朱日祥等，2015；范宏瑞等，2021)。然而，在胶东地区广泛发育的晚侏罗世花岗岩却与金矿有着密切的空间关系(Guoetal., 2013; Deng and Wang, 2016)，是该区富产金矿的重要地质背景。例如，胶北玲珑型花岗岩是该区最主要的赋矿围岩之一(Wenetal., 2015; Yangetal., 2016a, b; Zhangetal., 2020)。前人对这些花岗岩的年代学、矿物学、地球化学等方面进行了大量的研究，对该期花岗岩的特征及成因进行了广泛深入的讨论(Wangetal., 1998; 苗来成等，1998; 郭敬辉等，2005; Houetal., 2007; Zhangetal., 2010; Jiangetal., 2012; Yangetal., 2012, 2014, 2018; Maetal., 2013; 黄涛等，2014; 刘跃等，2014; Zhaoetal., 2016; Dengetal., 2018; Heetal., 2021)。已有的研究显示，这些岩体可能形成于相同的构造体制。胶东东部晚侏罗世花岗岩体中赋存的金矿数量明显少于胶北地体，因此前人对于成矿作用的研究主要集中在胶北金矿区，而对胶东东部出露的同期花岗岩体的研究相对较少。显然，开展胶东东部及西北部同一时期岩体的差异性研究，可对研究该区的金成矿作用提供重要的助力。

本文在前人研究的基础上，对胶东地区东部出露的两个晚侏罗世花岗岩体进行了锆石U-Pb和Hf同位素分析，并对比了出露在胶东几个不同构造单元内的同时代花岗岩体的继承锆石和Hf同位素特征，探讨了该期花岗岩的源区及成因，为研究胶东地区在侏罗纪时期的岩浆活动、构造演化及其对金成矿的贡献提供一定的制约。

1 区域地质背景及样品采集

胶东地区位于华北克拉通的东南缘，西侧与鲁西地块界于郯庐断裂，南侧为苏鲁超高压变质带(图1)。根据岩石和构造特征，胶东地区可分为三个构造单元：苏鲁超高压变质地体、胶北地体和胶莱盆地(张增奇等，2014；潘素珍等，2015)。

图1 胶东地区地质简图(据李洪奎等，2011; Deng et al., 2015修改)Fig.1 Geological sketch map of Jiaodong Peninsula (modified after Li et al., 2011; Deng et al., 2015)

苏鲁超高压变质地体位于扬子板块北缘，形成于三叠纪华北板块与扬子板块陆-陆碰撞及后续的深俯冲作用。地体内广泛发育超高压变质岩，其原岩年龄在780～740Ma(Liuetal., 2006; 薛怀民等，2006; Zongetal., 2010; Chenetal., 2013; 高名迪等，2018)。胶北地体位于华北板块东南缘，基底岩石主要包括太古代TTG片麻岩和胶东岩群、古元古代荆山群和粉子山群以及新元古代蓬莱群(Tangetal., 2007; Tametal., 2011; Shan and Zhai, 2020)。TTG片麻岩和胶东岩群的原岩年龄主要为～2.9Ga、～2.7Ga、～2.5Ga(Jahnetal., 2008; Liuetal., 2013)，并普遍遭受～1.8Ga的角闪岩相区域变质作用(Tangetal., 2007)。荆山群和粉子山群，主要岩性为黑云斜长片麻岩、石榴夕线黑云片岩、黑云变粒岩、斜长角闪岩和大理岩等，直接覆盖在TTG片麻岩和胶东岩群之上(Wangetal., 1998)。蓬莱群为一套浅变质沉积岩，岩性主要为大理岩、板岩、石英岩、千枚岩和泥灰岩等，不整合覆盖在粉子山群之上(Lietal., 2007)。胶莱盆地是一个形成于中生代期间的伸展断陷盆地(张岳桥等，2008)，盆地基底主要为新太古代胶东岩群、古元古代荆山群和粉子山群等变质岩系，盖层则主要由白垩系莱阳群、青山群及王氏群等陆相沉积岩及部分火山岩、火山碎屑岩等构成(邹为雷等，2010)。

中生代花岗岩在胶东地区发育广泛，根据锆石U-Pb年代学研究可将其分为三期：晚三叠世225～205Ma(Chenetal., 2003; 郭敬辉等，2005；曾令森等，2007)、晚侏罗世160～150Ma(苗来成等，1998；郭敬辉等，2005；Zhangetal., 2010; Jiangetal., 2012; Maetal., 2013; Yangetal., 2014; Heetal., 2021)和早白垩世130～110Ma(Houetal., 2007; Gossetal., 2010; Lietal., 2012; Yangetal., 2012; 李增达等，2018; Gaoetal., 2019; 董学等，2020; 宋英昕等，2020)。晚三叠世花岗岩主要出露在苏鲁超高压变质地体中，岩性以碱性的正长花岗岩为主，是俯冲的扬子陆壳板片断离引发的岩石圈地幔部分融熔的结果(郭敬辉等，2005；曾令森等，2007；陈竟志和姜能，2011)。晚侏罗世花岗岩在该区最为发育，包括胶北地体的玲珑和栾家河岩体，胶东东部的昆嵛山、垛崮山和文登岩体。该期花岗岩主要包括偏铝质到过铝质黑云母花岗岩、二长花岗岩和花岗闪长岩，且含有中生代印支期变质锆石(200～250Ma)及少量前寒武纪继承锆石，属于地壳重熔花岗岩(郭敬辉等，2005; Zhangetal., 2010; Jiangetal., 2012; Yangetal., 2012, 2014)。早白垩世花岗岩主要包括郭家岭型似斑状花岗闪长岩-二长花岗岩，伟德山型(或称艾山型)似斑状二长花岗岩，崂山型带晶洞碱长花岗岩-二长花岗岩。岩体局部可见暗色包体或团块，地球化学特征显示具有壳幔岩浆混合成因的特点(Houetal., 2007; Gossetal., 2010; Yangetal., 2012; Yan and Shi, 2014; 李增达等，2018; Gaoetal., 2019)。

本次研究的样品为出露在胶东东部的晚侏罗世花岗岩，包括位于苏鲁超高压地体内的文登岩体和位于苏鲁超高压地体与胶北地体交界处的垛崮山岩体(图1b)，所有样品均采自新鲜露头，未经历明显风化。其中文登岩体包括何家店花岗岩和架子山花岗岩。何家店花岗岩样品(SD06-10-01)为粗粒二长花岗岩，取自文登市东部的何家店采石场(N37°12′6″、E122°8′55″)。架子山花岗岩包括含斑中粗粒二长花岗岩(SD06-11-01)和正长花岗岩(SD06-11-05)，取自文登市东南部的架子山(N37°9′5″、E122°6′17″)。垛崮山岩体包括弱片麻状花岗闪长岩(SD06-19-01)和呈脉状产出的正长花岗岩(SD06-19-04)，二者均取自乳山市东南部的窗龙山(N36°53′14″、E121°44′45″)。

2 分析方法

在野外观察和室内岩相学研究的基础上，选取新鲜的样品进行了锆石挑选和制靶，并在中国地质科学院离子探针中心进行锆石阴极发光(CL)图像的采集。锆石U-Pb定年是在中国科学院地质与地球物理研究所岩石圈演化国家重点实验室的激光烧蚀等离子体质谱仪(LA-ICPMS)上完成。激光剥蚀系统采用GeoLasPlus型193nm ArF准分子激光器，ICPMS为Agilent7500a型四极杆等离子体质谱仪。激光束斑直径为32μm，年龄计算时以标准锆石91500和GJ-1为外标进行同位素比值校正，元素含量以国际标样NIST610为外标，Si为内标计算。测试数据采用204Pb进行普通铅校正，单个测试点的同位素比值误差均为1σ。详细的仪器运行条件和分析流程见谢烈文等(2008)。数据处理采用Glitter软件利用标准锆石91500线性内插的方式对分析时间相关的U-Th-Pb进行同位素比值校正。最后锆石的U-Pb谐和图和年龄加权平均计算均采用Isoplot软件(Ludwig, 2003)完成。

锆石Hf同位素测定是在中国科学院地质与地球物理研究所岩石圈演化国家重点实验室的激光剥蚀多接收电感耦合等离子体质谱仪(LA-MC-ICPMS)上完成。激光剥蚀系统为GeolasPlus型193nm ArF准分子激光器，多接收等离子体质谱仪为Neptune Plus。实验采用的束斑大小为40～50μm，直接覆盖在锆石U-Pb定年的测定部位，本次测试采用锆石标样Mud-Tank和FM0411为分析标样，仪器运行条件及详细分析过程见Wuetal. (2006a)。

3 分析结果

3.1 锆石U-Pb年龄

3.1.1 文登岩体

何家店二长花岗岩样品(SD06-10-01)中的锆石大部分为短柱状晶型，粒径50～150μm。阴极发光图像(图2a)显示大部分锆石颗粒存在振荡环带，显示岩浆成因特点。部分锆石具有核-边结构，核部在形态和结构上与边部的新生岩浆锆石存在明显区别，为继承锆石。还有部分锆石的边部阴极发光为亮白色，无环带结构，显示出变质增生成因的特点。对样品中26颗锆石进行了31个点的测定，测定结果见表1和图3a。其中16个新生岩浆锆石分析点的U含量和Th/U比值分别为72×10-6～1413×10-6和0.26～1.13，谐和年龄范围在155～162Ma，加权平均为158±1Ma。5个继承锆石的U含量和Th/U 比值分别为54×10-6～721×10-6和0.36～1.44，谐和年龄在770～789Ma，加权平均为782±9Ma。2个变质增生锆石的Th/U比值为0.14，谐和年龄为225±5Ma和234±4Ma。剩余的锆石测试点年龄不谐和，但均位于780～160Ma的U-Pb不一致线附近，为新元古代继承锆石在晚侏罗世岩浆热事件中发生Pb丢失造成。

架子山二长花岗岩(SD06-11-01)样品中锆石大部分为短柱状晶型，粒径100～200μm。阴极发光图像(图2b)显示大部分锆石颗粒具有振荡环带或板状环带，具有岩浆成因的特点。与何家店花岗岩中的锆石类似，一些锆石颗粒的核部也含有继承锆石。对样品中25颗锆石进行了30个点的测定，结果见表1和图3b、c。其中20个新生岩浆锆石分析点的U含量和Th/U 比值分别为54×10-6～745×10-6和0.19～1.60，谐和年龄范围在147～167Ma，加权平均为159±2Ma。继承锆石的U含量和Th/U 比值分别为68×10-6～672×10-6和0.34～2.14。其中5个继承锆石给出了773±8Ma的加权平均年龄，另有一个继承锆石谐和年龄为2535±37Ma。剩余锆石测试点的年龄不谐和，位于780～160Ma的U-Pb不一致线附近，为新元古代继承锆石在晚侏罗世岩浆热事件中发生Pb丢失造成。

架子山正长花岗岩(SD06-11-05)样品中锆石大部分为短柱状-长柱状晶型，粒径100～200μm。阴极发光图像(图2c)显示大部分锆石颗粒具有振荡环带或板状环带，属于岩浆成因。部分新生岩浆锆石的内部也含有继承锆石核。对样品中的18颗锆石进行了24个点的测定，测定结果见表1和图3d。其中17个新生岩浆锆石分析点的U含量和Th/U比值为53×10-6～2091×10-6和0.19～1.13，谐和年龄范围在154～165Ma，加权平均为159±1Ma。3个继承锆石的U含量和Th/U比值为147×10-6～402×10-6和0.91～3.09，谐和年龄范围在769～784Ma，加权平均为776±11Ma。剩余锆石分析点的年龄不谐和，位于780～160Ma的U-Pb不一致线附近，是继承锆石在晚侏罗世岩浆热事件中发生Pb丢失造成。

续表1Continued Table 1Spot No.含量(×10-6)Pb∗ThUTh/U同位素比值年龄(Ma)207Pb206Pb1σ207Pb235U1σ206Pb238U1σ207Pb206Pb1σ207Pb235U1σ206Pb238U1σ1280.14864541.070.06490.00101.14130.01690.12750.00147713277387748137.42132220.960.05050.00410.17110.01320.02460.00072181771601115741413.51384860.280.05250.00220.17900.00710.02470.00043099216761573152.352740.700.05400.00590.18010.01880.02420.0008370229168161545164.41291330.970.05020.00450.17230.01490.02490.0007205197161131595173.9261370.190.05370.00620.18900.02090.02550.00103592431761816361818.77644771.600.05470.00250.19240.00840.02550.00044019917971623194.3121730.070.04910.00410.15840.01270.02340.0006150186149111494202.493721.300.05460.00530.17310.01600.02300.00073952041621414742130.01891651.140.06410.00201.14430.03390.12960.0019743647751678511225.6741930.380.05390.00300.18290.00960.02460.0005367119171815732321.52747450.370.04800.00190.16720.00630.02520.000410191157516122428.41211780.680.06680.00161.15350.02620.12530.001683049779127619251.641540.760.05300.00780.16840.02390.02300.0010331305158211476263.9941240.760.04960.00380.17050.01240.02490.00061781671601115942816.4772120.370.06020.00190.50940.01560.06140.00096096841810384529149.214376722.140.06470.00101.13700.01630.12740.00147663177187738302.045590.760.04930.00720.17320.02440.02550.00111623101622116273119.22376270.380.05590.00200.20150.00700.02620.000444679186616723233.31272610.490.06330.00140.92110.01890.10560.001371845663106477架子山正长花岗岩 SD06-11-0514.01321261.050.05160.00470.17190.01510.02420.0007269196161131544216.82395990.400.05070.00160.17710.00530.02540.0003226711665161231.742530.790.05070.00580.18040.02000.02580.000822824516817164552.474800.930.04940.00490.16960.01640.02490.000716521815914159487.31252480.500.05230.00310.18110.01020.02510.0005297129169916031011.1591650.360.05770.00340.32470.01840.04080.000952012628614258611106.78824120.040.05810.00090.32170.00470.04020.000453234283425431216.21186110.190.04930.00250.17270.00840.02540.0005161115162716231399.212424023.090.06730.00141.17570.02270.12670.0015848427891176991426.41691591.060.06600.00171.17030.02780.12860.00178065178713780101518.51201550.770.06280.00200.75550.02300.08730.0013700675711354081628.624510810.230.04990.00120.17360.00410.02520.000319057163416121713.23414100.830.05160.00270.17160.00850.02420.000526511416171543187.1942490.380.04930.00370.17250.01250.02540.00061611681621116241928.629410460.280.05340.00160.18440.00520.02510.000334566172416022023.01341470.910.06590.00221.17490.03790.12920.0020804697891878412224.31291291.000.05040.00480.17580.01610.02530.0007215207164141615233.2671080.620.04960.00670.17730.02290.02600.00111742871662016572417.8931870.500.05690.00330.33670.01860.04290.00094871242951427162624.93079100.340.04910.00150.16740.00500.02470.000315472157415722760.992520910.440.04860.00080.16860.00280.02520.00031294015821602293.91331181.130.04940.00450.17110.01500.02510.0007165202160131605

3.1.2 垛崮山岩体

垛崮山花岗闪长岩(SD06-19-01)中的锆石大部分为短柱状-长柱状，粒径150～250μm。阴极发光图像(图2d)显示大部分锆石具有核-边结构，少数锆石还可识别出核-幔-边结构。锆石的边部具有较宽的振荡环带，显示岩浆成因特点。核部具有熔蚀结构，部分核部仍保留原始的环带结构。对样品中的34颗锆石进行了50个点的测定，测定结果见表1和图3e-g。其中23个新生岩浆锆石分析点的U含量和Th/U 比值为77×10-6～3486×10-6和0.03～1.22，谐和年龄范围在160～169Ma，加权平均为164±2Ma。继承锆石的U含量和Th/U 比值为56×10-6～3170×10-6和0.08～1.58。其中5颗继承锆石给出743～784Ma的谐和年龄，3颗继承锆石给出1669～1820Ma的谐和年龄，2个继承锆石测点给出了两个最老的谐和年龄分别为2358±20Ma和2576±39Ma。6个变质成因锆石测点的Th/U比值为0.08～0.26，其谐和年龄范围在224～243Ma，加权平均为230±7Ma。剩余锆石分析点的年龄不谐和，均位于谐和线右侧，应该是在中生代变质事件或岩浆事件时发生Pb丢失造成的。

垛崮山正长花岗岩(SD06-19-04)中的锆石大部分为长柱状晶型，粒径150～300μm。阴极发光显示它们具有振荡环带或板状环带(图2e)，为岩浆成因。对15颗岩浆锆石进行了17个点的测定，测定结果见表1和图3h。17个新生岩浆锆石分析点的U含量和Th/U 比值分别为254×10-6～3048×10-6和0.01～1.01，谐和年龄范围在150～158Ma， 加权平均为153±1Ma。

3.2 锆石Hf同位素

Hf同位素分析数据见表2及图4。

3.2.1 文登岩体

何家店二长花岗岩(SD06-10-01)的新生岩浆锆石εHf(t)=-27.8～-24.1 (平均值-25.6，n=16)，tDM2=2728～2961Ma。13颗继承锆石εHf(t)=-10.9～-4.7，tDM2=1978～2367Ma。2颗变质锆石εHf(t)=-24.9～-21.7，tDM2=2630～2837Ma。

架子山二长花花岗岩(SD06-11-01)的新生岩浆锆石εHf(t)=-30.7～-20.3(平均值-25.9，n=20)，tDM2=2494～3145Ma。7颗继承锆石εHf(t)=-11.4～-0.7，tDM2=1626～2298Ma。另外1颗年龄为～2535Ma的继承锆石核测试点εHf(t)=5.7，tDM2=2675Ma。

架子山正长花岗岩(SD06-11-05)的新生岩浆锆石εHf(t)=-29.4～-24.1 (平均值-26.7，n=17)，tDM2=2730～3059Ma。7个继承锆石测点εHf(t)=-12.5～1.0，tDM2=1616～2467Ma。

3.2.2 垛崮山岩体

垛崮山花岗闪长岩(SD06-19-01)的新生岩浆锆石εHf(t)=-27.5～-17.9 (平均值-22.1，n=23)，tDM2=2344～2947Ma。18颗继承锆石εHf(t)=-16.1～11.3，tDM2=1654～3021Ma。6颗变质锆石εHf(t)=-21.9～-9.2，tDM2=1845～2653Ma。

垛崮山正长花岗岩(SD06-19-04)的新生岩浆锆石εHf(t)=-25.9～-22.2 (平均值-24.2，n=17)，tDM2=2606～2838Ma。

4 讨论

4.1 文登岩体和垛崮山岩体的形成时代

本次研究获得文登岩体3个样品及垛崮山岩体花岗闪长岩样品的岩浆锆石年龄分别为158±1Ma、159±2Ma、159±1Ma和164±2Ma，表明文登岩体和垛崮山岩体均形成于晚侏罗世(～160Ma)。侵入到垛崮山花岗闪长岩体的正长花岗岩脉体(SD06-19-04)获得了较年轻的年龄(153±1Ma)，与北部的昆嵛山岩体定年结果一致(郭敬辉等，2005；Zhaoetal., 2016; Heetal., 2021)，代表稍晚一期花岗岩浆活动。郭敬辉等(2005)通过SHRIMP锆石U-Pb定年获得文登岩体和垛崮山岩体的年龄分别为160±3Ma和161±1Ma。Zhaoetal.(2016)对文登东部花岗岩的锆石U-Pb定年获得的年龄较为年轻，为151±1Ma。本文的结果与郭敬辉等(2005)及 Zhaoetal.(2016)的结果在误差范围内一致，表明胶东地区东部在晚侏罗世也存在一期较强的岩浆作用。

4.2 晚侏罗世花岗岩的源区――继承锆石的U-Pb年龄证据

花岗岩通常产生于深部地壳物质的重熔，因此岩石中常残留有难熔的壳源物质，尤其是锆石。因为锆石具有极高的矿物学稳定性，即使在经历多期热事件后，仍能够以继承锆石的形式存在于新生的岩石中，同时能够保留寄主岩石的一些信息(黄涛等，2014)。本次研究的文登岩体和垛崮山岩体中的继承锆石主要集中在新元古代以及三叠纪两个时期，此外还有部分继承锆石具有1700～1800Ma、～2400Ma、～2500Ma热事件的年龄记录。这些继承锆石的年龄数据，不仅能为两个岩体的源区提供限定，还可以为华北、扬子两个陆块的碰撞及演化提供重要的依据。

华北板块和扬子板块在中生代碰撞汇聚以前，各自经历了不同的岩浆-热事件(Bruguieretal., 1997; Hackeretal., 1998, 2000; Chenetal., 2003; Lietal., 2005a; Zhengetal., 2005, 2006; Tangetal., 2007)。新元古代是扬子板块所在的华南板块发生生长和再造的重要时代，并且扬子板块北缘广泛发育新元古代(750～820Ma)的岩浆活动(Zhouetal., 2002, 2006; Lietal., 2003a, b; Zhengetal., 2004; Wangetal., 2006; Wuetal., 2006b)。而在华北板块南缘则没有新元古代岩浆活动的记录，因此新元古代的岩浆岩年龄也是扬子板块区别于华北板块的一个重要特征(Hackeretal., 1998; 郑永飞，2003; Zhengetal.，2004; Tangetal., 2007, 2008)。文登岩体和垛崮山岩体均含有较多新元古代(700～800Ma)的继承锆石，这表明两个岩体的岩浆源区主要来自扬子板块。除新元古代年龄外，这两个岩体还含有许多三叠纪变质成因锆石，这与大别-苏鲁超高压变质带的变质峰期时代(～230Ma)(Liuetal., 2004, 2006; Wuetal., 2006c)相吻合，而且大别-苏鲁超高压变质岩的原岩年龄也是新元古代(740～780Ma)(Zhengetal., 2004; Tangetal., 2008; Liuetal., 2010)。样品中有部分核-边结构的锆石，它们核部的继承锆石为新元古代(～780Ma)的谐和年龄，边部变质锆石为三叠纪(～230Ma)的谐和年龄(图2a)，这与大别-苏鲁超高压变质岩中发现的核-边结构的锆石特征(Zheng, 2008)一致。这些锆石年代学特征表明，文登岩体和垛崮山岩体的主要源区物质是经历了超高压变质作用改造的扬子陆块。

华北板块前寒武纪基底岩石主要包括太古代花岗质片麻岩以及古元古代变质火山-沉积岩(Zhaietal., 2005; Zhaoetal., 2005)。花岗质片麻岩的原岩年龄范围在2900～2500Ma，峰值范围是2900～2700Ma和2600～2500Ma(Zhai and Liu, 2003; Zhaoetal., 2005)，并主要经历了1900～1800Ma的区域角闪岩到麻粒岩相的变质作用(Zhai and Liu, 2003; Zhaietal., 2005; Santoshetal., 2007)。李江海等(2006)统计了华北板块近500件高精度的锆石U-Pb年龄数据，结果显示华北板块基底经历了多期次的构造热事件，大陆生长演化过程中重大的热事件幕为2700～2600Ma、2520Ma、2400Ma、2400～2250Ma、1950Ma、1840Ma，其中2520Ma、1840Ma分别对应于华北太古宙/元古宙和古元古代/中元古代地层界线，是最为显著的构造热事件。包括全岩Sm-Nd、锆石U-Pb等多种同位素年代学研究(Kröneretal., 1998; 金文山和管爱莲，1999; 李江海等，2000; 杨进辉等，2005)也证明～2500Ma和～1800Ma是华北板块地壳增生的两个主要时代，而且2400～2100Ma范围的热事件年龄，也仅在华北板块出现(路孝平等, 2004; Lietal., 2005b; Wanetal., 2006)。胶北地体的花岗质片麻岩同样显示有～2500Ma的热事件年龄(Tangetal., 2007; Jahnetal., 2008; Liuetal., 2013; 万渝生等，2017)。本次研究的文登岩体含有～2500Ma的继承锆石，垛崮山岩体含有～2500Ma和～2400Ma以及～1700Ma的继承锆石，表明文登和垛崮山两个岩体的岩浆源区有部分华北板块的物源贡献。

4.3 晚侏罗世花岗岩的源区――锆石Hf同位素证据

图5 胶东地区晚侏罗世花岗岩体中继承锆石年龄分布直方图数据来源：文登岩体和垛崮山岩体来自本文和Zhao et al., 2016; 玲珑岩体引自Jiang et al., 2012; Yang et al., 2012, 2014, 2018; Ma et al., 2013; 黄涛等，2014; Zhao et al., 2016Fig.5 Age distribution histograms of inherited zircons from Late Jurassic granites in Jiaodong areaData sources: the Wendong and Duogushan pluton from this study and Zhao et al., 2016; the Linglong pluton from Jiang et al., 2012; Yang et al., 2012, 2014, 2018; Ma et al., 2013; Huang et al., 2014; Zhao et al., 2016

图6 晚侏罗世岩体中新生岩浆锆石的εHf(t)频率直方图数据来源：文登岩体和垛崮山岩体数据来自本文；玲珑岩体数据引自Jiang et al., 2012; Yang et al., 2012, 2014, 2018; Ma et al., 2013; Zhao et al., 2016Fig.6 εHf(t) histogram of magmatic zircon from Late Jurassic granite plutonsData sources: Wendeng and Duogushan plutons from this study；Linglong pluton from Jiang et al., 2012; Yang et al., 2012, 2014, 2018; Ma et al., 2013; Zhao et al., 2016

图7 新生岩浆锆石的176Lu/177Hf-176Hf/177Hf比值关系图数据来源：文登岩体和垛崮山岩体数据来自本文;玲珑岩体数据引自Jiang et al., 2012; Yang et al., 2012, 2014, 2018; Ma et al., 2013; Zhao et al., 2016; 郭家岭岩体数据引自Yang et al., 2012, 2014; Wang et al., 2014Fig.7 176Lu/177Hf vs. 176Hf/177Hf diagram of newly generated magmatic zirconsData sources: Wendeng and Duogushan plutons from this study; Linglong pluton from Jiang et al., 2012; Yang et al., 2012, 2014, 2018; Ma et al., 2013; Zhao et al., 2016；Guojialing pluton from Yang et al., 2012, 2014; Wang et al., 2014

文登岩体的三个花岗岩样品中具有非常一致的锆石Hf同位素特征：εHf(t)值和Hf模式年龄均具有双峰式分布，3个样品的εHf(t)值的两个峰值位于-25和-8，Hf两阶段模式年龄的两个峰值范围在2200～2400Ma和2700～2900Ma(图4)，表明三者来源于相同岩浆源区。垛崮山岩体的两阶段Hf模式年龄与文登岩体有所差异，Hf模式年龄变化范围较大，从1800～3000Ma，主要峰期在～2600Ma。其εHf(t)值也呈双峰式分布，峰值位于-21和-8，与文登岩体基本一致，只是变化范围更大。从Hf同位素特征来看垛崮山岩体与文登岩体的主要源区具有相似性，这与它们都具有新元古代和三叠纪的继承锆石年龄相吻合。但垛崮山岩体的εHf(t)值及Hf模式年龄较文登花岗岩体变化更大，这反应了二者的源区并不完全相同，应该还有其他源区的贡献，这也与我们在垛崮山获得了更多的华北板块继承锆石年龄结论相符合。值得注意的，本次文登岩体样品(SD06-11-01)中含有谐和年龄～2500Ma的继承锆石，其εHf(t)值为5.7，一阶段模式年龄为2623±32Ma与U-Pb年龄基本一致，表明～2500Ma的年龄记录了壳幔分异的年龄。同样地，垛崮山岩体样品(SD06-19-01)也含有2576±39Ma和2358±20Ma的继承锆石，其εHf(t)值分别为10.0和11.3，一阶段模式年龄为～2496Ma和～2259Ma，均与U-Pb年龄基本一致。这些壳幔分异年龄与华北地壳增生的时代记录一致(路孝平等，2004；Lietal., 2005b；Wanetal., 2006；Tangetal., 2007；Jahnetal., 2008；Liuetal., 2013；万渝生等，2017)，也进一步证明了在两个岩体的源区中有华北板块的物源贡献。

4.4 晚侏罗世花岗岩继承锆石及Hf同位素对胶东地区深部构造特征的制约

为了更好地理解胶东地区晚侏罗世地质构造特征，本文统计了前人(Jiangetal., 2012; Yangetal., 2012, 2014, 2018; Maetal., 2013; 黄涛等, 2014; Zhaoetal., 2016)发表的胶北地体中晚侏罗世玲珑岩体的继承锆石数据(图5)，进行对比研究。玲珑岩体中继承锆石的年龄包括三叠纪(～230Ma)、新元古代(～780Ma)、1750～2200Ma、～2500Ma、2700～2900Ma，其中最主要的年龄峰值在三叠纪和新元古代，表明玲珑岩体主要源区物质也来自于扬子板块。与文登岩体和垛崮山岩体不同的，玲珑岩体中与华北板块物源有关的继承锆石(1750～2200Ma、～2500Ma、2700～2900Ma)数量更多。从该地区大地构造格局来看，文登岩体所在的苏鲁超高压变质地体属于扬子板块，垛崮山岩体位于昆嵛山边界杂岩带，属华北板块与扬子板块交界处(Zhaietal., 2000; Heetal., 2021)，玲珑岩体地处的胶北地体属于华北板块。从文登岩体→垛崮山岩体→玲珑岩体，代表华北板块物源信息的继承锆石(～2500Ma、～1800Ma)逐渐增多，这与三个岩体所处的大地构造单元位置相符合。

文登岩体处于苏鲁超高压变质地体内，其超高压变质岩的原岩是扬子板块的属性，因而在文登岩体中出现华北板块物源的记录有重要的构造意义。前人也在苏鲁超高压变质地体内发现多处太古代-古元古代的岩石具有华北板块的亲缘性(Zhangetal., 2006, 2014; Zhouetal., 2008)。目前普遍认为大别-苏鲁造山带的俯冲极性是扬子板块俯冲到华北板块的下部，如果在苏鲁超高压地体下部有华北板块的物质，而在胶北地体(传统上认为属于华北陆块)下部有扬子陆块的物质，那传统的单向俯冲构造样式就很难解释。而类似鳄鱼嘴式的重叠地壳构造模式(Li, 1994; 翟明国等，2000; Zhaoetal., 2016)或许能更好的解释胶东地区晚侏罗世岩体中的继承锆石特征。

将文登岩体、垛崮山岩体和玲珑岩体中新生岩浆锆石(～160Ma)的Hf同位素特征进行对比(图6)，结果显示：文登岩体εHf(t)值在-30.7～-20.3之间，峰值集中在-26; 垛崮山岩体εHf(t)值在-27.5～-17.9之间，峰值集中在-23；而玲珑岩体εHf(t)值变化范围较大，主要集中在-29.3～-9.3之间，而且明显呈现两个峰值，分别在-25和-17。从整体上来看，从文登岩体→垛崮山岩体→玲珑岩体，εHf(t)值有略微增大的趋势(图6)。在176Lu/177Hf对176Hf/177Hf图解中(图7)，玲珑岩体的Hf同位素组成也明显分成两个不同范围，即高εHf(t)值(玲珑岩体H)和低εHf(t)值(玲珑岩体L)的玲珑岩体能识别出两种不同的Hf同位素组成。文登岩体和垛崮山岩体的Hf同位素均与玲珑岩体中低εHf(t)值(玲珑岩体L)的特征一致，表明了三个岩体的源区具有一定的相似性。而玲珑岩体还具有高εHf(t)值(玲珑岩体H)，这反映了胶北地体的地壳深部存在多个源区物质的混合，除了扬子板块物质、华北板块物质的重熔外，可能有部分幔源物质的加入，使锆石的εHf(t)值变大。Sr、Nd等同位素数据也表明，胶东地区晚侏罗世的岩浆活动主要来源于古老地壳融熔并伴随有幔源物质的加入(Dengetal., 2018)。此外，同样位于胶北地体内的早白垩世郭家岭岩体为壳幔混合成因(Yangetal., 2012, 2014; Wangetal., 2014; 宋英昕等，2020)，高εHf(t)值的玲珑岩体与郭家岭岩体具有相似的锆石Hf同位素组成(图7)，这也表明这类玲珑岩体的形成可能具有幔源物质的参与。在胶东地区，从东到西(即文登岩体→垛崮山岩体→玲珑岩体)金矿床逐渐增多，而且大量的金矿集中产在玲珑岩体所处的胶北地体中。虽然从时间上来看，金成矿作用主要集中在130～110Ma(Yang and Zhou, 2001; 陈衍景等，2004；Dengetal., 2015, 2020; 朱日祥等，2015；范宏瑞等，2021)，与晚侏罗世岩浆活动相隔至少30Ma，金成矿不可能由该期岩浆活动直接作用形成。但是，胶北地区晚侏罗世岩浆活动的多源区特点，尤其是幔源物质的加入能够促使金在深部预富集，为之后金在早白垩世短期巨量成矿提供了一个重要的物质条件。

在晚侏罗世时期，郯庐断裂开始发生左行走滑活动(张长厚等，2001；朱光等，2009)，该深达岩石圈地幔的断裂为壳/幔物质和能量交换提供了的重要场所(Zhaoetal., 2012；杨立强等，2019)。与胶东地区东部相比，胶北地体更靠近郯庐断裂带，显然更容易受到其构造活动的影响。幔源岩浆底侵诱发加厚的下地壳拆沉并发生部分融熔，是晚侏罗世花岗岩形成的动力学机制之一(Houetal., 2007; Maetal., 2013; Yangetal., 2018)。这些幔源岩浆的入侵不仅促使了胶北地体太古宙-元古宙绿岩带地体强烈活化改造，使老地体中的金及成矿相关的元素被激活(李洪奎等，2011)，而且相对富水、富金的幔源岩浆(Xiaetal., 2019; Wangetal., 2020)还可以形成富含金及硫化物的角闪岩堆积体(Davidsonetal., 2007; Houetal., 2017)，这些过程使金元素在下地壳深部完成初始富集，为之后在早白垩世爆发式金成矿提供了重要的物质基础(范宏瑞等，2021)。而胶东东部地区可能缺乏这种金的预富集机制，因而金矿较为贫乏。

5 结论

(1)胶东地区东部的文登岩体和垛崮山岩体均形成于晚侏罗世(～160Ma)，与胶北金矿区的玲珑岩体形成时代一致。两个岩体中有较多的三叠纪(～230Ma)和新元古代(～780Ma)的继承锆石，表明它们的岩浆源区主要为经历了超高压变质作用的扬子板块。垛崮山岩体中还含有一定数量的华北板块物源的年龄记录，表明华北板块对其物质源区有重要的贡献。

(2)文登岩体和垛崮山岩体的新生岩浆锆石εHf(t)值范围分别为-30.7～-20.3和-27.5～-17.9，玲珑岩体与东部这两个岩体相比，其新生岩浆锆石εHf(t)值的变化范围更大，在-29.3～-9.3之间，显示出源区更加复杂，εHf(t)值明显偏高暗示可能有更多地幔物质的影响。

(3)胶东地区晚侏罗世花岗岩εHf(160Ma)值从东到西呈现增高的趋势，可能的成因是三叠纪扬子板块与华北板块碰撞后，加厚的地壳在晚侏罗世部分熔融，形成花岗岩。在此过程中胶北地体深部受到幔源物质的影响，从而对金及成矿相关元素产生一定的预富集作用。而同时期胶东东部地区未受到明显的地幔物质影响，可能是该区金矿产出较少的原因之一。

致谢成文过程中，邓军教授对本文提出了许多建设性意见；两位审稿人对本文提出了十分宝贵的修改意见，提高了文章的质量；在此一并表示衷心的感谢！