METHOD FOR DETERMINING THE RESISTANCE OF THE RESPIRATORY SYSTEM OF A PATIENT

• US 20080234595A1
• Filed: 03/14/2008
• Published: 09/25/2008
• Est. Priority Date: 03/20/2007
• Status: Active Grant

• Current Assignee: Dragerwerk AG & Company KGaA (FPS Vermögensverwaltung GmbH)
• Original Assignee: Drager Medical AG & Co KG (Dragerwerk AG)
• Inventors: RANIERI, Marco, MICELLI, Cosimo, APPENDINI, Lorenzo
• Granted Patent: US 8,160,817 B2
• Time in Patent Office: Days
• US Class Current: 600/538
• CPC Class Codes: A61B 5/085 Measuring impedance of respiratory organs or lung elasticity
• CPC Class Codes: A61M 16/0051 with alarm devices
• CPC Class Codes: A61M 16/024 including calculation means, e.g. using a processor
• CPC Class Codes: A61M 2016/0027 pressure meter
• CPC Class Codes: A61M 2016/0036 in the breathing tube and used in both inspiratory and expiratory phase
• CPC Class Codes: A61M 2230/46 Resistance or compliance of the lungs

Abstract

A method for automatically determining the resistance of the respiratory system of a spontaneously breathing intubated patient. The method includes measuring flow (V′) and airway pressure (P_aw) as a function of time; determining the airway pressure at the end of the occlusion period as end of occlusion airway pressure P_aw(t_c)=P_occl; analyzing a time dependence of the flow starting at the end of the occlusion period to determine a transition time t_d at which the slope of the flow changes from a first high value to a second lower value, and/or analyzing a time dependence of the airway pressure starting at the end of the occlusion period to determine a transition time t_d at which the slope at the airway pressure changes from a first high value to a second lower value, or by detecting an increase of the airway pressure over a predetermined threshold value to determine the transition time t_d; determining the flow at the transition time as transition time flow V′(t_d)=V′_peak and determining the airway pressure at the transition time as transition time pressure P_aw(t_d); and calculating the resistance on the basis of the ratio of an estimated driving pressure at the transition time P_drv(t_d) and the transition time flow V′_peak, wherein the value of the transition time driving pressure P_drv(t_d) is calculated as sum of the transition time airway pressure P_aw(t_d) and the muscular effort P_mus(t_d) at the transition time, which muscular effort P_mus(t_d) is extrapolated based on a predetermined time dependence of P_mus(t) during the occlusion period and based on value of the end of occlusion airway pressure P_occl.

Claims

1. A method for automatically determining the resistance of the respiratory system of a spontaneously breathing intubated patient, wherein the method is conducted with or without assisted ventilation and includes an occlusion period at the start of inspiration, the method comprising:

measuring flow (V′

) and airway pressure (P_aw) as a function of time;

determining the airway pressure at the end of the occlusion period as end of occlusion airway pressure (P_aw(t_c)=P_occl);

analyzing a time dependence of the flow starting at the end of the occlusion period to determine a transition time (t_d) at which the slope of the flow changes from a first high value to a second lower value, and/or analyzing a time dependence of the airway pressure starting at the end of the occlusion period to determine the transition time (t_d) at which the slope at the airway pressure changes from a first high value to a second lower value, or by detecting an increase of the airway pressure over a predetermined threshold value to determine the transition time (t_d);

determining the flow at the transition time as transition time flow (V′

(t_d)=V′

_peak) and determining the airway pressure at the transition time as transition time pressure (P_aw(t_d));

calculating the resistance on the basis of the ratio of an estimated driving pressure at the transition time (P_drv(t_d)) and the transition time flow (V′

_peak), wherein the value of the transition time driving pressure (P_drv(t_d)) is calculated as sum of the transition time airway pressure (P_aw(t_d)) and a muscular effort pressure (P_mus(t_d)) at the transition time, which muscular effort pressure (P_mus(t_d)) is extrapolated based on a predetermined time dependence of muscular effort pressure (P_mus(t)) during the occlusion period and based on value of the end of occlusion airway pressure (P_occl).

2. The method according to claim 1, wherein the predetermined time dependence is selected such that muscular effort pressure (P_mus(t)) is constant after the occlusion (P_mus(t_d)=−

P_occl), and thus the extrapolated driving pressure is determined to be the transition time airway pressure less the occlusion airway pressure (P_drv(t_d)=P_aw(t_d)−

P_occl).

3. A method according to claim 1, wherein the predetermined time dependence is selected such that muscular effort pressure (P_mus(t)) rises after the occlusion with constant slope (k) after the occlusion (P_mus(t_d)=−

P_occl+k*(t_d−

t_c)), and thus the extrapolated driving pressure is determined to be transition time airway pressure less the occlusion airway pressure plus the constant (P_drv(t_d)=P_aw(t_d)−

P_occl+k*(t_d−

t_c)).

4. The method according to claim 1, wherein the predetermined time dependence is selected such that muscular effort pressure P_mus(t) proceeds after the occlusion with a curvature that is described by an n^th order polynom (P_mus(t_d)=−

P_occl+k₁*(t_d−

t_c)+k₂*(t_d−

t_c)²+k₃*(t_d−

t_c)³+ . . . +k_n*(t_d−

t_c)ⁿ and thus P_drv(t_d)=P_aw(t_d)−

P_occl+k₁*(t_d−

t_c)+k₂*(t_d−

t_c)²+k₃*(t_d−

t_c)³+ . . . +k_n*(t_d−

t_c)ⁿ).

5. The method according to claim 3, wherein the P_mus(t_d) is corrected according to a-priori known pressure-flow characteristics of the respiratory muscles by applying a time dependent gain factor on the P_mus signal.

6. The method according to claim 4, wherein P_mus(t) is corrected according to the a-priori known pressure-flow characteristics of the respiratory muscles by applying a time dependent gain factor on the P_mus signal.

7. The method according to claim 1, wherein the step of analyzing the time dependence of the flow data or the airway pressure data includes fitting of a predetermined non-linear function to the respective time-dependent data, said predetermined non-linear function having a first essentially linear portion and a second essentially linear portion with a smooth transition in between, wherein the slopes and lengths and the location of the smooth transition are determined by parameters which are varied in the fitting procedure.

8. The method according to claim 1, wherein the step of analyzing the time dependence of the flow or of the airway pressure data includes fitting two straight lines to respective the time-dependent data and determining the intersection of the two straight lines as the transition time value t_a.

9. The method according to claim 1, wherein said step of analyzing the time dependence of the flow data or the airway pressure data includes a linear filtering procedure of the respective time-dependent data in order to determine the change of slope point as the transition time value t_d.

10. A method for automatically determining the elastance (E) of the respiratory system of a spontaneously breathing intubated patient, wherein said method is conducted with or without assisted ventilation and includes an occlusion period at the start of inspiration, wherein said method comprises:

a) determining the resistance (K) of the respiratory system with or without assisted ventilation and including an occlusion period at the start of inspiration by;

measuring flow (V′

) and airway pressure (P_aw) as a function of time;

determining the airway pressure at the end of the occlusion period as end of occlusion airway pressure (P_aw(t)=P_occl);

analyzing a time dependence of the flow starting at the end of the occlusion period to determine a transition time (t_d) at which the slope of the flow changes from a first high value to a second lower value, and/or analyzing a time dependence of the airway pressure starting at the end of the occlusion period to determine the transition time (t_d) at which the slope at the airway pressure changes from a first high value to a second lower value, or by detecting an increase of the airway pressure over a predetermined threshold value to determine the transition time (t_d);

determining the flow at the transition time as transition time flow (V′

(t_d)=V′

_peak) and determining the airway pressure at the transition time as transition time pressure (P_aw(t_d));

calculating the resistance on the basis of the ratio of an estimated driving pressure at the transition time (P_drv(t_d)) and the transition time flow (V′

_peak), wherein the value of the transition time driving pressure (P_drv(t_d)) is calculated as sum of the transition time airway pressure (P_aw(t_d)) and a muscular effort pressure (P_mus(t_d)) at the transition time, which muscular effort pressure (P_mus(t_d)) is extrapolated based on a predetermined time dependence of muscular effort pressure (P_mus(t)) during the occlusion period and based on value of the end of occlusion airway pressure (P_occl);

b) measuring the inspiratory volume as a function of time;

c) determining the airway pressure (P_aw(t)) during the occlusion period (t_b<

t<

t_c) as a function of time (P_aw(t));

d) extrapolating the airway pressure during the occlusion (P_aw(t)) with occlusion period (t_b<

t<

t_c) data to the time when the inspiratory volume reaches a predetermined volume (t_c) to obtain an extrapolated airway pressure (P_extr(t_e));

e) subtracting the extrapolated airway pressure (P_extr(t_e)) from the ventilator pressure (P_aw(t_e)) to obtain an estimated driving pressure (P_drv(t_e));

f) calculating the resistive pressure (P_res(t_e)) at the predetermined inspiratory volume (⅔

VT) using the measured flow at that inspiratory volume and the resistance value determined in step a), and subtracting the resistive pressure (P_res(t_e)) from the estimated driving pressure (P⁻

_extr(t_e)) to determine the elastic pressure (P_el(t_e)), and determining the elastance on the basis of the ratio of elastic pressure (P_el(t_e)) to the predetermined inspiratory volume (V(t_e)).

11. The method according to claim 10, wherein said extrapolation of the time dependent airway pressure data is performed by fitting a predetermined function to the airway pressure data of the occlusion period and obtaining the extrapolated pressure by calculating the fitted function at the extrapolation time t_e where the inspiratory volume reaches a predetermined value.

12. The method according to claim 11 wherein a predetermined linear convex or concave function is used in the fitting procedure for extrapolation.

13. The method according to claim 10, wherein:

steps d) to f) are repeated for a sequence of predetermined volumes (x_i−

VT, x₁ . . . , x_N≦

1);

a corresponding sequence of elastance values (E(x_i−

VT)) is determined; and

an elastance value is selected from a region of the sequence where the dependence of the elastance from the predetermined volume is lowest.

14. The method according to claim 11, wherein:

steps d) to f) are repeated for a sequence of predetermined volumes (x_i−

VT, x₁ . . . , x_N≦

1);

a corresponding sequence of elastance values (E(x_i−

VT)) is determined; and

an elastance value is selected from a region of the sequence where the dependence of the elastance from the predetermined volume is lowest.

15. The method according to claim 12, wherein:

steps d) to f) are repeated for a sequence of predetermined volumes (x_i−

VT, x₁ . . . , x_N≦

1);

a corresponding sequence of elastance values (E(x_i−

VT)) is determined; and

an elastance value is selected from a region of the sequence where the dependence of the elastance from the predetermined volume is lowest.

16. The method according to claim 10, wherein the values for resistance (R) and elastance (E) as determined in a previous breathing cycle are used to calculate the muscular pressure ((P_mus)), from the inspiratory volume (V) and flow (V′

) (P_mus+P_aw=E·

V+R·

V′

, P_aw=VA·

V+FA·

V′

) as a function of inspiratory volume to compare it to the extrapolated pressure value to check the extrapolation.

17. The method according to claim 11, wherein the values for resistance (R) and elastance (E) as determined in a previous breathing cycle are used to calculate the muscular pressure ((P_mus)), from the inspiratory volume (V) and flow (V′

)(P_mus+P_aw=E·

V+R·

V′

, P_aw=VA·

V+FA·

V′

) as a function of inspiratory volume to compare it to the extrapolated pressure value to check the extrapolation.

18. The method according to claim 15, wherein the values for resistance (R) and elastance (E) as determined in a previous breathing cycle are used to calculate the muscular pressure ((P_mus)), from the inspiratory volume (V) and flow (V′

)(P_mus+P_aw=E·

V+R·

V′

, P_aw=VA·

V+FA·

V′

) as a function of inspiratory volume to compare it to the extrapolated pressure value to check the extrapolation.

19. A method for determining a resistance of a respiratory system of a patient, the method comprising the steps of:

inhaling through a breathing apparatus;

blocking a flow through the breathing apparatus at a start of said inhaling, said blocking being performed during an occlusion period;

measuring flow through the breathing apparatus during said inhaling as a function of time;

measuring airway pressure in the breathing apparatus during said inhaling as a function of time;

determining the airway pressure at an end of the occlusion period (P_occl);

determining a transition time (t_d) during said inhaling and after the occlusion period, the transition time being determined by one of,

determining a time when a slope of the flow changes from a first high value to a second lower value,

determining a time when a slope of the airway pressure changes from a first high value to a second lower value, and

determining a time when the airway pressure increases above a predetermined threshold;

determining the flow at the transition time t_d as V′

_peak;

determining the airway pressure at the transition time t_d as P_aw(t_d);

determining muscular pressure at the transition time (P_mus(t_d)) by extrapolation of the airway pressure during the occlusion period;

calculating a driving pressure at the transition time (P_drv(t_d)) by adding P_aw(t_d) and P_mus(t_d);

calculating the resistance of the respiratory system based on a ratio of P_drv(t_d) to V′

_peak.