A Spectral-Domain Pseudo-Inverse Method for True 3D Gravity Inversion

技术摘要：一种用于真 3D 重力反演的谱域伪逆方法

问题陈述
本文探讨了 3D 重力反演中固有的基本挑战：欠定性、不稳定性和深度分辨率差。传统的地面观测仅提供 2D 平面（$z=0$）上的数据，缺乏解析 3D 地下密度分布 $\rho(x,y,z)$ 所需的垂直波数信息。这导致独立约束数量不足。此外，重力场随深度增加而迅速衰减，导致来自深部源的信号受到严重削弱，使得反问题呈现病态特征，且对噪声高度敏感。现有的方法（如空间域迭代正则化，例如 Tikhonov 正则化）计算成本高昂且对初始模型敏感，而频域方法通常局限于界面重建而非完整的体积反演。

方法论
作者提出了一种基于酉可对角化系统谱域伪逆理论的“真 3D”反演方法。该方法论分为两个不同的理论层面：

1. 通过解析延拓进行数据空间补全：
   为了满足 3D 输入数据对谱域伪逆的要求，本文利用了重力异常 $g_z$ 在无源空气层（$z < 0$）中满足拉普拉斯方程（$\nabla^2 g_z = 0$）的事实。通过对地表数据 $g_z(x,y,0)$ 应用二维傅里叶变换，作者推导出了一个解析向上延续公式。这使得构建一个完整的、物理自洽的位于观测平面上方的 3D 数据体成为可能。至关重要的是，本文推导出了半空间谱的一般解析表达式：
   $$\tilde{g}^{obs}_z(k_x, k_y, k_z) = \tilde{g}_z(k_x, k_y, 0) \frac{1}{a - i k_z}$$
   其中 $a = \sqrt{k_x^2 + k_y^2}$。该表达式将水平谱与通用的垂直传播核分离，在不产生数值插值误差的情况下，提供了 3D 傅里叶变换所需的垂直波数采样。

2. 谱域伪逆反演：
   3D 重力的正向算子被证明在傅里叶域中是酉可对角化的。密度谱 $\tilde{\rho}$ 与重力谱 $\tilde{g}_z$ 之间的关系表示为：
   $$\tilde{g}_z(\mathbf{k}) = \lambda(\mathbf{k}) \tilde{\rho}(\mathbf{k}), \quad \text{其中 } \lambda(\mathbf{k}) = \frac{-i 4\pi G k_z}{K^2}$$
   这里 $K^2 = k_x^2 + k_y^2 + k_z^2$。基于作者此前关于酉可对角化系统的理论，构造了一个稳定的广义逆滤波器：
   $$q_\alpha(\lambda) = \frac{\bar{\lambda}}{|\lambda|^2 + \alpha}$$
   其中 $\alpha > 0$ 是正则化参数。该滤波器应用于构建的 3D 数据谱以获得密度谱。最终的空间密度分布通过 3D 逆傅里叶变换恢复。由于数据谱的极点结构（特别是 $k_z = -ia$ 处的极点），$k_z$ 的逆变换是通过复平面内的围道积分实现的，对于 $z > 0$ 情况在下半平面闭合。

核心贡献

• 数据空间补全： 本文证明了通过拉普拉斯方程进行的向上延续可以从 2D 地表数据严谨地构建出 3D 数据体，从而缓解了问题的欠定性，并提供了 3D 谱分析所需的垂直波数采样。
• 闭式反演： 本文证明了 3D 重力正向算子是精确酉可对角化的。因此，反演被简化为一个确定性的闭式计算过程，包括一次正向 3D 傅里叶变换、使用推导出的滤波器进行谱缩放步骤，以及一次逆变换，从而消除了对迭代优化的需求。
• 稳定性与一致性： 所构造的逆滤波器被证明具有有界稳定性（$\|q_\alpha\|_\infty \leq 1/(2\sqrt{\alpha})$）和一致性（当 $\alpha \to 0^+$ 时收敛于 Moore–Penrose 广义逆）。
• 解析框架： 半空间谱的推导（公式 7）提供了一个统一的输入框架，它将水平分量与垂直分量分离，适用于任何复杂的地下模型下的任意地表重力数据。

结果
该方法使用均匀球体模型进行了验证。结果表明：

• 该方法能够正确恢复异常体的 3D 位置。
• 随着正则化参数 $\alpha \to 0^+$，该方法成功消除了格林函数扩散效应，恢复了源点处的奇异行为（狄拉克 $\delta$ 函数）。
• 对于有限的 $\alpha$，解是一个平滑的分布，且对整个空间的积分严格守恒源体的总质量。
• 该方法恢复的是扰动密度（相对于未知背景），其均值为零；绝对密度需要独立的先验信息。密度分布的定量恢复取决于 $\alpha$ 的选择：较小的 $\alpha$ 会产生更锐利的图像和更准确的内部值，趋近于源点处的发散极限。

意义与主张
本文声称建立了一个理论自洽的解析框架，用于“真 3D”重力反演。其意义在于将范式从迭代优化转向结构分解方法。通过证明该问题可以通过单次谱滤波操作来解决，该方法为卷积型问题的稳定反演提供了一个统一框架。作者指出，虽然该方法在数值上等同于零阶 Tikhonov 正则化，但在概念上有所不同，因为它直接源于对泊松方程酉可对角化结构的推导，而非源于优化原理。文章最后指出，尽管该方法在连续谱假设下达到了理想结果，但未来的工作将侧重于实际数据的应用、数值实现细节以及鲁棒性分析。