Identification and petrogenesis of Wuhema highly fractionated I-type granitoids in Qinghai Province from Songpan-Ganzi Terrane
-
摘要:
为探讨松潘−甘孜地体内青海吾和玛花岗质岩石的成因类型与岩石成因，对其开展了岩相学、地球化学、锆石U−Pb年代学和Lu-Hf同位素等研究。吾和玛花岗质岩石被确定为高分异的I型花岗岩，其结晶年龄为216.6±4.3 Ma，为正长花岗岩，具有弱过铝质高钾钙碱性特征。岩石轻、重稀土元素分馏不明显 ((La/Yb)N <7)， 轻稀土元素轻微富集，重稀土元素分布平缓，具有明显的负Eu异常(δEu=0.08 ~ 0.43)，明显亏损Ba、Sr和高场强元素Nb 、Ta 、P 、Ti ，富集大离子亲石元素Rb、Th、U、K，εHf(t)值为−2.54～1.32。吾和玛花岗质岩石的母岩浆是岩石圈拆沉作用引起软流圈上涌，诱发中元古代地壳物质（类似于变杂砂岩的物质）部分熔融形成长英质母岩浆，后期又受到斜长石、钾长石、磷灰石等矿物分离结晶作用的控制，经历了高程度的分异形成的。
Abstract:In order to study the genetic type and petrogenesis of the Wuhema granitoids in the Songpan Ganzi Terrane, the petrography, geochemistry, zircon U−Pb chronology and Lu-Hf isotope have been studied in this paper. The Wuhema granitic rocks are identified as highly fractionated I-type granitoids with crystallization age of 216.6±4.3 Ma. The rocks are syenogranite, and have the characteristics of weakly peraluminous, high-K calc-alkaline. The fractionation of REE in the rocks is not obvious (La/Yb)N <7), light rare earth elements are slightly enriched, and heavy rare earth elements are gently distributed. The rocks have obvious negative Eu anomalies (δEu= 0.08 ~ 0.43), with obvious depletion of Ba, Sr and high field strength elements (HFSEs, e.g., Nb, Ta, P and Ti), and enrichment of large ion lithophile elements (LILEs, e.g., Rb, Th and U), εHf(t) values of the rocks range from -2.54 to 1.32. The primary magma of the Wuhema granitic rocks were formed by the upwelling of asthenosphere caused by the lithosphere delamination, and induced the partial melting of Mesoproterozoic crustal materials (similar to metagreywacke), then the magma went through a high degree of fractionation of plagioclase, potassium feldspar, apatite and other minerals.
-
矿物识别是采矿工程、地质学等相关学科研究的基础。国内目前主要的矿物识别方法是依据人工经验，通常是依据其颜色、硬度和透明度，同时结合结晶习性、解理、断口等物理特征判断矿物种类（赵明，2010）。这对于从事矿物识别及相关工作人员的专业知识与经验要求较高，而且对矿石识别的准确度较低，难以大规模应用于实际场景。近年来，为了提高矿石识别效率及解决识别准确度低的问题，以机器学习为代表的矿物智能识别成为了研究的热点(周永章等，2018；郭艳军等，2020；郝慧珍等，2021）。
常规的矿物识别方法主要是依据矿物的属性，如外表形态、物理性质、化学成分等属性特征及相应的组合，使用能量散射光谱（EDS）(El Haddad et al., 2019)、电子探针(Tsuji et al., 2010)、激光诱导击穿光谱（LIBS）(Khajehzadeh et al., 2016)、显微光学观察(Maitre et al., 2019)、光谱分析(张旭等，2020)等方法实现矿物与岩石的识别。这些方法虽然准确率较高，但是需要用到专业仪器，过程复杂，整体效率偏低。
近年来，深度学习在矿物识别领域取得了良好效果（Trejbal et al., 2020；许振浩等，2021；周永章等，2021）。Baykan et al.（2010）利用多层感知机（MLP，Multilayer Perceptron）识别石英、白云母、黑云母等矿物。Izadi et al.（2013）利用两层神经网络，成功识别14种矿物。徐述腾等（2018）设计了Unet模型，建立对矿石矿物镜下照片的智能识别。Liu et al. （2019）基于深度学习和迁移学习算法，以12种岩石矿物为研究对象，建立了岩石矿物识别模型。彭伟航等（2019）基于改进InceptionV3模型对常见的16类矿物进行图像识别，准确率达到86%。李明超等（2020）基于迁移学习的方法与深度学习模型，根据19类矿物的颜色与纹理特征，建立了矿物图像识别耦合模型，并搭建矿物识别应用系统。刘艳鹏等（2020）利用卷积神经网络进行成矿预测。王李管等（2020）利用预训练的卷积神经网络对黑钨矿图像进行迁移学习实验，对黑钨原矿石的识别准确率达到99.6%。Liu et al.（2021）结合模型深度、模型结构和数据集大小，建立了小型矿石图像分类的深度学习模型。Zhou et al.（2022）在MobileNet网络的基础上，利用迁移学习与SENet机制相结合的方式，使矿石分类的准确率达到96%。通过上述文献不难看出，基于深度学习的矿物识别可以较好地解决传统方法中效率偏低的问题，并且准确率有所提高，但仍存在矿物识别种类较少，同时矿物形态特征发生变化时，已有的深度学习模型难以准确提取特征，以致于泛化能力不强的问题。
针对现有问题，本文以26类矿物的矿石图像作为研究对象，提出了一种基于改进的ConvNeXt网络的矿物智能识别模型，通过利用迁移学习策略，引入超强通道注意力机制，在加强矿物特征提取的同时，进一步提高特征之间的融合，以及矿物识别准确率与模型的泛化性，从而实现对不同矿物的高效准确识别。
1. 基本概念
在计算机视觉中，卷积神经网络凭借其对图像特征的快速提取与识别物体准确率高一直占据着重要地位。ConvNeXt网络是2022年由Facebook团队提出的纯卷积神经网络架构模型（Liu et al., 2022）。对于ConvNeXt网络，根据模型计算复杂度的不同，共有4个（T/S/B/L）版本。本文综合考虑矿物的类别与数据集的大小，选取了ConvNeXt-T版本。在ConvNeXt-T网络中，主要是由ConvNeXt块构成。对于ConvNeXt块，其输入特征图为h×w×dim（高度×宽度×特征图的通道深度），经过深度可分离卷积 (DepthwiseConv2D）及层标准化（Layer Normalization），对多个特征通道进行融合，之后经过普通卷积（Conv2D）等操作，对特征图进行一系列升维与降维处理，使其输出特征图为h×w×dim，具体变化过程如图1所示。对于整个ConvNeXt-T网络结构，首先输入大小为224×224×3（即高度与宽度大小为224，通道数为3）的图像，经过第一次卷积操作后特征图变为56×56×96；其次，经过一系列ConvNeXt块同时结合下采样（Downsample）操作，特征图变为7×7×768,最后，经全局平均池化（Global Avg Pooling）、层标准化与全连接（Linear）等操作后，特征图输出大小为1000。具体结构如表1所示。
表 1 ConvNeXt-T结构Table 1. Structure diagram of ConvNeXt-T层名 输入 ConvNeXt-T 输出 conv1 224×224×3 4×4,96,stride4
Layer Norm56×56×96 conv2_x 56×56×96 [d7×7,961×1,3841×1,96]×3 56×56×96 conv3_x 56×56×96 Downsample [d7×7,1921×1,7681×1,192]×3 28×28×192 conv4_x 28×28×192 Downsample[d7×7,3841×1,15361×1,384]×9 14×14×384 conv5_x 14×14×384 Downsample[d7×7,7681×1,30721×1,768]×3 7×7×768 7×7×768 Global Avg PoolingLayer NormalizationLinear 1000 2. 模型构建
2.1 迁移学习
理论上，深度学习可以提取图像中不同层次的特征，表征物体之间的细微差别。然而，这需要以足够多且高质量的已标注数据集为前提，进而训练得出已建立网络模型中的巨量参数。在实际应用场景中，研究对象（如矿物图像领域）往往没有建立大量带标签的数据集，通常会采用迁移学习的方式来解决（Zamir et al., 2018）。目的主要是实现模型之间的参数共享与相应的特征迁移，提高其泛化能力，同时减少训练成本。
考虑到研究的矿物对象种类及数据集大小，本文通过迁移学习的方式，实现ConvNeXt原模型与矿物识别模型之间参数共享，以降低其训练成本，提高矿物智能识别模型的泛化性。①利用原有的ConvNeXt网络对ImageNet-1k数据集进行预训练，得到ConvNeXt网络预训练模型与对应的模型权重文件；②通过微调的方式，删除上述预训练模型中原有的全连接层，并设计适合矿物研究对象的全连接层，从而适合本文研究的矿物识别领域；③基于迁移学习的方式，将以矿物种类的数量为标准而设计的全连接层与ConvNeXt模型相结合，以共享参数信息与特征权重，从而重新训练模型参数，加快训练速度。
2.2 超强通道注意力模块（ECA-Net）
现有研究表明，将注意力机制添加到已有的卷积神经网络框架或模型中可以提高其性能。比如，SENet通过对特征图中的通道维度进行信息权重分配，在ImageNet数据集上降低了0.4%~1.8%的错误率（Hu et al., 2018）；SKNet采用自适应的注意力编码机制，与ResNeXt50相比降低了1.44%的错误率（Li et al., 2019）；CBAM提出了通道与空间融合的注意力机制，较ResNet50降低了1.9%的错误率（Woo et al., 2018）。在复杂的环境变化下，矿物的特征发现多样化，本文在ConvNeXt网络基础上，在ConvNeXt块之后加入超强通道注意力机制ECA（Efficient Channel Attention）（Wang et al., 2020），加强了对矿物特征的提取。其主要流程为：首先，输入矿物图像特征图大小为C×H×W；在此基础上，进行GAP（全局平均池化）操作以减少模型参数量；其次，使用一维卷积（大小为k）获取矿物图像特征在每个空间通道的权重，其中，k表示局部通道交互的覆盖率，同时利用Sigmoid进行对应权重的归一化；最后，将已生成的各个通道权重与输入特征图相加权，进而得到矿物图像特征权重。超强通道注意力（ECA）模块如图2所示。
对于整个ECA模块来讲，主要是用Wk来表示学习到的通道注意力，其中，Wk涉及K×C个参数，并且Wk避免了不同通道的完全独立：
[w1,1⋯w1,k00⋯⋯00w2,2⋯w2,k+10⋯⋯0⋮⋮⋮⋮⋱⋮⋮⋮0⋯00⋯wC,C−k+1⋯wC,C] (1) 其中，ECA模块主要的部分是：利用卷积核大小为k的一维卷积实现每个通道之间的信息交互，从而生成每个通道对应的权重：
w=σ(C1Dk(y)) (2) 其中，CID为一维卷积，σ为Sigmoid函数。因此，ECA在保证了模型性能的同时，提高了模型效率（降低了模型复杂度）。
2.3 矿物智能识别模型
针对矿物图像识别困难及准确率较低的缺陷，本文在原有ConvNeXt网络的基础上，利用迁移学习策略与引入ECA模块，构建矿物智能识别网络模型，名称为ECA-ConvNeXt网络模型。具体来讲，主要是利用在ImageNet已经预训练的ConvNeXt模型，经过一定处理后，通过迁移学习加载到矿物识别模型框架内，而在现有ConvNeXt网络的基础上，为了实现矿物图像特征之间的融合，在每一个ConvNeXt块后（即在conv_2,3,4,5后），都添加一个ECA模块，实现了矿物图像局部跨通道交互融合的目的，最后经全局平均池化等操作，实现了对于矿物的识别。最终，本文构建的ECA-ConvNeXt网络模型，以不同种类矿物的矿石图像数据为研究对象，进行一系列实验后，矿物识别准确率及模型性能与现有主流网络相比都有较大提升。矿物智能识别模型整体结构如图3所示。
3. 实验过程
3.1 矿物图像数据采集
本文选取常见的26类矿物作为研究对象，包括斑铜矿、辰砂、赤铁矿、磁黄铁矿、磁铁矿、毒砂、方解石、方铅矿、橄榄石、铬铁矿、黑钨矿、褐铁矿、黄铁矿、黄铜矿、辉钼矿、辉锑矿、辉铜矿、孔雀石、蓝铜矿、铝土矿、软锰矿、闪锌矿、石英、铁钼矿、雄黄、萤石。
由于矿物种类不同，而且同类矿物产地来源也有所不同。不同的矿物经采集、运输等过程，不可避免地会沾染灰尘等，影响矿物图像的清晰度，以致于影响模型的精度。因此，首先对26类矿物的矿石进行水洗、晾晒等操作，将已经处理过的矿物的矿石样本放置于实验平台之上，然后使用佳能单反相机从不同的方向拍摄，拍摄距离约30 cm，其部分矿物的矿石样本图像如图版Ⅰ所示。
3.2 矿物图像数据增强
在深度学习中，是否有充足的数据样本是保证训练模型效果优劣与泛化能力强弱的基础。但是在实际应用中，往往没有建立统一的大型数据集。为了解决因数据集的不足造成识别准确率低等问题，通常采用数据增强的方式提升模型的效果。针对本文中不同矿物矿石图像的特点，对其进行图像裁剪、旋转，图像亮度调整和随机添加噪声操作。对于这些数据增强方法操作后，每类矿物的矿石图像原始数量和增强后的数量如表2所示，按照序列号进行排序，共计34576张图片。
表 2 矿物类别及图像数量统计Table 2. Mineral category and image quantity statistics序列号 矿物类别 原始图像数量/张 增强后的图像数量/张 序列号 矿物类别 原始图像数量/张 增强后的图像数量/张 1 斑铜矿 194 776 14 黄铜矿 360 1440 2 辰砂 214 856 15 辉钼矿 278 2224 3 赤铁矿 258 1032 16 辉锑矿 472 1888 4 磁黄铁矿 146 584 17 辉铜矿 300 1200 5 磁铁矿 286 1144 18 孔雀石 466 1864 6 毒砂 216 864 19 蓝铜矿 310 1240 7 方解石 374 1496 20 铝土矿 298 1192 8 方铅矿 516 2064 21 软锰矿 324 1296 9 橄榄石 142 568 22 闪锌矿 362 1448 10 铬铁矿 246 984 23 石英 257 1285 11 黑钨矿 324 1296 24 铁钼矿 121 968 12 褐铁矿 440 1760 25 雄黄 352 1408 13 黄铁矿 399 1995 26 萤石 426 1704 3.3 实验方法与环境配置
（1）实验设计
在对矿物识别模型实验之前，需将矿物图片按6∶2∶2的比例划分为矿物训练集、矿物验证集与矿物测试集。其中，矿物训练集中的图片用于训练本文构建模型中的未知参数，验证集目的是为了验证模型训练的结果，而其测试集用于检验矿物识别模型的效果。根据本文矿物数据集的大小，设置batch size大小为32，初始学习率为0.001，迭代200次，采用交叉熵函数作为模型的损失函数。
在矿物图片实验过程中使用Adam作为优化器，同时利用余弦退火对矿物图片训练过程中学习率进行不断优化，使其接近全局最优值点，进而得到矿物识别模型。
（2）环境配置
本实验基于Pytorch框架进行矿物识别模型的训练、验证与测试。其中，硬件环境：处理器为英特尔 i9-10900K（64GB），显卡为NVIDIA GeForce RTX 3080。软件环境：CUDA11.4，cuDNN8.2.4，Pytorch1.11，PyCharm2021专业版。
4. 结果分析
4.1 训练结果
本文利用训练准确率与验证准确率评估ECA-ConvNeXt网络矿物识别模型的训练效果。其中，训练准确率与验证准确率相差太大表示模型过拟合，导致网络泛化能力差。
为了验证本文模型的有效性，在一定条件下（数据集与处理方法相同，训练200次），分别在训练集与验证集上与5种典型图像分类网络进行横向的对比，其结果如图4所示。图4−a为训练集上不同模型准确率的对比，而图4−b为验证集上不同模型的准确率对比。在训练集上，6个网络的训练准确率均达到90%以上；在验证集上，VGG19准确率为92.22%，GoogLeNet准确率为92.43%，ResNet50准确率为95.04%，ResNeXt50准确率为95.43%，ConvNeXt准确率为97.53%，本文模型准确率为99.42%。从以上结果可知，本文构建的ECA-ConvNeXt网络模型与其他5种网络相比，取得了较优成绩，总体上性能最佳。
![]()
图 4 训练集（a）与验证集（b）不同网络的准确率对比Figure 4. Comparison of the accuracy of different networks in training sets (a) and validation sets (b)4.2 实验结果
4.2.1 模型性能
混淆矩阵是评判模型优劣的一种指标，常用于评判图像分类模型的优劣。本文基于改进的ConvNeXt矿物智能识别模型对于26类矿物的混淆矩阵如图5所示。其中，大多数准确率达97%及以上，对黄铁矿进行识别分类时只有94%，识别成毒砂与方解石有2%的概率，而识别成黄铜矿与闪锌矿有1%的概率。
对于图像分类任务，如何对图像所属类别进行正确分类是其主要的目的。在实际情况中，预测值与真实值的情况如表3所示。其中，TP为真实值为正，预测值为正；FP为真实值为正，预测值为负；FN为真实值为负，预测值为正；TN为真实值为负，预测值为负。
表 3 分类指标Table 3. Classification index预测值 正例（Positive) 反例（Negative) 真实值 正例（Positive) TP FP 反例（Negative) FN TN 为了验证本文构建的矿物智能识别模型的有效性，选用准确率（Accuracy）、精确率（Precision)与召回率（Recall）指标对6种矿物识别模型进行效果评估。计算公式如下所示：
Accuracy=TP+TNTP+TN+FP+FN (3) Precision=TPTP+FP (4) Recall=TPTP+FN (5) 在相同训练集上进行训练的6种矿物识别模型，在同一测试集上进行对应的测试，结果如表4所示。从表4可以看出，分类准确率与模型的复杂度呈正相关。其中，GoogLeNet优于VGG19，而ResNet50引入残差结构，ResNeXt50则是在ResNet的基础上结合分组卷积的优点，使其效果得到提升，而ConvNeXt是借鉴Swin-Transformer的结构，改进卷积神经网络的结构提升其性能。本文在ConvNeXt的基础上，引入注意力机制，使其准确率、精确率、召回率分别达到98.58%、98.62%、98.73%，与其他网络相比均是最佳。因此，对于ConvNeXt网络的改进较有效。
表 4 各个模型的测试评估结果Table 4. Test evaluation results of each model模型 准确率 精确率 召回率 VGG19 91.25% 91.62% 91.70% GoogLeNet 92.66% 93.36% 92.33% ResNet50 94.88% 95.04% 94.95% ResNeXt50 95.07% 95.45% 95.03% ConvNeXt 96.60% 96.62% 96.73% 本文模型 98.58% 98.62% 98.73% 4.2.2 消融实验分析
为了验证所提出的优化方法对于原ConvNeXt模型性能的影响，在相同实验条件下，采用不同优化方法对ConvNeXt模型进行测试，其结果如表5所示。在原ConvNeXt网络模型的基础上，利用ECA或迁移学习的方法，使其在准确率、精确率与召回率3项指标上都有1%左右的性能提升，而本文提出的改进ConvNeXt模型，只增加了极少的参数量与模型大小，约有2%的性能提升，取得了不错的效果。因此，本文提出的改进ConvNeXt模型更适合于对矿物的智能识别。
表 5 不同优化方法对于ConvNeXt模型的影响Table 5. The impact of different optimization methods on the ConvNeXt model模型 参数量 模型大小 准确率 精确率 召回率 ConvNeXt 27.80M 106.25M 96.60% 96.62% 96.73% ECA
+ConvNeXt27.82M 106.26M 97.62% 97.65% 97.61% 迁移学习+ConvNeXt 27.82M 106.26M 97.87% 97.90% 98.02% 本文模型 27.84M 106.27M 98.58% 98.62% 98.73% 4.3 矿物图像特征图可视化
卷积神经网络主要是通过提取图像中物体或其他本身所具有的特征，达到识别、检测、分割等目的，但对于其所提取的特征是什么，哪些区域对于识别起作用，以及由此得出的最终结果都无从得知。为了确定本文提出的矿物识别模型对于矿石图像中哪些区域识别的影响较大，本文利用Grad-CAM（Selvaraju et al., 2017）获得可视化矿石图像的特征图。主要是将本文提出的网络模型的最后一层feature maps（特征图），利用Grad-CAM进行可视化输出，可以观察矿物图像中的有效区域对于模型识别的影响。
图6是在不同模型的条件下，所提取特征的可视化热力图。其中，图6中从上到下的5种矿物分别是：斑铜矿、赤铁矿、橄榄石、辉钼矿和石英。而图6中从左到右依次是原图与依据VGG19、GoogLeNet、ResNet50、ResNeXt50、ConvNeXt和本文所提出的模型生成与之对应的特征图。其中，通过对于原图与不同模型之间特征图的对比（图中颜色越深，代表模型对于所在区域关注度越高，对识别的结果影响越大），本文提出的矿物识别模型优于传统模型，同时在对矿物识别的过程中，主要以矿物的边缘和颜色为感兴趣区域，也从侧面说明本文模型对于矿物特征的提取效果好，实现了对矿物的准确识别，证明了模型的有效性。因此，本文提出的基于改进的ConvNeXt网络矿物识别模型，提高了对于矿物的特征提取能力，较好地解决了矿物识别效率低的问题。
5. 结 论
在原有ConvNext网络的基础上，利用26类矿物34576张图像，通过迁移学习与引入ECA模块，构建了基于改进的ConvNeXt网络的矿物智能识别模型，实现了对于矿物分类的高效准确识别。得到了以下结论。
（1）通过在ImageNet上已预训练的ConvNeXt模型，经一定处理后利用迁移学习的方法，加载到矿物识别模型框架中，同时针对不同矿物图像自身特点，通过引入超强注意力通道（ECA）机制，构建了ECA-ConvNeXt网络矿物识别模型，提高了矿物图像特征融合的能力与训练速度。
（2）采用ECA-ConvNeXt矿物图像识别模型，将26类矿物图像与VGG19、GoogLeNet、ResNet50、ResNeXt50和ConvNeXt网络相比，其在训练过程收敛更快，同时在准确率、精确度与召回率3项数据上均是最佳，分别达到了98.58%、98.62%与98.73%；而消融实验，进一步证明了模型的有效性。
（3）利用Grad-CAM方法，在矿物原图的基础上，通过与不同模型之间矿物图像特征图的可视化结果对比，进一步证明了本文模型的优越性与在矿物识别领域的显著性。
本文研究成果对于矿物智能识别领域具有一定的借鉴意义，如何进一步提高矿物智能识别的速度与准确率，以及如何将其应用到实际场景等仍然是重要的研究方向。
-
![]()
图 1 松潘−甘孜地体及邻区大地构造单元简图(a, 据Yin et al., 2000)和吾和玛地区地质简图（b）
Figure 1. Tectonic units of the Songpan−Ganzi Terrane and its adjacent area (a) and simplified geological map of the Wuhema area (b)
![]()
图 2 吾和玛花岗质岩石代表性锆石阴极发光图像和LA-ICP-MS锆石U−Pb测年结果
Figure 2. Cathodoluminescence images of representative zircon grains and LA-ICP-MS zircon U−Pb dating results of the Wuhema granitoid rock
![]()
图 3 吾和玛花岗质岩石分类图解
a—A/CNK-A/NK图解（底图据Maniar et al.,1989）；b—TAS图解（底图据Middlemost,1994）；c—QAP图解；d—SiO2-K2O图解（底图据Rickwood et al.,1989）
Figure 3. Classification diagrams for the Wuhema granitoid rock
![]()
图 4 吾和玛花岗质岩石稀土元素球粒陨石标准化配分图(a)及微量元素原始地幔标准化蛛网图(b) (标准化值据McDonough et al.,1995)
Figure 4. Chondrite-normalized REE distribution patterns (a) and primative mantle normalized distribution patterns of trace elements(b) of the Wuhema granitoid rock
![]()
图 5 吾和玛花岗质岩石地球化学分类图解
a—(Zr+Nb+Ce+Y)−TFeO/ MgO图解；b—(Zr+Nb+Ce+Y)−(Na2O+K2O)/CaO图解；c，d—地球化学散点图。FG—分异的长英质花岗岩；OGT—未分异的 M、I、S 花岗岩 （a，b底图均据Whalen et al.,1987）
Figure 5. The geochemical classification diagrams of Wuhema granites rock
![]()
图 6 吾和玛花岗质岩石源区组成判别图（a底图据Stern et al., 1996，b ~ d底图据Altherr et al., 2002；Kaygusuz et al., 2008）
a—Mg#−SiO2图解；b—CaO/（MgO+TFeO）−K2O/Na2O图解；c—（Na2O+K2O+MgO+TFeO+TiO2）−（Na2O+K2O）/（MgO+TFeO+TiO2）图解；d—（CaO+MgO+TFeO+TiO2）−CaO/（MgO+TFeO+TiO2）图解.
Figure 6. Discrimination diagrams of magma source for Wuhema granites rock
![]()
图 7 吾和玛高分异I型花岗岩 Sr−Ba (a)和 La−(La/Yb)N (b)关系图及分离结晶趋势（矿物分离结晶趋势线据Rudnick et al., 1995）
Pl—斜长石；PlAn50—斜长石（An=50）；PlAn10—斜长石（An=10）；Kfs—钾长石；Amp—角闪石；Bi—黑云母；Zr—锆石；Sph—榍石；Ap—磷灰石；Mon—独居石；Allan—褐帘石
Figure 7. Sr−Ba (a) and La−(La/Yb)N (b) diagrams showing the fractional crystallization trends for the Wuhema highly fractionated I type granite
表 1 吾和玛花岗质岩石样品(P006-1-TW1)LA-ICP-MS锆石 U−Th−Pb分析结果
Table 1 LA-ICP-MS zircon U−Th−Pb isotopic analyses for the sample (P006-1-TW1) from Wuhema granitoid rock
测点号 同位素比值 年龄/Ma 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 208Pb/232Th 1σ 206Pb/238U 1σ 1 0.06784 0.00395 0.31801 0.01833 0.034 0.00066 0.01486 0.0006 215.6 4.1 2 0.05287 0.0019 0.26214 0.00959 0.03596 0.00058 0.01137 0.00031 227.7 3.58 3 0.05818 0.00189 0.29154 0.00969 0.03635 0.00058 0.01273 0.00029 230.2 3.58 4 0.05088 0.0023 0.22936 0.01042 0.0327 0.00055 0.00826 0.00035 207.4 3.46 5 0.05048 0.00162 0.24334 0.00801 0.03497 0.00055 0.00969 0.00022 221.5 3.4 6 0.05288 0.0022 0.24881 0.01043 0.03413 0.00057 0.01042 0.00028 216.3 3.53 7 0.05986 0.00209 0.27428 0.00971 0.03323 0.00054 0.01053 0.00028 210.8 3.36 8 0.05236 0.00302 0.24352 0.01396 0.03374 0.00062 0.01108 0.00039 213.9 3.86 9 0.04904 0.00201 0.22609 0.00933 0.03344 0.00055 0.00955 0.00027 212.1 3.42 10 0.05268 0.00206 0.24899 0.00984 0.03428 0.00056 0.01051 0.0003 217.3 3.49 11 0.06386 0.0023 0.27834 0.01013 0.03161 0.00052 0.01295 0.00033 200.6 3.23 12 0.09867 0.00235 0.48867 0.01222 0.03592 0.00056 0.01552 0.00029 227.5 3.48 13 0.05043 0.00351 0.2396 0.01657 0.03446 0.00066 0.0113 0.00036 218.4 4.11 14 0.05113 0.0023 0.24263 0.01099 0.03442 0.00058 0.01022 0.0003 218.2 3.59 15 0.06856 0.00256 0.2753 0.01036 0.02913 0.00048 0.01268 0.00035 185.1 3.02 16 0.05385 0.00278 0.25171 0.01297 0.03391 0.0006 0.01174 0.00037 215 3.73 17 0.05591 0.0023 0.24805 0.0103 0.03218 0.00053 0.01085 0.00032 204.2 3.31 18 0.05553 0.00258 0.2759 0.01287 0.03604 0.00062 0.01222 0.00046 228.2 3.84 19 0.05093 0.002 0.2454 0.00975 0.03495 0.00057 0.00982 0.00031 221.5 3.56 20 0.33927 0.00537 3.34562 0.05896 0.07153 0.00111 0.18423 0.00285 445.4 6.67 
下载: 导出CSV
表 2 吾和玛花岗质岩石锆石Lu-Hf同位素分析结果
Table 2 Zircon Lu-Hf isotopic data of Wuhema granitoid rock
测点号 176Yb/177Hf 176Lu/177Hf ±2σ (176Hf/177Hf )i εHf(t) ±1σ TDM2/Ga Hf-BD02 0.023024 0.000850 0.000018 0.282764 -0.27 0.03 1.10 Hf-BD03 0.029089 0.001110 0.000022 0.282736 -1.28 0.03 1.17 Hf-BD04 0.022269 0.000861 0.000022 0.282727 -1.59 0.03 1.19 Hf-BD05 0.037337 0.001418 0.000022 0.282807 1.24 0.03 1.01 Hf-BD06 0.021790 0.000849 0.000019 0.282700 -2.54 0.03 1.25 Hf-BD07 0.036377 0.001300 0.000021 0.282809 1.32 0.03 1.00 Hf-BD08 0.037459 0.001345 0.000022 0.282765 -0.26 0.03 1.10 Hf-BD09 0.024019 0.000934 0.000021 0.282734 -1.35 0.03 1.17 Hf-BD10 0.022585 0.000857 0.000021 0.282711 -2.15 0.03 1.22 Hf-BD11 0.032601 0.001179 0.000024 0.282753 -0.68 0.03 1.13 注: εHf(t)={［(176Hf/177Hf)样品-(176Lu/177Hf)样品×(eλt－1)］/［(176Hf/177Hf)CHUR－(176Lu/177Hf)CHUR×(eλt–1)］－1)}×10000，tDM1=1/λ×ln{1 +［(176Hf/177Hf)样品-(176Lu/177Hf)DM］/［(176Lu/177Hf)样品-(176Lu/177Hf)DM］}，TDM2 = TDM1-( TDM1－t)×［(fcc–f样品)/(fcc–fDM)］，(176Lu/177Hf) CHUR =0. 0332，(176Hf/177Hf) CHUR= 0.282772(据Blichert-Toft et al.,1999) ； (176Lu/177Hf) DM = 0. 0384, (176Hf/177Hf)DM = 0.28325 (据Griffin et al.,2000) , fcc, f样品和fDM分别代表陆壳、样品和亏损地幔的fLu/Hf值，t表示锆石结晶年龄，λ = 1. 867 × 10−11/a
下载: 导出CSV
表 3 吾和玛花岗质岩石全岩主量、微量和稀土元素组成分析结果
Table 3 Whole-rock major, trace and rare earth elements of the Wuhema granitoid rock
样品编号 PM006-2-YH2 PM006-2-YH1 PM006-3-YH1 PM006-3-YH2 PM006-4-YH1 PM006-4-YH2 PM006-5-YH1 PM006-5-YH2 SiO2 75.5 74.0 76.8 76.8 76.4 77.3 76.5 76.7 TiO2 0.075 0.22 0.10 0.10 0.08 0.06 0.11 0.07 Al2O3 12.26 13.4 12.22 12.31 12.23 12.27 12.31 12.43 Fe2O3 1.50 0.34 0.14 0.05 0.23 0.09 0.05 0.03 FeO 1.30 1.40 1.4 1.30 1.4 0.80 1.50 1.10 MnO 0.02 0.04 0.03 0.03 0.03 0.02 0.04 0.04 MgO 0.16 0.25 0.12 0.11 0.11 0.06 0.15 0.07 CaO 0.65 1.2 0.80 0.70 0.70 0.50 0.70 0.70 Na2O 3.49 3.5 3.41 3.14 3.13 3.34 3.4 3.36 K2O 4.43 4.32 4.44 4.81 4.82 4.88 4.59 4.82 P2O5 0.012 0.06 0.02 0.02 0.02 0.01 0.02 0.02 H2O 4.53 4.32 4.78 4.8 4.75 4.91 4.74 4.78 总计 99.99 100 100 100.01 100 99.98 100.01 99.98 A/CNK 1.05 1.07 1.03 1.06 1.06 1.05 1.04 1.04 烧失量 0.60 1.25 0.58 0.65 0.95 0.69 0.59 0.70 Mg# 10 21 12 13 11 11 15 10 DI 92 93 93 93 93 95 90 94 Ta 1.56 0.85 0.9 0.81 0.85 0.93 1.16 0.99 U 10.33 5.38 7.08 8.73 10.9 7.78 8.53 8.81 Ba 27.3 89 100 78.3 63.4 54.6 66.2 37.9 Cs 12.74 14.6 11.6 17.7 19.1 14.5 19.8 14.2 Zr 77.95 135 88.6 104 85.2 88 95.8 76.3 Hf 3.34 4.82 3.76 4.44 4.01 4.12 4.48 3.53 Rb 317.82 258 312 348 397 418 384 429 Th 33.63 28.1 30.7 28 31.9 21.2 35.2 31.3 Sr 26.47 29 39.2 33 28.4 21.9 32.8 21.6 Sn 3.72 7.34 6.12 14.7 9.8 6.04 10.4 9.54 Nb 7.84 12.6 12.4 12.5 12.6 13 14.2 13 La 13.43 23.4 15.4 7.66 16.2 7.12 17 12.5 Ce 74.91 46.2 34.7 19.2 38.6 15.5 37.6 30.2 Pr 4.03 5.26 4.16 2.25 4.51 1.61 4.42 3.66 Nd 16.08 19.8 16.5 8.49 17.4 5.85 17 13.9 Sm 4.27 4.42 4.55 2.18 4.46 1.49 4.87 3.44 Eu 0.16 0.61 0.21 0.17 0.12 0.067 0.17 0.1 Gd 4.39 4.11 4.56 2.61 4.41 2.12 5.12 3.57 Tb 0.73 0.65 0.76 0.48 0.75 0.4 0.9 0.58 Dy 4.93 3.89 4.73 3.3 4.82 2.75 5.75 3.54 Ho 1.05 0.83 0.98 0.73 0.99 0.61 1.24 0.76 Er 3.33 2.43 2.84 2.3 3.01 1.9 3.66 2.27 Tm 0.58 0.38 0.44 0.36 0.48 0.31 0.58 0.36 Yb 4.25 2.61 2.94 2.56 3.26 2.16 4.12 2.58 Lu 0.72 0.37 0.43 0.36 0.48 0.3 0.57 0.36 Y 31.41 22.4 26.4 20.8 28 18.4 34.5 21.5 δEu 0.11 0.43 0.14 0.22 0.08 0.12 0.1 0.09 Sr/Y 0.05 0.04 0.05 0.05 0.05 0.05 0.05 0.05 (La/Yb)N 2.27 6.43 3.76 2.15 3.56 2.36 2.96 3.48 La/Yb 3.16 8.97 5.24 2.99 4.97 3.3 4.13 4.84 ΣREE 132.86 114.96 93.2 52.65 99.49 42.19 103 77.82 注：主量元素含量单位为%，微量和稀土元素含量单位为10−6；Mg#=100×Mg2+/(Mg2++TFe2+)；A/CNK=Al2O3/(CaO+Na2O+K2O)；A/NK=Al2O3/(Na2O+K2O)；δEu=2×EuN/(SmN+GdN)；N表示球粒陨石标准化
下载: 导出CSV
-
Altherr R, Siebel W. 2002. I−type plutonism in a continental back−arc setting: Miocene granitoids and monzonites from the central Aegean Sea, Greece[J]. Contributions to Mineralogy and Petrology, 143(4): 397−415. doi: 10.1007/s00410-002-0352-y
Blichert−Toft J, Albarède F. 1999. The Lu−Hf isotope geochemistry of chondrites and the evolution of the mantle−crust system[J]. Earth and Planetary Science Letters, 148: 243−258.
Cai H M, Zhang H F, Xu W C, et al. 2010. Petrogenesis of indosinian volcanic rocks in Songpan−Garze fold belt of the northeastern Tibetan Plateau: New evidence for lithospheric delamination[J]. Science China (Earth Sciences ), 53(9): 1316−1328. doi: 10.1007/s11430-010-4033-9
Chappell B W. 1999. Aluminium saturation in I− and S−type granites and the characterization of fractionated haplogranites[J]. Lithos, 46: 535−551. doi: 10.1016/S0024-4937(98)00086-3
Griffin W L, Pearson N J, Belousova E, et al. 2000. The Hf isotope composition of cratonic mantle: LA−MC−ICPMS analysis of zircon megacrysts in kimberlites[J]. Geochimica et Cosmochimica Acta, 64(1): 133−147. doi: 10.1016/S0016-7037(99)00343-9
Kaygusuz A, Siebel W, En C, et al. 2008. Petrochemistry and petrology of I−type granitoids in an arc setting: the composite Torul pluton, Eastern Pontides, NE Turkey[J]. International Journal of Earth Sciences, 97(4): 739−764. doi: 10.1007/s00531-007-0188-9
King P L, Chappell B W, Allen C M et al. 2001. Are A−typegranites the high−temperature felsic granites? Evidence from fractionated granites of the Wangrah suite[J]. Australian Journal of Earth Sciences, 48(4): 501−514. doi: 10.1046/j.1440-0952.2001.00881.x
Li X H, Li Z X, Li W X, et al. 2007. U−Pb zircon, geochemical and Sr−Nd−Hf isotopic constraints on age and origin of Jurassic I− and A−type granites from central Guangdong, SE China: a major igneous event in response to foundering of a subducted flat−slab?[J]. Lithos, 96: 186−204.
Liu L, Chu G, Li Y G, et al. 2018. Mesozoic high− and low−SiO2 adakites and A−type granites in the Lower Yangtze River Belt, Eastern China: Implications for petrogenesis and metallogeny[J]. Minerals, 8(8): 328. doi: 10.3390/min8080328
Ludwig K R. 2003. ISOPLOT 3.00: A Geochronological Toolkit for Microsoft Excel [M]. Berkeley Geochronology Center, California.
Maniar P D, Piccoli P M. 1989. Tectonic discrimination of granitoids[J]. GSA Bulletin, 101(5): 635−643. doi: 10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
McDonough W F, Sun S S. 1995. The composition of the Earth[J]. Chemical Geology, 120(3/4): 223−253.
Middlemost E A K. 1994. Naming Materials in the Magma / igneous Rock System[J]. Earth−Science Reviews, 37(3/4): 215−224.
Rickwood P C. 1989. Boundary lines within petrologic diagrams which use oxides of major and minor elements[J]. Lithos, 22(4): 247−263. doi: 10.1016/0024-4937(89)90028-5
Roger F, Malavieille J, Leloup P H, et al. 2004. Timing of granite emplacement and cooling in the Songpan−Garze Fold Belt (eastern Tibetan Plateau) with tectonic implications[J]. Journal of Asian Earth Sciences, 22: 465−481. doi: 10.1016/S1367-9120(03)00089-0
Rudnick R L. 1995. Making continental crust[J]. Nature, 378(6557): 571−577. doi: 10.1038/378571a0
Stern C R, Kilian R. 1996. Role of the subducted slab, mantle wedge and continental crust in the generation of adakites from the Andean Austral Volcanic Zone[J]. Contributions to Mineralogy and Petrology, 123(3): 263−281. doi: 10.1007/s004100050155
Van Achterbergh E, Ryan C G, Jackson S E, et al. 2001. Data reduction software for LA−ICP−MS[C]//Sylvester P J. Laser−Ablation−ICPMS in the Earth Sciences. Principles and Applications. Mineralogical Society of Canada Short Course Series, 29: 239–243.
Wang Q, Li Z X, Chung S L, et al. 2011. Late Triassic high−Mg andesite/dacite suites from northern Hohxil, North Tibet: Geochronology, geochemical characteristics, petrogenetic processes and tectonic implications[J]. Lithos, 126(1/2): 54−67.
Watson E B, Harrison T M. 1983. Zircon saturation revisited: temperature and composition effects in a variety of crustal magma types[J]. Earth and Planetary Science Letters, 64: 295−304. doi: 10.1016/0012-821X(83)90211-X
Whalen J B, Currie K L, Chappell B W. 1987. A−type granites: geochemical characteristics, discrimination and petrogenesis[J]. Contributions to Mineralogy and Petrology, 95: 407−419. doi: 10.1007/BF00402202
Wu F Y, Liu X C, Ji W Q, et al. 2017. Highly fractionated granites: recognition and research[J]. Science China Earth Sciences, 60(7): 1201−1219. doi: 10.1007/s11430-016-5139-1
Xiao L, Zhang H F, Clemens J D, et al. 2007. Late Triassic granitoids of the eastern margin of the Tibetan Plateau: Geochronology, petrogenesis and implications for tectonic evolution[J]. Lithos, 96(3/4): 436−452.
Yan D P, Zhou M F, Song H L, et al. 2003. Origin and tectonic significance of a Mesozoic multi−layer over−thrust system within the Yangtze Block (Yangtze Block)[J]. Tectonophysics, 361: 239−254. doi: 10.1016/S0040-1951(02)00646-7
Yin A, Harrison T M. 2000. Geologic evolution of the Himalayan Tibetan orogen[J]. Annual Review of Earth and Planetary Sciences, 28: 211−280. doi: 10.1146/annurev.earth.28.1.211
Yuan C, Zhou M F, Sun M, et al. 2010. Triassic granitoids in the eastern Songpan−Ganzi Fold Belt, SW China: Magmatic response to geodynamics of the deep lithosphere[J]. Earth and Planetary Science Letters, 290(3/4): 481−492. doi: 10.1016/j.jpgl.2010.01.005
Zhang H F, Parrish R, Zhang L, et al. 2007. A−type granite and adakitic magmatism association in Songpan−Garze fold belt, eastern Tibetan Plateau: Implication for lithospheric delamination[J]. Lithos, 97(3/4): 323−335.
Zhang H F, Zhang L, Harris N, et al. 2006. U−Pb zircon ages, geochemical and isotopic compositions of granitoids in Songpan−Garze fold belt, eastern Tibetan Plateau: Constraints on petrogenesis and tectonic evolution of the basement[J]. Contributions to Mineralogy and Petrology, 152(1): 75−88. doi: 10.1007/s00410-006-0095-2
Zhang L Y, Ding L, Pullen A, et al. 2014. Age and geochemistry of western Hoh−Xil−Songpan−Ganzi granitoids, northern Tibet: Implications for the Mesozoic closure of the Paleo−Tethys ocean[J]. Lithos, 190: 328−348.
边飞, 吴柏林, 高永旺. 2013. 青海扎日加花岗岩地球化学、锆石 LA−ICP−MS U−Pb定年及地质意义[J]. 矿床地质, 32(3): 625−640. 崔加伟, 郑有业, 孙祥. 2016. 松潘−甘孜造山带北段岗龙乡花岗岩岩石成因及地球动力学特征[J]. 地质科技情报, 35(2): 129−139. 国家地质实验测试中心. 2010. (GB/T 14506.28−2010)硅酸盐岩石化学分析方法 第28部分: 16个主次成分量测定[S]. 北京:中国标准出版社. 胡健民, 孟庆任, 石玉若, 等. 2005. 松潘−甘孜构造带内花岗岩锆石SHRIMP U−Pb定年及其构造意义[J]. 岩石学报, 21(3): 867−880. 李献华, 李武显, 李正祥. 2007. 再论南论燕山早期花岗岩的成因类型与构造意义[J]. 科学通报, 52(9): 981−992. 李艳广, 汪双双, 刘民武, 等. 2015. 斜锆石LA−ICP−MS U−Pb定年方法及应用[J]. 地质学报, 89(12): 2400− 2418. 吴福元, 李献华, 郑永飞, 等. 2007. Lu−Hf 同位素体系及其岩石学应用[J]. 岩石学报, 23(2): 185−220. 夏磊, 闫全人, 向忠金, 等. 2017. 松潘−甘孜地体中部晚三叠世安山质增生弧的确定及其意义[J]. 岩石学报, 33(2): 579−604. 许志琴, 侯立炜, 王宗秀. 1992. 中国松潘带的造山过程 [M]. 北京: 地质出版社. 张雪亭. 2007. 1∶100万青海省地质图及说明书[M]. 北京:地质出版社. 赵永久. 2007. 松潘−甘孜东部中生代中酸性侵入体的地球化学特征、岩石成因及构造意义[D]. 中国科学院研究生院(广州地球化学研究所)博士学位论文. 朱弟成, 莫宣学, 王立全, 等. 2009. 西藏冈底斯东部察隅高分异I型花岗岩的成因: 锆石U−Pb年代学、地球化学和Sr−Nd−Hf同位素约束[J]. 中国科学( D辑), 39(7): 833−848. -
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 2778
- HTML全文浏览量: 1788
- PDF下载量: 250
- 被引次数: 1
