基于频域反向传播(Frequency-Domain Backpropagation)重建光学不可达环境中微气泡振荡动力学的重构方法

《Results in Engineering》：Reconstruction of microbubble oscillation dynamics with frequency-domain backpropagation in optically inaccessible environments

　　

摘要：日益增多的超声介导诊断与治疗技术（包括声孔效应/Sonoporation及血脑屏障/Blood-Brain Barrier, BBB调控）依赖于微气泡(Microbubble, MB)造影剂，其中微气泡动力学行为的精确控制决定了生物响应、效率与安全性。然而，在光学不透明及深部组织环境中，对处于稳定弱非线性(Weakly Nonlinear) regime的微气泡振荡进行定量监测仍然具有挑战性。本研究提出一种线性声波传播与叠加(Linear Acoustic Wave Propagation and Superposition, LAWPS)框架，可直接从被动记录的声发射(Acoustic Emission, AE)中重构微气泡半径-时间(Radius-Time)动力学。通过将弱非线性振荡的傅里叶级数(Fourier-Series)表示与线性单极子辐射理论(Monopole Radiation Theory)耦合，LAWPS将经典单极子模型扩展至建立多频声发射与潜在径向气泡动力学之间的可逆关系。将该框架推广至单分散(Monodisperse)微气泡群，研究人员推导了最优激励与接收配置及群体水平修正因子，使得能够定量重构群体平均微气泡动力学。利用同步光学与声学测量，研究人员证明在振荡幅度可达平衡半径R0约15%时，微气泡振荡的恢复相对误差约为5%。最后，研究人员表明在该框架工作区间（ΔR/R0≤ 20%）内的振荡可在生理相关条件下对小至10 μm的囊泡产生与声孔效应相关的机械应力（毛细管数CaK≥ 0.01）。综上，本工作建立了一种基于声发射监测光学不可达环境中弱非线性微气泡振荡的定量框架，可提升新兴生物医学超声应用及其他领域的监测与控制能力。

论文解读：基于频域反向传播重构光学不可达环境中微气泡振荡动力学

超声介导的诊断与治疗技术（如声孔效应Sonoporation、血脑屏障BBB调控及超声造影成像）广泛依赖微气泡MB造影剂。现有应用多利用惯性空化与MB破裂产生宽频发射来指示空化活动，但越来越多的治疗应用（如细胞追踪、BBB表型调控）依赖稳定非惯性（即稳定非线性）MB振荡，其精确监测与控制决定安全性与疗效。目前定量表征MB振荡半径-时间(R(t))主要依靠高速光学显微术，但该方法局限于体外(in vitro)或浅表透光组织，无法用于深部或临床相关的光学不透明环境。被动空化检测(Passive Cavitation Detection, PCD)虽可记录MB声发射AE，但传统频域分析仅提供定性信息（谐波/次谐波含量），无法定量还原R(t)。少数专用声学方法（如"声学相机"需>10 MHz受限衰减、时域反演法对噪声敏感且未拓展至多泡群）存在明显局限。因此，亟需一种能从被动记录的AE中定量重构弱非线性MB振荡动力学、适用于光学不可达环境的方法。

为此，Georgia Institute of Technology的研究人员Lee H, Zangabad RP, Kim C, Menezes V, Park J, Degertekin FL, Arvanitis C提出了线性声波传播与叠加(LAWPS)框架，通过傅里叶级数分解将经典线性单极子辐射理论推广至多频弱非线性振荡，建立AE各次谐波复傅里叶系数与MB表面速度复傅里叶系数的对应关系，实现从记录压力反推表面速度再积分得R(t)（频域反向传播），并推广至单分散MB群并引入破坏性干涉/建设性干涉(DI/CI)比校正。实验利用同步高速显微光学与PCD验证。结果表明LAWPS在ΔR/R₀< 20%弱非线性区可高精度重构R(t)，且该区间振荡可产生生理相关声孔效应级机械应力。论文发表于《Results in Engineering》。

主要关键技术方法：

研究人员采用理论推导（将MB振荡R(t)及远场辐射压力P(r,t)作傅里叶级数展开，基于单极子条件kR₀?1推导各次谐波AE与表面速度的频域传递关系，建立LAWPS正向传播与反向传播公式）；数值仿真验证（使用Hoff脂质壳MB模型及Marmottant模型生成R(t)与AE，对比Vokurka方法验证LAWPS精度；采用泊松圆盘采样Poisson Disk Sampling生成随机分布MB群，计算泡间次级声压及DI/CI比统计分布以推导群平均修正因子）；体外同步声-光实验验证（微流控ibidi通道中注入单分散Bracco BG9732 MBs（R₀≈3 μm，稀释至空泡分数Void Fraction β=0.015%），以0.5 MHz超声激励，同步采集5 Mfps高速显微图像与经校准PCD记录的AE，扣除背景后施以DI/CI比、MB数目、通道透反射损失校正，进行LAWPS反向传播并与光学测得R(t)比对）；理论力学评估（基于Rayleigh-Nyborg-Westervelt微流RNW Streaming及毛细管数Ca_K评估LAWPS工作区内振荡产生的膜剪切应力与声孔效应相关性）。

研究结果：

2.1. Derivation of relationship between single MB oscillation and AE（单MB振荡与AE关系的推导）： 假设MB作纯球形弱非线性振荡（kR₀?1, ΔR?λ），将R(t)与远场P(r,t)展开为傅里叶级数，结合线性单极子波动解导出第n次谐波AE复系数P_n与表面速度复系数?_n的关系式（含距离r、介质密度ρ₀、角频率ω₀、平衡半径R₀及arctan相位项）。由此得到LAWPS正向（由?(t)算P(r,t)）与反向（由P(r,t) FFT各阶P_n反算?_n再求和得?(t)，积分得R(t)）算法，框架与具体MB本构模型无关。

2.2. Validation of LAWPS propagation algorithm with single MB model（单MB模型验证LAWPS正向算法）： 用Hoff模型生成R₀=2~4 μm、f₀=0.5 MHz、激励10~200 kPa的MB振荡及AE，对比LAWPS正向传播与Vokurka经典AE公式。结果显示两者波形吻合（含声速传播延迟），n=7阶谐波近似时AE峰峰值误差约5%；高阶谐波项贡献随n增大衰减，大R₀收敛更快，验证了LAWPS正向公式正确性。

2.3. Validation of LAWPS backpropagation algorithm（LAWPS反向传播算法验证）： 对Vokurka模拟AE施加LAWPS反向传播，发现仅需至2阶谐波（n=2）即可重构原始R(t)，峰峰值误差<3%；高于3阶谐波压力幅对R(t)重构贡献可忽略（传递函数中√((nk₀R₀)²+1)/n ∝ 1/n）。当ΔR/R₀>30%（接近壳破裂、出现惯性响应）时傅里叶级数失效误差增大；若输入错误R₀（偏差30%）反向传播误差可达15%，说明需较准确的R₀先验知识。换用Marmottant模型得相同结论，证LAWPS不依赖特定MB模型。

2.4. Multiple MB scheme and phase-based correction（多MB方案与基于相位的校正）： 模拟随机分布单分散MB群（泊松圆盘采样，β=0.05~0.5%），在低浓度（β≤0.26%，泡间距>14R₀）下相邻泡产生的次级声压<激发压4%，对原振荡影响<2%，泡间动力学耦合可忽略。空间分布引起ae相位相消（破坏性干涉di），研究人员定义di ci比=>

2.5. Experimental validation of LAWPS algorithm using concurrent acoustic and optical monitoring of MB dynamics（同步声光监测实验验证LAWPS）： 微流控实验中用BG9732 MBs（平均直径6 μm，多分散指数7%），稀释至β=0.015%，光学确认泡间距>15R₀。同步采集0.5 MHz激励下AE与5 Mfps光学影像，AE经背景减除、DI/CI比校正（基频4.67%）、通道损失校正、MB数目归一化后做5阶LAWPS反向传播（取R₀=3 μm±0.43 μm容差）。结果与光学测量吻合：ΔR/R₀≤20%时相对误差约5%（n=69个MB），能捕捉谐波振荡；超此区间误差增大。证实LAWPS结合DI/CI校正可有效从被动AE定量重构低浓度单分散MB群平均R(t)。

2.6. Sonoporation-relevant mechanical effects emerge within the LAWPS operating regime（LAWPS工作区内出现声孔效应相关机械效应）： 理论计算R₀=3 μm、f₀=0.5 MHz、ΔR/R₀=10%时，气泡流雷诺数Re_bs≈0.05（满足Stokes流），近壁最大切应变率?_max≈31.4×10³s^-1，介质η=水时Ca_K较小，但在血液等效高黏度介质(η≥6倍水)下对D=10 μm囊泡Ca_K≥0.01（可逆膜破裂临界值）。说明LAWPS适用区间（ΔR/R₀<20%）的稳定弱非线性振荡在生理黏滞条件下足以产生声孔效应相关机械应力。

讨论与结论：

讨论指出LAWPS通过频域谐波分解建立AE与MB表面速度可逆关系，类比泰勒展开捕获弱非线性主导项，模型无关且仅要求球形弱非线性振荡；相较既往时域反演法，频域反向传播抗噪且可通过选带滤波提升实用性，并首次拓展至多泡群DI/CI校正。源-接收器面对面放置最小化DI，为PCD系统设计提供新依据。局限性包括：假设稳态同相泡群（低浓度下次要）、需已知单分散R₀（可结合高频背向散射估算）、未完全建模边界近场非球振荡及体内衰减/血管分布（需3D几何相关DI/CI及超声定位显微ULM辅助）。LAWPS填补了光学不可达介质中AE定量监测MB R(t)的空白，可实时反馈控制MB稳定在所需振荡区，助力安全有效的MB-FUS治疗，且具硬件嵌入式实时处理潜力。

结论：研究人员建立的LAWPS框架可从被动记录声发射定量重构光学不可达环境中弱非线性微气泡振荡的半径-时间动力学，在ΔR/R₀≤20%下单泡重构误差~5%、经DI/CI校正可用于低浓度单分散微气泡群，且该工作区振荡在生理黏滞条件下可产生声孔效应相关机械应力，为超声空化物理、声学计量及生物系统超声调控提供了定量非侵入监测新工具。