Metallogenesis of the main collisional period in mid-western Gangdise:Zir-con U-Pb geochronology of the granite porphyry in Dexin deposit, Tibet
-
摘要:
德新银铅锌矿床位于拉萨地块中部隆格尔-工布江达断隆带中段,是冈底斯银铅锌多金属成矿带上典型的热液型矿床之一。区内的铅锌矿化主要赋存于围岩裂隙和花岗斑岩中。利用LA-ICP-MS锆石U-Pb测年方法对区内含矿花岗斑岩进行年龄测定,获得206Pb/238U年龄加权平均值为57.7±0.5Ma(MSWD=1.6,n=15),代表花岗斑岩及与之相关的铅锌矿化形成年龄。综合区域地质背景,认为德新矿区与花岗斑岩有关的铅锌矿化形成于印度-亚洲大陆主碰撞时期,与纳如松多隐爆角砾岩型铅锌矿属同一岩浆成矿系统,冈底斯中西段存在(近)主碰撞阶段成矿事件。
Abstract:The Dexin Ag-Pb-Zn deposit, located in the Longgeer-Gongbujiangda fault-block uplift belt at the center of Lhasa Ter-rane, is one of the representative hydrothermal type deposits in the Gangdise Ag-Pb-Zn polymetallic ore-forming belt.The Pb-Zn mineralization dominantly occurs in the fractures of wall rocks and the granitic porphyries.Zircons from the ore-bearing granitic por-phyry in Dexin were measured by the LA-ICP-MS U-Pb dating method, which yielded a weighted average 206Pb/238U age of 57.7±0.5Ma (MSWD=1.6, n=15) that represents the formation age of the granitic porphyries and the Pb-Zn mineralization in the granitic porphyries.With a comprehensive analysis of geological background, it is concluded that the Pb-Zn mineralization associated with the granitic porphyries in Dexin and the cryptoexplosive breccia type Pb-Zn mineralization in Narusongduo were formed by the same magmatic metallogenic system, and there were many magmatic-metallogenic events near/in the main Indian-Asian collisional period in the mid-west part of Gangdise.
-
Keywords:
- granitic porphyry /
- zircon U-Pb age /
- Dexin /
- Tibet
-
近年来,国家整装勘查项目在湘西花垣地区铅锌矿找矿工作取得重大突破,新发现了花垣超大型铅锌矿床(田)。矿床由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同位素表现出逐渐降低的特征,矿区范围内成矿流体的迁移方向是由北向南流动,可根据这一规律指导下一步找矿勘探。
致谢: 野外工作得到西藏金泰矿业公司徐仲勋总工程师、西藏地勘局第五地质大队驻德新矿区工作人员及“西藏则学地区矿产远景调查”项目组成员的热情指导和帮助,锆石U-Pb定年工作得到中国地质大学(武汉)地质过程与矿产资源国家重点实验室胡兆初老师的帮助,在此一并表示诚挚的感谢。 -
图 2 德新矿区地质图①
Q—第四系;E2n—年波组火山岩;C2-P1l2—拉嘎组上段含砾砂板岩;C2-P1l1—拉嘎组下段粉砂质板岩;K2δ—晚白垩世闪长岩;δμ—闪长玢岩;1—正断层;2—逆断层;3—地质界线;4—含矿花岗斑岩;5—银铅锌多金属矿体;6—采样位置
Figure 2. Geological map of the Dexin Ag-Pb-Zn deposit, Tibet
表 1 西藏德新矿区花岗斑岩LA-ICP-MS锆石U-Th-Pb同位素分析结果
Table 1 Zircon LA-ICP-MS U-Th-Pb data of the granite porphyry in Dexin deposit, Tibet
样品编号 Th/10-6 U/10-6 Th/U 同位素比值 同位素年龄/Ma 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ DX1-QP01-01 3042 8204 0.37 0.04874 0.00207 0.06110 0.00252 0.00911 0.00012 135 91 60 2 58.5 0.7 DX1-QP01-02 3151 8263 0.38 0.04548 0.00170 0.05584 0.00213 0.00889 0.00011 -30 75 55 2 57.1 0.7 DX1-QP01-03 2542 6910 0.37 0.04798 0.00198 0.05962 0.00245 0.00898 0.00010 98 89 59 2 57.6 0.7 DX1-QP01-04 2174 7951 0.27 0.04678 0.00226 0.05612 0.00266 0.00866 0.00011 38 99 55 3 55.6 0.7 DX1-QP01-05 3394 9888 0.34 0.04785 0.00180 0.05968 0.00240 0.00895 0.00013 92 85 59 2 57.4 0.8 DX1-QP01-06 2928 7086 0.41 0.04798 0.00202 0.06019 0.00247 0.00911 0.00011 98 88 59 2 58.5 0.7 DX1-QP01-07 2849 6044 0.47 0.04740 0.00194 0.05883 0.00244 0.00895 0.00011 69 88 58 2 57.4 0.7 DX1-QP01-08 3307 8903 0.37 0.04898 0.00180 0.06064 0.00222 0.00895 0.00012 147 80 60 2 57.5 0.8 DX1-QP01-09 2781 5510 0.50 0.04846 0.00239 0.06107 0.00302 0.00915 0.00011 122 109 60 3 58.7 0.7 DX1-QP01-10 1567 1474 1.06 0.05111 0.00433 0.06390 0.00494 0.00933 0.00018 246 172 63 5 60.0 1.0 DX1-QP01-11 3641 8126 0.45 0.04732 0.00189 0.05905 0.00239 0.00900 0.00012 66 84 58 2 57.8 0.8 DX1-QP01-12 5313 12458 0.43 0.04611 0.00183 0.05683 0.00239 0.00885 0.00014 4 82 56 2 56.8 0.9 DX1-QP01-13 2698 8829 0.31 0.04500 0.00172 0.05562 0.00211 0.00893 0.00011 -19 74 55 2 57.3 0.7 DX1-QP01-14 3129 7990 0.39 0.04814 0.00186 0.06072 0.00235 0.00906 0.00010 106 84 60 2 58.2 0.7 DX1-QP01-15 1027 885 1.16 0.04923 0.00534 0.05967 0.00629 0.00908 0.00022 159 230 59 6 58.0 1.0 -
孟祥金. 西藏碰撞造山带冈底斯中新世斑岩铜矿成矿作用研究[D]. 中国地质科学院博士学位论文, 2004. http://www.oalib.com/references/17373755 孟祥金, 侯增谦, 叶培盛, 等.西藏冈底斯银多金属矿化带的基本特征与成矿远景分析[J].矿床地质, 2007, 26(2):153-162. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200702001.htm 侯增谦, 曲晓明, 杨竹森, 等.青藏高原碰撞造山带:Ⅲ.后碰撞伸展成矿作用[J].矿床地质, 2006, 25(6):629-651. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200606000.htm 孟祥金, 侯增谦, 高永丰, 等.西藏冈底斯东段斑岩铜钼铅锌成矿系统的发育时限:帮浦铜多金属矿床辉钼矿Re-Os年龄证据[J].矿床地质, 2003, 22(3):246-252. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200303003.htm 唐菊兴, 陈毓川, 王登红, 等.西藏工布江达县沙让斑岩钼矿床辉钼矿铼-锇同位素年龄及其地质意义[J].地质学报, 2009, 83(5):698-704. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200905010.htm 高一鸣, 陈毓川, 唐菊兴, 等.西藏工布江达县亚贵拉铅锌钼多金属矿床石英斑岩锆石SHRIMP定年及其地质意义[J].地质学报, 2009, 83(10):1436-1444. doi: 10.3321/j.issn:0001-5717.2009.10.008 Zheng Y C, Fu Q, Hou Z Q, et al. Metallogeny of the northeastern Gangdese Pb-Zn-Ag-Fe-Mo-W polymetallic belt in the Lhasa terrane, southern Tibet[J]. Ore Geology Reviews, 2015, 70:510-532. doi: 10.1016/j.oregeorev.2015.04.004
Zhao J X, Qin K Z, Li G M, et al. Collision-related genesis of the Sharang porphyry molybdenum deposit, Tibet:evidence from zircon U-Pb ages, Re-Os ages and Lu-Hf isotopes[J]. Ore Geology Re-views, 2014, 56:312-326. doi: 10.1016/j.oregeorev.2013.06.005
付强, 杨竹森, 郑远川, 等.西藏龙马拉Cu-Fe-Pb-Zn多金属矿床金云母Ar-Ar定年及其地球动力学意义[J].岩石矿物学杂志, 2014, 33(2):283-293. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW201402007.htm 纪现华, 孟祥金, 杨竹森, 等.西藏纳如松多隐爆角砾岩型铅锌矿床绢云母Ar-Ar定年及其地质意义[J].地质与勘探, 2014, 50(2):281-290. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201402008.htm 付强, 杨竹森, 郑远川, 等.加拉普铁矿区花岗闪长岩锆石U-Pb年龄、Hf同位素及地球化学研究[J].矿床地质, 2013, 32(3):564-578. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ201303008.htm 辛存林, 都卫东, 张育平, 等.西藏洞中松多铅锌多金属矿床地质特征及成因[J].兰州大学学报(自然科学版), 2013, 49(1):24-31. http://www.cnki.com.cn/Article/CJFDTOTAL-LDZK201301006.htm 黄克贤, 郑远川, 张松, 等.西藏亚贵拉矿区两期岩体LA-ICPMS锆石U-Pb定年及地质意义[J].岩石矿物学杂志, 2012, 31(3):348-360. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW201203006.htm 付强, 郑远川, 黄克贤.西藏龙马拉多金属矿床硫、铅同位素组成及其地质意义[J].有色金属(矿山部分), 2012, 64(4):26-30. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKU201204009.htm 纪现华, 杨竹森, 于玉帅, 等.西藏纳如松多铅锌矿床成矿岩体形成机制:岩浆锆石证据[J].矿床地质, 2012, 41(4):758-774. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ201204009.htm 于玉帅, 杨竹森, 多吉, 等.西藏加多捕勒铁铜矿成矿岩体时代与成因:锆石U-Pb年龄、Hf同位素与稀土元素证据[J].矿床地质, 2011, 30(3):420-434. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ201103005.htm 李应栩, 谢玉玲, 陈伟, 等.西藏恰功铁矿二长花岗斑岩锆石的U-Pb年代学与地球化学特征及意义[J].岩石学报, 2011, 27(7):2023-2033. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201107012.htm 高一鸣, 陈毓川, 唐菊兴, 等.西藏工布江达地区亚贵拉铅锌钼矿床辉钼矿Re-Os测年及其地质意义[J].地质通报, 2011, 30(7):1027-1036. http://dzhtb.cgs.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20110704&journal_id=gbc 费光春, 温春齐, 周雄, 等.西藏洞中拉铅锌矿床石英激光探针40Ar-39Ar定年及地质意义[J].矿物岩石, 2010, 30(3):38-43. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=kwys201003005&dbname=CJFD&dbcode=CJFQ 杨勇, 罗泰义, 杨竹森, 等.冈底斯造山带两套不同成矿体系的含矿斑岩对比研究[J].矿床地质, 2010, 29(2):195-206. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ201002002.htm 刘波, 董随亮, 李光明, 等. 西藏斯弄多铅锌矿岩体地球化学特征及年龄地质意义[J]. 2010, 29(增刊): 472-473. 李光明, 刘波, 董随亮, 等. 冈底斯成矿带中段铜铁铅锌矿集区的叠合成矿作用及意义[J]. 2010, 29(增刊): 217-218. 谢玉玲, 李应栩, Cooke D R, 等.西藏恰功铁矿岩浆演化序列及斑岩出溶流体特征[J].地质学报, 2009, 83(12):1869-1886. doi: 10.3321/j.issn:0001-5717.2009.12.006 张林奎, 范文玉, 高大发.西藏林周县勒青拉铅锌多金属矿床地质特征及成因[J].地质与勘探, 2008, 44(5):10-16 http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT200805004.htm 朱斯豹, 汪雄武, 彭惠娟, 等.西藏索波贡巴地区中新世花岗岩类与矽卡岩型铅锌矿的关系[J].矿床地质, 2010, 29(增刊):365-366. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ2010S1187.htm 刘英超, 纪现华, 侯增谦, 等.一个与岩浆作用有关的独立铅锌成矿系统的建立——以西藏纳如松多铅锌矿床为例[J].岩石矿物学杂志, 2015, 34(4):539-556. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW201504008.htm 臧文栓, 孟祥金, 杨竹森, 等.西藏冈底斯成矿带铅锌银矿床的S、Pb同位素组成及其地质意义[J].地质通报, 2007, 26(10):1393-1397. doi: 10.3969/j.issn.1671-2552.2007.10.017 Yin A, Harrison T M. Geologic evolution of the Himalayan Tibet-an Orogen[J]. Earth and Planetary Science Letters, 2000, 28:211-280. doi: 10.1146/annurev.earth.28.1.211
侯增谦, 孟祥金, 曲晓明.西藏冈底斯斑岩铜矿带埃达克质斑岩含矿性:源岩相变及深部过程约束[J].矿床地质, 2005, 24(2):108-121. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200502002.htm Liu Y S, Hu Z C, Gao S. In situ analysis of major and trace ele-ments of anhydrous minerals by LA-ICP-MS Without applying an internal standards[J]. Chemical Geology, 2008, 257:34-43. doi: 10.1016/j.chemgeo.2008.08.004
Liu Y S, Gao S, Hu Z C. Continental and oceanic crust recyclinginduced melt-peridotite interactions in the Trans-North China Orogen:U-Pb dating, Hf isotopes and trace elements in zircons of mantle xenoliths[J]. Journal of Petrology, 2010, 51:537-571. doi: 10.1093/petrology/egp082
Liu Y S, Hu Z C, Zong K Q, et al. Reappraisement and refine-ment of zircon U-Pb isotope and trace element analyses by LAICP-MS[J]. Chinese Science Bulletin, 2010, 55(15):1535-1546. doi: 10.1007/s11434-010-3052-4
Andersen T. Correction of common lead in U-Pb analyses that do not report 204Pb[J]. Chemical Geology, 2002, 192:59-79. doi: 10.1016/S0009-2541(02)00195-X
Ludwig K R. User's Manual for ISOPLOT 3.00:A Geochronolog-ical Toolkit for Microsoft Excel[M]. Berkeley C A:Berkeley Geo-chronology Centre, Special Publication 2008, 4:1-76.
Rowley D B, Xue F, Tucker R D, et al. Ages of ultrahigh pressure metamorphic and source orthognisses from the eastern Dabie Shan:U/Th zircon geochronology[J]. Earth and Planetary Science Let-ters. 1997, 151:191-203. doi: 10.1016/S0012-821X(97)81848-1
Crofu F, Hanchar J M, Hoskin P W, et al. Atlas of zircon textures[J]. Reviews in Mineralogy and Geochemistry, 2003, 53:469-495. doi: 10.2113/0530469
吴元保, 郑永飞.锆石成因矿物学研究及其对U-Pb年龄解释的制约[J].科学通报, 2004, 49(16):1589-1604. doi: 10.3321/j.issn:0023-074X.2004.16.002 吴元保, 郑永飞, 龚冰, 等.北淮阳新开岭地区花岗岩锆石U-Pb年龄和氧同位素组成[J].地球科学-中国地质大学学报, 2005, 30(6):659-672. http://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200506002.htm Rowley D B. Age of initiation of collision between India and Asia:A review of stratigraphic data[J]. Earth and Planetary Science Let-ters, 1996, 145(1/4):1-13. http://www.academia.edu/10329059/Age_of_initiation_of_collision_between_India_and_Asia_A_review_of_stratigraphic_data
莫宣学, 赵志丹, 邓晋福, 等.印度-亚洲大陆碰撞过程的火山响应[J].地学前缘, 2003, 10(3):135-148. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200304001016.htm 王成善, 李祥辉, 胡修棉.再论印度-亚洲大陆碰撞的启动时间[J].地质学报, 2003, 77(1):16-24. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200301005.htm 侯增谦, 杨竹森, 徐文艺, 等.青藏高原碰撞造山带:Ⅰ.主碰撞造山成矿作用[J].矿床地质, 2006, 25(4):337-358. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200604000.htm 侯增谦, 潘桂棠, 王安建, 等.青藏高原碰撞造山带:Ⅱ.晚碰撞转换成矿作用[J].矿床地质, 2006, 25(5):521-533. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200605000.htm 杨勇, 罗泰义, 黄智龙, 等.西藏纳如松多银铅矿S、Pb同位素组成:对成矿物质来源的指示[J].矿物学报, 2010b, 30(3):311-318. http://www.cnki.com.cn/Article/CJFDTOTAL-KWXB201003007.htm 中国地质调查局发展研究中心.国内外铅锌矿床成矿理论与找矿方法[M].北京:地质出版社, 2004. 习中军. 西藏自治区谢通门县德新矿区铅多金属矿普查2011年地质总结. 西藏自治区地质矿产勘查开发局第五地质大队, 2011.