一种基于剪切波幅度和相位检测的组织粘弹性测量方法

摘要

本发明提供一种基于剪切波幅度和相位检测的组织粘弹性测量方法，使用声辐射力激励和超声回波检测两个探头同时进行检测，声辐射力激励探头产生推动波束，使组织产生一定频率的剪切波，超声回波检测探头产生检测波束，检测剪切波在传播路径上的幅度衰减和相位改变，通过幅度衰减系数和速度检测值估计出组织的弹性模量和粘性系数，即可以同时定量检测出组织的弹性模量和粘性系数，为多种疾病的诊断提供了精确丰富的力学参数信息，有利于拓展超声弹性成像的临床应用范围与层次水平；同时本发明仅使用一种频率剪切波即可进行检测，降低了系统复杂性和设计成本，提高检测效率，还避免高频率剪切波引起的相位检测噪声，提高力学组织参数检测的准确性。

权利要求

1.一种基于剪切波幅度和相位检测的组织粘弹性测量方法，其特征在于，所述测量方法包括：S1、使用声辐射力激励探头产生推动波束，激励组织产生频率为ω的剪切波；

S2、在剪切波传播方向上选择相距为Δr的两个检测点r1和r2，且Δr小于剪切波的波长λ；

S3、对检测点r1，使用超声回波检测探头按设定帧率f产生检测波束，焦点位于r1上，通过回波分析获取r1质点随时间-振动曲线1；

S4、使用卡尔曼滤波从时间-振动曲线1抽取频率为ω的信号，该信号的相位和幅度即为检测点r1处相位φ1和幅度A1；

S5、对检测点r2，使用超声回波检测探头按同样设定帧率f产生检测波束，焦点位于r2上，通过回波分析获取r2质点随时间-振动曲线2；

S6、使用卡尔曼滤波从时间-振动曲线2抽取频率为ω的信号，该信号的相位和幅度即为检测点r2处相位φ2和幅度A2；

S7、使用下式计算得出剪切波速度c：

c＝ωΔr/(φ1-φ2)  式(2)

S8、使用下式计算得出剪切波幅度衰减系数α：

α＝ln(A1/A2)/Δr  式(3)

S9、根据所述剪切波幅度衰减系数α、剪切波频率ω、组织密度ρ、剪切弹性模量μ1和粘性系数μ2，建模成如下关系：S10、根据已知的计算剪切波速度的公式 和所述

式(4)求解可得：

在式(5)和式(6)中，剪切波频率ω和组织密度ρ是已知量，利用式(2)求解的剪切波速度c和式(3)求解的剪切波幅度衰减系数α，即可根据式(5)和式(6)分别计算出组织的剪切弹性模量μ1和粘性系数μ2。

2.根据权利要求1所述的基于剪切波幅度和相位检测的组织粘弹性测量方法，其特征在于，所述声辐射力激励探头采用3MHz聚焦探头，直径为45mm，组织重复推动频率为50Hz，所述超声回波检测探头采用5MHz聚焦探头，直径为12.5mm，脉冲回波检测频率为2KHz。

说明书

一种基于剪切波幅度和相位检测的组织粘弹性测量方法

技术领域

[0001] 本发明涉及医学超声工程技术领域，具体涉及一种基于剪切波幅度和相位检测的组织粘弹性测量方法。

背景技术

[0002] 超声弹性成像是近年来提出并已在临床诊断上显示了较大应用前景。组织的病理改变通常会引起弹性(硬度)变化，如肝纤维化、肝硬化、乳腺癌、前列腺癌等。因此组织的弹
性成像可成为一些重大疾病诊断的依据。

[0003] 然而组织是典型的粘弹性体，除了弹性以外，还具有粘性特征。粘性也与组织的病理状态相关，粘性的检测在医学中也具有重要意义。但是本发明的发明人发现，当前许多弹
性成像方法仅提供组织的弹性参数值，在估计组织的弹性参数时将组织假设成线弹性体，
忽略了组织的粘性。这种假设不仅无法提供组织的粘性参数，而且使组织的弹性参数估计
与实际产生偏差，精度有限。

[0004] 剪切波频散超声激励技术是最近提出的一种同时估计组织弹性和粘性参数的新技术，它利用剪切波频散特性，检测时需要激励多种频率的剪切波。原理如下：

[0005]

[0006] 其中，μ1和μ2分别是需要估计的剪切弹性模量与粘性系数，ρ是组织密度，ω是剪切波频率，c是剪切波速度。对于各种软组织来说，其密度ρ差别很小，可以近似成常数1000千
克/立方米。因此，只需要检测出不同频率ω剪切波的速度c，使用式(1)进行拟合即可求出
参数μ1和μ2。

[0007] 但是，本发明的发明人经过研究发现，当剪切波频率较高时，速度检测将产生较强噪声，致使参数μ1和μ2的求解精度受限；另外，此方法需要激励组织产生多种频率剪切波，因
而超声探头设计复杂，性能要求高，检测速度慢。

发明内容

[0008] 针对现有超声弹性成像仅提供组织的弹性参数信息，剪切波频散超声激励技术需要激励组织产生多种频率剪切波，当剪切波频率较高时检测精度受限，以及检测速度慢的
技术问题，本发明提供一种基于剪切波幅度和相位检测的组织粘弹性测量方法，本发明将
组织建模成与组织本身特性更接近的粘弹性模型，仅利用一种频率的剪切波能同时检测组
织的弹性和粘性两种特性参数。

[0009] 为了解决上述技术问题，本发明采用了如下的技术方案：

[0010] 一种基于剪切波幅度和相位检测的组织粘弹性测量方法，所述测量方法包括：

[0011] S1、使用声辐射力激励探头产生推动波束，激励组织产生频率为ω的剪切波；

[0012] S2、在剪切波传播方向上选择相距为Δr的两个检测点r1和r2，且Δr小于剪切波的波长λ；

[0013] S3、对检测点r1，使用超声回波检测探头按设定帧率f产生检测波束，焦点位于r1上，通过回波分析获取r1质点随时间振动曲线1；

[0014] S4、使用卡尔曼滤波从时间-振动曲线1抽取频率为ω的信号，该信号的相位和幅度即为检测点r1处相位φ1和幅度A1；

[0015] S5、对检测点r2，使用超声回波检测探头按同样设定帧率f产生检测波束，焦点位于r2上，通过回波分析获取r2质点随时间振动曲线2；

[0016] S6、使用卡尔曼滤波从时间-振动曲线2抽取频率为ω的信号，该信号的相位和幅度即为检测点r2处相位φ2和幅度A2；

[0017] S7、使用下式计算得出剪切波速度c：

[0018] c＝ωΔr/(φ1-φ2)  式(2)

[0019] S8、使用下式计算得出剪切波幅度衰减系数α：

[0020] α＝ln(A1/A2)/Δr  式(3)

[0021] S9、根据所述剪切波幅度衰减系数α、剪切波频率ω、组织密度ρ、剪切弹性模量μ1和粘性系数μ2，建模成如下关系：

[0022]

[0023] S10、根据已知的计算剪切波速度的公式 和所述式(4)求解可得：

[0024]

[0025]

[0026] 在式(5)和式(6)中，剪切波频率ω和组织密度ρ是已知量，利用式(2)求解的剪切波速度c和式(3)求解的剪切波幅度衰减系数α，即可根据式(5)和式(6)分别计算出组织的
剪切弹性模量μ1和粘性系数μ2。

[0027] 进一步，所述声辐射力激励探头采用3MHz聚焦探头，直径为45mm，组织重复推动频率为50Hz，所述超声回波检测探头采用5MHz聚焦探头，直径为12.5mm，脉冲回波检测频率为
2KHz。

[0028] 与现有技术相比，本发明提供的基于剪切波幅度和相位检测的组织粘弹性测量方法，使用声辐射力激励探头和超声回波检测探头两个探头同时进行检测，声辐射力激励探
头用于产生推动波束，使组织产生一定频率的剪切波，超声回波检测探头产生检测波束，检
测剪切波在传播路径上的幅度衰减和相位改变，通过幅度衰减系数和速度检测值估计出组
织的剪切弹性模量μ1和粘性系数μ2，即可以同时定量检测出组织的弹性模量和粘性系数，为
多种疾病的诊断提供了精确丰富的力学参数信息，有利于拓展超声弹性成像的临床应用范
围与层次水平；同时，本发明克服了现有技术中需要使用多种频率剪切波进行检测的缺点，
仅使用一种频率剪切波即可进行检测，降低了系统复杂性和设计成本，提高了检测效率；另
外，本发明还避免了高频率剪切波引起的相位检测噪声，提高了力学组织参数检测的准确
性。

附图说明

[0029] 图1是本发明提供的基于剪切波幅度和相位检测的组织粘弹性测量方法流程示意图。

[0030] 图2是本发明提供的剪切波相位和幅度检测原理示意图。

具体实施方式

[0031] 为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体图示，进一步阐述本发明。

[0032] 在本发明的描述中，需要理解的是，术语“纵向”、“径向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不
能理解为对本发明的限制。在本发明的描述中，除非另有说明，“多个”的含义是两个或两个
以上。

[0033] 在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可
以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是
两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本
发明中的具体含义。

[0034] 请参考图1和图2所示，本发明提供一种基于剪切波幅度和相位检测的组织粘弹性测量方法，所述测量方法包括：

[0035] S1、使用声辐射力激励探头产生推动波束，激励组织产生频率为ω的剪切波；

[0036] S2、在剪切波传播方向上选择相距为Δr的两个检测点r1和r2，且Δr小于剪切波的波长λ；

[0037] S3、对检测点r1，使用超声回波检测探头按设定帧率f产生检测波束，焦点位于r1上，通过回波分析获取r1质点随时间振动曲线1；

[0038] S4、使用卡尔曼滤波从时间-振动曲线1抽取频率为ω的信号，该信号的相位和幅度即为检测点r1处相位φ1和幅度A1；

[0039] S5、对检测点r2，使用超声回波检测探头按同样设定帧率f产生检测波束，焦点位于r2上，通过回波分析获取r2质点随时间振动曲线2；

[0040] S6、使用卡尔曼滤波从时间-振动曲线2抽取频率为ω的信号，该信号的相位和幅度即为检测点r2处相位φ2和幅度A2；

[0041] S7、使用下式计算得出剪切波速度c：

[0042] c＝ωΔr/(φ1-φ2)  式(2)

[0043] S8、使用下式计算得出剪切波幅度衰减系数α：

[0044] α＝ln(A1/A2)/Δr  式(3)

[0045] S9、根据所述剪切波幅度衰减系数α、剪切波频率ω、组织密度ρ、剪切弹性模量μ1和粘性系数μ2，建模成如下关系：

[0046]

[0047] S10、根据已知的计算剪切波速度的公式 和所述式(4)求解可得：

[0048]

[0049]

[0050] 在式(5)和式(6)中，剪切波频率ω和组织密度ρ是已知量，利用式(2)求解的剪切波速度c和式(3)求解的剪切波幅度衰减系数α，即可根据式(5)和式(6)分别计算出组织的
剪切弹性模量μ1和粘性系数μ2。

[0051] 作为具体实施例，所述声辐射力激励探头和超声回波检测探头可以采用本领域技术人员熟知的，或者现有技术中的探头来实现相应功能，例如所述声辐射力激励探头采用
3MHz聚焦探头，直径为45mm，组织重复推动频率为50Hz(HP33120A，Agilent，Santa Clara，
CA)，所述超声回波检测探头采用5MHz聚焦探头，直径为12.5mm，脉冲回波检测频率为2KHz
(i3-0508-R-SU，Olympus NDT，Waltham，MA)。

[0052] 与现有技术相比，本发明提供的基于剪切波幅度和相位检测的组织粘弹性测量方法，使用声辐射力激励探头和超声回波检测探头两个探头同时进行检测，声辐射力激励探
头用于产生推动波束，使组织产生一定频率的剪切波，超声回波检测探头产生检测波束，检
测剪切波在传播路径上的幅度衰减和相位改变，通过幅度衰减系数和速度检测值估计出组
织的剪切弹性模量μ1和粘性系数μ2，即可以同时定量检测出组织的弹性模量和粘性系数，为
多种疾病的诊断提供了精确丰富的力学参数信息，有利于拓展超声弹性成像的临床应用范
围与层次水平；同时，本发明克服了现有技术中需要使用多种频率剪切波进行检测的缺点，
仅使用一种频率剪切波即可进行检测，降低了系统复杂性和设计成本，提高了检测效率；另
外，本发明还避免了高频率剪切波引起的相位检测噪声，提高了力学组织参数检测的准确
性。

[0053] 最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技
术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本
发明的权利要求范围当中。