Age and sources of the ore-forming material for the Qibaoshan Cu-polymetallic deposit in Hu'nan Province:Evidence from quartz vein Rb-Sr isotopic dating and S-Pb iso-topes
-
摘要:
湘东北七宝山铜多金属矿床位于钦杭成矿带西段,是湘东北规模最大的铜多金属矿床。对七宝山矿床含矿石英脉中的石英矿物进行Rb-Sr同位素年龄测定,获得的Rb-Sr等时线年龄为153.4±2.0Ma(MSWD=1.8),87Sr/86Sr初始值为0.71849±0.00026,与区内石英斑岩形成年龄(153~155Ma)相同,说明七宝山铜多金属矿床成因与石英斑岩体密不可分,成岩成矿年龄均在燕山期。为分析七宝山矿床成矿物质来源,对矿区内的黄铁矿进行了S、Pb同位素分析,矿石δ34S为3.24‰~4.84‰,平均值为4.198‰;岩体δ34S为2.22‰~3.86‰,平均值为2.805‰。δ34S值总体变化较小,岩体δ34S值较矿石小,更趋近于0值,说明岩体中硫极可能源于地幔;而矿床中硫来源于主体地幔硫和少量地壳硫混熔的混合硫源。Pb同位素变化范围也较小,矿石206Pb/204Pb值变化范围为18.315~18.396,平均值为18.359;207Pb/204Pb值变化范围为15.629~15.737,平均值为15.675;208Pb/204Pb值变化范围为38.376~38.856,平均值为38.609。矿化岩体数值与之相似,结果显示,七宝山铜多金属矿床的Pb同位素组成具有下地壳富钍(铅)贫铀(铅)的特点。据此提出,矿床的成矿物质主要来自与含矿斑岩体有联系的深部岩浆分异演化而析出的含矿气-液流体。含矿斑岩体定位-结晶时,通过对周围受热地下水的对流循环作用,可以从围岩中萃取少量成矿物质加入成矿作用。
Abstract:The Qibaoshan Cu-polymetallic deposit, located in the west section of the Qinzhou-Hangzhou metallogenic belt, is the largest Cu-polymetallic deposit in northeastern Hu'nan Province.However, the data of the deposit age are very insufficient.Based on quartz minerals Rb-Sr isotopic dating of the Qibaoshan deposit, the authors tried to obtain the age.The results show that the Rb-Sr isochron age is 153.4±2.0Ma (MSWD=1.8), initial 87Sr/86Sr=0.71849±0.00026, nearly identical with the quartz porphyry (153Ma to 155Ma) in the ore district in age, which suggests that the genesis of the Qibaoshan deposit was closely related to the quartz porphy-ry.The diagenesis and mineralization were both in the Yanshan period.The authors analyzed S, Pb isotope data of the pyrite to iden-tify the mineral sources of the Qibaoshan Cu-polymetallic deposit and obtained the following data:δ34S=3.24‰~4.84‰ for the ore, with an average of 4.198‰; δ34S=2.22‰~3.86‰ for the rock, with an average of 2.805‰.The δ34S values change overall in a small range, and the δ34S values of the rock change even in a smaller range, more close to the value of 0, which suggests that the sul-fur in the rock was most possibly derived from the mantle; and the sulfur in the ore came from the main mantle mixed with a small amount of crustal sulfur.The changes of lead isotopes are also in small ranges:206Pb/204Pb change in the range of 18.315~18.396 for the ore, with an average of 18.359;207Pb/204Pb change in the range of 15.629~15.737, with an average of 15.675;208Pb/204Pb change in the range of 38.376~38.856, with an average of 38.609.The values for the rock change in a similar range.These data show that lead isotope composition of the Qibaoshan Cu-polymetallic deposit has a lower crust (thorium-rich and uranium-depleted) characteris-tics.A comprehensive analysis of the Sr-S-Pb isotopic system of the Qibaoshan deposit indicates that the ore-forming material was derived from the ore gas-liquid fluid, resulting from the deep magma evolution and differentiation in the same evolution process as the quartz porphyry.At the same time, the host rock might have extracted a small amount of minerals from the wall rock to take part in the mineralization process through the circulating convection effect of heating groundwater, resulting in quartz porphyry intrusion and crystallization.
-
近年来,国家整装勘查项目在湘西花垣地区铅锌矿找矿工作取得重大突破,新发现了花垣超大型铅锌矿床(田)。矿床由7个矿段组成,矿体厚度达30m,钻孔见矿率在80%以上。其主矿区铅锌资源储量在1000×104t以上,预计远景储量超过2000×104t。将花垣矿田提升为千万吨级的世界资源基地,潜在价值达数千亿元。有望成为中国最大的铅锌矿基地,此项目将使中国增加一个世界级铅锌矿,有效地缓解资源瓶颈[1]。花垣地区铅锌矿床作为整个湘西地区铅锌矿床的典型代表,其矿床成因或成矿流体来源受到众多地质学者的广泛关注[2-11]。
C、H、O同位素组成是示踪成矿流体及其碳来源的有效手段[12-18],在认识碳酸盐成岩过程中的流体性质方面具有不可替代的作用[19]。热液成因方解石是花垣地区铅锌矿床原生矿石中最主要的脉石矿物,其形成贯穿整个成矿过程,因而研究花垣矿区方解石和闪锌矿的C、H、O同位素组成尤为重要,可提供成矿流体来源方面的重要信息。前人已对花垣铅锌矿床的C、H、O同位素组成进行过研究。夏新阶等[20]对花垣地区李梅锌矿床成矿流体来源进行研究,认为该区成矿溶液中的成矿物质来自寒武纪的海相碳酸盐岩围岩,即容矿层灰岩及其上、下层位的灰岩和白云岩,区内成矿溶液中的水应主要为地层水(包括一些油田水),还应有大量的当地地表水加入。刘文均等[10]对花垣铅锌矿床矿石及脉石矿物的液态包裹体进行成矿流体的H、O同位素组成测定,认为成矿流体来源于建造水,后期可能有雨水渗入。杨绍祥等[21]通过对湘西北铅锌矿床C、H、O同位素特征及成矿环境分析认为,矿床成矿流体是主要来源于深部的热卤水、大气降水和少量变质水的混合流体,成矿流体中的C和O主要来自围岩。蔡应雄等[22]对湘西—黔东下寒武统铅锌矿床的流体包裹体和C、O同位素地球化学特征进行分析,认为花垣铅锌矿床成矿流体可能为区域迁移流体与地层封存水构成的混合流体,成矿物质大部分来源于碳酸盐岩围岩地层。李堃等[23]通过对湘西黔东地区铅锌矿床C、O同位素地球化学特征的研究,认为成矿流体是一种高盐度的低温热卤水。有关成矿流体的来源众说纷纭,是该地区矿床系统研究的薄弱环节,前人在该区采集的用于C、H、O同位素组成研究的样品零散,各矿段系统研究的程度不足,数据少,分析结果统计意义不强,导致得出的结论无法达成一致,在一定程度上制约了对花垣铅锌矿床成因的认识。
在前人研究成果的基础上,本次补充测定了近年新开采地段或新发现矿体的代表性样品,根据矿床的纵向分带特征,整体上系统地分析了花垣地区由北向南分布的团结、李梅、土地坪、蜂塘、大石沟铅锌矿床中方解石的C、O同位素组成及方解石和闪锌矿的H、O同位素组成,探讨了成矿流体中碳的来源,以及区域内矿床成矿流体基本特征的变化趋势和成矿流体来源对矿物C、O同位素组成的影响,以期进一步认识并明确该地区铅锌矿床成矿流体中碳和成矿流体的来源。
1. 区域地质背景
湘西位于上扬子地块东南缘,区内地层发育齐全,除缺失石炭系外,从中元古界冷家溪群,新元古界板溪群、震旦系至古—中生界及第四系均有出露。其中,主要赋矿层位为下寒武统上部、中寒武统下部及下奥陶统下部。储矿地层的岩性均为碳酸盐岩。寒武系—奥陶系碳酸盐岩台地相十分发育,厚约5000m,为铅锌主要含矿岩系[11, 24]。区内的地壳构造运动经历了武陵期、雪峰期—加里东期、海西期、印支期—喜马拉雅期4个发展阶段。总体以北东向褶皱变形和深大断裂为主,深大断裂则为以花垣-张家界断裂、吉首-古丈断裂、麻栗场断裂为主干所组成的断裂带,呈北北东—北东—北东东向弧形展布(图 1),并呈向南西方撒开,往北东方收敛的帚状[11]。这3条断裂为早期控相、后期控矿的深大断裂,花垣铅锌多金属矿集区即局限于这3条断裂及其所控制的下寒武统清虚洞组藻灰岩中[25]。区内大范围未见花岗岩出露,仅于古丈背斜龙鼻咀见基性超基性岩侵入板溪群地层。湘西北的铅锌矿总体呈北东向断续展布,省域内走向长230km,宽50~80km,面积万余平方千米[24]。
图 1 湘西北铅锌矿带区域构造略图(据参考文献[11]修改)1—省界;2—矿田界线;3—向斜及其编号;4—背斜及其编号;5—平移断层及其编号;6—正断层及其编号;7—逆断层及其编号;8—南区斜坡相;9—北区斜坡相;10—花垣地区铅锌矿床Figure 1. Regional tectonic sketch map of the lead-zinc deposits in northwest Hu'nan2. 矿床地质特征
湘西花垣地区团结、李梅、土地坪、蜂塘和大石沟铅锌矿床均为层控矿床,含矿围岩岩性均为灰岩,矿化与热液作用形成的方解石化关系十分密切。铅锌矿体赋存于下寒武统清虚洞组,清虚洞组按岩性及岩相可分为5段,自上而下依次为:白云岩段(1q5)、鲕粒灰岩段(1q4)、藻灰岩段(1q3)、泥晶灰岩段(1q2)和条带灰岩段(1q1)。其中,清虚洞组藻灰岩段(1q3)为区内主要含矿层位,岩性为深灰色厚层状藻灰岩夹粉晶灰岩,藻屑灰岩及斑块状云化灰岩,灰色-浅灰色厚层细-中砂屑灰岩、含藻云质灰岩;鲕粒灰岩段(1q4)为区内次要含矿层位,岩性为藻球粒不等晶灰岩、亮晶含鲕粒碎屑灰岩等。含矿层位为下寒武统清虚洞组下段(藻)灰岩。
矿体形态以层状、似层状、透镜状为主,次为脉状、网脉状。矿体规模大小差异较大,层状、似层状矿体与围岩产状基本一致,偶尔有穿层现象[23]。
铅锌矿矿石的矿物组成较简单,以闪锌矿和方铅矿为主,次为黄铁矿,偶见微量白铁矿。脉石矿物主要为方解石和重晶石,次为少量萤石和沥青。矿化蚀变类型主要为方解石化,常伴有弱黄铁矿化、萤石化、重晶石化和不同程度的褪色重结晶等[11]。
矿石矿物以淡黄色和棕黄色闪锌矿为主,其次为方铅矿和黄铁矿。脉石矿物以方解石为主,次为重晶石、极少量的萤石和微量沥青。方解石常呈条带状或团块状,与铅锌矿石密切共生(图版Ⅰ)。矿石结构以自形-他形晶粒结构为主,此外还有充填或填隙结构、交代结构。闪锌矿以细中粒为主,方铅矿结晶粗大,而黄铁矿主要呈细粒及微细浸染状分布于闪锌矿与围岩接触面附近。矿石构造有浸染状、细粒斑点状、环带状、斑块状及网脉状。围岩蚀变以方解石化、白云石化、硅化为主,次为重晶石化、黄铁矿化、萤石化、沥青化及褪色化[23]。
图版Ⅰa.团结铅锌矿床大脉状方解石化闪锌矿石,闪锌矿沿方解石脉与围岩接触带分布;b.团结铅锌矿床斑脉状方解石化闪锌矿石,闪锌矿沿方解石脉体边缘分布;c.李梅铅锌矿床粗脉状方解石化闪锌矿石,细粒黄铁矿与闪锌矿沿方解石脉与围岩接触带分布;d.蜂塘铅锌矿床斑脉状方解石化闪锌矿石,闪锌矿沿方解石脉边缘分布;e.蜂塘铅锌矿床斑脉状方解石化闪锌矿矿石,闪锌矿沿方解石脉边缘分布;f.土地坪铅锌矿床斑脉状方解石化黄铁矿矿石;g.大石沟铅锌矿床斑脉状方解石化闪锌矿石,闪锌矿和方铅矿沿方解石脉边缘分布;h.大石沟铅锌矿床斑脉状方解石化闪锌矿矿石,细粒黄铁矿与闪锌矿沿方解石脉与围岩接触带分布。Sp—闪锌矿;Gn—方铅矿;Cal—方解石;Py—黄铁矿图版Ⅰ.3. 样品采集及测试方法
方解石和闪锌矿分别是花垣地区铅锌矿床中最重要的脉石矿物和矿石矿物。本次采集的方解石和闪锌矿样品分别来自花垣地区由北而南依次分布的团结、李梅、土地坪、蜂塘和大石沟铅锌矿床(图 1),方解石均形成于主成矿期,与闪锌矿和方铅矿紧密共生,呈块状、粗脉状或斑脉状,闪锌矿多沿方解石脉体边缘分布,或呈斑状、浸染粒状分布于方解石脉体中(图版Ⅰ)。
方解石样品的C、O同位素组成在武汉地质调查中心同位素地球化学实验室测定。首先将方解石样品粉碎至100目,在双目镜下挑选方解石单矿物,重复二次挑选,确保方解石纯度高于98%。挑选的单矿物在研钵中研磨成粉末,过200目筛后在烘箱中烘干2h备用。分析采用100%磷酸法,质谱仪型号为MAT-251,分析误差范围约为0.2‰。
称取30mg试样置于反应管中,并注入4ml 100%磷酸,抽真空2h并稳定在1.0Pa,待试样与磷酸充分混合后,将反应管置于恒温25℃的水中24h,再用液氮吸收CO2气体,纯化后的CO2气体在MAT-251质谱仪上测定C、O同位素组成。分析过程采用标准样品GBW04417和NBS19进行质量监控,分析结果以相对V-PDB的值给出。δ13C以PDB为标准,δ18O以SMOW为标准。计算δ18OSMOW时,采用Friedman等[26]的平衡方程:δ18OSMOW=1.03086×δ18OPDB+30.86。
方解石和闪锌矿样品的H、O同位素测试利用MAT-253完成,H同位素采用加热爆裂法从样品中提取原生流体包裹体中的H2O,将提取的包裹体中的H2O与Zn在400℃条件下反应30min制取H2,在MAT-253质谱仪上测定H2的δD,测试误差约为2‰。O同位素组成分析采用BrF5法,在MAT-253质谱计上测定δ18O值,测定精度为0.2‰。
4. 测试结果
4.1 C、O同位素
表 1为湘西花垣地区铅锌矿床热液方解石C、O同位素组成分析结果,表 2为湘西花垣地区铅锌矿床围岩的C、O同位素组成,可见以下特征。
表 1 湘西花垣地区铅锌矿床主成矿期方解石C、O同位素组成Table 1. C and O isotopic compositions of calcites from the Huayuan Pb-Zn ore deposit in western Hu'nan矿区 原样品号 样品名称 δ13CPDB/‰ δ18OPDB/‰ δ18OSMOW/‰ 团结 11TJ-1B9 斑脉状方解石 0.04 -8.27 22.33 11TJ-1B9 斑脉状方解石 0.05 -8.23 22.38 11TJ-1B11 斑脉状方解石 0.13 -8.13 22.48 11TJ-1B12 斑脉状方解石 -0.24 -8.28 22.32 11TJ-1B14 斑脉状方解石 0.43 -8.75 21.84 13TJ-B1 大脉块状方解石 0.29 -8.54 22.06 13TJ-B10 斑脉状方解石 0.17 -9.18 21.40 13TJ-B13 桁脉状方解石 0.11 -8.46 22.14 13TJ-B18 桁脉状方解石 -0.61 -10.77 19.76 13TJ-B18 桁脉状方解石 -0.61 -10.76 19.77 李梅 13LM-B9 桁脉状方解石 0.88 -9.39 21.18 13LM-B11 大脉块状方解石 1.21 -9.66 20.90 13LM-B13-1 大脉块状方解石 1.12 -10.49 20.05 13LM-B31-1 大脉块状方解石 1.16 -10.28 20.26 13HYC-B5 斑脉状方解石 -2.71 -12.35 18.13 13HYC-B7 斑脉状方解石 -0.50 -10.32 20.22 土地坪 13TDP-B3 桁脉状方解石 -0.90 -12.35 18.13 13TDP-2B3 桁脉状方解石 -1.33 -13.15 17.30 13TDP-2B3 桁脉状方解石 -1.41 -13.23 17.22 蜂塘 13FT-B9 斑脉状方解石 -0.85 -11.17 19.35 13FT-B16 斑脉状方解石 -1.03 -11.93 18.56 13FT-B16 斑脉状方解石 -1.00 -11.85 18.64 13FT-B23 斑脉状方解石 -1.01 -11.18 19.33 13FT-B24 斑脉状方解石 -2.24 -12.47 18.01 大石沟 13DSG-B1 斑脉状方解石 -1.82 -13.05 17.41 13DSG-B3 斑脉状方解石 -1.55 -14.33 16.09 13DSG-B6 桁脉状方解石 -2.60 -13.81 16.62 13DSG-B7 斑脉状方解石 -1.29 -12.11 18.38 注:δ18OSMOW=1.03086×δ18OPDB+30.86[26] 表 2 湘西花垣地区铅锌矿床碳酸盐岩围岩的C、O同位素组成Table 2. C and O isotopic compositions of calcites from the Huayuan Pb-Zn ore deposit in western Hu'nan矿区 原样品号 样品名称 δ13CPDB/‰ δ18OPDB/‰ δ18OSMOW/‰ 数据来源 团结 13NZB-B17 灰岩 0.72 -8.52 22.08 本文 13NZB-B21-1 灰岩 0.92 -7.60 23.03 13NZB-B21-2 灰岩 0.97 -7.48 23.15 13HYC-B21 灰岩 0.15 -9.30 21.27 13HYC-B22 灰岩 0.57 -7.92 22.70 李梅 13LM-B23 灰岩 1.17 -8.83 21.76 本文 13LM-B28 灰岩 0.55 -9.44 21.13 13LM-B30 灰岩 0.62 -9.80 20.76 L16-4 无矿化藻灰岩 0.29 / 22.06 [20, 23] L16-6 无矿化藻灰岩 0.46 / 21.49 LM-11 灰岩 1.05 -9.11 21.47 [22] LM-16 灰岩 0.52 -9.24 21.33 蜂塘 SZS-14 灰岩 0.65 -8.95 21.63 SZS-21 灰岩 0.98 -8.99 21.59 [22] SZS-23 灰岩 0.97 -6.76 23.89 13FT-B30 灰岩 0.19 -7.60 23.03 13FT-B33-1 灰岩 0.69 -8.84 21.75 本文 13FT-B33-2 灰岩 0.65 -8.95 21.63 大石沟 13DSG-B20 灰岩 0.44 -8.47 22.13 13DSG-B21 灰岩 0.55 -10.74 19.79 本文 13DSG-B26 灰岩 0.77 -9.02 21.56 注:δ18OSMOW= 1.03086×δ18OPDB+30.86[26] (1)花垣地区铅锌矿床热液方解石的δ13CPDB值为-2.71‰~1.21‰,平均值为-0.58‰,δ18OSMOW值为16.09‰~22.48‰,平均值为19.72‰。其中,团结矿区方解石的δ13CPDB平均值为-0.024‰,δ18OSMOW平均值为21.648‰。李梅矿区方解石的δ13CPDB平均值为0.193‰,δ18OSMOW平均值为20.123‰。蜂塘矿区方解石的δ13CPDB平均值为-1.226‰,δ18OSMOW平均值为18.778‰。土地坪矿区方解石的δ13CPDB平均值为-1.213‰,δ18OSMOW平均值为17.55‰。大石沟矿区方解石的δ13CPDB平均值为-1.815‰,δ18OSMOW平均值为17.125‰。不同矿床的δ 13CPDB值和δ18OSMOW值较接近,变化范围也较小,组成非常集中,相对均一。花垣北东区域的团结、李梅矿区的δ 13CPDB值和δ18OSMOW值较西南区域的蜂塘、土地坪、大石沟矿区稍高。花垣地区铅锌矿床热液方解石的δ13CPDB值和δ18OSMOW值部分落入海相碳酸盐岩(δ13CPDB=0±4‰,δ18OSMOW=20‰~24‰)[27]范围内,部分介于原生碳酸盐岩和海相碳酸盐岩之间(图 2)。
图 2 湘西花垣地区铅锌矿床成矿期方解石C、O同位素图解(底图据参考文献[28])Figure 2. C-O isotope diagram of the calcite from Pb-Zn ore deposits in western Hu'nan(2)花垣地区铅锌矿床的围岩(藻)灰岩的δ13CPDB值为0.15‰~1.17‰,平均值为0.66‰,δ18OSMOW值为19.79‰~23.89‰,平均值为21.87‰。其中,团结矿区灰岩的δ13CPDB值为0.15‰~0.97‰,平均值为0.67‰,δ18OSMOW值为21.27‰~23.15‰,平均值为22.45‰。李梅矿区(藻)灰岩的δ13CPDB值为0.29‰~ 1.17‰,平均值为0.67‰,δ18OSMOW值为20.76‰ ~ 22.06‰,平均值为21.43‰。花垣蜂塘矿区围岩灰岩的δ13CPDB值为0.19‰~0.98‰,平均值为0.69‰,δ18OSMOW值为21.63‰~23.89‰,平均值为22.25‰。花垣大石沟矿区围岩灰岩的δ13CPDB值为0.44‰~ 0.77‰,平均值为0.59‰,δ18OSMOW值为19.79‰ ~ 22.13‰,平均值为21.16‰。
(3)与热液成因方解石相比,围岩具有更高的δ 13CPDB值和δ 18OSMOW值,围岩的δ 13CPDB值和δ 18OSMOW值落入海相碳酸盐岩(δ 13CPDB=0±4‰,δ 18OSMOW=20‰~24‰)[27]范围内,与围岩为沉积碳酸盐岩的地质特征一致。
4.2 H、O同位素
本次研究选取了27件样品进行H、O同位素分析,包括23件方解石和4件闪锌矿,均与成矿密切相关。测试结果见表 3。
表 3 湘西花垣地区铅锌矿主成矿期成矿流体的D-O同位素组成Table 3. δO-δD composition of ore-forming fluid from typical lead-zinc deposits in Huayuan area, western Hu'nan矿区 样号 矿物 δ18OSMOW/‰ δDSMOW/‰ δOfluid/‰ 换算温度/℃ 团结 13TJ-B1 褐色脉状闪锌矿 / -85.50 4.30 150 13TJ-B2 褐色脉状闪锌矿 / -78.70 6.10 150 13TJ-B8 褐色脉状闪锌矿 / -85.30 -4.10 150 13TJ-B11 黄色脉状闪锌矿 / -91.10 3.30 150 13TJ-B1 方解石 21 -50.6 8.35 150 13TJ-B2 方解石 21.9 -39.4 9.25 150 13TJ-B3 方解石 21.9 -33.4 9.25 150 13TJ-B4 方解石 21.1 -44.3 8.45 150 13TJ-B7 方解石 21.3 -57.6 8.65 150 13TJ-B8 方解石 21.7 -51.6 9.05 150 13TJ-B10 斑脉状方解石 21.4 -50.7 8.75 150 13TJ-B11 方解石 21.8 -41.8 9.15 150 LM-1 斑脉状方解石 / -23 0.25 150 LM-2 斑脉状方解石 / -40 1.54 150 李梅 LM-3 斑脉状方解石 / -15 1.93 150 LM-4 斑脉状方解石 / -46 2.43 150 LM-5 斑脉状方解石 / -33 1.72 150 土地坪 13TDP-B3 粗脉状方解石 17.7 -55.8 3.89 135 13TDP-2B3 粗脉状方解石 17.3 -48.2 3.49 135 蜂塘 13FT-B26 斑脉状方解石 19.3 -44.2 5.49 135 13DSG-B10 块状方解石 17.7 -64.1 3.89 135 13DSG-B11 方解石 18 -54.7 4.19 135 13DSG-B12 方解石 17.1 -48.1 3.29 135 大石沟 13DSG-B13 方解石 16.3 -65.8 2.49 135 13DSG-B15 方解石 17.5 -46.5 3.69 135 13DSG-B16 方解石 16.8 -51.7 2.99 135 13DSG-B17 方解石 17 -46.8 3.19 135 注:δ18OSMOW=1.03086 ×δ18OPDB + 30.86[26];方解石与流体的转换公式采用1000lnα方解石-水=2.78×106T-2-2.89[29];换算温度采用测温结果的峰值,T为绝对温度 23件方解石样品的δDSMOW值介于-65.8‰~ -15‰之间,平均值为-45.75‰,流体中的δ18Ofluid值为0.25‰~9.25‰,平均值为5.02‰。4件闪锌矿样品的δDSMOW值为-91.1‰ ~-78.7‰,平均值为-85.15‰,流体中的δ18Ofluid值为-4.1‰~6.1‰,平均值为2.4‰。
将本次和前人测定的H、O同位素组成(表 3)投在δ18O-δD图解(图 3)上,可见与雨水线斜交的线性关系。李梅矿床的方解石H、O同位素组成更靠近雨水线,矿床成矿流体主要来源于建造水和大气降水,还有部分变质水的混入,而在后期可能有雨水渗入并造成流体盐度的稀释。这种特征也是许多与建造水有关的密西西比型矿床的共同特点[10]。
图 3 湘西花垣地区铅锌矿床成矿流体O-D同位素组成图解(底图据参考文献[30])Figure 3. O-D isotope diagram of the ore-forming fluid from Pb-Zn ore deposits in western Hu'nan5. 讨论
5.1 成矿流体中碳的来源
碳酸盐矿物的C同位素组成对于认识成岩过程中碳酸盐矿物碳的来源至关重要,O同位素组成在示踪碳酸盐成岩过程中流体性质和成岩温度方面具有不可替代的作用[19]。湘西花垣矿区碳酸盐类脉石矿物仅发现方解石,其他碳酸盐类脉石矿物少见,甚至未见,故方解石的δ13C基本代表成矿流体的C同位素组成,可以用方解石中C同位素的组成来指示成矿流体的碳来源。
花垣地区铅锌矿床的围岩及热液方解石样品的δ18O-δ13C同位素图解(图 2)上,热液方解石的δ13CPDB和δ18OSMOW值呈近水平分布,略低于同样呈近水平分布的围岩灰岩。由表 2和图 2可见,花垣地区铅锌矿床中远矿围岩的δ13CPDB值和δ18OSMOW值落入海相碳酸盐岩范围内,表明这些矿床的远矿围岩(灰岩)为沉积成因碳酸盐岩,与其地质特征吻合。不同矿床的热液成因方解石的C、O同位素组成范围较宽,但均处于地幔多相体系与海相碳酸盐岩间,且具有靠近海相碳酸盐岩C同位素组成缓慢升高的趋势(图 3)。表明花垣地区铅锌矿成矿流体中的碳可能来源于:① 沉积有机物的脱羟基作用;② 地幔多相体系的沉积岩混染或高温效应;③ 海相碳酸盐岩的溶解作用。如果碳来源于沉积物中有机质的脱羟基作用,则C同位素组成较有机质升高,O同位素组成较有机质降低(图 2)。花垣地区铅锌矿床热液方解石的C同位素组成明显高于沉积物有机质,O同位素组成明显低于沉积物有机质,沉积有机物的脱羟基作用很难形成如此大的C-O同位素分馏[31]。因此沉积有机物的脱羟基作用应该不是成矿流体中C来源的主要机制。地幔多相体系的沉积岩混染或高温效应会使从其中形成的矿物具有比其更高的δ13C和δ18O值,而花垣地区铅锌矿床中热液方解石的C同位素组成与地幔多相体系相近,但O同位素组成明显高于地幔多相体系(图 3)。另外,该地区铅锌矿床中热液方解石的δ13CPDB值和δ18OSMOW值投点呈近水平分布,如果C、O同位素组成的线性关系是由沉积岩混染作用或高温分异作用所致,则该作用对流体O同位素组成的影响并不明显,而对C同位素组成的影响显著,从而导致从该溶液中沉淀的方解石,其C同位素组成变化显著[12, 32-34]。这与测试结果显然不相符。因此,沉积岩混染或高温效应也不应是成矿流体中碳来源的主要因素。
海相碳酸盐岩的溶解作用是通过流体与围岩之间的水/岩反应,造成δ13CPDB与δ18OSMOW呈相关趋势[32-34],如果成矿流体中的碳形成于海相碳酸盐岩的溶解作用,则其C同位素组成与海相碳酸盐岩相似,但其O同位素组成较海相碳酸盐岩亏损[31]。因此,花垣地区铅锌矿床C、O同位素组成在δ18O-δ13C图解(图 2)中总体呈近水平展布,即因流体与围岩的水/岩反应,脉石矿物热液方解石的沉淀应该主要是由水/岩反应、温度降低耦合等作用所致[28]。因而花垣铅锌矿床成矿流体中的碳可能来源于围岩地层,碳酸盐岩溶解作用起关键作用(亏18O的成矿流体与地层海相碳酸盐岩围岩的水/岩反应)[12]。热液方解石与围岩碳酸盐岩相比,具有明显低的δ18OSMOW值,表明亏18O的成矿流体与围岩发生了同位素交换。
5.2 成矿流体的来源
刘文均等[10]对花垣铅锌矿床中大量流体包裹体进行了详细研究,结果显示,该区成矿流体为高盐度的低温CO2-H2O-NaCl型热卤水,流体密度为1.0~1.1g/cm3,成矿压力为35~40MPa,成矿深度为1.3~1.4km。杨邵祥等[11]通过花垣地区矿床流体包裹体测温,获得成矿流体温度为99~190℃。蔡应雄等[22]对花垣李梅铅锌矿床流体包裹体测温获得成矿流体温度为98~130℃,盐度范围为18.01%~19.95% NaCl eq.,成矿压力为22~40MPa。周云等[1]对花垣李梅、耐子堡、狮子山等典型矿床流体成矿作用的研究认为,该区成矿流体温度主要分布于80~230℃之间,盐度范围为9%~21%NaCl eq.,成矿压力为25~ 45MPa。综合湘西典型铅锌矿流体包裹体特征及数据,该区成矿流体温度主要为100~180℃,总盐度为10%~23%NaCl eq.,多数大于15%NaCl eq.,密度多数大于1g/cm3,流体气相成分以H2O为主,成矿流体体系为NaCl-CaCl2-H2O,CaCl2-H2O,NaClKCl-H2O,CaCl2-MgCl2-H2O-卤水,成分主要为Cl-、Na+、Ca2+、K+、Mg2+,是以钠和钙氯化物为主的高浓度溶液,属于具有中-高盐度、中-高密度、低温度的地下热卤水性质的含矿热水溶液,成矿压力平均值为34~48MPa,成矿深度平均值为1.20~ 1.66km,属浅成低温矿床[1]。前人对全球典型MVT铅锌矿床统计发现,此类矿床流体包裹体均一温度为50~250℃,但多在90~150℃之间(表 4)[35-37],成矿深度多为几百至千余米,包裹体的盐度在10%~30% NaCl eq.之间[38],与油田水组分相似,主要成分从多到少依次为Cl-、Na+、Ca2+、K+、Mg2+,成矿流体是高密度盆地卤水。通过比较发现,湘西地区铅锌矿床成矿流体与美国典型MVT铅锌矿床具有相似性,成矿流体均为低温、中高盐度、高密度性质的盆地卤水,但湘西地区铅锌矿床成矿流体温度整体比美国典型MVT铅锌矿床成矿流体温度稍高。
表 4 中美MVT铅锌矿床成矿温度对比Table 4. Metallogenic temperatures contrast between Chinese and American MVT lead-zinc deposits国家 矿床(区) 主矿物 Th/℃ 资料来源 中
国
湖
南花垣 闪锌矿
方解石120~160
90~180[1] 茶田 闪锌矿
方解石96~170
92~169闪锌矿 113~219 打狗洞 方解石 92~152 石英 85~195 闪锌矿 108~148 董家河 方解石 128~164 石英 100~343 闪锌矿 106~129 唐家寨 石英 100~220 方解石 115~139 美
国
中
部维伯纳姆和老铅带 闪锌矿 90~120 阿肯色北部 闪锌矿 95~170 三州交界地区 闪锌矿 80~120 上密西西比河谷地区 闪锌矿 75~160, 个别达220 [35-37] 田纳西中部 闪锌矿 90~150 田纳西东部 闪锌矿 70~170 肯塔基中部 闪锌矿 70~130 湘西花垣地区铅锌矿床成矿流体温度主要为100~180℃,花垣团结、李梅、土地坪、蜂塘、大石沟铅锌矿床的成矿流体温度依次下降,均一温度具有由北向南降低的趋势,反映了成矿流体的运移方向。刘文均等[10]通过研究也认为,在花垣矿区范围内,流体的温度、CO2、CH4、K+、Na+、Ca2+的含量,出现从矿区北段向南段、从矿区到外围同步下降的特点,表明在花垣矿区范围内流体的迁移方向为由北向南流动,且主要途径是沿清虚洞组三、四段间顺层流动。成矿流体在运移和沉淀过程中依次形成团结、李梅、土地坪、蜂塘、大石沟铅锌矿床。
花垣地区铅锌矿床的H、O同位素组成表明,成矿流体的主要来源是建造水和大气降水。自然界中天然水富集12C、16O等轻同位素而贫重同位素,尤其是18O的含量特别低,而沉积碳酸盐岩则富含13C、18O等重同位素[23]。因此,主要来源为建造水和大气降水的成矿流体相对集中较多的轻同位素,亏13C、18O的成矿流体与围岩地层的水-岩反应中,反复与碳酸盐围岩进行同位素交换,导致成矿期沉淀的方解石中13C、18O逐渐降低。花垣矿区范围内流体的迁移方向是由北向南,这也是团结、李梅、土地坪、蜂塘、大石沟铅锌矿床中成矿期方解石的13C、18O同位素逐渐降低的原因(表 1;图 2)。成矿流体与围岩的水-岩反应是导致湘西花垣地区铅锌矿床中方解石矿物沉淀的主要机制。
6. 结论
(1)湘西花垣地区铅锌矿床围岩为海相沉积成因碳酸盐岩,成矿期方解石中的碳主要来源于碳酸盐岩围岩,由围岩海相碳酸盐岩的溶解作用形成。该地区铅锌矿床脉石矿物热液方解石和矿石矿物闪锌矿的沉淀主要是成矿流体与围岩的水/岩反应、温度降低耦合等作用所致,亏18O的成矿流体与围岩发生了同位素交换。
(2)花垣地区铅锌矿床成矿流体的主要来源是建造水和大气降水,以及少量变质水。亏13C、18O的成矿流体与围岩地层发生水-岩反应,由北向南分布的铅锌矿床中成矿期方解石的13C、18O同位素表现出逐渐降低的特征,矿区范围内成矿流体的迁移方向是由北向南流动,可根据这一规律指导下一步找矿勘探。
致谢: 审稿专家对稿件提出了宝贵的修改意见,在此表示诚挚的谢意。 -
图 4 七宝山铜多金属矿床Pb同位素模式图(底图据参考文献[43])
Figure 4. The Pb isotopic model for the Qibaoshan Cu-polymetallic deposit
表 1 七宝山铜多金属矿床含矿石英脉中的石英矿物Rb-Sr同位素测定结果
Table 1 Rb-Sr isotopic data of fluid inclusions in quartz mineral of the Qibaoshan Cu-polymetallic deposit
样号 样品名称 Rb/10-6 Sr/10-6 87Rb/86Sr 87Sr/86Sr 2σ (87Sr/86Sr)i QB22-1 石英 1.335 0.4996 7.724 0.73554 0.00020 0.71863 QB22-2 石英 0.622 0.3968 4.524 0.72829 0.00006 0.71839 QB22-3 石英 1.004 0.4501 6.450 0.73254 0.00018 0.71842 QB22-4 石英 1.143 0.3416 9.674 0.73945 0.00010 0.71827 QB22-5 石英 2.408 0.4854 14.360 0.74981 0.00016 0.71837 表 2 七宝山铜多金属矿床硫化物的S同位素组成
Table 2 Surfur isotopic composition of pyrite from the Qibaoshan Cu-polymetallic deposit
编号 样品描述 采样位置 分析矿物 δ34SCDT/‰ QB2 块状黄铁矿矿石 老虎口矿段40m中段之33m分段13线 黄铁矿 4.50 QB4 块状黄铁矿磁铁矿矿石 老虎口矿段40m中段之33m分段93线 黄铁矿 4.11 QB6 块状黄铁矿矿石 老虎口矿段40m中段之33m分段7线 黄铁矿 4.45 QB6 块状黄铁矿矿石 老虎口矿段40m中段之33m分段7线 黄铁矿 4.68 QB8 块状桁粒黄铁矿矿石 老虎口矿段40m中段5线 黄铁矿 4.29 QB9 块状细粒黄铁矿矿石 老虎口矿段40m中段5线 黄铁矿 4.84 QB12 含铜磁铁矿矿石 鸡公湾矿段80m中段12线 黄铁矿 3.24 QB19 含铜黄铁矿矿石 大七宝山矿段146m矿段28线 黄铁矿 3.47 QB1 黄铁矿化石英斑岩 老虎口矿段40m中段之33m分段13线 黄铁矿 2.84 QB17 细脉侵染状黄铁矿化石英斑岩 大七宝山矿段80m矿段26线 黄铁矿 3.86 QB20 闪锌矿化石英斑岩 大七宝山矿段146m矿段28线 黄铁矿 2.22 QB20 闪锌矿化石英斑岩 大七宝山矿段146m矿段28线 黄铁矿 2.30 注:主量元素含量单位为%,微量、稀土元素为10-6 表 3 七宝山铜多金属矿床硫化物的Pb同位素组成
Table 3 Lead isotopic composition of pyrite from the Qibaoshan Cu-polymetallic deposit
样品类型 分析矿物 分析结果 样号 同位素比值 表面年龄/Ma φ值 μ值 Th/U 208Pb/(206Pb+207Pb) 206Pb/204Pb 207Pb/204Pb 208Pb/204Pb QB2 黄铁矿 18.384±0.007 15.689±0.007 38.509±0.015 295 0.594 9.64 3.79 1.1302 QB4 黄铁矿 18.315±0.002 15.661±0.002 38.376±0.005 310 0.596 9.59 3.77 1.1295 QB6 矿石 黄铁矿 18.384±0.002 15.737±0.003 38.856±0.007 351 0.600 9.74 3.95 1.1388 QB8 黄铁矿 18.326±0.006 15.665±0.004 38.632±0.015 307 0.596 9.6 3.87 1.1365 QB12 黄铁矿 18.350±0.001 15.629±0.001 38.577±0.004 247 0.590 9.53 3.83 1.1353 QB19 黄铁矿 18.396±0.003 15.666±0.003 38.705±0.008 258 0.591 9.6 3.87 1.1363 QB1 石 黄铁矿 18.390±0.003 15.711±0.001 38.571±0.006 316 0.596 9.69 3.82 1.1311 QB17 英 黄铁矿 18.318±0.004 15.652±0.003 38.569±0.01 297 0.595 9.58 3.85 1.1354 QB20 斑 黄铁矿 18.412±0.006 15.717±0.006 38.734±0.015 308 0.596 9.69 3.88 1.1349 QB20 岩 黄铁矿 18.373±0.007 15.707±0.006 38.728±0.016 324 0.597 9.68 3.89 1.1364 表 4 钦杭结合带主要斑岩相关铜矿床及其成岩成矿年龄
Table 4 Porphyry-related copper deposits and their ages of diagenesis and mineralization in QHSZ, southern China
序号 矿床名称 矿石种类 侵人体类型 侵人体年龄 成矿年龄 资料来源 测试方法 年龄/Ma 测试方法 年龄/Ma 1 德兴 Cu-Au-Mo 花岗闪长斑岩 锆石U-Pb 171±3 辉钼矿Re-Os 170.4±1.8 [30-31] 2 永平 Cu-W-Mo 石英斑岩 锆石U-Pb 160±2.3
135±7.4辉钼矿Re-Os 156.7±2.8
155.7±3.6[32] [33] 3 村前 Cu多金属 花岗斑岩 锆石U-Pb 169±1.1 [34] 4 铜山岭 Cu多金属 花岗闪长斑岩 SHRIMP锆石U-Pb 149±4 [35] 5 宝山 Cu-Mo-W 花岗闪长斑岩 锆石U-Pb 173±1.9 黄铁矿Rb-Sr 174±7 [4, 36] 6 水口山 Cu多金属 花岗闪长斑岩 锆石U-Pb 172.3±1.6 [4] 7 七宝山 Cu多金属 石英斑岩 锆石U-Pb 155~153 石英Rb-Sr 153.4±2.0 [23], 本文 8 圆珠顶 Cu-Mo 花岗斑岩 SHRIMP锆石U-Pb 154±2
153.4±1.6辉钼矿Re-Os 155±5 [37-38] -
杨明桂, 黄水保, 楼法生, 等.中国东南陆区岩石圈结构与大规模成矿作用[J].中国地质, 2009, 36(3):528-543. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200903006.htm 毛景文, 谢桂青, 郭春丽, 等.华南地区中生代主要金属矿床时空分布规律和成矿环境[J].高校地质学报, 2008, 14(4):510-526. http://www.cnki.com.cn/Article/CJFDTOTAL-GXDX200804007.htm 毛景文, 陈懋弘, 袁顺达, 等.华南地区钦杭成矿带地质特征和矿床时空分布规律[J].地质学报, 2011, 85(5):636-658. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201105006.htm 王岳军, 范蔚茗, 郭峰, 等.湘东南中生代花岗闪长岩锆石U-Pb法定年及其成因指示[J].中国科学(D辑), 2001, 31(9):745-751. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200109005.htm 陆玉梅, 殷浩然, 沈瑞锦.七宝山多金属矿床成因模式[J].矿床地质, 1984, 3(4):53-60. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ198404006.htm 沈瑞锦, 陆玉梅. 湖南七宝山多金属矿床成因模式[C]//钨家达. 湖南地学新进展. 长沙: 湖南科学技术出版社, 1996: 29-36. 胡祥昭, 彭恩生, 孙振家.湘东北七宝山铜多金属矿床地质特征及成因探讨[J].大地构造与成矿学, 2000, 24(4):365-370. http://cdmd.cnki.com.cn/Article/CDMD-10533-2002080803.htm 胡祥昭, 肖宪国, 杨中宝.七宝山花岗斑岩的地质地球化学特征[J].中南工业大学学报, 2002, 33(6):551-554. http://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD200206001.htm 胡俊良, 徐德明, 张鲲.湖南七宝山石英斑岩地球化学特征及其与成矿的关系[J].华南地质与矿产, 2012, 28(4):298-306. http://www.cnki.com.cn/Article/CJFDTOTAL-HNKC201204002.htm 刘姤群, 金维群, 张录秀, 等.湘东北斑岩型和热液脉型铜矿成矿物质来源探讨[J].华南地质与矿产, 2001, 17(1):40-47. http://www.cnki.com.cn/Article/CJFDTOTAL-HNKC200101005.htm 李华芹, 陈富文, 蔡红.新疆西准噶尔地区不同类型金矿床RbSr同位素年代研究[J].地质学报, 2000, 74(2):181-192. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200002008.htm 李华芹, 王登红, 陈富文, 等.湖南雪峰山地区铲子坪和大坪金矿成矿作用年代学研究[J].地质学报, 2008, 82(7):900-905. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200807007.htm 蔡明海, 韩凤彬, 何龙清, 等.湘南新田岭白钨矿床He, Ar同位素特征及Rb-Sr测年[J].地球学报, 2008, 29(2):167-173. http://www.cnki.com.cn/Article/CJFDTOTAL-DQXB200802007.htm 付建明, 李华芹, 马丽艳, 等.粤北东昌市和尚田钨锡多金属矿成矿时代及其地质意义[J].地质学报, 2013, 87(9):1349-1358. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201309012.htm 蔺志永, 王登红, 张长青.四川宁南跑马铅锌矿床的成矿时代及其地质意义[J].中国地质, 2010, 37(2):488-494. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201002024.htm Faure G. Principle of Isotope Geology(2nd Edition)[M]. New York:John Wiley & Sons, 1986:183-199.
胡俊良, 徐德明, 张鲲.湖南七宝山矿床石英斑岩锆石U-Pb定年及Hf同位素地球化学[J].矿床地质, 2014, 33(增刊):201-202. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ2014S1103.htm 李华芹, 刘家齐, 魏琳.热液矿床流体包裹体年代学研究及其地质应用[M].北京:地质出版社, 1993:1-27. 李华芹.新疆北部有色金属矿床成矿作用年代学[M].北京:地质出版社, 1998. 蔡应雄, 杨红梅, 段瑞春, 等.湘西-黔东下寒武统铅锌矿床流体包裹体和硫、铅、碳同位素地球化学特征[J].现代地质, 2014, 28(1):29-41. http://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ201401004.htm Qiu X F, Yang H M, Lu S S, et al. Geochronology and geochemis-try of Grenville-aged(1063±16Ma) metabasalts in the Shennongjia district, Yangtze block:implications for tectonic evolution of the South China Craton[J]. International Geology Review, 2015, 57:76-96. doi: 10.1080/00206814.2014.991949
Ludwig K R. IsoplotEx v. 2.6[M]. Berkeley Geochronological Center Spec Publ. LA. 1999.
胡俊良, 徐德明, 张鲲, 等.湖南七宝山铜多金属矿床石英斑岩时代与成因:锆石U-Pb定年及Hf同位素与稀土元素证据[J].大地构造与成矿学, 2016, 40(6):1196-1211. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201606007.htm Norman D Z, Lands G P. Source of mineralizing components in hydrothermal ore fluids as evidenced by 87Sr/86Sr and stable isotope data from the Pasto Bueno Deposit, Peru[J]. Economical Geology, 1983, 78:451-456. doi: 10.2113/gsecongeo.78.3.451
Rossman G R, Weis D, Wasserburg G J. Rb, Sr, Nd and Sm con-centrations in quartz[J]. Geochimica et Cosmochimica Acta, 1987, 51:2325-2329. doi: 10.1016/0016-7037(87)90286-9
Changkakoti A, Gray J, Krstic D, et al. Determimations of radio-genic isotopes (Rb/Sr, Sm/Nd, and Pb/Pb) in fluid inclusion wa-ters:An example from the Bluebell Pb-Zn deposit, British Colum-bia, Canada[J]. Geochimica et Cosmochimica Acta, 1988, 52:961-967. doi: 10.1016/0016-7037(88)90251-7
张理刚.稳定同位素在地质科学中的应用[M].西安:陕西科学技术出版社, 1985:152-185. 郑永飞, 陈江峰.稳定同位素地球化学[M].北京:科学出版社, 2000:143-245. 张长青, 李厚民, 代军治, 等.铅锌矿床中矿石铅同位素研究[J].矿床地质, 2006, 25(增刊):213-216. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200611003057.htm 王强, 赵振华, 简平, 等.德兴花岗闪长斑岩SHRIMP锆石U-Pb年代学和Nd-Sr同位素地球化学[J].岩石学报, 2004, 20(2):315-324. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200402011.htm Lu J J, Hua R M, Yao C L. Re-Os age for molybdenite form the Dexing porphyry Cu-Au deposit in Jiangxi province, China[J]. Geochimica et Cosmochimica Acta, 2005, 69(Suppl. A):882. https://www.researchgate.net/publication/270455478_Multifractal_distribution_of_geochemical_elements_in_Dexing_porphyry_copper_deposit_Jiangxi_Province_China
丁昕, 蒋少涌, 倪培, 等.江西武山和永平铜矿含矿花岗质岩体锆石SIMS U-Pb年代学[J].高校地质学报, 2005, 11(3):383-389. http://www.cnki.com.cn/Article/CJFDTOTAL-GXDX200503009.htm 李晓峰, Yasushi W, 华仁民, 等.华南地区中生代Cu-(Mo)-WSn矿床成矿作用于洋岭/转换断层俯冲[J].地质学报, 2008, 82(5):625-640. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200805007.htm 王强, 孙燕, 张雪辉, 等.江西省村前铜多金属矿床斜长花岗斑岩LA-ICP-MS锆石U-Pb年龄及地质意义[J].中国地质, 2012, 39(5):1143-1150. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201205002.htm 魏道芳, 鲍征宇, 付建明.湖南铜山岭花岗岩体的地球化学特征及锆石SHRIMP定年[J].大地构造与成矿学, 2007, 31(4):482-489. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200704014.htm 姚军明, 华仁民, 林锦富.湘南宝山矿床REE、Pb-S同位素地球化学及黄铁矿Rb-Sr同位素定年[J].地质学报, 2006, 80(7):1045-1054. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200607010.htm 陈富文, 李华芹, 王登红, 等.粤西圆珠顶斑岩型铜钼矿床成矿地质特征及成岩成矿作用年代学研究[J].地质学报, 2012, 86(8):1298-1305. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201208013.htm 陈懋弘, 李忠阳, 李菁, 等.初论广西大瑶山地区多期次花岗质岩浆活动与成矿系列[J].地学前缘, 2015, 22(2):41-53. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201502005.htm 侯明兰, 蒋少涌, 姜耀辉, 等.胶东蓬莱金城矿区的S-Pb同位素地球化学和Rb-Sr同位素年代学研究[J].岩石学报, 2006, 22(10):2525-2533. doi: 10.3969/j.issn.1000-0569.2006.10.013 路远发. GeoKit:一个用VBA构建的地球化学工具软件包[J].地球化学, 2004, 33(5):459-464. http://www.cnki.com.cn/Article/CJFDTOTAL-DQHX200405003.htm 孙省利. 西秦岭泥盆系西城矿化集中区烃碱流体成矿系列研究[D]. 成都理工大学博士学位论文, 2001: 43-44. 刘姤群, 杨世义, 张秀兰, 等.粤北大宝山多金属矿床成因的初步探讨[J].地质学报, 1985, 61(1):47-61. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE198501005.htm Zartman R E, Doe B R. Plumbotectonics-the model[J]. Tectono-physics, 1981, 75:135-162. doi: 10.1016/0040-1951(81)90213-4
中国地质调查局资源评价部等. 钦杭成矿带重要矿产勘查部署方案. 2010.