Linearised Calderón problem: Reconstruction of unbounded perturbations in 3D

这是一份关于论文《线性化 Calderón 问题：三维无界扰动的重构》（Linearised Calderón Problem: Reconstruction of Unbounded Perturbations in 3D）的详细技术总结。

1. 研究背景与问题定义

核心问题：
本文研究的是线性化 Calderón 问题（Linearised Calderón Problem）在三维空间中的重构问题。

• 物理模型：考虑单位球 $B \subset \mathbb{R}^3$ 中的电导率问题。电导率系数为 $\gamma = 1 + \eta$，其中 $\eta$ 是待重构的扰动项。
• 线性化设定：假设背景电导率为常数 1，扰动 $\eta$ 较小。定义 Fréchet 导数算子 $F = D\Lambda(1; \eta)$，其中 $\Lambda(\gamma)$ 是诺伊曼 - 狄利克雷（Neumann-to-Dirichlet, ND）映射。
• 重构目标：从线性化数据 $F\eta$（即施加电流密度 $f$ 和 $g$ 时边界电势变化的线性响应）中精确重构扰动 $\eta$。
• 关键挑战：
  • 传统的重构方法通常假设扰动 $\eta$ 属于 $L^\infty$ 或 $L^2$ 空间。
  • 本文旨在处理无界扰动，即允许 $\eta \in L^3(B)$（在三维中，$L^d$ 空间是线性化算子连续扩展的自然空间）。
  • 需要解决从边界测量数据到内部扰动的直接、非迭代重构算法。

2. 方法论：基于 3D Zernike 基的前向代入法

本文提出了一种基于**三维 Zernike 基（3D Zernike Basis）**的直接重构算法，其核心思想是“前向代入”（Forward Substitution），类似于二维情况下的方法，但针对三维几何进行了推广。

2.1 数学工具

• 3D Zernike 基函数：
  定义在单位球 $B$ 上的正交基函数 $\psi_{\ell, k}^m(r, \theta, \phi) = R_\ell^k(r) Y_\ell^m(\theta, \phi)$。
  • $Y_\ell^m$：球谐函数（Spherical Harmonics），处理角向部分。
  • $R_\ell^k$：3D 径向 Zernike 多项式，处理径向部分。
  • 该基构成了 $L^2(B)$ 的正交基，且可以展开 $L^3(B)$ 中的扰动 $\eta$。
• 线性化算子的性质：
  利用 Fréchet 导数的积分形式：
  $$\langle (F\eta)f, g \rangle_{L^2(\partial B)} = -\int_B \eta \nabla u_f \cdot \nabla u_g \, dx$$
  其中 $u_f, u_g$ 是对应边界条件 $f, g$ 的调和函数。

2.2 核心算法逻辑

算法利用球谐函数作为边界电流密度进行测量，通过特定的三角结构实现系数的递推求解：

1. 扰动展开：将未知扰动 $\eta$ 展开为 Zernike 基的级数：
   $$\eta = \sum_{k, \ell, m} c_{\ell, k}^m \psi_{\ell, k}^m$$
   目标是求解系数 $c_{\ell, k}^m$。

2. 特定测量策略：
   为了提取特定系数，选择特定的边界电流密度。对于给定的径向阶数 $k$，施加的电流密度为 $Y_0^{k+1}$（仅涉及 $m=0$ 的球谐函数），接收端测量 $Y_m^{\ell+k+1}$。

   • 优势：在三维中，不需要像二维那样使用所有 $m$ 的傅里叶模式，仅需 $m=0$ 的特定模式即可，大大减少了数据冗余。

3. 递推公式（前向代入）：
   利用球谐函数乘积的有限展开性质（Gaunt 系数）和径向多项式的正交性，推导出系数 $c_{\ell, k}^m$ 的显式公式：
   $$c_{\ell, k}^m = (Q_{\ell, 0}^{k, m})^{-1} \left( \langle (F\eta)Y_0^{k+1}, Y_m^{\ell+k+1} \rangle - \sum_{q=0}^{k-1} \sum_{s=0}^{k-q} Q_{\ell, s}^{k, m, q} c_{\ell+2s, q}^m \right)$$

   • 三角结构：计算 $k$ 阶的系数时，仅依赖于 $k$ 阶的测量数据以及更低阶（$q < k$）已计算出的系数。
   • Gaunt 系数：公式中涉及球谐函数三重积分的 Gaunt 系数，这些系数是预先计算好的常数，与 $\eta$ 无关。

4. 重构步骤：

   • $k=0$：直接由数据 $(F\eta)Y_0^1$ 重构 $\eta_0$（调和函数部分）。
   • $k=1$：利用数据 $(F\eta)Y_0^2$ 和已知的 $\eta_0$ 重构 $\eta_1$。
   • 归纳：依次递增 $k$，利用 $(F\eta)Y_0^{k+1}$ 和之前的 $\eta_0, \dots, \eta_{k-1}$ 重构 $\eta_k$。

3. 主要贡献与结果

3.1 理论贡献

• 扩展到 $L^3$ 空间：证明了线性化算子 $F$ 可以连续扩展到 $L^3(B)$ 空间，从而允许重构无界扰动（Unbounded Perturbations）。这在理论上放宽了对扰动光滑性的要求。
• 精确重构定理：证明了对于任意 $\eta \in L^3(B)$，可以通过上述前向代入法精确重构其 Zernike 展开系数（定理 1.1）。
• 数据效率：指出在三维情况下，重构所需的边界测量数据量比全基测量要少。由于球谐函数乘积的性质，仅需 $m=0$ 的特定电流密度即可，避免了 $m$ 索引的冗余（$2\ell+1$ 增长）。

3.2 数值实验结果

• 理想情况：在 Mathematica 中模拟高精度测量数据（16 位有效数字），算法能够近乎完美地重构局部化的高斯型扰动，即使对径向阶数 $\ell$ 进行截断。
• 实际情况：使用粗糙的有限元（FE）网格（P1 单元，约 6000 节点）模拟带有噪声的测量数据。
  • 结果显示，尽管存在数值误差，重构结果依然保持了良好的稳定性。
  • 随着截断阶数 $K$ 的增加，重构图像逐渐逼近真实扰动。
  • 证明了该方法对测量噪声具有一定的鲁棒性，且计算出的系数在数值上是稳定的。

3.3 与二维方法的对比

• 相似性：都利用 Zernike 基和三角结构实现直接重构，无需迭代优化。
• 差异性：
  • 数据需求：三维方法利用 $m=0$ 的特定模式，数据冗余更少；二维方法需使用所有傅里叶模式。
  • 数值稳定性：三维方法由于球谐函数乘积展开涉及更高阶 $\ell$ 的项（求和项），理论上比二维方法多一层求和，可能数值稳定性稍差，但实验表明仍可获得良好近似。

4. 意义与展望

• 计算效率：该方法属于“前向代入”，计算复杂度低，非常适合数值实现，避免了传统迭代反演方法（如牛顿法）的高计算成本和局部极小值问题。
• 处理无界扰动：能够处理 $L^3$ 空间的扰动，意味着该方法可以应用于电导率分布中存在较大突变或奇异性的场景，扩展了 Calderón 问题的应用范围。
• 实际应用潜力：数值实验表明，即使在粗糙网格和含噪数据下，该方法仍能产生合理的重构结果。这为医学成像（如 EIT）和地球物理勘探中的三维电阻率成像提供了新的理论工具和算法基础。
• 未来方向：虽然目前针对的是单位球，但该方法的核心思想（利用正交基和线性化算子的三角结构）可能为其他几何形状或非线性问题的重构提供思路。

总结：
这篇论文成功地将二维线性化 Calderón 问题的直接重构算法推广到了三维，并引入了 3D Zernike 基。其核心创新在于利用特定的测量策略和球谐函数的代数性质，构建了一个高效的、非迭代的、能够处理无界扰动的直接重构公式。数值实验验证了该方法在理论上的精确性和在实际含噪环境下的可行性。