System and method for characterizing circulatory blood flow

• US 9,002,440 B2
• Filed: 03/15/2013
• Issued: 04/07/2015
• Est. Priority Date: 07/08/2010
• Status: Active Grant

• Current Assignee: Intelomed, Inc.
• Original Assignee: Intelomed, Inc.
• Inventors: Berkow, Jan K., Brumfield, Anne M.
• Primary Examiner(s): Kahelin, Michael
• Assistant Examiner(s): Alter, Mitchell E
• Application Number:
  US13/839,534
• Publication Number:
  US 20130267858A1
• Time in Patent Office: 753 Days
• Field of Search: 600/526
• US Class Current: 600/479
• CPC Class Codes: A61B 5/0002 Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network A61B5/07 takes precedence; transmission systems for measured values G08C; transmission H04B; transmission of digital information H04L; wireless communication networks H04W
• CPC Class Codes: A61B 5/02028 Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
• CPC Class Codes: A61B 5/02116 of pulse wave amplitude A61B5/02225 takes precedence
• CPC Class Codes: A61B 5/02416 using photoplethysmograph signals, e.g. generated by infrared radiation A61B5/14552 takes precedence
• CPC Class Codes: A61B 5/0261 using optical means, e.g. infrared light
• CPC Class Codes: A61B 5/1455 using optical sensors, e.g. spectral photometrical oximeters
• CPC Class Codes: A61B 5/7225 Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation A61B5/30, A61B5/316 take precedence; transformers for specific applications H01F38/00; amplifiers H03F
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: A61B 5/726 using Wavelet transforms
• CPC Class Codes: A61B 5/7278 Artificial waveform generation or derivation, e.g. synthesising signals from measured signals A61B5/327 takes precedence

Abstract

A computer-implemented method for characterizing circulatory blood volume and autoregulatory compensatory mechanisms to maintain circulatory blood volume is disclosed. A biological signal that emulates the arterial pulse wave is collected from a sensor. Three derived parameters are extrapolated from the biological signal. The first parameter, circulatory stress, reflects of the changes of the heart rate frequency. The second, circulatory blood volume, reflects the changes in the frequency strength of the heart rate frequency. The third, Pulse Volume Alteration (PVA) Index is a ratio of the sum of the strengths of the heart rate frequency harmonics to the strength of the heart rate frequency of the unprocessed biological signal. Each parameter is compared to a threshold value and assessed to determine an adequacy of circulatory blood volume and an appropriateness of the autoregulatory mechanisms used to maintain circulatory blood volume adequacy.

Claims

1. A method for characterizing a circulating blood volume, the method comprising:

receiving, by a computing device, a biological signal emulating an arterial pulse wave from a sensor associated with a human body;

calculating, by the computing device, a plurality of integer harmonics of a heart rate from the biological signal, wherein each of the plurality of integer harmonics is characterized by an integer harmonic amplitude;

calculating from the plurality of integer harmonic amplitudes of the heart rate, by the computing device, at least one derived parameter in a frequency domain comprising a measure of total harmonic distortion; and

characterizing the circulating blood volume by comparing, by the computing device, the at least one derived parameter to a baseline value.

2. The method of claim 1, wherein the sensor comprises one or more of the following:

a transmissive photo-optic sensor;

a reflective photo-optic sensor;

a pressure transducer;

a tonometry device;

a strain gauge;

an ultrasound device;

an electrical impedance device; and

a radar device.

3. The method of claim 1, further comprising:

conditioning the biological signal to form a conditioned biological signal; and

calculating, by the computing device, at least one second derived parameter in the frequency domain based on the conditioned biological signal.

4. The method of claim 3, wherein conditioning the biological signal comprises amplifying the biological signal or filtering the biological signal.

5. The method of claim 1, wherein calculating a plurality of integer harmonics of a heart rate from the biological signal comprises:

sampling, by the computing device, the biological signal into a plurality of discrete data within a time window, to form windowed discrete signal data; and

performing, by the computing device, a spectrum analysis of the windowed discrete signal data.

6. The method of claim 5, wherein performing a spectrum analysis of the windowed signal data comprises applying a Fast-Fourier Transform to the windowed signal data.

7. The method of claim 5, wherein performing a spectrum analysis of the windowed signal data comprises applying a wavelet transformation to the windowed signal data.

8. The method of claim 1, further comprising comparing a pattern of the at least one derived parameter in the frequency domain over time to a library of patterns of the derived parameter in the frequency domain over time.

9. The method of claim 1, wherein calculating the measure of the total harmonic distortion comprises summing together the integer harmonic amplitude of each of the plurality of integer harmonics of the heart rate.

10. The method of claim 1, wherein calculating the measure of the total harmonic distortion comprises:

summing together the integer harmonic amplitude of each of the plurality of integer harmonics of the heart rate to yield a numerator, and

dividing the numerator by an amplitude of a fundamental frequency of the heart rate.

11. The method of claim 1, wherein calculating the measure of the total harmonic distortion comprises:

setting a denominator equal to the integer harmonic amplitude of an integer harmonic having a largest integer harmonic amplitude value among the plurality of integer harmonics of the heart rate;

taking a square root of a sum of squares of the integer harmonic amplitudes of the plurality of integer harmonics of the heart rate except for the integer harmonic of the heart rate having the largest integer harmonic amplitude value among the plurality of integer harmonics of the heart rate to yield a numerator, and

dividing the numerator by the denominator.

12. A method for characterizing a circulating blood volume in an animal body in a response to at least one stress, the method comprising:

characterizing a first circulating blood volume by;

receiving, by a computing device, a first biological signal emulating an arterial pulse wave from a sensor associated with the animal body,

calculating, by the computing device, a first plurality of integer harmonics of a heart rate from the first biological signal wherein each of the first plurality of integer harmonics is characterized by a first integer harmonic amplitude,

calculating from the plurality of first integer harmonic amplitudes of the heart rate, by the computing device, at least one first derived parameter in a frequency domain comprising a first measure of total harmonic distortion, and

characterizing the first circulating blood volume by comparing, by the computing device, the at least one first derived parameter to a first baseline value;

applying at least one stress to the animal body;

determining a measure of a second circulating blood volume by;

receiving, by the computing device, a second biological signal emulating an arterial pulse wave from a sensor associated with the animal body,

calculating, by the computing device, a second plurality of integer harmonics of a heart rate from the second biological signal, wherein each of the second plurality of integer harmonics is characterized by a second integer harmonic amplitude,

calculating from the plurality of second integer harmonic amplitudes of the heart rate, by the computing device, at least one second derived parameter in a frequency domain comprising a second measure of total harmonic distortion, and

characterizing the second circulating blood volume by comparing, by the computing device, the at least one second derived parameter to a second baseline value; and

comparing, by the computing device, the characterization of the first circulating blood volume to the characterization of the second circulating blood volume to characterize the circulating blood volume as a response to the at least one stress.

13. The method of claim 12, wherein the animal is a human.

14. The method of claim 12, wherein applying at least one stress comprises reducing the animal'"'"'s circulatory volume during dialysis.

15. The method of claim 12, wherein applying at least one stress comprises reducing a blood volume of the animal.