川西扎乌龙-青海草陇花岗伟晶岩型稀有金属矿床磷灰石地球化学特征及地质意义*

严清高 李建康 李超 陈振宇 熊欣

1.自然资源部成矿作用与资源评价重点实验室，中国地质科学院矿产资源研究所，北京 100037 2.北京大学地球与空间科学院，造山带与地壳演化教育部重点实验室，北京 100871 3.中国地质科学院国家地质实验测试中心，北京 100037

锂铍等稀有金属作为地球低碳化发展的关键元素(Sovacooletal.,2020)，其相关的成矿理论及找矿勘探成为了当今世界各国的研究热点。花岗伟晶岩型稀有金属矿床作为稀有金属的传统来源(Shawetal.,2016)，提供了全球一半的锂资源(Bensonetal.,2017)，作为最重要的矿床类型之一，受到了国内外学者的广泛关注(李建康等,2007;Thomas and Davidson,2013;London,2018;王汝成等,2019)。我国松潘-甘孜造山带内，产出许多大型-超大型花岗伟晶岩相关的锂铍稀有金属矿床，如甲基卡(王登红等,2005)、可尔因(李建康等,2006)、李家沟(Feietal.,2018,2020)、党坝(费光春等,2020)等。长期以来，学界对松潘-甘孜造山带相关稀有金属矿床开展了大量的调查与研究，大部分学者认为该造山带内花岗岩是地壳重熔的结果(De Sigoyeretal.,2014;Zhaoetal.,2018a,2021)，而空间上与之关系密切的伟晶岩是该类花岗岩分异的产物，岩浆的分异伴随着锂元素的富集(李建康等,2014,2017,2019;王登红等,2017)。近来，部分学者则强调原岩性质和部分熔融过程对稀有金属元素富集成矿具有重要的制约作用(Fanetal.,2020;Zhaoetal.,2022)。值得指出的是锂铍等稀有金属在地球圈层内丰度极低，一般为n×10-6级以下(Rudnick and Gao,2003)，矿床却达百分级含量(如扎乌龙1.2%～1.5% Li2O，Yanetal.,2020)，富集百倍甚至千倍。然而，对于锂矿元素的超常富集，究竟是部分熔融还是岩浆分异过程起主导作用仍不甚明了；岩浆演化过程中的信息记录问题尚未完全解决，尤其是缺乏矿物学证据的制约。以往研究多以全岩地球化学分析手段来解析稀有金属成矿机制(Deschampsetal.,2017;李兴杰等,2018;李贤芳等,2020)，然而微弱的蚀变、矿化以及风化作用等均会给全岩成分带来显著的改变，同时伟晶岩晶体颗粒粗大，很难获得具有代表性的地球化学组分，因此花岗伟晶岩体的全岩成分实际上并不能准确代表相应岩浆熔体的组成，限制了对花岗伟晶岩型稀有金属成矿机制的理解。

最近研究表明，磷灰石晶体结构内包含大量REEs、Sr、Mn、U、Th、F、Cl等，是控制岩浆岩或熔(流)体中稀土元素、卤素含量的重要载体矿物(张荣伟等,2019)。因其在各种地质环境或过程中具有较强的稳定性(Bouzarietal.,2016)，能灵敏地记录和保存原始岩浆信息，常被用于示踪岩浆的物质来源(刘敏等,2021)、了解岩浆演化过程以及氧逸度变化等(Caoetal.,2021)。花岗伟晶岩型稀有金属矿床通常与过铝质的S型花岗岩相关，磷灰石在该类花岗质岩浆中的溶解度随铝饱和指数的增大而增大(Wolf and London,1994)，因此，与稀有金属成矿相关的岩体中含大量P元素，在成岩-成矿过程中形成大量磷灰石。然而，关于花岗伟晶岩成矿体系中磷灰石的系统研究仍较为薄弱。

扎乌龙-草陇锂矿床位于松潘-甘孜造山带中西部，是一个典型的花岗伟晶岩型稀有金属矿床，具有与造山带内甲基卡稀有金属矿床相似的地质特征(Xiongetal.,2019;熊欣等,2021)，为揭示造山带内稀有金属元素富集成矿机制研究提供了理想场所。本文以扎乌龙-草陇矿床中磷灰石为对象，在细致的结构分析基础上，开展了精细的原位主量元素、微量元素及Sr同位素分析测试。试图从矿物学角度揭示扎乌龙-草陇矿床的物质来源，成岩、成矿的岩浆演化过程，为花岗伟晶岩型稀有金属矿床成因机制研究提供借鉴依据。

1 区域地质背景

松潘-甘孜造山带位于青藏高原东北缘，其北部以东昆仑-阿尼玛卿缝合带，东南部以龙门山断裂带，西部以金沙江缝合带为界(图1a)，在三叠纪晚期古特提斯洋闭合过程中，由华北板块、羌塘-昌都板块、扬子板块相互俯冲、陆内碰撞作用下所形成(许志琴等,1992,2018)。造山带内出露地层主要有分布于造山带东部的新元古代-古生代结晶-变质杂岩组成的基底，上覆中-上三叠统为主的巨厚复理石沉积盖层(10～15km，许志琴等，1992；Weislogel,2010)，以及在带内零星分布的少量第四纪沉积物。巨厚三叠纪复理石沉积物是在古特提斯洋闭合过程中，经周缘地块风化剥蚀产物堆积而成(Rogeretal.,2010;Gongetal.,2021)，统一称为西康群变沉积岩。西康群进一步划分为中-上三叠统杂谷脑组(T2-3z)以及上三叠统侏倭组(T3zw)、新都桥组(T3x)和雅江组(T3y)，岩性主要由一系列灰黑色长英质砂岩、粉砂岩、绢云母板岩、含碳质页岩、泥质板岩及少量灰岩等组成(Weislogeletal.,2010)。受印支末期造山运动影响，区内发育大量北西和北西西向褶皱构造；而燕山期水平挤压作用使得区内北东向断裂构造发育(许志琴等,1992)。受印支末期到燕山早期构造运动影响，区内形成大量构造穹隆，包括甲基卡、马尔康、雅江等穹隆(De Sigoyeretal.,2014)。围绕构造穹隆，变沉积岩常具高温-中低压巴罗式变质分带特征，基于变质矿物组合，从穹隆幔部由内向外可划分为夕线石带、十字石带、红柱石带、石榴石带、黑云母带和绢云母-绿泥石带(许志琴等,2018)。

松潘-甘孜造山带广泛发育印支末期-燕山早期中酸性花岗岩，由于其显著的构造意义及稀有金属矿化潜力，受到学者广泛关注。以往研究基于年代学、矿物学及地球化学属性将造山带内花岗岩划分为强过铝质S型花岗岩(De Sigoyeretal.,2014;Deschampsetal.,2017;Yanetal.,2020)、低钾到高钾钙碱性I型花岗岩(Yuanetal.,2010;Deschampsetal.,2017)以及高钾碱性A型花岗岩(Wuetal.,2021)。这些花岗岩于～227Ma到196Ma间侵位至西康群地层内(Zhengetal.,2020)，与锂铍稀有金属矿化相关的主要为强过铝质S型花岗岩(Yanetal.,2020)。围绕花岗岩，常发育不同分带伟晶岩，多呈脉状产出于花岗岩体内外接触带。高度分异的伟晶岩常伴随强烈的锂铍稀有金属矿化，分布于花岗岩体外围。造山带内，已勘探出大量锂铍稀有金属矿化相关的伟晶岩体，包括甲基卡、党坝、李家沟、可尔因等(图1a)，松潘-甘孜造山带因此成为了全球著名的锂铍稀有金属成矿带(李建康等,2014,2017)。其中，扎乌龙-草陇矿床位于造山带中西部，因具有巨大的Li2O资源量，且随着找矿勘探的不断突破，有望成为超大型稀有金属矿床而受到地质工作者的普遍重视(李兴杰等,2018)。

2 矿床地质背景

川西扎乌龙-青海草陇稀有金属矿床横跨四川石渠县和青海玉树州，因矿床主体位于扎乌龙，故本文简称扎乌龙矿床，出露地层主要为三叠纪西康群及覆盖率超过50%的第四系(图1b)。矿区出露岩体为白云母花岗岩，沿扎乌龙背斜褶皱核部呈不规则纺锤状侵入至三叠纪西康群地层内，出露面积约58km2。岩体顶部局部残留西康群围岩盖层(图1c)，表明剥蚀程度较低。岩体与围岩界限清楚，白云母花岗岩形成于212Ma，具强过铝质S型花岗岩特征，为矿区成矿母岩(李兴杰等,2018;Lietal.,2019;Xiongetal.,2019)。岩体具一定的结构分带：中部显示中-细粒花岗结构(中间相，图2a)，边部为中粒花岗结构(边缘相，图2b)，同时边部局部地段演化出似伟晶岩相(图2c)，暗示花岗质岩浆的分异演化过程(Yanetal.,2020)。尽管岩体具有一定结构分带特征，但其矿物组成基本一致，主要由石英、长石和云母等矿物组成，中间相蚀变不明显，而边缘相轻微蚀变。石英呈他形粒状结构，含量40%～55%；斜长石为半自形板状，具轻微黏土化，占20%～30%；碱性长石含量15%～20%，呈半自形板状；白云母含量5%～10%，呈他形或半自形片状，局部见由构造应力作用而形成的塑性变形，具轻微绢云母化，偶尔与黑云母(<5%)共生(图3)。副矿物主要包括锆石、磷灰石、独居石、榍石、石榴石、电气石、钛铁矿、铌钽铁矿、辉钼矿、金红石、锡石等。

矿区伟晶岩成群分布于白云母花岗岩内外接触带，以外接触带最为发育，尤以白云母花岗岩体南侧最为集中(图1b)。目前，矿区已统计伟晶岩脉111条，矿化率为36%(李建康等,2007)。伟晶岩脉规模不等，一般长80～300m，宽2～10m，多呈长条脉状或少量透镜状产出。根据伟晶岩主要矿物组合，从花岗岩体向外依次可划分为：微斜长石伟晶岩、微斜长石-钠长石伟晶岩、钠长石伟晶岩、钠长石-锂辉石伟晶岩和锂云母(白云母)伟晶岩。其中钠长石-锂辉石型伟晶岩和锂云母(白云母)型伟晶岩为矿区主要的稀有金属成矿脉体，已探明二十多条该类型伟晶岩脉，Li2O品位为1.2%～1.5%。14号伟晶岩脉为典型的钠长石-锂辉石型伟晶岩，产出巨晶状锂辉石，全脉矿化，Li2O平均品位约1.2%，储量可达10万吨，占扎乌龙稀有金属矿床Li2O储量的60%(李兴杰等,2018;Yanetal.,2020)。14号伟晶岩脉侵入至岩体外围西康群地层内，延伸长约2km，地表出露厚度约5m，具良好分带性，自边缘向核部分别为：Ⅰ-云母石英带；Ⅱ-微斜长石带；Ⅲ-钠长石-锂辉石带和Ⅳ-石英-锂辉石带(图2d)。云母石英带出现于14号脉边部，厚度约0.5m，主要由石英(含量约85%，粒径1～2cm)、白云母(含量约10%，粒径0.1～0.5cm)和斜长石(含量约5%，粒径0.1～0.5cm)组成(图2e)。微斜长石带厚度约1m，主要由石英(含量约55%，粒径～2cm)、微斜长石(含量约35%，粒径0.5～1.5cm)和白云母(含量约5%，粒径～2cm)组成(图2f)。钠长石-锂辉石带厚度约1.5m，矿物主要有石英(含量30%～45%，粒径～0.5cm)、钠长石(含量40%～45%，粒径1～3cm)、锂辉石(含量2%～25%，粒径多数>5cm)及白云母(含量1%～5%，粒径0.2～0.5cm)(图2g)。石英-锂辉石带厚度约2m，分布于14号伟晶岩脉核部，主要由巨粒锂辉石(含量40%～75%，粒径>10cm)和石英(含量10%～30%，粒径2～5cm)及少量白云母(含量3%～5%，粒径0.5～1cm)和钠长石(含量<1%，粒径～0.5cm)组成(图2h)。副矿物主要有石榴石、电气石、萤石、方解石、锆石、磷灰石、独居石、榍石、锡石、钽锰矿、铌钽铁矿、磁铁矿等。

3 样品采集及分析方法

本次研究分别在白云母花岗岩中间相、边缘相以及14号伟晶岩脉各矿物结构带内采集代表性样品。部分样品制成薄片，进行磷灰石形态学观察。部分样品用于磷灰石单矿物挑选，挑选中，将样品粉碎至50～80目，通过重、磁选方法对样品进行筛选，后在双目镜下挑出晶型完好的磷灰石颗粒，其单矿物样品纯度大于99%。将磷灰石颗粒清洗后制成环氧树脂靶，抛光后用于CL拍照。在CL图像下，选择具有代表性的位置进行原位主、微量及Sr同位素分析测试。

3.1 主量元素分析

磷灰石主量元素分析在中国地质科学院矿产资源研究所电子探针实验室完成。采用JEOL JXA-8230设备，测试条件为加速电压15kV，激发电流20nA，电子束斑5μm。测试过程中使用天然样品和人工合成氧化物作为标准样品，包括磷灰石(Ca，P)、石英(Si)、斜长石(Na)、镁铝榴石(Mg)、Fe2O3(Fe)、MnTiO3(Mn)、金云母(F)和石盐(Cl)。Ca和P的特征峰测定时间为10s，其他Na，Mg，Si，Fe，Mn，F，Cl等的特征峰测定时间为20s。挥发性卤素元素放在优先测试位置，如F元素放在第一道第一个进行测试，减少扩散效应；同时，磷灰石测试前、测试过程中以及测试后，均测试标准物质，以监测设备稳定性。测试过程中，每测试50个分析点后测一次标准物质。标准物质数据显示设备运行状态稳定。元素检出限F为392×10-6，Na为80×10-6，Al为47×10-6，Mg为62×10-6，P为151×10-6，K为61×10-6，Fe为187×10-6，Mn为169×10-6，Ti为377×10-6，Ce为653×10-6，Cl为27×10-6，Si为76×10-6，S为82×10-6。精确度与成分含量相关，Ca、P元素优于0.5%，其他元素相对较高，其中F元素精确度为5%。所有数据经过ZAF校正。

3.2 微量元素分析

磷灰石LA-ICP-MS微量元素分析在中国地质科学院国家地质实验测试中心完成。测试在ASI J-200 343nm飞秒激光(Applied Spectra公司，美国)和X-Series电感耦合等离子体质普仪(Thermo Fisher公司，美国)联机系统上进行。测试过程中采用点方式剥蚀颗粒，束斑直径25μm，频率8Hz，能量密度6J/cm2。以He作为运移样品剥蚀颗粒的载气，样品信号采集时间为20s，空白背景值采集时间为30s。以人工合成的硅酸盐玻璃标准物质NIST SRM610和SRM612作为标样，每完成10个样品点测试一组标样以监测仪器运行状态。数据处理采用ICPMS DataCal软件完成(Liuetal.,2010)。微量元素检测限为0.05×10-6～0.1×10-6之间。标样结果表明微量元素分析结果精准度在10%以内。

3.3 原位Sr同位素分析

磷灰石原位Sr同位素分析在中国地质科学院国家地质实验测试中心完成。测试在ASI J-200 343nm飞秒激光(Applied Spectra公司，美国)和Neptune Plus多接收电感耦合等离子体质普仪(MC-ICP-MS,Thermo Fisher公司，德国)联机系统上进行。采用线扫描方式剥蚀样品，束斑直径为40μm，线扫描速度0.65μm/s，频率10Hz，能量密度10J/cm2。载气条件和信号采集方式如3.2节所述，具体设备参数及测试流程详见Lietal.(2018)。国际标样Durango作为磷灰石标准物质，获得87Sr/86Sr平均值为0.7067±0.003(n=8)，与国际地球化学环境标准物质推荐值在误差范围内一致(0.7064±0.0005，http://georem.mpch-mainz.gwdg.de/)。

4 分析结果

4.1 磷灰石产出特征

白云母花岗岩中磷灰石以矿物包裹体形式产出于寄主岩石的白云母、斜长石等矿物中(图3)，主要呈自形-半自形的长柱状-短柱状，代表磷灰石早于云母、长石类矿物结晶，记录了相对原始的岩浆信息。其中，白云母花岗岩中间相的磷灰石表面干净、成分均一，部分颗粒可见明显的核-边生长环带结构(图4a)，为典型的岩浆磷灰石；边缘相的磷灰石边部具轻微溶蚀现象，记录了岩浆演化中流体交代过程(图4b)。14号伟晶岩脉的磷灰石以副矿物形式与白云母、钠长石等共生，颗粒较大，呈半自形长柱状-短柱状，显示与长石、云母类矿物同时结晶或稍晚，记录了岩浆晚期的演化过程(图3)。整体上磷灰石表面较暗，相对不平整，发育较多孔洞、裂隙，局部含有矿物(如独居石)或流体包裹体。溶蚀现象进一步扩大，呈逐渐交代反应或港湾状溶蚀结构，为流体交代作用下的磷灰石，反映了流体交代阶段的元素组成(图3、图4)。

4.2 主量元素测试结果

磷灰石主量元素测试结果见电子版附表1，磷灰石整体上具有较为一致的CaO(51.54%～56.27%)和P2O5(39.62%～43.32%)含量，高F(2.84%～5.05%)、低Cl(<0.11%，绝大部分<0.02%)特征，为典型氟磷灰石。具体而言，白云母花岗岩边缘相(2.99%～4.76%，平均4.01%)相较于中间相磷灰石(3.03%～4.23%，平均3.63%)具有更高的F含量，显示了从岩体核部向边部F含量升高趋势(图5)。同时，就单颗粒磷灰石而言，其核部向边部同样显示F含量升高(图4)。而从花岗岩到伟晶岩，其F含量急剧下降(图5)，14号伟晶岩脉边部带伟晶岩磷灰石F含量变化于2.84%～4.14%，平均为3.65%；而从边部向核部整体具升高特征，其范围分别为3.38%～5.05%、3.28%～4.58%、3.26%～4.56%，平均值分别为4.06%、3.87%、3.94%。此外，MnO含量<2.99%，而其他Na2O、MgO、Al2O3、K2O、SiO2、FeO、TiO2等含量更低(<0.1%)。

表1 扎乌龙矿床磷灰石原位Sr同位素分析结果Table 1 In-situ Sr isotope analytical results for apatites from Zhawulong deposit

4.3 微量元素测试结果

磷灰石微量元素测试结果见电子版附表2，不同岩石磷灰石具有不同的稀土总量及配分模式(图6)。其中，白云母花岗岩中间相磷灰石∑REE变化在2268×10-6～3370×10-6之间，具强烈Eu负异常(Eu/Eu*=0.06～0.11)，显示中稀土富集而轻、重稀土亏损的典型“海鸥”式配分模式(图6a)；边缘相磷灰石呈右倾配分模式(图6b)，具轻稀土富集、重稀土相对亏损、中等-弱负异常(Eu/Eu*=0.09～0.92)特点，其∑REE为763×10-6～1800×10-6，略低于白云母花岗岩中间相。然而，伟晶岩不同分带磷灰石具有较低的稀土总量，以富集轻稀土、亏损重稀土，具强烈Eu正异常为特征(图6c-f)。其∑REE分别为2.10×10-6～339×10-6、5.61×10-6～140×10-6、37.1×10-6～202×10-6和4.01×10-6～137×10-6；Eu/Eu*分别为0.19～654、0.92～11.6、0.85～30.6和4.20～41.1。

微量元素组成上，其Th、U、Th/U比值、Y等元素在不同样品中具较大差异性(图7)。其中白云母花岗岩中间相磷灰石具有较低的Th、U、Th/U比值而高的Y含量值，其平均值分别为1.72×10-6、47.7×10-6、0.04和1803×10-6；边缘相磷灰石具有略高的Th、U、Th/U比值而略低的Y含量，平均值分别为8.96×10-6、73.2×10-6、0.12和172×10-6；伟晶岩磷灰石显示变化的但总体较高的Th、U、Th/U比值而低的Y含量，其范围分别为0.19×10-6～75.8×10-6、0.05×10-6～382×10-6、0.04～1.20和0.02×10-6～109×10-6。整体上，随着∑REE值的降低，Th、U及其Th/U比值具升高趋势，而Y显示降低趋势；从白云母花岗岩到伟晶岩，其Mn含量升高而∑REE+Y含量降低(图7)。

4.4 原位Sr同位素测试结果

磷灰石原位Sr同位素分析结果见表1。原位测试88Sr信号值与磷灰石Sr含量及其设备状态有关，为了获得较为可信的Sr同位素值，在多个测试点中选择88Sr信号大于0.5V的Sr同位素值作为本次磷灰石有效的Sr测试结果(Lietal.,2018)。中细粒白云母花岗岩磷灰石88Sr信号为0.54～1.43V，平均为0.79V，获得的87Sr/86Sr变化范围较大，为0.7211～0.7333，平均为0.7263；而14号伟晶岩脉磷灰石88Sr信号为0.57～3.63V，平均为1.39V，获得的87Sr/86Sr值较为集中，为0.7134～0.7153，平均值为0.7140。

5 讨论

5.1 成岩成矿物质来源

据地球物理资料证实，松潘-甘孜造山带内三叠纪复理石沉积物厚度可达15km(王海燕等,2016)，造山过程中，由于地壳增厚，在地热梯度、构造剪切以及放射性衰变等产生的热源参与下，将会发生部分熔融作用(Rogeretal.,2010)。如果岩浆熔体来自于深部三叠纪复理石沉积物的部分熔融，磷灰石将记录其一定的地球化学信息。磷灰石晶格内保留显著的卤素成分，可用于示踪岩浆的卤素行为(Websteretal.,2009)。结果显示，无论是白云母花岗岩还是伟晶岩磷灰石，均具有较高的F和异常低的Cl含量(<0.11%，许多低于检出限，附表1)。一般而言，Cl是流体相容性元素，在地壳风化过程中易于丢失，导致残余地壳沉积物具有高F低Cl的特征(Blevin and Chappell,1992)，而经此类地壳物质部分熔融产生的初始岩浆也将继承该卤素组成。因此，磷灰石高F低Cl特征有力地佐证了初始岩浆熔体与三叠纪复理石沉积物部分熔融作用具有紧密的相关性。造山过程中，地壳物质经部分熔融产生的岩浆具有异常低的Cl含量也时见报道，如澳大利亚Lachlan复合造山带内S型花岗岩中磷灰石记录的Cl含量普遍低于0.1%(Sha and Chappell,1999)。

附表1 扎乌龙矿床磷灰石电子探针分析结果(wt%)

续附表1Continued Appendix Table 1Spot No.FNa2OMgOAl2O3CaOP2O5K2OSO3SiO2FeOMnOLa2O3ClCe2O3TiO2Total194.020.070.0150.02354.0842.190.0000.0000.0000.3521.110.0000.0050.0080.000100.2203.920.000.0220.00055.4942.830.0000.0000.0170.0000.070.1650.0030.0000.000100.9ZWL14-1-1,14号脉云母-石英带13.270.020.0000.02753.7341.550.0000.0380.0000.0000.000.1100.0000.0830.00097.4524.110.030.0040.02253.9641.800.0000.0210.0000.1010.160.3290.0000.1300.00098.9434.140.010.0010.02354.1042.070.0000.0080.0170.1250.210.3840.0040.0550.03199.4443.860.000.0000.02454.5642.410.0000.0000.0170.0000.310.0550.0120.0160.00099.6454.080.000.0000.00054.0341.960.0000.0130.0000.0000.270.1100.0000.0900.00098.8464.120.010.0200.01654.6942.050.0170.0250.0150.1250.310.0000.0000.0000.02799.6872.840.010.0140.00454.9942.370.0160.0130.0000.0600.280.0000.0000.0000.00099.4083.140.000.0100.00054.1942.210.0200.0000.0070.0780.250.0000.0080.0000.02798.6293.420.030.0070.00055.2442.530.0080.0000.0000.0000.110.0000.0000.0000.00099.90104.090.000.0040.00054.7242.170.0030.0000.0100.0000.140.0000.0000.0160.00099.43113.780.000.0000.02253.9740.830.0000.0000.0220.0000.360.3290.0100.0000.00497.73123.110.000.0110.00053.9541.580.0000.0000.0070.1730.410.0000.0000.0310.06398.03133.690.000.0170.03854.5241.990.0000.0170.0200.0000.410.0000.0000.0000.00099.14144.060.000.0000.02054.3841.590.0000.0000.0000.0000.210.0000.0050.0000.00098.56154.060.000.0430.00054.1740.990.0000.0210.0350.0000.190.0550.0000.0430.00097.90163.850.000.0000.00054.1341.960.0000.0130.0420.0000.030.2200.0000.0000.00098.62173.860.000.0250.00755.2641.760.0160.0000.0000.0240.140.0000.0000.0000.00099.47183.040.000.0080.00054.8442.180.0160.0170.0000.0000.020.0000.0000.0080.01898.87193.570.000.0000.00054.0142.110.0040.0210.0220.0950.630.5490.0000.0000.01399.52203.460.000.0000.00053.5940.860.0040.0000.0320.0000.830.0000.0080.0000.00097.32ZWL14-2-3,14号脉微斜长石带13.570.010.0010.00052.1740.810.0000.0380.0150.0712.990.2730.0000.0000.00098.4524.030.000.0000.01853.6142.140.0360.0040.0000.0480.720.0550.0030.0270.06799.0633.920.000.0060.01452.6542.100.0040.0300.0000.2141.820.2730.0110.0000.00099.3944.400.000.0000.02253.7341.070.0000.0000.0000.1370.980.0000.0010.0590.00498.5553.800.000.0000.00052.9441.620.0010.0040.0320.1011.880.0000.0110.0000.02298.8164.330.000.0000.00054.8343.320.0000.0000.0220.0780.110.0000.0000.0000.004100.975.050.010.0000.00054.3241.820.0000.0000.0000.0000.110.0550.0000.0240.02799.2884.570.000.0080.02253.3942.210.0000.0130.0170.0361.270.0000.0100.0820.09499.7993.560.010.0000.01151.9242.770.0050.0080.0000.3932.610.1090.0000.0000.03199.92104.440.000.0000.01254.3841.880.0030.0420.0170.0540.810.0000.0010.0000.03199.80114.200.000.0130.01854.0242.040.0000.0250.0270.2270.920.0000.0110.0240.00099.76124.390.040.0000.00055.1041.900.0000.0000.0200.1190.150.1100.0000.0630.000100.0134.290.000.0110.00554.1841.770.0000.0460.0000.1850.820.2740.0070.0000.02299.80144.190.000.0160.00054.8141.670.0100.0000.0100.0000.290.3300.0020.0000.00099.56153.760.000.0040.00055.4041.610.0140.0210.0000.0060.140.0000.0100.0000.00499.39164.030.000.0250.00153.7941.420.0130.0170.0000.1431.010.0000.0170.0000.00998.77

续附表1Continued Appendix Table 1Spot No.FNa2OMgOAl2O3CaOP2O5K2OSO3SiO2FeOMnOLa2O3ClCe2O3TiO2Total174.160.000.0000.00054.1142.270.0070.0000.0000.2030.940.0000.0040.0470.00099.98183.960.000.0000.00054.5841.430.0000.0380.0200.0000.220.0000.0000.0160.00098.60194.170.000.0240.03351.5441.170.0000.0000.0070.2082.600.0000.0000.0000.01398.01203.870.000.0010.00054.0641.570.0000.0000.0000.0481.300.0000.0000.0120.01399.24214.230.000.0000.01755.3842.110.0000.0000.0070.0240.190.1650.0160.0430.022100.4223.700.000.0000.00254.8142.200.0190.0000.0050.0120.360.0000.0000.0080.00499.56233.380.010.0160.00054.8141.850.0130.0210.0100.1020.020.0000.0000.0470.02798.89244.010.000.0000.01054.6641.580.0120.0170.0220.0000.850.0000.0120.0000.00099.48253.880.000.0000.00355.0442.530.0280.0210.0000.0420.180.2750.0020.0710.045100.5ZWL14-3-2, 14号脉钠长石-锂辉石带14.180.000.0000.00052.9941.800.0000.0040.0570.3682.040.0540.0000.0390.02299.8023.750.000.0000.00053.4141.480.0000.0000.0000.2191.900.1630.0090.0350.00099.3833.310.000.0000.01852.9342.160.0110.0210.0300.3082.590.0540.0070.0350.000100.143.500.020.0000.00052.6641.330.0000.0000.0100.4152.130.3270.0000.0040.00098.9254.380.000.0020.00456.1242.010.0050.0000.0000.0590.180.0000.0070.0710.045101.064.310.000.0020.00053.8340.880.0000.0040.0170.2731.440.1090.0020.0590.00099.1173.790.000.0000.01952.8541.770.0000.0040.0300.3082.410.0000.0030.0620.00999.6683.280.030.0000.00053.3842.180.0000.0000.0320.3682.310.0000.0000.0000.000100.293.610.000.0000.00052.5642.100.0000.0130.0250.3852.470.2180.0000.0000.00099.86103.950.000.0000.01652.4340.250.0240.0330.0050.4802.350.0000.0000.0000.00497.88114.030.000.0000.00053.2640.960.0000.0210.0470.3622.430.0000.0090.0000.00099.41124.510.000.0000.00056.2740.220.0120.0000.0620.0000.060.4380.0000.0200.03199.73133.950.000.0000.00655.5341.010.0040.0000.0670.0000.090.0000.0030.0780.00499.07143.340.000.0000.00055.3641.800.0000.0000.0540.0000.000.0000.0040.0000.00099.16154.580.000.0240.01152.1140.760.0080.0340.0150.4562.120.0000.0000.0120.01398.21ZWL14-4, 14号脉石英-锂辉石带13.430.000.0000.00052.5740.840.0130.0220.0120.3222.050.0000.0000.0000.02197.8224.180.000.0160.00152.1841.210.0000.0160.0260.3822.190.0000.0980.0000.02198.5434.380.000.0140.00054.9239.810.0000.0000.0050.0000.100.0000.0210.0620.00097.4643.550.000.0210.00053.3940.550.0000.0000.0610.2231.710.0000.0030.0000.01198.0154.030.000.0000.01554.7139.890.0100.0140.0200.0210.110.0000.0540.0000.00097.1764.310.000.0030.00053.7540.260.0020.0200.0210.0861.280.0000.0330.0520.07998.0774.090.000.0030.02052.8240.440.0000.0000.0090.2832.180.0000.1100.0300.00098.2484.080.000.0180.00053.6740.540.0010.0000.0260.1411.500.0000.0000.0000.01398.2793.260.000.0000.01454.8439.620.0070.0000.0280.0090.390.0000.0130.0000.04496.85104.370.010.0000.00055.1040.880.0060.0000.0340.0000.090.0000.0250.0520.01698.73113.330.000.0080.00255.4939.940.0130.0000.0150.0000.220.0000.0000.0500.01397.68123.820.000.0140.00252.4040.540.0000.0120.0040.2922.390.0000.0480.0580.00097.96134.110.000.0000.01154.3740.410.0000.0100.0460.0000.020.0000.0200.0060.00397.27144.560.000.0000.00052.9541.060.0000.0200.0000.2091.510.0000.0910.0000.05798.53

此外，白云母花岗岩内见钛铁矿、金红石等还原性特征矿物；同时磷灰石显示强烈的Eu负异常和Ce正异常(图9a)，指示寄主岩石形成于较高还原性质的岩浆。事实上，岩浆还原性质还被松潘-甘孜造山带甲基卡、扎乌龙等典型花岗伟晶岩型稀有金属矿床中花岗岩石英内赋存的熔体包裹体里发现的甲烷、氢气和石墨等还原性物质所证实(Lietal.,2014;Li and Chou,2015;Xiongetal.,2019)。三叠纪复理石沉积物包含大量泥质岩(Gongetal.,2021)，甚至野外可见大量碳质泥岩和碳质页岩地层，显示强烈的还原性质特征。一般来说，这种还原性质是由页岩组成的泥质岩自身所具有的典型特性(Richards andengör,2017)。因此，岩浆的还原性质主要继承于三叠纪复理石泥质沉积物的源区。

综合上述证据，本文认为扎乌龙矿床的初始岩浆主要源自于造山过程中深部三叠纪复理石沉积物的部分熔融作用，这与基于Sr-Nd-Hf同位素证据得到造山带内S型花岗岩的主要来源为地壳沉积物的认识相吻合(De Sigoyeretal.,2014;Zhaoetal.,2018b;Yanetal.,2020)。值得一提的是，西康群全岩地球化学数据显示复理石沉积物具有较高的稀有金属含量，尤其是Li(63.6×10-6～147×10-6，Yanetal.,2020)高于地壳丰度数倍。因此，在高F低Cl，较高还原性质条件下，部分熔融导致Li等稀有金属元素从深部三叠纪复理石转移至岩浆中，复理石沉积物不但提供成岩物质，还提供稀有金属成矿物质。然而，考虑到成矿伟晶岩与白云母花岗岩具有紧密的时空与成因关系(李兴杰等,2018;Lietal.,2019;Yanetal.,2020)，其初始岩浆均来自相同的物质源区，因此，物质来源似乎不是扎乌龙稀有金属矿床形成的关键控制因素。

5.2 磷灰石对稀有金属矿床富集成矿的指示

磷灰石微量元素地球化学对岩浆的演化机制具有重要指示意义(张荣伟等,2019)。前文述及，白云母花岗岩与伟晶岩具有紧密的时空及成因关系，起源于相同的物质源区，因此磷灰石不同的稀土配分模式及稀土含量(图6)，以及Th、U、Y、Mn等元素的变化(图7)，指示了白云母花岗岩和伟晶岩的演化过程。磷灰石是稀土完全配分型矿物，具有与寄主岩石相一致的稀土元素配分模式(曾令森等,2012)。以往研究认为，榍石常富集LREEs(Panetal.,1993)，因此白云母花岗岩中间相的磷灰石可能受到榍石在早期或同时结晶的影响，导致其轻、重稀土分馏程度较低；同时强烈的Eu负异常，可能是由于具Eu正异常特征的斜长石早于或同时与磷灰石结晶所致(Xieetal.,2018)。而其边缘相显示轻稀土富集而重稀土亏损的配分模式，则指示富含MREEs和HREEs矿物的结晶，比如锆石和石榴石等，这与野外大量石榴石产于白云母花岗岩的现象相符。伟晶岩阶段，稀土总量进一步降低(图6)，暗示富含稀土的矿物，如独居石、锆石、石榴石、磷灰石、榍石等，在花岗岩熔体阶段进行大量结晶，消耗熔体中大量REEs，进而影响到伟晶岩阶段磷灰石的低REEs含量组成。岩浆的演化过程同样得到微量元素的支持，如岩浆的结晶分异导致残余熔体Th、U值的升高(Miller and Wooden,2004)，而石榴石的结晶导致HREEs和Y值的下降。另外，受熔体结构及成分的影响，高分异熔体结晶的磷灰石具有较高的Mn含量(O’Sullivanetal.,2020;Bromiley,2021)。这些元素的系统性变化(图7)，充分指示了花岗质岩浆的分异演化过程。

岩浆演化过程伴随氧逸度的变化。从白云母花岗岩到伟晶岩，磷灰石Eu异常由强烈负异常到中等-弱负异常再到伟晶岩阶段的强烈正异常(图6)。Eu异常反应了熔体的氧化还原性质或长石的结晶(Bea,1996)，然而在伟晶岩核部带，如石英-锂辉石带，长石含量很低(<1%)，却依然呈现极高的Eu正异常，说明磷灰石结晶体系中Eu呈Eu2+与Ca2+发生了显著的替换作用，因此Eu的正异常可能受氧化还原的影响。一般来说，Eu和Ce是两种差异较大的元素，在氧化状态发生变化的条件下，在磷灰石中显示相反的分配特征(韩丽等,2016)。随着Eu异常值的升高，Ce异常值降低，指示从花岗岩熔体向伟晶岩熔体演化过程中氧逸度逐渐升高(图9a)。全岩的Fe3+/Fe2+值也是氧逸度变化的重要指标(Piccoli and Candela,2002)，松潘-甘孜造山带内花岗岩具有相对较低的Fe3+/Fe2+值，基本处于中等还原区域，而伟晶岩则普遍偏高，处于中等还原至中等氧化区域(图9b)，扎乌龙白云母花岗岩中间相至边缘相同样显示升高的Fe3+/Fe2+值，反应了熔体从花岗岩演化至伟晶岩，氧逸度的增高。另外，花岗岩内普遍发育钛铁矿、金红石等还原性指示矿物，其流体包裹体常富含CH4、石墨或者H2，基本不含CO2；而伟晶岩则出现磁铁矿等氧化性特征矿物，其熔体包裹体常富含CO2而不含或含极少量的CH4(Lietal.,2014;Li and Chou,2015;Feietal.,2018;Xiongetal.,2019;熊欣等,2021)，再次有力地表明了氧逸度随岩浆演化而升高。氧逸度的变化可能是由于伟晶岩熔体自花岗岩分离迁移发生在构造相对封闭到相对开放的转换过程中，伴随H2等还原性组分的逃逸，造成了氧逸度的升高(Rivalta and Dahm,2006)。

值得注意的是，白云母花岗岩中间相磷灰石CL图像显示较为均一、干净的形貌特征，而其边缘相磷灰石具有微弱溶蚀交代结构，而至伟晶岩阶段，溶蚀交代现象进一步扩大，形成明显的港湾状甚至交代残余结构(图4)，暗示随着岩浆熔体的分异演化，晚期存在高度富H2O等挥发分的熔体或者流体，对早期结晶矿物发生溶蚀或交代作用。实验研究表明，REEs和Y元素在磷灰石和熔体中的分配系数随着熔体聚合度的增加而增加(Prowatke and Klemme,2006)。然而，岩浆中挥发分如H2O、F、P、B等的增加将降低岩浆聚合度(Liuetal.,2019)。因此，白云母花岗岩到伟晶岩阶段REEs+Y含量的降低指示了晚期熔体H2O等挥发分的增加(图7)。还需指出的是，随着挥发分的增加，熔体结构发生了变化，从而导致了F的差异性分配行为。相较于流体，F更易于分配进入熔体相。因此在岩浆磷灰石中，无论是单颗粒磷灰石核部至边部，还是从白云母花岗岩中间相至边缘相，F均呈上升趋势，反映了熔体的分异演化过程；而受流体交代的磷灰石，F含量相对降低，但仍显示较高的值，尤其是14号伟晶岩脉内部的分异演化(图5)。这可能暗示了岩浆演化晚期，F在残余熔体中强烈富集并导致相当量的F进入流体相，从而使流体交代的磷灰石F含量整体相对降低但在伟晶岩内部结构带仍有升高趋势。富F流体以氢氟酸形式对岩浆磷灰石进行溶蚀，从而形成港湾状甚至交代残余结构，而部分磷灰石则直接从流体中结晶，显示较暗的形貌特征(图4)。

如前文所述，花岗伟晶岩的初始熔体以富F和还原性质为特征，起源于深部三叠纪复理石的部分熔融作用。事实上，地壳性质的源岩因富含大量含水矿物(如云母等)将会降低岩浆的固相线，容易在造山过程中诱发云母脱水熔融而产生岩浆(Douce and Harris,1998)。值得一提的是，通过类质同象进入云母类含羟基矿物中的F在脱水熔融过程中广泛释放到熔体(Dycketal.,2020)，这种高含量F通过替代破坏熔体结构中的桥氧键深刻影响着硅酸盐熔体挥发分的扩散系数，并增大水在硅酸盐熔体中的溶解度(Dingwell,1989)。此外，原岩中有机质(含碳)的熔融将使岩浆表现为还原性质，同样有利于提高熔体中水的溶解度(Bezmenetal.,1991)。因此，熔体的富F和还原性质促使了初始岩浆含大量溶解水。溶解在岩浆中的水不但降低液相线，而且因解聚作用使熔体聚合度降低，黏度降低，从而使岩浆的结晶分异过程更容易和持续时间更久，这将提高不相容组分(如挥发分和稀有金属元素)在残余熔体的富集程度。据Thomas and Davidson (2012)研究认为花岗岩熔体中的水低于10%，而伟晶岩熔体却能达到10%～40%。因此，岩浆演化晚期的富水熔体易于迁移至花岗岩上部乃至侵入围岩而形成伟晶岩。同时，水含量的增加还减小熔体中微量元素活度，导致元素在矿物/熔体间的分配系数减小(Wood and Blundy,2002)，增强成矿元素的不相容性，因此Li、Be等稀有金属元素更倾向于集中分配至岩浆演化晚期的富挥发分熔体，乃至形成富稀有金属的伟晶岩熔体。还需指出的是，氧逸度的升高导致碳质组分可能主要以CO2/碳酸盐的形式存在(Li and Chou,2016,2017)，促使伟晶岩熔体产生偏碱性环境。这种偏碱性环境能够提高Li等元素的溶解度和促进锂辉石的结晶，形成锂辉石伟晶岩(刘永超,2021;Liuetal.,2022)。总而言之，伴随岩浆的高度分异演化，衍生出一系列有利的物理化学条件，引起熔体结构发生实质性变化，并促进了Li、Be等稀有金属元素在岩浆演化晚期熔体内饱和而富集成矿。因此，岩浆的高度分异演化是花岗伟晶岩中稀有金属元素成矿的关键。

6 结论

磷灰石是揭示花岗伟晶岩型稀有金属矿床物质来源及岩浆演化信息的有效手段。本文通过对扎乌龙-草陇矿床中磷灰石特征的精细分析，得出以下认识：

(1)磷灰石原位Sr同位素、高F低Cl含量，以及还原性质特征，指示扎乌龙稀有金属矿床的初始熔体源自于造山过程中深部地壳三叠纪复理石沉积物的部分熔融作用。

(2)从白云母花岗岩到伟晶岩，连续降低的REEs、Y以及逐渐升高的Th、U、Mn含量，指示了花岗质熔体的分异演化过程。伟晶岩显示负Ce异常和正Eu异常以及典型的溶蚀结构，表明岩浆演化晚期熔体氧逸度的升高以及高度富集H2O、F等挥发分。

(3)熔体富F及还原性特征，能提高初始熔体H2O溶解度，演化晚期较高的H2O含量促使熔体自花岗岩向外迁移而结晶形成伟晶岩。同时，在高氧逸度条件下，CO2营造的偏碱性条件适宜Li的富集和锂辉石结晶。

致谢磷灰石测试及数据处理得到中国地质科学院矿产资源研究所电子探针实验室王臻博士及中国地质科学院国家地质实验测试中心孙鹏程硕士的指导和帮助；同时，承蒙审稿人的认真审阅，提出了宝贵的修改意见。在此一并深表谢意！