岩石学报  2021, Vol. 37 Issue (10): 2923-2943, doi: 10.18654/1000-0569/2021.10.01   PDF    
喜马拉雅淡色花岗岩——关键金属Sn-Cs-Tl的富集机制
高利娥, 曾令森, 严立龙, 高家昊, 赵令浩     
中国地质科学院地质研究所, 北京 100037
摘要: 关键金属是全球高科技产业不可或缺的战略性资源,其富集机制和成矿作用是目前国际矿床学研究的热点之一。我们对喜马拉雅带吉隆和亚东地区淡色花岗岩开展系统的地球化学研究,发现侵入到藏南拆离系的淡色花岗岩含有较高的Sn、Cs、Tl、Be、W、B、Li和Bi。全岩元素地球化学分析表明,这些淡色花岗岩具有如下特征:(1)富集关键金属元素;(2)为原始岩浆经历斜长石、锆石、独居石、磷灰石、云母分离结晶作用后的残余熔体;(3)关键元素的富集和矿化与花岗岩高度分离结晶作用密切相关。随着分异程度的增强,岩浆变为富挥发分的高SiO2体系,关键金属元素在残余熔体中富集,并且最后可能形成具有工业价值的矿床。由于地球化学特征的相似性,Cs和Tl呈类质同象替代钾、铷进入云母中。富集关键金属元素的花岗岩在时间上和空间上属于与藏南拆离系相关的同构造侵位花岗岩,藏南拆离系的活动促使了原始岩浆的广泛分离结晶作用,以及后期的关键金属元素(如Rb、Cs和Tl)的富集。
关键词: 关键金属元素    Sn-Cs-Tl    分离结晶作用    淡色花岗岩    喜马拉雅造山带    
Enrichment mechanisms of Sn-Cs-Tl in the Himalaya leucogranite
GAO LiE, ZENG LingSen, YAN LiLong, GAO JiaHao, ZHAO LingHao     
Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China
Abstract: Critical elements are the indispensable strategic resource for global high-tech industries,and the mechanisms for the enrichment and mineralization for such elements is currently one of the key areas in the international mineral deposit studies. In the Gyirong and Yadong areas along the Himalayan orogen,a systematic geochemical study of leucogranites indicate that the leucogranites which invaded the Southern Tibet Detachment System (STDS) contain a relatively elevated content of Sn,Cs,Tl,Be,W,B,Li and Bi. Whole-rock element geochemical data indicate that (1) these leucogranites are enriched in critical elements; (2) they represent the residual melts,which resulted from fractional crystallization of plagioclase,zircon,monazite,apatite,and mica of more primitive magmas; and (3) enrichment and mineralization of critical elements are closely related to the evolution of such highly fractionated granites. As magmas evolved into high SiO2 and volatile magmas by intensive fractional differentiation,critical elements with a high fluid compatibility are enriched in the residual melt and might eventually form economic ore deposits. Due to geochemical similarity,Cs and Tl are incorporated into mica by replacement of K or Rb. These granites enriched with critical elements are syntectonic granites spatially and temporally associated with the STDS whose activities promote extensive fractional crystallization of the more primitive magmas and later enrichment of critical elements such as Rb,Cs,and Tl.
Key words: Critical element    Sn-Cs-Tl    Fractional crystallization    Leuocogranite    The Himalayan orogen    

关键金属是指现今社会必需的、对战略性新兴产业的发展至关重要的、但安全供应存在较高风险的一类金属元素及其矿床的总称(涂光炽等, 2004; 蒋少涌等, 2020; 侯增谦等, 2020; 翟明国等, 2019),包括三稀金属(稀有、稀土和稀散)和稀贵金属(铂族金属)(毛景文等, 2019a, b ; 王登红, 2019)。关键金属作为全球高科技产业不可或缺的战略性资源,尤其在航空航天、军事、芯片和新能源等领域的重要作用,其富集和成矿机制是目前国际矿床学研究的热点(毛景文等, 2019b)。随着科技的发展,未来对关键金属的需求将迅猛增长。但是,关键金属矿产是国际上最近提出的资源概念,针对此类矿床的专门研究刚刚起步,认知程度较低(侯增谦等, 2020)。

喜马拉雅造山带内分布着两条新生代淡色花岗岩带,部分花岗岩经历了结晶分异作用,形成高硅花岗岩甚至伟晶岩(Wu et al., 2017; Gao et al., 2021; 吴福元等, 2015; 刘志超等, 2020)。随着中国地质调查局成都地调中心李光明研究员团队在拿日雍错(也叫错那洞)发现了Be多金属工业矿体(李光明等, 2017),最近几年以中国科学院地质与地球物理研究所吴福元院士和南京大学王汝成教授团队为代表的学者对喜马拉雅淡色花岗岩稀有金属的成矿作用开展了较详细的研究,在佩枯措、夏如、吉隆、定结、告乌、然巴、康马、错那等地,发现了代表Be成矿的绿柱石,代表Nb-Ta成矿的钽铁矿、烧绿石-细晶石、铌铁金红石,代表Sn成矿的锡石,代表Li成矿的锂云母、锂辉石等(Wang et al., 2017; Wu et al., 2020)。在拿日雍错花岗岩中,相对富集Bi、Cs、Li、Sn、Be、Pb、B、W、Ta等稀有金属成矿元素(高利娥等, 2019)。已有研究认为,藏南喜马拉雅高分异淡色花岗岩形成过程中伴随着稀有金属元素的富集,很可能成为我国新的重要稀有金属成矿带(Wang et al., 2017; Wu et al., 2020; Gao et al., 2021; 李光明等, 2017; 曾令森和高利娥, 2017)。但在喜马拉雅淡色花岗岩稀有金属富集和成矿机理、影响因素等方面有待深入研究,形成新认识,揭示成矿潜力高的勘探靶区。

本文以吉隆和亚东地区淡色花岗岩为研究对象,通过详细的主量元素和微量元素测定,查明这些花岗岩中关键金属元素的异常情况,尤其是Sn-Cs-Tl富集情况,揭示喜马拉雅淡色花岗岩的Sn-Cs-Tl富集行为和机制。

1 地质背景与样品

喜马拉雅造山带由特提斯喜马拉雅带、高喜马拉雅结晶带、低喜马拉雅带和次喜马拉雅带组成(图 1a)。它们之间的界限分别为藏南拆离系、主中央逆冲断层、主边界逆冲断层。在喜马拉雅造山带内,新生代淡色花岗岩沿着特提斯喜马拉雅带和高喜马拉雅带分布。在特提斯喜马拉雅带,新生代花岗岩侵入片麻岩穹窿核部。在高喜马拉雅带中,淡色花岗岩侵入到高喜马拉雅结晶岩系中,或侵入藏南拆离系中,自西向东主要有Mandhata(纳木那尼)、Manaslu(马纳斯鲁)、Gyirong(吉隆)、Mount Qomolangma(珠穆朗玛)、Makalu(马卡鲁)、Dinggye(定结)、Dingga(定日)、Khula Kangri(库拉岗日)和Cuona(错那)等岩体。

图 1 藏南喜马拉雅造山带地质简图(a)和亚东地区地质简图(b) YTS-雅鲁藏布江缝合带;STDS-藏南拆离系;MCT-主中央逆冲推覆带;MBT-主边界逆冲推覆带;MFT-主前缘逆冲推覆带;TH-特提斯喜马拉雅;HH-高喜马拉雅;LH-低喜马拉雅;SH-次喜马拉雅 Fig. 1 Simplified geologic map of the Himalayan orogenic belt, southern Tibet (a) and simplified geological map of the Yadong area (b) YTS-Yarlung-Tsangpo suture; STDS-Southern Tibet Detachment System; MCT-Main Center Thrust; MBT-Main Boundary Thrust; MFT-Main Fronter Thrust; TH-Tethyan Himalaya; HH-High Himalaya; LH-Lower Himalaya; SH-Sub-Himalaya

亚东位于高喜马拉雅带中部,新生代淡色花岗岩出露面积近200km2,根据锆石U-Pb年龄可以解体为三个独立岩体(图 1b, Liu et al., 2017)。告乌淡色花岗岩形成于23Ma,为含电气石的淡色花岗岩,位于藏南拆离系和亚东谷露裂谷的交接部位,以岩床的形式侵位于藏南拆离系和剪切带上盘的特提斯喜马拉雅沉积岩系中,整体无变形,部分表现出弱定性特征,为同构造侵位或构造活动末期侵位岩浆。20~18Ma的顶嘎岩体和17~16Ma淡色花岗岩体,岩性为二云母花岗岩,侵入藏南拆离系和剪切带下盘的高喜马拉雅结晶岩系中,这两期淡色花岗岩切穿藏南拆离系,岩浆侵位应发生于韧性剪切活动停止后(刘志超等, 2020)。

吉隆位于高喜马拉雅带西部,新生代淡色花岗岩侵位于藏南拆离系和剪切带下盘的高喜马拉雅结晶岩系中(图 2)。藏南拆离系由古生代花岗质片麻岩和新生代面理化的淡色花岗岩组成,并被后期的淡色花岗岩体侵位。面理化淡色花岗岩具有明显的同构造性质(图 3a-c),主期侵位约为22~20Ma,并记录了早期36Ma和26Ma的岩浆热液活动(Yang et al., 2009; 高利娥等, 2016)。后期的淡色花岗岩形成于17.7Ma,未变形,岩浆侵入到藏南拆离系中,形成时代在藏南拆离系活动停止之后(高利娥等, 2016)。高喜马拉雅带内,淡色花岗岩形成年龄为22~16Ma(高利娥等, 2016)。其中一条变形的淡色花岗岩脉, 形成于21.0~18.7Ma,记录了早期向北的伸展作用和后期向南的逆冲缩短作用(Wang et al., 2013)。

图 2 吉隆地区地质简图 Fig. 2 Simplified geological map of the Gyirong area

图 3 吉隆地区同构造淡色花岗岩的野外特征 Fig. 3 Field occurrences of the syn-tectonic leucogranites in the Gyirong area

吉隆和亚东地区花岗岩主要由石英、斜长石、钾长石、白云母、电气石、石榴子石(图 3d),以及少量的黑云母、条纹长石、磷灰石、独居石、锆石等组成(图 4a-f)。条纹长石呈脉状、树枝状、网状、蝌蚪状等形态(图 4b)。斜长石具有核边结构,核部含有大量的白云母包裹体,边部干净,核部白云母包裹体呈90°或120°沿着斜长石的解理缝分布(图 4c, d)。电气石半自形-他形柱状,浅黄褐色为主(图 4c, e)。磷灰石呈粒状、卵状等形式成堆出现(图 4f)。手标本上,淡色花岗岩中白云母和黑云母定向排列形成面理(图 3b, c)。未发生塑性变形,表明花岗岩冷凝后期未受强烈变形。亚东地区出现少量的钠长花岗岩,主要由斜长石、石英,以及少量的钾长石、石榴子石、磷灰石等组成(图 4g, h)。

图 4 富集关键金属元素淡色花岗岩(a-f)和钠长花岗岩(g、h)的显微照片 Ap-磷灰石; Bt-黑云母; Gt-石榴子石; Mus-白云母; Pl-斜长石; Qtz-石英; Tour-电气石 Fig. 4 Photomicrographs showing the texture and mineral assemblage of leucogranites enriched with critical metal elements (a-f) and albite granite (g, h) Ap-apatite; Bt-biotite; Gt-garnet; Mus-muscovite; Pl-plagioclase; Qtz-quartz; Tour-tourmaline
2 分析方法

主量及微量元素的测试在自然资源部国家地质实验测试中心进行。主量元素通过XRF(X荧光光谱仪3080E)方法测试,分析精度为5%。微量元素和稀土元素(REE)通过等离子质谱仪(ICP-MS-Excell)分析,含量大于10×10-6的元素的测试精度为5%,而小于10×10-6的元素精度为10%。个别在样品中含量低的元素,测试误差大于10%。31件亚东花岗岩的全岩元素分析结果列在表 1中。

表 1 亚东淡色花岗岩的主量元素(wt%)和微量元素(×10-6)地球化学特征 Table 1 Major element (wt%) and trace element (×10-6) data for the Yadong leucogranites
3 数据及结果

根据岩石组成和关键金属的富集程度,本文数据分为三类:富集关键金属的淡色花岗岩,不富集关键金属的淡色花岗岩和钠长花岗岩。吉隆花岗岩包括马拉山-吉隆的所有淡色花岗岩,数据来自于Gao and Zeng (2014)Gao et al. (2017)。从主量元素含量来看(图 5),大多数花岗岩具有较高的SiO2(69.5%~78.8%),Al2O3(11.0%~17.0%),K2O(2.3%~8.4%),和Na2O(1.4%~6.1%),但较低的CaO(≤2.0%)、FeO(< 2.0%)、MgO、MnO、P2O5和TiO2,较高的A/CNK (>1.1)和K2O/Na2O(≥1.0)比值。与不富集关键金属元素的花岗岩相比,富集关键金属元素的花岗岩具有较高的Na2O(>3.4)和P2O5,但较低的CaO(< 0.9%)、FeO(< 1.4%)、MgO、MnO和TiO2。值得关注的是,亚东地区的两件钠长花岗岩的Na2O含量高达8.4%。

图 5 淡色花岗岩的主量元素特征 Fig. 5 Major element compositions of the leucogranites

在微量元素特征上(图 6),富集关键金属元素的花岗岩具有Ba和Ti显著的负异常,Th和Sr弱负异常。相比较,钠长花岗岩还显示了K和REE的负异常。与不富集关键金属元素的花岗岩相比(Rb=104×10-6~386×10-6,Sr=14×10-6~320×10-6,Rb/Sr=0.6~19.0),富集关键金属元素的花岗岩具有较高的Rb(>300×10-6)和Rb/Sr(>5.0),较低的Sr(< 70×10-6)和Ba,钠长花岗岩具有最低的Rb和Sr,即Rb=71×10-6~77×10-6,Sr=11×10-6~12×10-6,Rb/Sr=5.9~6.6(图 7)。在高场强元素含量和比值上,两类花岗岩差别不大,Zr=3×10-6~165×10-6,Hf=0.2×10-6~10×10-6,Nb=0.2×10-6~35×10-6,Ta=0.1×10-6~27×10-6,Zr/Hf=10.7~40.6,Nb/Ta=1.2~16.3(图 8图 9)。相对比,钠长花岗岩具有较高的Nb和Ta(Nb=35×10-6~50×10-6,Ta=10×10-6~15×10-6),中等含量的Zr和Hf(Zr=17×10-6~20×10-6,Hf=2×10-6)和Nb/Ta比值(3.3~3.4),较低的Zr/Hf比值(8.6~9.5)。

图 6 富集关键金属元素淡色花岗岩的原始地幔标准化微量元素蜘蛛网图(a、c)和球粒陨石标准化稀土元素配分图(b、d)(标准化值据Sun and McDonough, 1989) Fig. 6 Primitive mantle (PM)-normalized trace element (a, c) and chondrite-normalized rare earth element (b, d) distribution patterns for the leucogranites enriched with critical metal elements (normalization values after Sun and McDonough, 1989)

图 7 淡色花岗岩的微量元素特征 Fig. 7 Trace elements concentrations and ratios for the leucogranites (a) Rb-Sr;(b) Ba-Rb/Sr;(c) Eu/Eu*-Sr;(d) TZr-CaO

图 8 淡色花岗岩的其它微量元素特征和岩浆温度 Fig. 8 Additional trace element concentrations, ratios, and trace element-based thermometers for the leucogranites (a) Zr-A/NK;(b) TZr-TLREE;(c) TiO2-Ta/Nb;(d) Ta-Ta/Nb;(e) Nd/Nd*-P2O5;(f) Nd/Nd*-∑LREE

图 9 淡色花岗岩的高场强元素特征 Fig. 9 High field strength element concentrations and ratios for the leucogranites (a) Zr-Zr/Hf;(b) Hf-Zr/Hf;(c) Nb-Nb/Ta;(d) Ta-Nb/Ta

在稀土元素特征上(图 6),富集关键金属元素的花岗岩显示略微富集轻稀土和中稀土,亏损重稀土,(La/Yb)N=2.17~19.82,显著的负Eu异常(Eu/Eu*=0.15~0.66),微弱的负Nd异常(Nd/Nd*=NdN/[(PrN2×SmN)1/3=0.67~0.86)。钠长花岗岩具有轻稀土和中稀土亏损,重稀土略富集,负Eu异常和负Nd异常的特征(Eu/Eu*=0.56~0.80,Nd/Nd*=0.79~0.80)。

富集关键金属元素的花岗岩含有较高的W、Sn、B、Be、Cs、Tl等关键金属元素(图 10图 11)。具体表现在:全部样品富集Sn、Cs、Tl和Be元素,变化范围为1.69×10-6~55.20×10-6、14.6×10-6~126.0×10-6、1.05×10-6~4.97×10-6和6.71×10-6~184×10-6,平均值为18.92×10-6、35.48×10-6、2.13×10-6和16.12×10-6,约是维氏花岗岩平均值的6.3倍、7.1倍、1.4倍和2.9倍;部分样品富集W、B、Li、Bi、U和Pb,最高值可分别达9.76×10-6、1610×10-6、182×10-6、12.2×10-6、32.6×10-6和58.8×10-6,平均值分别为2.12×10-6、200.6×10-6、146.5×10-6、5.51×10-6、9.39×10-6和41.01×10-6;略亏损Nb、Ta元素,变化范围为1.90×10-6~78.80×10-6和0.34×10-6~15.2×10-6,平均为10.93×10-6、2.47×10-6。相比较,钠长花岗岩富集Nb、Ta、Sn、Be、Cs。不富集关键金属元素的花岗岩中这些元素的浓度大多数都低于维氏花岗岩平均值。

图 10 淡色花岗岩的关键金属W (a)、Sn (b)、Nb (c)和Ta (d)与Rb/Sr的关系图解 Fig. 10 Selected rare metal W (a), Sn (b), Nb (c) and Ta (d) plotted against Rb/Sr for the leucogranites

图 11 淡色花岗岩的关键金属B (a)、Be (b)、Cs (c)和Tl (d)与Rb/Sr的关系图解 Fig. 11 Selected rare metal B (a), Be (b), Cs (c) and Tl (d) plotted against Rb/Sr for the leucogranites
4 讨论

从上面地球化学数据可以看出,富集关键金属元素的淡色花岗岩具有:(1)较高的SiO2、Al2O3和A/CNK比值,为过铝质淡色花岗岩(图 5);(2)较高的Rb(>300×10-6)和Rb/Sr(>5.0),较低的Sr,高度变化的高场强元素的浓度和比值(图 7-图 9);(3)略微富集LREE和MREE,亏损HREE,显著的负Eu异常和微弱的负Nd异常(图 6);和(4)富集关键金属元素Sn、Cs、Tl等。与富集关键金属元素的淡色花岗岩相比,钠长花岗岩具有不一致的地球化学特征, 包括(1)Na2O含量高达8.4%;(2)非常低的Rb和Sr含量(Rb < 80×10-6,Sr < 15×10-6),较高的Nb和Ta含量(Nb>35×10-6,Ta>10×10-6),Nb/Ta比值小于5,Zr/Hf比值小于10;(3)LREE和MREE亏损,HREE略富集,更加显著的Eu和Nd负异常;和(4)富集关键金属元素Nb、Ta、Sn、Be、Cs等。

4.1 富集关键金属淡色花岗岩的形成过程:分离结晶作用

随着岩浆分离结晶作用的进行,当岩浆演变成高硅体系,矿物组成和元素地球化学行为会发生改变,如:(1)主要造岩矿物相溶解度的改变(Ren, 2012);(2)微量元素分配系数的变化(Shearer and Papike, 1987; Glazner et al., 2008; Gao et al., 2021)以及稀土元素四分组效应的出现(Liu et al., 2016);(3)副矿物溶解行为的变化(Wolf and London, 1994; Bea, 1996; Linnen and Keppler, 1997, 2002; Gao et al., 2021a, b ; 郭春丽等, 2017; 曾令森和高利娥, 2017);(3)钠长花岗岩的出现(Liu et al., 2014; Wu et al., 2020);和(4)结晶富Mn石榴子石、高Hf锆石、电气石和富集稀有金属元素的绿柱石、铌铁矿等(Liu et al., 2019; 高利娥等, 2012; 刘志超等, 2020)。

图 5-图 9显示了富集关键金属元素淡色花岗岩中重要主微量元素的变化趋势和相互关系。在花岗质岩浆分离结晶过着中,随着温度的降低,斜长石中的CaO含量降低(Shearer and Papike, 1987; Scaillet et al., 1995; Patiño Douce and Harris, 1998)。与不富集关键金属元素的花岗岩相比,富集关键金属元素的花岗岩具有较高的Na2O(图 5b),但较低的CaO(图 5c)和锆饱和岩浆温度(图 8b),CaO含量与岩浆温度呈正相关关系(图 7d),并出现钠长花岗岩。另外,富集关键金属元素花岗岩还具有较低的Sr、Ba,和显著的Eu负异常,Eu异常与Sr含量正相关(图 7c),以上特征表明富集关键金属元素的淡色花岗岩是较原始花岗质岩浆经历不同程度斜长石分离结晶作用的产物。

Zr和Hf具有相似的地球化学性质,岩浆演化过程中Zr/Hf比值应保持一致。越来越多的文献数据表明,在岩浆演化过程中,Zr/Hf比值可发生变化(Dupuy et al., 1992; David et al., 2000; Linnen and Keppler, 2002)。在SiO2>74.0%的高硅花岗岩中,较低Zr/Hf比值的形成机制有两种:(1)锆石的分离结晶作用(Bea, 1996; David et al., 2000; Claiborne et al., 2006);(2)花岗岩熔体结构的变化(Linnen and Keppler, 2002; 郭春丽等, 2017; 曾令森和高利娥, 2017)。在锆石中,Zr的溶解度是Hf的五倍(Linnen and Keppler, 2002),因此锆石的分离结晶作用可以降低花岗质熔体的Zr/Hf比,但同时降低残余熔体的Zr和Hf的含量。在富集关键金属元素的花岗岩中,Zr与A/NK呈正相关(图 8a),表明锆石发生了分离结晶作用。但是随着Zr/Hf比值降低,Zr和Hf含量先降低,后升高(图 9a, b)。当熔体的Zr/Hf比值超过20时,Zr-Zr/Hf和Hf-Zr/Hf系统关系发生了改变(Gao et al., 2021a)。随着岩浆分异作用的进行,花岗质岩浆的结构将改变,锆石和铪石的溶解度发生变化,Hf的溶解度将明显高于Zr的溶解度(Linnen and Keppler, 2002),导致熔体的Zr/Hf比值降低。因此,锆石的分离结晶作用和熔体结构的变化共同控制了花岗岩中Zr和Hf的含量和系统关系。

Nb和Ta也表现出相似的地球化学性质,在多数地质作用中一般不发生明显的分馏。最近一系列研究得出,Nb和Ta易溶于云母和富Ti矿物相中(Stepanov and Hermann, 2013; Stepanov et al., 2014; Gao et al., 2017)。在金红石和钛铁矿中DTa>DNb,而在云母中DTa < DNb(Klemme et al., 2005; Stepanov et al., 2014)。在富集关键金属元素的花岗岩中,TiO2-Ta/Nb和Ta-Ta/Nb的相关性(图 8c, d)表明花岗岩经历了云母的分离结晶作用。与Zr-Hf相似,当熔体的Nb/Ta比值超过5时,Nb-Nb/Ta和Ta-Nb/Ta系统关系发生了改变(图 9c, d),也与花岗质岩浆的结构变化相关。

另外,富集关键金属元素的花岗岩中具有明显的Nd负异常(图 6)。在淡色花岗岩中,稀土元素主要赋存于独居石和磷灰石中(Whittington and Treloar, 2002; Zeng et al., 2005)。富集关键金属元素的花岗岩中,∑LREE和Nd异常正相关(图 8f),P2O5和Nd异常(图 8e)也正相关,表明花岗岩经历了独居石和磷灰石的分离结晶作用。

以上分析可以表明,富集关键金属元素的花岗岩经历了斜长石、锆石、独居石、磷灰石、云母等的分离结晶作用,虽然具有高的10000Ga/Al值,但落入分离结晶作用形成的S型花岗岩区域(图 12),实质为高分异的S型花岗岩。

图 12 淡色花岗岩的A型花岗岩判别图(据Whalen et al., 1987) Fig. 12 Ga/Al against K2O+Na2O (a), FeOT/MgO (b), Nb (c) and Zr (d) for the leucogranites (after Whalen et al., 1987)
4.2 喜马拉雅淡色花岗岩中关键金属Sn-Cs-Tl富集机制初探

三稀金属包括稀有、稀土和稀散,是关键金属的重要组成部分。稀有金属包括锂(Li)、铍(Be)、铷(Rb)、铯(Cs)、铌(Nb)、钽(Ta)、锆(Zr)、铪(Hf)、钨(W)、锡(Sn)等。钨锡是我国传统的优势矿产资源,蒋少涌等(2020)指出我国原生钨锡矿床有5种主要类型:斑岩型、云英岩型、矽卡岩型、石英脉型和锡石硫化物型;3种次要类型:蚀变花岗岩型、热液角砾岩型和低温热液脉型。花岗岩的高度分异作用对钨锡成矿至关重要,与钨矿化有关的岩浆岩主要为陆壳重熔型花岗岩。与锡矿床相关的岩浆岩除了高分异花岗岩外,还包括A型花岗岩等多种岩石类型(Chen et al., 2013; 华仁民等, 2010; 蒋少涌等, 2020)。铯的地壳丰度为3×10-6,地球化学性质与钾和铊相似,常以类质同像形式存在于富钾矿物中,如天河石(含铷微斜长石)中Cs2O可达0.1%~0.3%。铯大多呈分散态存在,含铯矿物常见于花岗伟晶岩中,随着岩浆分异程度增强,铯含量增加,容易富集于岩浆晚期阶段、碱性岩和花岗伟晶岩中,与绿柱石、锂辉石、锂云母、电气石以及铌钽矿物共生。此外,铯还可富集在盐湖卤水中。所以,铯矿床可分为碱性长石花岗岩(花岗伟晶岩型)和卤水型。

稀散金属一般指在地壳中丰度很低(多为10-9级),分布极为分散的元素,包括镓(Ga)、锗(Ge)、硒(Se)、镉(Cd)、铟(In)、碲(Te)、铼(Re)、铊(Tl)8种元素。据最新统计资料表明,尽管中国镓、锗、硒、镉和铟的资源储量在全球第一,但是,稀散金属矿床的成矿成因等方面研究程度相对较低(谢桂青等, 2020)。铊的地壳丰度很低(0.45×10-6),地壳岩石中的w(Tl)和w(Pb)正相关,铊主要呈类质同象替代钾、铷进入云母和钾长石(Rader et al., 2018)。一直以来铊主要作为伴生组分从铜矿床和铅锌硫化物矿床中作为综合利用对象而被回收(温汉捷等, 2020)。目前为止铊成矿的研究相对比较薄弱,温汉捷等(2020)大致划分出两个具有工业意义的矿床类型, 即低温热液型铊矿床和块状硫化物型含铊矿床。热液型铊矿床一般和金、汞等矿化关系密切,铊的超常富集主要出现在成矿热液演化较晚阶段,最大的特点是能够形成铊的独立矿床(陈代演和邹振西, 2000)。硫化物型含铊矿床中铊以类质同象进入黄铁矿晶格, 不产出铊独立矿物(Murao and Itoh, 1992),如中国安徽和县香泉发现铊矿。

关键金属元素通常含量极低,在常具有“稀、细、伴”的特征,独立矿物不易被直接观察,并且,地球化学行为较为复杂,矿化类型和成矿机制多样,认知难度大(侯增谦等, 2020)。毛景文等(2019a)初步把关键矿产分为8种成因类型,即:(1)与花岗岩-伟晶岩有关的W、Sn、Nb、Ta、Li、Rb、Cs、Be矿床;(2)与碳酸岩-碱性岩有关的REE、Nb、U矿床;(3)与镁铁质-超镁铁质岩有关的Ni、Co、Cr、Pt族元素矿床;(4)低温热液型Tl、Te矿床;(5)多种类型热液矿床中的伴生组分,包括Re、Ge、In、Cd、Tl、Te、Se、Sc、Ga等;(6)表生沉积型矿床;(7)与表生盐湖有关的Li、Rb、Cs矿床;(8)与表生风化作用有关的稀土矿床。

地球化学数据表明,本文中富集关键金属元素的花岗岩富集Sn、Cs、Tl等元素。这些花岗岩经历了斜长石、锆石、独居石、磷灰石、云母等的分离结晶作用,为高分异的S型花岗岩。根据目前关键矿产的成因类型划分,喜马拉雅淡色花岗岩属于与花岗岩-伟晶岩有关的W、Sn、Nb、Ta、Li、Rb、Cs、Be矿床。实验研究表明,REE、F、Li、Rb、Cs、Be、Sn和W在分离结晶过程中更倾向于进入熔体中(Keppler and Wyllie, 1991; Linnen, 1998)。随着分离结晶作用的增强,花岗岩中的钾长石向微斜长石甚至天河石(即富Rb钾长石)演化,斜长石含量减少且逐渐向富Na端元演化,最后形成钠长花岗岩(Wu et al., 2003a, b )。本文中富集关键金属元素的花岗岩中K-Rb-Cs-Tl具有很好的正相关性(图 13),表明Cs和Tl主要呈类质同象替代钾、铷进入云母。同时,钨锡的富集与花岗岩的高度分异演化密切相关。上述分析表明,随着岩浆分异程度增强,Sn、Cs、Tl、Be、W、B、Li、Bi、U、Pb等关键金属元素在残余熔体中富集,Cs和Tl呈类质同象替代钾、铷进入云母。

图 13 淡色花岗岩的K-Rb-Cs-Tl的关系图解 Fig. 13 Diagrams showing K-Rb-Cs-Tl relationship for the leucogranites

喜马拉雅造山带被认为是稀有金属成矿的有利地区,错那洞发现了较大规模的钨锡铍矿床(李光明等, 2017)。在错那洞穹窿内,花岗质岩浆发生高度的分离结晶作用(高利娥等, 2017, 2019),使熔体中Be、W、Sn以及挥发分富集,一方面,在伟晶岩中形成绿柱石;另一方面,从岩浆中出溶的成矿热液侵入碳酸盐岩,发生水岩反应,沉淀出Be、W、Sn等稀有金属元素(何畅通等, 2020)。所以,错那洞稀有金属的成矿类型包括两类:(1)岩浆伟晶岩型;和(2)热液矽卡岩型。在岩浆分离结晶过程中,由于稀有金属元素在硅酸盐矿物和熔体之间的分配系数小于1,有利于稀有金属元素在熔体中富集,并为后期热液矽卡岩型矿化创造有利条件(Evensen and London, 2002)。因此,喜马拉雅造山带不仅是稀有金属成矿的有利地区,而且可以富集其它关键金属元素,岩浆结晶分离作用是最重要的成因机制。

4.3 关键金属富集的动力学背景

从地质图上看(图 1图 2),这些富集关键金属元素的淡色花岗岩都侵位于藏南拆离系内,吉隆淡色花岗岩形成于20~22Ma(高利娥等, 2016),亚东淡色花岗岩形成于23Ma(Liu et al., 2017),为同构造侵位花岗岩(Yang et al., 2009; 刘志超等, 2020)。Wu et al. (2020)提出,喜马拉雅地壳厚度为70km,藏南拆离系侵角为17.3±0.3°,当岩浆沿着藏南拆离系上升,需要运移150~200km。如此长距离的移动,促使了岩浆的广泛分离结晶作用。研究表明,变形作用会驱动压滤效应,促使岩浆分异(Bea et al., 2005)。因此,侵位过程中,富集关键金属元素的淡色花岗岩受到藏南拆离系的韧性剪切影响,促使岩浆发生广泛的结晶分异作用,以及后期的关键金属元素富集。

5 结论

(1) 吉隆和亚东淡色花岗岩可以分成三类:富集关键金属的淡色花岗岩、不富集关键金属的淡色花岗岩和钠长花岗岩。

(2) 全岩主微量元素地球化学特征表明,富集关键金属元素的花岗岩经历了斜长石、锆石、独居石、磷灰石、云母等的分离结晶作用,属于高分异的S型花岗岩。

(3) 富集关键金属元素的花岗岩含有较高的Sn、Cs、Tl、Be、W、B、Li、Bi等关键金属元素,属于与花岗岩-伟晶岩有关的矿床。

(4) 关键金属元素的富集与花岗岩的高度分异演化密切相关,随着岩浆分异程度增强,Sn、Cs、Tl、Be、W、B、Li、Bi、U、Pb等关键金属元素在熔体中富集,Cs和Tl呈类质同象替代钾、铷进入云母。

(5) 富集关键金属元素的花岗岩侵位于藏南拆离系内,为同构造侵位花岗岩。藏南拆离系的活动促使了岩浆的广泛分离结晶作用,以及后期的关键金属元素富集。

参考文献
Bau M. 1996. Controls on the fractionation of isovalent trace elements in magmatic and aqueous systems: Evidence from Y/Ho, Zr/Hf, and lanthanide tetrad effect. Contributions to Mineralogy and Petrology, 123(3): 323-333 DOI:10.1007/s004100050159
Bea F. 1996. Controls on the trace element composition of crustal melts. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 87(1-2): 33-41 DOI:10.1017/S0263593300006453
Bea F, Fershtater GB, Montero P, Smirnov VN and Molina JF. 2005. Deformation-driven differentiation of granitic magma: The Stepninsk pluton of the Uralides, Russia. Lithos, 81(1-4): 209-233 DOI:10.1016/j.lithos.2004.10.004
Chen DY and Zou ZX. 2000. Studying situation about Lanmuchang-type thallium (mercury) ore deposits in southwestern Guizhou. Guizhou Geology, 17(4): 236-241 (in Chinese with English abstract)
Chen J, Wang RC, Zhu JC, Lu JJ and Ma DS. 2013. Multiple-aged granitoids and related tungsten-tin mineralization in the Nanling Range, South China. Science China (Earth Sciences), 56(12): 2045-2055 DOI:10.1007/s11430-013-4736-9
Claiborne LL, Miller CF, Walker BA, Wooden JL, Mazdab FK and Bea F. 2006. Tracking magmatic processes through Zr/Hf ratios in rocks and Hf and Ti zoning in zircons: An example from the Spirit Mountain batholith, Nevada. Mineralogical Magazine, 70(5): 517-543 DOI:10.1180/0026461067050348
David K, Schiano P and Allègre CJ. 2000. Assessment of the Zr/Hf fractionation in oceanic basalts and continental materials during petrogenetic processes. Earth and Planetary Science Letters, 178(3-4): 285-301 DOI:10.1016/S0012-821X(00)00088-1
Dupuy C, Liotard JM and Dostal J. 1992. Zr/Hf fractionation in intraplate basaltic rocks: Carbonate metasomatism in the mantle source. Geochimica et Cosmochimica Acta, 56(6): 2417-2423 DOI:10.1016/0016-7037(92)90198-R
Evensen JM and London D. 2002. Experimental silicate mineral/melt partition coefficients for beryllium and the crustal Be cycle from migmatite to pegmatite. Geochimica rt Cosmochimica Acta, 66(12): 2239-2265 DOI:10.1016/S0016-7037(02)00889-X
Gao LE, Zeng LS, Shi WG, Chen ZY, Hu MY and Sun DY. 2012. Two types of garnets in the Cenozoic granites from the Himayalan Orogenic Belt: Geochemical characteristics and implications for crustal anatexis. Acta Petrologica Sinica, 28(9): 2963-2980 (in Chinese with English abstract)
Gao LE and Zeng LS. 2014. Fluxed melting of metapelite and the formation of Miocene high-CaO two-mica granites in the Malashan gneiss dome, southern Tibet. Geochimica et Cosmochimica Acta, 130: 136-155 DOI:10.1016/j.gca.2014.01.003
Gao LE, Zeng LS, Wang L, Hou KJ, Gao JH and Shang Z. 2016. Timing of different crustal partial melting in the Himalayan Orogenic Belt and its tectonic implications. Acta Geologica Sinica, 90(11): 3039-3059 (in Chinese with English abstract)
Gao LE, Zeng LS and Asimow PD. 2017. Contrasting geochemical signatures of fluid-absent versus fluid-fluxed melting of muscovite in metasedimentary sources: The Himalayan leucogranites. Geology, 45(1): 39-42 DOI:10.1130/G38336.1
Gao LE, Gao JH, Zhao LH, Hou KJ and Tang SH. 2017. The Miocene leucogranite in the Nariyongcuo gneiss dome, southern Tibet: Products from melting metapelite and fractional crystallization. Acta Petrologica Sinica, 33(8): 2395-2411 (in Chinese with English abstract)
Gao LE, Zeng LS, Hu GY, Gao JH, Zhao LH and Wang YY. 2019. Rare metal enrichment in leucogranite within Nariyongcuo gneiss dome, South Tibet. Earth Science, 44(6): 1860-1875 (in Chinese with English abstract)
Gao LE, Zeng LS, Zhao LH, Hou KJ, Guo CL, Gao JH and Wang YY. 2021a. Geochemical behavior of rare metals and high field strength elements during granitic magma differentiation: A record from the Borong and Malashan Gneiss Domes, Tethyan Himalaya, southern Tibet. Lithos, 398-399: 106344 DOI:10.1016/j.lithos.2021.106344
Gao LE, Zeng LS, Zhao LH, Gao JH and Shang Z. 2021b. Behavior of apatite in granitic melts derived from partial melting of muscovite of metasedimentary sources. China Geology, doi: 10.31035/cg2021009
Glazner AF, Coleman DS and Bartley JM. 2008. The tenuous connection between high-silica rhyolites and granodiorite plutons. Geology, 36(2): 183-186 DOI:10.1130/G24496A.1
Guo CL, Zeng LS, Gao LE, Su HZ, Ma XH and Yin B. 2017. Highly fractionated grantitic minerals and whole-rock geochemistry prospecting markers in Hetian, Fujian Province. Acta Geologica Sinica, 91(8): 1796-1817 (in Chinese with English abstract)
He CT, Qin KZ, Li JX, Zhou QF, Zhao JX and Li GM. 2020. Preliminary study on occurrence status of beryllium and genetic mechanism in Cuonadong tungsten-tin-beryllium deposit, eastern Himalaya. Acta Petrologica Sinica, 36(12): 3593-3606 (in Chinese with English abstract) DOI:10.18654/1000-0569/2020.12.03
Hou ZQ, Chen J and Zhai MG. 2020. Current status and frontiers of research on critical mineral resources. Science Bulletin, 65: 3651-3652 (in Chinese)
Hua RM, Li GL, Zhang WL, Hu DQ, Chen PR, Chen WF and Wang XD. 2010. A tentative discussion on differences between large-scale tungsten and tin mineralizations in South China. Mineral Deposits, 29(1): 9-23 (in Chinese with English abstract)
Jiang SY, Zhao KD, Jiang H, Su HM, Xiong SF, Xiong YQ, Xu YM, Zhang W and Zhu LY. 2020. Spatiotemporal distribution, geological characteristics and metallogenic mechanism of tungsten and tin deposits in China: An overview. Science Bulletin, 65(33): 3730-3745 (in Chinese)
Keppler H and Wyllie PJ. 1991. Partitioning of Cu, Sn, Mo, W, U, and Th between melt and aqueous fluid in the systems haplogranite-H2O-HCl and haplogranite-H2O-HF. Contributions to Mineralogy and Petrology, 109(2): 139-150 DOI:10.1007/BF00306474
Klemme S, Prowatke S, Hametner K and Günther D. 2005. Partitioning of trace elements between rutile and silicate melts: Implications for subduction zones. Geochimica et Cosmochimica Acta, 69(9): 2361-2371 DOI:10.1016/j.gca.2004.11.015
Li GM, Zhang LK, Jiao YJ, Xia XB, Dong SL, Fu JG, Liang W, Zhang Z, Wu JY, Dong L and Huang Y. 2017. First discovery and implications of Cuonadong superlarge Be-W-Sn polymetallic deposit in Himalayan metallogenic belt, southern Tibet. Mineral Deposits, 36(4): 1003-1008 (in Chinese with English abstract)
Linnen RL and Keppler H. 1997. Columbite solubility in granitic melts: Consequences for the enrichment and fractionation of Nb and Ta in the Earth's crust. Contributions to Mineralogy and Petrology, 128(2): 213-227
Linnen RL. 1998. The solubility of Nb-Ta-Zr-Hf-W in granitic melts with Li and Li+F: Constraints for mineralization in rare metal granites and pegmatites. Economic Geology, 93(7): 1013-1025 DOI:10.2113/gsecongeo.93.7.1013
Linnen RL and Keppler H. 2002. Melt composition control of Zr/Hf fractionation in magmatic processes. Geochimica et Cosmochimica Acta, 66(18): 3293-3301 DOI:10.1016/S0016-7037(02)00924-9
Liu ZC, Wu FY, Ji WQ, Wang JG and Liu CZ. 2014. Petrogenesis of the Ramba leucogranite in the Tethyan Himalaya and constraints on the channel flow model. Lithos, 208-209: 118-136 DOI:10.1016/j.lithos.2014.08.022
Liu ZC, Wu FY, Ding L, Liu XC, Wang JG and Ji WQ. 2016. Highly fractionated Late Eocene (~35Ma) leucogranite in the Xiaru Dome, Tethyan Himalaya, South Tibet. Lithos, 240-243: 337-354 DOI:10.1016/j.lithos.2015.11.026
Liu ZC, Wu FY, Qiu ZL, Wang JG, Liu XC, Ji WQ and Liu CZ. 2017. Leucogranite geochronological constraints on the termination of the South Tibetan Detachment in eastern Himalaya. Tectonophysics, 721: 106-122 DOI:10.1016/j.tecto.2017.08.019
Liu ZC, Wu FY, Liu XC, Wang JG, Yin R, Qiu ZL, Ji WQ and Yang L. 2019. Mineralogical evidence for fractionation processes in the Himalayan leucogranites of the Ramba Dome, southern Tibet. Lithos, 340-341: 71-86 DOI:10.1016/j.lithos.2019.05.004
Liu ZC, Wu FY, Liu XC and Wang JG. 2020. The mechanisms of fractional crystallization for the Himalayan leucogranites. Acta Petrologica Sinica, 36(12): 3551-3571 (in Chinese with English abstract) DOI:10.18654/1000-0569/2020.12.01
Mao JW, Yuan SD, Xie GQ, Song SW, Zhou Q, Gao YB, Liu X, Fu XF, Cao J, Zeng ZL, Li TG and Fan XY. 2019a. New advances on metallogenic studies and exploration on critical minerals of China in 21st century. Mineral Deposits, 38(5): 935-969 (in Chinese with English abstract)
Mao JW, Yang ZX, Xie GQ, Yuan SD and Zhou ZH. 2019b. Critical minerals: International trends and thinking. Mineral Deposits, 38(4): 689-698 (in Chinese with English abstract)
Murao S and Itoh S. 1992. High thallium content in Kuroko-type ore. Journal of Geochemical Exploration, 43(3): 223-231 DOI:10.1016/0375-6742(92)90108-K
Patiño Douce AE and Harris N. 1998. Experimental constraints on Himalayan Anatexis. Journal of Petrology, 39(4): 689-710 DOI:10.1093/petroj/39.4.689
Rader ST, Mazdab FK and Barton MD. 2018. Mineralogical thallium geochemistry and isotope variations from igneous, metamorphic, and metasomatic systems. Geochimica et Cosmochimica Acta, 243: 42-65 DOI:10.1016/j.gca.2018.09.019
Ren MH. 2004. Partitioning of Sr, Ba, Rb, Y, and LREE between alkali feldspar and peraluminous silicic magma. American Mineralogist, 89(8-9): 1290-1303 DOI:10.2138/am-2004-8-918
Scaillet B, Pichavant M and Roux J. 1995. Experimental crystallization of leucogranite magmas. Journal of Petrology, 36(3): 663-705 DOI:10.1093/petrology/36.3.663
Shearer CK, Papike JJ and Laul JC. 1987. Mineralogical and chemical evolution of a rare-element granite-pegmatite system: Harney Peak Granite, Black Hills, South Dakota. Geochimica et Cosmochimica Acta, 51(3): 473-486 DOI:10.1016/0016-7037(87)90062-7
Stepanov A, Mavrogenes JA, Meffre S and Davidson P. 2014. The key role of mica during igneous concentration of tantalum. Contributions to Mineralogy and Petrology, 167: 1009 DOI:10.1007/s00410-014-1009-3
Stepanov AS and Hermann J. 2013. Fractionation of Nb and Ta by biotite and phengite: Implications for the "missing Nb paradox". Geology, 41(3): 303-306 DOI:10.1130/G33781.1
Sun SS and McDonough WF. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. In: Saunders AD and Norry MJ (eds. ). Magmatism in the Ocean Basins. Geological Society, London, Special Publications, 42(1): 313-345
Tu GC, Gao ZM, Hu RZ, Zhang Q, Li CY, Zhao ZH and Zhang BG. 2004. The Geochemistry and Ore-Forming Mechanism of the Dispersed Elements. Beijing: Geological Publishing House, 1-424 (in Chinese with English abstract)
Wang DH. 2019. Study on critical mineral resources: Significance of research, determination of types, attributes of resources, progress of prospecting, problems of utilization, and direction of exploitation. Acta Geologica Sinica, 93(6): 1189-1209 (in Chinese with English abstract)
Wang RC, Wu FY, Xie L, Liu XC, Wang JM, Yang L, Lai W and Liu C. 2017. A preliminary study of rare-metal mineralization in the Himalayan leucogranite belts, South Tibet. Science China (Earth Sciences), 60(9): 1655-1663 DOI:10.1007/s11430-017-9075-8
Wang XX, Zhang JJ, Liu J, Yan SY and Wang JM. 2013. Middle-Miocene transformation of tectonic regime in the Himalayan orogen. Chinese Science Bulletin, 58(1): 108-117 DOI:10.1007/s11434-012-5414-6
Wen HJ, Zhu CW, Du SJ, Fan Y and Luo CG. 2020. Gallium (Ga), germanium (Ge), thallium (Tl) and cadmium (Cd) resources in China. Science Bulletin, 65(33): 3688-3699 (in Chinese)
Whalen JB, Currie KL and Chappell BW. 1987. A-type granites: Geochemical characteristics, discrimination and petrogenesis. Contributions to Mineralogy and Petrology, 95(4): 407-419 DOI:10.1007/BF00402202
Whittington AG and Treloar PJ. 2002. Crustal anatexis and its relation to the exhumation of collisional orogenic belts, with particular reference to the Himalaya. Mineralogical Magazine, 66(1): 53-91 DOI:10.1180/0026461026610015
Wolf MB and London D. 1994. Apatite dissolution into peraluminous haplogranitic melts: An experimental study of solubilities and mechanisms. Geochimica et Cosmochimica Acta, 58(19): 4127-4145 DOI:10.1016/0016-7037(94)90269-0
Wu FY, Jahn BM, Wilde SA, Lo CH, Yui TF, Lin Q, Ge WC and Sun DY. 2003a. Highly fractionated I-type granites in NE China (Ⅰ): Geochronology and petrogenesis. Lithos, 66(3-4): 241-273 DOI:10.1016/S0024-4937(02)00222-0
Wu FY, Jahn BM, Wilde SA, Lo CH, Yui TF, Lin Q, Ge WC and Sun DY. 2003b. Highly fractionated I-type granites in NE China (Ⅱ): Isotopic geochemistry and implications for crustal growth in the Phanerozoic. Lithos, 67(3-4): 191-204 DOI:10.1016/S0024-4937(03)00015-X
Wu FY, Liu ZC, Liu XC and Ji WQ. 2015. Himalayan leucogranite: Petrogenesis and implications to orogenesis and plateau uplift. Acta Petrologica Sinica, 31(1): 1-36 (in Chinese with English abstract)
Wu FY, Liu XC, Ji WQ, Wang JM and Yang L. 2017. Highly fractionated granites: Recognition and research. Science China (Earth Sciences), 60(7): 1201-1219 DOI:10.1007/s11430-016-5139-1
Wu FY, Liu XC, Liu ZC, Wang RC, Xie L, Wang JM, Ji WQ, Yang L, Liu C, Khanal GP and He SX. 2020. Highly fractionated Himalayan leucogranites and associated rare-metal mineralization. Lithos, 352-353: 105319 DOI:10.1016/j.lithos.2019.105319
Xie GQ, Han YX and Li XH. 2019. A preliminary study of characteristics of dispersed metal-bearing deposits in Middle-Lower Yangtze River metallogenic belt. Mineral Deposits, 38(4): 729-738 (in Chinese with English abstract)
Xie GQ, Li XH, Han YX, Zhu QQ, Li W, Ye H and Song SW. 2020. Recent progress in study of enrichment mechanism of tellurium, selenium and thallium from oxidized gold-rich porphyry-skarn deposits. Mineral Deposits, 39(4): 559-567 (in Chinese with English abstract)
Yang XY, Zhang JJ, Qi GW, Wang DC, Guo L, Li PY and Liu J. 2009. Structure and deformation around the Gyirong basin, North Himalaya, and onset of the South Tibetan Detachment System. Science in China (Series D), 52(8): 1046-1058 DOI:10.1007/s11430-009-0111-2
Zeng LS, Saleeby JB and Asimow P. 2005. Nd isotope disequilibrium during crustal anatexis: A record from the Goat Ranch migmatite complex, southern Sierra Nevada batholith, California. Geology, 33(1): 53-56 DOI:10.1130/G20831.1
Zeng LS and Gao LE. 2017. Cenozoic crustal anatexis and the leucogranites in the Himalayan collisional orogenic belt. Acta Petrologica Sinica, 33(5): 1420-1444 (in Chinese with English abstract)
Zhai MG, Wu FY, Hu RZ, Jiang SY, Li WC, Wang RC, Wang DH, Qi T, Qin KZ and Wen HJ. 2019. Critical metal mineral resources: Current research status and scientific issues. Bulletin of National Natural Science Foundation of China, 33(2): 106-111 (in Chinese with English abstract)
陈代演, 邹振西. 2000. 贵州西南部滥木厂式铊(汞)矿床研究. 贵州地质, 17(4): 236-241. DOI:10.3969/j.issn.1000-5943.2000.04.003
高利娥, 曾令森, 石卫刚, 陈振宇, 胡明月, 孙东阳. 2012. 喜马拉雅造山带新生代花岗岩中两类石榴石的地球化学特征及其在地壳深熔作用中的意义. 岩石学报, 28(9): 2963-2980.
高利娥, 曾令森, 王莉, 侯可军, 高家昊, 尚振. 2016. 喜马拉雅碰撞造山带不同类型部分熔融作用的时限及其构造动力学意义. 地质学报, 90(11): 3039-3059. DOI:10.3969/j.issn.0001-5717.2016.11.006
高利娥, 高家昊, 赵令浩, 侯可军, 唐索寒. 2017. 藏南拿日雍错片麻岩穹窿中新世淡色花岗岩的形成过程: 变泥质岩部分熔融与分离结晶作用. 岩石学报, 33(8): 2395-2411.
高利娥, 曾令森, 胡古月, 高家昊, 赵令浩, 王亚莹. 2019. 藏南拿日雍措片麻岩穹窿淡色花岗岩稀有金属的富集. 地球科学, 44(6): 1860-1875.
郭春丽, 曾令森, 高利娥, 苏红中, 马星华, 尹冰. 2017. 福建河田高分异花岗岩的矿物和全岩地球化学找矿标志研究. 地质学报, 91(8): 1796-1817. DOI:10.3969/j.issn.0001-5717.2017.08.010
何畅通, 秦克章, 李金祥, 周起凤, 赵俊兴, 李光明. 2020. 喜马拉雅东段错那洞钨-锡-铍矿床中铍的赋存状态及成因机制初探. 岩石学报, 36(12): 3593-3606.
侯增谦, 陈骏, 翟明国. 2020. 战略性关键矿产研究现状与科学前沿. 科学通报, 65(33): 3651-3652.
华仁民, 李光来, 张文兰, 胡东泉, 陈培荣, 陈卫锋, 王旭东. 2010. 华南钨和锡大规模成矿作用的差异及其原因初探. 矿床地质, 29(1): 9-23. DOI:10.3969/j.issn.0258-7106.2010.01.003
蒋少涌, 赵葵东, 姜海, 苏慧敏, 熊索菲, 熊伊曲, 徐耀明, 章伟, 朱律运. 2020. 中国钨锡矿床时空分布规律、地质特征与成矿机制研究进展. 科学通报, 65(33): 3730-3745.
李光明, 张林奎, 焦彦杰, 夏祥标, 董随亮, 付建刚, 梁维, 张志, 吴建阳, 董磊, 黄勇. 2017. 西藏喜马拉雅成矿带错那洞超大型铍锡钨多金属矿床的发现及意义. 矿床地质, 36(4): 1003-1008.
刘志超, 吴福元, 刘小驰, 王建刚. 2020. 喜马拉雅淡色花岗岩结晶分异机制概述. 岩石学报, 36(12): 3551-3571. DOI:10.18654/1000-0569/2020.12.01
毛景文, 袁顺达, 谢桂青, 宋世伟, 周琦, 高永宝, 刘翔, 付小方, 曹晶, 曾载淋, 李通国, 樊锡银. 2019a. 21世纪以来中国关键金属矿产找矿勘查与研究新进展. 矿床地质, 38(5): 935-969.
毛景文, 杨宗喜, 谢桂青, 袁顺达, 周振华. 2019b. 关键矿产——国际动向与思考. 矿床地质, 38(4): 689-698.
涂光炽, 高振敏, 胡瑞忠, 张乾, 李朝阳, 赵振华, 张宝贵. 2004. 分散元素地球化学及成矿机制. 北京: 地质出版社, 1-424.
王登红. 2019. 关键矿产的研究意义、矿种厘定、资源属性、找矿进展、存在问题及主攻方向. 地质学报, 93(6): 1189-1209. DOI:10.3969/j.issn.0001-5717.2019.06.003
温汉捷, 朱传威, 杜胜江, 范裕, 罗重光. 2020. 中国镓锗铊镉资源. 科学通报, 65(33): 3688-3699.
吴福元, 刘志超, 刘小驰, 纪伟强. 2015. 喜马拉雅淡色花岗岩. 岩石学报, 31(1): 1-36.
谢桂青, 韩颖霄, 李新昊. 2019. 长江中下游成矿带含稀散金属矿床特征初探. 矿床地质, 38(4): 729-738.
谢桂青, 李新昊, 韩颖霄, 朱乔乔, 李伟, 叶晖, 宋世伟. 2020. 氧化性富金斑岩-矽卡岩矿床中碲、硒、铊富集机制的研究进展. 矿床地质, 39(4): 559-567.
曾令森, 高利娥. 2017. 喜马拉雅碰撞造山带新生代地壳深熔作用与淡色花岗岩. 岩石学报, 33(5): 1420-1444.
翟明国, 吴福元, 胡瑞忠, 蒋少涌, 李文昌, 王汝成, 王登红, 齐涛, 秦克章, 温汉捷. 2019. 战略性关键金属矿产资源: 现状与问题. 中国科学基金, 33(2): 106-111.