US10828444B2 - Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters - Google Patents
Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters Download PDFInfo
- Publication number
- US10828444B2 US10828444B2 US15/546,774 US201615546774A US10828444B2 US 10828444 B2 US10828444 B2 US 10828444B2 US 201615546774 A US201615546774 A US 201615546774A US 10828444 B2 US10828444 B2 US 10828444B2
- Authority
- US
- United States
- Prior art keywords
- mus
- fitting
- region
- respiratory
- pressure
- 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, expires
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/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/085—Measuring impedance of respiratory organs or lung elasticity
-
- 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/087—Measuring breath flow
-
- 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/4836—Diagnosis combined with treatment in closed-loop systems or methods
-
- 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
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/0051—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes with alarm devices
-
- 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
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/021—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes operated by electrical means
- A61M16/022—Control means therefor
- A61M16/024—Control means therefor including calculation means, e.g. using a processor
-
- 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
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/021—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes operated by electrical means
- A61M16/022—Control means therefor
- A61M16/024—Control means therefor including calculation means, e.g. using a processor
- A61M16/026—Control means therefor including calculation means, e.g. using a processor specially adapted for predicting, e.g. for determining an information representative of a flow limitation during a ventilation cycle by using a root square technique or a regression analysis
-
- 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
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/0057—Pumps therefor
- A61M16/0063—Compressors
-
- 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
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/0003—Accessories therefor, e.g. sensors, vibrators, negative pressure
- A61M2016/0027—Accessories therefor, e.g. sensors, vibrators, negative pressure pressure meter
-
- 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
- A61M16/00—Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
- A61M16/0003—Accessories therefor, e.g. sensors, vibrators, negative pressure
- A61M2016/003—Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter
- A61M2016/0033—Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter electrical
- A61M2016/0036—Accessories therefor, e.g. sensors, vibrators, negative pressure with a flowmeter electrical in the breathing tube and used in both inspiratory and expiratory phase
-
- 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
-
- 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/50—General characteristics of the apparatus with microprocessors or computers
-
- 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
- A61M2230/00—Measuring parameters of the user
- A61M2230/40—Respiratory characteristics
- A61M2230/46—Resistance or compliance of the lungs
-
- G06F19/3481—
-
- 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
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
-
- 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
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
- G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
Definitions
- the following relates to the respiratory therapy arts, respiratory monitoring arts, and related arts.
- Various types of respiratory therapy employ a mechanical ventilator.
- the ventilator In passive patient therapy, the patient is unable to breathe, and the ventilator operates in a pressure control mode in which the ventilator performs the entire work of breathing (WoB).
- WoB work of breathing
- active patient therapy the patient can perform some of the necessary work but cannot meet the respiratory demands independently.
- ventilator operates in a pressure support mode to provide sufficient pressure to overcome any deficiency in the patient's ability to breathe.
- Volume control modes of ventilator operation are also known, in which flow rate or volume is the controlled parameter, rather than controlling pressure (although pressure limit settings may also be applied to guard against pulmonary barotrauma), and are mainly used for passive patient therapy.
- PSV pressure support mode ventilation
- assessment of the patient's work of breathing which is the clinical parameter commonly used to infer the patient's effort per breath, is facilitated by evaluating the respiratory muscle pressure, P mus (t), over the breathing cycle. More specifically, the WoB is computed by integrating P mus (t) over the inhaled volume.
- Respiratory parameters such as respiratory resistance (R) and compliance (C) may also be of interest, or may need to be determined in order to assess other parameters.
- P mus (t) in support modalities of mechanical ventilation enables the ventilator to be set such that the patient and ventilator share the mechanical work performed on the respiratory system. Quantitative assessment of P mus can be used to select the appropriate level of ventilation support in order to prevent both atrophy and fatigue of patient's respiratory muscles.
- the respiratory muscle pressure P mus (t) is typically assessed by measuring the esophageal pressure (P es ) via insertion of a balloon-tipped catheter in the patient's esophagus.
- the measured P es (t) is assumed to be a good proxy for the pleural pressure (P pt ) and can be used, in conjunction with an estimate of chest wall compliance C cw , to compute the WoB via the so-called Campbell diagram or, equivalently, via explicit computation of P mus (t) and then of WoB.
- R and C are important per se, as they provide quantitative information to the physician about the mechanical properties of the patient's respiratory system and can be used to diagnose respiratory diseases and to select the appropriate ventilation modalities and therapeutic paths.
- R and C can also be used to estimate P mus (t) as a non-invasive alternative to the use of the esophageal catheter. Assuming R and C are known, P mus (t) is suitably calculated by the following equation (known as the Equation of Motion of the Lungs):
- P y ⁇ ( t ) R ⁇ ⁇ V . ⁇ ( t ) + V ⁇ ( t ) C + P mus ⁇ ( t ) + P 0 ( 1 )
- P y (t) is the pressure measured at the Y-piece of the ventilator (also known as pressure at the mouth of the patient)
- ⁇ dot over (V) ⁇ (t) is the flow of air into and out of the patient respiratory system (measured again at the Y-piece)
- V(t) is the net volume of air delivered to the patient (measured by integrating the flow signal ⁇ dot over (V) ⁇ (t) over time)
- P 0 is a constant term to account for the pressure at the end of expiration.
- P y ⁇ ( t ) R ⁇ ⁇ V . ⁇ ( t ) + V ⁇ ( t ) C + P 0 .
- P y (t), ⁇ dot over (V) ⁇ (t), and V(t) waveforms are fully determined by the selected ventilator settings and directly measurable, so that it is straightforward to generate a sufficient data set to determine R and C.
- the value of P mus (t) varies with time over the breath cycle, and Equation (1) is not easily solved.
- Equation (1) For active patients, Equation (1) has generally been applied to non-invasively estimate P mus (t) using a two-step approach, where R and C are estimated first and then Equation (1) is applied to compute P mus (t) using the estimated values of R and C.
- Estimation of R and C may be performed by applying the flow-interrupter technique (also called End Inspiratory Pause, EIP).
- EIP End Inspiratory Pause
- the flow-interrupter technique has the disadvantage of interfering with the ventilation pattern supplied to the patient.
- the patient's respiratory muscles ought to be fully relaxed during the EIP maneuver in order for the computation of R and C to be valid, which may not always be the case.
- R and C assessed via the EIP maneuver may be different from the R and C values attained during the ventilation pattern for which P mus (t) is to be determined.
- the EIP maneuver is performed in a specific ventilation mode (Volume Assisted Control, VAC) and the resulting R and C values might not be representative of the corresponding values that determine the dynamics of the lung mechanics under other ventilation modes, such as PSV, potentially leading to error in the subsequently computed P mus (t).
- VAC Volume Assisted Control
- Equation (1) Another approach for estimating R and C in the case of an active patient is to apply least-squares fitting of Equation (1) to flow and pressure measurements under specific conditions for which the term P mus (t) is assumed to be zero.
- Some conditions for which P mus (t) could be assumed to be close to zero include: (1) periods of patient paralysis while Continuous Mandatory Ventilation (CMV) is applied; (2) periods of high Pressure Support Ventilation (PSV) levels; (3) specific portions of every pressure-supported breath that extend both during the inhalation and the exhalation phases; and (4) exhalation portions of pressure-supported breaths, where the flow signal satisfies specific conditions that are indicative of the absence of patient inspiratory efforts.
- CMV Continuous Mandatory Ventilation
- PSV Pressure Support Ventilation
- Condition (1) and (2) are undesirable clinical states that cannot be properly induced as an expedient for measuring R and C.
- the assumption of P mus (t) ⁇ 0 for Condition (3) is questionable, especially during the inhalation phase.
- Condition (4) provides only a limited amount of data for the least squares fitting procedure. In sum, it has been difficult to attain a clinically useful period of sufficient time duration for which P mus (t) ⁇ 0 is reliably achieved in an active patient in order to estimate R and C.
- a medical ventilator performs a method including: receiving measurements of pressure of air inspired by or expired from a ventilated patient operatively connected with the medical ventilator; receiving measurements of air flow into or out of the ventilated patient operatively connected with the medical ventilator; dividing a breath time interval into a plurality of fitting regions; and simultaneously estimating respiratory system's resistance and compliance or elastance, and respiratory muscle pressure by fitting to a time series of pressure and air flow samples in each region.
- the fitting includes parameterizing the respiratory muscle pressure by a continuous differentiable function, such as a polynomial function, over the fitting region.
- the fitting is to an equation of motion of the lungs in each region with monotonicity constraints of the respiratory muscle pressure and domain constraints of the respiratory parameters applied in each region.
- One advantage resides in providing non-invasive estimation of respiratory parameters including resistance, compliance, and respiratory muscle pressure.
- Another advantage resides in providing a ventilator with improved data analysis computational robustness.
- the invention may take form in various components and arrangements of components, and in various steps and arrangements of steps.
- the drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
- FIG. 1 diagrammatically shows a ventilation system.
- FIG. 2 diagrammatically shows a data analysis algorithm disclosed herein which simultaneously estimates multiple respiration parameters by approximating the respiratory muscle pressure P mus (t) by a low-order polynomial function.
- FIG. 3 plots simulated respiration waveforms over about three breaths, with sensitivity of the parameter matrix plotted in the lowermost plot of FIG. 3 .
- FIG. 4 plots simulated respiration waveforms over about three breaths with a small amplitude, high frequency pressure signal ⁇ P(t) superimposed on the ventilator-applied pressure, with sensitivity of the parameter matrix plotted in the lowermost plot of FIG. 4 .
- FIG. 5 plots normal interaction between the ventilator and a patient emulated using a computer-simulated Lung Emulator employing an ideal R, C circuit and no noise.
- FIG. 6 plots the output of the disclosed constraint optimization algorithm (top plot) and the error (bottom plot) for the data of FIG. 5 .
- FIG. 7 plots normal interaction between the ventilator and a patient emulated using a computer-simulated Lung Emulator employing an ideal R, C circuit with numerically added noise.
- FIG. 8 plots the output of the disclosed constraint optimization algorithm (top plot) and the error (bottom plot) for the data of FIG. 7 .
- PSV pressure support ventilation
- a medical ventilator system includes a medical ventilator 10 that delivers air flow at positive pressure to a patient 12 via an inlet air hose 14 . Exhaled air returns to the ventilator 10 via an exhalation air hose 16 .
- a Y-piece 20 of the ventilator system serves to couple air from the discharge end of the inlet air hose 14 to the patient during inhalation and serves to couple exhaled air from the patient into the exhalation air hose 16 during exhalation. Note that the Y-piece is sometimes referred to by other nomenclatures, such as a T-piece 20 . Not shown in FIG. 1 are numerous other ancillary components that may be provided depending upon the respiratory therapy being received by the patient 12 .
- Such ancillary components may include, by way of illustration: an oxygen bottle or other medical-grade oxygen source for delivering a controlled level of oxygen to the air flow (usually controlled by the Fraction of Inspired Oxygen (FiO 2 ) ventilator parameter set by the physician or other medical personnel); a humidifier plumbed into the inlet line 14 ; a nasogastric tube to provide the patient 12 with nourishment; and so forth.
- the ventilator 10 includes a user interface including, in the illustrative example, a touch-sensitive display component 22 via which the physician, respiratory specialist, or other medical personnel can configure ventilator operation and monitor measured physiological signals and operating parameters of the ventilator 10 . Additionally or alternatively, the user interface may include physical user input controls (buttons, dials, switches, et cetera), a keyboard, a mouse, audible alarm device(s), indicator light(s), or so forth.
- FIG. 1 in an upper portion some additional salient aspects of the ventilator system are diagrammatically illustrated in a block diagram format, including the ventilator 10 represented as a simplified block diagram, and the Y-piece 20 as a diagrammatic box with operative connections indicated by connecting arrows.
- the ventilator 10 is operating in a pressure support ventilation (PSV) mode as implemented by a controller 30 .
- PSV is an appropriate ventilation mode for an active patient who is capable for expending at least some Work of Breathing (WoB) that is, whose diaphragm and other chest muscles are acting to at least assist in operating the lungs to perform breathing.
- WoB Work of Breathing
- the controller 30 may implement various ventilation modes depending on the patient's condition and the therapy to be delivered. For example, in the case of a passive patient who is providing no WoB, the controller 30 may operate the ventilator 10 in a Pressure Control Ventilation (PCV) mode.
- PCV Pressure Control Ventilation
- the controller 30 may operate the ventilator 10 in a Pressure Control Ventilation (PCV) mode.
- PCV Pressure Control Ventilation
- volume control ventilation modes are also sometimes used, although pressure limit settings may also be applied in volume control ventilation to guard against pulmonary barotrauma.
- the ventilation controller 30 is implemented as a microprocessor with ancillary electronics such as read only memory (ROM), electronically erasable read only memory (EEPROM), flash memory, or another non-volatile memory component storing software or firmware executed by the microprocessor, random access memory (RAM) chip(s) to provide working memory, and so forth. If EEPROM, flash memory, or other updatable memory is used to store the software or firmware, then capabilities of the ventilator 10 can advantageously be updated (within the limits of its hardware components) by updating the software or firmware.
- ROM read only memory
- EEPROM electronically erasable read only memory
- flash memory or another non-volatile memory component storing software or firmware executed by the microprocessor, random access memory (RAM) chip(s) to provide working memory, and so forth.
- RAM random access memory
- the PSV controller 30 outputs a desired pressure control signal as a function of time, which is used to control a ventilator compressor 32 (e.g. a pneumatic pump, turbopump, or so forth) that generates air flow at the controlled positive pressure that is applied to the Y-piece 20 via the inlet air hose 14 .
- a ventilator compressor 32 e.g. a pneumatic pump, turbopump, or so forth
- an oxygen regulator 34 may add a controlled fraction of oxygen to the air flow to achieve a Fraction of inspired Oxygen (FiO 2 ) set by the physician, respiratory specialist, or other medical personnel who sets the configuration of the ventilator 10 for the patient 12 .
- the pressure of the ventilatory pattern may vary during the breathing cycle to provide pressure-driven or pressure-assisted inhalation and to reduce pressure to facilitate exhalation.
- the ventilator system typically further includes physiological monitoring sensors, such as an illustrative pressure sensor 40 and an illustrative flowmeter 42 .
- the pressure sensor 40 measures the pressure at the Y-piece 20 (also known as pressure at the mouth of the patient), which is denoted here as P y (t).
- the flowmeter 42 measures the air flow rate into and out of the Y-piece 20 , denoted herein as ⁇ dot over (V) ⁇ (t).
- the flowmeter 42 also directly or indirectly provides the net volume of air delivered to the patient, denoted herein as V(t), which may be directly measured or may be derived by integrating the flow rate ⁇ dot over (V) ⁇ (t) over time.
- These measured values P y (t), ⁇ dot over (V) ⁇ (t), V(t), optionally along with other information such as the ventilator settings (e.g. FiO 2 , the pressure profile delivered by the PSV control, et cetera) may be variously used by a ventilator monitor 44 to efficacy of the mechanical ventilation, to detect any deterioration of the state of the patient 12 , to detect any malfunction of the ventilator 10 , or so forth.
- the ventilator monitor 44 is implemented as a microprocessor with ancillary electronics, and may be updateable by updating the software or firmware.
- the ventilator controller 30 and the ventilator monitor 44 may be implemented by a common microprocessor, and the controller and monitor functions may be integrated at various levels for example, it is contemplated to provide feedback-based ventilation control based on the measured values P y (t), ⁇ dot over (V) ⁇ (t), V(t) or parameters derived therefrom.
- Such software or firmware may be provided in the form of a non-transitory storage medium storing instructions readable and executable by the microprocessor of the ventilator monitor 44 to perform the monitoring functionality.
- the non-transitory storage medium may, for example, comprise a flash memory, optical disk, hard disk drive, or other storage medium.
- Equation (2) the system of equations represented by Equation (2) has more unknowns (N+2 unknowns) than equations (N equations), and hence is an underdetermined problem that cannot be solved because it has an infinite number of solutions, only one of which is the true “physical” solution.
- Equation (2) Due to being underdetermined, the set of equations represented by matrix Equation (2) is very sensitive to measurement noise, unknown disturbances and unmodeled dynamics. Problematically, the noise is on the same time scale as the variations in the measured signals P y (t), ⁇ dot over (V) ⁇ (t), V(t) and in the fitted respiratory muscle pressure P mus (t). Thus, even if the underdetermined nature of the simultaneous estimation problem is somehow overcome, the resulting parameter values tend to be noisy and hence of limited clinical value.
- ⁇ P(t) a relatively high-frequency and small-amplitude pressure signal
- ⁇ P(t) a relatively high-frequency and small-amplitude pressure signal
- this can be done by adding a small-amplitude sinusoidal ⁇ P(t) to the controlled pressure signal output by the controller 30 using a signal combiner 52 prior to its input to the ventilator compressor 32 .
- the amplitude of ⁇ P(t) is preferably chosen to be low enough to not appreciably impact the therapeutic value of the PSV signal output by the controller 30 .
- the frequency of ⁇ P(t) is preferably high enough to be significantly higher than the breath frequency (e.g. typically a few breaths per minute corresponding to a frequency of, e.g., about 0.2 Hz for 5-sec breaths).
- Equation (2) the underdeterminancy of the set of equations represented by matrix Equation (2) is addressed in embodiments disclosed herein (either with or without the optional superimposed ⁇ P(t)) by solving Equation (2) in fitting region(s) 60 for which the respiratory muscle pressure P mus (t), although not assumed to be zero, is reasonably assumed to have some constraining characteristic(s) that enable the number of parameters to be reduced sufficiently to make Equation (2) overdetermined.
- Equation (2) The simultaneous estimation of the R, C and P mus (t) characterizing one breath (represented, without loss of generality, by N recorded time samples) by Equation (2) is an underdetermined problem, since it requires the computation of N+2 unknowns (N values for the N time samples of P mus (t), plus an additional unknown for R, and an additional unknown for C) from N equations corresponding to the N time samples.
- N equations are not independent. Rather, it can be expected that the value of P mus (t) for neighboring samples should be continuous. In some regions, it may be reasonably assumed that P mus (t) is monotonically increasing, being flat, or monotonically decreasing.
- This approximation is used to construct a least squares (LS) problem over a time window of s samples (where s ⁇ N) in which the unknowns are R, C, and a 0 , . . . , a n .
- LS least squares
- the time interval covered by the N samples is divided into fitting regions in which P mus (t) is monotonically increasing, monotonically decreasing, or being flat over the entire fitting region.
- a quadratic program can be constructed leveraging the known monotonicity in the region. This ensures efficient determination of a unique solution.
- the fitting regions 60 are chosen to be small enough for P mus (t) to be well-fitted by a polynomial approximation.
- the estimation of R, C, and P mus (t) at each time 1, . . . ,N is obtained by solving a LS problem over a window of lengths. For the usual case in which s ⁇ N, the window slides forward in time (that is, the window of width s is applied to successive increments of width s in the time series of samples 1, . . . ,N).
- V [V( 1 ) V( 2 ) . . . V(s)] T
- V [V( 1 ) V( 2 ) . . . V(s)] T
- matrix ⁇ is an s ⁇ (n+3) matrix given by:
- an iterative least squares approximation approach such as gradient descent or Levenberg-Marquardt can be used to solve Equation (3) for the parameters ⁇ .
- the illustrative approach employs a polynomial approximation of order n of P mus (t) over the time window of width s.
- the approach can be generalized to approximating P mus (t) over the time window of width s by any continuous function that is smooth over the window of width s (i.e. that is differentiable over the window of width s).
- Other contemplated continuous and smooth approximation functions include spline functions, e.g. cubic spline functions.
- the normal interaction between the ventilator and a patient is emulated using a computer-simulated Lung Emulator.
- This normal interaction gives rise to waveforms of flow and volume, plotted in FIG. 3 , that make the data matrix ill-conditioned.
- the parameters estimated via Equation (3) are sensitive to noise or error in the measured data.
- the sensitivity of the parameter matrix is plotted in the lowermost plot of FIG. 3 . Note in this plot that the sensitivity ordinate ranges [0, 200,000].
- this noise is optionally counteracted by superimposition of the low amplitude, high frequency component ⁇ P(t). To illustrate, FIG.
- the signal generator 50 and signal combiner 52 can be implemented in software or firmware as part of the software or firmware of the ventilator controller 30 , or the signal generator 50 and signal combiner 52 can be components separate from the ventilator controller 30 , e.g. a voltage-controlled oscillator (VCO) circuit outputting the signal ⁇ P(t), and an op-amp-based signal combiner or other signal combiner implemented in hardware.
- VCO voltage-controlled oscillator
- Equation (1) can be modified to have quadratic characteristics as follows:
- Equation (4) ( R 0 + R 1 ⁇ ⁇ V . ⁇ ( t ) ⁇ ) ⁇ V . ⁇ ( t ) + ( 1 C 0 + V ⁇ ( t ) C 1 ) ⁇ V ⁇ ( t ) + P mus ⁇ ( t ) + P 0 ( 4 ) Equation (4) is characterized by a flow-dependent resistance R 0 +R 1 ⁇
- ⁇ dot over (V) ⁇ 0 [ ⁇ dot over (V) ⁇ ( 1 ) ⁇ dot over (V) ⁇ ( 2 ) . . . ⁇ dot over (V) ⁇ (s)] T
- ⁇ dot over (V) ⁇ 1 [ ⁇ dot over (V) ⁇ ( 1 )
- V 0 [V( 1 ) V( 2 ) . . . V(s)] T
- V 1 [V 2 ( 1 ) V 2 ( 2 ) . . . V 2 (s)] T
- matrix ⁇ is an s ⁇ (n+5) matrix given by:
- an iterative least squares approximation approach such as gradient descent or Levenberg-Marquardt can be used to solve Equation (5) for the parameters ⁇ NL .
- each fitting region 60 is chosen so that P mus (t) is monotonic (either monotonically increasing, or monotonically decreasing) in the entire fitting region.
- inequality constraints on the possible values of P mus (t), and domain constraints that R and C can take are introduced based on physiological considerations, such that the least squares (LS) solution becomes unique.
- the constraints are cast in linear form and define an objective function to be minimized of LS type, so that the mathematical formulation of the optimization problem to be solved falls into the category of quadratic programming. Not only is the uniqueness of solution now guaranteed, but also the routine to solve the program can be very efficient since quadratic programming is a mature mathematical technology.
- robustness of the resulting method to estimate R, C, and P mus (t) is further improved by the introduction of equality constraints.
- Robustness is advantageous for practical applications due to uncertainties and non-ideal factors that can affect the application (measurement noise, unknown disturbances, nonlinearities, unmodeled dynamics).
- Equality constraints on the values of P mus (t) are used to reduce the number of unknowns to describe P mus (t), hence making the overall estimation more robust.
- the normal interaction between the ventilator and a patient is emulated using a computer simulated Lung Emulator.
- Equation (7) the respiratory system's compliance C is replaced by the elastance E according to the relationship
- Eigenvalue decomposition of the quadratic matrix that can be constructed from the objective function using real data demonstrates that the problem is fully determined under Constraints (8)-(10). All the eigenvalues are negative but two, which are zero. For the quadratic problem to have a unique solution all the eigenvalues should be strictly negative. The eigenvectors associated with the zero eigenvalues are, however, minimizing directions forbidden by the given constraints, so that the underdeterminacy of the LS simultaneous estimation of R, C, and P mus (t) is overcome.
- Constraints in addition to, or instead of, the Constraints (8)-(10) are contemplated.
- the input to the algorithm is the set of measured P y (t), ⁇ dot over (V) ⁇ (t), and V(t) over a complete breath, where again V(t) is suitably obtained by integration of ⁇ dot over (V) ⁇ (t).
- the output includes a value for each of R, C (or E), and the waveform P mus (t) for the entire breath.
- FIGS. 5-8 experiments on simulated respiratory data indicate that the above-described quadratic program algorithm provides suitable estimates of R, C and P mus (t) when the pressure and flow data come from an ideal R, C circuit with and without additive noise corrupting the measurements.
- FIG. 5 shows the simulated respiration data for the ideal R, C circuit simulator with no noise in the signals
- FIG. 6 shows the output of the quadratic program algorithm (top plot) and the error (bottom plot; note that the ordinate of the error plot has a range [0,10 ⁇ 13 ] so that negligible error is observed throughout).
- FIGS. 7 and 8 show the same experiment as FIGS. 5 and 6 , but this time with numerically-generated noise added. While some error is observed due to the noise, the fit is still fairly accurate.
- the disclosed techniques can be combined, e.g. the quadratic program (Equation (7) with Constraints (8)-(10)) can be performed in conjunction with a parameterization of P mus (t), for example as described with reference to FIG. 2 or some other parameterization.
- the solid lines are the LS estimates of R and C obtained by measuring P es (t), which requires an invasive catheter (knowledge of P es (t) permits estimation of resistance and compliance of the respiratory system via LS, with no underdeterminacy issues).
- the quadratic program algorithm is therefore able to provide non-invasively the same R and C estimates that the current state-of-the-art obtains invasively.
Abstract
Description
where Py(t) is the pressure measured at the Y-piece of the ventilator (also known as pressure at the mouth of the patient), {dot over (V)}(t) is the flow of air into and out of the patient respiratory system (measured again at the Y-piece), V(t) is the net volume of air delivered to the patient (measured by integrating the flow signal {dot over (V)}(t) over time), and P0 is a constant term to account for the pressure at the end of expiration.
For the passive patient, Py(t), {dot over (V)}(t), and V(t) waveforms are fully determined by the selected ventilator settings and directly measurable, so that it is straightforward to generate a sufficient data set to determine R and C. By contrast, in the case of an active patient who is providing some WoB, the value of Pmus(t) varies with time over the breath cycle, and Equation (1) is not easily solved.
Y=Xθ (2)
where:
Y=[Py(1) Py(2) . . . Py(N)]T Pressure at Y-piece at
{dot over (V)}=[{dot over (V)}(1) {dot over (V)}(2) . . . {dot over (V)}(N)]T Flow rate at
V=[V(1) V(2) . . . V(N)]T Net air volume at
θ=[
and matrix X is an (N+2)×N matrix given by X=[{dot over (V)} V IN], where IN is an N×N identity matrix. By solving the system of equations Y=Xθ for the parameter vector θ, the resistance R, compliance C, and respiratory muscle pressure Pmus(t) can be obtained. However, the system of equations represented by Equation (2) has more unknowns (N+2 unknowns) than equations (N equations), and hence is an underdetermined problem that cannot be solved because it has an infinite number of solutions, only one of which is the true “physical” solution.
Y=χϕ (2a)
where:
In the above notation, the first sample in the window of width s is designated without loss of generality as sample t=1, so that the last sample in the window is designated as sample t=s. Matrix Equation (2a) thus represents a set of s equations with n+3 unknowns, and is overdetermined so long as s>(n+3). More typically, s>>n. For example, in one illustrative example n=2 (quadratic approximation for Pmus(t)), the sampling rate is 100 Hz, and the window is 0.6 sec long corresponding to s=60.
ϕ=(χTχ)−1χT Y (3)
Alternatively, an iterative least squares approximation approach such as gradient descent or Levenberg-Marquardt can be used to solve Equation (3) for the parameters ϕ.
Equation (4) is characterized by a flow-dependent resistance R0+R1·|{dot over (V)}(t)| and a volume-dependent elastance
The parameters to be estimated are now R0, R1, C0, C1, and Pmus(t). The least squares (LS) problem to be solved (with polynomial approximation of Pmus(t), i.e. corresponding to Equation (2a)) becomes:
Y=χ NLϕNL (5)
where:
Assuming an overdetermined set of equations, the matrix Equation (5) can be solved in the least squares sense according to:
ϕNL=(χNL TχNL)−1χNL T Y (6)
Alternatively, an iterative least squares approximation approach such as gradient descent or Levenberg-Marquardt can be used to solve Equation (5) for the parameters ϕNL.
J=Σ t=1 N(P y(t)−R{dot over (V)}(t)−EV(t)−P mus(t))2 (7)
In the objective J of Equation (7), the respiratory system's compliance C is replaced by the elastance E according to the relationship
The objective function J is minimized with respect to the parameters R, C (or E), and Pmus(1), . . . , Pmus(N), subject to inequality constraints capturing the known monotonic regions of Pmus(t). This problem can be cast as a quadradic program by minimizing J subject to the following inequality constraints:
where the time t=m is the “turning point”, that is, the point at which Pmus(t) goes from being monotonically decreasing (for t=1, . . . , m) to being monotonically increasing (for t=m+1, . . . ,N). Said another way, Pmus(m) is the time at which Pmus(t) reaches its minimum value. Optionally, the quadratic program can include additional constraints based on physiological knowledge. For example, if there is some known minimum respiratory muscle pressure Pmin and/or some known maximum respiratory muscle pressure Pmax (for example, it may in some instances be assumed that Pmax=0 as the diaphragm and chest muscles cannot act to apply positive pressure to the lungs), then the following inequalities can be added:
Similar limit (domain) constraints may optionally be placed on R and C:
R min ≤R≤R max (10)
E min ≤E≤E max
Claims (18)
P mus(t)=a 0 +a 1 t+ . . . +a n t n
P mus(t)=a 0 +a 1 t+ . . . +a n t n.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/546,774 US10828444B2 (en) | 2015-02-12 | 2016-01-12 | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562115182P | 2015-02-12 | 2015-02-12 | |
US15/546,774 US10828444B2 (en) | 2015-02-12 | 2016-01-12 | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters |
PCT/IB2016/050115 WO2016128846A1 (en) | 2015-02-12 | 2016-01-12 | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters |
Publications (2)
Publication Number | Publication Date |
---|---|
US20180001041A1 US20180001041A1 (en) | 2018-01-04 |
US10828444B2 true US10828444B2 (en) | 2020-11-10 |
Family
ID=55272529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/546,774 Active 2037-05-14 US10828444B2 (en) | 2015-02-12 | 2016-01-12 | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters |
Country Status (6)
Country | Link |
---|---|
US (1) | US10828444B2 (en) |
EP (1) | EP3256197B1 (en) |
JP (1) | JP6717843B2 (en) |
CN (1) | CN107205694B (en) |
RU (1) | RU2712040C2 (en) |
WO (1) | WO2016128846A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11247008B1 (en) | 2020-08-05 | 2022-02-15 | Effortless Oxygen, Llc | Flow triggered gas delivery |
US11318276B2 (en) | 2020-08-05 | 2022-05-03 | Effortless Oxygen, Llc | Flow triggered gas delivery |
US11420007B2 (en) | 2020-08-05 | 2022-08-23 | Effortless Oxygen, Llc | Flow triggered gas delivery |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017068464A1 (en) | 2015-10-19 | 2017-04-27 | Koninklijke Philips N.V. | Anomaly detection device and method for respiratory mechanics parameter estimation |
JP6831377B2 (en) | 2015-11-02 | 2021-02-17 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Breath-by-breath reassessment of patient lung parameters to improve estimated performance |
CN109069030B (en) | 2016-02-18 | 2021-10-19 | 皇家飞利浦有限公司 | Enhancing respiratory parameter estimation and asynchronous detection algorithms via use of central venous pressure manometry |
EP3544502A1 (en) * | 2016-10-26 | 2019-10-02 | Koninklijke Philips N.V. | System and methods for estimation of respiratory muscle pressure and respiratory mechanics using p0.1 maneuver |
EP3558431A1 (en) | 2016-12-20 | 2019-10-30 | Koninklijke Philips N.V. | Determining respiratory mechanic parameters in the presence of intrinsic positive end-expiratory pressure |
CN109350063B (en) * | 2018-12-03 | 2021-03-30 | 北京航空航天大学 | Respiratory mechanics parameter detection device and method suitable for chronic obstructive pulmonary disease monitoring |
JP6981551B2 (en) * | 2019-02-19 | 2021-12-15 | Jfeスチール株式会社 | Continuous casting machine control method, continuous casting machine control device, and slab manufacturing method |
WO2020170563A1 (en) * | 2019-02-19 | 2020-08-27 | Jfeスチール株式会社 | Control method for continuous casting machine, control device for continuous casting machine, and method for manufacturing slab |
CN110537917A (en) * | 2019-08-15 | 2019-12-06 | 郴州市第一人民医院 | mechanical ventilation intelligent monitoring system and monitoring method based on respiratory mechanics |
DE102020133460A1 (en) | 2020-01-07 | 2021-07-08 | Drägerwerk AG & Co. KGaA | Method and signal processing unit for determining a pneumatic measure using a lung mechanical model and a progression model |
CN113546260B (en) * | 2021-06-07 | 2023-08-22 | 中国人民解放军总医院第一医学中心 | Method, apparatus and storage medium for monitoring for use anomalies in transnasal high flow |
CN116796656B (en) * | 2023-05-17 | 2024-03-19 | 中国人民解放军总医院第一医学中心 | Method for estimating lumped parameter respiratory system model parameters |
CN117323525B (en) * | 2023-12-01 | 2024-02-23 | 南京沪家医疗科技有限公司 | Pressure control method and device of breathing machine |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030010339A1 (en) * | 1999-02-03 | 2003-01-16 | University Of Florida | Method and apparatus for nullifying the imposed work of breathing |
US20030159695A1 (en) * | 2000-04-26 | 2003-08-28 | Magdy Younes | Method and apparatus for determining respiratory system resistance during assisted ventilation |
JP2008000436A (en) | 2006-06-23 | 2008-01-10 | Air Water Safety Service Inc | Lung compliance estimation apparatus and estimation method and artificial respiratory apparatus provided with estimation apparatus |
US7322937B2 (en) * | 2002-11-20 | 2008-01-29 | Maquet Critical Care Ab | Method and breathing apparatus for assessing pulmonary stress |
US20090022192A1 (en) | 2004-02-13 | 2009-01-22 | Tadashi Itoh | Laser Device and Lasing Method |
US7484508B2 (en) | 2002-06-27 | 2009-02-03 | Yrt Limited | Method and device for monitoring and improving patient-ventilator interaction |
US20090221926A1 (en) | 2006-05-12 | 2009-09-03 | Magdy Younes | Method and Device for Generating of a Signal that Reflects Respiratory Efforts in Patients on Ventilatory Support |
US20100071696A1 (en) | 2008-09-25 | 2010-03-25 | Nellcor Puritan Bennett Llc | Model-predictive online identification of patient respiratory effort dynamics in medical ventilators |
US20100081119A1 (en) * | 2008-09-30 | 2010-04-01 | Nellcor Puritan Bennett Llc | Configurable respiratory muscle pressure generator |
US20110237970A1 (en) | 2008-12-10 | 2011-09-29 | Koninklijke Philips Electronics N.V. | Determining elastance and resistance |
US20130152934A1 (en) * | 2010-08-27 | 2013-06-20 | Resmed Limited | Adaptive cycling for respiratory treatment apparatus |
US20150314090A1 (en) * | 2014-04-30 | 2015-11-05 | Landseed Hospital | Alarming method for ventilator and ventilator alarm system |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5803066A (en) * | 1992-05-07 | 1998-09-08 | New York University | Method and apparatus for optimizing the continuous positive airway pressure for treating obstructive sleep apnea |
GB9525994D0 (en) * | 1995-12-20 | 1996-02-21 | Univ Manitoba | Improvements in delivery of assist modes of mechanical ventilation |
US6257234B1 (en) * | 1998-08-21 | 2001-07-10 | Respironics, Inc. | Apparatus and method for determining respiratory mechanics of a patient and for controlling a ventilator based thereon |
JP2000175886A (en) * | 1998-12-14 | 2000-06-27 | Nippon Koden Corp | Method and apparatus for processing ventilation data |
US7520279B2 (en) * | 2001-07-19 | 2009-04-21 | Resmed Limited | Pressure support ventilation of patients |
JP3860803B2 (en) * | 2003-06-03 | 2006-12-20 | エア・ウォーター防災株式会社 | Control device for gas supply servomechanism for ventilator |
US20070221221A1 (en) * | 2006-02-23 | 2007-09-27 | Cooke Richard H | Ventilator for Rapid Response to Respiratory Disease Conditions |
EP2091429B1 (en) * | 2006-11-16 | 2010-11-10 | Hamilton Medical AG | Method and device for determining the peep during the respiration of a patient |
EP1972274B1 (en) * | 2007-03-20 | 2015-12-30 | Drägerwerk AG & Co. KGaA | Method and apparatus for determining the resistance of the respiratory system of a patient |
EP2349423B1 (en) * | 2008-10-24 | 2014-02-26 | Hospitech Respiration Ltd. | Ventilation system |
WO2010070498A1 (en) * | 2008-12-19 | 2010-06-24 | Koninklijke Philips Electronics, N.V. | System and method for treating lung disease using positive pressure airway support |
JP5665869B2 (en) * | 2009-09-01 | 2015-02-04 | コーニンクレッカ フィリップス エヌ ヴェ | System and method for quantifying lung compliance in a self-ventilating subject |
CN102770070B (en) * | 2009-12-28 | 2015-11-25 | 佛罗里达大学研究基金会有限公司 | For the system and method for real-time assessment mechanics of lung |
MX2013001557A (en) * | 2010-08-13 | 2013-06-28 | Respiratory Motion Inc | Devices and methods for respiratory variation monitoring by measurement of respiratory volumes, motion and variability. |
WO2013098714A1 (en) * | 2011-12-30 | 2013-07-04 | Koninklijke Philips Electronics N.V. | Airway impedance measurement integrated with respiratory treatment devices |
CN104203093B (en) * | 2012-03-30 | 2017-10-03 | 皇家飞利浦有限公司 | For breathing system and method and the closed loop controller that energy is assessed in real time |
US10165966B2 (en) * | 2013-03-14 | 2019-01-01 | University Of Florida Research Foundation, Incorporated | Methods and systems for monitoring resistance and work of breathing for ventilator-dependent patients |
-
2016
- 2016-01-12 WO PCT/IB2016/050115 patent/WO2016128846A1/en active Application Filing
- 2016-01-12 EP EP16702201.1A patent/EP3256197B1/en active Active
- 2016-01-12 RU RU2017131740A patent/RU2712040C2/en active
- 2016-01-12 US US15/546,774 patent/US10828444B2/en active Active
- 2016-01-12 JP JP2017542069A patent/JP6717843B2/en active Active
- 2016-01-12 CN CN201680009858.1A patent/CN107205694B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030010339A1 (en) * | 1999-02-03 | 2003-01-16 | University Of Florida | Method and apparatus for nullifying the imposed work of breathing |
US20030159695A1 (en) * | 2000-04-26 | 2003-08-28 | Magdy Younes | Method and apparatus for determining respiratory system resistance during assisted ventilation |
US7484508B2 (en) | 2002-06-27 | 2009-02-03 | Yrt Limited | Method and device for monitoring and improving patient-ventilator interaction |
US7322937B2 (en) * | 2002-11-20 | 2008-01-29 | Maquet Critical Care Ab | Method and breathing apparatus for assessing pulmonary stress |
US20090022192A1 (en) | 2004-02-13 | 2009-01-22 | Tadashi Itoh | Laser Device and Lasing Method |
US20090221926A1 (en) | 2006-05-12 | 2009-09-03 | Magdy Younes | Method and Device for Generating of a Signal that Reflects Respiratory Efforts in Patients on Ventilatory Support |
JP2008000436A (en) | 2006-06-23 | 2008-01-10 | Air Water Safety Service Inc | Lung compliance estimation apparatus and estimation method and artificial respiratory apparatus provided with estimation apparatus |
US20100071696A1 (en) | 2008-09-25 | 2010-03-25 | Nellcor Puritan Bennett Llc | Model-predictive online identification of patient respiratory effort dynamics in medical ventilators |
US20100081119A1 (en) * | 2008-09-30 | 2010-04-01 | Nellcor Puritan Bennett Llc | Configurable respiratory muscle pressure generator |
US20110237970A1 (en) | 2008-12-10 | 2011-09-29 | Koninklijke Philips Electronics N.V. | Determining elastance and resistance |
US20130152934A1 (en) * | 2010-08-27 | 2013-06-20 | Resmed Limited | Adaptive cycling for respiratory treatment apparatus |
US20150314090A1 (en) * | 2014-04-30 | 2015-11-05 | Landseed Hospital | Alarming method for ventilator and ventilator alarm system |
Non-Patent Citations (11)
Title |
---|
Antonio Albanese: "Physiology-based Mathematical Models for the Intensive Care Unit: Application to Mechanical Ventilation", Jan. 1, 2014 XP055332137, ISBN: 978-1-303-94555-7 (Year: 2014) * |
ANTONIO ALBANESE: "Physiology-based Mathematical Models for the Intensive Care Unit: Application to Mechanical Ventilation", PROQUEST DISSERTATIONS PUBLISHING, 1 January 2014 (2014-01-01), XP055332137, ISBN: 978-1-303-94555-7, Retrieved from the Internet <URL:https://academiccommons.columbia.edu/download/fedora_content/download/ac:199653/CONTENT/Albanese_columbia_0054D_12039.pdf> * |
Brochard et al: "Insiratory Pressure Support Prevents Diaphragmatic Fatigue During Weaning From Mechanical Ventilation"; Am. Rev. Respier. Dis. vol. 139, pp. 513-521, 1989. |
Heyer et al: "Non-Invasive Detection of Respiratory Muscles Activity During Assisted Ventilation"; C.R. Biologies 325 (4), pp. 383-391, 2002. |
Iotti et al:"Noninvasive Evaluation of Instantaneous Total Mechanical Activity of the Respiratory Activity of the Respiratory Muscles During Pressure Support Ventilation"; Chest, vol. 108 (1):208-215, 1995. |
Kondili et al: "Estimatin of Inspiratory Muscle Pressure in Critically Ill Patients"; Intensive Care Med (2010), vol. 36(4), pp. 648-655, 2010. |
Lutchen et al: "Optimal Ventilation Waveforms for Estimating Low-Frequency Respiratory Impedance"; J. Appl. Physiology, vol. 75, pp. 478-488, 1993. |
Schott et al: "Real-Time Computation of a Patient's Respiratory Effort During Ventilation"; Journal of Clinical Monitoring and Computing (2006), vol. 20, pp. 193-200. |
Werner Nikischin, Tilo Gerhardt, Ruth Everett, and Eduardo Bancalari,"A New Method to Analyze Lung Compliance When Pressure—Volume Relationship is Nonlinea", 1998 (Year: 1998). * |
Yamada et al: "Respiratory Muscle Pressure Analysis in Pressure-Support Ventilation"; J. Appl, Physiol, vol. 77(5), pp. 2237-2243, 1994. |
Younes et al: "A Method for Monitoring and Improving Patient:Ventilator Interaction"; Intensive Care Med, (2007), vol. 33, pp. 1337-1346. |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11247008B1 (en) | 2020-08-05 | 2022-02-15 | Effortless Oxygen, Llc | Flow triggered gas delivery |
US11318276B2 (en) | 2020-08-05 | 2022-05-03 | Effortless Oxygen, Llc | Flow triggered gas delivery |
US11420007B2 (en) | 2020-08-05 | 2022-08-23 | Effortless Oxygen, Llc | Flow triggered gas delivery |
Also Published As
Publication number | Publication date |
---|---|
RU2017131740A3 (en) | 2019-06-14 |
EP3256197A1 (en) | 2017-12-20 |
RU2712040C2 (en) | 2020-01-24 |
EP3256197B1 (en) | 2019-07-24 |
US20180001041A1 (en) | 2018-01-04 |
JP6717843B2 (en) | 2020-07-08 |
WO2016128846A1 (en) | 2016-08-18 |
RU2017131740A (en) | 2019-03-12 |
CN107205694B (en) | 2021-04-06 |
CN107205694A (en) | 2017-09-26 |
JP2018506354A (en) | 2018-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10828444B2 (en) | Simultaneous estimation of respiratory parameters by regional fitting of respiratory parameters | |
RU2737295C2 (en) | Apparatus for mechanic artificial pulmonary ventilation and respiratory monitoring | |
CN109803708B (en) | Breathing apparatus and ventilator apparatus | |
EP2799008B1 (en) | Method for continuous and non-invasive determination of effective lung volume and cardiac output | |
US20130006133A1 (en) | Methods and systems for monitoring volumetric carbon dioxide | |
US20130006134A1 (en) | Methods and systems for monitoring volumetric carbon dioxide | |
JP6487425B2 (en) | Calculation of respiratory work based on noninvasive estimation of intrathoracic pressure and / or noninvasive estimation of intrathoracic pressure | |
US11191447B2 (en) | Breath by breath reassessment of patient lung parameters to improve estimation performance | |
US8876728B2 (en) | System and method for quantifying lung compliance in a self-ventilating subject | |
US20170367617A1 (en) | Probabilistic non-invasive assessment of respiratory mechanics for different patient classes | |
JP2018536510A (en) | Simultaneous estimation of respiratory mechanics and patient effort by parameter optimization | |
US11738161B2 (en) | System and method for adaptive scheduling of pause maneuvers used for estimation of compliance and/or resistance during mechanical ventilation | |
JP7168560B2 (en) | System and method for estimation of respiratory muscle pressure and ventilation dynamics using P0.1 maneuver | |
CN107690310B (en) | Non-invasive method for monitoring the respiratory state of a patient via continuous parameter estimation | |
WO2016132279A1 (en) | Enhancement of simultaneous estimation of respiratory parameters via superimposed pressure signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: KONINKLIJKE PHILIPS N.V., NETHERLANDS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ALBANESE, ANTONIO;VICARIO, FRANCESCO;WANG, DONG;AND OTHERS;SIGNING DATES FROM 20160111 TO 20170725;REEL/FRAME:043117/0123 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE AFTER FINAL ACTION FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: AWAITING TC RESP., ISSUE FEE NOT PAID |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |