Today, a number of techniques is available for human cardiac output monitoring. They can be based on invasive or non-invasive measurement techniques and are developed for continuous or non-continuous monitoring purposes. One continuous monitoring method is described as “Modelflow” in U.S. Pat. No. 5,400,793. This method couples the pulse-type blood-stream pressure signal derived from the aorta and calculates the flow where the aorta is regarded as a transmission line supplemented by a windkessel compliance.
Cardiac output monitoring, for instance, of patients who undergo major surgery, or are at intensive care, provides valuable information on patient status. Non-continuous methods are reliable, but are often invasive and require considerable skill in conducting the measurement. Continuous cardiac output monitoring is therefore becoming more and more popular. Therefore, methods and techniques for continuous and non-invasive cardiac output monitoring are important for future developments in the field of cardiac output monitoring.
For non-invasive situations, the above mentioned “Modelflow” technique requires the availability of a non-invasive blood pressure signal which is not always available at the operation- or intensive care units. Furthermore, the Modelflow technique must be calibrated using other complex cardiac output measuring techniques. The invention aims at providing a reliable continuous measurement of cardiac output, at the same time obviating a necessity for performing complex or invasive measurements or calibration. The invention aims at providing a reliable monitor of cardiac output, both for output changes and absolute output measurements.
To this end, the invention provides a method according to the features of claim 1. Specifically, by providing a circulation model wherein a measured CO2 partial pressure and oxygen uptake to a heart stroke volume per breath are related to the cardiac output, minimal attributes are necessary for providing a reliable measurement. In another aspect, the invention provides a ventilation unit for actively or passively ventilating the lungs and comprising a breathing mask piece according to claim 22. In particular, the invention provides a ventilation unit comprising:
a respiratory detector for measuring an expiratory tidal volume and respiratory rate;
an oxygen sensor for measuring an oxygen uptake from the expiratory tidal volume;
a CO2 sensor for measuring a CO2 partial pressure in the expiratory tidal volume; and
a processor programmed in consistency with a circulation model, for relating a measured CO2 partial pressure, expiratory tidal volume and oxygen uptake to a heart stroke volume per breath;
wherein the respiratory detector, the oxygen sensor and the CO2 sensor are arranged to be coupled to said processor for inputting a measured expiratory tidal volume; a measured oxygen uptake VO2 and a measured CO2 partial pressure, and the processor arranged to output a heart stroke volume per breath to an output unit consistent with said circulation model.
The inventors have found that variations of the CO2 partial pressures can be modelled using a mathematical model of the circulation system. Accordingly, by the invention, using capnography and a breath-to breath model, the cardiac output of the heart can be monitored non-invasively in a continuous manner.
In a preferred embodiment, said circulation model defines a distribution of apical and basal lung segments, each segment defining a predetermined ventilation perfusion ratio (V/Q), and where the heart stroke volume per breath n SVn is calculated to be consistent with an estimated fraction b of CO2 in air in each segment, derived from a measured end tidal CO2 partial pressure and pulmonary oxygen uptake. The ventilation and perfusion rates that the model takes into consideration may incorporate physical or pathological conditions of the circulation in a subject that is monitored.
On standing up, the distribution of blood flow over the lung changes due to hydrostatic pressure differences from apex to base in the upright position. Distribution of tidal volume also changes, due to an altered pressure/volume relationship of the lung compartments. Consequently, the overall ventilation/perfusion ratios changes. Results indicate that the model accurately estimates end tidal variations of the CO2 partial pressures during posture change and tracks spontaneous variations in partial CO2 pressure (PCO2) in a fixed body position, substantially ascribed to changes in cardiac output. It is noted that parts of this model were disclosed in a scientific paper of the inventor: Janneke Gisolf et al, “Tidal Volume, Cardiac output and functional residual capacity determine end-tidal CO2 transient during standing up in humans”, J Physiol. 554.2 p 579-590.
Further features and benefits will become apparent from the description read in conjunction with the figures.
In the figures:
In
The processor 10 is programmed in consistency with a circulation model which defines a distribution of apical and basal lung segments, each segment defining a predetermined ventilation perfusion ratio (V/Q), and where the heart stroke volume per breath n is calculated to be consistent with an estimated fraction of CO2 in air in each segment, derived from the measured end tidal CO2 partial pressure and pulmonary oxygen uptake.
The ventilation and perfusion rates that the model takes into consideration may incorporate physical or pathological conditions of the circulation in a subject that is monitored. The model is further detailed in
As illustrated in
As illustrated in
The respiratory quotient (RQ), defined as the ratio of carbon dioxide production (VCO2) to VO2, normally between 0.7 and 1.0, was set fixed at 0.9. Alternatively, variable values of the RQ value can be inputted in the model using known measurement techniques. The consequences of variations in oxygen uptake are 2-fold:
1. oxygen uptake is related to basal metabolism, and is related to CO2 production. For a resting, supine measurement there will be little variation in oxygen uptake.
2. the level of oxygenation of blood determines its ability to carry CO2 (known as the ‘Haldane effect’. For a normal, resting measurement the oxygenation will be optimal and the CO2 uptake will be defined by the Equation 1 detailed below.
The lung capillary volume and the small venule volume are lumped together, as gas exchange occurs in both. The major arteries of the lung are included in the venous compartment; the major veins of the lung are included in the arterial compartment. In a practical example, the total blood volume of 5.5 l is distributed over Vv (4.0 l), Va (1.3 l) and Vcap (0.2 l). The segmented model may include the effects of gravity to gravity-induced blood perfusion gradient in the lung. In the supine position, SV and VT are distributed equally over all compartments. With nine compartments, in the supine position each lung compartment receives one-ninth of the breath-to-breath SV and VT. In the upright position there is an apical-to-basal perfusion and ventilation gradient, with increased perfusion and ventilation at the lung base. The perfusion gradient is steeper than the ventilation gradient, resulting in a 7.9-0.8 apical-to-basal V/Q gradient. Furthermore, different values for anatomical dead space VD may be used going from supine to upright respiratory positions, for instance in a range varying from +53 ml (anatomical) to +81 ml (physiological).
Table 1 defines a distribution of stroke volume (SV), tidal volume (VT), functional residual capacity (FRC) and lung capillary blood volume (Vcap) per lung segment k, in the supine and standing position, that can be included in the model. Upright distributions are based on measurements previously performed by West J B (1962): “Regional differences in gas exchange in the lung of erect man”. J Appl Physiol 17, 893-898 West (1962). For the upright position, the FRC is increased with respect to the supine position.
For the purpose of tracking short-term end tidal PCO2 variations with posture change, data was selected starting 150 s prior to standing up and ending 150 s after standing up. Mean arterial blood pressure was measured with a Finapres (Model 5; Netherlands Organization for Applied Scientific Research, Biomedical Instrumentation, TNOBMI). The cuff was applied to the midphalanx of the middle finger of the dominant arm, which was placed at heart level. Beat-to-beat changes in SV were estimated by modeling flow from arterial pressure (Modelflow, TNOBMI). This method computes an aortic waveform from a peripheral arterial pressure signal using a non-linear 3-element model of the aortic impedance (Jellema et al. 1999; Harms et al. 1999). Cardiac output was the product of SV and HR. To obtain absolute values of Q to calibrate Modelflow Q, a Fick-determined Q was obtained from arterial and central venous O2 content and the VO2 in the supine and in the standing position. Absolute values of Q were used to calibrate Modelflow Q, averaged during 150 s in the supine position, and during 150 s of standing. Breath-to-breath online gas analysis was performed using a Medical-Graphics CPX/D metabolic cart. Respiratory gas was sampled continuously from a mouthpiece and partial gas pressures were obtained from a Zirkonia oxygen analyser (accuracy±0.03% O2) and a non-dispersive infrared sensor for CO2 (accuracy±0.05% CO2) that thus delivered VO2, VCO2and PETCO2.
For the definitions of the various quantities used in the model equations reference is made to Table 2.
The CO2 equilibrium curve relating blood CO2 content ([CO2]) to blood partial CO2 pressure (PCO2) is described as [CO2]=f (PCO2),with
f(x)=0.53(1.266−exp(−0.0257x)) Equation 1
To compute PCO2 from [CO2] in blood, we use the inverse function
f−1(x)=−ln(1.266−(x/0.53))/0.0257 Equation 2
To convert PCO2 in air (mmHg) to [CO2] (%), we use the conversion factor c, which amounts to 0.1316% mmHg−1. The distribution of SV and VT over each lung compartment k (k=1 . . . 9) is described by functions g and h, respectively. These functions, which are different for the supine and upright positions and yield the fractions for SV and VT listed in Table 2, are given by
g(k)= 1/9(in the supine position)−0.0205+0.0263k (in the upright position) Equation 3
and
h(l)= 1/9 (in the supine position) 0.226(1.102−exp(−0.1063k) (in the upright position)
Each lung compartment's share of FRC, Vcap and VD is given by the weight function
w(k)=0.10055(1.36708−exp(−0.3393k)) Equation 5
which yields the fractions for FRC and Vcap listed in Table 2.
For each breath n, a variation in CO2 in said venous compartment Vv is expressed by the amount A that arrives from the arterial compartment Va plus the amount B of CO2 created by the basal metabolism minus the amount C that exits the venous compartment; where the sum of CO2 created by the basal metabolism is expressed as a function of the oxygen uptake VO2 per breath. For each breath a, the venous CO2 content ([CO2]v,n) is calculated from its previous value [CO2]v,n−1 according to Equations 6-9. The amount of CO2 in the venous compartment increases by the amount that arrives from the arterial compartment (A) and the amount created by the basal metabolism (B), and decreases by the amount that leaves the compartment (C). Thus, we have
[CO2]v,n−[CO2]v,n−1+(A+B−C)/Vv
where
C=[CO2]v,n−1SVn Equation 7
A=[CO2]a,n−1SVn Equation 8
with [CO2]a denoting the arterial CO2 content, and
B=VO2, nRQ(TRESP,n/60) Equation 9
where VO2,n is the oxygen extraction for breath n (in ml min−1) and RQ is the respiratory quotient, which is set at 0.9 (the average as approximated from subject data, by dividing VCO2by VO2). The term is multiplied by the breath duration (in min) (TRESP,n/60) to estimate the CO2 produced per breath.
For each breath n, a variation in CO2 in said venous compartment Vv is expressed by the amount A that arrives from the arterial compartment Va plus the amount B of CO2 created by the basal metabolism minus the amount C that exits the venous compartment; where the sum of CO2 created by the basal metabolism is expressed as a function of the oxygen uptake VO2 per breath.
The arterial blood CO2 content for breath n ([CO2]a,n) is calculated from its previous value [CO2]a,n−1 according to Equations 10-12. The amount of CO2 in the arterial compartment increases by the amount of CO2 arriving from the lungs (D) and decreases by the amount of CO2 leaving the arterial compartment (E)
[CO2]a,n=[CO2]a,n−1+(D−E)/Va Equation 10
The amount D can be estimated from the end-tidal partial CO2 pressure in each lung compartment k (PkCO2k=1 . . . 9) through
Where f is the above function that relates blood CO2 content to the blood partial CO2 pressure and g is the above function that defines the distribution of SV over the nine lung compartments. The amount E is given by
E=[CO2]a,n−1SVn Equation 12
For each breath n, a CO2 amount in each segment k of said segmented lung model is expressed as the amount F of CO2 in the lung capillaries Vcap, in the functional residual capacity FRC, and in the anatomical dead space VD, as a function of an estimated CO2 partial pressure PkCO2n in the segments k; plus the amount G of CO2 carried to the lungs from the venous compartment by the heart stroke volume SVn; and where the estimated CO2 partial pressure PkCO2n in the segments is expressed in relation to an estimated fraction b of CO2 in air in each segment k.
The PCO2 of blood draining the lungs (PtcCO2) is dependent on the gravity-induced perfusion and ventilation gradients, as described by the above functions g and h. For each breath, the PCO2 in each lung segment k (PkCO2,n) is calculated according to Equations 13-18. At FRC, the amount of CO2 in lung segment k (F) is determined by the CO2 content in the lung capillaries, in the FRC and in the VD
F=f(PkCO2,n−1)w(k)Vcap+cPkCO2,n−1w(k)FRC+cPETCO2,n−1w(k)VD
with the weight function w and conversion factor c as described above. The contribution of CO2 in dead space (the right-most term) is computed noting that end-tidal air from the previous breath is returned to the lungs from dead space. The amount of CO2 carried to the lungs from the venous compartment (G) is given by
G=[CO2]v,n−1SVng(k) Equation 14
The ratio a of [CO2] in blood and [CO2] in air is approximated from the previous breath, n−1, according to
a=f(PETCO2,n−1)/(cPETCO2,n−1) Equation 16
The ratio b of the end-tidal amount of CO2 in air and the total amount of CO2 is given by
b=(w(k)FRC+h(k)VTn)/(a(w(k)Vcap+g(k)SVn)+w(k)FRC+h(k)VTn) Equation 16
The end-tidal [CO2] in each lung compartment k is determined by the total amount of CO2 (F+G), which is distributed over air and blood with ratio b, and the end tidal volume of air in compartment k
[CO2]k,n=b(F+G)/(w(k)FRC+h(k)VTn) Equation 17
A simple conversion using the above constant c then yields PkCO2,n. The PETCO2 depends on the distribution of tidal volume, which is given by the fraction h(k),k=1 . . . 9, and differs between the supine and the standing position, and is computed as
In conjunction with the above described model, through direct calculation, or via iterative testing, or via look up tables the heart stroke volume SVn can be determined in consistency with the model, for each breath.
Application of the Model for the Supine Position
For the supine position, measured CO2 partial pressure in the expiration can be related to the CO2 partial pressure in the lung-compartments. CO2 pressure as a percentage of total air pressure corresponds with CO2 concentration (also as percentage). With respect to CO2 equilibrium curve the CO2 pressure in blood corresponds with a much greater CO2 concentration (percentage), which can be calculated as according to the function of Equation 1 here above. For the supine position, the partial CO2 pressure in the lungs is equal for all segments and also to the measured CO2 partial pressure in the expiratory volume. Hence, the CO2 concentration in blood can be obtained straightforwardly. The blood volume responsible for producing the CO2 can be expressed as the sum of the cardiac output SV plus the capillary volume Vcap. The amount of dissolved CO2 is determined by the CO2 production, which can be estimated. Thus, the cardiac output can be derived by measuring expired CO2 air partial pressure, relating this to a CO2 concentration in the blood using Equation 1, and deriving a cardiac output by dividing the CO2 production by the CO2 concentration in blood, and subtracting an estimated capillary volume of the lungs:
SV=((CO2 production per breath/[CO2]blood)−Vcap) Equation 19
Using a measured respiratory frequency and the heart stroke rate, this value can easily be converted to a cardiac output value per heart stroke. It will be clear to those skilled in the art that the invention is not limited to the exemplary embodiments described with reference to the drawings but may comprise all kinds of variations thereof. Such variations are deemed to fall within the scope of protection of the appended claims.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/NL2004/000783 | 11/5/2004 | WO | 00 | 2/15/2008 |