岩石学报  2013, Vol. 29 Issue (11): 3767-3778   PDF    
西藏邦铺矿区辉绿玢岩成因及对区域构造岩浆演化的指示
赵晓燕1,2, 杨竹森2, 侯增谦3, 郑远川1, 刘英超3, 田世洪2, 付强3, 费凡1     
1. 中国地质大学地球科学与资源学院,北京 100083;
2. 中国地质科学院矿产资源研究所,北京 100037;
3. 中国地质科学院地质研究所,北京 100037
摘要: 邦铺是继驱龙、甲玛之后在冈底斯斑岩成矿带发现的又一大型的斑岩型矿床。该矿床的成矿作用以斑岩Mo-Cu矿化为主,伴生矽卡岩型Pb-Zn矿化。矿区出露岩体面积较广,种类较多,约占矿区总面积的50%。辉绿玢岩位于矿区的北侧,呈独立脉体产出。其SiO2含量变化于48.39%~50.08%之间,(K2O+ Na2O)平均为3.50%,铝饱和指数A/CNK为1.23,属于过铝质系列,MgO=7.92%~9.12%,Mg#=61.72~64.79。稀土总量∑REE变化于44.93×10-6~67.58×10-6之间,轻重稀土间的分馏不明显,富集Sr、Ba、U,强烈亏损Rb、K等元素。相容元素Ni、Cr含量平均值分别为89×10-6和256×10-6,相对于原生玄武岩浆范围较小,说明邦铺辉绿玢岩经历了一定程度的橄榄石、单斜辉石等镁铁质矿物的分离结晶,但程度不高,可以近似认为辉绿玢岩是岩浆源区平衡部分熔融的产物。锆石U-Pb年龄为14.46±0.38Ma,与成矿年龄14.92Ma相近,锆石原位εHf(t)平均值为4.7。邦铺辉绿玢岩岩石地球化学特征与达孜玄武岩具有高度的一致性,表明二者具有相同的源区,即由于印度板片发生断离,软流圈物质上涌并部分熔融而形成。其成岩年龄的差距是印度板片向北俯冲到不同位置的结果。邦铺辉绿玢岩在成岩过程中受到一定程度地壳物质的混染。它的发现可以作为后碰撞期岩浆活动的一个基性端元,从而为限定地幔性质、岩浆活动时期、以及探索地幔物质对区域岩浆活动及对成矿的影响提供了机会。
关键词: 软流圈地幔     辉绿玢岩     邦铺     冈底斯斑岩成矿带     西藏    
Petrogenesis of diabase porphyrite in Bangpu deposit of Tibet and its instructions to the regional tectonic and magmatic evolution
ZHAO XiaoYan1,2, YANG ZhuSen2, HOU ZengQian3, ZHENG YuanChuan1, LIU YingChao3, TIAN ShiHong2, FU Qiang3, FEI Fan1     
1. School of Geosciences and Resources, China University of Geosciences, Beijing 100083, China;
2. Institute of Mineral Resources, CAGS, Beijing 100037, China;
3. Institute of Geology, CAGS, Beijing 100037, China
Abstract: Bangpu is found as a large-scale porphyry deposit after Qulong and Jiama in Gangdese porphyry copper metallogenic belt which contains mainly porphyry Mo-Cu with skarn Pb-Zn accompanying. There are different species of rock masses in Bangpu mining area, which possess about half of the total area. The diabase porphyrite locates in the north of the mining area as independent vein. The contents of SiO2 are between 48.39%~50.08%, the average content of (K2O+Na2O) is 3.50%, A/CNK is 1.23 and it is belonging to the peraluminous series. MgO=7.92%~9.12%, Mg#=61.72~64.79. The REEs vary from 44.93×10-6 to 67.58×10-6 without detense differentiation between LREE and HREE. The average contents of Ni and Cr are 89×10-6 and 256×10-6 respectively which are lower to native basaltic magma, indicating the diabase porphyrite undergoes mafic minerals such as olivine and clinopyroxene fractional crystallization but not strong, so the diabase porphyrite can be considered as product of equilibrium partial melting. The zircon U-Pb age is 14.46±0.38Ma which is similar to the mineralization age. The average εHf(t) value in situ is 4.7. There is a high degree of consistency between the diabase porphyrite and Dazi basalt in geochemical characteristics which showing the same source region they share. That is the upwelling and partial melting of the asthenospheric because of the break off of Indian plate. The gap between their diagenetic ages is the result of continued northward subduction of Indian plate. The diabase porphyrite has been contaminated by crustal material to a certain extent. It is discovered as a magma-mafic-end member in post collision stage which gives us the opportunity to limit the character of mantle-derived magma and its effect to regional magmatism and mineralization.
Key words: Asthenosphere     Diabase porphyrite     Bangpu     Gangdese porphyry copper metallogenic belt     Tibet    
1 引言

青藏高原,这一世界屋脊,不仅被誉为探索地球的天然实验室,而且更是我国矿产资源的重要产地。其中冈底斯斑岩成矿带即驱龙-厅宫-邦铺斑岩铜钼成矿带位于南拉萨地体(图 1a),长达350km,宽约50km(侯增谦等,2001)。众多学者对冈底斯斑岩成矿带中典型矿床进行的详细解剖表明(孟祥金等, 2005, 2006侯增谦等,2007杨志明等, 2008a, b2011杨志明和侯增谦, 2009; 唐菊兴等,2010王焕等,2011罗茂澄等,2012;Yang et al., 2009Hou et al., 2009, 2011; Hou and Cook, 2009),区内矿床主要形成于印度亚洲大陆处于后碰撞伸展的环境,成岩成矿年龄集中于12~18Ma。近来Zheng et al.(2012a)在冈底斯南部明则矿区发现了形成于30Ma的斑岩Mo-Cu成矿作用,该矿床的成矿地质背景、岩浆源区和成矿作用均与中新世冈底斯斑岩铜矿带相似。

图 1 青藏高原造山带成矿区带图(a, 据Zheng et al., 2012a修改)和邦铺矿区地质简图(b, 底图据周雄,2012) Fig. 1 Geological map of metallogenic belt in Tibetan orogenic belt (a, modified after Zheng et al., 2012a) and of Bangpu deposit (b, base map after Zhou, 2012)

系统的岩石地球化学研究业已证明(曲晓明等,2001Hou et al., 2012a, b ),冈底斯成矿带内的含矿斑岩以碱性-钙碱性为主,且均具埃达克质岩属性,并认为岩浆应源自加厚下地壳的部分熔融(侯增谦等, 2003, 2012Qu et al., 2004, 2007)。与含矿岩浆同期的埃达克质岩浆遍布整条冈底斯带(Hou et al., 2012a, bZheng et al., 2012a, b ),下地壳发生如此广泛和大规模的部分熔融,若仅通过地壳内部机制实现是极难理解的,因而对于产于87°以西的埃达克质岩浆,大多数研究者认为同期源自地幔的钾质-超钾质岩浆可能为其形成提供了大量的热能;然而,钾质-超钾质岩浆在87°以东则尚未见有报道(胡文洁,2012孙晨光,2007Zhao et al., 2009),冈底斯东段后碰撞阶段广泛的埃达克质岩浆的热源问题仍未很好解决。

目前冈底斯斑岩铜矿带东段在后碰撞环境下尚缺乏独立基性端元岩浆存在的报道,仅在局部岩体内发现镁铁质包体。Zheng et al.(2012a)在明则含矿斑岩中发现了大量镁铁质包体,认为其源自拉萨地体岩石圈地幔的部分熔融,并为埃达克质岩浆的形成提供了热源。虽然前人针对赋存于中酸性岩浆里的镁铁质包体做过大量研究,然而由于镁铁质包体与寄主岩石通常会发生强烈的化学交换作用,通过包体来限定地幔性质存在一定的局限性。独立的基性岩则不存在这个问题,它受到中酸性岩浆的混染程度小,更好地保存了基性岩浆的原始特征,因而能更准确地反映地幔端元性质,使我们更深刻地了解区域岩浆的形成机制,为探讨深部动力学演化过程提供依据。本文报道的邦铺辉绿玢岩呈独立脉体产出,其发现为限定地幔性质、岩浆活动时期、以及探索地幔物质对区域岩浆活动及成矿的影响提供了机会。

2 矿区地质

邦铺矿床位于西藏特提斯构造域冈底斯念青唐古拉地体中南部。矿区岀露地层较为简单,除第四系外,主要有下二叠统洛巴堆组(P1l)灰岩和古近纪典中组(E1d)火山岩。典中组地层与洛巴堆组地层之间以断层接触(周雄,2009)。矿区侵入岩比较发育,约占矿区总面积的50%以上。岩体多数呈大规模的侵入体产出,也有的呈小型的岩枝、岩脉产出。岩浆岩主要包括:二长花岗斑岩、花岗斑岩、花岗闪长斑岩、黑云母二长花岗岩、闪长玢岩及辉绿玢岩等(图 1b)。Mo-Cu含矿斑岩主要为二长花岗斑岩及闪长玢岩。通过对邦铺矿区11个钻孔的详细编录,发现邦铺斑岩钼铜矿床具有典型斑岩矿床的蚀变分带特征,从中心向外围依次为:黑云母化带、硅化带、绢云母化带以及青磐岩化带,粘土化多呈补丁状分布于硅化带及青磐岩化带中。

辉绿玢岩主要分布在矿区的东北部,并呈岩脉状产出于典中组火山岩中(图 1b),出露面积约0.23km2,未见与其它岩体有穿插关系。岩体总体较为新鲜,未发生明显的热液蚀变,说明辉绿玢岩成岩作用发生于邦铺斑岩成矿作用之后。

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

本文用于测试的样品采自邦铺矿区东北部,经纬度坐标分别为29°53′52.7″N,91°56′12.9″E,具体位置如图 1b所示。共采集辉绿玢岩10块,挑选5块蚀变较弱的进行全岩化学分析,其余的均用于挑选锆石。

辉绿玢岩:风化面呈灰绿-墨绿色,新鲜面呈绿色(图 2a),斑状结构,块状构造,主要组成矿物为辉石和斜长石。斑晶以斜长石为主,约占10%。基质主要为辉石(45%)、斜长石(40%),以及少量隐晶质矿物(5%)。在显微镜下,板状的斜长石及板柱状的辉石较为杂乱地分布(图 2b),仅在局部见典型的辉绿结构。部分样品有轻微的蚀变,蚀变矿物主要为阳起石、绿泥石、方解石。

图 2 邦铺辉绿玢岩岩石学特征 Pl-斜长石;Pyr-辉石 Fig. 2 Petrological characteristics of the diabase porphyrite in Bangpu deposit Pl-plagioclase; Pyr-pyroxene
3.2 全岩主微量分析

测试的岩石样品通过手标本及镜下鉴定,选取较新鲜无污染的样品碎至200目,送至国家地质实验测试中心进行主微量元素测定。主量元素的分析是在荷兰帕纳科公司Axios波长色散X射线荧光光谱仪(XRF)上完成,微量元素的分析则在ICP-MS上完成。

3.3 锆石U-Pb定年

锆石样品的分离和挑选由廊坊市地岩矿物分选有限公司完成。锆石制靶及阴极发光图像由北京锆石领航科技有限公司完成。将待测锆石置于玻璃板上用无色透明环氧树脂固定,待环氧树脂充分固化后,将靶上的锆石粗磨、细磨至约一半使锆石中心部位暴露并抛光;样品经反射光和透射光照相后,用阴极发光(CL)扫描电镜进行图像分析,然后清洗镀金待测。根据阴极发光图像分析确定单颗粒锆石的成因类型,选择待测点。测定时选择颗粒较大、较自形、清晰的锆石,尽量避开裂纹和包裹体。

LA-MC-ICP-MS锆石U-Pb定年测试分析在中国地质科学院矿产资源研究所MC-ICP-MS实验室完成,锆石定年分析所用仪器为Finnigan Neptune型MC-ICP-MS及与之配套的Newwave UP 213激光剥蚀系统。激光剥蚀所用束斑直径为25μm,频率为10Hz,能量密度约为2.5J/cm2,以He为载气。数据处理采用ICPMSDataCal程序(Liu et al., 2010),测量过程中绝大多数分析点206Pb/204Pb>1000, 未进行普通铅校正,204Pb由离子计数器检测,204Pb含量异常高的分析点可能受包体等普通Pb的影响,对204Pb含量异常高的分析点在计算时剔除,锆石年龄谐和图用Isoplot 3.0程序获得。详细实验测试过程可参见侯可军等(2009) 。样品分析过程中,Plesovice标样作为未知样品的分析结果为336.5±1.1Ma(n=3, 2σ),对应的年龄推荐值为337.13±0.37(2σ)(Sláma et al., 2008),两者在误差范围内完全一致。

3.4 锆石原位Hf同位素分析

锆石Hf同位素测试同样是在中国地质科学院矿产资源研究所MC-ICP-MS实验室进行。所用仪器为Finnigan Neptune型MC-ICP-MS及与之配套的Newwave UP 213激光剥蚀系统,分析时激光束斑直径为55μm,分析点为已经打过U-Pb年龄并且粒径较大的锆石,相关分析流程见侯可军等(2007) 。Hf亏损地幔模式年龄(tDM)的计算采用现今的亏损地幔176Hf/177Hf=0.28325和176Lu/177Hf=0.0384(Vervoort and Blichert-Toft, 1999)。

4 分析测试结果 4.1 常量元素组成

邦铺辉绿玢岩的主微量元素分析结果如表 1所示。在辉绿玢岩五个样品中,BP12-3-1-2烧失量大于6%,因而常量元素结果不予考虑。其余四块样品SiO2含量变化于48.39%~50.08%之间,(K2O+Na2O)为3.27%~3.82%,平均为3.50%;Al2O3含量为14.89%~15.46%,平均为15.16%。MgO的含量为7.92%~9.12%,Mg#值为61.72~64.79,TiO2变化于0.79%~1.14%之间。在TAS全碱图(图 3a)和不活泼元素比值的Zr/TiO2-Nb/Y图解(图 3b)中,邦铺辉绿玢岩样品均落于亚碱性系列玄武岩区域内。

表 1 邦铺辉绿玢岩主量元素(wt%)及微量元素(×10-6)分析结果 Table 1 Major (wt%) and trace element (×10-6) analyze results of the diabase porphyrite in Bangpu deposit

图 3 邦铺辉绿玢岩SiO2-(K2O+Na2O)图解(a, 底图据Irvine and Baragar, 1971)和Nb/Y-Zr/TiO2图解(b, 底图据Cox et al., 1979) Fig. 3 SiO2-(K2O+Na2O) diagram (a, after Irvine and Baragar, 1971) and Nb/Y-Zr/TiO2 diagram (b, after Cox et al., 1979) of the diabase porphyrite in Bangpu deposit
4.2 微量元素组成

邦铺辉绿玢岩稀土总量∑REE较为均一,其值为44.93×10-6~67.58×10-6,平均53.24×10-6。轻重稀土分馏不明显,LREE/HREE=2.79~3.12,(La/Yb)N=2.03~2.33,显示为弱右倾型的稀土配分模式(图 4a)。δEu变化范围为1.00~1.12,显示微弱的正Eu异常。在原始地幔标准化的微量元素比值蛛网图(图 4b)上,邦铺辉绿玢岩富集Sr、Ba、U,强烈亏损Rb、K等元素。相容元素Ni、Cr含量平均值分别为89×10-6和256×10-6,相对于原生玄武岩浆范围(Ni=300×10-6~400×10-6,Cr=300×10-6~500×10-6Frey et al., 1978)较低。

图 4 邦铺辉绿玢岩及达孜玄武岩、驱龙闪长质包体、赛利普超钾质岩和查加寺钾质岩平均值的稀土元素球粒陨石标准化图解(a)和微量元素原始地幔标准化图解(b)(球粒陨石和原始地幔值据Sun and McDonough, 1989) 达孜玄武岩数据引自Gao et al. (2008) ;驱龙闪长质包体数据引自杨志明 (2008) ;赛利普数据孙晨光等(2007) ;查加寺数据胡文洁等(2012) Fig. 4 Chondrite-normalized rare earth element patterns (a) and primitive-mantle-normalized trace element patterns (b) of the Bangpu diabase porphyrite and of the average values of the Dazi basalt, Qulong diorite enclave, Sailipu ultrapotassic rock and Chajiasi potassic rock (chondrite and primitive mantle values after Sun and McDonough, 1989) Data of Dazi basalt from Gao et al. (2008) ; Qulong diorite enclaves from Yang (2008) ; Sailipu from Sun et al. (2007) ; Chajiasi from Hu et al. (2012)
4.3 锆石U-Pb定年

用于测试的锆石多数呈长柱状或短柱状,长宽比介于1:1~2:1之间,粒径变化范围为40~100μm。锆石自形程度较好,大部分锆石呈现白色或浅灰色,也有少数锆石呈黑色。如阴极发光(CL)图像(图 5a)所示,绝大多数锆石发育密集震荡环带,具有岩浆成因锆石的典型特征(吴元保和郑永飞,2004),还有一些锆石呈黑色,不见震荡环带。

图 5 邦铺辉绿玢岩锆石阴极发光图像(a)和U-Pb年龄和谐图(b) Fig. 5 Cathodoluminescence images (a) and U-Pb concordia diagram (b) of zircons from the diabase porphyrite in Bangpu deposit

邦铺辉绿玢岩共完成了18个锆石点的测试,其中10个点为锆石有效点。岩体LA-ICP-MS锆石U-Pb同位素分析结果见表 2。锆石中U含量变化于304.9×10-6~2435×10-6之间,平均值为1799×10-6;Th则变化于822.6×10-6~3580×10-6之间,平均值为1210×10-6;Th/U比值变化于1.07~3.46之间,远大于0.1,属于岩浆成因锆石(Belousova et al., 2002吴元保和郑永飞,2004)。利用Isoplot3(Ludwig,2003)程序进行了谐和曲线投影和206Pb/238U加权平均年龄的计算。结果表明,在207Pb/235U-206Pb/238U谐和图上(图 5b),所有分析点都集中在一致线及其附近的一个很小区域内,等时线年龄为14.46±0.38Ma(其中MSWD=2.0,n=10),表明辉绿玢岩岩脉形成于新生代,与成矿年龄14.92Ma(周雄,2012)近于一致。

表 2 邦铺矿床辉绿玢岩LA-ICP-MS锆石U-Pb年龄测试数据 Table 2 LA-ICP-MS zircon dating result of the diabase porphyrite in Bangpu deposit
4.4 锆石Hf同位素特征

对邦铺辉绿玢岩中晶形较大的10颗已经进行了U-Pb定年的锆石进行了Lu-Hf同位素分析,分析结果如表 3所示。10个测点的176Lu/177Hf比值基本都小于0.002,显示锆石在形成后具有极少的放射成因Hf的积累,因此所测定的176Hf/177Hf比值可以代表锆石结晶时体系的Hf同位素组成(Amelin et al., 1999)。176Hf/177Hf比值分布于0.282851~0.282960之间,平均值为0.282897,εHf(t)值为3.1~7.0,平均值为4.7,单阶段Hf模式年龄tDM为512Ma,二阶段Hf模式年龄tDMC为797Ma,fLu/Hf平均值为-0.95。

表 3 邦铺矿床辉绿玢岩锆石Lu-Hf同位素原位分析结果 Table 3 LA-ICP-MS in situ zircon Lu-Hf isotopic analytical results of the diabase porphyrite in Bangpu deposit
5 讨论 5.1 成岩年龄及意义

邦铺矿区地表岩体面积出露较广,约占矿区总面积的50%左右。含矿岩体为二长花岗斑岩及闪长玢岩。二长花岗斑岩锆石U-Pb年龄为16.23±0.19Ma(罗茂澄等,2011),闪长玢岩锆石U-Pb年龄为14.96±0.16Ma(Wang et al., 2012),二者与成矿年龄即辉钼矿Re-Os年龄15.32±0.79Ma(孟祥金等,2003)在误差范围内一致。邦铺辉绿玢岩成岩年龄为14.46±0.38Ma,与成矿年龄近于一致,说明辉绿玢岩的成岩作用与邦铺矿床成矿存在着密不可分的联系。邦铺辉绿玢岩年龄的精确测定表明14Ma左右在青藏高原存在有幔源岩浆的活动,并产生了独立的基性端元,这也表明幔源物质对冈底斯斑岩成矿带矿床的成矿发挥了不可替代的作用。

5.2 岩浆作用及其源区探讨

邦铺辉绿玢岩SiO2含量较低,介于48.39%~50.08%之间,Mg#变化于61.72~64.79之间,接近于原生玄武质岩浆的范围。但其Cr、Ni含量相对于原生玄武质岩浆较低,说明邦铺辉绿玢岩经历了一定程度的橄榄石、单斜辉石等镁铁质矿物的分离结晶。岩石稀土元素分配曲线呈平缓的右倾形式,轻重稀土分异不明显,无明显的负Eu异常,表明岩浆分离结晶的程度并不强,可以近似认为辉绿玢岩是岩浆源区平衡部分熔融的产物。所以,我们可以根据辉绿玢岩的岩石地球化学特征来示踪岩浆源区。

在YbN-(La/Yb)N(图 6a)及Y-Sr/Y(图 6b)图解中,邦铺辉绿玢岩五个样品都不具备埃达克岩的属性,而冈底斯斑岩铜矿带的含矿斑岩均位于埃达克岩的区域内(图 6),这都说明辉绿玢岩的源区与中新世同斑岩铜矿成矿相关的埃达克质斑岩不一致,不可能起源于碰撞加厚的新生下地壳。驱龙与邦铺同属于冈底斯斑岩铜矿带,空间距离上相距不到40km,是一超大型的斑岩Cu-Mo矿床。在驱龙出露一些闪长质的包体,其锆石U-Pb年龄为19.5±0.4Ma(杨志明,2008),如此相近的距离及相似的年龄使我们更加去关注它们在成因上的联系,以及是否起源于相同的源区?在微量元素配分图解(图 4b)上,邦铺辉绿玢岩与驱龙闪长质包体均显示出弧岩浆的特征,富集Ba、Sr、U等大离子亲石元素,亏损Nb、Ta、Ti等高场强元素,但驱龙闪长质包体强烈的稀土分异特征(图 4a)与邦铺辉绿玢岩明显不同。驱龙闪长质包体较高的Sr/Y及La/Yb比值,以及相对明显的MREE和HREE分异特征表明发生弱石榴石角闪岩化的新生下地壳是其形成的最理想源区(杨志明,2008)。这也再次证明邦铺辉绿玢岩不可能由新生下地壳部分熔融而形成。

图 6 邦铺辉绿玢岩YbN-(La/Yb)N (a)和Y-Sr/Y (b)图解(底图据侯增谦等,2007) 冈底斯斑岩成矿带含矿斑岩数据引自Hou et al., 2009 Fig. 6 YbN vs. (La/Yb)N (a) and Y vs. Sr/Y (b) diagrams of the diabase porphyrite in Bangpu deposit(after Hou et al., 2007) Data of ore-bearing porphyries in the Gangdese porphyry copper metallogenic belt are from Hou et al., 2009

邦铺辉绿玢岩稀土总量较低,平均值为53.24×10-6,(Ce/Yb)N平均值为2.08,与地球平均值(1.15)接近;(La/Yb)N变化于2.03~2.33之间,平均值为2.24,反映幔源原始岩浆的特征。由于Nb、Zr、Y这三种元素在地幔-玄武质熔浆间的总分配系数近于1,又由于它们的不活动性,不受交代、蚀变或风化作用的影响,因而可以基本上代表源区地幔的成分。在Nb-Zr-Y(图 7a, b)图解中,邦铺辉绿玢岩落在亏损地幔的端元内,显示出亏损地幔的源区特征。邦铺辉绿玢岩属于过铝质岩石,Th/Nb、Ba/Th、Zr/Nb及Ba/La的平均值分别为0.72、75.59、53.89、9.48,显示出壳幔源混合的特征,说明源区岩浆在成岩过程中受到一定程度地壳物质的混染。同时邦铺辉绿玢岩176Hf/177Hf比值平均值为0.282897,并且具有较高的εHf(t)值,其平均值为4.7,与矿区含矿的二长花岗斑岩(εHf(t)平均值为4.5)相似(图 8),接近亏损地幔Hf同位素组成,进一步表明辉绿玢岩具有亏损地幔的源区特征。

图 7 邦铺辉绿玢岩Zr-Nb (a)和Zr-Y (b)图解(底图据赵振华, 1997) Fig. 7 Zr vs. Nb (a) and Zr vs. Y (b) diagrams of the diabase porphyrite in Bangpu deposit (after Zhao, 1997)

图 8 邦铺辉绿玢岩的Hf同位素组成和锆石U-Pb年龄作图(底图据杨志明,2008) Fig. 8 Plot of εHf(t) versus U-Pb ages of zircons from the diabase porphyrite in Bangpu deposit (after Yang, 2008)

为了更好地限定岩浆的源区特征,把邦铺辉绿玢岩与达孜玄武岩、后碰撞期的钾质-超钾质岩(赛利普、查加寺)做以对比研究。达孜玄武岩具有较高的Mg含量,Mg#较高,介于53~77之间,K2O含量低,绝大部分小于0.5%,Al2O3变化于14.38%~17.99%之间(Gao et al., 2008),这些均与邦铺辉绿玢岩高Mg、Fe,低K的特征一致。邦铺辉绿玢岩的稀土配分模式与达孜玄武岩较为一致(图 4a),曲线平缓,轻重稀土分异不明显。微量元素原始地幔标准化曲线相差不大(图 4b),相较于达孜火山岩,邦铺辉绿玢岩Th亏损不明显,Nb、Ta分离,这是地壳物质混染的结果。二者Cr、Ni含量均较高,Pb较富集,具备岛弧玄武岩的特征。现在普遍被接受的观点是,达孜玄武岩源区为软流圈地幔,并受到俯冲沉积物的混染。邦铺辉绿玢岩与达孜玄武岩(40~38Ma)虽然成岩年龄相差很大,但相似的地球化学特征不能否认它们源区的一致性。所以邦铺辉绿玢岩应来自于上涌软流圈地幔的熔融,并受到了一定地壳物质的混染。

现在普遍认为,后碰撞的钾质-超钾质岩多显示富集地幔的特征,其源区可能是交代富集上地幔低度部分熔融的产物(胡文洁等,2012孙晨光等,2007)。在图 4中可见,赛利普和查加寺岩石具有相似的稀土和微量元素配分模式,类似于埃达克岩的成分特征,较强的轻重稀土分异特征及高的Sr/Y比值、低的MgO、Y含量明显不同于邦铺辉绿玢岩的岩石地化特征。钾质-超钾质岩浆绝大部分分布在87°以西的区域,这也证明在后碰撞伸展阶段东西地幔性质的不一致性。

5.3 成岩构造背景及区域构造岩浆演化

利用微量元素进行构造岩浆判别是构造背景分析中的一种有效手段,但构造背景存在多种制约因素,所以我们选取多种元素综合判别成岩的构造环境,希望可以真实还原当时的构造演化过程。在Hf-Th-Ta(图 9a)和Zr-Ti图解(图 9b)中,五个点均落入岛弧拉斑玄武岩的范围内;在Zr-Zr/Y图解(图 9c)中,岩石也全部落在火山弧玄武岩中,结合岩石地球化学特征,说明邦铺辉绿玢岩虽形成于大陆碰撞后伸展的环境,但其源区岩浆地球化学特征与岛弧岩浆岩类似。

图 9 邦铺辉绿玢岩Hf/3-Th-Ta图解(a)、Zr-Ti图解(b)和Zr-Zr/Y图解(c) (底图据赵振华, 1997) Fig. 9 Hf/3-Th-Ta diagram (a), Zr-Ti diagram (b) and Zr-Zr/Y diagram (c) of the diabase porphyrite in Bangpu deposit (after Zhao, 1997)

65Ma之前,新特提斯洋壳板片向欧亚大陆下俯冲,大约在65Ma左右,新特提斯洋结束俯冲,印度大陆与欧亚大陆开始强烈碰撞(图 10a)。65~40Ma之间,新特提斯洋片与印度岩石圈地幔一起向拉萨地体下俯冲,在重力或其它的作用力下,发生俯冲洋壳板片的回卷及印度岩石圈地幔的断离与回返,软流圈物质上涌,带来的热量导致特提斯洋片部分熔融,最终形成至今大面积分布的冈底斯花岗岩基和林子宗火山岩,强烈的碰撞也使地壳大幅度增厚(图 10b)。40~25Ma,俯冲的印度板片发生拆沉(40Ma),软流圈物质上涌,一部分软流圈物质部分熔融直接上升至地表,最终形成达孜火山岩;另一部分则底垫于加厚的岩石圈之下,传递软流圈的热量,使地壳及岩石圈地幔(新生下地壳)的地热梯度增高,最终部分熔融形成冈底斯南亚带的具埃达克性质的岩石、镁铁质的包体及钾质超钾质岩(图 10c)。25~8Ma,印度板片继续向北俯冲,并且推动已经拆沉的岩石圈地幔继续向北迁移,到达了冈底斯中亚带的位置;当上涌的软流圈部分熔融并直接上升至近地表时,形成邦铺辉绿玢岩,上升过程中受到来自地壳物质的混染;而当热的软流圈提供热量,导致岩石圈地幔及新生的下地壳部分熔融时,则形成冈 底斯斑岩铜矿带广泛分布的埃达克质斑岩、驱龙的闪长质包体以及钾质 超钾质岩等(图 10d)。~8Ma,印度板片继续向北俯冲,直至班公湖-怒江缝合带,不再有上涌的软流圈提供热量,后碰撞期的岩浆活动也即终止(图 10e, Zheng et al., 2012a, b )。

图 10 印度-亚洲碰撞造山带构造岩浆演化示意图 Fig. 10 Schematic illustrations of magma-structure evolution of Indian-Asian convergent zone
6 结论

(1) 邦铺辉绿玢岩富Mg、Fe,亏损K,轻重稀土分异不明显。锆石U-Pb年龄为14.46±0.38Ma,εHf(t)的平均值为4.7,接近于亏损地幔Hf同位素组成,表明邦铺辉绿玢岩主要源自具亏损性质地幔的部分熔融。

(2) 邦铺辉绿玢岩属于典型的弧岩浆岩,其岩石地球化学特征与达孜玄武岩具有高度的一致性,说明二者具有相同的源区,即由于印度板片发生断离,软流圈物质上涌并部分熔融而形成。其成岩年龄的差距是印度板片向北俯冲到不同位置的结果。邦铺辉绿玢岩在成岩过程中受到一定程度地壳物质的混染。

(3) 邦铺辉绿玢岩呈独立脉体产出,其较好地保存了岩浆的原始特征,因而岩浆可以作为后碰撞期岩浆演化的基性端元,从而为限定地幔性质、岩浆活动时期、以及探索地幔物质对区域岩浆活动及对成矿的影响提供了机会。

致谢 感谢国家测试中心邓月金老师在主微量测试中提供的帮助。中国地质科学院矿产资源研究所MC-ICP-MS实验室侯可军、郭春丽老师在锆石U-Pb及Hf同位素测试中给予了指导。梁维、张雄、徐玉涛等参与了野外的工作。感谢赵志丹老师和董国臣老师提出的宝贵建议!
参考文献
[] Amelin Y, Lee DC, Halliday AN, Robert TP. 1999. Nature of the Earth's earliest crust from hafnium isotopes in single detrital zircons. Nature, 399(6733): 252–255. DOI:10.1038/20426
[] Belousova E, Griffin WL, O'reilly SY, Fisher N. 2002. Igneous zircon: Trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143(5): 602–622. DOI:10.1007/s00410-002-0364-7
[] Cox KG, Bell JD, Pankhurst RJ. 1979. The Interpretation of Igneous Rocks. London: George, Allen and Unwin: 2-35.
[] Frey FA, Green DH, Roy SD. 1978. Integrated models of basalt petrogenesis: A study of quartz tholeiites to olivine melilitites from south eastern Australia utilizing geochemical and experimental petrological data. Journal of Petrology, 19(3): 463–513. DOI:10.1093/petrology/19.3.463
[] Gao YF, Wei RH, Hou ZQ, Tian SH, Zhao RS. 2008. Eocene high-MgO volcanism in southern Tibet: New constraints for mantle source characteristics and deep processes. Lithos, 105(1-2): 63–72. DOI:10.1016/j.lithos.2008.02.008
[] Hou KJ, Li YH, Zou TR, Qu XM, Shi YR, Xie GQ. 2007. Laser ablation-MC-ICP-MS technique for Hf isotope microanalysis of zircon and its geological applications. Acta Petrologica Sinica, 23(10): 2595–2604.
[] Hou KJ, Li YH, Tian YY. 2009. In situ U-Pb zircon dating using laser ablation-multion couting-ICP-MS. Mineral Deposits, 28(4): 481–492.
[] Hou ZQ, Qu XM, Huang W, Gao YF. 2001. Gangdise porphyry copper metallogenic belt: The possible second “Yulong” copper belt. Geology in China, 28(10): 27–29.
[] Hou ZQ, Mo XX, Gao YF, Qu XM, Meng XJ. 2003. Adakite, a possible host rock for porphyry copper deposits: Case studies of porphyry copper belts in Tibetan Plateau and in northern Chile. Mineral Deposits, 22(1): 1–12.
[] Hou ZQ, Pan XF, Yang ZM, Qu XM. 2007. Porphyry Cu-(Mo-Au) deposits no related to oceanic-slab subduction: Examples from Chinese porphyry deposits in continental settings. Geoscience, 21(2): 332–351.
[] Hou ZQ, Cook NJ. 2009. Metallogenesis of the Tibetan collisional orogen: A review and introduction to the special issue. Ore Geology Reviews, 36(1-3): 2–24. DOI:10.1016/j.oregeorev.2009.05.001
[] 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
[] Hou ZQ, Zhang HR, Pan XF, Yang ZM. 2011. Porphyry Cu (-Mo-Au) deposits related to melting of thickened mafic lower crust: Examples from the eastern Tethyan metallogenic domain. Ore Geology Reviews, 39(1-2): 21–45. DOI:10.1016/j.oregeorev.2010.09.002
[] Hou ZQ, Zheng YC, Yang ZM, Rui ZY, Zhao ZD, Qu XM, Sun QZ. 2012a. Contribution of mantle components within juvenile lower-crust to collisional zone porphyry Cu systems in Tibet. Mineralium Deposita, 48(2): 173–192. DOI:10.1007/S00126-012-0415-6
[] Hou ZQ, Zheng YC, Zeng LS, Gao LE, Huang KX, L iW, Li QY, Liang W, Sun QZ. 2012b. Eocene-Oligocene granitoids of southern Tibet: Constraints on crustal anatexis and tectonic evolution of the Himalayan orogen. Earth and Panletary Science Letters, 349-350: 38–52. DOI:10.1016/j.epsl.2012.06.030
[] Hou ZQ, Zheng YC, Yang ZM, Yang ZS. 2012. Metallogenesis of continental collision setting: Part Ⅰ. Gangdese cenozoic porphyry Cu-Mo systems in Tibet. Mineral Deposits, 31(4): 647–670.
[] Hu WJ, Tian SH, Yang ZS, Zhang ZQ. 2012. Petrogenesis of Miocene Chajiasi potassic rocks in western Lhasa block, Tibetan Plateau: Constraints from lithogeochemistry, Geochronology and Sr-Nd isotopes. Mineral Deposits, 31(4): 813–830.
[] Irvine TN, Baragar WRA. 1971. A guide to the chemical classification of the common volcanic rocks. Canadian Journal of Earth Sciences, 8(5): 523–548. DOI:10.1139/e71-055
[] Liu YS, Gao S, Hu ZC, Gao CG, Zong KQ, Wang DB. 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen: U-Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths. Journal of Petrology, 51(1-2): 537–571. DOI:10.1093/petrology/egp082
[] Ludwig KR. 2003. User's Manual for Isoplot 3.00: A Geochronological Toolkit for Microsoft Excel. Berkeley: Kenneth R. Ludwig
[] Luo MC, Wang LQ, Leng QF, Chen W. 2011. Zircon Hf Isotope and Ce4+/Ce3+ ratio of the monzogranite porphyry and biotite monzonitic granite in Bangpu Mo (Cu) deposit, Tibet. Mineral Deposits, 30(2): 266–278.
[] Luo MC, Mao JW, Wang LQ, Leng QF, Chen W. 2012. Fluid inclusion evidence for magmatic-hydrothermal evolution in the Bangpu porphyry molybdenum-copper deposit. Acta Geoscientica Sinica, 33(4): 471–484.
[] Meng XJ, Hou ZQ, Gao YF, Huang W, Qu XM, Qu WJ. 2003. Development of porphyry copper-molybdenum-lead-zinc ore-forming system in East Gangdese belt, Tibet: Evidence from Re-Os age of molybdenite in Bangpu copper polymetallic deposit. Mineral deposits, 22(3): 246–252.
[] Meng XJ, Hou ZQ, Li ZQ. 2005. Fluid inclusions and ore-forming processes of three porphyry copper deposits in Gangdese belt, Tibet. Mineral Deposits, 24(4): 398–408.
[] Meng XJ, Hou ZQ, Li ZQ. 2006. Sulfur and lead isotope compositions of the Qulong porphyry copper deposit, Tibet: Implications for the sources of plutons and metals in the deposit. Acta Geologica Sinica, 80(4): 554–560.
[] Qu XM, Hou ZQ, Huang W. 2001. Is Gangdese porphyry copper belt the second “Yulong” copper belt. Mineral Deposits, 20(4): 355–366.
[] Qu XM, Hou ZQ, Li YG. 2004. Melt components derived from a subducted slab in late orogenic ore-bearing porphyries in the Gangdese copper belt, southern Tibetan plateau. Lithos, 74(3-4): 131–148. DOI:10.1016/j.lithos.2004.01.003
[] 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
[] Sláma J, Kosler J, Condon DJ, Crowley JL, Gerdes A, Hanchar JM, Horstwood MSA, Morris GA, Nasdala L, Norberg N, Schaltegger U, Schoene B, Tubrett MN, Whitehouse MJ. 2008. Plesovice zircon: A new natural reference material for U-Pb and Hf isotopic microanalysis. Chemical Geology, 249(1-2): 1–35. DOI:10.1016/j.chemgeo.2007.11.005
[] Sun CG, Zhao ZD, Mo XX, Zhu DC, Dong GC, Zhou S, Dong S, Xie GG. 2007. Geochemistry and origin of the Miocene Sailipu ultrapotassic rocks in western Lhasa block, Tibetan Plateau. Acta Petrologica Sinica, 23(11): 2715–2726.
[] Sun SS, Mc Doaongh 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 Publications, 42(1): 313–345. DOI:10.1144/GSL.SP.1989.042.01.19
[] Tang JX, Wang DH, Wang XW, Zhong KH, Ying LJ, Zheng WB, Li FJ, Guo N, Qin ZP, Yao XF, Li L, Wang Y, Tang XQ. 2010. Geological features and metallogenic model of the Jiama copper-polymetallic deposit in Tibet. Acta Geoscientia Sinica, 31(4): 495–506.
[] Vervoort JD, Blichert-Toft J. 1999. Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochimica et Cosmochimica Acta, 63(3-4): 533–556. DOI:10.1016/S0016-7037(98)00274-9
[] Wang H, Tang JX, Ying LJ, Wang LQ, Qin ZP. 2011. Characteristics of the main ore minerals in the Jiama Cu-polymetallic ore deposit, Tibet. Journal of Chengdu University of technology (Science and Technology Edition), 38(1): 103–112.
[] Wang ZH, Liu YL, Liu HF, Guo LS, Zhang JS, Xu KF. 2012. Geochronology and geochemistry of the Bangpu Mo-Cu porphyry ore deposit, Tibet. Ore Geology Reviews, 46: 95–105. DOI:10.1016/j.oregeorev.2012.02.004
[] Wu YB, Zheng YF. 2004. Research of zircon minerageny and its restriction about the interpretation of U-Pb dating. Chinese Science Bulletin, 49(16): 1589–1604.
[] Yang ZM. 2008. The Qulong giant porphyry copper deposit in Tibet: Magmatism and Mineralization. Ph. D. Dissertation. Beijing: Chinese Academy of Geological Sciences: 1–54.
[] Yang ZM, Hou ZQ, Xia DX, Song YC, Li Z. 2008a. Relationship between western porphyry and mineralization in Qulong copper deposit of Tibet and its enlightenment to further exploration. Mineral Deposits, 27(1): 28–36.
[] Yang ZM, Hou ZQ, Song YC, Li ZQ, Xia DX, Pan FC. 2008b. Qulong superlarge porphyry Cu deposit in Tibet: Geology, alteration and mineralization. Mineral Deposits, 27(3): 279–318.
[] Yang ZM, Hou ZQ. 2009. Genesis of giant porphyry Cu deposit at Qulong, Tibet: Constraints from fluid inclusions and H-O isotopes. Acta Geologica Sinica, 83(12): 1838–1859.
[] Yang ZM, Hou ZQ, White NC, Chang ZS, Li ZQ, Song YC. 2009. Geology of the post-collisional porphyry copper-molybdenum deposit at Qulong, Tibet. Ore Geology Reviews, 36(1-3): 133–159. DOI:10.1016/j.oregeorev.2009.03.003
[] Yang ZM, Hou ZQ, Jiang YF, Zhang HR, Song YC. 2011. Sr-Nd-Pb and zircon Hf isotopic constraints on petrogenesis of the Late Jurassic granitic porphyry at Qulong, Tibet. Acta Petrologica Sinica, 27(7): 2003–2010.
[] Zhao ZD, Mo XX, Dilek Y, Niu YL, DePaolo DJ, Robinson P, Zhu DC, Sun CG, Dong GC, Zhou S, Luo ZH, Hou ZQ. 2009. Geochemical and Sr-Nd-Pb-O isotopic compositions of the post-collisional ultrapotassic magmatism in SW Tibet: Petrogenesis and implications for India intra-continental subduction beneath southern Tibet. Lithos, 113(1-2): 190–212. DOI:10.1016/j.lithos.2009.02.004
[] Zhao ZH. 1997. Principles of Trace Element Geochemistry. Beijing: Science Press.
[] Zheng YC, Hou ZQ, Li W, Liang W, Huang KX, Li QY, Sun QZ, Fu Q, Zhang S. 2012a. Petrogenesis and geological implications of the Oligocene Chongmuda-Mingze adakite-like intrusions and their mafic enclaves, southern Tibet. Chicago Journals, 120(6): 647–669.
[] Zheng YC, Hou ZQ, Li QY, Sun QZ, Liang W, Fu Q, Li W, Huang KX. 2012b. Origin of Late Oligocene adakitic intrusives in the southeastern Lhasa terrane: Evidence from in situ zircon U-Pb dating, Hf-O isotopes, and whole-rock geochemistry. Lithos, 148: 296–311. DOI:10.1016/j.lithos.2012.05.026
[] Zhou X. 2009. Study on the fluid inclusions of the Bangpu Mo-Cu polymetallic deposit, Tibet. Master Degree Thesis. Chengdu: Chengdu University of Technology, 1-90 (in Chinese with English summary)
[] Zhou X. 2012. Genesis study of Mo-Cu polymetallic deposit from the Bangpu, Tibet. Ph. D. Dissertation. Chengdu: Chengdu University of Technology, 1-193 (in Chinese with English summary)
[] 侯可军, 李延河, 邹天人, 曲晓明, 石玉若, 谢桂青. 2007. LA-MC-ICP-MS锆石Hf同位素的分析方法及地质应用. 岩石学报, 23(10): 2595–2604.
[] 侯可军, 李延河, 田有荣. 2009. LA-MC-ICP-MS锆石微区原位U-Pb定年技术. 矿床地质, 28(4): 481–492.
[] 侯增谦, 曲晓明, 黄卫, 高永丰. 2001. 冈底斯斑岩铜矿成矿带有望成为西藏第二条“玉龙”铜矿带. 中国地质, 28(10): 27–29.
[] 侯增谦, 莫宣学, 高永丰, 曲晓明, 孟祥金. 2003. 埃达克岩: 斑岩铜矿的一种可能的重要含矿母岩——以西藏和智利斑岩铜矿为例. 矿床地质, 22(1): 1–12.
[] 侯增谦, 潘小菲, 杨志明, 曲晓明. 2007. 初论大陆环境斑岩铜矿. 现代地质, 21(2): 332–351.
[] 侯增谦, 郑远川, 杨志明, 杨竹森. 2012. 大陆碰撞成矿作用: I. 冈底斯新生代斑岩成矿系统. 矿床地质, 31(4): 647–670.
[] 胡文洁, 田世洪, 杨竹森, 张兆卿. 2012. 拉萨地块西段中新世查加寺钾质火山岩岩石成因——岩石地球化学、年代学和Sr-Nd同位素约束. 矿床地质, 31(4): 813–830.
[] 罗茂澄, 王立强, 冷秋峰, 陈伟. 2011. 邦铺钼(铜)矿床二长花岗斑岩、黑云二长花岗岩锆石Hf同位素和Ce4+/Ce3+比值. 矿床地质, 30(2): 266–278.
[] 罗茂澄, 毛景文, 王立强, 冷秋峰, 陈伟. 2012. 西藏邦铺斑岩钼铜矿床岩浆-热液流体演化: 流体包裹体研究. 地球学报, 33(4): 471–484.
[] 孟祥金, 侯增谦, 高永丰, 黄卫, 曲晓明屈文俊. 2003. 西藏冈底斯东段斑岩铜钼铅锌成矿系统的发育时限: 帮浦铜多金属矿床辉钼矿Re-Os年龄证据. 矿床地质, 22(3): 246–252.
[] 孟祥金, 侯增谦, 李振清. 2005. 西藏冈底斯三处斑岩铜矿床流体包裹体及成矿作用研究. 矿床地质, 24(4): 398–408.
[] 孟祥金, 侯增谦, 李振清. 2006. 西藏驱龙斑岩铜矿S、Pb同位素组成: 对含矿斑岩与成矿物质来源的指示. 地质学报, 80(4): 554–560.
[] 曲晓明, 侯増谦, 黄卫. 2001. 冈底斯斑岩铜矿(化)带: 西藏第二条: “玉龙” 铜矿带. 矿床地质, 20(4): 355–366.
[] 孙晨光, 赵志丹, 莫宣学, 朱弟成, 董国臣, 周肃, 董昕, 谢国刚. 2007. 青藏高原拉萨地块西部中新世赛利普超钾质岩石的地球化学与岩石成因. 岩石学报, 23(11): 2715–2726.
[] 唐菊兴, 王登红, 汪雄武, 钟康惠, 应立娟, 郑文宝, 黎枫佶, 郭娜, 秦志鹏, 姚晓峰, 李磊, 王友, 唐晓倩. 2010. 西藏甲玛铜多金属矿矿床地质特征及其矿床模型. 地球学报, 31(4): 495–506.
[] 王焕, 唐菊兴, 应立娟, 王立强, 秦志鹏. 2011. 西藏甲玛铜多金属矿床主要矿石矿物特征. 成都理工大学学报 (自然科学版), 38(1): 103–112.
[] 吴元保, 郑永飞. 2004. 锆石成因矿物学研究及其对U-Pb年龄解释的制约. 科学通报, 49(16): 1589–1604.
[] 杨志明. 2008. 西藏驱龙超大型斑岩铜矿床——岩浆作用与矿床成因. 博士学位论文. 北京: 中国地质科学院, 1-154
[] 杨志明, 侯增谦, 夏代祥, 宋玉财, 李政. 2008a. 西藏驱龙铜矿西部斑岩与成矿关系的厘定: 对矿床未来勘探方向的重要启示. 矿床地质, 27(1): 28–36.
[] 杨志明, 侯增谦, 宋玉财, 李振清, 夏代详, 潘凤雏. 2008b. 西藏驱龙超大型斑岩铜矿床: 地质、蚀变与成矿. 矿床地质, 27(3): 279–318.
[] 杨志明, 侯增谦. 2009. 西藏驱龙超大型斑岩铜矿的成因: 流体包裹体及H、O同位素证据. 地质学报, 83(12): 1838–1859.
[] 杨志明, 侯增谦, 江迎飞, 张洪瑞, 宋玉财. 2011. 西藏驱龙矿区早侏罗世斑岩的Sr-Nd-Pb及锆石Hf同位素研究. 岩石学报, 27(7): 2003–2010.
[] 赵振华. 1997. 微量元素地球化学原理. 北京: 科学出版社.
[] 周雄. 2009. 西藏邦铺钼铜多金属矿床流体包裹体研究. 硕士学位论文. 成都: 成都理工大学, 1-90
[] 周雄. 2012. 西藏邦铺铜钼多金属矿床成因研究. 博士学位论文. 成都: 成都理工大学, 1-193