CN108937896B - 一种中心动脉血压测量方法及设备 - Google Patents
一种中心动脉血压测量方法及设备 Download PDFInfo
- Publication number
- CN108937896B CN108937896B CN201710369140.XA CN201710369140A CN108937896B CN 108937896 B CN108937896 B CN 108937896B CN 201710369140 A CN201710369140 A CN 201710369140A CN 108937896 B CN108937896 B CN 108937896B
- Authority
- CN
- China
- Prior art keywords
- blood pressure
- artery
- peripheral
- time
- central
- 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.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
- A61B5/0225—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers the pressure being controlled by electric signals, e.g. derived from Korotkoff sounds
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
- A61B5/02125—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Vascular Medicine (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Ophthalmology & Optometry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明实施例公开了一种中心动脉血压测量方法及设备,该方法包括:测量外周动脉的随时间变化的血压波形;获取血流从中心动脉到所述外周动脉的血压测量位置的传播时间;利用所述外周动脉的血压波形和所述传播时间计算所述中心动脉的血压波形。通过上述方式测量中心动脉血压波形计算量较少,且只需要测量单外周动脉的血压,测量难度低,使用方便,能够降低测量成本,且能够体现个体差异性。
Description
技术领域
本发明实施例涉及生物医学工程领域,特别是涉及一种中心动脉血压测量方法及设备。
背景技术
中心动脉血压(CAP)是独立于肱动脉血压的疾病预测因素,它对于高血压导致的动脉功能退化情况起着直接的预测作用,近年研究发现,中心动脉血压评价高血压治疗效果、疾病预后情况、动脉结构功能等方面都比传统手臂肱动脉血压测量更准确,其临床价值更高。
本发明的发明人在长期研究中发现,现有技术所采用的双外周中心动脉血压测量方法需要测量双外周动脉血压,计算量较大,使用不便且设备昂贵。
发明内容
本发明实施例主要解决的技术问题是提供一种中心动脉血压测量方法及设备,能够解决现有双外周中心动脉血压测量方法计算量大、使用不便的问题。
为解决上述技术问题,本发明实施例采用的一个技术方案是:提供一种中心动脉血压测量方法,包括:测量外周动脉的随时间变化的血压波形;获取血流从中心动脉到外周动脉的血压测量位置的传播时间;利用外周动脉的血压波形和传播时间计算中心动脉的血压波形。
为解决上述技术问题,本发明实施例采用的另一个技术方案是:提供一种中心动脉血压测量设备,包括:外周动脉测量仪,用于测量外周动脉随时间变化的血压波形;处理器,连接外周动脉测量仪,用于获取血流从中心动脉到外周动脉的血压测量位置的传播时间,利用外周动脉的血压波形和传播时间计算中心动脉的血压波形。
本发明实施例的有益效果是:在本发明实施例的方法中,通过测量外周动脉的随时间变化的血压波形,获取血流从中心动脉到外周动脉的血压测量位置的传播时间,最终利用外周动脉的血压波形和传播时间可计算得到中心动脉的血压波形,该方法中根据单外周动脉血压波形和血流传播时间即可计算得到中心动脉的血压波形,计算量较少,且只需要测量单外周动脉的血压,测量难度低,使用方便,且能够降低测量成本。
附图说明
图1是本发明中心动脉血压测量方法第一实施例的流程图;
图2是图1所示的方法所采用的单一动脉血管模型的示意图;
图3是手指血管关闭实验示意图;
图4是本发明中心动脉血压测量方法第二实施例的流程图;
图5是本发明中心动脉血压测量方法第三实施例的流程图;
图6是本发明中心动脉血压测量设备第一实施例的结构示意图;
图7是本发明中心动脉血压测量设备第二实施例的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,图1是本发明中心动脉血压测量方法第一实施例的流程图。如图1所示,本发明中心动脉血压测量方法包括:
步骤S11:测量外周动脉的随时间变化的血压波形;
具体地,基于血液动力学的基本原理,首先建立单一动脉血管模型20,如图2所示,假设中心动脉2012到外周动脉血管具有纯弹性无损的特性,并且时域独立的信号可以写成傅里叶序列,血管201各位置的平稳血流和平均血压相等,其特征阻抗以ZR表示。在血管末端,与之连接的是一个Windkessel模型202,用以描述远端血管子***特性,该子***可由一个等效特征阻抗Z(ω)表示。在血管末端,由于阻抗不匹配,动脉***存在反射现象,从心脏流向外周动脉的血流定义为前进波,而各处外周血管反射部分的叠加,由外周动脉流向心脏的血流定义为反射波。其中,前进波和反射波均为压力波,实际测得的血压波形是由前进波和反射波叠加而成的。
此外,人体在静息状态时,典型的手指血流仅为0.05ml/min,如此慢的血流速度表示单位时间血流量很少,暗示着手指动脉的终端阻力很大,则手指动脉血管可以近似看作为处于关闭状态,因此,手指动脉血压的前进波完全反射变成反射波,手指动脉测量的血压波形即为前进波的二倍。
为了验证手指动脉血管可以近似看作为处于关闭状态,如图3所示,在手指端设置两个袖带,其中袖带1连接手指动脉血压测量装置,用于测量手指动脉的血压波形(DigitalArterial Pressure,DAP),袖带2连接充气装置,用于在充气时给手指动脉血管施加压力,以使得手指动脉血管处于关闭状态。其中,采用容积补偿法测量该手指动脉的血压波形,当然,在其他实施例中,也可以采用脉搏波传导速度法等其他方法测量手指动脉的血压波形,此处不做具体限定。
通过测量30个志愿者手指动脉的血压波形,包括手指动脉血管处于非关闭状态和关闭状态的血压波形后,统计得到如下表1中的数据:
表1 实验结果
其中,n表示实验测量的血压波形中选取的样本数量,p-value是假定概率,表示对原假设的支持程度,通常采用显著性水平0.05与p-value比较,当p-value大于0.05时,支持原始假设,否则,拒绝原始假设。因此,有上述实验结果可知,p-value均大于0.05,即手指动脉血管处于非关闭状态和关闭状态时,手指动脉血压波形无明显变化,也就是说,手指动脉血管可以近似看作为处于关闭状态。
由此,在一个应用例中,基于如图2所示的单一动脉血管模型20,在手指端2011设置测量装置,例如血压测量袖带,测量随时间变化的手指动脉的血压波形,其中,手指动脉的血压波形测量方法可以采用容积补偿法。当然,在其他应用例中,手指动脉的血压波形也可以采用其他测量方法,若其他外周动脉的血管也可以近似看作处于关闭状态,也可以测量其他外周动脉的血压波形,此处不做具体限定。
其中,步骤S11进一步包括:
步骤S111:同步测量随时间变化的心电图信号;
具体地,在上述应用例中,在心脏部位2012设置心电测量仪,通过心电测量仪可以同步测量随时间变化的心电图信号。
步骤S12:获取血流从中心动脉到外周动脉的血压测量位置的传播时间;
其中,步骤S12进一步包括:
步骤S121:利用心电图信号和外周动脉的血压波形之间的相位差计算血流从中心动脉到外周动脉血压测量位置的传播时间。
具体地,心脏周期性地射血,血流沿着动脉流动产生周期性的血压脉搏波,这些波形随着离心脏距离的不同会有时间的延迟,即血流传播时间(PTT)。在一个应用例中,基于光电容积测量脉搏波传播时间的原理,在手指动脉部位2011设置一光电传感器,记录光电容积脉搏波(PPG),在心脏部位2012设置一心电测量仪,记录心电图信号(ECG),通过在一个周期中,计算ECG信号的R波波峰位置与同步测量的PPG信号之间的相位差,即可以得到血流从中心动脉到手指动脉血压测量位置的传播时间。当然,在其他应用例中,该传播时间也可以采用其他方法计算,此处不做具体限定。
步骤S13:利用手指动脉的血压波形和传播时间计算中心动脉的血压波形。
其中,步骤S13进一步包括:
步骤S131:利用外周动脉的血压波形分解外周动脉的血压测量位置的前进波和反射波;
其中,步骤S131进一步包括:
步骤S1311:通过如下公式分解出前进波和反射波:
其中,pDf(t)为外周动脉血压测量位置的前进波在t时刻的值,pDb(t)为外周动脉血压测量位置的反射波在t时刻的值,PD(t)为外周动脉血压测量位置在t时刻的血压。
具体地,基于图2所示的单一动脉血管模型,血流从中心动脉2012到手指动脉的血压测量位置2011之间的血管201中流动到远端血管202时,在血管201末端,由于阻抗不匹配,血流存在反射现象,从而出现前进波和反射波,同时,由于手指动脉血管可近似看作处于关闭状态,因此,在血管201末端,前进波完全反射为反射波,即如图2所示,在频域上,在血管201末端的反射系数
其中,Γ(ω)为外周动脉血压测量位置的反射系数在ω频点的值,Pf(ω)为外周动脉血压测量位置的前进波在ω频点的值,Pb(ω)为外周动脉血压测量位置的反射波在ω频点的值。同样的,时域上,该反射系数
因此,通过上述公式,可以利用测量到的手指动脉的血压波形以及步骤S1311计算出的手指动脉的血压测量位置的前进波和反射波。
步骤S132:利用手指动脉血压测量位置的前进波和反射波以及传播时间计算中心动脉的血压波形。
其中,步骤S132进一步包括:
步骤S1321:通过如下公式计算中心动脉的血压波形:
PC(t)=pDf(t+Δt)+pDb(t-Δt);
其中,PC(t)为中心动脉在t时刻的血压,pDf(t+Δt)和pDb(t-Δt)为外周动脉的血压检测位置的前进波和反射波分别在t+Δt和t-Δt时刻的值,Δt为血流从中心动脉到外周动脉的血压测量位置的传播时间。
具体地,基于图2所示的单一动脉血管模型,血管201中手指动脉的血压测量位置2011测得的前进波是由中心动脉2012的前进波经过Δt时间传播得到的,而中心动脉2012的反射波则是由手指动脉的血压测量位置2011测得的反射波经过Δt时间传播得到的,由此可以通过手指动脉的血压测量位置2011的前进波和反射波获得中心动脉2012的前进波和反射波,并且由于实际测得的血压波形是由前进波和反射波的压力叠加而成的,因此利用上述公式将手指动脉的血压测量位置2011的前进波和反射波变换后的得到的中心动脉2012的前进波和反射波相加,即可得到中心动脉的血压波形。
进一步参阅图2,在频域上,该单一动脉血管模型中,手指动脉的血压波形PD(ω)和中心动脉的血压波形PC(ω)之间的关系为
其中,H(ω)为中心动脉2012和手指动脉血压测量位置2011之间的传递函数在ω频点的值,PC(ω)为中心动脉血压在ω频点的值,PD(ω)为手指动脉血压测量位置的血压在ω频点的值。
此外,如图2所示,中心动脉2012的前进波和反射波与手指动脉的血压测量位置2011的前进波和反射波之间的关系,以及中心动脉2012的血压波形与中心动脉2012的前进波和反射波之间的关系分别为
PCf(ω)=Pf(ω)ejωΔt;
PCb(ω)=Pb(ω)e-jωΔt;
PC(ω)=PCf(ω)+PCb(ω);
其中,PC(ω)为中心动脉血压在ω频点的值,PCf(ω)为中心动脉血压的前进波在ω频点的值,PCb(ω)为中心动脉血压的发射波在ω频点的值,Pf(ω)为手指动脉血压测量位置的前进波在ω频点的值,Pb(ω)为手指动脉血压测量位置的反射波在ω频点的值,Δt为血流从中心动脉到手指动脉的血压测量位置的传播时间。
通过上述公式以及反射系数的公式可知,中心动脉2012和手指动脉血压测量位置2011之间的传递函数可以通过如下公式得到:
其中,H(ω)为中心动脉2012和手指动脉血压测量位置2011之间的传递函数在ω频点的值,Γ(ω)为手指动脉血压测量位置的反射系数在ω频点的值。
由上述步骤S1311可知,手指动脉血压测量位置的反射系数Γ(ω)的值为1,因此,上述公式可简化为
由此,中心动脉血压波形可以通过如下公式得到:
PC(ω)=PD(ω)cos(ωΔt);
其中,PC(ω)为中心动脉血压在ω频点的值,PD(ω)为手指动脉血压测量位置的血压在ω频点的值。
由此,在时域上,中心动脉血压波形则可以通过如下公式得到:
其中,PC(t)为中心动脉在t时刻的血压,PD(t+Δt)和PD(t-Δt)为手指动脉血压测量位置的血压分别在t+Δt和t-Δt时刻的值,Δt为血流从中心动脉到外周动脉的血压测量位置的传播时间。
上述实施例中,通过测量手指动脉随时间变化的血压波形,获取血流从中心动脉到手指动脉的血压测量位置的传播时间,最终利用手指动脉的血压波形和血流传播时间计算中心动脉的血压波形,只需根据单外周动脉血压波形和血流传播时间即可计算得到中心动脉的血压波形,计算量较少,且只需要测量单外周动脉的血压,测量难度低,使用方便,能够降低测量成本,并且由于不同人体实时测量的手指动脉的血压波形和血流传播时间均存在差异,该方法测量的中心动脉血压波形能够体现个体差异性。
当然,在其他实施例中,也可以测量其他外周动脉随时间变化的血压波形,此处不做具体限定。
请参阅图4,图4是本发明中心动脉血压测量方法第二实施例的流程图,本发明中心动脉血压测量方法第二实施例是在本发明中心动脉血压测量方法第一实施例的基础上,步骤S11进一步包括:
步骤S112:测量外周动脉关闭状态时的随时间变化的血压波形;
步骤S113:测量外周动脉非关闭状态时的随时间变化的血压波形;
由于外周动脉处于关闭状态时,外周动脉血压的前进波完全反射变成反射波,该外周动脉测量的血压波形即为前进波的二倍,并且利用该外周动脉处于非关闭状态时测量的血压波形和关闭状态时血压的一半,即可以计算得到该外周动脉血压的反射波。
具体地,在一个应用例中,基于如图2所示的单一动脉血管模型20,在外周动脉端2011设置测量装置,例如图3所示的两个血压测量袖带,测量外周动脉关闭状态和非关闭状态时的随时间变化的血压波形,其中,外周动脉的血压波形测量方法可以根据实际情况采用张力法、容积补偿法、脉搏波传导速度法等任意一种或其他测量方法,此处不做具体限定。
由此,进一步地,步骤S131包括:
步骤S1312:通过如下公式分解出外周动脉血压测量位置的前进波和反射波:
pDb(t)=PD2(t)-pDf(t);
其中,pDf(t)为外周动脉血压测量位置的前进波在t时刻的值,pDb(t)为外周动脉血压测量位置的反射波在t时刻的值,PD1(t)为外周动脉关闭状态时血压测量位置在t时刻的血压;PD2(t)为外周动脉非关闭状态时血压测量位置在t时刻的血压。
具体地,当该外周动脉血管不能近似看作处于关闭状态时,可以通过在外周动脉血压测量位置设置的测量装置,测量得到该外周动脉处于关闭状态和非关闭状态的血压波形,从而根据上述公式可以分解出该外周动脉血压测量位置的前进波和反射波,以便后续计算中心动脉的血压波形。当该外周动脉血管可以近似看作处于关闭状态时,可参考步骤S1311,此处不再重复。
本实施例可以与本发明中心动脉血压测量方法第一实施例相结合。通过上述方法,可以利用外周动脉血管处于关闭状态和非关闭状态时测量的血压波形,分解出前进波和反射波,以便后续计算中心动脉血压波形,该方法计算简单,测量方便,能够降低测量成本,并且由于不同人体实时测量的手指动脉的血压波形和血流传播时间均存在差异,该方法测量的中心动脉血压波形能够体现个体差异性。
请参阅图5,图5是本发明中心动脉血压测量方法第三实施例的流程示意图。如图5所示,本发明中心动脉血压测量方法第三实施例是在本发明中心动脉血压测量方法第一实施例的基础上,步骤S12包括:
步骤S122:利用先验知识估计血流从中心动脉到外周动脉的血压测量位置的传播时间。
具体地,在一个应用例中,针对每个测量对象,测量其外周动脉随时间变化的血压波形,提取舒张期的血压值,根据如图2所示的单一动脉血管模型及从心脏到该外周动脉血压测量位置的血管平均长度等先验知识,可以估计得到血流从中心动脉到外周动脉的血压测量位置的传播时间。根据多个测量对象的平均数据,即可得到最终估计的传播时间。当然,在其他应用例中,也可以提取外周动脉整个周期的血压波形或收缩期的血压值进行估计,此处不做具体限定。
本实施例的步骤可以替代步骤S121,本实施例可以与本发明中心动脉血压测量方法第一、第二实施例相结合。
请参阅图6,图6是本发明中心动脉血压测量设备第一实施例的结构示意图。如图6所示,本发明中心动脉血压测量设备60第一实施例的包括:外周动脉测量仪601和处理器602。
外周动脉测量仪601,用于测量外周动脉随时间变化的血压波形;
处理器602,连接外周动脉测量仪601,用于获取血流从中心动脉到外周动脉的血压测量位置的传播时间,利用外周动脉的血压波形和传播时间计算中心动脉的血压波形。
具体地,在一个应用例中,在待测试者手指部位设置外周动脉测量仪,例如图3所示的装置,利用容积补偿法测量出手指动脉随时间变化的血压波形,可以参考步骤S11,此处不再重复;然后将测量的手指动脉随时间变化的血压波形传输给处理器602,处理器602根据获得的血压波形,利用先验知识估计血流从中心动脉到外周动脉的血压测量位置的传播时间,具体可参考步骤S122,此处不再重复;最后利用外周动脉的血压波形和传播时间计算出中心动脉的血压波形,具体可以参考步骤S13,此处不再重复。
当然,在其他应用例中,可以测量其他外周动脉的血压波形,血压波形也可以通过其它方法测量,外周动脉测量仪也可以是其他装置,血流从中心动脉到外周动脉的血压测量位置的传播时间也可以通过其他方式获取,此处不做具体限定。
当然,在其他实施例中,中心动脉血压测量设备视具体需求还可以包括显示器、按键等,此处不做具体限定。
在上述实施例中,本发明中心动脉血压测量设备利用外周动脉测量仪测量外周动脉随时间变化的血压波形,利用处理器获取血流从中心动脉到外周动脉的血压测量位置的传播时间,最终利用外周动脉的血压波形和血流传播时间计算出中心动脉的血压波形,计算量较少,只需测量单外周动脉血压波形,测量难度低,使用方便,能够降低测量成本,并且由于不同人体实时测量的外周动脉的血压波形和血流传播时间均存在差异,因此该设备测量的中心动脉血压波形能够体现个体差异性。
请参阅图7,图7是本发明中心动脉血压测量设备第二实施例的结构示意图。图7与图6结构类似,此处不再赘述,不同之处在于,本发明中心动脉血压测量设备70进一步包括:心电测量仪703,连接于处理器702,用于同步测量随时间变化的心电图信号。
处理器702进一步用于利用心电图信号和外周动脉的血压波形之间的相位差计算血流从中心动脉到外周动脉血压测量位置的传播时间。
具体地,在一个应用例中,在待测试者的心脏部位设置心电测量仪703,在待测试者手指部位设置外周动脉测量仪701,同步测量随时间变化的心电图信号和外周动脉的血压波形,处理器703接收上述随时间变化的血压波形和心电图信号,基于光电容积测量脉搏波传播时间的原理,计算心电图信号和手指动脉的血压波形之间的相位差,从而获得血流从中心动脉到手指动脉血压测量位置的传输时间,具体计算方法可以参考步骤S121的方法,此处不再重复。当然,在其他应用例中,也可以测量其他外周动脉的血压波形,此处不做具体限定。
在上述实施例中,本发明中心动脉血压测量设备进一步利用测量的心电图信号和血压波形计算血流从中心动脉到外周动脉的血压测量位置的传播时间,最终利用外周动脉的血压波形和血流传播时间计算出中心动脉的血压波形,由于不同人体实时测量的外周动脉的血压波形和血流传播时间均存在差异,因此该设备测量的中心动脉血压波形能够体现个体差异性,并且计算量较少,只需测量单外周动脉血压波形,测量难度低,使用方便,能够降低测量成本。
以上所述仅为本发明的实施方式,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (2)
1.一种中心动脉血压测量设备,其特征在于,所述设备包括:
外周动脉测量仪,用于测量外周动脉关闭状态时以及非关闭状态时随时间变化的血压波形;
处理器,连接所述外周动脉测量仪,所述处理器利用先验知识估计血流从中心动脉到所述外周动脉的血压测量位置的传播时间;
利用所述外周动脉的血压波形分解出所述外周动脉的血压测量位置的前进波和反射波;
当所述外周动脉处于关闭状态时,通过如下公式分解出所述前进波和反射波:
其中,pDf(t)为所述外周动脉血压测量位置的前进波在t时刻的值,pDb(t)为所述外周动脉血压测量位置的反射波在t时刻的值,PD(t)为所述外周动脉血压测量位置在t时刻的血压;
当所述外周动脉处于非关闭状态时,通过如下公式分解出所述前进波和反射波:
pDb(t)=PD2(t)-pDf(t);
其中,pDf(t)为所述外周动脉血压测量位置的前进波在t时刻的值,pDb(t)为所述外周动脉血压测量位置的反射波在t时刻的值,PD1(t)为所述外周动脉关闭状态时血压测量位置在t时刻的血压;PD2(t)为所述外周动脉非关闭状态时血压测量位置在t时刻的血压;利用所述外周动脉血压测量位置的前进波和反射波以及所述传播时间计算所述中心动脉的血压波形。
2.根据权利要求1的设备,其特征在于,进一步包括:心电测量仪,用于同步测量随时间变化的心电图信号;
所述处理器连接所述心电测量仪,进一步用于利用所述心电图信号和所述外周动脉的血压波形之间的相位差计算血流从所述中心动脉到所述动脉血压测量位置的传播时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710369140.XA CN108937896B (zh) | 2017-05-23 | 2017-05-23 | 一种中心动脉血压测量方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710369140.XA CN108937896B (zh) | 2017-05-23 | 2017-05-23 | 一种中心动脉血压测量方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108937896A CN108937896A (zh) | 2018-12-07 |
CN108937896B true CN108937896B (zh) | 2021-08-03 |
Family
ID=64493672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710369140.XA Active CN108937896B (zh) | 2017-05-23 | 2017-05-23 | 一种中心动脉血压测量方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108937896B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113143230B (zh) * | 2021-05-11 | 2022-05-20 | 重庆理工大学 | 一种外周动脉血压波形重构*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103479343A (zh) * | 2013-09-27 | 2014-01-01 | 上海交通大学 | 基于振荡式血压计信号的中心动脉压检测***及方法 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL124787A0 (en) * | 1998-06-07 | 1999-01-26 | Itamar Medical C M 1997 Ltd | Pressure applicator devices particularly useful for non-invasive detection of medical conditions |
JP4470868B2 (ja) * | 2005-11-22 | 2010-06-02 | パナソニック電工株式会社 | 生体情報計測装置 |
US8313439B2 (en) * | 2009-03-20 | 2012-11-20 | Massachusetts Institute Of Technology | Calibration of pulse transit time measurements to arterial blood pressure using external arterial pressure applied along the pulse transit path |
KR101069489B1 (ko) * | 2009-04-30 | 2011-09-30 | 연세대학교 산학협력단 | 말초 혈관 및 혈류상태 진단장치 및 방법 |
US20110021929A1 (en) * | 2009-07-27 | 2011-01-27 | Nellcor Puritan Bennett Ireland | Systems and methods for continuous non-invasive blood pressure monitoring |
CN102197996B (zh) * | 2010-11-23 | 2013-04-17 | 中国科学院合肥物质科学研究院 | 利用桡动脉脉搏波获取主动脉脉搏波速度的方法 |
SG182872A1 (en) * | 2011-01-21 | 2012-08-30 | Healthstats Int Pte Ltd | Apparatus and method for altering the arterial pulse waveform of a body |
CN102499658B (zh) * | 2011-11-08 | 2013-12-11 | 中国科学院深圳先进技术研究院 | 中心血压波形重建模型及重建方法 |
CN103284703B (zh) * | 2013-06-18 | 2015-05-13 | 中国科学院合肥物质科学研究院 | 基于上肢动脉信息的主动脉脉搏波传递时间测量方法 |
CN104414679A (zh) * | 2013-09-04 | 2015-03-18 | 北京航空航天大学 | 一种中心动脉压重建中外周压力波延迟时间的确定方法 |
TWI535416B (zh) * | 2015-06-02 | 2016-06-01 | 國立中央大學 | 非侵入且非加壓式血壓波量測裝置及方法 |
WO2017024457A1 (zh) * | 2015-08-08 | 2017-02-16 | 深圳先进技术研究院 | 连续血压测量装置、测量模型建立方法和*** |
CN105595979A (zh) * | 2016-01-21 | 2016-05-25 | 中山大学 | 一种基于脉搏波传播时间的无创连续血压监测方法及装置 |
CN105725999A (zh) * | 2016-01-26 | 2016-07-06 | 四川长虹电器股份有限公司 | 一种基于脉搏波传导时间的血压测量方法 |
CN105748051B (zh) * | 2016-02-18 | 2018-10-09 | 京东方科技集团股份有限公司 | 一种血压测量装置 |
CN106343990A (zh) * | 2016-09-21 | 2017-01-25 | 中国矿业大学 | 一种指压式血压计及其测量方法 |
CN106333663A (zh) * | 2016-10-20 | 2017-01-18 | 深圳欧德蒙科技有限公司 | 一种血压监测方法和装置 |
CN106618537B (zh) * | 2016-12-21 | 2020-09-01 | 天津普仁万合信息技术有限公司 | 一种基于脉搏波传导的连续动态血压监测装置和方法 |
CN108236459A (zh) * | 2016-12-26 | 2018-07-03 | 深圳先进技术研究院 | 一种中心动脉血压测量方法及*** |
-
2017
- 2017-05-23 CN CN201710369140.XA patent/CN108937896B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103479343A (zh) * | 2013-09-27 | 2014-01-01 | 上海交通大学 | 基于振荡式血压计信号的中心动脉压检测***及方法 |
Non-Patent Citations (2)
Title |
---|
基于两个外周血压重建中心动脉压的方法及***关键技术的研究;蔡优飞;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140315;I140-649 * |
蔡优飞.基于两个外周血压重建中心动脉压的方法及***关键技术的研究.《中国优秀硕士学位论文全文数据库 信息科技辑》.2014,I140-649. * |
Also Published As
Publication number | Publication date |
---|---|
CN108937896A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10349838B2 (en) | Methods and apparatus for determining arterial pulse wave velocity | |
Kurylyak et al. | A Neural Network-based method for continuous blood pressure estimation from a PPG signal | |
US10709424B2 (en) | Method and system for cuff-less blood pressure (BP) measurement of a subject | |
US10398323B2 (en) | Methods and apparatus for determining pulse transit time as a function of blood pressure | |
EP2939592B1 (en) | Detection of parameters in cardiac output related waveforms | |
US10136823B2 (en) | Methods and apparatus for determining cuff blood pressure | |
Hahn et al. | Subject-specific estimation of central aortic blood pressure using an individualized transfer function: a preliminary feasibility study | |
JP2002253519A5 (zh) | ||
Arza et al. | Pulse transit time and pulse width as potential measure for estimating beat-to-beat systolic and diastolic blood pressure | |
CN107785081B (zh) | 计算中心血流动力学指标的方法、装置、存储介质及设备 | |
TW200307806A (en) | Systems and methods for calibrating a distorted signal with another signal of known calibration | |
Karamanoglu et al. | On-line synthesis of the human ascending aortic pressure pulse from the finger pulse | |
CN114652351B (zh) | 基于超声多普勒的连续血压测量方法、装置和电子设备 | |
US20160058409A1 (en) | Ultrasonic blood pressure measurement apparatus and blood pressure measurement method | |
US20040088123A1 (en) | Method for real-time monitoring of cardiac output and blood flow in arteries and apparatus for implementing the same | |
US11723543B2 (en) | Non-invasive system and method for measuring blood pressure variability | |
CN108937896B (zh) | 一种中心动脉血压测量方法及设备 | |
Zhang et al. | Pulse transit time-based blood pressure estimation using hilbert-huang transform | |
Nabeel et al. | Deep learning for blood pressure estimation: an approach using local measure of arterial dual diameter waveforms | |
CN105726001B (zh) | 一种血压测量方法及装置 | |
RU2644299C1 (ru) | Осциллографический способ измерения артериального давления | |
CN108236458A (zh) | 基于双外周血压重建肱动脉血压的方法 | |
Hassan et al. | Measuring blood pressure using a photoplethysmography approach | |
Wang et al. | An improved algorithm for noninvasive blood pressure measurement | |
Xu et al. | A calibration method for cuffless continue blood pressure measurement using Gaussian normalized pulse transit time |
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 |