CN105251140A - 一种聚焦声透镜的设计方法 - Google Patents
一种聚焦声透镜的设计方法 Download PDFInfo
- Publication number
- CN105251140A CN105251140A CN201510816714.4A CN201510816714A CN105251140A CN 105251140 A CN105251140 A CN 105251140A CN 201510816714 A CN201510816714 A CN 201510816714A CN 105251140 A CN105251140 A CN 105251140A
- Authority
- CN
- China
- Prior art keywords
- groove
- focusing
- depth
- width
- focus
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013461 design Methods 0.000 title abstract description 15
- 230000000694 effects Effects 0.000 claims description 30
- 230000005540 biological transmission Effects 0.000 claims description 17
- 230000002547 anomalous effect Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 230000000737 periodic effect Effects 0.000 claims description 7
- 230000005764 inhibitory process Effects 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 2
- 238000002604 ultrasonography Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000007423 decrease Effects 0.000 abstract description 2
- 238000009210 therapy by ultrasound Methods 0.000 abstract 1
- 238000012285 ultrasound imaging Methods 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 12
- 230000001225 therapeutic effect Effects 0.000 description 12
- 238000012360 testing method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000010606 normalization Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 206010028851 Necrosis Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010897 surface acoustic wave method Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- 239000003871 white petrolatum Substances 0.000 description 1
Landscapes
- Transducers For Ultrasonic Waves (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种聚焦声透镜的设计方法,属于聚焦超声领域。本发明通过在传统聚焦声透镜的声辐射面构建声学人工结构,即在传统聚焦声透镜的声辐射面设置周期性分布的同心环状凹槽结构,并采用有限元模拟计算对凹槽的结构参数进行优化设计,从而能够在特定频率范围内抑制声焦点处径向旁瓣的幅度,同时提高声能的聚焦增益。采用本发明的设计方法设计出的声透镜可以应用于聚焦超声成像和聚焦超声治疗领域,相比于现有的声透镜聚焦装置,能够有效降低径向旁瓣幅度,提高聚焦增益,减小焦点位置的偏移。
Description
技术领域
本发明属于超声治疗技术领域,更具体地说,涉及一种聚焦声透镜的设计方法。
背景技术
高强度聚焦超声(HIFU)技术是一种有效的无损创伤治疗技术,对体内深层、固体肿瘤组织的治疗具有很好的效果,其主要是通过一定的方式使超声波透过体表后聚集在病灶组织上,经过细胞与超声的相互作用,在病灶组织内产生不可逆转的凝固性坏死,从而达到治疗的目的。目前,该技术以其无害、安全、有效等优势得到了国内外许多学者和医生的广泛关注,并已在泌尿学、肿瘤学、神经外科、妇科、眼科等医学领域中被广泛应用和推广。
超声聚焦换能器是HIFU技术的核心器件,由于HIFU技术的应用对象是有生命的活体,治疗过程中要利用高强度超声能量消除病灶,但更重要的是要保证治疗对象的安全性,因此实现超声能量精准聚焦,即超声换能器的聚焦特性的好坏是HIFU技术治疗的关键。根据聚焦方式的不同,常见的超声聚焦换能器主要有球面自聚焦换能器和声透镜聚焦换能器。其中,声透镜聚焦是利用声波在弯曲界面上的折射来达到聚焦的目的,结构简单,成本低廉,能够通过改变透镜曲率来改变焦距,适合应用于超声成像和超声治疗。然而由于声透镜聚焦换能器的聚焦效果受到声波衍射效应的影响,焦点的旁瓣幅度相对较大,在超声成像中容易导致伪影和斑点等不利效果,在超声治疗中容易对目标治疗区域周围的人体组织造成一定的损伤。
Durnin在自由空间的波动方程中发现了非衍射解[Phys.Rev.Lett.58(15)1987],利用该理论实现的声束能够在轴向上保持稳定的幅度值。非衍射声束具有良好的声场特性,在医学成像,组织定征,多普勒流测量,材料无损检测等许多领域具有良好的应用前景,因而得到广泛研究。但这类解的算法复杂,能量激发效率不高,需要较大的空间尺度,从而限制了其推广应用。
近年来关于利用声学反常透射现象的声学人工结构的研究也取得了许多进展,在声学人工结构中利用声子晶体与基体界面上的表面波都可以获得非衍射声束(J.Appl.Phys.106(4),044512-044513,2009)。Christensen等人在声学刚性板上加上周期性凹槽结构后,利用在结构中激发的声表面波获得了声波反常透射和非衍射的效果[Nat.Phys.3(12),851-852,2007],该理论提供了一种高效简便的激发准直声束的方法。但由于非衍射声束理论属于较新的声学原理,截止于目前,国内外学者对非衍射声束理论的研究还处于比较浅显的阶段,还没有研究能够将非衍射声束理论与聚焦超声设备相结合,也没有研究能够将声准直效应应用于超声聚焦理论。
经检索,关于提高聚焦效果的聚焦设备已有相关专利公开。
如,中国专利申请号:201510010200.X,申请日为:2015年01月08日,发明创造名称:一种基于位移本征模式的声波聚焦器件,该申请案主要是通过将金属环状结构浸没在水中作为声聚焦器件,从而将焦点的形状由椭球形压缩到一个平面附近范围内,在一定程度上可以提高声波的聚焦效果,但该申请案并没有涉及声透镜的设计和制备,其原理与本发明不同,且其对径向旁瓣幅度的抑制效果及对焦点增益的提升效果也没有提及。
又如,中国专利申请号:201210142895.3,申请日为:2012年05月07日,发明创造名称为:一种聚焦方法和装置,该申请案利用时间反转理论能够实现自适应的高精度聚焦,提高了聚焦超声的精度和超声治疗的准确性。但该申请案也没有涉及声透镜的设计和制备,其主要是针对电子相控阵列聚焦方式进行的改进,虽能够在一定程度上增强焦点的超声强度,但不能够抑制焦点处径向旁瓣的幅度,也不能减小声焦点与几何焦点之间的偏移,从而不可避免地会对目标治疗区域周围的人体组织造成一定的损伤。
发明内容
1.发明要解决的技术问题
本发明的目的在于克服超声聚焦换能器采用现有声透镜进行聚焦时,焦点附近的径向旁瓣幅度较大,在超声治疗中会对目标区域周边的人体组织造成一定损伤的不足,提供了一种聚焦声透镜的设计方法。通过使用本发明的设计方法设计的聚焦声透镜,能够利用声波反常透射现象激发声准直效应,从而将能量更好地聚集到焦点区域,大大降低了焦点径向旁瓣的幅度,且提高了焦点的声波聚焦增益,从而提高了超声成像和超声治疗的准确性与安全性。
2.技术方案
为达到上述目的,本发明提供的技术方案为:
本发明的一种聚焦声透镜的设计方法,其步骤为:
步骤一、选定聚焦声透镜的初始目标工作频率f0,根据该初始目标工作频率f0计算出声波在工作介质中的纵波波长λ;在聚焦声透镜表面的声辐射面,即凹球面上设置周期性分布的同心环状凹槽结构,并设置该凹槽结构的初始分布参数:相邻凹槽的间隔为L0,凹槽的深度为h0,凹槽的宽度为d0;
步骤二、在有限元模型计算中循环调节凹槽的宽度d和深度值h,并采用频率扫描的方法分别计算焦点声压及透镜的声透射率相对于频率f的变化关系,当产生焦点声压最大的工作频率f1max与声波反常透射的峰值频率f2max的差值|f1max-f2max|达到预期要求时即停止循环,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1,并以此时的凹槽深度和宽度分别作为本次循环中优化的凹槽深度h1和宽度d1;
步骤三、在目标工作频率f1附近分析焦点处的径向旁瓣相对于主瓣的幅度比值,在有限元模型计算中保持优化后的凹槽深度h1和宽度d1不变,不断调节相邻凹槽的间隔L使上述幅度比值达到最小,从而取得最佳的径向旁瓣幅度抑制效果,并以此时相邻凹槽的间隔L作为本次循环中优化的相邻凹槽的间隔L1;
步骤四、评价经步骤一至步骤三中优化后的聚焦声透镜在目标工作频率f1处的聚焦效果,聚焦效果包括以下三项指标:(a)焦点处的径向旁瓣幅度相对于主瓣幅度的比值;(b)焦点处的聚焦增益,即焦点声压相对于平面活塞超声换能器发射声压的增益;(c)声压焦点位置与透镜几何焦点位置之间的偏移量;
步骤五、如果步骤四中聚焦声透镜在目标工作频率f1处的聚焦效果未达到设定目标,则重复进行步骤二至步骤四,直至步骤四中聚焦声透镜在目标工作频率f1处的聚焦效果达到设定目标;
步骤六、根据步骤五中最终得到的凹槽宽度d、凹槽深度h、相邻凹槽的间隔L来制作聚焦声透镜。
更进一步地,步骤一中所有凹槽的圆心均位于聚焦声透镜的轴线上,且每个凹槽的底部与顶部的宽度保持弧度值一致。
更进一步地,步骤一中设置相邻凹槽间隔的初始值L0以及凹槽深度的初始值h0均为波长λ,凹槽宽度的初始值d0为半波长λ/2。
更进一步地,步骤二中凹槽的深度h与宽度d的具体优化过程为:
(1)先保持凹槽的宽度d不变,以初始深度h0的1%为步长增大凹槽的深度h,若|f1max-f2max|减小,则继续以上述步长增大凹槽深度h,反之则以上述步长持续减小凹槽的深度h,当|f1max-f2max|开始变大时则停止调整,以|f1max-f2max|开始变大前的凹槽深度作为新的凹槽深度h;
(2)使新的凹槽深度h保持不变,以初始宽度d0的1%为步长增大凹槽的宽度d,若|f1max-f2max|减小,则继续以上述步长增大凹槽宽度d,反之则以上述步长持续减小凹槽宽度d,当|f1max-f2max|开始变大时则停止调整,以|f1max-f2max|开始变大前的凹槽宽度作为新的凹槽宽度d;
(3)循环执行步骤(1)、步骤(2),当|f1max-f2max|达到预期要求时即停止循环。
更进一步地,在步骤(3)中,当|f1max-f2max|小于声波反常透射的峰值频率f2max的1%时停止循环,此时得到的凹槽深度h和宽度d值即为本次循环中优化的凹槽深度h1和宽度d1,且以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1;若上述条件始终无法满足,则以循环过程中使|f1max-f2max|的值最小时的凹槽深度h和宽度d值作为本次循环中优化的凹槽深度h1和宽度d1,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1。
更进一步地,步骤三中相邻凹槽间隔L的具体优化过程为:以凹槽初始间隔L0的0.5%为步长增大相邻凹槽的间隔L,如果焦点处的径向旁瓣与主瓣的幅度比值变小,则继续以上述步长增大相邻凹槽的间隔L;如上述比值变大,则以凹槽初始间隔L0的1%为步长持续减小相邻凹槽的间隔L;当调整至上述幅度比值开始变大时则停止调整,并以上述幅度比值开始变大之前的相邻凹槽间隔作为优化的相邻凹槽间隔L1。
3.有益效果
采用本发明提供的技术方案,与现有技术相比,具有如下显著效果:
(1)本发明的一种聚焦声透镜的设计方法,通过在传统声透镜的辐射表面设置周期性分布的同心环形凹槽,并通过有限元计算对环形凹槽的分布参数,即相邻凹槽的间隔L,凹槽的深度h及凹槽的宽度d进行优化设计,能够在目标频率范围引入声准直效应以抑制衍射效应的发生。通过采用本发明的设计方法得到的聚焦声透镜,能够使聚焦声场焦点处的径向旁瓣幅度得到明显抑制,大大减少了超声成像中的斑点和伪影现象以及超声治疗中对目标区域周边组织产生的副作用。使用本发明得到的聚焦声透镜能够将更多的声能量聚焦到主瓣区域,从而提高聚焦效率,有利于提高超声成像的对比度和超声治疗的效率,且能够减小声场焦点和透镜几何焦点之间的偏移量,使超声成像和超声治疗在深度方向上更加准确。
(2)本发明的一种聚焦声透镜的设计方法,其设计方案不需要对声透镜的结构作较大改变,只须用周期结构来修饰透镜的辐射表面,具有加工难度低、容易实现的优点。
附图说明
图1为本发明的一种聚焦声透镜的设计方法的设计流程图;
图2中的(a)为传统聚焦声透镜的聚焦示意图;
图2中的(b)为本发明的聚焦声透镜的聚焦示意图;
图2中的(c)为传统聚焦声透镜的立体结构示意图;
图2中的(d)为本发明的聚焦声透镜的立体结构示意图;
图3中的(a)为利用传统聚焦声透镜聚焦时焦平面上的归一化声压分布的有限元模拟结果示意图;
图3中的(b)为利用本发明的聚焦声透镜聚焦时焦平面上的归一化声压分布的有限元模拟结果示意图;
图3中的(c)为利用传统聚焦声透镜聚焦时焦平面上的归一化声压分布的实际测量结果示意图;
图3中的(d)利用本发明的聚焦声透镜聚焦时焦平面上的归一化声压分布的实际测量结果示意图;
图4为分别利用传统的聚焦声透镜和本发明的聚焦声透镜聚焦时焦点处的径向声压分布相对于主瓣幅度的归一化值的对比图;
图5为分别利用传统的聚焦声透镜和本发明的聚焦声透镜聚焦时焦点处径向第一旁瓣幅度相对于主瓣幅度的归一化值随频率变化的对比图;
图6为分别利用传统的聚焦声透镜和本发明的聚焦声透镜聚焦时声轴线上的声压分布相对于焦点最大声压的归一化值的分布对比图;
图7为分别利用传统的聚焦声透镜和本发明的聚焦声透镜聚焦时声焦点与透镜几何焦点之间的距离随工作频率变化的分布对比图;
图8为分别利用传统的聚焦声透镜和本发明的聚焦声透镜聚焦时焦点增益的分布对比图。
具体实施方式
为进一步了解本发明的内容,现结合附图和实施例对本发明作详细描述。
实施例1
本实施例的一种聚焦声透镜的设计方法,其设计流程如图1所示,其具体步骤为:
步骤一、选定聚焦声透镜的初始目标工作频率f0=547kHz,根据该初始目标工作频率f0计算出声波在工作介质(声传播介质)中的纵波波长λ=c/f0,c为声波在工作介质中的纵波声速。在聚焦声透镜表面的声辐射面,即凹球面上设置周期性分布的同心环状凹槽结构,所有凹槽均指向透镜的几何焦点,凹槽圆心均位于聚焦声透镜的轴线上,且每个凹槽的底部与顶部的宽度保持弧度值一致。设置该凹槽结构的初始分布参数:相邻凹槽的间隔(即凹槽的环向分布周期长度)为L0,凹槽的深度为h0,凹槽的宽度为d0,且有L0=h0=λ,d0=λ/2。
如图2(a)和图2(b)所示分别为传统聚焦声透镜和本实施例所设计的聚焦声透镜的聚焦示意图,其中d为凹槽的宽度,L为凹槽的环向分布周期长度,即相邻凹槽的间隔,h为凹槽的深度,r为透镜辐射面的半径,θ为透镜的张角;声透镜的几何半径为75mm,透镜张角40°,声透镜材料选用铝,其密度和纵波声速c分别为2700千克/立方米和6400米/秒;声传播介质为水,其密度和纵波声速分别为1000千克/立方米和1500米/秒,此时传播介质中的纵波波长λ约为2.74毫米。因此,本实施例中设计L0=h0=2.74mm,d0=λ/2=1.37。
步骤二、在有限元模型计算中循环调节凹槽的宽度d和深度值h,并采用频率扫描的方法分别计算焦点声压及透镜的声透射率相对于频率f的变化关系,当产生焦点声压最大的工作频率f1max与声波反常透射的峰值频率f2max的差值|f1max-f2max|达到预期要求时即停止循环,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1,并以此时的凹槽深度和宽度分别作为本次循环中优化的凹槽深度h1和宽度d1。本实施例中凹槽的深度h与宽度d的具体优化过程为:
(1)先保持凹槽的宽度d不变,以初始深度h0的1%为步长增大凹槽的深度h,若|f1max-f2max|减小,则保持上述步长不变继续增大凹槽深度h,反之则以上述步长持续减小凹槽的深度h,当|f1max-f2max|开始变大时则停止调整,以|f1max-f2max|开始变大前的凹槽深度作为新的凹槽深度h。
(2)使新的凹槽深度h保持不变,以初始宽度d0的1%为步长增大凹槽的宽度d,若|f1max-f2max|减小,则保持上述步长不变继续增大凹槽宽度d,反之则以上述步长持续减小凹槽宽度d,当|f1max-f2max|开始变大时则停止调整,以|f1max-f2max|开始变大前的凹槽宽度作为新的凹槽宽度d。
(3)循环执行步骤(1)、步骤(2),当|f1max-f2max|小于声波反常透射的峰值频率f2max的1%时停止循环,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1,并以此时的凹槽深度和宽度分别作为本次循环中优化的凹槽深度h1和宽度d1。若上述条件始终无法满足,则以循环过程中使|f1max-f2max|的值最小时的凹槽深度h和宽度d值作为本次循环中优化的凹槽深度h1和宽度d1,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1。
步骤三、在目标工作频率f1附近(本实施例是在0.9f1~1.1f1范围内)分析焦点处的径向旁瓣相对于主瓣的幅度比值,在有限元模型计算中保持优化后的凹槽深度h1和宽度d1不变,不断调节相邻凹槽的间隔L使上述幅度比值达到最小,从而取得最佳的径向旁瓣幅度抑制效果,并以此时相邻凹槽的间隔L作为本次循环中优化的相邻凹槽的间隔L1。本实施例中相邻凹槽间隔L的具体优化过程为:以凹槽初始间隔L0的0.5%为步长增大相邻凹槽的间隔L,如果焦点处的径向旁瓣与主瓣的幅度比值变小,则继续以上述步长增大相邻凹槽的间隔L;如上述比值变大,则以凹槽初始间隔L0的1%为步长持续减小相邻凹槽的间隔L;当调整至上述幅度比值开始变大时则停止调整,并以上述幅度比值开始变大之前的相邻凹槽间隔作为优化的相邻凹槽间隔L1。
步骤四、评价经步骤一至步骤三中优化后的聚焦声透镜在目标工作频率f1处的聚焦效果,聚焦效果包括以下三项指标:(a)焦点处的径向旁瓣幅度相对于主瓣幅度的比值;(b)焦点处的聚焦增益,即焦点声压相对于平面活塞超声换能器发射声压的增益;(c)声压焦点位置与透镜几何焦点位置之间的偏移量。具体而言,在有限元模型中计算出上述三项指标,并与设定目标相比较,本实施例中三项指标的设定目标分别为:(a)焦点处的径向旁瓣幅度相对于主瓣幅度的比值相较于使用传统透镜时至少下降3dB;(b)焦点处的聚焦增益不低于使用传统透镜的情形;(c)声压焦点位置与透镜几何焦点位置之间的偏移量不超过0.5mm。
步骤五、如果步骤四中聚焦声透镜在目标工作频率f1处的聚焦效果未达到设定目标,则重复进行步骤二至步骤四,直至步骤四中聚焦声透镜在目标工作频率f1处的聚焦效果达到设定目标。本实施例中经有限元模拟计算并优化设计后,得到凹槽的最终设计尺寸为:凹槽宽度d为1.05mm,深度h为2.07mm,相邻凹槽的间隔L为2.72mm。
步骤六、根据步骤五中最终得到的凹槽宽度d、凹槽深度h、相邻凹槽的间隔L来制作聚焦声透镜,即得本实施例的聚焦声透镜。
为了证明本实施例所得聚焦声透镜的聚焦效果,同时制备一个传统聚焦声透镜,其材质及尺寸均与本实施例的聚焦声透镜完全一致,唯一的区别是其辐射表面没有设置周期分布的环形凹槽,两只透镜的立体示意图如图2(c)和图2(d)所示。将上述两种聚焦声透镜分别与一只平面活塞超声换能器(美国Panametrics公司,型号V389-SU,中心频率500kHz)利用白凡士林作为耦合剂进行贴合,浸没在水槽中(长60厘米,宽45厘米,高40厘米)进行实验。将一只针式水听器(美国ONDA公司,型号HNC-0400)装配在三维移动控制系统上进行声压扫描测量,信号由一台数字示波器进行采集和数字化。
如图3所示为两种不同透镜焦平面上(焦点处垂直于声传播方向的平面)的归一化声压分布,由图3可以看出,本实施例所提供聚焦声透镜相比传统聚焦声透镜具有更好的聚焦效果,其旁瓣幅度得到明显降低。
如图4所示为采用传统声透镜和本实施例的声透镜聚焦时焦点处径向声压分布相对于主瓣幅度的归一化值。其中,实线为利用传统聚焦声透镜的有限元模拟结果,虚线为利用本实施例的聚焦声透镜的有限元模拟结果,实心圆为利用传统聚焦声透镜进行实验的测量结果,空心圆为利用本实施例的聚焦声透镜进行实验的测量结果。由图4可以看出,使用传统聚焦声透镜时第一个旁瓣的归一化值的理论计算值为0.3,实验测量值为0.32;使用本实施例的聚焦声透镜时第一个旁瓣的归一化值的理论计算值为0.2,实验测量值为0.24,相对于传统聚焦声透镜其第一旁瓣的相对幅度可下降3分贝左右。因此,采用本实施例的聚焦声透镜能够有效抑制焦点处径向旁瓣的幅度。
如图5所示为采用传统聚焦声透镜和本实施例的聚焦声透镜时焦点处径向第一旁瓣幅度相对于主瓣幅度的归一化值随频率的变化,其中,实线为利用传统聚焦声透镜的有限元模拟结果,虚线为利用本实施例的聚焦声透镜的有限元模拟结果,实心圆为利用传统聚焦声透镜进行实验的测量结果;空心圆为利用本实施例的聚焦声透镜进行实验的测量结果。由图5可以看出本实施例的聚焦声透镜相比传统声透镜,可在520kHz到560kHz范围内有效降低焦点处的径向旁瓣幅度值。
如图6所示为声轴线上的声压分布相对于焦点最大声压的归一化值,其中,实线为利用传统聚焦声透镜的有限元模拟结果,虚线为利用本实施例的聚焦声透镜的有限元模拟结果,实心圆为利用传统聚焦声透镜进行实验的测量结果;空心圆为利用本实施例的聚焦声透镜进行实验的测量结果。由图6可以看出,利用本实施例的聚焦声透镜后,焦点声压幅度最大点相比于利用传统聚焦声透镜时更靠近声源。
如图7所示为声焦点与透镜几何焦点之间的距离随工作频率的变化,其中,实线为利用传统聚焦声透镜的有限元模拟结果,虚线为利用本实施例的聚焦声透镜的有限元模拟结果,实心圆为利用传统聚焦声透镜进行实验的测量结果;空心圆为利用本实施例的聚焦声透镜进行实验的测量结果。由图7可以看出,在实际目标工作频率541.5kHz处,利用本实施例的聚焦声透镜时声焦点和透镜几何焦点之间的距离约为0.3mm,而利用传统聚焦声透镜时该距离约为4.7mm。因此,使用本实施例的声透镜能够显著降低声压焦点位置与透镜几何焦点位置之间的偏移量。需要指出的是,虽然在450kHz频率处,两类透镜都能够使声焦点与透镜几何焦点之间的距离接近0,但是从图8可以发现,此时透镜的声透射率非常低,并不适合用于超声激励。
如图8所示为采用传统聚焦声透镜和本实施例的聚焦声透镜时的焦点增益(焦点最大声压值相对于平面活塞超声换能器发射声压的归一化值,单位为分贝)。其中,实线为利用传统聚焦声透镜的有限元模拟结果,虚线为利用本实施例的聚焦声透镜的有限元模拟结果,实心圆为利用传统聚焦声透镜进行实验的测量结果,空心圆为利用本实施例的聚焦声透镜进行实验的测量结果。由图8可以看出,在实际目标工作频率541.5kHz处,利用本实施例的聚焦声透镜时相比利用传统聚焦声透镜,其焦点增益可提升约5.5分贝。
综上所述,本实施例的一种聚焦声透镜的设计方法,通过在传统声透镜的辐射表面设置周期性分布的同心环形凹槽,并通过大量的研究,采用有限元模型对凹槽的分布参数进行了优化设计,从而能够在目标频率范围引入声准直效应以抑制衍射效应的发生,使聚焦声场焦点处的径向旁瓣幅度得到明显抑制,大大减少了超声成像中的斑点和伪影现象以及超声治疗中对目标治疗区域周边组织产生的副作用。使用本实施例得到的聚焦声透镜同时能够将更多的声能量聚焦到主瓣区域,从而提高聚焦效率,有利于提高超声成像的对比度和超声治疗的效率,且能够减小声场焦点和透镜几何焦点之间的偏移量,使超声成像和超声治疗在深度方向上更加准确。
Claims (6)
1.一种聚焦声透镜的设计方法,其步骤为:
步骤一、选定聚焦声透镜的初始目标工作频率f0,根据该初始目标工作频率f0计算出声波在工作介质中的纵波波长λ;在聚焦声透镜表面的声辐射面,即凹球面上设置周期性分布的同心环状凹槽结构,并设置该凹槽结构的初始分布参数:相邻凹槽的间隔为L0,凹槽的深度为h0,凹槽的宽度为d0;
步骤二、在有限元模型计算中循环调节凹槽的宽度d和深度值h,并采用频率扫描的方法分别计算焦点声压及透镜的声透射率相对于频率f的变化关系,当产生焦点声压最大的工作频率f1max与声波反常透射的峰值频率f2max的差值|f1max-f2max|达到预期要求时即停止循环,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1,并以此时的凹槽深度和宽度分别作为本次循环中优化的凹槽深度h1和宽度d1;
步骤三、在目标工作频率f1附近分析焦点处的径向旁瓣相对于主瓣的幅度比值,在有限元模型计算中保持优化后的凹槽深度h1和宽度d1不变,不断调节相邻凹槽的间隔L使上述幅度比值达到最小,从而取得最佳的径向旁瓣幅度抑制效果,并以此时相邻凹槽的间隔L作为本次循环中优化的相邻凹槽的间隔L1;
步骤四、评价经步骤一至步骤三中优化后的聚焦声透镜在目标工作频率f1处的聚焦效果,聚焦效果包括以下三项指标:(a)焦点处的径向旁瓣幅度相对于主瓣幅度的比值;(b)焦点处的聚焦增益,即焦点声压相对于平面活塞超声换能器发射声压的增益;(c)声压焦点位置与透镜几何焦点位置之间的偏移量;
步骤五、如果步骤四中聚焦声透镜在目标工作频率f1处的聚焦效果未达到设定目标,则重复进行步骤二至步骤四,直至步骤四中聚焦声透镜在目标工作频率f1处的聚焦效果达到设定目标;
步骤六、根据步骤五中最终得到的凹槽宽度d、凹槽深度h、相邻凹槽的间隔L来制作聚焦声透镜。
2.根据权利要求1所述的一种聚焦声透镜的设计方法,其特征在于:步骤一中所有凹槽的圆心均位于聚焦声透镜的轴线上,且每个凹槽的底部与顶部的宽度保持弧度值一致。
3.根据权利要求1或2所述的一种聚焦声透镜的设计方法,其特征在于:步骤一中设置相邻凹槽间隔的初始值L0以及凹槽深度的初始值h0均为波长λ,凹槽宽度的初始值d0为半波长λ/2。
4.根据权利要求1-3中任一项所述的一种聚焦声透镜的设计方法,其特征在于:步骤二中凹槽的深度h与宽度d的具体优化过程为:
(1)先保持凹槽的宽度d不变,以初始深度h0的1%为步长增大凹槽的深度h,若|f1max-f2max|减小,则保持上述步长不变继续增大凹槽深度h,反之则以上述步长持续减小凹槽的深度h,当|f1max-f2max|开始变大时则停止调整,以|f1max-f2max|开始变大前的凹槽深度作为新的凹槽深度h;
(2)使新的凹槽深度h保持不变,以初始宽度d0的1%为步长增大凹槽的宽度d,若|f1max-f2max|减小,则保持上述步长不变继续增大凹槽宽度d,反之则以上述步长持续减小凹槽宽度d,当|f1max-f2max|开始变大时则停止调整,以|f1max-f2max|开始变大前的凹槽宽度作为新的凹槽宽度d;
(3)循环执行步骤(1)、步骤(2),当|f1max-f2max|达到预期要求时即停止循环。
5.根据权利要求4所述的一种聚焦声透镜的设计方法,其特征在于:在步骤(3)中,当|f1max-f2max|小于声波反常透射的峰值频率f2max的1%时停止循环,此时得到的凹槽深度h和宽度d值即为本次循环中优化的凹槽深度h1和宽度d1,且以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1;若上述条件始终无法满足,则以循环过程中使|f1max-f2max|的值最小时的凹槽深度h和宽度d值作为本次循环中优化的凹槽深度h1和宽度d1,以此时产生焦点增益最大的工作频率f1max作为目标工作频率f1。
6.根据权利要求5所述的一种聚焦声透镜的设计方法,其特征在于:步骤三中相邻凹槽间隔L的具体优化过程为:以凹槽初始间隔L0的0.5%为步长增大相邻凹槽的间隔L,如果焦点处的径向旁瓣与主瓣的幅度比值变小,则继续以上述步长增大相邻凹槽的间隔L;如上述比值变大,则以凹槽初始间隔L0的1%为步长持续减小相邻凹槽的间隔L;当调整至上述幅度比值开始变大时则停止调整,并以上述幅度比值开始变大之前的相邻凹槽间隔作为优化的相邻凹槽间隔L1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510816714.4A CN105251140B (zh) | 2015-11-23 | 2015-11-23 | 一种聚焦声透镜的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510816714.4A CN105251140B (zh) | 2015-11-23 | 2015-11-23 | 一种聚焦声透镜的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105251140A true CN105251140A (zh) | 2016-01-20 |
CN105251140B CN105251140B (zh) | 2018-03-27 |
Family
ID=55091247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510816714.4A Active CN105251140B (zh) | 2015-11-23 | 2015-11-23 | 一种聚焦声透镜的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105251140B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462983A (zh) * | 2016-06-03 | 2017-12-12 | 三星电子株式会社 | 电磁波聚焦装置、光学设备和显微镜 |
CN107644636A (zh) * | 2017-09-01 | 2018-01-30 | 深圳先进技术研究院 | 一种基于平面人工结构的超声聚焦透镜 |
CN107812328A (zh) * | 2017-05-12 | 2018-03-20 | 菱可思株式会社 | 超声波装置和超声波单元 |
CN107870381A (zh) * | 2017-12-22 | 2018-04-03 | 深圳先进技术研究院 | 平面透镜和平面透镜的制作方法 |
CN108838747A (zh) * | 2018-08-10 | 2018-11-20 | 天津大学 | 一种基于声透镜的超声聚焦流体振动抛光系统 |
WO2019041299A1 (zh) * | 2017-09-01 | 2019-03-07 | 深圳先进技术研究院 | 一种基于平面人工结构的超声聚焦透镜 |
CN109431543A (zh) * | 2018-12-20 | 2019-03-08 | 深圳先进技术研究院 | 声波聚焦透镜、超声成像装置及方法 |
CN111112037A (zh) * | 2020-01-20 | 2020-05-08 | 重庆医科大学 | 透镜式多频聚焦超声换能器、换能系统及其声焦域轴向长度的确定方法 |
CN111326135A (zh) * | 2020-04-03 | 2020-06-23 | 青岛大学 | 一种宽频消色差声聚焦透镜 |
CN111407313A (zh) * | 2020-04-03 | 2020-07-14 | 深圳先进技术研究院 | 一种超声换能器和超声成像装置 |
CN112214926A (zh) * | 2020-09-21 | 2021-01-12 | 南京大学 | 具有超构表面声透镜的环形超声换能器及优化设计方法 |
WO2022179325A1 (zh) * | 2021-02-25 | 2022-09-01 | 西安邮电大学 | 一种用于超声声场测量系统的手动声轴对齐方法 |
CN117598728A (zh) * | 2024-01-23 | 2024-02-27 | 浙江大学 | 一种用于医学诊疗的平面型多焦点声透镜及声透镜装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102573983A (zh) * | 2009-08-14 | 2012-07-11 | 伊西康内外科公司 | 超声外科设备和硅波导以及其使用方法 |
US20120289827A1 (en) * | 2008-05-06 | 2012-11-15 | Ultrawave Labs, Inc. | Multi-Modality Ultrasound and Radio Frequency Methodology for Imaging Tissue |
WO2014176483A1 (en) * | 2013-04-26 | 2014-10-30 | Thync, Inc. | Focused transcranial ultrasound systems and methods for using them |
-
2015
- 2015-11-23 CN CN201510816714.4A patent/CN105251140B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120289827A1 (en) * | 2008-05-06 | 2012-11-15 | Ultrawave Labs, Inc. | Multi-Modality Ultrasound and Radio Frequency Methodology for Imaging Tissue |
CN102573983A (zh) * | 2009-08-14 | 2012-07-11 | 伊西康内外科公司 | 超声外科设备和硅波导以及其使用方法 |
WO2014176483A1 (en) * | 2013-04-26 | 2014-10-30 | Thync, Inc. | Focused transcranial ultrasound systems and methods for using them |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462983A (zh) * | 2016-06-03 | 2017-12-12 | 三星电子株式会社 | 电磁波聚焦装置、光学设备和显微镜 |
CN107812328A (zh) * | 2017-05-12 | 2018-03-20 | 菱可思株式会社 | 超声波装置和超声波单元 |
CN107644636A (zh) * | 2017-09-01 | 2018-01-30 | 深圳先进技术研究院 | 一种基于平面人工结构的超声聚焦透镜 |
WO2019041299A1 (zh) * | 2017-09-01 | 2019-03-07 | 深圳先进技术研究院 | 一种基于平面人工结构的超声聚焦透镜 |
US11650355B2 (en) | 2017-12-22 | 2023-05-16 | Shenzhen Institutes Of Advanced Technology | Planar lens and manufacturing method for planar lens |
CN107870381A (zh) * | 2017-12-22 | 2018-04-03 | 深圳先进技术研究院 | 平面透镜和平面透镜的制作方法 |
CN108838747A (zh) * | 2018-08-10 | 2018-11-20 | 天津大学 | 一种基于声透镜的超声聚焦流体振动抛光系统 |
CN109431543A (zh) * | 2018-12-20 | 2019-03-08 | 深圳先进技术研究院 | 声波聚焦透镜、超声成像装置及方法 |
CN111112037A (zh) * | 2020-01-20 | 2020-05-08 | 重庆医科大学 | 透镜式多频聚焦超声换能器、换能系统及其声焦域轴向长度的确定方法 |
CN111326135A (zh) * | 2020-04-03 | 2020-06-23 | 青岛大学 | 一种宽频消色差声聚焦透镜 |
CN111407313A (zh) * | 2020-04-03 | 2020-07-14 | 深圳先进技术研究院 | 一种超声换能器和超声成像装置 |
CN111326135B (zh) * | 2020-04-03 | 2023-07-07 | 青岛大学 | 一种宽频消色差声聚焦透镜 |
CN112214926A (zh) * | 2020-09-21 | 2021-01-12 | 南京大学 | 具有超构表面声透镜的环形超声换能器及优化设计方法 |
WO2022179325A1 (zh) * | 2021-02-25 | 2022-09-01 | 西安邮电大学 | 一种用于超声声场测量系统的手动声轴对齐方法 |
CN117598728A (zh) * | 2024-01-23 | 2024-02-27 | 浙江大学 | 一种用于医学诊疗的平面型多焦点声透镜及声透镜装置 |
CN117598728B (zh) * | 2024-01-23 | 2024-05-03 | 浙江大学 | 一种用于医学诊疗的平面型多焦点声透镜及声透镜装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105251140B (zh) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105251140B (zh) | 一种聚焦声透镜的设计方法 | |
Bessonova et al. | Focusing of high power ultrasound beams and limiting values of shock wave parameters | |
CN106267592B (zh) | 一种具有超构表面的凹面聚焦换能器设计方法 | |
US20120143100A1 (en) | Extended depth-of-focus high intensity ultrasonic transducer | |
CN109745077B (zh) | 基于聚焦超声声振信号的弹性特性检测方法 | |
Khokhlova et al. | Design of HIFU transducers to generate specific nonlinear ultrasound fields | |
Stephens et al. | Design aspects of focal beams from high-intensity arrays | |
Firouzi et al. | Efficient transcranial ultrasound delivery via excitation of lamb waves: Concept and preliminary results | |
Lu et al. | Enhanced-cavitation heating protocols in focused ultrasound surgery with broadband split-focus approach | |
CN103164588B (zh) | 一种大功率聚能反射罩的设计方法 | |
Qiao et al. | Transcostal high-intensity focused ultrasound treatment using phased array with geometric correction | |
Hunt | Principles of ultrasound used for hyperthermia | |
CN112214926B (zh) | 具有超构表面声透镜的环形超声换能器及优化设计方法 | |
Fan et al. | Acoustic characterization of high intensity focused ultrasound field generated from a transmitter with large aperture | |
Karzova et al. | Nonlinear effects in ultrasound fields of diagnostic-type transducers used for kidney stone propulsion: Characterization in water | |
Ming-He et al. | Nonlinear effect on focusing gain of a focusing transducer with a wide aperture angle | |
Li et al. | Adjustable acoustic field controlled by “ultrasonic projector” on ultrasound application | |
CN117598728B (zh) | 一种用于医学诊疗的平面型多焦点声透镜及声透镜装置 | |
Baker | Finite amplitude propagation of focused ultrasonic waves in water | |
Okada et al. | Influence of tough hydrophone shapes with titanium front plate and hydrothermal PZT thick film on distribution of acoustic bubbles around focal point of HIFU transducer | |
Igarashi et al. | Control of ultrasonic acoustic fields by multiple acoustic waveguides and piezoelectric transducers | |
Karaböce et al. | Investigation of ultrasonic fields produced by HIFU transducers used in cancer therapy | |
Ting-Bo et al. | Effects of vapour bubbles on acoustic and temperature distributions of therapeutic ultrasound | |
WU et al. | Numrical simulation of nonlinear focused ultrasound | |
JP2005005923A (ja) | 超音波送波器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220119 Address after: 210000 floor 21, block B, Zhongdan Life Science Industrial Park, No. 3-1, xinjinhulu, Jiangbei new area, Nanjing, Jiangsu Province Patentee after: JIANGSU TINGSN TECHNOLOGY Co.,Ltd. Address before: 210093 No. 22, Hankou Road, Gulou District, Jiangsu, Nanjing Patentee before: NANJING University |