Geochemical characteristics and chronological significance of Middle Ordovician granite in Halahu area of southern Qilian Mountain
-
摘要:
哈拉湖地区奥陶纪花岗岩相对富钠(Na2O>K2O),具有中钾-高钾钙碱性系列特征,岩石类型以辉长闪长岩和花岗闪长岩为主。岩石稀土及微量元素特征均显示物质来源有地幔物质参与,构造环境判别图显示形成于大洋俯冲的环境。利用LA-ICP-MS进行锆石U-Pb同位素年龄测定,样品锆石206Pb/238U年龄加权平均值为466±4 Ma,限定该花岗岩的形成时代为中奥陶世,该时期北祁连洋发生双向扩张俯冲,南部中祁连陆块遭受俯冲挤压,导致地幔物质上涌,同时陆壳部分熔融,岩体为两者混合形成的产物,该岩体年龄的确定为祁连洋的多阶段俯冲提供了新的证据。
Abstract:The Ordovician granite in the Hala Lake area is characterized by relatively rich sodium (Na2O>K2O) and medium potassium-high potassium calc-alkaline series.Rock types include gabbroic diorite and granodiorite.The whole-rock rare earth and trace element characteristics indicate that the material source was derived from mantle material, and the tectonic environment was most likely the island arc granite formed during the ocean subduction.The zircon U-Pb isotope ages were determined by LA-ICP-MS.These zircons yielded a 206Pb/238U weighted mean age of 466±4 Ma, suggesting that the granite was formed in Middle Ordovician.In the two-way expansion and subduction stage of the North Qilian Ocean, the southern Zhongqilian block was subjected to subduction and extrusion, which caused the mantle material to migrate upward.At the same time, the continental crust was partially melted and the intrusions were formed by the mixture of the two melts.The age of the intrusion provides new evidence for multi-stage subduction of the Qilian Ocean.
-
Keywords:
- South Qilian /
- Ordovician /
- island arc /
- granite
-
地质数据是地质工作中形成的重要信息资源,多年累积形成的海量地质数据信息为国民经济和社会发展提供了重要支撑[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。
经实验发现,矢量图层存储时并非以各地区的“图件”为存储单元存储到一个至几个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)本文提出的方案是有效的,为地质数据的存储与管理提供了一种可行的技术方案。
致谢: 野外工作得到哈拉湖西项目组全体成员的大力帮助,写作过程中北京大学宋述光老师给予了指导,审稿专家提出了宝贵的修改意见,在此表示诚挚的谢意。 -
图 1 哈拉湖地区大地构造位置图(a, 据参考文献[35]修改)和地质简图(b)
Figure 1. The tectonic position (a) and simplified geological map(b) of Halahu area
图 2 哈拉湖中奥陶世侵入岩TAS图解[36]
Figure 2. TAS diagram of Middle Ordovician
表 1 哈拉湖地区中奥陶世侵入岩地球化学分析结果
Table 1 Geochemical compositions of Middle Ordovician intrusions in Halahu area
样品号 ⅠGs2003 IGsC3097
-1QⅣPM21GS5
-1ⅠPm103Gs2
-1ⅠPm103Gs16
-1ⅠGS0809 岩性 石英闪长岩 英云闪长岩 辉长闪长岩 SiO2 56.89 57.45 57.99 50.44 50.45 51.64 TiO2 0.80 0.80 0.77 0.83 0.90 0.74 Al2O3 14.87 15.53 16.32 12.19 12.27 12.26 Fe2O3 3.47 3.97 3.78 2.23 2.14 1.36 FeO 4.12 3.73 4.75 7.65 7.67 8.80 MnO 0.12 0.12 0.16 0.17 0.18 0.20 MgO 4.71 4.41 3.79 9.07 8.99 9.46 CaO 5.14 4.48 3.30 10.08 10.11 6.85 Na2O 2.64 2.81 5.56 2.10 2.34 2.35 K2O 2.64 2.72 0.90 1.59 1.43 2.04 P2O5 0.26 0.27 0.14 0.17 0.23 0.24 烧失量 3.62 3.01 1.84 2.40 2.27 3.24 总量 99.71 99.70 99.73 99.76 99.74 99.74 Mg# 60.32 59.57 51.86 65.41 65.25 64.33 A/NK 2.07 2.05 1.61 2.36 2.27 2.02 A/CNK 0.90 0.99 1.01 0.52 0.52 0.66 La 23.70 22.80 10.77 16.06 18.06 21.46 Ce 45.54 46.30 23.36 32.07 36.55 42.58 Pr 5.72 5.95 3.10 4.13 5.14 5.27 Nd 22.12 22.91 12.81 16.84 21.57 20.34 Sm 4.38 4.44 3.12 3.62 4.71 4.37 Eu 1.09 1.15 0.92 0.96 1.21 1.09 Gd 4.08 3.97 3.24 3.62 4.42 3.79 Tb 0.70 0.63 0.59 0.59 0.71 0.64 Dy 3.82 3.87 3.77 3.27 3.73 3.69 Ho 0.81 0.77 0.84 0.68 0.76 0.76 Er 2.27 2.15 2.46 1.85 2.09 2.07 Tm 0.35 0.34 0.38 0.30 0.32 0.34 Yb 2.32 2.07 2.50 2.00 2.04 2.04 Lu 0.35 0.31 0.37 0.32 0.30 0.30 Sr 474.99 484.95 300.81 428.41 448.96 277.45 Rb 87.91 87.36 20.10 35.85 28.86 61.12 Ba 671.66 760.40 246.84 379.85 439.84 405.40 Th 8.05 9.21 5.41 4.96 5.68 7.13 Ta 0.70 0.69 0.47 1.28 0.82 0.54 Nb 6.31 7.46 6.16 5.53 6.12 7.09 Zr 124.90 129.00 82.20 67.70 75.00 96.40 Hf 4.31 4.52 2.43 2.04 2.36 2.57 Cr 115.81 91.19 30.35 274.80 279.40 369.60 V 207.94 183.50 272.83 252.36 220.00 U 1.57 1.73 1.14 1.03 1.24 1.71 Y 20.22 19.43 20.56 17.40 19.09 17.77 Sc 19.20 18.18 38.23 41.17 30.16 注:主量元素含量单位为%,微量和稀土元素含量单位为10-6 表 2 哈拉湖地区辉长闪长岩锆石U-Th-Pb同位素测定结果(IU-Pb0803)
Table 2 Zircon U-Th-Pb isotope determination results for gabbroic diorite in Halahu area
点号 Pb
/10-6U
/10-6206Pb/238U 1σ 207Pb/235U 1σ 207Pb/206Pb 1σ 208Pb/232Th 1σ 232Th/238U 1σ 206Pb/238U
年龄/Ma1σ 207Pb/235U
年龄/Ma1σ 1 59 777 0.0758 0.0008 0.602 0.012 0.0575 0.0009 0.0303 0.0010 0.0247 0.0005 471 5 478 10 2 98 714 0.1424 0.0014 1.342 0.019 0.0683 0.0008 0.0472 0.0013 0.137 0.001 858 8 864 12 3 41 502 0.0780 0.0008 0.620 0.013 0.0576 0.0011 0.0215 0.0006 0.572 0.003 484 5 490 10 4 44 574 0.0758 0.0007 0.598 0.009 0.0572 0.0008 0.0261 0.0006 0.318 0.002 471 5 476 7 5 23 294 0.0757 0.0007 0.590 0.014 0.0565 0.0012 0.0230 0.0005 0.3880 0.0004 470 5 471 11 6 38 474 0.0758 0.0007 0.592 0.012 0.0566 0.0011 0.0224 0.0005 0.528 0.001 471 5 472 10 7 41 570 0.0716 0.0007 0.565 0.010 0.0572 0.0009 0.0331 0.0008 0.2261 0.0001 446 5 455 8 8 49 590 0.0797 0.0008 0.640 0.009 0.0582 0.0008 0.0130 0.0004 0.924 0.002 494 4 502 7 9 91 1181 0.0762 0.0008 0.610 0.010 0.0581 0.0008 0.0224 0.0006 0.369 0.003 473 5 484 8 10 34 449 0.0743 0.0007 0.590 0.011 0.0576 0.0010 0.0242 0.0006 0.3392 0.0007 462 5 471 9 11 94 1275 0.0727 0.0007 0.576 0.008 0.0575 0.0007 0.0230 0.0005 0.3392 0.0004 453 4 462 7 12 53 676 0.0759 0.0007 0.598 0.014 0.0572 0.0013 0.0241 0.0006 0.3989 0.0004 472 4 476 11 13 19 244 0.0750 0.0008 0.599 0.018 0.0579 0.0016 0.0125 0.0003 0.946 0.005 466 5 477 14 14 43 568 0.0737 0.0007 0.577 0.011 0.0568 0.0010 0.0215 0.0005 0.4411 0.0003 458 5 462 9 15 73 995 0.0736 0.0007 0.578 0.009 0.0570 0.0009 0.0165 0.0004 0.450 0.003 458 4 463 7 16 46 762 0.0587 0.0006 0.446 0.007 0.0551 0.0008 0.0100 0.0003 0.7160 0.0005 368 5 374 6 17 60 821 0.0749 0.0007 0.589 0.010 0.0570 0.0008 0.0160 0.0005 0.317 0.002 465 4 470 8 18 44 564 0.0748 0.0007 0.583 0.012 0.0565 0.0011 0.0219 0.0006 0.490 0.001 465 4 466 9 20 28 405 0.0735 0.0007 0.576 0.013 0.0569 0.0013 0.0210 0.0005 0.0628 0.0000 457 4 462 11 21 35 471 0.0742 0.0007 0.585 0.011 0.0572 0.0009 0.0217 0.0005 0.3703 0.0004 462 4 468 8 22 150 1821 0.0732 0.0007 0.596 0.009 0.0590 0.0007 0.0210 0.0005 0.819 0.002 456 5 475 7 23 122 1841 0.0647 0.0006 0.498 0.007 0.0559 0.0007 0.0160 0.0004 0.497 0.002 404 4 411 6 24 95 1320 0.0722 0.0007 0.560 0.008 0.0563 0.0007 0.0226 0.0006 0.3214 0.0007 449 4 452 7 -
Song S G, Niu Y L, Zhang L F, et al.Time constraints on orogenesis from oceanic subduction to continental subduction, collision, and exhumation:an example from North Qilian and North Qaidam HP-UHP belts[J].Acta Petrologica Sinica, 2009, 25:2067-2077. http://en.cnki.com.cn/Article_en/CJFDTOTAL-YSXB200909003.htm
许志琴, 徐惠芬, 张建新, 等.北祁连走廊南山加里东俯冲杂岩增生地体及其动力学[J].地质学报, 1994, 68:1-15. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400804917 张雪亭, 杨生德, 杨站军, 等.青海省板块构造研究1:100万青海省大地构造说明书[M].北京:地质出版社, 2007:1-221. Xiao W J, Windley B F, Yong Y, et al.Early Paleozoic to Devonian multiple-accretionary model for the Qilian Shan, NW China[J].Journal of Asian Earth Sciences, 2009, 35(3/4):323-333. doi: 10.1016-j.jseaes.2008.10.001/
Yan Z, Aitchison J C, Fu C L, et al.Hualong Complex, South Qilian terrane:U-Pb and Lu-Hf constraints on Neoproterozoic micro-continental fragments accreted to the northern Proto-Tethyan margin[J].Precambrian Research, 2015, 266:65-85.
Xia L Q, Li X M, Yu J Y, et al.Mid-late Neoproterozoic to Early Paleozoic volcanism and tectonic evolution of the Qilianshan, NW China[J].Geo.Res.J., 2016, 9/12:1-41. https://www.sciencedirect.com/science/article/abs/pii/S2214242816300092
Zhang J X, Yu S Y, Mattinson C G..Early Paleozoic polyphase metamorphism in northern Tibet, China[J].Gondwana Research, 2017, 41:267-289. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=7367864fc4cf6d175a90dfd73ec90132
Fu C L, Yan Z, Guo X Q, et al.Assembly and dispersal history of continental blocks within the Altun-Qilian-North Qaidam mountain belt, NW China[J].International Geology Review, 2018:1-24. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1080/00206814.2018.1428831
Lu S N, Li H K, Zhang C L, et al.Geological and geo-chronologicai evidence for the Precambrian evolution of the tarim Craton and surrounding continental fragments[J].Precambrian Research, 2008, 160:94-107. http://www.sciencedirect.com/science/article/pii/S0301926807001581
陈宣华, 邵兆刚, 熊小松, 等.祁连山带断裂构造体系、深部结构与构造演化[J].中国地质, 2019, 46(5):995-1020. Wang Q Y, Pan Y M, Chen N S.Proterozoic polymetamorphism in the Quanji Block, northwestern china:evidence from microtextures, garnet compositions and monazite CHIME ages[J].Journal of Asian Earth Sciences, 2009, 34:686-698. http://www.sciencedirect.com/science/article/pii/S1367912008001600
Song S G, Niu Y L, Zhang L F, et al.Tectonic evolution of early Paleozoic HP metamorphic rocks in the North Qilian Mountains, NW China:New perspectives[J].Journal of Asian Earth Sciences, 2009, 35:334-353. doi: 10.1016-j.jseaes.2008.11.005/
Xu Y J, Du Y S, Cawood P A, et al.Detrital zircon record of continental collision:assembly of the Qilian Orogen, China[J].Sedimentary Geology, 2010, 230(1/2):35-45. doi: 10.1016-j.sedgeo.2010.06.020/
Zhang Z W, Li W Y, Gao Y B, et al.Sulfide mineralization associated with arc magmatism in the Qilian Block, western china:zircon U-Pb age and Sr-Nd-Os-S isotope constraints from the Yulonggou and Yaqu gabbroic intrusions[J].Mineralium Deposita, 2014, 49(2):279-292. doi: 10.1007/s00126-013-0488-x
李文渊.古亚洲洋与古特提斯洋关系初探[J].岩石学报, 2018, 34(8):2201-2210. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201808001 Song S G, Niu Y L, Su L, et al.Tectonics of the North Qilian orogeny, NW China[J].Gondwana Research, 2013, 23:1378-1401. http://www.sciencedirect.com/science/article/pii/S1342937X12000445
许志琴, 张建新, 许惠芬, 等.中国主要大陆山链韧性剪切带及动力学[M].北京:地质出版社, 1997:1-99. 冯益民, 何世平.祁连山大地构造与造山作用[M].北京:地质出版社, 1996:1-101. 张建新, 许志琴, 陈文, 等.北祁连中段增生杂岩/火山弧的时代探讨[J].岩石矿物杂志, 1997, 16(2):112-119. http://www.cqvip.com/Main/Detail.aspx?id=2511814 夏林圻, 夏祖春, 徐学义.北祁连山奥陶纪弧后盆地火山岩浆成因[J].中国地质, 2003, 30(1):48-60. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi200301006 Song S G, Yang L M, Zhang Y Q, et al.Qi-Qin accretionary belt in Central China Orogen:Accretion by trench jam of oceanic plateau and formation of intra-oceanic arc in the Early Paleozoic Qin-Qi-Kun Ocean[J].Science Bulletin, 2017, 62:1035-1038. http://www.cqvip.com/QK/86894X/201715/673136310.html
吴才来, 徐学义, 高前明, 等.北祁连早古生代花岗质岩浆作用及构造演化[J].岩石学报, 2010, 26(4):1027-1044. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201004003 宋述光, 张立飞, Niu Y, 等.青藏高原北缘早古生代板块构造演化和大陆深俯冲[J].地质通报, 2004, 23(9/10):918-925. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=200409164&flag=1 丁仨平, 裴先治, 李佐臣, 等.祁连造山带中段雾宿山群的构造属性探讨[J].中国地质, 2008, 35(4):577-589. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi200804002 刘志武, 王崇礼, 石小虎.南祁连党河南山花岗岩类特征及其构造环境[J].现代地质, 2006, 20(4):545-554. http://d.old.wanfangdata.com.cn/Periodical/xddz200604004 师江朋, 霍腾飞, 来强, 等.南祁连北缘东段早志留世刚察大寺花岗岩的成因:锆石LA-ICP-MS U-Pb年代学和岩石地球化学制约[J].地球学报, 2015, 36(6):781-789. http://d.old.wanfangdata.com.cn/Periodical/dqxb201506013 罗志文, 张志诚, 李建锋, 等.中南祁连西缘肃北三个洼塘地区古生代两类花岗质侵入岩年代学及其地质意义[J].岩石学报, 2015, 31(1):176-188. http://d.old.wanfangdata.com.cn/Periodical/ysxb98201501013 郭周平, 李文渊, 张照伟, 等.南祁连化隆地区鲁满山花岗岩的岩石成因:地球化学、锆石U-Pb年代学及Hf同位素约束[J].中国地质, 2015, 42(4):864-880. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi201504006 Yang H, Zhang H F, Luo B J, et al.Early Pale0zoic intrusive rocks from the eastern Qilian Orogen NE Tibetan plateau:petrogenesis and tectonic significance[J].Lithos, 2015, 224/225:13-31.
李春昱, 刘仰文, 朱宝清, 等.秦岭及祁连山构造发展史[C]//国际交流地质学术论文集(1).北京: 地质出版社, 1978: 174-187. 李五福, 张新远, 王春涛, 等.祁连山哈拉湖地区奥陶纪岛弧火山岩及其构造意义[J].地质通报, 2019, (388):1287-1296. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20190804&flag=1 左国朝, 张淑玲, 程建生, 等.祁连地区蛇绿岩带划分及其构造意义[C]//张旗.蛇绿岩与地球动力学研究.北京: 地质出版社, 1996: 129-134. Wan Y S, Zhang J X, Yang J S, et al.Geochemistry of highgrade metamorphic rocks of the North Qaidam Mountains and their geological significance[J].Journal of Asian Earth Sciences, 2006, 28:174-184. http://www.sciencedirect.com/science/article/pii/S1367912006000307
何世平, 李荣社, 王超, 等.祁连山西段甘肃肃北地区北大河岩群片麻状斜长角闪岩的形成时代[J].地质通报, 2010, 29(9):1275-1280. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20100903&flag=1 Song S G, Niu Y L, Su L, et al.Continental orogenesis from ocean subduction, continent collision/subduction, to orogen collapse, and orogen recycling:The example of the North Qaidam UHPM belt, NW China[J].Earth-Science Reviews, 2014, 129:59-84. http://www.sciencedirect.com/science/article/pii/S0012825213002055
Middlemost E A K.Naming materials in the magma/igneous rock system[J].Earth Science Reviews, 1994, 37(3/4):215-224. http://www.sciencedirect.com/science/article/pii/0012825294900299
Peccerillo R, Taylor S R.Geochemistry of eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey[J].Contributions to Mineralogy & Petrology, 1976, 58(1):63-81. doi: 10.1007-BF00384745/
Irvine T N, Baragar W R A.A guide to the chemical classification ofthe common volcanic rocks[J].Canadian Journal of Earth Sciences, 1971, 8(5):523-548. doi: 10.1007%2FBF02596878
Boynton W V.Cosmochemistry of the Rare Earth Elements:Meteoric Studies[J].Rare Earth Element Geochemistry, 1984, 2:63-114. https://www.sciencedirect.com/science/article/pii/B9780444421487500083
Sun S S, Mcdonough W F.Chemical and lsotopic Systematics of Oceanic Basalts:Implications for Mantle Composition and Processes[J].Geological Society London Special Publications, 1989, 42(1):313-345. https://www.researchgate.net/publication/231575101_Chemical_and_isotopic_systematics_of_oceanic_basalts_Implications_for_mantle_composition_and_processes
冯益民.祁连山大地构造与造山作用[M].北京:地质出版社, 1996. 夏林圻, 夏祖春.北祁连山古洋壳火山作用基本特征[C]//中国矿物学岩石学地球化学研究新进展学术交流会, 1994. 夏林圻, 夏祖春, 徐学义.北祁连山奥陶纪弧后盆地火山岩浆成因[J].中国地质, 2003, 30(1):48-60. http://d.old.wanfangdata.com.cn/Periodical/zgdizhi200301006 张建新, 孟繁聪.北祁连和北阿尔金硬柱石榴辉岩:冷洋壳俯冲作用的证据[J].科学通报, 2006, 51(14):1683-1688. http://d.wanfangdata.com.cn/Periodical/kxtb200614011 张建新, 孟繁聪, Mattinson C G.南阿尔金-柴北缘HP/UHP变质带研究进展、问题及挑战[J].高校地质学报, 2007, 13(3):526-545. 张建新, 许志琴.北祁连中段加里东俯冲-增生杂岩/火山弧带及其变形特征[J].地球学报, 1995, (2):36-46. http://www.cnki.com.cn/Article/CJFDTotal-DQXB502.003.htm 董顺利, 李忠, 高剑, 等.阿尔金-祁连-昆仑造山带早古生代构造格架及结晶岩年代研究进展[J].地质论评.2013, 59(4):731-747. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=HYC201405070000001439 Zhang Y Q, Song S G, Yang L M, et al.Basalts and picrites from a plume-type ophiolite in the South Qilian Accretionary Belt, Qilian Orogen:Accretion of a Cambrian Oceanic Plateau?[J].Lithos, 2017, 278:97-110.
宋述光, 苏犁, 张贵宾, 等.秦-祁-昆地区典型高压变质带与蛇绿混杂岩带关系研究.北京大学, 2015: 1-505.