岩石学报  2016, Vol. 32 Issue (8): 2281-2302   PDF    
扬子西缘乡城-丽江结合带燕山期斑岩Mo多金属矿床成矿系统
刘学龙1, 李文昌1,2, 张娜1, 赖安琦1, 李喆1, 杨富成1     
1. 昆明理工大学, 昆明 650093 ;
2. 云南省地质调查局, 昆明 650051
摘要: 位于扬子陆块西缘与西南三江造山带结合部位的乡城-丽江地区,北起四川乡城经云南格咱南至丽江地区,在印支期斑岩铜矿带上新发现叠加了燕山期斑岩Mo多金属成矿作用,形成乡城-丽江斑岩Mo矿带。受印支期古特提斯洋盆闭合后地壳缩短与加厚的影响,燕山晚期下地壳发生拆沉作用,导致后碰撞型花岗质岩浆侵位,发育了Mo多金属成矿作用,并构成斑岩成矿系统。岩石地球化学特征表明,燕山晚期含矿斑岩具高硅(SiO2=66.29%~79.36%)、高碱(K2O+Na2O=5.07%~9.24%)、富钾(K2O/Na2O=0.71~2.13)的特点,属于高钾钙碱性岩石系列。岩石富集轻稀土元素(LREE),具有负δEu异常;微量元素具富集大离子亲石元素K、U、Th、Rb,亏损高场强元素Nb、Ta、P、Ti的特征;大离子亲石元素的富集和Nb、Ta等高场强元素的亏损,表明形成这些岩体的岩浆主要来自地壳,具有造山带花岗岩的地球化学特征。区内代表性含矿斑岩相似的地球化学组成及分布特征表明,乡城-丽江结合带燕山晚期的成矿斑岩是由同源岩浆分异演化而来。成矿系统的分析表明,燕山晚期主要含矿斑岩成矿物质以壳源为主,但具有少量地幔物质的加入,成岩成矿物质来源具有相似或一致的源区特征。在成矿元素的组合上,由成矿斑岩体向外带表现出W、Mo→W、Mo、Cu→Pb、Zn、Ag的演化趋势。研究表明,本区燕山晚期花岗岩浆的侵入及Mo多金属成矿作用并不是独立或个别的成矿事件,而是纵跨义敦岛弧、甘孜-理塘结合带及扬子西缘的带状成矿活动,属于区内与燕山晚期岩浆侵入作用相关的统一斑岩成矿系统。
关键词: 斑岩     Mo多金属矿床     成矿系统     成矿物质来源     燕山晚期     乡城-丽江结合带    
Metallogenic system of the Yanshanian porphyry Mo polymetallic deposit in the Xiangcheng-Lijiang suture zone, western margin of Yangtze block, SW China
LIU XueLong1, LI WenChang1,2, ZHANG Na1, LAI AnQi1, LI Zhe1, YANG FuCheng1     
1. Kunming University of Science and Technology, Kunming 650093, China ;
2. Yunnan Geological Survey, Kunming 650051, China
Abstract: Xiangcheng-Lijiang suture zone located in combined position of the western margin of Yangtze block and Sanjiang orogenic belt, it north to Xiangcheng of Sichuan, through Geza of Yunnan, south to Lijiang area, in which newly discovered that porphyry copper belt of Late Triassic superimposed porphyry Mo-polymetallic mineralization in Late Cretaceous and formed Xiangcheng-Lijiang porphyry molybdenum belt. In Late Cretaceous, by Tethys basin closed and crustal shortening and thickening, lower crustal delamination and post-collision-type granitic magma intrusion developed Mo-polymetallic mineralization, which formed the porphyry metallogenic system. The lithogeochemical characteristics shown that Late Cretaceous ore-bearing porphyry with the characteristics of high silicon (SiO2=66.29%~79.36%), high alkali (K2O+Na2O=5.07%~9.24%), richen potassium (K2O/Na2O=0.71~2.13), which belongs to the K-high calc-alkali rock series. Rock enriched LREE and LILEs (K, U, Th, Rb) and depleted in HFSEs (Nb, Ta, P, Ti); enriched in LILEs and depleted in HFSEs (Nb, Ta) shown that the magma which formed these rocks are mainly from the crust and has the geochemical characteristics of orogenic granitoids. Similar geochemical composition and distribution characteristics of representative ore-bearing porphyry in this area shown that ore-forming porphyry in Xiangcheng-Lijiang suture zone was derived from the differentiation and evolution of the same source magma in Late Cretaceous. By analysis of metallogenic system, the metallogenic materials of ore-bearing porphyry are mainly derived from the crust, but with a small amount of mantle material added, rock-forming and metallogenic meterials with the similar or consistent source region characteristics. In the combination of ore-forming elements, from the metallogenic porphyry body to the outer zone shown that the evolutional trend of W, Mo→W, Mo, Cu→PbZn, Ag. This study shown the intrusion of the Late Cretaceous granite magma and the Mo-polymetallic mineralization in the area are not independent or individual mineralization events, but spanning the banded metallogenic activities of Yidun island arc and Ganzi-Litang suture zone and the western margin of the Yangtze and belong to a part of the unified porphyry metallogenic system related to magmatic intrusion in Late Cretaceous in the region.
Key words: Porphyry     Mo-polymetallic deposits     Metallogenic system     Source of ore-forming materials     Late Cretaceous     Xiangcheng-Lijiang suture zone    
1 引言

多年来的地质找矿和研究工作表明,乡城-丽江地区Cu-Mo多金属矿床的成矿作用与中酸性侵入岩体存在密切的成因联系,因此众多地质学家对区内典型成矿斑岩体开展了许多富有创新的研究工作。研究表明,区内Cu多金属矿床形成与洋壳俯冲造山作用相关,成岩成矿作用主要发生在晚三叠世甘孜-理塘洋壳西向俯冲背景下(杨岳清等,2002; 李文昌,2007; Li et al.,2011; 李文昌等,2011; 曹殿华等,2009; 庞振山等,2009; Deng et al.,2014)。岩石地球化学特征表明,斑岩体具有与洋壳俯冲成因埃达克岩相似的地球化学性质,主要属钙碱性系列(曾普胜等,2006; 冷成彪等,2007; 庞振山等,2009; 任江波等,2011; 刘学龙等,2012a2013)。近年来,不少研究者对成矿斑岩体开展了多方面研究工作,特别是取得了一些新的测年数据和同位素资料(曾普胜等,2006; 李建康等,2007; 李文昌,2007; 李文昌等,20092010; 杜杨松等,2007; 王守旭等,2008; 曹殿华等,2009; 尹光候等,2009; 黄肖潇等,2012; 刘学龙等,2013; 金灿海等,2013),多元测年方法对区内代表性斑岩体所获得年龄数据存在较大的差异,但主要集中于220~200Ma,指示成岩成矿作用主要发生在印支期。然而,前人的研究工作多涉及印支期研究较为成熟的成矿岩体(普朗、红山、雪鸡坪)的研究。燕山期花岗岩和斑岩Mo(Cu)矿床的研究成果相对较少,主要是获得了一些零星的年代学方面的证据(李建康等,2007; 李文昌等,2012; 黄肖潇等,2012; 彭惠娟等,2014; 余海军等,2015),岩浆活动的时限为96.4~65Ma,Mo多金属成矿作用的时限为86~77Ma(Li et al.,2014; 杨立强等,2015)。进入喜山期以来(65~15Ma),受青藏高原的碰撞隆起的影响,区内表现为逆冲-推覆构造和大规模走滑平移活动,形成了喜山期花岗岩侵位(李文昌,2007)。

杜杨松,庞振山,王功文,郭欣,李青. 2007. 云南省普朗铜矿成矿规律与成矿预测研究. 北京:中国地质大学,20-100

最近,李文昌团队研究工作对区内斑岩Cu多金属矿床成矿系统的研究取得了一定的成果。乡城-丽江结合带成矿作用以三期构造-岩浆活动为基础划分为三套斑岩成矿系统,即印支期斑岩Cu多金属成矿系统、燕山期斑岩Mo-Cu多金属成矿系统和喜山期富碱斑岩Au-Mo-Cu多金属成矿系统(李文昌等,2013)。印支期,由于岩浆作用强度、构造环境的局部差异,导致在空间上形成了不同规模、不同类型的矿床。主要表现在岩体内或外接触带产出斑岩型Cu(Au)矿床(普朗、雪鸡坪等)、矽卡岩型Cu矿床(浪都)及矽卡岩-斑岩复合型Cu(Mo-Pb-Zn)矿床(红山),在斑岩活动的外围形成产于构造角砾岩中的热液型含金富银的Pb-Zn(Cu)矿床(恩卡、亚杂)。印支期形成的斑岩型、矽卡岩型及热液脉型矿床在时间、空间上的共生,构成了一个较为完整的斑岩Cu矿成矿系统(刘学龙等,2012b)。燕山期,该区进入陆内碰撞地壳增厚和后碰撞造山阶段,发生地壳挤压;燕山晚期,软流圈上涌发生地壳熔融,形成了本区另外一条重要的构造-岩浆成矿带(李文昌,2007)。区内燕山晚期的花岗岩体呈近SN向带状叠加于印支期斑(玢)岩带上,该斑岩带向北延至四川乡城、稻城,南延至扬子陆块西缘的香格里拉铜厂沟(拉巴)、东炉房地区,成矿岩体的空间分布从北往南呈现出规模出露-半隐伏-隐伏的分布规律,由粗粒花岗岩-似斑状花岗岩-斑岩的分布规律。燕山晚期,受地壳缩短与加厚的影响,下地壳最终发生拆沉作用,软流圈热上涌,导致地壳重熔和花岗岩浆的侵位,形成以斑岩型Mo为主的斑岩-矽卡岩型矿床(刘学龙等,2013; Li et al.,2014)。喜山期,受青藏高原碰撞隆升和大型走滑作用的影响,发育有幔源岩浆的侵入作用,并形成了斑岩型Au(Cu-Mo)矿床(李文昌,2007; Deng and Wang,2016)。

区内燕山期花岗岩的大面积分布表明,在陆-陆碰撞造山作用过程中也存在构造-岩浆-成矿作用的发生。据目前勘查和研究的工作程度,乡城-丽江结合带燕山晚期的Mo多金属成矿作用的代表性矿床主要有铜厂沟斑岩型Mo-Cu矿床、休瓦促W-Mo-Cu多金属矿床、热林Cu-Mo矿床、红山Cu-Mo矿床。本文在前人研究工作基础上,对区内燕山晚期的构造-岩浆活动及成矿事件进行了系统厘定,建立了成矿系统模式,以期为该区地质找矿工作提供科学参考。

2 区域成矿背景

乡城-丽江结合带位于特提斯-喜马拉雅成矿域义敦岛弧带南段,属于义敦-香格里拉Cu-Au-Ag-Pb-Zn-Sn-Hg-Sb-W-Be成矿带(图 1)。区内褶皱、断裂发育,由一系列NW向紧密线性褶皱和NW-NNW向断裂组成,是控制沉积建造、变质作用、岩浆活动及其矿床产出的主要构造,而次级同向断裂及东西向断层则为容岩(矿)构造,晚期发育规模较小的北东向断层,切错了早期断裂及褶皱(李文昌和刘学龙,2015)。断裂构造早期为NNW向或近SN走向,NNW向断裂多为压性,切穿和改造了SN向断裂;NE向或近EW向断裂形成时间最晚,但对早期构造具改造作用。格咱河断裂是区内主要的断裂构造,呈近SN向展布,控制了印支期、燕山晚期出露的中酸性侵入岩的空间展布。印支期、燕山晚期出露的中酸性斑(玢)岩体仅在格咱河断裂以东研究区,主要呈岩株、岩墙、岩枝、岩脉状分布,并严格受NNW向断裂的控制,且具有带状分布的特点。喜山期走滑剪切和拉张构造环境形成的富碱斑岩,主要分布于格咱河断裂东侧的甭哥地区和三江口-剑川断裂西侧的东炉房等地区,断裂构造控制明显。伴随洋壳的俯冲消减,区内形成了与中酸性岩浆侵入作用密切相关的斑岩型(矽卡岩型)Cu矿床;燕山晚期,受地壳缩短与加厚的影响,下地壳发生拆沉作用,最终导致该区发生强烈的后碰撞型花岗岩浆侵位,并伴有Mo(Cu)多金属成矿作用的形成;喜山期,受青藏高原碰撞隆升和侧向走滑剪切作用的影响,伴随新生地壳的部分熔融发育有富碱斑岩的侵入作用,并形成了斑岩型Au(Cu-Mo)多金属矿床(杨岳清等,2002; 曾普胜等,2003; 李文昌,2007; 李文昌等,2010; 刘学龙等,2012a; Liu et al.,2013; Deng et al.,2015)。

3 区域典型矿床特征

燕山期发育的酸性(斑)岩带呈近南北向展布,纵跨了义敦岛弧带、甘孜-理塘结合带、扬子陆块西缘3个印支期构造单元。本次LA-ICP-MS锆石测年、岩石地球化学特征研究所选岩石样品采自滇西北乡城-丽江结合带燕山期的代表性成矿斑岩体,主要包括铜厂沟花岗闪长斑岩、休瓦促二长花岗岩、花岗闪长岩、热林二长花岗岩、红山深部(隐伏)花岗闪长斑岩(图 1)。采样位置主要为钻孔岩芯内及其新鲜的地表露头,样品质量大于3kg。

图 1 乡城-丽江结合带构造-岩浆-成矿分布图(据云南省地质调查院,2013修编) Ⅰ-扬子板块;Ⅱ-甘孜-理塘板块结合带;Ⅲ-义敦岛弧带;Ⅳ-中咱微陆块;Ⅴ-金沙江结合带;Ⅵ-江达-维西火山弧;Ⅶ-昌都-兰坪陆块;Ⅷ-三达山-景洪火山弧;Ⅸ-澜沧江结合带;Ⅹ-保山地块 Ⅰ-Yangtze plate; Ⅱ-Ganze-Litang suture; Ⅲ-Yidun Island-arc belt; Ⅳ-Zhongza massif; Ⅴ-Jingshajiang suture; Ⅵ-Jiangda-Weixi volcanic; Ⅶ-Changdu-Langping block; Ⅷ-Sandashan-Jinghong volcanic arc; Ⅸ-Langcangjiang suture; Ⅹ-Baoshan block Fig. 1 Tectonic-magmatic-metallogenic system of Xiangcheng-Lijiang suture zone

① 云南省地质调查院.2013.云南香格里拉县格咱地区铜多金属矿整装勘查成果报告.230-460

3.1 分析测试方法

锆石U-Pb同位素定年在中国科学院地质与地球物理研究所完成。测试采用Finnigan Neptune型多接收电感耦合等离子体质谱仪(MC-ICP-MS)及Newwave UP 213 Nd:YAG激光剥蚀系统建立的锆石微区U-Pb定年,采样方式为单点剥蚀,数据采集采用所有信号同时静态方式接收。锆石年龄采用锆石91500或GJ1或TEM作为外标,标准和样品测试前先对空白进行测量。锆石Hf同位素测试所用仪器为Finnigan Neptune多接收等离子质谱和Newwave UP213紫外激光剥蚀系统(LA-MC-ICP-MS),实验过程中采用He作为剥蚀物质载气,根据锆石大小,剥蚀直径采用55μm或40μm,时使用国际上通用的锆石标样GJ1作为参考物质,分析过程中锆石标准GJ1的176Hf/177Hf测试加权平均值为0.282015±8(2σ,n=10)。相关仪器运行条件及详细分析流程见文献(吴福元等,2007)。

Mo、W成矿元素分析在国土资源部昆明矿产资源监督检测中心完成,分析测试仪器为Finnigan MAT制造的Element Ⅰ型Element Ⅱ高分辨率电感藕合等离子质谱仪(ICP-MS)。全岩的主和微量元素测试在国家地质实验测试中心完成,主量元素的测定采用等离子质谱仪(PE300D)测定,精度优于5%。微量元素和稀土元素的测定采用等离子光谱仪(PE8300)测定,检出限为0.001。

硫化物硫同位素组成分析在核工业北京地质研究院分析测试中心完成。首先对矿石、岩石中结晶态的硫化物进行分选 、挑纯处理,硫同位素测试以Cu2O做氧化剂制备SO2,并用冷冻法收集,然后用MAT253稳定同位素比质谱仪分析硫同位素组成,测量结果以V-CDT为标准,分析精度优于±2‰。Re-Os同位素分析与测定在国家地质实验测试中心完成,详细分析流程见文献(Du et al.,2004)。

3.2 铜厂沟斑岩型Mo-Cu矿床

铜厂沟斑岩型Mo-Cu矿床位于乡城-丽江结合带南缘,与扬子地台边缘坳陷带及甘孜-理塘结合带相接。区内构造运动复杂,岩浆作用强烈,除广泛发育华力西期及印支期中基性火山岩外,超基性、基性火山岩(时代?)也很发育。印支期和燕山期发育的浅成-超浅成中酸性侵入岩与该区Cu-Mo矿床的形成密切相关。矿区出露地层主要为上二叠统黑泥哨组(P2h)玄武岩,中三叠统北衙组(T2b)的碳酸盐建造及少量第四系更新统(Qp)、全新统(Qh)。铜厂沟Mo-Cu矿床主要产出于NE向的背斜构造之上,由于受后期构造的破坏,导致背斜出露残缺不全(图 2)。矿区主干断裂属铜厂断裂(F1)的北段,铜厂断裂沿背斜轴部通过,两侧有多条次级断裂,主要有近EW向、NE向、NW向三组。矿区岩浆岩发育,火山岩主要以玄武岩为主,分布于黑泥哨组(P2h)内。侵入岩有燕山期中酸性浅成岩体沿铜厂断裂附近分布,侵入岩体的出露面积均很小(0.1km2)。

图 2 铜厂沟斑岩型Mo-Cu矿床地质简图(实测) Fig. 2 Geological sketch map of the Tongchanggou porphyry Mo-Cu deposit

铜厂沟Mo-Cu矿主要工业矿体有3个。其中,KT1矿体赋存于黑泥哨组(P2h)玄武岩与北衙组一段(T2b1)灰岩接触带的矽卡岩中矿体呈似层状分布,总体向北西陡倾,局部倒转而倾向南东,倾角50°~84°,走向35°。矿体表现为Mo-Cu共生,元素垂直分带性明显,矿体浅部以Cu为主,中部为Mo-Cu共生矿,深部则以Mo为主(图 3)。 矿体厚0.91~33.83m,

图 3 铜厂沟Mo-Cu矿床0号勘探线剖面图(据云南省地质调查院,2013) Fig. 3 Profile map of No.0 exploration line of Tongchanggou Mo-Cu deposit

Cu品位0.11%~7.32%,Mo品位0.03%~1.29%。矿体呈似层状,总体向NW陡倾,倾角45°~79°,走向35°。本文研究工作对铜厂沟辉钼矿化花岗闪长斑岩进行矿物鉴定和Mo、W成矿元素含量分析表明,花岗闪长斑岩中具有较好的白钨矿化,W品位为0.01%~0.15%(表 1)。矿区共探获Mo金属量25.9万吨,矿床规模达大型。

表 1 铜厂沟斑岩型Mo-Cu矿床Mo、W成矿元素分析(wt%) Table 1 Elements analysis of Mo,W of Tongchnaggou porphyry Mo-Cu deposit(wt%)

铜厂沟Mo-Cu矿床是乡城-丽江结合带内近年来新的找矿重大成果,是燕山期斑岩型Mo多金属成矿作用的典型代表,目前矿床规模已达大型,找矿潜力巨大。李文昌等(2012)测定了铜厂沟Mo-Cu矿床辉钼矿Re-Os年龄为85±2Ma,余海军等(2015)测定铜厂沟花岗闪长斑岩LA-ICP-MS锆石U-Pb年龄为84.57±0.29Ma(MSWD=0.73)。本文研究获得黄铜矿辉钼矿化花岗闪长斑岩的锆石U-Pb年龄为84±1.7Ma(测试数据见电子版附表),辉钼矿化花岗闪长斑岩的锆石U-Pb年龄为83.6±1.1Ma(数据见电子版附表、图 4),表明铜厂沟斑岩型Mo-Cu矿床成岩成矿的时限基本一致,都属于燕山晚期构造-岩浆活动的产物。

图 4 铜厂沟花岗闪长斑岩锆石206Pb/238U-207Pb/235U谐和图 Fig. 4 Concordia diagram of zircon 206Pb/238U-207Pb/235U of Tongchanggou granodiorite porphyry
3.3 休瓦促W-Mo矿床

休瓦促W-Mo矿区位于格咱-乡城断裂以东,休瓦促-热林Cu-W-Mo-PbZn矿带北段。矿区出露地层为三叠系上统喇嘛亚组一段(T3lm1)岩性为灰色砂岩、粉砂岩、粉砂质绢云板岩(图 5)。矿区内发育一条NW-SE向断裂,地貌上表现呈串珠状高原湖泊及分水岭、深沟等特征,平硐中见断层泥化带,断层产状55°∠65°。矿区东南角有大面积出露的二长花岗岩体,岩性绝大多数为黑云母二长花岗岩,少数为花岗闪长岩,岩石具硅化及绢云母化。

图 5 休瓦促W-Mo矿区地质图(据云南省地质调查院,2013修编) Fig. 5 Geological map of Xiuwacu W-Mo mining area

辉钼矿主要赋存于二长花岗岩体内及构造裂隙或破碎带石英脉中,矿(化)体产出严格受断裂或裂隙控制。矿区分为东西两矿段,共圈出7个矿体(KT1~KT7)。Mo品位0.43%~6.37%,平均品位1.89%,矿床有用组分变化较均匀,矿床规模达中型。矿石矿物成分主要为辉钼矿、白钨矿,偶见少量黄铁矿、黄铜矿,未见氧化矿石。李建康等(2007)采用辉钼矿Re-Os测年技术获得休瓦促Mo矿床的形成时代为83±1Ma,本文LA-ICP-MS锆石U-Pb测年结果表明,辉钼矿化二长花岗岩的年龄为83.3±1.7Ma(MWSD=0.28)(数据见电子版附表、图 6),属于燕山晚期的构造-岩浆活动产物。值得一提的是,本文获得了休瓦促东矿段花岗闪长岩的锆石U-Pb年龄为202±3.5Ma(MWSD=0.64)(数据见电子版附表),这表明休瓦促W-Mo矿区存在印支晚期的岩浆活动及侵入事件。

图 6 休瓦促二长花岗岩、花岗闪长岩锆石206Pb/238U-207Pb/235U谐和图 Fig. 6 Concordia diagram of zircon 206Pb/238U-207Pb/235U of Xiuwacu monzonitic granite and granodiorite
3.4 热林Cu-Mo矿床

矿区位于乡城-丽江结合带北段,热林背斜轴部附近。出露地层主要为曲嘎寺组三段(T3q3),分布于热林背斜轴部,主要为薄-中厚层状灰岩夹绢云板岩;图姆沟组一段(T3t1),主要分布于热林背斜的两翼,岩性为板岩夹灰岩,底部为砾岩;图姆沟组二段(T3t2),分布于矿区两侧,热林背斜两翼,为绢云板岩和中性及中酸性火山岩类,火山岩组合为安山岩-流纹岩-英安岩-火山角砾岩,总体为一套碎屑岩-火山岩建造(图 7)。矿区断裂及褶皱发育,构造活动强烈,断裂及褶皱控制了地层的分布及花岗岩体的侵位。侵入岩以浅成-超浅成中酸性侵入岩为主,形成时代主要为印支晚期-燕山晚期,岩石由中酸性二长花岗岩、二长花岗斑岩和石英闪长玢岩组成。

图 7 中甸地区热林-亚杂矿区地质矿产略图(据尹光侯等,2009修改) 1-第四系,2-图姆沟组,3-曲嘎寺组,4-晚白垩世二长花岗岩,5-晚白垩世二长花岗斑岩,6-三叠纪石英闪长玢岩,7-三叠纪英安斑岩,8-花岗脉岩类,9-逆断层,10-中性/基性火山岩,11-Ⅰ-细粒二长花岗岩,12-Ⅱ-中粒二长花岗岩,13-Ⅲ-二长花岗斑岩,14-地质界线、岩性岩相界线 1-Quaternary System; 2-Tumugou Foramtion; 3-Qugasi Foramtion; 4-monzonitic granite of Late Cretaceous; 5-monzogranite porphyry of Late Cretaceous; 6-quartz diorite porphyrite of Triassic; 7-dacite porphyry of Triassic; 8-granite vein rocks; 9-reverse fault; 10-middle/basic volcanic rock; 11-Ⅰ-fine grained monzonitic granite; 12-Ⅱ-medium grained monzonitic granite; 13-Ⅲ-monzogranite porphyry; 14-geological boundary,lithological facies boundaries Fig. 7 The geological mineral sketch of Relin-Yaza mining area in Zhongdian(after Yin et al.,2009)

地表工程控制多个矿体,矿体均分布于二长花岗岩内部,呈脉状分布,矿体厚度较小,品位较低。其中,KT1分布于热林丫口附近,呈脉状分布,矿体走向345°,倾向NE。矿体厚度0.8~10m,Mo品位0.019%~0.086%,平均0.049%。钻孔深部揭露到二长花岗岩,说明深部存在含矿斑岩体(图 8)。从空间上看,矿区矿化具有明显分带性,即岩体的内外接触带附近具Mo-W矿化,离接触带500m左右地段具黄铜矿、弱辉钼矿化,离接触带约1800~2000m地段具PbZn矿化。推测深部可能具有较好的Mo-W找矿潜力。矿石以原生硫化矿为主,矿石矿物包括黄铁矿、方铅矿、闪锌矿、黄铜矿、锡石、白钨矿、磁铁矿、辉钼矿等;脉石矿物为石英、长石、绿泥石、黑云母等。

图 8 热林Cu-Mo多金属矿区0号勘探线剖面图(据云南省地质调查院,2013修改) Fig. 8 Profile map of No.0 exploration line of copper molybdenum polymetallic mining area in Relin

李建康等(2007)获得热林矿床辉钼矿的Re-Os等时线年龄为81.2±2.3Ma;尹光侯等(2009)采用黑云母Ar-Ar定年技术测得热林黑云母二长花岗斑岩的黑云母Ar-Ar坪年龄为82.01±0.86Ma,等时线年龄为81.7±1.1Ma。本文新获取2件二长花岗岩样品的锆石U-Pb年龄为79.4±2.0Ma和78±1.6Ma(数据见电子版附表、图 9),表明热林矿区成岩成矿事件同属燕山晚期的构造-岩浆活动产物。

图 9 热林花岗闪长斑岩锆石206Pb/238U-207Pb/235U谐和图 Fig. 9 Concordia diagram of zircon 206Pb/238U-207Pb/235U of Relin granodiorite porphyry
3.5 红山Cu-Mo矿床

红山Cu-Mo矿床位于香格里拉县格咱乡红山一带,与雪鸡坪铜矿毗邻,处于普朗-红山Cu多金属矿亚带西侧。矿区内次级断裂发育,多数与区域构造线一致,呈NW-SE向(图 10)。特别是成矿前的层间裂隙、节理发育,是产生矽卡岩化及矿化最有利的构造条件。矿区内火山岩很不发育,侵入岩均为浅成闪长玢岩或二长斑岩,数量较少。玢岩或斑岩均呈NW向延伸,侵入于曲嘎寺组砂泥质板岩中,岩石的蚀变较强。经工程揭露隐伏的花岗闪长斑岩具较好的Mo矿化(图 11),外接触带矽卡岩型层状Cu矿化发育。角岩中可见沿裂隙发育的细脉、网脉状Mo、Cu矿化。

图 10 红山Cu-Mo矿矿区地质图(据云南省地质调查院,2013修编) Fig. 10 Geological map of Hongshan Cu-Mo mining area)
图 11 红山Cu-Mo多金属矿区9号勘探线剖面图(据云南省地质调查院,2013修改) Fig. 11 Profile map of No.9 exploration line of Cu-Mo polymetallic mining area in Hongshan

红山Cu-Mo矿床的成矿时代问题一直存在较大的争议,石英二长斑岩隐伏岩体的全岩Rb-Sr等时年龄为216Ma,徐兴旺等(2006)获得辉钼矿的Re-Os年龄为77±2Ma;李建康等(2007)获得夕卡岩型矿石中的黄铁矿Re-Os年龄为77±18Ma;李文昌等(2011)获得红山深部隐伏岩体中辉钼矿矿Re-Os年龄为80.2Ma;黄肖潇等(2012)获得了红山成矿斑岩印支期、燕山期两期研究获得的锆石U-Pb年龄为216±3.3Ma和75.8±1.3Ma,孟健寅等(2013)获得红山夕卡岩型矿石中辉钼矿Re-Os等时线年龄年龄为79.32±0.87Ma;彭惠娟等(2014)认为是红山矿床是形成于晚燕山期的夕卡岩型铜矿床。本文获得的红山深部隐伏的花岗闪长斑岩的锆石U-Pb年龄为80±1.1Ma(数据见电子版附表、图 12),属于燕山晚期构造-岩浆活动产物。上述表明,红山Cu-Mo矿床存在两期成岩成矿作用,与石英二长斑岩、闪长玢岩有关的Cu金属成矿作用形成于印支期,而与深部隐伏花岗闪长斑岩有关的Cu-Mo成矿作用形成于燕山晚期,晚期对早期具有明显的叠加和改造作用。

4 燕山期斑岩Mo多金属矿床成矿系统 4.1 岩石系统

燕山晚期Mo多金属矿床含矿斑岩体岩性主要为花岗闪长斑岩和二长花岗岩,岩石常含斜长石、钾长石、石英、黑云母为主的斑晶。岩石地球化学分析结果表明(表 2),铜厂沟Mo-Cu矿床花岗岩闪长斑岩SiO2含量为63.02%~74.06%,全碱(Na2O+K2O)含量为6.97%~8.79%。岩石里特曼指数(σ)平均为2.30,属高钾钙碱性岩石系列(图 13)。岩石具较高的稀土元素总量∑REE=168.4×10-6~258.7×10-6。轻稀土较为富集(LREE/HREE为12.2~15.8),δEu为0.89~0.97,表现出无明显的δEu异常。休瓦促二长花岗岩SiO2含量为66.29%~77.36%,全碱(K2O+Na2O)含量为5.07%~9.24%,K2O含量普遍高于Na2O,属高钾钙碱性系列,少数为钾玄岩系列(图 13),显示准铝质-弱过铝质特征。稀土元素总量∑REE=155.0×10-6~356.2×10-6。δEu值为0.26~0.74,δEu表现具明显的负异常(图 14a)。负Eu异常的存在与Sr的强烈亏损一致,说明斜长石和角闪石等富Ca矿物分离结晶或者源区存留。热林二长花岗岩中的SiO2含量69.88%~74.58%,全碱(K2O+Na2O)含量为7.23%~8.62%,属高钾钙碱性系列,少数为钾玄岩系列(图 13),具准铝质-弱过铝质特征。岩石稀土元素总量∑REE=80.72×10-6~284.9×10-6。(La/Yb)N值为7.63~32.8,具较高的(La/Yb)N和LREE/HREE值,表明轻稀土元素分馏较强,且岩浆分异程度相对较高,δEu为0.46~0.83,具负Eu异常(图 14a)。负Eu异常可能与岩体中斜长石较多,剥蚀较深存在一定关系。红山隐伏花岗闪长斑岩中SiO2含量为67.59%~71.84%,全碱(K2O+Na2O)含量为7.6%~9.55%,属高钾钙碱性-钾玄岩系列(图 13),显示准铝质特征。稀土元素特征表明,花岗闪长斑岩具较高的稀土总量,∑REE平均为197.8×10-6。在球粒陨石标准化曲线中,呈现典型的右倾模式(图 14a),LREE/HREE值为7.61~24.7,为轻稀土元素(LREE)富集型。(La/Yb)N值为8.29~50.7,较高的(La/Yb)N和LREE/HREE值,表明轻稀土元素分馏较强,且岩浆分异程度相对较高。δEu值为0.58~0.82,δEu表现出弱的负异常。

表 2 乡城-丽江结合带燕山期花岗岩主量(wt%)和微量元素特征(×10-6) Table 2 The characteristics of major elements (wt%) and tracee lements (×10-6) of Yanshanian granite in xiangcheng-Lijiang suture zone
图 12 红山花岗闪长斑岩锆石206Pb/238U-207Pb/235U谐和图 Fig. 12 Concordia diagram of zircon 206Pb/238U-207Pb/235U of Hongshan granodiorite porphyry
图 13 花岗岩TAS分类图解及花岗岩SiO2-K2O图 Fig. 13 TAS classification diagram and SiO2 vs. K2O diagram of granite
图 14 稀球粒陨石标准化土元素配分图(a,标准化值据Sun and McDonough,1989)和原始地幔标准化微量元素蛛网图(b,标准化值据Taylor and McClennan,1985) Fig. 14 Chlondrite-normalized REE patterns(a,normalizing values after Sun and McDonough,1989)and primitive mantle-normalized trace element spider diagrams(b,normalizing values after Taylor and McClennan,1985)of Pulang composite pluton

综上所述,燕山晚期含矿斑岩具有较高的SiO2含量,岩石具高硅、富碱、高钾的特征,属于高钾钙碱性-钾玄岩的岩石系列,表明具有陆-陆碰撞带含矿斑岩的地球化学特征(侯增谦等,2012)。在微量元素和稀土元素的地球化学组成上,本区碰撞造山带花岗岩具明显的轻稀土(LREE)元素富集的特点,具有一定的负Eu异常,其中休瓦促含矿斑岩表现出较为明显的负Eu异常,热林、红山、铜厂沟含矿斑岩的δEu值表现出相对弱的负异常。本区燕山晚期含矿斑岩中微量元素的分布表现出明显的多峰谷模式,富集大离子亲石元素(LILE)Rb、U、Ba、K,亏损高场强元素(HFSE)Nb、Ta、P、Ti(图 14b),Th具有明显的正异常,表明含Mo斑岩体主要系大陆地壳深熔作用产物。此外,本区燕山期的含矿斑岩的Sr值变化于15.6×10-6~907×10-6,平均563.3×10-6;Y值变化于8.58×10-6~28.50×10-6,平均14.7×10-6,Sr/Y比值为1.2~64.3,具高Sr,低Y和高Sr/Y比值的特征,显示出典型的埃达克岩的地球化学属性。

4.2 流体系统

与长英质岩浆相关的成矿系统成矿流体是与成矿作用相关的特殊地地质流体,富含挥发性组分和卤素组分的岩浆热液,是研究矿床形成过程中物质来源、运移、汇聚成矿的关键(Hedenquist and Lowenstern,1994)。已有的研究工作表明,休瓦促W-Mo矿床的形成过程可分为硅酸盐-氧化物、石英-硫化物、碳酸盐岩3个阶段。其中,石英-硫化物阶段是主要的成矿阶段,以石英脉及金属硫化物的生成为代表。该阶段原生包裹体主要有气液两相水溶液包裹体(L+V)、富CO2包裹体(L+VCO2或L+LCO2+VCO2)以及纯液相包裹体(L)。成矿流体为含CO2的中高温(146~550℃),中低盐度(3.15%~12.51% Nacleqv)的H2O-NaCl热液体系,可能主要来自于岩浆期后热液(王新松等,2015)。红山Cu-Mo矿床流体包裹体的研究表明,热液成矿作用可划分为6个阶段,其中深部隐伏斑岩侵入引起的矿化主要为辉钼矿±黄铜矿+黄铁矿-石英脉。该类脉体中发育含NaCl子矿物三相(SL)、气液两相(VL)、富气相(LV)及含CO2三相(LC)4种类型的流体包裹体。包裹体均一温度变化范围为251.2~321.3℃,盐度为4.4%~7.4% Nacleqv(含子矿物包裹体盐度为31.96%~41.1% Nacleqv),成矿流体为中高温、低盐度NaCl-CO2-H2O体系热液(李文昌等,2013)。铜厂沟Mo-Cu矿床流体包裹体的研究表明,石英-硫化物阶段发育的辉钼矿±黄铜矿+黄铁矿-石英脉中流体包裹体的类型为气液两相包裹体(L+V)和气相包裹体(V),包裹体较为细小,以气液两相包裹体为主。包裹体的完全均一温度变化范围为265~365℃,盐度为4.5%~9.0% Nacleqv(项目组余海军提供),属于中高温,中低盐度的H2O-NaCl热液体系。

上述表明,燕山期形成的Mo多金属矿床主要成矿阶段流体包裹体多为气液两相水溶液包裹体(L+V),并发育富CO2包裹体(L+LCO2+VCO2),流体的均一温度变化于200~400℃,属中高温热液矿床,成矿流体形成于中低盐度的H2O-NaCl热液体系(图 15)。上述特征表明,燕山期Mo多金属矿床的形成具有相似的流体特征,也指示成矿流体的源区可能具有一致或同源性。

图 15 流体包裹体均一温度、均一盐度直方图 Fig. 15 Uniform temperature and salinity uniform histogram of fluid inclusions
4.3 物质系统 4.3.4 硫同位素特征

硫同位素示踪研究是探讨成矿物质来源的主要依据,本文对燕山晚期主要成矿斑岩的研究表明,铜厂沟Mo-Cu矿床中辉钼矿、黄铜矿、黄铁矿的金属硫化物的硫同位素组成δ34S的变化范围为-0.7‰~3.8‰,平均为0.79‰(表 3),δ34S值的分布呈现出明显的塔式分布特征(图 16),具有幔源硫的特征,说明矿床中辉钼矿、黄铜矿、黄铁矿的硫同位素组成与岩浆热液流体所形成的矿床硫同位素值相一致。休瓦促W-Mo矿中辉钼矿、黄铜矿、黄铁矿的硫同位素组成表明,δ34S值的变化范围为1.4‰~4.4‰,平均为2.65‰,变化范围较窄,具有高度均一的硫源,说明在硫化物的结晶沉淀过程中没有引起成矿流体中硫同位素的强烈分馏,呈现出岩浆硫的特征。红山Cu-Mo矿中黄铜矿、黄铁矿的硫同位素组成非常稳定,δ34S的变化范围为3.52‰~6.11‰,平均值为4.4‰,说明矿床的硫源主要为岩浆硫。俎波等(2013)的研究表明,红山Cu-Mo矿石硫化物环带的δ34S集中分布于3.8‰~5.23‰之间,极差仅1.42‰,均值为4.56‰; 集中分布且较低的δ34S值可能指示矿石硫统一来源于深源岩浆,这与本文的结果基本一致。本次研究所获得的燕山期代表性Mo多金属矿床金属硫化物δ34S值主要集中分布在0值附近,这与岩浆热液流体所形成的矿床硫同位素值分布相一致,表明在矿床形成过程中岩浆热液流体发挥了重要作用。

表 3 燕山期斑岩型Mo多金属矿床硫化物硫同位素组成 Table 3 Sulfur isotopic composition of sulfide of Tongchanggou porphyry molybdenum deposit
图 16 硫同位素分布直方图 资料来源:铜厂沟数据为本文实测值;休瓦促数据杨岳清等,2002;王新松等,2015;项目组;红山数据孟健寅等,2013;俎波等,2013;项目组 Fig. 16 Distribution histogram of sulfur isotopic
4.3.2 铅同位素特征

铅同位素的示踪研究表明,铜厂沟Mo-Cu矿床矿石矿物的铅同位素组成:206Pb/204Pb为18.3325~18.694,207Pb/204Pb为15.588~15.663,208Pb/204Pb为38.454~39.008,铅同位素组成较为稳定,显示正常铅的特征。铅同位素组成与特征参数(△β与△γ、V1与V2)之间具有明显的正相关性,依据铅构造模式判别和成因分类的综合分析,铜厂沟斑岩Mo-Cu矿床的铅主要来源于深部,并显示壳幔混合来源的特征(刘学龙等,2016)。红山Cu-Mo矿床的铅同位素分析表明,206Pb/204Pb为18.561~18.712,207Pb/204Pb为15.621~15.695,208Pb/204Pb为38.971~39.332,铅同位素组成变化稳定,主要为正常铅(刘学龙等,2013)。休瓦促Cu-Mo矿床金属硫化物铅同位素组成表明,206Pb/204Pb为18.646~18.936,207Pb/204Pb为15.650~15.710,208Pb/204Pb为38.857~39.404,各种金属硫化物也具有较为相似的Pb同位素组成,表明具有相同的来源;且与成矿相关的花岗岩的Pb同位素初始值具有相似的变化范围,说明具有相似或一致的物质来源(王新松等,2015)。本文综合区内3个Mo多金属矿矿床的铅同位素进行分析,样品主要落入上地壳和造山带范围,明显以壳源为主(图 17)。成矿物质来源与碰撞造山后的构造演化可能存在密切关系。说明铅同位素来源于造山带构造-岩浆事件中岩浆热液的带入。此外,铅的模式年龄说明下部存在较老的地质体,与区内岩浆弧是在古生代地壳基础上发育起来的推断相吻合(李文昌等,2012)。

图 17 铅同位素207Pb/204Pb-206Pb/204Pb构造环境判别图解(据Zartman and Doe,1981)与△β-△γ成因分类图解(据朱炳泉,1998) 数据来源:铜厂沟数据为本文实测值;休瓦促数据王新松等,2015;红山数据宋保昌等,2006;刘学龙等,2013 Fig. 17 Tectonic environment discrimination diagram of lead isotope 207Pb/204Pb vs. 206Pb/204Pb(after Zartman and Doe,1981)and genetic classification diagram of △β-△γ(after Zhu,1998)
4.3.3 铼-锇同位素特征

研究表明,Re-Os同位素体系不仅可以精确测定金属矿床中硫化物的形成时间,而且可以示踪成矿物质来源以及指示成矿过程中不同来源物质的混入程度(Foster et al.,1996)。一般金属硫化物的来源从幔源→壳幔混合源→壳源,辉钼矿中Re的含量会逐渐降低。李文昌等(2012)对前人的研究成果总结表明,成矿物质来源于地幔或以地幔物质为主的矿床,其辉钼矿的Re含量变化于10×10-6~1000×10-6;具有壳幔混合源的矿床,每克辉钼矿中的Re含量多在(n×10)×10-6;成矿物质完全来自于地壳(上地壳),其辉钼矿Re含量明显较低,为(1~n)×10-6。本文分别统计了燕山晚期铜厂沟、休瓦促、热林、红山Mo多金属矿床中金属矿物中Re含量变化,见表 4所示。

表 4 燕山晚期Mo多金属矿床Re-Os测年Re含量 Table 4 Re content of Re-Os dating of porphyry copper molybdenum deposit in Late Yanshanian

铜厂沟Mo-Cu矿床中辉钼矿中Re含量为16.41×10-6~43.53×10-6,表明其成矿物质源区具有壳幔混合来源的特征,且更靠近壳源说明成矿物质来源以壳源为主有幔源组分的加入。休瓦促W-Mo矿中辉钼矿中Re含量为4.3×10-6~23.72×10-6,与壳幔混源的物质源区的辉钼矿中Re的含量相一致,表明成矿物质的来源具有壳幔混合来源的特征。热林铜钼Cu-Mo矿中辉钼矿中Re含量为7.2×10-6~23.7×10-6,源区特征仍表现出壳幔混源的特点。红山Cu-Mo矿床中辉钼矿的Re含量为2.5×10-6~94×10-6,与壳源和壳幔混源岩浆热液矿床中Re的含量相近,说明红山Cu-Mo矿床成矿物质来源以壳源为主,并混入少量幔源物质,具有壳幔混源的特征。因此,乡城-丽江结合带燕山晚期的斑岩型Mo多金属矿床中Re含量变化范围较为一致,指示成矿物质来源具有相同或相似的物质源区,即以壳源物质源区为主且有少量幔源物质加入的特征。综上所述,燕山晚期主要Mo多金属矿床S、Pb同位素示踪及Re-Os测年Re含量的分析结果表明,燕山晚期的含矿斑岩成矿物质以壳源为主,但具有明显的壳幔混源的特征。

5 成矿系统分析 5.1 矿床分布特征

本区燕山期成矿作用与中酸性岩浆活动密切相关。从含矿岩石的类型来看,以花岗闪长斑岩、二长花岗岩为主,Cu-Mo-W多金属矿化主要发育在斑岩体内及接触带附近。区内代表性矿床,在成矿元素的组合上表现出较好的相似性。休瓦促W-Mo矿床在斑岩体及接触带内矿化产出以Mo-W-Cu为主,外带产出有Pb-Zn-Ag及Sn矿化。热林Cu-Mo矿床成矿元素的组合类型以Cu-Mo为主,但近期的研究工作也发现了二长花岗斑岩中具有较好的W矿化,因此其成矿元素的组合类型仍以Cu-Mo-W为主。红山Cu-Mo矿以往的研究工作主要以Cu为主,侵入岩体在地表的出露面积较小,仅见有闪长玢岩的岩枝出露,地表分布有大量的角岩及矽卡岩,矿化元素以Cu为主,外带分布有Pb-Zn-Ag多金属化。但深部的花岗闪长斑岩的揭露也显示具有Mo的矿化元素组合,但矿区深部花岗斑岩中尚未见有W的矿化显示。铜厂沟Mo-Cu矿床作为本区燕山晚期成矿的典型代表,成矿元素的组合类型具较好的分带性,主要表现为矿体浅部以Cu为主,中部表现为Cu-Mo共生,内接触附件发育有W矿化,深部矿化以Mo为主的特征。

从矿化类型与斑岩体的空间分布情况来看,燕山晚期Mo多金属矿床矿化类型自斑岩体向外依次为斑岩型Mo-Cu-W矿化→矽卡岩型Cu-Mo矿化→热液型Cu-Pb-Zn-Ag矿化。矿化体主要受斑岩体中构造裂隙系统的控制,构造裂隙主要发育于岩体的顶部或边部,而岩体中心没有明显的发育。含矿裂隙多呈脉状、网脉状分布,裂隙越是发育矿化也越富集,如铜厂沟含矿斑岩的边部发育的构造裂隙中分布有富集的辉钼矿化,受构造裂隙控制的矿体多呈透镜状或似层状产出。休瓦促W-Mo矿中,辉钼矿的发育多与石英脉的产出密切相关,在较宽的石英脉或细脉状石英脉体的边缘均有辉钼矿、黄铜矿、白钨矿的产出。同时,Cu-Mo-W矿化也表现出斑岩体本身就是矿体的特征。

另外,燕山晚期的含矿斑岩(休瓦促、热林、铜厂沟)中都发现了白钨矿,W矿化的分布同样受构造裂隙的控制,与石英脉的发育密切关联,多呈星点状、细脉状分布,在局部矿化富集地段可见有稠密浸染状的分布。由于热液矿床中成矿元素的沉淀随温度的变化呈现出不同的空间分带性。成矿元素的组合上由岩体中心向外表现出Mo→Mo-Cu(W)→Cu(Mo)→Pb-Zn-Ag由高温→中高温→中低温的变化规律,显示出岩体→内接触带→外接触带→围岩的成矿效应。

5.2 成矿机制

已有的研究成果表明,格咱岛弧地区Cu-Mo多金属矿床主要形成于印支期、燕山期中酸性岩浆侵入作用,区内成矿作用显示了叠加改造的普遍性、物质来源的多样性和成矿元素的复杂性(李文昌等,2015)。印支期,自晚三叠世甘孜-理塘洋壳向西发生俯冲以来,源自俯冲板片的脱水流体随地幔的交代诱发源区岩石熔融,导致区内钙碱性岩浆发生浅成-超浅成侵位和斑岩型Cu矿(普朗、雪鸡坪、春都)及矽卡岩型Cu多金属矿床(红山、浪都)的形成。进入燕山期以来,受区域造山后伸展作用的影响和制约,该区发育有陆壳重融主导型花岗岩浆的侵入作用(李文昌,2007),同时伴随有Mo、W、Cu多金属矿床的成矿作用的发生,形成了与二长花岗岩有关的Mo、W斑岩型矿床及花岗闪长斑岩侵入有关的Mo多金属成矿作用,其空间上形成了本区另一条重要的构造-岩浆成矿带。燕山晚期,陆壳重融型岩浆的侵位揭示了本区区域应力场发生了重大的转变。造山作用的P-T-t轨迹表明,在挤压应力场向伸展转变的过程会引起物质的熔融和成矿流体的形成,从而导致强烈的岩浆作用和流体活动,是成岩成矿作用发育形成的有利阶段(陈衍景,2006)。由于地壳的增厚,山根温度相对低于软流圈温度,这种热-物质结构产生潜在的重力失衡,导致去根作用或下地壳的拆沉作用,其结果是下部热的软流圈大幅上涌并导致地壳部分熔融,造成大规模的岩浆侵入事件,即S型花岗岩的大规模发育。

区内燕山期成矿作用与造山期后花岗岩的岩浆作用关系密切,由于岩浆作用形成的热液可以继承岩浆成分,因此岩浆中成矿元素富集特征可以反映到岩浆热液成矿作用中,燕山期岩浆富集Mo、W、Cu、Ag、Au等金属成矿元素,因此,为W、Mo多金属成矿作用提供了必要的物质基础,并形成了与二长花岗岩有关的石英脉型Mo(W)矿(休瓦促、沙都格勒)及沿构造破碎带发育的含Au-Ag-W-Sn的石英脉型多金属矿床(热林)。中部的红山深部发现隐伏花岗闪长斑岩体,岩体内发育有较好的Mo(Cu)矿化。南部铜厂沟经工程揭露也发现了深部的隐伏岩体,岩性主要为花岗闪长斑岩,岩体内Mo矿化发育较好,地质找矿工作已获得重大的突破。综上所述,乡城-丽江结合带燕山期侵入岩带呈南北向叠加于印支期岛弧岩浆岩带上,矿床多产于复式岩体内,以同源早期岩浆为围岩,成矿作用以接触交代和热液充填为主,形成了高温热液交代型,高温成矿元素组合为主的Mo-Cu-W为主的多金属成矿系统。

自晚三叠世以来,格咱地区大规模的构造-岩浆作用使得岩浆与地层中的成矿元素相互叠加富集,导致大量Cu、Mo、Au等成矿元素预富集形成初时矿源层之一,为之后的斑岩Mo多金属矿床的形成奠定了物质基础。至燕山晚期,区域性的碰撞后至伸展作用的构造转换阶段伴随有大规模的构造-岩浆-成矿事件。在造山环境,陆壳相对较薄,在早期挤压构造的基础上形成了大型的张性断裂系统,这为后期岩浆-流体-成矿作用的形成提供了良好的空间。燕山早期受印支期洋壳俯冲造山作用及压性构造的影响,造山带内发育的钾玄岩岩石系列的侵入岩可视为印支期洋壳俯冲消减作用的结束;燕山期与碰撞造山作用有关火山-侵入岩主要为高钾钙碱性的壳幔型花岗岩,燕山晚期受碰撞造山向造山后伸展构造转换的影响,区内形成了大量高钾钙碱性-钾玄岩系列的花岗岩,是碰撞造山作用过程"山根"拆沉作用的有效记录,这种构造环境的转换是区内燕山晚期斑岩Mo多金属矿床形成的主要机制。

6 成矿系统模式

乡城-丽江结合带燕山晚期斑岩型Cu-Mo-W成矿系统是燕山期碰撞造山作用过程的地质记录,成矿系统的分析与成矿模式的建立对于理解和深化对该区燕山期陆-陆碰撞造山过程的演化及相关的成岩成矿作用具有重要的地质意义。本区燕山晚期花岗岩浆的侵入及伴生的斑岩型矿床和大规模成矿流体的活动表明,滇西北地区晚白垩世存在区域性的构造-岩浆-流体-成矿作用,各斑岩型Cu-Mo矿床并不是独立或个别的成矿事件,而是属于区内与晚白垩世岩浆侵入作用相关的斑岩成矿系统的一部分。

区内Cu-Mo多金属成矿时代主要为印支期、燕山期及喜马拉雅期。通过各类型矿床查明资源储量占总资源储量的比例得出,Cu矿成矿时代以印支期为主,成矿系列为斑岩成矿系列;Mo矿床的形成时代主要为燕山晚期,也属于斑岩成矿系列。同时代、同类型的矿床,由于各成矿地质条件和控矿因素差异性,燕山晚期的斑岩成矿(系列)产出部位、类型各异。本文在综合区域地质背景,区域构造演化、典型矿床特征、成岩成矿物质系统的研究基础上,建立了乡城-丽江结合带燕山期斑岩型Mo多金属矿床成矿系统模式与预测模型,见图 18所示。

图 18 燕山晚期斑岩型Mo多金属矿床成矿系统模式 Fig. 18 The metallogenic system model of Mo-polymetallic deposit in Late Cretaceous
7 结论

(1) 乡城-丽江结合带燕山期含矿斑岩具高硅、富碱、高钾的特征,属高钾钙碱性-钾玄岩岩石系列,岩石富集轻稀土(LREE),具有一定的负δEu异常,微量元素分布表现出明显的多峰谷模式,明显富集大离子亲石元素(LILE)Rb、U、K,亏损高场强元素(HFSE)Nb、Ta、P、Ti,Th具有明显的正异常,表明含Mo斑岩体的形成具有大陆地壳的成因。燕山晚期含矿斑岩具有相似的稀土元素和微量元素地球化学组成特征,表明该区成矿斑岩是由同源岩浆分异演化而来。在成岩成矿时代上,区内4个代表性矿床含矿斑岩锆石U-Pb年龄值都非常接近,近乎同时期形成,属于燕山晚期构造-岩浆活动的产物。

(2) 燕山晚期的斑岩Mo多金属矿床的形成与区内壳源岩浆侵入密切相关,成矿物质以壳源为主,但具有壳幔混源的特征,含矿斑岩成岩成矿物质来源具有相似或一致的源区特征。在成矿元素的组合上,由斑岩体内向外带具有W、Mo→W、Mo、Cu→Pb、Zn、Ag的演化趋势,而红山矿床具有上部Cu-PbZn下部Cu-Mo的矿化元素分带和铜厂沟上部以Cu为主,中部Cu-Mo,深部Mo-W的垂直矿化元素分带模式。

(3) 本区燕山晚期花岗岩浆的侵入及大规模成矿流体的活动表明,滇西北地区晚白垩世存在区域性的构造-岩浆-流体-成矿作用,各斑岩型Mo多金属矿床并不是独立或个别的成矿事件,而是属于区内与晚白垩世岩浆侵入作用相关的斑岩成矿系统的一部分。

致谢 本文是综合多个项目的研究成果。野外地质调查工作得到云南省地质调查院矿产地质调查所、云南鼎立矿业公司的帮助和支持,在此表示感谢。感谢审稿专家对本文提出的宝贵意见与建议!
参考文献
[1] Cao DH, Wang AJ, Huang YF, Zhang W, Hou KJ, Li RP, Li YY. 2009. SHRIMP geochronology and Hf isotope composition of zircons from Xuejiping porphyry copper deposit, Yunnan Province. Acta Geologica Sinica , 83 (10) :1430–1435.
[2] Chen YJ. 2006. Orogenic-type deposits and their metallogenic model and exploration potential. Geology in China , 33 (6) :1181–1196.
[3] Deng J, Wang QF, Li GJ, Li CS, Wang CM. 2014. Tethys tectonicevolution and its bearing on the distribution of important mineraldeposits in the Sanjiang region, SW China. Gondwana Research , 26 (2) :419–437. DOI:10.1016/j.gr.2013.08.002
[4] Deng J, Wang QF, Li GJ, Hou ZQ, Jiang CZ, Danyushevsky L. 2015. Geology and genesis of the giant Beiya porphyry-skarn gold deposit, northwestern Yangtze Block, China. Ore Geology Reviews , 70 :457–485. DOI:10.1016/j.oregeorev.2015.02.015
[5] Deng J, Wang QF. 2016. Gold mineralization in China:Metallogenic provinces, deposit types and tectonic framework. Gondwana Research , 36 :219–274. DOI:10.1016/j.gr.2015.10.003
[6] Du AD, Wu SQ, Sun DZ, Wang SX, Qu WJ, Richard MHS, John M, Dmitry M. 2004. Preparation and certification of Re-Os dating reference materials:Molybdenite HLP and JDC. Geostandard and Geoanalytical Research , 28 (1) :41–52. DOI:10.1111/ggr.2004.28.issue-1
[7] Foster JG, Lambert DD, Frick LR, Maas R. 1996. Re-Os isotopic evidence for genesis of Archaean nickel ores from uncontaminated komatiites. Nature , 382 (6593) :703–706. DOI:10.1038/382703a0
[8] Hedenquist JW, Lowenstern JB. 1994. The role of magmas in the formation of hydrothermal ore deposits. Nature , 370 (6490) :519–527. DOI:10.1038/370519a0
[9] 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.
[10] Huang XX, Xu JF, Chen JL, Ren JB. 2012. Geochronology, geochemistry and petrogenesis of two periods of intermediate-acid intrusive rocks from Hongshan area in Zhongdian arc. Acta Petrologica Sinica , 28 (5) :1493–1506.
[11] Jin CH, Fan WY, Zhang H, Shen ZW, Zhang Y, Gao JH. 2013. Geochemical characteristics and genesis of the Zhuoma lead-zinc deposit in Yunnan Province. Geology in China , 40 (6) :1902–1911.
[12] Leng CB, Zhang XC, Chen YJ, Wang SX, Gou TZ, Chen W. 2007. Discussion on the relationship between Chinese porphyry copper deposits and adakitic rocks. Earth Science Frontiers , 14 (5) :199–210.
[13] Li JK, Li WC, Wang DH, Lu YX, Yin GH, Xue SR. 2007. Re-Os dating for ore-forming event in the late of Yanshan Epoch and research of ore-forming regularity in Zhongdian arc. Acta Petrologica Sinica , 23 (10) :2415–2422.
[14] Li WC.2007. The tectonic evolution of the Yidun island arc and the metallogenic model of the Pulang porphyry copper deposit, Yunnan, SW China. Ph. D. Dissertation. Beijing: China University of Geosciences .
[15] Li WC, Yin GH, Lu YX, Liu XL, Xu D, Zhang SQ, Zhang N. 2009. The evolution and 40Ar-39Ar isotopic evidence of the Pulang complex in Zhongdian. Acta Geologica Sinica , 83 (10) :1421–1429.
[16] Li WC, Yin GH, Liu XL, Lu YX, Xu D, Cao XM, Zhang SQ, Yang SR. 2010. 40Ar-39Ar ages and antimony mineralization of the Pananniu chang porphyry in the Zhongdian area of Yunnan Province. Geology and Exploration , 46 (2) :267–271.
[17] Li WC, Zeng PS, Hou ZQ, White NC. 2011. The Pulang porphyry copper deposit and associated felsic intrusions in Yunnan Province, Southwest China. Economic Geology , 106 (1) :79–92. DOI:10.2113/econgeo.106.1.79
[18] Li WC, Yin GH, Yu HJ, Lu YX, Liu XL. 2011. The porphyry metallogenesis of Geza volcanic magmatic arc in NW Yunnan. Acta Petrologica Sinica , 27 (9) :2541–2552.
[19] Li WC, Yu HJ, Yin GH, Cao XM, Huang DZ, Dong T. 2012. Re-Os dating of molybdenite from Tongchanggou Mo-polymetallic deposit in Northwest Yunnan and its metallogenic environment. Mineral Deposits , 31 (2) :282–292.
[20] Li WC, Yu HJ, Yin GH. 2013. Porphyry metallogenic system of Geza arc in the Sanjiang region, southwestern China. Acta Petrologica Sinica , 29 (4) :1129–1144.
[21] Li WC, Yin GH, Yu HJ, Liu XL. 2014. The Yanshanian granites and associated Mo-polymetallic mineralization in the Xiangcheng-Luoji area of the Sanjiang-Yangtze conjunction zone in Southwest China. Acta Geologica Sinica , 88 (6) :1742–1756. DOI:10.1111/1755-6724.12341
[22] Li WC, Liu XL. 2015. The metallogenic regularity related to the tectonic and petrographic features of Pulang porphyry copper orefield, Yunnan, and its ore-controlling characteristics. Earth Science Frontiers , 22 (4) :53–66.
[23] Liu XL, Li WC, Yin GH. 2012a. The metallogenic system of porphyry-skarn type Cu, Mo(Au) polymetallic deposits in Geza island arc of Yunnan province. Geology in China , 39 (4) :1007–1022.
[24] Liu XL, Li WC, Yin GH. 2012b. Lead isotope characteristics and tracing significance of ore metallogenic material in Geza arc metallogenic belt, Yunnan. Geoscience , 26 (3) :445–452.
[25] Liu XL, Li WC, Yin GH, Zhang N. 2013. The geochronology, mineralogy and geochemistry study of the Pulang porphyry copper deposits in Geza arc of Yunnan Province. Acta Petrologica Sinica , 29 (9) :3049–3064.
[26] Liu XL, Li WC, Zhang N, Yang FC, Kang J, Zhang B. 2016. Characteristics of sulfur and lead Isotopes and tracing of mineral sources in the Tongchanggou porphyry Mo(Cu) deposit at the southern edge of Geza arc belt, Yunnan. Geology in China , 43 (1) :209–220.
[27] Liu XL, Li WC, Yin GH, Zhang N. 2013. Distribution of magmatic rocks and geochemical characteristics in Geza Arc, western Yunnan. Acta Geologica Sinica , 87 (Suppl.1) :736–738.
[28] Meng JY, Yang LQ, Lü L, Gao X, Li JX, Luo YZ. 2013. Re-Os dating of molybdenite from the Hongshan Cu-Mo deposit in Northwest Yunnan and its implications for mineralization. Acta Petrologica Sinica , 29 (4) :1214–1222.
[29] Pang ZS, Du YS, Wang GW, Cao Y. 2009. Geological and geochemical feature and petrogenesis of Pulang complex, Yunan Province, China. Geological Bulletin of China , 28 (4) :531–537.
[30] Peng HJ, Li HY, Pei RF, Zhang CQ, Zhou YM, Tian G, Li JX, Long F. 2014. Mineralogical characteristics and metallogeny of the Hongniu-Hongshan copper deposit in Zhongdian area, Yunnan Province, China. Acta Petrologica Sinica , 30 (1) :237–256.
[31] Ren JB, Xu JF, Chen JL. 2011. Zircon geochronology and geological implications of ore-bearing porphyries from Zhongdian arc. Acta Petrologica Sinica , 27 (9) :2591–2599.
[32] Song BC, Cai XP, Xu XW, Wang J, Qin KZ, Zhang BL. 2006. Cenozoic hot-spring exhalative mineralization in Hongshan copper-polymetallic deposit in Zhongdian, NW Yunnan. Chinese Journal of Geology , 41 (4) :700–710.
[33] Sun SS, 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. DOI:10.1144/GSL.SP.1989.042.01.19
[34] Taylor SR, McClennan SM.1985. The Continental Crust:Its Composition and Evolution. Palo Alto, CA: Blackwell Scientific Publications .
[35] Wang SX, Zhang XC, Leng CB, Qin CJ, Ma DY, Wang WQ. 2008. Zircon SHRIMP U-Pb dating of the Pulang porphyry copper deposit, northwestern Yunnan, China:The ore-forming time limitation and geological significance. Acta Petrologica Sinica , 24 (10) :2313–2321.
[36] Wang XS, Bi XW, Hu RZ, Leng CB, Yu HJ, Yin GH. 2015. S-Pb isotopic geochemistry of Xiuwacu magmatic hydrothermal Mo-W deposit in Zhongdian area, NW Yunnan:Constrains on the sources of metal. Acta Petrologica Sinica , 31 (11) :3171–3188.
[37] Wu FY, Li XH, Zheng YF, Gao S. 2007. Lu-Hf isotopic systematics and their applications in petrology. Acta Petrologica Sinica , 23 (2) :185–220.
[38] Xu XW, Cai XP, Qu WJ, Song BC, Qin KZ, Zhang BL. 2006. Later Cretaceous granitic porphyritic Cu-Mo mineralization system in the Hongshan area, northwestern Yunnan and its significances for tectonics. Acta Geologica Sinica , 80 (9) :1422–1433.
[39] Yang LQ, Gao X, He WY. 2015. Late Cretaceous porphyry metallogenic system of the Yidun arc, SW China. Acta Petrologica Sinica , 31 (11) :3155–3170.
[40] Yang YQ, Hou ZQ, Huang DH, Qu XM. 2002. Collision orogenic process and magmatic metallogenic system in Zhongdian arc. Acta Geoscientia Sinica , 23 (1) :17–24.
[41] Yin GH, Li WC, Jiang CX, Xu D, Li JK, Yang SR. 2009. The evolution of Relin duplex rock masses in Yanshan phase and Ar-Ar dating age and copper-molybdenum mineralization characteristics of Zhongdian volcanic-magma arc. Geology and Exploration , 45 (4) :385–394.
[42] Yu HJ, Li WC, Yin GH, Wang JH, Jiang WT, Wu S, Tang Z. 2015. Geochronology, geochemistry and geological significance of the intrusion from the Tongchanggou Mo-Cu deposit, northwestern Yunnan. Acta Petrologica Sinica , 31 (11) :3217–3233.
[43] Zartman R E, Doe BR. 1981. Plumbotectonics-the model. Tectonophysics , 175 (1-2) :135–162.
[44] Zeng PS, Mo XX, Yu XH, Hou ZQ, Xu QD, Wang HP, Li H, Yang CZ. 2003. Porphyries and porphyry copper deposits in Zhongdian area, northwestern Yunnan. Mineral Deposits , 22 (4) :393–400.
[45] Zeng PS, Li WC, Wang HP, Li H. 2006. The Indosinian Pulang superlarge porphyry copper deposit in Yunnan, China:Petrology and chronology. Acta Petrologica Sinica , 22 (4) :989–1000.
[46] Zhu BQ.1998. The Theory and Application of Isotope System in the Earth Science:On the Crust Mantle Evolution in China. Beijing: Science Press .
[47] Zu B, Xue CJ, Ya XE, Wang QF, Liang HY, Zhao Y, Liu MT. 2013. Sulfide zonal texture and its geological significance of ores from the Hongshan copper deposit in Shangri-la, Yunnan Province, China. Acta Petrologica Sinica , 29 (4) :1203–1213.
[48] 曹殿华, 王安建, 黄玉凤, 张维, 侯可军, 李瑞萍, 李以科.2009. 中甸弧雪鸡坪斑岩铜矿含矿斑岩锆石SHRIMP U-Pb年代学及Hf同位素组成. 地质学报 , 83 (10) :1430–1435.
[49] 陈衍景.2006. 造山型矿床、成矿模式及找矿潜力. 中国地质 , 33 (6) :1181–1196.
[50] 侯增谦, 郑远川, 杨志明, 杨竹森.2012. 大陆碰撞成矿作用:Ⅰ. 冈底斯新生代斑岩成矿系统. 矿床地质 , 31 (4) :647–670.
[51] 黄肖潇, 许继峰, 陈建林, 任江波.2012. 中甸岛弧红山地区两期中酸性侵入岩的年代学、地球化学特征及其成因. 岩石学报 , 28 (5) :1493–1506.
[52] 金灿海, 范文玉, 张海, 沈战武, 张玙, 高建华.2013. 云南卓玛铅锌矿床地球化学特征及矿床成因. 中国地质 , 40 (6) :1902–1911.
[53] 冷成彪, 张兴春, 陈衍景, 王守旭, 苟体忠, 陈伟.2007. 中国斑岩铜矿与埃达克(质)岩关系探讨. 地学前缘 , 14 (5) :199–210.
[54] 李建康, 李文昌, 王登红, 卢映祥, 尹光侯, 薛顺荣.2007. 中甸弧燕山晚期成矿事件的Re-Os定年及成矿规律研究. 岩石学报 , 23 (10) :2415–2422.
[] 李文昌. 2007. 义敦岛弧构造演化与普朗超大型斑岩铜矿成矿模型. 博士学位论文. 北京:中国地质大学, 15-80 http://cdmd.cnki.com.cn/article/cdmd-82501-2007066872.htm
[56] 李文昌, 尹光候, 卢映祥, 刘学龙, 许东, 张世权, 张娜.2009. 中甸普朗复式斑岩体演化及40Ar-39Ar同位素依据. 地质学报 , 83 (10) :1421–1429.
[57] 李文昌, 尹光候, 刘学龙, 卢映祥, 许东, 曹晓民, 张世权, 杨舒然.2010. 中甸矿集区普朗-红山铜多金属成矿亚带北段帕纳牛场斑岩体40Ar-39Ar年龄及锑矿化. 地质与勘探 , 46 (2) :267–271.
[58] 李文昌, 尹光侯, 余海军, 卢映祥, 刘学龙.2011. 滇西北格咱火山-岩浆弧斑岩成矿作用. 岩石学报 , 27 (9) :2541–2552.
[59] 李文昌, 余海军, 尹光侯, 曹晓民, 黄定柱, 董涛.2012. 滇西北铜厂沟钼多金属矿床辉钼矿Re-Os同位素年龄及其成矿环境. 矿床地质 , 31 (2) :282–292.
[60] 李文昌, 余海军, 尹光候.2013. 西南"三江"格咱岛弧斑岩成矿系统. 岩石学报 , 29 (4) :1129–1144.
[61] 李文昌, 刘学龙.2015. 云南普朗斑岩型铜矿田构造岩相成矿规律与控矿特征. 地学前缘 , 22 (4) :53–66.
[62] 刘学龙, 李文昌, 尹光侯.2012a. 云南格咱岛弧斑岩-矽卡岩铜、钼(金)矿床成矿系统. 中国地质 , 39 (4) :1007–1022.
[63] 刘学龙, 李文昌, 尹光侯.2012b. 云南格咱岛弧岩浆成矿带铅同位素特征及成矿物质来源示踪. 现代地质 , 26 (3) :445–452.
[64] 刘学龙, 李文昌, 尹光侯, 张娜.2013. 云南格咱岛弧普朗斑岩型铜矿年代学、岩石矿物学及地球化学研究. 岩石学报 , 29 (9) :3049–3064.
[65] 刘学龙, 李文昌, 张娜, 杨富成, 康健, 张彪.2016. 云南格咱岛弧带南缘铜厂沟斑岩型铜钼矿床硫铅同位素特征与成矿物质来源示踪. 中国地质 , 43 (1) :209–220.
[66] 孟健寅, 杨立强, 吕亮, 高雪, 李建新, 罗跃中.2013. 滇西北红山铜钼矿床辉钼矿Re-Os同位素测年及其成矿意义. 岩石学报 , 29 (4) :1214–1222.
[67] 庞振山, 杜杨松, 王功文, 曹毅.2009. 云南普朗复式岩体地质地球化学特征及成因. 地质通报 , 28 (4) :531–537.
[68] 彭惠娟, 李洪英, 裴荣富, 张长青, 周云满, 田广, 李建新, 龙飞.2014. 云南中甸红牛-红山矽卡岩型铜矿床矿物学特征与成矿作用. 岩石学报 , 30 (1) :237–256.
[69] 任江波, 许继峰, 陈建林.2011. 中甸岛弧成矿斑岩的锆石年代学及其意义. 岩石学报 , 27 (9) :2591–2599.
[70] 宋保昌, 蔡新平, 徐兴旺, 王杰, 秦克章, 张宝林.2006. 云南中甸红山铜-多金属矿床新生代热泉喷流沉积型矿床. 地质科学 , 41 (4) :700–710.
[71] 王守旭, 张兴春, 冷成彪, 秦朝建, 马德云, 王外全.2008. 滇西北普朗斑岩铜矿锆石离子探针U-Pb年龄:成矿时限及地质意义. 岩石学报 , 24 (10) :2313–2321.
[72] 王新松, 毕献武, 胡瑞忠, 冷成彪, 余海军, 尹光侯.2015. 滇西北中甸地区休瓦促岩浆热液型Mo-W矿床S、Pb同位素对成矿物质来源的约束. 岩石学报 , 31 (11) :3171–3188.
[73] 吴福元, 李献华, 郑永飞, 高山.2007. Lu-Hf同位素体系及其岩石学应用. 岩石学报 , 23 (2) :185–220.
[74] 徐兴旺, 蔡新平, 屈文俊, 宋保昌, 秦克章, 张宝林.2006. 滇西北红山晚白垩世花岗斑岩型Cu-Mo成矿系统及其大地构造学意义. 地质学报 , 80 (9) :1422–1433.
[75] 杨立强, 高雪, 和文言.2015. 义敦岛弧晚白垩世斑岩成矿系统. 岩石学报 , 31 (11) :3155–3170.
[76] 杨岳清, 侯增谦, 黄典豪, 曲晓明.2002. 中甸弧碰撞造山作用和岩浆成矿系统. 地球学报 , 23 (1) :17–24.
[77] 尹光候, 李文昌, 蒋成兴, 许东, 李建康, 杨舒然.2009. 中甸火山-岩浆弧燕山期热林复式岩体演化与Ar-Ar定年及铜钼矿化. 地质与勘探 , 45 (4) :385–394.
[78] 余海军, 李文昌, 尹光候, 王建华, 姜文涛, 吴松, 唐忠.2015. 滇西北铜厂沟Mo-Cu矿床岩体年代学、地球化学及其地质意义. 岩石学报 , 31 (11) :3217–3233.
[79] 曾普胜, 莫宣学, 喻学惠, 侯增谦, 徐启东, 王海平, 李红, 杨朝志.2003. 滇西北中甸斑岩及斑岩铜矿. 矿床地质 , 22 (4) :393–400.
[80] 曾普胜, 李文昌, 王海平, 李红.2006. 云南普朗印支期超大型斑岩铜矿床:岩石学及年代学特征. 岩石学报 , 22 (4) :989–1000.
[81] 朱炳泉. 1998. 地球科学中同位素体系理论与应用:兼论中国大陆壳幔演化. 北京: 科学出版社 .
[82] 俎波, 薛春纪, 亚夏尔, 王庆飞, 梁华英, 赵毅, 刘铭涛.2013. 云南香格里拉红山铜矿石硫化物环带及地质意义. 岩石学报 , 29 (4) :1203–1213.