# A Fault Diagnosis Approach for Gas Turbine Exhaust Gas Temperature Based on Fuzzy C-Means Clustering and Support Vector Machine

Authors: Zhi-tao Wang; Ning-bo Zhao; Wei-ying Wang; Rui Tang; Shu-ying Li
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

As an important gas path performance parameter of gas turbine, exhaust gas temperature (EGT) can represent the thermal health condition of gas turbine. In order to monitor and diagnose the EGT effectively, a fusion approach based on fuzzy C-means (FCM) clustering algorithm and support vector machine (SVM) classification model is proposed in this paper. Considering the distribution characteristics of gas turbine EGT, FCM clustering algorithm is used to realize clustering analysis and obtain the state pattern, on the basis of which the preclassification of EGT is completed. Then, SVM multiclassification model is designed to carry out the state pattern recognition and fault diagnosis. As an example, the historical monitoring data of EGT from an industrial gas turbine is analyzed and used to verify the performance of the fusion fault diagnosis approach presented in this paper. The results show that this approach can make full use of the unsupervised feature extraction ability of FCM clustering algorithm and the sample classification generalization properties of SVM multiclassification model, which offers an effective way to realize the online condition recognition and fault diagnosis of gas turbine EGT.


---

## Body

## 1. Introduction

With the development of high efficiency and clean energy, gas turbine plays an increasingly significant role in different domains, such as aviation and marine propulsion systems, electric power stations, and natural gas transportation petroleum [1]. With the increasing demand of security operation for gas turbine, the traditional regular maintenance technology has been unable to fully keep up with the actual demand and gas turbine health management technology has gradually become one of the most problems concerned by researchers and users in recent years [2]. In order to guarantee the gas turbine to run efficiently under the safe reliable condition, many sensors are often used to monitor the health state of gas turbine in the practical application. Massive amounts of data gathered by these sensors are easy to make difficulties in data analysis and affect the maintenance decision. Therefore, the choices of appropriate monitor parameters, signal processing methods and data mining techniques are very important to realize the health management of gas turbine.Exhaust gas temperature (EGT) is an important gas path performance parameter of gas turbine, which can represent the thermal health condition of gas turbine [3, 4]. Considering the characteristics of different gas path performance monitoring parameters, the multiple linear regression models for analyzing the relationship between EGT and other parameters were established by Song et al. [5]. Their results showed that there were strong linear correlations between different gas path performance parameters and all the low turbine outlet pressure, high rotational speed, high pressure compressor outlet temperature, low rotational speed, and high pressure compressor outlet pressure could be reflected through the change of EGT. Yilmaz [6] also found the similar results by analyzing the relationship between EGT and other engine operational parameters at two different power settings, including maximum continuous and take-off, in the CFM56-7B turbofan engine. Hence, EGT is often used as an important parameter to evaluate the health state of gas turbine and determine the maintenance policy [7].In the past half century, different methods have been developed to monitor and diagnose the EGT of gas turbine. Wang and Yang [8] analyzed many faults of PG6551B industrial gas turbine, such as turbine ablation, combustion component, and fuel system failure. They found that the uniformity of EGT could effectively reflect the feature of above fault. Chen et al. [9] proposed a general regression neural network (GRNN) approach to construct an autodetection network for EGT sensors, on the basis that they also studied the optimizing design of network and error controlling and developed the method of threshold for sensor detection. Based on the advantage of artificial neural networks (ANN), Muthuraman et al. [10] developed an autoassociative neural network approach to detect combustor-related damage by monitoring EGT. Błachnio and Pawlak [11] established a nonlinear observer and chose EGT as the important parameter to evaluate the health state of turbine blades. Korczewski [12–14] analyzed the change rules of EGT for a naval gas turbine engine under steady and unsteady operation conduction in detail. And they proposed an effective approach for detecting and evaluating the failures of the flow section and supply system of gas turbine by using EGT. Kenyon et al. [15] developed an intelligent system for detection of EGT anomalies in gas turbines by using the strong nonlinear mapping ability of ANN. Considering the characteristics of gas turbine operation control based on thermocouple measured exhaust temperatures, Xia et al. [16] discussed the application of Fiber-Bragg-grating-based sensing technology in the EGT measuring of gas turbine. Their results demonstrated that the fiber sensing method was more valuable for the monitoring and fault diagnosis of gas turbine because it could well reflect the changing of EGT. In order to increase the operational availability of industrial gas turbines, Yang et al. [17] presented a generalization of multidimensional linear regression to facilitate multisensor fault detection and signal reconstruction through the use of analytical optimization. Gülen et al. [18] discussed the relationship between EGT and other gas turbine performance parameters, on the basis of which an important diagnostic parameter named profile factor that was the ratio of the maximum exhaust thermocouple and the average of all exhaust temperature thermocouples was used to evaluate the performance of combustor and the whole gas turbine in their paper.From the reviews discussed above, it is noted that the average EGT is often used to evaluate the health state of gas turbine in most researches. However, it may be more important and valuable to extract the relationship among different EGT sensors in order to realize condition monitor and fault diagnosis of gas turbine effectively. Although many studies have presented the effects of EGT distribution characteristics on health state of gas turbine or its hot sections and many analyses were discussed in detail, there was still a lack of systematic research in the area of online automatic identification and fault diagnosis for gas turbine EGT. Besides, EGT can be affected by many uncertain factors in the practical applications, which make it difficult to realize fault diagnosis quickly by using traditional model-driven approach. Therefore, it is very useful to develop a data-driven approach based on artificial intelligence technology in order to improve fault diagnosis accuracy.The fault diagnosis of gas turbine EGT based on data-driven approach essentially is the cluster and classification of fault information. In the concrete implementation process, the training samples including normal and fault information need to be obtained firstly. Then the fault diagnosis model based on artificial intelligence algorithm can be established and trained by using training samples. The feature information collected from sensors will be inputted to the well trained fault diagnosis model and we can get the diagnostic results finally. Obviously, the establishment of fault information features space and design of artificial intelligence algorithm are two key steps to realize accurate fault diagnosis of gas turbine EGT. As mentioned above, many specific state patterns or fault types of gas turbine EGT cannot be determined directly based on experience in the practical applications due to the effects of many uncertain factors. Besides, enough prior knowledge including specific fault types is indispensable for the supervised artificial intelligence algorithms (such as ANN [19] and SVM [20]).Considering the distribution characteristics of gas turbine EGT and the deficiencies of present literatures, a fusion approach based on FCM clustering algorithm and SVM classification model (FCM-SVM) is proposed in this paper. Firstly, FCM clustering algorithm is used to realize clustering analysis and obtain the state patterns of EGT, which means that the preclassification of EGT is completed. Then, SVM multiclassification model is designed to carry out the state pattern recognition and fault diagnosis of EGT. As an example, the historical monitoring data of EGT from an industrial gas turbine is analyzed to verify the effectiveness of the FCM-SVM approach finally.The rest of this paper is organized as follows. In Section2, the distribution characteristics of gas turbine EGT are described briefly. Section 3 introduces the basic theory of FCM clustering algorithm and SVM classification model in detail. The fusion fault diagnosis approach which combines FCM clustering with SVM is discussed in Section 4. Application examples and discussion are included in Section 5. Finally, Section 6 presents some conclusions.


## 2. Signal Feature of Gas Turbine EGT

As mentioned above, it is very important to choose the appropriate measured parameters to monitor and diagnose the health state of gas turbine. Gas path and vibration parameters are two main types in the practical applications [2, 21]. Theoretically, as the most important gas path performance parameter for gas turbine, the outlet temperature of combustor chamber can not only affect the overall performance of engine, but also directly determine the ultimate strength of turbine blade. For example, the creep life of hot channel components can reduce the order of magnitude when the outlet temperature of combustor chamber increases 50°C [4], which may cause major fault and incur great maintenance costs. However, the outlet temperature of combustor chamber is usually so high that it cannot be measured directly by using conventional sensors. According to the well-defined Brayton thermodynamic cycle, there is a consistent relationship between the outlet temperature of combustor chamber and EGT. Therefore, EGT, as a measured parameter, is often used for gas engine control, condition monitoring, fault diagnosis, and maintenance decisions.Compared with the average EGT, EGT profile can contain more information about the health state of gas turbine. Figures1 and 2 show two EGT profiles with the same average EGT of an industrial gas turbine that has 12 EGT sensors. From Figures 1 and 2, it is easy to see that the EGT profiles of normal and fault condition are different although the average EGT are the same. This means that the fault information is incorrect or incomplete if only the average EGT is used to monitor the health state of gas turbine. Besides, Figures 1 and 2 also show that all the sensors should give similar outputs when gas turbine operates in normal condition. If the component of gas turbine is failure, different temperatures will be observed. Therefore, the uniformity of EGT can more effectively reflect the health state of gas turbine, especially for steady state condition.Figure 1
EGT profile of gas turbine with normal condition.Figure 2
EGT profile of gas turbine with fault condition.In order to quantitatively describe the uniformity of EGT, Mao [22] presented three indexes which can be calculated by the following functions. Assume that there are n sensors and their original outputs are Ti, i=1,…,n, respectively:(1)H1=T1′-Tn′,H2=T1′-Tn-1′,H3=T1′-Tn-2′,where Ti′ is the transform value of Ti and T1′>T2′>⋯>Tn-2′>Tn-1′>Tn′.According to above indexes, it is obvious that all the values ofH1, H2, and H3 are smaller when gas turbine operates in normal condition. If a sensor fails, it usually causes H1 or three indexes to increase. But the values of H1 and H2 or all three indexes can increase due to hot sections failure. However, it is worth pointing out that the above indexes only consider 4 EGT and others are ignored. Besides, it is difficult to diagnose the specific cause of failure when the sensor is fault because the above indexes ignore the adjacent information between different sensors. Therefore, there are some limitations to evaluate the uniformity of EGT only using the above three indexes. In order to solve this problem, all the measured EGT are used to realize cluster and fault diagnosis in this paper.


## 3. Basic Theory

### 3.1. Fuzzy C-Means Clustering Algorithm

As an unsupervised machine learning method, FCM clustering algorithm was improved by Bezdek [23] in 1981 in order to solve the hard clustering problem by using fuzzy set theory. In the FCM clustering algorithm, membership degree function is used to indicate the extent to which each data point belongs to each cluster, and this information is also used to update the values of cluster centers [24]. Based on the concept of fuzzy C-partition, FCM clustering algorithm has been applied successfully in a wide variety of applications, such as image segmentation [25], data mining [26], thermal system monitoring [27], and fault diagnosis [28].For the sample setX=x1,x2,…,xn, the object of FCM clustering algorithm is to divide the sample set into c groups and obtain the cluster centers by minimizing the following dissimilarity function [29]:(2)minJmU,V=∑i=1c∑k=1nuikmdik2,where m is the fuzzy weighting parameter varying in the range [1,∞]. The bigger the m, the more fuzzy the final cluster result. U is fuzzy partition matrix, V is cluster center matrix, and n and c are the number of samples and cluster centers, respectively. uik is the fuzzy membership degree of the kth sample in the ith cluster and it should be meeting the following three constraints [30]:(3)uik∈0,1,1≤i≤c,1≤k≤n,∑i=1cuik=1,1≤k≤n,∑k=1nuik∈0,n,1≤i≤c.For the distancedik between kth sample xk and the centre of ith cluster Vi, it can be calculated by using Euclidean distance as follows:(4)dik=xk-Vi=∑j=1mxkj-Vij2,1≤i≤c,1≤k≤n.V
i
j can be calculated by utilizing the following formulation:(5)Vij=∑k=1nuikmxkj∑k=1nuikm,1≤i≤c,1≤j≤m.In essence, fuzzy cluster is performed through an iterative optimization by updating fuzzy membership degree [29]:(6)uiks+1=1∑j=1cdiks/djks2/m-1,1≤i≤c,1≤k≤n,where s is the iterative step.When the below requirement is met, we can stop iteration and obtain the cluster result:(7)Us+1-Us≤ε,where ε is the iterative threshold in the range [0,1].Based on above method, the cluster process of FCM clustering algorithm is virtually to determine the fuzzy membership degree and cluster centers through continuous iteration, which is shown in Figure3.Figure 3
Calculation flowchart of fuzzy C-means clustering algorithm.


### 3.2. Support Vector Machine Classification Model

Compared with the conventional classifiers, support vector machine (SVM), developed by Vapnik [31], can effectively solve the classification problem by implementing the structure risk minimization based on statistical learning theory. Nowadays, SVM has been widely and successfully applied to detection and diagnosis of machine conditions due to its high accuracy and good generalization for a smaller number of samples [32, 33].SVM is initially used to deal with binary classification problems. Its core idea is to transform the sample data from original space to a higher dimensional feature space through some nonlinear mapping functions and then find the optimal separating hyperplane in this feature space to realize linear classification. Figure4 shows the classification principle based on SVM for the nonlinear classification problem.Figure 4
Classification of two classes using SVM.For the nonlinear training sample data set including two classesxi,yi, i=1,2,…,n, x∈Rn, y∈-1,+1, n is the number of samples. The nonlinear mapping function φ(x) will be used to transform the sample data from original space to a higher dimensional feature space and the optimal separating hyperplane can be constructed to maximize the margin between the two classes by the following linear function:(8)fx=ω·φx+b,where ω is the normal vector of optimal separating hyperplane and b is a scalar.In essence, the solution of optimal separating hyperplane is the corresponding constraint optimization problem:(9)min12ω2+C∑i=1nζis.t.yiω·φxi+b+ζi≥1,1≤i≤nζi≥0,1≤i≤n,where C is the penalty factor that can realize the trade-off between empirical risk and confidence interval. ζi is slack factor.Combining the method of Lagrange multipliers, the above convex optimization problem can be simplified into the dual quadratic optimization problem:(10)maxLa=∑i=1nai-12∑i,j=1naiajyiyjφxi·φxjs.t.ai≥0,1≤i≤n∑i=1naiyi=0,where a is Lagrangian multiplier.Then, the nonlinear decision function is described as(11)fx=sign∑i,j=1naiyiφxi·φxj+b.In order to calculate the value ofφ(xi)·φ(xj), the kernel function K(xi,xj) is used and the above function can be expressed as(12)fx=sign∑i,j=1naiyiKxi,xj+b.For the SVM, there are many kinds of kernel function, such as linear kernel, polynomial kernel, polynomial kernel, and radial basis function (RBF) kernel. Compared with other kernel functions, the RBF kernel can obtain the higher classification accuracy in many practical applications [34]. Therefore, the RBF kernel is used in this study.As previously mentioned, SVM is initially designed for binary classification. However, there are often many faults in the practical applications, which mean that it is necessary to develop a method to deal with a multiclassification problem. Currently, different methods have been developed for the multiclassification based on SVM, such as “one-against-one,” “one-against-all,” and directed acyclic graph (DAG). According to the comparison results obtained by Hsu and Lin [35], the “one-against-one” method is more suitable for practical use than other methods. For the sample set including c class, c(c-1)/2 SVM classifiers can be constructed by using “one-against-one” method and every SVM classifier is trained.


## 3.1. Fuzzy C-Means Clustering Algorithm

As an unsupervised machine learning method, FCM clustering algorithm was improved by Bezdek [23] in 1981 in order to solve the hard clustering problem by using fuzzy set theory. In the FCM clustering algorithm, membership degree function is used to indicate the extent to which each data point belongs to each cluster, and this information is also used to update the values of cluster centers [24]. Based on the concept of fuzzy C-partition, FCM clustering algorithm has been applied successfully in a wide variety of applications, such as image segmentation [25], data mining [26], thermal system monitoring [27], and fault diagnosis [28].For the sample setX=x1,x2,…,xn, the object of FCM clustering algorithm is to divide the sample set into c groups and obtain the cluster centers by minimizing the following dissimilarity function [29]:(2)minJmU,V=∑i=1c∑k=1nuikmdik2,where m is the fuzzy weighting parameter varying in the range [1,∞]. The bigger the m, the more fuzzy the final cluster result. U is fuzzy partition matrix, V is cluster center matrix, and n and c are the number of samples and cluster centers, respectively. uik is the fuzzy membership degree of the kth sample in the ith cluster and it should be meeting the following three constraints [30]:(3)uik∈0,1,1≤i≤c,1≤k≤n,∑i=1cuik=1,1≤k≤n,∑k=1nuik∈0,n,1≤i≤c.For the distancedik between kth sample xk and the centre of ith cluster Vi, it can be calculated by using Euclidean distance as follows:(4)dik=xk-Vi=∑j=1mxkj-Vij2,1≤i≤c,1≤k≤n.V
i
j can be calculated by utilizing the following formulation:(5)Vij=∑k=1nuikmxkj∑k=1nuikm,1≤i≤c,1≤j≤m.In essence, fuzzy cluster is performed through an iterative optimization by updating fuzzy membership degree [29]:(6)uiks+1=1∑j=1cdiks/djks2/m-1,1≤i≤c,1≤k≤n,where s is the iterative step.When the below requirement is met, we can stop iteration and obtain the cluster result:(7)Us+1-Us≤ε,where ε is the iterative threshold in the range [0,1].Based on above method, the cluster process of FCM clustering algorithm is virtually to determine the fuzzy membership degree and cluster centers through continuous iteration, which is shown in Figure3.Figure 3
Calculation flowchart of fuzzy C-means clustering algorithm.


## 3.2. Support Vector Machine Classification Model

Compared with the conventional classifiers, support vector machine (SVM), developed by Vapnik [31], can effectively solve the classification problem by implementing the structure risk minimization based on statistical learning theory. Nowadays, SVM has been widely and successfully applied to detection and diagnosis of machine conditions due to its high accuracy and good generalization for a smaller number of samples [32, 33].SVM is initially used to deal with binary classification problems. Its core idea is to transform the sample data from original space to a higher dimensional feature space through some nonlinear mapping functions and then find the optimal separating hyperplane in this feature space to realize linear classification. Figure4 shows the classification principle based on SVM for the nonlinear classification problem.Figure 4
Classification of two classes using SVM.For the nonlinear training sample data set including two classesxi,yi, i=1,2,…,n, x∈Rn, y∈-1,+1, n is the number of samples. The nonlinear mapping function φ(x) will be used to transform the sample data from original space to a higher dimensional feature space and the optimal separating hyperplane can be constructed to maximize the margin between the two classes by the following linear function:(8)fx=ω·φx+b,where ω is the normal vector of optimal separating hyperplane and b is a scalar.In essence, the solution of optimal separating hyperplane is the corresponding constraint optimization problem:(9)min12ω2+C∑i=1nζis.t.yiω·φxi+b+ζi≥1,1≤i≤nζi≥0,1≤i≤n,where C is the penalty factor that can realize the trade-off between empirical risk and confidence interval. ζi is slack factor.Combining the method of Lagrange multipliers, the above convex optimization problem can be simplified into the dual quadratic optimization problem:(10)maxLa=∑i=1nai-12∑i,j=1naiajyiyjφxi·φxjs.t.ai≥0,1≤i≤n∑i=1naiyi=0,where a is Lagrangian multiplier.Then, the nonlinear decision function is described as(11)fx=sign∑i,j=1naiyiφxi·φxj+b.In order to calculate the value ofφ(xi)·φ(xj), the kernel function K(xi,xj) is used and the above function can be expressed as(12)fx=sign∑i,j=1naiyiKxi,xj+b.For the SVM, there are many kinds of kernel function, such as linear kernel, polynomial kernel, polynomial kernel, and radial basis function (RBF) kernel. Compared with other kernel functions, the RBF kernel can obtain the higher classification accuracy in many practical applications [34]. Therefore, the RBF kernel is used in this study.As previously mentioned, SVM is initially designed for binary classification. However, there are often many faults in the practical applications, which mean that it is necessary to develop a method to deal with a multiclassification problem. Currently, different methods have been developed for the multiclassification based on SVM, such as “one-against-one,” “one-against-all,” and directed acyclic graph (DAG). According to the comparison results obtained by Hsu and Lin [35], the “one-against-one” method is more suitable for practical use than other methods. For the sample set including c class, c(c-1)/2 SVM classifiers can be constructed by using “one-against-one” method and every SVM classifier is trained.


## 4. Fusion Fault Diagnosis Modeling of EGT Based on FCM-SVM Approach

As the supervised artificial intelligence method, enough prior knowledge including specific fault types is necessary for SVM classification model. However, the fault types of gas turbine EGT cannot be determined directly based on experience in the practical applications due to the effects of many uncertain factors. In order to achieve the automatic monitor and diagnosis of EGT effectively, a fusion approach based on FCM clustering algorithm and SVM classification model (FCM-SVM) is proposed. Firstly, FCM clustering algorithm is used to realize clustering analysis and obtain the state patterns, which means that the preclassification of EGT is completed. Then, SVM multiclassification model is designed and used to carry out the online state pattern recognition and fault diagnosis of gas turbine EGT.Figure5 shows the fusion fault diagnosis framework of gas turbine EGT based on FCM-SVM approach. The detailed modeling processes are as follows.Figure 5
Fusion fault diagnosis framework of gas turbine EGT based on FCM-SVM approach.Step 1.
Generated sufficient EGT samples from the historical database and the essential preprocessing upon EGT data are carried out before data analysis, such as supplementary data, eliminating noise and outliers.Step 2.
According to the cluster process which is shown in Figure3, FCM clustering algorithm is used to obtain the initial clustering results of gas turbine EGT.Step 3.
Cluster validity indexλ(c) is used to evaluate the validity of clustering and determine the number of clusters. The c is optimum when λ(c) reaches its maximum value:(13)λc=∑i=1c∑k=1nuikmVi-x¯2/c-1∑i=1c∑k=1nuikmxk-Vi2/n-c,(14)x¯=∑i=1c∑k=1nuikmxkn.Step 4.
After obtaining the optimal clustering results, the fault diagnosis sample set including specific fault types can be established.Step 5.
SVM multiclassification model will be designed based on “one-against-one” method and trained by using fault diagnosis sample set.Step 6.
The measured EGT obtained from real gas turbine are preprocessed and inputted to the well trained SVM multiclassification model. Then we can get the final diagnostic results.Step 7.
The measured EGT also are stored into the historical database and used for later analysis.


## 5. Case Study and Discussions

In order to demonstrate the effectiveness of FCM-SVM approach introduced in this paper, the historical monitoring data of EGT from one industrial single shaft gas turbine will be analyzed as a case study in this section.


### 5.1. Sample Data

As an industrial single shaft gas turbine, Taurus70 is made in solar turbines incorporated and used for power generation. 12 thermocouple temperature sensors are used to measure the EGT and the average EGT is about 505°C when gas turbine operates in a normal state.Figure6 shows the changing curves of 12 EGT varying with time under normal running state of gas turbine. And the EGT profiles can be seen in Figure 7. From Figures 6 and 7, it is clear that there is significant difference between the measured outputs of different thermocouple temperature sensors at the same time even when the gas turbine is running in a normal state. Therefore, much feature information will be ignored which can decrease the fault diagnosis accuracy if only the average EGT is used to evaluate and analyse the health state of EGT. Considering the operating conditions of gas turbine, 490-group data including 4 classes are taken to establish the original sample set. 470 samples are selected randomly as training samples and the remaining 20 samples are selected as testing samples.Figure 6
Real-time measured EGT curve of gas turbine with normal condition.Figure 7
Real-time measured EGT profile of gas turbine with normal condition.


### 5.2. Optimal Clustering of EGT Based on FCM Clustering Algorithm

For the FCM clustering algorithm, it is very important to determine an appropriate number of clusters, which is called cluster validity problem. In this study, the cluster number is decided automatically by using the introduced cluster validity index which is shown in (13). Considering the computation complexity and accuracy, the scope of the number of clusters is commonly [2,n] in practical process and n is the number of samples. For the 470 training samples shown in Table 1, the FCM clustering algorithm can stepwise iterate from 2 to 21 clusters. Figure 8 shows the changing trends of cluster validity index λ(c) as a function of the number of clusters. Based on the result shown in Figure 8, it is clear that λ(c) increases at first and then drops down with the increases of the number of clusters. And λ(c) can reach its maximum value when the number of clusters is 4, which is in agreement with the real samples class. Therefore, it is concluded that the FCM clustering algorithm is suitable for optimal clustering of gas turbine EGT.Table 1
The original sample set of gas turbine EGT.
Number
T
1/°C
T
2/°C
T
3/°C
T
4/°C
T
5/°C
T
6/°C
T
7/°C
T
8/°C
T
9/°C
T
10/°C
T
11/°C
T
12/°C
1
493.6847
519.3781
510.2112
521.7342
488.1656
471.3032
475.5255
493.422
486.8506
513.879
518.3308
514.9266
2
497.8859
524.3516
513.879
524.09
493.422
475.5255
480.0065
497.8859
489.743
517.5453
523.043
519.1163
3
496.3108
521.996
512.3073
524.09
490.5316
473.151
478.9526
497.0984
488.1656
515.7123
520.9489
516.7598
4
500.2478
524.6133
513.355
526.9688
493.6847
477.1075
481.8502
499.723
491.8456
519.1163
523.043
520.6871
5
497.8859
521.7342
511.7833
522.7812
491.3201
473.9427
478.6891
496.5733
488.9544
515.9742
520.6871
516.7598
6
499.723
524.8751
514.1408
526.4454
494.21
476.3166
481.3235
499.1981
492.1084
518.8545
523.8282
520.1635
7
495.7857
520.6871
512.3073
523.3047
490.7945
473.4149
477.8983
495.2605
488.9544
514.1408
520.9489
516.7598
8
499.4606
525.1368
513.0931
526.4454
493.422
476.0529
481.5869
498.4109
491.5829
518.069
522.5195
519.3781
9
496.5733
523.043
513.355
524.09
491.8456
475.5255
480.5334
497.8859
489.2173
517.2835
521.7342
518.069
10
497.8859
523.043
513.879
525.1368
492.6339
475.2617
480.5334
497.6234
490.2688
517.2835
523.5664
519.9018
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
486
498.6733
525.3986
514.9266
525.922
493.9474
476.3166
482.1135
498.1484
490.5316
518.5927
522.7812
519.1163
487
498.6733
523.5664
514.9266
526.7071
493.1593
476.5803
482.1135
497.6234
490.7945
518.5927
523.043
519.64
488
498.6733
523.043
512.5692
524.8751
492.6339
475.7892
480.0065
497.0984
488.9544
518.069
521.7342
518.069
489
498.4109
522.5195
513.0931
524.09
491.8456
474.9979
480.2699
497.3609
489.4802
517.0217
521.2107
518.8545
490
499.4606
523.3047
514.6647
525.6602
492.6339
475.7892
481.0601
498.1484
491.0573
517.0217
521.7342
519.9018Figure 8
The effect of the number of clusters on cluster validity indexλ(c).Considering the high dimension characteristics of samples, it is difficult to realize graphical analysis directly. In this paper, three temperature uniformity indexes described by Mao [22] are used to further analyze and evaluate the cluster results of gas turbine EGT by using graphic approach. Figures 9–11 show the cluster results of gas turbine EGT based on FCM clustering algorithm. It may be clearly observed in Figures 9–11 that all the three temperature uniformity indexes of F1 class are relatively small (0°C≤H1≤30°C, 0°C≤H2≤30°C, and 0°C≤H3≤30°C). This means that 12 thermocouple temperature sensors give the almost same outputs, which belongs to normal state. Compared with F1 class, F2 class has the following characteristics: 40°C≤H1≤70°C, 30°C≤H2≤45°C, and 30°C≤H3≤45°C. The actual experimental results show that the fundamental reason for this phenomenon is turbine blade wear which can cause a difference of enthalpy drop between different turbine blade passages. For the F3 class, all the three temperature uniformity indexes are very large (140°C≤H1, 110°C≤H2, and 110°C≤H3) due to the effects of gas turbine load rejection. In addition, a careful inspection of Figures 9 and 10 reveals that the temperature uniformity index H1 is significantly larger than the other two indexes (500°C≤H1, 0°C≤H2≤30°C, and 0°C≤H3≤30°C). It means that one of the 12 thermocouple temperature sensors is fault which can result in a smaller output.Figure 9
Graphic clustering result by usingH1 and H2.Figure 10
Graphic clustering result by usingH1 and H3.Figure 11
Graphic clustering result by usingH2 and H3.


### 5.3. EGT Fault Diagnosis Based on SVM Classification Model

According to the optimal clustering results, the fault diagnosis training sample set including fault types can be established, which is shown in Table2. Then on this basis, we can develop 6 SVM classifiers based on “one-against-one” method. Table 3 shows the fault diagnosis performance of SVM multiclassification model for the training samples. From Table 3, it is clear that the fault diagnosis accuracy rate of trained SVM multiclassification model is 100% for the training samples, which means that the SVM multiclassification model has been well trained for fault diagnosis of gas turbine EGT.Table 2
Fault diagnosis training sample set of EGT.
Type
T
1/°C
T
2/°C
T
3/°C
T
4/°C
T
5/°C
T
6/°C
T
7/°C
T
8/°C
T
9/°C
T
10/°C
T
11/°C
T
12/°C
F1
493.6847
519.3781
510.2112
521.7342
488.1656
471.3032
475.5255
493.422
486.8506
513.879
518.3308
514.9266
497.8859
524.3516
513.879
524.09
493.422
475.5255
480.0065
497.8859
489.743
517.5453
523.043
519.1163
496.3108
521.996
512.3073
524.09
490.5316
473.151
478.9526
497.0984
488.1656
515.7123
520.9489
516.7598
500.2478
524.6133
513.355
526.9688
493.6847
477.1075
481.8502
499.723
491.8456
519.1163
523.043
520.6871
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
F2
489.743
511.5213
502.3465
512.0453
481.0601
464.9604
474.2065
490.2688
480.7968
502.8711
511.2593
506.2799
489.4802
512.0453
502.8711
513.617
482.1135
466.5472
473.9427
490.7945
482.1135
504.4445
512.5692
507.5906
487.9026
510.9973
501.0349
512.3073
480.0065
465.2249
472.3592
490.2688
480.5334
502.3465
510.9973
506.5421
488.4285
512.0453
503.6579
512.8312
480.7968
466.8116
473.151
490.7945
481.8502
502.8711
512.8312
507.5906
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
F3
444.8924
448.348
487.1684
474.5367
461.6006
427.5448
244.4364
333.8674
236.1834
323.3517
431.559
419.7637
444.0943
448.348
485.8542
472.955
461.3361
427.009
244.7409
334.1505
236.7968
322.4955
430.7567
420.3012
444.8924
448.6136
485.8542
474.0095
461.0715
427.5448
244.7409
334.4335
236.1834
322.781
431.559
420.5699
443.0299
448.0823
485.8542
473.2187
460.807
427.009
243.8271
333.8674
236.1834
321.9244
430.2217
420.3012
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
F4
472.3592
463.3729
454.6267
468.3976
−17.7778
471.3032
481.3235
476.8439
456.2188
462.5787
468.926
479.2161
472.3592
464.6959
454.3612
468.6618
−17.7778
471.0391
482.6401
477.1075
455.4228
462.0493
469.4544
477.3712
470.511
461.5197
451.174
466.8116
−17.7778
469.4544
480.7968
474.2065
454.3612
460.4602
467.869
478.4255
470.7751
462.314
451.7055
466.8116
−17.7778
469.7186
481.3235
475.5255
454.892
460.1953
467.869
477.1075
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮Table 3
Fault diagnosis performance of SVM multiclassification model for training samples.
Type
Number of training samples
Number of accurate diagnoses
Accuracy rate/%
F1
200
200
100
F2
200
200
100
F3
20
20
100
F4
50
50
100Then the testing samples shown in Table1 are used to further demonstrate the effectiveness of SVM multiclassification model. Table 4 shows the comparison between actual results and fault diagnosis results by using the well trained SVM multiclassification model for testing samples. Based on the results shown in Table 4, it is demonstrated that the well trained SVM multiclassification model can effectively diagnose the fault of gas turbine EGT with a 95% accuracy rate for the testing samples. Besides, the reason of misclassification is that the sample data is obtained when the turbine blade wear or corrosion is not severe. In order to compare with other models, backpropagation (BP) neural network model is also employed to make the same fault diagnosis and the results are also listed in Table 4. The comparative analysis shows that SVM classification model can improve the fault diagnosis accuracy of gas turbine EGT significantly compared with BP neural network model. All these indicate that SVM is more suitable for fault diagnosis of gas turbine EGT.Table 4
The comparison results of different fault diagnosis models for testing samples.
Number
SVM1
SVM2
SVM3
SVM4
SVM5
SVM6
SVM model
BP model
Actual results
1
F1
F1
F1
F2
F2
F3
F1
F1
F1
2
F1
F1
F1
F2
F2
F3
F1
F1
F1
3
F1
F1
F1
F2
F2
F3
F1
F1
F1
4
F1
F1
F1
F2
F2
F3
F1
F2
F1
5
F1
F1
F1
F2
F2
F3
F1
F1
F1
6
F2
F1
F1
F2
F2
F3
F2
F2
F2
7
F1
F1
F1
F2
F2
F3
F1
F1
F2
8
F2
F1
F1
F2
F2
F3
F2
F2
F2
9
F2
F1
F1
F2
F2
F3
F2
F2
F2
10
F2
F1
F1
F2
F2
F3
F2
F1
F2
11
F2
F3
F1
F3
F2
F3
F3
F3
F3
12
F2
F3
F1
F3
F2
F3
F3
F3
F3
13
F2
F3
F1
F3
F2
F3
F3
F3
F3
14
F2
F3
F1
F3
F2
F3
F3
F3
F3
15
F2
F3
F1
F3
F2
F3
F3
F3
F3
16
F1
F3
F4
F3
F4
F4
F4
F4
F4
17
F1
F3
F4
F3
F4
F4
F4
F4
F4
18
F1
F3
F4
F3
F4
F4
F4
F4
F4
19
F1
F3
F4
F3
F4
F4
F4
F4
F4
20
F1
F3
F4
F3
F4
F4
F4
F4
F4


## 5.1. Sample Data

As an industrial single shaft gas turbine, Taurus70 is made in solar turbines incorporated and used for power generation. 12 thermocouple temperature sensors are used to measure the EGT and the average EGT is about 505°C when gas turbine operates in a normal state.Figure6 shows the changing curves of 12 EGT varying with time under normal running state of gas turbine. And the EGT profiles can be seen in Figure 7. From Figures 6 and 7, it is clear that there is significant difference between the measured outputs of different thermocouple temperature sensors at the same time even when the gas turbine is running in a normal state. Therefore, much feature information will be ignored which can decrease the fault diagnosis accuracy if only the average EGT is used to evaluate and analyse the health state of EGT. Considering the operating conditions of gas turbine, 490-group data including 4 classes are taken to establish the original sample set. 470 samples are selected randomly as training samples and the remaining 20 samples are selected as testing samples.Figure 6
Real-time measured EGT curve of gas turbine with normal condition.Figure 7
Real-time measured EGT profile of gas turbine with normal condition.


## 5.2. Optimal Clustering of EGT Based on FCM Clustering Algorithm

For the FCM clustering algorithm, it is very important to determine an appropriate number of clusters, which is called cluster validity problem. In this study, the cluster number is decided automatically by using the introduced cluster validity index which is shown in (13). Considering the computation complexity and accuracy, the scope of the number of clusters is commonly [2,n] in practical process and n is the number of samples. For the 470 training samples shown in Table 1, the FCM clustering algorithm can stepwise iterate from 2 to 21 clusters. Figure 8 shows the changing trends of cluster validity index λ(c) as a function of the number of clusters. Based on the result shown in Figure 8, it is clear that λ(c) increases at first and then drops down with the increases of the number of clusters. And λ(c) can reach its maximum value when the number of clusters is 4, which is in agreement with the real samples class. Therefore, it is concluded that the FCM clustering algorithm is suitable for optimal clustering of gas turbine EGT.Table 1
The original sample set of gas turbine EGT.
Number
T
1/°C
T
2/°C
T
3/°C
T
4/°C
T
5/°C
T
6/°C
T
7/°C
T
8/°C
T
9/°C
T
10/°C
T
11/°C
T
12/°C
1
493.6847
519.3781
510.2112
521.7342
488.1656
471.3032
475.5255
493.422
486.8506
513.879
518.3308
514.9266
2
497.8859
524.3516
513.879
524.09
493.422
475.5255
480.0065
497.8859
489.743
517.5453
523.043
519.1163
3
496.3108
521.996
512.3073
524.09
490.5316
473.151
478.9526
497.0984
488.1656
515.7123
520.9489
516.7598
4
500.2478
524.6133
513.355
526.9688
493.6847
477.1075
481.8502
499.723
491.8456
519.1163
523.043
520.6871
5
497.8859
521.7342
511.7833
522.7812
491.3201
473.9427
478.6891
496.5733
488.9544
515.9742
520.6871
516.7598
6
499.723
524.8751
514.1408
526.4454
494.21
476.3166
481.3235
499.1981
492.1084
518.8545
523.8282
520.1635
7
495.7857
520.6871
512.3073
523.3047
490.7945
473.4149
477.8983
495.2605
488.9544
514.1408
520.9489
516.7598
8
499.4606
525.1368
513.0931
526.4454
493.422
476.0529
481.5869
498.4109
491.5829
518.069
522.5195
519.3781
9
496.5733
523.043
513.355
524.09
491.8456
475.5255
480.5334
497.8859
489.2173
517.2835
521.7342
518.069
10
497.8859
523.043
513.879
525.1368
492.6339
475.2617
480.5334
497.6234
490.2688
517.2835
523.5664
519.9018
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
486
498.6733
525.3986
514.9266
525.922
493.9474
476.3166
482.1135
498.1484
490.5316
518.5927
522.7812
519.1163
487
498.6733
523.5664
514.9266
526.7071
493.1593
476.5803
482.1135
497.6234
490.7945
518.5927
523.043
519.64
488
498.6733
523.043
512.5692
524.8751
492.6339
475.7892
480.0065
497.0984
488.9544
518.069
521.7342
518.069
489
498.4109
522.5195
513.0931
524.09
491.8456
474.9979
480.2699
497.3609
489.4802
517.0217
521.2107
518.8545
490
499.4606
523.3047
514.6647
525.6602
492.6339
475.7892
481.0601
498.1484
491.0573
517.0217
521.7342
519.9018Figure 8
The effect of the number of clusters on cluster validity indexλ(c).Considering the high dimension characteristics of samples, it is difficult to realize graphical analysis directly. In this paper, three temperature uniformity indexes described by Mao [22] are used to further analyze and evaluate the cluster results of gas turbine EGT by using graphic approach. Figures 9–11 show the cluster results of gas turbine EGT based on FCM clustering algorithm. It may be clearly observed in Figures 9–11 that all the three temperature uniformity indexes of F1 class are relatively small (0°C≤H1≤30°C, 0°C≤H2≤30°C, and 0°C≤H3≤30°C). This means that 12 thermocouple temperature sensors give the almost same outputs, which belongs to normal state. Compared with F1 class, F2 class has the following characteristics: 40°C≤H1≤70°C, 30°C≤H2≤45°C, and 30°C≤H3≤45°C. The actual experimental results show that the fundamental reason for this phenomenon is turbine blade wear which can cause a difference of enthalpy drop between different turbine blade passages. For the F3 class, all the three temperature uniformity indexes are very large (140°C≤H1, 110°C≤H2, and 110°C≤H3) due to the effects of gas turbine load rejection. In addition, a careful inspection of Figures 9 and 10 reveals that the temperature uniformity index H1 is significantly larger than the other two indexes (500°C≤H1, 0°C≤H2≤30°C, and 0°C≤H3≤30°C). It means that one of the 12 thermocouple temperature sensors is fault which can result in a smaller output.Figure 9
Graphic clustering result by usingH1 and H2.Figure 10
Graphic clustering result by usingH1 and H3.Figure 11
Graphic clustering result by usingH2 and H3.


## 5.3. EGT Fault Diagnosis Based on SVM Classification Model

According to the optimal clustering results, the fault diagnosis training sample set including fault types can be established, which is shown in Table2. Then on this basis, we can develop 6 SVM classifiers based on “one-against-one” method. Table 3 shows the fault diagnosis performance of SVM multiclassification model for the training samples. From Table 3, it is clear that the fault diagnosis accuracy rate of trained SVM multiclassification model is 100% for the training samples, which means that the SVM multiclassification model has been well trained for fault diagnosis of gas turbine EGT.Table 2
Fault diagnosis training sample set of EGT.
Type
T
1/°C
T
2/°C
T
3/°C
T
4/°C
T
5/°C
T
6/°C
T
7/°C
T
8/°C
T
9/°C
T
10/°C
T
11/°C
T
12/°C
F1
493.6847
519.3781
510.2112
521.7342
488.1656
471.3032
475.5255
493.422
486.8506
513.879
518.3308
514.9266
497.8859
524.3516
513.879
524.09
493.422
475.5255
480.0065
497.8859
489.743
517.5453
523.043
519.1163
496.3108
521.996
512.3073
524.09
490.5316
473.151
478.9526
497.0984
488.1656
515.7123
520.9489
516.7598
500.2478
524.6133
513.355
526.9688
493.6847
477.1075
481.8502
499.723
491.8456
519.1163
523.043
520.6871
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
F2
489.743
511.5213
502.3465
512.0453
481.0601
464.9604
474.2065
490.2688
480.7968
502.8711
511.2593
506.2799
489.4802
512.0453
502.8711
513.617
482.1135
466.5472
473.9427
490.7945
482.1135
504.4445
512.5692
507.5906
487.9026
510.9973
501.0349
512.3073
480.0065
465.2249
472.3592
490.2688
480.5334
502.3465
510.9973
506.5421
488.4285
512.0453
503.6579
512.8312
480.7968
466.8116
473.151
490.7945
481.8502
502.8711
512.8312
507.5906
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
F3
444.8924
448.348
487.1684
474.5367
461.6006
427.5448
244.4364
333.8674
236.1834
323.3517
431.559
419.7637
444.0943
448.348
485.8542
472.955
461.3361
427.009
244.7409
334.1505
236.7968
322.4955
430.7567
420.3012
444.8924
448.6136
485.8542
474.0095
461.0715
427.5448
244.7409
334.4335
236.1834
322.781
431.559
420.5699
443.0299
448.0823
485.8542
473.2187
460.807
427.009
243.8271
333.8674
236.1834
321.9244
430.2217
420.3012
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
F4
472.3592
463.3729
454.6267
468.3976
−17.7778
471.3032
481.3235
476.8439
456.2188
462.5787
468.926
479.2161
472.3592
464.6959
454.3612
468.6618
−17.7778
471.0391
482.6401
477.1075
455.4228
462.0493
469.4544
477.3712
470.511
461.5197
451.174
466.8116
−17.7778
469.4544
480.7968
474.2065
454.3612
460.4602
467.869
478.4255
470.7751
462.314
451.7055
466.8116
−17.7778
469.7186
481.3235
475.5255
454.892
460.1953
467.869
477.1075
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮Table 3
Fault diagnosis performance of SVM multiclassification model for training samples.
Type
Number of training samples
Number of accurate diagnoses
Accuracy rate/%
F1
200
200
100
F2
200
200
100
F3
20
20
100
F4
50
50
100Then the testing samples shown in Table1 are used to further demonstrate the effectiveness of SVM multiclassification model. Table 4 shows the comparison between actual results and fault diagnosis results by using the well trained SVM multiclassification model for testing samples. Based on the results shown in Table 4, it is demonstrated that the well trained SVM multiclassification model can effectively diagnose the fault of gas turbine EGT with a 95% accuracy rate for the testing samples. Besides, the reason of misclassification is that the sample data is obtained when the turbine blade wear or corrosion is not severe. In order to compare with other models, backpropagation (BP) neural network model is also employed to make the same fault diagnosis and the results are also listed in Table 4. The comparative analysis shows that SVM classification model can improve the fault diagnosis accuracy of gas turbine EGT significantly compared with BP neural network model. All these indicate that SVM is more suitable for fault diagnosis of gas turbine EGT.Table 4
The comparison results of different fault diagnosis models for testing samples.
Number
SVM1
SVM2
SVM3
SVM4
SVM5
SVM6
SVM model
BP model
Actual results
1
F1
F1
F1
F2
F2
F3
F1
F1
F1
2
F1
F1
F1
F2
F2
F3
F1
F1
F1
3
F1
F1
F1
F2
F2
F3
F1
F1
F1
4
F1
F1
F1
F2
F2
F3
F1
F2
F1
5
F1
F1
F1
F2
F2
F3
F1
F1
F1
6
F2
F1
F1
F2
F2
F3
F2
F2
F2
7
F1
F1
F1
F2
F2
F3
F1
F1
F2
8
F2
F1
F1
F2
F2
F3
F2
F2
F2
9
F2
F1
F1
F2
F2
F3
F2
F2
F2
10
F2
F1
F1
F2
F2
F3
F2
F1
F2
11
F2
F3
F1
F3
F2
F3
F3
F3
F3
12
F2
F3
F1
F3
F2
F3
F3
F3
F3
13
F2
F3
F1
F3
F2
F3
F3
F3
F3
14
F2
F3
F1
F3
F2
F3
F3
F3
F3
15
F2
F3
F1
F3
F2
F3
F3
F3
F3
16
F1
F3
F4
F3
F4
F4
F4
F4
F4
17
F1
F3
F4
F3
F4
F4
F4
F4
F4
18
F1
F3
F4
F3
F4
F4
F4
F4
F4
19
F1
F3
F4
F3
F4
F4
F4
F4
F4
20
F1
F3
F4
F3
F4
F4
F4
F4
F4


## 6. Conclusions

Considering the distribution characteristics of gas turbine EGT and its effect on the health state of gas turbine, a fusion approach based on FCM clustering algorithm and SVM classification model (FCM-SVM) is proposed and successfully applied to an industrial gas turbine in this paper. In the analysis presented in this study, it is demonstrated that FCM-SVM based approach can make full use of the unsupervised feature extraction ability of FCM clustering algorithm and the sample classification generalization properties of SVM multiclassification model, which offers an effective way to realize the online condition recognition and fault diagnosis of gas turbine EGT. In the concrete implementation process, the introduced FCM clustering algorithm is a good alternative to achieve automatic identification of the fault types of gas turbine EGT. In other words, it is effective to overcome the influence of experience judgment on fault types. Besides, the introduction of SVM multiclassification model has a great potential to improve the fault diagnosis performance of gas turbine EGT. It is worth noticing that the study of this paper is only focused on researching the artificial intelligence approach for the condition recognition and fault diagnosis of gas turbine EGT but ignores the effects of many other parameters such as inlet temperature of gas turbine. Therefore, more studies and improvement about the application of this approach are needed further.



---

*Source: 240267-2015-06-08.xml*