0 引言
辽宁青城子地区铅锌矿自明朝至今已发现400多年,是我国北方重要的铅锌多金属矿产地之一,富产铅、锌、金、银等多种金属。榛子沟铅锌矿床是青城子矿田内代表性矿床之一。自20世纪50年代,众多学者从岩石学、同位素地球化学和沉积学等方面对该区矿床成因进行了大量研究工作,并取得一系列成果。张朋[1]通过对青城子矿田内岩浆活动的研究,提出区内出露的黑云母花岗岩为成矿母岩,矿床属中温热液充填交代成因;刘国平[2]通过对辽东裂谷所经历的主要构造热事件研究,发现每次岩浆活动在空间上虽不重合,但是彼此相邻,进而提出同位成矿作用的观点;刘红霞等[3]通过对青城子矿石组构特征的研究,发现该矿床的形成与沉积作用、变质变形及热液叠加作用有关,提出该矿床为变质热液改造成因类型;刘志远等[4]通过对青城子地区硅质岩的成因研究以及区内同生断裂的发现,提出矿床早期存在海底喷流沉积成岩成矿作用;王秀福等[5]通过对青城子矿田地质、同位素和成矿流体特征等方面的系统研究,提出区内铅锌矿床属海相火山喷流沉积-变质-岩浆热液(侵入-陆相火山)叠加的再造型成因。目前多数研究者认为青城子地区铅锌矿床为沉积-变质-热液叠加成因,其产出受地层、岩浆和构造联合控制。本文在前人研究工作基础上,重点分析热液叠加期铅锌矿化特征,并通过流体包裹体和碳、氢、氧稳定同位素分析,对成矿流体来源进行分析和讨论。
1 区域地质背景
青城子铅锌矿田地处华北陆台辽东——吉南古元古代裂谷增生体内。该裂谷于古元古代期间经历了拉张裂解——快速沉降——回返挤压——隆升拆离——消亡全部演化过程,并于中生代期间经历了濒西太平洋构造带的叠加作用[6]。区内主要发育一套属绿片岩相-角闪岩相的巨厚辽河群变质岩系地层,由老到新地层依次为浪子山组、里尔峪组、高家峪组、大石桥组和盖县组。浪子山组岩性以千枚岩、绿泥绢云片岩为主,原岩为一套陆源碎屑岩建造;里尔峪组岩性以变质流纹岩、变粒岩、大理岩为主,原岩为一套基性——酸性火山-沉积岩建造;高家峪组岩性以变粒岩、板岩为主,原岩为一套陆源碎屑岩黏土建造;大石桥组岩性以大理岩为主,原岩为一套碳酸盐岩建造;盖县组岩性以二云片岩为主,原岩为一套黏土岩建造。地层呈东西展布,整合接触。区内构造活动强烈,以褶皱和断裂构造为主。早期矿田经历辽吉裂谷带的发展演化过程,形成了大型翻卷褶皱和同生断裂[4](荒甸子——新岭——朱家堡子断裂)、韧性断裂等构造;晚期经历变质作用及印支期和燕山期的构造岩浆活动,形成了一系列的叠加褶皱和韧性、脆性断裂。叠加褶皱构造由轴向近EW、NW和NE的褶皱构造组成,其中EW向褶皱构造为主体。断裂构造根据走向方向可划分为NW、NE、SN和EW向四组断裂。NW和NE向断裂构成区内主要菱形断裂构造格架(图 1),规模较大,是脉状矿体的主要导矿及容矿构造;近SN和EW向断裂为次生断裂构造,是矿田内层状铅锌矿体和金银矿体的主要容矿构造[8, 9]。
区内岩浆活动频繁,其出露面积占全区1/3以上,依据岩石类型及同位素年代学将其划分为四期,分别为:古元古代岩浆活动,发育钾质花岗岩(条痕状花岗岩),岩浆活动剧烈,主要分布在南部,形成于2 000 Ma年左右[9];中元古代岩浆活动,发育钠质花岗岩,岩浆活动较弱,呈稀散点状分布,以东南部的大顶子钠质花岗岩岩群为代表,见铅锌矿化、铁锰碳酸盐化和高岭土矿化,形成于1 621~1 741 Ma[2];印支期岩浆活动,发育似斑状花岗岩-似斑状二长花岗岩,岩浆侵入活动强烈,规模大,分布在青城子矿区北部和西南部,以新岭花岗岩和双顶沟花岗岩为代表,形成于200 Ma左右[9];燕山期岩浆活动,发育安山质火山岩及伟晶岩脉、煌斑岩脉、闪长岩脉、花岗斑岩脉等脉岩,规模小,岩脉附近常发育脉状铅锌矿化或切割脉体,分布在矿区东南部钾质花岗岩和钠质花岗岩之间,形成于137 Ma左右[10](图 1)。 2 矿床地质特征 2.1 矿区地质 2.1.1 地层
矿区出露地层为太古宙鞍山群和古元古代辽河群变质岩系的高家峪组(Pt1g)、大石桥组(Pt1d)和盖县组(Pt1gx)[11]。其中,高家峪组和大石桥组为主要赋矿地层。高家峪组岩石主要由片岩、变粒岩、透闪岩、大理岩及斜长角闪岩组成。石墨、黄铁矿的含量较高,尤其是在层间滑动面附近,石墨体积分数可达15%以上。该组是矿区内层状铅锌矿体的主要赋矿层位。
大石桥组岩石主要由大理岩、变粒岩、片岩及片麻岩等组成。根据岩性组合划分为3个岩性段:一段主要为白云质大理岩夹厚层菱镁矿,沿走向分布稳定,富铅、锌、金、银等多金属矿和菱镁矿,是榛子沟区2号矿体的主要赋矿层位;二段由片岩和变粒岩等组成,沿地层延伸方向稳定,以含石榴石为主要特征,富含金矿;三段占矿区总面积的一半以上,以钙质大理岩为主,富铅、锌、金、银等多种金属和菱镁矿,其中该层位脉状矿化发育,为矿区内主要的含铅锌矿层之一。
2.1.2 构造
榛子沟——甸南倾没背斜是矿床的主要褶皱构造,褶皱两翼宽缓,轴向近EW,倾角70°,地层倒转,核部出露辽吉花岗岩(混合变粒岩),两翼岩性由核部向外依次为高家峪组角闪岩、大石桥组大理岩和盖县组片岩。区内断层发育,以脆性断裂为主,可划分为NW、NE和SN三组断裂。NW和NE断裂主要位于高家峪组与花岗岩、大石桥组、盖县组的接触部位。其中:NW向断裂倾向NE,倾角20°~60°,层状矿体发育;SN向断裂切穿地层,并切割NW、NE向断裂,构成脉状、囊状矿体的容矿空间[6]。2.1.3 岩浆岩
区内经历多期岩浆活动,主要发育古元古代钾质花岗岩、印支期似斑状花岗岩-似斑状二长花岗岩和燕山期岩脉,其中以古元古代钾质花岗岩和燕山期的脉岩为主。钾质花岗岩发育榛子沟——甸南倾没背斜核部,与两翼地层呈整合接触,岩石呈肉红色,斑状-似斑状结构和块状构造;煌斑岩、闪长岩脉主要发育在断裂内,与脉状矿体关系密切。2.2 矿化特征
矿体主要赋存于高家峪组和大石桥组及其与辽吉花岗岩接触附近的层间断裂内,少量矿体赋存于二者附近的切层断裂内。依据矿体的产状和形态,矿石类型可划分为层状和脉状矿体两类。层状矿体呈层状、似层状和透镜状产出,与容矿围岩的产状一致或近于一致,主要为喷流沉积成因类型,具有规模大和沉积结构特征,受地层控制;少量的层状矿体是后期热液顺片理贯入形成的,属后期热液成因,规模小且明显受断裂控制。本文主要研究喷流沉积成因的层状矿体。脉状矿体位于层状矿体上盘,呈脉状、囊状产出,产状与围岩不协调,穿切围岩的层理。矿石矿物以黄铁矿、闪锌矿和方铅矿为主,次为黄铜矿、毒砂、磁黄铁矿;脉石矿物以方解石、石英和白云石为主,次为云母类矿物。常见矿石结构为交代残余结构、压裂结构、乳滴状结构等(图 2),构造为块状构造、浸染状构造、变余构造和层纹状构造等。
2.3 围岩蚀变特征 围岩蚀变较弱,规模很小,分布不均匀,分带性不明显,主要有白云石化、矽卡岩化、硅化、碳酸盐化、黄铁矿化、绢云母化和绿泥石化等。其中,碳酸盐岩类型围岩蚀变以白云石化和碳酸盐化为主,云母片岩类围岩常见硅化、绢云母化和绿泥石化。层状铅锌矿体相对脉状铅锌矿体蚀变较弱。2.4 热液叠加期成矿作用特征
已有研究表明,榛子沟铅锌矿床的形成经历了海底喷流、变质变形和热液叠加三期成矿作用[1, 2, 3, 4, 5, 6, 7, 8, 9],其中热液叠加成矿作用对矿床的形成起了重要作用。根据矿石矿物组合、结构构造和矿脉穿切关系将其划分为I黄铁矿-方铅矿-闪锌矿-石英和Ⅱ黄铁矿-方铅矿-石英-方解石两个阶段。热液叠加成矿作用的最显著特征是脉状矿体的形成和层状矿体的局部热液改造(图 3)。本区脉状矿体多赋存于层状——似层状矿体上盘大理岩、各类片岩地层之中,沿不同方向、规模及性质的断裂构造产出(图 4)。矿体形态多为脉状、树枝状或N字型等,走向为NW(300°~330°)或近SN,倾角变化较大,矿体走向延长60~200 m,倾斜延深30~60 m。I黄铁矿-方铅矿-闪锌矿-石英阶段矿石矿物以黄铁矿、闪锌矿和方铅矿为主,次为黄铜矿,脉石矿物以石英为主。Ⅱ黄铁矿-方铅矿-石英-方解石阶段矿石矿物以方铅矿、黄铁矿为主,脉石矿物以方解石、石英为主(图 2)。品位数据(表 1)表明,脉状矿体w(Pb+Zn)为14.75%,似层状矿体w(Pb+Zn)为4.34%~16.16%,大部分小于14.75%,只有2号矿体的w(Pb+Zn)大于脉状矿体,但总体上脉状矿体w(Pb+Zn)值大于似层状矿体;综合推断热液叠加期明显对成矿元素起到富集作用。
矿体号 | 矿化类型 | 品位/% | 平均品位/% | Pb/Zn | w(Pb+Zn)/% | ||
Pb | Zn | Pb | Zn | ||||
289号 | 似层状 | 0.45~20.34 | 0.95~12.50 | 2.33 | 2.01 | 1.16 | 4.34 |
2号 | 似层状 | 0.43~40.90 | 0.87~69.71 | 6.04 | 10.12 | 0.60 | 16.16 |
321号 | 似层状 | 0.78~11.23 | 0.58~5.12 | 5.23 | 5.45 | 0.96 | 10.68 |
320号 | 似层状 | 0.56~20.15 | 0.78~28.00 | 3.54 | 7.65 | 0.46 | 11.19 |
282号 | 脉状 | 0.87~48.76 | 0.56~28.54 | 5.78 | 8.97 | 0.64 | 14.75 |
1号 | 脉状 | 0.67~23.54 | - | 6.43 | - | - | - |
注:资料来源于文献[6]。 |
显微测温实验和拉曼光谱成分分析在吉林大学地球科学学院流体包裹体实验室完成。测温仪器包括德国Carl Zeiss Axiolab型显微镜(10×50)和英国Linkam THMS-600型冷热两用台,温度低于0 ℃时分析精度为±0.1 ℃,温度高于200 ℃时为±2 ℃;成分分析仪器采用分析精度为1 cm-1的Renishaw System-100 型激光拉曼光谱仪和扫描范围为4 500~8 501 cm-1的514 nm Ar+离子激光器。
碳、氢、氧同位素测试在核工业地质矿产研究所实验室中心完成,测量仪器为MAT253 EM型质谱计,氧同位素分析精度为±0.2,碳同位素分析精度为±0.2,氢同位素分析精度为±2。BrF5法制取氧,氢同位素采用爆裂和锌法制取氢,碳同位素采用爆裂和冷阱法制取碳。
3.2 流体包裹体岩相学特征
Ⅰ阶段石英颗粒中见4种类型流体包裹体:纯液相包裹体、气液两相包裹体及少量CO2包裹体和富气相包裹体,其中纯液相包裹体分布极少,不参与本文讨论。岩相特征(图 5a——f)如下:气液两相包裹体(VL):常温下由H2O的气液两项组成,气相占包裹体整体的5%~15%;呈近圆状、椭圆状和少量不规则状,长轴多为5~8 μm,偶见10~15 μm,在石英颗粒中多为成群分布,少量随机分布。
富气相包裹体(LV):常温下由H2O的气液两项组成,气液比为85%~90%;形态多为近圆状,长轴多见为5 μm,在石英颗粒中随机分布且少见。
CO2三相包裹体(LC):常温下由H2O的液相和CO2的气液两项组成,CO2气液比范围较大,为40%~90%;形态近椭圆,长轴为5~10 μm,在石英颗粒中随机分布且少见。
Ⅱ阶段石英颗粒中流体包裹体以液相包裹体为主,气液两相包裹体次之。气液两相包裹体(VL)岩相特征(图 5g——i)表现如下:常温下由H2O的气液两项组成,气相占包裹体的5%~15%,呈近圆状、近椭圆状和长方形状及少量不规则状,长轴为5~10 μm,成群分布。 3.3 流体包裹体显微测温特征
岩相学研究发现Ⅰ阶段铅锌矿体发育气液两相、富气相和CO2三相流体包裹体,Ⅱ阶段矿体发育气液两相包裹体。本次对VL、LV、LC包裹体进行显微测温研究,结果如表 2所示。成矿阶段 | 包裹体类型 | 大小/μm | 气液比/% | 固相CO2熔化温度/℃ | 笼形物消失温度/℃ | CO2部分均一温度/℃ | 冰点温度/℃ | 包裹体完全均一温度/℃ | w(NaCl)/% | 密度/(g·cm-3) |
Ⅰ | VL(65) | 5~15 | 5~15 | -6.3~-2.6 | 160~310 | 4.3~9.7 | 0.80~0.96 | |||
LV(12) | 5~10 | 85~90 | -2.1~-1.8 | 270~320 | 3.0~3.7 | 0.65~0.80 | ||||
LC(8) | 5~10 | 40~90 | -58.8~-58.1 | 6.9~7.3 | 23.1~24.2 | 300~360 | 5.2~5.9 | 0.72~0.89 | ||
Ⅱ | VL(40) | 5~10 | 5~15 | -0.9~-0.8 | 120~270 | 1.4~9.4 | 0.87~1.01 | |||
注:括号内数字为包裹体数。气液两相包裹体盐度计算据Hall(1988);二氧化碳三相包裹体盐度计算据Roedder(1984)。 |
Ⅰ阶段气液两相包裹体盐度为4.3%~9.7%,均一至液相温度为160~310 ℃,偶见410 ℃烧爆温度,峰值为210~260 ℃;富气相包裹体盐度为3.0%~3.7%,均一至液相温度为270~320 ℃,峰值为290~310 ℃;CO2三相包裹体的固相CO2熔化温度为-58.8~-58.1 ℃,盐度为5.2%~5.9%,均一至液相温度为300~360 ℃(图 6a、b)。
Ⅱ阶段气液两相包裹体盐度为1.4%~9.4%,均一至液相温度为120~270 ℃,峰值为180~220 ℃(图 6c、d)。
3.4 流体包裹体成分特征 激光拉曼光谱分析研究结果(图 7)表明:Ⅰ黄铁矿-方铅矿-闪锌矿-石英阶段CO2三相包裹体液相成分为H2O和CO2,气相成分为H2O和CO2、CH4和N2,其中CO2、CH4和N2的含量与CO2气液比正相关;Ⅰ阶段气液两相包裹体液相成分为H2O,气相成分以H2O为主,CO2次之,CH4和N2再次之。 3.5 流体包裹体碳、氢、氧同位素特征 测试矿物石英中13C和D来自石英中流体包裹体封存的挥发分流体,如H2O、CO2、CH4;由于水岩反应的存在而测试值18O偏高,因此测试值13C和D的范围更能够反映成矿流体的来源,所测石英的18O同位素值需要基于流体温度进行转换(表 3)。样品号 | 成矿阶段 | δ18Oq-SMOW/‰ | δDq-SMOW/‰ | δ13C/‰ | 均一温度/℃ | δ18OH2O-SMOW/‰ |
QCZ-4 | Ⅰ黄铁矿-方铅矿-闪锌矿-石英阶段 | 10.0 | -95.4 | -4.4 | 227 | -0.62 |
QCZ-5 | 9.5 | -96.5 | 250 | 0.04 | ||
QCZ-6 | 10.0 | -4.8 | ||||
QCZ-15a | Ⅱ黄铁矿-方铅矿-石英-方解石阶段 | 7.7 | -80.9 | -7.9 | 160 | -7.43 |
QCZ-14 | 8.0 | -86.5 | 180 | -5.57 | ||
QCZ-16 | 7.2 | -87.5 | -10.8 | 160 | -7.93 | |
QCZ-18 | 7.0 | -88.4 | 180 | -6.57 | ||
QCZ-15b | 7.5 | -82.9 | -12.6 | 160 | -7.63 | |
QCZ-17 | 7.5 | -80.0 | 180 | -6.07 | ||
注:据文献[12, 13] 。δ18Oq-SMOW-δ18OH2O-SMOW=3.38×106/T2-2.9(T为均一温度绝对温度)。 |
Ⅰ阶段石英的δ18Oq-SMOW值为9.5‰~10.0‰,δDq-SMOW值为-96.5‰、-95.4‰,依据公式[12, 13]计算相应的δ18OH2O-SMOW值为-0.62‰、0.04‰,在体系δDq-SMOW-δ18OH2O-SMOW中投点落于大气降水线与岩浆水之间,偏于岩浆水一侧(图 8),表明其成矿流体具有岩浆水与大气降水的混合热液特点[16, 17, 18];δ13C值为-4.8‰、-4.4‰,落于花岗岩(δ13C为-7‰~-3‰)范围内(图 9),指示成矿流体中碳主要来源于岩浆水。
Ⅱ阶段石英的δ18Oq-SMOW值为7.0‰~8.0‰,δDq-SMOW值为-88.4‰~-80.0‰,相应δ18OH2O-SMOW为-7.93‰~-5.57‰,在体系δDq-SMOW-δ18OH2O-SMOW中投点落于大气降水线右侧(图 8),推断成矿流体可能主要来源于大气降水,但不排除有少量的岩浆水参与;δ13C值为-12.6‰~-7.9‰,较低于花岗岩碳值(图 9),推断碳具有多源性,在围岩蚀变作用(碳酸盐化)、含矿岩层中有机质(石墨)的氧化作用和岩浆作用不同程度的混合作用下表现出略低于岩浆碳的现象。
4 成矿流体性质及来源
流体包裹体岩相、显微测温及碳氢氧同位素特征表明,热液叠加期Ⅰ和Ⅱ两个阶段的成矿流体来源存在差异(图 6、8、9),故其成矿流体来源不同且具有不同的地球化学性质。Ⅰ阶段铅锌矿体的流体包裹体研究发现,富气相包裹体、CO2三相包裹体和气液两相包裹体共存于矿体中,气相成分多样富含H2O、CO2、CH4、N2等。Ⅱ阶段的铅锌矿体发育较为单一的气液两相包裹体,气相成分以H2O为主。Ⅰ阶段矿体成矿流体属中低温、低盐度、低密度的CO2-H2O-NaCl体系热液(图 10a),Ⅱ阶段成矿流体由于温度压力的降低,属低温、低盐度、低密度的H2O-NaCl体系(图 10b)。
Ⅰ阶段碳同位素值投点落于花岗岩内和氧同位素值向大气降水方向漂移现象,均反映成矿流体来源于岩浆水与大气降水的混合热液。Ⅱ阶段包裹体的盐度与温度呈现某种正相关关系(图 10b),即盐度随着温度的降低而降低,推断成矿后期有大气降水加入。同位素分析结果从另一方面证实大气降水参与成矿作用,并萃取地层中有用物质。K-Ar测年结果表明,多处切割脉状矿体的煌斑岩脉年龄为84~100 Ma,多处被脉状富矿化切割的闪长岩脉年龄为140 Ma[6],即脉状铅锌矿体应该形成于84~140 Ma;推断该成矿期的成矿作用与第四期岩浆活动,即形成与137 Ma的燕山期岩浆构造活动密切相关。综合认为热液叠加成矿期成矿流体来源于燕山期岩浆热液和大气降水,大气降水在成矿后期参与的比重增大。 5 结论
1)榛子沟铅锌矿床热液叠加期成矿作用明显对成矿元素起到富集作用,且依据矿物共生组合、矿石组构及穿切关系,热液叠加期可划分为Ⅰ方铅矿-闪锌矿-黄铁矿-石英阶段和Ⅱ黄铁矿-方铅矿-石英-方解石阶段。2)Ⅰ阶段铅锌矿体石英颗粒中发育气液两相、富气相和CO2三相包裹体,成矿流体属中低温、低盐度、低密度的CO2-H2O-NaCl体系热液,气相成分以H2O为主,次为CH4、N2、CO2;Ⅱ阶段铅锌矿体石英颗粒中发育较单一的气液两相包裹体,成矿流体属低温、低盐度、低密度的H2O-NaCl体系热液,气相成分以H2O为主。
3)成矿流体来源于燕山期的岩浆水与大气降水的混合热液,且在成矿后期大气降水的混入比例增加。
[1] | 张朋.浅谈岩浆热液活动在青城子矿床中的成矿作用[J].中国高新技术企业,2009(7):127-128. Zhang Peng.Mineralization of Magmatic Hydrothermal in Qingchengzi Ore Deposits[J].Chinese Hi-Tech Enterprises,2009(7):127-128. |
[2] | 刘国平.辽宁青城子矿田的同位成矿作用[J].有色金属矿产与勘查,1999,8(5):277-282. Liu Guoping.Isospatial Metallogenesis in Qingchengzi Ore Fild,Liaoning[J]. Geological Exploration for Non Ferrous Metals,1999,8(5):277-282. |
[3] | 刘红霞,孔含泉,杨言辰.辽宁小佟家堡子金矿床地质特征及成因研究[J].黄金,2006,27(5):13-16. Liu Hongxia,Kong Hanquan,Yang Yanchen.Geologic Characteristics and Genesis of Xiaotongjiapuzi Gold Deposit,Liaoning Province[J].Gold,2006,27(5):13-16. |
[4] | 刘志远,徐学纯,田豫才,等.辽东青城子地区喷流-沉积成矿作用与金银多金属矿化的关系[J]. 矿床地质, 2007, 26(5):563-571. Liu Zhiyuan,Xu Xuechun,Tian Yucai,et al.Relationship Between Sedimentation-Exhalation Ore-Forming Process and Gold-Silver Polymetallic Mineralization in Qingchengzi Area,Liaoning Province[J].Mineral Depodits,2007,26(5):563-571. |
[5] | 王秀福,刘培栋,杨桂莲,等.青城子铅锌矿田地质特征[J].有色矿冶,2010,26(2):2-21. Wang Xiufu,Liu Peidong,Yang Guilian,et al.Study on Process Mineralogy for a Certain Slag[J]. Non-Ferrous Mining and Metallurgy,2010,26(2):2-21. |
[6] | 董存杰.青城子锌金银多金属矿田矿床地质特征及成矿系统分析[D]. 北京:中国地质大学,2012. Dong Cunjie.Geological Characteristics of the Deposits and Analysis of the Mineralization System of Qingchengzi Pb-Zn-Au-Ag Polymetallic Ore Field[D]. Beijing:China University of Geosciences,2012. |
[7] | 张秋生.中国早前寒武纪地质及成矿作用[M]. 长春:吉林人民出版社,1984. Zhang Qiusheng.Chinese Early Precambrian Geology and Mineralization[M].Changchun:Jilin People's Publishing House,1984. |
[8] | 孙立民.青城子铅锌矿田闪锌矿特征及意义[J].辽宁地质,1997(3):210-217. Sun Limin.Characteristics and Significance of Sphalerite Qingchengzi Pb-Zn Ore-Field[J].Liaoning Geology,1997(3):210-217. |
[9] | 李基宏.辽宁青城子铅锌金银矿集区成矿条件与成矿预测[D].长春:吉林大学,2005. Li Jihong.Study on Ore-Forming Conditions and Mineral Resource Assessment of Lead-Zinc-Silver-Gold Metallogenic Belt in Qingchenzi,Liaoning Province[D].Changchun:Jilin University,2005. |
[10] | 代军治.辽宁青城子地区金、银矿床成矿流体特征及成因探讨[D].长春:吉林大学,2005. Dai Junzhi. Characteristics of Ore-Forming and Discussion on the of Au, Ag Deposits in Qingchengzi Region,Liaoning Province[D].Changchun:Jilin University,2005. |
[11] | 刘洪津. 青城子铅锌矿田及其外围金银矿床成生关系研究[J]. 矿产与地质,2012(6):476-479. Liu Hongjin. The Relationship Between Students and Research into the Qingchengzi Pb-Zn Orefield and Its Peripheral Gold Silver Deposit[J].Mineral Resources and Geology,2012(6):476-479. |
[12] | Taylor H P. The Application of Oxygen and Hydrogen Isotope Studies to Problems of Hydrothermal Alteration and Ore Deposition[J].Economic Geology,1997,69:843-883. |
[13] | Clayton R N, O'neil J R, Mayeda T K. Oxygen Isotope Exchange Between Quartz and Water[J].Journal of Geophysical Research,1972,77:3057-3067. |
[14] | 张理刚.稳定同位素在地质科学中的应用[M].西安:陕西科学技术出版社,1985. Zhang Ligang.The Application of the Stable Isotope to Geology[M].Xi'an:Shaanxi Science and Technology Press,1985. |
[15] | 郑永飞,陈江峰.稳定同位素地球化学[M].北京:科学出版社,2000. Zheng Yongfei,Chen Jiangfeng. The Application of the Stable Isotope to Geology[M].Beijing:Science Press,2000. |
[16] | Bodnar R J.Revised Equation and Table for Determining the Freezing Point Depression of H2O-NaCl Solution[J].Geochimica et Cosmochimica Acta,1993,57(3):683-684. |
[17] | 李文昌,尹光侯,余海军,等. 云南普朗斑岩型铜矿床成矿流体特征及矿床成因[J]. 吉林大学学报(地球科学版),2013,43(5):1436-1447. Li Wenchang,Yin Guanghou,Yu Haijun,et al.Cha-racteristics of the Ore-Forming Fluid and Genesis of the Pulang Copper Deposit in Yunnan Province[J].Journal of Jilin University(Earth Science Edition),2013,43(5):1436-1447. |
[18] | 王承洋,王可勇,周向斌,等. 内蒙古东山湾钨钼多金属矿床成矿流体地球化学特征及成因[J]. 吉林大学学报(地球科学版),2015,45(3):759-771. Wang Chengyang,Wang Keyong,Zhou Xiangbin,et al.Geochemical Characterists of Ore-Forming Fluid and Genesis of Dongshanwan Tungesten-Molybdenum Polymetallic Deposit in Inner Mongolia[J]. Journal of Jilin University(Earth Science Edition),2015,45(3):759-771. |