黄层 袁顺达 赵盼捞
钨锡是洁净能源,是信息产业和国防工业等领域不可或缺的战略性金属,其成矿作用和找矿勘查是近年来国际矿床学研究的焦点(毛景文等, 2019; 翟明国等, 2019; 蒋少涌等, 2020; 袁顺达等, 2020)。我国的南岭成矿带是全球最重要的钨锡成矿带之一(Lehmann, 1990; 袁顺达, 2017; Yuanetal., 2019),其钨和锡的储量分别超过全球探明储量的50%和20% (Sunetal., 2012)。前人对南岭钨锡成矿带内钨锡矿床开展了大量的研究,并对成矿带内巨量钨锡超常富集成矿的关键控制因素提出了一系列重要认识,其中主要包括:(1)南岭地区基底为富含钨锡等成矿金属的变沉积岩(Liuetal., 1983; Romer and Kroner, 2016; Yuanetal., 2019; Zhaoetal., 2022a),为成矿提供了钨锡的预富集;(2)受古太平洋板块俯冲影响(Maoetal., 2021),钨锡花岗岩源区部分熔融温度条件适宜,促进钨锡从岩浆源区有效释放到初始岩浆中(Yuanetal., 2019; Songetal., 2022; Zhaoetal., 2022a, d);(3)成矿花岗岩具有高分异、还原的特征,有利于钨锡在岩浆演化过程中进一步富集(Zhaoetal., 2016, 2022c; Songetal., 2018; Yuanetal., 2018; 毛景文等, 2023; 吴福元等, 2023)。上述研究表明南岭地区钨锡成矿条件优渥,是开展进一步钨锡找矿勘查的理想地区。然而,如何开展钨锡找矿工作,尤其是在已有矿床周围高效圈定钨锡成矿远景区,是钨锡成矿作用研究以及找矿勘查部署工作面临的重要科学问题。
湘南地区的柿竹园和荷花坪矿田是南岭成矿带内两个具有代表性的大型-超大型锡钨多金属矿田,前人对其成岩成矿年龄、成矿机制以及成矿构造背景开展了一系列的研究(刘义茂等, 1995; 赵振华等, 2000; 龚庆杰等, 2004; 蔡明海等, 2006; 章荣清等, 2010; 刘晓菲等, 2012; 袁顺达等, 2012a, b; Maoetal., 2013; Zhaoetal., 2018),并取得了丰富的成果。然而,尽管两个矿田空间相邻(约5km),成矿时代相近,并且中间发育一系列钨锡矿床和矿化点,但是由于以往的研究认为荷花坪锡矿形成于晚三叠世(蔡明海等, 2006),与柿竹园地区晚侏罗世的大规模成矿作用存在明显时差,因而长期以来鲜有人关注两个矿田之间的成因联系以及两矿田之间的区域是否具有较大钨锡找矿潜力。
湘南双园冲锡矿位于荷花坪与柿竹园锡钨矿田中间部位,与两大型矿田处于同一北东向构造体系。目前对该矿床的研究程度较低,尤其是成岩成矿时代鲜有人涉及,制约了对该矿床成因的研究以及该区的找矿勘查部署工作。本文对双园冲锡矿和成矿花岗岩开展了系统的锆石、独居石和锡石U-Pb年代学研究,并在厘定双园冲锡矿成岩成矿时代的基础上结合矿床地质特征及成矿花岗岩岩石学特征的分析,提出荷花坪锡多金属矿田与柿竹园锡多金属矿田具有密切成因联系,两者中间具有较大锡钨找矿潜力。
华南陆块是扬子陆块和华夏陆块于新元古代沿江南造山带拼合而成,南岭地区位于华南腹地,处于两个陆块缝合带中段地区。此后,南岭地区构造格局先后受加里东造山运动、印支造山运动以及古太平洋板块俯冲-弧后伸展的影响,发育了大量花岗岩(舒良树, 2012; Yuanetal., 2015; Zhaoetal., 2019; 赵正等, 2022)。
双园冲矿床地处湖南省南部,位于南岭成矿带中西段。矿区出露地层主要有中泥盆统跳马涧组、棋梓桥组,岩性以砂岩、灰岩及白云质灰岩为主。区内断裂构造复杂,主要发育有NE向、NNE向、NW向、近SN向、EW向五组断裂,其中NE-NNE向断裂与成矿关系密切,为导矿、容矿构造(钟江临和李楚平, 2013),控制脉状锡矿的产出(图2);区内向斜呈NE走向,核部地层为棋梓桥组,两翼地层为跳马涧组。矿区出露的岩浆岩主要为花岗斑岩脉,呈NE向展布。深部隐伏岩体为中细粒云英岩化花岗岩,岩石整体呈灰绿色,变余花岗结构,块状构造,主要由长石、石英、白云母、绢云母及少量金属矿物和副矿物组成,其中长石蚀变强烈,生成大量细小的绢云母和粘土矿物(图3a、b)。
图1 南岭成矿带柿竹园-荷花坪地区地质简图(据章荣清等, 2016修改)Fig.1 Sketch geological map of the Shizhuyuan-Hehuaping area in the Nanling W-Sn metallogenic belt (modified after Zhang et al., 2016)
图2 双园冲矿床63线剖面(据钟江临和李楚平,2013修改)Fig.2 The Line 63 profile of the Shuangyuanchong deposit (modified after Zhong and Li, 2013)
图3 双园冲锡矿区云英岩化花岗岩(a、b)、云英岩型矿石(c、d)、石榴子石矽卡岩(e、f)和绿泥石硫化物矿石(g、h)手标本(左栏)和显微(右栏)照片图(b)中长石蚀变强烈,生成大量细小的绢云母和粘土矿物. (b、d)为正交偏光下照片,(f)为单偏光下照片,(h)为反射光下照片. Fsp-长石; Grt-石榴子石; Ms-白云母; Py-黄铁矿; Pyh-磁黄铁矿; Qz-石英; Ser-绢云母; Sp-闪锌矿Fig.3 Representative hand specimen photographs (left column) and microphotographs (right column) of greisenized granite (a, b), greisen-type Sn ore (c, d), garnet skarn (e, f) and chlorite sulfide ore (g, h) from the Shuangyuanchong Sn depositIn Fig.3b, feldspar alteration is strong, resulting in a large number of fine sericite and clay minerals. (b, d) are cross-polarized microphotographs, (f) is single-polarized microphotograph, and (h) is reflected light microphotograph. Fsp-feldspar; Grt-garnet; Ms-muscovite; Py-pyrite; Pyh-pyrrhotite; Qz-quartz; Ser-sericite; Sp-sphalerite
双园冲矿床矿体主要受NE-NNE向断裂带控制,倾向SE,赋存于泥盆系棋梓桥组下段灰岩与跳马涧组砂岩过渡层以及棋梓桥组下段下部层位的层间破碎带中,呈脉状、透镜状、似层状和层状产出。矿体品位Sn介于0.2%~1.52%,Pb介于0.7%~16.24%,Zn介于0.02%~6.06% (钟江临和李楚平, 2013)。矿化类型主要为云英岩型和矽卡岩型。据钻孔剖面显示,由深部至浅部,依次为云英岩型、矽卡岩型脉状Sn矿、矽卡岩型脉状Pb-Zn矿、Fe-Mn矿,呈现清晰的矿化蚀变分带(图2),其中Fe-Mn矿露头较差,风化较为严重。云英岩型矿石呈浅绿色-灰黑色,细粒变晶结构,块状构造,主要由半自形粒状石英(75%~80%),片状白云母(15%~20%),柱状长石(3%~5%)及锡石和黄铁矿等矿物组成(图3c, d)。矽卡岩型矿石主要包括进变质石榴子石矽卡岩和退变质的绿泥石硫化物矿石。前者呈灰绿色,细粒变晶结构,块状构造,主要由石榴子石、透辉石、黄铁矿等矿物组成。显微镜下石榴子石呈自形-半自形粒状,单偏光条件下呈浅粉色至黄褐色,糙面明显,具不规则裂纹,边缘和裂隙强烈蚀变为绿泥石和粘土矿物(图3e, f);绿泥石硫化物矿石呈灰黑色,块状构造,主要由黄铁矿、磁黄铁矿、闪锌矿和方铅矿等金属矿物和少量绿泥石组成(图3g, h)。
本次研究的花岗岩样品和矿石样品均采于双园冲矿区井下。锆石和独居石样品来自云英岩化花岗岩,锡石样品来自云英岩型矿石。首先将样品粉碎至80~100目后用常规浮选和电磁选的方法进行分选,然后在双目镜下挑选出晶形较好的锡石、独居石和锆石颗粒。
将无夹杂、无裂纹、无明显蚀变的矿物颗粒与外标在双目显微镜下用双面胶布粘在载玻片上,并在载玻片上放置PVC环,然后将环氧树脂和固化剂的混合物放入PVC环中,直至其完全凝固。最后,将样品靶取下来打磨抛光直至露出大部分矿物颗粒的内部部分。上机分析前,用乙醇清洗样品表面,以消除表面污染。
2.2.1 LA-ICP-MS锆石和独居石U-Pb定年
锆石和独居石LA-ICP-MS U-Pb同位素定年在北京中科矿研检测技术有限公司完成。激光剥蚀系统为ESI公司生产的NWR 193nm激光器,ICP-MS型号为Agilent 7500。激光剥蚀过程中采用氦气作载气,使用氩气为补偿气以调节灵敏度,二种气体在进入ICP之前通过一个Y型接头混合。单个样品测试过程包括大约15~20s的空白信号和45s的样品信号。测试获得的数据采用ICPMSDataCal软件进行离线处理,包括对样品和空白信号的选择、仪器灵敏度漂移校正、元素含量以及U-Th-Pb同位素比值和年龄计算,详细流程参考Liuetal. (2010)和侯可军等(2009)。
2.2.2 LA-ICP-MS锡石U-Pb定年
锡石测定采用非基体匹配校正方法,使用锆石标准物GJ-1作为外标进行同位素分馏校正,并采用Tera-Wasserburg图解法和207Pb法对普通铅进行校正(陈靖等,2021)。分析时采用的激光束斑直径为50μm,剥蚀频率为6Hz,能量密度为5J/cm2。分析后将实测数据结果进行线性拟合,得到上交点即为样品的普通Pb组成,下交点近似代表样品的形成年龄(Tera and Wasserburg, 1972)。然后根据上交点所获得的初始207Pb/206Pb同位素组成,对含有较高普通铅的样品进行207Pb校正,获得校正后的206Pb/238U年龄。锡石样品的T-W图绘制和加权平均年龄计算均采用Isoplot完成(Ludwig, 2001)。
双园冲花岗岩锆石测年结果见表1。锆石阴极发光图像(图4a)显示,双园冲花岗岩中的锆石多为半自形-自形柱状,长度20~120μm,长宽比为1.5~2,颜色为黑色-灰白色,发育有密集的振荡环带,大部分锆石由于颜色较深,环带不清晰。测试过程挑选晶形较好、颜色较浅、环带清晰、无包裹体和裂隙的锆石开展测试。测试结果显示,锆石206Pb/238U值变化范围为0.0258~0.0419,207Pb/235U值变化范围为0.2128~1.6537,207Pb/206Pb值变化范围为0.0602~0.2898。锆石207Pb/206Pb-238U/206Pb反等时线年龄为161.2±2.5Ma (MSWD=1.8),加权平均年龄为161.0±2.9Ma (MSWD=1.2)(图5a)。
图4 双园冲矿区云英岩化花岗岩中代表性锆石(a)、独居石(c)以及云英岩型矿石中代表性锡石(b)的CL图像Fig.4 Representative cathodoluminescence (CL) of zircons (a), monazite (c) for the greisenized granite and cassiterite (b) for the greisen-type ore from the Shuangyuanchong Sn deposit
表1 双园冲锡矿云英岩化花岗岩锆石LA-ICP-MS U-Pb定年分析结果
独居石的测年结果见表2。本次所采花岗岩样品中的独居石主要呈自形-半自形,长度为50~110μm,长宽比为1~3,整体颜色为均一的灰白色,独居石表面局部发育有裂隙和孔洞,边部偶有熔蚀现象(图4c)。对环带规则、无包裹体和裂隙的独居石测年结果显示,独居石Pb、Th、U总量分别为3.8×10-4~8.8×10-4、5.5×10-2~1.3×10-1、1.1×10-3~3.3×10-3,获得207Pb/206Pb-238U/206Pb的反等时线年龄为157.1±1.8Ma (MSWD=0.61),加权平均年龄为157.1±1.1Ma (MSWD=0.26)(图5b)。
图5 双园冲矿区云英岩化花岗岩中锆石(a)、独居石(b)以及云英岩型矿石中锡石(c)LA-ICP-MS U-Pb年龄结果Fig.5 LA-ICP-MS U-Pb ages of zircon (a), monazite (b) for the greisenized granite and cassiterite (c) for the greisen-type ore from the Shuangyuanchong Sn deposit
表2 双园冲锡矿云英岩化花岗岩独居石LA-ICP-MS U-Pb定年分析结果
续表2
锡石测年结果见表3。云英岩型矿石样品中锡石为深灰色-黑色,半自形-自形粒状结构,长度为50~190μm,长宽比为1~2.1,生长环带发育,表面无裂隙(图4b)。锡石测试结果显示,206Pb/238U值变化范围为0.0243~0.0799,207Pb/206Pb值变化范围为0.0513~0.5824,得到207Pb/206Pb-238U/206Pb反等时线年龄为158.9±2.9Ma (MSWD=3.9),207Pb法普通铅校正后的206Pb/238U加权平均年龄为158.6±2.6Ma (MSWD=3.3)(图5c)。
表3 双园冲锡矿锡石LA-ICP-MS U-Pb定年分析结果
锡石是锡矿床主要的矿石矿物,对其U-Pb同位素测年是确定锡矿床成矿时代的有效手段,已经被应用于各种含锡石热液矿床的年代学研究中(Yuanetal., 2008, 2011; Zhangetal., 2015; Dengetal., 2018; Liuetal., 2018)。然而,由于基体效应的影响和标样的缺乏,在一定程度上制约该测年方法的广泛应用。为了进一步拓展锡石U-Pb测年方法的应用,在使用He做载气以及加入不同辅助气(N2、水蒸气)条件下,利用锆石作为外部标准物质进行同位素分馏校准,通过采用大的激光剥蚀斑束(50μm)及浅的剥蚀深度克服基体效应,并采用Tera-Wasserburg图解法和207Pb法对普通铅进行校正,获得了准确的锡石206Pb/238U年龄,显示采用锆石作为标准物质对锡石样品进行原位LA-ICP-MS U-Pb定年是可行的(陈靖等, 2021)。本次研究采用陈靖等(2021)的方法以锆石为标样获得双园冲矿床锡石Tera-Wasserburg U-Pb下交点年龄为158.9±2.9Ma,与207Pb法普通铅校正后的206Pb/238U加权平均年龄(158.6±2.6Ma)一致,指示锡矿化年龄为159Ma左右,该年龄与本次研究获得的与锡矿密切共生花岗岩的独居石U-Pb年龄(157.1±1.8Ma)一致,表明花岗岩侵位和锡矿化同时形成于晚侏罗世。此外,本次研究获得的花岗岩锆石U-Pb年龄为161.2±2.5Ma,尽管由于锆石U含量过高导致测年结果误差较大(李秋立, 2016),但仍与上述成岩成矿年龄在误差范围内基本一致。结合双园冲花岗岩与锡矿在空间上的密切接触关系以及矿化蚀变分带特征(图2),表明花岗岩与锡矿化具有密切的成因联系。
从整个区域上来看,柿竹园-荷花坪地区分布有一系列大型-超大型锡钨多金属矿床,是南岭钨锡成矿带锡多金属矿床集中产出地区。前人对该区锡钨矿形成时代开展了大量的研究,例如Zhaoetal. (2018)对柿竹园地区锡钨多金属矿床的研究显示柿竹园W-Sn-Mo-Bi矿床形成于151±1Ma,野鸡尾Sn-Cu矿床形成于152.7±0.9Ma,玛瑙山W-Sn-Fe-Mn-Pb-Zn矿形成于154~158Ma,枞树坂Pb-Zn矿形成于152±2Ma,上述年龄与前人获得的该矿区红旗岭Sn-W-Pb-Zn矿年龄(153.3±1.0Ma, 袁顺达等, 2012a)、金船塘Sn-Bi矿年龄(158.8±6.6Ma, 刘晓菲等, 2012)一致,并且与荷花坪地区荷花坪锡矿(154~156Ma, Zhangetal., 2015; Zhaoetal., 2022b)、打鼓岭锡矿(157.8±4.1Ma, 章荣清, 2014)一致,指示双园冲锡矿与柿竹园-荷花坪地区的Sn-W-Mo-Bi-Pb-Zn多金属矿床集中形成于160~150Ma。此外,南岭成矿带西段香花岭Sn矿(154~158Ma, Yuanetal., 2007, 2008; Wuetal., 2022)、黄沙坪Pb-Zn-Sn-W多金属矿(155~160Ma, 雷泽恒等, 2010a; 原垭斌等, 2014; 赵盼捞等, 2018)、芙蓉Sn矿(156~160Ma, 毛景文等, 2004; Yuanetal., 2011)、新田岭W矿(159~162Ma,袁顺达等, 2012b)、魏家W矿(151~162Ma, Zhaoetal., 2016; Huangetal., 2017)等锡(钨)多金属矿床均形成于160~150Ma,该时期被称为南岭地区锡钨爆发式成矿期(毛景文等, 2007; Maoetal., 2013, 2019; 袁顺达等, 2020)。Maoetal. (2021)进一步通过对比南岭地区锡钨矿和东南沿海斑岩-矽卡岩型铜矿的时空分布特征提出南岭地区晚侏罗世锡钨矿床形成于古太平洋板块向华南板块俯冲形成的弧后伸展构造背景。
荷花坪与柿竹园钨锡多金属矿田是南岭地区空间相邻的两个大型-超大型矿田,近年来找矿成果显著,尤其是在荷花坪东南缘与千里山的西南缘发现了一系列钨锡矿床、矿化点,如柿竹园矿田西南缘的玛瑙山、横山岭和蛇形坪Pb-Zn矿深部发现了一系列W-Sn矿体等(蔡新华等, 2006; 雷泽恒等, 2010b; Zhaoetal., 2018)。然而,由于荷花坪Sn矿长期一直被认为形成于晚三叠世(蔡明海等, 2006),柿竹园与荷花坪Sn多金属矿田一直被视作两个独立的矿田开展研究。而对于玛瑙山矿床,尽管其与柿竹园Sn-W多金属矿床空间距离较远,以往一直被认为是与柿竹园Sn-W多金属矿有关的远端脉状铅锌矿。然而,千里山岩体为典型的富氟花岗岩,富氟岩浆的固相线可低至525℃(Websteretal., 1987; 熊小林等, 1998),此时远端Pb-Zn矿可能在离岩体较近的部位就位,而玛瑙山Pb-Zn-Mn矿床离柿竹园钨锡矿超过4km,因而前期的研究认为二者来源于独立演化的热中心,玛瑙山深部具有进一步的钨锡找矿潜力(袁顺达等, 2012a; 袁顺达, 2017),这一认识已经得到了最近玛瑙山Pb-Zn-Mn矿区深部钨锡找矿突破工作的证实。这样看来,双园冲锡矿与玛瑙山Sn(W)-Pb-Zn-Mn矿空间距离更为接近,二者正好位于荷花坪和柿竹园矿田中间(图1),说明柿竹园矿田与荷花坪矿田之间是相连一体的。为了进一步明确双园冲锡矿与柿竹园和荷花坪两个大型锡多金属矿田之间的联系,本文就三个矿床的成岩成矿特征进行了系统地对比,并对该区的成矿潜力进行了初步的研究。
首先,从成岩成矿时代来看,双园冲矿床与荷花坪锡多金属矿田(154~161Ma, Zhangetal., 2015)和柿竹园锡多金属矿田(151~155Ma, Lietal., 2004; Zhaoetal., 2018)形成时代基本一致。从矿体产出位置来看,锡矿体均主要产于花岗岩与中-上泥盆统棋梓桥组、佘田桥组或中泥盆统跳马涧组灰岩中间部位,或者产于中泥盆统跳马涧组中上部砂岩与中泥盆统棋梓桥组灰岩的过渡带(钟江临和李楚平, 2013; 王勇, 2021)。其次,从构造方面,成矿岩体和矿体受NE向断裂的控制,均呈北东向展布,北东向断裂是成矿岩浆和热液活动的重要通道。从成矿岩体岩石特征上来看,三个矿区成矿岩体均为黑云母花岗岩,且位于双园冲矿床两侧的柿竹园和荷花坪矿区成矿岩体具有一致的全岩Nd和锆石Hf同位素特征(表4、图6),显示相同的岩浆源区组成,表明三者的成矿花岗岩可能来自同一深部的岩浆房。此外,宋才见等(2020)通过重力等地球物理数据也提出柿竹园和荷花坪矿区深部花岗岩可能是相连的。综合以上信息,本文提出如下柿竹园-荷花坪地区钨锡矿床成矿模型(图7):该地区深部存在一个大的岩浆房,岩浆结晶分异形成富含钨锡等成矿金属和挥发分(氟、氯、水)的岩浆,由于这种高分异的岩浆密度较小,在岩浆房顶部突起部位聚集上侵,尤其是沿着NE向断裂等地壳薄弱部位上升侵位形成岩枝, 岩枝出溶富含成矿金属的流体在岩体顶部由于自交代作用形成云英岩型蚀变和相应的钨锡矿床,在碳酸盐岩围岩接触带部位形成矽卡岩型蚀变和相应钨锡矿床,在距离岩体更远的浅部形成脉状Pb-Zn-Fe-Mn矿化。因此,柿竹园-荷花坪钨锡矿田中间区域具有较大的钨锡找矿潜力,尤其是隐伏的晚侏罗世高分异花岗岩体与泥盆系碳酸盐岩围岩的接触带部位。
表4 荷花坪锡矿与柿竹园钨锡矿成矿岩体Nd同位素和锆石Hf同位素组成
图6 荷花坪锡矿与柿竹园钨锡矿成矿岩体锆石εHf(t)值(a)和全岩εNd(t)值(b)频数直方图Fig.6 The frequency histograms of zircon εHf(t) (a) and whole-rock εNd(t) (b) values for ore-related granite from the Hehuaping Sn deposit and Shizhuyuan W-Sn deposit
图7 南岭地区柿竹园-荷花坪地区钨锡矿床成矿模型Fig.7 The genesis model of W-Sn deposit in the Shizhuyuan-Hehuaping area, Nanling region
由于矽卡岩-云英岩型锡钨矿体体积和蚀变范围较小,并且多为隐伏矿体,如何有效的开展锡钨找矿勘查是广受关注但未较好解决的一个难题。袁顺达等(2012a),袁顺达(2017)和Zhaoetal. (2018)通过研究湘南地区钨锡多金属矿床特征发现,深部隐伏的云英岩-矽卡岩型W-Sn矿床在浅部往往发育同时代的脉状铅锌矿化和少量的锰铁矿化,两种矿化为同一套岩浆热液体系演化的结果,并据此进一步提出湘南地区浅部脉状铅锌矿和锰铁矿化深部具有较大的寻找同时期的钨锡矿的潜力。双园冲锡矿浅部发育较多的脉状铅锌矿以及铁锰矿化(图2),与上述找矿模型一致,并且近期在千里山岩体周缘的玛瑙山锰铅锌矿、妹子垄铅锌矿、横山岭-蛇形坪铅锌矿等深部发现一系列大中型W-Sn矿体及矿化点(蔡新华等, 2006; 雷泽恒等, 2010b; 郭海等, 2011; Zhaoetal., 2018),进一步证实了此找矿模型的有效性。因此南岭地区铅锌矿体深部,尤其是目前距离已知钨锡多金属矿床较远的铅锌矿体深部是钨锡多金属矿床找矿勘查的有利部位。
(1)南岭成矿带双园冲锡矿的锡矿化年龄(158.9±2.9Ma)与其成锡花岗岩的侵位年龄(157.1±1.8Ma)基本一致,均形成于晚侏罗世,结合矿床地质特征,指示两者具有密切成因联系。
(2)对比整个区域,双园冲锡矿与荷花坪-柿竹园地区Sn-W多金属矿床及南岭成矿带西段的锡多金属矿床集中形成于160~150Ma,形成于古太平洋板块向华南板块俯冲形成的弧后伸展构造背景。
(3)双园冲、荷花坪和柿竹园三个锡多金属矿床可能形成于同一次岩浆热液活动,其成矿岩体可能来自同一个大岩浆房,因此柿竹园和荷花坪之间的区域具有较大发育晚侏罗世花岗岩体及相关钨锡矿的潜力,研究区乃至整个南岭地区浅部脉状Pb-Zn-Fe-Mn矿化是深部Sn-W找矿勘查的有利部位。
致谢感谢北京中科矿研检测技术有限公司刘力老师对本文LA-ICP-MS U-Pb定年测试的支持、中国地质大学(北京)宋世伟副教授对样品薄片显微观察的指导;感谢审稿专家和编辑部老师提出的建议性意见和建议。