体外冲击波碎石焦点区域声压分布时域有限差分法.pdfVIP

体外冲击波碎石焦点区域声压分布的时域有限差分法 仿真研究 1 1 1 2 2 菅喜岐 , 李文龙 , 刘东顺 , 史启铎 , 谭中 （1.天津医科大学生物医学工程系，天津 300070；2.天津医科大学第二医院，天津 300211） 摘要：目的 研究体外冲击波碎石(Extracorporeal shock wave lithotropisy ，ESWL)治疗中不同厚度人体组织对 ESWL 焦点区 域大小、位置和声压分布的影响，为 ESWL 治疗计划的制定提供理论依据。方法 以 Reichenberger 的 ESWL 水中实验为参 照例，分别建立三维单一人体组织和多层复合人体组织仿真模型，用时域有限差分（Finite Difference Time Domain ，FDTD ） 数值仿真法，数值仿真 ESWL 治疗中形成的焦点区域。结果 高强度聚焦超声波在人体脂肪、肌肉层中形成的实际焦点区域 大小、位置和焦点区域声压分布，随着人体脂肪、肌肉层厚度的变化而变化。结论 ESWL 治疗过程中由于人体组织厚度的 不同，高强度聚焦超声非线性传播而形成的焦距和最大声压不同，对于不同体态的患者在 ESWL 治疗时，应该采用不同的 焦距设定方法。 关键词：体外冲击波碎石；FDTD；脉冲超声波的非线性传播 中图分类号：Q681 文献标识码：A 引言 目前，4mm 以上的泌尿外科结石临床治疗的首选治疗方法是ESWL ，被应用于临床的ESWL 装置焦 点区域最高声压一般设定在 40~120MPa 之间。用现有 ESWL 装置对结石症进行治疗时，目前临床上只 能根据临床经验来选择发生器发出的超声波或冲击波强度的按键,不能根据患者的体形、结石大小和种 类，相应调整控制发生器发出的能量、聚焦超声波形成的焦区体积和最高声压等；在用现有实际临床治 疗时，即使定位系统的焦点位置与人体内结石位置一致设定，有时可能碎石效果不佳，有时也可能会对 人体组织造成伤害，如出现尿血、肾血肿等并发症，重者出现肾破裂等[1,2] 。为了提高治疗效果、减少对 人体组织的伤害，有必要了解、掌握 ESWL 治疗过程中，实际形成的焦点位置、大小、焦点处的声压分 布等。通过实验测定法，存在现有高强度冲击波感知器价格非常昂贵，人体内不易检测等问题。有必要 用理论或者数值解析的方法分析焦区声压的分布，目前国内外做的相关研究中，我国西安交通大学的陆 明珠等[3]利用积分算法的超声传播公式结合遗传算法对 256 个阵元的内球壳式相控阵聚焦换能器的声场 进行了仿真研究，并根据不同的换能器阵元划分得到了单焦点、多焦点等不同的聚焦方式。但是其研究 采用的是积分算法无法对高强度超声波传播过程的非线性特性对声场的影响进行计算。重庆大学的陈景 秋等[4]对三种不同参数的椭圆球面聚焦形成的焦点位置进行了仿真研究，该研究由于采用了 Tait 状态方 程式进行计算，其仿真计算量非常大，一般采用巨型计算机才可进行仿真计算。美国华盛顿大学的 Ibrahim M. Hallaj 和波士顿大学的 Robin O. Cleveland[5]利用 KZK 方程、Pennes 生物热传导方程和 FDTD 仿真法 比较了不同声压时非线性效应对焦点处产生的温度的影响，在该研究中由于使用KZK 方程进行仿真， *基金项目：教育部科学研究重点项目（205006 ）; 教育部留学回国人员科研启动基金 作者简介：菅喜岐，（1962-），男，工学博士，副教授，主要研究方向：超声医学 通讯作者：菅喜岐，电话：022E-mail ：jianxiqi@ 1 需要对二阶微分进行仿真计算，其计算复杂计算量大。笔者等利用声波的运动方程和连续方程和 1966 年由 Yee 提出应用于电磁波模拟计算的 FDTD 方法，提出了在体积弹性系数内包含超声非线性量的算法 [6-10