2. 云南冶金集团股份有限公司, 昆明 650224;
3. 云南省国土资源规划设计研究院, 昆明 650216
2. Yunnan Metallurgical Group Co., Ltd., Kunming 650224, China;
3. Yunnan Land Resources Planning and Design Institute, Kunming 650216, China
0 前言
老君山成矿区位于云南省的东南部,以盛产锡、铅、锌、钨、银著称。在1 380 km2的范围内已探明都龙超大型锡矿床1处,新寨(锡)、南秧田(钨)大型矿床2处,戈岭(锡)、老寨(锡)中型矿床2处,以及小型多金属矿床(点)20余处,是滇东南多金属成矿带[1, 2](个旧——薄竹山——老君山)的重要组成部分。都龙锡多金属矿床曾被认为是一个典型的岩浆热液矿床[3, 4, 5]。随着研究的进行,越来越多的学者认为老君山地区乃至滇东南地区超大型、大型块状硫化物矿床的成因模式,不能简单地归结为花岗岩成因;在老君山成矿区围绕都龙、新寨、南秧田矿床先后出现了沉积-变质-岩浆热液叠加改造成因[6]、沉积-变质热液改造成因[7]、海底(火山)喷流热水沉积成因[8, 9, 10, 11, 12, 13]等。纵观都龙锡矿床以往研究成果,该矿床应为多期复合叠加成矿模式,与成矿相关的地质作用主要包括早加里东期(火山)喷流沉积成矿作用、印支期区域变质成矿作用和燕山期花岗岩叠加改造成矿作用。
老君山成矿区经历了多期成矿地质作用,不同地质作用对成矿的影响如何?形成的各类矿石、矿物在同位素组成上有何区别与联系?将是本文研究的重点。作者将运用同位素地球化学原理及研究方法,分析不同矿化类型、不同矿物同位素组成特征,并结合已有研究成果及矿区地质情况探讨上述问题,以深化该区成矿理论,指导找矿实践。
1 区域地质背景
老君山成矿区大地构造处于扬子地块、江南造山带、三江造山带等几大构造单元接合部位(图 1)。研究区主要出露寒武纪、泥盆纪、三叠纪地层。寒武系为扬子地块前缘海盆环境沉积的碎屑岩、泥质岩、碳酸盐岩,岩石普遍变质为片岩类、变粒岩类、大理岩类和矽卡岩类等,下、中寒武统为主要赋矿层位。泥盆系为内陆湖沼-浅海相碎屑岩、泥质岩和碳酸盐岩,部分层位具千枚岩化。三叠系为千枚岩化板岩、细砂岩。研究区中部为老君山燕山期S型复式花岗岩体,出露面积约153 km2。紧邻老君山花岗岩体形成了一套穹窿状变形变质岩系,被称为“老君山变质核杂岩”[14, 15, 16, 17],变形-变质作用的峰期为印支期(约235 Ma)[18, 19]。成矿区受到多期地质作用影响(早加里东期、印支期、燕山期),成矿过程十分复杂。 2 岩石类型与产状 老君山成矿区主要变质岩类型包括:①大理岩类,有方解石大理岩、白云石大理岩、含炭质大理岩、含砾钙泥质大理岩,与其他变质岩呈互层状产出。②片岩类,主要有绿泥石片岩、角闪石片岩、云母片岩、云母石英片岩、石英电气石片岩、长石云母片岩等,在研究区普遍发育。③混合片麻岩类,包括眼球状花岗片麻岩、条痕状花岗片麻岩,在南秧田地区发育较广。④矽卡岩是成矿区内主要的赋矿岩石,根据矿物组合及复杂程度,可分为简单矽卡岩及复杂矽卡岩两类。简单矽卡岩矿物组合简单,一般不超过4种非金属矿物,主要矿物为透辉石、绿帘石、阳起石、普通角闪石,含少许磁铁矿、磁黄铁矿、黄铁矿等,层位稳定,呈似层状、层状分布,厚度小(一般小于40 m),延长小于40 km。复杂矽卡岩由10余种硅酸盐矿物组成,主要有透辉石、钙铁辉石、斜黝帘石、绿帘石、绿泥石、阳起石、石榴子石等,其次有石英、云母、方解石、电气石、钠长石、萤石,普遍含金属硫化物,呈囊状、透镜状、似层状产出,厚度一般小于60 m,走向延伸超过10 km以上。复杂矽卡岩是成矿区最重要的赋矿围岩,都龙、新寨、南秧田等超大型、大型矿床均产出于复杂矽卡岩及其附近变质岩层位。老君山燕山期花岗岩体按照其演化顺序大致可分为三个亚期:①第一亚期占岩体出露总面积的2/3,为灰白色含斑中粗粒二云母花岗岩,组成矿物有斜长石、微斜长石、黑云母、白云母、石英,为似斑状花岗结构,块状构造。②第二亚期约占该岩体总面积的1/3,为灰白色中细粒二云母花岗岩,呈岩株状侵入于第一期岩体中,主要矿物为微斜长石、更长石、石英、黑云母、白云母,具细粒花岗结构,块状构造。③第三亚期为灰白色花岗斑岩,呈岩枝及岩脉穿插于早期花岗岩及变质岩系中,斑晶主要为钾长石、石英,次为黑云母、斜长石,具斑状结构,块状构造。经野外实地调研及综合分析,都龙、南秧田矿区出露的第三期花岗斑岩脉对锡钨等元素的富集具有积极意义[20]。
3 矿床类型及特征
根据赋矿围岩的不同,老君山成矿区矿床类型可划分为与花岗岩有关的石英脉、伟晶岩型钨锡矿床,简单矽卡岩型锡矿床,片岩、复杂矽卡岩型锡钨多金属矿床,以及大理岩型银铅锌矿床等,含矿围岩主要为花岗岩及变质岩类。与花岗岩有关的石英脉、伟晶岩型钨锡矿床位于老君山成矿区中心近花岗岩体部位,典型矿点有花石头、戈岭、老寨等。矿体呈脉状产出于花岗岩体内部及其与围岩的接触带中,组成矿脉的主要矿物有锡石、黑钨矿、绿柱石、黄铁矿、黄铜矿、石英、长石、云母、电气石、石榴子石、绿泥石。
简单矽卡岩型锡矿床分布于成矿区的北东部。矿体产出地层为中寒武统田蓬组(∈2t),层状含矿简单矽卡岩广泛分布,长达30余km(巷香寨——茅草坪——三水),普遍发育锡矿化,矿物组合简单,一般不超过4种非金属矿物,主要金属矿物有锡石、黄铁矿,脉石矿物有透辉石、绿泥石、阳起石等,典型矿化点有茅草坪、三水等。
片岩、复杂矽卡岩型矿床可分为锡石多金属硫化物矿床及白钨矿床两类,是老君山成矿区的主要矿床类型。锡石多金属硫化物矿床主要分布于花岗岩体南西部的都龙及北东部的新寨,赋矿地层为中寒武统田蓬组,矿体产出于碎屑岩与碳酸盐岩过渡带,矿体走向长几十米至2 000余m,厚度数米至数十米,倾斜延伸几十米至数百米;主要金属矿物有锡石、铁闪锌矿、磁铁矿、磁黄铁矿、黄铜矿、黄铁矿、毒砂;脉石矿物有钙铁辉石、透辉石、绿泥石、阳起石、石榴石、透闪石、绿帘石、石英、白云母、萤石、方解石等。白钨矿矿床分布于花岗岩体东部的南秧田,赋矿地层为下寒武统冲庄组(∈1ch)浅灰绿色云母石英片岩、透辉石片岩、阳起石片岩,下部矿层多为似层状、层状矿体展布,走向长10余km,厚度数米至10余m,倾斜延伸百余米至千余米,白钨矿颗粒细小,上部矿层有脉状、透镜状矿体穿插,脉状矿体中白钨矿颗粒粗大;主要金属矿物有白钨矿、磁黄铁矿、黄铁矿、毒砂、黄铜矿及少量锡石,脉石矿物有石英、长石、云母、透辉石、透闪石、阳起石、斜黝帘石。
大理岩型银铅锌矿床为位于老君山花岗岩体北部新寨附近的三保银铅锌多金属矿床和都龙南部的南当厂银铅锌多金属矿床,矿体产出地层为中寒武统田蓬组,矿体呈透镜状、似层状产于碎裂似层状大理岩中,主要金属矿物有方铅矿、铁闪锌矿、黄铁矿、毒砂、辉银矿、自然银,脉石矿物有方解石、白云石、石英等。
为了研究矿床成因,作者所在研究团队曾进行了岩、矿石的光薄片鉴定工作,发现了大量早加里东期海底喷流沉积成矿作用证据及其后期叠加改造成矿作用证据[11, 21]。如大量纹层条带状构造,鲕状、结核状构造,变余玻屑填间结构、变余斑状结构、变余间粒结构等,表明早期喷流沉积成矿作用的存在,不同金属矿物间的穿插交代,亦表明后期变质、叠加改造成矿作用对成矿影响强烈。
本文围绕不同种类单矿物的生成顺序进行了综合鉴定工作,对60余件矿石光片进行鉴定和拍照,典型金属矿物组合特征如下:
脉状构造:毒砂强烈破裂变形,破裂缝隙间充填磁黄铁矿,并具明显交代现象,晚期细脉状黄铁矿穿插上述矿物,呈脉状构造(图 2a);铁闪锌矿、黄铜矿穿插交代碎裂黄铁矿,呈脉状构造(图 2b);锡石沿铁闪锌矿裂隙充填,形成脉状构造(图 2c)。
网脉状构造:方铅矿、铁闪锌矿,沿碎裂黄铁矿、毒砂的缝隙充填,形成网脉状构造(图 2d)。
残余状结构:毒砂、磁黄铁矿、铁闪锌矿被方铅矿交代呈残余结构(图 2e);铁闪锌矿及磁黄铁矿被黄铜矿穿插交代,呈残余状结构(图 2f)。
老君山成矿区金属元素单矿物的总体生成顺序由早及晚为:毒砂→磁黄铁矿→黄铁矿→铁闪锌矿→黄铜矿、方铅矿和锡石,较早形成的毒砂、磁黄铁矿是层纹条带状、块状硫化物矿石的重要组成部分,后期形成的铁闪锌矿、黄铜矿、方铅矿等多构成脉状、网脉状及细脉浸染状矿石。
4 样品采集和测试方法
本文围绕研究区不同矿化类型的典型矿床进行了系统取样工作。在花石头、老寨、箐角等矿点采集花岗岩型、斑岩型及接触带矽卡岩型矿石和围岩样品;在茅草坪矿化点采集远离老君山花岗岩体的简单矽卡岩型矿石和相关围岩样品;在都龙(铜街矿段、曼家寨矿段、万隆山矿段)、新寨、南秧田等矿床采集外接触带片岩型、复杂矽卡岩型矿石和围岩样品;在南当厂、坝脚采集大理岩型矿石和围岩样品。挑选50余件矿石通过破碎、挑选单矿物后用于硫、铅、氢、氧同位素测试工作。测试矿物包括磁黄铁矿、毒砂、方铅矿、黄铁矿、黄铜矿、铁闪锌矿、萤石、石英等,测试工作在国土资源部宜昌地质矿产研究所完成。22件单矿物硫同位素的测试方法是:将挑选的单矿物与氧化亚铜在真空状态下加热,进行氧化反应,生成二氧化硫,用气体同位素质谱仪分析硫同位素组成(质谱仪为:德国产MAT251)。得到的δ34S测定值重复样品分析的相对双差小于±0.5‰,其准确度由标准物质的分析结果间接给出,为±0.2‰。
29件单矿物铅同位素的测试方法是:采用静态模式,用标准样SRM981控制质谱计中的质量分馏,监控Pb同位素比值测定的NBS981国际标准物质的测定值207Pb/206Pb为0.914 2±0.000 5,检测温度为20 ℃,检测湿度为30%,总体误差小于0.026%。
10件单矿物氢、氧同位素的测试方法和依据为:DZ/T0184.19-1997天然水中氢同位素锌还原法测定,及DZ/T0184.13-1997硅酸盐和氧化物矿物氧同位素组成的五氟化溴法测定。测试结果以SMOW 为标准,分析精度优于±0.2‰。
5 结果与讨论 5.1 硫同位素
通过综合整理分析本文硫同位素测试结果及前人已有测试资料(表 1、图 3),老君山成矿区δ34S分布范围为-1.50‰~8.61‰,均值为2.06‰,方差为5.33‰,频率分布直方图(图 3a)显示,在其主体-2‰~6‰区间呈塔式分布。由于幔源硫δ34S约为0±3‰[22],总体分析结果与该范围较接近(略为偏高),结合成矿区地质特征,推测矿床硫主要为幔源硫,在多期地质作用过程中有一定量的壳源硫混入,且硫同位素发生了一定程度的均一化。有3个样品δ34S值分布在7‰~9‰,推测与海水硫的混入有关。所有统计样品不同单矿物δ34S均值由小到大排序为:毒砂、磁黄铁矿、黄铜矿、铁闪锌矿、方铅矿、黄铁矿。通过与前节所述矿物生成顺序(毒砂→磁黄铁矿→黄铁矿→铁闪锌矿→黄铜矿、方铅矿、锡石)对比,发现早期形成的毒砂(0.14‰)和磁黄铁矿(1.25‰)δ34S均值较低,且接近于0,幔源硫特征明显。晚期形成的铁闪锌矿(1.72‰)、黄铜矿(1.68‰)和方铅矿(2.91‰)δ34S均值略大,反映了壳源硫的混入,推测与后期叠加改造成矿作用有关。黄铁矿δ34S值最大,是因为相对其他几种金属硫化物,黄铁矿具有较高与H2S的反应平衡同位素分馏系数[23],比其他金属硫化物易于从外界环境中俘获34S。
样号 | 取样地点 | 矿石类型 | 测试矿物 | δ34S/‰ | 资料来源 | 样号 | 取样地点 | 测试矿物 | δ34S/‰ | 资料来源 |
MJZ-03 | 都龙 | 穿层石英脉型 | 磁黄铁矿 | 0.27 | 本文 | T2-1 | 铜街 | 磁黄铁矿 | 0.60 | 文献 |
DL01 | 都龙 | 穿层石英脉型 | 磁黄铁矿 | -0.07 | 都硫5 | 铜街 | 铁闪锌矿 | 2.60 | [3, 10] | |
DL09 | 都龙 | 层状矽卡岩型 | 磁黄铁矿 | -1.03 | 都硫6 | 铜街 | 铁闪锌矿 | 1.60 | ||
WLS-04 | 都龙 | 层状矽卡岩型 | 黄铁矿 | 1.56 | A1005-1 | 铜街 | 铁闪锌矿 | 1.20 | ||
HST-01 | 花石头 | 花岗岩型 | 毒砂 | 1.24 | A1010 | 铜街 | 黄铁矿 | 2.20 | ||
DL23 | 南当厂 | 大理岩型 | 黄铁矿 | 4.93 | A1028 | 铜街 | 黄铁矿 | 3.80 | ||
DL24 | 南当厂 | 大理岩型 | 黄铁矿 | 4.78 | A1005-3 | 铜街 | 铁闪锌矿 | 1.00 | ||
DL25 | 南当厂 | 大理岩型 | 方铅矿 | 1.08 | A1005-2-2 | 铜街 | 磁黄铁矿 | 1.00 | ||
NYT6 | 南秧田 | 石英脉型 | 磁黄铁矿 | 5.78 | 5901-1 | 曼家寨 | 磁黄铁矿 | 0.80 | ||
NYT8-1 | 南秧田 | 石英脉型 | 磁黄铁矿 | 4.95 | 2701-1 | 曼家寨 | 磁黄铁矿 | 0.40 | ||
XZKD02 | 新寨 | 石英脉型 | 毒砂 | -0.56 | 1101-1 | 曼家寨 | 磁黄铁矿 | -0.70 | ||
N-05 | 新寨 | 石英脉型 | 毒砂 | -1.18 | 都硫7 | 曼家寨 | 铁闪锌矿 | -1.50 | ||
XZ15 | 新寨 | 层状矽卡岩型 | 毒砂 | 1.05 | 都硫8 | 曼家寨 | 铁闪锌矿 | 1.00 | ||
QJ1(h) | 箐脚 | 层状矽卡岩型 | 黄铁矿 | 2.58 | 都硫10 | 曼家寨 | 铁闪锌矿 | 2.70 | ||
QJ2(h) | 箐脚 | 层状矽卡岩型 | 方铅矿 | 1.37 | 都硫11 | 曼家寨 | 黄铜矿 | 0.60 | ||
QJ9 | 箐脚 | 层状矽卡岩型 | 黄铁矿 | 3.90 | LDL-16 | 曼家寨 | 铁闪锌矿 | 0.50 | ||
LZ1(h) | 老寨 | 层状矽卡岩型 | 黄铜矿 | 0.61 | LDL-25 | 曼家寨 | 黄铜矿 | 0.90 | ||
LZ2 | 老寨 | 层状矽卡岩型 | 黄铁矿 | 2.18 | 都硫S | 辣子寨 | 铁闪锌矿 | 0.40 | ||
LZ7 | 老寨 | 石英脉型 | 黄铜矿 | 3.94 | L1-1 | 辣子寨 | 铁闪锌矿 | 1.80 | ||
QJ7 | 箐脚 | 石英脉型 | 黄铁矿 | 4.15 | L1-2 | 辣子寨 | 磁黄铁矿 | 1.80 | ||
QJ10 | 箐脚 | 斑岩型 | 黄铜矿 | 2.37 | LDL-95 | 辣子寨 | 铁闪锌矿 | 0.70 | ||
QJ11 | 箐脚 | 斑岩型 | 黄铁矿 | 8.61 | LDL-102 | 水洞厂 | 铁闪锌矿 | 7.50 | ||
不同矿石类 | 花岗岩及花岗斑岩型3件 | 4.07 | 本文 | LDL-104 | 水洞厂 | 方铅矿 | 8.30 | |||
型均值统计 | 层状矽卡岩型8件 | 1.53 | LDL-64(1) | 岩冲 | 黄铁矿 | 3.50 | ||||
切穿层状矽卡岩的石英脉型8件 | 2.16 | LDL-64(2) | 岩冲 | 铁闪锌矿 | 2.80 | |||||
大理岩型3件 | 3.60 | LDL-64(3) | 岩冲 | 方铅矿 | 0.90 | |||||
不同矿物 | 磁黄铁矿11件 | 1.25 | 本文 | |||||||
均值统计 | 毒砂4件 | 0.14 | ||||||||
方铅矿4件 | 2.91 | |||||||||
黄铁矿11件 | 3.83 | |||||||||
黄铜矿5件 | 1.68 | |||||||||
铁闪锌矿13件 | 1.72 |
本文测试的22件样品(图 3b)显示,层状矽卡岩型矿石中硫化物δ34S分布范围为-1.03‰~3.90‰,均值为1.53‰,方差为2.11‰,分布范围相对集中,硫同位素组成偏幔源。切穿矽卡岩层石英脉型矿石中硫化物δ34S分布范围为-1.18‰~5.78‰,均值为2.16‰,方差为7.87‰;δ34S分布可分为两类:一类分布在-2‰~1‰,反映了对幔源硫的继承;另一类为3‰~6‰,为壳源硫。花岗岩及斑岩型矿石中的硫化物δ34S均值为4.07‰,为壳源硫。大理岩型矿石中的硫化物δ34S均值为3.60‰,以壳源硫为主。
5.2 铅同位素 由本文铅同位素测试结果与前人资料(表 2)经综合统计分析可知,老君山成矿区岩石、矿石铅同位素μ值分布范围较广,为9.18~10.11,均值为9.65。具有高μ值(大于9.58)的铅通常被认为是来自U、Th 相对富集的上部地壳,低μ值(≤9.58)铅源自于下地壳或地幔[24, 25, 26]。统计结果显示,高μ值铅占66.7%,低μ值铅占33.3%。因此,据μ值判断,研究区岩、矿石铅源主要为上地壳。样品编号 | 取样位置 | 矿石类型 | 样品名称 | 206Pb/204Pb | 207Pb/204Pb | 208Pb/204Pb | 模式年龄/Ma | μ | 资料来源 |
JBJ1 | 坝脚 | 大理岩型矿石 | 方铅矿 | 18.080 | 15.641 | 38.187 | 453.1 | 9.58 | 本文 |
JBJ2 | 坝脚 | 大理岩型矿石 | 方铅矿 | 18.100 | 15.654 | 38.230 | 454.0 | 9.61 | |
WLS-04 | 都龙 | 层状矽卡岩型硫化矿石 | 黄铁矿 | 18.745 | 15.899 | 39.472 | 289.4 | 10.02 | |
HST-01 | 花石头 | 花岗岩裂隙石英脉型硫化矿石 | 毒砂 | 18.533 | 15.766 | 38.964 | 281.0 | 9.78 | |
LZ1(h) | 老寨 | 层状矽卡岩型硫化矿石 | 黄铜矿 | 18.921 | 15.654 | 39.613 | -140.0 | 9.52 | |
LZ2 | 老寨 | 层状矽卡岩型硫化矿石 | 黄铁矿 | 18.719 | 15.736 | 39.051 | 113.2 | 9.70 | |
LZ4 | 老寨 | 层状矽卡岩型硫化矿石 | 黄铁矿 | 18.890 | 15.845 | 39.300 | 125.4 | 9.89 | |
LZ7 | 老寨 | 穿层状矽卡岩石英脉型硫化矿石 | 黄铜矿 | 18.736 | 15.735 | 39.086 | 100.0 | 9.70 | |
MJZ-03 | 曼家寨 | 穿层状矽卡岩石英脉型硫化矿石 | 磁黄铁矿 | 18.413 | 15.633 | 38.768 | 206.0 | 9.53 | |
DL01 | 曼家寨 | 穿层状矽卡岩石英脉型硫化矿石 | 磁黄铁矿 | 18.471 | 15.653 | 38.943 | 189.0 | 9.56 | |
MXZ 1 | 茅草坪 | 条带状石榴石矽卡岩 | 矽卡岩 | 20.662 | 15.823 | 38.807 | -1 000.0 | 10.03 | |
MXZ 4 | 茅草坪 | 矽卡岩化角岩 | 矽卡岩 | 18.358 | 15.782 | 38.732 | 421.0 | 9.83 | |
MXZ 5 | 茅草坪 | 阳起石矽卡岩 | 矽卡岩 | 18.539 | 15.780 | 38.817 | 293.3 | 9.80 | |
MXZ 6 | 茅草坪 | 含角闪石矽卡岩 | 矽卡岩 | 20.209 | 15.720 | 38.969 | -1000.0 | 9.63 | |
DL23 | 南当厂 | 大理岩型矿石 | 黄铁矿 | 18.466 | 15.668 | 38.913 | 210.8 | 9.59 | |
DL24 | 南当厂 | 大理岩型矿石 | 黄铁矿 | 18.568 | 15.778 | 39.272 | 270.7 | 9.80 | |
DL25 | 南当厂 | 大理岩型矿石 | 方铅矿 | 18.503 | 15.682 | 38.951 | 201.4 | 9.62 | |
NYT6(一) | 南秧田 | 穿层状矽卡岩石英脉型硫化矿石 | 磁黄铁矿 | 18.526 | 15.537 | 38.614 | 0.0 | 9.33 | |
NYT8-1(一) | 南秧田 | 穿层状矽卡岩石英脉型硫化矿石 | 磁黄铁矿 | 19.037 | 15.961 | 39.566 | 162.0 | 10.11 | |
QJ1(h) | 箐脚 | 层状矽卡岩型硫化矿石 | 黄铁矿 | 18.521 | 15.758 | 39.171 | 280.0 | 9.76 | |
QJ2(h) | 箐脚 | 层状矽卡岩型硫化矿石 | 方铅矿 | 18.426 | 15.636 | 38.726 | 200.0 | 9.53 | |
QJ7 | 箐脚 | 穿层状矽卡岩石英脉型硫化矿石 | 黄铁矿 | 18.489 | 15.725 | 39.027 | 263.3 | 9.70 | |
QJ9 | 箐脚 | 层状矽卡岩型硫化矿石 | 黄铁矿 | 18.520 | 15.794 | 39.097 | 322.9 | 9.83 | |
QJ10 | 箐脚 | 花岗岩型、花岗斑岩型硫化矿石 | 黄铜矿 | 18.600 | 15.721 | 39.092 | 180.0 | 9.68 | |
QJ11 | 箐脚 | 花岗岩斑岩型硫化矿石 | 黄铁矿 | 18.620 | 15.718 | 39.058 | 161.8 | 9.67 | |
DL09 | 铜街 | 层状矽卡岩型硫化矿石 | 磁黄铁矿 | 18.465 | 15.670 | 38.917 | 214.0 | 9.60 | |
XZKD02 | 新寨 | 穿层状矽卡岩石英脉型硫化矿石 | 毒砂 | 19.432 | 15.542 | 43.101 | -700.9 | 9.28 | |
N-05 | 新寨 | 穿层状矽卡岩石英脉型硫化矿石 | 毒砂 | 18.545 | 15.637 | 39.153 | 115.2 | 9.52 | |
XZ15 | 新寨 | 层状矽卡岩型硫化矿石 | 毒砂 | 19.560 | 15.651 | 41.133 | -635.5 | 9.48 | |
ZK83049-56 | 都龙 | 花岗斑岩 | 黄铁矿 | 18.625 | 15.704 | 38.986 | 140.9 | 9.65 | 文献 |
LDL-76 | 都龙 | 矿石 | 磁黄铁矿 | 18.524 | 15.679 | 38.951 | 182.6 | 9.61 | [10] |
LDL-76 | 都龙 | 矿石 | 铁闪锌矿 | 18.527 | 15.694 | 39.043 | 198.9 | 9.64 | |
LDL-90 | 都龙 | 矿石 | 铁闪锌矿 | 18.533 | 15.685 | 39.126 | 183.6 | 9.62 | |
LDL-95 | 都龙 | 矿石 | 铁闪锌矿 | 18.474 | 15.652 | 38.848 | 185.4 | 9.56 | |
LDL-25 | 都龙 | 矿石 | 黄铜矿 | 18.560 | 15.771 | 39.226 | 268.0 | 9.78 | |
LDL-102 | 都龙 | 矿石 | 铁闪锌矿 | 18.515 | 15.717 | 39.082 | 235.3 | 9.68 | |
LDL-104 | 都龙 | 矿石 | 方铅矿 | 18.473 | 15.655 | 38.873 | 189.8 | 9.57 | |
LDL-105 | 都龙 | 矿石 | 黄铁矿 | 18.495 | 15.679 | 38.946 | 203.4 | 9.61 | |
Stp-1 | 都龙 | 花岗片麻岩 | 片麻岩 | 18.763 | 15.693 | 38.832 | 27.4 | 9.61 | |
DLG-30 | 都龙 | 花岗片麻岩 | 片麻岩 | 18.434 | 15.644 | 38.514 | 204.4 | 9.55 | |
LDL-57 | 都龙 | 花岗岩 | 花岗岩 | 18.404 | 15.674 | 38.869 | 262.4 | 9.61 | |
LDL-56 | 都龙 | 花岗岩 | 花岗岩 | 18.533 | 15.738 | 39.127 | 247.8 | 9.72 | |
LJS-5 | 都龙 | 花岗岩 | 花岗岩 | 18.662 | 15.708 | 39.025 | 119.3 | 9.65 | |
YX-4 | 都龙 | 花岗斑岩 | 花岗斑岩 | 18.623 | 15.760 | 39.187 | 210.6 | 9.76 | |
ZK83049-56 | 都龙 | 花岗斑岩 | 花岗斑岩 | 18.538 | 15.675 | 38.902 | 167.6 | 9.60 | |
ZP-8 | 都龙 | 片岩 | 片岩 | 19.032 | 15.760 | 39.931 | -82.3 | 9.72 | |
ZP-29 | 都龙 | 片岩 | 片岩 | 20.665 | 15.782 | 42.688 | -1 000.0 | 10.03 | |
P13-1 | 都龙 | 片岩 | 片岩 | 18.918 | 15.670 | 39.667 | -116.5 | 9.55 | |
ZP-21 | 都龙 | 大理岩 | 大理岩 | 18.330 | 15.627 | 38.605 | 258.5 | 9.53 | |
P12-1# | 都龙 | 大理岩 | 大理岩 | 17.744 | 15.513 | 38.057 | 543.6 | 9.37 | |
P12-2# | 都龙 | 大理岩 | 大理岩 | 18.658 | 15.681 | 39.219 | 88.5 | 9.60 | |
P13-2# | 都龙 | 大理岩 | 大理岩 | 17.314 | 15.496 | 37.808 | 829.2 | 9.41 | |
LDL-82# | 都龙 | 大理岩 | 大理岩 | 18.349 | 15.449 | 38.016 | 17.2 | 9.18 | |
LDL-103 | 都龙 | 大理岩 | 大理岩 | 18.481 | 15.653 | 38.865 | 181.6 | 9.56 |
研究区铅同位素组成总体变化性较大,其中206Pb/204Pb=17.314~20.665,方差为0.315,207Pb/204Pb=15.449~15.961,方差为0.008,208Pb/204Pb=37.808~43.101,方差为0.796。其中,206Pb/204Pb与208Pb/204Pb的相关系数为0.650,在铅同位素构造模式图(图 4a、b)上,样品点投影线性分布特征明显;206Pb/204Pb与207Pb/204Pb的相关系数为0.401,在构造模式图(图 4c、d)上,样品投影点较分散。
在208Pb/204Pb-206Pb/204Pb铅同位素构造模式图(图 4a、b)上,岩、矿石样品主要投影在上地壳和下地壳之间的造山带附近;而在207Pb/204Pb-206Pb/204Pb铅同位素构造模式图(图 4c、d)上,岩、矿石大部分样品投影在上地壳附近,少量样品分布在造山带附近。与前述据μ值特征推测一致,研究区铅来源以上地壳为主,同时存在部分上地幔或下地壳来源。
由铅同位素构造模式图(图 4)可知,穿层石英脉型矿石、矽卡岩、片岩、片麻岩及其相关矿石,一部分样品投影靠近下地壳和造山带区域,反映了对早期深源铅的继承,其他样品铅源投影在上地壳,反映了后期变质含矿热液或花岗岩含矿热液的叠加。花岗岩、花岗斑岩及其相关矿石以壳源铅为主。大理岩及其相关矿石样品多投影在造山带靠近下地壳的位置,显示了深源铅的特征。
按照单阶段演化模式计算的模式年龄,结果中除去8个负值,其余的46件样品的模式年龄分布范围为829.2~0 Ma。在铅同位素年龄分布直方图(图 5)上,各类岩、矿石年龄分布区间与影响研究区的多期地质事件相吻合。
460~420 Ma代表了早加里东期地质事件,研究区加里东期夹持在弥勒——师宗断裂和越北古陆之间,为扬子地台前缘海盆环境,伴随地幔上隆[15],海底地壳发生拉张、断陷,深部高温、高浓度地幔物质上升,引发该区海底喷流沉积成矿作用的发生,形成了众多的块状硫化物矿床或矿源层。
300~160 Ma与海西——印支期地质作用相对应,年龄统计结果(图 5)可分为300~240 Ma(海西期)和220~160 Ma(印支期)两个地质作用峰期,两期地质作用引起了文山——麻栗坡断裂、马关都龙断裂等断裂的压扭走滑,地壳上隆,并引发了大范围的地层变形变质。老君山变质核杂岩已有研究[18]表明,不同采样地点、不同变质岩类采用SHRIMP变质锆石U-Pb法测得年龄为(233.5±1.4) Ma、(236.4±5.8) Ma,IMS榍石U-Pb法测得年龄为(231.4±4.2) Ma,与海西期地质作用时期较为接近;采用Ar-Ar法测得不同变质岩坪年龄为(170.8±0.7) Ma、(205.2±1.2) Ma、(220.4±2.7) Ma和(195.7±6.6) Ma,与印支期地质作用时期基本吻合。
据新寨锡矿详查报告[27],采用K-Ar法获得与锡石共生的云母年龄为197.1~217.0 Ma;近年研究采用Ar-Ar法和Re-Os同位素定年方法,对新寨锡矿床中的金云母及南秧田钨矿床中的辉钼矿进行的成矿时代定年,金云母Ar-Ar同位素坪年龄为(209.5±1.1)Ma,辉钼矿的Re-Os同位素模式年龄为(209.1±3.3)~(214.1±4.3) Ma,说明新寨锡矿床和南秧田钨矿床形成于印支期[19]。可见,近年测得的成矿精确年龄集中在本文统计结果 印支期变质成矿作用(220~160 Ma)的早期。
120~80 Ma与燕山期花岗岩地质作用相对应。近年老君山花岗岩SHRIMP 锆石U-Pb定年,得到的岩体形成年龄为83.3~92.9 Ma[28, 29],分布于本文统计结果的晚期。老君山花岗岩体的侵位,在叠加改造早期矿床的同时,形成了部分花岗岩型、伟晶岩型及接触带型矿床或矿化点。
研究区矿石测年数据在印支期分布最为集中,说明印支期变质作用对成矿影响最为强烈。每期地质作用年龄分布都表现为先高频出现之后逐渐降低的趋势,推测与该期地质作用的先强烈影响而后逐渐减弱的演化特征相对应。
同位素测年统计结果在区域上可以与地质事件一一对应,但值得注意的是,研究区不同类型单矿物年龄均值的分布,与矿物生成先后顺序并无明显关联。同位素测年结果代表的是样品中同位素体系形成时的时间,也可以是原同位素组成体系被打破新体系形成的时间,由于研究区存在多期地质作用,可能有些地质作用仅仅造成相关矿物或岩石的同位素体系变更,改变了测年结果,但未造成大规模矿液的运移或成矿物质的堆积。因此,同位素的测年时间能否作为矿床形成时间,还需考虑该期地质作用对成矿物质的实地影响,就如对变质的寒武系岩石进行测年,得到的年龄可以有多种,但地层形成时代仍属寒武纪。受多期多阶段成矿作用影响的矿床,对成矿物质来源的判别应更注重物质而不是时间,更不能以较新的测年时间推翻较早期成矿作用的存在。
5.3 氢氧同位素
10件样品氢氧同位素测试结果(表 3)显示,δDSMOW值变化范围为-88.9‰~-59.5‰,均值为-73.42‰,有8件样品δDSMOW值范围为-40‰~-80‰,即分布于Ohmoto[30]和Sheppard[31]界定的岩浆水氢同位素分布范围,仅2件样品分布值略低于该范围,说明研究区成矿流体主要源自于岩浆。δ18OSMOW变化范围为-9.99‰~8.91‰,均值为-2.449‰,在氢氧同位素分布图(图 6)上,有1件样品落入岩浆水范围内,其余样品δ18OSMOW组成有向大气降水线漂移的现象,这说明成矿流体有大气降水的混入。野外编号 | 采样位置 | 矿石类型 | 测试矿物 | δDSMOW/‰ | δ18OSMOW/‰ |
JBJ1 | 坝脚铅矿 | 大理岩中的石英脉型硫化矿 | 石英 | -74.6 | -7.54 |
JBJ2 | 坝脚铅矿 | 大理岩中的石英脉型硫化矿 | 石英 | -86.1 | 3.42 |
DL01 | 曼家寨 | 层状矽卡岩型硫化矿 | 铁闪锌矿 | -62.0 | 1.92 |
MJZ03 | 曼家寨 | 切穿层状矽卡石英脉型硫化矿 | 石英 | -76.9 | -9.26 |
DL09 | 曼家寨 | 层状矽卡岩型硫化矿 | 萤石 | -68.4 | -6.47 |
B03 | 曼家寨 | 层状矽卡岩型硫化矿 | 铁闪锌矿 | -70.8 | -0.70 |
TJ02 | 铜街 | 层状矽卡岩型硫化矿 | 铁闪锌矿 | -75.3 | -1.00 |
WLS01 | 万隆山 | 层状矽卡岩型硫化矿 | 铁闪锌矿 | -59.5 | 8.91 |
HST01 | 花石头 | 花岗岩石英脉型硫化矿 | 石英 | -88.9 | -3.78 |
LZ7 | 老寨 | 花岗岩石英脉型硫化矿 | 石英 | -71.7 | -9.99 |
众多单矿物测年数据统计结果可以与各期地质事件一一对应,反映了各期地质作用对该区成矿的影响。但对于变质和叠加改造地质事件频发区域,测年结果不代表该时期必然发生成矿物质的巨量堆积,有关矿床成因的研究应更注重其他方面的证据,如矿体产状特征,矿石结构、构造特征,矿石的同位素、稀土、微量元素组成特征等。
通过对老君山成矿区矿床地质及硫、铅、氢、氧同位素地球化学方面的研究,综合分析前人研究成果,成矿区矿床的形成是主要受早加里东期海底喷流沉积、印支期区域变质和燕山期花岗岩侵位三期地质作用的制约。
硫同位素δ34S分布范围-1.50‰~8.61‰,均值为2.06‰,属壳幔混合来源;矿物生成顺序较早的单矿物硫同位素组成偏幔源,晚期为壳源;层状矽卡岩型矿石硫同位素组成偏幔源,花岗岩和花岗斑岩型为壳源硫,穿层石英脉型可分为幔源硫和壳源硫两类。
铅同位素分析结果属壳幔混合来源,以壳源为主。穿层石英脉型矿石、矽卡岩、片岩、片麻岩及其相关矿石,一部分样品投影靠近下地壳和造山带区域,反映了对早期深源铅的继承,其他样品铅源投影在上地壳,反映了后期变质热液或花岗岩热液的叠加。花岗岩、花岗斑岩及其相关矿石以壳源铅为主。大理岩及其相关矿石样品多投影在造山带靠近下地壳的位置,显示了地层沉积时深源铅的特征。
氢、氧同位素组成指示成矿流体以岩浆水为主,并有大气降水的混入。
[1] | 周建平, 徐克勤, 华仁民, 等. 滇东南锡多金属矿床成因商榷[J]. 云南地质, 1997, 16(4):309-349. Zhou Jianping, Xu Keqin, Hua Renmin, et al. A Discussion on Genesis of the Tin Polymetallic Sulfide Deposits of Southeastern Yunnan[J]. Yunnan Geology, 1997, 16(4):309-349. |
[2] | 丁星妤, 戴塔根, 刘晓玮. 滇东南南部多金属成矿带找矿潜力分析[J]. 吉林大学学报(地球科学版), 2012, 42(6):1730-1739. Ding Xingyu, Dai Tagen, Liu Xiaowei. Ore Prospecting Potential of the Polymetallic Metallogenic Belt in the South of SE Yunnan[J]. Journal of Jilin University(Earth Science Edition), 2012, 42(6):1730-1739. |
[3] | 宋焕斌. 云南东南部都龙锡石-硫化物型矿床的成矿特征[J]. 矿床地质, 1989, 8(4):29-38. Song Huanbin. Ore Froming Characteristics of the Dulong Cassiterite Sulfide Type Ore Deposit in Southeastern Yunnan[J]. Mineral Deposits, 1989, 8(4):29-38. |
[4] | 安保华. 老君山岩体特征、成因及其找矿意义探讨[J]. 西南矿产地质, 1990, 4(1):30-35. An Baohua. Discussion on the Characteristics, Genesis and Exploration Significance of the Laojunshan Granites[J]. Geology of Ore Deposits in Southwest China, 1990, 4(1):30-35. |
[5] | 官容生. 滇东南构造岩浆带花岗岩体的含矿性探讨[J]. 矿物岩石, 1991,11(1):92-101. Guan Rongsheng. An Approach of Mineralization of Granite Mass in the Structure Tectonic Magmatic Belt in the Southeast of Yunnan[J]. Minerals and Rocks, 1991,11(1):92-101. |
[6] | 付国辉. 云南都龙锡多金属矿床地质特征及成矿规律[J]. 西南矿产地质, 1992, 6(2):25-37. Fu Guohui. Geological Characteristics and Mineralization Regularity of the Dulong Sn Polymetallic Deposit, Yunnan Province[J]. Geology of Ore Deposits in Southwest China, 1992, 6(2):25-37. |
[7] | 王学焜. 麻栗坡新寨锡矿床地质地球化学特征[J]. 云南地质, 1994, 13(1):1-16. Wang Xuekun. Geological-Geochemical Characteristics of Xinzhai Tin Deposit in Malipo[J]. Yunnan Geology, 1994, 13(1):1-16. |
[8] | 周建平, 徐克勤, 华仁民,等. 滇东南喷流沉积块状硫化物特征与矿床成因[J]. 矿物学报, 1998, 18(2):158-168. Zhou Jianping, Xu Keqin, Hua Renmin, et al. Characteristics and Genesis of Exhalative Sedimentary Massive Sulfides in Southeastern Yunnan Province[J]. Acta Mineralogica Sinica, 1998, 18(2):158-168. |
[9] | 曾志刚, 李朝阳, 刘玉平, 等. 老君山成矿区变质成因矽卡岩的地质地球化学特征[J].矿物学报, 1999, 19(1):48-55. Zeng Zhigang, Li Chaoyang, Liu Yuping, et al. Geology and Geochemistry of Metamorphogenic Skarn from Laojunshan Metallogenic Province[J]. Acta Mineralogica Sinica, 1999, 19(1):48-55. |
[10] | 刘玉平, 李朝阳, 刘家军. 都龙矿床含矿层状矽卡岩成因的地质地球化学证据[J]. 矿物学报, 2000, 20(4):378-383. Liu Yuping, Li Chaoyang, Liu Jiajun. Characteristics and Genesis of Stratiform Skarn from Dulong Tin-Zinc Polymetallic Deposit[J]. Acta Mineralogica Sinica, 2000, 20(4):378-383. |
[11] | 秦德先. 近年在我国三个大型锡矿床中相继发现的火山(成矿)作用及其科学意义[J]. 昆明理工大学学报:理工版, 2008, 33(1):Ⅰ-Ⅱ. Qin Dexian. Discovery of Volcanic Mineralization of Three Giant Sn Deposits in Our Country in Recent Years and Its Science Significance[J]. Journal of Kunming University of Science and Technolog:Natural Science Edition, 2008, 33(1):Ⅰ-Ⅱ. |
[12] | 戴婕, 张林奎, 潘晓东, 等. 滇东南南秧田白钨矿矿床矽卡岩矿物学特征及成因探讨[J].岩矿测试, 2011, 30(3):269-275. Dai Jie, Zhang Linkui, Pan Xiaodong, et al. Mineralogical Characteristics and Genesis Discussion of Skarn in Nanyangtian Sheelite Deposite of Southestern Yunnan[J]. Rock And Mineral Analysis, 2011, 30(3):269-275. |
[13] | 石洪召, 张林奎, 任光明, 等. 云南麻栗坡南秧田白钨矿床层控似矽卡岩成因探讨[J]. 中国地质, 2011, 38(3):67-69. Shi Hongzhao, Zhang Linkui, Ren Guangming, et al. The Genesis of Skarnoid from the Nanyangtian Stratabound Scheelite Deposit in Malipo, Yunnan Province[J]. Geology in China, 2011, 38(3):67-69. |
[14] | 张世涛,冯明刚,吕伟. 滇东南南温河变质核杂岩解析[J]. 中国区域地质, 1998, 17(4):390-397. Zhang Shitao, Feng Minggang, Lü Wei. Analysis of the Nanwenhe Metamorphic Core Complex Southeastern Yunnan[J]. Chinese Regional Geology, 1998, 17(4):390-397. |
[15] | 李东旭, 许顺山. 变质核杂岩的旋扭成因:滇东南老君山变质核杂岩的构造解析[J].地质论评, 2000, 46(2):113-119. Li Dongxu, Xu Shunshan. Rotation-Shearing of Metamorphic Core Complex:Structure Analysis of Metam Orphic More Momplex in Laojunshan, Southeastern Yunnan Province[J]. Geological Review, 2000, 46(2):113-119. |
[16] | Yan D P, Zhou M F, Wang C Y, et al. Structural Styles and Chronological Evidences from Dulong-Song Chay Tectonic Dome:Earlier Spreading of South China Sea Basin due to Late Mesozoic to Early Cenozoic Extension of South China Block[J]. Earth Science-Journal of China University of Geosciences, 2005, 30(4):402-412. |
[17] | Yan D P, Zhou M F, Wang C Y, et al. Structural and Geochronological Constraints on the Tectonic Evolution of the Dulong-Song Chay Tectonic Dome in Yunnan Province, SW China[J]. Journal of Asian Earth Sciences, 2006, 28(4):332-353. |
[18] | 刘玉平,徐伟,廖震,等. 老君山变质核杂岩隆升的热历史解析与动力学机制探讨[J]. 矿物岩石地球化学通报, 2007, 26(增刊1):87-88. Liu Yuping,Xu Wei,Liao Zhen, et al. Thermal History Analysis and Dynamic Mechanism of Metamorphic Core Complex in Laojunshan Bulletin of Mineralogy[J]. Petrology and Geochemistry, 2007, 26(Sup.1):87-88. |
[19] | 冯佳睿, 毛景文, 裴荣富, 等. 滇东南老君山地区印支期成矿事件初探:以新寨锡矿床和南秧田钨矿床为例[J]. 矿床地质,2011, 30(1):57-73. Feng Jiarui, Mao Jingwen, Pei Rongfu, et al. Atentative Discussion on Indosinian Ore-Forming Events in Laojunshan Area of Southeastern Yunnan:A Case Study of Xinzhai Tin Deposit and Nanyangtian Tungsten Deposit[J]. Mineral Deposits, 2011, 30(1):57-73. |
[20] | 贾福聚. 云南老君山成矿区成矿系列及成矿规律研究[D]. 昆明:昆明理工大学, 2010:1-126. Jia Fuju. Researches on Metallogenic Series and Metallogenic Regularities in Laojunshan Metallogenic Belt, Yunnan Province[D]. Kunming:Kunming University of Science and Technology, 2010:1-126. |
[21] | 贾福聚, 高建国, 念红良,等. 滇东南老君山锡多金属成矿区含矿岩系稀土元素地球化学[J]. 矿物学报, 2013, 33(4):658-664. Jia Fuju, Gao Jianguo, Nian Hongliang, et al. REE Geochemistry of Ore-Bearing Rocks in the Laojunshan Tin Poly-Metallic Metallogenic Region, Southeastern Yunnan Province, China[J]. Acta Mineral Sinica, 2013, 33(4):658-664. |
[22] | Chaussidon M, Lorand J P. Sulphur Isotope Composition of Orogenic Spinel Herzolite Massifs from Ariege (North-Eastern Pyrenees, France):An Ion Microprobe Study[J]. Geochimica et Cosmochimica Acta, 1990, 54(10):2835-2846. |
[23] | 张理刚. 稳定同位素在地质科学中的应用:金属活化热液成矿作用及找矿[M]. 西安:陕西科学技术出版社, 1985:16-250. Zhang Ligang. Stable Lisotope Applications in Geological Science:Hydrothermal Mineralization of Metal Activation and Prospecting[M]. Xi'an:Shaanxi Science and Technology Press, 1985:16-250. |
[24] | Doe B R, Zartman R E. Plumbotectonics, the Pha-nerozoic[J]. Geochemistry of Hydrothermal Ore Deposits, 1979, 2:22-70. |
[25] | Zartman R E, Doe B R. Plumbotectonics:The Model[J]. Tectonophysics, 1981, 75(1):135-162. |
[26] | 吴开兴, 胡瑞忠, 毕献武, 等. 矿石铅同位素示踪成矿物质来源综述[J]. 地质地球化学, 2002, 30(3):73-81. Wu Kaixing, Hu Ruizhong, Bi Xianwu,et al. Ore Lead Isotopes as a Tracer for Ore-Forming Material Sources:A Review[J]. Geology-Geochemistry, 2002, 30(3):73-81. |
[27] | 西南地质勘探局306队. 新寨锡矿详查报告[R]. 昆明:云南省地质资料馆, 1988. Team 306 of Southwest Geological Exploration Bureau. Detailed Survey Report of Xinzhai Tin Deposit[R]. Kunming:Geological Archives of Yunnan Province, 1988. |
[28] | 刘玉平, 李正祥, 李惠民, 等. 都龙锡锌矿床锡石U-Pb年代学:滇东南白垩纪大规模花岗岩-成矿事件[J]. 岩石学报, 2007, 23(5):967-976. Liu Yuping, Li Zhengxiang, Li Huimin, et al. U-Pb Geochronology of Cassiterite and Zircon from the Dulong Sn-Zn Deposit:Evidence for Cretaceous Large-Scale Granitic Magmatism and Mineralization Events in Southeastern Yunnan Province, China[J]. Acta Petrologica Sinica, 2007, 23(5):967-976. |
[29] | 冯佳睿, 毛景文, 裴荣富, 等. 云南瓦渣钨矿区老君山花岗岩体的SHRIMP锆石U-Pb定年、地球化学特征及成因探讨[J]. 岩石学报,2010, 26(3):845-857. Feng Jiarui, Mao Jingwen, Pei Rongfu, et al. SHRIMP Zircon U-Pb Dating and Geochemical Characteristics of Laojunshan Granite Intrusion from the Wazha Tungsten Deposit, Yunman Province and Their Implications for Petrogenesis[J]. Acta Petrologica Sinica, 2010, 26(3):845-857. |
[30] | Ohmoto H. Stable Isotope Geochemistry of Ore Deposits[J]. Reviews in Mineralogy and Geochemistry, 1986, 16(1):491-559. |
[31] | Sheppard S M F. Characterization and Isotopic Variations in Natural Waters[J]. Reviews in Mineralogy and Geochemistry, 1986, 16(1):165-183. |
[32] | Sheppard S M F, Taylor H P. Hydrogen and Oxygen Isotope Evidence for the Origins of Water in the Boulder Batholith and the Butte Ore Deposits, Montana[J]. Economic Geology, 1974, 69(6):926-946. |