[论文解读] 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），在相邻层次之间平衡逼近误差与收敛半径，确保整个层次结构中的收敛性。
• 分析直接问题的弗雷chet可微性，并推导出与频率和划分方式相关的相关常数的界，该界在远离狄利克雷本征频率时成立。

实验结果

研究问题

• RQ1 在具有间断波速的赫姆霍兹反问题中，如何缓解随着子域数量增加而引起的稳定性常数指数增长？
• RQ2 压缩逼近速率需满足何种条件，才能在不稳定性加剧的情况下确保多层级迭代重建格式的收敛性？
• RQ3 即使在频率数据有限的情况下，多层级投影最速下降方法是否仍能从粗略逼近出发收敛到真实解？
• RQ4 当频率趋近于零时，收敛半径的行为如何？该行为是否可用于改善低频区域的稳定性？
• RQ5 在分段常数波速的分层重建中，逼近误差与稳定性常数之间最优权衡是什么？

主要发现

• 在 $L^2(\theta)$ 中建立了赫姆霍兹反问题的条件性利普希茨稳定性估计，稳定性常数随分段常数划分中子域数量的增加而呈指数增长。
• 通过分层压缩逼近，最优地界定了稳定性常数，并推导出一个压缩速率条件，确保稳定性常数保持受控。
• 只要压缩速率满足所推导条件，多层级迭代格式即收敛至真实解，当逼近误差 $\eta_n$ 趋近于零时成立。
• 迭代格式的收敛半径 $\rho$ 随频率 $\omega$ 趋近于零而趋于无穷大，从而在低频下实现稳定重建。
• 该方法对有限频率数据（包括任意低频）仍有效，只要数据避开狄利克雷本征频率。
• 该格式在区域划分细化下具有鲁棒性，通过公式 (4.13) 平衡逼近误差与收敛半径，确保在各层次间均能保证收敛。