A CBA-KELM-Based Recognition Method for Fault Diagnosis of Wind Turbines with Time-Domain Analysis and Multisensor Data Fusion

. Fault diagnosis technology (FDT) is an eﬀective tool to ensure stability and reliable operation in wind turbines. In this paper, a novel fault diagnosis methodology based on a cloud bat algorithm (CBA)-kernel extreme learning machines (KELM) approach for wind turbines is proposed via combination of the multisensor data fusion technique and time-domain analysis. First, the derived method calculates the time-domain indices of raw signals, and the fused time-domain indexes dataset are obtained by the multisensor data fusion. Then, the CBA-based KELM recognition model that can identify fault patterns of a wind turbine gearbox (WTB) is automatically established with the fused dataset. The dataset includes a large number of samples involving 6 fault types under diﬀerent operational conditions by 5 accelerometers. The eﬀectiveness and feasibility of this proposed method are proved by adopting the datasets originated from the test rig, and it achieves a diagnostic accuracy of 96.25%. Finally, compared with the other peer-to-peer methods, the experimental classiﬁcation results show that the proposed CBA-KELM technique has the best performances.


Introduction
e wind turbines [1,2], usually working under rough and complex environmental conditions, are very vulnerable to failure throughout their lifetime [3][4][5].Consequently, condition monitoring (CM) and fault diagnosis technology (FDT) for wind turbines are of vital significance not only to ensure the unit in safety running but also to reduce the operation and maintenance (O&M) costs [6].However, the varying loads and changing environmental conditions make the CM of the wind turbines a difficult task.In order to solve these problems, many in-depth and detailed studies have been carried out recently [7][8][9][10][11].It can be divided into two different kinds of classifications including qualitative diagnosis and quantitative diagnosis methods.
Among the qualitative experience approaches, in [12], it provides a method using RMS and extreme values of vibration signals for the CM of the wind turbine gearbox (WTB), but the statistical features of vibrations that characterize key condition indicators for WTB have no uniform standard.In [13], a fault diagnosis expert system of a wind turbine is established, but this method is highly dependent on expert experience.To address the aforementioned issues, in [14], a fault diagnosis technology based on the expert system and fault tree analysis is proposed, but the drawback of this approach is that it only applies to a simple wind turbine system.In [15], a signed directed graph method is adopted for fault diagnosis of wind turbines.However, it has been verified and approved only for systems with fewer loops.
Quantitative fault diagnosis, as the most common method for wind turbine fault diagnosis, includes the analytical model-based methods and the data-driven methods which are the main analytical methods.e method based on signal processing has been presented in [8,11,[16][17][18][19][20].In [8], a complex wavelet transform is proposed for multifault detection using vibration signals collected from a real wind turbine.However, the time-varying characteristics of periodic signals are not fully considered and utilized.In [11], the Hilbert transform and ensemble empirical mode decomposition method are applied to carry out fault analysis and diagnosis of a direct-drive wind turbine.However, the high computational costs and the noise pollution in the signal reconstruction by adding white noise increase the difficulty of fault diagnosis.A single senor-based blind source separation method presented in [16] has been utilized for fault diagnosis of the WTB.However, it neglects the effect of the time epoch number on the diagnostic results.In [17], bispectrum and hidden Markov model technology are applied to detect wind turbine faults, but it is shown that the complexity of this algorithm is too high.
In recent years, as the intelligent fault diagnosis algorithms, the machine learning technique has achieved enormous success.In [21,22], the artificial neural network approach is applied to fault diagnosis.Nevertheless, the success of this method highly relies on the selection of initial threshold, weight values, and the database.Wang and Liu [23] developed a fault diagnosis strategy for wind turbine brake systems using the evidence theory, wavelet theory, and fuzzy neural network.However, how to measure the information's uncertainty affects the output of results.Bechkaoui et al. [24] and Merabet et al. [25] present the monitoring strategy of wind turbines by the fuzzy logic technique, but the accuracy of the complex model simulation system seriously affects the diagnostic effect.In [26], the extreme learning machine technique is used for fault classification of wind turbine equipment.In order to increase the accuracy of prediction, a kernel extreme learning machine (KELM) is proposed in [27][28][29].KELM has a number of unique advantages, such as high learning speed, more generalized performance, and high robustness.However, it is sensitive to the parameters such as a kernel parameter and cost parameter, which may result in a low diagnostic rate.
In this paper, with time-domain analysis and multisensor data fusion technology, a variant of the basic KELM, called the cloud bat algorithm-(CBA-) based KELM, is applied for the first time to fault diagnosis of the WTB.e structure of the paper is presented as follows.Section 2 succinctly introduces the wind turbine drive train and gearbox.In Section 3, the time-domain analysis, the multisensor data fusion, and the improved KELM technology are described, and the implementation process of the proposed failure diagnosis technology for WTB is also detailed in this section.In Section 4, the validity of the proposed approach is verified using the dataset collected from the test bench.Diagnostic accuracies of the proposed approach are also compared with that of the other peer-to-peer methods.Finally, conclusions are summarized in Section 5.

Wind Turbine Drive Train and Gearbox
2.1.Drive Train of Wind Turbine.At present, wind power has been widely used as one of the most important promising environmentally clean energy sources.
e wind turbine drive train (shown in Figure 1) is the load-bearing component in the process of energy conversion and transmission, including a rotor hub, main bearings, gearbox, and generator.e wind energy captured by the rotor hub with blades is transformed to the mechanical energy, then the main bearings transfers the mechanical energy to the generator via the gearbox that increases the lower rotor speed to a higher speed for the generator.
e generator converts mechanical energy into electrical energy.

Wind Turbine Gearbox.
e gearbox is a key component of the wind turbine.On the one hand, the ever-increasing capacity of wind turbines has caused the gearbox structure to become more and more complex, in the event of failure will directly lose nearly $1 million [30].On the other hand, when it is still in the early stage of failure, due to the influence of environmental noise and other factors, the resulting fault feature components are relatively weak, resulting in the difficulty of fault diagnosis [31,32].
Compared to any other application, the wind turbine gearboxes tend to fail early.Although their life expectancy is 20 years, it is far from meeting the requirement in the actual operation process [33].e construction of the typical WTB which consists of two-stage bevel wheel and one-stage planetary and its transmission diagram is shown in Fig- ure 2. e multispeed transmission is represented by three parts, namely, low, intermediate, and high speed.
Observing the pictures from Figures 3(a)-3(d), the worn tooth, tooth fracture, tooth surface deformation, and bearing faults of the WTB are demonstrated according to the results of the field investigation.Due to long-term operation of WTB under variable loads and variable speed conditions, it will lead to frequent faults and difficulty in diagnosis.In order to meet the actual situation and ensure the accuracy of the diagnosis, different kinds of faults and the same fault under different working conditions are considered in detail, including complex faults which are difficult to diagnose in this research.Shock and Vibration on different types of wind turbines.To improve these conditions, the fault diagnosis of WTB is studied using the CBA-KELM and information fusion technology in this paper.

Multisensor Data Fusion Technology and CBA-KELM-Based Fault Monitoring Technique
After collecting vibration information on the wind turbine gearbox, it can preprocess the original vibration signals by calculating the 13 different kinds of time-domain indicators (TDI), such as where x rms represents the root mean square value, σ 2 denotes the variance, x q represents the kurtosis index, C represents the peak index, I denotes the impulse index, x i represents the mean value, and x max , x min , x p−p , x r , K, L, x ′ , and N denote the maximum value, the minimum value, the peak to peak value, the square root amplitude, the waveform index, the abundance index, the average amplitude, and the length calculation of each index, respectively.

Multisensor Data Fusion Technology.
With the rapid development of many different types of sensors, more and more data can be used for scientific researches.Data fusion technology is an effective way to optimize the use of multisource massive data, and it can achieve inferences that are not feasible from a source or single sensor.ree different kinds of fusion methods, namely, decision-level fusion, feature-level fusion, and data-level fusion, are introduced in detail in [36,37].
As shown in Figure 4, in order to reduce the computational complexity and improve the fault diagnosis rate, multisensor datasets are obtained in combination with feature-level fusion and decision-level fusion technology in parallel superposition in this research.In the feature-level fusion, representative features are extracted from multisensors observation data and combined into a single connection feature vector, which is input into a pattern recognition model.
e decision-level fusion involves merging information at a higher level of abstraction, combining the results of multiple sensor information to generate the final fusion decision.
As shown in Figure 5, the specific process is described as follows: (1) Collect the vibration signals from the different acceleration sensors, and it can be defined by x t , where n denotes the number of sensors: (2) Assuming that the finite length of discrete time series information collected from each sensor is expressed as , where c n i denotes the fault classification of n th sensor and M represents the data length of each sensor.
(3) According to equations (1)-( 5), the original vibration signals are pretreated and analyzed, and it can get the representative feature-level information by calculating the 13 different kinds of TDI.In addition, decision-level information can also be obtained in this step.(4) After the above steps, it obtains the dataset for each sensor, where x 1 (t) denotes the dataset for the first sensor.By analogy, datasets (x t ) of all sensors can be obtained through time-domain statistical analysis: (5) Combining the feature-level fusion and data-level fusion technologies, we can get the fusion dataset χ x in parallel superposition, where χ x ∈ R (M/N)×(13×n+1) , s � M/N.us, it can be defined as 4 Shock and Vibration

Kernel Extreme Learning Machine.
ELM is a kind of new learning algorithm for the single-hidden layer feedforward networks (SLFNs) developed by Huang et al. [38].Compared with the other traditional gradient-descent algorithms such as support vector machine and back propagation (BP) algorithm, the ELM tends to have better scalability and generalization performance.e architecture of an SLFN is shown in Figure 6.
Give N arbitrary distinct samples (x i , t i ), where us, it can be known that In this architecture, where w i � [w i1 , w i2 , . . ., w ik ] T represents the weight vector between the i th hidden node and the input nodes, β i � [β i1 , β i2 , . . ., β im ] T denotes the weight vector between the i th hidden neuron and the output neurons, y j represents the outputs of the network, b i is the threshold of neurons in the i th hidden layer, n, k, and m indicates the number of nodes in the input layer (hidden layer/output layer), respectively, and g(x) denotes the sigmoid function.
e ELM classifier can be expressed as the following equation: With zero error means, the SLFNs with k hidden neurons can approximate these N samples that  k j�1 ‖y j − t j ‖ � 0. us, it can know that there exist b i , w i , and β i , which can be described as en, its compact version can be fully expressed as

Shock and Vibration 5
Hβ � T, (11) where us, the output layer weight vector can be denoted by equation (13) according to the basic theory of ELM, where H * represents the Moore-Penrose generalized inverse matrix, which can be described in equation ( 14) and C is called the cost parameter, related to the stability and generalized performance: Because the feature mapping H do not need to know to users, a kernel-based ELM is applied for improving learning speed, generalized performance, and accuracy of judgment in this paper.e KELM algorithm is an improved algorithm in combination with the kernel function on the proposed ELM algorithm.us, the kernel matrix and the output function of the ELM classifier can be denoted by equations (15) and (16), respectively.And the radial basis function kernel is chosen as follows (equation ( 17)), where c is the kernel parameter:

Kernel Extreme Learning Machine Optimized Cloud Bat
Algorithm.Due to the existence of the kernel function, this algorithm is sensitive to parameter setting.e parameters c and C have a significant influence on the performance of the algorithm when executing the KELM learning algorithm.e CBA algorithm [39] not only has the advantages of the bat algorithm (BA) but also has other advantages such as reduction of the loudness, increasing of the pulse rate emission, and strong stability.erefore, c and C are optimized using the CBA algorithm to create a CBA-optimized KELM prediction model.e CBA algorithm which takes the bat echolocation predation mechanism as the starting point and takes the universality of the normal cloud model as the basis is proposed in this research.Suppose the KELM classification accuracy as pre(C, c).us, within a given range, it needs to be guaranteed to get the maximum classification accuracy as follows: R � max pre(C, c).(18) e proposed algorithm of the CBA-KELM can be described as follows: (1) Randomly initialize the relevant parameters, evaluate the initial positions using the fitness function, and search for a population elite x pbest .( 2 where β is a random digit in the range [0, 1], the pulse frequency range is defined by [f min , f max ], X t i denotes the position of the bat, and V t i represents the velocity of the bat.( 5) Steps ( 2)-( 4) are executed until the predefined maximum number of optimizing generation is satisfied.
Compared with the KELM method, the CBA-KELM is verified more effectively by using the same sample in this section.
e programs are run in the MATLAB R2017a environment, and the results are shown in Figure 7. Obviously, the CBA-KELM has a better performance.

CBA-KELM-Based Fault Monitoring Technique for Wind
Turbines.According to the previous discussion, the 13n parameters are chosen to form the feature vector, including the signals collected from n different sensors.In this paper, the raw signals originated from five sensors in different locations containing six types of faults.e KELM-based SLFN fault diagnosis model is established and illustrated in Figure 8, which has 65 inputs and six outputs.Based on the labeled fused dataset, the SLFN neural network is trained according to the learning algorithm supported by the CBA-KELM algorithm.As shown in Figure 9, the application of the CBA-KELM method for the WTB fault diagnosis can be summarized as follows: Step 1. Collect the vibration signals x t from the sensors mounted on the WTB, and the fault classification is defined according to different faults.
Step 2. In order to obtain the fused dataset, the TDI index of x n (t) is calculated.us, we can get the χ x matrix from the original vibration signal using equations ( 1)-( 6), and the feature level fusion and data level fusion technology are merged in a way of parallel superposition.
Step 3. Normalize the fused datasets.It can obtain the preprocessed data samples to test or train the fault diagnosis model after normalizing the fused datasets.
Step 4. Divide the data samples.e ratio of the testing dataset and training dataset is randomly set.e testing samples are applied to validate the fault diagnosis model, and the training samples are adopted to train this model.
e KELM parameters c and C play an important role in accuracy and generalization, so they are optimized by using the CBA algorithm in this step.
Step 6. Obtain the fault diagnosis model using the CBA-KELM approach.It can judge the fault type when employing the fault diagnosis model.

Experimental Verification
In order to correctly verify the proposed method, the raw vibration signals of the faulty gearbox originated from Qian Peng Company's rotating machinery vibration analysis and fault diagnosis test platform system [40].e wind power gearbox runs in a harsh environment for a long time and works under different loads and rotational speeds.erefore, the experimental device takes these conditions into full consideration and simulates the gearbox faults under different loads and rotating speeds.e collected data reflect the various fault types of the gearbox, including not only single faults such as pitting tooth, worn tooth, and tooth fracture but also comprehensive composite faults that are difficult to diagnose, which is consistent with the actual operation of the gearbox in a wind farm.
e presented CBA-KELM method is then used in the fault diagnosis, and different kinds of methods are also shown for the comparative analysis.

Data Acquisition and Preprocessing.
e test platform system, which has the 5.12 kHz sampling rate, consists of a variable-speed drive motor, bearings, gearbox, shafts, and so on.Various faults can be quickly simulated by adjusting the counterweight, the installation position of the parts, and the organic combination of components.e map diagram of this test platform system is shown in Figure 10.It can complete the collection of on-site vibration signals and has the following functions: signal acquisition, real-time display, signal preservation, etc. e gearbox used for the experiment is a decelerator.An input gear (N Z 1 � 55, modulus � 2) is installed on the input shaft, and an output gear (N Z 2 � 75, modulus � 2) is mounted on the output shaft, where Z 1 and Z 2 represent the pinion gear and the large gear, respectively.And N Z 1 and N Z 2 are the number of teeth of the pinion and the large gear.
In this research, six gearbox failures (i.e., c n i � 1, 2, 3, 4, 5, 6 { }) were simulated at different rotational speeds and loads, including complex faults which are difficult to diagnose.e detailed explanation of the dataset is shown in Table 1.e original vibration signals originate from five accelerometers (sensor 1-5) mounted at different positions on the gearbox, namely, input shaft motor-side bearing, the output shaft motor-side bearing, the input shaft load-side bearing, the output shaft load-side bearing Y, and the output shaft load-side bearing X.In addition to the five accelerometers, a sensor (sensor 6) for measuring the gearbox input shaft speed is installed.e magnetic powder torque device shown on the test bench can automatically adjust the torque.Corresponding to six kinds of fault classification, 6,064,640 vibratory signals are collected.us, the total datasets are obtained from the gearbox fault diagnosis experiments.
Using the proposed approach, the raw signals are preprocessed.According to the time-domain analysis, the initial matrix set is converted to another matrix.e time-domain parameters are calculated in rows of 512.In parallel superposition way, the feature level and data level are merged Shock and Vibration at the same time.Hence, a fused domain matrix with 2369 lines and 66 columns is obtained.e fused time-domain matrix is used to train the KELM model.

Experimental Results.
e CBA-KELM-based method is used to identify the fault classification of gearboxes with the multisensor fusion.
e CBA-KELM model is developed with the following parameters: 13 different kinds of timedomain indicators are calculated for fused datasets; the number of the neurons in the input layer is 65; the number of the neurons in the output layer is six; the kernel function is the radial basis function (RBF); the number of iterations is determined to be 100; the acceleration coefficient is identified as c1 � c2 � 1.49445; the population number is set as 10; and the parameter ranges of c and C are [0.1, 10] and [0, 1000], respectively.
To train the CBA-KELM model, trials are repeated 20 times, where 80% of the data samples are randomly exploited to be as the training dataset, and the remaining 20% data samples are used as the testing dataset.When applying the CBA method for KELM approach, the KELM parameters c and C should be determined after optimization.e fault Among all these trails, the average training accuracy is 100%, and the average testing accuracy is 96.25%, which means that the fault diagnosis and classification have higher accuracy based on the method shown in Table 2.   e proposed BA-KELM method utilizing the fused dataset is introduced in this paper.All the parameter settings are selected the same as the CBA algorithm except that the cloud theory is not used to optimize the weight of the BA algorithm.To train the BA-KELM model, the experiment is also repeated 20 times, with 80% of the data samples being used randomly as the training datasets and the remaining being used as the test datasets.Using the BA-KELM algorithm, the fault classification results of these 20 experiments are shown in Figure 11(b).Among these experiments, the average training accuracy rate is 100%.And the average testing accuracy rate is 95.66%, which is lower than the CBA-KELM approach.

Shock and Vibration
e FA technique is a new type of bionic intelligent algorithm, which originates from the group behavior of fireflies in the natural world through information exchange.
e FA-KELM model has the same parameter settings as the CBA-KELM method and the BA-KELM technique.From Figure 11(c), it can be known that the average training accuracy is 100% and the average testing accuracy is 95.52%.e KELM algorithm, which is implemented in a single learning step, and the parameters are usually assigned based on experience or trajectory [41].Its parameter settings are the same as those in the KELM model in BA-KELM and CBA-KELM algorithms.After 20 trials, the average training accuracy reaches 90.99% and the average testing accuracy reaches 87.90%, as shown in Figure 11(d).
e ELM model is developed with the following parameters: 13 different kinds of time-domain indicators are calculated for fused datasets; the number of the neurons in the input layer is 65; the number of the neurons in the output layer is six; and the number of the neurons in the hidden layer is 10.In order to train the ELM model, the experiment is repeated 20 times, 80% of the data samples are randomly used as training datasets and the rest are used as test datasets.After 20 trials, the fault classification results are shown in Figure 11(e).
e average training accuracy and average testing accuracy are 69.06% and 68.46%, respectively.e parameters of the BPNN-based technique are set as follows: the number of the input layer and output layer are set as one; the number of hidden layers is set as two; the number of input and output neurons is determined as 65 and 6, respectively; the number of the first and the second hidden layers is set as 14.Trails also are repeated 20 times, and the  is denotes that the fused signals of these multisensors are more meaningful, which can fully reflect the fault condition.
e classification results prove that information fusion technology can improve fault diagnosis accuracies.

Conclusions
is paper presents a novel approach for FDT of a wind turbine gearbox based on CBA-KELM with time-domain analysis and multisensor data fusion technology.In combination with feature-evel fusion and data-level fusion technology in parallel superposition, the original vibration signal is preprocessed in order to obtain the fused datasets.
e fused datasets contain a large number of samples involving six classifications by five accelerometers, including complex faults which are difficult to diagnose.And a fault classification method of KELM is optimized by using the CBA algorithm, which can obtain a high accuracy.e validity of the proposed FDT method is proved by using the fused data collected from the test platform system from the Qian Peng Company.Compared with a variety of other peer-to-peer methods, the proposed approach has the highest accuracy and stability.

Figure 2 :
Figure 2: Structure of a wind turbine gearbox.

Figure 7 : 2 nFigure 8 :
Figure 7: e diagnostic accuracy of the KELM technique with (a) the training sample and (b) the test sample.e diagnostic accuracy of the CBA-KELM method with (c) the training sample and (d) the test sample.
) Generate the drops of the cloud based on the information provided by the current optimal solution.pbc i (E x , E n , H e ) represents the population bats cloud, here H e � E n /c2, E x � x pbest , and E n �

Table 1 :
Description of the dataset.

Table 2 :
Comparisons with peer-to-peer classification.usingthismethodareshown in Figure11(f ).It is shown that the average training accuracies are 77.78% and the average testing accuracies are 68.97%.In Figures11(g)-11(l), the CBA-KELM, BA-KELM, FA-KELM, KELM, ELM, and BPNN algorithms based on unfused datasets are applied in this paper.e average training accuracies and training accuracies are very poor, no more than 55%.