论文信息：Li, D., Liang, Q., Du, J., Sun, S., Zhang, Y., & Chen, C. (2020). Transforming Total‐Field Magnetic Anomalies Into Three Components Using Dual‐Layer Equivalent Sources. Geophysical Research Letters, 47, e2019GL084607. https://doi.org/10.1029/2019GL084607

磁法勘探中大都是单分量、标量的观测，而地下磁性体的信息其实藏在更完整的矢量场里。本文由中国地质大学（武汉）与华中科技大学、浙江大学等单位的同行完成，本文作者之一张壹参与了其中浙江大学一方的合作。我们提出的双层等效源方法，用"浅层+深层"两层等效源结合预条件反演，把日常测量得到的总场磁异常 ΔB 稳定地转换为三分量磁场（Bx、By、Bz），尤其在长波长场上精度得到明显提升。下面围绕这篇 GRL 研究短文，梳理它的动机、方法与主要结论。

背景与动机：标量 ΔB 的局限

野外磁测（地面、航空、海洋）最常见到的输出是总场磁异常 ΔB——扣除背景地磁场 B0 后，反映地下磁性体分布的那个标量。ΔB 虽然对磁性体敏感，却同时强烈依赖于岩石磁化强度、磁性体的产状方向以及所处的纬度。换句话说，同一个地下模型在不同位置、不同方向的测量下，ΔB 的形态会有很大差异，这给数据处理与地质解释带来麻烦。

一个行之有效的破解办法，是把标量的 ΔB 转换为磁场的三个方向分量 Bx、By、Bz。矢量化的磁场一方面能压低"磁化方向对异常形态的影响"，一方面能压低反演与解释中的非唯一性，为后续正演与反演提供更干净的输入。

问题在于：传统转换手段（如快速傅里叶变换 FFT）要求观测位于水平平面且数据点在规则网格上。可实际野外数据往往落在起伏地形、不规则分布的测点上，FFT 类方法就"水土不服"了。而等效源方法正好能在任意曲面、任意分布上处理位场数据——这是它被广泛采用的根本原因。

方法：从单层等效源到双层等效源

等效源的思想很简单：在观测面之下布设一层假想的离散小体（等效源单元），通过反演确定每个单元的物性参数（这里为磁化强度），使得这一层假想源正演出的场能拟合观测数据；此后在这个等效源上做任意正演，即可得到上延、下延、分量转换等结果。

然而传统的单层等效源有一个致命短板：经验上（Dampney, 1969）等效源层应放在观测点网格间距的 2–6 倍深度处，它主要用来刻画短波长信息；而位场异常（含磁异常）总是同时包含长、短波长成分。一个贴近地面的单层等效源，无法充分捕获由深部源产生的长波长异常，于是当把 ΔB 转换为三分量、尤其做向上延拓或提取长波长分量时，误差就会显著放大，边界效应也难以压制。

为此本文引入双层等效源（图1）：

  • 浅层：紧贴观测面下方、随地形起伏的曲面/层面，由小而密的矩形棱柱单元组成，负责恢复短波长场；
  • 深层：位于地壳较深位置、置于一个平面上，单元尺度更大，负责捕获长波长信号。深层深度是一个关键参数，可通过磁数据的功率谱分析估计出一个平均深度（本此研究实测中由谱分析定在 5 km）。
三种方法重建三分量误差对比
图1　合成数据上三种方法重建三分量的绝对误差对比（引自论文 Figure 4，为论文核心结果图）

第1–4列分别为 ΔB、By、Bx、Bz 的重建误差分布；(a)–(d) 单层等效源、(e)–(h) 双层等效源（无 PCG）、(i)–(l) 双层等效源（含 PCG）。可见单层源在三分量上误差显著，加入深层源后明显改善，再引入预条件（PCG）后边界效应进一步被压制。

反演的目标是线性系统 d = G·m（d 为 ΔB 观测、m 为等效源磁化参数、G 为核函数矩阵）。难点在于：磁场随与源的距离呈快速衰减，若不加以处理，反演出的物性会过度集中于浅层，从而低估深层源的贡献——这正是深层的"存在感"被淹没的根源。本文借鉴磁法三维反演中的深度加权思想（Li & Oldenburg, 1996；Pilkington, 1997），构造预条件矩阵 S = z^β I，引入预条件共轭梯度（PCG）求解。对于棱柱单元，由于衰减不仅取决于深度还取决于单元尺度与长宽比，采用 z^3（偶极近似）并不合适，数值实验表明 β = 2 与棱柱核函数的衰减曲线最吻合（论文补充图 S3）。预条件矩阵只作用在与深层单元对应的那部分核函数矩阵块上，从而有针对性地"抬升"深层的有效权重。

主要结果：三分量与长波长场精度显著提升

在合成数据上（三个埋深不同的磁性体，观测位于起伏的任意曲面上），对三种方案做了系统的误差对比（图1）：

  • 单层等效源在重建观测面附近的 ΔB 时尚可，但无法高精度地重建三分量（图1 a–d 各分量误差显著）；
  • 加入深层等效源后，三分量重建精度明显改善（图1 e–h）；
  • 再叠加预条件（PCG）后表现最佳：边界效应被显著削减（论文补充图 S4 显示上延结果的边界伪差几乎消失），数据区内最大的误差（位于西南角）经本转换方法后可压低到 1% 以内。

在实测数据（广西社桐的隐伏岩体与地下磁性基底）上，方法同样稳定：用有限的地面观测即可重建出不同高度面上的 ΔB、By、Bx、Bz 四个量（图2），更好地刻画了当地磁场随高度的变化特征。这说明双层等效源+预条件反演对真实、起伏地形上的数据同样有效。

实测数据重建三分量磁场
图2　实测数据（广西社桐）重建的 ΔB 与三分量磁场（引自论文 Figure 6）

(a)–(d) 地面上的 ΔB、By、Bx、Bz；(e)–(h) 1 km 高度平面上的对应结果。该区下方存在隐伏岩体与磁性基底，重建结果更好地还原了磁场随高度变化的特征。

结论与讨论

这项工作把"浅层+深层"的双层等效源与深度加权/PCG 预条件反演结合起来，一举解决了两个纠缠的问题：一是单层等效源天然缺失长波长能力，二是反演中物性向浅层聚集而低估深层。加入深层等效源并配合预条件矩阵后，三分量磁场的重建精度大幅提高，尤其是长波长场。

当然，方法也有明确的代价：计算时间与内存开销较大。当数据规模很大或要求拟合精度非常高时，需要高性能计算资源与较长运行时间——这也是作者指出的未来改进方向（加速反演算法）。同时，由于假设不存在剩余磁化、忽略退磁效应且每个单元磁化强度恒定，方法适用的是感应磁化为主的场景。

值得说明的是，这一方法的适用范围并不仅限于地面磁测：它天然适合任意分布、任意地形的数据，因此可推广到位场建模、多源数据融合，以及填补缺失数据区（如卫星与航空磁测、重力数据的空白带），具有较好的工程与应用价值。

延伸阅读与引用

本文为对该论文的中文解读，图片均引自论文原文。