澳大利亚三维地质填图进展与实例

陈应军，严加永

(1.中铁二局第二工程有限公司，四川成都 610091;2.中国地质科学院矿产资源研究所，国土资源部成矿作用与资源评价重点实验室，北京 100037;3.东华理工大学，核工程与地球物理学院，江西南昌 330013)

1 引言

近年来，矿集区层次的综合地球物理探测和三维地质填图，在国际上已经逐渐成为解决成矿学基础地质问题和深部找矿的主要手段。三维地质填图可将不同学科、不同尺度的数据，进行基于三维空间的分析与对比，将矿产勘查由传统的二维定性解释向三维定性与定量相结合方向转变，深化理解地下地质体和构造之间的空间、成因和演化关系和成矿系统，结合流体模拟和对成矿过程的三维模拟和解释，预测深部找矿靶区，提高深部找矿的效率。

澳大利亚联邦科学与工业研究组织(简称CSIRO)开展了一项雄心勃勃的计划，称为“玻璃地球计划”，希望通过三维可视化和地质模拟等技术，使大陆表层一千米“像玻璃一样透明”，目的是使澳大利亚大陆地表以下一千米内所发生的深层过程变得透明，以便可以发现澳大利亚下一代巨型矿床，其实质就是进行深部三维填图。在该计划中，重磁电震等地球物理方法在确定深部控矿构造和围岩分层，建立初始模型，并通过三维重磁模拟和多尺度分析来约束和修改初始模型，确定地质界线、构造和岩体分布等信息的过程中起到了重要作用(Graham，2002;严加永等，2008)。根据这项计划，澳大利亚科学家已经进行了许多在大陆板块、陆块、成矿带、矿集区和矿田不同尺度的三维地质填图工作，并取得了良好的效果。

在我国，随着地质找矿工作的不断发展和深入，许多成矿带中易于开采的露天矿及埋藏较浅的矿床大部分都已查明，持续不断的开采造成该地区各类型矿山的保有储量急剧下降，一批大中型矿山陷入资源枯竭的危机之中，开展接替资源的找矿工作已经迫在眉睫。拓展深部找矿空间，挖掘深部的资源潜力已经成为实现我国资源可持续发展的战略选择。然而，深部找矿却面临巨大的挑战。从技术层面上说，深部找矿勘查要求所采用的探测技术探测深度大、精度高、抗干扰能力强、能适应复杂地形条件，目前任何一种单一的方法都难以满足进行500 m以下的矿产勘查的需求(吕庆田等，2005)。在这种情况下，吕庆田等(2013)在铜陵狮子山矿田进行了三维地质填图工作，成功预测了西湖、焦冲、严冲、五贵桥、天马山－章木山、观音冲－朱村－团山7个深部成矿靶区。祁光等(2013)结合泥河矿区的物性研究，也成功对泥河进行了三维地质填图工作。

本文以澳大利亚三维填图为例，从三维地质填图的概念出发，介绍了三维地质填图的工作流程和常用的工作平台。结合澳大利亚三个三维地质填图实例，介绍了三维地质填图在成矿带、油气和水资源领域中的应用和进展。通过对澳大利亚三维地质填图的综合分析，提出了对我国开展三维地质填图工作的启示和建议。

2 三维地质填图的概念

三维地质填图(Three－dimensional geological mapping)，又称为立体地质填图，它是指在进行地面地质填图的同时，利用各种物探(如重力、磁法、电法、地震勘探等)、钻探、遥感和地球化学等综合探测技术，按照一定的精度获取地表及地下相应深度的地质信息，并利用地球物理反演技术和三维可视化技术，揭示研究区地下一定深度(2～5km)的构造、岩体、地层、断裂的空间形态，制作区域三维空间地质图的方法。

2.1 三维地质填图要素

(1)确定研究区域;

(2)确定三维地质填图需要解决的问题;

(3)在(2)的基础上，确定填图精度;

(4)在满足精度要求的基础上，收集地质资料、物探数据，并根据需要补充进行物探、化探等综合探测;

(5)构建三维地质图。

2.2 三维地质填图的目的和意义

通过创建一个包含所有地质信息的虚拟三维空间，人们可以更加直观地理解地质体与构造之间的空间关系、构造成因和演化关系等，理解成矿系统，开展流体模拟，预测深部找矿靶区，并帮助研究人员理解区域3D结构框架，检验过去的地质模型或对区域地质的认识，可为地质找矿、地质灾害监测预警提供深部地质信息，保障国家地质生态安全、促进地球科学发展提供强有力的支撑。

3 三维地质填图工作流程及软件

3.1 三维地质填图流程

三维地质图的构建需要整合各种数据，包括地质资料(包括地质图、构造地质资料、地层层序资料等)、地球物理数据(地震资料、重磁数据、电磁探测数据等)和岩石物理数据(磁化率、密度)等，这就需要有一套合理高效的工作流程。虽然不同学者提出的流程的细节上各有不同，但基本上都包括以下三个方面:构建初始模型、二维或三维的重磁反演模拟、三维显示与地质解释(祁光，2013;吕庆田，2013)。流程图见图1。

图1 三维地质建模工作流程图(据吕庆田等，2013)Fig.1 Diagram showing workflow of 3D geological modeling(after Lü et al.，2013)

3.1.1 资料收集和处理

所收集资料包括简化地质图(其实是对地表岩性单元或地质单元在合适的范围内简化)、钻孔数据、地层年代学数据，还要进行岩石物性测量、岩性与物性对应关系分析(Kowalczyk et al.，2008)、重磁数据处理(如编辑、网格化、滤波和局部场分离等)和地震剖面解释等。在地质构造、岩性变化复杂的地区，对岩性单元的适当简化尤为重要，可以降低反演模拟的难度。在重磁反演中，深钻孔提供的深部地层单元的边界深度信息作为重要的约束，在反演中保持不变。区域场和局部异常分离出的局部异常将作为考量模型是否合理的依据。

3.1.2 二维地质模型的构建

确定的剖面间距，结合已有地质、钻孔资料的分析，依次推断、绘制建模区域的相应二维地质剖面，这些二维地质剖面应由若干紧密关联的模型体(也即地质体)构成，要求剖面大致反映对剖面区域的地层、构造、岩体和矿体空间分布的认识。

在这一步的工作中，地质学家对研究区域的地质理解应作为构建模型的重要参考。

3.1.3 2.5D/3D 反演

反演主要包括2.5D和3D重磁反演模拟。先前构建的二维地质模型作为2.5D重磁模拟的初始模型。使用人机交互“试错法”对2D剖面上的模型进行修改，直到获得合理的地质模型和满意的数据拟合为止。模型体的物性和空间形态应由物性数据和地质合理性决定。

完成研究区所有2.5D剖面的重磁模拟之后，按照剖面的空间顺序依次将2.5D模型拼合成3D模型进行计算，得出3D模型的理论异常，并与实际异常对比。拟合误差较大的地方，应返回到3D剖面进行修改，直到获得满意的结果为止。

在整个模拟过程中，物性与岩性的对应关系保持不变。

3.1.4 三维可视化与地质解释

将三维模型输出到3D可视化平台(如Encom PATM)。通过对三维模型的解释，可提取深部成矿信息，结合成矿模型开展深部成矿预测，还可作为储量计算、矿山设计和预测深部或边部矿体的参考。

3.2 三维地质填图常用软件组合

目前还没有任何一套软件平台可以同时完成所有的填图(建模)工作，只能使用多套软件进行填图工作。

2D剖面模拟工具(地质体在剖面上为多边形):GM －SYS Profile(NGA);ModelVision(Encom)。

3D多面体离散建模工具:ModelVision(Encom);QuickMag(Encom)

界面建模(深度、厚度与物性)工具:

3D物性反演工具:UBC MAG3D and GRAV3D(Mira Geoscience)

3D岩性反演工具:VPmg(Fullagar Geophysics);3D GeoModeller(Intrepid Geophysics)

GM－SYS软件是一种交互式、重磁协同地震的综合解释软件，软件包含的重磁数据2D剖面模拟功能，可模拟任意形状的二维体，并赋予其各类物性参数(密度、感磁、剩磁和地震波速等)，其中的3D处理模块，可进行重磁数据三维正反演模拟，可以建立由一个或多个三维层面的栅格图定义的地质构造模型(王林飞等，2011)。

ModelVision是一款专业的重磁处理解释软件，适用于找矿、油气勘探、工程和环境地球物理等行业。利用ModelVision可进行人机交互式2D剖面正演模拟，通过对已知的地质地球物理和钻孔资料的分析，依次推断、绘制区域的所有2D地质剖面，每条2D剖面可大致反映出对剖面穿过区域的地层、构造、岩体和矿体空间分布状况。ModelVision提供一个交互式界面，可进行重磁3D地质反演①。

QuickMag是一款先进和高效的三维磁力地质建模软件②。

UBC MAG3D and GRAV3D是一款领先的重磁反演软件，可进行重磁三维反演，在反演中，可加入先验信息进行约束反演，从而提高反演精度③。

VPmg是一款重力，TMI，重力梯度和磁场梯度的建模和反演软件，既可以进行地质约束(如钻孔、物性等约束条件)反演，也可进行无约束反演，VPmg 3D重力正向模拟和反演三维建模和标准重力数据反演。VPmg 3D用户可以选择反演任一重力(磁力、TMI)数据，或多个或单个重力(磁力、TMI)梯度张量的分量④。

GeoModeller是一款三维地质建模软件。包括地质编辑模块，钻孔及网格/网格管理，地球物理正演和反演模块等。GeoModeller可根据反演结果进一步细化模型⑤。

软件的选择需要综合考虑价格、功能、兼容性和软件的显示功能等因素。

4 澳大利亚三维地质填图实例

4.1 Eastern Yilgarn Craton成矿带3D填图

Eastern Yilgarn Craton(EYC)是一个包括了众多世界一流金矿和镍矿的重要成矿带。然而在EYC成矿带内部，各部分都存在各自不同成矿系统都有着不同的地球动力学过程，研究各个部分的地球动力学演化过程及成矿系统，意义十分重大。描绘EYC地区的构造和岩性结构的3D地质图是澳大利亚地球科学中心(GA)、西澳大利亚地质调查中心(GSWA)和西澳大利亚大学(UWA)的一项合作项目的主要成果之一。此3D地质图是地质和地球物理数据在三维上的处理和集成的成果的显示(Blewett et al.，2006)。

4.1.1 构建3D地质图

构建EYC地区的三维地质图，需遵循由2D到2.5D，再到3D的过程来整合地下空间形态的原则。本研究使用了来自从澳大利亚联邦政府、州政府，再到工业界和学术界人士的数据。值得注意的是，收集早期研究所采集的相关数据十分必要，因为在整合过程中再对后续重要数据进行重新评估虽然可以提高3D图的质量，但却会延迟构建3D图的工作进度。工作流程:

根据前期收集的跨学科的数据资料，按照工作流程(如图1)，进行构建三维地质图的工作。地质和地球物理相结合的方法在研究中已经被证明是非常适合三维空间的。地下地层、岩体单元和断层的及其几何形态的解析方法(如图2):

图2 ModelVision 二维正演(据Blewett et al.，2006)Fig.2 ModelVision 2D forward modeling(after Blewett et al.，2006)

确定地层和岩石的关系;

整合来自各种资料的信息，确定地层、断层要素等;

构建一系列二维剖面(截面);

使用位场正演验证二维剖面的合理性(使用Encom’s ModelVision软件)，根据地质资料调整二维剖面上的地层、断层等等要素，使二维剖面与实际最为合理。这一个过程可能需要反复进行才能达到地质资料、事实与地球物理数据符合最合理、最接近地下实际构造的情况。

将前面所得模型输入到三维软件中，如Gocad等等，并进行调整、完善。在构建三维地质图的过程中，二维剖面(截面)是约束三维曲面的基础。

4.1.2 地质与地球物理约束条件

包括断层、钻孔、地层接触关系、2.5维位场正演模型、三维反演模型，反射地震资料等。

4.1.3 所得三维地质模型及其意义

构建EYC地区的3D地质图(如图3)，突出了以下与找矿及科研相关的问题:

(1)模型构建前，对地质的良好认识是必不可少的;

(2)高品质高分辨率的地球物理数据集是至关重要的;

(3)为所有不同方法不同来源的地质与地球物理数据找到一致并且合理的地质解释;

(4)使用逆向建模的时候，需要一个良好的初始模型约束;

(5)在三维空间中对各种数据进行可视化、集成和分析十分重要;

(6)建模使用的2D和3D软件包的数据格式应能够互相兼容，这可以节省大量工作;

(7)在三维空间中进行研究，可使成矿系统和成矿过程的研究更为直观;

(8)三维建模技术可对成矿靶区优选提供显著支撑。

4.2 Capel and Faust盆地三维填图与油气远景预测

4.2.1 研究概况

澳大利亚地球科学中心(Geoscience Australia)研究人员在澳大利亚东部近海的Capel and Faust basins完成了一项通过的三维地质模型来预测油气的工作(Higgins et al，2011)。由于缺少先验信息和相关研究，对该地区前沿沉积盆地的解释和分析存在较高的不确定性。因此，澳大利亚地球科学中心收集了2006～2007年间该地区相关研究数据，包括岩石样本、高分辨率反射地震数据、位场数据和测深数据。

4.2.2 该区的三维建模与油气勘探

对这些2D数据集成的分析，难以获得令人满意的地质解释，原因在于:

图3 EYC三维模型(据Blewett et al，2006)Fig.3 Three-dimensional model of EYC(after Blewett et al.，2006)

地震测线间距较宽(15～40km);

缺乏盆地沉积物组成和年龄的直接证据(如深钻孔);

广泛存在的火山岩和侵入岩掩盖了地震剖面的特征;

高度区域化的盆地结构。

正是因为这些问题的存在，对这些数据进行有效的整合和分析十分必要。通过使用由GOCAD所提供的三维可视化和建模环境，各种有用数据被整合起来，通过处理和解释来创建一个集成的三维地质模型。这个模型可以使关键地质和油气远景问题得到回答。在盆地尺度上进行三维地质分析是澳大利亚近海研究的一个较新领域，通常，油气勘探公司通过进行油气成藏和流体建模来分析盆地勘探前景和规模。

4.2.3 构建集成盆地三维模型

Capel and Faust盆地的三维地质模型有两部分组成:水平面和断层面。基于二维地震解释，在水平面上定义了沉积和侵蚀的问题⑥，七个关键的水平面被用来约束三维地质模型，这些水平面约束了可能成为油气来源、储存和运移的主要沉积单位。断层面主要指示构造运动，也便于盆地可能产生的有机物的富集、运移和渗透。

地震数据是用于构建模型的主要数据参考，但地震测线相对较宽的间距要求整合其他数据类型。因此，重力、磁法以及水下地形数据都在建模过程中进行了整合。重磁梯度数据、边缘检测、向上延拓的处理结果，显示出了地下地形和区域结构。这种最大程度上的整合已有信息和资料的工作方法，使得三维地质模型得以构建。

4.2.4 油气远景评估

在项目过程中对三维层面和断层进行反复修改，最终得到的三维层位和断层整合较好，得到较为完美的三维地质模型(图4)。可以看到，该地区有一系列被断层切割的层面，断层面的几何形态也显示了断层随随即按的演变过程。

盆地内的三维地质和地质演化模拟表明，盆地中部和西北部最大沉积区的沉积物厚度足够生烃(Funnell et al.，2011;Hashimoto et al.，2011)。地震特征分析、区域构造重建和模拟表明，潜在的烃源岩在断裂前的沉积层、断裂时形成的沉积物和断裂形成的沉积物底层是连续的，储层在断裂上部连续，上断裂后段存在区域性盖层(Colwell et al.，2006;Hashimoto et al.，2011)。

图4 Capel and Faust basin盆地的3D结构和断层(红色)Fig.4 Three-dimensional basement structures and faults(red)of the Capel and Faust basins

三维建模也提高了研究区内可能的油气运移途径和富集构造的可视化程度。盆地模拟表明，最有生油潜力的烃源岩可能存在于该地区晚白垩纪，油气生成和运移大部分发生在白垩纪和早新生代。三维地质模拟还显示，一些背斜构造可能未被侵蚀，从而保留一些早期生成的有机物。

这些推断还需要进一步的详细调查。

4.3 Spring Hill地区三维填图与地下水资源评价

应用于水文地质学研究的三维(3D)地质填图法得到了越来越广泛的应用，越来越多的用户正在探讨开发如何将矿物质和烃类3D计算和可视化的方法应用于地下水。世界各地的许多团体已经从事地质信息三维可视化开发多年(Turner，1991;Berg et al.，2009)。

4.3.1 研究背景

在过去十年中，澳大利亚南部大部分地区的降雨量有所下降，造成许多地区用水紧张，澳大利亚维多利亚州政府制定了地下水管理计划。基于3D的水文地质方法能否产生更强大的、可重复利用的以及可以使地下水资源更好管理的水文地质结果，尤其是在用水紧张的时期更是如此，为更准确地确定总的水资源管理框架，迫切需要了解地表水与地下水之间的相互关系，并将3D地质和水文地质概念模型与数值模拟软件包相综合，提高预测能力，以应对未来气候变化的挑战，成为了一个研究热点。

4.3.2 研究区概况

Spring Hill地下水供水保护区(GSPA)面积为253km2，位于澳大利亚维多利亚州中部的Loddon河流域南端。该区由抬升和切割的古生代基岩地形构成，被古河谷填充，上覆第三纪沉积物，或被厚的晚第三纪的玄武岩流层序和火山喷发点中断。

研究区的地下水特性是不断变化，地下水位下降最大的区域与地下水开采强度最大的区域相吻合。一些地区的地下水位在冬季后就会恢复到抽水前的水位，而一些地区呈现长期下降趋势。

4.3.3 3D 地质模型的构建

本研究采用ArcGIS(ESRI)和GOCAD的套装软件，这是Paradigm地球物理勘探有限公司的软件产品。GOCAD是专门为石油工业而开发的软件，地质学家可以根据各种可用的数据集构建地下3D模型。

为构建Spring Hill区的3D地质模型，需要收集整理详细的地表地质GIS图、钻孔记录、地表地形的数字高程模型以及地球物理数据，也需要对钻孔记录进行准确而详细的核查，以进一步了解地层，然后在记录中标出具体的、重要的界线。在构建地下玄武岩结构的3D模型时只用到了220个钻孔数据。

同时，对异常表面进行校正，如果其他数据(如地表地质学或地球物理学)显示另一种解释可能性更大，则需要对其进行修正。一旦构建了地表地形、基岩古地形和任何中间层，就可以采用GOCAD建模工具来构建各个元素(地质单元或含水层)的体素模型(3D块)。之后在GOCAD中计算这些单元的体积和截面积，在此基础上计算地下水资源，也可用于实现可视化，例如，可以用不同的颜色表示地下水矿化度、地下水位或构造特征。在GOCAD中，可以很容易地对基岩古河谷内含水层的大小和形状进行可视化，结果表明，主要含水层单元是火山喷发点(或火山渣锥)和玄武岩。

4.3.4 主要研究成果

根据研究区的3D地质图(图5)，评价地下水资源，包括:

(1)根据估算的孔隙度，计算含水层的体积和该区最大可能的地下水量;

(2)根据露头面积、降雨量和水位图得出年补给量;

(3)通过玄武岩和深层砂矿含水层流出该区域的地下水流量;

图5 Spring Hill区的3D地质模型(据 Gill et al，2011)Fig.5 Three-dimensional geological models of the Spring Hill area(after Gill et al.，2011)

(4)到达主要排水线过程中地下水的损失，即基流。

本研究中开发的方法主要用于提高地下水资源评价能力，减少含水层系统的物理尺寸和结构形式的不确定性。此外，还可以提高将含水层系统进行概念化的能力，并确定其主要特征。三维地质软件，如GOCAD，能可靠且重复地计算含水层体积、截面积、含水层表面积、以及水位降低量。与传统的地下水资源评价方法相比，可以实现含水层系统的可视化以及建立基于物理的简单质量平衡模型，进而使分析工作更加透明和公开。

尽管数据并不丰富，但3D地质建模方法和质量平衡模型可以很好地利用现有的数据，采用明确和可重复的方式，对研究区的地下水资源进行量化。采用模型中的一些参数，可以验证重要假设是否合理，并保证结果可靠(Maimone，2004)。采用这一方法，生成了数字、3D地质和水文地质数据，可以存储在公共数据库内，为将来的用户和数据更新所用。

5 澳大利亚三维地质对我国的启示

随着经济的高速发展，对各种能源和矿产的需求飞速上升，环境也在急剧的恶化当中。开展深部探测，尤其是石油产区和大型矿集区的三维地质填图也成为了解决我国面临的重大资源和环境问题的迫切需求。除了推进深部找矿外，三维地质填图在重大地质灾害区进行地质灾害预测和水资源勘查方面，也将起到十分重要的作用，对可持续发展、维护国家安全和稳定，具有特别重要的科学意义和现实意义。三维地质填图是我国当今社会发展的需要，更是地质调查工作的一次重大创新。对我国三维地质填图工作提出几点建议:

(1)在三维地质填图试点工作的基础上，编制三维地质填图技术标准和规范，指导进一步开展可以满足不同需求、适合不同地质背景、不同比例尺的三维地质填图研究工作;

(2)加强遥感、深部地球物理探测、大陆科学钻探、地球物理反演及三维可视化技术等现代信息采集和处理技术的研究，在三维地质填图工作中，应用高精度、高分辨率的地球物理探测技术，提高探测精度，从而提高地质填图水平;

(3)加强三维可视化表达软件系统的开发和研究;

(4)不仅要在石油产区和大型矿集区开展三维地质填图，同时，应在我国主要经济区、重大工程建设区、重大地质灾害高发区进行三维地质填图，查明地下地质结构和工程地质条件，为经济区规划建设、地质灾害预防和治理提供基础资料。

(5)建立不同层次的深部综合地质信息管理与服务系统，加强不同层次的三维地质填图信息管理，提高成果的应用和社会化服务水平。

［注释］

① pbEncom.Modelvision［ED］.［2014 －4 － 30］.http://www.pbencom.com/products/Geophysics/encom －modelvision.shtml

② pbEncom.QuickMag［ED］.［2014 － 4 － 30］.http://www.pbencom.com/products/Geophysics/encom-quickmag.shtml

③ pbEncom.UBCInversion［ED］.［2014 －4 －30］.http://www.pbencom.com/products/Geophysics/UBC-inversion.shtml

④ Mira Geoscience.Fullagar Geophysics VPmg/VPem1D［ED］.［2014－4 － 30］.http://www.mirageoscience.com/our-products/software-product/fullagar-geophysics-vpmg-vpem1d

⑤ Intrepid-geophysics.GeoModeller［ED］.［2014 －4 －30］.http://www.intrepid-geophysics.com/ig/index.php?page=geomodeller

⑥ Colwell，J.B.，Hashimoto，T.，Rollet，N.，Higgins，K.，Bernardel，G.and McGiveron，S.，2010.Interpretation of Seismic Data，Capel and Faust Basins，Australia’s Remote Offshore Eastern Frontier.Geoscience Australia Record 2010/06，58pp.

Berg RC，Russel HAJ.2009.Thorleifson LH.Introduction:Three dimensional geologic mapping—An international perspective［C］.Three dimensional Geologic Mapping Workshop GSA Portland Oregon.2010/06，58

Blewett R S，Hitchman A P.2006.3D Geological Models of the Eastern Yilgarn Craton，Project Y2，Final Report［R］.Geoscience Australia Record 2006:272

Colwell J.，Foucher J.-P.，Logan G.，Balut，Y.2006.Partie 2，Programme AUSFAIR(Australia-Aairway basin bathymetry and sampling survey)Cruise Report［C］.344

Funnell R.，Stagpoole V.2011.Petroleum systems models for Capel and Faust basins，eastern offshore Australia［A］.GNS Science Report prepared for Geoscience Australia［C］.298

Gill B，Cherry D，Adelana M.2011.Using three-dimensional geological mapping methods to inform sustainable groundwater development in a volcanic landscape，Victoria，Australia［J］.Hydrogeology Journal，19(7):1349－1365

Graham Carr.2002.Glass earth-the geochemical front of penet rateovercast layer［J］.Land and Resources Inf ormat ion，.12:46 － 50

Hashimoto，T.，Higgins，K.L.，Rollet，N.，Stagpoole，V.，Petkovic，P.，Colwell，J.B.，Hackney，R.，Logan，G.A.，Funnell，R.，Bernardel，G.，2011.Geology and prospectivity of the Capel and Faust basins in the deepwater Tasman Sea region［A］.APPEA Journal and Conference Proceedings［C］，51，Extended Abstracts.

Hashimoto，T.，Rollet，N.，Earl，K.，Bernardel，G.，2008.Capel and Faust basins—new information from the offshore frontier between Australia，New Zealand and New Caledonia.［A］New Zealand Petroleum Conference Proceedings，Crown Minerals，Ministry of Economic Development［C］，Wellington.

Higgins K.L.，Hashimoto T.，Hackney R.，Petkovic P.，Milligan，P.，2011.3D Geological Modelling and Petroleum Prospectivity Assessment in Offshore Frontier Basins using GOCAD:Capel and Faust Basins，Lord Howe Rise［J］.Geoscience Australia Record，(01):98－107

Kowalczyk P.D.，Oldenburg N.Phillips，T.H.Nguyen，V.Thomson.2010.Acquisition and analysis of the 2007－2009 geoscience BC airborne data.In:Lane RJ(ed.).Expanded Abstracts of ASEG －PESA Airborne Gravity Workshop.ASEG －PESA Airborne Gravity Workshop，Canberra.Canberra:Geoscience Australia，115－124

Lü Qing-tian ，Shi Da-nian，Zhao Jin-hua，Yan Jia-yong，Xu Ming-cai.2005.Seismic method f or deep mineral exploration:problems and prospects-A case study of the Tongling ore district［J］.Geological Bulletin of China，24(3):211－218(in Chinese with English abstract)

Lü Qing-tian，Shi Da-nian，Tang Jing-tian，WU Ming-an，Chang Yin-fo，SinoProbe－03－CJ Team.2011.Probing on deep structure of middle and lower reaches of the yangtze metallogenic belt and typical ore concentration area:A review of annual progress of SinoProbe－03［J］.Acta Geoscientica Sinica，32(3):257 － 268(in Chinese with English abstract)

Lü Qing-tian，Chang Yin-fo，SinoProbe － 03 － CJ Team.2011.Crustal Structure and Three-dimensional Deep Exploration Technology for Mineral Resources:An Introduction to SinoProbe［J］.Acta Geoscientica Sinica，32(Supp.1):49 － 64(in Chinese with English abstract)

Lü，Qing-tian，Guang Qi，Jiayong Yan.2012."3D geologic model of Shizishan ore field constrained by gravity and magnetic interactive modeling:A case history［J］.Geophysics，78(1):B25 － B35

Maimone M.Defining and managing sustainable yield［J］.Ground Water 2004.42(6):809 －814

Qi Guang，Lü Qing-tian，Yan Jia-yong，Wu Min-an，Liu Yan.2012.Geologic constrained 3D Gravity and magnetic modeling of Nihe deposit—A case study［J］.Chinese Journal of Geophysics，55(12):4194－4206(in Chinese with English abstract)

Turner AK［Editor］.1991.Three-dimensional modelling with geoscientific information systems，NATO ASI series C:mathematical and physical sciences［M］，vol 354，Kluwer，Dordrecht，The Netherlands，443

Yan Jia-yong，Teng Ji-wen，Lü Qing-tian.2008.Geophysical exploration and application of deep metallic ore resources［J］.Progress in Geophysics，23(3):871 －891(in Chinese with English abstract)

［附中文参考文献］

吕庆田，史大年，赵金花，严加永，徐明才.2005.深部矿产勘查的地震学方法:问题与前景)铜陵矿集区的应用实例［J］，地质通报，24(3):211－218

吕庆田，史大年，汤井田，吴明安，常印佛，SinoProbe－03－CJ项目组.2011.长江中下游成矿带及典型矿集区深部结构探测——SinoProbe－03年度进展综述［J］.地球学报，03:257－268

吕庆田，常印佛，SinoProbe－3项目组.2011.地壳结构与深部矿产资源立体探测技术实验——SinoProbe－03项目介绍［J］:地球学报，32(Supp.1):49 －64

祁 光，吕庆田，严加永，吴明安，刘 彦.2012.先验地质信息约束下的三维重磁反演建模研究——以安徽泥河铁矿为例［J］.地球物理学报，12:4194－4206

王林飞，熊盛青，何 辉.2012.非地震地球物理软件发展现状与趋势［J］.物探与化探，35(6):837－844

严加永，滕吉文，吕庆田.2008.深部金属矿产资源地球物理勘查与应用［J］.地球物理学进展，23(3):871－891