CN115568876A - 血流速度测量值的修正方法及多普勒血流检测装置 - Google Patents

Abstract

本申请实施例提供了一种血流速度测量值的修正方法及多普勒血流检测装置，其中，血流速度测量值的修正方法应用于多普勒血流检测装置；方法包括：获取经由第一晶元组和第二晶元组对待检测血管进行检测得到的第一晶元组血流速度测量值和第二晶元组血流速度测量值；确定探头在检测状态下的滚转角，滚转角为探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向之间的夹角，发射晶元为第一晶元组和/或第二晶元组中的发射晶元；根据滚转角对第一晶元组血流速度测量值和第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值；根据第一修正值和第二修正值确定待输出的血流速度测量值。

Description

血流速度测量值的修正方法及多普勒血流检测装置

技术领域

本发明涉及医疗检测领域，特别是涉及一种血流速度测量值的修正方法及多普勒血流检测装置。

背景技术

随着科技的发展，对于人体出现的各类疾病，可以通过医疗检测手段获取人体的相关参数来帮助医生进行诊断，其中血流速度是一项重要的生理参数，许多疾病都能在血流速度上得到反映。

利用多普勒技术进行血流速度的检测，其基本原理是超声探头发出一定强度的超声波，当超声波遇到运动的血液时，由于超声波的波长大于血液中的红细胞的直径，超声波遇到红细胞时将发生散射。同时由于红细胞在运动，所以散射的回波信号发生了多普勒频移，速度越大，频移也就越大。通过分析回波信号中的频移量，结合多普勒技术的数学公式推算出红细胞的移动速度。

目前，使用多普勒血流检测装置对患者进行血流速度检测在临床上已经得到了广泛应用，但是在实际使用中，测量结果往往存在着一定的误差，尤其是在医护人员操作不标准，或者患者在长时间监测过程中出现移动时，测量结果的误差会更大，这导致血流速度测量结果的可参考性被降低，严重时可能会造成医护人员对患者病情的误判。

发明内容

有鉴于此，本申请实施例为解决背景技术中存在的至少一个问题而提供一种血流速度测量值的修正方法及多普勒血流检测装置。

第一方面，本申请一实施例提供了一种血流速度测量值的修正方法，应用于多普勒血流检测装置，所述多普勒血流检测装置包括探头，所述探头包括用于放置在待检测血管不同位置处的第一晶元组和第二晶元组；所述方法包括：

获取经由所述第一晶元组对所述待检测血管进行检测得到的第一晶元组血流速度测量值；

获取经由所述第二晶元组对所述待检测血管进行检测得到的第二晶元组血流速度测量值；

确定所述探头在检测状态下的滚转角，所述滚转角为所述探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向之间的夹角，所述发射晶元为所述第一晶元组和/或所述第二晶元组中的发射晶元；

根据所述滚转角对所述第一晶元组血流速度测量值和所述第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值；

根据所述第一修正值和所述第二修正值确定待输出的血流速度测量值。

结合本申请的第一方面，在一可选实施方式中，所述多普勒血流检测装置还包括设置在探头内的姿态传感器；所述确定所述探头在检测状态下的滚转角，包括：

利用所述姿态传感器确定所述探头在检测状态下的滚转角。

结合本申请的第一方面，在一可选实施方式中，所述姿态传感器包括三轴加速度传感器；所述确定所述探头在检测状态下的滚转角，包括：

基于所述三轴加速度传感器获得在检测状态下重力加速度在相互垂直的三个坐标轴方向上的加速度分量；

根据所述加速度分量确定所述滚转角。

结合本申请的第一方面，在一可选实施方式中，所述姿态传感器包括陀螺仪；所述确定所述探头在检测状态下的滚转角，包括：

基于所述陀螺仪获得在检测状态下的所述滚转角。

结合本申请的第一方面，在一可选实施方式中，所述根据所述滚转角对所述第一晶元组血流速度测量值和所述第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值，包括采用以下公式得到第一修正值和第二修正值：

V1 = V / sin

；其中，V1表示第一修正值或第二修正值，V表示第一晶元组血流速度测量值或第二晶元组血流速度测量值，

表示滚转角。

结合本申请的第一方面，在一可选实施方式中，所述根据所述第一修正值和所述第二修正值确定待输出的血流速度测量值，包括：

将所述第一修正值和所述第二修正值的平均值确定为所述待输出的血流速度测量值。

结合本申请的第一方面，在一可选实施方式中，所述获取第一晶元组血流速度测量值和/或所述获取第二晶元组血流速度测量值，包括：

获取发射角，所述发射角为所述多普勒血流检测装置中发射晶元发射信号的方向与待检测血管中血流运动方向之间的夹角；

获取所述多普勒血流检测装置中位于同一晶元组中的所述发射晶元和接收晶元之间的距离以及所述发射晶元与所述待检测血管之间的距离；

根据所述发射角、所述发射晶元和接收晶元之间的距离、以及所述发射晶元与所述待检测血管之间的距离确定接收角，所述接收角为所述多普勒血流检测装置中接收晶元接收信号的方向与血流运动方向之间的夹角；

根据所述发射角和所述接收角确定所述第一晶元组血流速度测量值和/或所述第二晶元组血流速度测量值。

结合本申请的第一方面，在一可选实施方式中，所述发射晶元与所述待检测血管之间的距离包括第一距离，所述第一距离为沿所述发射晶元发射信号的方向上所述发射晶元到所述待检测血管的距离；

所述根据所述发射角、所述发射晶元和接收晶元之间的距离、以及所述发射晶元与所述待检测血管之间的距离确定接收角，包括采用以下公式确定接收角：

| β - α | = arctan ( D / L )；其中，β为接收角，α为发射角，| β - α | 为接收角与发射角之间的差值的绝对值，D为发射晶元和接收晶元之间的距离，L为第一距离。

结合本申请的第一方面，在一可选实施方式中，获取所述发射晶元与所述待检测血管之间的距离，包括：获取所述发射晶元发射信号至所述接收晶元接收所述信号所需的时长；根据所述所需的时长以及所述信号的传播速度确定所述信号在所述所需的时长内传播的距离；根据所述传播的距离确定所述第一距离。

第二方面，本申请一实施例提供了一种多普勒血流检测装置，用于实现上述第一方面中任意一项所述的血流速度测量值的修正方法。

本申请实施例所提供的血流速度测量值的修正方法及多普勒血流检测装置，其中，血流速度测量值的修正方法应用于多普勒血流检测装置，多普勒血流检测装置包括探头，探头包括用于放置在待检测血管不同位置处的第一晶元组和第二晶元组；方法包括：获取经由第一晶元组对待检测血管进行检测得到的第一晶元组血流速度测量值；获取经由第二晶元组对待检测血管进行检测得到的第二晶元组血流速度测量值；确定探头在检测状态下的滚转角，滚转角为探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向之间的夹角，发射晶元为第一晶元组和/或第二晶元组中的发射晶元；根据滚转角对第一晶元组血流速度测量值和第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值；根据第一修正值和第二修正值确定待输出的血流速度测量值；如此，采用放置在待检测血管不同位置处的至少两个晶元组分别进行测量，即便出现实际测量中的探头与待检测血管之间的夹角与装置中设定的多普勒夹角不一致的情况，也能够根据各晶元组的血流速度测量值综合确定测量结果，减小上述情况造成的误差，并且通过对第一晶元组和第二晶元组在测量过程中由于发射晶元发射的声场方向出现偏移导致的误差进行修正，进一步提高了各晶元组测量结果的准确性，最终能够输出可靠性更高的血流速度测量值，为医护人员对患者病情的判断提供更为可靠的依据。

本申请附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本申请的实践了解到。

附图说明

此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：

图1为本申请一实施例提供的多普勒血流检测装置的结构框图；

图2为本申请一实施例提供的血流速度测量值的修正方法的流程示意图；

图3为滚转角不等于90°时三轴加速度传感器与待检测血管的位置关系示意图；

图4为偏航角、俯仰角和滚转角的示意图；

图5为在有滚转角误差时和在无滚转角误差时对待检测血管进行检测的位置示意图；

图6为一具体应用中采用多普勒血流检测装置对待检测血管进行检测时的位置示意图；

图7为采用多普勒血流检测装置对浅表的血管和深层的血管进行检测时的接收角对比示意图；

图8为采用一种多普勒血流检测装置对待检测血管进行检测的示意图；

图9为采用另一种多普勒血流检测装置对待检测血管进行检测的示意图；

图10为本申请一具体示例提供的利用超声波确定第一距离的原理示意图；

图11为本申请一具体示例提供的利用超声波确定第一距离的电路框图；

图12为本申请一具体示例提供的利用超声波确定第一距离的流程示意图；

图13为本申请一具体示例提供的采用包括第一接收晶元和第二接收晶元的第一晶元组和/或第二晶元组对待检测血管进行检测时的位置示意图。

具体实施方式

为使本发明的技术方案和有益效果能够更加明显易懂，下面通过列举具体实施例的方式进行详细说明。其中，附图不一定是按比例绘制的，局部特征可以被放大或缩小，以更加清楚的显示局部特征的细节；除非另有定义，本文所使用的技术和科学术语与本申请所属的技术领域中的技术和科学术语的含义相同。

除另作定义外，本申请所涉及的技术术语或者科学术语应具有本申请所属技术领域具备一般技能的人所理解的一般含义。在本申请中的“一”、“一个”、“一种”、“该”、“这些”等类似的词并不表示数量上的限制，它们可以是单数或者复数。在本申请中所涉及的术语“包括”、“包含”、“具有”及其任何变体，其目的是涵盖不排他的包含；例如，包含一系列步骤或模块（单元）的过程、方法和系统、产品或设备并未限定于列出的步骤或模块（单元），而可包括未列出的步骤或模块（单元），或者可包括这些过程、方法、产品或设备固有的其他步骤或模块（单元）。在本申请中所涉及的“多个”是指两个或两个以上。“和/或”描述关联对象的关联关系，表示可以存在三种关系，例如，“A和/或B”可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。在本申请中所涉及的术语“第一”、“第二”、“第三”等，只是对相似对象进行区分，并不代表针对对象的特定排序。

本申请中，患者指的是正在接受或者可能接受多普勒血流检测装置的检测的任何人，患者也可被称为病人。医护人员，例如包括使用多普勒血流检测装置对患者进行监测的人，因而，也作为使用多普勒血流检测装置的用户之一；但，用户不仅可以包括医护人员，还可以包括其他直接或间接使用多普勒血流检测装置的人员。

本申请实施例提供了一种血流速度测量值的修正方法，应用于多普勒血流检测装置。首先，请参考图1，图1示出了本申请一实施例提供的多普勒血流检测装置的结构。如图所示，多普勒血流检测装置，包括：探头，信号处理模块，以及装置主体模块。具体的，探头可以包括第一晶元组、第二晶元组和姿态传感器。其中，第一晶元组和第二晶元组分别与信号处理模块连接，从而将探测到的信号传输至信号处理模块；信号处理模块在接收到由各晶元组传输的信号后，对上述信号进行处理。具体地，信号处理模块包括滤除电路1、滤除电路2、信号解调1、信号解调2、ADC1、以及ADC2；第一晶元组与信号处理模块中的滤除电路1连接，第二晶元组与信号处理模块中的滤除电路2连接，滤除电路1和滤除电路2分别对第一晶元组和第二晶元组探测到的信号进行滤波处理；滤除电路1与信号解调1连接，滤除电路2与信号解调2连接，信号解调1和信号解调2分别对经滤除电路1和滤除电路2滤波处理后的信号进行解调处理；信号解调1与ADC1连接，信号解调2与ADC2连接，ADC1和ADC2分别对经信号解调1和信号解调2解调处理后的信号进行模数转换处理。经信号处理模块处理后的信号传输至处理器。姿态传感器与处理器连接，将探测到的传感器信号传输至处理器。为了存储接收到的检测结果，以及存储运算过程中的中间量和算法程序，该装置中还包括与处理器连接的存储模块，如内存、存储器等。当然，该装置中还可以包括网络口。

通讯接口为和其它智能设备的通讯提供接口，以进行数据的传输，实现反馈。通讯接口可以是USB/UART/网口/蓝牙/WIFI/CAN等常用的通讯硬件接口，本申请对其不做具体限制。

此外，该装置还可以包括音频功放以及与音频功放连接的喇叭。喇叭可以用于对处理器确定的结果进行音频输出以及对需要提醒用户的信息进行音频输出。

本具体示例中，该装置还可以包括输入设备；其中，输入设备包含但不限于以下至少之一：键盘、按键、声控输入装置、触摸屏。可以理解的是，其他能够实现信号输入的设备均包含在本申请的含义中。此外，该装置还可以包括输出设备；其中，输出设备包含但不限于LCD显示屏。当然，其他能够实现信号输出的设备也均包含在本申请的含义中。如此，该装置能够实现与用户的信息交互。

本申请实施例提供的血流速度测量值的修正方法，应用于多普勒血流检测装置。多普勒血流检测装置包括探头，探头包括用于放置在待检测血管不同位置处的第一晶元组和第二晶元组。其中，多普勒血流检测装置例如为图1所示的多普勒血流检测装置。

请参考图2，该血流速度测量值的修正方法包括：

步骤201、获取经由第一晶元组对待检测血管进行检测得到的第一晶元组血流速度测量值；

步骤202、获取经由第二晶元组对待检测血管进行检测得到的第二晶元组血流速度测量值；

步骤203、确定探头在检测状态下的滚转角；

其中，滚转角为探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向之间的夹角，发射晶元为第一晶元组和/或第二晶元组中的发射晶元；

步骤204、根据滚转角对第一晶元组血流速度测量值和第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值；

步骤205、根据第一修正值和第二修正值确定待输出的血流速度测量值。

可以理解的，本申请实施例提供的血流速度测量值的修正方法，采用放置在待检测血管不同位置处的至少两个晶元组分别进行测量，即便出现实际测量中的探头与待检测血管之间的夹角与装置中设定的多普勒夹角不一致的情况，也能够根据各晶元组的血流速度测量值综合确定测量结果，减小上述情况造成的误差，并且通过对第一晶元组和第二晶元组在测量过程中由于发射晶元发射的声场方向出现偏移导致的误差进行修正，进一步提高了各晶元组测量结果的准确性，最终能够输出可靠性更高的血流速度测量值，为医护人员对患者病情的判断提供更为可靠的依据。

其中，第一晶元组和第二晶元组用于放置在待检测血管不同位置处，例如，第一晶元组和第二晶元组中的一者放置在待检测血管的第一位置处，而另一者放置在待检测血管的第二位置处，第一位置和第二位置分别为沿血流运动方向上待检测血管的前后两个位置，即血流运动方向为由第一位置至第二位置的方向。应当理解，待检测血管被皮肤包裹，第一晶元组和第二晶元组显然并不是直接与待检测血管接触，本文所述的第一晶元组和第二晶元组放置在待检测血管不同位置处，指的是第一晶元组和第二晶元组隔着皮肤放置在与待检测血管对应的不同位置处。

请继续参考图1，多普勒血流检测装置还包括设置在探头内的姿态传感器；确定探头在检测状态下的滚转角，包括：利用姿态传感器确定探头在检测状态下的滚转角。

姿态传感器是一种高性能的三维运动姿态测量仪器。姿态传感器可以包括以下至少之一：三轴加速度传感器、陀螺仪、三轴电子罗盘等运动传感器。

作为一种可选的实施方式，姿态传感器包括三轴加速度传感器；确定探头在检测状态下的滚转角，包括：基于三轴加速度传感器获得在检测状态下重力加速度在相互垂直的三个坐标轴方向上的加速度分量；根据加速度分量确定滚转角。

在具体应用中，开始进行血流速度的测量时，用户将探头按照预先规定的使用方式放置在测量位置，其中，预先规定的使用方式例如为保证探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向相垂直，不产生离散角，从而减小测量误差。在测量的过程中，可以认为病人是静止的。此时可以记录三轴加速度传感器的初始输出值，为后续计算得到的滚转角做参考。

可以理解的，以待检测血管位于水平面内为例，如果三轴加速度传感器的芯片是水平放置的，并且以三轴加速度传感器的方向建立包括XYZ三个坐标轴的直角坐标系，那么重力加速度g在X轴和Y轴方向的分量均为0，在Z轴方向的分量为g。

而在实际测量时，特别是长时间的监测中，很可能会出现由于病人移动导致探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向不垂直的情况，从而造成测量误差。

图3示出了滚转角不等于90°时三轴加速度传感器与待检测血管的位置关系，由于三轴加速度传感器位于探头内，因此图3也即探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向不垂直的情况三轴加速度传感器与待检测血管的位置关系。如图所示，重力加速度g在X轴方向的分量为Ax，Ax与待检测血管101之间的夹角为a1，与重力加速度g的夹角为a；重力加速度g在Y轴方向的分量为Ay，Ay与待检测血管101之间的夹角为b1，与重力加速度g的夹角为b；重力加速度g在Z轴方向的分量为Az，Az与待检测血管101之间的夹角为γ1，与重力加速度g的夹角为γ。其中，各个角度之间的关系为：a =90°- a1，b = 90°- b1，γ = 90°- γ1。

重力加速度g在各轴方向上的分量为：Ax = gcosa，Ay = gcosb，Az = gcosγ；或者：Ax = gsina1，Ay = gsinb1，Az = gsinγ1。

进一步的：

g*g = Ax*Ax + gcosa1*gcosa1；

gcosa1 = squr(g*g - Ax*Ax)；

同理可知，

gcosb1 = squr (g*g - Ay*Ay)；

gcosγ1 = squr(g*g - Az*Az)。

在立体几何中，g相当于立方体的对角线，Ax、Ay、Az相当于三条边，所以根据勾股定理：Ax*Ax + Ay*Ay + Az*Az = g*g。

以X轴为例，sina1 = Ax/g，cosa1 = squr(g*g - Ax*Ax) / g ；

tana1 = Ax / squr(g*g - Ax*Ax) = Ax / squr(Ay*Ay + Az*Az) 。

同理可知，

tanb1 = Ay / squr(Ax*Ax+Az*Az)；

tanγ1 = Az / squr(Ax*Ax +Ay*Ay)。

由此得出了三轴加速度传感器的测量值与a1、b1和γ1之间的关系。这里，a1、b1和γ1的值与对应的弧度值相等。

接下来，计算各轴与水平线（参考图中待检测血管101）之间的弧度值，其中，采用a11 、b11 、γ11分别代表X轴方向、Y轴方向、Z轴方向与水平线之间的弧度值；则：

a11= arctan(Ax / squr(Ay*Ay + Az*Az))；

b11= arctan(Ay / squr(Ax*Ax+Az*Az))；

γ11= arctan(Az / squr(Ax*Ax +Ay*Ay))。

接下来使用以下公式进行计算：弧度 = 角度值*πR/180；即，角度值 = 弧度*180/πR，其中R取1。最后得到各轴的角度值分别为：

；

；

；

其中，

为弧度值b11对应的角度值，

为弧度值γ11对应的角度值，

为弧度值a11对应的角度值；

即为滚转角。

可以理解的，相互垂直的三个坐标轴方向根据三轴加速度传感器在测量状态下的方向建立，三轴加速度传感器与发射晶元共同位于探头内，因此，三轴加速度传感器所处的方向可以代表发射晶元所处的方向；故，根据加速度分量能够确定出发射晶元发射的声场在水平面上的投影方向与血流运动方向之间的夹角。

作为另一种可选的实施方式，姿态传感器包括陀螺仪；确定探头在检测状态下的滚转角，包括：基于陀螺仪获得在检测状态下的滚转角。

图4为偏航角、俯仰角和滚转角的示意图；其中，（a）示出了偏航角

，（b）示出了俯仰角

，（c）示出了滚转角

。将运动物体的空间轨迹分解为沿Z、Y、X三轴的独立运动，运动矢量在每个轴向的分量由该矢量与各轴之间夹角的余弦表示。其中，运动物体例如为姿态传感器。方向余弦共九个参数，可以用一个3阶方向余弦矩阵

来表示。

分量形式如下：

。

运动物体的3次的欧拉旋转可以用三个方向余弦矩阵表示：

；

；

。

余弦矩阵的乘积

用来表示从地理坐标系到运动物体坐标系的变换：

；

；

。

在得到运动物体坐标系相对于地理坐标系的姿态更新矩阵

后，就可以将姿态传感器获取的角速率数据转换为姿态角数据。

，

表示加速度计在各轴方向上感测的重力分量，

表示地理坐标系下在加速度计各轴测量重力场的分量。

利用反三角函数可以得到倾斜角（滚转角

和俯仰角

）：

；

。

在实际使用中，根据人体的特性，可认为Z轴方向和待检测血管方向的相对角度是不变化的，即γ和γ1保持不变，即俯仰角不变，只需要利用滚转角

进行角度修正。

可以理解的，当人体变化位置时，姿态传感器和待检测血管的位置是同步变化的。在多普勒效应中，起作用的是和待检测血液中血流运动方向所在的平面平行的夹角。

根据滚转角对第一晶元组血流速度测量值和第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值，包括采用以下公式得到第一修正值和第二修正值：V1 =V / sin

；其中，V1表示第一修正值或第二修正值，V表示第一晶元组血流速度测量值或第二晶元组血流速度测量值，

表示滚转角。

请参考图5，在按照预先规定的使用方式将探头放置在测量位置时，滚转角

等于90°，此时测量结果无滚转角误差。而在有滚转角误差时，即探头发生转动后，按照速度分量的原理，此时获得的血流速度测量值是无滚转角误差应时应当获得的血流速度测量值在转动后方向上的速度分量，因此通过公式V1 = V / sin

可以获得修正后的血流速度，即或者理论上等于无滚转角误差应时应当获得的血流速度测量值的修正值。

可以理解的，在滚转角

等于90°时，根据上述公式V1 = V / sin

，计算可得V1 =V，进一步说明此时测量结果无滚转角误差。

通过对第一晶元组血流速度测量值和第二晶元组血流速度测量值进行修正，可以得到第一修正值和第二修正值。如前所述，第一晶元组血流速度测量值和第二晶元组血流速度测量值分别为经由第一晶元组和第二晶元组对待检测血管进行检测得到的测量值，而第一晶元组和第二晶元组用于放置在待检测血管不同位置处。请参考图6，在一具体示例中，为了便于理解，可认为第一晶元组和第二晶元组分别放置在左右两个位置处，即可将第一晶元组和第二晶元组分别称为左晶元组和右晶元组。图中虚线框示出了探头与待检测血管之间的夹角与多普勒血流检测装置中设定的多普勒夹角一致的情况，即实际多普勒夹角等于设定的多普勒夹角，二者均等于发射角α。此时左晶元组和右晶元组在图中分别以左0和右0示出。然而在实际测量中，往往难以保证实际多普勒夹角与设定的多普勒夹角完全一致，即会产生测量误差。图中实线框示出了实际多普勒夹角与多普勒血流检测装置中设定的多普勒夹角不一致的情况，此时实际多普勒夹角为α’，左晶元组和右晶元组在图中分别以左1和右1示出。可以理解的，由于获得的第一晶元组血流速度测量值或第二晶元组血流速度测量值V为实际血流速度V0在相应的晶元组发射或接收信号的方向上的垂直分量，因此V=V0/cosα。那么，如果α和α’之间的偏差角度r为1°，则测量误差△V = V0/cosα - V0/cos(α’) = V0/cosα - V0/cos(α±1)。由此可见，设定的多普勒夹角越大，产生的误差就越大。

为了减小误差，本申请实施例所应用的多普勒血流检测装置的探头中包括上述第一晶元组和第二晶元组，通过第一晶元组和第二晶元组实现两路多普勒通道采集，由于第一晶元组和第二晶元组在待检测血管上的放置位置不同，从而两路多普勒通道采集具有两个多普勒夹角，第一晶元组血流速度测量值和第二晶元组血流速度测量值分别基于两个不同的多普勒夹角得到，并且分别进行了上述修正，如此，通过对修正后的第一修正值和第二修正值进行结合计算，可以获得更为准确的血流速度测量值。

其中，第一晶元组中包括至少一个发射晶元和至少一个接收晶元，第二晶元组中同样包括至少一个发射晶元和至少一个接收晶元。第一晶元组和第二晶元组同步发射和采样。第一晶元组血流速度测量值和第二晶元组血流速度测量值对应的测量时间相同。左晶元组中的发射晶元例如从左至右发射信号，右晶元组中的发射晶元例如从右至左发射信号。

在具体应用中，第一晶元组和第二晶元组中发射晶元发射的信号的频率不同，从而多普勒血流检测装置为双频双通道异常探头检测装置。可选的，第一晶元组和第二晶元组发射的信号的频率在1MHz-16MHz的范围内，以满足双频双通道检测的需求。进一步地，第一晶元组和第二晶元组发射的信号的频率在2MHz-8MHz的范围内，检测结果的准确性更好。第一晶元组中发射晶元发射的信号的频率和第二晶元组中发射晶元发射的信号的频率之一为F1，另一为F2，为了提高检测结果的准确性，F1和F2满足如下条件：F1 < F2；F2取F1的余数大于0.5MHz，F2%F1>0.5MHz；且F2 < 3F1。取余数就是取出不能整除的多出来的部分；具体地，如果F2 = nF1+x，其中，n为正整数，x < F1，则F2取F1的余数等于x。

作为一种较佳的示例，F1 = 2.1MHz - 2.7MHz，F2 = 3.8MHz - 4.3MHz。

在一具体示例中，第一晶元组和第二晶元组的信号滤波中做限波滤波器，其中之一的带宽的范围为F1 - 0.5MHz ~ F1 + 0.5MHz；其中另一的带宽范围为F2 - 0.5MHz ~F2 + 0.5MHz。如此，不同的通道设置不同的滤波带宽。在双通道同时发射时，两侧的接收晶元都能够接收到两个发射晶元的信号，也就是互相会有干扰，在本具体示例中，选择了对应通道对应频段的滤波，可以把另外一个频段的信号抑制，达到更好的信噪比。

通过滤波和解调技术去除两个通道中的超声波干扰。

作为一种可选的实施方式，根据第一修正值和第二修正值确定待输出的血流速度测量值，包括：将第一修正值和第二修正值的平均值确定为待输出的血流速度测量值。

下面结合图6，以α = 60°，r = 5°为例，分析通过将第一修正值和第二修正值的平均值确定为待输出的血流速度测量值对减小误差的作用情况：

令△f*c / 2f = A；

实际应当检测得到的血流速度V11 = △f * c / (2f * cosα) = A / cosα = A/ 0.5 = 2A；

如果不采用本实施方式进行修正而检测得到的血流速度V12 = △f * c / (2f *cosα’) = △f * c / (2f * cos(α+r)) = A / cos65° = 2.3663A；

那么，误差为：（V12 – V11）/ V11 = （2.3663A - 2A）/2A = 18.3%。

如果采用本实施方式进行修正，那么图中左1对应的晶元组检测得到的血流速度即为上述V12；

图中右1对应的晶元组检测得到的血流速度V22 = △f * c / (2f * cos(α-r))= A / cos55° = 1.7434A；

为了便于计算，这里先不考虑根据滚转角带来的误差，即认为第一修正值和第二修正值分别与V12和V22相等，则待输出的血流速度测量值等于：（V12 + V22）/ 2 =（2.3663A + 1.7434A）/ 2 = 2.0548A；

采用本实施方式得到的待输出的血流速度测量值的误差为：（2.0548A – 2A）/ 2A= 2.74%。

由此可见，采用本实施方式可以使误差从18.3%减小到2.74%。

如此，采用双通道采集的修正方法极大地减小了由于使用时探头与待检测血管之间的夹角与多普勒血流检测装置中设定的多普勒夹角不一致的问题引起的误差。

作为一种可选的实施方式，在步骤201和步骤202中，获取第一晶元组血流速度测量值和/或获取第二晶元组血流速度测量值，包括：获取发射角，发射角为多普勒血流检测装置中发射晶元发射信号的方向与待检测血管中血流运动方向之间的夹角；获取多普勒血流检测装置中位于同一晶元组中的发射晶元和接收晶元之间的距离以及发射晶元与待检测血管之间的距离；根据发射角、发射晶元和接收晶元之间的距离、以及发射晶元与待检测血管之间的距离确定接收角，接收角为多普勒血流检测装置中接收晶元接收信号的方向与血流运动方向之间的夹角；根据发射角和接收角确定第一晶元组血流速度测量值和/或第二晶元组血流速度测量值。

如此，通过先根据发射角、发射晶元和接收晶元之间的距离、以及发射晶元与待检测血管之间的距离确定接收角，再利用发射角和接收角确定血流速度测量值，可以获得更为准确的血流速度测量值。尤其针对待检测血管为浅表的血管，比如人体颈动脉等情形，采用本实施方式可以获得更贴近真实情况的结果。

具体地，根据发射角和接收角确定第一晶元组血流速度测量值和/或第二晶元组血流速度测量值，包括采用以下公式确定第一晶元组血流速度测量值和/或第二晶元组血流速度测量值：V = Δf * c / ( f * ( cosα + cosβ ))；其中，V表示第一晶元组血流速度测量值或第二晶元组血流速度测量值，Δf表示采用多普勒血流检测装置测量血流速度获得的频率差值，c表示发射晶元发射的信号的传播速度，f为发射晶元发射的信号的频率，α为发射角，β为接收角。

其中，c具体可以为声音在人体组织中的传播速度，一般视组织不同取1570m/s或者1540m/s。

下面，请参考图7，图中1011表示浅表的血管，1012表示深层的血管，T为发射晶元，R为接收晶元，T和R位于同一晶元组中。其中，深层的血管1012距离皮肤表面的深度大于浅表的血管1011距离皮肤表面的深度，相应的，深层的血管1012与晶元组的距离大于浅表的血管1011与晶元组的距离。如图所示，在测量深度较深的深层的血管1012时，发射角α与对应的第二接收角β2之间的差距较小，在计算血流速度测量值时可认为二者近似相等。然而，在测量深度较浅的浅表的血管1011时，超声都还是在近声场，发射角α与对应的第一接收角β1之间的差距较大，一般由于深度问题可能会产生2°~10°的误差。如果不考虑实际的接收角，显然会导致测量结果与实际情况存在较大偏差。需要说明的是，图中仅示意性的说明血管与晶元组的距离不同时，β1和β2会相差较大，图中的角度、尺寸以及比例关系并不严格代表实际情况。

以图8所示的情况为例，其中，发射晶元T相比于接收晶元R更靠近待检测血管101，假设：α=60°，D=5mm，正常人颈动脉离皮肤深度约20mm，由于设计传感器时，晶元组无法完全贴在皮肤表面，加上晶元组安装时的角度，H≈25mm；那么，沿发射晶元T发射信号的方向上发射晶元T到待检测血管101的距离L = 25/sin60°= 25/0.866 = 28.867 mm。(β - α) =arctan(D/L) = arctan(5/28.867) = 9.823°；β = α + (β - α) = 60°+9.823°=69.823°；则cosα+cosβ=cos60°+cos69.823°= 0.5+0.345 = 0.845。

如果按照发射角与接收角近似相等的方法计算血流速度测量值，则V’=Δf * c /( f * 2cosα )。

而实际上，V = Δf * c / ( f * ( cosα + cosβ ))。

误差百分比为：|V’-V|/V = |V’ /V -1| = | ( cosα + cosβ )/2cosα -1| =15.5%。

图9为采用另一种多普勒血流检测装置对待检测血管进行检测的示意图，在该多普勒血流检测装置中，发射晶元T相比于接收晶元R更远离待检测血管101。可以理解的，此种情况下采用公式V = Δf * c / ( f * ( cosα + cosβ ))确定血流速度测量值，同样可以减小误差。具体的误差率，这里不展开计算。

可以理解的，虽然本实施方式中描述了获取第一晶元组血流速度测量值和/或获取第二晶元组血流速度测量值的具体步骤，但是上述具体步骤显然也适用于多普勒血流检测装置中只有一个晶元组的情况。

其中，获取发射角具体可以包括：获取多普勒血流检测装置中预先存储的发射角α。可以理解的，发射角α的具体数值可以是多普勒血流检测装置在设计时就确定下来的固定角度值。

获取多普勒血流检测装置中位于同一晶元组中的发射晶元T和接收晶元R之间的距离D，具体可以包括：获取多普勒血流检测装置中预先存储的位于同一晶元组中的发射晶元T和接收晶元R之间的距离D。该D值也可以是多普勒血流检测装置在设计时就确定下来的固定值。

发射晶元T与待检测血管101之间的距离包括第一距离L，第一距离L为沿发射晶元T发射信号的方向上发射晶元T到待检测血管101的距离；根据发射角α、发射晶元T和接收晶元R之间的距离D、以及发射晶元T与待检测血管101之间的距离L确定接收角β，包括采用以下公式确定接收角β：| β - α | = arctan ( D / L )；其中，β为接收角，α为发射角，| β -α | 为接收角与发射角之间的差值的绝对值，D为发射晶元和接收晶元之间的距离，L为第一距离。

请参考图10，示例性的，获取发射晶元T与待检测血管101之间的距离，包括：获取发射晶元T发射信号至接收晶元R接收信号所需的时长t；根据所需的时长t以及信号的传播速度c确定信号在所需的时长内传播的距离；根据传播的距离确定第一距离L。

实际应用中，获取发射晶元T与待检测血管101之间的距离可以通过A超、B超等技术测量实现。请继续参考图10，发射器T例如为超声发射器，接收器R例如为超声接收器，障碍物例如为待检测血管；发射器T向障碍物发射一小段固定频率的超声波，并在开始发射时进行计时，超声波在物体中有传播的速度c和时间，被障碍物反射后进入到接收器R，从发射到接收的时间可测得为t，超声波在物体中的传播速度为c，则可以计算得到测量距离，根据测量距离可以估算出实际距离。可以通过以下公式确定第一距离：L = c*t/2。

具体的，可以使用原采集电路和传感器的组成，采用PW脉冲波多普勒的方式进行多普勒的采集，并在每个发生周期记录发射到接收的时间差，从而计算测量距离和角度。此方法使用时域数据进行分析计算。

由于超声波在人体组织中的传播速度c = 1570m/s = 1570000mm/s，颈部的血管深度在5mm-30mm之间（即H在5mm-30mm之间），L的范围在7mm-42mm之间。T的范围在7/1570000~42/1570000s，7/1570~42/1570 ms，4.5μs~27μs之间。使用的超声频率可以在2MHz~8MHz，具体地，超声频率为4MHz。以超声频率为4MHz为例，每个发射周期的时长为1/4MHz =0.25μs；那么，发射5个超声周期，总计发射时长为1.25μs，这能够满足采样间隔时间。

图11为本申请一具体示例提供的利用超声波确定第一距离的电路框图；图12为本申请一具体示例提供的利用超声波确定第一距离的流程示意图。

首先，请参考图11。图中电路为总体电路中信号调理部分中的一部分，专门用来进行距离测试。其中，发生控制电路、信号解调和信号放大部分共用。在进行距离测试时，发生控制电路控制发射晶元发出信号，接收晶元接收到信号后，由信号解调和信号放大分别对信号进行解调和放大处理；其中，接收到信号会被输入到ADC（图中未示出）中进行采集。在进行距离测试时，信号进入模拟比较器进行比较，模拟比较器输出为数字信号，直接由处理器的I/O管脚采集。

接下来，请参考图12。利用超声波确定第一距离的流程包括：在开始距离测试后，发射晶元发射超声波信号，并且开始计时，关闭发射；超声波信号通过组织传输到血液相撞，发生多普勒效应后反射至接收晶元，接收晶元对超声波信号进行i/q解调，解调信号输出为微弱的模拟信号，信号放大电路把信号放大至设定的范围。当血流速度越大时，信号会越大。设定一个比较电压V0，当信号大于V0时，模拟比较器输出为1；当信号小于V0时，模拟比较器输出为0。

确定每一个测距周期的方法为：从发射晶元发射超声波信号并开始计时，如果模拟比较器输出从0转为1，则认为是该周期测得的超声传输时间。具体地，判断比较器触发的结果为是还是否，如果结果为是，则记录时间N；如果结果为否，则等待时间到，直至总采集时间到。

由于在一个脉动周期中，动态血流的速度是从小到大再转小，从设计电路中，血流速度越大，信号越大，会越早的触发模拟比较器，这种情况下测得的时间越准确。所以，可选的，设定一个时间段，如，2s-5s的测量时间，优选为3s，保证在一次测量时间中能够测得最大的血流状态下的信号。3s的情况下，被测患者的脉搏是20bpm，能够满足到临床的要求。

如图12所示，从发射超声波到总采集时间到，记录的时间从小到大排序，去除异常点，选择三个最小的时间取平均，距离测试结束。

根据公式L = c*t/2，确定第一距离。再通过如下公式计算得到接收角：| β - α |= arctan ( D / L )。

作为另一种可选的实施方式，获取发射晶元T与待检测血管101之间的距离，包括：获取由用户基于输入设备输入的第二距离H，第二距离H为沿垂直于皮肤表面的方向上发射晶元T到待检测血管101的距离；根据第二距离H确定第一距离L，其中采用以下公式确定第一距离L：L = H / sinα；其中，H为第二距离。

其中，第二距离H例如由医护人员输入至多普勒血流检测装置中。

接下来，请参考图13。可选的，位于第一晶元组和/或第二晶元组中的接收晶元包括第一接收晶元R1和第二接收晶元R2，第一接收晶元R1和第二接收晶元R2接收同一发射晶元T发射的信号，发射晶元T和第一接收晶元R1之间的距离D1与发射晶元T和第二接收晶元R2之间的距离D2相同；根据发射角、发射晶元和接收晶元之间的距离、以及发射晶元与待检测血管之间的距离确定接收角，包括：确定第一接收角β1和第二接收角β2；其中，第一接收角β1为第一接收晶元R1接收信号的方向与血流运动方向之间的夹角，第二接收角β2为第二接收晶元R2接收信号的方向与血流运动方向之间的夹角；根据发射角和接收角确定血流速度测量值，包括：根据发射角α和第一接收角β1确定第一血流速度测量值V10；根据发射角α和第二接收角β2确定第二血流速度测量值V20；根据第一血流速度测量值V10和第二血流速度测量值V20确定血流速度测量值V。

其中，发射晶元T和第一接收晶元R1之间的距离D1与发射晶元T和第二接收晶元R2之间的距离D2相同，所以α-β2 ≈ β1-α，设此角度差值为ε。发射晶元T发射信号时，第一接收晶元R1和第二接收晶元R2能够同时接收到信号。那么，第一接收晶元R1和第二接收晶元R2对应的速度分量相等，Δf1 = Δf2 = Δf*cosε。

假设α = 60°，通过ε = 9.823°，因此β1 = α+ε = 69.823°，β2 = α-ε = 50.177°。

令△f*c / 2f = A；

则如果不考虑发射角和接收角之间的差异，实际应当检测得到的血流速度V11 =△f * c / (2f * cosα) = A / cosα = A / 0.5 = 2A；

而在本方法中，根据第一血流速度测量值V10和第二血流速度测量值V20确定血流速度测量值V。具体地，例如将第一血流速度测量值V10和第二血流速度测量值V20的平均值确定为血流速度测量值V。那么，

V10 = △f * c / (f * (cosα + cosβ1)) = 2A / (cos60° + cos69.823°) =2.332A；

V20 = △f * c / (f * (cosα + cosβ2)) = 2A / (cos60°+ cos50.177°) =1.748A；

则V = (V1+V2)/2 =2.04A。

如此，使用本方法可以减小误差。且β越小时，误差越小。

在此基础上，本申请实施例还提供了一种多普勒血流检测装置，用于实现上述实施例中任意一项的血流速度测量值的修正方法。

应当理解，关于多普勒血流检测装置的各细节技术特征可以参考上述血流速度测量值的修正方法实施例中的相关描述，此处不再赘述。

应该理解，虽然本申请各流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。

应当理解，以上实施例均为示例性的，不用于包含权利要求所包含的所有可能的实施方式。在不脱离本公开的范围的情况下，还可以在以上实施例的基础上做出各种变形和改变。同样的，也可以对以上实施例的各个技术特征进行任意组合，以形成可能没有被明确描述的本发明的另外的实施例。因此，上述实施例仅表达了本发明的几种实施方式，不对本发明专利的保护范围进行限制。

Claims (10)

1.一种血流速度测量值的修正方法，应用于多普勒血流检测装置，其特征在于，所述多普勒血流检测装置包括探头，所述探头包括用于放置在待检测血管不同位置处的第一晶元组和第二晶元组；所述方法包括：

获取经由所述第一晶元组对所述待检测血管进行检测得到的第一晶元组血流速度测量值；

获取经由所述第二晶元组对所述待检测血管进行检测得到的第二晶元组血流速度测量值；

确定所述探头在检测状态下的滚转角，所述滚转角为所述探头中发射晶元发射的声场在水平面上的投影方向与待检测血管中血流运动方向之间的夹角，所述发射晶元为所述第一晶元组和/或所述第二晶元组中的发射晶元；

根据所述滚转角对所述第一晶元组血流速度测量值和所述第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值；

根据所述第一修正值和所述第二修正值确定待输出的血流速度测量值。

2.根据权利要求1所述的血流速度测量值的修正方法，其特征在于，所述多普勒血流检测装置还包括设置在探头内的姿态传感器；所述确定所述探头在检测状态下的滚转角，包括：

利用所述姿态传感器确定所述探头在检测状态下的滚转角。

3.根据权利要求2所述的血流速度测量值的修正方法，其特征在于，所述姿态传感器包括三轴加速度传感器；所述确定所述探头在检测状态下的滚转角，包括：

基于所述三轴加速度传感器获得在检测状态下重力加速度在相互垂直的三个坐标轴方向上的加速度分量；

根据所述加速度分量确定所述滚转角。

4.根据权利要求2所述的血流速度测量值的修正方法，其特征在于，所述姿态传感器包括陀螺仪；所述确定所述探头在检测状态下的滚转角，包括：

基于所述陀螺仪获得在检测状态下的所述滚转角。

5.根据权利要求1所述的血流速度测量值的修正方法，其特征在于，所述根据所述滚转角对所述第一晶元组血流速度测量值和所述第二晶元组血流速度测量值进行修正，得到第一修正值和第二修正值，包括采用以下公式得到第一修正值和第二修正值：

V1 = V / sin

；其中，V1表示第一修正值或第二修正值，V表示第一晶元组血流速度测量值或第二晶元组血流速度测量值，

表示滚转角。

6.根据权利要求1所述的血流速度测量值的修正方法，其特征在于，所述根据所述第一修正值和所述第二修正值确定待输出的血流速度测量值，包括：

将所述第一修正值和所述第二修正值的平均值确定为所述待输出的血流速度测量值。

7.根据权利要求1所述的血流速度测量值的修正方法，其特征在于，所述获取第一晶元组血流速度测量值和/或所述获取第二晶元组血流速度测量值，包括：

获取发射角，所述发射角为所述多普勒血流检测装置中发射晶元发射信号的方向与待检测血管中血流运动方向之间的夹角；

获取所述多普勒血流检测装置中位于同一晶元组中的所述发射晶元和接收晶元之间的距离以及所述发射晶元与所述待检测血管之间的距离；

根据所述发射角、所述发射晶元和接收晶元之间的距离、以及所述发射晶元与所述待检测血管之间的距离确定接收角，所述接收角为所述多普勒血流检测装置中接收晶元接收信号的方向与血流运动方向之间的夹角；

根据所述发射角和所述接收角确定所述第一晶元组血流速度测量值和/或所述第二晶元组血流速度测量值。

8.根据权利要求7所述的血流速度测量值的修正方法，其特征在于，所述发射晶元与所述待检测血管之间的距离包括第一距离，所述第一距离为沿所述发射晶元发射信号的方向上所述发射晶元到所述待检测血管的距离；

所述根据所述发射角、所述发射晶元和接收晶元之间的距离、以及所述发射晶元与所述待检测血管之间的距离确定接收角，包括采用以下公式确定接收角：

| β - α | = arctan ( D / L )；其中，β为接收角，α为发射角，| β - α | 为接收角与发射角之间的差值的绝对值，D为发射晶元和接收晶元之间的距离，L为第一距离。

9.根据权利要求8所述的血流速度测量值的修正方法，其特征在于，获取所述发射晶元与所述待检测血管之间的距离，包括：获取所述发射晶元发射信号至所述接收晶元接收所述信号所需的时长；根据所述所需的时长以及所述信号的传播速度确定所述信号在所述所需的时长内传播的距离；根据所述传播的距离确定所述第一距离。

10.一种多普勒血流检测装置，其特征在于，用于实现上述权利要求1至9中任意一项所述的血流速度测量值的修正方法。

Images