Methods and apparatus for determining cardiac output

• US 9,289,133 B2
• Filed: 10/13/2010
• Issued: 03/22/2016
• Est. Priority Date: 02/10/2003
• Status: Active Grant

• Current Assignee: Massachusetts Institute of Technology, The Board of Trustees of Michigan State University (Michigan State University)
• Original Assignee: Massachusetts Institute of Technology, The Board of Trustees of Michigan State University (Michigan State University)
• Inventors: Cohen, Richard J., Mukkamala, Ramakrishna, Sherman, Derin A.
• Primary Examiner(s): BOCKELMAN, MARK
• Application Number:
  US12/903,915
• Publication Number:
  US 20110034813A1
• Time in Patent Office: 1,987 Days
• Field of Search: 600/301, 600/483, 600/500, 600/526
• US Class Current: 1/1
• CPC Class Codes: A61B 5/02007 Evaluating blood vessel condition, e.g. elasticity, compliance
• CPC Class Codes: A61B 5/02255 the pressure being controlled by plethysmographic signals, e.g. derived from optical sensors
• CPC Class Codes: A61B 5/029 Measuring or recording blood output from the heart, e.g. minute volume A61B8/065 takes precedence
• CPC Class Codes: A61B 5/352 Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
• CPC Class Codes: A61B 5/726 using Wavelet transforms

Abstract

The present invention provides methods and apparatus for determining a dynamical property of the systemic or pulmonary arterial tree using long time scale information, i.e., information obtained from measurements over time scales greater than a single cardiac cycle. In one aspect, the invention provides a method and apparatus for monitoring cardiac output (CO) from a single blood pressure signal measurement obtained at any site in the systemic or pulmonary arterial tree or from any related measurement including, for example, fingertip photoplethysmography.

According to the method the time constant of the arterial tree, defined to be the product of the total peripheral resistance (TPR) and the nearly constant arterial compliance, is determined by analyzing the long time scale variations (greater than a single cardiac cycle) in any of these blood pressure signals. Then, according to Ohm'"'"'s law, a value proportional to CO may be determined from the ratio of the blood pressure signal to the estimated time constant. The proportional CO values derived from this method may be calibrated to absolute CO, if desired, with a single, absolute measure of CO (e.g., thermodilution). The present invention may be applied to invasive radial arterial blood pressure or pulmonary arterial blood pressure signals which are routinely measured in intensive care units and surgical suites or to noninvasively measured peripheral arterial blood pressure signals or related noninvasively measured signals in order to facilitate the clinical monitoring of CO as well as TPR.

Claims

1. A method for determining cardiac output to within a constant scale factor, comprising:

measuring a physiological signal of a patient using a sensor, the physiological signal indicative of cardiovascular system activity and measured over a plurality of cardiac cycles;

determining, from the physiological signal, beat-to-beat variability of the physiological signal across multiple cardiac beats;

determining cardiac output to within a constant scale factor using the beat-to-beat variability information extracted from the physiological signal; and

presenting the cardiac output on a display of a computing device;

wherein determining beat-to-beat variability includes

constructing a cardiac contractions timing signal from the physiological signal by forming an impulse train in which each impulse is located at the onset of upstroke of each pulse;

determining an impulse response, which when convolved with the cardiac contractions timing signal, fits the physiological signal;

determining a time constant by fitting an exponential function to a tail end of the impulse response; and

determining cardiac output to within a constant scale factor using the time constant.

2. The method of claim 1 further comprises measuring the physiological signal invasively or non-invasively at any site in the systemic or pulmonary arterial tree.

3. The method of claim 1 wherein the physiological signal is arterial blood pressure.

4. The method of claim 1 wherein the physiological signal is pulmonary arterial blood pressure.

5. The method of claim 1 wherein the physiological signal is a signal related to arterial blood pressure.

6. The method of claim 1 wherein determining beat-to-beat variability further comprises estimating a function that represents the response of the physiological signal to a single cardiac contraction over a time period greater than a single cardiac cycle.

7. The method of claim 6 further comprises determining a characteristic time of the function and determining cardiac output to within a constant scale factor by dividing the magnitude of the physiological signal by the characteristic time.

8. The method of claim 1 further comprises determining vascular resistance to within a constant scale factor using the beat-to-beat variability information.

9. A method for determining a dynamical property of the systemic or pulmonary arterial tree, comprising:

measuring a physiological signal of a patient using a sensor, the physiological signal indicative of cardiovascular system activity over a plurality of cardiac cycles;

determining the timing of cardiac contractions in relation to each other, which occur in a sequence of cardiac contractions, where the determination is made from the physiological signal;

determining, from the timing of cardiac contractions and the physiological signal, a mathematical relationship between the timing of cardiac contractions in the sequence and the physiological signal occurring over the plurality of cardiac cycles;

determining a dynamical property of the systemic or pulmonary arterial tree using the relationship;

determining at least one of cardiac output and vascular resistance from the dynamical property; and

presenting the at least one of cardiac output and vascular resistance on a display of a computing device.

10. The method of claim 9 wherein the physiological signal is arterial blood pressure.

11. The method of claim 9 wherein the physiological signal is a signal related to arterial blood pressure.

12. The method of claim 9 further comprises determining the timing of cardiac contractions by measuring and analyzing a surface electrocardiogram (ECG).

13. The method of claim 9 further comprises determining at least one of cardiac output and vascular resistance to within a constant scale factor from the dynamical property.

14. The method of claim 9 further comprises fitting an exponential-like function to a portion of the relationship, estimating a time constant of the function and determining cardiac output to within a constant scale factor by dividing arterial blood pressure by the time constant.

15. A method for determining cardiac output to within a constant scale factor, comprising:

measuring a physiological signal of a patient using a sensor, the physiological signal indicative of cardiovascular system activity and measured over a plurality of cardiac cycles;

modeling the physiological signal as an output of a mathematical function;

determining an input to the mathematical function from the physiological signal, where the input is indicative of cardiac contractions;

identifying the mathematical function, which when applied to the input, fits the measured physiological signal over a plurality of cardiac cycles;

determining cardiac output to within a constant scale factor from the identified function; and

presenting the cardiac output on a display of a computing device.

16. The method of claim 15 further comprises low pass filtering the physiological signal prior to the step of determining an input.

17. The method of claim 15 wherein the physiological signal is arterial blood pressure.

18. The method of claim 15 wherein the mathematical function is a convolutional equation governed by an impulse response.

19. The method of claim 18 wherein the input to the mathematical function is an impulse train or a train of finite width pulses.

20. The method of claim 19 further comprises identifying the mathematical function from the input using autoregressive exogenous input identification.

21. The method of claim 15 further comprises estimating a time constant of the mathematical function and determining cardiac output to within a constant scale factor by dividing arterial blood pressure by the time constant.