Geochemical characteristics of dacite from the Bayinsuhetu deposit in South Mongolia and its tectonic implications
-
摘要:
在南蒙古巴音苏赫图矿区古生代地层中新发现一套火山岩,其Rb-Sr同位素年龄为409±19Ma,与前人划分的地层年代(志留纪-泥盆纪)较吻合。英安岩的SiO2含量为63.23%~73.36%,Al2O3含量为10.25%~15.56%。英安岩轻稀土元素富集,重稀土元素亏损,具有较强的负Eu异常。岩石学和地球化学综合研究表明,该套火山岩基本同期同源演化,且岩性变化较小,以流纹岩-英安岩类为主。根据钾含量,该套火山岩为中钾钙碱性系列。英安岩落在靠近盆岭构造域的位置,显示出略亏损的特征,与原始地幔PM端元接近。推测该套火山岩形成于古亚洲洋闭合过程中的岛弧环境,表明火山岩形成时古亚洲洋可能还没有闭合。
Abstract:A set of volcanic rocks were recently found in Palaeozoic strata in the Bayinsukhtuore district of South Mongolia. The RbSr isotopic age of the rocks is 409±19Ma, which is consistent with the chronological age (Silurian-Devonian) obtained from previous work. The SiO2 content of dacite is 63.23%~73.36%, while the Al2O3 content is 10.25%~15.56%. The dacite is LREE enrichment, HREE loss, and has strong Eu negative anomaly. Based on petrology and geochemistry, it is held that this set of volcanic rocks evolved in the same period and came from the same source. The lithologic character of the rocks is mainly of dellenite-dacite type, with little change. According to the Kalium analysis, the volcanic rocks were derived from the potassium calc-alkaline series. It is suggested that the melted mantle source area was reformed by the subducting plate. During this process the original magma was generated from the liquid mantle of the subducting zone. At the same time crustal contamination happened. The magmatic evolution would be the main reason for generating the combined set of the volcanic rocks. It also implies that this set of volcanic rocks was formed in the island arc during the closing of pal-Asia, which indicates that the closing of pal-Asia was later than the formation of this set of volcanic rocks.
-
Keywords:
- South Mongolia /
- volcanic rock /
- geochemical characteristics /
- mantle /
- origin
-
地质数据是地质工作中形成的重要信息资源,多年累积形成的海量地质数据信息为国民经济和社会发展提供了重要支撑[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 中蒙边境钨-钼矿化集中区区域地质简图(据参考文献[31]修改)
Figure 1. Simplified regional geological map of the W-Mo mineralized area along the Sino-Mongolian border
图 7 Y-Sr/Y关系图解[37]
Figure 7. Sr/Y versus Y diagram
图 9 研究区火山岩Isr-εNd(t)判别图(底图据参考文献[17])
DM—亏损地幔;EMⅠ、EMⅡ—富集地幔;PM—普通地幔;MORB—洋中脊玄武岩;OIB—洋岛玄武岩
Figure 9. Isr-εNd(t) diagram for volcanic rocks in the study area
表 1 研究区火山岩主量、微量和稀土元素分析结果
Table 1 The analyzed data of major, trace and rare earth elements of volcanic rocks in the study area
样品 HS02 HS03 HS04 HS11 HS12 HS13 HS14 HS15 SiO2 71.39 69.41 67.83 69.32 69.63 72.97 73.36 63.23 TiO2 0.43 0.54 0.58 0.54 0.47 0.48 0.33 0.53 Al2O3 14.47 13.56 14 12.85 15.56 12.63 10.25 12.94 A/CNK 1.35 1.18 1.5 1.09 1.65 1.49 0.94 1.52 A/NK 1.53 1.83 1.88 1.82 1.95 1.91 1.45 1.85 MnO 0.03 0.05 0.04 0.06 0.03 0.04 0.04 0.05 MgO 1.39 1.72 1.97 2.2 1.69 1.8 1.22 1.91 CaO 0.67 2.23 1.04 2.59 0.81 1.02 2.13 0.82 Na2O 4.56 3.05 2.96 3.14 2.9 2.57 2.33 3.3 K2O 1.82 2.2 2.4 1.75 2.97 2.2 2.99 1.46 Na2O+K2O 6.38 5.25 5.36 4.89 5.87 4.77 5.32 4.76 P2O5 0.13 0.22 0.23 0.2 0.14 0.2 0.12 0.21 烧失量 1.9 3.56 2.82 3.96 2.96 2.58 2.86 0.76 合计 99.99 100 100 100 100 100 100 100.01 Ba 418 424.78 400.54 313.54 613.7 389.84 616.79 372.47 Cr 42.41 60.02 63.69 64.15 41.89 47.67 26.68 62.15 Ni 18.52 15.45 20.81 29.47 18.37 19.43 10.44 23.39 Sr 358.02 172.12 236.2 271.6 348.43 241.36 139.79 260.96 V 49.25 58.62 69.42 65.65 48.37 60.54 39.4 69.67 Zr 180.59 203.91 207.43 152 211.53 155.09 180.69 184.73 TFe2O3 3.2 3.94 4.28 4.31 3.51 3.91 2.59 4.11 Ga 16.68 16.9 16.09 15.35 20.17 14.37 14.03 17.8 Rb* 75.09 89.44 99.98 80.53 116.27 89.61 117.44 61.51 Sr 356.36 184.82 220.06 308.6 303.38 240.04 171.03 275.75 Y 47.26 38.31 38.37 30.27 49.28 32.4 24.44 32.94 Zr 180.24 244.92 197.59 170.66 201.97 169.89 210.03 221.7 Nb 11.33 10.7 10.53 9.91 13 8.1 8.18 10.47 Cs 5.98 6.35 9.82 5.16 6.89 6.68 3.36 4.86 Ba* 431.13 483.26 422.32 358.66 561.68 403.22 752.71 461.75 La 39.42 37.43 41.19 37.45 52.03 36.32 31.96 38.86 Ce 86.33 60.06 75.73 70.49 87.18 67.72 58.45 74.2 Pr 10.34 9.51 11.06 9.2 13.54 9.09 7.41 9.64 Nd 41.17 37.99 43.32 37 53.71 36.57 27.99 39.17 Sm 8.69 7.56 8.64 7.24 10.86 7.32 5.29 7.75 Eu 1.08 1.52 1.76 1.41 1.33 1.58 1.05 1.62 Gd 8.45 6.96 7.99 6.58 9.76 6.82 4.72 7.15 Tb 1.49 1.15 1.27 1.03 1.61 1.09 0.74 1.12 Dy 9.06 7.03 7.43 5.84 9.56 6.28 4.35 6.44 Ho 1.89 1.51 1.53 1.2 2 1.3 0.89 1.28 Er 5.1 4.18 4.16 3.25 5.51 3.44 2.57 3.63 Tm 0.77 0.64 0.62 0.49 0.84 0.51 0.4 0.54 Yb 5.03 4.19 4.08 3.17 5.57 3.27 2.58 3.49 Lu 0.77 0.64 0.62 0.49 0.86 0.51 0.39 0.53 Hf 5.82 7.14 5.86 4.88 6.74 4.9 6.22 6.56 Ta* 1.11 0.88 0.83 0.75 1.29 0.63 0.65 1.03 Tl* 0.51 0.58 0.7 0.48 0.7 0.51 0.68 0.38 Pb 19.08 29.49 18.45 14.33 18.18 23.95 18.23 17.94 Bi* 0.35 0.22 0.15 0.29 0.46 0.19 0.1 0.21 Th 18.43 14.42 12.97 11.13 20.65 10.14 13.03 12.67 U 3.98 2.83 3.2 2.04 3.55 2.23 2.11 2.2 LREE 187.04 154.06 181.69 162.78 218.65 158.6 132.15 171.24 HREE 32.55 26.31 27.69 22.05 35.71 23.22 16.65 24.19 ∑REE 219.59 180.37 209.38 184.83 254.36 181.82 148.79 195.42 LREE/HREE 5.75 5.85 6.56 7.38 6.12 6.83 7.94 7.08 注:主量与微量元素在中国科学院地质与地球物理研究所分别采用ICP-AES和ICP-MS分析;球粒陨石标准化值据参考文献[35];主量元素含量单位为%,微量和稀土元素为10-6 表 2 研究区火山岩Sr-Nd同位素分析结果
Table 2 Sr and Nd isotope composition of volcanic rocks in study area
样号 HS02 HS03 HS04 HS11 HS13 HS15 Rb/10-6 85.28 90.04 110.1 69.31 87.67 59.3 Sr/10-6 364 173.3 246.1 278.9 243.8 275.8 87Rb/86Sr 0.6783 1.506 1.296 0.7196 1.041 0.6226 87Sr/86Sr(2σ) 0.711279±15 0.716905±11 0.715226±14 0.712401±14 0.713748±9 0.711880±12 Sm(10-6) 7.654 6.898 8.177 6.425 7.105 7.151 Nd(10-6) 34.09 32.7 38.75 31.44 33.69 33.91 147Sm/144Nd 0.1359 0.1277 0.1277 0.1237 0.1277 0.1277 143Nd/144Nd(2σ) 0.512362±6 0.512323±7 0.512334±13 0.512332±8 0.512338±7 0.512329±6 TDM 1545 1467 1448 1386 1439 1454 Isr 0.707329 0.708135 0.707678 0.708209 0.707682 0.708254 εSr(t) 47.02 58.48 51.99 59.53 52.04 60.16 εNd(t) -2.2 -2.5 -2.3 -2.2 -2.2 -2.4 fSm/Nd -0.31 -0.35 -0.35 -0.37 -0.35 -0.35 T2DM 1333 1360 1343 1328 1335 1349 注:2σ为测量误差(10-6)。研究区火山岩(87Sr /86Sr)i和εNd(t)由409Ma推算;所用参数143Nd /144Nd(CHUR)0= 0.512638,147Sm /144Nd(CHUR)0= 0.1967,λRb87= 1.42×10-12a-1,λSm147= 6.57×10-12a-1 -
Dobrelsov N L, Berzin N A, Buslov M M. Opening and tectonic evolution of the Paleo-Asian Ocean. International[J]. Geo. Rev., 1995, 37:335-360. doi: 10.1080/00206819509465407
Yin A, Mark H. A Phanerozoic palinspastic reconstruction of China and its neighboring regions[C]//Yin A, Mark Harrison. The tectonic evolution of Asia. Cambridge University Press, 1996: 422-485.
董树文, 吴锡浩, 吴珍汉, 等.论东亚大陆的构造壳变——燕山运动的全球意义[J].地质论评, 2000, 46(1):8-13. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzlp200001001&dbname=CJFD&dbcode=CJFQ 张招崇, 王永强.冀北印支期碱性岩浆活动及其地球动力学意义[J].矿物岩石地球化学通报, 1997, 16(4):6-9. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=kydh704.001&dbname=CJFD&dbcode=CJFQ 张理刚.东亚岩石圈块体地质——上地幔、基底和花岗岩同位素地球化学及其动力学[M].北京:科学出版社, 1995:1-249. 陈义贤, 陈文寄.辽西及邻区中生代火山岩——年代学地球化学和构造背景[M].北京:地震出版社, 1995:1-279. 许保良, 王式洸, 韩宝福, 等.富集型与亏损型A型花岗岩——以燕山地区和乌伦古河地区岩石为例[J].北京大学学报(自然科学版), 1998, 34(2/3):352-362. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=bjdz8z1.028&dbname=CJFD&dbcode=CJFQ 李伍平, 路凤香, 孙善平, 等.北京西山东岭台组(J3d)火山岩的成因及其构造环境探讨[J].岩石学报, 2000, 16(3):345-352. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200003004&dbname=CJFD&dbcode=CJFQ 李伍平, 路凤香, 李献华, 等.北京西山髫髻山组火山岩的地球化学特征与岩浆起源[J].岩石矿物学杂志, 2001, 2:123-133. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yskw200102002&dbname=CJFD&dbcode=CJFQ 李伍平, 路凤香, 李献华, 等.辽西地区义县组火山岩的起源及壳幔相互作用[J].矿物岩石, 2001, 21(4):1-6. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=kwys200104000&dbname=CJFD&dbcode=CJFQ 李伍平, 李献华, 路凤香.辽西中侏罗世高Sr低Y型火山岩的成因及其地质意义[J].岩石学报, 2001, 17(4):523-532. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200104002&dbname=CJFD&dbcode=CJFQ 李伍平, 路凤香, 李献华.北京西山髫髻组粗安岩的成因及其地质意义[J].岩石矿物学杂志, 2001, 20(3):247-254. http://www.oalib.com/paper/4337617 李伍平, 李献华, 路凤香, 等.辽西早白垩世义县组火山岩的地质特征及其构造背景[J].岩石学报, 2002, 2:193-204. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200202006&dbname=CJFD&dbcode=CJFQ 李伍平, 李献华.燕山造山带中段中晚侏罗世中酸性火山岩的成因及其意义[J].岩石学报, 2004, 3:501-510. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200403014&dbname=CJFD&dbcode=CJFQ 张旗, 王焰, 王元龙.燕山期中国东部高原下地壳组成初探:埃达克质岩Sr、Nd同位素制约[J].岩石学报, 2001, 17(4):505-513. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200104000&dbname=CJFD&dbcode=CJFQ 韩宝福, 何国琦, 王式洸.后碰撞幔源岩浆活动、底垫作用及准噶尔盆地基底的性质[J].中国科学(D辑), 1999, 1:16-21. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jdxk199901002&dbname=CJFD&dbcode=CJFQ 洪大卫, 王式洸, 谢锡林, 等.兴蒙造山带正ε(Nd, t)值花岗岩的成因和大陆地壳生长[J].地学前缘, 2000, 7(2):441-456. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dxqy200002016&dbname=CJFD&dbcode=CJFQ 阎国翰, 牟保磊, 许保良, 等.燕辽-阴山三叠纪碱性侵入岩年代学和Sr, Nd, Pb同位素特征及意义[J].中国科学(D辑), 2000, 30(4):380-387. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jdxk200004005&dbname=CJFD&dbcode=CJFQ Jahn B M, Wu F Y, Capdevila R, et al. Higly evolved juvenile granites with tetrad REE patterns:the Woduhe and Baerzhe granites from the Great Xing'an Mountains in NE China[J]. Lithos, 2001, 59(4):171-198. doi: 10.1016/S0024-4937(01)00066-4
Wu F Y, Sun D Y, Li H, et al. A-type granites in northeastern China:age and geochemical constraints on their petrogenesis[J]. Chemical Geology, 2002, 187(1/2):143-173. http://www.sciencedirect.com/science/article/pii/S0009254102000189
郭锋, 范蔚茗, 李超文, 等.早古生代古亚洲洋俯冲作用:来自内蒙古大石寨玄武岩的年代学与地球化学证据[J].中国科学(D辑), 2009, 5:569-579. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jdxk200905004&dbname=CJFD&dbcode=CJFQ 刘建峰. 内蒙古林西-东乌旗地区晚古生代岩浆作用及其对区域构造演化的制约[D]. 吉林大学, 2009. http://cdmd.cnki.com.cn/Article/CDMD-10183-2009093598.htm Blight J H S, Petterson M G, Crowley Q G, et al. The Oyut Ulaan Volcanic Group:stratigraphy, magmatic evolution and timing of Carboniferous arc development in SE Mongolia[J]. Journal of the Geological Society, 2010, 167(3):491-509. doi: 10.1144/0016-76492009-094
Badarch G, Cunningham W D, Windley B F. A New Terrane Subdivision for Mongolia; Impications for the Phanerozoic Crustal Growth of Central Asia[J]. Journa of Asian Earth Sciences, 2002, 21:87-110. doi: 10.1016/S1367-9120(02)00017-2
Ren X D, Kiosses W B, Alexander S M. Regulation of the small GTP-binding protein Rho by cell adhesion and the cytoskeleton[J]. The EMBO Journal, 1999, 18:578-585. doi: 10.1093/emboj/18.3.578
Zeng Q, Liu J, Zhang Z, et al. Geology and geochronology of the Xilamulun molybdenum metallogenic belt in eastern Inner Mongolia, China[J]. International Journal of Earth Sciences, 2011, 100:1791-1809. doi: 10.1007/s00531-010-0617-z
林强, 葛文春, 曹林, 等.大兴安岭中生代双峰式火山岩的地球化学特征[J].地球化学, 2003, 32(3):208-222. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqhx200303001&dbname=CJFD&dbcode=CJFQ 聂凤军, 张万益, 杜安道, 等.内蒙古朝不楞矽卡岩型铁多金属矿床辉钼矿铼-锇同位素年龄及地质意义[J].地球学报, 2007, 28:315-323. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqxb200704000&dbname=CJFD&dbcode=CJFQ 张德全.大兴安岭南段不同构造环境中的两类花岗岩[J].岩石矿物学杂志, 1993, 1:1-11. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yskw199301000&dbname=CJFD&dbcode=CJFQ 周建波, 张兴洲, 马志红, 等.中国东北地区的构造格局与盆地演化[J].石油与天然气地质, 2009, 30:530-538. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=syyt200905002&dbname=CJFD&dbcode=CJFQ 聂凤军.中蒙边境中东段金属矿床成矿规律和找矿方向[M].北京:地质出版社, 2007. Chen F K, Hegner E, Todt W. Zircon Ages, Nd Isotopic and Chemical Compositions of Orthogneisses from the Black Forest, Germany:Evidence for a Cambrian Mag-matic Arc[J]. International Journal of Earth Sciences, 2000, 88(4):791-802. doi: 10.1007/s005310050306
Gao S, Rudnick R L, Yuan H L, et al. Recycling Lower Continental Crust in the North China Craton[J]. Nature, 2004, 432:892-897. doi: 10.1038/nature03162
Chen F K, Siebel W, Satir M, et al. Geochronoloy of the Karadere Basement (NW Turkey) and Implications for the Geological Evolution of the Istanbul Zone[J].International Journal of Earth Sciences, 2002, 91:469-481. doi: 10.1007/s00531-001-0239-6
Boynton, Rose W, Dunn, et al. Manual of Ambulatory PEDIATRICS[J]. Ajn American Journal of Nursing, 1984, 84(12):1486. https://www.libreriauniverso.it/articolo/41387/0781788749/Boynton-Manual-of-Ambulatory-Pediatrics
Sun S S, McDonough W S, Chemical and isotopic systematics of oceanic basalts:implications for mantle composition and processes[J]. Geological Society, London, Special Publications, 1989, 42(1):313-345. doi: 10.1144/GSL.SP.1989.042.01.19
Drummond M S, Defant M J. A model for Trondhjemite-Tonalite-Dacite Genesis and crustal growth via slab melting:Archean to modern comparisons[J]. Journal of Geophysical Research Solid Earth, 1990, 95(B13):21503-21521. doi: 10.1029/JB095iB13p21503
吴华英, 张连昌, 周新华, 等.大兴安岭中段晚中生代安山岩年代学和地球化学特征及成因分析[J].岩石学报, 2008, 6:1339-1352. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200806018&dbname=CJFD&dbcode=CJFQ Manya S, Makenya A H M, Nakamurab E. The geochemistry of High-Mg andesite and associated adakitic rocks in the musoma-mara greenstone belt, northern tanzania:possible evidence for neoarchaean ridge subduction[J]. Precambrian Research, 2007, 159:241-259. doi: 10.1016/j.precamres.2007.07.002
Shinio R. Geochemistry of High-Mg andesite and the tectonic evidence of the okinawa trough-ryukyu arc system[J]. Chemical Geology, 1999, 157:69-88. doi: 10.1016/S0009-2541(98)00199-5
Sorokin A A, Sorokin A P, Ponomarchuk V A, et al. The Age and Geochemistry of Volcanic Rocks on the Eastern Flank of the Umlekan-Ogodzha Volcanoplutonic Belt (Amur Region)[J]. Russian Geology and Geophysics, 2010, 51:369-379. doi: 10.1016/j.rgg.2010.03.004
邵济安, 路凤香, 李伍平.辽西中生代陆内底侵作用背景下形成的安山岩[J].岩石学报, 2007, 23(4):701-708. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200704003&dbname=CJFD&dbcode=CJFQ 路凤香, 桑隆康.岩石学[M].北京:地质出版社, 2002:15-26. 李承东, 张福勤, 苗来成, 等.吉林色洛河晚二叠世高镁安山岩SHRIMP锆石年代学及其地球化学特征[J].岩石学报, 2007, 23(4):767-776. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb200704009&dbname=CJFD&dbcode=CJFQ