岩石学报  2017, Vol. 33 Issue (2): 476-494   PDF    
西藏多龙矿集区尕尔勤斑岩铜矿床年代学及地球化学——兼论硅帽的识别与可能的浅成低温热液矿床
张志1, 方向1, 唐菊兴2, 王勤1, 杨超2, 王艺云1, 丁帅1, 杨欢欢1     
1. 中国地质调查局成都地质调查中心, 成都 610081;
2. 中国地质科学院矿产资源研究所, 国土资源部成矿作用与资源评价重点实验室, 北京 100037
摘要: 尕尔勤斑岩铜矿床是多龙超大型铜金矿集区内具有较大找矿潜力的矿床之一。本文以尕尔勤斑岩铜矿区内成矿花岗闪长斑岩及新发现的硅帽为研究对象,进行了矿床年代学、地球化学及找矿方向研究工作。尕尔勤矿床花岗闪长斑岩SiO2含量在60.01%~62.81%之间,K2O含量在1.86%~2.16%之间,Al2O3含量在15.12%~16.51%之间,属于钙碱性岩石系列。斑岩体锆石U-Pb年龄为124.4±0.4Ma(MSWD=0.42,n=25),属于早白垩世侵位;锆石εHft)值主要集中于1.15~9.71之间,tDMc模式年龄主要集中于559~1105Ma之间,表明岩浆岩区起源于具有幔源特征的深成熔体,而个别负εHft)值(-18.65,-19.75)的出现表明其受到了古老成熟地壳的轻微混染。斑岩体稀土元素显示为重稀土亏损的右倾分布型式,在原始地幔标准化图解中微量元素显示为Rb、Th、U、Sr等大离子亲石元素相对富集,而Nb、Ta、Ti、P、Zr等高场强元素相对亏损特征,体现出岛弧岩浆特有的地球化学特征,综合研究表明矿床形成与班公湖-怒江洋北向俯冲密切相关。硅帽中角砾状样品与层纹状(致密块状)样品地球化学特征具有明显区别,发育较少的层纹状硅帽样品可能为热水沉积与生物化学沉积共同作用形成的硅质岩,而大量存在的角砾状硅帽体现出明显热液充填交代成因特点,经地球化学剖面测量显示硅帽展布区域Cu、Au、Ag、Pb、Zn、As、Sb等成矿指示元素异常明显,结合高光谱矿物识别表明该区域存在浅成低温热液型矿床的巨大可能性,具有重要的勘查指示意义。
关键词: 斑岩铜矿床     浅成低温热液矿床     年代学     地球化学     尕尔勤     多龙矿集区     西藏    
Chronology, geochemical characteristics of the Gaerqin porphyry copper deposit in the Duolong ore concentration area in Tibet and discussion about the identification of the lithoscaps and the possible epithermal deposit
ZHANG Zhi1, FANG Xiang1, TANG JuXing2, WANG Qin1, YANG Chao2, WANG YiYun1, DING Shuai1, YANG HuanHuan1     
1. Chengdu Institute of Geology and Mineral Resources, Chinese Geological Survey, Chengdu 610081, China;
2. MLR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, CAGS, Beijing 100037, China
Abstract: Gaerqin porphyry copper deposit is one of the most prospecting deposits in Duolong large copper-gold ore concentration area. This paper has a research for the chronology, geochemistry and prospecting direction of the deposit based on the ore-bearing porphyry-granodiorite porphyry and new discovered lithoscaps. The result shows that granodiorite porphyry (SiO2 of 60.01%~62.81%, K2O of 1.86%~2.16%, Al2O3 of 15.12%~16.51%) belong to the calc-alkaline series. Zircon U-Pb age of granodiorite porphyry is 124.4±0.4Ma (MSWD=0.42, n=25), which suggest Early Cretaceous emplacement. Zircons from ore-forming porphries main show positive εHf(t) values of 1.15~9.71 and younger tDMc model ages of 559~1105Ma, while few negative εHf(t) values (-18.65, -19.75) have been found, it indicates that the granodiorite porphyry derived from mantle-derived magmas with minor ancient crustal contamination. The granodiorite porphyry show a right-sloping Rare earth elements (REE) pattern, and enrichment of large-ion lithophile elements (LILE:Rb, Th, U, Sr), depletion of high field strength elements (HFSE:Nb, Ta, Ti, P, Zr), consistent with geochemical feature of arc magmas. Comprehensive studies show formation of the deposit is closely relate to northward subduction of the Bangong-Nujiang ocean. The geochemical characteristics of brecciated lithoscaps samples show obvious hydrothermal replacement and are different from that of lamellar lithoscaps (dense block) samples which may be the products of the bio-genesis and hydrothermal sedimentation. Mineralization indicator elements Cu, Au, Ag, Pb, Zn, As, Sb are anomaly obvious according to geochemical profile, suggesting enormous possibilities of epithermal deposits in this deposit and important exploration significance.
Key words: Porphyry copper deposit     Epithermal Epithermal deposit     Chronology     Geochemistry     Gaerqin     Duolong ore concentration area     Tibet    
1 引言

青藏高原是世界的屋脊,其不仅发育有与经典洋壳俯冲模式吻合的成岩成矿作用 (Sillitoe, 1972; Zhu et al., 2009, 2010; Ji et al., 2009; 曲晓明等, 2009; 唐菊兴等, 2010; Li et al., 2011a, 2013),也发育着一系列与大陆碰撞作用相关的成矿事件 (侯增谦等, 2004; Hou et al., 2003, 2005, 2009; Qu et al., 2007, 2009; Yang et al., 2009; Li et al., 2006, 2011a; 陈衍景, 2013),是陆-陆碰撞造山动力学过程及其成矿作用探讨的天然实验室 (Yin and Harrison, 2000; Xu et al., 2012)。横贯青藏高原中部的班公湖-怒江缝合带作为冈瓦纳大陆的北界 (潘桂棠等, 2004; Yin and Harrison, 2000; Zhu et al., 2011, 2013),是冈底斯陆块与羌塘陆块的分界线 (任纪舜和肖黎薇, 2004),是西藏继玉龙、冈底斯成矿带后的第三条斑岩-矽卡岩铜矿带 (秦克章等, 2006; Li et al., 2011b, 2013)。该带虽找矿潜力巨大,但研究程度却十分薄弱,其从洋盆开启至俯冲与碰撞的构造演化过程也一直存在争论 (Pearce and Deng, 1988; Kapp et al., 2003; 李金祥等, 2008; 李光明等, 2011; 耿全如等, 2011; 江军华等, 2011; 曲晓明等, 2012, 2013; Li et al., 2014),制约着找矿工作的深入进行。作为缝合带构造演化研究的反演指示针与终点站,更多的成岩-成矿作用研究势在必行,特别是对于大型-超大型矿床乃至矿集区成岩成矿作用及其构造背景研究对于深入开展缝合带区域找矿工作来说尤为重要。

多龙矿集区是近年来藏北取得重大找矿突破的超大型铜金矿集区,已陆续发现了一批大型-超大型斑岩、斑岩-浅成低温热液型矿床 (李光明等, 2007; 李金祥等, 2008; 李玉彬等, 2012; Li et al., 2012a, b; 唐菊兴等, 2014),具有重大的科研与经济价值。截止目前为止,矿集区内已勘查出大型斑岩铜金矿床3个 (多不杂、波龙、拿若)(铜资源量>700万吨),超大型斑岩-浅成低温热液铜金矿床1个 (铁格隆南)(铜资源量>800万吨),此外还有一系列有待进一步进行勘查工作的矿床 (点)(图 1)。尕尔勤斑岩铜金矿位于多龙矿集区东南部,是矿集区内最具找矿潜力的矿床之一,虽已陆续进行了一些探矿工作,但至今未能取得大的找矿突破,且研究程度十分薄弱,对该矿床成岩成矿作用的研究对认识多龙矿集区的成因和进一步勘查具有重要的意义。因此,本文在详细野外地质调查的基础上,对尕尔勤矿床成矿斑岩年代学、岩石地球化学及锆石Hf同位素组成进行了详细研究,探讨了岩石成因及矿床形成构造背景。与此同时,对野外地质调查过程中于矿床北东部发现的具有强烈Cu、Au、Pb等成矿元素异常的硅帽进行了地化剖面、岩石地球化学及高光谱研究,探讨了存在浅成低温热液型铜金矿的可能性,以期为矿集区成岩成矿作用及其构造背景研究提供新素材的同时,为矿床找矿工作提供新的方向。

图 1 多龙矿集区地质简图 (据王勤等, 2015) 1-第四系;2-新近系康拖组;3-下白垩统美日切错组;4-中下侏罗统色哇组;5-下侏罗统曲色组;6-上三叠统日干配错组;7-二长花岗斑岩;8-花岗斑岩;9-石英斑岩;10-花岗闪长斑岩;11-辉长岩;12-橄榄岩;13-地层整合接触界线;14-地层不整合接触界线;15-断层;16-地表蚀变范围;17-矿床 (点);18-硅帽位置 Fig. 1 Geological map of the Duolong ore concentrated area (modified after Wang et al., 2015)
2 区域与矿床地质

多龙矿集区位于班公湖-怒江缝合带北缘的多不杂岩浆弧内,矿集区内已发现有多不杂、波龙、铁格隆南 (荣那)、拿若等大型-超大型斑岩或斑岩-浅成低温热液矿床,另还分布有尕尔勤、色那、拿厅、铁格山、地堡那木岗等具有较大找矿潜力的矿床点 (图 1)。

矿集区内出露地层主要可见有上三叠统日干配错组、下侏罗统曲色组 (J1q)、中-下侏罗统色哇组 (J1-2s)、下白垩统美日切错组 (K1m),另可见大面积分布的第四系 (Q)。其中下侏罗统曲色组及中-下侏罗统色哇组是矿集区内重要的赋矿围岩地层。矿集区内主要发育一套中酸性侵入岩,岩性主要可见花岗闪长斑岩、石英闪长玢岩、花岗斑岩、闪长玢岩、闪长岩等,其中以岩株或岩脉形式产出的花岗闪长斑岩、石英闪长玢岩为矿集区内重要的成矿斑岩。矿集区内构造主要可见断裂构造,可见有EW向、NE向、NW向断裂构造发育,其中地堡那木岗、波龙、多不杂、铁格隆南及拿若等矿床均沿北东向断裂构造分布,而尕尔勤、色那、赛角及拿若则沿北西向断裂构造分布。

尕尔勤斑岩铜矿床位于多龙矿集区东南部 (图 1),是矿集区内最具找矿潜力的矿床之一。前人已陆续进行了一些探矿工作,但未能取得大的突破。矿区内出露地层主要可见第四系 (Q)、下白垩统美日切组错下段 (K1m)、下侏罗统曲色组第一、二岩性段 (J1q1、J1q2) 以及上三叠统日干配错组 (T3r)(图 2)。曲色组第一岩性段主要由浅灰白色、浅黄褐色薄至中厚层状弱褐铁矿化、弱硅化长石石英砂岩组成,该套地层在矿区南西部大量出露;曲色组第二岩性段主要由浅灰色、浅黄褐色薄至中厚层状弱褐铁矿化、弱硅化长石石英砂岩组成,曲色组第二岩性段与下覆第一岩性地层呈整合接触关系。下白垩统美日切错组主要可见一套紫红色安山质火山碎屑岩、安山玢岩,主要分布在矿区的北部,分布范围较窄。上三叠统日干配错组展布于矿区北东部,范围较小,岩性为一套单一的灰黑色结晶灰岩,与曲色组第二岩性段呈断层接触关系。

图 2 尕尔勤矿床地质简图 (据西藏自治区地质调查院, 2011修改) Fig. 2 Geological map of the Gaerqin deposit

①  西藏自治区地质调查院. 2011.西藏班公湖-怒江成矿带西段铜多金属资源调查报告

矿区岩浆岩岀露较少,主要可见呈小岩株或岩脉形态展布的花岗闪长斑岩。构造主要可见断裂构造,矿区内可见有F1及F2两条断裂构造。F1断层展布于矿区中北部,近EW向,走向285°±,断层面北倾,倾角35°~70°之间。F2断层范围较窄,展布于矿区北东部,为日干配错组灰黑色结晶灰岩与曲色组第二岩性段的接触界线。

截止目前为止,矿区揭露矿体较小,均赋存于花岗闪长斑岩体中,呈分枝状,厚度不大,矿化不连续,矿体品位不高,属于斑岩型矿床中边部的枝状贫矿体。矿石矿物主要可见黄铜矿、黄铁矿,脉石矿物主要可见石英、长石、云母及粘土矿物。矿石构造主要可见细脉浸染状及块状构造,矿石结构主要可见自形粒状结构。矿区范围内热液蚀变较为明显,可见有黑云母化、硅化、钾长石化、绢云母化、高岭土化及绿泥石化等,鉴于探矿工程以及目前探获岩体形态的制约,没有明显的蚀变分带现象。

3 样品及分析方法 3.1 采样位置及样品特征

本次用于锆石定年及Hf同位素研究的样品采集于矿区中部的花岗闪长斑岩矿化体中 (图 2),斑岩体岩石地球化学样品分别采集于钻孔ZK2308(3件)、ZK2312(2件)。花岗闪长斑岩呈灰白色,斑状结构,块状构造,主要由斑晶 (约占40%) 与基质 (约占60%) 组成。斑晶主要由长石、石英、角闪石及少量黑云母组成,基质主要由长英质矿物组成。可见有硅化、钾化及粘土化蚀变 (图 3a-f)。

图 3 尕尔勤矿区斑岩体及硅帽照片 (a-c) 绢英岩化及钾化花岗闪长斑岩;(d-f) 为 (a-c) 对应的镜下照片;(g) 层纹状 (致密块状) 硅帽;(h、i) 角砾状硅帽,可见明显充填交代现象;(j) 为 (g) 的对应镜下照片;(k、l) 为 (h、i) 对应镜下照片 Fig. 3 Photos of porphyry and lithoscapss in Gaerqin deposit

本次研究所用硅帽样品采集于矿区北东部强硅化区,采样点见图 2。紫红色、红色、灰红色,可见有层纹状 (致密块状) 及角砾状构造两种 (图 3g-l)。层纹状硅帽发育较少,可见明显弧形层理,块状构造,质地坚硬。角砾状硅帽广泛发育,其硅质角砾呈棱角-次棱角状,磨圆度较差,粒径多在0.5~15mm之间,由铁质、硅质胶结而成 (图 3h, i, k, l),角砾状硅帽样品多发育孔洞,指示经历过强酸淋滤。

3.2 锆石U-Pb定年及Lu-Hf同位素测定

本次样品破碎和锆石挑选由河北省廊坊区域地质矿产调查研究所实验室完成,破碎后样品经淘洗除去比重轻的矿物,采用浮选和磁选分选出锆石,在双目镜下挑选晶形、色泽较好,透明度高的锆石颗粒进行制靶 (宋彪等, 2002),然后利用阴极发光 (CL) 显微照相观察锆石的内部结构,避开包裹体、裂隙及残留核,最后选择合适的测点进行分析。锆石U-Pb分析及Lu-Hf同位素测定均在中国地质科学院矿产资源研究所LA-ICP-MS实验室完成。锆石年龄计算以国际标准锆石91500为外标。测试结果用GLITTER (ver4.0) 软件计算得出,并按照Andersen (2002)的方法进行普通铅校正,之后用Isoplot程序完成年龄计算及谐和图绘制。具体分析步骤和数据处理过程参见文献 (Gao et al., 2003; 柳小明等, 2002; 袁洪林等, 2003; 侯可军等, 2009)。

Lu-Hf同位素分析采用配有193nm激光的Neptune多接收电感耦合等离子质谱仪进行分析,分析过程中采用6~8Hz的激光频率、100mJ的激光强度和50~60μm的激光束斑。试验中以He作为剥蚀物质的载气,采用91500作为外部标准,测试过程中,插入的91500标样的176Hf/177Hf比值变化于0.282310~0.282319之间,完全在仪器的正常测试范围内,具体测试过程见文献 (Wu et al., 2006)。

3.3 主量、微量元素分析及高光谱测定

主量元素和微量元素分析均在西南冶金地质测试所完成。主量元素采用X荧光光谱法 (XRF) 分析,分析精度优于2%;微量元素采用等离子质谱 (ICP-MS) 分析方法,测试精度优于5%。

高光谱利用美国ASD光谱分析仪进行,识别项目主要为扫面矿物识别,可获得混合矿物光谱曲线及各种矿物含量比重。

4 花岗闪长斑岩测试结果 4.1 锆石U-Pb定年结果

尕尔勤花岗闪长斑岩内锆石无色透明,长轴与短轴之比在1:1~2.5:1之间,长轴多大于100μm,颗粒大者长轴可达300μm,主要呈长柱状或短柱状,个别呈浑圆粒状,部分颗粒内部可见细小包裹体。阴极发光图像显示锆石成分含量比较均匀,发育密集震荡环带 (图 4),具岩浆成因锆石特征 (吴元保和郑永飞, 2004)。

图 4 尕尔勤矿区斑岩体锆石阴极发光图像及锆石U-Pb定年 (实线圈) 和Hf同位素 (虚线圈) 测点位置 Fig. 4 Zircon cathodoluminescence (CL) images and analysis spots of U-Pb (real coil) and Hf isotope (virtual coil) from porphyry of Gaerqin deposit

本次对花岗闪长岩样品挑选出的锆石晶体共进行了25个点的测试,锆石U-Pb同位素分析结果见表 1。样品中Th含量在69×10-6~270×10-6之间,平均值为148.3×10-6,U含量在289×10-6~784×10-6之间,平均值为477.7×10-6,Th/U比值在0.20~0.40之间,具岩浆成因锆石特征。样品测点年龄的计算是利用Isoplot软件进行处理,得到的锆石206Pb/238U加权平均年龄为124.4±0.4Ma (MSWD=0.42,n=25),该年龄代表了花岗闪长岩的结晶年龄,属于早白垩世侵位。在206Pb/238U-207Pb/235U谐和图上可以看出 (图 5),所有分析数据点均分布在谐和曲线上,拟合度较高,得到的年龄值真实可靠。

表 1 尕尔勤矿区斑岩体锆石定年数据结果 Table 1 U-Pb isotopes analyzed data of the zircons of porphyry in the Gaerqin deposit

图 5 尕尔勤矿区斑岩体锆石年龄 Fig. 5 Concordia diagrams of zircon U-Pb ages for porphyry of Gaerqin deposit
4.2 锆石Lu-Hf同位素组成

本次Lu-Hf同位素测试共进行了20个点的测试 (表 2),测试结果表明,花岗闪长斑岩内锆石176Hf/177Hf分布于之间0.282140~0.2829720之间,176Lu/177Hf变化于0.000631~0.002242之间。根据对应锆石年龄计算的εHf(t) 值变化于-19.75~9.71之间,主要集中于1.15~9.71之间;εHf(0) 主要分布于-22.36~7.07之间;地壳tDM模式年龄变化于398~1569Ma之间,tDMc模式年龄变化于559~2424Ma之间, 主要集中于559~1105Ma之间。

表 2 尕尔勤矿区斑岩体锆石Hf同位素数据结果 Table 2 Hf isotopes analyzed data of the zircons of porphyry in the Gaerqin deposit
4.3 地球化学结果 4.3.1 主量元素

本次测试所得花岗闪长斑岩主量元素结果见表 3。其SiO2含量在60.01%~62.81%之间;Al2O3含量普遍较高,在15.12%~16.51%之间,平均值16.17%;MgO含量在1.79%~2.14%之间,平均值1.95%;CaO含量在2.57%~3.85%之间,平均值3.28%;Na2O含量普遍高于K2O含量,Na2O/K2O在1.46~2.43之间,AR (碱度率)([Al2O3+CaO+(Na2O+K2O)]/[Al2O3+CaO-(Na2O+K2O)]) 在1.7~1.85之间。在SiO2-K2O图解及AR-SiO2图解中样品均落入了钙碱性区域 (图 6a, b),表明尕尔勤矿床成矿斑岩体属于钙碱性岩石系列。

表 3 尕尔勤矿区花岗闪长斑岩与硅帽主量元素 (wt%) 与微量元素 (×10-6) 分析结果表 Table 3 Major elements (wt%) and trace elements (×10-6) of porphyry in Gaerqin deposit

图 6 尕尔勤矿区斑岩体SiO2-K2O图解 (a, 据Peccerillo and Taylor, 1976)、SiO2-AR图解 (b, 底图据Wright, 1969)、球粒陨石标准化稀土元素配分图 (c, 标准化值据Boynton, 1984) 及原始地幔标准化微量元素蛛网图 (d, 标准化值据Sun and McDonough, 1989) 多不杂成矿斑岩数据据李金祥, 2008;波龙成矿斑岩数据据Li et al., 2013; 陈华安等, 2013 Fig. 6 SiO2-K2O plot (a, map after Peccerillo and Taylor, 1976), SiO2-AR plot (b, base map after wright, 1969), chondrite-normalized REE patterns (c, normalization values after Boynton, 1984) and primitive mantle-normalized trace element spidergrams (d, normalization values after Sun and McDonough, 1989) of porphyry in Gaerqin deposit
4.3.2 微量元素

花岗闪长斑岩微量元素分析结果见表 3,可以看出,花岗闪长斑岩微量元素具有以下特征:稀土总量偏低 (∑REE在54.3×10-6~72.07×10-6之间);LREE/HREE在6.3~10.2之间,(La/Yb)N介于6.27~14.38之间,指示轻、重稀土元素分馏较明显;δEu在0.82~1.20之间,无明显Eu异常,δCe在0.99~1.07之间,无明显Ce异常;稀土配分模式曲线显示为向右倾的轻稀土富集型 (图 6c)。

从微量元素蛛网图可以看出 (图 6d),Rb、Th、U、Sr等大离子亲石元素相对富集,而Nb、Ta、Ti、P、Zr等高场强元素则表现为相对亏损特征,体现出岛弧岩浆的特有的地球化学特征。

5 硅帽测试结果 5.1 主量元素

本次测试所得硅帽样品主量元素结果见表 3。层纹状 (致密块状) 硅帽样品SiO2含量在93.16%~95.99%之间;Al2O3含量在0.41%~0.61%之间,Si/Al比值在138~207之间;CaO含量在0.29%~0.67%之间;Na2O含量与K2O含量相当且都比较低;TiO2含量在0.54%~1.02%之间;MnO含量在0.031%~0.14%之间;P2O5含量较低,在0.027%~0.052%之间;Al/(Al+Fe+Mn) 在0.06~0.48之间,U/Th比值在0.33~0.92之间,Ni/Co比值在0.38~2.67之间。角砾状硅帽样品SiO2含量明显低于层纹状样品,其SiO2含量在79.12%~89.77%之间,Al2O3含量在0.45%~0.48%之间,Si/Al比值在146~176之间;CaO含量在0.58%~0.59%之间;TiO2含量普遍低于层纹状硅帽,含量在0.16%~0.24%之间;MnO含量在0.074%~0.31%之间;P2O5含量较低,在0.051%~0.079%之间。

5.2 微量元素

本次研究硅帽样品微量与稀土元素结果见表 3。从稀土元素结果可以看出,两种硅帽样品的∑REE含量具有明显差别,层纹状硅帽样品∑REE含量在109.6×10-6~182.3×10-6之间,平均值139.5×10-6,明显高于角砾状硅帽样品 (∑REE含量在29.31×10-6~35.93×10-6之间,平均值31.71×10-6)。两种硅帽样品因稀土含量的不同其标准化曲线存在着截然的差别,从球粒陨石标准化曲线结果可以看出 (图 7b),两类硅帽样品均表现为轻、重稀土分异明显的右倾型,均具有弱负Ce异常,但层纹状硅帽样品表现为弱负Eu异常,而角砾状硅帽则表现出弱正Eu异常。在北美页岩稀土元素标准化图上可以看出 (图 7a),两种硅帽样品均表现出轻、重稀土分异不明显的平坦分布模式,层纹状硅帽样品表现出弱正Eu异常,而角砾状硅帽样品则表现出强烈正Eu异常,Ce异常不明显。综合显示两类硅帽样品可能是两个不同时期不同环境的产物。

图 7 尕尔勤硅帽样品北美页岩标准化 (a) 与球粒陨石标准化 (b) 稀土元素分布型式图 (标准化值据Boynton, 1984) Fig. 7 North American shale-normalized (a) and chondrite-normalized (b) REE patterns of lithoscaps in Gaerqin deposit (normalization values after Boynton, 1984)

矿区硅帽样品总体表现出微量元素低的特点,相对于地壳元素丰度而言 (表 3),Nb、Ta、Th、Ba等不相容元素亏损较为明显,而Cu、Pb、Zn、Sb、Y、As等亲铜元素则表现为富集特征,含量均在克拉克值的1.5倍以上,特别是角砾岩型硅帽中Cu、Zn、As、Sb等元素表现为强烈富集 (图 8)。此外Ni、Co、Rb、V、Ba等也表现为亏损特征。

图 8 两种硅帽样品部分微量元素含量与地壳丰度对比图 Fig. 8 Contrast graphic of some trace elements and Clarke values of two kinds of lithoscaps
5.3 高光谱

本次利用美国ASD光谱仪对矿区角砾状硅帽样品进行了高光谱测试,重点针对粘土矿物类矿物的识别。本次共进行2个样品的分析测试,每件样品测试3个面,获得光谱解混结果,可半定量获得矿物组成。本次测试结果识别出了埃洛石、伊利石、高岭石、方解石等矿物 (图 9),主要为泥化带/氧化带蚀变产物。

图 9 尕尔勤矿区角砾状硅帽样品高光谱分析结果 Fig. 9 Hyperspectral analysis of lithoscaps in Gaerqin deposit
6 讨论 6.1 斑岩体岩石成因及其源区探讨

岩石中如存在H2O、CO2等挥发份会致使其在烧灼时逃逸形成烧失量 (LOI),其中CO2常以碳酸盐矿物的形式出现,H2O则多与含水矿物有关。蚀变是斑岩型矿床的必备条件,其发育的绢英岩化及青磐岩化蚀变带中均具有大量含水矿物及一定的碳酸盐矿物,因此斑岩型矿床成矿斑岩在分析过程中往往会出现相对较高的烧失量。本次用于测试样品均选择蚀变相对较弱的进行,但烧失量依然较高。从Harker图解及岩石分析数据表中可以看出 (图 10表 3),LOI与CaO的含量呈负相关关系,表明岩石烧失量主要是与含水蚀变矿物有关;LOI与Na2O及Sr的含量呈微弱的负相关关系,而与K2O的含量则成微弱的正相关关系,表明岩石具有一定程度的钾化及绢英岩化,从而造成K的相对升高而Na、Sr的相对降低,但相关性并不明显,表明蚀变程度并不高;LOI与Y、Yb、La等稀土元素,Rb、Ba等多数大离子亲石元素元素以及高场强元素并无明显的相关性 (图 10表 3),表明稀土元素及大部分微量元素可信度均较高。

图 10 尕尔勤矿床花岗闪长斑岩的Harker图解 Fig. 10 Harker diagrams for granodiorite porphyry from Gaerqin deposit

尕尔勤矿区与成矿关系密切的花岗闪长斑岩具有以下地球化学特征:SiO2含量均大于60%,Al2O3含量均大于15%,MgO含量较低 (1.79%~2.14%之间),富钠而贫钾 (Na2O含量在3.03%~3.82%之间,K2O含量在1.57%~2.16%之间),高Sr (312×10-6~503×10-6之间)、低Y (10.4×10-6~11.2×10-6之间) 和Yb (0.76×10-6~1.34×10-6之间),Sr/Y比值在29~49之间,Rb、Th、U、Sr等大离子亲石原始元素相对富集,而Nb、Ta、Ti、P、Zr等高场强元素则表现为相对亏损特征,与典型的埃达克岩地球化学特征相似 (Defant and Drummond, 1990; 王焰等, 2000; 王强等, 2001; 张旗等, 2001; Castillo, 2006),因此可将尕尔勤矿区成矿斑岩体定为埃达克岩 (图 11a)。

图 11 尕尔勤矿区斑岩体Sr/Y-Y图解 (a, 据Defant and Drummond, 1990) 和锆石U-Pb年龄-εHf(t) 图解 (b) 多不杂数据据李金祥, 2008; 波龙数据据陈华安等, 2013; Li et al., 2013 Fig. 11 Sr/Y vs. Y plot (a, base map after Defant and Drummond, 1990) and U-Pb ages vs. εHf(t) diagram (b) of porphyry in Gaerqin deposit

已有研究表明,埃达克岩存在着O型埃达克岩与C型埃达克岩之分。O型埃达克岩形成于俯冲玄武质洋壳的部分熔融 (Defant and Drummond, 1990; 张旗等, 2001),而C型埃达克岩则形成于下地壳底侵玄武质岩石的部分熔融或加厚下地壳的部分熔融 (Atherton and Petford, 1993; Wang et al., 2005, 2007)。尕尔勤矿区斑岩体富Na而贫K,且具有低的Th含量 (3.93×10-6~4.28×10-6之间) 及Th/Ce比值 (0.12~0.20之间),与O型埃达克岩地球化学特征相似而区别于C型埃达克岩。上述特点与矿集区内临近的多不杂、波龙的大型斑岩矿床成矿斑岩体特征十分相似 (李金祥等, 2008; 陈华安等, 2013),表明尕尔勤花岗闪长斑岩来源于俯冲洋壳的部分熔融 (Defant and Drummond, 1990)。

矿区斑岩体轻重稀土分异明显,重稀土相对亏损,具有高Sr,低Y、Yb的特点,暗示源区内有角闪石和石榴石的残余 (Defant and Drummond, 1990; Condie, 2005; Castillo, 2006; Davidson et al., 2007)。锆石Hf同位素体现出总体分布范围较大而相对集中的特点,其εHf(t) 值主要集中于1.15~9.71之间,仅3个样品εHf(t) 值小于0,在εHf(t) 与锆石-Pb年龄图解中可以看出 (图 11b),矿区花岗闪长斑岩体Hf同位素组成绝大多数靠近于亏损地幔线,表明其为起源于具有幔源性质的深成熔体,与俯冲玄武质洋壳部分熔融形成埃达克岩一致 (Stern and Kilian, 1996),但个别较大负εHf(t) 值的出现表明其具有轻微的古老地壳的混染,这一特点与邻近的多不杂、波龙等矿床成矿斑岩体有所区别 (李金祥, 2008; Li et al., 2013; 陈华安等, 2013)。

已有研究资料表明,多龙矿集区内与成矿有关的斑岩体多具有埃达克岩地球化学特征。李金祥等 (2008)对多不杂矿区含矿斑岩体地球化学特征研究结果表明其属于O型埃达克岩,起源于俯冲洋壳的部分熔融;Li et al. (2013)陈华安等 (2013)对波龙矿区含矿斑岩体地球化学特征研究结果指出其与多不杂矿床具有相似的岩浆源区。尕尔勤矿床与邻近的波龙、多不杂矿床相比,与其成矿有关的斑岩体微量元素与稀土元素地球化学特征 (图 6c, d) 及Hf同位素组成均非常相似 (图 11b),成岩年龄相当 (尕尔勤矿床成矿斑岩成岩年龄为124.34±0.36Ma;多不杂矿床成矿斑岩成岩年龄为121.6±1.9Ma;波龙矿床成矿斑岩成岩年龄为118~120Ma)(李金祥等, 2008; Li et al., 2013; 陈华安等, 2013),表明三者可能具有相似的岩浆源区与岩石成因,即矿区成矿斑岩的侵位与班公湖-怒江特提斯洋壳的俯冲消减密切相关。

6.2 矿床形成的构造背景探讨

班公湖-怒江缝合带作为冈瓦纳大陆的北界 (潘桂棠等, 2004),是拉萨陆块与冈底斯陆块的分界线 (任纪舜和肖黎薇, 2004),其构造演化问题一直是西藏地质学界争论的焦点。缝合带SSZ型蛇绿岩的研究已表明班公湖-怒江特提斯洋至少在中侏罗世已开始由扩张转为俯冲消减 (史仁灯, 2004; 邱瑞照等, 2004; 周涛等, 2014),但对于班公湖-怒江特提斯洋的俯冲消减极性以及碰撞闭合的时限一直存在较大的争论。前期研究者认为俯冲形式为班公湖-怒江特提斯洋洋壳向羌塘板块下的单向俯冲 (Pearce and Deng, 1988; Kapp et al., 2003),但随着该构造带研究资料的丰富及研究的深入,越来越多的学者认为该构造带具有向南、向北双重俯冲的特征 (秦克章等, 2006; 李光明等, 2007; 许荣科等, 2007; 李金祥等, 2008; 康志强等, 2010; 杜德道等, 2011; 耿全如等, 2011)。目前“双向俯冲”的观点已基本被大多数研究者所认同,但关于俯冲→碰撞缝合阶段的演化却又存在不同的观点,主要集中在晚侏罗世末-早白垩世初缝合带是否已经闭合的问题上。Kapp et al. (2003)根据缝合带沉积相的分析认为班公湖-怒江特提斯洋在145Ma左右已经闭合消失,该观点得到了后来申扎-班戈一带碰撞后A型花岗岩时代的佐证 (曲晓明等, 2012, 2013)。另一方面,也有众多研究表明班公湖-怒江特提斯洋在早白垩世初还没消亡,洋盆闭合时限可能推迟到早白垩世中晚期,主要研究资料包括:王忠恒等 (2005)于该构造带带中段扎加藏布下游塔仁本发现了早白垩世 (110Ma左右) 近百余平方千米的洋岛,表明在班公湖-怒江特提斯洋中段在早白垩世中晚期还未消亡;朱弟成等 (2006)通过塔仁本和多玛OIB型洋岛玄武岩的研究,指出班公湖-怒江洋壳在大约110Ma时尚未彻底消亡;李金祥等 (2008)Li et al.(2011a, 2013) 通过对藏北多龙矿集区内多不杂及波龙斑岩矿床成矿斑岩的年代学及地球化学特征研究,指出班公湖-怒江特提斯洋在120Ma左右仍在向北俯冲;康志强等 (2009)通过对拉萨地块北部于早白垩世中期形成的多尼组火山岩的地球化学及锆石年代学研究,认为该套火山岩更大的可能性形成于班公湖-怒江特提斯洋的南向俯冲环境;高顺宝等 (2011)对班戈地区侵入岩进行了锆石U-Pb年代学及地球化学研究,指出班公湖-怒江特提斯洋向南俯冲可能早于中侏罗世,而在早白垩世中期才碰撞闭合;耿全如等 (2011)在综合总结前人研究资料及野外实地调查基础上,提出了班公湖-怒江特提斯洋在中侏罗世至早白垩世向北、向南双向俯冲,在晚白垩世碰撞缝合的观点;赵元艺等 (2011)通过对舍索铜矿成矿母岩的年代学及地球化学研究,指出成岩于116Ma的花岗闪长岩形成于班公湖-怒江洋南向俯冲闭合之前的岛弧环境;Zhu et al. (2011)通过对拉萨地体全面的酸性岩浆岩Hf同位素研究,指出113Ma左右班公湖-怒江洋盆已消失,深俯冲的洋壳发生了板片断离,表明此时已进入碰撞阶段;Li et al. (2014)通过对南羌塘中酸性侵入岩进行地球化学、锆石年代学及Hf同位素研究表明班公湖-怒江特提斯洋壳从俯冲至碰撞时限为170~110Ma。综上所述研究表明,班公湖-怒江缝合带构造演化是一个极其复杂的问题,并不是一个简单的年龄所能判别,鉴于已有众多资料表明区域上存在着大量的与班公湖-怒江特提斯洋俯冲相关的早白垩世中晚期岛弧岩浆岩,因此笔者更倾向于认为班公湖-怒江特提斯洋在早白垩世末仍未消亡,至早白垩世中晚期方才碰撞闭合。

多龙矿集区作为班公湖-怒江成矿带上最为重要的超大型铜金矿集区,其构造背景问题同样存在争论,目前有俯冲背景下的陆缘弧或增生楔环境 (李金祥等, 2008; 李光明等, 2011)、洋壳消亡后陆-陆碰撞过程中的晚阶段-碰撞后伸展环境 (曲晓明等, 2009, 2012, 2013) 两种不同的观点。本次对矿集区内重要的矿床—尕尔勤矿床成矿斑岩研究结果表明,其具有相对富集Rb、Th、U、Sr等大离子亲石元素,而相对亏损Nb、Ta、Ti、P、Zr等高场强元素的岛弧岩浆特征,在Rb-Y+Nb及Rb-Y+Ta图解中所有样品均落入了岛弧花岗岩区域 (图 12a, b)。同时,斑岩体显示出明显的高Sr和低Y、Yb的O型埃达克岩特点,Hf同位素组成也表现为起源于深成熔体的源区特征,综合表明尕尔勤矿区成矿斑岩体形成于岛弧环境,矿床形成于班公湖-怒江洋壳向北俯冲阶段,其成岩成矿作用与洋壳俯冲密切相关,该认识与区内多不杂、波龙等矿床已有研究结果一致 (李金祥等, 2008; 李光明等, 2011; Li et al., 2011; 陈华安等, 2013),表明多龙矿集区成岩-成矿作用总体上可能形成于岛弧俯冲环境。尕尔勤斑岩体成岩年龄为124.4±0.4Ma,多龙矿集区内成矿斑岩体基本集中在120Ma左右,结合前述区域构造演化研究,笔者认为至少在早白垩世中期班公湖-怒江特提斯洋壳还处于向北俯冲阶段。

图 12 尕尔勤矿区斑岩体Rb-Y+Nb (a) 和Rb-Yb+Ta (b) 图解 (底图据Pearce et al., 1984) Fig. 12 Rb vs. Y+Nb (a) and Rb vs. Yb+Ta (b) diagrams of porphyry in Gaerqin deposit (base map after Pearce et al., 1984)
6.3 硅帽的识别与可能的浅成低温热液铜金矿床

科学研究的最终落脚点在于生产与社会经济效益,因此正确的科学理论应用于野外地质调查对于圈定找矿靶区来说至关重要。硅帽 (或岩帽) 是指在地质环境上处于古地表与浅成中酸性侵入体之间,富硅的,经受了高级泥化蚀变的斑岩型热液活动/蚀变系统的顶部 (Sillitoe, 1973, 1983, 2010),硅帽的形成与识别往往对浅成低温热液型及隐伏斑岩型矿床具有重要的勘探意义 (Sillitoe et al., 1998; 徐庆生等, 2010),但却往往被人忽视。早在2006年,秦克章等 (2006)便提出了西藏班公湖-怒江成矿带具备有产出斑岩-浅成低温热液矿床的地质条件,但鉴于藏北高原总体工作程度的制约,该类矿床找矿工作一直未能取得突破,直到铁格隆南超大型斑岩-浅成低温热液矿床的发现 (唐菊兴等, 2014),才使人们意识到班公湖-怒江成矿带在该类型矿床中的巨大找矿潜力,而多龙矿集区便是其中的重点突破地区。

本次研究在野外地质调查过程中,于尕尔勤矿床北东角突起的山脊处发现了高硅质岩石,前人一直未进行工作,本次发现的高硅质岩石具有层纹状 (致密块状) 及角砾状两种产出,角砾状高桂质岩石大量分布,而层纹状则相对较少,且呈凌乱分布。角砾状高硅质岩石可见明显的孔洞,表明其经受了强烈的风化淋滤,岩石化学特征表明,角砾状高硅质岩石相对贫SiO2(79.12%~89.77%之间),稀土元素含量较低 (∑REE含量在29.31×10-6~35.93×10-6之间);镜下鉴定结果表明,角砾状高硅质岩石具有明显充填交代现象 (图 3K, L),角砾磨圆度差,泥化蚀变矿物普遍存在,手标本、岩石化学特征及镜下特征均表明角砾状高硅质岩石明显不同于生物沉积成因或热水沉积成因硅质岩 (Herzig et al., 1988; Yin and Harrison, 2000; 冯彩霞, 2011), 而属于原岩经受风化淋滤后经热液充填交代成因的硅化岩 (即硅帽)。层纹状高硅质岩石稀土元素特征与区域上已报到的改则北硅质岩具有一定的相似之处 (冯彩霞, 2011),但其具有更为宽泛的Al/(Al+Fe+Mn) 值、低的U/Th比值以及低的Ni/Co比值,表明层纹状硅化岩可能是生物化学沉积于热水沉积共同作用的产物 (Rona, 1978; Crerar et al., 1982),鉴于资料的约束本文不作深入讨论。角砾状硅帽应晚于层纹状硅帽形成,这种类型高硅质岩石的形成往往与断裂构造活动或火山活动有关,鉴于硅帽展布区域与日干配错组灰岩呈断层接触,同时硅帽边部也发现了零星凝灰岩分布,推测该区域大量分布的角砾状硅帽为火山活动与断裂构造双重影响的结果。

前述研究已表明,角砾状硅帽样品强烈富集Cu、Pb、Zn、Sb、Y、As、Ga等亲铜元素,表明其具有一定的勘查指示意义。本次研究对识别出的硅帽地区进行了十字地化剖面测量 (图 13),结果显示异常较为明显。在A-A’剖面中,Cu、Au、Ag、Pb、Zn、As、Sb、Hg等成矿指示性元素总体含量值总体均较大,而且高值异常区非常明显,且异常叠合完好,其中Cu、Au、Ag、Pb、Zn中心异常基本完全套合,而As、Sb、Hg异常则紧靠前述元素外围。A-A’剖面异常位置靠近矿区F2断层,暗示该异常可能与该断裂活动有关。在B-B’剖面中,Cu、Au、Ag、Pb、Zn、As、Sb、Hg等成矿指示性元素都具有高含量特征,但无明显高值异常区域,与A-A’剖面中西段类似,表明尽管该区域成矿指示性元素具有一定程度的富集,但在靠近F2断层的硅帽区域富集程度更强,是探矿工程布设的首选位置。

图 13 尕尔勤矿床硅帽地区A-A’及B-B’剖面 Fig. 13 The A-A' and B-B' sections of the lithoscapss area in Gaerqin deposit

本次高光谱测试识别出了埃洛石、伊利石、高岭石、方解石等泥化带矿物,前述矿物与强硅化组成了顶部高级泥化带,结合前述地化剖面的异常特征显示其具有存在浅成低温热液或斑岩-浅成低温热矿床的巨大可能性。目前并未识别出明矾石、绢云母、冰长石等更具有进一步指示意义的标志性矿物,是否属于高硫或低硫型还需要进一步工作去判别,但总体表明硅帽区域具有重要的勘查意义。值得注意的是,在与强烈富集成矿指示元素的硅帽邻近的日干配错组灰岩中发现了大量的由铁锰碳酸盐矿物组成的网脉 (图 14a),经后期氧化形成团块状褐铁矿、软锰矿等矿物 (图 14b),表明该区域流体活动强烈,属于典型的流体逃逸现象,这也是硅帽区域值得进一步进行勘查工作的侧面证明。

图 14 尕尔勤日干配错组灰岩中的流体逃逸现象 Fig. 14 Phenomenon of the fluidescaping of the Riganpeicuo Formation in Gaerqin deposit
7 结论

(1) 利用锆石U-Pb年代学及Hf同位素组成研究,查明了尕尔勤矿区内成矿斑岩体成岩年龄为124.4±0.4Ma,侵位于早白垩世;锆石εHf(t) 值变化于-19.75~9.71之间,tDMc模式年龄变化于559~2424Ma之间,显示为起源于具有幔源性质的深成熔体。矿床的形成与班公湖-怒江洋北向俯冲密切相关。

(2) 区内成矿斑岩主量元素显示为钙碱性岩石系列,稀土元素表现为向右倾的轻稀土富集型分布型式;微量元素显示为Rb、Th、U、Sr等大离子亲石元素相对富集,而Nb、Ta、Ti、P、Zr等高场强元素相对亏损特征,体现出岛弧岩浆特有的地球化学特征,表明区内成矿斑岩体形成于岛弧环境。

(3) 区内硅帽样品存在着两种地球化学及镜下特征迥异的类型,其中大面积分布的角砾岩型硅帽显示为热液充填交代成因。高光谱测试结果识别出高级泥化带矿物组合,地球化学剖面测量显示硅帽区域存在着明显的Cu、Au、Ag、Pb、Zn、As、Sb等成矿指示元素异常。综合研究表明该区域存在浅成低温热液或斑岩-浅成低温热液型矿床的巨大可能性,具有重要的勘查指示意义。

致谢 参加野外工作的还有中国地质科学院矿产资源研究所宋扬、刘志博助理研究员,成都理工大学郭娜副教授以及高轲、韦少岗、唐楠、汪重午等同学;野外工作过程中得到了中国科学院地质与地球物理研究所李光明副研究员及赵俊兴博士等在思路上的建设性启发;在实验过程中得到了中国地质科学院矿产资源研究所LA-ICP-MS实验室侯可军老师的帮助;李金祥副研究员和赵俊兴副研究员悉心审稿,提出中肯的意见与建议;在此一并致以衷心的感谢!
参考文献
[] Andersen T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb. Chemical Geology, 192(1-2): 59–79. DOI:10.1016/S0009-2541(02)00195-X
[] Atherton MP, Petford N. 1993. Generation of sodium-rich magmas from newly underplated basaltic crust. Nature, 362(6416): 144–146. DOI:10.1038/362144a0
[] Boynton WV. 1984. Cosmochemistry of the rare earth elements:Meteorite studies. In:Henderson P (ed.). Developments in Geochemistry. Amsterdam:Elsevier, 2:63-114
[] Castillo PR. 2006. An overview of adakite petrogenesis. Chinese Science Bulletin, 51(3): 257–268. DOI:10.1007/s11434-006-0257-7
[] Chen HA, Zhu XP, Ma DF, Huang HX, Li GM, Li YB, Li YC, Wei LJ, Liu CQ. 2013. Geochronology and geochemistry of the Bolong porphyry Cu-Au deposit, Tibet and its mineralizing significance. Acta Geologica Sinica, 87(10): 1593–1611.
[] Chen YJ. 2013. The development of continental collision metallogeny and its application. Acta Petrologica Sinica, 29(1): 1–17.
[] Condie KC. 2005. TTGs and adakites:Are they both slab melts?. Lithos, 80(1-4): 33–44. DOI:10.1016/j.lithos.2003.11.001
[] Crerar DA, Namson J, So Chyi M, Williams L, Feigenson MD. 1982. Manganiferous cherts of the Franciscan assemblage:Ⅰ. General geology, ancient and modern analogues, and implications for hydrothermal convection at oceanic spreading centers. Economic Geology, 77(3): 519–540. DOI:10.2113/gsecongeo.77.3.519
[] Davidson J, Turner S, Handley H, Macpherson C, Dosseto A. 2007. Amphibole "sponge" in arc crust?. Geology, 35(9): 787–790. DOI:10.1130/G23637A.1
[] Defant MJ, Drummond MS. 1990. Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature, 347(6294): 662–665. DOI:10.1038/347662a0
[] Du DD, Qu XM, Wang GH, Xin HB, Liu ZB. 2011. Bidirectional subduction of the Middle Tethys oceanic basin in the west segment of Bangonghu-Nujiang suture, Tibet:Evidence from zircon U-Pb LAICPMS dating and petrogeochemistry of arc granites. Acta Petrologica Sinica, 27(7): 1993–2002.
[] Feng CX. 2011. Geochemical characteristics and sedimentary environment of siliceous rocks of Gerze area in western Bangong Co-Nujiang suture zone, Tibet. Mineral Deposits, 30(5): 773–786.
[] Gao S, Liu XM, Yuan HL, Hattendorf B, Günther D, Chen L, Hu SH. 2003. Analysis of forty-two major and trace elements in USGS and NIST SRM Glasses by LA-ICPMS. Geochimica et Cosmochimica Acta, 67(Suppl.18): A116.
[] Gao SB, Zheng YY, Wang JS, Zhang Z, Yang C. 2011. The geochronology and geochemistry of intrusive rocks in Bange area:Constraints on the evolution time of the Bangong Lake-Nujiang ocean basin. Acta Petrologica Sinica, 27(7): 1973–1982.
[] Geng QR, Pan GT, Wang LQ, Peng ZM, Zhang Z. 2011. Tethyan evolution and metallogenic geological background of the Bangong Co-Nujiang belt and the Qiangtang massif in Tibet. Geological Bulletin of China, 30(8): 1261–1274.
[] Herzig PM, Becker KP, Stoffers P, Bäcker H, Blum N. 1988. Hydrothermal silica chimney fields in the Galapagos Spreading Center at 86°W. Earth and Planetary Science Letters, 89(3-4): 261–272. DOI:10.1016/0012-821X(88)90115-X
[] Hou KJ, Li YH, Tian YR. 2009. In situ U-Pb zircon dating using laser ablation-multi ion counting-ICP-MS. Mineral Deposits, 28(4): 481–492.
[] Hou ZQ, Ma HW, Zaw K, Zhang YQ, Wang MJ, Wang Z, Pan GT, Tang RL. 2003. The Himalayan Yulong porphyry copper belt:Product of large-scale strike-slip faulting in eastern Tibet. Economic Geology, 98(1): 125–145.
[] Hou ZQ, Gao YF, Meng XJ, Qu XM, Huang W. 2004. Genesis of adakitic porphyry and tectonic controls on the Gangdese Miocene porphyry copper belt in the Tibetan orogen. Acta Petrologica Sinica, 20(2): 239–248.
[] Hou ZQ, Zhong DL, Deng WM and Zaw K. 2005. A tectonic model for porphyry copper-molybdenum-gold deposits in the eastern Indo-Asian collision zone. In:Porter TM (ed.). Super Porphyry Copper & Gold Deposit:A Global Perspective. Adelaide:PGC Publishing, 423-440
[] Hou ZQ, Yang ZM, Qu XM, Meng XJ, Li ZQ, Beaudoin G, Rui ZY, Gao YF, Zaw K. 2009. The Miocene Gangdese porphyry copper belt generated during post-collisional extension in the Tibetan orogen. Ore Geology Reviews, 36(1-3): 25–51. DOI:10.1016/j.oregeorev.2008.09.006
[] Ji WQ, Wu FY, Chung SL, Li JX, Liu CZ. 2009. Zircon U-Pb geochronology and Hf isotopic constraints on petrogenesis of the Gangdese batholith, southern Tibet. Chemical Geology, 262(3-4): 229–245. DOI:10.1016/j.chemgeo.2009.01.020
[] Jiang JH, Wang RJ, Qu XM, Xin HB, Wang ZZ. 2011. Crustal extension of the Bangong Lake arc zone, western Tibetan Plateau, after the closure of the Tethys Oceanic Basin. Earth Science, 36(6): 1021–1032.
[] Kang ZQ, Xu JF, Wang BD, Dong YH, Wang SQ, Chen JL. 2009. Geochemistry of Cretaceous volcanic rocks of Duoni Formation in northern Lhasa block:Discussion of tectonic setting. Earth Science, 34(1): 89–104.
[] Kang ZQ, Xu JF, Wang BD, Chen JL. 2010. Qushenla Formation volcanic rocks in north Lhasa block:Products of Bangong Co-Nujiang Tethy's southward subduction. Acta Petrologica Sinica, 26(10): 3106–3116.
[] Kapp P, Murphy MA, Yin A, Harrison TM, Ding L, Gou JH. 2003. Mesozoic and Cenozoic tectonic evolution of the Shiquanhe area of western Tibet. Tectonics, 22(4): 1029.
[] Li GM, Qin KZ, Ding KS, Liu TB, Li JX, Wang SH, Jiang SY, Zhang XC. 2006. Geology, Ar-Ar age and mineral assemblage of Eocene skarn Cu-Au±Mo deposits in the southeastern Gangdese arc, Southern Tibet:Implications for deep exploration. Resource Geology, 56(3): 315–336. DOI:10.1111/rge.2006.56.issue-3
[] Li GM, Li JX, Qin KZ, Zhang TP, Xiao B. 2007. High temperature, salinity and strong oxidation ore-forming fluid at Duobuza gold-rich porphyry copper deposit in the Bangonghu tectonic belt, Tibet:Evidence from fluid inclusions. Acta Petrologica Sinica, 23(5): 935–952.
[] Li GM, Duan ZM, Liu B, Zhang H, Dong SL, Zhang L. 2011. The discovery of Jurassic accretionary complexes in Duolong area, northern Bangong Co-Nujiang suture zone, Tibet, and its geologic significance. Geological Bulletin of China, 30(8): 1256–1260.
[] Li GM, Li JX, Qin KZ, Duo J, Zhang TP, Xiao B, Zhao JX. 2012a. Geology and hydrothermal alteration of the Duobuza gold-rich porphyry copper district in the Bangongco metallogenetic belt, Northwestern Tibet. Resource Geology, 62(1): 99–118. DOI:10.1111/rge.2012.62.issue-1
[] Li JX. 2008. Geochronology, petrology and metallogeneses of high oxidized magma-hydrothermal fluid of Duobuza gold-rich porphyry copper deposit in Bangonghu belt, Northern Tibet. Ph. D. Dissertation. Beijing:Institute of Geology and Geophysics, Chinese Academy of Sciences, 1-224 (in Chinese)
[] Li JX, Li GM, Qin KZ, Xiao B. 2008. Geochemistry of porphyries and volcanic rocks and ore-forming geochronology of Duobuza gold-rich porphyry copper deposit in Bangonghu belt, Tibet:Constraints on metallogenic tectonic settings. Acta Petrologica Sinica, 24(3): 531–543.
[] Li JX, Qin KZ, Li GM, Xiao B, Zhao JX, Chen L. 2011a. Magmatic-hydrothermal evolution of the Cretaceous Duolong gold-rich porphyry copper deposit in the Bangongco metallogenic belt, Tibet:Evidence from U-Pb and 40Ar/39Ar geochronology. Journal of Asian Earth Sciences, 41(6): 525–536. DOI:10.1016/j.jseaes.2011.03.008
[] Li JX, Qin KZ, Li GM, Xiao B, Chen L, Zhao JX. 2011b. Post-collisional ore-bearing adakitic porphyries from Gangdese porphyry copper belt, southern Tibet:Melting of thickened juvenile arc lower crust. Lithos, 126(3-4): 265–277. DOI:10.1016/j.lithos.2011.07.018
[] Li JX, Li GM, Qin KZ, Xiao B, Chen L, Zhao JX. 2012b. Mineralogy and mineral chemistry of the Cretaceous Duolong gold-rich porphyry copper deposit in the Bangongco arc, Northern Tibet. Resource Geology, 62(1): 19–41. DOI:10.1111/rge.2012.62.issue-1
[] Li JX, Qin KZ, Li GM, Xiao B, Zhao JX, Cao MJ, Chen L. 2013. Petrogenesis of ore-bearing porphyries from the Duolong porphyry Cu-Au deposit, central Tibet:Evidence from U-Pb geochronology, petrochemistry and Sr-Nd-Hf-O isotope characteristics. Lithos, 160-161: 216–227. DOI:10.1016/j.lithos.2012.12.015
[] Li JX, Qin KZ, Li GM, Richards JP, Zhao JX, Cao MJ. 2014. Geochronology, geochemistry, and zircon Hf isotopic compositions of Mesozoic intermediate-felsic intrusions in central Tibet:Petrogenetic and tectonic implications. Lithos, 198-199: 77–91. DOI:10.1016/j.lithos.2014.03.025
[] Li T. 1992. The statistical characteristics of the abundance of chemical elements in the earth's crust. Geology and Prospecting, 28(10): 1–7.
[] Li YB, Duo J, Zhong WT, Li YC, Qiangba WD, Chen HQ, Liu HF, Zhang JS, Zhang TP, Xu ZZ, Fan AH, Suolang WQ. 2012. An exploration model of the Duobuza porphyry Cu-Au deposit in Gaize Country, northern Tibet. Geology and Exploration, 48(2): 274–287.
[] Liu XM, Gao S, Yuan HL, Hattendorf B, Günther D, Chen L, Hu SH. 2002. Analysis of 42 major and trace elements in glass standard reference materials by 193nm La-ICPMS. Acta Petrologica Sinica, 18(3): 408–418.
[] Pan GT, Wang LQ, Zhu DC. 2004. Thoughts on some important scientific problems in regional geological survey of the Qinghai-Tibet Plateau. Geological Bulletin of China, 23(1): 12–19.
[] Pearce JA, Harris NBW, Tindle AG. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. Journal of Petrology, 25(4): 956–983. DOI:10.1093/petrology/25.4.956
[] Pearce JA, Deng WM. 1988. The ophiolites of the Tibetan Geotraverses, Lhasa to Golmud (1985) and Lhasa to Kathmandu (1986). Philosophical Transactions of the Royal Society A:Mathematical, Physical and Engineering Sciences, 327(1594): 215–238. DOI:10.1098/rsta.1988.0127
[] Peccerillo A, Taylor SR. 1976. Geochemistry of eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey. Contributions to Mineralogy and Petrology, 58(1): 63–81. DOI:10.1007/BF00384745
[] Qin KZ, Li GM, Zhang Q, Li JX, Miao Y, Xiao B, Zhang TP, Duo J, Li JG and Lu Y. 2006. Metallogenic conditions and possible occurrences for epithermal gold mineralizations in Gangdese and Bangonghu belts, Tibet:In view of porphyry-epithermal Cu-Au metallogenetic systematic. In:Proceedings of 8th State Mineral Deposit Conference. Beijing:Geological Publishing House, 660-670 (in Chinese)
[] Qiu RZ, Zhou S, Deng JF, Li JF, Xiao QH, Cai ZY. 2004. Dating of gabbro in the Shemalagou ophiolite in the western segment of the Bangong Co-Nujiang ophiolite belt, Tibet:With a discussion of the age of the Bangong Co-Nujiang ophiolite belt. Geology in China, 31(3): 262–268.
[] Qu XM, Hou ZQ, Zaw K, Li YG. 2007. Characteristics and genesis of Gangdese porphyry copper deposits in the southern Tibetan Plateau:Preliminary geochemical and geochronological results. Ore Geology Reviews, 31(1-4): 205–223. DOI:10.1016/j.oregeorev.2005.03.012
[] Qu XM, Hou ZQ, Zaw K, Mo XX, Xu WY, Xin HB. 2009. A large-scale copper ore-forming event accompanying the rapid uplift of the southern Tibetan Plateau:Evidence from zircon SHRIMP U-Pb dating and LA ICP-MS analysis. Ore Geology Reviews, 36(1-3): 52–64. DOI:10.1016/j.oregeorev.2009.04.001
[] Qu XM, Wang RJ, Xin HB, Zhao YY, Fan XT. 2009. Geochronology and geochemistry of igneous rocks related to the subduction of the Tethys oceanic plate along the Bangong Lake arc zone, the western Tibetan Plateau. Geochimica, 38(6): 523–535.
[] Qu XM, Xin HB, Du DD, Chen H. 2012. Ages of post-collisional A-type granite and constraints on the closure of the oceanic basin in the middle segment of the Bangonghu-Nujiang suture, the Tibetan plateau. Geochimica, 41(1): 1–14.
[] Qu XM, Xin HB, Du DD, Chen H. 2013. Magma source of the A-type granite and slab break-off in the middle segment of the Bangonghu-Nujiang suture, Tibet plateau. Acta Geologica Sinica, 87(6): 759–772.
[] Ren JS, Xiao LW. 2004. Lifting the mysterious veil of the tectonics of the Qinghai-Tibet Plateau by 1:250000 geological mapping. Geological Bulletin of China, 23(1): 1–11.
[] Rona PA. 1978. Criteria for recognition of hydrothermal mineral deposits in oceanic crust. Economic Geology, 73(2): 135–160. DOI:10.2113/gsecongeo.73.2.135
[] Shi RD. 2004. The time limits of the Bangong Co-Nujiang Ocean:SSZ-type ophiolite age constraints. Chinese Science Bulletin, 52(2): 223–227.
[] Sillitoe RH. 1972. A plate tectonic model for the origin of porphyry copper deposits. Economic Geology, 67(2): 184–197. DOI:10.2113/gsecongeo.67.2.184
[] Sillitoe RH. 1973. The tops and bottoms of porphyry copper deposits. Economic Geology, 68(6): 799–815. DOI:10.2113/gsecongeo.68.6.799
[] Sillitoe RH. 1983. Enargite-bearing massive sulfide deposits high in porphyry copper systems. Economic Geology, 78(2): 348–352. DOI:10.2113/gsecongeo.78.2.348
[] Sillitoe RH, Steele GB, Thompson JFH, Lang JR. 1998. Advanced argillic lithocaps in the Bolivian tin-silver belt. Mineralium Deposita, 33(6): 539–546. DOI:10.1007/s001260050170
[] Sillitoe RH. 2010. Porphyry copper systems. Economic Geology, 105(1): 3–41. DOI:10.2113/gsecongeo.105.1.3
[] Song B, Zhang YH, Wan YS, Jian P. 2002. Mount making and procedure of the SHRIMP dating. Geological Review, 48(Suppl.): 26–30.
[] Stern CR, Kilian R. 1996. Role of the subducted slab, mantle wedge and continental crust in the generation of adakites from the Andean Austral volcanic zone. Contributions to Mineralogy and Petrology, 123(3): 263–281. DOI:10.1007/s004100050155
[] Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts:Implications for mantle composition and processes. In:Saunders AD and Norry MJ (eds.). Magmatism in the Ocean Basins. Geological Society, London, Special Publication, 42(1):313-345
[] Tang JX, Li FJ, Li ZJ, Zhang L, Tang XQ, Deng Q, Lang XH, Huang Y, Yao XF, Wang Y. 2010. Time limit for formation of main geological bodies in Xiongcun copper-gold deposit, Xietongmen County, Tibet:Evidence from zircon U-Pb ages and Re-Os age of molybdenite. Mineral Deposits, 29(3): 461–475.
[] Tang JX, Sun XG, Ding S, Wang Q, Wang YY, Yang C, Chen HQ, Li YB, Li YB, Wei LJ, Zhang Z, Song JL, Yang HH, Duan JL, Gao K, Fang X, Tan JY. 2014. Discovery of the epithermal deposit of Cu (Au-Ag) in the Duolong ore concentrating area, Tibet. Acta Geoscientica Sinica, 35(1): 6–10.
[] Wang Q, Xu JF, Zhao ZH. 2001. The summary and comment on research on a new kind of igneous rock:Adakite. Advance in Earth Sciences, 16(2): 201–208.
[] Wang Q, McDermott F, Xu JF, Bellon H, Zhu YT. 2005. Cenozoic K-rich adakitic volcanic rocks in the Hohxil area, northern Tibet:Lower-crustal melting in an intracontinental setting. Geology, 33(6): 465–468. DOI:10.1130/G21522.1
[] Wang Q, Wyman DA, Xu JF, Jian P, Zhao ZH, Li CF, Xu W, Ma JL, He B. 2007. Early Cretaceous adakitic granites in the Northern Dabie complex, Central China:Implications for partial melting and delamination of thickened lower crust. Geochimica et Cosmochimica Acta, 71(10): 2609–2636. DOI:10.1016/j.gca.2007.03.008
[] Wang Q, Tang JX, Fang X, Lin B, Song Y, Wang YY, Yang HH, Yang C, Li YB, Wei LJ, Feng J, Li L. 2015. Petrogenetic setting of andsites in Rongna ore block, Tiegelong Cu (Au-Ag) deposit, Duolong ore concentration area, Tibet:Evidence from zircon U-Pb LA-ICP-MS dating and petrogeochemistry of andsites. Geology in China, 42(5): 1324–1336.
[] Wang Y, Zhang Q, Qian Q. 2000. Adakite:Geochemical characteristics and tectonic significances. Scientia Geologica Sinica, 35(2): 251–256.
[] Wang ZH, Wang YS, Xie YH, Sun ZG, Lu ZL, Qu YG, Li CZ, Jiang XF. 2005. The Tarenben oceanic-island basalts in the middle part of the Bangong-Nujiang suture zone, Xizang and their geological implications. Sedimentary Geology and Tethyan Geology, 25(1-2): 155–162.
[] Wright JB. 1969. A simple alkalinity ratio and its application to questions of non-orogenic granite genesis. Geological Magazine, 106(4): 370–384. DOI:10.1017/S0016756800058222
[] Wu FY, Yang YH, Xie LW, Yang JH, Xu P. 2006. Hf isotopic compositions of the standard zircons and baddeleyites used in U-Pb geochronology. Chemical Geology, 234(1-2): 105–126. DOI:10.1016/j.chemgeo.2006.05.003
[] Wu YB, Zheng YF. 2004. Genesis of zircon and its constraints on interpretation of U-Pb age. Chinese Science Bulletin, 49(15): 1554–1569. DOI:10.1007/BF03184122
[] Xu QS, Qin F, Liu Y, Yuan B, Sun H, Chen XF, Zheng J, Niu XL. 2010. Lithocaps:Geological characteristics and implication to exploration of epithermal and porphyry-style deposits. Geology and Exploration, 46(1): 20–23.
[] Xu RK, Zheng YY, Zhao PJ, Shan L, Zhang YL, Cao L, Qi JH, Zhang GY, Dai FH. 2007. Definition and geological significance of the Gacangjian volcanic arc north of Dongqiao, Tibet. Geology in China, 34(5): 768–777.
[] Xu ZQ, Ji SC, Cai ZH, Zeng LS, Geng QR, Cao H. 2012. Kinematics and dynamics of the Namche Barwa Syntaxis, eastern Himalaya:Constraints from deformation, fabrics and geochronology. Gondwana Research, 21(1): 19–36. DOI:10.1016/j.gr.2011.06.010
[] Yang ZS, Hou ZQ, Meng XJ, Liu YC, Fei HC, Tian SH, Li ZQ, Gao W. 2009. Post-collisional Sb and Au mineralization related to the South Tibetan detachment system, Himalayan orogen. Ore Geology Reviews, 36(1-3): 194–212. DOI:10.1016/j.oregeorev.2009.03.005
[] Yin A, Harrison TM. 2000. Geologic evolution of the Himalayan-Tibetan Orogen. Annual Review of Earth and Planetary Sciences, 28: 211–280. DOI:10.1146/annurev.earth.28.1.211
[] Yuan HL, Wu FY, Gao S, Liu XM, Xu P, Sun DY. 2003. Determination of U-Pb age and rare earth element concentrations of zircons from Cenozoic intrusions in northeastern China by laser ablation ICP-MS. Chinese Science Bulletin, 48(22): 2411–2421.
[] Zhang Q, Wang Y, Qian Q, Yang JH, Wang YL, Zhao TP, Guo GJ. 2001. The characteristics and tectonic-metallogenic significances of the adakites in Yanshan period from eastern China. Acta Petrologica Sinica, 17(2): 236–244.
[] Zhao YY, Cui YB, Lü LN, Shi DH. 2011. Chronology, geochemical characteristics and the significance of Shesuo copper polymetallic deposit, Tibet. Acta Petrologica Sinica, 27(7): 2132–2142.
[] Zhou T, Chen C, Liang S, Chen GF, Li HL, Li DW. 2014. Zircon U-Pb geochronology and geochemical characteristics of volcanic rocks in the ophiolite mélange at the Bangong Lake, Tibet. Geotectonica et Metallogenia, 38(1): 157–167.
[] Zhu DC, Pan GT, Mo XX, Wang LQ, Zhao ZD, Liao ZL, Geng QR, Dong GC. 2006. Identification for the Mesozoic OIB-type basalts in central Qinghai-Tibetan Plateau:Geochronology, geochemistry and their tectonic setting. Acta Geologica Sinica, 80(9): 1312–1328.
[] Zhu DC, Zhao ZD, Pan GT, Lee HY, Kang ZQ, Liao ZL, Wang LQ, Li GM, Dong GC, Liu B. 2009. Early cretaceous subduction-related adakite-like rocks of the Gangdese Belt, southern Tibet:Products of slab melting and subsequent melt-peridotite interaction?. Journal of Asian Earth Sciences, 34(3): 298–309. DOI:10.1016/j.jseaes.2008.05.003
[] Zhu DC, Mo XX, Zhao ZD, Niu YL, Wang LQ, Chu QH, Pan GT, Xu JF, Zhou CY. 2010. Presence of Permian extension-and arc-type magmatism in southern Tibet:Paleogeographic implications. Geological Society of America Bulletin, 122(7-8): 979–993. DOI:10.1130/B30062.1
[] Zhu DC, Zhao ZD, Niu YL, Mo XX, Chung SL, Hou ZQ, Wang LQ, Wu FY. 2011. The Lhasa Terrane:Record of a microcontinent and its histories of drift and growth. Earth and Planetary Science Letters, 301(1-2): 241–255. DOI:10.1016/j.epsl.2010.11.005
[] Zhu DC, Zhao ZD, Niu YL, Dilek Y, Hou ZQ, Mo XX. 2013. The origin and pre-Cenozoic evolution of the Tibetan Plateau. Gondwana Research, 23(4): 1429–1454. DOI:10.1016/j.gr.2012.02.002
[] 陈华安, 祝向平, 马东方, 黄瀚霄, 李光明, 李玉彬, 李玉昌, 卫鲁杰, 刘朝强. 2013. 西藏波龙斑岩铜金矿床成矿斑岩年代学、岩石化学特征及其成矿意义. 地质学报, 87(10): 1593–1611.
[] 陈衍景. 2013. 大陆碰撞成矿理论的创建及应用. 岩石学报, 29(1): 1–17.
[] 杜德道, 曲晓明, 王根厚, 辛洪波, 刘治博. 2011. 西藏班公湖-怒江缝合带西段中特提斯洋盆的双向俯冲:来自岛弧型花岗岩锆石U-Pb年龄和元素地球化学的证据. 岩石学报, 27(7): 1993–2002.
[] 冯彩霞. 2011. 班公湖-怒江缝合带西段改则硅质岩地球化学特征及沉积环境. 矿床地质, 30(5): 773–786.
[] 高顺宝, 郑有业, 王进寿, 张众, 杨成. 2011. 西藏班戈地区侵入岩年代学和地球化学:对班公湖-怒江洋盆演化时限的制约. 岩石学报, 27(7): 1973–1982.
[] 耿全如, 潘桂棠, 王立全, 彭智敏, 张璋. 2011. 班公湖-怒江带、羌塘地块特提斯演化与成矿地质背景. 地质通报, 30(8): 1261–1274.
[] 侯可军, 李延河, 田有荣. 2009. LA-MC-ICP-MS锆石微区原位U-Pb定年技术. 矿床地质, 28(4): 481–492.
[] 侯增谦, 高永丰, 孟祥金, 曲晓明, 黄卫. 2004. 西藏冈底斯中新世斑岩铜矿带:埃达克质斑岩成因与构造控制. 岩石学报, 20(2): 239–248.
[] 江军华, 王瑞江, 曲晓明, 辛洪波, 王振中. 2011. 青藏高原西部班公湖岛弧带特提斯洋盆闭合后的地壳伸展作用. 地球科学, 36(6): 1021–1032.
[] 康志强, 许继峰, 王保弟, 董彦辉, 王树庆, 陈建林. 2009. 拉萨地块北部白垩纪多尼组火山岩的地球化学:形成的构造环境. 地球科学, 34(1): 89–104.
[] 康志强, 许继峰, 王保弟, 陈建林. 2010. 拉萨地块北部去申拉组火山岩:班公湖-怒江特提斯洋南向俯冲的产物. 岩石学报, 26(10): 3106–3116.
[] 李光明, 李金祥, 秦克章, 张天平, 肖波. 2007. 西藏班公湖带多不杂超大型富金斑岩铜矿的高温高盐高氧化成矿流体:流体包裹体证据. 岩石学报, 23(5): 935–952.
[] 李光明, 段志明, 刘波, 张晖, 董随亮, 张丽. 2011. 西藏班公湖-怒江结合带北缘多龙地区侏罗纪增生杂岩的特征及意义. 地质通报, 30(8): 1256–1260.
[] 李金祥, 李光明, 秦克章, 肖波. 2008. 班公湖带多不杂富金斑岩铜矿床斑岩-火山岩的地球化学特征与时代:对成矿构造背景的制约. 岩石学报, 24(3): 531–543.
[] 李金祥. 2008.班公湖带多不杂超大型富金斑岩铜矿床的成岩成矿年代学、岩石学及高氧化岩浆-流体-成矿作用.博士学位论文.北京:中国科学院地质与地球物理研究所, 1-224
[] 黎彤. 1992. 地壳元素丰度的若干统计特征. 地质与勘探, 28(10): 1–7.
[] 李玉彬, 多吉, 钟婉婷, 李玉昌, 强巴旺堆, 陈红旗, 刘鸿飞, 张金树, 张天平, 徐志忠, 范安辉, 索朗旺钦. 2012. 西藏改则县多不杂斑岩型铜金矿床勘查模型. 地质与勘探, 48(2): 274–287.
[] 柳小明, 高山, 袁洪林, HattendorfB, GüntherD, 陈亮, 胡圣红. 2002. 193nm LA-ICPMS对国际地质标准参考物质中42种主量和微量元素的分析. 岩石学报, 18(3): 408–418.
[] 潘桂棠, 王立全, 朱弟成. 2004. 青藏高原区域地质调查中几个重大科学问题的思考. 地质通报, 23(1): 12–19.
[] 秦克章, 李光明, 张旗, 李金祥, 缪宇, 肖波, 张天平, 多吉, 李金高, 陆彦. 2006.西藏浅成低温金-银矿的成矿条件与可能产出区分析——从斑岩浅成低温铜金成矿系统的角度.见:第八届全国矿床会议论文集.北京:地质出版社, 660-670
[] 邱瑞照, 周肃, 邓晋福, 李金发, 肖庆辉, 蔡志勇. 2004. 西藏班公湖-怒江西段舍马拉沟蛇绿岩中辉长岩年龄测定——兼论班公湖-怒江蛇绿岩带形成时代. 中国地质, 31(3): 262–268.
[] 曲晓明, 王瑞江, 辛洪波, 赵元艺, 樊兴涛. 2009. 西藏西部与班公湖特提斯洋盆俯冲相关的火成岩年代学和地球化学. 地球化学, 38(6): 523–535.
[] 曲晓明, 辛洪波, 杜德道, 陈华. 2012. 西藏班公湖-怒江缝合带中段碰撞后A型花岗岩的时代及其对洋盆闭合时间的约束. 地球化学, 41(1): 1–14.
[] 曲晓明, 辛洪波, 杜德道, 陈华. 2013. 西藏班公湖-怒江缝合带中段A-型花岗岩的岩浆源区与板片断离. 地质学报, 87(6): 759–772.
[] 任纪舜, 肖黎薇. 2004. 1:25万地质填图进一步揭开了青藏高原大地构造的神秘面纱. 地质通报, 23(1): 1–11.
[] 史仁灯. 2004. 班公湖SSZ型蛇绿岩年龄对班-怒洋时限的制约. 科学通报, 52(2): 223–227.
[] 宋彪, 张玉海, 万渝生, 简平. 2002. 锆石SHRIMP样品靶制作、年龄测定及有关现象讨论. 地质论评, 48(增): 26–30.
[] 唐菊兴, 黎凤佶, 李志军, 张丽, 唐晓倩, 邓起, 郎兴海, 黄勇, 姚晓峰, 王友. 2010. 西藏谢通门县雄村铜金矿主要地质体形成的时限:锆石U-Pb、辉钼矿Re-Os年龄的证据. 矿床地质, 29(3): 461–475.
[] 唐菊兴, 孙兴国, 丁帅, 王勤, 王艺云, 杨超, 陈红旗, 李彦波, 李玉彬, 卫鲁杰, 张志, 宋俊龙, 杨欢欢, 段吉琳, 高轲, 方向, 谭江云. 2014. 西藏多龙矿集区发现浅成低温热液型铜 (金银) 矿床. 地球学报, 35(1): 6–10.
[] 王强, 许继锋, 赵振华. 2001. 一种新的火成岩——埃达克岩的研究综述. 地球科学进展, 16(2): 201–208.
[] 王勤, 唐菊兴, 方向, 林彬, 宋扬, 王艺云, 杨欢欢, 杨超, 李彦波, 卫鲁杰, 冯军, 李力. 2015. 西藏多龙矿集区铁格隆南铜 (金银) 矿床荣那矿段安山岩成岩背景:来自锆石U-Pb年代学、岩石地球化学的证据. 中国地质, 42(5): 1324–1336.
[] 王焰, 张旗, 钱青. 2000. 埃达克岩 (Adakite) 的地球化学特征及其构造意义. 地质科学, 35(2): 251–256.
[] 王忠恒, 王永胜, 谢元和, 孙忠纲, 鲁宗林, 曲永贵, 李存直, 姜雪飞. 2005. 西藏班公湖-怒江缝合带中段塔仁本洋岛型玄武岩的发现及地质意义. 沉积与特提斯, 25(1-2): 155–162.
[] 吴元保, 郑永飞. 2004. 锆石成因矿物学研究及其对U-Pb年龄解释的制约. 科学通报, 49(16): 1589–1602.
[] 徐庆生, 覃锋, 刘阳, 袁波, 孙赫, 陈喜峰, 郑杰, 牛向龙. 2010. 岩帽:地质特征及找矿意义. 地质与勘探, 46(1): 20–23.
[] 许荣科, 郑有业, 赵平甲, 陕亮, 张雨莲, 曹亮, 齐建宏, 张刚阳, 代芳华. 2007. 西藏东巧北尕苍见岛弧的厘定及地质意义. 中国地质, 34(5): 768–777.
[] 袁洪林, 吴福元, 高山, 柳小明, 徐平, 孙德有. 2003. 东北地区新生代侵入体的锆石激光探针U-Pb年龄测定与稀土元素成分分析. 科学通报, 48(14): 1511–1520.
[] 张旗, 王焰, 钱青, 杨进辉, 王元龙, 赵太平, 郭光军. 2001. 中国东部燕山期埃达克岩的特征及其构造-成矿意义. 岩石学报, 17(2): 236–244.
[] 赵元艺, 崔玉斌, 吕立娜, 石登华. 2011. 西藏舍索矽卡岩型铜多金属矿床年代学与地球化学特征及意义. 岩石学报, 27(7): 2132–2142.
[] 周涛, 陈超, 梁桑, 陈桂凡, 李华亮, 李德威. 2014. 西藏班公湖蛇绿混杂岩中火山岩锆石U-Pb年代学及地球化学特征. 大地构造与成矿学, 38(1): 157–167.
[] 朱弟成, 潘桂棠, 莫宣学, 王立全, 赵志丹, 廖忠礼, 耿全如, 董国臣. 2006. 青藏高原中部中生代OIB型玄武岩的识别:年代学、地球化学及其构造环境. 地质学报, 80(9): 1312–1328.