西南三江普朗铜矿岩浆混合作用:矿物学和地球化学证据*

2015-03-15刘欢张长青贾福东周云满娄德波

岩石学报 2015年11期

刘欢张长青＊＊贾福东周云满娄德波

LIU Huan1，ZHANG ChangQing1＊＊，JIA FuDong1，ZHOU YunMan2 and LOU DeBo1

1. 中国地质科学院矿产资源研究所，国土资源部成矿作用与资源评价重点实验室，北京 100037

2. 云南黄金矿业集团股份有限公司，昆明 650224

1. MLR Key Laboratory of Metallogeny and Mineral Assessment，Institute of Mineral Resources，Chinese Academy of Geological Sciences，Beijing 100037，China

2. Yunnan Gold ＆ Mineral Group Co. ，Ltd，Kunming 650224，China

2015-02-10 收稿，2015-05-31 改回.

1 引言

普朗铜矿床位于西南三江义敦岛弧南段的格咱岛弧，是晚三叠世甘孜-理塘洋壳向中咱地块俯冲的产物(李文昌等，2010;邓军等，2011，2012;毛景文等，2012;Deng et al.，2014a，b;Hou et al.，2007)(图1)。矿床产于普朗复式侵入体内，储量巨大，Cu 资源量达650 万吨，伴生Au、Ag、Mo、Pb及Zn 等资源(刘学龙等，2013)。前人从岩石学、地球化学及年代学等方面对普朗复式岩体开展了大量研究，但对岩体的成因研究相对薄弱。曹殿华等(2009)、刘学龙和李文昌(2013)分别从矿物学、岩石地球化学及同位素地球化学等方面提出了普朗复式岩体具有岩浆混合成因的特点，Liu and Zhang (2014)从同位素年代学方面阐述了岩浆混合作用与成矿关系。岩浆混合作用作为岩浆演化的重要表现形式，不仅是形成火成岩多样性的主要机制(Didier and Barbarin，1991;Jarrar et al.，2013;Renjith et al.，2014;Ma et al.，2013;莫宣学等，2002)，也与Cu、Pb、Zn、Mo、Au、Sn、W 等成矿作用有着密切的联系(Wallace and Carmicheal，1992;Hou et al.，2011)，其研究对于探讨地球深部动力学过程、认识成岩成矿过程、分析岩浆演化等有着重要意义。然而，对于普朗铜矿岩浆混合作用，尤其是岩浆混合程度的研究相对薄弱。

图1 普朗斑岩铜矿地质图(据李文昌等，2013;云南迪庆有色金属有限责任公司，2009①修改)Fig.1 Geological mapping of Pulang porphyry copper deposit (modified after Li et al.，2013)

岩浆混合程度作为岩浆混合作用的重要研究内容之一，是指混合岩浆中基性岩浆所占的比例，它在数值上是确定的，但在实际应用时难以得到。传统的岩浆混合程度的研究主要基于以下三种方法:①云南迪庆有色金属有限责任公司. 2009. 云南省迪庆普朗铜矿区勘探地质报告质量平衡原理计算(董国臣等，2006);②定量模型估算(O’Hara，1977;Perugini et al.，2013);③二元混合方程确定。上述方法多涉及岩浆岩中微量元素的浓度和分配系数等参数，并且使用这些方法的前提条件是:研究的岩浆混合体系在能量和成分上达到平衡;然而，在实际的岩浆混合杂岩体中，均匀和非均匀的岩浆混合岩往往伴生在一起(Perugini and Poli，2012)，因此难以准确确定其岩浆混合程度。近年来的研究发现，岩浆混合的结构是一种混沌的性质，导致岩浆混沌混合的基本动力学机理是拉伸和褶曲作用，这两种作用都是非线性的耦合过程，分形是刻画这一过程的有效工具(Perugini and Poli，2000，2012;Perugini et al.，2003，2013)。

分形是研究自然界空间结构复杂性的一门学科，是复杂性科学的重要组成部分。分形理论最早由Mandelbrot(1983)提出，随后被广泛地应用于各个科学领域。运用分形理论的思想和方法，能从无序中发现有序，揭示杂乱、破碎、混沌等极不规则的复杂现象内部所蕴含的规律。分形理论主要包括自相似分形、自仿射分形与多重分形三类。多重分形作为分形理论的重要分支，其在形式上包括将自相似分形分解为多个相关联的简单分形，每个简单分形均有各自的奇异指数及相应的分维值，这两者构成了多重分形谱。地质现象在时间和空间上均表现出极强的复杂性和不规则性，借助多重分形理论能真正有效地解决地质现象的复杂性问题，对其内部结构进行精细解析(Cheng，1995，1999;Wang et al.，2010a，2011a，b)。多重分形理论已被广泛应用于刻画地质体的不规则分布(Agterberg et al.，1996;Deng et al.，2010，2011b;Liu et al.，2012;Wang et al.，2010b，2012;Zuo et al.，2013，2015)，然而该应用主要集中于宏观尺度，微观尺度上的应用相对薄弱。

文章以普朗铜矿的复式岩体及其内发育的暗色微粒包体为研究对象，基于系统的野外观测与显微镜观察，结合电子探针分析与主微量元素测试工作，从矿物学与地球化学方面阐述了包体的岩浆混合成因证据，证实了研究区存在岩浆混合作用;继而采用多重分形理论定量刻画了包体中主要元素的分布特征，厘定了岩浆混合程度与元素分布的多重分形参数间的关系，为岩浆混合程度的解析提供了定量化依据。

2 区域地质与矿床地质特征

格咱岛弧位于西南三江义敦岛弧南段，是古特提斯演化形成的重要地质构造单元(Wang et al.，2014)，区内构造发育、岩浆活动频繁，其印支期的俯冲造山作用诱发了大规模的火山岩浆活动，形成了以普朗为代表的一系列大型、超大型矿床(邓军等，2013，2014;Mao et al.，2014)。区域出露地层主要为晚三叠世曲嘎寺组(T3q)、图姆沟组(T3t)和喇嘛哑组(T3lm)(图1b)。区内断裂、褶皱主要由一系列NW 向线性褶皱与同向断裂及近EW 向断裂组成，晚期发育规模较小的NE 向断层;区内岩浆岩总体以岛弧安山岩及中酸性浅成-超浅成侵入岩为主。

3 包体的岩相学特征

普朗复式岩体中的石英闪长玢岩与石英二长斑岩内均发育大量随机分布的暗色微粒包体。包体大小不一(粒径在2 ～15cm 之间)，粒度较细，一般为灰黑色，明显暗于寄主岩，呈椭球状、长条状、不规则状等产出，少数具有塑性形变，与寄主岩多呈截然接触关系，少数呈过渡接触关系(图3a、b)。在包体与寄主岩的接触面，偶见斜长石巨晶(图3b)。显微镜观察表明:包体具岩浆结构，为岩浆结晶形成(图3c-f)。包体与寄主岩具有相近的矿物组成，但与寄主岩相比，包体中明显富集角闪石、黑云母等暗色矿物。

包体与寄主岩中均发育大量针状磷灰石(图3e、f)，为淬冷结晶作用的重要标志。包体与寄主岩中均可见具有暗色矿物镶边的眼球状石英，具有次生加大边的斜长石(图3g)。此外，还可见粗大的斜长石斑晶包裹角闪石、斜长石构成嵌晶结构(图3h)。包体与寄主岩中的斜长石具有复杂的成分与结构不平衡现象。在寄主岩中偶见斜长石具有富钙的核(An=51 ～42)和富钠的幔部(An=16)，且两者变化截然(图3g);包体中的斜长石斑晶，核部成分均一(An =36 ～34)，核部被富钠的斜长石(An =28)包裹，后者又被富钙的斜长石(An=35)包裹(图3h)。

图2 普朗斑岩铜矿1 线剖面图(云南迪庆有色金属有限责任公司，2009 修改)Fig.2 Vertical cross-section along exploration line 1 in Pulang porphyry copper deposit

4 样品的选取与分析方法

在系统的岩相学观察与鉴定的基础上，研究选取了10组相对新鲜的包体与寄主岩对进行了主微量元素组成测定。同时，研究选取了两件产于石英闪长玢岩中的暗色微粒包体，借助电子探针面分析技术，获得了暗色微粒包体内Mg、Fe、Ca、Al、K、Ti、P 及Ba 等8 种元素含量的趋势图;进而采用多重分形理论解析了这8 种元素的分布特征。基于篇幅，测试数据省略。

4.1 全岩主微量元素测试

主微量元素测试在北京大学造山带与地壳演化重点实验室完成。主量元素测试采用X 射线荧光光谱法(XRF)，所用仪器为美国热电公司生产的ARL ADVANT XP +，精度优于＜1%。微量元素及稀土元素测试采用电感耦合等离子体质谱法(ICP-MS)在美国安捷伦科技有限公司生产的Agilent ICPMS 7500ce 型质谱仪上完成，精度优于＜5%。

4.2 电子探针分析

电子探针面分析测试在中国地质科学院矿产资源研究所电子探针实验室完成，测试仪器为日本JOEL 公司生产的JXA-823V，加速器电压为15kV，束电流为20nA，束斑直径为1μm，扫描面积为7.08mm2。使用硬玉(Si、Na 和Al)，镁橄榄石(Mg)，正长石(K)，磷灰石(P)，硅灰石(Ca)，金红石(Ti)，以及合成氧化物(Cr、Mn、Fe、Ni)作为标准。

4.3 多重分形计算过程

基于数据的特点，本文采用数盒子法对暗色微粒包体中Mg、Fe、Ca、Al、K、Ti、P 及Ba 元素的分布特征进行了多重分形解析。数盒子法是将数据集划分为多个互不重叠的小盒子以构成多尺度模式。用n(ε)个边长为ε 的盒子不重叠地覆盖数据集(Halsey et al.，1986)，则每个盒子中的测度(μj)可表示为:

式中，mj表示第j 个盒子中元素的含量和，mT表示所有盒子元素的总含量，n(ε)=N/ε。随后，配分函数χ(q，ε)可由下式求得:

图3 普朗斑岩铜矿寄主岩与暗色微粒包体的野外照片(a-b)与典型岩相学显微照片(c-h)MME-暗色微粒包体;Qtz-石英;Pl-斜长石;Bt-黑云母;Amp-角闪石;Ap-磷灰石;Px-辉石;Pnn-叶绿泥石. (g)与(h)中白色圈内数字表示钙长石分子(An)Fig.3 Field photographs (a，b)and photomicrographs (c-h)of host rocks and mafic microgranular enclaves from Pulang porphyry copper depositQtz-quartz;Pl-plagioclase;Bt-biotite;Amp-amphibole;Ap-apatite;Px-pyroxene;Pnn-penninite. The An contents are shown in circles

其中，

权重q 为测量不同位置的奇异指数提供了更准确更细微的方法。q ＞1，μq主要指示高含量的变化，q ＜1 则主要反映低含量的分布，q=1 时，μ 指示原始数据的分布。

色谱条件：Shimadzu-GL ODS-2色谱柱（250 mm×4.6 mm，5 μm）；流动相为甲醇-水（75∶25）；体积流量1.0 mL/min；检测波长228 nm；进样量20 μL，柱温30 ℃。

权重q 的质量值数τ(q)可表示为(Hentschel and Procaccia，1983):

广义维Dq可由以下幂率关系产生(Hentschel and Procaccia，1983):

权重q 与相应的广义维Dq可构成Rényi 谱。对于给定的q 的范围，Dq的变化越大，表明相应数据集的空间不均匀程度越高(Arias et al.，2011)。当q =2 时，D2为关联维，且与τ(2)相等。D2值越高，表明元素含量序列中高值聚集分布，反之，则指示元素含量的高值离散分布(Liu et al.，2012)。

依据公式(4)得到的质量指数τ(q)，通过勒让德变换可求得多重分形谱的奇异指数α(q)及其相应的分维值f(α)。

α(q)与f(α)构成了多重分形谱，其形状为倒钟形。当q为最大值、最小值时，可分别求得amin与f(amin)及amax与f(amax)。分形谱的宽度及其分形谱两端的高差可表示为:

Δα 与Δf 是描述多重分形谱两个必不可少的参数。随着Δα 增加，元素的分布由均匀(随机)过渡为不均匀(有序、复杂、聚集)。Δf 为正值时，分形谱的形状为右钩，数据集中低值占主导地位;Δf 为负值时，分形谱呈左钩，数据集中以高值为主。若高值与低值的比例相当，则分形谱呈左右对称的钟形(Deng et al.，2011b;Wang et al.，2011a，b)。

5 结果分析

5.1 元素地球化学特征

主量元素结果显示，暗色微粒包体的SiO2含量为53.67% ～61.50%，均值为58.35%，比寄主岩(SiO2含量为60.74% ～66.13%，均值为64.78%)偏基性。包体中CaO、MgO 和Fe2OT3等含量分别为3.14% ～6.69%、3.12% ～5.40%与3.38% ～9.00%，均值分别为4.71%、4.39% 与6.48%，比寄主岩中CaO(2.36% ～6.08%)、MgO(1.78% ～3.59%)和Fe2OT3(2.33% ～6.06%)的含量高(图4)。在SiO2与主要氧化物的二元图解中，暗色微粒包体与寄主岩SiO2含量与TiO2、Fe2OT3与MgO 的含量具有良好的线性关系，均表现为随着SiO2含量的升高而降低;K2O 与Na2O 含量随SiO2含量的增加变化不明显，表明暗色微粒包体与寄主岩的K2O 与Na2O 含量相似。

稀土元素球粒陨石标准化配分曲线显示，暗色微粒包体与其寄主岩具有相似的稀土元素配分模式，呈明显的右倾斜型，均表现为轻稀土元素富集，重稀土元素亏损，表明两者经历了相似的岩浆演化过程，具有密切的成因联系(图5a、b)。然而，与寄主岩相比，暗色微粒包体中稀土元素总量普遍较高(119.0 ×10-6～308.9 ×10-6)，并具有较明显的Eu 负异常(δEu=0.56 ～0.95，均值为0.72)。

在原始地幔标准化的微量元素蛛网图中，暗色微粒包体与其寄主岩具有相似的微量元素组成，均表现为富集大离子亲石元素(K、Rb、Ba、Sr)，亏损高场强元素(Nb、Hf、P、Ti)，无明显的δCe 异常(图5c、d)，具有典型的岛弧岩浆岩的微量元素特征。然而与寄主岩相比，暗色微粒包体的富集程度更高。此外，在主微量元素的二元图解(SiO2-V)中，暗色微粒包体与寄主岩显示出良好的线性关系(图6a)，表明两者经历了相似的演化过程，在V-Ba/V 图中，两者显示了曲线变化关系(图6b)，在一定程度上指示了岩浆混合过程。

表1 包体中元素分布的多重分形参数一览表Table 1 Multifractal parameters elemental distribution for mafic microgranular enclaves

5.2 包体内元素分布多重分形特征

借助电子探针面分析技术，获得了暗色微粒包体(样品编号为PL1-MME 与PL2-MME)内Mg、Fe、Ca、Al、K、Ti、P 及Ba 等8 种元素含量的趋势图，其中Ca 元素的分布图如图7所示。采用数盒子法计算得到了PL1-MME 与PL2-MME 内元素分布的多重分形谱(图8、图9)，与多重分形参数(表1)。PL1-MME 中，Al、Ba、Ca 及K 元素分布的多重分形谱呈左钩，Fe、Mg、P 及Ti 元素分布的多重分形谱呈右钩，表明该暗色微粒包体中Al、Ba、Ca 及K 元素的含量以高值为主，Fe、Mg、P 及Ti 元素的含量以低值居多。此外，Al、Ba 及K 元素分布的关联维D2相对于其他元素较高，指示这三种元素的高值聚集分布;Al、Ca、Fe 及K 元素的多重分形谱相对于其他元素较宽，表明相应元素分布的多重分形特征较明显。PL2-MME 中，除元素P 与Ti 之外，其他元素分布的多重分形谱均呈左钩。Al 与Ba 元素分布的关联维D2相对于其他元素较高，Ca、Fe、K 与Mg 元素分布的多重分形谱相对于其他元素较宽。不同暗色微粒包体中同一元素的分布的多重分形参数对比显示，与PL2-MME 相比，在PL1-MME 中Ba、K、P及Ti 元素分布的关联维较高，Al 与K 元素分布的多重分形谱较宽。Ti、K 及P 在不同暗色微粒包体中分布的关联维与多重分形谱的宽度相差较大(图10)。

图4 普朗铜矿暗色微粒包体与寄主岩(石英闪长玢岩或石英二长斑岩)的哈克图解Fig.4 Harker plots of mafic microgranular enclaves (MME)and their host rocks (quartz diorite porphyry or quartz monzonite porphyry)from Pulang porphyry copper deposit

6 讨论

6.1 普朗复式岩体的岩浆混合成因

普朗铜矿的石英闪长玢岩与石英二长斑岩中发育的暗色微粒包体具有典型的岩浆岩矿物组合和结构构造特点，未见石榴子石、堇青石、红柱石等典型的富铝矿物，亦不具有变晶结构、交代结构和面理构造等变质岩常见的特征，因此不是地壳深部变质岩的难熔残留体。研究区的暗色微粒包体与捕掳体在形态上有较大差别，围岩捕掳体多为棱角状，但暗色包体形态多为浑圆状和椭球状，并具有塑性流变特点。以上特点可排除暗色微粒包体的围岩捕掳体成因。同源岩浆早期结晶形成的析离体或堆积体成因的包体，其同种矿物的粒度与寄主岩相近，但暗色微粒包体的粒度却明显小于其寄主岩。此外，与寄主岩相比，包体具有更高的稀土元素总量，表明包体不是寄主岩浆早期的堆晶体。因为REE 为强不相容元素，倘若暗色微粒包体是花岗质岩浆早期结晶分异产物的堆积体，则其REE 含量应比寄主岩低，REE 配分曲线应当位于寄主岩的下方，故又可以排除包体的析离体成因。综合以上特点，研究区发育的暗色微粒包体为岩浆混合成因的铁镁质微粒包体，这为该区岩浆混合作用的存在提供了有效证据。

图5 普朗斑岩铜矿中暗色微粒包体与寄主岩的球粒陨石标准化稀土元素配分曲线(a、b)及原始地幔标准化微量元素蛛网图(c、d)(标准化值据Sun and McDonough，1989)Fig.5 Chondrite-normalized REE (a，b)and primitive-mantle-normalized multi-element (c，d)patterns for mafic microgranular enclaves (MME)and their host rocks (quartz diorite porphyry or quartz monzonite porphyry)from Pulang porphyry copper deposit(normalizing values after Sun and McDonough，1989)

图6 普朗铜矿暗色微粒包体与寄主岩(石英闪长玢岩或石英二长斑岩)的主微量二元图解Fig.6 Binary diagrams of mafic microgranular enclaves (MME)and their host rocks (quartz diorite porphyry or quartz monzonite porphyry)from Pulang porphyry copper deposit

图7 普朗斑岩铜矿暗色微粒包体的电子探针面分析范围的背散射图像(a、b)与其内Ca 元素分布图(c、d)Fig.7 Backscattered electron images (a，b)and the distribution of element Ca (c，d)for mafic microgranular enclaves (MME)from Pulang porphyry copper deposit

暗色微粒包体具聚塑性流动特征，表明包体经历了与寄主岩共同的液态共存阶段，暗示发生过液态物质交换-混合作用(mixing)和混和作用(mingling)，混合作用改变了暗色微粒包体的成分，混和作用造成了包体呈不同的大小和形态散布到寄主岩中。包体中存在的具细粒暗色矿物(多数为角闪石)镶边的眼球状石英，可能是寄主花岗岩浆中的斑晶，在铁镁质岩浆与花岗质岩浆发生混合作用时，被铁镁质岩浆捕获，并发生反应。石英在基性岩浆中是不稳定的矿物相，当被基性岩浆捕获时，其边角将被熔化。熔融作用的吸热效应在紧靠石英边缘的一圈熔体中形成了局部过冷条件，导致细粒暗色矿物集合体围绕石英捕掳晶晶出(莫宣学等，2002)。粗大的斜长石斑晶包裹角闪石、斜长石构成嵌晶结构。这是由于偏基性的岩浆或正在混合的岩浆体系骤冷，产生很多小颗粒角闪石和斜长石，而后缓慢冷却至外部环境温度，平衡熔体冷却到接近液相线温度时斜长石晶核缓慢结晶成粗大晶体，形成包裹着早期形成的角闪石和斜长石的嵌晶结构(Zhao et al.，2012)。在包体和寄主岩石中，均可见到部分斜长石具有富钙的残留核被富钠的斜长石包裹，且两者变化截然。富钙的残留斜长石可能代表早期从基性岩浆中结晶的斜长石，边部相对富钠的斜长石可能代表晚期从混浆中结晶或从壳源花岗质岩浆中结晶的斜长石。斜长石的这种成分与结构不平衡现象，有效地记录了基性岩浆与酸性岩浆发生混合的过程。暗色微粒包体与寄主岩具有相似的微量元素组成与稀土元素配分模式，在SiO2与主要氧化物、部分微量元素相关图解中，两者呈明显的线性相关，表明两者经历了相似的岩浆演化过程，具有一定的亲缘性，这种亲缘性反映包体与寄主岩可能具有同源性，包体可能是同源岩浆早期结晶的析离体或堆积体，也有可能是岩浆混合作用过程中物质交换的结果。但是在正常的岩浆分异演化过程中，随着结晶分异的加强，残余熔体的δEu 负异常将会愈加强烈，同源岩浆早期结晶产物的δEu 负异常程度将显著小于晚期产物。然而，与寄主岩相比，包体具有较明显的δEu 负异常。因此，我们认为，是岩浆混合过程中物质交换导致了包体与寄主岩的亲缘性。

6.2 多重分形参数与岩浆混合程度间的关系

暗色包体中CaO 含量的变化可指示岩浆混合程度(Perugini and Poli，2000)，随着岩浆混合作用的进行，暗色微粒包体中CaO 含量逐渐降低。PL1-MME 中CaO 含量为6.69%，PL2-MME 中CaO 含量为4.92%，表明PL2-MME 中岩浆混合程度更高。基于此，我们对比分析了两件包体样品中Mg、Fe、Ca、Al、K、Ti、P 及Ba 元素分布的多重分形参数与CaO 含量的关系。结果显示，多数元素分布的多重分形谱宽度△a 与CaO 含量成反相关，分形谱两端的高差△f(a)与CaO 含量成正相关，关联维D2与CaO 含量关系较为复杂，主体呈正相关关系。这表明，暗色微粒包体中主量元素分布的多重分形谱越宽，数据集中高值与低值所占的比例相差越小，高值聚集程度越低，则岩浆混合程度越高。基性与酸性岩浆混合过程中的拉伸与褶曲导致了暗色微粒包体的形态具有多重分形特征，由于包体的形态难以定量刻画，因此采用包体中主要元素的分布图来确定包体的形态分形特征(Perugini and Poli，2000;Perugini et al.，2013)。随着岩浆混合作用的进行，包体的体积越来越小，形态越来越不规则，这种不规则表现为多重分形特征越来越明显，对应于包体中主要元素分布的多重分形谱宽度增大。在混合过程中，酸性岩浆与基性岩浆之间存在物质(斑晶与元素)交换，从而导致两者间元素的差异性相对减小，表现为主要元素高值的聚集分布程度降低，数据集中高值与低值所占的比例的差异性相比于混合前降低。

图8 暗色微粒包体PL1-MME 中元素分布的多重分形谱(a)Al;(b)Ba;(c)Ca;(d)Fe;(e)K;(f)Mg;(g)P;(h)TiFig.8 Multifractal spectrum of elemental distribution for PL1-MME

不同暗色微粒包体间，Ti、K 与P 元素的多重分形参数相对于其他元素变化较大，表明随着岩浆混合作用的进行，Ti、K 与P 元素分布的不均一变化较大，这三种元素在指示岩浆混合程度上更灵敏。这可能与基性岩浆的组成有关，因为Ti、K 与P 是组成基性岩浆的重要元素。

图9 暗色微粒包体PL2-MME 中元素分布的多重分形谱(a)Al;(b)Ba;(c)Ca;(d)Fe;(e)K;(f)Mg;(g)P;(h)TiFig.9 Multifractal spectrum of elemental distribution for PL2-MME

7 结论

(1)普朗斑岩铜矿石英闪长玢岩与石英二长斑岩中发育大量随机分布的暗色微粒包体，包体具岩浆结构，粒度较细，呈椭圆状、不规则状产出，偶见塑性变形，富集暗色矿物;包体与寄主岩中可见具暗色矿物镶边的眼球状石英、针状磷灰石以及具有结构与成分不平衡现象的斜长石斑晶，包体与寄主岩具有相似的微量元素组成与稀土元素配分模式，这些现象共同佐证了研究区岩浆混合作用的存在。

图10 包体中元素分布的多重分形参数折线图Fig.10 Curves of multifractal parameters elemental distribution for mafic microgranular enclaves

(2)证实了多重分形在岩浆混合作用研究中的有效性。研究结果表明，多重分形参数能有效刻画岩浆混合程度。包体中Al、Ca 及K 等主要元素分布的多重分形谱越宽，关联维越小，则岩浆混合程度越高，反之亦然。暗色微粒包体中Ti、P 与K 元素的分布能更有效的反应岩浆混合程度。

(3)文章丰富了岩浆混合作用及多重分形在微观尺度的研究内容。

致谢野外与室内研究工作得到了毛景文老师的悉心指导和彭惠娟、肖昌浩同学的无私帮助;实验测试工作得到了中国地质科学院矿产资源研究所电子探针实验室的陈振宇老师及陈小丹硕士以及北京大学造山带与地壳演化重点实验室的大力支持与帮助;两名审稿人与主编对文章提出了宝贵的意见;在此一并表达最诚挚的感谢!

Agterberg FP，Cheng QM，Brown A and Good D. 1996. Multifractal modeling of fractures in the Lac du Bonnet Batholiths，Manitoba.Computer ＆ Geosciences，22(5):497 -507

Arias M，Gumiel P，Sanderson DJ and Martin-Izard A. 2011. A multifractal simulation model for the distribution of VMS deposits in the Spanish segment of the Iberian pyrite belt. Computers ＆Geosciences，37(12):1917 -1927

Cao DH，Wang AJ，Li WC，Wang GS，Li RP and Li YK. 2009. Magma mixing in the Pulang porphyry copper deposit:Evidence from petrology and element geochemistry. Acta Geologica Sinica，83(2):166 -175 (in Chinese with English abstract)

Cheng QM. 1995. The perimeter-area fractal model and its application to geology. Mathematical Geology，27(1):69 -82

Cheng QM. 1999. The gliding box method for multifractal modeling.Computers ＆ Geosciences，25(9):1073 -1079

Deng J，Wang QF，Yang LQ，Wang YR，Gong QJ and Liu H. 2010.Delineation and explanation of geochemical anomalies using fractal models in the Heqing area，Yunnan Province，China. Journal of Geochemical Exploration，105(3):95 -105

Deng J，Wang QF，Wan L，Liu H，Yang LQ and Zhang J. 2011. A multifractal analysis of mineralization characteristics of the Dayingezhuang disseminated-veinlet gold deposit in the Jiaodong gold province of China. Ore Geology Reviews，40(1):54 -64

Deng J，Yang LQ and Wang CM. 2011. Research advances of superimposed orogenesis and metallogenesis in the Sanjiang Tethys.Acta Petrologica Sinica，27(9):2501 - 2509 (in Chinese with English abstract)

Deng J，Wang CM and Li GJ. 2012. Style and process of the superimposed mineralization in the Sanjiang Tethys. Acta Petrologica Sinica，28(5):1349 -1361 (in Chinese with English abstract)

Deng J，Ge LS and Yang LQ. 2013. Tectonic dynamic system and compound orogeny:Additionally discussing the temporal-spatial evolution of Sanjiang orogeny，Southwest China. Acta Petrologica Sinica，29(4):1099 -1114 (in Chinese with English abstract)

Deng J，Wang QF，Li GJ，Li CS and Wang CM. 2014a. Tethys tectonic evolution and its bearing on the distribution of important mineral deposits in the Sanjiang region，SW China. Gondwana Research，26(2):419 -437

Deng J，Wang QF，Li GJ and Santosh M. 2014b. Cenozoic tectonomagmatic and metallogenic processes in the Sanjiang region，southwestern China. Earth-Science Reviews，138:268 -299

Deng J，Wang CM，Li WC，Yang LQ and Wang QF. 2014. The situation and enlightenment of research of the tectonic evolution and metallogenesis in the Sanjiang Tethys. Earth Science Frontiers，21(1):52 -64 (in Chinese with English abstract)

Deng J，Wang QF，Li GJ and Zhao Y. 2015a. Structural control and genesis of the Oligocene Zhenyuan orogenic gold deposit，SW China.Ore Geology Reviews，65:42 -54

Deng J，Wang QF，Li GJ，Hou ZQ，Jiang CZ and Danyushevsky L.2015b. Geology and genesis of the giant Beiya porphyry-skarn gold deposit，northwestern Yangtze Block，China. Ore Geology Reviews，70:457 -485

Didier J and Barbarin B. 1991. Enclaves and Granite Petrology.Amsterdam:Elsevier，625

Dong GC，Mo XX，Zhao ZD，Zhu DC，Wang LL，Chen T and Li B.2006. Magma mixing in middle part of Gangdise magma belt:Evidences from granitoid complex. Acta Petrologica Sinica，22(4):835 -844 (in Chinese with English abstract)

Halsey TC，Jensen MH，Kadanoff LP，Procaccia I and Shraiman BI.1986. Fractal measures and their singularities，the characterization of strange sets. Physical Review (A)，33(2):1141 -1151

Hentschel HGE and Procaccia L. 1983. The infinite number of generalized dimensions of fractals and strange attractors. Physica(Series D)，8(3):435 -444

Hou ZQ，Zaw K，Pan GT，Mo XX，Xu Q，Hu YZ and Li XZ. 2007.Sanjiang Tethyan metallogenesis in S. W. China:Tectonic setting，metallogenic epochs and deposit types. Ore Geology Reviews，31(1-4):48 -87

Hou ZQ，Zhang HR，Pan XF and Yang ZM. 2011. Porphyry (Cu-Mo-Au)deposits related to melting of thickened mafic lower crust:Examples from the eastern Tethyan metallogenic domain. Ore Geology Reviews，39(1 -2):21 -45

Jarrar GH，Yaseen N and Theye T. 2013. A hybrid composite dike suite from the northern Arabian Nubian Shield， southwest Jordan:Implications for magma mixing and partial melting of granite by mafic magma. Journal of Volcanology and Geothermal Research，254:80-93

Li WC，Pan GT，Hou ZQ，Mo XX，Wang LQ，Ding J and Xu Q. 2010.Archipelagic-Basin， Forming Collision Theory and Prospecting Techniques along the Nujiang-Lancangjiang-Jinshajiang Area in Southwestern China. Beijing:Geological Publishing House (in Chinese)

Li WC，Zeng PS，Hou ZQ and White NC. 2011. The Pulang porphyry copper deposit and associated felsic intrusions in Yunnan Province，Southwest China. Economic Geology，106(1):79 -92

Li WC，Yu HJ and Yin GH. 2013. Porphyry metallogenic system of Geza arc in the Sanjiang region，southwestern China. Acta Petrologica Sinica，29(4):1129 -1144 (in Chinese with English abstract)

Liu H，Wang QF，Li GJ and Wan L. 2012. Characterization of multitype mineralizations in the Wandongshan gold poly-metallic deposit，Yunnan (China)，by fractal analysis. Journal of Geochemical Exploration，122:20 -33

Liu H and Zhang CQ. 2014. Magma mixing and its genetic relationship with mineralization in the Pulang porphyry copper deposit. Acta Geologica Sinica，88(Suppl.2):565 -567

Liu XL and Li WC. 2013. The Indo-Chinese epoch magmatism in Geza arc of Yunnan:Evidences from zircon U-Pb dating and Hf isotopic composition. Earth Science Frontiers，20(5):57 -74 (in Chinese with English abstract)

Liu XL，Li WC，Yin GH and Zhang N. 2013. The geochronology，mineralogy and geochemistry study of the Pulang porphyry copper deposits in Geza arc of Yunnan Province. Acta Petrologica Sinica，29(9):3049 -3064 (in Chinese with English abstract)

Ma XH，Chen B and Yang MC. 2013. Magma mixing origin for the Aolunhua porphyry related to Mo-Cu mineralization，eastern Central Asian Orogenic Belt. Gondwana Research，24(3 -4):1152 -1171

Mandelbrot BB. 1983. The Fractal Geometry of Nature. San Francisco:Freeman

Mao JW，Xie GQ，Duan C，Pirajno F，Ishiyama D and Chen YC.2011a. A tectono-genetic model for porphyry-skarn-stratabound Cu-Au-Mo-Fe and magnetiteapatite deposits along the Middle-Lower Yangtze River Valley，Eastern China. Ore Geology Reviews，43(1):294 -314

Mao JW，Zhang JD，Pirajno F，Ishiyama D，Su HM，Guo CL and Chen YC. 2011b. Porphyry Cu-Au-Mo-epithermal Ag-Pb-Zn-distal hydrothermal Au deposits in the Dexing area，Jiangxi Province，East China:A linked ore system. Ore Geology Reviews，43(1):203-216

Mao JW，Zhou ZH，Feng CY，Wang YT，Zhang CQ，Peng HJ and Yu M. 2012. A preliminary study of the Triassic large-scale mineralization in China and its geodynamic setting. Geology in China，39(6):1437 -1471 (in Chinese with English abstract)

Mao JW，Pirajno F，Lehmann B，Luo MC and Berzina A. 2014.Distribution of porphyry deposits in the Eurasian continent and their corresponding tectonic settings. Journal of Asian Earth Sciences，79:576 -584

Mo XX，Luo ZH，Xiao QH，Yu XH，Liu CD，Zhao ZD and Zhou S.2002. Evidence of Magma mixing in granitoids plutons and the way of investigation on granitoids. In:Xiao QH (ed. ). The Ways of Investigation on Granitoids. Beijing:Geological Publishing House，53 -70 (in Chinese)

O’Hara MJ. 1977. Geochemical evolution during fractional crystallization of a periodically refilled magma chamber. Nature，266(5602):503-507

Perugini D and Poli G. 2000. Chaotic dynamics and fractals in magmatic interaction processes:A different approach to the interpretation of mafic microgranular enclaves. Earth and Planetary Science Letters，175(1 -2):93 -103

Perugini D，Poli G and Mazzuoli R. 2003. Chaotic advection，fractals and diffusion during mixing of magmas:Evidence from lava flows.Journal of Volcanology and Geothermal Research，124(3 -4):255-279

Perugini D and Poli G. 2012. The mixing of magmas in plutonic and volcanic environments:Analogies and differences. Lithos，153:261-277

Perugini D，De Campos CP，Dingwell DB and Dorfman A. 2013.Relaxation of concentration variance:A new tool to measure chemical element mobility during mixing of magmas. Chemical Geology，335:8 -23

Renjith ML，Charan SN，Subbarao DV，Babu EVSSK and Rajashekhar VB. 2014. Grain to outcrop-scale frozen moments of dynamic magma mixing in the syenite magma chamber，Yelagiri alkaline complex，South India. Geoscience Frontiers，5(6):801 -820

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 Ocean Basins. Geological Society，London，Special Publications，42(1):313 -345

Wallace P and Carmicheal ISE. 1992. Sulfur in basaltic magmas.Geochimica et Cosmochimica Acta，56(5):1863 -1874

Wang QF，Deng J，Liu H，Yang LQ，Wan L and Zhang RZ. 2010a.Fractal models for ore reserve estimation. Ore Geology Reviews，37(1):2 -14

Wang QF，Deng J，Zhao J，Wan L，Gong QJ，Yang LQ and Liu H.2010b. Tonnage-cutoff model and average grade-cutoff model for a single ore deposit. Ore Geology Reviews，38(1 -2):113 -120

Wang QF，Deng J，Liu H，Wan L and Zhang ZJ. 2011a. Fractal analysis of the ore-forming process in a skarn deposit:A case study in the Shizishan area，China. In:Sial AN，Bettencourt JS，De Campos CP and Ferreira VP (eds. ). Granite-Related Ore Deposits. Geological Society，London，Special Publications，350(1):89 -104

Wang QF，Deng J，Zhang QZ，Liu H，Liu XF，Wan L，Li N，Wang YR，Jiang CZ and Feng YW. 2011b. Orebody vertical structure and implications for ore-forming processes in the Xinxu bauxite deposit，western Guangxi，China. Ore Geology Reviews，39(4):230 -244 Wang QF，Deng J，Zhao JC，Li N and Wan L. 2012. The fractal relationship between orebody tonnage and thickness. Journal of Geochemical Exploration，122:4 -8

Wang QF，Deng J，Li CS，Li GJ，Yu L and Qiao L. 2014. The boundary between the Simao and Yangtze blocks and their locations in Gondwana and Rodinia:Constraints from detrital and inherited zircons. Gondwana Research，26(2):438 -448

Zhao KD，Jiang SY，Yang SY，Dai BZ and Lu JJ. 2012. Mineral chemistry，trace elements and Sr-Nd-Hf isotope geochemistry and petrogenesis of Cailing and Furong granites and mafic enclaves from the Qitianling batholith in the Shi-Hang zone，South China.Gondwana Research，22(1):310 -324 Zuo RG，Xia QL and Zhang DJ. 2013. A comparison study of the C-A

and S-A models with singularity analysis to identify geochemical

anomalies in covered areas. Applied Geochemistry，33:165 -172 Zuo RG，Wang J，Chen GX and Yang MG. 2015. Identification of weak anomalies:A multifractal perspective. Journal of Geochemical Exploration，148:12 -24

附中文参考文献

曹殿华，王安建，李文昌，王高尚，李瑞萍，李以科. 2009. 普朗斑岩铜矿岩浆混合作用:岩石学及元素地球化学证据. 地质学报，83(2):166 -175

邓军，杨立强，王长明. 2011. 三江特提斯复合造山与成矿作用研究进展. 岩石学报，27(9):2501 -2509

邓军，王长明，李龚健. 2012. 三江特提斯叠加成矿作用样式及过程. 岩石学报，28(5):1349 -1361

邓军，葛良胜，杨立强. 2013. 构造动力体制与复合造山作用——兼论三江复合造山带时空演化. 岩石学报，29(4):1099 -1114

邓军，王长明，李文昌，杨立强，王庆飞. 2014. 三江特提斯复合造山与成矿作用研究态势及启示. 地学前缘，21(1):52 -64

董国臣，莫宣学，赵志丹，朱弟成，王亮亮，陈涛，李冰. 2006. 冈底斯岩浆带中段岩浆混合作用——来自花岗杂岩的证据. 岩石学报，22(4):835 -844

李文昌，潘桂棠，侯增谦，莫宣学，王立全，丁俊，徐强. 2010. 西南“三江”多岛弧盆-碰撞造山成矿理论与勘查技术. 北京:地质出版社

李文昌，于海军，尹光候. 2013. 西南"三江"格咱岛弧斑岩成矿系统. 岩石学报，29(4):1129 -1144

刘学龙，李文昌. 2013. 云南格咱岛弧印支期岩浆作用的锆石年龄和铪同位素证据. 地学前缘，20(5):57 -74

刘学龙，李文昌，尹光侯，张娜. 2013. 云南格咱岛弧普朗斑岩型铜矿年代学、岩石矿物学及地球化学研究. 岩石学报，29(9):3049 -3064

毛景文，周振华，丰成友，王义天，张长青，彭惠娟，于淼. 2012.初论中国三叠纪大规模成矿作用及其动力学背景. 中国地质，39(6):1437 -1471