CN113823409A - 一种透析中低血压事件风险的评估方法和*** - Google Patents
一种透析中低血压事件风险的评估方法和*** Download PDFInfo
- Publication number
- CN113823409A CN113823409A CN202111112325.5A CN202111112325A CN113823409A CN 113823409 A CN113823409 A CN 113823409A CN 202111112325 A CN202111112325 A CN 202111112325A CN 113823409 A CN113823409 A CN 113823409A
- Authority
- CN
- China
- Prior art keywords
- rbv
- data
- rate
- threshold
- risk
- 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
Links
- 208000001953 Hypotension Diseases 0.000 title claims abstract description 31
- 238000000502 dialysis Methods 0.000 title claims abstract description 23
- 208000021822 hypotensive Diseases 0.000 title claims abstract description 21
- 230000001077 hypotensive effect Effects 0.000 title claims abstract description 21
- 238000011156 evaluation Methods 0.000 title claims description 25
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000001914 filtration Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 14
- 230000009467 reduction Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 33
- 230000007423 decrease Effects 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000011157 data evaluation Methods 0.000 claims description 14
- 230000007704 transition Effects 0.000 claims description 12
- 238000005094 computer simulation Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 230000036543 hypotension Effects 0.000 abstract description 7
- 238000012545 processing Methods 0.000 abstract description 4
- 239000008280 blood Substances 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 5
- 230000003247 decreasing effect Effects 0.000 description 5
- 238000012502 risk assessment Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 208000012866 low blood pressure Diseases 0.000 description 3
- 238000000108 ultra-filtration Methods 0.000 description 3
- 230000036772 blood pressure Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- 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
-
- 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/026—Measuring blood flow
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Public Health (AREA)
- Data Mining & Analysis (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Mathematical Optimization (AREA)
- Cardiology (AREA)
- Veterinary Medicine (AREA)
- Databases & Information Systems (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Hematology (AREA)
- Operations Research (AREA)
- Vascular Medicine (AREA)
- External Artificial Organs (AREA)
Abstract
本发明公开了一种透析中低血压事件风险的评估方法和***。确定理想RBV曲线的一阶导数曲线。以预定间隔时长采集透析过程中的RBV数据,每次采集的数据使用卡尔曼滤波器进行滤波,计算RBV数据的下降率,使用卡尔曼滤波器对该下降率进行滤波。通过一阶导数曲线计算当前时刻的安全线和禁止线,根据比对下降率与安全线和禁止线的大小,确定风险系数。本发明能滤除采集的RBV数据中的干扰,提高使用RBV评估患者透析过程中低血压事件风险的准确性。同时,使用动态阈值进行低血压事件风险的计算,使其能够更好地适应不同时间段患者对失衡程度的耐受度变化。滤波算法简单,处理速度快,可以实现对RBV数据的实时处理。
Description
技术领域
本发明涉及医疗管理技术领域,尤其是一种透析中低血压事件风险的评估方法和***。
背景技术
透析中低血压是慢性肾脏病患者透析治疗中最常见的并发症,出现此症状的主要原因之一是透析中的超滤速率与患者体内的再注入速率失衡导致患者体内核心区域的血容量迅速减少,超出患者的耐受能力。通过分析患者在治疗中的相对血容量(RBV)变化情况,可以推断出患者体内的失衡情况,从而评估患者发生低血压事件的风险。此风险值可应用于透析过程中的生物反馈闭环控制,可以提前干预患者体内的失衡情况,降低患者出现低血压事件的几率。具体地,通常通过RBV的绝对值、变化量、变化率(一阶导数)、二阶导数是否超过某一个或多个阈值来判断患者低血压事件是否或发生者发生的风险程度。
CN110151153A公开了一种现有技术方案,其通过每5分钟采集患者RBV数据,对RBV的二阶导数SDRBV进行阈值判断,并可选地结合血压数据,从而判断患者是否会发生低血压事件。另外,在该文献中,还采用了另一种模糊逻辑技术,使用至少两个模糊模块接收血液动力学参数(例如相对血容量和血压),通过对各个模糊模块的输出进行加权,最终输出至少一个变量,用以评价患者发生低血压事件的危险程度。
CN104346521A公开了另一种现有技术,其采用利用学习算法或神经网络对患者个体透析参数进行学习后对治疗中患者的血液动力学参数进行预测的装置和方法。
现有技术存在以下的不足:
1.RBV的噪声会很大程度干扰RBV的导数变化情况,从而降低对低血压事件判断的准确性。
2.未考虑患者在治疗中对失衡程度的耐受度变化。
3.非实时,可能错过最佳干预时机。
4.判断结果为逻辑值,不利于后续精确的反馈控制。
发明内容
本发明的发明目的在于:针对上述存在的全部或部分问题,提供一种透析中低血压事件风险的评估方法,以及一种透析中低血压事件风险的评估***,以提高使用RBV评估患者透析过程中低血压事件风险准确性。
本发明采用的技术方案如下:
一种透析中低血压事件风险的评估方法,包括:
以预定间隔时长采集透析过程中的RBV数据,对于每次采集的RBV数据,执行以下流程:
使用预构建的卡尔曼滤波器对所述RBV数据进行滤波;计算观测到的RBV数据的下降率,使用所述卡尔曼滤波器对计算的所述下降率进行滤波;
基于对当前时刻设定的评估门限,对经滤波后的所述下降率进行评估。
进一步的,所述对当前时刻设定的评估门限,包括第一门限、以及高于所述第一门限的第二门限。
进一步的,所述基于对当前时刻设定的评估门限,对经滤波后的所述下降率进行评估,为基于所述第一门限和第二门限,计算经滤波后的所述下降率的风险系数。
进一步的,所述基于对当前时刻设定的评估门限,对经滤波后的所述下降率进行评估,包括:
在经滤波后的所述下降率处于所述第一门限与第二门限之间时,采用线性插值法或单调递增曲线对经滤波后的所述下降率进行评估。
进一步的,所述第一门限为设定的理想RBV曲线的一阶导数在当前时刻的值,所述第二门限为所述第一门限的N倍,N为大于1的常数。
进一步的,所述理想RBV曲线的一阶导数在第i时刻的值FDRBVi的计算方法为:
FDRBVi=n×EDRBV×ai,
式中,DRBV为预期RBV下降量,T为治疗总时间,n为放大倍数,a为衰减参数。
进一步的,所述卡尔曼滤波器的***状态预测方程为:
更新方程为:
式中,Xk为当前时刻的***状态,Xk-1为上一时刻的***状态,Zk为当前时刻观测到的***状态,A为状态转移矩阵,H为观测矩阵,Q、R分别表示***预测模型和传感器的噪声协方差矩阵,Pk为Xk的协方差矩阵,Pk-1为Xk-1的协方差矩阵,I为单位矩阵,Kk为卡尔曼增益矩阵,是一个中间计算量,带有“^”表示是估计值,带有“-”表示是预测值。
进一步的,所述状态转移矩阵A的计算方法包括:
构建RBV在t时刻的计算模型;
将所述计算模型离散化表示;
计算RBV的变化速率;
使用离散化的计算模型和所述变化速率表示***状态,与***状态标准表示方式类比,得到状态转移矩阵。
进一步的,所述观测矩阵H的计算方法包括:
计算观测到的RBV数据的下降率;
使用离散化的计算模型和所述下降率表示观测到的***状态,与观测到的***状态标准表示方式类比,得到观测矩阵。
本发明提供的一种透析中低血压事件风险的评估***,包括数据采集部、数据评估部和数据输出部,所述数据评估部预构建有卡尔曼滤波器,设定有对应各时刻的评估门限;
所述数据采集部,根据配置的采集周期,采集透析过程中的RBV数据;
所述数据评估部,使用所述卡尔曼滤波器对所述RBV数据进行滤波,计算观测到的RBV数据的下降率,使用所述卡尔曼滤波器对计算的所述下降率进行滤波;基于当前时刻的评估门限,对经滤波后的所述下降率进行评估;
所述数据输出部,输出所述数据评估部的评估结果。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1、本发明采用卡尔曼滤波对相对血容量数据进行预处理,利用简单的算法滤除干扰信号的影响,使其更好地反应出患者实际的相对血容量变化情况,提高了使用RBV评估患者透析过程中低血压事件风险的准确性。并且,由于卡尔曼滤波算法简单,处理速度快,可以实现对RBV数据的实时处理,减少修正措施的反应时间。
2、本发明使用动态阈值进行低血压事件风险的计算,使其能够更好地适应不同时间段患者对失衡程度的耐受度变化。
3、本发明方法和***最终的评估结果为对下降率评估得到的归一化风险系数,对风险等级判断更加精细,方便后续的反馈控制或修正措施的参数确定。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是本发评估方法的一个实施方式的流程图。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书(包括任何附加权利要求、摘要)中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
实施例一
如图1所示,透析中低血压事件风险的评估方法包括以下流程:
第一步、确定透析过程的理想RBV曲线的一阶导数(FDRBV)曲线
首先确定本次治疗预期RBV下降量(DRBV),以及治疗总时长T,默认DRBV设置为24%,治疗总时长T为4小时。当然,两个数据可根据患者实际情况确定。则整个治疗过程中的RBV平均下降率(EDRBV)为:
取EDRBV的5倍作为初始FDRBV,治疗中第i时刻的FDRBV表示为FDRBVi,则FDRBVi使用下式计算:
FDRBVi=5×EDRBV×0.985i
第二步、建立卡尔曼滤波器方程和确定滤波参数
卡尔曼滤波是一种通过***输入输出观测数据对***状态进行最优估计的算法,观测数据中包括***中的噪声和干扰的影响,所以最优估计也可看作是滤波过程。卡尔曼滤波器(下文简称滤波器)的工作依据两个方程:***状态转移方程和***观测方程。
以Xk表示当前时刻的***状态,***的状态转移方程可以表示为:
Xk=AXk-1+Buk+ωk
式中,A为状态转移矩阵;Xk-1表示上一时刻的***状态;B为控制矩阵;uk表示控制向量;ωk~N(0,Qk)表示过程激励噪声,是均值为0,方差为Qk的高斯白噪声。
同时,以Zk表示使用测量***观测到的当前时刻的***状态,其可表示为:
Zk=HXk+vk
式中,H称为观测矩阵;vk~N(0,Rk),表示观测噪声,是均值为0,方差为Rk的高斯白噪声。
透析治疗期间患者的RBV由于超滤是不断变化(下降)的,记为rbv,下降速率记为m。那么t时刻的rbv可以表示为:
rbvt=rbvt-Δt+mt-Δt×Δt
式中,Δt表示变化时长,rbvt-Δt表示t-Δt时刻的rbv,mt-Δt表示t-Δt的变化速率。
将上式的rbv离散化表示为:
rbvk=rbvk-1+mk-1
从收集的rbv数据分析得出,在没有意外干扰(患者稳定治疗)的情况下,下降速率m在治疗期间是不断减小的,经验公式为:
mk=a×mk-1
则***的预测方程可以归纳为:
将***的状态统一表示为Xk,即:
那么有:
即:
B=0
a的取值在本实施例中为0.985,同第一步中计算FDRBVi公式中的参数。
对于直接观测到的治疗中的rbv,记为rbv′。定义rbv′的下降率m′为:
其中,C为一常数,linearFit(C)表示对最近的C个rbv′数据进行线性拟合得出的斜率。C的取值决定了对rbv′的短期还是长期变化的估计,在本实施例中,需要对短期变化进行估计,则C取较小值,具体取值取决于rbv′数据的采集周期,在周期为1min的情况下,C取5。
则观测到的***状态统一表示为Zk:
由于观测到的参数rbv′、m′与***状态rbv、m是相同的量,无需进行转化,所以Zk又可以表示为:
也即:
滤波需要用到以下5个公式:
***状态预测方程:
式中,带有“^”表示是估计值,带有“-”表示是预测值(下同)。
更新方程:
式中,Pk为Xk的协方差矩阵,Pk-1为Xk-1的协方差矩阵,Q、R分别表示***预测模型和传感器的噪声协方差矩阵,需要根据实际情况确定,Kk为卡尔曼增益矩阵,是中间计算量,I为单位矩阵。在本实施例中,优选取值为:
第三步、采集RBV数据,进行滤波
以预定时间间隔(如每分钟)采集一次RBV数据。
将采集的RBV数据作为滤波器的第一个输入;并依照第二步中的相关算法得出斜率(linearFit(C)),作为滤波器的第二个输入,对RBV数据及其斜率进行滤波,去除噪声,输出滤波后的RBV和斜率。当然,可以仅对斜率进行滤波。
第四步、生成低血压风险评估系数
根据第一步中的方法计算出当前时刻的FDRBVi,将其作为RBV斜率的安全线,并将其进行一定程度的放大(例如n×FDRBVi,n为大于1的正整数),作为RBV斜率的禁止线。基于设定的安全线及禁止线,对采集的RBV数据进行风险评估。具体的,为根据滤波后的RVB的斜率与安全线和禁止线间的大小关系进行风险评估。
1.如果滤波后的RVB的斜率超过FDRBVi(即安全线),则认为极有可能已(或将)发生低血压事件,其风险系数记为1。
2.如果滤波后的RBV的斜率小于FDRBVi,则认为没有发生低血压事件的风险,其风险系数记为0。
3.如果滤波后的RBV的斜率在FDRBVi与n×FDRBVi之间,则采用线性插值法或其它单调递增曲线计算风险系数。
第五部、判断治疗过程是否结束,若是,则结束评估,否则,跳转到第三步,继续对下一次采集的RBV数据进行评估。
实施例二
本实施例公开了一种透析中低血压事件风险的评估***,包括数据采集部、数据评估部和数据输出部。数据评估部预构建有卡尔曼滤波器,并且设定有对应各时刻的评估门限。
评估门限是基于理想RBV曲线的一阶导数曲线(FDRBV)设定的。首先确定本次治疗预期RBV下降量(DRBV),以及治疗总时长T,默认DRBV设置为24%,治疗总时长T为4小时。当然,两个数据可根据患者实际情况确定。则整个治疗过程中的RBV平均下降率(EDRBV)为:
取EDRBV的n倍作为初始FDRBV,治疗中第i时刻的FDRBV表示为FDRBVi,则FDRBVi使用下式计算:
FDRBVi=n×EDRBV×ai;
式中,n为放大倍数,a为衰减参数,a<1;在本实例中n取5,a取0.985。
评估阈值包括安全线和禁止线,以当前时刻的FDRBVi作为RBV斜率的安全线,并将其进行一定程度的放大(例如N×FDRBVi,N为大于1的常数),作为RBV斜率的禁止线。
预构建的卡尔曼滤波器具体如下:
卡尔曼滤波器(下文简称滤波器)的工作依据两个方程:***状态转移方程和***观测方程。
以Xk表示***状态,***的状态转移方程可以表示为:
Xk=AXk-1+Buk+ωk
式中,A为状态转移矩阵;Xk-1表示上一时刻的***状态;B为控制矩阵;uk表示控制向量;ωk~N(0,Qk)表示过程激励噪声,是均值为0,方差为Qk的高斯白噪声。
同时,以Zk表示使用测量***观测到的***当前状态,其可表示为:
Zk=HXk+vk
式中,H称为观测矩阵;vk~N(0,Rk),表示观测噪声,是均值为0,方差为Rk的高斯白噪声。
透析治疗期间患者的RBV由于超滤是不断变化(下降)的,记为rbv,下降速率记为m。那么t时刻的rbv可以表示为:
rbvt=rbvt-Δt+mt-Δt×Δt
式中,Δt表示变化时长,rbvt-Δt表示t-Δt时刻的rbv,mt-Δt表示t-Δt的变化速率。
将上式的rbv离散化表示为:
rbvk=rbvk-1+mk-1
从收集的rbv数据分析得出,在没有意外干扰(患者稳定治疗)的情况下,下降速率m在治疗期间是不断减小的,经验公式为:
mk=a×mk-1
则***的预测方程可以归纳为:
将***的状态统一表示为Xk,即:
那么有:
即:
B=0
a的取值在本方案中为0.985,同第一步中计算FDRBVi公式中的参数。
对于直接观测到的治疗中的rbv,记为rbv′。定义rbv′的下降率m′为:
其中,C为一常数,linearFit(C)表示对最近的C个rbv′数据进行线性拟合得出的斜率。C的取值决定了对rbv′的短期还是长期变化的估计,在本实施例中,需要对短期变化进行估计,则C取较小值,具体取值取决于rbv′数据的采集周期,在周期为1min的情况下,C取5。
则观测到的***状态统一表示为Zk:
由于观测到的参数rbv′、m′与***状态rbv、m是相同的量,无需进行转化,所以Zk又可以表示为:
也即:
滤波需要用到以下5个公式:
***状态预测方程:
式中,带有“^”表示是估计值(下同);带有“-”表示是预测值(下同)。
更新方程:
式中,Q、R分别表示***预测模型和传感器的噪声协方差矩阵,需要根据实际情况确定,Kk为卡尔曼增益矩阵,是中间计算量,I为单位矩阵。在本实施例中,优选取值为:
数据采集部以预定时间间隔(如每分钟)采集一次RBV数据,每次采集的RBV数据均传递给数据评估部。
数据评估部将接收的RBV数据作为滤波器的第一个输入;并依照第二步中的相关算法得出斜率(linearFit(C)),作为滤波器的第二个输入,对RBV数据及其斜率进行滤波,去除噪声,输出滤波后的RBV和斜率。当然,可以仅对斜率进行滤波。
对于经滤波后的斜率,根据其与安全线和禁止线间的大小关系进行风险评估:
1.如果滤波后的RVB的斜率超过FDRBVi(即安全线),则认为极有可能已(或将)发生低血压事件,其风险系数记为1。
2.如果滤波后的RBV的斜率小于FDRBVi,则认为没有发生低血压事件的风险,其风险系数记为0。
3.如果滤波后的RBV的斜率在FDRBVi与n×FDRBVi之间,则采用线性插值法或其它单调递增曲线计算风险系数。
数据输出部连接数据评估部,输出数据评估部每次对数据采集部采集的RBV数据的评估结果(即风险系数)。数据输出部输出的数据可作为反馈控制或修正措施的参数。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。
Claims (10)
1.一种透析中低血压事件风险的评估方法,其特征在于,包括:
以预定间隔时长采集透析过程中的RBV数据,对于每次采集的RBV数据,执行以下流程:
使用预构建的卡尔曼滤波器对所述RBV数据进行滤波;计算观测到的RBV数据的下降率,使用所述卡尔曼滤波器对计算的所述下降率进行滤波;
基于对当前时刻设定的评估门限,对经滤波后的所述下降率进行评估。
2.如权利要求1所述的透析中低血压事件风险的评估方法,其特征在于,所述对当前时刻设定的评估门限,包括第一门限、以及高于所述第一门限的第二门限。
3.如权利要求2所述的透析中低血压事件风险的评估方法,其特征在于,所述基于对当前时刻设定的评估门限,对经滤波后的所述下降率进行评估,为基于所述第一门限和第二门限,计算经滤波后的所述下降率的风险系数。
4.如权利要求3所述的透析中低血压事件风险的评估方法,其特征在于,所述基于对当前时刻设定的评估门限,对经滤波后的所述下降率进行评估,包括:
在经滤波后的所述下降率处于所述第一门限与第二门限之间时,采用线性插值法或单调递增曲线对经滤波后的所述下降率进行评估。
5.如权利要求2~4任一所述的透析中低血压事件风险的评估方法,其特征在于,所述第一门限为设定的理想RBV曲线的一阶导数在当前时刻的值,所述第二门限为所述第一门限的N倍,N为大于1的常数。
8.如权利要求7所述的透析中低血压事件风险的评估方法,其特征在于,所述状态转移矩阵A的计算方法包括:
构建RBV在t时刻的计算模型;
将所述计算模型离散化表示;
计算RBV的变化速率;
使用离散化的计算模型和所述变化速率表示***状态,与***状态标准表示方式类别,得到状态转移矩阵。
9.如权利要求8所述的透析中低血压事件风险的评估方法,其特征在于,所述观测矩阵H的计算方法包括:
计算观测到的RBV数据的下降率;
使用离散化的计算模型和所述下降率表示观测到的***状态,与观测到的***状态标准表示方式类比,得到观测矩阵。
10.一种透析中低血压事件风险的评估***,其特征在于,包括数据采集部、数据评估部和数据输出部,所述数据评估部预构建有卡尔曼滤波器,设定有对应各时刻的评估门限;
所述数据采集部,根据配置的采集周期,采集透析过程中的RBV数据;
所述数据评估部,使用所述卡尔曼滤波器对所述RBV数据进行滤波,计算观测到的RBV数据的下降率,使用所述卡尔曼滤波器对计算的所述下降率进行滤波;基于当前时刻的评估门限,对经滤波后的所述下降率进行评估;
所述数据输出部,输出所述数据评估部的评估结果。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111112325.5A CN113823409B (zh) | 2021-09-23 | 2021-09-23 | 一种透析中低血压事件风险的评估方法和*** |
PCT/CN2022/099129 WO2023045425A1 (zh) | 2021-09-23 | 2022-06-16 | 一种透析中低血压事件风险的评估方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111112325.5A CN113823409B (zh) | 2021-09-23 | 2021-09-23 | 一种透析中低血压事件风险的评估方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113823409A true CN113823409A (zh) | 2021-12-21 |
CN113823409B CN113823409B (zh) | 2023-12-29 |
Family
ID=78921048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111112325.5A Active CN113823409B (zh) | 2021-09-23 | 2021-09-23 | 一种透析中低血压事件风险的评估方法和*** |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113823409B (zh) |
WO (1) | WO2023045425A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023045425A1 (zh) * | 2021-09-23 | 2023-03-30 | 重庆山外山血液净化技术股份有限公司 | 一种透析中低血压事件风险的评估方法和*** |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117009831B (zh) * | 2023-10-07 | 2023-12-08 | 山东世纪阳光科技有限公司 | 一种精细化工事故风险预测评估方法 |
CN117798744B (zh) * | 2024-02-29 | 2024-05-10 | 茌平县汇通机械制造有限公司 | 一种数控机床运行状态监测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109195508A (zh) * | 2016-07-22 | 2019-01-11 | 爱德华兹生命科学公司 | 平均动脉压(map)导出的对未来低血压的预测 |
CN109195507A (zh) * | 2016-07-08 | 2019-01-11 | 爱德华兹生命科学公司 | 低血压分析参数的预测加权 |
CN110648755A (zh) * | 2019-09-10 | 2020-01-03 | 云南博亚医院有限公司 | 一种血液透析质量评估及管理*** |
CN111939353A (zh) * | 2019-05-14 | 2020-11-17 | 吴元昊 | 血液透析中低血压事件的预测模型的构建方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9002440B2 (en) * | 2010-07-08 | 2015-04-07 | Intelomed, Inc. | System and method for characterizing circulatory blood flow |
CN104703534B (zh) * | 2012-08-28 | 2017-11-03 | 弗雷泽纽斯医疗保健控股有限公司 | 通过测量相对血容量、血压和心率来检测透析中病态事件的指示的方法 |
DE102013101989A1 (de) * | 2013-02-28 | 2014-09-18 | B. Braun Avitum Ag | Fuzzy-Logik |
CN108986419B (zh) * | 2018-10-17 | 2021-02-26 | 暨南大学 | 一种用于血液透析的数据报警方法 |
CN113823409B (zh) * | 2021-09-23 | 2023-12-29 | 重庆山外山血液净化技术股份有限公司 | 一种透析中低血压事件风险的评估方法和*** |
-
2021
- 2021-09-23 CN CN202111112325.5A patent/CN113823409B/zh active Active
-
2022
- 2022-06-16 WO PCT/CN2022/099129 patent/WO2023045425A1/zh unknown
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109195507A (zh) * | 2016-07-08 | 2019-01-11 | 爱德华兹生命科学公司 | 低血压分析参数的预测加权 |
CN109195508A (zh) * | 2016-07-22 | 2019-01-11 | 爱德华兹生命科学公司 | 平均动脉压(map)导出的对未来低血压的预测 |
CN111939353A (zh) * | 2019-05-14 | 2020-11-17 | 吴元昊 | 血液透析中低血压事件的预测模型的构建方法 |
CN110648755A (zh) * | 2019-09-10 | 2020-01-03 | 云南博亚医院有限公司 | 一种血液透析质量评估及管理*** |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023045425A1 (zh) * | 2021-09-23 | 2023-03-30 | 重庆山外山血液净化技术股份有限公司 | 一种透析中低血压事件风险的评估方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN113823409B (zh) | 2023-12-29 |
WO2023045425A1 (zh) | 2023-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113823409A (zh) | 一种透析中低血压事件风险的评估方法和*** | |
CN109222949B (zh) | 心率检测方法和心率检测装置 | |
JP4322227B2 (ja) | 血液力学パラメータの測定装置 | |
CN107809948B (zh) | 用于获知血压走向的方法和设备 | |
JP5019457B2 (ja) | 急激な血圧変化の検知 | |
KR101630887B1 (ko) | 시간 종속 측정 신호를 처리하기 위한 방법 및 장치와 컴퓨터 판독 가능한 기록 매체 | |
CA2411030C (en) | Apparatus for controlling heart assist devices | |
Fetanat et al. | A physiological control system for an implantable heart pump that accommodates for interpatient and intrapatient variations | |
EP1896725B1 (en) | An apparatus for controlling blood flow in an extracorporeal circuit | |
AU2011213041B2 (en) | System and method for diagnosing sleep apnea based on results of multiple approaches to sleep apnea identification | |
WO2001097059A2 (en) | Combining measurements from different sensors | |
CN104921736B (zh) | 一种包含参数估计功能滤波模块的连续血糖监测设备 | |
AU2008281616A1 (en) | Method and apparatus for measuring breathing rate | |
ITTO950560A1 (it) | Metodo e apparecchiatura di dialisi automatica | |
KR20150113700A (ko) | 진단 시스템 및 방법 | |
AU2007255448A1 (en) | Prediction of rapid symptomatic blood pressure decrease | |
US20230122618A1 (en) | Method and system for postdialytic determination of dry weight | |
CN204261114U (zh) | 在体外血液处理期间判定内部过滤的方法和装置 | |
EP3197517A2 (en) | Method and apparatus for determining aortic valve opening | |
CN116682566B (zh) | 一种血液透析的数据处理方法及*** | |
CN111860102A (zh) | 用于在含噪环境中分析***的状态的设备和方法 | |
AU2019231059A1 (en) | Heartbeat rate calculation device and method | |
US20080167562A1 (en) | Device For Determining the Transition Between Systole and Diastole | |
WO2020220403A1 (zh) | 疲劳状态检测方法、装置、设备及存储介质 | |
CN113365553A (zh) | 用于确定患者的血压值的方法、血压测量设备和透析*** |
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 |