岩石学报  2012, Vol. 28 Issue (10): 3334-3346   PDF    
浙江安吉铅锌多金属矿床岩浆侵位序列与成矿控制
谢玉玲, 唐燕文, 李应栩, 刘保顺, 邱立明, 张欣欣, 姜妍岑     
北京科技大学土木与环境工程学院, 北京 100083
摘要: 本文在野外地质调研基础上, 通过室内岩矿相鉴定和岩石化学分析, 结合锆石U-Pb年代学研究结果, 对安吉矿区主要侵入岩单元的岩浆侵位序列、岩浆演化、构造背景及矿化控制进行了探讨。矿区岩浆活动时限为141~117Ma, 处于晚侏罗-早白垩世中国东部构造活化和构造机制转化阶段。矿区不同侵入岩单元的侵位顺序依次为黑云母二长花岗岩、正长花岗岩、花岗闪长岩、细粒花岗岩及双峰式脉岩组合, 除双峰式脉岩组合外, 其余岩石单元构成坞山关杂岩体。杂岩体不同岩石单元间具有一定的地球化学亲缘关系, 其成分变化是先存断裂的多次张开和焊结、深源物质注入和结晶分异联合作用的结果。双峰式脉岩组合是区域岩石圈减薄高峰、深部物质快速上涌的表现, 也是中国东部应力场进入引张阶段的标志。矿区主要经历了两次矿化事件。第一期与花岗闪长岩(137Ma) 侵位有关, 形成以石榴子石矽卡岩为特征的矽卡岩型铁铜矿化, 其主要分布于矿区西部;第二期为与细粒花岗岩(134~135Ma) 侵位有关的铅锌(银)、钼矿化, 形成了以绿帘石矽卡岩为主的矽卡岩型铅锌(银) 矿化, 同时在矿区北部的石灰岩地层中形成了碳酸盐交代型铅锌矿化、在黑云母二长花岗岩中形成脉状、网脉状铅锌矿化及网脉状-浸染状辉钼矿化。主要侵入岩单元的岩石化学表明, 与成矿有关的花岗闪长岩显示了I型花岗岩的地球化学特征, 而细粒花岗岩的岩石化学特征界于I型与S型之间, 但均表现有深源物质加入, 细粒花岗岩相对于花岗闪长岩显示相对氧化条件和较强的分异特征。深源物质加入和进一步的结晶分异是造成熔体中成矿金属元素富集的重要机制。矿区内与成矿有关的花岗闪长岩和细粒花岗岩主要发育于矿区东西两侧, 因此矿区西侧为铁铜勘查的有利地区, 而东部的铅锌(银) 钼矿化值得进一步工作, 但矿区后期NE和近EW向断裂构造对矿体的保存条件影响较大, 因此在勘查工作中应加强矿区构造地质方面的研究工作。
关键词: 岩石学     岩石化学     成矿控制     岩浆侵位序列     铅锌多金属矿区     浙江安吉    
Magmatic intrusive series and their implication for the ore prospecting in Anji exploration area, Zhejiang Province
XIE YuLing, TANG YanWen, LI YingXu, LIU BaoShun, QIU LiMing, ZHANG XinXin, JIANG YanCen     
Civil and Environmental Engineering Institute, University of Science and Technology Beijing, Beijing 100083, China
Abstract: Based on field observation, petrographic and lithogeochemical results, together with zircon U-Pb chronology results, this paper discusses the intrusive sequence in Anji deposit and their evolution link, tectonic setting, as well as the mineralization implication of varied intrusive units. The magmatic intrusions in this area occur during Late Jurassic-Early Cretaceous (141~117Ma), and they are the products of tectonic reactivation and tectonic conversion in eastern China during Late Mesozoic. The intrusive sequences are from early biotite monzonitic granite followed by syenogranite, granodiorite, fine-crystalline granite and bimodal dikes. Except biomodal dikes, the other four units compose the Wushanguan complex. Geochemical affinity and compositions difference of the units in the complex may result from mantle material injection and crystallization differentiation of magma during faults activation at Late Jurassic-Early Cretaceous. The bimodal dikes marked lithospheric thinning peak and mantle material upwelling, as well as the stress field conversion into extension in eastern China. Two metallogenic events are recognized in Anji area including early Fe-Cu metallogenic event which genetically related to granodiorite and late Pb-Zn (Ag)-Mo metallogenic event which is associated with fine-crystalline granite. Geochemically, the granodiorite shows Ⅰ-type granite character, whereas the fine-grained granite shows the combined character of S-type and Ⅰ-type granite, and both show more mantle materials input than the other units in Wushanguan complex. Mantle materials involving and strong differentiation may be the reason for metal elements (e.g. Fe, Cu, Mo, Pb, Zn and Ag) enrichment and related mineralization. The causative intrusive units occurs in the west and east part of the deposit respectively, so, the west of the deposit is the favorable Fe-Cu exploration area, and the east of the deposit is the favorable Pb-Zn (Ag) as well as Mo exploration area. Well-developed post-mineralization NE and EW faults have seriously influence to the preservation of ores, so the structure geology studies are very important to further exploration in Anji area.
Key words: Petrography     Petrogeochemistry     Ore-forming control     Intrusive sequences     Pb-Zn polymetallotic deposit     Anji, Zhejiang Province    

安吉港口铅锌多金属矿位于钦杭成矿带东段(图 1a)(杨明桂和梅勇文, 1997杨明桂等, 2009) 北缘。钦杭成矿带是我国新近确立的又一重要成矿带, 其大体呈NE向沿钦杭(钦州湾-杭州湾) 古板块结合带及其两侧展布。钦杭结合带位于扬子和华夏两板块之间, 南起广西钦州湾, 经梧州-郴州-萍乡-德兴, 北至浙江杭州。杨明桂和梅勇文(1997)将绍兴-北海断裂带作为其东南界线, 而将凭祥-歙县断裂带作为其西北界线。钦杭结合带在长沙-萍乡附近走向明显变化, 以此为界可将该成矿带分为东、西两段, 安吉矿区即位于该成矿带的东段的东北缘。德兴大型矿集区是该成矿带中较早发现、经济意义巨大的已知矿集区, 除此之外, 在该成矿带的东段和西段近年来找矿工作取得重大进展, 在江西、皖南等地相继发现数十个大中型铜、金、钨、钼和铅锌矿床。在浙西发现的淳安银山铅锌多金属矿(何国锦等, 2011) 以及安吉铅锌多金属矿代表了钦杭成矿带东缘的找矿突破, 填补了钦杭成矿带浙西地区找矿工作的空白。

图 1 钦杭成矿带东段简图(a, 据杨明桂和梅勇文, 1997修改) 及安吉矿区地质图(b, 据南京地质矿产研究所, 2010, 修绘) 1-第四系;2-印渚埠组; 3-西阳山组;4-华严寺组;5-杨柳岗组;6-大陈岭组;7-荷塘组;8-正长花岗岩9-花岗闪长岩;10-黑云母二长花岗岩;11-细粒花岗岩;12-花岗斑岩;13-正长花岗斑岩(或霏细岩) 脉;14-辉绿岩脉;15-石英二长斑岩脉;16-石英正长斑岩;17-硅化破碎带;18-断层;19-矿化带;20-矿体;21-沿硅化破碎带发育的矿化;22-村庄 Fig. 1 The sketch map of east part of Qin-Hang metallogenic belt (a, after Yang and Mei, 1997) and geological map of the Anji deposit (b) 1-Quaternary; 2-Ordovician Yinzhifu Formation; 3-Cambrian Xiyangshan Formation; 4-Cambrian Huayansi Formation; 5-Cambrian Yangliugang Formation; 6-Cambrian Dachenling Formation; 7-Cambrian Hetang Formation; 8-syenogranite; 9-granodiorite; 10-biotite monzonitic granite; 11-fine crystalline granite; 12-granite porphyry; 13-syenogranite porphyry (felsite); 14-diabase; 15-beschtauite; 16-quartz syenite porphyry; 17-silicatized fracture zone; 18-fault; 19-mineralized zone; 20-ore-body; 21-mineralized zone and fracture zone; 22-village

①南京地质矿产研究所.2010.浙江省安吉县港口地区金多金属矿地质报告.1-100(内部资料)

安吉港口矿区位于浙江省安吉县港口乡, 矿区内矿化类型多样, 包括矽卡岩型铁铜矿化、脉状铁铜矿化、矽卡岩型铅锌矿化、碳酸盐地层中的交代型铅锌矿化、热液脉型铅锌银矿化和萤石矿化及产于岩体中呈细网脉状和浸染状的辉钼矿化。矿区岩浆活动强烈, 侵位序列复杂, 但对矿区内岩浆侵位序列及与成矿的关系研究尚相当薄弱。开展安吉多金属矿床详细的地质和地球化学研究, 将对认识该矿床的成因, 建立成矿模型具有重要意义, 同时也对研究钦杭成矿带东缘的构造演化、岩浆活动及相关的成矿过程具有重要意义。本文在详细的野外地质调研、矿床地质、主要侵入岩单元的岩石学、岩石化学分析基础上, 结合锆石U-Pb定年结果, 对矿区内主要侵入岩的单元的侵位序列、岩浆演化、构造背景及成矿控制进行了探讨。结果表明, 矿区岩浆活动主要以燕山期中酸性侵入活动为主, 另外还发育有基性、酸性和碱性脉岩, 目前控制的岩浆活动时限从141Ma延续至117Ma。岩浆侵位序列从早到晚依次为:黑云母二长花岗岩、正长花岗岩、花岗闪长(斑) 岩、细粒花岗岩及双峰式脉岩组合。区内矽卡岩型铁(铜、锌) 矿化主要与花岗闪长岩有成因联系, 而铅锌(银)、钼矿化主要与区内细粒花岗岩有关。花岗闪长岩与古生界灰岩地层接触部位控制了矿区矽卡岩型铁(铜、锌) 矿化, 细粒花岗岩与灰岩接触部位控制了矽卡岩型铅锌矿化。矿区经历了两次主要的成矿事件, 与成矿有关的岩体侵位年龄约为134Ma和137Ma, 总体介于晚侏罗世和早白垩世, 与中国东部大规模成岩成矿事件时间(毛景文和王志良, 2000) 和长江中下游成矿事件(周涛发等, 2008) 时间吻合。笔者认为, 晚中生代中国东部大陆岩石圈活化和构造机制从压扭向张扭转化使先存断裂的多次张开、焊结, 以及与之相伴的深源物质注入和岩浆的结晶分异是造成杂岩体内不同侵入岩单元成分变异的主要原因。

1 矿床地质特征 1.1 矿区地质概况

位于浙江西北部的安吉矿区夹持于天目山、莫干山两个火山盆地之间。矿区大地构造上属扬子地块皖南-浙西地向斜(杨松生, 1991)。矿区内除广泛分布的第四系外, 主要出露寒武系和奥陶系(图 1b), 从老至新分别为下寒武统荷塘组(€1h) 炭质、泥质、硅质页岩;下寒武统大陈岭组(€1d) 灰岩, 中寒武统杨柳岗组(€2y) 泥质白云质灰岩、灰岩、硅质泥岩、炭质页岩, 上寒武统华严寺组(€3h) 灰岩、泥质白云质灰岩, 上寒武统西阳山组(€3x) 白云质、泥质灰岩和灰岩, 下奥陶统印渚埠组(O1y) 粉砂质泥岩和泥岩, 矿区内地层间均为整合接触(南京地质矿产研究所, 2010)。矿区地处北东向于潜-三桥埠复向斜南翼, 地层整体NE倾, 构造以断裂构造为主, 主要发育NEE向断裂、NW向断裂和NE向断裂, 以及NNW、NWW、NEE等方向的次级断裂。区内岩浆活动强烈, 以燕山期侵入岩最为发育, 主要表现为坞山关杂岩体和一系列酸性、基性、碱性岩脉。

1.2 矿区的蚀变、矿化类型及空间分布特征

矿区蚀变作用广泛发育, 包括岩体内的绿泥石化、粘土矿化、黄铁矿化、硅化、绢云母化、钾化、绿帘石化以及灰岩地层的大理岩化、灰岩与岩体接触带附近的矽卡岩化、泥岩的角岩化等, 其中绿泥石化、黄铁矿化、硅化、绢云母化、矽卡岩化是矿区最为常见的蚀变类型, 以绿泥石化、黄铁矿化分布范围最大。

矿区目前勘查发现的金属矿化类型多样, 包括矽卡岩型铁铜矿化、矽卡岩型铅锌矿化、脉状-网脉状铅锌矿化、碳酸盐地层交代型铅锌矿化及岩体中的脉状-网脉状-浸染状辉钼矿化。矽卡岩型铁铜(锌) 矿化发育于矿区西部观音堂-五庙桥一带, 主要产于花岗闪长岩与寒武系灰岩地层的接触部位(图 1b), 以铁、锌为主, 局部含铜较高, 矽卡岩类型以石榴子石矽卡岩为主, 另发育少量透辉石矽卡岩和绿帘石矽卡岩, 主要矿石矿物为磁铁矿、铁闪锌矿、黄铜矿等, 脉石矿物主要为石榴子石、透辉石、石英、绿帘石等。矽卡岩型铅锌矿化主要见于矿区东部, 主要产于细粒花岗岩与寒武系灰岩地层的接触部位, 以绿帘石矽卡岩为主, 主要金属矿物组合为方铅矿、闪锌矿、黄铁矿、黄铜矿等, 脉石矿物主要为石英、方解石、绿帘石、阳起石等。与铅锌矿化有关的矽卡岩类型明显不同于与铁铜矿化有关的石榴子石矽卡岩。碳酸岩地层交代型铅锌矿化主要发育于矿区东部的北侧寒武系地层中, 呈囊团状或不规则状, 可见灰岩地层的透闪石化、大理岩化。脉状-网脉状-浸染状辉钼矿化主要产于矿区东部桐坞里一带和俞家坞村口, 以辉钼矿-石英脉(脉宽 < 1~2cm)、石英-绿泥石-辉钼矿细脉或网脉、辉钼矿细网脉及浸染状辉钼矿化为主, 产于黑云母二长花岗岩中。

2 矿区主要侵入岩单元的产状、分布及岩石学特征

通过野外地质调研和室内岩矿相鉴定, 并结合锆石U-Pb定年结果(谢玉玲等, 2012唐燕文等, 2012a, b;谢玉玲等, 未发表资料), 矿区出露的侵入岩单元包括(按岩体侵位先后顺序):

黑云母二长花岗岩  呈岩株产于矿区中部, 在东部和西部也有零星出露。岩石呈浅肉红色, 等粒-似斑状结构, 块状构造, 斑晶由自形-半自形钾长石组成, 基体主要由半自形-他形斜长石、石英、钾长石以及黑云母组成。副矿物包括锆石、磷灰石、磁铁矿。矿区范围内岩体普通发育蚀变, 主要表现为斜长石的绢云母化、泥化和黑云母的绿泥石-黄铁矿化及脉状、网脉状硅化、绿泥石化、黄铁矿化等。锆石U-Pb定年结果表明其侵位时间为141.0±1.4Ma。

正长花岗岩  呈岩株出露于矿区中南部和东南部。岩石呈肉红色, 等粒结构, 中等粒度, 块状构造。矿物组成为钾长石、斜长石、石英;暗色矿物主要为黑云母。副矿物有锆石、磁铁矿等。镜下可见斜长石的双晶和斜长石的绢云母化蚀变。锆石U-Pb定年结果表明其侵位时间为138.1±1.0Ma。

花岗闪长(斑) 岩  浅灰到灰白色, 不等粒-斑状结构(图 2a), 块状构造, 主要出露于矿区西部, 出露面积约1.8km2。组成矿物主要为斜长石, 钾长石, 石英;暗色矿物为角闪石和黑云母, 副矿物主要为锆石、磁铁矿等。镜下可见斜长石的聚片双晶、卡钠复合双晶和环带结构及角闪石的双晶等。矿区西北部主要为花岗闪长斑岩, 具有多斑的特征, 斑晶含量在60%~80%, 斑晶主要为斜长石、钾长石和角闪石。局部蚀变强烈, 原岩结构不清, 主要可见斜长石的绢云母化、角闪石的碳酸盐化、绿泥石化和黑云母的绢云母、绿泥石化、黄铁矿化等。锆石U-Pb定年结果表明其侵位时间为137.0±1.8Ma。

图 2 不同岩石单元岩石薄片的显微镜下照片 (a)-花岗闪长斑岩(正交偏光);(b)-石英二长斑岩(正交偏光);(c)-石英二长斑岩中角闪石的碳酸盐化(单偏光);(d)-细粒花岗岩(单偏光);(e)-细粒花岗岩(正交偏光);(f)-强绢云母化细粒花岗岩(正交偏光);(g)-正长花岗斑岩中浑圆状石英和长石斑晶及钾长石的环斑结构(正交偏光);(h)-霏细岩的球粒结构(正交偏光);(i)-辉绿岩的辉绿结构(正交偏光).Kf-钾长石;Q-石英;Hb-角闪石;Pl-斜长石;Bi-黑云母;Ser-绢云母;Px-辉石;Me-基质;R-硅酸盐球粒;Py-黄铁矿;Cal-方解石 Fig. 2 Micrographies polarization photos of intrusions in Anji mining area (a)-porphyritic granodiorite (CPL); (b)-beschtauite (CPL); (c)-carbonized hornblende in beschtauite (PPL); (d)-fine crystalline granite (PPL); (e)-fine crystalline granite (CPL); (f)-sericitized fine crystalline granite (CPL); (g)-round quartz and K-feldspar phenocryst in syenogranitic porphyry (CPL); (h)-eurite with centric texture (CPL); (i)-diabase with diabasic texture (CPL). Kf-K-feldspar; Q-quartz; Hb-hornblende; Pl-plagioclase; Bi-biotite; Ser-sericite; Px-pyroxene; Me-metrix; R-silicate chondri; Py-pyrite; Cal-calcite

石英二长斑岩  呈岩脉产出于寒武系地层和黑云母二长花岗岩中, 见于矿区中北部, 全长大于1km。岩石呈灰白色, 斑状结构(图 2b, c), 块状构造。斑晶包括钾长石、斜长石、石英、角闪石, 偶见黑云母斑晶, 岩石普遍发育强度不一的蚀变, 可见晚期张性石英脉穿插其中。靠近石英脉处岩石蚀变较强, 远离石英脉时蚀变较弱。常见长石的星点状绢云母化、角闪石的碳酸盐化及绿泥石和黄铁矿交代角闪石现象;副矿物主要为锆石、磁铁矿等。锆石U-Pb定年结果表明其侵位时间为137.3±1.6Ma。

细粒花岗岩  呈小岩株出露于矿区东侧, 侵位于黑云母二长花岗岩中, 出露面积约0.4km2。岩石呈肉红色, 不等粒-斑状结构(图 2d, e), 结晶粒度在 < 0.1~0.5mm, 中要集中在细粒范围, 个别达中粒。主要矿物组成为石英、钾长石, 少量斜长石和暗色矿物。岩石蚀变常见, 且局部蚀变较强, 主要为石英绢云母化(图 2f), 获得的两组锆石U-Pb定年结果分别为133.9±1.3Ma和134.5±1.6Ma (谢玉玲等, 2012)。

双峰式脉岩组合  矿区常见酸性和基性脉岩相伴产出, 成对出现, 矿区范围内可见正长花岗斑岩-辉绿岩组合和霏细岩(霏细斑岩)-辉绿岩组合, 呈岩墙或岩脉产出。

酸性脉岩在矿区东南部主要表现为正长花岗斑岩, 岩石呈暗红色至肉红色, 斑状结构、块状构造, 斑晶主要为钾长石、石英, 另可见少量黑云母斑晶。常见钾长石的环斑结构(图 2g) 及钾长石和石英的聚斑, 基质隐晶, 无论是长石斑晶还是石英斑晶均呈浑圆状, 溶蚀现象明显。锆石U-Pb定年结果表明其侵位时间为122.1±2.0Ma (谢玉玲等, 未发表资料);在矿区西北部主要表现为霏细岩和霏细斑岩, 岩石具霏细结构、球粒结构(图 2h) 或斑状结构, 岩石呈白色, 磁状断口, 具有斑状结构时, 斑晶为钾长石, 霏细岩锆石U-Pb定年结果表明其侵位时间为125.9±1.5Ma。

基性脉岩在空间上与正长花岗斑岩(或霏细岩、霏细斑岩) 相伴产出, 有时可见两者共同构成同一脉体的两侧或辉绿岩分布在正长花岗斑岩两侧。岩石呈灰黑-灰绿色, 具典型的辉绿结构(图 2i)、块状构造, 矿物组成为斜长石、辉石、角闪石、黑云母和钛铁矿。副矿物主要为锆石等。在矿区南部有时可相变为辉石闪长玢岩。锆石U-Pb定年结果表明其侵位时间约为117.0±2.0Ma (谢玉玲等, 未发表资料)。

除上述侵入岩单元外, 矿区还发现有花岗斑岩和石英正长斑岩, 其与其它岩石单元的关系尚不清楚。

3 主要侵入岩单元的主量和稀土、微量元素组成特征 3.1 样品及测试方法

岩石主量、微量和稀土元素在核工业北京地质研究院分析测试中心完成。主量元素含量分析使用Philips PW2404型X荧光光谱仪(XRF) 完成, 分析精度优于1%;微量元素和稀土元素分析使用Finnigan MAT ElementⅠ型电感耦合等离子体质谱仪(ICP-MS) 完成, 分析精度多优于3%。

本文在唐燕文等(2012a, b) 和谢玉玲等(2012)基础上, 补充了双峰式脉岩组合的岩石化学分析结果(表 1)。

表 1 双峰式脉岩组合主量(wt%)、微量和稀土元素(×10-6) 元素分析结果 Table 1 Major (wt%), trace and rare earth (×10-6) element contents of bimodal dikes
3.2 主量元素分析结果

黑云母二长花岗岩  SiO2含量变化于69.42%~70.66%之间, Al2O3含量为13.45%~13.81%, K2O+Na2O为7.93%~8.33%, K2O/Na2O主要变化于1.44~1.62。属钾玄岩系列;铝饱和指数A/CNK介于0.93~0.97之间, A/NK介于1.25~1.31, 属于碱性-准铝质岩石。

正长花岗岩  SiO2含量变化于73.76%~75.27%之间, K2O+Na2O为8.06%~8.84%, K2O/Na2O主要变化于1.62~1.8, 属钾玄岩系列;Al2O3含量为12.04%~12.84%, 铝饱和指数A/CNK介于0.97~1.06之间, A/NK介于1.12~1.19, 其投点位于准铝质与过铝质范围交界附近。

花岗闪长岩  SiO2含量变化于66.73%~67.79%之间, K2O+Na2O为6.95%~7.25%, K2O/Na2O主要变化于1.04~1.21, 属高钾钙碱性-钾玄岩系列;Al2O3含量为13.73%~13.93%, 铝饱和指数A/CNK介于0.87~0.89之间, A/NK介于1.43~1.48, 属于准铝质岩石。

石英二长斑岩   SiO2含量变化于65.91%~66.76%之间, K2O/Na2O主要变化于1.58~1.82, 平均1.71;属钾玄岩系列;Al2O3含量为13.28%~13.72%, 平均13.5%, 铝饱和指数A/CNK介于0.86~0.91之间, A/NK介于1.25~1.32, 属于准铝质岩石。

细粒花岗岩  SiO2含量74.48%~75.66%, Al2O3含量12.52%~13.08%, Na2O+K2O含量8.14%~8.46%, 铝饱合指数A/CNK介于1.07~1.13, A/NK介于1.14~1.24, 属高钾钙碱性、过铝质岩石。

霏细岩  SiO2含量变化于76.07%~76.86%之间, K2O+Na2O为9.18%~9.46%, K2O/Na2O主要变化于1.47~1.80, 为钾玄岩系列;Al2O3含量为11.22%~11.42%, 铝饱和指数A/CNK介于0.82~0.90之间, A/NK介于0.92~0.96, 属于过碱质岩石。

正长花岗斑岩  SiO2含量变化于72.39%~73.07%之间, K2O+Na2O为8.64%~8.86%, 平均8.75%;K2O/Na2O主要变化于1.92~2.12, 属钾玄岩系列;Al2O3含量为12.28%~12.93%, 平均12.58%;铝饱和指数A/CNK介于0.92~0.95之间, A/NK介于1.12~1.17, 属于准铝质岩石。

辉绿岩  SiO2含量变化于50.36%~51.41%之间, K2O+Na2O为5.06%~5.36%;K2O/Na2O主要变化于0.28~0.45, 为钙碱性-高钾钙碱性系列;Al2O3含量为15.76%~15.99%, 平均15.85%, 铝饱和指数A/CNK介于1.11~1.16之间, A/NK介于2.02~2.11, 属于过铝质岩石。

岩石化学分析表明, 花岗闪长岩属Ⅰ型花岗岩, 而细粒花岗岩的地球化学特征界于Ⅰ型和S型之间。SiO2-K2O图解(图 3a) 中花岗闪长岩落在钙碱性区域, 而石英二长斑岩和细粒花岗岩落在碱性区域。在A/NK-A/CNK图解(图 3b) 中均落入过铝质-准铝质区域。

图 3 矿区主要侵入岩单元的K2O-SiO2 (a) 和A/CNK-A/NK (b) 图解(采用路远发, 2004 GeoKit软件绘制) Fig. 3 The SiO2-K2O and A/CNK-A/NK diagrams (Using GeoKit software by Lu, 2004)
3.3 稀土和微量元素分析结果

黑云母二长花岗岩  稀土总量为164.2×10-6~216.5×10-6, (La/Yb)N=21.21~23.93, 显示轻稀土富集型稀土配分曲线(图 4a);δEu=0.50~0.67。微量元素蛛网图(图 4b) 显示, Rb、K、U、Th, Pb明显富集, 相对亏损Ba、Nb、Ta、P。

图 4 矿区主要侵入岩单元的稀土元素球粒陨石标准化配分曲线(a, 标准值据Boynton, 1984) 和原始地幔标准化微量元素蛛网图(b, 标准值据McDonough and Sun, 1995) Fig. 4 Chondrite-normalized REE distribution patterns (a, normalization values after Boynton, 1984) and primitive mantle normalized trace elements spider diagram (b, normalization values after McDonough and Sun, 1995) of intrusions in Anji deposit

正长花岗岩  稀土总量为124.6×10-6~163.8×10-6, (La/Yb)N=7.47~9.90, 轻稀土相对富集, 重稀土相对轻稀土亏损(图 4a);δEu=0.38~0.50, Eu负异常明显。微量元素蛛网图(图 4b) 总体显示为右倾型, 即相对富集不相容元素, Rb、K、U、Th、Pb明显富集, 但Ba、P、Ti相对亏损。

花岗闪长(斑) 岩  稀土总量为145.5×10-6~151.5×10-6, LREE/HREE为2.96~3.97, 平均3.65;(La/Yb)N=9.69~13.45, 平均11.59;δEu=0.80~0.86, 平均0.81。稀土配分曲线显示为右倾的轻稀土富集型(图 4a), 有微弱的负Eu异常。微量元素蛛网图(图 4b) 中Rb、K、U、Th、K、Pb和稀土明显富集, 相对亏损高场强元素Nb、Ta、Zr和P。

石英二长斑岩  稀土总量为128.2×10-6~148.9×10-6, 平均137.1×10-6, ∑LREE/∑HREE=3.19~3.46, 平均3.38, (La/Yb)N=9.99~10.82, 平均10.32, 显示轻稀土富集型分配曲线(图 4a);δEu=0.74~0.80, 平均0.76, Eu负异常较弱。微量元素蛛网图(图 4b) 显示, U、Th、K、Pb明显富集。

细粒花岗岩  稀土总量在76.8×10-6~123.7×10-6, (La/Yb)N=13.74~21.41。显示轻稀土富集型分配曲线(图 4a)。重稀土相对轻稀土亏损, 轻稀土内部分异明显, (La/Sm)N=7.81~12.94, 而重稀土内部分异不明显, (Gd/Lu)N=0.77~1.06。δEu=0.65~0.83, 铕异常不明显。微量元素蛛网图(图 4b) 显示总体右倾型, ) Rb、K、U、Th和Pb明显富集, 相对亏损Nb、Ti、P。

正长花岗斑岩  稀土总量为398.2×10-6~412.9×10-6, 平均403.2×10-6, ∑LREE/∑HREE=2.96~3.04, 平均3.15, (La/Yb)N=8.27~9.05, 平均8.50, 重稀土相对轻稀土亏损(图 4a);δEu=0.20~0.24, 平均0.22, Eu负异常明显。微量元素蛛网图(图 4b) 显示, 总体右倾, 即富集不相容元素, 且Rb、K、Th、U、Pb明显正异常, 而Ti、P呈明显负异常。

霏细岩  稀土总量为292.7×10-6~315.1×10-6, 平均302.8×10-6, ∑LREE/∑HREE=1.76~1.96, 平均1.88, (La/Yb)N=4.04~4.39, 平均4.23, 轻重稀土分异较弱(图 4a)。δEu=0.04, Eu负异常明显(图 4a)。微量元素蛛网图(图 4b) 显示明显亏损Ba、P和Ti。

辉绿岩  稀土总量∑REE=154.5×10-6~163.5×10-6, 平均值156.0×10-6, ∑LREE/∑HREE=3.19~3.26, 平均3.23, (La/Yb)N=9.09~9.80, 平均9.38, 轻稀土相对富集(图 4a);重稀土相对轻稀土亏损;δEu=0.97~1.07, 平均1.01, 无Eu异常明显。微量元素蛛网(图 4b) 图较为平缓, 但U明显正异常, 而Nb和Ta有明显的负异常。

4 讨论 4.1 岩浆侵位序列、岩浆演化及构造背景

野外观察表明, 花岗闪长岩、石英二长斑岩、细粒花岗岩及双峰式脉岩组合均可见侵位于黑云母二长花岗岩中, 侵位关系清楚, 但正长花岗岩与黑云母二长花岗岩呈断层接触, 其间可见双峰式脉岩侵位, 另可见双峰式脉岩侵入花岗闪长岩和细粒花岗岩中。细粒花岗岩与花岗闪长岩未见侵位关系。由野外不同岩石组合、岩石学和岩石化学特征看, 矿区花岗闪长(斑) 岩与石英二长斑岩、正长花岗斑岩和霏细岩分别为同期岩浆活动不同侵位深度的表现。矿区的NS和EW向断裂可能是造成其剥蚀深度差异的主要原因。根据野外岩体侵位关系, 结合不同侵入岩单元的锆石U-Pb定年结果(据谢玉玲等, 2012唐燕文等, 2012a;谢玉玲等未发表资料), 矿区内岩浆侵位序列依次为黑云母二长花岗岩、正长花岗岩、花岗闪长(斑) 岩和石英二长斑岩、细粒花岗岩、双峰式脉岩组合。黑云母二长花岗岩、正长花岗岩、花岗闪长(斑) 岩和细粒花岗岩共同构成坞山关杂岩体。

岩石化学结果表明, 除双峰式脉岩组合外, 其余均属准铝质-过铝质岩石系列、钙碱性-钾玄岩, 侵位年龄相近, 且在SiO2-Al2O3、SiO2-Fe2O3+MgO+TiO2、Rb/Sr-Rb/Ba图解(图 5a-c;) 中均显示很好线性演化趋势, 表明其为同源演化的产物, 而在Sr/Y-Ta/Nb图解(图 5d) 中明显分为三个区间, 显示其成分变异特征。双峰式脉岩组合中的正长花岗斑岩(或霏细岩) 显示了与杂岩体不同的岩石化学特性。杂岩体总体表现为Rb、U、Th、K等大离子亲石元素明显富集, 相对亏损高场强元素Nb、Ti、P, 表明其源区可能有角闪石、金红石、斜长石的难容残留相;从黑云母二长花岗岩到正长花岗岩, 其中SiO2含量升高, 而Al2O3、CaO、TiO2、Fe2O3和MgO降低, 并伴有Sr/Y、Nb/Ta、La/Yb明显降低和Eu负异常的加强, 但Sr/Ba比值变化不大, 表明可能有金红石、长石、磁铁矿等的分离结晶。正长花岗岩到花岗闪长岩(石英二长斑岩), SiO2含量降低, 而Fe2O3+MgO+TiO2随之呈线性升高, 并伴有Sr、Ba、Sr/Y、Ta/Nb、La/Yb明显升高, 但Sr/Ba比值变化不大, 可能代表了一次深源物质的贯入过程。由于深源物质的加入也造成Cr、V、Ni、Cu、Zn等指示深源的成矿金属元素含量明显增高。花岗闪长岩到细粒花岗岩SiO2、Al2O3、Fe2O3+MgO+TiO2含量均降低, 且Rb/Sr、Fe2O3/FeO、Zr/TiO2和La/Yb升高, Sr、Ba含量降低, 但Sr/Ba比值变化不大, 反映了角闪石和斜长石的分离结晶作用。双峰式脉岩组合中, 正长花岗斑岩(霏细岩) 相对于辉绿岩富集大离子亲石元素(LILE) 高场强元素(HFSE), 如Zr、Hf、Nb、Ta、Cr和REE (特别是LREE), 同时Rb、Li、Be等亲流体元素和Pb富集, 而辉绿岩中相对富集Cu、Zn、Co、Ni、V等, 特别是Zn可达地壳丰度的50倍以上。

图 5 矿区花岗质岩石单元的主微量元素图解 (a)-SiO2-Al2O3图解;(b)-SiO2-Fe2O3+MgO+TiO2图解;(c)-Rb/Sr-Rb/Ba图解;(d)-Sr/Y-Ta/Nb图解 Fig. 5 Major and trace elements correlation diagrams of granitic rocks in Anji deposit

双峰式脉岩组合中的正长花岗斑岩和霏细岩, 其岩石化学、岩石定年结果均表明其与杂岩体具有不同的岩浆源区性质和不同的构造背景, 具有高硅、低铁、低Sr/Y、Ta/Nb和高Rb/Sr、Rb/Ba的特征。正长花岗斑岩和霏细岩与辉绿岩具有互补的稀土和微量元素组成式样, 霏细岩相对于辉绿岩具有强的Eu负异常。从野外侵位关系和同位素定年结果看, 正长花岗斑岩(霏细岩) 先于辉绿岩侵位, 因此正长花岗斑岩不可能是基性岩浆进一步结晶分异的产物, 而可能代表了硅酸岩熔体的不混溶过程, 互补的稀土和微元素元素组成是元素在不混溶两相之间分配造成的。

对中国东部中生代以来的构造演化前人进行了大量的工作(任纪舜等, 1984许志琴等, 1982任建业等, 1998李思田等, 1987环文林等, 1982), 结果表明, 印支运动使亚洲东部三个不同陆块(扬子、中朝、西伯利亚) 进一步叠接之后(任纪舜, 1984), 到晚中生代(J3-K1) 中国大陆东部进入又一个构造活化期, 表现为以大规模的火山喷发、花岗岩浆侵位、盆岭式断陷盆地系和变质核杂岩等为特征的大陆裂陷作用(任建业等, 1998)。中国东部构造形变和构造应力场自中生代以来经历了三个重要的演化阶段, 各阶段的构造动力与太平洋内各时期板块运动的方向和边界条件有着密切关系(环文林等, 1982)。矿区岩浆活动时限为141~117Ma, 代表了中国东部晚中生代构造活化和构造体制转换的产物。中国东部岩石圈减薄发生在晚中生代, 且在120~130Ma的早白垩世达到高峰(吴福元等, 2003), 以其间发育的双峰式岩浆作用为代表(许文良等, 2004)。

在构造环境判别R1-R2图解中(图 6a) 中黑云母二长花岗岩落入同碰撞-晚造山区域, 而正长花岗岩落入造山晚期-非造山, 花岗闪长岩落入碰撞前-碰撞后隆起, 细粒花岗岩落入晚造山-非造山边界区域, 而晚期侵位的正长花岗斑岩落入同碰撞-晚造山边界区域、霏细岩落在晚造山-非造山区域且以落入非造山区域为主;在Yb+Ta-Rb构造环境图解(图 6b), 黑云母二长花岗岩落入同碰撞区域, 石英二长斑岩落入火山弧花岗岩和同碰撞花岗岩边界区域, 花岗闪长岩则落入火山弧花岗岩区域, 细粒花岗岩和正长花岗岩落入同碰撞区域, 正长花岗斑岩、霏细岩落在板内环境区域, 结合岩石定年数据, 其总体反映从碰撞造山环境向板内(非造山) 环境的转化。晚中生代构造活化导致一系列北北东向断裂的再活动, 在燕山主期(J2-J3), 该断裂以左旋压扭作用为特征, 到J3末期, 中国东部构造应力场实现了由左旋压扭到右旋张扭的巨大反转(许志琴等, 1982)。构造应力场从左旋压扭到右旋压扭应力场向右旋张扭应力场的转化, 致使早先压扭条件下紧闭的深断裂引张开裂, 导致围压降低, 高热流背景下积聚的能量释放, 发生大规模的热事件和岩浆侵位, 同时也触发了规模巨大的裂陷作用(李思田等, 1987)。在矿区这一构造机制转化诱发了地幔软流圈上涌和区域的高热流值, 并引起下地壳的部分熔融。岩浆沿先期NE向断裂上侵形成区内黑云母二长花岗岩侵位并使断裂焊结, 深部岩浆房被封存, 在近等压条件下演化, 并发生了金红石、长石、磁铁矿等矿物的结晶分离, 并使岩浆中Fe、Mg、Ti等含量降低, 而SiO2含量增高, 断裂再次张开形成正长花岗岩侵位和断裂的再次焊结;其后由于断裂继续活动, 并向深部扩展, 幔源物质注入, 岩浆中SiO2含量相对降低, 而显示深源的元素含量增加, 花岗闪长岩和石英二长斑岩中的暗色包体是深源物质注入的证据之一。其后断裂再次被焊结, 岩浆再次处于等压降温过程, 此时由于岩浆房上移, 压力相对较低, 因此在等压降温过程中造成大量长石的结晶分离, 其后断裂再次张开形成细粒花岗岩。其后区域构造活动减弱。在120~130Ma, 由于岩石圈的持续减薄、地幔物质快速上涌, 并沿先期形成的压性构造上侵形成了双峰式脉岩组合。因此, 双峰式脉岩组合代表了区域应力从张扭进入到张性的标志, 自此, 造山作用结束, 矿区进入板内岩浆活动时期。从双峰式脉岩的形态和产状看, 其产状和厚度较为稳定, 显示原压性构造再次张开的特点。

图 6 矿区主要侵入岩单元构造环境图解 (a)-R1-R2图解, 底图据Batchelor and Bowden, 1985;(b)-Yb+Ta-Rb图解, 底图据Pearce et al., 1984.①-幔斜长花岗岩;②-破坏性活动板块边缘(板块碰撞前) 花岗岩;③-板块碰撞后隆起期花岗岩;④-晚造山期花岗岩;⑤-非造山A型花岗岩;⑥-同碰撞(S型) 花岗岩;⑦-造山期后A型花岗岩;syn-COLG-同碰撞花岗岩;VAG-火山弧花岗岩;WPG-板内花岗岩;ORG-洋岛花岗岩 Fig. 6 Discriminant diagrams of the tectonic setting of R1-R2 (a, after Batchelor and Bowden, 1985) and Rb-Yb+Ta (b, after Pearce et al., 1984) for intrusions in Anji deposit
4.2 岩浆活动与成矿控制

花岗岩与成矿关系密切, 特别是斑岩型和热液型矿床, 前人对花岗岩与成矿的关系进行了大量的工作。大地构造环境是控制花岗岩成矿元素组合的重要因素(汪雄武和王晓地, 2002), Pirajno et al.(1992)曾探讨不同大地构造背景下的与花岗岩有关的矿床(GRMD) 类型, 决定花岗岩成矿类型的重要因素包括花岗岩的类型、演化、分异程度及氧化还原条件等。Blevin and Chappell (1992, 1995), Blevin et al.(1996)在系统总结澳大利亚东部1000多个与花岗岩有关的成矿类型基础上提出, Sn矿化与S型花岗岩或还原的、并经过分离结晶的Ⅰ型花岗岩有关;铜和金与磁铁矿型或/和含榍石的、氧化的、中性Ⅰ型花岗岩有关;Mo与铜金相似, 但相关的岩浆分异更强、更氧化;W与多种花岗岩类型有关, 与岩浆的氧化还原状态的相关性不明显。此外, 具有偏向于高硅含量的花岗岩省一般与锡、钼、钨和金矿化有关, 而偏向于中等或低硅含量的花岗岩一般与铜矿化有关(Blevin and Chappell, 1992)。从矿区主要侵入岩单元的岩石化学看, 总体具有高硅的特征。

野外观察表明, 矿区西北部的矽卡岩型铁铜矿化主要产在花岗闪长(斑) 岩与围岩的接触带处, 矽卡岩类型以石榴子石矽卡岩为主, 而矿区东部的矽卡岩铅锌矿化主要产于细粒花岗岩与围岩的接触带, 以绿帘石矽卡岩为主, 岩体与围岩接触部位可见少量内矽卡岩发育, 为细粒花岗岩与铅锌矿化的成因联系提供了可靠的野外证据。另外, 从矿区范围内已知矿化的空间分布上看, 古磁铁矿采坑或采洞及目前勘查证实的磁铁矿化(含铜、锌) 主要分布于矿区西侧, 花岗岩闪长岩分布区内或外围, 而古采银遗迹和勘查证实的铅锌(银)、钼矿化主要集中在矿区东部, 在细粒花岗岩体内、细粒花岗岩与围接触带或其围岩碳酸盐地层内, 可见岩体与碳酸盐围岩接触带的矽卡岩型和碳酸盐交代型铅锌矿化及黑云母二长花岗岩中的脉状、网脉状及浸染状铅锌矿化。区内已知的钼矿化多呈细脉状、网脉状、浸染状, 主要发现于细粒花岗岩西侧的黑云母二长花岗岩中。花岗闪长岩常具斑状结构, 而细粒花岗岩为微细粒-斑状结构, 均显示了浅成侵位的特征。从矿区蚀变特征看, 矿区范围内以绿泥石化、黄铁矿化、硅化最为发育, 在细粒花岗岩及围岩黑云母二长花岗岩中发育大规模绿泥石、黄铁矿化蚀变, 局部发育强的呈线状分布的石英绢云母化, 目前钻孔控制的绿泥石、黄铁矿化延深达500余米, 大规模发育的热液蚀变表明岩浆富挥发分。

岩石微量元素组成表明, 花岗闪长岩和石英二长斑岩具有最高的Cu+Zn及Sr和Ba含量, 而细粒花岗岩具有最高的W和Cr含量和较高的Sr、Ba含量。细粒花岗岩中Cu、Zn变化范围较大, 其中一个特高值可能是由于样品中轻微矿化造成的。细粒花岗岩相对于花岗闪长岩分异程度和氧化程度均更高, 在Rb/Sr-Fe2O3/FeO图解(图 7) 其投影点即不在典型的成铜金的强氧化花岗岩区域内, 也与典型的成钨锡的壳源花岗岩不同, 而与世界范围内成W、W-Mo的花岗岩投影区域相近。矿区方铅矿铅同位素组成与细粒花岗岩的铅同位素组成接近(李应栩等, 2012), 也表明二者具有成因联系。

图 7 矿区花岗质岩石单元的Rb/Sr-Fe2O3/FeO图解 Fig. 7 Rb/Sr-Fe2O3/FeO diagram of granitic rocks in Anji deposit

上述地质和地球化学特征表明, 矿区主要经历了两期的成矿事件, 第一期是与花岗闪长岩侵位有关的铁(铜、锌) 矿化, 其成矿时代应为晚侏罗世晚期, 年龄在137Ma左右。第二期是与细粒花岗岩侵位有关的铅锌(银、钼) 矿化, 其成矿时代在晚侏罗世-早白垩世之间, 约134Ma左右。这些成矿年龄明显不同于钦杭结合带中段的德兴矿集区等。斑岩铜矿德兴矿区与成矿有关的花岗闪长斑岩锆石SHRIMP U-Pb年龄为171±3Ma (王强等, 2004);江西熊家山钼矿床的辉钼矿Re-Os等时线年龄为152±2.0Ma (孟祥金等, 2007);皖南东源含钨、钼花岗闪长斑岩的锆石U-Pb年龄为146±1Ma (周翔等, 2011), 因此其代表了钦杭成矿带上又一重要的成矿事件。从其成矿时代上看, 其与长江中下游第一期和第二期成矿作用时间(周涛发等, 2008, 2010) 较为一致。

4.3 对矿区找矿勘查的指示意义

由区内主要侵入岩单元的空间分布及与矿化的关系可以看出, 矿区与成矿有关的岩浆侵位主要分布于矿区的东西两侧, 总体呈不规则的椭圆状, 长轴方向为NE, 表明矿区内主要的找矿方向应重点放在矿区东西两侧。矿区内矿化类型多样, 但总体分布于西部花岗闪长岩岩体内及其与岩体接触带以及矿区东部细粒花岗岩岩体内、岩体与碳酸岩地层接触带、岩体周边围岩中。从初步的围岩蚀变研究结果看, 矿区西部岩体以石榴子石矽卡岩化、硅化和黑云母化为主, 岩体多为中-粗粒结构, 表明其剥蚀深度已较大, 进一步的找矿空间较小, 与西部石榴石矽卡岩铁锌铜矿化具有类似矿物组合的五庙桥西部, 上部主要为寒武系地层覆盖, 局部出露小岩株, 剥蚀程度较浅, 在河边已发现高品位铁闪锌矿露头, 往深部寻找类似隐伏矿体具有一定潜力。矿区东部银山附近发育较强的石英绢云母化, 其外围为大面积的绿泥石化、黄铁矿化、硅化, 其中绿泥石化分布范围大, 初步控制约3~5km2, 钻孔控制的蚀变延深约500m, 说明该细粒花岗岩造成的蚀变强度和范围均较大, 且剥蚀程度较浅, 深部仍有一定的找矿空间, 若能进一步对各种蚀变、矿化的空间展布进行详细工作, 确定其热液活动中心, 并在其附近实施进一步的勘查工程, 将有望取得找矿突破。另外, 矿区东部和东北部有大面积的寒武系碳酸岩地层, 在碳酸盐地层中加强碳酸岩交代型铅锌矿化的勘查也是很重要的。

5 结论

(1) 矿区岩浆活动强烈, 活动时限从141~117Ma, 岩浆侵位的顺序为黑云母二长花岗岩、正长花岗岩、花岗闪长岩、细粒花岗岩、双峰式脉岩组合。构造体制转换所造成先期断裂活化、压力降低、地幔物质上涌和高热流值是控制岩浆形成和侵位的重要机制, 而断裂被岩浆侵入焊结、岩浆封存演化、再次开裂和焊结的多次循环, 造成了矿区杂岩体内多期次、成分不同的岩石单元形成。

(2) 矿区岩浆活动是中国东部岩石圈活化及大规模成岩成矿作用的产物, 其构造背景代表了中国东部从压扭性向张扭性转化的过程, 双峰式脉岩组合可以代表区域应力场进入引张的标志。

(3) 矿区主要经历了两次重要的成矿事件, 其对应的成矿时代分别为137Ma和134Ma左右, 两者相隔约3Ma左右, 分别形成了与花岗闪长岩有关的铁(铜、锌) 矿化和与细粒花岗岩有关的铅锌(银)、钼矿化。

(4) 矿化的空间分布主要受与成矿有关的侵入体及构造的双重控制。

(5) 矿区西部花岗闪长岩体内部和外围是寻找铁铜(锌) 矿化的有利地区, 而东部细粒花岗岩体内及其与寒武系灰岩接触部分是寻找铅锌矿化的有利地区, 构造对矿化蚀变具有明显的控制作用, 且成矿后构造对矿体保存影响较大, 应进一步加强构造地质方面的研究工作。

致谢 在野外地质调研过程中得到南京地质矿产研究所王爱国、贺菊瑞等, 铜陵金九矿业集团公司、浙江省第一地质大队及港口铅锌矿普查项目组全体同仁的大力帮助;岩石化学分析过程中得到核工业地质研究院刘牧等人的大力帮助;在此一并致谢!
参考文献
[] Boynton WV. 1984. Geochemistry of the rare earth elements: Meteorite studies. In: Henderson P (ed.). Rare Earth Element Geochemistry. Amsterdam: Elsevier, 63-114
[] Batchelor RA, Bowden P. 1985. Petrogenetic interpretation of granitoid rock series using multicationic parameters. Chemical Geology, 48(1-4): 43–55. DOI:10.1016/0009-2541(85)90034-8
[] Blevin PL, Chappell BW. 1992. The role of magma sources, oxidation states and fractionation in determining the granite metallogeny of eastern Australia. Transactions of the Royal Society of Edinburgh: Earth Sciences, 83(1-2): 305–316. DOI:10.1017/S0263593300007987
[] Blevin PL, Chappell BW. 1995. Chemistry, origin, and evolution of mineralized granites in the Lachlan fold belt, Australia: The metallogeny of Ⅰ-and S-type granites. Economic Geology, 90(6): 1604–1619. DOI:10.2113/gsecongeo.90.6.1604
[] Blevin PL, Chappell BW, Allen CM. 1996. Intrusive metallogenic provinces in eastern Australia based on granite source and composition. Transactions of the Royal Society of Edinburgh: Earth Science, 87(1-2): 281–290. DOI:10.1017/S0263593300006684
[] He GJ, Yang XC, Wu GM, Zhang GF, Cai XX, Zheng J. 2011. A study of ore mineral characteristics and metallogenic stages of the Yinshan Ag-Pb-Zn polymetallic ore deposit, Northwest Zhejiang Province. Acta Geoscientica Sinica, 32(3): 304–312.
[] Huan WL, Shi ZL, Yan JQ. 1982. Meso-Cenozoic tectonic evolution of eastern China and adjacent area and movement of the Pacific plate. Chinese Journal of Geology(2): 179–190.
[] Li ST, Yang SG, Wu CL, Huang JF, Cheng ST, Xia WC, Zhao GR. 1987. Mesozoic taphrogeny in northeastern China and rift basin system in northeastern Asia. Science in China (Series B)(2): 185–195.
[] Li YX, Xie YL, Tang YW, Wang AG, Qiu LM, Zhang XX. 2012. Primary study on the source of lead in Gangkou polymetallic deposit, Anji, Zhejiang Province, China: Lead-isotope evidence. Mineral Deposits, 31(4): 917–929.
[] Lu YF. 2004. GeoKit-A geochemical toolkit for Microsoft Excel. Geochimica, 33(5): 460–464.
[] Mao JW, Wang ZL. 2000. A preliminary study on time limits and geodynamic setting of large-scale metallogeny in East China. Mineral Deposits, 19(4): 289–296.
[] Meng JX, Hou ZQ, Dong GY, Liu JG, Qu WJ, Yang ZS, Zuo LY, Wan LJ, Xiao MZ. 2007. The geological characteristics and Re-Os isotope age of molybdenite of the Xiongjiashan molybdenum deposit, Jiangxi Province. Acta Geologica Sinica, 81(7): 946–951.
[] Mcdonough WF, Sun SS. 1995. The composition of the Earth. Chemical Geology, 120(3-4): 223–253. DOI:10.1016/0009-2541(94)00140-4
[] Pearce JA, Harris NBW, Tindle AG. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. Journal of Petrology, 25(4): 956–983. DOI:10.1093/petrology/25.4.956
[] Pirajno F. 1992. Hydrothermal Mineral Deposits: Principles and Fundamental Concepts for the Exploration Geologist. Berlin, Heidelberg: Springer-Verlag: 1-709.
[] Ren JS. 1984. The indosinian orogeny and its significance in the tectonic evolution of China. Bulletin of the Chinese Academy of Geological Sciences(2): 31–44.
[] Ren JY, Hu XY, Zhang JX. 1998. The Late Mesozoic tectonic activation in the eastern Chinese continent and its evolution process. Geotectonica et Metallogenia, 22(2): 89–96.
[] Tang YW, Xie YL, Li YX, Qiu LM, Liu BS, Li Y, Zhang XX, Jiang YC, Han YD. 2012a. Petrogeochemistry, petrography characteristics and petrogenesis of Wushanguan complex body in Anji mining area, Zhejiang Province. Mineral Deposits, 31(4): 903–916.
[] Tang YW, Xie YL, Li YX, Qiu LM, Liu BS, Li Y, Zhang XX, Jiang YC, Han YD, Zhou JJ. 2012b. Geochemical characteristics and zircons LA-ICP-MS U-Pb age of quartz-monzonite porphyry in Anji polymetallic deposit of Zhejiang and their geological implications. Geoscience, 26(4): 647–655.
[] Wang Q, Zhao ZH, Jian P, Xu JF, Bao ZW, Ma JL. 2004. SHRIMP zircon geochronology and Nd-Sr isotopic geochemistry of the Dexing granodiorite porphyries. Acta Petrologica Sinica, 20(2): 315–324.
[] Wang WX, Wang XD. 2002. Some diagnostic criteria for mineralized granite. Acta Petrologica et Mineralogica, 21(2): 119–130.
[] Wu FY, Ge WC, Sun DY, Guo CL. 2003. Discussions on the lithospheric thinning in eastern China. Earth Science Frontiers, 10(3): 51–60.
[] Xie YL, Tang YW, Li YX, Qiu LM, Liu BS, Li Y, Zhang XX, Han YD, Jiang YC. 2012. Petrogeochemistry characteristics and LA-ICP-MS zircons U-Pb dating of fine grained granite and its geological significance in Anji polymetallic deposit, Zhejiang Province. Mineral Deposits, 31(4): 891–902.
[] Xu WL, Wang QH, Wang DY, Pei FP, Gao S. 2004. Processes and mechanism of Mesozoic lithospheric thinning in eastern North China Craton: Evidence from Mesozoic igneous rocks and deep-seated xenoliths. Earth Science Frontiers, 11(3): 309–317.
[] Xu ZQ, Zhang QD, Zhao M. 1982. The main characteristics of the paleorift in the middle section of the Tancheng-Lujiang fracture zone. Bulletin of the Chinese Academy of Geological Sciences, 4(1): 17–44.
[] Yang MG, Mei YW. 1997. Characteristics of geology and metallization in the Qinzhou-Hangzhou paleoplate juncture. Geology and Mineral Resources of South China(3): 52–59.
[] Yang MG, Huang SB, Lou FS, Tang WX, Mao SB. 2009. Lithospheric structure and large-scale metallogenic process in Southeast China continental area. Geology in China, 36(3): 528–543.
[] Yang SS. 1991. A preliminary analysis on metallogenetic conditions of polymetallic ore deposits in southern Anhui-western Zhejiang area. Geology and Prospecting(5): 13–17.
[] Zhou TF, Fan Y, Yuan F. 2008. Advances on petrogensis and metallogeny study of the mineralization belt of the Middle and Lower Reaches of the Yangtze River area. Acta Petrologica Sinica, 24(8): 1665–1678.
[] Zhou TF, Fan Y, Yuan F, Song CZ, Zhang LJ, Qian CC, Lu SM, David RC. 2010. Temporal-spatial framework of magmatic intrusions in Luzong volcanic basin in East China and their constrain to mineralizations. Acta Petrologica Sinica, 26(9): 2694–2714.
[] Zhou X, Yu XQ, Wang DN, Zhang DH, Li CL, Fu JZ, Dong HM. 2010. Characteristics and geochronology of the W, Mo-bearing granodiorite porphyry in Dongyuan, southern Anhui. Geoscience, 25(4): 201–210.
[] 何国锦, 杨晓春, 吴光明, 张国防, 蔡雄翔, 郑剑. 2011. 浙江西北银山银铅锌多金属矿床矿石矿物特征及成矿期次初步研究. 地球学报, 32(3): 304–312.
[] 环文林, 时振梁, 鄢家全. 1982. 中国东部及邻区中新生代构造演化与太平洋板块运动. 地质科学(2): 179–190.
[] 李思田, 杨士恭, 吴冲龙, 黄家福, 程守田, 夏文臣, 赵根榕. 1987. 中国东北部晚中生代裂陷作用及东北亚断陷盆地系. 中国东北部晚中生代裂陷作用及东北亚断陷盆地系.中国科学(B辑)(2): 185–195.
[] 李应栩, 谢玉玲, 唐燕文, 王爱国, 邱立明, 张欣欣. 2012. 浙江安吉港口多金属矿床铅物质来源初探--铅同位素证据. 矿床地质, 31(4): 917–929.
[] 路远发. 2004. GeoKit:一个用VBA构建的地球化学工具软件包. 地球化学, 33(5): 460–464.
[] 毛景文, 王志良. 2000. 中国东部大规模成矿时限及其动力学背景的初步探讨. 矿床地质, 19(4): 289–296.
[] 孟祥金, 侯增谦, 董光裕, 刘建光, 屈文俊, 杨竹森, 左力艳, 万禄进, 肖茂章. 2007. 江西金溪熊家山钼矿床特征及其Re-Os年龄. 地质学报, 81(7): 946–951.
[] 任纪舜. 1984. 印支运动及其在中国大地构造演化中的意义. 中国地质科学院院报(2): 31–44.
[] 任建业, 胡祥云, 张俊霞. 1998. 中国大陆东部晚中生代构造活化及其演化过程. 大地构造与成矿学, 22(2): 89–96.
[] 唐燕文, 谢玉玲, 李应栩, 邱立明, 刘保顺, 李媛, 张欣欣, 韩宇达, 姜妍岑. 2012a. 浙江安吉坞山关杂岩体岩相学、岩石地球化学特征及岩石成因. 矿床地质, 31(4): 903–916.
[] 唐燕文, 谢玉玲, 李应栩, 邱立明, 刘保顺, 李媛, 张欣欣, 姜妍岑, 韩宇达, 周俊杰. 2012b. 浙江安吉多金属矿区石英二长斑岩岩石化学特征和锆石U-Pb年龄及其地质意义. 现代地质, 26(4): 647–655.
[] 王强, 赵振华, 简平, 许继峰, 包志伟, 马金龙. 2004. 德兴花岗闪长斑岩SHRIMP锆石U-Pb年代学和Nd-Sr同位素地球化学. 岩石学报, 20(2): 315–324.
[] 汪雄武, 王晓地. 2002. 花岗岩成矿的几个判别标志. 岩石矿物学杂志, 21(2): 119–130.
[] 吴福元, 葛文春, 孙德有, 郭春丽. 2003. 中国东部岩石圈减薄研究中的几个问题. 地学前缘, 10(3): 51–60.
[] 谢玉玲, 唐燕文, 李应栩, 邱立明, 刘保顺, 李媛, 张欣欣, 韩宇达, 姜妍岑. 2012. 浙江安吉多金属矿床细粒花岗岩岩石化学特征、锆石U-Pb定年及其与成矿关系探讨. 矿床地质, 31(4): 891–902.
[] 许文良, 王清海, 王冬艳, 裴福萍, 高山. 2004. 华北克拉通东部中生代岩石圈减薄的过程与机制:中生代火成岩和深源捕虏体证据. 地学前缘, 11(3): 309–317.
[] 许志琴, 张巧大, 赵民. 1982. 郯庐断裂中段古裂谷的基本特征. 中国地质科学院院报, 4(1): 17–44.
[] 杨明桂, 梅勇文. 1997. 钦-杭古板块结合带与成矿带的主要特征. 华南地质与矿产(3): 52–59.
[] 杨明桂, 黄水保, 楼法生, 唐维新, 毛素斌. 2009. 中国东南陆区岩石圈结构与大规模成矿作用. 中国地质, 36(3): 528–543.
[] 杨松生. 1991. 皖南浙西地区多金属矿成矿条件浅析. 地质与勘探(5): 13–17.
[] 周涛发, 范裕, 袁峰. 2008. 长江中下游成矿带成岩成矿作用研究进展. 岩石学报, 24(8): 1665–1678.
[] 周涛发, 范裕, 袁峰, 宋传中, 张乐骏, 钱存超, 陆三明, DavidRC. 2010. 庐枞盆地侵入岩的时空格架及其对成矿的制约. 岩石学报, 26(9): 2694–2714.
[] 周翔, 余心起, 王德恩, 张德会, 李春麟, 傅建真, 董会明. 2011. 皖南东源含W、Mo花岗闪长斑岩及成矿年代学研究. 现代地质, 25(4): 201–210.