通过分布式源反演重建有效超声换能器模型

《Ultrasonics》：Reconstructing effective ultrasound transducer models via distributed source inversion

　　

对超声波传播的精确建模对于超声波检测（Ultrasonic Testing，UT）中的高保真仿真和成像至关重要。一个主要的挑战在于表征激励源，特别是对于孔径相对于声波长较大的换能器。在这种情况下，非均匀激励和空间干涉显著影响产生的辐射方向图。本文提出了一种分布式源反演（Distributed Source Inversion，DSI）策略，以重建一个有效的时空换能器模型，该模型能够重现实验测量的波场。重建的源模型捕获了孔径相关的相位和振幅变化，而不需要详细了解换能器结构。该方法通过铝半圆柱上的方向性测量进行了验证，其中结合重建源模型的仿真与实验方向性图和波形形状显示出良好的一致性。最后，对逆时偏移（Reverse Time Migration，RTM）和全波形反演（Full-Waveform Inversion，FWI）的合成研究表明，精确的换能器建模对于基于仿真的成像和反演工作流的成功至关重要，并显著提高了重建质量。

**论文解读：分布源反演重建有效超声换能器模型及其在仿真成像中的关键作用**

**研究背景与问题**

在超声波检测（Ultrasonic Testing，UT）中，高保真仿真与成像高度依赖对波传播的精确建模。其中，激励源的表征是核心挑战，尤其当换能器孔径与声波长相当甚至更大时，非均匀激励和空间干涉会显著改变辐射方向图（radiation pattern）。传统射线类方法（如合成孔径聚焦技术SAFT、全聚焦法TFM）基于直线射线近似，忽略源特性；而基于仿真的成像方法，如逆时偏移（Reverse Time Migration，RTM）和全波形反演（Full-Waveform Inversion，FWI），虽能利用波场的完整信息，但要求换能器模型高度匹配真实物理。然而，实验条件下换能器的具体结构、压电晶振的固有模态以及耦合条件常未知，导致源特性估计不准确，严重制约仿真反演的收敛性和可靠性。现有源估计方法（如从头波提取、反卷积、水声标定等）或假设空间均匀性，或仅恢复时间波形，难以同时捕获孔径上相位和振幅的时空变化。因此，需发展一种无需详细换能器内部知识、能直接从实测数据重建有效时空源模型的方法。

**研究内容与意义**

为应对上述挑战，研究人员提出了一种基于伴随的分布式源反演（Distributed Source Inversion，DSI）框架，旨在从实验测量数据中重建超声换能器的有效时空模型。该方法将换能器建模为分布在孔径上的点源集合，通过迭代优化使模拟波场与观测波场的L2范数最小化，从而同时恢复空间和时间的源特征。研究利用铝半圆柱的方向性测量数据集验证了DSI的有效性，并进一步通过合成RTM和FWI实验评估精确源建模对成像与反演质量的影响。该研究发表于《Ultrasonics》，其核心意义在于：为仿真成像方法（尤其是FWI）提供了获取精确换能器模型的实用途径，显著提升了基于仿真的缺陷检测的轴向和横向分辨率，并揭示了源模型不准确是限制仿真反演成功的关键瓶颈。

**主要关键技术方法**

研究人员采用了以下关键技术方法：1）基于伴随方法的梯度计算：通过一次正向模拟和一次伴随（adjoint）模拟，高效计算目标泛函对源系数的梯度，伴随力由数据残差加权得到；2）有限记忆Broyden–Fletcher–Goldfarb–Shanno（L-BFGS）优化算法：迭代更新源系数，结合最优步长计算，实现快速收敛；3）谱元法（Spectral Element Method，SEM）作为正演求解器，采用Gauss–Lobatto–Legendre（GLL）点进行空间离散和中心差分法（CDM）进行时间积分；4）源参数化：将有效源表示为空间基函数（点源位置）和时间基函数（线性插值）的乘积组合，系数为待求优化变量。实验数据来源于慕尼黑工业大学文献[21]的铝半圆柱数据集，换能器为A V103-RB型P波换能器（直径12.7mm，标称频率1MHz），接收器为17个等间距分布在曲面上的相同换能器。

**研究结果**

**3. 源校准（Source calibration）**
**3.1 实验设置与模拟配置**：利用铝半圆柱（直径0.3m）上的方向性测量数据，将分布式源建模为20个沿x∈[-0.01m,0.01m]均匀分布的点源，每个点源独立波形。比较两种接收加权策略：均匀加权（DSI 1）和振幅补偿加权（DSI 2，权重为1/A_r，其中A_r为接收信号最大振幅）。
**3.2 分布式源反演**：从零源模型开始，依次进行10次最速下降和40次L-BFGS迭代。结果表明：DSI 2在较高孔径角上显著改善了模拟与实验波形之间的相关系数（Correlation Coefficient，CC）和归一化均方误差（Mean Squared Error，MSE），而DSI 1仅在小角度准确。重建的源振幅分布呈现轻微右偏，但仍在名义孔径12.7mm内，且有效捕获了圆形孔径的三维效应。DSI 2重建的波形在各角度与实验数据吻合良好，而点源模型无法产生正确的方向性。
**3.3 源采样（Source sampling）**：改变点源数量Ns=5,10,20,40,80,160。发现当Ns≥10（即每P波波长至少4个源点）时，重建质量几乎不受影响；Ns=5时无法恢复有效方向性。
**3.4 接收采样（Receiver sampling）**：改变接收器数量Nr=3,5,9,17。结果表明，接收器间距≤20°（Nr≥9）可成功重建整个方向性，而稀疏配置（80°或40°间距）无法内插中间角度的辐射特征。

**4. 对RTM和FWI的影响**
**4.1 研究设置**：以DSI 2为真实源模型生成合成参考数据（含三个不同直径的孔缺陷），分别用点源（PS）、均匀源（US）、DSI 1和DSI 2四种源模型执行RTM和FWI。
**4.2 结果**：RTM中，PS无法识别任何缺陷；US、DSI 1和DSI 2能检测上部缺陷，但均不能分辨底部孔，且异质源模型（DSI）伪影较少。FWI中，PS和US失效；DSI 1与US类似，重建质量差；只有DSI 2成功重建了全部三个缺陷，且F1分数、AUPRC和AUROC显著优于其他模型，表明精确的高孔径源建模是FWI利用长偏移数据提高分辨率的关键。

**总结讨论与结论**

讨论部分指出：DSI成功重建了有效换能器模型，对空间离散化鲁棒，源采样的最小要求为每P波波长4个点；接收器间距20°可捕获完整方向性；振幅补偿加权（1/A_r）提供了平衡的重建效果。仿真成像研究表明：RTM对源模型敏感性较低，但FWI高度依赖源准确性，只有DSI 2允许一致利用低和高孔径数据，从而实现高质量损伤识别。计算成本方面，2D情况下单次SEM模拟约5秒（NVIDIA GeForce RTX 3090 GPU），完整源校准<15分钟；3D扩展需并行策略。局限性包括：忽略三维几何扩散和衰减、源模型转移性和接收器不确定性需进一步研究。

结论部分（翻译自“6. Conclusions and outlook”）：
本研究引入了一种基于伴随的分布式源反演（DSI）框架，用于直接从超声测量中重建有效换能器模型。该方法避免了压电元件的显式建模，恢复了一个空间分布的、随时间变化的源，该源隐式捕捉了实验不确定性，如压电缺陷、耦合条件和定位误差。采用振幅补偿接收加权策略重建的源模型，在整个孔径角度范围内准确重现了振幅和波形特征。此外，该方法对空间源离散化和接收采样密度表现出强鲁棒性。基于这些结果，使用RTM和FWI的合成研究凸显了精确源建模对基于仿真的成像的重要性。RTM相对较不敏感，因为它主要依赖主导波场分量；相反，FWI通过利用长偏移数据，提高了轴向和横向分辨率，尤其是对于靠近背壁的深层缺陷。在这种场景下，精确表征源特性（尤其在高孔径角）对于充分利用可用波形信息并实现高质量重建至关重要。尽管存在若干限制（如三维效应、转移性和接收器不确定性），DSI在基于仿真的成像之外还有显著潜力，例如生成与物理高度一致的合成数据，用于验证常规方法（如SAFT、TFM）和训练数据驱动方法，从而为实验测量、数值建模和射线或数据驱动方法之间提供实用桥梁，支持更鲁棒、高分辨率的检测技术发展。