[论文解读] Stability estimates in stationary inverse transport

Guillaume Bal, Alexandre Jollivet|ArXiv.org|Apr 8, 2008

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

一句话总结

该论文通过将边界数据映射到与速度向量正交的平面的几何变换，建立了从阿尔贝多算子出发的静止线性输运方程中空间依赖散射和吸收系数重建的 $L^1$-稳定性估计。关键贡献是在 $n \geq 3$ 维下实现了严格稳定性结果，扩展了以往仅限于空间独立散射系数的工作。

ABSTRACT

We study the stability of the reconstruction of the scattering and absorption coefficients in a stationary linear transport equation from knowledge of the full albedo operator in dimension $n\geq3$. The albedo operator is defined as the mapping from the incoming boundary conditions to the outgoing transport solution at the boundary of a compact and convex domain. The uniqueness of the reconstruction was proved in [M. Choulli-P. Stefanov, 1996 and 1999] and partial stability estimates were obtained in [J.-N. Wang, 1999] for spatially independent scattering coefficients. We generalize these results and prove an $L^1$-stability estimate for spatially dependent scattering coefficients.

研究动机与目标

• 证明静止线性输运方程中空间依赖散射和吸收系数重建的 $L^1$-稳定性估计。
• 扩展先前仅对空间独立散射系数建立部分稳定性的结果。
• 构建一个几何框架，其中阿尔贝多算子将一个平面上的入射条件映射到另一个平面上的出射测量值，从而简化分析。
• 在次临界性条件下，建立输运算子在 $L^1$ 空间中的适定性和有界性。
• 为在实际成像几何中，从边界测量值稳定重建反向散射和吸收系数提供基础。

提出的方法

• 通过径向延拓变换域，将原始边界 $\Gamma_\pm$ 映射到与速度向量正交的新平面 $F_\pm$，从而简化单次散射贡献的分析。
• 将阿尔贝多算子 $\mathcal{A}$ 定义为从 $L^1(F_-)$ 到 $L^1(F_+)$ 的映射，表示从入射数据到出射数据的正向输运解。
• 采用泛函分析方法将阿尔贝多算子分解为奇异分量，通过使用支撑缩小的测试函数分离出单次散射贡献。
• 基于算子 $\mathbf{T} = \mathbf{T}_1 + A_2$ 的扰动论证，其中 $\mathbf{T}_1$ 为自由输运算子，$A_2$ 表示散射，证明其有界可逆性。
• 应用弗雷德霍姆选择性与诺伊曼级数，证明 $\mathbf{T}^{-1}$ 在 $L^1(O, |v|\,dxdv)$ 中存在且有界，从而确保输运方程的适定性。
• 通过 $\mathcal{A}$ 在 $\mathcal{L}(L^1(F_-), L^1(F_+))$ 中的算子范数有界，建立系数差值的稳定性。

实验结果

研究问题

• RQ1 能否为静止线性输运方程中空间依赖散射和吸收系数的重建建立 $L^1$-稳定性估计？
• RQ2 测量设置的几何结构如何影响反向输运重建的稳定性？
• RQ3 单次散射贡献在稳定性分析中起什么作用，如何通过使用支撑缩小的测试函数将其分离？
• RQ4 在何种条件下，输运算子 $\mathbf{T}$ 在基于 $L^1$ 的函数空间中可逆？
• RQ5 是否可以在一种变换后的几何中分析阿尔贝多算子，使得入射和出射数据定义在与速度向量正交的平面上？

主要发现

• 该论文在空间依赖散射和吸收系数的情况下，建立了 $L^1$-稳定性估计，推广了以往仅限于空间独立散射系数的成果。
• 在系数的适当次临界性假设下，阿尔贝多算子 $\mathcal{A}$ 被证明是从 $L^1(F_-)$ 到 $L^1(F_+)$ 的有界算子。
• 输运方程的解在 $L^1(O, |v|\,dxdv)$ 中唯一可解，且解算子有界，其范数受 $\|\sigma\|_{\infty}$ 和 $\|\sigma_p\|_{\infty}$ 控制。
• 输运算子 $\mathbf{T}$ 的逆在 $L^1(O, |v|\,dxdv)$ 中有界，这确保了正向问题的稳定性，并为推导稳定性估计提供了可能。
• 稳定性估计通过 $\mathbf{T}^{-1} = \mathbf{T}_1^{-1}(I + A_2\mathbf{T}_1^{-1})^{-1}$ 的诺伊曼级数展开获得，其中算子 $A_2\mathbf{T}_1^{-1}$ 的范数有界于 $1 - e^{-2R\mathrm{v}_0^{-1}\|\sigma_p\|_{\infty}} < 1$，从而保证了级数收敛。
• 最终的稳定性界表示为 $\|f\|_{L^1(O,|v|dxdv)} \leq C_0 \|f_-–\|_{F_-}$，其中 $C_0$ 依赖于 $R$、$\mathrm{v}_0^{-1}$、$\|\sigma\|_{\infty}$ 以及扰动算子逆的范数。