基于K-wave的经颅超声聚焦仿真研究

作者:杨施成; 刘逍逸; 王欢; 刘春泽; 陈友元; 周红生
来源:声学技术, 2023, 42(01): 39-45.
DOI:10.16300/j.cnki.1000-3630.2023.01.007

摘要

由于颅骨形状不规则,声速、密度等声学参数分布不均匀,使用传统相控聚焦的方法进行经颅超声精确聚焦时会出现焦点偏移、焦点形状畸变、散焦的现象,时间反转法被认为是能够克服颅骨非均匀特性,实现经颅聚焦的有效手段。而利用传统有限元,有限差分法进行数值仿真时,往往需要划分十分精密的网格进行计算,耗费大量的时间。K空间伪谱法通过在空间域上进行傅里叶变换,时间域上进行有限差分的方法求解声波方程,能够在降低计算网格密度的同时保证计算精度。为了克服颅骨对聚焦超声造成的焦移影响,同时避免声学模型计算量过大的问题,本文使用Matlab K-wave开源工具箱,对食蟹猴颅骨CT数据进行三维声学建模,实现了基于时间反转法的经颅超声聚焦,讨论了时间反转法对于颅骨不规则形状及非均匀声速造成的焦移的补偿。仿真结果表明,K-wave在提高计算速度的同时能够服务于精准经颅超声聚焦,且相较于传统基于声时差法的相控聚焦方法,时间反转法在预设焦点位置,焦点强度更高,焦斑形状更为规则,聚焦效果更好。

全文