DOi:10.16111/j.0258_7106.2017.01.012
 海南岛高通岭钼矿床赋矿岩体LA-ICP-MS锆石U_Pb年龄、Hf同位素和Ce4+ /Ce3+特征 
 朱昱桦1,2,于亮亮3,4**,于得水1,2,许德如1,付王伟5

(1 中国科学院广州地球化学研究所, 中国科学院矿物学与成矿学重点实验室, 广东 广州 510640; 2 中国科学院大学, 北京100049; 3 中山大学地球科学与地质工程学 院, 广东 广州510275; 4 广东省地质过程与矿产资源探查重点实 验室, 广东 广州 510275; 5 哈密地区安全生产监督管理局, 新疆 哈密839000)

第一作者简介朱昱桦, 男, 1989年生, 硕士研究生, 矿物学、岩石学、矿床学专业。 Email: zhuyuhua23@163.com
**通讯作者于亮亮, 男, 1983年生, 博士, 主要从事构造与成矿学研究。 Email: sh ushengbeibei@163.com

收稿日期2015_11_26

本文受到国家重点研发计划“深地资源勘查开采”(编号: 2016YFC0600401)重点专项和 国土资源部全国矿产资源潜力评价项目(编号: 1212010881625)资助

摘要:高通岭钼矿床位于海南岛中部的早白垩世—晚白垩世高通岭岩体中, 是一个小型热 液脉型钼矿床。作者对赋矿围岩黑云母钾长花岗岩进行了岩相学观察、LA-ICP-MS锆 石U_Pb 定年、锆石Lu_Hf同位素测试和锆石微量元素分析。结果显示:赋矿围岩矿物组合为钾长石 (45%~60 %)、石英(20%~30%)、斜长石(10%~20%)及黑云母(10%~15%);矿石矿物以辉钼矿 为主 ,矿石类型主要为石英脉型,次为蚀变黑云母钾长花岗岩型和碎裂岩型;高通岭岩体中的锆 石,具有清晰的振荡环带,显示其岩浆成因;锆石U_Pb年代学研究获得206U/ 238Pb加权平均年龄为(102.5±1.8) Ma,锆石的176Hf/17 7Hf比值为0.282 349~0.282 663, εHf(t)为-12.8~-1.7 ,二阶段地壳模式年龄(TCDM)为1976~1271 Ma,峰值为1450~1400 Ma, 锆石稀土元素配分模式显示高通岭岩 体的陆壳源区属性。地表取样的黑云母钾长花岗岩(GTL_03和GTL_04)具有较低Ce4+ /Ce3+比值(160左右),反映了较低氧逸度,钻孔取样的黑云母钾 长花岗岩具有较高Ce4+/Ce3+比值(平均值377)。高通岭岩体总体上氧逸度较 低,暗示其大规模成矿潜力较低。
关键词: 地球化学;锆石U_Pb定年;Hf同位素;Ce4+/Ce3+比值; 高通岭钼矿床
文章编号: 0258_7106 (2017) 01_0185_15  中图分类号: P618.65 文献标志码:A
 LA-ICP-MS zircon U_Pb dating, Hf isotopic composition and Ce4+/Ce3+ 
    characteristics of Gaotongling molybdenum ore deposit in Hainan Provinc e 
    and their implications for metallogeny 
ZHU YuHua1,2, YU LiangLiang3,4, YU DeShui1,2, XU DeRu1 and FU WangWei5 

(1 CAS Key Laboratory of Mineralogy and Metallogeny, Guangzhou Institute of Geoc hemistry, Chinese Academy of Sciences, Guangzhou 510640, Guangdong, China; 2 Uni versity of Chinese Academy of Sciences, Beijing 100049, China; 3 School of Earth Science and Geological Engineering, Sun Yat_sen University, Guangzhou 510275, G uangdong, China; 4 Guangdong Key Laboratory on Geological Process and Mineral Re sources Exploration, Guangzhou 510275, Guangdong, China; 5 Hami Administration o f Work Safety, Hami 83900, Xinjiang, China)

2015_11_26

Abstract:The Gaotongling molybdenum ore deposit classified as the quartz vein type or e deposit is located in the central area of the Hainan Island. It is hosted in t he Cretaceous granite named Gaotongling terrace. In order to investigate minera l composition, metallogenic age and metallogenic source of the Gaotongling molyb denum deposit, the authors conducted petrographic observations, LA-ICP-MS zirc on U_Pb dating, zircon in situ Lu_Hf isotope analysis and trace elements analysis of the ore_hosting rocks. The ore_hosting rocks are medium_coarse grained biotit e moyite, and the ore types can be divided into quartz vein type molybdenite ore , altered moyite type molybdenite ore and cataclastic rock type molybdenite ore. LA-ICP-MS zircon dating of biotite moyite rock yielded 206U/ 238Pb ages varying from 105.10 Ma to 99.89 Ma with weighted average age of (102.5±1.8) Ma. LA -ICP- MS zircon Hf isotopic analyses show εHf(t) values of biotite moyite r ocks from -12.8 to -1.7. The two_stage crust model ages (TCDM) of Gaotonglin g terrace vary from 1976 Ma to 1271 Ma, with the peaks from 1450 Ma to 1400 Ma. Combined with REE characteristics of zircons, the authors hold that the Gaotongl ing rocks belong to continental crust. Early Gaotongling granite rocks have low Ce4+/Ce3+ ratios (around 160), reflecting a low oxygen fugacity, whe reas Gaotongling granite samples collected from the drill hole have high Ce4 +/Ce3+ ratios (377 on average), showing a higher oxygen fugacity. Gaoton gling granites possess confined metallogenic capacity. 
Key words: geochemistry, zircon U_Pb dating, Hf isotope, Ce4+/Ce 3+ ratio, Gaotongling molybdenum deposit
         近年来,随着找矿勘查工作的开展,已在海南岛发现9个钼矿床(陈沐龙等,2002),分别 是罗 葵洞大型钼矿床(刘君等,2010; 王国君等,2010; 于立红等,2010; 于立红, 2014)、石 门山中型钼铅锌多金属矿床、新村中型钼矿床、龙门岭中型钼矿床、高通岭小型钼矿床(廖 香俊等,2008)、红门岭小型钼钨矿床(付王伟等,2013)、报告村小型钼矿床、梅岭小型 铜钼矿床(龙文国等,2003)和红岭钼矿点(图1)。这一系列矿床时代相近,成矿动力学 背景相同,矿床类型相似,且大多数分布在九所-陵水深大断裂两侧,因此,前人将这些矿 床统 一归为石门山式钼矿床(廖香俊等,2008)或海南岛早白垩世晚期至晚白垩世早期(120~7 0 Ma)钼成矿作用的产物(付王伟, 2013; 李孙雄等,2014)。地理位置上,高通岭钼矿床 位 于昌江-琼海深大断裂带,梅岭铜钼矿床分布在海南岛东北部,其余矿床均分布在九所-陵 水 断裂两侧的岩体中。矿床类型上,红门岭钼钨矿床、龙门岭钼矿床和高通岭钼矿床均属热液 脉 型钼矿床,其余钼矿床均为斑岩型钼矿床(陈沐龙等,2002; 付王伟等,2011;李孙雄等, 2014)。由此可见,高通岭钼矿床在地理位置上和矿床类型上与其他钼矿床存在着差异,这 些钼矿床是否能统一划成一类仍需要更多证据的支持。
        前人对高通岭钼矿床的成矿时代进行过初步的研究。廖香俊等(2008)对高通岭钼矿床进行 了 辉钼矿Re_Os同位素测年,获得了(98.4±2.5) Ma的等时线年龄。陈沐龙等(2002)对 高通岭岩体进行锆石U_Pb定年,获得了(95.6±1.6) Ma的年龄,并认为这一年龄代表 该矿床形成的年龄。贾小辉等(2010)和Wang等(2012)对高通岭岩体进行LA-ICP -MS锆石U_Pb定年,获得(108.1±1.7) Ma的加权平均年龄。付王伟等(2014)对高 通岭岩体进行定年,获得了加权平均年龄(100.8±2.4) Ma和(100.2±1.7) Ma。那 么,这些年龄的意义分别是什么?到底哪一年龄可以代表高通岭钼矿的形成年龄?这些问题 尚需进一步探讨。
        高通岭钼矿床成因的研究工作尚处于起步阶段,对其矿床成因和成矿物质来源相关研究不多 。廖香俊等(2008)通过矿石石英中包裹体的C、H、O同位素研究,推测该矿床为燕山晚期 岩 浆活动形成的热液脉型矿床,并推测高通岭钼矿成矿流体来自于高通岭岩体。付王伟等( 20 13;2014)根据高通岭岩体年龄与辉钼矿年龄非常相近,推测高通岭钼矿的成矿流体直接来 自于高通岭岩体。高通岭的岩浆性质、氧逸度和成矿过程等尚不明确,仍有待进一步研究。 
        本文对高通岭钼矿床赋矿围岩黑云母钾长花岗岩进行了LA-ICP-MS锆石U_Pb定年、锆石 微 量元素含量和锆石Hf同位素组成分析,以探讨高通岭岩体及钼矿床形成年龄、高通岭岩体的 源区、岩体的氧逸度和含矿能力,从而为研究高通岭钼矿床的成因提供制约。
图 1海南岛区域地质简图(据Xu et al., 2013)
     Fig. 1Simplified map showing the main stratigraphic and magmatic units, and or e deposits in the Hainan Island, 
    south China (after Xu et al., 2013)        
1成矿地质背景及矿床地质特征
        海南岛是一个大陆性岛屿,通过琼州海峡与华南大陆相连(Xu et al., 2013; 2014)。海 南岛地处欧亚板块、印度-澳大利亚板块和太平洋板块的交汇地带,具有复杂的大地构造特 征和构造演化历史,经历了晋宁、加里东、海西-印支期、燕山期和喜马拉 雅构造运动(付王伟等,2014;Xu et al., 2013; 2014)。海南岛出露地层主要为古生界 ,其次为元古宇和新生界(图1)。最老地层为中元古代抱板群,但是出露较为零星且 被后期侵入 岩体切割(陈哲培等,1997)。海南岛大量发育印支期和燕山期侵入岩,约占整个海南岛面 积的60%(陈哲培等,1997)。海南岛主要发育EW向和NE向构造带(图1)。自北向南, 海南岛分布4条EW向构造带,分别是王五-文教断裂、昌江-琼海断裂、尖峰-吊罗断裂和 九所-陵水断裂(图1)。NE向构造带主要包括白沙断裂和著名的含金矿带——戈枕 韧性 剪切带。以九所-陵水深大断裂为界,海南岛的北部是华南地层大区的五指山地体,南部是 南海大区的三亚地体(陈哲培等,1997)。
        高通岭钼矿床位于昌江-琼海断裂中段燕山晚期的高通岭岩体内(付王伟等,2011) 。高通岭黑云母钾长花岗岩呈岩株状分布在屯昌复式岩体之中,出露面积约为3. 0 km2。矿区内发育大量燕山晚期岩脉,主要为花岗斑岩脉、煌斑岩脉(图2a),其次有伟 晶岩脉、 辉石闪长岩脉、花岗细晶岩脉和闪长岩脉等。岩脉规模及产状变化较大,常呈平行排列,具 分 叉复合现象。高通岭岩体内发育近SN向、NW向断裂和破碎带,矿体严格受控于NW向断裂破碎 带。
        高通岭钼矿床矿体呈脉状分布在高通岭岩体中(图2b)。矿化带以高通岭为中心,地表宽度 约200 m、长度约600 m,走向约115~123°,倾向NE,倾角约49~85°(图2a、b)。矿化 带 中分布14个钼矿体和矿化体(图2a),其中9个矿体可达到最低可采厚度和最低工业品位的 要求。矿体NW向平行分布,间距一般10~20 m(图2b)。矿体中心一般为含辉钼矿的石英- 多金属硫化物脉,矿体两侧发育浸染状矿化现象。矿体地表长度一般在200~400 m,最长约 529 m(V9矿体),厚一般1~2 m,最厚3.61 m(V2矿体),倾向延深一般100~200 m,最深约390 m(V2矿体)。矿体的钼品位为0.06%~0.28%,最高达2.70%,全 矿床平均钼品位为0.1%。截止2006年,该矿床已查明钼金属储量1697.5 t。 
     图 2高通岭钼矿床矿区地质简图(a,据付王伟等,2014)及矿床地质剖面图(b,据付王 伟,2013)
     Fig. 2Geological map of the Gaotongling molybdenum ore deposit(a, after Fu et al., 2014)and geological 
    section of the Gaotongling molybdenum ore deposit( b, after Fu, 2013)  
图 3高通岭钼矿区赋矿围岩和矿石照片
     a. 黑云母钾长花岗岩手标本; b. 正交偏光下显微照片; c. 石英脉型钼矿石; d. 自形 叶片状、板片状辉钼矿;
    e. 蚀变钾长花岗岩型钼矿石; f. 碎裂岩型钼矿石
     Kf—钾长石; Q—石英; Pl—斜长石; Bi—黑云母; Ser—绢云母; Mo—辉钼矿
     Fig. 3The photographs for the rocks and the ore of the Gaotongling molybdenum ore deposit
     a. Biotite moyite sample; b. Microscopic photograph(crossed nicols); c. Quartz vein type molybdenum ore; d. Euhedral foliated and plate_like
     molybdenite; e. Altered biotite moyite type molybdenum ore; f. Cataclastic rock type molybdenum ore
     Kf—K_feldspar; Q—Quartz; Pl—Plagioclase; Bi—Biotite; Ser—Sericite; Mo—Moly bdenum ore        
2赋矿围岩和矿石矿物组成
        高通岭钼矿床赋矿围岩为黑云母钾长花岗岩 (图3a),呈肉红色-灰白色,花岗结构,局部可见钾长 石斑晶,块状构造。主要矿物成分 为钾长石(45%~60%,微斜长石为主)、石英(20%~30%)、斜长石 (10%~20%)及黑云母(10%~15%),副矿物有锆石、磷灰石、磁铁矿 等, 次生矿物有高岭石、绢云母等(图3b)。按照矿物颗粒大小,可分成细粒黑云母钾长花岗岩 与粗粒黑云母钾长花岗岩,两者呈渐变关系。
        矿石类型主要为石英脉型(图3c),次为蚀变 黑云母钾长花岗岩型和碎裂岩型。石英脉型矿石通常是在石英脉两侧发育厚层状辉钼矿,辉 钼矿通常以粗粒板片状和叶片状辉钼矿集合体形式存在(图3d)。矿石矿物以辉钼矿为主, 含少量黄铁矿、黄铜矿、闪锌矿、黑钨矿等,脉石矿物为石英、斜长石、钾长石、绢云母、 绿泥石等。蚀变黑云母钾长花岗岩型钼矿石,通常产于石英脉的外侧,强烈钾化,辉钼矿则 呈细粒浸染状分布其中(图3e)。碎裂岩型钼矿石主要发育在断裂破碎带中,石英脉 中含石英脉碎块、断层泥以及细粒状和团块状辉钼矿(图3f)。
3样品及分析方法
        本文对1件样品(GTL_03)进行LA-ICP-MS锆石U_Pb定年、锆石微量元素含量和锆石Hf同位 素 组成分析,2件样品(GTL_04和GTL_09)进行锆石微量元素含量和锆石Hf同位素组成分析。 这3件样品均采自高通岭钼矿区的高通岭岩体。GTL_03样品取自V2矿脉坑道口附近的细粒 黑云母钾长花岗岩(图2b);GTL_04样品取自V2矿脉坑道口南 16 m处的近矿细粒黑云母钾长花岗岩;GTL_09样品取自钻孔ZK22702的127回次265 m 处(图2a、b),为粗粒黑云母钾长花岗岩。
        本次所采样品新鲜,基本无矿化,部分样品 发生轻微蚀变。先对样品清洗后,挑选一部分新鲜的岩石,粉碎至120目左右 ,经电磁选和重选富集锆石。在双目镜下挑选出晶形完好的锆石颗粒,将锆石制靶并拍摄阴 极 发光图像(CL)、透射光图像和反射光图像。锆石的分选、阴极发光图像、透射光图像和反 射光图像拍摄均在北京锆年领航科技公司完成。在阴极发光图像上选取适合的位置,首先进 行LA-ICP-MS锆石U_Pb定年和锆石微量元素测试,之后在对应的测点上进行Lu_Hf同位素测 试 ,上述测试均在中国科学院广州地球化学研究所同位素地球化学国家重点实验室完成。LA- ICP-MS锆石U_Pb定年、锆石微量元素测定和Lu_Hf同位素测定均采用Resolution M_50型激 光 剥蚀系统和Agilent 7500a型的ICP_MS联机。锆石U_Pb定年、锆石微量元素含量测定时,采 用标准锆石TEMORA (417 Ma)作为测年外标,测试时剥蚀激光斑束直径为31 μm,频率为8 Hz,分析过程见涂湘林等(2011)。锆石Lu_Hf同位素组成测试时,剥蚀激光斑束直径为44 μm,频率为8 Hz。锆石U_Pb定年和锆石微量元素数据处理使用ICPMSDataCal 8.6软件(L iu et al., 2008),谐和图的绘制采用Isoplot v3.0软件完成。
4数据分析
4.1锆石U_Pb定年
        LA-ICP-MS锆石U_Pb同位素数据见表1和图4。图4所显示的是GTL_03样品的锆石颗粒阴极发 光 图像(CL),测点位置及锆石U_Pb年龄均进行了标注。大部分锆石为长柱状,晶形完好,颗 粒较大,锆石颗粒的长度大多在70~280 μm,长/宽比在1.5~4之间。其中较大的锆石颗 粒反 映较饱和的锆含量。绝大部分锆石颗粒具有清晰明显的震荡环带,表明这些锆石颗粒属于岩 浆成因锆石(Hoskin et al., 2003)。锆石Th/U比值在0.55~0.86之间, 也反映出典型的岩浆锆石的Th/U比值特征(吴元保等,2004;高晓英等,2011)。从震荡环带的交切关系可发 现,大部分锆石颗粒都具有核,这些处于核部锆石同样具有明显的震荡环带。在锆相对饱和 的岩浆中,这种结构是由于锆石沿着尚未溶解的锆石周 围继续结晶而形成(Hoskin et al., 2003)。 
 表 1高通岭钼矿花岗岩围岩锆石LA-ICP-MS U_Pb年龄定年结果
     Table 1LA-ICP-MS zircon U_Pb data of granite in the Gaotongling molybdenum o re deposit   
图 4高通岭钼矿赋矿花岗岩样品GTL_03的锆石CL图像
     实线圆为LA-ICP-MS U_Pb年龄测点,虚线圆为Hf同位素测试点
     Fig. 4Cathodoluminescence (CL) images of zircon from granite of the Gaotonglin g molybdenum ore deposit
     White solid circles indicate locations of LA-ICP-MS U_Pb dating and dotted cir cles indicate LA_MC_ICP_MS Hf analytical 
    positions, the U_Pb age is also shown    
       图 5高通岭钼矿区赋矿花岗岩的锆石U_Pb年龄谐和图和加权平均年龄图(GTL_04和GTL_09 引自付王伟等,2014)
     Fig. 5Zircon concordia diagrams for the granite of the Gaotongling molybdenum ore deposit and its average age
     (GTL_04 and GTL_09 after Fu et al., 2014)       
        在锆石U_Pb年龄分析时,对于谐和度小于90%的点予以舍弃。GTL_03样品的206Pb /238U年龄分布在99.89~105.97 Ma之间,获得了(102.5± 1.8) Ma的 206Pb/238U加权平均年龄(图5)。这一年龄与付王伟等(2014)采用LA -ICP-MS方法获得的高通岭岩体锆石206Pb/238U加权平均年龄(100 .8±2.4) Ma和(100.2±1.7) Ma在误差范围内一致。
4.2锆石Lu_Hf同位素组成
        对GTL_03、GTL_04和GTL_09三件黑云母钾长花岗岩样品进行了锆石原位Hf同位素组成测试。 Hf同位素测试结果如表2所示。GTL_03黑云母钾长花岗岩锆石的176Hf/1 77Hf比值范围分布在0.282 44~0.282 61, 176Lu/177 H f比值范围分布在0.000 42~0.001 38。 GTL_04黑云母钾长花岗岩锆 石的 176Hf/177Hf比值范围分布在0.282 38~0.282 66,比值范围分布在0.000 52~0.002 10。GTL_09黑云母钾长 花岗岩锆石的176Hf/177Hf比值范围分布在 0.282 35~0.282 6 6,176Lu/177Hf比值范围分布在0.000 57~0.001 05 。
表 2高通岭钼矿赋矿黑云母钾长花岗岩锆石Hf同位素组成
     Table 2Zircon Hf isotope data for the ore_bearing granite of the Gaotongling m olybdenum ore deposit
    
图 6高通岭黑云母钾长花岗岩锆石的地壳Hf同位素模式年龄(TCDM)分布图( a)和εHf(t)_年龄图(b)
     Fig. 6Distribution of zircon Hf isotope crust model ages of Gaotongling granit e (TCDM) (a)and the plots of the εHf(t
    versus U_ Pb age diagram for the zircon of ore_bearing granite (b)176 Lu/177Hf    
        高通岭赋矿围岩黑云母钾长花岗岩的εHf(t)分布在-12.8~-1.7,T CDM主要分布在1976~1271 Ma,主要峰值落在1450~1400 Ma范围内(图6a)。图6 b可以发现高通岭钼矿εHf(t)大部分落在1.5~ 1.3 Ga演化线之上,但很 少一部分锆石的εHf(t)落在1.5~2.0 Ga上地壳演化线之间,这说明高通 岭岩体母岩主要由1.5~1.3 Ga从亏损地幔分离出来的物质加上小部分更老的上地壳物质 混合形成。这一小部分稍老的上地壳物质形成时间不晚于2.0 Ga。
4.3锆石微量元素特征与Ce4+/Ce3+比值
        锆石中,Hf、Y、P、Th、U、REE等元素常会以类质同象替换锆石中的Zr(雷玮琰等,2013) 。磷钇矿(化学式Y(PO4))具有与锆石相类似的结构,在锆石中磷钇矿、磷灰石等是常见 的 包体。从阴极发光图像可以发现,高通岭锆石中也含有包体(图4),虽然在实验过程中尽 量避开包体,但一些微小的包体是难以避免的。那些具有异常高的P(>500×10 -6)、Y(>2000×10-6)和Ca含量的测点,很可能含有包体,因此在分析锆石 微量元素时需进行舍弃。表3为黑云母钾长岗岩样品GTL_03、GTL_04和GTL_09中锆石微量元 素含量。
        锆石的球粒陨石标准化稀土元素配分模式表明(图7),GTL_03、GTL_04和GTL_09样品中的锆 石明显 的富集重稀土元素,具有Ce正异常和Eu负异常等特征,并落在大陆地壳锆石成因范围内 (Grimes et al., 2007)。GTL_03样品中锆石的Eu/Eu*比值范围在0.31~0.64(表3),平 均值为0.45;GTL_04样品中锆石的Eu/Eu*比值为0.19~0.55(表3),平均值为0.32;G TL_09样品中锆石的Eu/Eu*比值为0.23~0.51(表3),平均值为0.41。 
     岩浆中Ce3+离子可以被氧化成Ce4+,Ce4+离子具有与Zr4+离子相 同的电荷数以及相近的离子 半径,所以会更容易取代锆石中的锆离子,因此锆石会呈现出强烈的Ce正异常(Liang et al., 2006)。Ce4+/Ce3+主要受氧逸度的控制,受温度和压力影 响较 小(Ballard et al., 2002; Liang et al., 2006),因此,锆石中Ce4+/Ce3+比值可以反映岩浆的氧逸 度。目前Ce 4+/Ce3+尚不能通过光谱分析直接测得,但是可以通过相关理论计算求得(Ball ard et al., 2002),公式为:
     (Ce(Ⅳ)/Ce(Ⅲ))zircon=(Cemelt-Cezircon/Dzircon/melt Ce(Ⅲ))/
    (Cezircon/Dzircon/meltCe(Ⅳ) -Cemelt)(1)
        根据Blundy等(1994)提出的晶格扭曲模型,推导出分配系数与半径的关系,即: 
lnDi=lnD0-4πENA(ri/3+r0/6)(ri-r 0)2/RT(2) 
        同一岩浆中温度和压力一定,Di为某元素在锆石-熔体相分配系数,D0为假定 的无应变时离子的理想分配系数,R为气体常数,E为杨氏模量,NA为阿伏加 德罗常数,ri为某元素的离子半径,r0为最优离子半径,即Zr离子八 配位的离 子半径0.84 ?。由式(2)可知,lnDi与(ri/3+r0/6)( ri-r0)2线性相关。
        取用三价稀土元素和四价Hf、Th和U的lnDi与相应(ri/3+r0/6)( ri-r0)2作图(图8,以GTL_03_11为例)可以计算出Ce3+和Ce4+ 在锆石-熔体相的分配系数(表4),更多细节可以参考(Ballard et al., 2002; 辛洪波 等,2008;张俊杰等,2012)。将Ce3+和Ce4+的锆石-熔体相分配系数,岩石 样品的Ce含量(Cemelt)和锆石Ce含量(Cezircon)带入式(1),即可 计算出Ce4+/Ce 3+比值。通过上述方法计算出GTL_03样品的Ce4+/Ce3+比值范围在17~322 ,平均值为166;GTL_04的Ce4+/Ce3+比值范围在20~412,平均值为155;G TL_09 的Ce4+/Ce3+比值比前两者明显增加,范围在87~1053之间,平均值为377。
5讨论
5.1高通岭岩体年龄及源区
        目前对于高通岭矿床的年龄有以下4种结果:
        (1) 陈沐龙等(2002)给出高通岭钼矿的年龄为(95.6±1.6) Ma;
        (2) 廖香俊等(2008)的辉钼矿Re_Os定年结果为(98.4±2.5) Ma;
        (3) 贾小辉等(2010)和Wang等 (2012)获得高通岭岩体锆石U_Pb年龄为(108.1±1. 7) Ma(20个测点);
        (4) 付王伟等(2014)对高通岭钼矿围岩进行定年,获得加权平均年龄(100.8±2.4) Ma和(100.2±1.7) Ma。
表 3黑云母钾长花岗岩锆石微量元素含量
     Table 3Trace elements data for the zircon of the ore_bearing granitoid    
 图 7高通岭黑云母钾长花岗岩锆石的REE球粒陨石标
    准化图(Grimes et al., 20 07)(球粒陨石据Sun et al., 
    1989)
     Fig. 7Chondrite_normalized REE patterns of zircon from Gaotongling granite (af ter Grimes et al., 2007) (chondrite 
    values after Sun et al., 1989)   
         这些年龄都曾被解释为高通岭钼矿床的年龄。其中(95.6±1.6) Ma这一年龄数据在发表 时,具体样品个数、测试点数、数据和采样地点均没有列出,这一年 龄明显偏小,可能是后期侵入脉体的年龄(陈沐龙等,2002)。廖香俊等(2008)得到的辉钼 矿Re_ Os年龄为(98.4±2.5) Ma,样品采自含矿石英脉,代表高通岭钼矿的成矿年龄。贾小辉 等(2010)和Wang等(2012)在对高通岭岩体进行锆石U_Pb定年的同时,还对屯昌 岩体进行了锆石U_Pb定年,获得了(107.2±1.3) Ma加权平均年龄。高通岭岩体作为屯 昌岩体的侵入岩体,其年龄应小于屯昌岩体的年龄。相比之下,高通岭岩体年龄(108.1± 1.7) Ma偏 大,且其测点偏少。本文获得的高通岭赋矿岩体黑云母钾长花岗岩中锆石U_Pb定年结果(( 102.5±1.8) Ma)与付王伟等(2014)定年结果((100.8±2.4) Ma和(100.2±1 .7) Ma)在误差范围内一致,可以代表高通岭钼矿赋矿围岩的年龄。
        目前,在海南岛的白 垩 纪钼矿床中,只有高通岭钼矿床、石门山 钼矿床((100±1) Ma;陈沐龙等,2014)和罗葵洞钼 
    表 4微量元素八倍配位的离子半径值和含矿斑岩微量元素在锆石_熔体相的分配系数
     Table 4Ionic radii of trace elements as cations in 8_fold coordination and the distribution coefficients for trace elements
     between individual zircon and wh ole rock for the ore_bearing granitoid
    图 8分配系数与离子半径函数的线性拟合图
     Fig. 8The liner fitting between lnDi and function of ri
        矿床((99.7±0.4) Ma,李孙雄等,2014)具有可靠的年 龄数据。在华南,90~100 Ma这一时期,古太平洋板块(或称伊泽乃奇板块)向欧亚板块俯 冲(毛景文等,2008)而产生大规模的成矿事件(孙卫东等,2008)。海南岛钼矿床的成矿 时间与华南这一成矿事件的时间吻合,因此,推测海南岛白垩纪钼矿床成矿带很可能是华南 白垩纪成矿带的南缘。
        锆石Hf同位素模式年龄主要峰值为1450~1400 Ma。在海南岛,峰值为1477 Ma的碎 屑锆石(Zhou et al., 2015)和(1439±9) Ma的花岗岩体(Li et al., 2008)都已经 被 确认,这些碎屑锆石和岩体的模式年龄大多分布在1.6~1.8 Ga范围内(梁新权,1995) 。虽 然尚不能确定高通岭岩体是不是这些岩石重熔的产物,但可以确定1400~1450 Ma在海南岛 有 大规模岩浆事件。高通岭花岗岩锆石的εHft)值均小于0,表明这些 岩石的源区以地壳为主。锆 石的稀土元素配分模式与壳源锆石的稀土元素分配模式相似(图7),进一步证实了高通岭 岩体的源区为地壳岩石。
5.2高通岭岩体与钼矿床成矿过程
        目前对于高通岭钼矿床属于石英热液脉型矿床已经取得共识(陈沐龙等 ,2002; 廖香俊等,2008;付王伟等,2011;2014)。前人研究表明,高通岭钼矿床的成 矿 流体源于岩浆流体(廖香俊等,2005;2008;付王伟等,2014;李孙雄等,2014)。岩浆流 体形成的热液脉型硫化物矿床与斑岩型矿床有着重要联系,二者存在类似的成矿机制(江思 宏等,2004;凌洪飞,2011)。氧逸度可以影响进而控制长英质岩浆的矿化作用( Mungall, 2002; Sun et al., 2004)。在岩浆中,氧逸度控制着硫的价态,在低的氧逸度 时,硫会以S2-的形式存在;当氧逸度较高时,硫则以SO或SO2的形式存在(Liang et al., 2006)。SO和SO2与Cu、Mo等离子更具亲和性,且有较高的溶解度,随着岩浆分 异,流体-熔体中会逐渐富集Cu、Mo等元素,故在高氧逸度时,流体中会富集更多Cu、 M o金属物质,从而有利于成矿。因此,高的氧逸度常作为斑岩型矿床成矿的标志之一(Balla rd et al., 2002; Liang et al., 2006; 辛洪波等,2008)。Ce4+/Ce3+比值 大小能反映出岩浆氧逸度高低(Trail et al., 2012; 赵振华,2010),国内外学者根据Ce 4+/Ce3+比值和Eu/Eu*界定斑岩型矿床的成矿岩体与不成矿岩体,Ce4+ /Ce3+临界值一般分布在260~300之间,Eu/Eu*为0.4,大于这2个临界值才具 有 成矿潜力(Liang et al., 2006; 辛洪波等,2008;赵振华,2010)。本文的Ce4+/C e3+比值计算结果表明,黑云母钾长花岗岩样品GTL_03、GTL_04的 Ce4+/Ce3+比值较低,说明氧逸度偏低;而黑云母钾长花岗岩样品GTL_09的Ce 4+/Ce3+比值高,具有较高的氧逸度。表4定年数据 显示GTL_03、GTL_04样品年龄数据稍大于GTL_09样品,且GTL_09样品采样位置位于265 m深 处钻孔内, 因此认为GTL_09粗粒黑云母钾长花岗岩样品代表晚期岩浆结晶的产物。
        高通岭岩体早期结晶岩石具有较低氧逸度,在岩浆演化晚期,氧逸度增加,此时岩浆流体得 以携带更多成矿物质。高通岭岩体不同部位氧逸度存在差异的原因,可能是岩浆演化的结果 ,也可能是局部有高氧逸度物质混染(凌洪飞,2011)。总体上,高通岭岩体氧逸度较低, 不具大规模成矿的潜力。
6结论
        通过对高通岭钼矿床赋矿围岩和矿石进行岩相学观察、赋矿围岩的LA-ICP-MS锆石U_Pb定 年、锆石Lu_Hf同位素组成测试和锆石微量元素分析得出以下结论:
        (1) 高通岭钼矿赋矿围岩是黑云母钾长花岗岩。矿石类型可划分为石英脉型、蚀变黑云母 钾长花岗岩型和碎裂岩型。
        (2) 高通岭岩体的锆石206U/238Pb加权平均年龄为(102.5±1.8 ) Ma。
        (3) 高通岭岩体的锆石εHf(t)为-12.8~-1.7,二阶段地壳模式年 龄 (TCDM)分布在1976~1271 Ma,峰值为1450~1400 Ma。结合球粒陨石标准化锆 石稀土元素配分模式图,显示高通岭岩体的陆壳源区属性。
        (4) GTL_03、GTL_04黑云母钾长花岗岩样品的锆石Ce4+/Ce3+比值平均分别 为166、155,GTL_09黑云母钾长花岗岩样品的锆石Ce4+/Ce3+比值高,平均 值为377,显示高通岭岩体 总体上氧逸度偏低,大规模钼成矿作用受限。    
        志谢本文的野外工作得到海南省地矿局的支持,在此表示感谢。同时感谢 匿名审 稿专家和编辑老师对本文的审阅及宝贵意见,这使得本文质量有了提高。           
参考文献
References    
     Ballard J R, Palin M J and Campbell I H. 2002. Relative oxidation states of magm as inferred from Ce(Ⅳ)/Ce(Ⅲ) in zircon: Application to porphyry copper deposi ts of northern Chile[J]. Contributions to Mineralogy and Petrology, 144: 347_3 64.
    Blundy J and Wood B. 1994. Prediction of crystal_melt partition coefficient s fro m elastic moduli[J]. Nature, 372:452_454.
     Chen M L, Fu C R, Zhuang Y G and Yun P. 2002. First discovery of quartz hydrothe rmal vein type molybdenum deposit in QiongZhong region[J]. Geology and Mineral Resources of South China, (1): 9 (in Chinese).
     Chen M L, Ma C Q, Lv Z Y, Yun P and Liu Y Y. 2014. Zircon U_Pb chronology of Qia njia rock mass in Hainan Island and its geological significance[J]. Geological Science and Technology Information, 33(6):1_10+27 (in Chinese).
     Chen Z P. 1997. Stratigraphy (Lithostratic) of Hainan Province[M]. Wuhan: Chin a University of Geosciences Press. 1_100 (in Chinese). 
     Fu W W, Xu D R, Fu Ya R, Yang C S, Zhou Y C, Yang D S, Wu C J and Wang Z L. 2011 . Regional metallogenic regularity of molybdenum deposit in Hainan Island[J]. Journal of Minerals, (S1):772_773 (in Chinese with English abstract).
     Fu W W. 2013. An research on ore deposit_types and geodynamic background of the Mo ore deposits in Hainan Province, South China[D]. Supervisor: Xu D R. Guangz ho u: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences. 184p (in Ch inese with English abstract).
     Fu W W, Xu D R, Fu Y R, Wu C J, Yang C S, Zhou Y C and Wang Z L. 2013. Molybdeni te Re_Os isotopic dating of Hongmenling Mo_W deposit in Hainan Province and its geological implications[J]. Journal of East China Institute of Technology, 35( 2):135_142 (in Chinese with English abstract).
     Fu W W, Xu D R, Wu C J, Fu Y R, Zhou Y C, Zhou Y Q, Wang Z L and Lin G. 2014. LA -ICP-MS zircon U_Pb dating of syenogranites hosting Gaotongling Mo deposit in Ha inan Province: Implications for metallogenesis[J]. Mineral Deposits, 33(2): 41 9_427(in Chinese with English abstract).
     Gao X Y and Zheng Y F. 2011. On the Zr_in_rutile and Ti_in_zircon geothermometer s[J]. Acta Petrologica Sinica, 27(2): 417_432(in Chinese with English abstract ).
     Grimes C B, John B E, Kelemen P B, Mazdab F K, Wooden J L, Cheadle M J, Hangh?j K and Schwartz J J. 2007. Trace element chemistry of zircons from oceanic cr ust: A method for distinguishing detrital zircon provenance[J]. Geology, 35: 643. 
     Hoskin P W O and Schaltegger U. 2003. The composition of zircon and igneous and metamorphic petrogenesis[J]. Reviews in Mineralogy and Geochemistry, 53:27_62. 
     Jia X H, Wang Q, Tang G J, Jiang Z Q, Zhao Z H, Yang Y H, Wang X D and Zhao W Q. 2010. Zircon U_Pb geochronology, geochemistry and petrogenesis of the Late Earl y Cretaceous adakitic intrusive rocks in the Tunchang area, Hainan Province[J] . Geochimica, 39(6): 497_519 (in Chinese with English abstract).
     Jiang S H, Nie F J, Zhang Y and Hu P. 2004. The latest advances in the research of epithermal deposits[J]. Earth Science Frontiers (China University of Geosci ences, Beijing), 11(2):401_411(in Chinese with English abstract).
     Lei W Y, Shi G H and Liu Y X. 2013. Research progress on trace element character istics of zircon different origins[J]. Earth Science Frontiers, 20(4): 273_284 (in Chinese with English abstract).
     Li S X, Chen M L, Yang D S and Wang Y H. 2014. The molybdenite Re_Os age and ana lysis of geodynamic background in Hainan island[J]. Geology and Mineral Re sources of South China, 30(3): 272_279(in Chinese with English abstract).
     Li Z X, Li X H, Li W X and Ding S J. 2008. Was Cathaysia part of Proterozoic La urentia? _ new data from Hainan Island, South China[J]. Terra Nova, 20(2): 154_164. 
     Liang X Q. 1995. Sm_Nd Ages of the Precambrian granitic_greenstone series in Hai nan Island and their geological significance[J]. Acta Petrroologica Sinica, 11 (1):71_76 (in Chinese with English abstract).
     Liao X J, Wang P A, Ding S J, Huang X D, Dong F X, Liu X C and Lei W Z. 2005. Ma i n minerogenetic series and metallogenic characteristics on Hainan Island[J]. J ournal of Geomechanics, 11(2): 187_194 (in Chinese with English abstract).
     Liao X J, Wang P A, Qin H C, Lu X K, Dong F X, Liu X C and Shu B. 2008. Geology, geochemistry and ore_forming age of the Gaotongling molybdenum deposit, Tunch ang area, Hainan, China [J]. Geological Bulletin of China, 27(4):560_570 (in Chinese with English abstract).
     Liang H Y, Campbell I H, Allen C, Sun W D, Liu C Q, Yu H X, Xie Y W and Zhang Y Q. 2006. Zircon Ce4+/Ce3+ ratios and ages for Yulong ore_bearing por phyries in eastern Tibet[J]. Mineralium Deposita, 41:152_159.
     Ling H F. 2011. Theory of granite type uranium deposit hydrothermal origin _ fro m oxygen fugacity condition restriction[J]. Geological Review, 57(2):193_206 (in Chinese with English abstract).
     Liu J, Cao Y L, Yu L H, Wang G J and Wang Y J. 2010. Wallrock alteration and pro specting criteria of the luokuidong molybdenum deposit in Hainan Province[J]. Geology and Resources, 19(1): 63_66 (in Chinese with English abstract).
     Liu Y, Hu Z, Gao S, Günther D, Xu J, Gao C and Chen H. 2008. In situ analysis o f major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard[J]. Chemical Geology, 257:34_43.
     Long W G, Ling Q Y, Fan Q H, Xie Z S, Zhang Y G and Yang Z Q. 2003. Geological f eatures of Yantang Cu_Mo deposit in eastern Hainan Island[J]. Geology and Mine ral Resources of South China, 30(3): 30_36(in Chinese with English abstract).
     Mao J W, Xie G Q, Guo C L, Yuan H D, Cheng Y B and Chen Y C. 2008. Spatial_tempo ral distribution of Mesozoic ore deposits in South China and their metallogenic settings [J]. Geological Journal of China Universities, 14(4): 510_526 (in Chi nese with English abstract).
     Mungall J E. 2002. Roasting the mantle: Slab melting and the genesis of major Au and Au_rich Cu deposits[J]. Geology, 30(10): 915_918.
     Sun S S and McDonough W F. 1989. Chemical and isotopic systematics of oceanic ba salts: implications for mantle composition and processes[J]. Geological So ciety, London, Special Publications, 42(1): 313_345.
     Sun W D, Arculus R J, Kamenetsky V S and Binns R A. 2004. Release of gold_bearin g fluids in convergent margin magmas prompted by magnetite crystallization[J]. Nature, 431(7011): 975_978.
     Sun W D, Ling M X, Wang F Y, Ding X, Hu Y H, Zhou J B and Yang X Y. 2008. Pacifi c Plate subduction and Mesozoic geological event in eastern China[J]. Bulletin o f Mineralogy, Petrology and Geochemistry, 27(3): 218_225 (in Chinese with En glish abstract).
     Trail D, Bruce W E and Tailby N D. 2012. Ce and Eu anomalies in zircon as proxie s for the oxidation state of magmas[J]. Geochimica et Cosmochimica Acta, 97:70 _87.
     Tu X L, Zhang H, Deng W F, Ling M X, Liang H Y, Liu Y and Sun W D. 2011. Applica tion of resolution in_situ laser ablation ICP_MS in trace element analyses[J]. Geochimica, 40(1): 83_98 (in Chinese with English abstract).
     Wang G J, Liu J, Cao Y L and Zhang X N. 2010. Metallogenic characteristics of Lu okuidong porphyry molybdenum deposit in the Baoting County, Hainan Province[J] . Mineral Exploration, 1(5): 453_457 (in Chinese with English abstract).
     Wang Q, Li X H, Jia X H, Wyman D, Tang G J, Li Z X, Ma L, Yang Y H, Jiang Z Q an d Gou G N. 2012. Late Early Cretaceous adakitic granitoids and associated magnes ian and potassium_rich mafic enclaves and dikes in the Tunchang_Fengmu area, Hai nan Province (South China): Partial melting of lower crust and mantle, and magma hybridization[J]. Chemical Geology, 328: 222_243.
     Wu Y B and Zheng Y F. 2004. Genesis of zircon and its constraint on interpretati on of U_Pb age[J]. Chinese Science Bulletin,49(16) : 1554_1569 (in Chinese). 
     Xin H B and Qu X M. 2008. Relative oxidation states of ore_bearing porphyries in ferred from Ce (Ⅳ) /Ce (Ⅲ) ratio in zircon: Application to the porphyry c opper belt at Gangdese, Tibet[J]. Acta Mineralogica Sinica,28(2) : 152_16 0 ( in Chinese with English abstract).
     Xu D R, Wang Z L, Cai J X, Wu C J, Nonna B C, Wang L, Chen H Y, Baker M J and Ku siak M A. 2013. Geological characteristics and metallogenesis of the shilu Fe_or e deposit in Hainan Province, South China[J]. Ore Geology Reviews, 53:318_342. 
     Xu D R, Wang Z L, Chen H Y, Hollings P, Jansen N H, Zhang Z C and Wu C J. 2014. Petrography and geochemistry of the Shilu Fe_Co_Cu ore district, South China: Im plications for the origin of a Neoproterozoic BIF system [J]. Ore Geology Revi ews, 57:322_350.
     Yu L H, Wang G J and Wang Y J. 2010. Study on molybdenum wall_rock alteration an d mineralization of Luokuidong Mo ore in Hainan Province [J]. Non_ferrousminin g and Metallurgy, 26(01):9_12 (in Chinese with English abstract).
     Yu L H. 2014. Molybdenum ore body geological characteristics and prospecting mar ks of Luokuidong Mo ore in Hainan Province [J]. Heilongjiang Science and Techn ology Information, (28):90_92 (in Chinese).
     Zhang J J, Wang G J, Yang X Y, Sun W D and Dai S Q. 2012. The petrogenesis of th e Jingde granodiorite and its MMEs: Constraints from geochemistry,zircon U_Pb d ating and Hf isotopic compositions[J]. Acta Petrologica Sinica,28(12 ): 4 047_4063(in Chinese with English abstract).
     Zhou Y, Liang X, Liang X, Jiang Y, Wang C, Fu J and Shao T. 2015. U_Pb geochrono logy and Hf_isotopes on detrital zircons of Lower Paleozoic strata from Hainan I sland: New clues for the early crustal evolution of southeastern South China[J ]. Gondwana Research, 27(4):1586_1598.
     Zhao Z H. 2010. Trace element geochemistry of accessory minerals and its applica tions in petrogenesis and metallogenesis[J]. Earth Science Frontiers, 17(1): 267_286 (in Chinese with English abstract).    
     附中文参考文献    
     陈沐龙, 符策锐, 庄有光, 云平. 2002. 琼中地区首次发现热液石英脉型钼矿床 [J]. 华南地质与矿产, (1): 9.
     陈沐龙, 马昌前, 吕昭英, 云平, 刘园园. 2014. 海南岛千家复式岩体锆石U_Pb年代学及其 地质意义[J]. 地质科技情报, (6): 1_10+27.
     陈哲培. 1997. 海南省岩石地层[M]. 中国地质大学出版社. 1_100.
     付王伟, 许德如, 傅杨荣, 杨昌松, 周迎春, 杨东生, 吴传军, 王智琳. 2011. 海南岛钼矿 床区域成矿规律初探[J]. 矿物学报, (S1): 772_773.
     付王伟. 2013. 海南岛钼矿床成因类型及成矿动力学背景研究[D]. 导师:许德如. 中国 科学院大学. 184页.
     付王伟, 许德如, 傅杨荣, 吴传军, 杨昌松, 周迎春, 王智琳. 2013. 海南省红门岭钼钨矿 床辉钼矿Re_Os同位素定年及地质意义[J]. 东华理工大学学报(自然科学版), (2): 135_1 42.
     付王伟, 许德如, 吴传军, 傅杨荣, 周迎春, 周岳强, 王智琳, 林舸. 2014. 海南省高通岭 钼矿床赋矿岩体LA-ICP-MS锆石U_Pb定年及成矿意义[J]. 矿床地质, (2): 419_427.
     高晓英, 郑永飞. 2011. 金红石Zr和锆石Ti含量地质温度计[J]. 岩石学报, (2): 417_43 2.
     贾小辉, 王强, 唐功建, 姜子琦, 赵振华, 杨岳衡, 王晓地, 赵武强. 2010. 海南屯昌早白 垩世晚期埃达克质侵入岩的锆石U_Pb年代学、地球化学与岩石成因[J]. 地球化学, (6): 497_519.
     江思宏, 聂凤军, 张义, 胡朋. 2004. 浅成低温热液型金矿床研究最新进展[J]. 地学前 缘, (2): 401_411.
     雷玮琰, 施光海, 刘迎新. 2013. 不同成因锆石的微量元素特征研究进展[J]. 地学前缘, (4): 273_284.
     李孙雄, 陈沐龙, 杨东生, 汪焰华. 2014. 海南岛钼矿床Re_Os年龄及其成矿地球动力学背 景探讨[J]. 华南地质与矿产, (3): 272_279.
     梁新权. 1995. 海南岛前寒武纪花岗岩_绿岩系Sm_Nd同位素年龄及其地质意义[J]. 岩石 学报, (1): 71_76.
     廖香俊, 王平安, 丁式江, 黄香定, 董法先, 刘晓春, 雷伟志. 2005. 海南岛主要成矿系列 与矿床成矿规律研究[J]. 地质力学学报, (2): 187_194.
     廖香俊, 王平安, 覃海灿, 路西坤, 董法先, 刘晓春, 舒斌. 2008. 海南屯昌地区高通岭钼 矿床的地质、地球化学特征及成矿时代[J]. 地质通报, (4): 560_570.
     凌洪飞. 2011. 论花岗岩型铀矿床热液来源——来自氧逸度条件的制约[J]. 地质论评, ( 2): 193_206.
     刘君, 曹玉莲, 于立红, 王国君, 王英杰. 2010. 海南省罗葵洞钼矿围岩蚀变找矿标志[J ]. 地质与资源, (1): 63_66.
     龙文国, 林起玉, 范庆贺, 谢自, 庄有光, 杨志强. 2003. 海南琼海烟塘梅岭铜钼矿床地 质特征[J]. 华南地质与矿产, (3): 30_36.
     毛景文, 谢桂青, 郭春丽, 袁顺达, 程彦博, 陈毓川. 2008. 华南地区中生代主要金属矿床 时空分布规律和成矿环境[J]. 高校地质学报, (4): 510_526.
     孙卫东, 凌明星, 汪方跃, 丁兴, 胡艳华, 周继彬, 杨晓勇. 2008. 太平洋板块俯冲与中国 东部中生代地质事件[J]. 矿物岩石地球化学通报, (3): 218_225.
     涂湘林, 张红, 邓文峰, 凌明星, 梁华英, 刘颖, 孙卫东. 2011. RESOlution激光剥蚀系统 在微量元素原位微区分析中的应用[J]. 地球化学, (1): 83_98.
     王国君, 刘君, 曹玉莲, 张馨宁. 2010. 海南罗葵洞斑岩型钼矿地质特征及矿床成因[J]. 矿产勘查, (05): 453_457.
     吴元保, 郑永飞. 2004. 锆石成因矿物学研究及其对U_Pb年龄解释的制约[J]. 科学通报, (16): 1589_1604.
     辛洪波, 曲晓明. 2008. 西藏冈底斯斑岩铜矿带含矿岩体的相对氧化状态:来自锆石Ce(Ⅳ)/ Ce(Ⅲ)比值的约束[J]. 矿物学报, (2): 152_160.
     于立红, 王国君, 王英杰. 2010. 海南省罗葵洞钼矿围岩蚀变与矿化研究[J]. 有色矿冶, (1): 9_12.
     于立红. 2014. 海南罗葵洞钼矿矿体地质特征及找矿标志[J]. 黑龙江科技信息, (28): 9 0_92.
     张俊杰, 王光杰, 杨晓勇, 孙卫东, 戴圣潜. 2012. 皖南旌德花岗闪长岩与暗色包体的成因 :地球化学、锆石U_Pb年代学与Hf同位素制约[J]. 岩石学报, (12): 4047_4063.
     赵振华. 2010. 副矿物微量元素地球化学特征在成岩成矿作用研究中的应用[J]. 地学前 缘, (1): 267_286.