The invention relates generally to a system for performing a system-specific analysis of turbomachinery such as detecting and isolating faults in turbomachinery and, more particularly, to a system for detecting and isolating faults in an aircraft engine.
Various types of turbomachinery are known and are generally in use in a range of applications, such as jet engines, industrial gas turbines, steam turbines and so forth. Typically, the components of the turbomachinery may be subjected to general wear and tear during their lifetime. In addition, the components may be exposed to abnormal conditions while in operation. As a result, the components of the turbomachinery can deteriorate, fail and lead to faults and inefficient operation of the turbomachinery. Consequently, it may be desirable to detect and isolate such faults in the engine for tracking engine health to ensure efficient operation of the engine.
Many specific techniques have been developed for detecting faults in the components of the engines and other systems. For example, in some systems data related to specific engine parameters is collected over a period of time and such data is analyzed and used for predicting engine outputs. Typically, models of jet engine performance are used to predict engine outputs based upon the collected input data. In addition, sensors are employed for measuring engine outputs for various components of the engine. Further, the differences between the predicted engine outputs and measured engine outputs are used for detecting and isolating faults in the engine.
In general, the engine model is based upon a sample of fleet of engines and there may be a significant variation in such model from one engine to another in the fleet of engines. In addition, there is variation in the input data across the fleet of engines. Such variation in the input data across the fleet of engines and the model errors causes scatter in the differences between the predicted and measured engine outputs that limits the performance of fault detection and isolation of faults.
Accordingly, it would be desirable to develop a system to detect and isolate faults in turbomachinery in a more efficient manner. More specifically, it would be desirable to have an efficient fault detection system for an aircraft engine that reduces the scatter in the engine outputs across a fleet of engines and thereby provides efficient fault detection and isolation along with a relatively accurate estimate of an engine component health for a particular engine over a period of time.
Briefly, in accordance with one aspect of the present invention a method for system-specific analysis of an engine includes applying control inputs to the engine and an engine model and estimating outputs from the engine model based upon the control inputs. The method includes sensing outputs from the engine and analyzing residuals between estimated and sensed outputs. The method also includes customizing the engine model to reduce residuals for a particular engine and detecting the faults in the engine based upon the residuals for the particular engine. Computer-readable medium that afford functionality of the type defined by this method is also provided by the present technique.
In accordance with another aspect of the present invention a system for detecting faults in an engine includes an engine model configured to receive control inputs corresponding to the engine control inputs and sensed inputs and to estimate outputs based upon the control inputs and the sensed inputs. The system also includes a plurality of sensors configured to sense outputs from the engine and an estimator configured to customize the engine model to reduce residuals between the estimated and sensed outputs.
These and other features, aspects, and advantages of the present invention will become better understood when the following detailed description is read with reference to the accompanying drawings in which like characters represent like parts throughout the drawings, wherein:
As discussed in detail below, embodiments of the present technique function to detect and isolate faults in turbomachinery such as an aircraft engine, an industrial gas turbine and a steam turbine. Turning now to drawings and referring first to
In this embodiment, the inputs from the PLA 18 are processed through a closed loop control 20 to generate control inputs 22 (u) for the engine 12 and the engine model 14. In one embodiment, the engine model 14 includes a physics-based model. In another embodiment, the engine model 14 includes an empirical model. Further, the engine model 14 may include a steady state model. Alternatively, the engine model 14 may include a transient model. In the illustrated embodiment, the control inputs 22 may include, but are necessarily not limited to, a fuel flow, an active clearance control, variable geometry, power extraction and combinations thereof for components of the engine 12. Typically, the components of the engine 12 include a fan, a booster, a high-pressure compressor, a low-pressure compressor, a high-pressure turbine, a low-pressure turbine and a combustor, among others. Other or different components and parameters may, of course, be monitored and controlled by the present techniques, depending upon the aircraft type, its equipment, and the control regimes envisaged. In addition, the engine model 14 may receive sensed inputs 23 such as, but not necessarily limited to, temperature, pressure, altitude, Mach number, or combinations thereof. Further, as will be appreciated by one skilled in the art such sensed inputs 23 will typically be experienced by the engine 12 in operation. Of course, these sensed inputs 23 are considered inputs for the engine model 14.
In the illustrated embodiment, the components of the engine 12 operate based upon the control inputs 22. Further, a plurality of sensors (not shown) may be coupled to each of the components of the engine 12 for sensing outputs from the engine 12. In certain embodiments, the sensed outputs from engine may include noise components due to factors such as random variation 24 (w) and sensor errors 26 (v) for the plurality of sensors coupled to the components of the engine 12. In operation, the plurality of sensors measure sensed outputs 28 (y) of the components. Examples of sensed outputs 28 include temperature, pressure, rotor speed, efficiency, flow capacity, inter-component temperature and so forth.
In a presently contemplated configuration, the engine model 14 generates predicted sensor outputs 30 (ŷ) based upon the control inputs 22. The predicted sensor outputs 30 from the engine model 14 do not include any noise components due to the random variation 24 and the sensor errors 26. The sensed outputs 28 from the plurality of sensors are combined and compared with the predicted senor outputs 30 as represented by reference numeral 32 to estimate residuals 34 (v). Further, the estimated residuals 34 may be analyzed to detect and isolate faults in the engine 12 by comparing the estimated residuals 34 with fault signatures via a fault diagnostics system 36 that will be described below. In a presently contemplated configuration the fault diagnostics system 36 is a part of the FADEC 16. However, those skilled in the art will appreciate that the fault diagnostics system 36 may be isolated from the FADEC 16 or other control systems. This is particularly true in aircraft where no FADEC 16 is present. In certain embodiments, the fault diagnostics system 36 may be partitioned within the FADEC 16 from the other control modules. In one embodiment, the estimated residuals 34 may be analyzed in real time on-wing. Alternatively, the estimated residuals 34 may be analyzed at a diagnostic location on ground, either in real time, near real time, or at a later time.
In accordance with the present techniques, the engine model 14 may be customized to reduce residuals 34 between the sensed outputs 28 and the predicted sensor outputs 30 for the particular engine 12. Further, the residuals 34 for a customized engine model 14 may function to reduce model errors and errors due to noise such as random variation 24 and sensor errors 26, thereby providing a substantially accurate detection of faults in the engine 12. It should be noted that it is the ability to correct for model errors that permits system-specific customization of the engine model 14. The engine model 14 may be customized to match the particular engine 12 by coupling an estimator to the engine model 14 that will be described below with reference to
In operation, the estimated parameters from the estimator 42 are applied to the engine model 14 to update the parameters of the engine model 14 for reducing the residuals 34. It should be noted that in a presently contemplated embodiment, the parameters of the engine model 14 are updated at a bandwidth sufficiently fast to track changes in the engine 12 and sufficiently slow to avoid masking faults occurring in the engine 12 to avoid the engine model 14 adapting to faults in the engine 12 that are otherwise required to be detected for an efficient operation of the engine 12 (i.e., customizing the model to undesired conditions). Advantageously, the residuals 34 from the customized engine model 14 may be utilized for detecting and isolating the faults in the engine 12 via the fault diagnostics system 36. In this embodiment, the tracking filter 46 analyzes the residuals (sometimes referred to as “innovations”) 34 to estimate the engine parameters that will be described below with reference to
In the illustrated mathematical representation, a non-linear function 70 (h(xk)) relates the estimated state 68 of the system to measurements from the system 54. In this embodiment, measurement noise 72 (vk) such as due to sensor errors may be incorporated into the function 70 of the system 54 as represented by reference numeral 74. As a result, the system 54 generates outputs or measurements 76 (Zk) at the given time step k that may be utilized by the filter 52 for estimation of residuals as will be described below.
As described above, the filter 52 receives the deterministic inputs 56 that are employed by the system 54 for generating the measurements 76. In the illustrated embodiment, the deterministic inputs 56 are applied to a function 78 (f({circumflex over (x)}k|k, uk)) for generating an a priori estimate of the state of the system 54 for the next time step from a previous time step k as represented by reference numerals 80 and 82. Again, a delay operator 84 (z−1) may be employed to incorporate any changes to the state due to any delays between subsequent time steps to generate an updated estimate 86 ({circumflex over (X)}k|k-1). Such updates may be incorporated into the existing state 80 as computed by summers 82 and 88. Further, a non-linear function 92 (h({circumflex over (x)}k|k-1)) may utilize the updated estimate 86 to estimate predicted measurements 94 ({circumflex over (z)}k).
In a presently contemplated configuration, the measurements 76 from the system 54 are collated and compared with the predicted outputs 94 as represented by reference numeral 96 to generate residuals (innovations) 98 (vk). Moreover, the generated residuals 98 are multiplied by an observer gain 100 (Kk({circumflex over (x)}k|k-1)) that may be employed to reduce the residuals 96 to customize the existing state 86. In this embodiment, the observer gain 100 includes a Kalman gain Kk that is given by the following equations:
Kk=Pk|k-1HkT(HkPk|k-1HkT+Rk) (1)
Pk|-1=AkPk-1|k-1AkT+GkQkGkT (2)
Pk|k=(I−KkHk)Pk|k-1 (3)
where:
Referring now to
In the illustrated embodiment, the residuals 118 are provided to a tracking filter 120 for generating “personalized” (i.e., system-specific) parameter estimates 122 for the engine, that is, estimates that are adapted to the particular aircraft and equipment rather than for a generic fleet of aircraft and equipment. In this embodiment, the tracking filter 120 includes an extended Kalman filter. The tracking filter 120 analyzes the residuals 118 and generates the personalized parameter estimates 122 for reducing the residuals 118 between the generated outputs and the sensed engine outputs 114.
The estimated personalized parameters 122 from the engine model 106 are applied to an engine model 124 for prediction of residuals. In the illustrated embodiment, the parameters of the engine model 124 are updated based upon the personalized parameters 122 to match the current state of the engine model 124 with the current state of the engine model 106. Moreover, the engine model 124 receives engine inputs 126 from a pre-determined number of prediction flights NP. Again, the engine inputs 126 may include engine control inputs and sensed inputs as described earlier. In this embodiment, the engine model 124 generates outputs based upon the engine inputs 126. Examples of generated outputs include temperature, pressure, rotor speed, efficiency, flow capacity, inter-component temperature and so forth.
Additionally, the engine model 124 receives engine sensed outputs 128 from the pre-determined number of prediction flights NP. Subsequently, the generated outputs from the engine model 124 are collated and compared with the engine sensed outputs 128 as represented by reference numeral 130. As a result, residuals 132 for the engine model 124 are estimated based upon the generated outputs and the engine sensed outputs 128. Advantageously, the residuals 132 from the personalized engine model 124 may be employed for diagnosing the faults in the engine via a multiple model hypothesis test 134.
In a presently contemplated configuration, the residuals 132 from the personalized engine model 124 are then compared with a set of faults or fault signatures 136 via the multiple model hypothesis test 134 for detecting and isolating the faults in the engine. In certain embodiments, fault probabilities 138 may be computed by the multiple model hypothesis test 134. In certain other embodiments, a severity estimate for the detected faults may be generated. In such embodiments, the severity estimate is calculated based upon the fault probabilities 138 and a magnitude of the fault signatures. Further, the estimated parameters 122 may be employed for generating a trend over time for detecting abnormal deterioration of the components of the engine.
As described above, the personalized parameter estimates 122 may be generated from the pre-determined number of training flights NT and the residuals 132 from the engine model 124 may be obtained from the pre-determined number of prediction flights NP.
In the illustrated embodiment, the personalized parameters from the training windows 142, 144 and 146 are utilized by the engine model for predicting engine outputs in the prediction windows as represented by reference numerals 154, 156 and 158. In addition, for each of the prediction windows 154, 156 and 158 the predicted engine outputs are compared with sensed outputs from the engine to generate residuals. As noted above, the residuals may be further utilized for detecting and isolating faults in the engine.
In the illustrated embodiment, the residuals 172 generated from the plurality of Kalman filters 164 are applied to a probability density function such as a Gaussian probability density function 174 for detecting faults based upon a likelihood of the residuals 172. In certain embodiments, Bayes rule 176 along with a hidden Markov model (HMM) 178 may be employed for determining fault probabilities 180 (P(fault′i′|z)) from the residuals 172. In the illustrated embodiment, a fault in the system 162 may be detected based upon the fault probabilities 180 and pre-determined thresholds.
In this embodiment, the fault probability 180 for each of the faults may be estimated between time updates by employing a probability transition matrix CP. The probability of ith fault at a given time k is estimated based upon the measurements up to time k−1 and is given by the following equation:
P(fi|tk,vk-1,vk-2, . . . )=CP·P(fi|tk-1,vk-1,vk-2, . . . ) (4)
where:
Further, a likelihood of the residual vk for each fault fi may be estimated by using a Gaussian distribution as given by the following equation:
where:
Next, at step 192 a sum of Gaussian PDF is calculated for a number of samples over a period of time to estimate a final exponent (Zfault) for the number of samples as represented by step 194. At step 196 Bayes rule may be applied to the residuals to determine fault probabilities for each of the fault as shown at step 198. In certain embodiments, the final vector of exponent (Zfault) may be augmented with an exponent for an unknown fault (Zunknown) to detect and isolate a fault other than the set of faults (step 200). In one embodiment, the unknown fault may function as a threshold for a detected fault that does not match any of the fault signatures. Thus, the augmentation of the final exponent with the unknown fault facilitates substantially accurate prediction of the fault probabilities. Further, the probability of occurrence of the unknown fault may be separated from the set of faults as represented by step 202.
As noted above, the present technique detects and isolates faults in the engine by analyzing the residuals between outputs estimated from an engine model and outputs measured from the engine. It should be noted that the estimated residuals may be analyzed in real-time on wing. Alternatively, the estimated residuals may be analyzed on a diagnostic location on ground. That is, parameter data, either raw or processed, may be transmitted from the aircraft to a ground location for computation of the derived parameters and residuals, and for analysis of the residuals as described above. This may be done in real-time, near real time, or even at a later time (e.g., following a flight).
In operation, the engine 206 receives input control inputs 220. Examples of such inputs include fuel flow, an active clearance control, variable geometry, power extraction, or combinations thereof for components of the engine 206. In addition, the control inputs 220 also include sensed inputs such as temperature, pressure, altitude, Mach number and combinations thereof. Further, the control inputs 220 are applied to the engine models 212 and 216 for the on-wing and the remote diagnostics units 208 and 210 for predicting outputs from the engine models 212 and 216 based upon the control inputs 220. The tracking filters 214 and 218 are configured to analyze the residuals between the predicted and sensed outputs from the engine 206 for generating personalized parameter estimates 222 for the particular engine 206.
In one embodiment, the engine model 212 of the on-wing diagnostic unit 208 may be employed for estimating outputs. Subsequently, the tracking filter 214 of the on-wing diagnostic unit may be employed for analyzing the residuals and for generating personalized parameter estimates 222. In another embodiment, the engine model 212 of the on-wing diagnostic unit 208 may be employed for estimating outputs and the estimated residuals between the estimated outputs and sensed outputs may be analyzed at the remote diagnostic location 210 via the tracking filter 218. Thus, a combination of the engine models 212 and 216 along with the tracking filters 214 and 218 may be employed for analyzing the residuals thereby facilitating the detection of faults in the engine 206.
The estimated personalized parameters 222 for the particular engine 206 are utilized for fault detection and isolation through a fault detection system 224. The fault detection system 224 analyzes the residuals and detects faults in the engine 206 by comparing the residuals with fault signatures via a fault detection module. Additionally, based upon the estimated personalized parameters 222 a trend of deterioration of the engine 206 may be generated via a trending module 228 to detect any abnormal deterioration of the components of the engine 206. In this embodiment, the parameters corresponding to the faults detected by the fault detection system 224 may be made available to a user through an output 230. Examples of such parameters include fault probabilities and severity estimates for the detected faults.
The various aspects of the technique described hereinabove have utility in turbomachinery for example, an aircraft engine, an industrial gas turbine and a steam turbine. As will be appreciated by those skilled in the art, the present technique provides an efficient fault detection system for an aircraft engine that personalizes an engine model to match the individual engine. In addition, the technique provides a mechanism to reduce the scatter in the engine outputs across a fleet of engines and thereby provides a relatively accurate estimate of an engine component health for the particular engine over a period of time.
While only certain features of the invention have been illustrated and described herein, many modifications and changes will occur to those skilled in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and changes as fall within the true spirit of the invention.