Arterial CO.sub.2 Monitor and Closed Loop Controller

This invention related to a method and apparatus for continuously and non-invasively monitoring arterial blood CO.sub.2 partial pressure (PaCO.sub.2) of artificially ventilated patients.

Mechanical ventilation is required by patients in an intensive care unit who are unable to control their own respiration. The rate of ventilation must be adjusted so that arterial CO.sub.2 is within a desirable range. Conventionally clinicians adjust the ventilator settings based on periodically drawn blood samples. In order to monitor rapidly changing PaCO.sub.2 (for monitoring or closed loop control purposes), a continuous and non-invasive monitor is desirable. Known transcutaneous transducers are non-invasive but require heating of the patient's skin to 44.degree. C. and a long stabilization time of 30 minutes which renders them unsatisfactory for continuous monitoring. The known method of assuming a constant arterial to end-tidal CO.sub.2 difference is not reliable during ventilation/perfusion changes, and attempts to implement closed loop ventilation control have failed largely due to the inability to continuously and non-invasively observe the variable to be controlled, that is, the PaCO.sub.2.

Thus, the direct methods of monitoring PaCO.sub.2 are invasive, and indirect methods are not reliable, particularly because end-tidal CO.sub.2 is influenced by deadspace, which is an unmeasurable quality.

It would thus be desirable to provide a method and apparatus for providing a continuous and substantially non-invasive PaCO.sub.2 estimation.

Briefly, and in general terms the invention provides a method and apparatus for continuously and non-invasively monitoring arterial blood CO.sub.2 partial pressure (PaCO.sub.2) of artificially ventilated patients, by monitoring a patient's breath and determining PaCO.sub.2 based upon a determination of a deadspace ratio, which is the ratio of the alveolar deadspace to alveolar tidal volume. The method generally comprises the steps of continuously monitoring, measurable parameters of a patient's breath; obtaining an input value of PaCO.sub.2 from a blood sample of the patient and using the patients breath parameters and the input value to calculate the deadspace ratio; and continuously determining PaCO.sub.2 based on the assumption that the deadspace ratio subsequently remains constant.

Decision rules obtained from other measurable data are preferably also used to identify the onset of changes in the deadspace ratio, and a new deadspace ratio is then determined from the patient's breath parameters and a further input value of PaCO.sub.2 from the patient's blood sample.

The determination of PaCO.sub.2 is preferably based upon the equation ##EQU1##

where V.sub.D.sup.alv is the alveolar deadspace,

V.sub.T.sup.alv is the alveolar tidal volume,

PE.CO.sub.2 is the mixed-expired CO.sub.2 from the alveolar tidal volume, and

PiCO.sub.2 is inspired CO.sub.2.

The mixed-expired CO.sub.2, inspired CO.sub.2, alveolar tidal volume and the alveolar deadspace are the measurable parameters of the patient's breath.

The other measurable data used to determine decision rules for identifying changes in the deadspace ratio are preferably related to lung mechanics and trends in CO.sub.2 production.

The method preferably further involves adjusting patient ventilation based on the determined value of PaCO.sub.2.

In another aspect of the invention, an apparatus is provided for continuously and non-invasively monitoring arterial blood CO.sub.2 partial pressure (PaCO.sub.2) of artificially ventilated patients. The apparatus preferably includes a capnograph for monitoring continuously measurable parameters relevant to a patient's breath and providing data relating thereto, and means for determining a deadspace ratio connected to the capnograph to receive the breath parameter data and adapted to receive information relating to the PaCO.sub.2 of a blood sample of the patient based upon the PaCO.sub.2 information and the breath parameter data. Means are also preferably provided for continuously determining PaCO.sub.2 based on the deadspace ratio, and the assumption that the deadspace ratio remains subsequently constant. However, means are also preferably provided for further receiving decision rules enabling identification of the onset of changes in the deadspace ratio to thereby signal the need for a further blood sample to re-calculate the deadspace ratio.

The apparatus is preferably connected to a mechanical ventilator to control operation of the ventilator based on the values of PaCO.sub.2 determined by the apparatus.

These and other aspects and advantages of the invention will become apparent from the following detailed description, and the accompanying drawings, which illustrate by way of example the features of the invention .

FIG. 1 is a basic functional block diagram of the apparatus of the invention;

FIG. 2 is a block diagram in the form of a software flow chart for the apparatus of FIG. 1;

FIG. 3 is a graph of airwave CO.sub.2 partial pressure versus expired volume for each breath of the patient;

FIG. 4 is a schematic diagram showing the main method steps according to the invention;

FIG. 5 is a graph of test results showing PaCO.sub.2 estimation; and

FIGS. 6(A) and 6(B) show a detailed functional block diagram of the apparatus of the invention.

In FIG. 1 the outlet port of capnograph 10 is connected, via connection 11, to analogue to digital (A/D) converter 12. The capnograph 10 may, for example, be a HP-78356 and the A/D converter 12 may comprise analogue devices such as an RT1-815.

The A/D converter 12 is connected via connections 13 and 14 to computer 15 which may, for example, be an IBM compatible PC-AT286. The connection 13 is to the interrupt port of the computer 15 and the connection 14 is to the data port.

A mechanical ventilator 16 such as a Puritan-Bennett 7200 is typically connected both to the A/D converter 12 and the computer 15 as shown. The connections 17 and 18 are to the analogue ports of the ventilator and carry signals related to the pressure and flow respectively. The connection 19 is to the digital port of the ventilator 16 and provides an RS-232 link to the computer 15.

Referring now to FIG. 4, the method in fact involves two major modules, the first of which is the PaCO.sub.2 estimator 20 and the second of which is the decision module 21. The estimation module 20 receives information on input 22 related to measurements taken from each breath of patient 23, that is, by a breath-by-breath analysis. The decision module 21 receives information on input 24 obtained from other measurable quantities relevant to the patient such as lung mechanics and trends in CO.sub.2 production for example, and contains decision rules obtained by previous experimentation. The rules are implemented as the rule-base of an Expert system.

The PaCO.sub.2 estimator 20 is described by equation (1). ##EQU2##

where V.sub.D.sup.alv is the alveolar deadspace,

V.sub.T.sup.alv is the alveolar tidal volume,

PE.CO.sub.2 is the mixed-expired CO.sub.2 from the alveolar tidal volume, and

PiCO.sub.2 is inspired CO.sub.2.

The deadspace ratio is the ratio of alveolar deadspace to alveolar tidal volume.

The various parameters may be obtained from the plot of airwave CO.sub.2 partial pressure versus expired volume for each breath a shown in the single breath test graph of FIG. 3 On the graph PE.sup.1 CO.sub.2 is end tidal CO.sub.2, and V.sub.T.sup.alv is tidal volume involved in gas exchange. The alveolar deadspace (represented by area Y) is that part of inspired gas which reaches the alveoli but does not take part in gas exchange. VD.sup.aw, the airway deadspace, is the point of maximum inflection of the plot. From FIG. 1, PE.CO.sub.2 and PiCO.sub.2 may be found using Equation (2)-(5). ##EQU3## where V is the respiratory rate. VCO.sub.2 is carbon dioxide production which can be calculated each minute by integrating the CO.sub.2 fraction (FCO.sub.2 ) and the flow signal, as shown in Equation (6)-(7). ##EQU4## where V is flow, FCO.sub.2 is fraction of CO.sub.2, PCO.sub.2 is the capnograph signal, P.sub.airway is airway pressure, P.sub.B is barometric pressure and P.sub.H2O is vapor pressure.

Breath-by-breath processing yields the mixed expired CO.sub.2, inspired CO.sub.2, alveolar tidal volume and the airway deadspace. The only unknowns being PaCO.sub.2 obtained after a blood sample analysis, the deadspace ratio can be calculated. Assuming that the deadspace ratio remains subsequently constant, further PaCO.sub.2 can be calculated using Equation (1). A PaCO.sub.2 estimate is calculated once every minute based on the average of the breaths in the minute.

Certain corrections are needed when implementing the system as set forth below:

(a) Flow signals have to be corrected form BTPS (Body Temperature Pressure Saturated) to STPD (Standard Temperature Pressure Dry).

(b) Correction must be made for delay-time between the flow signals and capnograph signals. Delay time is found by simple breath-holding and rapid expiration through the airway tubing, and lining up the start of flow and capnograph signals.

(c) Correction for compliance of the airway tubing. Flow due to compliance volume is subtracted form the analog flow signal, using ##EQU5## where P=airway pressure, C=tubing compliance,

V=flow due to compliance.

(d) Correction for rebreathing is done by continuously integrating flow and CO.sub.2 fraction, using Equation (6).

(e) Correction of capnograph signal for vapor pressure and airway pressure is done by Equation (7).

As mentioned above the decision rules are obtained by experimentations to determine rules which indicate a change in the deadspace ratio whereby the system may signal that a new blood test is required.

The following are the decision rules derived:

(i) If Alveolar minute volume increases and CO.sub.2 production decreases, deadspace ratio may have changed.

(ii) If Alveolar minute decreases and CO.sub.2 production increases, deadspace ratio may have changed.

(iii) If Alveolar minute volume increases and arterial or end-tidal CO.sub.2 increases, deadspace may have changed.

(iv) If Alveolar minute volume decreases and arterial or end-tidal CO.sub.2 decreases, deadspace may have changed.

It is possible to derive further rules to indicate a change in deadspace ratio. For example, changes to airway resistance, peak airway pressure (PAP), peak flow, SaO.sub.2, inspiratory to expiratory ratio, and positive end-expiratory pressure (PEEP) should indicate a change in the deadspace ratio. By automatically recording these parameters during a clinical trial, including blood test results, correlation between the change in parameters and change in deadspace ratio can be performed.

If desired, a closed-loop control of ventilation may be implemented based on the predicated PaCO.sub.2. The controller in this case is a set of rules which decides on the tidal volume and respiratory rate settings for the mechanical ventilator, to achieve and maintain PaCO.sub.2 at a set-point. The controller rules are based on existing clinical protocol for ventilator settings.

The control algorithm is presented below:

First, a PaCO.sub.2 setpoint has to be determined as follows:

Given the pH value from the most recent blood gas result, if pH is between 7.36 and 7.44, the PaCO.sub.2 setpoint is 40 mm Hg, the default value. If pH exceeded the limits, the following equation is used to calculate a new PaCO.sub.2 setpoint.

where pH=7.4, and HCO.sub.3.sup.- is from the most recent blood gas analysis results. The PaCO.sub.2 setpoint can also be set by the clinician, overriding the above calculations.

Next, ventilation settings for the next 5 minutes can be set by the following equation.

where PaCO.sub.2 is the setpoint (V.sub.T.sup.alv * f) is the alveolar minute ventilation needed to achieve the setpoint, PaCO.sub.2 ' is the latest estimation and (V.sub.T.sup.alv '* f') is the latest minute ventilation.

To decide on the actual V.sub.T and f, from the minute volume, the following procedures are followed:

Increase V.sub.T and keep f constant, so that

(1) V.sub.T is not smaller then 500 ml.

(2) V.sub.T is not bigger than 1000 ml.

(3) Peak airway pressure (PAP) is not greater than 40.

If (2) or (3) cannot be satisfied, keep V.sub.T constant and increase f such that

(a) f is not smaller than 10 bpm.

(b) f is not bigger than 20 bpm. (and not less than 10 bpm.)

(c) Inspiratory to expiratory ratio (I:E) is smaller than 1:2.

If (b) or (c) is exceeded, Peak air flow (PAF) should be increased.

If the required minute volume cannot be achieved without exceeding the limits in any one of VT/f/PAP/mean pressure/PAF/I:E, a warning message will be displayed to alert the clinician.

As is evident above, the system of the invention is implemented, according to one embodiment, on a PC-AT computer. Analogue flow and pressure signals from the Puritan-Bennett 7200 ventilator and airwave CO.sub.2 signals from the HP-78356 capnograph are analogue-to-digital converted and processed in real time using the computer. Each breath is checked to reject unphysiological waveforms, such as incomplete spontaneous breaths, before further processing. A PaCO.sub.2 estimate is calculated every minute and the decision rules are invoked.

Tests have been performed to verify all calculations made in formulating the decision rules. CO.sub.2 production and mixed-expired CO.sub.2 were tested against the Douglas bag method. Air way deadspace was checked by introducing a known deadspace into the ventilator circuit. Corrections had to be made to account for airway tubing compliance, time delay between the flow and capnograph signal, and rebreathing. To test the accuracy of PaCO.sub.2 estimation, clinical trials were carried out on ICU patients. For each patient, tidal volume and ventilation rate were changed in different combinations to a maximum of 30% of the initial settings. After each change and a stabilization period of about 10 minutes, a blood sample was drawn to check the estimation. Each clinical experiment involved 6-8 manipulations, over 90 minutes.

To identify the factors and the degree that they affect the alveolar deadspace ratio, specific procedures are incorporated into the clinical experiments to change the deadspace. Posture of the patient is changed by tuning the patient or inclining of the bed. Another procedure has been to administer bronchodilators. Various measurable parameters are recorded during the experiments for correlation with deadspace changes. These include airway compliance and resistance, peak airway pressure, peak air flow, inspiratory time, positive end-expiratory pressure (PEEP), inspiratory to expiratory ratio (I:E), slope of the CO.sub.2 versus expired-volume waveform, end-tidal CO.sub.2 and SaO.sub.2.

The test results show reliable estimation (.+-.5 mmHg) of PaCO.sub.2, even when deadspace ratio changed by up to 30% of the initial value. This indicates that the estimator is robust to some changes in the deadspace ratio.

Test results also showed that deadspace ratio change can be expected when alveolar tidal volume and frequency changes are not: followed by expected changes in end-tidal CO.sub.2, estimated PaCO.sub.2 or CO.sub.2 production. Results from a trial are presented in FIG. 5. At point A, the increase in alveolar tidal volume and ventilation rate product (alveolar minute-volume) is not followed by a drop in both end-tidal CO.sub.2 and estimated PaCO.sub.2, indicating a blood test is needed. At point B, alveolar minute-volume decrease is not followed by an increase in end-tidal CO.sub.2. In each case, the new estimation system correctly identifies the deadspace ratio change and estimates PaCO.sub.2 reliably, compared to using the traditional method based on a constant arterial--end tidal difference.

Changes to compliance and resistance, peak airway pressure, peak flow, SaO.sub.2 and I:E should indicate a change in the deadspace ratio but more results are needed before these relationships can be quantified.

The results show that the PaCO.sub.2 estimator is sufficiently rebust to permit continuous estimation for a wide range of ventilator settings. For large deadspace changes, the experimentally derived rules can be relied upon to signal for a blood gas test. Nevertheless, further clinical runs are necessary to cover all possible cases of deadspace ratio changes. By incorporating an Expert System, the knowledge base may be easily extended as more clinical data becomes available.

For the purpose of an even fuller understanding of the invention, the following description provides the pseudo code for programming the apparatus according to the invention. The program should be read in conjunction with FIG. 6 which is a self-explanatory functional block diagram of the apparatus.

System setup: Hardware setup and parameters initialization.

Repeat

Timer: Keeping track of real-time, using each interrupt service form fixed-frequency AtoD conversion as time counter. Sets flags after each minute ("onemin" flag) and five minutes ("fivemin" flag). If "one min" flag is true,

Ventrequest: Request data from ventilator by sending request codes to serial port. The first request is "SPD" send patient data

FndCO.sub.2 prod: Request calculate carbon dioxide production and mixed expired carbon dioxide for the minute.

Checkcomq: checking comm. queue. If queue is not empty,

Readvent: Read characters received from ventilator.

If calculation above are satisfactory,

Sbt.sub.-- CO.sub.2 calc: Request calculation of physiological parameters and PaCO.sub.2 prediction.

If "fivemin" flag is true and calc. above are satisfactory,

InvokeExpert: Test if deadspace ratio changes.

InvokeControl: Control rules.

Checkkey: Check keyboard for keyhit. If Keyhit,

Processkeyhit: "M" for marking of blood taken.

Entering of PaCO.sub.2, pH, HCO.sub.3.

"Q" to exit from program.

If hardware test failed,

Safe-exit.: Disables all interrupts, turn off

A/D conversion, and close all files before exiting.

If 0.5-4.0 seconds have elapsed, process data from queue:

Toscrn: Display on screen.

Tostore: Store in file.

FndCO.sub.2 prod: ongoing CO.sub.2 production calc.

Sbt.sub.-- CO.sub.2 calc: Ongoing SBT calc.

Until exit from system

Initialize variables.

Initialize graphics (axes for plotting, etc.)

Hardware setup: Setting up interrupts for digital and analog interfacing.

Setup a circular buffer (AtoDqueue) for storing real-time data.

Set up a circular buffer (Commqueue) for storing character strings received from serial port.

Program Analog to Digital Conversion card to convert at set frequency (100Hz), and to generate an interrupt on completion of each conversion, for the first channel of three. (The three channels are the airway flow, pressure and capnograph signals.)

Set and enable interrupt vector (Interrupt Service Routine: SetAtoD) for A/D card, to read analog flow and pressure from the ventilator, and capnograph waveform from the capnograph.

Set and enable interrupt vector (Interrupt Service Routine: SetComm) for serial port, to read digital data from ventilator.

Increment timing counter, for "Timer" routine.

Read converted data from port.

Request and read the other two channels. (The three channels are flow, pressure and capnograph waveform).

Put delay time between flow, pressure waveforms and the capnograph waveform to synchronize the signals.

Put data into the circular buffer (AtoDqueue).

Put characters received in circular buffer (Commqueue).

If "onemin" or "fivemin" flag has been set, clear it (Token removal to ensure that the token is passed around the real-time loop only once).

Checks timing counts. (Generated by ISR, SetAtoD).

If one minute ha elapsed, set "onemin" flag.

If five minutes have elapsed, set "fivemin" flag.

During ongoing calculations (by the data count),

Repeat from queue head to tail,

Read CO.sub.2 (mmHg), airway flow. (Lpm) and pressure signals (mmHg).

Calculate CO.sub.2 fraction: CO.sub.2 fraction=CO.sub.2 /(PB-P.sub.H20 +pressure), where PB=760, P.sub.H20 =47.

Correction for flow:

Flow(Lpm)=flow:

Compliance * (pressure-previous pressure) * 60/PTTOMS, where

Compliance=5.17 ml/mmHg, PTTOMS=10.

If airway flow >+1.2 1pm (expiration phase),

CO.sub.2 flow=CO, flow+flow*CO.sub.2 fraction.

(Integration of airway flow and CO.sub.2 fraction)

Expvol=expvol+flow. (Calculation of expired volume in the breath)

If request for a minute's CO.sub.2 production,

CO.sub.2 production=CO.sub.2 flow * BTPStoSTPD * PTTOMS/60.

where BTPStoSTPD=0.8262.

Expvol (ml)=expvol * BTPStoSTPD * PTTOMS/60.

Mixed-expired

CO.sub.2 (mmHg)=CO.sub.2 production,863*0.8262/Expvol.

Return

Read characters from circular buffer (Commqueue). If carriage return is encountered (ie. end of message received),

Check whether the heading of message is SPD (send patient data), SLM (send lung mechanics) or SVS (send ventilator status).

For each message received, error checking by size of message, and whether various parameters read are within physiological range.

If SPD is received, SLM will be requested.

If SLM is received, SVS will be requested.

If SVS is received, and if all data are acceptable, a flag will be set to indicate completion of ventilator requests.

(Parameters read from ventilator include: RR (rate), MV (minute volume) , MAP (mean air pressure), IE (inspiratory to expiratory ratio), VT (tidal vol), SMV (spontaneous MV), PAP (peak airway pressure), SRR (set rate), SVT (set VT), PIF (peak inspiratory pressure), PEEP (positive end expiratory pressure), DMC (Dynamic Compliance), DMR (Dynamic Resistance)).

If ongoing calculation,

For each breath, if flow <-1.2Lpm (inspiration phase),

Tinsp (inspiratory time) is counted.

For each breath, if flow >+1.2Lpm (expiration),

Correct for compliance flow.

Calculate tidal volume(Vt):

Expvol(ml)=expvol+(flow*PTTOMS/1000).

Plot Single Breath Co.sub.2 Expiration curve (airway

CO.sub.2 (mmHg) vs expvol(ml)).

Calculate gradient of airway CO.sub.2 vs expvol.

Record the maximum inflection point of CO.sub.2 vs expvol;

This is the airway deadspace, VD.sub.airway.

At end of expiration (flow<-1.2Lpm),

Reject data from the breath if waveform is unphysiological (Less than 50 data points, or maximum inflexion point is not found).

Determine end-tidal CO.sub.2, the maximum airway CO.sub.2.

Calculate the plateau slope of CO.sub.2 vs expvol.

Calculate the area under the Single Breath CO.sub.2 (area X), by summation.

Parameters calculated for this breath are summed with parameters from previous breaths within a minute. The parameters are: the number of breaths, areaX, VD.sub.airway, tidal volume (V.sub.T) , end-tidal CO.sub.2 (ETCO.sub.2), plateau phase slope, Tinsp and number of rejected waveforms.

If request for a CO.sub.2 prediction,

If arterial CO.sub.2 sample is drawn (marked) but not available yet,

The minute's parameters are kept in a buffer.

If arterial CO.sub.2 is available,

An average for a 11 parameters over the last minute is calculated. (The average value of each parameter is the summation/total number of breaths.)

From the average, the following calculations are done:

V.sub.t.sup.alv (alveolar tidal vol)=V.sub.T -VD.sub.airway.

PE*CO.sub.2 (bag mixed expired CO.sub.2)=area X/V.sub.T.sup.alv.

PE.sup.bag CO.sub.2 (bag mixed expired CO.sub.2)=

CO.sub.2 prod,863/(V.sub.T *f)

PE.sup.bag CO.sub.2 =PE.sup.bag CO.sub.2 * V.sub.T /V.sub.T.sup.alv.

PiCO.sub.2 -PECO.sub.2 =PE.sup.bag CO.sub.2.

If PaCO.sub.2 is new,

Calculate new deadspace ratio: V.sub.D.sup.alv (alveolar deadspace)=VT.sup.alv (PaCO.sub.2)/(PaCO.sub.2 -PiCO.sub.2). Alvcon (deadspace ratio)=V.sub.D.sup.alv /V.sub.T.sup.alv.

VDphy(physiological deadspace)=V.sub.D.sup.alv /V.sub.T.sup.airway.

Estimated PCO.sub.2 (EPaCO.sub.2)=(PECO.sub.2 -AlvconPiCO.sub.2)/(1-Alvcon).

InvokeExpert: Rules to check if deadspace ratio has changed.

(Parameters needed are: MV, ETCO.sub.2, EPaCO.sub.2, CO.sub.2 prod, slope, DMC, DMR, Tinsp, PEEP, PAP).

Whenever the rules below are triggered three times consecutively, a warning is generated (Note that the parameters are compared with the values obtained from the most recent blood gas results):

If MV increases by 800 ml or more and ETCO.sub.2 does not drop by at least 2 mmHg.

If MV decreases by 800 ml or more and ETCO.sub.2 does not increase by at least 2 mmHg.

If MV increases by 800 ml or more and EPaCO.sub.2 does not drop by at least 2 mmHg.

If MV decreases by 800 ml or more and EPaCO.sub.2 does not increase by at least 2 mmHg.

If MV increases by 800 ml or more and CO.sub.2 prod does not increase by at least 20 ml.

If MVdecreases by 800 ml or more and CO.sub.2 prod does not decrease by at least 20 ml.

If slope changes by more than 0.5 mmHg/ml.

IF DMR changes by more than 10 cmH2O/L/s.

IF DMC changes by more than 10 ml/cmH2O.

If Tinsp changes by more than 0.5 s.

If PAP changes by more than 10 cmH2O.

If previous PEEP is less than 10 cmH2O and changes by more than 5 cmH2O.

If previous PEEP is equal or more than 10 cmH2O and changes by more than 2 cmH2O.

InvokeControl: Rules for ventilator control. (Parameters needed are:

Entered from blood gas results: pH, HCO.sub.3, settings limits (I:E ratio minimum, rate limits, volume limits, peak pressure limit).

From ventilator requests: RR, MV, MAP, IE, VT, SMV, PAP, SRR, SVT, PIF, PEEP, DMC, DMR.

From Calculations: Vdphy, EPaCO.sub.2.)

If new blood test result is available, calculate new CO.sub.2 setpoint:

If (pH<7.36) or (pH>7.44)

CO.sub.2 setpoint=1.6706,HCO.sub.3

Calculate alveolar minute volume needed (AMVneeded);

AMV (Alveolar minute volume(ml))=(MV,1000-RR*VDphy).

AMVneeded=EPaCO.sub.2 , AMV/CO.sub.2 setpoint.

Correct for spontaneous breathing; If (SMV>0) and (RR>SRR)

RRspont (Spontaneous rate)=SRR-RR.

AMVneeded=AMVneeded--(SMV, 1000-RRspont,VDphy).

Control.

Calculate an initial tidal volume needed using current frequency, to produce the required minute volume:

newVT=(AMVneeded/RR)+VDphy.

If PAP>PAP limit,

Repeat

Check the proposed VT & RR to rest if their limits have been exceeded (CheckVT, CheckRR);

Results of these checks are entered into the look-up table. Calculate VT & RR using look-up table 1; Results from look-up table decides whether calculated settings are acceptable.

Until the result from look-up table is either "Implement" or "Impossible". else

Repeat

Check VT & RR if limits exceeded.

Calculate VT & RR using look-up table 2.

Until the result is either "Implement" or "Impossible".

If VT>maximumVT, result is "Not Increase".

If VT<minimumVT, result is "Not Decrease". Otherwise result is "OK".

If RR>maximumRR, result is "Not Increase RR".

If RR<minimumRR, result is "Not Decrease RR".

Otherwise result is "OK".

After doing the checks for VT & RR, the VT & RR results are used with the appropriate look-up table to determine the next calculation.

If result is "Inc V" increase VT by 50 ml

If result is "Drop V", decrease VT by 50 ml.

If VT has changed, calculate new RR for the required minute volume (AMVneeded):

new RR=AMVneeded/(VT-VDphy).

If result is "Inc RR" increase RR by 1 BPM

If result is "Drop RR" decrease RR by 1 BPM

If RR has changed, calculate new VT for the required AMVneeded:

new VT=(AMVneeded/RR)+VDphy.

If result from the look-up tables is "Implement",

If (IE<minimumIE)

If (new RR is bigger or equal to the current setting)

Increase PAF by 10 LPM. (Ensure IE ratio is above the minimum by increasing peak air flow).

Implement the new VT & RR.

If result is "Impossible",

A warning alarm is generated to indicate inability to implement required settings.

It will be apparent from the foregoing that while particular forms of the invention have been illustrated and described, various modifications can be made without departing from the spirit and scope of the invention. Accordingly, it is not intended that the invention be limited, except as by the appended claims.