This application claims the priority benefit of China application serial no. 202010134735.9, filed on Mar. 2, 2020. The entirety of the above-mentioned patent application is hereby incorporated by reference herein and made a part of this specification.
The invention relates to fault diagnosis, and more particularly relates to a method and a system for fault diagnosis of a gearbox of a wind turbine generator based on stacked denoising autoencoders.
In recent years, tremendous progress has been made in wind power generation. However, compared with conventional power systems such as those using coal, natural gas, etc., the operation and maintenance cost for a wind power system is higher. Therefore, it is necessary to adopt a proactive strategy for the maintenance of a wind power project, and, as a consequence, it is crucial to monitor the state, diagnose, predict, and implement health management to reduce the maintenance cost of a wind turbine. A majority of gearbox faults result from faulted gears. The maintenance required for a gearbox fault is complicated, and the cost for assembling/dissembling, transporting, and maintaining the gearbox is high. Therefore, to ensure normal operation of a wind turbine, further studies on the faults of gearboxes of wind power generators are required.
Chinese patent publication no. CN104792520A discloses a fault diagnosis method for a gearbox of a wind turbine generator system, and proposes a fault diagnosis method for a gearbox of a wind turbine generator system based on local mean decomposition and an optimized K mean value clustering algorithm, in which original vibration signals are decomposed and reconstructed, and an analysis is carried out on the reconstructed signals. Since the working conditions of a gearbox are complicated, decomposing the original signal may result in loss of high-dimensionality features, making it unable to obtain favorable fault diagnosis performance. Chinese patent publication no. CN108256556A discloses a fault diagnosis method for a wind generating set gearbox on the basis of a deep belief network, according to which a waveform database of working conditions of the gearbox is directly constructed, and original signals are input into a trained deep faith network to generate a waveform for comparison with the waveforms of different working conditions in the database. However, such an experience-based deep structure according to this method is unable to optimally extract features, and the efficacy of the fault diagnosis so rendered is significantly affected.
In view of the issues in the conventional art, the invention proposes a method and a system for fault diagnosis of a gearbox of a wind turbine generator based on stacked denoising autoencoders capable of facilitating the diagnosis efficacy based on the current diagnosis method.
For the above objectives, an aspect of the invention provides a method for fault diagnosis of a gearbox of a wind turbine generator based on stacked denoising autoencoders. The method includes:
(1) respectively obtaining a plurality of sets of original vibration signals under respective fault conditions, performing a Fourier transformation process and a normalization process on each of the original vibration signals to obtain a spectrum signal corresponding to each of the original vibration signals, and forming training data from all the spectrum signals;
(2) performing unsupervised training on a plurality of denoising autoencoders by using the training data;
(3) stacking together hidden layers of the respective denoising autoencoders that are trained, and adding the hidden layers to a logic regression layer to form the stacked denoising autoencoders; and
(4) optimizing the stacked denoising autoencoders by performing supervised training using a quantum-behaved particle swarm optimization method to obtain optimized stacked denoising autoencoders, so as to perform fault diagnosis by using the optimized stacked denoising autoencoders.
According to an embodiment of the invention, Step (2) includes:
(2.1) obtaining respective mapped signals by performing random mapping to the spectrum signals in the training data;
(2.2) adding non-masking noise to each of the mapped signals to obtain noise-contaminated signals, and mapping each of the noise-contaminated signals to the hidden layer; and
(2.3) obtaining respective reconstruction signals through reconstruction by a decoder in the hidden layer, and obtaining an optimal parameter of the denoising autoencoder through obtaining a minimum value of squared reconstruction errors according to the respective reconstruction signals and the respective spectrum signals.
According to an embodiment of the invention, an optimal parameter {θf, θg} of the denoising autoencoder is obtained by obtaining a minimum value of
wherein θf represents a parameter set {W1, bg}, θg represents a parameter set {W2, b2}, X2 represents the spectrum signal, X5 represents the reconstruction signal, and X5=σ(W2h+b2), h represents the hidden layer, h=σ(W1X4+b1), σ is a sigmoid function for realizing non-linear deterministic mapping, W1 represents a weight after mapping conversion of the hidden layer, b1 represents an offset after mapping of the hidden layer, X4 represents the noise-contaminated signal, W2 represents a weight after reconstruction, b2 represents an offset for reconstruction, X2
According to an embodiment of the invention, before Step (4), the method further includes:
initializing parameters of the stacked denoising autoencoders by using optimal parameters of the respective denoising encoders obtained in the unsupervised training, and updating weight values of the stacked denoising autoencoders by using a stochastic gradient descent method.
According to an embodiment of the invention, Step (4) includes:
(4.1) mapping a learning rate and a hidden layer number of the stacked denoising autoencoders as particle positions;
(4.2) obtaining an optimal individual position of each particle and a global optimal position of a swarm according to an adaption value of each particle in the swarm;
(4.3) obtaining a global optimal position of a corresponding particle according to the optimal individual position of each particle, and updating the particle positions according to the global optimal positions of the respective particles;
(4.4) repeating Steps (4.1) to (4.3) until an iteration stop condition is met, and using a swarm global optimal position that is obtained as the learning rate and the hidden layer number of the stacked denoising autoencoder.
According to an embodiment of the invention, an adaption value fitness (Nh,lr) of each particle in the swarm is obtained according to fitness
wherein lr represents the learning rate of the stacked denoising autoencoders, Nh represents the hidden layer number of the stacked denoising autoencoders, M represents a swarm size, xi represents actual values of the learning rate and the hidden layer number of the stacked denoising autoencoders, and yi represents predicted values of the learning rate and the hidden layer number of the stacked denoising autoencoders.
According to an embodiment of the invention, Step (4.3) includes: updating the particle positions according to
wherein mbest represents global optimal positions of all individuals, mbest j represents a center of optimal current positions in a jth dimension, Pi represents an optimal current position of an ith particle, Pij represents an optimal position of the ith particle in the jth dimension, Pgj represents an optimal position of a gth particle in the jth dimension, Pc
According to an embodiment of the invention, performing the fault diagnosis by using the optimized stacked denoising autoencoders includes:
obtaining a target vibration signal of a to-be-diagnosed gearbox of a wind turbine generator, and performing the Fourier transformation process and the normalization process on the target vibration signal to obtain a target spectrum signal;
extracting a fault feature signal by using the stacked denoising autoencoders, and identifying the fault feature signal by using a least squares support vector machine to obtain a fault type.
According to another aspect of the invention, a system for fault diagnosis of a gearbox of a wind turbine generator based on stacked denoising autoencoders is provided. The system includes:
a data processing module, configured to respectively obtain a plurality of sets of original vibration signals under respective fault conditions, perform a Fourier transformation process and a normalization process on each of the original vibration signals to obtain a spectrum signal corresponding to each of the original vibration signals, and form training data from all the spectrum signals;
a first training module, configured to perform unsupervised training on a plurality of denoising autoencoders by using the training data;
a stacked denoising autoencoder constructing module, configured to stack together hidden layers of the respective denoising autoencoders that are trained, and add the hidden layers to a logic regression layer to form stacked denoising autoencoders; and
a second training module, configured to optimize the stacked denoising autoencoders by performing supervised training using a quantum-behaved particle swarm optimization method to obtain optimized stacked denoising autoencoders, so as to perform fault diagnosis by using the optimized stacked denoising autoencoders.
According to yet another aspect of the invention, a computer-readable storage medium is provided. The computer-readable storage medium stores a program command, and when the program command is executed by a processor, the method for the fault diagnosis of the gearbox of the wind turbine generator based on the stacked denoising autoencoders according to the embodiments of the invention is realized.
In general, compared with the conventional art, the technical solution of the invention is capable of constructing the stacked denoising autoencoders by using the collected gearbox fault signals. Since the signals used in the construction process are obtained from the same gearbox, the optimized stacked denoising autoencoders are able to effectively extract the fault features in the signals of the same gearbox. The extracted fault features include high-dimensionality information of the original vibration signals. By inputting the feature signals into the least squares support vector machine, the fault types can be effectively identified. Accordingly, a strong basis is provided for pinning down the location of the fault of the gearbox and maintaining the gearbox. Accordingly, the stability and reliability of the operation of the equipment are ensured.
The accompanying drawings are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification. The drawings illustrate embodiments of the invention and, together with the description, serve to explain the principles of the invention.
Reference will now be made in detail to the present preferred embodiments of the invention, examples of which are illustrated in the accompanying drawings. Wherever possible, the same reference numbers are used in the drawings and the description to refer to the same or like parts.
In order to make the objectives, technical solutions, and advantages of the invention clearer, the following further describes the invention in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein only serve to explain the invention, but not to limit the invention. In addition, the technical features involved in the various embodiments of the invention described below may be combined with each other as long as such technical features do not conflict with each other.
In the examples of the invention, terms such as “first”, “second”, etc. are used to distinguish different objects, and do not necessarily imply a specific order or sequence.
S1: Collecting, by using an acceleration sensor, original vibration signals X under different fault conditions, performing a Fourier transformation process on the original vibration signals X to obtain spectrum signals X1, and normalizing the spectrum signals X1 to obtain normalized spectrum signals X2: A normalization formula is as follows:
wherein xnorm represents a normalized value, xmin and xmax respectively represent the minimum value and the maximum value in the spectrum signals X1.
S2: Constructing a stacked denoising autoencoder structure:
As shown in
X
3
=q
D(X3|X2) (2),
wherein D represents an original data set.
Non-masking noise is added to the mapped signals X3 to obtain noise-contaminated signals X4. Then, the obtain noise-contaminated signals X4 are mapped to a hidden layer h represented as follows:
h=f(X4,θf)=σ(W1X4+b1) (3),
wherein θf represents a parameter set {W1,b1}, W1 represents a weight matrix for mapping of the hidden layer, b1 is an offset vector for mapping of the hidden layer, and σ is a sigmoid function for realizing non-linear deterministic mapping, whose formula is represented as follows:
Then, reconstructed data X5 is obtained by the hidden layer h through decoder reconstruction. An expression thereof is as follows:
X
5
=g(h,θg)=σ(W2h+b2) (5),
wherein θg represents a parameter set {W2,b2}, W2 represents a weight matrix upon reconstruction, b2 represents an offset vector upon reconstruction, and an optimal parameter set {θf,θg} is obtained by obtaining the minimum value of squared reconstruction errors. A formula for obtaining the minimum value of the squared reconstruction errors is as follows:
wherein n is the number of the normalized spectrum signals X2.
As shown in
By using the corresponding parameters obtained during the process of unsupervised training, the parameters of the stacked denoising autoencoders are initialized. Then, a back propagation algorithm is adopted to perform supervised training on annotated information of the parameters, i.e., updating weight values by using a stochastic gradient descent method. The initial stacked denoising autoencoder structure is still unable to render the optimal efficacy in extracting features of gearbox faults. Here, a quantum-behaved particle swarm optimization algorithm is introduced to provide a stacked denoising autoencoder structure with a favorable feature extraction efficacy. Here, the annotated information of the parameters refers to a corresponding parameter set {θf,θg} corresponding to a corresponding individual denoising autoencoder obtained during the stepwise training process of the individual denoising autoencoder.
An iterative optimization formula for particles in the quantum-behaved particle swarm optimization algorithm is represented as follows:
wherein mbest and mbest j respectively represent centers of optimal current positions of all individuals and a jth dimension, Pi is an optimal current position of an ith particle, Pij and Pgj are respectively optimal positions of the ith and gth particles in the jth dimension, Pc
α=0.5+0.5×(tmax−t)/tmax (8),
wherein tmax represents a maximum number of iterations, and t represents the number of iterations.
The steps of iterative optimization are as follows:
(1) initializing the quantum-behaved particle swarm optimization algorithm including particle positions, an optimization range, compression/expansion factors, and the number of iterations, wherein a learning rate lr and a hidden layer number Nh of the stacked denoising autoencoders to be optimized are mapped as particle positions;
(2) calculating an adaption function of each particle in a swarm to obtain an optimal individual position of each particle and a global optimal position of the swarm, wherein the adaption function formula is represented as follows:
wherein M represents a swarm size, xi represents actual values of the learning rate and the hidden layer number of the stacked denoising autoencoder, yi represents predicted values of the learning rate and the hidden layer number of the stacked denoising autoencoder obtained from Formula (7), and an optimization objective is to obtain the minimum value of fitness(Nh, lr);
(3) calculating an optimal mean of the individual positions of all the particles in the swarm, i.e., a particle global optimal position, and updating the particle positions according to Formula (7);
(4) repeating (1) to (3) of the iterative process for the quantum-behaved particle swarm optimization algorithm until an iteration stop condition is met, wherein output optimization results are the learning rate lr and the hidden layer number Nh of the stacked denoising autoencoder.
Thus far, the stacked denoising autoencoder structure determined according to the original gearbox vibration signals under different fault states is obtained. Firstly, the pre-processed gearbox vibration signals are input into a single layer of the denoising autoencoder in the structure for unsupervised training, so as to obtain a denoising autoencoder capable of effectively extracting fault features from gearbox signals. Then, the individual denoising autoencoders are stacked and added to the logic regression layer to form a deep structure. By using the corresponding parameters obtained in the process of unsupervised training, the parameters of the deep structure are initialized, and a process of supervised reverse fine-tuning is performed. Then, the quantum-behaved particle swarm optimization algorithm is introduced to optimize the learning rate and the hidden layer number of the initial stacked denoising autoencoder structure, so as to obtain the stacked denoising autoencoders with a favorable efficacy in extracting the features of the gearbox vibration signals.
Step 3: Performing fault diagnosis on the currently input gearbox vibration signals:
After collecting a plurality of sets of original vibration signals, the sets of original vibration signals are respectively subjected to pre-processes such as Fourier transformation and normalization. The pre-processed signals are input into the optimized stacked denoising autoencoders, and a feature able to specify a fault type thereof is extracted from each signal. Then, the extracted fault features are input into a least squares support vector machine for fault classification. Accordingly, the state of the current signal is diagnosed.
In the embodiment of the invention, a least squares support vector machine using a Gaussian radial basis function as the core function may be adopted. The core function is represented as follows:
wherein σ represents a core parameter, and xi and xj respectively represent an ith sampling value and a jth sampling value. The least squares support vector machine is represented as follows:
wherein αi represents a Lagrange multiplier, yi represents −1 or 1 of a class, β represents a compensating parameter, and l represents a sampling number.
Compared with the conventional art, the embodiment of the invention constructs the stacked denoising autoencoders by using the collected gearbox fault signals. Since the signals used in the construction process are from the same gearbox, the optimized stacked denoising autoencoders are able to effectively extract the fault features in the signals of the same gearbox. The extracted fault features include high-dimensionality information of the original vibration signals. By inputting the feature signals into the least squares support vector machine, the fault types can be effectively identified. Accordingly, a strong basis is provided for pinning down the location of the fault of the gearbox and maintaining the gearbox. Thus, the stability and reliability of the operation of the equipment are ensured.
a data processing module 201, configured to respectively obtain a plurality of sets of original vibration signals under respective fault conditions, perform a Fourier transformation process and a normalization process on each of the original vibration signals to obtain a spectrum signal corresponding to each of the original vibration signals, and form training data from the spectrum signals;
a first training module 202, configured to perform unsupervised training on a plurality of denoising autoencoders by using the training data;
a stacked denoising autoencoder constructing module 203, configured to stack together hidden layers of the respective denoising autoencoders that are trained, and add the hidden layers to a logic regression layer to form stacked denoising autoencoders; and
a second training module 204, configured to optimize the stacked denoising autoencoders by performing supervised training using a quantum-behaved particle swarm optimization method to obtain optimized stacked denoising autoencoders, so as to perform fault diagnosis by using the optimized stacked denoising autoencoders.
Regarding the details of the respective modules, reference is made to the above descriptions made for the embodiment of the method, and the same details will not be repeated in the following for the embodiment of the invention.
In another embodiment of the invention, a computer-readable storage medium is provided. The computer-readable storage medium stores a program command. By executing the program command by a processor, the method for fault diagnosis of the gearbox of the wind turbine generator based on the stacked denoising autoencoders according to the embodiment is realized.
Analysis on Experimental Results
Owing to the limitations of objective conditions, it is difficult to collect a large amount of fault data for research within a short period of time. Therefore, the wind turbine generator gearbox fault data adopted in the embodiment is provided by a motor-driven planetary gearbox fault simulation system. The system includes a motor, a parallel shaft gearbox, a planetary gear, a low-speed bearing, a high-speed bearing, and a magnetic brake, and is capable of simulating various different gearbox fault conditions. As shown in
After performing the Fourier transformation process on 80 sets of training data, the normalization process was performed on the data according to Formula (1). Then, the processed training data were input into the individual denoising autoencoders to perform the unsupervised training according to Formulae (2) to (6). The initial stacked denoising autoencoders were formed through stacking, and the quantum-behaved particle swarm optimization algorithm was initialized. The stacked denoising autoencoders were optimized through supervised training according to Formulae (7) to (9), and a t-SNE non-linear dimensionality reduction algorithm was adopted to reduce the dimensionality of the extracted high-dimensionality features to two, so as to validate the feature extraction efficacy of the stacked denoising autoencoders. The result of a single trial of feature extraction is as shown in
After the stacked denoising autoencoders were optimized, 120 sets of test data were subjected to the Fourier transformation process and the normalization process. Then, the processed data were input into the stacked denoising autoencoders to extract fault features. Then, the fault feature signals were input into the least squares support vector machine for fault type identification, so as to obtain the percentage of the correctly diagnosed samples with respect to the total samples. Accordingly, the diagnosis accuracy according to the method of the invention was obtained. The diagnosis result is as shown in Table 1.
100%
According to Table 1, among the four different working conditions of the gear box, the lowest diagnosis accuracy was as high as 93.33%, and the mean diagnosis accuracy was as high as 96.67%. The results suggest the method for fault diagnosis of the gearbox of the wind turbine generator based on the stacked denoising autoencoders proposed in the invention renders a favorable diagnosis efficacy, and provides a different line of thinking as well as a novel approach for the fault diagnosis of gearboxes of wind power generators.
It is noted that, based on practical needs, the respective steps/parts in the invention may be further divided into more steps/parts, or two more steps/parts or portions of the steps/parts may be combined to form new steps/parts to realize the objectives of the invention.
The method according to the invention may be implemented in hardware or firmware, or realized as software or computer codes in a recording medium (e.g., CD ROM, RAM, soft disk, hard disk, or magneto-optical disc), or computer codes originally stored in a remote recording medium or non-transitory machine-readable medium and downloadable via a network to be stored in a local recording medium and processed by software stored in a recording medium of a general-purpose computer, a dedicated processor, or a programmable or dedicated hardware component (e.g., ASIC or FPGA). It is understood that a computer, a processor, a microprocessor controller, or a programmable hardware component include a storage component (e.g., RAM, ROM, flash memory, etc.) capable of storing or receiving software or computer codes. When the software or computer codes are accessed and executed by a computer, a processor, or a hardware component, the processing method described herein is realized. In addition, when a general-purpose computer accesses the codes for realizing the processes described herein, by executing the codes, the general-purpose computer is turned into a dedicated computer for executing the processes described herein.
It will be apparent to those skilled in the art that various modifications and variations can be made to the structure of the present invention without departing from the scope or spirit of the invention. In view of the foregoing, it is intended that the present invention cover modifications and variations of this invention provided they fall within the scope of the following claims and their equivalents.
Number | Date | Country | Kind |
---|---|---|---|
202010134735.9 | Mar 2020 | CN | national |