1. 中国科学院大学地球与行星科学学院 中国科学院计算地球动力学重点实验室, 北京 100049;
2. 北京科技大学土木与资源工程学院, 北京 100083;
3. 中国科学院地质与地球物理研究所, 北京 100029
2018年5月23日 收稿; 2018年7月11日 收修改稿
基金项目: 全国危机矿山接替资源综合研究项目(20089931)和国家科技支撑计划项目(2011BAB03B09)资助
通信作者: 李永兵,E-mail:yongbingli@ucas.ac.cn
摘要: 成矿热力学参数是矿床成因研究的重要内容,共生的矿物必然能共同反映其形成时的成矿热力学环境。基于此,通过分析兰家沟钼矿不同成矿阶段的矿物共生组合,利用流体包裹体测定不同成矿阶段的成矿温度和压力范围,并基于热力学计算获得不同成矿阶段的矿物共生组合相图,进而得到兰家沟钼矿床3个成矿阶段的成矿热力学参数。在石英斑晶阶段,氧逸度的区间范围为10-40~10-25,硫逸度的区间范围为10-40~10-5。在斑岩型矿化阶段,pH值为6.9~6.1,氧化还原电位在-1.3~0.3 V。在石英大脉型矿化阶段,氧逸度应小于10-13,硫逸度为10-14~10-6,矿床的成矿环境属于弱还原环境。这些研究结果有助于对兰家沟钼矿成矿环境的认识。
关键词: 兰家沟钼矿pH值硫逸度氧逸度矿物共生组合流体包裹体
Study on metallogenic thermodynamic parameters of Lanjiagou molybdenum deposit
GUO Yang1, LIU Lidong2, CHU Shaoxiong3, LI Yongbing1, YAN Quanren1, SHI Yaolin1
1. Key Laboratory of Computational Geodynamics of CAS, College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China;
2. School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China;
3. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
Abstract: The thermodynamic parameters of ore-forming process are the important content of the genesis of the deposits, and the paragenetic minerals inevitably show the metallogenic thermodynamic environment during their formation. Based on this, the mineral assemblages at different metallogenic stages in Lanjiagou molybdenum deposit are evaluated and three mineralization stages are determined. Then the ore-forming temperatures and pressure ranges of the three mineralization stages are obtained by measuring the temperatures of the inclusions, the phase diagrams of the paragenetic compositional minerals at the different mineralization stages are calculated, and the metallogenic thermodynamic parameters of the three metallogenic stages in Lanjiagou molybdenum deposit are estimated. At the stage of quartz phenocryst, the oxygen fugacity ranged from 10-40 to 10-25, while the sulfur fugacity ranged from 10-40 to 10-5. At the porphyry mineralization stage, the pH value of the ore deposits was about 6.9-6.1, the oxidation-reduction potential was between -1.3 and 0.3 V. At the quartz large vein mineralization stage, the oxygen fugacity was less than 10-13, the sulfur fugacity ranged from 10-14 to 10-6, and the ore deposit was in a weakly reducing environment. All these results help to understand the metallogenic environment of Lanjiagou molybdenum deposit.
Keywords: Lanjiagou molybdenum depositpH valuesulfur fugacityoxygen fugacitymineral assemblagesfluid inclusions
兰家沟钼矿位于华北克拉通北缘,辽宁省葫芦岛市境内,是大型斑岩型钼矿床,由下兰家沟、中兰家沟、上兰家沟、元宝山、西山、小马沟等钼矿床组成[1-2]。流体包裹体测试表明兰家沟矿床成矿属中-高温成矿系列,矿质在弱酸-酸性介质中、还原条件下运移、沉淀而成矿[3],兰家沟钼矿的成矿流体主要是来自混合的岩浆水与大气降水,成分主要为H2O,以及微量的CO2、CO32-[4]。高温时,钼主要以H2MoO4形式运移,在低温或含极高HCl条件下,钼能够以少量MoO2Cl2形式存在。在温压条件变化引起的流体演化过程中,钼络合物可能产生选择性运移和沉淀[5]。目前认为:兰家沟钼矿正是由于成矿金属元素以络合物的形式迁移到有利的成矿空间,并在适当物理化学条件下沉淀,而形成矿床[6]。
地球演化过程中的成矿成岩作用是非常复杂的物理化学过程,对这些过程的理解需要借助成矿成岩过程中形成矿物的物理化学条件,包括成矿成岩时的介质浓度、温度、压力、组分的活(逸)度及Eh-pH值等[7-12],从而限定成矿元素迁移沉淀的物理化学条件,阐明矿床的形成机理。因此,特定物理化学条件下形成的矿床,其共生矿物组合无疑是矿床形成和演化过程的最直接的指示标记,这些共生矿物形成时的物理化学条件的交集最有可能指示成矿的物理化学条件[13]。这使得基于矿物共生组合研究成矿热力学参数成为可能。矿床共生矿物热力学参数的相关研究不仅可以揭示成矿流体中成矿元素在迁移、沉淀过程中的物理化学条件以及矿物的演化规律,而且能够为利用矿物共生组合找矿提供依据[14-16]。
为此,以兰家沟钼矿3个成矿阶段的矿物共生组合划分为基础,测试并估算不同成矿阶段石英流体包裹体温度和压力,计算不同成矿阶段矿物平衡共生体系的热力学参数,得到Eh-pH相图和fO2-fS2相图,获得兰家沟斑岩型钼矿床成矿的物理化学条件。本文的讨论可为今后研究该类型矿床的成矿热力学条件和成矿机制起到一定的借鉴作用。
1 兰家沟钼矿地质特征及成矿阶段划分兰家沟钼矿位于华北克拉通北缘东段,燕辽沉降带的东端,华北克拉通北缘断裂(赤峰—开源断裂)以南,辽西火山岩凹陷与山海关隆起的过渡带[17-18],女儿河断裂与下兰家沟—老虎洞断裂的交汇部位[19-20]。本区在燕山期发生强烈的构造活化,主要表现为中、晚元古代至三叠纪地层发生强烈的褶皱和断裂,形成大背岭背斜、寺儿堡向斜等褶皱构造以及女儿河断裂、阜新—锦州等NE向断裂、上奈林—药王庙(南段)等NNE向断裂,白马石—碱厂等SN向断裂。兰家沟矿床岩体岩性主要为粗粒花岗岩、细粒花岗岩、石英斑岩,脉岩主要为角闪斜长煌斑岩、二长斑岩、闪长玢岩、辉长岩辉绿岩[21],矿区出露地层为中元古界蓟县系雾迷山组白云质灰岩、含燧石条带白云岩,寒武系—奥陶系灰岩、页岩和白垩系义县组火山碎屑岩不整合覆盖于其上[22]。
兰家沟钼矿的成矿阶段划分见表 1。从表 1可以看到兰家沟钼矿的成矿阶段从高到低依次为石英斑晶阶段、斑岩型矿化阶段和石英大脉型矿化阶段。这3个成矿阶段的温度覆盖区间范围都很大,并且有重叠。
Table 1
表 1 兰家沟钼矿成矿阶段划分Table 1 The metallogenic stages of Lanjiagou molybdenum deposit
| 表 1 兰家沟钼矿成矿阶段划分Table 1 The metallogenic stages of Lanjiagou molybdenum deposit |
1) 石英斑晶阶段,矿物组合为黑云母+钾长石+斜长石+石英+磁铁矿,包裹体均一温度为242~346 ℃,频数最高温度为297 ℃,获得压力为67.42 MPa。
2) 斑岩型矿化阶段,矿物组合为绢云母+石英±钾长石+辉钼矿±黄铁矿,包裹体均一温度为231~457 ℃,频数最高温度为268 ℃,获得压力为41.1 MPa。
3) 石英大脉型矿化阶段,矿物组合为绢云母+石英+辉钼矿+黄铁矿,包裹体均一温度为150~335 ℃,频数最高温度为205 ℃,获得压力为25 MPa。
图 1是兰家沟钼矿床不同成矿阶段流体包裹体的均一温度直方图。兰家沟钼矿的包裹体类型主要有:1)CO2+H2O-NaCl A型包裹体;2)CO2+H2O-NaCl B型包裹体;3)L-V+S含石盐子晶的NaCl-H2O体系包裹体;4)L-V型富液相NaCl-H2O溶液包裹体;5)富气相型包裹体;6)Ⅴ型纯气相包裹体;7)Ⅵ型纯液相水溶液包裹体。各成矿阶段热液石英中发育不同的包裹体组合。石英斑晶阶段的流体包裹体组合为类型1、2、3、4;斑岩型矿化阶段的流体包裹体组合为类型4、5,偶见类型2、3;石英大脉型矿化阶段主要以类型4包裹体为主体。为对比不同成矿阶段的包裹体温压数据,选择在3个成矿阶段出现频数都很高的类型4的包裹体温压数据作为热力学参数计算的温压条件。对于石英大脉型成矿阶段,虽然240~260 ℃间的频数很高,但考虑到220~360 ℃间的温度频数低于220 ℃以下的温度频数,因此选择频数最多的温度205 ℃作为这个成矿阶段的代表温度。
Fig. 1
Download: JPG larger image | |
图 1 不同成矿阶段流体包裹体均一温度直方图 Fig. 1 Homogeneous temperature histogram of fluid inclusions at different ore-forming stages 图 1 不同成矿阶段流体包裹体均一温度直方图 Fig. 1 Homogeneous temperature histogram of fluid inclusions at different ore-forming stages --> |
包裹体岩相学和显微测温分析在中国地质大学(北京)流体包裹体重点实验室和北京科技大学土木与环境工程学院完成。岩相学观察所用仪器为ZEISS Axioskop 40型正交偏光、反光显微镜和Axiolab型正交偏光显微镜,放大倍数为100~800倍。显微测温用冷热台为Linkam THMSG600型,并利用美国FLUIDINC公司提供的人工合成流体包裹体标准样品对冷热台进行温度标定,该冷热台可测温范围为-196~600 ℃, 冰点温度误差为±0.1 ℃,均一温度误差小于±2 ℃。流体包裹体测试过程中,升、降温的速率为(0.2~5) ℃/min,相变点附近升、降温速率均<1 ℃/min。
2 兰家沟钼矿不同成矿阶段的Eh-pH成矿流体的Eh-pH有很多种方法可以直接测量, 例如水岩体系法等,但是由于受实验条件及空气中氧等多种因素的影响,所测数据重现性差[22]。此外,很多****通过包裹体法[13]、去极化法[23]、“戴塔根”法[24]进一步探讨了测定的氧化还原电位的准确度。考虑到成岩成矿流体是非常复杂的多相组分系统,必然存在矿物共生组合,而矿物必须具有同样的热力学参数才能共生。因此可以利用兰家沟钼矿不同成矿阶段的矿物共生组合和能斯特公式大致估算不同温度下Eh-pH之间的关系并绘制相图,确定共生矿物的Eh-pH范围。
2.1 石英斑晶阶段这一成矿阶段没有辉钼矿产出, 矿物组合为黑云母+钾长石+斜长石+石英+磁铁矿,由于表生氧化作用,磁铁矿被氧化为褐铁矿。石英包裹体均一温度为242~346 ℃,频数最高的包裹体温度为297 ℃,磁铁矿的出现说明形成矿体的环境处于氧化环境。在热液系统中磁铁矿可以有如下化学过程形成:
$3{\rm{F}}{{\rm{e}}_2}{\rm{Si}}{{\rm{O}}_4} + {{\rm{O}}_2} = 2{\rm{F}}{{\rm{e}}_3}{{\rm{O}}_4} + 3{\rm{Si}}{{\rm{O}}_2}.$ | (1) |
$Eh = {E^0} + \frac{{RT}}{{nF}}\ln \frac{{\left[ {氧化态} \right]}}{{\left[ {还原态} \right]}},$ | (2) |
${E^0} = \frac{{{\Delta _r}G_m^0}}{{nF}}.$ | (3) |
Table 2
表 2 基于反应(1)利用热力学原理估算的Eh值Table 2 Estimated Eh values based on the thermodynamic principles and reaction (1)
| 表 2 基于反应(1)利用热力学原理估算的Eh值Table 2 Estimated Eh values based on the thermodynamic principles and reaction (1) |
2.2 斑岩型矿化阶段本成矿阶段的矿石矿物为辉钼矿、黄铁矿,脉石矿物为绢云母、石英、钾长石。包裹体均一温度为231~457 ℃,频数最高的包裹体温度为268 ℃,与成矿密切相关的围岩蚀变主要为硅化、绢云母化,两种蚀变主要发生于细粒斑状花岗岩和粗粒花岗岩中。这个成矿阶段主要的矿物共生平衡反应如下:
${\rm{FeS}} + 0.5{{\rm{S}}_2} = {\rm{Fe}}{{\rm{S}}_2},$ | (4) |
${\rm{Mo}}{{\rm{O}}_4} + {{\rm{S}}_2} = {\rm{Mo}}{{\rm{S}}_2} + 2{{\rm{O}}_2},$ | (5) |
$\begin{array}{l}\;\\\;\;\;\;\;\;2{\rm{CaN}}{{\rm{a}}_4}{\rm{A}}{{\rm{l}}_6}{\rm{S}}{{\rm{i}}_{14}}{{\rm{O}}_{40}} + {{\rm{K}}_2}{\rm{O}} + 2{{\rm{H}}_2}{\rm{O}} = \\2{\rm{KA}}{{\rm{l}}_3}{\rm{S}}{{\rm{i}}_3}{{\rm{O}}_{10}}{({\rm{OH}})_2} + 2{\rm{N}}{{\rm{a}}_2}{\rm{O}} + {\rm{CaO}} + 8{\rm{Si}}{{\rm{O}}_2},\end{array}$ | (6) |
$\begin{array}{l}3{\rm{KAlS}}{{\rm{i}}_3}{{\rm{O}}_8} + 2{{\rm{H}}^ + } = \\{\rm{KA}}{{\rm{l}}_2}\left[ {{\rm{AlS}}{{\rm{i}}_3}{{\rm{O}}_8}} \right]{({\rm{OH}})_2} + 6{\rm{Si}}{{\rm{O}}_2} + 2{{\rm{K}}^ + },\end{array}$ | (7) |
$\begin{array}{*{20}{l}}{{\rm{K}}{{({\rm{Mg}},{\rm{Fe}})}_3}\left[ {{\rm{AlS}}{{\rm{i}}_3}{{\rm{O}}_8}} \right]{{({\rm{OH}})}_2} + 2{\rm{AlC}}{{\rm{l}}_3} = }\\{{\rm{KA}}{{\rm{l}}_2}\left[ {{\rm{AlS}}{{\rm{i}}_3}{{\rm{O}}_8}} \right]{{({\rm{OH}})}_2} + 3({\rm{Mg}},{\rm{Fe}}){\rm{C}}{{\rm{l}}_2}.}\end{array}$ | (8) |
Fig. 2
Download: JPG larger image | |
图 2 231, 268, 457 ℃下不同体系的共生相图 Fig. 2 The phase diagrams at temperatures of 231, 268, and 457 ℃ 图 2 231, 268, 457 ℃下不同体系的共生相图 Fig. 2 The phase diagrams at temperatures of 231, 268, and 457 ℃ --> |
利用化学反应6、7、8和标准热力学参数计算得到Al-K-Si-H2O体系相图。根据图 2(b)可以观测到本成矿阶段在3个不同温度下,pH变动范围为6.9~6.1。
综合以上分析,本成矿阶段可以认为pH为6.9~6.1,而氧化还原电位为-1.3~0.3,处于弱氧化到弱还原环境条件下,并且处于弱还原环境的概率更大。
2.3 石英大脉型矿化阶段本成矿阶段的矿石矿物为辉钼矿、黄铁矿,脉石矿物为绢云母、石英。包裹体均一温度为150~335 ℃,频数最高的包裹体温度为205 ℃。由于辉钼矿比黄铁矿的稳定区将更大,只需研究黄铁矿在不同温度下与绢云母共生的实际情况,就可以得到此阶段的成矿热力学参数。
利用化学反应6、7、8和标准热力学参数计算得到Al-K-Si-H2O体系相图。从图 3可以大致推断本成矿阶段从低温到高温Eh都处于Eh=0附近,变化范围很小,大约0.5 V,pH值变动范围为6.9~5.2。整体看本成矿阶段矿床的成矿环境属于弱还原环境。
Fig. 3
Download: JPG larger image | |
图 3 不同温度条件的Al-K-Si-H2O体系 Fig. 3 Al-K-Si-H2O system at different temperatures 图 3 不同温度条件的Al-K-Si-H2O体系 Fig. 3 Al-K-Si-H2O system at different temperatures --> |
3 兰家沟钼矿不同成矿阶段的fO2-fS2氧逸度和硫逸度对于矿床的形成具有非常重要的意义,氧逸度在一定程度上控制了岩浆与热液成矿过程中成矿元素的活化迁移和富集成矿[27],硫在岩浆中的溶解度与岩浆熔体的Fe离子含量有关[28]。
3.1 石英斑晶阶段以化学反应(1)为例,本成矿阶段的氧逸度计算公式[29]如下:
${K_{\rm{P}}} = \frac{{\alpha _{{\rm{F}}{{\rm{e}}_3}{0_4}}^2 \cdot {\alpha _{{{\rm{S}}_{\rm{i}}}{0_2}}}}}{{\alpha _{{\rm{F}}{{\rm{e}}_2}{\rm{s}}{{\rm{i}}_4}}^3 \cdot {f_{{0_2}}}}},$ | (9) |
$\Delta G_{反应}^0 = - RT\ln {K_p},$ | (10) |
Fig. 4
Download: JPG larger image | |
图 4 Fe-O-S体系中磁铁矿氧逸度与硫逸度综合图 Fig. 4 Oxygen fugacity and sulfur fugacity of magnetite in Fe-O-S system 图 4 Fe-O-S体系中磁铁矿氧逸度与硫逸度综合图 Fig. 4 Oxygen fugacity and sulfur fugacity of magnetite in Fe-O-S system --> |
3.2 斑岩型矿化阶段本成矿阶段的矿石矿物为辉钼矿、黄铁矿,可以利用黄铁矿、辉钼矿的矿物共生关系,估计本成矿阶段的硫逸度,计算结果见图 5。结合图 5可以看出黄铁矿和辉钼矿的硫逸度共生变化范围广,本成矿阶段的硫逸度范围应大于10-15,同前一阶段相比,硫逸度下限变大。
Fig. 5
Download: JPG larger image | |
图 5 斑岩型矿化阶段辉钼矿和黄铁矿的氧逸度与硫逸度综合图 Fig. 5 Oxygen fugacity and sulfur fugacity of molybdenite and pyrite in the prophyry mineralization stage 图 5 斑岩型矿化阶段辉钼矿和黄铁矿的氧逸度与硫逸度综合图 Fig. 5 Oxygen fugacity and sulfur fugacity of molybdenite and pyrite in the prophyry mineralization stage --> |
综上,大致可以确定本成矿阶段的氧逸度应小于10-17, 硫逸度为10-15~10-3。
3.3 石英大脉型矿化阶段本成矿阶段的矿石矿物为辉钼矿、黄铁矿,图 6(a)、6(b)分别是辉钼矿和黄铁矿在3个温度下的氧逸度、硫逸度综合图。与上一阶段相比,硫逸度变化范围变小,在150 ℃条件下,硫逸度下限为10-14,随着温度增高,硫逸度逐渐增加。同前一成矿阶段相比,氧逸度变化不大。综上,根据共生的辉钼矿、黄铁矿氧逸度和硫逸度变化图,大致可以确定本成矿阶段的氧逸度应小于10-13,硫逸度变化的区间为10-14~10-6。
Fig. 6
Download: JPG larger image | |
图 6 石英大脉型矿化阶段辉钼矿和黄铁矿的氧逸度与硫逸度综合图 Fig. 6 Oxygen fugacity and sulfur fugacity of molybdenite and pyrite in the quartz large vein mineralization stage 图 6 石英大脉型矿化阶段辉钼矿和黄铁矿的氧逸度与硫逸度综合图 Fig. 6 Oxygen fugacity and sulfur fugacity of molybdenite and pyrite in the quartz large vein mineralization stage --> |
4 结论与讨论本文基于兰家沟钼矿3个成矿阶段的矿物共生组合,通过不同成矿阶段形成的石英包裹体估算不同成矿阶段的温度与压力,进而估算兰家沟矿床的成矿热力学参数,计算结果与已有研究对比,基本一致。估算的3个成矿阶段的热力学参数如表 3所示。
Table 3
表 3 兰家沟钼矿不同成矿阶段的热力学参数Table 3 Thermodynamic parameters at different metallogenic stages of Lanjiagou molybdenum deposit
| 表 3 兰家沟钼矿不同成矿阶段的热力学参数Table 3 Thermodynamic parameters at different metallogenic stages of Lanjiagou molybdenum deposit |
石英斑晶阶段形成矿体的环境为氧化环境,氧逸度的区间范围为10-40~10-25,硫逸度的区间为10-40~10-5;斑岩型矿化阶段,其形成环境的pH在6.9~6.1之间,氧化还原电位在-1.3~0.3之间,氧逸度小于10-17,硫逸度为10 -15~10-3,成矿环境为弱还原环境;石英大脉型矿化阶段的氧逸度小于10-13,硫逸度为10-14~10-6,成矿环境为弱还原环境。
矿物共生组合是成矿成岩地质过程的必然结果,因此基于矿物共生组合估算成矿热力学条件应能准确获得成矿成岩的热力学条件,也为定量估计成矿热力学条件提供了可行的方法。这为查明多金属矿床的成矿物理化学性质、演化、成矿模式和勘查模型以及矿区的深部和外围潜在资源评价和找矿提供依据。但是使用时要注意:1)尽可能准确界定不同成矿阶段的矿物共生组合;2)尽可能采用同一来源的矿物热力学基础数据。
参考文献
[1] | 李友权. 兰家沟钼矿田成因类型及找矿标志[J]. 中国科技博览, 2012(20): 244-244. |
[2] | 马永昌, 王长刚, 冯国清, 等. 杨家杖子矿区再找矿[J]. 矿床地质, 2002, 21(S1): 435-438. |
[3] | 代连铎, 张玉平, 宋雨春. 兰家沟钼矿床地质特征及成矿模式[J]. 有色矿冶, 2008, 24(6): 4-7. Doi:10.3969/j.issn.1007-967X.2008.06.003 |
[4] | 代军治, 毛景文, 谢桂青, 等. 辽西兰家沟钼矿床成矿流体特征及成因探讨[J]. 矿床地质, 2007, 26(4): 443-454. Doi:10.3969/j.issn.0258-7106.2007.04.008 |
[5] | 孙燕, 刘建明, 曾庆栋. 斑岩型铜(钼)矿床和斑岩型钼(铜)矿床的形成机制探讨:流体演化及构造背景的影响[J]. 地学前缘, 2012, 19(6): 179-193. |
[6] | Sangster D. Mississippi Valley-type lead-zinc[J]. Geology of Canadian Mineral Deposit Types Edited by OR Eckstrand, WD Sinclair and RI Thorpe Geological Survey of Canada, 1996, 78(8): 253-261. |
[7] | 吴湘滨, 戴塔根. 湘西南石英脉型金矿矿物流体包裹体pH值和Eh值的估算[J]. 大地构造与成矿学, 1999, 23(3): 248-253. Doi:10.3969/j.issn.1001-1552.1999.03.007 |
[8] | 梁鹤.华南白垩纪岩背斑岩锡矿成矿斑岩的成因及其成矿意义[D].北京: 中国科学院大学, 2017. |
[9] | 邓金灿, 李蘅, 徐文炘, 等. 广西巴里-龙头山100号锡多金属矿床成矿物理化学环境研究:氧逸度(fO2)、CO2逸度(fCO2)、硫逸度(fS2)、mΣs[J]. 矿产与地质, 2001, 15(3): 167-171. Doi:10.3969/j.issn.1001-5663.2001.03.004 |
[10] | 张聚全, 李胜荣, 卢静. 中酸性侵入岩的氧逸度计算[J]. 矿物学报, 2018, 38(1): 1-14. |
[11] | 张荣伟, 彭建堂, 邓起东, 等. 滇东北茂租热液型铅锌矿床矿物组合共生分异的热力学Eh-pH相图[J]. 地质论评, 2017, 63(5): 1 391-1 400. |
[12] | 弋嘉喜.察尔汗盐湖矿物组合特征及其成因指示[D].北京: 中国科学院大学, 2017. http://www.cnki.com.cn/Article/CJFDTOTAL-YHYJ201702007.htm |
[13] | 艾永富, 马瑞志. 下兰家沟钼矿中菱锰矿-方解石共生组分特征及其对钼矿化的指示意义[J]. 地质与勘探, 1983(2): 32-36. |
[14] | Worley B, Powell R. Making movies:phase diagrams changing in pressure, temperature, composition and time[J]. Geological Society, London, Special Publications, 1998, 138(1): 269-280. Doi:10.1144/GSL.SP.1996.138.01.15 |
[15] | 刘斌. 利用流体包裹体及其主矿物共生平衡的热力学方程计算形成温度和压力[J]. 中国科学(B辑), 1987(3): 81-88. |
[16] | 尹淑苹, 李永兵, 田会全, 等. 黑龙江乌拉嘎金矿成矿热力学参数研究[J]. 岩石矿物学杂志, 2014, 33(4): 720-725. Doi:10.3969/j.issn.1000-6524.2014.04.010 |
[17] | 申永治. 兰家沟钼矿床同位素年代学研究[J]. 辽宁地质, 1989(4): 312-320. |
[18] | Zeng Q, Liu J, Qin K, et al. Types, characteristics, and time-space distribution of molybdenum deposits in China[J]. International Geology Review, 2013, 55(11): 1 311-1 358. Doi:10.1080/00206814.2013.774195 |
[19] | 田豫才. 辽西兰家沟钼矿区成矿构造、岩浆演化及成矿作用[J]. 矿产与地质, 1999, 13(3): 135-140. Doi:10.3969/j.issn.1001-5663.1999.03.002 |
[20] | 梁正宇. 兰家沟含钼花岗岩岩石化学及其应用[J]. 有色矿冶, 2011, 27(1): 1-4. Doi:10.3969/j.issn.1007-967X.2011.01.001 |
[21] | 代军治, 毛景文, 赵财胜, 等. 辽西兰家沟钼矿床花岗岩SHRIMP锆石U-Pb年龄及岩石化学特征[J]. 地质学报, 2008, 82(11): 1 555-1 564. |
[22] | 刘金辉, 孙占学. 确定砂岩型铀矿体定位新方法:水岩体系Eh-pH法[J]. 吉林大学学报(地球科学版), 2004, 34(1): 44-48. |
[23] | 王爱军, 石超英. 去极化方法测定海水的氧化还原电位初探[J]. 海洋技术学报, 2010, 29(2): 120-122. |
[24] | 戴塔根. 地球化学样品的Eh-pH值测定及应用[J]. 湖南地质, 1992(S1): 72-74. |
[25] | 傅献彩, 沈文霞, 姚天扬. 物理化学[M]. 北京: 高等教育出版社, 1990. |
[26] | 林传仙, 白正华, 张哲儒. 矿物及有关化合物热力学数据手册[M]. 北京: 科学出版社, 1985. |
[27] | Mavrogenes J A, O'Neill HSC. The relative effects of pressure, temperature and oxygen fugacity on the solubility of sulfide in mafic magmas[J]. Geochimica et Cosmochimica Acta, 1999, 63(7/8): 1 173-1 180. |
[28] | 刘善琪, 李永兵, 曾庆栋, 等. 利用矿物共生关系确定成矿条件:以海沟金矿为例[J]. 岩石矿物学杂志, 2013, 32(1): 99-105. Doi:10.3969/j.issn.1000-6524.2013.01.008 |
[29] | 印永嘉, 奚正楷, 张树永. 物理化学简明教程[M]. 4版. 北京: 教育出版社, 2007. |