斑岩铜矿蚀变分带与成矿机制：玉龙矿床水-岩反应热力学平衡模拟例析
张少颖1,2，和文言2，高雪2，田成华2，肖仪武1
1 矿冶科技集团有限公司
2 中国地质大学(北京)地质过程与矿产资源国家重点实验室，深时数字地球前沿科学中心
作者简介：张少颖，博士，矿物学、岩石学、矿床学专业。
通讯作者：和文言，副教授，主要从事矿床学与矿产勘查学工作。
斑岩铜矿具有蚀变分带特征，准确厘定蚀变分带，是分析成矿机制、建立斑岩成矿模型和找矿预测的关键，对指导找矿具有积极意义。
Porphyry copper deposits exhibit characteristics of alteration zoning. Accurately defining the alteration zoning is crucial for analyzing the mineralization mechanism, establishing a porphyry mineralization model, and predicting ore deposits, which has positive significance in guiding ore prospecting.
通常认为青磐岩化和绢英岩化蚀变的形成是由岩浆热液与地下水混合形成的。然而，也有研究强调其形成主要与岩浆热液温度下降有关，而无需地下水的混合，本文研究亦如此。即仅仅依靠温度的下降，单一岩浆流体与斑岩体的水岩反应过程即可产生在斑岩系统中观察到的所有蚀变和矿化特征。
本文以金沙江-哀牢山成矿带玉龙斑岩铜矿为例，引入热力学平衡的基本概念和两种常用的算法模型，通过热力学平衡模拟定量化岩浆热液矿床蚀变和矿化过程，开展了蚀变分带成因和金属成矿机制研究。结果表明，与钾硅酸盐化蚀变相关的流体具有较高pH值(5.0~7.0)和氧逸度(ΔFMQ=+2.9~+3.6)特征；而与青磐岩化阶段有关的流体pH值(5.0~6.4)和氧逸度(ΔFMQ=+1.1)均有所下降；与绢英岩化阶段有关的流体pH值(5.0~5.7)进一步小幅下降，而氧逸度则(ΔFMQ=+1.7)略有回升。
文中通过利用LMA反演及GEM正演相结合定量化地刻画了玉龙斑岩铜矿水岩反应过程，由此揭示了斑岩矿床蚀变分带是逐渐冷却的单一岩浆热液与斑岩体不断反应的结果，且不同温度窗口对应着钾硅酸盐化(450~360℃)、青磐岩化(360~320℃)和绢英岩化(320~200℃)蚀变矿物的形成，故含矿流体温度的快速下降可能是玉龙铜矿蚀变叠加的重要因素。此外，Cu、Mo络合离子溶解度对温度变化的差异响应，导致了Mo矿化主要发育于靠近斑岩体的高温区域，而Cu则以网脉状-浸染状叠加到Mo矿化之上，并广泛分布于斑岩体周边的高-中温区域。
基金项目：本文受国家自然科学基金项目(42302082、42372098)、国家重点研发计划项目(2022YFF0800902)、高等学校学科创新引智计划(BP0719021)、中国地质大学地质过程与矿产资源国家重点实验室专项基金(MSFGPMR201804)和深时数字地球前沿科学中心“深时数字地球”中央高校科技领军人才团队项目(2652023001)联合资助
1 热力学平衡的概念及计算方法
1.1 质量作用定律
1.2 吉布斯能最小化
2 玉龙铜矿地质特征
3 玉龙铜矿水岩反应热力学平衡模拟
3.1 成矿初始流体物化条件模拟
3.2 GEM正演流体演化轨迹
3.3 LMA反演流体演化轨迹
4 玉龙铜矿蚀变分带与金属成矿
4.1 流体演化轨迹与蚀变分带成因
4.2 Cu、Mo沉淀机制
5 结论
斑岩型矿床中所具有的相对一致的蚀变分带和脉体共生及穿插关系，如钾硅酸盐化、青磐岩化及绢英岩化等同心环状蚀变分带和与之对应的A、B、D脉体及其特征矿物组合，可能反映了岩浆-热液流体演化中一个重要的共性过程。这一过程通常被归因于岩浆持续结晶所控制的流体出溶性质的变化，这些流体的组成进而决定了由流体引起的热液蚀变特征。然而，也有学者提出，斑岩型矿床特征的蚀变分带和脉体类型是由于温度下降而非岩浆流体的周期性脉冲所造成的出溶流体组成变化，且大气降水的参与对于斑岩矿床蚀变分带的贡献也存在争议。在这种背景下，铜、钼等金属多被认为主要沉淀于早期的钾硅酸盐化阶段。然而，一些研究表明，许多斑岩矿床矿化却主要位于青磐岩化甚至绢英岩化蚀变带中。同时，辉钼矿和黄铜矿等矿物通常赋存于不同脉体中，暗示Cu和Mo元素的分离和沉淀可能是不同期次发生的。这促使Mo和Cu元素在斑岩型矿床中通常表现出不同的空间分布模式，即Mo富集在靠近岩体的中心矿壳中，而Cu富集在外矿壳中。对于上述斑岩矿床蚀变分带与金属成矿关系和Cu-Mo空间分离现象引起的争议，通常可通过流体包裹体和水岩反应实验等进行验证。然而，由于成矿流体在地质演化过程中难以有效保存，导致部分研究结果相互矛盾且难以解释。同时，复杂的水岩反应实验研究面临诸多方法和技术上的问题，例如开放系统的构建和矿物缓冲对的选择。此外，将实验研究在矿物尺度上观察的微观现象与热液矿床在岩石尺度上观察的宏观现象相匹配也是一个挑战。
近年来，计算机模拟技术的成熟为热液矿床成矿过程研究提供了新的可能性。特别是，结合了流体-矿物平衡热力学动态传输模型的建立，对于理解和预测热液矿床的形成过程起到了关键作用。以斑岩矿床为例，热力学模拟探讨了岩浆-热液流体性质的瞬时性，以及由压力波动引起的相分离、高盐度流体包裹体和富气相流体包裹体共存等现象的成因机制。利用金属的实验气/液分配系数以及火山气/熔体分配系数，一些热力学模型还探讨了气体运移、沉淀金属的能力，以及大气水的混合对多相混溶H2O-NaCl流体中金属沉淀的影响。以上研究显示，热力学模拟为衡量热液系统中流体相平衡和热力学性质(PVTX)提供了新的思路和解决方案，对分析成矿流体演化和成矿作用发挥着重要作用。
我国斑岩型矿床受独特的构造演化历程所导致的更快的抬升速率或复杂的构造叠加影响，蚀变叠加现象往往较为发育，这一特点使得我国斑岩型矿床的蚀变分带与金属成矿机制变得更加复杂。其中，玉龙矿床是青藏高原东缘金沙江-哀牢山成矿带中最大的斑岩型铜(钼)矿床。与后碰撞型斑岩铜矿的蚀变分带特征相似，玉龙铜矿晚期的绢英岩化强烈叠加在早期的钾硅酸盐化带之上，这模糊了蚀变分带特征及其与矿化之间的关系，阻碍了对成矿过程的精细刻画及其成因机制的深入理解。目前，随着对玉龙铜矿蚀变分带及矿物组合特征、流体包裹体、矿物-全岩地球化学及同位素的深入研究，为构建水岩反应热力学模型、厘定流体演化过程及金属成矿机制奠定了良好的基础。为此，本文引入热力学平衡的基本概念和两种常用的算法模型，并以玉龙铜矿为例，通过热力学平衡模拟定量化岩浆热液矿床蚀变和矿化过程，明确斑岩型矿床蚀变分带成因及金属成矿机制。
任意化学系统均应向总吉布斯自由能最小的状态发展，即热力学平衡状态。因此，寻找系统化学平衡状态即在保证系统吉布斯自由能最小且满足系统平衡约束条件下获得系统中各物质的含量。然而，水岩反应总是涉及多相-多组分的过程，这就需要采用高效且鲁棒性强的计算方法来处理复杂化学反应过程。Smith and Missen(1982)最开始将化学平衡算法划分为化学计量数和非化学计量数两类方法。其中，前者主要利用质量平衡和质量作用方程(LMA)，而后者则通过系统吉布斯自由能最小化(GEM)来约束平衡条件。然而，在Van Zeggeren and Storey(2011)看来，LMA和GEM的化学平衡算法是互补的，二者之间并无本质区别。两种计算方法的特征详见表1。
表1 LMA与GEM平衡计算方法比较
质量作用定律(LMA)算法通过求解一个结合质量平衡和质量作用表达式的非线性方程系统来确定化学反应平衡常数K。
其中，A、B、C、D表示化学物质，a、b、c、d为化学计量数。当该反应达到化学平衡时，可知：
式(3)括号中内容为生成物活度的化学计量数次方乘积除去反应物活度的化学计量数次方乘积，通常定义为Q，故上式一般表现形式为：
其中，Δrμo指的是纯的、处于标准状态的生成物与反应物之间的吉布斯自由能差值。Δrμ则指的是实际反应中达到或未达到平衡状态的生成物与反应物之间的自由能差值。Q代表的活度(如物质的浓度、逸度以及摩尔分数等)会随着反应的进行而不断变化，进而导致Δrμ的不断变化。因此，Δrμo对于一个可逆反应来说是一个确定值，不随反应物与生成物的活度而变，Δrμ和Q则会随着反应物与生成物的活度随时变化。
K称为热力学平衡常数，即平衡时生成物活度与反应物活度的化学计量数之比，不随系统物质组成而变化。任意温度、压力下的K值可由相关计算软件获得，之后即可计算平衡状态下各物质的活度。目前，已开发的用于计算K值的软件包括WATEQ、MINTEQA2、Eq.3/6、CHESS、PHREEQC，以及Supcrt92、SupcrtBL和Geochemist's Workbench等。一般来说，这些基于化学计量数方法的模型只考虑化学反应的平衡条件，而无法获得系统的吉布斯自由能。此外，Leal et al.(2013)提出上述地球化学计算模型也存在着一些数学方法上的不连续性，即系统中主要物质含量是通过质量平衡方程确定的，而次要物质含量是通过利用质量作用方程的连续替代方法计算的，这种数学方法的不连续性将导致一系列的算法效率及函数收敛问题。
吉布斯能最小化模型(GEM)通过指定化学元素的总量以及电荷平衡来确保系统质量平衡，通过最小化系统的吉布斯自由能来寻找平衡时物质相组合，即：
其中，质量平衡方程An=b确保了封闭系统下元素守恒，不等式n≥0则保证了物质的量在平衡时是非负的。基于吉布斯能最小化模型(GEM)的地球化学建模软件包括Chem Sage、THERIAC)、HCh、Fact Sage，以及GEM-Selektor和Reaktoro等。早期的软件包如Chem Sage首次构建了基于物质非理想性化学势的吉布斯能最小化算法，然而其导数却对应于理想模型。为了解决这种做法带来的函数收敛问题，Eriksson and Hack(1990)通过人为确定初始值来提高算法效率。然而，这种改进并没有改变该方法的收敛性问题。其后，GEM-Selektor软件包既未采用对数障碍函数，也未采用常规的KKT扰动方法，而是采用了内插点算法来确定多相系统的最小吉布斯自由能。在每次迭代中，类似于线搜索策略的内插点算法寻找一个可以使系统吉布斯能下降的方向，同时满足质量平衡限制。在此基础上，基于C++/Python编译平台的地球化学建模软件Reaktoro还考虑到可能的非线性平衡约束，如pH值、活度或逸度等。
玉龙铜矿位于青藏高原东缘，是新生代金沙江-哀牢山成矿带的组成部分(图1a)。目前，玉龙铜矿已探明的铜和钼金属资源量分别为650万t(0.62%)和41万t(0.042%)，是该成矿带内最大的后碰撞斑岩矿床。矿床以始新世复式斑岩体为中心，呈岩株状侵位于恒星错背斜轴部的由三叠纪海相碎屑岩和碳酸盐岩层序组成的上三叠统地层中(图1b)。该复式斑岩体以二长花岗斑岩(MGP)为主，钻孔中可见与二长花岗斑岩呈侵入接触的花岗斑岩及钠长斑岩脉。矿区地表发育大面积粘土化蚀变，地层中多可见褐铁矿化或孔雀石化。
图1 玉龙斑岩铜矿大地构造位置图(a)、地质简图(b)及剖面图(c)(b,c)
矿区MGP蚀变强烈，从中心向外依次发育有钾化蚀变带、钾硅酸盐化蚀变带(钾化及黑云母化)、绢英岩化蚀变带及浅部的粘土化蚀变带(图1c)。钾化蚀变带位于岩体深部，以发育大量钾长石为主要特征。钾硅酸盐化蚀变带靠近岩体中心，可见钾长石+黑云母+石英矿物组合，斜长石常见有钾长石次生加大边(图2a)，伴随着石英+金属硫化物脉体的分布(图2b)，及部分钾长石斑晶呈弥漫状分布在斑岩体中(图2c)。青磐岩化蚀变带在矿区不甚发育，仅可见部分黑云母发生绿泥石化。绢英岩化在矿区显著发育，包括基质中呈放射鳞片状的绢云母集合体，以及长石斑晶表面强烈的绢云母化(图2d)。粘土化蚀变在矿区较为广泛，常与绢英岩化共同叠加在早期的钾硅酸盐化带内，X射线衍射分析和短波红外光谱表明，其主要由高岭石、蒙脱石和伊利石组成。
图2 玉龙铜矿典型蚀变类型特征
(a)斜长石边部的钾长石次生加大边；(b)钾硅酸盐化脉体穿切斑岩体；(c)弥漫状钾长石化；(d)斜长石斑晶及基质的绢云母化
钻孔ZK0908编录结果显示，赋矿二长花岗斑岩体蚀变类型有钾长石化、黑云母化、绿泥+绿帘石化，以及浅部发育的粘土化(图3)。其中，以钾长石和黑云母为代表的钾硅酸盐化蚀变主要位于钻孔300m以下深部，而以绢云母为代表的绢英岩化则全孔均较为发育。此外，以绿泥石+绿帘石为代表的青磐岩化蚀变主要叠加于黑云母化蚀变带之上，但仅局部可见。该钻孔400m以上Cu矿化较好，且由于浅部粘土化蚀变的强烈发育，这促进了黄铜矿等金属硫化物的次生富集，使得该钻孔浅部Cu品位更高。
目前，有关于玉龙流体包裹体的研究业已约束其部分物理化学性质，如初始单相中等密度流体包裹体的均一温度(400~550℃)、盐度((9.1±3.0)%NaCleqv)和主要金属含量(~1000×10-6Cu和(22±19)×10-6Mo)。然而，某些关键参数，如流体pH值、氧逸度和挥发分含量等却难以获得其定量结果。经典热液成矿模型显示，含矿流体受温度、压力等强度性质变化及水岩反应过程影响，溶液金属络合物溶解度逐渐下降而导致金属逐渐沉淀。因此，为获得流体出溶时其物化参数，可逆向化处理上述矿物沉淀过程，即利用LMA方法构建含挥发份流体与硅酸盐矿物和金属硫化物三端元之间的相平衡过程，获得的平衡后流体即可代表该状态下的流体组成。模拟使用的矿物和离子的热力学数据来源于MINES2023数据库，并导入到GEM-Selektor软件包进行数据库管理。该数据库包括Holland and Powell(1998)的矿物热力学数据、SUPCRT92中的溶液离子热力学数据、含Al溶液离子的最新实验数据，并且包含使用GEMSFITS软件优化的含K、Na和Si溶液离子的具有内部一致性的热力学数据集。
首先，参考玉龙流体包裹体结果，本文将模拟温度和压力设置为450℃和0.1GPa。其次，为了在模拟过程中能够限制流体的主要成分浓度范围，我们搜集了斑岩矿床中富气相和中等密度流体包裹体的组成信息(表2)。二者S(SO2和H2S)、Cl(HCl和NaCl)和CO2含量相似，分别为0.1%~2.0%S、0.6%~9.6%Cl和1.1%~5.3%CO2。虽然目前缺乏直接测量流体包裹体中对氧化还原敏感的离子浓度的方法，但考虑到斑岩矿床深部贫矿的网脉状区域通常可见有赤铁矿+磁铁矿+硬石膏+金红石等表征高氧逸度的矿物组合，成矿初始流体应处于较高氧逸度状态，如ΔFMQ=+2。在450℃条件下，获得的系统氧逸度为logfO2=-23.8。此外，由于流体中挥发份S的组成及价态较为复杂，为利用LMA方法模拟溶液中S浓度的可行性，本文将溶液中的S简化为H2S和SO2两种缓冲端元，其余S价态的化合物均可由上述两种端元反应而来，即H2Saq+3O2,aq=2SO2,aq+2H2O。本文利用MINES2023热力学数据库及Reaktoro中的Thermofun模块通过LMA方法计算的上述矿物反应式平衡常数为logK=72.6，由此获得的流体中SO2/H2S浓度比例为3.5。结合典型斑岩矿床中流体包裹体S的浓度范围，本文设置流体中含3.5%SO2和1.0%H2S。已有资料显示，玉龙铜矿初始单相中等密度流体包裹体盐度范围为3.5%~9.1%NaCleqv，但对于流体中HCl和NaCl比例关系尚未确定。考虑到斑岩矿床中钾化阶段矿物组合包括钾长石、硬石膏、磁铁矿等，其限定的流体pH值应为近中性。由此获得的流体中电离的H+离子浓度应为10-7mol/L，进而限定的溶液中NaCl与HCl浓度之比为8左右。结合玉龙铜矿流体包裹体等效盐度范围，本文设置流体中含8.0%NaCl和1.0%HCl。另外，部分斑岩矿床成矿流体还可见有CO2的存在，且其存在可能抑制了主要含矿络合物(CuHSaq、CuCl2-、Au(HS)2-、FeCl42-)的稳定性而促进了Cu、Au、Fe的沉淀。为此，模型流体设置CO2含量为固定值2.0%。综上所述，本文构建了一个相应的含挥发份流体、钾化阶段矿物组合和金属硫化物三端元相平衡模拟：(1)1kg含8.0%NaCl、1.0%HCl、2.0%CO2、3.5%SO2和1.0%H2S的水溶液；(2)代表钾化阶段的矿物组合：钾长石KAlSi3O8、钠长石NaAlSi3O8、黑云母K(Mg，Fe)3AlSi3O10(OH)2、石英SiO2、磁铁矿Fe3O4、硬石膏CaSO4、金红石TiO2，并保证其过量；(3)代表金属来源的矿物组合：黄铜矿CuFeS2、辉钼矿MoS2，并保证其过量。
表2典型斑岩矿床中富气相和中等密度流体包裹体的组成信息
模拟结果显示，平衡后的流体pH值为4.7，logfO2=-23.0(ΔFMQ=+2.7)，为一酸性强氧化性流体。此时，该溶液金属饱和时溶解Cu含量为1138×10-6，Mo为1.2×10-6。模拟的初始流体Cu浓度接近LA-ICP-MS测试的玉龙初始单相中等密度流体包裹体中Cu浓度(~1000×10-6)，而Mo浓度也处于该类包裹体中Mo浓度范围内((20±19)×10-6)。
为了解水岩反应过程对于玉龙斑岩铜矿蚀变与矿化的影响，我们利用GEM-Selektor软件包和MINES2023热力学数据库建立了基于GEM算法的水岩反应动态传输模型。模型以3.1节获得的平衡流体组成为水岩反应过程的流体端元，并选择赋矿二长花岗斑岩体(MGP)的全岩地球化学数据为岩石端元，具体见表3。
动态传输模型共设计25个岩石单元，含矿流体从下至上依次流经该岩石序列。由于温度下降和水岩反应的影响，矿物会不断从流体中沉淀。具体模型设置如下：(1)每个单元中含有100g二长花岗斑岩。流体从岩石单元1开始与之反应并达到平衡状态，反应后的流体继续流经岩石单元2、3、4等，并与之反应后持续与其余岩石单元反应；(2)初始流体组成为含挥发份流体、钾化阶段矿物组合和金属硫化物三端元相平衡后的流体组成(节3.1)，质量设置为1000g；(3)模拟温度设置为温度为450℃到200℃区间范围，且随流体与岩石单元平衡后阶梯式下降，即每个岩石单元降温10℃；(4)由于压力变化对于斑岩成矿系统金属元素溶解度影响较小，故本文模拟暂不考虑压力变化情况。模拟压力依据斑岩体侵位深度，将其固定为0.1GPa。
模拟结果(图4)显示，在450~410℃温度范围内，当初始酸性及强氧化性流体流入MGP时，受控于初始阶段高水岩比(10/1)的影响，流体的pH值和氧逸度均受流体控制从而表现出酸性(pH=~4.8)和高氧逸度(ΔFMQ=+2.9~+2.5)特征。当温度处于410~360℃范围内，流体pH值从4.8升高到5.1附近，而流体的氧逸度ΔFMQ则逐渐从+2.5降低到+1.1，伴随着钾长石、钠长石、硬石膏、黑云母和赤铁矿、磁铁矿等典型钾硅酸盐化蚀变矿物的沉淀。当温度在360~320℃范围内，代表青磐岩化蚀变阶段的典型矿物如绿帘石、铁绿泥石和斜绿泥石等逐渐形成，而流体的pH值从5.0缓慢下降为4.8附近，流体氧逸度ΔFMQ从+0.1小幅回升到+1.1。当温度进一步从320℃下降到200℃时，流体pH值从4.8缓慢升高为5.4左右，氧逸度ΔFMQ也持续升高。在此期间，绢云母和方解石等开始沉淀并形成典型的绢英岩化蚀变。
图4 玉龙铜矿水岩反应热力学模拟
(a)水岩反应过程示意图；(b)流体pH值变化；(c)流体氧逸度ΔFMQ值变化；(d)流体中溶解的SO2/H2S比值；(e)钾硅酸盐化蚀变过程中蚀变矿物组合；(f)青磐岩化蚀变过程中蚀变矿物组合；(g)绢英岩化蚀变过程中蚀变矿物组合；(h)水岩反应过程中流体溶解金属离子浓度变化
温度的下降和水岩反应过程均降低了Cu、Mo络合离子在溶液中的溶解度，从而导致了相关金属的沉淀，但上述过程对各络合离子稳定性的影响不一。以HMoO4-和MoO42-为载体的Mo元素主要沉淀在高温区间如450~370℃范围内，Mo在溶液中的溶解度从1.2×10-6下降到0.01×10-6。主要以CuCl(CuCl43-、CuCl2-、CuCl)、少量以CuS(CuHS、Cu(HS)2-)为载体的Cu元素沉淀也受到了温度的控制作用。在450~370℃范围内，Cu的溶解度从1041.9×10-6下降到4.7×10-6，且随着温度下降其持续降低到1.0×10-6左右。然而，在青磐岩化阶段，由于Cu、Mo溶解度对温度的响应差异，溶解的Cu/Mo离子比值急剧增大，这导致了该阶段以Cu为主的金属沉淀。此外，由于斑岩型矿化为多阶段流体叠加的综合反映，故本文分别模拟了1、10和50轮流体与MGP反应后的金属沉淀含量。结果显示，当1轮流体与MGP反应时，黄铜矿和辉钼矿分别沉淀在450~300℃和450~330℃区间内，其沉淀的质量分别为8.9×10-4kg和1.1×10-6kg(图5)。当与MGP反应的流体轮数增加为10和50时，黄铜矿在450~300℃范围内沉淀的质量增长5~10倍，且在300~200℃范围内也开始有少量黄铜矿沉淀。另外，随着补充流体从10轮升高为50轮，辉钼矿在450~330℃内的沉淀量进一步加大，其沉淀质量平均增加3~5倍；而在300~200℃内涨幅较小，仅增长1倍有余。
图5 水岩反应过程中金属沉淀量
(a)水岩反应过程中黄铜矿沉淀质量；(b)水岩反应过程中辉钼矿沉淀质量.图中10和50分别代表有10和50轮含矿流体与二长花岗斑岩反应
为验证GEM模拟水岩反应过程下的成矿流体演化过程，本文根据玉龙铜矿蚀变矿物共生关系并依托LMA方法构建了相应的蚀变矿物pH-fO2相图(图6)。综合矿物温压计及流体包裹体研究结果，本文分别将玉龙铜矿的钾硅酸盐化、青磐岩化及绢英岩化蚀变阶段温度和压力限定在400℃、300℃、300℃和0.1GPa条件下。各温度、压力条件下的离子浓度参考GEM模型中流体离子浓度。黑云母固溶体中铁云母与金云母的比例为3：7左右，而绿泥石固溶体中斜绿泥石与铁绿泥石的比例则根据GEM模型结果确定为6：1.2(表4)。在黑云母蚀变为绿泥石过程中，GEM模型显示的流体Mg2+、K+离子分别为0.00021mol/kg、0.1067mol/kg。在绿泥石蚀变为绢云母过程中，GEM模型显示流体Mg2+、K+离子的浓度分别为0.006mol/kg、0.08mol/kg。
图中矿物平衡反应线主要依据表4反应平衡常数logK计算而来;圈1代表钾硅酸盐化蚀变矿物组合限定的pH和fO2范围，圈2代表青磐岩化蚀变，圈3代表绢英岩化蚀变
为保证热力学参数的一致性，本文利用MINES2023数据库，并导入到Reaktoro软件包进行数据库管理，进而利用其中的ThermoFun模块计算获得表4中各矿物反应式的平衡常数logK值，最终依托该值获得蚀变矿物pH-fO2相图。
模拟结果显示，受钾长石、黑云母、硬石膏及磁铁矿和赤铁矿共存区限制，与钾硅酸盐化蚀变相关的流体pH值为7左右，而氧逸度ΔFMQ为+3.6左右，属于强氧化、偏中性流体(图6a点1)。当温度从400℃降低为300℃时，斑铜矿+黄铁矿的稳定相区将大幅扩大，故钾硅酸盐化阶段的黑云母+黄铜矿+斑铜矿脉体的出现应是温度下降导致的而非氧逸度变化导致的(图6b点2)，这与前人实验结果一致。当温度下降到300℃附近时，由黑云母绿泥石化、角闪石绿帘石化等典型青磐岩化蚀变矿物限定的流体pH值为6.4，而受钾硅酸盐化阶段大量赤铁矿和磁铁矿结晶影响，流体氧逸度略有下降。在绢英岩化阶段，伴随着绢云母+方解石组合的出现，流体的pH值进下降至绢云母和方解石稳定域，即pH值位于5.2~5.7之间。此外，受黄铁矿沉淀导致的溶液还原性S浓度下降影响，与绢英岩化阶段有关的流体氧逸度略有回升(图6b点3)。
综上所述，玉龙斑岩铜矿的流体从早期的高温钾硅酸盐化、中温青磐岩化到晚期的绢英岩化，成矿流体的pH值经历了持续的小幅度下降；流体氧逸度则由初始阶段钾硅酸盐化的强氧化特征，进而在青磐岩化阶段经历了小幅度下降，以及最终绢英岩化阶段的缓慢回升。
LMA反演与GEM正演模拟结果均反映了玉龙铜矿类似的流体pH值和氧逸度演化轨迹(图7)：(1)与钾硅酸盐化蚀变相关的流体具有较高pH值(5.0~7.0)和氧逸度(ΔFMQ=+2.9~+3.6)特征；而与青磐岩化阶段有关的流体pH值(5.0~6.4)和氧逸度(ΔFMQ=+1.1)均有所下降；与绢英岩化阶段有关的流体pH值(5.0~5.7)进一步小幅下降，而氧逸度则(ΔFMQ=+1.7)略有回升。GEM模拟结果显示，受流体中SO2的歧化作用(4SO2+4H2O=H2S+3H2SO4；图4d)的影响，SO2/H2S比值从初始的0.86快速下降为0.01，并伴随着H2SO4的不断形成。这导致流体pH值逐渐降低，并蚀变早期生成的黑云母和长石类矿物，形成了绿泥石和绢云母等矿物(图4f，g)。这一假设得到了玉龙铜矿钾化阶段(-4.5‰)和绢英岩化阶段(+2.1‰)黄铁矿S同位素差异性变化的支持，二者均显示出SO2的歧化作用是控制流体pH值下降的重要因素。此外，模拟显示初始含矿流体具有高氧逸度特征，这与来自岩浆角闪石和黑云母氧逸度计算结果相一致。同时，由于高氧逸度、富含SO2的高温流体可以促进相对难溶的高场强元素和贱金属进入金红石晶体结构中。因此，玉龙矿化脉体中富钨热液金红石的发现，也进一步验证了初始含矿流体的高氧逸度特征。之后，在水岩反应的驱动下，高温钾化阶段典型矿物如磁铁矿或赤铁矿开始结晶沉淀，这促进了流体氧逸度在钾化阶段的持续下降(图4e)。之后，在青磐岩化和绢英岩化阶段，受金属硫化物如黄铜矿、黄铁矿和辉钼矿不断沉淀所控制，流体氧逸度开始缓慢回升至ΔFMQ=+1.0~+1.5范围内。
模拟形成蚀变矿物组合与玉龙铜矿各蚀变分带矿物组合特征较为吻合，即早期高温钾化阶段钾长石、钠长石、硬石膏、黑云母和赤铁矿、磁铁矿等矿物组合的形成，和中温青磐岩化阶段绿帘石、铁绿泥石和斜绿泥石的出现，以及晚期低温绢英岩化阶段绢云母和方解石的形成(图4)。通常认为青磐岩化和绢英岩化蚀变的形成是由岩浆热液与地下水混合形成的。然而，也有研究强调其形成主要与岩浆热液温度下降有关，而无需地下水的混合。本次模拟显示，仅仅依靠温度的下降，单一岩浆流体与斑岩体的水岩反应过程即可产生在斑岩系统中观察到的所有蚀变和矿化特征。早期高温钾化蚀变和同阶段低温青磐岩化蚀变，以及更晚期低温的绢英岩化等蚀变分带作为斑岩矿床广泛存在的一大特征，说明斑岩矿床蚀变分带的存在可能与某些特殊的外部条件关系不大，而更可能是由于一些具有广泛意义(强度性质如温度、压力等)指标变化的结果。因此，斑岩矿床中普遍存在同心环状蚀变分带现象及A、B、D脉体穿切关系序列，应是单一岩浆流体冷却的必然结果。
然而，由于我国斑岩矿床，特别是后碰撞型斑岩矿床，受独特的构造演化历程所导致的更快的抬升速率或复杂的构造叠加影响，其矿化与绿泥石或绢云母等青磐岩化和绢英岩化蚀变常具有紧密的伴生关系。目前，玉龙铜矿石英Ti扩散速率研究显示，玉龙铜矿的形成过程经历了复杂的多阶段流体网脉阶段和较短的形成时间。结合本次模拟结果，温度变化是控制蚀变分带中各矿物出现顺序的重要因素。考虑到在斑岩成矿系统中，相比于高温岩浆热液阶段形成的钾硅酸盐化蚀变，青磐岩化蚀变多被认为形成于岩浆热液相对低温阶段，且青磐岩化蚀变可因岩浆热液系统的坍塌而随着大气降水“倒灌”侵入先期形成的钾化蚀变带，造成局部叠加于钾化阶段上，本文认为碰撞环境下较高的同矿化抬升剥蚀速率导致的岩浆热液体系快速冷却可能是控制玉龙铜矿绢英岩化异常发育且强烈叠加到钾硅酸盐化带之上的重要因素。
热液流体中Cu和Mo的沉淀受到多种因素影响，其中包括降温、压力下降及水岩反应等。这些因素可能导致成矿物质在流体中的溶解度降低。然而，不同于浅成低温热液矿床，斑岩环境下Cu和Mo沉淀过程不会引起流体中fO2和pH值的重大变化，这显示水岩反应导致的流体物化指标变化似乎无法使Cu或Mo在以矿化斑岩为代表的小岩体中高效富集沉淀。此外，有关卤水/蒸汽Cu和Mo分配系数研究以及斑岩型Cu-Mo矿床中普遍存在的富Cu或Mo卤水包裹体显示，Cu和Mo优先分配到卤水相中。因此，压力下降导致的相分离也不太可能是斑岩矿床中Mo沉淀的主要原因。因此，在斑岩环境中，温度下降导致的含矿卤水凝结可能是将Cu和Mo富集并沉淀矿化在小岩体中的关键。
已有研究显示，玉龙铜矿钾硅酸岩化蚀变带内脉体中可见辉钼矿，而与绿泥石或绢云母相关脉体则主要是黄铜矿和黄铁矿，且不同脉体石英阴极发光结果显示辉钼矿沉淀温度高于黄铜矿。与之类似，模拟结果显示，以HMoO4-和MoO42-为载体的Mo元素在狭窄高温区间内(450~370℃)溶解度快速下降，而以CuCl(CuCl43-、CuCl2-、CuCl)为主要络合离子的Cu在中-高温阶段(450~300℃)平稳下降。此外，随着模拟过程中含矿流体反应轮数的不断增加，黄铜矿在450~200℃整个温度区间内沉淀质量均稳定增长，而Mo在450~330℃内的沉淀量进一步加大，但在300~200℃内涨幅较小。这与斑岩矿床中常见的Mo矿化多发育在靠近岩体侵入部位，Cu呈网脉状-浸染状叠加到Mo矿化之上等现象一致(图7)，均显示出温度变化对于斑岩矿床矿化分带的重要控制作用。由此可见，Cu和Mo络合离子稳定性对温度的差异响应是与玉龙铜矿富含辉钼矿脉体和以Cu-Fe硫化物为主脉体解耦分布所引起的局部空间上Cu、Mo分离密切相关。
(1)与钾硅酸盐化蚀变相关的流体具有较高pH值(5.0~7.0)和氧逸度(ΔFMQ=+2.9~+3.6)特征；而与青磐岩化阶段有关的流体pH值(5.0~6.4)和氧逸度(ΔFMQ=+1.1)均有所下降；与绢英岩化阶段有关的流体pH值(5.0~5.7)进一步小幅下降，而氧逸度则(ΔFMQ=+1.7)略有回升。
(2)斑岩矿床蚀变分带是逐渐冷却的单一岩浆热液与斑岩体不断反应的结果，且不同温度窗口对应着钾硅酸盐化(450~360℃)、青磐岩化(360~320℃)和绢英岩化(320~200℃)蚀变矿物的形成。
(3)以HMoO4-和MoO42-为载体的Mo元素在狭窄高温区间内(450~370℃)溶解度快速下降，而以CuCl(CuCl43-、CuCl2-、CuCl)为主要络合离子的Cu在中-高温阶段(450~300℃)平稳下降。
阳光创译介绍
北京阳光创译语言翻译有限公司(Suntrans)是一家聚焦于矿业和能源领域的翻译和咨询服务提供商。阳光创译于2008年2月成立于北京，美国纽约设有分公司，并在乌干达和巴基斯坦设立有办事处。历时18年来，在董事长吕国博士的带领下，阳光创译快速稳健发展，核心定位已经由“打造中国地质矿业翻译领军品牌”，逐渐延伸扩展成“中国国际矿业能源服务大平台”，涵盖矿业能源领域翻译、会展、咨询、猎头、“一带一路”矿业能源商会、矿业能源媒体等国际服务板块。
阳光创译是国际领先的专业领域多语言服务提供商，是中国翻译协会成员、中国语言服务产业技术创新联盟成员和中国矿业联合会全球地质信息共享委员会理事会成员。阳光创译为中国地质、矿业以及石油领域企业国际化和本地化提供整体语言解决方案，提供矿业和能源咨询服务，助推中国矿业企业的国际化进程。
矿业和能源咨询：矿产资源评估、尽职调查、矿业投融资对接、境外矿业营销和销售咨询、大数据分析客户挖掘。
矿产资源评估：是指对矿产资源进行综合评价的过程，包括矿产储量、品位、开采条件等方面的评估，以确定矿产资源的开发价值和潜力。
尽职调查：是指在进行商业交易或投资前对相关企业或项目进行全面的调查和分析，以确认其真实性、可行性和风险，确保投资方能做出明智的决策。
矿业投融资对接：是指为矿业企业提供投融资服务，帮助企业寻找合适的投资方或融资渠道，促成投资交易的达成。
境外矿业营销和销售咨询：是指为矿业企业提供在海外市场进行营销和销售的咨询服务，帮助企业拓展海外业务，提升产品销售水平。
大数据分析客户挖掘：是指利用大数据技术和分析方法，对市场数据进行挖掘和分析，帮助企业发现潜在客户群体，制定营销策略和提升销售效率。
为国内外个业提供地矿与能源产业链的信息服务、投融资项目对接、技术咨询、法律法规咨询等业务。专家顾问来自于加州大学、迈阿密大学、中国地质大学、中国石油大学、中国矿业大学、北京大学等院与相关科研机构，为海内外的客户提供全方位的地质矿业和能源咨询服务。
矿业人才猎头：阳光创译发展历时18年，建立了庞大的地质矿业能源人才库，拥有上千名具有海外背景的地质矿业能源专业人才。可为国内企事业单位推荐推荐地质矿业和能源及其它领域的知名专家、学者和教授等高端专业人才。
矿业会展：阳光创译已成功举办多届论坛，对境内外矿业投资企业的发展和转型等起到了积极促进作用，已受到广泛关注与好评。企业可以利用各种信息渠道宣传自己的产品，推介自己的品牌、形象，了解本行业最新产品动态和行业
如果您对我们的矿业咨询服务感兴趣或有任何疑问，欢迎随时与我们联系。我们将竭诚为您提供帮助和解答。让我们携手合作，共同开创矿业领域的新篇章
阳光创译的核心定位已经由最初的“阳光创译=矿业翻译”——中国地质矿业翻译领军品牌逐渐延伸扩展成“阳光创译=中国国际矿业服务大平台”，涵盖矿业翻译、矿业会展、矿业咨询、矿业猎头、“一带一路”矿业商会、矿业媒体等国际矿业服务板。
热文推荐：
● 阳光创译吕国博士赴哈萨克斯坦矿业考察
● 阳光创译吕国博士考察蒙古国：矿业新政激活中蒙合作新机遇
● 紫金矿业豪掷 12 亿美元，收购哈萨克斯坦金矿
● 蒙古国矿产资源：地大物博的矿业宝藏
● 西芒杜铁矿项目迎重大突破：2025 年 7 月钢轨落位，全球铁矿格局加速重塑
●《中国矿业报》专访阳光创译吕国博士
● 阳光创译 | 矿业翻译服务
● 专业矿业能源翻译服务，帮您开启全球合作之门
相关业务需求即可致电下方相关责任人
| 吕国博士 |
董事长 |
1127203572@qq.com |
| 电话：13401184506 |
微信：suntrans-mining |
| 孙玉梅 |
总经理 |
853495100@qq.com |
| 电话：010-82865216 |
微信：853495100 |
| Gary |
矿业外语猎头 |
582934973@qq.com |
| 电话：13311203631 |
微信：suntrans222 |
| Bruce |
英语翻译 |
1935014417@qq.com |
| 电话：13693672413 |
微信：13693672413 |
| 刘鑫 |
论文翻译 |
1749173384@qq.com |
| 电话：13683660873 |
微信：lvy_09260312 |
| 苗菁 |
地矿报告翻译 |
2534271235@qq.com |
| 电话：13683661813 |
微信：zyy-suntrans |
| Jerry |
口译 |
2268356581@qq.com |
| 电话：13681439112 |
微信：miningtranslation |
| 范瑞瑞 |
媒体合作 |
1985130191@qq.com |
| 电话：17310253590 |
微信：suntrans2014 |
| 张逸远 |
总经理助理 |
347562171@qq.com |
| 电话：15210897852 |
微信：15210897852 |
| 赵宏 |
矿业咨询部主任 |
1583289628@qq.com |
| 电话：13683660097 |
微信：suntrans_1543 |
| 程波 |
矿业咨询部副主任 |
490148863@qq.com |
| 电话：13683660397 |
微信：13683660397 |
| Mary |
矿产贸易总监 |
1529169438@qq.com |
| 电话：17310657953 |
微信：ha15313128283 |
| 吕彦锦 |
翻译 |
3602237436@qq.com |
| 电话：15701297623 |
微信:suntrans_6787
|
阳光创译的核心定位已经由最初的“阳光创译=矿业翻译”——中国地质矿业翻译领军品牌逐渐延伸扩展成“阳光创译=中国国际矿业服务大平台”，涵盖矿业翻译、矿业会展、矿业咨询、矿业猎头、“一带一路”矿业商会、矿业媒体等国际矿业服务板块。
文章版权声明：本文来源 ：张少颖，和文言，高雪，田成华，肖仪武.2024.斑岩铜矿蚀变分带与成矿机制：玉龙矿床水-岩反应热力学平衡模拟例析.岩石学报，40(6)：1837-1852.，版权归原作者所有，本文不代表阳光创译立场，并对文中观点保持中立，仅供各位阅读者交流参考之目的。本号所转载内容没有任何商业宣传目的，仅供交流，如有侵权，请联系主编删除（主编微信：suntrans2008），另外图片版权归原作者所有，如有侵权请联系我们，我们将会立刻删除！给您带来的不便，尽请谅解！
商务合作：请加编辑微信suntrans2008（添加时请您标明“公众号商务合作”意向）；联系电话：15313128281
图片版权声明：部分图片来源于网络，版权归原作者所有。如涉及侵权或者对版权有疑问，请联系邮件1943581358@qq.com，我们会尽快处理！
