CN101692976A - 基于微波、射频多普勒效应的血液动力学检测方法及其装置 - Google Patents

Abstract

本发明公开了一种基于微波、射频多普勒效应的血液动力学检测方法及其装置，特点是采用异形开口同轴天线，保证收发端天线的匹配；通过在微波和射频频段优化信号频率，以获得能兼顾微波或射频信号对人体的穿透深度和多普勒信号的分辨率的性能；通过设计超窄带带通滤波器，以抑制设备测试中的强零频干扰，并采用了弱信号处理手段以减低对探头位置的依赖和减少外界干扰的影响，提高血液流速检测的实时性，精确度和方便性；通过应用先进的时频联合分析方法，能够对心脏压缩和扩张时血液最大速度和最小速度进行跟踪，并可提供一个心脏周期中的相关血管中的血流速度变化信息，为医生提供可靠的诊断信息，可有效降低血液检测过程中意外感染的危险性。

Description

基于微波、射频多普勒效应的血液动力学检测方法及其装置

技术领域

本发明涉及一种检测方法及装置，尤其是利用微波、射频多普勒效应，对视力可及或浅表的血管内血液的流速进行实时无损检测的方法及装置。

背景技术

血流速测量技术在医疗上有非常重要的实用价值，可以为医生提供有效的诊断信息，也可在手术中为医生的校正操作提供及时、精确和有效的信息。血流速度测量可极大地减少由于早期技术性问题而引起的术后并发症的发生，比如移植手术中血管缝合的不畅通问题。在医学技术快速发展的时代，这种技术为手术成功提供了必要条件：保证血管的通畅性，同时也可以作为血压测量的一种补充。对医生而言，使用该项技术可以提高对病人的监护标准，为手术后病人的护理方案提供依据，降低医生误诊的心理负担。对医院而言，该项技术除了可以改善医疗质量、减少医疗事故、避免再次手术，减轻患者的痛苦，还可以减少医疗诉讼案及医院承担的责任。血动力检测的信息还可以作为医生对心血管系统疾病诊断的一种潜在的辅助手段。

现有的血流速测量技术主要有三种，分别利用了超声多普勒效应、电磁感应和激光多普勒效应。三项技术分别具有自己的一些优势，但是同时也都存在着一些不足。比如超声多普勒或电磁感应技术对探头的位置，角度非常敏感，可控性差、操作要求比较高。激光多普勒技术则容易受到外界光的干扰，并且不适合测量单个血管的血流速。

目前有如下一些专利涉及到了微波多普勒测量流速的研究：

专利一(美国专利号：4167736)设计了基本的微波多普勒测量所需的硬件，包括天线，微波电路等，属于较早利用该方法测量流速的专利之一。

专利二(美国专利号：4091385)给出了一种校准方法，能够在没有第三方测速手段定标的情况，独立给出有效的速度测量值。

专利三(WO专利号：2001061283A1)能够测量混合流体的平均速度，而且给出了一套算法，能够利用反射功率来计算不同流体所占的百分比，在后续信号处理中引入了数字信号算法。

专利四(美国专利号：5315880)除了利用微波多普勒效应来测量流速外，还利用了一个超声波发生器测量导管的横截面积，从而利用完全无损的方法求得单位时间内流量。

但是，这些专利均没有直接用于血流的测量。专利五(美国专利号：US7263398B2)给出了一种利用微波脉冲测量血管中的血流量的方法和装置，但其需要将隔开一定距离的同轴收发天线插入血管，通过测量微波脉冲在血管中的时延来测量血管中的血流量。这是一种有损伤、侵入式的测量方法，易造成血液检测过程中的意外感染。

发明内容

本发明的目的是针对现有技术的不足而给出的一种基于微波、射频多普勒效应，对视力可及或浅表的血管内血液的流速进行实时无损检测的方法及装置，本发明利用血流对微波、射频信号的互作用产生的多普勒效应，通过相干检测获取信号准确的多普勒频移，从而测得血管中血流的速度信息。

本发明的目的是这样实现的：

一种基于微波、射频多普勒效应的血液动力学检测方法，其特征在于该方法包括以下具体步骤：

a)、将一发射天线及一接收天线分别紧贴在待测者血管外壁上或血管外的皮肤上，发射电磁波的方向和血管轴向呈一个夹角，该夹角的范围在0～30°；接收天线的位置与发射天线的位置沿血管或紧靠血管两侧，并隔开一定距离；发射天线连接微波或射频收发模块的一信号放大器；接收天线连接微波或射频收发模块的一带通滤波器；

b)、根据血流速度、多普勒信号的频率分辨率、电磁波的传播衰减等因素，设定微波或射频信号的频率；

c)、将微波或射频信号的一部分经功率分配器进入混频器的本振端，作为本振信号，微波或射频信号的另一部分经放大后由发射天线发射，经血流多普勒频移调制的已调微波或射频信号通过接收天线接收，进入混频器的射频信号端，作为测血流的接收信号；本振信号与已调微波或射频信号经混频器相干差频作用后在混频器的低频输出端输出多普勒低频信号；并滤除其中的零频直流信号干扰信号与高频噪声；

d)、经滤除干扰和噪声后的多普勒低频信号经放大后转换为数字信号提供给微处理器；

e)、采用时频分析获得血流速度随时间的变化关系，对心脏压缩和扩张时血流最大速度和最小速度进行跟踪，给出血管中血流最大速度和最小速度值，提供一个心脏周期血管中的血流速度变化信息图。

实施上述方法的装置包括天线、微波或射频收发模块、信号处理电路及电源，所述天线包括发射天线及接收天线，或称为收发天线探头；所述微波或射频收发模块由微波或射频信号发生器、功率分配器、混频器、第一信号放大器、第二信号放大器、第一带通滤波器及第二带通滤波器组成，微波或射频信号发生器连接功率分配器，功率分配器分为两路，一路经第二信号放大器连接发射天线，另一路连接混频器的本振端；混频器的信号端经第一信号放大器、第二带通滤波器连接接收天线，其输出端经第一带通滤波器与信号处理电路连接；所述电源分别连接微波或射频收发模块与信号处理电路。

所述天线为一异型开口同轴天线，由硬同轴电缆制成，其一端沿垂直于轴线方向切成平面，并在其侧面切割出一个槽，以保证辐射方向图满足测试要求。

所述信号处理电路由模数变换器、微处理器及第三信号放大器组成，所述微处理器经模数变换器连接第三信号放大器，第三信号放大器与微波或射频收发模块中的第一带通滤波器连接。

本发明中设计了带匹配网络的异形开口同轴天线，保证收发端天线的匹配和测试所需的辐射方向图；通过在微波和射频频段优化信号频率，以获得能兼顾微波或射频信号对人体的穿透深度和多普勒信号的分辨率的性能；通过设计超窄带带通滤波器，以抑制设备测试中的强零频干扰，并采用了弱信号处理手段以减低对探头位置的依赖和减少外界干扰的影响，提高血液流速检测的实时性，精确度和方便性；通过应用先进的时频联合分析方法，能够对心脏压缩和扩张时血液最大速度和最小速度进行跟踪，并可根据需要，提供一个心脏周期中的相关血管中的血流速度变化信息，为医生提供可靠的诊断信息，可有效降低血液检测过程中意外感染的危险性。

附图说明

图1为本发明检测装置结构框图

图2为本发明检测装置中天线的结构示意图

图3为一个心脏活动周期的模拟时域多普勒信号图

图4为模拟时域多普勒信号的伪WV分布重排结果图

具体实施方式

参阅图1、图2，本发明的检测装置包括天线1、微波或射频收发模块2、信号处理电路3及电源4，所述天线1包括发射天线12及接收天线11；所述微波或射频收发模块2由微波或射频信号发生器21、功率分配器22、混频器23、第一信号放大器24、第二信号放大器25、第一带通滤波器26及第二带通滤波器28组成，微波或射频信号发生器21连接功率分配器23，功率分配器23分为两路，一路经第二信号放大器25连接发射天线12，另一路连接混频器23的本振端；混频器23的射频信号端经第一信号放大器24、第二带通滤波器28连接接收天线11，混频器23的低频输出端经第一带通滤波器26与信号处理电路3连接；所述电源4分别连接微波或射频收发模块2及信号处理电路3。所述天线1由硬同轴电缆制成，其一端沿垂直于轴线方向切成平面14，并在其侧面切割出一个槽15。所述信号处理电路3由模数变换器31、微处理器32及第三信号放大器33组成，微处理器32经模数变换器31连接第三信号放大器33，第三信号放大器33与微波或射频收发模块2中的第一带通滤波器26连接。

下面对本发明进行详细说明：

本发明天线包括收、发共两个，也可采用单天线加收发双工器，考虑实际应用中对天线探头的尺寸要求，天线可设计为小尺寸的开口同轴天线。开口同轴天线由硬同轴电缆加工优化而成。优化是指根据所需要的天线辐射、阻抗特性在同轴线前端优选开槽、切割的形状与位置，来满足系统对收发微波或射频信号的性能要求。由于人体血管是由高介电常数的有耗材料组成，与自由空间有较大的差别，收、发天线的输入端设计了匹配网络，以改善天线与系统的匹配特性。收、发天线安装在探头中，探头的材料要求对微波或射频信号透明，并能满足医学消毒要求。使用时探头根据需要，或紧贴在血管外壁上(如手术中)，或紧贴在血管外的皮肤上(如浅表血管中血流流速探测的应用)，电磁波发射方向和血管轴向呈一个角度θ，θ的范围在0到30°之间。发射天线发射的部分电磁波穿透皮肤、血管壁，进入血液，其散射波被接收天线所接收到。接收天线的位置与发射天线的位置沿血管或紧靠血管两侧排列，并隔开一定距离。由于血液的流动导致的多普勒效应，使得接收天线所收到的电磁波频率与发射信号的频率相比有一个频率偏移，且偏移大小与血流速度成正比。血流速度与多普勒频率的关系为

$f_d=\frac{f_0}{c}\·v\·\cos \mathrm{\θ}---\left(I\right)$

(I)式中的v为血流速度，f_d为多普勒频率，f₀为微波或射频信号的频率，θ为波的传播方向与血流方向的夹角，c为真空中的光速，当环境媒质非真空时，c改为电磁波在该媒质中传播的速度。f₀可在射频频段或微波频段中选择，f₀的频率越低，电磁波对血管壁和血液的穿透能力越强，但从(I)式中可以看到，测到的多普勒频率会很低，频率分辨力将下降。另外收发天线的尺寸将增加，对构筑测试设备不利。f₀提高的结果正好相反。f₀的选择还需考虑的是减小外界干扰信号的影响，以及满足国家无线电管理委员会对频谱管理的要求。本发明中通过优选f₀来保证设备的综合性能。

微波或射频收发模块2由微波或射频信号发生器21，放大器24、25，功率分配器22，带通滤波器26、28，混频器23和电源4组成，微波或射频信号发生器21能够产生稳定的、频率可调的微波或射频信号，放大器分为微波或射频放大器25和低频信号放大器24，分别将信号源产生的信号放大到足够的电平至发射天线输出，以及将接收到的经血流多普勒频移调制的微波或射频信号放大至混频器23的射频信号端，为了保证检测设备检测到的信号的频率稳定性，本发明采用了相干接收技术，即将微波或射频信号发生器21的输出通过功率分配器22把微波或射频信号一分为二，一路送发射天线12发射，另一路送混频器23作为本振信号。经发射天线12发射的部分微波或射频信号，与血管中的血流相互作用之后，产生信号频率f₁＝f₀±f_d的新信号，该信号与发射信号f₀同时被接收天线11接收，通过放大器24放大之后，进入混频器23的射频信号端，与混频器23本振端的本振信号f₀差频作用之后，在混频器23输出端输出含零频直流成分的多普勒低频信号。由于直射信号比散射信号大得多，零频直流信号对多普勒信号会产生很大的干扰，本发明中采用一个超窄带的带通滤波器26，以滤除零频干扰信号与高频噪声。

多普勒低频信号的采集由可变增益低频放大器33与模数(A/D)变换器31完成，可变增益低频放大器33将带通滤波器26输出的信号放大至模数变换器31所需的合适的输入电平，模数变换器31将模拟多普勒信号转换成数字信号，以便进行血流速度的数字信号处理。

上述数字信号的处理可由通用微处理机或高性能嵌入式微处理器完成，由于人体血流速度因心脏的运动具有时变特性，本发明对测得的多普勒信号首先进行模数变换，然后再采用时频分析技术获得血流速度随时间的变化关系，能够对心脏压缩和扩张时血流最大速度和最小速度进行跟踪，提供血管中血流最大速度和最小速度值，并可根据需要，提供一个心脏周期中的相关血管中的血流速度变化信息。

为了帮助更好的理解本发明，将结合附图给出进一步的说明，并给出一个实施例。本发明的保护范围不仅限于所提供的实施例方案。

实施例

参阅图2，将一段特性阻抗为50欧姆的硬同轴电缆1的一端沿垂直轴线方向切成平面14，并在其侧面切割出一个槽15；这里切割的角度、切割后所获得平面的大小、开槽的位置与长度可以利用电磁波理论计算获得。目的是使得辐射电磁波能量能够尽可能多地进入血管，并按照预期的方向传播。然后把两个这样的天线分别安装在如图1所示的微波或射频收发模块2的发射端和接收端。图1中的微波或射频收发模块2能够产生所需要频率稳定的正弦波信号，比如可采用石英晶体稳频的锁相微波源，信号频率为430MHz，输出功率为10dBm。该信号经过功率分配器22分为两路，功率分配器22采用10dB的定向耦合器，一路约0dBm的信号经过一个10dB增益的信号放大器25放大输送约10dBm的微波功率给发射端通过发射天线12发射出去，另外一路约9.5dBm的信号直接输送给混频器23作为本振信号，实现相干检测。混频器23可由一个能够输出零中频的混频器和相应的滤波器来实现，经血流散射有多普勒频移的微波或射频信号被接收天线11接收，信号放大器24把从接收天线11接收而来的微弱电磁波信号放大后送到混频器23，混频器23比较出该信号与本振信号频率的差异，并将这个差值输出到数据采集和模数变换器31做进一步的数字信号处理。本发明以检测直径为数毫米级的动脉和静脉中的血流为主，设计的流速范围为1cm/s到100cm/s，根据式(I)和上述采用的信号源的频率，并考虑到测量中的人体媒质参数，多普勒频率范围大约在0.1Hz～12Hz。由于在测量中除受血流调制的微波或射频散射信号外还有从发射天线直接到达接收天线的直射信号，直射信号要远远大于散射信号，因此混频器23输出除所需的多普勒信号外，还会有一个很强的零频直流电平，其会影响数据的处理。因此混频器23的输出电路中需要一个中心频率低，带宽窄，矩形系数好的滤波器消除零频干扰和高频带外噪声，本实施例采用美信公司的开关电容滤波器芯片MAX7490。这里的中心频率和带宽可根据信号的特点来确定。数据采集和模数变换器31可通过NI公司的带USB接口的数据采集卡，这样接下来的信号处理和数据显示可应用通用微处理机进行，也可采用高性能的嵌入式微处理器芯片，同时实现模数变换和数据处理。

由于人体血流速度因心脏的运动具有时变特性，本实施例对测得的多普勒信号采用时频分析技术，获得血流速度随时间的变化关系，能够对心脏压缩和扩张时血流最大速度和最小速度进行跟踪，提供血管中血流最大速度和最小速度值，并可根据需要，提供一个心脏周期中的相关血管中的血流速度变化信息。

为了分析和处理非平稳信号，人们提出了信号的时频分析法，即利用时间和频率的联合函数来表示信号，时频分析分为线性和二次型两种。二次型的时频分析可用来描述信号的能量密度分布，典型的是Wigner-Ville分布。Wigner-Ville分布可被看作信号能量在时域和频域中的分布，它属于Cohen类时频分布。信号s(t)的Wigner-Ville分布定义为：

$W_z(t,f)={\∫}_{-\∞}^{\∞}z(t+\frac{\mathrm{\τ}}{2})z^{*}(t-\frac{\mathrm{\τ}}{2})e^{-j2\mathrm{\π\τf}}\mathrm{d\τ}---\left(\mathrm{II}\right)$

式中，z(t)是s(t)的解析信号，*表示取共轭。

虽然Wigner-Ville分布具有好的时频聚集性，但是对于多分量信号，其Wigner-Ville分布会出现交叉项，产生“虚假信号”。交叉项是二次型时频分布的固有结果，它来自于多分量信号中不同信号分量之间的交叉作用，而交叉项的抑制又主要通过核函数的设计来实现。常用的加核函数后的Wigner-Ville分布有伪Wigner-Ville分布(PWD)：

${\mathrm{PWD}}_z(t,f)={\∫}_{-\∞}^{\∞}z(t+\frac{\mathrm{\τ}}{2})z^{*}(t-\frac{\mathrm{\τ}}{2})h\left(\mathrm{\τ}\right)e^{-j2\mathrm{\π\τf}}\mathrm{d\τ}---\left(\mathrm{III}\right)$

式中h(τ)是窗函数。

本实施例应用了伪Wigner-Ville分布(PWD)来对血流的多普勒信号进行时频分析处理。为了了解时频信号处理的效果，图3显示了一个心脏活动周期的模拟时域多普勒信号，图4给出了模拟时域多普勒信号的伪WV分布重排结果。从此结果可以获得一个心脏活动周期的血管中血流最大、最小速度所对应的频率，以此可计算出最大、最小流速；也可以从中得到一个心脏活动周期的血流速度随时间的变化关系。

Claims (4)

1.一种基于微波、射频多普勒效应的血液动力学检测方法，其特征在于该方法包括以下具体步骤：

a)、将一发射天线探头及一接收天线探头分别紧贴在待测者血管外壁上或血管外的皮肤上，发射电磁波的传播方向和血管轴向呈一个夹角，该夹角的范围在0～30°；接收天线的位置与发射天线的位置沿血管或紧靠血管两侧，并隔开一定距离；发射天线连接微波或射频收发模块的一信号放大器；接收天线连接微波或射频收发模块的一带通滤波器；

b)、根据血流速度、多普勒信号的频率分辨率、电磁波的传播衰减因素，设定微波或射频信号的频率；

c)、将微波或射频信号的一部分经功率分配器进入混频器的本振端，作为本振信号，微波或射频信号的另一部分经放大后由发射天线发射，经血流多普勒频移调制的、已调微波或射频信号通过接收天线探头接收，进入混频器的信号端，作为测血流的接收信号；本振信号与已调微波或射频信号经混频器相干差频作用后在混频器的低频输出端输出多普勒低频信号；并滤除其中的零频直流信号、干扰信号与高频噪声；

d)、经滤除干扰和噪声后的多普勒低频信号经放大后转换为数字信号提供给微处理器；

e)、采用时频分析获得血流速度随时间的变化关系，对心脏压缩和扩张时血流最大速度和最小速度进行跟踪，给出血管中血流最大速度和最小速度值，提供一个心脏周期血管中的血流速度变化信息图。

2.一种实施权利要求1所述方法的装置，其特征在于该装置包括天线、微波或射频收发模块、信号处理电路及电源，所述天线包括发射天线及接收天线，或称为收发天线探头；所述微波或射频收发模块由微波或射频信号发生器、功率分配器、混频器、第一信号放大器、第二信号放大器、第一带通滤波器及第二带通滤波器组成，微波或射频信号发生器连接功率分配器，功率分配器分为两路，一路经第二信号放大器连接发射天线，另一路连接混频器的本振端；混频器的信号端经第一信号放大器、第二带通滤波器连接接收天线，其输出端经第一带通滤波器与信号处理电路连接；所述电源分别连接微波或射频收发模块与信号处理电路。

3.根据权利要求2所述的装置，其特征在于所述天线为一异型开口同轴天线，其由硬同轴电缆制成，其一端沿垂直于轴线方向切成平面，并在其侧面切割出一个槽，以保证辐射方向图满足测试要求。

4.根据权利要求2所述的装置，其特征在于所述信号处理电路由模数变换器、微处理器及第三信号放大器组成，所述微处理器经模数变换器连接第三信号放大器，第三信号放大器与微波或射频收发模块中的第一带通滤波器连接。

Images