岩石学报  2017, Vol. 33 Issue (3): 896-906   PDF    
赣杭构造带灵山花岗岩体黑云母的矿物化学特征及其对岩石成因的指示意义
叶茂1, 赵赫2, 赵沔1, 舒珣1, 张若曦1, 杨水源1     
1. 中国地质大学地质过程与矿产资源国家重点实验室, 武汉 430074;
2. 成都理工大学地球科学学院, 成都 610059
摘要: 灵山花岗岩体在平面上为一环状分布的侵入体,中心为角闪石黑云母花岗岩,外围为黑云母花岗岩。在角闪石黑云母花岗岩中分布有大量的暗色镁铁质微粒包体。黑云母是大多数中酸性火成岩中比较重要的一种镁铁质矿物,它能很好地反映寄主岩浆的属性和成岩时的物理、化学条件,因此,本文对这两种花岗岩及镁铁质微粒包体中的黑云母开展了系统的岩相学观察和电子探针化学组成研究,探讨灵山岩体的物质来源、成岩条件和岩浆的混合作用过程。研究结果表明两种花岗岩体的黑云母具有不同化学成分,而暗色镁铁质微粒包体中黑云母的化学成分则变化较大。三种黑云母均在低氧逸度条件下晶出。两种花岗岩中的黑云母均富Fe贫Mg,属于铁质黑云母,含铁系数[(Fe3++Fe2+)/(Fe3++Fe2++Mg2+)]分别为0.65~0.70,0.72~0.78,FeOT/MgO均接近7.04。MF值[2×Mg/(Fe2++Mg+Mn)]分别为0.64~0.76和0.48~0.60,指示两种花岗岩的物质来源都是以壳源为主。镁铁质微粒包体中黑云母的MF值变化范围比较大,为0.63~1.06,为铁质黑云母到镁质黑云母,暗示包体岩浆经历过不同程度的岩浆混合作用。镁铁质微粒包体中部分黑云母与角闪石黑云母花岗岩中黑云母的结晶条件相似,而部分则有明显差异,推测是由于基性的镁铁质包体岩浆注入到酸性的花岗岩浆是一个连续多阶段的过程。
关键词: 黑云母     矿物化学     结晶条件     岩浆混合作用     灵山花岗岩    
Mineral chemistry of biotite and its petrogenesis implication in Lingshan granite pluton, Gan-Hang Belt, SE China
YE Mao1, ZHAO He2, ZHAO Mian1, SHU Xun1, ZHANG RuoXi1, YANG ShuiYuan1     
1. State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan 430074, China;
2. School of Earth Science, Chengdu University of Technology, Chengdu 610059, China
Abstract: The Lingshan composite granite pluton is composed of an outer zone of biotite granite, and a central part of amphibole-biotite granite which contains many mafic microgranular enclaves (MME). Biotite is an important ferromagnesian mineral in most of intermediate-felsic igneous rocks, and its mineral chemistry data can record the properties of magma and the petrogenetic physicochemical conditions. In this study, we carried out a detailed petrographic study together with electric micro-probe analysis on biotites for the biotite granite, the amphibole-biotite granite and the MME, to discuss the source, physicochemical conditions, and the magma mingling/mixing processes in the Lingshan composite granite pluton. The results show significantly different chemical compositions for biotites from the biotite granite and the amphibole-biotite granite, while chemical compositions of biotites in the MME are variable. The crystallization of these biotites was all in a low oxygen fugacity. The biotites in amphibole-biotite granite and biotite granite are both rich in Fe, poor in Mg, and belong to ferri-biotite. The ratio of (Fe3++Fe2+)/(Fe3++Fe2++Mg2+) are 0.65~0.70 and 0.72~0.78, respectively. The ratios of FeOT/MgO of biotites in both granites are close to 7.04. The MF values[2×Mg/(Fe2++Mg+Mn)] are 0.64 to 0.76 for biotite from the amphibole-biotite granite, 0.48 to 0.60 for the biotite from the biotite granite, suggesting these tow granites are both derived from crustal source rocks. The biotites in MME vary from ferri-biotite to magnesio-biotite, and their MF values range between 0.63 and 1.06, which imply that magma of MME experienced magma mixing/mingling in various degrees. Amphibole-biotite granite and parts of MME have a similar crystallization condition, while other MME are obviously different with amphibole-biotite granite. The contrast between amphibole-biotite granite and MME may be a consequence of continuous and multistage interaction between felsic and mafic magmas.
Key words: Biotite     Mineral chemistry     Crystallization condition     Magma mixing     Lingshan pluton    

灵山花岗岩体由中心部分的角闪石黑云母花岗岩和外围的黑云母花岗岩及少部分岩瘤状钠长石黑云母花岗岩和花岗斑岩四种岩性组成。前两种花岗岩构成了灵山岩体的主体部分,且中心部分的角闪石黑云母花岗岩中普遍存在有暗色的镁铁质微粒包体。灵山花岗岩体因产出有Nb-Ta等相关的多金属矿床而受到广泛关注,前人对灵山花岗岩体的研究也多侧重于成矿方面的研究 (陈志中, 1984; 章崇真等, 1985; 袁忠信等, 1987, 1988; 黄定堂, 2003; 章平和田邦生, 2005; 田邦生和周贤旭, 2008; 刘友华等, 2011; Zhu et al., 2015)。对灵山岩体中各类花岗岩的物质来源、岩浆结晶的物理化学条件及演化历史等研究则相对较少。

在很多花岗岩体中常有暗色包体的出现,被广泛地认为指示了两种或以上的不同性质的岩浆在成岩时共存且有混合作用的存在。这些暗色包体为镁铁质岩浆在长英质岩浆演化中起重要作用提供了证据,岩浆的混合作用 (Mingling/Mixing) 研究已经成为了国内外花岗岩研究的热点问题 (Perugini et al., 2003; Barbarin, 2005; Słaby and Martin, 2007; Yang and Jiang, 2013; Temizel et al., 2014; Weidendorfer et al., 2014; 张旗等, 2007)。前人对灵山花岗岩中的镁铁质微粒包体进行过一定程度的研究,且存在不同的成因认识。例如,袁忠信等 (1988)认为灵山花岗岩中的镁铁质微粒包体是寄主角闪石黑云母花岗岩母源岩石经过部分熔融作用后的残留物,而并非岩浆混合作用的产物。郑建平等 (1996)通过镁铁质微粒包体的岩相学和岩石地球化学探讨了包体的成因及其在造山带岩石圈均衡调整中的壳幔相互作用,认为灵山花岗岩中普遍存在岩浆混合作用。Zhou et al. (2013)在研究江南造山带中的A型花岗岩时,对灵山花岗岩的中心和外围两种花岗岩展开了岩石地球化学研究工作,包括锆石U-Pb定年,主量、微量、Sr-Nd-Hf同位素组成的研究,认为灵山花岗岩为A型花岗岩,并且灵山花岗岩具有相对较高的全岩εNd(t) 值 (外围的黑云母花岗岩εNd(t) 值为-0.96,中心部分的角闪石黑云母花岗岩εNd(t) 值为-3.28) 和锆石εNd(t) 值 (外围的黑云母花岗岩εNd(t) 平均值为-0.47,中心部分的角闪石黑云母花岗岩εNd(t) 平均值为-2.02),指示了这些花岗岩形成时可能有幔源组分的加入,存在壳幔混合作用。

岩浆岩中的矿物地球化学数据可以为岩浆固结的物理化学条件、岩浆房中发生的各种地质过程 (如分离结晶、岩浆的混合作用和岩浆的消耗与补充等演化过程) 提供直接的证据。矿物的化学成分可以反映出矿物结晶的物理化学参数 (压力、温度和氧逸度等)、结晶生长历史和岩浆物质来源等信息 (Speer, 1987; Temizel et al., 2014)。黑云母是花岗岩中最常见的暗色矿物之一,其矿物化学成分中各元素存在广泛的类质同像替换,如Fe2+-Mg-Mn、Si-Al-Fe-Ti及OH-F-Cl等 (Rieder et al., 1999; 孙世华和于洁, 1989; 吕志成等, 2003)。黑云母化学成分很好地记录了其寄主岩浆的属性,它的成分变化是寄主岩浆固结成岩时不同的地质环境和物理化学条件下所对应的结果,近年来已经被应用到岩浆混合作用的研究领域 (Pignatelli et al., 2016; Gao et al., 2016)。通过对黑云母矿物化学的研究可以反演出寄主岩体的成岩条件、演化过程及岩浆物质来源等。因此,对黑云母的研究具有重要的地质意义并已广泛的应用于岩浆岩成因的研究之中 (Burkhard, 1991; El Sheshtawi et al., 1993; Abdel-Rahman, 1994, 1996; Stone, 2000; Zhao et al., 2005; Zhu et al., 2014; Xu et al., 2015; 彭花明, 1997; 张遵忠等, 2005; 蒋少涌等, 2006, 2008; 东前等, 2011; 徐耀明等, 2013; 赵沔等, 2015)。

因此,为了进一步研究灵山花岗岩及镁铁质微粒包体的成因,特别是为了揭示该岩体形成过程中是否存在岩浆混合作用,本文拟以灵山岩体中黑云母为研究对象,开展详细的岩相学研究和矿物化学特征研究,通过对比灵山岩体外围和中心的花岗岩及镁铁质微粒包体中的黑云母,结合前人已有的岩石地球化学数据 (锆石U-Pb定年,主量、微量、Sr-Nd-Hf同位素组成的研究),探讨灵山花岗岩的岩石成因及其与镁铁质包体岩浆混合作用的过程。

1 地质概况

灵山岩体位于江西省的东北角,大地构造上处于两个大地一级构造单元的结合部位,以宜春-绍兴断裂为界,北区为扬子地块,南区为华夏地块。区域上处于赣杭构造带上,距德兴地体与怀玉地体间的赣东北深断裂仅约10km (章平和田邦生, 2005)。灵山岩体出露面积约200km2,平面呈一近似圆形的环状体,如图 1所示,主要由四个期次的侵入体组成。中心位置为第一期次的角闪石黑云母花岗岩 (约65%),构成灵山岩体的主体部分。第二期次的黑云母花岗岩 (约30%) 分布在角闪石黑云母花岗岩的外围。两者被第三、第四期的小型细粒钠长石化黑云母花岗岩和花岗斑岩在四周零星地侵入穿切。

图 1 灵山复式岩体地质简图 (据Zhu et al., 2015修改) Fig. 1 Geological sketch map of Lingshan granitic pluton (modified after Zhu et al., 2015)
2 样品描述

野外采样沿着茗洋关水库东岸,从岩体外围到中心地区 (图 1),依次采取新鲜的黑云母花岗岩 (图 2a)、角闪石黑云母花岗岩 (图 2b) 和角闪石黑云母花岗岩中不同色率的暗色包体样品 (图 2c)。

图 2 灵山黑云母花岗岩 (a、d)、角闪石黑云母花岗岩 (b、e) 和镁铁质微粒包体 (MME)(c、f) 的手标本照片及单偏光镜下照片 Bt-黑云母;Amp-角闪石;Pl-斜长石;Kfs-钾长石;Qt-石英 Fig. 2 The hand specimen photographs and microphotographs of biotite granite (a, d), amphibole-biotite granite (b, e) and mafic microgranular enclaves (c, f) from Lingshan Bt-biotite; Amp-amphibole; Pl-plagioclase; Kfs-K-feldspar; Qt-quartz

黑云母花岗岩为灵山岩体的环状外围部分,岩石呈肉红色-黄褐色,中至粗粒,块状构造。矿物组成为钾长石 (40%~60%),斜长石 (15%~25%),石英 (25%~35%) 及黑云母 (4%~8%)。黑云母呈半自形生长于其他颗粒之间 (图 2d),但是并没有发生挤压变形。黑云母中常包含锆石、磷灰石、磁铁矿等副矿物。角闪石黑云母花岗岩为灵山岩体的中心主体部分,岩石呈灰白色,粗粒至巨粒,似斑状结构,块状构造。主要矿物组分为钾长石 (30%~50%),斜长石 (25%~35%),石英 (25%~30%),角闪石 (2%~5%) 及黑云母 (4%~8%),副矿物见有磁铁矿、榍石、锆石、磷灰石、钛铁矿等。黑云母呈半自形生长于其他颗粒之间 (图 2e),并没有发生挤压变形。在灵山岩体中,仅在中心部分的角闪石黑云母花岗岩中发现有暗色镁铁质微粒包体存在。镁铁质微粒包体呈卵圆形、椭球形等,大小不一,小者直径3~5cm,大者直径则可达1.5~2m。包体颜色深浅也不一样,根据色率的高低可以形成一个较为完整的连续系列。包体总体呈灰黑色,有时为似斑状结构,有时为细粒-微粒半自形等粒结构,块状构造。主要矿物成分为云母 (15%~30%),角闪石 (10%~25%),斜长石 (20%~30%),钾长石 (5%~15%) 及石英 (5%~15%),副矿物主要是呈长柱状的磷灰石,以及榍石、锆石、磁铁矿、钛铁矿等 (图 2f)。

3 分析方法及结果

电子探针测试工作在中国地质大学 (武汉) 地质过程与矿产资源国家重点实验室电子探针室完成,电子探针型号为JEOL JXA-8100。样品在上机测试之前先按照Zhang and Yang (2016)提供的办法进行了镀碳,将样品镀上尽量均匀的厚度约20nm的碳膜。测试时加速电压为15kV,加速电流20nA,电子束斑为1μm,测试结果采用ZAF方法校正,元素标样全部采用美国SPI提供的矿物标样,分别为橄榄石 (Si)、金红石 (Ti)、镁铝榴石 (Al, Fe)、蔷薇辉石 (Mn)、透辉石 (Mg, Ca)、硬玉 (Na)、透长石 (K) 及萤石 (F)。以22个氧为基础 (Rieder et al., 1999) 计算黑云母的结构式,并采用林文蔚和彭丽君 (1994)的计算方法计算黑云母中的Fe2+、Fe3+含量。计算结果见表 1

表 1 灵山黑云母花岗岩、角闪石黑云母花岗岩及包体中的黑云母的化学组成 (wt%) Table 1 Electron microprobe analysis of the biotite from the biotite granite, amphibole-biotite granite and the mafic microgranular enclaves (wt%)

表 1的电子探针分析结果及图 3黑云母成分变化趋势图来看,两种花岗岩均体现富Fe、Ti、F,贫Mg的特征。角闪石黑云母花岗岩中黑云母的FeOT(全铁,下同) 质量分数为25.83%~28.71%,MgO的质量分数为6.89%~8.09%,含铁系数[(Fe3++Fe2+)/(Fe3++Fe2++Mg)]为0.65~0.70。黑云母花岗岩中的黑云母的FeOT质量分数为28.12%~30.14%,MgO的质量分数为4.87%~6.28%,含铁系数[(Fe3++Fe2+)/(Fe3++Fe2++Mg)]为0.72~0.78。两者相对比,黑云母花岗岩要比角闪石黑云母花岗岩中的黑云母更富Fe,而贫Mg。镁铁质微粒包体的黑云母的化学成分情况与前两者有明显差异,FeOT质量分数为21.63%~28.69%,MgO为6.78%~12.42%,含铁系数[(Fe3++Fe2+)/(Fe3++Fe2++Mg)]为0.49~0.70。与花岗岩中的黑云母相比,Mg、Fe含量变化范围大,一部分镁铁质微粒包体中黑云母的成分与角闪石黑云母花岗岩中的黑云母相似,而一部分则相对富Mg贫Fe。由镁铁质微粒包体→角闪石黑云母花岗岩→黑云母花岗岩,MF[2×Mg/(Fe2++Mg+Mn)]平均值不断减小 (0.75→0.70→0.55)。

图 3 灵山黑云母花岗岩、角闪石黑云母花岗岩及镁铁质微粒包体中黑云母成分变化趋势图 Fig. 3 Composition variations of biotites in biotite granite, amphibole-biotite granite and mafic microgranular enclaves from Lingshan

角闪石黑云母花岗岩、黑云母花岗岩和包体中黑云母的TiO2和Al2O3含量变化不大,三者的TiO2变化范围分别为2.79%~4.05%,3.07%~3.72%,2.03%~3.95%;Al2O3变化范围分别为11.28%~12.32%,11.26%~12.54%,11.24%~13.08%,但总的来说镁铁质微粒包体中的黑云母的Al和Ti含量变化范围相对较大 (图 3),这与包体岩浆与其寄主岩浆的混合程度及黑云母的结晶环境变化有关 (Temizel et al., 2014)。

在Mg-(Al+Fe3++Ti)-(Fe2++Mn) 黑云母分类图 (图 4) 中,角闪石黑云母花岗岩中的黑云母落在黑云母花岗岩中的黑云母的上方,但都落在铁质黑云母区域;镁铁质微粒包体中的黑云母分布范围最广,既有镁质黑云母又有铁质黑云母。

图 4 灵山黑云母花岗岩、角闪石黑云母花岗岩及镁铁质微粒包体中黑云母分类图 (据Foster, 1960) A-金云母;B-镁质黑云母;C-铁质黑云母;D-铁叶云母;E-铁白云母;F-白云母.图例同图 3,后图同 Fig. 4 Classification diagram for biotites in biotite granite, amphibole-biotite granite and mafic microgranular enclaves from Lingshan (after Foster, 1960) A-phlogopite; B-magnesian biotite; C-ferrobiotite; D-sidero-phyllite; E-ferro-muscovite; F-muscovite. See Fig. 3 for legend in Fig. 4-Fig. 7
4 讨论 4.1 岩石成因

岩浆岩中黑云母的成分特征很好地记录了寄主岩浆的物理化学信息 (Abdel-Rahman, 1994, 1996)。Wones and Eugster (1965)利用低Al,且与磁铁矿和钾长石共生的黑云母,实验研究出黑云母中Fe3+、Fe2+及Mg三者的含量可估算出云母结晶时的氧逸度。根据镜下观察和电子探针数据,本区的黑云母符合上述估算方法的条件。在Fe3+-Fe2+-Mg三角图解 (图 5) 中,本区三种黑云母样品落在Ni-NiO缓冲线附近,表明三者都是在较低氧逸度条件下结晶形成的。

图 5 灵山黑云母花岗岩、角闪石黑云母花岗岩及镁铁质微粒包体中黑云母的Fe3+-Fe2+-Mg图解 (据Wones and Eugster, 1965) Fig. 5 Fe3+-Fe2+-Mg ternary diagram for biotites in biotite granite, amphibole-biotite granite and mafic microgranular enclaves from Lingshan (after Wones and Eugster, 1965)

本区的角闪石黑云母花岗岩与黑云母花岗岩中黑云母的MF值的变化范围分别为0.64~0.76、0.48~0.60。根据谢应雯和张玉泉 (1987)对花岗岩类岩石成因类型的总结,本区两种花岗岩中黑云母的MF值在0.5~1.0范围内,均属富铁黑云母,是壳型花岗岩特有的标型矿物。周作侠 (1988)利用FeOT/(FeOT+MgO)-MgO图解 (图 6),来区分三种不同物质来源地质体中的云母。图 6也指示出角闪石黑云母花岗岩和黑云母花岗岩均具壳源型花岗岩的特征。从图 6中还可以看到角闪石黑云母花岗岩接近壳幔混合源与壳源花岗岩的边界,这可能是镁铁质包体岩浆与其发生岩浆混合作用的体现。暗色镁铁质微粒包体的黑云母MF值变化范围为0.63~1.06,为铁黑云母到镁黑云母,谢应雯和张玉泉 (1987)指出镁黑云母 (MF=1~1.3) 是壳幔混合型花岗岩的标型特征,在图 6中也可看到部分镁铁质微粒包体中的黑云母落在了壳幔混合区,可以看作是较接近原始包体岩浆中结晶的黑云母;而其它落在壳源区的部分,则是包体岩浆与花岗岩岩浆混合作用后的岩浆结晶出的黑云母。

图 6 灵山黑云母花岗岩、角闪石黑云母花岗岩及镁铁质微粒包体中黑云母的FeOT/(FeOT+MgO)-MgO (%) 图解 (据周作侠, 1988) Fig. 6 FeOT/(FeOT+MgO)-MgO diagram on material sources of biotites in biotite granite, amphibole-biotite granite and mafic microgranular enclaves from Lingshan (after Zhou, 1988)

Abdel-Rahman (1994)曾对造山和非造山岩系中的黑云母成分有过系统研究,提出了利用黑云母成分判别构造环境的判别图解 (图 7),指出造山钙碱性岩系 (I型花岗岩,C区) 中的黑云母相对富Mg;过铝岩系 (S型花岗岩,P区) 中的黑云母富Al,为铁叶云母质的花岗岩;而非造山碱性岩系 (A型花岗岩,A区) 中的黑云母则相对富Fe,近铁云母。本区角闪石黑云母花岗岩与黑云母花岗岩中的黑云母均表现富Fe,在Mg-(Al+Fe3++Ti)-(Fe2++Mn) 黑云母分类图 (图 4) 中都落在铁质黑云母区域,两者FeOT/MgO均接近非造山碱性岩系的值 (7.04; Abdel-Rahman, 1994),在图 7中也指示两种花岗岩的构造环境为非造山碱性岩系。因此,本区的角闪石黑云母花岗岩和黑云母花岗岩具有A型花岗岩的地球化学特征。Zhou et al. (2013)开展灵山花岗岩的岩石地球化学研究结果显示,灵山两种花岗岩都具有富稀土元素,Nb、Zr、Hf及其它高场强元素,亏损Ba和Sr,且高Ga/Al比值,在Zr+Nb+Ce+Y vs. FeOT/MgO和Zr+Nd+Ce+Y vs. 104×Ga/Al图解 (Whalen et al., 1987) 中位于A型花岗岩区域,表明灵山花岗岩属于A型花岗岩。我们根据黑云母的矿物化学对花岗岩成因类型的判断与根据岩石地球化学特征的判断是相符合的。灵山花岗岩体的形成年龄在130Ma左右 (Zhou et al., 2013),与赣杭构造带中的相山盆地火山侵入杂岩 (Yang et al., 2010, 2011; 杨水源等, 2012, 2013)、大茅山花岗岩 (Jiang et al., 2011)、铜山花岗岩 (Jiang et al., 2011)、白菊花尖花岗岩 (Wong et al., 2009)、新路盆地中的火山侵入杂岩 (Yang et al., 2012; 杨水源等, 2015)、大洲流纹岩 (Yang et al., 2013) 等岩体的形成年龄大致相同,都属于早白垩世岩浆活动的产物。岩石地球化学研究表明这些岩体具有A型花岗岩的地球化学特征,并且,前人认为该时期的赣杭构造带处于太平洋板块与华南板块俯冲碰撞后撤所引起的弧后拉张环境 (Yang et al., 2012; Li et al., 2014)。

图 7 灵山黑云母花岗岩、角闪石黑云母花岗岩及镁铁质微粒包体中黑云母构造环境判别图 (据Abdel-Rahman, 1994) A-A型花岗岩;C-I型花岗岩;P-S型花岗岩 Fig. 7 Discrimination diagram of tectonic setting for biotites in biotite granite, amphibole-biotite granite and mafic microgranular enclaves from Lingshan (after Abdel-Rahman, 1994) A-A-type granite; C-I-type granite; P-S-type granite
4.2 岩浆混合作用

花岗岩中的包体可以是多成因的,可以是围岩的捕掳体,与花岗岩同源的晚期析离体,以及不同性质成分的岩浆不完全混合的产物。对于花岗岩中包体的成因类型一直以来都存在严重分歧,但最新国内外研究表明岩浆的混合成因是最重要和被广泛接受的成因模式 (Perugini et al., 2003; Barbarin, 2005; Słaby and Martin, 2007; Zhao et al., 2010; Temizel et al., 2014; Weidendorfer et al., 2014; 张旗等, 2007)。

灵山花岗岩中的镁铁质微粒包体也属于岩浆混合成因。从镁铁质微粒包体的岩相学特征来看,镁铁质微粒包体与寄主花岗岩之间常存在岩浆快速淬冷形成的反应边,也有表现为一个色率较浅的过渡带,而并没有截然的界限。镁铁质微粒包体具有斑状结构或者微粒半自形等粒结构等典型的岩浆结构。此外,镁铁质微粒包体中发育大量长宽比10~40或更高的针状、长柱状磷灰石 (图 2f)。这些都说明镁铁质微粒包体为岩浆成因,发生过淬冷,与花岗岩岩浆发生过混合作用。

另一方面,前人对镁铁质微粒包体及其寄主花岗岩的岩石化学分析结果列于表 2。镁铁质微粒包体的SiO2含量低,且变化范围较大,在58.53%~66.90%之间,平均值为62.12%,为中性岩,其寄主角闪石黑云母花岗岩的SiO2含量较高,变化较小,在67.75%~69.89%之间,平均值为69.22%。镁铁质微粒包体中Fe、Mg,Ti、Ca、P含量远高于寄主花岗岩,其变化范围也很大。综合上述镁铁质微粒包体及其寄主花岗岩的岩石地球化学特征,我们推测由于镁铁质包体岩浆受到了寄主花岗岩浆不同程度的岩浆混合作用,才导致镁铁质微粒包体的成分属于中性岩。Zhou et al. (2013)的测试结果 (表 2) 显示,灵山内、外两种花岗岩具有较高的εNd(t) 值 (分别为-3.28,-0.96) 和较高的锆石εNd(t) 值 (分别为-3.7~-0.4,-2.2~+0.8),暗示灵山花岗岩可能有地幔物质的加入。Zhou et al. (2013)对灵山花岗岩的锆石U-Pb定年为130Ma左右,但并未对其中的镁铁质微粒包体进行U-Pb定年,但根据镁铁质微粒包体的岩相学特征,灵山的镁铁质微粒包体应是与其寄主花岗岩同时代岩浆作用的产物。综上,我们认为灵山花岗岩中的镁铁质微粒包体是壳幔混合作用直接的证据。从黑云母的矿物化学特征来看,包体中黑云母的Mg含量随着Fe的增加逐渐降低,呈现明显的线性负相关 (图 3e)。包体中的黑云母在Mg-(Al+Fe3++Ti)-(Fe2++Mn) 分类图 (图 4) 中则既有镁质黑云母又有铁质黑云母,其MF值 (表 1) 为0.63~1.06,变化范围大。这些特征表明随着花岗岩岩浆与包体岩浆不同程度的物质交换,混合作用进行的程度有深有浅,故导致镁铁质微粒包体中结晶出来的黑云母成分变化范围比较大。这与手标本中包体的色率呈现出一个由深到浅的系列是对应的。通过对比研究,色率较深的包体中的黑云母具有较高的Mg含量,成分上更接近或位于壳幔混合的范围里面 (图 6)。这些色率较高的包体样品少有钾长石、斜长石和石英的捕掳斑晶,表明这部分包体岩浆与花岗岩浆的混合作用相对较弱。

表 2 灵山镁铁质微粒包体及其寄主花岗岩的岩石化学成分 (wt%) Table 2 Chemical composition for mafic microgranular enclaves and their host rocks from Lingshan (wt%)

据上所述,包体岩浆注入到花岗质岩浆应是一个连续多阶段的过程。在花岗岩岩浆演化初期,基性的镁铁质包体岩浆注入到花岗岩岩浆房中。由于包体岩浆注入的扰动及花岗岩浆的热对流作用,包体岩浆被扯碎为一个个小岩浆团,由于此时花岗岩岩浆的粘度较低,包体岩浆受到花岗岩岩浆的混合作用较强,包体岩浆的成分受花岗岩岩浆成分的影响较大,导致包体中结晶出来的黑云母与寄主花岗岩中的黑云母成分相近,富Fe贫Mg,并而形成色率较浅的包体。而在花岗岩岩浆演化的晚期,由于花岗岩岩浆粘度增大,导致此时注入到花岗岩岩浆中的包体岩浆与花岗岩岩浆的混合作用也逐渐变弱,包体岩浆的成分受花岗岩岩浆的影响越来越小,因此包体中晶出的黑云母就保留了原始相对富Mg的特征。由于混合作用的程度由强到弱,就形成了色率逐渐增高的各种类型的包体。因此灵山花岗岩中的暗色镁铁质微粒包体可能是拉张伸展的动力学环境下镁铁质基性岩浆与壳源长英质酸性岩浆不同程度混合作用的产物。

5 结论

本文通过对灵山角闪石黑云母花岗岩、黑云母花岗岩及角闪石黑云母花岗岩中镁铁质微粒包体三者中黑云母的化学成分分析,探讨了灵山花岗岩的岩石成因和壳幔混合作用的过程,得出以下两点结论:

(1) 角闪石黑云母花岗岩和黑云母花岗岩中的黑云母富Fe贫Mg,均为铁质黑云母,后者更富Fe,两者的其它化学成分差异不大。角闪石黑云母花岗岩和黑云母花岗岩中的黑云母指示两种花岗岩的物质来源以壳源为主。镁铁质微粒包体的黑云母相对富Mg,成分变化较大,为铁质黑云母到镁质黑云母。三种黑云母都是在较低氧逸度的环境下形成的。结合前人的研究认为本区两种花岗岩属非造山碱性岩系,具有A型花岗岩的地球化学特征。

(2) 灵山花岗岩存在壳幔混合作用,角闪石黑云母花岗岩中镁铁质微粒包体的黑云母记录了岩浆混合作用过程,镁铁质的包体岩浆注入到花岗岩岩浆中是一个连续多阶段的过程,随着花岗岩岩浆的演化,镁铁质包体岩浆受到花岗质岩浆混合作用的强度逐渐降低。

致谢 本文的野外采样得到了江西省核工业地质局二六五大队徐成勇工程师和中国地质大学 (武汉) 刘国奇的帮助和支持,在此表示感谢。审稿过程中两位匿名审稿人对本文提出了非常具有建设性的意见,作者谨致谢意!
参考文献
[] Abdel-Rahman AFM. 1994. Nature of biotites from alkaline, caloalkaline and peraluminous magmas. Journal of Petrology, 35(2): 525–541. DOI:10.1093/petrology/35.2.525
[] Abdel-Rahman AFM. 1996. Discussion on the comment on nature of biotites in alkaline, calc-alkaline and peraluminous magmas. Journal of Petrology, 37(5): 1031–1035. DOI:10.1093/petrology/37.5.1031
[] Barbarin B. 2005. Mafic magmatic enclaves and mafic rocks associated with some granitoids of the central Sierra Nevada batholith, California:Nature, origin, and relations with the hosts. Lithos, 80(1-4): 155–177. DOI:10.1016/j.lithos.2004.05.010
[] Burkhard DJM. 1991. Temperature and redox path of biotite-bearing intrusives:A method of estimation applied to S-and I-type granites from Australia. Earth and Planetary Science Letters, 104(1): 89–98. DOI:10.1016/0012-821X(91)90240-I
[] Chen ZZ. 1984. A study on division of stage and their metallogenic relationship with tungsten and tantalum in Linshan, northeastern Jiangxi. Geotectonica et Metallogenia, 8(1): 59–70.
[] Dong Q, Du YS, Cao Y, Pang ZS, Song LX, Zheng Z. 2011. Compositional characteristics of biotites in Wushan granodiorite, Jiangxi Province:Implications for petrogenesis and mineralization. Journal of Mineralogy and Petrology, 31(2): 1–6.
[] El Sheshtawi YA, Salem AKA, Aly MM. 1993. The geochemistry of ferrous biotite and petrogenesis of Wadi-El-Sheikh granitoid rocks Southwestern Sinai, Egypt. Journal of African Earth Sciences (and the Middle East), 16(4): 489–498. DOI:10.1016/0899-5362(93)90106-Z
[] Foster MD. 1960. Interpretation of the composition of trioctahedral micas. Geological Survey Professional Paper, 354B:1-49
[] Gao P, Zhao ZF, Zheng YF. 2016. Magma mixing in granite petrogenesis:Insights from biotite inclusions in quartz and feldspar of Mesozoic granites from South China. Journal of Asian Earth Sciences, 123: 142–161. DOI:10.1016/j.jseaes.2016.04.003
[] Huang DT. 2003. Evolving characteristics and related rare-metal metallogenesis of Lingshan rock body. Geology and Prospecting, 39(4): 35–40.
[] Jiang SY, Zhao KD, Jiang YH, Ling HF, Ni P. 2006. New type of tin mineralization related to granite in South China:Evidence from mineral chemistry, element and isotope geochemistry. Acta Petrologica Sinica, 22(10): 2509–2516.
[] Jiang SY, Li L, Zhu B, Ding X, Jiang YH, Gu LX, Ni P. 2008. Geochemical and Sr-Nd-Hf isotopic compositions of granodiorite from the Wushan copper deposit, Jiangxi Province and their implications for petrogenesis. Acta Petrologica Sinica, 24(8): 1679–1690.
[] Jiang YH, Zhao P, Zhou Q, Liao SY, Jin GD. 2011. Petrogenesis and tectonic implications of Early Cretaceous S-and A-type granites in the northwest of the Gan-Hang rift, SE China. Lithos, 121(1-4): 55–73. DOI:10.1016/j.lithos.2010.10.001
[] Li JH, Zhang YQ, Dong SW, Johnston ST. 2014. Cretaceous tectonic evolution of South China:A preliminary synthesis. Earth-Science Reviews, 134: 98–136. DOI:10.1016/j.earscirev.2014.03.008
[] Lin WW, Peng LJ. 1994. The estimation of Fe3+ and Fe2+ contents in amphibole and biotite from EMPA data. Journal of Changchun University of Earth Sciences, 24(2): 155–162.
[] Liu YH, Li KD, Tu JF. 2011. Geological characteristics and genesis of Huangshan niobium (tantalum) deposit in Geyuan, Jiangxi Province. Resources Survey & Environment, 32(4): 291–298.
[] Lü ZC, Duan GZ, Dong GH. 2003. Mineral chemistry of biotite from granites associated with different mineralization in three stages of Yanshanina Period in the southern-middle parts of the Da Hinggan Ling mountains and its petrogenetic and metallogenic significance. Acta Mineralogica Sinica, 23(2): 177–184.
[] Peng HM. 1997. Geological characteristics of biotite from Yangxi granite body and their geological implications. Acta Petrologica et Mineralogica, 16(3): 271–281.
[] Perugini D, Poli G, Christofides G, Eleftheriadis G. 2003. Magma mixing in the Sithonia Plutonic Complex, Greece:Evidence from mafic microgranular enclaves. Mineralogy and Petrology, 78(3-4): 173–200. DOI:10.1007/s00710-002-0225-0
[] Pignatelli I, Faure F, Mosser-Ruck R. 2016. Self-mixing magma in the Ruiz Peak rhyodacite (New Mexico, USA):A mechanism explaining the formation of long period polytypes of mica. Lithos, 266-267: 332–347. DOI:10.1016/j.lithos.2016.10.024
[] Rieder M, Cavazzini G, D'Yakonov YS, Frank-Kamenetskii VA, Gottardi G, Guggenheim S, Koval PV, Müller G, Neiva AMR, Radoslovich EW, Robert JL, Sassi FP, Takeda H, Weiss Z, Wones DR. 1999. Nomenclature of the micas. Mineralogical Magazine, 63(2): 267–296. DOI:10.1180/minmag
[] Słaby E, Martin H. 2007. Mafic and felsic magma interaction in granites:The Hercynian Karkonosze pluton (Sudetes, Bohemian Massif). Journal of Petrology, 49(2): 353–391. DOI:10.1093/petrology/egm085
[] Speer JA. 1987. Evolution of magmatic AFM mineral assemblages in granitoid rocks:The hornblende+melt=biotite reaction in the Liberty Hill pluton, South Carolina. American Mineralogist, 72: 863–878.
[] Stone D. 2000. Temperature and pressure variations in suites of Archean felsic plutonic rocks, Berens River area, Northwest Superior Province, Ontario, Canada. The Canadian Mineralogist, 38(2): 455–470. DOI:10.2113/gscanmin.38.2.455
[] Sun SH, Yu J. 1989. Interpretation of chemical composition and subdivision of Mg-Fe micas, Part B:The natural subdivision of Mg-Fe micas, Part B:The natural subdivision of Mg-Fe micas. Scientia Geologica Sinica(2): 176–189.
[] Temizel I·, Arslan M, Abdioğlu E, Yücel C. 2014. Mineral chemistry and thermobarometry of Eocene monzogabbroic stocks from the Bafra (Samsun) area in Turkey:Implications for disequilibrium crystallization and emplacement conditions. International Geology Review, 56(10): 1226–1245. DOI:10.1080/00206814.2014.933363
[] Tian BS, Zhou XX. 2008. Gold ore belt division and the metallogenic models in northeast Jiangxi trialangle area. Contributions to Geology and Mineral Resources Research, 23(2): 130–134.
[] Weidendorfer D, Mattsson HB, Ulmer P. 2014. Dynamics of magma mixing in partially crystallized magma chambers:Textural and petrological constraints from the basal complex of the austurhorn intrusion (SE Iceland). Journal of Petrology, 55(9): 1865–1903. DOI:10.1093/petrology/egu044
[] Whalen JB, Currie KL, Chappell BW. 1987. A-type granites:Geochemical characteristics, discrimination and petrogenesis. Contributions to Mineralogy and Petrology, 95(4): 407–419. DOI:10.1007/BF00402202
[] Wones DR, Eugster HP. 1965. Stability of biotite:Experiment, theory, and application. The American Mineralogist, 50: 1228–1272.
[] Wong J, Sun M, Xing GF, Li XH, Zhao GC, Wong K, Yuan C, Xia XP, Li LM, Wu FY. 2009. Geochemical and zircon U-Pb and Hf isotopic study of the Baijuhuajian metaluminous A-type granite:Extension at 125~100Ma and its tectonic significance for South China. Lithos, 112(3-4): 289–305. DOI:10.1016/j.lithos.2009.03.009
[] Xie YW, Zhang YQ. 1987. Typomorphic peculiarities of biotites from different genetic types of granite in the Hengduanshan region. Acta Mineralogica Sinica, 7(3): 245–254.
[] Xu B, Jiang SY, Wang R, Ma L, Zhao KD, Yan X. 2015. Late Cretaceous granites from the giant Dulong Sn-polymetallic ore district in Yunnan Province, South China:Geochronology, geochemistry, mineral chemistry and Nd-Hf isotopic compositions. Lithos, 218-219: 54–72. DOI:10.1016/j.lithos.2015.01.004
[] Xu YM, Jiang SY, Zhu ZY, Zhou W, Kong FB, Sun MZ, Xiong YG. 2013. Geochronology, geochemistry and mineralogy of ore-bearing and ore-barren intermediate-acid intrusive rocks from the Jiurui ore district, Jiangxi Province and their geological implications. Acta Petrologica Sinica, 29(12): 4291–4310.
[] Yang SY, Jiang SY, Jiang YH, Zhao KD, Fan HH. 2010. Zircon U-Pb geochronology, Hf isotopic composition and geological implications of the rhyodacite and rhyodacitic porphyry in the Xiangshan, Jiangxi Province, China. Science China (Earth Sciences), 53(10): 1411–1426. DOI:10.1007/s11430-010-4058-0
[] Yang SY, Jiang SY, Jiang YH, Zhao KD, Fan HH. 2011. Geochemical, zircon U-Pb dating and Sr-Nd-Hf isotopic constraints on the age and petrogenesis of an Early Cretaceous volcanic-intrusive complex at Xiangshan, Southeast China. Mineralogy and Petrology, 101(1-2): 21–48. DOI:10.1007/s00710-010-0136-4
[] Yang SY, Jiang SY, Zhao KD, Jiang YH, Ling HF, Luo L. 2012. Geochronology, geochemistry and tectonic significance of two Early Cretaceous A-type granites in the Gan-Hang Belt, Southeast China. Lithos, 150: 155–170. DOI:10.1016/j.lithos.2012.01.028
[] Yang SY, Jiang SY, Zhao KD, Jiang YH, Fan HH. 2012. Zircon U-Pb geochronology, geochemistry and Sr-Nd-Hf isotopic compositions of the rhyolite porphyry from the Zhoujiaoshan deposit in Xiangshan uranium ore field, Jiangxi Province, SE China. Acta Petrologica Sinica, 28(12): 3915–3928.
[] Yang SY, Jiang SY. 2013. Occurrence and significance of a quartz-amphibole schist xenolith within a mafic microgranular enclave in the Xiangshan volcanic-intrusive complex, SE China. International Geology Review, 55(7): 894–903. DOI:10.1080/00206814.2012.752662
[] Yang SY, Jiang SY, Zhao KD, Jiang YH. 2013. Petrogenesis and tectonic significance of Early Cretaceous high-Zr rhyolite in the Dazhou uranium district, Gan-Hang Belt, Southeast China. Journal of Asian Earth Sciences, 74: 303–315. DOI:10.1016/j.jseaes.2012.12.024
[] Yang SY, Jiang SY, Zhao KD, Jiang YH, Ling HF, Chen PR. 2013. Timing and geological implications of volcanic rocks from the Ruyiting section, Xiangshan uranium ore field, Jiangxi Province, SE China. Acta Petrologica Sinica, 29(12): 4362–4372.
[] Yang SY, Wen X, Zhao KD, Jiang YH, Ling HF, Chen PR. 2015. Geochronology, petrogenesis and geological significance of the trachyte porphyry from Xinlu Basin, Gan-Hang Belt, SE China. Acta Petrologica Sinica, 31(3): 747–758.
[] Yuan ZX, Bai G, Li SB, Guo YQ, Wang PH. 1987. Characteristics of feldspars from the Lingshan granite body and their petrogenetic significance. Acta Petrologica et Mineralogica, 6(4): 331–338.
[] Yuan ZX, Bai G, Yu SM, Li SB, Guo YQ, Liu GH. 1988. Geological Characteristics, Petrogenesis and Mineralization in Lingshan Granites, Jiangxi Province. Beijing: Beijing Science & Technology Press.
[] Zhang CZ, Li JT, Huang DT, Zou JR, Ouyang HD. 1985. The evolution and mineralization of Lingshan granite in Jiangxi. Geotectonica et Metallogenia, 9(4): 323–340.
[] Zhang P, Tian BS. 2005. Primary discussion on the characteristics and metallogenesis of Lingshan composite massif in northeastern Jiangxi Province. World Nuclear Geoscience, 22(1): 31–38.
[] Zhang Q, Pan GQ, Li CD, Jin WJ, Jia XQ. 2007. Granitic magma mixing versus basaltic magma mixing:New viewpoints on granitic magma mixing process:Some crucial questions on granite study (1). Acta Petrologica Sinica, 23(5): 1141–1152.
[] Zhang RX, Yang SY. 2016. A mathematical model for determining carbon coating thickness and its application in electron probe microanalysis. Microscopy and Microanalysis, 22(6): 1374–1380. DOI:10.1017/S143192761601182X
[] Zhang ZZ, Gu LX, Wu CZ, San JZ, Tang XQ, Xi AH, Wang S. 2005. Weiya complex, eastern Tianshan:Single-sourced or diverse-sourced? Evidence from biotite. Geochimica, 34(4): 328–338.
[] Zhao KD, Jiang SY, Jiang YH, Wang RC. 2005. Mineral chemistry of the Qitianling granitoid and the Furong tin ore deposit in Hunan Province, South China:Implication for the genesis of granite and related tin mineralization. European Journal of Mineralogy, 17(4): 635–648. DOI:10.1127/0935-1221/2005/0017-0635
[] Zhao KD, Jiang SY, Zhu JC, Li L, Dai BZ, Jiang YH, Ling HF. 2010. Hf isotopic composition of zircons from the Huashan-Guposhan intrusive complex and their mafic enclaves in northeastern Guangxi:Implication for petrogenesis. Chinese Science Bulletin, 55(6): 509–519. DOI:10.1007/s11434-009-0314-0
[] Zhao M, Yang SY, Zuo RG, Zhao KD, Jiang YH, Ling HF, Chen PR. 2015. Magmatic evolution characteristics of Xiangshan volcanic-intrusive complex from the Gan-Hang Belt:Studies on the mineral chemistry of plagioclase and biotite. Acta Petrologica Sinica, 31(3): 759–768.
[] Zheng JP, Li CN, Xue CS, Ye XC. 1996. Genesis of basalt inclusion in Lingshan granites, Jiangxi. Geological Science and Technology Information, 15(1): 19–24.
[] Zhou J, Jiang YH, Xing GF, Zeng Y, Ge WY. 2013. Geochronology and petrogenesis of Cretaceous A-type granites from the NE Jiangnan Orogen, SE China. International Geology Review, 55(11): 1359–1383. DOI:10.1080/00206814.2013.774199
[] Zhou ZX. 1988. Chemical characteristics of mafic mica in intrusive rocks and its geological meaning. Acta Petrologica Sinica, 4(3): 63–73.
[] Zhu ZY, Jiang SY, Hu J, Gu LX, Li JW. 2014. Geochronology, geochemistry, and mineralization of the granodiorite porphyry hosting the Matou Cu-Mo (±W) deposit, Lower Yangtze River metallogenic belt, eastern China. Journal of Asian Earth Sciences, 79: 623–640. DOI:10.1016/j.jseaes.2013.07.033
[] Zhu ZY, Wang RC, Che XD, Zhu JC, Wei XL, Huang XE. 2015. Magmatic-hydrothermal rare-element mineralization in the Songshugang granite (northeastern Jiangxi, China):Insights from an electron-microprobe study of Nb-Ta-Zr minerals. Ore Geology Reviews, 65: 749–760. DOI:10.1016/j.oregeorev.2014.07.021
[] 陈志中. 1984. 赣东北灵山岩体的阶段划分与钨钽成矿关系的研究. 大地构造与成矿学, 8(1): 59–70.
[] 东前, 杜杨松, 曹毅, 庞振山, 宋林旭, 郑震. 2011. 江西武山花岗闪长斑岩中黑云母成分特征及其成岩成矿意义. 矿物岩石, 31(2): 1–6.
[] 黄定堂. 2003. 灵山岩体演化特征及其与稀有金属的成矿关系. 地质与勘探, 39(4): 35–40.
[] 蒋少涌, 赵葵东, 姜耀辉, 凌洪飞, 倪培. 2006. 华南与花岗岩有关的一种新类型的锡成矿作用:矿物化学、元素和同位素地球化学证据. 岩石学报, 22(10): 2509–2516.
[] 蒋少涌, 李亮, 朱碧, 丁昕, 姜耀辉, 顾连兴, 倪培. 2008. 江西武山铜矿区花岗闪长斑岩的地球化学和Sr-Nd-Hf同位素组成及成因探讨. 岩石学报, 24(8): 1679–1690.
[] 林文蔚, 彭丽君. 1994. 由电子探针分析数据估算角闪石、黑云母中的Fe3+、Fe2+. 长春地质学院学报, 24(2): 155–162.
[] 刘友华, 李康东, 涂金飞. 2011. 江西葛源黄山铌 (钽) 矿床地质特征及成因分析. 资源调查与环境, 32(4): 291–298.
[] 吕志成, 段国正, 董广华. 2003. 大兴安岭中南段燕山期三类不同成矿花岗岩中黑云母的化学成分特征及其成岩成矿意义. 矿物学报, 23(2): 177–184.
[] 彭花明. 1997. 杨溪岩体中黑云母的特征及其地质意义. 岩石矿物学杂志, 16(3): 271–281.
[] 孙世华, 于洁. 1989. Mg-Fe云母化学成分的解释和分类 (Ⅱ)——Mg-Fe云母的自然分类. 地质科学(2): 176–189.
[] 田邦生, 周贤旭. 2008. 赣东北"金三角区"成矿带的划分及成矿模式. 地质找矿论丛, 23(2): 130–134.
[] 谢应雯, 张玉泉. 1987. 横断山不同成因类型花岗岩类岩石中黑云母的标型特征. 矿物学报, 7(3): 245–254.
[] 徐耀明, 蒋少涌, 朱志勇, 周巍, 孔凡斌, 孙明志, 熊永根. 2013. 江西九瑞矿集区成矿与未成矿中酸性侵入岩年代学、岩石化学、矿物化学特征的异同及地质意义. 岩石学报, 29(12): 4291–4310.
[] 杨水源, 蒋少涌, 赵葵东, 姜耀辉, 范洪海. 2012. 江西相山铀矿田邹家山矿床中流纹斑岩的锆石U-Pb年代学、岩石地球化学与Sr-Nd-Hf同位素组成. 岩石学报, 28(12): 3915–3928.
[] 杨水源, 蒋少涌, 赵葵东, 姜耀辉, 凌洪飞, 陈培荣. 2013. 江西相山铀矿田如意亭剖面火山岩的年代学格架及其地质意义. 岩石学报, 29(12): 4362–4372.
[] 杨水源, 文霞, 赵葵东, 姜耀辉, 凌洪飞, 陈培荣. 2015. 赣杭构造带新路火山盆地粗面斑岩的年代学、岩石成因及其地质意义. 岩石学报, 31(3): 747–758.
[] 袁忠信, 白鸽, 李双保, 郭永泉, 王佩华. 1987. 灵山花岗岩的长石特征及其成因意义. 岩石矿物学杂志, 6(4): 331–338.
[] 袁忠信, 白鸽, 余时美, 李双保, 郭永泉, 刘光海. 1988. 江西灵山花岗岩地质特征及其成岩成矿作用. 北京: 北京科学技术出版社.
[] 章崇真, 李吉涛, 黄定堂, 邹继容, 欧阳辉东. 1985. 江西灵山花岗岩的演化和成矿. 大地构造与成矿学, 9(4): 323–340.
[] 章平, 田邦生. 2005. 赣东北灵山复式岩体特征及其成矿作用初探. 世界核地质科学, 22(1): 31–38.
[] 张旗, 潘国强, 李承东, 金惟俊, 贾秀勤. 2007. 花岗岩混合问题:与玄武岩对比的启示——关于花岗岩研究的思考之一. 岩石学报, 23(5): 1141–1152.
[] 张遵忠, 顾连兴, 吴昌志, 三金柱, 汤晓茜, 郗爱华, 王硕. 2005. 东天山尾亚杂岩体:同源还是异源?——来自黑云母的证据. 地球化学, 34(4): 328–338.
[] 赵沔, 杨水源, 左仁广, 赵葵东, 姜耀辉, 凌洪飞, 陈培荣. 2015. 赣杭构造带相山火山侵入杂岩的岩浆演化特征——来自斜长石和黑云母的化学成分研究. 岩石学报, 31(3): 759–768.
[] 郑建平, 李昌年, 薛重生, 叶新才. 1996. 江西灵山花岗岩中玄武岩包体的成因. 地质科技情报, 15(1): 19–24.
[] 周作侠. 1988. 侵入岩的镁铁云母化学成分特征及其地质意义. 岩石学报, 4(3): 63–73.