花岗伟晶岩型稀有金属矿床流体成矿机制研究进展

郑范博王国光倪 培

内生金属矿床成矿机制研究国家重点实验室, 南京大学地球科学与工程学院, 江苏 南京 210023

0 引言

战略性关键金属对于新能源领域、高性能材料、生物技术、信息技术、高端装备、绿色环保以及航空航天等战略新兴产业具有不可替代的用途,属于高科技发展的核心资源(Linnen et al.,2012;陈骏,2019;蒋少涌等,2019;毛景文等,2019;翟明国等,2019;侯增谦等,2020)。战略性关键金属包括稀有金属、稀土金属、稀散金属和部分稀贵金属。稀有金属包括锂、铍、铌、钽、铷、铯、锆、铪、钨、锡等;稀土金属为镧系元素以及性质相近的钪、钇;稀散元素为镓、锗、硒、镉、铟、蹄、铼和铊;部分稀贵金属包括铂族元素、铬、钴等。文中主要关注与花岗伟晶岩关系密切的稀有金属矿产资源。

在地球科学领域中,流体包括气体、液体、熔体和地球中受构造应力作用而移动的物体(卢焕章,2019)。流体具有流动性、压缩性和压力传递作用,这些力学属性为成矿作用研究提供了基础(徐兴旺等,2019)。成矿流体是金属富集、迁移和沉淀过程中的重要介质。流体包裹体作为保存在晶体中的记录原始成矿流体信息最直接的证据,是获取成矿过程中温度、压力和成分等信息的关键媒介(Roedder, 1984;Bodnar, 2003; 卢焕章等,2004;倪培等,2018;王国光等,2020)。与稀有金属矿产有关的花岗伟晶岩型矿床成矿流体来源、成分、温压条件和演化过程仍存在争议(Norton, 1973;Stewart, 1978; 王联魁等,1999,2000;李健康,2006;Simmons and Webber, 2008;London,2014,2015,2018;Shaw et al., 2016;Thomas and Davidson, 2016; Lv et al., 2018;Thomas et al., 2019;张辉等,2019)。文章以成矿流体为核心内容,回顾花岗岩伟晶岩型稀有金属矿床流体性质和成因机制的研究进展,指出目前研究存在的问题,并且进一步提出相关研究展望。

1 花岗伟晶岩型稀有金属矿床基本特征

1.1 定义及分类

花岗伟晶岩是指具有骨架、文象或其他矿物晶体定向生长的粗粒花岗岩。伟晶岩内的矿物颗粒往往比较粗大,并且常呈现从边缘向核部颗粒逐渐粗大的趋势。伟晶岩中标志性的结构为普遍发育特征性文象结构(图1),并且矿物晶体具有定向生长的习性,而且伟晶岩内常常发育显著的矿物组合空间分带(London,2018)。尽管稀有金属矿床常与伟晶岩有关,值得指出的是,伟晶岩中仅有极少部分发育稀有金属矿化,矿化伟晶岩仅占伟晶岩总量的1%～2% (Cˇerny',1991a;London and Kontak,2012)。伟晶岩有多种分类方案,其中London(2008)将伟晶岩划分为简单伟晶岩和稀有金属伟晶岩,前者主要与长石和石英等非金属矿产有关,后者为稀有金属的重要母岩。应用最为广泛的为Cˇerny'(1991b) 和Cˇerny' and Ercit(2005)的方案,他们将伟晶岩分为LCT(锂铯钽)族和NYF(铌钇氟)族。LCT族伟晶岩富集Li、Be、Ta>Nb、Cs、B、F、P、Mn、Ga、Rb、Sn和 Hf(Cˇerny',1991b; Bradley et al.,2017),主要与S型花岗岩有关,少量与Ⅰ型花岗岩 有 关 (Cˇerny' and Ercit, 2005; Simmons and Webber,2008);NYF族伟晶岩富集Nb > Ta、Y、F、Be、 Sn、 B、 Ti、 REE、 Zr、 Th、 U 和 Sc(Cˇerny',1991b;Rakovan,2008),与A型花岗岩有关(Cˇerny' and Ercit,2005)。LCT族伟晶岩通常是板块汇聚造山过程中岩浆活动的产物(Cˇerny',1991b;Tkachev, 2011; McCauley and Bradley,2014;Tkachev et al.,2018)。此外,LCT族伟晶岩年龄分布的主峰与超大陆聚合的时代一致(McCauley and Bradley,2014;王汝成等,2021)。如2368 Ma、1800 Ma、962 Ma、529～485 Ma和371～274 Ma的LCT族伟晶岩出现的主峰分别与Sclavia and Superia、 Nuna、 Rodinia、 Gondwana、Pangea超大陆聚合时代一致(McCauley and Bradley,2014)。LCT族伟晶岩提供了全世界约二分之一的Li(Benson et al.,2017),十分之一的Be,大量的Ta和几乎全部的Cs(Bradley et al.,2017)。

1.2 产状

花岗伟晶岩脉通常产于母岩花岗岩内部或顶部片麻岩、角闪岩、片岩中,大多数分布在与母岩花岗岩相距数千米范围内(London,2015,2018),有时可达10 km(Baker, 1998)。伟晶岩群的演化和矿化程度往往具有空间分带特征,靠近母岩的伟晶岩演化程度较低,矿化较弱,仅发育造岩矿物;远端的伟晶岩演化程度更高,矿化更强,更倾向于含有锂辉石、绿柱石等稀有金属独立矿物,形成具有重要经济意义的稀有金属矿体(London,2018)。另外,随着距母岩花岗岩距离的增加,伟晶岩脉数目一般会逐渐减少(London,2018)。不仅伟晶岩脉群在整个区域上具有带状分布特征,而且单个伟晶岩脉内也常见带状分布特征(Cameron,1949)。伟晶岩脉通常受围岩构造、组构和层理控制,大小不一,呈脉状、囊状、透镜状、管状、不规则状产出,一般与围岩接触关系明显(Simmons and Webber,2008)。

Qtz—石英;Kfs—钾长石图1 典型稀有金属花岗伟晶岩的文象结构(样品为新疆大红柳滩锂矿手标本)Fig.1 Representative photo of rare metal pegmatites showing graphic texture (The sample is a hand specimen from the Dahongliutan lithium deposit, Xinjiang)Qtz-quartz; Kfs-K-feldspar

1.3 单个花岗伟晶岩脉内部分带

通常单个花岗伟晶岩脉具有显著的内部分带特征(图2)。Cameron(1949)总结了大量资料,提出了广为接受的伟晶岩空间分带模型。从外部到中心,依次发育边缘带、外部带、中间带和内核。边缘带的厚度不大,矿物颗粒细小,主要为长石和石英。外侧带矿物颗粒较粗,主要由长石和石英组成,具有文象结构或花岗结构。中间带主要矿物颗粒粗大,一般由微斜长石组成,有时由长石和石英组成,厚度较大。内核带由石英块体组成。此外,伟晶岩往往结构复杂,常发育后期的裂隙充填体和交代体。Norton(1983)修改了上述模型,提出最后形成的内核带为富钠长石和锂云母或锂辉石。London and Morgan(2017)实验模拟验证了在纯石英结晶之后确实可以形成钠长石-锂云母组合。LCT族伟晶岩矿床中加拿大Tanco矿床发育良好的空间分带(Cˇerny' and Ercit,2005)。Li矿体主要发育在上过渡带,Cs矿体产出在铯镏石带,Ta矿体产出在细晶钠长石带和中过渡带。相似地,超大型Greenbushes伟晶岩矿床中Ta矿化和Li矿化发生在不同位置(Linnen et al.,2012)。中国可可托海3号脉由外向内根据矿物组合特征划分为9个结构带(邹天人和李庆昌,2006),依次为文象、变文象伟晶岩带,糖晶状钠长石带,块体微斜长石带,白云母-石英带,叶钠长石-锂辉石带,石英-锂辉石带,白云母-薄片钠长石带,锂云母-薄片钠长石带,块体微斜长石核和块体石英核,构成同心环状结构分带。但不是所有的LCT族伟晶岩矿床都有良好的内部分带。Swanson(2012)报道了美国北卡罗来纳州和南卡罗来纳州的锂辉石伟晶岩矿床,发现该区域的伟晶岩含有中粗粒锂辉石、斜长石、钾长石、石英、白云母,以及少量的副矿物如绿柱石、锡石等,但是没有显著的空间分带性。Barros et al.(2020)发现爱尔兰东南部的Leinster钠长石-锂辉石伟晶岩为弱分带至均质化。另外,中国四川甲基卡也归类为非显著空间分带性的伟晶岩型矿床(London,2018)。

UST(unidirectional solidification texture)—单向固结组构,位于中间带图2 美国加利福尼亚州圣迭戈县帕洛马山附近一条完整的伟晶岩岩脉部分(29 cm厚;London,2018)Fig.2 A complete section of a pegmatite dike from near Palomar Mountain, San Diego County, California, USA (29 cm-thickness; London, 2018)UST-unidirectional solidification texture, located in the intermediate zone

2 花岗伟晶岩型稀有金属矿床的成矿流体特征

花岗伟晶岩型矿床的成矿流体成分和物理性质是正确理解成矿机制的关键要素。成矿流体中的挥发分(B、P、F和H2O)是稀有金属在溶液中迁移的络合剂,控制着元素的迁移和沉淀;成矿元素含量的变化对于精确理解成矿机制至关重要;成矿流体的物理性质如黏度是理解伟晶岩型矿床成矿机制非常重要的参数。

2.1 成矿流体挥发分

稀有金属倾向于与F和Cl等形成络合物在岩浆热液体系迁移 (Salvi et al., 2000; Veksler,2004; Thomas et al., 2009, 2011a; Zaraisky et al.,2010; Veksler et al., 2012; Timofeev et al., 2015;

Harlaux et al., 2017)。然而,花岗伟晶岩型稀有金属矿床成矿流体挥发分(B、P、F和H2O)含量仍存在争议。部分学者通过全岩分析、模拟计算和实验岩石学研究指出,挥发分含量相对较低,只含有3%～6%H2O(London, 1992;Nabelek et al., 2010)。Morgan and London(1987)通过对加拿大Tanco稀有金属花岗伟晶岩研究,估算了成矿流体中含有655×10-6的B,1419×10-6的F。Stilling et al.(2006)基于对Tanco稀有金属花岗伟晶岩的整体轮廓进行了计算机3D重建,计算了其整体成分,提出Tanco稀有金属花岗伟晶岩含有0.07%的B2O3、0.12%的F、0.86%的P2O5和0.28%的H2O。London (2014) 使用成分类似伟晶岩的富锂流纹质玻璃进行结晶实验,发现产物存在类似于花岗伟晶岩的结构,因此,认为伟晶岩中挥发分(B、P和F)含量<2%,H2O含量为4.5%～5.1%。

然而,另外一些学者通过包裹体研究发现花岗伟晶岩型稀有金属矿床成矿流体可能含有极高的挥发分。如Badanina et al.(2004)对俄罗斯Khangilay杂岩体中Li-F稀有金属花岗岩开展了熔体包裹体成分分析工作,报道最高程度演化的熔体中挥发分质量百分数可以达到11% (H2O:8.6%;F: 1.6%;B2O3: 1.5%)。Thomas and Davidson(2012)总结了一些伟晶岩矿床中的熔体包裹体数据,表明成矿流体中H2O的含量主要在8.1%～26.2%之间。Borisova et al.(2012)使用LA-ICPQMS对Ehrenfriedersdorf伟晶岩的石英中的富晶体包裹体进行的初步分析研究发现B含量可达16500×10-6。熊欣等(2019)运用扫描电镜-能谱仪对甲基卡花岗伟晶岩型锂矿锂辉石中的包裹体(图3)固相矿物进行鉴定发现存在磷灰石和萤石,表明成矿流体中富P和F。Li et al.(2019)通过稀有金属矿物在热液金刚石压腔(HDAC)水溶液中的结晶实验研究表明,成矿流体可能富含H2O。Nabelek et al. (2010)根据实验研究表明,溶解的H2O可以提高成核所需的自由能、降低黏度,在伟晶岩形成中起着关键作用。Maneta and Baker(2014)在水饱和、含锂条件下实验性地产生了石英-长石共生的文象结构。

除了挥发分B、P、F和H2O外,成矿流体中也常见碳酸盐矿物、CO2等挥发分。如London(1986)在Tanco伟晶岩中发现扎布耶石(Li2CO3),Thomas and Davidson(2016)在德国Ehrenfriedersdorf伟晶岩石英中的熔体包裹体中发现气相成分为CO2。Li and Chou(2017)运用激光拉曼对甲基卡锂辉石中的富晶体包裹体研究发现,包裹体中的固相矿物除了硅酸盐矿物(方石英、锂辉石、锂绿泥石)外还含有扎布耶石(Li2CO3)和方解石,室温下,流体相中含10%～80%体积百分比的CO2,这些结果表明,富晶体包裹体代表了一种捕获富含碳酸盐和硅酸盐的流体。因此,成矿流体可能是富碳酸盐和硅酸盐的流体。

花岗伟晶岩型稀有金属矿床成矿流体中的挥发分(B、P、F和H2O)可显著降低花岗伟晶岩结晶温度、成核速率和黏度,并且明显提高了成矿流体中稀有金属元素和大离子亲石元素的浓度(Manning, 1981;Dingwell et al., 1996;Maneta et al.,2015)。London et al.(1989) 研究表明在200 MPa条件下,挥发分可以使成矿流体的固相线降低到450℃,比简单含水花岗岩固相线低200℃以上。Thomas et al.通过对花岗伟晶岩中的包裹体研究表明成矿流体具有超强稀有金属元素溶解能力和强迁移性 (Thomas and Davidson, 2016;Thomas et al.,2019)。

2.2 成矿流体金属含量

通过稀有金属含量的研究可以准确了解花岗伟晶岩型稀有金属矿床成矿机制,但其含量目前仍存在争议,一些学者认为成矿流体中具有极高的金属含量。如Webster et al. (1997)针对德国Ehrenfriedersdorf Sn-W矿床的伟晶岩中石英斑晶开展了熔体包裹体电子探针和二次离子质谱分析,发现熔体包裹体有高含量的Sn、F、P、Li、Rb、Cs、Nb、Ta、Be,部分熔体包裹体中Sn含量高达1000×10-6～2000×10-6,比无矿火成岩高出2个数量级。Thomas et al.(2011b)在德国Ehrenfriedersdorf伟晶岩中的熔体包裹体中测得Be含量为1234×10-6～11025×10-6。Borisova et al.(2012)使用激光烧蚀等离子体四级杆质谱仪对Ehrenfriedersdorf伟晶岩的石英中的富晶体包裹体进行的初步分析研究发现Li 含量可达28400×10-6,Ta含量可达1720×10-6。Thomas et al.基于熔体包裹体的研究表明,成矿流体包含大量的稀有金属元素,具有超强的元素溶解和迁移能力,Be含量高达11600×10-6,Sn含量高达6865×10-6,Ta含量高达3236×10-6,Li含量高达28405×10-6(图4a、4b;Thomas and Davidson, 2016; Thomas et al.,2019)。Li and Chou(2017)运用激光拉曼对甲基卡锂辉石中的富晶体包裹体研究发现,包裹体中的固相矿物含有扎布耶石(Li2CO3)和锂绿泥石,表明成矿流体可能含有极高浓度的Li等金属。

然而,部分学者认为成矿流体中稀有金属含量相对较低。如Morgan and London(1987)通过对加拿大Tanco稀有金属花岗伟晶岩研究,估计成矿流体中含有Li 3523×10-6、Be 169×10-6、Rb 5525×10-6、Cs 3036×10-6、Nb 56×10-6、Ta 300×10-6。London (2014)根据熔体组成为10×10-6Be的花岗岩熔体和矿物-熔体分配系数DBe=0.25,对Be富集进行了瑞利分馏测试计算,结果显示在结晶度达到99.9%时Be丰度仅达到1778×10-6;同时发现在14种最富含绿柱石的伟晶岩中,Be的平均丰度仅为205×10-6。

VCO2—气相二氧化碳;LCO2—液相二氧化碳;LH2O—液相水;S—固相;V—气相a—甲基卡锂辉石中的富晶体包裹体;b—甲基卡锂辉石中富CO2 FIA(Fluid Inclusion Assemblage,流体包裹体组合);c—甲基卡锂辉石中的富晶体FIA;d—大红柳滩锂辉石中的富晶体FIA图3 花岗伟晶岩型锂矿锂辉石中典型的包裹体照片Fig.3 Typical inclusions in the spodumene from the granitic pegmatite-type lithium-bearing ore. (a) Crystal-rich inclusion hosted in the spodumene from Jiajika; (b) CO2-rich Fluid Inclusion Assemblage hosted in the spodumene from Jiajika; (c) Crystal-rich Fluid Inclusion Assemblage hosted in the spodumene from Jiajika; (d) Crystal-rich Fluid Inclusion Assemblage hosted in the spodumene from DahongliutanVCO2-Gas phase carbon dioxide; LCO2-Liquid phase carbon dioxide; LH2O-Liquid phase water; S-Solid phase; V-Gas phase

2.3 成矿流体的黏度

流体的黏度对成岩成矿流体迁移非常重要,流体的黏度越大越难迁移。London(2009)认为成矿流体H2O含量较低,黏度在700℃时大约为105Pa·s,在过冷到450℃时黏度会更大(108Pa·s),这种黏度相当于常温下的沥青。然而,Thomas et al.(2011c)在距离母岩花岗岩很远的地方观察到了只有几厘米宽的伟晶岩(例如挪威Froland伟晶岩),认为高度黏稠、H2O不饱和、过冷、迅速结晶的成矿流体从花岗岩源区长距离(高达10 km)的迁移几乎不可能,与距离母岩较远的地方出现的细脉伟晶岩事实不符。Thomas and Davidson(2012)研究表明,花岗伟晶岩稀有金属矿床成矿流体黏度较低,在任何给定的温度条件下,花岗岩熔体大约比形成的相应伟晶岩熔体的黏度高1万倍。Li and Chou(2017)对甲基卡伟晶岩锂辉石中的富晶体包裹体研究发现,其成分看起来接近于具有低黏度、高元素扩散性和传质能力的含水硅酸盐液体,使得稀有金属元素的浓度极高。

影响成矿流体黏度的因素很多,包括流体组成、温度、压力等。Audétat and Keppler (2004) 通过实验表明流体的黏度只随着硅酸盐含量的增加而缓慢增加。含20%硅酸盐的流体在800℃时的黏度等于纯水在室温下的黏度(10-3Pa·s)。含有50%硅酸盐的流体在相同条件下的黏度为0.8×10-1Pa·s,与橄榄油在室温下的黏度相当。Nabelek et al.(2010)研究表明H2O可显著降低流体黏度。Thomas and Davidson(2012)通过对含水量不同的熔体黏度计算发现,含水越多熔体黏度越低。此外,Sowerby and Keppler(2002)认为大量的F、B、P和碱(Li、Na) 可显著降低流体黏度。

a—Ehrenfriedersdorf花岗岩-伟晶岩系统石英中熔体包裹体中Be质量浓度-H2O浓度图;b—Ehrenfriedersdorf花岗岩-伟晶岩系统石英中熔体包裹体中CA浓度-H 2O浓度图(该图表明在临界条件下,某些元素在超临界流体或熔体中的溶解度非常高);c—假二元硅酸盐熔体-H 2O系统中A型和B型熔体包裹体温度-H2 O浓度图;d—5个不同伟晶岩石英中熔体包裹体的结果绘制的假二元溶线CA代表Be、Sn、As、P、Cl、Ta;CA-crit代表在临界H2O浓度下的CA浓度;H2O-crit代表临界H2O浓度);T C代表临界温度图4 假二元硅酸盐熔体-H 2O系统的温度-H2 O浓度图(Thomas and Davidson,2016)Fig.4 Temperature versus H2O concentration plot of the pseudo-binary silicate melt-H2 O system (Thomas and Davidson,2016).(a) Be concentration versus H2O concentration plot in melt inclusions in the Ehrenfriedersdorf granite-pegmatite system. (b)CAversus H2O concentration plot in melt inclusions hosted in quartz in the Ehrenfriedersdorf granite-pegmatite system. (c)Relationship of type-A and type-B melt in clusions in a temperature versus H2O concentration plot of the pseudo-binary silicate melt-H2 O system. (d) Results of melt inclusions in quartz of five different pegmatites plot a pseudo-binary solvus(The figure b shows that certain elements have very high solubility in supercritical fluids or melts under critical conditions; CA represents Be, Sn, As, P, Cl,Ta; CA-crit represents the concentration of CA at the critical H2O concentration; H2O-critre represents critical H2O concentration;T C represents critical temperature)

2.4 成矿流体的温压条件

成矿流体的温压条件是理解成矿流体演化的重要因素。不同学者运用矿物温度计法、稳定同位素法和流体包裹体温压计法等方法获得的花岗伟晶岩稀有金属矿床成矿温压范围较广,争议较大(图5)。传统观点认为在锂铝硅酸盐相图中,假设伟晶岩形成于平衡条件,锂辉石形成于温度为600℃压力超过350 MPa的高温高压条件(Linnen et al., 2012)。吴长年等(1994) 运用包裹体温度计法对阿尔泰伟晶岩锂辉石中的包裹体研究表明,其均一温度为480～550℃,据富CO2流体包裹体等容线估算的形成压力为200～360 MPa。Kontak et al. (2002)根据已获得的加拿大新斯科舍省Peggys Cove伟晶岩形成压力330～350 MPa,对石英宿主的流体包裹体均一温度进行压力校准获得伟晶岩流体包裹体捕获温度为600～650℃。Ackerman et al. (2007)基于对捷克斯洛伐克Vlastějovice伟晶岩中包裹体等容线计算,结合细晶花岗岩的水饱和固相线和长石温度计重建了伟晶岩型锂矿的结晶温压条件:温度为500～570℃,压力为310～430 MPa。Simmons et al. (2016)运用石榴子石-黑云母温压计法和矽线石、石英、白云母、黑云母和碱性长石矿物组合推断了美国Mt Mica伟晶岩形成的温压条件:温度为630～650℃,压力为300 MPa。Li and Chou(2017)运用热液金刚石压腔(在接近锂辉石形成条件的温度和压力下)对甲基卡锂辉石中富晶体包裹体均一实验发现,均一温度为500～720℃;通过富晶体包裹体均一时热液金刚石压腔水介质的压力估计包裹体内压为300～500 MPa。Li et al.(2019)也通过稀有金属矿物在热液金刚石压腔水溶液中的结晶实验研究表明,锂辉石结晶温度为530～700℃,结晶压力为300～500 MPa。黄永胜等(2016)运用英国Linkam公司的THMSG600冷热台对新疆阿尔泰伟晶岩研究表明,该地区三叠纪伟晶岩早期结构带中的CO2-NaCl-H2O流体包裹体均一温度为400～581℃,利用含CO2包裹体等容线相交法确定流体捕获压力为235～308 MPa。Mulja and Williams-Jones(2018)通过流体包裹体研究结合斜长石-石榴子石-白云母-黑云母地质温压计和接触变质组合的多平衡模型对加拿大魁北克稀有金属花岗伟晶岩研究表明,锂辉石形成温度为450～480℃,压力为350 MPa。Xiong et al.(2019)通过Linkam THMS600冷热台对扎乌龙花岗伟晶岩锂矿中的富晶体包裹体均一实验表明,均一温度为500～580℃,结合PVT相图估算压力为310～480 MPa。

图5 世界主要伟晶岩矿床形成的温压条件Fig.5 Temperature and pressure conditions of main pegmatite deposits in the world

然而,部分学者认为花岗伟晶岩型锂矿形成温度较低。如London(1986)对加拿大Tanco伟晶岩中的富晶体流体包裹体和相平衡研究表明,包裹体均一于375～450℃,根据CO2包裹体和H2O包裹体等容线交点估算压力为260～290 MPa。Chakoumakos and Lumpkin (1990)通过伟晶岩液相线、固相线和矿物组合,包裹体等容线,围岩变质条件,伟晶岩矿物平衡综合评估了美国新墨西哥州Harding伟晶岩的P-T路径,表明其形成温度为300～650℃,压力为300～350 MPa。Morgan VI and London (1999)运用长石温度计法对美国加利福尼亚州Ramona伟晶岩研究表明,其形成于375～425℃。Sirbescu and Nabelek(2003)对美国南达科他州Tin Mountain伟晶岩中流体包裹体研究表明,其均一温度为260～420℃,同时据包裹体等容线估算压力在120～420 MPa之间。

3 花岗伟晶岩型稀有金属矿床成矿机制

3.1 花岗伟晶岩成因

关于花岗伟晶岩型稀有金属矿床的成因存在花岗岩浆高程度结晶分异(图6)和富集稀有金属源区小比例深熔两种流行的观点。一些学者认为花岗伟晶岩型稀有金属矿床成矿流体是花岗岩浆高度分异结晶形成的残余熔体(Norton, 1973;Raimbault et al., 1995;ernand Ercit, 2005;Hulsbosch et al., 2014; London, 2014;吴福元等,2015; Knoll et al., 2018; Roda-Robles et al.,2018)。花岗伟晶岩熔体从源区上升侵位结晶分异过程中,长石中K/Rb比值逐渐降低,钾长石和白云母中的Li、Rb、Cs增高,绿柱石中Na/Li比值降低,但是Cs含量升高,铌铁矿族矿物的Nb/Ta比值升高 (ern, 1989)。 Roda-Robles et al.(2012)通过对西班牙Pinilla de Fermoselle花岗伟晶岩系统的研究也发现在岩浆分异演化过程中,长石、云母和电气石中Li、Be和Sr都有系统性变化,Li和Be随着演化程度的增加而增加,Sr随着演化程度增加而降低;同时发现云母和钾长石中的Rb、Cs和Ba也有系统性变化,Rb和Cs随着演化程度的增加而增加,Ba随着演化程度增加而降低。电气石中的Nb和Ta随着演化程度的增加而增加,Zn则显示相反的变化趋势。

图6 花岗岩和花岗伟晶岩关系示意图(London,2008)Fig.6 Schematic diagram of the relationship between granite and granitic pegmatite (London, 2008)

但是,有学者指出在全球范围内发现的几个伟晶岩矿田和花岗岩侵入体之间缺乏空间、时间和成分的关系。如加拿大Tanco伟晶岩、澳大利亚Greenbush伟晶岩和中国可可托海3号伟晶岩,其成矿母岩花岗岩至今未发现和证实(Lv et al.,2018;张辉等,2019)。另外有学者提出了源区富集的观点,该观点强调地表风化作用导致源区富集稀有金属和B、Li等助熔成分(Stewart, 1978;Kontak et al., 2005; Simmons and Webber, 2008;Melleton et al., 2012; Deveaud et al., 2015; Shaw et al., 2016; Müller et al., 2017; Fuchsloch et al.,2018; Konzett et al., 2018; Gourcerol et al., 2019)。例如,在麻粒岩相变质过程中,区域剪切带可能引导富含稀有金属和助熔元素的流体。这种流体可能会促进地壳深部的部分熔融,从而形成富含稀有金属的熔体(Cuney and Barbey, 2014)。

3.2 稀有金属富集作用

3.2.1 流体不混溶作用

许多学者通过研究发现流体不混溶作用对稀有金属的富集至关重要(Jahns and Burnham,1969; Veksler and Thomas, 2002; Thomas and Davidson, 2016; Fan et al.,2020)。 Warren and Pincus (1940)提出成网组分和变网组分之间为了电荷平衡而竞争非桥氧可引起不混溶作用。传统观点认为伟晶岩是缓慢冷却的结果,而且伟晶岩形成过程中可以出现熔体和流体共存(Jahns and Burnham, 1969)。Roedder (1992)综述了自然界岩浆环境的熔融包裹体,提出岩浆演化过程中存在硅酸盐熔体到富水硅酸盐熔体再到晚期水质流体的转变过程,并且进一步提出不混溶作用可能在岩浆演化过程中普遍存在。Veksler and Thomas(2002) 利用冷封快速淬火高压釜开展了0.1～0.2 GPa压力水饱和条件下的实验研究,发现在高B、P和F(每种质量含量为5%)条件下硅酸盐熔体、低密度溶液和高盐度熔体共存。Veksler(2004)实验研究也证实花岗伟晶岩形成过程存在液态不混溶。Badanina et al.(2004)对俄罗斯Khangilay杂岩体中Li-F稀有金属花岗岩开展了熔体包裹体研究,发现高演化伟晶岩的成分与低演化熔体成分无法用分离结晶作用解释,而可能是晚阶段残留熔体不混溶的结果。在不混溶作用中,铝硅酸盐熔体相富含K,超高盐度的流体相富含Na,这也与钠化花岗岩位于上部的地质事实相吻合。Borisova et al. (2012)认为伟晶岩形成过程中发生了熔体和流体相的不混溶。Thomas et al.发现德国Ehrenfriedersdorf伟晶岩石英中存在成分互补的两种类型包裹体(Thormes et al.,2011a,2019;Thomas and Davidson,2016):A型熔体包裹体(贫H2O-过铝质包裹体)和B型熔体包裹体(富H2O-过碱质包裹体),对其进行加压均一化实验,发现它们沿硅酸盐熔体-H2O体系的假二元溶解曲线分布(图4c),因此,认为它们是由熔体-熔体不混溶产生的共轭熔体。此外,Thomas and Davidson(2016)对5个花岗伟晶岩稀有金属矿床系 统 (Ehrenfriedersdorf、 Zinnwald、 Königshain、Malkhan和Tanco)进行综合比较,表明其都具有相似的熔体-熔体不混溶模式(图4d)。Ballouard et al.(2020)对南非纳马夸兰中元古代奥兰治河伟晶岩带中的稀有金属伟晶岩研究表明,在伟晶岩相关稀有金属矿床的形成过程中,熔体-熔体和熔体-流体不混溶和交代作用起着重要作用。Fan et al.(2020)研究了西藏西昆仑白龙山花岗伟晶岩的全岩和矿物(如白云母和锂辉石)的主量、微量元素和锂同位素组成,表明贫锂和富锂伟晶岩是花岗岩浆演化后期熔体-流体分离的产物。贫锂伟晶岩形成于贫H2O富硅酸盐熔体系统,而富锂伟晶岩产生于富H2O贫硅酸盐熔体(超临界流体)系统。熔体-流体分离可导致显著的锂同位素分馏,7Li富集在强结合的残余硅酸盐熔体中,6Li倾向于流体中较弱的水合键;或者,富含6Li的富锂伟晶岩可能是由6Li在超临界或近临界流体中的优先富集引起的。超临界流体出溶在锂同位素分馏和稀有金属成矿中起着重要作用。

此外,某些自然矿床在构造、地质和地球化学方面也显示出流体不混溶。陈毓川等(2003)对新疆阿尔泰成矿带地质特征和成矿规律进行总结时发现,发生流体不混溶的岩体的围岩都具有良好的封闭作用。李建康(2006)对甲基卡矿区研究发现,甲基卡二云母花岗岩内部的伟晶岩脉与岩体接触关系明显,呈现出贯入式特点,二云母花岗岩和伟晶岩微量和稀土元素组成表现出一定的突变性,指示甲基卡矿床为岩浆液态不混溶成因。廖芝华等(2019)对可尔因二云母花岗岩体和伟晶岩地球化学特征研究,结果显示二云母花岗岩与伟晶岩间存在Al2O3、Na2O和Si2O、K2O相分离的现象;二云母花岗岩与伟晶岩Li、Be、Nb、Ta等稀有元素含量,F、B等挥发分含量,ΣREE、ΣCe/ΣY、(La/Yb)N、(Gd/Yb)N、δEu等具有突变现象。这些特征表明伟晶岩为二云母花岗岩液态不混溶形成。

3.2.2 过冷却条件下的组成带纯化模型

London (1992)提出的新成因模型认为,伟晶岩形成过程未必极其富水,强调过冷却作用和非平衡结晶是伟晶岩重要机制。由于过冷却效应的存在,相对应的富挥发分硅酸盐熔体可以维持低至500℃,甚至350℃。花岗伟晶岩型稀有金属矿床成矿过程中需要稀有金属矿物达到饱和。如Stewart (1978)和Maneta et al. (2015)确定,含水花岗岩熔体在Li2O含量大约为1.5%时会饱和形成透锂长石、锂辉石和锂霞石。然而,达到经济意义的伟晶岩Li2O平均含量接近0.5% (Stewart,1978),Maneta et al. (2015) 用该值模拟伟晶岩中锂铝硅酸盐的结晶,发现当与饱和值比较时,很明显锂铝硅酸盐未达到饱和,不应该在结晶开始就存在,但是在熔体充分过冷却时,锂铝硅酸盐会在熔体结晶开始时结晶。London(2018)总结发现形成稀有金属花岗伟晶岩的熔体中稀有金属的丰度大多低于含水花岗岩系统固相线处的饱和值,因此,形成花岗伟晶岩型矿床的熔体在结晶开始前通过明显的液相线过冷在稀有元素矿物中达到饱和。

London(2014,2018)提出了组成带纯化富集(constitutional zone refining;图7)来解释花岗岩伟晶岩的内部结构和矿物学演化,其认为在熔体过冷到远低于其液相线和固相线温度时,组成带纯化富集作用开始时,在伟晶岩从脉体边缘向核部结晶,结晶前缘会形成一层边界层流体,亲石元素通过边界层流体先结晶,稀有金属物质在边界层流体中逐渐聚集,直至伟晶岩脉结晶最后阶段,稀有金属沉淀形成花岗伟晶岩型稀有金属矿床。边界层中诸如H2O、B、P和F等挥发分的积累增强了元素的横向扩散。富挥发分边界层被认为有利于伟晶岩内部大块晶体的生成。过冷却条件下组成带纯化模型主要表现在伟晶岩组构、矿物生长习性和矿物组合分带等伟晶岩显著特征方面。矿物学分带普遍表现为斜长石主要在外部带,钾长石主要在中间带,石英集中于最后形成的内核,稀有金属矿物集中于中间带和内核。

图7 花岗伟晶岩浆组成带纯化示意图(London,2014,2018)Fig.7 Schematic diagram of the constitutional zone refining of granitic pegmatite magma(London, 2014, 2018). (a) During the constitutional zone refining, the compatible components dissolve from the bulk melt and attach to the surface of the diagenetic mineral through the boundary layer. (b) Because the volatiles decrease the solid temperature and enhance the miscibility of the components,the excluded rare metal components are enriched in the boundary layer liquid. (c) Once the composition of the melt is exhausted, the boundary layer liquid will crystallize, resulting in a mutation in the composition of the growing minerals (mica, tourmaline)

4 甲基卡花岗伟晶岩型锂矿

稀有金属锂是最重要的关键金属之一,被广泛应用于新能源汽车、电池、航空航天、可控核聚变、玻璃陶瓷、医药、润滑脂等领域(Linnen et al.,2012)。随着战略性新兴产业的发展,尤其是新能源汽车需求量的增加,极大地刺激了全球各国对锂资源的需求。中国国务院办公厅印发的《新能源汽车产业发展规划(2021—2035)》提及:“到2025年,中国新能源汽车销售总量达到汽车新车销售总量的20%左右,到2035年,纯电动汽车成为新销售车辆的主流”。作为新能源电动汽车的核心材料,中国对锂的需求量也将随着新能源汽车量的增多而加大。

锂资源主要来源有花岗伟晶岩型锂矿、盐湖卤水型锂矿和沉积岩型锂矿(Kesler et al.,2012;Li et al.,2015;刘丽君等,2017;Benson et al.,2017;Bowell et al.,2020)。各类型锂矿在全球分布广泛(图8),盐湖卤水型锂矿主要分布智利、玻利维亚、阿根廷和中国等国家;花岗伟晶岩型锂矿主要分布在澳大利亚、中国、南非和加拿大等国家;沉积型锂矿全球分布较少,可开采的主要矿物为贾达尔石(羟硼硅钠锂石),主要在塞尔维亚。

1—甲基卡;2—可尔因;3—阿尔泰;4—大红柳滩;5—Zavitskoye;6—Goltsovoer;7—Tastyq;8—Vishnvakovskoe;9—Lakha;10—Ural mining;11—Ullava lantta;12—Minade Barroso;13—Guarda;14—Zinnwald;15—Winneba;16—Manono-Kitolo;17—Kamativi;18—Bikita;19—Cape Cross-Brandberg-Uis;20—Rubicon Mine;21—Greenbushes;22—MountCaitlin;23—Mount Marion;24—Kemerton;25—Tanco;26—Quebec;27—Kings Mtn;28—Aracuai;29—Sao Joaodel Rei;30—扎布耶碳酸盐型盐湖;31—西台吉乃尔硫酸盐型盐湖;32—东台吉乃尔硫酸盐型盐湖;33—Salton sea;34—Silver Peak;35—Searles;36—Uyuni;37—Atacama;38—Jadar;39—Mcdermitt图8 世界主要锂矿床分布图Fig.8 Distribution of major lithium deposits in the world1-Jiajika; 2-Keryin; 3-Altay; 4-Dahongliutan; 5-Zavitskoye; 6-Goltsovoer; 7-Tastyq; 8-Vishnvakovskoe; 9-Lakha; 10-Ural mining; 11-Ullava lantta; 12-Minade Barroso; 13-Guarda; 14-Zinnwald; 15-Winneba; 16-Manono-Kitolo; 17-Kamativi; 18-Bikita; 19-Cape Cross-Brandberg-Uis; 20-Rubicon Mine; 21-Greenbushes; 22-MountCaitlin; 23-Mount Marion; 24-Kemerton;25-Tanco; 26-Quebec; 27-Kings Mtn; 28-Aracuai; 29-Sao Joaodel Rei; 30-Zabuye carbonate-type salt lake; 31-West Taiji’naier sulfate-type salt lake; 32-East Taiji’naier sulfate-type salt lake; 33-Salton sea; 34-Silver Peak; 35-Searles; 36-Uyuni; 37-Atacama; 38-Jadar; 39-Mcdermitt

盐湖卤水型锂矿虽然储量巨大,但其品位低。中国盐湖卤水型锂矿近80%分布在青藏高原的盐湖中(王秋舒,2016),开发条件差,且大多数Mg/Li值较高,提锂技术尚未完全成熟,开发利用困难。沉积岩型锂矿储量较小,品位中等,锂主要赋存于锂蒙脱石等黏土矿物中(Bowell et al.,2020)。同时国内的沉积型锂矿多数与铝土矿、煤伴生,尚没有独立开发利用(于沨等,2019)。花岗伟晶岩型锂矿品位高、易开采,锂主要赋存于锂辉石、透锂长石、锂云母中。中国花岗伟晶岩型锂矿主要分布在松潘-甘孜巨型硬岩型锂矿带和阿尔泰造山带,这是锂资源的最重要来源。目前新疆阿尔泰锂资源面临枯竭问题,松潘-甘孜造山带内的花岗伟晶岩型锂矿将成为国内今后锂资源开发的重要基地。甲基卡花岗伟晶岩型锂矿位于中国西部大型松潘-甘孜锂矿带内(许志琴等,2018),是目前亚洲规模最大的硬岩型锂矿,具有规模大、品位高、矿种多、埋藏浅、选矿性能好等特点(王登红和付小方,2013;付小方等,2014)。该矿是由S型花岗岩高度结晶分异使得锂元素富集的产物(梁斌等,2016;代鸿章等,2018;李贤芳等,2020;Zhang et al.,2021)。因此,文中以甲基卡花岗伟晶岩型锂矿为例,探讨花岗伟晶岩型稀有金属矿床流体成矿机制研究目前存在的问题。

甲基卡花岗伟晶岩型锂矿矿区整体呈片麻岩穹窿构造(许志琴等,2019,2020)。矿区出露地层为三叠系西康群砂页岩经多阶段区域变质和接触变质作用而形成的黑云母石英片岩、二云母石英片岩和红柱石十字石石英片岩等中浅变质岩系(熊欣等,2019)。矿区出露的唯一侵入岩体为马颈子二云母花岗岩岩体,位于矿区南部,该岩体的LA-ICP-MS锆石U-Pb法年龄为223±1 Ma (郝雪峰等,2015),属于印支晚期的产物(熊欣等,2019)。矿区具有一定规模的伟晶岩共509条,其中矿(化)脉占62%(付小方等,2017)。空间上,花岗伟晶岩脉自侵入体接触带向外大致出现微斜长石型伟晶岩带(Ⅰ)→微斜长石-钠长石型伟晶岩带(Ⅱ)→钠长石型伟晶岩带(Ⅲ)→钠长石-锂辉石型伟晶岩带(Ⅳ)→锂(白)云母型伟晶岩带(Ⅴ)→石英脉带等不规则分带,稀有金属在空间上大致具有Be→Li→Nb+Ta→Cs+Rb的分带性变化(付小方等,2017;Huang et al.,2020;图9)。

图9 甲基卡稀有金属矿田地质简图(Huang et al.,2020)Fig.9 Geological sketch of the Jiajika rare metal ore field(Huang et al., 2020)

尽管已对甲基卡锂矿开展了较多的研究,但是甲基卡锂矿的流体成矿机制仍存在如下问题,①与经典伟晶岩型锂矿如加拿大Tanco、新疆可可托海3号脉花岗伟晶岩型锂矿相比,甲基卡花岗伟晶岩型锂矿单个矿脉内部分带不明显(付小方等,2017)。因此,甲基卡花岗伟晶岩型了锂矿被London (2018)归类为矿脉内部无显著分带的伟晶岩矿床。②甲基卡花岗伟晶岩型锂矿发育大量<5 mm的细粒度的锂辉石,被称为细晶锂辉石(付小方等,2017;刘善宝等,2019)。甲基卡花岗伟晶岩型锂矿大量的细晶锂辉石与粗晶、巨晶锂辉石的成因联系不清。③相关学者对甲基卡花岗伟晶岩型锂矿做了较多的成矿流体研究,揭示出甲基卡花岗伟晶岩型锂矿发育丰富的熔体和流体包裹体,是开展岩浆-热液过渡研究的理想对象。已有报道中甲基卡花岗伟晶岩型锂矿锂辉石中的富晶体包裹体均一温度较高(Li and Chou,2017),这与甲基卡花岗伟晶岩锂矿富含挥发分,可能形成于过冷却条件相矛盾。④甲基卡花岗伟晶岩型锂矿成因机制目前仍不明确。一些学者认为甲基卡花岗伟晶岩型锂矿是液态不混溶作用形成(李建康,2006;付小方等,2017)。但是,由于伟晶岩矿床中大量发育次生流体包裹体,流体包裹体与熔体包裹体是否为流体不混溶作用产物仍存在巨大争议。因此,甲基卡这种特殊的花岗伟晶岩型锂矿的成矿控制因素、成矿过程、成矿机制有待深入研究。

5 总结和展望

稀有金属是一类重要的关键金属。花岗伟晶岩是稀有金属赋存的重要岩石类型之一。花岗伟晶岩型稀有金属矿床成矿流体富集挥发分(B、P、F和H2O)和稀有金属元素,这使得该类矿床成矿流体黏度降低,迁移性增强。此外,挥发分可与稀有金属结合形成络合物,增强稀有金属元素的溶解能力。花岗伟晶岩型稀有金属矿床成矿流体形成的温压条件变化较大,可能形成于高温高压条件,也可能形成于过冷却条件。花岗伟晶岩型稀有金属矿床成矿流体形成机制主要有两种,分别为花岗质岩浆高度结晶分异演化和富成矿元素地壳物质小比例深熔。流体不混溶和组成带纯化作用可能促进稀有金属超常富集。

目前对于花岗伟晶岩型稀有金属矿床的相关研究所取得的认识方面还存在不足,需要进一步研究来完善。以下为将来稀有金属伟晶岩型矿床成矿流体研究的几个重点方向:

(1)天然样品在形成后会受到多期热液改造事件,造成了原生包裹体识别困难。阴极发光成像技术可以清晰地识别出不同热液事件的叠加关系,进而可帮助判断包裹体是否为原生包裹体还是后期叠加的次生包裹体。因此,基于阴极发光成像技术的详细显微岩相学是开展伟晶岩流体包裹体研究的重要手段。

(2)花岗伟晶岩型稀有金属矿床矿物捕获的包裹体挥发分含量高,内压较高,常温下进行均一化容易爆裂,因此需要在施加合适外压的热液金刚石压腔中进行。

(3)以往对花岗伟晶岩型稀有金属矿床成矿中的包裹体研究主要是以显微激光拉曼探针定性和半定量研究,或者针对子矿物开展电子探针分析,缺乏开展系统的单个包裹体LA-ICP-MS成分分析。通过对甲基卡花岗伟晶岩型锂矿开展单个包裹体LA-ICP-MS成分分析,可快捷、高效、准确获得多种元素信息,为成矿过程示踪带来一系列关键证据。

(4)目前相关学者主要集中于明显分带的伟晶岩研究,对于分带性不明显的锂辉石伟晶岩矿床研究较少,需要深入研究成矿流体与成矿机制。

致谢:感谢吕古贤教授和张宝林教授约稿撰写本文。感谢四川省地质调查院付小方教授级高级工程师在野外采样中提供的帮助。感谢南京大学许志琴院士、李广伟教授、郑碧海博士和中国科技大学李万才博士在研究中提供的帮助。感谢审稿专家和编辑部提供的指导和修改意见。