岩石学报  2021, Vol. 37 Issue (3): 886-910, doi: 10.18654/1000-0569/2021.03.15   PDF    
新疆阿尔泰巴斯铁列克钨多金属矿区花岗质岩体年代学和地球化学特征及其对成矿的启示
包长甲1, 柴凤梅1, 李永2, 祁冬梅1, 张雪冰1, 田猛2     
1. 新疆大学新疆中亚造山带大陆动力学与成矿预测实验室, 乌鲁木齐 830046;
2. 新疆维吾尔自治区有色地质勘查局地质矿产勘查研究院, 乌鲁木齐 830000
摘要: 巴斯铁列克钨多金属矿床位于新疆阿尔泰造山带南缘,是近年来在区内发现的首例中型钨多金属矿床。矿体主要产于二叠纪花岗岩与上志留统-下泥盆统康布铁堡组火山-沉积岩接触带的矽卡岩中。钨矿化与矿区花岗质岩石有明显的空间关系。然而,与钨矿化有关的花岗质岩石成因尚不清楚。本文对矿区出露的与矿化关系密切的黑云母花岗岩、二长花岗岩和二云母花岗岩进行了锆石U-Pb年代学和岩石地球化学研究。3个样品的LA-ICP-MS锆石U-Pb加权平均年龄分别为282.3±3.2Ma、284.3±2.2Ma和284.8±2.3Ma,属早二叠世岩浆活动的产物,与成矿年龄一致。所有岩石具有高硅(SiO2=73.6%~78.3%)、富碱(K2O+Na2O=5.15%~9.62%)富钾(K2O/Na2O>1.1)、贫钙(CaO=0.19%~0.75%)和钛(TiO2=0.04%~0.24%)、弱过铝-强过铝质(A/CNK=1.01~1.39)特征。这些岩石稀土元素总量(∑REE)变化较大(变化于20.3×10-6~328×10-6),但二云母花岗岩显示轻重稀土元素分异不明显((La/Yb)N=0.96~2.06)、Eu强烈负异常(δEu=0.07~0.41)的深"V"型稀土元素分布特征,黑云母花岗岩和二长花岗岩显示轻稀土略富集((La/Yb)N分别为2.8~5.5和4.8~7.4)且Eu负异常(δEu=0.33~0.39和0.34~0.63)明显的右倾型稀土元素分布特征。所有样品均显示相对富集Rb、Th、U、Pb元素和相对亏损Nb、Ti、P、Sr、Ba元素,但二云母花岗岩中W含量(4.6×10-6~9.4×10-6)相对低于黑云母花岗岩和二长花岗岩中W含量(分别为15.1×10-6~168×10-6和8.4×10-6~16.0×10-6)。所有样品的锆石具有正的高εHft)值(+3.8~+11)和相对年轻的亏损地幔模式年龄(Hf的tDM2为0.60~1.0Ga)。以上特征说明,这些岩石属高钾钙碱性分异I-A过渡型花岗岩。结合区域地质背景,认为这些岩体是二叠纪时期后碰撞伸展环境下两个独立岩浆事件的产物,母岩浆均来源于新生地壳熔体与幔源岩浆,经过高度分异演化后结晶形成矿区岩石。花岗质岩浆活动为巴斯铁列克钨矿床提供了成矿物质,岩浆演化过程(结晶分异与熔体-流体作用)对成矿元素有富集作用。
关键词: 锆石U-Pb定年    Hf同位素    岩石地球化学    花岗岩    巴斯铁列克    
Geochronology and geochemical characteristics of granitoids in the Bastielieke tungsten polymetallic deposit in the southern margin of Altay: Implications for tungsten mineralization
BAO ChangJia1, CHAI FengMei1, LI Yong2, QI DongMei1, ZHANG XueBing1, TIAN Meng2     
1. Xinjiang Key Laboratory for Geodynamic Processes and Metallogenic Prognosis of the Central Asian Orogenic Belt, Xinjiang University, Urumqi 830046, China;
2. Institute of Geology and Mineral Resources Exploration, Xinjiang Nonferrous Metals Geoexploration Bureau, Urumqi 830000, China
Abstract: The medium-sized Bastielieke tungsten polymetallic deposit is a newly discovered deposit in the Kelan Basin at the southern margin of Altay Orogenic Belt. The skarn hosted the orebody which is at the outer contact zone between the granitoids and volcanic-sedimentary rocks of the Upper Silurian to Lower Devonian Kangbutiebao Formation. Two-mica granite, biotite granite and monzogranite outcrop in the Bastielieke ore district showing obvious spacial relation to the hosted ore deposits. The petrogenesis of these rocks associated with the tungsten mineralization is indefinitive. LA-ICP-MS zircon U-Pb geochronology and geochemical analysis of major and trace elements and zircon Hf isotope are performed to discuss their petrogenesis. These granites show consistent weighted mean 206Pb/238U ages of 282.3±3.2Ma, 284.3±2.2Ma and 284.8±2.3Ma, respectively, which indicate that the Permian magmatic events at the southern margin of Altay is consistent to the ages of mineralization. These granitoids exhibit high SiO2 concentration (73.6%~78.3%), high total alkali (Na2O+K2O) (5.15%~9.62%) concentration and low CaO (0.19%~0.75%) and TiO2 (0.04%~0.24%) concentrations. In addition, they are enriched in potassium (K2O/Na2O=1.08~2.06) with high A/CNK values (1.01~1.39). The total REE contents of these granites vary from 20.3×10-6 to 328×10-6 with different chondrite-normalized REE patterns. The two-mica granite is characterized by low (La/Yb)N ration (0.96~2.06) with pronounced negative Eu anomaly from 0.07 to 0.41, and showing a deep V-type pattern. The biotite granite and monzogranite display enrichment of LREE with (La/Yb)N ration ranging from 2.8 to 5.5 and 4.8 to 7.4, and moderate negative Eu anomalies ranging from 0.33 to 0.39 and 0.34 to 0.63, respectively. All samples display noticeable negative anomalies of Nb, Ti, P, Sr, Ba and strong positive anomalies of Rb, Th, U and Pb. The two-mica granite has lower W concentrations varying from 4.6×10-6 to 9.4×10-6 than those in the biotite granite and monzogranite from 15.1×10-6 to 168×10-6 and 8.4×10-6 to 16.0×10-6, respectively. In addition, zircon εHf(t) values range from +3.8 to +11, two stage model ages (tDM2) ranging from 0.6Ga to 1.0Ga. The petrologic and geochemical data indicate that these rocks are high-K calc-alkaline and weak peraluminious to strongly peraluminous, which is similar to the geochemical behavior of the fractionated I-A type granitoids. Combined with regional geology, we suggest that they are derived from melting of the juvenile crust and mantle, and followed by subsequent crystallization fractionation in deeper magma chamber in a post-collision environment by two independent magmatic events. These granites provide key ore components for the formation of Bastielieke tungsten polymetallic deposit and preconcentrate tungsten by crystallization fractionation and melts-fluid exsolution.
Key words: Zircon U-Pb dating    Hf isotope    Geochemistry    Granitoids    Bastielieke    

阿尔泰造山带是中亚造山带的重要组成部分,因其独特的地质构造特征和丰富的矿产资源一直是人们关注的热点地区之一。前人研究结果表明,二叠纪是新疆阿尔泰一次重要的构造-岩浆-成矿作用时期(王涛等, 2005, 2010; 肖文交等, 2006; 杨富全等, 2018),区内出露大量二叠纪侵入岩(王涛等, 2010; Gao and Zhou, 2013; Tong et al., 2014; Liu et al., 2017; Tang et al., 2017),发育有一系列不同成因类型的金属矿床,包括造山型金矿(多纳拉萨依金矿和赛都金矿, 闫升好等, 2004)、矽卡岩型铁矿(加尔巴斯套铁矿, 刘锋等, 2012; 库额尔齐斯铁矿, 张超等, 2013)、岩浆期后热液脉型铜矿(克因布拉克铜矿, 杨富全等, 2010; 张志欣等, 2014)、岩浆型铜镍硫化物矿床(喀拉通克铜镍矿, 张作衡等, 2005)、伟晶岩型稀有金属矿床(大喀拉苏Be-Nb-Ta矿, 任宝琴等, 2011)和矽卡岩型钨矿床(巴斯铁列克钨多金属矿床, 李永, 2018; 杨富全等, 2019)等。这些二叠纪侵入岩与成矿作用关系密切,不同类型侵入岩与矿床种属有明显的对应关系,岩浆热液型矿床最为发育。前人对区内与岩浆热液矿床有关的二叠纪花岗质岩石开展过大量研究,也获得了一些重要认识,如这些花岗质岩石主要分布于阿尔泰造山带南缘与额尔齐斯一带,多呈未变形圆形或不规则状产出,沿NW-SE向断裂分布(王涛等, 2010; Cai et al., 2011),持续时间在290~270Ma,峰期约在281~271Ma(王涛等, 2010; Tong et al., 2014);岩石类型主要包括碱性长石花岗岩(Liu et al., 2017)、黑云母花岗闪长岩、二长花岗岩(王涛等, 2005);岩石普遍富碱、贫硅和铁,多属亚碱性系列岩石,部分为钾玄质系列,铝饱和指数变化较大,多为准铝质-过铝质(王涛等, 2010),主要为I-A过渡型花岗岩(王涛等, 2005; Tong et al., 2014; 李强等, 2019)和A2型花岗岩(Shen et al., 2011; Tong et al., 2014; Liu et al., 2017)。然而,对于这些二叠纪时期侵入岩形成的构造背景,尚存在有后碰撞伸展环境(Tong et al., 2014; Liu et al., 2017)或俯冲环境(Yuan et al., 2007)的争议。

巴斯铁列克钨多金属矿床是2013年新疆有色地质勘查局地质矿产勘查研究院在克兰盆地内实施找矿勘查任务时发现的,是新疆阿尔泰境内首次发现的富钨矿床,目前矿床的钨储量可达中型(李永, 2018)。巴斯铁列克钨多金属矿床矿体赋存于二叠纪花岗岩与晚志留世-早泥盆世康布铁堡组外接触带的矽卡岩内,形成于二叠纪(284.4Ma, 杨富全等, 2019), 张国锋等(2019)对矿床的矿物学特征开展了研究,认为矿区矽卡岩是岩浆热液交代大理岩的产物。矿区内出露的花岗质岩石类型丰富,但是与钨相关的花岗质岩石特征尚不清楚,针对与钨有关的花岗质岩石开展岩相学、地球化学以及岩石成因研究,将有助于区内找矿勘查工作部署。

本文对巴斯铁列克钨多金属矿区发育的花岗质岩石开展系统的岩石地球化学及同位素年代学研究,利用LA-ICP-MS锆石U-Pb定年法约束岩体形成时代,结合岩石地球化学和锆石Hf同位素分析,探讨岩浆来源和演化。在此基础上,结合区域地质演化历史的研究,限定了该时期岩浆活动与矿床形成的构造背景。这将为巴斯铁列克钨多金属矿床成因研究提供重要信息,也为揭示阿尔泰造山带南缘构造演化提供依据。

1 区域地质概况

阿尔泰造山带是中亚造山带的重要组成部分,呈北西-南东向横贯中、俄、哈、蒙四国,构造上处于西伯利亚板块和哈萨克斯坦-准噶尔板块接合部位。阿尔泰造山带被北西-南东向断裂大致分为北、中、南3个块体(Xiao et al., 2004)。巴斯铁列克钨多金属矿床位于南阿尔泰块体内的克兰盆地(图 1)。盆地内出露的地层有中-上志留统库鲁木提群变质岩;上志留-下泥盆统康布铁堡组变质中酸性火山沉积岩;中泥盆统阿勒泰镇组浅变质火山沉积岩。盆地内已发现有阿巴宫铁矿、铁木尔特铅锌铜矿、塔拉特铅锌铁矿、恰夏铁铜矿、萨热阔布金矿、红墩铅锌矿等(图 1)。花岗质岩浆侵入活动以奥陶纪和二叠纪最为剧烈,其中奥陶纪花岗岩出露面积不大,多已发生变形,岩石类型有黑云母花岗岩、二云母花岗岩和白云母花岗岩,属强过铝质高钾钙碱性岩石(刘锋等, 2009; 柴凤梅等, 2010),是俯冲作用的产物。二叠纪侵入岩出露面积相对较广,集中形成于286~261Ma(王涛等, 2005; 刘锋等, 2009; 李强等, 2019),岩体基本不变形,多呈不规则圆形、椭圆形,多切割区域构造线(王涛等, 2005),主要岩石类型为黑云母二长花岗岩、二云母花岗岩和钾长花岗岩,具有高钾钙碱性岩石特征,属I-A过渡型和A型花岗岩(王涛等, 2005; 童英等, 2010; 刘锋等, 2012; Liu et al., 2017; 李强等, 2019)。

图 1 阿尔泰造山带南缘区域地质简图及克兰盆地地质略图(据Windley et al., 2002; Li et al., 2003; 有色地质勘查局706队, 2003修改) 矿床:1-萨热阔布金矿;2-恰夏铁铜矿;3-铁木尔特铅锌铜矿;4-阿巴宫铁矿;5-塔拉特铅锌铁矿;6-红墩铅锌矿;7-巴斯铁列克钨多金属矿L1矿段;8-巴斯铁列克钨多金属矿L2矿段;9-巴斯铁列克钨多金属矿L3矿段 Fig. 1 Regional geological sketch maps of the southern margin of the Altay orogenic belt and geological sketch map of the Kelang Basin (modified after Windley et al., 2002; Li et al., 2003) Names of ore deposits: 1-Sarekuobu Au deposit; 2-Qiaxia Fe-Cu deposit; 3-Tuomoerte Pb-Zn-Cu deposit; 4-Abagong Fe deposit; 5-Talate Pb-Zn-Fe deposit; 6-Hongdun Pb-Zn deposit; 7-L1 oreblock of Bastielieke tungsten polymetallic deposit; 8-L2 oreblock of Bastielieke tungsten polymetallic deposit; 9-L3 oreblock of Bastielieke tungsten polymetallic deposit

①   新疆维吾尔自治区有色地质勘查局706队. 2003. 新疆阿尔泰市铁木尔特地区区域地质报告

2 矿区地质及岩体地质

巴斯铁列克矿床位于克兰盆地东南(图 2)。矿区出露地层为上志留-下泥盆统康布铁堡组第二岩性段浅粒岩、变粒岩、大理岩和黑云石英片岩。矿体产于花岗岩体与康布铁堡组接触带的矽卡岩中。地表出露的矿化带长约16km,宽50~150m,共圈出了3条钨矿体、10条钨铜(锌)矿体、4条铜矿体和1条锌矿体。目前矿化带由西北向东南划分为3个矿段,分别为L1矿段(原为克兹尔尕因W-Cu矿段)、L2矿段(原为巴斯铁列克Cu-Zn-W矿段)和L3矿段(原为巴斯铁列克东W矿段)。

图 2 巴斯铁列克钨多金属矿区地质图(据新疆维吾尔自治区有色地质勘查局地质矿产勘查研究院, 2019修改) Fig. 2 Geological sketch maps of the Bastielieke tungsten polymetallic deposit

②  新疆维吾尔自治区有色地质勘查局地质矿产勘查研究院. 2019. 新疆福海县巴斯铁列克铜多金属矿预查报告

矿区出露的花岗质岩石在地表呈不规则状岩株和岩枝产出,岩体边部呈港湾状,面积约7.74km2,主要分布于阿巴宫断裂北侧,沿NW-SE向延伸,侵位于晚志留-早泥盆世火山-沉积岩中,与地层呈断层接触。在巴斯铁列克不同矿段出露的岩石类型略有差异,其中L1矿段出露有二云母花岗岩,呈岩株状;L2矿段出露有中-粗粒黑云母花岗岩和二长花岗岩,其中黑云母花岗岩与地层接触带发育矽卡岩化蚀变和钨铜矿化,在二长花岗岩中见细脉状钨矿化。各类型岩石特征如下。

二云母花岗岩呈灰白色,花岗结构,块状构造(图 3a, b)。主要由石英(30%)、钾长石(约35%)、酸性斜长石(20%)、少量黑云母和白云母(15%)组成。其中石英约0.2 ~0.8mm,呈他形粒状,偶见微弱波状消光。钾长石主要为微斜长石(约占总体积的25%)和条纹长石(约占总体积的10%),呈半自形-他形板状,分别发育格子状双晶和条纹结构;酸性斜长石主要为钠长石(牌号为8~10,约占总体积的15%)和更长石(约5%)为主,呈长板状,自形程度与钾长石相当,长约0.3~0.6mm,聚片双晶发育;白云母和黑云母含量相当,主要呈片状分布于石英和长石间,片径大小约0.1~0.5mm。副矿物有锆石、磷灰石和磁铁矿(< 1%),样品偶见石榴子石和黑电气石。

图 3 巴斯铁列克钨多金属矿区花岗质岩石照片 (a、b)二云母花岗岩;(c、d)黑云母花岗岩;(e、f)二长花岗岩. Q-石英;Kfs-钾长石;Pl-斜长石;Bt-黑云母;Ms-白云母;Gt-石榴子石 Fig. 3 Photographs and microphotographs of different intrusions in the Bastielieke tungsten polymetallic deposit (a, b) two mica granite; (c, d) biotite granite; (e, f) monzongranite. Q-quartz; Kfs-K-feldspar; Pl-plagioclase; Bt-biotite; Ms-muscovite; Gt-garnet

黑云母花岗岩呈灰白色,中粗粒花岗结构,块状构造(图 3c, d)。主要由碱性长石、斜长石和石英组成。其中石英含量约35%,多呈2~6mm他形粒状;钾长石含量约35%,多呈半自形板状-他形粒状,介于2~7mm之间,以微斜长石为主,也有少量正长石和条纹长石;酸性斜长石有钠长石和更长石,含量约15%,呈半自形板状,多已发生绢云母化蚀变;黑云母呈片状,片径大小约0.1~1mm,含量约13%,解理发育,局部蚀变为绿泥石。副矿物主要为锆石、磷灰石和钛铁矿等。显微镜下可见约2%的自形粒状石榴子石和他形白钨矿。

二长花岗岩呈灰白色,中粒结构,块状构造(图 3e, f)。石英(30%)、碱性长石(35%)、酸性斜长石(25%)和黑云母(8%)组成。其中石英均呈他形粒状分布于其它颗粒间,大小约2~3mm;钾长石主要为格子状双晶发育的微斜长石,呈半自形板状-他形粒状,介于3~5mm之间,偶见包裹他形粒状石英和板条状斜长石,多发生高岭土化,偶见微文象结构;酸性斜长石主要为钠长石,呈半自形板状,约3~4mm,聚片双晶发育, 多已发生绢云母化蚀变;黑云母呈他形片状,大小约0.5~1mm左右,偶见边缘和解理缝蚀变。副矿物主要为锆石、磷灰石、榍石和钛铁矿,显微镜下见自形石榴子石。

3 样品采集和分析方法 3.1 样品采集

本次研究所用样品采自巴斯铁列克矿区地表出露的新鲜-弱蚀变的岩石。主要为二云母花岗岩(采样坐标:47°39′13.86″N、88°34′57.83″E)、黑云母花岗岩(采样坐标:47°36′39.25″N、88°38′17.92″E)和二长花岗岩(采样坐标:47°36′43.22″N、88°38′16.44″E)。其中二云母花岗岩(5件)、黑云母花岗岩(6件)、二长花岗岩(5件)用于全岩主微量元素成分分析,不同类型岩石各一件用于锆石定年分析和Hf同位素分析。具体采样位置见图 2

3.2 分析方法

全岩主微量元素成分分析在澳实分析检测(广州)有限公司澳实矿物实验室完成。主量元素采用P61-XRF26FS方法,所用仪器为X荧光光谱仪(型号PANalytical AXIOS),分析精度优于0.01%;FeO含量用滴定法测定,按照GB/T14506.14—2010标准执行,分析精度优于0.1%。微量和稀土元素分析采用M61-MS81(等离子质谱/光谱定量综合分析)方法,所用仪器为电感耦合等离子体质仪(型号ELAN9000), 不同元素的分析精度和准确度略有不同,其中Th、U、Ge、Sb为0.05×10-6,Ga、Sr、Nd、Ta、Sc、Co和W为0.1×10-6,Rb、Nb、Hf、Li、Cu、Ni和Sn为0.2×10-6,Ba和Pb为0.5×10-6,Cr、V和Zn为1×10-6;稀土元素中La、Ce和Y的分析精度和准确度为0.1×10-6,其余元素优于0.03×10-6

锆石单矿物挑选在河北省地质测绘院岩矿实验测试中心完成,经破碎、重选和磁选后,在双目镜下挑选出晶形和透明度较好的颗粒。锆石的制靶、透反射、阴极发光(CL)图像等工作由南京宏创地质勘查技术服务有限公司完成。定年测试前圈出内部无包体、表面无裂纹、环带发育的晶体颗粒及最佳分析部位。

锆石LA-ICP-MS U-Pb同位素分析在北京科融恩科技有限责任公司完成。所用仪器为Agilent 7900型ICP-MS等离子体质谱仪和与之配套的ESI NWR UP213激光剥蚀系统。激光剥蚀过程中采用氦气作载气、氩气为补偿气以调节灵敏度,激光剥蚀斑束直径为40μm,频率为10Hz。每测试9个样品点测定2次标准锆石91500,锆石标准样品Plesǒvice作为外标,用于监测仪器状态和分析数据质量。本次测试分析质量监控样品Plešovice锆石共计11点,测定年龄为336.9±2.1Ma。样品的同位素比值和元素含量计算采用GLITTER(Version 4.0)处理,锆石的U-Pb年龄计算及图的绘制采用Isoplot(Version 3.7)完成。单次测量数据点的误差为1σ,年龄值选用206Pb/238U年龄,加权平均值有95%的置信度。

锆石的原位Hf同位素分析在北京科融恩科技有限责任公司完成。所用仪器为ThermoFisher Neptune Plus电感耦合等离子体质谱仪(LA-MC-ICP-MS)与之配套的ASI J200飞秒激光剥蚀系统。激光剥蚀过程采用氦气作载气,氩气作补偿气,激光剥蚀脉冲频率为10Hz,激光束斑直径为40μm,每个分析点的数据采集时间为40秒。采用标准锆石(91500、Plesǒvice、Temora或GJ-1)作监控标样,每分析10个样品点,分析1次标准样品,以监控仪器状态。此次锆石原位Hf同位素比值测定选择了锆石U-Pb年龄谐和度好的颗粒,在靠近U-Pb年龄分析位置圈点测定。εHf(t)值根据测定的锆石U-Pb年龄计算获得,采用的参数为:现今球粒陨石的176Lu/177Hf为0.0332,176Hf/177Hf为0.282772(Blichert-Toft and Albarède, 1997),亏损地幔的176Lu/177Hf为0.0384,176Hf/177Hf为0.28325 (Vervoot and Blichert-Toft, 1999),大陆平均地壳的176Lu/177Hf为0.015(Griffin et al., 2002),Lu衰变常数λ=1.867×10-11/y(Söderlund et al., 2004)。

4 分析结果 4.1 全岩主量、微量元素

巴斯铁列克钨多金属矿区二云母花岗岩、黑云母花岗岩和二长花岗岩的主量和微量元素分析结果列于表 1。这些岩石的元素成分略有不同。它们具有高SiO2含量(均高于73%)和全碱含量,但二云母花岗岩的全碱含量(K2O+Na2O=8.53%~9.62%)明显高于黑云母花岗岩和二长花岗岩的全碱含量(分别为7.11%~7.70%和5.15%~7.77%),除个别样品外,大部分样品明显富钾((K2O/Na2O=1.08~2.06)。在SiO2-K2O图解上(图 4a),除1个二长花岗岩样品位于钙碱性岩区外,其余所有样品均落于高钾钙碱性岩区。二云母花岗岩的TiO2、MgO、CaO含量明显低于黑云母花岗岩和二长花岗岩的相应元素含量。所有岩石属过铝质岩(图 4b),但二云母花岗岩的铝饱和指数(A/CNK介于1.01~1.15)明显小于黑云母花岗岩和二长花岗岩的铝饱和指数(A/CNK)(分别为1.20~1.30和1.19~1.39),且它们的Al2O3含量较钙碱性强过铝质花岗岩低(周刚, 2006)。前人研究结果表明,利用花岗岩锆石饱和温度可以获得花岗岩初始岩浆温度。巴斯铁列克矿区花岗质岩石中见有捕获/继承锆石,说明母岩浆中Zr含量已达到饱和,利用锆石饱和温度计获得660~712℃、741~788℃和814~855℃分别代表了二云母花岗岩、黑云母花岗岩和二长花岗岩初始岩浆温度的上限(Watson and Harrison, 1983; Miller et al., 2003)。在SiO2与其他氧化物的Harker图解上(图 5),黑云母花岗岩与二长花岗岩有一定的演化关系,SiO2与Al2O3、CaO、TiO2、FeOT、MgO有明显的负相关,与P2O5、K2O和Na2O关系不很明显,表明不同类型花岗岩为花岗岩浆经历不同程度结晶分异过程的产物。

表 1 巴斯铁列克钨多金属矿区花岗质岩石的主量(wt%)和微量(×10-6)元素组成 Table 1 Major (wt%) and trace (×10-6) element data for granitoids in the Bastielieke tungsten polymetallic deposit

图 4 巴斯铁列克钨多金属矿区花岗质岩石SiO2-K2O图(a, 据Peccerillo and Taylor, 1976)和A/CNK-A/NK图(b, 据Maniar and Piccoli, 1989) Fig. 4 Diagrams of SiO2 vs. K2O diagram (a, after Peccerillo and Taylor, 1976) and A/CNK vs. A/NK (b, after Maniar and Piccoli, 1989) for granitoids in the Bastielieke tungsten polymetallic deposit

图 5 巴斯铁列克钨多金属矿区花岗质岩石哈克图解 Fig. 5 Harker plots for granitoids in the Bastielieke tungsten polymetallic deposit

在微量元素成分上,不同类型岩石也显示了不同的特征。二云母花岗岩的稀土元素总量(∑REE)(20.32×10-6~33.54×10-6)明显低于黑云母花岗岩(69.7×10-6~123.3×10-6)和二长花岗岩(153.8×10-6~327.6×10-6)的稀土元素总量。在球粒陨石标准化稀土元素图解上(图 6a),二云母花岗岩显示了轻重稀土元素分异不明显((La/Yb)N=0.96~2.06),Eu强烈负异常(δEu=0.07~0.41)的深“V”型稀土元素分布型式;黑云母花岗岩(图 6c)和二长花岗岩(图 6e)显示了相似的轻稀土略富集((La/Yb)N分别为2.8~5.5和4.79~7.41)的Eu负异常(δEu分别为0.33~0.39和0.34~0.63)明显的右倾型稀土元素分布特征。在原始地幔标准化微量元素蛛网图上(图 6b, d, f), 三者均显示了Ti、P、Sr、Ba的明显负异常,Rb、Th、U、K、Pb、Zr、Hf的正异常,但元素含量上有明显差别。二云母花岗岩的Sr含量(6.0×10-6~13.8×10-6)和Ba含量(2.5×10-6~29.9×10-6)低于黑云母花岗岩与二长花岗岩的Sr含量(分别为32.0×10-6~57.6×10-6和55.8×10-6~85.3×10-6)和Ba含量(分别为88.2×10-6~166×10-6和121.5×10-6~338×10-6)。此外,二云母花岗岩有相对低的W含量(4.6×10-6~9.4×10-6),二长花岗岩具有中等的W含量(8.4×10-6~16.0×10-6),黑云母花岗岩有相对高的W(15.1×10-6~168×10-6)含量。

图 6 巴斯铁列克钨多金属矿区花岗质岩石球粒陨石标准化稀土元素配分图(a, c, e)和原始地幔标准化微量元素蛛网图(b, d, f)(原标准化值据Sun and McDonough, 1989) Fig. 6 Chondrite-normalized REE patterns (a, c, e) and primitive mantle-normalized trace element patterns (b, d, f) of the whole rock samples from granitoids in the Bastielieke tungsten polymetallic deposit (normalization values after Sun and McDonough, 1989)
4.2 锆石U-Pb年龄

用于锆石LA-ICP-MS U-Pb同位素定年的3件样品中,所有锆石晶形较好,均呈半自形-自形柱状及双锥状,晶棱及晶面清楚,晶体大小变化较大,其中二云母花岗岩中锆石长轴变化于100~150μm之间,长短轴比一般为2:1~3:1;黑云母花岗岩中锆石长轴变化于100~200μm之间,长短轴比一般为3:1~4:1左右;二长花岗岩中锆石长轴变化于80~200μm之间,长短轴比一般为2:1~3:1。大部分锆石颗粒具有岩浆锆石的典型振荡环带,个别锆石可见核幔结构(图 7)。测试结果列于表 2

图 7 巴斯铁列克钨多金属矿区花岗质岩石锆石阴极发光图,示测试位置 Fig. 7 Cathodoluminescence images of zircon from granitoids in the Bastielieke tungsten polymetallic deposit, showing locations of analysis

表 2 巴斯铁列克钨多金属矿区花岗质岩石LA-ICP-MS锆石U-Pb年龄测定结果 Table 2 LA-ICP-MS U-Pb dating data of zircons from granitoids in the Bastielieke tungsten polymetallic deposit

二云母花岗岩(样品号BKL1-32)的30个分析结果中,23粒锆石分析点的Th/U比值介于0.13~0.67之间,它们的206Pb/238U表面年龄变化范围不大,介于273.5~288.6Ma之间,在误差范围内有一致的206Pb/238U、207Pb/235U和207Pb/206U值,数据点聚集在年龄谐和线附近较小的范围内,其206Pb/238U加权平均年龄为282.3±3.2Ma(MSWD=0.22)(图 8a),与谐和年龄值281.8±1.6Ma(MSWD=1.5)一致,可以代表二云母花岗岩侵位年龄。3个分析点的较老年龄(959.2Ma、489.8Ma、393.7Ma)分别与阿尔泰古老的陆块年龄值(Wang et al., 2006; 刘伟等, 2010)、晚寒武世-早奥陶世地层年龄值(481~495Ma)(杨富全等, 2017)和泥盆纪岩浆岩年龄值(381~400Ma)(杨富全等, 2010; 单强等, 2011; 柴凤梅等, 2012)一致,应该代表了继承或捕获的早期岩浆锆石的年龄。4个置信度较低的分析点(07点274.0Ma, 10点288.6Ma, 16点276.9Ma, 27点281.7Ma)与其他23个分析点的年龄在误差范围内重合,但它们的207Pb/235U较大,可能主要是207Pb难以测准或者Pb丢失导致的。

图 8 巴斯铁列克钨多金属矿区花岗质岩石LA-ICP-MS锆石U-Pb年龄谐和线和加权平均年龄分布图 Fig. 8 Concordia diagrams and distributions of LA-ICP-MS U-Pb weighted mean age from granitoids in the Bastielieke tungsten polymetallic deposit

黑云母花岗岩(样品号BKL2-7)中30个分析点中,除1个年龄较老的样品点(24点322.7Ma)和年龄值偏低的样品点(12点254.6Ma)之外,其余28个分析点的Th/U比值介于0.13~1.25之间(表 2),并且它们的206Pb/238U年龄介于275.4~289.9Ma之间,206Pb/238U、207Pb/235U和207Pb/206U值在误差范围内一致,数据点聚集在谐和线附近较小的范围内,其206Pb/238U加权平均值为284.3±2.2Ma(MSWD=0.29),与谐和年龄值284.2±1.1Ma(MSWD=0.9)一致(图 8b),可以代表岩体的形成年龄。

二长花岗岩(样品号BKL2-32)30个分析点中,除1个年龄值偏年轻的样品点(10点249.9Ma)和4个年龄值偏老的样品点(07点333.7Ma, 15点307.0Ma, 26点336.7Ma, 28点476.8Ma)之外,其余25个分析点的Th/U比值介于0.13~1.17之间(表 2),它们的206Pb/238U年龄结果变化范围较小,介于281.3~290.7Ma之间,206Pb/238U、207Pb/235U和207Pb/206U值在误差范围内一致,数据点聚集在谐和线附近较小的范围内,谐和年龄值(284.7±1.2Ma, MSWD=0.12)与206Pb/238U加权平均年龄(284.8±2.3Ma, MSWD=0.31)(图 8c)一致,可以代表二长花岗岩的侵位年龄。其中较老的年龄值(476.8Ma),可能代表了捕获或继承的早期岩浆锆石的年龄,较小的年龄值可能是由于Pb丢失造成的。

综上所述,本次研究的3个样品的年龄数据可靠,谐和年龄值与加权平均年龄值接近(图 8),表明这些锆石形成后U-Pb体系保持封闭,没有明显的U或Pb同位素的丢失和加入。这些年龄可以代表巴斯铁列克钨多金属矿区花岗质岩体的结晶年龄。

4.3 锆石Hf同位素组成

二云母花岗岩、黑云母花岗岩和二长花岗岩中已定年锆石的21个点、23个点和24个点的Hf同位素分析结果列于表 3

表 3 巴斯铁列克钨多金属矿区花岗质岩石中锆石Lu-Hf同位素组成 Table 3 Zircon Hf isotopic composition of the granitoids in the Bastielieke tungsten polymetallic deposit

所有测试锆石的176Hf/177Hf误差值在0.00003以内,表明Hf同位素数据质量较好。二云母花岗岩的锆石176Lu/177Hf比值介于0.001172~0.003907之间,锆石具有一致的176Hf/177Hf初始比值,21个点的(176Hf/177Hf)i比值变化于0.282754~0.282882,εHf(t)相对集中且较高,介于+5.6~+10之间,fLu/Hf变化于-1.0~-0.88,对应的两阶段Hf模式年龄(tDM2)范围为0.66~0.95Ga。黑云母花岗岩的锆石176Lu/177Hf比值主要介于0.001730~0.004719之间,23个点的(176Hf/177Hf)i比值变化于0.282716~0.282905之间,εHf(t)介于+4.2~+11之间,fLu/Hf变化于-0.95~-0.86,其两阶段Hf模式年龄(tDM2)范围为0.60~1.0Ga。二长花岗岩的锆石176Lu/177Hf比值主要介于0.001720~0.005233之间,24个点的(176Hf/177Hf)i比值变化于0.282702~0.282903之间,εHf(t)介于+3.8~+11之间,fLu/Hf变化于-0.92~-0.86,其两阶段Hf模式年龄(tDM2)范围为0.61~1.0Ga(表 3)。

总之,3个岩体的锆石Hf同位素组成稍有差异,但它们的εHf(t)均为大于零的正值,与新疆阿尔泰造山带二叠纪花岗岩锆石具有高(正)εHf(t)值特征相同,如大桥南岩体(+6.1~+8.0)、锡伯渡岩体(+10~+12.9)、布尔津岩体(+5.6~+12.9)(Tong et al., 2014)、恰库尔特岩体(+10~+18.7)(Shen et al., 2011)。fLu/Hf均明显小于硅铝质地壳的fLu/Hf(-0.72)值(Vervoort et al., 1996),两阶段Hf模式年龄(tDM2)远高于岩体的结晶年龄(284Ma),可以代表各岩体的源岩脱离亏损地幔的时间。

5 讨论 5.1 矿区花岗岩年代学及指示意义

阿尔泰造山带南缘古生代花岗岩广为发育,随着近年来对阿尔泰造山带内不同时期、成因有别的花岗岩的深入研究,认为古生代花岗质岩浆活动主要发生在470~440Ma(中晚奥陶世)、425~360Ma(晚志留世-晚泥盆世)、355~318Ma(晚石炭世)、290~270Ma(早二叠世)、245~190Ma(中晚三叠世-早侏罗世)共5个期次,其中晚志留世-晚泥盆世和早二叠世花岗岩最为发育(王涛等, 2010)。

本研究获得巴斯铁列克钨多金属矿区二云母花岗岩、黑云母花岗岩与二长花岗岩的锆石LA-ICP-MS U-Pb年龄结果接近,分别为282.3±3.2Ma,284.3±2.2Ma和284.8±2.3Ma,与前人获得的阿尔泰造山带内喇嘛昭岩体(276±9Ma, 王涛等, 2005)、将军山岩体和大喀拉苏岩体(268.3±1.9Ma, 270.4±1.9Ma和261.4±2.1Ma, Liu et al., 2017; 李强等, 2019)、阿苇滩岩体和大桥南岩体(271±2Ma和267±5Ma, Tong et al., 2014)、沙尔布拉克南岩体(275.1±1.7Ma, 孙桂华等, 2009)、可斯尔鬼岩体(286.6±2.6Ma, 刘锋等, 2012)、克因布拉克岩体(278.6±3.5Ma, 李香仁等, 2012)、克孜勒岩体(289.4±4.3Ma, 赵玉梅等, 2016)等年龄结果一致,并且矿区花岗质岩石的锆石U-Pb年龄与矿区内辉钼矿Re-Os同位素年龄结果(284.4Ma,杨富全等, 2019)在误差范围内一致,指示这些花岗质岩石侵位与矿区内钨多金属成矿作用均发生于早二叠世,与区域内不同成因类型的金属矿床成矿时限一致(闫升好等, 2004; 张作衡等, 2005; 杨富全等, 2010; 刘锋等, 2012; 张超等, 2013),是阿尔泰造山带南缘二叠纪成岩成矿作用的组成部分。

5.2 岩石类型及关系

巴斯铁列克钨多金属矿区内二云母花岗岩、黑云母花岗岩以及二长花岗岩之间未见到明确的接触关系,不同岩石类型间的侵位序列尚不清楚。不同类型花岗质岩石有一致的锆石U-Pb年龄,表明它们为同一时期岩浆作用的产物。这些花岗质岩石出现白云母和石榴子石矿物组合,SiO2含量高,贫FeO、MgO、TiO2和CaO含量,富碱高钾(K2O>Na2O),属高钾钙碱性过铝质岩石(A/CNK>1.0),明显富集Rb、Th、U和Pb,相对亏损Sr、Ba、Eu、P、Ti和Zr等元素,有相近的Hf同位素组成,暗示它们可能存在某种成因联系。

巴斯铁列克矿区所有花岗岩样品具有低的K/Rb比值(88~131),Zr/Hf比值(16.8~36.6)低于正常花岗岩的相应值(38),暗示这些岩石是母岩浆经过分异形成。有研究表明,低的K/Rb比值是高程度岩浆演化体系的特征(Irber, 1999; 苏慧敏和蒋少涌, 2017);在准铝质或过铝质花岗岩中,随着矿物的分离会导致Zr/Hf比值的降低(Linnen and Keppler, 2002; 郭春丽等, 2017; Yuan et al., 2018, 2019)。在(Zr+Nb+Ce+Y)-(K2O+Na2O)/CaO和(Zr+Nb+Ce+Y)-FeOT/MgO图上(图 9),大部分样品位于分异花岗岩区。

图 9 巴斯铁列克钨多金属矿区花岗质岩石类型分类图 (a) (Zr+Nb+Ce+Y)-(K2O+Na2O)/CaO图;(b) (Zr+Nb+Ce+Y)-FeOT/MgO图(底图据Whalen et al., 1987Eby, 1990). FG-分异I&S型花岗岩;OGT-未分异I & S型花岗岩 Fig. 9 Type discrimination diagrams of granitoids in the Bastielieke tungsten polymetallic deposit (a) (Zr+Nb+Ce+Y) vs. (K2O+Na2O)/CaO; (b) (Zr+Nb+Ce+Y) vs. FeOT/MgO (after Whalen et al., 1987; Eby, 1990). FG-field for fractioned I- and S-type granitoids, OGT-field for unfractioned I- and S-type granitoids

二云母花岗岩与黑云母花岗岩和二长花岗岩在岩石学和地球化学上有明显不同:(1)二云母花岗岩中副矿物为磁铁矿,也见少量黑电气石,黑云母花岗岩与二长花岗岩副矿物中发育钛铁矿;(2)二云母花岗岩有较低的稀土总量(∑REE=20.32×10-6~33.54×10-6),呈现轻重稀土元素分馏不明显((La/Yb)N=0.96~2.06)和Eu负异常强烈的(δEu=0.07~0.41)深“V”型稀土元素配分模式,Ba、Sr和Ti亏损程度更高;(3)二云母花岗岩的A/CNK值(1.01~1.15)较黑云母花岗岩与二长花岗岩的A/CNK值(1.19~1.39)低,但Al2O3、Na2O和全碱(K2O+Na2O)含量明显高于后两者的相应元素含量;(4)在哈克图解上(图 5),二云母花岗岩与黑云母花岗岩多表现了相互平行的元素变化趋势。这些特征表明,二云母花岗岩与黑云母花岗岩和二长花岗岩有不同的成因。

二云母花岗岩具有极高的Rb/Sr比值(29.1~55.2),极低的Nb/Ta比值(4.3~5.6)和Zr/Hf比值(16.8~20.7),这与高分异花岗岩特征一致(郭春丽等, 2017; 吴福元等, 2017; Yuan et al., 2018, 2019)。岩石中电气石的存在、较低的锆石饱和温度、Al2O3与SiO2明显的负相关、低的Zr含量等说明母岩浆经历了较充分的长石和锆石的分离结晶过程,因为B等挥发分的存在可以导致固相线温度下降, 进而延长岩浆的分离结晶过程(陈伟等, 2018)。这可能与岩浆结晶晚期熔体-流体的强烈相互作用有关(Jahn et al., 2001; 郭春丽等, 2017)。黑云母花岗岩与二长花岗岩有相同的矿物组合,有一致的微量元素和稀土元素配分模式以及Hf同位素组成,元素间也存在一定的演化关系(图 5),表明它们系同源岩浆演化的产物。与黑云母花岗岩相比,二长花岗岩有相对低的SiO2含量和Rb/Sr(3.8~4.7)比值(黑云母花岗岩Rb/Sr=5.6~9.9),有高的锆石饱和温度(814~855℃),相对高的稀土元素总量、LREE/HREE比值和δEu值,暗示黑云母花岗岩是二长花岗岩分异演化的结果。有研究表明,当花岗质岩浆向更酸性的方向发生分异演化时,它们的Rb/Sr比值会逐渐升高,稀土元素总量会减少,LREE/HREE比值和δEu会降低(Ballouard et al., 2016; 郭春丽等, 2017)。Al2O3、CaO、TiO2和FeOT与SiO2呈负相关,说明岩浆在演化过程中经过了锆石、钛铁矿以及长石的分离结晶。

巴斯铁列克矿区花岗质岩石均出现石榴子石和白云母等S型花岗岩中典型的铝过饱和矿物。然而,有研究指出,经过分异演化后的I型、S型和A型花岗岩往往有较为相近的地球化学特征(Chappell, 1999; 吴福元等, 2007);在岩浆演化过程中,由于受初始岩浆成分、岩浆结晶温压条件和岩浆-热液作用影响,石榴子石与白云母这种富铝矿物组合也可以出现在分异I型花岗岩中(Chappell, 1999; 吴福元等, 2007),并且斜长石、普通角闪石和黑云母的分馏作用可以使I型花岗岩演化为弱过铝质-强过铝质特征岩石(Zen et al., 1986; 李世超等, 2016)。矿区花岗质岩石显示了高10000×Ga/Al值(2.5~3.3)和FeOT/(FeOT+MgO) 值(>0.85),强烈亏损Ba等A型花岗岩特征,但岩石中无碱性暗色矿物,低的Zr+Ce+Nb+Y含量(< 350×10-6)和(K2O+Na2O)/CaO比值(图 9)(Whalen et al., 1987),低的岩浆温度(A型花岗岩通常>880℃)(刘昌实等, 2003)等,与典型A型花岗岩不同。矿区花岗质岩石中不含角闪石,有钠长石和黑云母,除个别样品中P2O5含量低(< 0.1%)之外,大多数样品中P2O5含量高(>0.1%),且与SiO2含量没有明显的相关性,Eu亏损明显的稀土元素配分模式、高的Rb/Sr比值(3.8~55.2)和Rb/Ba比值(1.08~132)(王涛等, 2005)等与I型花岗岩不同。因此,巴斯铁列克矿区花岗质岩石兼具I型、S型和A型花岗岩特征,应该是分异的钙碱性I-A过渡型花岗岩。这与新疆阿尔泰造山带二叠纪花岗岩特征一致(王涛等, 2010; 李强等, 2019),与中国境内中亚造山带高分异花岗岩特征一致(Jahn et al., 2001; Wu et al., 2003)。

综上所述,巴斯铁列克矿区二云母花岗岩、黑云母花岗岩和二长花岗岩同属于二叠纪岩浆作用形成的分异花岗岩,它们代表了同一时期两个独立岩浆事件。二云母花岗岩与黑云母花岗岩和二长花岗岩不具演化关系,黑云母花岗岩和二长花岗岩系同源岩浆不同演化阶段的产物,黑云母花岗岩是二长花岗岩进一步演化的结果。它们可能是具有相同同位素组成、不同化学成分和地球化学特征的源区岩石形成熔体后,经过不同程度的结晶分异形成。

5.3 岩石成因

前人对阿尔泰造山带构造演化进行过大量研究, 认为阿尔泰造山带大致从早古生代(奥陶纪-泥盆纪)开始发生俯冲、碰撞、增生,至石炭纪基本奠定了构造格架(Windley et al., 2002; Xiao et al., 2004; Wang et al., 2006; 王涛等, 2010)。然而,阿尔泰造山带俯冲与碰撞构造转换的具体时间(Gao and Zhou, 2013; 沈晓明等, 2013; Zhang et al., 2018; Liu et al., 2019; Yang et al., 2020)、阿尔泰造山带南缘后碰撞的起始时间(王京彬和徐新, 2006; 龙晓平等, 2006; 肖文交等, 2006; Xiao et al., 2009; 沈晓明等, 2013)以及晚石炭世-早二叠世岩浆作用的动力学机制等一直是争议的焦点。

巴斯铁列克钨多金属矿区二云母花岗岩、黑云母花岗岩与二长花岗岩属二叠纪岩浆活动的产物,它们应该与区域上广泛发育的同时代岩体具有相同的构造背景。尽管对于晚石炭世-早二叠世岩浆作用的动力学机制存在有板片窗(Gao and Zhou, 2013; Zhang et al., 2018)、岩石圈拆沉和地幔柱的争议(Pirajno, et al., 2008; Tang et al., 2017)。越来越多的研究成果表明,晚石炭世-早二叠世阿尔泰造山带南缘处于后碰撞伸展环境(Shen et al., 2011; Gao and Zhou, 2013; Liu et al., 2017)。如晚石炭世-早二叠世陆相火山-沉积岩角度不整合覆盖于早期海相沉积岩之上(周刚等, 2006; Tong et al., 2014);晚石炭世发育淡色花岗岩(沈晓明等, 2013);二叠纪时期存在高温-超高温低压变质作用(仝来喜, 2014; 刘兆和仝来喜, 2015; Liu et al., 2017)、出露大面积未变形的I-A型和A型花岗岩(王涛等, 2005, 2010)和大量镁铁-超镁铁质岩石(张作衡等, 2005; Tang et al., 2017; Gao and Zhou, 2013; Liu et al., 2017),并且这些岩石均具有正的高εNd(t)和εHf(t)值(王涛等, 2005, 2010; 童英等,2010; Gao and Zhou, 2013; Tang et al., 2017)。后碰撞时期是板块构造活动中的一个重要阶段,这一时期的俯冲板块拆离、岩石圈拆沉、伸展跨塌及岩浆底侵作用可导致大量后碰撞花岗质岩石形成(杨坤光和杨巍然, 1997; 夏祖春等, 2005)。因此,巴斯铁列克钨多金属矿区出露的花岗质岩体是后碰撞伸展构造环境的产物。

通常,高硅富碱过铝质花岗岩可以由地壳物质部分熔融形成(Kapp et al., 2002),也可以由壳幔岩浆混合形成(Korhonen et al., 2015; 张永明等, 2019)。巴斯铁列克矿区花岗质岩石中的锆石εHf(t)值分别为+5.6~+10、+4.2~+11和+3.8~+11,这些锆石的两阶段Hf同位素模式年龄(0.6~1.0Ga)大于锆石结晶年龄(284Ma),表明其成岩过程有多种组分参与。这些变化大的εHf(t)值介于亏损地幔和球粒陨石演化线之间,远高于古老地壳Hf演化线,暗示可能来源于新增生地壳或亏损地幔。花岗质岩浆不可能直接由地幔部分熔融形成,并且大量研究成果表明,阿尔泰造山带是典型的增生造山带,存在着一定规模的显生宙陆壳生长(王涛等, 2005),发育大量古生代具有正的高εNd(t)和εHf(t)值的花岗岩(王涛等, 2005, 2010; 童英等, 2010; Gao and Zhou, 2013; Tang et al., 2017),被认为是陆壳生长的结果(王涛等, 2005; Tang et al., 2017),同时这些花岗岩中有捕获锆石存在,岩石中有一定量的Cr和Ni等元素,因此这些花岗岩很可能来源于新生地壳物质与幔源岩浆,在上升过程中混入了少量古老地壳物质。二叠纪时期阿尔泰造山带南缘处于后碰撞伸展环境,伴随同时期大规模基性岩(Gao and Zhou, 2013; Tang et al., 2017),说明幔源岩浆为巴斯铁列克矿区花岗质岩石的形成不仅提供了重要的热源,也伴随有少量混入。考虑到岩石中有古老锆石颗粒和强烈的分离结晶,说明存在深部岩浆房。

综上所述,巴斯铁列克钨多金属矿区花岗质岩石形成于后碰撞伸展环境,底侵的幔源岩浆加热新生成地壳并导致部分熔融和混合,产生的混合岩浆在深部岩浆房发生强烈的结晶分异,最终上升侵位形成矿区各类岩石,表现出高钾钙碱性等壳源岩石特点,这些岩石类型的不同与母岩浆成分和结晶分异程度有关。

5.4 矿区花岗岩对钨矿化的启示

本研究获得的巴斯铁列克钨多金属矿区出露的二云母花岗岩、黑云母花岗岩和二长花岗岩的成岩年龄为284~282Ma,与矿区辉钼矿Re-Os同位素年龄结果(284.4Ma,杨富全等, 2019)一致,表明这些岩体在时空上与钨多金属矿化密切相关。

大量研究资料表明,与钨铜多金属矿化有关的花岗质岩体可以是S型也可以是I型(Soloviev and Kryazhev, 2018; Mao et al., 2019),含矿岩石主要为高钾钙碱性或碱性系列花岗闪长岩、二长花岗岩、二云母花岗岩和黑云母花岗岩(Hu et al., 2012; Mao et al., 2019),母岩浆性质可以是氧化的也可以是还原的,多数学者认为,富挥发分(富F)高分异氧化性岩浆有利于钨矿化。实验岩石学研究结果表明,钨在岩浆岩中表现为强不相容性(Adam and Green, 2006; Fulmer et al., 2010),在部分熔融过程中易于进入熔体,在结晶分异中易残留于熔体(Yuan et al., 2018, 2019)。钨也是亲氧元素,在岩浆岩中以[WO4]2-形式存在(Che et al., 2013; Wade et al., 2013),WO3在700~800℃花岗质熔体中的饱和溶解度超过1000×10-6(Štemprok, 1990),在岩浆结晶过程中很少形成钨的独立矿物,仅在钨含量较高的花岗岩中出现少量黑钨矿或白钨矿(Che et al., 2013; 祝新友等, 2016),主要以类质同象方式进入白云母、锂云母、黑云母和斜长石晶格中。但是W在岩浆熔体中的饱和溶解度与熔体成分强烈相关(Štemprok, 1990),与岩浆熔体的氧化还原状态相关性弱,富挥发份亚铝质或过铝质岩浆在低温时可以直接结晶出黑钨矿(Che et al., 2013; 丁腾等, 2017)。此外,岩浆挥发分的存在(H2O、HCl、HF、CO2、F)可以增加W的溶解度(Wood and Samson, 2000),尤其是F可以大幅降低熔体的粘度、密度和固相线温度,改变熔体结构(Agangi et al., 2012),使W在残余熔体中富集(Audétat et al., 2000)。这些成果表明花岗岩的分异演化过程也是钨的高度富集过程,分异花岗岩中高的金属元素W(Sn)或者由于岩浆结晶分异导致或者与岩浆-热液流体作用有关(毛景文等, 2008; Yuan et al., 2011)。

巴斯铁列克钨多金属矿区的花岗岩中钨的含量高,黑云母花岗岩为15.1×10-6~168×10-6,二长花岗岩为8.4×10-6~16×10-6,二云母花岗岩为4.6×10-6~9.4×10-6,均明显高于钨在下地壳(0.6×10-6)、上地壳(1.9×10-6)(Rudnick and Gao, 2004)和低钙型花岗岩中的平均含量(2×10-6Turekian and Wedepohl, 1961),可以为成矿作用提供丰富的物质来源。前已述及,巴斯铁列克钨多金属矿区花岗质岩体属分异花岗岩,母岩浆来源于新生地壳和地幔物质部分熔融形成的混合熔体,经过了一定程度的演化。黑云母花岗岩的钨含量明显高于二长花岗岩的钨含量,造岩矿物黑云母中的F和钨含量明显高于二长花岗岩中黑云母的F和钨的含量(未发表资料),黑云母花岗岩中出现白钨矿,与我国华南地区与钨矿化相关的多个岩体中(新田岭岩体, 章荣清, 2014;黄沙坪岩体, 丁腾等, 2017;瑶岗仙岩体, 董少花等, 2014)出现的白钨矿类似,被造岩矿物包裹也与未蚀变黑云母共生。尽管黑云母花岗岩中钨含量变化大(15.1×10-6~168×10-6),并且花岗岩中钨矿物常被认为是岩浆期后热液交代早期结晶花岗岩有关(Kempe et al., 1999)。以上这些特征为岩浆结晶分异对钨富集提供了重要证据。

前人研究结果表明,碱性长石的结晶分异可以导致Ba、Sr和Eu的亏损,但岩浆分异晚期流体作用可以使这些元素亏损更明显(Jahn et al., 2001);在花岗质岩浆发生高度分异演化形成流体时,由于Ta较Nb难进入流体,导致花岗岩体系中Nb/Ta比值显著降低(Zarasky et al., 2010)。钨在相对还原的条件下表现为亲铁性,并且随着还原作用越强,钨更容易从岩浆迁移到热液中(Sylvester, 1998),使钨含量迅速增加。二云母花岗岩具有深“V”型稀土配分模式,较黑云母花岗岩和二长花岗岩有更加亏损的Ba、Sr和Eu,有极低的Nb/Ta比值(4.3~5.6),采用Qiu et al. (2013)的方法计算二云母花岗岩中锆石的氧逸度(logfo2)值低(-24.6~-14.4),岩石中黑电气石的存在可以延长岩浆房出溶流体的时间并降低出溶流体的最低温度(Chang and Meinert, 2008)。以上特征应该与岩浆结晶晚期发生了熔体-流体强烈相互作用有关。

综上所述,巴斯铁列克钨多金属矿区二云母花岗岩、黑云母花岗岩和二长花岗岩是原生岩浆经过高程度演化形成高硅富碱、富含挥发分并富集成矿金属元素的母岩浆结晶形成。花岗质岩浆为钨矿床的形成提供了物质基础,岩浆演化(结晶分异作用和晚期熔体-流体的过程)对钨有初始富集作用。

前已述及,巴斯铁列克钨多金属矿区花岗质岩石具有I-A过渡型特征,锆石εHf(t)值均为高的正值,我国钨锡成矿作用极为发育的华南地区与钨矿化相关花岗岩多为S型(黄兰椿和蒋少勇等, 2013; 原垭斌等, 2014; Zhao et al., 2016; Yuan et al., 2018, 2019),锆石εHf(t)值均为高的负值;巴斯铁列克含矿岩石母岩浆来源于新生地壳和幔源岩浆,华南地区含矿岩体母岩浆主要来源于古老地壳(黄兰椿和蒋少勇等, 2013; 原垭斌等, 2014; Zhao et al., 2016),这些原生岩浆有高的钨初始含量,高度分异和熔体-流体出溶等对钨的富集(黄兰椿和蒋少勇等, 2013; 郭春丽等, 2017; 陈伟等, 2018)是两地区岩石的共同特征。由此可见,原生岩浆的分异演化过程是导致该类型矿床成矿的关键控制因素之一。

6 结论

(1) 巴斯铁列克钨多金属矿区二云母花岗岩、黑云母花岗岩和二长花岗岩结晶锆石的加权平均年龄分别为282.3±3.2Ma、284.3±2.2Ma和284.8±2.3Ma,是早二叠世岩浆活动的产物,与矿区辉钼矿成矿年龄一致。

(2) 二云母花岗岩、黑云母花岗岩和二长花岗岩是高钾钙碱性过铝质系列岩石,属分异I-A过渡型花岗岩。具有富集Th、U、K、Pb,亏损Ba、Sr、P、Ti和高的正εHf(t)值,是由新生地壳熔体和幔源岩浆混合形成的熔体经分异演化形成。它们代表了同一时期两个不同的岩浆事件。黑云母花岗岩和二长花岗岩是同一源区不同演化阶段的产物,二云母花岗岩来自不同源区且有更高分异演化程度。

(3) 二云母花岗岩、黑云母花岗岩和二长花岗岩是巴斯铁列克钨多金属矿床的成矿母岩,岩浆的演化分异(结晶分异与熔体-流体)对成矿元素有一定的初始富集作用,为矿床的形成提供了成矿元素和成矿流体。

致谢      野外工作期间得到新疆维吾尔自治区有色地质勘查局706队郭旭吉、王永强等领导的支持,中国地质科学院矿产资源研究所的杨富全研究员、李强和杨成栋副研究员的指导和帮助;项目实施过程中得到了新疆维吾尔自治区有色地质勘查局地质矿产勘查研究院领导和技术人员的帮助;样品的测试工作在实验人员的帮助下完成;审稿专家提出了宝贵的修改意见。在此一并表示感谢。

参考文献
Adam J and Green T. 2006. Trace element partitioning between mica-and amphibole-bearing garnet lherzolite and hydrous basanitic melt: 1. Experimental results and the investigation of controls on partitioning behaviour. Contributions to Mineralogy and Petrology, 152(1): 1-17
Agangi A, Kamenetsky VS and McPhie J. 2012. Evolution and emplacement of high fluorine rhyolites in the Mesoproterozoic Gawler silicic large igneous province, South Australia. Precambrian Research, 208-211: 124-144 DOI:10.1016/j.precamres.2012.03.011
Audťat A, D and Heinrich CA. 2000. Magmatic-hydrothermal evolution in A fractionating granite: A microchemical study of the Sn-W-F-mineralized Mole granite (Australia). Geochimica Et Cosmochimica Acta, 64(19): 3373-3393 DOI:10.1016/S0016-7037(00)00428-2
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 DOI:10.1130/G37475.1
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
Cai KD, Sun M, Yuan C, Long XP and Xiao WJ. 2011. Geological framework and Paleozoic tectonic history of the Chinese Altai, NW China: A review. Russian Geology and Geophysics, 52(12): 1619-1633 DOI:10.1016/j.rgg.2011.11.014
Chai FM, Dong LH, Yang FQ, Liu F, Geng XX and Huang CK. 2010. Age, geochemistry and petrogenesis of Tiemierte granites in the Kelang basin at the southern margin of Altay, Xinjiang. Acta petrologica Sinica, 26(2): 377-386 (in Chinese with English abstract)
Chai FM, Yang FQ, Liu F, Geng XX, Jiang LP, Lü SJ, Guo XJ and Chen B. 2012. Geochronology and genesis of meta-felsic volcanic rocks from the Kangbutiebao Formation in Chonghuer basin on southern margin of Altay, Xinjiang. Geological review, 58(6): 1023-1037 (in Chinese with English abstract)
Chappell BW. 1999. Aluminium saturation in I- and S-type granites and the characterization of fractionated haplogranites. Lithos, 46(3): 535-551 DOI:10.1016/S0024-4937(98)00086-3
Chang ZS and Meinert LD. 2008. The Empire Cu-Zn mine, Idaho: Exploration implications of unusual skarn features related to high fluorine activity. Economic Geology, 103(5): 909-938 DOI:10.2113/gsecongeo.103.5.909
Che XD, Linnen RL, Wang RC, Aseri A and Thibault Y. 2013. Tungsten solubility in evolved granitic melts: An evaluation of magmatic wolframite. Geochimica et Cosmochimica Acta, 106: 84-98 DOI:10.1016/j.gca.2012.12.007
Chen W, Chen B and Sun KK. 2018. Petrogenesis of the Maogongdong highly differentiated granite in the Dahutang tungstenore field, Jiangxi Province. Acta Petrologica Sinica, 34(6): 1704-1724 (in Chinese with English abstract)
Ding T, Ma DS, Lu JJ, Zhang RQ and Zhang ST. 2017. Mineral geochemistry of granite porphyry in Huangshaping pollymetallic deposit, southern Hunan Province, and its implications for metallogensis of skarn scheelite mineralization. Acta Petrologica Sinica, 33(3): 716-728 (in Chinese with English abstract)
Dong SH, Bi XW, Hu RZ and Chen YW. 2014. Petrogenesis of the Yaogangxian granites and implications for W mineralization, Hunan Province. Acta Petrologica Sinica, 30(9): 2749-2765 (in Chinese with English abstract)
Eby GN. 1990. The A-type granitoids: A review of their occurrence and chemical characteristics and speculations on their petrogenesis. Lithos, 26(1-2): 115-134 DOI:10.1016/0024-4937(90)90043-Z
Fulmer EC, Nebel O and Westrenen WV. 2010. High-precision high field strength element partitioning between garnet, amphibole and alkaline melt from Kakanui, New Zealand. Geochimica et Cosmochimica Acta, 74(9): 2741-2759 DOI:10.1016/j.gca.2010.02.020
Gao JF and Zhou MF. 2013. Magma mixing in the genesis of the Kalatongke dioritic intrusion: Implications for the tectonic switch from subduction to post-collision, Chinese Altay, NW China. Lithos, 162-163: 236-250 DOI:10.1016/j.lithos.2013.01.007
Griffin WL, Wang X, Jackson SE, Pearson NJ, O'Reilly SY, Xu XS and Zhou XM. 2002. Zircon chemistry and magma mixing, SE China: In-situ analysis of Hf isotopes, Tonglu and Pingtan igneous complexes. Lithos, 61(3-4): 237-269 DOI:10.1016/S0024-4937(02)00082-8
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)
Hu RZ and Zhou M. 2012. Multiple Mesozoic mineralization events in South China: An introduction to the thematic issue. Mineralium Deposita, 47(6): 579-588 DOI:10.1007/s00126-012-0431-6
Huang LC and Jiang SY. 2013. Geochronology, geochemistry and petrogenesis of the tungsten bearing porphyritic granite in the Dahutang tungsten deposit, Jiangxi Province. Acta Geologica Sinica, 29(12): 4323-4335 (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
Jahn BM, Wu FY, Capdevila R, Martineau F, Zhao ZH and Wang YX. 2001. Highly evolved juvenile granites with tetrad REE patterns: the Woduhe and Baerzhe granites from the Great Xing'an Mountains in NE China. Lithos, 59(4): 171-198 DOI:10.1016/S0024-4937(01)00066-4
Kapp JDA, Miller CF and Miller JS. 2002. Ireteba pluton, Eldorado Mountains, Nevada: Late, deep-source, peraluminous magmatism in the Cordilleran Interior. The Journal of Geology, 110(6): 649-669 DOI:10.1086/342864
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
Korhonen F, Brown M, Clark C, Foden JD and Taylor R. 2015. Are granites and granulites consanguineous?. Geology, 43(11): 991-994 DOI:10.1130/G37164.1
Li Q, Yang FQ and Yang CD. 2019. Geochronology and geochemical characteristics of Dahalasu granite in Altay of Xinjiang, China and their tectonic significance. Journal of Earch Sciences and Environment, 41(4): 396-413 (in Chinese with English abstract)
Li SC, Li YF, Wang XA and Zhu K. 2016. Delineation of the Late Triassic granitic pluton from the middle part of Greater Xing'an Mountains showing tetrad REE patterns and its geological implications. Acta Petrologica Sinica, 32(9): 2793-2806 (in Chinese with English abstract)
Li XR, Liu F and Yang FQ. 2012. Geological times and its signification of the two mica syenogranite in the Keyinblak Cu-Zn deposit area in Altay, Xinjiang. Xinjiang Geology, 30(1): 5-11 (in Chinese with English abstract)
Li Y. 2018. Geological characteristics and prospecting significance of Basitielieke scheelite deposit in southern margin of the Altai, Xinjiang. Xinjiang Nonferrous Metals, 41(3): 14-15 (in Chinese)
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 CS, Chen XM, Chen PR and Hu H. 2003. Subdivision, discrimination criteria and genesis for A-type rock suites. Geological Journal of China Universities, 9(4): 573-591 (in Chinese with English abstract)
Liu F, Yang FQ, Mao JW, Chai FM and Geng XX. 2009. Study on chronology and geochemistry for Abagong granite in Altay orogen. Acta Petrologica Sinica, 25(6): 1416-1425 (in Chinese with English abstract)
Liu F, Zhang C and Yang FQ. 2012. A study of chronology and ore-forming process of Jialbasto iron deposit on southern margin of Altay. Mineral Deposits, 31(6): 1277-1288 (in Chinese with English abstract)
Liu W, Liu LJ, Liu XJ, Shang HJ and Zhou G. 2010. Age of the Early Devonian Kangbutiebao formation along the southern Altay Mountains and its northeastern extension. Acta Petrologica Sinica, 26(2): 387-400 (in Chinese with English abstract)
Liu XJ, Liu W and Si CQ. 2019. Petrogenesis and source rocks of the high-K calc-alkaline and shoshonitic I type granitoids in the northwestern part of East Junggar, NW China. Lithos, 326-327: 298-312 DOI:10.1016/j.lithos.2018.12.033
Liu YL, Zhang H, Tang Y, Zhang X, Lv ZH and Zhao JY. 2017. Petrogenesis and tectonic setting of the Middle Permian A-type granites in Altay, northwestern China: Evidences from geochronological, geochemical, and Hf isotopic studies. Geological Journal, 53(2): 527-546
Liu Z and Tong LX. 2015. The petrogenetic link between the Late Paleozoic high-temperature metamorphism in the Altay orogen and the Tarim mantle plume: Evidence from metapelitic and mafic granulites. Acta Petrologica Sinica, 31(6): 1761-1773 (in Chinese with English abstract)
Long XP, Sun M, Yuan C, Xiao WJ, Chen HL, Zhao YJ, Cai KD and Li JL. 2006. Genesis of Carboniferous volcanic rocks in the eastern Junggar: Constraints on the closure of the Junggar Ocean. Acta Petrologica Sinica, 22(1): 31-40 (in Chinese with English abstract)
Maniar PD and Piccoli PM. 1989. Tectonic discrimination of granitoids. Geological Society of American Bulletin, 101: 635-643 DOI:10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
Mao JW, Xie GQ, Guo CL, Yuan SD, Cheng YB and Chen YC. 2008. Spatial-temporal distribution of mesozoic ore deposits in South China and their metallogenic settings. Geological Journal of China Universities, 14(4): 510-526 (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. SEG Special Publications, 22: 411-482
Miller CF, McDowell SM and Mapes RW. 2003. Hot and cold granite? Implications of zircon saturation temperatures and preservation of inheritance. Geology, 31(6): 529-532 DOI:10.1130/0091-7613(2003)031<0529:HACGIO>2.0.CO;2
Peccerillo R 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
Pirajno F, Mao JW, Zhang ZC, Zhang ZH and Chai FM. 2008. The association of mafic-ultramafic intrusions and A-type magmatism in the Tian Shan and Altay orogens, NW China: Implications for geodynamic evolution and potential for the discovery of new ore deposits. Journal of Asian Earth Sciences, 32(2-4): 165-183 DOI:10.1016/j.jseaes.2007.10.012
Qiu JT, Yu XQ, Santosh M, Zhang DH, Chen SQ and Li PJ. 2013. Geochronology and magmatic oxygen fugacity of the Tongcun molybdenum deposit, Northwest Zhejiang, SE China. Mineralium Deposita, 48(5): 545-556 DOI:10.1007/s00126-013-0456-5
Ren BQ, Zhang H, Tang Y and Lü ZH. 2011. LA-ICP-MS U-Pb zircon geochronology of the Altai pegmatites and its geological significance. Acta Mineralogica Sinica, 31(3): 587-596 (in Chinese with English abstract)
Rudnick RL and Gao S. 2004. Composition of the continental crust. Treatise on Geochemistry, 3: 1-64
Shan Q, Zeng QS, Luo Y, Yang WB, Zhang H, Qiu YZ and Yu XY. 2011. SHRIMP U-Pb ages and petrology studies on the potassic and podic rhyolites in Altai, North Xinjiang. Acta Petrologica Sinica, 27(12): 3653-3665 (in Chinese with English abstract)
Shen XM, Zhang HX, Wang Q, Wyman DA and Yang YH. 2011. Late Devonian-Early Permian A-type granites in the southern Altay Range, Northwest China: Petrogenesis and implications for tectonic setting of "A2-type" granites. Journal of Asian Earth Sciences, 42(5): 986-1007 DOI:10.1016/j.jseaes.2010.10.004
Shen XM, Zhang HX and Ma L. 2013. Discovery of the Late Carboniferous leucogranite in the southern Altay Range and its tectonic implications. Geotectonica et Metallogenia, 37(139): 721-729 (in Chinese with English abstract)
Söderlund Ulf, 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
Soloviev SG and Kryazhev SG. 2018. Geology, mineralization, and fluid inclusion characteristics of the Kashkasu W-Mo-Cu skarn deposit associated with a high-potassic to shoshonitic igneous suite in Kyrgyzstan, Tien Shan: Toward a diversity of W mineralization in Central Asia. Journal of Asian earth sciences, 153(3): 425-449
Štemprok M. 1990. Solubility of tin, tungsten and molybdenum oxides in felsic magmas. Mineralium Deposita, 25(3): 205-212 DOI:10.1007/BF00190382
Su HM and Jiang SY. 2017. A comparison study of tungsten-bearing granite and related mineralization in the northern Jiangxi-southern Anhui provinces and southern Jiangxi Province in South China. Science China (Earth Sciences), 60(11): 1942-1958 DOI:10.1007/s11430-016-9071-6
Sun GH, Li JY, Yang TN, Li YP, Zhu ZX and Yang ZQ. 2009. Zircon SHRIMP U-Pb dating of two linear granite plutons in southern Altay Mountains and its tectonic implications. Geology in China, 36(05): 976-987 (in Chinese with English abstract)
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 Special Publication, 42(1): 313-345
Sylvester PJ. 1998. Post-collisoonal strongly peraluminius granites. Lithos, 45(1-4): 29-44 DOI:10.1016/S0024-4937(98)00024-3
Tang DM, Qin KZ, Xue SC, Mao YJ, Evans NJ, Niu YJ and Chen JL. 2017. Genesis of the Permian Kemozibayi sulfide-bearing mafic-ultramafic intrusion in Altay, NW China: Evidence from zircon geochronology, Hf and O isotopes and mineral chemistry. Lithos, 292-293: 49-68 DOI:10.1016/j.lithos.2017.08.021
Tong LX, Chen YB and Chen LL. 2014. Confirmation of ultrahigh-temperature metapelitic granulite in the Altay orogen and its geological significance. Chinese Science Bulletin, 59(20): 1984-1991 (in Chinese) DOI:10.1360/csb2014-59-20-1984
Tong Y, Wang T, Hong DW, Han BF, Zhang JJ, Shi XJ and Wang C. 2010. Spatial and temporal distribution of the Carboniferous-Permian granitoids in northern Xinjiang and its adjacent areas, and its tectonic significance. Acta Petrologica et Mineralogica, 29(6): 619-641 (in Chinese with English abstract)
Tong Y, Wang T, Jahn BM, Sun M, Hong DW and Gao JF. 2014. Post-accretionary Permian granitoids in the Chinese Altai orogen: Geochronology, petrogenesis and tectonic implications. American Journal of Science, 314(1): 80-109 DOI:10.2475/01.2014.03
Turekian KK and Wedepohl KH. 1961. Distribution of the elements in some major units of the earth's crust. Geological Society of America Bulletin, 72(2): 175-192 DOI:10.1130/0016-7606(1961)72[175:DOTEIS]2.0.CO;2
Vervoot JD, Pachelt PJ, Gehrels GE and Nutman AP. 1996. Constraints on early earth differentiation from hafnium and neodymium isotopes. Nature, 379(6566): 624-627 DOI:10.1038/379624a0
Vervoot JD and Blichert-Toft J. 1999. Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochimica et Cosmochimica Acta, 63: 533-556 DOI:10.1016/S0016-7037(98)00274-9
Wade J, Wood BJ and Norris CA. 2013. The oxidation state of tungsten in silicate melt at high pressures and temperatures. Chemical Geology, 335: 189-193 DOI:10.1016/j.chemgeo.2012.10.011
Wang JB and Xu X. 2006. Post-collisional tectonic evolution and metallogenesis in northern Xinjiang, China. Acta Geologica Sinica, 80(1): 23-31 (in Chinese with English abstract)
Wang T, Hong DW, Tong Y, Han BF and Shi YR. 2005. Zircon U-Pb SHRIMP age and origin of post-orogenic Lamazhao granitic pluton from Altai orogen: Its implications for vertical continental growth. Acta Petrologica Sinica, 21(3): 640-650 (in Chinese with English abstract)
Wang T, Hong DW, Jahn BM, Tong Y, Wang YB, Han BF and Wang XX. 2006. Timing, petrogenesis, and setting of Paleozoic synorogenic intrusions from the Altai Mountains, Northwest China: Implications for the tectonic evolution of an accretionary orogen. The Journal of Geology, 114: 735-751 DOI:10.1086/507617
Wang T, Tong Y, Zhang JJ, Li S, Shi XJ, Li JY, Han BF and Hong DW. 2010. Spatial and temporal variations of granitoids in the Altay orogen and their implications for tectonic setting and crustal growth: Perspectives from Chinese Altay. Acta Petrologica et Mineralogica, 29(6): 595-618 (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
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
Windley BF, Kröener A, Guo JH, Qu GS, Li YY and Zhang C. 2002. Neoproterozoic to Paleozoic geology of the Altai Orogen, NW China: New zircon age data and tectonic evolution. The Journal of Geology, 110(6): 719-737 DOI:10.1086/342866
Wood SA and Samson IM. 2000. The hydrothermal geochemistry of tungsten in granitoid environments: Ⅰ. relative solubilities of ferberite and scheelite as a function of T, P, pH, and mNaCl. Economic Geology, 95(1): 143-182
Wu FY, Jahn BM, Wilde SA, Lo CH, Yui TF, Lin Q, Ge WC and Sun DY. 2003. Highly fractionated I-type granites in NE China (Ⅱ): Isotopic geochemistry and implications for crustal growth in the Phanerozoic. Lithos, 67: 191-204 DOI:10.1016/S0024-4937(03)00015-X
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, 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
Xia ZC, Xu XY, Xia LQ, Li XM, Ma ZP and Wang LS. 2005. Geochemistry of the Carboniferous-Permian post-collisional granitic rocks from Tianshan. Northwestern Geology, 38(1): 1-14 (in Chinese with English abstract)
Xiao WJ, Windley BF, Badarch G, Sun SL, Li JY, Qin KZ and Wang ZH. 2004. Palaeozoic accretionary and convergent tectonics of the southern Altaids: Implications for the growth of Central Asia. Journal of the Geological Society, 161(3): 339-342 DOI:10.1144/0016-764903-165
Xiao WJ, Han CM, Yuan C, Chen HL, Sun M, Lin SF, Li ZL, Mao QG, Zhang JE, Sun S and Li JL. 2006. Unique Carboniferous-Permian tectonic-metallogenic framework of Northern Xinjiang (NW China): Constraints for the tectonics of the southern Paleoasian domain. Acta Petrologica Sinica, 22(5): 1062-1076 (in Chinese with English abstract)
Xiao WJ, Windley BF and Yuan C. 2009. Paleozoic multiple subduction-accretion processes of the Southern Altaids. American Journal of Science, 309: 221-270 DOI:10.2475/03.2009.02
Yan SH, Chen W, Wang YT, Zhang ZC and Chen BL. 2004. 40Ar/39Ar dating and its significance of the Ertix gold metallogenic belt in the Altay Oregon, Xinjiang. Acta Geologica Sinica, 78(4): 500-506 (in Chinese with English abstract)
Yang FQ, Mao JW, Liu F, Dong YG, Chai FM and Geng XX. 2010. Geological characteristics and metallogenesis of the Keyinbulake copper-zinc deposit in Altay, Xinjiang. Acta Petrologica Sinica, 26(2): 361-376 (in Chinese with English abstract)
Yang FQ, Chai FM, Li Q, Gu GL and Yang CD. 2017. The discovery of the Late Cambrian-Early Ordovician strata in the basement of the Kelan basin, Altay, Xinjiang: Constrains from zircon LA-MC-ICP-MS U-Pb dating. Acta Geoscientica Sinica, 38(6): 859-871 (in Chinese with English abstract)
Yang FQ, Wang YQ, Yang CD, Li Q and Guo XJ. 2018. Metallic ore deposits in the Kelan Basin, Altay, Xinjiang: Geological characteristics and metallogenesis. Mineral Deposits, 37(3): 441-462 (in Chinese with English abstract)
Yang FQ, Li Y, Chen GM, Zhang GF, Li C, Tian M, Li P and Li Q. 2019. Ore-forming age of the Basitielieke tungsten polymetallic deposit in the Altay, Xinjiang: Constraits of zircon U-Pb and molybdenite Re-Os geochronology. Geotectonica et Metallogenia, 43(6): 1144-1154 (in Chinese with English abstract)
Yang FQ, Zhang B, Yang CD, Li Q and Wang YQ. 2020. Geology and geochronology of the volcanogenic massive sulphide polymetallic deposits in Altay Orogenic Belt, Xinjiang, Northwest China: Examples from the Kelan Basin. International Geology Review, 1-16
Yang KG and Yang WR. 1997. Post-collision orogeny process and origin of huge quantity granites in orogenic belt. Geological Science and Technology Information, 16(4): 16-22 (in Chinese with English abstract)
Yuan C, Sun M, Xiao WJ, Li XH, Chen HL, Lin SF, Xia XP and Long XP. 2007. Accretionary orogenesis of the Chinese Altai: Insights from Paleozoic granitoids. Chemical Geology, 242(1-2): 22-39 DOI:10.1016/j.chemgeo.2007.02.013
Yuan SD, Peng JT, Hao S, Li HM, Geng JZ and Zhang DL. 2011. In situ LA-MC-ICP-MS and ID-TIMS U-Pb geochronology of cassiterite in the giant Furong tin deposit, Hunan Province, South China: New constraints on the timing of tin-polymetallic mineralization. Ore Geology Reviews, 43(1): 235-242 DOI:10.1016/j.oregeorev.2011.08.002
Yuan SD, Williams-Jones AE, Mao JW, Zhao PL, Yan C and Zhang DL. 2018. The origin of the Zhangjialong tungsten deposit, South China: Implications for W-Sn mineralization in large granite batholiths. . Economic Geology, 113: 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: 1005-1012 DOI:10.5382/econgeo.4669
Yuan YB, Yuan SD, Chen CJ and Huo R. 2014. Zircon U-Pb ages and Hf isotopes of the granitoids in the Huangshaping mining area and their geological significance. Acta Petrologica Sinica, 30(1): 64-78 (in Chinese with English abstract)
Zarasky GP, Korzhinskaya V and Kotova N. 2010. Experimental studies of Ta2O5 and columbite-tantalite solubility in Fluoride solutions from 300 to 550℃ and 50 to 100MPa. Mineralogy and Petrology, 99(3/4): 287-300
Zen EA. 1986. Aluminum enrichment in silicate melts by fractional crystallization: Some mineralogic and petrographic constraints. Journal of Petrology, 27(5): 1095-1117 DOI:10.1093/petrology/27.5.1095
Zhang C, Liu F and Yang FQ. 2013. A discussion on genetic mechanism of the Kuerqis iron deposit in Fuyun County, Xinjiang. Rock and Mineral Analysis, 32(1): 157-165 (in Chinese with English abstract)
Zhang C, Liu DD, Luo Q, Liu LF, Zhang YZ, Zhu DY, Wang PF and Dai QQ. 2018. An evolving tectonic environment of Late Carboniferous to Early Permian granitic plutons in the Chinese Altai and Eastern Junggar terranes, Central Asian Orogenic Belt, NW China. Journal of Asian Earth Sciences, 159: 185-208 DOI:10.1016/j.jseaes.2017.08.008
Zhang GF, Li Y, Yang FQ, Li Q, Yang CD and Chai FM. 2019. Mineral characteristics of Basitielieke W-polymetallic deposit in Altay of Xinjiang and their geological significance. Mineral Deposits, 38(1): 101-119 (in Chinese with English abstract)
Zhang RQ. 2014. Petrogenesis and metallogeny of the W and Sn-bearing granites in southern Hunan province: case study from wangxianling and Xintianling. Ph. D. Dissertation. Nanjing: Nanjing University (in Chinese with English summary)
Zhang YM, Pei XZ, Li ZC, Li RB, Liu CJ, Pei L, Chen YX and Wang M. 2019. Zircon U-Pb age, geochemical characteristics and geological significance of the Caledonian strongly peraluminous granites in the Nanshan area, Qinghai Province. Geological Bulletin of China, 38(5): 742-756 (in Chinese with English abstract)
Zhang ZH, Chai FM, Du AD, Zhang ZC, Yan SH, Yang JM, Qu WJ and Wang ZL. 2005. Re-Os dating and ore-forming material tracing of the Karatungk Cu-Ni sulfide deposit in northern Xinjiang. Acta Petrologica et Mineralogica, 24(4): 285-293 (in Chinese with English abstract)
Zhang ZX, Chai FM, Wang LF, Liu F and Geng XX. 2014. Geochemical characteristics and its significance of biotite monzogranite in the Keyinbulake copper-zinc deposit at the southern margin of Altay. Xinjiang. Geoscience, 28(3): 461-471 (in Chinese with English abstract)
Zhao PL, Yuan SD, Mao JW, Santosh M, Li C and Hou KJ. 2016. Geochronological and petrogeochemical constraints on the skarn deposits in Tongshanling ore district, southern Hunan Province: Implications for Jurassic Cu and W metallogenic events in South China. Ore Geology Reviews, 78: 120-137 DOI:10.1016/j.oregeorev.2016.03.004
Zhao YM, Peng G, Zhao DL and You YC. 2016. LA-ICP-MS zircon U-Pb age and geochemistry of the Kezile granitic pluton in Altay of China and their geological implications. Geology and Exploration, 52(2): 271-282 (in Chinese with English abstract)
Zhou G, Zhang ZC, Gao GZ, Yang WP, He B, Zhang XL, Luo SB, Wang X and He YK. 2006. Dating of granitic plutons and its geological implications in the lower reaches of the Qinggeli River in the northern part of East Junggar, Xinjiang. Geoscience, 20(1): 141-150 (in Chinese with English abstract)
Zhu XY, Wang YL, Cheng XY and Tian Y. 2016. Geological and geochemical characteristics of xenoliths in Yaogangxian granite, Hunan Province. Acta Petrologica et Mineralogica, 35(1): 16-32 (in Chinese with English abstract)
柴凤梅, 董连慧, 杨富全, 刘锋, 耿新霞, 黄承科. 2010. 阿尔泰南缘克朗盆地铁木尔特花岗岩体年龄、地球化学特征及成因. 岩石学报, 26(2): 377-386.
柴凤梅, 杨富全, 刘锋, 耿新霞, 姜丽萍, 吕书君, 郭旭吉, 陈斌. 2012. 阿尔泰南缘冲乎尔盆地康布铁堡组变质酸性火山岩年龄及岩石成因. 地质论评, 58(6): 1023-1037. DOI:10.3969/j.issn.0371-5736.2012.06.003
陈伟, 陈斌, 孙克克. 2018. 江西大湖塘矿集区茅公洞钨矿有关的高分异花岗岩成因研究. 岩石学报, 34(6): 1704-1724.
丁腾, 马东升, 陆建军, 章荣清, 张世涛. 2017. 湘南黄沙坪多金属矿床花岗斑岩的矿物化学及其对矽卡岩白钨矿成矿的指示意义. 岩石学报, 33(3): 716-768.
董少花, 毕献武, 胡瑞忠, 陈佑纬. 2014. 湖南瑶岗仙复式花岗岩岩石成因及与钨成矿关系. 岩石学报, 30(9): 2749-2765.
郭春丽, 曾令森, 高利娥, 苏红中, 马星华, 尹冰. 2017. 福建河田高分异花岗岩的矿物和全岩地球化学找矿标志研究. 地质学报, 91(8): 1796-1817. DOI:10.3969/j.issn.0001-5717.2017.08.010
黄兰椿, 蒋少涌. 2013. 江西大湖塘富钨花岗斑岩年代学、地球化学特征及成因研究. 岩石学报, 29(12): 4323-4335.
李强, 杨富全, 杨成栋. 2019. 新疆阿尔泰大喀拉苏花岗岩年代学、地球化学特征及其构造意义. 地球科学与环境学报, 41(4): 396-413. DOI:10.3969/j.issn.1672-6561.2019.04.002
李世超, 李永飞, 王兴安, 朱凯. 2016. 大兴安岭中段晚三叠世四分组效应花岗岩的厘定及其地质意义. 岩石学报, 32(9): 2793-2806.
李香仁, 刘锋, 杨富全. 2012. 阿尔泰克因布拉克铜锌矿区二云母正长花岗岩成岩时代及地质意义. 新疆地质, 30(1): 5-11. DOI:10.3969/j.issn.1000-8845.2012.01.002
李永. 2018. 新疆阿尔泰南缘巴斯铁列克白钨矿床地质特征及找矿意义. 新疆有色金属, 41(3): 14-15.
刘昌实, 陈小明, 陈培荣, 王汝成, 胡欢. 2003. A型岩套的分类、判别标志和成因. 高校地质学报, 9(4): 573-591. DOI:10.3969/j.issn.1006-7493.2003.04.011
刘锋, 杨富全, 毛景文, 柴凤梅, 耿新霞. 2009. 阿尔泰造山带阿巴宫花岗岩体年代学及地球化学研究. 岩石学报, 25(6): 1416-1425.
刘锋, 张超, 杨富全. 2012. 新疆阿尔泰南缘加尔巴斯套铁矿床成矿时代及成矿作用研究. 矿床地质, 31(6): 1277-1288. DOI:10.3969/j.issn.0258-7106.2012.06.012
刘伟, 刘丽娟, 刘秀金, 尚海军, 周刚. 2010. 阿尔泰南缘早泥盆世康布铁堡组的SIMS锆石U-Pb年龄及其向东向北延伸的范围. 岩石学报, 26(2): 387-400.
刘兆, 仝来喜. 2015. 阿尔泰造山带晚古生代高温变质作用与塔里木地幔柱活动的成因联系: 来自泥质和镁铁质麻粒岩的证据. 岩石学报, 31(6): 1761-1773.
龙晓平, 孙敏, 袁超, 肖文交, 陈汉林, 赵永久, 蔡克大, 李继亮. 2006. 东准噶尔石炭系火山岩的形成机制及其对准噶尔洋盆闭合时限的制约. 岩石学报, 22(1): 31-40.
毛景文, 谢桂青, 郭春丽, 袁顺达, 程彦博, 陈毓川. 2008. 华南地区中生代主要金属矿床时空分布规律和成矿环境. 高校地质学报, 14(4): 510-526. DOI:10.3969/j.issn.1006-7493.2008.04.005
任宝琴, 张辉, 唐勇, 吕正航. 2011. 阿尔泰造山带伟晶岩年代学及其地质意义. 矿物学报, 31(3): 587-596.
单强, 曾乔松, 罗勇, 杨武斌, 张红, 裘瑜卓, 于学元. 2011. 新疆阿尔泰康布铁堡组钾质和钠质流纹岩的成因及同位素年代学研究. 岩石学报, 27(12): 3653-3665.
沈晓明, 张海祥, 马林. 2013. 阿尔泰南缘晚石炭世淡色花岗岩的发现及其构造意义. 大地构造与成矿学, 37(139): 721-729.
苏慧敏, 蒋少涌. 2017. 赣南和赣北-皖南钨成矿带含钨花岗岩及其成矿作用特征对比研究. 中国科学(地球科学), 47(11): 1292-1308.
孙桂华, 李锦轶, 杨天南, 李亚萍, 朱志新, 杨之青. 2009. 阿尔泰山脉南部线性花岗岩锆石SHRIMP U-Pb定年及其地质意义. 中国地质, 36(05): 976-987. DOI:10.3969/j.issn.1000-3657.2009.05.003
仝来喜, 陈义兵, 陈林丽. 2014. 阿尔泰造山带超高温变泥质麻粒岩的确定及其地质意义. 科学通报, 59(20): 1984-1991.
童英, 王涛, 洪大卫, 韩宝福, 张建军, 史兴俊, 王超. 2010. 北疆及邻区石炭-二叠纪花岗岩时空分布特征及其构造意义. 岩石矿物学杂志, 29(6): 619-641. DOI:10.3969/j.issn.1000-6524.2010.06.003
王京彬, 徐新. 2006. 新疆北部后碰撞构造演化与成矿. 地质学报, 80(1): 23-31.
王涛, 洪大卫, 童英, 韩宝福, 石玉若. 2005. 中国阿尔泰造山带后造山喇嘛昭花岗岩体锆石SHRIMP年龄、成因及陆壳垂向生长意义. 岩石学报, 21(3): 640-650.
王涛, 童英, 李舢, 张建军, 史兴俊, 李锦轶, 韩宝福, 洪大卫. 2010. 阿尔泰造山带花岗岩时空演变、构造环境及地壳生长意义——以中国阿尔泰为例. 岩石矿物学杂志, 29(6): 595-618. DOI:10.3969/j.issn.1000-6524.2010.06.002
吴福元, 李献华, 杨进辉, 郑永飞. 2007. 花岗岩成因研究的若干问题. 岩石学报, 23(6): 1217-1238. DOI:10.3969/j.issn.1000-0569.2007.06.001
吴福元, 刘小驰, 纪伟强, 王佳敏, 杨雷. 2017. 高分异花岗岩的识别与研究. 中国科学(地球科学), 47(7): 745-765.
夏祖春, 徐学义, 夏林圻, 李向民, 马中平, 王立社. 2005. 天山石炭-二叠纪后碰撞花岗质岩石地球化学研究. 西北地质, 38(1): 1-14. DOI:10.3969/j.issn.1009-6248.2005.01.001
肖文交, 韩春明, 袁超, 陈汉林, 孙敏, 林寿发, 厉子龙, 毛启贵, 张继恩, 孙枢, 李继亮. 2006. 新疆北部石炭纪-二叠纪独特的构造-成矿作用: 对古亚洲洋构造域南部大地构造演化的制约. 岩石学报, 22(5): 1062-1076.
闫升好, 陈文, 王义天, 张招崇, 陈柏林. 2004. 新疆额尔齐斯金成矿带的40Ar/39Ar年龄及其地质意义. 地质学报, 78(4): 500-506.
杨富全, 毛景文, 刘锋, 董永观, 柴凤梅, 耿新霞. 2010. 新疆阿尔泰克因布拉克铜锌矿床地质特征及成矿作用. 岩石学报, 26(2): 361-376.
杨富全, 柴凤梅, 李强, 谷高灵, 杨成栋. 2017. 新疆阿尔泰克兰盆地基底发现晚寒武世-早奥陶世地层-来自LA-MC-ICP-MS锆石U-Pb年龄证据. 地球学报, 38(6): 859-871.
杨富全, 王永强, 杨成栋, 李强, 郭旭吉. 2018. 新疆阿尔泰克兰盆地金属矿床地质特征及成矿作用. 矿床地质, 37(3): 441-462.
杨富全, 李永, 陈贵民, 张国锋, 李超, 田猛, 李鹏, 李强. 2019. 新疆阿尔泰巴斯铁列克钨多金属矿成矿时代-来自锆石U-Pb年龄和辉钼矿Re-Os年龄约束. 大地构造与成矿学, 43(6): 1144-1154.
杨坤光, 杨巍然. 1997. 碰撞后的造山过程及造山带巨量花岗岩的成因. 地质科技情报, 16(04): 16-22.
原垭斌, 袁顺达, 陈长江, 霍然. 2014. 黄沙坪矿区花岗岩类的锆石U-Pb年龄、Hf同位素组成及其地质意义. 岩石学报, 30(1): 64-78.
张超, 刘锋, 杨富全. 2013. 新疆富蕴县库额尔齐斯铁矿成因机制探讨研究. 岩矿测试, 32(1): 157-165. DOI:10.3969/j.issn.0254-5357.2013.01.026
张国锋, 李永, 杨富全, 李强, 杨成栋, 柴凤梅. 2019. 新疆阿尔泰巴斯铁列克钨多金属矿矿物特征及其地质意义. 矿床地质, 38(1): 101-119.
章荣清. 2014. 湘南含钨和含锡花岗岩成因与成矿作用: 以王仙岭和新田岭为例. 博士学位论文. 南京: 南京大学
张永明, 裴先治, 李佐臣, 李瑞保, 刘成军, 裴磊, 陈有炘, 王盟. 2019. 青海南山地区加里东期强过铝质花岗岩锆石U-Pb年龄、地球化学特征及其地质意义. 地质通报, 38(5): 742-756.
张志欣, 柴凤梅, 王丽芬, 刘锋, 耿新霞. 2014. 新疆阿尔泰南缘克因布拉克铜锌矿区黑云母二长花岗岩地球化学特征及意义. 现代地质, 28(3): 461-471. DOI:10.3969/j.issn.1000-8527.2014.03.002
张作衡, 柴凤梅, 杜安道, 张招崇, 闫升好, 杨建民, 屈文俊, 王志良. 2005. 新疆喀拉通克铜镍硫化物矿床Re-Os同位素测年及成矿物质来源示踪. 岩石矿物学杂志, 24(4): 285-293. DOI:10.3969/j.issn.1000-6524.2005.04.005
赵玉梅, 彭戈, 赵得龙, 尤悦程. 2016. 中国阿尔泰克孜勒花岗岩体LA-ICP-MS锆石U-Pb定年、岩石地球化学特征及其地质意义. 地质与勘探, 52(2): 271-282.
周刚, 张招崇, 谷高中, 杨文平, 何斌, 张小林, 罗世宾, 王祥, 贺永康. 2006. 新疆东准噶尔北部青格里河下游花岗岩类的时代及地质意义. 现代地质, 20(1): 141-150. DOI:10.3969/j.issn.1000-8527.2006.01.017
祝新友, 王艳丽, 程细音, 田野. 2016. 湖南瑶岗仙花岗岩体中包体的地质地球化学特征与岩浆演化. 岩石矿物学杂志, 35(1): 16-32. DOI:10.3969/j.issn.1000-6524.2016.01.002