岩石学报  2016, Vol. 32 Issue (8): 2248-2264   PDF    
滇西芦子园Pb-Zn-Fe多金属矿床蔷薇辉石地球化学特征及其成因意义
邓明国1, 徐荣1, 王朋1, 孙柏东2, 曾磊3, 余海军2, 王涛2, 沙建泽2     
1. 昆明理工大学国土资源工程学院, 昆明 650093 ;
2. 云南省地质调查局, 昆明 650051 ;
3. 四川省煤田地质局137队, 达州 635000
摘要: 芦子园超大型Pb-Zn-Fe多金属矿床是保山地块上系列同类矿床的典型代表,矿体呈似层状顺层产出,含矿围岩以矽卡岩为主;地质、物探、化探、遥感等综合推断深部存在隐伏岩体,但已知矿床内未见到隐伏岩体且未获得直接的成矿年龄等,矿床成因一直存在较大分歧。矿床深部发现蔷薇辉石与Fe、Pb、Zn、Cu矿化具密切的共伴生关系以及往深部蔷薇辉石含量增加、颗粒增大等特征,具有重要的矿床成因指示意义。本文对蔷薇辉石的主量元素、微量元素、稀土元素以及H、O同位素进行了研究,结果表明:蔷薇辉石主要化学组成SiO2 47.01%~48.66%,MnO 34.42%~40.19%,CaO 7.52%~8.51%,FeO 2.85%~8.76%,MgO 0.22%~0.28%,TiO2和Al2O3含量均较低;微量元素Li、Be、Sc、Zr、Rb、Sr等不相容元素不亏损或弱亏损,相容元素Cr、Co、Ni明显亏损;稀土元素总量较低;包裹体δDV-SMOW为-81.7‰~-73.3‰,δ18OH2O为5.4‰~6.1‰,在δD-δ18OH2O图解中投点均落入岩浆水区域;矿石矿物Y/Ho值与周边燕山期志本山、漕涧花岗岩较为一致。上述均表明芦子园矿床成矿物质主要来源于深部岩浆岩,综合分析认为,芦子园Pb-Zn-Fe多金属矿床为与燕山期隐伏中酸性岩浆岩有关的远程矽卡岩型矿床。
关键词: 蔷薇辉石     主量元素     微量稀土元素     H、O同位素     Pb-Zn-Fe多金属矿床     滇西芦子园    
Geochemistry of the rhodonite in the Luziyuan Pb-Zn-Fe polymetallic deposit in West Yunnan and their genesis significance
DENG MingGuo1, XU Rong1, WANG Peng1, SUN BaiDong2, ZENG Lei3, YU HaiJun2, WANG Tao2, SHA JianZe2     
1. Faculty of Land Resource Engineering Kunming University of Science and Technology, Kunming 650093, China ;
2. Yunnan Bureau of Geological Survey, Kunming 650051, China ;
3. Geological Team 137 of Sichuan Coalfield Geology Bureau, Dazhou 635000, China
Abstract: Luziyuan super-large Pb-Zn-Fe polymetallic deposit is a typical representative of the similar deposits located in Baoshan block. The ore bodies occur as statiform-like along the strata, and the ore-bearing wall-rocks are mainly skarns. The geological, geophysical, geochemical prospecting and remote sensing interpretation evidences all showed that the concealed rock body existed in the deep. While the inferred hidden body and the accurate mineralization age hasn't been found and obtained in the known deposit, thus the ore genesis still remains controversial. This research discovers that there is an obvious Symbiotic relationship between rhodonite and iron, lead-zinc, and copper ores. In addition, the amount and grain size of rhodonite gradually increase down the deep. All these phenomena observed in the deposit are of great significance for us to understand the origin of it. Major, trace, rare earth elements and H-O isotope compositions of the rhodonite are studied in this paper, with the results that the rhodonite mainly comprises SiO2 of 47.01%~48.66%, MnO of 34.42%~40.19%, CaO of 7.52%~8.51%, FeO of 2.85%~8.76%, MgO of 0.22%~0.28%, and lower TiO2 and Al2O3. The incompatible elements (such as Li, Be, Sc, Zr, Rb, Sr) are not depleted, or weakly depleted in the rhodonite; while some compatible elements including Cr, Co, and Ni are obviously depleted; the REEs' total contents of both rhodonite and sphalerite are lower; δDV-SMOW and δ18OH2O values of the rhodonite range from -81.7‰ to -73.3‰, 5.4‰ to 6.1‰, respectively, which are all plotted within magmatic waters' area on the δDV-SMOW vs. δ18OH2O diagram. Y/Ho value of ore minerals is relatively identical to the surrounding Yanshanian granites such as Zhibenshan and Caojian. All of the above-mentioned features demonstrate that the metallogenic fluids of Luziyuan deposit are characterized by the deep magmatic hydrothermal solutions. Therefore, it is suggested that the Luziyuan Pb-Zn-Fe polymetallic deposit is a distal-skarn type Pb-Zn deposit closely related to the concealed intermediate and acidic igneous rocks of Yanshan Period.
Key words: Rhodonite     Major elements     Trace elements     H,O isotope     Pb-Zn deposit     Luziyuan in West Yunnan    

滇西保山地块是我国西南三江成矿带中重要的Fe-Cu-Pb-Zn-Ag-W-Sn及稀有金属富集区,区内已发现保山核桃坪大型Pb-Zn-Cu多金属矿床、西邑大型Pb-Zn矿床、施甸东山、龙陵勐兴中大型Pb-Zn矿床和镇康芦子园超大型Pb-Zn-Fe(Cu)多金属矿床以及数十处矿(化)点(图 1b)。其中,镇康芦子园Pb-Zn-Fe(Cu)多金属矿床因其矿产资源量大、围岩蚀变强烈、矿物组合多样和矿化元素及伴生元素多等特点而备受矿床学界关注,已成为研究保山地块系列同类热液矿床的典型代表。前人从矿床地质特征、矿床地球化学、成岩成矿年代学及同位素地球化学等方面做了诸多研究(邓必方,1995; 夏庆霖等,2005; 董文伟和陈少玲,2007; 陶琰等,2010; 朱飞霖等,2011; 杨玉龙等,2012; 林冰霞,2013),但迄今对该矿床的成因认识仍众说纷纭,其中比较主流的成因观点为沉积热液改造矿床(邓必方,1995; 李文昌和莫宣学,2001; 陈永清等,2005; 董文伟和陈少玲,2007; 卢映祥,2011)以及与深部隐伏岩体有关的热液矿床(夏庆霖等,2005; 陶琰等,2010; 朱飞霖等,2011; 杨玉龙等,2012)。产生上述分歧的可能原因:一是矿体本身的顺层产出特点与层控矿床的地质特征相似;二是矿区周边及深部在时空上至今尚未发现与矿床有直接成因联系的岩体;三是赋矿地层的Pb、Zn等成矿元素背景值较高(邓必方,1995; 董文伟和陈少玲,2007),但未有充分的研究证据显示地层具"矿源层"特征。综上认为,对成矿物质来源的深入探讨是研究芦子园矿床成因的关键所在。

本次工作对芦子园矿床中发现的蔷薇辉石进行了系统调查,野外观察发现蔷薇辉石的产出与Pb、Zn、Fe矿化有密切的共伴生关系,在蔷薇辉石分布广的地段矿体规模变大、品位变富,且普遍存在铅锌矿直接交代(或充填)蔷薇辉石的现象,这些地质事实暗示蔷薇辉石在该矿床中具有重要的成因指示意义。本文在前人研究的基础上,开展了蔷薇辉石的主量元素、微量元素、稀土元素和H、O同位素研究,为探讨芦子园矿床的成矿物质来源、矿床形成环境及成因提供直接证据。

1 区域地质背景及矿床地质特征 1.1 区域地质背景

三江特提斯造山带是东特提斯构造带的重要组成部分,因金沙江、澜沧江、怒江纵贯全区而得名(图 1a; Deng et al.,2014ab; Wang et al.,2016; Li et al.,20152016)。保山地块位于三江造山带南段,东界澜沧江断裂,西界怒江断裂,北部由于澜沧江断裂和怒江断裂汇拢而消失,南部延伸至缅甸,构成滇缅泰马地体的北延部分(图 1a,b; 陶琰等,2010; Deng et al.,2014b)。早二叠世之前保山地块与腾冲、羌塘、拉萨地块附属于冈瓦纳大陆北东部,于早二叠世萨克马尔期末保山地块和腾冲地块开始从冈瓦纳大陆裂离向北漂移,导致中特提斯洋的开启;晚三叠世至早白垩世,拉萨地块向北漂移,中特提斯洋进入闭合阶段,导致腾冲地块和保山地块碰撞而形成高黎贡碰撞构造带(Metcalfe,2011; 莫宣学和潘桂堂,2006; Deng et al.,2014b)。保山和思茅地块之间的昌宁-孟连缝合带曾为古特提斯主洋,洋盆开启于中泥盆世,中三叠世闭合(Wang et al.,2000; 邓军等,2014)。保山地块基底为新元古代-寒武系公养河群,系一套类复理石砂岩、板岩及硅质岩,沉积盖层由古生代和中生代的台地浅海相碎屑岩、碳酸盐组成,仅上石炭统、上三叠统和中侏罗统地层中夹有中基性火山岩。

图 1 芦子园铅锌矿床地质简图(据云南省地质调查院,2012修改) (a)三江造山带大地构造位置图(Deng et al.,2014b);(b)保山地块区域简图(据Deng et al.,2014b修绘);(c)镇康地区矿产分布图(夏庆霖等,2005):①岩脚山铅锌矿;②枇杷水铅锌矿;③草坝寨铅锌矿;④天生桥铁矿;⑤放羊山铅锌矿;⑥小河边铁矿;⑦芦子园铅锌铁多金属矿;⑧水头山铅锌矿;⑨乌木兰锡矿;⑩小干沟金矿;(d)31号勘探线剖面示意图(据云南省地质调查局,2012修绘).1-寒武系保山组第一段粉砂质粘板岩夹泥质条纹结晶灰岩;2-寒武系沙河厂组第三段大理岩、大理岩化灰岩;3-寒武系沙河厂组第二段大理岩化灰岩;4-辉绿岩脉;5-铅锌矿体;6-铁矿体;7-铜矿体;8-实测断层及编号;9-推测断层及编号;10-矽卡岩;11-蔷薇辉石 Fig. 1 Geological sketch map of Luziyuan Pb-Zn deposit

① 云南省地质调查院.2012.云南省镇康县芦子园铅锌铁多金属矿勘查核实报告

保山地块内主要经历了早古生代和中-新生代花岗岩侵入活动,具幕式、阶段性演化特征。早古生代以地块南部的平河花岗岩基为代表,锆石U-Pb年龄为502~466Ma(董美玲等,2012; Wang et al.,2013; 禹丽等,2014)。古-中生代印支期花岗岩体主要为木厂A型花岗岩,是由钠闪英碱正长岩、钠闪花岗岩、霓石花岗岩组成的复式岩体,曲晶石、霓霞石U-Pb年龄为253~241Ma,该期岩体具Sn、Zn、∑Ce等综合矿化,具体成矿作用尚不清楚(顾影渠等,1988)。中-新生代岩浆活动形成的岩体主要有志本山、柯街及曹涧酸性花岗岩体。其中,燕山晚期岩浆活动主要形成志本山、柯街花岗岩,由黑云母花岗岩、二云母花岗岩及浅色花岗岩组成,锆石U-Pb年龄分别为126.7±1.6Ma和93±13Ma,相关的成矿作用为Pb、Zn、Fe、Cu等矿化,形成如保山核桃坪、镇康芦子园等Pb-Zn-Fe(Cu)多金属矿床(陶琰等,2010);喜山期岩浆活动形成的岩体主要为曹涧二云母花岗岩,锆石U-Pb年龄为73±0.19Ma,可能与云龙铁厂锡矿等成矿相关(陶琰等,2010; 禹丽等,2014)。除此,区内还发育大量辉绿岩脉及岩墙。

1.2 矿床地质特征

镇康芦子园矿区是保山地块上典型的Fe-Cu-Pb-Zn-Ag-W-Sn多金属矿化集中区,区内分布有乌木兰锡矿、小河边铁矿、天生桥铁矿、芦子园铅锌矿、放羊山铅锌矿、水头山铅锌矿等一系列高-中-低温多金属矿床组合(图 1c),其中芦子园矿床是本区热液矿床的典型代表。该矿床位于保山-镇康地块南部镇康复背斜之芦子园次级背斜之上,主要赋矿地层为上寒武统沙河厂组二段大理岩、大理岩化灰岩及板岩等。矿床具以下显著地质特征:

① 矿体沿NE向层间断裂破碎带呈层状、似层状整合产出于矽卡岩及大理岩中,受NE向区域构造及层间断裂构造控制明显(图 1);

② 矿床规模达超大型,Pb+Zn资源量达4.2Mt,Fe矿石资源量达300Mt;

③ 矿体走向及倾向延伸大(走向延长大于3.5km,倾向延伸大于1.5km),深部找矿潜力大;

④ 矿化元素(Fe、Cu、Zn、Pb、Ag)及伴生元素(W、Sn、Cd、Co、Mn)多,综合利用价值高;

⑤ 蚀变矿化分带明显:上部大理岩型Pb-Zn(Ag)→中部大理岩-矽卡岩型Pb-Zn-Fe-Cu→下部矽卡岩型Fe-Cu(图 1d);

⑥ 围岩蚀变强烈(矽卡岩化、大理岩化、硅化、碳酸盐化、黄铁矿化)、金属矿物组合多(闪锌矿、方铅矿、磁铁矿、黄铜矿、赤铁矿、黄铁矿、磁黄铁矿、斑铜矿、辉铜矿、毒砂等),钙矽卡岩矿物(钙铁榴石-钙铝榴石、透辉石、钙铁辉石、钠长石、斜长石、透闪石等)和锰矽卡岩矿物(蔷薇辉石、锰铝榴石、锰阳起石等)发育,以锰质矽卡岩为特征;

⑦ 由浅表至中深部,伴随蔷薇辉石的增多,铅锌矿体品位变富、规模增大,与铅锌矿化紧密伴生;

根据野外地质调查、手标本及显微镜观察等综合研究,可将芦子园矿床划分为矽卡岩期(主要为磁铁矿化)、热液成矿期(Pb、Zn、Cu等矿化)和表生期三个主要成矿期(表 1)。

表 1 芦子园铅锌矿床矿物生成顺序表 Table 1 The mineral arisen sequence of Luziyuan lead-zinc deposit
2 蔷薇辉石岩相学特征

芦子园矿床蔷薇辉石主要产于矿床中深部的断裂构造破碎带,在Ⅰ、Ⅱ号主矿体及附近广泛分布。新鲜的蔷薇辉石为玫瑰红至桃红色,常见其包裹残余角砾状大理岩(图 2a,b),主要呈团块状(图 2b,e,f)、角砾状(图 2b,c)和脉状(图 2d,e)产出,与铅锌矿共伴生关系密切。推测蔷薇辉石为深部热液沿构造滑脱部位交代大理岩而成,受后期持续构造应力作用,团块状蔷薇辉石被压碎为角砾状,以阳起石为主的矽卡岩矿物对其交代形成灰绿色放射状环带、锯齿状交代边,显示出明显的热液交带特征(图 2a,b)。

图 2 芦子园铅锌矿床蔷薇辉石野外及手标本照片 (a)1350中段蔷薇辉石(LZY11-4、LZY11-5、LZY13-1、LZY13-2、LZY26采样中段);(b)1220中段蔷薇辉石(LZY21、LZY24、LZY31、LZY40采样中段);(c)角砾状蔷薇辉石;(d)脉状蔷薇辉石;(e)阳起石蔷薇辉石矽卡岩,蔷薇辉石呈团块状和脉状产出;(f)团块状蔷薇辉石,见方解石脉穿插蔷薇辉石而未切穿绿帘石和阳起石 Fig. 2 Field photographs of rhodonite in Luziyuan Pb-Zn deposit

在手标本上,蔷薇辉石为粉红-桃红色,结晶颗粒较大(多为0.5~3mm),可见两组解理发育,常被铅锌矿及湿矽卡岩阶段矿物(阳起石、绿帘石等)充填交代成角砾状、碎粒状等(图 2b-d);据野外地质特征及手标本判断,蔷薇辉石的形成先于阳起石和绿帘石等含水硅酸盐矿物(图 2b,f)。在显微镜下,蔷薇辉石为淡红色,呈粒状、板柱状连晶,因有尘点状铁锰质分布,其表面多为粗糙不洁状,正高突起。正交下具一级橙干涉色,呈厚板柱状、粒状产出,简单双晶发育,可见条纹叶片状透辉石嵌于蔷薇辉石解理中呈连生现象,为含闪锌矿的蔷薇辉石矽卡岩(图 3c,d);金属矿物以闪锌矿为主,含少量黄铜矿、方铅矿等,沿蔷薇辉石的脉壁呈不规则细脉分布,部分充填于蔷薇辉石及围岩的裂隙中或矿物颗粒间形成细脉,表明金属矿物的形成晚于蔷薇辉石(图 3a,b)。

图 3 芦子园铅锌矿床蔷薇辉石镜下照片 (a)闪锌矿及黄铜矿呈细脉穿入蔷薇辉石(反射光);(b)闪锌矿分布于蔷薇辉石脉边缘(反射光);(c)蔷薇辉石呈板柱状平行连晶(单偏光);(d)蔷薇辉石呈平行连晶生长,与闪锌矿及透辉石等共生(正交偏光) Fig. 3 Photomicrographs of rhodonite in Luziyuan Pb-Zn deposit

矿区的蔷薇辉石主要有二种:一是团块状、角砾状蔷薇辉石,为热液交代构造破碎带中的团块状、角砾状及透镜状大理岩而成,主要形成于渗滤交代作用,常见蔷薇辉石中包裹未交代完全的大理岩残块;二是脉状蔷薇辉石,为热液沿构造裂隙充填交代大理岩的过程中形成。两种产状的蔷薇辉石颜色及共伴生矿物组合都较为一致,团块状、角砾状蔷薇辉石的结晶粒度略大于脉状蔷薇辉石。

3 样品采集及分析方法

全岩样品采自芦子园矿床的典型钻孔及坑道,主要针对大理岩、矽卡岩、构造岩及矿石分别进行采样。单矿物样品采自矿床深部蔷薇辉石分布广泛的1220、1350中段Ⅰ、Ⅱ号主矿体,采集的蔷薇辉石、阳起石及闪锌矿样品均具有明显的共伴生关系。样品分析方法如下所述:

(1) 主量元素分析:挑选新鲜未风化的岩(矿)石破碎后在陶瓷研钵中研磨至200目,岩石粉末的主量元素测试在西北有色金属地研测试中心采用等离子质谱(ICP-MS)定量测试完成,测试误差小于5%。

(2) 电子探针成分分析:蔷薇辉石电子探针片的磨制在昆明理工大学地球科学系完成,将磨制好的探针片先在显微镜下观察、拍照,并圈出典型的蔷薇辉石矿物。电子探针成分分析在武汉地调中心(宜昌地质矿产研究所)完成,仪器型号为JXA-8100型。

(3) 微量、稀土元素分析:蔷薇辉石、阳起石、闪锌矿单矿物微量、稀土元素测试在中国科学技术大学地球与空间科学学院固体同位素地球化学实验室完成。将含蔷薇辉石、阳起石和闪锌矿的样品磨至60~80目,在双目镜下挑选出纯度达99%的单矿物样品,用玛瑙研钵研磨至200目,采用等离子质谱(ICP-MS)进行微量、稀土元素定量测试,分析误差小于5%,具体分析方法和流程见Liang et al.(2000)

(4) H、O同位素分析:蔷薇辉石H、O同位素测试在北京核工业地质研究院完成。将挑选好的蔷薇辉石在装有玻璃碳的陶瓷管里与C发生充分还原反应,生成供测试用的CO及H2气体。CO及H2在高纯氦气流的带动下,经过色谱柱分离后分别进入质谱MAT253进行分析。测量结果以SMOW为标准,记为δ18OV-SMOW及δDV-SMOW,分析精度分别优于±0.2‰及±1‰,具体分析方法见刘汉彬等(2013)

4 分析结果 4.1 岩(矿)石主量元素组成

芦子园矿床中岩(矿)石主量元素测试结果见表 2。CaO从大理岩→矽卡岩→构造岩→矿石含量逐渐降低,与其它主量元素主要表现为负相关关系。这一特征表明,蚀变过程是一Si、Fe、Mn以及Al等元素带入、Ca强烈带出的过程,从而使CaO的含量逐渐降低,形成以阳起石、绿帘石、透辉石、绿泥石、蔷薇辉石及石榴子石为主的一套锰质矽卡岩建造。锰质矽卡岩仅产于矽卡岩型Pb-Zn-Ag矿床或Pb-Zn矿床中,为矽卡岩型Pb-Zn(Ag)矿床的重要找矿标志(赵一鸣和李大新,2003)。Mn作为蔷薇辉石的主要阳离子成分,在芦子园矿区具有重要的找矿意义,MnO从大理岩→矽卡岩→矿石→蔷薇辉石含量逐渐增加。大理岩中MnO的含量最低,变化范围为0.08%~0.29%,平均为0.22%;在矽卡岩中为0.46%~2.31%,平均1.42%,是大理岩含量的6.5倍;在矿石中明显增高,为0.17%~3.71%,平均2.25%,是大理岩含量的10倍(表 2);在蔷薇辉石中含量最高,为34.42%~40.19%(表 3)。

表 2 芦子园矿区岩(矿)石主量元素成分分析(wt%) Table 2 Major element compositions of the rock from Luziyuan Pb-Zn deposit(wt%)
表 3 蔷薇辉石及主要矽卡岩矿物电子探针分析结果(wt%) Table 3 Electron microprobe analyses of the rhodonite and altered rocks from Luziyuan Pb-Zn deposit(wt%)
4.2 蔷薇辉石电子探针分析结果

表 3可以看出芦子园矿床蔷薇辉石的主要化学成分SiO2为47.01%~48.66%,平均47.60%,其含量介于石榴子石和阳起石之间;MnO含量最高,为34.42%~40.19%,平均37.14%;CaO为7.52%~8.51%,平均8.13%;FeO变化范围较大,为2.85%~8.76%,平均6.03%;MgO为0.22%~0.28%,平均0.24%;TiO2和Al2O3含量均较低。与产于福建马坑铁矿床外接触带矽卡岩中的钙蔷薇辉石(韩发等,1980)相比,芦子园矿床的蔷薇辉石更富集MnSiO3而亏损CaSiO3分子(图 4),SiO2和FeO含量相近。1350中段蔷薇辉石(LZY26)的MnO含量略高于1220中段(LZY24)(表 3),与湖南柿竹园矽卡岩型钨矿床中蔷薇辉石的MnO含量从近接触带至远接触带略有增高的变化特点一致(毛景文等,1994)。根据蔷薇辉石成分平均含量计算其晶体化学式为(Ca0.1857Mn0.6703Fe0.0671Mg0.0077)Si1.0145O3,符合蔷薇辉石的结晶化学式中Ca不超过0.44,Mn不超过0.68,Fe不超过0.10的规律(韩发等,1980)。将本区蔷薇辉石按MnSiO3、FeSiO3、CaSiO3计算其分子百分比,投于MnSiO3-FeSiO3-CaSiO3三角图解上,投点均落入蔷薇辉石区内(图 4)。

图 4 MnSiO3-FeSiO3-CaSiO3成分关系图解(底图据Sundius,1931) Fig. 4 Component relationship diagram of MnSiO3-FeSiO3-CaSiO3(after Sundius,1931)
表 4 芦子园矿床蔷薇辉石、阳起石及闪锌矿单矿物微量元素含量(×10-6) Table 4 Trace element compositions of the rhodonite,actinolite and sphalerite from Luziyuan Pb-Zn deposit(×10-6)
4.3 微量元素组成

芦子园Pb-Zn-Fe矿床蔷薇辉石、闪锌矿及阳起石单矿物微量元素含量列于表 4。与原始地幔微量元素含量(Taylor and Mclennan,1985)相比的蛛网图见图 5

图 5 芦子园矿床单矿物原始地幔标准化微量元素蛛网图(标准化值据Taylor and Mclennan,1985) Fig. 5 Trace element enrichment coefficient in Luziyuan Pb-Zn deposit(normalization values after Taylor and Mclennan,1985)

表 4图 5可以看出,芦子园矿床中9件蔷薇辉石样品不相容元素Li、Be、Sc、Zr、Rb、Sr不亏损或弱亏损,相容元素Cr、Co、Ni明显亏损,这一特征可能暗示蔷薇辉石的形成与深部岩浆岩有关,在岩浆结晶分异过程中,相容元素优先进入岩浆分异早期的结晶矿物中,不相容元素留在残余岩浆中运移至浅部交代大理岩等碳酸盐岩,显示出相对富集的特征,这与前述矿区岩石主量元素变化特征相对应。与保山地块上燕山期漕涧花岗岩相对比,除成矿元素Zn含量高,亲石元素Rb偏低外,蔷薇辉石其他微量元素组成与漕涧花岗岩具有较好的一致性,进一步表明蔷薇辉石的形成与岩浆岩有关。

表 5 芦子园矿床蔷薇辉石、阳起石及闪锌矿稀土元素含量(×10-6)及特征参数 Table 5 REE compositions(×10-6)of the rhodonite,actinolite and sphalerite from Luziyuan Pb-Zn deposit

为了解微量元素在蔷薇辉石、阳起石和闪锌矿中赋存的异同点,本次还分析了4件闪锌矿和2件阳起石的微量元素。从蛛网图(图 5)中可以看出,芦子园矿床中闪锌矿的微量元素赋存特征与蔷薇辉石、阳起石表现出较为一致的微量元素分配模式,进一步表明闪锌矿和蔷薇辉石具有相似的物质来源,可能为同一成矿流体的产物。Cu和Co在蔷薇辉石和阳起石中的含量明显低于闪锌矿,与Zn具有良好的正相关关系,反映蔷薇辉石和阳起石的形成早于闪锌矿,早期高温流体交代大理岩形成石榴子石、透辉石以及蔷薇辉石等干矽卡岩阶段矿物,随着流体的演化,中晚期含矿流体交代大理岩及干矽卡岩阶段矿物形成阳起石、绿帘石、绿泥石等湿矽卡岩阶段矿物,由于温度、pH、Eh等物理化学条件的变化,金属硫化物开始沉淀。

4.4 稀土元素组成

芦子园矿床中蔷薇辉石、阳起石及闪锌矿的稀土元素含量见表 5。以Taylor and Mclennan(1985)的球粒陨石作为标准,作稀土元素配分曲线如图 6。从表 5图 6可以看出:

图 6 芦子园矿床单矿物及区域岩浆岩球粒陨石标准化稀土元素配分曲线(标准化值据Taylor and Mclennan,1985)志本山、柯街花岗岩据陶琰等,2010; 漕涧花岗岩据禹丽等,2014 Fig. 6 REE patterns of the single mineral seperates in Luziyuan deposit andthe granite in Baoshan Land mass(normalization values after Taylor and Mclennan,1985)

(1) 蔷薇辉石稀土元素配分模式可分为以下2种类型。团块状、角砾状蔷薇辉石(LZY11-4、LZY13-1、LZY24、LZY40)表现为平坦型(图 6a),∑REE为1.33×10-6~3.75×10-6,平均2.29×10-6,LREE/HREE为1.53~3.70,平均2.22,轻重稀土分异程度很低,δEu为0.75~2.07,为弱负铕异常到正铕异常,δCe为0.89~1.06,几乎不显示铈异常。脉状蔷薇辉石(LZY11-5、LZY13-2、LZY21、LZY31)为右倾型(图 6b),∑REE为4.65×10-6~5.75×10-6,平均5.11×10-6,比团块状蔷薇辉石含量高,LREE/HREE为3.77~6.29,平均4.64,轻重稀土分异明显,δEu为0.99~1.53,表现为弱正铕异常,δCe为0.96~1.08,无明显铈异常。分析结果与前述蔷薇辉石地质特征相对应,表明芦子园铅锌矿床存在2种类型的蔷薇辉石:第一种以稀土元素配分模式呈平坦型的团块状、角砾状蔷薇辉石为主;第二种主要为脉状蔷薇辉石,其稀土配分模式为轻稀土富集右倾型。

(2) 2件阳起石样品稀土元素配分模式均为右倾型(图 6c)。∑REE为2.13×10-6~12.11×10-6,平均7.21×10-6,LREE/HREE为9.81~11.8,平均10.8,δEu为0.50~0.81,为铕负异常,δCe为0.90~0.94,几乎不显示铈异常。与蔷薇辉石相比,阳起石轻重稀土明显分异,铕从正异常向负异常转变,铈无明显变化,表明随蚀变作用的进行,深部热液流体对阳起石的叠加改造更强烈,使其稀土元素特征与蔷薇辉石明显不同,而与闪锌矿相似。

(3) 闪锌矿稀土元素配分模式为右倾型(图 6d)。∑REE为0.41×10-6~10.72×10-6,平均3.43×10-6,LREE/HREE为5.40~40.2,平均17.2,δEu为0.41~0.75,为铕负异常,δCe为0.60~1.03,表现为负铈异常。闪锌矿轻重稀土分馏程度比阳起石及蔷薇辉石高,具明显铕负异常。

(4) 总体上,本次测试样品轻稀土相对富集,LREE/HREE均大于1,随着蚀变作用的进行,轻重稀土分异程度升高,稀土元素配分模式逐渐向区域岩浆岩的配分模式过渡(图 6)。

4.5 蔷薇辉石H、O同位素组成

芦子园矿床蔷薇辉石H、O同位素测试结果见表 6,根据公式1000lnαRhodonite-Water=3.85×106/T2+(-8.71)×103/T+2.46计算获得δ18OH2O值(Zheng,1993)。蔷薇辉石流体包裹体中的δDV-SMOW为-81.7‰~-73.3‰,δ18OH2O为5.4‰~6.1‰。在δD-δ18OH2O同位素图解上(图 7),样品投点均落入原生岩浆水左下侧。

表 6 芦子园矿床蔷薇辉石H、O同位素分析结果 Table 6 Hydrogen and Oxygen isotopic compositions of the rhodonite in Luziyuan Pb-Zn deposit
5 讨论 5.1 成矿物质来源 5.1.1 稀土元素示踪

因各个稀土元素之间具有相近的地球化学行为,常作为一个整体参与地质地球化学演化过程(龙汉生等,2011)。稀土元素地球化学研究在矿床成因(Schade et al.,1989)及热液体系中的水/岩反应机理(Bau,1991; Klinkhammer et al.,1994)等方面得到了广泛应用。REE3+的离子半径在0.086~0.103之间(8次配位)(王中刚等,1989),与Ca2+的离子半径(0.099,6~8次配位)(牟保磊,1999)较为相似,REE3+容易以类质同象的形式进入岩石晶格中替代岩石中的Ca2+等阳离子成分;与Zn2+的离子半径(0.074,4~6次配位)(刘英俊等,1984)差异较大,REE3+很难以类质同象的形式替代Zn2+进入矿物晶格,表明闪锌矿主要是在矿质沉淀的同时继承了成矿流体中的REE,其REE组成不容易被后期地质作用改造,更能代表成矿期流体的REE组成特征。

在芦子园矿床中,岩浆热液沿断裂构造上升运移的同时,与碳酸盐岩发生水岩反应,大量Mn2+、Fe2+、Mg2+等离子进入围岩中形成石榴子石、蔷薇辉石和阳起石等矽卡岩矿物。同时,REE以类质同象的形式进入碳酸盐岩晶格中,对其REE组成叠加改造,从而使蔷薇辉石的REE组成与该区岩浆岩差异较大。杨富全等(2007)对蒙库矽卡岩型铁矿床的研究认为:矽卡岩矿物的稀土组成与周边花岗岩差异较大时,表明其不是接触交代成因,这与矿区矽卡岩呈层状、似层状顺层产出而不产于岩体接触带的地质事实相符,而蔷薇辉石存在两种稀土配分模式与受岩浆热液改造的程度有关。随着深部岩浆热液上涌,成矿流体温度降低,有利于贱金属元素沉淀(李剑等,2005),此时,大规模的金属硫化物形成,闪锌矿的REE主要赋存于其包裹体或粒间裂隙中,其REE配分模式为轻稀土富集右倾型,具明显的Eu负异常,与周边燕山期志本山、柯街花岗岩(陶琰等,2010)以及漕涧花岗岩(禹丽等,2014)稀土配分模式相似(图 6),反映了矿区成矿物质主要来自深部岩浆的特征。

图 7 芦子园矿床蔷薇辉石δD-δ18OH2O图解(底图据Sheppard,1977) 图中保山核桃坪矽卡岩型铅锌多金属矿δD-δ18OH2O范围据薛传东等(2008) Fig. 7 Hydrogen and Oxygen isotopic diagram of the rhodonite from Luziyuan Pb-Zn deposit(after Sheppard,1977)

Y和Ho具有相同的价态和离子半径而被认为其地球化学行为相近,在许多地质过程中Y/Ho值保持不变(Nozaki et al.,1997)。近年来,Y/Ho在矿床成因(Krupenin,2004)、成矿流体来源(宋史刚等,2008; 章永梅等,2014)、构造背景演化(Huang et al.,2013)及成矿环境判别(黄虎等,2012)等方面得到了广泛应用。原始岩浆岩及碎屑物质(主要为玄武质和长英质地壳)Y/Ho值与球粒陨石相似,变化范围为24~34,平均28(Taylor and Mclennan,1985; Bau et al.,1996; Bolhar et al.,2005),而现代海水Y/Ho值较高,变化范围大,为44~74(Bau and Dulski,1999)。芦子园矿床矿石Y/Ho值为21.1~32.4,平均27.6(林冰霞,2013),蔷薇辉石Y/Ho值变化范围为22.0~33.2,平均28.4,阳起石样品Y/Ho值相对偏低,平均19.4,闪锌矿为24.7(表 5)。本文将该矿床矿石及各类矿物的Y/Ho值与周边燕山期岩浆岩作了对比(图 8),从图 8可以看出,芦子园铅锌矿区的Y/Ho值与志本山、漕涧花岗岩较为相似,这一特征进一步表明芦子园矿床的成矿物质主要来自深部岩浆岩。

图 8 芦子园铅锌矿区与区域花岗岩Y/Ho值对比 EPR(东太平洋洋脊)、MAR(中大西洋洋脊)、BAB弧后盆地和现代海水的Y/Ho值数据引自Douville et al.(1999)Bau et al.(1997)Bau and Dulski(1999);芦子园矿区矿石Y/Ho值林冰霞(2013);漕涧地区花岗岩Y/Ho禹丽等(2014);志本山花岗岩Y/Ho陶琰等(2010) Fig. 8 Compared the Y/Ho ratios of mineral from Luziyuan and regional granite
5.1.2 氢、氧同位素示踪

一种矿物的成因或起源常可由稳定同位素分析来推断(郑永飞等,2000)。水是成矿流体的重要组分,其H、O同位素组成是鉴别水来源的重要示踪剂(刘英超等,2009)。芦子园矿床蔷薇辉石流体包裹体中H、O同位素分析结果(表 6图 7)显示,蔷薇辉石的δ18OH2O变化范围较窄,δDV-SMOW变化范围相对较宽。在图 7中,蔷薇辉石H、O同位素投点均落于原生岩浆水左下侧,与保山核桃坪矽卡岩型Pb-Zn-Cu多金属矿床成矿期石英的δD-δ18OH2O(薛传东等,2008)组成较为接近,其H同位素未发生漂移可能是因为芦子园矿床形成早期成矿流体主要来自深部岩浆热液,尚未受到大气降水等的混入。对芦子园矿床蔷薇辉石包裹体研究表明,本区蔷薇辉石包裹体主要为气液两相包裹体,据拉曼特征峰值可知,液相成分主要为H2O,气相成分为CO2和CH4,均一温度变化范围为230~320℃,平均259℃,盐度变化范围为10.61%~17.60%NaCleqv)(未发表数据),为H2O-NaCl-CO2三元体系,其流体性质与张德会等(2011)对成矿热液系统分类中的岩浆流体类似。上述表明形成芦子园矿床蔷薇辉石的热液流体主要来自深部岩浆热液。

5.2 蔷薇辉石Eu正异常的形成及成矿环境

目前对Eu正异常的形成机理尚存较大争议,包括岩浆分异、水岩反应、岩石对流体中Eu2+的吸附作用以及对承载Eu2+络合物的吸附作用等,同时成岩成矿的氧化还原环境、温度和pH值等也是出现Eu异常的重要条件(Bau,1991; Klinkhammer et al.,1994; Hass et al.,1995; 杨富全等,2007; 任云生等,2010; 张承帅,2012)。芦子园矿床蔷薇辉石Eu以正异常为主,样品LZY11-4、LZY24显示Eu弱负异常,LZY11-5、LZY21未显示Eu异常,其他样品均显示Eu正异常,δEu为1.16~2.07,变化范围较小,Eu含量变化与Sm、Gd含量变化较为一致。Ghaderi et al.(1999)认为若Eu正异常以Eu3+的形式出现时,Eu与Sm、Gd的含量变化同步,且Eu异常的变化范围小;若Eu正异常以Eu2+的形式出现,因其离子半径较大,相对其它REE3+的地球化学行为是独立的,Eu异常强度取决于矿物中Sm、Gd的富集程度,二者呈负相关关系,因此认为蔷薇辉石中Eu正异常以Eu3+的形式出现。据Ce和Eu在不同氧化还原环境下的价态变化,Eu从正异常向负异常转变指示成矿环境向还原环境演化的过程,即Eu正异常越大,氧化程度越强,Eu负异常越大,还原程度越强,且高温环境是Eu正异常形成的重要条件(赵振华,1997; 赵劲松等,2007; 丁振举等,2013)。由表 5图 6可以看出,Eu异常从蔷薇辉石→阳起石→闪锌矿呈递减趋势,由正异常向负异常转变,表明蔷薇辉石形成于高温氧化环境,且成矿环境由氧化环境向还原环境转变。

矽卡岩的矿物组合也能为成矿环境判别提供有力证据,以钙铁榴石和富镁辉石矿物为主的矽卡岩指示氧化的环境,而钙铁辉石、钙铝榴石则主要形成于还原环境(Kwak,1994; Meinert et al.,2005);且石榴石中的Fe2+/Fe3+可以有效判别矽卡岩形成时的氧化还原环境,当Fe2+/Fe3+值较低时,主要为氧化环境;当Fe2+/Fe3+值较高时,以还原环境为主(Zaw and Singoyi,2000)。据电子探针分析结果,石榴子石的主要端元组分为钙铁榴石,含量高达84%,Fe2+/Fe3+值较低,为0.03~0.22,平均0.12;辉石中MgO含量高达9%,Mg/Fe值大多>2,顽火辉石端元组分达27%,为富Mg辉石(曾磊,2014),表明矿床矽卡岩矿物形成于氧化环境。同时在矽卡岩早期生成大量的氧化物(磁铁矿、赤铁矿)而几乎未有硫化矿生成也支持了上述观点。

结合对蔷薇辉石包裹体均一温度(230~320℃,平均259℃)和石英包裹体均一温度(108~250℃,平均136℃)(林冰霞,2013)的研究认为:芦子园矿床成矿环境从矽卡岩期高温氧化环境向石英-硫化物期低温还原环境转变,更有利于Pb、Zn等成矿元素的沉淀富集(张艳等,2015)。这也与芦子园矿床客观地质事实相吻合,早期在矿床下部的高温氧化环境下形成厚大的磁铁矿氧化物矿体,后期转变为低温还原条件并在浅部形成Pb-Zn(Cu)多金属硫化物矿体。

5.3 成因意义

尽管前人对芦子园Pb-Zn-Fe多金属矿床作了大量研究,但对矿床成因方面的认识还存在较大争议(邓必方,1995; 陈永清等,2005; 陶琰等,2010; 朱飞霖等,2011; 卢映祥,2011),其争论的焦点在于成矿物质和流体来源、地层与成矿关系的不确定。基于矿体赋存于具有较高Pb、Zn背景值的寒武系沙河厂组地层,产状与围岩一致,具一定的沉积韵律结构(如矿体呈层状、条带状、纹层状等),且矿区内尚未发现与成矿有关的岩浆岩体,仅见少量辉绿岩脉出露,加之赋矿围岩具有蚀变类型多样等基本地质特征,部分学者因此提出沉积改造成因的观点,认为该矿床是经早期沉积并受到后期动力变质、区域变质及岩浆热液的改造、富集(邓必方,1995; 李文昌和莫宣学,2001; 陈永清等,2005; 董文伟和陈少玲,2007)而成。而卢映祥(2011)认为该矿床为多期成矿作用叠加的产物,分别为早古生代寒武系地层中的Pb、Zn、Cu预富集成矿作用→与晚古生代、中生代基性岩浆岩侵入、喷发有关的铁铜成矿作用→与新生代(隐伏)酸性岩体有关的接触、交代变质成矿作用。也有学者认为所谓的"矿源层"应具有相应较大范围内的矿化和其他地质作用现象,而这些特征在芦子园矿区尚未得到进一步证实,因此矿区地层不一定是矿床成矿物质的来源(卢映祥,2011)。

因矿区赋矿围岩主要为热液蚀变形成的矽卡岩(图 1d),多数学者认为该矿床为沉积改造成因的可能性小。夏庆霖等(2005)通过流体包裹体研究认为芦子园矿区经历了中低温(160~280℃)和中高温(280~420℃)两个主要成矿阶段,测得的S-Pb同位素组成变化范围窄、相对均一,显示成矿作用曾受岩浆活动的影响,初步认为芦子园Pb-Zn-Fe多金属矿床为与上寒武统局部地层和隐伏岩体有关的热液矿床。朱飞霖等(2011)通过闪锌矿Rb-Sr测年获得等时线年龄为141.9±2.6Ma,Sr同位素初始比值(0.714497)与核桃坪矽卡岩型铅锌矿床(0.711850)(陶琰等,2010)相似,进一步认为芦子园矿床为与燕山期隐伏花岗岩有关的岩浆热液矿床,这与杨玉龙等(2012)通过S、Pb、H、O同位素所获得结果相一致。本文对蔷薇辉石的H、O同位素及微量稀土研究表明,芦子园矿区成矿物质主要来源于深部岩浆岩,这与地质事实及以往认识相吻合。

锰质矽卡岩总是产于离侵入体有一定距离的碳酸盐岩中,多与中酸性岩浆岩有关,并和铅锌矿化关系密切(赵一鸣等,19862002)。与锰质矽卡岩有关的代表性矿床有辽宁八家子(夏学惠等,1999)、湖南柿竹园(毛景文等,1994)、福建大排和马坑(丰成友等,2006)等Pb-Zn-Fe多金属矿床,该类矿床具有以下共同特征:①矿体主要产于与中酸性中小岩体有一定距离的碳酸盐岩过渡带,容矿或控矿构造常为层间断裂破碎带或推覆构造;②矿体以层状或似层状为主,矿化受层状矽卡岩控制;③围岩蚀变强烈,以锰质矽卡岩化为主,常伴有绿泥石化和黄铁矿化等;④除铅锌矿化外,常伴有铜、银、磁铁矿化,在空间上显示出一定的蚀变矿化分带特征。 芦子园铅锌矿床矿体均呈层状、似层状或透镜状沿NE向层间断裂顺层产出,在空间上与矽卡岩密切伴生,由深部至浅部具磁铁矿化→磁铁、铅锌(铜)矿化→铅锌矿化的分带特征(图 1d)。矿区围岩蚀变强烈,主要为矽卡岩化,蚀变岩锰质含量高,以蔷薇辉石为代表,具典型的锰质矽卡岩建造;赋矿围岩主要为矽卡岩和大理岩,在蔷薇辉石产出部位矿体规模变大,品位较富(邓明国等,2013);矿石主要为热液充填块状、脉状及热液交代浸染状矿石(图 2c-e),具明显的热液成矿特点,符合远源带特征(Meinert,1997)。李开毕等(2012)通过高精度磁测推断芦子园矿区深部存在隐伏的中酸性岩体;吾守艾力·肉孜等(2015)Liang et al.(2015)通过航磁、重力及地质资料综合研究认为芦子园矿区存在隐伏酸性岩体,并通过AMT二维反演出矿区深部地质构造及隐伏岩体的埋深及形态。隐伏岩体可能为矿区大规模矽卡岩化及铅锌矿化提供流体、热源及成矿物质。上述表明,芦子园Pb-Zn-Fe多金属矿床具有典型远程矽卡岩矿床的特征(赵一鸣等,198619972002)。

综上认为,芦子园矿床为与燕山期隐伏中酸性岩浆岩有关的远程矽卡岩型Pb-Zn-Fe多金属矿床。矿床矿体受NE向层间断裂控制明显,矿石矿物及围岩蚀变等显示出明显的热液成因,H、O同位素及微量稀土数据均表明成矿流体来自深部岩浆岩。基于上述认识,未来在芦子园矿区深部及外围的勘查找矿工作中,应将矽卡岩型铁矿及与区域上中酸性岩浆成矿系统有关的矿化元素(如Cu、Mo、Sn和Hg等)列为重点勘查对象。

6 结论

(1) 芦子园Pb-Zn-Fe多金属矿床中深部发现的蔷薇辉石具有二种产出特征:第一种蔷薇辉石主要为团块状、角砾状,稀土配分模式为平坦型;第二种蔷薇辉石为脉状,稀土配分模式为轻稀土富集右倾型。二种蔷薇辉石在空间上与铅锌矿紧密共(伴)生,在该区及该类型矿床中具重要的找矿指示意义及矿床成因意义。其H、O同位素投点在δD-δ18OH2O图解中均落入岩浆水区域,显示蔷薇辉石为岩浆热液成因的特征。

(2) Eu异常由蔷薇辉石→闪锌矿的变化特征,结合蚀变矿物组合特点及包裹体均一温度结果,认为矿床成矿环境为一高温氧化环境向低温还原环境转化的过程,更有利于后期金属硫化物的沉淀。

(3) H、O同位素、微量、稀土分析结果均表明芦子园矿区成矿物质主要来自深部岩浆岩,综合矿床地质特征及前人研究结果认为芦子园矿床为与燕山期隐伏中酸性岩浆岩有关的远程矽卡岩型Pb-Zn-Fe多金属矿床。

致谢 本研究相关野外地质调查及样品采集工作得到镇康县鸿骏矿业开发有限公司领导和云南省地质调查院驻镇康芦子园项目组工作人员的大力帮助和支持;匿名审稿人和编辑对论文的修改完善提出了宝贵意见;谨此致谢!
参考文献
[1] Bau M. 1991. Rare-earth element mobility during hydrothermal and metamorphic fluid-rock interaction and the significance of the oxidation state of europium. Chemical Geology , 93 (3-4) :219–230. DOI:10.1016/0009-2541(91)90115-8
[2] Bau M, Koschinsky A, Dulski P, Hein JR. 1996. Comparison of the partitioning behaviours of yttrium, rare earth elements, and titanium between hydrogenetic marine ferromanganese crusts and seawater. Geochimica et Cosmochimica Acta , 60 (10) :1709–1725. DOI:10.1016/0016-7037(96)00063-4
[3] Bau M, M ller P, Dulski P. 1997. Yttrium and lanthanides in eastern Mediterranean seawater and their fractionation during redox-cycling. Marine Chemistry , 56 (1-2) :123–131. DOI:10.1016/S0304-4203(96)00091-6
[4] Bau M, Dulski P. 1999. Comparing yttrium and rare earths in hydrothermal fluids from the Mid-Atlantic Ridge:Implications for Y and REE behaviour during near-vent mixing and for the Y/Ho ratio of Proterozoic seawater. Chemical Geology , 155 (1-2) :77–90. DOI:10.1016/S0009-2541(98)00142-9
[5] Bolhar R, Van Kranendonl MJ, Kamber BS. 2005. A trace element study of siderite-jasper banded iron formation in the 3. 45Ga Warrawoona Group, Pilbara Graton-formation from hydrothermal fluids and shallow seawater. Precambrian Research , 137 (1-2) :93–114.
[6] Chen YQ, Lu YX, Xia QL, Jiang CX, Liu HG, Lü ZC. 2005. Geochemical characteristics of the Hetaoping Pb-Zn deposit, Baoshan, Yunnan, and its genetic model and ore prospecting model pattern. Geology in China , 32 (1) :90–99.
[7] Deng BF. 1995. The metallogenic model of Hg-Pb-Zn in Zhenkang area, Baoshan, Yunnan. Yunnan Geology , 14 (4) :355–364.
[8] Deng J, Wang QF, Li GJ, Santosh M. 2014a. Cenozoic tectono-magmatic and metallogenic processes in the Sanjiang region, southwestern China. Earth-Science Reviews , 138 :268–299. DOI:10.1016/j.earscirev.2014.05.015
[9] Deng J, Wang QF, Li GJ, Li CS, Wang CM. 2014b. Tethys tectonic evolution and its bearing on the distribution of important mineral deposits in the Sanjiang region, SW China. Gondwana Research , 26 (2) :419–437. DOI:10.1016/j.gr.2013.08.002
[10] Deng J, Wang CM, Li WC, Yang LQ, Wang QF. 2014. The situation and enlightenment of the research of the tectonic evolution and metallogenesis in the Sanjiang Tethys. Earth Science Frontiers , 21 (1) :52–64.
[11] Deng MG, Li WC, Wen HJ, Cang FB, Wang P, Lü CL, Xu R, Zeng L, Liu XL, Liu W, Yin GH. 2013. The finding and prospecting significance of rhodonite in the Luziyuan lead-zinc polymetallic ore deposit in Zhenkang, western Yunnan. Geological Bulletin of China , 32 (11) :1867–1869.
[12] Ding ZJ, Yao SZ, Liu CQ, Zhou ZG, Yang MG. 2013. The characteristics of exhalation-sedimentary deposit of Donggouba polymetal deposit:Evidence from ore's REE composition. Acta Petrologica Sinica , 19 (4) :792–798.
[13] Dong ML, Dong GC, Mo XX, Zhu DC, Nie F, Xie XF, Wang X, Hu ZC. 2012. Geochronology and geochemistry of the Early Palaeozoic granitoids in Baoshan block, western Yunnan and their implications. Acta Petrologica Sinica , 28 (5) :1453–1464.
[14] Dong WW, Chen SL. 2007. The characteristics & genesis of Luziyuan Pb-Zn deposit, Zhenkang. Yunnan Geology , 26 (4) :404–410.
[15] Douville E, Bienvenu P, Charlou JI, Donval JP, Fouquet Y, Appriou P, Gamo T. 1999. Yttrium and rare earth elements in fluids from various deep-sea hydrothermal systems. Geochimica et Cosmochimica Acta , 63 (5) :627–643. DOI:10.1016/S0016-7037(99)00024-1
[16] Feng CY, Zhang DQ, She HQ, Li DX, Li JW, Wu JS. 2006. Major types of lead-zinc deposits in Fujian Province in relation to manganese skarnization. Mineral Deposits , 25 (Suppl.1) :341–344.
[17] Ghaderi M, Palin M, Campbell IH, Sylvester PJ. 1999. Rare earth element systematics in scheelite from hydrothermal gold deposits in the Kalgoorlie-Norseman region, Western Australia. Economic Geology , 94 (3) :423–437. DOI:10.2113/gsecongeo.94.3.423
[18] Gu YQ, Qan TH, Ye ZS, Yan YB. 1988. The Petrographic and geochemical characteristics of the Muchang A-type granite in Zhenkang, western Yunnan. Acta Petrologica et Mineralogica , 7 (1) :38–48.
[19] Han F, Ge CH, Tan HJ, Yuan RG, Zhao YM. 1980. The rhodonite and manganese clinopyroxene of Makeng, south of Fujian Provence. Geological Reviews , 26 (3) :208–216.
[20] Hass JR, Shock EL, Sassani DC. 1995. Rare earth elements in hydrothermal systems:Estimates of standard partial molal thermodynamic properties of aqueous complexes of the rare earth elements at high pressures and temperatures. Geochimica et Cosmochimica Acta , 59 (21) :4329–4350. DOI:10.1016/0016-7037(95)00314-P
[21] Huang H, Du YS, Yang JH, Tao P, Huang HW, Huang ZQ, Xie CX, Hu LS. 2012. Geochemical features of siliceous sediments of the Shuicheng-Ziyun-Nandan rift basin in the Late Paleozoic and their tectonic implication. Acta Geologica Sinica , 86 (12) :1994–2010.
[22] Huang H, Du YS, Huang ZQ, Yang JH, Huang HW, Xie CX, Hu LS. 2013. Depositional chemistry of chert during late Paleozoic from western Guangxi and its implication for the tectonic evolution of the Youjiang Basin. Science China (Earth Sciences) , 56 (3) :479–493. DOI:10.1007/s11430-012-4496-y
[23] Klinkhammer GP, Elderfield H, Edmond JM, Mitra A. 1994. Geochemical implications of rare earth element patterns in hydrothermal fluids from mid-ocean ridges. Geochimica et Cosmochimica Acta , 58 (23) :5105–5113. DOI:10.1016/0016-7037(94)90297-6
[24] Krupenin MT. 2004. Y/Ho ratio as genetic indicator of sparry magnesites in South Urals, Russia. Acta Petrologica Sinica , 20 (4) :803–816.
[25] Kwak TAP. 1994. Hydrothermal alteration in carbonate-replacement deposits:Ore skarns and distal equivalents. In:Lentz DR (ed.). Alteration and Alteration Processes Associated with Ore-Forming Systems. Geological Association of Canada, Short Course Notes , 11 :381–402.
[26] Li GJ, Wang QF, Huang YH, Chen FC, Dong P. 2015. Discovery of Hadean-Mesoarchean crustal materials in the northern Sibumasu block and its significance for Gondwana reconstruction. Precambrian Research , 271 :118–137. DOI:10.1016/j.precamres.2015.10.003
[27] Li GJ, Wang QF, Huang YH, Gao L, Yu L. 2016. Petrogenesis of middle Ordovician peraluminous granites in the Baoshan block:Implications for the Early Paleozoic tectonic evolution along East Gondwana. Lithos , 245 :76–92. DOI:10.1016/j.lithos.2015.10.012
[28] Li J, Lü ZC, Dong FL. 2005. Ore-forming fluids related to the intermediate-acid magmatism system and their trace element geochemical characteristics. Geology in China , 32 (4) :614–624.
[29] Li KB, Yang SS, Cai X. 2012. The effect and result of high-precision magnetic survey on prospecting in Luziyuan mine in Zhenkang Country. Journal of Yunnan University (Natural Sciences) , 34 (Suppl.2) :157–162.
[30] Li WC, Mo XX. 2001. The Cenozoic tectonics and metallogenesis in the "Three-River" area of Southwest China. Yunnan Geology , 14 (4) :333–346.
[31] Liang Q, Jing H, Gregoire DC. 2000. Determination of trace elements in granites by inductively coupled plasma mass spectrometry. Talanta , 51 (3) :507–513. DOI:10.1016/S0039-9140(99)00318-5
[32] Liang SX, Jiao YJ, Guo J. 2015. Prediction of hidden granites in the Luziyuan area of Yunnan Province and the prospecting direction. Acta Geologica Sinica , 89 (5) :1781–1782. DOI:10.1111/acgs.2015.89.issue-5
[33] Lin BX. 2013. Research on geological and geochemical features of Pb-Zn deposit in Luziyuan, Zhenkang, Yunnan Province. Master Degree Thesis. Kunming:Kunming University of Science and Technology (in Chinese with English summary)
[34] Liu HB, Jin GS, Li JJ, Han J, Zhang JF, Zhang J, Zhong FW, Guo DQ. 2013. Determination of stable isotope composition in uranium geological samples. World Nuclear Geoscience , 30 (3) :174–179.
[35] Liu YC, Yang ZS, Hou ZQ, Tian SH, Wang ZL, Song YC, Xue WW, Lu HF, Wang FC, Zhang YB, Zhu T, Yu CJ, Su AN, Li ZZ, Yu YS. 2009. Geology and hydrogen, oxygen and carbon isotope geochemistry of Dongmozhazhua Pb-Zn ore deposit, Yushu area, Qinghai Province. Mineral Deposits , 28 (6) :770–784.
[36] Liu YJ, Cao LM, Li ZL, Wang HN, Chu TQ, Zhang JR.1984. Element Geochemistry. Beijing: Science Press .
[37] Long HS, Luo TY, Huang ZL, Zhou MZ, Yang Y, Qian ZK. 2011. Rare Earth element and trace element geochemistry of pyrite ores in the Laochang large size silver ploymetallic deposit of Lancang, Yunnan Province, China. Acta Mineralogica Sinica , 31 (3) :462–473.
[38] Lu YX. 2011. Pb-Zn-Cu-Fe superimposed metallogenesis and regional ore-prospecting model in Baoshan-Zhenkang continent block. In:Academic-Salon Collected Works for Neodoxy 55:The Significant Influence of Plates Convergence and Mantle Plume on Yunnan Regional Metall-Ogenesis. Kunming:China Association of Science and Technology Academic Department, 78-83 (in Chinese)
[39] Mao JW, Li YH, Wang AP. 1994. Manganoan skarn in the Shizhuyuan polymetallic tungsten deposit, Hunan Province. Mineral Deposits , 13 (1) :38–47.
[40] Meinert LD. 1997. Application of skarn deposit zonation models to mineral exploration. Exploration and Mining Geology , 6 (2) :185–208.
[41] Meinert LD, Dipple GM, Nicolescu S. 2005. World skarn deposits , 100 :299–236.
[42] Metcalfe I. 2011. Tectonic framework and phanerozoic evolution of Sundaland. Gondwana Research , 19 (1) :3–21. DOI:10.1016/j.gr.2010.02.016
[43] Mo XX, Pan GT. 2006. From the Tethys to the formation of the Qinghai-Tibet Plateau:Constrained by tectono-magmatic events. Earth Science Frontiers , 13 (6) :43–51.
[44] Mou BL.1999. Element Geochemistry. Beijing: Peking University Press .
[45] Nozaki Y, Zhang J, Amakawa H. 1997. The fractionation between Y and Ho in the marine environment. Earth and Planetary Science Letters , 148 (1-2) :329–340. DOI:10.1016/S0012-821X(97)00034-4
[46] Ren YS, Zhao HL, Lei E, Wang H, Ju N, Wu CZ. 2010. Trace element and rare earth element geochemistry of the scheelite and ore genesis of the Yangjingou large scheelite deposit in Yanbian area, northeastern China. Acta Petrologica Sinica , 26 (12) :3720–3726.
[47] Rouzi WSAL, Liang SX, Zou GF, Luo MJ. 2015. The application of AMT and 3D gravity methods to the prospecting for concealed rock body in the Luziyuan area. Geophysical and Geochemical Exploration , 39 (3) :525–529.
[48] Schade J, Cornell DH, Theart HFJ. 1989. Rare earth element and isotopic evidence for the genesis of the Prieska massive sulfide deposit, South Africa. Economic Geology , 84 (1) :49–63. DOI:10.2113/gsecongeo.84.1.49
[49] Sheppard SMF. 1977. The Cornubian batholith, SW England:D/H and 18O/16O studies of kaolinite and other alteration minerals. Journal of the Geological Society , 133 (6) :573–591. DOI:10.1144/gsjgs.133.6.0573
[50] Song SG, Ding ZJ, Yao SZ, Zhou ZG. 2008. REE and Nd-Sr Isotopic compositions of the iron-rich siliceous rock in Bikou terrane:Implication to ancient sedimentary environment. Journal of Mineralogy and Petrology , 28 (3) :57–63.
[51] Sundius N. 1931. On the triclinic manganiferous pyroxenes. Am. Mineralogist, 16:411-429, 488-518 http://cn.bing.com/academic/profile?id=2145919621&encoded=0&v=paper_preview&mkt=zh-cn
[52] Tao Y, Hu RZ, Zhu FL, Ma YS, Ye L, Cheng ZT. 2010. Ore-forming age and the geodynamic background of the Hetaoping lead-zinc deposit in Baoshan, Yunnan. Acta Petrologica Sinica , 26 (6) :1760–1772.
[53] Taylor SR and Mclennan SM. 1985. The Continental Crust:Its Composition and Evolution. Oxford:Blackwell Scientific Publication http://www.oalib.com/references/19233886
[54] Wang CM, Bagas L, Lu YJ, Santosh M, Du B, McCuaig TC. 2016. Terrane boundary and spatio-temporal distribution of ore deposits in the Sanjiang Tethyan Orogen:Insights from zircon Hf-isotopic mapping. Earth-Science Reviews , 156 :39–65. DOI:10.1016/j.earscirev.2016.02.008
[55] Wang XF, Metcalfe I, Jian P, He LQ, Wang CS. 2000. The Jinshajiang-Ailaoshan suture zone, China:Tectonostratigraphy, age and evolution. Journal of Asian Earth Sciences , 18 (6) :675–690. DOI:10.1016/S1367-9120(00)00039-0
[56] Wang YJ, Xing XW, Cawood PA, Lai SC, Xia XP, Fan WM, Liu HC, Zhang FF. 2013. Petrogenesis of Early Paleozoic peraluminous granite in the Sibumasu Block of SW Yunnan and diachronous accretionary orogenesis along the northern margin of Gondwana. Lithos , 182-183 :67–85. DOI:10.1016/j.lithos.2013.09.010
[57] Wang ZG, Yu XY, Zhao ZH.1989. REE Geochemistry. Beijing: Science Press .
[58] Xia QL, Chen YQ, Lu YX, Jiang CX, Liu HG, Lü ZC. 2005. Geochemistry, fluid inclusion, and stable isotope studies of Luziyuan Pb-Zn deposit in Yunnan Province, southwestern China. Earth Science , 30 (2) :177–186.
[59] Xia XH, Liu CT, Yan F. 1999. The discovery of pyroxmangite in Bajiazi area of Liaoning Province and its ore_hunting implication. Acta Petrologica et Mineralogica , 18 (2) :150–156.
[60] Xue CD, Han RS, Yang HL, Yang ZM, Tian SH, Liu YQ, Hao BW. 2008. Isotope geochemical evidence for ore-forming fluid resources in Hetaoping Pb-Zn deposit, Baoshan, northwestern Yunnan. Mineral Deposits , 27 (2) :243–252.
[61] Yang FQ, Mao JW, Xu LG, Zhang y, Liu F, Huang CL, Zhou G, Liu GR, Dai ZJ. 2007. REE geochemistry of the Mengku iron deposit, Xinjiang, and its indication for iron mineralization. Acta Petrologica Sinica , 23 (10) :2443–2456.
[62] Yang YL, Ye L, Cheng ZT, Bao T, Gao W. 2012. A tentative discussion on the genesis of skarn Pb-Zn deposits in the Baoshan-Zhenkang terrane. Acta Petrologica et Mineralogica , 31 (4) :554–564.
[63] Yu L, Li GJ, Wang QF, Liu XF. 2014. Petrogenesis and tectonic significance of the Late Cretaceous magmatism in the northern part of the Baoshan block:Constraints from bulk geochemistry, zircon U-Pb geochronology and Hf isotopic compositions. Acta Petrologica Sinica , 30 (9) :2709–2724.
[64] Zaw K, Singoyi B. 2000. Formation of magnetite-scheelite skarn mineralization at Kara, northwestern Tasmania:Evidence from mineral chemistry and stable isotopes. Economic Geology , 95 (6) :1215–1230. DOI:10.2113/gsecongeo.95.6.1215
[65] Zeng L. 2014. Geochemical characteristics of Skran minerals and genesis research of Luziyuan Pb-Zn deposit, Wester Yunnan. Master Degree Thesis. Kunming:Kunming University of Science and Technology, 55-62 (in Chinese with English summary)
[66] Zhang CS.2012. Geology and geochemistry of Makeng Fe-Mo deposit Fujian. Ph. D. Dissertation. Beijing: China University of Geosciences .
[67] Zhang DH, Jin XD, Mao SD, Wang LL. 2011. The classification of ore-forming fluid and the efficiency of ore formation of magmatic hydrothermal solution. Earth Science Frontiers , 18 (5) :90–102.
[68] Zhang Y, Han RS, Wei PT, Qiu WL. 2015. pH-logfO2 and pH-loga for Pb-Zn paragenesis and separation in the Zhaotong lead-zinc deposit. Geology in China , 42 (2) :607–620.
[69] Zhang YM, Gu XX, Peng YW, Zheng L, Zhang Y, Gao HJ, Dong SY. 2014. Geology, geochemistry and genesis of the Bonga carbonatite-type niobium deposit, Angola. Earth Science Frontiers , 21 (5) :50–68.
[70] Zhao JS, Qiu XL, Zhao B, Tu XL, Yu J, Lu TS. 2007. REE geochemistry of mineralized skarns from Daye to Wushan region, China. Geochimica , 36 (4) :400–412.
[71] Zhao YM, Lin WW, Bi CS, Li DX. 1986. Basic geological characteristics of Skarn deposits of China. Bulletin of the Chinese Academy of Geological Sciences (14) :59–87.
[72] Zhao YM, Zhang YN, Lin WW. 1997. Characteristics of pyroxenes and pyroxenoids in Skarn deposits of China and their relationship with metallization. Mineral Deposits , 16 (4) :318–329.
[73] Zhao YM, Li DX, Dong YG. 2002. Further discussion on Pb-Zn-Ag-bearing manganoan skarn formation. Mineral Deposits , 21 (S1) :548–552.
[74] Zhao YM, Li DX. 2003. Amphiboles in skarn deposits of China. Mineral Deposits , 22 (4) :345–359.
[75] Zhao ZH.1997. Trace Element Geochemistry. Beijing: Science Press .
[76] Zheng YF. 1993. Calculation of oxygen isotope fractionation in anhydrous silicate minerals. Geochimica et Cosmochimica Acta , 57 (5) :1079–1091. DOI:10.1016/0016-7037(93)90042-U
[77] Zheng YF, Xu BL, Zhou GT. 2000. Geochemical studies of stable isotopes minerals. Earth Science Frontiers , 7 (2) :299–320.
[78] Zhu FL, Tao Y, Hu RZ, Liao MY, Wang YX, Li YB. 2011. Formation age of the Luziyuan lead-zinc deposit in Zhengkang, Yunnan. Bulletin of Mineralogy, Petrology and Geochemistry , 30 (1) :73–79.
[79] 陈永清, 卢映祥, 夏庆霖, 蒋成兴, 刘红光, 吕志成.2005. 云南保山核桃坪铅锌矿床地球化学特征及其成矿模式与找矿模型. 中国地质 , 32 (1) :90–99.
[80] 邓必方.1995. 保山-镇康地区汞、铅锌矿床的成矿模式. 云南地质 , 14 (4) :355–364.
[81] 邓军, 王长明, 李文昌, 杨力强, 王庆飞.2014. 三江特提斯复合造山与成矿作用研究态势及启示. 地学前缘 , 21 (1) :52–64.
[82] 邓明国, 李文昌, 温汉捷, 仓福彬, 王朋, 吕昶良, 徐荣, 曾磊, 刘学龙, 刘伟, 尹光侯.2013. 滇西镇康芦子园铅锌多金属矿床中蔷薇辉石的发现及找矿意义. 地质通报 , 32 (11) :1867–1869.
[83] 丁振举, 姚书振, 刘丛强, 周宗桂, 杨明国.2013. 东沟坝多金属矿床喷流沉积成矿特征的稀土元素地球化学示踪. 岩石学报 , 19 (4) :792–798.
[84] 董美玲, 董国臣, 莫宣学, 朱弟成, 聂飞, 谢许峰, 王霞, 胡兆初.2012. 滇西保山地块早古生代花岗岩类的年代学、地球化学及意义. 岩石学报 , 28 (5) :1453–1464.
[85] 董文伟, 陈少玲.2007. 镇康芦子园铅锌矿床特征及成因. 云南地质 , 26 (4) :404–410.
[86] 丰成友, 张德全, 佘宏全, 李大新, 李进文, 吴建设.2006. 福建省主要铅锌矿床类型及其与锰质矽卡岩化的关系. 矿床地质 , 25 (S1) :341–344.
[87] 顾影渠, 钱天宏, 叶芝珊, 颜以彬.1988. 滇西镇康木厂A型花岗岩岩石学及地球化学特征. 岩石矿物学杂志 , 7 (1) :38–48.
[88] 韩发, 葛朝华, 谭惠静, 袁润广, 赵一鸣.1980. 闽南马坑等地钙蔷薇辉石及锰三斜辉石. 地质论评 , 26 (3) :208–216.
[89] 黄虎, 杜远生, 杨江海, 陶平, 黄宏伟, 黄志强, 谢春霞, 胡丽沙.2012. 水城-紫云-南丹裂陷盆地晚古生代硅质沉积物地球化学特征及其地质意义. 地质学报 , 86 (12) :1994–2010.
[90] 李剑, 吕志成, 董方浏.2005. 中酸性岩浆体系成矿流体及微量元素地球化学特征. 中国地质 , 32 (4) :614–624.
[91] 李开毕, 杨淑胜, 蔡旭.2012. 高精度磁测在镇康芦子园矿区勘查中的作用及效果. 云南大学学报(自然科学版) , 34 (S2) :157–162.
[92] 李文昌, 莫宣学.2001. 西南"三江"地区新生代构造及其成矿作用. 云南地质 , 20 (4) :333–346.
[93] 林冰霞. 2013. 云南省镇康芦子园铅锌矿床地质地球化学特征研究. 硕士学位论文. 昆明:昆明理工大学 http://cdmd.cnki.com.cn/article/cdmd-10674-1014177293.htm
[94] 刘汉彬, 金贵善, 李军杰, 韩娟, 张建锋, 张佳, 钟芳文, 郭东侨.2013. 铀矿地质样品的稳定同位素组成测试方法. 世界核地质科学 , 30 (3) :174–179.
[95] 刘英超, 杨竹森, 侯增谦, 田世洪, 王召林, 宋玉财, 薛万文, 鲁海峰, 王富春, 张玉宝, 朱田, 俞长捷, 苏嫒娜, 李真真, 于玉帅.2009. 青海玉树东莫扎抓铅锌矿床地质特征及碳氢氧同位素地球化学研究. 矿床地质 , 28 (6) :770–784.
[96] 刘英俊, 曹励明, 李兆麟, 王鹤年, 储同庆, 张景荣. 1984. 元素地球化学. 北京: 科学出版社 .
[97] 龙汉生, 罗泰义, 黄智龙, 周明忠, 杨勇, 钱志宽.2011. 云南澜沧老厂大型银多金属矿床黄铁矿稀土和微量元素地球化学. 矿物学报 , 31 (3) :462–473.
[98] 卢映祥. 2011. 保山-镇康陆块铅锌铜铁叠加成矿作用与区域找矿模型. 见:新观点新学说学术沙龙文集55:板块汇聚、地幔柱对云南区域成矿作用的重大影响. 昆明:中国科学技术协会学会学术部, 78-83
[99] 毛景文, 李红艳, 王平安.1994. 湖南柿竹园钨多金属矿床中的锰质矽卡岩. 矿床地质 , 13 (1) :38–47.
[100] 莫宣学, 潘桂棠.2006. 从特提斯到青藏高原形成:构造-岩浆事件的约束. 地学前缘 , 13 (6) :43–51.
[101] 牟保磊. 1999. 元素地球化学. 北京: 北京大学出版社 .
[102] 任云生, 赵华雷, 雷恩, 王辉, 鞠楠, 吴昌志.2010. 延边杨金沟大型钨矿床白钨矿的微量和稀土元素地球化学特征与矿床成因. 岩石学报 , 26 (12) :3720–3726.
[103] 吾守艾力·肉孜, 梁生贤, 邹光富, 罗茂金.2015. AMT与重力方法在云南芦子园地区隐伏岩体勘查中的应用. 物探与化探 , 39 (3) :525–529.
[104] 宋史刚, 丁振举, 姚书振, 周宗桂.2008. 碧口地块富铁硅岩REE及Nd-Sr同位素组成及其古环境意义. 矿物岩石 , 28 (3) :57–63.
[105] 陶琰, 胡瑞忠, 朱飞霖, 马言胜, 叶霖, 程增涛.2010. 云南保山核桃坪铅锌矿成矿年龄及动力学背景分析. 岩石学报 , 26 (6) :1760–1772.
[106] 王中刚, 于学元, 赵振华. 1989. 稀土元素地球化学. 北京: 科学出版社 .
[107] 夏庆霖, 陈永清, 卢映祥, 蒋成兴, 刘红光, 吕志成.2005. 云南芦子园铅锌矿床地球化学、流体包裹体及稳定同位素特征. 地球科学 , 30 (2) :177–186.
[108] 夏学惠, 刘昌涛, 闫飞.1999. 辽宁八家子三斜锰辉石的发现及其找矿意义. 岩石矿物学杂志 , 18 (2) :150–156.
[109] 薛传东, 韩润生, 杨海林, 杨志明, 田世洪, 刘勇强, 郝百武.2008. 滇西北保山核桃坪铅锌矿床成矿流体来源的同位素地球化学证据. 矿床地质 , 27 (2) :243–252.
[110] 杨富全, 毛景文, 徐林刚, 张岩, 刘锋, 黄成林, 周刚, 刘国仁, 代军治.2007. 新疆蒙库铁矿床稀土元素地球化学及对铁成矿作用的指示. 岩石学报 , 23 (10) :2443–2456.
[111] 杨玉龙, 叶霖, 程增涛, 鲍谈, 高伟.2012. 保山镇康地块矽卡岩型铅锌矿床成因初探. 岩石矿物学杂志 , 31 (4) :554–564.
[112] 禹丽, 李龚健, 王庆飞, 刘学飞.2014. 保山地块北部晚白垩世岩浆岩成因及其构造指示:全岩地球化学、锆石U-Pb年代学和Hf同位素制约. 岩石学报 , 30 (9) :2709–2724.
[113] 曾磊. 2014. 滇西芦子园铅锌矿床矽卡岩矿物地球化学特征及矿床成因研究. 硕士学位论文. 昆明:昆明理工大学, 55-62 http://cdmd.cnki.com.cn/article/cdmd-10674-1015547950.htm
[114] 张承帅. 2012. 福建马坑铁矿床地质与地球化学. 博士学位论文. 北京:中国地质大学 http://cdmd.cnki.com.cn/article/cdmd-11415-1012364394.htm
[115] 张德会, 金旭东, 毛世德, 王丽丽.2011. 成矿热液分类兼论岩浆热液的成矿效率. 地学前缘 , 18 (5) :90–102.
[116] 张艳, 韩润生, 魏平堂, 邱文龙.2015. 云南昭通铅锌矿pH-logfO2和pH-loga相图对铅锌共生分异的制约. 中国地质 , 42 (2) :607–620.
[117] 章永梅, 顾雪祥, 彭义伟, 郑硌, 张岩, 高海军, 董树义.2014. 安哥拉Bonga碳酸岩型铌矿床地质地球化学特征及成因. 地学前缘 , 21 (5) :50–68.
[118] 赵劲松, 邱学林, 赵斌, 涂湘林, 虞珏, 芦铁山.2007. 大冶-武山矿化夕卡岩的稀土元素地球化学研究. 地球化学 , 36 (4) :400–412.
[119]
[120] 赵一鸣, 张轶男, 林文蔚.1997. 我国夕卡岩矿床中的辉石和似辉石特征及其与金属矿化的关系. 矿床地质 , 16 (4) :318–329.
[121] 赵一鸣, 李大新, 董永观.2002. 再论含Pb-Zn(Ag)锰质夕卡岩建造. 矿床地质 , 21 (S1) :548–552.
[122] 赵一鸣, 李大新.2003. 中国夕卡岩矿床中的角闪石. 矿床地质 , 22 (4) :345–359.
[123] 赵振华. 1997. 微量元素地球化学原理. 北京: 科学出版社 .
[124] 郑永飞, 徐宝龙, 周根陶.2000. 矿物稳定同位素地球化学研究. 地学前缘 , 7 (2) :299–320.
[125] 朱飞霖, 陶琰, 胡瑞忠, 廖名扬, 王银喜, 李玉帮.2011. 云南镇康芦子园铅-锌矿的成矿年龄. 矿物岩石地球化学通报 , 30 (1) :73–79.