[論文レビュー] 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 の場合に a = 1/2 よりも C(δ) が大きくなることが観察された。
• この結果は、磁気ポテンシャルデータのみからでは、磁化されたリソスフィア層の厚さは一意に特定できない（系理3.8による）。