岩石学报  2012, Vol. 28 Issue (2): 557-570   PDF    
黑龙江省乌拉嘎金矿赋矿花岗闪长斑岩锆石 U-Pb年龄、岩石成因及其地质意义
王永彬1, 刘建明1, 孙守恪1, 李艳2, 李凤友3, 胡海涛3     
1. 中国科学院地质与地球物理研究所矿产资源重点实验室,北京 100029;
2. 中国石油渤海钻探工程有限公司第二录井分公司,沧州 062552;
3. 黑龙江乌拉嘎黄金矿业有限责任公司,伊春 153221
摘要: 黑龙江省是我国著名金矿产区之一,发育多个浅成低温热液型金矿,它们与中生代陆相火山-次火山岩有密切的关系,但对这些陆相火山-次火山岩的成因缺乏系统研究,制约了本区金矿的成因认识和矿床勘查。本次研究采用先进的LA-ICP-MS锆石定年法、地球化学Sr-Nd-Pb同位素示踪等方法,对与乌拉嘎浅成低温热液金矿存在密切成因联系的次火山岩(花岗闪长斑岩)进行深入研究。LA-ICP-MS锆石定年法获得乌拉嘎矿区葡萄沟岩体及其南部含矿岩枝的成岩年龄分别为108.2±1.2Ma和106±1.1Ma,与区内宁远村组火山岩成岩时间基本相近,推断金矿成矿时代为早白垩世晚期,与东安金矿和高松山金矿为同期。岩石地球化学确定该岩体为高钾钙碱性(σ=1.83~2.18)偏铝质I型花岗岩特征,结合微量元素和Sr-Nd-Pb研究显示其具有活动陆缘弧岩浆岩特点,进一步得出岩浆源区和成矿物质具有来源于新元古代形成的镁铁质下地壳的部分熔融的属性。早白垩世时,中国东部处于伸展应力体制下,尤其黑龙江构造活动强烈,起源于新生下地壳重融的岩浆活动频繁,结合已有的浅成低温热液金矿床资料显示,推断区内具有巨大的成矿潜力。
关键词: 花岗闪长斑岩     锆石U-Pb     岩石成因     乌拉嘎金矿     黑龙江省    
Zircon U-Pb geochronology, petrogenesis and geological implication of ore-bearing granodiorite porphyry in the Wulaga gold deposit, Heilongjiang Province
WANY YongBin1, LIU JianMing1, SUN ShouKe1, LI Yan2, LI FengYou3, HU HaiTao3     
1. Key Laboratory of Mineral Resources,Institute of Geology and Geophysics,Chinese Academy of Sciences,Beijing 100029, China;
2. Second Mud Logging Company of CNPC Bohai Drilling Engineering Company Limited,Cangzhou 062552, China;
3. Wulaga Gold Mining Corporation,Yichun 153221, China
Abstract: Heilongjiang is one of China's most famous gold producing districts. Past production in the area was partially from epithermal gold mineralization, which was associated with the Mesozoic subaerial volcanic to subvolcanic intrusions. However, a little work has been done for the petrogenesis of these igneous rocks. In this paper, LA-ICP-MS zircon U-Pb ages, coupled with whole-rock major, trace elements and Sr-Nd-Pb isotopic data are present to constrain the geochronology and the petrogenesis of ore-bearing granodiorite in the Wulaga gold deposit, also discuss the epithermal gold mineralization potential of Heilongjiang Province. Two U-Pb zircon ages from Putaogou granodiorite intrusion and its southern ore-bearing granodiorite porphyry apophysis that occur in the Wulaga gold district are respectively 108.2±1.2Ma and 106±1.1Ma, which are similar to the age of surrounding volcanic rocks. Thus it is inferred that the gold mineralization hosted within the granodiorite porphyry apophysis was also formed in the Early Cretaceous, which was contemporaneous with the Dong'an gold deposit and Gaosongshan gold deposits. The granodiorite porphyry is the typically I-type, calc-alkaline (σ=1.83~2.18) granitoid with meta-aluminous attributes (A/NCK=0.97~1.16). Geochemical and Sr-Nd-Pb isotopic features suggest that the granodiorite porphyry was likely derived from partial melting of the Neoproterozoic juvenile mafic low-crust, which is inferred to be the source of partial mineral-forming components by comparing with other gold deposits. Early Cretaceous magmatism favored an extensional regime in NE China, when complex tectonic activities occurred, at the same time as intense low crust-sourced magmatism. This result suggests that there is likely huge mineralization potential for epithermal gold deposits in the Heilongjiang Province, through comparison with regional mining epithermal gold deposits.
Key words: Granodiorite porphyry     Zircon U-Pb     Petrogenesis     Wulaga gold deposit     Heilongjiang Province    
1 引言

自二十世纪七八十年代,浅成低温热液金矿成矿模型的建立,带动了全世界新一轮寻找同类型金矿的热潮,在中国东部中生代陆相火山-次火山岩中陆续发现了多个同类型的金矿床(如乌拉嘎金矿等),这为深入研究我国发育的浅成热液成矿系统提供了良好的研究实例。大量研究表明,浅成低温热液金矿与火山-次火山岩具有密切时空联系(Sillitoe and Bonham,1984)。岩浆活动不仅为金矿形成提供主要的热源,而且还提供一定量的成矿流体和成矿物质(Hedenquist and Houghton,1987White and Hedenquist,1990Giggenbach,1996;)。因此,深入研究与金矿化作用密切相关的火山-次火山岩,对深化中国东部中生代浅成低温热液成矿系统的认识与指导区域找矿都具有重要意义。

乌拉嘎金矿(又称团结沟金矿)作为我国浅成低温热液型金矿典型代表之一,受到众多学者关注,在矿床地质特征(吴尚全,1995)、成矿流体(孙丰月等,2008)和矿床成因(来又东,2008聂喜涛,2010)等方面获得大量成果,但前人对与成矿关系密切的次火山岩缺乏深入系统的岩石成因研究,制约了成矿物质来源及其矿床勘探的研究。因此,本次研究选择乌拉嘎矿区赋矿次火山岩体作为研究对象,在详细野外地质观察基础上,通过对采集样品的精细年代学、岩石地球化学、同位素示踪等研究,结合区域资料与国内外相关研究进展,深入探讨本区次火山岩的成因类型、源区性质及其与金矿化之间的联系。

2 成矿地质背景和矿床地质特征

乌拉嘎金矿区大地构造位置独特,位于中亚造山带(天山-兴蒙造山带)的东段,佳木斯地块与松嫩地块的衔接处北端(图 1a),同时又处于滨太平洋构造域(Nei,1991; Şengör AMC and Natal'In,1996; Jahn et al.,2000a; Li,2006)。自古生代到中生代先后受到了古亚洲洋和古太平洋演化的影响,具有多块体、多阶段的演化特征(唐克东等,1995; 李锦轶,1998; Li,2006; Wu et al.,2007; Zeng et al.,20092011ab)。

图 1 区域构造示意图与矿区地质图 (a)-东北地区构造格架图(据Gu et al.,2007);(b)-乌拉嘎金矿矿区地质图(据聂喜涛,2010);(c)-乌拉嘎金矿102线地质剖面图 Fig. 1 Regional tectonic sketch and local geological map (a)-tectonic framework of northeastern China(after Gu et al.,2007);(b)-simplified geological map of Wulaga gold deposit(after Nie,2010);(c)-geological cross-section of No.102 exploration line of Wulaga gold deposit

矿区构造主要表现为脆性断裂,呈现多期活动的特点,古老的东西向断裂可能在中生代复活并被改造,NNE向的乌拉嘎深断裂旁侧发育着一系列NWW、NNW、NEE向次级断裂,不同程度地控制本区的成岩成矿作用。

矿区地层简单,呈双层结构,下伏的元古代黑龙江群,为一套绿片岩相到角闪岩相变质岩,主要由斜长角闪岩、绿泥片岩、云母石英片岩等各种中深成相变质岩组成。90年代,黑龙江群中断断续续发现蓝片岩、含放射虫化石的硅质岩、含几丁虫化石的千枚岩等,认为黑龙江群为一套增生碰撞杂岩(Wu et al.,2007),并识别出洋壳残片,对其定年研究,确认存在777±18Ma的结晶年龄,并受到437±7Ma的变质作用影响(颉颃强等,2008)。

上覆地层为晚中生代宁远组火山岩,为陆相火山碎屑岩建造和火山熔岩堆积,其岩性为中酸性(安山质-英安质-流纹质)火山熔岩和碎屑岩为主,并以夹有珍珠岩、黑耀岩为特征(来又东,2008);宁远村组(K1n)的前人测得全岩K-Ar同位素年龄分别为118Ma、109Ma及100.1Ma(曲关生,2008)。

矿区浅成侵入岩分布于乌拉嘎断陷东侧的葡萄沟、张才沟、团结沟一带,有数个岩体出露,主要呈岩株、岩舌状,具有燕山晚期的火山岩浆沿嘉荫断陷东缘断裂乌拉嘎段的火山管道上升至近地表实现侵位的次火山岩相特征。乌拉嘎含矿岩体呈岩株状,矿区钻孔推断其可能是葡萄沟岩体南部的舌状分枝。本次研究确认两者岩性相似,均为花岗闪长斑岩。岩相学研究确认:主要造岩矿物由斜长石(35%~55%)、石英(20%~25%)、正长石(10%~20%)及少量黑云母(3%~5%)组成,副矿物为锆石、磷灰石和磁铁矿等。以含有宽大的斜长石和石英斑晶的斑状结构为特征,其中斜长石斑晶为中长石,具有典型的环带结构,普遍发育聚片双晶,石英斑晶多被熔蚀成港湾结构(如图 2d),明显的高温浅成相产物。

图 2 乌拉嘎金矿床矿石标本(a)、花岗闪长斑岩(b)及花岗闪长斑岩显微照片(c、d) Qz-石英;Or-正长石;Pl-斜斜长石;Bt-黑云母 Fig. 2 Photographs of ore samples(a),granodiorite porphyry sample(b),microphotographs of granodiorite porphyry(c,d)in Wulaga gold deposit Qz-quartz; Or-orthoclase; Pl-plagioclase; Bt-biotite

乌拉嘎金矿化范围波及早白垩世火山岩、花岗闪长斑岩及中元古界黑龙江地层,但是具有工业意义的矿体则主要赋存于花岗闪长斑岩上部的隐爆角砾岩带(如图 1b,c)。因金矿化不均匀,变化性较大,矿体与围岩无明显界线,矿体形态和大小只能根据边界品位圈定,总体呈NWW-NW走向的雁形排列,倾向北北东,倾角20°~45°。金属矿物以黄铁矿、白铁矿和辉锑矿为主,硅化和绢云母化等蚀变与矿体紧密共生。

前人(孙凤兴等,1996魏仪方和刘春华,1996)研究中,矿化与岩浆活动具有密切的成因关系,其成矿物质和成矿流体的来源,均有岩浆作用的迹象,且岩体成矿元素分析,发现中生代火山岩-次火山岩本身金的丰度较高,近乎数倍于地壳克拉克值,且其标准离差和变化系数较小,反映其中的金分布比较均匀且集中,是可能的矿源岩,因此深入探讨含矿次火山岩,对金矿的研究具有重要意义。

3 测试方法 3.1 LA-MC-ICP锆石定年

本次研究所用锆石,来自两个样品,即:WLG-W-1(乌拉嘎含矿岩枝)和PTG-1(北部葡萄沟岩体)。锆石的分选工作由河北省区域地质调查大队地质实验室完成:样品粉碎后,通过浮选和磁选获得锆石(>1000粒);样品靶的制备按照SHRIMP定年锆石样品的制备方法进行(宋彪等,2002),在中国科学院地质与地球物理研究所IMS-1280实验室制靶中心完成:在双目镜下挑选出晶形和透明度较好的锆石颗粒,然后将它们粘贴在环氧树脂表面,抛光后将锆石进行透射光、反射光和阴极发光显微照相。锆石阴极发光(CL)显微照相在中国科学院地质与地球物理研究所的电子探针与电镜实验室,通过德国LEO1450VP扫描电子显微镜(SEM)完成(电压15kV;电流1.1nA);锆石定年和微量元素分析在中国科学院地质与地球物理研究所多接受等离子体质谱仪实验室(MC-ICPMS Laboratory)完成:激光剥蚀系统为GeolasPlus型193nm ArF准分子激光器,ICP-MS设备为Agilent7500a型Q-ICPMS。实验室的激光输出能量密度在15J/cm2的情况下,激光束斑大小为4~160μm调整,激光脉冲速率为2~10Hz。实验室采用He气作为剥蚀物质的载气,Ar气用来增加传输速率,由于采用高纯度的Ar和He气,204Pb和202Hg的背景值大多<100cps,数据采集模式为Time-resolved Analysis,选用一个质量峰采集一点的跳峰方式(peak hopping,one point per peak),每8个样品分析点分别测定一组标样(两个标准锆石91500和GJ-1,一个微量元素标样NIST SRM 610)。每个分析点的气体背景采集时间为30s,信号采集时间为40s,207Pb/206Pb,206Pb/238U,207Pb/235U(235U=238U/137.88),208Pb/232Th的比值采用标准锆石GJ-1为外部标准进行校正(除了考虑样品和外标锆石GJ-1同位素比值在测定过程中产生的标准偏差外,GJ-1相对偏差设定为2%)。分馏校正和结果计算采用GLITTER4.0软件处理,普通铅校正采用Anderson(2002) 的方法,年龄的处理采用Ludwig(2003) 的ISOPLOT3.0程序。元素浓度采用GLITTER(4.0)程序计算,以NIST SRM610作外标,以Si作为内标。详细的分析流程参考Xie et al.(2008)

3.2 主微量测试

挑选新鲜的岩石样品,由河北省区域地质调查大队地质实验室粉碎成200目粉末。岩石地球化学主量元素、微量元素及稀土元素的分析在核工业北京地质研究院分析测试:主量元素分析采用PHILIPS PW-2404 X-荧光光谱仪,详细测试流程见Norrish and Hutton(1969) ,样品先按照1∶5的比例放入Li2B4O7溶液中,在1050~1250℃温度下熔化,然后将熔化样品制成玻璃薄片进行分析,分析精度估计1%(SiO2)和2%(其它氧化物)。微量元素和稀土元素的分析采用ELEMENT-2型等离子质谱仪,详细的测试流程参考Qi et al.(2000) ,将50mg全岩粉末(200目)置于Teflon烧瓶中分别用HNO3和HF溶解2d后加入HClO4进一步溶解。蒸干后用5%的HNO3溶液将样品稀释到50mL,最后做微量元素分析。分析过程中同时测定两个GBW系列标样:一个标样校正Li,Ba,V,Cr,Co,Ni,Cu,Zn,Ga,Rb,Sr,Cs,Ba,Pb,Th,U,Sc,Y及稀土元素;另一个标样校正W,Mo,Nb,Ta,Zr和Hf等元素。分析误差<5%。

3.3 Sr-Nd-Pb同位素测试

Sr-Nd-Pb同位素分析在核工业北京地质研究院分析测试研究中心完成。根据微量元素中Rb-Sr-Sm-Nd含量称取适量样品于Teflon 闷罐中,加入87Rb-84Sr和149Sm-150Nd 混合稀释剂并用HF、HNO3和HClO4充分溶解后用离子交换树脂分离出Rb、Sr、Sm和Nd,最后在ISOPROBE-T热电离质谱仪(TIMS)上测试。整个分析流程实验本底为:Rb,Sr<100×10-12;Sm,Nd<50×10-12

称取约120mg样品于Teflon闷罐中,加入HF和HNO3,在100℃电热板上加热一周充分溶解后用离子交换树脂将Pb分离出来,同样在ISOPROBE-T热电离质谱仪上测试。Pb的全流程本底小于50×10-12

4 实验测试结果 4.1 锆石LA-ICP-MS 定年 4.1.1 锆石特征

锆石广泛存在于各类岩石中,富含U和Th,低普通Pb以及非常高的矿物稳定性,使得锆石U-Pb定年成为同位素年代学研究中最常用和最有效的方法之一(吴元保和郑永飞,2004),但锆石可以形成于岩浆结晶过程、岩浆期后热液活动过程中,并在高温高压变质过程中可以增生,因此对其成因的判断尤其重要。岩浆锆石的判断标准包括:(1) 阴极发光照像(CL)显示典型振荡环带结构(简平等,2001);(2) Th/U比值多为0.1~1.0(Belousova et al.,2002),一般大于0.5(Hoskin and Schaltegger,2003);(3) 岩浆锆石REE含量和(La/Yb)N更低,δCe更高(Hoskin et al.,2005)。

本次研究的锆石,来自两个样品(图 3图 4),但是其特征具有相似性。锆石无色透明,晶型较好,多为长柱状或短柱状,长100~200μm,长宽比介于1.5∶1~3∶1。阴极发光下具有清晰均匀的震荡环带,未发现增长边或核幔结构。样品中锆石Th、U含量分别为:55.01×10-6~286.4×10-6;206.2×10-6~581.4×10-6,Th/U比值均大于0.1(0.2~0.54之间),Th、U之间正相关性较好,数据见表 1,同时利用LA-ICP-MS测试的锆石稀土数据作图(图 3图 4),两个样品表现具有一致性,其REE含量较高(526×10-6~833×10-6),富集重稀土,强烈亏损轻稀土,其(La/Yb)N比值近乎0,但铈异常明显较高,δCe值达到4.33~156。综合以上特征本次选测的锆石为典型的岩浆锆石,其年龄可以代表其结晶年龄。

图 3 乌拉嘎金矿含矿花岗闪长斑岩岩枝中锆石的形态、分析点位图及锆石稀土模式图(球粒陨石标准值据Sun and Mcdonough,1989) Fig. 3 Cathodoluminescence(CL)image and chondrite-normalized REE diagrams of zircon from ore-bearing granodiorite porphyry apophysis in Wulaga gold deposit(chondrite REE values after Sun and Mcdonough,1989)

图 4 乌拉嘎金矿葡萄沟岗闪长斑岩中锆石的形态、分析点位图及锆石稀土模式图(球粒陨石标准值据Sun and Mcdonough,1989) Fig. 4 Cathodoluminescence(CL)image and chondrite-normalized REE diagrams of zircon from ore-baring granodiorite porphyry in Wulaga gold deposit(chondrite REE values after Sun and Mcdonough,1989)
4.1.2 锆石年龄

样品锆石LA-ICP-MS U-Pb年龄数据列于表 1。本次测定有效外标GJ-1锆石数据共20个点,其中样品WLG-W-1,附带标样6个,表面年龄206Pb/238U: 595.8~610Ma,加权平均值为603.7Ma;样品PTG-1,附带标样7个,表面年龄206Pb/238U :597.9~612.4Ma,加权平均值为607.6Ma。上述标样年龄值,均在标准值603.12±2.4Ma(柳小明等,2007)的误差(±5%)范围内,本次测量数据真实可靠。

表 1 乌拉嘎金矿花岗闪长斑岩LAI-CPMS锆石U-Pb分析结果 Table 1 Zircon U-Pb LA-ICPMS data for granodiorite porphyry Wulaga gold deposit

样品WLG-W-1,测定了19个有效单颗粒锆石,19个测点锆石U-Pb同位素组成在误差范围内非常谐和,均落在谐和线上,206Pb/238U加权平均年龄为106±1.1Ma(可信度95%,MSWD=0.61)(图 5)。

图 5 乌拉嘎金矿含矿花岗闪长斑岩岩枝(WLG-W-1)锆石年龄谐和图 Fig. 5 Concordia diagram of the ore-bearing granodiorite porphyry apophysis(WLG-W-1)in Wulaga gold district

样品PTG-1,测定了19个有效单颗粒锆石,19个测点锆石U-Pb同位素组成在误差范围内非常谐和,均落在谐和线上,206Pb/238U加权平均年龄为108.2±1.2Ma(可信度95%,MSWD=0.82)(图 6)。

图 6 葡萄沟花岗闪长斑岩(PTG-1)锆石年龄谐和图 Fig. 6 Concordia diagram of granodiorite porphyry(PTG-1)in the Putaogou area of the Wulaga gold district
4.2 岩石地球化学特征 4.2.1 主量元素特征

本次研究选择了6个新鲜样品进行分析,数据列表见表 2。地球化学特征具有相似性:SiO2的含量为64.34%~65.62%,属中酸性岩范畴;Al2O3=13.85%~16.21%,A/NCK为0.97~1.16,多数值落在1~1.1之间,为偏铝质-弱过铝质花岗岩特征;Na2O=3.29%~3.95%,K2O=2.57%~3.28%,Na2O+K2O含量高,平均为6.62%,富钠(Na2O/K2O>1),在SiO-K2O联合图解上,样品落入高钾钙碱性区域,与计算的里特曼指数[(=(Na2O+K2O)2/(SiO2-43)]一致,为钙碱性((=1.85~2);相对富钙(2.48%~3.12%)和富镁(1.37%~2.01%),镁指数较高[Mg#=Mg2+/(Mg2++0.9×Fe2O3T/72)=0.49~0.52]。另外,P2O5含量低且变化小(0.11%~0.15%),低TiO2(0.46%~0.76%),低FeOT。岩石主量数据在TAS和R1-R2投图(图 7),结果相似,位于亚碱性的花岗闪长岩区间,与镜下岩相学定名具有统一性。

表 2 乌拉嘎金矿含矿花岗闪长斑岩主量元素(wt%)、稀土元素及微量元素(×10-6)分析结果 Table 2 The analyzed data of major(wt%),rare earth and trace elements(×10-6)for ore-bearing granodiorite porphyry in Wulaga gold deposit

图 7 乌拉嘎金矿含矿花岗闪长斑岩主量元素特征 a)-R1-R2图解;(b)-TAS图解(底图据Middlemost,1985);(c)-硅钾图(底图据Rickwood,1989);(d)-铝饱和指数图(底图据 Maniar and Piccoli,1989) Fig. 7 Characteristics of major elements of granodiorite porphyry in Wulaga gold deposit (a)-R1-R2 diagram;(b)-TAS diagram(after Middlemost,1985);(c)-SiO2-K2O diagrams(after Rickwood,1989);(d)-discriminant diagram of A/NK-A/CNK(after Maniar and Piccoli,1989)
4.2.2 稀土及微量元素特征

乌拉嘎金矿含矿岩体的稀土含量变化不大,配分模式几乎平行(表 2图 8)。岩体的稀土含量整体较低,变化为78.8×10-6~134.3×10-6,但样品的稀土配分模式相似的表现为轻稀土富集(∑LREE=72.18×10-6~122.8×10-6),重稀土亏损(∑HREE=6.62×10-6~11.8×10-6)的右倾模式,轻重稀土分馏明显((La/Yb)N=13.19~16.72),无明显铕异常(δEu=0.82~1.05),暗示源区可能有石榴石或角闪石矿物残留。微量元素(图 8)组成上,亏损Nb、Ta、P、Zr、Hf和Ti等高场强元素,富集Rb、K和Sr等大离子亲石元素,强烈富集U和Pb,弱亏损Th,高Sr含量(278×10-6~323×10-6),低Y含量(9.04×10-6~17.4×10-6),中等Sr/Y比值(16.73~35.73,均值24.61);Rb/Sr比值(0.24~0.44)变化不大;稀土配分模式和微量元素蛛网图综合分析,具有活动陆缘弧花岗岩特点。

图 8 乌拉嘎金矿含矿花岗闪长斑岩稀土元素配分曲线(a)及微量元素蛛网图(b)(标准化值来自Sun and Mcdonough,1989) Fig. 8 Chondrite-normalized REE(a)and chondrite-normalized trace element spider diagrams(b)of ore-barren granodiorite porphyry in Wulaga gold deposit(normalization values after Sun and Mcdonough,1989)
4.2.3 Sr-Nd-Pb同位素特征

本次研究中乌拉嘎金矿成矿母岩-葡萄沟岩体(编号PTG)的花岗闪长斑岩的Sr-Nd-Pb同位素组成见表 3表 4图 9图 10。基于先前的同位素年龄数据,取108Ma的平均年龄代表乌拉嘎岩体的侵位时代,并据此对有关参数进行计算。在计算Nd模式年龄时,采用两阶段模式。所采用的参数见表 3,计算公式引自Li and McCulloch(1996) 以及沈渭洲等(1999)

表 3 乌拉嘎金矿含矿花岗闪长斑岩Sr-Nd同位素分析结果及主要参数 Table 3 The Sr-Nd isotopic determination data and some major parameters for ore-bearing granodiorite porphyry in Wulaga gold deposit

表 4 乌拉嘎金矿含矿花岗闪长斑岩Pb同位素分析结果及主要参数 Table 4 The Pb isotopic determination data and some major parameters for ore-bearing granodiorite porphyry in Wulaga gold deposit

其(87Sr/86Sr)i变化不大,介于0.705209~0.706178,平均值为0.705586;143Nd/144Nd初始值较高(介于0.512505~0.512533),平均为0.512526,εNd(t)绝对值较小,为0.1~1.1,该值远高于华北克拉通古老下地壳的εNd(t)(-44~-32;Jahn et al.,1999),与兴蒙造山带显生宙花岗岩的εNd(t)值(-2.22~7.11,平均值为+2.0;吴福元等,1999洪大卫等,20002003)一致。亏损地幔模式年龄tDM2较为年轻,其值介于800~897Ma之间。研究区花岗岩的fSm/Nd均为负值,数据比较集中(介于-0.46~-0.41),表明源区的稀土元素Sm、Nd分馏不明显,因此运用上述测试数据计算的钕二阶段同位素模式年龄是可靠的。

本次测试选择花岗岩具有较高的放射性成因铅组成和较低的非放射性成因铅204Pb。(206Pb/204Pb)i=18.232~18.387,平均为18.315;(207Pb/204Pb)i=15.51~15.537,平均值为15.519;( 208Pb/204Pb)i=38.08~38.254,平均值为38.155。由铅同位素比值计算出的μ值介于9.3~9.34之间,低于μ=9.74的陆壳演化线。上述Pb同位素组成和大兴安岭构造同位素地球化学省Pb同位素特征(张理刚,1995)相似。

综上所述,乌拉嘎金矿岩体具有低锶初始值,较高的钕初始值,εNd(t)略高于0,钕同位素亏损地幔二阶段模式年龄较年轻,岩石富集放射性铅组分,而非放射性成因铅204Pb含量较低,总体特征显示兴蒙造山带显生宙花岗岩特点。

5 讨论 5.1 成岩时代

乌拉嘎金矿矿体赋存于燕山晚期的花岗闪长斑岩与中元古界黑龙江群下亚群接触带附近的近东西向展布的隐爆角砾岩中,金矿化主要受隐爆角砾岩带控制。成矿与矿区内花岗闪长斑岩有密切的联系,斑岩年龄的确定是对成矿时代有效的制约。前人曾对此进行了各种方法的定年尝试:斑岩锆石同位素年龄为105Ma、114Ma(江雄新,1986);全岩K/Ar同位素年龄为100Ma、102Ma、112.6Ma(吴尚全,1995);锆石LA-ICP-MS同位素年龄109Ma(聂喜涛,2010),基本确定岩体为早白垩世岩浆浅成侵入产物。但是受限于当时的定年方法和测试精度,使得斑岩年龄波动较大(100~114Ma)。

本次研究对矿区出露的两个岩体(WLG-W-1;PTG-1)分别采样,选择锆石进行系统定年,通过精确的LA-ICP-MS定年方法,得到准确的年龄分别为(106±1.1)和(108.2±1.2)Ma,其中含矿岩体(WLG-W-1)略晚于矿区北部的无矿岩体(PTG-1),与聂喜涛(2010) 所获得的年龄几乎一致。由此验证了乌拉嘎含矿岩体为葡萄沟岩体在南部的舌状分支,且推断含矿岩体为岩浆演化晚期产物,可能富集更多的气液成分,为成矿提供了必要条件。

前人对区内存在的白垩纪宁远村组火山岩研究成果显示,其成岩时代区间为100~118Ma(全岩K-Ar同位素年龄),而聂喜涛(2010) 通过对宁远村组英安质岩屑晶屑熔结凝灰岩进行了单颗粒锆石年代学研究,获得单颗粒锆石的年龄在102~121Ma,其加权平均值为113.3±1.2Ma,为早白垩世晚期火山作用产物。而区内花岗闪长斑岩年龄为106~108Ma,比火山岩稍晚,基本为一期火山-岩浆作用产物,结合前人(孙凤兴等,1996魏仪方和刘春华,1996)研究,乌拉嘎金矿显示的亲岩浆特征,推断成矿时代为早白垩世晚期。

黑龙江省内分布大量同期的早白垩世火山-次火山岩,且伴生一定量的浅成低温热液金矿,例如产于光华组的东安金矿,其成岩和成矿年龄分别为:112Ma和108Ma(敖贵武等,2004);产于甘河组(美丰组)的高松山金矿(刘桂阁等,2006边红业等,2009尹西君等,2010)等。早白垩世火山岩中浅成低温热液金矿的发育,尤其突显早白垩世火山-侵入作用的重要,本次研究的岩体正可作为其中典型代表,深入探讨其背景,有助于加深区域矿床的研究和勘查。

5.2 成因类型

花岗岩的分类历来都是学者关注和研究的热点,当下最常用的花岗岩成因分类方案是MISA(即M、I、S和A型)(吴福元等,2007)。A型和M型具有典型成因特征,本次研究涉及的岩体明显不符,所以在此仅做I型和S型探讨即可。Chappell and White(2001) 在I型和S型花岗岩分类提出25周年之际,再次总结了两类花岗岩的分类特征:相比而言,前者具有富钠、富钙、富暗色矿物(角闪石等),且铝饱和指数低(A/NCK<1.1),少过铝质矿物等。

通过对乌拉嘎金矿有关的花岗闪长斑岩的地球化学特征研究发现,该岩体为中酸性(SiO2=64.34%~65.62%),钙碱性岩石系列(σ=1.85~2),全碱含量较高(6.52%~7.31%):富钠(Na2O/K2O>1),相对富钙(2.48%~3.12%)和富镁(1.37%~2.01%),镁指数略高(Mg#=0.49~0.52),但P2O5含量低且变化小(0.11%~0.15%),铝饱和指数(A/NCK)不高,为0.97~1.16,多数值落在1~1.1之间,为偏铝质-弱过铝质花岗岩特征。以上地球化学特征均显示其为典型的I型花岗岩。

同期含矿的火山岩(宁远村组、光华组和甘河组)以中酸性安山-流纹岩为主,杂生基性火山岩(尹西君等,2010),地球化学研究显示多为偏铝-弱过铝质的钙碱性火山岩。稀土元素分析,配分曲线与乌拉嘎花岗闪长斑岩相似,均呈重稀土富集,轻稀土亏损的右倾模式,其中宁远村组和甘河组安山质火山岩没有明显的铕异常,而光华组与矿有关的潜流文岩具有明显铕异常(于建波,2006),铕异常的差别可能源于样品的不同,潜流纹岩形成于火山活动晚期,源区残留斜长石的缘故。黑龙江省白垩世火山岩地球化学性质总体类似活动大陆活动边缘火山岩,众所周知,此类火山岩石环太平洋浅成低温热液金矿的重要赋矿围岩(Mitchell and Garson,1981White and Hedenquist,1990)。

5.3 岩浆源区

前文讨论认为本次研究的花岗闪长斑岩为高钾钙碱性I型花岗岩,且轻重稀土明显分馏,亏损高场强元素(Nb、Ta、Ti、P),富集大离子亲石元素(Rb、Th、U和Pb)和微弱的Eu异常,其源岩应为壳源火成岩。但其Mg#值(均值为0.51)略高于壳源中酸性岩石的最高值(约0.5),推断其可能为下地壳底部部分熔融,或有部分地幔物质的加入。同时花岗闪长斑岩具有较低的(87Sr/86Sr)i(0.705209~0.706178)和弱富集的Nd同位素组成[εNd(t)=0.1~1.1],在εNd(t)-(87Sr/86Sr)i图解(图 9)中,样品点均位于亏损地幔和EMI型富集地幔之间,处于兴蒙造山带显生宙花岗岩的εNd(t)值区间内,明显高于世界典型壳源花岗岩(洪大卫等,2000),Wu et al.(2000) 认为εNd(t)>0的花岗质岩浆源自年轻的新生的镁铁质下地壳。初始Pb-Pb同位素组成图解(图 10)上,花岗闪长斑岩样品点投在MORB与下地壳区域交界处。一般认为由地幔部分熔融直接形成中酸性岩浆的可能性极小(Green,1980Martin,1999),因此乌拉嘎金矿赋矿围岩的源区物质应该是新生的镁铁质地壳的部分熔融,这种镁铁质下地壳可能是早期幔源岩浆底侵作用的产物,结合其具有相近且年轻的亏损地幔Nd模式年龄(tDM2=800~897Ma),推测这种年轻的地壳主要形成于新元古代之初。

图 9 乌拉嘎金矿含矿花岗闪长斑岩εNd(t)-(87Sr/86Sr)i图解(据Wu et al.,2000) Fig. 9 The εNd(t)-(87Sr/86Sr)i diagram of the granodiorite porphyry in Wulaga gold deposit(after Wu et al.,2000)

图 10 乌拉嘎金矿含矿花岗闪长斑岩初始208Pb/204Pb-206Pb/204Pb和207Pb/204Pb- 206Pb/204Pb图解(据Zindler and Hart,1986) Fig. 10 The 208Pb/204Pb-206Pb/204Pb and 207Pb/204Pb-206Pb/204Pb diagrams of granodiorite porphyry in Wulaga gold deposit(after Zindler and Hart,1986)
5.4 成矿岩浆背景

晚中生代,中国东部受控于太平洋俯冲作用,并发生了由东西向转变为北东向-北北东向构造体制的转折。至早白垩世,中国东部处在伸展应力状态下,产生了一些列的NE-NNE向的断层和左旋拉伸盆地,乌拉嘎断陷盆地就是其中之一。在西太平洋板块的俯冲作用下,牡丹江-乌拉嘎等断裂复活,诱发了大规模具有活动陆缘性质的陆相火山岩(宁远组、光华组和甘河组等)喷发和岩浆侵入(Wu et al.,2005)同时伴生了一批浅成低温热液型金矿(乌拉嘎金矿、东安金矿和高松山金矿等)。

本次研究的乌拉嘎金矿赋矿花岗闪长斑岩,其地球化学特征和Sr-Nd-Pb同位素显示,其来源可能为晚元古代新生的镁铁质下地壳,其形成多有幔源岩浆的贡献。Taylor et al.(1985)统计大陆地壳中金丰度时发现,下部大陆地壳金丰度(3.4×10-9)明显高于上部地壳(1.8×10-9)。而近年来胶东金矿研究显示,地幔作用在成矿过程中提供了部分成矿物质(毛景文等,20022005刘建明等,2003),侯增谦(2004) 认为初生的埃达克质熔体与地幔橄榄岩/幔源熔体的相互作用,可能是埃达克质熔体获取金属和硫的重要途径。由此推断与地幔作用具有紧密关系的下地壳可能是潜在的矿源层,而乌拉嘎、东安金矿和高松山金矿赋矿火山-次火山岩金的高含量(魏仪方和刘春华,1996郭继海等,2004),也侧面验证了这一推断。

从成因学角度来说,热液型矿床的定义为成矿物质被成矿热液搬运到适当的空间富集成具有经济价值的地质体。矿床的形成条件:充足的矿源层提供成矿物质,强烈的岩浆活动提供热源和热液,复杂且强烈的构造活动提供运移通道。而早白垩纪世伸展体制下,黑龙江构造活动强烈,起源于新生下地壳重融的岩浆活动频繁,结合已有的矿床资料显示,区内具有巨大的成矿潜力。

6 结论

(1) 锆石LA-ICP-MS定年获得乌拉嘎矿区葡萄沟岩体(PTG)及其南部含矿分枝(WLG)的成岩年龄分别为108.2±1.2Ma和106±1.1Ma,与区内宁远村组火山岩成岩时间基本相近,推断金矿成矿时代为早白垩世晚期,与东安金矿和高松山金矿为同期。

(2) 地球化学研究表明赋矿花岗闪长斑岩为高钾钙碱性、偏铝质I型花岗岩,其微量元素和Sr-Nd-Pb研究显示其具有活动陆缘弧岩浆特点,推断岩浆源区和成矿物质具有来源于新元古代形成的镁铁质下地壳的部分熔融的属性。

(3) 早白垩世时,中国东部处于北西向伸展应力体制下,黑龙江构造活动强烈,起源于新生下地壳(潜在矿源层)重熔的岩浆活动频繁,具有巨大的浅成热液金矿床成矿潜力。

致谢 感谢黑龙江省乌拉嘎金矿公司有关领导和技术人员为野外工作提供的帮助;锆石U-Pb定年得到了中国科学院地质与地球物理研究所杨岳衡老师的帮助;地球化学(主量、微量、Sr-Nd-Pb)测试得到了北京核工业第三研究院刘牧老师的帮助;论文写作过程中与曾庆栋老师、褚少雄、孙燕等进行了有益的交流;论文修改过程中,匿名审稿专家和编辑部俞良军老师提供了宝贵的建设性意见,使作者收益匪浅;在此一并表示衷心的谢忱!
参考文献
[] 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
[] Ao GW, Liu ZM, Xue MX, Pei ZX, Wang CS, Chen HM and Jiang YJ. 2004. Metallogenic mechanism of the Dong'an gold deposit in Heilongjiang. Gold Geology , 10 (1) :15–19.
[] Belousova EA, Griffin WL, O'Reilly SY and Fisher NL. 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
[] Bian HY, Chen M, Liu HL and Zhao CR. 2009. Geology and genesis of the Gaosongshan gold deposit in Xunke County, Heilongjiang Province. Geology and Resources , 18 (2) :91–95.
[] Chappell BW and White AJR. 2001. Two contrasting granite types: 25 years later. Australian Journal of Earth Sciences , 48 :489–499. DOI:10.1046/j.1440-0952.2001.00882.x
[] Giggenbach WF. 1986. The use of gas chemistry in delineating the origin of fluids discharged over the Taupo volcanic zone: A review. International Volcanological Congress, New Zealand. Proc. Symposium , 5 :47–50.
[] Green TH. 1980. Island arc and continent-building magmatism: A review of petrogenic models based on experimental petrology and geochemistry. Tectonophysics , 63 :367–385. DOI:10.1016/0040-1951(80)90121-3
[] Gu L, Zheng Y, Tang X, Zaw K, Delle-Pasque F, Wu C, Tian Z, Lu J, Ni P, Li X, Yang F and Wang X. 2007. Copper, gold and silver enrichment in ore mylonites within massive sulphide orebodies at Hongtoushan VHMS deposit, NE China. Ore Geology Review , 30 :1–29. DOI:10.1016/j.oregeorev.2005.09.001
[] Guo JH, Wang CS and Shi YJ. 2004. Geological and geochemical features of the Dongan gold deposit in Heilongjiang. Geology and Prospecting , 40 (4) :37–41.
[] Hedenquist JW and Houghton BF. 1987. Epithermal gold mineralization and its molcanic environments. Earth Resources Foundation, the University of Sydney :422–423.
[] Hong DW, Wang SG, Xie XL and Zhang JS. 2000. Genesis of positive εNd(t) granitoids in the Da Hinggan Mountains-Mongolia orogenic belt and growth of continental crust. Earth Science Frontiers , 7 (2) :441–456.
[] Hong DW, Wang SG, Xie XL, Zhang JS and Wang T. 2003. Correlation between continental growth and the supercontinental cycle: Evidence from the granites with positive εNd(t) in the Central Asian orogenic belt. Acta Geologica Sinica , 77 (2) :203–209.
[] Hoskin PWO and Schaltegger U. 2003. The composition of zircon and igneous and metamorphic petrogenesis. Reviews in Mineralogy and Geochemistry , 53 (1) :27–62. DOI:10.2113/0530027
[] Hoskin PWO. 2005. Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the Jack Hills, Australia. Geochimica et Cosmochimica Acta , 69 (3) :637–648. DOI:10.1016/j.gca.2004.07.006
[] Hou ZQ. 2004. Porphyry Cu-Mo-Au deposit: Some new insights and advances. Earth Science Frontiers , 11 (1) :131–144.
[] Jahn BM, Wu FY, Lo CH and Tsai CH. 1999. Crust-mantle interaction induced by deep subduction of the continental crust: Geochemical and Sr-Nd isotopic evidence from post-collisional mafic-ultramafic intrusion of the northern Dabie complex, central China. Chemical Geology , 157 :119–146. DOI:10.1016/S0009-2541(98)00197-1
[] Jahn BM, Wu FY and Chen B. 2000a. Granitoids of central Asian orogenic belt and continental growth in the Phanerozoic. Transactions of the Royal Society of Edinburgh: Earth Sciences , 91 :181–193. DOI:10.1017/S0263593300007367
[] Jahn BM, Griffin WL and Windley R. 2000b. Continental growth in the Phanerozoic: Evidence from Central Asia. Tectonophysics , 328 :7–10.
[] Jian P, Cheng YQ and Liu DY. 2001. Petrographical study of metamorphic zircon: Basic roles in interpretation of U-Pb age of high grade metamorphic rocks. Earth Science Frontiers , 8 (3) :183–191.
[] Jiang XX. 1986. The main geological characteristics and mineralization of Tuanjiegou gold deposits. In: Shenyang Institute of Geology and Mineral Resources (ed.). Collection of Gold Geological Research. Beijing: Geological Publishing House, 71-83(in Chinese)
[] Lai YD. 2008. Study on geological characteristics and enrichment regularities of gold mineralization in Wulaga gold deposit, Heilongjiang Province. Master Degree Thesis. Changchun: Jilin University :35–39.
[] Li JY. 1998. Some new ideas on tectonics of NE China and its neighboring areas. Geological Review , 44 (4) :339–347.
[] Li JY. 2006. Permian geodynamic setting of Northeast China and adjacent regions: Closure of the Paleo-Asian Ocean and subduction of the Paleo-Pacific Plate. Journal of Asian Earth Sciences , 26 :207–224. DOI:10.1016/j.jseaes.2005.09.001
[] Li XH and McCulloch MT. 1996. Secular variation in the Nd isotopic composition of Neoproterozoic sediments from the southern margin of Yangtze Block: Evidence for a Proterozoic continental collision in southeast China. Precambrian Research , 76 :67–76. DOI:10.1016/0301-9268(95)00024-0
[] Liu GG, Wang ED, Chang CJ, Jin BY and Yu WQ. 2006. Discussion on metallogeny of Gaosongshan gold deposit, Xunke County, Heilongjiang Province. Non-ferrous Mining and Metallurgy , 22 (4) :1–4.
[] Liu JM, Ye J, Xu JH, Sun JG and Shen K. 2003. C-O and Sr-Nd isotopic geochemistry of carbonate minerals from gold deposits in East Shandong, China. Acta Petrologica Sinica , 19 (4) :775–785.
[] Liu XM, Gao S, Diwu CR, Yuan HL and Hu ZC. 2007. Simultaneous in-Situ determination of U-Pb age and trace elements in zircon by LA-ICP-MS in 20μm spot size. Chinese Science Bulletin , 52 (9) :1257–1264. DOI:10.1007/s11434-007-0160-x
[] Ludwig KR. 2003. User's Manual for Isoplot 3. 00, a geochronological Toolkit for Microsoft Excel. Berkeley Geochronological Center Special Publication , 4 :25–32.
[] Lugmair GW and Harti K. 1978. Lunar initial 143Nd/144Nd: Differential evolution of the lunar crust and mantle. Earth and Planetary Science Letters , 39 :349–357. DOI:10.1016/0012-821X(78)90021-3
[] Maniar PD and Piccoli PM. 1989. Tectonic discrimination of granitoids. Geological Society of America Bulletin , 101 :635–643. DOI:10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
[] Mao JW, He Y and Ding TP. 2002. Mantle fluids involved in metallogenesis of Jiaodong (East Shandong) gold district: Evidence of C, O and H isotopes. Mineral Deposits , 21 (2) :121–128.
[] Mao JW, Li HM, Wang YT, Zhang CQ and Wang RT. 2005. The relationship between mantle-derived fluid and gold ore-formation in the eastern Shandong Peninsula: Evidences from D-O-C-S isotopes. Acta Geologica Sinica , 79 (6) :839–857.
[] Martin H. 1999. Adakitic magmams: Modern analogues of Archean granitoids. Lithos , 46 :411–429. DOI:10.1016/S0024-4937(98)00076-0
[] Middlemost EAK. 1985. Magmas and Magmatic rocks. Essex: Longman Group Ltd .
[] Mitchell AHG and Garson MS. 1981. Mineral Deposits and Global Tectonic Settings. Academic Press, London, 405-406
[] Nei S. 1991. Paleoclimatic and paleomagnetic constraints on the Paleozoic reconstructions of South China, North China and Tarim. Tectonophysics , 196 :279–308. DOI:10.1016/0040-1951(91)90327-O
[] Nie XT. 2010. Study on metallogenic geological backgroun and metallogenic model in Tuanjiegou gold deposit, Heilongjiang Province. Master Degree Thesis. Changchun: Jilin University :34–50.
[] Norrish K and Hutton JT. 1969. An accurate X-ray spectrographic method for the analysis of a wide range of geological samples. Geochimica et Cosmochimica Acta , 33 (4) :431–453. DOI:10.1016/0016-7037(69)90126-4
[] Qi L, Hu J and Conrad GD. 2000. Determination of trace elements in granites by inductively coupled plasma mass spectrometry. Talanta , 51 (3) :507–513. DOI:10.1016/S0039-9140(99)00318-5
[] Qu GS.2008. Petrology and Formation of Heilongjiang. Wuhan: Chinese Geology University Press : 1 -300.
[] Rickwood PC. 1989. Boundary lines within petrologic diagrams which use oxides of minor elements. Lithos , 22 :247–263. DOI:10.1016/0024-4937(89)90028-5
[] Shen WZ, Ling HF, Li WX, Huang XL and Wang DZ. 1999. Study on the Nd-Sr isotopic composition of granitoids in SE China. Geological Journal of China Universities , 5 (1) :22–25.
[] Şengör AMC and Natal'In BA. 1996. Paleotectonics of Asia: Fragments of a synthesis. In: Yin A and Harrison M (eds.). The Tectonic Evolution of Asia. Cambridge: Cambridge University Press, 486-640
[] Sillitoe RH and Bonham HF. 1984. Volcanic landforms and ore deposits. Econ. Geol. , 79 :1286–1298. DOI:10.2113/gsecongeo.79.6.1286
[] Song B, Zhang YH, Wan YS and Jian P. 2002. Sampling dating and discussion of the zircon SHRIMP. Geological Review , 48 (Suppl) :26–30.
[] Sun FX, Wu GX and Yang P. 1996. Geological features and metallogenetic model from Tuanjiegou gold deposit, Heilongjiang Province, China. Jilin Geology , 15 (2) :52–60.
[] Sun FY, Wang L, Huo L and Wang KY. 2008. Fluid inclusion study on Wulaga gold deposit in Heilongjiang Province and implications for ore genesis. Geology in China , 35 (6) :1267–1273.
[] 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 Oceanic Basalts. Geological Society Special Publication :313–345.
[] Tang KD, Wang Y, He GQ and Shao JA. 1995. Continental-margin structure of Northeast China and its adjacent areas. Acta Geologica Sinica , 69 (1) :16–28.
[] Wei YF and Liu CH. 1996. The geological and geochemical model of Tuanjiegou gold deposit. Gold , 2 :3–11.
[] White NC and Hedenquist JW. 1990. Epithermal environments and styles of mineralization: Variations and their causes, and guidelines for exploration. In: Hedenquist JW, White NC and Siddeley G (eds.). Epithermal Gold Mineralization of the Circum-Pacific: Geology, Geochemistry, Origin and Exploration, I. J. Geochem. Explor., 36: 445-468
[] Wu FY, Jahn BM and Lin Q. 1997. Isotopic characteristics of the postorogenic granites in orogenic belt of northern China and their implications in crustal growth. Chinese Science Bulletin , 43 (5) :420–424.
[] Wu FY, Sun DY and Lin Q. 1999. Petrogenesis of the Phanerozoic granites and crustal growth in Northeast China. Acta Petrologica Sinica , 15 (2) :181–189.
[] Wu FY, Jahn BM, Wilde SA and Sun DY.2000. Phanerozoic crustal growth: U-Pb and Sr-Nd isotopic evidence from the granites in northeastern China. : 89 -113.
[] Wu FY, Yang JH, Lo CH Wilde SA, Sun DY and John BM.2007. The Heilongjiang Group: A Jurassic accretionary complex in the Jiamusi Massif at the western Pacific margin of northeastern China. : 156 -172.
[] Wu FY, Li XH, Yang JH and Zheng YF. 2007. Discussions on the petrogenesis of granites. Acta Petrologica Sinica , 23 (6) :1217–1238.
[] Wu SQ.1995. Geology of Tuanjiegou Porphyry Gold Deposit in Heilongjiang. Beijing: Seismological Press : 5 -44.
[] Wu YB and Zheng YF. 2004. Genesis of zircon and its constraints on interpretation of U-Pb age. Chinese Science Bulletin , 49 (15) :1554–1569. DOI:10.1007/BF03184122
[] Xie HQ, Zhang FQ, Miao LC, Chen FK and Liu DY. 2008. Zircon SHRIMP U-Pb dating of the amphibolite from “Heilongjiang Group” and the granite in Mudanjiang area, NE China, and its geological significance. Acta Petrologica Sinica , 24 (6) :1237–1250.
[] Xie LW, Zhang YB, Zhang HH, Sun JF and Wu FY. 2008. In situ simultaneous determination of trace elements, U-Pb and Lu-Hf isotopes in zircon and baddeleyite. Chinese Science Bulletin , 53 (10) :1565–1573.
[] Yin XJ, Lian YL and Pan C. 2010. Rock geochemical characteristics of Gaosongshan gold district of Xunke, Heilongjiang Province. Gold , 31 (10) :22–26.
[] Yu JB, Su RK and Liu ZM. 2005. Ore-controlled facer and ore-forming material originate of the Donggan gold deposit, Heilongjiang Province. Gold Science and Technology , 13 (6) :8–11.
[] Zeng QD, Liu JM, Zhang ZL, Jia CS, Yu CM, Ye J and Liu HT. 2009. Geology and lead-isotope study of the Baiyinnuoer Zn-Pb-Ag deposit, south segment of the Da Hinggan Mountains, Northeastern China. Resource Geology , 59 (2) :170–180. DOI:10.1111/rge.2009.59.issue-2
[] Zeng QD, Liu JM, Yu CM, Ye J and Liu HT. 2011a. Metal Deposits in the Da Hinggan Mountains, NE China: Styles, characteristics, and exploration potential. International Geology Review , 53 :846–878. DOI:10.1080/00206810903211492
[] Zeng QD, Liu JM, Zhang ZL, Chen WJ and Zhang WQ. 2011b. Geology and geochronology of the Xilamulun molybdenum metallogenic belt in eastern Inner Mongolia, China. International Journal of Earth Sciences , 100 :1791–1809. DOI:10.1007/s00531-010-0617-z
[] Zhang LG.1995. Block-Geology of Eastern Asia Lithosphere: Isotope Geochemistry and Dynamics of Upper Mantle, Basement and Granite. Beijing: Science Press : 153 -173.
[] Zindler A and Hart S. 1986. Chemical geodynamics. Ann. Rev. Earth Plant. Sci. , 14 :493–571. DOI:10.1146/annurev.ea.14.050186.002425
[] 敖贵武, 刘志明, 薛明轩, 裴志霞, 汪长生, 陈海明, 姜迎久.2004. 黑龙江东安金矿床成矿机理. 黄金地质 , 10 (1) :15–19.
[] 边红业, 陈满, 刘洪利, 赵春荣.2009. 黑龙江省逊克县高松山金矿床地质特征及成因分析. 地质与资源 , 18 (2) :91–95.
[] 郭继海, 汪长生, 石耀军.2004. 黑龙江东安金矿地质及地球化学特征. 地质与勘探 , 40 (4) :37–41.
[] 洪大卫, 王式洸, 谢锡林, 张季生.2000. 兴蒙造山带正εNd(t)值花岗岩的成因和大陆地壳生长. 地学前缘 , 7 (2) :441–456.
[] 洪大卫, 王式洸, 谢锡林, 张季生, 王涛.2003. 从中亚正εNd(t)值花岗岩看超大陆演化和大陆地壳生长的关系. 地质学报 , 77 (2) :203–207.
[] 侯增谦.2004. 斑岩Cu-Mo-Au矿床: 新认识与新进展. 地学前缘 , 11 (1) :131–144.
[] 简平, 程裕淇, 刘敦一. 2001. 变质锆石成因的岩相学研究——高级变质岩U-Pb年龄解释的基本依据. 地学前缘, 8(3): 183-191
[] 江雄新. 1986. 团结沟金矿床主要地质特征及其成矿作用. 北京: 地质出版社 : 71 -83.
[] 来又东. 2008. 黑龙江省乌拉嘎金矿矿床地质特征及矿化富集规律研究. 硕士学位论文. 长春: 吉林大学, 35-59
[] 李锦轶.1998. 中国东北及邻区若干地质构造问题的新认识. 地质论评 , 44 (4) :339–347.
[] 刘桂阁, 王恩德, 常春郊, 金宝义, 喻万强.2006. 黑龙江省逊克县高松山金矿成因探讨. 有色矿冶 , 22 (4) :1–4.
[] 刘建明, 叶杰, 徐九华, 孙景贵, 沈昆.2003. 胶东金矿床碳酸盐矿物的碳-氧和锶-钕同位素地球化学研究. 岩石学报 , 19 (4) :775–785.
[] 柳小明, 高山, 第五春荣, 袁洪林, 胡兆初.2007. 单颗粒锆石的20μm小斑束原位微区LA-ICP-MS U-Pb年龄和微量元素的同时测定. 科学通报 , 52 (4) :228–235.
[] 聂喜涛. 2010. 黑龙江省团结沟金矿床成矿地质背景与成矿模式研究. 硕士学位论文. 长春: 吉林大学, 34-50
[] 毛景文, 赫英, 丁悌平.2002. 胶东金矿形成期间地幔流体参与成矿过程的碳氧氢同位素证据. 矿床地质 , 21 (2) :121–128.
[] 毛景文, 李厚民, 王义天, 张长青, 王瑞廷.2005. 地幔流体参与胶东金矿成矿作用的氢氧碳硫同位素证据. 地质学报 , 79 (6) :839–857.
[] 曲关生. 2008. 黑龙江省岩石地层. 武汉: 中国地质大学出版社 : 1 -300.
[] 沈渭洲, 凌洪飞, 李武显, 黄小龙, 王德滋.1999. 中国东南部花岗岩类Nd-Sr同位素研究. 高校地质学报 , 5 (1) :22–25.
[] 宋彪, 张玉梅, 万渝生, 简平.2002. 锆石SHRIMP样品靶制作、年龄测定及有关现象讨论. 地质论评 , 48 (增刊) :26–30.
[] 孙凤兴, 吴国学, 杨鹏.1996. 团结沟金矿床地质模型. 吉林地质 , 15 (2) :52–60.
[] 孙丰月, 王力, 霍亮, 王可勇.2008. 黑龙江乌拉嘎大型金矿床流体包裹体特征及矿床成因研究. 中国地质 , 35 (6) :1267–1273.
[] 魏仪方, 刘春华.1996. 团结沟金矿床地质地球化学模型. 黄金 , 2 :3–11.
[] 唐克东, 王莹, 何国琦, 邵济安.1995. 中国东北及邻区大陆边缘构造. 地质学报 , 69 (1) :16–28.
[] 吴福元, 江博明, 林强.1997. 中国北方造山带造山后花岗岩的同位素特点与地壳生长意义. 科学通报 , 42 (20) :2188–2192.
[] 吴福元, 孙德有, 林强.1999. 东北地区显生宙花岗岩的成因与地壳增生. 岩石学报 , 15 (2) :181–189.
[] 吴福元, 李献华, 杨进辉, 郑永飞.2007. 花岗岩成因研究的若干问题. 岩石学报 , 23 (6) :1217–1238.
[] 吴尚全. 1995. 黑龙江省团结沟斑岩金矿地质. 北京: 地震出版社 : 5 -44.
[] 吴元保, 郑永飞.2004. 锆石成因矿物学研究及其对U-Pb年龄解释的制约. 科学通报 , 49 (16) :1589–1604.
[] 颉颃强, 张福勤, 苗来成, 陈福坤, 刘敦一.2008. 东北牡丹江地区“黑龙江群”中斜长角闪岩与花岗岩的锆石SHRIMP U-Pb定年及其地质学意义. 岩石学报 , 24 (6) :1237–1250.
[] 谢烈文, 张艳斌, 张辉煌, 孙金凤, 吴福元.2008. 锆石/斜锆石U-Pb和Lu-Hf同位素以及微量元素成分的同时原位测定. 科学通报 , 53 (2) :220–228.
[] 于建波, 苏仁奎, 刘志明.2005. 东安金矿床控矿因素及成矿物质来源浅析. 黄金科学技术 , 13 (6) :8–11.
[] 尹西君, 连永牢, 潘超.2010. 黑龙江省逊克县高松山金矿区岩石地球化学特征. 黄金 , 31 (10) :22–26.
[] 张理刚.1995. 东亚岩石圈块体地质-上地幔、基底和花岗岩同位素地球化学及其动力学 :153–173.