Advanced Mathematical and Numerical Methods in Control and Optimization for Smart Grids
View this Special IssueResearch Article  Open Access
Abnormal Detection of Wind Turbine Based on SCADA Data Mining
Abstract
In order to reduce the curse of dimensionality of massive data from SCADA (Supervisory Control and Data Acquisition) system and remove data redundancy, the grey correlation algorithm is used to extract the eigenvectors of monitoring data. The eigenvectors are used as input vectors and the monitoring variables related to the unit state as output vectors. The genetic algorithm and cross validation method are used to optimize the parameters of the support vector regression (SVR) model. A high precision prediction is carried out, and a reasonable threshold is set up to alarm the fault. The condition monitoring of the wind turbine is realized. The effectiveness of the method is verified by using the actual fault data of a wind farm.
1. Introduction
SCADA system changes the operation mode of wind farm systems with healthy working environment and reduces the costs of operation and maintenance. However, a large set of high dimensions and many types of data are not fully utilized or developed; only stay on realtime data and historical data reporting statistics are typically monitored or gathered. Therefore, it is important to make full use of the data collected by wind power centralized control center to collect the data of the massive wind turbines, conduct state monitoring of the turbines, and predict their condition and life [1, 2]. Several surveys of WT failures have been conducted in the last two decades to identify failure rates and associated downtime for different subassemblies. However, the different taxonomies used by different turbine manufacturers, wind farm operators, and researchers make comparisons between these surveys challenging. The evaluation of 15 years of data from the German “250 MW Wind” programme [3] and >95% of all the turbines operating between 1997 and 2005 in Sweden [4] gave first insights into the reliability of the first onshore WTs. The German turbines had an average availability of about 98%. An average failure rate of 0.4 failures per turbine per year resulted in an average downtime of 130 hours per turbine per year for the Swedish turbines. A distinctive difference between failure rate and downtime distribution in subassembly groups was identified. The electrical and electronic control systems were identified as the most failureprone, but gearbox and generator failures caused the longest downtime.
Many scholars have researched large power wind turbine monitoring and fault diagnoses [5], based on statistical learning method to detect abnormal situations through the wind turbine response model of the weighted least squares support vectorbased wind power generator and external regression conditions [6]. The results show that the model is better than conventional forecasting methods. Pandit and Infield [7] used an indepth analysis of commonly used stationary covariance functions in which wind turbine power curve was used, where GPbased power curve was constructed using different stationary covariance functions, and after that, a comparative analysis was carried out in order to identify the most effective covariance function. The commonly used squared exponential covariance function is taken as the benchmark, against which other covariance functions are assessed. The results show that the performance (in terms of model accuracy and uncertainty) of GP fitted power curve models based on rational quadratic covariance functions is almost the same as for the most commonly used squared exponential function. The studies of Astolfi et al.[8, 9] are a catalog of generalizable methods for studying wind turbine power curve upgrades. In particular, from the study of the selected test cases, it arises that complex wind conditions might affect wind turbine operation such that the production improvement is nonnegligibly different from what can be estimated under the hypothesis of ideal wind conditions. Wan et al. [10] proposed wind form using wavelet based on energy function for asymmetrical fault detection in doubly fed induction generator. The proposed method not only detects the fault within one and half cycle of fundamental wave but also reveals the effectiveness under timevarying conditions. Turbine condition monitoring (TCM) through vibration analysis has pros and cons: basically high diagnostic power against high cost and high complexity for elaborating the information [11] from the data stream into knowledge. Chunyu et al. [12] put forward a dynamic prediction model of wind turbine blade failure based on the grey theory. The relative error between prediction and field investigation data is less than 5%, meeting the actual needs of engineering and verifying the effectiveness and applicability of the proposed algorithm. The main contribution of Chakkor et al. [13] is designing an intelligent wireless remote monitoring and control system according to features and requirements of wind turbines. This system based on IP communication combines Web and database client/server technology to copy data measurements received from the different sensors installed in the wind turbine machines. Eggers et al. [14] used Hotelling T2 control chart and an automatic relevance neural network to analyze the wind turbine power to identify wind turbine detection faults. Zhang [15], combined with AHP and variable weight theory, used a wind turbine performance evaluation model based on Grey Theory and a variable weight fuzzy comprehensive evaluation. However, these studies did not consider the correlation and coupling between the components of the unit, which makes the model inaccurate. Zhang et al. [16], based on SVR prediction model, helped to establish a prediction model, which takes the amount of SCADA systems as input and the active power of the unit as output. The disadvantage of this model is that the feature extraction of the highdimensional input vectors is not easy, and the power is used as the only standard to diagnose the state of the unit. BP neural network is used to model and predict gearbox and generator [17], and multiagent method is used to synthetically analyze the diagnosis results of different components, giving the overall operation status of the unit. However, the use of neural network modelling requires timeconsuming learning process, and the selection of learning samples lacks basis.
Being based on statistical analysis, it commonly requires vast datasets for providing meaningful indications: the most common opinion therefore is that SCADA can detect incipient faults at a late stage. Astolfi et al. [18] employed artificial neural networks, for their capability in reconstructing nonlinear dependency between inputs and outputs, and formulated simple models for the diagnosis of occurring faults at the level of gearbox. The datasets employed have the 10minute sampling time of the common SCADA control systems; the gearbox vibrations and the gearbox temperatures are selected as target output to model. It will be shown that the time resolution of SCADA is too coarse for reliable vibration analysis, which should be rather observed at its proper time scale (several Hz). At present, data mining methods such as clustering and statistical model are widely used in domestic and foreign enterprises, but their cleaning process is complicated and the cleaning conditions are harsh [19]. Therefore, in order to make a reliable analysis of the power generation performance of wind turbines, an efficient and versatile cleaning method is urgently needed.
In view of this, this paper firstly extracts the features from the massive and highdimensional data collected by the SCADA system, removes the irrelevant and redundant parameters of the operation state of the unit, and improves the monitoring accuracy of the wind turbine by improving the model input. The reasonable threshold is selected to alarm the abnormal state of wind turbine to avoid false alarm and untimely alarm.
The paper is organized in three sections. Section 2 discusses feature selection and sparse learning technology to reduce the dimension of the operation parameters of the SCADA system, remove the independent and redundant parameters of the operating state of the wind turbine, and retain the related characteristic parameters. In Section 3, the multiinput and multioutput SVR model, which takes the active power, the speed of the blade, and the pitch 1 angle as the output vector and the characteristic parameter as the input vector, is established. Cross validation (CV) is combined with a genetic algorithm (GA) for parameter optimization. In Sections 4 and 5, the proposed method is applied to the industrial data. Performance of the proposed model is also discussed. Section 6 concludes the paper.
2. Data Mining of Characteristic Parameters for Wind Turbines
The data collected and recorded by the SCADA system of the wind turbine has highdimensional characteristics. In this paper, 74 digits of the wind turbine components are selected. The method of feature parameter data mining reduces the number of features and dimension disaster so that the generalization ability of the model is stronger and the overfitting phenomenon is reduced. The commonly used methods for selecting characteristic parameters include principal component analysis (PCA) [20], the Pearson correlation coefficient [14], and the random forest method [21]. When the data are highdimensional vectors, the calculation of PCA is complicated and it is most suitable for linear data. Pearson correlation coefficient is only sensitive to the disadvantages of the most obvious linear relationship. The random forest method is prone to an overfitting phenomenon. Therefore, in this paper, a data mining algorithm based on the grey correlation degree [22] is proposed to overcome the above shortcomings and to improve the accuracy and effectiveness of wind turbine operation state assessments.
2.1. Extraction of the Characteristic Parameters Based on Grey Relational Grade
There are 74 variables in wind turbine information recorded by the SCADA system. The acquisition interval is 10 minutes, as shown in Figure 1. Figure 1 shows the monitoring variables collected by the wind turbine SCADA system and its response code. The object of study in this paper is that the wind turbines are in the condition of unlimited power and healthy operation. They have some monitoring quantities such as the control mode and alarm of some parameters recorded in the SCADA system, speed mode, state of shaft 1, shaft 2. and shaft 3 converters, etc. Variables in an invariant state can be ignored. Table 1 is part of the parameter alarm information of GE wind turbine manufacturer. To preprocess these eigenvectors, we must remove these eigenvectors to avoid the disaster of dimensionality caused by too many features.

Wind turbine operation is mainly reflected in active power, rotor speed, and pitch angle; these three parameters are used as input vectors. We take the pitch angle 1 as an example to monitor the pitch angle. The grey correlation with other variables is calculated to reduce the wind turbine data dimension while ensuring the smallest loss of information. The concrete steps of extracting the operating characteristic parameters using the grey correlation degree are as follows:(1)The characteristic set of the wind turbine operating state is .(2)According to the parameters of the primary wind turbine , the corresponding parameters are extracted from the SCADA system as the sample set of the grey correlation degree :where is the parameter of the samples. The degree of correlation between the calculated parameters is found as follows.
Determine the reference sequence and the comparison sequence according to the training sample .
The absolute difference between and sequences is calculated from the sample set as
The absolute difference is used to calculate the maximum difference and minimum difference of level two, respectively, as
is calculated, and the sequence in the moment of for correlation coefficient is compared by
2.2. Mining of Characteristic Parameter Data of Wind Turbines
The grey correlation analysis mentioned above is applied to 1.6 MW wind turbines. By preprocessing the eigenvectors, 23 variables related to the operating state of the wind turbine are selected, and the values of these variables are extracted from the SCADA system. Based on operational data from January 1, 2015, to January 1, 2017, the sample set features selected by grey correlation analysis and the characteristic parameters of the grey correlation matrix for color map are as shown in Figure 2. According to Figure 1, the grey correlation degree between each primary parameter of the power, rotor speed, and pitch angle is different, which makes it feasible to excavate the characteristic parameters of the generator set. In this paper, variables with correlation degree greater than 0.5 are selected as input of monitoring variables. The set of characteristic parameters is shown in Tables 2–4.



3. Prediction Model Based on Support Vector Regression
After the data effectiveness analysis and dimensionality reduction are conducted, the parameters of the wind turbine are regressed. SVR [23] algorithm of structural risk minimization criterion solves the practical problems of small sample, nonlinearity, and high dimension and overcomes the shortcomings of the indetermination of the network structure and local minima, over learning and under learning. Therefore, this paper chooses SVR algorithm to build regression prediction model. The specific algorithm is as follows.
Set a given sample training set for , is the input vector, and is the output vector. Nonlinear mapping should be used for the nonlinear SVR model . The mapping sample sets are used to feature spaces . The optimal decision function is as follows:where is the characteristic space weight coefficient vector and is the bias. It is assumed that all training samples can be in precision with linear functions at accuracy. According to the principle of structural risk minimization, the problem can be formulated aswhere is the relaxation factor. Introducing the Lagrange function, the optimization problem in the dual space is used to obtain the following formula:where , and are the Lagrange multipliers and is the penalty factor respective to , and . Find the partial derivative and make it equal to 0 and bring the derivative into the Lagrange function:
Using the positive definite matrix theorem, the inner product is replaced by a kernel function . Therefore, the SVR function can be obtained as follows:
In the kernel function, the structure of the radial basis function (RBF) kernel is simple, and its generalization ability is better. Based on this, the kernel function of the model is selected as the radial basis function. , where is the kernel width.
3.1. Genetic Algorithm (GA)
In this model, the penalty coefficient and parameter of the kernel function affect the SVR precision. Therefore, GA is used to optimize the parameters of the SVR model, which is based on the natural selection and genetic mechanism of the theory of biological evolution by Darwin [24] to find the optimal solution. The main process is to encode the solution to the problem. There are two ways to code the solution individually, including binary coding and real number coding, which essentially maps the solution space to the chromosome space. Then, a reasonable initial population is generated in these solution spaces, and individuals are selected according to fitness function, genetic selection, crossover, and mutation operation. The individual with high fitness value is kept and vice versa. This new generation of offspring retains the advantages of the previous generation, whereas the last generation did not. This process is iterated many times until the optimal solution is obtained.
3.2. Cross Validation (CV)
In machine learning, CV is mainly used for model performance evaluation and learning. The basic principle is that the original sample is divided into a training set and a validation set, and then the training set is used to train the model. The model validated by the test set is obtained from the training model. As a performance index evaluation model, CV considers the training error as well as the generalization error. The most common CV method is kfolding cross validation (Kfold CV), and the specific algorithms are as follows:(1)The sample are divided into subsets that are not intersected, and the number of samples is . are remembered.(2)For each model , do the following. Training set and model . Get the corresponding hypothesis function . Validation set is used to calculate the generalization error .(3)The average generalization error of each model is calculated, and the model with the least generalization error is selected.
3.3. SVR Parameter Selection Based on CV and GA
The selection of the parameters of the SVR model is essentially the optimization of the model. The algorithm of Kfold CV and GA is used to optimize the parameters of SVR, and . The model as shown in Figure 3 is as follows:(1)The SVR parameters are coded to form the initial population.(2)For population decoding, we calculate the fitness of individuals based on the Kfold CV method. In this paper, , the minimum mean square error MSE of samples , is chosen as the fitness function value of the GA algorithm.(3)Judge whether or not to meet the terminating condition if it is satisfied to turn (5); otherwise, proceed to (4).(4)Update the population by selection, cross, and variation; then, return to (2).(5)The optimal and optimal model is output.
4. Condition Monitoring Based on SVR Parameter Optimization
4.1. Data Processing
In this paper, power, rotor speed, and pitch angle are taken as the output vectors and other feature parameters are taken as input vectors, and then multipleinput and multipleoutput SVR model is established. The accuracy of the proposed model is verified by running data of the wind turbine for four months.(1)According to the fault information recorded by the SCADA system, the samples of the maintenance shutdown due to the failure of the wind turbine and the samples of the less power operation are eliminated.(2)Consider that the cutting wind speed is 3 m/s, the rated wind speed is 12 m/s, and the cutting wind speed is 25 m/s. According to the actual power curve of wind turbine, the wind speed range selected in this paper is 3 to 25.(3)To eliminate the magnitude of interference between the parameters, the parameters are normalized to [0, 1], according to the dimensions.
4.2. Model Establishment
In order to prove the validity of the model, this paper selects four months effective data of wind turbine to predict. In this model, the principle of cross validation selection first considers the minimum error MSE. According to the errors of MSE, to avoid the occurrence of the learning state, a group of smaller penalty parameters is selected as the best parameter. From Figure 4, we can see that 5CV and GA model are the best. The average fitness curve in Figure 4 indicates the average fitness of all the individuals in each generation. The best fitness curve represents the maximum fitness of all individuals in each generation. The convergence of the fitness curve is very fast, and the convergence level of the final fitness curve is relatively consistent, which reflects the optimization of SVR parameters. When the power, rotor speed, and pitch angle are output, the best parameter is applied to the SVR model. The comparison between the actual and predicted values of the SVR model is shown in Figure 5. Table 5 lists the power, rotor speed, and pitch error. The mean relative error values indicating the good prediction accuracy and stability of the SVR model CVGA algorithm is shown in Table 5.
(a)
(b)
(c)
(d)
(a)
(b)
(c)

4.3. Discriminating the Health Status of Wind Turbine Based on Threshold
According to the trained SVR model, the observed values of power, rotor speed, and pitch angle can be obtained from the current input vector, and the distance between the measured value and the observed value extracted from SCADA system is . According to the distance between the measured value and the observed value, the threshold is compared with the set threshold to distinguish the operating state of the wind turbine. The distance from the observed value to the measured value in the SCADA system is defined as
If is selected too low, the algorithm is too sensitive to the change of the operating state of the wind turbine and it is prone to misjudge the results. If selection is too large, the prediction time will be reduced and the detection rate of the abnormal operating state will be affected. To solve this problem, it is necessary to select the appropriate threshold . In this paper, twoyear SCADA data of the 1.6 MW wind turbine is analyzed, and the appropriate threshold is determined. From Table 5, it can be seen that threshold selection will cause the operating state of the wind turbine to not be normally recognized. To find a lower detection rate and the misjudgment rate, the threshold is selected as (in Table 6).

5. Example Analysis
At 5:32 on August 7, 2016, a wind turbine in Hebei Province went into shutdown due to the SCADA system failure alarm. After checking the pitch gear of blade 1 of the wind turbine, the wind turbine went into shutdown due to failure. 950 sets of data were extracted from the SCADA system before the blower alarm stopped. 950 sets of collected data are input into the model, as shown in Figure 6. As can be seen from the figure, close to 150 data points were detected as abnormal points before the wind turbine shut down. The situation of power restriction indicates that some of the wind turbines have begun to deteriorate at these times, and the unit has issued abnormal alarms. The reduced power generated by the unit indicates that the model can give a hint before the failure occurs. Therefore, the proposed model is effective for the state monitoring and fault prediction of the wind turbine, and it can avoid the continued deterioration of the fault and the influence on the safe operation of the power grid.
6. Conclusion
In order to extract relevant state from massive and highdimensional data of the SCADA system, realizing the monitoring of the state of the wind turbine, a grey correlation degree is proposed based on data mining technology to extract characteristic parameters of the wind turbine’s operating state, which reduces the data dimensions and computation. To improve the precision, GA and CV are combined to optimize the parameters of the regression model. To verify the validity of the model, the threshold of the SVR model is analyzed, and the model is applied to wind farm. The results show the following:(1)By establishing a data mining model of the characteristic parameters based on the grey correlation analysis, we extract parameters that are more related to the power, rotor speed, and pitch angle, effectively avoiding “dimension disaster.”(2)By comparing the power, rotor speed, pitch angle regression model, and the measured values, the results showed that the average relative error of the SVR model is low. The regression model has high accuracy and generalization ability; it can be applied to wind turbine anomaly distinguishing analysis.(3)Applying the model to practice, the analysis of the model results and the SCADA system can be used to record the measured values. The results show that when using the distance threshold to choose the appropriate conditions, the wind turbine condition monitoring can reflect the operating status of the wind turbine to provide technical references for online monitoring of wind turbines.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Liang Tao provided data and ideas as well as optimization algorithm; Qian Siqi mined and analyzed data; Zhang Yingjuan used GAoptimized SVR algorithm to predict monitoring variables; and Shi Huan analyzed and summarized the experimental results.
References
 T. Liang, Z. Yuan, C. Mei et al., “Wind power remote centralized intelligent monitoring system design and management,” Automation and Instrumentation, vol. 30, no. 11, pp. 50–53, 2015. View at: Google Scholar
 P. Tchakoua, R. Wamkeue, M. Ouhrouche, F. SlaouiHasnaoui, T. Tameghe, and G. Ekemb, “Wind turbine condition monitoring: stateoftheart review, new trends, and future challenges,” Energies, vol. 7, no. 4, pp. 2595–2630, 2014. View at: Publisher Site  Google Scholar
 B. Hahn, M. Durstewitz, and K. Rohrig, “Reliability of wind turbines experience of 15 year with 1500 WTs,” in Wind Energy, J. Peinke, P. Schaumann, and S. Barth, Eds., pp. 329–332, Springer, Berlin, Germany, 2007. View at: Google Scholar
 J. Ribrant and L. M. Bertling, “Survey of failures in wind power systems with focus on Swedish wind power plants during 1997–2005,” IEEE Transactions on Energy Conversion, vol. 22, no. 1, pp. 167–173, 2007. View at: Publisher Site  Google Scholar
 Lingwu Energy, H. Naishan, and R. Li, “Shortterm wind power prediction based on cloud support vector machine model,” Power Automation Equipment, vol. 33, no. 7, pp. 34–38, 2013. View at: Google Scholar
 H. H. Yang, M. L. Huang, and P. C. Huang, “Detection of wind turbine faults using a data mining approach,” Journal of Energy Engineering, vol. 142, no. 3, Article ID 04015021, 2016. View at: Publisher Site  Google Scholar
 R. K. Pandit and D. Infield, “Comparative analysis of Gaussian process (GP) power curve models based on different stationary covariance functions for the purpose of improving model accuracy,” Renewable Energy, vol. 140, pp. 190–202, 2019. View at: Publisher Site  Google Scholar
 D. Astolfi, F. Castellani, and L. Terzi, “Wind turbine power curve upgrades,” Energies, vol. 11, no. 5, p. 1300, 2018. View at: Publisher Site  Google Scholar
 D. Astolfi and F. Castellani, “Wind turbine power curve upgrades: part II,” Energies, vol. 12, no. 8, p. 1503, 2019. View at: Publisher Site  Google Scholar
 S. Wan, J. Wan, and C. Zhang, “Comprehensive evaluation of wind power unit performance evaluation based on grey theory and variable weight fuzzy mathematics,” Acta Energiae Solaris Sinica, vol. 36, no. 9, pp. 2285–2291, 2015. View at: Google Scholar
 K. Bassett, R. Carriveau, and D. S.K. Ting, “Vibration analysis of 2.3 MW wind turbine operation using the discrete wavelet transform,” Wind Engineering, vol. 34, no. 4, pp. 375–388, 2010. View at: Publisher Site  Google Scholar
 Y. Chunyu, G. Jianying, and X. Shiguang, “Failure mechanism analysis and failure number prediction of wind turbine blades,” TELKOMNIKA (Telecommunication Computing Electronics and Control), vol. 12, no. 3, pp. 533–540, 2014. View at: Publisher Site  Google Scholar
 S. Chakkor, M. Baghouri, and A. Hajraoui, “Wind turbine fault detection system in real time remote monitoring,” International Journal of Electrical and Computer Engineering (IJECE), vol. 4, no. 6, 2014. View at: Publisher Site  Google Scholar
 J. J. Eggers, R. Bauml, R. Tzschoppe, and B. Girod, “Scalar costa scheme for information embedding,” IEEE Transactions on Signal Processing, vol. 51, no. 4, pp. 1003–1019, 2003. View at: Publisher Site  Google Scholar
 C. Zhang, “To evaluate the performance of Waanje, wind turbine fuzzy comprehensive evaluation based on Grey theory and variable weight,” Journal of Solar Energy, vol. 36, no. 9, pp. 2285–2291, 2015. View at: Google Scholar
 X. Zhang, S. Hao, J. Li et al., “Grey correlation based analysis on impacting factors of maximum power point tracking control of wind power generating unit,” Power System Technology, vol. 39, no. 2, pp. 445–449, 2015. View at: Google Scholar
 A. Zaher, S. D. J. McArther, D. G. Infield, and Y. Patel, “Online wind turbine fault detection through automated SCADA data analysis,” Wind Energy, vol. 12, no. 6, pp. 574–593, 2009. View at: Publisher Site  Google Scholar
 D. Astolfi, L. Scappaticci, and L. Terzi, “Fault diagnosis of wind turbine gearboxes through temperature and vibration data,” International Journal of Renewable Energy Research (IJRER), vol. 7, no. 2, pp. 965–976, 2017. View at: Google Scholar
 Y. Feng, Y. Qiu, C. J. Crabtree, H. Long, and P. J. Tavner, “Monitoring wind turbine gearboxes,” Wind Energy, vol. 16, no. 5, pp. 728–740, 2013. View at: Publisher Site  Google Scholar
 W. Liao, R. Zhang, W. Yu et al., “PV output power forecast based on similar samples and PCA,” Solar Journal, vol. 37, no. 9, pp. 2377–2385, 2016. View at: Google Scholar
 L. Breiman, “Random forests,” Machine Learning, vol. 45, no. 1, pp. 5–32, 2001. View at: Google Scholar
 X. Zhang, H. Sipeng, J. Li et al., “Analysis of influence factors of MPPT control based on grey correlation degree,” Power Grid Technology, vol. 39, no. 2, pp. 445–449, 2015. View at: Google Scholar
 Wind Turbineg Ruiming, Support Vector Machine Analysis Theory and Application, China Power Press, Beijing, China, 2007.
 Y.Y. Hong and Y. Y. Chen, “Placement of power quality monitors using enhanced genetic algorithm and wavelet transform,” IET Generation, Transmission & Distribution, vol. 5, no. 4, pp. 461–466, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Liang Tao 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.