关于在傅里叶域使用非均匀快速傅里叶变换进行二维磁异常高效建模研究的学术报告
本研究题为 “Efficient 2D modeling of magnetic anomalies using NUFFT in the Fourier domain”，由来自桂林理工大学电子工程与自动化学院的Dongdong Zhao（通讯作者）、Central South University地球科学与信息物理学院的Xulong Wang与Jianxin Liu，以及桂林理工大学地球科学学院的Qianjiang Zhang共同完成。该研究论文发表于Springer Nature旗下的学术期刊 *Pure and Applied Geophysics*。
一、 研究的学术背景 该研究隶属于应用地球物理学领域，特别是针对磁法勘探这一成熟的地球物理勘探方法。磁法勘探已广泛应用于矿产勘查、区域地质填图、管线探测、考古测绘等多个方面。二维正演建模在磁法勘探的直接解释中扮演着核心角色，其目标是根据已知的地下地质体形态和磁化率分布，计算出其在地表产生的磁异常。传统上，有多种方法（如解析法、积分方程法）可用于计算简单规则地质体引起的磁异常。然而，对于具有任意形状和非均匀磁化率分布的复杂地质体模型，通常难以获得解析表达式，因此需要高效的数值模拟方法。
傅里叶域方法是解决大规模数值模拟问题的一种高效途径。传统方法通常采用标准快速傅里叶变换（Fast Fourier Transform, FFT）进行计算。但FFT在应用时存在局限性：它要求空间域和傅里叶域数据必须在均匀分布的规则网格上进行采样，导致在实际建模时为了减小边界截断效应（即边缘效应）和强加的周期性影响，往往需要对模型网格进行大幅扩展。这虽然能提升精度，但极大地增加了计算负担和内存需求。虽然存在如Gauss-FFT等方法试图在精度和效率间取得平衡，但仍存在计算时间较长或难以兼顾的问题。非均匀快速傅里叶变换（Nonuniform Fast Fourier Transform, NUFFT）作为一种能够处理非均匀采样数据的算法，为改善上述问题提供了可能。
因此，本研究的主要目标是：提出并验证一种在傅里叶域中利用NUFFT技术进行二维磁异常正演建模的高效、高精度方法。研究旨在解决标准FFT方法因强制周期性和均匀网格要求带来的精度限制，并通过引入一种针对起伏地形的快速算法，进一步提升复杂地形条件下的计算效率。
二、 研究的详细工作流程 本研究主要包括理论推导、算法设计、数值实验验证三个主要流程，涉及模型构建、算法实现和性能测试等多个步骤。
1. 理论推导与算法框架建立 首先，研究从空间域的磁位积分表达式出发。对于一个任意形状、具有任意磁化率分布的二维地质体，将其剖分为多个具有矩形截面的单元，并假设每个单元内磁化率为常数。通过在水平方向（x方向）对积分表达式进行一维傅里叶变换，将空间域的二维积分问题转换为傅里叶域中一系列关于深度（z方向）的一维积分问题。这一转换使得建模网格可以在垂直方向进行非均匀设置（例如，浅部加密，深部稀疏），从而更灵活地适应地质模型。对于每个单元的一维垂直积分，研究推导出了其解析解，进一步保证了计算精度。
在这一理论框架下，傅里叶域的磁位（以及通过微分关系导出的磁力场和梯度张量）表达式最终归结为两项关键计算：一是磁化强度的傅里叶变换，二是与深度相关的核函数的求和。传统的标准FFT方法要求均匀网格采样。本研究则引入1D类型3的NUFFT算法来计算磁化强度的傅里叶变换。NUFFT算法的核心思想是结合局部插值和标准FFT，通过高斯卷积与解卷积操作，允许在空间域或傅里叶域进行非均匀采样，从而克服传统FFT的网格限制。
2. 关键技术设计与优化 波数选择方案： 由于傅里叶域磁位谱随波数增大呈指数衰减（并非严格单调），为了在保证精度的同时减少计算量，研究提出了一种在对数轴上等间隔选取波数的优化方案。即在高波数区（谱值较小）采样较稀疏，在低波数区（谱值较大）采样较密集。通过设定一个最小波数阈值（例如10⁻⁴）和奈奎斯特频率，在对数坐标上将波数范围划分为若干区间进行采样。这种策略有效平衡了精度与效率。
起伏地形快速算法： 针对起伏地形观测，传统方法需要在地形最高点与最低点之间离散出多个高程面，并对每个高程面重复计算整个模型空间的贡献，计算成本随高程面数量线性增加。研究团队提出了一种快速算法。通过分析傅里叶域表达式中垂直积分项的数学结构，他们发现可以将总积分区域巧妙分解为三个部分：当前观测面以上的部分、当前观测面以下至模型区域顶部的部分、以及模型区域以下的部分。其中，模型区域以下部分的贡献对所有观测面是相同的，只需计算一次并存储结果。当计算从地形最低点向最高点逐层进行时，当前高程面以上的积分区域可以从上一高程面的计算结果中继承或快速更新，无需重复进行完整的数值积分。最后，通过三次样条插值获得实际起伏地形上的磁异常值。该算法显著减少了重复计算量，其效率优势随着地形离散高程面数量的增加而愈加明显。
3. 数值实验与性能验证 研究设计了四个不同的二维模型，对NUFFT方法的可靠性、精度、效率及适用性进行了系统测试。所有计算代码均使用Fortran 95编写，并在标准工作站上运行以比较计算时间。同时引入了相对均方根误差（Relative Root Mean Square Error, RRMS）作为精度评价标准。
3.1 算法验证（模型一：简单矩形棱柱体） 为了验证新算法的正确性，研究首先设计了一个简单的二维矩形棱柱体模型，并计算其在平面观测条件下的磁力场和梯度张量。将NUFFT方法（分别采用51和201个波数）和标准FFT网格扩展法（采用不同的扩展系数ke）的计算结果与解析解进行对比。结果显示，NUFFT方法（n=201）与解析解吻合度极高，相对误差普遍低于0.3%。而标准FFT方法即使进行了网格扩展，在异常零点附近仍存在较大相对误差，并且强制周期性导致的边界效应明显，随着扩展系数ke增大（计算量剧增），该效应才逐渐减弱。此实验初步证实了NUFFT方法在精度上的优越性。
3.2 精度与效率平衡分析（模型一续用） 研究进一步深入比较了两种方法在达到相同精度（RRMS < 0.1%）时所需的计算时间。对于磁力场分量，NUFFT方法仅需13毫秒（波数n=251）即可满足精度要求，而标准FFT网格扩展法需要236毫秒（ke=15），NUFFT方法效率高出约17.1倍。对于磁梯度张量，标准FFT方法在ke=3时需69毫秒达到精度要求，NUFFT方法效率优势为4.3倍。这清晰地表明，在处理磁力场正演问题时，NUFFT方法在计算效率上具有显著优势。
3.3 起伏地形快速算法测试（模型二：复杂山谷地形） 为测试所提出的起伏地形快速算法的性能，研究构建了一个具有复杂山谷地形的模型，其下埋藏一个矩形磁性体。通过改变离散地形的插值平面数量，对比了常规累加算法与快速算法的计算时间。当插值平面数量较少（如5个）时，两种算法效率相当。但当插值平面数量增加到71个时，常规算法耗时167毫秒，而快速算法仅需58毫秒，效率提升显著。同时，研究表明当插值平面数量超过41个后，磁异常的计算精度（RRMS）趋于稳定且小于0.1%，验证了快速算法在保证精度的同时大幅提升了起伏地形建模的计算效率。
3.4 复杂地形影响研究（模型三：半空间山谷地形模型） 研究设计了一个半空间山谷地形模型，该模型的特点在于地形本身具有磁化率，且地形下伏背景介质磁化率不为零。通过计算不同观测高度（从地形表面向上延拓）的总场磁异常，研究直观展示了地形本身磁化引起的显著异常。结果强调，如果不考虑地形影响或仅简单地进行向上延拓，地形引起的异常会与目标地质体的异常叠加，可能导致错误的解释结论。这凸显了在磁异常反演解释中必须仔细考虑地形效应的重要性。
3.5 真实地形应用测试（模型四：基于数字高程模型的真实地形） 最后，研究将方法应用于中国云南-贵州高原区域的真实数字高程模型（DEM）。在复杂的真实地形下，设置了两个尺寸相同但埋深不同的磁性体。计算得到的总场磁异常清晰地显示：距离起伏观测点更近的右侧磁性体产生的异常幅度更大；地形的剧烈变化导致了磁异常的局部波动。异常的正负值边界与磁性体的边界对应良好，符合正演建模的基本规律。该实例有力地证明了该方法对实际复杂地形具有很强的适用性和稳健性。
三、 研究的主要结果 1. 算法正确性验证结果： 对于简单矩形模型，NUFFT方法（n=201）的计算结果与解析解高度一致，相对误差极小，证实了该算法理论推导和程序实现的正确性。标准FFT方法则表现出明显的边界效应和零点附近误差。 2. 精度与效率对比结果： 在达到相同计算精度（RRMS<0.1%）的前提下，对于磁力场正演，NUFFT方法比标准FFT网格扩展法快17.1倍；对于磁梯度张量正演，快4.3倍。这证明了NUFFT方法在实现高精度模拟方面具有更高的计算效率。 3. **起伏地形快速算法性能结果：** 快速算法在处理复杂地形时效率优势明显，尤其当离散高程面数量较多时，计算时间远少于常规累加算法，且能保持高精度（插值平面>41时RRMS稳定<0.1%）。 4. 地形效应影响结果： 数值模型清楚地表明，地形（尤其是具有磁性的地形）本身会产生可观的磁异常，这种异常无法通过简单的数据向上延拓来消除。这为磁法勘探数据的处理和解释提供了重要警示。 5. 真实地形建模结果： 基于真实DEM数据的建模成功获得了合理的磁异常分布，异常特征与模型设置（磁性体位置、地形起伏）相符，证明了该方法处理实际复杂地质问题的能力。
这些结果层层递进，从验证基础算法，到量化对比性能优势，再到展示对复杂及真实地形问题的处理能力，系统地支持了研究提出的核心论点：基于NUFFT的傅里叶域正演方法是一种高效、高精度且适用于复杂场景的磁异常建模工具。
四、 研究的结论与价值 本研究提出并成功验证了一种在傅里叶域中利用非均匀快速傅里叶变换进行二维磁异常正演建模的新方法。研究得出结论：与传统的标准FFT网格扩展法相比，NUFFT方法通过允许傅里叶域的非均匀采样，能够有效最小化边缘效应的影响，在保证数值精度的同时，显著提高了计算效率。同时，针对起伏地形设计的快速算法进一步降低了复杂地形条件下的计算成本。
该研究的科学价值在于：1）为地球物理磁法正演模拟提供了一种新的、更优的数值计算框架；2）深入分析了傅里叶域方法的误差来源，并通过波数优化和算法创新提出了解决方案；3）明确揭示了地形效应在磁异常中的重要性，对数据解释具有指导意义。其应用价值体现在：该方法能够高效、精确地模拟复杂形状、非均匀磁化地质体在起伏地形下产生的磁异常，可为矿产勘查、地质填图等实际工作中的高精度正反演解释提供强有力的工具。
五、 研究的亮点 1. 方法创新性： 首次将1D类型3的NUFFT系统性地应用于二维磁异常傅里叶域正演计算中，解决了传统FFT方法因强制周期性和均匀网格要求带来的精度与效率矛盾。 2. 关键技术优化： 提出了基于对数间隔的波数优化选择方案，以及针对起伏地形观测的快速积分算法，这两项优化显著提升了整体方法的实用性和计算性能。 3. 系统性的验证体系： 通过从简单模型到复杂模型、从理论对比到真实数据应用的一系列精心设计的数值实验，全面、量化地验证了方法的可靠性、优越性和稳健性。 4. 对实际问题的关注： 研究不仅关注算法本身，还通过模型实验强调了地形效应这一在实际勘探中至关重要但易被忽视的因素，提升了研究的现实意义。