Gas Turbine Transient Performance Tracking Using Data Fusion Based on an Adaptive Particle Filter

This paper considers the problem of gas turbine transient performance tracking in a cluttered environment. To increase the accuracy and robustness of state estimation, a data-fusion nonlinear estimation method based on an adaptive particle filter (PF) is proposed. This method needs local estimates transmitted to a central filtering unit for data fusion, and then global data feedback to the local PF for consensus propagation. The computational burden is shared by the local PF and central filtering unit in the data-fusion architecture. Furthermore, the PF algorithm used for the data fusion is embedded with the prior knowledge of engine health condition and adaptive to the measurement noise, and hence is called the adaptive PF. The heuristic information of state variables represented by inequality constraints tunes the local estimates by a probability density truncation method. The covariance of measurement noise is calculated by wavelet transform and utilized to update the particle importance function of the real time PF. The performance improvements of the proposed method are indicated through extensive experiments for gradual and abrupt shift performance tracking under conditions of gas turbine transient operation.


Introduction
Gas turbine engines provide the power for airplanes, and their reliability is vital to flight safety and performance. Nevertheless, the engine working conditions are terrible and they usually must endure high speeds, extreme temperatures and strong vibrations. Erosion and fouling of major components are unavoidable and result in a gradual deterioration of engine performance during its lifetime. Besides, foreign/domestic object damage will cause engine performance to sharply shift, and it is called an abrupt fault [1]. Engine data are periodically collected by airlines for engine health evaluation. Maintenance schedules adapted to reliable health tracking results leads to safe operation and reduced costs [2]. Hence, how to get reliable information about engine health conditions in time has drawn a lot of attention.
Engine health parameters, which are correction factors for the efficiency of the major components, are introduced to quantify the engine performance changes [3]. The engine performance tracking problem can be regarded as calculating the health parameters. Various methods such as Kalman filters (KFs), neural networks, fuzzy logic, genetic algorithms and expert systems have been proposed to obtain health parameters for engine health monitoring [3][4][5][6][7]. KF-based methods seem to be the most common ones for gas turbine health estimation, but these techniques are mainly focused on engine steady state, e.g., under cruise or average conditions [8].
The engine transient condition represents how a gas turbine operates from one steady state to another steady state when the input changes. The rotating components of the engine (fan, compressor, high-pressure turbine (HPT) and low-pressure turbines (LPT)) run close to their surge boundaries during transient behavior, and some performance parameters exceed their extreme value in a short time. Generally speaking, a low-bypass engine performance anomaly more easily occurs during a transient process, e.g., acceleration and deceleration. That is to say, gas turbine transient performance tracking is more urgent compared to steady state monitoring, but up to this point in time no studies of dynamic behavior monitoring for gas turbines have been presented. This paper proposes an adaptive particle filter (PF) based data fusion approach to engine transient health monitoring.
In this paper we emphasize the problem of monitoring the health performance of key gas turbine rotating components, including fans, compressors (high-pressure compressors (HPCs)), high-pressure turbines (HPTs) and low-pressure turbines (LPTs). Failures of engine actuators, sensors, and other components are not considered. The engine performance tracking of gradual drifts and abrupt shifts during transient operation is mainly addressed, but the technique can also be applied to engine steady behavior. A data fusion architecture to monitor engine health conditions is designed, and prior knowledge of engine health state and measurement noise is utilized for a Monte Carlo simulation to tune the state estimates. The present contribution is a derivation of a data fusion method based on adaptive PF for gas turbine transient performance tracking at the theoretical, implementation and performance level.
This paper is organized as follows: Section 2 presents a review of the basic PF and the problem formulation; the adaptive PF algorithm and data fusion architecture implementation for gas turbine health monitoring is given in Section 3; the performance comparisons of the proposed method for gradual and sharp shifts in different operation conditions are discussed in Section 4; and Section 5 concludes this paper.

Problem Formulation
A gas turbine engine is a low-bypass turbofan engine, see Figure 1. A single inlet supplies airflow to the fan. Air leaving the fan is separated into two streams: one stream passes through the engine core, and the other stream passes through the annular bypass duct and then leaves. The fan is driven by the low pressure turbine. The air passing through the engine core moves through the HPC, which is driven by the HPT. Fuel is injected in the combustor and burned to produce hot gas for driving the turbines. The gas leaves the LPT and is mixed with the air from the bypass duct through the convergent nozzle, which has a variable cross section area.
Energies 2015, 8, page-page 2 be the most common ones for gas turbine health estimation, but these techniques are mainly focused on engine steady state, e.g., under cruise or average conditions [8].
The engine transient condition represents how a gas turbine operates from one steady state to another steady state when the input changes. The rotating components of the engine (fan, compressor, high-pressure turbine (HPT) and low-pressure turbines (LPT)) run close to their surge boundaries during transient behavior, and some performance parameters exceed their extreme value in a short time. Generally speaking, a low-bypass engine performance anomaly more easily occurs during a transient process, e.g., acceleration and deceleration. That is to say, gas turbine transient performance tracking is more urgent compared to steady state monitoring, but up to this point in time no studies of dynamic behavior monitoring for gas turbines have been presented. This paper proposes an adaptive particle filter (PF) based data fusion approach to engine transient health monitoring.
In this paper we emphasize the problem of monitoring the health performance of key gas turbine rotating components, including fans, compressors (high-pressure compressors (HPCs)), high-pressure turbines (HPTs) and low-pressure turbines (LPTs). Failures of engine actuators, sensors, and other components are not considered. The engine performance tracking of gradual drifts and abrupt shifts during transient operation is mainly addressed, but the technique can also be applied to engine steady behavior. A data fusion architecture to monitor engine health conditions is designed, and prior knowledge of engine health state and measurement noise is utilized for a Monte Carlo simulation to tune the state estimates. The present contribution is a derivation of a data fusion method based on adaptive PF for gas turbine transient performance tracking at the theoretical, implementation and performance level. This paper is organized as follows: Section 2 presents a review of the basic PF and the problem formulation; the adaptive PF algorithm and data fusion architecture implementation for gas turbine health monitoring is given in Section 3; the performance comparisons of the proposed method for gradual and sharp shifts in different operation conditions are discussed in Section 4; and Section 5 concludes this paper.

Problem Formulation
A gas turbine engine is a low-bypass turbofan engine, see Figure 1. A single inlet supplies airflow to the fan. Air leaving the fan is separated into two streams: one stream passes through the engine core, and the other stream passes through the annular bypass duct and then leaves. The fan is driven by the low pressure turbine. The air passing through the engine core moves through the HPC, which is driven by the HPT. Fuel is injected in the combustor and burned to produce hot gas for driving the turbines. The gas leaves the LPT and is mixed with the air from the bypass duct through the convergent nozzle, which has a variable cross section area.  Considering air flow mass, power and momentum conservation laws [9,10], a nonlinear aero-thermodynamic model of a gas turbine engine is given by: where k is the time index, y is the 8-element measured output, x is the 6-element augmented state, and u is the 2-element control input. The noise terms w k and v k represent the process inaccuracies and measurement inaccuracies in the model. The sensor measurements are low-pressure spool speed N L , high-pressure spool speed N H , fan outlet pressure P 22 , HPC outlet pressure P 3 , fan outlet temperature T 22 , HPC outlet temperature T 3 , HPT outlet temperature T 43 and LPT outlet temperature T 6 . The augmented state x includes the 2-element original state x o (N L and N H ) and health parameter vector p (fan efficiency SE1, HPC efficiency SE2, HPT efficiency SE3 and LPT efficiency SE4, where section efficiency is defined SE). The elements of the control vector in the model are fuel flow W f and nozzle area A 8 , which defines the engine operating point. There are two information entropy definitions used to select the system parameters. Auto Information Entropy is utilized to select the measured parameters, while the Cross Information Entropy to analyze the correlations between control, measured and health parameters.
In the framework of gas turbine performance monitoring, the quantities of interest are the differences between the estimated engine health parameters status and their reference ones. The actual engine performance is represented by the estimated values of health parameter, the prior value of which is adapted to the current measurements in a recursive approach. The block diagram of the model-based approach to monitoring the engine health condition is shown in Figure 2, and it is a closed-loop state estimator correcting structure. The KF, especially the linear KF (LKF), is an optimal state estimator for linear systems with noisy and inaccurate measurements, and is widely used for gas turbine engine health monitoring [11]. Considering air flow mass, power and momentum conservation laws [9,10], a nonlinear aero-thermodynamic model of a gas turbine engine is given by: where k is the time index, y is the 8-element measured output, x is the 6-element augmented state, and u is the 2-element control input. The noise terms k w and k v represent the process inaccuracies and measurement inaccuracies in the model. In the framework of gas turbine performance monitoring, the quantities of interest are the differences between the estimated engine health parameters status and their reference ones. The actual engine performance is represented by the estimated values of health parameter, the prior value of which is adapted to the current measurements in a recursive approach. The block diagram of the model-based approach to monitoring the engine health condition is shown in Figure 2, and it is a closed-loop state estimator correcting structure. The KF, especially the linear KF (LKF), is an optimal state estimator for linear systems with noisy and inaccurate measurements, and is widely used for gas turbine engine health monitoring [11]. The state variable model (SVM) is a piecewise linear representation of the engine and needs to be established before the LKF implementation. It should be pointed out that the SVM obtained at a steady operating point is only a representation of the engine near this operating point. The state estimation accuracy by the LKF varies with the distance from the actual operating point to the design point, and also with the magnitude of the state deviation [12]. That means that the strong nonlinearity of gas turbine makes the LKF effective only in a small working range around the nominal state, and it can't do well in cases of transient behavior, e.g., during acceleration and deceleration. The state variable model (SVM) is a piecewise linear representation of the engine and needs to be established before the LKF implementation. It should be pointed out that the SVM obtained at a steady operating point is only a representation of the engine near this operating point. The state estimation accuracy by the LKF varies with the distance from the actual operating point to the design point, and also with the magnitude of the state deviation [12]. That means that the strong nonlinearity of gas turbine makes the LKF effective only in a small working range around the nominal state, and it can't do well in cases of transient behavior, e.g., during acceleration and deceleration.
Note that state variables change together in a situation of normal gradual engine deterioration, and the SVM is derived by the partial derivative calculation in the linearization process. Thus it is impossible for the SVM to correctly describe simultaneous shifts of all states, thus the tracking accuracy of gradual shifts by the LKF is not satisfactory. In addition, the assumption that the derivative of health parameter h is approximated by zero is utilized in the LKF equations, which is because engine deterioration generally occurs very slowly relative to the dynamics of the original state variables x 0 . Therefore any sharp shift of health parameters is not well estimated by the LKF due to the contradiction with the hypothesis . h " 0. Hence, the generic LKF-based engine monitoring approach to the abrupt and multi-state gradual shift has a theoretical drawback.
To aim at the above issues, various nonlinear filtering methods are proposed and applied to state estimation, especially the extended Kalman filter (EKF), unscented Kalman filter (UKF) and PF [13,14]. Nevertheless, the EKF method is often used for weak nonlinear Gaussian systems due to a Taylor Series expansion truncated to the first order [8], and the performance closely depends on how often Jacobians are updated. Compared to the EKF and UKF, the PF is based on sequential Monte Carlo sampling theory, and it does not necessitate simplification of nonlinearity or any hypothesis of specific distributions [15,16]. Therefore, the PF-based approach is the one of the best ways that we address strongly nonlinear engine health monitoring issues in the following section.

The Particle Filter
The nonlinear model of a gas turbine is given by Equtaion (1), and let x 0:k fi tx 0 ,¨¨¨, x k u and y 1:k fi ty 1 ,¨¨¨, y k u denote the series of the state and measurement. Assume the probability density function (PDF) of the prior condition is ppx 0 q, and the posterior PDF ppx 0:k |y 1:k q is characterized by a set of weighted random samples x i 0:k , w i k ( N i"1 , wherein N s is the particle number. The particle set x i 0:k ; i " 0,¨¨¨, N ( is associated to the weights w i k ; i " 0,¨¨¨, N ( , and the PDF at time k can be approximated by: The case that the particles for Monte Carlo sampling are directly generated from the posterior PDF ppx 0:k |y 1:k q is expected, but it is generally unavailable. Thus the importance sampling distribution function qpx 0:k |y 1:k q is defined before sampling: Then the i th particle weight w i k can be approximated by: ppx 0:k |y 1:k q qpx 0:k |y 1:k q 9 With this choice and normalization, the importance weights can be computed: One problem of the basic PF algorithm is that more particles have negligible weights after a few recursive steps. This implies that particle degeneracy occurs and a large computational effort for updating particle is meaningless. A typical method for solving this issue is importance re-sampling, and each particle is assigned by equal weight w i k " 1{N whenever the effective number N eff of particles becomes less than a threshold value N th .
Once N th is close to value N eff,k , all particles have almost the same significance. The architecture of a conventional PF mentioned above is that the measured data from different sensors are sent to one central PF to process together, and this is so-called the central architecture [15,17]. The advantage of this architecture is minimal information loss, but it also raises the problems that all measurements are identically treated at a time and the central filter bears a heavy computational burden, especially in the framework of the PF [18][19][20]. With the development of information fusion, the fusion PF structure employing a bank of local PFs and one master filter is presented and the computational effort is then shared by several filters [21,22].
In the case of the model-based approach to gas turbine health monitoring, the heuristic health knowledge is usually ignored and not considered in estimate algorithms. The typical anomaly modes and gradual deterioration rules of the engine have been summarized, and the magnitude ranges of health parameters are then determined. This prior knowledge about health parameters is represented by the constraints, and it can be used for state estimation in the PF algorithm. Furthermore, measurement noise levels vary with the order of the sensed value during the engine's dynamic operation. Generally speaking, the magnitude of sensed noise increases as the engine works at a larger operating point during the dynamic operation. Hence, it is very important for the PF to tune the estimate with the noise level in the situation of the transient behavior of the engine.

Particle Filter with Inequality Constraints
The prior state information can be described by equality or inequality constraints in the state estimation algorithm. These constraint approaches are mainly concentrated and applied to linear systems, especially combined with a KF, which has been proved to increase performance estimation accuracy. KFs with equality constraint approaches include the moving horizon estimation and smoothly constrained KF [23]. The estimate projection and truncation methods are incorporated to inequality constraints on state estimates, and the latter one has better performance for deterioration tracking due to its handling of two-sided constraints. The nonlinear state estimation with inequality constraints is more useful to gas turbine engine health monitoring, but there are not enough theoretical studies. For a nonlinear dynamic system, the heuristic knowledge represented by inequality constraints is introduced to the PF algorithm. This idea based on a previously published method [24] has been extended to the nonlinear PF approach for abrupt shift tracking.
The prior information we used in this paper is derived from the gradual deterioration rule and abrupt fault feature. Engine performance degrades with use, and the health parameters never improve and change in one way. For example, we know that the health parameter representing engine efficiency declines over time, and the parameter usually varies within the´10% magnitude due to abrupt faults. With this information, we can determine the bound of the state in advance and it is represented by the inequality constraint of the state estimator. The probability density function (PDF) of the PF estimate at the prior inequality constraint is truncated and the constrained filter estimate as the mean of the truncated PDF is calculated. Now consider the nonlinear system of Equation (1) where the s scalar constraints are added: where a k,m ď b k,m , a k,m and b k,m are the two sided constraint of the m th health parameter at time k. The health parameter estimatex k and covariance P k is derived by the unconstraint PF algorithm: The PF algorithm with inequality constraints is to truncate the unconstrained PDF, N (x k , P k ). Once the mean r x k and covariance r P k of the truncated PDF are computed, we can obtain the constrained health parameter estimate. The state estimate r x k,m and covariance r P k,m are defined after the first m scalar constraints enforced, and then the new transformation is performed: where T m is orthogonal, W m is diagonal (these two quantities, T m and W m , can be derived from the Jordan canonical decomposition of r P k,m ) and S m is obtained using Gram-Schmidt orthogonalisation. The first m inequality constraints are normalized: Since z k,m is an identity covariance with statistically independent element, only the first element is constrained in Equation (10), and the PDF truncation reduces to a one-dimensional PDF. The part outside of the constraints is removed due to the fact z k,m lays between c k,m and d k,m . The truncated PDF is normalized after computing the area of the inside portion of the PDF. The mean µ m and variance σ 2 m of the first element of r z k,m with the constraint enforcement are expressed by: where α is a magnification factor. The inverse transformation of the Equation (9) is taken, and the mean and variance of the constrained state estimate are therefore given as: We repeat the process of Equations (9)- (12) to enforce the next constraint to the state estimate and jump out of the iteration until m = s. Hence, we obtain the state estimate and covariance at time k as x k " r x k,s , P k " r P k,s by the constrained PF algorithm.

Measure Noise Tuned Particle Filter
As we know, particle importance weight w i is closely dependent on the distance between the real sensed value and its estimate, and the smaller the distance, the larger the weight. In the basic PF algorithm, the particle importance weight w i is defined as follows: where the quantity σ is consistent with the covariance of measurement noise due to the fact that the sensed noise usually is generally recognized as Gaussian white noise. We set the quantity σ by the heuristic experience of the noise and it remains constant in the conventional PF algorithm. However, the probability distribution of the sensed noise changes in engine transient behavior, and it raises the problem of measurement noise uncertainty with regards to state estimation. The deviation from the actual measurement noise will result in a decline in estimation accuracy. In the paper the on-line tuning quantity σ of the PF with the measurement noise variance with the help of the wavelet transform is designed. The measured stream y ω pkq could be approximated by a low-order polynomial or a piecewise low-order polynomial in an observation interval [25]. This interval size is corresponding to sampling number. Suppose the polynomial of the sensed stream y ω pkq is described as follows: Then the wavelet transform is implemented to the sensed stream m w pkq: W y ω ps, τq " y ω pkq˚ϕ s,τ pkq " W y ps, τq`W ω ps, τq where φpkq is a wavelet function, "*" is the convolution operation, W y ps, τq is an approximate part of the wavelet coefficients and W ω ps, τq represents a detailed part of the wavelet coefficients. We define a vanishing moment α of the function φpkq, and the measurement noise part is extracted from the sensed stream part once the quantity α is greater than the highest-order of the polynomial, namely α > M. Then the standard deviation of the noisy term in an interval can be calculated by the following expression: where the scale s equals to 0.5, and the series W y ω ps, t h q is the K/2 wavelet coefficients of ty ω pkq |k " 0, 1,¨¨¨K u, wherein 0 ď h ď K{2. Normalizing parameter M is usually set at 0.6745, and the function Med() represents the middle value calculation of a series. In order to tune the quantity σ of the PF algorithm in real-time during the engine dynamic behavior, the measurement noise of a series in an interval is computed by a wavelet transformation. The sensed stream y ω is segmented by a window with fixed width L, and it slides forward along the sampling time. The variance of the sensed noise part in the sliding window is estimated, therefore, the framework of the adaptive PF algorithm is established as shown in Figure 3.
,   Figure 3. The framework of the adaptive particle filter (PF) algorithm.

Data Fusion Based on Adaptive Particle Filter
The fusing PF based on the integration of adaptive PF algorithm and information fusion theory is proposed for gas turbine transient performance monitoring. The sensor used to detect engine health condition is divided into several teams, and each team is applied to estimate the local estimates by the local PF. The engine component layout and thermodynamics of operation are taken into account for the partitions of local filters in the fusion architecture. The fan of the engine is a cold and low pressure (LP) component, the compressor is a cold and high pressure (HP) one, HPT is a hot and HP one, and LPT is a hot and LP one [26,27]. Two different ways to partitioning are performed, namely, the cold-hot component partition and the LP-HP component partition. Therefore, there are four combinations of the four rotating components above, which are the cold group, the hot group, the HP group and the LP group. For example, the cold partition includes two components (fan and compressor). Likewise, the hot partition consists of HPT and LPT, the HP partition of compressor and HPT, and the LP partition of fan and LPT.
For the cold group, the measured parameters are T22, P22, T3 and P3. In a similar manner the measurements T43 and T6 are for the hot group, the sensors T22, P22 and T6 for the LP group, and the sensors T3, P3, and T43 for the HP group. In addition, two spool speeds (NL and NH) are important quantities representing engine operation status, which are communal measurements and utilized in each part. Hence, the measurements of four local filter can be denoted by The data fusion estimation based on the adaptive PF mainly includes three stages. First, several local filters perform in parallel to obtain individual sensor-based estimates. Second, all local estimates are combined in a master filter, where a global state estimate is yielded. Third, the global state and covariance is fed back to local filters with an information-sharing strategy for next cycle. The procedure of data fusion for the engine transient performance estimation is summarized as follows: Step 1: Initialization Given the initial values of global state 0 x , estimation error covariance 0 P , and process noise covariance 0 Q in the master filter. The four local filters are initialized with the information allocation strategy: Figure 3. The framework of the adaptive particle filter (PF) algorithm.

Data Fusion Based on Adaptive Particle Filter
The fusing PF based on the integration of adaptive PF algorithm and information fusion theory is proposed for gas turbine transient performance monitoring. The sensor used to detect engine health condition is divided into several teams, and each team is applied to estimate the local estimates by the local PF. The engine component layout and thermodynamics of operation are taken into account for the partitions of local filters in the fusion architecture. The fan of the engine is a cold and low pressure (LP) component, the compressor is a cold and high pressure (HP) one, HPT is a hot and HP one, and LPT is a hot and LP one [26,27]. Two different ways to partitioning are performed, namely, the cold-hot component partition and the LP-HP component partition. Therefore, there are four combinations of the four rotating components above, which are the cold group, the hot group, the HP group and the LP group. For example, the cold partition includes two components (fan and compressor). Likewise, the hot partition consists of HPT and LPT, the HP partition of compressor and HPT, and the LP partition of fan and LPT.
For the cold group, the measured parameters are T 22 , P 22 , T 3 and P 3 . In a similar manner the measurements T 43 and T 6 are for the hot group, the sensors T 22 , P 22 and T 6 for the LP group, and the sensors T 3 , P 3 , and T 43 for the HP group. In addition, two spool speeds (N L and N H ) are important quantities representing engine operation status, which are communal measurements and utilized in each part. Hence, the measurements of four local filter can be denoted by y The data fusion estimation based on the adaptive PF mainly includes three stages. First, several local filters perform in parallel to obtain individual sensor-based estimates. Second, all local estimates are combined in a master filter, where a global state estimate is yielded. Third, the global state and covariance is fed back to local filters with an information-sharing strategy for next cycle. The procedure of data fusion for the engine transient performance estimation is summarized as follows: Step 1: Initialization Given the initial values of global state x 0 , estimation error covariance P 0 , and process noise covariance Q 0 in the master filter. The four local filters are initialized with the information allocation strategy: Q j,0 " β´1 j Q m,0 P j,0 " β´1 j P m,0 X j,0 " β´1 j X m,0 j " 1,¨¨¨, 4 where the information distribution factor β j follows: Step 2: The adaptive PF performs in the local filter.
The particles tx i j,0:k´1 , w i j,k´1 u N i"1 are generated based on the prior distribution NpX j,k´1 , P j,k´1 q, and are propagated through the nonlinear model Equation (1). The numerical characteristics of the measured noise are computed by wavelet transformation, and the quantity σ is applied to the importance weight calculation. The known health information is imposed to the PF to produce the local constrained estimates of the statex j,k , error covarianceP j,k and noise covariance Q j,k .
Step 3: Information fusion implements in the master filter.
The local estimatesx j,k ,P j,k , and Q j,k are sent to the master filter to fulfill the information fusion for the global optimal estimate x m,k . The estimate error covarianceP j,k is an important parameter representing the performance of the local filter, and it is used to calculate the fusing weight of the local filter. The larger the covarianceP j,k is, the smaller the fusing weight in the global state is in the paper: Step 4: Information distribution strategy The state estimate calculated by the master filter is transmitted back to each local filter with an information assignment strategy: Q j,k " β´1 j Q m,k P j,k " β´1 j P m,k X j,k " X m,k j " 1¨¨¨N (20) Due to the fact that information distribution factor β j has no effect on estimation accuracy, it is set by β j " 1 N .
Steps (2)-(4) present the fusing PF algorithm at iteration k. For iteration k + 1, the state and covariance delay a time index and Steps (2)-(4) are repeated. The fusion filter architecture is shown in Figure 4 for gas turbine transient performance estimation.
As shown in the figure, we can see that calculation loads are shared both by the local filters, and the master filter no longer undertakes the whole process like in the basic PF algorithm. Since the time update and measurement update are carried out independently in every local filter, the individual estimate by the local filter is not immediately affected by others. In the gas turbine engine health monitoring application, the data fusion filter architecture has a more efficient capability to deal with state estimation in cases of sensor fault due to fusing weight adaptive to estimation accuracy of the local filter.

Simulation and Analysis
The data fusion based on the adaptive PF approach is evaluated for the engine performance monitoring using the Matlab software. In the simulation environment, we use the component-level-model (CLM) engine to take the place of the actual engine, and the sampling rate equals to 50 Hz. The hardware of computer used for simulation is configured as follows: CPU i3-2100 @ 3.10 GHz and RAM 2GB. The standard deviations of measurement and health parameter are shown in Tables 1 and 2. Gaussian noise v with magnitude specified in Table 1 is added to the simulated measured values, and the independent system noise and initial measured noise separately follow Gas turbine health condition is represented by health parameters as mentioned in the previous section. Table 2 shows four abrupt faults representative of possible situations expected to be encountered in practice, and the health parameter deviation in each case refers to lab record of the Rolls-Royce Company (London, UK) [3,28]. For example, Case 1 is a fan abrupt fault with its efficiency SE1 deviating 1%. Assume that there are no fault on sensor measurements in the following experiments.

Simulation and Analysis
The data fusion based on the adaptive PF approach is evaluated for the engine performance monitoring using the Matlab software. In the simulation environment, we use the component-level-model (CLM) engine to take the place of the actual engine, and the sampling rate equals to 50 Hz. The hardware of computer used for simulation is configured as follows: CPU i3-2100 @ 3.10 GHz and RAM 2GB. The standard deviations of measurement and health parameter are shown in Tables 1 and 2. Gaussian noise v with magnitude specified in Table 1 is added to the simulated measured values, and the independent system noise and initial measured noise separately follow ω " Np0, Qq and v " Np0, Rq, wherein Q " 0.16ˆ10´4 I 6ˆ6 . Gas turbine health condition is represented by health parameters as mentioned in the previous section. Table 2 shows four abrupt faults representative of possible situations expected to be encountered in practice, and the health parameter deviation in each case refers to lab record of the Rolls-Royce Company (London, UK) [3,28]. For example, Case 1 is a fan abrupt fault with its efficiency SE1 deviating 1%. Assume that there are no fault on sensor measurements in the following experiments. Engine gradual deterioration due to normal usage is simulated by linear drift of four health parameters, beginning from a healthy engine (the four parameters equal to 1) at cycle number n = 0 and with the degeneration at the end of the sequence at n = 6000:´2.18% on SE1,´6.71% on SE2, 3.22% on SE3 and´0.81% on SE4. Considering the magnitude of both the engine abrupt fault and deterioration, the bounds of the health parameter representing the PF inequality constraints are separately set by a = [1.005,1.005,1.005,1.005] T and b = [0.97,0.90,0.96,0.98] T . The performance of the engine anomaly detection is assessed by three indices, namely, root-mean-square error (RMSE), convergence time and root-mean-square deviation (RMSD). The RMSE and RMSD are separately defined by: where S is the sampling step andx i the mean of estimate value. The convergence time T c is used to indicate the delay in fault recognition. We define this time index in this paper that is from the starting deviation to the estimate steady state within˘0.02% range and no longer out of this range in two consecutive steps.

Abrupt Fault Diagnosis in Steady Operation Conditions
The tests on gas turbine abrupt fault diagnosis are first performed at ground steady conditions (H = 0 m, Ma = 0, W f = 2.48 kg/s). The abrupt faults depicted in Table 2 are simulated, and the noise is not changed in the engine steady behavior. Given the stochastic character of the measurement noise, each test-case has been run five times. Then the engine fault diagnostic performances of basic KF, basic PF, fusion particle filter (F-PF) and fusion adaptive particle filter (FA-PF) are given in Table 3. The particle number of the PF is 60, and that of both the F-PF and FA-PF is 30.
The RMSEs and RMSDs of the two fusion PF approaches shown in Table 3 are smaller than those of the conventional PF and KF in the cases of the four abrupt fault modes, and the FA-PF one is superior to the others. The convergence time of the three PF approaches are nearly the same, and vary clearly in different cases. The KF is a linear estimator and it takes less time to reach the steady state. Two speed measurements are repeatedly utilized in each local filter and there is the reason for the fact the importance of their weights is increasing compared to the remaining measurements in the fusion filter structure. The fusion PF approaches have more satisfactory estimation accuracy due to sufficient extraction and information-sharing of key measurements such as the speeds for health monitoring. Because of the heuristic health knowledge enforced through inequality constraints, the FA-PF has the best estimation accuracy of fault diagnosis in the engine steady operation. Table 3. The engine fault diagnostic performances of four filtering approaches in steady behavior. Kalman filter: KF; fusion particle filter: F-PF; fusion adaptive particle filter: FA-PF.

Abrupt Fault Diagnosis in Dynamic Operation
In order to further evaluate the proposed method performance in engine transient performance tracking, more simulation is carried out in the case of mixed gradual deterioration and abrupt faults. Gradual deterioration refers to all health parameters degrading linearly from the healthy condition to the end of 3000 cycles, and a simulated abrupt fault of magnitude´1% is added to SE2 at 5 s. The engine operates in the dynamic behavior mode from W f = 2.48 kg/s to 1.98 kg/s under three operation conditions: ground (H = 0 m, Ma = 0), high-altitude 1 (H = 8000 m, Ma = 0.5), and high-altitude 2 (H = 11000 m, Ma = 0.8). The particle number of the three PFs is set the same as previously, and the measurement noises remain unchanged. A comparison of the three approaches for the engine transient performance tracking in the cases of mixing gradual deterioration and abrupt fault at ground is depicted in Figure 5, where the dotted line and solid line are the real and estimated values of health parameters, respectively. It can readily be found from the Figure 5 that the PF working in fusion architecture outperforms the basic PF. Table 4 further presents the estimation performance by the three PF methods in terms of number and data at different operating conditions. As can be seen from Table 4, the RMSE of the F-PF and FA-PF is smaller than that of PF, among which the FA-PF in three operation conditions are below 0.007. Hence, the estimation accuracy of fusion PF architecture outperforms that of the generic PF structure.

Performance Estimation with Uncertain Noise in Dynamic Operation
The stochastic feature of the engine measured noise changes at different operation points. An experiment of gradual deterioration tracking with uncertain measurement noise is performed to It can readily be found from the Figure 5 that the PF working in fusion architecture outperforms the basic PF. Table 4 further presents the estimation performance by the three PF methods in terms of number and data at different operating conditions. As can be seen from Table 4, the RMSE of the F-PF and FA-PF is smaller than that of PF, among which the FA-PF in three operation conditions are below 0.007. Hence, the estimation accuracy of fusion PF architecture outperforms that of the generic PF structure.

Performance Estimation with Uncertain Noise in Dynamic Operation
The stochastic feature of the engine measured noise changes at different operation points. An experiment of gradual deterioration tracking with uncertain measurement noise is performed to assess the FA-PF algorithm. The engine experiences dynamic operation in the case of gradual deterioration as the same as the Section 4.2, but no abrupt fault is added. A series of the tests are implemented with ground conditions, including the change of only one measurement noise and all measurement noises simultaneously. Noise generated in the engine core, by sources such as the HPC, combustor, HPT and LPT, plays the most important roles to the overall noise under low-power conditions. While jet and fan noises have dominated over core noise at high engine power during takeoff [29].
The sensor noise varies with engine power condition in transient process. The uncertain noise of one sensor P 22 in Figure 6  assess the FA-PF algorithm. The engine experiences dynamic operation in the case of gradual deterioration as the same as the Section 4.2, but no abrupt fault is added. A series of the tests are implemented with ground conditions, including the change of only one measurement noise and all measurement noises simultaneously. Noise generated in the engine core, by sources such as the HPC, combustor, HPT and LPT, plays the most important roles to the overall noise under low-power conditions. While jet and fan noises have dominated over core noise at high engine power during takeoff [29].
The sensor noise varies with engine power condition in transient process. The uncertain noise of one sensor P22 in Figure 6  As shown in Figure 6, the engine transient performance estimates by the fusion PF approaches seem to deviate from their real values once the noise of sensor P22 varies, no matter the quantity of measurement noise σ is equal to R0, R1 or R2. The fusion adaptive PF has sound tracking performance due to the quantity σ adaptive to the real measurement noise real time in Figure 6d. We can obtain similar results in the case that all sensors for transient performance monitoring have uncertain measured noise in Figure 7. Table 5 summarizes the RMSE of the engine transient health estimation performance by the fusion PF approach with different the variable σ. The covariance of the measurement noise is estimated and real time tunes the quantity  of the FA-PF algorithm as presented previously, and the RMSEs of the FA-PF are almost the same as those of the fusion PF where quantity σ is equal to the As shown in Figure 6, the engine transient performance estimates by the fusion PF approaches seem to deviate from their real values once the noise of sensor P 22 varies, no matter the quantity of measurement noise σ is equal to R 0 , R 1 or R 2 . The fusion adaptive PF has sound tracking performance due to the quantity σ adaptive to the real measurement noise real time in Figure 6d. We can obtain similar results in the case that all sensors for transient performance monitoring have uncertain measured noise in Figure 7. Table 5 summarizes the RMSE of the engine transient health estimation performance by the fusion PF approach with different the variable σ. The covariance of the measurement noise is estimated and real time tunes the quantity σ of the FA-PF algorithm as presented previously, and the RMSEs of the FA-PF are almost the same as those of the fusion PF where quantity σ is equal to the true value of measurement noise covariance in Table 5. Nevertheless, if the quantity σ is not tuned to the covariance of the real measurement noise, the fusion PF will produce a larger RMSE in terms of the engine transient performance estimation.  Table 5. Nevertheless, if the quantity  is not tuned to the covariance of the real measurement noise, the fusion PF will produce a larger RMSE in terms of the engine transient performance estimation.

Engine Health Monitoring Test
Finally a test of the engine health monitoring is carried out to evaluate the proposed method at ground level. The engine input variables, Wf and A8, are fed into the engine as shown in Figure 8. The engine NH representing the engine operation varies as follows: about 0.91 before 2.7 s, increasing from 0.91 to 1.0, then decreasing from 1.0 to 0.91, and about 0.91 to the end. During this process, the engine thrust increases from 0.819 to 1.0, and then back to 0.819. Four abrupt faults depicted in Table 2 are separately injected into the engine at 2 s. Figure 9 depicts the variance estimates of measurement noise, which are applied to tune the quantity  in the FA-PF method. The noise pollutes the true measurement, and it changes along the engine operation condition (the larger the power condition, the more noise enforcement). Table 6 summarizes the performance of engine health monitoring by four PF methods in the abrupt fault cases. The fusion PF with inequality constraints is defined by the FC-PF. The detection results in

Engine Health Monitoring Test
Finally a test of the engine health monitoring is carried out to evaluate the proposed method at ground level. The engine input variables, W f and A 8 , are fed into the engine as shown in Figure 8. The engine NH representing the engine operation varies as follows: about 0.91 before 2.7 s, increasing from 0.91 to 1.0, then decreasing from 1.0 to 0.91, and about 0.91 to the end. During this process, the engine thrust increases from 0.819 to 1.0, and then back to 0.819. Four abrupt faults depicted in Table 2 are separately injected into the engine at 2 s. Figure 9 depicts the variance estimates of measurement noise, which are applied to tune the quantity σ in the FA-PF method. The noise pollutes the true measurement, and it changes along the engine operation condition (the larger the power condition, the more noise enforcement). Table 6 summarizes the performance of engine health monitoring by four PF methods in the abrupt fault cases. The fusion PF with inequality constraints is defined by the FC-PF. The detection results in Table 6 show that the FC-PF and FA-PF have the less estimation errors and provide more stable estimates than the other PF methods because of the prior knowledge used. Furthermore, the performance of the FA-PF for engine health monitoring is the best due to the quantity of the importance function σ adaptive to the measurement noise real-time.  Table 6 show that the FC-PF and FA-PF have the less estimation errors and provide more stable estimates than the other PF methods because of the prior knowledge used. Furthermore, the performance of the FA-PF for engine health monitoring is the best due to the quantity of the importance function  adaptive to the measurement noise real-time.

Conclusions
This paper describes the use of data fusion based on the adaptive PF for gas turbine dynamic performance monitoring. The state estimation in the fusion estimator architecture includes three steps: several local filters working in parallel to obtain individual sensor-based estimates, one master  Table 6 show that the FC-PF and FA-PF have the less estimation errors and provide more stable estimates than the other PF methods because of the prior knowledge used. Furthermore, the performance of the FA-PF for engine health monitoring is the best due to the quantity of the importance function  adaptive to the measurement noise real-time.

Conclusions
This paper describes the use of data fusion based on the adaptive PF for gas turbine dynamic performance monitoring. The state estimation in the fusion estimator architecture includes three steps: several local filters working in parallel to obtain individual sensor-based estimates, one master

Conclusions
This paper describes the use of data fusion based on the adaptive PF for gas turbine dynamic performance monitoring. The state estimation in the fusion estimator architecture includes three steps: several local filters working in parallel to obtain individual sensor-based estimates, one master filter fusing these local estimates to yield a global state estimate, and the global estimates serving as the feedback to each local filter with information-sharing strategy. A systematic comparison of the fusion PF methods for transient performance estimation is presented. Gradual deterioration, abrupt faults and their mixtures are typically considered as the engine anomaly scenarios in the test. The fusion PF architecture has better estimation accuracy and less convergence time than the conventional PF architecture. The implementation of the fusion PF is quite straightforward and involves only basic matrix operations. The convergence time by the fusion PF is similar to that by the basic PF, yet the computational burden of the master filter is reduced, because it is shared by local filters in the fusion PF structure.
Moreover, an adaptive PF algorithm is proposed to sufficiently utilize prior information for the engine transient performance detection with measurement noise uncertainty. The heuristic health knowledge is usually neglected in model-based engine diagnoses due to the complex mathematics application to the conventional PF. In this paper, the engine prior health information represented by the inequality constraints is enforced in the fusion PF algorithm. In addition, the uncertainty of the measurement noise in the engine dynamic operation is considered in the fusion PF. The covariance of the sensed noise is estimated and then applied to tune the importance weights of the PF. The improvements brought by the data fusion based on the adaptive PF have been illustrated on the application of the engine transient performance monitoring with noise uncertainty. The experiments show that the proposed method leads to the more reliable assessments of the engine health conditions no matter whether the cases involve gradual engine deterioration or abrupt faults. With information fusion, a priori knowledge and measurement noise adaption in mind, the data fusion based on the adaptive PF approach seems therefore to be more promising. The present work by the authors has shown the advantages of data fusion based on an adaptive PF for gas turbine health monitoring. It is not considered that the phase errors between different measurements combined with the senor lag will affect the shift in engine performance, and it would be interesting to discuss the proposed approach of this paper with addition of measurement differences.