Cardiovascular Measurements by Sound Velocity Dilution

The present invention relates, in general, to the measurement of hemodynamic parameters by indicator dilution, and more particularly to measurements of blood flow, blood volume, cardiac output, protein concentration, lung water, and solution clearance through the measurement of the sound velocity dilution effects of an indicator in the bloodstream of a patient.

In a large number of medical procedures, it is very important to monitor with precision hemodynamic parameters in a patient. For example, in procedures such as hemodialysis, a filter is used to remove selected particles and liquids from the patient's bloodstream, but if this material is removed too quickly, the vascular system of the patient may collapse.

Medical parameters such as protein concentration in the blood, cardiac output, lung water content, blood circulating volume, and the like also require careful measurement and monitoring, and improved accuracy in such measurements is highly desirable.

Difficulties in making such measurements have resulted from the fact that such procedures usually involve extracorporeal circulation of the blood from a patient through, for example a blood treatment system, and in many cases the effects of the system itself on the blood flow or on the measurement devices is unknown. For example, if blood is directed to a dialysis filter through plastic tubing, the effect of the plastic material on measuring equipment using ultrasound waves may not be known with any certainty, since the characteristics of the material can vary from one tube to another. Difficulties are also encountered in perivascular measurements, i.e., measurements made using sensors mounted on a blood vessel for the thickness of a vessel wall can adversely affect the measurement. Even when measurements are being made in the patient's body, as by an ultrasonic blood flow meter clamped around a blood vessel, the thickness and morphology of the vein/artery will affect the measurement.

It is, therefore, an object of the present invention to provide a method and apparatus for accurately and reliably measuring blood parameters in situ within a patient's vein or artery, externally of the vein or artery, or in tubing leading to a blood treatment system which carries the blood exteriorly of the patient's body.

It is another object of the invention to accurately and reliably measure a patient's blood to determine blood parameters such as blood volume changes, protein concentration or cardiac output, to determine lung water, and to determine solution clearance of a blood testing or monitoring system, through measurement, over a period of time, of blood dilution.

It is another object of the present invention to measure parameters of a patient's blood by measuring sound velocity dilution in the blood resulting from the injection of known quantities of a diluting indicator material.

Briefly, the present invention relates to a method and apparatus for measuring blood parameters by dilution. More specifically, the invention relates to a method and apparatus for measuring various blood parameters utilizing a dilution process wherein an indicator material such as a saline solution is injected into a bloodstream and its downstream effects on measured quantities such as the velocity of ultrasound in the bloodstream and electrical impedance of the blood provide information about a variety of blood flow parameters.

It is known that the velocity of ultrasound in blood is a function of, among other things, the proteins and ions contained in the blood, with the sound velocity increasing with an increase in protein concentration. Accordingly, the velocity of sound through a blood sample can be varied by diluting the blood with an indicator material having different acoustical characteristics than those of the blood; for example, through the use of a saline solution that has no proteins. By injecting a known volume of such an indicator into a known blood flow, the diluting effect of the indicator over a period of time can be accurately determined by a sensor responsive to changes of sound velocity in the blood. The sensor is positioned downstream of the injection so that the indicator passes the sensor, with the measured diluting effect being used to determine various blood parameters.

The measurements can be made in situ in a patient's vessel through the use of a catheter to provide intravascular measurements, or can be done through the use of measuring sensors clamped onto a patient's blood vessel, to provide perivascular measurements. In another alternative, the measurements may be made in an extracorporeal blood system in which clamp-on measuring sensors are secured for example, to tubing leading to exterior blood treatment equipment such as a hemodialysis machine, or the like. In general, then, in accordance with the present invention, ultrasound measurements are made of the characteristics of a patient's blood. Thereafter, a bolus, or known volume, of an indicator material having different characteristics than those of the blood is injected into the bloodstream and measurements are made of changes in sound velocity to determine the passage of the bolus past the sensor. The changes in such characteristics can then be plotted and used to determine parameters such as blood flow, blood volume, cardiac output, protein concentration, lung water, solution clearance, blood recirculation within a treatment system, and the like. Dilution measurements are principally by sound velocity, with increased accuracy being provided by measuring both sound velocity and amplitude. Additional sensors can be used to measure various other blood characteristics such as temperature or electrical impedance to increase the accuracy of the determinations made by the present invention.

More particularly, the blood parameter measurements in accordance with one embodiment of the present invention can be carried out by injecting a known volume of an indicator isotonic saline solution at blood temperature through a catheter into an artery, for example, or into an arterial tube in a blood treatment system. This solution, which is carried downstream by the arterial blood flow, dilutes the blood protein concentration, thereby producing a change in the measured velocity of sound in the blood when it passes through an arterial sensor located downstream from the injection point. This change in sound velocity permits a rapid and accurate intravascular measurement of, for example, cardiac output (volume per unit time), as well as measurement of blood protein concentration.

In cases where the sensor does not make direct contact with the blood, as in perivascular or clamp-on applications, the ultrasonic wave of the sensor passes through a vessel or tube as well as through the blood. The measured sound velocity then is a function of the material and the geometry of the vessel or tube and their acoustical properties, as well as of the acoustical properties of the blood. In such cases, the effects of the vessel or tube must be determined and then canceled out of later measurements so that the results of later blood parameter measurements will be independent of the sonic and geometric properties of the tube or the artery/vein.

Sound velocity is also a function of blood temperature, and accordingly in another embodiment of the invention the indicator may be a volume of a liquid having a different temperature than that of the bloodstream. In this case, the indicator material can be detected not only by a sound velocity sensor, but by a temperature sensor, as well, to provide improved accuracy. In similar manner, indicators having a different optical density or electrical impedance than that of the bloodstream may be used, with appropriate sensors, to detect changes in blood characteristics.

The intravascular sensors described above can be located in either an artery or a vein of a patient, or both. When in an artery, including the pulmonary artery, measurements of cardiac output (CO) may be made, but when a vein, the measurement is of the flow through the patient's organs. In the latter case, the measurements can be used to determine solution clearance; that is, the ability of an organ to respond to, or to recover, different solutes. Thus, for example, a special solution can be injected into an artery upstream of an organ to be tested and measurements made in a vein downstream of the organ to measure the organ's reaction to that solution and from this to determine specific characteristics of the organ. For example, a patient's liver can be tested by injecting a solution containing urea.

Perivascular measurements may also be made on either an artery or a vein, or both, but these measurements are made externally of the vessel, as noted above. Sensors for this purpose may be ultrasonic flow meters such as those manufactured by Transonic Systems, Inc., Ithaca, N.Y. or other suitable flow meters utilizing ultrasonic waves for measurement of flow. In perivascular systems, however, the material and geometric characteristics of the vessel itself affect the measurements, so compensation is required, as will be described hereinbelow.

In another form of the invention, referred to herein as a clamp-on measurement system, measurements of blood are made outside the patient's body; for example, in tubing leading from the patient's vascular system to a blood treatment system such as a dialysis machine. In such a system, blood is drawn from the patient, passed through suitable tubing to a dialysis filter and is then returned through tubing to the patient, usually to the same vessel from which the blood being treated was initially drawn, but downstream of the withdrawal site. Clamp-on sensors such as Transonic Systems' ultrasonic flow meters are suitable for measuring blood flow through the tubing, but the measured parameters are functions not only of blood characteristics, but of tube material and geometry. In such a case, calibration of the system is required to eliminate the effects of the tubing. As with intravascular and perivascular systems, arterial or venous measurements, or both, may be made in the clamp-on system, with injections of indicator material in either the artery or the vein, or both, to permit determination of various blood characteristics as well as permitting measurement of the effectiveness of a treatment system.

Ultrasonic sensors measure sound velocity dilution as the indicator material is carried past the sensor by the bloodstream, and changes in sound velocity are plotted to permit calculation of various blood parameters. The time at which the indicator material reaches the sensor after injection, the area under the plotted curve representing the changes in sound velocity at the sensor, and the amplitude of the measurement all provide information concerning the blood characteristics, in accordance with the invention.

The foregoing, and additional objects, features, and advantages of the present invention will become apparent to those of skill in the art from the following detailed description of preferred embodiments thereof, taken in conjunction with the accompanying drawings, in which:

FIG. 1 is a diagrammatic illustration of an intravascular blood measurement system utilizing sound velocity measurements;

FIG. 2 is a diagrammatic illustration of a catheter for use in the system of FIG. 1;

FIG. 3 is a dilution curve showing the change in blood sound velocity following injection of an indicator material;

FIG. 4 is a temperature dilution curve for the system of FIG. 1;

FIG. 5 is a curve representing calculated sound velocity changes due to temperature;

FIG. 6 is a sound velocity dilution curve corrected for temperature variations;

FIGS. 7-9 are diagrammatic illustrations of sensors suitable for use with the catheter of FIG. 1;

FIG. 10 is a diagrammatic illustration of a perivascular blood measurement system utilizing sound velocity measurements;

FIG. 11 is a diagrammatic illustration of a clamp-on sensor system utilizing sound velocity measurements;

FIGS. 12-15 are diagrammatic illustrations of various blood measurement sensors which may be utilized in the systems of FIGS. 10 and 11;

FIG. 16 is a graphical illustration of a change in blood temperature over time in a tube;

FIG. 17 is a graphical illustration of a change in tube sound velocity due to the temperature change of FIG. 16;

FIG. 18 is a graphical illustration of a change in amplitude of received sound waves due to the tube temperature change of FIG. 16;

FIG. 19 is a graphical illustration of sound velocity changes produced by injections of indicator material;

FIG. 20 is a graphical illustration of changes in venous and arterial sound velocities following an intravenous injection of an isotonic solution in the system of FIG. 11 for determining cardiac output;

FIG. 21 is a graphical illustration of changes in venous and arterial sound velocities following an arterial injection of a solution for determining organ or filter clearance;

FIG. 22 is a graphical illustration of changes in sound velocity in the arterial and venous lines of a blood treatment system for use in determining blood clearance with an arterial and a venous injection;

FIG. 23 is a graphical illustration of changes in sound velocity in the arterial and venous lines of a blood treatment system for use in determining blood circulation, using arterial and venous injections;

FIG. 24 is a graphical illustration of changes in sound velocity in the arterial and venous lines of a blood treatment system for use in determining blood recirculation, using only an intravenous injection;

FIGS. 25-27 are graphical illustrations of the transfer of water and solute between a patient's lungs and bloodstream as a bolus of hypertonic solution passes through the lungs;

FIGS. 28 and 29 are graphical illustrations of sound velocity and blood impedance changes for use in calibrating the measurement of lung water;

FIGS. 30 and 31 are graphical illustrations of arterial sound velocity and blood impedance changes responsive to a calibrating intravenous injection of hypertonic saline solution for determination of lung water;

FIGS. 32-35 are graphical illustrations of arterial blood sound velocity and impedance measurements in response to an intravenous injection of isotonic saline solution for lung water volume measurement;

FIGS. 36-39 are graphical illustrations of venous or arterial blood sound velocity dilution and electrical impedance dilution measurements in response to a venous injection of hypertonic saline solution; and

FIGS. 40-42 are graphical illustrations of arterial sound velocity changes due to temperature for lung water measurement.

Turning now to a more detailed description of the present invention, there is illustrated in FIG. 1 in diagrammatic form an intravascular blood measurement system 10 utilizing sound velocity measurements for determining blood characteristics. The ultrasound velocity in blood, C.sub.b, can be described as a function of total protein concentration P, and concentration of ions I, (D. Schneditz and T. Kenher, J.Acoust.Soc.Am. 86 (6), pp. 2073-2081, (1989) "A sound-speed sensor for the measurement of total protein concentration in disposable, blood-perfused tubes") and by V. A. Del Grosso and C. W. Mader, J.Acoust.Soc.Am. 52, 961-974 (1972) "Speed of sound in sea water samples"), and may be expressed as follows:

where A.sub.0, A.sub.1, A.sub.2, and A.sub.3 are known coefficients that are functions of temperature; N is the amount of protein in blood volume V.sub.b ; and I is the average concentration of ions in plasma and erythrocytes.

The cardiovascular system of a patient is indicated in FIG. 1 by a vessel 12 which may be an artery leading to a patient's organs or to a selected organ indicated by block 14. A vessel 16 is a vein carrying blood away from the organ 14 to other parts of the cardiovascular system indicated by block 18. Additional arteries or veins indicated in phantom at 20 may also lead to the organ 14. An arterial catheter 22 is inserted into the artery 12, in accordance with one embodiment of the present invention, to provide a passageway for injecting an indicator material such as a saline solution into the artery at a selected location. The catheter thus includes an inlet port 24 accessible outside the patient's body and an outlet port 26 located on the catheter for positioning within a vessel for this purpose. The catheter also includes a plurality of sensors, illustrated in greater detail in FIG. 2, to which reference is now made.

As illustrated, the catheter preferably incorporates a conventional temperature sensor 30 for measuring blood temperature, and a sound velocity sensor 32 for determining the velocity of ultrasonic waves in the blood flowing past the catheter. These two sensors are located downstream of the catheter outlet 26 so that indicator material injected into the bloodstream will pass by the sensors. The direction of blood flow in the artery 12 and past the catheter 22 is indicated by arrows 34 in FIGS. 1 and 2. It is noted that a balloon 36 may be provided on the catheter for securing it at a desired location, the expansion of the balloon being controlled by a source 38 of fluid under controlled pressure connected to the balloon internally of the catheter as by way of line 40. The temperature sensor 30 is connected by way of lead 42 to a suitable temperature detecting circuit 44 while the sound velocity sensor 32 is connected by way of lead 46 to a suitable sound velocity detecting circuit 48. In addition to these sensors, a blood pressure sensor 50 may also be provided on the catheter 22 and connected by way of lead 52 to a blood pressure detector 54. If desired, an additional measurement of sound amplitude may be made by means of sound amplitude detector 56 also connected by way of lead 46 to sensor 32.

As illustrated in FIG. 1, each of the detectors 44, 48, 54, and 56 is connected by a corresponding one of leads 60, 62, 64,and 66 to an evaluating device 70, which may be a personal computer for carrying out the curve plotting and calculations described herein.

The catheter with its sound velocity sensor 32, when placed in the blood vessel 12, provides a direct measure of the velocity of ultrasound in the blood carried by the artery, and permits a determination of blood parameters such as the total protein concentration in the blood and its changes based on sound velocity changes resulting from the injection of an indicator material and the resultant dilution of the blood. The dilution of flowing blood by the injection of an isotonic saline solution, for example, injected at body temperature, decreases the blood sound velocity in accordance with the following equation:

where P is the total protein concentration in the blood before the injection. The change in protein concentration with time .DELTA.P(t) can be calculated from equation 2 and gives the following result: ##EQU1## where V.sub.b is the volume of blood that mixes with the indicator and V.sub.isot. (t) is the amount of indicator at time (t) in the blood. The resulting change in blood sound velocity can then be rewritten as follows: ##EQU2## Summation over a period of time such as one minute will give the average flow through the system during this time, as follows: ##EQU3## where S.sub.b.isot.sound is the area under a dilution curve representing changes in blood sound velocity over time, Q is the average flow through the system during the given period of time, and V.sub.isot. is the volume of the injected indicator solution. FIG. 3 illustrates such a dilution curve 72, showing changes in sound velocity due to an injection of an indicator medium at time t.sub.1 from injection outlet 26 of catheter 22. The change in blood sound velocity is measured at sensor 32 (FIG. 2) downstream of the injection site. The cross-hatched area 74 under the curve 72 is determined in evaluating device 70 to permit calculation of the average flow through the system in accordance with Eq. 6.

For the sound velocity sensor located in an artery within the blood circulating system of a patient; for example, in a pulmonary artery, and for an intravenous injection V.sub.isot. at time t.sub.1, the measured value Q from equation 6 is the cardiac output CO, so that Q=CO. For the intraarterial injection, this measured flow is also the blood flow through the patient's organ or organs. The thermal sensor 30 on the catheter increases the accuracy of measurement because the velocity of sound in blood varies with temperature. By directly measuring blood temperature, therefore, more accurate calculations of the "A" coefficients in equation 1 for protein and for CO calculations in equation 6 are obtained.

If the temperature of the injected indicator medium is not equal to body temperature; for example is colder, than measurements to provide the curves of FIGS. 4, 5, and 6 may be made to correct the recorded sound velocity to obtain highly accurate blood parameter calculations. Thus, variations in blood temperature at sensor 30 following the intravenous injection at port 26 are measured and plotted, as illustrated by curve 76 in FIG. 4, to provide a temperature dilution curve .DELTA.T(t). This dilution curve 76 is used to calculate sound velocity changes due to the influence of temperature changes, the calculations being indicated by curve 78 in FIG. 5. This curve provides a calculated temperature correction C.sub.T (t)=A.DELTA.T(t), where A is known, which is used to correct a recorded sound velocity dilution curve 80 (FIG. 6), indicated by C.sub.b (t), to obtain a highly accurate cardiac output value indicated by curve 82 and providing the value C.sub.b.correct (t).

Because an isotonic solution leaves an intravascular space slowly, as indicated by curve 72 in FIG. 3, this curve can also be used to calculate blood circulating volume (BCV), which can be calculated from FIG. 3 as follows:

where .DELTA.C(t.sub.mix) is the decrease of sound velocity produced when there is a full mixing of the indicator material with the blood in the artery 12.

As further illustrated in FIG. 1, a second catheter 22' may also be incorporated in the blood measurement system 10, this catheter being inserted in vein 16 for use in measuring blood flow therein. The catheter 22' is substantially identical to catheter 22 with similar elements being similarly numbered, but primed. The measurements available with the use of catheter 22' will be described hereinbelow.

Sensors usable with the catheter 22 are illustrated diagrammatically in FIGS. 7, 8, and 9. FIG. 7 illustrates a sound velocity sensor 32 having an ultrasonic sound wave generating crystal 90 and a similar receiving crystal 92 located on opposite ends of a cavity 94 formed in the surface of the catheter. Sound waves emitted by crystal 90 are received by crystal 92, the sound waves traveling through the blood surrounding the catheter and filling the cavity 94. The blood flow carries the indicator medium, as described above, which changes the velocity of the sound waves crossing from one end to the other of cavity 94. The output from sensor 32 is supplied to the sound velocity detector 48 which measures the time required for the waves to travel from transmitter 90 to receiver 92, thereby providing the data for plotting the curve of FIG. 3.

FIG. 8 similarly illustrates an electrical impedance sensor 96 located on the surface of catheter 22, the sensor incorporating a plurality of electrically conductive strips 98, 100 which measure the conductivity of the blood 34 passing the sensor. This conductivity measurement allows detection of the impedance of the blood by a corresponding detector (not shown). FIG. 9 illustrates an optical sensor 102 located in a cavity 104 on a surface of catheter 22. The sensor includes a light source 106 at one end of the cavity and a light sensor 108 at the other end of the cavity for detecting the optical density of the blood 34 passing the catheter 22. A suitable detector connected to sensor 108 may be provided to produce measurements of the optical characteristics of the blood. These and other known sensors may be used with catheter 22 to provide the required blood characteristics.

A modification of the system of FIG. 1 is illustrated in FIG. 10, wherein the system utilizes perivascular measurements instead of the intravascular measurements obtained by the catheter of FIG. 1. Elements of FIG. 10 common to FIG. 1 carry similar identifying numerals. In the perivascular system of FIG. 10, clamp-on sensors 110, 112, 114, and 116 are provided on artery 12, with corresponding sensors 110', 112', 114', and 116' being provided on the exterior of vein 16. Sensors 110 and 110' are blood flow measurement sensors which provide output signals by way of lines 118 and 118' to corresponding blood flow measurement detectors 120 and 120'. Each of the blood flow sensors may be a Bypass Flow Meter model HT 109 or model T106 produced by Transonic Systems, Inc., Ithaca, N.Y., for example. Sensors 112 and 112' measure sound velocity in the bloodstream of the artery 12 and the vein 16, respectively, and produce outputs on their corresponding lines 122 and 122' to corresponding sound velocity measurement detectors 124 and 124'. The sensors 112 and 112' each may be a sensor such as model No. DT 106, also available from Transonic Systems, Inc.

Sensors 114 and 114' are blood dilution sensors which are connected by way of lines 126 and 126' to dilution measurement detectors 128 and 128. Similarly, sensors 116 and 116' are additional blood dilution sensors and produce output signals on lines 130 and 130' to corresponding dilution detectors 132 and 132'. The outputs of detectors 120, 124, 128, and 132 are connected by way of lines 134, 136, 138, and 140, respectively, to the evaluating device 70, and in similar manner detectors 120', 124', 128', and 132' are also connected to the evaluating device. The dilution sensors 114, 114' and 116, 116' may be any of a variety of sensors for detecting various blood characteristics. For example, sensors 114 and 114' may be electrical impedance sensors, and sensors 116, 116' may be optical sensors, the particular sensors selected being dependent on the blood characteristics of interest.

A still further modification of the invention is illustrated in FIG. 11, wherein the exterior, or clamp-on sensors utilized in the system of FIG. 10 are used to measure blood flow through a conduit, or tubing, 150 leading from a vessel 152 of a patient to external, i.e., extracoporeal, equipment such as a hemodialysis filter 154. The vessel 152 may be an artery of a patient or a shunt placed in the cardiovascular system of a patient. As indicated by arrows 156, a portion of the blood flowing through the vessel 152 is drawn through inlet, or arterial tubing 150 for delivery to the hemodialysis filter 154. The blood which passes through the filter is then returned to vessel 152 by way of outlet tubing 158 which may be referred to as the venous tubing. Surrounding, or clamped onto the arterial tubing 150 are a plurality of blood sensors such as a blood flow sensor 160 connected by way of line 162 a blood flow detector 164. Sensor 160 is similar to the blood flow sensor 110 of FIG. 10 and provides an output on line 162 which is a measure of the volume of blood flow flowing through tubing 150. The signal from detector 164 is supplied by way of line 166 to the evaluating device 70 which carries out the determinations of blood parameters in accordance with the relationships described herein.

Also clamped around arterial tubing 150 is a temperature sensor 170 which provides an output by way of line 172 to a detector 174, the output of which is supplied by way of line 176 to the evaluating device 70. The temperature sensor measures the temperature of the blood flowing through tubing 150. Alternatively, the sensor 170 may be a sound amplitude sensor and detector 174 may provide a measure of variations of sound amplitude changes.

A third sensor 180 is clamped around tube 150 for dilution measurement. Although this sensor may be any one of the dilution sensors discussed above, a sound velocity sensor for changes in sound velocity in the bloodstream is preferred. The output from sensor 180 is supplied by way of line 182 to a sound velocity detector 184, the output of which is supplied by way of line 186 to the evaluating device 70.

A fourth sensor 190 is also provide on arterial tubing 150. This sensor is another dilution sensor, and may be an optical sensor, an impedance sensor, or the like for measuring additional blood parameters, with the output from the sensor being supplied by way of line 192 to a corresponding detector 194, the output of which is supplied by way of line 196 to the evaluating device 70.

The venous tubing 158 may also carry similar sensors indicated at 160', 170', 180', and 190', the outputs of which are supplied to corresponding detectors 164', 174', 184' and 194', with the detected signals being supplied to the evaluating device 70.

FIGS. 12-15 illustrate various clamp-on sensors which may be utilized in the systems of FIGS. 10 and 11. Thus, FIG. 12 illustrates an ultrasonic sensor, which may be used as sensor 180 or 190, for example, having an ultrasonic wave generating crystal 200 positioned on one side of tubing 150 and a corresponding receiving crystal 202 located on the opposite side of the tubing. Sound waves 204 generated by crystal 200 travel through the wall of tubing 150, through the blood 156 flowing within the central passageway 206 of the tubing, back through the wall of the tubing and are received by crystal 202. As is known, the received ultrasound varies with changes in the acoustical characteristics of the blood flow.

FIG. 13 illustrates a second type of ultrasonic sensor, also suitable for use as the sensor 180 or 190 in FIG. 11, for example, but differing in that the transmitting crystal 200 and the receiving crystal 202 are located on the same side of the tube 150. In this device, the wave 204 emitted by crystal 200 is transmitted through the tubing and the bloodstream 156, and is reflected by a mirror 208 back through the tubing and the bloodstream to the receiving crystal 202.

FIG. 14 illustrates an electrical impedance sensor which also may be used as the sensor 180 or 190 in FIG. 11. This sensor includes a pair of electrodes 210 and 212 located, for example, on opposite sides of tubing 150 for measuring the electrical conductivity of the bloodstream 156. Still another sensor suitable for use in the system of FIG. 11 is a clamp-on optical sensor such as that illustrated in FIG. 15, which may include a light source 214 on one side of the tubing 150 and a light sensor 216 on the opposite side for measuring the optical density of the bloodstream.

In the embodiments of both FIG. 10 and FIG. 11, the ultrasonic waves of the sound sensors, such as sensors 110 and 112 in FIG. 10, and sensors 160 and 180 in FIG. 11, pass through the conduit, which is either the vessels 12, 16, or the tubing 150,158 as well as the blood. In such cases, the measured sound velocity is a function of the geometry of the tube/vessel and its acoustical properties, as well as the acoustical properties of the blood. The relationship between the measured sound velocity C and an unknown blood sound velocity C.sub.b can be stated as:

where C.sub.t is the average velocity of ultrasound through the material of the tube/vessel, where l.sub.t is the equivalent path length of ultrasound through the tube/vessel, and where l.sub.b is the equivalent path length through the blood.

In order to take into account the effects of the tube/vessel acoustical properties, it is necessary to calibrate the systems of FIGS. 10 and 11. For the system of FIG. 11, the tube system can be calibrated by filling it with a solution of known sound properties, such as a saline solution, and measuring the sound velocity. This measurement can then be compared with the velocity of sound in the blood, and the relationship between C and C.sub.b can be found. This relationship will be valid for a constant temperature.

If, during dialysis for example, the blood temperature changes, it will change the sound velocity not only in the blood, but in the tube, as well. This will decrease the accuracy of the measurements because of the strong influence of tube temperature on sound velocity, the velocity changing in the range of 3 to 7 m/s/.degree. C. for different materials. Because of this, the real change in sound velocity due, for example, to a variation in protein concentration .DELTA.C(t) will be the difference between the measured changes in ultrasound velocity .DELTA.C.sub.m (t) and the change due to the influence of temperature .DELTA.C.sub.T (t), which may be expressed as follows:

Temperature changes in the bloodstream, such as the changes illustrated in FIG. 16 by curve 220, may occur during treatment. As noted above, a temperature change will also influence both the ultrasound velocity in the tube itself, and the attenuation of sound in the tube itself, as illustrated by curve 222 in FIG. 17 and curve 224 in FIG. 18. Thus, for example, an increase in blood temperature T.sub.b (curve 220 in FIG. 16) heats the tube, making it softer and changing the conditions for sound wave propagation through the tube material. As illustrated by curve 222 in FIG. 17 the ultrasound velocity drops with increased temperature and, as illustrated by curve 224 in FIG. 18, attenuation of the signal increases, causing the amplitude A of the received signal to drop. Experiments have shown that because these changes arise from the same cause, they behave synchronously and may be expressed as C.sub.T (t), the change in sound velocity with temperature. The recorded changes in amplitude of the received acoustical signals are only due to tube temperature changes, because blood does not attenuate sound waves up to about 10 to 15 MHz; accordingly, up to these frequencies, changes in blood parameters such as protein content don't influence sound attenuation. This makes it possible to compensate for the influence of temperature on the measured ultrasound velocity by recording the changes in the amplitude of the received sound signal .DELTA.A(t) due to changes in temperature in accordance with the following equation:

where K.sub.T is a coefficient which can be calculated for the type of tubing that is used in the blood treatment system. Such calculations may be made prior to use of the tubing or may be calculated during use in, for example, dialysis as a ratio of the slope of the sound velocity and the slope of the sound amplitude at the beginning of the dialysis procedure when the tube is being heated or cooled by the blood, as expressed in the following equation:

It is noted that the perivascular system of FIG. 10 cannot be precalibrated because the vessel wall size and its properties are unknown. Accordingly, the perivascular probe must be calibrated during dialysis or other treatment.

Another way to calibrate the perivascular sensors or the clamp-on tube sensors of FIGS. 10 and 11 is to make a calibration injection of a known indicator material directly prior to the location of the sound velocity sensor while simultaneously measuring the blood flow through the tube/vessel. A calibration injection into the arterial line, for example by way of arterial inlet ports 230 or 232 in the systems of FIGS. 10 and 11, respectively, will change the measured sound velocity in the arterial sensors 112 or 180 as follows:

or

where S.sub.A.cal.isot.sound. is the measured sound velocity dilution area generated by blood sound velocity dilution changes (S.sub.A.b.cal.isot.sound). From equations 6 and 13 the average flow through the tube/vessel during the time when the indicator passes is expressed as follows: ##EQU4## where Q.sub.A.cal. is the average flow through the tube/vessel and V.sub.A.cal.isot. is the volume of the calibrating injection.

The sensors carried by the vein 16 or the venous tube 158 in FIGS. 10 and 11, respectively, may also be used to calculate various blood flow parameters. Thus, for example, by analogy with equation 6, the cardiac output CO can be calculated after an injection at body temperature of an isotonic saline indicator solution into the venous port (by way of catheter 234 in FIG. 10 or 236 in FIG. 11) and detecting the indicator at sensors 112' or 180'. Thus, from equations 6, 13, and 14, the cardiac output is calculated as follows: ##EQU5## where S.sub.A.isot.sound. is the measured sound velocity dilution area under a curve 240 generated by the measurement of dilution produced by an intravenous injection of a volume V.sub.V.isot. at t.sub.1 illustrated in FIG. 19.

Prior to the foregoing measurement, a measured sound velocity dilution area S.sub.A.cal.isot.sound. is produced by curve 241 in FIG. 19 and, as illustrated, follows closely after a calibration injection upstream of the sensor.

The cardiac output can be calculated from a single intravenous injection at port 236, when the dilution sensors are located on both the arterial and venous parts of the system and are closely matched, as illustrated in FIG. 11.

The cardiac output can be calculated from a single intravenous injection (Q.sub.A.cal.isot. .fwdarw.Q.sub.V.isot. and S.sub.A.cal.isot. .fwdarw.S.sub.V.isot.) by obtaining the measurements illustrated in FIG. 20. As there illustrated, an intravenous injection is made at a time t.sub.1 indicated by vertical line 250, and shortly thereafter sound velocity sensor 180' detects in the venous line the dilution illustrated by curve 252. A portion of the indicator material is then recirculated to the arterial line 150 where it is detected by sensor 180, as indicated at curve 254. The cardiac output CO can then be calculated as follows: ##EQU6## where Q.sub.V.isot. is the blood flow in the venous line during the isotonic injection and S.sub.V.isot. is the dilution area generated in the venous sensor. These calculations of cardiac output are independent of the properties of the tubing/vessel and the blood.

By analogy with equation 5, relatively small changes of circulating blood volume can be measured in the systems illustrated in FIGS. 1, 10, and 11. In the system of FIG. 1, blood volume can be measured by intravascular sound sensors on catheter 22 as follows: ##EQU7## where .DELTA.V.sub.b is the change in blood volume. In similar manner, the change of blood volume can be calculated for the configurations of FIGS. 10 and 11 from equations 12, 14, and 17, as follows: ##EQU8##

The foregoing illustrates that the monitoring of changes in ultrasound velocity .DELTA.C(t) due to indicator dilution permits the monitoring of total protein changes (equation 3) and relative changes of blood volume not only with intravascular catheters (equation 17), but noninvasively, as during a dialysis procedure (equation 18). In each case, this is done with a single calibration injection and a single measurement injection. This is important in dialysis procedures to prevent the unpleasant effects of fluid unbalancing in such a procedure.

The systems of FIGS. 1, 10, and 11 may also be used to determine the clearance K.sub.clear of different solutes passing through a patient's organ or organs 14 or through blood treatment equipment such as the hemodialysis filter 154. For this purpose, a suitable solution is injected into the appropriate arterial input port 24, 230, or 232, and the dilution signals and flows are recorded in both the arterial and venous lines, as indicated by curves 260 and 262 in FIG. 21. The changes that take place in the solution as it passes through the hemofilter/organ is obtained by the following equation:

where Q.sub.A and Q.sub.V are the corresponding blood flows in the arterial and venous lines averaged over the time the bolus of the indicator material (or medium) is passing through the sensors and where S.sub.A and S.sub.V are the corresponding dilution areas before and after the hemofilter/organ (FIG. 21).

In the case where there is only a single sound velocity sensor located on the venous line in any of the systems of FIGS. 1, 10, and 11, the clearance K.sub.clear can be calculated by two injections of an indicator material, one before the filter/organ; for example, by way of inlet port 232 in FIG. 11, and the other after the filter/organ; for example, by way of inlet port 236 in FIG. 11. The resulting curves 264 and 266 are illustrated in FIG. 22, where the intraarterial injection at inlet port 232 is made at the time indicated by dotted line 268 to produce the dilution area S.sub.AV measured at the venous sound velocity sensor 180'. FIG. 22 also illustrates the later intravenous injection made at inlet port 236 at a time indicated by dotted line 270, which produces the dilution curve 266 and the resulting area S.sub.vv. The intravenous injection is made at inlet port 236 and is detected at sound velocity sensor 180'. The value K.sub.clear can be calculated as follows:

where S.sub.AV and S.sub.VV are the corresponding dilution areas from the injections before and after the filter/organ.

As noted above, some of the blood returned to the vascular system 152 from the hemodialysis filter in the system of FIG. 11 is recirculated back to the arterial line 150, as indicated at 246. The recirculation coefficient K.sub.rec. for such a system can be calculated using intravenous and intraarterial injections of isotonic indicator material at the temperature of the blood flowing in the system, and can be calculated as follows: ##EQU9## where S.sub.A.rec. is the recirculation dilution area under curve 272, and where Q.sub.A.rec. is the blood flow in arterial line 150 during the recirculation, as illustrated in FIG. 23.

When two sensors are used, one located in the arterial and one in the venous tubing lines; for example, sensors 180 and 180' in FIG. 11, only one intravenous injection at port 236 (for FIG. 11) is needed, allowing the recirculation coefficient to be determined as follows: ##EQU10## where S.sub.A. and S.sub.V. represent the dilution areas beneath the curves obtained from the measurements in the arterial line 150 by sensor 180 and within the venous line 158 by sensor 180. These are indicated in FIG. 24 by curves 274 and 276, respectively. It is noted that when the systems of FIGS. 1 and 10 are used, the dilution curves are recorded both before and after organ deconvolution analyses are done.

The system of the present invention may also be used to measure extravascular lung water. The osmotically active part of extravascular lung water V.sub.1 can be measured from the water volume and the hypertonic solutes moving between the lungs and the vascular system when a bolus of hypertonic indicator solution passes through the vasculature of the lung (N. M. Krivitski, V. V. Kislukhin, "Determination of the Volume of Extravascular Fluid in the Lungs by Blood Electrical Resistance Using a Dilution Method", Med. Tekn. Jan-Feb; (1), pp. 6-9 (1987); N. M. Krivitski et al, "Determination of the Volume of Extravascular Fluid in the Lungs by Electroimpedance Indicators", Fisiol. Zh SSSR, Apr., 75(4), pp. 589-595 (1989). This lung water is expressed by the following equation:

where V.sub.w (t.sub.e) is the volume of water transferred to the blood from the lung until time t.sub.e (see FIGS. 25-27); t.sub.e is the time when the osmotic gradient between lung tissue and blood becomes zero (see FIG. 26); T.sub.b is the osmolality of blood before the injection of an indicator; and .DELTA.T.sub.b (t.sub.e) is the blood osmolality increase at time t.sub.e (see FIG. 25). The change of osmolality is indicated by curve 280 in FIG. 25 following the injection of an indicator at time 282. For the case where a hypertonic solution (for example NaCl), illustrated by curve 283 in FIG. 27, also moves from blood to lungs, then equation 23 will be expressed as follows:

where G.sub.NaCl (t.sub.e) is the amount of NaCl transferred to the lungs from the blood until t.sub.e, and where I is the average concentration of ion plasma and erythrocytes.

Eq. 24 is based on the following considerations. As a hypertonic bolus passes through the lung vasculature, water initially moves from the lung tissue into the blood and NaCl moves from the blood into lung tissue down its osmotic gradient, as illustrated by curve 284 in FIG. 26. As the bolus reaches the lungs at time 286, the curve 284 starts to move in the +direction toward a peak 288. This water movement increases the osmolality of the lung fluids, as indicated by curve 290 in FIG. 25 and after the peak 280 is reached curve 280 in FIG. 25) gradually decreases the osmolality of the blood. As the bolus moves out of the lung vasculature, the water transfer decreases, as indicated by curve 284 after peak 288, so that at a time t.sub.e the osmolality of the lung tissue becomes equal to that of the blood, as indicated at point 292 in FIG. 25, and net water movement ceases. After time t.sub.e the osmolality of the blood falls further, and the osmotic gradient then drives water from the blood into the lung tissue, as indicated at region 294 of curve 284, and drives NaCl from the lung tissue to the blood until osmotic equilibrium is reached at time 296.

The osmolality of blood before injection, T.sub.b, and the ion concentration I can be measured routinely. The two values of the are close, if indicated in the same units. Then, to calculate the extravascular lung water from Eq. 24 it is necessary to determine V.sub.w (t.sub.e), G.sub.NaCl (t.sub.e), and T.sub.b (t.sub.e). This can be done utilizing the apparatus of FIGS. 1, 10 and 11 when two dilution sensors are located in the arterial line. For example, in addition to the sound velocity measurement sensor 180 and its detector 184, a second dilution measurement sensor 190 is used. This may, for example, be an electrical impedance sensor with a corresponding electrical impedance detector 194. These sensors are used simultaneously with the blood flow measurement probe utilizing sensor 160 and detector 164, also in the arterial part of the system of FIG. 11.

Lung water can then be measured by injecting a volume V.sub.A.inj.cal.isot. of isotonic saline solution at body temperature at a time 298 into the arterial port 232 for calibration. The ultrasound, dilution calibration area S.sub.A.cal.isot.sound under curve 300 in FIG. 28 and the electrical impedance dilution area S.sub.A.cal.isot.imp. under curve 302 in FIG. 29 are obtained from recording the outputs of detectors 184 and 194, respectively. To obtain curves 300 and 302, an intraarterial isotonic injection of a volume V.sub.A.inj.cal.hyp. is made into the arterial port at the time indicated by dotted line 298 and the resulting changes in the arterial line of the blood sound velocity represented by curve 300 in FIG. 28 and the electrical impedance represented by curve 302 in FIG. 29 are recorded. An intravenous injection of isotonic indicator having a volume V.sub.v.inj.isot. is then made into port 236 at time 304 to produce dilution curves 306 and 308 and to corresponding dilution areas S.sub.A.isot.sound (t.sub.e) and S.sub.A.isot.imp (t.sub.e) illustrated in FIGS. 28 and 29.

Thereafter, a volume V.sub.A.inj.cal.hyp. of hypertonic saline solution having an osmolality T.sub.inj.osm. at body temperature is injected at time 309 (FIGS. 30 and 31) into the arterial port 232 for further calibration of the system. A resulting ultrasound dilution calibration area S.sub.A.cal.hyp.sound, illustrated at curve 310 in FIG. 30, and an electrical impedance solution calibration area S.sub.A.cal.hyp.imp. illustrated at 312 in FIG. 31 are then recorded. An intravenous injection of volume V.sub.v.inj.hyp. and having the same concentration as the volume V.sub.A.inj.cal.hyp. is made into the venous port 236 and the changes in the arterial line of the blood sound velocity C.sub.A.hyp. (t) represented by curve 3 14 in FIG. 30 and the impedance of the hypertonic solution Z.sub.A.hyp. (t) illustrated by curve 316 in FIG. 31 are recorded.

During the time a hypertonic solution passes through the lungs, water moves from the lungs to the blood and NaCl leaves the blood and travels to the lungs, as described above. This decreases the blood sound velocity as illustrated in FIG. 30 by the curve 314 and as illustrated by the corresponding area S.sub.A.trans.sound (t.sub.e) represented by the cross-hatched area 320. This area 320 is the difference between the total area S.sub.A.hyp.sound (t.sub.e) under the dilution curve 314 less the area S.sub.A.Or.hyp.sound (t.sub.e), illustrated by cross-hatched area 322, which is the area generated by the dilution of the blood, produced by the indicator solution, which would occur if no osmotic movement occurred, as expressed by the following equation:

By analogy, the water moving from the lungs to the blood and the NaCl leaving the blood and moving to the lungs increases the blood electrical impedance, as is illustrated in FIG. 31. In this figure, curve 324 represents the change in measured impedance upon injection of a hypertonic saline solution which would occur if no osmotic movement occurred. The area under this curve is expressed as S.sub.A.Or.hyp.imp. (t.sub.e). Curve 316 represents the change in electrical impedance in the Z.sub.A.hyp (t) measured as a result of the intravenous injection of hypertonic saline solution, with the area under this curve up to t.sub.e being S.sub.A.hyp.imp. (t.sub.e). The difference between the area under curve 324 and the area under curve 316 is indicated in FIG. 31 at area 326, which is expressed as S.sub.A.trans.imp. (t.sub.e). The area 326 may also be expressed as follows:

Because both of the areas S.sub.A.trans.sound (t.sub.e) and S.sub.A.trans.imp. (t.sub.e) are generated by the same volumes of water V.sub.w. (t.sub.e) and salt G.sub.NaCl (t.sub.e), two equations with two unknowns make it possible to calculate these volumes: ##EQU11## where K.sub.w.sound is the sound velocity dilution area from 1 ml of water, K.sub.NaCl.sound is the sound velocity dilution area from 1 gram of NaCl, K.sub.W.imp. is the electrical impedance dilution area from 1 ml of water, and K.sub.NaCl imp. is the electrical impedance dilution area from 1 gram of NaCl. These coefficients are calculated from two pairs of equations, as follows: ##EQU12## where C.sub.isot. it is the concentration of isotonic NaCl and C.sub.hyp. it is the concentration of hypertonic NaCl.

Since changes in the blood electrical impedance Z.sub.A.hyp (t) are proportional to the injected salt concentration in the blood, and since that in turn is proportional to changes in osmolality, then the change in osmolality at t.sub.e, that is, .DELTA.T.sub.b (t.sub.e) it is given by:

where K.sub.p is an unknown coefficient of proportionality. The coefficient K.sub.p can be evaluated by the summation of changes in osmolality and blood electrical impedance over a selected period such as one minute: ##EQU13## This may be restated as follows:

where CO* (1-P) is the water content of the cardiac output.

Saline equations 30 and 32 for .DELTA.T.sub.b (t.sub.e) gives the following result: ##EQU14##

Substituting equation 33 and the calculated volumes of water V.sub.w. (t.sub.e) and salt G.sub.NaCl (t.sub.e) in equation 24, it becomes possible to derive the osmotically active extravascular lung water and the extravascular lung volume from which the water was transferred; that is, the mostly intracellular lung volume, and to derive the extravascular volume to which the NaCl was transferred, which is mostly interstitial lung volume.

In a system configuration such as that of FIG. 11, where an electrical impedance sensor 190, a sound dilution sensor 180, and a flow probe 160 are present in both the arterial and venous parts of the system, lung water may be measured by injecting a volume V.sub.v.inj.isot. of isotonic saline solution at body temperature into the venous port; for example, the port 236. Thereafter, the ultrasound dilution area 330, indicated as S.sub.v.isot.sound in FIG. 32 and the electrical impedance dilution area 332, indicated as S.sub.v.isot. Imp. in FIG. 33, are recorded. In addition, the changes in the arterial line of the blood sound velocity C.sub.A.isot (t) illustrated by curve 334 in FIG. 34 and the changes in electrical impedance Z.sub.A.isot. (t) illustrated by line 336 in FIG. 35 are also recorded. Thereafter, the V.sub.v.inj.hyp. of a hypertonic saline solution having an osmolality T.sub.inj.osm. at body temperature is injected into the venous port 236. The resulting ultrasound dilution area S.sub.V.hyp.sound is illustrated at curve 340 in FIG. 36 and the corresponding electrical impedance solution area S.sub.V.hyp.imp. is illustrated at curve 342 in FIG. 37. Both are in the venous line, and are based on recorded dilution measurements.

The resulting changes in the arterial line of the blood sound velocity Z.sub.A.hyp. (t) illustrated by curve 344 in FIG. 38 and of the arterial line blood impedance Z.sub.A.hyp. (t), illustrated by line 346 in FIG. 39 are also recorded. In this embodiment, the same formulas as those described above with respect to equations 25-33 are used by making the following substitutions:

S.sub.V.isot.sound for S.sub.A.cal.isot.sound

V.sub.v.isot. for V.sub.A.cal.isot.

S.sub.v.hyp.sound for S.sub.A.cal.hyp.sound

V.sub.v.hyp. for V.sub.A.cal.hyp.

S.sub.v.isot.imp. for S.sub.A.cal.isot.imp.

S.sub.v.hyp.imp. for S.sub.A.cal.hyp.imp.

to permit a determination of lung water.

Although the foregoing description utilizes a sound velocity measurement in combination with an electrical impedance measurement, it will be understood that the sensors 180 and 190, as well 180' and 190' can be other types of dilution measurement sensors and detectors. For example, an optical dilution sensor may be combined with a sound dilution sensor or an optical dilution sensor can be combined with an electrical impedance dilution sensor to calculate extravascular lung water in an analogous way. By using a single type of sensor in multiple measurements, relative changes in lung water can be calculated.

Another way to calculate lung water is to inject intravenously a diffusible indicator; for example, cold water, at a temperature other than body temperature, and a non-diffusible indicator such as, for example, a green dye (Oppenheimer et al, "Thermal-dye lung water measurements: effects of edema and embolism", J. Surge. Res. 26; pp. 504-512 (1979)). Extravascular lung water can then be calculated as follows:

where V.sub.therm. is the volume of the diffusible and non-diffusible indicators and MTT.sub.DIF. -MTT.sub.Nondif. represents the mean transit times of the diffusible indicator and the non-diffusible indicator, respectively. The mean transit time is calculated from the following equation, as is known: ##EQU15## where F(t) is the dilution curve. For this variant of lung water calculation, two sensors are needed, one for the diffusible and one for the non-diffusible indicator.

In another embodiment of the invention, lung water can be determined using only a single sound velocity sensor. As indicated in equation 1, blood sound velocity is a function of temperature, and this property is used to record thermal changes caused by blood dilution. Two intravenous isotonic injections of an indicator solution are made, one at body temperature, and the other at a different temperature; for example, colder than body temperature. The sound velocity dilution curve recorded in the arterial line for the intravenous isotonic injection is illustrated by curve 350 in FIG. 40. The sound velocity dilution curve due to the cold injection is illustrated by curve 352 in FIG. 41. The difference between the curve 352 representing the cold injection C.sub.A.isot.sound.cold (t) and the sound velocity dilution curve 350 representing the body temperature C.sub.A.isot.sound (t) gives the sound velocity changes that are related only to temperature changes C.sub.A.sound.thermal (t), as indicated in the following equation:

Using C.sub.A.sound.thermal (t) as the curve of the diffusible indicator and C.sub.A.isot.sound (t) as the curve of the non-diffusible indicator in the formulas of equations 34 and 35 results in the extravascular lung water volume measured without temperature measurements. In the same way, an electrical impedance sensor may be used to obtain lung water volume.

Although the present invention has been described in terms of preferred embodiments, various modifications may be made without departing from the true spirit and scope thereof, as set forth in the following claims: