心血管动力学参数的检测方法.pdf

本发明为一种心血管动力学参数的无创伤检测方法。其步骤是，先采集寸口关部的脉搏波及血压信息，标定压力脉波，确定其中的10个特征点，提取14个特征量，由此计算获得心血管动力学参数，包括：脉率、左室射血时间、心室射血压、动脉内射流压、平均脉压、平均动脉压、心舒末容量、射血分数、每搏血量、心输出量、外周血管阻力、动脉顺应性等。这些参数可描述血压生理变化规律，并判断血压个性化分类，是高血压等疾病个性化诊断和治疗的辅助决策的重要客观依据。本发明方法的模型科学性强，与动物实验和人体实验符合程度高，参数敏感性强，显著性好。

1、  一种心血管动力学参数的检测方法，其特征在于具体步骤如下：
(1)用压力传感器按标准方法采集寸口关部脉搏波，同时准确采集相应部位的血压信息；
(2)将上述信息相结合和放大，低通滤波进行模数转换，标定为压力脉波P(t)，再作数据处理；
(3)用斜率法和模式识别法对压力脉波识别10个特征点，提取14个特征量；
(4)对心血管生理学参数进行计算；
其中，10个特征点为：
(1)a点：心房收缩开始点，
(2)b点：主动脉辨开放点，
(3)c点：主动脉最高压力点，
(4)d点：主动脉弹性振动点，
(5)e点：左心室停止射血点，
(6)e₁点：左心室舒张开始点，
(7)e₂点：主动脉辨关闭点，
(8)f点：二尖辨开放点，
(9)g点：反潮波起点，
(10)b’点：下一个周期的主动脉辨开放点。
14个特征量为：
(1)bc：射血时主动脉压上升时间，
(2)be₁：左心室快速射血时间，
(3)be₂：左心室射血时间，
(4)bf：左心室收缩时间，
(5)bg：动脉内主潮波时间，
(6)bb’：每搏周期时间，
(7)c高：主脉波高度，
(8)e₁高：潮波高度，
(9)e₂高：心舒点高度，
(10)f高：波谷点高度，
(11)g高：反潮波高度，
( 12 ) - - - ∫ b e 1 AP ( t ) · dt = APE 1 , ]]>射血期动脉压力积分，
( 13 ) - - - ∫ e 1 b ′ AP ( t ) · dt = APE 2 , ]]>心舒期动脉压力积分，
( 14 ) - - - ∫ b b ′ AP ( t ) · dt = APE 3 , ]]>每搏期动脉压力积分。

2、  根据权利要求1所述的检测方法，其特征在于
(1)脉率：PR＝60/bb’，单位为每分钟搏动次数               (1)
(2)左室射血时间：LVET＝be₂(ms)，单位为秒                (2)

3、  根据权利要求1所述的检测方法，其特征在于
(1)心室射血压：VEP＝DP+VPE+AEP                           (3)
其中，DP为舒张压，VPE为管道内存留压，
VPE＝e₁高×APE3/APE1×(1-e₁高/APE2                     (4)
AEP为动脉内射压，
AEP＝APE3×ATC                                           (5)
这里ATC为心室舒期后动脉内压稳定下降所表现的血液输运系数，
DP(t)＝e₁高×Exp(-ATC·t)
(2)平均脉压：MPP＝APE3/bb’单位为kPa                     (6)
(3)平均动脉压：MAP＝MPP+DP                               (7)
(4)动脉血容量：ABV＝SV×MAP/MPP单位为ml                  (8)
这里SV为每搏血量。

4、  根据权利要求1所述的检测方法，其特征在于
(1)心舒末容量：
EDV＝(-1/α)×Ln(1-β/K(0)×Ln(VEP×r))单位为毫升，      (9)
其中，VEP为心室射血压力，1/α＝80.76±7.27，β/K(0)＝0.146±0.06，r＝0.985/(1-Exp(vwk×t))，VWK＝0.378±0.037；
(2)射血分数：EF＝AEP/VEP                                 (10)
(3)每搏血量：SV＝EDV×EF单位L/min，                      (11)
(4)心输出量：CO＝SV×PR/1000单位L/min，                  (12)
(5)总外周血管阻力：TPR＝79.96×MAP/CO单位kPa·s/t，      (13)
(6)动脉顺应性：AC＝SV/(SP-DP)单位ml/kPa，                (14)。

心血管动力学参数的检测方法
技术领域
本发明属人体生理学参数检测技术领域，具体涉及一种心血管动力学参数的检测方法。
背景技术
贯彻个性化治疗原则，做到“量体裁衣”或“按锁配钥匙”，是提高高血压等疾病的疗效和控制率的关键措施。目前，个性化决策的常用方法，仍然是根据病史、体检、化验、心电、超声、X光等多方面的结果，凭专科医生的知识与经验进行决策。台湾邱绪新发明的DynaPluse200M型仪器，是在给患者肱部以袖袋式血压计测压时，拾取局部的震动信号来检测患者的心血管弹性参数，它不可能对人体整体的病理生理学状态进行系统分析。
本发明人曾在长期的心血管监测中发现：脉搏与血压和心血管疾病的关系密切，而且中医三指切脉也是采集脉搏与血压信息同步分析。血压相同的患者，其脉搏信息不同，则具有不同的病理生理学状态，可作为高血压等疾病个性化辨证论治的客观依据。发明人和合作者在脉诊现代化、数字化早期研究中，曾先后建立过三套脉诊参数算法，即以改进的弹性腔模型为基础的SS参数算法(1977年)和EZ参数算法(1980年)，二者在全国性办班推广已充分公开，并被改头换面使用(如安徽的ZXG型、XXG型心血管功能检测仪等)，这些仪器所用的参数，数学模型比较简单，没有考虑到人体神经调节、体液调节及心血管自身调节等因素，更没有经过动物实验和人体对照试验的验证与修正，所以不很符合人体病理生理学规律，难以用作患者病情判断和治疗的依据。
1986年后，发明人和合作者又改用状态空间模型为基础的ES参数算法，但其算法仍不能准确描述高血压等疾病的病理生理学变化规律，如按1999WHO/ISH标准分类，3级高血压的心输出量不应高于1级高血压。
发明内容
本发明的目的在于通过建立能真实反映血液循环动力学的数学模型，提出一种能对心血管动力学参数进行检测、分析的方法，为建立相应的无创、简便、快速检测的智能仪器打下基础。
本发明人长期在心脏直视手术围术期监测的过程中，研究了压力脉波与血液循环动力学状态的定量关系，发现有14项心血管动力学参数，可描述心血管生理学和病理生理变化规律，并判断其个性状态分类，以之作为高血压等病情个性化治疗辅助决策的客观依据。于是建立了这些心血管动力学参数无损伤检测和定量判断的方法，具体步骤如下：
1、用压力传感器按标准方法采集寸口关部脉搏波，同时准确采集相应部位的血压信息；
2、将上述信息相结合和放大，低通滤波后进行模数转换，标定为压力脉波P(t)，再作数据处理；
3、用斜率法和模式识别法对压力脉波识别10个特征点，提取14个特征量；
4、对心血管生理学参数进行计算。
下面就分析模型和计算方法作具体描述：
1、对压力脉波进行数据处理，识别特征点和提取特征量(图1)。
特征点10个：
(1)a点：心房收缩开始点，
(2)b点：主动脉辨开放点，
(3)c点：主动脉最高压力点，
(4)d点：主动脉弹性振动点，
(5)e点：左心室停止射血点，
(6)e₁点：左心室舒张开始点，
(7)e₂点：主动脉辨关闭点，
(8)f点：二尖辨开放点，
(9)g点：反潮波起点，
(10)b’点：下一个周期的主动脉辨开放点。
图1中，a-b波  (a’-b’波)为心房收缩波，m-n波为主动脉弹性振动波，L-a’波为主动脉静压排空段。
特征量：共14个
(1)bc：射血时主动脉压上升时间，
(2)be₁：左心室快速射血时间，
(3)be₂：左心室射血时间，
(4)bf：左心室收缩时间，
(5)bg：动脉内主潮波时间，
(6)bb’：每搏周期时间，
(7)c高：主脉波高度，
(8)e₁高：潮波高度，
(9)e₂高：心舒点高度，
(10)f高：波谷点高度，
(11)g高：反潮波高度，
( 12 ) - - - ∫ b e 1 AP ( t ) · dt = APE 1 , ]]>射血期动脉压力积分，
( 13 ) - - - ∫ e 1 b ′ AP ( t ) · dt = APE 2 , ]]>心舒期动脉压力积分，
( 14 ) - - - ∫ b b ′ AP ( t ) · dt = APE 3 , ]]>每搏期动脉压力积分。
其中，特征量(1)-(6)为时值，特征量(7)-(11)为幅值，特征量(12)-(14)为3个包络面积值。AP(t)为(Arterial Pressure)动脉压力函数。
根据图1测量所得的特征量可以计算下述参数：
(1)计算脉率(PR-Pulse Rate)，从心动周期bb’得：
PR＝60/bb’，bb’单位秒                                   (1)
(2)计算左室射血时间(LVET——L.V.EjectionTime)
LVET＝be₂(ms)，单位为秒                                   (2)
2、将脉搏波与血压相结合，分析各压力值(P)与流量值(V)的相互关系(图2)，
图2中，SP为收缩压，DP为舒张压，PP为脉压，VEP为心室射血压，AEP为(动脉内)射流压，ESV为心缩末血容量，EDV为心舒末血容量，SV为每搏血量，ABV为动脉血容量，MAP为平均动脉压，PEP为射血前期，LVET为左室射血时间，LVDT为左室舒张时间。
(1)计算心室射血压(VEP——Ventricular Ejection Pressure)
左室收缩全期造成的压力，可分解为舒张压(DP)，动脉内射流压(AEP)，管道内存留压(VPE)三部分，其计算式为：
VEP＝DP+VPE+AEP，单位为kPa                                 (3)
式中，
VPE＝e₁高×APE3/APE1×(1-e₁高/APE2)                        (4)
(2)计算动脉内射流压AEP——(Arterial ejection Pressure)。左室射血时，动脉内血流的压力一面上升并部分由动脉壁吸收为势能，一面使血液向末梢动脉输运，从脉搏波的起伏可以看出其动量与压力的变化过程。
AEP＝APE3×ATC                                             (5)
式中，ATC(Arterial Transpotation Coefficiency)为心室舒期后动脉内压稳定下降所表现的血液输运系数，如图1的L~a段，曲线为一下降的指数曲线。
其曲线方程为：
DP(t)＝e₁高×Exp(-ATC×t)
上式中，e₁后的压力积分和t可从压力脉波测量，DP(t)为心室舒张期动脉内压。
(3)计算平均脉压(MPP——Mean Pulse Pressure)
MPP＝APE3/bb’，单位为kPa                               (6)
(4)计算平均动脉压(MAP——Mean Arterial Pressure)
MAP＝MPP+DP，单位为kPa                                  (7)
(5)计算动脉血容量(ABV——Arterial Blood Volume)
ABV＝SV×MAP/MPP，单位为ml                              (8)
式中，SV为每搏血量(见后)。
3、按生物控制论中的状态空间法，建立与压力脉波特征量相应的微分方程，以描述心血管系统生理学状态。
如在心脏生理学中有著名的Starling定律(见周衍椒、张镜如主编：生理学，人民卫生出版社，1989)：在心率一定的条件下，心舒末左室容量EDV和心室射血压力VEP或每搏血量SV成函数关系(图3)。试将曲线放大，并忽略其下降支。则形成心室射血压力、容积、流量曲线(图4)，其曲线方程表达式为：
VEP(t)＝VEP(0)×(Exp(K(0)/α(1-Exp(-α×EDV)))
式中，VEP为心室射血压力，VEP(o)为压力的边界值，K(0)/α为心肌变力系数，EDV为心舒末(左室)容量。
通过一批动物实验和人体(心脏手术中)心导管压力曲线的描记和相关数据的提取，进行曲线拟合相关分析，其相关系数r＝0.99，表明所建曲线方程有高精度的可靠性。
下面根据这一曲线方程，对几个函数(或参数)求解。
(1)计算心舒末容量(EDV——End Diastolic Volume)，根据心脏手术前后的记录，得知心脏舒缩运动时体积变化(图5)及其方程表达式为：
V(t)＝V(0)×Exp(-VWK×t)
表述方程中；V(o)为体积缩小的初始量，即心舒末(左室)容量EDV。VWK为心室体积缩小曲线的时间常数，t为心缩时间(图2中“PEP+LVET”时间)。
根据动物(犬)试验和100例人体心脏手术中有关测值的统计，VWK＝0.378±0.037(sec-1)，1/α＝80.76±7.27，β/K(0)＝0.146±0.06。由于从压力脉波测得的VWK和t值可以求得：
r＝0.985/(1-Exp(vwk×t))。
根据图5中的曲线表述方程，移项可得：
EDV＝(-1/α)×Ln(1-β/K(0)×Ln(VEP×r))    单位为毫升，    (9)
在测量VEP之后，便可根据上述已知条件，求解EDV值。
(2)计算射血分数(EF——Ejection Fraction)已知一搏期间心室射血压为VEP，射血产生的动脉内射流压为AEP，二者的比值即射血分数EF。
EF＝AEP(kPa)/VEP(kPa)，                                    (10)
即：EF＝SV(ml)/EDV(ml)，
(3)计算每搏血量(SV——Stroke Volume)从方程组(14)移项，得SV＝EDV×EF单位ml，
(4)计算心输出量(CO——Cadiac Output)                       (11)
CO＝SV(ml)×PR/1000，单位L/min，                           (12)
(5)计算总外周血管阻力(TPR——Total Peripheral Resistance)
TPR＝79.96×MAP/CO  单位kPa·s/t，                         (13)
根据公认的知识，如果不知道右房压或由于右房压通常很低，计算时可以忽略，则可以通过上式(13)计算。
(6)计算动脉顺应性(AC——Arterial Compliance)。
动脉顺应性是指单位时间内动脉管段血液体积变量与压力变量的比值(即dv/dp)，故可表述为：
AC＝SV/(SP-DP)  单位ml/kPa，                               (14)
上述参数若按体重与体表面积计算，还可以求得反映个性化差别的各种指数。这些生理学参数具有如图6所示的相互关系，也表明了产生高血压的病理生理学机理，某些参数的异常特性，可表示患者的病理生理学个性。
附图说明
图1为典型脉波特征点和特征量测量图示。
图2为压力脉波波形与各压力值、流量值的关系示意图。
图2中，自a点左室开始收缩，经等容收缩期(或射血前期-PEP)使动脉中压力超过舒张压(b点)，压力从b点加速上升至c点，动脉升压因素和降压因素达到平衡。C点与e点地波动是由于主动脉壁弹性振动，e点后左心室舒张，并引起动脉中的反潮波(g波)，其后动脉中压力呈指数曲线下降。
图3为Frank-Starling关系曲线。
图3(A)表明EDV(心舒未容量)与VEP(心室射血压)的关系，图3(B)表明EDV与SV(每搏血量)的关系。
图4为左室压力、容积、流量关系曲线的上升段。
图5为心室舒缩运动时的体积变化曲线。
图6为所检测14项心血管动力学参数相关的关系图。
本发明提出的血液动力学参数的无损伤性检测方法是建立在严格的血液循环动力学基础之上的。所建立的分析、计算模型科学性强，与动物实验及人体试验符合程度高，所获得的病理生理学参数能较好地显示出高血压等疾病的病理生理学特性，参数的敏感性强，显著性好，可用以区分不同高血压等疾病患者的个性类型。从而有针对性提供合理而恰当的防治建议，对患者和专业医生均有实用价值。依据本发明方法提供的仪器(装置)，用嵌入式软硬件建造，结构简单，使用方便，可实现无损伤检测，所得结果可靠。实践表明、本发明有助于对高血压等疾病患者进行个性化的辨证论治和随访，研究药效与机理，进而可提高高血压等疾病的控制率。
具体实施方式
本发明采集人体左寸口关部(桡动脉)的体表脉搏信息，在标准体位与规范操作下，可取得典型的拐点清晰的脉波(图1)，结合用血压值标定为压力脉波，并可矫正予置压力的影响。由于采用超低频频响(0.01-50Hz)的压力传感器、线性放大器、低通滤波器，使失真度＜0.2％。得到图1和图2的曲线图。根据这2个图获得10个特征值和14个特征量。再根据计算式(1)-(14)，获得各项心血管动力学参数。
由本发明方法计算获得的动力学参数准确可靠。
(1)对所建立的数学模型分解的状态空间方程组进行曲线拟合相关分析，14个方程组的拟合相差系数均在0.99以上。证实了数学模型的可行性及其精度。在此基础上建立的XXG参数算法的理论误差甚小。
(2)经动物(犬)试验。用双盲法进行放血、输血等处理，证明本发明测定的血容量可靠(17次试验中16次一致，1次基本符合。也证明人工失血性休克及用抗休克液的结果和所监测的参数变化规律吻合。
(3)经动物(犬)试验和临床直接对照验证，结果小输出量(CO)测定的相关系数的心导管(Fick)法r＝0.805，热稀释法r＝0.868，色素法r＝0.891，电磁流量计法r＝0.953。心室射血压(VEP)测定54例，r＝0.963，总外周血管阻力(TPR)测定的结果，符合病理生理学和药物作用应有的变化。
本发明可以定量分析和判断一些心血管疾病的病理生理学状态，如对479例高血压患者的调查发现，患者虽在用药治疗，但查出种种参数的改变，有的属于2-3级高血压，根据特征性参数的异常，可将高血压区分为高动力型(VEP、VPE、SV升高)、高阻力型(TPR升高)、高容量型(ABV、EDV升高)，或高模量型(MPP升高、AC降低)，因而可用以辅助医生进行个体化治疗决策，也可用以随访和比较病情。如对49列高血压患者用美托洛尔(Betaloc)治疗前后中检测比较，有统计学意义的结果为MPP与MAP均明显降低的同时，VEP、CO及PR亦明显降低，而EDV、TPR、ABV、AC则无明显改变，这不仅定量阐明了该药降压作用的机理，也表明该药仅适用于有VEP、CO、PR升高的高动力型高血压。
根据上述病理生理学参数的检测和人工智能判断的方法，本发明设计了相应的专用检测仪器，该仪器包括检测系统，予处理系统和计算分析系统。其中检测系统用压力传感器配放大器采集信息，人工输入某些人体生理参数：性别、年龄、身高、体重等。予处理系统以单片机为核心，采用EC3-1451，586型或兼容CPU，配有模数转换器(12bitA/D)、ROM、RAM及与计算机系统连接的通讯接口。计算分析系统采用基于“和欣”操作系统的专用程序，该程序固化在专用ROM或FLASH中，另外，还配有显示器、打印机和电源等。
计算机分析系统包含有三个模块：①压力脉波自动识别取值模块(PRMVP程序)，②模型分解、参数计算模块(XXGP程序)③面向临床应用的专家咨询模块。AICEXP程序用C语言编程，包括①知识库，(知识规则与参数关系)；②数据库(正常标准与升降判定)；③推理机制与状态判式；