磁异常

论文信息：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）后边界效应进一步被压制。 ...