MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 832851 10.1155/2014/832851 832851 Research Article A GM (1, 1) Markov Chain-Based Aeroengine Performance Degradation Forecast Approach Using Exhaust Gas Temperature Zhao Ning-bo Yang Jia-long Li Shu-ying Sun Yue-wu Xu Qingsong College of Power and Energy Engineering Harbin Engineering University Harbin 150001 China hrbeu.edu.cn 2014 642014 2014 15 12 2013 09 03 2014 16 03 2014 6 4 2014 2014 Copyright © 2014 Ning-bo Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Performance degradation forecast technology for quantitatively assessing degradation states of aeroengine using exhaust gas temperature is an important technology in the aeroengine health management. In this paper, a GM (1, 1) Markov chain-based approach is introduced to forecast exhaust gas temperature by taking the advantages of GM (1, 1) model in time series and the advantages of Markov chain model in dealing with highly nonlinear and stochastic data caused by uncertain factors. In this approach, firstly, the GM (1, 1) model is used to forecast the trend by using limited data samples. Then, Markov chain model is integrated into GM (1, 1) model in order to enhance the forecast performance, which can solve the influence of random fluctuation data on forecasting accuracy and achieving an accurate estimate of the nonlinear forecast. As an example, the historical monitoring data of exhaust gas temperature from CFM56 aeroengine of China Southern is used to verify the forecast performance of the GM (1, 1) Markov chain model. The results show that the GM (1, 1) Markov chain model is able to forecast exhaust gas temperature accurately, which can effectively reflect the random fluctuation characteristics of exhaust gas temperature changes over time.

1. Introduction

Prognostic and health management for aeroengine are the main concerns for many researchers and users in order to provide more useful information for the safe operation [1, 2]. Performance degradation forecast technology for assessing degradation states quantitatively based on gas path performance parameters is one of the most important technologies, which can improve the safety, reliability, and maintenance of aeroengine [3, 4]. Therefore, how to improve the forecast precision by the advanced condition monitoring techniques to extract hidden, unknown, and the useful information from large amounts of monitoring data is emphasized on in the study of performance degradation forecast.

In the past half-century, different methods have been developed to analyze aeroengine gas path performance parameters for performance degradation [5, 6], fault diagnosis , remaining service life , and reliability . Li developed a novel adaptive gas path analysis (adaptive GPA) approach to estimate actual engine performance and gas path component health status by using gas path measurements . Chen established the Artificial Neural Network (ANN) model of aeroengine performance trend forecasting by using the strong nonlinear mapping ability of ANN and the phase space reconstruction theory . Jiang took complex engine gas path system as a grey system and evaluated the samples utilizing the grey relationship degree theory to achieve the gas path fault diagnosis of aeroengine . Wu et al. proposed Support Vector Machines (SVM) to predict the residual life of aeroengine based on the data of the actual gas path parameters monitoring information and failure event report from the aeroengine . Ren and Zuo developed a residual life prediction model based on Bayesian updating methods and particle swarm optimization with immunity algorithms through analyzing the performance degradation process of aeroengine . On the basis that the gas path performance parameters of aeroengine were analyzed, the aeroengine residual life prediction method based on performance deterioration pattern was proposed by Fu et al. . Wang and Jiang provided a performance reliability prediction method based on Support Vector Machines (SVM) for aeroengine by using condition monitoring information .

Several approaches have been introduced to forecast the gas path performance parameters of aeroengine. Based on the research of aeroengine performance parameters relativity, with the condition of small samples and variables with multiple correlations, Shi et al. proposed a partial least-squares regression method to build short time forecasting model of aeroengine performance parameter under the condition of small samples . In order to improve the forecasting accuracy of aeroengine performance parameters, Li et al. decomposed the original sequence by using wavelet transform and forecasted the subsequences in different frequency bands by using Auto Regressive Moving Average (ARMA) or Auto Regressive Integrated Moving Average (ARIMA) . Zhong et al. adopted the process neural networks to predict aeroengine performance parameters . Zhang and Wang adopted Support Vector Machine (SVM) regression approach to monitor an aeroengine health and condition by building monitoring models of main aeroengine performance parameters . In order to predict the change tendency of aeroengine performance parameters effectively, a novel exhaust gas temperature (EGT) prediction method named process support vector machine (PSVM) was proposed by Fu and Zhong . Ilbas and Turkmen dealt with the estimation of exhaust gas temperature (EGT) of a CFM56-7B turbofan engine using Artificial Neural Network (ANN) at two different power settings, maximum continuous and take-off .

From the literature described above, statistical and artificial intelligence based approaches are the two main techniques. Auto Regressive (AR), Moving Average (MA), Auto Regressive Moving Average (ARMA), and Auto Regressive Integrated Moving Average (ARIMA) can be mentioned as statistical models, while Artificial Neural Network (ANN) and Support Vector Machines (SVM) have been most widely used as artificial intelligence approaches. The essences of the above approaches are establishing the appropriate time series model by analyzing historical data. The modeling processes of the statistical based approaches are relatively simple [12, 13]. However, these approaches for accurately forecasting trends depend on the law of the distribution of historical data as well as large amounts of observed data. ANN-based approaches eliminate the limitations of traditional regression methods and accurately establish mapping between input and output variables [19, 20]. It can approximate an arbitrary nonlinear function with satisfactory precision. A great deal of training data and relatively long training period for robust generalization can ensure the success of ANN-based approaches. Although the support vector machines-based approach has many special advantages that can resolve problems, such as the small-sample set, nonlinear, and high dimensions. The kernel functions and penalty factor are difficult to determine, which can also influence the forecast accuracy [17, 18]. Conventional methods of accurate forecast time series trends depend on the sufficiency and completeness of the information obtained.

In practical applications, it is difficult to obtain the complete information because of many reasons. Besides, the aeroengine gas path performance parameters are often highly nonlinear, stochastic, and nonstationary. Therefore, not only the conventional statistical models are not as accurate as the artificial neural network-based approaches for aeroengine gas path performance parameters trend forecast problems, but the traditional methods may also be too complex to be used in forecasting future values of time series.

Grey system theory, proposed by Deng in 1982 , is a mathematical analysis of systems with incomplete information and discrete data. As a simple forecast model, GM (1, 1) model has been widely and successfully applied to various systems such as economic and industrial and so forth . In aeroengine, Wu introduced a wear fault forecast model based on grey system theory for aeroengine . The results showed that compared with the traditional prediction methods of time sequence AR( n ) model, the forecast method of the GM (1, 1) model has the advantages of higher precision of forecast and lower sensitivity to the unequal intervals among the original data sequences for the modeling. However, the forecast accuracy of GM (1, 1) model is unsatisfactory when original data shows great randomness because GM (1, 1) model is only a just order single variable grey model .

The goal of the paper is to introduce the time series forecast based on grey system theory to the forecast modeling of aeroengine gas path performance parameters. Firstly, a type of time series forecast method based on GM (1, 1) model is introduced for aeroengine exhaust gas temperature (EGT). This method can effectively solve the trend forecast problems of EGT under incomplete information and discrete small sample data. Then, Markov chain model is integrated into GM (1, 1) model in order to enhance the forecast performance, which can solve the influence of random fluctuation data on forecast accuracy and achieve an accurate estimate of the nonlinear EGT. A real case of aeroengine EGT from CFM56 aeroengine of China Southern is used to test the capability of the proposed improved model.

The rest of this paper is organized as follows. In Section 2, the degradation signature of aeroengine is analyzed. Section 3 briefly describes the modeling methodologies of aeroengine based on GM (1, 1) model. The hybrid model which combines Markov chain model with GM (1, 1) model is discussed in Section 4. The application and discussion are illustrated and detailed in Section 5. Finally, some conclusions are presented in Section 6.

2. Aeroengine Degradation Signatures and Problem Description

The performance of an aeroengine will deteriorate over the time due to different gas path component degradations such as fouling, erosion, corrosion, and foreign object damage . There are many gas path performance monitoring parameters for civil aeroengine. Due to the fact that performance degradation modes of different parts are different with the increase of aeroengine service time, it is very important to select the appropriate measured parameters that can reflect the performance degradation of aeroengine to realize aeroengine performance degradation forecast.

Outlet temperature of combustor chamber is the most important performance parameter for aeroengine. Not only does it affect the overall performance of the engine, but also it directly determines the ultimate strength of turbine blade. For example, the creep life of hot channel components can reduce the order of magnitude when outlet temperature of combustor chamber increases 50°C , which may cause major fault and incur great maintenance costs. However, outlet temperature of combustor chamber is usually too high to be measured with available instrumentation and temperature distribution is extremely uneven. According to the well-defined Brayton thermodynamic cycle, there is a consistent relationship between outlet temperature of combustor chamber and exhaust gas temperature (EGT). Thus, EGT, as a measured parameter, is often used for engine control, condition monitoring, fault diagnosis, and maintenance decisions. When other conditions remain the same, the higher the EGT is, the more serious the performance degradation of aeroengine is.

Considering the gas path performance monitoring parameters, the multiple linear regression models for the relationship between EGT and other parameters were established by Song et al. . The results showed that there were a strong linear correlation between the performances parameters, such as low turbine outlet pressure, high rotational speed, high pressure compressor outlet temperature, low rotational speed, and high pressure compressor outlet pressure can be reflected through the change of EGT. Yilmaz  found similar results by analyzing the relationship between EGT and engine operational parameters at two different power settings, including maximum continuous and take-off, in the CFM56-7B turbofan engine. Hence, the EGT is often used to evaluate the health states of aeroengine and determine the maintenance policy .

Aeroengine EGT can be divided into take-off EGT and cruise EGT in accordance with different data acquisition stages during flight . Take-off EGT means the exhaust gas temperature in the take-off stage and maximum thrust. The aeroengine can exceed the normal temperature most easily in this stage which can bring the destructive effect to the engine. In the actual process of engineering application, the take-off EGT margin (EGTM) is often used instead of the take-off EGT parameter to performance analysis. EGTM is the deviation between the actual value and EGT red line value when plane takes off with full power at sea level pressure, inflection point temperature conditions. Among them, EGT red line value refers to the allowed maximum EGT that is given by the manufacturers. The calculation formula of EGTM is defined as (1) EGTM ( C ° ) = EG T B - EG T M , where subscript B represents a redline value and subscript M represents actual monitoring value.

Figure 1 shows the sequence that denotes the EGTM procured from six CFM56 aeroengines of China Southern . The sampling interval is about a 200-flight cycle, which can be approximately considered continuous equal interval sampling after data preprocessing based on multiple interpolation method.

EGTM sequence and trend analysis for CFM56 aeroengine .

From Figure 1, it is easy to see that the entire EGTM data sequence has obvious downtrend over time for the six aeroengines. The EGTM can decrease 72°C in the 12400 cycle. However, the degradation process of different engines is not quite the same. For a single engine, the changing process of EGTM sequence is a complex nonlinear process. Most of the differences between adjacent sampling points are less than 10°C and the maximum reached 20°C.

Based on the above analysis, aeroengine performance degradation forecast can be solved as EGTM time series forecasting problem. However, it is difficult to establish a precise mathematical model to describe EGTM that can be affected by many uncertain factors. Therefore, the key problem lies in how to establish precise forecast model under incomplete information and discrete small sample data in order to achieve an accurate estimate of the nonlinear EGTM parameters.

3. EGTM Forecast Modeling Based on GM (1, 1) Model

Based on the temporal variation characteristics of aeroengine EGTM mentioned above, the system of aeroengine performance parameter EGTM can be regarded as a grey dynamic system. This section briefly describes the modeling methodologies about GM (1, 1) model and provides an EGTM trend forecasting framework based on GM (1, 1) model.

3.1. Modeling Methodologies Based on GM (1, 1) Model

In grey systems theory, the most commonly used grey forecast model is GM (1, 1) model, which is successfully employed in time series forecast applications with the uncertain problems under discrete data and incomplete information .

Generally speaking, forecast based on GM (1, 1) model can be regarded as curve fitting analysis in time series . In order to effectively reduce the discreteness of the original time series and reveal the hidden regular pattern in the system development, the first order accumulation generating operation (1-AGO) is used firstly before the first order differential equation is adopted to match the data. Then, whitening equation can be solved by the ordinary least square method to the time response sequences at time n . Finally, the first order inverse accumulating generation operation (1-IAGO) is employed to establish the GM (1, 1) forecast model and obtain the predicted value. The detailed procedure is shown as in Figure 2.

The forecast system based on GM (1, 1) model .

Consider the following nonnegative EGTM time sequence X ( 0 ) : (2) X ( 0 ) = ( x ( 0 ) ( 1 ) , x ( 0 ) ( 2 ) , , x ( 0 ) ( n ) ) , n 4 , where n is the sample size of the data.

To increase the forecast of the GM (1, 1) model, the first order accumulating generation operation (1-AGO) is derived from the original sequence X ( 0 ) . The 1-AGO is derived as (3) x ( 1 ) ( k ) = i = 1 k x ( 0 ) ( i ) , k = 1,2 , , n .

Then the accumulated sequence X ( 1 ) that is obtained using the 1-AGO formation can be expressed as (4) X ( 1 ) = ( x ( 1 ) ( 1 ) , x ( 1 ) ( 2 ) , , x ( 1 ) ( n ) ) , n 4 .

The generated mean sequence Z ( 1 ) of X ( 1 ) is defined as (5) Z ( 1 ) = ( z ( 1 ) ( 1 ) , z ( 1 ) ( 2 ) , , z ( 1 ) ( n ) ) , where z ( 1 ) ( k ) is the mean value of adjacent data; that is, (6) z ( 1 ) ( k ) = 1 2 x ( 1 ) ( k ) + 1 2 x ( 1 ) ( k - 1 ) , k = 2,3 , , n .

The grey model GM (1, 1) can be expressed by one variable, and the grey difference equation is defined as (7) x ( 0 ) ( k ) + a z ( 1 ) ( k ) = b .

And its whitening equation is (8) d x ( 1 ) ( t ) d t + a x ( 1 ) ( t ) = b , where coefficients a and b are called developing and grey input coefficients, respectively. By the ordinary least square method, the coefficients a and b can be obtained as (9) [ a , b ] T = ( B T B ) - 1 B T Y , where (10) Y = [ x ( 0 ) ( 2 ) , x ( 0 ) ( 3 ) , , x ( 0 ) ( n ) ] T , B = [ - z ( 1 ) ( 2 ) 1 - z ( 1 ) ( 3 ) 1 · · · · · · - z ( 1 ) ( n ) 1 ] .

The solutions x ( 1 ) ( t ) of (8) can be obtained using the ordinary least squares method as follows: (11) x ( 1 ) ( t ) = [ x ( 1 ) ( 1 ) - b a ] e - a ( t - 1 ) + b a , where x ( 1 ) ( 1 ) = x ( 0 ) ( 1 ) .

Hence, the time response sequences of (8) at time ( k + 1 ) are (12) x p ( 1 ) ( k + 1 ) = [ x ( 0 ) ( 1 ) - b a ] e - a k + b a .

To obtain the forecast value of the primitive data at time ( k + 1 ) , the first order inverse accumulating generation operation (1-IAGO) is employed to establish the following grey model: (13) x p ( 0 ) ( k + 1 ) = [ x ( 0 ) ( 1 ) - b a ] e - a k ( 1 - e a ) .

And the predicted value of the primitive data at time ( k + H ) is (14) x p ( 0 ) ( k + H ) = [ x ( 0 ) ( 1 ) - b a ] e - a ( k + H - 1 ) ( 1 - e a ) .

Compared with the statistical models, GM (1, 1) model need not find the statistics features of original series. So GM (1, 1) model gets rid of the shadow of large-sample statistics in terms of information availability degree . Besides, only two coefficients are required to be identified in (8), which means that the number of data sample used in GM (1, 1) model is rather small. In other words, GM (1, 1) model can realize the forecast by using only sample data sequence and is often used as a short term forecast. This is the biggest advantage over the artificial intelligence method.

3.2. EGTM Forecast Framework Based on GM (1, 1) Model

According to the above method, a GM (1, 1) model based approach for aeroengine performance degradation forecast using EGTM signatures is illustrated in Figure 3. The following shows the details of the forecast process.

Flowchart of EGTM forecast based on GM (1, 1) model approach.

Step 1.

Generated sufficient EGTM samples from the historical database and the essential preprocessing upon EGTM data are carried out before data analysis, such as supplementary data, eliminating noise and outliers. After that, the samples can be divided into training samples and testing samples.

Step 2.

The EGTM forecast model based on GM (1, 1) method is established by using the training samples.

Step 3.

The testing samples are used to verify the forecast performance of the GM (1, 1) model. Step 2 will be repeated if the forecast model accuracies do not meet the requirements.

Step 4.

Apply the GM (1, 1) model that meets the accurate requirement to EGTM measured signatures obtained from real aeroengine to forecast.

4. EGTM Forecast Modeling Based on GM (1, 1)-Markov Chain Model

As a first order single variable grey model, GM (1, 1) model provides an excellent approach to forecast uncertain systems . However, the forecast accuracy of GM (1, 1) model for EGTM series with large random fluctuations is lower, which cannot satisfy the engineering requirement. In order to enhance the forecast performance, Markov chain model is integrated into GM (1, 1) model to extract the random fluctuation of experimental data and solve the influence of random fluctuation data on forecast accuracy. The improved model is defined as GM (1, 1) Markov chain model.

A GM (1, 1) Markov chain based approach for aeroengine performance degradation forecast using EGTM signatures is illustrated in Figure 4. To achieve the aim of this study, the forecast has two stages. The original data are modeled by the GM (1, 1) model firstly. Then the residual errors between the fitting values and the actual values for all previous time steps can be obtained. After that, the transition behavior of those residual errors by Markov transition matrices is established and the possible correction for the forecast value can be made from those Markov matrices. The following shows the details of the forecast process.

Flowchart of EGTM forecast based on GM (1, 1) Markov chain approach.

4.1. Residual Errors

According to the forecast values x p ( 0 ) ( k ) obtained from GM (1, 1) model by (13) and the real values x ( 0 ) ( k ) , the residual errors series can be obtained as (15) ε ( k ) = x ( 0 ) ( k ) - x p ( 0 ) ( k ) , where k respects the time step.

4.2. Division State

The real values of x ( 0 ) ( k ) are distributed in the region of the forecast value x p ( 0 ) ( k ) which may be divided into a convenient number of contiguous intervals. When x ( 0 ) ( k ) falls in interval i , one of n such intervals, it may be regarded as corresponding to a state Q i which can be denoted as follows: (16) Q i = [ Q 1 i , Q 2 i ] , i = 1,2 , , n , Q 1 i = x p ( 0 ) ( k ) + A i , Q 2 i = x p ( 0 ) ( k ) + B i , where Q 1 i and Q 2 i are the lower boundary and upper boundary of the Q i state. A i and B i are the deviation between Q 1 i and Q 2 i with forecast value x p ( 0 ) ( k ) , respectively. Hence, Q 1 i and Q 2 i can reflect the dynamic characteristics of the error residual series.

4.3. Transition Probability and Matrix

Let the state space of a Markov chain { x m } be I , the current state be i , and the next state be j ; then, the transition probability is written as (17) P i j = P { x m + 1 = j | x m = i | } , ( i , j ) I , where the P i j is independent of m .

The matrix P , formed by placing P i j in row i and column j , for all i and j , is called the transition probability matrix or chain matrix. Note that the elements of the matrix P satisfy the following two properties: (18) P i j 0 , i , j I , j I P i j = 1 , i I .

The transition probability of state is written as (19) P i j ( r ) = M i j ( r ) M i , i , j = 1,2 , , n , where P i j ( r ) is the probability of transition from state i to j by r steps. M i j ( r ) is the transition times from state i to j by r steps and M i is the number of data belonging to the i th state. Because the transition for the last r entries of the series is indefinable, M i should be counted by the first as l - r entries. l is the quantity of entries of the original series. Then, the transition probability matrix of state can be written as (20) P ( r ) = [ P 11 ( r ) P 12 ( r ) P 1 n ( r ) P 21 ( r ) P 22 ( r ) P 2 n ( r ) P n 1 ( r ) P n 2 ( r ) P n n ( r ) ] .

The transition probability matrix of states P ( r ) reflects the transition rules of the system. The transition probability of states P i j ( r ) reflects the probability of transition from initial state i to probable state j by r steps. It is the foundation of forecast by the Markov probability matrix.

Generally speaking, consider r = 1 and the maximum transition step is 1. Then, P ( 1 ) can be obtained. If the forecast original data is located in the Q i state, the state of next step is determined by the i th row vector of transition probability states P i j ( 1 ) . If max P i j ( 1 ) = P i 1 ( 1 ) , the state of next step is Q 1 . When the state of next step cannot be determined by 1-step transition probability, the 2-step transition probability will be selected.

4.4. Obtaining the Forecast Value

When the possibility of a certain state of the next step is determined by the probabilities in k - 1 row vectors, Q 1 i and Q 2 i can also be obtained. The median in [ Q 1 i , Q 2 i ] is selected as the forecast value, so forecast value of original data sequence is obtained according to the above explanation. (21) x ^ p ( 0 ) ( k ) = x p ( 0 ) ( k ) + ( A i + B i 2 ) .

The main assumption in a Markov chain model is that knowledge of the current state occupied by the process can be sufficient to describe the future probabilistic behavior of the process. Another unique property of this Markov chain model is the existence of a steady state matrix.

5. Experimental Study

In this section, the forecast approach based on GM (1, 1) Markov chain introduced in this paper will be applied to EGTM forecast of CFM56 aeroengine to demonstrate the potential capability of the new approach. The comparisons between the EGTM forecast capabilities using the GM (1, 1) model, GM (1, 1) Markov chain model, and other traditional methods are adopted. Four performance measures are used to examine the forecast accuracy of forecast models in this paper. The relative percentage error (RPE), mean square error (MSE), absolute mean error (AME), and absolute mean percentage error (AMPE) are calculated using the following functions, respectively: (22) RPE = | ε ( k ) | x ( 0 ) ( k ) × 100 % , MSE = 1 n k = 1 n ε ( k ) 2 , AME = 1 n k = 1 n | ε ( k ) | , AMPE = 1 n k = 1 n | ε ( k ) | x ( 0 ) ( k ) × 100 % .

5.1. Sample Data

For this study, the investigators gathered data samples from CFM56 aeroengine of China Southern, which have been described in Figure 1. In order to highlight the forecast performance of the forecast methods under incomplete information and discrete small sample data, only 45 samples are taken to construct a time series { EGT M m }    ( m = 1,2 , , 45 ) from aeroengine 2. Figure 5 shows the distribution of training and testing sample data. As the training sample, the first 20 sample data are first used to establish forecast model. Then the next 25 sample data are used to test the effectiveness of forecast model.

Sample sets partition for EGTM forecast.

5.2. Forecast Results Based on GM (1, 1) Model

After the generation of samples, GM (1, 1) model is established by (13) to forecast the trend from the generated training samples. Specific modeling steps and methods have been given in Figure 3.

The forecast results of EGTM based on GM (1, 1) model and the original series are plotted in Figure 6. From Figure 6, it is clear that the trend of EGTM can be forecasted based on GM (1, 1) model. However, the forecast accuracies of GM (1, 1) model are unstable at some points. And because the accumulated sequence obtained using the 1-AGO formation is monotonically increasing, which is seen in Figure 7, GM (1, 1) model cannot extract random fluctuations of EGTM.

Forecast values of EGTM based on GM (1, 1) model.

The accumulated data of EGTM using the 1-AGO formation.

In order to compare other models, the linear regression model y = - 0.0064 x + 99.968 , nonlinear regression model y = 100.19 exp ( - 0.00007 x ) , ARIMA (1, 0, 0) model, and Radial Basis Function (RBF) Neural Networks are also employed to make the same forecasts. The comparative results are listed in Table 1. We can see that the forecast performance based on nonlinear model is closely approximate to the forecast performance based on GM (1, 1) model. Besides, GM (1, 1) model is better than linear model or RBF model that often requires a large amount of training samples to get a higher forecast precision. Furthermore, the differences between GM (1, 1) model and ARIMA (1, 0, 0) model are that GM (1, 1) model does not have a requirement to the distribution characteristics of training samples while ARIMA (1, 0, 0) model not only requires a large amount of training samples, but also requires the law for the distribution of training samples. Considering the above factors, GM (1, 1) model is more suitable for EGTM forecast, especially for small sample data.

The comparison results of GM (1, 1) model and GM (1, 1) Markov chain model.

Cycle Forecast model
Original values Linear model Nonlinear model ARIMA model RBF model GM (1, 1) model
4200 86 73.088 74.66925168 76.24820454 83.893778 74.850179
4400 86 71.808 73.63116571 85.04659158 75.345291 73.899849
4600 76 70.528 72.60751169 85.09212012 66.440678 72.961583
4800 72 69.248 71.59808897 79.58775994 77.213138 72.035231
5000 64 67.968 70.60269971 69.85924741 75.907728 71.12064
5200 66 66.688 69.6211488 69.85924741 51.433499 70.217661
5400 72 65.408 68.65324386 70.37022346 77.51105 69.326147
5600 70 64.128 67.69879518 74.57638112 68.512595 68.445951
5800 58 62.848 66.75761567 72.84613016 52.139635 67.576931
6000 62 61.568 65.82952088 62.6255194 75.584126 66.718945
6200 66 60.288 64.91432888 65.48119394 63.496734 65.871852
6400 58 59.008 64.01186029 68.51302315 64.808649 65.035514
6600 60 57.728 63.12193824 61.40259349 54.990815 64.209794
6800 60 56.448 62.24438829 62.78391435 60.518705 63.394558
7000 60 55.168 61.37903843 62.49368451 62.428861 62.589673
7200 64 53.888 60.52571907 62.25825584 61.24122 61.795007
7400 60 52.608 59.68426294 65.60171099 57.527047 61.01043
7600 68 51.328 58.85450511 62.02015608 58.92226 60.235815
7800 62 50.048 58.03628296 68.8277359 53.706033 59.471034
8000 68 48.768 57.22943611 63.72946594 60.035962 58.715964
8200 66.5 47.488 56.4338064 68.75093548 60.237787 57.97048
8400 65 46.208 55.64923791 67.47264457 67.031388 57.234461
8600 55 44.928 54.87557684 66.17459205 53.436919 56.507787
8800 55 43.648 54.11267156 57.5174732 55.602918 55.79034
9000 57 42.368 53.36037254 57.31506291 58.04245 55.082001
5.3. Forecast Results Based on GM (1, 1) Markov Chain Model 5.3.1. Division State

In the previous literature, there is no unified standard to determine number of state and the boundary of the Q i state, and state division is performed on experience, which is called the hard division approach, which might result in the weakness of forecast accuracy and algorithm application, due to the fact that the different divisions would be made by different persons. Generally speaking, the value ranges of n are 3–6, and they ensure that each state interval has data.

According to the absolute error series ε ( k ) , the states are partitioned by establishing five contiguous intervals for Markov-chain forecast model. Figure 8 shows the results of state division. The blue curve shows the obtained trend curve by GM (1, 1) model and the red curves show the state lines.

State division for the time series of EGTM.

5.3.2. Transition Probability and Matrix

From Figure 8, it can be seen that the number of samples in each state is as follows: M 1 ( 1 ) = 2 , M 2 ( 1 ) = 1 , M 3 ( 1 ) = 7 , M 4 ( 1 ) = 7 , M 2 ( 1 ) = 2 . Then one-step transition probability for every state can be calculated by using (19) and (20). Consider the following: (23) P ( 1 ) = [ 0 0 1 2 1 2 0 0 0 0 1 0 2 7 1 7 0 2 7 2 7 0 1 7 3 7 3 7 0 0 0 1 0 0 ] .

5.3.3. Forecast Values and Analysis

When the possibility of a certain state of the next step is determined according to 1-step probability transition matrices P ( 1 ) , Q 1 i and Q 2 i can be obtained. The forecast values can be obtained through the same computations. The forecast values by GM (1, 1) model and GM (1, 1) Markov chain model are shown in Figure 9. What is not worthy is that the dynamic change characteristics of EGTM are reflected by calculating the state transfer probability. Figures 10, 11, and 12 show the RPE comparative analysis of different forecast models. From the values of RPE comparison, we find that the forecast accuracies of GM (1, 1) Markov chain model are significantly higher than other models.

: The forecast values by GM (1, 1) model and GM (1, 1) Markov chain model for EGTM.

The forecast RPE by linear regression model and nonlinear regression model for EGTM.

The forecast values by ARMIA (1, 0, 0) model and RBF model for EGTM.

The forecast values by GM (1, 1) model and GM (1, 1) Markov chain model for EGTM.

The comparison of forecast accuracy of the different models is listed in Table 2. From Table 2, it is easy to see that the forecast performance of GM (1, 1) Markov chain model outperforms other models. When all of the factors are considered together, the modeling processes of GM (1, 1) Markov chain model are much easier than ARIMA (1, 0, 0) model and RBF model. Through the test values, we find that GM (1, 1) Markov chain model gives more satisfactory performances in MSE, AME, and AMPE than other models. The AMPE of GM (1, 1) Markov chain model is only 2.727393%, which can meet the demand of engineering application. The results validate the effectiveness of GM (1, 1) Markov chain model.

The comparative analysis results of different models.

Models Evaluating standards
MSE AME AMPE/%
Linear model 108.1934 8.572960000 13.07414
Nonlinear model 37.67832 4.861358116 7.209246
ARIMA model 37.36459 4.785827391 7.447413
RBF model 48.05668 5.59155023 8.443975
GM (1, 1) model 34.2345 4.673887228 6.982612
GM (1, 1) Markov chain model 11.75861 1.981937814 2.727393
6. Conclusions

According to the characteristics of the aeroengine gas path performance parameters, EGTM is used to realize aeroengine performance degradation forecast. Based on the change law of aeroengine EGTM, EGTM forecast is solved as a grey system forecast problem. However, it is shown that GM (1, 1) model is only able to accurately forecast the trend of EGTM and the forecast accuracy of GM (1, 1) model is not satisfactory when the EGTM data show great randomness. In order to enhance the forecast performance of GM (1, 1) model, Markov chain model is integrated into GM (1, 1) model. The comparison results show that the forecast accuracy of the improved model named GM (1, 1) Markov chain model is better than other models, especially for the small samples. GM (1, 1) Markov chain model can solve the influence of random fluctuation data on forecast accuracy and achieve an accurate estimate of the nonlinear EGTM.

Conflict of Interests

The authors declared that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The present work is supported by the Fundamental Research Funds for the Central Universities of China (no. HEUCFZ1005).

Bytington C. S. Roemer M. J. Galie T. Prognostic enhancements to diagnostic systems for improved condition-based maintenance 6 Proceedings of the IEEE Aerospace Conference 2002 2815 2824 10.1109/AERO.2002.1036120 Marinai L. Probert D. Singh R. Prospects for aero gas-turbine diagnostics: a review Applied Energy 2004 79 1 109 126 2-s2.0-2642561086 10.1016/j.apenergy.2003.10.005 Litt J. S. Simon D. L. Garg S. A survey of intelligent control and health management technologies for aircraft propulsion systems 2005 NASA-TM-2005-213622 Melcher K. J. Maul W. A. Garg S. Propulsion health management system development for affordable and reliable operation of space exploration systems AIAA 2007 2007-6237 Li Y. G. Gas turbine performance and health status estimation using adaptive gas path analysis Journal of Engineering for Gas Turbines and Power 2010 132 4 9 2-s2.0-78149394220 10.1115/1.3159378 041701 Chen G. Forecasting engine performance trend by using structure self-adaptive neural network Acta Aeronautica et Astronautica Sinica 2007 28 3 535 539 2-s2.0-34347273726 Jiang L. P. Gas path fault diagnosis system of aero-engine based on grey relationship degree Procedia Engineering 2011 15 4774 4779 2-s2.0-84055213144 10.1016/j.proeng.2011.08.893 Wu X. H. Wen D. Li R. G. Chen Z. Z. Huang H. Z. Lv Z. Q. A model to predict the residual life of aero-engine based upon Support Vector Machine Proceedings of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering (QR2MSE '13) July 2013 Chengdu, China 1880 1882 10.1109/QR2MSE.2013.6625946 Ren S. H. Zuo H. F. A combined prediction method for the residual-life of civil aviation engines based on performance degradation Mechanical Science and Technology For Aerospace Engineering 2011 30 1 23 29 Fu X. Y. Zhong S. S. Zhang M. N. An improved method for aeroengine residual life prediction Journal of Harbin Institute of Technology 2013 45 5 51 55 Wang H. X. Jiang Y. Performance reliability prediction of complex system based on the condition monitoring information Mathematical Problems in Engineering 2013 2013 7 10.1155/2013/836517 836517 Shi Y. S. Li Y. Z. Song Y. X. Aeroengine performance parameter forecast based on partial least-squares regression method under condition of small samples Aircraft Design 2012 32 6 47 51 Li X. B. Cui X. L. Lang R. L. Forecasting method for aeroengine performance parameters Journal of Beijing University of Aeronautics and Astronautics 2008 34 3 253 256 2-s2.0-43249083781 Zhong S. S. Lei D. Ding G. Convolution sum discrete process neural network and its application in aeroengine exhausted gas temperature prediction Acta Aeronautica et Astronautica Sinica 2012 33 3 438 445 2-s2.0-84860296710 Zhong S. S. Luan S. G. Ding G. Aeroengine condition monitoring based on process neural networks and aerothermo-parameters Journal of Nanjing University of Science and Technology 2006 30 5 533 536 2-s2.0-33845544944 Ding G. Xu M. Q. Hou L. G. Prediction of aeroengine exhaust gas temperature using process neural network Journal of Aerospace Power 2009 24 5 1035 1039 2-s2.0-66349117138 Zhang C. X. Wang N. Aero-engine condition monitoring based on support vector machine Physics Procedia 2012 24, part B 1546 1552 10.1016/j.phpro.2012.02.228 Fu X. Y. Zhong S. S. Aeroengine turbine exhaust gas temperature prediction using process support vector machines Advances in Neural Networks 2013 7951 Berlin, Germany Springer 300 310 10.1007/978-3-642-39065-4_37 Ilbas M. Turkmen M. Estimation of exhaust gas temperature using artificial neural network in turbofan engines Journal of Thermal Science and Technology 2012 32 2 11 18 Flores J. J. Graff M. Rodriguez H. Evolutive design of ARMA and ANN models for time series forecasting Renewable Energy 2012 44 225 230 2-s2.0-84858289782 10.1016/j.renene.2012.01.084 Deng J. L. Control problems of grey systems Systems & Control Letters 1982 1 5 288 294 2-s2.0-50849151631 Kayacan E. Ulutas B. Kaynak O. Grey system theory-based models in time series prediction Expert Systems with Applications 2010 37 2 1784 1789 2-s2.0-71749116803 10.1016/j.eswa.2009.07.064 Wu Z. F. Guo L. Zuo H. F. The wear fault prediction model of aero-engine based on the gray system theory Proceedings of the IEEE International Conference on Grey Systems and Intelligent Services (GSIS '07) November 2007 Nanjing, China 528 532 2-s2.0-50249154901 10.1109/GSIS.2007.4443330 Li G. D. Masuda S. Nagai M. The prediction model for electrical power system using an improved hybrid optimization model International Journal of Electrical Power & Energy Systems 2013 44 1 981 987 10.1016/j.ijepes.2012.08.047 Wang L. Li Y. G. Abdul Ghafir M. F. Swingler A. A Rough Set-based gas turbine fault classification approach using enhanced fault signatures Proceedings of the Institution of Mechanical Engineers A 2011 225 8 1052 1065 2-s2.0-84856265734 10.1177/0957650911416913 Hamilton P. E. P. Ha D. Exhaust gas temperature capabilities now in system 1 software Product Update 2005 25 1 88 89 Song Y. X. Zhang K. X. Shi Y. S. Research on aeroengine performance parameters forecast based on multiple linear regression forecasting method Journal of Aerospace Power 2009 24 2 427 431 2-s2.0-62549145803 Yilmaz I. Evaluation of the relationship between exhaust gas temperature and operational parameters in CFM56-7B engines Proceedings of the Institution of Mechanical Engineers G 2009 223 4 433 440 2-s2.0-69249092923 10.1243/09544100JAERO474 He E. M. Song L. T. Analysis of EGT and measures to increase the EGT margin Aviation Engineering & Maintenance 1999 6 20 21 Bengtsson M. Condition based maintenance system technology-where is development heading? Proceedings of the 17th European Maintenance Congress May 2004 Barcelona, Spain Zhao Y. Z. Research on methods of life on wing prediction for civil aeroengine [M.S. thesis] 2010 Nanjing, China Nanjing University of Aeronautics and Astronautics Deng J. L. Introduction to grey system theory The Journal of Grey System 1989 1 1 1 24 MR1025897 ZBL0701.90057 He Y. Bao Y. D. Grey-Markov forecasting model and its application System Engineering: Theories & Practice 1992 9 4 59 63 Huang M. He Y. Cen H. Predictive analysis on electric-power supply and demand in China Renewable Energy 2007 32 7 1165 1174 2-s2.0-37849186439 10.1016/j.renene.2006.04.005 Kumar U. Jain V. K. Time series models (Grey-Markov, Grey Model with rolling mechanism and singular spectrum analysis) to forecast energy consumption in India Energy 2010 35 4 1709 1716 2-s2.0-77949266106 10.1016/j.energy.2009.12.021 Li G. D. Yamaguchi D. Nagai M. A GM(1,1)-Markov chain combined model with an application to predict the number of Chinese international airlines Technological Forecasting and Social Change 2007 74 8 1465 1481 2-s2.0-34548301313 10.1016/j.techfore.2006.07.010 Hsu Y. T. Liu M. C. Yeh J. Hung H. F. Forecasting the turning time of stock market based on Markov-Fourier grey model Expert Systems with Applications 2009 36 4 8597 8603 2-s2.0-60249093964 10.1016/j.eswa.2008.10.075