Next Article in Journal
Effect of Clay Mineral Composition on Low-Salinity Water Flooding
Previous Article in Journal
Experimental and Analytical Simulation Analyses on the Electrical Performance of Thermoelectric Generator Modules for Direct and Concentrated Quartz-Halogen Heat Harvesting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Gas Path Fault Contribution Matrix for Marine Gas Turbine Diagnosis Based on a Multiple Model Fault Detection and Isolation Approach

1
College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China
2
Centre for Efficiency and Performance Engineering, University of Huddersfield, Huddersfield HD1 3DH, UK
*
Authors to whom correspondence should be addressed.
Energies 2018, 11(12), 3316; https://doi.org/10.3390/en11123316
Submission received: 16 October 2018 / Revised: 9 November 2018 / Accepted: 13 November 2018 / Published: 27 November 2018
(This article belongs to the Section J: Thermal Management)

Abstract

:
To ensure reliable and efficient operation of gas turbines, multiple model (MM) approaches have been extensively studied for online fault detection and isolation (FDI). However, current MM-FDI approaches are difficult to directly apply to gas path FDI, which is one of the common faults in gas turbines and is understood to mainly be due to the high complexity and computation in updating hypothetical gas path faults for online applications. In this paper, a fault contribution matrix (FCM) based MM-FDI approach is proposed to implement gas path FDI over a wide operating range. As the FCM is realized via an additive term of the healthy model set, the hypothetical models for various gas path faults can be easily established and updated online. In addition, a gap metric analysis method for operating points selection is also proposed, which yields the healthy model set from the equal intervals linearized models to approximate the nonlinearity of the gas turbine over a wide range of operating conditions with specified accuracy and computational efficiency. Simulation case studies conducted on a two-shaft marine gas turbine demonstrated the proposed approach is capable of adaptively updating hypothetical model sets to accurately differentiate both single and multiple faults of various gas path faults.

1. Introduction

With increasing gas turbine complexity, requirements for ensuring reliability and security are growing significantly, and diagnosis of gas path component faults is becoming a major issue [1]. The main causes of gas path faults include fouling, erosion, corrosion, and foreign object damage (FOD), which degrade the performance of the gas turbine, further reduce its safety and stability, and lead to decreased fuel economy and increases in operation and maintenance costs [2]. Gas path fault diagnosis approaches aim to detect and locate a fault in time to prevent a catastrophic consequence, and therefore contribute to the development of a maintenance schedule. A large number of gas path fault diagnosis approaches have been proposed and achieved good results in practical applications, such as neural networks [3,4,5,6], gas path analysis derivatives [7], genetic algorithms [8], and adaptive estimation [9,10,11]. Fault detection and isolation (FDI), as the core of the gas turbine fault diagnosis, has received extensive attention from researchers [12,13,14,15].
The multiple model (MM) approach is highly regarded due to its ability to simultaneously ensure high levels of diagnosis accuracy whilst maintain acceptable computational costs. It is based on hypothesis testing and conditional probability, which transforms a complex problem into several simple problems to detect and isolate the fault [16]. The MM approach was pioneered by Magill [17]. Up to now, three generations of the MM approach have been designed, including autonomous MM, interacting MM, and variable structure MM. A detailed description of these approaches can be found in References [16,18]. Maybeck applied the MM approach to a gas turbine to detect and isolate sensor and actuator faults [19,20,21,22], and Khorasani applied the MM approach to gas path fault detection and isolation for jet engines [23,24,25]. In Reference [23], a nonlinear MM-FDI approach based on the extended Kalman filter (EKF) and the unscented Kalman filter (UKF) was proposed for a jet engine. Then, the MM-FDI approach was used to detect and isolate the gas path fault, and a hierarchical architecture was developed that enabled the detection and isolation of both single and multiple faults in the jet engine [24].
Gas path fault, actuator fault, and sensor fault are three types fault of a gas turbine. Application of the MM-FDI approach in gas path fault identification is still under development. One of the main problems is the design of hypothetical fault models [26], which have a significant impact on the performance of the MM-FDI approach. The gas turbine healthy model set can be described by a set of linear models in Equation (1):
x k + 1 = A i x k + B i u k + w k y k = C i x k + v k
where i is an operating point, i = 1, 2, …, M; x n × 1 , y l × 1 , and u t × 1 denote the state vector, measurement vector, and control vector, respectively; and n, l, and t denote the dimension of the state vector, measurement vector, and control vector, respectively. w k n × 1 ( 0 , Q ) and v k l × 1 ( 0 , R ) are independent discrete-time process noises and measurement noises, respectively; and A i n × n , B i n × t , and C i l × n denote the state matrix, control matrix, and output matrix of the healthy condition at the operating point i, respectively.
Faults in the sensors cause changes in the measurement vector y, and faults in the actuators cause changes in the control vector u. Therefore, hypothetical models for sensor and actuator faults can be established by directly changing the corresponding measurement vector and control vector of the healthy model in Equation (1) [18,22,24]. However, the consequences of gas path faults due to component defects cause changes in the system matrices A, B, and C, hence the traditional gas path fault model can be described in the form of Equation (2):
x k + 1 = A i x k + B i u k + w k y k = C i x k + v k
where A i n × n , B i n × t , and C i l × n are the system matrices of a gas path fault condition at the operating point i .
According to Equation (2), each gas path fault at the operating point i is matched to a hypothetical fault model, which means that the system matrices corresponding to each gas path fault need to be established separately in advance for applying the MM-FDI approach. Therefore, establishing the hypothetical fault models is a very time-consuming task, and they are also not easy to update automatically online. Although the proposed hierarchical architecture can reduce the number of models of each level [24], due to the stochasticity and diversity of fault occurrences, the establishment of the most possible hypothetical fault models will face combinatorial explosion problems, which make the MM-FDI approach difficult to implement practically.
In this paper, a gas path fault contribution matrix (FCM) was introduced to improve the performance of the MM-FDI approach to address the above problem, and a systemic method to select the operating points that form the healthy model set was carried out based on gap metric analysis to achieve MM-FDI over a wide range of operating conditions. To evaluate the effectiveness of the FCM based MM-FDI approach, several simulation case studies on a two-shaft marine gas turbine were conducted. Performance of the proposed approach for FDI and the influence of the number of available measurements and measurement outliers was analyzed for single and multiple fault cases.
The remainder of this paper is organized as follows: Details pertaining to the FCM based MM-FDI approach are introduced in Section 2. In Section 3, a gap metric analysis method is introduced to yield the healthy model set. Section 4 shows the application of the proposed approach to a two-shaft marine gas turbine, and the results of the simulation case studies are analyzed and discussed. Finally, the main conclusions are presented in Section 5.

2. The FCM Based MM-FDI Approach

2.1. Model Set Design Based on FCM

Inspired by the estimation algorithm used in adaptive performance optimization of turbofan engines [27], this paper introduces the gas path fault contribution matrices E and D into Equation (1) as additive terms of the healthy model, so that they can be adjusted easily to balance the deviations caused by gas path faults. This gas turbine model set design based FCM can be represented in Equation (3):
{ x k + 1 = A x k + [ B E ] [ u k f k ] y k = C x k + [ 0 D ] [ u k f k ]
where the matrices E n × q and D l × q denote FCMs which represent the effect of gas a path fault on the healthy model, q represents the number of gas path faults in the gas turbine, f k = i = 1 q b i z i is the gas path fault vector, b i is the fault amplitude, and z i is the location of i - th fault. In particular, f k is a q × 1 zero-valued vector for the healthy condition, and it becomes a nonzero value when a gas path fault has occurred.
In Equation (3), the j-th hypothetical model at the operating point i can be represented by
{ x k + 1 j = A i x k j + [ B i E i ] [ u k j f k j ] + w k y k j = C i x k j + [ 0 D i ] [ u k j f k j ] + v k
As shown in Equation (4), different gas path hypothetical fault models can be established by directly quantifying the value of the corresponding gas path fault matrix f k , which is similar to the effect of the control vector u k on the actuator fault model.

2.2. Model Conditional Filtering

The Kalman filter algorithm for the j - th hypothetical model at the i - th operating point is shown in Equation (5).
T i m e u p d a t e : x ^ k + 1 | k ( j ) = A i x ^ k | k ( j ) + [ B i E i ] [ u k ( j ) f k ( j ) ] T P ^ k + 1 | k ( j ) = A i P k | k ( j ) A i T + Q M e a s u r m e n t u p d a t e : K k ( j ) = P k + 1 | k ( j ) C i T [ C i P k + 1 | k ( j ) C i T + R ] 1 x ^ k + 1 | k + 1 ( j ) = x ^ k + 1 | k ( j ) + K k ( j ) [ y k ( C i x ^ k + 1 | k ( j ) + [ 0 D i ] [ u k ( j ) f k ( j ) ] T ) ] P k + 1 | k + 1 ( j ) = P ^ k + 1 | k ( j ) K k ( j ) C i P ^ k + 1 | k ( j )
where the Kalman filter algorithm contains two steps, time update and measurement update; therein, x ^ k | k and x ^ k + 1 | k + 1 are the system state estimate at time k and k + 1 , respectively; x ^ k + 1 | k is the system state predict at time k + 1 ; P k | k and P k + 1 | k + 1 are the state estimate error covariance at time k and k + 1 , respectively; P ^ k + 1 | k is the state predict error covariance at time k + 1 ; and K k is the Kalman gain at time k .

2.3. Model Probability Update

The conditional probability μ i ( k ) is used to represent the approximation of each hypothetical model to the gas turbine actual condition. The gas turbine actual conditions A C match H m i when a given measurement vector y k is at a discrete time k ; that is
μ i ( k ) = Pr [ A C = H m i | y ( t k ) = y k ]
The conditional probability at a given time k of the hypothetical model set can be recursively calculated by the Bayesian Law according to the conditional probability value of the previous moment and the Gaussian probability density corresponding to the current filter residual, as shown in Equation (7).
μ i ( k ) = f ( y k | H m i , y k 1 ) μ i ( k 1 ) j = 1 W f ( y k | H m j , y k 1 ) μ j ( k 1 )
where f ( y k | H m i , y k 1 ) denotes the conditional probability density of the i - th hypothetical model when the measurements are y k , and μ i ( k 1 ) denotes the conditional probability of the i - th hypothetical model at the time k 1 . The Gaussian conditional probability density function is shown in Equation (8).
f ( y k | H m j , y k 1 ) = 1 ( 2 π ) l / 2 | S k j | 1 / 2 exp [ 1 2 ( γ k j ) T S k j γ k j ]
where l is the dimension of the measurement parameters, and S k j and γ k j are the innovation covariance and the innovation of the filter corresponding to the j - th hypothetical model, respectively, which can be obtained in the filtering process.

2.4. Fault Detection and Isolation

The conditional probability of the hypothetical model set can be used as an indication of a fault, as it provides a meaningful measure of how likely each fault is at a given time [28]. Therefore, the actual fault can be detected and isolated by using the maximum probability criteria, as shown in Equation (9):
j   =   arg   max i = 1 W μ i ( k ) > μ T H
where j represents that the j - th model is closest to the actual condition, and μ T H is a preset detect and isolate threshold. To ensure a fast response to changes due to fault effects, a minimum μ i = 0.001 was used for each model [27]. Additionally, the preset detection and isolation threshold μ T was set to 0.98 to reduce the amount of warnings [22].

2.5. Hypothetical Model Update

In practice, multiple faults exist in which a second fault occurs after the first fault. In this scenario, no hypothetical model in the initial hypothetical model set matches the actual fault condition. In addition, a gas path abrupt fault with large amplitude may also result in incorrect detection. Therefore, to maintain the performance of the proposed approach, it is a key issue to update the hypothetical model set to detect multiple faults after detecting the first fault, or to detect faults with larger amplitude.
Assuming the initial hypothetical model set is constant and has W hypothetical models, the demonstration of the hypothetical model set update process is shown in Figure 1; the initial hypothetical model set is shown in the left of the figure, where the subscript j 1 denotes the ( j 1 ) - th element of f k . The initial fault amplitude is assumed to have a deviation of −1% (−3% is considered to be severe in practice [29,30]) for each hypothetical fault model. If the hypothetical model j has the largest conditional probability, and it is in excess of the preset threshold μ T H , then the first fault is detected and isolated as the hypothetical model j . Subsequently, the initial hypothetical model set will be updated to detect whether or not the current condition is multiple faults or a single fault that is more serious. In the hypothetical model set update process, each hypothetical model should contain the first detected fault; the j - th hypothetical model represents a more serious fault (−2%), as shown in the right of Figure 1. In the case of more faults, the update process of the hypothetical model set is in the same manner.
The detailed hypothetical model set update algorithm is shown in Table 1. First, the detected hypothetical model is determined to obtain the current fault type and the number of faults, and determine the matrix F u p d a t a corresponding to the updated hypothetical model set. Then, each column of F u p d a t a is substituted into Equation (4) to obtain the updated hypothetical model set. After the hypothetical model set is updated, the conditional probability of each model and the filter will be reinitialized.

2.6. Implementation of FCM Based MM-FDI

The overall flow diagram of FCM based MM-FDI is shown in Figure 2, which consists of following steps: Establishing a gas turbine model set based on FCM that can represent system structure at the healthy condition, and where all parameters use the percentage change in the model set. Then, based on a hypothetical model set update algorithm, a bank of Kalman filters corresponding to the hypothetical model set which represents hypothetical conditions of the gas turbine is designed. The conditional probability of each hypothetical model is recursively calculated according to the filter residuals and Bayesian law. Finally, the model that most closely aligns with the actual condition of the gas turbine can be detected and isolated using the maximum probability criteria. Compared with the MM-FDI algorithm in References [23,24], the FCM based MM-FDI approach proposed in this paper solves the problem of gas path fault model set establishment and automatic update of the hypothetical model set.

3. Selecting Operating Points Based on Gap Metric Analysis

Selecting the few operating points to form the healthy model set over a wide operating condition is another issue of the MM approach. It usually uses equal intervals or control command [24], which may lead to the redundant or incomplete selection of the linear model. To enhance the computational efficiency, the trajectory piecewise-linear method [31] and range reduction techniques [32] have been proposed for the optimal selecting of operating points.
In this paper, a systemic approach for selecting the operating points of the gas turbine linearized models was carried out based on gap metric analysis. EI-Sakkary [33] pioneered a gap metric that can be used as an indication of the approximation of two linear systems. For two linear models L 1 and L 2 with the same input and output, the gap between two linear models is [34]
δ ( L 1 , L 2 ) = max ( δ ( L 1 , L 2 ) , δ ( L 2 , L 1 ) )
where δ ( L i , L j ) is
δ ( L i , L j ) = inf Q H ( M i N i ) ( M j N j ) Q
where i , j = { 1 , 2 } , i j . Q H , and M i and N i , are normalized right coprime factorizations of L i , and defined as
L i = N i M i 1 , i = 1 , 2
According to the gap metric, the range of the gap metric for two linear models is [0, 1], and the smaller the gap metric, the closer the dynamic response between the two linear models. Therefore, when the gap metric of the two linear models is less than the preset threshold δ T H , it means that just one of the models is sufficient, as the other one has become redundant. Otherwise, the two models are not redundant and need to be selected simultaneously.
In this paper, a set of operating points with equal intervals and number of z are selected in the range of the control variable u , and the interval φ between any two adjacent operating points u p and u p + 1 ( p = 1 , 2 , , z 1 ) satisfies Equation (13):
φ = | u p u p + 1 | δ u
where δ u is a sufficiently small value such that the gap metric of the linearization models corresponding to any two adjacent points is less than the preset threshold δ T H .
Correspondingly, a set of linearized models L i ( i = 1 , , z ) corresponding to z operating points was established. Moreover, the gap metric between any two linearized models can be calculated. Finally, the gap metric matrix G of all linearized models is obtained in Equation (14):
G = [ δ ( L 1 , L z ) δ ( L 2 , L z ) δ ( L z , L z ) δ ( L 1 , L 2 ) δ ( L 2 , L 2 ) δ ( L z , L 2 ) δ ( L 1 , L 1 ) δ ( L 2 , L 1 ) δ ( L z , L 1 ) ]
where for the matrix G z × z , δ ( L i , L j ) = δ ( L j , L i ) , i , j = 1 , 2 , , z , and δ ( L i , L j ) = 0 , when i = j . i = 1 defines the first operating point.
According to the preset threshold δ T and the obtained gap metric matrix G, a set of operating points, from which corresponding linearized models form the healthy model set, can be selected from the z operating points. It can be seen from Equation (14) that the i - th column in matrix G is the gap metric between the i - th linearized model and z linearized models. We keep i unchanged and j increases sequentially, then compare δ ( L i , L j ) with the preset threshold δ T H . If δ ( L i , L j ) δ T H , then the linearized model at an operating point j can be replaced by the linearized model at the operating point i . Repeating the above process until δ ( L i , L j ) > δ T H or L j = L z , then the j 1 i linearized models from the operating point i to j 1 can be represented by the linearized model at the operating point i . Meanwhile, the linearized model at the operating point j is the next selected operating point.

4. Simulation Results and Discussion

In this section, the proposed FCM based MM-FDI approach is applied to a two-shaft marine gas turbine. Several case studies were conducted to evaluate the performance of the proposed approach, which include single fault detection and isolation, multiple fault detection and isolation, and the influence of the number of available measurements and the measurements’ outliers.
The general layout of the two-shaft marine gas turbine is shown in Figure 3 [35]. The gas path components include a compressor, a combustion chamber, a compressor turbine, and a power turbine, wherein the power turbine drives the propeller through the reduction gearbox, which is simplified as a load, as shown in Figure 3.
The nonlinear dynamic model of the two-shaft marine gas turbine was taken as a real engine, which was established according to the available literature [36], as shown in Equation (15); for detailed information refer to authors’ previous work [35].
N ˙ 1 = ( 30 π ) 2 1 J 1 N 1 [ m C T η m c p g ( T 3 T 4 ) m C c p ( T 2 T 1 ) ] N ˙ 2 = ( 30 π ) 2 1 J 2 N 2 [ m P T c p g ( T 4 T 5 ) Φ N 2 3 ] T ˙ 3 = R g T 3 P 3 V 1 c v g [ k ( c p g T 2 m C + L H V η C C w f c p g T 3 m C T ) c p g T 3 ( m C + w f m C T ) ] P ˙ 3 = P 3 T 3 T ˙ 3 + R g T 3 ( m C + w f m C T ) V 1 P ˙ 4 = ( m C T m P T ) R g T 4 V 2 P ˙ 5 = ( m P T Γ P 5 ) R g T 5 V 3
where J 1 and J 2 are the inertia of the compressor shaft and power turbine shaft, respectively; Φ is the relationship coefficient between N 1 and the load; L H V is the fuel low heating value; and V 1 , V 2 , and V 3 are the component volumes. m C , m C T , and m P T represent the mass inside the compressor, compressor turbine, and power turbine, respectively. c p and c p g represent the heat capacity of air and gas at constant pressure, respectively; c v g represents the volumetric heat capacity of gas; and R g denotes the gas constant. η m denotes the mechanical efficiency, and Γ is the relationship coefficient between m P T and P 5 .
Next, the performance of the FCM based FDI approach was tested in a simulation environment, as shown in Figure 4. According to Equation (15), there are six component performance parameters in a two-shaft marine gas turbine: compressor mass flow m C , compressor efficiency η C , compressor turbine mass flow m C T , compressor turbine efficiency η C T , power turbine mass flow m P T , and power turbine efficiency η P T . When the gas turbine is in a healthy condition, the compressor mass flow and efficiency are denoted m C , H and η C , H , respectively; compressor turbine mass flow and efficiency are denoted m C T , H and η C T , H , respectively; and power turbine mass flow and efficiency are denoted m P T , H and η P T , H , respectively.
In this paper, gas path faults were simulated via changing component performance parameters. Typical single gas path faults used in this paper are shown in Table 2. When two or more single faults occurred simultaneously they were considered multiple faults. The sensor noise used in this paper was Gaussian noise, where the mean value was zero, and the relative standard deviation of each measurement variable is shown in Table 3.

4.1. FCM Based Model Set Testing

In the two-shaft marine gas turbine model set, as shown in Equation (3), the state vector was x = [ N 1 , N 2 , T 3 , P 3 , P 4 ] T , the measurement vector y = [ N 1 , N 2 , T 2 , P 2 , T 4 , P 4 , T 5 ] T , and the control vector u = w f . Sample time was 0.02 s.

4.1.1. Mode Set Accuracy Testing

When f k = 0 , the FCM based model set is a set of healthy linearized models. In this paper, the fitting approach based on perturbation [37] was used to establish the gas turbine linearized model, as shown in Equation (4), and gap metric analysis was used to select the operating point.
The range of the control vector w f was [0.3, 1]. With an equal interval of 0.01 for setting the operating points, the operating range was divided into 70 subranges, and a total of 71 gas turbine linearization models were established. According to Equation (14), the matrix of the gap metrics between the gas turbine linearization models was calculated and is shown in Figure 5.
The gap metric matrix is a symmetric matrix, and the diagonal elements are the gap metrics of each model with itself, hence the value is zero. The gas turbine operating range was segmented by comparing the pre-set threshold δ T H = 0.25 starting at the design operating point w f = 1 (namely L 1 ), and six operating points were selected, as shown in Table 4, to yield the FCM based model set.
The dynamic response of each measurement parameter was compared between the gas turbine nonlinear model [35] and FCM based model set. Figure 6 presents the corresponding differences between two models at a given fuel flow schedule. In this case, the fuel flow w f decreases gradually from 1.0 to 0.3 in the period from 0 to 40 s. It can be seen that errors at each selected operating point appear to be the largest for each subrange, which agrees with the operating point selection using the gap metric method. However, the absolute value of the maximum error for each measured parameter is less than 4%; that is, it is in an acceptable error range. At the expense of certain model simulation accuracy, the model set covering the wide operating condition only uses 8.5% of the operating points of gas turbine linearization models, which greatly reduces the storage space and calculation requirements of the model set.

4.1.2. Hypothetical Fault Simulation

To verify whether the FCM based model set can be used to establish the hypothetical models, a comparative study was conducted with a nonlinear model in Equation (15), the traditional gas path fault model shown in Equation (2), and the hypothetical fault model auto-generated from the FCM based mode set, as shown in Equation (4). In this section, the traditional gas path model F m c at the operating point w f = 1.0 is selected. Figure 7 shows the comparison dynamic tracking response in artificial simulation scenarios. The gas turbine was abruptly changed from operating point w f = 1.0 to operating point w f = 0.99 at t = 15 s; a F m c decrease of 1% occurs at t = 15 s and recovery to the healthy condition at t = 25 s; after 10 s, a F m c t decrease of 1% occurs until the end.
It can be seen from Figure 7 that the hypothetical fault model auto-generated from the FCM based mode set has the same fault dynamic tracking response with a nonlinear model under a steady state condition and transient condition. The traditional gas fault model can only simulate the specified fault (here is F m c ). The hypothetical fault model auto-generated from FCM based mode set has adaptive fault tracking capability through adjusting f k , which makes computation simpler than the traditional gas path fault model.

4.2. Single Fault Scenarios

4.2.1. FDI Results of a Single Fault under Different Operating Conditions

In this section, the gas turbine was worked under the steady condition at the operating point w f = 1.0 until t = 360 s, and then the transient condition while the operating point dropped gradually to w f = 0.3 at t = 720 s. Six gas path faults occurred separately in sequence under the steady condition and transient condition with an amplitude of −1%, and the diagnosis results are shown in Figure 8.
It can be seen that the FCM based MM-FDI approach can accurately detect and isolate single gas path faults under both the steady and transient conditions at expected detection times and isolation periods, respectively.

4.2.2. Performance under Different Single Fault Amplitudes

In this section, 100 simulations of single gas path faults were performed each time at the operating point, and repeated 10 times for a total of 1000 simulations based on the gas turbine nonlinear model. For each simulation, the type of fault was randomly generated, and the range of the fault amplitude bi was [−2.7%, −0.7%]. Keeping the hypothetical model set in an initial model sequence (fault amplitude bi is −1%), the correct detections (CD), the incorrect detections (ID), and the missed detections (MD) of the FCM based MM-FDI approach under a single fault were obtained from this Monte Carlo test, and are listed in Table 5.
It can be seen that the proposed approach allows 95.5% faults to be detected and isolated under the initial hypothetical mode set, which demonstrates that the proposed approach produces high diagnostic accuracy under different single fault scenarios.
All the incorrect detections are on Fηpt, which incorrectly detected as Fmc, as shown in Table 5. For an in-depth understanding, all the fault amplitudes of Fηpt faults that were simulated in the case study and the corresponding FDI results are presented in Figure 9. It can be seen that all incorrect detection occurs when the fault amplitude is too large. Therefore, the cause of the incorrect detection is that the amplitude of the initial hypothetical Fηpt model is small, and it does not match the actual fault condition. Fortunately, the large amplitude fault can be accurately detected and isolated when the hypothetical model set is updated, as advocated in this method.

4.2.3. Performance under Different Numbers of Available Sensors

In practice, sensors often fail to provide effective measurement information, which results in the number of available sensors decreasing. Therefore, the influence of the number of available sensors on the performance of the proposed approach was studied at the operating point w f = 1.0 . Three cases were considered, seven sensors ( N 1 , N 2 , T 2 , P 2 , T 4 , P 4 , and T 5 ), four sensors ( N 1 , N 2 , P 2 , and T 4 ), and two sensors ( N 1 and N 2 ), and each fault had a 1% decrease in its output. The detection time and isolation time of each fault are shown in Table 6.
As shown in Table 6, the proposed approach can detect and isolate the fault with four sensors. However, several faults were not able to be differentiated, which are denoted as either MD or ID in the scenario of two available sensors. In addition, the fault detection time and the isolation time increases as the number of sensors decreases, which is mainly due to the reduction in the available information to distinguish different faults that the sensors can provide.

4.2.4. Performance under Measurement Outliers

It is inevitable to have outliers in measurements due to various interferences, such as failed communication, hardware malfunction, or environment change, which can result in a false alarm. Therefore, some gas path diagnosis methods have been proposed for pre-processing the measurement signals to remove these outliers and noise before FDI. In this case, the outliers were artificially added to the measurement signals and the outliers were selected at the 2% level of the nominal value of these measurements. Suppose a F m c decrease of 1% occurs at t = 20 s, the variation of each measurement parameter and the probability of each hypothetical model are shown in Figure 10. It can be seen that the proposed approach can accurately detect and isolate faults when there are outliers, and the outliers only slightly disturb the conditional probability of the model. Therefore, the FCM-based FDI approach proposed in this paper is robust in the measurement of outliers.

4.3. Multiple Fault Scenarios

4.3.1. FDI Results of Multiple Faults in the Gas Path

Figure 11 shows the FDI results for the scenarios of a −2% change in F η p t at t = 25 s and a 1% change in F η c at t = 35 s under operating point w f = 1.0 . In this case, the hypothetical model set update approach was used for detecting and isolating these faults. It can be seen that when a F η p t decreasing by 2% fault occurs, the conditional probability of the 7# hypothetical model in Figure 11a is the largest and exceeds μ T H , which means this hypothetical model is closest to the actual fault condition. Then, the hypothetical model set is updated and the conditional probability is reinitialized correspondingly. With the updated hypothetical model set, the 7# hypothetical model, as shown in Figure 11b, is closest to the actual fault condition, so the conditional probability of the model is largest and over μ T H . Then, the hypothetical model set will be updated again and the conditional probability will be reinitialized again. Finally, as shown in Figure 11c, the F η p t and the F η c faults are accurately detected and isolated.

4.3.2. Performance under Multiple Faults

A simulation was also conducted for the scenario of multiple faults at the operating point w f = 1.0 , and the detection and isolation results are shown in Figure 12 and Table 7. In this scenario, the type of fault was random and the range of the fault amplitude b i was [−2.7%, −0.7%]. In addition, the fault amplitude of the gas path fault caused by fouling, corrosion, and erosion may be gradual, so both abrupt and gradual changes in fault amplitude were considered. In the case of an abrupt fault, the first fault occurs at t = 20 s and the second fault occurs at t = 70 s. In the case of a gradual fault, the fault amplitude of both faults changes from t = 0 s, with the first fault changing to a given amplitude at t = 20 s, and the second fault changing to a given amplitude at t = 70 s.
It can be seen that, regardless of whether the fault amplitude changes abruptly or gradually, the first fault was 100% accurately detected and isolated using the proposed approach. In the case of an abrupt fault, the second fault shows five faults were incorrectly identified as three or more faults, and one fault was missed due to its very small amplitude. In the case of a gradual fault, the second fault shows that five faults were incorrectly identified to be three faults. Therefore, the proposed approach can effectively improve the FDI performance of large amplitude single faults, and has high FDI performance with 94% and 95% accuracy in the multiple abrupt faults and gradual faults scenarios, respectively.

4.3.3. FDI Results of Multiple Faults in both the Actuator and Gas Path

In Equation (3), gas path faults are converted into a type of fault in the model set that is similar to the actuator fault. Due to the good scalability of the FCM based FDI approach, it is easy to apply the proposed approach to gas path and actuator multiple faults.
In this section, the actuator fault corresponds to the fuel flow ( F w f ). The fuel actuator gets stuck at some value [38] a 0 = ( w f / w f , H 1 ) × 100 % , which is the fuel actuator fault F w f ( a 0 % ) ; u k = m k + a 0 , where m k is controller fuel mass output. Considering the healthy condition and six gas path faults, as shown in Table 2, there will be eight hypothetical models in the model set. Although the number of hypothetical models has increased, the process of the FDI of multiple faults in both the gas path and actuator is consistent with the process of the FDI of gas path faults.
Figure 13 shows the FDI results in the scenarios of the actuator fault, where the Fwf decreases by 5% fault occurs at t = 15 s, the F m c decreases by 1% fault occurs at t = 35 s, and the F η p t decreases by 1% fault occurs at t = 45 s. In these scenarios, the gas turbine engine works at the operating point w f = 1.0 . Assume that the initial fault amplitude of Fwf is −5% and the gas path fault amplitude is −1%. The FDI results show that both the gas path fault and actuator fault are accurately detected and isolated, which indicates that the occurrence of an actuator fault will not affect the detection and isolation of a gas path fault.

5. Conclusions

In this paper, an FCM based MM-FDI approach was developed to detect and isolate the gas path faults of gas turbines over a wide operating condition. It shows that the FCM can be realized via an additive term of the healthy model set, allowing hypothetical models for various gas path faults to be easily established and updated online. A gap metric analysis for operating points selection converts the nonlinearity of the gas turbine into multiple linearized models over a wide range of operating conditions with specified accuracy and computational efficiency. Several simulation case studies on a two-shaft marine gas turbine were conducted, and the performance of the proposed approach for the FDI of a single fault and multiple faults, and the influence of measurement outliers, was analyzed. The simulation results show that the proposed FCM based MM-FDI approach realizes the adaptive generation of the fault model set by taking the gas path fault as an additive term and extending the gas path fault matrix into the control vector, which makes MM-FDI approaches achievable and efficient for online gas path fault detection and isolation application. The gap metric analysis for selecting operating points of the healthy model set greatly reduces the storage space and calculation requirement of the model set.

Author Contributions

This research was completed through a collaboration of all authors. Q.Y. and Y.C. proposed the FCM based MM-FDI approach and wrote the original manuscript; S.L. was the team leader of this work and responsible for coordination; F.G., A.S., and Y.C. designed the case studies and analyzed the results.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities [HEUCFP201722].

Acknowledgments

Gratitude is extended to Harbin Engineering University for supporting Y.C. to carry out collaborative research in the Centre for Efficiency and Performance Engineering at the University of Huddersfield.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Symbols
A i state matrix of the healthy condition at the operating point i
A i state matrix of a gas path fault condition at the operating point i
B i control matrix of the healthy condition at the operating point i
B i control matrix of a gas path fault condition at the operating point i
C i output matrix of the healthy condition at the operating point i
C i output matrix of a gas path fault condition at the operating point i
D gas path fault contribution matrix for the output vector
E gas path fault contribution matrix for the state vector
G the gap metric matrix of all linearized models
H m i the i-th hypothetical model
J inertia of shaft
K Kalman gain
L i the i-th linear model
M i , N i the normalized right coprime factorizations of linear model Li
N rotational speed
P pressure
Q process noise covariance
R measurement noise covariance for FDI
R g gas constant
S innovation covariance
T temperature
V component volume
W the number of hypothetical models
b fault amplitude
c p air constant pressure specific heat
c p g gas constant pressure specific heat
c v g gas constant volume specific heat
f gas path fault vector
k adiabatic index
l the dimension of the measurement vector
m gas/air mass flow of component
n the dimension of the state vector
q number of gas path faults
s fault amplitude changes during model set updating
t the dimension of the control vector
u control vector
v k measurement noise
w k process noise
w f fuel mass flow
x state vector
y measurements vector
z fault location vector
Greek
Γcoefficient between mass flow and pressure
Φ coefficient between N2 and load
γ filter innovation
δ ( L i , L j ) gap metric between linear models Li and Lj
δ T H the preset threshold of the gap metric
η isentropic efficiency
μ T H the preset threshold of the conditional probability
μ i the conditional probability of the i-th hypothetical model
φ the interval between any two adjacent operating points
Subscript
C compressor
C C combustion chamber
C T compressor turbine
H health condition
P T power turbine
k discrete time k

References

  1. Volponi, A.J. Gas Turbine Engine Health Management: Past, Present, and Future Trends. J. Eng. Gas Turbines Power 2014, 136. [Google Scholar] [CrossRef]
  2. Zwebek, A.; Pilidis, P. Degradation Effects on Combined Cycle Power Plant Performance—Part I: Gas Turbine Cycle Component Degradation Effects. J. Eng. Gas Turbines Power 2003, 125, 651. [Google Scholar] [CrossRef]
  3. Lu, P.-J.; Zhang, M.-C.; Hsu, T.-C.; Zhang, J. An Evaluation of Engine Faults Diagnostics Using Artificial Neural Networks. J. Eng. Gas Turbines Power 2001, 123, 340–346. [Google Scholar] [CrossRef]
  4. Ogaji, S.; Sampath, S.; Singh, R.; Probert, D. Novel approach for improving power-plant availability using advanced engine diagnostics. Appl. Energy 2002, 72, 389–407. [Google Scholar] [CrossRef]
  5. Ogaji, S.O.T.; Singh, R. Gas path fault diagnosis framework for a three-shaft gas turbine. Proc. Inst. Mech. Eng. Part A J. Power Energy 2003, 217, 149–157. [Google Scholar] [CrossRef]
  6. Campora, U.; Cravero, C.; Zaccone, R. Marine gas turbine monitoring and diagnostics by simulation and pattern recognition. Int. J. Naval Archit. Ocean Eng. 2018, 10, 617–628. [Google Scholar] [CrossRef]
  7. Li, Y.G.; Korakiantis, T. Nonlinear Weighted-Least-Squares Estimation Approach for Gas-Turbine Diagnostic Applications. J. Propuls. Power 2011, 27, 337–345. [Google Scholar] [CrossRef]
  8. Zedda, M.; Singh, R. Gas Turbine Engine and Sensor Fault Diagnosis Using Optimization Techniques. J. Propuls. Power 2002, 18, 1019–1025. [Google Scholar] [CrossRef]
  9. Pu, X.; Liu, S.; Jiang, H.; Yu, D. Adaptive gas path diagnostics using strong tracking filter. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2014, 228, 577–585. [Google Scholar] [CrossRef]
  10. Lu, F.; Huang, J.; Lv, Y. Gas Path Health Monitoring for a Turbofan Engine Based on a Nonlinear Filtering Approach. Energies 2013, 6, 492–513. [Google Scholar] [CrossRef] [Green Version]
  11. Chang, X.; Huang, J.; Lu, F.; Sun, H. Gas-Path Health Estimation for an Aircraft Engine Based on a Sliding Mode Observer. Energies 2016, 9, 598. [Google Scholar] [CrossRef]
  12. Gupta, S.; Ray, A.; Sarkar, S.; Yasar, M. Fault detection and isolation in aircraft gas turbine engines. Part 1: Underlying concept. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2008, 222, 307–318. [Google Scholar] [CrossRef] [Green Version]
  13. Abbasi Nozari, H.; Aliyari Shoorehdeli, M.; Simani, S.; Dehghan Banadaki, H. Model-based robust fault detection and isolation of an industrial gas turbine prototype using soft computing techniques. Neurocomputing 2012, 91, 29–47. [Google Scholar] [CrossRef]
  14. Sadough Vanini, Z.N.; Khorasani, K.; Meskin, N. Fault detection and isolation of a dual spool gas turbine engine using dynamic neural networks and multiple model approach. Inf. Sci. 2014, 259, 234–251. [Google Scholar] [CrossRef]
  15. Sina Tayarani-Bathaie, S.; Khorasani, K. Fault detection and isolation of gas turbine engines using a bank of neural networks. J. Process Control 2015, 36, 22–41. [Google Scholar] [CrossRef]
  16. Li, X.R.; Jilkov, V.P. Survey of maneuvering target tracking. Part V: Multiple-model methods. IEEE Trans. Aerosp. Electron. Syst. 2005, 41, 1255–1321. [Google Scholar] [CrossRef]
  17. Magill, D.T. Optimal adaptive estimation of sampled stochastic processes. IEEE Trans. Autom. Control 1965, 10, 434–439. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, Y.; Li, X.R. Detection and diagnosis of sensor and actuator failures using IMM estimator. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 1293–1313. [Google Scholar] [CrossRef]
  19. Menke, T.E.; Maybeck, P.S. Sensor/actuator failure detection in the Vista F-16 by multiple model adaptive estimation. IEEE Trans. Aerosp. Electron. Syst. 1995, 31, 1218–1229. [Google Scholar] [CrossRef] [Green Version]
  20. Eide, P.; Maybeck, P. An MMAE failure detection system for the F-16. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 1125–1136. [Google Scholar] [CrossRef] [Green Version]
  21. Hanlon, P.D.; Maybeck, P.S. Multiple-model adaptive estimation using a residual correlation Kalman filter bank. IEEE Trans. Aerosp. Electron. Syst. 2000, 36, 393–406. [Google Scholar] [CrossRef]
  22. Maybeck, P.S. Multiple model adaptive algorithms for detecting and compensating sensor and actuator/surface failures in aircraft flight control systems. Int. J. Robust Nonlinear Control 1999, 9, 1051–1070. [Google Scholar] [CrossRef]
  23. Naderi, E.; Meskin, N.; Khorasani, K. Nonlinear Fault Diagnosis of Jet Engines by Using a Multiple Model-Based Approach. J. Eng. Gas Turbines Power 2011, 134, 011602. [Google Scholar] [CrossRef]
  24. Meskin, N.; Naderi, E.; Khorasani, K. A Multiple Model-based Approach for Fault Diagnosis of Jet Engines. IEEE Trans. Control Syst. Technol. 2013, 21, 254–262. [Google Scholar] [CrossRef]
  25. Pourbabaee, B.; Meskin, N.; Khorasani, K. Sensor Fault Detection, Isolation, and Identification Using Multiple-Model-Based Hybrid Kalman Filter for Gas Turbine Engines. IEEE Trans. Control Syst. Technol. 2016, 24, 1184–1200. [Google Scholar] [CrossRef]
  26. Pagán Rubio, J.A.; Vera-García, F.; Hernandez Grau, J.; Muñoz Cámara, J.; Albaladejo Hernandez, D. Marine diesel engine failure simulator based on thermodynamic model. Appl. Therm. Eng. 2018, 144, 982–995. [Google Scholar] [CrossRef]
  27. Espana, M. On the estimation algorithm for adaptive performance optimization ofturbofan engines. In 29th Joint Propulsion Conference and Exhibit; Joint Propulsion Conferences; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 1993. [Google Scholar]
  28. Ru, J.; Rong Li, X. Variable-Structure Multiple-Model Approach to Fault Detection, Identification, and Estimation. IEEE Trans. Control Syst. Technol. 2008, 16, 1029–1038. [Google Scholar] [CrossRef]
  29. Syverud, E.; Brekke, O.; Bakken, L.E. Axial compressor deterioration caused by saltwater ingestion. J. Turbomach. 2005, 129, 119–126. [Google Scholar] [CrossRef]
  30. Kurz, R.; Brun, K.; Wollie, M. Degradation Effects on Industrial Gas Turbines. J. Eng. Gas Turbines Power 2009, 131, 062401. [Google Scholar] [CrossRef]
  31. Rewienski, M.; White, J. A Trajectory Piecewise-Linear Approach to Model Order Reduction and Fast Simulation of Nonlinear Circuits and Micromachined Devices. IEEE Trans. Comput. Aided Des. Integr. Circuits Syst. 2003, 22, 155–170. [Google Scholar] [CrossRef]
  32. Lin, M.-H.; Tsai, J.-F. Range reduction techniques for improving computational efficiency in global optimization of signomial geometric programming problems. Eur. J. Oper. Res. 2012, 216, 17–25. [Google Scholar] [CrossRef]
  33. El-Sakkary, A. The gap metric: Robustness of stabilization of feedback systems. IEEE Trans. Autom. Control 1985, 30, 240–247. [Google Scholar] [CrossRef]
  34. Georgiou, T.T.; Smith, M.C. Optimal Robustness in the Gap Metric. IEEE Trans. Autom. Control 1990, 35, 673–686. [Google Scholar] [CrossRef]
  35. Yang, Q.; Li, S.; Cao, Y. A new component map generation method for gas turbine adaptation performance simulation. J. Mech. Sci. Technol. 2017, 31, 1947–1957. [Google Scholar] [CrossRef]
  36. Camporeale, S.M.; Fortunato, B.; Mastrovito, M. A Modular Code for Real Time Dynamic Simulation of Gas Turbines in Simulink. J. Eng. Gas Turbines Power 2006, 128, 506. [Google Scholar] [CrossRef]
  37. Yu, H.; Yuecheng, Y.; Shiying, Z.; Zhensheng, S. Comparison of linear models for gas turbine performance. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2014, 228, 1291–1301. [Google Scholar] [CrossRef]
  38. Deshpande, A.P.; Zamad, U.; Patwardhan, S.C. Online Sensor/Actuator Failure Isolation and Reconfigurable Control Using the Generalized Likelihood Ratio Method. Ind. Eng. Chem. Res. 2009, 48, 1522–1535. [Google Scholar] [CrossRef]
Figure 1. Demonstration of the hypothetical model set update.
Figure 1. Demonstration of the hypothetical model set update.
Energies 11 03316 g001
Figure 2. Diagram of the fault contribution matrix (FCM) based multiple model fault detection and isolation (MM-FDI) approach.
Figure 2. Diagram of the fault contribution matrix (FCM) based multiple model fault detection and isolation (MM-FDI) approach.
Energies 11 03316 g002
Figure 3. The layout of a two-shaft marine gas turbine. T 0 is the inlet temperature of the compressor and P 0 is the ambient pressure; T 2 and P 2 are the compressor outlet temperature and outlet pressure, respectively, and N 1 is the compressor speed; w f is the fuel mass flow which sprays into the combustion chamber to produce outlet temperature T 3 and outlet pressure P 3 ; T 4 and P 4 are the exhaust gas temperature and pressure of the compressor turbine, respectively; T 5 and P 5 are the outlet temperature and outlet pressure of the power turbine, respectively, and N 2 is the power turbine speed.
Figure 3. The layout of a two-shaft marine gas turbine. T 0 is the inlet temperature of the compressor and P 0 is the ambient pressure; T 2 and P 2 are the compressor outlet temperature and outlet pressure, respectively, and N 1 is the compressor speed; w f is the fuel mass flow which sprays into the combustion chamber to produce outlet temperature T 3 and outlet pressure P 3 ; T 4 and P 4 are the exhaust gas temperature and pressure of the compressor turbine, respectively; T 5 and P 5 are the outlet temperature and outlet pressure of the power turbine, respectively, and N 2 is the power turbine speed.
Energies 11 03316 g003
Figure 4. FCM based FDI approach simulation testing.
Figure 4. FCM based FDI approach simulation testing.
Energies 11 03316 g004
Figure 5. The matrix of the gap metric between the gas turbine linearized models.
Figure 5. The matrix of the gap metric between the gas turbine linearized models.
Energies 11 03316 g005
Figure 6. Comparison of the response of each measurement parameter between the nonlinear model and the FCM based model set at a given fuel flow schedule, and the corresponding error.
Figure 6. Comparison of the response of each measurement parameter between the nonlinear model and the FCM based model set at a given fuel flow schedule, and the corresponding error.
Energies 11 03316 g006
Figure 7. Performance comparison between the traditional Fmc fault linear model at design operating points and the FCM based model set.
Figure 7. Performance comparison between the traditional Fmc fault linear model at design operating points and the FCM based model set.
Energies 11 03316 g007
Figure 8. FDI result of a single gas path fault under the steady state condition and transient condition.
Figure 8. FDI result of a single gas path fault under the steady state condition and transient condition.
Energies 11 03316 g008
Figure 9. The fault amplitude of the Fηpt fault and the corresponding FDI results of 145 simulations.
Figure 9. The fault amplitude of the Fηpt fault and the corresponding FDI results of 145 simulations.
Energies 11 03316 g009
Figure 10. Performance of the FCM based MM-FDI approach under outliers. (a) The measurement parameters in the scenarios of Fmc decreasing by 1% at t = 20 s; (b) The conditional probability µi of each hypothetical model in the scenarios of Fmc decreasing by 1% at t = 20 s along with outliers.
Figure 10. Performance of the FCM based MM-FDI approach under outliers. (a) The measurement parameters in the scenarios of Fmc decreasing by 1% at t = 20 s; (b) The conditional probability µi of each hypothetical model in the scenarios of Fmc decreasing by 1% at t = 20 s along with outliers.
Energies 11 03316 g010
Figure 11. The multiple fault detection and isolation results for the scenarios of a 2% decrease in Fηpt at t = 15 s and 1% decrease in Fηc at t = 35 s. (a) The conditional probability of the initial model set; (b) The conditional probability of the updated model set 1; (c) The conditional probability of the updated model set 2.
Figure 11. The multiple fault detection and isolation results for the scenarios of a 2% decrease in Fηpt at t = 15 s and 1% decrease in Fηc at t = 35 s. (a) The conditional probability of the initial model set; (b) The conditional probability of the updated model set 1; (c) The conditional probability of the updated model set 2.
Energies 11 03316 g011
Figure 12. The type and amplitude of multiple faults and the corresponding FDI results of 100 stochastic simulations. (a) The type and amplitude of the first abrupt fault and corresponding FDI results; (b) The type and amplitude of the first gradual fault and corresponding FDI results; (c) The type and amplitude of the second abrupt fault and corresponding FDI results; (d) The type and amplitude of the second gradual fault and corresponding FDI results.
Figure 12. The type and amplitude of multiple faults and the corresponding FDI results of 100 stochastic simulations. (a) The type and amplitude of the first abrupt fault and corresponding FDI results; (b) The type and amplitude of the first gradual fault and corresponding FDI results; (c) The type and amplitude of the second abrupt fault and corresponding FDI results; (d) The type and amplitude of the second gradual fault and corresponding FDI results.
Energies 11 03316 g012
Figure 13. The multiple faults detection and isolation results in the scenarios of F w f decreases by 5% fault occurs at t = 15 s, F m c decreases by 2% fault occurs at t = 35 s, and F η p t decreases by 1% fault occurs at t = 45 s. (a) The conditional probability of the initial model set; (b) The conditional probability of the updated model set 1; (c) The conditional probability of the updated model set 2.
Figure 13. The multiple faults detection and isolation results in the scenarios of F w f decreases by 5% fault occurs at t = 15 s, F m c decreases by 2% fault occurs at t = 35 s, and F η p t decreases by 1% fault occurs at t = 45 s. (a) The conditional probability of the initial model set; (b) The conditional probability of the updated model set 1; (c) The conditional probability of the updated model set 2.
Energies 11 03316 g013
Table 1. The hypothetical mode set update algorithm.
Table 1. The hypothetical mode set update algorithm.
Step 1: Determine the type and number of the detected faults
m i = f i n d ( μ > 0.98 ) z i = f i n d ( F ( : , m i ) < 0 ) n = n u m e l ( z i )
Step 2: Determine the matrix Fupdate corresponding to the updated hypothetical model set
F u p d a t e = F ; f o r   k = 1 : n       F u p d a t e ( z i ( k ) , : ) = F ( z i ( k ) , m i ) ;       F u p d a t e ( z i ( k ) , z i ( k ) + 1 ) = F ( z i ( k ) , m i ) s i ; e n d
Step 3: Substituting each column of Fupdate into Equation (4) to obtain the updated mode set
{ x k + 1 = A x k + [ B E ] [ u k f k ] y k = C x k + [ 0 D ] [ u k f k ]
Notation
μ denotes the conditional probability of the W hypothetical models, a 1 × W vector;
F   denotes   a   matrix   which   is   composed   of   f k   corresponding   to   each   hypothetical   model   in   the   current   hypothetical   model   set ,   F = [ f k 1 # f k 2 # f k W # ] ,   a   q × W matrix;
F u p d a t a   denotes   a   matrix   which   is   composed   of   f k   corresponding   to   each   hypothetical   model   in   the   updated   model   set ,   a   q × W matrix;
s denotes the fault amplitude changes during model set updating, s = −1.
Table 2. Typical single gas path faults.
Table 2. Typical single gas path faults.
Fault DescriptionSymbolFault Vector
Healthy condition H f k = [ 0   0   0   0   0   0 ] T
Decrease   in   m C F m c ( b 1 % ) f k = [ b 1   0   0   0   0   0 ] T ,   b 1 = ( m C / m C , H 1 ) × 100 %
Decrease   in   η C F η c ( b 2 % ) f k = [ 0   b 2   0   0   0   0 ] T ,   b 2 = ( η C / η C , H 1 ) × 100 %
Decrease   in   m C T F m c t ( b 3 % ) f k = [ 0   0   b 3   0   0   0 ] T ,   b 3 = ( m C T / m C T , H 1 ) × 100 %
Decrease   in   η C T F η c t ( b 4 % ) f k = [ 0   0   0   b 4   0   0 ] T ,   b 4 = ( η C T / η C T , H 1 ) × 100 %
Decrease   in   m P T F m p t ( b 5 % ) f k = [ 0   0   0   0   b 5   0 ] T ,   b 5 = ( m P T / m P T , H 1 ) × 100 %
Decrease   in   η P T F η p t ( b 6 % ) f k = [ 0   0   0   0   0   b 6 ] T ,   b 6 = ( η P T / η P T , H 1 ) × 100 %
Table 3. Sensor noise.
Table 3. Sensor noise.
Measurement ParametersN1N2T2P2T4P4T5
Standard deviation (%)0.0510.0510.230.1640.0970.1640.097
Table 4. The result of operating points selection based on the gap metric analysis.
Table 4. The result of operating points selection based on the gap metric analysis.
Models1#2#3#4#5#6#
i1217426371
wf10.990.840.590.380.3
Table 5. The fault detection and isolation (FDI) results of the 1000 stochastic simulations.
Table 5. The fault detection and isolation (FDI) results of the 1000 stochastic simulations.
Confusion Matrix of the FDI ResultFinal Result
HFmcFηcFmctFηctFmptFηptCDIDMD
H135000000
Fmc01230000095.5%0.45%0
Fηc001520000
Fmct000121000
Fηct000015300
Fmpt000001710
Fηpt0450000145
Table 6. The detection time and the isolation time for different numbers of measurements available.
Table 6. The detection time and the isolation time for different numbers of measurements available.
Fault TypeSeven SensorsFour SensorsTwo Sensors
td (s)ti (s)td (s)ti (s)td (s)ti (s)
Fmc1.465.024.817.652231009
Fηc0.93.684.7851.4MD/IDMD/ID
Fmct0.280.90.75.029123739
Fηct4.3217.92.9262.8MD/IDMD/ID
Fmpt0.341.25.1225.7MD/IDMD/ID
Fηpt0.541.82.088.522.5116.2
Table 7. The FDI results of 100 stochastic simulations in the scenarios of multiple faults.
Table 7. The FDI results of 100 stochastic simulations in the scenarios of multiple faults.
Fault ChangeFaultsSimulation TimesFDI Result
CDIDMD
AbruptFirst fault10010000
Second fault1009451
Total1009451
GradualFirst fault10010000
Second fault1009550
Total1009550

Share and Cite

MDPI and ACS Style

Yang, Q.; Li, S.; Cao, Y.; Gu, F.; Smith, A. A Gas Path Fault Contribution Matrix for Marine Gas Turbine Diagnosis Based on a Multiple Model Fault Detection and Isolation Approach. Energies 2018, 11, 3316. https://doi.org/10.3390/en11123316

AMA Style

Yang Q, Li S, Cao Y, Gu F, Smith A. A Gas Path Fault Contribution Matrix for Marine Gas Turbine Diagnosis Based on a Multiple Model Fault Detection and Isolation Approach. Energies. 2018; 11(12):3316. https://doi.org/10.3390/en11123316

Chicago/Turabian Style

Yang, Qingcai, Shuying Li, Yunpeng Cao, Fengshou Gu, and Ann Smith. 2018. "A Gas Path Fault Contribution Matrix for Marine Gas Turbine Diagnosis Based on a Multiple Model Fault Detection and Isolation Approach" Energies 11, no. 12: 3316. https://doi.org/10.3390/en11123316

APA Style

Yang, Q., Li, S., Cao, Y., Gu, F., & Smith, A. (2018). A Gas Path Fault Contribution Matrix for Marine Gas Turbine Diagnosis Based on a Multiple Model Fault Detection and Isolation Approach. Energies, 11(12), 3316. https://doi.org/10.3390/en11123316

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop