[论文解读] A Lipschitz stable reconstruction formula for the inverse problem for the wave equation

Shitao Liu, Lauri Oksanen|arXiv (Cornell University)|Oct 3, 2012

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

一句话总结

本文提出了一种全局可重构、局部Lipschitz稳定的反波速问题公式，结合了边界控制法与复几何光学解。研究证明，在黎曼流形 $ (M, c^{-2}dx^2) $ 上存在一个无临界点的严格凸函数时，$ c^{-2} $ 的低频分量可在满足几何条件的前提下实现Lipschitz稳定性重构。

ABSTRACT

We consider the problem to reconstruct a wave speed $c \in C^\infty(M)$ in a domain $M \subset \R^n$ from acoustic boundary measurements modelled by the hyperbolic Dirichlet-to-Neumann map $Λ$. We introduce a reconstruction formula for $c$ that is based on the Boundary Control method and incorporates features also from the complex geometric optics solutions approach. Moreover, we show that the reconstruction formula is locally Lipschitz stable for a low frequency component of $c^{-2}$ under the assumption that the Riemannian manifold $(M, c^{-2} dx^2)$ has a strictly convex function with no critical points. That is, we show that for all bounded $C^2$ neighborhoods $U$ of $c$, there is a $C^1$ neighborhood $V$ of $c$ and constants $C, R > 0$ such that |\F\ll( ilde c^{-2} - c^{-2} r)(ξ)| \le C e^{2R |ξ|} orm{ ilde Λ- Λ}_*, \quad ξ\in \R^n, for all $ ilde c \in U \cap V$, where $ ilde Λ$ is the Dirichlet-to-Neumann map corresponding to the wave speed $ ilde c$ and $ orm{\cdot}_*$ is a norm capturing certain regularity properties of the Dirichlet-to-Neumann maps.

研究动机与目标

• 开发一种用于从Dirichlet-to-Neumann数据反演波速 $ c $ 的全局可重构方法。
• 在反波速问题中实现 $ c^{-2} $ 的局部Lipschitz稳定性，优于通常的对数或Hölder型稳定性。
• 在几何条件——稳定可观测性——下建立稳定性，该条件由黎曼流形 $ (M, c^{-2}dx^2) $ 上存在无临界点的严格凸函数所表征。
• 弥合边界控制法与复几何光学方法之间的差距，以获得更强的稳定性估计。
• 为通过非线性Landweber迭代实现高效数值反演提供基础，该方法要求Lipschitz稳定性以保证线性收敛。

提出的方法

• 结合边界控制法与形式为 $ e^{i( ho + i heta)ullet x/2} $ 的复几何光学（CGO）解，其中 $ | ho| = | heta| $，$ ho ot heta $，构造重构公式。
• 使用调和指数函数作为测试解，利用其振荡结构提取关于 $ c^{-2} $ 的信息。
• 应用改进的边界控制法，结合可观测性估计与紧致流形上的Friedrichs不等式进行能量控制。
• 在Dirichlet-to-Neumann算子上定义范数 $ orm{ullet}_{*} $，以捕捉正则性与稳定性，该范数结合了 $ orm{ullet}_{L^2} $、$ orm{ullet}_{H^1} $ 与因果反射算子 $ ilde{ ho} $。
• 通过傅里叶域中重构误差的先验界建立稳定性：$ | ilde{F}( ilde{c}^{-2} - c^{-2})( ho)| riangleq ext{const} imes e^{2R| ho|} orm{ ilde{ ho} - ho}_{*} $，其中 $ ilde{c} $ 属于 $ C^2 $-邻域。
• 依赖于在 $ (M, c^{-2}dx^2) $ 上存在一个无临界点的严格凸函数 $ ho $，以确保稳定可观测性，并在 $ c $ 的邻域内对几何量保持一致有界。

实验结果

研究问题

• RQ1 能否构造一个波速 $ c $ 的全局重构公式，使其在Dirichlet-to-Neumann映射的小扰动下保持稳定？
• RQ2 是否可能在反波速问题中实现局部Lipschitz稳定性，而非仅对数或Hölder型稳定性？
• RQ3 黎曼流形 $ (M, c^{-2}dx^2) $ 上的何种几何条件足以保证此类Lipschitz稳定性？
• RQ4 如何通过引入复几何光学解来增强边界控制法，以提升稳定性与重构质量？
• RQ5 所提出的重构公式是否允许在非线性Landweber迭代等迭代求解器中实现收敛性保证？

主要发现

• 本文通过边界控制法与复几何光学解的结合，建立了 $ c^{-2} $ 的全局可重构公式。
• 在稳定可观测性假设下（即在 $ (M, c^{-2}dx^2) $ 上存在无临界点的严格凸函数），重构具有局部Lipschitz稳定性。
• 具体而言，对于所有属于有界 $ C^2 $-邻域 $ U $ 的 $ ilde{c} $ 以及 $ c $ 的 $ C^1 $-邻域 $ V $，对所有 $ ho eq 0 $，不等式 $ | ilde{F}( ilde{c}^{-2} - c^{-2})( ho)| riangleq Ce^{2R| ho|} orm{ ilde{ ho} - ho}_{*} $ 成立，其中 $ C, R > 0 $。
• 范数 $ orm{ullet}_{*} $ 捕捉了Dirichlet-to-Neumann映射的正则性，是量化稳定性的关键。
• 该结果意味着反问题可通过非线性Landweber迭代实现高效数值求解，该方法在Lipschitz稳定性下可实现线性收敛。
• 该稳定性结果是紧致的，即若无几何假设，则Lipschitz稳定性不成立，如Stefanov与Uhlmann通过存在截留测地线的反例所证明。