 Research
 Open Access
 Published:
On the practical identifiability of a twoparameter model of pulmonary gas exchange
BioMedical Engineering OnLine volume 14, Article number: 82 (2015)
Abstract
Background
Successful application of mechanical ventilation as a lifesaving therapy implies appropriate ventilator settings. Decision making is based on clinicians’ knowledge, but can be enhanced by mathematical models that determine the individual patient state by calculating parameters that are not directly measurable. Evaluation of models may support the clinician to reach a defined treatment goal. Bedside applicability of mathematical models for decision support requires a robust identification of the model parameters with a minimum of measuring effort. The influence of appropriate data selection on the identification of a twoparameter model of pulmonary gas exchange was analyzed.
Methods
The model considers a shunt as well as ventilationperfusionmismatch to simulate a variety of pathologic pulmonary gas exchange states, i.e. different severities of pulmonary impairment. Synthetic patient data were generated by model simulation. To incorporate more realistic effects of measurement errors, the simulated data were corrupted with additive noise. In addition, real patient data retrieved from a patient data management system were used retrospectively to confirm the obtained findings. The model was identified to a wide range of different FiO _{2} settings. Just one single measurement was used for parameter identification. Subsequently prediction performance was obtained by comparing the identified model predicted oxygen level in arterial blood either to exact data taken from simulations or patients measurements.
Results
Structural identifiability of the model using one single measurement for the identification process could be demonstrated. Minimum prediction error of blood oxygenation depends on blood gas level at the time of system identification i.e. the measurement situation. For severe pulmonary impairment, higher FiO _{2} settings were required to achieve a better prediction capability compared to less impaired pulmonary states. Plausibility analysis with real patient data could confirm this finding.
Discussion and conclusions
Dependent on patients’ pulmonary state, the influence of ventilator settings (here FiO _{2}) on model identification of the gas exchange model could be demonstrated. To maximize prediction accuracy i.e. to find the best individualized model with as few data as possible, best ranges of FiO _{2}settings for parameter identification were obtained. A less effort identification process, which depends on the pulmonary state, can be deduced from the results of this identifiability analysis.
Background
Mechanical ventilation is a lifesaving intervention in intensive care, maintaining pulmonary function in critically ill patients. Appropriate ventilator settings need to be found by the clinician to ensure both sufficient oxygenation and carbon dioxide removal. Target values for arterial partial pressures of oxygen (PaO _{2}) and carbon dioxide (PaCO _{2}) can be reached by changing inspired oxygen fraction (FiO _{2}) and minute ventilation (MV). Removal of CO_{2} and therefore PaCO _{2} in the patient is mainly regulated by adjusting MV. In critically ill patients, e.g. patients suffering from acute respiratory distress syndrome (ARDS), high levels of FiO _{2} and appropriate PEEP are usually necessary to ensure sufficient oxygenation. Finding the appropriate FiO _{2} setting follows a trialanderror approach that may not only be tedious but also exposes the patient to the potential risk of hypoxia and hyperoxia [1–4]. Pulse oximetry allows a continuous measurement of peripheral oxygen saturation (SpO _{2}), however this method has limitations in sensitivity and accuracy due to calibration assumptions, optical interference, and signal artifact [5]. Therefore, an invasive blood gas analysis is required at the end of each trial to evaluate the individual effect of a change in FiO _{2} accurately.
In mechanical ventilation therapy, both the risk of ventilator induced lung injury (VILI) and the effort to find adequate settings may be reduced if medical decision support would provide recommendations on how to adjust a patient’s settings to reach a prescribed treatment goal. In general, decision support can be divided into knowledge based (KDSS) and model based systems (MDSS). KDSS builds on rules of typical i.e. average patient behavior to represent reactions to changes of ventilator settings. In contrast, MDSS that are adapted to patient specific physiologic properties can simulate the individual reaction to changes in therapy settings. Using the inverse model, MDSS therefore suggests individualized ventilator settings by evaluating the approximated physiology of the patient.
Parameters of a model contain compact information about the individual patient state and dynamics once they are quantified in a parameter identification process (PIP). Parameter identification requires information from patient measurements often obtained during certain clinical maneuvers. Success and robustness of the PIP strongly depends on the properties of the model to reflect the required dynamics of the patient, the signal quality and amount of data available at the bedside. As the identified parameters are used for forward calculation of the model equations, they directly influence prediction performance of the model. While using multiple and even redundant measurements helps compensating noise induced errors, measurement efforts and applying the necessary maneuvers should not interfere with clinical processes. Thus, those measurements should be kept to the minimum necessary to ensure a robust PIP.
Models of pulmonary gas exchange are able to predict the effect of FiO _{2} and MV on PaO _{2} and PaCO _{2} in the patient. Oneparameter models [6, 7] usually only consider shunt, i.e. the amount of venous blood that is mixed with the oxygenated blood, to describe a patient’s oxygenation status and to predict the effect of an increase of FiO _{2} on PaO _{2}. However, using only one parameter to describe gas exchange impairments fails at low FiO _{2} when mismatches between alveolar ventilation (V̇) and perfusion (Q) occur. Several studies [8, 9] have come to the conclusion that PaO _{2} /FiO _{2} ratio, usually used to categorize lung impairment, changes with FiO _{2}. Thus, besides shunt, mathematical models of gas exchange should either include a parameter to describe oxygen diffusion limitation [10, 11] or a parameter to characterize \(\dot{V}/Q\) mismatch [9, 12–14]. Latter have shown to reproduce measurements at different oxygenation levels with sufficient accuracy compared to more complex models or MIGET measurements [15]. A twoparameter model of pulmonary gas exchange including shunt and \(\dot{V}/Q\) mismatch has previously been published by Kjaergaard et al. [12]. Karbing et al. [16] evaluated this model with data of severely ill intensive care patients. The model has been found to be identifiable with four pulse oxymetric (SpO _{2}) measurements at different levels of FiO _{2} and one blood gas analysis (BGA) providing PaO _{2} and PaCO _{2} together with the acid–base parameters pH, base excess (BE) and the hemoglobin concentration (cHb) as well as the endtidal gas fractions of oxygen (FetO _{2}) and carbon dioxide (FetCO _{2}). Although systems have been built to perform the necessary measurements in 10–15 min [17], lowering the number of measurements required for identification and therefore minimizing the required time and effort is highly relevant. Therefore we investigated if the number of measurements that are necessary to identify the model can be reduced to one FiO _{2}setting. Additionally, we evaluated the influence of the chosen level of FiO _{2} during model identification.
Methods
Gas exchange model with \(\dot{V}/Q\)mismatch and shunt
The mathematical model of human pulmonary gas exchange consists of two alveolar compartments that are perfused and ventilated and one shunt compartment that is perfused but not ventilated. The alveolar compartments are separated into a compartment with high \(\dot{V}/Q\)ratio and a compartment with low \(\dot{V}/Q\)ratio. This allows the consideration and simulation of limitations in gas exchange for both oxygen (O_{2}) and carbon dioxide (CO_{2}) concentrations in blood. Shunt, i.e. the fraction of venous blood not participating in gas exchange, is quantified by model parameter f _{ s } multiplied with blood flow Q. 90 % of the nonshunted blood ((1 − f _{ s })*Q) is distributed to the low \(\dot{V}/Q\) compartment, 10 % of the nonshunted blood is delivered to the high \(\dot{V}/Q\) compartment. Model parameter f _{ A } represents the fraction of alveolar ventilation V̇ _{ A } that reaches the low \(\dot{V}/Q\) compartment. Figure 1 shows the model structure of the pulmonary gas exchange model.
The model assumes equilibrium in blood gas concentrations as well as constant alveolar ventilation and perfusion without separating ventilation into phases of inspiration and expiration. Model inputs are FiO _{2} as well as endtidal blood gas fractions of oxygen (FetO _{2}) and carbon dioxide (FetCO _{2}), respectively. Inspired carbon dioxide is set to 0. Tidal volume V _{ tid } and respiratory frequency f _{ R } are assumed to be constant during simulation and are provided as additional model inputs. Model outputs are the resulting arterial blood gas parameters PaO _{2} and PaCO _{2}.
Alveolar ventilation V̇ _{ A } is calculated from f _{ R } and the difference between tidal volume V _{ tid } and the anatomic dead space volume V _{ ds }:
FetO _{2} and FetCO _{2} are composed of alveolar gas fractions FAO _{2} and FACO _{2} in both compartments, such that
Index x represents O_{2} and CO_{2} in Eq. 2 and in all following equations. Index 1 refers to the alveolar compartment with high \(\dot{V}/Q\), while index 2 denotes the low \(\dot{V}/Q\) compartment. Oxygen consumption V̇ _{ O2} and carbon dioxide production V̇ _{ CO2} are derived from alveolar air flow to each of the compartments and the difference between inspired and alveolar gas fractions as described in Eqs. (3) and (4):
Capillary blood gas concentrations Cc _{ x } are derived from alveolar gas fractions using O_{2} and CO_{2} dissociation curves [18, 19] (T—temperature):
Venous blood gas concentrations are then calculated as described in Eqs. (6) and (7).
Cardiac output Q is measured at the bedside or estimated from the patient’s body surface area.
FAO _{2} and FACO _{2} are solved numerically for a given measured FetO _{2} and FetCO _{2} with the condition that venous concentration in both compartments is equal. Finally, arterial blood gas concentrations CaO _{2} and CaCO _{2} are calculated as:
Arterial partial pressures of oxygen and carbon dioxide are then calculated from the reversed dissociation curves:
Model simulation
Forward calculation of the model equations is termed as model simulation. The flowchart of the model simulation process M _{ Θ } is depicted in Fig. 2 on the left. Vector Ψ summarizes physical constants for measurements needed for model simulation:
Here, minute ventilation MV is calculated as:
Vector Θ includes model parameters f _{ s } and f _{ A } as well as Ψ:
Output values PaO _{2,sim } and PaCO _{2,sim } are calculated depending on Θ and FiO _{2}:
Model identification
Parameters that need to be identified in the presented model are shunt fraction (f _{ s }) and the fraction of alveolar ventilation that is distributed to the alveolar compartment with low \(\dot{V}/Q\)ratio (f _{ A }). The process of model identification \((\tilde{M}_{a})\) is shown on the right of Fig. 2. \((\tilde{M}_{a})\) is a minimization process of an objective function. \(\widetilde{FiO}_{2}\), i.e. the level of inspired oxygen at the time of the measurement, as well as the other constant physiological values required during identification, are represented in vector α:
\(PaO_{2,meas}\) and \(PaCO_{2,meas}\) are the measured blood gases obtained at a specific condition described by α. They are used to determine f ^{*}_{ s } and f ^{*}_{ A } that best reproduce the measurements in the forward model:
Parameter identification was performed by minimizing the sum of the squared error (SSE) between measured (meas) and predicted (pred) partial blood gas pressures in arterial blood:
The weighting factor of 3 for PaCO _{2} was chosen to avoid imbalanced influence of PaO _{2} data on the identification process, as dimension of PaCO _{2} is approximately three times smaller than PaO _{2}. Minimization of the above described objective function was carried out using fminsearchbnd in MATLAB (R2012a, The Mathworks, Natick, MA, USA). fminsearchbnd is distributed under the BSD license and is based on fminsearch, the MATLAB function that employs the NelderMead simplex search method [20]. According to [21] a shunt of 50 % and above leads to increases in FiO _{2} having no effect on PaO _{2}. Additionally, f _{ A } values above 0.9 lead to a swap of the high \(\dot{V}/Q\)compartment with the low \(\dot{V}/Q\)compartment, essentially mirroring \(\dot{V}/Q\)values of f _{ A } below 0.9. Thus parameter constraints for {f _{ s }, f _{ A }} were set as nonnegative lower boundaries LB = {0, 0} and as upper boundaries UB = {0.5, 0.9}.
Structural identifiability of the model using multiple measuring points was shown in a previous report [22]. Initial fs was set to 0.2, which showed to lead to the global minimum of the objective function in all test cases. Initial f _{ A } was set arbitrarily to 0.5 to start the minimization process at a certain initial mismatch between ventilation and perfusion. Constant patient state for the time of model prediction is assumed.
Model prediction
Prediction of blood gas levels depending on FiO _{2} is done using forward calculation of the model M_{ Θ* } with Θ ^{*} = (f ^{*}_{ s } , f ^{*}_{ A } , Ψ):
Data
We have employed both simulated and recorded real patient data to evaluate how well the described model is identifiable with data obtained at one single FiO _{2} level.
Simulated data The same twoparameter model of pulmonary gas exchange was used to create experimental data. It allows calculating the impact of noise in the data because the correct results for parameter identification are known a priori. Twelve classes of patient data sets have been generated, that differ in the parameter combinations of f _{ s } and f _{ A }. Those have been chosen to represent different stages of pulmonary disease. Model parameters used for data generation, resulting FiO _{2} /PaO _{2}ratios and the classifications of pulmonary impairment [24] are listed in Table 1. BGA and physiological standard values of an adult man were used for data generation. These physical constants are listed and explained in Table 2.
More formally, we define twelve patient classes by the model parameters
For each of the twelve patient classes, 1000 simulated measurements equidistant between FiO _{2} of 21 % and 100 % were determined. Depending on FiO _{2} settings, model simulation led to
Measuring PaO _{2} and PaCO _{2} via blood gas samples drawn from the arterial line is the current gold standard in clinical practice [25, 26], while measuring arterial oxygen saturation via pulse oximetry is accurate within ±2 % of the true value [27]. Thus, to account for measurement noise that would be present in a real setting, both PaO _{2,sim } and PaCO _{2,sim } data were superimposed with uniformly distributed noise in a range of ±5 %:
The quality of system identification was assessed with a test set (ts) of 17 distinct \(FiO_{2}^{ts}\) values ranging from 21 to 100 % in steps of 5 %:
Real patient data: Two real patient data sets were used for the plausibility check of the results obtained from the theoretical analysis. Real patient data including at least four blood gas measurements at different FiO _{2} settings in mandatory ventilation mode were retrieved from a patient data management system of the university medical centre in Kiel [7]. Two data sets, a mild (Pat R1) and a critically ill patient (Pat R2), met those demands. The recorded levels of FiO _{2} were applied on a therapeutical basis, thus not systematically in the context of a clinical trial. The data sets included measurements of PaO _{2}, PaCO _{2}, f _{ R }, Q, V _{ tid }, V _{ ds }, cHb, pH, T and FetCO _{2} at each of the applied FiO _{2} levels. Patient data did not include FetO _{2} measurements, thus FetO _{2} was approximated from:
Here, the respiratory quotient RQ was assumed to be 0.8. As with the simulated data sets, initial conditions of {fs, fA} for parameter identification were set to {0.2, 0.5}.
Analysis of structural identifiability
To verify structural identifiability of a model system, uniqueness of the solution of parameter identification has to be proven. The simplicity of the twoparameter model allows a numerical calculation and twodimensional visualization of the objective function. The SSE is calculated and plotted for different parameter combinations to visualize the contour of the error surface. A single global minimum of the objective function indicates structural identifiability of the model.
Structural identifiability of the model using one measurement point at one \(\widetilde{FiO}_{2}\) level was analyzed with the synthetic as well as the two real patient data sets. The error surfaces (SSE) were plotted as a function of model parameters f _{ s } and f _{ A } with a resolution of 90 × 90.
Evaluation of quality of fit
Besides quantity, quality of measurements used for model identification is essential for the accuracy of parameter identification. To verify practical identifiability, the influence of measuring errors on identification behavior of the model system was evaluated using the 1000 virtual measurements (\(PaO_{2,meas}^{j} , PaCO_{2,meas}^{j}\)) in each of the eight classes of virtual patients.
For every patient class the gas exchange model was identified with only one single of the 1000 noisy measurements. Please note that each of those measurements belongs to one FiO _{2} setting:
With the identified parameters f ^{*}_{ s } and f ^{*}_{ A } , the evaluations \(PaO_{{2,pred_{i} }}^{ts}\) as well as \(PaCO_{{2,pred_{i} }}^{ts}\) were calculated for all 1 ≤ i ≤ 17 FiO _{2}values in the test set:
Predictive performance for blood oxygenation as well as partial pressure of carbon dioxide was evaluated by comparing \(PaO_{2,sim}^{ts}\) and \(PaCO_{2,sim}^{ts}\), i.e. the values of the original simulation set, with \(PaO_{2,pred}^{ts}\) and \(PaCO_{2,pred}^{ts}\) respectively. Mean deviations were calculated with
and
For statistical evaluation, the 1000 \(\widetilde{FiO}_{2}\) levels, ranging from 21 % to 100 %, were divided into eight clusters (21–30, 30–40,…, 90–100 %), each cluster containing 125 values. Mean and standard deviations of \(\overline{\Delta PaO}_{2}\) for each cluster were calculated.
Verification of results with patient data
Two real patient data sets were used to confirm the findings obtained with simulated data. Identification was conducted at each FiO _{2} value that was recorded in the particular patient. Predictive performance was evaluated by comparing measured and predicted PaO _{2} and PaCO _{2} at all four recorded FiO _{2} levels:
Results
Visualizing the objective function
Figure 3 visualizes the contour of the objective function evaluated at one single measuring point. Figure 3a shows synthetic patient data, while Fig. 3b is devoted to real patient data. The contour lines (SSE) are scaled logarithmically to improve visibility of the minimum. For all analyzed data sets, a single global minimum in the error surface as a function of the model parameters f _{ s } and f _{ A } could be detected. The parameter combination leading to the global minimum was in agreement with the parameters used for data generation. Within accuracy of numerical representation (double), the SSE value was zero at the minimum. The global minimum is inside a narrow and flat valley, parallel to the axis of parameter f _{ A }. A similar type of shape of the objective function could be shown for the analysis of both simulated and real data.
Prediction of PaCO _{2} and PaO _{2} depending on FiO _{2}
Over all tested data sets, mean \(\overline{\Delta PaCO}_{2}\) was 2.4 % (1.3 mmHg) of the true value with a standard deviation of ±1.6 % (±0.8 mmHg). Figure 4a shows mean deviation of \(\overline{\Delta PaO}_{2}\) with respect to α. PaO _{2} prediction was less accurate when identification data were recorded at low FiO _{2} levels, especially for data sets representing pulmonary impairment. For data set 1 (healthy lung) minimal \(\overline{\Delta PaO}_{2}\) was achieved for 40 % < FiO _{2} < 50 % whereas model identification with data set 12 (severely impaired lung) shows the smallest \(\overline{\Delta PaO}_{2}\) for α = (100 %, Ψ). In the minima, the model was able to reproduce PaO _{2} of all of the simulated patient data sets with a mean deviation of less than 2.5 % (<2.5 mmHg) of the true value with a standard deviation of less than 3 % (<2 mmHg). In Fig. 4a, the minima are marked with vertical lines and the respective patient numbers.
Mean deviation of ΔPaO _{2} for real patient data sets is shown in Fig. 4b. Predictions with model parameters being identified at low and high FiO _{2} show higher deviations from measured values than for identification at medium FiO _{2} levels. Mean deviation of ΔPaO _{2} was less than 10 % (8 % or 5.8 mmHg for Patient R1 and 6 % or 4.6 mmHg for Patient R2) at the minimum. The best performance of PaO _{2} prediction was found for α = (40 %, Ψ) (Pat R1) and α = (80 %, Ψ) (Pat R2) respectively.
Figure 5 summarizes the FiO _{2} levels leading to a minimum of mean deviation of \(\overline{\Delta PaO}_{2}\) for simulated data sets (colored markers). 32 additional simulated data sets (black markers) were generated to illustrate the relation between optimal FiO _{2} and pulmonary impairment more precisely. Location of the minimum was shifted to a higher FiO _{2} cluster with increasing pulmonary impairment, i.e. decreasing PaO _{2}/FiO _{2}ratio.
Discussion
Using mathematical models for decision support in clinical practice requires high level of safety and accuracy of the model predictions. Furthermore, measuring effort for model identification, i.e. identification of patient specific parameters, should be kept to a minimum.
Identification of the twoparameter model of human pulmonary gas exchange requires data from arterial blood gas analysis and photoplethysmographic saturation measurement. The model has previously been applied with four measurements at different levels of FiO _{2}. To reduce the effort required in clinical practice to identify gas exchange models it was investigated if a reduced number of FiO _{2} levels for data acquisition is sufficient.
Structural identifiability of the model applying a single identification data point was examined in this study with both simulated and real patient data sets. One single global minimum of the objective function is an indication for structural identifiability of the model. Without noise in the data, one single blood gas measurement is sufficient for robust model identification.
The effect of a change in FiO _{2} on the concentration of carbon dioxide in arterial blood is negligible. PaCO _{2} data could be reproduced by the model with high accuracy. Prediction error of CO _{2} data was below the noise level of ±5 % for all data sets. However, measuring errors may decrease the accuracy of parameter identification and therefore model prediction of PaO _{2}, especially if identification is based on only one measurement. Results show that the gas exchange model with shunt and \(\dot{V}/Q\)mismatch is able to fit both the synthetic and the real patient data with good accuracy, as already presented in former work [16]. Oxygenation data of all data sets could be reproduced by the model with a mean deviation below 10 % in spite of measuring errors in the identification data. However, when identifying the model with noisy data, the FiO _{2} setting influences the prediction accuracy of PaO _{2}. This influence was therefore examined to find a guideline how to choose an appropriate FiO _{2} level for data acquisition in the identification process.
The identification processes of both simulated and real patient data sets representing a variety of different disease states showed surprisingly similar results. It could be pointed out that accuracy of model prediction of blood gas concentration is related to the FiO _{2} setting when recording identification data. The optimum FiO _{2} level depends on the level of pulmonary impairment whereupon FiO _{2} should be increased with increasing severity of pulmonary impairment. In severely ill patients, oxygenation of arterial blood is inhibited, thus higher FiO _{2} levels are required in order to achieve adequate PaO _{2} levels.
Figure 6a shows the mean deviation of \(\overline{\Delta PaO}_{2}\) with respect to \(PaO_{2,meas}^{j}\). It could be observed that a minimal prediction error is achieved for \(PaO_{2,meas}^{j}\) in the range of 150–200 mmHg for the entire simulated data sets. Identification at higher levels leads to a small increase of mean deviation. However, mean deviation of \(\overline{\Delta PaO}_{2}\) was found to be still below 5 % for all tested patient cases. Identifying the model with PaO _{2} levels of less than 100 mmHg is potentially leading to an increase in both mean and standard deviation of model prediction. Because of the high severity of pulmonary impairment, data sets 7 and 8 do not achieve a PaO _{2} of 100 mmHg, even when a FiO _{2} of 100 % is applied. Patients with such high pulmonary impairment have to be ventilated using the highest FiO _{2} possible to achieve a sufficient oxygenation.
Mean deviation of ΔPaO _{2} over \(PaO_{2,meas}^{real}\) is shown in Fig. 6b. Both curves representing real data confirm the results of the analysis with simulated data, but minimum mean deviations were found for a PaO _{2,meas } of 73 mmHg and 81 mmHg respectively.
When only one measuring point is used for model identification, success of parameter identification obviously depends on the level of PaO _{2} at the particular time of the measurement. Low PaO _{2} levels bear the risk of large influence of measuring errors in identification data. In low PaO _{2} levels, even small changes in PaO _{2} used for identification may lead to an overestimation of shunt fraction f _{ s } and therefore to a smaller increase of PaO _{2} for higher FiO _{2} levels compared to correct data. In severely ill patients, this effect is more prominent than in patients with less pulmonary distress.
PaO _{2} levels suggested to be optimal for identification of the model may be not achievable in patients with severe pulmonary impairment. Here, identification at an FiO _{2} of 100 % was shown to achieve predictions with highest accuracy.
In the presented work, we have investigated a FiO _{2} range 21–100 %. To separate effects of shunt from those caused by low \(\dot{V}/Q\)ratio, subatmospheric oxygen levels should be considered as well. However, the intended use of the applied model and the presented routine of identifying the model parameters with only one blood gas measurement are in a critical care environment where such oxygen levels are not applied.
Karbing et al. [16] have previously presented a threeparameter extension of the model used in this work which uses an adjustable distribution of nonshunted blood among the alveolar compartments. This model shows to be superior in terms of reproducing PaCO _{2} especially in patients with \(\dot{V}/Q\)ratios below 0.27. The presented work however focuses mainly on finding the optimal calibration point of FiO _{2}, which has only a minor effect on PaCO _{2}. Nevertheless, investigating the structural identifiability and the possibility of using only one measurement set to also identify the threeparameter model should be considered in future work.
The model of gas exchange applied in the presented study includes assumptions such as steady state conditions of blood gases and constant alveolar ventilation and perfusion. Thus, tidal breathing as is the reality in humans is not considered. Those assumptions present shortcomings compared to the reality those models try to reproduce [28]. Several models including tidal breathing have been presented in the past [28–30] however continuous measurement of blood gases in combination with continuous measurements of gases in inspired and expired air is not routinely available at the bedside at this moment. Thus clinical applications are currently limited to models assuming continuous ventilation and perfusion as well as equilibrated blood gases. Still, Karbing et al. [16] have shown previously that the applied model is well capable of reproducing patient data for a wide range of lung impairments and that the model can thus be used in a clinical environment as a prediction tool.
Model based decision support in clinical practice implies that the mathematical model is identifiable with a minimum of measuring effort. We could show that the twoparameter gas exchange model with shunt and \(\dot{V}/Q\)mismatch is structural identifiable with only one single blood gas measurement. Using only one single measurement, possible measuring errors are not averaged. However, simulation results show that model based prediction of blood gases for different FiO _{2} is possible with a mean prediction error below 10 % for a maximum measuring error of 5 %. Simultaneously, we could determine the range of PaO _{2} level where prediction error is minimized for data representing a wide range of different pulmonary states. Our work provides scientific findings in developing a robust parameter identification process for the gas exchange model with low measuring effort. For a given accuracy of the blood gas measurements used for identification it will be possible to estimate the accuracy of the model prediction of blood gases.
This study is faced with the limitation that only two real patient data sets were available to confirm the findings from the study with simulated data. Furthermore, the real data were not from a systematic patient study, but retrieved from a patient data management system, giving no information about the interventions of the clinicians between the measurements. Therefore, a change in the patients’ pulmonary state cannot be excluded. A systematic study with a higher number of mechanically ventilated patients is necessary to consolidate our findings.
Conclusions
The study showed that the identification point has a significant impact on the predictive performance of the presented gas exchange model. Measuring errors, i.e. noise in identification data, could lead to prediction errors when only one measurement is applied. A combination of simulated and real patient data provides a valuable tool in determining the optimal identification point where influence of measurement errors is minimal.
Abbreviations
 ARDS:

acute respiratory distress syndrome
 BE :

base excess
 BGA:

blood gas analysis
 Cc :

capillary concentration
 cHB :

concentration of hemoglobin
 CO _{2} :

carbon dioxide
 Cv :

venous concentration
 f _{ A } :

fraction of ventilation distribution
 FA :

alveolar fraction
 Fet :

endtidal fraction
 FiO _{2} :

inspired oxygen fraction
 f _{ Q } :

fraction of perfusion distribution
 f _{ R } :

respiratory frequency
 f _{ s } :

shunt fraction
 KDSS:

knowledge based decision support system
 LB:

lower boundaries
 MBSS:

model based decision support system
 MIGET:

multiple inert gas elimination technique
 MV:

minute ventilation
 O _{2} :

oxygen
 Pa :

arterial partial pressure
 PIP:

parameter identification process
 Q :

cardiac output
 RQ :

respiratory quotient
 SpO _{2} :

peripheral oxygen saturation
 SSE :

sum of the squared errors
 T :

patient temperature
 UB:

upper boundaries
 V̇ :

ventilation
 V̇ _{ A } :

alveolar ventilation
 V̇ _{ CO2} :

carbon dioxide production
 V _{ ds } :

dead space volume
 VILI:

ventilator induced lung injury
 V̇ _{ O2} :

oxygen consumption
 V _{ tid } :

tidal volume
References
 1.
Edmark L, KostovaAherdan K, Enlund M, Hedenstierna G. Optimal oxygen concentration during induction of general anesthesia. Anesthesiology. 2003;98(1):28–33.
 2.
Nash G, Blennerhassett JB, Pontoppidan H. Pulmonary lesions associated with oxygen therapy and artificial ventilation. N Engl J Med. 1967;276(7):368–74.
 3.
Rinaldo JE, Rogers RM. Adult respiratorydistress syndrome: changing concepts of lung injury and repair. N Engl J Med. 1982;306(15):900–9.
 4.
Altemeier WA, Sinclair SE. Hyperoxia in the intensive care unit: why more is not always better. Curr Opin Crit Care. 2007;13(1):73–8.
 5.
Sinex JE. Pulse oximetry: principles and limitations. The American Journal of Emergency Medicine. 1999;17(1):59–66.
 6.
Benatar SR, Hewlett AM, Nunn JF. The use of isoshunt lines for control of oxygen therapy. Br J Anaesth. 1973;45(7):711–8.
 7.
Kretschmer J, Becher T, Riedlinger A, Schadler D, Weiler N, Moller K. A simple gas exchange model predicting arterial oxygen content for various FiO_{2} levels. Conf Proc IEEE Eng Med Biol Soc. 2013;2013:465–8.
 8.
Karbing DS, Kjaergaard S, Smith BW, Espersen K, Allerod C, Andreassen S, Rees SE. Variation in the PaO_{2}/FiO_{2} ratio with FiO_{2}: mathematical and experimental description, and clinical relevance. Crit Care. 2007;11(6):R118.
 9.
Whiteley JP, Gavaghan DJ, Hahn CE. Variation of venous admixture, SF6 shunt, PaO_{2}, and the PaO_{2}/FiO_{2} ratio with FiO_{2}. Br J Anaesth. 2002;88(6):771–8.
 10.
Andreassen S, Egeberg J, Schroter MP, Andersen PT. Estimation of pulmonary diffusion resistance and shunt in an oxygen status model. Comput Methods Programs Biomed. 1996;51(1–2):95–105.
 11.
Andreassen S, Rees SE, Kjaergaard S, Thorgaard P, Winter SM, Morgan CJ, Alstrup P, Toft E. Hypoxemia after coronary bypass surgery modeled by resistance to oxygen diffusion. Crit Care Med. 1999;27(11):2445–53.
 12.
Kjaergaard S, Rees SE, Nielsen JA, Freundlich M, Thorgaard P, Andreassen S. Modelling of hypoxaemia after gynaecological laparotomy. Acta Anaesthesiol Scand. 2001;45(3):349–56.
 13.
Petros AJ, Dore CJ, Nunn JF. Modification of the isoshunt lines for low inspired oxygen concentrations. Br J Anaesth. 1994;72(5):515–22.
 14.
Riley RL, Cournand A. Ideal alveolar air and the analysis of ventilationperfusion relationships in the lungs. J Appl Physiol. 1949;1(12):825–47.
 15.
Loeppky JA, Caprihan A, Altobelli SA, Icenogle MV, Scotto P, Vidal Melo MF. Validation of a twocompartment model of ventilation/perfusion distribution. Respir Physiol Neurobiol. 2006;151(1):74–92.
 16.
Karbing DS, Kjaergaard S, Andreassen S, Espersen K, Rees SE. Minimal model quantification of pulmonary gas exchange in intensive care patients. Med Eng Phys. 2011;33(2):240–8.
 17.
Thomsen LP, Karbing DS, Smith BW, Murley D, Weinreich UM, Kjaergaard S, Toft E, Thorgaard P, Andreassen S, Rees SE. Clinical refinement of the automatic lung parameter estimator (ALPE). J Clin Monit Comput. 2013;27(3):341–50.
 18.
Kelman GR. Digital computer subroutine for the conversion of oxygen tension into saturation. J Appl Physiol. 1966;21(4):1375–6.
 19.
Sharan M, Singh MP, Aminataei A. A mathematical model for the computation of the oxygen dissociation curve in human blood. Biosystems. 1989;22(3):249–60.
 20.
Lagarias JC, Reeds JA, Wright MH, Wright PE. Convergence properties of the Nelder–Mead simplex method in low dimensions. SIAM J Optim. 1998;9(1):112–47.
 21.
Lumb AB. Distribution of pulmonary ventilation and perfusion. In: Nunn’s applied respiratory physiology. 7th ed. Edinburgh: Elsevier; 2010.
 22.
Riedlinger A, Schranz C, Möller K: Robustness analysis of a mathematical gas exchange model. Biomed Tech (Berl) 2013; 58 (Suppl 1).
 23.
Schranz C, Knobel C, Kretschmer J, Zhao Z, Moller K. Hierarchical parameter identification in models of respiratory mechanics. IEEE Trans Biomed Eng. 2011;58(11):3234–41.
 24.
Ranieri VM, Rubenfeld GD, Thompson BT, Ferguson ND, Caldwell E, Fan E, Camporota L, Slutsky AS. Acute respiratory distress syndrome: the Berlin definition. JAMA. 2012;307(23):2526–33.
 25.
Byrne AL, Bennett M, Chatterji R, Symons R, Pace NL, Thomas PS. Peripheral venous and arterial blood gas analysis in adults: are they comparable? A systematic review and metaanalysis. Respirology. 2014;19(2):168–75.
 26.
Honarmand A, Safavi M. Prediction of arterial blood gas values from arterialized earlobe blood gas values in patients treated with mechanical ventilation. Indian J Crit Care Med. 2008;12(3):96–101.
 27.
Criticare Systems inc: Gebrauchsanweisung 504/504US, 504P/504USP. 1993.
 28.
Hahn CE, Farmery AD. Gas exchange modelling: no more gills, please. Br J Anaesth. 2003;91(1):2–15.
 29.
Benallal H, Denis C, Prieur F, Busso T. Modeling of endtidal and arterial PCO_{2} gradient: comparison with experimental data. Med Sci Sports Exerc. 2002;34(4):622–9.
 30.
Whiteley JP, Farmery AD, Gavaghan DJ, Hahn CE. A tidal ventilation model for oxygenation in respiratory failure. Respir Physiol Neurobiol. 2003;136(1):77–88.
Authors’ contributions
AR carried out the analysis and drafted the manuscript. JK participated in the implementation of the model as well as representation and appraisal of the results. KM conceived of the study, and participated in the design of the study. All authors read and approved the final manuscript.
Acknowledgements
This work was partially supported by the Bundesministerium für Bildung und Forschung under Grant 01IB10002D (WiMVent) as well as the FP7IRSES 318943 (eTime). The authors express their gratitude to Dr. Schädler and Dr. Becher from the Department of Anesthesiology and Intensive Care Medicine at the University Medical Center of Schleswig–Holstein, Campus Kiel, Germany, for providing the anonymous made patient data for this study.
Compliance with ethical guidelines
Competing interests There are no significant competing interests that might have influenced the results of the work described in this manuscript.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Riedlinger, A., Kretschmer, J. & Möller, K. On the practical identifiability of a twoparameter model of pulmonary gas exchange. BioMed Eng OnLine 14, 82 (2015). https://doi.org/10.1186/s1293801500776
Received:
Accepted:
Published:
Keywords
 Physiological modeling
 Gas exchange model
 Parameter identification