江西香炉山钨矿床矽卡岩矿物和成矿花岗岩地球化学特征及其指示意义*

赵文 张怀瑾

中国地质大学资源学院，武汉 430074

长江中下游成矿带是我国华南地区重要的多金属成矿带，以广泛分布的中生代岩浆热液型铜金钼铁矿床而闻名(Xieetal.,2011;Maoetal.,2013)。近年来，在该成矿带南侧相继发现了朱溪、大湖塘等大型-超大型钨矿床(陈国华等,2012;项新葵等,2012)，加上已往发现的北东向分布的一系列钨多金属矿床，共同构成了一条在南岭之外的新的重要钨矿带(Zhaoetal.,2017,2021)。前人对于南岭成矿带已经有了较深入的研究，取得了包括“五层楼+地下室”勘查模型和“九龙脑模式”等一系列理论成果(王登红等,2010;赵正等,2017a,b;Zhaoetal.,2018b;尹政等,2021)。而对于长江中下游钨矿床的研究主要集中于朱溪、大湖塘、阳储岭等中大型矿床(刘善宝等,2017;李宏伟等,2021)，显示主要成矿时间为早燕山期150～135Ma(Zhaoetal.,2017,2021)。

香炉山钨矿是长江中下游南侧钨(铜钼)矿带中一个大型矽卡岩型钨矿床，赋存达29.6万吨的WO3(田邦生和袁步云,2008;张家菁等,2008;吴胜华等,2014)。前人针对该矿床的地质特征、成矿年代学、矿物学和流体包裹体等方面进行了探讨(陈耿炎,1990;田邦生和袁步云,2008;张家菁等,2008;刘勇和周贤旭,2010;陈波和周贤旭,2012;吴胜华等,2014;Daietal.,2018;Wuetal.,2019;Lietal.,2020;吴胜华等,2020)。这些成果为该矿床的进一步研究及区域成矿规律奠定了一定基础。研究结果显示，香炉山钨矿形成时代为晚燕山期(～127Ma)(陈耿炎,1990;张家菁等,2008;Zhaoetal.,2017;Daietal.,2018;Lietal.,2020)，远远年轻于该区域其他典型钨矿床。然而，对其矽卡岩特征与成矿过程仍缺乏详细探讨。Wuetal.(2019)、Lietal.(2020)和吴胜华等(2020)通过对香炉山矿床中的白钨矿、锆石和硫化物微量元素和同位素的研究，显示出其相较于该带早燕山期大湖塘钨矿和南岭的早燕山期柿竹园钨矿，具有相对单一的硫化物组成和花岗岩低氧逸度的特征，指示其具有相对独特的成矿过程和背景。因此，对该矿床特征及花岗岩地球化学特征更进一步的刻画和限定，将有助于长江中下游成矿带的总体认识。

矽卡岩钨矿床作为最重要的钨矿床类型之一，通常分为氧化型矽卡岩钨矿床和还原型矽卡岩钨矿床(Newberry and Einaudi,1981;Newberry and Layer,1998)。一般来讲，还原型矽卡岩钨矿床相较于氧化型矽卡岩钨矿床规模更大(Meinertetal.,2005)。因此，限定矽卡岩型钨矿的类型，将有助于对该矿床的成矿机理研究和资源潜力评估。氧化型和还原型矽卡岩钨矿床分别在矽卡岩矿物组合和成分，以及成矿花岗岩特征上均具有显著的差别(Newberry and Layer,1998;Meinertetal.,2005)。本文将通过研究香炉山钨矿床中矽卡岩矿物和成矿花岗岩的成分特征，限定香炉山钨矿的矽卡岩氧化-还原类型，结合前人研究成果，为成矿流体及物质来源和花岗岩地球化学特征提供指示，同时为区域成矿规律和成矿预测提供依据。

1 区域地质背景

香炉山钨矿位于华南扬子板块北部长江中下游成矿带的西南端，北邻秦岭-大别造山带和华北板块(图1)。该成矿带以斑岩-矽卡岩矿床为主，包含北侧的铜金钼铁矿带(Xieetal.,2011;Maoetal.,2013)，和南侧的钨(铜钼)矿带(Maoetal.,2013;Zhaoetal.,2017,2021)。该地区基底为一套中元古界双桥山群浅变质岩，上覆新元古界-志留系碎屑岩、中泥盆统-下三叠统碳酸盐岩、中三叠统-下侏罗统陆源碎屑岩和早白垩世火山岩(Maoetal.,2013)。该地区主要出露新元古代九岭花岗岩和众多中生代花岗岩。九岭花岗岩呈巨大岩基展布，其岩性主要为黑云母花岗闪长岩、黑云母花岗岩和二云母花岗岩，其年龄为～830Ma(钟玉芳等,2005)。中生代花岗岩侵位于九岭花岗岩及区域地层中，其岩性主要为黑云母花岗岩、花岗闪长岩和花岗斑岩，其年龄为158～118Ma(袁媛等,2012)。中生代花岗岩与区域钨多金属成矿作用密切相关(Maoetal.,2013;Zhaoetal.,2017,2021)。区域的基底构造以南北向的褶皱系和近东西向构造系为主(章泽军等,2003)，盖层构造以北东东向和北北东向褶皱和断裂构造系为主(黄修保,2001)。中生代花岗岩和矿体主要分布于北北东向和北东东向断裂的复合部位(丰成友等,2012)。

2 矿床地质特征

香炉山矿区地处九江坳陷与九岭隆起的交界地带，位于香炉山-观音堂背斜西南倾伏端(图2)。矿区出露地层主要为中寒武统杨柳岗组含炭质灰岩和上寒武统华严寺组深灰色条带状泥质灰岩。香炉山地区的花岗岩出露于矿区东北部，进一步由西向东分为任家山岩体、高湖岩体和太阳山岩体，但三者在深部为一整体并无明显区分(吴胜华等,2014)。

香炉山花岗岩呈细粒到中粗粒结构，主要由石英(55%～60%)、钾长石(～20%)、斜长石(10%～15%)和黑云母(5%～10%)组成，副矿物以钛铁矿为主，为黑云母花岗岩(Daietal.,2018)。香炉山白钨矿体主要赋存于杨柳岗组与香炉山花岗岩接触带部位，矿化带长度>4km，宽1.5km(田邦生和袁步云,2008)。矽卡岩矿体呈似层状(图3a)与地层产状近乎平行，倾角较缓。矿体整体呈北东向展布，向接触带两侧尖灭，以背斜核部为界，分别北西和南东向倾斜。由接触中心向外矽卡岩化逐渐减弱。在外接触带中岩石常呈现细粒变晶结构或条纹条带状构造，并主要形成了辉石、黑云母、长石、石榴石、红柱石、符山石、绿帘石、阳起石、透闪石等矽卡岩矿物；而在内接触带中，热液变质作用主要表现为云英岩化、硅化、绿泥石化、萤石化和高岭土化等(田邦生和袁步云,2008)。其中矽卡岩化和云英岩化是钨矿成矿阶段主要蚀变，而硅化、绿泥石化、萤石化等蚀变则是后期石英-硫化物阶段的典型产物(田邦生和袁步云,2008)。在接触带外侧还发育一些扁豆状和透镜状矿体(吴胜华等,2014)。

矽卡岩型矿石呈浸染状-细脉状构造(图3b)，矽卡岩矿物主要为石榴石和辉石，矿石矿物主要为白钨矿和黄铜矿(图3c)。石榴石和辉石呈嵌晶结构。石榴石矽卡岩中，石榴石呈一级灰干涉色并见双晶(图3d)，主要为钙铝榴石。辉石矽卡岩中，辉石呈半自形到他形粒状，空隙充填萤石(图3e)。白钨矿呈半自形或他形粒状为主，石榴石和石英矿物缝隙，部分与黄铜矿和磁黄铁矿共生(图3f)。整个成矿过程可分为两个成矿期，包括矽卡岩期和石英硫化物期。矽卡岩期可分成三个成矿阶段，主要形成各种钙、铁、镁、铝的硅酸盐矿物。在早期矽卡岩阶段主要生成以岛状和链状的无水硅酸盐。似层状白钨矿主要在这一阶段中的钙矽卡岩中产出。晚期矽卡岩阶段主要特征为形成的角闪石、绿帘石、阳起石、透闪石等含水硅酸盐矿物交代早期矽卡岩矿物。最后的氧化物阶段介于矽卡岩期与石英硫化物期之间，具有过渡性质。在这一阶段可形成长石类以及云母类矿物。在这一阶段生成的白钨矿呈网脉状穿插在早期矽卡岩矿物中。在石英硫化物期中二氧化硅独立的形成大量的石英，并伴有大量金属硫化物的形成，可划分为两个成矿阶段。在早期硫化物阶段中主要生成黄铜矿、黄铁矿、毒砂等铁铜硫化物。而晚期硫化物阶段中主要形成碳酸盐矿物与方铅矿、闪锌矿等铅锌硫化物。

3 样品和测试方法

本次测试样品均为新鲜岩石，采自香炉山钨矿矿区。其中香炉山花岗岩样品6件，制作成粉末后，用于全岩主微量元素分析。矽卡岩样品14件，制作成光片后，分别在光学显微镜和电子显微镜下进行观察和矿物鉴别，随后选择进行原位主微量分析和原位硫同位素分析。

3.1 全岩主微量分析

香炉山花岗岩全岩主微量元素分析在广州澳实分析检测中心完成。全岩主量元素分析仪器使用PANalytical PW2424 X射线荧光光谱仪。主量各元素分析谱线均为Kα，数据校正采用理论α系数法，测试相对标准偏差(RSD)<5%。全岩微量元素分析仪器使用Agilent 7900电感耦合等离子体质谱。样品测试前加入硼酸锂(LiBO2/Li2B4O7)熔剂，混合均匀，在熔炉中于1025℃熔融。待熔融液冷却后，用硝酸、盐酸和氢氟酸消解并定容。

3.2 矿物原位主微量分析

石榴石、辉石和白钨矿原位主量元素分析在香港大学地球科学系完成，仪器型号为日本电子JXA8230，测试加速电压为15kV，电流20nA，束斑大小为1μm。白钨矿原位微量元素分析在南京聚谱检测科技有限公司完成，仪器型号为Agilent 7700X ICP-MS连接Photon Machines Excite 193nm激光剥蚀系统。激光剥蚀束斑大小为25μm。NIST 610作为外标，GSE-1G、BHVO-2G、BCR-2G和NIST 612作为监控标样。Ca作为白钨矿分析内标元素，数据处理时用相应点位的电子探针数据作校正。

3.3 矿物原位硫同位素分析

硫化物原位硫同位素分析在芬兰地质调查局完成，仪器型号为Nu Plasma HR MC-ICP-MS连接Photon Machine Analyte G2激光剥蚀系统。激光分析束斑直径为50μm，频率为5Hz。氦气作为载气运送剥蚀物质(Mülleretal.,2009)。空白背景信号采集时间为20s，样品剥蚀成分信号采集时间为60s。校正标样为黄铁矿国际标样PPP-1(Gilbertetal.,2014)，同时内部标样Py2作为监控标样。具体实验步骤和分析流程参见Wongetal.(2017)。

4 测试结果

4.1 岩石全岩主微量成分

表1 香炉山花岗岩主量元素(wt%)和微量元素(×10-6)成分Table 1 Major element (wt%)and trace element (×10-6)compositions of the Xianglushan granite

香炉山花岗岩稀土元素总含量(∑REE)为187×10-6～237×10-6，其中轻稀土 (LREE)为133×10-6～175×10-6，重稀土(HREE)为54.2×10-6～68.0×10-6，轻重稀土比值(LREE/HREE)为2.2～2.8，和(La/Yb)N=5.4～7.4。稀土元素总体表现为轻稀土富集。球粒陨石标准化稀土元素配分图呈现右倾形式(图4b)，δEu为0.17～0.20，具有强负异常；δCe为0.94～1.04，异常不明显。香炉山花岗岩E-MORB标准化微量元素蛛网图为右倾形式(图4a)，呈现出Rb、U等大离子亲石元素相对富集，Nb、Zr等高场强元素相对亏损。

4.2 矿物原位主微量成分

石榴石、辉石和白钨矿原位主量成分分析结果见表2。石榴石的SiO2、Al2O3和FeO较高且变化不大，分别为37.25%～38.64%、18.06%～19.98%和6.72%～9.39%，而CaO和MnO含量变化较大，分别为21.35%～34.02%和1.30%～11.27%。TiO2和MgO的含量较低，部分低于电子探针检测限，最高值分别为0.75%和0.06%。石榴石端元组分成分跨度较大，其中钙铝榴石为54%～81%，钙铁榴石为5%～14%，以及锰铝-铁铝榴石为5%～40%(图5a)，其中矽卡岩晚期石榴石往往呈现富锰铝-铁铝榴石组分的特征，接近于锰质矽卡岩。

辉石的SiO2和CaO含量较高且变化较小，分别为49.22%～52.21%和22.15%～22.91%，而FeO和MgO含量变化较大，分别为15.73%～26.06%和0.50%～8.18%。MnO和Al2O3的含量较低，分别为1.35%～3.95%和0.12%～0.53%。TiO2的含量可低于电子探针检测限，最高仅可达0.05%。辉石端元成分主要为钙铁辉石(50%～86%)，其次为透辉石(3%～46%)和锰钙辉石(4%～13%)(图5b)，其中早期矽卡岩辉石相对于晚期矽卡岩辉石有相对富透辉石组分的趋势，并且单颗粒辉石中有时呈现核部富透辉石组分和边部富钙铁辉石组分的特征。

白钨矿的WO3和CaO含量较高且变化不大，分别为79.83%～81.16%和18.55%～19.90%。SiO2和FeO的含量较低，分别为0.10%～0.29%和0.02%～0.04%。MgO、MnO和Al2O3的含量部分低于电子探针检测限，最高值分别为0.04%、0.02%和0.12%，可能受显微矿物包裹体的影响。

白钨矿原位微量元素成分分析结果见表3。白钨矿的Mo含量较低，为31.2×10-6～453×10-6，同时含有一定量的Ga(0.07×10-6～1.73×10-6)、Sr(13.4×10-6～46.5×10-6)、Nb(12.3×10-6～90.5×10-6)、Ta(0.11×10-6～2.74×10-6)和Pb(1.01×10-6～7.62×10-6)。Na、Cu、Zn、Rb、Sn、Th和U部分低于ICP-MS检测限，但最高值分别可达61.1×10-6、0.61×10-6、2.87×10-6、1.87×10-6、0.80×10-6、29.5×10-6和59.7×10-6。白钨矿稀土元素的含量较高且变化范围大，为66×10-6～1001×10-6，其在球粒陨石标准化图中呈现右倾的特征，且具有显著的Eu正异常(图6a)。根据Zhaoetal.(2018a)获得的矽卡岩白钨矿晶格应变模型(Lattice Strain Model)参数(r0=1.072Å,EM=4565kbar)，以及吴胜华等(2014)通过流体包裹体获得的温度(T=500K)，通过晶格应变模型变形方程：log(xiD0/Di)=logxi+[(4πNAloge)/R]×(EM/T)×[r0(ri-r0)2/2+(ri-r0)3/3]，获得白钨矿沉淀时成矿热液流体中稀土元素配分形式(图6b)。如图6b所示，呈现出右倾且重稀土相对于轻稀土升高的特征。

表3 香炉山矽卡岩矿床白钨矿LA-ICP-MS分析结果(×10-6)Table 3 LA-ICP-MS results of scheelite from the Xianglushan deposit (×10-6)

4.3 硫化物原位硫同位素成分

黄铜矿和磁黄铁矿原位硫同位素成分分析结果见表4。黄铜矿的δ34SCDT值为1.0‰～4.2‰，磁黄铁矿的δ34SCDT值为3.0‰～5.2‰。可见磁黄铁矿相较于黄铜矿富集重硫同位素。

表4 香炉山矽卡岩矿床黄铜矿和磁黄铁矿原位硫同位素分析结果Table 4 In-situ S isotopic compositions of chalcopyrite and pyrrhotite from Xianglushan

5 讨论

5.1 香炉山花岗岩和矽卡岩类型

矽卡岩型钨矿的成矿岩体通常具有过铝质的特征(Meinertetal.,2005)。香炉山花岗岩成分在TAS图解(图7a)中，投点于高K2O+Na2O和高SiO2花岗岩类范围，具有亚碱性的特征。同时，花岗岩成分在K2O-SiO2图解(图7b)中，投点于钾玄岩系列范围。香炉山花岗岩铝饱和指数为1.31～1.35，属于强过铝质(图7c)。香炉山花岗岩成分特征与全球典型矽卡岩型钨矿成矿花岗岩特征接近(Meinertetal.,2005)。在Zr+Nb+Ce+Y-(K2O+Na2O)/CaO图解(图7d)中，香炉山花岗岩成分落入分异的S型和分异的I型重合区域。此外，香炉山花岗岩成分在A/NK-A/CNK图解中，投点也位于S型花岗岩一侧。在La/Sm-La图解(图8a)中，香炉山花岗岩成分显示岩石形成过程和分离结晶过程密切相关。因此，香炉山花岗岩属于分异的S型花岗岩(Sylvester,1998)，这也与华南其他与钨成矿相关花岗岩的成分特征相一致(图7;Zhaoetal.,2021)。

根据矽卡岩矿物特征，矽卡岩型钨矿被划分为氧化型和还原型：氧化型的矽卡岩型钨矿一般以富钙铁榴石为特征，而还原型的矽卡岩型钨矿则主要以矽卡岩早期富钙铁辉石及矽卡岩晚期富脉状贫Mo白钨矿和富锰铝-铁铝榴石为特征(Newberry and Einaudi,1981;Newberry and Layer,1998)。香炉山矽卡岩中，石榴石和辉石的成分在三角图中均落入典型的矽卡岩型钨矿区域(图5)。其中，辉石主要以钙铁辉石为主，而石榴石以钙铝榴石为主并有少部分锰铝-铁铝榴石。白钨矿主要以大颗粒出现在矽卡岩晚期脉中，呈现贫Mo的特征(表3)。因此，香炉山钨矿属于典型的还原型矽卡岩钨矿。矽卡岩型钨矿的氧化还原特征与成矿岩体和围岩的氧化还原状态密切相关(Newberry and Layer,1998)。前人对香炉山花岗岩的研究显示，其呈现出低氧逸度(fO2)的特征(Lietal.,2020)，这也与岩体中赋存钛铁矿的特征相符，指示其还原的特征。同时，香炉山矽卡岩的赋存围岩杨柳岗组主要为含炭质灰岩(田邦生和袁步云,2008)，同样具有还原特征。因此，可以推断香炉山钨矿的还原型特征是花岗岩和围岩性质共同作用的结果，使其成为大型矽卡岩型钨矿床。

5.2 成矿流体及物质来源

香炉山矽卡岩中晚期黄铜矿和磁黄铁矿的δ34SCDT值分别为1.0‰～4.2‰和3.0‰～5.2‰，均呈现出典型的岩浆硫同位素的特征(Bowman,1998)。此外，通过晶格应变模型及白钨矿稀土数据计算得出的后期流体配分模式，显示出与花岗岩类似的右倾配分模式(图4a、图6b)，指示继承的特征。这表明，香炉山晚期成矿流体仍然以岩浆流体为主。前人针对香炉山矽卡岩不同阶段的流体包裹体研究也显示，从早期到晚期的成矿流体均主要来自于岩浆，随着温度的降低和物质组成的变化形成不同阶段的矽卡岩和矿化(吴胜华等,2014)。因此，岩浆源为成矿流体和物质的主要来源。

香炉山花岗岩为S型花岗岩，多为沉积岩经部分熔融后形成的产物(Chappell and White,1974)。同源岩浆的Nb/Ta比值一般相同(Foley,1984;Barthetal.,2000)，可用于指示部分熔融过程中外来物质的加入。香炉山花岗岩的Nb/Ta比值为6.7～11.4，变化范围较宽，但远小于地幔平均值，而与地壳平均值接近(Wedepohl,1995)，暗示壳源物质为原始岩浆的物源，可能有极少量的幔源物质加入。这与锆石Hf同位素(εHf(t)=-5.7～-3.1，Daietal.,2018；εHf(t)=-6.9～-4.1，Lietal.,2020)和辉钼矿Re元素(Re=12.12×10-6～22.77×10-6，Daietal.,2018)特征，所获得的结论一致。在Rb/Ba-Rb/Sr图解(图8b)中，香炉山花岗岩成分与大湖塘和朱溪成矿花岗岩投点一致落入富粘土原岩区域，与计算的泥质来源熔体较为接近。Romer and Kroner (2016)的研究也显示，钨矿的成矿花岗岩源区与强烈风化沉积物有关。成矿物质在岩浆源区已经有预富集，通过岩浆演化分异再次富集；成矿热液流体从岩浆中的析出，又一次富集成矿物质；最后，经过温度降低及与碳酸质围岩的反应，形成矽卡岩钨矿。区域上，中元古界双桥山群浅变质岩具有较高的钨含量(陈波和周贤旭,2012)，可能为香炉山花岗岩和该区域其他钨成矿花岗岩的物源。

5.3 成矿地球动力学背景

前人针对香炉山钨矿床的年代学研究显示，其形成时代为早白垩世(～127Ma)(陈耿炎,1990;张家菁等,2008;Zhaoetal.,2017;Daietal.,2018;Lietal.,2020)。在白垩纪，香炉山钨矿所在的长江中下游成矿带中的岩浆作用和成矿作用，主要受控于古太平洋板块的俯冲，体现为145～136Ma和136～120Ma两阶段的成矿特征(Sunetal.,2012;Wuetal.,2012;Maoetal.,2013;Zhaoetal.,2017,2021)。对于该时期的成矿动力学背景仍存在争议。Maoetal.(2013)和Zhaoetal.(2021)认为古太平洋板片(伊扎纳吉板片)在135Ma前为斜向大陆边缘俯冲，之后平行大陆边缘俯冲。而Zhaoetal.(2017)和Lingetal.(2009)认为长江中下游成矿带受控于古太平洋板块的洋脊俯冲，并且在～136Ma后，南侧太平洋板片单侧南东向回撤(Zhaoetal.,2017)。长江中下游成矿带成矿规律呈非对称展布(图1)。在该成矿带北侧的铜金钼铁矿带中，分布有一系列的早白垩世火山岩盆地，以及A型花岗岩(Maoetal.,2013;Xueetal.,2015)，指示拉伸环境，及强烈地幔物质的参与。而在成矿带南侧的钨(铜钼)矿带中，包括香炉山花岗岩在内，阳储岭矿床、大湖塘矿床和朱溪矿床的成矿花岗岩均显示出地壳物质部分熔融，及少部分地幔物质加入的特征(Huang and Jiang,2014;李岩等,2014;Maoetal.,2015;李宏伟等,2021)。此外，整个华南在136Ma之后形成北东向展布的拉张盆地及向俯冲边界逐渐年轻的火成岩形成特征(Zhouetal.,2006;Li and Li,2007)。因此，136Ma之后，华南板块下部很可能是伴随着太平洋板片的单向回撤。香炉山花岗岩成分在Rb-(Y+Nb)图解(图9a)和Nb-Y图解(图9b)中，投点也均接近于板内花岗岩区域(Pearceetal.,1984)，并与全球典型矽卡岩型钨矿成矿花岗岩成分投点接近(Meinertetal.,2005)。值得注意的是，长江中下游成矿带中钨成矿相关花岗岩自早燕山期到晚燕山期投点从火山弧-碰撞型花岗岩区域转变为板内花岗岩区域(图9)，这也伴随着花岗岩氧逸度由高向低的转变(Lietal.,2020)，进一步指示了136Ma之后地球动力学背景的转变。综上，香炉山钨矿形成于板内的拉张环境，可能与古太平洋板块中伊扎纳吉板片与太平洋板片间的洋脊俯冲，及太平洋板片回撤有关。

6 结论

香炉山钨矿成矿花岗岩为黑云母花岗岩，具有强过铝质特征，属于分异S型花岗岩。香炉山矽卡岩矿物主要为钙铁辉石为主，而石榴石以钙铝榴石为主并有少部分锰铝-铁铝榴石，同时主要的白钨矿分布于矽卡岩晚期网脉中呈贫Mo的特征。因此，香炉山钨矿属于典型的还原型矽卡岩钨矿。还原型的岩浆热液流体和还原型的炭质围岩共同造成了还原型的矽卡岩特征。香炉山花岗岩的物质来源主要为壳源富粘土物质，可能伴有极少量的幔源物质加入。香炉山花岗岩具有板内花岗岩的特征，指示形成于板内的拉张环境，这可能与古太平洋板块中伊扎纳吉板片与太平洋板片间的洋脊俯冲，及太平洋板片回撤有关。

致谢感谢江西省修水香炉山钨业有限公司为矿区野外工作提供的便利。感谢付晓和黄锦洪在实验测试方面提供的全力协助。感谢周美夫教授、高剑峰研究员、陈伟研究员、赵正研究员和刘泽瑞教授在不同阶段对华南钨锡矿床系列研究的大力支持。同时感谢主编及匿名审稿人对论文提出的宝贵修改意见。