CN103990233A - 高强度超声经颅聚焦的相控方法 - Google Patents

高强度超声经颅聚焦的相控方法 Download PDF

Info

Publication number
CN103990233A
CN103990233A CN201410218215.0A CN201410218215A CN103990233A CN 103990233 A CN103990233 A CN 103990233A CN 201410218215 A CN201410218215 A CN 201410218215A CN 103990233 A CN103990233 A CN 103990233A
Authority
CN
China
Prior art keywords
focus
skull
transducer
signal
array element
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
Application number
CN201410218215.0A
Other languages
English (en)
Other versions
CN103990233B (zh
Inventor
菅喜岐
王意喆
张千
周文征
丁鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Medical University
Original Assignee
Tianjin Medical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin Medical University filed Critical Tianjin Medical University
Priority to CN201410218215.0A priority Critical patent/CN103990233B/zh
Publication of CN103990233A publication Critical patent/CN103990233A/zh
Application granted granted Critical
Publication of CN103990233B publication Critical patent/CN103990233B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Surgical Instruments (AREA)
  • Thermotherapy And Cooling Therapy Devices (AREA)

Abstract

本发明涉及一种高强度超声经颅聚焦的相控方法,特别是对相控换能器不同触发信号组的幅值、初始相位及组间时间延迟进行控制,解决经颅超声对颅骨及周边软组织损伤的问题。它包括以下步骤:S1:确定颅内靶区位置。S2:对聚焦颅内组织的相控换能器各阵元控制信号的触发声强、初始相位调控的阶段。S3:根据聚焦在颅内靶区形成的声压场与温度场分布,确定颅骨内声压和温升较高的热点位置的阶段。S4:计算消除颅骨热点的相控换能器各个阵元的控制信号的声强及时间延迟的阶段。S5:激发阵元,验证颅骨内有无热点的阶段。本发明可实现经颅聚焦,降低经颅超声对颅骨及其周边脑组织和头皮组织造成的损伤。

Description

高强度超声经颅聚焦的相控方法
技术领域
本发明涉及一种高强度超声经颅聚焦的相控方法,通过调控相控阵列换能器每个阵元激励信号的方法,根据颅骨内高声压处即热点处的稳态声压幅值及振动相位,设定用于消除热点的每个阵元激励信号及不同触发信号组间的时间延迟,降低颅骨内热点处温升,达到避免热点对颅骨及周边软组织造成损伤的目的。
背景技术
目前,聚焦超声已经普遍用于工程、农业、医学等各个领域。高强度聚焦超声技术具有穿透性能好、穿透深度大和指向性强等优势,特别是在医学领域,作为一种无创或微创的治疗方法,受到普遍关注。高强度聚焦超声技术已成功用于如前列腺癌、乳腺癌、肝癌、子宫肌瘤等疾病的临床治疗中。然而,在颅脑肿瘤治疗中,由于颅骨具有三层孔隙物理结构,其声阻抗远远大于软组织,使超声经颅聚焦形成的声场与温度场往往不能形成可治疗焦域或形成的可治疗焦域偏离目标区域。
在经颅聚焦超声治疗领域中,主要用到两种换能器,一种是半球形相控换能器,另一种是球冠状凹球面相控换能器。其中,以色列Insightec公司(Insightec Ltd.)生产的Exablate 4000高强度聚焦超声治疗系统(Exablate 4000 High-Intensity Focused Ultrasound Therapy System)中,包含1000阵元的一个可用于经颅治疗的半球形相控阵列换能器,其较大的开口直径和曲率半径可扩大超声透过颅骨的声窗,降低颅骨单位面积内的超声能量。但该设备采用1000个阵元的相控换能器,换能器的驱动及控制电路复杂、成本高且其仅限于可调控范围小的脑深部组织治疗;对颅内浅表组织治疗时,在颅骨内形成温升高的热点,热点会对颅骨及周边软组织造成损伤。对于球冠状凹球面相控换能器,虽然其相控换能器所含阵元较少,驱动及控制电路相对简单,但受其几何形状的影响,该类型换能器发射的声束透过颅骨声窗较小,易造成超声能量沉积,引起颅骨及周边损伤。目前,市场中尚无可用于经颅治疗的球冠状凹球面相控换能器及相关治疗系统。
发明内容
本发明为解决上述现有技术存在问题而提供了一种高强度聚焦超声经颅聚焦的相控方法。
图2所示为球冠状凹球面相控换能器经颅聚焦的空间示意图,该方法利用相控换能器每个阵元触发声强和初始相位可控的特性及聚焦超声线性传播的特性,叠加具有时间延迟的多组相控换能器激励信号,降低颅骨处高声压,消除颅骨内部存在的热点,实现聚焦超声经颅对脑组织的聚焦及治疗。
本发明所采取的技术方案是:
一种高强度超声经颅聚焦的相控方法,包括如下步骤:
S1:确定颅内靶区位置。
S2:根据颅内靶区位置及相控换能器可变焦范围,选定治疗目标深度与颅骨和相控换能器的相对位置,在目标焦点处设置点声源SF,在一段时间内形成的声压场分布中,记录处于接收模式的相控换能器阵元n的声压信号SFn(t),对相控换能器每个阵元记录的声压信号按时间序列进行反转处理,利用函数逼近法,根据SFn(t)计算可聚焦在颅内靶区的相控换能器阵元n的激励信号的初始相位延迟△tn,根据治疗剂量设定可聚焦在颅内靶区的相控换能器的输入声强I0和相控换能器阵元激励信号的声压最大幅值P0
S3:以信号                                                激励相控换能器的阵元n,根据其形成的声压场和一定时间的温度场分布,确定颅骨内需进行消除的热点R位置及热点R处稳态声压PR(S0n(t))的最大振幅P0R和初始相位φ0R阶段。
S4:在热点R处设置点声源SR,在一段时间内形成的声压场分布中,记录处于接收模式的相控换能器的阵元n的声压信号记录SRn(t),根据SRn(t)计算可聚焦在热点处的相控换能器阵元n激励信号的初始相位延迟△tiRn(i为热点消除次数)。根据公式(1),计算可消除热点的相控换能器的激励声强Ii和对应的相控换能器阵元激励信号的最大幅值Pi
                                 (1)
其中,图1为公式(1)在该模型热点R处的系统输入声强与输出声压的线性关系。当输入声强分别设定为Ia、Ib,通过计算得到在热点R处稳态时的声压最大值PaR和PbR。以信号激励相控换能器的阵元n,同上述S3,根据其形成的声压场和一定时间的温度场分布,确定消除热点信号组在热点位置的稳态声压PiRn(Si(t))的初始相位φiR。根据公式(2),计算S0n(t)与SiRn(t)两换能器激励信号组间的时间延迟△ti。记阵元n热点R处消除激励信号为
                                      (2)
S5:以计算公式(3)为激励信号驱动相控换能器,当一次热点消除时,叠加触发信号S0n(t)和SiRn(t),重新触发换能器,验证形成的声压场与温度场超声经颅聚焦目标靶区颅骨内原有热点是否消除及是否产生新热点。若仍存在热点,重复S4,进行第Q次热点消除。
           (3)
本发明具有的优点和积极效果是:
a.本发明可在相控阵列换能器可变焦范围内选择性聚焦,克服颅骨对颅内形成焦点产生偏移或散焦的作用。
b.本发明涉及的热点消除方法,可针对颅骨内部选择的热点位置,降低颅骨内声压及温升,达到消除颅内温升高的热点,使颅骨及其周边脑组织和脑组织免于受到损伤。在声压场中,图4为热点消除前后声轴上声压最大值特性曲线,采用热点消除方法后颅骨内部最大声压与焦点最大声压比值为0.45,相比未进行热点消除降低49.3%。在温度场中,图5为热点消除前后声轴上温升特性曲线,采用热点消除后颅骨内部最高温升降低至10.31℃,相比未进行热点消除降低了61.1%。 
c.热点消除方法主要作用于热点区域,对比图4与图5中,可知消除热点前后焦点(x=80mm)处的声压幅值及最大温升的变化趋势小于热点(x=60.5mm)处声压与温升的变化趋势。而热点消除方法对焦域的影响较小,对比图8与图9,可知消除热点前后焦点区域的大小、形状相近。
d.采用本发明涉及的高强度超声经颅聚焦的相控方法,如图5所示,热点消除后轴上焦点的最高温升为31.4℃,大于温升18℃对脑组织造成凝固性坏死的温升阈值,可对目标靶区组织产生损伤,有治疗作用。
e.本发明可针对颅骨内热点,对消除颅骨热点的相控换能器每个阵元的正弦激励信号的幅值及初始相位延迟进行计算,并根据聚焦在颅内脑组织的相控换能器激励信号,设定激励信号组间的时间延迟。本方法缩短了制定高强度治疗超声经颅治疗治疗计划的时间,提高了制定效率。
f.本方法利用制作成本较少及驱动控制电路较简单的凹球面相控阵列换能器,实现高强度聚焦超声经颅聚焦的同时降低颅骨内热点对颅骨及周边软组织造成的损伤。
附图说明
图1为模型热点处的系统输入声强与输出声压的线性关系图;
图2为球冠状凹球面相控换能器经颅聚焦空间示意图;
图3为球冠状凹球面相控换能器经颅聚焦模型图;
图4为热点消除前后声轴上声压最大值特性曲线;
图5为热点消除前后声轴上15s最高温升特性曲线;
图6为热点消除前采用时间反转法X-Z平面声压场最大声压分布图;
图7为热点消除后X-Z平面声压场最大声压分布图;
图8为热点消除前采用时间反转法X-Z平面15s温度场最大温升分布图;
图9为热点消除后X-Z平面15s温度场最大温升分布图;
图10为本发明的高强度超声经颅聚焦的相控方法流程图。
图中主要标号说明:
1—凹球面64阵元相控阵列超声换能器 
2—人体颅骨模型
3—目标焦点   
4—脱气水
5—脑组织
6—换能器单阵元
Φ1—换能器开口直径     
Φ2—阵元直径
Φ3—换能器中心开孔直径  
θ—奇偶环层相邻阵元中心与换能器中心连线的夹角。
具体实施方式
以下结合附图和实例对本发明的技术方案作进一步说明描述。本发明的优选实例不构成对本发明的限制。
本发明提出的一种强度超声经颅聚焦的相控方法。如图3所示为实例模型,其中换能器以一球冠状凹球面64阵元相控阵列换能器在频率为0.8MHz条件下辐照为例,按如下步骤实施。
S1:确定颅内靶区位置;
S2:根据颅内靶区位置及相控换能器可变焦范围,选定治疗深度L2为25mm,颅骨与相控换能器的相对距离L1为55mm。在目标焦点处设置点声源,处于接收模式的相控换能器的每个阵元在一段时间内记录点声源激发的声压场。对换能器各个阵元记录的声压信号按时间序列进行反转处理,计算聚焦在颅内靶区的相控换能器阵元n的激励信号的初始相位延迟△tn,并根据治疗剂量设定相控换能器的输入声强I0,本实例中I0为2.5W/cm2
S3:采用信号激励相控换能器阵元n,根据其形成的声压场及一定时间的温度场分布,确定颅骨内热点位置为声轴上x=60.5mm处。
S4:在轴上热点x=60.5mm设置点声源,其在一段时间内形成的声压场分布由处于接收模式的相控换能器的各个阵元记录,计算可聚焦在热点处的相控换能器阵元n的激励信号的初始相位延迟△tni,根据公式(1),推算可消除热点的相控换能器激励信号的统一激励声强Ii为0.29W/cm2,记该组激励信号为
                                   (1)
在信号触发相控换能器形成的治疗声场和信号触发相控换能器形成消除热点的声场中,分别记录颅骨热点处记录声压稳态信号P(S0n(t))和PiRn(Si(t)),两组声压的正弦初始相位分别为,根据公式(2),计算可消除热点的相控换能器激励信号的时间延迟为1.78μs。
                                        (2)
S5:根据计算公式(3),叠加触发信号S0n(t)和SiRn(t),重新触发换能器,验证形成的声压场与温度场超声经颅聚焦目标靶区颅骨内原有热点是否消除及是否产生新热点。在本实例中,消除热点次数Q为1,即经一次消除后热点不再存在且无新热点形成,则利用该信号组进行经颅治疗。
 
        (3)
 在上述超声经颅聚焦消除颅骨热点的实例中,图4为热点消除前后声轴上声压最大值特性曲线,图6、7为对应的声压场分布图。如图4所示,热点消除前颅骨内最大声压与焦点最大声压比值(骨焦比)为0.89,热点消除后骨焦比降低到0.45,颅骨内最大声压降低幅度大于焦点最大声压变化幅度。图5为热点消除前后声轴上温升特性曲线,图8、9为对应的温度场分布图。如图5所示,热点消除前颅骨内最大温升为26.65℃,焦点最高温升为28.30℃,热点消除后颅骨内最大温升降低至为10.31℃,焦点最高温升为31.40℃。通过比较以上相关参数,验证本发明可提高经颅超声在焦点区域的聚焦效果的同时,消除颅骨内热点,降低对颅骨及周边组织损伤。

Claims (3)

1.一种高强度超声经颅聚焦的相控方法,其特征在于,该方法包括以下步骤:
S1:确定颅内靶区位置;
S2:根据颅内靶区位置及相控换能器可变焦范围,选定治疗目标深度与颅骨和相控换能器的相对位置,在目标焦点处设置点声源SF,在一段时间内形成的声压场分布中,记录处于接收模式的相控换能器阵元n的声压信号SFn(t),对相控换能器每个阵元记录的声压信号按时间序列进行反转处理,利用函数逼近法,根据SFn(t)计算可聚焦在颅内靶区的相控换能器阵元n的激励信号的初始相位延迟△tn,根据治疗剂量设定可聚焦在颅内靶区的相控换能器的输入声强I0和相控换能器阵元激励信号的声压最大幅值P0阶段;
S3:以信号                                                激励相控换能器的阵元n,根据其形成的声压场和一定时间的温度场分布,确定颅骨内需进行消除的热点R位置及热点R处稳态声压PR(S0n(t))的最大振幅P0R和初始相位φ0R阶段;
S4:在热点R处设置点声源SR,在一段时间内形成的声压场分布中,记录处于接收模式的相控换能器的阵元n的声压信号记录SRn(t),根据SRn(t)计算可聚焦在热点处的相控换能器阵元n激励信号的初始相位延迟△tiRn,根据公式(1)即,计算可消除热点的相控换能器的激励声强Ii和对应的相控换能器阵元激励信号的最大幅值Pi,以信号激励相控换能器的阵元n,同上述S3,根据其形成的声压场和一定时间的温度场分布,确定消除热点信号组在热点位置的稳态声压PiRn(Si(t))的初始相位φiR,根据公式(2)即 ,计算S0n(t)与SiRn(t)两换能器激励信号组间的时间延迟△ti,记阵元n热点R处消除的激励信号为
S5:以计算公式(3)即 为激励信号驱动相控换能器,当一次热点消除时,叠加触发信号S0n(t)和SiRn(t),重新触发换能器,验证形成的声压场与温度场超声经颅聚焦目标靶区颅骨内原有热点是否消除及是否产生新热点;若仍存在热点,重复S4,进行第Q次热点消除阶段。
2.根据权利要求1所述的高强度超声经颅聚焦脑组织的相控方法,其特征在于:基于线性系统输入与输出的特性,根据公式(1)即,计算消除颅骨内热点的相控换能器激励信号在热点处的稳态声压幅值近似为热点声压峰值Pi的输入声强Ii
3.根据权利要求1所述的高强度超声经颅聚焦的相控方法,其特征在于:在触发的治疗声场和触发的消除声场中,分别记录颅骨最大声压处P(S0n(t))和PiRn(Si(t))的声压稳态信号,用来近似描述两组声压的正弦初始相位分别为φ0R和φiR,根据公式(2)即,计算S0n(t)与SiRn(t)两换能器激励信号组间的时间延迟△ti
CN201410218215.0A 2014-05-22 2014-05-22 高强度超声经颅聚焦的相控方法 Active CN103990233B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410218215.0A CN103990233B (zh) 2014-05-22 2014-05-22 高强度超声经颅聚焦的相控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410218215.0A CN103990233B (zh) 2014-05-22 2014-05-22 高强度超声经颅聚焦的相控方法

Publications (2)

Publication Number Publication Date
CN103990233A true CN103990233A (zh) 2014-08-20
CN103990233B CN103990233B (zh) 2017-02-01

Family

ID=51304723

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410218215.0A Active CN103990233B (zh) 2014-05-22 2014-05-22 高强度超声经颅聚焦的相控方法

Country Status (1)

Country Link
CN (1) CN103990233B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105233430A (zh) * 2015-11-18 2016-01-13 天津医科大学 高强度超声经颅聚焦的颅外驻波抑制方法
CN107019853A (zh) * 2017-05-19 2017-08-08 天津医科大学 高强度超声聚焦焦域的信号叠加相控方法
CN107929960A (zh) * 2017-12-12 2018-04-20 天津医科大学 双激励信号叠加调控聚焦超声焦域平台式温度分布的方法
CN110191740A (zh) * 2016-12-22 2019-08-30 新宁研究院 用于执行经颅超声治疗和成像程序的系统和方法
CN111589001A (zh) * 2020-06-10 2020-08-28 浙江大学 一种阵列式超声治疗系统
CN112967709A (zh) * 2021-02-20 2021-06-15 河北工业大学 一种基于无损脑神经刺激聚焦声场的相控方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6023977A (en) * 1997-08-01 2000-02-15 Acuson Corporation Ultrasonic imaging aberration correction system and method
US6419648B1 (en) * 2000-04-21 2002-07-16 Insightec-Txsonics Ltd. Systems and methods for reducing secondary hot spots in a phased array focused ultrasound system
CN100475289C (zh) * 2006-01-26 2009-04-08 上海交通大学 相控阵聚焦超声的大焦域形成方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105233430A (zh) * 2015-11-18 2016-01-13 天津医科大学 高强度超声经颅聚焦的颅外驻波抑制方法
CN105233430B (zh) * 2015-11-18 2018-01-19 天津医科大学 高强度超声经颅聚焦的颅外驻波抑制方法
CN110191740A (zh) * 2016-12-22 2019-08-30 新宁研究院 用于执行经颅超声治疗和成像程序的系统和方法
CN110191740B (zh) * 2016-12-22 2021-12-21 新宁研究院 用于执行经颅超声治疗和成像程序的系统和方法
CN107019853A (zh) * 2017-05-19 2017-08-08 天津医科大学 高强度超声聚焦焦域的信号叠加相控方法
CN107929960A (zh) * 2017-12-12 2018-04-20 天津医科大学 双激励信号叠加调控聚焦超声焦域平台式温度分布的方法
CN107929960B (zh) * 2017-12-12 2019-06-21 天津医科大学 双激励信号叠加调控聚焦超声焦域平台式温度分布的方法
CN111589001A (zh) * 2020-06-10 2020-08-28 浙江大学 一种阵列式超声治疗系统
CN111589001B (zh) * 2020-06-10 2022-03-01 浙江大学 一种阵列式超声治疗系统
CN112967709A (zh) * 2021-02-20 2021-06-15 河北工业大学 一种基于无损脑神经刺激聚焦声场的相控方法

Also Published As

Publication number Publication date
CN103990233B (zh) 2017-02-01

Similar Documents

Publication Publication Date Title
US20240033542A1 (en) Histotripsy using very short ultrasound pulses
CN103990233A (zh) 高强度超声经颅聚焦的相控方法
US20220219019A1 (en) Histotripsy therapy systems and methods for the treatment of brain tissue
US10576304B2 (en) Thermal therapy apparatus and method using focused ultrasonic sound fields
US10806952B2 (en) Therapeutic ultrasound apparatus and method
CN107913477B (zh) 一种阵列超声换能器的激励方法、装置、设备及存储介质
CN104093452B (zh) 超声波像系统及超声波像生成方法
EP3016594A1 (en) Histotripsy excitation sequences optimized for bubble cloud formation using shock scattering
Nakamura et al. Large volume coagulation utilizing multiple cavitation clouds generated by array transducer driven by 32 channel drive circuits
US11896853B2 (en) Transrectal ultrasound probe for boiling histotripsy ablation of prostate, and associated systems and methods
CN104353196B (zh) 高强度超声经颅聚焦的幅值调制方法
US7645235B2 (en) Method for effecting local increases in temperature inside materials, particularly body tissue
CN104001277A (zh) 一种复频超声波肿瘤治疗头
CN203874298U (zh) 一种复频超声波肿瘤治疗头
CN108523923B (zh) 基于随机分布82阵元相控超声换能器多频分区激励方法
CN205127175U (zh) 高能聚焦超声肿瘤治疗累积热损伤自动保护避让系统
Liu et al. Close-loop lesion formation control using multiple-focus dual mode ultrasound array
CN105233430B (zh) 高强度超声经颅聚焦的颅外驻波抑制方法
Yeats et al. In Vivo Cavitation-based Aberration Correction of Histotripsy in Porcine Liver
Li et al. Dual-focus scanning in volumetric HIFU ablation: Preliminary simulation study
WO2019126991A1 (zh) 一种阵列超声换能器的激励方法、装置、设备及存储介质
Gateau et al. Reaching the optimal focusing and steering capabilities of transcranial HIFU arrays based on time reversal of acoustically induced cavitation bubble signature
Bigelow Improved heating efficiency with High-Intensity Focused Ultrasound using a new ultrasound source excitation
Cai et al. Treatment Planning of Ultrasound Intensity and Heating Time for Phased HIFU Surgery

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
OL01 Intention to license declared
OL01 Intention to license declared