岩石学报  2021, Vol. 37 Issue (3): 637-664, doi: 10.18654/1000-0569/2021.03.01   PDF    
内蒙古维拉斯托稀有金属-锡多金属矿床的成矿作用: 来自花岗质岩浆结晶分异的启示
武广1, 刘瑞麟1,2, 陈公正1, 李铁刚1, 李睿华1,2, 李英雷1,2, 杨飞1,2, 张彤3     
1. 中国地质科学院矿产资源研究所, 自然资源部成矿作用与资源评价重点实验室, 北京 100037;
2. 北京大学, 造山带与地壳演化教育部重点实验室, 北京 100871;
3. 内蒙古自治区地质调查院, 呼和浩特 010020
摘要: 内蒙古维拉斯托稀有金属-锡多金属矿床位于大兴安岭南段西坡,是一个以锡为主,共伴生锌、钨、铜、钼、铌、钽、锂和铷的大型矿床。该矿床以发育铌、钽、铷、锂矿化有别于大兴安岭南段的其他锡多金属矿床。本文对维拉斯托地区的北大山岩体和维拉斯托岩体开展了LA-ICP-MS锆石U-Pb定年、全岩地球化学研究和锆石Lu-Hf同位素分析。北大山岩体和维拉斯托岩体的结晶年龄分别为140±2Ma和137±1Ma~138±1Ma。2个岩体均具有高硅、富碱、贫铁镁钙特征,主体为高钾钙碱性的弱过铝质碱长花岗岩,轻重稀土分馏不明显,呈现出显著的负铕异常。2个岩体均富集Rb、Th、U、K、Ta、Hf和轻稀土元素,强烈亏损Ba、Sr、P、Eu和Ti等元素。与北大山岩体相比,维拉斯托岩体有更低的稀土总量,更明显的Eu负异常,展示了明显的稀土元素四分组效应;Ba、Sr、Eu、Ti和Y亏损更强烈。二者的时空关系和地球化学特征表明2个岩体为同一期岩浆活动不同演化阶段的产物,岩浆演化过程经历了钾长石、斜长石、黑云母、独居石、榍石、褐帘石、磷灰石、磷钇矿、锆石和钛铁矿的分离结晶,维拉斯托岩体还发生了熔体-流体相互作用。北大山岩体和维拉斯托岩体锆石的εHft)值分别为5.4~8.6和4.2~9.9,两阶段Hf同位素模式年龄分别为643~847Ma和556~921Ma,表明它们的岩浆源区为起源于亏损地幔的新元古代新生地壳物质。北大山岩体属于高分异I型花岗岩,维拉斯托岩体发育钠长石、天河石、黄玉、锂云母和萤石,属于超分异的Li-F花岗岩。花岗质岩浆的高程度分离结晶和熔体-流体相互作用是形成维拉斯托超分异花岗岩并发生稀有金属-锡多金属成矿作用的重要控制因素。维拉斯托稀有金属-锡多金属矿床属于岩浆-热液矿床。维拉斯托地区的早白垩世花岗质岩石形成于伸展环境。
关键词: 锆石U-Pb定年    地球化学    锆石Lu-Hf同位素    高分异I型花岗岩    北大山岩体和维拉斯托岩体    维拉斯托稀有金属-锡多金属矿床    大兴安岭南段    
Mineralization of the Weilasituo rare metal-tin-polymetallic ore deposit in Inner Mongolia: Insights from fractional crystallization of granitic magmas
WU Guang1, LIU RuiLin1,2, CHEN GongZheng1, LI TieGang1, LI RuiHua1,2, LI YingLei1,2, YANG Fei1,2, ZHANG Tong3     
1. MNR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 100037, China;
2. MOE Key Laboratory of Orogenic Belts and Crustal Evolution, Peking University, Beijing 100871, China;
3. Inner Mongolia Institute of Geological Survey, Hohhot 010020, China
Abstract: The Weilasituo rare metal-tin-polymetallic deposit of Inner Mongolia, located on the west slope of the southern Great Xing'an Range, NE China, is a large-sized Sn deposit associated with Zn, W, Cu, Mo, Nb, Ta, Li, and Rb. The deposit is different from other Sn-polymetallic deposits in the southern Great Xing'an Range by developing Nb, Ta, Rb, and Li mineralization. In this paper, LA-ICP-MS zircon U-Pb dating, whole-rock geochemistry, and zircon Lu-Hf isotopic analysis for the Beidashan and Weilasituo intrusions in the Weilasituo area were carried out. The crystallization ages of the Beidashan and Weilasituo intrusions are 140±2Ma and 137±1Ma~138±1Ma, respectively. Both of them have high SiO2 and total alkali contents and low Fe2O3, FeO, MgO, and CaO contents, mainly belonging to weakly peraluminous alkali-feldspar granite with a high-K calc-alkaline affinity. The two intrusions are characterized by enrichment of Rb, Th, U, K, Ta, Hf, and LREE and strong depletion of Ba, Sr, P, Eu, and Ti, with extremely negative Eu anomalies and weak fractionation between LREE and HREE. Compared with the Beidashan pluton, the Weilasituo stock has lower total REE content, more obvious negative Eu anomalies, and more strong depletion of Ba, Sr, Eu, Ti, and Y, showing an obvious tetrad effect of rare earth elements. The temporal-spatial relationship and geochemical characteristics of the two intrusions indicate that they were the products of different evolution stages of the same magmatic activity. The magma experienced the fractional crystallization of K-feldspar, plagioclase, biotite, monazite, sphene, allanite, apatite, xenotime, zircon, and ilmenite during its evolution. In addition, the interaction between melt and fluid occurred during the formation of the Weilasituo stock. The zircon εHf(t) values of the Beidashan pluton and the Weilasituo stock are 5.4~8.6 and 4.2~9.9, respectively, and the two-stage Hf isotopic model ages are 643~847Ma for the Beidashan pluton and 556~921Ma for the Weilasituo stock, indicating that their magma sources were Neoproterozoic juvenile crustal materials originated from depleted mantle. The Beidashan pluton is highly fractionated Ⅰ-type granite, and the Weilasituo stock belongs to ultra fractionated Li-F granite, which is characterized by albite, amazonite, topaz, lepidolite, and fluorite typomorphic minerals. The high degree fractional crystallization of granitic magma as well as fluid-melt interaction is the important controlling factors for the formation of the Weilasituo ultra fractionated granite and the rare metal-Sn-polymetallic deposit. The Weilasituo rare metal-Sn-polymetallic deposit genetically belongs to magmatic-hydrothermal deposit. The Early Cretaceous granitic rocks in the Weilasituo area formed in an extensional environment.
Key words: Zircon U-Pb dating    Geochemistry    Zircon Lu-Hf isotope    Highly fractionated Ⅰ-type granite    Beidashan pluton and Weilasituo stock    Weilasituo rare metal-tin-polymetallic deposit    Southern Great Xing'an Range    

大兴安岭南段处于古亚洲成矿域和滨西太平洋成矿域的交汇部位,是中国北方最重要的锡银多金属成矿带(Wang et al., 2001; 毛景文等, 2013; Ouyang et al., 2015)。这些锡多金属矿床主要表现为锡银铅锌(铜钨钼)组合(Liu et al., 2001; Zhou et al., 2012; Zhai et al., 2014, 2017; Zeng et al., 2015; Wang et al., 2016a; Zhang et al., 2017, 2019; 陈公正, 2018; 李睿华, 2019; Mi et al., 2020)。维拉斯托稀有金属-锡多金属矿床为Sn-Zn-W-Mo-Li组合(刘瑞麟等, 2018a, b),以深部的岩体顶部花岗岩中的Sn-Zn(-Rb-Nb-Ta)矿化和浅部隐爆角砾岩筒内的Sn-Zn-Li(-Rb-Be)矿化区别于该区的其他锡多金属矿床。前人对维拉斯托稀有金属-锡多金属矿床的研究主要集中在成岩成矿时代的厘定(郭贵娟, 2016; Liu et al., 2016; 翟德高等, 2016; 祝新友等, 2016; Wang et al., 2017; 刘瑞麟等, 2018a)和流体包裹体研究(刘瑞麟等, 2018b)。维拉斯托稀有金属-锡多金属矿床研究存在的主要问题是:缺少成岩成矿关系的系统研究,基本上未涉及岩浆演化对矿床的矿化类型和矿种的制约,导致矿床成因类型一直存在争议。多数研究者认为该矿床属于斑岩成矿系统(刘翼飞等, 2014; Liu et al., 2016; 翟德高等, 2016; 樊志勇等, 2017; Wang et al., 2017),但亦有部分研究者认为属于与花岗岩有关的岩浆-热液成矿系统(祝新友等, 2016; 刘瑞麟等, 2018b);此外,该区不同的锡多金属矿床在矿种组合上存在较大区别,其原因一直未得到应有的重视。鉴于上述原因,本次选取维拉斯托稀有金属-锡多金属矿床为研究对象,在详细的野外地质调查基础上,对矿区外围的北大山岩体和矿区深部隐伏的维拉斯托岩体开展系统的岩相学、LA-ICP-MS锆石U-Pb定年、全岩地球化学和锆石Lu-Hf同位素分析,目的是厘定维拉斯托岩体与北大山岩体的成因联系,查明岩体的源区组成和岩石成因,揭示岩浆演化与成矿关系,厘定矿床成因类型,建立维拉斯托稀有金属-锡多金属矿床成矿模式,结合前人对大兴安岭南段典型锡多金属矿床及其成矿岩体的研究成果,探讨岩浆源区组成、岩浆演化过程对锡多金属成矿的制约,为大兴安岭南段锡银多金属找矿提供理论支撑,为该区成矿规律研究提供依据。

1 区域地质背景

大兴安岭南段位于中亚造山带东段(图 1a)。二连-贺根山断裂带为其北界与兴安地块相分隔;西拉木伦断裂带为南界与华北板块北缘的西拉木伦造山带相隔;东以嫩江-八里罕断裂为界,与松辽盆地分隔;西界大致位于锡林浩特-东乌珠穆沁旗一线(图 1b)。前中生代该区受古亚洲洋构造体制控制,中生代期间主要受蒙古-鄂霍茨克洋和古太平洋两大构造体制控制(Wang et al., 2012; Xu et al., 2013)。大兴安岭南段出露的地层有中元古界、下古生界、上古生界、侏罗系-白垩系和新生界(图 1b)。中元古代锡林郭勒杂岩(本文称宝音图群)主要岩石类型为斜长角闪片麻岩和黑云斜长片麻岩。长期以来锡林郭勒杂岩的时代归属一直存在较大争议,内蒙古自治区地质矿产局(1991)认为其形成于古元古代;一些研究者对锡林郭勒杂岩中片麻岩测年获得了早古生代的年龄,认为锡林郭勒杂岩并非古老陆块(薛怀民等, 2009);但更多的年龄数据表明原锡林郭勒杂岩既包括中元古代形成的古老基底,又包括新元古代基性-超基性岩和古生代花岗岩类,解体后的中元古代基底被称为锡林浩特岩群(葛梦春等, 2011; 于洋等, 2012)。下古生界为海相火山岩-碎屑岩,上古生界为海相火山岩-碎屑岩和陆相碎屑岩,侏罗系-白垩系为陆相火山岩-碎屑岩,新生代包括玄武岩以及松散沉积物(内蒙古自治区地质矿产局, 1991)。晚古生代期间,伴随着古亚洲洋的闭合,形成了NEE-NE向褶皱和断裂构造,奠定了本区的基本构造格架;中生代则以NE-NNE向断裂为主。区内主要发育燕山期、印支期和华力西期3期侵入岩,其中燕山期侵入岩规模最为宏大(图 1b)。

图 1 大兴安岭南段大地构造位置图(a, 据Ge et al., 2007修改)及主要的锡多金属矿床分布图(b, 地质底图据徐志刚等, 2008修编) 锡多金属矿床编号:1-白音查干东山;2-毛登-小孤山;3-花敖包特;4-查木罕;5-维拉斯托;6-安乐;7-黄岗;8-老盘道背后;9-道伦达坝;10-东山湾;11-宝盖沟;12-沙布楞山;13-莫古吐;14-边家大院;15-大井;16-白音诺尔;17-浩布高;18-敖脑达坝;19-孟恩陶勒盖 Fig. 1 Sketch geological map of the southern Great Xing'an Range and its neighboring areas, showing geotectonic units (a, modified after Ge et al., 2007) and locations of major Sn-polymetallic and Sn-bearing polymetallic deposits (b, geological base map modified after Xu et al., 2008) Names of numbered Sn-(Ag)-polymetallic deposits: 1-Baiyinchagandongshan; 2-Maodeng-Xiaogushan; 3-Huaaobaote; 4-Chamuhan; 5-Weilasituo; 6-Anle; 7-Huanggang; 8-Laopandaobeihou; 9-Daolundaba; 10-Dongshanwan; 11-Baogaigou; 12-Shabulengshan; 13-Mogutu; 14-Bianjiadayuan; 15-Dajing; 16-Baiyinnuoer; 17-Haobugao; 18-Aonaodaba; 19-Mengentaolegai

北大山岩体和维拉斯托岩体所在的维拉斯托地区位于大兴安岭南段主脊断裂与西坡断裂之间的米生庙复背斜东南翼,北距西乌珠穆沁旗政府所在地约60km(图 1b)。区域上出露地层主要有中元古界宝音图群、上石炭统本巴图组和阿木山组、下-中二叠统大石寨组、上二叠统林西组、中侏罗统万宝组、上侏罗统满克头鄂博组和第四系(图 2)。宝音图群岩性为斜长角闪片麻岩和黑云斜长片麻岩;本巴图组以海相粉砂岩和杂砂岩为主,夹安山质凝灰岩、安山岩及灰岩;阿木山组为一套海相碎屑岩、碳酸盐岩建造;大石寨组为一套海相碎屑岩与火山碎屑岩;林西组为一套陆相碎屑岩建造;万宝组主要是陆相含碳泥岩、粉砂岩及砾岩夹煤层;满克头鄂博组为一套以流纹岩、流纹质熔结凝灰岩为主的陆相酸性火山岩;第四系冲积层及风成砂土广泛分布(刘瑞麟等, 2018b)。区域内主要的构造由贯穿全区的米生庙复背斜和断裂组成。复背斜和断裂以NEE-NE向为主,NW向断裂及近EW向断裂也较发育,构成了该区近网格状的构造格局。NEE-NE向的米生庙复背斜宽约60km,由中元古界宝音图群组成复背斜轴部,石炭系、二叠系组成翼部(图 2)。NE向断裂控制着区内岩体及脉岩的分布,其中NE向的华力西期断裂控制着区内石英闪长岩及大多数脉岩的分布;NE向的燕山期断裂控制着区内花岗岩的分布。NE向的燕山期断裂为主要的导矿和容矿构造,而EW向的燕山期断裂仅为容矿构造。区域侵入岩岩浆活动主要发生在华力西期和燕山期,受NE向断裂控制,呈岩基或岩株状分布。华力西期花岗岩类以石英闪长岩为主,亦有花岗岩和花岗闪长岩,该期岩体侵入中元古界宝音图群及上石炭统本巴图组,并在下-中二叠统砂砾岩内见其角砾。王瑾(2009)王新宇等(2013)对维拉斯托地区的花岗岩类进行了SHRIMP锆石U-Pb测年,获得了320.5±4.1Ma的黑云母花岗岩的结晶年龄和313.9±3.4Ma黑云母二长花岗岩的结晶年龄(王瑾, 2009)、298.0±2.5Ma的花岗闪长岩结晶年龄和308.3±4.2Ma的石英闪长岩结晶年龄(王新宇等, 2013),表明维拉斯托矿区在晚石炭世先有一期花岗质岩浆侵入,后有一期闪长质岩浆侵入。薛怀民等(2010)报道了维拉斯托矿区石英闪长岩和闪长岩的SHRIMP锆石U-Pb年龄,其形成年龄分别为310±2Ma和311±2Ma,亦表明在晚石炭世期间维拉斯托地区存在闪长质岩浆活动。燕山期有3种类型的岩浆侵入活动:(1)肉红色花岗岩,分布于矿区东北约7.5km处的达青牧场一带,呈几个小岩株零星出露;(2)浅灰色斑状花岗岩,主要分布于矿区东南部约8~9km的北大山地区,刘翼飞(2009)报到了该岩体的SHRIMP锆石U-Pb年龄为140±3Ma,其边缘相斑状二长花岗岩形成年龄为138±2Ma;(3)浅黄色中细粒花岗岩,分布于矿区西北部约13~14km的巴音高勒苏木西部,岩体局部和部分伟晶岩内可见天蓝色天河石。

图 2 维拉斯托地区地质图(据内蒙古维拉斯托矿业有限责任公司和内蒙古地质勘查有限责任公司, 2015修改) Fig. 2 Geological map of the Weilasituo area, showing locations of ore deposits

① 内蒙古维拉斯托矿业有限责任公司, 内蒙古地质勘查有限责任公司. 2015. 内蒙古自治区克什克腾旗维拉斯托矿区锡多金属矿勘探报告. 1-235

2 维拉斯托稀有金属-锡多金属矿床和岩体地质特征 2.1 维拉斯托稀有金属-锡多金属矿床地质特征

维拉斯托大型稀有金属-锡多金属矿床位于内蒙古克什克腾旗政府所在地北约100km处(图 1中的5号矿床),矿区地理坐标:44°04′00″~44°05′45″N、117°28′15″~117°30′45″E。目前已控制锡金属量8.98万吨,平均品位0.80%,锌金属量8.00万吨,平均品位0.72%,三氧化钨金属量1.33万吨,平均品位0.44%,钼金属量0.03万吨,平均品位0.13%(刘瑞麟等, 2018b; 内蒙古自治区地质调查院,2018);另外,在隐爆角砾岩中发现大量的含锂云母,其Li2O品位在0.8%~3.6%之间,平均1.25%,估算Li2O金属量超过60万吨,在斑状细粒碱长花岗岩体的顶部还存在铷、铌、钽等成矿元素(刘瑞麟等,2018b)。

① 内蒙古自治区地质调查院. 2018. 中国矿产地质志内蒙古卷(有色金属矿产篇)报告. 1-678

维拉斯托矿床主体属于隐伏矿床,地表仅出露少量的NNE向和近EW向脉状矿体,大量的矿体均隐伏在地表以下50~750m范围内(图 3)。该矿床发育3种类型的矿体:(1)岩体顶部的浸染状锡锌(铷铌钽)矿体;(2)隐爆角砾岩筒中的浸染状-细脉状锡锌锂(铷铍)矿体;(3)断裂构造控制的石英脉型锡锌(钨钼铜)矿体(图 3b)。浸染状锡锌(铷铌钽)矿体位于矿床深部,矿体呈帽状产于斑状细粒碱长花岗岩体的顶部,矿体在岩体顶部厚度较大,在两侧逐渐变薄(图 3b)。矿石矿物主要为锡石和闪锌矿,除了锡和锌主成矿元素外,尚伴生铷、铌和钽,矿体中的锡石和闪锌矿常常一起产出形成团块,多呈稀疏浸染状分布于岩体顶部,团块大小不一,直径多集中在0.3~1.5cm,最大者可达10cm;脉石矿物主要为钠长石、天河石和黄玉。隐爆角砾岩型锡锌锂(铷铍)矿体位于斑状细粒碱长花岗岩岩体上方,呈陡倾筒状,地表出露范围约200m×300m(图 3)。角砾主要由黑云母斜长片麻岩、斜长角闪片麻岩、石英闪长岩组成,大小不一,棱角分明;胶结物为富锂云母云英岩,主要为锂云母,含少量的黄玉、石英和萤石。石英脉型锡锌(钨钼铜)矿体位于矿床浅部(图 3b),是维拉斯托矿床最重要的工业锡多金属矿体。该类型矿体的矿石品位较高,锡平均含量约1%。矿体总体呈舒缓的波状,倾向SE110°,倾角约40°(图 3b)。目前控制矿脉最大长度约800m,倾向延伸约700m,矿脉宽约1m,最大宽度可达8m。石英脉型矿体以锡为主,共伴生锌、钼、铜等,除了共伴生锡锌(钨、钼、铜)矿体外,在局部地段亦可圈出少量的独立锌矿体和钼矿体(图 3b)。该类型矿体的主要矿石矿物为锡石、闪锌矿、毒砂,少量辉钼矿、黑钨矿、磁黄铁矿、黄铁矿、黄铜矿、黝锡矿等;脉石矿物为石英,少量萤石、白云母和方解石等。

图 3 维拉斯托矿床地质图(a)和A-A′勘探线剖面图(b)(据内蒙古维拉斯托矿业有限责任公司和内蒙古地质勘查有限责任公司, 2015修改) Fig. 3 Geological map (a) and geological cross-section along the exploration line A-A′ (b) of the Weilasituo ore deposit
2.2 北大山岩体和维拉斯托岩体地质特征 2.2.1 北大山岩体

北大山岩体出露于维拉斯托稀有金属-锡多金属矿区东南约5km处,呈岩基产出,侵入中上侏罗统地层中(图 2)。岩体由中粗粒、中粒和细粒花岗岩组成,岩石具花岗结构、块状构造。花岗岩的主要矿物为碱性长石(40%~45%)、石英(40%~45%)、斜长石(5%~10%)和黑云母(5%~7%),另有少量白云母(图 4a-c),副矿物主要为锆石和磷灰石。斜长石一般呈自形板状或长板状,有的具有微弱的环带构造;钾长石多呈半自形板状-他形粒状,局部见碳酸盐化;石英他形粒状,部分石英发育裂纹;黑云母呈叶片状。

图 4 北大山岩体中-粗粒花岗岩(a-c)和维拉斯托岩体斑状细粒碱长花岗岩(d-i)显微镜下照片 (a)粗粒花岗岩,主要由石英和钾长石组成,其次为黑云母和斜长石;(b)中粒花岗岩,主要由石英和钾长石组成,其次为黑云母和斜长石;(c)中粒花岗岩,见铁叶云母;(d)斑状细粒碱长花岗岩,发育石英、微斜长石、钠长石、黄玉斑晶;(e)斑状细粒碱长花岗岩,发育石英、钾长石、钠长石、锂云母斑晶;(f)斑状细粒碱长花岗岩,发育钾长石、微斜长石、钠长石、黄玉、白云母斑晶,锂云母中见锆石副矿物;(g)斑状细粒碱长花岗岩的雪球结构,石英斑晶具有明显的两圈距离较远矿物生长带;(h)斑状细粒碱长花岗岩的雪球结构,石英斑晶具有单圈矿物生长带;(i)斑状细粒碱长花岗岩的雪球结构,石英斑晶雪球中包含钠长石矿物颗粒. Ab-钠长石;Bi-黑云母;Cbn-碳酸盐;Ep-绿帘石;Kf-钾长石;Lep-锂云母;Mic-微斜长石;Ms-白云母;Pl-斜长石;Qz-石英;Spl-铁叶云母;Toz-黄玉;Zrn-锆石. (c)为单偏光下照片,其他均为正交偏光下照片 Fig. 4 Representative photomicrographs for medium-coarse grained granite of the Beidashan pluton (a-c) and porphyritic fine-grained alkali-feldspar granite of the Weilasituo stock (d-i) (a) coarse-grained granite, mainly consisting of quartz and K-feldspar, followed by biotite and plagioclase; (b) medium-grained granite, mainly consisting quartz and K-feldspar, followed by biotite and plagioclase; (c) medium-grained granite, developing siderophyllite; (d) porphyritic fine-grained alkali-feldspar granite, developing quartz, microcline, albite, and topaz phenocrysts; (e) porphyritic fine-grained alkali-feldspar granite, developing quartz, K-feldspar, albite, and lepidolite phenocrysts; (f) porphyritic fine-grained alkali-feldspar granite, developing K-feldspar, microcline, albite, topaz, and muscovite phenocrysts, and zircon occurring in lepidolite; (g) snowball texture of porphyritic fine-grained alkali-feldspar granite, and two distinct mineral growth zones with a long distance distributing in quartz phenocrysts; (h) snowball texture of porphyritic fine-grained alkali-feldspar granite, and quartz phenocrysts developing monocyclic mineral growth zone; (i) snowball texture of porphyritic fine-grained alkali-feldspar granite, and quartz phenocryst snowball containing albite mineral particles. Ab-albite; Bi-biotite; Cbn-carbonate; Ep-epidote; Kf-K-feldspar; Lep-lepidolite; Mic-microcline; Ms-muscovite; Pl-plagioclase; Qz-quartz; Spl-siderophyllite; Toz-topaz; Zrn-zircon. (c) under plane polarized light, others under orthogonal polarized light
2.2.2 维拉斯托岩体

维拉斯托岩体产于维拉斯托稀有金属-锡多金属矿床深部,呈岩株状侵入中元古界宝音图群。斑状细粒碱长花岗岩为浅灰白色,局部地段因富天河石而整体呈浅天蓝色。岩石具似斑状结构,块状构造。斑晶包括石英、钾长石、钠长石,斑晶呈多级,粒度差别极大;基质为细粒结构,主要矿物为钠长石、钾长石、石英、黄玉、锂云母和萤石(图 4d-f)。石英斑晶呈短柱状,边部因熔蚀呈浑圆状,粒度2~8mm,含量40%,石英斑晶内部包含早期结晶的细粒钠长石,石英斑晶边缘含大量细粒半自形钠长石、他形钾长石,沿结晶面生长,构成一圈或多圈环带,形成典型的“雪球结构(snowball texture)”(图 4g-i)。钾长石斑晶粒度稍细,粒径1~2mm,含量30%,自形晶为主,多表现为微斜长石,少量条纹长石,包裹少量钠长石,发育有格子双晶或条纹结构(图 4d),与石英边部的似文象结构相似。钠长石,除早阶段包裹于石英、钾长石中的钠长石外,大量钠长石呈细粒状分布于石英、钾长石斑晶之间,粒径0.05~0.5mm,粒度大体介于斑晶矿物与基质矿物之间,含量20%,绝大部分钠长石呈他形晶。基质包括细粒石英、长石、黄玉、萤石、锂云母等,总含量小于20%。其中,黄玉、石英分两类,一部分呈微细粒状分布于基质中,另一类局部集中呈不规则状。

3 样品和分析方法

用于主量、微量和稀土元素分析的样品包括9件北大山岩体中粗粒花岗岩(样品BD1-1~BD1-9)和10件维拉斯托岩体斑状细粒碱长花岗岩(样品ZK16、ZK17、ZK19、ZK21、ZK37、ZK38、ZK43、ZK44、ZK46和ZK47);选取1件中粗粒花岗岩(样品BD1-1)和2件斑状细粒碱长花岗岩(样品ZK44和ZK46)进行LA-ICP-MS锆石U-Pb定年和Lu-Hf同位素分析。

北大山岩体的9件样品全部取自维拉斯托矿床东南侧5km处岩体地表露头,采样剖面见图 2,样品特征见前文北大山岩体岩相学描述。

维拉斯托岩体10件斑状细粒碱长花岗岩样品全部取自维拉斯托稀有金属-锡多金属矿床的钻孔中。样品ZK16、ZK17、ZK19、ZK21、ZK37、ZK38、ZK43、ZK44、ZK46和ZK47分别采自钻孔IZK02307的679~682m处、ZK03104的590~593m处、IZK03105的521~524m处、IZK02306的500~503m处、IZK03901的586~589m处、IZK03901的553~563m处、ZK03101的700~703m处、IZK2305的802~799m处、IZK1509的1075~1078m处和IZK01505的514~517m处。部分样品的采样钻孔见图 3b,样品特征见前文的维拉斯托岩体岩相学特征描述。

3.1 LA-ICP-MS锆石U-Pb定年

锆石由河北省廊坊诚信地质服务公司采用常规重力和磁选方法分选,为保证各粒级锆石颗粒的完整性,将重约5~10kg的岩石经过粗碎及中碎,将粉碎后的样品分粒级进行淘洗,然后在双目镜下挑选出锆石,最后,将锆石置于环氧树脂中,磨制约一半大小,使锆石内部暴露,用于阴极发光和LA-ICP-MS U-Pb分析。锆石阴极发光研究在北京锆年领航科技有限公司完成,锆石LA-ICP-MS U-Pb分析在北京燕都中实测试技术有限公司Finnigan Neptune型多接收等离子体质谱仪上完成,该仪器连接New Wave UP型213nm激光剥蚀系统。样品分析流程及原理参见Griffin et al. (2004)侯可军等(2009)。用参考锆石GEMOC GJ-1(207Pb/206Pb年龄为608.5±1.5Ma,Jackson et al., 2004)进行元素间的分馏校正,用标准样品Mud Tank(截距年龄为732±5Ma, Black and Gulson, 1978)控制测试精度,年龄数据处理采用ISOPLOT 3.0程序(Ludwig, 2003)。因为低的信号强度和204Hg气体的干扰,不能测试204Pb值。普通Pb校正采用EXCEL程序ComPbCorr#3 15G(Andersen, 2002)。单个点的分析误差是1σ水平,整体年龄误差在2σ水平上是95%。

3.2 锆石Lu-Hf同位素分析

锆石原位Lu-Hf同位素分析在中国地质科学院矿产资源研究所自然资源部成矿作用与资源评价重点实验室进行,使用的仪器为Neptune型多接收等离子体质谱仪,配以ArF激发激光剥蚀系统(193nm)。仪器工作条件和数据采集见Yuan et al.(2004, 2018a)。测试时,激光束斑直径为44μm,激光重复率为10Hz,能量为100mJ。标准锆石Mud的176Hf/177Hf比值为0.282500±0.000030。用测试的176Hf/177Hf和176Lu/177Hf比值计算初始的176Hf/177Hf比值,采用1.867×10-11/y作为176Lu的衰变常数(Söderlund et al., 2004)。εHf(t)值计算时采用的球粒陨石现代176Hf/177Hf和176Lu/177Hf比值分别0.282772和0.0332(Blichert-Toft and Albarède, 1997)。计算Hf同位素模式年龄时采用如下参数:亏损地幔的176Hf/177Hf和176Lu/177Hf现代比值分别为0.28325和0.0384(Blichert-Toft and Albarède, 1997);大陆上地壳平均的176Lu/177Hf值为0.0093(Vervoort and Blichert-Toft, 1999)。具体的计算过程见Nowell et al. (1998)Griffin et al. (2000)描述。

3.3 全岩主量、稀土和微量元素分析

样品粉末样磨制在河北省廊坊诚信地质服务有限公司完成。将样品粗碎至厘米级的块体,用纯化水冲洗干净,烘干并粉碎至200目以备使用。主量、稀土和微量元素分析在核工业北京地质研究院实验室完成。主量元素分析采用AxiosmAX型X射线荧光光谱仪,配以高色散光栅系统,分析误差根据美国地质调查局岩石标样BCR-1和AVG-2及中国国家岩石标样GSR-3确定,TiO2和P2O5的分析精度分别约为1.5%和2.0%,其他氧化物的分析精度优于1%。稀土和微量元素分析采用Agilent-7500a型电感耦合等离子体质谱仪,分析数据的质量用美国地质调查局的岩石标样BCR-1和BHVO-1监控,绝大多数稀土和微量元素的分析精度优于5%。

4 测试结果 4.1 LA-ICP-MS锆石U-Pb定年

对1件北大山岩体中粒花岗岩样品中的9粒锆石和2件维拉斯托稀有金属-锡多金属矿床内的斑状细粒碱长花岗岩样品中的23粒锆石开展了LA-ICP-MS锆石U-Pb定年。3件样品的锆石分析点位和它们的CL图像见图 5,其分析结果见表 1,锆石年龄谐和图见图 6

图 5 北大山岩体中粗粒花岗岩(a)和维拉斯托岩体斑状细粒碱长花岗岩(b、c)锆石阴极发光图,示测点位置及对应年龄值 Fig. 5 Representative cathodoluminescence (CL) images of zircons from medium-coarse grained granite of the Beidashan pluton (a) and porphyritic fine-grained alkali-feldspar granite of the Weilasituo stock (b, c), showing U-Pb analytical spots and corresponding 206Pb/238U ages

表 1 北大山岩体(BD1-1)和维拉斯托岩体(ZK44和ZK46)LA-ICP-MS锆石U-Pb分析结果 Table 1 LA-ICP-MS zircon U-Pb data for the Beidashan pluton (Sample BD1-1) and the Weilasituo stock (samples ZK44 and ZK46) in the Weilasituo area

图 6 北大山岩体(a、b)和维拉斯托岩体(c-f)的锆石U-Pb年龄谐和图及加权平均年龄 Fig. 6 Zircon U-Pb concordia diagrams and age histograms for the Beidashan pluton (a, b) and the Weilasituo stock (c-f)
4.1.1 北大山岩体花岗岩

北大山岩体中粗粒花岗岩样品BD1-1中的锆石呈无色-浅黄色。锆石阴极发光图像显示,该样品中的锆石晶体形态较好,主要呈单锥或双锥状,多数锆石具生长韵律环带(图 5a),表明其岩浆成因。表 1显示,锆石的U和Th含量分别介于804×10-6~7714×10-6和356×10-6~2413×10-6。经计算,锆石的Th/U比值介于0.25~0.82之间,亦表明其为岩浆成因(Pupin, 1980)。9粒锆石的10个分析点集中分布在谐和线上(图 6a),其206Pb/238U年龄的加权平均值为140±2Ma(MSWD=0.10, n=10; 图 6b)。表明北大山岩体的结晶年龄为140±2Ma。

4.1.2 维拉斯托岩体斑状细粒碱长花岗岩

锆石的阴极发光显示,2件样品的锆石内部结构不清楚,颜色普遍偏暗,个别锆石出现明亮白色(如ZK46-3.1),无法观察到振荡环带(图 5b, c),具有明显的高U特征。锆石的U含量极高(9732×10-6~50853×10-6),这符合多数与钨锡矿成矿相关的高分异花岗岩岩体普遍具有高U含量的特征(Wang et al., 2016b)。在单偏光镜下,该岩体锆石呈淡黄色,透明,少量半透明,裂纹不发育,表面光滑,具金刚光泽,边界平整,为柱状、长柱状自形晶,锆石大小100~300μm,晶体长宽比为1~4。一般情况下,典型的岩浆成因锆石Th/U比值一般为0.1~1.0,且U和Th之间具有明显的正相关关系(Belousova et al., 2002)。该岩体锆石的Th/U比为0.03~0.28,多数小于0.1,可能暗示锆石受到了流体改造。

样品ZK-44中锆石的Th、U含量分别为875×10-6~5270×10-6(平均2585×10-6)和11687×10-6~50853×10-6(平均27568×10-6),Th/U值介于0.03~0.24之间,绝大多数小于0.1,平均0.11。11粒锆石的12个分析点集中分布在谐和线上(图 6c),其206Pb/238U年龄的加权平均值为138±1Ma,MSWD值为0.76(图 6d)。样品ZK-46锆石的Th、U含量分别为835×10-6~3031×10-6(平均1772×10-6)和9732×10-6~45926×10-6(平均19365×10-6),Th/U值介于0.03~0.28之间,绝大多数小于0.1,平均0.11。12粒锆石的16个分析点集中分布在谐和线上及其附近(图 6e),其206Pb/238U年龄的加权平均值为137±1Ma,MSWD值为0.32(图 6f)。表明维拉斯托岩体的结晶年龄为138~137Ma。

4.2 全岩主量、稀土和微量元素分析

北大山岩体9件中粗粒花岗岩和维拉斯托岩体10件斑状细粒碱长花岗岩样品的主量、微量和稀土元素分析结果见表 2

表 2 北大山和维拉斯托岩体全岩主量元素(wt%)和微量元素(×10-6)分析结果 Table 2 Major (wt%) and trace (×10-6) elements analytical results for the Beidashan pluton and the Weilasituo stock in the Weilasituo area
4.2.1 北大山岩体

北大山岩体9件中粗粒花岗岩样品的SiO2含量为74.40%~76.47%,Al2O3含量为12.42%~13.44%,Fe2O3含量为0.18%~0.43%,FeO含量为0.89%~1.23%;岩石具有低的TiO2、MnO、MgO、CaO和P2O5含量;Na2O含量为3.81~4.37%,K2O含量为4.57%~5.45%,Na2O/K2O比值为0.72~0.96;铝饱和指数为0.99~1.06(表 2)。在花岗岩分类的QAP图解中(Streckeisen, 1976),这些中粗粒花岗岩属于碱长花岗岩(图 7a);在R1-R2图解(De La Roche et al., 1980)中,所有样品均落入碱性花岗岩区域(图 7b);在Maniar and Piccoli (1989)的A/CNK-A/NK图解中,几乎所有的样品均落入弱过铝质区域(图 7c);在Peccerillo and Taylor (1976)的SiO2-K2O图解中,所有的样品均落入高钾钙碱性系列区域(图 7d)。

图 7 北大山和维拉斯托岩体岩石分类图解 (a)QAP图解(Streckeisen, 1976);(b)R1-R2图解(De La Roche et al., 1980);(c)A/CNK-A/NK图解(Maniar and Piccoli, 1989);(d)SiO2-K2O图解(Peccerillo and Taylor, 1976) Fig. 7 Diagrams of petrographic classification for the Beidashan pluton and the Weilasituo stock (a) quartz-alkali feldspar-plagioclase (QAP) diagram (Streckeisen, 1976); (b) R1 vs. R2 diagram (De La Roche et al., 1980); (c) A/CNK vs. A/NK diagram (Maniar and Piccoli, 1989); (d) SiO2 vs. K2O diagram (Peccerillo and Taylor, 1976)

表 2可知,北大山岩体中粗粒花岗岩具有中-高的稀土总量(206×10-6~436×10-6),相对高的Yb含量(5.36×10-6~20.3×10-6),强烈的铕负异常(Eu/Eu*值为0.01~0.04)。轻重稀土分馏程度较低,其(La/Yb)N值介于2.24~5.17之间。在球粒陨石标准化的稀土元素配分曲线上,这些岩石相对富集轻稀土、亏损重稀土,轻稀土分馏相对弱,重稀土分馏不明显,稀土配分曲线为轻稀土略富集型(图 8a)。

图 8 北大山和维拉斯托岩体球粒陨石标准化稀土元素配分图(a, 标准化值据Boynton, 1984)和原始地幔标准化微量元素蛛网图(b,标准化值据Sun and McDonough, 1989) Fig. 8 Chondrite-normalized REE patterns (a, normalization values after Boynton, 1984) and primitive mantle-normalized trace element spider grams (b, normalization values after Sun and McDonough, 1989) of the Beidashan pluton and the Weilasituo stock

在原始地幔标准化的微量元素蛛网图上,北大山岩体明显富集Rb、Th、U、Ta、Hf、Zr、Y和稀土元素,强烈亏损Ba、Sr、P、Eu和Ti,弱亏损K和Nb(图 8b)。另外,这些岩石具有极低的Sr含量(5.73×10-6~18.3×10-6,平均13.5×10-6)和高的Y含量(37.2×10-6~127×10-6,平均84.0×10-6)。

4.2.2 维拉斯托岩体

维拉斯托岩体10件斑状细粒碱长花岗岩样品的SiO2含量为70.77%~73.87%,TiO2含量为0.01%~0.03%,Al2O3含量为15.08%~17.07%;岩石具有低的Fe2O3、FeO、MnO、MgO、CaO和P2O5含量,相对低的K2O含量(2.90%~4.01%),高的Na2O含量(6.11%~7.59%)。Na2O/K2O比值为1.61~2.58,铝饱和指数为0.97~1.16(平均1.06)。在QAP分类图解中,绝大多数样品落入碱长花岗岩区域,少数样品落入石英碱长正长岩区域(图 7a)。在R1-R2图解中,绝大多数落入石英正长岩区域,少数样品落入碱性花岗岩区域(图 7b)。在A/CNK-A/NK图解中,绝大多数样品落入弱过铝质区域,个别样品落入准铝质和强过铝质区域(图 7c)。在SiO2-K2O图解中,绝大多数样品落入高钾钙碱性系列区域,有少部分样品落入钙碱性系列区域(图 7d)。

表 2可知,早白垩世斑状细粒碱长花岗岩具有低的稀土总量(28.8×10-6~50.2×10-6),相对高的Yb含量(3.71×10-6~6.17×10-6),强烈的铕负异常(Eu/Eu*值为0.01~0.07)。轻重稀土分馏弱,其(La/Yb)N值介于0.37~0.85之间。稀土元素四分组效应明显,展示出典型的“M”型分布特征,即(La、Ce、Pr、Nd)、(Pm、Sm、Eu、Gd)、(Gd、Tb、Dy、Ho)和(Er、Tm、Yb、Lu)构成四条上凸曲线,其中Ce、Sm、Dy、Yb分别为四条曲线的最高点,其中Eu由于强烈亏损而偏离曲线(图 8a)。

在原始地幔标准化的微量元素蛛网图上,这些早白垩世斑状细粒碱长花岗岩明显富集Rb、U、Ta、La、Ce、Hf、Gd、Tb、Dy、Yb和Lu,强烈亏损Ba、Sr、Eu、Ti和Y,弱亏损K和Nb(图 8b)。另外,这些岩石有极低的Sr含量(1.54×10-6~13.8×10-6,平均4.38×10-6)和Y含量(1.44×10-6~4.17×10-6,平均2.83×10-6)。

4.3 锆石Lu-Hf同位素分析

在已经定年的锆石点上开展锆石原位Lu-Hf同位素分析,其分析结果见表 3图 9

表 3 北大山岩体(BD1-1)和维拉斯托岩体(ZK44和ZK46)锆石的Lu-Hf同位素分析结果 Table 3 Zircon Lu-Hf isotopic data for the Beidashan pluton (Sample BD1-1) and the Weilasituo stock (samples ZK 44 and ZK46)

图 9 北大山和维拉斯托岩体锆石的Hf同位素图解 中亚造山带东段显生宙火成岩Hf同位素组成据Xiao et al.(2004)Chen et al.(2009);燕山褶皱逆冲带显生宙火成岩Hf同位素组成据Yang et al.(2006) Fig. 9 Correlation diagrams between Hf isotopic compositions and ages of zircons from the Beidashan pluton and the Weilssituo stock, respectively Hf isotopic composition of Phanerozic igneous rocks in the eastern Central Asian orogenic belt after Xiao et al. (2004) and Chen et al. (2009), and Hf isotopic composition of Phanerozic igneous rocks in the Yanshan fold and thrust belt after Yang et al. (2006)
4.3.1 北大山岩体

北大山碱长花岗岩样品BD1-1中9粒锆石的9个分析点的Hf同位素初始176Hf/177Hf值介于0.282839~0.282929之间,平均0.282875。它们的εHf(t)值介于5.4~8.6之间,平均6.7。Hf同位素单阶段和两阶段模式年龄分别为467~579Ma和643~847Ma,均值分别为534Ma和765Ma(表 3)。

4.3.2 维拉斯托岩体

维拉斯托斑状细粒碱长花岗岩2件样品中22粒锆石的25个分析点的初始176Hf/177Hf值介于0.282807~0.282968之间,平均0.282860。它们的εHf(t)值介于4.2~9.9之间,平均6.1。Hf同位素单阶段和两阶段模式年龄分别为402~647Ma和556~921Ma,均值分别为578Ma和799Ma(表 3)。

5 讨论 5.1 维拉斯托岩体与北大山岩体成因联系

维拉斯托地区广泛发育早白垩世花岗岩,除了矿区东南部的北大山岩体外,维拉斯托矿区的北东、南西和北西部均出露早白垩世花岗岩(图 2)。这些岩石的结构变化较大,既有中粗粒结构、中粒结构,也有中细粒结构,似斑状花岗岩亦很常见,这些岩体在深部可能相连,构成早白垩世花岗岩基。刘翼飞(2009)获得北大山岩体中心部位的中粗粒花岗岩的SHRIMP锆石U-Pb年龄为140±3Ma,其边缘相斑状细粒花岗岩形成年龄为138±2Ma;祝新友等(2016)报到了维拉斯托矿区深部细粒碱长花岗岩的LA-ICP-MS锆石U-Pb年龄为139.5±1.2Ma。本次研究获得北大山岩体中粒花岗岩的LA-ICP-MS锆石U-Pb年龄为140±2Ma,维拉斯托岩体斑状细粒碱长花岗岩的锆石U-Pb年龄为137~138Ma。本文的定年结果和前人获得的年龄数据均表明北大山岩体和维拉斯托岩体为同一期岩浆活动的产物。维拉斯托岩体与北大山岩体密切的时空关系暗示两者可能存在成因联系。在以SiO2含量为横坐标的一系列哈克图解中(图 10),2个岩体的主要氧化物和代表岩浆演化的Rb、Sr和Ba与SiO2具有明显的线性关系。从北大山岩体到维拉斯托岩体,TiO2、FeOT、CaO、K2O、Sr和Ba含量降低,Al2O3、MgO、Na2O、全碱、P2O5和Rb含量增加,但2个岩体之间亦具有明显的成分间断,暗示形成北大山岩体的岩浆经历了不止一次的分离结晶作用,北大山岩体与维拉斯托岩体之间可能还存在分异程度介于两者之间的岩体,只是由于剥蚀和保存状况及矿床勘探深度,尚不能确定中间分异岩体的位置。综上,我们认为维拉斯托岩体和北大山岩体是同一期岩浆活动的产物,北大山岩体结晶略早,形成北大山岩体的残余岩浆经过不止一次的结晶分异,最终形成了具有斑状结构的维拉斯托细粒碱长花岗岩。

图 10 北大山和维拉斯托岩体哈克图解 Fig. 10 Harker diagrams of the Beidashan pluton and the Weilasituo stock
5.2 岩石成因类型和岩浆源区组成 5.2.1 岩石成因类型

北大山岩体和维拉斯托岩体的岩石类型均为碱长花岗岩,两者均具有高硅,富碱,贫铁、镁、钙、钛和磷的特征,主体均属高钾钙碱性的弱过铝质岩石。岩石展示了强烈的Eu负异常,轻重稀土分馏均较弱。与北大山岩体相比,维拉斯托岩体具有更低的稀土含量,更为明显的Eu负异常,Ba、Sr、Ti和Y亏损更强烈,另外,维拉斯托岩体展示了明显的稀土元素四分组效应(图 8)。岩相学研究显示,2个岩体均不存在典型的S型花岗岩中特有的富铝矿物(如堇青石和石榴石),结合其非常低的P2O5含量和弱过铝质岩石特征,我们首先排除它们为S型花岗岩的可能性。尽管2个岩体均展示了A型花岗岩的大部分地球化学特征,但其Zr含量明显低于A型花岗岩(>250×10-6, Whalen et al., 1987),它们的锆饱和温度介于568~726℃之间,明显低于A型花岗岩的锆饱和温度(平均839℃, King et al., 1997),亦低于正常I型花岗岩的锆饱和温度(平均781℃, King et al., 1997);此外,2个岩体均具有高的分异指数(DI=94.5~98.4)。因此,我们认为北大山岩体和维拉斯托岩体均不是A型花岗岩,而是分异的I型花岗岩。Whalen et al. (1987)提出了一系列花岗岩成因类型判别图,在10000Ga/Al-(Na2O+K2O)判别图中,2个岩体样品全部投影在A型花岗岩区域(图 11a)。在10000Ga/Al-Zr判别图中,这2个岩体样品亦全部落在A型花岗岩区域(图 11b)。吴福元等(2017)指出,尽管A型花岗岩也存在分异作用,但它在10000Ga/Al-Zr判别图上的趋势是从A型花岗岩区向高分异花岗岩区演化;相反,I型或S型花岗岩在分异的过程中,其10000Ga/Al比值逐渐升高,显示与A型花岗岩相反的趋势(图 11b)。北大山岩体与维拉斯托岩体的演化关系明显不同于A型花岗岩的分异趋势。在(Zr+Nb+Ce+Y)-10000Ga/Al判别图中,北大山岩体投影在未分异的I、S、M型花岗岩、分异型花岗岩和A型花岗岩三者的过渡区;维拉斯托岩体的所有样品均落入分异型花岗岩区域(图 11c)。在(Zr+Nb+Ce+Y)-(Na2O+K2O)/CaO判别图中,北大山岩体投影在分异型花岗岩与A型花岗岩的过渡区;维拉斯托岩体投影在分异型花岗岩区域的上方(图 11d)。综上,我们认为,北大山岩体和维拉斯托岩体均为分异的I型花岗岩。根据吴福元等(2017)厘定的高分异花岗岩和超分异花岗岩划分标准,我们认为北大山岩体属于高分异花岗岩,维拉斯托岩体属于超分异花岗岩。

图 11 花岗岩成因类型判别图解(据Whalen et al., 1987; b中的演化趋势线据吴福元等, 2017) Fig. 11 Chemical discrimination diagrams for types of granites (after Whalen et al., 1987; the evolution trend line in Fig. 11b after Wu et al., 2017)
5.2.2 岩浆源区

北大山岩体的锆石εHf(t)值介于5.4~8.6之间,平均6.7;两阶段Hf模式年龄介于643~847Ma之间,均值为765Ma。维拉斯托岩体的锆石εHf(t)值介于4.2~9.9,平均6.1;两阶段Hf模式年龄介于556~921Ma之间,均值为799Ma。维拉斯托岩体与北大山岩体具有相似的锆石Hf同位素组成,在εHf(t)-年龄图解中(图 9),2个岩体的样品均落于球粒陨石和亏损地幔演化线之间,且位于Xiao et al. (2004)Chen et al. (2009)限定的中亚造山带东段显生宙火成岩Hf同位素组成区域。由于幔源岩浆结晶分异难以形成大规模花岗岩(吴福元等, 2007),故上述特征表明维拉斯托地区的早白垩世岩体的源区主要为起源于亏损地幔的新元古代新生地壳物质。结合大兴安岭南段黄岗(Zhou et al., 2012)、东山湾(Zeng et al., 2015)、道伦达坝(陈公正, 2018)、白音查干东山(李睿华, 2019)、毛登-小孤山(李睿华, 2019)、宝盖沟(Mi et al., 2020)、查木罕(Zhang et al., 2019)等锡多金属矿床的成矿岩体锆石Hf同位素组成,我们推测新元古代期间大兴安岭南段存在一次重要的地壳生长事件,新元古代新生地壳物质部分熔融形成的花岗质岩浆是大兴安岭南段锡多金属矿床形成的重要因素之一。

5.3 成岩构造背景

北大山岩体和维拉斯托岩体的岩石类型主要为碱长花岗岩,岩石均属高钾钙碱性系列。在Pearce et al. (1984)的Yb-Ta判别图中,北大山岩体和维拉斯托岩体均落入板内花岗岩区域(图 12a)。在Pearce et al. (1984)的(Yb+Ta)-Rb判别图中,2个岩体的样品均落入板内花岗岩与同碰撞花岗岩过渡区(图 12b)。在Brown et al. (1984)的SiO2-log[CaO/(Na2O+K2O)]判别图中,北大山岩体样品全部落入伸展环境形成的花岗岩区域,维拉斯托岩体落入伸展型的下方,但靠近伸展环境形成的花岗岩区域(图 12c)。在Batchelor and Bowden (1985)推荐的R1-R2判别图中,北大山岩体落在晚造山花岗岩和非造山花岗岩过渡区域,维拉斯托岩体投影在非造山花岗岩区域,但亦靠近晚造山花岗岩区域(图 12d)。北大山岩体和维拉斯托岩体的岩石类型、地球化学特征及构造环境判别图表明它们形成于伸展环境。

图 12 北大山和维拉斯托岩体构造环境判别图 (a)Yb-Ta判别图(Pearce et al., 1984);(b)(Yb+Ta)-Rb判别图(Pearce et al., 1984);(c)SiO2-log[CaO/(Na2O+K2O)]判别图(Brown et al., 1984);(d)R1-R2判别图(Batchelor and Bowden, 1985). ORG-大洋中脊花岗岩;syn-COLG-同碰撞花岗岩;VAG-火山弧花岗岩;WPG-板内花岗岩;Ⅰ-地幔分异的花岗岩;Ⅱ-板块碰撞前的花岗岩;Ⅲ-碰撞后隆升的花岗岩;Ⅳ-晚造山花岗岩;Ⅴ-非造山花岗岩;Ⅵ-同碰撞花岗岩;Ⅶ-后造山花岗岩 Fig. 12 Tectonic discrimination diagrams for the Beidashan pluton and the Weilasituo stock (a) Yb vs. Ta diagram (Pearce et al., 1984); (b) Yb+Ta vs. Rb diagram (Pearce et al., 1984); (c) SiO2 vs. log[CaO/(Na2O+K2O)] diagram (Brown et al., 1984); (d) R1 vs. R2 diagram (Batchelor and Bowden, 1985). ORG-oceanic ridge granites; syn-COLG-syn-collisional granites; VAG-volcanic arc granites; WPG-within-plate granites; Ⅰ-mantle differentiated granites; Ⅱ-pre-collisional granites; Ⅲ-post-collisional uplifting granites; Ⅳ-late-orogenic granites; Ⅴ-anorogenic granites; Ⅵ-syn-collisional granites; Ⅶ-post-orogenic granites

晚侏罗世-早白垩世期间,包括大兴安岭南段在内的中国东北地区主要受东部的古太平洋构造体制(Wu et al., 2011)和北部的蒙古-鄂霍茨克洋构造体制(Ying et al., 2010; Wang et al., 2012)控制。关于早白垩世期间伸展的地球动力学背景,目前主要存在以下几种观点:(1)受古太平洋板块向欧亚大陆下俯冲的远程效应影响(毛景文等, 2005);(2)受蒙古-鄂霍茨克海闭合后造山带垮塌形成的伸展环境影响(Davis et al., 1998; Meng, 2003);(3)受蒙古-鄂霍茨克海闭合后的造山后伸展环境及古太平洋板块向欧亚大陆下俯冲形成的弧后伸展环境的叠加影响(Wang et al., 2012; Xu et al., 2013; Chen et al., 2017; Wu et al., 2017);(4)受由EW向演变为NNE向的中国东部构造体制转换影响(翟明国等, 2004; 毛景文等, 2013);(5)受板内非造山环境影响(Pirajno et al., 2009)。大兴安岭南段晚侏罗世-早白垩世火成岩主要呈NNE向展布,其展布方向与NEE向的蒙古-鄂霍茨克造山带斜交,因此,蒙古-鄂霍茨克造山带垮塌不可能是维拉斯托地区早白垩世花岗岩形成的主要原因。190Ma开始,小兴安岭-张广才岭地区形成了呈NE向展布的弧岩浆岩,表明古太平洋板块于早侏罗世开始向欧亚大陆之下俯冲(Yang and Li, 2008; Yu et al., 2012; Xu et al., 2013);145~130Ma,由于古太平洋板块的俯冲方向改变和/或俯冲板片后撤,中国东部整体进入伸展环境,发生大规模的岩浆活动与成矿事件(Wu et al., 2002; Liu et al., 2005; Wang et al., 2012; Mao et al., 2013, 2019, 2021)。因此,我们认为,古太平洋板块的俯冲方向改变和/或俯冲板片后撤导致大兴安岭地区处于伸展环境,在伸展背景下软流圈物质上涌,底侵于下地壳底部,加热新元古代期间形成的地壳,使其发生部分熔融,岩浆经过不同程度的演化,分别形成北大山岩体和维拉斯托岩体。

5.4 岩浆演化与成矿 5.4.1 岩浆演化过程

前已述及,北大山岩体和维拉斯托岩体都属于分异型花岗岩。随着岩浆演化,花岗岩通常会经历一个由Zr不饱和向Zr饱和演化过程(Watson and Harrison, 1983, 1984)。图 13显示,从北大山岩体到维拉斯托岩体,Zr含量逐渐降低,这一过程与岩浆演化趋势一致。随着Zr含量的降低,FeOT、TiO2、Sr含量和Zr/Hf、Nb/Ta、Y/Ho和Ba/Sr比值降低,Hf、Nb、Ta、Rb含量和Rb/Sr比值增加(图 13)。Ba倾向于进入钾长石,钾长石的分离结晶将导致残余岩浆中Ba含量和Ba/Sr比值降低(Li et al., 2007a, b)。Sr倾向于进入斜长石,斜长石的分离结晶导致残余岩浆中Sr含量降低,Ba/Sr比值升高(Li et al., 2007a, b)。Zr和Hf主要富集在锆石中,Zr、Hf含量和Zr/Hf比值的降低往往和Zr饱和之后的锆石分离结晶有关(Watson and Harrison, 1983, 1984)。黑云母的分离结晶会使残余岩浆富集Nb和Ta,而且会使岩浆的Ba/Sr和Nb/Ta比值降低(Stepanov et al., 2014)。角闪石分离结晶将增加残余岩浆中的Nb和Ta含量,并降低熔体的Nb/Ta比值(Tiepolo et al., 2000; Pfänder et al., 2007)。金红石、钛铁矿、榍石的分离结晶导致残余岩浆中Ti、Nb和Ta含量降低,Nb/Ta比值升高(Linnen and Keppler, 1997)。与北大山岩体相比,维拉斯托岩体的稀土元素含量降低,而且显示出明显的稀土四分组效应,暗示岩浆演化过程中发生了独居石、榍石、褐帘石和磷灰石等副矿物的分离结晶作用(Bau, 1996; Irber, 1999; Monecke et al., 2002; Gelman et al., 2014)。由于地幔物质具有较高的Zr/Hf和Nb/Ta比值(Zr/Hf约为34.3,Nb/Ta约为19.9,Münker et al., 2003)和较低的Zr和Nb含量,地幔物质的参与会使北大山岩体和维拉斯托岩体具有较高的Zr/Hf和Nb/Ta比值和低的Zr和Nb含量,而不会产生随着Zr含量降低Zr/Hf比值降低的趋势,也无法解释随着Nb/Ta比值降低,Nb和Ta含量升高的趋势。因此,北大山岩体和维拉斯托岩体Zr/Hf和Nb/Ta的分异不是地幔物质参与的结果。上地壳物质的Zr/Hf比值约为36.7(Rudnick and Gao, 2004)、Nb/Ta比值约为13.4(Rudnick and Gao, 2004),且具有低的Nb和Ta含量,上地壳物质的混染,也无法解释Nb和Ta含量随着Nb/Ta比值降低而增加的趋势。因此,岩浆的分离结晶是Zr/Hf和Nb/Ta分异的根本原因。随着岩浆分离结晶程度的增加,岩浆演化晚期的流体将导致Zr和Hf的置换比例增大,Hf含量升高,表明Zr/Hf比值的分异也与流体的参与有关(Zhang et al., 2008)。北大山岩体和维拉斯托岩体变化较小的锆石εHf(t)值和Hf同位素两阶段模式年龄,亦表明岩体的源区组成较为单一,而且未受到地幔和/或上地壳物质的混染。北大山岩体和维拉斯托岩体地球化学特征及岩浆演化趋势表明,从北大山岩体到维拉斯托岩体,岩浆演化过程中经历了钾长石、斜长石、黑云母、角闪石、金红石、钛铁矿、独居石、榍石、褐帘石、磷灰石和锆石的分离结晶;维拉斯托岩体显著的稀土四分组效应指示岩浆演化晚期发生了熔体-流体相互作用。

图 13 北大山和维拉斯托岩体中Zr含量对FeOT、TiO2、Hf、Nb、Ta、Rb、Sr、Zr/Hf、Nb/Ta、Y/Ho、Rb/Sr、Ba/Sr相关性图解 Fig. 13 Variations of Zr against FeOT, TiO2, Hf, Nb, Ta, Rb, Sr, Zr/Hf, Nb/Ta, Y/Ho, Rb/Sr, and Ba/Sr of the rocks from the Beidashan pluton and the Weilasituo stock

基于瑞利分馏模型,结合岩体中的微量元素含量及实验测定的相关元素在不同矿物相和熔体之间的分配系数进行定量模拟,模拟结果见图 14图 14显示,从北大山岩体到维拉斯托岩体,岩浆主要发生了斜长石、钾长石、黑云母、磷灰石、独居石和褐帘石的分离结晶作用,斜长石的分离结晶导致维拉斯托岩体具有更强的Eu和Sr的负异常;钾长石的分离结晶明显降低了维拉斯托岩体的钡含量;黑云母的分离结晶将降低残余熔体中的铁、钛、钡含量;富含轻稀土的副矿物独居石、榍石、褐帘石和磷灰石的分离结晶导致维拉斯托岩体轻稀土含量明显低于北大山岩体,维拉斯托岩体重稀土和钇含量明显低于北大山岩体,指示锆石和磷钇矿发生了分离结晶作用。另外,北大山岩体和维拉斯托岩体本身也发生了钾长石、斜长石、黑云母、磷灰石和独居石的分离结晶(图 14)。2个岩体均表现出Ti的负异常,可能与钛铁氧化物的分离结晶作用有关。

图 14 北大山和维拉斯托岩体微量元素定量模拟图解(据Xu et al., 2015修改) Ba、Sr和Eu分配系数据Rollinson (1993);La和Yb分配系数据Mahood and Hildreth (1983)Fujimaki (1986)Yurimoto et al. (1990). Aln-褐帘石;Amp-角闪石;Ap-磷灰石;Bt-黑云母;Grt-石榴石;Kf-钾长石;Mnz-独居石;Ms-白云母;Pl-斜长石;Zrn-锆石 Fig. 14 Numerical modelling diagrams for trace elements of the Beidashan pluton and the Weilasituo stock (modified after Xu et al., 2015) Partition coefficients of Ba, Sr, and Eu cited from Rollinson (1993); partition coefficients of La and Yb cited from Mahood and Hildreth (1983), Fujimaki (1986), and Yurimoto et al. (1990). Aln-allanite; Amp-amphibole; Ap-apatite; Bt-biotite; Grt-garnet; Kf-K-feldspar; Mnz-monozite; Ms-muscovite; Pl-plagioclase; Zrn-zircon

总之,北大山岩体和维拉斯托岩体地球化学特征及岩浆演化趋势表明,它们均属于分异型花岗岩,岩浆演化过程中经历了钾长石、斜长石、黑云母、独居石、榍石、褐帘石、磷灰石、锆石和钛铁矿的分离结晶;除了岩浆分离结晶,形成维拉斯托岩体的残余岩浆还发生了熔体-流体相互作用,展示了超分异花岗岩特征。

5.4.2 岩浆演化与锡多金属成矿作用

Zr和Hf、Nb和Ta及Y和Ho等元素的地球化学行为一致,被称为“孪生元素”,Zr/Hf、Nb/Ta及Y/Ho比值在一般的岩浆体系中并不发生数值的变化(Green, 1995),但当岩浆由于分异而发生性质上的明显改变时,这些比值都将显著变小(Bau, 1996; Linnen and Keppler, 2002; Deering and Bachmann, 2010; Dostal et al., 2015; Ballouard et al., 2016)。因而,Zr/Hf和Nb/Ta比值可视为花岗岩浆结晶分异程度的标志(Pérez-Soba and Villaseca, 2010; 陈璟元和杨进辉, 2015; Ballouard et al., 2016; 袁顺达, 2017; Yuan et al., 2018b)。如Bau(1996)就提出将Zr/Hf=26作为花岗岩体系岩浆-热液的分界,Ballouard et al. (2016)以Nb/Ta=5将花岗岩划分为正常结晶分异成因和岩浆-热液相互作用成因。吴福元等(2017)认为这些都可能是高分异花岗岩与超分异花岗岩的地球化学界限,只不过超分异花岗岩很可能与岩浆热液作用关系密切。按照上述标准,北大山岩体属于高分异花岗岩,而维拉斯托岩体属于超分异花岗岩(图 15a, b)。Irber (1999)以1.0和1.1的TE1, 3值为界限,认为TE1, 3值小于1.0的岩石不存在稀土元素四分组效应,TE1, 3值大于1.1的岩石显示了明显的稀土四分组效应。北大山岩体的TE1, 3值介于1.02~1.08之间,维拉斯托岩体的TE1, 3值介于1.48~1.61之间(表 2),指示维拉斯托岩体存在明显的稀土四分组效应,而北大山岩体仅有轻微的稀土四分组效应。从北大山岩体到维拉斯托岩体,岩石的稀土四分组效应明显增强(图 15c, d),Y/Ho比值降低(图 15d)。通常认为,岩浆发生高度分离结晶后,熔体中出现大量的流体,而熔体与流体的相互作用,将改变岩浆中稀土元素的地球化学行为,从而使它们呈现四分组效应(Bau, 1996; Irber, 1999; Monecke et al., 2002)。赵振华等(1992)系统讨论了中国华南稀有金属花岗岩的稀土元素四分组效应,认为花岗质熔体的高度分离结晶和与富含F和Cl的挥发分流体的相互作用是形成稀土四分组效应的重要控制因素。岩浆的分离结晶作用只能导致Y/Ho比值发生小的变化,通常在24~34之间,但熔体-流体的相互作用可以显著改变Y/Ho比值。在图 15d中,北大山岩体的Y/Ho比值基本落在24~34区域,但维拉斯托岩体的Y/Ho比值明显降低,表明北大山岩体仅经历了分离结晶作用,而维拉斯托岩体除了分离结晶,还发生了显著的熔体-流体相互作用。

图 15 花岗岩演化程度判别图 Fig. 15 Discrimination diagrams for evolution degree of granites

除了维拉斯托稀有金属-锡多金属矿床,大兴安岭南段还发育众多的锡多金属矿床。我们收集了该区典型锡多金属矿床成矿岩体地球化学数据,并计算了相关参数。白音查干东山超大型Sn-Ag-Pb-Zn-Cu-Sb矿床成矿岩体的Nb/Ta比值介于6.6~6.9之间,Zr/Hf比值介于20.6~22.0之间,TE1, 3值介于1.08~1.10之间(李睿华, 2019),岩体的分异程度较高,接近超分异花岗岩,显示了相对明显的稀土四分组效应。黄岗超大型Sn-Fe-W矿床成矿岩体的Nb/Ta比值介于4.6~9.1之间,Zr/Hf比值介于20.2~27.8之间,TE1, 3值介于1.02~1.13之间(Zhou et al., 2012),岩体的总体分异程度亦较高,显示了相对明显的稀土四分组效应。大井大型Sn-Zn-Cu-Pb-Ag矿床成矿岩体的Nb/Ta比值介于10.1~10.7之间,Zr/Hf比值介于31.6~32.3之间,TE1, 3值介于0.99~1.00之间(Liu et al., 2001),岩体的分异程度相对较低,不显示稀土四分组效应。道伦达坝中型Cu-W-Sn矿床成矿岩体的Nb/Ta比值介于8.3~11.7之间,Zr/Hf比值介于17.5~27.1之间,TE1, 3值介于1.19~1.37之间(陈公正, 2018),岩体的分异程度相对较高,显示了明显的稀土四分组效应。毛登-小孤山中型Sn-Zn-Cu矿床成矿花岗斑岩的Nb/Ta比值介于12.3~14.7之间,Zr/Hf比值介于29.1~31.8之间,TE1, 3值介于0.97~1.06之间(李睿华, 2019),岩体的分异程度较低,且稀土四分组效应不明显。东山湾小型Sn-W-Mo-Cu矿床成矿岩体的Nb/Ta比值介于6.8~10.7之间,Zr/Hf比值介于15.9~24.2之间,TE1, 3值介于0.86~1.12之间(Zeng et al., 2015; Zhang et al., 2017),岩体的分异程度相对较高,总体上,岩体的稀土四分组效应不明显。宝盖沟小型Sn矿床成矿岩体的Nb/Ta比值介于3.1~11.1之间,Zr/Hf比值介于27.8~30.3之间,TE1, 3值介于0.99~1.02之间(Mi et al., 2020),岩体的分异程度相对较低,稀土四分组效应不明显。查木罕小型W-Mo-Be矿床成矿岩体的Nb/Ta比值介于5.9~14.1之间,Zr/Hf比值介于13.5 ~27.5之间,TE1, 3值介于1.08~1.21之间(Zhang et al., 2019),总体上,岩体的分异程度较高,且显示了明显的稀土四分组效应。本文和前人的研究结果表明,大兴安岭南段锡多金属矿床的成矿岩体均为分异型花岗岩,岩浆演化程度和熔体-流体相互作用是导致大兴安岭南段锡多金属矿床不同矿种组合的控制因素。稀有金属-锡多金属矿床为超分异花岗岩,具有明显的稀土元素四分组效应;Sn-W组合矿床的岩体分异程度亦较高,显示较明显的稀土元素四分组效应;独立的Sn矿床岩体的分异程度明显低于稀有金属-锡多金属组合和Sn-W组合矿床,且通常不显示稀土四分组效应;Cu-Sn-Pb-Zn-Ag矿床的成矿岩体分异程度最低,不显示稀土四分组效应,如:大井矿床、毛登-小孤山矿床。

5.4.3 维拉斯托矿床成因类型和成矿模式

关于维拉斯托稀有金属-锡多金属矿床的成因类型,尚有不同认识,多数研究者认为维拉斯托稀有金属-锡多金属矿床属于斑岩-热液脉复合型成矿系统(刘翼飞等, 2014; Liu et al., 2016; 翟德高等, 2016; 樊志勇等, 2017; Wang et al., 2017);但亦有部分研究者认为属于岩浆-热液型锡多金属矿床(祝新友等, 2016)。矿床成因争议的根本原因是对矿区隐伏的花岗质岩株岩石类型的不同认识和对岩体中钠长石形成机制的理解。刘翼飞等(2014)Liu et al. (2016)翟德高等(2016)Wang et al. (2017)认为该岩株的岩石类型为石英斑岩,而祝新友等(2016)将该隐伏岩株厘定为细粒斑状碱长花岗岩;多数研究者认为岩体中的钠长石是热液成因的钠长石化所致(Liu et al., 2016; 刘瑞麟等, 2018b),部分研究者认为钠长石是岩浆结晶产物(祝新友等, 2016)。维拉斯托岩体顶部岩石发育似斑状结构,向下逐渐变为中粒花岗岩,因此维拉斯托岩体不是石英斑岩。岩体顶部斑状花岗岩的斑晶主要为石英、钾长石、钠长石,斑晶粒度差别极大;基质为细粒结构,主要矿物为钠长石、天河石、石英、黄玉、锂云母和萤石。石英斑晶边缘含大量细粒半自形钠长石、他形钾长石,它们沿石英结晶面生长,构成一圈或多圈环带,形成明显的“雪球结构”。“雪球结构”是黄玉-锂云母-钠长石花岗岩中一种特有结构,是岩浆作用所致,即在富F、Na、Al的花岗质岩浆结晶初期,由于高成核密度和低生长速率而形成了大量的钠长石晶核(也可能同时形成了少量石英晶核),随着岩浆的演化,体系逐渐转变为适合于钠长石、石英和钾长石同时生长,在H2O含量较高的条件下,石英的生长速率高于钠长石,相对快速生长的石英晶体在向外生长的过程中逐渐包容了钠长石晶体,最终形成了按斑晶晶面排列的“雪球结构”(Pollard, 1989; Kempe et al., 1999; 朱金初等, 2002)。Kovalenko and Kovalenko (1976)朱金初等(2002)认为稀有金属花岗岩的成矿作用发生在结晶分异的过程中,极高的分异程度造成高场强元素在残余熔体中高度富集,最终达到饱和并发生矿化。综上所述,维拉斯托岩体不是石英斑岩,而是超分异的Li-F花岗岩;钠长石主要形成于岩浆阶段;维拉斯托矿床不是斑岩成矿系统,而是与超分异花岗岩有关的岩浆-热液成矿系统。其成矿模式如图 16

图 16 维拉斯托稀有金属-锡多金属成矿系统的成矿模式 Fig. 16 Ore deposit model of the Weilasituo rare metal-tin-polymetallic metallogenic system

早白垩世期间,大兴安岭南段处于伸展环境,导致岩石圈减薄、软流圈上涌,幔源岩浆以底侵方式进入下地壳,触发了新元古代期间增生的来自亏损地幔的新生下地壳的熔融,形成了早白垩世I型花岗岩质岩浆。岩浆上侵过程中发生分异,形成了不同演化程度的杂岩体。花岗质熔体的高度分离结晶和与富含F和Cl的挥发分流体的相互作用可以使稀有金属元素在岩浆阶段富集成矿(赵振华等, 1992; Yuan et al., 2019)。其成矿过程为:富氟富稀有金属元素的岩浆先通过结晶分异作用,结晶出石英、钠长石、钾长石和云母等造岩矿物,导致残余熔体逐渐富集不相容元素和挥发分;当稀有金属成矿元素和挥发分富集到一定程度后,黄玉开始晶出,导致熔体中的氟含量明显降低,从而使稀有金属矿物(如铌铁矿、钽铁矿等)开始沉淀成矿(Cuney et al., 1992),形成维拉斯托岩体顶部的蚀变花岗岩型Sn-Zn(-Rb-Nb-Ta)矿体。随着体系中晶体相所占比例的逐渐增多,导致熔体中的水含量在一定时间达到饱和,进而独立的流体相出溶,进入岩浆-热液过渡阶段,该阶段因含有大量的挥发分,发生隐爆,在隐爆角砾岩筒中发育广泛云英岩化,形成以含锂白云母和锂云母为主的富云母云英岩,出溶流体中还含有较多的Rb和Be,导致维拉斯托隐爆角砾岩型Sn-Zn-Li(-Rb-Be)矿体的形成。当岩浆体系的温度下降至固相线时,熔体相消失,此时体系中只存在固体相和流体相,成矿系统进入热液成矿阶段,高温成矿流体运移到断裂-裂隙系统,形成受断裂构造控制的石英脉型锡多金属矿体;成矿流体的温度继续降低,演化到中温阶段,形成了少量的钼多金属矿体(图 16)。

6 结论

(1) 维拉斯托地区的北大山岩体和维拉斯托岩体的LA-ICP-MS锆石U-Pb年龄分别为140Ma和137~138Ma,两者为同一期岩浆活动不同演化阶段的产物。

(2) 北大山岩体和维拉斯托岩体都经历了钾长石、斜长石、黑云母、独居石、榍石、褐帘石、磷灰石、锆石和钛铁矿的分离结晶;除了岩浆分离结晶,维拉斯托岩体还发生了熔体-流体相互作用,表现出明显的稀土元素四分组效应。

(3) 花岗质岩浆的高程度分离结晶和熔体-流体相互作用是形成维拉斯托超分异花岗岩的重要控制因素,并发生稀有金属-锡多金属成矿作用。

(4) 北大山岩体和维拉斯托岩体的岩浆源区为起源于亏损地幔的新元古代新生地壳物质,岩体形成于伸展环境。

致谢      野外工作得到了内蒙古维拉斯托矿业有限公司各级领导的大力支持;内蒙古自治区地质调查院的章培春、李雪娇、赵静、闫洁、张婷婷和魏雅玲参加了部分野外工作;两位审稿人和终审主编认真审阅了本文并提出了详细的修改意见,这些意见和建议使本文整体质量得到很大的提升。在此一并致以诚挚的感谢。

参考文献
Andersen T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb. Chemical Geology, 192(1-2): 59-79 DOI:10.1016/S0009-2541(02)00195-X
Ballouard C, Poujol M, Boulvais P, Branquet Y, Tartèse R and Vigneresse JL. 2016. Nb-Ta fractionation in peraluminous granites: A marker of the magmatic-hydrothermal transition. Geology, 44(3): 231-234 DOI:10.1130/G37475.1
Batchelor RA and Bowden P. 1985. Petrogenetic interpretation of granitoid rock series using multicationic parameters. Chemical Geology, 48(1-4): 43-55 DOI:10.1016/0009-2541(85)90034-8
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
Belousova E, Griffin W, O'Reilly SY and Fisher N. 2002. Igneous zircon: Trace element composition as an indicator of source rock type. Contributions to Mineralogy and Petrology, 143(5): 602-622 DOI:10.1007/s00410-002-0364-7
Black LP and Gulson BL. 1978. The age of the mud tank carbonatite, strangways range, northern territory. BMR Journal of Australian Geology and Geophysics, 3(3): 227-232
Blichert-Toft J and Albarède F. 1997. The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle-crust system. Earth and Planetary Science Letters, 148(1-2): 243-258 DOI:10.1016/S0012-821X(97)00040-X
Boynton WV. 1984. Cosmochemistry of the rare earth elements: Meteorite studies. Developments in Geochemistry, 2: 63-114
Brown GC, Thorpe RS and Webb BC. 1984. The geochemical characteristics of granitoids in contrasting arcs and comments on magma sources. Journal of the Geological Society, 141(3): 413-426 DOI:10.1144/gsjgs.141.3.0413
Bureau of Geology and Mineral Resources of Inner Mongolian Autonomous Region. 1991. Regional Geology of Inner Mongolian Autonomous Region. Beijing: Geological Publishing House, 1-498 (in Chinese)
Chen B, Jahn BM and Tian W. 2009. Evolution of the Solonker suture zone: Constraints from zircon U-Pb ages, Hf isotopic ratios and whole-rock Nd-Sr isotope compositions of subduction- and collision-related magmas and forearc sediments. Journal of Asian Earth Sciences, 34(3): 245-257 DOI:10.1016/j.jseaes.2008.05.007
Chen GZ. 2018. Metallogenesis of the Daolundaba Cu-polymetallic deposit in the southern Graet Xing'an Range. Master Degree Thesis. Beijing: China University of Geosciences, 1-131 (in Chinese with English summary)
Chen JY and Yang JH. 2015. Petrogenesis of the Fogang highly fractionated Ⅰ-type granitoids: Constraints from Nb, Ta, Zr and Hf. Acta Petrologica Sinica, 31(3): 846-854 (in Chinese with English abstract)
Chen YJ, Zhang C, Wang P, Pirajno F and Li N. 2017. The Mo deposits of Northeast China: A powerful indicator of tectonic settings and associated evolutionary trends. Ore Geology Reviews, 81: 602-640 DOI:10.1016/j.oregeorev.2016.04.017
Cuney M, Marignac C and Weisbrod A. 1992. The Beauvoir topaz-lepidolite albite granite (Massif Central, France): The disseminated magmatic Sn-Li-Ta-Nb-Be mineralization. Economic Geology, 87(7): 1766-1794 DOI:10.2113/gsecongeo.87.7.1766
Davis GA, Cong W, Zheng YD, Zhang JJ, Zhang CH and Gehrels GE. 1998. The enigmatic Yinshan fold-and-thrust belt of northern China: New views on its intraplate contractional styles. Geology, 26(1): 43-46 DOI:10.1130/0091-7613(1998)026<0043:TEYFAT>2.3.CO;2
De La Roche H, Leterrier J, Grandclaude P and Marchal M. 1980. A classification of volcanic and plutonic rocks using R1R2-diagram and major-element analyses: Its relationships with current nomenclature. Chemical Geology, 29(1-4): 183-210 DOI:10.1016/0009-2541(80)90020-0
Deering CD and Bachmann O. 2010. Trace element indicators of crystal accumulation in silicic igneous rocks. Earth and Planetary Science Letters, 297(1-2): 324-331 DOI:10.1016/j.epsl.2010.06.034
Dostal J, Kontak DJ, Gerel O, Gregory SJ and Fayek M. 2015. Cretaceous ongonites (topaz-bearing albite-rich microleucogranites) from Ongon Khairkhan, Central Mongolia: Products of extreme magmatic fractionation and pervasive metasomatic fluid: Rock interaction. Lithos, 236-237: 173-189 DOI:10.1016/j.lithos.2015.08.003
Fan ZY, Qiu HY, Fu X, Wang KX and Hugejiletu. 2017. Discovery and exploration of Weilasituo large porphyry-type tin-polymetal deposit in Inner Mongolia and its geological significances. Gold Science and Technology, 25(1): 9-17 (in Chinese with English abstract)
Fujimaki H. 1986. Partition coefficients of Hf, Zr, and REE between zircon, apatite, and liquid. Contributions to Mineralogy and Petrology, 94(1): 42-45 DOI:10.1007/BF00371224
Ge MC, Zhou WX, Yu Y, Sun JJ, Bao JQ and Wang SH. 2011. Dissolution and supracrustal rocks dating of Xilin Gol complex, Inner Mongolia, China. Earth Science Frontiers, 18(5): 182-195 (in Chinese with English abstract)
Ge WC, Wu FY, Zhou CY and Zhang JH. 2007. Porphyry Cu-Mo deposits in the eastern Xing'an-Mongolian Orogenic Belt: Mineralization ages and their geodynamic implications. Chinese Science Bulletin, 52(24): 3416-3427 DOI:10.1007/s11434-007-0466-8
Gelman SE, Deering CD, Bachmann O, Huber C and Gutiérrez FJ. 2014. Identifying the crystal graveyards remaining after large silicic eruptions. Earth and Planetary Science Letters, 403: 299-306 DOI:10.1016/j.epsl.2014.07.005
Gervasoni F, Klemme S, Rocha-Júnior ERV and Berndt J. 2016. Zircon saturation in silicate melts: A new and improved model for aluminous and alkaline melts. Contributions to Mineralogy and Petrology, 171(3): 21 DOI:10.1007/s00410-016-1227-y
Green TH. 1995. Significance of Nb/Ta as an indicator of geochemical processes in the crust-mantle system. Chemical Geology, 120(3-4): 347-359 DOI:10.1016/0009-2541(94)00145-X
Griffin WL, Pearson NJ, Belousova E, Jackson SE, Van Achterbergh E, O'Reilly SY and Shee SR. 2000. The Hf isotope composition of cratonic mantle: LAM-MC-ICPMS analysis of zircon megacrysts in kimberlites. Geochimica et Cosmochimica Acta, 64(1): 133-147 DOI:10.1016/S0016-7037(99)00343-9
Griffin WL, Belousova EA, Shee SR, Pearson NJ and O'Reilly SY. 2004. Archean crustal evolution in the northern Yilgarn Craton: U-Pb and Hf-isotope evidence from detrital zircons. Precambrian Research, 131(3-4): 231-282 DOI:10.1016/j.precamres.2003.12.011
Guo GJ. 2016. Discussion on geological characteristics and metallogenic origin of Weilasituo Sn polymetal deposit in Inner Mongolia. Master Degree Thesis. Beijing: China University of Geosciences, 1-78 (in Chinese with English summary)
Hou KJ, Li YH and Tian YR. 2009. In situ U-Pb zircon dating using laser ablation-multi ion counting-ICP-MS. Mineral Deposits, 28(4): 481-492 (in Chinese with English abstract)
Irber W. 1999. The lanthanide tetrad effect and its correlation with K/Rb, Eu/Eu*, Sr/Eu, Y/Ho, and Zr/Hf of evolving peraluminous granite suites. Geochimica et Cosmochimica Acta, 63(3-4): 489-508 DOI:10.1016/S0016-7037(99)00027-7
Jackson SE, Pearson NJ, Griffin WL and Belousova EA. 2004. The application of laser ablation-inductively coupled plasma-mass spectrometry to in situ U-Pb zircon geochronology. Chemical Geology, 211(1-2): 47-69 DOI:10.1016/j.chemgeo.2004.06.017
Kempe U, Götze J, Dandar S and Habermann D. 1999. Magmatic and metasomatic processes during formation of the Nb-Zr-REE deposits Khaldzan Buregte and Tsakhir (Mongolian Altai): Indications from a combined CL-SEM study. Mineralogical Magazine, 63(2): 165-177 DOI:10.1180/002646199548402
King PL, White AJR, Chappell BW and Allen CM. 1997. Characterization and origin of aluminous A-type granites from the Lachlan Fold Belt, Southeastern Australia. Journal of Petrology, 38(3): 371-391 DOI:10.1093/petroj/38.3.371
Kovalenko VI and Kovalenko NI. 1976. Ongonites (Topaz-Bearing Quartz Keratophyre) Subvolcanic Analogues of Rare Metal Li-F Granites. Moscow: Nauka Press, 1-124
Li RH. 2019. Metallogenesis of tin-copper-polymetallic ore deposits in the Xilinhot area, the southern Great Xing'an Range, China. Ph. D. Dissertation. Beijing: Peking University, 1-239 (in Chinese with English summary)
Li XH, Li WX and Li ZX. 2007a. On the genetic classification and tectonic implications of the Early Yanshanian granitoids in the Nanling Range, South China. Chinese Science Bulletin, 52(14): 1873-1885 DOI:10.1007/s11434-007-0259-0
Li XH, Li ZX, Li WX, Liu Y, Yuan C, Wei GJ and Qi CS. 2007b. U-Pb zircon, geochemical and Sr-Nd-Hf isotopic constraints on age and origin of Jurassic I- and A-type granites from central Guangdong, SE China: A major igneous event in response to foundering of a subducted flat-slab?. Lithos, 96(1-2): 186-204 DOI:10.1016/j.lithos.2006.09.018
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-3): 213-227 DOI:10.1007/s004100050304
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 RL, Wu G, Li TG, Chen GZ, Wu LW, Zhang PC, Zhang T, Jiang B and Liu WY. 2018a. LA-ICP-MS cassiterite and zircon U-Pb ages of the Weilasituo tin-polymetallic deposit in the southern Great Xing'an Range and their geological significance. Earth Science Frontiers, 25(5): 183-201 (in Chinese with English abstract)
Liu RL, Wu G, Chen GZ, Li TG, Jiang B, Wu LW, Zhang PC, Zhang T and Chen YC. 2018b. Characteristics of fluid inclusions and H-O-C-S-Pb isotopes of Weilasituo Sn-polymetallic deposit in southern Da Hinggan Mountains. Mineral Deposits, 37(2): 199-224 (in Chinese with English abstract)
Liu W, Li XJ and Tan J. 2001. Petrogenetic and metallogenetic background of the Dajing Cu-Sn-Ag-Pb-Zn ore deposit, Inner Mongolia, and characteristics of the mineralizing fluid. Resource Geology, 51(4): 321-331 DOI:10.1111/j.1751-3928.2001.tb00105.x
Liu W, Siebel W, Li XJ and Pan XF. 2005. Petrogenesis of the Linxi granitoids, northern Inner Mongolia of China: Constraints on basaltic underplating. Chemical Geology, 219(1-4): 5-35 DOI:10.1016/j.chemgeo.2005.01.013
Liu YF. 2009. Metallogenic study of Bairendaba Ag polymetallic deposit in Hexigten Banner, Inner Mongolia. Master Degree Thesis. Beijing: Chinese Academy of Geological Sciences, 1-98 (in Chinese with English summary)
Liu YF, Fan ZY, Jiang HC, Nie FJ, Jiang SH, Ding CW and Wang FX. 2014. Genesis of the Weilasituo-Bairendaba porphyry-hydrothermal vein type system in Inner Mongolia, China. Acta Geologica Sinica, 88(12): 2373-2385 (in Chinese with English abstract)
Liu YF, Jiang SH and Bagas L. 2016. The genesis of metal zonation in the Weilasituo and Bairendaba Ag-Zn-Pb-Cu-(Sn-W) deposits in the shallow part of a porphyry Sn-W-Rb system, Inner Mongolia, China. Ore Geology Reviews, 75: 150-173 DOI:10.1016/j.oregeorev.2015.12.006
Ludwig KR. 2003. User's Manual for Isoplot 3.0: A Geochronological Toolkit for Microsoft Excel. Berkeley: Berkeley Geochronology Center Special Publication, 1-71
Mahood G and Hildreth W. 1983. Large partition coefficients for trace elements in high-silica rhyolites. Geochimica et Cosmochimica Acta, 47(1): 11-30 DOI:10.1016/0016-7037(83)90087-X
Maniar PD and Piccoli PM. 1989. Tectonic discrimination of granitoids. GSA Bulletin, 101(5): 635-643 DOI:10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
Mao JW, Xie GQ, Zhang ZH, Li XF, Wang YT, Zhang CQ and Li YF. 2005. Mesozoic large-scale metallogenic pulses in North China and corresponding geodynamic settings. Acta Petrologica Sinica, 21(1): 169-188 (in Chinese with English abstract)
Mao JW, Cheng YB, Chen MH and Pirajno F. 2013. Major types and time-space distribution of Mesozoic ore deposits in South China and their geodynamic settings. Mineralium Deposita, 48(3): 267-294 DOI:10.1007/s00126-012-0446-z
Mao JW, Zhou ZH, Wu G, Jiang SH, Liu CL, Li HM, Ouyang HG and Liu J. 2013. Metallogenic regularity and minerogenetic series of ore deposits in Inner Mongolia and adjacent areas. Mineral Deposits, 32(4): 715-729 (in Chinese with English abstract)
Mao JW, Ouyang HG, Song SW, Santosh M, Yuan SD, Zhou ZH, Zheng W, Liu H, Liu P, Cheng YB and Chen MH. 2019. Geology and metallogeny of tungsten and tin deposits in China. Society of Economic Geologists Special Publication, 22: 411-482
Mao JW, Zheng W, Xie GQ, Lehmann B and Goldfarb R. 2021. Recognition of a Middle-Late Jurassic arc-related porphyry copper belt along the Southeast China coast: Geological characteristics and metallogenic implications. Geology: 49 DOI:10.1130/G48615.1
Meng QR. 2003. What drove Late Mesozoic extension of the northern China-Mongolia tract?. Tectonophysics, 369(3-4): 155-174 DOI:10.1016/S0040-1951(03)00195-1
Mi KF, Lü ZC, Yan TJ, Yao XF, Ma YX and Lin CG. 2020. Zircon geochronological and geochemical study of the Baogaigou tin deposits, southern Great Xing'an Range, Northeast China: Implications for the timing of mineralization and ore genesis. Geological Journal, 55(7): 5062-5081 DOI:10.1002/gj.3729
Monecke T, Kempe U, Monecke J, Sala M and Wolf D. 2002. Tetrad effect in rare earth element distribution patterns: A method of quantification with application to rock and mineral samples from granite-related rare metal deposits. Geochimica et Cosmochimica Acta, 66(7): 1185-1196 DOI:10.1016/S0016-7037(01)00849-3
Münker C, Pfänder JA, Weyer S, Büchl, Kleine T and Mezger K. 2003. Evolution of planetary cores and the earth-moon system from Nb/Ta systematics. Science, 301(5629): 84-87 DOI:10.1126/science.1084662
Nowell GM, Kempton PD, Noble SR, Fitton JG, Saunders AD, Mahoney JJ and Taylor RN. 1998. High precision Hf isotope measurements of MORB and OIB by thermal ionisation mass spectrometry: Insights into the depleted mantle. Chemical Geology, 149(3-4): 211-233 DOI:10.1016/S0009-2541(98)00036-9
Ouyang HG, Mao JW, Zhou ZH and Su HM. 2015. Late Mesozoic metallogeny and intracontinental magmatism, southern Great Xing'an Range, northeastern China. Gondwana Research, 27(3): 1153-1172 DOI:10.1016/j.gr.2014.08.010
Pearce JA, Harris NBW and Tindle AG. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. Journal of Petrology, 25(4): 956-983 DOI:10.1093/petrology/25.4.956
Peccerillo A and Taylor SR. 1976. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey. Contributions to Mineralogy and Petrology, 58(1): 63-81 DOI:10.1007/BF00384745
Pérez-Soba C and Villaseca C. 2010. Petrogenesis of highly fractionated Ⅰ-type peraluminous granites: La Pedriza pluton (Spanish Central System). Geologica Acta, 8(2): 131-149
Pfänder JA, Münker C, Stracke A and Mezger K. 2007. Nb/Ta and Zr/Hf in ocean island basalts-Implications for crust-mantle differentiation and the fate of niobium. Earth and Planetary Science Letters, 254(1-2): 158-172 DOI:10.1016/j.epsl.2006.11.027
Pirajno F, Ernst RE, Borisenko AS, Fedoseev G and Naumov EA. 2009. Intraplate magmatism in Central Asia and China and associated metallogeny. Ore Geology Reviews, 35(2): 114-136 DOI:10.1016/j.oregeorev.2008.10.003
Pollard PJ. 1989. Geochenmistry of Granite Associated with Tantalum and Niobium Mineraliazation Lanthanides, Tantalum and Niobium. Berlin: Springer, 145-168
Pupin JP. 1980. Zircon and granite petrology. Contributions to Mineralogy and Petrology, 73(3): 207-220 DOI:10.1007/BF00381441
Rollinson HR. 1993. Using Geochemical Data: Evaluation, Presentation, Interpretation. Harlow: Longman, 1-325
Rudnick RL and Gao S. 2004. Composition of the continental crust. Treatise on Geochemistry, 4: 1-51
Söderlund U, Patchett PJ, Vervoort JD and Isachsen CE. 2004. The 176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions. Earth and Planetary Science Letters, 219(3-4): 311-324 DOI:10.1016/S0012-821X(04)00012-3
Stepanov AS, Mavrogenes JA, Meffre S and Davidson P. 2014. The key role of mica during igneous concentration of tantalum. Contributions to Mineralogy and Petrology, 167(6): 1009 DOI:10.1007/s00410-014-1009-3
Streckeisen A. 1976. To each plutonic rock its proper name. Earth-Science Reviews, 12(1): 1-33 DOI:10.1016/0012-8252(76)90052-0
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 Publication, 42(1): 313-345
Tiepolo M, Vannucci R, Oberti R, Foley S and Zanetti A. 2000. Nb and Ta incorporation and fractionation in titanian pargasite and kaersutite: Crystal-chemical constraints and implications for natural systems. Earth and Planetary Science Letters, 176(2): 185-201 DOI:10.1016/S0012-821X(00)00004-2
Vervoort JD and Blichert-Toft J. 1999. Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochimica et Cosmochimica Acta, 63(3-4): 533-556 DOI:10.1016/S0016-7037(98)00274-9
Wang CG, Sun FY, Sun GS, Sun JD, Li Y and Feng HD. 2016a. Geochronology, geochemical and isotopic constraints on petrogenesis of intrusive complex associated with Bianjiadayuan polymetallic deposit on the southern margin of the Greater Khingan, China. Arabian Journal of Geosciences, 9(5): 334 DOI:10.1007/s12517-016-2333-x
Wang FX, Bagas L, Jiang SH and Liu YF. 2017. Geological, geochemical, and geochronological characteristics of Weilasituo Sn-polymetal deposit, Inner Mongolia, China. Ore Geology Reviews, 80: 1206-1229 DOI:10.1016/j.oregeorev.2016.09.021
Wang J. 2009. Chronology and geochemistry of granitoid for the Weilasituo copper polymetal deposit in Inner Mongolia. Master Degree Thesis. Beijing: China University of Geosciences (Beijing), 1-68 (in Chinese with English summary)
Wang JB, Wang YW, Wang LJ and Uemoto T. 2001. Tin-polymetallic mineralization in the southern part of the Da Hinggan Mountains, China. Resource Geology, 51(4): 283-291 DOI:10.1111/j.1751-3928.2001.tb00102.x
Wang T, Guo L, Zheng YD, Donskaya T, Gladkochub D, Zeng LS, Li JB, Wang YB and Mazukabzov A. 2012. Timing and processes of Late Mesozoic mid-lower-crustal extension in continental NE Asia and implications for the tectonic setting of the destruction of the North China Craton: Mainly constrained by zircon U-Pb ages from metamorphic core complexes. Lithos, 154: 315-345 DOI:10.1016/j.lithos.2012.07.020
Wang X, Chen J and Ren MH. 2016b. Hydrothermal zircon geochronology: Age constraint on Nanling Range tungsten mineralization (Southeast China). Ore Geology Reviews, 74: 63-75 DOI:10.1016/j.oregeorev.2015.10.034
Wang XY, Hou QY, Wang J, Chen YL, Liu JB, Wang Z and Li DP. 2013. SHRIMP geochronology and Hf isotope of zircons from granitoids of the Weilasituo deposit in Inner Mongolia. Geoscience, 27(1): 67-78 (in Chinese with English abstract)
Watson EB and Harrison TM. 1983. Zircon saturation revisited: Temperature and composition effects in a variety of crustal magma types. Earth and Planetary Science Letters, 64(2): 295-304 DOI:10.1016/0012-821X(83)90211-X
Watson EB and Harrison TM. 1984. Accessory minerals and the geochemical evolution of crustal magmatic systems: A summary and prospectus of experimental approaches. Physics of the Earth and Planetary Interiors, 35(1-3): 19-30 DOI:10.1016/0031-9201(84)90031-1
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
Wu FY, Sun DY, Li H, Jahn BM and Wilde S. 2002. A-type granites in northeastern China: Age and geochemical constraints on their petrogenesis. Chemical Geology, 187(1-2): 143-173 DOI:10.1016/S0009-2541(02)00018-9
Wu FY, Li XH, Yang JH and Zheng YF. 2007. Discussions on the petrogenesis of granites. Acta Petrologica Sinica, 23(6): 1217-1238 (in Chinese with English abstract)
Wu FY, Sun DY, Ge WC, Zhang YB, Grant ML, Wilde SA and Jahn BM. 2011. Geochronology of the Phanerozoic granitoids in Northeastern China. Journal of Asian Earth Sciences, 41(1): 1-30 DOI:10.1016/j.jseaes.2010.11.014
Wu FY, Liu XC, Ji WQ, Wang JM and Yang L. 2017. Highly fractionated granites: Recognition and research. Scientia Sinica (Terrae), 47(7): 745-765 (in Chinese) DOI:10.1360/N072016-00139
Wu G, Li XZ, Xu LQ, Wang GR, Liu J, Zhang T, Quan ZX, Wu H, Li TG, Zeng QT and Chen YC. 2017. Age, geochemistry, and Sr-Nd-Hf-Pb isotopes of the Caosiyao porphyry Mo deposit in Inner Mongolia, China. Ore Geology Reviews, 81: 706-727 DOI:10.1016/j.oregeorev.2016.06.002
Xiao WJ, Zhang LC, Qin KZ, Sun S and Li JL. 2004. Paleozoic accretionary and collisional tectonics of the eastern Tianshan (China): Implications for the continental growth of central Asia. American Journal of Science, 304(4): 370-395 DOI:10.2475/ajs.304.4.370
Xu B, Jiang SY, Wang R, Ma L, Zhao KD and Yan X. 2015. Late Cretaceous granites from the giant Dulong Sn-polymetallic ore district in Yunnan Province, South China: Geochronology, geochemistry, mineral chemistry and Nd-Hf isotopic compositions. Lithos, 218-219: 54-72 DOI:10.1016/j.lithos.2015.01.004
Xu WL, Pei FP, Wang F, Meng E, Ji WQ, Yang DB and Wang W. 2013. Spatial-temporal relationships of Mesozoic volcanic rocks in NE China: Constraints on tectonic overprinting and transformations between multiple tectonic regimes. Journal of Asian Earth Sciences, 74: 167-193 DOI:10.1016/j.jseaes.2013.04.003
Xu ZG, Chen YC, Wang DH, Chen ZH and Li HM. 2008. Scheme of the Classification of the Minerogenetic Units in China. Beijing: Geological Publishing House, 1-103 (in Chinese)
Xue HM, Guo LJ, Hou ZQ, Zhou XW, Tong Y and Pan XF. 2009. The Xilingele complex from the eastern part of the Central Asian-Mongolia Orogenic Belt, China: Products of Early Variscan orogeny other than ancient block: Evidence from zircon SHRIMP U-Pb ages. Acta Petrologica Sinica, 25(8): 2001-2010 (in Chinese with English abstract)
Xue HM, Guo LJ, Hou ZQ, Tong Y, Pan XF and Zhou XW. 2010. SHRIMP zircon U-Pb ages of the Middle Neopaleozoic unmetamorphosed magmatic rocks in the southwestern slope of the Da Hinggan Mountains, Inner Mongolia. Acta Petrologica et Mineralogica, 29(6): 811-823 (in Chinese with English abstract)
Yang JH, Wu FY, Shao JA, Wilde SA, Xie LW and Liu XM. 2006. Constraints on the timing of uplift of the Yanshan Fold and Thrust Belt, North China. Earth and Planetary Science Letters, 246(3-4): 336-352 DOI:10.1016/j.epsl.2006.04.029
Yang W and Li SG. 2008. Geochronology and geochemistry of the Mesozoic volcanic rocks in western Liaoning: Implications for lithospheric thinning of the North China Craton. Lithos, 102(1-2): 88-117 DOI:10.1016/j.lithos.2007.09.018
Ying JF, Zhou XH, Zhang LC and Wang F. 2010. Geochronological framework of Mesozoic volcanic rocks in the Great Xing'an Range, NE China, and their geodynamic implications. Journal of Asian Earth Sciences, 39(6): 786-793 DOI:10.1016/j.jseaes.2010.04.035
Yu JJ, Wang F, Xu WL, Gao FH and Pei FP. 2012. Early Jurassic mafic magmatism in the Lesser Xing'an-Zhangguangcai Range, NE China, and its tectonic implications: Constraints from zircon U-Pb chronology and geochemistry. Lithos, 142-143: 256-266 DOI:10.1016/j.lithos.2012.03.016
Yu Y, Ge MC, Zhou WX, Sun JJ and Liu ZR. 2012. Petrology and metamorphic temperature-pressure conditions of Xilinhot Group, Inner Mongolia, China. Earth Science Frontiers, 19(5): 136-143 (in Chinese with English abstract)
Yuan HL, Gao S, Liu XM, Li HM, Günther D and Wu FY. 2004. Accurate U-Pb age and trace element determinations of zircon by laser ablation-inductively coupled plasma-mass spectrometry. Geostandards and Geoanalytical Research, 28(3): 353-370 DOI:10.1111/j.1751-908X.2004.tb00755.x
Yuan SD. 2017. Several crucial scientific issues related to the W-Sn metallogenesis in the Nanling Range and their implications for regional exploration: A review. Bulletin of Mineralogy, Petrology and Geochemistry, 36(5): 736-749 (in Chinese with English abstract)
Yuan SD, Mao JW, Zhao PL and Yuan YB. 2018a. Geochronology and petrogenesis of the Qibaoshan Cu-polymetallic deposit, northeastern Hunan Province: Implications for the metal source and metallogenic evolution of the intracontinental Qinhang Cu-polymetallic belt, South China. Lithos, 302-303: 519-534 DOI:10.1016/j.lithos.2018.01.017
Yuan SD, Williams-Jones AE, Mao JW, Zhao PL, Yan C and Zhang DL. 2018b. The origin of the Zhangjialong tungsten deposit, south China: Implications for W-Sn mineralization in large granite batholiths. Economic Geology, 113(5): 1193-1208 DOI:10.5382/econgeo.2018.4587
Yuan SD, Williams-Jones AE, Romer RL, Zhao PL and Mao JW. 2019. Protolith-related thermal controls on the decoupling of Sn and W in Sn-W metallogenic provinces: Insights from the Nanling region, China. Economic Geology, 114(5): 1005-1012 DOI:10.5382/econgeo.4669
Yurimoto H, Duke EF, Papike JJ and Shearer CK. 1990. Are discontinuous chondrite-normalized REE patterns in pegmatitic granite systems the results of monazite fractionation?. Geochimica et Cosmochimica Acta, 54(7): 2141-2145 DOI:10.1016/0016-7037(90)90277-R
Zeng QD, Sun Y, Chu SX, Duan XX and Liu JM. 2015. Geochemistry and geochronology of the Dongshanwan porphyry Mo-W deposit, Northeast China: Implications for the Late Jurassic tectonic setting. Journal of Asian Earth Sciences, 97: 472-485 DOI:10.1016/j.jseaes.2014.07.027
Zhai DG, Liu JJ, Zhang HY, Yao MJ, Wang JP and Yang YQ. 2014. S-Pb isotopic geochemistry, U-Pb and Re-Os geochronology of the Huanggangliang Fe-Sn deposit, Inner Mongolia, NE China. Ore Geology Reviews, 59: 109-122 DOI:10.1016/j.oregeorev.2013.12.005
Zhai DG, Liu JJ, Li MJ, Zhang M, Li BY, Fu X, Jiang HC, Ma LJ and Qi L. 2016. Geochronological study of Weilasituo porphyry type Sn deposit in Inner Mongolia and its geological significance. Mineral Deposits, 35(5): 1011-1022 (in Chinese with English abstract)
Zhai DG, Liu JJ, Zhang AL and Sun YQ. 2017. U-Pb, Re-Os, and 40Ar/39Ar geochronology of porphyry Sn ±Cu ±Mo and polymetallic (Ag-Pb-Zn-Cu) vein mineralization at Bianjiadayuan, Inner Mongolia, northeast China: Implications for discrete mineralization events. Economic Geology, 112(8): 2041-2059 DOI:10.5382/econgeo.2017.4540
Zhai MG, Meng QR, Liu JM, Hou QL, Hu SB, Li Z, Zhang HF, Liu W, Shao J and Zhu RX. 2004. Geological features of Mesozoic tectonic regime inversion in eastern North China and implication for geodynamics. Earth Science Frontiers, 11(3): 285-297 (in Chinese with English abstract)
Zhang AC, Wang RC, Li YL, Hu H, Lu XC, Ji JF and Zhang H. 2008. Tourmalines from the Koktokay No.3 pegmatite, Altai, NW China: Spectroscopic characterization and relationships with the pegmatite evolution. European Journal of Mineralogy, 20(1): 143-154 DOI:10.1127/0935-1221/2008/0020-1779
Zhang LL, Jiang SH, Bagas L, Han N, Liu Y and Liu YF. 2019. Element behaviour during interaction of magma and fluid: A case study of Chamuhan granite, and implications on the genesis of W-Mo mineralisation. Lithos, 342-343: 31-44 DOI:10.1016/j.lithos.2019.05.025
Zhang XB, Wang KY, Wang CY, Li W, Yu Q, Wang YC, Li JF, Wan D and Huang GH. 2017. Age, genesis, and tectonic setting of the Mo-W mineralized Dongshanwan granite porphyry from the Xilamulun metallogenic belt, NE China. Journal of Earth Science, 28(3): 433-446 DOI:10.1007/s12583-016-0934-1
Zhao ZH, Masuda A and Shabani MB. 1992. Tetrad effects of rare-earth elements in rare-metal granites. Geochimica, 28(3): 221-233 (in Chinese with English abstract)
Zhou ZH, Mao JW and Lyckberg P. 2012. Geochronology and isotopic geochemistry of the A-type granites from the Huanggang Sn-Fe deposit, southern Great Hinggan Range, NE China: Implication for their origin and tectonic setting. Journal of Asian Earth Sciences, 49: 272-286 DOI:10.1016/j.jseaes.2012.01.015
Zhu JC, Rao B, Xiong XL, Li FC and Zhang PH. 2002. Comparison and genetic interpretation of Li-F rich, rare-metal bearing granitic rocks. Geochimica, 31(2): 141-152 (in Chinese with English abstract)
Zhu XY, Zhang ZH, Fu X, Li BY, Wang YL, Jiao ST and Sun YL. 2016. Geological and geochemical characteristics of the Weilasito Sn-Zn deposit, Inner Mongolia. Geology in China, 43(1): 188-208 (in Chinese with English abstract)
陈公正. 2018. 大兴安岭南段道伦达坝铜多金属矿床成矿作用研究. 硕士学位论文. 北京: 中国地质大学(北京), 1-131
陈璟元, 杨进辉. 2015. 佛冈高分异I型花岗岩的成因: 来自Nb-Ta-Zr-Hf等元素的制约. 岩石学报, 31(3): 846-854.
樊志勇, 邱慧远, 付旭, 王可祥, 胡格吉乐吐. 2017. 内蒙古维拉斯托大型斑岩型锡多金属找矿勘查及启示. 黄金科学技术, 25(1): 9-17.
葛梦春, 周文孝, 于洋, 孙俊俊, 鲍建泉, 王世海. 2011. 内蒙古锡林郭勒杂岩解体及表壳岩系年代确定. 地学前缘, 18(5): 182-195.
郭贵娟. 2016. 内蒙古维拉斯托锡多金属矿床地质特征及成因探讨. 硕士学位论文. 北京: 中国地质大学(北京), 1-78
侯可军, 李延河, 田有荣. 2009. LA-MC-ICP-MS锆石微区原位U-Pb定年技术. 矿床地质, 28(4): 481-492. DOI:10.3969/j.issn.0258-7106.2009.04.010
李睿华. 2019. 大兴安岭南段锡林浩特地区锡铜多金属矿床的成矿作用. 博士学位论文. 北京: 北京大学, 1-239
刘瑞麟, 武广, 李铁刚, 陈公正, 武利文, 章培春, 张彤, 江彪, 刘文元. 2018a. 大兴安岭南段维拉斯托锡多金属矿床LA-ICP-MS锡石和锆石U-Pb年龄及其地质意义. 地学前缘, 25(5): 183-201.
刘瑞麟, 武广, 陈公正, 李铁刚, 江彪, 武利文, 章培春, 张彤, 陈毓川. 2018b. 大兴安岭南段维拉斯托锡多金属矿床流体包裹体和同位素特征. 矿床地质, 37(2): 199-224.
刘翼飞. 2009. 内蒙古克什克腾旗拜仁达坝银多金属矿床成因研究. 硕士学位论文. 北京: 中国地质科学院, 1-98
刘翼飞, 樊志勇, 蒋胡灿, 聂凤军, 江思宏, 丁成武, 王丰翔. 2014. 内蒙古维拉斯托-拜仁达坝斑岩-热液脉状成矿体系研究. 地质学报, 88(12): 2373-2385.
毛景文, 谢桂青, 张作衡, 李晓峰, 王义天, 张长青, 李永峰. 2005. 中国北方中生代大规模成矿作用的期次及其地球动力学背景. 岩石学报, 21(1): 169-188.
毛景文, 周振华, 武广, 江思宏, 刘成林, 李厚民, 欧阳荷根, 刘军. 2013. 内蒙古及邻区矿床成矿规律与成矿系列. 矿床地质, 32(4): 715-729. DOI:10.3969/j.issn.0258-7106.2013.04.006
内蒙古自治区地质矿产局. 1991. 内蒙古自治区区域地质志. 北京: 地质出版社, 1-498.
王瑾. 2009. 内蒙古维拉斯托铜多金属矿床矿区花岗岩类年代学与地球化学. 硕士学位论文. 北京: 中国地质大学(北京), 1-68
王新宇, 侯青叶, 王瑾, 陈岳龙, 刘金宝, 王忠, 李大鹏. 2013. 内蒙古维拉斯托矿床花岗岩类SHRIMP年代学及Hf同位素研究. 现代地质, 27(1): 67-78. DOI:10.3969/j.issn.1000-8527.2013.01.007
吴福元, 李献华, 杨进辉, 郑永飞. 2007. 花岗岩成因研究的若干问题. 岩石学报, 23(6): 1217-1238. DOI:10.3969/j.issn.1000-0569.2007.06.001
吴福元, 刘小驰, 纪伟强, 王佳敏, 杨雷. 2017. 高分异花岗岩的识别与研究. 中国科学(地球科学), 47(7): 745-765.
徐志刚, 陈毓川, 王登红, 陈郑辉, 李厚民. 2008. 中国成矿区带划分方案. 北京: 地质出版社, 1-103.
薛怀民, 郭利军, 侯增谦, 周喜文, 童英, 潘晓菲. 2009. 中亚-蒙古造山带东段的锡林郭勒杂岩: 早华力西期造山作用的产物而非古老陆块?锆石SHRIMP U-Pb年代学证据. 岩石学报, 25(8): 2001-2010.
薛怀民, 郭利军, 侯增谦, 童英, 潘晓菲, 周喜文. 2010. 大兴安岭西南坡成矿带晚古生代中期未变质岩浆岩的SHRIMP锆石U-Pb年代学. 岩石矿物学杂志, 29(6): 811-823. DOI:10.3969/j.issn.1000-6524.2010.06.016
于洋, 葛梦春, 周文孝, 孙俊俊, 刘泽瑞. 2012. 内蒙古锡林浩特岩群岩石学特征及变质温压条件. 地学前缘, 19(5): 136-143.
袁顺达. 2017. 南岭钨锡成矿作用几个关键科学问题及其对区域找矿勘查的启示. 矿物岩石地球化学通报, 36(5): 736-749. DOI:10.3969/j.issn.1007-2802.2017.05.004
翟德高, 刘家军, 李俊明, 张梅, 李泊洋, 付旭, 蒋胡灿, 马立军, 漆亮. 2016. 内蒙古维拉斯托斑岩型锡矿床成岩、成矿时代及其地质意义. 矿床地质, 35(5): 1011-1022.
翟明国, 孟庆任, 刘建明, 侯泉林, 胡圣标, 李忠, 张宏福, 刘伟, 邵济安, 朱日祥. 2004. 华北东部中生代构造体制转折峰期的主要地质效应和形成动力学探讨. 地学前缘, 11(3): 285-297. DOI:10.3321/j.issn:1005-2321.2004.03.027
赵振华, 增田彰正, 夏巴尼 MB. 1992. 稀有金属花岗岩的稀土元素四分组效应. 地球化学, 28(3): 221-233. DOI:10.3321/j.issn:0379-1726.1992.03.003
朱金初, 饶冰, 熊小林, 李福春, 张佩华. 2002. 富锂氟含稀有矿化花岗质岩石的对比和成因思考. 地球化学, 31(2): 141-152. DOI:10.3321/j.issn:0379-1726.2002.02.005
祝新友, 张志辉, 付旭, 李柏阳, 王艳丽, 焦守涛, 孙雅琳. 2016. 内蒙古赤峰维拉斯托大型锡多金属矿的地质地球化学特征. 中国地质, 43(1): 188-208. DOI:10.3969/j.issn.1000-3657.2016.01.014