CN106037816A - 一种超声粘弹性测量方法及系统 - Google Patents

一种超声粘弹性测量方法及系统 Download PDF

Info

Publication number
CN106037816A
CN106037816A CN201610334175.5A CN201610334175A CN106037816A CN 106037816 A CN106037816 A CN 106037816A CN 201610334175 A CN201610334175 A CN 201610334175A CN 106037816 A CN106037816 A CN 106037816A
Authority
CN
China
Prior art keywords
ultrasonic
tissue
signal
ultrasonic transducer
original
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.)
Pending
Application number
CN201610334175.5A
Other languages
English (en)
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.)
Shenzhen University
Original Assignee
Shenzhen 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 Shenzhen University filed Critical Shenzhen University
Priority to CN201610334175.5A priority Critical patent/CN106037816A/zh
Publication of CN106037816A publication Critical patent/CN106037816A/zh
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明提供了一种超声粘弹性测量方法及系统，方法包括：任意信号发生器同时产生原始组织激励信号和超声换能器激励信号；所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。本发明将编码检测技术应用到超声瞬时弹性成像系统，以增加检测深度。

Description

一种超声粘弹性测量方法及系统
技术领域
本发明涉及粘弹性测量技术领域，尤其涉及一种超声粘弹性测量方法及系统。
背景技术
通常组织发生病变会导致组织粘弹性的变化。因此，对组织的弹性、粘性进行定量的测量对疾病的临床诊断有一定的参考价值，弹性成像技术应运而生。弹性成像的基本原理是：对组织施加一个内部(包括自身的)或外部的动态或静态/准静态的激励，在弹性力学、生物力学等物理规律作用下，组织将产生一个响应，例如位移、应变、速度的分布有一定的差异。弹性模量较大即较硬的组织应变较小，或者振动的幅度较小或速度较大。利用超声成像、磁共振成像或者光学成像等方法，结合数字信号处理技术，可以估算出组织内部的响应情况，从而反映组织内部的弹性模量。
基于外加振子的瞬时弹性成像方法是一种安全、易于集成到现有商用超声系统上的技术。它主要采用机械振动器在体表，或者将振动器和超声探头集成在一起，通过振动器向体内传播低频振动，用脉冲回波检测振动引起的剪切波传播时的幅度、相位和速度信息，由此得出组织的弹性模量。但该方法存在检测深度有限的问题。
现有的外加振子的瞬时弹性成像方法主要有种方式：透射式超声瞬时弹性成像和反射式超声瞬时弹性成像。透射式超声瞬时弹性成像是将超声换能器固定在介质的一侧，用于发射和接收超声波，低频振荡器固定在介质的另一侧，两者都垂直于介质的表面，超声换能器固定在低频振荡器的轴向上，其作用类似于一个活塞式振荡器(振荡器的位移垂直于活塞表面)，用于产生低频剪切波，这种方法检测较为不方便，被测物体不能太厚。反射式超声瞬时弹性成像，超声换能器固定在低频振荡器的轴向上，低频振荡器带动超声换能器振动，在介质内部产生低频剪切波，超声换能器用于发射和接收超声波，同时起到活塞式振荡器的作用，这种方法剪切波所产生的绝对位移等于实验所测得的位移减去超声换能器自身振动的位移上。
上述方法均采用传统的超声脉冲检测回波信号的方法，检测深度存在一定的限制。
因此，现有技术还有待改进和发展。
发明内容
鉴于上述现有技术的不足之处，本发明的目的在于提供一种超声粘弹性测量方法及系统，旨在解决现有技术中超声脉冲检测回波信号的方法，检测深度存在限制的问题。
为了达到上述目的，本发明采取了以下技术方案：
一种超声粘弹性测量方法，其中，所述方法包括以下步骤：
A、任意信号发生器同时产生原始组织激励信号和超声换能器激励信号；
B、所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；
C、所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；所述步骤C与所述步骤B同步进行；
D、根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；
E、根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。
所述超声粘弹性测量方法，其中，所述步骤B具体包括：
B1、所述原始组织激励信号通过功率放大器放大，得到组织激励信号；
B2、所述组织激励信号加载在微型激励器上，通过微型激励器驱动振动杆对所述组织施加指定波形的振动，使组织内产生剪切波；其中，深度为d0的剪切波d(t)＝Asin(2πfvt+φv)+d0；其中，A是振动幅值，fv是振动频率，φv是起始相位，t是时间参数。
所述超声粘弹性测量方法，其中，所述超声换能器激励信号包括由任意信号发生器产生的第一超声换能驱动子信号，和与第一超声换能驱动子信号相位相差180度的第二超声换能驱动子信号。
所述超声粘弹性测量方法，其中，所述步骤C具体包括：
C1、所述第一超声换能驱动子信号和所述第二超声换能驱动子信号经超声收/发系统耦合后，得到超声换能器驱动信号；
C2、超声换能器接收所述超声换能器驱动信号，并发出原始超声波p(t)；其中，p(t)＝sin(2πfot)，f0是超声换能器的发射频率，t是时间参数；
C3、所述原始超声波经所述组织的反射和散射后，得到原始回波信号r(t，k)，并被所述超声换能器接收；其中，r(t，k)＝p(t-Δtk)，Δtk表示组织振动产生的延时，t是时间参数，k取任意正整数，θ是原始超声波的超声波束与组织振动方向的夹角，c是超声波声速；
C4、所述原始回波信号r(t，k)经超声收/发系统放大滤波滤掉高频后，得到回波信号H(t，k)；其中，
C5、采集卡接收所述回波信号H(t，k)。
所述超声粘弹性测量方法，其中，所述步骤D中脉冲回波处理算法为Voigt模型，且根据Voigt模型和回波信号得到剪切波在组织中的传播速度cs；其中，剪切波在组织中的频率为ωs，ρ是组织密度，μ1是弹性模量，μ2是粘性系数。
所述超声粘弹性测量方法，其中，所述步骤E中组织的粘弹性数据模型为其中E是弹性模量，G是杨氏模量，σ是泊松比，ρ是组织密度。
一种超声粘弹性测量系统，其中，包括：
任意信号发生器，用于同时产生原始组织激励信号和超声换能器激励信号；
组织激励模块，用于所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；
脉冲回波检测模块，用于所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；所述脉冲回波检测模块与所述组织激励模块同时启动；
剪切波速度计算模块，用于根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；
弹性参数获取模块，用于根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。
所述超声粘弹性测量系统，其中，所述组织激励模块具体包括：
功率放大器，用于将所述原始组织激励信号进行放大，得到组织激励信号；
微型激励器，用于接收所述组织激励信号并驱动振动杆；
振动杆，用于在微型激励器的驱动下对所述组织施加指定波形的振动，使组织内产生剪切波。
所述超声粘弹性测量系统，其中，所述超声换能器激励信号包括由任意信号发生器产生的第一超声换能驱动子信号，和与第一超声换能驱动子信号相位相差180度的第二超声换能驱动子信号。
所述超声粘弹性测量系统，其中，所述脉冲回波检测模块具体包括：
耦合单元，用于将所述第一超声换能驱动子信号和所述第二超声换能驱动子信号进行耦合，得到超声换能器驱动信号；
原始超声波发送单元，用于接收所述超声换能器驱动信号，并发出原始超声波p(t)；其中，p(t)＝sin(2πfot)，f0是超声换能器的发射频率，t是时间参数；
超声波接收单元，用于接收所述原始超声波经所述组织的反射和散射后而得到原始回波信号r(t，k)；其中，r(t，k)＝p(t-Δtk)，Δtk表示组织振动产生的延时，t是时间参数，k取任意正整数，θ是原始超声波的超声波束与组织振动方向的夹角，c是超声波声速；
高频滤波单元，用于将所述原始回波信号r(t，k)滤波滤掉高频后，得到回波信号H(t，k)；其中，
采集单元，用于接收所述回波信号H(t，k)；
其中，所述耦合单元、及所述高频滤波单元均设置在超声收/发系统中，所述原始超声波发送单元、及所述超声波接收单元均设置在超声换能器中，所述采集单元为采集卡。
本发明所述的超声粘弹性测量方法及系统，方法包括：任意信号发生器同时产生原始组织激励信号和超声换能器激励信号；所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。本发明将编码检测技术应用到超声瞬时弹性成像系统，以增加检测深度。
附图说明
图1为本发明所述超声粘弹性测量方法较佳实施例的流程图。
图2为本发明所述超声粘弹性测量系统较佳实施例的结构框图。
图3为本发明所述超声粘弹性测量系统具体实施例的示意框图。
具体实施方式
本发明提供一种超声粘弹性测量方法及系统，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
请参考图1，其为本发明所述超声粘弹性测量方法较佳实施例的流程图。如图1所示，所述超声粘弹性测量方法，包括以下步骤：
步骤S100、任意信号发生器同时产生原始组织激励信号和超声换能器激励信号；
步骤S200、所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；
步骤S300、所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；所述步骤S300与所述步骤S200同步进行；
步骤S400、根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；
步骤S500、根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。
本发明的实施例中，由任意信号发生器同时产生原始组织激励信号和超声换能器激励信号。之所以需要同时产生原始组织激励信号和超声换能器激励信号，是因为要确保步骤S200和步骤S300同步进行。步骤S200中振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波，但由于剪切波存在时间短且很快就会消失，故需要超声换能器同步发出用于检测的超声波，以实现对剪切波成分的检测和分析，故步骤S200和步骤S300同步进行，也即每次超声波检测剪切波时要保证组织中的激励波形(也即剪切波的波形)处在同样的状态。
具体实施时，所述步骤S100中，任意信号发生器中输出信号具有多样性，也即任意信号发生器提供了组织激励方式的多样性。例如，可以采用不同脉冲宽度的标准正弦激励组织(类似于突然敲打一下组织的作用)，也可以是N个周期的正弦脉冲，或者是某一频率的连续振动波形等。通过任意信号发生器产生不同形式的振动波形激励组织，使组织中的散射子做不同形式的振动以产生剪切波，有利于对组织粘弹性的测量。
所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器的检测脉冲，可以通过控制任意信号发生器灵活地产生。本发明不仅可以实现传统的脉冲回波检测，还可以实现编码检测，且该发明仅适合相位编码方式。
进一步的，所述步骤S200具体包括：
步骤S201、所述原始组织激励信号通过功率放大器放大，得到组织激励信号；
步骤S202、所述组织激励信号加载在微型激励器上，通过微型激励器驱动振动杆对所述组织施加指定波形的振动，使组织内产生剪切波；其中，深度为d0的剪切波d(t)＝Asin(2πfvt+φv)+d0；其中，A是振动幅值，fv是振动频率，φv是起始相位，t是时间参数。
d(t)还与外加机械振动有关，如果只是采用单独一个正弦波的振动形式，A将随时间衰减；如果采用连续正弦振动形式，A近似为常数，φv是起始相位。
进一步的，所述步骤S300具体包括：
步骤S301、所述第一超声换能驱动子信号和所述第二超声换能驱动子信号经超声收/发系统耦合后，得到超声换能器驱动信号。
步骤S302、超声换能器接收所述超声换能器驱动信号，并发出原始超声波p(t)；其中，p(t)＝sin(2πfot)，f0是超声换能器的发射频率，t是时间参数。
步骤S303、所述原始超声波经所述组织的反射和散射后，得到原始回波信号r(t，k)，并被所述超声换能器接收；其中，r(t，k)＝p(t-Δtk)，Δtk表示组织振动产生的延时，t是时间参数，k取任意正整数，θ是原始超声波的超声波束与组织振动方向的夹角，c是超声波声速。
步骤S304、所述原始回波信号r(t，k)经超声收/发系统放大滤波滤掉高频后，得到回波信号H(t，k)；其中，
对原始回波信号r(t，k)进行正交分解，得到信号z(t，k)，其中：应用一个线性相位低通滤波器将信号z(t,k)中的高频滤掉，得到回波信号H(t，k)；
其中，
H(t，k)的相位为
可知，组织的振动信号为当t＝t0＝2d0/c时(其中t0＝2d0/c表示超声波往返超声换能器与深度为d0处的散射子所需的时间)，得到深度为d0处的散射子的慢时振动信号d(t0+kT)。当采用根据不同深度下的组织振动慢时信号峰值出现的时刻不一样，即达峰时间，即可计算出振动在组织中的传播速度。
这里假设散射子做正弦振动，当散射子为其他形式的振动波形时，振动检测算法原理依然可以采用该过程。这一脉冲回波处理算法采用了传统的脉冲回波检测，即探头发射M个正弦波的形式，当采用编码检测方式时，在对回波信号进行正交解调之后，需要进行解码，后面的步骤不变。
步骤S305、采集卡接收所述回波信号H(t，k)。
进一步的，所述步骤S400中脉冲回波处理算法为Voigt模型，且根据Voigt模型和回波信号得到剪切波在组织中的传播速度cs；其中，剪切波在组织中的频率为ωs，ρ是组织密度，μ1是弹性模量，μ2是粘性系数。
进一步的，所述步骤S500中组织的粘弹性数据模型为其中E是弹性模量，G是杨氏模量，σ是泊松比，ρ是组织密度。
基于上述方法实施例，本发明还提供了一种超声粘弹性测量系统。如图2和图3所示，所述超声粘弹性测量系统包括：
任意信号发生器1-1，用于同时产生原始组织激励信号和超声换能器激励信号；
组织激励模块200，用于所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；
脉冲回波检测模块300，用于所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；所述脉冲回波检测模块与所述组织激励模块同时启动；
剪切波速度计算模块400，用于根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；
弹性参数获取模块500，用于根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。
具体实施时，所述剪切波速度计算模块400及弹性参数获取模块500设置在电脑5-2中。
进一步的，在所述超声粘弹性测量系统中，所述组织激励模块200具体包括：
功率放大器2-1，用于将所述原始组织激励信号进行放大，得到组织激励信号；
微型激励器3-1，用于接收所述组织激励信号并驱动振动杆3-2；
振动杆3-2，用于在微型激励器3-1的驱动下对所述组织施加指定波形的振动，使组织内产生剪切波。
进一步的，在所述超声粘弹性测量系统中，所述超声换能器激励信号包括由任意信号发生器1-1产生的第一超声换能驱动子信号，和与第一超声换能驱动子信号相位相差180度的第二超声换能驱动子信号。
进一步的，在所述超声粘弹性测量系统中，所述脉冲回波检测模块300具体包括：
耦合单元，用于将所述第一超声换能驱动子信号和所述第二超声换能驱动子信号进行耦合，得到超声换能器驱动信号；
原始超声波发送单元，用于接收所述超声换能器驱动信号，并发出原始超声波p(t)；其中，p(t)＝sin(2πfot)，f0是超声换能器的发射频率，t是时间参数；
超声波接收单元，用于接收所述原始超声波经所述组织的反射和散射后而得到原始回波信号r(t，k)；其中，r(t，k)＝p(t-Δtk)，Δtk表示组织振动产生的延时，t是时间参数，k取任意正整数，θ是原始超声波的超声波束与组织振动方向的夹角，c是超声波声速；
高频滤波单元，用于将所述原始回波信号r(t，k)滤波滤掉高频后，得到回波信号H(t，k)；其中，
采集单元，用于接收所述回波信号H(t，k)；
其中，所述耦合单元、及所述高频滤波单元均设置在超声收/发系统4-1中，所述原始超声波发送单元、及所述超声波接收单元均设置在超声换能器3-3中，所述采集单元为采集卡5-1。
具体实施时，本发明所述超声粘弹性测量系统由任意信号发生器1-1、功率放大器2-1，微型激励器3-1，振动杆3-2，超声换能器3-3，仿体3-4，超声收/发系统4-1，采集卡5-1，电脑5-2组成。其中任意信号发生器1-1在外部时钟7的控制下产生系统内所有需要的信号：原始组织激励信号1，组织激励信号2，第一超声换能驱动子信号3，第二超声换能驱动子信号4，超声换能器驱动信号5，采集卡的触发信号6，原始回波信号8，回波信号9。
由功率放大器2-1放大任意信号发生器1-1产生的原始组织激励信号1，得到组织激励信号2；由微型激励器3-1接收组织激励信号2，并通过带动振动杆3-2、超声换能器3-3作用于仿体3-4；由超声收/发系统4-1耦合任意信号发生器1-1产生的第一超声换能驱动子信号3，第二超声换能驱动子信号4，，将耦合后的超声换能器驱动信号5作用于超声换能器3-3并将其原始回波信号8放大和滤波；采集卡5-1通过采集经过超声收/发系统4-1处理过的回波信号9，处理之后在计算机5-2上存储、处理、显示。
其中：任意信号发生器1-1选用美国的LeCroy Arbstudio 1104，采用外部时钟控制方式，在模拟口产生特定形式的组织激励信号，在数字口产生采集卡的触发信号、相位相差180度的第一超声换能驱动子信号3及第二超声换能驱动子信号4，该第一超声换能驱动子信号3及第二超声换能驱动子信号4可以通过编码进行调制得到超声换能器驱动信号5；
功率放大器2-1和微型激励器3-1是选用丹麦Brüel&Denmark的PowerAmplifier Type 2718和Mini-shaker Type 4810，配套使用，理论上微型激励器可以产生频率低于18kHz的任意波形的振动信号；超声收/发系统4-1和超声换能器3-3是采用德力凯TCD EMS-9A系列的超声收/发模块和配套的超声换能器，其中心频率是2MHz；采集卡5-1是采用ADLINK的PCI-9820，它是一个2通道14位65MS/s采样数字化仪，同时可实现外部触发和外部时钟控制。
在超声收/发系统4-1提供的时钟信号控制下，意信号发生器1-1产生系统所有的信号，且均由该外部时钟7分频得到，这保证了这些信号具有公共的周期，从而可以实现信号之间的严格同步。该外部时钟信号为8MHz，选择超声收/发系统4-1的脉冲重复频率8e6/1080≈7.4kHz，f＝8e6/69120≈115Hz，用f，2f或者更高的f的倍频激励组织使组织产生剪切波，并用超声换能器3-3得到传统检测的结果和编码检测的结果。
本发明的实施例中，设置组织激励模块200、脉冲回波检测模块300和采集卡5-1的时钟源为同一个外部时钟，组织激励模块200、脉冲回波检测模块300和采集卡5-1的触发信号统一由任意信号发生器1-1产生，且这三个触发信号的频率均由外部时钟信号分频得到，保证了这些信号具有公共的周期，从而可以实现信号之间的严格同步。
将本发明所述的超声粘弹性测量系统对标准弹性仿体进行实验时，将标准弹性仿体Elasticity QA Phantom Model 049A放置在超声换能器3-3下方，在超声换能器3-3和仿体3-4之间涂上一层医用耦合剂，用频率f＝8e6/69120≈115Hz去激励组织，脉冲重复频率为8e6/1080≈7.4kHz，采集卡5-1的采样率为8MHz，用传统4个波的脉冲方法获得回波信号，对回波数据进行上述方法处理得到标准弹性仿体的剪切波速度值是2.7m/s，从而得到弹性值是7.3kPa。而后，在相同位置，用7位的Barker码，每个位中有4个波的方式调制超声换能器驱动信号之后接收回波信号，计算得到标准弹性仿体的剪切波速度为2.74m/s，从而得到弹性值是7.5kPa。标准弹性仿体的弹性标称值是7-9.7kPa，编码与传统检测都能在标准弹性仿体上得到在标称值之内的值，且在相同条件下，二者得到的结果很接近。应当理解的是本实施例中的编码是为了解释本发明方法和装置，本发明系统不局限于只有该编码，所有的相位编码均可以用于本发明系统中。
综上所述，本发明提供了一种超声粘弹性测量方法及系统，方法包括：任意信号发生器同时产生原始组织激励信号和超声换能器激励信号；所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。本发明将编码检测技术应用到超声瞬时弹性成像系统，以增加检测深度。
可以理解的是，对本领域普通技术人员来说，可以根据本发明的技术方案及本发明构思加以等同替换或改变，而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。

Claims (10)

1.一种超声粘弹性测量方法，其特征在于，所述方法包括以下步骤：
A、任意信号发生器同时产生原始组织激励信号和超声换能器激励信号；
B、所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，组织内产生剪切波；
C、所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；所述步骤C与所述步骤B同步进行；
D、根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；
E、根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。
2.根据权利要求1所述超声粘弹性测量方法，其特征在于，所述步骤B具体包括：
B1、所述原始组织激励信号通过功率放大器放大，得到组织激励信号；
B2、所述组织激励信号加载在微型激励器上，通过微型激励器驱动振动杆对所述组织施加指定波形的振动，使组织内产生剪切波；其中，深度为d0的剪切波d(t)＝Asin(2πfvt+φv)+d0；其中，A是振动幅值，fv是振动频率，φv是起始相位，t是时间参数。
3.根据权利要求1所述超声粘弹性测量方法，其特征在于，所述超声换能器激励信号包括由任意信号发生器产生的第一超声换能驱动子信号，和与第一超声换能驱动子信号相位相差180度的第二超声换能驱动子信号。
4.根据权利要求3所述超声粘弹性测量方法，其特征在于，所述步骤C具体包括：
C1、所述第一超声换能驱动子信号和所述第二超声换能驱动子信号经超声收/发系统耦合后，得到超声换能器驱动信号；
C2、超声换能器接收所述超声换能器驱动信号，并发出原始超声波p(t)；其中，p(t)＝sin(2πfot)，f0是超声换能器的发射频率，t是时间参数；
C3、所述原始超声波经所述组织的反射和散射后，得到原始回波信号r(t，k)，并被所述超声换能器接收；其中，r(t，k)＝p(t-Δtk)，Δtk表示组织振动产生的延时，t是时间参数，k取任意正整数，θ是原始超声波的超声波束与组织振动方向的夹角，c是超声波声速；
C4、所述原始回波信号r(t，k)经超声收/发系统放大滤波滤掉高频后，得到回波信号H(t，k)；其中，
C5、采集卡接收所述回波信号H(t，k)。
5.根据权利要求1所述超声粘弹性测量方法，其特征在于，所述步骤D中脉冲回波处理算法为Voigt模型，且根据Voigt模型和回波信号得到剪切波在组织中的传播速度cs；其中，剪切波在组织中的频率为ωs，ρ是组织密度，μ1是弹性模量，μ2是粘性系数。
6.根据权利要求1所述超声粘弹性测量方法，其特征在于，所述步骤E中组织的粘弹性数据模型为其中E是弹性模量，G是杨氏模量，σ是泊松比，ρ是组织密度。
7.一种超声粘弹性测量系统，其特征在于，包括：
任意信号发生器，用于同时产生原始组织激励信号和超声换能器激励信号；
组织激励模块，用于所述原始组织激励信号经放大后加载在微型激励器上，并通过振动杆对指定的组织施加指定波形的振动，使组织内产生剪切波；
脉冲回波检测模块，用于所述超声换能器激励信号经超声收/发系统发出以驱动超声换能器，超声换能器发射出的超声波经所述组织反射和散射得到的原始回波信号并由超声换能器接收，并经超声收/发系统传输至采集卡；所述脉冲回波检测模块与所述组织激励模块同时启动；
剪切波速度计算模块，用于根据脉冲回波处理算法，及所述回波信号，得到剪切波在组织中的传播速度；
弹性参数获取模块，用于根据组织的粘弹性数据模型及剪切波在组织中的传播速度，得到组织的弹性数据和粘性数据。
8.根据权利要求7所述超声粘弹性测量系统，其特征在于，所述组织激励模块具体包括：
功率放大器，用于将所述原始组织激励信号进行放大，得到组织激励信号；
微型激励器，用于接收所述组织激励信号并驱动振动杆；
振动杆，用于在微型激励器的驱动下对所述组织施加指定波形的振动，使组织内产生剪切波。
9.根据权利要求7所述超声粘弹性测量系统，其特征在于，所述超声换能器激励信号包括由任意信号发生器产生的第一超声换能驱动子信号，和与第一超声换能驱动子信号相位相差180度的第二超声换能驱动子信号。
10.根据权利要求9所述超声粘弹性测量系统，其特征在于，所述脉冲回波检测模块具体包括：
耦合单元，用于将所述第一超声换能驱动子信号和所述第二超声换能驱动子信号进行耦合，得到超声换能器驱动信号；
原始超声波发送单元，用于接收所述超声换能器驱动信号，并发出原始超声波p(t)；其中，p(t)＝sin(2πfot)，f0是超声换能器的发射频率，t是时间参数；
超声波接收单元，用于接收所述原始超声波经所述组织的反射和散射后而得到原始回波信号r(t，k)；其中，r(t，k)＝p(t-Δtk)，Δtk表示组织振动产生的延时，t是时间参数，k取任意正整数，θ是原始超声波的超声波束与组织振动方向的夹角，c是超声波声速；
高频滤波单元，用于将所述原始回波信号r(t，k)滤波滤掉高频后，得到回波信号H(t，k)；其中，
采集单元，用于接收所述回波信号H(t，k)；
其中，所述耦合单元、及所述高频滤波单元均设置在超声收/发系统中，所述原始超声波发送单元、及所述超声波接收单元均设置在超声换能器中，所述采集单元为采集卡。
CN201610334175.5A 2016-05-19 2016-05-19 一种超声粘弹性测量方法及系统 Pending CN106037816A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610334175.5A CN106037816A (zh) 2016-05-19 2016-05-19 一种超声粘弹性测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610334175.5A CN106037816A (zh) 2016-05-19 2016-05-19 一种超声粘弹性测量方法及系统

Publications (1)

Publication Number Publication Date
CN106037816A true CN106037816A (zh) 2016-10-26

Family

ID=57177167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610334175.5A Pending CN106037816A (zh) 2016-05-19 2016-05-19 一种超声粘弹性测量方法及系统

Country Status (1)

Country Link
CN (1) CN106037816A (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106840362A (zh) * 2017-03-20 2017-06-13 西安交通大学 基于声辐射力脉冲响应的激光测振监测hifu损伤粘弹性的方法
CN107440740A (zh) * 2017-07-21 2017-12-08 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
CN107505393A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 介质粘弹性的测量方法和装置
CN108459082A (zh) * 2018-03-27 2018-08-28 润电能源科学技术有限公司 一种材料检测方法及材料检测装置
CN108685596A (zh) * 2017-04-06 2018-10-23 美国西门子医疗解决公司 利用超声医疗成像的组织性质估计
CN108872903A (zh) * 2017-05-10 2018-11-23 冯原 头部磁共振弹性成像检测方法及成像驱动装置
CN108982674A (zh) * 2017-06-05 2018-12-11 郑翊 使用单个激励脉冲和具有不同宽度和形状的多个激励脉冲测量组织剪切波性质
CN109199381A (zh) * 2018-09-11 2019-01-15 合肥工业大学 一种全息微波弹性成像系统及其成像方法
CN110251159A (zh) * 2019-05-31 2019-09-20 深圳先进技术研究院 一种组织检测方法、装置和存储介质
CN111886496A (zh) * 2017-12-01 2020-11-03 海德拉运动有限公司 通过振动波的反射来测量流体性质
US20210038195A1 (en) * 2018-04-28 2021-02-11 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Ultrasound transient elasticity measurement device and method
WO2022141631A1 (zh) * 2021-01-04 2022-07-07 深圳迈瑞生物医疗电子股份有限公司 粘弹性测量方法和超声成像系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202051729U (zh) * 2011-03-01 2011-11-30 深圳市一体医疗科技股份有限公司 一种用于测量粘弹性介质弹性的测量探头、系统
CN103054552A (zh) * 2012-12-24 2013-04-24 深圳先进技术研究院 生物组织粘弹性测量方法和系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202051729U (zh) * 2011-03-01 2011-11-30 深圳市一体医疗科技股份有限公司 一种用于测量粘弹性介质弹性的测量探头、系统
CN103054552A (zh) * 2012-12-24 2013-04-24 深圳先进技术研究院 生物组织粘弹性测量方法和系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YAN-RONG GUO,ET AL.: "In-Vitro Quantification of Rat Liver Viscoelasticity with Shear Wave Dispersion Ultrasound Vibrometry", 《35TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE EMBS》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106840362A (zh) * 2017-03-20 2017-06-13 西安交通大学 基于声辐射力脉冲响应的激光测振监测hifu损伤粘弹性的方法
CN106840362B (zh) * 2017-03-20 2019-08-23 西安交通大学 基于声辐射力脉冲响应的激光测振监测hifu损伤粘弹性方法
CN108685596A (zh) * 2017-04-06 2018-10-23 美国西门子医疗解决公司 利用超声医疗成像的组织性质估计
CN108872903A (zh) * 2017-05-10 2018-11-23 冯原 头部磁共振弹性成像检测方法及成像驱动装置
CN108982674A (zh) * 2017-06-05 2018-12-11 郑翊 使用单个激励脉冲和具有不同宽度和形状的多个激励脉冲测量组织剪切波性质
CN107505393A (zh) * 2017-07-21 2017-12-22 无锡海斯凯尔医学技术有限公司 介质粘弹性的测量方法和装置
CN107440740A (zh) * 2017-07-21 2017-12-08 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
US11719613B2 (en) 2017-07-21 2023-08-08 Wuxi Hisky Medical Technologies Co., Ltd. Method and device for quantifying viscoelasticity of a medium
CN111886496B (zh) * 2017-12-01 2024-04-30 海德拉运动有限公司 通过振动波的反射来测量流体性质
CN111886496A (zh) * 2017-12-01 2020-11-03 海德拉运动有限公司 通过振动波的反射来测量流体性质
CN108459082A (zh) * 2018-03-27 2018-08-28 润电能源科学技术有限公司 一种材料检测方法及材料检测装置
CN108459082B (zh) * 2018-03-27 2021-01-05 润电能源科学技术有限公司 一种材料检测方法及材料检测装置
US20210038195A1 (en) * 2018-04-28 2021-02-11 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Ultrasound transient elasticity measurement device and method
CN109199381B (zh) * 2018-09-11 2021-11-02 合肥工业大学 一种全息微波弹性成像系统及其成像方法
CN109199381A (zh) * 2018-09-11 2019-01-15 合肥工业大学 一种全息微波弹性成像系统及其成像方法
CN110251159A (zh) * 2019-05-31 2019-09-20 深圳先进技术研究院 一种组织检测方法、装置和存储介质
WO2022141631A1 (zh) * 2021-01-04 2022-07-07 深圳迈瑞生物医疗电子股份有限公司 粘弹性测量方法和超声成像系统

Similar Documents

Publication Publication Date Title
CN106037816A (zh) 一种超声粘弹性测量方法及系统
CN103026257B (zh) 使用剪切波成像的方法和装置
Palmeri et al. A finite-element method model of soft tissue response to impulsive acoustic radiation force
Amador et al. Shear elastic modulus estimation from indentation and SDUV on gelatin phantoms
Gallot et al. Passive elastography: shear-wave tomography from physiological-noise correlation in soft tissues
JP4405920B2 (ja) 人間または動物組織の弾性測定装置および方法と、この弾性の二次元または三次元表示の設定
Viola et al. Radiation force imaging of viscoelastic properties with reduced artifacts
IL167172A (en) Imaging method and device employing shear waves
CN102724917A (zh) 通过传播剪切波来测量哺乳动物软组织物理参数的方法和装置
CN107015230A (zh) 一种超声波测距方法
CN107405133A (zh) 适用于各向异性介质成像的剪切波弹性成像方法及装置
Orescanin et al. 3-D FDTD simulation of shear waves for evaluation of complex modulus imaging
CN112367918A (zh) 用于弹性成像和粘弹性成像的系统和方法
CN103054552A (zh) 生物组织粘弹性测量方法和系统
JP4221555B2 (ja) 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法
EP2622379B1 (en) Device and method for imaging of non-linear and linear properties of formations surrounding a borehole
Zheng et al. Shear Wave Propagation in Soft Tissue with Ultrasound Vibrometry
JP4389091B2 (ja) 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法
Yang et al. A simple viscoelastic model for soft tissues the frequency range 6-20 MHz
Emmer et al. Radial modulation of single microbubbles
Hossain et al. On the quantitative potential of Viscoelastic Response (VisR) ultrasound using matrix array transducers: In silico demonstration
JP5435455B2 (ja) 集束型加振装置
Zhao et al. Elasticity and viscosity estimation from shear wave velocity and attenuation: a simulation study
WO2022213949A1 (zh) 弹性成像方法、装置、电子设备及存储介质
Ganor et al. Time and pressure dependence of acoustic signals radiated from microbubbles

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20161026