[论文解读] Recovery of a source term or a speed with one measurement and applications

Plamen Stefanov, Günther Uhlmann|arXiv (Cornell University)|Mar 6, 2011

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

一句话总结

本文在各向异性介质中，通过单个边界测量，建立了恢复源项或声速的精确唯一性与稳定性条件。结合Carleman估计与几何论证——特别是严格凸曲面的叶状结构——证明了在最小几何与正则性假设下，声速可被唯一且稳定地重构，直接应用于光声断层扫描。

ABSTRACT

We study the problem of recovery the source $a(t,x)F(x)$ in the wave equation in anisotropic medium with $a$ known so that $a(0,x) ot=0$ with a single measurement. We use Carleman estimates combined with geometric arguments and give sharp conditions for uniqueness. We also study the non-linear problem of recovery the sound speed in the equation $u_{tt} -c^2(x)Δu =0$ with one measurement. We give sharp conditions for stability, as well. An application to thermoacoustic tomography is also presented.

研究动机与目标

• 在各向异性介质中，通过单个边界测量，建立对波方程中源项 $ a(t,x)F(x) $ 的精确唯一恢复条件。
• 通过几何术语（特别是严格凸曲面的叶状结构）表述条件，拓展基于Carleman估计的反问题方法的应用范围。
• 证明从单个测量中恢复声速 $ c(x) $ 的反问题的局部与全局稳定性估计，且在有界 $ C^k $-范数下具有统一的界。
• 提供一个几何框架，以改进超bolic反问题中先前非精确的拟凸性与凸性条件。
• 将结果应用于光声断层扫描，表明声速可从初始压力数据的边界测量中稳定恢复。

提出的方法

• 使用时间依赖权函数的Carleman估计，适应问题的几何结构，避免依赖于单一全局拟凸函数。
• 通过一族严格凸的等值面 $ ext{dist}(x, ar{ heta}) = s $（其中 $ s o 0 $）实施逐步唯一延拓，将信息从边界传播至内部。
• 采用几何叶状结构论证：若边界 $ ar{ heta} $ 严格凸，且等值面在距离 $ T $ 内保持严格凸，则 $ F = 0 $ 在 $ ar{ heta} $ 的邻域内成立。
• 通过线性化与波方程的结构，将恢复 $ c(x) $ 的非线性反问题转化为恢复 $ F = ilde{c}^2 - c^2 $ 的线性问题，其中 $ a = abla ilde{u} $。
• 利用Dirichlet-to-Neumann映射 $ B $ 及其参数主 $ Q $，证明在条件 (51) 下，$ B $ 是阶数为 0 的经典拟微分算子，且主符号为椭圆。
• 通过紧致扰动论证建立稳定性：$ F $ 的解算子有有界性，且含紧致误差项，从而导出Fredholm型估计，表明在数据小扰动下具有稳定性。

实验结果

研究问题

• RQ1 在何种几何与分析条件下，可从单个边界测量唯一恢复波方程 $ u_{tt} - c^2(x) abla abla u = a(t,x)F(x) $ 中的源项 $ F(x) $？
• RQ2 声速 $ c(x) $ 能否从 $ [0,T] imes ar{ heta} $ 上波解的单个测量中稳定重构？$ T $ 与几何结构的精确条件是什么？
• RQ3 如何用严格凸曲面的几何叶状结构条件替代Carleman估计中的拟凸性条件，以改进先前非精确的假设？
• RQ4 反问题中 $ c(x) $ 的稳定性如何？稳定性常数如何依赖于 $ c $ 与 $ ilde{c} $ 的正则性与有界性？
• RQ5 结果能否应用于光声断层扫描？从初始压力数据中重构 $ c(x) $ 的稳定性如何？

主要发现

• 当 $ ext{dist}(x, ar{ heta}) o s o 0 $ 时，若等值面严格凸，且 $ T $ 足够小，则 $ F $ 在 $ ar{ heta} $ 的邻域内唯一确定。
• 条件 $ T o 0 $ 是精确的：若当 $ s o 0 $ 时等值面 $ ext{dist}(x, ar{ heta}) = s $ 不严格凸，则唯一性可能不成立。
• 通过紧致扰动论证证明了 $ c(x) $ 反问题的稳定性：若第二时间导数 $ ilde{u}_{tt} $ 接近 $ u_{tt} $，则 $ ilde{c} $ 在 $ L^2(K) $ 中接近 $ c $。
• 稳定性估计 $ orm{ ilde{c} - c}_{L^2(K)} o C orm{ ilde{u}_{tt} - u_{tt}}_{L^2([0,T] imes ar{ heta})} $ 成立，且 $ C $ 在 $ c $ 的 $ C^k $-有界性与椭圆性条件下一致，其中 $ k > n/2 $。
• 当 $ f o H^{k+1} $ 且 $ c o C^k $，$ k > n/2 $ 时，解 $ u $ 属于 $ C^2([0,T]; C^2(ar{ heta})) $，确保边界数据充分正则，定义良好。
• 该方法适用于光声断层扫描：在相同几何与正则性条件下，声速 $ c(x) $ 可从初始压力的单个测量中稳定恢复。