CN105592868B - 流体分析仪和相关联方法 - Google Patents

流体分析仪和相关联方法 Download PDF

Info

Publication number
CN105592868B
CN105592868B CN201480024298.8A CN201480024298A CN105592868B CN 105592868 B CN105592868 B CN 105592868B CN 201480024298 A CN201480024298 A CN 201480024298A CN 105592868 B CN105592868 B CN 105592868B
Authority
CN
China
Prior art keywords
mrow
msub
theta
pressure
msup
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
Application number
CN201480024298.8A
Other languages
English (en)
Other versions
CN105592868A (zh
Inventor
格雷格·皮特菲尔德
西奥多·卡尔沃斯基
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Atrium Medical Corp
Original Assignee
Atrium Medical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Atrium Medical Corp filed Critical Atrium Medical Corp
Priority claimed from PCT/US2014/029319 external-priority patent/WO2014144769A1/en
Publication of CN105592868A publication Critical patent/CN105592868A/zh
Application granted granted Critical
Publication of CN105592868B publication Critical patent/CN105592868B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/03—Detecting, measuring or recording fluid pressure within the body other than blood pressure, e.g. cerebral pressure; Measuring pressure in body tissues or organs
    • A61B5/036—Detecting, measuring or recording fluid pressure within the body other than blood pressure, e.g. cerebral pressure; Measuring pressure in body tissues or organs by means introduced into body tracts
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M1/00—Suction or pumping devices for medical purposes; Devices for carrying-off, for treatment of, or for carrying-over, body-liquids; Drainage systems
    • A61M1/71—Suction drainage systems
    • A61M1/74—Suction control
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816—Measuring devices for examining respiratory frequency
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4848—Monitoring or testing the effects of treatment, e.g. of medication
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6846—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive
    • A61B5/6847—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive mounted on an invasive device
    • A61B5/6852—Catheters
    • 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00—General characteristics of the apparatus
    • A61M2205/33—Controlling, regulating or measuring
    • A61M2205/3331—Pressure; Flow
    • A61M2205/3334—Measuring or controlling the flow rate
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00—General characteristics of the apparatus
    • A61M2205/33—Controlling, regulating or measuring
    • A61M2205/3331—Pressure; Flow
    • A61M2205/3344—Measuring or controlling pressure at the body treatment site
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00—General characteristics of the apparatus
    • A61M2205/33—Controlling, regulating or measuring
    • A61M2205/3379—Masses, volumes, levels of fluids in reservoirs, flow rates
    • A61M2205/3389—Continuous level detection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2205/00—General characteristics of the apparatus
    • A61M2205/50—General characteristics of the apparatus with microprocessors or computers
    • A61M2205/52—General characteristics of the apparatus with microprocessors or computers with memories providing a history of measured variating parameters of apparatus or patient
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2210/00—Anatomical parts of the body
    • A61M2210/10—Trunk
    • A61M2210/101—Pleural cavity

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Hematology (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Vascular Medicine (AREA)
  • Anesthesiology (AREA)
  • External Artificial Organs (AREA)

Abstract

设备、计算机可读程序和方法通过使用与流体收集罐相关联的引流管的端处的至少一个压力值和远离所述收集罐的位置处的引流管内的至少一个压力值确定患者参数，包括通过引流管从患者的胸腔引流的流体的体积和/或流速。用非线性求解器对所述压力值进行处理以确定所述患者参数。

Description

流体分析仪和相关联方法
相关申请
本申请要求享有于2013年3月15日递交的美国临时专利申请序列号61/794,712和于2014年3月14日递交的美国非临时专利申请序列号14/212,336的优先权，其全部内容通过引用并入本文。
关于联邦资助的声明
没有设想或缩减本发明的实施例以通过联邦资助或补助来实施。
技术领域
本发明的实施例涉及其中期望知道从流体源诸如个体的胸腔移除的流体的量的治疗装备和方法。
背景技术
有几种疾病状态和损伤，需要从患病个体的胸腔移除流体。临床医生监测所移除的流体的量，以便确定愈合过程的进展、疾病的进展和用于实现这样的引流的设备的功能。通常，患者具有***到胸腔中的引流管并且该管引流到称为胸腔引流装置(CDU)的特殊密闭容器中。存在用于测量密封容器中的流体的各种装置并且每个装置与成本相关联，无论是在作出观察或测量的临床医生的时间方面，还是在提供用于作出测量并且使临床医生或医院信息***注意到测量数据的传感器和信令设备的装备成本方面。
例如，但不限于，一些设备的特征在于用于测量容器的重量或改变的重量的一个或多个应变计。随着应变计的精度增加，应变计变得更昂贵并且更容易受到损坏。应变计特别容易被损坏，并且如果设计得较可靠稳定则对重量改变具有较小的分辨率。
记录患者的流体引流体积水平所需的、临床人员对CDU上的患者的监测是棘手、费时并且因此昂贵的，并且不像原本可能需要的那样频繁地测量或检查。引流设备通常被保持在地板上，并且要求临床医生四肢着地趴下来读取流体液位，并且例如在使用期间CDU不小心被碰倒则还要记得调整误差。
发明内容
本发明的实施例旨在提供用于在具有引流管的患者中确定至少一个患者参数的方法、计算机程序和设备，所述引流管具有第一端和第二端，所述第一端装配到患者胸腔中，并且第二端接纳在收集容器中。如本文所用，术语“患者参数”是指与引流管相关联的患者参数，以举例而不是限制的方式，包括：从胸腔移除的流体(气体和/或液体)的体积、移除流体的流率、保持所移除的流体的液体部分的容器的本体积、过滤设备和管的本功能以及患者呼吸状态。
本发明的方法、计算机程序和设备的实施例的特征在于，所述收集容器和引流管中的至少一者具有：第一压力传感器，所述第一压力传感器感测(i)接近所述收集容器的位置中的所述引流管和(ii)所述收集容器中的至少一者内的压力；以及第二压力传感器，所述第二压力传感器感测远离所述收集容器的位置处的所述引流管内的压力。也就是说，所述第一压力传感器与所述收集容器紧密地相关联，不管是位于所述收集容器内还是靠近所述收集容器的引流管中。所述第二压力传感器位于远离所述收集容器的引流管中，更接近患者。
所述方法包括从所述第一压力传感器获得至少第一压力值和从所述第二压力传感器获得至少第一压力值的步骤。所述方法进一步包括通过非线性求解器由所述压力值计算至少一个患者参数的步骤。
本发明方法的一个实施例的特征在于，使每个压力值与时间相关联并且随时间获得多个连续压力值。所述连续压力值优选在时间上以定时间隔分开；也就是说，它是周期性的。
本发明方法的实施例非常适合由计算机装置执行。如本文所用，术语“计算机装置”是指在设备内部或与设备分开和不同地承载的计算机或中央处理单元(CPU)，诸如在大型计算机上承载的CPU、服务器、传统的台式计算机和膝上型计算机、微控制器、车载计算机诸如树莓派(RaspberryPi)以及手持式通信设备诸如手机和平板设备、通过互联网和/或无线通信设备链接的场外计算机等。优选地，将第一压力值和连续压力值存储在计算机存储器中。如本文所用，术语“计算机存储器”是指数据存储，例如与CPU信号通信并且可以由这样的CPU处理的芯片、盘、存储器驱动器、硬盘驱动器、闪存驱动器等。例如，CPU从存储器存取这种第一时间压力值和连续压力值并且执行对这种数据进行组织、计算至少一个患者参数并且产生显示所述一个患者参数的输出的步骤。
如本文所用，术语“非线性求解器”是指假设/模型数学计算机程序意义上的机器学习过程。用已知的输入和输出对非线性求解器进行训练，并且导出或编程经训练的非线性求解器，以接收第一压力值和第二压力值以及所有其他连续压力值，通过假设/模型对这样的值进行处理，以获得患者参数。
实施例的特征在于以从由支持向量机方案、人工神经网络、遗传算法和遗传编程组成的组中选择的一种或多种方案来实施的非线性求解器。详细的讨论将以线性回归假设和机器学习过程的形式描述非线性方程。
例如，一个非线性方程是线性回归假设：
如上所用：
T＝表示当前时间(t＝0)，或在数据之后的固定时间的点
t＝整数时间样本，0最新，N最老
f(x,y)＝取决于x和y的未知函数
h(T)＝流体罐中的流体的体积
PV(t)＝出口处的作为时间的函数的压力，
PC(t)＝患者的CDU侧的作为时间的函数的压力，
N＝计算所需的时间样本的数
Theta＝下面列出的所有Theta的集
Thetabias＝针对假设的任意DC偏差常数
Thetac,t＝针对在时间t的CDU压力的任意参数系数
ThetaV,t＝针对在时间t的出口压力的任意参数系数。
线性回归方程具有用“m”个值的训练集表示的成本函数J，表示如下：
用“m”个项的训练集来确定未知的“Theta”值，如下：
例如，但不限于，通过从由梯度下降、粒子群优化和人工蜂群组成的组中选择的一种或多种优化算法来求解非线性方程及其导数。
所述一个方法的特征在于如下推导的Theta：
J(Theta)的偏导数表示为：
其中Xx是与Thetax相关联的偏差(Bias)、PV或PC参数。
将体现本发明的特征的另一个非线性方程表示为：
其中：
T＝表示当前时间(t＝0)，或在数据之后的固定时间点
t＝整数时间样本，0最新，N最老
f(x,y)＝取决于x和y的未知函数
h(T)＝流体罐中的流体的体积
Pv(t)＝出口处的作为时间的函数的压力，
PC(t)＝患者的CDU侧的作为时间的函数的压力，
N＝计算所需时间样本数
Theta＝以下列出的所有Theta的集
Thetabias＝针对假设的任意DC偏差常数
Thetac,t＝针对在时间t的CDU压力的任意参数系数
ThetaV,t＝针对在时间t的出口压力的任意参数系数。
本发明的进一步的实施例旨在提供用于确定至少一个患者参数的计算机可读程序。所述计算机可读程序包括用于使具有存储器的计算机处理单元从第一压力传感器获得至少一个第一压力值并且从第二压力传感器获得至少一个第一压力值的指令。所述计算机可读程序进一步指令计算机装置通过非线性求解器由所述压力值来计算至少一个患者参数。
计算机可读程序在计算机装置上操作或运行，以接收所述至少一个第一压力值和例如采用定时间隔的随时间的任何连续压力值，并且将这样的数据存储在存储器中。组织和使用所述数据以计算用于显示的一个或多个患者参数。根据上述非线性方程和方法对数据进行处理。
本发明的进一步的实施例涉及一种用于确定至少一个患者参数的设备。所述设备包括：第一压力传感器，所述第一压力传感器被构造和布置为感测(i)接近收集容器的位置中的引流管或(ii)收集容器本身中的至少一者内的压力；以及第二压力传感器，所述第二压力传感器被构造和布置为感测远离收集容器的位置处的引流管内的压力。所述第一压力传感器和所述第二压力传感器被构造和布置用于放置为与具有存储器的计算装置信号通信，以从所述第一压力传感器获得第一压力值并且从所述第二压力传感器获得第一压力值。所述计算机装置通过非线性求解器由所述压力值来计算至少一个患者参数。
优选的设备包括计算机装置。所述计算机装置优选接收多个第一压力值和例如但不限于采用定时间隔的随时间的任何连续压力值，所述多个第一压力值和例如但不限于采用定时间隔的随时间的任何连续压力值被存储在存储器中、被组织和处理以计算一个患者参数。
组织和使用所述数据以计算用于显示的一个或多个患者参数。根据上述非线性方程和方法对数据进行处理并且在下面的详细描述中予以更全面地描述。
本发明的方法、计算机可读程序和设备允许监测从患者的胸腔移除的流体。这种流体包括任何生物流体，诸如气体、水、血清、血液、血浆、淋巴液、脓、治疗性输注流体及其混合物。术语“气体”用于指空气和施予以用于治疗目的的气体，诸如氧气、氮气、氦气和一氧化氮。本发明方法、计算机可读程序和设备的实施例允许计算气体体积和液体体积以及流体的相应流率。
在查看下面简单描述的附图和研究下面的本发明的详细描述后，这些和其他特征和优点对于本领域技术人员将是显而易见的。
附图说明
图1示出了根据示例实施例的用于确定从患者的胸腔引流的流体的体积和/或流量的一个示例性流体分析仪。
图2是根据示例实施例以进一步的示例性细节示出图1的电子控制模块的示意图。
图3是根据示例实施例图示用于收集和存储时间序列压力值的一个示例性方法的流程图。
图4是根据示例实施例图示用于处理图3的所存储的时间序列压力值以测量胸腔引流流体的一个示例性方法的流程图。
图5是根据示例实施例图示将图1的引流设备模型化的一个示例性电子电路的示意图。
图6描绘了根据示例实施例的机器学习过程。
具体实施方式
将相对于用于测量来自外部源诸如患者的胸腔的流体体积的设备、***和方法对本发明的实施例进行详细描述。现在转向图1，描绘了体现本发明的特征的设备，在本文中有时称为流体分析仪，并且总的用标号100表示。
图1示出了用于确定从源诸如患者110的胸腔引流的流体的体积和/或流量的一个示例性流体分析仪100。在图1的示例中，在引流设备120内实施流体分析器100，所述引流设备120操作例如以使用导管或管道诸如双腔导管109收集从患者110的胸腔引流的流体。在图1的示例中，引流设备120被示为耦接到真空源112，以促进从患者110的流体引流。
引流设备120被示为包括电子控制模块(ECM)102，所述电子控制模块(ECM)102操作以控制对流体通过导管109到流体罐108中的引流。如所描绘，ECM 102包括用户界面106，用于与设备120的用户(例如，临床人员)交互。流体罐108是例如收集从患者的胸腔引流的流体的胸腔引流单元(CDU)。ECM 102从(a)感测导管109的出口通路的患者连接附近的压力的出口压力传感器126和(b)感测源自相同患者点的引流管和收集***的端处的压力的CDU压力传感器128接收压力值。出口压力传感器126可以位于导管109的排放管上，尤其在导管109不是双腔型的情况下。与引流设备120和/或导管109一起布置的额外压力传感器可以与***100一起使用，而不脱离本发明的范围。
引流设备120可以包括其他组件(例如，阀、过滤器、传感器、转换器、接口等)，而不脱离本发明的范围。引流设备120处理来自压力传感器126和128的时间序列的压力值，以确定从患者110引流的流体的流体体积和/或流率。ECM 102可以确定其他度量，诸如，但不限于：CDU过滤器堵塞、患者呼吸、患者吸入、患者呼出、每分钟的呼吸、患者呼吸困难(例如，咳嗽、呛咳、笑、哭等)、管堵塞和闭塞或堵塞的程度、线路中的流体以及流体的患者倾倒。
图2是以进一步的示例性细节示出ECM 102的示意图。最好结合下面的描述查看图1和图2。ECM 102被示为包括用户界面106、处理器206、存储器208和传感器接口210。存储器208被说明性地示为存储软件104，所述软件104包括机器可读指令，所述机器可读指令当由处理器206执行时实施用于确定从患者110引流的流体的体积和/或流率的引流设备120内的功能。传感器接口210与出口压力传感器126和CDU压力传感器128两者耦接并且包括模拟-数字转换器，以将感测的压力值提供到软件104。软件104包括数据采集器212和计算器214。数据采集器212操作以随时间从压力传感器126和128收集压力值并且将这些值作为时间序列压力数据220存储在存储器208内。时间序列压力数据220包含经由传感器接口210从压力传感器126和128周期性地读取的压力值。计算器214操作以处理时间序列压力数据220，以确定罐108内的流体的体积并且可选地确定流体进入罐108的流率。计算器214可以包括学习算法，所述学习算法基于从传感器接口210读取的压力值和定义的体积自动学习这些压力值和流体体积和/或流率之间的关系。计算器214基于每个压力传感器126和128处的识别的压力变化生成压力特征数据222。压力特征数据222可在下文中称为出口压力特征和/或CDU压力特征。计算器214还生成流体测量结果224以包括罐108内的流体体积和进入罐108的流体流率中的一者或两者。数据采集器212和计算器214可以同时操作。
图3是图示用于读出和存储时间序列压力值的一个示例性方法300的流程图。图4是图示用于处理所存储的时间序列压力数据以测量胸腔引流流体的一个示例性方法的流程图。例如，方法300和400是例如在ECM102的软件104内实施的在计算机上运行的单独处理线程。最好结合下面的描述查看图1到图4。
在步骤302中，方法300从出口压力传感器读取第一压力值。在步骤302的一个示例中，由处理器206执行的软件104使用传感器接口210从出口压力传感器126读取第一压力值。在步骤304中，方法300从CDU压力传感器读取第二压力值。在步骤304的一个示例中，由处理器206执行的软件104使用传感器接口210从CDU压力传感器128读取第二压力值。
在步骤306中，方法300存储步骤302和304的压力值作为时间序列压力数据。在步骤306的一个示例中，数据采集器212将步骤302和304的压力值存储在存储器208内作为时间序列压力数据220。
方法300周期性地重复以收集和存储压力值并且基于第一压力值和第二压力值计算流体参数。
在步骤402中，方法400从存储器读取时间序列压力数据。在步骤402的一个示例中，计算器214从存储器208读取时间序列压力数据220。在步骤404中，方法400从时间序列压力数据提取特征。在步骤404的一个示例中，计算器214从时间序列压力数据220提取压力特征并且将这些压力特征作为压力特征数据222存储在存储器208内。
步骤406是判定。在步骤406中，方法400确定在步骤404中确定的特征是否有效，并且如果这些特征有效，则方法400以步骤408继续；否则，方法400以步骤402继续。
在步骤408中，方法400通过非线性求解器对特征进行处理。在步骤408的一个示例中，计算器214通过非线性方程处理压力特征数据222，以生成包括流体体积和流体流率中的一者或两者的流体测量结果224。在步骤408的另一个示例中，计算器214通过预先训练的支持向量机对压力特征数据222进行处理。在步骤408的另一个示例中，计算器214使用人工神经网络、遗传算法和遗传编程中的一者或多者来对压力特征数据222进行处理。例如，计算器214可以基于所收集的训练、验证和测试数据集中的一者或多者来利用数据。
步骤410是判定。在步骤410中，方法400确定在步骤408中生成的数据是否有效，并且如果该数据有效，则方法400以步骤412继续；否则，方法400以步骤402继续。
在步骤412中，方法400报告所确定的体积和其他值。在步骤412的一个示例中，计算器214将流体测量结果224输出到用户界面106。
方法400周期性地重复以基于时间序列压力数据220确定和输出流体测量结果224。
图5示出了表示引流设备120的不同部分的电阻、电感和电容的一个电子模拟500。因此，根据示例实施例，电子模拟500将引流设备120模型化。然而，可以替代地以其他方式将引流设备120模型化，而不脱离本发明的范围。电子模拟500对两个分支是固定的，其中唯一的例外是流体罐108，所述流体罐108由可变电容器模型化。给定由患者创建的动态压力信号和由每个压力传感器126、128“看见”的“电网”差异，该定时将取决于流体罐108内的流体体积而非线性地变化。流体罐108中的流体体积的增加将减小影响压力定时的电路臂中的空气电容。***内过滤器(图5的Rf)可以在电阻上变化(随时间增加)。该电阻变化应该是线性影响并且随时间逐渐发生。
图6描绘***或设备100的机器学习过程。将步骤600、610、620和640表征为训练。步骤600收集用于训练的输入和输出。将在受控设定中通过实验获得这些输入和输出，以产生训练步骤610中的输入和输出的矩阵。训练步骤610产生一个或多个输出，所述一个或多个输出被评估接受性或适应性。遗传算法和程序基于期望的输出并且基于被编译成一个值的训练期间的所有预期输出在适应性的基础上对输入和输出进行评估，其中在步骤640对模型进行可选的改变。如果模型产生可接受的结果，则在步骤630将模型导出到设备100。
在步骤650将模型集成到设备100中。而且，在步骤660使用模型计算针对模型的输入参数。在步骤670使用模型对输入参数进行处理，并且在步骤680，设备100根据输出结果行动。
示例性方案：
在某些实施例中，使用非线性方程假设来计算流体罐108中的流体的体积。可以使用固定模型或动态和自校正模型。通常期望固定模型，使得常数(重量或质量)不改变并且可以在软件104内原样使用，而不用再训练。假设对于固定模型或动态和自校正模型将相同，但可能需要针对模型对方程进行优化。
假设：
令：
·T＝表示当前时间(t＝0)，或在数据之后的固定时间点
·t＝整数时间样本，0最新，N最老
·f(x,y)＝取决于x和y的未知函数
·J(Theta)＝针对Theta的假设vs.实际的成本
·V(T)＝流体罐108中的流体的体积(实际)
·h(T)＝流体罐108中的流体的体积(假设)
·PV(t)＝例如从压力传感器126确定的出口处的作为时间的函数的压力(或出口压力特征，以及其他参数)
·PC(t)＝例如从压力传感器128确定的患者的CDU侧的作为时间的函数的压力(或CDU压力特征，以及其他参数)
·m＝训练集中的样本数
·N＝计算所需时间样本数
·Theta＝下面列出所有Theta的集
·Thetabias＝针对假设的任意DC偏差常数8
·Thetac,t＝针对在时间t的CDU压力的任意参数系数
·ThetaV,t＝针对在时间t的出口压力的任意参数系数
为了近似V(T)，设置以下线性回归假设：
可以用“m”个项的训练集如所示表示成本函数J，如下：
可以用“m”个项的训练集找到未知的“Theta”值，如下：
有多种方式可以求解或修改该方程。求解或修改方程的一种可能的方案是通过使用“梯度下降”方案/算法。梯度下降方案提供对方程解的良好近似。替代方法包括最近和更快的优化***，诸如粒子群优化***。使用人工蜂群(ABC)能够以低开销快速求解方程。另一种方案是使用外部数学包诸如MatLab、R或Octave来求解“Theta”。这些求解方案是现有技术中已知的。例如，诸如以软件104的形式将Theta的这些值输入或以其他方式并入ECM102中。这种方案对于任何静态“Theta”组通常是足够的。
梯度下降需要关于“Theta”的成本函数的偏导数向下移动“Theta”，如下：
J(Theta)的偏导数可以表示为：
其中Xx是与Thetax相关联的偏差、PV或PC参数。
如果上述方程假设不收敛，可以添加二次幂项(second power term)和可能甚至三次幂项(third power term)，以帮助实现收敛。如果解展现高偏差情况(欠拟合)，则将需要并入额外的幂项。如果存在太多的参数，这称为高方差(过拟合)，则可能会出现相反的问题。可以通过一些技术来解决过拟合，诸如正规化和添加更多的样本。
下面是包括较高次项的第二方程假设，例如如果一阶模型没有令人满意地收敛，则使用该第二方程假设：
通过使用推导参数，可以改善上述数学模型的精度。例如，可以从压力波形推导两个波形之间的最大值和时间差。使用推导的参数数据集而不是大的时间序列上的分析将减小成本计算中的“误差”项。在不对原始数据进行参数提取的情况下，将需要处理更多的样本，以消除误差。
对方程的更复杂的求解是通过将时间信号分解成它的固有模态函数或IMF。固有模态函数(IMF)可以用于以增加的复杂性/计算为代价提高精度。最近的工作由NordenE.Huang以希尔伯特-黄变换(HHT)的形式完成。该方案使用经验模态分解(EMD)和希尔伯特谱分析(HAS)来提高精度，缺点是计算费用增加。在特定实施例中，测试上面讨论的数学方程假设，并且仅在数学假设不令人满意的情况下实施HHT技术。
引流设备120能够至少部分地使用来自压力传感器126、128的压力数据确定除流体罐108中的流体的体积以外的信息。另外，某些实施例进一步包括一个或多个阀传感器(未示出)或阀的编程定时以检测引流设备120中的一个或多个阀的状态。在这些实施例中，ECM 102使用来自阀传感器的阀数据连同来自压力传感器126、128的数据例如以确定额外的信息。ECM 102例如使用线性回归技术，诸如类似于上面讨论的用于确定流体罐108流体体积的技术确定该额外信息。作为另一示例，在一些实施例中，ECM 102使用逻辑回归(也称为分类、概率输出/类呈现)确定额外的信息。
如上面所讨论，引流设备120能够估计流体罐108中的流体的近似体积。由于流体罐108中的***的“空气电容(air capacity)”的改变，这样的估计的精度与流体体积成反比。另外，某些实施例能够使用表征随时间的过滤器改变的信息来检测细微的体积改变。引流设备120例如通过检测由于过滤器的润湿导致的可变电阻器的电阻的改变来确定过滤器改变。此外，一些实施例能够基于流体罐108体积可如何随时间改变来检测诸如由于断开的缝线或缝合导致的重倾倒(heavily dumping)患者110。
一些实施例能够经由在引流设备120内包括来自患者泄压阀130的操作的阀数据和来自压力传感器128的CDU压力数据的线性回归模型来确定通过流体罐108的空气流(也被认为是流体)。例如通过使用来自***的反馈和表征当前空气通过的改变的信息，精度增加。此外，一些实施例能够允许空气的潮流进入流体罐108并且保持阀以近似释放前的空气压力，以进一步提升精度。在这些实施例中，引流设备120通常仅保持压力一瞬间以实现更好的临床数据，而不会显著影响治疗。将围绕特征诸如这一个特征设置安全限制。压力时间衰变数据可以用于帮助近似用以在一段时间内将被相加的通过的体积/团块，以报告空气流率。
可以使用类似的方案找到呼吸速率，但可能需要大的时间数据集(或来自上述数据集的样本)。
引流设备120的某些实施例能够使用一个或多个逻辑回归模型确定以下事件中的一者或多者的存在或存在的概率：
·来自线性回归的体积数据有效
°压力信号与训练示例不足够相似
°报告体积中的可能的误差
·流体罐108过滤器堵塞，请更换CDU
·患者正在呼吸
·患者正在吸入
·患者正在呼出
°基于上述事件计算的每分钟的呼吸次数
·患者情况，诸如患者呼吸困难(例如咳嗽、呛咳、笑、哭等)
·管堵塞(例如，完全闭塞)
·管闭塞程度
·线路中的流体
·患者倾倒流体
本文中所公开的***和方法的某些实施例实现以下优点中的一个或多个：
·更好的患者恢复时间；
·可以将管线清除设备(TLC)嵌入作为引流设备120的部分，以在需要时清除引流管，并且检测依赖循环、监测流体收集和速率、空气流和速率；
·当以高置信度存在依赖循环和流体收集时，可以激活TLC。这将延长电池寿命，不需要的TLC起动(firing)对治疗的干扰较小。在没有堵塞的情况下激活TLC时，患者风险较小。
·通过更好地监测改善患者安全，例如，***可以检测患者是否正确连接到***和/或患者是否正在呼吸。
在这种极端的情况下可以召唤临床帮助。
·***可以确定引流管是否闭塞或具有阻塞路径的流体，对流体阻塞起动TLC，并且警告临床管扭结等。
·***在数据中提供更多的值以记录在患者记录中或用于研究或者发送到中心护士站
·联网到中心护士站的***可以向护理人员提供关于以下的信息：
°呼吸率(患者管断开、睡眠或非睡眠相关呼吸暂停)
°患者痛苦
°流体的体积
°空气流率
°管通畅
°患者积极倾倒流体。
在某些实施例中，在捕获数据的训练集和然后从ECM记录的测试集之后实施本发明的方法、设备和软件。数据的训练集包含例如以10ml的体积增量的许多样本。通过额外的训练样本以及对成本函数和/或参数集的修改来解决计算中的高方差(内插的样本可具有比预期误差更高的误差)。
在各种示例性实施例中，样本数据集可以包含在不同的时间/体积的“坏”数据。“坏”数据将被视为其中***模型已经被违反并且不再以足够的预测确定性发挥作用的测量。一个“坏”数据示例是当流体被截留在引流线路或安装在线路上的夹具中时获取的数据。其他“坏”示例可以包括呼吸的突然改变，诸如咳嗽和没有呼吸。另一个示例可以包括以不同的体积的空气泄漏。
本领域技术人员将认识到根据需要将被用于分析数据以找到“theta”并且然后求解出体积和其他参数的各种工具诸如Matlab和称为Octave的免费工具的可用性。
逻辑回归
在某些实施例中，应用逻辑回归以开发生成输入的分类(例如，是/否为闭塞状态、数据有效/数据不是有效的、患者呼吸等)的***。逻辑回归非常类似于上述线性回归，其中方程假设由其输出是从0到1的Sigmoid函数包裹。
被表示为逻辑假设，方程可以采取以下形式：
改变成本函数，以保持凸解，以帮助避免局部最小值。成本函数的新形式，其中如果y属于该类则为1并且如果y不属于该类则为0：
偏导数仅在分类器y(0或1)方面类似：
例如，在一个vs.所有结构中完成分类。也就是说，对于每一个分类，将对神经网络进行训练，使得仅被训练的类是“在类中”或1。所有其他样本将是“在类外”或以其他方式设为0。这意味着将存在针对所期望的每一个分类设置的不同的“theta”。
支持向量机方案
支持向量机对于分类是非常适合和专业的。通过将所有上述逻辑回归分成经训练的支持向量机(SVM)，得到生成存在于输入数据中的类的概率的计算上高效的***。
SVM与逻辑回归的不同之处在于SVM在其超平面(n维空间)内“存储”正例。然后，其报告样本与所学***面中的已知正相关的样本的概率。根据各种示例实施例，被称为“libSVM”的开源和优化的C++库可以被用来设计原型和/或集成到引流设备120中，用于类预测。
为了让机器学习***泛化(generalize)到没有教导的情况，特征本身需要是一般化的。例如，使用压力传感器126、128两者的绝对压力量值可能无法良好泛化。施加到患者110的压力和由患者施加的压力将与(i)规定的负压和由患者的呼吸***施加的力，以及(ii)患者的胸膜间隙中的空气死体积的量相关。
从压力传感器126、128接收到的数据取决于流体罐108流体体积，并且可以形成更一般的特征。以下是可以在某些实施例中使用的特征的部分列表：
·峰值延迟(例如，以微秒或毫秒为单位的时间)-使用最慢的压力传感器波形(最有可能的流体罐108压力)，报告从其峰值压力回到出口压力峰值的峰值压力的时间。
·出口的周期最大-最小(压力)
·出口压力峰值之间的距离(以时间为单位，诸如毫秒)
·流体罐108的周期最大-最小(压力)
·流体罐108压力峰值之间的距离(以时间为单位，诸如毫秒)
·在出口压力峰值，流体罐108和出口之间的压力差
·在CDU压力峰值，流体罐108和出口之间的压力差
·出口压力的斜率(导致峰值)
·流体罐108压力的斜率(导致峰值)
·压力曲线下的面积
·阀门 定时(开/关时间)
本发明的一个实施例的特征在于通过遗传编程(GP)的机器学习。不同于优化假设的解决方案的遗传算法和机器学习方案，GP创建“程序”来解决问题。学习过程的结果是解决应用问题的“公式/源代码”。
有几种不同的方案来遗传编程，包括以下：
-GP—基于指令的功能树的标准遗传编程
-LGP—线性遗传编程，线性处理的固定指令长度(快速)
-CGP—笛卡尔遗传编程，随机连线的随机指令点的固定网格
可以生成C代码以容易地嵌入软件104内的市售LGP应用程序的示例是Discipulus5遗传编程预测建模(RML Technologies公司，科罗拉多州利特尔顿，http://www.rmltech.com/)。
在某些实施例中，可以省略出口压力传感器126。这可以在例如外部施加的刺激的特性是已知的情况下发生。仍然可以基于已知的患者源使用流体罐108中的CDU压力传感器128来确定体积和其他特征。例如，在患者在心肺机或呼吸机上的情况下，这样的外部刺激是已知的。
可以在不脱离本发明的范围的情况下在上述方法和***中做出改变。因此应当指出的是，包含在上述说明中或在附图中示出的事项应当被解释为说明性的而不是限制性的。所附权利要求旨在覆盖本文中所描述的所有一般的和具体的特征以及因为语言的关系可以说落在它们之间的本发明方法和***的范围的所有陈述。

Claims (30)

1.一种用于在具有引流管的患者中确定至少一个患者参数的方法，所述引流管具有第一端和第二端，所述第一端装配到所述患者的胸腔中，并且所述第二端与收集容器流体连通，所述收集容器和所述引流管中的至少一者具有：第一压力传感器，所述第一压力传感器感测(i)接近所述收集容器的位置中的所述引流管和(ii)所述收集容器中的至少一者内的压力；以及第二压力传感器，所述第二压力传感器感测远离所述收集容器的位置处的所述引流管内的压力，所述方法包括以下步骤：
a.分别从所述第一压力传感器和所述第二压力传感器获得连续的压力值；
b.把所述压力值作为时间序列压力数据存储在计算装置的存储器中；
c.通过使用非线性方程的非线性求解器由所述时间序列压力数据利用所述计算装置计算至少一个患者参数；
d.响应于所述患者参数，激活管线清除设备，以便清除所述引流管。
2.根据权利要求1所述的方法，其中，在定时间隔处获得所述压力值。
3.根据权利要求2所述的方法，其中，所述压力值与从由阀状态、时间、患者吸入状态和管状态组成的组中选择的一个或多个事件、状态或值相关联。
4.根据权利要求1所述的方法，其中，所述计算装置产生显示所述至少一个患者参数的输出。
5.根据权利要求1所述的方法，其中，利用从由支持向量机方案、人工神经网络、遗传算法和遗传编程组成的组中选择的一种或多种方案来实施所述非线性方程。
6.根据权利要求1所述的方法，其中，所述非线性方程是线性回归假设：
<mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>,</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>Theta</mi> <mrow> <mi>b</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>N</mi> </msubsup> <mo>&amp;lsqb;</mo> <msub> <mi>Theta</mi> <mrow> <mi>c</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>v</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
其中，
N＝计算所需的时间样本的数，
T＝表示当前时间(t＝0)，或在数据之后的固定时间的点，
t＝整数时间样本，0最新，N最老，
j是时间样本标号，取值0至N，
f(x,y)＝取决于x和y的未知函数，
h(T)＝流体罐中的流体的假设体积，
Pv(t)＝出口处的作为时间t的函数的压力，
Pc(t)＝患者的胸腔引流装置CDU侧的作为时间t的函数的压力，
Theta＝下面列出的所有Theta的集：
Thetabias＝针对假设的任意偏差常数，
Thetac,j＝针对时间样本j的胸腔引流装置CDU压力的任意参数系数，
Thetav,j＝针对时间样本j的出口压力的任意参数系数。
7.根据权利要求6所述的方法，其中，所述线性回归假设具有成本函数J，表示如下：
<mrow> <mi>J</mi> <mrow> <mo>(</mo> <mi>T</mi> <mi>h</mi> <mi>e</mi> <mi>t</mi> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mo>*</mo> <mi>m</mi> </mrow> </mfrac> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </msubsup> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>h</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>-</mo> <mi>V</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> </mrow>
其中，
m＝训练集中的样本数，
i是取值1至m的整数，
V(T)＝流体罐中的流体的实际体积。
8.根据权利要求7所述的方法，其中，用m个项的训练集来确定未知的Theta值，如下：
<mrow> <msub> <mi>min</mi> <mrow> <mi>T</mi> <mi>h</mi> <mi>e</mi> <mi>t</mi> <mi>a</mi> </mrow> </msub> <mi>J</mi> <mrow> <mo>(</mo> <mi>T</mi> <mi>h</mi> <mi>e</mi> <mi>t</mi> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mo>*</mo> <mi>m</mi> </mrow> </mfrac> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </msubsup> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>h</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>-</mo> <mi>V</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>.</mo> </mrow>
9.根据权利要求8所述的方法，其中，通过从由梯度下降、粒子群优化和人工蜂群组成的组中选择的一种或多种优化算法来求解所述非线性方程及其导数。
10.根据权利要求1所述的方法，其中，将所述非线性方程表示为：
<mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>,</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>Theta</mi> <mrow> <mi>b</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>N</mi> </msubsup> <mo>&amp;lsqb;</mo> <msub> <mi>Theta</mi> <mrow> <mi>c</mi> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>C</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>c</mi> <mn>2</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>v</mi> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>v</mi> <mn>2</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
其中，
N＝计算所需的时间样本的数，
T＝表示当前时间(t＝0)，或在数据之后的固定时间的点，
t＝整数时间样本，0最新，N最老，
j是时间样本标号，取值0至N，
f(x,y)＝取决于x和y的未知函数，
h(T)＝流体罐中的流体的体积，
Pv(t)＝出口处的作为时间t的函数的压力，
Pc(t)＝患者的胸腔引流装置CDU侧的作为时间t的函数的压力，
Theta＝下面列出的所有Theta的集：
Thetabias＝针对假设的任意偏差常数
Thetac1,j，Thetac2,j＝针对时间样本j的胸腔引流装置CDU压力的任意参数系数
Thetav1,j，Thetav2,j＝针对时间样本j的出口压力的任意参数系数。
11.根据权利要求1所述的方法，其中，所述流体是包括由气体、水、血清、血液、血浆、淋巴液、脓、治疗性输注流体及其混合物组成的组中的任一者的生物流体。
12.根据权利要求1所述的方法，其中，所述患者参数是被收集在收集容器中的液体流体的估计。
13.根据权利要求1-12中任一项所述的方法，其中，对所述非线性求解器进行训练。
14.根据权利要求1所述的方法，其中，所述患者参数指的是所述收集容器中的流体的体积。
15.一种用于在具有引流管的患者中确定至少一个患者参数的设备，所述引流管具有第一端和第二端，所述第一端装配到所述患者的胸腔中和/或建立与所述患者的胸腔的流体连通，并且所述第二端接纳在收集容器中，所述设备包括：
a.第一压力传感器，所述第一压力传感器被构造和布置为感测(i)接近所述收集容器的位置中的所述引流管和(ii)所述收集容器中的至少一者内的压力；
b.第二压力传感器，所述第二压力传感器被构造和布置为感测远离所述收集容器的位置处的所述引流管内的压力；
c.所述第一压力传感器和所述第二压力传感器被构造和布置用于放置为与具有存储器的计算装置信号通信，以使得：
i.分别从所述第一压力传感器和所述第二压力传感器获得连续的压力值；
ii.将所述压力值作为时间序列压力数据存储在所述存储器中；
iii.通过使用非线性方程的非线性求解器由所述时间序列压力数据利用所述计算装置计算至少一个患者参数；
iv.响应于所述患者参数，激活管线清除设备，以便清除所述引流管。
16.根据权利要求15所述的设备，其中，以定时间隔得出所述压力值。
17.根据权利要求15所述的设备，其中，所述计算装置产生显示所述至少一个患者参数的输出。
18.根据权利要求15所述的设备，其中，利用从由支持向量机方案、人工神经网络、遗传算法和遗传编程组成的组中选择的一种或多种方案来实施所述非线性方程。
19.根据权利要求15所述的设备，其中，所述非线性方程是线性回归假设：
<mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>,</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>Theta</mi> <mrow> <mi>b</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>N</mi> </msubsup> <mo>&amp;lsqb;</mo> <msub> <mi>Theta</mi> <mrow> <mi>c</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>v</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
其中，
N＝计算所需的时间样本的数，
T＝表示当前时间(t＝0)，或在数据之后的固定时间的点，
t＝整数时间样本，0最新，N最老，
j是时间样本标号，取值0至N，
f(x,y)＝取决于x和y的未知函数，
h(T)＝流体罐中的流体的假设体积，
Pv(t)＝出口处的作为时间t的函数的压力，
Pc(t)＝患者的胸腔引流装置CDU侧的作为时间t的函数的压力，
Theta＝下面列出的所有Theta的集：
Thetabias＝针对假设的任意偏差常数，
Thetac,j＝针对时间样本j的胸腔引流装置CDU压力的任意参数系数，
Thetav,j＝针对时间样本j的出口压力的任意参数系数。
20.根据权利要求19所述的设备，其中，所述线性回归假设具有成本函数J，表示如下：
<mrow> <mi>J</mi> <mrow> <mo>(</mo> <mi>T</mi> <mi>h</mi> <mi>e</mi> <mi>t</mi> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mo>*</mo> <mi>m</mi> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>h</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>-</mo> <mi>V</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> </mrow>
其中，
m＝训练集中的样本数，
i是取值1至m的整数，
V(T)＝流体罐中的流体的实际体积。
21.根据权利要求20所述的设备，其中，用m个项的训练集来确定未知的Theta值，如下：
<mrow> <msub> <mi>min</mi> <mrow> <mi>T</mi> <mi>h</mi> <mi>e</mi> <mi>t</mi> <mi>a</mi> </mrow> </msub> <mi>J</mi> <mrow> <mo>(</mo> <mi>T</mi> <mi>h</mi> <mi>e</mi> <mi>t</mi> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mo>*</mo> <mi>m</mi> </mrow> </mfrac> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </msubsup> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>h</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>-</mo> <mi>V</mi> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mi>i</mi> </msup> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> <mo>.</mo> </mrow>
22.根据权利要求19所述的设备，其中，通过梯度下降、粒子群优化和人工蜂群方案来求解所述非线性方程及其导数。
23.根据权利要求15所述的设备，其中，将所述非线性方程表示为：
<mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>,</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>Theta</mi> <mrow> <mi>b</mi> <mi>i</mi> <mi>a</mi> <mi>s</mi> </mrow> </msub> <mo>+</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>N</mi> </msubsup> <mo>&amp;lsqb;</mo> <msub> <mi>Theta</mi> <mrow> <mi>c</mi> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>C</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>c</mi> <mn>2</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>c</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>v</mi> <mn>1</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Theta</mi> <mrow> <mi>v</mi> <mn>2</mn> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>*</mo> <msub> <mi>P</mi> <mi>v</mi> </msub> <msup> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
其中，
N＝计算所需的时间样本的数，
T＝表示当前时间(t＝0)，或在数据之后的固定时间的点，
t＝整数时间样本，0最新，N最老，
j是时间样本标号，取值0至N，
f(x,y)＝取决于x和y的未知函数，
h(T)＝流体罐中的流体的体积，
Pv(t)＝出口处的作为时间t的函数的压力，
Pc(t)＝患者的胸腔引流装置CDU侧的作为时间t的函数的压力，
Theta＝下面列出的所有Theta的集：
Thetabias＝针对假设的任意偏差常数，
Thetac1,j，Thetac2,t＝针对时间样本j的CDU压力的任意参数系数，
Thetav1,j，Thetav2,t＝针对时间样本j的出口压力的任意参数系数。
24.根据权利要求15所述的设备，其中，所述流体是包括由气体、水、血清、血液、血浆、淋巴液、脓、治疗性输注流体及其混合物组成的组中的任一者的生物流体。
25.根据权利要求15所述的设备，进一步包括收集容器。
26.根据权利要求15所述的设备，进一步包括引流管。
27.根据权利要求15所述的设备，其中所述计算装置监测从由体积数据中的误差、流动数据中的误差、过滤器闭塞、呼吸状态、引流管闭塞、线路中的液体和过度流体移动组成的组中选择的至少一个额外患者参数。
28.根据权利要求15所述的设备，其中所述患者参数是被收集在收集容器中的液体流体的估计。
29.根据权利要求15-28中任一项所述的设备，其中，对所述非线性求解器进行训练。
30.根据权利要求15所述的设备，其中，所述患者参数指的是所述收集容器中的流体的体积。
CN201480024298.8A 2013-03-15 2014-03-14 流体分析仪和相关联方法 Active CN105592868B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361794712P 2013-03-15 2013-03-15
US61/794,712 2013-03-15
PCT/US2014/029319 WO2014144769A1 (en) 2013-03-15 2014-03-14 Fluid analyzer and associated methods

Publications (2)

Publication Number Publication Date
CN105592868A CN105592868A (zh) 2016-05-18
CN105592868B true CN105592868B (zh) 2018-03-09

Family

ID=54842924

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480024298.8A Active CN105592868B (zh) 2013-03-15 2014-03-14 流体分析仪和相关联方法

Country Status (2)

Country Link
EP (1) EP2968707B1 (zh)
CN (1) CN105592868B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107261219A (zh) * 2017-06-23 2017-10-20 彭殿松 一种用于胸外科胸腔防堵引流装置
CN109172888A (zh) * 2018-10-12 2019-01-11 谭泓 医用智能调压吸痰吸引负压治疗仪
WO2023209483A1 (en) * 2022-04-26 2023-11-02 Ethicon, Inc. Prolonged air leak perception

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4654029A (en) * 1982-12-13 1987-03-31 Howmedica, Inc. Electronic drainage system
US5195995A (en) * 1989-05-09 1993-03-23 Sherwood Medical Company Control vacuum to measure gas flow
US6977140B1 (en) * 1998-09-29 2005-12-20 Organ Recovery Systems, Inc. Method for maintaining and/or restoring viability of organs
CN1266179A (zh) * 2000-04-12 2000-09-13 姜家玉 测定与调控密封容器内液体流入和流出速度的方法及设备
US7976519B2 (en) * 2002-12-31 2011-07-12 Kci Licensing, Inc. Externally-applied patient interface system and method
JP4590352B2 (ja) * 2003-12-22 2010-12-01 メデラ ホールディング アーゲー 排液装置および方法
US7599894B2 (en) * 2005-03-04 2009-10-06 Hrl Laboratories, Llc Object recognition using a cognitive swarm vision framework with attention mechanisms
US8882678B2 (en) * 2009-03-13 2014-11-11 Atrium Medical Corporation Pleural drainage system and method of use
WO2010126625A1 (en) * 2009-04-30 2010-11-04 Medtronic, Inc. Patient state detection based on support vector machine based algorithm
CH705248A1 (de) * 2011-07-07 2013-01-15 Medela Holding Ag Thoraxdrainagevorrichtung mit reduziertem Gegendruck.

Also Published As

Publication number Publication date
EP2968707A1 (en) 2016-01-20
EP2968707A4 (en) 2016-11-09
CN105592868A (zh) 2016-05-18
EP2968707B1 (en) 2020-05-06

Similar Documents

Publication Publication Date Title
AU2014228778B2 (en) Fluid analyzer and associated methods
CN116705337B (zh) 一种健康数据采集及智能分析方法
Zhang et al. A context-aware mhealth system for online physiological monitoring in remote healthcare
CN105592868B (zh) 流体分析仪和相关联方法
WO2020008339A1 (en) Sensing system and method for monitoring time-dependent processes
Khan et al. Automated system design for classification of chronic lung viruses using non-linear dynamic system features and k-nearest neighbour
CN107506781A (zh) 一种基于bp神经网络的人体行为识别方法
CN109842682B (zh) 一种基于物联网的分布式环境安全学习及预警方法
WO2023097780A1 (zh) 机械通气过程中人机异步现象的分类方法和分类装置
Yan et al. Topological descriptors of gait nonlinear dynamics toward freezing-of-gait episodes recognition in parkinson’s disease
Masoudi et al. Early detection of apnea-bradycardia episodes in preterm infants based on coupled hidden Markov model
CN109119155A (zh) 基于深度学习的icu死亡危险性评估***
CN110841143A (zh) 一种用于预测输液管路状态的方法及***
CN110598789A (zh) 一种基于模糊感知器的人体疲劳状态预测方法及***
CN104980442A (zh) 一种基于元样本稀疏表示的网络入侵检测方法
CN207070264U (zh) 一种心率耳机
CN105266790A (zh) 一种心率测量的数字信号处理方法、***及设备
CN108775924A (zh) 一种基于机器人的室内空气质量监测***
Javeed et al. Physical Locomotion Detection using pre-classification for Muscles monitoring
CN117766152B (zh) 一种智能检测全身麻醉静脉管道链接脱落的方法
Dusarlapudi et al. COVID-19 patient breath monitoring and assessment with MEMS accelerometer-based DAQ-a Machine Learning Approach
Duffy et al. An investigation into smartphone based weakly supervised activity recognition systems
CN109065152A (zh) 一种健康数据分析***及方法
CN118141377B (zh) 患者的负性情绪监测***及方法
Wu et al. An Intelligent Method for Upper Limb Posture Recognition Based on Limited MEMS Data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant