METHOD AND APPARATUS FOR OPTIMISATION OF MECHANICAL VENTILATION

Abstract: The physiological outcome of a mode of ventilation accomplished by a ventilator, which mode constitutes an alternative to the mode by which the patient is being ventilated is studied through computer simulation. The simulation is performed on the basis of data, which in mathematical terms describes both the mechanics of the patient's respiratory system and mode of ventilation accomplished by the ventilator. The simulation starts from one of two principle points of departure. According to one of them the operator indicates the single or multiple modes of ventilation, which should be analysed through simulation. According to the second one the operator indicates the intended therapeutic goal. Then the computer performs a process of simulation with varying modes of ventilation until a mode has been identified, which leads close to the goal indicated by the operator.

Description: Method and apparatus for optimisation of mechanical ventilation The present invention pertains to a method for optimisation of ventilator function aiming at that desired therapeutic goals should be reached. The optimisation is achieved on the basis of investigated physiological properties of the respiratory system. Mechanical ventilation may cause lung damage. High pressures lead to damage denoted barotrauma. Ventilation at low pressures and volumes may lead to collapse and re-expansion of lung compartments during expirations and inspirations, respectively. This phenomenon may be denoted Recorex. Recorex may lead to lung trauma because of shear forces between zones of collapsed and aerated lung parenchyma. Studies of lung mechanics and adaptation of mechanical ventilation to lung physiology allow reduction of barotrauma and Recorex. The elastic pressure volume diagram, the Pei/V curve, can be recorded, for example with an electronically controlled ventilator. With one or several P e i/V curves recorded under different circumstances one can judge whether the tidal volume, Vt, as a whole or in part falls within the pressure/volume range within which barotrauma and Recorex are minimal. Under guidance of this information one may change the pattern of ventilation in such a way that an optimal P/V range is used to ventilate the lungs. One may for example increase the positive end expiratory pressure, PEEP, if the lung volume is to low. An alternative is to increase the frequency or to reduce the time for expiration in relation to the time for inspiration so that the lungs during expiration do not empty to a degree such that elastic forces exerted by the thoracic cage and the lungs are equilibrated. Thereby, one accomplishes that the alveolar pressure at the end of expiration is positive, which is denoted auto-PEEP. This implies that the lung volume is increased so that Recorex is avoided. If one rather would find that the pressure should be diminished during inspiration, one may decrease ventilation CONFIRMATION COPY and thereby lower the peak airway pressure in order to decrease the risk for barotrauma. Either an increase of PEEP or a lowering of the peak airway pressure frequently lead to a decrease of pulmonary gas exchange which may be deleterious. A lowering of the oxygenation of the blood can be counteracted by an increase in the fraction of oxygen inspired gas. This is, however, associated with risks. Another effect is that CO 2 -elimination is decreased which implies that carbon dioxide is retained in the body. Lately one has frequently accepted this effect, denoted permissive hypoventilation, however, without improved clinical results. Several physiological mechanisms are dependent on pH and thereby on the partial pressure of CO 2 in arterial blood, PaCO 2 . Accordingly it is important, not only to exercise control over airway pressures but also over gas exchange so that PaCO 2 is maintained within suitable limits. By recording of a so called single breath test for CO 2 one may estimate to what extent a change of Vt and minute ventilation, Vmin, will lead to a change of CO 2 elimination expressed in ml/min. By measuring, or estimating, how much CO 2 elimination changes at a change of breathing pattern one can during the following breaths estimate to what extent the change will lead to a change of PaCO 2 . In order to improve the results of mechanical ventilation at grave lung disease both barotrauma and Recorex should be avoided. In order to get around that the gas exchange is unduly affected so that CO 2 retention or hypoxia develops one may better use the Vt. by flushing the connecting tubes, particularly the tracheal tube with unspoiled gas during the later part of expiration. Thereby one may decrease the dead space so that one may decrease Vt without the risk for reduced gas exchange. The pulmonary exchange of carbon dioxide and oxygen are coupled to one another because for each volume of oxygen taken up and consumed through metabolism a thereto proportional volume of carbon dioxide is produced and eliminated. What has been said above about carbon dioxide is accordingly paralleled by a corresponding phenomenon for oxygen. Technically it is more difficult to make fast and accurate determination of oxygen concentration than of carbon dioxide concentration. Particularly, at high inspired oxygen concentrations it is very difficult to accurately measure oxygen uptake. For these reasons carbon dioxide is below focussed upon in discussions about gas exchange although most aspects are relevant also with respect to oxygen. Strategies for avoidance of lung damage following mechanical ventilation, and even to protect the lungs against damage, is denoted protective lung ventilation, PLV. One objective of the present invention concerns new principles for determination of how PLV should be performed on the basis of the physiological characteristics of the individual patient and in consideration of the demands of an adequate gas exchange. The invention is founded on measurements of respiratory mechanics and also of gas exchange. The measurements are followed by a mathematical characterisation of these physiological characteristics. Thereafter a computer performs an analysis of the mode of ventilation by the ventilator with the goal that the therapeutic goal defined by the responsible operator (doctor or therapist) is reached. In a typical case the goal is defined as a combination of PLV and adequate gas exchange. Methods for studies of the mechanics of the respiratory system are previously known, for example through the Swedish patent C2 506521 and reference 1 and 2. Accordingly, one may use a computer-controlled ventilator to study the mechanics of the respiratory system. A pressure and volume range extending beyond the range of the tidal volume at the current setting of the ventilator, Vt, may be studied. A modification of one expiration by prolongation of the time for expiration and/or by reducing PEEP allows a volume range below the Vt to be studied. A volume range above the Vt may be studied by increasing of the volume of gas that is insufflated during the inspiration following the modified expiration. The three segments of the sigmoid P e iA curve are described by an equation in which the elastic recoil pressure is related to the volume, V. The reference volume for V is the lowest observed volume during the measurement. The intermediate segment which begins at the volume Vlip and the pressure Plip is linear and has, accordingly, a constant slope corresponding to a value of compliance denoted Clin. The lower and the upper segments are non-linear. The slopes of these segments (compliance) approach asymptotically zero when the curves are extrapolated towards a low and a high volume, denoted Vmin and Vmax, respectively. The equation that describes Pel as a function of V can be as follows: For V < Vlip: (1) Pel = Plip - (Vlip-Vmin) / Clin * ln((Vmin - Vlip)/(Vmin - V)) For Vlip ≤ V ≤ Vuip: Pel = Plip + (V - Vlip) / Clin For V > Vuip: Pel = Puip + (Vmax-Vuip) / Clin * ln((Vmax - Vuip) / (Vmax -V)) The pressure that drives flow through the airways can in analogy with Ohm's law be calculated as resistance multiplied by flow rate. Inspiratory resistance, R(i), often varies to a small extent during the breath. It may decrease with increasing lung volume. The pressure that drives flow during insipiration Pres(ϊ), can defined as a function of flow rate, F, and volume, V: Presrø = (rO + M * V) * F (2) rO denotes the resistance at zero volume, while r1 , usually a negative number, gives the variation of R(i) with volume. Equation 1 and 2 describe the mechanical behaviour of the respiratory system during inspiration under various circumstances. Their coefficients can be determined with known methods. However, simpler equations may often be used. Within a limited volume range a linear equation may with adequate accuracy describe the P^IM curve. The inspiratory resistance can often be represented by a constant. This invention can be embodied on the basis of other equations, which in the actual situation adequately represent the elastic and resistive properties of the respiratory system. Furthermore, the invention is not dependent upon which method is used for measurement and analysis leading to a mathematical description of the elastic and resistive properties of the respiratory system. With the intent to avoid the limitations of previously known methods for establishment of PLV method and apparatus have been given the characterisation according to claims 1 and 25, respectively. Preferred embodiments of the invention will below be described with reference to figures 1-2, which represent printouts of the computer screen. According to a preferred embodiment of the invention the function of the respiratory system i studied with a known system based upon a computer- controlled ventilator complemented with a CO 2 analyser, (CO 2 analyzer 930, Siemens-Elema, Solna). During a sequence of breaths analogue signals corresponding to flow rate, pressure and CO 2 concentration are converted to a digital format and stored. The flow signal is corrected for compression/decomplression of the gas in the tubings and is then integrated to a volume signal. The reference level of the latter is the minimum volume during the sequence of breaths. From the pressure measured in the ventilator the pressure in the trachea, Ptr, is calculated by subtraction of the pressure driving flow through connecting tubing. The latter pressure is calculated as the product between flow rate and tube resistance. The mechanical properties of the respiratory system can be studied and defined as described above with known methods. In order to fulfil make possible a computer analysis of different modes of ventilation the mechanics during expiration must also be mathematically defined. This is done as described below. The tracheal pressure has during mechanical ventilation of the sedated patient that has no own muscular activity two components, the elastic recoil pressure, Pel and the the resistive pressure compnent, Pres: Ptr = Pres + Pel (3) According to a preferred embodiment measurements are performed during insufflation of a modulated sinusoidal flow pattern. Utilising a known iterative technique the parameters of equation 1 and 2 are calculated. Thereby the elastic and resistive properties of the respiratory system under inspiration become known. It is known from studies at health and in disease that the P^IM curve displays no significant hysteresis over ordinary or modestly increased tidal volumes. Accordingly, Eq. 1 can often be used to describe Pel as a function of volume both for inspirations and expirations. Resistance differs often between inspiration and expiration. Expiratory resistance, R(e) can be determined according to the following: For each sampling interval during one or more expirations R(e) is calculated first by calculation of Pres from Eq. 3, and then by dividing Pres by flow rate, F. F and Ptr are measured, while Pel is calculated from volume according to Eq. 1. The expiratory conductance G(e) is then calculated as 1/Rfe . G(e) is approximately linearly correlated to volume : Gte) = g0 + g1 * V (4) In order to avoid the influence by fast accelerations and decelerations at the flow transitions between inspiration and expiration the calculation of the coefficients gO and g1 are based on a volume interval which does not include the start and end of the expiration that is studied. In order to provide data for gas exchange the single breath test for CO 2 is analysed in a special format. The volume of CO 2 , which is expired during the course of a breath, Vt C o 2 is calculated from the integral of the product between flow rate and CO 2 concentration. Vt C o2 is related to expired volume, Fig. 1 , panel 1. The upper part of the curve, indicated in heavy line, can with a high precision be described as a 2 nd degree equation: Vt C 02 = a C o2+ co2- Vt + Cco2- Vt 2 (5) Equations 1-5 constitute the mathematical description of the respiratory physiology from which analyses and simulation of lung function is performed. The coefficients of the equations, presented in panel 3, Fig. 1 and 2, are automatically transferred from the Excel work sheet in which they are calculated to the work sheet in which further analysis is made. The description of the physiological properties of the respiratory system can be expressed in equations which are alternatives to eq.1-5. The mathematical description of respiratory mechanics and, according to a preferred embodiment also gas exchange, serve according to the invention as a basis for simulation of the course of physiological events that would follow a change in mode of ventilation. The description of the physiological properties of the respiratory system may be in the format of tables as well as equations. The methods employed for the description of the respiratory physiology may vary as long as they fulfil the demand the simulation process raises. Simulations of mechanics and gas exchange, as they would appear after the establishment of an alternative mode of ventilation, can with the present system be performed from two different starting points. 1. The operator defines an alternative mode of ventilation, which he wants to analyse with respect to mechanical consequences and gas exchange. The computer then performs the analysis and presents the results as shown in Fig. 1. Then the operator considers if the analysed mode of ventilation accords with the immediate therapeutic goal he has defined. He could also analyse further modes in order to find one that is deemed preferable with respect to the therapeutic goal. 2. The operator defines the therapeutic goal. The computer then performs an iterative search for a mode of ventilation, which leads as close as possible to the therapeutic goal. In this search limits of the parameters defining the mode of ventilation defined by the operator are respected. Independent upon which starting points the invention is practised, the mathematical simulation of the events following a change of mode of ventilation follows the same principle. For each tested mode of ventilation six breaths are mathematically simulated in order to provide a steady state with respect to mechanical events, notably the volume and the pressure from which a breath starts. Simulation can be done for volume or pressure controlled ventilation. In this context volume controlled ventilation is further described. The analysis starts from a mathematical description of the alternative mode of ventilation in terms of a number of parameters depicted in panel 3, Fig. 1. 1. Minute volume, MIN_VOL. 2. Breathing frequency, Frekv. MIN_VOL and Frekv define together tidal volume, Vt: Vt = MIN_VOL/Frekv 3. Inspiratory time in % of the breathing cycle, Ti%. 4. Flow ramp during inspiration, RAMP. 0 implies a constant inspiratory flow rate. 100 implies that flow rate increases from zero to twice mean inspiratory flow. -100 implies that flow decreases from twice mean inspiratory flow to zero. 5. The pressure in the ventilator at start of expiration, Pe st art . 6. The pressure in the ventilator at end of expiration, Pe end - The simulation is performed by partitioning of each simulated breath into short time intervals. Most intervals are 1 % of he breathing cycle. At the start of each phase of the breathing cycle the intervals are 1 0/00 of the breathing cycle in order to increase the precision of the simulation of transient events. During each interval, with the duration Δt, the process is simulated by a series of calculations, which are based upon the properties of the respiratory system and the function of the ventilator during the actual time interval or the immediately preceding interval. The resulting phase shift between different steps in the calculation, which may correspond to one time interval, does not lead to significant inconveniences as the intervals are short, in particular those around the transitions between inspirations and expirations. For each time interval during the inspiration the flow rate from the ventilator, F, is first calculated from the frequency, Vt, Ti% and RAMP, as well as from the time of the interval in relation to the start of the inspiration. According to an alternative embodiment the flow pattern during inspiration is not simulated as a linear ramp but in accordance with an elective equation or a table. For example, the flow rate may initially during the inspiration increase linearly to reach an early maximum and then be given the form of an exponential decline. Independent upon which flow pattern is applied, F is corrected for the volume that does not leave the tubing in order to maintain the tracheal flow rate equal to Ftr in accordance to previously known principles. The volume, V, is during the complete simulation process by consecutively summing the product Ftr*Δt from the moment and from the volume of the lungs at the beginning of the measurement procedure of lung mechanics. Then, Pel is calculated from V and the parameters defining the P^/V curve according to Eq. 1. Pres(i) is calculated according to Eq. 2. Ptr is then calculated from Eq. 3. The pressure in the ventilator, Pvent, is calculated by adding to Ptr the pressure drop from the Y-piece to the trachea calculated by multiplication Ftr with the resistance of the tube at the actual value of Ftr. The resistance of the tube, Rtube, at a certain flow is calculated according to Rohrer's formula from the coefficients K_1 and K_2 shown in panel 5 in FIG. 1 and 2. During simulated expirations and for each time interval the pressure controlled by the ventilator, Pe, is calculated. Pe is according to the preferred embodiment allowed to vary as a linear ramp from Pe s tart to Pβend- Then, the expiratory flow is calculated as: 10 ((Pel-Pe)/R(e)+Rtube). Then, the volume is calculated as has been described for inspiration and also Pel according to Eq. 1. A new value for G(e) is calculated according to Eq. 4. G(e) is inverted to obtain R(e). For each time interval Rtube and Ptr are then calculated as has been described for the inspiratory phase. After the simulation of one breath an additional number of breaths are simulated until the system comprising the ventilator and the respiratory system have come to a steady state with respect to pressure, flow and volume during succeeding breaths. The events during the last simulated breath constitute the result of the simulation. The result is presented in the format of curves over flow 6, of the pressure in the ventilator 7, and of Ptr 8, all plotted against time during the breath. Another presentation 9 shows Presp, Ptr and Pel plotted as a loop against volume changes during the breath. From the column Alternative in a panel 10 the following values constituting the outcome of the simulation are shown: Tidal volume Vt, peak and mean Ptr, peak and mean Pel and PEEPtot which is the pressure in the alveoli corresponding to Pel at the end of expiration. One may further read the value of PaCO 2 , which is estimated to result at a steady state as described below. A mode of ventilation that is considered to constitute an alternative to the one applied at the time t1 is below denoted alternative mode. Immediately after the institution of an alternative mode, at the time t2 Vmin ∞ 2 will change from the value Vmin C02 ( , ) to the value \Jm\n C o 2 (t2). At a new steady state after 20 to 40 minutes, at the time t3 PaC0 2 (t1) will have changed to PaC0 (t3). The change of PaCO 2 within the period t1 to t3 is according to known principles indirectly proportional to the change of Vmin C o 2 from t1 to t2. The values of Vt at t1 and t2 are entered into equation 5. VtCO 2 can be calculated at these moments and then by multiplication with the frequency F at t1 and t2 also \fm\r\ C0 2(t1) and Vmin C02 (?2) are calculated. Then, PaC0 2 (t3) is calculated: PaCO 2 ff3J = PaCO 2 ^j * Vmin C0 2(» ) I Vmin C o2(?2) (6) 11 If PaC0 2 (t1) is unknown its value may be set to 1. Then, Eq. 6 will give PaC0 2 (t3) in relationship to PaC0 2 (t1). The simulation of mechanics, and according to a preferred embodiment also of gas exchange, may according to the invention and according to the principle for calculation described above be performed from two perspectives. From one perspective the operator determines such alternative modes of ventilation that he intends to be studied. The simulation by the computer is followed by the operator's judgement if the alternative mode of ventilation should be implemented or if further simulations should test other alternative modes. The computer simulation according to the invention is performed with a program, which according to one embodiment is realised as a sheet for calculation of type Excel. A corresponding function may, however, be obtained with most programming languages. Figure 1 depicts the computer screen as it appears after a simulation. The content of the screen can easily be varied in order to be in harmony with the value an operator attributes to various partitions of the total information. According to the example of Fig. 1 the operator has applied the perspective that he determines the mode of ventilation that he desires to analyse through simulation. He introduces values describing an alternative mode of ventilation in the column Alternative within the panel Mode of Ventilation 4, which mode he wants to analyse through simulation. Thereby he brings the program to analyse that mode. From Fig. 1 it appears that the operator, as a comparison to the actual mode of ventilation, desired to analyse the outcome of a change of minute ventilation from 11 to 14 and a change of the frequency from 20 to 30. From the panel OUTCOME 10 he can among other observations note that the alternative mode is predicted to increase PEEPtot from 4.6 to 6.6 and to decrease PaCO 2 from 40 to 36. If the operator considers that this result corresponds to the therapeutic goals he intends to reach with an alternative mode of ventilation he may initiate such a mode. He might alternatively choose to enter further alternative in panel 4 and thereby 12 simulate any number of alternative modes until he finds a mode, which with respect to expected result corresponds to his expectations. Simulation of the events at pressure controlled ventilation is performed according to the same principles as at volume controlled ventilation. The difference is that for each time interval during inspiration the pressure generated by the ventilator, Presp, is first calculated from an equation describing how that pressure varies during inspiration. Presp minus Pel gives the pressure driving flow during the time interval. By dividing this pressure by the resistance in the tracheal tube and airways the inspiratory flow rate is obtained. The integral of the flow rate gives the volume change during the interval. From the volume Pel is then calculated according to eq. 1. Ptr is then calculated followed by calculation of a new value of Rtub. During expiration the simulation is performed as for volume controlled ventilation. Going out from the second of above mentioned starting points simulation is performed within and iterative procedure serving to find that mode of ventilation, which as closely as possible corresponds to a predefined therapeutic goal. This goal may for a certain category of patients be defined by a programmed so called default. According to a preferred embodiment it is, however, the operator who defines the desired therapeutic GOAL. The GOAL may be of purely mechanical nature such as to achieve the lowest possible pressure in the patient's airways or alveoli. The GOAL may alternatively be a combination of results with respect to both mechanics and gas exchange. The importance of each component described by the GOAL may be weighted according to principles described below. The lung is for gas exchange. Disregarding that an operator according to the preferred embodiment may choose and weigh different components comprising the GOAL, an embodiment based upon that the gas exchange should be assured to an extent defined by the operator is further described. From this prerequisite the simulation process seeks a mode in order to reach other components of the GOAL. In a patient ventilated by a ventilator the alveolar 13 pressure at the end of expiration corresponds to Pel. This pressure is denoted total PEEP, PEEPtot. In order to maintain PLV PEEPtot should not fall below the pressure at which lung collapse occurs. Therefore a common component of the GOAL is to maintain a defined PEEPtot. Another component of the GOAL is that tracheal or alveolar peak pressure should be as low as possible in order to avoid barotrauma. A combination that maintains a defined PEEPtot and minimises peak pressure is important for PLV. In patients with severely compromised circulation it may, however, be important to decrease the mean tracheal or alveolar mean pressure. The preferred embodiment gives the operator a large degree of freedom to define the GOAL in terms of several components. Definition of the GOAL and the following simulation may follow the example illustrated in FIG. 2 showing a printout of the computer screen after a complete process of simulation. The operator has defined the GOAL by filling in the column GOAL in panel 11 and by defining the weight of each partial goal in panel 12. Panel 11 shows that the operator according to the example desires to let PaCO 2 rise from 40 to 50 mm Hg. Furthermore, he wants that PEEPtot should increase from 4.6 to 10. Panel 12 shows that the operator attributed the partial goal to reach the desired PEEPtot a weight of 2 and to minimise peak tracheal pressure a weight of 1. Remaining pressures in panel 12 have been attributing the weight of zero. This implies that the values of these pressures do not influence the process of simulation. The significance of the concept weight is further explained below. When the operator according to the preferred embodiment has noted that a partial goal is to achieve a certain PaCO 2 the degrees of freedom of the iterative process is limited to such alternative which satisfies the condition that the given PaCO 2 is achieved. The program attributes accordingly the given value PaCO 2 a weight that in principle is infinite. This effect is achieved in the following way. The value of PaCO 2 , defined as goal by the operator corresponds to PaC0 2 (t3) in equation 6. In the same equation PaC0 2 (t1) and \Zm C o 2 (t1) correspond to actual known values. Vmin C o 2 ^2 is the only unknown entity in the equation and may accordingly be calculated. In the iterative process alternative 14 values of Vt are tested in the search for an alternative mode of ventilation. The examination of a certain value of Vt starts with a calculation of VtCO 2 at the moment (t2) according to equation 5. The Frequency is calculated as the quotient \Jm\x\coz(t2)Ntco 2 (t2). The minute ventilation is then calculated as WFrequency. In the described manner the given value of PaCO 2 is at a tested value of Vt coupled to given values of Frequency and minute ventilation. This implies that only such combinations of Vt, Frequency and minute ventilation, which satisfies the demand that PaCO 2 , after initiation of an alternative mode and the establishment of a steady state, will reach the value PaCO 2 , G oAL- The values of Frequency and minute ventilation, which are calculated from the target PaCO 2 and the value of Vt found by simulation, are presented in panel 13. The analytical method used for calculation of Frequency and minute volume from an alternative value of Vt implies that the number of parameters describing the mode of ventilation analysed in the iterative process is reduced by one. This gives advantages in the form of less time consumption and improved security with respect to convergence of the iterative procedure toward a solution satisfying the goal determined by the operator. Similar results may in principle be achieved according to other embodiments allowing the iterative process to vary independently two of the factors Vt, Frequency and minute ventilation. Convergence of the iterative process toward a goal implying that a certain PaCO 2 should be reached may be achieved by entering PaCO 2 among other partial goals in panel 12 and by attributing this partial goal a high weight. Apart from defining the goal as described the operator may enter the parameters describing the mode of ventilation, which are allowed to be changed. There are furthermore reasons for limitation of the change of the parameters. One reason is that the validity of description of the function of the respiratory system according to panel 3 is greater within a limited range of modes of ventilation than within a larger range. Another reason may be that the ventilator used may have technical limitation with respect to variability of the mode of 15 ventilation. A third reason may be that tradition and experience are restricted to a limited range of mode of ventilation. For example, a negative experience from a particular category of patients may merit limitation of the allowed mode of ventilation by entering limits for changes of each parameter describing the mode of ventilation. It should also be mentioned that a limitation of allowed changes of the parameters is further merited by the existence of so called "local minima of the error signal" Such minima are particularly present at extreme combinations of the parameters which are allowed to vary. Risks associated with convergence of the numerical process toward a non-relevant minimum may be practically eliminated by limiting the range of allowed variation of the parameters describing the mode of ventilation. According to the embodiment depicted in Fig. 2, parameters describing the mode of ventilation are described in panel 14 together with the limits within which the parameters are allowed to vary. The operator may change these limits presented as their default values. Accordingly, he needs only to enter aberrations from the default values. The system based upon default may according to the invention be used to facilitate for the operator to use the system, by calling for menus specific for particular groups of patients and situations to be presented on the computer screen. One may for example call upon a menu intended for ARDS patients. In this menu a partial goal might be a value of PEEPtot calculated from the coefficients in panel 3 describing the Pel/V curve of the patient. This goal may be attributed a high weight. One may follow a tradition aiming at limitation of peak Ptr, Ptrp eak . Starting from a menu of this character, which is intended for ARDS patients, the operator only needs to determine if PaCO 2 should be maintained at the current value or if it should be changed. Accordingly, by calling upon a menu for ARDS patients, which is defined in the treatment protocol of the establishment, it is possible even for a modestly experienced operator to make use of the system. 16 If the patient suffers from severe obstructive lung disease the goal may be to reduce PaCO 2 while, on the same time one seeks a setting, which gives the lowest possible pressure in the alveoli at the end of expiration, i.e. PEEPtot. In that particular situation a menu may be called upon in which by default values such goals are predefined and weighted. In a preferred embodiment of the invention default values for different groups of patients are combined with freedom for the trained operator to modify default values in accordance with his experience and knowledge. When the operator has defined the GOAL and the limits for alternative modes of ventilation he starts the iterative estimation of such a mode that in an optimal way leads towards the GOAL. Within the iterative process such a combination of parameters describing the mode of ventilation, which leads to a minimum of a parameter denoted ERROR, is sought. ERROR is calculated from deviations of those partial goals, which the operator has defined or accepted as default. The outcome is read from the last breath of the simulation. ERROR is the sum comprising components from each partial goal. For a partial goal, which reflects that a certain value of a parameter describing the outcome should be reached, a primary error is calculated as the deviation between the defined partial goal in panel 12, which is denoted Ideal, and the value of the parameter, which represents the outcome of the simulation, denoted Simulated in panel 12. This primary error is attributed a positive sign either by taking the absolute value of the error or by raising the error to the second power. According to a preferred embodiment, in a situation in which PEEPtot should attain a certain value the error of PEEPtot is calculated by raising the deviation to the second power before the resulting error is multiplied by the factor WEIGHT. Thereby one achieves that small deviations contribute little to ERROR while large deviations contribute much. This is advantageous within the iterative process as this process more efficiently converges towards an optimal solution. When the operator wants that the target PEEPtot should be 17 accomplished with a high precision he attributes the partial goal PEEPtot a relatively high weight according to the example in panel 12, Fig. 2. In order to achieve that a partial goal in panel 12, which reflects that a parameter describing the outcome should be brought to a minimum, the primary error is calculated as the deviation between the defined partial goal and the simulated value of the parameter. Thereafter, this primary error is multiplied by the weight the operator has attributed the error. The contribution by each partial error to the total error is presented in the column ERROR in panel 12. This is also the case with respect to their sum, which as described is denoted ERROR. As is shown by the embodiment illustrated in Fig. 2 panel 12, the partial goal with respect to PEEPtot may be to maintain a certain target value or to minimise the value. Other values for which minimal values may be sought are peak and mean values of Ptr and Pel. The operator may as shown be free to choose a combination of parameters, which together define the GOAL. He does so by attributing the error for different parameters a weight greater than zero. In order to make different errors equally important and to make them independent upon which denomination they are expressed in, it is according to the invention possible to normalise the primary error for each parameter in such a way that they in a typical situation will contribute to the error to the same extent if the operator has attributed them the same weight. When, as described, the operator has defined the partial therapeutic goals as well as the limits for alternative modes of ventilation he starts an iterative simulation process. In a numerical process the computer seeks the mode of ventilation, which leads to the lowest value of ERROR. Several methods for iterative estimation of parameters are known and commonly used. According to experience, the iterative estimation of parameters describing the mode of ventilation, which minimises ERROR, is in the preferred embodiment based upon a principle commonly denoted Newton's. The principle is applied by successive approximations based upon linear tangents in the multidimensional sphere of data. Computer tools for such estimations are available in several commercially 18 available programs. In a preferred embodiment the tool Solver (EXCEL 5.0 Microsoft) was used. The results of the simulation are presented for the operator both in terms of parameters describing the mode of ventilation, which according to the simulation leads as closely as possible to the by the operator defined GOAL. These parameters are shown in Fig. 2, panel 14, column Alternative. According to the preferred embodiment also the predicted course of a breath is displayed in diagrams over time showing flow rate 15, Pvent 16 and Ptr 17. Furthermore, a diagram 18 is shown, depicting the pressure in the ventilator, Ptr and Pel presented as loops against the volume V during the breath. The operator analyses the outcome of the simulation and considers if the outcome with respect to security and the medical reasons is such that a change of mode of ventilation in accordance with the simulation can be expected to lead to the defined therapeutic goals. The operator further judges if this mode of ventilation can be expected to benefit the patient from a comprehensive perspective. If that is the case he may initiate a change of the mode of ventilation in accordance with the simulation. If he finds that further alternative need to be explored he may after modification of the GOAL of the limits for parameters describing the mode of ventilation repeat the simulation and on the basis of the new simulation again consider if an alternative mode of ventilation can be expected to benefit the patient. The operator may also manually modify the parameters describing the mode of ventilation according to the simulation. Such modification is automatically followed by a new simulation and presentation of the outcome of such a modified mode of ventilation. The operator may again judge if the latter outcome can be used as a basis for an alternative mode of ventilation. If so, he may initiate that mode. At pressure controlled ventilation Vt is a result of the course of Presp during the breathing cycle. Accordingly, Vt is unknown at the start of a simulation process. At simulation of pressure controlled ventilation with the goal to maintain a certain value of PaCO 2 , this cannot be reached by limiting the iterative process to such modes of ventilation, which leads to the desired value of PaCO 2 . In this 19 situation PaCO 2 is dealt with according to principals corresponding to those above described for PEEPtot according to the preferred embodiment. This implies that the difference between desired value for PaCO 2 and the value obtained through simulation is raised to the 2 nd power and is attributed a relatively high weight. Thereby, to a corresponding degree, all deviations will give an important positive contribution to ERROR. According to a preferred embodiment the computer, which is used for studies of the respiratory system physiology, is also used for the simulation of alternative modes of ventilation. The system used for such purposes has been described 1,2 . An advantage of this embodiment is that it takes no system for transfer of data describing the physiology of the respiratory system from a measurement system to the computer used for simulations. Other advantages may, however, be reached through alternative embodiments. The simulation may be independent upon the measurement system, which another apparatus and at another place. After transmission of data describing the situation of the patient and his physiology experts in respiratory care may perform simulation and on the basis of these contribute to improved respiratory care without being present at bedside. The characterisation of the mechanics of the respiratory system as described above may cover a volume and elastic pressure range far above and below the volume range of the current tidal volume. As an alternative mode of ventilation may reach into these zones this may be an advantage, as the physiology over the alternative volume range is then known. However, in some conditions like ARDS, the physiology may immediately change from that prevailing during normal tidal volume when a wider range of volume is studied 3 . The elastic pressure volume curve will change as a result of a reduction of lung volume at which airway closure and alveolar collapse occurs. Under such conditions, the P e i/V curve studied with a starting point at low lung volumes will not be valid during tidal ventilation at higher lung volumes. An alternative embodiment of the invention has been developed in order to avoid these problems. In this 20 alternative embodiment the physiology of the respiratory system is studied during ventilation with volumes which do not cover such a volume range that the physiology is significantly changed by the measurement procedure. Preferably, in the alternative embodiment respiratory physiology is studied just within the current tidal volume range. In the alternative embodiment flow rate, airway pressure and carbon dioxide the computer records signals for a certain number of consecutive breaths. For at least one breath the pressure following a postexpiratory pause and a postinspiratory pause is recorded. The analysis of the sequence of breaths, most of which may be ordinary breaths, starts with a computer identification of the pauses mentioned and reading of the postinspiratory pressure, denoted plateau pressure, Pplat, and the postexpiratory alveolar pressure denoted PEEPtot. The elastic pressure volume curve is considered to be linear within the limited volume range corresponding to the tidal volume. This is known to be a reasonable approximation. The course of the tracheal pressure during the breaths is calculated as described above. The conductance of the respiratory system, G, is for each moment calculated as: G = Flow/(Ptr-Pel) (7) Conductance and its variation with lung volume during inspiration and during expiration are calculated as separate linear equations. The perturbed data at transitions between inspiration and expiration are not included in this calculation. The exchange of CO 2 is analysed as described above. The data describing the physiology of the ordinary breath are transferred to the simulation program. An advantage of a more simple mathematical description of the physiology of the respiratory system is that the simulation process may be simplified. The use of continuous and often linear equations makes simulations possible with analytic rather than by numerical techniques. When numerical techniques nevertheless are employed the iterative simulations more easily converges toward a true minimum. An other advantage of the alternative embodiment is that the measurement procedure is not associated with any significant deviations from the current breathing pattern. This eliminates potential 21 hazards associated with particularly low and high airway pressures. The measurement procedure can be repeated at frequent intervals and be completely automated. The validity of data is increased with respect to the current volume and pressure range. However, extrapolation outside this range is associated by potential errors. Accordingly, the initiation of modes of ventilation, which are very different from the current mode based upon simulations, may lead to unpredictable results. Repeated measurements and simulation procedures as further described below may, however, compensate for this potential drawback. The simulation program of the alternative embodiment based upon tidal volume analysis follows the principles described above. However, after execution of the simulation according to either of the two starting points, the initiation of an alternative mode of ventilation may if so decided by the operator be realised in a stepwise manner. The operator may enter the number of steps and the time interval between the steps by which the alternative mode should be realised. Then, the computer suggests a new mode of ventilation by interpolation in between the parameters describing the current ventilator setting and the alternative obtained through simulation. In a preferred embodiment each step is preceded by a new study of the respiratory physiology and a repeated simulation procedure. The operator may initiate and define a process through which the alternative mode is introduced stepwise in a manual, or alternatively in fully or partially automated procedure. A typical setting would be that the initial measurement and simulation procedure is followed by a change in mode of ventilation such that for each parameter describing the mode the gap between current setting and the setting suggested from simulations is reduced by for example 25 or 50%. After an interval defined by the operator the measurement and simulation procedure is repeated. Thereby a successive approach toward a mode of ventilation that in an optimal way leads to the goal defined by the operator is achieved. The computer is equipped with means for controlling the function of the ventilator as is previously described 1,2 . By emitting signals to the 22 ventilator the computer can override the settings on the front panel of the ventilator. These means makes it possible to let the computer effectuate changes of mode of ventilation within wide technical limits. However, as described above the operator defines the limits within which he allows changes of parameters describing the mode of ventilation. The procedure may according to the operator's decision be automated with or without operator surveillance. 23 References: 1. Respiratory mechanics during mechanical ventilation in health and in disease. Cecilia Svantesson, These 1997, Department of Clinical Physiology, Lund University, ISBN 91-682-2766-9 2. A single computer-controlled mechanical insufflation allows determination of the pressure-volume relationship of the respiratory system. Svantesson, C, Drefeldt, B., Sigurdsson, S., Larsson A., Brochard, L. and Jonson, B. J Clinical Monitoring and Computing 00: In press 1999.

Claims: 24 Claims 1. Method for optimisation of the mode of ventilation by a ventilator, characterised in that the course with respect to flow rate in the airway, the lung volume changes and/or the pressure at one or more levels in the airways and the lungs are calculated through computer simulation on the basis of in mathematical terms expressed properties of the elastic and resistive properties of the respiratory system and from in mathematical terms expressed modes of ventilation, which modes of ventilation may constitute alternative to the current mode, whereby the course, which according to the simulation is predicted to follow an initiation of the alternative mode of ventilation, is calculated to serve as a basis of analysis of the alternative mode of ventilation with respect to goals which according to the intentions should be reached through treatment with a ventilator. 2. Method according to claim 1 characterised in that the calculation and analysis of the course during one or more consecutive breaths apart from flow rate, volume or pressure also comprises the exchange of carbon dioxide or oxygen, which calculation is performed on the basis of in mathematical terms expressed relationship between the extent of gas exchange and parameters describing the course of lung ventilation, which results from such mode of ventilation, which may constitute alternative to the current mode of ventilation 3. Method according to claim 1 or 2 characterised in that calculation through computer simulation of the course at an alternative mode of ventilation is performed with numerical technique. 4. Method according to claim 1 or 2 characterised in that the operator indicates the single or those multiple alternative modes of ventilation, which through computer simulation should be analysed with respect to the resulting course, after which the operator performs the analysis of whether some 25 alternative mode of ventilation leads to the goal, which according to the intentions should be achieved. 5. Method according to claim 1 or 2 characterised in that the computer through simulations performs simulations of such modes of ventilation, which are programmed as defaults, after which the operator performs the analysis of whether some alternative leads to the goal, which according to the intentions should be achieved. 6. Method according to claim 1 or 2 cha racterised in that the operator indicates the goal , which is intended to be achieved through an alternative mode of ventilation, after which the computer performs a mathematical analysis in which the in mathematical terms expressed mode of ventilation is varied so as to identify through simulation, which mode leads close to the indicated goal. 7. Method according to claim 1 or 2 and some combination of claim 4, 5 and 6 characterised in that simulated modes of ventilation comprise a combination of such modes of ventilation which are analysed in an iterative numerical process, such modes which have been programmed as defaults and such modes which the operator has indicated. 8. Method according to claim 6 or 7 characterised in that alternative modes of ventilation, which are analysed in a mathematical process, are limited with respect to which parameters describing the mode of ventilation are allowed to be varied to achieve the indicated goal and/or with respect to the range within which these parameters are allowed to vary. 9. Method according to claim 2 characterised in that the characterisation of the extent of gas exchange in the lungs is achieved through an equation describing the volume of exchanged carbon dioxide or oxygen in relation to the tidal volume. 10. Method according to claim 2 characterised in that the characterisation of the gas exchange in the lungs is achieved through an equation describing the dead space in relation to the tidal volume. 26 11. Method according to claim 1 or 2 cha racterised in that the characterisation of the elastic and resistive properties of the respiratory system is achieved through equations describing the pressure gradient conditioned by resistance within an airway segment as a function of flow rate and/or lung volume, and the pressure conditioned by the elastic recoil of the respiratory system as a function of lung volume. 12. Method according to any of claims 1-11 characterised in that separate equations are used to express the mechanical properties of the respiratory system during inspiration and expiration. 13. Method according to claim 1 or 2 characterised in that the elastic and resistive properties of the respiratory system is expressed in tables, which indicate these properties in relation to changing lung volume or other parameter which bears relation to lung volume. 14. Method according to claim 1 characterised in that alternative modes of ventilation are defined through equations describing pressure or flow rate achieved by the ventilator during the respiratory cycle. 15. Method according to claim 1 characterised in that alternative modes of ventilation are defined through tables describing pressure of flow rate achieved by the ventilator during the respiratory cycle. 16. Method according to any of claim 1-15 characterised in that a goal, which is intended to be achieved through treatment with a ventilator is to achieve a defined level of pressure during the full or partial phase within the respiratory cycle at some level within airways or lungs. 17. Method according to any of claim 1-16 characterised in that a goal, which according to the intentions of an operator should reached through treatment with a ventilator, is to achieve a defined extent of gas exchange. 18. Method according to any of claim 1-17 characterised in that a goal, which is intended to be reached through treatment with a ventilator is to achieve a 27 mathematically defined combination of pressure levels during full or indicated phases within one or more levels of airways and lungs. 19. Method according to any of claim 1-18 characterised in that a goal, which is intended to be reached through treatment with a ventilator is to achieve a mathematically defined combination of pressure levels and extent of gas exchange. 20. Method according to any of claim 1-19 characterised in that the alternative mode of ventilation, which according to the simulation process leads close to the goal, which is intended to be reached through treatment with a ventilator, is approached stepwise by changing the mode of ventilation in steps in the direction indicated by the simulation process. 21. Method according to claim 20 cha racterised in that in between the stepwise changes of the mode of ventilation the physiological properties of the respiratory system are studied and a new simulation process is performed. 22. Method according to claim 20 or 21 characterised in that the computer executes stepwise changes in a process supervised by an operator. 23. Method according to claim 20 or 21 characterised in that the computer executes the stepwise changes in an automated process. 24. Method according to any of claims 1-20 characterised in that the computer reads the physiological result of a change in mode of ventilation and issues a warning if the results do not indicate that the goal is being approached. 25. Apparatus for optimisation of the mode of ventilation comprising a ventilator, sensors for flow rate and pressure and a computer, which is programmed so as to on the basis from in mathematical terms expressed properties of elastic and resistive properties of the respiratory and from in mathematical terms expressed modes of ventilation of the ventilator, which modes may constitute alternative to the current mode, calculate the course through one or more consecutive breaths with respect to airway flow rate, lung volume change and/or pressure at one or more levels in airways and lungs. 28 26. Apparatus according to claim 25 cha racterised in that sensors for oxygen or carbons dioxide connected to the computer, which on the basis of in mathematical terms expressed relationship between the extent of gas exchange and parameters describing the course relating to lung ventilation, is programmed also to calculate through simulation the extent of gas exchange, which results from such mode of ventilation, which may constitute an alternative to the current mode of ventilation. 27. Apparatus according to claim 25 characterised in that the computer program used for simulations constitutes a calculus sheet. (Swedish: Kalkylblad) 28. Apparatus according to claim 25 characterised in that the computer is equipped with means so as to receive indications from an operator, which indications define which single or multiple alternative modes of ventilation that should be simulated. 29. Apparatus according to claim 25 cha racterised in that the computer is programmed so as to by mathematical analysis search for such a mode of ventilation, which leads near to the by the operator indicated goal of treatment with a ventilator. 30. Apparatus according to claim 29 characterised in that the computer allows mathematical analysis with iterative technique. 31. Apparatus according to any of claims 25-30 characterised in that the computer is equipped with means sending signals to the ventilator, which signals influences the ventilator so that it delivers a mode of ventilation, which is in conformity with a mode of ventilation that according to the simulation would lead toward the goal of treatment indicated by the operator. 32. Apparatus according to claim 31 characterised in that the alternative mode of ventilation, which according to the simulation process leads close to the goal, which is intended to be achieved through treatment with a ventilator, is approached stepwise by changing the mode of ventilation in steps in the direction indicated by the simulation process.

Inventors: JONSON BJOEM, JONSON, BJOEM
Applicants: INNOTEK AB, JONSON BJOEM, JONSON, BJOEM
Application year: 1999
Main classification: A61M16
IPC codes: A61M16