斑岩矿床成矿时间尺度的研究进展：以藏东玉龙斑岩铜(钼)矿床为例 *

陈奇 王长明 祝佳萱 杜斌 段泓羽 石康兴 钱金龙 刘俐君

1. 中国地质大学(北京)地球科学与资源学院，北京 100083 2. 中国地质大学地质过程与矿产资源国家重点实验室，北京 100083 3. 有色金属矿产地质调查中心，北京 100012

斑岩矿床是世界上最重要的岩浆-热液矿床之一，常产出于与板片俯冲密切相关的岛弧和大陆边缘弧环境(Sillitoe, 1972, 2010; Richards, 2003; Cookeetal., 2005, 2014; Wilkinson, 2013; 陈华勇和吴超, 2020)，以及与非俯冲有关的陆内环境和陆陆碰撞环境(侯增谦和杨志明, 2009; Hou and Cook, 2009; Houetal., 2015; Yangetal., 2015, 2016; Wangetal., 2017, 2018c; Dengetal., 2018a; Zhengetal., 2019; 侯增谦等, 2020)。典型的斑岩矿床主要是由位于中-上地壳的岩浆房(通常≥100km3)释放大量热液(通常≥10km3；Richards, 2003, 2011)，岩浆热能驱动的岩浆-热液流体在上地壳发生超常元素富集，从而形成具有经济价值的斑岩矿床(Sillitoe, 2010)。成矿过程的持续时间是控制超大型斑岩矿床形成的重要因素(Chelle-Michouetal., 2017; Chiaradia and Caricchi, 2017; Richards, 2018; Korgesetal., 2020)。然而，对于斑岩成矿系统的多期次岩浆-热液活动，仍缺少详细而精确的时间尺度来限定这些过程。

斑岩矿床的矿化事件最终会导致金属矿物和脉石矿物直接从热液中沉淀，或通过流体-岩石反应代替先前存在的矿物相(Chiaradiaetal., 2014)。与成矿作用有关的矿石和脉石矿物在各个矿化阶段普遍存在，这就为确定斑岩系统中岩浆-热液活动的持续时间提供了最直接的测试对象。精细和准确地限定斑岩成矿系统中岩浆-热液活动的持续时间一直是矿床学研究的热点和难点(Chelle-Michouetal., 2015)。借助于传统的原位锆石U-Pb定年(LA-ICP-MS或SHRIMP)和辉钼矿Re-Os等时线定年(ICP-MS)，一些斑岩矿床成矿的持续时间高达几个甚至十几个百万年(m.y.; Sillitoe and Mortensen, 2010)。这种“超长”的时间尺度，既可能是岩浆多期次侵入或矿化周期性叠加的结果(Chiaradiaetal., 2013)；也可能是由于定年方法有限的测试精度(～2%)，导致无法获得精确限定成矿作用的时间尺度(von Quadtetal., 2011; Schalteggeretal., 2015; Lietal., 2017)。此外，对于某些短时间快速形成的斑岩矿床(例如在几十万年尺度内)，传统的测试精度便无法限定其成矿过程的持续时间。越来越多的研究表明，许多超大型斑岩矿床的金属聚集和沉淀是多期次、周期性的岩浆-热液“脉冲”(pulses)的结果，且单次矿化事件很可能小于一个百万年(Chiaradiaetal., 2014; Lietal., 2017)。因此，要解决成矿作用的时间尺度，就需要更为精确的年代学方法。

现今，借助高精度的同位素定年、热力学数值模拟以及石英的钛扩散模型等方法，斑岩矿床中岩浆-热液活动的时间尺度可以被限定在几千年到几十万年之间(Lawley and Selby, 2012; Weisetal., 2012; Merceretal., 2015)。例如，将石英的钛扩散模型应用于斑岩矿床，岩浆和热液活动的持续时间被限定在几万年的尺度内(Merceretal., 2015; Cernuschietal., 2018; Wangetal., 2022; Zhangetal., 2021)。这一模型通常使用恒定的最高初始温度来模拟石英中钛扩散的时间，获得的时间尺度可以代表该温压条件下最短的扩散时间。这些最短的时间尺度为探索斑岩系统中岩浆停留、岩脉注入、热液活动冷却等过程的时间提供了一个起点。

精确限定成矿过程的时间尺度对于更好地理解斑岩矿床的形成和演化具有极其重要意义。本文在回顾斑岩矿床中各种定年方法限定成矿时间尺度的基础上，以三江特提斯成矿带超大型玉龙斑岩铜(钼)矿床为研究对象，构建石英的钛扩散模型以确定岩浆-热液活动的时间尺度，并与矿物的高精度同位素定年对比。本次研究将对限定玉龙矿床多期次岩浆-热液活动以及斑岩矿床成矿过程的时间尺度等内容提供重要的参考。

1 斑岩成矿系统时间尺度的研究进展

通过矿物系统定年建立的年代学框架，进而确定成矿的时代和持续时间，对理解和完善斑岩矿床的成矿模型具有重要意义(Chiaradiaetal., 2013, 2014)。目前，斑岩系统中确定成矿过程时间尺度的方法主要包括：(1)通过放射性同位素定年直接测定矿物系统的时间尺度；(2)对与成矿相关的热扩散过程进行数值模拟(Cathles, 1997; Weisetal., 2012; Korgesetal., 2020)，并结合热年代学方法(Hickeyetal., 2014)；(3)通过在热液成矿系统中测量流体流速和金属浓度来推断所需的时间(Simmons and Brown, 2006)。

1.1 高精度同位素定年

随着矿物原位微区定年技术的广泛应用，直接测定矿床中代表岩浆-热液活动最早和最新的矿物，就可以从理论上确定斑岩系统的成矿时间(Chiaradiaetal., 2014)。常见的矿物年代学应用包括锆石及热液矿物的U-Pb定年、辉钼矿的Re-Os定年以及部分脉石矿物(如黑云母、白云母或绢云母、明矾等)的40Ar/39Ar定年等(Chiaradiaetal., 2009; Chelle-Michouetal., 2015; Hart-Madiganetal., 2020)。将不同的定年技术结合在一起，可以限定斑岩系统中从岩浆侵位到低温热液蚀变和剥露过程的时间框架(图1)。矿物的同位素定年解决斑岩系统中时间尺度的能力，取决于所用方法相对于成矿过程持续时间的可实现精度。如果岩浆-热液活动的持续时间短于定年方法的不确定度，则无法限定成矿的时间尺度。

图1 斑岩系统中各类同位素定年的封闭温度范围(据Chiaradia et al., 2014修改)Fig.1 Range of closure temperatures for various geochronometers (modified after Chiaradia et al., 2014)

利用化学磨损-同位素稀释-热电离质谱法(CA-ID-TIMS: Isotope Dilution-Thermal Ionization Mass Spectrometry)，锆石U-Pb定年的分析精度和准确度提高了一个数量级，从约±0.4%提高到约±0.1%(Schalteggeretal., 2015)。最新的高精度锆石U-Pb地质年代学研究表明，短时间的斑岩侵位和热液活动常发生在大型岩浆房驱动成矿演化的末期(von Quadtetal., 2011; Chelle-Michouetal., 2015; Buretetal., 2016; Tapsteretal., 2016)，并且已经能够在万年(～10k.y.)的时间尺度内解析母岩浆在几十万年里停留和结晶的历史(Buretetal., 2016; Tapsteretal., 2016; Largeetal., 2018; Rottieretal., 2020)。同位素稀释-负热电离质谱法(ID-NTIMS: Isotope Dilution-Negative Thermal Ionization Mass Spectrometric)辉钼矿Re-Os定年凭借其可以直接获得各阶段热液脉的高精度年龄(不确定度<0.2%)的特点，成为解析多期矿化事件时间尺度的有力工具(Selbyetal., 2007; Lawley and Selby, 2012; Chiaradiaetal., 2014; Changetal., 2017; Lietal., 2017)。通过对智利El Salvador和El Teniente两个超大型斑岩矿床使用高精度辉钼矿Re-Os定年，获得其成矿持续时间分别为0.6m.y.和1.7m.y.，其单次热液矿化的时间尺度可以小于100k.y.(Zimmermanetal., 2014; Spenceretal., 2015)。国内学者也利用该方法成功确定了多个超大型斑岩铜(钼)矿床的时间尺度。例如，Lietal. (2017)通过对西藏超大型驱龙斑岩铜矿使用高精度辉钼矿Re-Os定年，发现矿床在266±13k.y.的成矿时间内，可能至少存在两个短期的A-B-D脉旋回(时间尺度分别为38±11k.y.和59±10k.y.)和一个不完整(B-D脉)旋回。由此表明，周期性的岩浆-热液活动既可以出现于持续时间较长的斑岩矿床(例如几个百万年; Chiaradiaetal., 2014)，也可以在短期形成的斑岩矿床中循环发生。与大多数斑岩型铜矿不同，Zhaoetal. (2021)通过对我国东北超大型岔路口斑岩钼矿的年代学研究发现大部分钼金属发生沉淀的时间尺度小于650k.y.，即单次的、短暂的岩浆-热液活动同样具有形成超大型钼矿的潜力。

目前，已经有研究认为矿物的高精度同位素定年更适合于解决中新世和较年轻矿床中的几十万年尺度的热液事件；而对于较老的斑岩矿床，其定年不确定性增加会使年龄的分辨率变差(Lietal., 2017)。因此，斑岩系统中放射性同位素定年方法的选择需要结合实际情况，谨慎选择和考虑。

1.2 热力学模拟与钛扩散模型

通过对斑岩岩体的冷却速率研究发现，受热力学约束与成矿相关的热液系统的最大寿命仅仅为数万年(Cathles, 1997)。热力学数值模拟还表明，斑岩系统的形成可能发生在5～10万年的范围内(Weisetal., 2012)。从活跃地热系统获得的数据(Simmons and Brown, 2006)以及对年轻热液系统的U-Th不平衡定年也支持了这些模型约束(Grimesetal., 1998)。这些数据表明理论上仅需要几万年就能形成超大型斑岩矿床(Simmons and Brown, 2006)。此外，通过模拟斑岩系统中岩浆的幕式注入和多相流体活动，Korgesetal. (2020)认为与斑岩成矿有关的岩浆房是通过岩浆快速、间歇性注入形成的。

图2 三江特提斯地质构造简图(a,据Deng et al., 2017, 2018b; Wang et al., 2015a, b绘制)及始新世玉龙斑岩铜矿带地质简图(b,据Hou et al., 2003绘制)图b中年龄数据来源：Hou et al., 2003; Jiang et al., 2006; Chang et al., 2017; Huang et al., 2019b; 梁华英等, 2008Fig.2 Simplified tectonic framework of the Sanjiang Tethys metallogenic belt and its adjacent areas (a, modified after Deng et al., 2017, 2018b; Wang et al., 2015a, b) and simplified geological map of the Eocene Yulong porphyry Cu belt (b, modified after Hou et al., 2003.)Data sources in Fig.2b: Hou et al., 2003; Jiang et al., 2006; Liang et al., 2008; Chang et al., 2017; Huang et al., 2019b

图3 玉龙斑岩矿床地质图(据Chang et al., 2018绘制)Fig.3 Geological map of the Yulong deposit (modified after Chang et al., 2018)

图4 玉龙矿床矿化、脉体类型、斑岩侵入体及热液蚀变关系示意图Fig.4 Schematic space-time plot for mineralization, veins, and alteration relative to the intrusive history at Yulong deposit

图5 玉龙矿床含矿热液石英脉的手标本、镜下以及石英的阴极发光图(a、b)石英-黄铜矿-黄铁矿(A3)脉切穿早期黑云母(EB)脉和石英-钾长石-黑云母(A1)脉的手标本和透射光照片；(c)A3脉主要由半自形、具有弱振荡分带的石英构成；(d)石英-黄铜矿脉(B3)被晚期黄铁矿-石英脉(D脉)切穿的手标本照片；(e)石英黄铜矿脉(A3)的透射光照片；(f)B3脉石英显示出良好的生长环带且充填黄铜矿. 白色点位为使用EMPA测定的钛元素含量(×10-6).矿物代号：Bt-黑云母；Ccp-黄铜矿；Kfs-钾长石；Py-黄铁矿；Qz-石英；Rt-金红石Fig.5 Hand specimen photographs, photomicrographs and CL textures of the hydrothermal quartz veins at Yulong deposit(a, b) early biotite-quartz (EB) veins and quartz-K-feldspar-biotite (A1) veins cut by quartz-K-feldspar-pyrite-chalcopyrite (A3) veins; (c) A3 veins are dominated by subhedral quartz with oscillatory zoning; (d) quartz-chalcopyrite (B3) veins cut by late pyrite-quartz (D) veins; (e) B3 veins shown under plane polarised light; (f) bright-CL euhedral growth zones and filled with chalcopyrite in B3 veins. Representative EMPA analyses are shown with Ti concentrations given in ×10-6 (white font). Abbreviations: Bt-biotite; Ccp-chalcopyrite; Kfs-K-feldspar; Py-pyrite; Qz-quartz; Rt-rutile

基于脉石矿物中流体-岩石平衡的扩散模型可以为各种地质过程的时间尺度提供物理约束，并且已经成为一种可以量化几年至几万年时间尺度的方法(Matthewsetal., 2012)。斑岩矿床中广泛分布不同阶段的热液脉体，石英凭借其物理化学性质的稳定性，成为记录岩浆-热液活动的重要载体。前人的研究表明，在高温条件下，斑岩矿床中石英生长结构的阴极发光(CL)亮度与钛元素的浓度密切相关(Rusketal., 2008)。利用石英中钛的浓度、钛石英温度计(Huang and Audétat, 2012)和钛扩散速率(Cherniaketal., 2007; Jollandsetal., 2020)相结合的方法，可以建立石英生长结构与钛浓度梯度的扩散模型。该扩散模型已经成功应用于确定火山系统喷发的时间尺度(Chamberlainetal., 2014; Seitzetal., 2016; Ackersonetal., 2018)和控制变质事件的时间尺度(Spearetal., 2012)。最新的研究显示，基于扩散年代学的钛扩散模型已经作为一种有效的工具应用于斑岩系统，约束成矿事件的持续时间以及不同阶段热液活动的时间尺度(Merceretal., 2015; Cernuschietal., 2018; Chenetal., 2021b; Zhangetal., 2021; Wangetal., 2022)。

2 藏东玉龙斑岩铜(钼)矿床应用实例

2.1 矿床地质背景

三江特提斯经历了增生造山到碰撞造山的复杂构造演化过程(邓军等, 2012, 2019; 杜斌等, 2016; 王长明等, 2017; 陈奇等, 2019; Dengetal., 2014, 2021; Wangetal., 2014b, 2018b, 2020; Chenetal., 2021a)。该研究区发育有一系列斑岩-矽卡岩型矿床，前人对其进行了综合分析和研究(邓军等, 2020; 毕献武等, 2019; Wangetal., 2014a, 2018a; Yangetal., 2019)。玉龙斑岩铜矿带位于三江特提斯构造带中部的东羌塘地块，与新生代碰撞造山作用相关(图2a; Wangetal., 2016)。这条具有重要的研究价值和经济意义成矿带南北延长达300km，东西宽15～30km，分布有多个中型到大型的矿床(图2b; 唐仁鲤和罗怀松, 1995; Houetal., 2003)。

玉龙斑岩铜(钼)矿床位于该矿带的北端，蕴藏的Cu金属量超过6.24Mt(@0.62%；杨志明等, 2020; Yang and Cooke, 2019)，同时伴生有大量的Mo、Au等金属，是该矿带储量最高、规模最大的矿床。研究区沉积地层以三叠纪碳酸盐岩和碎屑沉积岩为主，并被始新世多期斑岩体侵入(马鸿文, 1990; 唐仁鲤和罗怀松, 1995; Houetal., 2003; 姜耀辉等, 2006a, b)。岩浆侵入体和矿体均受多条环形断裂和北西向恒星错背斜控制(图3)。玉龙铜矿的矿化类型主要包括：斑岩侵入体内部的细脉浸染状的石英-黄铁矿-黄铜矿矿化，蚀变角岩带中的黄铁矿-黄铜矿-辉钼矿矿化，以及表生富集带内层状或透镜状的氧化矿化(唐仁鲤和罗怀松, 1995; Houetal., 2003)。热液蚀变具有呈同心状叠加发育的特征，在斑岩侵入体的深部和中部发育钾硅酸盐蚀变，自内向外逐渐叠加了绢云母和泥化蚀变，外围发育广泛的青磐岩化蚀变(唐仁鲤和罗怀松, 1995; Houetal., 2003)。

玉龙斑岩铜(钼)矿床发育多种类型的含矿热液石英脉(Changetal., 2017; Sunetal., 2021; 孙茂妤等, 2015)。根据脉体的横切关系、空间分布和矿物组合特征，将这些热液脉划分为三个阶段，即早阶段EB脉(黑云母+石英±钾长石±磁铁矿)和A脉(石英+钾长石±黑云母±辉钼矿±黄铜矿)；主成矿阶段B脉(石英+黄铜矿+辉钼矿+黄铁矿±辉铜矿±黝铜矿)；以及晚阶段D脉(黄铁矿+石英+黄铜矿±斑铜矿±铜蓝)，并伴随有不同类型的热液蚀变(图4)。根据矿物组合及相互关系，A脉和B脉均可以继续细分为三种不同类型(Chenetal., 2021b)。本文重点对早期成矿阶段(A3脉)和主成矿阶段(B3脉)的含矿热液石英脉进行研究。A3脉切穿早期的EB脉和A1脉(石英-钾长石-黑云母脉)，主要由石英、黄铜矿和钾长石组成(图5a, b)；B3脉被晚期D脉(石英-黄铁矿脉)切穿，脉体主要由石英和黄铜矿组成，部分黄铜矿边缘可见辉铜矿(图5d, e)。这些热液石英在阴极发光(CL)图中均可见亮暗交替的生长环带，且较亮的位置显示出较高的钛含量(图5c, f)，表明石英CL强度与钛含量具有一定的正相关关系。此外，热液金红石与硫化物伴生或单独出现(图5f)，表明热液脉中钛浓度可能处于饱和状态。

前人已对玉龙矿床的岩浆活动和铜钼矿化开展了大量的年代学研究(Houetal., 2003, 2006; Jiangetal., 2006; Liangetal., 2006; Lietal., 2012; 梁华英等, 2008; 王成辉等, 2009)，成岩成矿时代得到很好限定(Changetal., 2017; Huangetal., 2019b)，这些工作为进一步解析斑岩矿床多阶段岩浆-热液过程的时间尺度奠定了坚实基础。

2.2 实验方法及数据处理

扫描电子显微镜-阴极发光技术(SEM-CL)揭示了石英中独特的生长环带结构，这些结构受控于矿物生长过程中压力、温度以及元素扩散速率的影响(Landtwing and Pettke, 2005; Rusketal., 2006; Frelingeretal., 2015; Yuguchietal., 2020)。已有研究表明，在高温条件下(>400℃)，热液石英中具有亮暗交替特征的环带往往与钛元素的分布具有较强的相关关系(Rusketal., 2008; Rusk, 2012)。扩散年代学的原理是通过一维空间的扩散模型模拟矿物内相邻环带间元素的剖面变化，以确定与元素成分变化相关的矿物生长所需的时间。石英中相邻的亮暗环带具有不同的钛含量，其交界位置的钛浓度梯度允许对不同石英世代之间的钛扩散进行建模。本次研究使用电子探针分析技术(EMPA)在更精细的尺度上(<10μm)测定石英中不同阶段生长环带的钛含量。为构建石英生长结构和钛含量的联系，本文采用ImageJ V1.8获取CL图像中石英亮度的灰度值曲线，并以EMPA测定的元素含量来校正相同位置的灰度值，由此得到的灰度曲线量化了相对CL亮度和钛浓度的关系(Merceretal., 2015)。

实验研究发现，石英中钛的浓度和扩散速率与压力和温度密切相关。随着校正方法的不断改进，地质学家们已经提出了多种适用于不同环境的石英钛温度计(TitaniQ: Thomasetal., 2010; Huang and Audétat, 2012)。在较低的压力条件下(1～10kbar)，Huang and Audétat (2012) 校准得到更适用于地壳浅部岩浆房或热液环境的石英钛温度计：

(1)

式中T的单位为开尔文(K)，P的单位为千巴(kbar)，该温度计已经在其他斑岩矿床中得到广泛应用(Mercer and Reed, 2013; Maoetal., 2017; Cernuschietal., 2018)。斑岩系统的深度和压力估测一直是矿床学领域的研究难点。针对玉龙斑岩矿床的多期次岩浆-热液流体活动，本文采用了多种方法估测其各阶段的压力。根据斑岩体中角闪石的Al地质温度计(Huangetal., 2019a)，估测石英斑晶的形成压力为1.5kbar。热液脉中石英的形成压力则参考最新的流体包裹体研究结果(Changetal., 2018; Sunetal., 2021)，早阶段脉体(A3脉)、主成矿脉体(B3脉)分别形成于1.1kbar和0.8kbar(Chenetal., 2021b)。

考虑到斑岩系统中石英生长的温度和压力条件，假设石英中平行于c轴(001)的扩散活化能为273±12kJ/mol，在700～1150℃的温度范围内，Cherniaketal. (2007)通过使用合成石英和天然石英以及TiO2粉末作为钛源的扩散实验，得到了平行于(001)扩散的阿伦尼乌斯关系式：

(2)

图6 不同钛扩散速率的对比及其在典型斑岩矿床的外推结果Cherniak et al. (2007)和Jollands et al. (2020)实验的扩散方向均平行于矿物结晶的c轴Fig.6 Comparison of different Ti-in-quartz diffusivities and their extrapolation results in typical porphyry depositsThe Ti-in-quartz diffusion experiments of Cherniak et al. (2007) and Jollands et al. (2020) are parallel to the c-axis

式中，R是通用气体常数，T是扩散的起始温度(K)，DTi为钛的扩散速率。此外，Jollandsetal. (2020)最新的研究使用TiO2和SiO2粉末混合作为钛源，提出了在900～1490℃范围内的钛扩散速率公式。在相同温度和压力条件下，本文发现两种方法计算的扩散速率相差可达二到三个数量级。当处于低压(<3kbar)、中低温(<700℃)等类似斑岩系统的条件时，通过Jollandsetal. (2020)计算的扩散速率显得不合理的“漫长”(图6)。相比之下，Cherniaketal. (2007)的钛扩散速率计算结果更适用于热液环境，并已广泛应用于斑岩矿床的研究(Merceretal., 2015; Cernuschietal., 2018)。假设石英中钛的扩散模型建立在一个一维的、与浓度无关的半无限介质中的扩散，本文使用Carslaw and Jaeger (1946)和Crank (1975)计算扩散时间的方程：

(3)

式中C表示钛含量不同的两代石英中钛的浓度沿梯度分布；cmin和cmax是最小和最大钛含量(×10-6)；x是扩散起始到扩散边缘的距离(μm)；t是扩散时间(s)；D是石英中钛的扩散速率(m2/s，公式2)。其中，钛扩散速率的选择是基于最高的温度估计(Merceretal., 2015)。使用EXCEL的互补误差函数确定方程(3)的最佳拟合解，该函数使用最小二乘法解决非线性数据拟合问题，以确定大多数扩散边界的最佳拟合时间。最后，使用卡方检验(χ2)来确定扩散模型中最佳的拟合时间尺度，避免了由于主观的视觉判断造成的人为误差(Borradaile, 2003)。

2.3 数据结果及分析

利用钛扩散模型，本文选取了矿床早阶段(A2脉)和主成矿阶段(B3脉)中石英的典型结构，确定热液活动的时间尺度。对石英CL图中亮暗交替或变化明显的位置进行EMPA分析(图7a、图8a；表1)，并在同一位置使用ImageJ获取灰度值。对于每次线分析，细灰虚线表示原始的使用钛含量校正的灰度值，粗灰线表示使用ImageJ中“平滑”(smooth)处理后的灰度值(图7b、图8b)。石英中钛的扩散模拟以校准后的粗灰线为拟合对象，忽略了原始灰度曲线的细微干扰。粗黑色曲线及其对应的时间尺度表示每个扩散曲线的最佳拟合结果，边界模拟以浅灰色细线表示，时间尺度以年为单位(图7c、图8c)。

钛扩散模拟显示，早阶段(A3脉)热液脉中石英形成、充填和冷却过程所需的时间范围为32000～210000年(图7)；而主成矿阶段(B3脉)的扩散时间与早阶段石英类似，范围为110000～870000年(图8)。这些时间尺度被解释为从石英边界形成到温度降到有效闭合温度以下(约400℃)发生钛扩散的总时间，也可以代表热液活动中离散的、短暂的加热事件所需时间的总和(Merceretal., 2015)。通过模拟结果的对比可以发现(表2)，使用钛扩散模型估算的时间尺度常受到钛浓度梯度、扩散速率(或扩散初始温度)以及扩散距离的影响。例如，在相同的扩散速率，相近的钛浓度梯度条件下，扩散时间随扩散距离的增加而增加(A3脉：A和B)；当扩散距离相同，钛浓度梯度相近时，扩散时间往往又与扩散速率或扩散起始温度呈正相关(B3脉：A和C；B和D)。在其他条件相似时，钛浓度梯度的增加往往会使扩散时间增长(A3脉：A和D)。

表1 热液石英脉的EMPA元素含量(×10-6)及钛温度计的计算结果

表2 玉龙斑岩铜(钼)矿床不同类型热液石英脉(A3脉和B3脉)的时间尺度

图8 玉龙矿床热液脉(B3脉)中石英的钛扩散模拟图Fig.8 Examples of modeling Ti diffusion in quartz from B3 veins (17-08B10) at Yulong deposit

Changetal. (2017)通过对不同阶段的脉体使用的高精度辉钼矿Re-Os定年(ID-NTIMS)，发现玉龙斑岩铜(钼)矿床的成矿过程总体持续时间长达约5.1m.y.，但铜钼矿化集中发生在约1.3m.y.范围内，其中约80%的铜沉淀和富集在约820000年的时间尺度内。这一结果与使用扩散模型确定的最长时间尺度类似，即在几十万年的尺度范围内。此外，使用钛扩散模型可以将时间尺度最短限定在几万年内(如A3脉中D点：32000年)。由此表明，玉龙斑岩铜(钼)矿床热液活动的时间尺度可以被限定在几万年至几十万年内。

图9 超大型斑岩铜-钼矿床的储量与成矿持续时间和时间尺度的对比图使用石英的钛扩散模型、高精度同位素定年(ID-TIMS辉钼矿Re-Os或锆石U-Pb)以及传统定年方法(LA-ICP-MS锆石U-Pb).数据来源：Sillitoe and Mortensen (2010), von Quadt et al. (2011), Mercer et al. (2015), Spencer et al. (2015), Buret et al. (2016), Chang et al. (2017), Li et al. (2017), Cernuschi et al. (2018), Large et al. (2018), Zhang et al. (2021), Zhao et al. (2021), Wang et al. (2022)以及本文Fig.9 The tonnage of selected giant porphyry Cu-Mo deposits versus the duration and time scales of main mineralization processesAge constrained by Ti diffusion model in quartz, high-precision dating (ID-TIMS molybdenite Re-Os or zircon U-Pb), and traditional dating (LA-ICP-MS zircon U-Pb). Data from Sillitoe and Mortensen (2010), von Quadt et al. (2011), Mercer et al. (2015), Spencer et al. (2015), Buret et al. (2016), Chang et al. (2017), Li et al. (2017), Cernuschi et al. (2018), Large et al. (2018), Zhang et al. (2021), Zhao et al. (2021), Wang et al. (2022), and this study

2.4 斑岩矿床时间尺度的约束

最新的研究显示，石英的钛扩散模型逐渐成为约束斑岩系统成矿持续时间以及不同阶段热液活动时间尺度有效方法。例如，Cernuschietal. (2018)使用钛扩散模型成功估测了秘鲁Haquira East斑岩矿床最大岩浆-热液活动持续时间为170000年，并认为矿石可以在不到35000年的时期内发生快速地沉淀。Merceretal. (2015)则使用该方法限定了美国Butte斑岩矿床的时间尺度，认为斑岩岩浆停留和石英晶体形成的时间尺度为50～6000年，而各期次热液石英脉形成和冷却的时间尺度为10～60000年。因此，利用石英中的钛扩散模型，斑岩矿床中岩浆和热液活动的时间尺度被精确限定在几万年之内。

值得注意的是，与国外具有较短时间尺度(最短可达几百年)的斑岩矿床相比(例如，Butte：Merceretal., 2015; Haquira East: Cernuschietal., 2018)，玉龙斑岩铜矿则表现出较长的时间尺度。通过分析这些矿床的地质特征发现，Butte和Haquira East矿床往往形成于更大的深度(约8～10km；Merceretal., 2015; Cernuschietal., 2018)，属于深成斑岩型矿床。这种深度所估测的压力条件会明显大于中浅成的玉龙矿床(深度约2～5km；Changetal., 2018; Sunetal., 2021)，因而会使用更高的起始扩散温度和扩散速率。此外，根据前人的钛扩散实验(Cherniaketal., 2007)，Merceretal. (2015)所外推的钛扩散速率值也发生明显“过高”的偏移。虽然尚不清楚其偏大的原因，但是更快钛扩散速率一定程度上会获得更小的时间尺度。因此，对于钛扩散模型的使用仍需小心谨慎，特别是在精确获得钛含量的基础上，需要结合矿床地质背景或其他实验方法合理估测温度和压力条件。

尽管钛扩散模型不能直接提供斑岩系统成矿过程的绝对时间，但是该模型仍然能够为某一阶段的岩浆-热液活动提供更精确的时间尺度。石英钛扩散模型的广泛使用，使我们能够从矿物的细微结构入手，基于不同阶段或期次的温压条件，来分析岩浆-热液流体活动的时间尺度。同时，精确限定斑岩矿床成矿过程的时间尺度离不开各种方法的协同合作。矿物高精度同位素定年所获得的绝对时间，与使用元素的扩散年代学获得的相对时间结合，可以在更为精细的尺度上完善斑岩矿床岩浆-热液活动的时间框架。

随着人们对斑岩成矿过程的研究程度不断提高，超大型斑岩矿床时间尺度的精度得到明显提升(Schalteggeretal., 2015; Spenceretal., 2015; Cernuschietal., 2018; Zhaoetal., 2021)。整体而言，传统定年方法确定了百万年尺度的成岩成矿持续时间；高精度同位素定年随测试精度的提升而获得几十万年尺度的持续时间；而石英的钛扩散模型从矿物结构的元素扩散年代学角度可以解析几万年甚至更小级别的时间尺度(图9)。因此，斑岩系统的母岩浆房和成矿过程可能整体经历了百万年尺度的演化，但相对短暂的、几万到几十万年甚至更短时间尺度的岩浆-热液活动对金属超常富集过程具有重要影响(Changetal., 2017; Lietal., 2017; Zhaoetal., 2021)。通过对岩浆成因的矿物学研究也表明，下伏岩浆房的生长受原始岩浆多次幕式注入的影响，每次注入补充的关键成矿物质可以有效地释放到溶出流体中(Lietal., 2018; Zhengetal., 2020)，这可能也决定了斑岩矿床的成矿体量。

3 结论及展望

超大型斑岩矿床的形成往往与多期次的岩浆-热液成矿事件发生叠加有关。岩浆-热液活动的时间尺度对限定和量化斑岩矿床的形成过程具有重要意义。本文以玉龙超大型斑岩铜(钼)矿床为对象，重点识别和剖析热液脉中普遍存在的石英，利用钛元素的扩散年代学方法，来精确限定斑岩矿床中多期岩浆-热液活动的时间尺度。研究结果显示玉龙矿床中热液活动的时间尺度为32000～870000年，有力支持了超大型斑岩矿床可以在几万至几十万年甚至更短时间内形成的观点，同时也为斑岩矿床成矿时间尺度的研究提供了成功范例。这些短暂的成矿事件可能包括岩浆-热液活动的单次性或周期性补充，从而使金属在较长的成矿期内不断发生富集和沉淀，最终导致超大型斑岩矿床的形成。

目前，使用石英钛扩散模型的案例仍然较少且多集中于超大型斑岩矿床，尚无法确定中小型斑岩矿床的岩浆-热液活动是否也具有类似的矿物生长条件和时间尺度。随着人们对矿床学研究的不断深入，与矿物生长密切相关的元素扩散年代学不仅有助于限定斑岩成矿的时间尺度，而且可尝试和拓展到其他类型矿床的研究。

致谢论文的完成得益于邓军教授等老师的指导。野外工作期间西藏玉龙铜业股份有限公司刘申态主任给予了大力支持；实验过程中得到了中国地质科学院矿产资源研究所陈振宇研究员和核工业北京地质研究院邓刘敏硕士的大力帮助。两位匿名审稿人提出了宝贵的意见和建议，使文章得以完善。在此一并致以诚挚的谢意。