# Patent

## Title
Distinguishing closed and open respiratory airway apneas by complex admittance values

## Abstract

Methods and apparatus are disclosed for determining the occurrence of a closed or open apnea. Respiratory air flow from a patient is measured to give an air flow signal. The determination of an apnea is performed by applying an oscillatory pressure waveform of known frequency to a patient's airway, calculating a complex quantity representing a patient admittance (12) and comparing its value with ranges (14,16) indicative of open or closed apneas. The method distinguishes open from closed apneas even when the model used to calculate admittance is not based on details of the respiratory apparatus. In addition the patient admittance may be compared with admittance during normal breathing to avoid having to characterize the airway.

## Background
RELATED APPLICATION
This application claims priority of U.S. Provisional Application 60/823,973, filed Aug. 30, 2006, and U.S. Provisional Application 60/916,147, filed May 4, 2007, the specifications and drawings of which are incorporated herein by reference.


FIELD OF THE INVENTION
This invention relates to the discrimination of open and closed apneas (i.e. the complete cessation of breathing). In an open apnea the airway is patent, while in a closed apnea there is a total obstruction of the airway. The discrimination between such apneas is advantageous in the diagnosis and treatment of respiratory conditions that have adverse effects on a person's wellbeing.
BACKGROUND OF THE INVENTION
The expression “airway” as used herein is to be understood as the anatomical portion of the respiratory system between the nares and the bronchi, including the trachea. The expression “respiration” is to be understood as the continually repeating events of inspiration (inhaling) followed by expiration (exhaling).
In the Sleep Apnea syndrome a person stops breathing during sleep. Cessation of airflow for more than 10 seconds is called an “apnea”. Apneas lead to decreased blood oxygenation and thus to disruption of sleep. Apneas are traditionally categorized as either central, where there is no respiratory effort, or obstructive, where there is respiratory effort. With some central apneas, the airway is patent, and the subject is merely not attempting to breathe. Conversely, with other central apneas and all obstructive apneas, the airway is not patent (i.e. occluded). The occlusion is usually at the level of the tongue or soft palate. The airway may also be partially obstructed (i.e. narrowed or partially patent). This also leads to decreased ventilation (hypopnea), decreased blood oxygenation and disturbed sleep.
The dangers of obstructed breathing during sleep are well known in relation to the Obstructive Sleep Apnea (OSA) syndrome. Apnea, hypopnea and heavy snoring are recognized as causes of sleep disruption and risk factors in certain types of heart disease. Increased upper airway resistance (Upper Airway Resistance syndrome) during sleep without snoring or sleep apnea also can cause sleep fragmentation and daytime sleepiness.
The common form of treatment of these syndromes is the administering of Continuous Positive Airway Pressure (CPAP). Briefly stated, CPAP treatment acts as a pneumatic splint of the airway by the provision of a positive pressure, usually in the range 4-20 cm H2O. The air is supplied to the airway by a motor driven blower or other flow generator (FG) whose outlet passes via an air delivery hose to a nose (or nose and/or mouth) mask sealingly engaged to a patient's face. An exhaust port is provided in the delivery tube proximate to the mask. More sophisticated forms of CPAP, such as bi-level CPAP and autosetting CPAP, are described in U.S. Pat. Nos. 5,148,802 and 5,245,995 respectively.
As noted, central apneas need not involve an obstruction of the airway, and often occur during very light sleep and also in patients with various cardiac, cerebrovascular and endocrine conditions unrelated to the state of the upper airway. In those cases where the apnea occurs without obstruction of the airway, there is little benefit in treating the condition by techniques such as CPAP. In automated CPAP systems, it is important to accurately distinguish apneas with an open airway from apneas with a closed airway, in order to avoid inappropriately increasing the CPAP splinting air pressure. Such unnecessary increases in pressure reflexly inhibit breathing, further aggravating the breathing disorder.
U.S. Pat. No. 5,245,995 describes how snoring and abnormal breathing patterns can be detected by inspiration and expiration pressure measurements while sleeping, thereby leading to early indication of preobstructive episodes or other forms of breathing disorder. Particularly, patterns of respiratory parameters are monitored, and CPAP pressure is raised on the detection of pre-defined patterns to provide increased airway pressure to, ideally, prevent the occurrence of the obstructive episodes and the other forms of breathing disorder.
Prior Use of the Forced Oscillation Technique
U.S. Pat. No. 5,704,345, entitled “Detection Of Apnea And Obstruction Of The Airway In The Respiratory System” describes various techniques for sensing and detecting abnormal breathing patterns indicative of obstructed breathing, including the determination of airway patency by a forced oscillation technique (FOT) in which an oscillatory pressure waveform of known frequency is applied to a patient's airway and the magnitude of the component of an airflow signal at the known frequency induced by the oscillatory pressure waveform is calculated and compared with a threshold value. The present invention is an improvement of the methods and apparatus disclosed in the '345 patent.
The use of pressure oscillations at frequencies of the order of 4 Hz to determine airway patency was used in the ResMed AutoSet Clinical automatic CPAP device and the PII Plus. In these machines which used FOT, the pressure was measured at the mask, and the flow was measured very close to the mask, on the patient side of the mask vent. The present invention finds an advantage in measuring pressure and flow at or near the flow generatory at least for analyzing the effect of the forced oscillation technique.
The prior art implementations of FOT are less accurate in distinguishing between closed and open apneas when there is present moderate leak and moderate “resistance” in the airpath between the flow generator and the patient. For example, a passive patient simulation consisting of a 3-4 cm H2O/(l/s) resistance, with an adjustable leak, would indicate an open airway at a leak of 15 l/min and a closed airway at a leak of 20 l/min. What is needed is a system that more accurately distinguishes between open and closed apneas. In particular what is needed is a system that goes beyond treating the components of the airpath simply as nonlinear resistances and which utilizes an algorithm that takes into account the capacitive and inductive components of the airpath impedance.
BRIEF DESCRIPTION OF THE INVENTION
In a prior filed U.S. Provisional Patent Application Ser. No. 60/823,973, filed Aug. 30, 2006, it was disclosed that an effective method for discriminating between closed and open respiratory airway apneas was by determining the complex admittance of a patient airway and comparing the absolute value of the complex admittance to threshold values. That disclosure employed an algorithm that takes into account the patient circuit and models each component (patient circuit, vent flow and leak) in order to determine the patient pressure and flow. In particular, the small-signal hose pressure drop was modeled as a two port network in which the parameters were a function of mean flow. The present invention includes embodiments that are an improvement on that technique by including more information from that complex quantity, namely the phase angle of the admittance, or both the real and imaginary components of the complex quantity, in discriminating between open and closed apneas. In particular, it has been surprisingly discovered by analyzing data previously identified with open and closed apneas, that there is significantly better separation of the data if the phase angle of the admittance is taken into account. In addition, the previous system is simplified by modeling the pressure drop without having to measure the AC impedance of the patient circuit across a range of flows.
A plot of the magnitude of admittance against the phase of admittance provides distinct areas of the plot that indicate open or closed airways. Alternatively, the real part of the admittance may be plotted against the imaginary part in an Argand diagram, resulting again in distinct areas of the diagram being associated with open and closed apneas.
The invention provides an improved method and apparatus for treating a patient subject to apneas. In particular it relies upon improved methods for determining whether an apnea is closed or open, by identifying sets of values of the complex patient admittance that are respectively characteristic of open or closed apneas. The admittance is the reciprocal of the complex impedance of the apparatus parameterized for example as a two port model of flow impedance or determining resistive and inertial components from a theoretical model employing empirically derived resistive constants and theoretically determined flow inertance constants.
A sinusoidal (e.g. 4 Hz) pressure oscillation is applied at the input to the airpath, while flow and pressure is measured both at the input and output of the airpath. Based upon the two port model or the more theoretical hose drop model, the admittance is determined from the AC component of the patient airflow (found by subtracting the AC components of vent flow and leak) and the AC component of the mask pressure. The vent flow, in turn, is determined from an improved quadratic relationship to mask pressure in the case of the two port model, or a linearized calculation in the more theoretical hose drop model. The leak is determined from an estimated leak coefficient. The calculation of AC components of these quantities is improved over prior estimations by using Fourier sine and cosine components at the input oscillation frequency rather than by approximating an orthonormal set of functions by square waves.
The invention discloses a method for determining patency of the airway of a patient, the method comprising the steps of:
applying an oscillatory pressure waveform of known frequency to the patient's airway;
measuring respiratory air flow and pressure at the flow generator;
calculating the AC values of flow and pressure at the mask from a 2-port impedance model,
determining that the airway is patent by determining whether the complex admittance is in a region characteristic of patency.
Advantageously the admittance is determined from the ratio of AC values of patient flow and mask pressure, and there is the step of comparing the value of the complex admittance with ranges of values for which the airway is declared patent.
The invention yet further discloses a method for determining the degree of patency of the airway of a patient, the method comprising the steps of:
applying an oscillatory pressure waveform of known frequency, magnitude and phase at an entrance to the patient's airway;
measuring respiratory air flow from the patient;
determining the magnitude and phase of the component of said air flow at said known frequency induced by said oscillatory pressure waveform; and
determining the degree of patency as the location of the complex admittance within the region indicative of patency.
To minimize the effect of the delivery system not supplying a predetermined waveform in the aforementioned method, the pressure waveform actually produced is measured. In this technique we initially apply a waveform of some amplitude at the flow generator, calculate or observe the magnitude of the pressure waveform at the mask, then adjust (typically increasing) the amplitude at the flow generator in order to produce a desired amplitude at the mask. Given that the system is approximately linear for these small signals, the calculation as to how much to increase the driving waveform is a calculation of the radio of the desired to the actual pressure magnitude at the mask. The advantage of this approach is an improvement in the signal to noise ratio.
The invention yet further discloses a method for controlling the administration of CPAP treatment to the airway of a patient by means controllable to supply breathable air to the patient's airway continually at a selectable pressure elevated above atmospheric pressure, the method comprising the step of:
commencing or increasing CPAP treatment pressure if an apnea is occurring, determined by the steps of:
measuring respiratory air flow from the patient as a function of time; and
determining the deviation of said admittance from the centroid of the region of patency.
The invention yet further discloses apparatus for determining patency of the airway of a patient, the apparatus comprising:
means for applying an oscillatory pressure waveform of known frequency to the patient's airway;
means for measuring respiratory air flow from the patient; and
means for determining that the airway is patent if there is an admittance within a patency region at said known frequency induced by said oscillatory pressure waveform.
The invention yet further discloses apparatus for determining the degree of patency of the airway of a patient, the apparatus comprising:
means for applying an oscillatory pressure waveform of known frequency and magnitude to the patient's airway;
means for determining the complex admittance of respiratory air flow from the patient; and
means for determining the degree of patency as the deviation of said induced admittance from the centroid of a patency region.
The invention yet further discloses a method of distinguishing between open and closed airway apneas of a patient comprising the steps of:
(i) connecting a respiratory device to a patient via an air delivery tube and a patient interface;
(ii) delivering an alternating pressure waveform to the patient from the respiratory device to the patient via the air delivery tube;
(iii) measuring a flow rate and pressure of air at the respiratory device;
(iv) determining a capacitive component of an air delivery tube impedance;
(v) determining a patient admittance from said measured flow and pressure of air and said capacitive component;
(vi) distinguishing between an open and closed airway apnea on the basis of said patient admittance.
The invention yet further discloses method of distinguishing between open and closed airway apneas of a patient comprising the steps of:
(i) connecting a respiratory device to a patient via an air delivery tube and a patient interface;
(ii) delivering an alternating pressure waveform to the patient from the respiratory device to the patient via the air delivery tube;
(iii) measuring a flow rate and pressure of air at the respiratory device;
(iv) determining an inductive component of an air delivery tube impedance;
(v) determining a patient admittance from said measured flow and pressure of air and said capacitive component;
(vi) distinguishing between an open and closed airway apnea on the basis of said patient admittance.
In another form of the invention patency of the airway is determined to be unknown, or in a “don't know” region. Another aspect of the invention is that patency is defined as having a gradual scale from open to closed. Another aspect is that apnea discrimination regions are defined in the complex plane, the regions not necessarily being contiguous. A further refinement is that the discrimination regions may be functions of the average leak level; in particular, the don't-know regions probably increase in size with increasing leak
Further forms of the invention are as set out in the claims.



BRIEF DESCRIPTION OF THE DRAWINGS
Embodiments of the invention will now be described with reference to the accompanying drawings, in which:
FIG. 1 shows a flow diagram of the basic methodology of an embodiment;
FIG. 2 shows, in diagrammatic form, apparatus embodying the invention;
FIG. 3 shows an electronic analogue of the 2-port analysis of the present invention.
FIG. 4 shows values of the complex admittance based on a detailed algorithm.
FIG. 5 shows values of the complex admittance based on a simplified algorithm.



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS AND BEST MODE
A System Implementing the Invention
The present invention is an improvement upon the embodiments disclosed in U.S. Pat. No. 5,704,345, which is incorporated by reference in its entirety. FIG. 1 is a flow diagram of the basic methodology of one embodiment of the present invention. The first step 10 of the present invention is the measurement of respiratory flow and pressure at points near the flow generator over time where an apnea is occurring. This information is processed in step 12 to generate admittance values to be used as qualitative measures for subsequent processing. Steps 14-16 detect whether a closed, open or mixed apnea is occurring by comparison of the complex admittance value in a time window with a patency region.
If an apnea is in progress there then follows a determination whether the apnea is open or closed. If an apnea with an open airway is occurring, and, if desired, the event is logged in step 18. If the result of step 16 is that an apnea with a closed airway is occurring, an increase in CPAP treatment pressure occurs in step 20. If desired, step 20 may include the optional logging of the detected abnormality.
In the instance of an apnea with an open airway the CPAP treatment pressure is reduced, in accordance with usual methodologies that seek to set the minimal pressure required to obviate, or at least reduce, the occurrence of apneas. The amount of reduction in step 22 may, if desired, be zero.
The methodology represented in FIG. 1 is of a clinical embodiment, where patient CPAP pressure is controlled over time as appropriate. A purely diagnostic embodiment operates in the same manner except it omits the CPAP pressure increase and pressure decrease actions of step 20 and step 22 respectively.
FIG. 2 shows, in diagrammatic form, clinical CPAP apparatus in accordance with one embodiment for implementing the methodology of FIG. 1. A mask 30, whether either a nose mask and/or a face mask, is sealingly fitted to a patient's face. Fresh air, or oxygen enriched air, enters the mask 30 by flexible tubing 32 which, in turn, is connected with a motor driven turbine (flow generator) 34 to which there is provided an air inlet 36. The motor 38 for the turbine is controlled by a motor-servo unit 40 to either increase or decrease the pressure of air supplied to the mask 30 as CPAP treatment. The mask 30 also includes an exhaust port 42 that is close to the junction of the tubing 34 with the mask 30.
Adjacent to the flow generator 34 is a flow-resistive element 44. This can take the form of an iris across which a differential pressure exists. The mask side of the flow-resistive element 44 is connected by a small bore tube 46 to a pressure transducer 48 and to an input of a differential pressure transducer 60. Pressure at the other side of the flow-resistive element 44 is conveyed to the other input of the differential pressure transducer 50 by another small bore tube 52.
The pressure transducer 48 generates an electrical signal in proportion to the flow pressure, which is amplified by amplifier 56 and passed both to a multiplexer/ADC unit 58 and to the motor-servo unit 40. The function of the signal provided to the motor-servo unit 40 is as a form of feedback to ensure that the static pressure is controlled to be closely approximate to the set point pressure.
The differential pressure sensed across the flow-resistive element 44 is output as an electrical signal from the differential pressure transducer 50, and amplified by another amplifier 60. The output signal from the amplifier 56 therefore represents a measure of the mask or respiratory airflow rate. The controller 62 is programmed to perform a number of processing functions.
The pressure and flow may be considered to be composed of steady state values and AC values, the latter reflecting the effect of a imposed oscillatory signal on the pressure having a frequency of 4 Hz. In what follows, a “steady-state” quantity is either (a) one from which the oscillatory component has been deliberately removed, for example by a filtering operation, or (b) is the result of a calculation based partly or wholly on quantities from which the oscillatory component has been removed, or (c) is a quantity which is calculated based on instantaneous quantities, such as pressure and flow measured at the flow generator (which thus include an oscillatory component), and a model of the airpath and patient leak which either partly or wholly ignores the reactive components of the system and treats it simply as a (possibly nonlinear) system of resistances. The steady state pressure loss along tubing 32 is calculated from the flow through the tube, and knowledge of the static pressure-flow characteristic of the tubing, for example by table lookup. The steady state pressure at the mask is then calculated by subtracting the tube pressure loss. The pressure loss along tube 32 is then added to the desired set pressure at the mask to yield the desired instantaneous pressure at the pressure generator 34. The flow through the exhaust 42 is calculated from the pressure at the mask from the pressure-flow characteristic of the exhaust, for example by table lookup. The steady state mask flow is calculated by subtracting the flow through the exhaust 42 from the flow through the tubing 32. The steady state patient flow is then calculated by subtracting the steady-state estimated leak, which may be determined for example in a CPAP device by a 1st order low pass filter with a time constant of 10 seconds whose input is the instantaneous mask flow, from the steady-state mask flow.
The methodology put into place by the controller 62 will now be described with reference to the apparatus of FIG. 2. If the patient respiratory flow is very low or zero (note that the mask flow will not cease when the patient is apnoeic if there is any leak), a determination of airway patency (steps 14-16) is made by using an externally induced oscillation technique. If the airway is open, but the respiratory muscles are relaxed (i.e. a central apnea with open airway), then small externally originating fluctuations in the mask pressure will induce a small respiratory airflow by inflating and deflating the lungs, and by compressing and decompressing the gas in the lungs. Conversely, if the airway is closed, no airflow will be induced. This is quantified as follows:
Discriminating Apneas by Admittance Thresholds
The admittance Y is given by

Y=G+iB 

where G is conductance and B is susceptance.
In order to decide whether the airway is open, the complex value of the patient admittance, Y, is compared with a region of values. The value of this threshold may be selected on the basis of the following observations.
An explanation for the angle being a better classifier than the magnitude is possibly related to the interaction of the complex impedances in the circuit. The leak and vent flow are essentially non-linear resistances with no complex (j) component. The patient is basically a resistance and compliance (capacitance) in series to ground. The hose drop also has a complex component in the form of the inertance of the flow. So the magnitude of the overall impedance is influenced by all components (hose, leak, vent flow, patient airway) as they all have real parts. However the imaginary component is the interaction between the inertial hose drop and the capacitive lungs. It is comparing the relative “strengths” of the inertial and capacitive reactances.
It is also possible to use the real and imaginary parts of the impedance rather than the angle and magnitude. This has an advantage of reduced CPU load. The magnitude calculation requires a square root and the angle calculation involves an inverse tan. (Both of these can be done to easily the required accuracy by techniques involving lookup tables, which are computationally fairly cheap.) Using both the real and imaginary or magnitude and angle gives a more robust classification. The grouping of the open and closed values becomes much more apparent.
As determined from the values in FIG. 4, a possible set of thresholds would be














angle = ARG(patient_admittance)



if ( angle < −2 OR angle > 1.2 )









AirwayState = CLOSED;









else if ( angle < 0.9 AND angle > −1.5 )









AirwayState = OPEN;









else









AirwayState = UNKNOWN;









end if









    
    
        where angle is in radians. The acceptable values for the thresholds may be any that permits the separation of values representing open and closed apneas determined from FIG. 4.
    
    


The patient resistance at 4 Hz also indicates the state of the airway. This resistance is not the reciprocal of the conductance in the above equation. The patient impedance at 4 Hz is

Z=R+iX 

where




  Y
  =
  
    
      1
      Z
    
    .
  




The patient resistance R is given by





  R
  =
  
    
      Re
      ⁡
      
        (
        Z
        )
      
    
    =
    
      Re
      ⁡
      
        (
        
          1
          Y
        
        )
      
    
  




In an alternative embodiment R may be used to characterize the state of the airway rather than |Y|.
Calculating Admittance
Patient admittance is calculated by the equation




  
    
      Y
      pat
    
    =
    
      
        Q
        
          pat
          ,
          
            
          
          ⁢
          AC
        
      
      
        P
        
          mask
          ,
          AC
        
      
    
  
  ,




    
    
        where the numerator is the magnitude of the AC patient flow, which is in effect a measure of the differential (with respect to time) of patient flow. The denominator is the same differential of the mask pressure.
    
    


Correcting for Vent Flow and Leak
The patient airflow is determined by subtracting flow through the vent and leakage flow from the inflow to the mask.

Qpat,AC=Qmaskin,AC−Qvent,AC−Qleak,AC 

The mean vent flow is determined either from mean mask pressure during the period of measurement, or more preferably from the square of the mean of the square root of the mask pressure during this time (which is also used to make leak calculations). The mask pressure itself is calculated in the conventional way, ignoring AC behavior of the airpath.
AC vent flow is calculated using a linear approximation about the operating point:




  
    Q
    
      vent
      ,
      AC
    
  
  =
  
    
      P
      
        mask
        ,
        AC
      
    
    ⁢
    
      
        d
        
          Q
          vent
        
      
      
        d
        
          P
          mask
        
      
    
  




Modeling Mask Pressure to Determine AC Flow at Mask
The mask pressure is determined by evaluating the coefficients k1 and k2 in the equation which gives the pressure drop across the vent as a function of vent flow:

Pmask=k1Qvent+k2Qvent2.

    
    
        Thus to calculate the vent flow at a particular mask pressure, one solves the quadratic equation for vent flow at that mask pressure.
    
    






  
    d
    
      Q
      vent
    
  
  
    d
    
      P
      mask
    
  




can be obtained by differentiation of the previous equation, which gives





  
    
      d
      
        P
        mask
      
    
    
      d
      Q
    
  
  =
  
    
      k
      1
    
    +
    
      2
      ⁢
      
        
      
      ⁢
      
        k
        2
      
      ⁢
      
        Q
        vent
      
    
  




Since




  
    
      
        Δ
        ⁢
        
          
        
        ⁢
        Q
      
      
        Δ
        ⁢
        
          
        
        ⁢
        P
      
    
    ≈
    
      
        d
        Q
      
      
        d
        P
      
    
  
  ,




a small-signal approximation yields





  
    
      Δ
      ⁢
      
        
      
      ⁢
      Q
    
    ≈
    
      Δ
      ⁢
      
        
      
      ⁢
      P
      ⁢
      
        
          d
          Q
        
        
          d
          P
        
      
    
  
  ,




from which it follows as a reasonable approximation that





  
    Q
    
      vent
      ,
      AC
    
  
  =
  
    
      
        P
        
          mask
          ,
          AC
        
      
      
        
          k
          1
        
        +
        
          2
          ⁢
          
            
          
          ⁢
          
            k
            2
          
          ⁢
          
            Q
            vent
          
        
      
    
    .
  




Calculating Leakage Flow
To complete the calculation of the patient airflow it is necessary to also calculate the leakage flow from the inflow to the mask. The actual leak is calculated over the time period during which the patient 4 Hz admittance is calculated. Typical leak coefficient estimates relate to what happened some time ago, and one may, as in the prior art AutoSet CS and AutoVPAP devices estimate the leak coefficient Kleak from




  
    K
    leak
  
  =
  
    
      
        Q
        nonvent
      
      _
    
    
      
        
          P
          mask
        
      
      _
    
  




    
    
        where the overbars indicate the mean over the measurement period. During apnea the leak is equal to the total non-vent flow

Qnonvent=QFGTotal−Qvent 

        and during breathing, the leak is (to a good approximation) equal to the average non-vent flow. In the above estimation of Kleak, all quantities are calculated without reference to the AC characteristics of the airpath, based on instantaneous values, treating the airpath as a nonlinear resistor.
    
    


The model for instantaneous leak, both DC (on which the above formula for Kleak is based) and AC is

Qleak=Kleak√{square root over (Pmask)}

    
    
        and, as with vent flow, the AC component (at 4 Hz) of leak flow is found using the small-signal approximation
    
    






  
    Q
    
      leak
      ,
      AC
    
  
  =
  
    
      P
      
        mask
        ,
        AC
      
    
    ⁢
    
      
        d
        
          Q
          leak
        
      
      
        d
        
          P
          mask
        
      
    
  




    
    
        where we may calculate
    
    






  
    d
    
      Q
      leak
    
  
  
    d
    
      P
      mask
    
  




by differentiating the above equation directly, giving





  
    
      d
      
        Q
        leak
      
    
    
      d
      
        P
        mask
      
    
  
  =
  
    
      
        K
        leak
      
      
        2
        ⁢
        
          
            P
            mask
          
        
      
    
    .
  




The mean for √{square root over (Rmask)}, may be either √{square root over ( Pmask or √{square root over (Pmask)}. The latter is preferred (because it gives an unbiased estimate of the mean value of




  
    d
    
      Q
      leak
    
  
  
    d
    
      P
      mask
    
  




when Pmask is not constant, and we have already calculated it in order to estimate Kleak), yielding





  
    
      d
      
        Q
        leak
      
    
    
      d
      
        P
        mask
      
    
  
  =
  
    
      
        
          Q
          nonvent
        
        _
      
      
        2
        ⁢
        
          
            (
            
              
                
                  
                    P
                    mask
                  
                  ⁢
                  
                    
                  
                
              
              _
            
            )
          
          2
        
      
    
    .
  




Determination of AC Values
Throughout this description, the AC values referred to are determined from measurements as trigonometric Fourier components at the exciting frequency. This is an improvement over the prior art use of the forced oscillation technique, which used square waves as an orthonormal set of functions.
All AC quantities are calculated as complex numbers. In particular, the AC pressure and flow at the mask are calculated by finding the sine and cosine components of instantaneous pressure and flow oven the period of measurement using the inner products with the sine and cosine function respectively, yielding component coefficients cs and cc, then writing either cs+icc or cc+ics as the AC value. (Either of these two forms is used, the two forms differing merely in a phase shift, but the one form is used consistently throughout). Alternatively other standard methods of estimating the amplitude and phase of the sinusoidal component of the pressure and flow at the exciting frequency may be used, such as least-squares fitting of a sine and a cosine.
Explicitly we have



  
    
      
        c
        s
      
      =
      
        
          1
          T
        
        ⁢
        
          
            ∫
            0
            
              2
              ⁢
              
                π
                /
                w
              
            
          
          ⁢
          
            
              f
              ⁡
              
                (
                t
                )
              
            
            ⁢
            
              sin
              ⁡
              
                (
                
                  ω
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
                )
              
            
            ⁢
            
              d
              t
            
          
        
      
    
    ;
  




  
    
      
        c
        C
      
      =
      
        
          1
          T
        
        ⁢
        
          
            ∫
            0
            
              2
              ⁢
              
                π
                /
                w
              
            
          
          ⁢
          
            
              f
              ⁡
              
                (
                t
                )
              
            
            ⁢
            
              cos
              ⁡
              
                (
                
                  ω
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
                )
              
            
            ⁢
            
              d
              t
            
          
        
      
    
    ;
  



    
    
        where the exciting frequency is T/2B, T is a normalizing factor, and the integral may be replaced by a summation using discrete sampled values of sine and cosine. f(t) is any function whose Fourier coefficients are required.
    
    


All quantities are determined over a 6 second sliding window (the “admittance window”) symbolically indicated as t=0. An admittance calculation can in principle be made at the algorithmic sampling frequency (say 50 Hz), but for clinical purposes this is not essential, and the calculation may be efficiently performed at the end of every 4 Hz cycle. Alternatively, calculation of admittance at 2 Hz or even 1 Hz frequencies is reasonable for clinical purposes.
Due to motor controller delays, the first cycle (250 ms) of the pressure waveform may not be sinusoidal, and there will be some delay in setting up a steady state in the airpath, so the first 250 ms of data may be ignored.
The sine and cosine values would be stored in a table generated at startup or specified as constants in the source code. Sufficient accuracy is provided if multiplications are fixed point, 16 bit*16 bit with a 32 bit result, assuming of course that overflow of the sum does not occur, and there are no numerical stability issues with this approach. Thus the computational cost is relatively small.
Modeling Impedance on a Two Port Network
The calculation of AC pressure at the mask and AC flow entering the mask is in one embodiment determined from a 2-port electrical network analogue of the flow in the system at a particular frequency.
From FIG. 3, routine circuit analysis gives



  
    
      V
      out
    
    =
    
      
        V
        
          i
          ⁢
          
            
          
          ⁢
          n
        
      
      -
      
        
          I
          
            i
            ⁢
            
              
            
            ⁢
            n
          
        
        ⁢
        
          Z
          inout
        
      
    
  




  
    
      I
      out
    
    =
    
      
        I
        
          i
          ⁢
          
            
          
          ⁢
          n
        
      
      -
      
        
          V
          out
        
        
          Z
          outgnd
        
      
    
  



    
    
        and of course these correspond to pressure and the mask and flow entering the mask respectively. Explicitly:
    
    





  
    
      P
      
        mask
        ,
        
          A
          ⁢
          
            
          
          ⁢
          C
        
      
    
    =
    
      
        P
        
          FG
          ,
          
            A
            ⁢
            
              
            
            ⁢
            C
          
        
      
      -
      
        
          Q
          
            FG
            ,
            
              A
              ⁢
              
                
              
              ⁢
              C
            
          
        
        ⁢
        
          Z
          inout
        
      
    
  




  
    
      Q
      
        maskin
        ,
        
          A
          ⁢
          
            
          
          ⁢
          C
        
      
    
    =
    
      
        Q
        
          FG
          ,
          
            A
            ⁢
            
              
            
            ⁢
            C
          
        
      
      -
      
        
          P
          
            mask
            ,
            
              A
              ⁢
              
                
              
              ⁢
              C
            
          
        
        
          Z
          outgnd
        
      
    
  



The impedances Zinout and Zoutgnd are taken to be functions of average flow. These impedances are measured for a particular airpath at a number of average flow levels (see below for details). The average total flow generator flow during the admittance window is determined, and for each of Zinout and Zoutgnd, is used to linearly interpolate between the measured impedances to determine the impedance at that average flow level.
The inlet and outlet AC flows and AC pressures are determined over a period of 30 seconds (to reduce noise) using the standard inner product method described above for the calculation of Fourier coefficients, yielding complex values. These values are used to calculate the 2-port parameters (refer to FIG. 1, where voltages and currents in that figure correspond to pressures and flows respectively in the following), by:



  
    
      Z
      inout
    
    =
    
      
        
          P
          
            i
            ⁢
            
              
            
            ⁢
            n
          
        
        -
        
          P
          out
        
      
      
        Q
        
          
            
          
          ⁢
          
            i
            ⁢
            
              
            
            ⁢
            n
          
        
      
    
  




  
    
      Z
      outgnd
    
    =
    
      
        P
        out
      
      
        
          Q
          
            
              
            
            ⁢
            
              i
              ⁢
              
                
              
              ⁢
              n
            
          
        
        -
        
          Q
          out
        
      
    
  



    
    
        If the denominators are very small or zero, the impedances are taken to have some numerically very large value in relation to typical airpath impedances.
    
    


The Hose Pressure Drop Theoretical Model
The hose drop is the pressure difference between the internal mask pressure and the FG pressure. It is the pressure loss across all components that are placed before the mask. This includes mufflers, humidifiers, AB filters and the mask connection hose.
The complete hose drop is made up of the resistive hose drop as well as the inertial pressure drop




  
    P
    hosedrop
  
  =
  
    
      
        K
        1
      
      ⁢
      
        Q
        FG
        2
      
    
    +
    
      
        K
        2
      
      ⁢
      
        Q
        FG
      
    
    +
    
      
        K
        L
      
      ⁢
      
        
          d
          
            Q
            FG
          
        
        
          d
          t
        
      
    
  




K1,K2=empirically derived constants
KL=theoretically determined flow inertance constant
QFG=FG Flow




  
    
      K
      L
    
    =
    
      
        ρ
        ⁢
        
          
        
        ⁢
        l
      
      A
    
  
  ,




ρ=air density (1.19 kg/m3)
l=tube length (2 or 3 m)
A=cross sectional area




  (
  
    
      
        π
        ⁢
        
          
        
        ⁢
        
          d
          2
        
      
      4
    
    ,
    
      d
      =
      
        0.019
        ⁢
        
          
        
        ⁢
        m
      
    
  
  )




In the case of masks such as the ResMed Activa and Swift masks there is another tube between the main hose and the mask. This must also be taken into account in the inertial constant.
The mean hose drop is made up of purely the resistive part (the over bar indicates mean):

 Phosedrop=K1 QFG+K2 QFG 

K1,K2=empirically derived constants
 QFG=Mean FG Flow
In order to calculate the AC hose drop some linearization is necessary. The resistive and inertial components are separated. The resistive component is derived by linearizing the above mean hose drop formula about the operating point.




  
    
      d
      
        P
        hosedrop
      
    
    
      d
      
        Q
        FG
      
    
  
  =
  
    
      2
      ⁢
      
        K
        
          1
          ⁢
          
            
          
        
      
      ⁢
      
        Q
        FG
      
    
    +
    
      K
      2
    
  




Phosedrop resisitive,AC=(2K1 QFG+K2)QFG,AC 

The AC hose drop is thus the AC flow multiplied by gradient of the hose non-linear resistance at the location around which the small oscillation occurs. For the derivation it is assumed that the oscillations are of a small magnitude compared to the change in the gradient of the hose resistance.
The inertial component is a constant multiplied by the time derivative of the FG flow. By standard linear circuit theory, an inductance L has impedance sL, which for sinusoidal signals at the frequency ω is jωL. The inertial hose drop component then becomes:

Phosedrop iner tan ce,AC=KLjωQFG,AC 

ω=2πf

Combining the two gives

Pdrop,AC=(2K1 QFG+K2+KLjω)QFG,AC 

Mask Pressure
The mask pressure is the FG pressure minus the hose pressure drop:

Pmask=PFG−Phose drop 

Pmask,AC=PFG,AC−Phose drop,AC 

EXAMPLES
In FIG. 4 a plot is presented in which the thick lines represent open apneas and the thin lines represent closed apneas. The horizontal axis represents the absolute magnitude of the admittance, while the vertical axis represents the phase of the complex admittance. As may be seen from the figure, there is a complete separation of the open apneas, represented by the thick curves in the center of the graph, from the closed apneas, represented by the two thinner curve regions, one near the origin and the other at the extreme upper value of the y axis. If one looks only at the absolute value of the admittance, in effect by projecting all of the data onto the x axis, the separation of the data is lost, with the projection of the thin (closed apnea) curves meeting the projection of the thick (open apnea) data.
FIG. 4 shows the result of plotting patient admittance data calculated according to the formulas above. The x axis values are the average absolute value of the patient admittance. The y axis values are the phase angle of the patient admittance. As may be seen from FIG. 4 there is a clean separation of the values for closed apneas and for open apneas.
FIG. 5 shows the result of plotting patient admittance where instead of the previous calculation of admittance a value for the complex admittance is determined from a model that is essentially independent of the modelling details of any particular configuration of mask or ventilation circuit. In particular, the data for FIG. 5 were determined by using a simple leak orifice model P=(KQ)2. It may be seen from FIG. 5 that the open and closed apneas result in separated regions of the plot of the complex admittance.
Moderate (Approximated Patient Admittance) Algorithm
The third algorithm attempts to approximate the patient admittance by generically modeling the patient circuit and combining the leak and vent flows.
The patient circuit is modeled in the same way as the Complex algorithm above, however the constants are fixed at those for a 2 m hose with AB filter.
The leak and vent flow modeling are combined into one parameter which is the non patient flow. This is modeled as
Leak+Vent Flow
The leak and vent flow are combined and modeled as using the following equations:
The leak orifice can be approximated at runtime using:




  
    K
    leak
  
  =
  
    
      
        Q
        FG
      
      _
    
    
      
        
          P
          mask
        
      
      _
    
  




The instantaneous leak is then calculated using

Qleak=Kleak√{square root over (Pmask)}

Using small signal approximation we can say that




  
    Q
    
      leak
      ,
      
        A
        ⁢
        
          
        
        ⁢
        C
      
    
  
  =
  
    
      P
      
        mask
        ,
        
          A
          ⁢
          
            
          
          ⁢
          C
        
      
    
    ⁢
    
      
        d
        
          Q
          leak
        
      
      
        d
        
          P
          mask
        
      
    
  




where





  
    
      d
      
        Q
        leak
      
    
    
      
        d
        
          P
          mask
        
      
      ⁢
      
        
      
    
  
  =
  
    
      K
      leak
    
    
      2
      ⁢
      
        
          P
          mask
        
      
    
  




(from instantaneous leak equation)

substituting in Kleak we get




  
    
      Q
      
        leak
        ,
        
          A
          ⁢
          
            
          
          ⁢
          C
        
      
    
    =
    
      
        P
        
          mask
          ,
          
            A
            ⁢
            
              
            
            ⁢
            C
          
        
      
      ⁢
      
        
          
            Q
            FG
          
          _
        
        
          2
          ⁢
          
            
              (
              
                
                  
                    P
                    mask
                  
                
                _
              
              )
            
            2
          
        
      
    
  
  ,




as we are dealing with small mean Pmask changes we can approximate ( √{square root over (Pmask)})2 with Pmask so the AC leak becomes





  
    Q
    
      leak
      ,
      
        A
        ⁢
        
          
        
        ⁢
        C
      
    
  
  =
  
    
      P
      
        mask
        ,
        
          A
          ⁢
          
            
          
          ⁢
          C
        
      
    
    ⁢
    
      
        
          Q
          FG
        
        _
      
      
        2
        ⁢
        
          
            P
            mask
          
          _
        
      
    
  




Approximated Patient Flow
The patient flow is then

Qpatient,AC=QFG,AC−Qleak,AC 

Qleak is the leak flow combined with the vent flow (all non patient flow).
Classification
There are a number of ways to use the approximated admittance to classify the airway state. Good separation in just the angle can be seen when plotting the angle vs magnitude of the admittance. Therefore a simple classification would be to have two thresholds, but purely on angle. However the calculation of the angle involves a division and an inverse tan, requiring more calculation. It is suggested that it is preferable to have the real and imaginary parts of the admittance used instead.
Thresholds
The simplest thresholds which provide good classification accuracy are simple straight diagonal lines. These are of the form y>x+b and y<x+c. In this case when plotting the imaginary part (y) vs the real part (x) of the admittance, CLOSED is y>x+b, and OPEN is y<x+c. Based on the floating-point algorithm the values for the constants were b=0, c=−0.03. This gave no misclassification. Any other curves that separate the data regions could be used, but the linear separation reduces the complication of the threshold calculations and is for that reason preferred.

## Claims

1. A method of classifying airway apneas of a patient as open and closed with a positive airway pressure respiratory apparatus connectable to a patient via an air delivery tube and a patient interface, the method comprising:
(i) controlling delivery of a pressure waveform to a patient from the respiratory apparatus via the air delivery tube;
(ii) measuring a flow rate and pressure of air at the respiratory apparatus;
(iii) determining a capacitive component of an air delivery tube impedance with a processor of the respiratory apparatus;
(iv) determining a patient admittance from said measured flow and pressure of air and said capacitive component with the processor; and
(v) classifying the airway apnea on a basis of said patient admittance with the processor.

2. The method of claim 1 further comprising correcting, with the processor, said flow rate for leak.

3. The method of claim 1 further comprising correcting, with the processor, the said flow rate for vent flow.

4. A method of classifying airway apneas of a patient as open and closed comprising:
(i) delivering a pressure waveform to the patient from a positive airway pressure respiratory apparatus to a patient via an air delivery tube;
(ii) measuring a flow rate and pressure of air at the respiratory apparatus;
(iii) determining an inductive component of an air delivery tube impedance;
(iv) determining a patient admittance from said measured flow and pressure of air and said inductive component;
(v) classifying the airway apnea on the basis of said patient admittance in the respiratory apparatus.

5. The method of claim 4 further comprising the step of correcting the said flow rate for leak.

6. The method of claim 4 further comprising the step of correcting the said flow rate for vent flow.

7. The method of claim 1 further comprising determining, with the processor, the phase angle of the admittance and classifying the airway apnea based on the phase angle.

8. The method of claim 1 further comprising determining the magnitude of the admittance and classifying the airway apnea based on the magnitude of the admittance.

9. The method of claim 1 further comprising determining the real component of the admittance and classifying the airway apnea based on the real component of the admittance.

10. The method of claim 1 further comprising determining the imaginary component of the admittance and classifying the airway apnea based on the imaginary component of the admittance.

11. The method of claim 1 whereby said pressure waveform is an alternating pressure waveform having a frequency of approximately 4 Hz.

12. The method of claim 1 whereby said pressure waveform is an alternating pressure waveform having a frequency in the range of approximately 2-8 Hz.

13. The method of claim 1 whereby said pressure waveform is an alternating pressure waveform having a frequency in the range of approximately 1-16 Hz.

14. The method of claim 1 whereby approximately the first 250 ms of flow or pressure data following an application of the pressure waveform are not measured.

15. The method of claim 1 wherein classifying the airway apnea comprises plotting a magnitude of the admittance against a phase angle of the admittance.