[论文解读] Increasing stability in acoustic and elastic inverse source problems

Mozhgan Nora Entrekhabi, Victor Isakov|arXiv (Cornell University)|Aug 30, 2018

Numerical methods in inverse problems参考文献 16被引用 3

一句话总结

本文通过利用多频数据和波数域中的傅里叶分析，建立了声波与弹性波散射中反源问题的递增稳定性估计。通过应用傅里叶变换、解析延拓的界以及双曲系统的精确可观测性，作者推导出近乎利普希茨的稳定性界，且随着波数区间增大而显著改善，从而在亥姆霍兹系统与拉梅系统中显著提升了源恢复的分辨率与鲁棒性。

ABSTRACT

We study increasing stability in the inverse source problems for the Helmholtz equation and the classical Lame system from (minimal) boundary data at multiple wave numbers. By using the Fourier transform with respect to wave numbers, explicit bounds for analytic continuation of the data to larger wave numbers, the Hyugens' principle, and sharp bounds in the corresponding dynamical initial boundary value problems, increasing (with growing wave numbers interval) stability estimates for source terms are obtained.

研究动机与目标

• 为解决椭圆方程与方程组反源问题的不适定性，此类问题的标准对数稳定性会限制分辨率。
• 通过利用一系列波数下的边界测量，而非单一频率，以改进稳定性估计。
• 建立从近场数据恢复声源与弹性源的近乎利普希茨稳定性界。
• 将现有递增稳定性结果扩展至更复杂的弹性系统，该系统包含多种波型（P波与S波）。
• 为生物医学成像、地球物理学及天线合成等实际应用提供理论框架。

提出的方法

• 对波数进行傅里叶变换，将多频反问题转化为频域中的问题。
• 利用解析解延拓的精确界，将稳定性估计扩展至更高波数。
• 利用惠更斯原理与关联双曲初值边值问题的精确可观测性估计。
• 通过帕塞瓦尔恒等式与索伯列夫空间中的能量估计，推导稳定性估计，适用于源项。
• 采用亥姆霍兹分解，将弹性系统中的压强波（P波）与剪切波（S波）分量分离。
• 引入随波数区间 $K$ 增大而减小的对数稳定性项，从而实现近乎利普希茨行为。

实验结果

研究问题

• RQ1 能否利用多频边界数据，为亥姆霍兹方程的反源问题建立递增稳定性？
• RQ2 随着波数区间 $K$ 增大，源恢复的稳定性如何改善？
• RQ3 类似递增稳定性结果能否扩展至包含P波与S波的更复杂的动力弹性系统？
• RQ4 傅里叶分析与解析延拓在转换稳定性估计中起到何种作用？
• RQ5 所推导的界与单频情况下的标准对数稳定性相比如何？

主要发现

• 对于亥姆霍兹方程，$f_0$ 与 $f_1$ 的稳定性估计随 $K$ 改进，达到形式为 $C\big(\frac{M_1^2}{1+K^{4/3}|E_0|^{1/2}}\big)$ 的近乎利普希茨界，其中 $E_0 = -\ln \varepsilon_0$。
• 源项的 $H^1$-范数的稳定性界以 $K^{-4/3}$ 速率衰减，表明波数区间越大，稳定性显著提升。
• 对于弹性系统，类似递增稳定性已证明，其界的形式为 $C\big(\frac{M_{3}^2}{1+K^{4/3}|E_e|^{1/2}}\big)$，随着 $K$ 增大，表现出近乎利普希茨行为。
• 对数稳定性成分随 $K$ 增大而减小，剩余项成为主导，从而在数值重建中实现更高分辨率。
• 结果通过傅里叶分析、解析延拓与双曲系统精确可观测性的结合推导得出，提供了稳健的理论框架。
• 该方法适用于全边界数据与实际测量（如仅压力或仅位移），增强了在真实世界反问题中的适用性。