CN111637962A - 一种剪切波衰减系数测量方法与*** - Google Patents

一种剪切波衰减系数测量方法与*** Download PDF

Info

Publication number
CN111637962A
CN111637962A CN202010508029.6A CN202010508029A CN111637962A CN 111637962 A CN111637962 A CN 111637962A CN 202010508029 A CN202010508029 A CN 202010508029A CN 111637962 A CN111637962 A CN 111637962A
Authority
CN
China
Prior art keywords
magnetic resonance
nuclear magnetic
vibration
attenuation coefficient
radio frequency
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
CN202010508029.6A
Other languages
English (en)
Other versions
CN111637962B (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.)
Wuxi Marvel Stone Healthcare Co Ltd
Original Assignee
Wuxi Marvel Stone Healthcare Co Ltd
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 Wuxi Marvel Stone Healthcare Co Ltd filed Critical Wuxi Marvel Stone Healthcare Co Ltd
Priority to CN202010508029.6A priority Critical patent/CN111637962B/zh
Publication of CN111637962A publication Critical patent/CN111637962A/zh
Priority to PCT/CN2020/130245 priority patent/WO2021243975A1/zh
Application granted granted Critical
Publication of CN111637962B publication Critical patent/CN111637962B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
    • G01H11/02—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by magnetic means, e.g. reluctance

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施例涉及一种剪切波衰减系数测量方法与***，解决了超声检测法不能穿透某些介质导致测量失败的问题；剪切波衰减系数测量方法包括：包括：向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，分析处理所述核磁共振回波信号得到剪切波衰减系数；所述***包括磁共振***、机械振动激励装置和核磁共振控制台；本发明实施例基于低场核磁共振以非侵入式、非破坏式的方式测量固体或者半固体的衰减系数。

Description

一种剪切波衰减系数测量方法与***
技术领域
本发明实施例涉及一种剪切波衰减系数测量方法与***。
背景技术
衰减系数又称衰减常数。是传播系数的实数部分。它包括两部分：经典吸收和分子吸收。经典吸收是由于空气的粘滞性、热传导效应以及空气分子转动等所产生的声能耗散，其大小与声波频率的平方成正比例，并且与空气温度和气压有关，这种吸收一般可以不考虑。分子吸收主要是空气中氧和氮分子振动弛豫效应所引起的，它与空气的温度和湿度密切相关，也随声波频率的增减而变化，但变化规律较为复杂。
衰减系数描述了波在物质中传播的衰减速率，是反映材料吸收性能的重要指标。在医学上，衰减系数可以反映组织的弹性、脂肪含量、含水量的变化等与生理病理状态密切相关的参数指标。
对于衰减系数的计算往往需要利用剪切波的振幅数据来进行。剪切波在介质中的衰减可以用如下公式描述
Figure BDA0002527250080000011
其中，A0为剪切波的初始振幅，r为被测物体所在位置与振源的距离，α为衰减系数。
现有衰减系数的测量技术大多基于该衰减模型进行改进。该测量技术需要在介质内两点采集振幅，再根据式(1)算出衰减系数：
Figure BDA0002527250080000012
其中，A1，A2分别为波在被测物体内距振源不同距离两点的振幅，Δr为两点间距离，a为衰减系数。
剪切波的振幅数据需要经过测量得到。根据采集剪切波振幅的方式分类，常见的有驻波法和超声检测法。
驻波法为最直接的测量方法，适用于多种机械波。如图1所示，由发生器产生在固定长度介质内传播的驻波，其在介质内损耗大致遵循A＝A0(1+R)e-αr，其中，R为反射系数，α是介质的衰减系数，A和A0分别为距振源r处与振源处振幅。因为发生器和接收器为同一材料制成，输出电压(U)和接收电压(U0)与振幅成正比。可得
Figure BDA0002527250080000021
上述方法需将截取物质装盛在测量设备容器里，或将发生/接收器深入介质内测量传播信号振幅大小，其适用于气体和液体。
超声检测法是基于超声成像的测量方法；适用于固体内衰减系数的测量，即在施加一个外部激励振动，在组织内产生频率为ω的剪切波，并在波传播方向上选择两个相距Δr的点，分别距激励振源r1，r2。在选定两点用超声探头以同样的设定帧率激发采集到两条质点随时间振动曲线，后分别从两条振动曲线中利用卡尔曼滤波抽取出振动频率ω的信号，所抽取出信号幅度即为对应位置的振幅。
驻波法适用于液体和气体，测量时需要将探测器浸入物质内部或截取物质样本，适用范围较窄；超声检测衰减系数具有低成本、快速定量等优点，但超声波穿透性易受介质影响，容易导致测量失败。
发明内容
本发明实施例提供一种剪切波衰减系数测量方法与***，解决了超声检测法不能穿透某些介质导致测量失败的问题。
第一方面，一种剪切波衰减系数测量方法，包括：向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，分析处理所述核磁共振回波信号得到剪切波衰减系数。
结合第一方面，在第一种可能的实现方式中，所述向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，包括：向作简谐振动的、处于梯度磁场中的被测物体施加90°射频脉冲，根据振动源的频率确定至少一个180°射频脉冲的施加时刻并于该时刻向所述振动状态稳定的被测物体施加所述至少一个180°射频脉冲。
结合第一方面的第一种可能的实现方式，在第二种可能的实现方式中，所述测量所述被测物体至少两个不同位置的核磁共振回波信号，包括检测所述被测物体不同中心频率的核磁共振回波信号；
或包括检测磁体和/或探头与所述被测物体的不同位置处的核磁共振回波信号。
结合第一方面的第二种可能的实现方式，在第三种可能的实现方式中，所述测量所述被测物体至少两个不同位置的核磁共振回波信号，包括对每一个位置等间隔的调整振动源的初始相位用以获取该位置一系列的核磁共振回波信号。
结合第一方面或第一方面的第三种可能的实现方式，在第四种可能的实现方式中，所述对每一个位置等间隔的调整振动源的初始相位，包括等间隔的调整所述振动源的初始相位，所述间隔大于2π。
结合第一方面或第一方面的上述任一种可能的实现方式，在第五种可能的实施方式中，所述分析处理所述核磁共振回波信号得到剪切波衰减系数，包括：
设S1(θ,n)为在第一位置采集到的核磁共振回波信号，S2(θ,n)为在第二位置采集到的核磁共振回波信号，θ表示振动源的初始相位，第二维n表示回波的点数；
分别对S1(θ,n)和S2(θ,n)的第二维做一维傅里叶变换，取最低频部分的相位，分别得到
Figure BDA0002527250080000041
和
Figure BDA0002527250080000042
分别对
Figure BDA0002527250080000043
和
Figure BDA0002527250080000044
做一维相位反卷折运算；
分别计算
Figure BDA0002527250080000045
和
Figure BDA0002527250080000046
即分别计算
Figure BDA0002527250080000047
和
Figure BDA0002527250080000048
的峰值，得到
Figure BDA0002527250080000049
Figure BDA00025272500800000410
计算第一位置和第二位置处剪切波振幅
Figure BDA00025272500800000411
Figure BDA00025272500800000412
式(5)和(6)中，G为运动敏感梯度，γ为拉莫尔频率，T为剪切波振动周期，N为在整个运动编码持续时间内剪切波的周期数目；
按以下公式(7)运算出剪切波在组织内传播的衰减系数值
Figure BDA00025272500800000413
式(7)中α是衰减系数，r1是第一位置距振源距离，r2是第二位置距振源距离。
第二方面，一种用于实现所述的剪切波衰减系数测量方法的测量***，包括
磁共振***，用于按照核磁共振脉冲序列指令发射射频脉冲，定位被测物体空间位置，用于接收核磁共振回波信号；
机械振动激励装置，用于接收所述磁共振***发出的射频脉冲信号并根据所述信号使被测物体产生简谐振动；
核磁共振控制台，控制磁共振***运行核磁共振脉冲序列指令，接收磁共振***采集到的核磁共振回波信号，对所述核磁共振回波信号分析处理。
结合第二方面，在第一种可能的实现方式中，所述磁共振***包括磁共振谱仪、射频功放、前置放大器、收发切换模块、磁体模块以及射频探头；
磁共振谱仪，用于向机械振动激励装置发送触发信号，所述磁共振谱仪通过射频功放与收发切换模块连接；
收发切换模块，用于切换所述磁共振***的发射状态和接收状态；
射频探头，所述射频探头与收发切换模块连接；发射状态下，所述射频探头用于向被测物体发射射频脉冲；接收状态下，所述射频探头用于接收被测物体检测靶位受激发后产生的核磁共振回波信号；
磁体模块，用于产生静态梯度磁场对被测物体进行空间定位。
结合第二方面的第一种可能的实现方式，在第二种可能的实现方式中，所述机械振动激励装置包括信号发生器、功率放大器和振动发生器；
信号发生器，用于接收磁共振谱仪发出的触发信号并根据所述触发信号驱动振动发生器产生简谐振动；
振动发生器，所述振动发生器与传动杆连接；所述振动发生器通过传动杆作用于被测物体使所述被测物体产生简谐振动。
结合第二方面的第二种可能的实现方式，在第三种可能的实现方式中，收发切换模块通过前置放大器与磁共振谱仪连接；所述振动发生器通过功率放大器与信号发生器连接。
本发明实施例的一种剪切波衰减系数测量方法与***，基于磁共振原理以非侵入式，非破坏式的方式测量固体或者半固体的衰减系数；显然，避免了超声检测法不能穿透某些介质导致测量失败的问题。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面对实施例描述中所需要使用的附图作简要介绍。
图1是驻波法测量衰减系数的原理示意图。
图2是本发明实施例的测量方法流程示意图。
图3是本发明实施例的***结构示意图。
图4是本发明实施例增加梯度***的结构示意图。
图5是本发明实施例的磁体模块磁场分布示意图。
图6是本发明实施例的核磁共振脉冲序列示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施例作进一步地详细描述。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明实施例提供一种剪切波衰减系数测量方法与***，用以解决现有技术中超声检测法不能穿透某些介质导致测量失败的问题。其中，方法和***是基于同一发明构思的，由于方法及装置解决问题的原理相似，因此装置与方法的实施可以相互参见，重复之处不再赘述。
本发明中所涉及的多个，是指两个或两个以上。另外，在本发明的描述中，“第一”、“第二”等词汇，仅用于区分描述的目的，而不能理解为指示或暗示相对重要性，也不能理解为指示或暗示顺序。
本发明实施例提供了一种剪切波衰减系数测量方法与***，该方法与***基于低场核磁共振原理实现，以非侵入式，非破坏式的方式测量固体或者半固体的衰减系数。本发明实施例的测量方法还具有低成本、测量时间短、定位准确、测量结果可重复性高的优点。具有既可以用于工业检测，也可以用于医学检测的特点。下面结合附图对本发明优选的实施方式进行详细说明。
参阅图2所示，本发明实施例提供一种剪切波衰减系数测量方法，包括：向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，分析处理所述核磁共振回波信号得到剪切波衰减系数。
通过对被测物体施加核磁共振脉冲序列，对被测物体进行运动编码，定位被测物体的空间位置，采集至少两个不同位置的磁共振信号，最终由磁共振信号的相位信息和所测点与振动源的距离等信息直接计算出衰减系数。
可选地，所述向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，包括：向作简谐振动的、处于梯度磁场中的被测物体施加90°射频脉冲，根据振动源的频率确定至少一个180°射频脉冲的施加时刻并于该时刻向所述振动状态稳定的被测物体施加所述至少一个180°射频脉冲。
如图5所示，在振动状态稳定后，施加一个90°射频脉冲，由于磁场梯度的存在，ROI范围内一定厚度内的自旋将被激发，根据振动源的频率可以确定180°脉冲的施加时刻。在180°脉冲之后，自旋聚相，在物质内不同两点检测回波信号并传回控制台完成计算。在实际中，可以通过部分运动编码缩短回波时间或多次信号累加来提高信噪比。此外，为了提高相位信噪比，可以施加多个180°脉冲。
定位被测物体的空间位置，可以使用梯度磁场作用与被测物体的方式来实现。梯度磁场主要用于空间定位，包括相位编码及频率编码，可以通过梯度磁场明确空间上的任意位置。本发明实施例的***采用了特殊磁场分布的单边磁体，一般地，采用传统的低场磁共振***或其他形式的磁体设计方案也能达到相同的目的。
可选地，测量所述被测物体至少两个不同位置的核磁共振回波信号，包括检测所述被测物体不同中心频率的核磁共振回波信号；或包括检测磁体和/或探头与所述被测物体的不同位置处的核磁共振回波信号。
本发明实施例需要测量两个位置的磁共振信号，可以通过改变中心频率的方式或者改变磁体、探头与被测物体的相对位置的方式实现。
可选地，所述测量所述被测物体至少两个不同位置的核磁共振回波信号，包括对每一个位置等间隔的调整振动源的初始相位用以获取该位置一系列的核磁共振回波信号。
优选地，对每个位置采集一组核磁共振回波信号；对每个位置的一组核磁共振回波信号的采集，可以通过对每一个位置等间隔的调整振动源的初始相位来获得，所述等间隔指的是各个振动状态下各个初始相位的差值固定；优选地，各个初始相位的等间隔的值大于2π。
可选地，所述分析处理所述核磁共振回波信号得到剪切波衰减系数，包括：
设S1(θ,n)为在第一位置采集到的核磁共振回波信号，S2(θ,n)为在第二位置采集到的核磁共振回波信号，θ表示振动源的初始相位，第二维n表示回波的点数；
分别对S1(θ,n)和S2(θ,n)的第二维做一维傅里叶变换，取最低频部分的相位，分别得到
Figure BDA0002527250080000081
和
Figure BDA0002527250080000082
分别对
Figure BDA0002527250080000083
和
Figure BDA0002527250080000084
做一维相位反卷折运算；
分别计算
Figure BDA0002527250080000091
和
Figure BDA0002527250080000092
即分别计算
Figure BDA0002527250080000093
和
Figure BDA0002527250080000094
的峰值，得到
Figure BDA0002527250080000095
Figure BDA0002527250080000096
计算第一位置和第二位置处剪切波振幅
Figure BDA0002527250080000097
Figure BDA0002527250080000098
式(5)和(6)中，G为运动敏感梯度，γ为拉莫尔频率，T为剪切波振动周期，N为在整个运动编码持续时间内剪切波的周期数目；
按以下公式(7)运算出剪切波在组织内传播的衰减系数值
Figure BDA0002527250080000099
式(7)中α是衰减系数，r1是第一位置距振源距离，r2是第二位置距振源距离。
其中，当测量的位置为两个时，第一位置和第二位置与测量的两个位置一一对应；当测量的位置为3个及以上时，第一位置和第二位置为3个及以上的测量位置中的任意两个。
本发明实施例仅阐述采集了两个位置进行衰减系数计算的方法。实际上，可以采集两个以上的位置，用与本发明公开方法类似的方法测量衰减系数，都应属于本发明的保护范围。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤，而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
如图3所示，本发明实施例还提供一种用于实现所述的剪切波衰减系数测量方法的测量***，包括磁共振***，用于按照核磁共振脉冲序列指令发射射频脉冲，定位被测物体空间位置，用于接收核磁共振回波信号；机械振动激励装置，用于接收所述磁共振***发出的射频脉冲信号并根据所述信号使被测物体产生简谐振动；核磁共振控制台，控制磁共振***运行核磁共振脉冲序列指令，接收磁共振***采集到的核磁共振回波信号，对所述核磁共振回波信号分析处理。
可选地，为便于携带所述磁共振谱仪和磁体模块均为便携式。所述磁共振***为便携式低场磁共振***，从而，使整体***更轻巧便捷，成本更低。
核磁共振控制台发送核磁共振脉冲序列指令至机械振动激励装置，机械振动激励装置作用于被测物体使其产生简谐振动，被测物体在磁共振***中通过梯度磁场进行空间位置定位，磁共振***将被测物体的核磁共振回波信号接收后返回至核磁共振控制台；核磁共振控制台通过指令控制机械振动激励装置对每一个位置等间隔的调整振动源的初始相位来获得一系列核磁共振回波信号返回至核磁共振控制台；核磁共振控制台也可以通过指令以及相应的执行机构以改变中心频率的方式或者改变磁体、探头与被测物体的相对位置的方式来产生不同位置的核磁共振回波信号并将其返回至核磁共振控制台；核磁共振控制台接收到不同位置的一系列核磁共振回波信号后对所述核磁共振回波信号分析处理，得到衰减系数。
可选地，所述磁共振***包括磁共振谱仪、射频功放、前置放大器、收发切换模块、磁体模块以及射频探头；磁共振谱仪，用于向机械振动激励装置发送触发信号，所述磁共振谱仪通过射频功放与收发切换模块连接；收发切换模块，用于切换所述磁共振***的发射状态和接收状态；射频探头，所述射频探头与收发切换模块连接；发射状态下，所述射频探头用于向被测物体发射射频脉冲；接收状态下，所述射频探头用于接收被测物体检测靶位受激发后产生的核磁共振回波信号；磁体模块，用于产生静态梯度磁场对被测物体进行空间定位。
优选地，收发切换模块通过前置放大器与磁共振谱仪连接；所述振动发生器通过功率放大器与信号发生器连接。
具体地，如图3所示，磁共振谱仪通过射频功放与收发切换开关连接；磁共振谱仪通过前置放大器与收发切换开关连接；磁共振谱仪通过收发切换开关与射频探头连接；所述收发切换模块用于切换所述便携式磁共振***处于发射状态或接收状态，发射状态下，所述探头发射射频脉冲；接收状态下，所述射频探头用于接收被测物体靶位受激发后产生的核磁共振回波信号。磁共振谱仪与机械振动激励装置通讯；磁共振控制台通过射频功放分别与磁共振谱仪和收发切换开关连接；所述磁体模块设于射频探头后侧；被测物体设于射频探头与机械振动激励装置的传动杆之间；传动杆将机械振动激励装置的振动传给被测物体，使被测物体在磁体模块的梯度磁场中作简谐振动。
优选地，若使用静态磁场，所述磁共振谱仪不设梯度控制模块，因此也无需设置相应的梯度线圈。磁体模块设有单边磁体，背面连接设有磁轭，背面磁场强度迅速衰减，所述单边磁体产生一个静态磁场，所述静态磁场在ROI范围内存在AP-LR平面内的线性或近似线性梯度，在ROI范围之外，磁场均匀地迅速衰减；磁体模块的磁场分布如图6所示。
可选地，若使用可控梯度磁场，如图4所示，所述磁共振***还包括梯度***；所述梯度***包括梯度放大器和梯度线圈；所述梯度放大器与磁共振谱仪通讯，梯度放大器与梯度线圈连接；所述梯度线圈设于射频探头和磁体模块之间，用于对被测物体产生可控梯度磁场。
梯度***由梯度放大器和梯度线圈组成，由磁共振谱仪控制，经由梯度放大器将信号放大后由梯度线圈在被测物体所在空间内形成梯度磁场；本实施方式中除脉冲序列和静态梯度磁场配合进行运动编码外，脉冲序列也可与梯度***配合，利用梯度线圈产生的可控梯度磁场来进行运动编码。
优选地，所述磁共振谱仪设有发射和门控单向信号通路，所述磁共振谱仪通过发射和门控单向信号通路与射频功放连接；所述磁共振谱仪通过收发转换门控与收发切换开关连接，所述射频功放将发射信号放大后与收发切换开关连接，所述收发切换开关与所述射频探头连接。
可选地，机械振动激励装置包括信号发生器、功率放大器和振动发生器；信号发生器，用于接收磁共振谱仪发出的触发信号并根据所述触发信号驱动振动发生器产生简谐振动；振动发生器，所述振动发生器与传动杆连接；所述振动发生器通过传动杆作用于被测物体使所述被测物体产生简谐振动。
优选地，所述信号发生器为波形发生器，波形发生器向功率放大器发送50HZ正弦波进行放大，振动功率放大器与信号发生器连接，信号发生器与所述磁共振谱仪连接，接收磁共振谱仪发出的触发信号由信号发生器传递至振动发生器，振动发生器驱动传动杆产生简谐振动；所述简谐振动器经由传动杆紧贴在被测物质邻近表面上产生振动后在其内部产生剪切波，所述剪切波在AP方向上传播，引起组织内质点在LR方向上做简谐振动，所述核磁共振控制台与所述磁共振谱仪连接，控制运行磁共振脉冲序列指令，并接收磁共振谱仪采集到的核磁共振回波信号，完成实时数据处理。
本领域普通技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明实施例进行各种改动和变型而不脱离本发明实施例的精神和范围。这样，倘若本发明实施例的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (10)

1.一种剪切波衰减系数测量方法，其特征在于：包括：向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，分析处理所述核磁共振回波信号得到剪切波衰减系数。
2.根据权利要求1所述的一种剪切波衰减系数测量方法，其特征在于：所述向振动状态稳定的被测物体施加核磁共振脉冲序列，并配合静态或可控梯度磁场进行运动编码，以检测所述被测物体至少两个不同位置的核磁共振回波信号，
包括：向作简谐振动的、处于梯度磁场中的被测物体施加90°射频脉冲，根据振动源的频率确定至少一个180°射频脉冲的施加时刻并于该时刻向所述振动状态稳定的被测物体施加所述至少一个180°射频脉冲。
3.根据权利要求2所述的一种剪切波衰减系数测量方法，其特征在于：所述测量所述被测物体至少两个不同位置的核磁共振回波信号，包括检测所述被测物体不同中心频率的核磁共振回波信号；
或包括检测磁体和/或探头与所述被测物体的不同位置处的核磁共振回波信号。
4.根据权利要求3所述的一种剪切波衰减系数测量方法，其特征在于：所述测量所述被测物体至少两个不同位置的核磁共振回波信号，包括对每一个位置等间隔的调整振动源的初始相位用以获取该位置一系列的核磁共振回波信号。
5.根据权利要求4所述的一种剪切波衰减系数测量方法，其特征在于：所述对每一个位置等间隔的调整振动源的初始相位，包括等间隔的调整所述振动源的初始相位，所述间隔大于2π。
6.根据权利要求1-5任意一项所述的一种剪切波衰减系数测量方法，其特征在于：所述分析处理所述核磁共振回波信号得到剪切波衰减系数，包括：
设S1(θ,n)为在第一位置采集到的核磁共振回波信号，S2(θ,n)为在第二位置采集到的核磁共振回波信号，θ表示振动源的初始相位，第二维n表示回波的点数；
分别对S1(θ,n)和S2(θ,n)的第二维做一维傅里叶变换，取最低频部分的相位，分别得到
Figure FDA0002527250070000021
和
Figure FDA0002527250070000022
分别对
Figure FDA0002527250070000023
和
Figure FDA0002527250070000024
做一维相位反卷折运算；
分别计算
Figure FDA0002527250070000025
和
Figure FDA0002527250070000026
即分别计算
Figure FDA0002527250070000027
和
Figure FDA0002527250070000028
的峰值，得到
Figure FDA0002527250070000029
Figure FDA00025272500700000210
计算第一位置和第二位置处剪切波振幅
Figure FDA00025272500700000211
Figure FDA00025272500700000212
式(5)和(6)中，G为运动敏感梯度，γ为拉莫尔频率，T为剪切波振动周期，N为在整个运动编码持续时间内剪切波的周期数目；
按以下公式(7)运算出剪切波在组织内传播的衰减系数值
Figure FDA00025272500700000213
式(7)中α是衰减系数，r1是第一位置距振源距离，r2是第二位置距振源距离。
7.一种用于实现权利要求1-6任意一项所述的剪切波衰减系数测量方法的测量***，其特征在于：包括
磁共振***，用于按照核磁共振脉冲序列指令发射射频脉冲，定位被测物体空间位置，用于接收核磁共振回波信号；
机械振动激励装置，用于接收所述磁共振***发出的射频脉冲信号并根据所述信号使被测物体产生简谐振动；
核磁共振控制台，控制磁共振***运行核磁共振脉冲序列指令，接收磁共振***采集到的核磁共振回波信号，对所述核磁共振回波信号分析处理。
8.根据权利要求7所述的测量***，其特征在于：所述磁共振***包括磁共振谱仪、射频功放、前置放大器、收发切换模块、磁体模块以及射频探头；
磁共振谱仪，用于向机械振动激励装置发送触发信号，所述磁共振谱仪通过射频功放与收发切换模块连接；
收发切换模块，用于切换所述磁共振***的发射状态和接收状态；
射频探头，所述射频探头与收发切换模块连接；发射状态下，所述射频探头用于向被测物体发射射频脉冲；接收状态下，所述射频探头用于接收被测物体检测靶位受激发后产生的核磁共振回波信号；
磁体模块，用于产生静态梯度磁场对被测物体进行空间定位。
9.根据权利要求8所述的测量***，其特征在于：所述机械振动激励装置包括信号发生器、功率放大器和振动发生器；
信号发生器，用于接收磁共振谱仪发出的触发信号并根据所述触发信号驱动振动发生器产生简谐振动；
振动发生器，所述振动发生器与传动杆连接；所述振动发生器通过传动杆作用于被测物体使所述被测物体产生简谐振动。
10.根据权利要求9所述的测量***，其特征在于：所述收发切换模块通过前置放大器与磁共振谱仪连接；所述振动发生器通过功率放大器与信号发生器连接。
CN202010508029.6A 2020-06-05 2020-06-05 一种剪切波衰减系数测量方法与*** Active CN111637962B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010508029.6A CN111637962B (zh) 2020-06-05 2020-06-05 一种剪切波衰减系数测量方法与***
PCT/CN2020/130245 WO2021243975A1 (zh) 2020-06-05 2020-11-19 剪切波衰减系数测量方法与***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010508029.6A CN111637962B (zh) 2020-06-05 2020-06-05 一种剪切波衰减系数测量方法与***

Publications (2)

Publication Number Publication Date
CN111637962A true CN111637962A (zh) 2020-09-08
CN111637962B CN111637962B (zh) 2021-04-20

Family

ID=72329812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010508029.6A Active CN111637962B (zh) 2020-06-05 2020-06-05 一种剪切波衰减系数测量方法与***

Country Status (2)

Country Link
CN (1) CN111637962B (zh)
WO (1) WO2021243975A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112244813A (zh) * 2020-10-22 2021-01-22 无锡鸣石峻致医疗科技有限公司 一种低场核磁共振弹性测量方法及***
WO2021243975A1 (zh) * 2020-06-05 2021-12-09 无锡鸣石峻致医疗科技有限公司 剪切波衰减系数测量方法与***

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5810731A (en) * 1995-11-13 1998-09-22 Artann Laboratories Method and apparatus for elasticity imaging using remotely induced shear wave
JP2005245775A (ja) * 2004-03-04 2005-09-15 Toshiba Corp 磁気共鳴イメージング装置
JP2005270583A (ja) * 2004-03-26 2005-10-06 Toshiba Corp 磁気共鳴装置
CN1707290A (zh) * 2004-06-04 2005-12-14 施卢默格海外有限公司 使用脉冲场梯度nmr测量确定流体取样测井工具中流体性质的方法和装置
US20100074501A1 (en) * 2008-09-22 2010-03-25 Siemens Medical Solutions Usa, Inc. Co-Registering Attenuation Data and Emission Data in Combined Magnetic Resonance/Positron Emission Tomography (MR/PET) Imaging Apparatus
CN101732048A (zh) * 2008-11-21 2010-06-16 西门子公司 处理组合的磁共振-发射断层造影图像的方法和装置
US20130285657A1 (en) * 2012-04-25 2013-10-31 Los Alamos National Security, Llc Classification of Materials Using Nuclear Magnetic Resonance Dispersion and/or X-Ray Absorption
KR20130130555A (ko) * 2012-05-22 2013-12-02 서울대학교산학협력단 광자감쇠지도 생성 방법 및 ｍｒｉ 융합영상시스템
DE102012214012A1 (de) * 2012-08-07 2014-02-13 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Identifizierung zumindest eines strahlungsschwächenden Implantats für eine medizinische MR-PET-Bildgebung
CN103792585A (zh) * 2014-01-28 2014-05-14 淮南矿业(集团)有限责任公司 多介质结构测量方法及装置
CN104013388A (zh) * 2014-06-13 2014-09-03 中国医学科学院生物医学工程研究所 基于低频连续波的磁声耦合成像激励与检测方法及装置
CN104142326A (zh) * 2014-06-27 2014-11-12 中国石油化工股份有限公司 一种衰减系数检测方法
CN106290564A (zh) * 2016-10-08 2017-01-04 沈阳工业大学 粮食声衰减系数测量***及方法
CN107656223A (zh) * 2017-09-29 2018-02-02 福建加谱新科科技有限公司 基于傅里叶变换核磁共振谱叠加峰形的核磁共振谱获取方法
CN108375746A (zh) * 2017-01-18 2018-08-07 上海联影医疗科技有限公司 一种相位反卷折方法与设备
CN110161126A (zh) * 2019-06-27 2019-08-23 中国工程物理研究院化工材料研究所 适用于宽温域的固体应力波幅值衰减系数测试装置及方法
CN110726775A (zh) * 2019-11-14 2020-01-24 四川省食品药品检验检测院(四川省药品质量研究所、四川省医疗器械检测中心) 声速与声衰减系数测量装置及方法
CN110780247A (zh) * 2019-11-12 2020-02-11 无锡鸣石峻致医疗科技有限公司 一种基于磁共振原理的器官脂肪无创定量检测方法
CN110916662A (zh) * 2019-12-05 2020-03-27 无锡鸣石峻致医疗科技有限公司 一种便携式核磁共振器官弹性无创定量检测***

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6528995B1 (en) * 2001-09-10 2003-03-04 Schlumberger Technology Corporation Methods and apparatus for measuring flow velocity in a wellbore using NMR and applications using same
WO2016159979A1 (en) * 2015-03-31 2016-10-06 Halliburton Energy Services, Inc. Fiber optic nuclear magnetic resonance sensor
CN107389794B (zh) * 2017-08-14 2020-01-10 中国石油大学(北京) 测量岩石衰减系数的方法与***
CN110927252B (zh) * 2019-10-30 2023-05-12 深圳大学 一种靶向剪切波弹性成像检测方法
CN111637962B (zh) * 2020-06-05 2021-04-20 无锡鸣石峻致医疗科技有限公司 一种剪切波衰减系数测量方法与***

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5810731A (en) * 1995-11-13 1998-09-22 Artann Laboratories Method and apparatus for elasticity imaging using remotely induced shear wave
JP2005245775A (ja) * 2004-03-04 2005-09-15 Toshiba Corp 磁気共鳴イメージング装置
JP2005270583A (ja) * 2004-03-26 2005-10-06 Toshiba Corp 磁気共鳴装置
CN1707290A (zh) * 2004-06-04 2005-12-14 施卢默格海外有限公司 使用脉冲场梯度nmr测量确定流体取样测井工具中流体性质的方法和装置
US20100074501A1 (en) * 2008-09-22 2010-03-25 Siemens Medical Solutions Usa, Inc. Co-Registering Attenuation Data and Emission Data in Combined Magnetic Resonance/Positron Emission Tomography (MR/PET) Imaging Apparatus
CN101732048A (zh) * 2008-11-21 2010-06-16 西门子公司 处理组合的磁共振-发射断层造影图像的方法和装置
US20130285657A1 (en) * 2012-04-25 2013-10-31 Los Alamos National Security, Llc Classification of Materials Using Nuclear Magnetic Resonance Dispersion and/or X-Ray Absorption
KR20130130555A (ko) * 2012-05-22 2013-12-02 서울대학교산학협력단 광자감쇠지도 생성 방법 및 ｍｒｉ 융합영상시스템
DE102012214012A1 (de) * 2012-08-07 2014-02-13 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Identifizierung zumindest eines strahlungsschwächenden Implantats für eine medizinische MR-PET-Bildgebung
CN103792585A (zh) * 2014-01-28 2014-05-14 淮南矿业(集团)有限责任公司 多介质结构测量方法及装置
CN104013388A (zh) * 2014-06-13 2014-09-03 中国医学科学院生物医学工程研究所 基于低频连续波的磁声耦合成像激励与检测方法及装置
CN104142326A (zh) * 2014-06-27 2014-11-12 中国石油化工股份有限公司 一种衰减系数检测方法
CN106290564A (zh) * 2016-10-08 2017-01-04 沈阳工业大学 粮食声衰减系数测量***及方法
CN108375746A (zh) * 2017-01-18 2018-08-07 上海联影医疗科技有限公司 一种相位反卷折方法与设备
CN107656223A (zh) * 2017-09-29 2018-02-02 福建加谱新科科技有限公司 基于傅里叶变换核磁共振谱叠加峰形的核磁共振谱获取方法
CN110161126A (zh) * 2019-06-27 2019-08-23 中国工程物理研究院化工材料研究所 适用于宽温域的固体应力波幅值衰减系数测试装置及方法
CN110780247A (zh) * 2019-11-12 2020-02-11 无锡鸣石峻致医疗科技有限公司 一种基于磁共振原理的器官脂肪无创定量检测方法
CN110726775A (zh) * 2019-11-14 2020-01-24 四川省食品药品检验检测院(四川省药品质量研究所、四川省医疗器械检测中心) 声速与声衰减系数测量装置及方法
CN110916662A (zh) * 2019-12-05 2020-03-27 无锡鸣石峻致医疗科技有限公司 一种便携式核磁共振器官弹性无创定量检测***

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
KEITH D PAULSEN 等: ""Magnetic resonance elastography in the breast: initial reconstructions of damping coefficient"", 《IEEE》 *
冯若: "《超声手册》", 31 December 1999 *
吴水才: "《医学超声及应用》", 31 December 2012 *
尹昌 等: ""钢材料中声速与声衰减系数的测量"", 《声学技术》 *
张钊骞 等: ""衰减系数测量技术的发展和展望"", 《无损探伤》 *
李传辉: ""含水合物地层地震波衰减特性及衰减系数估计方法研究"", 《中优秀博士学位论文全文数据库》 *
李水 等: ""水声材料横波速度和衰减系数参量源法测量***"", 《声学学报》 *
牛金海: "《超声原理及生物医学工程应用》", 31 December 2017 *
程昌林: ""固体横波衰减系数的测定"", 《电子科技大学学报》 *
郝亚平 等: ""超声测定肝脏衰减系数与糖脂代谢的关系"", 《上海医学》 *
雷斯尼克: "《骨与关节疾病诊断学》", 31 December 2009 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021243975A1 (zh) * 2020-06-05 2021-12-09 无锡鸣石峻致医疗科技有限公司 剪切波衰减系数测量方法与***
CN112244813A (zh) * 2020-10-22 2021-01-22 无锡鸣石峻致医疗科技有限公司 一种低场核磁共振弹性测量方法及***
CN112244813B (zh) * 2020-10-22 2022-08-26 无锡鸣石峻致医疗科技有限公司 一种低场核磁共振弹性测量方法及***

Also Published As

Publication number Publication date
WO2021243975A1 (zh) 2021-12-09
CN111637962B (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
US6246895B1 (en) Imaging of ultrasonic fields with MRI
KR101398948B1 (ko) 진폭―위상 변조된 초음파를 이용한 점탄성 측정
US5099848A (en) Method and apparatus for breast imaging and tumor detection using modal vibration analysis
FI64282B (fi) Diagnosapparatur foer bestaemmande av vaevnadernas struktur oc sammansaettning
US6566873B1 (en) Method of and apparatus for nuclear quadrupole resonance testing a sample
US8508229B2 (en) Shear mode pressure-activated driver for magnetic resonance elastography
US8558547B2 (en) System and method for magnetic resonance radio-frequency field mapping
US9588205B2 (en) Simultaneous and dynamic determination of longitudinal and transversal relaxation times of a nuclear spin system
CN111637962B (zh) 一种剪切波衰减系数测量方法与***
JP2003225223A5 (zh)
CN205338953U (zh) 一种骨密度检测设备
CN103857331A (zh) 用于量化人体中肝脏脂肪的方法和***
CN110916663B (zh) 一种便携式核磁共振器官弹性无创定量检测方法
CN112244813B (zh) 一种低场核磁共振弹性测量方法及***
US5194809A (en) Nmr-scope
Selby et al. Dynamic mechanical analysis with portable NMR
JP3837180B2 (ja) 磁気共鳴で被検体の弾性画像を発生する方法
JP2009025093A (ja) 電磁超音波測定装置および電磁超音波を用いた板厚および応力の測定方法
Aindow et al. Quantitative investigation of disk ultrasonic sources
CN110916662B (zh) 一种便携式核磁共振器官弹性无创定量检测***
WO2006119511A2 (en) Ultrasonic method to determine bone parameters
JP4125134B2 (ja) 磁気共鳴アコーストグラフィ
Battaglini et al. Frequency modulated continuous wave ultrasonic radar
JP4609975B2 (ja) 磁気共鳴イメージング装置
JP2664189B2 (ja) 核磁気共鳴を用いた検査装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant