2. 中国科学院矿物学与成矿学重点实验室, 广州 510640;
3. 西藏自治区地质矿产勘查开发局, 拉萨 850000
2. Key Laboratory of Mineralogy and Metallogeny, Chinese Academy of Sciences, Guangzhou 510640, China;
3. Geological Mining Exploration and Development Bureau, Tibet Autonomous Region, Lhasa 850000, China
0 引言
冈底斯—念青唐古拉地块(也称拉萨地块)位于亚洲大陆的最南缘，是青藏高原重要的组成部分，其记录了多期次碰撞造山运动，因此，一直以来备受国内外地学学者的关注和青睐。前人[1-6]研究结果表明，冈底斯成矿带中岩浆岩带相继经历了主碰撞(65~41 Ma)、晚碰撞(40~26 Ma)和后碰撞(25~0 Ma)演化过程，尤其以始新世主碰撞造山(52~41 Ma，主要集中于52~47 Ma[6])发生的岩浆作用最具特征，形成壳/幔混源的花岗岩-辉长岩组合，产于冈底斯弧花岗岩基内部，呈复式岩体断续分布于其南缘[3]。冈底斯斑岩铜钼成矿带位于冈底斯岩基的南部，是中国重要的铜矿资源富集区，已探明的驱龙、甲玛、厅宫、雄村和朱诺等超大型—大型斑岩铜矿床[7-9]铜矿资源量超过3 000万t。花岗质侵入岩的研究不仅对阐明冈底斯造山带的形成演化具有重要意义，而且对论证岩浆活动与地壳演化、成矿作用的关系均具有重要意义。
多仁则—桑阿卡地区发育大面积的侵入岩，以中酸性岩为主。前人开展亚模等4幅区域地质调查[10]时，在该地区厘定出不同的侵入岩填图单元，并做了一定量的岩石地球化学样品测试工作，为多仁则—桑阿卡地区侵入岩形成、演化提供了扎实的基础资料。但缺少岩浆作用与成矿作用的研究，不利于指导勘查找矿工作。作者依托“西藏昂仁—谢通门朱诺铜矿整装勘查区矿产调查与找矿预测”项目，对朱诺矿区南部的多仁则—桑阿卡地区进行了1:50 000专项地质填图，在该区的黑云花岗闪长岩体中发现8条铜矿化体和1条铜金多金属矿化体。在前人工作基础上，选择研究区内与成矿密切相关的侵入岩为研究对象，对研究区内的侵入岩进行详细的野外地质调查、典型样品采集，并开展了岩相学、锆石U-Pb年代学、岩石地球化学及Sr-Nd同位素研究，确定其岩石类型、形成时代、地球化学特征，探讨其岩石成因及地质意义，以期为朱诺整装勘查区区域成矿规律研究和区域矿产预测提供新资料，扩大找矿方向。
1 区域地质概况冈底斯成矿带呈近东西向展布于亚洲大陆南缘，夹持于班公湖—怒江缝合带与印度河—雅鲁藏布江缝合带之间[5]。多仁则—桑阿卡地区位于冈底斯—念青唐古拉复合岩浆弧(Ⅳ)之南冈底斯火山-岩浆弧(Ⅳ-1)中段南部[4](图 1a)，同时处于朱诺整装勘查区南部始新世火山-岩浆成岩巴热拉铅多金属成矿亚(次级)带。研究区出露地层为古新统林子宗群典中组(E1d)火山岩、始新统秋乌组(E2q)沉积岩及少量的第四系(Q)冲洪积物(图 1b)。地层产状总体倾向南西。
|
| 1.第四系；2.秋乌组二段；3.秋乌组一段；4.典中组；5.始新世黑云母二长花岗岩；6.始新世中粒二长花岗岩；7.始新世黑云母花岗闪长岩；8.始新世细粒花岗闪长岩；9.脉岩(花岗岩脉、石英脉、闪长玢岩脉、辉绿岩脉)；10.黑云母花岗闪长岩；11.断层；12.地质界线；13.岩相界线；14.角度不整合界线；15.铜(金)矿(化)体及代号；16.铜矿点、铜金矿化点；17.亚模地区区域地质调查项目的采样位置及样品编号；18.本次采集样品位置及编号；19.同位素年龄；20.地名。a底图据文献[11]修编。 图 1 多仁则—桑阿卡地区大地构造位置(a)及研究区地质图和采样位置(b) Fig. 1 Tectonic sketch map of Duorenze-Sangaka area (a) and geological map showing sampling position(b) of the study area |
|
|
研究区内岩浆活动强烈，岩石类型较多，分布在研究区北部，岩性主要为花岗闪长岩、黑云母花岗闪长岩、二长花岗岩、黑云母二长花岗岩。在研究区南部、西南部，岩体与林子宗群的典中组火山岩地层、始新统的秋乌组沉积岩地层呈侵入接触。岩体内部发育较多的晚期脉岩，如花岗岩脉、石英脉、闪长玢岩脉、辉绿岩脉。研究区内发育北西向、北东向两组断裂构造，控制着该区的火山-岩浆活动和热液蚀变及矿化。
研究区内及邻区始新世岩浆作用与铜多金属成矿作用具有十分密切的空间和成因关系。例如：巴热拉矿床辉钼矿Re-Os年龄为48.6 Ma(同吴松副教授口头交流)；本次调查中新发现两处矿床(点)，分别为曲嘎甫铜钼多金属矿床(距巴热拉矿区东边界2.0~2.5 km、距多仁则—桑阿卡铜多金属矿点西边界12~13 km，与成矿密切的细粒花岗岩形成时代为始新世)、多仁则—桑阿卡铜多金属矿点(与成矿密切的黑云母花岗闪长岩形成时代为始新世)；距朱诺斑岩铜矿区南约28 km的罗布真铅锌矿床矿区花岗闪长斑岩锆石U-Pb年龄为(52.7±1.8)Ma[12]，赋矿火山岩锆石U-Pb年龄为(49.14±0.86)Ma[13]。
2 矿区地质及矿化特征 2.1 矿区地质多仁则—桑阿卡铜多金属成矿区位于昂仁县秋窝乡北的多仁则—查嘎一带，主要出露林子宗群。其中：典中组火山岩地层分布在研究区西部，岩性为英安岩、英安质(晶屑)凝灰岩。秋乌组分为上下两段：下段岩性为紫红色复成分砾岩夹中—厚层状含砾不等粒岩屑长石砂岩和炭质页岩，局部形成煤线和薄煤层(如秋乌煤矿)；上段岩性为长石砂岩夹泥质粉砂岩，偶夹页岩。
矿区内岩浆岩发育。火山岩主要分布在林子宗群典中组中。侵入岩主要为始新世花岗闪长岩、黑云母花岗闪长岩、二长花岗岩、黑云母二长花岗岩；次为花岗闪长斑岩，岩体中可见到暗色微粒包体(MME)，暗示岩浆发生混染作用。其中，黑云母花岗闪长岩及晚期的岩脉(辉绿岩脉和闪长玢岩脉)构成含矿岩体，以岩基的形式产出，地表形态呈左躺卧向上开口的“喇叭状”，岩基长轴总体北西向展布，长7.0~8.0 km，短轴长1.6~4.0 km，面积约17 km2。矿区内主要发育北东向和北西向断裂，断裂构造夹持区矿化较好，是主要的赋矿地段。区内已知的矿(化)点有多仁则铜矿点、桑阿卡铜矿点及查嘎金铜多金属矿化点。因此，该地区具备很好的成矿地质条件和很大的找矿潜力。
2.2 矿化特征多仁则—桑阿卡地区铜多金属矿(化)体主要产于黑云母花岗闪长岩体中，受区域性北东向和北西向断裂构造所产生的一系列小裂隙、破碎蚀变带控制明显。本次工作共圈定出9个矿(化)体(图 1b)。各矿(化)体特征如下：
S-CuⅠ：该矿脉赋矿围岩为黑云母花岗闪长岩。地表出露长120~130 m，宽20~100 cm，脉内可见孔雀石化，未见原生金属矿化，发育硅化、钾化和绢英岩化等蚀变现象。Cu品位为0.31%~4.39%，平均品位为1.61%。
S-CuⅡ：该矿脉赋矿围岩为黑云花岗闪长岩，其中北东部矿化较南西部矿脉弱。地表出露长160~170 m，宽20~80 cm，脉内有孔雀石化，未见原生金属矿化，发育硅化、钾化等蚀变。Cu品位为0.50%~1.54%，平均品位为0.80%。
S-CuⅢ：该矿脉赋矿围岩为灰白色黑云母花岗闪长岩，该矿北西为一条构造破碎带，推测该构造与成矿关系密切。矿化好，目前控制规模较小，长125~130 m，宽2~3 m，矿脉中见孔雀石、黄铜矿、斑铜矿(图 2a，b)，蚀变主要为硅化和钾化。Cu品位为0.41%~0.61%，平均品位为0.53%。
|
| a.岩体露头中铜矿化；b.含黄铜矿、斑铜矿、孔雀石、蓝铜矿黑云母花岗闪长岩；c，d.采样露头岩性特征；e.黑云母花岗闪长岩手标本；f，g.黑云母花岗闪长岩镜下特征。Mal.孔雀石；Az.蓝铜矿；Cp.黄铜矿；Bn.斑铜矿；Qtz.石英；Pl.斜长石；Kfs.钾长石；Hbl.角闪石；Bt.黑云母；Gra.文象结构。 图 2 多仁则—桑阿卡地区黑云母花岗闪长岩体露头及岩相学特征 Fig. 2 Occurrence and petrographic characteristics of the biotite granodiorite in Duorenze-Sangaka area |
|
|
S-CuⅣ：矿化产于花岗闪长岩体中，规模不大，孔雀石化蚀变岩带长3~5 m，宽10~30 cm，蚀变为钾化、硅化和褐铁矿化等。Cu品位为0.71%。
S-CuⅤ：产于花岗闪长岩体中，规模不大，断续长3~4 m，宽20~50 cm，矿脉内见孔雀石化，发育钾化、硅化和褐铁矿化等蚀变。Cu品位为0.23%。
S-CuⅥ：产于断裂破碎蚀变带中，孔雀石矿化脉长6~7 m，宽20~30 cm，北东向展布，脉内可见孔雀石化，并发育硅化、钾化及黏土化等蚀变。Cu品位为0.18%。
S-CuⅦ：矿(化)体产于花岗闪长岩体中，规模不大，连续长20~30 m，宽20~100 cm，矿石矿物见孔雀石化，发育钾化、硅化和褐铁矿化等蚀变。Cu品位为1.97%。
S-CuⅧ：矿(化)体产于花岗闪长岩体中的蚀变带内，规模大，连续性好，矿化带走向40°~43°，倾向北西，延长270~300 m，宽3~6 m，矿脉内见孔雀石化，发育钾化、硅化和褐铁矿化等蚀变。Cu品位为0.62%~1.97%，平均品位为1.16%。
C-AuCuⅠ：该矿脉赋矿围岩为黑云花岗闪长岩，该矿化脉明显受近东西向的构造控制。规模较大，地表长500~550 m，宽2~3 m，脉内可见有孔雀石化，未见原生金属矿化，发育硅化、钾化和绢英岩化等蚀变。3件样品分析结果显示：Au品位为0.27×10-6~12.45×10-6，平均品位为4.44 ×10-6；Cu品位为0.15%~2.04%，平均品位为1.21%。其中1件样品Ag品位大于1 583×10-6。
3 样品采集及分析测试 3.1 取样位置及岩性特征本文所测样品(2018DRZ-1TW/1GS/1SN、2018DRZ-2GS/2SN、2018DRZ-3GS)均为中细粒黑云母花岗闪长岩，采自多仁则铜矿点南部，采样坐标位置分别为(87°25′45″ E、29°24′32″ N)(87°25′43″ E、29°24′30″ N)(87°25′42″ E、29°24′28″ N)，其余样品为亚模地区区域地质调查项目采集的样品。其中，2018DRZ-1TW做了年代学测试，2018DRZ-1SN和2018DRZ-2SN做了Sr-Nd同位素测试，3件样品均做了地球化学分析。样品位置见图 1b。
中细粒黑云母花岗闪长岩新鲜面为灰白色(图 2c，d，e)，具中细粒结构(图 2f，g)，块状构造。岩石由主要矿物斜长石、石英、钾长石及次要矿物黑云母、角闪石等组成。斜长石为半自形柱状，发育有聚片双晶，具轻微绢云母化，粒径为0.2~5.0 mm，体积分数约为50%；钾长石为半自形柱状、板状，主要为正长石，部分与石英呈有规律的镶嵌形成显微文象结构(图 2f，g)，粒径为0.2~2.0 mm，体积分数约20%；石英为他形粒状，粒径为0.2~2.0 mm，体积分数约为20%；黑云母为自形片状，部分不同程度发育绿泥石化，片径为0.2~1.0 mm，体积分数为5%~8%；角闪石为自形柱状，粒径为0.2~1.0 mm，体积分数为2%~5%。
3.2 分析测试主、微量和稀土元素测试在华北有色地质勘查局燕郊中心实验室完成。主量元素采用XRF仪器测定，仪器型号为AB-104L，微量、稀土元素采用ICP-MS等离子体质谱分析测试，准确度误差均小于5%。Sr-Nd同位素测试在核工业北京地质研究院分析测试研究中心完成，测试方法为同位素稀释法，测试仪器为ISOPROBE-T，实验过程中相对湿度50%，温度为20 ℃，同位素比值误差以2σ计。相关仪器运行条件及详细分析流程见文献[14]。
锆石分选在河北廊坊区域地质调查研究所实验室完成。首先利用无污染标准技术对锆石进行分选，并在双目显微镜下尽量挑选出无裂缝、无包体、透明干净的自形锆石颗粒，将其与标样一起置于环氧树脂中制靶，磨至一半，使锆石内部暴露；然后对锆石进行透射光和反射光显微照相以及阴极发光图像分析，观察锆石的内部结构，选择无裂纹、无包裹体的部位测试点位。LA-ICP-MS锆石U-Pb测年在中国地质大学(武汉)地质过程与矿产资源国家重点实验室完成，测试仪器型号为Agilent7700型ICP-MS，激光斑束直径为32 μm，剥蚀深度为20~40 μm。采用He作为剥蚀物质的载气，测试过程中采用每隔5个测点中插入1个标准样品点，测试样品之前和之后分别插入2个91500和GJ-1测点。运用中国地质大学(武汉)开发的ICPMSDataCal 8.3处理信号[15]，数据处理过程中采用91500标准锆石作为外标，29Si作为内标。样品谐和图解和加权平均年龄采用Isoplot 3.0完成[16]。
4 测试结果 4.1 LA-ICP-MS锆石U-Pb测年中细粒黑云母花岗闪长岩LA-ICP-MS锆石U-Pb测年结果见表 1。锆石为浅色或浅灰白色，长柱状或柱状，自形程度较好。颗粒大小差别较大，长为150 ~200 μm，宽为60 ~130 μm，长宽比为1:1~2:1。CL图像特征显示有清晰的内部结构和振荡环带(图 3)；锆石的Th和U质量分数分别为(96.81~804.47)×10-6和(110.68~325.79)×10-6，Th/U值为0.81~2.47(表 1)，具有典型岩浆锆石特征[17]。测点年龄投点落在谐和线上(图 4a)，表现为锆石群的特征，说明测点具有一致的207Pb/235U、206Pb/238U表面年龄和锆石在岩浆演化中没有发生明显的Pb丢失。20颗锆石20个数据点的206Pb/238U年龄在(49.8~47.6)Ma之间，206Pb/238U加权平均年龄为(49.0±0.7)Ma(MSWD=0.27，n=20)(图 4b)。因此，该年龄代表中细粒黑云母花岗闪长岩的结晶年龄，其形成时代为始新世。
| 测点号 | wB/10-6 | Th/U | 同位素比值 | 年龄/Ma | |||||||||||
| Pb | Th | U | 207Pb/206Pb | σ | 207Pb/235U | σ | 206Pb/238U | σ | 207Pb/235U | σ | 206Pb/238U | σ | |||
| 01 | 1.20 | 96.81 | 117.97 | 0.82 | 0.051 000 | 0.004 988 | 0.049 586 | 0.004 296 | 0.007 483 | 0.000 192 | 48.6 | 3.8 | 48.1 | 1.2 | |
| 02 | 1.74 | 147.65 | 164.77 | 0.90 | 0.050 271 | 0.006 992 | 0.051 107 | 0.006 658 | 0.007 430 | 0.000 272 | 51.6 | 3.9 | 47.7 | 1.7 | |
| 03 | 1.97 | 192.58 | 181.41 | 1.06 | 0.058 929 | 0.010 130 | 0.058 993 | 0.008 613 | 0.007 626 | 0.000 391 | 50.7 | 3.9 | 49.0 | 2.5 | |
| 04 | 2.79 | 250.25 | 268.94 | 0.93 | 0.049 461 | 0.006 697 | 0.050 442 | 0.005 562 | 0.007 759 | 0.000 267 | 49.1 | 4.2 | 49.8 | 1.7 | |
| 05 | 1.86 | 204.76 | 169.19 | 1.21 | 0.052 701 | 0.005 017 | 0.053 068 | 0.004 401 | 0.007 463 | 0.000 225 | 50.6 | 6.4 | 47.9 | 1.4 | |
| 06 | 1.82 | 198.93 | 174.65 | 1.14 | 0.051 053 | 0.003 838 | 0.052 697 | 0.003 261 | 0.007 709 | 0.000 192 | 58.2 | 8.3 | 49.5 | 1.2 | |
| 07 | 1.64 | 144.93 | 154.52 | 0.94 | 0.045 048 | 0.005 344 | 0.045 200 | 0.004 897 | 0.007 526 | 0.000 289 | 50.0 | 5.4 | 48.3 | 1.9 | |
| 08 | 1.90 | 160.32 | 181.35 | 0.88 | 0.051 147 | 0.005 162 | 0.048 915 | 0.004 265 | 0.007 546 | 0.000 221 | 52.5 | 4.2 | 48.5 | 1.4 | |
| 09 | 2.13 | 196.15 | 211.56 | 0.93 | 0.051 698 | 0.006 611 | 0.050 377 | 0.005 282 | 0.007 754 | 0.000 337 | 52.1 | 3.1 | 49.8 | 2.2 | |
| 10 | 2.13 | 181.52 | 222.75 | 0.81 | 0.052 345 | 0.005 348 | 0.054 786 | 0.004 941 | 0.007 728 | 0.000 197 | 44.9 | 4.8 | 49.6 | 1.3 | |
| 11 | 2.23 | 219.09 | 213.26 | 1.03 | 0.051 012 | 0.005 038 | 0.051 582 | 0.004 455 | 0.007 701 | 0.000 232 | 48.5 | 4.1 | 49.5 | 1.5 | |
| 12 | 1.80 | 175.09 | 173.79 | 1.01 | 0.056 157 | 0.005 980 | 0.054 550 | 0.004 640 | 0.007 628 | 0.000 232 | 49.9 | 5.1 | 49.0 | 1.5 | |
| 13 | 3.00 | 278.49 | 286.86 | 0.97 | 0.051 982 | 0.004 147 | 0.054 859 | 0.004 213 | 0.007 736 | 0.000 214 | 54.2 | 4.8 | 49.7 | 1.4 | |
| 14 | 2.56 | 235.47 | 250.43 | 0.94 | 0.057 533 | 0.005 699 | 0.055 127 | 0.004 281 | 0.007 742 | 0.000 254 | 51.1 | 4.3 | 49.7 | 1.6 | |
| 15 | 2.14 | 188.58 | 218.34 | 0.86 | 0.057 048 | 0.008 032 | 0.054 411 | 0.006 146 | 0.007 677 | 0.000 316 | 53.9 | 4.5 | 49.3 | 2.0 | |
| 16 | 4.48 | 804.47 | 325.79 | 2.47 | 0.057 540 | 0.007 134 | 0.053 544 | 0.005 610 | 0.007 661 | 0.000 256 | 54.2 | 4.1 | 49.2 | 1.6 | |
| 17 | 2.50 | 320.28 | 216.39 | 1.48 | 0.055 214 | 0.010 920 | 0.046 134 | 0.009 906 | 0.007 408 | 0.000 369 | 54.5 | 4.1 | 47.6 | 2.4 | |
| 18 | 2.38 | 217.91 | 221.58 | 0.98 | 0.051 000 | 0.004 988 | 0.049 586 | 0.004 296 | 0.007 483 | 0.000 192 | 53.8 | 5.9 | 48.1 | 1.2 | |
| 19 | 1.74 | 151.44 | 170.90 | 0.89 | 0.050 271 | 0.006 992 | 0.051 107 | 0.006 658 | 0.007 430 | 0.000 272 | 53.0 | 5.4 | 47.7 | 1.7 | |
| 20 | 1.12 | 133.84 | 110.68 | 1.21 | 0.058 929 | 0.010 130 | 0.058 993 | 0.008 613 | 0.007 626 | 0.000 391 | 45.8 | 9.6 | 49.0 | 2.5 | |
|
| 图 3 多仁则—桑阿卡地区黑云母花岗闪长岩体锆石CL图像 Fig. 3 Cathodoluminescence(CL) of representative zircon grains from the biotite granodiorite in Duorenze-Sangaka area |
|
|
|
| 图 4 多仁则—桑阿卡地区黑云母花岗闪长岩体锆石U-Pb谐和图(a)和加权平均年龄图(b) Fig. 4 Zircon U-Pb age concordia diagram(a) and weighted mean age diagram(b) of the biotite granodiorite in Duorenze-Sangaka area |
|
|
中细粒黑云母花岗闪长岩的主量元素、微量元素及稀土元素分析结果见表 2。
| 样品号 | SiO2 | TiO2 | Al2O3 | Fe2O3 | FeO | MnO | MgO | CaO | Na2O | K2O | P2O5 | CO2 | H2O- | ||
| 2018DRZ-1GS | 66.68 | 0.51 | 15.06 | 1.92 | 2.28 | 0.07 | 1.25 | 2.32 | 3.36 | 3.68 | 0.16 | 1.16 | — | ||
| 2018DRZ-2GS | 67.00 | 0.51 | 15.03 | 1.97 | 2.33 | 0.07 | 1.28 | 2.28 | 4.19 | 3.28 | 0.14 | 0.62 | — | ||
| 2018DRZ-3GS | 67.32 | 0.50 | 14.98 | 2.18 | 2.38 | 0.07 | 1.14 | 2.17 | 3.87 | 3.53 | 0.15 | 0.52 | — | ||
| D6048-1* | 67.39 | 0.53 | 15.52 | 2.08 | 1.76 | 0.08 | 1.12 | 3.04 | 3.88 | 3.59 | 0.16 | — | 0.24 | ||
| D6049-1* | 67.06 | 0.54 | 15.40 | 1.64 | 1.93 | 0.06 | 1.21 | 3.13 | 4.19 | 3.63 | 0.17 | — | 0.27 | ||
| D7143-1* | 68.04 | 0.45 | 15.04 | 2.00 | 1.75 | 0.07 | 1.21 | 2.88 | 3.50 | 3.74 | 0.15 | — | 0.22 | ||
| D7146-1* | 68.61 | 0.42 | 15.21 | 1.40 | 1.79 | 0.07 | 1.12 | 2.97 | 3.55 | 3.86 | 0.12 | — | 0.22 | ||
| D7148-1* | 66.35 | 0.54 | 15.54 | 1.95 | 2.10 | 0.08 | 1.43 | 3.48 | 3.46 | 3.86 | 0.15 | — | 0.21 | ||
| D6047-1* | 66.25 | 0.58 | 15.18 | 1.90 | 2.43 | 0.07 | 1.65 | 3.47 | 3.43 | 3.85 | 0.17 | — | 0.20 | ||
| D6050-1* | 67.15 | 0.53 | 15.23 | 2.30 | 1.76 | 0.06 | 1.12 | 1.83 | 4.02 | 3.90 | 0.16 | — | 0.37 | ||
| D6052-1* | 66.58 | 0.52 | 15.22 | 2.18 | 1.53 | 0.05 | 1.09 | 2.34 | 3.91 | 3.98 | 0.16 | — | 0.40 | ||
| Cγδ** | 62.17 | 0.60 | 16.54 | 2.17 | 1.94 | 0.11 | 2.54 | 4.51 | 3.21 | 1.72 | 0.18 | 0.44 | — | ||
| 样品号 | 烧失量 | 总和 | ID | A/NK | A/CNK | AR | R1 | R2 | A/MF | C/MF | Cr | Co | Ni | Cu | |
| 2018DRZ-1GS | 2.66 | 101.11 | 80.21 | 1.583 | 1.097 | 2.26 | 2 299 | 615 | 1.70 | 0.48 | 5.67 | 7.42 | 3.37 | 12.01 | |
| 2018DRZ-2GS | 1.85 | 100.55 | 80.23 | 1.439 | 1.030 | 2.52 | 2 108 | 610 | 1.66 | 0.46 | 5.16 | 4.66 | 3.35 | 7.55 | |
| 2018DRZ-3GS | 1.63 | 100.44 | 80.25 | 1.470 | 1.060 | 2.52 | 2 177 | 590 | 1.66 | 0.44 | 6.82 | 5.12 | 3.27 | 9.03 | |
| D6048-1* | 0.73 | 100.12 | 77.22 | 1.511 | 0.983 | 2.35 | 2 176 | 691 | 1.94 | 0.69 | 21.60 | 5.45 | 1.44 | 6.28 | |
| D6049-1* | 0.85 | 100.08 | 78.18 | 1.423 | 0.933 | 2.46 | 2 043 | 704 | 1.95 | 0.72 | 13.40 | 4.08 | 2.34 | 6.11 | |
| D7143-1* | 0.95 | 100.00 | 77.55 | 1.534 | 1.000 | 2.28 | 2 332 | 671 | 1.86 | 0.65 | 13.50 | 7.65 | 5.30 | 5.93 | |
| D7146-1* | 0.70 | 100.04 | 78.29 | 1.518 | 0.986 | 2.28 | 2 331 | 678 | 2.12 | 0.75 | 7.84 | 4.87 | 3.78 | 4.50 | |
| D7148-1* | 0.85 | 100.00 | 74.41 | 1.574 | 0.959 | 2.14 | 2 191 | 756 | 1.71 | 0.70 | 17.40 | 9.68 | 6.51 | 9.72 | |
| D6047-1* | 0.82 | 100.00 | 73.99 | 1.547 | 0.942 | 2.16 | 2 187 | 759 | 1.51 | 0.63 | 17.30 | 6.10 | 5.66 | 8.39 | |
| D6050-1* | 1.67 | 100.10 | 81.54 | 1.406 | 1.075 | 2.73 | 2 054 | 561 | 1.84 | 0.40 | 12.80 | 5.16 | 2.32 | 7.66 | |
| D6052-1* | 2.25 | 100.21 | 80.27 | 1.417 | 1.015 | 2.61 | 2 056 | 618 | 1.97 | 0.55 | 12.70 | 6.14 | 1.26 | 6.04 | |
| Cγδ** | — | 96.13 | 62.76 | 2.316 | 1.078 | 1.61 | 2 588 | 976 | 1.38 | 0.69 | 33.50 | 12.20 | 13.60 | 13.00 | |
| 样品号 | Zn | Mo | W | Sn | Ag | Au | Bi | Rb | Ba | Th | U | Nb | Ta | Pb | Sr |
| 2018DRZ-1GS | 42.70 | 0.61 | 2.87 | 2.72 | 0.06 | 2.02 | 0.21 | 122.00 | 679 | 16.90 | 2.57 | 12.10 | 0.97 | 19.50 | 345.00 |
| 2018DRZ-2GS | 39.60 | 0.55 | 1.63 | 2.59 | 0.06 | 3.32 | 0.55 | 82.50 | 660 | 15.20 | 2.59 | 13.10 | 0.88 | 19.70 | 441.00 |
| 2018DRZ-3GS | 38.10 | 0.85 | 1.35 | 3.15 | 0.05 | 1.23 | 0.14 | 95.50 | 651 | 16.90 | 2.81 | 13.10 | 1.01 | 17.00 | 421.00 |
| D6048-1* | 54.70 | 1.10 | 1.12 | 2.24 | 0.04 | 0.48 | 0.18 | — | 662 | 14.20 | 3.14 | 11.60 | 1.00 | 23.40 | 390.00 |
| D6049-1* | 42.30 | 0.78 | 1.22 | 1.61 | 0.07 | 0.47 | 0.12 | — | 599 | 14.10 | 2.89 | 12.00 | 1.08 | 21.00 | 349.00 |
| D7143-1* | 38.40 | 0.83 | 1.10 | 2.50 | 0.04 | 0.48 | 0.07 | — | 410 | 25.40 | 7.66 | 16.00 | 1.85 | 16.60 | 401.00 |
| D7146-1* | 29.20 | 0.83 | 0.60 | 3.05 | 0.03 | 0.76 | 0.07 | — | 354 | 24.40 | 9.45 | 23.20 | 2.33 | 16.40 | 357.00 |
| D7148-1* | 52.60 | 0.64 | 0.64 | 1.96 | 0.04 | 0.60 | 0.06 | — | 485 | 24.20 | 5.86 | 12.80 | 1.40 | 17.60 | 450.00 |
| D6047-1* | 39.90 | 0.98 | 1.17 | 2.40 | 0.04 | 0.71 | 0.12 | — | 484 | 21.10 | 6.46 | 12.00 | 1.40 | 19.60 | 430.00 |
| D6050-1* | 38.00 | 1.10 | 1.85 | 1.86 | 0.07 | 0.71 | 0.16 | — | 676 | 15.80 | 3.63 | 11.90 | 1.06 | 18.90 | 362.00 |
| D6052-1* | 39.80 | 0.87 | 1.60 | 1.64 | 0.07 | 0.42 | 0.10 | — | 580 | 13.40 | 2.86 | 10.40 | 0.88 | 23.80 | 302.00 |
| Cγδ** | 55.00 | 0.46 | 0.33 | 1.10 | 67.00 | 0.64 | 0.07 | 40.00 | 476 | 4.90 | 1.24 | 7.10 | 0.44 | 12.00 | 437.00 |
| 样品号 | Zr | Hf | Y | La | Ce | Pr | Nd | Sm | Eu | Gd | Tb | Dy | Ho | Er | |
| 2018DRZ-1GS | 172.00 | 8.56 | 15.70 | 30.40 | 63.90 | 7.42 | 26.10 | 5.06 | 1.19 | 4.23 | 0.69 | 3.70 | 0.68 | 2.08 | |
| 2018DRZ-2GS | 178.00 | 8.97 | 18.50 | 32.40 | 65.00 | 7.48 | 28.40 | 5.22 | 1.13 | 4.73 | 0.75 | 4.16 | 0.78 | 2.34 | |
| 2018DRZ-3GS | 185.00 | 9.22 | 19.80 | 35.60 | 70.20 | 8.33 | 31.60 | 5.80 | 1.22 | 5.05 | 0.76 | 4.55 | 0.80 | 2.46 | |
| D6048-1* | 179.00 | 5.40 | 32.60 | 37.90 | 73.00 | 7.95 | 28.60 | 5.20 | 1.32 | 4.87 | 0.81 | 4.26 | 0.95 | 3.16 | |
| D6049-1* | 168.00 | 5.21 | 32.80 | 35.30 | 65.70 | 7.32 | 26.80 | 5.00 | 1.14 | 4.68 | 0.80 | 4.19 | 0.95 | 3.17 | |
| D7143-1* | 142.00 | 4.68 | 36.10 | 103.00 | 162.00 | 14.20 | 44.30 | 6.90 | 1.17 | 6.64 | 0.98 | 4.84 | 1.05 | 3.52 | |
| D7146-1* | 130.00 | 4.84 | 47.80 | 31.10 | 79.60 | 9.82 | 35.00 | 7.10 | 1.08 | 6.35 | 1.11 | 5.89 | 1.36 | 4.68 | |
| D7148-1* | 158.00 | 4.60 | 24.80 | 33.60 | 69.20 | 7.59 | 26.50 | 4.70 | 1.07 | 4.20 | 0.67 | 3.39 | 0.75 | 2.53 | |
| D6047-1* | 154.00 | 5.00 | 26.00 | 33.00 | 66.50 | 7.23 | 25.60 | 4.53 | 1.04 | 4.15 | 0.65 | 3.33 | 0.74 | 2.49 | |
| D6050-1* | 160.00 | 5.20 | 34.80 | 37.00 | 72.20 | 7.98 | 28.40 | 5.24 | 1.22 | 4.94 | 0.85 | 4.50 | 1.01 | 3.34 | |
| D6052-1* | 161.00 | 4.92 | 27.50 | 33.20 | 64.80 | 7.41 | 26.60 | 4.89 | 1.16 | 4.39 | 0.73 | 3.64 | 0.82 | 2.70 | |
| Cγδ** | 112.00 | 4.80 | 14.10 | 26.00 | 50.00 | 6.37 | 26.60 | 5.30 | 1.13 | 2.90 | 0.62 | 3.70 | 0.75 | 2.30 | |
| 样品号 | Tm | Yb | Lu | ∑REE | LREE/HREE | (La/Yb)N | (La/Sm)N | (Gd/Yb)N | δEu | ||||||
| 2018DRZ-1GS | 0.36 | 2.12 | 0.32 | 537.33 | 9.46 | 10.29 | 3.88 | 1.11 | 0.76 | ||||||
| 2018DRZ-2GS | 0.38 | 2.39 | 0.34 | 565.43 | 8.81 | 9.72 | 4.01 | 1.15 | 0.68 | ||||||
| 2018DRZ-3GS | 0.39 | 2.35 | 0.36 | 614.32 | 9.13 | 10.87 | 3.96 | 1.21 | 0.67 | ||||||
| D6048-1* | 0.53 | 3.33 | 0.50 | 638.93 | 8.36 | 8.16 | 4.71 | 1.09 | 0.79 | ||||||
| D6049-1* | 0.53 | 3.31 | 0.49 | 599.22 | 7.80 | 7.65 | 4.56 | 1.06 | 0.71 | ||||||
| D7143-1* | 0.60 | 3.82 | 0.55 | 1 193.95 | 15.07 | 19.34 | 9.64 | 1.23 | 0.52 | ||||||
| D7146-1* | 0.86 | 5.76 | 0.87 | 742.56 | 6.09 | 3.87 | 2.83 | 1.04 | 0.48 | ||||||
| D7148-1* | 0.42 | 2.85 | 0.40 | 569.87 | 9.38 | 8.46 | 4.62 | 1.14 | 0.72 | ||||||
| D6047-1* | 0.44 | 2.84 | 0.41 | 555.27 | 9.16 | 8.33 | 4.70 | 1.16 | 0.72 | ||||||
| D6050-1* | 0.58 | 3.52 | 0.54 | 641.41 | 7.89 | 7.54 | 4.56 | 1.06 | 0.72 | ||||||
| D6052-1* | 0.44 | 2.85 | 0.41 | 569.06 | 8.64 | 8.36 | 4.38 | 1.09 | 0.75 | ||||||
| Cγδ** | 0.39 | 1.88 | 0.31 | 480.50 | 8.98 | 9.92 | 3.17 | 0.85 | 0.80 | ||||||
| 注：上角*为亚模等4幅区域地质调查报告[10]数据；上角**为中国花岗岩类化学元素丰度[18]数据；Cγδ为中国花岗闪长岩。ID.分异指数；R1=4Si-11(Na+K)-2(Fe+Ti)(原子数)；R2=6Ca+2Mg+Al(原子数)。主量元素质量分数单位为%；微量、稀土元素质量分数单位为10-6；Au元素质量分数单位为10-9。 | |||||||||||||||
岩石地球化学分析主量元素结果具有如下特征：w(SiO2)=66.25%~68.61%，均值为67.13%, w(K2O)=3.28%~3.98%，均值为3.72%，全碱(w(K2O+Na2O))=7.04%~7.92%，均值为7.48%，均高于中国花岗闪长岩(分别为62.17%、1.72%、4.93%[18])；在侵入岩w(Na2O+K2O)-w(SiO2)分类图(图 5a)上，样品点落到花岗闪长岩区域；碱度率指数(AR)为2.14~2.73(均值2.32)；在w(K2O)- w(SiO2)关系图(图 5b)上，样品均落在高钾钙碱性系列范围；w(Al2O3)=14.98%~15.54%，均值为15.22%，低于中国花岗闪长岩16.64%；铝饱和指数(A/CNK)为0.933~1.097(均值为1.013)，在A/NK-A/CNK图解(图 5c)上，部分样品点位于准铝质区域，部分位于过铝质区域，其中本次采集的3件样品出现较高质量分数的刚玉分子(2.25~4.48)，高于中国花岗闪长岩刚玉分子(1.7)，而收集样品的刚玉分子(0~1.49)低于中国花岗闪长岩[18]，表明岩体不同部位源区存在差异性；岩石具较高的分异指数(ID=73.99~81.54)；w(MgO)=1.09%~1.65%，均值为1.34%，小于3%。上述分析表明该岩体具有高钾钙碱性铝质高分异花岗岩地球化学特征。
稀土元素总量(w(∑REE))=537.33×10-6~1 193.95×10-6(均值657.00×10-6)，高于中国花岗闪长岩(480.50×10-6)[18]。在球粒陨石数据标准化的稀土配分模式图(图 6a)中[22]，呈明显的右倾型。(La/Yb)N为3.87~19.34，轻稀土较重稀土分馏明显，具有负铕异常(δEu=0.48~0.79，均值为0.70)。
|
| 图 6 多仁则—桑阿卡地区黑云母花岗闪长岩体稀土元素配分模式图(a)、微量元素蛛网图(b) Fig. 6 Chondrite-normalized REE distribution patterns diagram(a) and primitive mantle-normalized trace element spider diagram (b) of the biotite granodiorite in Duorenze-Sangaka area |
|
|
微量元素方面，在原始地幔数据标准化的微量元素蛛网图(图 6b)中[22]，岩石表现出富集强不相容元素Th、U、K、Nd、Zr、Hf和亏损高场强元素Nb、Ta、Sr、Ti、P等特征。
4.3 Sr-Nd同位素特征在壳-幔地球化学和深部地质作用过程时，Sr、Nd等同位素在示踪源区方面有突出的优点，但放射性成因子体同位素的积累程度与岩石的年龄有关[12, 23]。首先要对已测出的Sr、Nd同位素进行初始值计算，相关参数计算采用岩石结晶年龄(49.0 Ma)，测试结果及相关计算参数见表 3。
| 样品号 | 87Rb/86Sr | 87Sr/86Sr | 2 σ | (87Sr/86Sr)i | εSr(0) | εSr(t) | fRb/Sr | u/% |
| 2018DRZ-1SN | 1.023 373 | 0.705 988 | 0.000 018 | 0.705 280 | 21.1 | 11.9 | 11.37 | 60.41 |
| 2018DRZ-2SN | 0.541 388 | 0.705 902 | 0.000 016 | 0.705 530 | 19.9 | 15.4 | 5.55 | 62.09 |
| 样品号 | 147Sm/144Nd | 143Nd/144Nd | 2 σ | (143Nd/144Nd)i | εNd(t) | fSm/Nd | T2DM | INC |
| 2018DRZ-1SN | 0.113 712 | 0.512 501 | 0.000 012 | 0.512 465 | -2.2 | -0.4 | 1 034 | 0.46 |
| 2018DRZ-2SN | 0.107 808 | 0.512 526 | 0.000 008 | 0.512 492 | -1.6 | -0.5 | 992 | 0.43 |
| 注：为减少147Sm/144Nd变化对Nd模式年龄计算产生的影响，表中所列T2DM年龄统一采用二阶段模式计算，计算公式为T2DM=(1/ λSm)ln(1+ A)。式中：λSm=6.54×10-12a-1；A={(143Nd/144Nd)m-(143Nd/144Nd)DM-[(147Sm/144Nd)m-(147Sm/144Nd)c](e λt -1)}/[(147Sm/144Nd)c-(147Sm/144Nd)DM]，其中，(143Nd/144Nd)DM=0.513151，(147Sm/144Nd)c=0. 118，(147Sm/144Nd)DM=0.2136。εNd(t)=[(143Nd/144Nd)m /(143Nd/144Nd)CHUR-1]×104。式中，(143Nd/144Nd)CHUR=0.512630。以上式中，下标m代表样品现今测定值，下标DM代表亏损地幔值[24]，下标CHUR代表球粒陨石值[25]，下标c代表大陆壳平均值[26]。u为花岗岩中幔源成分百分比，u=(Ic- Ii)/(Ic- Im)×100%。式中：Ic为花岗岩形成时地壳的锶初始比值，Ic=0.7092-0.00000247× t；Ii为花岗岩样品的锶初始值；Im为花岗岩形成时上地幔的锶初始值(0.704)，Im =0.704-0.0000011× Ii× t([27]。INC为钕地壳指数[28-30]，代表岩石形成演化中地壳物质的贡献，INC=[ εNd(Rock)- εNd(M)]/ [ εNd(C)- εNd(M)]。式中：εNd(Rock)为测试样品值；εNd(M)为地幔值(-8.4)；εNd(C)为地壳值(-14.7)。 | ||||||||
Sr、Nd同位素具如下特征：1)(87Sr/86Sr)i分别为0.705 280、0.705 530，锶初始值小于大陆地壳范围(0.706~0.718)[12, 23]，通过邢凤鸣[27]花岗岩幔源成分百分比计算获知，幔源成分百分比(u)分别为60.41%、62.09%；2)(143Nd/144Nd)i分别为0.512 465、512 492，εNd(t)=-2.2~-1.6，fSm/Nd= -0.5~-0.4，均在-0.6~ -0.2之间，说明岩浆源区的Sm、Nd分馏不明显，故给出的模式年龄是有效的，可以参与有关讨论[31]，其对应的亏损地幔Nd的模式年龄T2DM=1 034~992 Ma(均值1 013 Ma)，通过D. J. Depaolo[28]钕地壳指数计算获得INC为0.43、0.46(均值为0.445)。以上Sr、Nd同位素特征表明岩浆源区为壳幔混源。
5 讨论 5.1 成岩时代大量研究[1-3]表明，与印度板块和欧亚板块碰撞相关的岩浆作用发生在早古近纪(始于65 Ma)。冈底斯成矿带中的岩浆岩带内岩浆作用具有两个明显的峰期，分别为早始新世(52~47 Ma，峰值约50 Ma)和早—中中新世(18~14 Ma，峰值约16 Ma)[32-33]。其中，早始新世形成的产于冈底斯南部(昂仁—谢通门—南木林—尼木—曲水)一带的大型花岗岩基是冈底斯岩浆作用的爆发期产物(52~47 Ma)，构成了冈底斯岩基的主体，这些岩体中含有大量的暗色微粒包体(MME)，暗示壳源花岗岩浆与幔源基性岩浆之间发生岩浆混合[5-6, 12, 33]。这与上述元素、同位素地球化学判别结果吻合，更能充分地说明岩浆源区具壳幔混源特征。
本次获得多仁则—桑阿卡地区含矿岩体成岩年龄为(49.0±0.7)Ma，形成时代为始新世，在冈底斯成矿带早始新世岩浆大爆发时期(52~47 Ma)范围内，是该时期的产物。与冈底斯成矿带上始新世成矿地质体成岩成矿年龄一致(表 4)，属于冈底斯成矿带早始新世大规模的铜多金属成岩成矿作用产物。
| 矿(化)区 | 坐标位置 | 成因类型 | 成矿时代 | 岩石/矿石测试对象 | 测试方法 | 年龄/Ma | 资料来源 |
| 沙让钼矿 | 92°41.5′E 30°11′N |
斑岩型 | 始新世 | 含矿花岗斑岩(锆石) 矿石(辉钼矿) |
U-Pb Re-Os等时线 |
53.0±1.0 51.0±1.0 |
[34] [35] |
| 热玛铁矿 | 90°59′E 29°53′N |
矽卡岩型 | 始新世 | 花岗质岩石(锆石) | MC-ICP-MS U-Pb | 49.7±0.1 | [36] |
| 吉如铜矿 | 88°53′E 29°39′N |
斑岩型 | 始新世 | 黑云母二长花岗岩(锆石) 矿石(辉钼矿) |
SHRIMP U-Pb Re-Os等时线 |
48.68±0.49 50.8~48.3 |
[37] |
| 加多捕勒铁铜矿 | 88°33.5′E 29°56′N |
矽卡岩型 | 始新世 | 黑云母二长花岗岩(锆石) | SHRIMP U-Pb | 50.9±1.8 | [38] |
| 郎达铜多金属矿 | 矽卡岩型 | 始新世 | 花岗质岩体(黑云母) | Ar-Ar等时线 | 47.6 | [37] | |
| 吉容铜多金属矿 | 斑岩-矽卡岩型 | 始新世 | 花岗闪长岩(锆石) | MC-ICP-MS U-Pb |
53.2±0.6 53.2±0.7 |
[39] | |
| 多仁则—桑阿卡铜多金属矿点 | 87°27′E 29°23′N |
斑岩型、热液型 | 始新世 | 黑云母花岗闪长岩(锆石) | LA-ICP-MS U-Pb | 49.0±0.7 | 本次 |
| 巴热拉铅锌多金属矿 | 87°13′E 29°27′N |
斑岩型、热液脉型 | 始新世 | 酸性侵入岩(锆石) 矿石(辉钼矿) |
U-PbRe-Os等时线 | 49.848.6 | 2019年6月同吴松口头交流 |
| 次玛班硕铜矿 | 87°21′E 29°35′N |
斑岩型 | 始新世中新世 | 二长花岗斑岩(锆石) 花岗闪长斑岩(锆石) |
LA-ICP-MS U-Pb | 50.0±0.8 15.2 |
赵亚云，待刊数据 [23] |
| 罗布真红山矿段 | 87°18′E 29°31′N |
斑岩型 | 始新世 | 花岗闪长斑岩(锆石) 英安岩(锆石) |
LA-ICP-MS U-Pb LA-ICP-MS U-Pb |
52.7±1.8 49.1±0.8 |
[12-13] |
| 曲水铜多金属矿化 | 始新世 | 花岗闪长岩(锆石) | SHRIMP U-Pb | 51.2±1.1 | [40] |
前人对冈底斯成矿带古新世—始新世岩浆起源研究较多，主要存在两种观点：Mo等[41]认为起源于碰撞过程中残留的新特提斯洋壳部分熔融；Lee等[42]认为地球化学特征的宽泛变化(如w(SiO2)=45%~80%)，源于壳幔多层次的相互作用。因此，有必要结合岩体特征、成岩年龄及构造背景对研究区含矿岩体进行深入解析，探讨岩石成因，有利于进一步丰富和加深对冈底斯成矿带内与壳幔花岗岩有关的Cu-Au-Mo-Fe-Pb-Zn成矿系统(52~47 Ma)[3]的认识。
研究区岩体具有高硅(w(SiO2)=67.13%)、高钾(w(K2O)=3.72%)、富碱(w(K2O+Na2O)=7.48%)、贫MgO(w(MgO)=1.34%，小于3%)的特征。在w(K2O)-w(SiO2)图解(图 5b)上，样品均落在高钾钙碱性系列范围。(La/Yb)N为3.87~19.34，轻稀土较重稀土分馏明显，具有负铕异常(δEu均值0.70，表 2、图 6a)，表明岩浆源区可能存在石榴石、角闪石的残留，岩浆发生斜长石的分离结晶作用，w(P2O5)与w(SiO2)呈负相关关系(图略)，表明多仁则—桑阿卡地区含矿岩体为Ⅰ型花岗岩[41]。
Th、U、K、Nd、Zr、Hf的富集和Nb、Ta、Sr、Ti、P的亏损，暗示岩浆部分熔融或结晶分异过程中具明显的斜长石分离[42]，但在La/Sm-w(La)图解(图 7a)中，具明显的部分熔融特征。MgO质量分数低(1.34%)，明显低于起源于俯冲洋壳的部分熔融，更可能来自于加厚下地壳的部分熔融[44]。在A/MF-C/MF图解(图 7b)中，除1个样品外，其余样品均落入到基性岩的部分熔融与变质砂岩的部分熔融交汇区域，但位于基性岩的部分熔融区域，表明物质源区更可能来源于基性下地壳熔融的贡献。
|
| a底图据文献[43]；b底图据文献[45]；c底图据文献[48]；d底图据文献[50]。MORB.大洋中脊玄武岩。 图 7 多仁则—桑阿卡地区黑云母花岗闪长岩体La/Sm-w(La)图(a)、A /MF- C/MF图(b)、εNd(t)-(87Sr/86Sr)i图(c)、w(Nb)-w(Y)构造环境判别图(d) Fig. 7 La/Sm vs. w(La) (a) and A/MF vs. C /MF (b) and εNd(t) vs. (87Sr/86Sr)i (c) and w(Nb) vs. w(Y) diagrams of the biotite granodiorite in Duorenze-Sangaka area |
|
|
Nb与Ta、Zr与Hf为强不相容元素，在蚀变和变质作用过程中比较稳定，因此，其可以示踪原始岩浆源区的特征[46]。样品的Nb/Ta值为8.57~14.82(均值为11.13)，接近大陆地壳平均值11[46]，Zr/Hf值为19.84~34.35(均值为28.29)，低于大陆地壳平均值(33)[33]，暗示岩浆源区主要来自地幔物质，但具明显壳幔混源特征；样品Sm/Nd值为0.16~0.20(均值为0.18)，高于上地壳(0.17)，表明岩浆源区相对较深[47]。
含矿岩体的同位素(87Sr/86Sr)i分别为0.705 280、0.705 530，小于大陆地壳范围(0.706~0.718)[14, 23]。通常情况下，εNd(t)值以0为界，正值代表岩浆来源于亏损地幔，负值代表源区来源于地壳或富集地幔[14]，本次计算获得含矿岩体的εNd(t)=-2.2~-1.6。在εNd(t) - (87Sr/86Sr)i图解(图 7c)中，样品投点分布在雅江MORB向下地壳演化趋势线和中上地壳演化之间，更靠近下地壳趋势演化线[48]，表明岩浆源区与下地壳更为密切相关。花岗岩的Nd模式年龄可以帮助判断其物质来源和岩石成因[49]。研究区含矿岩体((49.0±0.7)Ma)具有古老的Nd二阶段模式年龄(T2DM=1 013 Ma)，表明岩体岩浆形成后向上侵入的过程中遭受到具有中元古代古老地壳物质的混染。这一认识与花岗岩幔源成分比计算结果(61.25%)及钕地壳指数计算结果(44%)在误差范围内基本吻合。
综上各方面研究表明，多仁则—桑阿卡含矿岩体岩浆源区具壳幔混源特征。
5.3 成岩构造背景在w(Nb)-w(Y)构造环境判别图解(图 7d)中，本次研究样品落入火山弧花岗岩+同碰撞花岗岩区域，结合前人对冈底斯成矿带始新世花岗质岩基的研究[1-6, 10]和本次研究所得出的认识，岩浆形成于同碰撞构造环境。结合区域地质背景分析认为，岩浆源区与印度-欧亚大陆俯冲-碰撞后板片断离构造背景有关，岩浆起源于新生基性下地壳/地幔的部分熔融。其形成动力学机制是印度板块与欧亚板块碰撞过程中，印度大陆在57~53 Ma发生陡深俯冲[51]，使西藏大陆岩石圈加厚，深部板片发生断离，热的软流圈地幔物质通过板片断离而成的断离窗向壳幔界面上涌并注入，使岩石圈地幔熔融，继至诱发下地壳熔融。这些幔源岩浆在下地壳底部经历底侵和MASH(熔融、混染、储存和混合过程)，在应力松弛或地壳伸展背景下[3]，形成该地区乃至冈底斯成矿带壳幔混源岩浆-热液成矿作用。岩浆就位于不同的构造位置控制着不同类型的矿床(点)。
5.4 岩浆作用与成矿作用的关系现有资料[3, 33, 51-52]分析认为，多仁则—桑阿卡地区含矿岩体与铜金多金属成矿关系密切。冈底斯岩浆带的南部早始新世壳幔混源花岗岩，其形成普遍被认为与印度-欧亚板块碰撞后的板片断离有关。板片断离形成的断离窗是深部地幔流体上升的有利通道，岩浆-热液流体本身就载有成矿物质，并且在上升过程中将通道围岩中的成矿元素萃取出来，最后经岩浆混合作用形成具有更高背景值成矿(如W、Sn、Mo、Au)黑云母花岗闪长岩，为该地区成矿提供有利矿源和热源等条件。
结合冈底斯成矿带上早始新世成岩成矿事实(表 4)，如冈底斯东段北亚带沙让斑岩型钼矿床，冈底斯东段南亚带朗达矽卡岩型铜多金属矿床、吉容斑岩-矽卡岩型铜多金属矿床、热玛矽卡岩型铁矿床及曲水铜多金属矿化，冈底斯中段吉如斑岩型铜矿床、加多捕勒矽卡岩型铁铜矿床，冈底斯中段西南部的巴热拉热液脉型-斑岩型铅锌多金属矿床、多仁则—桑阿卡斑岩型-热液型铜多金属矿点、次玛班硕早期的斑岩型铜矿床及罗布真红山矿段的斑岩型铜多金属矿床，矿床成因类型主要为斑岩型和矽卡岩型，矿种以铜、铁、钼为主。因此，笔者认为多仁则—桑阿卡地区斑岩型-热液型铜多金属矿化是冈底斯成矿带中段中南部早始新世岩浆活动大爆发滞后的成矿响应，是与冈底斯成矿带内的壳幔花岗岩有关的Cu-Au-Mo-Fe-Pb-Zn成矿系统(52~47 Ma)的重要组成部分。上述认识为冈底斯成矿带上与始新世岩浆作用有成因联系的矿床研究提供了理论基础和资料，对朱诺整装勘查区矿产地质调查和找矿预测提供了很好的找矿方向，对下一步勘查工作部署和找矿突破具有重要现实意义。
6 结论1) 多仁则—桑阿卡地区含矿中细粒黑云母花岗闪长岩LA-ICP-MS锆石U-Pb年龄为(49.0±0.7)Ma，其形成时代为始新世。
2) 多仁则—桑阿卡地区含矿岩体具有高硅、高钾、富碱、贫MgO的特征，为高钾钙碱性系列的高分异Ⅰ型花岗岩；Eu具有负异常和Sr的亏损，暗示岩浆发生斜长石的分离结晶作用；微量元素表现出Th、U、K、Nd、Zr、Hf的富集和Nb、Ta、Sr、Ti、P的亏损特征。
3) 全岩Sr-Nd同位素特征、微量元素及元素比值揭示，岩浆源区为以地幔源区为主的壳幔混源，是在印度板块与亚洲大陆俯冲-碰撞后板片断离构造背景下，热的软流圈地幔物质通过板片断离窗向壳幔界面上涌并注入，使岩石圈地幔熔融，继而诱发下地壳部分熔融，在应力松弛或地壳伸展背景下，形成该地区乃至冈底斯成矿带壳幔混源岩浆-热液成矿作用。
4) 多仁则—桑阿卡地区斑岩型-热液型铜多金属矿化是冈底斯成矿带中段中南部早始新世岩浆活动大爆发滞后的成矿响应，是与冈底斯成矿带内的壳幔花岗岩有关的Cu-Au-Mo-Fe-Pb-Zn成矿系统(52~47 Ma)的重要组成部分。
| [1] |
莫宣学, 董国臣, 赵志丹, 等. 西藏冈底斯带花岗岩的时空分布特征及地壳生长演化[J]. 高校地质学报, 2005, 11(3): 281-290. Mo Xuanxue, Dong Guochen, Zhao Zhidan, et al. Spatial and Temporal Distribution and Characteristics of Granitoids in the Gangdese, Tibet and Implication for Crustal Growth and Evolution[J]. Geological Journal of China Universities, 2005, 11(3): 281-290. |
| [2] |
潘桂堂, 莫宣学, 侯增谦, 等. 冈底斯造山带的时空结构及演化[J]. 岩石学报, 2006, 22(3): 521-533. Pan Guitang, Mo Xuanxue, Hou Zengqian, et al. Spatial-Temporal Framework of the Gangdese Orogenic Belt and Its Evolution[J]. Acta Petrologica Sinica, 2006, 22(3): 521-533. |
| [3] |
侯增谦, 王二七, 莫宣学, 等. 青藏高原碰撞造山与成矿作用[M]. 北京: 地质出版社, 2008: 1-980. Hou Zengqian, Wang Erqi, Mo Xuanxue, et al. Collisional Orogeny and Metallogenesis of the Tibetan Plateau[M]. Beijing: Geological Publishing House, 2008: 1-980. |
| [4] |
纪伟强, 吴福元, 钟孙霖, 等. 西藏南部冈底斯岩基花岗岩时代与岩石成因[J]. 中国科学:D辑:地球科学, 2009, 39(7): 849-871. Ji Weiqiang, Wu Fuyuan, Zhong Sunlin, et al. Geochronology and Petrogenesis of Granitic Rocks in Gangdese Batholith, Southern Tibet[J]. Science China:Series D:Earth Science, 2009, 39(7): 849-871. |
| [5] |
孟元库, 许志琴, 高存山, 等. 藏南冈底斯带中段始新世岩浆作用的厘定及其大地构造意义[J]. 岩石学报, 2018, 34(3): 513-546. Meng Yuanku, Xu Zhiqin, Gao Cunshan, et al. The Identification of the Eocene Magmatism and Tectonic Significance in the Middle Gangdese Magmatic Belt, Southern Tibet[J]. Acta Petrologica Sinica, 2018, 34(3): 513-546. |
| [6] |
Mo X X, Dong G C, Zhao Z D, et al. Timing of Magma Mixing in the Gangdese Magmatic Belt During the India-Asian Collision:Zircon SHRIMP U-Pb Dating[J]. Acta Geologica Sinica, 2005, 79: 66-76. DOI:10.1111/j.1755-6724.2005.tb00868.x |
| [7] |
Hou Z Q, Cook N J. Metallogenesis of the Tibetan Collisional Orogen:A Review and Introduction to the Special Issue[J]. Ore Geology Reviews, 2009, 36(1/2/3): 2-24. |
| [8] |
Wang R, Richards J P, Hou Z Q, et al. Extent of Underthrusting of the Indian Plate Beneath Tibet Controlled of Miocene Porphyry Cu-Mo±Au Deposits[J]. Mineralium Deposita, 2014, 49: 165-173. DOI:10.1007/s00126-013-0507-y |
| [9] |
郑有业, 多吉, 张刚阳, 等. 西藏吉如斑岩铜矿床的发现过程及意义[J]. 矿床地质, 2007, 26(3): 317-321. Zheng Youye, Duo Ji, Zhang Gangyang, et al. Discovery of Jiru Porphyry Copper Deposit in Tibet and Its Significance[J]. Mineral Deposits, 2007, 26(3): 317-321. |
| [10] |
西藏自治区地质矿产勘查开发局区域地质调查大队.西藏亚模地区1: 50000区域地质调查报告[R].北京: 中国地质调查局, 2013: 1-299. Regional Geological Survey Team, Tibet Autonomous Region Geological Mining Exploration and Development Bureau.1: 50000 Regional Geological Survey Report of the Yamo Area, Tibet[R]. Beijing: China Geological Survey, 2013: 1-299. |
| [11] |
赵亚云, 刘晓峰, 刘远超, 等. 西藏朱诺矿区外围次玛班硕地区铜成矿有利条件分析[J]. 甘肃地质, 2017, 26(4): 28-36. Zhao Yayun, Liu Xiaofeng, Liu Yuanchao, et al. Copper Metallogenic Condition of Cimabanshuo Area Around Zhunuo Copper Mine in Tibet[J]. Gansu Geology, 2017, 26(4): 28-36. |
| [12] |
赵亚云, 刘晓峰, 刘波, 等. 西藏罗布真矿区花岗闪长斑岩锆石U-Pb年龄、地球化学特征[J]. 西藏地质, 2018, 45(2): 15-25. Zhao Yayun, Liu Xiaofeng, Liu Bo, et al. Geochemical Characteristics, Zircon U-Pb Age of the Granodiorite Porphyry in Luobuzhen Orefield, Tibet[J]. Tibet Geology, 2018, 45(2): 15-25. |
| [13] |
赵亚云, 杨春四, 吕金梁, 等. 西藏罗布真矿区林子宗群火山岩锆石U-Pb年龄、地球化学特征及意义[J]. 现代地质, 2019, 33(1): 73-85. Zhao Yayun, Yang Chunsi, Lü Jinliang, et al. Its Significance and Geochemical Characteristics, Zircon U-Pb Age of the Linzizong Group Volcanic Rocks in Luobuzhen Orefield, Tibet[J]. Geoscience, 2019, 33(1): 73-85. |
| [14] |
赵亚云, 张树明, 汤琳, 等. 龙首山中段芨岭花岗岩体Sr-Nd-Pb同位素特征及意义[J]. 地球科学, 2016, 41(6): 1016-1031. Zhao Yayun, Zhang Shuming, Tang Lin, et al. Sr-Nd-Pb Isotopic Characteristics and Its Geological Significance of the Jiling Grantic Pluton in the Middle Longshou Mountains Areas[J]. Earth Science, 2016, 41(6): 1016-1031. |
| [15] |
Liu Y S, Hu Z C, Gao S, et al. In Situ Analysis of Major and Trace Elements of Anhydrous Minerals by LA-ICP-MS Without Applying an Internal Standard[J]. Chemical Geology, 2008, 257: 34-43. DOI:10.1016/j.chemgeo.2008.08.004 |
| [16] |
Ludwig K R. User's Manual for Isoplot 3.0: A Geochronological Toolkit for Microsoft Excel[J]. Berkeley Geochronology Center, Special Publication, 2003, 4: 1-71.
|
| [17] |
姚晓峰, 唐菊兴, 李志军, 等. 班怒带西段尕尔穷铜金矿两套侵入岩源区及其地质意义:来自Hf同位素特征的指示[J]. 吉林大学学报(地球科学版), 2012, 42(增刊2): 188-198. Yao Xiaofeng, Tang Juxing, Li Zhijun, et al. Magma Origin of Two Plutons from Gaerqiong Copper-Gold Deposit and It's Geological Significance, Western Bangonghu-Nujiang Metallogenic Belt, Tibet:Implication from Hf Isotope Characteristics[J]. Journal of Jilin University(Earth Science Edition), 2012, 42(Sup.2): 188-198. |
| [18] |
史长义, 鄢明才, 迟清华, 等. 中国花岗岩类化学元素丰度[M]. 北京: 地质出版社, 2008: 1-124. Shi Changyi, Yan Mingcai, Chi Qinghua, et al. China Granite Chemical Element Abundance[M]. Beijing: Geological Publishing House, 2008: 1-124. |
| [19] |
Le Maitre R W. Igneous Rocks:A Classification and Glossary of Terms:Recommendations of the International Union of Geological Sciences, Sub-Commission on the Systematics of Igneous Rocks[M]. Cambridge: Cambridge University Press, 2002.
|
| [20] |
Peccerillo A, Taylor S R. Geochemistry of Eocene Calc-Alkaline Volcanic Rocks from the Kastamonu Area, Northern Turkey[J]. Contributions to Mineralogy and Petrology, 1976, 58(1): 63-81. DOI:10.1007/BF00384745 |
| [21] |
Maniar P D, Piccoli P M. Tectonic Discrimination of Granitoids[J]. Geological Society of America Bulletin, 1989, 101(5): 635-643. DOI:10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2 |
| [22] |
Sun S S, McDonough W F. Chemical and Isotopic Systematics of Oceanic Basalts: Implications for Mantle Composition and Processes[C]//Saunders A D, Norry M J. Magmatism in the Ocean Basins.London: Geological Society Special Publication, 1989: 313-345.
|
| [23] |
赵亚云, 刘晓峰, 刘远超, 等. 西藏次玛班硕地区由秋米斑岩体锆石U-Pb年龄、地球化学特征[J]. 地球科学, 2018, 43(12): 4551-4565. Zhao Yayun, Liu Xiaofeng, Liu Yuanchao, et al. Zircon U-Pb Ages and Geochemical Characteristics of Youqiumi Porphyry Pluton in Cimabanshuo Area, in Tibet[J]. Earth Science, 2018, 43(12): 4551-4565. |
| [24] |
Miller R G, O'Nions R K. Source of Precambrian Chemical and Clastic Sediments[J]. Nature, 1985, 314: 325-330. DOI:10.1038/314325a0 |
| [25] |
Bouvier A, Vervoort J D, Patchett P J. Lu-Hf and Sm-Nd Istope Composition of CHUR:Constraints from Unequilibrated Chlndrites and Implications for the Bulk Composition of Terrestrial Planets[J]. Earth and Plannetary Science Letters, 2008, 273(1/2): 48-57. |
| [26] |
Jahn B M, Condie K C. Volution of the Kaapvaal Craton as Viewed from Geochemical and Sm-Nd Isotopic Analyses of Intracratonic Pelites[J]. Geochimica et Cosmochimica Acta, 1995, 59(11): 2239-2258. DOI:10.1016/0016-7037(95)00103-7 |
| [27] |
邢凤鸣. 锶同位素初始比值在划分花岗岩成因类型上的应用探讨[J]. 岩石学报, 1987, 2(5): 75-80. Xing Fengming. On the Application of Initial Strontium Isotope Rations in Division of Granitic Types[J]. Acta Petrologica Sinica, 1987, 2(5): 75-80. |
| [28] |
Depaolo D J. Isotope Studies of Processes in Mafic Magma Chambers:Ⅰ:Kiglapait Intrusion, Labrador[J]. Journal of Petrology, 1985, 26: 925-951. DOI:10.1093/petrology/26.4.925 |
| [29] |
Depaolo D J, Perry F V, Baldridge W S. Crustal Versus Mantle Sources of Granitic Magmas:A Two-Parameter Model Based on Nd Isotopic Studies[J]. Royal Society of Edinburgh Transactions:Earth Sciences, 1992, 83: 439-446. DOI:10.1017/S0263593300008117 |
| [30] |
Mo X X, Hou Z Q, Niu Y L, et al. Mantle Contributions to Crustal Thickening During Continental Collision Evidence from Cenozoic Igneous Rocks in Southern Tibet[J]. Lithos, 2007, 96(1/2): 225-242. |
| [31] |
吴福元, 林强, 江博明. 造山带造山后花岗岩的同位素特点与地壳生长意义[J]. 科学通报, 1997, 42(20): 2188-2192. Wu Fuyuan, Lin Qiang, Jahn Boming. Northern China Orogenic Belt Post-Tectonic of Granitic Isotope Characteristics and Crustal Growth's Significance[J]. Chinese Science Bulletin, 1997, 42(20): 2188-2192. |
| [32] |
Zhu D C, Wang Q, Cawood P A, et al. Raising the Gangdese Mountains in Southern Tibet[J]. Journal of Geoophysical Research (Solid Earth), 2017, 122(1): 214-223. DOI:10.1002/2016JB013508 |
| [33] |
张泽明, 丁慧霞, 董昕, 等. 喜马拉雅造山带两种不同类型榴辉岩与印度大陆差异性俯冲[J]. 地球科学, 2019, 44(5): 1602-1620. Zhang Zeming, Ding Huixia, Dong Xin, et al. Two Contrasting Eclogite Types in the Himalayan Orogen and Differential Subduction of Indian Continent[J]. Earth Science, 2019, 44(5): 1602-1620. |
| [34] |
秦克章, 李光明, 赵俊兴, 等. 西藏首例独立钼矿:冈底斯沙让大型斑岩钼矿的发现及其意义[J]. 中国地质, 2008, 35(6): 1101-1112. Qin Kezhang, Li Guangming, Zhao Junxing, et al. Discovery of Sharang Large-Scale Porphyry Molybdenum Deposit, the First Single Mo deposit in Tibet and Its Significance[J]. Geology in China, 2008, 35(6): 1101-1112. |
| [35] |
唐菊兴, 陈毓川, 王登红, 等. 西藏工布江达县沙让斑岩钼矿床辉钼矿铼-锇同位素年龄及其地质意义[J]. 地质学报, 2009, 83(5): 698-704. Tang Juxing, Chen Yuchuan, Wang Denghong, et al. Re-Os Dating of Molybdenite from the Sharang Porphyry Molybdenum Deposit in Gongbo'gyamda County, Tibet and Its Geological Significance[J]. Acta Geologica Sinica, 2009, 83(5): 698-704. |
| [36] |
Zheng Y C, Fu Q, Hou Z Q, et al. Metallogeny of the Northeastern Gangdese Pb-Zn-Ag-Fe-Mo-W Polymetallic Belt in the Lhasa Terrane, Southern Tibet[J]. Ore Geology Reviews, 2015, 70: 510-532. DOI:10.1016/j.oregeorev.2015.04.004 |
| [37] |
张刚阳, 郑有业, 龚福志, 等. 西藏吉如斑岩铜矿:与陆陆碰撞过程相关的斑岩成岩成矿时代约束[J]. 岩石学报, 2008, 24(3): 473-479. Zhang Gangyang, Zheng Youye, Gong Fuzhi, et al. Geochronologic Constraints on Magmatic Intrusions and Mineralization of the Jiru Porphyry Copper Deposit, Tibet, Associated with Continent-Continent Collisional Process[J]. Acta Petrologica Sinica, 2008, 24(3): 473-479. |
| [38] |
于玉帅, 杨竹森, 多吉, 等. 西藏加多捕勒铁铜矿成矿岩体时代与成因:锆石U-Pb年龄、Hf同位素与稀土元素证据[J]. 矿床地质, 2011, 30(3): 420-435. Yu Yushuai, Yang Zhusen, Duo Ji, et al. Age and Petrogenesis of Magmatic Rocks from Jiaduobule Skarn Fe-Cu Deposit in Tibet:Evidence from Zircon SHRIMP U-Pb Dating, Hf Isotope and REE[J]. Mineral Deposits, 2011, 30(3): 420-435. |
| [39] |
Jiang Z Q, Wang Q, Derek A, et al. Zircon U-Pb Geochronology and Geochemistry of Late Cretaceous-Early Eocene Granodiorites in the Southern Gangdese Batholith of Tibet:Petrogenesis and Implications for Geodynamics and Cu±Au±Mo Mineralization[J]. International Geology Review, 2015, 57(3): 373-392. DOI:10.1080/00206814.2015.1009503 |
| [40] |
Mo X X, Dong G C, Zhao Z D, et al. Mantle Input to the Crust in Southern Gangdese, Tibet, During the Cenozoic:Zircon Hf Isotopic Evidence[J]. Journal of Earth Science, 2009, 20(2): 241-249. |
| [41] |
Mo X X, Niu Y L, Dong G C, et al. Contribution of Syncollisional Felsic Magmatism to Continental Crust Growth:A Case Study of the Paleogene Linzizong Volcanic Succession in Southern Tibet[J]. Chemical Geology, 2008, 250(1/2/3/4): 49-67. |
| [42] |
Lee H Y, Chung S L, Lo C H. Eocene Neotethyan Slab Break off in Southern Tibet Inferred from the Linzizong Volcanic Record[J]. Tectonophysics, 2009, 477: 20-37. DOI:10.1016/j.tecto.2009.02.031 |
| [43] |
赵亚云.龙首山中段古生代花岗岩岩石学、地球化学特征及地质意义[D].南昌: 东华理工大学, 2016: 1-79. Zhao Yayun. Petrology, Geochemistry Characteristics and Geological Significance of Paleozoic Granites in Middle Longshou Moutains, Gansu Province[D]. Nanchang: East China University of Technology, 2016: 1-79. http://cdmd.cnki.com.cn/Article/CDMD-10405-1016758454.htm |
| [44] |
陈希节, 许志琴, 孟元库, 等. 冈底斯带中段中新世埃达克质岩浆作用的年代学、地球化学及Sr-Nd-Hf同位素制约[J]. 岩石学报, 2014, 30(8): 2253-2268. Chen Xijie, Xu Zhiqin, Meng Yuanku, et al. Petrogenesis of Miocene Adakitic Diorite-Porphyrite in Middle Gangdese Batholith, Southern Tibet:Constraints from Geochemistry, Geochronology and Sr-Nd-Hf Isotopes[J]. Acta Petrologica Sinica, 2014, 30(8): 2253-2268. |
| [45] |
Alther R, Holl A, Hegner E. High-Potassium, Calc-Alkaline Ⅰ-Type Plutonism in the European Variscides:Northern Vosges(France) and Northern Schwarzwald(Germany)[J]. Lithos, 2000, 50: 51-73. DOI:10.1016/S0024-4937(99)00052-3 |
| [46] |
郭春丽, 曾令森, 高利娥, 等. 福建河田高分异花岗岩的矿物和全岩地球化学找矿标志研究[J]. 地质学报, 2017, 91(8): 1796-1817. Guo Chunli, Zeng Lingsen, Gao Li'e, et al. Highly Fractionated Grantitic Minerals and Whole-Rock Geochemistry Prospecting Markers in Hetian, Fujian Province[J]. Acta Geologica Sinica, 2017, 91(8): 1796-1817. |
| [47] |
李奋其, 李益多, 张士贞, 等. 西藏朗县地区增生楔杂岩带90 Ma岛弧型深成岩浆活动和意义[J]. 中国地质, 2016, 43(1): 142-152. Li Fenqi, Li Yiduo, Zhang Shizhen, et al. The 90 Ma Island-Arc Type Plutonism in the Subduction-Accretionary Complex in Langxian County Area, Tibet[J]. Geology in China, 2016, 43(1): 142-152. |
| [48] |
Hou Z Q, Yang Z M, Lu Y J, et al. A Genetic Linkage Between Subduction and Collision-Related Porphyry Cu Deposits in Continental Collision Zones[J]. Geology, 2015, 43(3): 247-250. |
| [49] |
凌洪飞, 沈渭洲, 孙涛, 等. 广东省22个燕山期花岗岩的源区特征及成因:元素及Nd-Sr同位素研究[J]. 岩石学报, 2006, 22(12): 2687-2703. Ling Hongfei, Shen Weizhou, Sun Tao, et al. Genesis and Source Characteristics of 22 Yanshanian Granites in Guangdong Province:Study of Element and Nd-Sr Isotopes[J]. Acta Petrologica Sinica, 2006, 22(12): 2687-2703. |
| [50] |
Pearce J A, Harris N B W, Tindle A G. Trace Element Discrimination Diagram for the Tectonic Interpretation of Granitic Rocks[J]. Journal of Petrology, 1984, 25(4): 956-983. DOI:10.1093/petrology/25.4.956 |
| [51] |
Zhu D C, Wang Q, Zhao Z D. Magmatic Record of India-Asia Collision[R].Beijing: Scientific Reports, 2015: 14289.
|
| [52] |
Leech M L, Singh S, Jain A K, et al. The Onset of India-Asia Continental Collision:Early, Steep Subduction Required by the Timing of HUP Metamorphism in the Western Himalaya[J]. Earth and Planetary Science Letters, 2005, 234: 83-97. DOI:10.1016/j.epsl.2005.02.038 |