A 1.5D Spectral Kurtosis-Guided TQWT Method and Its Application in Bearing Fault Detection

Bearings are the key parts of rotating machinery, and their operation status is related to the operation safety of the whole equipment. Vibration signals often contain periodic impulse components which can reflect the fault state of bearings. However, due to the interference of signal transmission path and the influence of operating environment noise, the periodic impulse components in the signal are often submerged by the nonperiodic transient impulse components, modulation harmonic components, and noise components.,erefore, the core problem of bearing fault diagnosis theory is used to accurately extract the frequency band of bearing fault state information and suppress the frequency band of interference information. In this paper, the signal is processed by the tunable Q-factor wavelet transform (TQWT), the midfrequency band of the signal is tightly divided by selecting different Q-values, and the 1.5D spectral kurtosis defined in frequency domain is used to select the optimal subband. Simulated analysis shows that this method can avoid low-frequency harmonic interference, nonperiodic transient impulse components, and strong noise components in the time domain.,erefore, it can effectively realize the selection of the subbands of periodic impulse components and effectively extract fault feature information. ,rough experimental signal analysis, TQWT has good sparsity decomposition characteristics and can reasonably divide the signal frequency band, so as to separate the useful fault characteristic frequency band and interference frequency band. At the same time, compared with the kurtosis index, 1.5D spectral kurtosis has better robustness and resolution for low signal-to-noise ratio signals, which can achieve the purpose of fault characteristic band extraction.


Introduction
Rolling bearing is an important component of rotating machinery, and its operation status is directly related to whether the whole equipment can run safely and smoothly. erefore, the fault analysis of rolling bearing has become a hot issue in the field of signal analysis [1][2][3]. Considering the problems of signal transmission path, signal acquisition convenience, and cost, vibration signal is a widely used analysis medium. So many analysis and processing methods based on vibration signal are derived [4][5][6]. Especially, in recent years, with the development of mathematical model construction and parameter optimization method of processing algorithm [7][8][9], signal processing algorithm has been greatly improved in terms of accuracy, efficiency, and adaptability.
e key link of bearing vibration signal processing is filtering process. rough certain analysis methods, the fault information within the specific frequency band can be retained to the maximum extent, and the interference components outside the frequency band can be suppressed to the maximum extent, so as to realize the signal filtering and fault information extraction.
ere are two main problems: one is to use appropriate methods to segment the frequency band, which makes the subband after segmentation have certain sparsity, and the other is to use appropriate index to evaluate the subband after segmentation, so that the subband with fault impulse component has certain specificity. e kurtogram algorithm [10] proposed by Antoni is considered to be a typical signal processing method, which can divide the signal plane into subbands with different bandwidth and center frequency through filter banks and select the subbands by the kurtosis index. Subsequently, many scholars began to improve the frequency band segmentation process and index selection of this kind of algorithm. Lei et al. [11] used wavelet packet decomposition to replace the equivalent filter banks in the kurtogram algorithm. Wang and Tsui [12] used dynamic Bayes to optimize the wavelet parameters and adjust the center frequency and bandwidth of the equivalent filtering characteristics, so that it can more accurately locate the fault impact subband. Qin et al. [13] constructed a new M-band flexible wavelet transform to identify the potential fault features in the measured signal. is transform provides more intensive time-frequency plane sampling and retains the tunable filter parameters and dilation factors. e simulated and experimental results show that this transform has better denoising performance than analytic flexible wavelet transform, orthogonal wavelet transform, and biorthogonal wavelet transform. Meanwhile, Qin et al. [14] proposed a novel K-SVD algorithm to get the optimized transient dictionary for the one-dimensional signal. Simulated and experimental results show that this method has advantages over the traditional K-SVD method and infogram method in periodic transient impulse components detection. e kurtosis index is too sensitive to impulse components and is easily disturbed by nonperiodic transient impulses and noise. Based on these considerations, Antoni further defines the square envelope entropy index and the square envelope spectral entropy index in the time domain and frequency domain, respectively, and constructs the infogram algorithm [15]. Li et al. [16] proposed a multiscale clustering gray infogram, which combined two negative entropies with the gray method. Moshrefzadeh and Fasana [17] proposed an autogram method using square envelope unbiased autocorrelation instead of time-domain signal. rough this operation, uncorrelated interference components can be discarded.
Selesnick designed the tunable Q-factor wavelet transform (TQWT) from the perspective of frequency domain filtering [18]. TQWT can make the sparsity of frequency band division more diversified by adjusting its Q-factor. e key problem of TQWT is how to determine its parameter value and make the selection of frequency band more accurate by adjusting the oscillation characteristics of wavelet. Zhao et al. [19] used the particle swarm optimization algorithm to optimize TQWT parameters, combined with the spectral kurtosis index to select the optimal subband. Ma et al. [20] also used particle swarm optimization to optimize the parameters of TQWT and further combined with frequency slice wavelet transform (FSWT) to select the frequency band. Ding et al. [21] constructed a filter bank with almost constant bandwidth by adjusting the Q-factor and decomposition level of TQWT.
In this paper, a 1.5D spectral kurtosis-guided TQWT method is proposed. In this method, the 1.5D spectral kurtosis of the square envelope sequence of wavelet subband signals with different Q-values and scales is calculated to replace the traditional kurtosis index to select the optimal subband.
is method can suppress the noise and the nonperiodic transient impulse interference without coupling relationship. e outline of this paper is as follows. Section 2 introduces the 1.5D spectral kurtosis-guided TQWT method and its superiority in signal periodic feature extraction. Sections 3 and 4 verify the effectiveness of the proposed method with the aid of simulated and experimental data.

1.5D Spectral Kurtosis-Guided TQWT
2.1. TQWT. Compared with the traditional discrete wavelet transform, TQWT can determine the wavelet basis function by designing the Q-factor (the Q-factor is the ratio of signal center frequency to bandwidth, Q � f c /B ω ) and the redundancy factor according to the signal oscillation characteristics. is design makes the choice of wavelet function more flexible, makes the signal adaptability stronger, and makes it easier to control the sparsity of frequency band division. TQWT uses band-pass filter banks to realize signal decomposition and reconstruction in an iterative manner. e three-level decomposition and reconstruction is taken as an example. e corresponding dual channel filter banks are shown in Figure 1.
In Figure 1, H 1 (ω) and H 0 (ω) are the frequency response functions of the high-pass filter and the low-pass filter (H 1 (ω) and H 0 (ω) can be expressed as equations (1) and (2)), H * 1 (ω) and H * 0 (ω) are the conjugate complex numbers of H 1 (ω) and H 0 (ω), respectively, β � 2/(Q + 1) is the high-pass scaling factor, and α � 1 − β/r is the low-pass scaling factor, where r is the redundancy factor: e maximum number of decomposition layers of TQWT can be obtained by e scale factors α and β can be substituted into equation In summary, the filtering characteristics of TQWT change with the Q-factor and the redundancy factor r. As the number of decomposition layers increases, the center frequency of the filter gradually decreases and the bandwidth also decreases. When the sampling frequency of the signal is f s , the center frequency f c and bandwidth B ω can be expressed as equations (5) and (6), respectively:

1.5D Spectral Kurtosis.
For a discrete sequence of signal x(t), the system response is represented by X(t), which can be expressed as where H(t, f) is the transfer function and dx(f) is the spectral increment. e second-order spectral distance of X(t) can be expressed as When 2n � 2, equation (8) can be abbreviated as Furthermore, the fourth-order spectral cumulant can be expressed as e spectral kurtosis can be defined by fourth-order normalized cumulants, as shown in For a signal x(t), the diagonal slice of the third-order where E · { } represents mathematical expectation. e 1.5D spectrum of the original signal x(t) can be obtained by performing Fourier transform on the three diagonal slices as e spectral distance of the spectrum discrete sequence is calculated and brought into equation (11), and the 1.5D spectral kurtosis value 1.5 D K is obtained as In order to compare the robustness of the kurtosis index defined in the time domain with the 1.5D spectral kurtosis index defined in the frequency domain against noise, different degrees of noise are added to the periodic impulse simulated signal. e mathematical model of the simulated signal can be expressed by Figure 1: e corresponding dual channel filter banks of (a) decomposition process and (b) reconstruction process.
Shock and Vibration 3 (15) where A is the amplitude, ξ is the damping ratio, f o is the natural frequency of the system, f outer is the characteristic frequency, and n(t) is the noise. e parameters of the simulated signal are listed in Table 1.
Figures 2(a)-2(p) are the time-domain waveforms and 1.5D spectrum of simulated signals with different SNR. Figure 2(q) shows the variation trend of the kurtosis index in the time domain and the 1.5D spectral kurtosis index in the frequency domain with the SNR. It can be seen from the figure that when the periodic impulse in the time-domain signal is gradually overwhelmed by noise, the spectral impact component composed of the fault characteristic frequency and its harmonics in the 1.5D spectrum is still obvious. is shows that the 1.5D spectrum has a relatively strong ability to suppress noise. By comparing the variation trends of the time-domain kurtosis index and the 1.5D spectral kurtosis index, it can also be found that when the SNR ratio is low, the 1.5D spectral kurtosis has a higher slope, indicating that it has a higher resolution. en, the following simulated signals (Figures 3(a) and 3(b)) are used to compare the robustness of the kurtosis index and the 1.5D spectral kurtosis index against nonperiodic transient impulse components. e time-domain kurtosis and 1.5D spectral kurtosis are calculated for Gaussian noise, and the results are 2.9587 and 12.1610, respectively. en, a set of nonperiodic transient impulse are added at 0.3s of the simulated signal, and the kurtosis and 1.5D spectral kurtosis values are 9.3163 and 11.9932, respectively. e results show that the kurtosis index defined in the time domain is easily interfered by nonperiodic transient impulse, and the 1.5D spectral kurtosis value is more robust to nonperiodic transient impulse.

Proposed Method
(1) Firstly, the parameters of TQWT are selected.
According to the characteristics of the frequency band where the bearing fault information is located, the number of signal decomposition layers J is selected as 6, which is beneficial to abandon the lowfrequency modulation harmonic characteristics and the high-frequency interference influence. e Q-factor is obtained through traversal optimization within a certain range, and this paper uses 0.5 as the step size to obtain the optimization between 2.5 and 5.
(2) e subbands under different Q-factors and decomposition scales are demodulated by Hilbert envelope, and the 1.5D spectrum of the demodulated signal is calculated, and then the 1.5D spectral kurtosis value is calculated for each discrete spectrum sequence. (3) e optimal subband is selected, and its time-domain signal and the corresponding 1.5D spectrum are analyzed, which can extract fault feature information.

Simulated Analysis
In order to demonstrate the process of the proposed method and verify the effectiveness of the proposed method, a group of simulated signals are used for analysis. e mathematical model of simulated signal can be expressed by · sin 2πf o t , e simulated signal consists of four parts: y 1 represents periodic impulse components, y 2 represents nonperiodic transient impulse components, y 3 represents modulation harmonic components, n(t) represents a certain degree of noise, and Y is the composite signal of four components. Where A � 2 is the amplitude of the periodic impulse components, A 1 � 3 and A 2 � 6 are the amplitudes of the nonperiodic transient impulse components, A 3 � A 4 � 0.6 are the amplitudes of the modulation harmonic components. ξ is the damping ratio, f o is the natural frequency of the system, and f outer is the characteristic frequency of the periodic impulse components. b 1 � 260, b 2 � 360, t 1 � 0.2, t 2 � 0.3, and f 1 � f 2 � 3600 Hz are damping parameters, time parameters, and natural frequencies of nonperiodic transient impulse components, respectively. n(t) is the noise (SNR � − 12 dB). e parameters of the simulated signal are listed in Table 2.  e simulated signal is analyzed using the method proposed in this paper. By adjusting the adjustable factor Q,    Figure 5. e horizontal axis represents different scales, the color of the line represents the Q-factor, and the vertical axis represents the 1.5D spectral kurtosis value. It can be seen that when the Q-factor is 5 and the scale is 5, the corresponding subband index value reaches the maximum, indicating that the subband is the optimal subband. e process of finding the optimal Q-factor and scale is equivalent to finding the optimal center frequency and bandwidth of the filter bank, which can be represented by the following equivalent filtering process ( Figure 6). e selected subband with Q-factor of 5 and scale of 5 is analyzed, and its time domain and 1.5D spectrum are shown in Figures 7(a) and 7(b). And the fault characteristic frequency and its multiplication component can be extracted.
For comparison, the kurtosis-guided TQWT method is used to select subbands. e optimal subband is corresponding to the Q-factor of 5 and scale of 1 (Figure 8(a)). e subband signal is analyzed (Figures 8(b) and 8(c)), and it can be found that the subband signal is the nonperiodic transient impulse components, and no fault characteristic information can be obtained.

Experimental Data with Nonperiodic Transient Impulse
Components. In order to further verify the effectiveness of the proposed method in processing bearing fault signals containing abnormal impulse components, the bearing vibration signal data of CWRU are analyzed. e experimental platform is shown in Figure 9. e bearing model used is JEMSKF6203-2RS. e selected experimental data are the rolling element fault data with relatively small fault size and Shock and Vibration 7 weak fault degree. e driving motor speed is 1750 rpm, and the sampling frequency is 12000 Hz. rough the above information, combined with the calculation equations of bearing fault characteristic frequency, the theoretical fault characteristic frequency under different types of bearing faults can be obtained, as shown in Table 3. e time-domain waveform of the intercepted experimental data is shown in Figure 10. It can be seen from the figure that this segment of data has obvious nonperiodic transient impulse components. e original signal is processed by the method proposed in this paper, and the results are shown in Figure 11. It can be seen from Figure 11(a) that the optimal subband is selected when the Q-factor is 5 and the scale is 2.
en, the subband is analyzed. e time-domain waveform is shown in Figure 11(b), which contains obvious periodic impulse components. e 1.5D spectrum is shown in Figure 11(c). e characteristic frequency of rolling element fault is extracted. is shows that the proposed method is robust to the nonperiodic transient impulse components.
As a comparison, the kurtosis index is used to select the subband signals under different Q-factors and scales. e optimal subband is shown in Figure 12(a). e time-domain waveform of the subband is shown in Figure 12(b). It can be seen that its main component is nonperiodic transient impulse, and its 1.5D spectrum is shown in Figure 12(c), and no obvious fault feature information can be found.  Shock and Vibration

Experimental Data with Low Signal-to-Noise Ratio.
In order to verify the processing ability of the proposed method for low signal-to-noise ratio signal, the data of QPZZ-II bearing test-bed are used for analysis and verification. e structure of the experimental platform is shown in Figure 13. e power source of the test-bed is the variable frequency drive motor. e motor is connected with the rotor through the coupling. e two ends of the rotor are bearing seats. e normal bearing is in the left bearing seat, and the fault bearing is in the right bearing seat. e structural parameters of the bearing are shown in Table 4. e bearing model is 6205E. e speed of driving motor is 1470 rpm. e theoretical fault characteristic frequency under different types of bearing faults can be obtained, as shown in Table 5.
e time-domain waveform of the experimental signal is shown in Figure 14. It can be seen from the time-domain waveform that the signal-to-noise ratio of the vibration    signal is low, and the periodic impulse characteristics are submerged. e original signal is analyzed by the kurtogram algorithm, and the analysis results are shown in Figure 15. It can be seen from Figure 15(a) that the optimal frequency band is in the low-frequency band (the frequency band is actually noise band, and the frequency band where the periodic impulse characteristics of bearing fault is located is about 2000 Hz, which is submerged). e time-domain waveform and envelope spectrum of the optimal subband signal are shown in Figures 15(b) and 15(c), respectively. e method fails to extract the effective fault information. e original signal is processed by the method proposed in this paper, and the results are shown in Figure 16. It can be seen from Figure 16(a) that the optimal subband is selected when the Q-factor is 3 and the scale is 5. en, the subband is analyzed. e time-domain waveform is shown in Figure 16(b). e 1.5D     spectrum is shown in Figure 16(c). e characteristic frequency of rolling element fault is extracted. is shows that the proposed method is robust to low signal-to-noise ratio signal.

Conclusions
e core problem in bearing signal processing is how to correctly select the frequency band of the fault information, so as to extract the information within the frequency band and filter out the interference information outside the frequency band. TQWT can provide a signal frequency band division method, and because its Q-factor is adjustable, it can adjust the wavelet oscillation characteristics to adapt to various signal characteristics. In order to improve the accuracy of subband selection and avoid the noise and nonperiodic transient impulse interference in the environment, this paper proposes a Q-factor and scale optimization method of TQWT with 1.5D spectral kurtosis as the evaluation index. Simulated and experimental results show that the proposed method has good analysis effect when SNR is low and when there is nonperiodic transient impulse interference.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.