[論文レビュー] Inverse boundary value problem for the Helmholtz equation: Multi-level approach and iterative reconstruction

Elena Beretta, Maarten V. de Hoop+2 more|arXiv (Cornell University)|Jun 10, 2014

Numerical methods in inverse problems参考文献 17被引用数 8

ひとこと要約

本稿では、ヘルムホルツ方程式の逆問題において、ディリクレ・トゥ・ノイマンデータから不連続な波動速度を反復的に再構成するための、バナッハ空間における多層非線形射影勾配降下法を提示する。階層的圧縮近似を通じて、安定性定数の最適な上限を得る条件的リプシッツ安定性推定を確立し、部分領域数の増加に伴う安定性定数の指数的増大に対しても収束を可能にする。

ABSTRACT

We study the inverse boundary value problem for the Helmholtz equation using the Dirichlet-to-Neumann map at selected frequency as the data. We develop an explicit reconstruction of the wavespeed using a multi-level nonlinear projected steepest descent iterative scheme in Banach spaces. We consider wavespeeds containing discontinuities. A conditional Lipschitz stability estimate for the inverse problem holds for wavespeeds of the form of a linear combination of piecewise constant functions with an underlying domain partitioning, and gives a framework in which the scheme converges. The stability constant grows exponentially as the number of subdomains in the domain partitioning increases. To mitigate this growth of the stability constant, we introduce hierarchical compressive approximations of the solution to the inverse problem with piecewise constant functions. We establish an optimal bound of the stability constant, which leads to a condition on the compression rate pertaining to these approximations.

研究の動機と目的

• 周波数領域のディリクレ・トゥ・ノイマンデータを用いて、不連続な波動速度を伴うヘルムホルツ方程式の逆境界値問題に対処する。
• 区分的定数波動速度再構成における領域分割の増加に伴う安定性定数の指数的増大を克服する。
• 階層的レベル間の近似誤差と収束半径のバランスを保つことで、収束を保証する多層反復スキームを開発する。
• 安定性定数を最適化し、頑健な再構成を可能にする圧縮近似レートに関する条件を確立する。
• 周波数が低いかゼロに近い場合でも、ディリクレ固有周波数を避けることで、安定かつ収束する反復再構成の枠組みを提供する。

提案手法

• 波動速度から $L^2(\theta)$ 内のディリクレ・トゥ・ノイマン写像への非線形作用素として逆問題を定式化し、有限周波数の境界測定値からのデータを用いる。
• 複素幾何的オプティクス解を用いて、$L^2(\theta)$ 内で条件的リプシッツ安定性推定を導出する。安定性定数は、区分的定数分割における部分領域数に依存する。
• より細分化された領域分割における区分的定数関数を用いた、圧縮近似の多スケール階層を導入する。
• バナッハ空間における射影勾配降下法を適用し、各細分化レベルに対応する有限次元部分空間に反復的更新を射影する。
• 隣接するレベル間の近似誤差と収束半径のバランスを保つ十分条件（式 (4.13)）を導出し、階層全体での収束を保証する。
• 直接問題のフレシェ微分可能性を分析し、周波数および分割に依存する定数の上限を導出する。これはディリクレ固有周波数を除く範囲で有効である。

実験結果

リサーチクエスチョン

• RQ1 不連続な波動速度を伴う逆ヘルムホルツ問題において、部分領域数の増加に伴う安定性定数の指数的増大は、どのように軽減可能か？
• RQ2 圧縮近似レートにどのような条件を課すことで、不安定性が増す中でも多層反復再構成スキームの収束を保証できるか？
• RQ3 初期近似が粗い場合でも、有限周波数データに制限がある中で、多層射影勾配降下法が真の解に収束するか？
• RQ4 周波数がゼロに近づく際、収束半径 $\rho$ はどのように振る舞い、低周波数領域での安定性向上にどのように利用可能か？
• RQ5 区分的定数波動速度の階層的再構成において、近似誤差と安定性定数の最適なトレードオフは何か？

主な発見

• 区分的定数分割における部分領域数の増加に伴い、安定性定数が指数的に増大するが、$L^2(\theta)$ 内で逆ヘルムホルツ問題の条件的リプシッツ安定性推定が確立された。
• 階層的圧縮近似を通じて安定性定数が最適に制御され、安定性定数が制御されたままとなる圧縮レートに関する条件が導出された。
• 圧縮レートが導出された条件を満たす限り、近似誤差 $\eta_n$ がゼロに近づくにつれて、多層反復スキームは真の解に収束する。
• 反復スキームの収束半径 $\rho$ は周波数 $\omega$ がゼロに近づくにつれて無限大に近づき、低周波数領域での安定再構成を可能にする。
• データがディリクレ固有周波数を避ける限り、有限周波数データ（任意に低い周波数を含む）に対しても、本手法は有効である。
• 式 (4.13) を用いた近似誤差と収束半径のバランスにより、領域分割の細分化に対する耐性があり、全レベルにわたり収束が保証される。