滕吉文^1,2，薛国强^1,2,3，宋明春⁴

  1 中国科学院地质与地球物理研究所

  2 中国科学院地球科学研究院

     3 中国科学院大学地球与行星科学学院

  4 河北地质大学地球科学学院

     第一作者：滕吉文，研究员，博士生导师，主要研究方向：深部动力学

  通讯作者：薛国强，研究员，博士生导师，主要研究方向：瞬变电磁探测理论方法研究。

导读：

我国新一找矿突破战略行动正在加速推进，第二空间找矿正是主战场。21世纪以来，我国提出了向地球深部进军的口号，积极开展深部探测和地质找矿，开拓第二找矿空间，不断取得找矿新成果。

随着勘查工作的不断深入，第二找矿空间难度也越来越大，而地球物理电磁法勘探具有探测深度大、分辨率高的优势，特别是近年来，我国自主研发的电磁法勘探新理论、新方法和新仪器在第二找矿空间发挥了越来越重要的作用。

滕吉文院士长期倡导和并积极推动第二深度空间找矿，他在本文中不仅介绍了第二深度空间找矿理念及其成矿模型，评述了矿产资源电磁探测高分辨方法技术发展现状，还提供了2处电磁法找矿实例。本文研究成果对指导找矿具有重要参考。

------内容提纲------

摘要

0  引言

1  第二深度空间矿产资源探查科学理念

2  矿产资源电磁探测高分辨方法技术

3  深部矿产资源地球物理探测实例

3.1  河北围场银窝沟铜多金属矿SOTEM探测

3.2  胶东焦家金矿田CSAMT深部探测

4  结论

---------------

摘要：我国第一深度空间(0~500m)矿产资源量已不能满足国家战略需求，基于我国深部矿产资源形成的空间和条件的优越性及复杂性，提出了第二深度空间(500~2000m)探查大型与超大型矿产资源的科学理念。本文通过两个电磁探测实例，论证这一理念的正确性。首先，分析了近源电磁探测和基于飞行平台的电磁探测新方法，以及噪声去除、反演成像和综合解释等精细数据处理技术。利用近年来新发展的短偏移距瞬变电磁法，在河北围场铜多金属矿区进行了示范探测，圈定6处深部矿化异常，全部得到钻孔验证；采用可控源音频大地电磁法，在山东省胶东金矿区精细刻画出深部构造，并圈定了有利赋矿部位，揭示了焦家、马塘、寺庄等独立分布的金矿床，其主矿体均在第二深度找矿空间相互连接，实际属于同一金矿床，是世界罕见的超巨型金矿床。两个地球物理探测成果，是第二深度空间矿产资源成矿理论和找矿实践的有效范例。

关键词：第二深度空间；矿产资源；成矿理论；地球物理；找矿实践

0   引言

第一找矿、勘探与开发深度空间，是指由地表到深度为500m处；第二深度空间，是指500~2000m的深度范围；第三深度空间，系指2000m以下的找矿、勘探与开发空间。当前我国矿产资源勘查主要在浅部，公布的各类矿产资源量，大部分为浅部第一深度空间探查的矿产资源量。从理论上看，深部空间是地壳物质、能量交换和相互作用的集中区域，适合成矿元素的聚集和矿床的产生。在全球范围内，对于第二深度空间的矿产资源的探测和开发近年取得了明显的效果。我国第二深度空间的找矿潜力巨大，但需要借助高科技地球物理勘探来获取深部信息，以有效识别大型和超大型矿床。

地球物理方法基于物理场在地球介质中的响应特征，获取地下目标体的物性参数属性。在矿产资源探测中，当含矿岩体的物理参数属性与背景围岩存在差别时，即可能借助地球物理场实现对矿体的识别和定位，比如地震方法可以给出地下介质的弹性波速度属性，电磁法可获得电阻率和极化率属性，重力和磁法可以分别得到密度和磁化率属性。这些地球物理方法是寻找和圈定深部隐伏矿床的重要工具。随着探测目标体埋深的增大，地球物理方法和技术的革新成为推动第二深度空间矿产资源探测的重要力量。高分辨电磁探测技术、高精度重力及磁法、高精度地壳宽角反射／折射技术、高分辨率地震反射技术和小孔径精细三维地震层析成像技术等新技术的出现，为深部隐伏矿床和矿集区的探查提供了可靠的地球物理手段。

电磁法是矿产资源探测的核心方法之一。根据场源和装置形式的不同，电磁法可分为：大地电磁法(MT)、音频大地电磁法(AMT)。可控源音频大地电磁法(CSAMT)、广域电磁法(WFEM)、长偏移瞬变电磁法(LOTEM)(Strack，1992；严良俊等，2001)、时频电磁法(TFEM)和短偏移电磁法(SOTEM)等。上述方法中，除AMT和MT方法外均属人工源方法，其中CSAMT和WFEM为频率域方法，LOTEM和SOTEM为时间域方法，TFEM是在CSAMT和LOTEM基础上发展起来的频率域与时间域混合的双模方法。电磁法能够实现地下由浅及深的电阻率成像，对深部矿产资源圈定起到至关重要的作用。如何兼顾探测深度和精度，保障第二深度空间的探测效果，是当前电磁新方法和新技术研究的重要方向。

中国陆地由多个块体镶嵌而成，成矿条件良好，资源潜力巨大。世界三大成矿域在我国均有分布，且最新的成矿理论研究和深部预测结果均表明：我国陆地叠加有多期次、大规模成矿作用，在2000m深部蕴藏着巨量矿产资源。在第二深度空间矿产资源探查理念指导下，有效发挥地球物理探测技术，对推动深部矿产资源勘探具有重要的作用。本文以电磁法为例，介绍了近年新发展的、适于第二深度空间矿产资源探测的电磁法新技术，并结合典型矿区的电磁探测结果，给出了第二深度空间矿产资源探测的有效范例。

1  第二深度空间矿产资源探查科学理念

从地球动力学的角度来看，大型、超大型矿床的时空分布在本质上受到深部物质运移和深层过程的制约，而壳、幔物质的结构和构造则是控制地壳内部大规模流体的分异、调整、运移、聚集与形成的最主要的深部因素。深入探讨岩石圈及与相邻圈层间物质和能量交换及其深层动力过程，有益于认识成矿元素的离散与聚集的机制及时空格局，以及预测金属矿床出现的类型和聚集的空间部位，从而为大型、超大型矿床及矿集区的探查和预测提供导向。

矿产资源的形成与聚集，涉及地球深部物质与能量的交换和深层动力过程（图1）。因此，在第二深度空间金属矿产资源探查前，必须从机理上给出新的认识和建立新的运动学和动力学的模型。地面浅表处的金属矿产资源，包括大型、超大型矿床和多金属矿集区的形成，均是由地球内部物质与能量的交换所致（图1），涉及到地球深处壳、幔介质与构造格局、物质状态、运移行为、物质属性和其空间展布的深层动力过程。这表明：当今第一深度空间的已有矿产资源，并非是地下矿产资源的整体或全部，在第二深度空间必然存在着矿体的继续伸展和聚集。

图1  第二深度空间流体成矿模型

2  矿产资源电磁探测高分辨方法技术

随着地球物理方法和技术的革新，近年来多种地球物理方法已在第二深度空间探测中发挥重要作用，比如高精度重力和磁法、电磁法、高精度反射／折射地震法等。新一代的高精度磁力仪，如超导磁力仪、光泵磁力仪和旋进式磁力仪，也在近年来得到发展。分型、遗传算法、贝叶斯方法和神经网络算法等新方法应用于数据处理也取得了重要进展。结合地震高精度探测技术的发展和对成矿规律的认识，开展矿产资源勘查和成矿动力学研究是目前重要的研究方向。地震宽角发射／折射技术，能够获得深部的速度结构、泊松比和深部各向异性信息等，从而描绘出深部成矿带的空间构造和环境，对定位矿体或含矿构造具有重要作用。

相比于浅部矿产资源探测，深部探测难以采用浅表获取的数据对结果进行约束，因此具有更大的多解性和不确定性在解释地球物理反演结果时，浅表的地质与地球物理信息往往作为先验模型，以获取合理的球物理解释而在深部探测时，这种做法往往会导致数据严重偏离真实模型地球内部结构复杂，物理-化学环境状态异变，深部矿产资源勘探和开发的难度更大，此时，地球物理方法的探测精度成为关键获取由浅及深高精度探测结果，成为探测结果可靠性的有效保障，也是矿产资源探测获得成功的关键。近年来，在人工源电磁法基本理论的指导下，国内先后研发出兼具大深度和高精度的电磁探测方法，比如广域电磁法和电性源瞬变电磁法在这些方法上所涌现的数据处理和反演的新技术，进一步提升了方法的应用效果。

广域电磁法由何继善院士（2010）提出，它是一种人工源频率域电磁探测新技术。广域电磁法的观测方式类似于可控源音频大地电磁法，但观测区域较后者更大。通过采用更精确的电磁场响应定义“广义视电阻率”，广域电磁法能够获得比CSAMT更为宽广的观测区域、更高的探测效率和更大的勘探深度。目前，广域电磁法已在国内多个矿区开展了应用示范研究，包括隐伏金矿、多金属矿区等。高精度三维反演将是该方法的研究和应用重点。瞬变电磁法是利用不接地回线或接地导线向地下发送一次脉冲磁场，在一次脉冲磁场的间隙期间，观测二次涡流场，以进行地下矿产资源勘探的一种电磁勘探方法。电性源瞬变电磁法采用接地导线作为发射源，因其发射功率较大，能够开展大深度探测以满足第二空间矿产资源的勘探需求。根据观测区域与接地导线源之间的相互关系，可分为长偏移距瞬变电磁法，短偏移距瞬变电磁法和多道瞬变电磁法(MTEM)。LOTEM和MTEM主要应用于地壳构造探测、油气资源勘探和地热探查。SOTEM方法由国内学者薛国强等(2013)提出，在近源区域观测电磁场信号，偏移距通常为0.3~2倍探测深度。目前SOTEM已在包括煤矿、铁矿和银铅锌矿等矿产资源探测中获得成功应用。

电磁法高分辨探测技术研究主要还包括：

（1）基于飞行平台的探测技术。航空电磁法和半航空电磁法是近年得到快速发展的电磁法技术。通过将全部或者部分数据采集设备置于飞行平台上，这两类方法能够在复杂地形和地貌条件下开展探测，尤其在沙漠、沼泽和高原等难以以人力开展数据采集的复杂环境下具有明显优势此外，该方法具备探测效率高、采样密集和覆盖范围广等特点，在矿产资源的普查和详查方面均能发挥重要作用。经过几十年的发展和测试，目前时间域航空电磁法已是国内外的研发重点，发展了VTEM，SkyTEM，HeliTEM等商用探测系统，并已在矿产资源探测中发挥重要作用。近年来，国内也发展了包括CHTEM-I和CHTEM-II的航空瞬变电磁探测系统，中国地质调查局航遥中心和中国科学院地质与地球物理所也在相关科研项目的资助下，开展航空瞬变电磁探测系统的研发与测试工作。半航空电磁法仅需将接收装备搭载在飞行平台上，故可采用小型无人机开展相关研发和测试工作，目前国内包括成都理工大学、骄鹏公司和吉林大学等，都发展了半航空电磁探测系统，并开展了实际测试与探测。

（2）噪声去除技术。当今电磁法勘探大多要处于干扰环境中实施，即须在强干扰的背景下提取微弱有效信号，电磁法噪声去除技术成为探测取得成功的前提条件之一。近年来，基于原始数据的噪声去除技术得以发展，这其中既包括了在硬件系统设计时在噪声压制方面的考虑和设计，也包括对采集的原始时间序列的各种噪声压制技术。在硬件噪声压制方面，低噪声磁场传感器、等值反磁通传感器和磁通门的研发和设计能在一定程度上提高信噪比，同时也包括空心线圈传感器的优化设计等。在基于算法的噪声压制技术方面，近年来发展了基于数字递归陷波的噪声周期噪声压制、基于Hilbert-Huang的时频域噪声去除和基于正交多项式法的航空电磁运动噪声去除等新技术，基于深度学习和神经网络的相关技术开始被应用于瞬变电磁数据的噪声压制。尽管新的噪声压制技术不断涌现，但是目前还没有一种对于各类噪声均适用的去除技术在矿产资源探测的实际应用中，噪声源的定位和分析对于规避可能的噪声干扰至关重要。同时，记录并分析测区的纯噪声信息也能够为后续的数据处理提供必要的基础。

（3）反演成像技术。随着仪器装备的发展，近年来在航空和地空电磁方法的数值模拟、数据预处理方法和正反演方法的研究方面也取得了一系列研究成果。其中，作为从数据到电阻率映射的核心环节，反演和成像研究是近年来的研究热点，这其中包括视电阻率成像、横向约束反演、贝叶斯反演和2D/3D反演等技术。针对航空和半航空电磁探测系统的特殊性，对地面电磁法中得到有效应用的相关技术进行改造和拓展，发展了小波去噪技术、视电阻率成像和合成孔径成像等数据处理和成像方法。在反演方面，基于一维正演算子，横向约束反演在航空电磁法中已成为数据反演的主流方法。基于随机采样的贝叶斯方法近年来也得到有效应用，在恢复模型电阻率的基础上，更进一步的提供模型参数的不确定性，能够为数据的解释提供更坚实的支撑。在数据的高维反演方面，基于Footprint的基本概念所发展的高维反演方法，已完成模型测试和部分实测数据的反演。电磁数据的高维反演一般采用高斯-牛顿法(GN)、拟牛顿法(QN)、L-BFGS或非线性共轭梯度法(NLCG)。其中L-BFGS和NLCG每次迭代仅需开展一次正演模拟和一次伴随正演，无需计算海森矩阵，在电磁数据的高维反演中应用较多。但是，L-BFGS和NLCG方法的收敛速度较慢，当具有高效的数值模拟算子时，GN和QN方法是更好的选择。

（4）综合解释技术。反演能够从采集的电磁数据中获取电阻率信息，但是为了进一步推测地下地质结构或地质目标体，仍然需要结合地质信息和其他地球物理结果。数据融合和模式识别技术的发展，使得电磁数据的自动化或者解释成为一种新的趋势这种解释方法的速度更快，且更具客观性，能够对探测目标体进行快速、定量化预测。Campanyd等（2019）采用统计分析方法将电测井分类结果映射到电阻率成像剖面进行融合。Gulbrandsen等（2017）融合电磁法和测井数据，并将其直接用于地质建模。数据融合算法和不断增加的多源数据，使得相关研究逐步成为可能。此外，基于不同数据的误差范围，给出数据解释的可靠性评价，也对数据的综合解释具有重要意义。随着重点测区多源地球物理场、地质构造和物理参数测井数据的不断增加，使得综合解释成为一个重要的研究方向。

3  深部矿产资源地球物理探测实例

3.1 河北围场银窝沟铜多金属矿SOTEM探测

河北围场地处华北地台北缘的内蒙地轴中段，区内太古宙-古元古代发生了角闪岩相和麻粒岩相区域变质作用，形成了良好的金多金属矿源层中生代处于燕山陆内造山带北部，经历了印支运动、燕山运动和多期强烈的构造岩浆活动，具有非常有利的金、银、多金属成矿的地质条件。测区属燕山山系，海拔标高一般为1000~1300m，地形起伏严重，以往探测深度300m。

测区出露主要地层为侏罗系张家口组火山岩地层，局部为白垩系花吉营组张家口组分布在区内大部分地域，由于风化强烈，产状不甚清楚，大致沿火山口呈环状向外倾斜的低角度产出。地表局部覆盖浅层残坡积物和洪冲积物。测区位于乌龙沟-上黄旗“Y”型断裂东支西侧，受区域构造影响，区内主要发育两组断裂构造，主要是北东向，其次是北西向。两组断裂构造共同控制着区内六楞沟火山机构，环状及放射状次级断裂也较为发育，为矿液的运移、沉淀提供了有利的构造空间。

以往资料表明，不含金属矿物类岩石的电阻率均值多在2500Ωm以上，而含金属矿物类岩石的电阻率视矿物的含量而定，一般小于不含矿岩石，如含铜品位3％左右的浸染状铜矿石电阻率约为200Ωm。因此，矿区内矿化异常体表现为低电阻率特性。

为调查区内控矿、赋矿构造，圈定找矿靶区，在矿区范围内共布设10条SOTEM测线，每条测线长度为800m，测线间距为50m，测点间距20m，共有410个测点。发射源沿测区西南向村庄附近的一条道路布设，发射源方位大致为东偏南46°，长度为750m。测线L1距离发射源最近，线中心距发射源中心距离为645m，测线L10距离发射源最远，中心点距离为1090m（图2）。

图2  测线布置示意图

对10条测线的数据进行数据滤波与反演成像，得到每条线的电阻率-高程断面图。反演中，模型最大深度为1500m。图3给出了L4和L7线的反演结果。处理结果显示测区内电阻率由浅及深逐渐增高，由浅部的几十欧姆增大到深部的几千欧姆横向上，一定深度范围内的电阻率出现几处高低阻的突变。结合地质工作查明的断裂情况，推断电阻率-高程断面图中出现的横向电阻率突变位置对应于断裂带穿过的地方。因此，在图中可以定性地圈定出几条断裂，如图3中红线所示。

将所有测线处理结果汇总成三维立体图，如图4所示。根据成矿理论及矿床富集规律，特别是针对该区的岩浆成矿机制，断裂带附近尤其是多个断裂交汇处是成矿元素容易随岩浆运移、汇集的区域，是找矿的主要目标区。图5所示的电阻率平面图中，均可以圈定几处明显的断裂交汇导致的低阻区域，结合地表踏勘及已有地质资料，认为该区域可能是矿体的富集区。

图3  L4、L7线反演电阻率-高程断面图

图4 三维立体图与平面等值线图

图5 电阻率平面分布图与建议钻孔位置

根据SOTEM电阻率-高程断面图及不同深度的电阻率平面图，我们推断出6处为可能的赋矿位置，如图5所示，建议实施钻探验证，建议钻探深度600m，并有可能对钻孔1和钻孔6钻进至1000m。钻孔结果证实深部发现了铅、锌、铜、铁等矿化岩石（表1）。

表1  钻孔揭露含矿岩石信息

3.2  胶东焦家金矿田CSAMT深部探测

胶东位于华北克拉通东南边缘，是中国最重要的金成矿区，已探明金资源量超过5000t，占全国的近1/3，其中，近年探明的500~2000m第二深度空间金资源量已达第一深度空间的2倍。地球物理方法在第二深度空间找矿中发挥了重要作用。

焦家金矿田处于胶东的西北部，其西侧邻近渤海湾。该区主要由早前寒武纪变质岩系和侏罗—白垩纪花岗岩类组成，地表被第四系松散沉积物浅覆盖。焦家断裂是焦家金矿田的控矿断裂，该断裂长约60km，宽50~500m。断裂总体走向NNE，但变化较大，在10°~75°之间呈S型变化；断裂倾角10°~78°，呈浅部缓、深部陡，陡、缓交替的铲式阶梯状特征。断裂下盘发育较多与走向平行或呈“入”字型相交的次级断裂。断裂总体发育于早前寒武纪变质岩系与侏罗—白垩纪花岗岩类接触部位。该断裂控制了焦家、新城、河东、河西等20余处金矿床，已探明的第一深度空间金资源量约500t，第二深度空间金资源量超过1000t。

研究表明，焦家断裂为浅部倾角陡向深部变缓的铲式断裂，而且由浅至深显示陡、缓交替变化的台阶式或坡坪式特点，矿体厚大部位赋存于台阶的陡缓转折和平缓部分，在2000m垂向深度以浅出现浅部和深部两个大规模的赋矿台阶（或称矿化富集带，其中各有多个小规模的赋矿台阶），二者之间为无矿间隔或弱矿化带。由此建立了胶东深部金矿阶梯成矿模式，金矿控矿断裂沿倾向出现若干个倾角由陡变缓的变化台阶，破碎带蚀变岩型矿体主要沿台阶的陡、缓转折和平缓部位富集，构成阶梯式分布型式（图6）。根据阶梯成矿模式和与成矿有关地质体的物性特征，首先采用频率域电磁探测深部构造和有利赋矿部位，识别控矿断裂向深部倾角变缓的台阶，即在地表通过高精度地球物理探测，查明控矿断裂或成矿地质体向深部的结构变化，根据阶梯模式推断深部矿的位置、规模，这一找矿方法被称为深部金矿阶梯找矿方法。

图6  胶东深部金矿阶梯成矿模式

在焦家金矿田的寺庄矿区360勘探线施工的CSAMT剖面，揭示了焦家断裂的深部特征和深部矿体位置。剖面长度3.2km，探测深度2.0km。探测结果表明，视电阻率曲线可分为中上部的低阻带、底部的高阻带和之间的过渡梯级带（图7）。根据该区主要地质体的物性特征推断，剖面中上部的低电阻率电性层与新太古代变质岩系相对应，底部的高阻电性层与侏罗纪玲珑花岗岩相对应。浅部矿体已得到钻孔控制，过渡梯级带等值线自上而下呈舒缓波状定向延深，与焦家断裂带相对应。在2200号测点以东的已知钻孔控制区域，2300—2750号点之间见有多层金矿体，该段的视电阻率等值线显示间距宽大、同步向下弯曲的特征，为断裂倾角较缓段，符合阶梯成矿特征。在2200—0号测点之间的未知地段，视电阻率等值线显示有4段向下弯曲变化相对较大、等值线间距相对较大的区段，分别位于1900—2200号点间、1200—1500号点间、700—1000号点间和0—300号点间，推断为断裂倾角变化的台阶位置，即为深部成矿的有利部位。

图7 寺庄矿区360线CSAMT法剖面和推断的深部地质成矿特征

1-已控制断裂蚀变带；2-推测断裂蚀变带；3-已控制金矿体；4-推断深部矿体；5-已施工钻孔位置

根据阶梯成矿模式的电磁探测结果，经过系统的深部钻探，在940~2030m深度范围内探获金金属量389t，探获矿体的分布范围和位置与本文上述的描述大致一致，验证了阶梯找矿方法在胶东深部找矿中的有效性。通过深部金矿勘探还发现，以往认为独立分布的焦家、马塘、寺庄、寺庄深部、马塘深部、焦家深部、纱岭-前陈等金矿床，其主矿体在第二深度找矿空间相互连接，实际属于同一金矿床，这一矿床已累计探明金资源量超过1200t，为一世界罕见的超巨型金矿床。

4  结论

我国处于快速工业化和经济高速发展时期，亟需巨量矿产资源作为支撑，矿产资源的紧缺已成为现实，对外依存度不断增高，实施第二深度空间金属矿产资源的找矿、勘探和开发已迫不及待。在深部金属矿产资源的探查和开发中，必须依靠新的认识和新的理念，采用先进的、高精度的地球物理探测方法和技术，并建立起有自主特色的深部矿产资源探查理论、方法和技术体系。本文所举的两个实例证实了在第二深度空间探查与发现大型和超大型矿床有着极大的远景，有力地支撑了第二深度空间矿产资源成矿理论的发展，今后的地球物理找矿研究还应该关注以下几个方向：

（1）加大第三深度空间地球物理精细探查工作。开展第二深度空间探测时，厘清深部矿产资源的形成、演化和分布规律，开展成矿作用和地壳深部结构以及深层动力学研究，是确保探测结果可靠性和延续性的重要手段。为此，开展第二深度空间矿产资源探测，不应将目光仅局限在500-2000m的第二深度赋矿空间，而应探索更深部的地球物理参数分布，解译深部的成矿和构造环境，揭示矿产资源的垂向分布和演化。因此，地球物理方法和技术的创新，应不断向更大深度的探测空间延伸，着力发展具有更大探测深度潜力的地球物理方法。

（2）研发立体化探测新技术体系。地球物理数据采集通常仅能在地表有限区域开展，随着近年钻探能力的提升和航空飞行器等装备的飞速发展，在地下深处和空中也能观测到地球物理数据。这些数据与地面观测的数据能够形成交叉互补，会极大的减弱观测数据反演结果的多解性。尽可能利用潜在的观测空间，开展井中-地面-空中立体探测，成为获取深部可靠信息的重要方向。当前立体化探测仍然出于分散研发和粗放应用阶段，如何发展一体化探测方法和数据处理技术，是开展大规模、集成化深部探测的重要途经。

（3）综合利用多源地学信息。矿产资源的形成与地球内部物质和能量交换，及其深部过程和动力学机制密切相关。加强第二深度空间的找矿工作，离不开地球物理学与相关学科的交叉和互促。尤其将地球科学作为整体开展研究时，如何有效整合多学科数据和信息，开展基础和理论性研究，从多方位、多角度解析多源信息，地球科学从定性分析向定量刻画挺进凝练交叉学科的核心科学问题，推动地球科学学科交叉融合，将成为新时代找矿勘查的重要方向。

致谢 感谢中国科学院地质与地球物理研究所李海和陈卫营副研究员在本文成文过程中提供的帮助。

  原文来源：滕吉文，薛国强，宋明春.2022.第二深度空间矿产资源探查理念与电磁法找矿实践.地球物理学报，65(10)：975-3985,doi：10. 6038/cjg2022P0762.