Methods and apparatus for determining arterial pulse wave velocity

• US 10,349,838 B2
• Filed: 08/12/2011
• Issued: 07/16/2019
• Est. Priority Date: 08/12/2010
• Status: Active Grant

• Current Assignee: The Board of Trustees of Michigan State University (Michigan State University)
• Original Assignee: The Board of Trustees of Michigan State University (Michigan State University)
• Inventors: Mukkamala, Ramakrishna, Xu, Da, Zhang, Guanqun, Gao, Mingwu, Moslehpour, Mohsen
• Primary Examiner(s): Mallari, Patricia
• Assistant Examiner(s): Portillo, Jairo H
• Application Number:
  US13/816,427
• Publication Number:
  US 20130184595A1
• Time in Patent Office: 2,895 Days
• Field of Search: None
• CPC Class Codes: A61B 5/02125 of pulse wave propagation time
• CPC Class Codes: A61B 5/0285 Measuring or recording phase velocity of blood waves
• 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/7235 Details of waveform analysis A61B5/349 takes precedence
• CPC Class Codes: A61B 5/7278 Artificial waveform generation or derivation, e.g. synthesising signals from measured signals A61B5/327 takes precedence

Abstract

Methods are presented for determining pulse transit time (PTT) and/or pulse wave velocity (PWV) of a subject by application of parametric system identification to proximal and distal arterial waveforms. The two waveforms are measured from the subject. A system is defined that relates the proximal arterial waveform to the distal arterial waveform (or vice versa) in terms of the unknown parameters of a parametric mathematical model. The model parameters are determined from the measured waveforms using system identification. PTT between the proximal and distal arterial sites is then determined from the system model. PWV may also be determined by dividing the distance between measurement sites (D) by PTT.

Claims

1. A method for determining an arterial pulse transit time of a subject, comprising:

measuring a proximal waveform indicative of blood pressure at a proximal site of the subject using a first sensor;

measuring a distal waveform indicative of blood pressure at a distal site of the subject using a second sensor;

defining a dynamic system that relates one of the waveforms to the other waveform in terms of unknown parameters of a parametric mathematical model that assumes a specific structure for the dynamic system, where one of the unknown parameters accounts for wave reflection and another of the unknown parameters accounts for pulse transit time;

determining the unknown parameters of the mathematical model from the measured proximal waveform and the measured distal waveform; and

determining a pulse transit time between the proximal site and the distal site using the determined parameters from the mathematical model, where the steps of determining the unknown parameters of the mathematical model and determining a pulse transit time are implemented by a computer processor of a computing device.

2. The method of claim 1 further comprises determining a distance between the proximal and distal sites and determining a pulse wave velocity by dividing the distance by the pulse transit time.

3. The method of claim 1 further comprises determining the unknown parameters by fitting one of the proximal waveform or the distal waveform to the other waveform.

4. The method of claim 1 wherein the parametric mathematical model is a parametric physical tube-load model transfer function in which the pulse transit time is an explicit parameter of the model.

5. The method of claim 4 further comprises determining the parameters by finding the transfer function, which when applied to one of the waveforms, best fits the other waveform in a least squares sense.

6. The method of claim 1 wherein the mathematical model is a parametric black-box model that is characterized by non-physical parameters.

7. The method of claim 6 wherein the parametric black-box model is further defined as an impulse response.

8. The method of claim 7 further comprises determining the non-physical parameters by finding the impulse response, which when convolved with the proximal arterial waveform, best fits the distal arterial waveform in the least squares sense.

9. The method of claim 8 wherein the pulse transit time is determined as a time delay of the impulse response.

10. The method of claim 7 wherein the mathematical model is as an autoregressive exogenous input equation:

11. The method of claim 10 further comprises estimating the unknown parameters using a least squares method.

12. The method of claim 11 further comprises, for low heart rate, dividing a beat of the waveforms into two segments, determining the parameters of the transfer function for each of the two segments, and taking an average of the pulse transit time parameters for each of the two segments.

13. A method for determining an arterial pulse wave velocity of a subject, comprising:

measuring a proximal waveform indicative of the arterial pulse at a proximal site of the subject using a first sensor;

measuring a distal waveform indicative of the arterial pulse at a distal site of the subject using a second sensor;

defining a dynamic system for the subject that relates one of the waveforms to the other waveform in terms of unknown parameters of a physical tube-load model transfer function in which pulse transit time is an explicit parameter of the model;

determining the parameters of the transfer function from the measured proximal and the measured distal waveforms; and

determining a pulse wave velocity from the pulse transit time parameter, where the steps of determining the parameters of the transfer function and determining a pulse transit velocity are implemented by a computer processor of a computing device.

14. The method of claim 13 wherein the waveforms are indicative of at least one of blood pressure, blood velocity, or blood flow rate.

15. The method of claim 13 further comprises measuring the proximal waveform from a carotid artery and measuring the distal waveform from a femoral artery.

16. The method of claim 13 wherein the physical tube-load model having a tube which is uniform and frictionless and represents a wave travel path between the proximal and distal arterial measurement sites, and a load distal to the distal measurement site, where the load has a frequency-dependence.

17. The method of claim 13 wherein the parameters are determined by finding the transfer function, which when applied to one of the waveforms, best fits the other waveform using a least square technique.

18. The method of claim 13 wherein determining pulse wave velocity further comprises determining a distance between the two measurement sites of the waveforms and determining the pulse wave velocity by dividing the distance by the pulse transit time parameter.

19. A method for tracking blood pressure of a subject, comprising:

measuring a proximal waveform indicative of the arterial pulse at a proximal site of the subject using a first sensor;

measuring a distal waveform indicative of the arterial pulse at a distal site of the subject using a second sensor;

defining a dynamic system that relates the proximal arterial waveform to the distal arterial waveform in terms of a mathematical model, where one parameter of the mathematical model accounts for wave reflection and another parameter of the mathematical model accounts for pulse transit time;

determining the mathematical model for the subject from the measured waveforms;

determining a pulse transit time between the proximal site and the distal site from the mathematical model, where the steps of determining the mathematical model and determining a pulse transit time are implemented by a computer processor of a computing device; and

tracking blood pressure by mapping the pulse transit time through a calibration curve that relates pulse transit time to blood pressure of the subject.

20. The method of claim 19 wherein the mathematical model is a parametric model that assumes a specific structure for the dynamic system.

21. A method for improving a cardiac output measure of a subject, comprising:

measuring a waveform indicative of the arterial pulse at a proximal site of the subject using a first sensor;

measuring a waveform indicative of the arterial pulse at a distal site of the subject using a second sensor;

defining a dynamic system that relates the proximal arterial waveform to the distal arterial waveform in terms of a mathematical model, where one parameter of the mathematical model accounts for wave reflection and another parameter of the mathematical model accounts for pulse transit time;

determining the mathematical model from the measured waveforms;

determining a pulse transit time between the proximal site and the distal site from the mathematical model where the steps of determining the mathematical model and determining a pulse transit time are implemented by a computer processor of a computing device; and

improving a cardiac output measure for the subject using the pulse transit time.

22. The method of claim 21 wherein the mathematical model is a parametric model that assumes a specific structure for the dynamic system.

23. The method of claim 21 further comprises improving a cardiac output measure from a blood pressure waveform by scaling it by an exponential of the pulse transit time so as to correct it for arterial compliance changes.

24. The method of claim 21 further comprises improving a cardiac output measure by calibrating it to absolute cardiac output using the pulse transit time and anthropomorphic information of the subject.