等效源

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

论文信息：Zhang, Y., Xu, Y., Mooney, W. D., & Chen, C. (2021). Local separation of potential field anomalies using equivalent sources: application for the 3-D structure of mantle uplift beneath Von Kármán crater, the Moon. Geophysical Journal International, 227(3), 1612–1623. https://doi.org/10.1093/gji/ggab307 位场（重力与磁力）观测是研究地下密度与磁性结构的重要途径，但在实际资料中，浅部小尺度目标的异常总是叠加在深部大尺度场的背景之上，如何把它们可靠地分离开来一直是位场资料处理与解释中的经典难题。本文以一作张壹为代表的研究，提出了一种基于等效源技术 + 迭代反演的位场异常局部分离新方法，并将其应用于月球 Von Kármán 撞击坑（VKC，嫦娥四号着陆区）下方地幔隆起的 3D 结构重建。下面围绕这篇 Geophysical Journal International 论文，梳理其动机、方法要点与主要结论。 为什么要做"局部分离"：区域—剩余异常分离的局限 观测位场数据是多个互相叠加的地下源体所产生的重力或磁效应之和，各源体的异常往往难以单独提取。传统上，我们把反映深部大尺度背景的异常称为区域异常，把浅部小尺度目标所对应的异常称为剩余异常，二者的分离（regional–residual separation）是位场解释中极有价值的一步。几十年来发展出了大量方法，大体可分两类： 基于数据的方法：利用区域与剩余异常在图形特征或频率上的差异来分离，例如多项式拟合（Beltrao et al., 1991）、最小曲率法（Mickus et al., 1991）、有限元（Mallick & Sharma, 1999）、三维主成分与纹理分析（Zhang et al., 2009）、双维经验模态分解（Hou et al., 2012）、维纳滤波（Pawlowski & Hansen, 1990）、独立成分分析（Forootan & Kusche, 2012）以及各类小波分析（Fedi & Quarta, 1998；Xu et al., 2009）等。 基于模型的方法：利用位场的固有特性来分离，例如向上延拓作为标准分离滤波（Jacobsen, 1987；Zeng et al., 2007）、优化滤波参数（Pilkington & Cowan, 2006）、以及基于**等效层（equivalent layer）**概念的带通滤波等（Pawlowski, 1994；Guo et al., 2013）。 这些方法虽多，却普遍存在几个不足：其一，对被分离异常的质量缺乏可量化的评价指标，选择何种分离结果仍带有较强的主观性；其二，无论是小波分析还是球谐分析本质上都是信号处理手段，所生成异常的空间物理合理性（physical plausibility）并不总能得到保证。更重要的是，正如 Li & Oldenburg (1998) 指出的，以往绝大多数分离方法针对的是埋深明显不同的源体所产生、波长差异较大的异常；而当我们想要从观测中分离出波长相近的局部异常时（例如几处相邻源体的异常叠加），基于频率差异的方法便难以奏效，这正构成了本文要解决的"局部分离"问题。 ...