Zircon U-Pb geochronology and trace elements geochemistry of Hiltaba Suite in the north of Moonta-Wallaroo mining district, Yorke Peninsula, Gawler Craton, Australia
-
摘要:
澳大利亚南澳州高勒克拉通广泛发育的Hiltaba岩套与该地区的铁氧化物铜金矿床关系密切，查明该岩套的成岩时代、成岩温度及构造背景，对铁氧化物铜金矿床的研究至关重要。对约克半岛Moonta-Wallaroo矿集区北部的眼球状花岗岩进行锆石U-Pb测年，测得年龄为1589±43 Ma，与前人确定的Hiltaba岩套岩浆事件（1595~1575 Ma）时间一致。通过锆石Ti温度计计算出花岗岩成岩温度为848~971℃，平均878℃，接近A型花岗岩的成岩温度。锆石的稀土元素配分模式呈现轻稀土元素亏损、重稀土元素富集及向左倾斜的特征，具有明显的正Ce异常和负Eu异常，δCe值变化范围较大（4.85~107.63），δEu值变化范围较小（0.17~0.36）。在锆石Y-Yb/Dy等判别图解中，样品多落入板内构造环境范围及附近，少数落入火山弧构造环境范围及附近；在锆石Hf-U/Yb等判别图解中，样品均落入陆壳环境范围。综合分析认为，该花岗岩形成于活动大陆边缘向陆内环境演变的构造背景。
Abstract:The Hiltaba Suite is widely developed in the Gawler Craton, South Australia, which is closely related to the IOCG deposits in the area. It is very important to find out the diagenetic age, diagenetic temperature and tectonic setting of the suite for the study of IOCG deposits. Zircon U-Pb dating of the augen granite collected from the north of the Moonta-Wallaroo mining district in the York Peninsula shows that the age is 1589 ±43 Ma, which is consistent with the magmatic event (1595~1575 Ma) of the Hiltaba Suite determined by previous geologists. According to zircon Ti thermometer, the diagenetic temperature of granite is 848~971℃, with an average of 878℃, which is close to that of A-type granite. The rare earth element patterns of zircon show the characteristics of LREE loss, HREE enrichment and left-dip, with obvious positive Ce anomaly and negative Eu anomaly, δCe varies widely (4.85~107.63) and δEu varies small (0.17~0.36). In the zircon Y-Yb/Dy and other discriminant diagrams, most of the samples fall within and near the within plate tectonic environment, and a few fall within and near the volcanic arc tectonic environment. In the zircon Hf-U/Yb and Y-Yb/Dy discriminant diagrams, the samples fall within the continental crust environment. Comprehensive analysis shows that the granite was formed in the tectonic setting of the evolution from the active continental margin to the intracontinental environment.
-
Keywords:
- Gawler Craton /
- IOCG deposits /
- Hiltaba Suite /
- diagenetic age /
- diagenetic temperature /
- tectonic setting
-
进入新世纪以来，中国面临城镇化进程的跨越。党的十九大报告提出实施区域协调发展战略，要求“以城市群为主体构建大中小城市和小城镇协调发展的城镇格局”，新型城镇化建设进入新时期。城市不断建设和扩张，建设用地、工业用地增长的同时，侵蚀着耕地、水域、草地等。安庆市是快速发展的长三角区域的城市之一，多年来城市地区人类活动剧烈，原始地貌改变众多，地表覆被严重，给城市勘察带来众多困难。其中一个典型问题为，原有的地表河浜可能由于各种原因被填埋形成暗浜。暗浜底部多以淤泥为主，部分地段还含有大量生活垃圾，土性极差。随着时间的推移，这些埋藏于地下的不良地质条件逐渐成为后期市政道路、桥梁工程建设中的质量隐患，造成工程开挖难度加大，尤其给河网地区浅表的工程建设带来较多的困难。因此，提取城市地区的暗浜分布范围，为城市工程建设与城市规划提供决策支持，具有十分重要的意义。
暗浜反映出地表覆被类型的变化。目前暗浜的探测多以物探方法为主[1-2]，如二维微动剖面探测法、瞬态瑞雷波法、高密度电法, 这些探测方法一般是在对工作区地质条件较熟悉的条件下开展的，且适合小范围的勘察，成本较高。随着卫星对地观测技术的发展，遥感影像的时空分辨率在不断提高，这为监测和刻画土地利用变化时空过程提供了丰富、可靠的数据支撑。变化检测是通过观察一个对象或一个现象在不同时间的状态，来确定其发生的差异过程[3-6]。遥感影像变化检测是利用不同时期覆盖同一地表区域的多源遥感影像和相关地理空间数据，结合相应地物特性和遥感成像机理，采用图像、图形处理理论及数理模型方法，确定和分析该地域地物的变化，包括地物位置、范围的变化和地物性质、状态的变化[7]。其研究的目的是找出感兴趣的变化信息，滤除作为干扰因素的不相干的变化信息。
多时相高分辨率遥感影像具有空间信息丰富、地物细节清晰的特点, 可以利用空间信息提高变化检测精度，已越来越多地应用于城市监测等领域[8-12]。面向对象的处理是高分辨率遥感影像解译领域的重要思想。面向对象的变化检测不仅可以利用基于像素的方法所使用的光谱、纹理和特征变换等作为变化的测度，还可以利用有关对象的形状特征和空间关系的附加信息。面向对象的遥感图像变化检测首先要解决的问题是如何根据多时相影像获取适合变化检测的图像对象，以图像分割技术为基础，基于获取的图像对象进行变化检测。Desclé等[13]对SPOT多期影像进行分割，计算每期影像的光谱值，利用统计方法提取出变化的对象。Conchedda等[14]为检测红树林的变化，利用分割方法中的区域生长法多时相影像进行分割。Stow等[15]应用变化检测技术提取植被的变化信息，其分割方法选择试错法确定分割参数。Duveiller等[16]将多时相遥感影像叠加，同时分割，用以提取森林覆盖的变化。Wang等[17]在此分割方法基础上，对多时相图像对象提取了多种特征，选择最优的阈值界定“变”与“不变”，提取土地覆盖类型的变化信息。很多学者利用基于对象的分类后比较方法对地图或GIS图层进行更新[18-20]，能够提取出“from-to”信息，在土地覆盖/利用变化检测中应用广泛。就遥感影像分类方法而言，SVM的独特优势在于可以处理小样本、高维数、非线性数据问题，将SVM算法应用于遥感影像分类时，可以通过少量的训练样本获得比传统分类方法更高的分类精度。由于SVM分类算法应用广泛且分类精度较高，为进行暗浜自动化监测提供了有效途径，故本文选择SVM算法为遥感影像分类方法。
暗浜反映出地面以下地质体与周围物体结构的差异，如水塘填埋物与其周围土体结构必然是不同的。而传统的光学遥感卫星难以检测地下土体或岩石结构特征。微动测深是一种便捷有效的物探手段，基于微动信号估算地层速度结构，抗干扰能力强，探测深度大，适用于人口密度大、有振动干扰的城区房屋密集区环境[21-24]，广泛应用于地铁工程、道路工程、城市地质等工程领域。为了验证利用遥感变化检测提取的暗浜范围是否准确可靠，将遥感技术与微动探测方法相结合，基于暗浜提取结果，选取典型区域，采用微动探测进行实地勘察。
本文以提取安庆研究区(主要分布在主城区)的暗浜分布范围为目的，利用面向对象的分类后比较的变化监测方法，分别对2期遥感影像T1与T2进行面向对象的监督分类，进而对分类结果进行叠加分析，提取出前期土地利用类型为水体但目前变为其他土地覆盖类型的区域，即为暗浜的分布区域。在此基础上，将检测出来的变化区域作为靶区，选择典型区域，利用物探工作方法中的微动测量方法，探测地下是否存在暗浜，验证变化检测结果的可靠性。本次工作旨在利用遥感影像自动识别暗浜方法和技术流程，为城市工程建设与城市规划提供决策支持。
1. 研究区与数据介绍
安庆市位于安徽省西南部，长江下游北岸，皖河入江处，河网密集，是一个快速发展的城市。多年来安庆市采取了“东进、北扩、西拓”的扩张方式，对于城市的改造力度很大，由此导致了城市地表覆被类型的变化。在城市规划和工程建设中，暗浜是一个较典型的不良工程地质现象，对工程施工、尤其对基础工程建设构成隐患。为避免暗浜对工程建设尤其是基坑工程建设带来不利影响，有必要查明暗浜的特征和空间分布范围，为实施地球物理勘探方法圈定靶区。
暗浜调查区域布置为矩阵区域，以安庆市辖区为主，包括宜秀区、十里铺、老峰镇、大龙山镇等，数据选取2004年的Quickbird影像与2018年的worldview-2影像。调查区域与数据选择见图 1。
![]()
图 1 研究区概况及遥感影像数据a—研究区位置图；b—T2(2018年)影像；c—T1(2004年)影像Figure 1. Illustration of the study area and remote sensing images2. 研究方法
暗浜提取采用以遥感图像识别为主，微动探测为辅的思路，首先利用遥感图像变化检测方法提取暗浜范围，进而结合微动探测技术，选择典型区域，验证结果的可靠性。遥感图像变化检测技术通过处理相同区域在不同时相的遥感影像，得到该区域的地物类型变化，已经在城市建设、森林保护、土地监测、灾害评估等领域发挥了重要的作用。在变化检测中，分类后比较法是一种根据分类比较方法得到变化细节的判别方法，目前使用非常广泛。分类后比较具体的做法是分别对不同时相的遥感影像进行分类，得到不同时期的专题图，通过叠加这些专题图像，判别随着时间跃迁，类别的具体变化。微动探测方法可以勘察地表以下地层结构，根据横波波速结构或视横波波速剖面图，结合钻孔钻探资料及其他地质资料加以解释，给出地质现象解译，判断暗浜是否存在。
技术路线图见图 2，首先针对2期高分辨率遥感影像T1和T2，对数据预处理，进而采用分类后比较的变化检测方法，其中分类时以图像分割得到的影像对象为基本处理单元，不仅利用遥感影像的光谱信息，而且充分利用影像的形状、纹理、上下文等信息，提高分类精度。根据分类结果，提取暗浜范围，进而选择典型的区域，辅以微动探测，布设试验，验证暗浜提取结果的可靠性。
2.1 图像预处理
影像预处理是面向对象的变化检测的前提，主要包括辐射校正和几何校正。首先对2期影像进行相对辐射校正，以降低因图像辐射差异对变化检测结果的影响。然后对影像进行几何校正，使2期影像精确配准，控制误差小于0.5像素。
2.2 图像分割
由于高分辨率遥感影像单个地物特征内反射率的变异性及图像的复杂性，传统的逐像素分析方法难以适应于高分影像的应用需求。研究表明，基于对象的影像分析(OBIA)技术可减少地理参考及光谱变异性等带来的影响。图像分割是OBIA的核心，分割过程将图像划分为光谱上相似且在空间上相邻的同质对象。理想的分割结果是，对象内部的方差最小，对象与对象之间的异质性较大。通过分割生成图像对象，即分割体，获得更丰富的信息，包括纹理、形状和与相邻对象的空间关系，从而利用空间上下文信息进一步分类。
目前的分割方法可以分为2类：①基于边缘的影像分割；②基于区域的影像分割。基于边界的分割算法是按照影像中各类地物像素灰度值不连续的特点、边界特征差异性的原则实现影像分割[19]。目前常用的边缘检测算子有Krisch算子、Robert算子等。基于边缘的方法更符合人们的认知，当影像各地物类型之间有显著差异时，该方法能够取得良好的效果。基于区域分割的方法主要分为2种：①区域增长；②区域分裂与合并。区域分割的依据是影像内像素点或像素块的特征相似性，包括灰度、颜色、纹理、空间结构等特征信息[25]。
影像分割的核心问题在于如何选择最优的影像分割尺度，在最优分割尺度下获得的影像对象能够更好地表达实际地物的特征信息，有利于变化检测结果精度的提高[26]。本文影像分割步骤采用ENVI 5.3的Feature Extraction模块。该模块采用基于边缘的分割方法，可以根据遥感影像内邻近像素的色度、亮度、纹理等信息对影像进行分割，同时通过分割效果预览框可以实时地进行分割尺度的调整，操作简单方便，应用性较强。本文反复试验对比，通过目视判断确定研究区的最优分割尺度。
2.3 特征提取
利用影像的原始特征光谱，提取指数特征及纹理特征作为特征组合。指数特征有归一化水体指数NDWI及归一化植被指数NDVI。
NDWI=ρgreen−ρnirρgreen+ρnir (1) (1) 式中，ρnir与ρgreen分别表示影像近红外波段、绿色波段的反射率，突出水体信息。
NDVI=ρnir−ρredρnir+ρred (2) (2) 式中，ρred表示影像红色波段的反射率，增强植被信息。
纹理特征使用灰度共生矩阵计算，通过8种基于二阶矩阵的纹理滤波的结果对比，选择均值滤波的结果作为纹理特征。利用PCA提取特征组合。通过协方差矩阵计算主成分，选择包含95%的信息量的前3个特征分量。
2.4 图像分类
基于分割结果，利用支持向量机(SVM) 算法，对图像进行监督分类。经过对研究区影像进行判读分析，确定建设用地(包括城镇用地、农村居民点、其他建设用地)、林地、草地、耕地、水域、未利用土地等类别，选择各类地物的训练样本，训练样本的位置随机选择。在此基础上，利用SVM算法对影像进行分类。
SVM算法的基本思想是，基于训练样本集D在样本空间找到具有“最大间隔”的划分超平面，将不同类别的样本分开。SVM在线性分类的基础上引入了结构风险最小化原理，提高了算法的泛化能力，在小样本的情况下也能达到较好的分类精度。SVM与核函数方法结合，解决了非线性分类问题。由于此算法的上述优势，在遥感影像分类中广泛应用。
具体模型如下：
min (3) 其中, w=(w1; w2; …; wd)，为超平面的法向量，b为位移项，xi∈D为训练样本集中的一个样本。
对式(3)使用拉格朗日乘子可得到其“对偶问题”：
\begin{gathered} \max \limits_{\alpha} \sum\limits_{i=1}^{m} \alpha_{i}-\frac{1}{2} \sum\limits_{i=1}^{m} \sum\limits_{j=1}^{m} \alpha_{i} \alpha_{j} y_{i} y_{j} x_{i}^{T} x_{j} \\ { s.t. }\ \sum\limits_{i=1}^{m} \alpha_{i} y_{j}=0 , \\ \alpha_{i} \geqslant 0, i=1, 2, \cdots, m \end{gathered} (4) 对于线性不可分问题，将核技巧应用到支持向量机，基本思想是通过非线性变换将输入空间对应一个特征空间，使得在输入空间Rn中的超曲面模型对应于特征空间H中的超平面模型。这样，分类问题的学习任务通过在特征空间中求解线性支持向量机就可以完成。此时，对偶问题则为：
\begin{aligned} &\max \limits_{\alpha} \sum\limits_{i=1}^{m} \alpha_{i}-\frac{1}{2} \sum\limits_{i=1}^{m} \sum\limits_{j=1}^{m} \alpha_{i} \alpha_{j} y_{i} y_{j} K\left(x_{i}, x_{j}\right) \\ & { s.t. }\ \sum\limits_{i=1}^{m} \alpha_{i} y_{j}=0, 0 \leq \alpha_{i} \leq C, i=1, 2, \cdots, m \end{aligned} (5) 其中, K(xi, xj)为核函数，C为惩罚参数。常用的SVM核函数有多项式核函数，高斯核函数，sigmoid核函数等。已有研究表明，高斯核函数更适用于土地覆被的分类[27]，因而采用高斯核函数作为SVM分类器中的核函数，其数学形式为：
K(x, z)=exp \left(-\frac{|| x-z||^{2}}{2 \sigma^{2}}\right) (6) 2.5 变化检测
变化检测的流程如图 3所示。变化检测是基于SVM算法得到地物分类结果，将2期图像的结果进行对比，提取时相T1(2004年)影像上的水体在时相T2(2018年)影像的变化情况。首先将2期影像的分类结果图转为分类矢量图，对2个图层进行叠加分析，得到包含变化信息的多边形，同时赋予这些多边形类别变化的属性信息，完成变化检测，并统计各类别的增加、减少情况。
2.6 地球物理特征分析
基于变化检测提取的暗浜范围，选择前期是坑塘目前已变为耕地或道路的区域作为采样点，布设微动测线，根据现场的探测结果解译，判断地下是否有存在暗浜的可能。图 4为在安庆市迎江区布设的3条微动测线(底图 4，为2004年遥感影像)，可知在2004年此区域地表覆盖类型为池塘，而现场测量时发现，此处已在开发建设，坑塘早已被填埋。若暗浜填充物的波阻抗大于周围地层介质波阻抗，则微动速度-深度剖面上暗浜区域应表现为相对低速异常，反之则为相对高速异常。由此，推断地质断面图并判断暗浜的位置。
![]()
图 4 微动测量工作布置示意图(2004影像)a—测线1和测线2；b—测线3Figure 4. Layout of microtremor survey (images of 2004)3. 结果
3.1 变化检测实验结果
对2期影像分别进行裁剪、拼接、相对辐射校正、地理配准，得到T1与T2两期影像。本文目的在于提取暗浜的空间分布，即以前地表为水体包括沟塘河渠的地区，现在变化为非水体。因此，实验中，对T1影像即2004年QuickBird影像分为水体和非水体2类，将T2影像即2018年WorldView-2影像分为建设用地(包括城镇用地、农村居民点、其他建设用地)、林地、草地、耕地、水域、未利用土地等类别，分类结果如图 5所示。
![]()
图 5 影像分类结果a—研究区位置图；b—T2影像分类结果；c—T1影像水体提取结果Figure 5. The classification results of remote sensing images最后对分类结果进行精度评价，发现Kappa系数很高，结果准确度较高(表 1)。
表 1 分类精度Table 1. Accuracy of classification影像类型 总体精度 Kappa系数 T2影像 90.49% 0.8506 T1影像 96.12% 0.94 对2期影像分类结果进行叠加分析，得到T1影像表现为水体而T2影像表现为非水体的图斑，作为暗浜空间范围(图 6)。
对暗浜图斑的属性进行赋值，分析得到研究区暗浜幼前期的水体目前转变为何种地物类别(图 7)。
3.2 微动探测试验结果
基于变化检测结果，选取2个区域，布设微动探测剖面3条，共计140个点，采用直线型布阵方式进行测量，台阵半径1 m-2 m-4 m，采样频率250 Hz，测量时间20 min，点距8 m。采用标准化流程进行数据处理，并根据收集的钻探成果对其过程参数进行修正，在总体与钻探成果较为吻合后确定最终的微动速度-深度剖面图，并推断地质断面图，如图 8所示。根据图 8-a并结合遥感图像可知，1~5号点为深色暗浜或水田区域，11~60号点为暗浜区域，61~71号点为深色暗浜或水田区域。图 8-b为测线1推断地质断面图，10 m以浅的低速异常推测为暗浜分布区域，速度范围在100~300 m/s之间，变化检测结果基本一致。同时，按剖面图 8-c中100~300 m/s速度范围区分暗浜，发现15 m以浅的低速异常为暗浜分布区域，由于该微动测线布设靠近暗浜中心，故而低速异常深度相比图 8-a稍大。
![]()
图 8 不同测线微动速度-深度剖面和推断地质断面图(a、c、e分别为测线1、测线2、测线3微动速度-深度剖面图，b、d、f分别为测线1、测线2、测线3推断地质断面图)Velocity-depth profile of microtremor survey and inferred geographical profile of different line由图 8-e分析，发现在14号点、22号点、26号点存在局部的低速异常，根据变化检测结果可知，在2004年时15~22号点都属于暗浜范围(图 8-f)，至2018年地物地貌变化极大，暗浜基本被回填。
4. 讨论
利用本文的算法，提取出研究区暗浜面积9.86 km2，从空间分布看，大部分暗浜位于安庆市辖区南部沿江地区，包括宜秀区、十里铺乡、老峰镇、白泽胡乡、山口乡，暗浜表面主要地物类型为建筑物(表 2)。整体看，有89.15%的暗浜区域，其地表覆盖类型目前为建设用地，如道路、桥梁、房屋等，其中居民点占42.35%，其他建设用地占46.80%。这跟安庆市近十五年的发展建设一致，城市的扩张使得不透水层面积急剧增加，以前的很多小池塘沟渠被填埋，用以城市开发建设。此外，暗浜区域地表覆盖类型为耕地的面积为0.69 km2，占暗浜总面积的6.97%，一些湖泊被人为造田，使水域面积减少。
表 2 暗浜空间范围分布面积Table 2. Area statistics of underground siltkm2 名称 建设用地 耕地 林地 草地 未利用土地 暗浜总计 宜秀区 2.3368 0.0486 0.0053 0.1844 — 2.5751 十里铺 2.0097 0.0064 0.0510 0.0008 — 2.0679 老峰镇 1.6146 0.3560 0.0010 0.0050 0.0093 1.9859 白泽湖 0.5841 0.0622 — — — 0.6463 山口乡 0.5852 — 0.0186 0.0010 — 0.6048 长风乡 0.5028 0.0607 0.0016 0.0067 — 0.5718 龙狮乡 0.5120 0.0496 — — — 0.5616 大龙山 0.3972 0.0021 0.0055 0.0159 — 0.4207 大观区 0.0932 — 0.0007 0.0000 — 0.0939 月山镇 0.0716 — 0.0072 0.0035 0.0007 0.0830 海口镇 0.0003 0.0659 0.0000 0.0000 — 0.0662 迎江区 0.0184 0.0268 0.0045 0.0138 — 0.0635 菱北街 0.0306 — — 0.0010 — 0.0316 杨桥镇 0.0049 — 0.0180 0.0024 — 0.0253 鲟鱼镇 0.0225 — 0.0009 0.0000 — 0.0234 五横乡 — — 0.0202 0.0015 — 0.0217 新洲乡 0.0090 0.0091 — 0.0018 — 0.0199 总计 8.7929 0.6874 0.1345 0.2378 0.0100 9.8626 基于变化检测结果及微动探测实地验证，推断地质断面，能够勘察得到暗浜的深度和边界。经实地验证，发现布设的3处测线均存在暗浜，填充物较复杂，填充不均匀，且2处由于地表城市规划已开始施工作业，导致微动剖面速度特征杂乱；在填充物简单、地表未进行人类工程活动时，微动探测应用效果更好。
5. 结论
遥感技术的发展为地质调查提供了极大的便利。本文将遥感与物探手段结合，首先利用遥感变化检测技术监测暗浜分布范围，为实施物探划定出靶区或重点区域，这样既能极大程度节约成本，又能提高暗浜检测的准确性。首先从遥感数据入手，利用变化检测方法自动提取安庆市城区暗浜的空间位置与范围，共提取暗浜面积9.86 km2，且发现约89%的区域目前已为建设用地，且建设用地中超过一半的区域为道路、桥梁等工矿用地，留有工程隐患的可能性很大。然而，利用遥感影像提取到的暗浜是二维的空间分布形态，无法得到暗浜位置准确的深度信息。本文基于变化检测技术监测暗浜分布范围，布设3条测线进行了实地勘察，验证了变化检测结果的可靠性。本文的暗浜提取方法与结果能够为城市工程建设与规划提供决策支持，具有十分重要的意义。
致谢: 本次样品采集和数据处理工作由中国地质调查局南京地质调查中心和澳大利亚南澳州地质调查局共同完成，南澳地调局安东尼教授和侯宝宏教授给予了大力支持，审稿专家在成文过程中提出了诸多宝贵意见，在此一并表示衷心的感谢。 -
![]()
图 1 高勒克拉通地质矿产简图(据参考文献[16]修改)
Figure 1. Geological and mineral map of the Gawler Craton
![]()
图 2 约克半岛Moonta-Wallaroo矿集区北部地区地质简图(据参考文献[5]修改)
Figure 2. Geological map of the region in the north of Moonta-Wallaroo District, Yorke Peninsula
![]()
图 4 钻孔PBD30岩石特征
a—眼球状花岗岩岩心照片；b—中细粒花岗岩岩心照片；c—样品岩心照片；d—样品镜下照片，正交偏光。Kfs—钾长石；Qtz—石英；Bt—黑云母
Figure 4. Characteristics of rocks from borehole PBD30
![]()
图 5 眼球状花岗岩中代表性锆石阴极发光图像、分析点位置及年龄值
Figure 5. Cathodoluminesence images, location of U-Pb spot analyses and ages of representative zircon
![]()
图 6 眼球状花岗岩中代表性锆石球粒陨石标准化稀土元素配分模式(球粒陨石标准值据参考文献[26])
Figure 6. Chondrite-normalized rare earth element patterns of representative zircons in the augen granite
![]()
图 7 眼球状花岗岩中代表性锆石La-(Sm/La)N判别图(a)和(Sm/La)N-Ce/Ce*判别图(b)(底图据参考文献[30])
Figure 7. La-(Sm/La)N (a) and (Sm/La)N-Ce/Ce* (b) diagrams of representative zircons in the augen granite
![]()
图 8 眼球状花岗岩中代表性锆石U-Pb谐和图(a) 和207Pb/206Pb年龄谱图(b、c)
Figure 8. U-Pb concordia diagram (a) and 207Pb/206Pb age spectrums (b, c) of representative zircons in the augen granite
![]()
图 9 不同构造背景下锆石微量元素判别图(底图据参考文献[40])
a—Y-Yb/Dy判别图；b—Y-Ce判别图；c—Cd-Th判别图；d—Er-U判别图；e—Nb-Ta判别图；f—Y-Lu/Hf判别图；g—Y-Th/U判别图
Figure 9. Trace element discriminant diagrams of zircons in different tectonic setting
![]()
图 10 陆壳和洋壳环境中锆石Hf-U/Yb (a) 与Y-U/Yb (b) 判别图(底图据参考文献[41])
Figure 10. Hf-U/Yb (a) and Y-U/Yb (b) discriminant diagrams of zircon in continental and ocean crust
表 1 眼球状花岗岩中代表性锆石LA-ICP-MS U-Th-Pb测试结果
Table 1 LA-ICP-MS U-Th-Pb data of representative zircons in the augen granite
点号 含量/10-6 Th/U 同位素比值 年龄/Ma Th U 207Pb/206Pb 1 σ 207Pb/235U 1 σ 206Pb/238U 1 σ 207Pb/206Pb 1 σ 207Pb/235U 1 σ 206Pb/238U 1 σ 1 342 621 0.55 0.0966 0.0023 3.3944 0.0826 0.2516 0.0033 1559 44 1503 19 1447 17 2 307 488 0.63 0.0979 0.0022 3.7373 0.0832 0.2739 0.0033 1584 43 1579 18 1560 17 3 451 974 0.46 0.0953 0.0020 3.3399 0.0694 0.2509 0.0029 1544 40 1490 16 1443 15 4 131 150 0.88 0.0970 0.0025 4.0733 0.1018 0.3018 0.0041 1566 48 1649 20 1700 20 5 239 490 0.49 0.0979 0.0023 3.6270 0.0793 0.2654 0.0031 1585 43 1555 17 1517 16 6 81 96 0.84 0.0959 0.0028 3.9585 0.1151 0.2967 0.0049 1546 49 1626 24 1675 24 7 454 535 0.85 0.0986 0.0028 3.7424 0.0969 0.2721 0.0037 1598 54 1580 21 1552 19 8 466 742 0.63 0.0976 0.0026 3.5643 0.1082 0.2611 0.0057 1589 51 1542 24 1495 29 9 401 881 0.46 0.0992 0.0028 3.7575 0.0941 0.2715 0.0038 1610 47 1584 20 1549 19 10 254 400 0.63 0.0981 0.0025 4.2550 0.1022 0.3117 0.0046 1588 48 1685 20 1749 22 11 279 429 0.65 0.0978 0.0023 4.1256 0.0965 0.3024 0.0041 1583 44 1659 19 1703 20 12 478 979 0.49 0.0971 0.0022 3.2652 0.0769 0.2411 0.0033 1569 42 1473 18 1392 17 13 366 607 0.60 0.0980 0.0021 3.7823 0.0774 0.2774 0.0030 1587 41 1589 16 1578 15 14 307 682 0.45 0.0980 0.0023 3.4856 0.0956 0.2558 0.0051 1587 43 1524 22 1469 26 15 373 570 0.65 0.0986 0.0023 4.0140 0.0912 0.2930 0.0033 1598 43 1637 19 1656 17 
下载: 导出CSV
表 2 眼球状花岗岩中代表性锆石微量元素含量
Table 2 Trace element compositions of representative zircons in the augen granite
点号 La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Y Hf Nb Ta ΣREE ΣLREE ΣHREE δEu δCe 1 1.84 61.3 2.83 20.7 33.9 8.00 113 33.4 311 90.9 351 68.2 627 113 2676 12228 8.98 5.26 1837 129 1708 0.36 5.33 2 0.97 51.6 0.89 6.36 8.05 1.77 32.7 11.8 128 45.0 196 42.2 395 76.8 1300 11434 7.87 4.32 997 70 927 0.29 12.58 3 1.77 85.6 2.17 11.3 11.8 2.57 35.6 12.5 133 45.8 209 46.9 463 92.4 1361 12850 10.7 5.59 1154 115 1039 0.35 9.22 4 0.68 29.4 0.29 3.10 5.64 1.06 24.4 8.09 92.8 34.9 153 31.4 287 56.0 996 10298 3.66 1.69 727 40 687 0.24 16.35 5 0.83 49.5 1.45 12.5 17.8 4.74 72.7 21.0 186 54.3 206 41.3 383 73.9 1619 12009 7.74 4.03 1125 87 1039 0.35 8.61 6 2.01 23.1 0.67 5.77 7.18 1.60 28.3 8.73 95.4 34.5 149 29.8 267 52.5 996 9958 3.39 1.67 705 40 665 0.30 4.85 7 1.66 64.2 2.12 19.1 27.6 7.18 114 33.4 314 92.0 341 64.2 550 102 2732 10079 6.93 3.29 1732 122 1610 0.34 7.14 8 0.88 104 1.08 7.16 8.26 1.54 29.4 10.6 126 46.4 217 47.4 462 93.2 1374 11539 10.3 5.43 1154 123 1032 0.27 22.50 9 0.85 56.0 1.09 8.45 12.2 3.06 44.8 14.1 140 44.1 190 41.4 410 83.0 1332 13084 9.24 5.37 1050 82 968 0.36 12.13 10 0.02 37.4 0.10 1.48 2.95 0.48 15.6 5.61 69.1 27.6 134 30.5 301 62.8 845 11482 6.93 4.11 688 42 646 0.17 107.63 11 0.07 47.9 0.23 2.71 4.26 0.98 23.1 7.77 89.9 34.3 160 35.1 339 69.3 1026 11477 7.20 3.90 814 56 758 0.24 58.59 12 2.10 108 1.64 8.47 7.76 1.83 27.7 9.43 109 39.9 192 44.0 445 94.1 1246 12913 10.1 6.35 1090 130 960 0.34 13.54 13 0.22 52.3 0.37 2.98 4.26 0.86 19.7 7.27 86.1 32.7 156 35.7 353 74.3 985 11777 8.36 4.64 825 61 764 0.24 35.87 14 0.99 52.7 1.52 11.1 17.8 4.42 72.2 21.0 201 60.9 257 58.2 545 105 1821 12174 11.3 6.25 1408 89 1319 0.33 8.55 15 0.15 43.9 0.30 2.75 4.54 0.86 23.8 8.49 99.1 39.2 181 39.2 374 76.3 1158 11651 9.11 4.71 894 53 842 0.20 38.41 注：微量元素含量单位为10-6
下载: 导出CSV
表 3 锆石Ti温度计计算结果
Table 3 Calculation results of zircon Ti thermometer
点号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 T/℃ 870 867 859 927 859 971 899 872 848 874 875 858 865 858 866
下载: 导出CSV
-
Perring C S, Pollard P J, Dong G, et al. The Lingtning Creek sill complex, Cloncurry district, Northwest Queensland: A source of fluids for Fe oxide Cu-Au mineralization and sodic-calcic alteration[J]. Economic Geology, 2000, 95(5): 1067-1089. doi: 10.2113/gsecongeo.95.5.1067
Pollard P J. Evidence of a magmatic fluid and metal source for Fe-oxide Cu-Aumineralization[C]//Porter T M. Hydrothermal iron oxide copper-gold and related deposits: A global perspective. Adelaide: Australian Mineral Foundation, 2000: 27-41.
Pollard P J. An intrusion-related origin for Cu-Au mineralization in iron oxide-copper-gold(IOCG) provinces[J]. Mineralium Deposita, 2006, 41(2): 179-187. doi: 10.1007/s00126-006-0054-x
Sillitoe R H. Iron oxide-copper-gold deposits: an Andean view[J]. Mineralium Deposita, 2003, 38(7): 787-812. doi: 10.1007/s00126-003-0379-7
毛景文, 余金杰, 袁顺达, 等. 铁氧化物-铜-金(IOCG)型矿床: 基本特征、研究现状与找矿勘查[J]. 矿床地质, 2008, 27(3): 267-278. doi: 10.3969/j.issn.0258-7106.2008.03.001 聂凤军, 江思宏, 路彦明. 氧化铁型铜-金(IOCG)型矿床的地质特征、成因机理与找矿模型[J]. 中国地质, 2008, 35(6): 1074-1087. doi: 10.3969/j.issn.1000-3657.2008.06.005 赵宇浩, 毛大华, 王天刚, 等. 智利Candelaria-Punta del Cobre铁氧化物铜金矿床地质和成矿作用简[J]. 地质通报, 2017, 36(12): 2296-2307. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20171221&flag=1 Reid A J. The Olympic Cu-Au Province, Gawler Craton: A Review of the Lithospheric Architecture, Geodynamic Setting, Alteration Systems, Cover Successions and Prospectivity[J]. Minerals, 2019, 9: 1-37. http://www.researchgate.net/publication/333908454_The_Olympic_Cu-Au_Province_Gawler_Craton_A_Review_of_the_Lithospheric_Architecture_Geodynamic_Setting_Alteration_Systems_Cover_Successions_and_Prospectivity/download
Creaser R A. Petrogenesis of a Mesoproterozoic quartz latite-granitoid suite from the Roxby Downs area, South Australia[J]. Precambrian Research, 1996, 79(3): 371-394. http://www.sciencedirect.com/science/article/pii/S0301926896000022
Budd A R. A- and I-type subdivision of the Gawler Ranges-Hiltaba Volcano-Plutonic Association[J]. Geochimica et Cosmochimica Acta, 2006, 70(18-supp-S): A72. http://www.onacademic.com/detail/journal_1000035388115810_638f.html
Chapman N, Ferguson M, Meffre S, et al. Pb-isotopic constraints on the source of A-type Suites: Insights from the Hiltaba Suite-Gawler Range Volcanics Magmatic Event, Gawler Craton, South Australia[J]. Lithos, 2019, 346/347: 1-18. http://www.sciencedirect.com/science/article/pii/S002449371930307X
Brotodewo A, Tiddy C J, Reid A, et al. Relationships between magmatism and deformation in northern Yorke Peninsula and southeastern Proterozoic Australia[J]. Australian Journal of Earth Sciences, 2018, 65: 619-641. doi: 10.1080/08120099.2018.1470573
姚仲友, 王天刚, 王国平, 等. 大洋洲地区优势矿产资源潜力评价[M]. 科学出版社, 2015: 99-107. 姚仲友, 王天刚, 傅朝义, 等. 大洋洲地区大地构造格架与优势矿产资源[J]. 地质通报, 2014, 33(2/3): 143-158. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=2014020301&flag=1 Reid A J, Fabris A. Influence of pre-existing low metamorphic grade sedimentary sucessions on the distribution of iron oxide-copper-gold mineralization in the Olympic Cu-Au province, Gawler craton[J]. Economic Geology, 2015, 110: 2147-2157 doi: 10.2113/econgeo.110.8.2147
丁绍磊. 澳大利亚南澳省Hillside与Olympic Dam铁氧化物铜-金矿床(IOCG)对比研究[D]. 成都理工大学硕士学位论文, 2016. Reid A J. Geology and metallogeny of the Gawler Craton[C]//Phillips G N. Australian Ore Deposits. Melbourne: The Australian Institute of Mining and Metallrugy, 2017: 589-594.
Allen S R, McPhie J, Ferris G, et al. Evolution and architecture of a large felsic igneous province in western Laurentia: the 1.6 Ga Gawler Range Volcanics, South Australia[J]. Journal of Volcanology and Geothermal Research, 2008, 172: 132-147. doi: 10.1016/j.jvolgeores.2005.09.027
Haynes D W, Cross K C, Bills R T, et al. Olympic Dam ore genesis; a fluid-mixing model[J]. Economic Geology, 1995, 90(2): 281-307. doi: 10.2113/gsecongeo.90.2.281
Belperio A, Flint R, Freeman H. Prominent Hill-a hematite-dominated, iron-oxide copper-gold system[J]. Economic Geology, 2007, 102(8): 1499-1510. doi: 10.2113/gsecongeo.102.8.1499
McPhie J, Kamenetsky V S, Chambefort I, et al. Origin of the supergiant Olympic Dam Cu-U-Au-Ag deposit, South Australia: was a sedimentary basin involved?[J]. Geology, 2011, 39(8): 795-798. http://www.researchgate.net/profile/Vadim_Kamenetsky/publication/263891547_Origin_of_the_supergiant_Olympic_Dam_Cu-U-Au-Ag_deposit_South_Australia_Was_a_sedimentary_basin_involved/links/551a22a80cf244e9a45854a7.pdf
Ferris G M, Schwarz M P. Proterozoic gold province of the Central Gawler Craton[J]. MESA Journal, 2003, 30: 4-12. http://www.researchgate.net/publication/284033530_Proterozoic_gold_province_of_the_central_Gawler_Craton
Mortimer G E, Cooper J A, Oliver R L. The geochemical evolution of Proterozoic granitoids near Port Lincoln in the Gawler orogenic domain of South Australia[J]. Precambrian Research, 1988, 40/41: 387-406. doi: 10.1016/0301-9268(88)90077-0
Reid A, Hand M, Jagodzinski E, et al. Paleoproterozoic orogenesis in the southeastern Gawler craton, South Australia[J]. Australian Journal of Earth Sciences, 2008, 55: 449-471. doi: 10.1080/08120090801888594
Slama J, Kosler J, Condon D J, et al. Plešovice zircon-A new natural reference material for U-Pb and Hf isotopic microanalysis[J]. Chemical Geology, 2008, 249(1/2): 1-35.
Sun S S, McDonough W F. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes[C]//Saunders A D, Norry M J. Magmatism in the Ocean Basins. The Geological Society, London, Special Publications, 1989, 42(1): 313-345.
Dubińska E, Bylinab P, Kozowskia A, et al. U-Pb dating of serpentinization: Hydrothermal zircon from ametasomatic rodingite shell(Sudetic ophiolite, SW Poland)[J]. Chemical Geology, 2004, 203(3/4): 183-203.
Schaltegger U. Hydrothermal zircon[J]. Elements, 2007, 3(1): 51-79. doi: 10.2113/gselements.3.1.51
赵振华. 副矿物微量元素地球化学特征在成岩成矿作用研究中的应用[J]. 地学前缘, 2010, 17(1): 267-286. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201001027.htm Hoskin P W O. Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the Jack Hills, Australia[J]. Geochimica et Cosmochimica Acta, 2005, 69(3): 637-648. doi: 10.1016/j.gca.2004.07.006
范飞鹏, 陈乐柱, 李海立, 等. 南岭东段枫树洞稀土矿中辉长质包体锆石U-Pb年代学, 地球化学特征及其成岩作用[J]. 华东地质, 2020, 41(4): 325-338. https://www.cnki.com.cn/Article/CJFDTOTAL-HSDZ202004003.htm Watson E B, Harrison T M. Zircon thermometer reveals minimum melting conditions on earliest Earth[J]. Science, 2005, 308(5723): 841-844. doi: 10.1126/science.1110873
Watson E B, Wark D A, Thomas J B. Crystallization thermometers for zircon and rutile[J]. Contributions to Mineralogy & Petrology, 2006, 151(4): 413-433. doi: 10.1007%2Fs00410-006-0068-5.pdf
Ferry J M, Watson E B. New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers[J]. Contributions to Mineralogy & Petrology, 2007, 154(4): 429-437. http://www.researchgate.net/profile/Edward_Watson2/publication/216832207_New_thermodynamic_models_and_revised_calibrations_for_the_Ti-in-zircon_and_Zr-in-rutile_thermometers/links/54ff52ad0cf2eaf210b91ff2.pdf
高晓英, 郑永飞. 金红石Zr和锆石Ti含量地质温度计[J]. 岩石学报, 2011, 27(2): 417-431. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201102006.htm Hayden L A, Watson E B. Rutile saturation in hydrous siliceous melts and its bearing on Ti-thermometry of quartz and zircon[J]. Earth and Planetary Science Letters, 2007, 258(3/4): 561-568. http://www.onacademic.com/detail/journal_1000035380473610_d5fe.html
King P L, White A J R, Chappell B W, et al. Characterization and Origin of Aluminous A-type Granites from the Lachlan Fold Belt, Southeastern Australia[J]. Journal of Petrology, 1997, 38(3): 371-391. doi: 10.1093/petroj/38.3.371
刘昌实, 陈小明, 陈培荣, 等. A型岩套的分类、判别标志和成因[J]. 高校地质学报, 2003, 9(4): 573-591. doi: 10.3969/j.issn.1006-7493.2003.04.011 蔡杨, 吴维平, 马涛, 等. 安徽九华山柯村花岗斑岩锆石U-Pb年代学、地球化学特征及地质意义[J]. 华东地质, 2018, 39(3): 177-186. https://www.cnki.com.cn/Article/CJFDTOTAL-HSDZ201803003.htm 陈超, 滕学建, 潘志龙, 等. 内蒙古北山造山带中段石板井地区A型花岗岩锆石U-Pb年龄对北山洋闭合时间的限定[J]. 地质通报, 2020, 39(9): 1448-1460. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20200910&flag=1 Schulz B, Klemd R, Bratz H. Host rock compositional controls on zircon trace element signatures in metabasites from the Austroalpine basement[J]. Geochimica et Cosmochimica Acta, 2006, 70(3): 697-710. doi: 10.1016/j.gca.2005.10.001
Grimes C B, John B E, Kelemen P B, et al. 2007. Trace element chemistry of zircons from oceanic crust: A method for distinguishing detrital zircon provenance[J]. Geology, 2007, 35(7): 643-646. doi: 10.1130/G23603A.1
Budd A R, Wyborn L A I, Bastrakova I V. The Metallogenic Potential of Australian Proterozoic Granites. 2001.
Stewart K P, Foden J. Mesoproterozoic granites of South Australia. 2003.
Coles B, Tassell S W, Street G J, et al. Progress reports, annual reports and final report to licence expiry/full surrender. 1991.
Conor C H H. Moonta-Wallaroo region: an interpretation of the geology of the Maitland and Wallaroo 1: 100 000 sheet areas. 1995.
Department of Primary industries and Resources SA. WHYALLA 1: 250, 000 geological map. 2006.
-
期刊类型引用(11)
1. 刘娟娟，林一杰，张新英. 农田土水稻系统中硒与镉的相互关系研究. 农业与技术. 2025(01): 1-5 . 
百度学术
2. 谭卓贤，杜建军，孙星，易琼，徐培智，张木. 石灰、磷酸盐及硅酸盐对土壤硒有效性及水稻累积硒的影响. 江苏农业学报. 2024(03): 450-456 . 百度学术
3. 刘瑞，周卫军，彭素华，商贵铎，周雨舟，李敏. 湖南省石门县耕地土壤硒含量特征及影响因素解析. 土壤通报. 2023(04): 840-847 . 百度学术
4. Xiu-jin Liu，Ke Yang，Fei Guo，Shi-qi Tang，Ying-han Liu，Li Zhang，Hang-xin Cheng，Fei Liu. Effects and mechanism of igneous rock on selenium in the tropical soil-rice system in Hainan Province, South China. China Geology. 2022(01): 1-11 . 必应学术
5. 王蕊，陈楠，张二喜. 基于总量与形态的矿区周边土壤重金属生态风险与健康风险评估. 环境科学. 2022(03): 1546-1557 . 百度学术
6. 吴超，孙彬彬，成晓梦，周国华，贺灵，曾道明，梁倍源. 丘陵山区多目标区域地球化学调查不同成因表层土壤代表性研究——以浙江绍兴地区为例. 地质通报. 2022(09): 1539-1549 . 本站查看
7. 袁宏伟，陈江均，郭腾达，吴艳君，杨敏. 巴彦淖尔市临河区狼山镇和新华镇一带富硒土壤地球化学特征及影响因素. 地质与勘探. 2022(05): 1027-1041 . 百度学术
8. 刘冰权，沙珉，谢长瑜，周强强，魏星星，周梵. 江西赣县清溪地区土壤硒地球化学特征和水稻根系土硒生物有效性影响因素. 岩矿测试. 2021(05): 740-750 . 百度学术
9. 杨泽，刘国栋，戴慧敏，张一鹤，肖红叶，阴雨超. 黑龙江兴凯湖平原土壤硒地球化学特征及富硒土地开发潜力. 地质通报. 2021(10): 1773-1782 . 本站查看
10. 郭军，刘明，汤恒佳，廖琦，邢新丽. 湖南汨罗市范家园地区富硒土壤特征及其影响因素研究. 华南地质. 2021(04): 387-397 . 百度学术
11. 于炎炎. 安徽砀山县土壤锌地球化学特征及影响因素. 矿产与地质. 2021(06): 1147-1155+1170 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 2681
- HTML全文浏览量: 528
- PDF下载量: 1796
- 被引次数: 13
