经颞窗动态聚焦超声的校准仿真

《Cell Reports Methods》：Calibrated simulations for dynamic focusing of ultrasound through the temporal window

　　

聚焦超声可通过颞窗送达以调节异质分布的脑区。声学仿真（acoustic simulations）可在动态靶向脑结构时用于安全性评估，但由于对颅骨与换能器做了机械上不准确的假设，仿真与实测焦压之间的失配会随可转向范围变化。研究人员描述了利用轴对称投影与跨157名受试者深部脑靶点的三维可转向子空间中的稀疏采样进行仿真—测量校准的高效方法。为处理颅骨透射中仿真—现实失配，研究人员利用经八块人颞窗骨片测得的与预测的压强值推导出优化的骨衰减系数。综上，该校准框架与优化颞窗系数可广泛跨研究使用，以提升个性化神经调控治疗中报告准确性及依赖的安全性评估。

研究背景与问题提出

超声神经调控（ultrasound neuromodulation）作为新兴脑刺激技术，已在难治性抑郁、精神分裂症纹状环路干预等场景显现治疗潜力。多元件阵列可在二维平面驱动下于三维空间动态转向焦点，但仿真预测的焦压常与实测显著偏离，且偏差随靶点位置变化。现有校准方法要么依赖自由场单点标定、要么需用声学全息（acoustic holography）逐元反演，难以在大规模个性化部署中复用，且对转向空间位置“盲校准”。更关键的是，颅骨（尤其是颞窗骨小梁结构）的微观散射无法被常规CT/MRI分辨率捕获，导致以吸收为主的均匀衰减项低估了真实透射损失。为支持个性化神经调控的安全规划，亟需一套兼顾设备级与颅骨级失配、并能给出最坏情形误差界的实用校准框架。该研究发表于《Cell Reports Methods》。

主要关键技术方法

研究人员从157名受试者MRI中提取尾状核、杏仁核、丘脑中央中核、丘脑底核四类双侧深脑靶点相对单阵列的坐标，将圆对称覆盖展平为轴向—横向校准平面，定义覆盖95%靶点的截锥校准空间。在水中以6 mm等距稀疏采样（二维平面约50点、三维锥约20点）用水听器测峰压，与k-Wave/BabelBrain前向仿真峰压作差，用各向异性高斯函数拟合低频失配面；另取5具尸头以硅胶垫模拟头皮，定位佩戴式器件后切取左右颞窗骨片（8片有效），行零回波时间（ZTE）MRI提取虚拟颞窗，以逐步均匀赋值骨衰减系数α₀做前向仿真，对照骨片后向时间反转聚焦的实测峰压，用均方误差最小化搜优α₀。

研究结果

引言部分指出，仿真—测量失配分设备相关（元串扰、指向性）与患者相关（颅骨微结构）两类，设备级表征是个性化建模前提；多元件动态转向若无差分校准将产生大幅误差。

结果部分首先展示：基于157受试者靶坐标定义的校准空间内，声学全息反演元源后在可转向范围峰压失配超20%，主因是隔离测量未捕获大串扰；改用焦压图差异的各向异性高斯拟合，二维平面随机子采样至50点（6 mm间距）即达R²=0.83，三维全空间6 mm采样下约20点即收敛，贝叶斯模型选择确认三维各向异性高斯最优，且在另5只同构阵列复现。该焦校法在500 kHz稀疏阵列与900 kHz费马螺旋碗形换能器均适用。

继而研究人员用8片颞窗骨片做透射测量：时间反转靶向40点/片，平均靶偏1.68±0.96 mm；扫描均匀骨衰减α₀发现MSE最小点在15.3 dB/cm/MHz，与文献13.3–16.6 dB/cm（1 MHz）吻合；以此优化值计，仿真焦压的5th百分位高估33%、95th百分位低估44%。热例（5 Hz、占空比10%、30 min，吸收取衰减16.8%）显示骨温升界+2°C内可被该误差带覆盖。

讨论与结论翻译

研究人员提出一套跨转向空间校准阵列焦点仿真的方法，并优化颅骨衰减以提升峰值压误差界准确性，从而改进神经调控报告与安全参数理解。无转向或轴对换能器可简化校准；颅骨衰减以集总参数给出，其值依附于所用仿真框架（网格、幂律指数、黏弹求解器），换框架需重标。校准空间定义亦可借人群生物计量或大模型MRI数据集而定。焦校法仅对齐峰压失配，旁瓣/栅瓣失配未校；8片尸颞窗未覆盖活体头皮灌注与 demographic 微结构变异，故5th–95th误差带为保守安全上限。综上结论：稀疏采样校准优于声学全息；优化颞窗衰减系数α₀=15.3 dB/cm/MHz；误差界支持个性化治疗最坏情形评估；该框架须作患者特异性颅骨传输建模的前置基线而非替代。