2012, Vol. 28
2. 武警黄金地质研究所，廊坊 065000;
3. 武警黄金一支队，牡丹江 157000
, Yang LiQiang1, Chen YongFu2, Qing Min2, Yan JiaPan2, Ge LiangSheng2, Guo XiaoDong2, Wang JingRui3
2. Gold Geological Institute of CAPF, Langfang 065000, China;
3. No.1 Party of Gold Geology, CAPF, Mudanjiang, 157000, China
埃达克岩最早由Kay(1978) 在美国阿留申群岛中的Adak岛发现，由Defant and Drummod(1990) 命名，是具有特定地球化学特征的一套中酸性火山岩和侵入岩组合。钱青(2001) 综合研究前人有关资料后，认为它是一类高铝的中酸性火山岩或侵入岩。张旗等(2002) 指出，埃达克岩不是指某一种具体的岩石类型，而是指具有特定地球化学性质的一套中酸性火山岩和侵入岩组合。张旗等(2004) 依埃达克岩的多样性进一步将埃达克岩分为两类6种类型。
黑龙江东宁金厂铜金矿床位于吉黑东部的延边-东宁地区，该区兴蒙造山运动在古生代末期(250Ma)结束，印支晚期(210Ma)进入环太平洋运动阶段(Wu et al., 2011)。受太平洋板块俯冲影响，区域为挤压造山构造环境，I型花岗岩岩浆活动发育。燕山晚期(110~120Ma)，区域构造环境变为弧后(陆缘弧)伸展(Wu et al., 2011)。据地震资料证实(高立新，2011；Zhao et al., 2011；黄金莉，2010)，此时俯冲板片前缘已越过壳幔边界，从而引起壳幔物质熔融，导致区域上中酸性岩脉广泛发育，其中有一类具埃达克质岩石的特点(张旗等，2001；李印等，2009)。该类岩石在东北地区分布较广，且与成矿作用关系密切。张炯飞等(2004) 统计我国东北地区与埃达克岩有关的矿床已有10个之多，包括矿床类型有斑岩型、岩浆热液型、浅成低温型，矿种主要为Au和Cu，次为Mo、Ag、Pb和Zn，成矿环境多为火山弧环境。
东宁金厂铜金矿床中J-0号矿体深部闪长玢岩特征与埃达克岩相似，闪长玢岩与J-0号矿体时空关系密切，是J-0号矿体的成矿地质体，可以作为找矿的标志。为此，本文通过岩相学、岩石地球化学和U-Pb同位素定年对闪长玢岩的成因及成矿时代进行分析与探讨，以期对金厂矿区的成矿地质体的物质来源、岩浆演化、成岩时代及成矿地质背景进行初步探讨。
2 地质背景金厂铜金矿床位于古亚洲构造域与环太平洋构造域复合部位，兴凯地块西南部，位于Ⅲ级构造单元张广才岭-太平岭边缘隆起带中的太平岭隆起与老黑山断陷结合部位，受双重大地构造单元的制约和影响，成矿地质背景复杂。
金厂铜金矿床矿区内大面积分布中生代侵入岩，主要为印支期-燕山期中酸性侵入岩，地层分布很少，仅在矿区外围出露上元古界黄松群变质岩系，矿区东部及外围零星出露中生代火山岩系，岩性为流纹质-英安质岩屑晶屑凝灰岩-流纹质火山角砾岩(图 1)。
|
图 1 金厂矿区地质简图 1-第四系；2-第三系；3-三叠系；4-燕山晚期花岗斑岩；5-燕山晚期闪长玢岩；6-燕山早期闪长岩；7-印支晚期-燕山早期粗粒花岗岩；8-印支晚期-燕山早期中细粒花岗岩；9-印支晚期-燕山早期文象花岗岩；10-环形构造；11-断裂；12-推断断裂；13-环状矿体；14-放射状矿体；15-角砾岩型矿体 Fig. 1 Sketch geological map of Jinchang deposit 1 -Quaternary; 2-Tertiary; 3-Triassic; 4-granite porphory of Late Yanshannia; 5-diorite porphory of Late Yanshannia; 6-diorite of Early Yanshannia;7-Late Indonisia and Early Yanshannia coarse grain granites; 8-Late Indonisia and Early Yanshannia median to fine grain granites; 9- Late Indonisia and Early Yanshannia graphic granites; 10-circularity structure; 11-faults; 12-sespected faults; 13-circularity ore body; 14-actimophous ore body; 15-breccia ore body |
区内构造运动和岩浆活动的频繁，构造主要为断裂和穹隆。矿区北西侧与北东向绥阳断裂相邻，该断裂控制了区域地层和印支-燕山期侵入岩体的展布，是太平岭隆起带与老黑山断陷的分界断裂。绥阳断裂的长期活动派生出一系列次级断裂，使区内断裂构造十分发育，形成了北东向、东西向、南北向、北西向的线形断裂，以及由这些断裂构造控制的隐爆角砾岩筒。矿区内环线构造及放射状构造发育，其形成与深部侵入体有关。
矿区岩浆岩主要有两期岩浆活动：印支晚期-燕山早期岩浆活动和燕山晚期岩浆活动。印支晚期-燕山早期岩浆岩主要为花岗岩、花岗闪长岩，呈岩株状，为矿体的主要围岩；燕山晚期岩浆岩主要为闪长玢岩和花岗斑岩岩脉，呈脉状和小岩枝状，与矿体时空关系密切。
金厂矿区现已发现大小金矿体二十几个，金矿体矿化类型存在多样性。矿体类型划分为三类五亚类：角砾岩型、斑岩型和岩浆穹隆型；其中角砾岩型矿体又可划分为隐爆角砾岩型、热液充填角砾岩型和塌陷角砾岩型矿体三亚类；岩浆穹隆型矿体划分为环状裂隙充填型和放射状裂隙充填型两亚类。
J-0号矿体为复合型矿体，上部为热液充填角砾岩型矿体(呈脉状产于隐爆角砾岩筒中)，中部为坍塌角砾岩型矿体(产于隐爆角砾岩筒与闪长玢岩过渡带中)，深部为斑岩型矿体(产于闪长玢岩岩体中)。本文主要对J-0号矿体深部的闪长玢岩岩体进行分析讨论。
3 闪长玢岩地质-地球化学特征矿区内闪长玢岩位于高丽沟J-0号矿体深部，地表仅零星出露。闪长玢岩为一小岩枝，目前控制岩枝最大直径约100m，控制延伸200m。小岩枝侵位于印支晚期正长花岗岩中，从闪长玢岩岩体中心通过角砾岩带至正长花岗岩，总体变化趋势为：未蚀变或弱蚀变闪长玢岩→蚀变闪长玢岩→闪长玢岩角砾岩带→正长花岗岩角砾岩带→蚀变正长花岗岩→正长花岗岩。
3.1 岩相学特征本次工作对J-0号矿体的一至五中段闪长玢岩进行了野外地质调查及室内的岩相学分析及研究。闪长玢岩自岩体中心向外蚀变依次变强，矿化强度与蚀变成正相关。本文依其蚀变强度进行了闪长玢岩的岩相学的相关分析(表 1、图 2)。
|
表 1 闪长玢岩岩相学特征表 Table 1 Petrography characterstics of diorite porphyry |
|
图 2 闪长玢岩显微特征照片 Fig. 2 Microphotograph characteristics of diorite porphyry |
笔者依闪长玢岩的蚀变程度不同，共采取主微量分析样品22件，未蚀变闪长玢岩6件，弱蚀变闪长玢岩9件，强蚀变闪长玢岩7件。样品由河北省区域地质矿产调查研究所实验室测试。
金厂闪长玢岩的岩石化学成分如表 2所示。高丽沟未蚀变闪长玢岩的主量元素特征显示：为中性，SiO2含量54.62%~58.62%，平均57.24%；富钠贫钾，Na2O含量3.86%~5.25%，平均4.79%，K2O含量0.37%~0.69%，平均0.51%；富钙，CaO含量6.93%~8.35%，平均7.32%；富铝，Al2O3含量15.08%~18.25%，平均17.18%；镁含量中等，MgO含量2.70%~5.68%，平均4.11%； 铁含量中等，Fe2O3含量0.84%~1.93%，平均1.27%；FeO含量2.28%~3.74%，平均3.01%。
|
|
表 2 闪长玢岩的主量元素(wt%)微量元素(×10-6)含量及相关参数 Table 2 Major elements (wt%) and trace elements (×10-6) compositions and related parameters of diorite porphyry |
弱蚀变闪长玢岩显示：中酸性，SiO2含量55.38%~66.38%，平均60.75%；Na2O含量3.95%~5.31%，平均4.90%，K2O含量0.30%~1.61%，平均0.83%；CaO含量3.84%~7.95%，平均5.68%；Al2O3含量15.54%~17.51%，平均16.63%；MgO含量1.87%~4.41%，平均3.44%；Fe2O3含量0.48%~2.45%，平均1.56%；FeO含量1.84%~3.52%，平均2.72%。
强蚀变闪长玢岩显示：偏酸性，SiO2含量55.99%~67.82%，平均62.42%；Na2O含量4.56%~5.35%，平均5.03%；K2O含量0.47%~1.11%，平均0.78%；CaO含量3.37%~7.09%，平均5.27%；Al2O3含量14.52%~16.76%，平均16.12%；MgO含量1.66%~5.55%，平均2.72%；Fe2O3含量0.26%~1.85%，平均0.86%；FeO含量1.58%~2.85%，平均2.21%。
从闪长玢岩的主量元素含量可以看出：SiO2含量随着玢岩蚀变程度的升高而增高，如硅-碱图所示(图 3)；Al2O3含量随着玢岩蚀变程度的升高而降低；MgO含量随着玢岩蚀变程度的升高而降低；Fe2O3含量随着玢岩蚀变程度的升高而先增高后降低；FeO含量随着玢岩蚀变程度的升高而降低。其原因有三点：(1) 闪长玢岩蚀变的强度与正长花岗岩角砾的含量正相关，酸性正长花岗岩角砾的加入使中性闪长玢岩的硅铝(SiO2、Al2O3、Na2O)含量增加，而镁铁(MgO、FeO、CaO)含量减少；(2) 闪长玢岩的蚀变使角闪石和斜长石生成绿泥石化、碳酸盐化、绢云母化和高岭土化等次生矿物，这些次生矿物易溶于水，从而使镁钙含量等减少；(3) 闪长玢岩的蚀变使FeO易被流体中的氧氧化，生成Fe2O3，使闪长玢岩中FeO的含量减少而Fe2O3的含量增加。
|
图 3 闪长玢岩硅-碱图解(据Middlmost, 1994) Ir-Irvine 分界线，上方为碱性，下方为亚碱性.1-橄榄辉长岩；2a-碱性辉长岩；2b-亚碱性辉长岩；3-辉长闪长岩；4-闪长岩；5-花岗闪长岩；6-花岗岩；7-硅英岩；8-二长辉长岩；9-二长闪长岩；10-二长岩；11-石英二长岩；12-正长岩；13-副长石辉长岩；14-副长石二长闪长岩；15-副长石二长正长岩；16-副长正长岩；17-副长深成岩；18-霓方钠岩/磷霞岩/粗白榴岩.第1组-未蚀变闪长玢岩；第2组-弱蚀变闪长玢岩；第3组-强蚀变闪长玢岩；后图同 Fig. 3 SiO2-(Na2O+K2O)diagram of diorite porphyry(after Middlmost, 1994) Ir-irvine line, the above is alkaline, below is subalkaline. 1-olivine gabbro; 2a-alkaline gabbro; 2b-subalkaline gabbro 3-gabbro diorite; 4-diorite; 5-granodiorite; 6-granite; 7-silicon rock; 8- monzonitic gabbro; 9-monzonitic diorite; 10-monzonite; 11-quartz monzonite; 12-syenite; 13-feldspathoid gabbro; 14-feldspathoid monzonitic diorite; 15-feldspathoid monzonite; 16-feldspathoid syenite; 17-feldspathoid plutonic; 18-Ni-sodalite/P-nepheline/Leucite; Group I-unaltered diorite porphyry;Group II-weakly altered diorite porphyry; Group III-strongly altered diorite porphyry |
在AR-SiO2图解中(图 4a)，闪长玢岩除两个点外所有样品均落在钙碱性区域内，而在亚碱性岩石进一步分类常用的K2O-SiO2(图 4b)图解上，样品点大多数落在低钾拉斑系列。
|
图 4 闪长玢岩AR-SiO2图解(a，据Wright, 1969)和K2O-SiO2图解(b，实线据 Peccerillo and Taylor, 1976；虚线据 Middlemost, 1985) Fig. 4 AR-SiO2 diagram(a, after Wright, 1969)and K2O-SiO2 diagram(b, solid line after Peccerillo and Taylor, 1976; dash line after Middlemost, 1985)of diorite porphyry |
闪长玢岩与Defant and Drummond(1990) 所定义的典型埃达克岩相比，具有很好的稳合性，如SiO2 ≥56%、Al2O3≥15%(很少低于此值)和MgO通常小于3%(极少大于6%)，K2O/Na2O比值＜0.5。
3.3 微量元素特征闪长玢岩的标准化顺序和数值均采用Sun and McDonough(1989) 标准化顺序和数值。
闪长玢岩的微量元素测试结果如表 2所示，微量元素(不相容元素)的原始地幔标准化蛛网图(图 5a)，从图中可以看出，未蚀变、弱蚀变和强蚀变闪长玢岩的蛛网图曲线一致，说明闪长玢岩的蚀变对微量元素的含量影响不大。
|
图 5 闪长玢岩原始地幔标准化的微量元素蛛网图(a)和稀土元素球粒陨石标准化配分图(b) Fig. 5 Primitive mantle-normalized trace element spider diagram(a)and chondrite-normalized REE distribution patterns of diorite porphyry(b) |
结果显示相容元素Co、Ni和Cr含量较低。而原始地幔标准化的微量元素蛛网图(图 5a)则表明研究区闪长玢岩以大离子亲石元素(LILE)相对富集，部分高场强元素(HSF)相对亏损为特征。其中Pb富集而Nb、Ta、P和Ti相对亏损，Zr、Hf异常不明显，U相对于Th富集；且Nb/Ta值主要变化于11~19之间，部分高于19，具明显的幔源特征(Sun and McDonough, 1989)。
新鲜样品具有高Sr(433.4×10-6~606.7×10-6)、低Y(9.69×10-6~11.72×10-6)特征，Sr含量高(很少小于400×10-6)，Sr/Y＞20~40，La/Yb≥10。说明闪长玢岩表现出显著的埃达克质岩石特性。
3.4 稀土元素特征闪长玢岩稀土元素测试结果如表 2所示，从球粒陨石标准化稀土配分曲线上看(图 5b)，明显具右倾趋势。轻重稀土分馏明显，轻稀土相对富集，重稀土相对亏损，LR/HR普遍较高，集中在5.06~14.24之间，平均9.16。稀土总量中等，∑REE在46.01×10-6~136.8×10-6之间，平均86.76×10-6，稀土总量具有从未蚀变闪长玢岩(∑REE在80.23×10-6~122.9×10-6之间，平均103.5×10-6)→弱蚀变闪长玢岩(∑REE在46.01×10-6~136.8×10-6之间，平均84.14×10-6)→强蚀变闪长玢岩(∑REE在48.63×10-6~102.6×10-6之间，平均74.36×10-6)逐渐降低的趋势，可能是岩石蚀变造成的稀土流失。除一个点为弱负异常(δEu为0.85)外其它样品都表现出Eu的正异常(δEu为1.0~1.46之间，平均1.21)，表明在岩浆分离结晶过程中斜长石不占据元素总分配系数的主导地位，Rollison(1993) 认为在中性熔浆体系中，角闪石对REE是相容的，并在Dy和Er之间具有最高的分配系数值，辉石和榍石也具有这种效应，闪长玢岩中角闪石和辉石能占到10%~15%，表明在闪长玢岩结晶分异过程中角闪石和辉石占据稀土元素总体系数的主导地位。轻稀土分馏程度大于重稀土，(La/Sm)N为3.06~7.0，(Gd/Yb)N为1.57~3.25。Yb含量1.19×10-6~1.41×10-6，较高的Sr/Y比值(41.30~59.30)，La/Yb比值主要介于(10.74~16.69)之间。Y和重稀土元素含量低(Y＜18×10-6、Yb≤1.9×10-6)，闪长玢岩表现出显著的埃达克质岩石特性。
4 锆石U-Pb年龄 4.1 样品采集及分析流程在野外地质观察与研究的基础上，选择了高丽沟J-0号矿体五中段的赋矿闪长玢岩——蚀变闪长玢岩和岩体中心相——未蚀变闪长玢岩样品为测试对象，样品编号分别为09JB400和09JB401，样品采集重量大于10kg。锆石的分选是在河北廊坊区域地质调查研究所完成，经过手工破碎大约10kg重的样品、淘洗、电磁选、重液分选，之后在双目镜下挑选，得到含包裹体少、无明显裂隙且晶形完好的锆石(图 6a，b)。在双目镜下将锆石样品粘在双面胶上，制成靶，备用。然后依据样品的反射光和透射光照片，排除多裂纹和抛光面不清晰的锆石；再根据锆石阴极发光(CL)图像，进一步完成锆石内部结构的分析。
|
图 6 闪长玢岩中锆石的显微镜照片 Fig. 6 Micrograph of zircons from diorite porphyry |
09JB400号样品的锆石阴极发光(CL)图像是在北京大学完成，锆石的制靶、光学显微镜照相及LA-ICP-MS锆石U-Th-Pb同位素年代学测试在中国地质大学(北京)地质过程与矿产资源国家重点实验室完成。
所用ICP-MS为美国Agilent科技有限公司生产的7500a型，剥蚀系统为美国New Wave贸易有限公司，UP193SS型。分析所用的光斑直径为36μm，年龄计算时以标准锆石91500为外标进行同位素比值校正，标准锆石TEMORA为监控盲样；元素含量以国际标样NIST610为外标，Si为内标计算，NIST612和NIST614为监控盲样。锆石U-Th-Pb年代学测试数据处理采用Isoplot软件进行。204Pb校正方法同Andersen(2002) ，实验室测试详细方法参考Song et al.(2010) 。
09JB401号样品锆石的阴极发光(CL)显微分析在中国科学技术大学物质结构中心实验室完成。锆石原位U-Pb 同位素年龄分析在中国科学技术大学壳幔物质与环境重点实验室的LA-ICP-MS 仪器上用标准测定程序进行。分析仪器为Elan DRC Ⅱ型四极杆质谱仪和GeoLas Pro 193μm 型激光剥蚀系统，激光器为193nm ArF 准分子激光器，激光剥蚀斑束直径为30μm，激光剥蚀样品的深度为20~40μm，并采用He作为剥蚀物质的载气。锆石年龄计算采用标准锆石91500 作为外标，TEM 作为内标，29Si 作为内标元素进行校正。采样方式为单点剥蚀，数据采集选用一个质量峰一点的跳峰方式(peak jumping)。详细的实验原理、流程及仪器参数见翟明国等(2004) 。测试结果通过LaDating@Zrn.xls ver 1.0软件计算得出，并按照Andersen(2002) 的ComPbCon# 32151 程序对其进行了普通铅校正。
4.2 测试结果09JB400样品锆石粒径多为100~300μm，多为长柱状或短柱状，自形程度较高，为自形到半自形。在锆石CL图像上，本次测试的锆石多具有清楚的振荡环带，为典型的岩浆成因锆石特征(图 7a)。部分锆石内部结构不均，可见内核及外边，但其发育振荡环带的主体部分岩浆成因特征明显，振荡环带部分从内到外具有一致的表面年龄。锆石U-Th-Pb年龄测试结果见表 3。
|
图 7 09JB400(a)和09JB401(b)部分锆石阴极发光图像 Fig. 7 Part of the zircon CL images of diorite porphyry sample 09JB400(a)and 09JB401(b) |
|
|
表 3 样品09JB400锆石U-Th-Pb 分析数据 Table 3 Zircon U-Pb analytical results of sample 09JB400 |
09JB401样品锆石粒径多为80~300μm，多为长柱状或短柱状，自形程度较高，为半自形。在锆石CL图像上，本次测试的锆石多具有振荡环带，为典型的岩浆成因锆石特征(图 7b)，其发育振荡环带的岩浆成因特征明显，振荡环带从内到外具有一致的表面年龄。锆石U-Th-Pb年龄测试结果见表 4。
|
|
表 4 样品09JB401锆石U-Th-Pb 分析数据 Table 4 Zircon U-Pb analytical results of sample 09JB401 |
在09JB400号样品测试的34个锆石数据中，样品锆石的Th/U比值，所有测试点都大于0.4，最高达到1.13，表明09JB400中的锆石为典型的岩浆成因锆石。经分析认为，有8个锆石测试数据的协和率较低，其余24颗锆石的年龄较为一致。
前人研究表明，采用激光探针进行U-Pb同位素定年时，需要进行普通铅的校正，对大于1Ga年龄的锆石采用207Pb/206Pb年龄合适，而对于年轻的锆石样品采用206Pb/238U年龄较为合适。由于闪长玢岩成岩年龄小于250Ma，故采用206Pb/238U的年龄，206Pb/238U表面年龄加权平均值为118.1±1.6Ma(MSWD＝0.14，概率为0.95)(图 8a)；26个样品点的表面年龄相关性较好，符合正态分布分布规律。26个样品点在206Pb/238U-207Pb/235U谐和图上表现为较好的谐和性(图 8a)，年龄值较为集中，谐和年龄为118.1±1.6Ma(MSWD＝0.14)。
|
图 8 09JB400(a)和09JB401(b)锆石U-Pb年龄 Fig. 8 Zircon concodia diagrams and U-Pb ages of diorite porphyry sample 09JB400(a)and 09JB401(b) |
在09JB401号样品测试的39个锆石数据中，样品锆石的Th/U比值，所有测试点都大于0.6，最高达到1.42，表明09JB401中的锆石为典型的岩浆成因锆石。经分析认为，有29个锆石测试数据的协和率较低，其余10颗锆石的年龄较为一致。
由于矿区火山岩成岩年龄小于250Ma，故采用206Pb/238U的年龄，206Pb/238U表面年龄加权平均值为115.7±2.0Ma(MSWD＝0.67)(图 8b)；10个样品点的表面年龄相关性较好，10个样品点在206Pb/238U-207Pb/235U谐和图上表现为较好的谐和性(图 8b)，年龄值较为集中，谐和年龄为115.7±2.0Ma(MSWD＝0.67)。
5 讨论 5.1 埃达克质闪长玢岩金厂矿区未蚀变闪长玢岩的SiO2含量54.62%~58.62%，平均57.24%，显示其为中性，与闪长岩的成分相当，结合其岩相及斑状结构特征，确定高丽沟岩体为闪长玢岩，属于低钾拉斑系列和钙碱性系列。闪长玢岩在主、微量和稀土等岩石化学方面均表现出与埃达克岩相似的特点。在闪长玢岩Sr/Y-Y图解(图 9a)和(La/Yb)N-YbN图解(图 9b)中，闪长玢岩样品均落入埃达克岩中，因此，金厂铜金矿床高丽沟闪长玢岩是一种埃达克岩或埃达克质岩。从判别图中(见后文)可以看出闪长玢岩的形成大地构造背景为火山弧环境，结合张旗等(2004) 对埃达克岩的分类方案，认为，金厂闪长玢岩属O型埃达克岩或埃达克质岩。
|
图 9 闪长玢岩Sr/Y-Y图解(a)和(La/Yb)N-YbN图解(b)(据 Defant and Drummond, 1990) Fig. 9 Sr/Y-Y diagram(a)and(La/Yb)N-YbN diagram(b)for diorite porphyry(after Defant and Drummond, 1990) |
前人(Thieblemont and Stein, 1997；Sajona et al., 1998；Bellon and Yumul, 2001；Oyarzun et al., 2001；张旗等，2002；侯增谦等，2003，2005；郭晓东等，2009)研究发现，埃达克岩或埃达克质岩与斑岩矿床和浅成低温热液矿床在时间上、空间上和成因上具有特殊的关系。世界主要的大型-超大型斑岩铜矿主要与板片俯冲形成的O型埃达克岩有关(Thieblemont and Stein, 1997；Oyarzun et al., 2001；张旗等，2009)。侯增谦等(2003) 认为埃达克岩是斑岩铜矿的重要的含矿母岩。Defant et al.(2002) 认为埃达克岩可作为找矿的标志来使用。
O型埃达克岩贫钾，说明它来源于贫钾的玄武质源岩，西太平洋许多O型埃达克岩形成于岛弧环境是人所共知的，那是因为有其它地质证据的佐证或有我们可以看得见的证据(张旗等，2009)。O型埃达克岩贫钾，说明它来源于贫钾的玄武质源岩，产于板块俯冲带(张旗等，2009)，与金厂矿区的埃达克质闪长玢岩产出的构造背景相吻合。
O型埃达克岩最重要的价值在于其与金铜成矿作用的密切关系(张旗等，2009)。金厂铜金矿区闪长玢岩为O型埃达克质岩或埃达克岩，是高丽沟J-0号矿体的容矿岩，J-0号矿体为一斑岩型铜金矿体。与前人提出的关于埃达克质岩与斑岩型铜金关系密切的结论是相符的，金厂矿区埃达克质闪长玢岩的确定对矿区的深部找矿具有重大的指导意义。
5.2 岩浆源区由于矿区及其外围没有与原始地幔性质相近的玄武岩(要求MgO含量大于7%)的产出，以其作为岩浆演化起源的端元组分进行研究无法进行。而结晶分异对母岩浆的不相容元素影响较小，可以作为判别岩浆地幔起源的有效手段。
闪长玢岩原始地幔标准化的微量元素蛛网图(图 5a)显示的Nb、Ta和Ti等亏损，一般认为主要有以下几个原因：①岩浆源区部分熔融过程中残留高场强元素(HFSE)富集的矿物(如金红石和钛铁矿等)；②岩浆上涌过程中受到地壳的混染；③与俯冲流体交代有关。由于岩浆源区部分熔融过程中若有金红石残留，这样起源的岩浆不仅会造成Nb、Ta的亏损，而且会引起明显的Zr、Hf和Ti的亏损。事实表明不能用岩浆源区有富集高场强元素(HFSE)的矿物残留来解释。虽然Pb的富集具有某些大陆地壳的性质，但单纯按地壳物质混染或者壳幔混合熔融很难解释Nb、Ta的低含量，因为金厂金矿区闪长玢岩的Nb、Ta含量甚至低于部分研究人员所计算的陆壳物质中的含量(Rudnick and Gao, 2010)。
金厂闪长玢岩δEu的变化范围为0.85~1.46，平均为1.21，(La/Yb)N平均值为10.18，表明金厂闪长玢岩岩体与壳幔型花岗岩地球化学特征相似。在(La/Yb)N-δEu变异图上，闪长玢岩样品均落入壳-幔型花岗岩范围内(图 10)，也说明本区岩浆主要来源于深部。
|
图 10 闪长玢岩(La/Yb)N-δEu变异图 Fig. 10 (La/Yb)N-δEu diagram for diorite porphyry |
因此，金厂金矿区闪长玢岩的这种岩石地球化学特征与岩浆在源区受俯冲流体交代有关。这种来自大洋地壳和大洋远洋沉积物的俯冲流体以富集大离子亲石元素(LILE)和轻稀土元素(LREE)而亏损高场强元素(HFSE)为特征(MacDonald et al., 2000)。其交代上覆地幔橄榄岩导致其富集大离子亲石元素(LILE)、轻稀土元素(LREE)和亏损高场强元素(HFSE)。同时，研究表明相对于Th而言，U更易于从俯冲板片迁移到地幔楔而相对富集(Ayers, 1998)，这在金厂金矿区闪长玢岩中也有显示。总之，金厂金矿区闪长玢岩可能起源于俯冲流体交代的地幔楔熔融岩浆，在其演化侵位过程中有壳源物质的加入。
闪长玢岩微量元素特征表明其大离子亲石元素(LILE)和轻稀土元素(LREE)相对富集，而高场强元素(HFSE)和重稀土元素(HREE)相对亏损，显示地壳混染特征(Jahn et al., 1999)。极度富集Pb而亏损Nb、Ta具有大陆下地壳的性质(McDonough and Sun, 1995; 高山等，2009；牛耀龄，2010)。因此，闪长玢岩的原始岩浆可能起源与大陆下地壳或者是岩浆上涌过程中受到了下地壳的混染。而所测样品Nb/Ta值主题接近于19且部分高于19，具明显的幔源特征(Sun and McDonough, 1989)。
5.3 岩浆演化闪长玢岩岩石化学特征表明金厂金矿区闪长玢岩应属于亚碱性岩石组合，具有钙碱性演化趋势，但钾含量较低。野外观察和镜下鉴定显示闪长玢岩的斑状结构非常明显，其斑晶主要以斜长石、角闪石为主，并见少量辉石斑晶。表明岩浆演化过程中经历了这些矿物的分离结晶作用。而稀土元素的Eu正异常表明斜长石在岩浆的结晶分异过程中不起主导作用，而与角闪石和辉石的分离结晶作用有关。此外，过度族元素Cr、Co、Ni等主要赋存在橄榄石、辉石等镁铁质矿物中，HREE则主要分配在单斜辉石和角闪石等铁镁质矿物中。金厂矿区闪长玢岩的低含量过度族元素和相对亏损的HREE特征进一步说明岩浆演化过程中结晶分异的存在，并以辉石和角闪石的分异为主。同时闪长玢岩中见有少量磁铁矿副矿物，其结晶分异也能引起HREE的亏损和过度族元素含量的降低，更可能是微量元素Ti弱亏损的主要因素。
在岩体的La/Sm-La关系图上，部分熔融表现为一条斜线，分离结晶过程表现为一条水平线。从图 11上看，金厂闪长玢岩岩体样品点按其蚀变强度不同而具有一定的水平线性关系，且总体趋势也是近于水平的，也说明金厂闪长玢岩岩体的形成以结晶分异作用为主。
|
图 11 闪长玢岩La/Sm-La关系图 Fig. 11 Ce/Yb-Eu/Yb diagram for diorite porphyry |
赵玉锁等(2011①)在研究金厂铜金矿床的成矿地质背景时，在总结前人的研究成果的基础上，认为延边-东宁成矿带的成矿地质背景为“活动大陆边缘的火山弧构造环境”，具备典型的沟-弧-盆体系。从金厂闪长玢岩的Rb-(Y+Nb)判别图和Ta-Yb判别图(图 12a，b)中可以看出，闪长玢岩的所有样品点均落在火山弧花岗岩区，说明闪长玢岩的成岩大地构造背景为大陆边缘的火山弧环境，与前人得出的结论是一致的。
① 赵玉锁，陈永福，阎家盼，肖力，王艳忠，王景瑞，张景海，崔雷，牛翠祎，段晓君，曹勇.2011. 延边-东宁成矿带“金厂式”金矿综合研究及勘查突破. 廊坊: 黄金地质研究所
|
图 12 闪长玢岩Rb-(Y +Nb)(a)和Ta-Yb(b)判别图解(据Pearce et al., 1984) Fig. 12 Rb-(Y+Nb)(a)and Ta-Yb(b)discrimination diagram for diorite porphyry(after Pearce et al., 1984) |
综上所述，金厂铜金矿区在太平洋板块的俯冲作用下，由大洋地壳和大洋远洋沉积物组成的俯冲板片持续向我国东北部的欧亚板块俯冲，在俯冲带的高温高压下，俯冲板片中的流体被挤出，被挤出的流体交代的地幔楔熔融形成深部岩浆，板片的继续俯冲造成岩石圈拆沉或板片折返造成深部岩浆上升侵位，在其演化侵位过程中持续有壳源物质的加入，使闪长玢岩具有壳幔混源的特点。
5.5 成岩时代及地质意义区域上小兴安岭甘河组火山岩(Zhang et al., 2010b)以及吉黑东部绥芬河组、屯田营组和金沟岭组火山岩(纪伟强，2007)都显示与金厂金矿区闪长玢岩相似的微量元素原始地幔标准化特征，且不管是火山岩还是侵入岩其成岩年代都集中于120~110Ma。说明这种受俯冲流体交代的地幔岩浆演化而来的岩浆岩在东北东部可能具有普遍性，其形成可以用中生代晚期太平洋板块的俯冲以及后期的岩石圈剪切拆沉解释(Wang et al., 2006；Zhang et al., 2010a；Wu et al., 2011)。
金厂金矿的成岩及成矿时代前人(陈锦荣等，2000①；张华锋，2007；门兰静，2008；鲁颖淮等，2009)已做过较为详细的研究，综合前人所述认为：①成岩时代集中于印支晚期190~210Ma和燕山中晚期100~130Ma；②金厂金矿的成矿时代集中于燕山中晚期110~125Ma之间。
① 陈锦荣，金宝义，王科强，武玉海. 2000. 黑龙江省东宁县金厂矿区及外围金矿成矿规律与深部预测. 廊坊: 武警黄金地质研究所，1-124
芮宗瑶等(1995) 以小西南岔、闹枝沟、五凤和刺猬沟四个典型矿床为例，建立了吉黑东部斑岩型-浅成热液型铜金矿床多重成矿模型；赵玉锁等(2011) 以延边-东宁成矿带典型铜金矿床建立了“金厂式”——浅成低温-斑岩型金矿成矿模式。笔者认为金厂金矿床并不具有典型斑岩型矿床的蚀变分带特征，通过综合分析金厂金矿床产出大地构造环境、矿床地质特征、矿床地球化学特征、成矿流体(王永，2006；秦江艳2008；门兰静，2008；肖力和赵玉锁，2009②；徐文喜，2009；王可勇等，2010)特征结合成矿年龄，认为金厂铜金矿床与典型斑岩型矿床即有相似性又存在明显的区别，故认为其属于浅成低温-斑岩型矿床的过渡的类型，属斑岩成矿系统的组成部分。
② 肖力，赵玉锁.2009.黑龙江省东宁县金厂矿区及外围金矿构造控矿规律与找矿预测研究. 廊坊: 黄金地质研究所
J-0号矿体的近矿强蚀变闪长玢岩成岩时代为118.1±1.6Ma，岩体中心未蚀变闪长玢岩的成岩时代为115.7±2.0Ma。蚀变闪长玢岩和未蚀变闪长岩的成岩时代存在不一致性，有2.4Ma的差距，本文认为：①锆石样品09JB401与09JB400的测试单位及测试仪器不同，可能为仪器的相对误差引起；②09JB401的误差为1.6Ma, 09JB400的误差为2.0Ma，与2.4Ma差距不大，可以认为其在误差允许范围之内；③蚀变闪长玢岩年龄较老，未蚀变闪长玢岩年龄相对较新，说明岩体边部结晶较早，而岩体中心部位结晶较晚，由于岩浆边部与早期围岩接触温度较低结晶较早，而中心岩浆温度较高，结晶较晚，这与岩体的地质特征反映的情况一致。故本文认为118.1±1.6Ma可以代表J-0号斑岩型矿体的成矿时代，其成矿时代与区域成矿时代有很好的吻合性，其成岩成矿时代与我国东北地区的燕山晚期太平洋板块的俯冲活动相对应。
6 结论黑龙江省东宁金厂铜金矿床为斑岩-浅成低温热液过渡类型矿床，成矿与闪长玢岩关系密切，有时闪长玢岩本身就是角砾岩型矿体(如J-0号矿体)的容矿岩体。J-0号矿体为金厂铜金矿床中较为典型的斑岩型矿体，其深部闪长玢岩岩体的发现及斑岩型矿体的厘定一方面为“金厂式”铜金矿矿床成矿模式提供了佐证，另一方面说明其它浅成低温型矿体的深部可能有存在斑岩型矿体，因此，金厂矿区的找矿潜力非常巨大。
岩相学和地球化学特征表明闪长玢岩岩浆演化过程中经历了结晶分异过程。同时，原始地幔标准化的微量元素蛛网图显示矿区闪长玢岩据富Pb、亏损Nb、Ta和Ti的特征，表明闪长玢岩岩浆可能起源于俯冲流体交代的地幔楔部分熔融，且岩浆上涌过程中经历了地壳混染。此外，闪长玢岩成岩年龄介于118.1±1.6~115.7±2.0Ma，代表金厂铜金矿床成矿的年龄为118.1±1.6Ma，与区域吉黑东部120~110Ma的岩浆活动时间一致，表明成矿作用与太平洋板块俯冲作用以及其后的岩石圈拆沉构造背景有关。
金厂铜金矿床闪长玢岩微量元素丰度呈低Y、高Sr且Sr/Y值大于40的特征，结合其它岩石化学特征显示其为埃达克质岩石；埃达克岩与大型-超大型斑岩型铜金和浅成低温热液型金矿关系密切。因此，埃达克质成矿闪长玢岩的确定可能暗示金厂铜金矿床具备成超大型矿床的条件，深部及外围找矿潜力巨大。
致谢 笔者在野外工作中得到了武警黄金一支队王艳忠总工程师、张景海工程师、程军高级工程师和曹勇助理工程师的大力帮助，在此表示衷心感谢！感谢审稿人对本文的评论和提出的修改意见。| [] | Andersen T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb. Chemical Geology , 192 :59–79. DOI:10.1016/S0009-2541(02)00195-X |
| [] | Ayers J. 1998. Trace element modeling of aqueous fluid-periotite interaction in the mantle wedge of subduction zones. Contributions to Mineralogy and Petrology , 132 :390–404. DOI:10.1007/s004100050431 |
| [] | Bellon H and Yumul Jt GP. 2001. Miocene to Quaternary adakites and related rocks in western Philippine arc sequences. Earth and Planetary Sciences Letters , 333 :343–350. |
| [] | Defant MJ and Drummond MS. 1990. Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature , 347 :662–665. DOI:10.1038/347662a0 |
| [] | Defant MJ, Xu JF, Kepezhinskas P, Wang Q, Zhang Q and Xiao L. 2002. Adakites: Some variations on a theme. Acta Petrologica Sinica , 18 (2) :129–142. |
| [] | Gao LX. 2011. Dynamic background and space-time characteristics analysis of earthquake activities in Northeast China. Earthquake , 31 (1) :41–51. |
| [] | Gao S, Zhang JF, Xu WL and Liu YS. 2009. Delimination and destruction of the North China Craton. Chinese Science Bulletin , 54 (1) :1962–1973. |
| [] | Hou ZQ, Mo XX, Gao YF, Qu XM and Meng XJ. 2003. Adakite, a possible host rock for porphyry copper deposits: Case studies of porphyry copper belts in Tibetan Plateau and in Northern Chile. Mineral Deposits , 22 (1) :1–11. |
| [] | Hou ZQ, Meng XJ and Qu XM. 2005. Copper ore potential of adakitic intrusives in Gangdese porphyry copper belt: Constrains from rock phase and deep melting process. Mineral Deposits , 24 (2) :108–121. |
| [] | Huang JL. 2010. Multi-scale seismic tomography and its tectonic significance of the mainland and key areas in China. The 26th annual meeting of the Chinese Geophysical Society and the 13th Academician Conference of the Seismological Society of China Proceedings,7-8(in Chinese) |
| [] | JahnBM, WuF, LoCH and Tsai CH.1999. Crust-mantle interaction induced by deep subduction of the continental crust: Geochemical and Sr-Nd isotopic evidence from post-colllisional mafic-ultramafic intrusions of the northern Dabie complex. Chemical Geology , 157 :119–146. |
| [] | Ji WQ. 2007. Chronology and geochemistry of Late Mesozoic volcanic rocks in eastern Jilin and Heilongjiang provinces. Ph. D. Dissertation. Changchun: Jilin University, 1-64(in Chinese with English summary) |
| [] | Kay RW. 1978. Aleutian magnesium andesites: Melts from subducted Pacific oceanic crust. J. Volcanol. Geotherm. Res. , 4 :117–132. DOI:10.1016/0377-0273(78)90032-X |
| [] | Li Y, Ling MX, Ding X, Liu J, Han F and Sun WD. 2009. Adakites or adakitic rocks and associated metallogenesis in eastern China. Geotectonica et Metallogenia , 122 (3) :448–464. |
| [] | Lu YH, Zhang Y, Lai Y and Wang YZ. 2009. LA-ICPMS zircon U-Pb dating of magmatism and mineralization in the Jinchang gold ore-field, Heilongjiang Province. Acta Petrologica Sinica , 25 (11) :2902–2912. |
| [] | MacDonald R, Hawkesworth CJ and Heath E. 2000. The lesser Antilles volcanic chain: A study in arc magmatism. Earth Science Reviews , 49 :1–76. DOI:10.1016/S0012-8252(99)00069-0 |
| [] | McDonough WF and Sun SS. 1995. The composition of the Earth. Chemical Geology , 120 :223–253. DOI:10.1016/0009-2541(94)00140-4 |
| [] | Men LJ. 2008. Geological, geochemical characteristics and metallogenic model of Jinchang super-large gold deposit, Dongning, Heilongjiang Province. Master Degree Thesis. Changchun: Jilin University, 1-62(in Chinese with English summary) |
| [] | Middlemost EAK.1985. Magmas and Magmatic Rocks. London: Longman : 1 -266. |
| [] | Middlemost EAK. 1994. Naming materials in the magma/igneous rock system. Earth-Science Reviews , 37 :215–224. DOI:10.1016/0012-8252(94)90029-9 |
| [] | Niu YL. 2010. Some basic concepts and problems on the petrogenesis of intra-plate ocean island basalts. Chinese Science Bulletin , 55 (2) :103–114. |
| [] | Oyarzun R, Morquez A, Lillo J, et al. 2001. Giant versus small porphyry copper deposits of Cenozoic age in northern Chile: Adakitie versus normal calc-alkaline magmatism. Mineral Deposit , 36 :794–798. DOI:10.1007/s001260100205 |
| [] | Pearce JA, Harris NBW and Tindle AG. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. J. Petrol. , 25 :956–983. DOI:10.1093/petrology/25.4.956 |
| [] | Peccerillo R and Taylor SR. 1976. Geochemistry of eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey. Contrib. Mineral. Petrol. , 58 :63–81. DOI:10.1007/BF00384745 |
| [] | Qian Q. 2001. Adakite: Geochemical characteristics and genesis. Acta Petrologica et Mineralogica , 20 (3) :297–306. |
| [] | Qin JY. 2008. Fluids geochemical characters of the magma domal structure gold deposit in Jinchang, Heilongjiang. Master Degree Thesis. Beijing: China University of Geoscience, 1-75(in Chinese with English summary) |
| [] | Rollison HR.1993. Using Geochemical Data: Evaluation, Presentation, Interpretation. New York: Longman Group UK Ltd. : 1 -352. |
| [] | Rudnick RL and Gao S. 2010. Composition of continental crust. In: Holland HD and Turekian KK (eds.). Reading from the Treatise On Geochemistry, 6: 131-195 |
| [] | Rui ZY, Zhang HT, Wang LS, Chen RY, Jin BL, Jin FZ, Wan YS, Zhou YC and Meng QL. 1995. Porphyry-epithermal copper-gold depositsin Yanbian area, Jilin Province. Mineral Deposits , 14 (2) :99–126. |
| [] | Sajona FG and Maury RC. 1998. Association of adakites with gold and copper mineralization in the Philippines. CRACAD SCI Ⅱ A , 326 (1) :27–34. |
| [] | Song SG, Niu YL, Wei CJ, Ji JQ and Su L. 2010. Metamorphism, anatexis, zircon ages and tectonic evolution of the Gongshan block in the northern Indochina continent: An eastern extension of the Lhasa Block. Lithos , 120 :327–34. DOI:10.1016/j.lithos.2010.08.021 |
| [] | Sun SS and McDonough WF. 1989. Chemical and isotope systematicas of oceanic basalts: Implications for mantle composition and processes. In: Saunders AD and Norry MJ (eds.). Magmatism in Oceanic Basins. Spec. Publ. Geol. Soc. Lond., 42: 313-345 |
| [] | Thieblemont D and Stein G. 1997. Lescuyer J-L. Gisements epithemaux et Prophyuiques: la connexion adakite. Earth Planet. Sci. Lett. , 325 :103–109. |
| [] | Wang F, Zhou XH, Zhang LC, Ying JF, Zhang YT, Wu FY and Zhu RX. 2006. Late Mesozoic volcanism in the Great Xing'an Range (NE China): Timing and implications for the dynamic setting of NE Asia. Earth and Planetary Science Letters , 251 :179–198. DOI:10.1016/j.epsl.2006.09.007 |
| [] | Wang KY, Qing M, Zhang XN, Wan D and Xiao L. 2011. Study on the characteristics of fluid inclusions and metallogenic evolution of Jinchang gold deposit, Heilongjiang Province. Acta Petrologica Sinica , 27 (5) :1275–1286. |
| [] | Wang Y. 2006. Geochemical characters of fluid inclusions of the ore-body in the magma domal structure and the argument on the genesis of the deposit in Jinchang gold deposit, Heilongjiang. Master Degree Thesis. Beijing: China University of Geoscience, 1-58 (in Chinese with English summary) |
| [] | Wright JB. 1969. A simple alkalinity ratio and its application to questions of non-orogenic granite genesis. Geological Magazine , 106 :370–384. DOI:10.1017/S0016756800058222 |
| [] | Wu FY, Sun DY, Ge WC, Zhang YB, Grant ML, Wilde SA and Jahn BM. 2011. Geochronology of the Phanerozoic granitoids in northeastern China. Journal of Asian Earth Sciences , 41 :1–30. DOI:10.1016/j.jseaes.2010.11.014 |
| [] | Xu WX. 2009. Geological characters, metallogenic regularity and model of the gold (copper) ore fields in Jinchang, Heilongjiang Province. Ph. D. Dissertation. Beijing: China University of Geoscience, 1-182(in Chinese with English summary) |
| [] | Zhai MG, Meng QR, Liu JM, Hou QL, Hu SB, Li Z, Zhang HF, Liu W, Shao JA and Zhu RX. 2004. Geological features of Mesozoic tectonic regime inversion in Eastern North China and implication for geodynamics. Earth Science Frontiers , 11 (3) :285–297. |
| [] | Zhang HF. 2007. Wall-rock alteration, mineralization time and deposite type of the Jinchang gold deposit from Dongning, Heilongjiang Province, China. Post-Doctor Research Report. Beijing: China University of Geoscience, 1-84(in Chinese with English summary) |
| [] | Zhang JF, Li ZT and Jin CZ. 2004. Adakites in northeastern China and their mineralized implications. Acta Petrologica Sinica , 20 (2) :361–368. |
| [] | Zhang JH, Gao S, Ge WC, Wu FY, Yang JH, Wilde SA and Li M. 2010a. Geochronology of the Mesozoic volcanic rocks in the Great Xing'an Range, northeastern china: Implications for subduction-induced delamination. Chemical Geology , 276 :144–165. DOI:10.1016/j.chemgeo.2010.05.013 |
| [] | Zhang Q, Wang Y, Qian Q, Yang JH, Wang YL, Zhao TP and Guo GJ. 2001. The characteristics and tectonic-metallogenic significances of the adakites in Yanshan Period from eastern China. Acta Petrologica Sinica , 17 (2) :236–244. |
| [] | Zhang Q, Wang Y, Liu W and Wang YL. 2002. Adakite: Its characteristics and implications. Geological Bulletin of China , 21 (7) :431–435. |
| [] | Zhang Q, Xu JF, Wang Y, Xiao L, Liu HT and Wang YL. 2004. Diversity of adakite. Geological Bulletin of China , 23 (9-10) :959–965. |
| [] | Zhang Q, Jin WJ, Xiong XL, Li CD and Wang YL. 2009. Characteristics and implication of O-type adakite in China during different geological periods. Geotectonica et Metallogenia , 33 (3) :432–447. |
| [] | Zhang ZC, Mao JW, Wang YB, Pirajno F, Liu JL and Zhao ZD. 2010b. Geochemistry and geochronology of the volcanic rocks associated with the Dong'an adularia-sericite epithermal gold deposit, Lesser Hinggan Range, Heilongjiang Province, NE China: Constraints on the metallogenesis. Ore Geology Reviews , 37 :158–174. DOI:10.1016/j.oregeorev.2010.03.001 |
| [] | Zhao DP, Yu S and Ohtani E. 2011. East Asia: Seismotectonics, magmatism and mantle dynamics. Journal of Asian Earth Sciences , 40 :689–709. DOI:10.1016/j.jseaes.2010.11.013 |
| [] | Defant MJ, 许 继峰, Kepezhinsksa P, 王 强, 张 旗, 肖 龙. 2002. 埃达克岩:关于其成因的一些不同观点. 岩石学报 , 18 (2) :129–142. |
| [] | 高立新.2011. 中国东北地区地震活动的动力背景及其时空特征分析. 地震 , 31 (1) :41–51. |
| [] | 高山, 章军锋, 许文良, 刘勇胜.2009. 拆沉作用与华北克拉通破坏. 科学通报 , 54 (1) :1962–1973. |
| [] | 郭晓东, 侯增谦, 陈祥, 王治华.2009. 云南省马厂箐富碱斑岩埃达克岩性质的厘定及其成矿意义. 岩石矿物学杂志 , 28 (4) :375–386. |
| [] | 纪伟强. 2007. 吉黑东部中生代晚期火山岩的年代学和地球化学. 博士学位论文. 长春: 吉林大学, 1-64 |
| [] | 黄金莉. 2010. 中国大陆及重点地区多尺度地震层析成像及其构造意义.中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会会议论文,7-8 |
| [] | 侯增谦, 莫宣学, 高永丰, 曲晓明, 孟祥金.2003. 埃达克岩: 斑岩铜矿的一种可能的重要含矿母岩——以西藏和智利斑岩铜矿为例. 矿床地质 , 22 (1) :1–11. |
| [] | 侯增谦, 孟祥金, 曲晓明, 高永丰.2005. 西藏冈底斯斑岩铜矿带埃达克质斑岩含矿性: 源岩相变及深部过程约束. 矿床地质 , 24 :108–121. |
| [] | 李印, 凌明星, 丁兴, 刘健, 韩峰, 孙卫东.2009. 中国东部埃达克岩及成矿作用. 大地构造与成矿学 , 122 (3) :448–464. |
| [] | 路远发.2004. GeoKit:一个用VBA构建的地球化学工具软件包. 地球化学 , 33 (5) :459–464. |
| [] | 鲁颖淮, 张宇, 赖勇, 王艳忠.2009. 黑龙江金厂金矿田岩浆和成矿作用的LA-ICPMS锆石定年. 岩石学报 , 25 (11) :2902–2912. |
| [] | 门兰静. 2008. 黑龙江东宁县金厂超大型金矿床的地质地球化学特征及成矿模式. 硕士学位论文.长春: 吉林大学, 1-62 |
| [] | 牛耀龄.2010. 板内洋岛玄武岩(OIB)成因的一些基本概念和存在的问题. 科学通报 , 55 (2) :103–114. |
| [] | 钱青.2001. Adakite的地球化学特征及成因. 岩石矿物学杂志 , 20 (3) :297–306. |
| [] | 秦江艳. 2008. 黑龙江金厂岩体官窿构造型金矿床的流体地球化学研究. 硕士学位论文. 北京: 中国地质大学, 1-75 |
| [] | 芮宗瑶, 张洪涛, 王龙生, 陈仁义, 金柏禄, 金逢洙, 万玉胜, 周永昶, 孟庆丽.1995. 吉林延边地区斑岩型-浅成热液型金铜矿床. 矿床地质 , 14 (2) :99–126. |
| [] | 王可勇, 卿敏, 张欣娜, 万多, 肖力.2011. 黑龙江金厂金矿床流体包裹体特征及成矿作用研究. 岩石学报 , 27 (5) :1275–1286. |
| [] | 王永. 2006. 黑龙江金厂金矿岩浆穹隆内矿体流体地球化学特征及矿床成因探讨.硕士学位论文.北京: 中国地质大学, 1-58 |
| [] | 徐文喜. 2009. 黑龙江金厂金(铜)矿田地质特征、成矿规律与成矿模式. 博士学位论文. 北京: 中国地质大学, 1-182 |
| [] | 翟明国, 孟庆任, 刘建明, 侯泉林, 胡圣标, 李忠, 张宏福, 刘伟, 邵济安, 朱日祥.2004. 华北东部中生代构造体制转折峰期的主要地质效应和形成动力学探讨. 地学前缘 , 11 (3) :285–297. |
| [] | 张华峰. 2007. 黑龙江省东宁县金厂金矿的围岩蚀变特征与成矿时代、类型研究. 博士后研究工作报告.北京: 中国地质大学, 1-84 |
| [] | 张炯飞, 李之彤, 金成洙.2004. 中国东北部地区埃达克岩及其成矿意义. 岩石学报 , 20 (2) :361–68. |
| [] | 张旗, 王焰, 钱青, 杨进辉, 王元龙, 赵太平, 郭光军.2001. 中国东部燕山期埃达克岩的特征及其构造-成矿意义. 岩石学报 , 17 (2) :236–244. |
| [] | 张旗, 王焰, 刘伟, 王元龙.2002. 埃达克岩的特征及其意义. 地质通报 , 21 (7) :431–435. |
| [] | 张旗, 许继峰, 王焰, 肖龙, 刘红涛, 王元龙.2004. 埃达克岩的多样性. 地质通报 , 23 (9-10) :959–965. |
| [] | 张旗, 金惟俊, 熊小林, 李承东, 王元龙.2009. 中国不同时代O型埃达克岩的特征及其意义. 大地构造与成矿学 , 33 (3) :432–447. |