体外冲击波碎石中组织厚度对焦点位置及其附近声压影响的时域有限差分数值解析.pdfVIP

体外冲击波碎石中组织厚度对焦点位置及其附近声 压影响的时域有限差分数值解析1 1 2 2 2 3 菅喜岐 ，王建民 ，史启铎 ，谭中 ，森田长吉 1 天津医科大学生物医学工程系，天津（300070 ） 2 天津医科大学第二医院，天津（300211 ） 3 千叶工业大学电气电子工学科，日本习志野（2750016 ） 摘 要：目的数值仿真体外冲击波碎石(Extracorporeal shock wave lithotropisy ，ESWL)治疗 中超声波非线性传播过程，研究不同厚度人体组织对ESWL 焦点区域声压及其实际形成焦 点位置的影响。方法 以实际ESWL 实验为例，建立仿真模型，用FDTD（Finite Difference Time Domain ）法进行数值仿真。结果 高强度聚焦超声波在人体脂肪、肌肉层非线性传播过程 中形成的实际焦点位置，随着人体脂肪、肌肉层厚度的不同而变化。结论 ESWL 治疗时， 随着患者体形的不同在人体内形成的ESWL 治疗焦点位置而移动。 关键词：ESWL ，脉冲超声波的非线性传播，FDTD 解析，实际焦点位置 目前，4mm 以上的泌尿外科结石临床治疗的首选疗法是 ESWL 。这种疗法是应用人 体外波源发出的低能量脉冲超声波或者冲击波通过水体传入人体内，在人体内聚焦形成高 能量、高强度的冲击波，用形成的这种冲击波破碎泌尿系统结石，结石碎片随着尿液排出 体外的疗法。ESWL 装置的波源主要有发低能量超声波的陶瓷压电式、电磁感应式和发低 能量冲击波的液电冲击式、微爆炸式发生器等。被应用于临床的 ESWL 装置的焦点区域 内最高声压一般设定在40-120MPa 之间，在 40-120MPa 声压之间设有 4 个声压高低不同 的超声波发射按键。 用现有ESWL装置对结石症进行治疗时，目前临床上只能根据临床经验来选择发生器 发出的超声波或者冲击波强度的按键,不能根据患者的体形和结石的大小和种类，相应调 整控制发生器发出的能量、聚焦超声波形成的焦区面积和最高声压等；在实际临床治疗时， 如果定位系统的焦点位置与人体内结石位置一致设定，有时临床碎石效果不佳。临床治疗 时只能根据临床经验设定最佳治疗焦点位置，这种凭经验的临床碎石治疗，可能会对人体 组织造成伤害，出现尿血、肾血肿等并发症，重者出现肾破裂等[1-4] 。为了提高治疗效果、 减少对人体组织的伤害，有必要掌握、了解ESWL治疗过程中，实际形成的焦点位置、焦 点处的声压分布等。通过实验测定及其数据分析的研究方法，存在现有高强度冲击波感知 器价格非常昂贵，几次实验测定之后感知器被破损，测量精度不能保证，人体内不易检测 等问题[1，5] 。有必要用理论或者数值解析的方法分析焦区声压的分布，但是到目前为止， 分析ESWL焦区的理论或者数值解析方法几乎没有[6-8] 。笔者等曾提出了超声波非线性传播 的FDTD仿真方法[9-12] ，并以Reichenberger 的电磁感应式ESWL水中实验为例[13，14],数值仿 真了焦点附近 4 个点的声压波形,得到与Reichenberger 的水中实验波形完全一致的仿真波 形[10] 。本研究仍以Reichenberger 的ESWL实验为例，研究高强度超声波传过不同厚度人体 肌肉组织时，不同的人体组织厚度对ESWL设备几何焦点附近点的声压、聚焦超声波形成 1本课题得到教育部科学技术研究重点项目基金（No.205006）和教育部留学回国人员科研启动基金的资助。 -1- 的实际焦点位置的影响。 1. 方法 1.1 基本方程式 当不考虑吸收衰减、声能密度（Lagrangian density ）和马赫(Mach)数远远小于时， [15,16] 声波非线性传播的方程式