青海玉树赵卡隆铁铜多金属矿床安山岩地球化学特征及其成因

李欢，奚小双，吴城明，张代斌

（1.中南大学 地球科学与信息物理学院，有色金属成矿预测教育部重点实验室，湖南 长沙，410083；2.九州大学工学府地球资源系统，日本 福冈，819-0395；3.广西国土资源规划院，广西 南宁，530022；4.青海诺贝尔矿业有限公司，青海 西宁，810008）

青海省玉树县赵卡隆铁铜多金属矿床是原玉树地质队于1959年发现的，有着长期找矿勘察的历史。但是，由于受到青海西部高原恶劣环境的限制，前人对赵卡隆矿床的地质研究较少。近年来，青海省有色地质矿产勘查局对矿区铁铜铅锌银矿进行了普查，初步认为矿区位于玉树—义敦岛弧带西段的火山沉积盆地中，含矿地层为一套中性火山岩系、碳酸盐岩及碎屑岩系，上部为铅锌，下部为铜，具有喷流沉积的“黑矿”成矿特征，矿床的形成与矿区岛弧安山岩的喷发密切相关，成因类型属火山喷流沉积−改造型铁铜铅锌银多金属矿床[1]。为了研究矿床的大地构造环境、成矿物质来源、控矿因素与成矿规律，本文作者根据矿区安山岩地球化学特征对赵卡隆矿床进行调查研究。岛弧和弧后安山岩与铁、铜等多种金属矿化密切相 关[2−7]，根据安山岩的岩石学[8−9]、同位素地质学[10]、地球化学[11−14]可以探知安山岩的形成背景、演化成因及其与金属矿化之间的关系[15−16]。

1 矿床简介

赵卡隆铁铜多金属矿床大地构造位置处于我国的“三江”成矿带羌塘—昌都地体江达岛弧隆起带北西段之车所—生达地堑中[17]。近年来在“三江”成矿带北段不断发现新矿床，因而备受关注[18−20]。在晚二叠纪，羌塘－昌都地区发生活化，到三叠纪转为太平洋型陆缘，形成青藏川滇印支期沟－弧－盆体系构造环境[21−22]。江达岛弧带是在早二叠纪开始发生差异性抬升、中晚三叠纪发育成熟的[23−24]。区域岩浆岩活动频繁，印支期主要为浅海海底中基性的多旋迴喷发，其岩石主要为安山岩、安山玄武岩类。岩浆岩主要分布于赵卡隆以东、想恩大增一带，呈层状、似层状及透镜状，顺层整合产于上三叠系巴塘群上部碎屑岩组（T3btc）地层内，南部结扎群紫红色碎屑岩组（T3jza）和含煤碎屑岩组（T3jzc）地层中亦有少量安山玄武岩分布。赵卡隆铁铜多金属矿床的成矿物质来源与该期晚期中性火山喷发活动有关。赵卡隆矿区地层简单，出露有上三叠系巴塘群（T3bt）上部碎屑岩组（T3btc）和顶部碳酸盐组（T3btd）及第四系（Q）（图1）。上部碎屑岩组（T3btc）按其岩性组合及段系岩性的含矿性划分为上、中、下3个岩性段。下岩性段（T3btc1）为一套中细粒砂岩、粉砂岩、泥质板岩夹层间石英砾岩、含砾鲕状灰岩组成；中岩性段（T3btc2）由一套粉砂岩、粉砂质板岩互层和中细粒长石石英砂岩组成，除长石石英砂岩中未发现矿体外，其他岩层中均发现有铁铜铅锌银矿矿化；上岩性段（T3btc3）岩性由灰岩、安山岩、泥质板岩、炭质板岩、长石石英砂岩、白云岩等组成，具轻度绿泥石化。赵卡隆铁铜多金属矿床主要产于上岩性段（T3btc3）内。该矿床是以铁、铜为主的多金属矿床，除主元素铁、铜外，还有能综合利用的共（伴）生铅、锌、银、金等有色金属和贵金属元素。矿体严格受地层和岩性控制，含矿带走向长约2 400 m，水平宽度200～800 m，倾斜延深大于700 m。根据矿体分布特点，划分为4个矿化带、7个矿群（图1）。其中以分布于矿区中部的Ⅰ，Ⅱ和Ⅲ号矿群规模最大[25]。

2 安山岩地质产状及矿物组成

图1 赵卡隆矿区地质图（据文献[25]）Fig.1 Geological map of Zhaokalong Deposit（according to Ref.[25]）

安山岩分布于矿区中部，呈层状、带状展布。根据层位不同划分为南、中、北3个岩带。南部安山岩带分布于Ⅵ矿体下部，走向长约 1 200 m，厚度为30～90 m。中部安山岩带规模最大，走向长为1 500 m，厚度为80～200 m，东端分布于Ⅰ矿体下部，与上部岩层呈整合接触，西端分布于Ⅶ矿体顶部，与上部岩组为断层接触。北部安山岩带规模较小，走向长约为1 km，厚度为20～80 m，分布于Ⅳ矿体顶部和Ⅲ矿体底部。成为Ⅳ矿体顶板围岩和Ⅲ矿体底板围岩。矿区各层安山岩岩性特征基本一致，呈灰绿－暗绿色，斑状结构，块状构造。岩石由斑晶和基质组成。斑晶主要为斜长石，少量的角闪石、斜长石聚斑，它们占岩石总量（质量分数，下同）的25%～40%。斑晶的粒径一般为1～2 mm，大者为3 mm。斜长石斑晶形态呈碎裂状、阶梯状、棱角状、少量具熔蚀现象。角闪石斑晶形态呈聚斑状，被绿泥石所交代，保留其外形，往往在斑晶周围形成一圈磁铁矿暗化边。基质主要由斜长石微晶、泥钙铁质和玻璃质组成，为细晶粒状结构，少量具长条状结构，含量达60%～75%。

3 岩石样品的处理及测试方法

岩石样品采自赵卡隆矿区出露的安山岩及部分钻孔岩芯，确保岩石样品新鲜、无风化。对岩石样品进行切片镜下观察之后，挑选出具有代表性的安山岩样品7个，进行送样前预处理。预处理过程为：将岩石样品进行初碎，然后，将初碎后的样品放入Tl−100型碳化钨体碎样机中，细碎至料径为75 μm。每次细碎前，同样要对钵体进行清理，用自来水冲洗干净后，用蒸馏水再次冲洗，接着用干净的纸巾擦拭干，后用蘸有酒精的棉花再次擦洗1遍，最后用吹风机将其吹干，这样确保样品之间不互相污染。将细碎好的样品装入袋中编号，置于烘箱中于105 ℃烘烤2 h，之后将烘干的样品转入干燥器中，以备后期实验送样使用。主量元素、稀土元素和微量元素分析在武汉综合岩矿测试中心（湖北省地质实验研究所）完成，主量元素采用X-Ray Fluorescence（XRF）（RIGAKU2100型）玻璃熔片法进行测定，FeO和Fe2O3含量单独用湿化学法分析进行测定。对主量元素首先测定烧失量，然后加入助熔剂及脱模剂等高温加热，最后将玻璃熔片放入 XRF自动进样系统进行测试。各元素测试结果相对标准偏差大多＜2%，实验结果精度达到5%。稀土及其他微量元素采用ICP-MS分析方法（美国Agileni公司生产的Agileni7500a仪器）分析。把经化学处理、加热、蒸干、加水稀释等加工程序后的样品，用ICPMS溶液法测定微量元素含量，标样为BCR-1，标准参考物测得的相对标准偏差（RSD）大多＜5%，大部分元素（如Rb，Th，U，REE，Nb，Ta，Y，Zr和Hf等）的分析精度高于 5%。稀土元素标准化采用 1982年Taylor球粒陨石标准值。同位素测试工作由武汉地质矿产研究所（原宜昌同位素测试中心）完成，由MAT系列固体及气体同位素质谱仪测定，分析精度为±0.2‰。硫同位素标准品为美国代阿布洛大峡谷（Canyon Diablo）铁陨石中的陨硫铁（CDT），处理后的结果用硫同位素 δ34SCDT表示。赵卡隆矿区安山岩显微照片见图2。所有测试数据及处理结果（质量分数）见表1。

图2 赵卡隆矿区安山岩显微照片Fig.2 Andesite micrographs of Zhaokalong deposit

表1 赵卡隆矿区安山岩主量元素、稀土元素、微量元素及硫同位素分析结果（质量分数）Table 1 Major elements,trace elements,rare earth elements and S isotopic values of andesites from Zhaokalong deposit

4 安山岩地球化学特征

4.1 安山岩主量元素特征

矿区安山岩的 SiO2含量（质量分数，下同）为47.79%～61.51%，平均含量为 55.41%；K2O含量为0.06%～2.90%，平均含量为 1.28%；Na2O含量为2.96%～6.26%，平均含量为 4.81%；全碱（K2O+Na2O）含量平均值为 6.09%。矿区安山岩具有低钾（Na2O含量高于K2O含量），富铝（Al2O3含量为9.59%～17.64%，平均为 15.44%）的特征。CIPW 标准矿物计算结果显示，Q含量为 6.35～21.8%，An含量较小，为0.8%～8.22%，Ab含量为32.92%～53.56%，Or含量为0.35%～18.08%；岩石分异指数（DI）为39.63～72.52。利用CIPW计算结果进行火山岩QAPF图投影（图3），投影结果显示矿区安山岩主要落入英安岩、安山岩区内，反映岩性较为均一。另外，在 Ol′-Q′-Ne′图中，样品全部落入亚碱性系列（图4（a））；在FAM图中，样品主要落入拉斑玄武岩系列（图 4（b））；在 Ab′-An-Or图上投影，样品主要为钠质类型和普通类型（图4（c））。矿区安山岩主要为低钾拉斑玄武岩系列（IAT），少量为钙碱性系列（图 4（b））。DI集中于 35～75，呈单峰式分布。在里特曼−戈蒂里图解中，样品也主要落在了消减带火山岩中（图4（d））。

图3 赵卡隆矿区火山岩Q-A-P-F图解Fig.3 Q-A-P-F diagram of volcanic rocks from Zhaokalong deposit

4.2 安山岩稀土元素特征

矿区安山岩的稀土总量（ΣREE）最低为71.45×10−6，最高为 131.63×10−6；轻稀土元素富集，w（LREE）/w（HREE）为 3.97～9.13，w（La）/w（Yb））N为3.50～8.54，（w（Ce）/w（Yb））N为 2.01～5.97，（w（Sm）/w（Eu））N为1.17～1.44。配分曲线呈右倾式，起伏不大，集中度较高。δEu为 0.92～1.10，Eu无异常；δCe为0.68～0.84，Ce呈弱的负异常，见图5。

4.3 安山岩微量元素特征

矿区安山岩的微量元素的分配型式为强不相容元素（大离子亲石元素）富集型，见图6。数据经标准化处理后，w（RbN）/w（YbN）为 0.44～26.64，平均 11.87；K*含量为0.08～4.70（“*”表示经原始地幔标准化），平均为 2.08；Th*含量为 0.41～32.36，Nb*为 0.21～0.83，Sr*为 0.09～0.89，P*为 0.35～0.63，Zr*为 1.35～2.21，Ti*为0.28～0.64；w（Zr）/w（Nb）为9.32～17.69，w（Hf）/w（Ta）为3.96～5.40。矿区安山岩富集Rb，Ba，Th和Zr等大离子亲石元素，同时亏损 Nb，Sr，P和Ti等高场强元素。蛛网图中自Ce至Lu，各元素的含量下降很快，并且Y和Yb元素有显著亏损，Y含量＜20 μg/g，YN含量为 2～4 μg/g，Yb 含量小于 1.93 μg/g，YbN含量小于 4 μg/g。

4.4 安山岩S同位素特征

火山岩全岩硫稳定同位素分析（表1）表明：δ34SCDT分布范围为−1.60‰～15.32‰；硫同位素比值较集中，除2个样品硫同位素含量比值低于10‰之外，其他高于10‰，平均值为10.53‰。

5 讨论

江达岛弧在区域上属于义敦岛弧带。义敦岛弧带是一个发育于张性薄陆壳基底上的三叠纪岛弧，以发育岛弧裂谷和成对的火山弧为特征[26]。江达构造带源自洋−陆作用的岛弧体制，从晚三叠世以前起经历了陆−陆碰撞焊合、陆内裂谷、陆内造山等多种构造体制的转换，主要的成矿作用发生在陆内及弧后裂谷期[27]。

图4 赵卡隆矿区安山岩地球化学图解Fig.4 Lithogeochemical diagrams of andesites from Zhaokalong deposit

图5 安山岩稀土元素球粒陨石标准化配分型式图Fig.5 Rare-earth elements distribution patterns of andesites

图6 安山岩微量元素原始地幔标准化蛛网图Fig.6 Spider diagram of trace elements from andesites

火山岩的系列、类型是分类命名、划分岩区、分析组合及研究岩浆成因来源、演化机理、构造环境、矿产分布的基础；而火山岩系列、类型的确定，主要依据其化学成分。从主元素分析数据得知，赵卡隆矿床安山岩主成分以贫硅、高钠、低钾含量为特征，属于亚碱性−拉斑玄武岩系列，具有岛弧火山岩性质。火山岩QAPF图显示其属于英安岩类，说明其构造位置可能偏于岛弧中的内弧火山活动[28−29]。矿区安山岩主要为低钾拉斑玄武岩系列（IAT），少量为钙碱性系列（图 4（b）），DI集中于 35～75，呈单峰式分布，反映了构造环境为板块边缘的消减带。在里特曼−戈蒂里图解中，样品也主要落在了消减带火山岩中（图 4（d））。消减带包括岛弧、活动大陆边缘（活动陆缘）2种环境。而与活动陆缘不同的是，岛弧火山岩 SiO2含量为50%～66%（大多＜56%），以低钾（K2O 含量平均为1.60%）安山岩为主[30−31]。因此，从主量元素上看，矿区安山岩主要表现为岛弧火山岩的性质。

从稀土元素特征值来看，赵卡隆安山岩富集轻稀土元素，δEu无明显异常，可能与岩石中斜长石分离结晶较弱有关，表明赵卡隆安山岩与玄武岩浆的分离结晶作用无关。稀土元素右倾曲线表明岩浆来源较为复杂，造成了轻、重稀土元素分异。此外，Ce呈弱的负异常，具有岛弧安山岩的特征，说明喷发环境为海相。与大陆型安山岩不同相比，大洋型安山岩中的稀土含量较低，w（ΣREE）为（25～178）×10−6（大陆型为（67～341）×10−6）[32]，矿区安山岩的各项稀土含量说明其为大洋型安山岩。

在微量元素方面，矿区安山岩富集 Rb，Ba，Th和Zr等大离子亲石元素，同时亏损Nb，Sr，P和Ti等高场强元素，具有与俯冲作用有关的火山岩特征，表现出岛弧火山岩的相关特性[33−34]，显示出俯冲板片流体在岩浆生成中的重要作用。而K的富集以及Nb，P和Ti的亏损，表明物源以壳源为主。同时，微量元素蛛网图中自Ce至Lu，各元素的含量下降很快，表现出裂谷火山岩的特征[35−36]。高场强元素的比值（w（Zr）/w（Nb），w（Hf）/w（Ta））具有类似MORB特性，并且Y和Yb元素有显著的亏损，Y含量＜20 μg/g，YN含量为2～4，Yb含量＜1.93 μg/g，YbN含量小于4 μg/g。微量元素以富集大离子亲石元素和亏损高场强元素为特征，显示其具有热水喷流沉积物的特征[37]，反映了其活动大陆边缘的性质[38]。基于 Th，Yb，Ta及 La分析结果，作出 w（Th）/w（Yb）−w（Ta）/w（Yb）图（图 7（a））和 w（La）/w（Yb）−w（Th）图（图 7（b））。在 w（Th）/w（Yb）−w（Ta）/w（Yb）图中，赵卡隆安山岩样品投点全部落于活动大陆边缘安山岩区；在w（La）/w（Yb）−w（Th）图中，样品主要落在大陆边缘弧和安山弧区域。w（Th）/w（Yb）−w（Ta）/w（Yb）图和 w（La）/w（Yb）−w（Th）图说明赵卡隆安山岩具有大陆边缘安山岩的性质，同时保留了部分岛弧安山岩的特征。此外，岩浆混染了大陆壳物质时，使Nb*含量小于1。而K*可用于指示大洋消减作用与岩浆成分的关系。K*含量大于 1 μg/g，表明其分布于岛弧区，且成分收消减作用的影响。其原因是消减作用能携带活动性元素K进入地幔楔形区并交代这种源区物质[39]。赵卡隆安山岩Nb*和K*的特征说明其形成与岛弧的一侧俯冲有关。

图7 w（Th）/w（Yb）−w（Ta）/w（Yb）[40]和 w（La）/w（Yb）−w（Th）[41]图解Fig.7 w（Th）/w（Yb）−w（Ta）/w（Yb）and w（La）/w（Yb）−w（Th）diagrams

S同位素对判别火山岩的形成环境有着重要作用[42−43]。赵卡隆安山岩δ34SCDT分布范围与海水中的硫同位素特征值相近，说明安山岩硫来源主要为海水，符合深海海底喷发模式。全岩硫稳定同位素分析数据佐证了赵卡隆安山岩为深海相喷发成因。

通过以上对安山岩主量元素、微量元素、稀土元素及硫同位素的研究得出：赵卡隆安山岩既有岛弧安山岩的性质，同时具有活动大陆边缘弧后盆地火山岩的特性。赵卡隆安山岩的成因应该比较复杂，一方面与大洋消减作用岛弧的形成过程有关，另一方面与岛弧形成后的弧后扩张有关。岛弧火山岩在岛弧演化的不同阶段具有不同性质。江达构造带晚石炭世开始岛弧体制的演化,中三叠世末岛弧体制终结，晚三叠世开始陆内弧后扩张活动，赵卡隆安山岩反映了整个岛弧活动过程。安山岩在一定程度上表现出继承了岛弧安山岩的性质，而实应为弧后扩张盆地的大地构造背景，为壳幔2种岩浆混合的产物。赵卡隆矿床岛弧火山岩的产生，一方面与江达岛弧东侧的大洋消减板块俯冲有关，另一方面与岛弧西侧的弧后扩张有关，其形成过程与岛弧的演化过程是一致的。

6 结论

（1）赵卡隆安山岩主量元素以低钾、富铝含量为特征，大多属于英安岩类、拉斑玄武岩系列；稀土元素特征表现为轻稀土元素富集，配分曲线呈右倾式，Eu含量无异常，Ce含量呈弱的负异常；安山岩富集Rb，Ba，Th和Zr等大离子亲石微量元素，同时亏损Nb，Sr，P和 Ti等高场强元素，δ34SCDT平均值为10.53‰。

（2）赵卡隆安山岩既有活动大陆边缘弧后盆地火山岩性质，同时具有岛弧安山岩的特性。赵卡隆安山岩的成因应该比较复杂，一方面与大洋消减作用岛弧的形成过程有关，另一方面与岛弧形成后的弧后扩张作用有关。

（3）江达构造带于晚石炭世开始岛弧体制的演化，于中三叠世末岛弧体制终结，于晚三叠世开始弧后扩张活动，赵卡隆安山岩反映了整个岛弧活动过程。安山岩在一定程度上表现出继承了岛弧安山岩的性质，而实应为弧后扩张盆地的大地构造背景，为壳幔2种岩浆混合的产物，其形成过程与岛弧的演化过程是一致的。

[1]杨武德,张代斌,曹守林.青海省玉树县赵卡隆铁铜多金属矿详查报告[R].西宁:青海省有色地质矿产勘查局七队,青海诺贝尔矿业有限公司,2009.YANG Wu-de,ZHANG Dai-bin,CAO Shou-lin.Detailed investigation report of Zhaokalong Fe-Cu polymetallic deposit,Yushu,Qinghai[R].Xining:Team 7,Geological Exploration Bureau for Nonferrous Metals of Qinghai Province,Qinghai Nobel Mining Ltd,2009.

[2]Donna M,Christopher J V.The Lerokis and Kali Kuning submarine exhalative gold-silver-barite deposits,Wetar Island,Maluku,Indonesia[J].Journal of Geochemical Exploration,1994,50（1/2/3）:351−370.

[3]Prokin V A,Buslaev F P.Massive copper-zinc sulphide deposits in the Urals[J].Ore Geology Reviews,1998,14（1）:1−6.

[4]彭省临,陈子龙.论华南古断拉谷与多因复成铜多金属矿床的关系[J].中南工业大学学报:自然科学版,1991,22（5）:493−499.PENG Sheng-lin,CHEN Zi-long.On the relationship between the paleo-aulacogens and polygenetic compound copper polymetallic ore deposits in south China[J].Journal of Central South University of Technology:Natural Science,1991,22（5）:493−499.

[5]Wright I C,Ronde C E J,Faure K,et al.Discovery of hydrothermal sulfide mineralization from southern Kermadec arc volcanoes（SW Pacific）[J].Earth and Planetary Science Letters,1998,164（1/2）:335−343.

[6]Sorin C,Marin S.Petrographic and geochemical interpretation of the Late Cretaceous volcaniclastic deposits from the Hateg Basin[J].Palaeogeography,Palaeoclimatology,Palaeoecology,2010,293（3/4）:306−318.

[7]Cole J W.Andesites of the Tongariro volcanic centre,North Island,New Zealand[J].Journal of Volcanology and Geothermal Research,1978,3（1/2）:121−153.

[8]Escuder J,Contreras F,Stein G,et al.Magmatic relationships and ages between adakites,magnesian andesites and Nb-enriched basalt-andesites from Hispaniola:Record of a major change in the Caribbean island arc magma sources[J].Lithos,2007,99（3/4）:151−177.

[9]Mathijs A,Ian E M,Jeffrey L,et al.Evolving volcanism at the tip of a propagating arc:The earliest high-Mg andesites in northern New Zealand[J].Journal of Volcanology and Geothermal Research,2010,195（2/3/4）:83−96.

[10]Thirlwall M F,Graham A M,Arculus R J,et al.Resolution of the effects of crustal assimilation,sediment subduction,and fluid transport in island arc magmas:Pb-Sr-Nd-O isotope geochemistry of Grenada,Lesser Antilles[J].Geochimica et Cosmochimica Acta,1996,60（23）:4785−4810.

[11]Ayer J A,Davis D W.Neoarchean evolution of differing convergent margin assemblages in the Wabigoon Subprovince:geochemical and geochronological evidence from the Lake of the Woods greenstone belt,Superior Province,Northwestern Ontario[J].Precambrian Research,1997,81（3/4）:155−178.

[12]James F,Dhanapati H.Barren Island Volcano（NE Indian Ocean）:Island-arc high-alumina basalts produced by troctolite contamination[J].Journal of Volcanology and Geothermal Research,2006,149（3/4）:177−212.

[13]ZHANG Zhao-chong,MAO Jing-wen,CAI Jin-hong,et al.Geochemistry of picrites and associated lavas of a Devonian island arc in the northern Junggar terrane,Xinjiang（NW China）:Implications for petrogenesis,arc mantle sources and tectonic setting[J].Lithos,2008,105（3/4）:379−395.

[14]Andrew J,David R,David P,et al.The nature of magmatism at Palinpinon geothermal field,Negros Island,Philippines:implications for geothermal activity and regional tectonics[J].Journal of Volcanology and Geothermal Research,2004,129（4）:321−342.

[15]Babaie H A,Ghazi A M,Babaei A,et al.Geochemistry of arc volcanic rocks of the Zagros Crush Zone,Neyriz,Iran[J].Journal of Asian Earth Sciences,2001,19（1/2）:61−76.

[16]Dampare S B,Shibata T,Asiedu D K,et al.Geochemistry of Paleoproterozoic metavolcanic rocks from the southern Ashanti volcanic belt,Ghana:petrogenetic and tectonic setting implications[J].Precambrian Research,2008,162（3/4）:403−423.

[17]陈建平,唐菊兴,付小芳,等.西南三江中段成矿规律与成矿预测研究[M].北京:地质出版社,2008:70−76.CHEN Jian-ping,TANG Ju-xing,FU Xiao-fang,et al.Study on metallogenic regularrity and metallogenic prediction of middle part of Sanjiang River,Southwestern China[M].Beijing:Geological Press,2008:70−76.

[18]王富春,王彦,李宝强,等.三江成矿带北西段主要矿床特征研究[J].西北地质,2006,39（2）:179−194.WANG Fu-chun,WANG Yan,LI Bao-qiang,et al.Characteristics of major ore deposits in the northwestern segment of the Sanjian metallogenic belt,southwest China[J].Northwestern Geology,2006,39（2）:179−194.

[19]杨钻云,施泽进,钟康慧,等.青海南部三江成矿带北缘治多地区铁铜铅锌多金属成矿带矿产资源潜力评价[J].地质与勘探,2010,46（3）:426−433.YANG Zuan-yun,SHI Ze-jin,ZHONG Kang-hui,et al.Potential evaluation of multimetal ore resource in the Zhidoi area,northern margin of the Sanjiang metallogenetic belt,southern Qinghai province[J].Geology and Exploration,2010,46（3）:426−433.

[20]董庆吉,肖克炎,陈建平,等.西南“三江”成矿带北段区域成矿断裂信息定量化分析[J].地质通报,2010,29（10）:1479−1485.DONG Qing-ji,XIAO Ke-yan,CHEN Jian-ping,et al.The quantitative analysis of regional metallogenic fault in the northern segment of the Sanjiang metallogenic belt,southwestern China[J].Geological Bulletin of China,2010,29（10）:1479−1485.

[21]罗建宁.羌塘—昌都陆块主要地质构造事件分析[C]//中国地质学会.青藏高原及邻区地质与资源环境学术讨论会论文摘要汇编.成都,2003:74−75.LUO Jian-ling.Analysis of major tectonic events in Qiangtang—Changdu area[C]//Geological Society of China.Abstract compilation of Geological Resources and Environment Symposium,Qinghai—Tibet Plateau and its adjacent areas.Chengdu,2003:74−75.

[22]潘桂堂,徐强.西南“三江”多岛弧造山过程成矿系统与资源评价[M].北京:地质出版社,2003:30−40.PAN Gui-tang,XU Qiang.Orogenic ore-forming system and the resource evaluation process of southwest “Sanjiang”multi-arc[M].Beijing:Geological Press,2003:30−40.

[23]Shellnutt J G,Zhou M F.Permian peralkaline,peraluminou and metaluminous A-type granites in the Panxi district,SW China:Their relationship to the Emeishan mantle plume[J].Chemical Geology,2007,243:286−316.

[24]杨夕辉.三江成矿带中甸地区斑岩型铜矿成矿综合信息研究及资源预测[D].昆明:昆明理工大学国土资源工程学院,2005:13−25.YANG Xi-hui.Ore forming information and resource projections of porphyry copper deposits in Zhongdian area,Sanjiang metallogenic belt[D].Kunming:Kunming University of Science and Technology.Faculty of Land Resource Engineering,2005:13−25.

[25]李欢,奚小双,吴城明,等.青海玉树赵卡隆铁铜多金属矿床地质特征及成因探讨[J].地质与勘探,2011,47（3）:380−387.LI Huan,XI Xiao-shuang,WU Cheng-ming,et al.Geology and metallogenesis of the Zhaokalong iron-copper-polymetallic ore deposit in Yushu,Qinghai Province[J].Geology and Exploration,2011,47（3）:380−387.

[26]侯增谦,莫宣学.义敦岛弧的形成演化及对“三江”地区块状硫化物矿床的控制作用[J].地球科学:中国地质大学学报,1991,16（2）:153−164.HOU Zeng-qian,MO Xuan-xue.The evolution of yidun island-arc and implications in the exploration of kuroko-type volcanogenic assive sulphide deposits in sanjiang area[J].China Earth Science:Journal of China University of Geosciences,1991,16（2）:153−164.

[27]周慧文,钟康惠,唐菊兴,等.藏东江达构造带形成演化与成矿[J].四川地质学报,2008,28（4）:276−282.ZHOU Hui-wen,ZHONG Kang-hui,TANG Ju-xing,et al.Evolution and metallogeny of the Jomda tectonic zone,in east Tibet[J].Acta Geologica Sichuan,2008,28（4）:276−282.

[28]Bayaraa B,Tsuchiya N,Bignall G.Magmatism of the shuteen complex and carboniferous subduction of the Gurvansaikhan terrane,South Mongolia[J].Journal of Asian Earth Sciences,2010,37（5/6）:399−411.

[29]Ian E M,Timothy J,Richard C,et al.Petrogenesis of dacite in an oceanic subduction environment:Raoul Island,Kermadec arc[J].Journal of Volcanology and Geothermal Research,2006,156（3/4）:252−265.

[30]Hiroaki I,Kaori D.Geochemical indications of provenance change as recorded in Miocene shales:Opening of the Japan Sea,San’in region,southwest Japan[J].Marine Geology,1997,144（1/2/3）:211−228.

[31]Michael B,Ulrich K.Significance of slab-derived partial melts and aqueous fluids for the genesis of tholeiitic and calc-alkaline island-arc basalts:Evidence from Mt.Arayat,Philippines[J].Chemical Geology,1993,105（4/5）:233−251.

[32]NIU Yao-ling,Michael J.MORB mantle hosts the missing Eu（Sr,Nb,Ta and Ti）in the continental crust:New perspectives on crustal growth,crust-mantle differentiation and chemical structure of oceanic upper mantle[J].Lithos,2009,112（1/2）:1−17.

[33]Pearce J A.Trace element characteristics of lavas from destructive plate boundaries[C]//Thrope R S.Andesites:Orogenic andesites and related rocks.New York:John Wiley and Sons,1982:525−548.

[34]Dostal J,Dupuy C,Carron J P,et al.Partition coefficients of trace elements:Application to volcanic rocks of St Vincent,West Indies[J].Geochimica et Cosmochimica Acta,1983,47（3）:525−533.

[35]赵振华.关于岩石微量元素构造环境判别图解使用的有关问题[J].大地构造与成矿学,2007,31（1）:92−103.ZHAO Zhen-hua.How to use the trace element diagrams to discriminate tectonic settings[J].Geotectonica et Metallogenia,2007,31（1）:92−103.

[36]赵振华.微量元素地球化学原理[M].北京:科学出版社,1997:24−55.ZHAO Zhen-hua.Principles of trace element geochemistry[M].Beijing:Science Press,1997:24−55.

[37]王登红,付小方,应汉龙.四川西部现代热泉沉积物地球化学特征及意义[J].吉林大学学报:地球科学版,2007,37（5）:878−883.WANG Deng-hong,FU Xiao-fang,YING Han-long.Geochemistry and significance of modern hot-spring sinters in western Sichuan[J].Journal of Jilin University:Earth Science Edition,2007,37（5）:878−883.

[38]Polat A,Peter W U,Brian F,et al.Trace element systematics of the Neoarchean Fiskenæsset anorthosite complex and associated meta-volcanic rocks,SW Greenland:Evidence for a magmatic arc origin[J].Precambrian Research,2009,175（1/2/3/4）:87−115.

[39]Randall A,Martin R,William M,et al.Isotopic and trace element constraints on mixing and melting models of marginal basin volcanism,Bransfield Strait,Antarctica[J].Earth and Planetary Science Letters,1992,111（2/3/4）:287−303.

[40]Pearce J A,Hatris N,Tlndle A G.Trace element discrimination diagrams for the tectonic interpretation of granitic rocks[J].Journal of Petrology,1984,25（4）:956−983.

[41]Condie K C.Geochemitry and tectonic setting of early proterozoic supercrustal rocks in the southwestern United States[J].The Journal of Geology,1986,94（6）:845−861.

[42]Hoog J C M,Taylor B E,Bergen M J.Sulfur isotope systematics of basaltic lavas from Indonesia:implications for the sulfur cycle in subduction zones[J].Earth and Planetary Science Letters,2001,189（3/4）:237−252.

[43]Downes H,Thirlwall M F,Trayhorn S C.Miocene subduction-related magmatism in southern Sardinia:Sr-Nd- and oxygen isotopic evidence for mantle source enrichment[J].Journal of Volcanology and Geothermal Research,2001,106（1/2）:1−22.