基于高效多重网格求解器的金属矿电磁数据三维反演研究学术报告
一、 研究作者、机构及发表信息 本项研究的主要作者包括：第一作者及通讯作者王永斐（清华大学数学科学系）、郭荣文（中南大学关键矿产研究与勘探国家重点实验室等多个单位）、柳卓（斯坦福大学地球与行星科学系）、杨顶辉（清华大学数学科学系）和王登康（中南大学）。该研究成果以“Three-dimensional inversion of metallic minerals electromagnetic data using efficient multigrid solver”为题，发表于《Transactions of Nonferrous Metals Society of China》期刊第36卷（2026年），页码为943-955。
二、 学术背景与研究目的 本研究属于地球物理勘探与计算地球物理学交叉领域，具体聚焦于大地电磁法（Electromagnetic, EM）三维反演这一前沿课题。矿产资源是国家经济发展和可持续增长的重要保障，而理解成矿过程、深部物质来源与运移路径对于提升资源潜力评估至关重要。电磁法作为资源勘探的利器，其三维反演相较于传统的一维和二维反演，能够提供更详细、更准确的地下电性结构信息。然而，三维反演面临的核心挑战在于其高昂的计算成本，包括巨大的内存需求和漫长的计算时间，这严重制约了其在实际勘探中的广泛应用。反演过程的核心是大量、反复的正演模拟计算，因此，发展高效、精确的正演求解器是突破三维电磁反演计算瓶颈的关键。
在此背景下，多重网格（Multigrid, MG）求解器因其优异的收敛特性以及在网格规模增大时计算时间呈线性增长的优势，近年来在地球物理应用中受到广泛关注。然而，关于MG求解器在反演中的数值性能，特别是将其集成到反演流程中能带来多大效益，相关文献记载甚少。因此，本研究旨在填补这一空白。研究的目标是开发一种基于高效MG正演求解器的三维电磁反演算法，显著提升反演计算效率，并通过合成模型测试和实际矿区数据应用，验证该算法在金属矿勘探中的有效性和实用性。
三、 详细研究流程与方法 本研究遵循“方法开发-数值验证-实际应用”的完整技术路线，具体流程如下：
1. 反演算法框架构建： 本研究采用非线性共轭梯度（Nonlinear Conjugate Gradient, NLCG）法作为反演核心优化算法。NLCG法的优势在于其无需直接存储庞大的雅可比矩阵，而是通过“伴随正演模拟”计算雅可比矩阵与向量的乘积，从而显著降低了大规模三维反演的内存需求。反演的目标函数定义为数据拟合差与模型粗糙度的加权和，通过最小化该目标函数来寻找最优的地下电阻率模型。
2. 高效正演与伴随求解器开发： 这是本研究的创新核心。正演模拟需要求解频域麦克斯韦方程，本研究采用有限差分法（Finite Difference Method, FDM）在交错网格上进行离散化。为高效求解由此产生的大型线性方程组，研究团队没有使用传统的基于Krylov子空间（如BiCGSTAB）的迭代求解器，而是自主研发并应用了一种基于双色平面高斯-赛德尔（Two-Color Plane Gauss-Seidel, GS）平滑器的多重网格（MG）求解器。 * MG求解器原理：该方法的核心思想是“粗细网格交替求解”。首先在细网格上进行平滑迭代，快速消除误差的高频（短波）分量；然后将剩余的平滑（低频、长波）误差限制（Restrict）到更粗的网格上。在粗网格上，求解问题规模变小，计算成本低廉。再将粗网格上的解通过延拓（Prolongation）插值回细网格，修正细网格上的解。这个过程可以以V-cycle、W-cycle或F-cycle等形式在多层网格间递归进行。 * 关键技术——双色平面GS平滑器：针对大地电磁正演问题中旋度-旋度方程存在零空间导致长周期收敛慢甚至发散的问题，以及网格拉伸带来的挑战，本研究采用了特定的平滑技术。将整个三维模型沿某个方向切割成一系列平行的平面（Plane）。将这些平面分为两种颜色（如红色和黑色），要求同种颜色的任意两个平面互不相邻。这样，属于同一种颜色的所有平面系统可以完全独立、并行地求解（如图4所示）。对于每个平面系统，其内部的未知电场分量（图4a中红色箭头）通过LU分解高效求解，而周围边界条件（灰色箭头）的值则来自前一次迭代的结果。这种设计极大地促进了并行计算，提升了计算效率。 * 在反演中的集成：在NLCG反演的每次迭代中，通常需要进行两次正演模拟（对应两种极化模式）和一次伴随方程求解。本研究的巧妙之处在于，将上述MG求解器同时用于正演方程和伴随方程的求解。并且，在计算一种极化模式的正演时进行的LU分解结果可以被存储并重复利用，用于求解另一种极化模式的正演及其伴随计算，从而进一步节省了计算时间（如图5流程图所示）。
3. 合成模型测试与数值性能评估： 为了验证算法的准确性和效率，研究设计了一个包含六个异常体的合成模型进行测试。模型背景为100 Ω·m的均匀半空间，其中嵌入了三个10 Ω·m和三个1000 Ω·m的异常块体，分为浅埋和深埋两组。正演模拟在0.1至100秒的8个周期上进行，设置了144个测点，并在数据中添加了3%的高斯随机噪声。 * 精度验证：将MG求解器的正演结果与业界广泛使用的ModEM软件（其使用块不完全LU分解预条件的BiCGSTAB求解器，简称blockILU-BiCGSTAB）的结果进行对比。在10秒周期沿一条测线的视电阻率和相位响应对比显示（图7），两者结果高度吻合，视电阻率相对差异小于0.021%，相位差异小于0.0001°，充分证明了MG求解器的计算精度。 * 效率对比：对比了所有周期上MG与blockILU-BiCGSTAB的总计算时间（表1）。结果显示，MG求解器将计算时间减少了69%至89%。特别值得注意的是，传统求解器的计算时间随周期增大而显著增加，而MG求解器的计算时间在不同周期下保持相对稳定，展现了其优异的鲁棒性。 * 反演测试：使用相同的合成数据，分别用本研究开发的MG反演代码（MG-Inversion）和ModEM进行三维反演。反演网格为128×128×128个单元格。两者均经过8次迭代后使均方根误差（RMS）从初始的3.56降至目标值1.05以下（MG-Inversion为1.01，ModEM为0.99）。从反演结果（图9）看，两种方法都能有效恢复异常体的位置和形状，但MG-Inversion恢复的异常体边界更为清晰。效率提升是最大亮点：ModEM完成反演总耗时42.27小时，每次迭代平均4.85小时；而MG-Inversion仅需11.05小时，每次迭代平均1.28小时，加速比达到3.83倍，这与摘要中提到的“提升约3倍”相符。
4. 实际数据应用——北衙金多金属矿区案例： 在验证算法有效性后，研究将其应用于中国云南著名的北衙金多金属矿区的实测大地电磁数据。该矿区位于金沙江-哀牢山富碱斑岩带，其成因存在走滑、大陆俯冲、岩石圈拆沉等多种理论争议。研究旨在通过三维电性结构成像，为理解成矿过程提供深部地球物理约束。 * 数据与反演参数：使用了44个宽频大地电磁测点，周期范围0.003-1370秒。相位张量分析（图12）表明该区地下结构具有强烈的三维特性。反演区域大小为205 km × 205 km × 104 km，离散化为64×64×64的网格。初始模型和参考模型均为100 Ω·m的均匀半空间。 * 反演结果：经过174次迭代，RMS从13.31降至2.92。最终获得的三维电阻率模型揭示了显著的电性结构（图14, 15, 16）：识别出三个主要的低阻体（C1, C2, C3）和一个覆盖在C1之上的高阻体。C1和C2从近地表延伸至约10公里深度，并在深部可能相连。在更深部（20-30公里），C2和C3可能合并为一个统一的低阻通道。
四、 主要研究结果 1. 算法精度得到验证：通过合成模型正演对比，MG求解器与成熟商业软件ModEM的结果高度一致，证明了其数值解的准确性。 2. 计算效率显著提升：合成模型测试表明，MG求解器本身比传统求解器快数倍。集成到NLCG反演框架后，完整的反演计算速度比ModEM快3.83倍，实现了“单次反演计算效率提升约3倍”的目标。这主要归功于MG算法本身的计算复杂度与网格规模的线性关系，以及双色平面GS平滑器带来的高效并行能力。 3. 反演效果良好：对合成数据的反演表明，MG-Inversion不仅能准确恢复地下异常，其模型边界分辨率甚至优于作为基准的ModEM。这可能是由于高效求解器允许在相同时间内进行更精细的迭代或采用不同的正则化策略。 4. 实际应用揭示深部结构：对北衙矿区的实际数据反演，获得了稳定的、拟合良好的三维电阻率模型。模型清晰地刻画了浅部与深部的电性异常体空间分布。结合地质背景，研究者对低阻体（C1, C2, C3）和高阻体的地质含义进行了解释：C2可能对应哀牢山-红河剪切带相关的构造水和部分熔融；C1可能对应古老的岩浆房及相关的金属硫化物和含盐流体，是北衙多金属矿的矿源；其上覆的高阻体是良好的储矿层。深部连通的低阻通道（C2/C3）可能代表了深部物质上涌的通道。这些结果为“矿床形成与哀牢山-红河剪切带密切相关，深部富金属流体沿断裂系统上升成矿”的观点提供了地球物理证据。
五、 研究结论与价值 本研究成功开发并验证了一种基于高效双色平面GS多重网格求解器的三维大地电磁反演方法。主要结论如下： 1. 方法学结论：将MG求解器集成到NLCG反演框架中是可行且高效的。该方法通过重用LU分解等技术，优化了正演与伴随问题的求解流程。 2. 数值性能结论：新算法在保证精度的前提下，将三维反演的计算效率提升了近4倍，有效缓解了计算瓶颈。 3. 应用结论：该算法能够稳定、有效地处理实际野外数据，获得具有地质意义的三维电性结构模型，为利用电磁法精细刻画金属矿集区深部结构提供了有力的工具。
本研究的价值体现在： * 科学价值：推动了计算地球物理中高效数值算法与反演理论的结合，为大规模三维电磁反演的实际应用提供了新的技术方案。 * 应用价值：显著降低了大面积、高精度三维电磁勘探的时间与计算资源成本，使得利用三维反演进行深部找矿和构造研究变得更加可行，对矿产资源勘探具有直接的促进作用。
六、 研究亮点 1. 算法创新性强：首次将一种特定的、针对地球物理问题优化的双色平面GS多重网格求解器系统性地集成到三维大地电磁NLCG反演中，并详细评估了其在完整反演流程中的性能增益。 2. 效率提升显著且经过严格验证：研究不仅提出了新算法，还通过与国际主流软件ModEM在相同硬件平台上的“头对头”对比测试，用详实的数据（计算时间、加速比、RMS收敛曲线、模型切片对比）定量化地证明了其效率优势（3.83倍加速）和至少同等的精度与效果。 3. 研究链条完整：从方法原理阐述、算法实现、合成数据验证到实际案例应用，构成了一个完整的方法学研究与应用示范闭环，增强了结论的说服力。 4. 解决了实际勘探痛点：直面三维电磁反演计算成本高的行业难题，提出的方案具有明确的工程实用价值，有助于推动三维电磁勘探技术的普及。
七、 其他有价值内容 本研究在引言部分对三维电磁反演的发展现状进行了精炼的综述，涵盖了正演方法（积分方程法、有限元法、有限体积法、有限差分法）、主流反演算法（高斯-牛顿法、Occam法、拟牛顿法、非线性共轭梯度法）以及一些提高分辨率与效率的新思路（如多值反演、子域分解与最小熵约束），为读者提供了清晰的学术背景定位。此外，文中对相位张量三维性指标（偏斜角β）的分析，是处理复杂三维数据时一个重要的定性分析步骤，体现了严谨的数据处理流程。