Repetitive Transient Extraction Using the Optimized SES Entropy Wavelet for Fault Diagnosis of Rotating Machinery

Repetitive transients are usually generated in the monitoring data when a fault occurs on the machinery. As a result, many methods such as kurtogram and optimized Morlet wavelet and kurtosis method are proposed to extract the repetitive transients for fault diagnosis. However, one shortcoming of these methods is that they are constructed based on the index of kurtosis and are sensitive to the impulsive noise, leading to failure in accurately diagnosing the fault of the machinery operating under harsh environment. To address this issue, an optimized SES entropy wavelet method is proposed. In the proposedmethod, the optimized parameters including bandwidth and central frequency of Morlet wavelets are selected. ,en, based on the wavelet coefficients decomposed using the optimized Morlet wavelet, the SES entropy is calculated to select the scales of wavelet coefficients. Finally, the repetitive transients are reconstructed based on the denoising wavelet coefficients of the selected scales. One simulation case and vibration data collected from the experimental setup are used to verify the effectiveness of the proposed method. ,e simulated and experimental analyses showed that the signal-to-noise ratio (SNR) of the proposed method has the largest value. Specifically, the SNR in the experimental analysis of the proposed method is 0.6, while that of the other three methods is 0.043, 0.0065, and 0.0045, respectively. ,erefore, the result shows that the proposed method is superior to the traditional methods for repetitive transient extraction from the vibration data suffered from impulsive noise.


Introduction
Fault diagnosis plays a vital role in ensuring long-term safe running of rotating machinery for avoiding huge economic loss and casualties.
us, many fault diagnosis methods based on the collected monitoring signals such as sound [1], infrared images [2], and current [3] are carried out. Of all the different types of monitoring signal, vibration signal contains abundant information of machinery health conditions and thus is the most analyzed signal, and many signal processing based methods have been proposed to process vibration signals for fault diagnosis [4,5]. e abundant information refers to the repetitive transients generated periodically when faults occur on the rotating machinery [6]. By analyzing the frequency of these transients, it is easy to know which part has faults or dynamic changes [7].
Although vibration signals containing repetitive transients can be used for diagnosing the fault, machinery usually operates under harsh environment and the transients are submerged in noise. As a result, it is difficult to extract useful information from vibration signals and inaccurate results probably obtained with ineffective methods. To solve this problem, Qiao et al. [8][9][10][11] applied stochastic resonance to fault diagnosis, making use of the enhancement of noise to periodicity signals generated due to the faults. Multistable stochastic resonance [12] can also be found for fault detection. However, selecting the suitable optimal parameters of stochastic resonance is not an easy thing; although the stochastic resonance can exhibit the energy of noise for transient extraction, these signal processing-based methods cannot be widely used in real applications. Moreover, blind deconvolution technique [13] and minimax concave regularized sparse representation [14] are also used for extracting repetitive transients for fault diagnosis, but they are not robust to heavy noise.
Time-frequency signal processing techniques including wavelet transform and empirical mode decomposition are also commonly used to process vibration signals for extracting useful fault information. e empirical mode decomposition and its variants [15] can be found to extract repetitive transients in fault diagnosis of rotating machines. ese methods decompose vibration signals into different models of intrinsic model functions and then the signals can be classified into nonstationary and nonlinear components for analysis. However, it should be noted that these methods have many problems such as endpoint effects, modal aliasing, and lacking mathematical theory, which limit the use of these methods. Wavelet transform is another commonly seen time-frequency analysis method and is performed by decomposing signals into different scales of frequency bands. Continuously, the frequency bands of repetitive transients can be selected by choosing scales for the reconstruction of repetitive transients. Zhang et al. [16] achieved the bearing fault diagnosis based on multiple Q-factor Gabor wavelets. Mishra et al. [17] used envelope analysis and wavelet denoising to diagnose rolling element bearing fault under slow-speed operation. In recent years, more improved wavelet-based diagnosis methods including M-band flexible wavelet transform [18], empirical wavelet transform [19,20], tunable Q-factor wavelet transform [21], and dynamic Bayesian wavelet transform [22] are proposed. But there exist some common issues for these wavelet-transform-based methods. First, it is difficult to choose a suitable wavelet function that should be similar to the transients. A random selection of wavelet function may lead to the wrong representation of transients and then no useful information is probably obtained. Second, there are no effective theories to guide the selection of wavelet scales. e wavelet scales representing the frequency of repetitive transients should be selected, but which scales should be selected should be guided with an effective index.
Kurtosis is sensitive to impulsive signals and exhibits large values when there are impulsive signals [23]. Based on these concepts, many useful kurtosis-based diagnosis methods have been proposed. Antoni [24] constructed a kurtogram method, which uses short-time Fourier transform to decompose signals into different frequency bands and select the optimal frequency band based on the index of kurtosis. Wang [25] proposed an enhanced kurtogram method for fault diagnosis of rolling element bearings based on the binary wavelet packet transform. Chen [26] constructed an improved fast spatialspectral ensemble kurtosis kurtogram and applied it to mechanical signature analysis of short duration data from rotating machinery. By combining minimum entropy deconvolution and spectral kurtosis, a new l0-norm embedded MED method was proposed for roller element bearing fault diagnosis at the early stage of damage [27]. A sub-band averaging kurtogram was proposed by incorporating with dual-tree complex wavelet packet transform [28]. Adaptive correlated kurtogram was constructed based on scale-space representation and empirical wavelet transform [29]. Qin [30] proposed an optimized Morlet wavelet and kurtosis method, in which the kurtosis is used to select the optimal scales. Although kurtosis is sensitive to repetitive transients, it is also more sensitive to the single impulsive noise. e single impulsive noise is widely seen because the machinery usually operates under harsh environment, and disturbances are inevitable. e disturbance for short time can lead to impulsive noise. As a result, these kurtosis methods probably fail to detect fault from the data containing impulsive noise. Infogram was proposed to overcome this issue by calculating square envelope spectrum entropy, but this method is not robust to heavy noise and there are usually several frequency bands appearing with large values [30]. Afterward, the square envelope spectrum entropy has attracted much attention and many variants of infogram have been constructed [31][32][33][34].
To effectively extract repetitive transients from signals suffered from heavy noise including the Gaussian white noise and the impulsive noise, a new fault diagnosis method is proposed based on the optimized Morlet wavelet and the square envelope spectrum entropy. In the proposed method, the Morlet wavelet function which is similar to the impulsive signals is used. e parameters of Morlet including bandwidth and the central wavelet frequency are optimized according to the Shannon entropy of wavelet coefficients. Based on the wavelet transform using the optimal Morlet wavelet function, square envelope spectrum entropy is calculated for guiding the selection of wavelet scales. Finally, the repetitive transients can be reconstructed based on the selected wavelet coefficients after denoising. e effectiveness of the proposed method is demonstrated using simulation signal and experimental signal, respectively. e proposed method can be used for the fault diagnosis of machinery suffered from both heavy Gaussian white noise and impulsive noise. Impulsive noise is widely seen in the collected data due to random disturbance, and traditional methods such as kurtogram cannot be used to process these data effectively. us, the proposed method is suitable for a wide range of applications in real engineering.
Based on the previously mentioned statements, the contributions of this paper are summarized as follows: (1) An optimized SES entropy wavelet is proposed to extract the repetitive transients for fault diagnosis, and the proposed method presents good performance even when impulsive noise is mixed in the diagnosed signal (2) e Morlet wavelet function which is similar to impulsive signals is used and its parameters are optimized and selected by calculating the Shannon entropy of wavelet coefficients (3) Square envelope spectrum entropy is introduced for guiding the selection of wavelet scales e rest of this paper is organized as follows. Section 2 is devoted to introduce the basic theory of wavelet transforms. In Section 3, the detailed procedure of the proposed method is present. In Section 4, the proposed method is demonstrated based on simulated and experimental signals, respectively. Finally, the conclusions are drawn in Section 5.

2
Shock and Vibration

The Theory of Wavelet Transform
Wavelet transform is widely used to analyze the signals in the time-frequency domain and is efficient for capturing the time of occurrence of a particular frequency. Accordingly, this method is robust for analyzing nonstationary signals and thus used in the proposed method. To help a better understanding of the proposed method, the theory of wavelet transform is introduced in detail. Wavelet can be considered as a short-duration wave growing and decaying over a limited period. e wavelet denoted using ψ(t) satisfies the following two conditions: (1) By scaling and shifting ψ(t) through a and b, respectively, a set of basis functions ψ a,b (t) can be obtained according to Wavelet transform is a powerful mathematic tool for obtaining a representation, decomposition, and reconstruction of signals. Specifically, wavelet transform is performed by decomposing the studied signals into a combination of time-shifted and dilated or compressed local basis functions according to where x(t) denotes the studied signal b. e asterisk represents the computation of complex conjugate, and thus, wavelet transform can be considered as a convolution of the signal x(t) with a window ψ a,b (t) that is dilated by a and shifted in time by b. As a result, each wavelet can be located at a different position around b along the time axis. Furthermore, with the help of a, wavelets of different frequencies can be obtained. a with large values indicating that the frequency of wavelet of frequency is low, while a with small values correspond to high frequency.

The Proposed Method
In this section, a repetitive transient extraction method based on an optimized SES entropy wavelet is constructed and introduced for machinery fault diagnosis. e flowchart of the proposed method is graphically shown in Figure 1, and its detailed procedures are described as follows: (1) First, the monitoring data are collected from machinery and used for fault diagnosis, and some parameters of the Morlet wavelet including f b and f c are initialized. f b and f c represent the bandwidth and the central wavelet frequency, respectively.
e Morlet wavelet which is similar to the impulsive signal generated due to faults is used in wavelet transform. e Morlet transform is shown as It can be found that the oscillation attenuation depends greatly on f b , specifically, the increase of f b enhances the oscillating of the Morlet wavelet, f c

Monitoring data
Initialize f b =x m and f c =R 1 for the Morlet wavelet Calculate wavelet transform using the above Morlet wavelet Select the optimal f c marked as f co according to the minimum E Shock and Vibration controls the oscillatory frequency, and the increase of f c leads to the decrease of the frequency resolution of the Morlet wavelet. To find the optimal f c , f b should be set as a large value. erefore, f b is set to be 30, while f c is set as 0.2. f c will be iterated from 0.2 to 1 with a step of 0.1, while f b remain unchanged. (2) Wavelet transform is calculated based on the Morlet wavelet according to equation (4), so the wavelet coefficients of different scales can be obtained.
where p i is the distribution sequence calculated based on wavelet coefficients according to (4) If f c does not exceed R2, iterate f c with an increase in T c and then compute the new wavelet Shannon entropy according to these above steps. f co and f bo , respectively. As a result, wavelet coefficients of different scales can be obtained and denoted using |W(a i , t)|. (8) e square envelope of these coefficients can be calculated using (9) Based on the discrete Fourier transform, the square envelope spectrum (SES) can be calculated using where α denotes the cyclic frequency and F s is the sampling frequency. (10) e SES negentropy of these Fourier coefficients is calculated according to the following equation: (11) e characteristic scales whose SES negentropy satisfies the following conditions are selected as the useful scales for reconstructing repetitive transients: where T is a threshold value and set to be 0.75. In this way, the wavelet coefficients with larger SES negentropy can be selected and they mostly contain repetitive transients. (12) e wavelet coefficients of the selected characteristics are denoised using soft thresholding to enhance the ability of transient extraction according to where θ is a threshold value.
In a continued manner, the reconstructed repetitive transients can be used for fault diagnosis by calculating the Hilbert envelop spectrum of these transients.

Transient Extraction for Simulated Signal of Bearing
Faults. In order to demonstrate the effectiveness of the proposed method for repetitive transient extraction, the simulated signals of the fault bearing are considered in this case study. e signal is simulated using the following model: where A m represents the amplitude value of the simulated signal. A m is randomly selected and varies from 0.6 to 1.6. u(t) denotes a unit step function, β w denotes structure damping ratio, and β w � 900. e resonance frequency f re and the sampling frequency f s are 3000 Hz and 12 kHz, respectively. e total sampling time is 0.5 s. Moreover, m denotes the number of impulses that periodically appear in the simulated signal. M is the total number of the simulated impulses and is set to 200. f o is the fault characteristic frequency and set to 50 Hz in this simulation.
In addition, to simulate the impulsive noise due to the harsh environment, I(t) is also added to the simulated signal, and its amplitude, resonance frequency, and bandwidth is 8 g, 5 kHz, and 800 Hz, respectively. N(t) denotes the additive Gaussian white noise with zero mean and satisfies the following equation: where D is the intensity of the Gaussian white noise and is selected as 0.125. e time domain of waveforms of the simulated signal is shown in Figure 2. From Figure 2(a), it can be found that transients are generated repetitively with a periodicity of 1/f o s. Unlike the repetitive transients, the impulsive noise appears with a large amplitude. is noise is commonly seen in the real monitoring data due to disturbance from the operational environment. In addition to the impulsive noise, Gaussian white noise is added to the signal shown in Figure 2(a) and the repetitive transients are seriously deteriorated in the signal shown in Figure 2(b). Obviously, it is difficult to extract the repetitive transients for fault diagnosis from the signal using traditional signal processing based fault diagnosis methods such as kurtogram. e proposed method is applied to processing this signal, and the time-domain waveform of the extracted repetitive transients is shown in Figure 3(a). It can be seen that the repetitive transients can be successfully extracted using the proposed method by removing both the Gaussian white noise and the impulsive noise. en, the Hilbert envelope spectrum of the signal shown in Figure 3(a) is calculated and presented in Figure 3(b). e amplitude of fault frequency and its harmonic frequency such as the second, third, and fourth are high and clearly seen. It can be inferred that a fault occurs on the monitoring bearing, which agrees with the simulated data. erefore, the proposed method is effective for repetitive transient extraction for fault diagnosis.
For comparison, the traditional methods including kurtogram, infogram, and the optimized Morlet wavelet and kurtosis method are used to process this simulated signal. Kurtogram is a widely used method to extract repetitive transients for fault diagnosis, so this method is considered for comparison. e result of the kurtogram is shown in Figure 4, and it can be found that the large value appears at the center frequency (f c ) of 4968.75 Hz and at the level of 7, which can be detected as the optimal frequency band and marked using a red color line. Consequently, the signal at this frequency band is filtered for extracting the repetitive transients. e extracted repetitive transients using kurtogram are shown in Figure 5(b). Obviously, the impulsive noise is contained in the repetitive transients, indicating the kurtogram cannot be used for repetitive transient extraction when there exists heavy noise in the collected signal. e envelop spectrum of the signal shown in Figure 5(b) is presented in Figure 5(c). ere is no useful information presented in the spectrum graph and the bearing fault cannot be successfully detected. Kurtogram is constructed by calculating the kurtosis index in different frequency bands decomposed using a short Fourier transform. e kurtosis index is very sensitive to impulsive noise and has a very large value under a sole impulsive noise. As a result, the frequency band of the impulsive noise will be wrongly selected for fault diagnosis.
Infogram is also proposed by Antoni for overcoming the shortcomings of the kurtogram, which is constructed based on Shannon entropy. e results of the infogram are shown in Figure 6. ere are several frequency bands with large values because the infogram is not robust to noise. e frequency band with the largest value is at the center frequency of 4500 Hz and at the second level, which is also marked using a red color line. e extracted repetitive transients are shown in Figure 7(b), and we can find that both the repetitive transients and the impulsive noise are contained in the repetitive transients. As shown in Figure 7(c), although the amplitude of the fault frequency can also be seen, there are also some disturbances in the envelope spectrum due to the impulsive noise. In sum, the infogram cannot be used to remove the impulsive noise because of its coarsely calculated procedure.
Like the kurtogram, the Morlet wavelet and kurtosis method is constructed based on the index of kurtosis and is sensitive to impulsive noise. e time-domain waveform of the extracted repetitive transients using the Morlet wavelet and kurtosis method is shown in Figure 8(a), in which the impulsive noise rather than the periodically appeared impulsive is clearly seen. As a result, the amplitude at the fault frequency has a small value and there are just low-frequency components in the envelop spectrum, as shown in Figure 8(b). It is because the scale of the signal containing impulsive noise has the largest kurtosis and thus is wrongly selected as the repetitive transients in the comparison method.
For quantitative advantages of the proposed method, the SNR is used for the envelope spectrum of the repetitive transients. e expression of SNR is generally written as Shock and Vibration    Shock and Vibration where i (0, 1, . . . N) denotes the order of the harmonics and N is specified as 4. A i represents the amplitude of the extracted ith harmonic and E is the sum of the amplitude of envelope spectrums. e SNR of the proposed method, infogram, kurtogram, and the Morlet wavelet and kurtosis method can be calculated as 0.135, 0.092, 0.011, and 0.0124, respectively, which verify the effectiveness of the proposed method.

Transient Extraction for Real Signal Collected from a Fault
Bearing. In this section, real signals collected from a fault bearing are used to verify the effectiveness of the proposed method. e data are provided by the bearing data center of Case Western Reserve University [35] and is commonly used for the demonstration of fault diagnosis methods. e  Shock and Vibration experimental setup is shown in Figure 9, which consists of a 2 HP motor (left), a torque transducer (center), and a dynamometer (right). e bearing of 6205-2Rs JEM SKF is used for testing, which is installed under the motor shaft. Some structural parameters of the bearing are presented and listed in Table 1. An accelerometer sensor, which is attached to the housing with magnetic bases, is used to collect vibration signals with a sampling of 12 kHz. e vibration signal of fault outer race was collected and processed using the proposed method. e fault is a single-point fault that is produced using electro-discharge machining and the fault diameter is 7 mils (1 mil � 0.001 inches). e rotation speed of the motor is 1797 r/min, and the corresponding outer fault frequency is 91.4 Hz. e time-domain waveform of the vibration signal is shown in Figure 10(a), and it can be found that the repetitive transients are submerged in the noise. Moreover, an impulsive noise whose amplitude value is 5 is introduced into the signal and the corresponding waveform is shown in Figure 10(b). e proposed method is applied to extract repetitive transients from this signal for fault diagnosis. e timedomain waveform of the repetitive transients is shown in Figure 11(a). It can be found that the simulated impulsive noise is not contained in the signal, indicating that the proposed method is not sensitive to the noise. To further verify the effectiveness of repetitive transients extracting and diagnosis the fault of bearing, the envelope spectrum of the signal is given in Figure 11(b). e amplitude of the outer    Shock and Vibration race fault frequency and its harmonic frequency can be clearly seen. It can be inferred that the bearing suffers from the outer fault, which agrees well with the fact. e traditional methods are also used for comparison. e results of the kurtogram show that the center frequency of the repetitive transients is 0 Hz at the 7.5 th level (seen in Figure 12), but the extracted signal of the corresponding frequency marked using the red line contains the impulsive noise (seen in Figure 13(b)), and as shown in Figure 13(c), there is no useful information in its Hilbert envelop spectrum. e center frequency at 0 Hz and 2 nd level is selected for repetitive transient extraction using infogram, as shown in Figure 14. e extracted repetitive transients are shown in Figure 15(b), and no useful information can also be found in its Hilbert envelop spectrum,    Figure 13: (a) e data after prewhitening, (b) the extracted repetitive transients using kurtogram, and (c) the Hilbert envelope spectrum of (a).
as shown in Figure 15(c). As shown in Figure 16, a similar result is found, when the optimized Morlet wavelet and kurtosis method is applied to the signal. e SNR of infogram, kurtogram, and the Morlet wavelet and kurtosis method can be calculated as 0.043, 0.0065, and 0.0045, respectively, while the SNR of the proposed method is 0.6. In sum, the proposed method is superior to the traditional methods and is robust to repetitive transient extraction from the signal suffering impulsive noise.

Conclusions
is paper presents a repetitive transient extraction method based on the optimized SES entropy wavelet for machinery fault diagnosis. In the proposed method, the parameters including bandwidth and central wavelet frequency are optimized by computing wavelet Shannon entropy. e wavelet transform of monitoring signal using the Morlet wavelet under these optimized parameters is performed.   Next, the square envelope of wavelet coefficients is calculated and then the SES entropy can be obtained. e repetitive transients can be reconstructed using the denoised wavelet coefficients for fault diagnosis. e effectiveness of the proposed method is verified by the simulated and experimental signals. Traditional methods including kurtogram, infogram, and the optimized Morlet wavelet and kurtosis method are used for comparison. e results show that the proposed method is robust to extract repetitive transients submerged in impulsive noise, while traditional methods usually fail to do that because they are sensitive to impulsive noise. erefore, the proposed method has a significant advantage over traditional methods and is suitable for fault diagnosis of machinery operating under harsh environment. In future, some useful optimal methods [36][37][38] will be introduced to make the proposed method more robust.
Data Availability e data that support the findings of this study are available on request from the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest.