地球早期大陆地壳的形成与演化一直是前寒武纪地质学研究的热点问题。在过去的20年中，地球早期演化历史研究主要集中在两个方面：（1）古大陆的增生与克拉通化过程；（2）板块构造起始过程[1]。华北克拉通作为地球早期陆壳的重要组成部分，与世界其他克拉通以～2.7 Ga的岩浆活动为主要活动期不同，华北克拉通～2.5 Ga 的岩浆活动非常强烈[2,3]。对华北克拉通独特和复杂的演化过程，不同学者存在不同认识，但大致可归纳为以下两种主要观点：一种观点认为华北克拉通于～2.5 Ga完成克拉通化，以微陆块拼合的“五台运动”为标志[4-6]；另一种观点则认为华北克拉通基底在太古宙末尚未完全固结，后期在古元古代存在古陆块或弧陆碰撞[7-9]，～1.85 Ga东部陆块和西部陆块沿中部造山带碰撞拼合，最终完成克拉通化[10]。显然，有关华北克拉通基底拼合时间及方式仍是前寒武纪地质研究的焦点问题，而～2.5 Ga的岩浆事件性质的研究无疑会提供华北克拉通早前寒武纪地质演化的重要信息。
华北克拉通～2.5 Ga的岩浆事件的一个显著特点是TTG岩浆活动的末期普遍伴有钾质花岗岩的侵位，如鞍山地区齐大山钾质花岗岩、冀东地区山海关钾质花岗岩、鲁西地区四海山及鲁山钾质花岗岩、嵩山地区路家沟钾质花岗岩、怀安地区钾质花岗岩以及赞皇杂岩中的菅等钾质花岗岩等[11，12]。邓晋福等[4]认为华北克拉通新太古代末期花岗岩类为T1T2G1G2组合，代表半成熟陆壳的组成。晋冀蒙交界一带，沿山西繁峙县北东一直至河北宣化县分布一条钾质花岗岩带，北东东向断续延伸近200 km，被认为是桑干构造带的重要组成部分。以往一般认为该条岩浆岩带具有碰撞和后碰撞花岗岩的特征[13]，是怀安古陆和恒山古陆之间碰撞带中岩石部分熔融的产物，在古陆的拼合或大陆增长中起着某种作用[14]。但也有研究者认为，所谓怀安古陆和恒山古陆原来可能曾是一个古陆，目前的差异只是它们出露的深度不同，代表了不同的地壳水平，而桑干构造带以及其中的高压变质岩石和花岗岩带是大陆地壳基底逆掩和大规模滑脱作用的产物，且岩浆活动可能从太古宙末（～2.5 Ga）一直持续到古元古代中期（～2.15 Ga）[15]。另外，张华锋[16]等则获得董家沟一带黑云母二长花岗岩岩与石榴正长花岗岩的锆石U-Pb 年龄分别为～2 440 Ma和～2 000 Ma，并认为这两期钾质花岗岩岩浆事件分别与新太古代末微陆块拼合和古元古代板内伸展或裂解有关。在怀安地区，近年来研究者们对～2.5 Ga 岩浆作用的研究更多的偏重于TTG 方面[17-19]，而对与TTG片麻岩密切伴生的二长花岗岩的研究则相对薄弱。为此，本文对天镇-怀安地区二长花岗岩进行了岩相学、地球化学、锆石U-Pb定年等工作，探讨其岩石成因、时代及构造环境，以期对晋冀蒙一带新太古代末期构造演化提供约束。
1 地质背景及岩石学特征
晋冀蒙交界地区是华北克拉通典型的高级变质岩区，早前寒武纪变质基底大致以大同-兴和一线为界划分为两套高级变质岩系，南东一侧为太古代怀安杂岩，北西一侧为古元古代孔兹岩系。研究区位于怀安杂岩中的天镇县一带，主体地质体为新太古代变质深成片麻岩，出露面积约占变质基底总面积的80%以上，主要由包括闪长质片麻岩、英云闪长质片麻岩、奥长花岗质片麻岩的TTG组合以及少量的二长花岗岩组成，原岩年龄介于2.45～2.55 Ga，变质年龄为1.82～1.85 Ga[16-18,20,21]。此外，区内出露少量新太古代硅铁建造、榴云片麻岩以及古元古代高压基性麻粒岩-大理岩、富铝片麻岩等表壳岩组合[22]。还见大量古元古代变质基性岩墙，侵入新太古代TTG片麻岩杂岩中（图1）。
图1 晋冀蒙交界地区（b）和天镇-怀安地区（c）早前寒武纪地质简图及华北克拉通构造分区图（a）（图b、c据参考文献[22]修改，图a据参考文献[20]修改）
Fig.1 Precambrain geological sketch of the border area of Shanxi,Heibei and Inner Mongolia Provinces(b)and Tianzhen-Huai’an area(c)and tectonic subdivision of the North China Craton（a）
1.新太古代TTG片麻岩；2.新太古代二长花岗岩；3.新太古代条带状硅铁建造；4.新太古代榴云片麻岩岩组；5.古元古代高压基性麻粒岩-大理岩组合；6.古元古代富铝片麻岩；7.中元古代盖层；8.古元古代二辉麻粒岩（变质基性岩墙）；9.中元古代基性岩墙；10.岩墙型高压基性麻粒岩出露点；11.片麻理/地层产状；12.构造接触界线；13.断层；14.同位素测年及岩石化学分析样品采样点；15.岩石化学分析样品采样点
区内的二长花岗岩出露面积有限，但分布较为广泛，以天镇东南石厂沟-白家烟一带和怀安县南部赵家窑-太平沟一带出露规模相对较大，多呈不规则小岩株状、脉状或枝脉状侵入TTG片麻岩中，片麻状构造不明显（图2a），内部见有基性麻粒岩呈“夹层”状产出（图2b）；在一些构造变形较强的区域，则见二长花岗岩与英云闪长质片麻岩彼此交互产出，二者接触界线平直，产状与英云闪长质片麻岩片麻理协调一致，表现出深熔混合岩化的特征。
二长花岗岩岩性整体均匀，片麻理不发育，岩石呈肉红色，中细粒花岗变晶结构，块状-弱片麻状构造（图2c、d）。主要由斜长石、钾长石、石英及少量暗色矿物黑云母组成。斜长石呈它形粒状，粒度一般为0.5～3.0 mm，杂乱分布，含量30%～35%。钾长石呈它形粒状，主要为条纹长石，部分岩石含微斜长石，粒度一般为1.0～3.0 mm，杂乱分布，少部分与石英呈文象结构，部分粒内嵌布少量斜长石、石英等颗粒，含量35%～40%。石英呈它形粒状，粒度一般为1.0～2.0 mm，杂乱分布，部分与钾长石呈文象结构，粒内具轻微波状消光，含量20%～25%。暗色矿物多被绿泥石、次闪石等交代呈假像，局部残留为黑云母，星散分布，含量一般在5%以下，局部含量较高，可达5%～10%。
2 分析方法
测年样品锆石分选工作由廊坊市宇能岩石矿物分选技术服务有限公司完成。样品按照常规粉碎淘洗，经磁选和重液分离，然后在双目镜下人工挑选纯度在99%以上的锆石。其中一件测年样品的锆石制靶和透射光、反射光、阴极发光照相及样品测年工作在中国地质调查局天津地质调查中心完成，采用LAMC-ICP-MS 进行锆石U-Pb 同位素定年测试，ICPMS为Agilent 7500a，分析中采用的激光束斑直径为35 μm，频率为8 Hz，能量为5 mJ。分析时采用GJ-1作为年龄外标，NIST610作为元素含量外标，分析流程见文献[23]。最终测试数据的离线处理采用软件ICPMASDataCal[24]，U-Pb年龄谐和图绘制和年龄权重平均计算均采用ISOPLOT3.0程序[25]；另外一件样品的锆石制靶和透射光、反射光、阴极发光照相在北京离子探针中心实验室完成，锆石U-Pb定年在北京离子探针中心SHRIMP Ⅱ上完成，详细分析方法见文献[26]。测试时一次流O2-强度为3～5 nA，束斑直径为25 μm。标样M257（U=840×10-6）[27]和TEM（年龄为417 Ma）[28]分别用于锆石U 含量和年龄校正。TEM与未知样品测定比例为1/3～1/4。对于样品和标准锆石数据点测定均由5组扫描给出。数据处理采用SQUID和ISOPLOT程序[25]。
样品主微量元素测试分析在中国地质调查局天津地质调查中心实验室完成。主量元素用X射线荧光光谱法（XRF）测试，FeO应用氢氟酸-硫酸溶样、重铬酸钾滴定的容量法，分析精度优于2%，微量元素使用ICP-MS测试，分析精度优于5%。
3 分析结果
3.1 锆石U-Pb年代学
本次工作选择两件二长花岗岩样品进行了年代学测定，样品岩性特征见图2a、b。其中，样品17DMTS-1采用LA-MC-ICPMS 锆石U-Pb测年；样品PM08TW53-1采用SHRIMP 锆石U-Pb定年，测试分析结果见表1、表2。
表1 二长花岗岩样品（17DMTS-1）LA-MC-ICPMS 锆石U-Pb 测年分析结果
Table 1 LA-MC-ICPMS Zircon U-Pb dating results of the monzonitic granite（17DMTS-1）
注：测试单位为中国地质调查局天津地质调查中心实验室
样品号17DMTS-1.1.1 17DMTS-1.1.2 17DMTS-1.2 17DMTS-1.3 17DMTS-1.4 17DMTS-1.5 17DMTS-1.6 17DMTS-1.7 17DMTS-1.8 17DMTS-1.9 17DMTS-1.10 17DMTS-1.11 17DMTS-1.12 17DMTS-1.13 17DMTS-1.14 17DMTS-1.15 17DMTS-1.16 17DMTS-1.17 17DMTS-1.18 17DMTS-1.19 17DMTS-1.20 17DMTS-1.21 17DMTS-1.22 17DMTS-1.23 17DMTS-1.24 17DMTS-1.25 17DMTS-1.26 17DMTS-1.27 17DMTS-1.28 17DMTS-1.29 17DMTS-1.30 17DMTS-1.31 17DMTS-1.32含量/（×10-6）Pb 27 10 32 10 41 11 11 22 22 17 20 35 9 30 9 16 29 18 37 21 12 935 0 32 86 107 16 13 4 27 99 27 U 225 286 169 130 45 124 23 440 6 339 334 78 56 564 40 88 8 152 123 373 169 78 44 247 9 525 184 261 148 52 249 130 8 Th 387 220 487 136 493 137 69 524 236 401 212 102 155 298 64 205 323 243 229 343 284 126 54 619 358 1039 299 289 178 90 333 618 303 Th/U 0.58 0.06 0.14 0.47 0.50 0.38 0.58 0.12 0.04 0.19 0.91 0.68 0.04 0.37 0.05 0.03 0.14 0.15 0.32 0.55 0.62 0.62 0.35 0.46 0.07 0.37 0.24 0.53 0.50 0.27 0.28 0.46 0.26同位素比值206Pb/238U 0.459 8 0.281 8 0.462 2 0.446 8 0.364 2 0.452 0 0.451 3 0.447 9 0.309 4 0.460 8 0.459 2 0.420 1 0.332 6 0.471 3 0.284 1 0.282 4 0.358 9 0.367 5 0.482 5 0.467 4 0.467 6 0.350 7 0.455 6 0.445 4 0.329 3 0.470 4 0.471 9 0.462 9 0.467 7 0.326 8 0.474 0 0.472 5 0.475 1 1σ 0.010 7 0.002 9 0.004 7 0.004 7 0.003 9 0.008 0 0.005 4 0.004 7 0.003 1 0.004 9 0.004 9 0.004 1 0.004 0 0.005 8 0.003 5 0.003 0 0.004 0 0.005 8 0.004 9 0.005 1 0.004 7 0.003 6 0.004 6 0.004 7 0.003 6 0.005 1 0.005 4 0.006 0 0.005 7 0.003 7 0.005 3 0.005 0 0.004 9 207Pb/235U 10.179 8 5.063 0 10.150 8 9.401 6 7.469 6 9.470 1 9.641 3 9.939 6 4.780 8 10.362 1 10.199 8 9.026 3 5.196 1 10.625 8 5.012 8 4.507 0 7.754 0 6.330 3 10.832 0 10.696 6 10.486 4 7.152 7 9.949 6 9.439 4 5.236 8 10.570 0 10.572 9 10.130 2 10.273 5 6.298 7 10.619 7 10.491 2 10.566 5 1σ 0.255 0 0.070 7 0.138 3 0.130 3 0.102 2 0.197 6 0.144 3 0.135 6 0.064 2 0.142 9 0.142 3 0.117 0 0.078 2 0.158 2 0.077 3 0.062 2 0.110 7 0.141 6 0.144 9 0.149 9 0.137 7 0.099 9 0.134 7 0.131 4 0.079 3 0.148 5 0.156 9 0.159 1 0.157 1 0.099 4 0.151 7 0.147 1 0.145 1 207Pb/206Pb 0.160 6 0.130 3 0.159 3 0.152 6 0.148 7 0.152 0 0.154 9 0.161 0 0.112 1 0.163 1 0.161 1 0.155 8 0.113 3 0.163 5 0.128 0 0.115 7 0.156 7 0.124 9 0.162 8 0.166 0 0.162 6 0.147 9 0.158 4 0.153 7 0.115 3 0.163 0 0.162 5 0.158 7 0.159 3 0.139 8 0.162 5 0.161 0 0.161 3 1σ 0.000 4 0.001 1 0.000 4 0.000 9 0.000 3 0.001 2 0.000 9 0.000 5 0.000 6 0.000 7 0.000 6 0.000 5 0.001 1 0.000 7 0.001 1 0.000 6 0.000 4 0.000 6 0.000 3 0.001 1 0.000 8 0.001 7 0.003 6 0.000 2 0.000 3 0.000 2 0.000 3 0.000 6 0.000 8 0.002 2 0.000 4 0.000 3 0.000 4年龄/Ma 206Pb/238U 2 439 1 600 2 449 2 381 2 002 2 404 2 401 2 386 1 738 2 443 2 436 2 261 1 851 2 489 1 612 1 604 1 977 2 017 2 538 2 472 2 473 1 938 2 420 2 375 1 835 2 485 2 492 2 452 2 473 1 823 2 501 2 495 2 506 1σ 57 17 25 25 22 43 29 25 18 26 26 22 22 30 20 17 22 32 26 27 25 20 24 25 20 27 28 32 30 20 28 26 26 207Pb/235U 2 451 1 830 2 449 2 378 2 169 2 385 2 401 2 429 1 782 2 468 2 453 2 341 1 852 2 491 1 821 1 732 2 203 2 023 2 509 2 497 2 479 2 131 2 430 2 382 1 859 2 486 2 486 2 447 2 460 2 018 2 490 2 479 2 486 1σ 61 26 33 33 30 50 36 33 24 34 34 30 28 37 28 24 31 45 34 35 33 30 33 33 28 35 37 38 38 32 36 35 34 207Pb/206Pb 2 462 2 102 2 448 2 375 2 332 2 368 2 401 2 466 1 833 2 488 2 467 2 411 1 853 2 492 2 070 1 891 2 420 2 028 2 485 2 517 2 483 2 322 2 439 2 387 1 885 2 487 2 482 2 442 2 448 2 224 2 482 2 466 2 469 1σ 21 22 21 21 21 21 21 21 23 21 21 21 22 21 22 22 21 24 21 21 21 22 21 21 23 21 21 21 21 22 21 21 22
表2 二长花岗岩样品（PM08TW53-1）锆石SHRIMP U-Pb 测年分析结果
Table 2 SHRIMP Zircon U-Pb dating results of the monzonitic granite（PM08TW53-1）
注：测试单位为北京离子探针中心实验室。
样品号PM08TW53-1-1.1 PM08TW53-1-2.1 PM08TW53-1-2.2 PM08TW53-1-3.1 PM08TW53-1-4.1 PM08TW53-1-4.2 PM08TW53-1-5.1 PM08TW53-1-6.1 PM08TW53-1-6.2 PM08TW53-1-7.1 PM08TW53-1-8.1 PM08TW53-1-9.1 PM08TW53-1-8.2 PM08TW53-1-10.1 PM08TW53-1-11.1 PM08TW53-1-10.2 PM08TW53-1-12.1 PM08TW53-1-13.1 PM08TW53-1-14.1 206Pbc×10-6 0.05 0.44 1.35 0.20 0.13 0.17 0.83 0.24 1.91 0.61 0.22 0.56 0.68 0.33 0.19 0.07 0.06 0.32 0.26 U×10-6 172 44 1 478 321 192 1 185 1 316 529 1 652 1 318 209 1 188 1 361 170 552 286 598 171 191 Th 179 17 69 106 110 89 465 49 97 649 111 592 105 52 351 60 64 118 174 Th/U 1.07 0.39 0.05 0.34 0.59 0.08 0.36 0.10 0.06 0.51 0.55 0.52 0.08 0.32 0.66 0.22 0.11 0.71 0.94同位素比值206Pb*/238U 0.487 0 0.480 9 0.254 5 0.467 6 0.460 0 0.268 7 0.168 4 0.454 5 0.304 2 0.367 8 0.482 0 0.210 9 0.358 6 0.504 1 0.451 5 0.510 3 0.466 8 0.348 9 0.348 5%1.5 1.8 1.3 1.5 1.4 1.3 1.3 1.4 1.3 1.4 1.5 1.3 1.3 1.5 1.4 1.4 1.4 1.5 1.5 207Pb*/235U 11.16 10.65 4.542 10.36 10.00 4.641 2.144 9.79 5.702 7.49 10.68 3.022 7.157 11.38 9.91 11.60 10.20 5.362 5.215%1.7 2.5 1.5 1.7 1.6 1.5 1.6 1.4 1.6 1.4 1.6 1.5 1.4 2.0 1.5 1.5 1.4 1.8 1.8 207Pb*/206Pb*0.166 2 0.160 7 0.129 5 0.160 6 0.157 7 0.125 28 0.092 34 0.156 17 0.136 0 0.147 83 0.160 79 0.103 95 0.144 80 0.163 7 0.159 28 0.164 87 0.158 44 0.111 5 0.108 6%0.80 1.7 0.77 0.71 0.66 0.60 0.93 0.45 0.89 0.42 0.61 0.69 0.44 1.3 0.56 0.50 0.47 1.1 1.1误差相关性0.880 0.741 0.876 0.907 0.909 0.912 0.824 0.948 0.854 0.959 0.923 0.893 0.952 0.738 0.925 0.943 0.945 0.811 0.791年龄/Ma 206Pb/238U 2 558 2 531 1 461 2 473 2 439 1 534 1 003 2 415 1 712 2 019 2 536 1 234 1 975 2 631 2 402 2 658 2 470 1 929 1 927 1σ±31±38±17±32±29±18±12±27±20±24±31±15±23±31±27±31±28±24±24 207Pb/206Pb 2 520 2 463 2 091 2 462 2 431 2 033 1 474 2 415 2 176 2 321 2 464 1 696 2 285 2 494 2 448 2 506 2 439 1 823 1 775 1σ±14±28±13±12±11±11±17±8±14±7±10±12±7±22±94±9±8±19±21 Disc/%-2-3 30 002 5 32 0 21 13-3 27 14-6 2-6-1-6-9
图2 天镇-怀安地区二长花岗岩野外露头及显微照片
Fig.2 Field outcrop and microphotographs of the monzogranites in the Tianzhen-Huai’an area
a.二长花岗岩露头特征；b.“夹层”状产于二长花岗岩中的基性麻粒岩；c.强变形带处二长花岗岩平行英云闪长质片麻岩片麻理并与之交互产出；d.二长花岗岩显微照片
样品17DMTS-1 分选锆石自形程度相对较高，多呈长柱状，长宽比一般在2/1～3/1之间。锆石CL影像多具核边结构，核部结晶环带保留较好，为岩浆成因，变质边相对较窄，颜色较为均匀（图3a）。样品共分析32 个数据点，其中，核部结晶环带较好的锆石207Pb/206Pb 表面年龄多在2 400～2 500 Ma 之间（Th/U比值均＞0.1），部分数据点可能受变质作用影响，测试数据为混合年龄，结果偏小；还有部分数据点存在Pb丢失，落在谐和线下方（图4a）。选择在谐和线上方相对集中的17个数据点进行加权平均，获得的成岩年龄为2 472±10 Ma。另外，样品中4个锆石变质边分析点的Th/U小于0.07，207Pb/206Pb表面年龄在1 891～1 833 Ma之间，代表古元古代变质事件的年龄。
样品PM08TW53-1分选锆石多呈长柱状，其长宽比一般在2/1～3/1之间，在CL图像上，锆石具有窄的变质边，但核部保留较好的环带构造，显示岩浆成因（图3b）。样品通过SHRIMP 锆石U-Pb定年，共分析19个数据点，其中的10个数据点相对集中的落在谐和线同一位置附近（图4b），锆石207Pb/206Pb表面年龄多在2 414～2 520 Ma之间（Th/U比值均＞0.1），选择其中6个谐和度较高的数据点的加权平均年龄为2 448±9 Ma，代表岩石形成时代。
图3 二长花岗岩锆石阴极发光图像
Fig.3 CL images of selected zircons of the monzogranites
图4 二长花岗岩锆石U-Pb年龄谐和图
Fig.4 Zircon U-Pb concordia diagrams of the monzogranites
3.2 岩石地球化学特征
对所采集的岩石样品进行了主量元素、稀土及微量元素的测试，详细分析结果及各类参数见表3。
表3 二长花岗岩主量元素（%）、微量元素（×10-6）分析结果
Table 3 Major（%）and trace（×10-6）elements result of the monzonitic granite
样品编号SiO2 TiO2 Al2O3 Fe2O3 FeO MnO MgO CaO Na2O K2O P2O5 LOI Total K2O/Na2O A/CNK σ Mg#DI Pb Cr Ni Co Rb Sr Ba Sc Nb Ta Zr Hf Ga U Th La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Y ΣREE LREE HREE LREE/HREE LaN/YbN PM08Tw53-1 73.31 0.2 13.78 0.8 0.8 0.021 0.42 1.32 3.04 5.74 0.066 0.41 99.50 1.89 1.01 2.54 0.33 90.2 18.8 5.88 5.28 3.03 184 289 1200 2.85 1.84 0.034 214 6.42 16.3 0.46 20.5 53 90.1 8.96 29 3.5 1.2 2.7 0.24 0.65 0.095 0.25 0.034 0.22 0.04 2.3 189.99 185.76 4.23 43.93 172.80 PM08Xt.Wl.Gs53-1 73.06 0.19 13.89 0.44 1.13 0.02 0.46 1.33 3.12 5.72 0.061 0.47 99.42 1.83 1.01 2.60 0.35 89.73 19 5.91 5.72 2.77 166 294 1120 3 2.36 0.034 176 5.33 15.1 0.44 19.3 44.4 76.4 7.66 25.2 3.21 1.12 2.5 0.23 0.7 0.1 0.25 0.035 0.24 0.039 2.55 162.08 157.99 4.09 38.59 132.70 TW8038-1 73.91 0.16 13.99 0.98 0.18 0.014 0.34 1.02 3.46 5.31 0.036 0.57 99.40 1.53 1.05 2.49 0.37 91.91 17.5 3.57 2.52 1.93 173 223 857 8.95 1.54 0.032 122 4.27 17.2 0.7 11.8 26 41 4.17 13.2 1.6 0.81 1.4 0.12 0.42 0.068 0.18 0.026 0.19 0.033 1.98 89.22 86.78 2.44 35.61 98.16 17DMTS-1 74.29 0.16 14.09 0.3 0.54 0.008 0.18 1.15 3.43 5 0.033 0.78 99.18 1.46 1.07 2.27 0.29 91.74 22.8 2.2 1.92 0.8 210 178 704 9.95 3.37 0.031 128 4.53 12.6 2.82 25.9 39.3 69.3 7.13 22.6 2.88 0.54 2.42 0.22 0.74 0.12 0.31 0.042 0.29 0.051 3.08 145.94 141.75 4.19 33.81 97.21 SGS-1 73.71 0.16 14.51 0.7 0.13 0.011 0.15 1.23 4.58 4.26 0.031 0.52 99.47 0.93 1.01 2.54 0.26 92.35 23.7 3.7 1.79 0.77 143 110 212 9.22 11.3 0.46 70.6 2.82 20.6 2.24 13.8 12.7 22.2 2.73 9.62 1.66 0.32 1.37 0.15 0.64 0.11 0.26 0.033 0.2 0.028 2.7 52.02 49.23 2.79 17.64 45.55 TW8060-1 73.46 0.23 14.02 1.23 0.38 0.03 0.31 1.23 3.45 5.12 0.06 0.42 99.51 1.48 1.04 2.41 0.27 90.49 23.8 3.28 2.14 2.37 213 210 958 9.9 8.7 0.48 174 5.49 17.5 1.88 22 49.3 86.4 8.82 27.6 3.66 0.74 3.24 0.34 1.48 0.26 0.74 0.11 0.8 0.12 7.36 183.61 176.52 7.09 24.90 44.20 YQ7260-1 70.55 0.32 15.35 1.69 0.38 0.02 0.47 1.37 3.68 5.11 0.10 0.93 99.03 1.39 1.09 2.80 0.31 88.19 17.6 7.3 4.8 4.1 139 382 1200 10.7 3 0.11 242 6.76 18.4 0.51 24 69.5 116 11.8 36.2 3.87 1 3.35 0.24 0.66 0.096 0.24 0.03 0.21 0.036 2.63 243.232 238.37 4.86 49.03 237.39
3.2.1 主量元素特征
二长花岗岩样品SiO2含量在70.55%～74.29%之间，平均为73.18%。Al2O3含量在13.78%～15.35%之间，平均为14.23%。Na2O 含量在3.04%～4.58%之间，平均为3.54%。K2O含量在4.26%～5.74%之间，平均为5.18%。除一个样品K2O/Na2O 小于1 外（0.93），其余样品K2O/Na2O在1.39～1.89之间，平均为1.60，属钾质岩石系列。CaO 含量在1.02%～1.37%之间，平均为1.24%。MgO 含量低，介于0.15%～0.47%之间，平均为0.33%，Mg#值介于26～37之间，平均值为31。TiO2含量在0.16%～0.32%之间，平均为0.20%。岩石整体相对富Si、Na、K、Al，贫Ca、Ti、Mg。样品的A/NCK在1.01～1.09之间，平均为1.04，属过铝质岩石；分异指数（DI）为88.19～92.35，分异程度高；里特曼指数（σ43）为2.27～2.8，为钙碱性岩。
在侵入岩TAS图解中（图5），样品投点均落于花岗岩区内，属亚碱性系列；在SiO2-K2O图解中（图6），样品主要落在高钾钙碱性系列中。
图5 二长花岗岩TAS图解（底图据参考文献[29]）
Fig.5 TAS diagram of the monzogranites
Ir.分界线，上方为碱性，下方为亚碱性；1.橄榄辉长岩；2a.碱性辉长岩；2b.亚碱性辉长岩；3.辉长闪长岩；4.闪长岩；5.花岗闪长岩；6.花岗岩；7.硅英岩；8.二长辉长岩；9.二长闪长岩；10.二长岩；11.石英二长岩；12.正长岩；13.副长石辉长岩；14.副长石二长闪长岩；15.副长石二长正长岩；16.副长正长岩；17.副长深成岩；18.霓方钠岩/磷霞岩/粗白榴岩
图6 二长花岗岩SiO2-K2O图解（底图据参考文献[30]）
Fig.6 SiO2-K2O diagram of the monzogranites
3.2.2 稀土、微量元素特征
二长花岗岩的稀土总量变化较大，ΣREE 在52.02×10-6～243.23×10-6之间，平均值为152.30×10-6，含量整体偏低。LREE 为49.23×10-6～238.37×10-6，平均值148.06×10-6，HREE 为2.44×10-6～7.09×10-6，平均值4.24×10-6。在球粒陨石标准化稀土元素配分模式图中（图7a），稀土配分曲线整体呈轻稀土相对富集、重稀土亏损的右倾分布形式，LREE/HREE为17.64 × 10- 6～49.03 × 10- 6，（La/Yb）N 为44.20～237.39。铕异常变化较大，其中的4件样品具弱的负铕异常，δEu介于0.63～0.85之间；另外3件样品为正铕异常，δEu介于1.19～1.65之间。在原始地幔标准化微量元素蛛网图中（图7b），样品整体显示富集Rb、Ba、K、Sr等大离子亲石元素，亏损U、Nb、Ta、P、Ti等高场强元素的特点。另外，样品Sr含量变化较大，在110×10-6～382×10-6之间（平均值240.86×10-6）；Y、Yb含量低，其中Y含量在1.98×10-6～7.36×10-6之间（平均值3.23×10-6）；Yb含量在0.19×10-6～0.8×10-6（平均值0.31×10-6）。过渡元素中Cr、Ni 含量低，分别为2.2×10-6～7.3×10-6（平均值4.55×10-6）、1.79×10-6～5.72×10-6（平均值3.45×10-6）。
图7 二长花岗岩球粒陨石标准化稀土配分模式图（a）和原始地幔标准化微量元素蛛网图（b）（标准化值据参考文献[31]）
Fig.7 Chondrite-normalized REE distribution patterns(a)and primitive mantle-normalized trace element patterns(b)of the monzogranites
4 讨论
4.1 成岩时代
前人对怀安杂岩中TTG片麻岩开展了大量同位素年代学工作，Zhao等[20]获得蔓菁沟一带英云闪长岩、奥长花岗岩的SHRIMP锆石U-Pb年龄分别为2 515±20 Ma、2 499±19 Ma；刘富等[17]获得天镇瓦窑口一带英云闪长质片麻岩的LA-ICP-MS锆石U-Pb年龄为2 502±11 Ma，怀安东洋河一带奥长花岗质片麻岩的LA-ICP-MS 锆石U-Pb 年龄为2 506±13 Ma；Wang 等[21]获得蔓菁沟一带英云闪长质片麻岩LAICP-MS 锆石U-Pb 年龄2 503±17 Ma。而有关怀安杂岩中新太古代末期钾质花岗岩锆石年龄的报道相对要少，Zhang 等[16]获得董家沟黑云母二长花岗岩LA-ICP-MS锆石U-Pb年龄为2 437±10 Ma，本次工作获得天镇大馒头山一带二长花岗岩LA-ICP-MS和SHRIMP锆石U-Pb年龄分别为2 472±10 Ma、2 448±9 Ma。由上述测年结果结合野外地质关系看，怀安地区新太古代末钾质花岗岩成岩年龄略晚于TTG片麻岩，二者应是同一地质事件不同阶段岩浆作用的产物。
4.2 岩石成因
研究区二长花岗岩富Si、高K，且样品K2O/Na2O基本都大于1，属钾质岩石，A/NCK介于1.01～1.09之间，属过铝质岩石。样品富集LREE且多数Th含量高，并具有与中、上部地壳非常相似的微量元素蛛网图曲线，显示出明显的地壳重熔特征。另外，样品低Mg（MgO 含量介于0.15%～0.47%之间，平均为0.33%，Mg#值介于26～37之间，平均值为31），Cr、Ni 含量低（分别为2.2×10- 6～7.3×10- 6、1.79×10-6～5.72×10-6），样品的Rb/Sr 比值≥0.36，接近或高于地壳平均值（0.35）[32]，而样品强烈亏损Nb、Ta，且Nb/Ta比值高，在18.13～108.71 之间，平均为50.05，高于地幔平均值（约17.5）[31]，说明岩浆源区以壳源物质为主，无明显地幔物质混入。同时，样品Sr含量变化大，在110×10-6～382×10-6，平均为240.86×10-6，但Y、Yb 含量很低，分别在1.98×10-6～7.36×10-6、0.19×10-6～0.8×10-6之间，在Sr-Yb判别图解中（图8），样品落在Ⅱ区（低Sr低Yb型）内，指示形成于中高压力条件下，相应的地壳深度大致在40～50 km。综合上述认识，区内二长花岗岩应是加厚下地壳部分熔融的产物，而无明显幔源组分的加入。
图8 二长花岗岩Sr-Yb分类图解（底图据参考文献[33]）
Fig.8 Sr-Yb diagram of the monzogranites
Ⅰ.高Sr 低Yb 型（埃达克型花岗岩岩）；Ⅱ.低Sr 低Yb 型（喜马拉雅型花岗岩）；Ⅲ.高Sr 高Yb 型（广西型花岗岩）；Ⅳ.低Sr 高Yb 型（浙闽型花岗岩）；Ⅴ.非常低Sr 高Yb 型（南岭型花岗岩）
区内二长花岗岩产出形态主要呈两种样式：一种产出形态是呈小型侵入体侵入于TTG片麻岩中，岩石一般具很弱的片麻理；另一种产出形态是在强变形带处呈顺层脉状与TTG片麻岩交互产出，这些顺层脉体宽度从几厘米到几米不等，或者在TTG 片麻岩中呈几毫米的红色细窄脉体与暗色基体相间分布，使TTG片麻岩表现出明显深熔混合岩化的特征。对比混合岩化TTG片麻岩中脉体与二长花岗岩矿物组成，基本可以确定区内二长花岗岩主要是由TTG 岩石重熔形成。另外，大量高温高压实验成果表明TTG的局部熔融产生钾质花岗岩的上地壳[34]，多数研究者也认为华北克拉通2.5～2.45 Ga 的钾质花岗岩可能是由先期侵入下地壳的TTG发生重熔形成[16,35-38]。在怀
安杂岩中获得钾质花岗岩的Nd同位素单阶段模式年龄为2 587～2 596 Ma[16]，与怀安杂岩中TTG片麻岩中获得的Nd同位素2.5～2.6 Ga的单阶段模式年龄一致[17]，这也反映二者Nd 同位素方面的继承性。因此，可以推测怀安地区新太古代末钾质二长花岗岩也是由TTG岩石重熔而成。
续表3
注：测试单位为中国地质调查局天津地质调查中心实验室。Mg#=Mg2+/(Mg2++Fe2+)；里特曼指数σ=[ω(K2O+Na2O)]2/[ω(SiO2-43)]
样品编号PM08Tw53-1 PM08Xt.Wl.Gs53-1 TW8038-1 17DMTS-1 SGS-1 TW8060-1 YQ7260-1 δEu δCe Rb/Sr Nd/Ta 1.19 1.01 0.64 54.12 1.21 1.02 0.56 69.41 1.65 0.97 0.78 48.13 0.63 1.02 1.18 108.71 0.65 0.92 1.30 24.57 0.66 1.02 1.01 18.13 0.85 0.99 0.36 27.27
4.3 构造环境及形成机制
TTG片麻岩与二长花岗岩共同构成了怀安杂岩的主体，二者相近的形成时代以及成因上的密切联系也就限定了需要在同一构造模式下讨论新太古代末期的这期岩浆事件。虽然目前对TTG的成因尚存争议，但多数研究者认为华北克拉通新太古代末期TTG 岩石是与板块俯冲有关的岩浆活动的产物[20，37，39，40]。目前，一些研究也认为怀安地区的TTG岩石系可能来自新太古代俯冲板块的熔融[16-19]，但对于区内的钾质花岗岩，则认为其形成机制是由俯冲板片断裂引起软流圈物质上涌，导致先期侵入下地壳的TTG发生重熔而成[16,19]。但研究区的二长花岗岩相对较高的Sr含量与极低的Yb、Y含量指示形成压力条件较高，属低Sr、低Yb 的喜马拉雅型花岗岩[33]，不太可能形成于与底侵有关的伸展背景下，而更可能是地壳加厚阶段下地壳TTG 岩石熔融的产物。同时，在Rb/10-Hf-Ta×3图解中，样品投点落在碰撞大地背景上的花岗岩区域内（图9a），在Rb-（Y+Nb）图解中，样品投点基本都落在同碰撞花岗岩（Syn-COLG）区域内（图9b）。
图9 二长花岗岩Rb-Hf-Ta图解（底图据参考文献[41]）与Rb-Y＋Nb图解（底图据参考文献[42]）
Fig.9 Rb-Hf-Ta（a）and Rb-Y＋Nb（b）diagram of the monzogranites
综合上述分析，区内新太古代TTG片麻岩形成与俯冲板片的熔融有关，而二长花岗岩形成则可能与碰撞导致的地壳加厚引发的下地壳TTG岩石重熔有关。由此本文认为，天镇-怀安地区新太古代末期经历了由洋壳俯冲到陆块碰撞的构造演化过程，也标志着华北克拉通新太古代末期（约2.5～2.45 Ga）克拉通化的完成。
5 结论
（1）天镇-怀安地区二长花岗岩LA-ICP-MS 锆石U-Pb年龄为2 472±10 Ma，SHRIMP锆石U-Pb年龄为2 448±9 Ma，是新太古代末期岩浆活动结束的标志。
（2）二长花岗岩具有富Si、高K、贫Fe、Mg等常量元素特征，属钾质花岗岩，微量元素上富集Rb、Ba、K、Sr等大离子亲石元素，亏损U、Nb、Ta、P、Ti等高场强元素，同时Cr、Ni含量低，具有明显的地壳重熔特征。
（3）天镇-怀安地区TTG片麻岩和二长花岗岩是新太古代末期板块俯冲-碰撞不同阶段的产物，代表华北克拉通一次重要的陆壳水平增生事件。
[1]王惠初，于海峰，苗培森，等.前寒武纪地质学研究进展与前景[J].地质调查与研究，2011，34（04）:241-252.
[2]耿元生，沈其韩，任留东，等.华北克拉通晚太古代末-古元古代初的岩浆事件及构造热体制[J].岩石学报，2010，26（7）:1945-1966.
[3]Yang C H, Du L L, Ren L D, et al.Delineation of the ca.2.7Ga TTG gneisses in the Zanhuang Complex, North China Craton and its geological implications[J].Journal of Asian Earth Sciences,2013,72(10):178-189.
[4]伍家善，耿元生，沈其韩，等.中朝古大陆太古宙地质特征及构造演化[M].北京:地质出版社，1998:1-104.
[5]邓晋福，吴宗絮，赵国春，等.华北地台前寒武花岗岩类、陆壳演化与克拉通形成[J].岩石学报，1999，15（2）:190-198.
[6]Zhai M G, Xiao W J, Kusky T M, et al.Tectonic evolution of China and adjacent crustal fragments[J].Gondwana Research,2007,12(1-2):1-3.
[7]马杏垣，白瑾，索书田，等.中国前寒武纪构造格架及研究方法[M].北京:地质出版社，1987.
[8]孙大中，李惠民，林源贤，等.中条山前寒武纪年代学、年代构造格架和年代地壳结构模式的研究[J].地质学报，1991，65（3）:216-231.
[9]吴昌华.华北克拉通的变质沉积岩及其克拉通的构造划分[J].高校地质学报，2007，13（3）:442-457.
[10]Zhao G C, Cawood P A, Wilde S, et al.Review of global 2.1-1.8 Ga orogens: Implications for a pre-Rodinia supercontinent[J].Earth Science Reviews,2002,59(1):125-162.
[11]杨崇辉，杜利林，任留东，等.赞皇杂岩中太古宙末期菅等钾质花岗岩的成因及动力学背景[J].地学前缘，2011，18（2）:62-78.
[12]Wan Y S, Dong C Y, Liu D Y, et al.Zircon ages and geochemistry of late Neoarchean syenogranites in the North China Craton: A review[J].Precambrian Research, 2012,222-223:265-289.
[13]李永刚，郭敬辉，翟明国.天镇-怀安地区早前寒武纪钾质花岗岩地球化学特征及构造环境[J].华北地质矿产杂志，1995，10（2）:223-235.
[14]翟明国，李永刚，郭敬辉，等.晋冀内蒙交界地区麻粒岩地体中两条花岗岩带及其对早前寒武纪地壳生长的意义[J].岩石学报，1996，12（2）:299-314.
[15]郭敬辉，翟明国，李江海，等.华北克拉通早前寒武纪桑干构造带的岩石组合特征和构造性质[J].岩石学报，1996，012（002）:193-207.
[16]Zhang H F , Zhai M G , Santosh M , et al.Geochronology and petrogenesis of Neoarchean potassic meta- granites from Huai’an Complex:Implications for the evolution of the North China Craton[J].Gondwana Research, 2011, 20(1):82-105.
[17]刘富，郭敬辉，路孝平，等.华北克拉通2.5 Ga 地壳生长事件的Nd-Hf 同位素证据: 以怀安片麻岩地体为例[J].科学通报，2009，（17）:2517-2526.
[18]Liu F,Guo J H,Peng P,et al.Zircon U-Pb ages and geochemistry of the Huai’an TTG gneisses terrane: Petrogenesis and implications for ～2.5 Ga crustal growth in the North China Craton[J].precambrian research, 2012, 212-213(225-244).
[19]张华锋，王浩铮，豆敬兆，等.华北克拉通怀安陆块新太古代低铝和高铝TTG 片麻岩的地球化学特征与成因[J].岩石学报，2015，31（6）:1518-1534.
[20]Zhao G C, Wilde S A, Sun M , et al.SHRIMP U-Pb zircon geochronology of the Huai'an Complex: Constraints on Late Archean to Paleoproterozoic magmatic and metamorphic events in the Trans-North China Orogen[J].American Journal Ofence,2008,308(3):270-303.
[21]Wang J, Wu Y B, Gao S, et al.Zircon U-Pb and trace element data from rocks of the Huai’an Complex:New insights into the Late Paleoproterozoic collision between the eastern and western blocks of the North China Craton[J].Precambrian Research,2010,178(1-4):59-71.
[22]张家辉，王惠初，田辉，等.华北克拉通怀安杂岩中“MORB”型高压基性麻粒岩的成因及其构造意义[J].岩石学报，2018，35（11）:3506-3528.
[23]Yuan, H L, Gao S, Liu X M, et al.Accurate U-Pb Age and Trace Element Determinations of Zircon by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry[J].Geostandards and Geoanalytical Research, 2004, 28(3):353-370.
[24]Liu Y S, Gao S, Hu Z C, et al.Continental and Oceanic Crust Recyclinginduced Melt-Peridotite Interactions in the Trans-North China Orogen:U-Pb Dating,Hf Isotopes and Trace Elements in Zircons from Mantle Xenoliths[J].Journal of Petrology,2010,51(1-2):537-571.
[25]Ludwig K R.User’s Manual for Isoplot 3.00:A geochronological toolkit for microsoft excel[M].Berkeley Geochronology Center,Special Publication,2003,4:1-70.
[26]Williams I S.U-Th-Pb geochronology by ion microprobe[C].In: McKibben M A, Shanks W C III, Ridley W I, eds.Applications of Microanalytical Techniques to Understanding Mineralizing Processes.Reviews of Economic Geology:Society of Economic Geologists,1998,7:1-35.
[27]Nasdala L, Hofmeister w, Norberg N, et al.Zircon M257: A homogeneous natural reference material for the ion microprobe U-Pb analysis of zircon[J].Geostandards and Geoanalytical Research,2008,32(3):247-265.
[28]Black L P, Kamo S L, Allen C M, et al.TEMORA 1: a new zircon standard for Phanerozoic U- Pb geochronology[J].Chemical Geology,2003,200(1):155-170.
[29]Middlemost E A K.Naming materials in the magma/igneous rock system[J].Earth Science Reviews, 1994, 37(3-4):215-224.
[30]Peccerillo A, Taylor S R.Geochemistry of eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey[J].Contributions to Mineralogy & Petrology, 1976,58(1):63-81.
[31]Sun S S,McDonough W F.Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes[J].Geological Society,London,Special Publications,1989,42(1):313-345.
[32]Taylor S R, Mclennan S M., The Geochemical Evolution of the Continental Crust[J].Review of Geophys, 1995, 33(2):241-265.
[33]张旗，王焰，李承东，等.花岗岩的Sr-Yb 分类及其地质意义[J].岩石学报，2006，22（9）:2249-2269.
[34]Johannes W, Holtz F.Petrogenesis and Experimental Petrology of Granitic Rocks[M].Berlin:Spring-verlag,1996:1-335.
[35]赵宗溥.中朝准地台前寒武纪地壳演化[M].北京:科学出版社，1993:284-330.
[36]伍家善，耿元生，沈其韩.中朝古大陆太古宙地质特征及构造演化[M].北京:地质出版社，1998:192-211.
[37]Zhai M G,Liu W J.Palaeoproterozoic tectonic history of the North China craton: a review[J].precambrian research,2003,122(1-4):183-199.
[38]刘树文，吕勇军，王伟，等.冀北太古代花岗质片麻岩的成因[J].岩石学报，2011，27（4）:909-921.
[39]李江海，牛向龙，程素华，等.大陆克拉通早期构造演化历史探讨:以华北为例[J].地球科学，2006，（03）:5-13.
[40]万渝生，刘敦一，王世炎，等.登封地区早前寒武纪地壳演化—地球化学和锆石SHRIMP U-Pb年代学制约[J].地质学报，2009，83（7）:982-999.
[41]Harris N B W,Marzouki F M H,Ali S.The Jabel Sayd complex, Arabian shield: Geochemistry constrains on the origin of peralkaline and related granites[J].Journal of the geological society,1986,143(2):287-295.
[42]Pearce J A,Harris N B W,Tindle A G.Trace element discrimination diagrams for the tectonic interpretation of granitic rocks[J].Journal of petrology,1984,25(4):956-983.