岩石学报  2012, Vol. 28 Issue (12): 3798-3808   PDF    
湖南黄沙坪W-Mo-Bi-Pb-Zn多金属矿床硫铅同位素地球化学研究
祝新友1,2, 王京彬1,2,3, 王艳丽3, 程细音1, 傅其斌4     
1. 中色地科矿产勘查有限公司,北京 100012;
2. 北京矿产地质研究院,北京 100012;
3. 有色金属矿产地质调查中心,北京 100012;
4. 昆明理工大学,昆明 650093
摘要: 湖南黄沙坪W-Mo-Bi-Pb-Zn多金属矿床规模大、矿种多、范围小、分带明显,是南岭有色金属成矿带的代表性矿床之一。成矿地质体为碱长花岗斑岩,与下石炭统灰岩接触带发生大规模矽卡岩化,形成大型钨、钼、铋、萤石以及铁的共生矿床。围绕矽卡岩向外,分布铜锌、铅锌、铅锌银的分带,对应的矿化组合分别为粗粒磁黄铁矿-闪锌矿-黄铜矿、中粗粒磁黄铁矿-闪锌矿-方铅矿、胶状黄铁矿-闪锌矿-方铅矿。围绕花岗斑岩,硫化物矿物的δ34S值呈带状分布,其δ34S总体变化为2.3‰~17.5‰,花岗斑岩中浸染状辉钼矿δ34S为17.1‰,矽卡岩中硫化物δ34S>15‰,矽卡岩附近及外侧的铅锌矿体10‰ < δ34S < 15‰,外围的铅锌银矿体δ34S < 10‰。下石炭统中代表沉积特点的细粒浸染状黄铁矿δ34S为-3.1‰~-22.6‰。铅同位素206Pb/204Pb为18.525~18.603,207Pb/204Pb为15.706~15.792,208Pb/204Pb为38.889~39.178。综合研究表明,黄沙坪矿床成矿物质硫和铅主要来自花岗斑岩。经硫同位素热力学平衡计算,引起δ34S值围绕花岗斑岩体分带的主要原因是随温度下降以及物理化学条件变化导致的热力学分馏作用,其次是沉积岩围岩中低δ34S值硫的加入。对南岭地区花岗岩、古生代地层等的δ34S值对比研究发现,引起花岗斑岩岩浆高δ34S值的主要原因是深部富含硫化物(δ34S值高)地层对富含挥发份(Li-F)的碱长花岗岩岩浆的混染作用,其次是成矿作用过程中地层与岩浆的相互作用(包括同化混染)。围绕黄沙坪矿床,湘南地区矽卡岩型钨多金属矿存在一个较高的δ34S值分布区。宝山矽卡岩型Cu-Mo-Pb-Zn矿床矿石硫化物δ34S为-1‰~3.6‰,206Pb/204Pb为18.602~18.672,207Pb/204Pb为15.693~15.780,208Pb/204Pb为38.901~39.186。因此,宝山矿床与黄沙坪矿床的物质来源和成矿机制不同,宝山矿床硫、铅同位素组成集中,分布范围不同于黄沙坪,成矿物质来自岩浆岩。黄沙坪、宝山矿床代表了南岭地区燕山早期存在两类不同性质的岩浆活动与成矿组合。
关键词: 硫同位素     铅同位素     矿床成因     成矿机制     钨矿     铅锌矿     湖南省    
Sulfur and lead isotope constraints on ore formation of the Huangshaping W-Mo-Bi-Pb-Zn polymetallic ore deposit, Hunan Province, South China
ZHU XinYou1,2, WANG JingBin1,2,3, WANG YanLi3, CHENG XiYin1, FU QiBin4     
1. Sinotech Mineral Exploration Co., Ltd, Beijing 100012, China;
2. Beijing Institute of Geology and Mineral Resources, Beijing 100012, China;
3. China Non-ferrous Metals Resource Geological Survey, Beijing 100012, China;
4. Kunming University of Science and Technology, Kunming 650093, China
Abstract: Huangshaping W-Mo-Bi-Pb-Zn polymetallic ore deposit, located at Yizhang County, Hunan Province, South China, is one of the most famous deposits in Nanling tungsten belt. This deposit shows characteristics of multiple elements mineralization, large ore reserve in a small area < 1km2, and typical zonation from skarn-type mineralization to distal hydrothermal veins. The skarn mineralization of W, Mo, Bi and fluorite occurred at the contact zone between alkali feldspar granite-porphyry and the Lower Carboniferous limestone. The Cu-Zn, Pb-Zn and Pb-Zn-Ag sulfide ore bodies distribute from the skarn zone to distal limestone strata, and the mineral assemblage are coarse-grained pyrrhotite (Po) -sphalerite (Sp) -chalcopyrite (Cp), Po-Sp-Galena (Ga), and colloidal pyrite-Sp-Ga, respectively. The δ34S values of the sulfide minerals are from 2.3‰ to 17.5‰, and a zonation of the δ34S value is observed around the granite-porphyry. The δ34S value of molybdenite disseminated at granite-porphyry is 17.1‰, the sulfide minerals from the skarn zone show δ34S values >15‰, the lead-zinc orebodies near the skarn show δ34S values between 10‰ and 15‰, and the distal Pb-Zn-Ag ore bodies form the limestone strata show δ34S values < 10‰. The disseminated pyrite in the Lower Carboniferous sediment rocks shows δ34S values from -3.1‰ to -22.6‰. The lead isotopic compositions of the sulfide minerals are as following: 206Pb/204Pb of 18.525~18.603, 207Pb/204Pb of 15.706~15.792, and 208Pb/204Pb of 38.889~39.178. According to these results, it is suggested that the sulfur and lead of the Huangshaping deposit may have originally came from the granite-porphyry, but also with minor contributions from the limestone strata. The δ34S zonation around the porphyry is probably caused by the thermodynamic isotope fractionation related to the temperature decreasing, varying of the redox conditions and pH values. Alternatively, it is also possible that some amounts of sulfur with lower δ34S value may have also been derived from the Lower Carboniferous sedimentary rocks. The high δ34S values in granite porphyry magma may have been caused by incorporation of the sulfides with high δ34S values from the sedimentary rocks at depth (i.e. Devonian carbonate rocks) and by contamination of the carbonate wall rocks during ore-forming processes. The Baoshan Cu-Mo-Pb-Zn skarn deposit is located about 9km north of the Huangshaping deposit, the sulfide minerals from Baoshan show δ34S values of -1‰ to 3.6‰, and 206Pb/204Pb of 18.602~18.672, 207Pb/204Pb of 15.693~15.780, and 208Pb/204Pb of 38.901~39.186. The Baoshan Cu-Mo-Pb-Zn polymetallic ore deposit show different origin and mineralization mechanism from the Huangshaping deposit. The Bashan deposit has a very narrow range of sulfur and lead isotopic compositions, and the hydrothermal fluid should have derived mainly from the granodiorite-porphyry magma. Therefore, these two deposits may represent the two types of magmatic activity and mineralization at Late Jurassic period in Nanling region, South China.
Key words: Sulfur isotope     Lead isotope     Ore genesis     Tungsten deposit     Lead-zinc deposit     Skarn deposit     Hunan Province    

南岭是世界著名的钨锡多金属矿产地,其成矿作用多与中细粒碱长花岗岩成因演化有关,成矿物质主要来自于岩浆岩的分异演化作用(图 1)。黄沙坪是南岭有色金属成矿带的代表性矿床之一(图 2),花岗斑岩为成矿地质体,未出露地表,成矿部位水平横截面积也仅有200m×300m,探明铅锌金属资源量超过200万吨,近年来探明的矽卡岩型钨、钼、铋、萤石矿均达到大型以上规模(共生)。不仅成矿地质体(花岗斑岩)与其它矽卡岩型矿床不同,铅锌矿的比例也非常大,规模远大于柿竹园、瑶岗仙等钨多金属矿。导致这些特点的主要机制是什么?本文通过对黄沙坪矿床花岗斑岩体及矽卡岩成矿系统的硫铅同位素组成研究,并与宝山矿床和南岭地区其它钨锡多金属矿床开展对比, 提出黄沙坪矿床花岗斑岩独特的成因与演化机制可能是导致黄沙坪矿床成矿特殊性的主要原因。

图 1 南岭钨锡多金属成矿带地质略图 图中标出了矽卡岩型钨锡矿δ34S平均值 Fig. 1 Sketch geological map of Nanling region δ34S value of skarn type W-Sn deposits in this figure

图 2 黄沙坪矿区地质简图(a)及16线剖面简图(b) Fig. 2 Geological map of surface (a) and section map (b) of the Huangshaping W-Mo-Pb-Zn polymetallic ore deposit
1 黄沙坪矿床地质概况

黄沙坪矿床地处华厦地块中部,骑田岭岩体的西北侧。矿区发育一套晚古生代浅海相台地沉积岩,包括上泥盆统锡矿山组(D3x)砂岩,下石炭统岩关阶石磴子组(C1s)厚层灰岩、测水组(C1c)细砂岩及泥质砂岩等,呈NNE向展布。由C1s灰岩与C1c细碎屑岩构成向东倾伏的上银山倒转背斜,沿背斜核部附近发育NNE纵向逆冲断裂F1、F3、F2,以及NWW向横向断裂F0,共同构成黄沙坪矿床“田”字型构造框架(图 1)。石英斑岩沿F1、F2侵入,形成51#和52#岩株(0.3km2)和NNE向岩脉。英安斑岩(54#)在矿区南部沿F0侵入。花斑岩(304#)和花岗斑岩(301#)分别位于矿区西北部和东南部,均为隐伏小岩株,岩体与C1s灰岩接触带出现大规模的矽卡岩。地质地球化学研究显示,花斑岩与花岗斑岩属同一类岩体,只是侵位高度和蚀变程度不同,二者沿NNW向断裂侵入,NNW与NNE向构造一同控制区内成矿作用(祝新友等,2010)。

黄沙坪矿床矿化类型包括钨钼多金属矿和铅锌矿,前者为矽卡岩型,分布于花岗斑岩、花斑岩与灰岩的接触带,后者呈囊状、不规则脉状分布于矽卡岩边部和外围,主要分散于花岗斑岩的北西部灰岩中,二者共同构成黄沙坪矽卡岩型矿床成矿系统。不同矿化类型围绕花岗斑岩构成良好的分带,矽卡岩型钨钼矿分布于花岗斑岩岩体与灰岩的接触带,主要矿物为石榴子石、透辉石、萤石、磁铁矿、白钨矿、辉钼矿、辉铋矿、黄铁矿等,发育碱交代脉(钾长石-萤石脉)。钨钼铋矿化主要形成于退变质-氧化物阶段,与碱交代脉两侧阳起石化、磁铁矿化相伴随。铅锌硫化物属矽卡岩演化的晚期阶段。分布于矽卡岩边部的硫化物矿体以囊状块状磁黄铁矿为主,其中含粗粒黑色闪锌矿(图 3a),方铅矿较少,局部含黄铜矿;外侧的铅锌矿呈不规则脉状分布于灰岩中,或石英斑岩边部,块状,含大量磁黄铁矿、闪锌矿和方铅矿,是铅锌矿的主要类型,在矿床深部这类脉体中富含黄铜矿,局部达到工业品位;最外侧的是铅锌银矿,分布于远离花岗斑岩的灰岩中,不规则脉状,矿石呈块状,矿物以方铅矿、闪锌矿、黄铁矿为主,主要呈胶状结构(图 3b),富含银(王静纯等,1993),局部见少量镜铁矿。

图 3 黄沙坪矿床主要铅锌矿石结构 (a)-矽卡岩边部的铅锌矿体,中粗粒状,闪锌矿中含乳滴状黝铜矿和黄铜矿;(b)-远离矽卡岩产于灰岩中的铅锌矿体,胶状结构,胶状白铁矿交代黄铁矿.Apy-毒砂; Cal-方解石; Cp-黄铜矿; Sp-闪锌矿; Tt-黝铜矿; Mc-白铁矿; Py-黄铁矿 Fig. 3 Photomicrographs of lead-zinc ore (a)-lead-zinc ore in skarn, texture of coarse-grain; (b)-lead-zinc ore in distal limestone rocks, pyrite alterated by colloidal marcacite. Apy-arsenopyrite; Cal-calcite; Cp-chalcopyrite; Sp-sphalerite; Tt-tetrahedrite; Mc-marcasite; Py-pyrite
2 硫、铅同位素分析结果 2.1 硫同位素分析结果

前人在黄沙坪矿床开展了大量的硫同位素分析研究工作,不同作者共测定硫化物硫同位素样品146件(童潜明,1986车勤建,2005息朝庄等,2009),其中,黄铁矿43件,闪锌矿50件,方铅矿44件、黄铜矿4件、磁黄铁矿3件、毒砂2件。除1件黄铁矿δ34S为-2.2‰外,其它所有数据均落在2.3‰~17.5‰范围内,大部分数据分布于6‰~17‰范围内(图 4)。黄铁矿2.8‰~16.8‰,平均10.5‰,闪锌矿4‰~17.5‰,平均12.3‰,方铅矿2.3‰~14.5‰,平均9.8‰、黄铜矿3.8‰~5.9‰,平均4.5‰,磁黄铁矿4.2‰~10.6‰,平均7.0‰、毒砂13‰~14.9‰,平均14.0‰。

图 4 黄沙坪和宝山矿床硫同位素组成直方图 (a)-黄沙坪,息朝庄(2009);(b)-黄沙坪,本文;(c)-宝山,本文 Fig. 4 Histogram of the δ34S values of sulfide minerals from the Huangshaping and Baoshan deposits

这些分析结果显示出黄沙坪具有与南岭地区大多数钨矿明显不同的硫同位素组成。早在20世纪80年代,童潜明(1986)对黄沙坪273中段作了系统的硫化物硫同位素分析工作,发现黄铁矿、闪锌矿、方铅矿的δ34S值均以花岗斑岩为中心向北、向西降低。该分带与硫化物矿物中微量元素含量的分带相对应,进而确定花岗斑岩为铅锌矿成矿母岩。息朝庄等(2009)总结出301矿段较304矿段δ34S值高,硫源可能是岩浆来源与海水(地层)硫混合作用形成。

为进一步探讨黄沙坪矿床成矿物质来源与成矿机制,在前人工作基础上,我们针对黄沙坪矿床中不同地质体、不同产状作了系统的采样分析。分析单位为核工业北京地质研究院分析测试研究中心,分析方法参照行标“DZ/T 0184.14-1997硫化物中硫同位素组成的测定”(中华人民共和国地质矿产行业标准, 1997)。硫同位素分布范围与前人分析结果一致(表 1),δ34S为3.0~17.4‰(除一件黄铁矿为24.2‰外)。其中,黄铁矿8.7‰~24.2‰,平均14.6‰;闪锌矿3.8‰~17.4‰,平均10.5‰;方铅矿3.0%~15.0‰,平均9.0‰;辉钼矿14.0‰~17.1‰,平均14.8‰。δ34S值呈现出辉钼矿>黄铁矿>闪锌矿>方铅矿的特征,表明这些矿物间硫同位素基本上达到平衡分馏状态(图 4)。

表 1 黄沙坪钨钼铅锌多金属矿和宝山铜钼铅锌矿床硫铅同位素分析结果 Table 1 The sulfur and lead isotope compositions of sulfide minerals from the Huangshaping W-Mo-Pb-Zn and Baoshan Cu-Mo-Pb-Zn ore deposits

①黄崇轲等.1995.南岭地区与钨锡(铅锌)相关银矿成矿规律与成矿预测研究报告.白银地质基金委员会研究报告(内部资料)

宝山矿区的硫同位素组成较集中,全部样品δ34S为-1‰~3.6‰,与车勤建(2005)测定的δ34S (4.4‰~4.5‰)范围略有不同。其中,黄铁矿-1.0‰~2.7‰,平均1.4‰;闪锌矿1.8‰~3.6‰,平均2.4‰;方铅矿-0.3%~2.1‰,平均0.94‰。呈塔式分布,其中黄铁矿主要取自矽卡岩,闪锌矿和方铅矿主要取自铅锌矿石,可能是黄铁矿与铅锌矿物的硫同位素显示不平衡状态的主要原因。

硫化物的硫同位素值与样品产出的空间位置密切相关,图 5是以黄沙坪矿床20m中段为研究对象,将20m中段及上、下-16m、56m中段的样品投影在20m中段后的结果,显然,硫化物的硫同位素δ34S值围绕301岩体呈现明显的分带现象。以301花岗斑岩体为中心,产于矽卡岩中的硫化物,包括辉钼矿、部分黄铁矿和磁黄铁矿,大部分样品的δ34S>15‰,外侧的铅锌矿体10‰ < δ34S < 15‰,更外侧,主要是铅锌银矿体,也包括部分334#花斑岩体有关矽卡岩中的铅锌矿,δ34S < 10‰, 甚至更低。这种围绕花岗斑岩体形成的分带,其总体走向呈NNW向,SSE方向开口,与矿化分带完全一致,花岗斑岩边部为矽卡岩,外侧为不规则脉状铅锌矿,更外侧为脉状铅锌银矿。

图 5 黄沙坪矿床20m中段硫化物硫同位素分带 Fig. 5 The zoning of the δ34S values of sulfide minerals at 20m level of the Huangshaping deposit

这种分带现象在垂向上同样明显,上部中段较下部中段具有明显不同的δ34S值。上部中段δ34S值高,往下逐渐减小。

依据20m中段硫同位素δ34S值的变化规律,NNW方向上,自矽卡岩往外,δ34S值的递减梯度大约为3‰/100m,这与成矿温度的递减梯度约60~80℃/100m相对应(矽卡岩温度约550℃,外带铅锌银矿成矿温度约150~200℃)。

2.2 铅同位素分析结果

前人也开展了大量的铅同位素研究工作,其中,息朝庄等(2009)测定的黄沙坪矿床共47件硫化物样品,测试矿物以方铅矿为主,这些资料表现在铅同位素图件上部分数据点非常分散,一部分甚至分布于地球增长线右侧。数据是由桂林矿产地质研究院分析完成,由于年代较早,且设备老化,测定结果误差较大。在本项目的研究中,共测定了29件铅同位素样品,其中黄沙坪16件,宝山13件(表 1)。分析单位为核工业北京地质研究院分析测试研究中心,分析方法参照国标GB/T 17672-1999《岩石中铅锶钕同位素测定方法》(中华人民共和国国家标准, 1999)。

本文测定的黄沙坪矿床铅同位素组成与息朝庄(2009)的数据平均值接近。这些数据共分3批测试,第二批测试结果误差相对较大,尤其是其中含Pb较低的辉钼矿,这3件辉钼矿样品的208Pb/204Pb明显高,相对富含放射性成因铅。如果不考虑这3件样品,黄沙坪矿床与宝山矿床的样品大体呈平行的线性分布关系,黄沙坪矿床分布于左侧,这种线性关系与204Pb的分析误差有关。在后面的分析讨论中,我们也扣除了这3件样品。

对比黄沙坪、宝山、瑶岗仙三个矿床,黄沙坪与瑶岗仙矿床铅同位素组成更接近,宝山矿床具有明显多的放射性成因铅。

3 讨论 3.1 硫的来源与演化

有关黄沙坪矿床硫同位素的研究较多,不同研究者获得的δ34S数据值范围大体一致, 但对硫的来源有不同认识。部分研究者认为成矿物质硫来自于花岗斑岩(童潜明,1986李裕祖, 1986),也有一些研究者认为有部分地层硫的加入(毕承思,1986),或地层与岩浆硫的混合(息朝庄等,2009)。

黄沙坪矿区矽卡岩成矿系统自矽卡岩中心至外围铅锌银矿床,硫化物δ34S呈现有规律的变化,形成以碱长花岗斑岩为中心的环带,因此,本文认为成矿物质硫主要来自于花岗斑岩体。由于矽卡岩外围铅锌矿成矿温度较低,其热力学分馏可引起δ34S值降低,而下石炭统围岩中低δ34S值硫的加入也可能导致外围铅锌矿体δ34S值降低。那么,黄沙坪矿床硫化物δ34S值围绕矽卡岩(或花岗斑岩)降低的主要机制是什么呢?

1)矽卡岩成矿系统自内向外降温引起的热力学分馏是引起δ34S值分带的重要原因

由于黄沙坪矿床中没有明显的有机质或生物参与成矿的证据,因此,硫同位素的分馏机制可能主要是通过热力学分馏完成的。反映成矿物理化学环境变化的矿物组合在退变质矽卡岩-氧化物阶段为磁铁矿-钾长石;硫化物阶段,内带铅锌矿为磁黄铁矿-黄铁矿-方解石,外带铅锌银矿为黄铁矿-方解石-绢云母-少量镜铁矿。显示出成矿作用自早至晚,成矿系统自内而外,成矿温度降低,成矿环境自偏还原态向偏氧化态、偏碱性向偏酸性发展。若来自岩浆的流体δ34SΣS值未发生显著变化,根据大本硫同位素分馏模型(Ohmoto and Lasaga, 1982),晚期以及外围沉淀的硫化物矿物δ34S值必然逐渐降低。因此,成矿系统自内向外环境变化引起的热力学分馏会导致δ34S值的逐渐降低。

我们大体将硫化物根据产状分为:花岗斑岩中浸染状硫化物、矽卡岩退变质阶段的硫化物、矽卡岩及近矽卡岩中产出的不规则状硫铁铅锌矿、远离矽卡岩的脉状铅锌矿。流体包裹体研究显示,矽卡岩中石榴子石均一温度>500℃,退变质氧化物阶段的温度为400℃±,矽卡岩中的辉钼矿、黄铁矿多形成于此阶段。硫化物阶段形成块状硫铁铅锌矿石,其中,矽卡岩内及附近富含磁黄铁矿,形成温度约为300~350℃,远离矽卡岩磁黄铁矿含量降低,逐渐变为铅锌矿,分布于外侧的铅锌银矿体温度更低,一般在200~250℃,矿石中常见胶状结构。据此,根据不同矿物间的硫同位素分馏关系,以及不同矿物产出关系,以辉钼矿、闪锌矿为例,计算该矿物形成时对应的体系δ34SH2S值(表 2)。计算结果显示,岩体及附近或高温阶段流体中δ34SH2S为13.2‰~16.2‰,岩体边部为16.6‰,外侧为9.8‰,远离矽卡岩产于灰岩中的脉状矿体,其成矿流体δ34SH2S为4.3‰。

表 2 根据辉钼矿、闪锌矿计算的体系H2S的δ34S值 Table 2 The calculated δ34S value of the H2S in ore-forming fluids based on the δ34S values of molybdenite and sphalerite

假定成矿过程中形成矽卡岩及铅锌矿的流体均直接来自碱长花岗岩岩浆,且热液δ34SΣS值未发生明显变化,可以根据大本模型(Ohomoto and Rye, 1979)计算出铅锌银矿形成时体系的x=SO42-/H2S摩尔比为0.52。如果考虑外围铅锌矿、铅锌银矿的成矿热液可能是内带矿物沉淀后的热液向外迁移的产物,形成外带铅锌银矿的流体δ34SΣS值可能更高,要达到δ34SH2S=4.3‰的水平,x=SO42-/H2S摩尔比会更高。虽然铅锌银矿体中出现微量镜铁矿,但主要矿物为闪锌矿、方铅矿、黄铁矿,这种矿物组合不大可能形成于氧化环境之中。

另外,宝山矿床的矽卡岩中也发育镜铁矿,发育程度甚至超过黄沙坪,但矿床δ34S分布于0值附近,无显著分馏。因此,黄沙坪晚阶段尤其是外围矿体中硫化物δ34S的降低(平均δ34S=4.7‰)不大可能全部因热力学分馏导致,应该还存在其它原因,这就是外来低δ34S硫的参与。

2)下石炭统围岩中低δ34S硫的参与也是外围铅锌矿体中δ34S值降低的重要原因

矽卡岩阶段→硫化物阶段硫同位素经历了复杂的演化,分布于矽卡岩中的硫化物δ34S值远高于外围的铅锌矿体。流体包裹体研究发现,位于矽卡岩阶段与硫化物阶段之间的退变质-氧化物阶段(钨钼成矿阶段),成矿流体性质发生迅速变化,不仅有氧化还原环境的变化,也存在大量大气降水的参与。伴随大气降水的外来物质尤其是来自围岩中的硫的加入也可能导致体系硫同位素组成发生变化。

黄沙坪容矿地层石磴子组(C1s)灰岩夹炭质薄层灰岩,局部含有浸染状细粒黄铁矿,前人测定了其中的黄铁矿δ34S为-3.1‰,而矿区测水组(C1c)炭质页岩中的黄铁矿结核两个样品δ34S值分别为-11.7‰和-22.6‰(王育民等, 1988, 转引自桂林冶金地质研究院),显示矿区容矿围岩δ34S值低,远低于花岗斑岩的δ34S值。这种来自围岩中的低δ34S值的硫的加入将引起黄沙坪矿床花岗斑岩外侧铅锌矿δ34S值的降低。Ishihara et al.(2003)在研究南岭钨矿时也认识到这一点,认为黄沙坪硫同位素分带与大量地层中低δ34S值硫加入有关。

位于湖南宜章县的界牌岭锡钨铍多金属矿亦产于石磴子组(C1s)灰岩顶部附近,产出层位与黄沙坪矿床相当。矿床中发育大范围的顺层萤石化,其中黄铁矿、闪锌矿、方铅矿等硫化物δ34S为-1.0‰~7.5‰,平均4.48‰,塔式分布(彭振安,1991)。矿区未出露成矿岩体,流体包裹体均一温度200~250℃,属远离岩体的矿化。虽然如此,硫同位素组成特点显示硫主要源自岩浆,围岩的影响有限。这也显示出黄沙坪矿床外围铅锌矿体的硫也主要来自于碱长花岗斑岩,围岩只提供了少量的低δ34S硫。也就是说,造成黄沙坪矿床硫同位素分带的主要原因可能是自内而外的物理化学条件和系统变化导致的热力学分馏引起,其次是围岩中低δ34S硫的加入。

综上所述,形成花岗斑岩的岩浆具有较高的δ34S值,围绕花岗斑岩和矽卡岩的硫同位素分带现象主要原因是由于体系成矿温度和氧化还原环境变化引发的热力学分馏引起,但外界尤其是来自容矿地层中低δ34S值的硫的参与也是重要原因。因此,黄沙坪矿床成矿硫主要来自于花岗斑岩体,少量来自矿区下石炭统。

宝山矿床硫化物矿物δ34S为-1‰~3.6‰,平均1.5‰,接近陨石硫的组成,硫主要源自于花岗闪长斑岩。

3.2 花岗斑岩岩浆高δ34S值的成因问题

锡钨多金属矿成矿作用与酸性侵入岩的岩浆分异与岩浆热液作用有关,一般认为,这些矿床成矿物质来自岩浆岩。国内外绝大多数钨矿床,尤其是南岭地区大部分的钨矿区,矿石硫同位素大多分布于0值附近,尤其是石英脉型钨矿,显示出岩浆硫的性质。如湖南瑶岗仙石英脉型钨矿床δ34S为-2.6‰~3.9‰(n=42,样品数,下同),平均-0.4‰(陈依壤,1988;祝新友等,2012);川口钨矿δ34S为0.6‰~2.2‰(n=7)(车勤建,2005);江西淘锡坑黄铁矿的δ34S为0.1‰~-2.1‰(n=20),平均-1.2‰(宋生琼等,2011);大吉山钨矿δ34S为-0.2‰~-3.7‰(n=24),平均-2..4‰(庄龙池等,1991);西华山钨矿δ34S为-1.1‰~2.2‰(n=33)(刘家齐,1989);黄砂脉型钨矿床δ34S为-1.3‰~2.2‰(n=59)(冯志文等,1989);广西栗木δ34S为-3.5‰~2.5‰(n=28)(徐文炘,1988)等。部分矽卡岩型钨锡矿δ34S值呈现不同程度的正值,如锡田矿区δ34S为3.2‰~3.7‰(n=3)(车勤建,2005)、杮竹园野鸡尾矿区δ34S为2.8‰~8.9‰(n=21)(刘悟辉等,2006)、香花岭矿区δ34S为-1.0‰~6.7‰(周涛等,2008)。湘南芙蓉矿区δ34S则分为两组,一组在零值附近并呈现不同程度的正值(-0.5‰~8.5‰, 刘悟辉等,2006; 蒋少涌等, 2006);另一组δ34S为较低的负值(-20.4‰~-10.0‰,蒋少涌等,2006), 指示硫可能来自围岩。

国外一些著名矽卡岩型钨矿床也有类似特点,如加拿大MacDame矽卡岩型钨矿平均δ34S为7.7‰(Cooke et al., 1984);韩国Sannae W-Mo矽卡岩型钨钼矿床中辉钼矿δ34S为5.2‰~6.0‰, 黄铁矿δ34S为6.0‰~7.1‰(Shelton et al., 1986)。Weolag矽卡岩型钨钼矿中,辉钼矿δ34S为3.0‰~5.7‰、黄铁矿δ34S为4.5‰~7.6‰、磁黄铁矿δ34S为4.8‰~5.1‰、黄铜矿δ34S为4.8‰、闪锌矿δ34S为3.0‰~5.7‰、方铅矿δ34S为3.9‰,水岩值也很高(So et al., 1983)。秘鲁San Cristobal矽卡岩型钨矿贱金属含量较高,δ34S较高,具有蒸发岩来源,而3个石英脉型钨矿Pasto Bueno, Panasqueira, 和Tungsten Queen矿床,δ34S在0值附近,具有岩浆源(Campbell, 1987)。石英脉型钨矿硫同位素更多显示出岩浆来源的特点,如美国北卡罗来那Queen钨矿床,δ34S在0值附近,具有深源特点(Casadevall and Rey, 1980),南非Namaqualand石英脉型钨矿,δ34S为3.6‰~4.5‰(Raith et al., 2000)。

黄沙坪钨多金属矿硫化物δ34S为3‰~17‰,明显不同于南岭地区大多数锡钨多金属矿δ34S接近0值幔源硫或岩浆岩硫的特点,以往的研究均将其简单地解释为地层来源或复杂来源(息朝庄等,2009),或回避这一问题。已有研究证明,碱长花岗斑岩体是矽卡岩及铅锌矿床的成矿母岩,为成矿地质体(童潜明,1986祝新友等,2010)。硫化物硫同位素δ34S值围绕花岗斑岩体发生分带说明,岩浆是成矿作用中硫的重要供体。由于高温下辉钼矿、黄铁矿与流体H2S的分馏很小,且高温还原条件下流体中S的主要形式为H2S (Ohmoto and Lasaga, 1982),因此,取自301花岗斑岩体附近的矽卡岩中辉钼矿、黄铁矿的δ34S值可大体代表岩浆的值,花岗斑岩岩浆的δ34S>15‰。取自花岗斑岩中网脉状矿化石英细脉中的辉钼矿δ34S为17.1‰,也说明形成花岗斑岩的岩浆具有较高的δ34S值。

本文的研究赞同成矿物质硫主要来自于花岗斑岩,这也意味着花岗斑岩岩浆具有较高的δ34S值,导致这一现象的原因是什么?有以下3种可能:

1)深部高δ34S值硫化物及膏盐层的混合

华南地区泥盆系-石炭系为连续沉积,上泥盆统中富含浸染状硫化物,其δ34S较高。在广东凡口铅锌矿、大宝山多金属矿区的中、上泥盆统碳酸盐岩中,已获得大量的浸染状细粒黄铁矿δ34S值,主要分布于10‰~25‰(张术根等,2009)。黄沙坪矿区含矿地层为下石炭统,深部存在泥盆系。形成钨锡多金属矿的花岗斑岩富含挥发份,尤其是其中Li-F含量较高,有利于萃取围岩中的硫,从而改变岩体的性质。正是由于岩体富硫的特性,形成规模巨大的铅锌矿床。相对于杮竹园的钨钼规模,黄沙坪矿床的铅锌矿规模大得多。

2)碳酸盐岩中硫酸盐的影响

一般情况下,海相沉积形成的碳酸盐岩中含有微量的硫酸根离子,其δ34S值较高,相当于海水硫酸盐的水平(Claypool et al., 1980;Strauss,1983)。湘南地区大多数矽卡岩型矿床的赋矿围岩为中、上泥盆系碳酸盐岩,富硫且δ34S值较高。矽卡岩型钨矿床δ34S值高于石英脉型钨矿,可能与灰岩中微量硫酸根离子的加入有关,同时,也与矽卡岩型钨矿床成矿机制有关,成矿过程中碳酸盐岩地层对岩浆产生了较大程度的同化混杂作用。研究发现,在石榴子石、透辉石矽卡岩形成之前,碳酸盐岩围岩与岩浆的混染作用已经进行,并改变了岩浆的成分,如岩体边部、岩枝或岩脉的钠/斜长石的An值升高(另文论述)。这一同化混染过程也会提高岩浆的硫的含量、导致岩浆δ34S值的升高。

3)岩浆演化的液态分离作用

纵观南岭地区两类钨矿床,脉型矿床δ34S值范围狭窄,且绝大多数局限在0值附近。矽卡岩型矿床δ34S值普遍较高。研究已经发现,脉型矿床形成于相对封闭系统中(另文发表),矽卡岩型矿床形成于相对开放系统中。开放系统的岩浆热液分异过程有大量气态物质的逸失,可能会导致残留的岩浆中相对富集34S,进而有利于矽卡岩型矿体中硫化物δ34S值的升高。有关这一点还需要寻找进一步的证据。

综合黄沙坪地区岩浆演化与地球化学特征,花岗斑岩岩浆侵入时受到高δ34S值的围岩混染应是岩浆δ34S值增高的主要原因。有关岩浆上侵过程中遭受围岩中硫混染的研究,多集中于超基性岩有关Cu-Ni-Pt矿床方面(Ripley et al., 2003王瑞廷等,2005),在一些矽卡岩型铜、铁矿床以及斑岩铜矿床的研究中,很多成果也发现围岩的混染作用,尤其是硫的加入对金属物质富集成矿的影响。在湘南地区,黄沙坪矿床矽卡岩的硫化物δ34S值最高,外围的柿竹园、新田岭、锡田等矽卡岩矿床硫化物δ34S值也较高,大体以黄沙坪为中心,周边的矽卡岩型矿床构成δ34S值的高值区(图 1)。瑶岗仙和尚滩矽卡岩矿床δ34S为-1.3‰~3.6‰,平均1.0‰(作者未发表数据),大体可作为其边界;江西崇义宝山矽卡岩型钨矿δ34S平均1.7‰(朱焱龄等,1981)。这一地区性的δ34S值增高与矽卡岩矿床分布相对应,显示这一区域的围岩对硫的混染作用很可能对矽卡岩演化及其成矿作用有着制约性影响,这种影响可能首先从改变岩浆岩的性质开始。

3.3 铅同位素组成

图 6可知,黄沙坪、宝山以及瑶岗仙、凡口矿床矿石的μ值(238U/204Pb)相近,大体在9.7左右(图 6)。凡口铅锌矿床206Pb/204Pb值比其他几个矿体稍低,这可能与该矿床的形成与盆地卤水作用有关(祝新友等,2013),容矿地层为中、上泥盆统碳酸盐岩,成矿物质来源与围岩相对应。宝山、黄沙坪、瑶岗仙矿床的206Pb/204Pb值变化基本相同,其成矿作用与燕山早期中酸性侵入岩有关,铅主要来自于岩浆岩。肖红全等(2003)对柿竹园矿田金船塘矿床的铅同位素分析结果也表明该矿床的铅来源于千里山花岗岩浆,他们还获得一条Pb-Pb等时线年龄,与千里山岩浆侵入年龄一致,进一步说明成岩成矿的一致性。

图 6 黄沙坪矿床207Pb /204Pb-206Pb /204Pb演化图解 虚线为基于Doe and Zartman (1979)的上地壳、造山带和地幔演化图解(图中也标出了宝山、瑶岗仙和凡口矿床的铅同位素数据,供对比,数据引自祝新友等, 2012) Fig. 6 207Pb/204Pb vs. 206Pb/204Pb diagram of the Huangshaping deposit The dash lines present the evolution of the upper crust, orogenic and mantle evolution curves adopted from Doe and Zartman (1979)

①祝新友,王艳丽,汤静如. 2012.湖南瑶岗仙钨矿典型矿床研究(湘南-粤北地区锡钨多金属矿床成矿规律总结研究专题).研究报告(内部资料)

Doe and Zartman (1979)的铅同位素来源演化图解上(图 6虚线),不同矿床具有相似的初始来源,样品点均落在上地壳演化曲线附近,这与岩浆岩的成因相对应,黄沙坪、瑶岗仙矿床的成矿花岗岩均属幔源物质上涌导致的地壳物质重熔的产物,宝山矿床花岗闪长岩被认为属I型花岗岩,有更多幔源物质的参与(路远发等,2006)。

黄沙坪矿床侵入岩成岩年龄为152±3Ma (SHRIMP U-Pb,雷泽恒等,2010),成矿年龄为158.4±1.3Ma (Re-Os等时线,齐钒宇等,2012)。宝山矿床花岗闪长斑岩成岩年龄SHRIMP U-Pb 158±2Ma,辉钼矿Re-Os等时线年龄为160±2Ma (路远发等,2006)。两个矿床侵入岩成岩年龄与成矿年龄大体相当,显然铅同位素的差异性主要不是年龄问题,而是来源与演化过程不同所致。

宝山与黄沙坪、瑶岗仙的差异更大程度上反映出成岩物质U/Th含量的差异,最终导致放射性成因铅含量的不同,尤其是206Pb/204Pb的不同。宝山矿床铅同位素单阶段模式年龄更加靠近岩体的年龄值,而钨矿有关花岗岩单阶段模式年龄远大于成岩成矿年龄,显示出宝山花岗闪长岩的铅同位素组成更加接近于地球铅同位素单阶段演化历史,这与宝山花岗闪长岩属I型花岗岩、有较多幔源物质参与的结论是一致的;而黄沙坪、瑶岗仙及其它一些钨矿区花岗岩成岩物质为古老变质岩的重熔,因此其铅同位素演化不符合单阶段演化模式的前提,从而其单阶段模式年龄值偏老。

3.4 硫化物的Pb-S同位素图解

以方铅矿、闪锌矿、黄铁矿的硫、铅同位素分析结果制作黄沙坪、宝山矿床的Pb-S同位素图解,如图 7。黄沙坪矿床铅同位素比值与δ34S值之间存在较显著的线性相关关系,δ34S值与206Pb/204Pb、207Pb/204Pb、208Pb/204Pb的相关系数分别为-0.720、-0.699、-0.724,为避免204Pb分析误差,图中也显示了δ34S值与206Pb/208Pb、207Pb/208Pb、207Pb/206Pb的相关关系,线性关系也很明显。这种线性关系一方面说明Pb来源与S的来源相同,另一方面也显示出2种来源物质混合的特性,即S和Pb均存在两种来源。

图 7 黄沙坪、宝山矿床矿石铅同位素值-δ34S图解 Fig. 7 δ34S value vs. lead isotope ratio of sulfide minerals in the Huangshaping and Baoshan deposits

结合黄沙坪矿床硫化物δ34S值围绕花岗斑岩分带的现象,主要物质应来自于花岗斑岩体,其具有较高的δ34S值,和相对低的放射性成因铅,另一类物质具有相对低的δ34S值,放射性成因铅含量相对高,对应的源区可能是矿区下石炭统容矿围岩,其中的浸染状黄铁矿内有较低的δ34S值。大多数情况下,沉积岩较同时代侵入岩往往具有更高含量的放射性成因铅。

宝山矿床δ34S、铅同位素值分布范围均很小,也不存在线性分布特点,显示S、Pb来源较单一,成矿物质来自于花岗闪长斑岩。

比较黄沙坪、宝山矿床的S、Pb同位素组成,在图 7中,二者分布特点、位置完全不同,也显示出两个矿床完全不同的岩石成因、成矿物质来源和矿床成因。黄沙坪矿床成矿岩体与南岭钨多金属矿区花岗岩相似,具有陆壳重熔花岗岩的特点,成矿矿种为W-Mo-Bi-Pb-Zn及萤石等,宝山矿床花岗闪长斑岩具有I型花岗岩特点,成矿矿种为Cu-Mo-Pb-Zn。

3.5 花岗斑岩岩浆演化及其对成矿的影响

在南岭钨矿床中,黄沙坪矿床成矿岩体小,矿床规模大,是南岭地区小岩体成大矿的典型。除此之外,与南岭钨多金属矿相比,成矿作用具有独特性:1)深部W-Mo-Bi规模巨大,铁品位较高,共生大型铁矿,而且外围发育大规模的铅锌矿化,其铅锌矿规模居南岭锡钨多金属矿床之首,超过杮竹园等矿床;2)成矿地质体为花岗斑岩,黄沙坪为南岭地区唯一与斑岩体有关的大型锡钨多金属矿床;3)岩体具有异常高的δ34S值。

这些特点也说明成矿地质体花岗斑岩具有独特的演化过程。花岗斑岩岩浆侵位深度更高,岩浆结晶温度和成矿温度相对更低,这就要求岩浆具有更高的Li-F等挥发份,使岩浆上侵至较高的高度(Robb, 2005)。在南岭地区,黄沙坪矿床的萤石含量不仅远高于石英脉型钨矿,也高于柿竹园矽卡岩型钨矿床。

岩体侵位高,赋存的围岩较新,黄沙坪地区位于下石炭统沉积岩之下更深部层位的泥盆系厚层碳酸盐岩富含浸染状黄铁矿,且具有较高的δ34S值。这种较大的迁移路径以及富F特点增加了岩浆与地层同化混染的机会,也可能是黄沙坪花岗斑岩岩浆具有高δ34S的重要原因。相比之下,杮竹园、瑶岗仙、锡田、香花岭与矿化有关的沉积岩为泥盆系,芙蓉、新田岭下部为大规模岩基,δ34S呈小的正值可能与此有关。

矽卡岩型矿床较石英脉型钨矿具有明显高的δ34S值,如果这些硫主要是深部及围岩地层同化混染的结果,将导致这类岩浆具有更高含量的硫,这也就可以进一步解释为什么富含硫的铅锌矿只出现于矽卡岩的外围,而石英脉型钨矿床则缺乏这类铅锌矿,虽然石英脉型钨矿对应的成矿花岗岩中Pb、Zn含量也较高。

宝山矿床的成矿特点不同于黄沙坪,而更相似于水口山、银山等矿床,属于南岭另一类矽卡岩成矿系列(胡永嘉等,1984)。由于岩浆性质不同,如F含量低,属正常I型花岗岩等,受围岩的同化混染程度低,成矿物质基本来自岩浆岩,受地层影响不明显。

4 结论

黄沙坪锡钨铅锌多金属矿床硫化物δ34S值围绕碱长花岗斑岩形成分带,花岗斑岩δ34S>14‰。铅同位素组成较均一,其比值与δ34S值呈线性分布,相关度较高。成矿的硫和铅主要来自花岗斑岩,少量的硫可能来自含矿地层,后者较前者更加富含放射性成因铅。引起δ34S值围绕花岗斑岩体分带的原因和机制主要是硫同位素的热力学分馏作用,以及成矿作用过程中来自地层的低δ34S值硫的加入。

黄沙坪矿床的花岗斑岩具有南岭地区钨矿床中最高的δ34S值,引起花岗斑岩岩浆高δ34S值的主要原因是深部富含硫化物(δ34S值高)的地层的混染作用,其次是成矿围岩碳酸盐岩的同化混染作用。富含挥发份(Li-F)的岩浆、赋矿地层较新且深部存在泥盆系等是黄沙坪花岗斑岩δ34S值增高的重要因素。矽卡岩型钨矿床成矿作用过程中地层与岩浆的强烈作用(包括同化混染)也会促进这一趋势,而石英脉型钨矿床这一特点不明显,δ34S值一般都在0值附近,成矿物质来自岩浆岩。

湘南地区围绕黄沙坪矿床存在一个δ34S值较高的区域矽卡岩型钨矿带。

宝山矿床硫、铅同位素组成更加集中,分布范围不同于黄沙坪,成矿物质来自岩浆岩。黄沙坪与宝山矿床均为矽卡岩型,但岩浆岩性质不同、成矿机制不同,二者代表了南岭地区燕山早期存在的两类不同性质的岩浆活动以及不同类型的成矿作用。

致谢 王静纯、李顺庭、唐攀科等参加了部分野外地质工作;同时,野外工作还得到黄沙坪铅锌矿江元成、王成发、汪林峰等人的大力支持和帮助;三位审稿人对本文提出了宝贵的修改意见;在此一并表示感谢。
参考文献
[] Bi CS. 1986. The metallogesis regularities of the skarn-type lead and zinc deposits in China. Bulletin of the Institute of Mineral Deposits, CAGS, 2: 155–165.
[] Bi CS, Wu JS, Wang MY, Ai SJ. 1988. Stable isotope geology of the Xintianling scheelite deposit. Mineral Deposits, 7(2): 39–48.
[] Campbell A R. 1987. A sulfur isotopic study of the San Cristobal tungsten-base metal mine, Peru. Mineralium Deposita, 22(1): 42–46.
[] Casadevall T, Rye RO. 1980. The tungsten queen deposit, Hamme district, Vance County, North Carolina: A stable isotope study of a metamorphosed quartz-huebnerite vein. Economic Geology, 75(4): 523–537. DOI:10.2113/gsecongeo.75.4.523
[] Che QJ. 2005. Study on Yanshanian magmatic to hydrothermal mineralization processes of major tin-polymetallic deposit-concentrated area in South Hunan. Ph. D. Dissertation. Wuhan: China University of Geoscience, 1-116 (in Chinese with English summary)
[] Chen YR. 1988. Geological and geochemical characteristics and diagenetic: Mineral genetic processes of Yaogangxian granite. Mineral Resources and Geology, 2(1): 62–72.
[] Claypool GE, Hoslter WT, Kaplan IR, Sakai H, Zak I. 1980. The age curves of sulfur and oxygen isotopes in marine sulfate and their mutual interpretation. Chemical Geology, 28: 199–260. DOI:10.1016/0009-2541(80)90047-9
[] Cooke BJ, Godwin CI. 1984. Geology, mineral equilibria, and isotopic studies of the McDame tungsten skarn prospect, north-central British Columbia. Economic Geology, 79(5): 826–847. DOI:10.2113/gsecongeo.79.5.826
[] Doe BR and Zartman HE. 1979. Plumbotectonic, the phanerozoic. In: Barnes HL (ed.). Geochemistry of Hydrothermal Ore Deposits. 2nd Edition. New York: John Wiley, 509-567
[] Feng ZW, Xia WH, Zhang JT, Chen ZY. 1989. Geological characteristics of the Huangsha vein-type tungsten deposit and a discussion on the property of its ore-forming fluid, Jiangxi. Earth Science, 14(4): 423–432.
[] Hu YG, Huang RM, Gong MY. 1984. The lead-zinc deposits related to two metallogenic series in the South Hunan. Hunan Geology, 3(1): 1–13.
[] Industry Standard of Geology and Mineral Resources of the People's Republic of China. 1997. Isotope analysis of geological samples. DZ/T 0184. 1-0184. 2 2-1997. Beijing: Standards Press of China (in Chinese)
[] Ishihara S, Wang PA, Kajiwara Y, Watanabe Y. 2003. Origin of sulfur in some magmatic-hydrothermal ore deposits of South China. Bulletin of the Geological Survey of Japan, 54(3-4): 161–169. DOI:10.9795/bullgsj.54.161
[] Jiang SY, Zhao KD, Jiang YH, Ling HF, Ni P. 2006. A new type of tin mineralization related to granite in South China: Evidence from mineral chemistry, element and isotope geochemistry. Acta Petrologica Sinica, 22(10): 2509–2516.
[] Lei CH, Chen FW, Chen ZH, Xu YM, Gong SQ, Li HQ, Mei YP, Qu WJ, Wang DH. 2010. Petrogenetic and metallogenic age determination of the Huangshaping lead-zind polymetallic deposit and its geological significance. Acta Geoscientia Sinica, 31(4): 532–540.
[] Li YZ. 1986. On the genesis of the Huangshaping Pb-Zn deposit, Guiyang, Hunan in special reference to fluid inclusion studies of its minerals. Acta Mineralogica Sinica, 6(1): 77–85.
[] Liu JQ. 1989. The Xihuashan granite and its mineralization. Bulletin of the Chinese Academy of Geological Sciences, 19: 83–104.
[] Liu WH, Xu WX, Dai TG, Li H. 2006. Isotope geochemistry of the Yejiwei deposit in the Shizhuyuan W-Sn orefield in Hunan Province. Acta Petrologica Sinica, 22(10): 2517–2524.
[] Lu YF, Ma YL, Qu WJ, Mei YP, Chen XQ. 2006. U-Pb and Re-Os isotope geochronology of Baoshan Cu-Mo Polymetallic ore deposit in Hunan Province. Acta Petrologica Sinica, 22(10): 2483–2492.
[] Ohmoto H and Rye RO. 1979. Isotopes of sulfur and carbon. In: Barnes HL (ed.). Geochemistry of Hydrothermal Ore Deposits. 2nd Edition. New York: John Wiley, 509-567
[] Ohmoto H, Lasaga AC. 1982. Kinetics of reactions between aqueous sulfates and sulfides in hydrothermal systems. Geochimica et Cosmochimica Acta, 46(10): 1727–1745. DOI:10.1016/0016-7037(82)90113-2
[] Pen ZA. 1991. Approach on genesis and alteration-mineralization zoning of Jiepailing tin deposit, Hunan. Mineral Resources and Geology, 5(4): 293–299.
[] Qi FL, Zhang Z, Zhu XY, Li YS, Zhen SM, Gong FY, Gong XD, He P. 2012. Skarn geochemistry of the Huangshaping W-Mo polymetallic deposit in Hunan and its geological significance. Geology in China, 39(2): 338–348.
[] Raith JG, Stein HJ. 2000. Re-Os dating and sulfur isotope composition of molybdenite from tungsten deposits in western Namaqualand, South Africa: Implications for ore genesis and the timing of metamorphism. Mnieralium Deposita, 35(8): 741–753. DOI:10.1007/s001260050276
[] Ripley EM, Lightfoot PC, Li CS, Elswick ER. 2003. Sulfur isotopic studies of continental flood basalts in the Noril'sk region: Implications for the association between lavas and ore-bearing intrusions. Geochimica et Cosmochimica Acta, 67(15): 2805–2817. DOI:10.1016/S0016-7037(03)00102-9
[] Robb L. 2005. Introduction to Ore-Forming Processes. MA: Blackwell Science Ltd..
[] Shelton KL, So CS, Rye DM, Park ME. 1986. Geologic, sulfur isotope, and fluid inclusion studies of the Sannae W-Mo Mine, Republic of Korea: Comparison of sulfur isotope systematics in Korean W deposits. Economic Geology, 81(2): 430–446. DOI:10.2113/gsecongeo.81.2.430
[] So CS, Rye DM, Shelton KL. 1983. Carbon, hydrogen, oxygen, and sulfur isotope and fluid inclusion study of the Weolag tungsten-molybdenum deposit, Republic of Korea: Fluid histories of metamorphic and ore-forming events. Economic Geology, 78: 1551–1572. DOI:10.2113/gsecongeo.78.8.1551
[] Song SQ, Hu RZ, Bi XW, Wei WF, Shi SH. 2011. Hydrogen, oxygen and sulfur isotope geochemical characteristics of Taoxikeng tungsten deposit in Chongyi County, southern Jiangxi Province. Mineral Deposits, 30(1): 1–10.
[] Strauss H. 1993. The sulfur isotopic record of Precambrian sulfates: New data and a critical evaluation of the existing record. Precambrian Research, 63(3-4): 225–246. DOI:10.1016/0301-9268(93)90035-Z
[] The National Standard of the People's Republic of China. 1999. Determinations for isotopes of lead、strontium and neodymium in rock samples (GB/T 17672-1999). Beijing: Standards Press of China (in Chinese)
[] Tong QM. 1986. The characteristics of metallogeny of the Huangshaping lead-zinc ore deposit of southern Hunan. Geological Review, 32(6): 565–577.
[] Wang JC, Yao Y, Wang KY, Jiang YC, Shi LC, Liu W. 1993. The characteristic of the silver mineralization in Huangshaping lead-zinc deposit. Geology and Exploration, 29(1): 20–27.
[] Wang RT, Mao JW, Ren XH, He Y, Tang ZL, Wang JY. 2005. Isotopic geochemical evidence of rock-forming and ore-forming processes in the Jianchaling sulfide nickel deposit, Lueyang County, Shaanxi Province. : 513-519.
[] Wang YM, Zhu JA, Yu QH. 1988. The Geology of Lead-zinc Deposits, Hunan Province. Beijing: Geological Publishing House.
[] Xi CZ, Dai TG, Liu WH. 2009. Lead and sulfur isotope geochemistry of the Huangshaping lead-zinc deposit, Hunan Province. Acta Geoscientica Sinica, 30(1): 89–94.
[] Xiao HQ, Zhao KD, Jiang SY, Jiang YH, Ling HF. 2003. Lead isotope geochemistry and ore-forming age of Jinchuantang Sn-Bi deposit in Dongpo ore field, Hunan Province. Mineral Deposits, 22(3): 264–270.
[] Xu WX. 1988. Stable isotope study of several tin deposits in South China. Mineral Resources and Geology, 2(Suppl.1): 128–136.
[] Zhang SG, Ding CG, Li MG, Liu SB. 2009. A genetic mineralogical study of sphalerite in the Fankou Pb-Zn ore district. : 364-374.
[] Zhou T, Liu WH, Li H, Xu WX, Dai TG. 2008. Isotope geochemistry of the Xianghualing tin-polymetallic deposit in Hunan Province. Acta Geoscientica Sinica, 29(6): 703–708.
[] Zhu XY, Wang JB, Zhang Z, Wang YL, Jiang YC. 2010. Identification of the NNW trending ore-controlling structure in the Huangshaping lead-zinc deposit, Hunan Province and its significance to exploration. Geology and Exploration, 46(4): 609–615.
[] Zhu XY, Wang JB, Liu SB, Wang YL, Han Y, Zhen SM and Guo LL. 2013. SHRIMP U-Pb zircon dating of diabase in Fankou lead-zinc deposit, Guangdong, and its geological significance. Acta Geologica Sinica, (1), in press (in Chinese with English abstract)
[] Zhu YL, Li CY, Lin YH. 1981. The Geology of Tungsten Deposits, Southern Jiangxi Province. Nanchang: Jiangxi People's Publishing House: 1-387.
[] Zhuang LC, Lin WS, Xie TH. 1991. Stable isotope geochemistry of the Dajieshan tungsten mine in Jiangxi Province. Bulletin of Yichang Institute of Geology and Mineral Resources, CAGS, 16: 109–120.
[] 毕承思. 1986. 中国矽卡岩铅锌矿床成矿特征的几个问题. 中国地质科学院矿床地质研究所文集, 2: 155–165.
[] 毕承思, 吴静淑, 王美玉, 艾上铮. 1988. 新田岭白钨矿床稳定同位素地质学研究. 矿床地质, 7(2): 39–48.
[] 车勤建. 2005.湘南锡多金属矿集区燕山期岩浆-流体-成矿过程研究.博士学位论文.武汉:中国地质大学, 1-116
[] 陈依壤. 1988. 瑶岗仙花岗岩地质地球化学特征与成岩成矿作用. 矿产与地质, 2(1): 62–72.
[] 冯志文, 夏卫华, 章锦统, 陈紫英. 1989. 江西黄沙脉钨矿床特征及成矿流体性质讨论. 地球科学, 14(4): 423–432.
[] 胡永嘉, 黄日明, 龚茂杨. 1984. 湘南两个内生金属成矿系列的铅锌矿床. 湖南地质, 3(1): 1–13.
[] 蒋少涌, 赵葵东, 姜耀辉, 凌洪飞, 倪培. 2006. 华南与花岗岩有关的一种新类型的锡成矿作用:矿物化学、元素和同位素地球化学证据. 岩石学报, 22(10): 2509–2516.
[] 雷泽恒, 陈富文, 陈郑辉, 许以明, 龚述清, 李华芹, 梅玉萍, 屈文俊, 王登红. 2010. 黄沙坪铅锌多金属矿成岩成矿年龄测定及地质意义. 地球学报, 31(4): 532–540.
[] 李裕祖. 1986. 湖南桂阳黄沙坪铅锌矿床矿物流体包裹体研究及矿床成因讨论. 矿物学报, 6(1): 77–85.
[] 刘家齐. 1989. 西华山花岗岩及其成矿作用. 中国地质科学院院报, 19: 83–104.
[] 刘悟辉, 徐文忻, 戴塔根, 李蘅. 2006. 湖南柿竹园钨锡多金属矿田野鸡尾矿床同位素地球化学研究. 岩石学报, 22(10): 2517–2524.
[] 路远发, 马丽艳, 屈文俊, 梅玉萍, 陈希清. 2006. 湖南宝山铜-钼多金属矿床成岩成矿的U-Pb和Re-Os同位素定年研究. 岩石学报, 22(10): 2483–2492.
[] 彭振安. 1991. 界牌岭锡矿床蚀变矿化分带及其成因探讨. 矿产与地质, 5(4): 293–299.
[] 齐钒宇, 张志, 祝新友, 李永胜, 甄世民, 公凡影, 巩小栋, 何鹏. 2012. 湖南黄沙坪钨钼多金属矿床矽卡岩地球化学特征及其地质意义. 中国地质, 39(2): 338–348.
[] 宋生琼, 胡瑞忠, 毕献武, 魏文凤, 石少华. 2011. 赣南崇义淘锡坑钨矿床氢、氧、硫同位素地球化学研究. 矿床地质, 30(1): 1–10.
[] 童潜明. 1986. 湘南黄沙坪铅-锌矿床的成矿作用特征. 地质论评, 32(6): 565–577.
[] 王静纯, 姚永, 王康遥, 汪元成, 施倪承, 刘卫. 1993. 黄沙坪铅锌矿伴生银矿化特征研究. 地质与勘探, 29(1): 20–27.
[] 王瑞廷, 毛景文, 任小华, 赫英, 汤中立, 汪军谊. 2005. 煎茶岭硫化镍矿床成岩成矿作用的同位素地球化学证据. 地球学报, 26(6): 513–519.
[] 王育民, 朱家鳌, 余琼华. 1988. 湖南铅锌矿地质. 北京: 地质出版社.
[] 息朝庄, 戴塔根, 刘悟辉. 2009. 湖南黄沙坪铅锌多金属矿床铅、硫同位素地球化学特征. 地球学报, 30(1): 89–94.
[] 肖红全, 赵葵东, 蒋少涌, 姜耀辉, 凌洪飞. 2003. 湖南东坡矿田金船塘锡铋矿床铅同位素地球化学及成矿年龄. 矿床地质, 22(3): 264–270.
[] 徐文炘. 1988. 我国南方若干锡矿床成矿物质来源的同位素证据. 矿产与地质, 2(增刊): 128–136.
[] 张术根, 丁存根, 李明高, 刘慎波. 2009. 凡口铅锌矿区闪锌矿的成因矿物学特征研究. 岩石矿物学杂志, 28(4): 364–374.
[] 中华人民共和国地质矿产行业标准. 1997.同位素地质样品分析方法, DZ/T 0184. 1-0184. 2 2-1997.北京:中国标准出版社
[] 中华人民共和国国家标准. 1999.岩石中铅、锶、钕同位素测定方法GB/T 17672-1999.北京:中国标准出版社
[] 周涛, 刘悟辉, 李蘅, 徐文炘, 戴塔根. 2008. 湖南香花岭锡多金属矿床同位素地球化学研究. 地球学报, 29(6): 703–708.
[] 祝新友, 王京彬, 张志, 王艳丽, 江元成. 2010. 湖南黄沙坪铅锌矿NNW向构造的识别及其找矿意义. 地质与勘探, 46(4): 609–615.
[] 祝新友, 王京彬, 刘慎波, 王艳丽, 韩英, 甄世民, 郭宁宁. 2013.广东凡口MVT铅锌矿床成矿年代--来自辉绿岩锆石SHRIMP定年证据.地质学报, (1), 待刊
[] 朱焱龄, 李崇佑, 林运淮. 1981. 赣南钨矿地质. 南昌: 江西人民出版社: 1-387.
[] 庄龙池, 林伟圣, 谢廷焕. 1991.大吉山钨矿的稳定同位素地球化学.中国地质科学院宜昌地质矿产研究所所刊, 第16号: 109-120