使用SAI‑Krylov的降阶谱元方法用于高精度且快速的三维瞬态电磁正演模拟

用更快的虚拟传感器窥探地下

寻找埋藏的矿体、地下水或隐伏断层，常常依赖地球电磁场中微弱的波动。要把这些波动转化为地下结构的图像，需要在超大型三维模型上进行大量数值计算，这往往非常缓慢。本文提出了一种新的电磁测量模拟方法，在保持现场勘探所需精度的同时，大幅降低计算时间和内存开销。

为什么电磁测量的模拟很难

瞬态电磁（TEM）方法通过在地表的线圈或接地电缆中注入电流脉冲工作。当电流断开时，地下会感应出涡流并随时间衰减，这些衰减携带有关岩石电导率和隐伏矿体的线索。为了预测接收器的观测值，科学家需要在三维网格上求解麦克斯韦方程。传统方法在每个单元内采用仅线性变化的低阶基函数，这类低阶近似易于实现，但难以在不爆炸性增加未知量的情况下捕捉细节，也会导致需要通过成千上万的小时间步进进行昂贵的时间推进。

用更少单元实现更精细的网格

作者采用了高阶谱元方法，用更灵活的数学形状丰富每个网格单元，而不是依赖直线变化。实际效果是，模型能够在相同网格下更平滑地表示电场变化，或以远少于传统方法的单元数达到给定精度。他们精心设计了沿边界放置未知量的方式、将物理单元映射到简单参考立方体的过程，以及如何使用特殊的高斯‑洛巴托点和权重。其中一项关键技巧称为约化积分，它在略微放松某些积分精确性的同时，使得关键的电导率矩阵严格对角化。这显著增加了稀疏性，降低了内存使用和线性系统求解的代价，同时在感兴趣的阶数范围内保持高精度。

在不丢失物理量纲的情况下压缩数学问题

即使使用更精细的单元，按传统方式逐步推进时间仍然代价高昂，因为时间步变化往往迫使对一个巨大矩阵重新分解。作者将问题改写为连续时间系统，并采用基于移位‑反演（shift‑and‑invert）Krylov子空间算法的模型降阶技术。通俗地说，他们将完整的大系统投影到一组更小的“代表性”模态上，这些模态能捕捉电磁场的实际演化。一个巧妙的谱变换将重要模态聚拢，使算法快速收敛。关键在于，方法能在运行时基于廉价的残差度量自适应决定降阶子空间的大小，作为停止准则。因此，求解器仅需一次主要的矩阵分解加上少量回代，就能在不逐步遍历每一时刻的情况下评估任意时刻的电磁响应。

将方法付诸测试

为了评估精度和效率，作者对若干基准情形进行了模拟。对于具有已知解析行为的简单半空间，他们比较了不同多项式阶数和积分策略与标准向后时间步进方案的结果。从二阶及以上，误差稳定在约百分之一或更低，并且在新的降阶模型下平均误差可降至约百分之零点五，同时在高阶运行中实现约16到20倍的加速，而内存需求并未显著增加。在具有强电阻率对比的分层模型中，谱元加Krylov方法在整个时间窗口内误差维持在约百分之一以内，并优于依赖低阶单元和更耗散时间步进的传统有限元基线。最后，在一个完全三维的硫化物矿体情形中，该方法能够追踪由接地线和线圈源激发的感应场如何传播、与复杂的导体体相互作用，并最终以高空间分辨率描绘出埋藏目标。

这对地下勘探意味着什么

对于依赖精细三维电磁建模的地球物理勘测，这项工作提供了一种兼顾高保真与高速度的途径。通过将高阶谱元与精心设计的模型降阶方案结合，作者展示了在实际TEM响应下可以实现亚百分比误差，同时将计算时间相比可信基线缩减一个量级。通俗地说，这意味着从原始测量设计到可解释的地下成像的周转更快，使得使用基于严格物理的模拟而非过度简化的捷径来勘探矿产、评估地质灾害或追踪地下水变得更可行。

引用: Fan, Y., Lu, K., Huang, Y. et al. Model-order-reduced spectral-element method for high-accuracy and fast 3-D transient electromagnetic forward modeling with SAI-Krylov. Sci Rep 16, 13356 (2026). https://doi.org/10.1038/s41598-026-44053-y

关键词: 瞬态电磁, 谱元方法, 模型降阶, 地球物探, 数值模拟