Geochemistry of the Yaochong granite in Dabie orogenic belt and its geological implications
-
摘要:
姚冲岩体位于秦岭-大别山造山带上，岩体的主要岩性为黑云母花岗岩和花岗斑岩。岩石地球化学分析结果显示，姚冲岩体SiO2含量为71.28%~77.24%，K2O含量为3.81%~5.28%，Na2O含量为3.76%~5.03%，Al2O3含量为11.61%~14.6%，具有高钾钙碱性和准铝质-过铝质 (A/CNK=0.96~1.03) 特征。岩石稀土元素总量较低 (11.82×10-6~123.65×10-6)，轻稀土元素相对富集，重稀土元素相对亏损，具有Eu负异常，岩体富集Rb、U、Hf和Y，亏损Ba、Nb、Ta和Ti，显示出I型花岗岩的特点。锆石Hf同位素分析结果表明，姚冲岩体的εHf(t) 值为-30~-22.8，位于地幔演化线之下，tDM2值变化于2.21~2.60Ga之间，指示其源于扬子板块北缘的古老地壳，其组成类似于扬子板块北缘TTG型岩浆岩。结合区域地质背景，姚冲岩体形成于扬子陆块与华北陆块碰撞造山后的陆内伸展环境，对应的地球动力学背景为晚侏罗世—早白垩世地壳从挤压收缩向区域性伸展的构造体制大转换过程。
Abstract:Located in the Qinling-Dabie orogenic belt, the Yaochong granite body mainly consists of biotite granite and granite porphyry. According to the geochemical analysis, Yaochong granites have the following data: SiO2=71.28%~77.24%, K2O=3.81%~5.28%, Na2O=3.76%~5.03%, Al2O3=11.61%~14.6%, and A/CNK=0.96~1.03. The rocks are enriched in LREE, depleted in HREE with negative Eu anomalies, enriched in Rb, U, Hf, Y and depleted in Ba, Nb, Ta, Ti, which indicates that Yaochong granites belong to metaluminous to peraluminous, high-K calc-alkaline I-type granite. According to zircon Hf isotope analysis, the εHf(t) values of the granite range from-30 to-22.8, lying below the depleted mantle line. The tDM2 values range from 2.21Ga to 2.60Ga, which indicates that the granite resulted from the lower crust of Yangtze block. The components of the granite porphyry are similar to those of Neoproterozoic TTG magmatic rocks. Yaochong granites were emplaced after the major collisional period of the North China block and Yangtze block, in a post-collisional tectonic setting under the regional extension mechanism. The geodynamic setting was the transformation of the tectonic regime of the late Jurassic-Early Cretaceous period.
-
Keywords:
- Dabie /
- Yaochong /
- granite /
- geochemistry /
- geological implications
-
“玻璃地球”（Glass Earth）计划最初由澳大利亚学者提出，指应用地质、地球物理勘探和三维可视化技术，使大陆表层1000m以浅“像玻璃一样透明”，并通过计算机网络技术为地质工作者及社会公众提供地学信息分析及决策支持服务[1]。它主要用于描述地壳浅层地质结构、成分及其空间拓扑关系，釆集、管理和处理基础地质调查、矿产地质勘查、矿产资源开釆、地下水管理、矿权管理与生产安全监控、工程地质、灾害地质勘查等信息。加拿大和法国相继提出了类似的计划，把目标提高到地下3000m以浅[2-3]。随后，其他欧美国家纷纷响应并制定了相应的规划，开展了关键技术研究[1, 4]，其中深部勘探技术和三维地质建模技术得到长足发展[5-7]。在资源与环境的双重压力下，让地球深部“透明化”已成为越来越多国家关注的热点。
中国地质找矿和地质灾害监测预警所需的深部地质信息还比较匮乏，迫切需要基础地质工作者转变工作模式。近年来，随着计算机三维技术的迅速发展，三维地质填图和建模作为描绘地质信息的技术手段也逐渐成熟[8]。因此，中国基础地质调查“从传统走向现代、从单一走向综合、从二维走向三维”的理念呼之欲出。2012年，中国地质调查局正式启动三维地质调查试点项目，设立了一个计划项目，在地壳结构、含油气盆地、成矿带、矿集区、经济区、技术方法等领域设置了14个工作项目。《相山火山盆地三维地质调查》是其中的一个工作项目，重点探索矿集区1：5万甚至更大比例尺的三维地质调查方法。本文简要介绍其中三维建模的实践过程，以及展望与建议。
1. 三维地质调查与建模发展历程
20世纪70年代末，美国率先启动的大陆反射地震探测计划（COCORP），在30个州采集超过10000km长的地震反射剖面，以揭示北美地壳深部精细结构[9]。受该计划影响和人类对矿产资源需求的扩大，80年代主要发达国家相继开展了大量深部地质探测计划及三维地质建模工作，如国际岩石圈委员会全球地学断面计划（GGT）、欧洲地球探测计划（EUROPROBE）、法国地壳探测计划（ECORP）、德国大陆反射地震计划（DEKORP）、英国反射地震计划（BIRPS）、瑞士地壳探测计划（NRP20）、意大利地壳探测计划（CROP）、加拿大岩石圈探测计划（LITHOPROBE）等[10-15]。90年代后，以澳大利亚为代表的西方国家发起了新一轮的三维地质调查计划，包括四维地球动力学计划（AGCRC）、玻璃地球计划（GlassEarth）、地球探测计划（AuScope）等[1]。一系列深部地质调查计划，极大地推动了三维地质建模技术的发展和三维地质建模软件的研发，如EarthVision、Micromine、DataMine、GOCAD、Surpac等现今国际主流建模软件都与这些计划密切相关。
1989年底首次召开了专门讨论3D地学建模问题的国际会议[16]。GOCAD软件创始人法国南希大学的Mallet教授分别于1989年、2004年提出了离散光滑插值技术（Discrete Smooth Interpolation，DSI）和（古地理）坐标转换技术（UVT），成为GOCAD、SKUA-GOCAD软件的核心技术[17-18]。众多国外学者基于地质填图数据建模[19-20]、地质剖面建模[21-26]、物探数据建模[27-29]、多源数据建模[30-31]等开展了大量研究，取得了丰硕的成果。
中国的三维地质调查与建模研究始于20世纪80年代，主要涉及地矿、石油、冶金、铁路等领域[32]。1985年以来，中国陆续在唐山、邢台、延庆-怀来盆地、东秦岭、燕山褶皱带等地开展了深地震反射探测。1998年，中国地震局实施了《长白山天池火山岩区岩浆系统、地壳结构的三维深地震测深研究》项目。2000年，中国科学院启动了《华北地区内部结构探测研究计划》项目。2001年，董树文等[33]完成了横贯大别山前陆褶皱带的深地震反射剖面探测，近年来实施的横穿造山带内部的深反射地震剖面揭示了造山带内的地壳结构。2008年，中国地质调查局先后启动了《长江中下游综合地球物理调查的立体地质填图应用试验》和《华南岩体形态圈定与研究》项目，初步尝试利用物探和钻探开展三维地质填图的方法组合技术[34]。2009年，中国地质科学院岩石圈中心在华北地区开展了一条550km的深反射剖面，南起怀来盆地，北至内蒙古二连盆地，揭示了华北北缘与兴蒙造山带的深部结构。21世纪以来，许多地矿与勘测部门在城市地质、工程地质、固体矿产及其可视化研究应用方面做了大量的尝试[35-38]，积累了一批数据和经验。2004年，上海市启动了三维城市地质调查工作，建立了三维基岩地质、第四纪地质、工程地质和水文地质结构模型，开启了地质工作社会化服务的新篇章[39]。
在三维技术领域方面，李德仁等[40]、李清泉等[41]、龚健雅等[42]提出了矢量与栅格集成的3D数据模型，指出空间实体由多种空间对象组成，这些对象可以由矢量或栅格数据表达。朱良峰等[43]提出了一种由工程钻孔数据构建三维地层模型的方法，能够将钻孔剖面图融入实际建模流程。2007年，张宝一等[44]提出了利用三维地质建模与可视化技术进行固体矿产储量估算的可行性方案，并提出了应用轮廓线进行三维矿体表面建模时尖灭与分支情况处理、带约束的三维矿体表面建模的解决方案。2013年，高阳等[45]利用国产3DMine矿业工程软件系统建立了广东石人嶂钨矿的地层模型、构造模型、矿体模型、地质工程模型等，实现了矿区深部数据的三维可视化。陈建平等[46-47]探讨了成矿带和含矿地质体三维模型建立并进行储量估算。
根据建模数据来源可将三维地质模型分为基于钻孔数据、基于地质剖面、基于地球物理测深数据、基于平面地质图、基于多源混合数据等类型[48-52]。受经济条件的制约，钻孔、勘探线剖面、物探资料等数据不易获取，数量有限，分布离散且不均匀。平面地质图成本相对低廉，容易获取且是能覆盖整个研究区的数据源，在一定程度上可缓解数据源带来的瓶颈[53-55]。因此，在其他地质数据匮乏的前提下，利用平面地质图构建区域三维地质模型是一种有效的解决方法[56]，既可以从整体上了解研究区域的地质构造，又可为后续增加钻孔、剖面数据后进一步精化模型作准备[30]。许多学者对以平面地质图为主要数据源的三维地质建模方法进行了研究[57-59]。笔者认为，在未开展地球物理测深工作、钻孔、坑道资料不多的地区，可以开展地质图三维建模，或称地质概念模型。
由于资金投入、地质构造复杂程度等因素的限制，很多三维数据模型尚处于探索和尝试阶段。目前，成熟的三维地质模拟技术只集中于简单层状地质体的三维结构与可视化表达分析上，而对结构复杂、物化属性分布不均匀的复杂地质体的三维可视化建模研究不够深入。当前，随着中国综合国力的增强，信息技术的快速进步，开展不同类型地区三维地质填图试点工作在客观上已经成为可能。
2011年3月，中国地质调查局在北京召开了三维地质填图学术交流与研讨会，论证了开展三维地质填图符合中国当今地质工作需求和未来发展趋势，是地质调查工作的一次重大创新，对实现地质调查由二维向三维的转变，显著提升地质工作的服务功能和支撑作用具有重大意义。
中国地质调查局从2012年开始实施三维地质调查项目。由笔者主持的《相山火山盆地三维地质调查》项目是在地表地质填图（1：5万）的基础上，综合运用已有的钻探与坑探资料，开展遥感解译、重磁三维反演、深孔综合测井和专题研究工作，根据深部目标地质体情况有针对性地部署MT、CSAMT测量工作，通过多参数交互解译确定目标地质体的三维空间展布特征，在GOCAD软件平台上构建不同类型的三维地质模型，为找矿勘探提供深部依据，并总结了矿集区三维地质调查方法[60-67]。
2. 三维地质建模方法
GOCAD是以工作流为核心的地质建模软件，三维建模、可视化、空间分析能力强大。它的建模是从面模型到实体模型，可以设计空间几何对象，也可以表现空间属性分布，具有真三维地质建模、快速建模、兼容性好等优点，拥有强大的地质统计分析功能。依托GOCAD软件平台，根据建模数据类型的不同，尝试了3种三维地质建模方法和技术。
2.1 数字地质填图建模
2.1.1 概念的提出
三维地质填图需要大量的物探、钻孔等数据支撑，成本昂贵。目前中国开展的数字地质调查积累的大量空间数据尚未得到进一步开发利用。因此，笔者尝试在GOCAD软件平台上，利用数字地质填图系统的野外路线PRB数据（P为地质点，R为分段路线，B为点间地质界线）直接进行浅表层三维地质建模，并提出“数字地质填图建模”概念。
数字地质填图三维模型可以对测区的地下地质情况进行一定深度的可视化，具有三维空间计算、动态更新、任意方向切制剖面图输出等功能，是区域地质调查成果的新型表达方式，也可作为今后深层次三维地质建模的基础。数字地质填图建模所需的基本数据源是填图过程中的野外路线数据和地形数据，比较容易获取，建模成本低，具有广阔的应用前景。
2.1.2 技术流程
数字地质填图建模可划分为3个关键步骤，即运用已有的PRB数据构建地质界面、将地质界面组合成面模型、将面模型生成体（网格）模型，其中以地质界面的构建最关键。按建面方法的不同，数字地质填图建模中的界面大致可以划分为模型边界面、DEM面、断层面、地层界面、第四系界面、岩体界面、残留顶盖界面、俘虏体界面等类型。地质界面的构建应遵循以下原则：①地质填图路线、产状数据为主要建模数据，其他数据为次要建模数据；②可信度高的数据约束、校正可信度较低的数据；③先建模型的边界面（DEM面、模型的底界面、模型的四周边界面），再建断层面，最后建其他地质界面。其他地质界面的构建，先建新地质体界面，后建老地质体界面；④由分段地质界面组合成更大范围的地质界面，界面构建应从小到大，并不断约束[64-65]。
地质界面的构建，利用PRB数据中的B数据及其对应的产状生成各个分段地质界面，再把同类分段地质界面组合生成更大范围的地质界面。对建立好的地质界面尽可能利用已有的钻孔等高精度数据进行约束、DSI等处理，生成更符合实际地质情况的地质界面（图 1）。
![]()
图 1 数字地质填图建模方法中的地质界面构建流程1—点间界线及编号；2—岩层产状；3—钻孔；4—地质界面；①—根据点间界线（B）及对应产状生成分段地质界面；②—根据分段地质界面生成地质界面；③—钻孔、点间界线、产状约束地质界面；④—裁剪地质界面Figure 1. Construction process of geological interface in Digital Geological Mapping Modeling2.2 地质剖面建模
2.2.1 地质剖面建模技巧
地质剖面数据是三维地质建模的重要数据源，运用地质剖面构建三维地质模型的方法应用较广泛[21-26, 47, 67-72]。地质剖面数据源主要是图切地质剖面、矿床勘探线剖面、中断平面图、物探解译的地质剖面等，建模流程主要由数据预处理及录入、地质界面构建、面模型构建和实体模型构建4个主要步骤组成（图 2）。
![]()
图 2 地质剖面建模方法流程图①—数据预处理及数据录入系统；②—地质界面的构建；③—面模型的构建；④—实体模型的构建Figure 2. Flow chart of Geological Profile Modeling层状地质体界面构建时，首先按照连接线生成面的方式构建一个初始地质界面。当地质剖面间距较大时，可以在相邻剖面线中根据地质体特征自动内插辅助地质界线。初始地质界面是较粗糙的界面，界面效果不佳且存在一些“坏”三角网、漏洞等错误，可以通过添加（删除、移动）三角网节点、增加（删除）三角网、打碎三角网、瓦解三角网、转换对角三角网等一系列简易操作解决。然后综合考虑地质体的属性与空间形态特征，以原始地质界线为约束条件，将初始地质界面数据在SWM（Structural Modeling Workflow）中重新生成地质界面，这个地质界面能够平滑过渡，效果更佳（图 3）。最后根据地表地质数据、剖面数据、钻孔数据等建模原始数据，采用约束-插值技术对构建的地质界面进行约束、插值、平滑处理，使生成的地质界面与所有原始建模数据吻合，且能平滑过渡。
![]()
图 3 剖面构建流程①—连接剖面线生成初始地质界面；②—以初始地质界面为数据，在SWM中生成地质界面Figure 3. Profile modeling process复杂岩体、褶皱构造、多值曲面等界面的构建，无法像构建层状地质体界面那样在SWM流程中自动、快速生成，而需要进行人机交互操作。如复杂岩体的构建，仅凭地质剖面中有限地质界线无法实现地质界面的准确构建，需要更多的人机交互过程：①创建一系列等间距的水平面，这些水平面需涵盖整个地质界线Z值范围；②提取不同高程上的地质界线与水平面的交点，分别生成不同高程的地质界线。根据少量交点数据生成一条穿过这些散点的包络线，加密包络线中的节点，并将包络线中交点所在位置的节点设置为控制点，通过DSI技术调整界线中的节点空间位置，使地质界线既能够平滑过渡，又能够穿过交点；③依次连接不同高程的地质界线生成地质界面；④运用剖面中的原始地质界线对地质界面进行约束、平滑处理，确保生成的地质界面与剖面地质界线完全吻合。
地质界面的构建过程中，空间数据的约束、插值是关键。当原始数据不够时，构建的地质界面会较粗糙。为了提高地质界面的可视化效果，需要在已有数据的基础上进行数据插值。GOCAD软件提供了DSI技术，可以在保持控制点空间位置不变的情况下，综合考虑地质体的属性与空间形态特征之间的关系对未知区域进行插值。通过不断调整插值点空间位置，使生成的面达到最优效果，最大限度地解决地质数据的不足。
建模过程中应用可信度较高的数据约束、校正可信度较低的数据。一般先根据主要建模数据构建初始地质界面，再用少量可信度更高的数据（如钻孔数据）约束和校正，达到构建的地质界面能够与所有地质数据尤其是精度较高的数据吻合。随着后期建模数据的不断增加，对已经构建的地质界面进行不断修改与更新。
2.2.2 地质图三维建模
应用地质剖面建模方法，根据平面地质图、地形数据可以直接构建三维地质概念模型，由地表模型和图切剖面模型两部分整合而成。
（1）建立地表模型。运用1：5万等高线数据构建DEM面。在DEM面的纹理属性中，设置地质图或遥感影像图，以便对其进行三维显示。
（2）制作系列图切剖面。根据DEM、地表地质图、钻孔、勘探线剖面等数据，制作一系列图切剖面。这些剖面在MapGIS软件中绘制而成，把绘制好的MapGIS软件格式的点、线、面数据转换成DXF格式数据，再将DXF格式的剖面数据导入建模软件中。运用剖面的4个角点坐标对剖面进行三维空间定位。
（3）地质概念模型构建。根据系列图切剖面，按由新到老的顺序依次构建地质界面。将DEM面、模型边界面和地质界面组合成为三维地质模型，即地质概念模型。
在建模之初，尚未有物探测深数据，笔者综合运用已有的平面地质图、地形图、钻孔、勘探线剖面等数据，制作系列图切剖面，建立了相山火山盆地地质概念模型。
2.3 多源数据融合建模
三维地质建模的原始数据来源具有多源性，涵盖地貌、地质、物探、化探、遥感、钻探等数据，数据类型和格式众多，难以有效利用。另一方面，建模数据获取难度大、成本高，又具有不均匀和不连续性，建模过程中总感觉数据不够。多源数据融合建模就是试图综合利用各类数据，更真实地反映地质体的三维空间特征。但建模数据之间的融合难度大，目前没有较完整的方法体系[47, 49]。笔者在相山盆地多源数据融合建模方面开展了一些探索，得出如下经验。
2.3.1 多源数据融合方法
（1）统一所有原始资料比例尺、坐标系统、投影参数，构建原始资料数据库。数据库采用平面直角坐标系，高斯克吕格（等角横切椭圆柱）投影，北京54坐标系或西安80坐标系，比例尺为1：1000，采用3°分带，坐标单位为m，Z轴正值方向向上。
（2）厘定建模单元，创建建模单元属性表，统一建模原始数据。浅表层三维地质模型的建模单元可以与地表填图单元相同，或稍作合并。深部三维地质建模可以依据主要建模数据源特点，在地表填图单元的基础上进行有效合并，重点突出目标地质体。对原始资料数据库中的点、线、面、体，按照建模单元分组归类。
（3）以DEM面为载体，实现地表所有建模数据的有效融合。将数字地质填图路线PRB数据投影到DEM面上，根据点间界线（B）及对应产状生成分段地质界面。地质点（P）、分段路线（R）、点间界线（B）、分段地质界面共同构成路线剖面[66]。用矢量地质图的地质界线切割DEM面，生成一系列小块DEM面，根据地质图分别对每个小块DEM面赋予不同颜色、属性、花纹等信息，生成一幅三维矢量地质图。而对于栅格地质图和遥感影像图，可以设置成为DEM面的纹理实现三维显示[66]。所有地表数据，以DEM为载体实现有机融合。
（4）根据地表填图路线和实测地质剖面构建浅表层三维地质模型，或根据地质图构建三维地质概念模型。对获取的深部物探数据进行反演与解译时，参考这些初始模型，以便更好地考虑地质体的三维空间展布特征。物探数据的地质解译和三维地质模型的修正反复交替进行，不断完善。
（5）根据覆盖全区的主要建模数据构建初始地质界面，再用高精度数据对已构建的初始地质界面进行约束，实现多种数据的有机融合。
2.3.2 约束-插值建模技术
约束-插值建模技术是以已知建模数据和地质属性为约束条件，对未知区进行插值处理，用以解决建模数据的不均匀性和不连续性。它可细分为点对线约束、点（线）对面约束、面对面约束3种。
（1）点对线约束。该方法类似于地质图连图，采用平滑的曲线将确定的空间点连接起来，确保了光滑曲线完全通过已知空间点，整条曲线又能平滑过渡。
（2）点（线）对面约束。当构建的初始地质界面与已有的高精度原始数据不吻合时，需要用后者约束地质界面。高精度原始数据可能是点、线、面等不同类型的数据，应将这些数据转换成离散点或线才能对面进行约束处理。
（3）面对面约束。面对面约束主要运用于地质界面组合生成面模型，主要功能是确保每个地质体能够完全密闭，并构建模型中不同面之间的空间拓扑关系。
3. 三维地质模型构建与分析
相山火山盆地位于江西省抚州市西南部，大地构造位置属于华夏板块湘桂赣地块的北东缘，钦杭结合带南侧。研究区面积约582km2，是一个火山-侵入杂岩区，主体岩石为早白垩世火山岩、火山-沉积岩系，下部打鼓顶组（K1d）岩性主要为流纹英安岩，上部鹅湖岭组（K1e）主要为碎斑熔岩，周边有浅成的早白垩世花岗斑岩（ηγπK1）呈环形岩墙状产出，西北角有上白垩统红层（K2）覆盖。盆地基底为青白口系（Qb）中浅变质岩。该盆地是中国目前最大的火山岩型铀矿田，素有“中国铀都”之美誉。从1956年至今，已有60余年的探、采历史，积累了丰富的勘探资料，目前找矿潜力仍然非常巨大。
在GOCAD软件平台上，依次构建了相山火山盆地三维地质概念模型、陀上幅三维地质模型、邹家山矿床三维模型、沙洲矿床三维模型、邹家山-居隆庵三维地质模型、相山火山盆地三维地质模型等不同范围、不同数据源的三维模型。
3.1 陀上幅三维地质模型
陀上幅三维地质模型采用数字地质填图建模方法构建而成，模型范围为1：5万陀上幅（G50E003008）标准图幅范围（东经115°45′~ 116°00′、北纬27°30′~27°40′），覆盖了相山火山盆地主体部分，面积约456km2，Z轴方向深度500m，底界面与DEM面平行（图 4）。建模单元与地表填图单元一致，共24个单元（表 1）。以野外数字地质填图路线PRB数据为主要建模数据源，地质图、地质剖面图、遥感影像数据等为辅助数据源，利用钻孔、勘探线剖面、中段平面图数据对地质界面进行约束（表 2）。
![]()
图 4 陀上幅三维地质模型a—三维遥感影像图（RGB（ETM7，ETM4，ETM1）+ALOS PAN）；b—矢量三维地质图；c—三维地质模型。1—第四系残坡积物；2—上白垩统莲荷组二段砾岩；3—上白垩统莲荷组一段砂岩；4—上白垩统塘边组三段粉砂岩；5—上白垩统塘边组二段砂岩；6—上白垩统塘边组一段粉砂岩；7—上白垩统河口组三段砾岩；8—上白垩统河口组二段复成分砾岩；9—上白垩统河口组一段砾岩；10—下白垩统鹅湖岭组二段中心相含花岗质团块碎斑熔岩；11—下白垩统鹅湖岭组二段过渡相碎斑熔岩；12—下白垩统鹅湖岭组二段边缘相含变质岩角砾碎斑熔岩；13—下白垩统鹅湖岭组一段砂岩、凝灰岩；14—下白垩统打鼓顶组二段流纹英安岩；15—下白垩统打鼓顶组一段砂岩、凝灰岩；16—青白口系上施组二段千枚岩；17—青白口系上施组一段片岩；18—青白口系库里组二段片岩；19—青白口系库里组一段片岩；20—早白垩世二长花岗斑岩；21—早泥盆世二长花岗岩；22—煌斑岩脉；23—硅化脉；24—地质界线；25—断层；26—河流Figure 4. 3D geological model of Tuoshang region表 1 各模型建模单元Table 1. The units of models模型 建模单元 模型 建模单元 模型 建模单元 陀上幅三维地质模型 Q K1e1 邹家山矿床三维模型 K1e2b 相山火山盆地三维地质模型 Q K2l2 K1d2 K1e2a K2 K2l1 K1d1 K1e1 K1e K2t3 Qbŝŝ2 K1d2 K1d K2t2 Qbŝŝ1 K1d1 T3z K3t1 Qbk2 Qb Qbŝŝ K3h3 Qbk1 邹家山-居隆庵三维地质模型 F Qbk K3h2 ηγπK1S Q Qbŝ K3h1 ηγD1J K1e ηγπK1S K1e2c q K1d ηγD1 K1e2b χ Qb F K1e2a F ηγπK1S 沙洲矿床三维模型 Q F ηγπK1S Qb F 注：Q—第四系；K2l2—莲荷组二段；K2l1—莲荷组一段；k2t3—塘边组三段；K2t2—塘边组二段；K2t1—塘边组一段；K2h3—河口组三段；K2h2—河口组二段；K2h1—河口组1段；K2—早白垩世红层；K1e2c—鹅湖岭组二段中心相；k1e2b—鹅湖岭组二段过渡相；K1e2a—鹅湖岭组二段边缘相；K1e1—鹅湖岭组一段；K1e—鹅湖岭组；K1d2—打鼓顶组二段；K1d1—打鼓顶组一段；K1d—打鼓顶组；T3z—紫家冲组；QBŝŝ2—上施组二段；QBŝŝ1—上施组一段；QBŝŝ—上施组；Qbk2—库里组二段；Qbk1—库里组一段；Qbk—库里组；QBŝ—神山组；Qb—青白口系变质岩；ηγπK1S—沙洲单元；ηγD1J—焦平单元；ηγD1—早泥盆世二长花岗岩；q—石英脉；x—煌斑岩脉；F—断层 表 2 数字地质填图建模数据Table 2. Data of digital geological mapping modeling数据类型 数据情况 标准图框 1:5万陀上幅标准图框 等高线数据 1:5万陀上幅标准图幅地形图 野外路线PRB数据 1:5万陀上幅数字地质调查野外路线232条 野外实测产状数据 主要是区调野外路线中的产状 钻孔数据 497个勘探钻孔 地质图 1:5万陀上幅地质图 遥感影像 TM/ETM+、ASTER、ALOS三种类型遥感影像 3.2 邹家山、沙洲矿床三维模型
从相山矿田的西部成矿区和北部成矿区各选一个典型铀矿床，构建三维地质模型，用于精细刻画不同地质体在三维空间的展布特征及其与铀矿体的空间关系。矿床三维模型的构建采用地质剖面建模方法，建模数据主要为勘探线剖面、钻孔编录资料、采矿中段平面图、等高线、地表地质图等（表 3）。矿床模型范围小、深度较大，与陀上幅浅表层、大范围的三维地质模型一起，为后期物探剖面地质解译、三维建模提供约束条件。
表 3 邹家山矿床建模数据源Table 3. The data sources of Zoujiashan deposit model数据类型 数量 数据源简介 勘探线剖面 26条 图片格式(JPEG)，勘探线间距50m 钻孔资料 415个钻孔 包含孔位表、测斜表、岩性表和矿化蚀变表等，孔间距25〜50m 采矿中段平面图 10幅 图片格式(JPEG)，中段之间的高差为40m 等高线数据 建模区域 矢量化数据 地质图 1幅 MapGIS格式的矢量化数据 邹家山矿床三维模型建模面积约1km2，建模深度与勘探剖面深度一致，最深处为1000m，建模单元共7个（表 1）。沙洲矿床三维模型建模面积约0.83km2，深度为360m（图 5）。用于建模的数据信息包括290个钻孔数据、5个采矿中段平面图和34条勘探线剖面数据。钻孔间距25~50m，采矿中段之间的高差为40m，勘探线间距50m。
![]()
图 5 沙洲矿床三维模型1—第四系；2—青白口系千枚岩、片岩；3—早白垩世二长花岗斑岩；4—断层Figure 5. 3D model of the Shazhou deposit3.3 邹家山-居隆庵三维地质模型
邹家山-居隆庵地区位于盆地西部，是重要的铀矿勘查区，区内有邹家山、居隆庵、牛头山等铀矿床，在牛头山地区深部存在丰富的铅、锌、银多金属矿化。对该成矿有利区域构建三维地质模型，可以揭示深部成矿地质要素的空间组合关系。在该区施测了14条相互平行的可控源大地电磁测深（CSAMT）剖面，剖面间距500m，单条剖面长4500m。以CSAMT解译的地质剖面为主要建模数据，地质图、钻孔、勘探线剖面图、中段平面图等数据为辅助数据。模型面积约32km2，深度3km，共有6个建模单元（表 1）。
3.4 相山火山盆地三维地质模型
运用多源数据融合建模方法对上述4个已构建的模型、深部探测等其他建模数据融合，构建了相山火山盆地三维地质模型。建模范围为东经115° 46′24″~116°03′30″、北纬27°27′1″~27°38′12″，涵盖了整个相山火山盆地，面积约582km2，深度达3km，有11个建模单元（图 6；表 1）。建模数据包括前期三维模型4个、均匀覆盖全区的大地电磁测深（MT）地质解译剖面19条、CSAMT地质解译剖面14条、地质图、勘探资料、地表等高线、遥感影像等。勘探资料有钻孔1459个、中段平面图和勘探线剖面图447幅。
![]()
图 6 山火山盆地三维地质模型1—第四系；2—上白垩统红层；3—下白垩统鹅湖岭组碎斑熔岩；4—鹅湖岭期火山通道；5—下白垩统打鼓顶组流纹英安岩；6—上三叠统紫家冲组砂岩；7—青白口系上施组千枚岩；8—青白口系库里组片岩；9—青白口系神山组千枚岩；10—早白垩世二长花岗斑岩；11—早泥盆世二长花岗岩；12—断层Figure 6. 3D geological model of Xiangshan volcanic basin3.5 模型分析
通过模型空间结构分析，对相山地区三维地质结构特征、构造演化及铀成矿有利因素取得了一些新认识：①盆地3层结构（基底岩系、火山岩系、红层盖层）特征清晰，青白口纪变质基底在三维空间上构成一个向北东东倾伏的复式背斜，相山火山盆地发育于核部倾伏端之上。盆地南西部有早泥盆世花岗岩侵入，由南西向北东倾伏，表明盆地西部具有变质岩-花岗岩双基底。②火山岩盆地下部打鼓顶组在西部厚度较大，呈东西向带状分布，与已发现的主要铀矿床在空间分布上关系密切。③鹅湖岭组在盆地中部厚度大，喷发中心的火山颈相位于相山主峰北西侧，总体上呈陡立管状，深部向南东倾伏，浅部向南东撒开，地表半径约2km。④早白垩世粗斑花岗斑岩呈向西开口的半环形分布。⑤格架性断裂构造主要为北东向，其次为北西向。新发现1条北东向断裂隐伏于白垩纪红层盆地之下，盆地北部发现1条弧形火山塌陷构造。2条北东向断裂夹持中部地堑块体，同时有东西走向的隆凹相间格局，使基底界面起伏较大。
相山铀矿田60余年来的地质勘查与开采，总结出“三盲”（盲断裂、盲岩体、盲矿床）、“三界面”（鹅湖岭组/打鼓顶组之间的组间界面、打鼓顶组/变质基底之间的基底界面、花岗斑岩/火山岩之间的侵入接触界面）找矿规律[73-75]。三维地质模型的建立，使这些找矿地质要素的三维空间特征变得清晰，已探明的铀矿化产于盆地内次级正负构造单元的转折部位，并与火山构造相联系。盆地西部变质岩-花岗岩双基底的鉴别，各地质体三维地质空间格局的系统解读，为相山地区新一轮铀矿勘查提供了重要依据。
4. 对三维地质建模的几点思考
4.1 完善地质思维的表达方式至关重要
三维地质建模关键在于地质表达方式的创新，需要地质学家从过去侧重地质思维（加上地史演化，属4D范畴）转变到兼顾重视现实表达方式上来。地质思维兼有逻辑思维与艺术思维，注重高度概括，忽略许多细节，在建模需要无缝连接时存在依据不足，同时用想象建立起来的各种模式难以用数学参数表达出来。编程员的信息技术思维难于理解地质体相互关系的深层次含义，导致在三维建模过程中，思维方式矛盾突显。
三维模型是在3D空间里对地质“事实”的记录、解释和简化，迫使研究者必须在3D空间里理解所有地质体、地质现象的几何形态和成因关系。因此，三维建模需要地质、物探、计算机人员通力合作，其中地质人员的全程参与特别重要，包括对地质现象的4D理解（含地史演化）、对地球物理场解译的主导作用。要重视培养掌握区域地质知识、物探解译常识和信息技术技能的复合型建模人员。
三维地质建模的技术方法仍欠成熟，缺乏研究者认可的统一流程。地质学家、地球物理学家、建模技术人员的水平，以及相互沟通能力是高质量建模的关键。不同类型地区三维地质建模的技术流程总结也势在必行。
4.2 实施分阶段建模，不断完善三维地质模型
由于当前的客观原因，三维地质调查无法以网格式部署物探、钻探工作量，所以目前的建模单位应该是目标地质体。矿集区三维地质建模要以控矿构造、含矿地质体为主要目标地质体。三维地质建模是一个不断完善的过程，当前可以依据现有资料构建地质概念模型，也就是地质图三维建模，然后有针对性地开展物探工作，不断完善深部地质结构解释。这样也便于实现多学科之间的数据融合，对已经建立的成矿模式在3D空间里接受检验。
4.3 数字地质填图建模值得探索并推广
数字地质填图建模值得进一步探索，可以作为地表区域地质填图的一种新型表达方式。过去区域地质调查图件的图切剖面、盆地模型、成矿模型等也试图揭示地下情况，但受条件限制只能用二维纸质图件表达。在1：5万图切剖面图上，1cm相当于500m深。在当今信息化时代，从地表填图数据直接进行三维建模（可暂定大致200m以浅），可以将这些图切剖面、立体模型图件实现实体化、可视化。同时，它也可以作为一种过渡性模型，用作钻孔设置的参考，或作为更深层次三维建模的约束条件。
4.4 三维地质模型的展示需要一种通用软件平台
三维地质建模专业软件价格昂贵，非专业人员使用也较困难。为了能够让更多人方便地使用构建好的模型，需要一种通用展示方式。笔者基于Adobe Acrobat软件平台开发了三维地质模型展示系统，实现了将Micromine、MapGIS、Surpac、GOCAD、Geosoft等众多专业地学软件构建的各类地质模型在通用平台上展示（转换为PDF文件），并拥有便利的人机交互功能。该展示系统主要有如下功能。
（1）多源地质模型导入。这是最基本的功能，通过U3D标准数据接口，可获取不同类型软件中的三维地质模型。目前，Adobe Acrobat提供2种获取三维模型的方式，一种是通过文件方式，直接读取含有模型的PDF文件，另一种是通过软件间的捕获功能实现。
（2）三维模型动态展示。可以利用Adobe Acrobat提供的快捷菜单，实现对模型的动态展示，如动态改变三维模型的位置、朝向，模型树展示、自定义剖面展示等功能。
（3）模型编辑与美化。利于Adobe Acrobat及Adobe 3D Reviewer工具相结合，可实现三维模型的编辑、材质设置、场景灯光设置等功能。
（4）多视图定义与转换。利于Adobe Acrobat视图管理功能，可实现对三维模型的任一视角、任一视图的操作与转换。
（5）灵活的人机交互。允许用户从表单域、书签、批处理序列、文档或页面动作调用JavaScript脚本，为PDF三维模型添加模型隐藏和显示、模型动画设置、模型渲染、模型缩放、剖面控制、模型指针设置等功能。
5. 结论
（1）三维地质建模应遵循由表及里、从粗到细的原则，分层次、分阶段建模，按照从地质概念模型到地质结构模型的顺序有效地开展工作。通过多种软件的遴选，笔者认为GOCAD软件比较适合于火山-侵入杂岩区的三维建模。
笔者尝试了在GOCAD软件平台上，利用数字地质填图数据直接进行浅表层三维地质建模，并提出“数字地质填图建模”概念。陀上幅三维地质模型由野外数字地质填图路线数据直接构建，钻孔数据作为约束条件。数字地质填图三维模型，可以对测区的地下地质情况进行一定深度的可视化，具有三维空间计算、动态更新、任意方向切制剖面图输出等功能。数字地质填图建模所需的基本源数据是填图过程中的野外路线数据和地形数据，比较容易获取，建模成本低，具有广阔的应用前景。
（2）相山火山盆地的其他4个三维地质模型的构建方法主要是地质剖面建模和多源数据融合建模，即以地质剖面数据为主要建模数据，钻孔、地质图等为约束条件。相山火山盆地三维地质模型和邹家山-居隆庵三维地质模型，主要根据大地电磁测深MT、CSAMT解译的地质剖面，结合地表填图数据进行建模，勘探线剖面图、采矿中段平面图、钻孔数据等作为约束条件；邹家山矿床三维模型和沙洲矿床三维模型，是根据勘探线剖面图、采矿中段平面图、钻孔数据、地质图等建立的。这些模型具有空间关系的确定性、三维可视化、可动态编辑与更新、分析和计算功能、成图与输出功能，使相山盆地深部信息可视化。
（3）三维地质调查工作强调多源数据的交互解译与有机融合，物探数据、勘探资料、地质图、地质演化规律相互印证。强调地质人员全程参与，与物探人员、计算机人员沟通交流，多学科融合，以便正确理解地球物理信息的地质含义，在3D空间里理解所有地质体、构造的几何形态和成因关系。
（4）相山地区三维地质模型的构建、地质体空间展布可视化，为该区地质结构及其与成矿关系的研究提供重要依据。三维模型揭示出火山盆地界面、火山口位置、火山通道形状、浅成斑岩体、基底褶皱和断裂的空间展布情况，使找矿地质要素的三维空间特征变得更清晰。相山盆地西部变质岩-花岗岩双基底的鉴别，也为本区铀成矿物源的研究提供了新思路。
致谢: 核工业北京地质研究院分析测试中心刘牧老师在全岩地球化学分析中提供帮助，中国地质大学（武汉）地质过程与矿产资源国家重点实验室在LA-MC-ICP-MS锆石Hf同位素实验测试中提供支持，中国地质大学（北京）朱弟成老师对本文提出指导性和建设性修改意见，在此一并表示感谢. -
![]()
图 1 大别造山带地质略图 (据参考文献[13]修改)
1—中-新生界 (K-E)；2—二郎坪群 (Pt3-Pz)；3—龟山岩组 (Pt2g)、南湾组 (Dn)；4—肖家庙岩组 (Z-Ox)；5—秦岭岩群 (Pt2qn)；6—桐柏-大别山变质杂岩 (Ar3-Pt1)；7—红安岩群 (Pt3)；8—榴辉岩；9—白垩纪火山岩 (K)；10—石炭系；11—燕山期花岗岩；12—晋宁期花岗岩；13—地质界线；14—断裂带及编号；15—大别造山带边界；16—钼矿床位置；Ⅰ—华北地块；Ⅱ—大别造山带；Ⅲ—扬子地块
Figure 1. Geological sketch map of the Dabie orogenic belt
![]()
图 2 大别山北麓姚冲钼矿区地质略图 (据参考文献[6]修改)
Q—第四系；PtDog—大别片麻杂岩；ArDb —太古宇大别岩群 (黑云斜长片麻岩)；K1zδο—早白垩世周湾独立单元（石英闪长岩）；Pz2γ—晚古生代含榴花岗岩；ec—榴闪岩及斜长角闪片麻岩；λπ—花岗岩脉；γπ—花岗斑岩脉；ζπ—正长斑岩脉；δμ—闪长玢岩脉；1—矿化体及编号；2—地质界线；3—不明性质断层；4—隐伏正断层/不明性质断层；5—钻孔及编号
Figure 2. Simplified geological map of the Yaochong Mo district
![]()
图 3 花岗岩岩相学特征
a—黑云母花岗岩标本；b—花岗斑岩标本；c—黑云母花岗岩显微结果；d—花岗斑岩的显微斑状结构.Qz—石英；Bt—黑云母；Pth—条纹长石；Pl—斜长石；Ms—白云母
Figure 3. Petrography of granite porphyry
![]()
图 4 姚冲花岗岩ANOR-Q′标准化矿物图解 (a, 据参考文献[19]) 和SiO2-A/CNK图解 (b)(Q′=100×Q/(Q+Or+Ab+An)；ANOR=100×An/(Or+An))
Figure 4. Q′-ANOR standardized mineral diagram (a) and SiO2-A/CNK diagram (b) of Yaochong granites
![]()
图 6 姚冲花岗岩原始地幔标准化微量元素蛛网图 (a) 和球粒陨石标准化稀土元素分布模式 (b)
（原始地幔标准值和球粒陨石标准化值据参考文献[20]）
Figure 6. PM-normalized incompatible element spider diagram (a) and chondrite-normalized REE patterns (b) of Yaochong granites
![]()
图 8 姚冲花岗岩（Zr+Nb+Ce+Y）-(Na2O+K2O)/CaO图解（底图据参考文献[32]）
A—A型花岗岩；FG—分异M+I+S型花岗岩；OGT—非分异M+I+S型花岗岩
Figure 8. Zr+Nb+Ce+Y versus (Na2O+K2O)/CaO diagram of Yaochong granites
![]()
图 9 姚冲花岗岩100*(TFeO+MgO+TiO2)/SiO2-(Al2O3+CaO)/(Na2O+K2O+TFeO) 图解（底图据参考文献[33]）
Figure 9. 100*(TFeO+MgO+TiO2)/SiO2 versus (Al2O3+CaO)/ (Na2O+K2O+TFeO) diagram of Yaochong granites
![]()
图 10 姚冲花岗岩造岩矿物结晶分异判别图解（底图据参考文献[35]）
PlAn50—斜长石 (An=50)；PlAn15—斜长石 (An=15)；Kfs—钾长石；Bt—黑云母；Ms—白云母；Grt—石榴子石；Amp—角闪石
Figure 10. Ba versus Sr and Rb versus Sr plots for Yaochong granites
![]()
![]()
表 1 姚冲矿区花岗岩主量、微量和稀土元素分析结果
Table 1 Major, trace and rare earth elements composition of Yaochong granites
样品号岩性 ZK0403-1 ZK0403-2 ZK0403-4 ZK0403-8 ZK0403-19 ZK0403-16 ZK0403-27 ZK0403-30 ZK0403-33 ZK0403-38 黑云母花岗岩 花岗斑岩 SiO2 76.64 77.24 71.28 72.41 72.61 76.15 75.35 75.01 75.49 74.70 TiO2 0.17 0.20 0.20 0.19 0.21 0.04 0.05 0.04 0.04 0.05 Al2O3 11.97 11.61 14.60 14.34 13.81 13.00 13.67 13.41 13.41 13.87 Fe2O3 1.09 1.17 1.22 1.08 1.25 0.31 0.42 0.51 0.57 0.54 FeO 0.39 0.48 0.43 0.41 0.49 0.46 0.30 0.46 0.22 0.33 MnO 0.03 0.04 0.03 0.03 0.05 0.00 0.02 0.01 0.03 0.03 MgO 0.23 0.32 0.27 0.38 0.45 0.25 0.10 0.16 0.07 0.08 CaO 0.64 0.73 1.15 0.99 1.07 0.55 0.71 0.77 0.77 0.58 Na2O 3.89 3.76 4.54 4.57 4.32 4.53 5.03 4.83 4.71 4.70 K2O 4.12 3.98 5.28 4.57 4.57 4.23 3.81 4.13 4.14 4.39 P2O5 0.04 0.05 0.07 0.06 0.05 0.02 0.01 0.02 0.01 0.03 LOI 0.71 0.62 1.08 0.94 1.10 0.48 0.53 0.68 0.49 0.78 Total 99.92 100.20 100.16 99.99 99.99 100.01 99.99 100.02 99.96 100.08 TFeO 1.37 1.53 1.38 1.53 1.61 0.74 0.68 0.92 0.74 0.82 A/CNK 1.00 0.98 1.00 0.96 0.99 1.00 1.00 0.97 0.98 1.03 ANK 1.10 1.11 1.15 1.11 1.15 1.08 1.10 1.08 1.10 1.11 Mg# 23.1 27.0 32.6 24.0 33.3 37.4 20.2 23.3 14.7 15.3 DI 94 94 92 93 92 95 96 95 95 95 Q 36.65 38.13 22.24 25.84 27.22 32.02 29.89 29.46 30.62 29.25 An 2.93 3.06 3.90 4.57 4.85 2.62 3.47 2.73 3.24 2.81 Ab 33.19 31.96 38.80 39.06 36.98 38.51 42.79 41.15 40.08 40.06 Or 24.55 23.63 31.51 27.28 27.32 25.11 22.64 24.57 24.60 26.14 C 0.05 - 0.22 - - 0.01 0.01 - - 0.37 Sc 5.35 4.97 5.62 5.47 6.46 3.39 6.09 6.23 3.99 2.74 V 21.3 16.5 20.9 14.7 18.7 10.2 7.89 12.7 5.69 5.29 Cr 7.01 3.45 5.23 1.07 2.21 2.15 0.26 2.19 0.34 0.05 Co 1.14 1.16 1.10 1.13 1.01 0.22 0.12 0.45 0.23 0.15 Ni 3.90 2.87 2.56 1.50 2.22 3.70 0.69 1.84 0.38 0.10 Ga 20.1 21.1 26.5 24.6 24 20.5 22.3 21.2 25.7 20 Rb 191 181 216 224 233 173 247 234 307 245 Sr 270 231 231 274 184 79.4 35.8 52.3 43.1 118 Y 4.56 5.44 5.86 5.83 5.76 4.54 5.89 6.18 3.9 2.16 Zr 116 134 143 136 175 67.6 72.1 70.9 60.6 44 Nb 17.5 19.7 18.2 15.3 17.1 21 21.5 23.6 20.8 10.9 Cs 2.45 2.13 2.7 2.5 2.72 1.63 4.2 10.9 2.71 2.41 Ba 642 818 1033 1053 840 107 76.5 63.7 65.9 82.1 La 19.6 28.1 30.9 32 33 6.03 6.94 6.67 5.17 3.43 Ce 35.9 48.6 55.8 54.8 56.6 6.03 6.94 6.67 5.17 3.43 Pr 3.57 5.26 5.72 6.03 5.93 1.02 1.21 1.16 1.19 0.73 Nd 12.9 18.6 20.4 21.6 20.3 2.99 3.53 3.53 4.2 2.56 Sm 1.87 2.81 2.99 3.30 2.74 0.47 0.48 0.53 0.76 0.43 Eu 0.47 0.54 0.59 0.68 0.63 0.09 0.09 0.10 0.10 0.12 Gd 1.43 1.66 1.99 2.14 1.70 0.34 0.58 0.65 0.48 0.47 Tb 0.18 0.28 0.24 0.30 0.20 0.05 0.10 0.08 0.10 0.07 Dy 0.71 1.05 1.10 1.14 0.74 0.45 0.65 0.58 0.41 0.15 Ho 0.20 0.16 0.19 0.21 0.16 0.11 0.12 0.15 0.09 0.06 Er 0.33 0.44 0.41 0.53 0.42 0.44 0.45 0.42 0.20 0.09 Tm 0.08 0.11 0.09 0.08 0.06 0.10 0.11 0.08 0.03 0.02 Yb 0.81 0.49 0.74 0.70 0.49 0.65 0.90 0.96 0.37 0.23 Lu 0.09 0.09 0.14 0.14 0.13 0.13 0.18 0.13 0.06 0.04 Hf 5.51 3.8 5.67 4.68 5.73 6.03 5.06 6.05 4.76 2.86 Ta 1.19 1.17 1.34 0.97 0.76 1.75 1.46 1.79 1.39 0.77 Pb 46.3 38.8 43.6 42.9 39 34.2 63.8 75.7 102 73.9 Th 15.2 18.3 16.3 14.9 15.6 26.8 28.5 24.8 21.9 8.46 U 7.26 9.33 7.5 7.02 4.55 22.4 24.7 24.9 28.9 12.5 Nb/Ta 14.71 16.84 13.58 15.85 22.56 12.00 14.73 13.18 14.96 14.23 Zr/Hf 21.05 35.26 25.22 29.06 30.54 11.21 14.25 11.72 12.73 15.38 ∑REE 78.13 108.20 121.30 123.65 123.10 18.91 22.28 21.71 18.34 11.82 LREE/HREE 19.43 24.24 23.77 22.62 30.55 7.31 6.22 6.11 9.51 9.55 (La/Yb)N 17.45 40.90 29.97 32.85 48.33 6.64 5.51 4.99 10.08 10.80 δEu=Eu/Eu* 0.84 0.71 0.70 0.74 0.83 0.66 0.53 0.52 0.49 0.84 TZr/℃ 804 816 820 811 838 757 762 759 747 723 注： N=球粒陨石标准化， 标准化值据参考文献[20]； A/CNK=molarAl/(Ca+Na+K)； A/NK=molarAl/( Na+K)； Mg#=100∗molarMg/(Mg+TFe2+)； δEu=Eu/Eu∗=2EuN/(SmN+GdN)； TZr=12900/[2.95+0.85M+ln(496000/Zrmelt)]， M=(Na+K+2Ca)/(Al∗Si)， Zrmelt为锆含量[21]； 主量元素含量单位为%， 微量和稀土元素含量为 10-6 
下载: 导出CSV
表 2 姚冲花岗岩锆石Hf同位素组成
Table 2 Zircon Hf isotope compositions of Yaochong granite
测试点 t/Ma 176Yb/177Hf 176Lu/177Hf 176Hf/177Hf 2σ ε (t) tDM1/Ma tDM2/Ma fLu/Hf ZK0403-8-1 137 0.020240 0.001064 0.282062 0.000020 -25.1 1677 2336 -0.97 ZK0403-8-4 138 0.013480 0.000517 0.282011 0.000010 -26.9 1724 2435 -0.98 ZK0403-8-4-1 138 0.014891 0.000564 0.281997 0.000010 -27.4 1744 2462 -0.98 ZK0403-8-5 136 0.019261 0.000951 0.282056 0.000016 -25.3 1680 2347 -0.97 ZK0403-8-8 142 0.016409 0.000697 0.282009 0.000012 -27.0 1735 2439 -0.98 ZK0403-8-12 141 0.024508 0.001021 0.281926 0.000144 -29.9 1864 2598 -0.97 ZK0403-8-13 137 0.016821 0.000753 0.282048 0.000016 -25.6 1683 2363 -0.98 ZK0403-8-15 138 0.014251 0.000667 0.281923 0.000062 -30.0 1852 2606 -0.98 ZK0403-8-19 136 0.007879 0.000383 0.282039 0.000015 -25.9 1680 2381 -0.99 ZK0403-8-20 142 0.018767 0.000897 0.282051 0.000010 -25.5 1685 2357 -0.97 ZK0403-8-22 137 0.014897 0.000691 0.282127 0.000026 -22.8 1572 2212 -0.98 ZK0403-8-23 137 0.018694 0.000855 0.282050 0.000014 -25.5 1686 2361 -0.97 ZK0403-8-24 141 0.022524 0.001021 0.282081 0.000018 -24.4 1649 2299 -0.97 ZK0403-16-1 140 0.015684 0.000730 0.282068 0.000014 -24.9 1655 2325 -0.98 ZK0403-16-4 141 0.018214 0.000887 0.282106 0.000023 -23.6 1610 2252 -0.97 ZK0403-16-7 139 0.023037 0.001115 0.282103 0.000052 -23.7 1623 2257 -0.97 ZK0403-16-8 137 0.018405 0.000888 0.282076 0.000013 -24.6 1651 2310 -0.97 ZK0403-16-10 138 0.025669 0.001237 0.282065 0.000015 -25.0 1681 2330 -0.96 ZK0403-16-11 142 0.018295 0.000806 0.282061 0.000017 -25.1 1667 2338 -0.98 ZK0403-16-12 137 0.020539 0.000951 0.282067 0.000015 -24.9 1666 2327 -0.97 ZK0403-16-13 137 0.018980 0.000999 0.282086 0.000014 -24.3 1642 2291 -0.97 ZK0403-16-14 140 0.014076 0.000607 0.282022 0.000011 -26.5 1713 2414 -0.98 ZK0403-16-17 137 0.020442 0.000972 0.282081 0.000021 -24.4 1647 2299 -0.97 ZK0403-16-20 137 0.014623 0.000627 0.282037 0.000014 -26.0 1694 2386 -0.98 ZK0403-16-21 135 0.018464 0.000887 0.282098 0.000013 -23.8 1620 2267 -0.97 ZK0403-16-22 138 0.020463 0.000977 0.282095 0.000017 -23.9 1628 2272 -0.97 ZK0403-16-23 141 0.016615 0.000836 0.282096 0.000015 -23.9 1621 2271 -0.97 ZK0403-16-25 142 0.015171 0.000650 0.282056 0.000016 -25.3 1668 2348 -0.98 注：εHf(t)=10000×﹛[(176Hf/177Hf)S-(176Lu/177Hf)S×(eλt-1)]/[ (176Hf/177Hf)CHUR.0-(176Lu/177Hf)CHUR×(eλt-1)]-1﹜；tDM1=1/λ×ln﹛1+[(176Hf/177Hf)S-(176Hf/177Hf)DM]/[ (176Lu/177Hf)S-(176Lu/177Hf)DM]﹜；tDM2=1/λ× ln﹛1+[(176Hf/177Hf)S.t-(176Hf/177Hf)DM.t]/[ (176Lu/177Hf)c-(176Lu/177Hf)DM]﹜+t；(176Hf/177Hf)S和 (176Lu/177Hf)S为样品测定值；(176Hf/177Hf)CHUR.0=0.282772，(176Lu/177Hf)CHUR=0.0332，(176Hf/177Hf)DM=0.28325, (176Lu/177Hf)DM=0.0384[24-25]；λ=1.867×10-11/a[26]；(176Lu/177Hf) c=0.015；t为锆石结晶时间
下载: 导出CSV
-
汤加富, 侯明金, 李怀坤, 等.扬子地块东北缘多期叠加变形及形成演化[J].大地构造与成矿学, 2003, 27(4): 313-326. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200304002.htm 王勇生, 朱光, 王道轩, 等.大别山东缘郯庐两期走滑剪切带形成的温压条件与造山带折返的关系[J].大地构造与成矿学, 2004, 28(3): 228-238. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200403001.htm 徐树桐, 江来利, 刘贻灿, 等.大别山区 (安徽部分) 的构造格局和演化过程[J].地质学报, 1992, 66(1): 1-14. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE199201000.htm 杨巍然, 杨坤光, 刘忠明, 等.桐柏-大别造山带加里东期构造热事件及其意义[J].地学前缘, 1999, 6(4): 247-253. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY199904011.htm 翟明国.华北克拉通破坏前的状态——对讨论华北克拉通破坏问题的一个建议[J].大地构造与成矿学, 2008, 32(4): 516-520. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200804017.htm Mao J W, Pirajno F, Xiang J F, et al. Mesozoic molybdenum deposits in the east Qinlin-Dabie orogenic belt: Characteristics and tectonic settings[J]. Ore Geology Reviews, 2011, 43: 264-293. doi: 10.1016/j.oregeorev.2011.07.009
毛景文, 华仁民, 李晓波.浅议大规模成矿作用与矿集区[J].矿床地质, 1999, 18(4): 291-299. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ199904000.htm 毛景文, 王志良.中国东部大规模成矿时限及其动力学背景的初步初探[J].矿床地质, 2000, 19(4): 289-296. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200004000.htm 李毅, 胡海珠, 陈丽娟, 等.大别山北麓姚冲钼矿床地质特征及找矿标志[J].地质与勘探, 2013, 49(2): 280-288. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201302013.htm 陶继华, 李武显, 李献华, 等.赣南龙源坝地区燕山期高分异花岗岩年代学、地球化学及锆石Hf-O同位素研究[J].中国科学 (D辑), 2013, 43(5): 760-778. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201305008.htm 高山, 张本仁, 金振民.秦岭-大别造山带下地壳拆沉作用[J].中国科学 (D辑), 1999, 29(6): 532-541. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK199906007.htm 王清晨, 丛柏林.大别山超高压变质带的大地构造框架[J].岩石学报, 1998, 14(4): 481-492. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB804.006.htm 杨泽强.河南商城县汤家坪钼矿辉钼矿铼锇同位素年龄及地质意义[J].矿床地质, 2007, 26(3): 289-295. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ200703004.htm 李俊平, 李永峰, 罗正传, 等.大别山北麓钼矿找矿重大进展及其矿床地质特征研究[J].大地构造与成矿学, 2011, 35(4): 576-586. http://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201104014.htm 刘清泉, 柳玉虎, 李永峰, 等.大别山北麓斑岩型钼矿床成矿地质条件及矿床成因[J].地质找矿论丛, 2013, 28(1): 27-33. http://www.cnki.com.cn/Article/CJFDTOTAL-DZZK201301004.htm Hu Z C, Liu Y S, Gao S, et al. Improved in situ Hf isotope ratio analysis of zircon using newly designed X skimmer cone and Jet sample cone in combination with the addition of nitrogen by laser ablation multiple collector ICP-MS[J]. Journal of Analytical Atomic Spectrometry, 2012, 27:1391-1399. doi: 10.1039/c2ja30078h
Geolaerts A, Mattielli N, Dejong J, et al. Hf and Lu isotopic reference values for zircon standard 91500 by MC-ICP-MS[J]. Chemical Geology, 2004, 206: 1-9. doi: 10.1016/j.chemgeo.2004.01.008
Liu Y S, Gao S, Hu Z C. Continental and oceanic crust recyclinginduced melt-peridotite interactions in the Trans-North China Orogen: U-Pb dating, Hf isotopes and trace elements in zicons from mantle xenoliths[J]. Journal of Petrology, 2010, 51(1/2): 537-571. https://www.researchgate.net/publication/268411794_Continental_and_Oceanic_Crust_Recycling-induced_MeltPeridotite_Interactions_in_the_Trans-North_China_Orogen_UPb_Dating_Hf_Isotopes_and_Trace_Elements_in_Zircons_from_Mantle_Xenoliths
Streckeisen A, Maitre R W. A chemical approximation to the modal QAPF classification of the igneous rocks[J]. Neues Jahrb Mineral Abh, 1979, 136: 169-206. http://www.oalib.com/references/9045201
Sun S S, Mcdonough W F. 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
Watson E B, Harrison T M. Zircon saturation revisited: Temperature and compositional effects in variety of crustal magma types[J]. Eeath and Planetaty Science Letters, 1983, 64(2): 295-304. doi: 10.1016/0012-821X(83)90211-X
张旗, 潘国强, 李承东, 等.花岗岩构造环境问题:关于花岗岩研究的思考之三[J].岩石学报, 2007, 23(11): 2683-2698. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200711003.htm Villaseca C, Orejana D, Paterson B A. Zr-LREE rich minerals in residual peraluminous granulites, another facter in the orgin of low Zr-LREE granitic melt[J]. Lithos, 2007, 96(3/4): 375-386. https://www.researchgate.net/profile/C_Villaseca/publication/229316595_Zr-LREE_rich_minerals_in_residual_peraluminous_granulites_another_factor_in_the_origin_of_low_Zr-LREE_granitic_melts/links/0912f509ccc82ae5c5000000/Zr-LREE-rich-minerals-in-residual-peraluminous-granulites-another-factor-in-the-origin-of-low-Zr-LREE-granitic-melts.pdf
Wu Y B, Zheng Y F, Zhang S B, et al. Zircon U-Pb ages and Hf isotope compositions of migmatie from the North Dabie terrane in China: Constrains on partial melting[J]. Journal of Metamorphic Geology, 2007, 25(9): 991-1009. doi: 10.1111/j.1525-1314.2007.00738.x
Rubatto D, Gebauer D. Use of cathodoluminescence for U-Pb zircon dating by IOM microprobe: Some examples from the western Alps[C]//Cathodoluminescence in Geoscience. Heidelberg, Germany: Springer-Verlag, 2000: 373-400.
Griffin W L, Pearson N J, Belousova E, et al. The Hf isotope composition of cratonic mantle: LA-MC-ICP-MS analysis of zircon megacrysts in kimberlites[J]. Geochimica et Cosmochimica Acta, 2000, 64(1): 133-147. doi: 10.1016/S0016-7037(99)00343-9
朱弟成, 莫宣学, 王立全, 等.西藏冈底斯东部察隅高分异I型花岗岩的成因:锆石U-Pb年代学、地球化学和Sr-Nd-Hf同位素约束[J].中国科学 (D辑), 2009, 39(7): 833-848. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200907001.htm 吴福元, 李献华, 杨进辉, 等.花岗岩成因研究的若干问题[J].岩石学报, 2007, 23(2): 185-220. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200706000.htm Chappell B W. Aluminium saturation in I-and S-type granites and the characterization of fractionated haplogranites[J]. Lithos, 46(3): 533-551. http://www.scirp.org/reference/ReferencesPapers.aspx?ReferenceID=1169996
Wu F Y, Wilder S A, Lo C H, et al. Highly fractionated I-type granites in NE China (Ⅰ): Geochronology and petrogenesis[J]. Lithos, 2003, 66(3/4): 241-273. http://www.academia.edu/13669525/Highly_fractionated_I-type_granites_in_NE_China_I_geochronology_and_petrogenesis
Li X H, Li Z X, Li W X, et al. U-Pb zicron, geochemical and Sr-Nd-Hf isotopic constrains on age and origin of Jurassic I-and Stype granites from central Guangdong, SE China: A major igneous event in response to foundering of a subducted flat-slab[J]. Lithos, 2007, 96(1/2): 186-204.
Whalen J B, Currie K L, Chappell B W, et al. A-type granites:Geochemical characteristics, discrimination and petrogenesis[J]. Contribution to Mineralogy and Petrology, 1987, 95: 407-419. doi: 10.1007/BF00402202
Sylvester P J. Post-collisional alkaline granites[J]. Geology, 1989, 97: 261-280. doi: 10.1086/629302
Xiong X L, Adam J, Green T H. Rutile stability and rutile/melt HFSE partitioning during partial of hydrous basalt: Implications for TTG genesis[J]. Chemical Geology, 2005, 218(3/4): 339-359. http://gemoc.mq.edu.au/Abstracts/Abs2006/XiongIAVCEI06.pdf
Janousek V Finger F, Roberts M, et al. Deciphering the petrogenesis of deeply buried granites: Whole-rock geochemical constraints on the origin of largely undepleted granulites from the Moldanubian Zone of the Bohemian Massif[J]. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 2004, 95: 141-159. doi: 10.1017/S0263593300000985
Arth J G. Behaviour of trace elements during magmatic processes: A summary of theoretical models and their applications[J]. Journal of Research of the Us Geological Survey, 1976, 4: 41-47. http://www.academia.edu/3420784/Behavior_of_trace_elements_during_magmatic_processes_a_summary_of_theoretical_models_and_their_applications
Mahood G, Hildreth W. Large partition coefficients for trace elements in high-silica rhyolithes[J]. Geochimica et Cosmochimica Acta, 1983, 47: 11-30. doi: 10.1016/0016-7037(83)90087-X
Green T H, Pearson N J. Rare earth element partitioning between titanite and coexisting silicate liquid at high pressure and temperature[J]. Chemical Geology, 1986, 74: 201-216. doi: 10.1007%2FBF00429424
Yurimoto H, Duck E F, Papike J J, et al. Are discontinuous chondrite-normalized REE pattern in pegmatitic granites systems the results of monazite fractionation?[J]. Geochimica et Cosmochimica Acta, 1990, 54: 2141-2145. doi: 10.1016/0016-7037(90)90277-R
Wedepohl K H. The composition of the continental crust[J]. Geochimica et Cosmochimica Acta, 1995, 59(7): 1217-1232. doi: 10.1016/0016-7037(95)00038-2
Blichert J, Albarede F. The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle-crust system[J]. Earth and Planetary Science Letters, 1997, 148: 243-258. doi: 10.1016/S0012-821X(97)00040-X
Patchett P J, Kouvo O, Hedge C E, et al. Evolution of continental crust and mantle heterogeneity: Evidence from Hf isotopes[J]. Contributions to Mineralogy Petrology, 1981, 78(3): 279-297. https://www.researchgate.net/publication/225258212_Evolution_of_continental_crust_and_mantle_heterogeneity_Evidence_from_Hf_isotopes
Nowell G M, Kempton P D, Noble S R, et al. High precision Hf isotope measurements of MORB and OIB by thermal ionization mass spectronmetry: Insights into the depleted mantle[J]. Chemical Geology, 1998, 149: 211-233. doi: 10.1016/S0009-2541(98)00036-9
Zhao Z F, Zheng Y F, Wei C S, et al., Zircon U-Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China[J]. Chemical Geology, 2008, 253: 222-242. doi: 10.1016/j.chemgeo.2008.05.011
Zhang S B, Zheng Y F, Wu Y B, et al. Zircon isotope evidence for≥3.5Ga continental crust in the Yangtze craton of China[J]. Precambrian Research, 2006, 146: 16-34. doi: 10.1016/j.precamres.2006.01.002
Zheng J P, Griffin W L, O'Reilly S Y, et al. Widespread Archean basement beneath the Yangtze craton[J]. Geology, 2006, 34: 417-420. doi: 10.1130/G22282.1
Xiong Q, Zheng J P, Yu C M, et al. Zircon U-Pb age and Hf isotope of Quanyishang A-type granite in Yichang: signification for the Yangtze continental cratonization in Paleoproterozoic[J]. Chinese Science Bulletin, 2008, 54: 436-446. http://www.oalib.com/paper/1673813
Zheng Y F, Fu B, Gong B, et al. Stable isotope geochemistry of ultrahigh pressure metamorphic rocks from the Dabie-Sulu orogen in China: Implications for geodynamics and fluid regime[J]. Earth-Science Reviews, 2003, 62(1/2): 105-161. https://www.researchgate.net/profile/Long_Li10/publication/223673324_Stable_isotope_geochemistry_of_ultrahigh_pressure_metamorphic_rocks_from_the_Dabie-Sulu_orogen_in_China_Implications_for_geodynamics_and_fluid_regime/links/0deec5254230e4dad5000000.pdf?origin=publication_list
Zheng Y F, Zhang S B. Formation and evolution of Precambrian continental crust in southern China[J]. Chinese Science Bulletin, 2007, 52(1): 1-10. doi: 10.1007/s11434-007-0015-5
Zhao Z F, Zheng Y F, Wei C S, et al. Post-collisional granitoids from the Dabie orogen in China: Zircon U-Pb age, element and O isotope evidence for recycling of subducted continental crust[J]. Lithos, 2007, 93(3/4): 248-272.
Zhao Z F, Zheng Y F. Remelting of subducted continental lithosphere: Petrogenesis of Mesozoic magmatic rocks in the Dabie-Sulu orogenic belt[J]. Science in China (Series D), 2009, 52(9): 1295-1318. doi: 10.1007/s11430-009-0134-8
Zhang S B, Zheng Y F, Zhao Z F, et al. Neoproterozoic anatexis of Archean lithosphere: Geochemical evidence from felsic to mafic instrusions at Xiaofeng in the Yangtze Gorge, South China[J]. Precambrian Research, 2008, 163: 210-238. doi: 10.1016/j.precamres.2007.12.003
Zhang S B, Zheng Y F, Zhao Z F, et al. Origin of TTG-like rocks from anatexis of ancient lower crust: Geochemical evidence from Neoproterozoic granitoids in South China[J]. Lithos, 2009, 113: 347-368. doi: 10.1016/j.lithos.2009.04.024
续海金, 叶凯, 马昌前.北大别早白垩世花岗岩类Sm-Nd和锆石Hf同位素及其构造意义[J].岩石学报, 2008, 24(1): 87-103. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200801008.htm 陈伟, 徐兆文, 李红超, 等.河南新县花岗岩岩基的岩石成因、来源及对西大别构造演化的启示[J].地质学报, 2013, 87(10): 1510-1524. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201310003.htm Zhao Z F, Zheng Y F, Wei C S, et al. Zircon U-Pb ages, Hf and O isotopes constrain the crustal architecrure of the ultrahigh-pressure Dabie orogen in China[J]. Chemical Geology, 2008, 253: 222-242. doi: 10.1016/j.chemgeo.2008.05.011
Liu F L, Xue H M. Review and prospect of SHRIMP U-Pb dating on zircons from Sulu-Dabie UHP metamorphic rocks[J]. Acta Petrologica Sinica, 2007, 23(11): 2737-2756. https://www.researchgate.net/publication/285632071_Review_and_prospect_of_SHRIMP_U-Pb_dating_on_zircons_from_Sulu-Dabie_UHP_metamorphic_rocks
Zheng Y F, Fu B, Gong B, et al. Stable isotope geochemistry of ultrahigh pressure metamorphic rocks from the Dabie-Sulu orogen in China: implications for geodynamics and fluid regime[J]. Earth-Science Reviews, 2003, 62: 105-161. doi: 10.1016/S0012-8252(02)00133-2
He Y S, Li S G, Hoefs J, et al. Post-collisional granitoids from the Dabie orogen: new evidence for partial melting of a thickened continental crust[J]. Geochimica et Cosmochimica Acta, 2011, 75: 3815-3838. doi: 10.1016/j.gca.2011.04.011
He, Y S, Li S G, Hoefs J, et al. Sr-Nd-Pb isotopic compositions of Early Cretaceous granitoids from the Dabie orogen: Constraints on the recycled lower continental crust[J]. Lithos, 2013, 156(2): 156-159.
Wang Q, Wyman D A, Xu J F, et al. Early Cretaceous adakitic granites in the Northern Dabie Complex, central China: Implications for partial melting and delamination of thickened lower crust[J]. Geochimica et Cosmochimica Acta, 2007, 71: 2609-2636. doi: 10.1016/j.gca.2007.03.008
Zhao Z F, Zheng Y F, Wei C S, et al. Origin of postcollisional magmatic rocks in the Dabie orogen: implications for crust-mantle interaction and crustal architecture[J]. Lithos, 2011, 126(1/2): 99-114. https://www.researchgate.net/profile/Yong-Fei_Zheng/publication/229421573_Origin_of_postcollisional_magmatic_rocks_in_the_Dabie_orogen_Implications_for_crust-mantle_interaction_and_crustal_architecture/links/0046352cbac0434632000000/Origin-of-postcollisional-magmatic-rocks-in-the-Dabie-orogen-Implications-for-crust-mantle-interaction-and-crustal-architecture.pdf
Ames L, Tilton G R, Zhou G. Timing of collision of the Sino-Korean and Yangtze cratons: U-Pb zircon dating of coesite-bearing eclogites[J]. Geology, 1993, 21: 339-342. doi: 10.1130/0091-7613(1993)021<0339:TOCOTS>2.3.CO;2
李锦轶.中朝地块与扬子地块碰撞的时限与方式——长江中下游地区震旦纪-侏罗纪沉积环境的演变[J].地质学报, 2001, 75(1): 25-34. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200101004.htm 张国伟, 张本仁, 袁学诚, 等.秦岭造山带与大陆动力学[M].北京:科学出版社, 2001: 1-855. Leech M L. Arrested orogenic development: Eclogitization, delamination, and tectonic collapse[J]. Earth and Planetary Science Letters, 2001, 185(1/2): 149-159. http://www.academia.edu/7686502/Arrested_orogenic_development_eclogitization_delamination_and_tectonic_collapse
Vanderhaeghe O, Teyssier C. Partial melting and flow of orogens[J]. Tectonophysics, 2001, 342(3/4): 451-472. https://www.researchgate.net/profile/Christian_Teyssier/publication/223623194_Partial_melting_and_flow_of_orogens/links/02e7e51624e80e25b4000000.pdf
Btyant D L, Ayers J C, Gao S. Geochemical, age, and isotopic constraints on the northern Dabie comples, east central China[J]. Geological Society of America Bulletin, 2004, 116(5/6): 698-717. https://www.researchgate.net/profile/John_Ayers/publication/230839344_Geochemical_Age_and_Isotopic_Constraints_on_the_Location_of_the_Sino-KoreanYangtze_Suture_and_Evolution_of_the_Northern_Dabie_Complex_East_Central_China/links/00b7d51840418141d7000000.pdf?origin=publication_list
Hacher B R, Ratachbacher L, Liu J G. Subduction, collision and exhumation in the ultrahigh-pressure Qin-ling-Dabie orogen[C]// Malpas J, Fletcher C, Ali J R. Aspects of the Tectonic Evolution of China. Geological Society, Special Publication, London, 2004, 226: 157-175.
Liu X C, Jahn B M, Liu D Y, et al. SHRIMP U-Pb dating on zircon of a metagabbro and eclogites from western Dabieshan (Hong' an Block), China, and its tectonic implications[J]. Tectonophysics, 2004, 394(3/4): 171-192. https://www.researchgate.net/publication/222696716_SHRIMP_U-Pb_zircon_dating_of_a_metagabbro_and_eclogites_from_western_Dabieshan_Hong%27an_Block_China_and_its_tectonic_implications
高阳, 叶会寿, 李永峰, 等.大别山千鹅冲钼矿区花岗岩的SHRIMP锆石U-Pb年龄、Hf同位素组成及微量元素特征[J].岩石学报, 2014, 30(1): 49-63. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201401004.htm 徐克勤, 胡受奚, 孙明志, 等.华南两个成因系列花岗岩及其成矿特征[J].矿床地质, 1982, 2: 1-14. http://www.cnki.com.cn/Article/CJFDTOTAL-KCDZ198202000.htm 卢欣祥.秦岭花岗岩揭示的秦岭构造演化过程[J].地球科学进展, 1998, 13(2): 213-214. http://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ802.018.htm 任纪舜.论中国大陆岩石圈构造的基本特征[J].中国区域地质, 1991, 2: 289-293. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD199104000.htm 李永峰. 豫西熊耳山地区中生代花岗岩类时空演化与钼 (金) 成矿作用[D]. 中国地质大学 (北京) 博士学位论文, 2005: 1-143. 河南省地质矿产局.河南省区域地质志[M].北京:地质出版社, 1989: 1-689. 毛景文, 张作衡, 余金杰, 等.华北及邻区中生代大规模成矿的地球动力学背景:从金属矿床年龄精测得到启示[J].中国科学 (D辑), 2003, 33(4): 289-299. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200304000.htm 刘清泉, 邵拥军, 张智慧, 等.大别山姚冲花岗岩锆石U-Pb年龄、Hf同位素及地质意义[J].中国有色金属学报, 2015, 25(2): 479-491. http://www.cnki.com.cn/Article/CJFDTOTAL-ZYXZ201502027.htm 刘清泉, 张智慧, 李永峰, 等.大别山北麓斑岩型钼矿床地质特征、成矿时代及其成矿构造背景[J].地质与勘探, 2014, 50(2): 199-215. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201402001.htm
计量
- 文章访问数: 2197
- HTML全文浏览量: 254
- PDF下载量: 374
