[论文解读] Unique reconstruction of simple magnetizations from their magnetic potential

Laurent Baratchart, Christian Gerhards+2 more|arXiv (Cornell University)|Mar 10, 2021

Numerical methods in inverse problems参考文献 42被引用 4

一句话总结

本文建立了在何种条件下，可通过磁势唯一重构体积体的磁化强度。证明了调和梯度与盒单纯向量场（在长方体上分段常数）可实现唯一反演，解决了真实地球物理模型中反磁静力学问题的非唯一性问题。

ABSTRACT

Inverse problems arising in (geo)magnetism are typically ill-posed, in particular {they exhibit non-uniqueness}. Nevertheless, there exist nontrivial model spaces on which the problem is uniquely solvable. Our goal is here to describe such spaces that accommodate constraints suited for applications. In this paper we treat the inverse magnetization problem on a Lipschitz domain with fairly general topology. We characterize the subspace of $L^{2}$-vector fields that causes non-uniqueness, and identify a subspace of harmonic gradients on which the inversion becomes unique. This classification has consequences for applications and we present some of them in the context of geo-sciences. In the second part of the paper, we discuss the space of piecewise constant vector fields. This vector space is too large to make the inversion unique. But as we show, it contains a dense subspace in $L^2$ on which the problem becomes uniquely solvable, i.e., magnetizations from this subspace are uniquely determined by their magnetic potential.

研究动机与目标

• 解决(地球)磁学中反磁化问题产生的非唯一性问题。
• 识别使反问题唯一可解的模型空间——特别是调和梯度与分段常数向量场。
• 提供磁化强度可从其磁势唯一重构的理论与实际条件。
• 探讨这些结果在实际地球物理场景（如岩石圈磁化与岩心样品分析）中的适用性。
• 证明尽管分段常数场在一般情况下不唯一，但其在稠密子空间（盒单纯场）中可确保唯一重构。

提出的方法

• 利用霍奇分解刻画不可见磁化强度的空间，识别出 H1₀(Ω) ⊕ Hdiv,0(Ω) 与调和梯度 ∇IΘ 为非唯一性的来源。
• 通过霍奇分解与调和函数的性质，证明只有磁化强度的调和梯度分量对域外磁势有贡献。
• 将盒单纯向量场定义为在长方体平行六面体上分段常数的场，其在 L²(Ω, ℝ³) 中构成稠密子空间。
• 通过定理 4.3 建立在盒单纯向量场空间上反问题的唯一可解性，即使事先未知盒子的数量、大小或位置。
• 利用离散化磁势算子的数值实验，验证稳定性常数 C(δ) 随分辨率 δ 呈指数增长。
• 使用谱范数与高精度矩阵求逆算法，计算离散化势算子 PΩc←Ω 的逆算子范数。

实验结果

研究问题

• RQ1 在利普希茨域中，体积磁化强度的反磁化问题在何种条件下是唯一可解的？
• RQ2 反问题的非唯一性是否可完全通过磁化强度场的霍奇分解来表征？
• RQ3 当假设磁化强度为分段常数时，是否可实现唯一重构？若可，是在何种子空间上？
• RQ4 磁化强度在域内局部化如何影响反问题的稳定性？
• RQ5 在离散逼近中，随着空间分辨率 δ 增大，稳定性常数 C(δ) 的行为如何？

主要发现

• 不可见磁化强度的空间被表征为 H1₀(Ω) ⊕ Hdiv,0(Ω)，可能还需补充调和梯度 ∇IΘ，具体取决于域的拓扑结构。
• 只有磁化强度的调和梯度分量对域外磁势有贡献，这意味着仅凭势数据无法探测非调和部分。
• 在盒单纯向量场空间——即在长方体平行六面体上分段常数的场——上，反问题具有唯一可解性，且无需事先知晓盒子的配置。
• 数值实验表明，稳定性常数 C(δ) 随分辨率 δ 呈指数增长，对于全支撑磁化强度，有 C(δ) ≈ exp(−7.933 + 4.562δ⁻⁰.⁸⁰⁴⁴)。
• 磁化强度在域的子集内局部化会增加不稳定性界，如 a = 1/4 时的 C(δ) 大于 a = 1/2 时。
• 该结果表明，根据推论 3.8，仅凭磁势数据无法唯一确定磁化岩石圈层的厚度。