Novel ConditionMonitoringMethod forWind Turbines Based on the Adaptive Multivariate Control Charts and SCADA Data

A novel condition monitoring method based on the adaptive multivariate control charts and the supervisory control and data acquisition (SCADA) system is developed. Two types of control charts are adopted: one is the adaptive exponential weighted moving average (AEWMA) control chart for abnormal state detection, and the other is the multivariate exponential weighted moving average (MEWMA) control chart for anomaly location determination. Optimization procedures for these control charts are implemented to achieve minimum out-of-control average running length. Multivariate regression analysis is utilized to obtain the normal condition prediction model of wind turbine with fault-free SCADA data. After comparing the regression accuracy of several popular algorithms in the MRA, the random forest is adopted for feature selection and regression prediction. Various tests on the wind turbine with normal and abnormal states are conducted. (e performance and robustness of various control charts are compared comprehensively. Compared with conventional control charts, the AEWMA control chart is more sensitive to the abnormal state and thus has a more effective anomaly identification ability and better robustness. It is shown that the MEWMA control chart combined with the out-of-limit number index can effectively locate and identify the abnormal component.


Introduction
With the increasing sustainable energy and environmental demands, wind energy has become one of the world's fastest growing renewable and green energy sources. Due to unstable and unpredictable wind speed characteristics and energy potentials, which are very sensitive to variations in topography and weather patterns, the cost ratios of the operation and maintenance (O&M) costs over the total energy costs per unit output electrical energy from wind turbine systems are considerately high, which is up to 20%∼ 25% [1]. Such high ratios of O&M costs may impede the applications of wind turbine systems compared to other renewable energy systems, such as solar photovoltaic or hot water systems. Consequently, effective condition monitoring (CM) methods for wind turbines are essential for maintenance decisions, which aim to reduce O&M cost [2]. Various signals, such as vibration [3], acoustic emission [4], and motor current [5,6], have been utilized for wind turbine CM systems. However, these approaches require the installation of additional sensors and data acquisition devices, which increase the capital cost and wiring complexity of wind turbine systems. Supervisory control and data acquisition (SCADA) systems have been installed in most modern wind turbines to monitor operational performances.
Currently, the SCADA signal has received a lot of attention owing to its application in wind speed-power forecasting [7][8][9], wind power assessment [10,11], and wind farm performance analysis [12]. A typical SCADA system records comprehensive wind turbine condition parameters, including temperatures (e.g., bearing temperature and oil temperature), wind parameters (e.g., wind speed and wind direction), and energy conversion parameters (e.g., output power, pitch angle, and rotor speed), which would be fault informative. Since no additional sensors or data acquisition devices are needed, the wind turbine CM method based on SCADA data is a cost-effective approach to improve the reliability of wind turbines [13].
Building a model to predict the normal behavior of SCADA parameters is the first issue of the wind turbine CM system. By using advanced SCADA data mining methods, various normal condition prediction models (NCPMs) have been developed to detect the significant changes in wind turbine behavior prior to anomaly occurrences. Kusiak et al. [14][15][16] first employed various data mining algorithms to construct NCPMs for wind turbine anomalies. After detailed comparisons based on the SCADA data collected at a large wind farm, they found that the random forest (RF) algorithm models provided the best accuracy [14]. Gill et al. [17] developed a probabilistic model of a power curve for CM purposes based on copula statistics. Its practical use was demonstrated on the SCADA data taken from a fleet of operational wind turbines. e adaptive neurofuzzy interference systems [18,19] and neural networks [20] have also been adopted to develop various NCPMs. Wang et al. [21] proposed a new NCPM based on heterogeneous signals and information collected from the SCADA system. A linear mixture self-organizing map classifier was applied to differentiate abnormal types. After simulations were carried out on the actual data from a wind farm in north China, the proposed technique was shown to be effective for abnormality detection and prediction. Recently, the Bayesian framework [22], spatiotemporal pattern network [23], and mathematical optimization models [24] were introduced for the early and unsupervised fault diagnosis of wind turbines using SCADA data.
For a given NCPM, the relationship between the input and output SCADA state variables of the wind turbine could be learned. Subsequently, the departure of the current turbine state from the NCPM could be measured online and yield a time series of residuals. e control chart from statistical process control is a time-honored tool to monitor the residuals [25]. If the residuals are statistically different from a normal (or fault-free) reference, the process is considered out of control, and an alarm would be raised accordingly. In recent years, the NCPM combined with control charts has been increasingly used in wind turbine CM systems. Most studies [26][27][28][29][30] used the Shewhart-type control charts, which have been proven to be very effective for detecting greater shifts [25]. However, they are slow in reacting to small and moderate shifts in the mean process. In that regard, the exponential weighted moving average (EWMA) control chart was developed to provide more sensitivity to small mean shifts [25]. Cambron et al. [31][32][33] first applied the EWMA control chart for the CM of wind turbines. Using several applications on the actual SCADA data, the results showed that a shift of 3.4% in annual energy production over a period of 5 years could be detected in time to plan proper maintenance. Helbing and Ritter [34] explained a straightforward method to incorporate nonconstant variance to construct a flexible EWMA control chart. Simulations showed that the FEWMA has lower false alarm rate than the EWMA. Wang et al. [35] deployed the EWMA control chart to derive the criteria for detecting the oil temperature shifts of wind turbine gearboxes. Yang et al. [36] proposed an approach combining data mining and control charts for fault detection in actual wind turbines. Both EWMA and multivariate EWMA (MEWMA) control charts were constructed for comparisons. eir observations showed that the MEWMA is more suitable for early detection and avoidance of errors.
Although the EWMA control chart can provide greater sensitivity to small shifts, it is not as effective as the Shewhart chart, where the shifts in the process mean level are relatively large due to the inertia problem [37]. In actual applications, such as monitoring of wind turbines, the shift of the residuals from the NCPM is unknown, which might cause the insufficiency of the EWMA control chart if the larger shift appears. To overcome the inertia problem, Capizzi and Masarotto [38] first presented an adaptive EWMA (AEWMA) by adaptively adjusting the weight on past observations according to a function of the prediction error. Later, Shu [39] extended the idea of the AEWMA chart on monitoring process locations to the case of monitoring process dispersion. e AEWMA chart is a smooth combination of the Shewhart and EWMA charts; thus, it can reduce the inertia problem. Using the examples on capsule weights and simulated data, both Capizzi and Masarotto [38] and Shu [39] showed that the AEWMA control chart is able to offer an overall good detection performance against shifts of different sizes. However, in the CM of wind turbines, the residual data would be more complicated, and the possibility of the AEWMA control chart holding a better performance than the EWMA chart is still unknown. To the authors' knowledge, the AEWMA control chart has not been used in the CM of wind turbines in the open literature.
In actual engineering, it is not only expected to alarm an abnormal state as early as possible, but determining the cause and location of the abnormal state is also expected. Since the SCADA system records condition parameters of the main components of wind turbines (e.g., the blade, gearbox, main bearing, and generator), the components with the abnormal state might be identified by modeling the control charts of these multivariate conditional parameters. Lately, Yang et al. [36] used the MEWMA to determine which components are likely to contribute to the fault. eir results showed that the MEWMA has a good potential in locating anomaly. e limitation of Yang et al.'s study [36] is that only specified values of MEWMA parameters were tested, indicating that the presented MEWMA might not be the optimal control chart. e optimal design of MEWMA should be conducted to fully realize the potential of the MEWMA in the CM of wind turbines.
A literature review indicated that only a few studies have used the multivariate control charts for the CM of wind turbines; this is particularly true for the abnormal state alarm of wind turbine using adaptive control charts. Moreover, there have been few attempts to comprehensively compare the performance and robustness of both EWMA and AEWMA control charts in monitoring the residuals from the NCPM of wind turbine SCADA data. erefore, the novelty and contributions of this study can be summarized as follows: (i) e framework for the CM of wind turbines is introduced based on the adaptive multivariate control charts (AMCCs). Two AMCCs (AEWMA and MEWMA) are introduced for abnormal state alarm and anomaly location of wind turbines, respectively. An optimal design is conducted to ensure that the obtained control charts are in the optimal state. e remainder of this paper is organized as follows. Section 2 introduces the proposed control charts. Section 3 provides the optimal design procedures. Section 4 describes feature selection and regression prediction on the SCADA data acquired from an operating wind turbine. Section 5 presents the flowchart of the AMCC-based CM method. Section 6 provides several CM examples and discusses the results. Finally, Section 7 lists the conclusions of the study.

AMCCs
Two AMCCs (AEWMA and MEWMA) are introduced for abnormal state alarm and anomaly location of wind turbines, respectively. e structures and procedures for these two control charts are derived in this section.

Abnormal State Alarm.
Monitoring data that obey the same distribution are represented by X i,1 , X i,1 , . . . , X i,1 , where i � 1, 2, . . . is the sampling time and n is the size of each sample. e mean and variance of the data are denoted by μ 0 and σ 2 0 , respectively. When the process is out of control, the mean of the data becomes μ 1 � μ 0 + δσ 0 , in which δ is the shift parameter. We define the mean of the sample data as X i � (1/n) n j�1 X ij , and thus, the EWMA statistics for monitoring mean shift of the sample data could be written as follows: where λ is the smoothing parameter, and 0 < λ < 1. Without loss of generality, we can let Y 0�0 . Lucas and Saccucci [40] pointed out that for smaller value of λ, the EWMA statics can detect a smaller mean shift faster. When takes λ a greater value, the EWMA statistic would have an accurate sensitivity to the larger mean shift. eoretically, the EWMA control charts can be customized to detect specific shifts in the process. However, for the actual wind turbine monitoring data, the mean shift is usually fluctuated in a certain range. e designed value of λ makes it difficult to adapt to the change in the actual mean shift. To overcome this inertia problem, the AEWMA statistic is proposed by [38] where e i � X i − Y i is the error term and ϕ(·) represents the score function. Note that for e i ≠ 0, and the AEWMA statistic can be rewritten as where w(e i ) � ((ϕ(e i ))/e i ) is the equivalent smoothing parameter. Evidently, the AEWMA statistic can adaptively adjust the weight of the estimate value at the previous time (Y i − 1) according to the prediction error at the current time. us, it can balance the requirements of various mean shifts to the smoothing parameters. Yashchin [37] suggested the Huber function as the score function, and its expression is given by where c is the error limit. e static Y i also obeys the same distribution with X i and has the same mean value with X i . When the sampling size n is large enough, the variance of Y i can be expressed as σ 2 Y � σ 2 X , where σ 2 X is the variance of X i , leaving us with σ 2 X � (σ 2 0 /n). erefore, the upper control limit (UCL) and lower control limit (LCL) of the AEWMA control chart could be expressed as follows: where k is the control limit parameter. From equations (3)-(5), it can be discovered that three parameters: λ, c, and k should be determined to obtain the control limits of AEWMA control charts. e determination of these parameters will be discussed in the following section. It is observed that for c ⟶ ∞, we have ϕ(e i ) � λe i and w(e i ) � λ. In this case, the AEWMA statistic degenerates into the EWMA statistic, and its control limits can be expressed as ,

Anomaly Location.
In the abnormal state alarm of wind turbines, the data monitored by the AEWMA control chart are univariate, i.e., the output power data of the wind turbine. In addition to the early warning of an abnormal state, we also expect this method to identify the cause and location of the anomaly state. Fortunately, the SCADA system records condition parameters of the main components of wind turbines (e.g., the blade, gearbox, main bearing, and generator). us, we introduce the MEWMA control charts to monitor these multivariate conditional parameters, and then the components with anomaly state might be identified. From the univariate EWMA control chart, Lowry et al. [41] proposed the MEWMA control chart, and its statistic can be expressed as Shock and Vibration 3 where Y i and X i are the p dimensional multivariate data vectors. We assume that Y 0 � 0. r denotes the smoothing parameter, leaving us with 0 < r < 1. e MEWMA control chart will sound an alarm if the following conditions are satisfied: (8) in which H is the given control limit and Y i denotes the covariance matrix of the Y i . Hence, we have where X i denotes the covariance matrix of the X i . To evaluate the contribution of different dimensional data to the MEWMA statistic, following variables are defined as in which j � 1, 2, . . . , p. e larger variation of Q ij indicates that the contribution of the j th dimensional variable to the MEWMA statistic is significant. e component corresponding to this variable is more likely to be in the abnormal state. erefore, the MEWMA control chart is adopted to monitor the multidimensional SCADA data and also to identify and locate the anomaly component of the wind turbine by analyzing the variation of Q ij .

Optimal Design of Control Charts
e average run length (ARL), which refers to the average number of extracted samples from the beginning of the control to the emission of alarm, is used to measure the performance of various control charts. Here, the ARL 0 is used as the in-control ARL and the ARL 1 as the out-ofcontrol ARL. Typically, ARL 0 is desired to be as large as possible and ARL 1 to be as small as possible. Hence, the designed control chart can raise an alarm on the existence of abnormal deviations as soon as possible on the premise of a lower false alarm rate. By satisfying the goal of minimizing the ARL 1 under certain ARL 0 , the optimal design procedures for the AEWMA and MEWMA control charts are presented, respectively, in the following subsections.

AEWMA.
ree parameters, including λ, c, and k, should be determined for the AEWMA control charts. Clearly, the selection of λ and c plays a key role in the performance of AEWMA control charts. Generally, the lower value of λ or greater value of c should be selected for small mean shift, while the greater value of λ or lower value of c would be favorable for detecting large mean shift. erefore, the design of AEWMA control charts is a multiobjective optimization problem. Capizzi and Masarotto [38] utilized the simulated annealing algorithm (SAA) for the parameter optimization of AEWMA control charts. However, the requirement for the initial value of SAA is relatively high. Once the initial value significantly deviates from the optimal value, it is difficult to converge to the optimal value. To improve the convergence speed of SAA, Shu [39] proposed a "two-step method." First, the AEWMA control chart is treated as a conventional EWMA control chart, and the optimal value of λ is obtained under certain ARL 0 using SAA. en, on the premise of given value of λ, the value of c is optimized. Figure 1 shows the flowchart for the optimal design of the AEWMA control chart. Detailed procedures are described as follows: (1) Sample size n and in-control ARL ARL 0 are selected.
(2) Typically, the range of parameter optimization is selected as λ < 0.8, 3 < λ < 7, and k < 10. (3) By setting c ⟶ ∞ (e.g., c � 1000), the AEWMA control chart is degenerated to the EWMA control chart. e parameter of the EWMA under the shift δ 1 is denoted by θ 1 � (λ, k). e optimal θ 1 should satisfy the following optimization problem: (4) A small positive number α (taken as α � 0.05 in this study) is selected to ensure that the control chart will not lose too much accuracy after the introduction of c. Based on the optimal parameter λ obtained in step (3), the optimal parameter θ 2 � (c, k) of the AEWMA control charts with mean shift δ 2 could be obtained by solving the following optimization problem: In the above steps, the calculation of ARL can be obtained using the Monte Carlo sampling method.

MEWMA.
e MEWMA control chart has two parameters: r and H. Similar to the AEWMA control chart, small values of the smoothing parameter r should be selected for small mean shifts, while large r values have advantages in detecting large mean shift. Runger and Prabhu [42] proposed a Markov chain algorithm (MCA) for designing a MEWMA control chart. For i, j � 0, 1, 2, . . . , m 1 , the transition probability from state i to state j is denoted by p(i, j), and its definition is given by 4 Shock and Vibration where χ 2 (p, c) is the noncentral chi-square random variable, p is the number of dimensions, c � (((1 − r)ig 1 )/r) 2 is noncentral parameter, and g 1 � (2/(2m 1 + 1)) ����������� (r/(r − 2))H. Based on the transition probability, the (m 1 + 1) dimension transition matrix P 0 could be constructed. us, the ARL 0 of MEWMA control charts could be calculated by where s is the initial probability vector, I is the unit vector, and 1 is a vector with all of its elements equal to 1. Similarly, the ARL 1 of the MEWMA control chart could also be obtained.
Based on the MCA [42], we use the partition method to obtain the optimal parameters of MEWMA control charts. e partition method generates a combination of a smoothing parameter r and a control limit H, satisfying a given ARL 0 , and finding the optimal smoothing parameter. Figure 2 presents the flowchart for the optimal design of MEWMA control charts. Detailed procedures are described as follows: (1) For a fixed smoothing parameter r, the method inspects the middle point of a lower control limit H low and a upper control limit H up such that ARL H low ≤ ARL 0 and ARL H up ≥ ARL 0 . (2) Once H mid , the middle point of two control limits is obtained, and ARL can be calculated by using the MCA. If the difference between ARL 0 and the newly computed ARL is less than a small number (i.e., ϵ < 10 − 3 ), the smoothing parameter r and the control limit H mid is a pair that can satisfy the given ARL 0 . Otherwise, keep following the previous procedures until a sought pair is found. (3) If this task is carried out until the method covers a whole range of smoothing parameter (0 < r ≤ 1), a number of combinations of r and H can be obtained.
With the combinations obtained, ARL 1 values can be calculated for a given shift δ. en, the smoothing parameter r for which ARL 1 is the smallest can be identified.

MRA on Fault-Free SCADA Data
In previous sections, both control charts of AEWMA and MEWMA have been introduced for the abnormal state alarm and anomaly location of wind turbines. e optimal design procedures for these control charts have been presented. e residuals monitored by these control charts are yielded by the departure of real-time SCADA data from the predictions of NCPM. In this section, we utilize the MRA to construct the NCPM of wind turbines with fault-free SCADA data. Several popular algorithms in MRA, including the RF, LASSO, and RFE, are used for feature selection and regression prediction.

Data Descriptions.
is study aims to monitor and diagnose doubly fed wind turbines with rated power of 2 MW. Typically, the SCADA data of the unit include output power, speed, torque, temperature, and pitch angle. e data record interval is 10 min. To correctly establish the NCPM of wind turbines, the anomaly data should be avoided as much as possible. By reading the record table of the SCADA system, it was found that no anomaly was reported in the time period from 12/26/2013 to 2/12/2014. e wind turbine unit was built and connected to the grid in early 2012. In this time period, the unit has passed the initial running stage and is in the stage of normal power generation. erefore, the data segment is ideal for MRA to construct the NCPM of wind turbines.
Shock and Vibration 5 points during the maintenance downtime, the total amount of data is 6135 points.

Feature Selection and Regression Prediction.
At the beginning of MRA on the fault-free SCADA data, to minimize the problem of model deviation due to the lack of important variables, we usually select as many argument variables as possible. In this study, we select the output power as the response variable and the remaining 44 variables as argument variables. However, in the process of actual modeling, it is necessary to select a variable subset (feature selection) which has the best ability to explain the response variable to improve the regression and prediction accuracy of the NCPM ( [43,44]). Before feature selection, the raw SCADA data should be standardized as follows: where n is the sample size and y i and y i are predicted and actual values of output power, respectively. e RF, LASSO, and RFE, which are popular algorithms in MRA, are utilized for feature selection and regression prediction. Basic ideas and characteristics for these algorithms are introduced. e RF is an integrated machine learning method [45]. It employs random resampling technology bootstrap and node random splitting technology to construct multiple decision trees, with the final classification results are obtained by voting. e RF has the ability to analyze the classification characteristics of complex interactions. It has a fine robustness for noise data and a faster learning speed. Its variable importance measure can be used as a feature selection tool for high dimensional data. e core algorithm uses the RF package in R software, in which the parameter takes the value of � n √ recommended by Breiman [45] (n is the number of features of the training data set). e number of trees is set to be n tree � 500. Out-of-control ARL for given shi ARL 1 Identification of optimal parameters r, H LASSO [46] is a linear model for estimating sparse parameters, especially for reducing the number of parameters. is method uses the l 1 norm to compress the coefficient of the model and directly makes the values of insignificant model parameters smaller (including zero).
is gives the LASSO the advantage of feature selection and ridge regression. Without changing the accuracy of the model test set, the dimension of the feature could be effectively reduced by using the LASSO regression model. e core algorithm adopts the LARS package in R software, and the cross validation is utilized to select the penalty parameter that controls the sparse parameter estimation. e main idea of RFE [47] is to build the model iteratively and then select the best (or worst) feature (which can be selected according to the value of coefficients). e iteration process on the remaining feature will be conducted until all the features have been traversed. e stability of RFE mostly depends on the type of iteration model.
In this study, the core algorithm is implemented through the CARET package in R software. After a series of tests, the decision tree model (treebagFuncs) is selected as the iteration model.
For the fault-free wind power SCADA data, the above three algorithms are used for regression prediction and feature selection. Table 1 displays the regression accuracy in the metrics of RMSE, MAPE, and MAE. After comparisons have been conducted, it can be observed that the RF has the best accuracy in regression. us, the feature selection based on the RF is carried out, and the top 15 features are shown in Table 2. According to the common sense, the parameters closely related to the output power of the wind turbine are wind speed, generator speed, torque, rotor speed, etc. ese features have been reflected in the feature ranking of RF. In addition, the rankings of generator phase current and phase voltage, as well as the parameters of several temperature measuring points (including the gearbox, bearing, generator, and even nacelle), were relatively higher. ese parameters are not easy to judge and select directly through common sense. Figure 3 also presents the comparison between the source SCADA data and regression prediction results. For the sake of simplicity, only four argument variables (i.e., the rotation torque, generator current, average wind speed in 10 min, and generator speed) on the response variable (output power of the wind turbine unit) are illustrated in the figure.

Wind Turbine CM System Based on AMCCs
Some key contents, including the structures of AMCCs, optimal design procedures of these control charts, and construction of NCPM with fault-free SCADA data, have been introduced in the previous sections, respectively.
How to implement these core algorithms needs to be explained for engineering applications. Figure 4 presents the flowchart for the wind turbine CM system based on AMCCs.
e entire process could be summarized as follows: (1) MRA is utilized to construct the NCPM of wind turbines with fault-free SCADA data. In this study, the RF shows better performance in feature selection and regression prediction. (2) Time-variable residuals of output power are produced by measuring the difference between the realtime SCADA data and the predictions of NCPM. (3) For the goal of minimum out-of-control ARL (see Figure 1), the optimal AEWMA control chart is constructed to monitor the output power residuals.
Steps (2) and (3) will be continued until the abnormal state is alarmed. (4) e optimal MEWMA control chart is established (see Figure 2) to model conditional parameters of main components, which are acquired from realtime SCADA data. e component under an anomaly state could then be located.
In the following, the effectiveness of the proposed CM method is shown by several examples. e performance and robustness of various control charts are compared in detail.

CM Examples
Based on the feature selection and regression prediction results, CM practice on the wind turbine unit is carried out. During the period from 12/1/2015 to 6/1/2016, there were three anomalies, namely, the generator brush worn, gearbox running hot in low generator stage, and shaft bearing overtemperature. e specific time of alarm log is shown in Table 3. For each anomaly, the number of monitored data points is 500, and the exact anomaly data point is also given in the table for comparisons.  Generator phase A current 3 Average wind speed in 10 min 4 Generator speed 5 Rotor speed 6 Blade yaw angle 7 Generator temperature 8 Gearbox bearing temperature 9 Nacelle angle 10 Generator phase A voltage 11 Nacelle revolution 12 Gearbox temperature 13 Ambient temperature 14 Nacelle temperature 15 Bearing temperature 6.1. Abnormal State Alarm. By using the NCPM model obtained in the previous section, the output power of the unit before and after the fault (500 data points in Table 3) is predicted, and then the residual is obtained by measuring the difference in the actual output power. e mean and variance of the predicted residuals for three fault data are all less than 0.05 and 0.08. Given ARL 0 � 500 and shift range (0.4-4), the optimal parameters of AEWMA control chart are then obtained, as shown in Table 4. For comparison, the parameters of the optimal EWMA control charts corresponding to different shifts (0.4, 2, 4) are also given in the table.
As mentioned before, the out-of-control ARL is an important index to evaluate the performance of control charts. Figure 5 shows the variation of ARL 1 with mean shift in the range of 0-4 for the designed control charts in Table 4. Evidently, when the shift is zero, the out-of-control ARL is equal to the in-control ARL, i.e., ARL 1 � ARL 0 . With an increase in a shift, the value of ARL 1 gradually decreases. Under small shifts (δ < 1.5), the value of ARL 1 for the AEWMA is lower than that of the EWMA control charts, especially for the EWMA control chart with larger smoothing parameters (EWMA-2 and EWMA-3). is means that the AEWMA behaves more sensitively and could give warnings of abnormal states earlier than the EWMA control charts.
When the shift becomes large enough (δ > 2), the difference in ARL 1 between the AEWMA and EWMA control charts is not significant, indicating that under large shift, the AEWMA can still maintain a performance comparable to the EWMA control charts. is is consistent with the theoretical expectation of the AEWMA. e AEWMA control charts are established for the output power residuals with anomaly A, B, and C, as shown in Figure 6. Figure 7 also presents the residuals monitored by various EWMA control charts for comparisons. It can be  7 h), respectively. us, one can say that the alarm time of AEWMA control charts could be several hours ahead of the SCADA system, and in this study, the maximum promptness appears in anomaly A (about 5.5 h). Compared with the AEWMA control chart, the EWMA control charts behave less sensitively to fault and have poor robustness. For the EWMA-3 of anomaly A (see Figure 7(a)), and EWMA-1, EWMA-2, and EWMA-3 of anomaly B (see Figure 7(b)), the abnormal state is not identified during the monitoring period. Although the faults    are signaled earlier than the SCADA system, the alarm time still lags behind the AEWMA control chart among the rest of EWMA control charts. For anomaly A, the EWMA-2 control chart (see Figure 7 From the above CM examples, one can say that compared with the EWMA control charts, the AEWMA control chart behaves more sensitively to the abnormal state. us, it can effectively identify the abnormal state and has better robustness. is is of great application value to the CM of practical wind turbine units.

Anomaly Location.
In the previous section, it is demonstrated that the AEWMA control chart can effectively identify the abnormal state. However, for complex electromechanical systems (i.e., the wind turbine), in addition to the early warning of abnormal state, it is also expected to identify the anomaly component, which is called the  Table 2, in addition to environmental features (such as the average wind speed and ambient temperature), there are also features characterizing the working conditions of main components, including the generator speed, gearbox temperature, and blade yaw angle.
is section uses the MEWMA control chart to model multidimensional data, studies the influence of various features on the monitoring statistics, and realizes the effective location and identification of the anomaly. e input parameter for MEWMA control chart should be determined by δ � μ T X i μ, where μ represents the mean vector of multidimensional data. After calculation, it is discovered that the values of δ for three samples are all lower than 4.5. us, by setting δ � 4.5, data dimension p � 15, and ARL 0 � 500, the optimal parameters of MEWMA control chart can be obtained as r � 0.48 and H � 29.65. For the 15dimensional monitoring data containing the anomaly A, B, and C, the MEWMA statistic could be evaluated. In Section 2.2, we have defined Q ij (see equation (10)) to evaluate the contribution of various dimensional data to the MEWMA statistic.
When the different dimension variables are excluded, the changes in Q ij are shown in Figure 8, respectively. As  can be identified from the figures, some specific dimension variables contribute greatly to the MEWMA statistic, while others seem to have little influence on the statistic. Nevertheless, from Figure 8, it is not easy to directly identify which dimension variables have remarkable contribution on Q ij . We can define the number (or frequencies) of MEWMA statistic beyond the control limit as an metric. en, the out-of-limit number (OLN) variation of various dimensional data is shown in Figure 9, respectively. Some observations could be summarized as follows: (1) For anomaly A (see Figure 9(a)), the variations of OLN after the removal of variable 7, 1, and 2 are greater than the removal of other variables. From Table 2, one can see that the variable 7 represents "generator temperature," the variable 1 "rotation torque," and variable 2 "generator phase A current." Consequently, it is estimated that the generator is more likely to have an anomaly. e alarm log of the SCADA system confirms that anomaly A does appear in the generator, which is described as the generator brush worn in Table 3. (2) For anomaly B, as shown in Figure 9(b), the maximum OLN variation occurs at variable 12 ("gearbox temperature"), indicating that the gearbox might be in an abnormal state. is is also consistent with the anomaly description "gearbox running hot in low generator stage" (see Table 3). (3) When anomaly C is considered (see Figure 9(c)), one can find that the maximum OLN variation appears at   Shock and Vibration 13 variable 15 ("bearing temperature"), and the main bearing is more likely to be in an abnormal state. is agrees well with the description of anomaly C "shaft bearing overtemperature" (see Table 3).
rough the accurate location of the three different anomalies, one can see that the MEWMA control chart combined with the OLN index can effectively locate and identify the abnormal component.

Conclusions
A novel CM method of wind turbines is introduced based on AMCCs and SCADA data. Two AMCCs (AEWMA and MEWMA) are proposed for abnormal state alarm and anomaly location of wind turbines, respectively. Optimization procedures for these control charts are implemented with the goal of minimum out-of-control ARL. MRA is utilized to obtain the NCPM of wind turbine with fault-free SCADA data. After conducting comparisons of the regression accuracy of several popular algorithms in the MRA, the RF is used for feature selection and regression prediction. Various tests on a wind turbine with normal and abnormal states are conducted. e performance and robustness of various control charts are compared comprehensively. Compared with the EWMA control charts, the AEWMA control chart behaves more sensitively to the abnormal state and thus has a more effective anomaly identification ability and better robustness. By accurately locating three different anomalies, it is demonstrated that the MEWMA control chart combined with the OLN index can effectively locate and identify the abnormal component.

Data Availability
e wind turbine data used to support the findings of this study were supplied by a wind power plant under license and so cannot be made freely available. Requests for access to these data should be made to (Qinkai Han, Email: han-qinkai@hotmail.com).  Variable no. 14 Shock and Vibration