Zircon U-Pb dating, geochemistry and geological significance of the Duobagou granite pluton on the southern margin of the Dunhuang landmass
-
摘要:
多坝沟黑云二长花岗岩体出露于敦煌地块南缘，呈近东西向带状分布，侵位于太古宙—古元古代敦煌岩群。LA-ICPMS锆石U-Pb年龄显示，该岩体形成于晚二叠世，206Pb/238U年龄加权平均值为252.1±2.3Ma。该岩体具高SiO2含量(68.32%~72.82%)和全碱含量(7.76%~8.00%)；A/CNK值为0.97~1.02，均小于1.1，为高钾钙碱性准铝质花岗岩。岩石具有明显的轻、重稀土元素分异和弱的负Eu异常，具高Sr/Yb和(La/Yb)N值富集大离子亲石元素(K、Rb、Ba、Cs)，显著亏损Nb、Ta、P和Ti，表明源岩具有经典岛弧岩石地球化学特征；(87Sr/86Sr)i值为0.707083，εNd(t)为-2.22，Nd模式年龄TDM2=1102Ma。元素和同位素地球化学特征显示，多坝沟花岗岩形成于加厚下地壳中岛弧玄武岩的部分熔融。结合岩体地质特征、年代学和区域地质背景，初步认为敦煌地块在晚古生代晚期依然受古亚洲洋汇聚碰撞相关的造山活动的影响，可能是中亚造山带的最南缘组成部分。
-
关键词:
- 花岗岩 /
- LA-ICP-MS锆石U-Pb年龄 /
- 二叠纪 /
- 敦煌地块 /
- 中亚造山带
Abstract:The Duobagou biotite monzonitic granitic pluton, distributed as a nearly EW-trending belt, intruded into Archean-Paleoproterozoic Dunhuang complexes in southern Dunhuang landmass. Zircon LA ICP-MS U-Pb data revealed the crystallization age of 252.1±2.3Ma, suggesting that Duobagou pluton was formed in Late Permian. Geochemical characteristics of monzonitic granite are similar to those of typical arc rocks. For example, the granites have high content of SiO2 (68.32%~72.82%) and total alkali (7.76%~8.00%), and their A/CNK ratios are 0.97 to 1.02, which are less than 1.1, indicating that Duobagou granites are high potassium calcium alkaline aluminium. They have obvious differentiation of REEs and weak Eu negative anomaly, and are characterized by strong enrichment of K, Rb, Ba and Cs and depletion of Nb, Ta, P, and Ti, with high ratios of Sr/Yb and (La/Yb)N. The sample PM04-1 shows high initial Sr isotopic ratio (0.707083), relatively low εNd(t) (-2.22) and Nd model ages (TDM2=1102Ma). The authors thus hold that the Duobagou biotite monzonitic granites may be the product of partial melting of island arc basalt in thickened low crust, and the Dunhuang landmass was likely to be involved into a series of orogenic events by the closing of the Palaeo-Asian Ocean at the end of Late Paleozoic, and possibly became the southernmost part of central Asian orogenic belt.
-
Keywords:
- granites /
- LA-ICP-MS zircon U-Pb age /
- Permian /
- Dunhuang landmass /
- Central Asian orogenic belt
-
地质数据是地质工作中形成的重要信息资源，多年累积形成的海量地质数据信息为国民经济和社会发展提供了重要支撑[1]。地质行业历史悠久，地质资料积累丰厚[2]，因此，实现海量地质数据存储具有重要的现实意义。中国的地质数据库建设工作自20世纪80年代中后期起步，时间上晚于西方国家。目前为止，建立的地质数据库主要包括油气及固体矿产资源勘探、地质环境调查、水文地质、矿产资源潜力评价等内容，数量也已经具有一定的规模[3]，采取的方法主要为关系数据库。张博[4]以ArcGIS为基础信息系统，通过ArcCatalog自带功能，建立基于Personal Geodatabase的陕西府谷县地质灾害空间数据库。蒲凯[5]在创建地质空间数据库系统时，把Oracle当做存储管理所有空间数据的基础工具，在ArcSDE的支持下，实现了空间数据的统一存储和管理，并进行应用示范。刘灿娟[6]在Microsoft Visual Studio.NET环境下，结合Oracle 10G数据库系统相关技术，根据应用要求，设计并完成地质数据的存储管理工作。然而，数据的不断积累，海量数据的出现，使得通用的关系型数据库无法很好地应对数据规模剧增时系统扩展性和性能问题，Hadoop的出现，则很好地解决了大数据存储的问题。由于Hadoop处理过程主要是通过Hadoop分布式文件系统（HDFS）来实现的[7], 因此本文采用Hadoop中的HDFS为底层存储架构，HBase为元数据存储机制的方式存储地质矿产大数据。实验证明，该存储方式明显优于传统关系数据存储方式。
1. Hadoop系统框架
地质数据具有多源、多元、异构、时空性、方向性、相关性、随机性、模糊性、非线性等特性[8]，因此地质大数据具有多样性，存在多维度、结构化与非结构化并用的特点。HDFS[9]作为Hadoop的分布式文件系统，可用于保存基本是顺序访问的海量结构化、非结构化等数据，并且适合运行在普通硬件上[10]，而HBase的主要优势为快速随机访问数据。因此，只有将两者结合起来才能更好地实现影像、矢量、文本等海量地质矿产数据的存储与检索。本文存储方法设计如下：将实际地质矿产数据存储到HDFS中（对于存储大的（基于MB、GB）影像地质数据条目直接存储到HDFS中，小的（KB）文本、矢量地质数据条目，先进行小文件合并等措施，再存储到HDFS中），HBase存储其索引数据。应用程序在一个新的HDFS文件中写入结果，同时更新基于HBase的元数据。
Hadoop为用户提供了一个开源的大数据分布式处理框架，包括分布式文件系统HDFS、分布式处理模型MapReduce、非关系型数据库HBase，该框架帮助用户在了解底层实现的基础上，通过函数编程和操作接口，完成对大数据的存储和处理[11]。其系统框架如图 1所示。
（1）HDFS：分布式文件系统，在整个软件生态系统中提供数据存储的功能，并提供数据的流式访问接口，相比集中式数据存储机制，它可以有效地提高系统数据访问的吞吐量。HDFS部署在一个由大量普通PC服务器组成的集群中，采用一次写入、多次读取的文件访问模型，并具有很强的容错处理能力和良好的平台可移植性。
（2）MapReduce：一个处理大数据的编程模型[12]，采用Map和Reduce两个简化过程将复杂的任务分成相互独立的子问题，并自动调度计算处理节点来处理这些子问题。利用数据/代码互定位技术有效地减少节点间的数据通信，还能够为应用开发者隐藏系统层的功能实现细节。
（3）HBase：一种提供列存储、实时读写，高性能的分布式数据库，是非关系型数据库NoSQL的一种具体实现，用于存放非结构化和半结构化的松散数据，并且可以通过主键（RowKey）和主键range来实现检索功能。
2. 地质矿产大数据存储模型
2.1 基于关系型数据库的存储模型
关系型数据库指采用关系模型来组织数据的数据库，模型由IBM研究员Codd在1970年首次提出[13]。地理信息系统（Geographic Information System, 简称GIS）以计算机科学为主要支撑技术，专门针对空间数据的应用问题[14]，能够用于地质数据的查询管理、空间分析、可视化等操作，将其与Oracle大型空间数据库结合并基于一定的标准进行数据库系统开发，是海量地质数据存储管理的现有存储方案。近半个世纪以来，中国针对国家级基础地质数据库的投资种类达到上百种，投资的数据总量也已超过100TB，并且目前数量还在急剧增长[15]。这些数据管理大多以GIS为基础实现平台，集合大型数据库存储技术、计算机技术进行组织、生产和展示，部分甚至采取三维建模等技术进行一体化管理。然而，这些管理应用中很少能充分利用数据库的功能，因此，面对海量数据的存储，将变得不堪重负。
2.2 基于非关系型数据库的存储模型
随着地质数据的快速增长，基于关系型数据库的存储模式已经很难高效地满足大量地质数据的存储要求。目前，云平台的出现及“大数据”的发展使得越来越多的学者倾向于运用Hadoop相关技术进行数据存储的研究。本文目标是结合HDFS与HBase的优点，设计一个复合式的地质大数据存储系统，将原始地质数据资料存储在HDFS中，将地质矿产元数据存储于HBase中，以实现大数据存储与高效率检索，系统架构如图 2所示。
HDFS主要由NameNode和一系列DataNode组成，其体系结构见图 3，其中，NameNode和Secondary NameNode协同合作管理HDFS的目录树和相关的元数据文件。NameNode的fsimage是元数据镜像文件，存储着整个文件系统的目录树，edits则是元数据的操作日志，且其大小影响NameNode的启动和运行速度，所以需要Secondary NameNode定期通过http get从NameNode中获取上述2个文件，并合并生成新的fsimage文件返回给NameNode，随后NameNode更新本地的fsimage文件，并清空edit操作日志内容，从而提升HDFS节点管理性能[16]。
2.2.1 基于HDFS的地质大数据存储
HDFS的文件操作主要包括读、写2个主要流程，两者均由HDFS客户端发起。
（1）文件写入：①客户端通过调用Distributed File System对象中的create（）方法，初始化写入文件实例。②NameNode首先会验证系统中没有要创建的文件，并确保客户端有创建文件的权限，否则会抛出I/O异常。③若校验通过，FSData Output Stream会将文件切割成多个packets，并协调NameNode和DataNode，通过write packet方法将对象写入DataNode。④最后一个文件写入成功后会返回一个ack packet给客户端，客户端则调用close方法关闭所有数据流，并将完成的信息汇报给NameNode。
本文选取重庆市矿产潜力评价数据中的“重庆市城口南大巴山式侵入岩型金矿预测工作区遥感影像图GEOTIFF”（230MB）为测试数据将其存储到HDFS中，核心代码如下：
FSDataOutputStream fsDataOutputStream = fs. create（new Path（“/opt/hadoop/dfs/Data/重庆市/金矿/重庆市城口南大巴山式侵入岩型金矿预测工作区遥感影像图GEOTIFF.tif”））；//在hdfs上打开一个文件输出流FileInputStream fileInputStream = new FileInputStream（“/home/hadoop/Desktop/重庆市城口南大巴山式侵入岩型金矿预测工作区遥感影像图GEOTIFF.tif”）；//在hdfs上打开一个文件输入流IOUtils.copy（file Input Stream，fs Data Output Stream）。
（2）文件读取：①客户端调用open函数，向远程的NameNode发起RPC请求，调用元数据节点，得到文件的数据块信息。②客户端在打开的文件流上调用read（）函数。DFSInputStream连接保存距离此文件最近的DataNode，并把要读的数据返回客户端。③读取完当前block的数据后，就会断开到DataNode的链接，然后选择下一个DataNode来获取下一个数据块。④所有数据读取完毕后，调用close函数来关闭数据流。
本文选取hdfs中data下的矢量数据为读取对象，核心代码如下：
RemoteIterator <LocatedFileStatus> listFiles = fs.listFiles（new Path（“/opt/hadoop/dfs/data”）, true）;
While（ListFiles.hasNext（））{
LocatedFileStatus file = ListFiles.next（）;
System.out.println（file.getPath（）.getName（））; }
2.2.2 基于HBase的地质元数据存储
HBase是一个分布式的、面向列的开源数据库，是基于Google论文《Bigtable:A distributed storage syatems for structured data》 [17]开源实现，HBase将HDFS作为其底层存储系统；以HadoopMapReduce进行海量数据分析；采用Zookeeper进行数据的统一管理。
以矢量地质数据为例，本文采用每个省份数据作为一个表，每个矿种对应一行，用行健表示各省份对应的矿种信息；用不同的列族代表“综合研究成果资料”、“省级矿产资料”；用不同的列代表不同的要素图（表 1）。
表 1 HBase列式存储设计Table 1. HBase column storage design table行健 时间戳 列族SPE 列族S_SPE 矿种ID 专题ID 图件ID 路径 其他属性 矿种ID 专题ID 图件ID 路径 其他属性 矿种RowKey TimeStamp MID SID FID LCT … MID SID FID LCT … 在表中，列族“SPE”代表“综合研究成果资料”，列族“S_SPE ”代表“省级矿产资料”。列族“SPE”下的‘MID’代表矿种的编码信息（编码信息由国家统一规定），如01代表“铁”，02代表“锰”等，后面命名规则参照前述方式。该事例表示的为“湖南省硫矿种”的编码信息和存储路径信息：
4319 column=SPE: MID, timestamp= 1478462443280, value=4319
4319 column=SPE: LCT, timestamp= 1478462443280, value=hdfs://192.168.20.51/opt/
hadoop/dfs/data/HUNAN
2.3 小文件存储问题
HDFS是一个具有高容错性、成本低廉性等特点的分布式文件系统，被设计用来对大文件进行流式存储，而在处理小文件时则会产生一些问题[18]。地质空间数据与其他空间数据相比，具有数据量大、专题内容多、信息量大等特点，且以矢量图层、遥感影像等大文件数据为主，非常适合于HDFS存储方法。以重庆市的金矿种为例，矿产潜力评价数据如图 4。
![]()
图 4 潜力评价数据（以重庆市金矿种为例）Figure 4. Potential evaluation data exemplified by a gold mine in Chongqing经实验发现，矢量图层存储时并非以各地区的“图件”为存储单元存储到一个至几个Block内存中，而是以WL、WT等文件单独进行存储，分别占据一个Block内存，形成无数小文件，造成HDFS存储瓶颈（图 5）。
图 4显示“重庆市城口石门口南大巴山岩浆热液型金矿区典型成矿要素图”中的一个矢量要素。该图只有15892个字节，却占据了一个Block节点，分配单独的split，给集群带来了存储压力。
2.4 Hadoop自带小文件优化方案
对于小文件问题，Hadoop系统自身提供了3个解决方案，分别是Hadoop Archive[19]，Sequence File[20]，CombineFileInputFormat。
HAR（Hadoop Archive）文件归档技术主要通过将HDFS上的小文件打包成HAR，缓解大量小文件消耗NameNode内存的问题。但是读取HAR中的文件需要读取两层index文件及文件本身数据，且HAR不支持存档文件的压缩，因此使用HAR处理小文件效率较低。
SequenceFile是HDFS提供的一种二进制文件技术，通过 <key, value>对序列化到SequenceFile文件实现多个小文件合并，同时支持数据块的压缩，显著减少了NameNode的内存及数据节点的磁盘空间。然而该方法没有建立相应的小文件到大文件的映射，读取效率低。
CombineFileInputFormat是一种新的输入格式，将HDFS上的多个文件合并成一个单独的split进而提高效率。该方法适用于HDFS中已经存在大量小文件的情况，并且CombineFileInputFormat是一个抽象类，需要用户自己实现实体类。
参考以上方法，CombineFileInputFormat方案更可行，然而，大量小文件数据尚未存到Hadoop平台，若先进行存储，再优化，效率则更低，所以最理想的策略是将“图层”下的各矢量因素合并为一个大文件，以此来提高小文件的存储性能。
2.5 IPutMerge解决方案
针对HDFS小文件存储，《hadoop实战》 [21]提供了一个PutMerge程序，用于将本地磁盘文件合并，然后再复制到HDFS。然而该方案仅面向本地单个文件夹下的数据，对于大量不同路径的数据存储，效率则较低，因此，本文设计了一种IputMerge（Improved PutMerge）程序，该方法可以实现数据的跨路径存储，操作更便捷。
管理split的总时间和构建map任务的总时间决定作业的执行时间，因此合理的分片大小决定负载平衡的质量，本文将分片大小设置为Block的大小128MB。其次，调用Configuration conf得到所需要的FileSystem实例。Configuration对象封装了客户端或服务器的配置，通过设置配置文件读取类路径实现core-site.xml。有了FileSystem实例后，调用factory方法FileSystem.getLocal（Configuration conf）。Hadoop文件API使用Path对象编制文件和目录名，使用FileStatus对象存储文件和目录的元数据。IPutMerge程序将调用递归算法合并本地不同目录中的所有文件。使用FileSystem的ListStatus（）方法得到本地的目录。FSDataInputStream是Java标准类java.io.DataInputStream的子类，可以提供随机访问的功能，本文用该方法将数据写入HDFS。为实现IPutMerge，创建一个for循环逐一读取inputFiles中的所有文件。部分代码如下所示：
Configuration conf = new Configuration（）;
......
Path localPath = new Path（（“/home/hadoop/Desktop/重庆市/金矿典型矿床成矿要素图”））;
Path hdfsPath = new Path（（“/opt/hadoop/dfs/data/重庆市/金矿种潜力评价图库/矿产及其预测/金矿典型矿床成矿要素图”））;
if（out==null）{
out = hdfsFS.create（hdfsPath）;
}
readLocalFileWriteToHDFS（localPath, localFS）;
......
FileStatus[] inputFiles = localFS.listStatus（localPath）;
for（int i =0; i < inputFiles.length; i++）{
boolean isDir = inputFiles[i].isDirectory（）;
if（isDir）{
readLocalFileWriteToHDFS（inputFiles[i].getPath（）, localFS）;
}
if（inputFiles[i].isFile（））{
……
}
}
}
最后，调用MapReduce程序，将数据上传到HDFS之上。Shell命令为：
[hadoop@hdmaster1 Desktop] $ hadoop jar HDFSpd.jar数据.HDFSpd2
3. 实验验证
3.1 小文件优化对比实验
为了验证优化后的方案更具优势，本文从2个方面进行对比测试：①实验数据在HDFS上的写入性能；②文件操作过程中HDFS的NameNode内存占用结果。
（1）实验环境
本文实验采用4个节点的Hadoop集群，其中包括1个NameNode和3个DataNode，具体配置如表 2。
表 2 集群配置情况Table 2. Cluster configuration组件 配置 Hadoop版本 Hadoop2.7.2 操作系统 Red Hat Enterprise Linux Server release 6.7(Santiago) Linux内核版本 2.6.32-573.e16.x86_64 JDK 1.7.0_91 网络带宽 100MB NameNode 八核2.4.0GHzCPU, 16G内存，600G硬盘，数量1 DataNode 八核2.4.0GHzCPU, 8G内存，600G硬盘，数量3 （2）写操作性能对比测试
本实验分别选取500, 1000, 10000, 20000, 30000, 40000个大小为1~1000KB的小文件为存储数据，分别进行数据存储，记录它们各自的存储时间，取3次实验结果的平均值，实验结果见图 6。
图 6表明，采用IPutMerge方法优化后的HDFS的写操作明显快于原始HDFS。这是因为，每写入一个文件到HDFS，NameNode都要对其进行创建和分配数据块的操作，而本文优化后的HDFS中，当n个小文件组成的大文件写入时，NameNode仅被调用一次。另外，Client对大文件使用的缓存机制，每128MB数据的请求被发送到NameNode，而不是发送每个小文件。
（3）内存占用对比测试
海量小文件会耗费主节点内存，造成NameNode瓶颈问题，进而严重影响HDFS的性能，为了分析基于HDFS的分布式文件存储系统中的NameNode占用情况并设计实验。
本实验分别选取500, 1000, 10000, 20000, 30000, 40000个大小为1~1000KB的小文件作为实验数据，在半小时内不断写入各文件到HDFS，并观察NameNode所占系统内存变化，取3次实验结果的平均值。内存使用量结果如图 7所示。
图 7显示了2种情况下的内存使用情况，从结果可以看出，在小文件数量小于1000个的情况下，优化效果并不明显，然而随着小文件数量的急剧增长，优化后的NameNode所占内存量远远低于优化前的HDFS中。这是因为，NameNode主要保存文件元数据及文件所在块的元数据。优化前的HDFS进行数据存储时，每个小文件占据了一个Block，则NameNode需占用大量内存来保存这些数据。而优化后的HDFS，NameNode保存了单个组合文件及Block的元数据，因此，系统内存占用量则会大幅度降低。
3.2 Hadoop与Oracle对比实验
本文实验通过对比2种空间数据存储方式：①传统空间数据存储方式，即采用Oracle与MapGis相结合进行空间数据存储。② Hadoop平台，通过HDFS进行数据存储。分析实验结果，验证本文方法更具优势。
（1）实验环境
表 3给出了测试用到的硬件环境，其中Oracle采用单节点，Hadoop则是在4台电脑上部署的集群。
表 3 测试数据库硬件环境Table 3. Test database hardware environment组件 Oracle Hadoop 操作系统 Windows7专业版 Red Hat Enterprise Linux Server release 6.7(Santiago) 服务器数量 1台(单节点) 4台(Hadoop集群） 网络带宽 100MB 100MB CPU性能 四核 八核 CPU主频 2.80GHz 2.40GHz 内存 8GB 8GB （2）数据导入对比实验
本实验以重庆市矿产资源潜力评价数据作为研究对象，选取100, 1000，10000, 50000, 100000个大小为1~1000KB的小文件为存储数据，分别采用Oracle与Hadoop先后进行数据存储，求取3次对比实验结果的平均值（图 8）。
从图 8可以看出，当文件数量较少时，2种存储方式效率持平，改进后的Hadoop存储方式则略占优势；当数据量很少时，较之传统的Hadoop存储，Oracle效率可能更高，这是因为使用Hadoop进行数据存储时，涉及到NameNode的元数据存储及DataNode的选取，耗费时间较长，优化后的Hadoop则仅需要调用少量的NameNode，大大节省了系统开销时间；当数据量较大时，尤其达到10000个以后，Hadoop在地质数据存储方面则显示出了巨大优势，写入数据所需时间大为缩短。而传统的Oracle数据库，数据存储由于主要集中在一台物理存储节点上，随着数据量不断增加大致呈线性增长。
3.3 数据导出对比实验
本实验以重庆市矿产资源潜力评价数据作为研究对象，选取100, 1000，10000, 50000, 100000个大小为1~1000KB的小文件为研究数据，将Oracle与Hadoop中存储的数据分别进行下载，记录它们各自的响应时间，对比3次实验结果，求取其平均值。
对比2种方式，从实验结果（图 9）可以看出，在数据量很小时，使用Hadoop进行数据读取并没有太大优势，比起传统Oracle仅略占优势。这是因为，当数据量少时，Hadoop本身的系统开销需要花费更多的数据处理时间，然而改进后的Hadoop将多个小文件合并为一个大文件，大大减少了访问数据量。因此，比起存储小数据量更占优势的Oracle来说，效率还要稍高。当数据量大幅度上升时，尤其数据量超过10000个时，Oracle的单节点数据输出则较吃力，系统响应时间大致呈线性增长，而Hadoop的多节点并行工作突显出了更大的优势。
4. 结语
大数据时代的到来带动了地质大数据的发展，而地质数据对城市规划、建设、安全等都具有重要的指导意义，如何有效存储地质数据成为了亟待解决的现实问题。本文针对海量地质矿产数据的存储问题，提出了一种全新的数据存储方案，对地质行业的发展具有一定的参考价值。实验表明：
（1）HDFS不适合海量地质小文件的存储，而优化后的HDFS数据存储则具有显著优势；
（2）当数据量很少时，基于Hadoop的存储效率优势并不明显，而当数据量急剧增长时，Hadoop则展示出了海量数据的存储优势；
（3）本文提出的方案是有效的，为地质数据的存储与管理提供了一种可行的技术方案。
致谢: 审稿专家提出宝贵的修改意见，野外工作中得到甘肃省地质调查院张渊工程师、马国杰助理工程师、张兵兵高级工程师的协助，中科院地球化学研究所张辉研究员对本文提出了建设性的意见，在此一并致谢。 -
![]()
![]()
图 2 多坝沟岩体地质简图和岩石学特征
1—第四系；2—敦煌岩群第四岩组；3—敦煌岩群第三岩组；4—敦煌岩群第二岩组；5—二长花岗岩；6—花岗岩脉；7—面理产状；8—片理产状；9—片麻理产状；10—采样点及编号。Kfs—钾长石；Bi—黑云母
Figure 2. Simplified geological map and petrographic features of Duobagou granites
![]()
图 3 多坝沟花岗岩锆石阴极发光图像
Figure 3. Cathodoluminescence images of zircons from Duobagou granite
![]()
图 4 多坝沟花岗岩U-Pb谐和图和206Pb/238U年龄加权平均值图
Figure 4. U-Pb concordia diagram and weighted average 206Pb/238U age diagram of zircons from the Duobagou granites
![]()
![]()
图 6 多坝沟花岗岩稀土元素配分曲线(a)和微量元素蛛网图(b)(标准化值据参考文献[30])
Figure 6. Chondrite-normalized REE patterns(a) and primitive mantle-normalized trace element spidergrams(b) of the Duobagou granite
![]()
图 7 多坝沟花岗岩10000Ga/Al-TFeO/MgO(a)[34]和Zr-TiO2(b)图解
Figure 7. 10000Ga/Al versus TFeO/MgO (a)and Zr versus TiO2 diagrams (b)for Duobagou granite
![]()
图 8 多坝沟花岗岩Y-Sr/Y(a)和YbN-(La/Yb)N(b)图解[35]
Figure 8. Sr/Y versus Y (a)and YbN versus (La/Yb)N (b) diagrams for Duobagou granite
![]()
图 9 多坝沟花岗岩源岩Ta/Yb-Th/Yb(a)和源区SiO2-Mg#(b)判别图
Figure 9. Ta/Yb versus Th/Yb (a)and SiO2 versus Mg#(b)diagrams for Duobagou granite
![]()
图 10 多坝沟花岗岩构造环境判别图解[31]
ORG—洋脊花岗岩；WPG—板内花岗岩；VAG—火山弧花岗岩；syn-COLG—同碰撞花岗岩
Figure 10. Discrimination diagrams of tectonic setting for Duobagou granites
表 1 多坝沟花岗岩LA-ICP-MS锆石U-Th-Pb年龄测定结果
Table 1 Zircon LA-ICP-MS U-Th-Pb ages of the granite from Duobagou region
测点号 含量/10-6 Th/U 同位素比值 表面年龄/Ma 232Th 238U 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 206Pb/207Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ dtw05-01 437 583 0.75 0.05156 0.0009 0.2823 0.0051 0.03972 0.0005 266 20 253 4 251 3 dtw05-02 672 810 0.83 0.05058 0.0013 0.2744 0.0068 0.03935 0.0005 222 33 246 5 249 3 dtw05-03 541 569 0.95 0.05684 0.0021 0.3591 0.0127 0.04582 0.0007 485 49 312 9 289 5 dtw05-04 447 438 1.02 0.05119 0.0009 0.2875 0.0054 0.04075 0.0005 249 21 257 4 257 3 dtw05-05 724 670 1.08 0.05504 0.0011 0.3075 0.0064 0.04053 0.0005 414 24 272 5 256 3 dtw05-06 389 505 0.77 0.05412 0.0011 0.2992 0.0061 0.04011 0.0006 376 23 266 5 254 3 dtw05-07 478 693 0.69 0.05119 0.0009 0.2734 0.0051 0.03875 0.0005 249 21 245 4 245 3 dtw05-08 651 700 0.93 0.05163 0.0011 0.2911 0.0064 0.04091 0.0006 269 27 259 5 258 3 dtw05-09 538 445 1.21 0.05057 0.0009 0.3003 0.0056 0.04311 0.0006 221 21 267 4 272 4 dtw05-10 441 573 0.77 0.05171 0.0011 0.2790 0.0060 0.03916 0.0005 273 26 250 5 248 3 dtw05-11 537 639 0.84 0.05204 0.0009 0.2853 0.0052 0.03979 0.0005 287 20 255 4 252 3 dtw05-12 620 608 1.02 0.05197 0.0011 0.2952 0.0062 0.04122 0.0006 284 25 263 5 260 3 dtw05-13 549 603 0.91 0.05201 0.0009 0.2781 0.0049 0.03879 0.0005 286 19 249 4 245 3 dtw05-14 296 352 0.84 0.05129 0.0008 0.2823 0.0048 0.03996 0.0005 254 18 253 4 253 3 dtw05-15 357 305 1.17 0.05080 0.0009 0.2824 0.0053 0.04034 0.0005 232 21 253 4 255 3 dtw05-16 484 701 0.69 0.05382 0.0010 0.2909 0.0055 0.03921 0.0005 364 21 259 4 248 3 dtw05-17 384 463 0.83 0.05194 0.0011 0.2916 0.0064 0.04074 0.0006 283 27 260 5 257 3 dtw05-18 466 542 0.86 0.05241 0.0009 0.2818 0.0051 0.03901 0.0005 303 20 252 4 247 3 dtw05-19 375 387 0.97 0.05192 0.0009 0.2879 0.0051 0.04023 0.0005 282 20 257 4 254 3 dtw05-20 422 570 0.74 0.05581 0.0011 0.3019 0.0063 0.03925 0.0005 444 24 268 5 248 3 
下载: 导出CSV
表 2 多坝沟花岗岩体主量、微量和稀土元素组成
Table 2 Geochemical composition of major, trace and rare earth elements of Duobagou granite
样品号 PM04-1 DYQ23 DYQ24 DYQ25 DYQ26 SiO2 69.75 68.32 68.69 70.59 72.82 Al2O3 14.66 14.78 14.74 14.49 13.73 TiO2 0.35 0.39 0.40 0.32 0.25 Fe2O3 0.92 1.14 1.01 0.94 0.85 FeO 1.77 1.89 2.01 1.56 1.17 CaO 2.55 2.13 2.55 2.25 1.57 MgO 0.91 1.06 1.10 0.85 0.59 K2O 4.30 4.14 4.27 4.55 4.58 Na2O 3.48 3.78 3.48 3.39 3.43 MnO 0.07 0.08 0.08 0.07 0.05 P2O5 0.11 0.14 0.14 0.12 0.08 H2O+ 0.71 1.24 1.03 0.52 0.66 H2O- 0.10 0.14 0.14 0.14 0.11 烧失量 0.99 2.01 1.34 0.72 0.76 总计 99.85 99.86 99.84 99.85 99.88 A/NK 1.41 1.38 1.42 1.38 1.29 A/CNK 0.97 1.01 0.98 0.99 1.02 Mg# 39 39 40 39 35 Li 27.5 14.4 27.7 19.3 17.9 Be 2.42 2.57 2.40 1.88 2.74 Sc 4.61 5.40 4.75 4.51 4.05 V 41.3 37.0 36.9 32.2 28.0 Cr 5.41 7.59 6.36 6.41 5.43 Co 18.2 5.54 5.83 4.65 3.65 Ni 3.44 4.15 3.17 2.83 2.30 Cu 12.1 5.36 6.10 16.0 24.6 Zn 54.8 51.8 62.6 40.7 48.0 Ga 15.7 16.2 15.6 14.0 15.6 Rb 167 125 99.2 147 163 Sr 348 289 334 308 230 Zr 184 199 180 144 145 Nb 10.1 7.83 8.75 6.58 7.98 In 0.03 0.03 0.04 0.02 0.02 Cs 6.66 2.66 5.83 5.34 6.27 Ba 812 536 622 602 581 Hf 6.26 6.71 5.71 4.31 4.51 Ta 0.70 0.96 0.96 0.69 1.08 W 190 0.83 0.37 0.69 2.87 Pb 19.0 23.7 19.9 20.8 24.8 Bi 0.06 0.16 0.17 0.03 0.09 Th 9.25 15.1 8.73 11.1 17.7 U 1.56 2.49 1.32 1.70 3.14 F 342 467 677 451 467 La 21.3 24.4 19.8 23.2 21.5 Ce 44.1 49.4 43.9 48.6 44.1 Pr 5.29 5.74 5.31 5.51 5.01 Nd 17.5 20.5 19.4 19.1 16.9 Sm 3.00 3.71 3.61 3.13 2.79 Eu 0.81 0.91 0.96 0.85 0.59 Gd 2.72 3.51 3.33 3.13 2.82 Tb 0.42 0.56 0.56 0.43 0.39 Dy 2.25 3.24 3.13 2.33 2.10 Ho 0.43 0.67 0.64 0.47 0.45 Er 1.29 1.84 1.77 1.33 1.28 Tm 0.21 0.34 0.33 0.24 0.26 Yb 1.60 2.10 1.98 1.47 1.55 Lu 0.26 0.33 0.32 0.25 0.29 Y 12.4 18.2 17.6 12.9 12.8 ΣREE 113.5 135.4 122.6 123.0 112.8 LREE 92.0 104.6 93.0 100.4 90.9 HREE 9.2 12.6 12.1 9.7 9.1 LREE/HREE 10.0 8.3 7.7 10.4 10.0 (La/Yb)N 9.58 8.32 7.17 11.27 9.96 δEu 0.87 0.77 0.84 0.83 0.65 注:主量元素含量单位为%，微量和稀土元素为10-6，其中Au元素单位为10-9
下载: 导出CSV
表 3 多坝沟花岗岩全岩Sr-Nd同位素数据
Table 3 Sr-Nd isotopic data of the Duobagou granite
样号 Rb Sr 87Rb/86Sr 87Sr/86Sr 2σ (87Sr/86Sr)i Nd Sm 147Sm/144Nd 143Nd/144Nd 2σ (143Nd/144Nd)i εNd(0) εNd(t) fSm/Nd TDM2 /10-6 /10-6 /Ma PM04-1 166.56 347.65 1.38734 0.712058 0.000007 0.707083 17.49 3 0.103323 0.51237 0.000004 0.5122 -5.23 -2.22 -0.47 1102
下载: 导出CSV
-
张建新, 李怀坤, 孟繁聪, 等.塔里木盆地东南缘(阿尔金山)"变质基底"记录的多期构造热事件:锆石U-Pb年代学的制约[J].岩石学报, 2011, 27(1):23-46. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201101002 Zhao Y, Sun Y, Diwu C R, et al. The Dunhuang block is a Paleozoic orogenic belt and part of the Central Asian Orogenic Belt (CAOB), NW China[J]. Gondwana Research, 2016, 30:207-223. doi: 10.1016/j.gr.2015.08.012
张志诚, 郭召杰, 邹冠群, 等.甘肃敦煌党河水库TTG地球化学特征, 锆石SHRIMP U-Pb定年及其构造意义[J].岩石学报, 2009, 25(3):495-505. http://d.old.wanfangdata.com.cn/Periodical/ysxb98200903002 朱涛, 王洪亮, 徐学义, 等.敦煌地块南缘石炭纪埃达克岩的发现及其地质意义[J].岩石学报, 2014, 30(2):491-502. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201402013 刘东晓, 王玉玺, 王金荣, 等.塔里木板块东南缘早二叠世角闪石岩的发现及大地构造意义[J].兰州大学学报(自然科学版), 2016, 52(2):167-178. http://www.cnki.com.cn/Article/CJFDTOTAL-LDZK201602004.htm 王玉玺, 第鹏飞, 陈万峰, 等.塔里木板块东南缘多坝沟金矿区花岗岩类年代学、地球化学特征及找矿意义[J].西北地质, 2017, 50(1):134-150. doi: 10.3969/j.issn.1009-6248.2017.01.013 王玉玺, 王金荣, 周小玲, 等. Columbia超大陆裂解:来自塔里木克拉通东南缘大红山A型花岗岩的证据[J].地质学报, 2017, 91(11):2369-2386. doi: 10.3969/j.issn.0001-5717.2017.11.001 许志琴, 杨经绥, 张建新, 等.阿尔金断裂两侧构造单元的对比及岩石圈剪切机制[J].地质学报, 1999, 73(3):193-205. doi: 10.3321/j.issn:0001-5717.1999.03.001 Lu S N, Zhao G C, Wang H C, et al. Precambrian metamorphic basement and sedimentary cover of the North China Craton:a review[J]. Precambrian Research, 2008, 160(1):77-93. doi: 10.1016-j.precamres.2007.04.017/
辛后田, 刘永顺, 罗照华, 等.塔里木盆地东南缘阿克塔什塔格地区新太古代陆壳增生:米兰岩群和TTG片麻岩的地球化学及年代学约束[J].地学前缘, 2013, (1):240-259. http://d.old.wanfangdata.com.cn/Periodical/dxqy201301020 赵燕, 朱涛, 王洪亮, 等.敦煌三危山地区晚泥盆世斜长花岗岩的发现及其地质意义[J].岩石学报, 2015, (7):1855-1869. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201507004 郑荣国, 吴泰然, 张文, 等.阿拉善地块北缘雅干花岗岩体地球化学, 地质年代学及其对区域构造演化制约[J].岩石学报, 2013, 29(8):2665-2675. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201308004 Zheng R G, Wu T R, Zhang W, et al. Late Paleozoic subduction system in the northern margin of the Alxa block, Altaids:Geochronological and geochemical evidence from ophiolites[J]. Gondwana Research, 2014, 25(2), 842-858 doi: 10.1016/j.gr.2013.05.011
Ge R F, Zhu W B, Wu H L, et al. The Paleozoic northern margin of the Tarim Craton:Passive or active?[J]. Lithos, 2012, 142/143:1-15. doi: 10.1016/j.lithos.2012.02.010
Zong K Q, Zhang Z M, He Z Y, et al. Early Palaeozoic highpressure granulites from the Dunhuang block, northeastern Tarim Craton:constraints on continental collision in the southern Central Asian Orogenic Belt[J]. Journal of Metamorphic Geology, 2012, 30(8):753-768. doi: 10.1111/jmg.2012.30.issue-8
He Z Y, Zhang Z M, Zong K Q, et al. Metamorphic p-T-t evolution of mafic HP granulites in the northeastern segment of the Tarim Craton (Dunhuang block):Evidence for early Paleozoic continental subduction[J]. Lithos, 2014, 196(5):1-13. http://www.sciencedirect.com/science/article/pii/S0024493714000760
王楠, 吴才来, 马昌前, 等.敦煌地块三危山地区花岗岩体地球化学, 锆石U-Pb定年及Hf同位素特征[J].地质学报, 2016, 90(10):2681-2705. doi: 10.3969/j.issn.0001-5717.2016.10.010 Zhang J X, Yu S Y, Gong J H, et al. The latest NeoarcheanPaleoproterozoic evolution of the Dunhuang block, eastern Tarim craton, northwestern China:Evidence from zircon U-Pb dating and Hf isotopic analyses[J]. Precambrian Research, 2013, 226:21-42. doi: 10.1016/j.precamres.2012.11.014
辛后田, 赵凤清, 罗照华, 等.塔里木盆地东南缘阿克塔什塔格地区古元古代精细年代格架的建立及其地质意义[J].地质学报, 2011, 85(12):1977-1993. http://d.old.wanfangdata.com.cn/Periodical/dizhixb201112002 陈柏林, 杨屹, 王小凤, 等.阿尔金北缘大平沟金矿床成因[J].矿床地质, 2005, 24(2):168-178. doi: 10.3969/j.issn.0258-7106.2005.02.009 李月臣, 陈柏林, 陈正乐, 等.阿尔金北缘红柳沟-拉配泉一带铜金矿床硫同位素特征及其意义[J].地质力学学报, 2007, 13(2):131-140. doi: 10.3969/j.issn.1006-6616.2007.02.006 孟繁聪, 张建新, 相振群, 等.塔里木盆地东北缘敦煌群的形成和演化:锆石U-Pb年代学和Lu-Hf同位素证据[J].岩石学报, 2011, 27(1):59-76. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201101004 朱江, 吕新彪, 彭三国, 等.甘肃北山小西弓金矿区石英正长斑岩LA-ICP-MS锆石U-Pb年龄和地球化学特征[J].地质通报, 2015, 34(8):1460-1469. doi: 10.3969/j.issn.1671-2552.2015.08.006 Jackson S E, Pearson N J, Griffin W L, et al. The application of laser ablation-inductively coupled plasma-mass spectrometry to in situ U-Pb zircon geochronology[J]. Chemical Geology, 2004, 211(1):47-69. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=17d96cdc3a96b2aa4c52b4ebb0b5da26
Van Achterbergh E, Ryan C G, Jackson S E, et al. Data reduction software for LA-ICP-MS in the Earth Sciences[J]. Mineralogical Association of Canada, 2001, 29:239-243. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=3861e59a4546c7a0b35d78f6e5f71b33
Anderson T. Correction of common lead in U-Pb analyses that do not report 204Pb[J]. Chemical Geology, 2002, 192(1):59-79. http://www.sciencedirect.com/science/article/pii/S000925410200195X
Ludwig K R. Mathematical-statistical treatment of data and errors for 230Th/U geochronology[J]. Reviews in Mineralogy and Geochemistry, 2003, 52(1):631-656. http://www.sciencedirect.com/science/article/pii/S0012821X0201107X
Maniar P, Piccoli P. Tectonic discrimination of granitoids[J]. Geol. Soc. Amer. Bull., 1989, 101:635-643. doi: 10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
Collins W J, Beams S D, White A J K, et al. Nature and origin of a type granites with articular reference to southeastern Australia[J]. Contrbutions to Mineralogy Petrology, 1982, 80:189-200. doi: 10.1007/BF00374895
Sun S S, McDonough W F. Chemical and isotopic systematic of oceanic basalts: Implications for mantle composition and processes[C]//Saunders A D, Norry M J. Magmatism in the Ocean Basins.Geol. Soc. Special Publication, 1989, 42: 313-345.
Pearce J A, Harris N B W, Tindle A G. Trance element discrimination diagrams for the tectonic interpretation of granitic rocks[J]. Journal of Petrology, 1984, 25:956-983. doi: 10.1093/petrology/25.4.956
吴福元, 李献华, 杨进辉, 等.花岗岩成因研究的若干问题[J].岩石学报, 2007, 23(6):1217-1254. doi: 10.3969/j.issn.1000-0569.2007.06.001 张旗, 李承东.花岗岩:地球动力学意义[M].青岛:海洋出版社, 2012:1-81. Whalen J B, Currie K L, Chappell B W. A-type granites:Geochemical characteristics, discrimination and Petrogenesis[J]. Contribution to Mineralogy and Petrology, 1987, 95(4):407-419. doi: 10.1007/BF00402202
Defant M J, Drummond M S. Derivation of some modern arc magmas by melting of young subducted lithosphere[J]. Nature, 1990, 347:662-665. doi: 10.1038/347662a0
张旗, 王焰, 熊小林, 等.埃达克岩和花岗岩:挑战与机遇[M].北京:中国大地出版社, 2008. Sajona F G, Maury R C, Pubellier M, et al. Magmatic source enrichment by slab-derived melts in a young post-collision setting, central Mindanao (Philippines)[J]. Lithos, 2000, 54(3):173-206. doi: 10.1016-S0024-4937(00)00019-0/
Reagan M K, Gill J B. Coexisting calcalkaline and high-niobium basalts from Turrialba Volcano, Costa Rica:Implications for residual titanates in arc magma sources[J]. Journal of Geophysical Research:Solid Earth, 1989, 94(B4):4619-4633. doi: 10.1029/JB094iB04p04619
Gao S, Rudnick R L, Yuan H L, et al. Recycling lower continental crust in the North China craton[J]. Nature, 2004, 432(7019):892-897. doi: 10.1038/nature03162
Rapp R P, Shimizu N, Norman M D, et al. Reaction between slab-derived melts and peridotite in the mantle wedge:experimental constraints at 3.8GPa[J]. Chemical Geology, 1999, 160(4):335-356. doi: 10.1016/S0009-2541(99)00106-0
徐学义, 何世平, 王洪亮, 等.早古生代北秦岭-北祁连结合部构造格局的地层及构造岩浆事件约束[J].西北地质, 2008, 41(1):1-21. doi: 10.3969/j.issn.1009-6248.2008.01.001 Batchelor R A, Bowden P. Petrogenetic interpretation of granitic rock series using multicationic parameters[J]. Chemical Geology, 1985, 48:43-55. doi: 10.1016/0009-2541(85)90034-8
Jian P, Liu D Y, Kröner A, et al. Evolution of a Permian intraoceanic arc-trench system in the Solonker suture zone, Central Asian Orogenic Belt, China and Mongolia[J]. Lithos, 2010, 118(1):169-190. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=40ff8aa320c272880289a116d19525f3
冯继成, 张文, 吴泰然, 等.甘肃北山桥湾北花岗岩体的年代学、地球化学及其地质意义[J].北京大学学报(自然科学版), 2012, 481(1):61-70. http://d.wanfangdata.com.cn/Periodical_bjdxxb201201010.aspx 夏林圻, 李向民, 夏祖春, 等.天山石炭-二叠纪大火成岩省裂谷火山作用与地幔柱[J].西北地质, 2006, 39(1):1-49. doi: 10.3969/j.issn.1009-6248.2006.01.001 赵泽辉, 郭召杰, 韩宝福, 等.新疆三塘湖盆地古生代晚期火山岩地球化学特征及其构造-岩浆演化意义[J].岩石学报, 2006, 22(1):199-214. http://d.old.wanfangdata.com.cn/Periodical/ysxb98200601020 Yang J H, Wu F Y, Wilde S A, et al. Tracing magma mixing in granite genesis:in situ U-Pb dating and Hf-isotope analysis of zircons[J]. Contrib. Mineral. Petrol., 2007, 153:177-190. http://www.tandfonline.com/servlet/linkout?suffix=CIT0084&dbid=16&doi=10.1080%2F00206814.2017.1375438&key=10.1007%2Fs00410-006-0139-7
计量
- 文章访问数: 3905
- HTML全文浏览量: 480
- PDF下载量: 2653
