四川乌斯河铅锌矿床成矿物质来源及矿床成因：来自原位S-Pb同位素证据*

韦晨 叶霖 李珍立 胡宇思 黄智龙 刘玉平 王皓宇

1. 中国科学院地球化学研究所,矿床地球化学国家重点实验室，贵阳 5500812. 中国科学院大学，北京 1000491.

扬子板块西南缘的四川、云南和贵州接壤地区铅锌矿床(点)成群成带分布，目前已发现不同规模的铅锌矿床超过400余处(柳贺昌和林文达，1999)，被称为川滇黔铅锌矿矿集区(图1；王奖臻等，2002；黄智龙等，2004；张长青，2005；Yeetal.，2011；Zhangetal.，2015)，累计探明的铅锌矿资源量超过2000万吨，约占全国总量的27%，是我国重要的Pb、Zn、Ag及多种分散元素的生产基地之一(黄智龙等，2011；Zhangetal.，2015)。这些矿床明显受构造控制，主要赋存于震旦系-下二叠统碳酸盐岩地层中，其中震旦-寒武系地层为该区最为重要的铅锌矿赋矿层位，产出多个大型铅锌矿床如四川大梁子、天宝山、赤普、云南金沙厂、茂租等。四川乌斯河铅锌矿床是在川滇黔接壤地区近年新发现的大型铅锌矿床之一，但与赋存于同一层位的其他铅锌矿床相比，该矿床的研究程度相对较低，其成矿物质来源不清，仅少数学者根据野外地质特征以及少量的地球化学数据提出不同的矿床成因类型，包括：(1)热水沉积成因(林方成，2005；Zhuetal.，2018)；(2)喷流沉积-后期热液叠加富集改造成因(郑绪忠，2012)；(3)MVT型铅锌矿(Xiongetal.，2018；Luoetal.，2020)。事实上，各种观点仍然缺少地球化学依据的支持，致使该矿床铅锌成矿作用的认识依然存在争议。对乌斯河大型铅锌矿开展深入的研究及比较矿床学分析，不仅有助于认识其成矿物质来源，建立碳酸盐岩为赋矿围岩的铅锌成矿模型，也将对该地区及外围地质勘探的深入具有重要的理论指导意义。

图1 川滇黔接壤地区大地构造位置(a，据Ye et al.，2011修改)和铅锌矿床分布图(b，据柳贺昌和林文达，1999修改)Fig.1 Tectonic setting of the Sichuan-Yunnan-Guizhou (SYG) Pb-Zn metallogenic province (a, modified after Ye et al.，2011) and Pb-Zn deposits distribution and regional structures in this area (b, modified after Liu and Lin，1999)

硫和铅同位素是示踪硫化物矿床成矿流体和物质来源的重要手段之一，已被广泛应用于各类硫化物矿床的成因和形成过程研究(李文博等，2006；Zhouetal.，2013，2015；Baoetal.，2016；金中国等，2016；Lietal.，2020)。长期以来，硫、铅等稳定同位素分析多采用整体化学溶样法后利用质谱对样品整体分析测定，其获得结果虽精度较高，但很难将不同期次或不同环带的硫化物完全区分开，获得的分析结果为整个颗粒不同部分的混合值，不能真正代表成矿物质或成矿流体系统中真实的同位素组成，从而掩盖了成矿过程中的一些重要信息(Chenetal.，2014；Yuanetal.，2015；Baoetal.，2017)。为此，本文借助激光剥蚀多接收器等离子体质谱(LA-MC-ICPMS)，对采自乌斯河铅锌矿床不同成矿阶段的硫化物(黄铁矿、闪锌矿和方铅矿)，开展原位硫和铅同位素地球化学研究，旨在获取精确同位素数据，从而查明该矿床的成矿物质来源、还原硫的形成机制以及示踪成矿过程，为深入认识该类矿床铅锌成矿作用提供地质地球化学依据。

1 矿区及矿床地质特征

1.1 区域地质

川滇黔铅锌矿集区位于扬子地台西南缘(图1a)，北侧与秦岭褶皱带、松潘-甘孜缝合带相接，西南侧紧邻三江造山带，东南临近华南褶皱带和右江盆地(柳贺昌和林文达，1999)。研究区具有基底和盖层双层结构特征(柳贺昌和林文达，1999)，其中基底由古太古代-古元古代的中-深变质杂岩、中元古代的变质细碎屑岩夹变质火山沉积岩以及新元古代浅变质碎屑岩和碳酸盐岩组成。沉积盖层地层出露较为齐全，除白垩系、第三系缺失外，自震旦系至第四系地层均有分布(Yanetal.，2003)，其中石炭系、二叠系和三叠系地层出露齐全、分布广且沉积厚度大。

区内构造极为复杂，主体构造呈NS向、NW向和NEE-NNE向展布，主要包括NS向安宁河断裂带、小江断裂带，NW向昭通-曲靖隐伏断裂带和垭都-紫云断裂带及NE师宗-水城断裂(图1b；柳贺昌和林文达，1999；张志斌，2006)。这些深大断裂带不仅控制着区内岩相分布，特别是赋矿的碳酸盐岩地层，同时控制了铅锌矿床(点)的展布(黄智龙等，2004；Zhangetal.，2015)。迄今为止，区内已发现铅锌矿床(点)400余处，这些矿床绝大多数分布于上三叠统-震旦系碳酸盐岩中，主要产于下二叠统阳新组-震旦系灯影组白云岩中，矿体明显受控于断裂构造，具有成矿温度较高(170～270℃)、中-低盐度(4.2%～13.5%NaCleqv)、成矿明显晚于围岩形成时代等特征(黄智龙等，2004；张长青，2005，2008；Hu and Zhou，2012；Zhangetal.，2015)。

1.2 矿区地质

近年来，在康滇地轴北段的四川汉源县乌斯河一带，发现多个(似)层状铅锌矿床(点)，其中包括牛心山、宝水溪、白塔、中溪坪、养善坪、核桃坪和乌斯河等铅锌矿床(点)，这些铅锌矿床(点)均稳定产于上震旦统-下寒武统的灯影组顶部麦地坪段含磷硅质白云岩中，平面上构成了长度大于50km的铅锌矿化带(林方成，2005；郑绪忠，2012)。乌斯河铅锌矿床位于川滇黔铅锌矿集区北部(图1b)，是该地区最具代表性大型铅锌矿床之一(林方成，2005)，据四川省地质矿产局207地质队在2003年的矿区普查估算结果，该矿床的铅锌资源量为370万吨，Pb+Zn平均品位～15.7% (四川省地质矿产局207地质队,2003(1)据四川省地质矿产局207地质队.2003.四川省汉源县乌斯河铅锌矿床普查报告)。

矿区出露的地层主要包括中元古界峨边群变质岩、南华系苏雄组陆相火山岩和碎屑岩、震旦统灯影组白云岩、寒武系碳酸盐岩和碎屑岩、奥陶系-二叠系碳酸盐岩和碎屑岩以及第四系沉积物(图2)。其中以上震旦统灯影组和寒武系地层沉积厚度最大(最大可达上千米)是该地区主要的铅锌矿化层位。自下而上震旦系灯影组可划分为三个岩性段：(1)第一段(Z3dn1)，以局部含燧石微晶白云岩为主，厚约100～200m；(2)第二段(Z3dn2)，主要岩性为硅质条带白云岩其底部发育一层厚约10～30m的紫色砂页岩；(3)第三段(Z3dn3)，即麦地坪段，岩性主要为含小壳动物化石及少量胶磷矿的硅质白云岩，厚约30～70m，其中，麦地坪段顶部的角砾状白云岩以及黑色硅质岩为该矿床主要的赋矿围岩。

矿区内构造极不发育，主要控矿构造为王帽山逆断层和万里村向斜。其中，王帽山逆断层位于矿区北西部，属近NE走向，倾向NW，倾角65°～70°，断层破碎带宽约15～30m。该断层上下盘分别为震旦系灯影组和寒武-奥陶系地层。万里村向斜，轴向近NS，核部依次出露奥陶系、志留系和二叠系地层，两翼分为寒武系及震旦系地层，表现出中间新两边老的向斜构造特征，其东翼岩层倾角为5°～16°，而西翼岩层倾角为4°～22°，该向斜属于轴部宽缓而两翼不对称的向斜，南端倾伏于乌斯河一带，倾伏角～10°，北端止于王帽山断层。此外，在矿区东部还发育一些小型断层，这些次级断层仅在矿山坑道中发育，其规模小且断距小，仅对局部铅锌矿体造成破坏。

图2 四川乌斯河矿区地质图(a)及19-19′剖面图(b)(据林方成，2005)Fig.2 Geological map (a) and 19-19′ section (b) of the Wusihe Pb-Zn deposit, Sichuan Province (modified after Lin，2005)

矿区铅锌矿体的空间展布受万里村向斜控制，矿体产状与地层产状基本一致，主要呈层状、似层状或透镜状顺层产出(图3a, b)。矿体横向延伸规模较大，且连续性较好，长度数米至数百米不等，地表露头连续可达6000m以上(林方成，2005；郑绪忠，2012)。前人依据铅锌矿化程度和矿物岩石组合，将麦地坪段矿化带分为上下两个矿化带：(1)上矿化带厚度变化不大，以层状的矿体稳定产出(平均厚度～13m)，赋矿岩性主要为中厚层状的深灰色白云岩及泥炭质白云岩，该矿化带内铅锌矿体为乌斯河矿床的主要工业矿体，矿体形态较规则，局部可见脉状矿体沿穿层裂隙充填；(2)下矿化带厚度变化较大，平均厚度为15m，是矿区次要赋矿层位，岩性为中厚层状的白云岩及硅质白云岩，其铅锌矿化较弱，表现为局部发育的透镜状矿体，围岩蚀变强烈，主要发育硅化和碳酸盐化。矿区内铅锌矿体多为单层，厚度变化为0.2～3m(平均1m)，矿石品位分别为Pb 0.025%～14.12%(平均3.15%)和Zn 0.59%～37.49%(平均12.50%)。在层间破碎带或者构造裂隙的交汇部位，矿体变厚，且品位明显变高。乌斯河铅锌矿床围岩蚀变较弱，仅在个别构造发育地段出现较小规模的围岩蚀变，并且主要发育弱硅化、白云石化、黄铁矿化、沥青化或萤石化等(林方成，2005；郑绪忠，2012)。

矿床中矿石矿物组合较为简单，主要金属矿物为闪锌矿和方铅矿，次为黄铁矿，偶见硫锑铅矿(图3、图4)。闪锌矿颜色较浅，以浅棕色为主，偶见红棕色(图3d，h)，脉石矿物主要为石英、白云石、沥青和方解石等(图3、图4)。矿石构造主要为块状(图3c，d)、条带状(图3e)、角砾状(图3f)、脉状(图3g)和浸染状(图3h)。主要结构包括自形、半自形粒状(图4a, d, f, i)，其次为交代残余结构(图4b)和交代溶蚀结构(图4c，g)和碎裂结构(图4h)。

图3 乌斯河铅锌矿床矿体产出特征及矿物组合 (a)白云岩中层状铅锌矿细脉；(b)层状铅锌矿体； (c)块状铅锌矿石；(d)硫化物与沥青紧密共生；(e)条带状矿石；(f)白云岩角砾被硫化物胶结；(g)脉状铅锌矿石；(h)浸染状铅锌矿石. 矿物缩写: Sp-闪锌矿；Gn-方铅矿；Py-黄铁矿；Bit-沥青；HD-热液白云石；Cal-方解石；Q-石英；Dol-白云岩Fig.3 Photos of orebodies and mineral assemblage in the Wusihe Pb-Zn deposit

在镜下可见到自形的黄铁矿颗粒(Py1)空隙被稍晚形成的闪锌矿(Sp1)和白云石充填(图4a)；被主阶段黄铁矿(Py2)包裹(图4c)，表明自形的黄铁矿为最早期形成的硫化物。闪锌矿为该矿床最为主要的硫化物，其形成可分为两个阶段：早阶段生成的闪锌矿(Sp1)交代方铅矿(Gn1)(图4b)或充填在Py1的边缘(图4a，c)；主阶段形成的闪锌矿(Sp2)充填在沥青的空隙(图4f)或被最晚形成的方铅矿(Gn2)交代(图4g)。此外，镜下可观察到大量的胶状闪锌矿产出(郑绪忠，2012)。早期形成的方铅矿(Gn1)被闪锌矿交代呈残余小颗粒状分布在Sp1中；晚生成的方铅矿交代早期形成的硫化物(图4c)或被沥青包括(图4e)。此外，石英为最早生成的脉石矿物，多孔状Sp1分布在其两侧，后被又被方解石脉溶蚀交代(图4d)，表明方解石脉最晚形成。

根据矿物在野外和镜下鉴定特征，初步拟定该矿床热液期可分为三个成矿阶段：早阶段、主阶段和晚阶段。具体的矿物生成顺序见图5。

2 分析方法

2.1 S同位素

硫化物原位S同位素测试在西北大学大陆动力学国家重点实验室完成。其中激光剥蚀系统为193nm 准分子激光剥蚀系统 (RESOlution M-50, ASI)，多接收等离子体质谱为为ASI 公司的Nu Plasma 1700 MC-ICP-MS。S同位素分析采用单点剥蚀方式，剥蚀斑径为25μm，激光频率为3Hz，能量密度为3.6J/cm2。测试过程中对不同类型硫化物一般选择相同基体作为标准样品，其中方铅矿标样采用CBI-3，黄铁矿标样Py-4、闪锌矿标样NBS123、最后相对于内标的δ34Sv-standard值全部转化为CDT值(δ34SV-CDT)，测试精度优于 ±0.2‰。详细分析过程参照Chenetal.(2014)和Baoetal. (2017)。

沥青的S同位素测试在中国科学院地球化学研究所环境地球化学国家重点实验室完成，采用英国GV公司生产的GC-IsoPrime连续流气体同位素质谱仪。根据沥青中S的含量(>1%)称取相应的样品量(使测量样品中S含量达到50μg)，并在CF-IsoPrime JB144连续流气体同位素质谱仪上测得硫同位素组成。本方法用国标GBW04414、GBW04415和IAEA-S3做内标，以CDT为标准，分析精度为0.2‰。

图4 乌斯河铅锌矿床硫化物矿石显微结构特征 (a) 晚生成闪锌矿(Sp1)充填的在自形黄铁矿(Py1)的孔隙中；(b)细粒的方铅矿(Gn1)分布在稍晚生成的闪锌矿(Sp1)中；后被更晚生成方铅矿(Gn2)交代；(c)早生成自形黄铁矿(Py1)和闪锌矿(Sp1)被晚生成黄铁矿(Py2)和方铅矿(Gn2)交代；(d)晚期方解石脉交代早期石英脉；(e)沥青分布在Sp1边缘；(f)沥青的空隙被晚生成Sp2充填；(g) 主阶段生成的闪锌矿(Sp2)被晚生成方铅矿(Gn2)交代；(h)方解石脉切穿晚生成闪锌矿(Sp2)晶体，形成碎裂结构；(i)主阶段生成方铅矿(Gn2)和硫锑铅矿分布在自形的方解石孔隙中，后被沥青充填. 矿物缩写: Jam-硫锑铅矿Fig.4 The micro-texture features of sulfide ores in the Wusihe Pb-Zn deposit

2.2 原位Pb同位素

硫化物原位Pb同位素分析在西北大学大陆动力学国家重点实验室完成。实验流程：将样品磨制成厚约80μm的探针片，原位铅同位素分析采用LA-MC-ICPMS和膜去溶装置联机进行。其中激光剥蚀系统为193nm RESOLoution M-50，多接收等离子体质谱为Nu Plasma II型MC-ICP-MS(Nu Ins，UK)，膜去溶系统是Cetac公司的Aridus II型 Desolvation Nebulizer System (Aridus II，Cetac，USA)。剥蚀时间约50s，剥蚀方铅矿的激光频率为3Hz，剥蚀斑束为9～100μm不等。仪器分析所得数据经质量监控标样NIST SRM 981归一化校正后获得最终Pb同位素比值。国际标准物质NIST SRM 610作为外标，其外部精度优于±0.01%。详细分析过程参见Chenetal.(2014)和Yuanetal.(2015)。

3 分析结果

3.1 硫同位素

本次测试了8件硫化物和9件沥青样品，共计45个硫同位素数据(表1)，可以看出不同矿物的δ34S值变化范围较宽(图6a)，介于+1.3‰～+24.3‰，均值+14.6‰(n=45)；其中9件沥青样品的δ34S值介于+14.8‰～+24.3‰，均值+20.0‰；8件硫化物中36个原位测点的δ34S值介于 +1.3‰～+23.3‰，均值+13.3‰；比前人通过溶液法获得S同位素数据范围更宽(+1.8‰～+15.5‰；林方成，2005；李同柱，2007；Xiongetal.，2018；Luoetal.，2020)。同时，不同成矿阶段的硫化物硫同位素明显不同(图6b)。成矿早阶段的Py1、Sp1和Gn1的δ34S值分别为+1.3‰～+9.8‰(均值+4.7‰；n=4)、+8.5‰～+14.2‰(均值+12.2‰；n=8)和+7.8‰～+10.9‰(均值9.8‰；n=3)；而主成矿阶段的硫化物硫同位素组成明显更高，Py2、Sp2和Gn2的δ34S值分别为+15.4‰～+23.3‰(均值20.6；n=5)、+13.0‰～+19.3‰(均值+15.4‰；n=10)和+11.6‰～+14.3‰(均值+12.5‰；n=6)。

3.2 Pb同位素

本次共测试9件硫化物样品 (闪锌矿和方铅矿)19个测点的原位铅同位素组成(表2、图7)。全部硫化物测点的206Pb/204Pb=18.079～18.140，均值18.108；207Pb/204Pb= 15.639～15.689，均值15.662；208Pb/204Pb=38.201～38.343，均值38.262；所有测点的变化范围明显比前人获得Pb同位素组成更为集中(Xiongetal.，2018；Luoetal.，2020)。但不同成矿阶段的Pb同位素组成存在明显差异，主成矿阶段的硫化物的206Pb/204Pb、207Pb/204Pb和208Pb/204Pb值(均值分别为18.114、15.664和38.267)略高于成矿早阶段硫化物样品的206Pb/204Pb、207Pb/204Pb和208Pb/204Pb值(均值分别为18.099、15.658和38.254)，暗示前者更加富集放射性成因Pb。

表1 乌斯河铅锌矿床硫化物LA-MC-ICPMS原位S同位素组成

图5 乌斯河铅锌矿床矿物生成顺序Fig.5 Mineral paragenesis of the Wusihe Pb-Zn deposit, Sichuan Province

4 讨论

4.1 硫源及硫酸盐还原机制

矿物学研究表明，乌斯河铅锌矿金属矿物组合简单且仅为硫化物，以闪锌矿为主，含少量方铅矿、其次为黄铁矿，未见硫酸盐类矿物(图3、图4)，表明硫化物沉淀时成矿流体的氧逸度相对较低，同时闪锌矿中发育大量富烃类的包裹体，进一步说明成矿热液流体中硫主要以H2S形成存在。在这种情况下，金属硫化物的硫同位素组成(δ34S硫化物)可以近似代表成矿流体中的总硫同位素组成(δ34S∑s)即δ34S硫化物≈δ34S∑s(Ohmoto，1972；Ohmoto and Goldhaber，1997；Seal II，2006)。因此，本次测试的硫化物的硫同位素平均值(+13.3‰)近似代表成矿流体中δ34S∑s的值。该矿床硫同位素值明显高于幔源岩浆硫的δ34S值(-3‰～+3‰；Chaussidonetal.，1989)和变质基底的岩石中硫同位素组成(～5‰，Zhangetal., 2019)，故排除了岩浆和基底岩石提供大量硫的可能性。值得注意的是，本区铅锌矿的赋矿地层及其上覆寒武系地层为一套蒸发相白云岩，含有石膏、重晶石等硫酸盐矿物。前人获得其δ34S值在28.0‰～34.8‰之间(白俊豪等，2013)，与震旦-寒武纪海水硫酸盐的δ34S值相近(28‰～34‰；Claypooletal.，1980)。因此，乌斯河铅锌矿床成矿流体中还原硫可能主要来源于赋矿地层的蒸发岩，这与邻区的赋存于石炭系碳酸盐岩中会泽、毛坪和富乐铅锌矿床基本一致(黄智龙等，2004；李文博等，2006；Lietal.，2020)此外，大量研究表明，硫酸盐(SO42-)还原产生还原硫(S2-)主要通过两种途径，细菌硫酸盐岩还原(BSR)和热化学硫酸盐还原(TSR)，但两种硫酸盐还原机制的触发的温度以及产生的硫同位素分馏明显不同(Seal II，2006；Basukietal.，2008)。通常来讲，BSR发生温度较低(适宜温度30～40℃；最高<110℃)，形成的S2-的硫同位素值多数为负值，且变化范围较宽，可达55‰(Jørgensenetal.，1992；Seal II，2006)。TSR则发生的温度相对较高(>120℃)，形成的S2-的硫同位素值相对稳定，与硫酸盐的分馏<15‰(Ohmoto，1972；Jørgensenetal.，1992)。乌斯河铅锌矿成矿温度介于120～260℃(熊索菲等，2016)，明显高于细菌存活的温度，且乌斯河铅锌矿硫化物的明显富集重的硫同位素，多数δ34S值>10‰，同时矿区发育大量的沥青，其δ34S值14.8‰～24.3‰，这些均表明TSR作用在乌斯河铅锌矿成矿流体的还原硫生成过程起到至关重要的作用。但本次研究发现乌斯河铅锌矿除了富集重的硫同位素(图6)外，成矿早阶段的硫化物δ34S值明显较低(<10%；图6)，比震旦系海水硫酸盐低14‰～32‰，不难发现单一TSR作用无法解释这一现象，其可能存在两种解释：(1)沉积成岩过程中地层圈闭一些BSR作用产生的H2S；(2)成矿早阶段有外部硫加入。考虑到BSR作用生成的δ34S值变化范围较宽，且多为负值(Ohmoto，1972；Seal II，2006)。因此，本文认为第二种解释更为合理。事实上，乌斯河矿床周缘发育了许多油气田，如四川盆地西南缘产于震旦系灯影组的威远大型气田，天然气中含0.8%～1.4%的H2S，其δ34S值多集中在11.5‰～14.4‰，同时四川盆地一些气田中沥青的硫化物δ34S值在 +5.80‰～+8.45‰之间(朱光有等，2006)。此外，矿床成矿早阶段发育了大量胶状闪锌矿(图8；郑绪忠，2012)，这种闪锌矿是由于充足的还原硫与携带金属的流体快速混合的产物，受限于还原硫生成的速率，BSR作用一般不会形成这类闪锌矿。因此，成矿早阶段的还原硫可能来源于古老油气藏圈闭的H2S。单个闪锌矿晶体形态和δ34S值变化特征也印证了这一点(图8)，如其核部为胶状而边缘为粗晶状，同时边缘的δ34S值明显高于核部，呈现先增加后降低的趋势。值得注意的是，边部呈现降低的硫同位素可能为快速硫化物沉淀产生强烈的瑞利分馏所致，全面深入理解其的分馏机制还需进一步研究。

表2 乌斯河铅锌矿床硫化物LA-MC-ICPMS原位Pb同位素组成

图6 乌斯河铅锌矿床不同矿物硫同位素直方图(a)和不同阶段硫化物硫同位素组成散点图(b) 湿化学法获得硫同位素数据引自林方成(2005)，李同柱(2007)，Xiong et al. (2018)和Luo et al.(2020).Py1-早阶段黄铁矿；Sp1-早阶段闪锌矿；Gn1-早阶段方铅矿；Bt-沥青；Py2-主阶段黄铁矿；Sp2-主阶段闪锌矿；Gn2-主阶段方铅矿Fig.6 Histogram of S isotope compositions of S-bearing minerals (a) and diagram of S-bearing minerals at different stages from the Wusihe Pb-Zn deposit (b)

图7 乌斯河铅锌矿床及围岩硫化物原位铅同位素组成图解 (a)乌斯河铅锌矿床硫化物与经年龄校正的赋矿围岩和基底地层的207Pb/204Pb-206Pb/204Pb对比图解 (数据来源于林方成，1995；常向阳，1997；白俊豪，2013);不同成矿阶段硫化物的 207Pb/204Pb-206Pb/204Pb(b)和208Pb/204Pb-206Pb/204Pb (c)Fig.7 The comparison plots of in situ S isotopic ratios of sulfides at different stages from the Wusihe Pb-Zn deposit

图8 单个闪锌矿晶体的透射光(a)和反射光(b)照片以及原位硫同位素组成变化趋势(c)Fig.8 The photographs of transmission light (a) and reflected light (b) in a single sphalerite crystal, and corresponding sulfur isotope changes from core to rim (c)

综上所述，还原硫主要来源于赋矿地层的蒸发岩，TSR作用在硫酸盐还原过程中起到举足轻重的作用，成矿早阶段可能有油气系统中先存的 H2S 加入。

4.2 成矿物质来源

以往只对乌斯河铅锌矿床的方铅矿进行了Pb同位素地球化学研究，闪锌矿中Pb同位素研究尚未报道(Xiongetal.，2018；Luoetal.，2020)。本次分析了方铅矿和闪锌矿的Pb同位素组成，旨在示踪铅锌等成矿物质来源。目前对川滇黔接壤地区铅锌矿床成矿物质来源研究基本上达成共识，潜在的成矿物质源区主要包括变质基底、赋矿围岩以及峨眉山玄武岩(柳贺昌和林文达，1999；黄智龙等，2004；李文博等，2006；Zhouetal.，2013，2015；Baoetal.，2016；Xuetal.，2020)。由于硫化物中U和Th含量较低，其Pb同位素组成无需进行年龄校正即可代表成矿流体的Pb同位素组成(Carretal., 1995)。本次研究发现乌斯河铅锌矿不同硫化物(闪锌矿和方铅矿)的原位Pb同位素变化范围较窄，表明成矿金属的来源较为单一或混合多个Pb同位素组成相似的源区(图7a)。考虑到上述三类岩石中普通Pb含量较高，在示踪成矿物质来源时，需要对其Pb同位素组成进行校正。近来，有学者报道乌斯河铅锌矿床Rb-Sr等时线年龄为411±10Ma(MSWD=1.4；Xiongetal.，2018)，其成矿年龄明显早于峨眉山玄武岩喷发(～260Ma；Heetal.，2007)，峨眉山玄武岩作为潜在的成矿物质源区基本可以排除。因此，本文采用411Ma作为该矿床的成矿年龄对变质基底、赋矿围岩石进行Pb同位素校正。

将乌斯河铅锌矿床硫化物原位Pb同位素组成投入经成矿年龄校正的变质基底和赋矿围岩图解上(图7a)，不难发现乌斯河铅锌矿床的硫化物位于变质基底范围内，且更靠近赋矿围岩(图7a)。因而，乌斯河铅锌矿成矿物质主要来源于变质基底，赋矿围岩可能贡献少量的成矿物质。值得注意的是，不同成矿阶段的硫化物的Pb同位素明显不同，主成矿阶段的硫化物更加富集放射性成因Pb，与赋矿围岩的Pb同位素组成范围较为相似(图7b，c)，暗示主成矿阶段赋矿围岩贡献成矿物质的比例更高，但考虑到碳酸盐岩地层中Pb、Zn含量较低，分别<6.1×10-6和<3.1×10-6(李同柱，2007)，即使发生强烈水岩反应，围岩对成矿金属的贡献也相对有限，这就导致观察到Pb同位素线性混合趋势并不明显，Zn同位素组成也证实这一点(Zhangetal.，2019)。可见，该矿床金属成矿物质主要来源来于变质基底，水岩反应可能使赋矿围岩贡献少量的成矿物质。

4.3 矿床成因

乌斯河铅锌矿硫化物矿体与地层整合产出，矿层延伸规模巨大，条带状、层纹状构造发育，致使部分学者将其归于热水喷流沉积成因(林方成，2005；郑绪忠，2012；Zhuetal.，2018)。事实上矿区中层状、似层状和陡倾斜的脉状矿体均有发育，显示明显的后生成矿特征(郑绪忠，2012)，因此，可以排除原生沉积成矿的可能性。

比较矿床学分析有助于我们进一步认识乌斯河铅锌矿床的类型。表3给出乌斯河铅锌矿与国内外典型的MVT铅锌矿的详细对比。就地质产出特征而言，本矿床中矿物组成相对简单，主要包括闪锌矿、方铅矿、白云石、沥青和方解石，含少量的黄铁矿和石英；矿石类型以浸染状、脉状、角砾状和块状为主；围岩蚀变较弱，主要包括方解石化、白云石化，偶见黄铁矿化；这些地质特征与典型的MVT型铅锌矿床(Leach and Sangster，1993；Leachetal.，2001)基本一致。此外，流体包裹体分析结果表明，其成矿流体以中-低温(130～260℃)、中高盐度(5.4%～21.8%NaCleqv)为特征(熊索菲等，2016)，也与MVT型铅锌矿床成矿温度(90～200℃)和盐度(10%～30 %NaCleqv)接近(Leachetal.，2001，2005)。近来，李珍立(2019)对乌斯河铅锌矿的闪锌矿开展了详细的LA-ICPMS研究，其结果表明闪锌矿中In、Co和Mn含量较低，明显不同于块状硫化物和岩浆热液矿床(In>50×10-6，Mn>1000×10-6；Cooketal.，2009；Yeetal.，2011；Weietal.，2018a)，也不同于远源矽卡岩型矿床(Co>200×10-6，Mn>1000×10-6；Yeetal.，2011)，矿床中闪锌矿相对富集Cd、Ge和Ga，这与典型的MVT型铅锌矿闪锌矿微量元素特征相似(Ge>45×10-6，Cd>2700×10-6；Yeetal.，2011；叶霖等，2016；Weietal.，2018b；胡宇思等，2019)。此外，本次研究发现乌斯河铅锌矿床硫化物的还原硫主要来源于赋矿地层的蒸发岩，铅锌等成矿物质主要来源来于变质基底。可见，上述地质地球化学特征均显示乌斯河铅锌矿床与典型MVT型矿床较为相似，故本文认为本矿床属于MVT型成因类型。

4.4 成矿过程

综合各类地质地化资料, 不难发现乌斯河铅锌矿床成矿金属主要来源于变质基底，而还原硫主要来自于赋矿地层的蒸发岩的热化学还原作用。此外，该矿床成矿流体温度较低，介于130～260℃之间(熊索菲等，2016)，已有的实验地球化学研究表明温度较低时，成矿流体很难同时携带大量的Pb、Zn金属离子和还原硫(Zhongetal.，2015)。因此，本文倾向于认为来自基底富Pb、Zn等成矿元素的流体与赋存于地层中富含还原硫的流体混合是乌斯河铅锌矿床金属硫化物沉淀的主要机制。该矿床的形成可能经历了如下过程：晚志留世(～410Ma)的华北板块与华南板块发生碰撞造山致使盆地卤水发生大规模迁移，淋滤和萃取基底岩石中Pb、Zn等成矿元素，形成富金属的成矿流体，沿区域性深大断裂向上迁移，在震旦系灯影组白云岩地层中成矿流体与少量先存的H2S快速结合到达过饱和状态，沉淀出胶状闪锌矿，然后TSR作用启动，还原地层中蒸发岩形成大量的还原硫，与富Pb、Zn 等成矿元素的流体混合引起大量的铅锌硫化物沉淀并释放出H+。生成的H+与围岩发生水/岩相互反应，进行元素和同位素交换，同时溶解碳酸盐岩产生大量CO2调节成矿流体中的物理化学条件 (如pH、 Eh)， 碳酸盐矿物反复溶解-再结晶过程维持成矿流体系统稳定(Warren，2000)，促进铅锌等硫化物的持续沉淀。这一过程也很好的解释了围岩可能贡献了少量的成矿物质。

表3 乌斯河铅锌矿床与国内外典型MVT铅锌矿床地质、地球化学特征对比

5 结论

通过乌斯河铅锌矿床地质特征及硫化物原位S和Pb同位素对比研究，本文获得以下认识：

(1)乌斯河铅锌矿床的还原硫主要来源于赋矿地层中的蒸发岩，TSR作用在硫酸盐还原过程中起到举足轻重的作用，成矿早阶段可能有油气系统中先存的 H2S 加入；

(2)乌斯河铅锌矿床的成矿物质主要来源来于变质基底，水岩反应致使围岩贡献少量的成矿物质；

(3)该矿床属于MVT型铅锌矿床，其铅锌成矿作用可能与加里东晚期碰撞造山运动引发的大规模流体运移有关。流体混合作用是乌斯河矿床铅锌等成矿物质主要富集沉淀机制。

致谢实验工作得到得到西北大学大陆动力学国家重点实验室包志安和陈开运工程师等的帮助；成文过程与中国科学院地球化学研究所罗泰义研究员、刘仕玉、吴涛和向震中等进行了有益的讨论；审稿专家和本刊编辑提出了诸多富有建设性的修改意见和建议；在此对他们表示衷心的感谢！