Next Article in Journal
Creation of Artificial Luciferase 60s from Sequential Insights and Their Applications to Bioassays
Next Article in Special Issue
Investigating Machine Learning and Control Theory Approaches for Process Fault Detection: A Comparative Study of KPCA and the Observer-Based Method
Previous Article in Journal
Concentration Retrieval in a Calibration-Free Wavelength Modulation Spectroscopy System Using Particle Swarm Optimization Algorithm
Previous Article in Special Issue
Real-Time Forecasting of Subsurface Inclusion Defects for Continuous Casting Slabs: A Data-Driven Comparative Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Denoising Method for Hydroelectric Unit Vibration Signal Based on ICEEMDAN–PE–SVD

1
Intelligent Power Equipment Technology Research Center, Wuhan University, Wuhan 430072, China
2
School of Power and Mechanical Engineering, Wuhan University, Wuhan 430072, China
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(14), 6368; https://doi.org/10.3390/s23146368
Submission received: 17 June 2023 / Revised: 7 July 2023 / Accepted: 11 July 2023 / Published: 13 July 2023

Abstract

:
Vibration monitoring and analysis play a crucial role in the fault diagnosis of hydroelectric units. However, accurate extraction and identification of fault features from vibration signals are challenging because of noise interference. To address this issue, this study proposes a novel denoising method for vibration signals based on improved complementary ensemble empirical mode decomposition with adaptive noise (ICEEMDAN), permutation entropy (PE), and singular value decomposition (SVD). The proposed method is applied for the analysis of hydroelectric unit sway monitoring. Firstly, the ICEEMDAN method is employed to process the signal and obtain several intrinsic mode functions (IMFs), and then the PE values of each IMF are calculated. Subsequently, based on a predefined threshold of PE, appropriate IMFs are selected for reconstruction, achieving the first denoising effect. Then, the SVD is applied to the signal after the first denoising effect, resulting in the SVD spectrum. Finally, according to the principle of the SVD spectrum and the variation in the singular value and its energy value, the signal is reconstructed by choosing the appropriate reconstruction order to achieve the secondary noise reduction effect. In the simulation and case analysis, the method is better than the commonly used wavelet threshold, SVD, CEEMDAN–PE, and ICEEMDAN–PE, with a signal-to-noise ratio (SNR) improvement of 6.9870 dB, 4.6789 dB, 8.9871 dB, and 4.3762 dB, respectively, and where the root-mean-square error (RMSE) is reduced by 0.1426, 0.0824, 0.2093 and 0.0756, respectively, meaning that our method has a better denoising effect and provides a new way for denoising the vibration signal of hydropower units.

1. Introduction

Water turbine generator units play a crucial role in the global energy transformation towards green, low-carbon, and sustainable development. With the development of the economy and society, the requirements for the safety, reliability, and stability of water turbine generator units continue to increase. However, these units operate in a complex environment and face various challenges that pose significant threats to their safe and stable operation. Therefore, effective monitoring and diagnosis of the unit’s operating condition have become a research focus in the industry [1]. During the long-term operation of hydropower units, various faults and abnormalities may occur due to the impact of water flow, load changes, mechanical wear and tear, and other factors. These faults and abnormalities often lead to changes in the vibration signal of the unit. By analyzing and processing the vibration signal, key features regarding the unit status and fault information can be extracted for status evaluation and fault diagnosis of the unit. Studies have shown that nearly 80% of faults in water turbine generator units can be reflected in vibration signals [2].
Vibration signals, as an important monitoring indicator for assessing the operating status of units [3], can reflect the internal operation and fault information of the units. Through the analysis and processing of vibration signals, equipment anomalies can be detected promptly, potential failures can be predicted, downtime can be reduced, and the availability and operational efficiency of the units can be improved, ensuring their safe and stable operation [4]. Therefore, vibration signal analysis holds significant value in the fault diagnosis of water turbine generator units. However, the analysis, processing, and feature extraction of vibration signals from these units face certain challenges. Firstly, the operating status of water turbine generator units is dynamic, and vibration signals exhibit nonlinearity and nonstationarity. The frequency spectrum and amplitude vary with time, presenting complex time–frequency characteristics. This complexity makes the accurate analysis and diagnosis of vibration signals more difficult. Secondly, due to the complexity of the operating environment and the limitations of sensors in water turbine generator units, vibration signals are often subjected to various noise interferences such as mechanical noise, electromagnetic interference, and environmental noise. The presence of noise can mask or distort useful fault features, reducing the reliability of fault information in vibration signals and posing challenges for feature extraction and identification.
The accurate and reliable extraction and identification of fault characteristics in vibration signals are essential for the normal operation and maintenance of hydroelectric units. Vibration signal denoising is of great significance in the condition evaluation and fault diagnosis of hydropower units. Firstly, through the denoising process, noise interference can be effectively reduced and the quality and availability of vibration signals can be improved. Secondly, the denoising of the vibration signal of the hydropower unit can extract the characteristic information related to the operation status and fault of the unit. The vibration signal contains rich information on mechanical features, dynamic characteristics, and harmonic components, which can be used to evaluate the working status of the unit, detect abnormal operations, and predict potential faults. By denoising the vibration signal, the characteristic information can be highlighted to help engineers and technicians more accurately determine the status of the unit and take the corresponding maintenance and repair measures in time to improve the reliability and operational efficiency of the unit. Therefore, researchers have carried out a lot of research work in the field of vibration signal denoising, aiming to filter out noise interference and improve the effective extraction of real signal features.
Currently, commonly used methods for vibration signal denoising include Fourier transform, wavelet transform, empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD), complementary ensemble empirical mode decomposition (CEEMD), complementary ensemble empirical mode decomposition with adaptive noise (CEEMDAN), variational mode decomposition (VMD), and singular value decomposition (SVD) [5,6,7,8,9]. Fourier transform is a classical linear signal analysis method suitable for analyzing stationary, regular linear signals. However, for nonstationary and nonlinear vibration signals, the Fourier transform performs poorly in extracting fault features. Wavelet transform is a method capable of analyzing nonstationary signals and allows for the analysis of signals in the time–frequency domain [6]. However, the parameter settings of the wavelet transform lack adaptability, as different signals require different wavelet basis functions, limiting its application. Empirical mode decomposition (EMD) is a relatively new time–frequency analysis algorithm that decomposes signals into several intrinsic mode functions (IMFs), enabling comprehensive analysis in the time–frequency domain. However, EMD tends to produce mode mixing and endpoint effects during signal decomposition, which limits its application in vibration signal denoising [7]. EEMD is an improved method based on EMD [8]. Although it can partially suppress mode mixing and endpoint effects caused by EMD, this method itself introduces new issues such as residual noise, decreasing the denoising effectiveness [10,11,12]. To address this, Yeh and Huang [13] proposed the CEEMD method, which effectively suppresses residual noise. However, the method has not fully addressed issues such as incomplete signal decomposition and low computational efficiency. VMD [14,15,16] effectively avoids endpoint effects seen in similar EMD methods but requires presetting the decomposition parameter K, making it difficult to achieve adaptive decomposition that meets the requirements of online automatic monitoring and analysis of vibration signals from water turbine generator units. CEEMDAN is an improved empirical mode decomposition method [17,18,19] that has certain advantages in denoising nonstationary signals. It can suppress mode mixing and endpoint effects and improve denoising effectiveness through an adaptive noise model [5,7,9]. However, it still faces issues such as residual noise and pseudo-modes [20,21,22]. Marcelo A. Colominas [23] proposed an improved denoising method called improved complementary ensemble empirical mode decomposition with adaptive noise (ICEEMDAN), which is based on CEEMDAN and addresses some of its limitations [20,21,22,24]. Unlike CEEMDAN, ICEEMDAN incorporates white noise as part of the complete noise ensemble instead of directly adding Gaussian white noise. SVD is a matrix decomposition method [25,26,27] that decomposes and transforms matrices, allowing the collected signals to be decomposed into a series of superimposed linear components. It can effectively detect subtle information variations in signals under complex backgrounds and is widely used in denoising and feature extraction [28,29,30,31,32]. In recent years, nonstationary vibration signal noise reduction via CEEMDAN–PE, ICEEMDAN–PE, and CEEMDAN–SVD has been used for the vibration feature extraction of hydropower units [7], underwater acoustic monitoring denoising [9], ball mill barrel vibration signal denoising [18], heart and brain electrical signal noise filtering [19], weld signal denoising [21], coal machine and gas signal denoising [33,34], building structure vibration signal denoising [26,35], parameter prediction [22,24,36,37] and gear-bearing fault diagnosis [31,38,39,40,41,42,43]. The initial applications of these methods have made good progress, but there is still a need to improve the vibration signal denoising of hydropower units.
To overcome the pseudo-modal problem in the CEEMDANethod and the limitations of the SVD method including matrix size limitation and information loss, as well as the unsatisfactory effect of the single denoising method of ICEEMDAN and SVD in the denoising of hydropower unit vibration signals, this study proposes a method that combines ICEEMDAN with PE and SVD for noise reduction in hydropower unit vibration signals, aiming to improve the accuracy and reliability of vibration signals and provide more effective technical support for fault diagnosis and predictions regarding hydropower units. The method processes the signal by ICEEMDAN, obtains the effective modal components (IMF), and calculates the PE value of each component. The IMF components are selected to be reconstructed according to the preset PE threshold to achieve the primary noise reduction effect [24,44,45]. Subsequently, the signal after noise reduction is processed using SVD decomposition to obtain the singular value difference spectrum, and the second noise reduction in the signal is carried out by selecting the appropriate reconstruction order; it is necessary to consider the variation in singular values and energy values, which can effectively avoid information loss to further improve the denoising effect. The results show that the proposed ICEEMDAN–PE–SVD denoising method has a higher signal-to-noise ratio and smaller root-mean-square error. Compared with the commonly used methods, this method can better retain the detailed information of the signal, effectively filter out the noise interference, and carry out more accurate and reliable fault feature extraction. Therefore, the research in this paper provides a new method for noise reduction in the vibration signal of hydropower units. In this paper, the theoretical basis, algorithm flow, simulation analysis, and example analysis of the method will be introduced in detail, which is followed by a discussion of the results.

2. Related Theory and Methods

2.1. ICEEMDAN Algorithm

The ICEEMDAN signal processing method is an improvement over CEEMDAN [20,21,22,23,24]. Unlike the traditional EMD method, ICEEMDAN incorporates adaptive noise and a complete ensemble strategy to enhance the stability and accuracy of the decomposition process. The improved method differs from CEEMDAN in that it selects the component of white noise obtained from the EMD decomposition instead of directly adding Gaussian white noise during the decomposition process. The flowchart of the ICEEMDAN algorithm is shown in Figure 1.
Let X ( t ) denote the original signal to be decomposed. E k ( ) represents the Kth I M F (intrinsic mode function) obtained through EMD decomposition. M ( ) denotes the calculation of the envelope of the reconstructed signal, which yields the local mean of the reconstructed signal. ω ( i ) ( t ) is the ith ( i = 1 , 2 , 3 , , I ) Gaussian white noise value following a standard normal distribution N ( 0 , 1 ) (mean = 0, unit variance). r k ( t ) represents the residual of the kth stage. The coefficient β k represents the signal-to-noise ratio of the kth stage. When k = 0 β 0 = ε 0 σ ( x ( t ) ) σ ( E 1 ( ω ( i ) ( t ) ) ) , and when k 1 , β k = ε 0 σ ( r k ( t ) ) (where ε 0 is a predetermined amplitude parameter and σ ( ) is the standard deviation operator). The notation denotes the operator for calculating the average value [21,22].
The steps of the ICEEMDAN decomposition algorithm are as follows [23]:
Step 1. Add white noise to the original signal X ( t ) according to Equation (1):
X ( i ) ( t ) = X ( t ) + β 0 E 1 ( ω ( i ) ( t ) )
Step 2. Decompose the obtained signal X ( i ) ( t ) using the EMD method, calculate the local mean M ( X ( i ) ( t ) ) of each component as described in Equation (2), and obtain r 1 ( t ) (the residual of the k = 1 stage):
r 1 ( t ) = M ( X ( i ) ( t ) )
Step 3. Calculate the k = 1 component of the original signal I M F 1 ( t ) as described in Equation (3):
I M F 1 ( t ) = X ( t ) r 1 ( t )
Step 4. Add white noise to the residual of the first stage r 1 ( t ) , calculate the local mean of the signal, and obtain r 2 ( t ) (the residual of the k = 2 stage) using Equation (4):
r 2 ( t ) = M ( r 1 ( t ) + β 1 E 2 ( ω ( i ) ( t ) ) )
Step 5. Calculate the k = 2 component of the original signal I M F 2 ( t ) as described in Equation (5):
I M F 2 ( t ) = r 1 ( t ) r 2 ( t )
Step 6. When k = 3 , 4 , 5 , K , calculate r k ( t ) (the residual of the kth stage) using Equation (6):
r k ( t ) = M ( r k 1 ( t ) + β k 1 E k ( ω ( i ) ( t ) ) )
Step 7. Calculate the kth component of the original signal as described in Equation (7):
I M F k ( t ) = r k - 1 ( t ) - r k ( t )
Step 8. Return to Step 6 and calculate the next k.
Step 9. The original signal, after being decomposed by ICEEMDAN, can be represented as:
X ( t ) = k = 1 K I M F k ( t ) + r k ( t )

2.2. PE Algorithm

PE is a nonlinear analysis method used to detect randomness and dynamic transitions in time series data [19,24,44,45]. It is an indicator commonly used to measure the complexity and nonlinear characteristics of signals and is widely employed in the analysis of time series and complex systems. The calculation of PE is based on the arrangement of the signal, where the signal sequence is sorted to obtain different permutations. The occurrence frequency of each permutation is then counted, and the probability distribution of each permutation is calculated. PE reflects the irregularity of the signal over time, where higher values indicate greater complexity and nonlinear features in the signal. It is highly sensitive to signal transitions and can quantitatively assess the presence of random noise in a signal sequence. The calculation process of PE is illustrated in Figure 2.
Let X ( t ) be a time series of length N: { X ( t ) , t = 1 , 2 , 3 . , N } , m is the embedding dimension,  τ is the time delay, and K = N ( m 1 ) τ is the number of reconstructed components. The calculation steps for PE are as follows [19,24,44,45]:
Step 1. Perform phase-space reconstruction on the time series X ( t ) to obtain the m-dimensional delay vector sequence, denoted as matrix Z:
Z = x ( 1 ) x ( 1 + τ ) x ( 1 + ( m 1 ) τ ) x ( j ) x ( j + τ ) x ( j + ( m 1 ) τ x ( K ) x ( K + τ ) x ( K + ( m 1 ) τ
Step 2. Each row of the matrix represents a reconstructed component. Rearrange each component ( x ( j ) , x ( j + τ ) , , x ( j + ( m 1 ) τ ) in matrix Z in ascending order according to their values:
x ( t + ( j 1 1 ) τ ) x ( t + ( j 2 1 ) τ ) x ( t + ( j m 1 ) τ )
where t represents the column index, and j 1 , j 2 , , j m represent the positions of each element in the reconstructed component.
Thus, for any time series X ( t ) , the matrix Z obtained through reconstruction will have a set of symbolic sequences.
S ( l ) = ( j 1 , j 2 , , j m )
where l = 1 , 2 , , k , and 1 k m ! , the m-dimensional phase space maps different sequences of symbols with a total of m ! possibilities of permutations, and S ( l ) is one of them.
Step 3. Calculate the probability of occurrence for each symbolic sequence, denoted as P 1 , P 2 , , P k . The probability is calculated as the count of occurrences of each symbolic sequence divided by m ! (the total count of different symbolic sequences). Define the PE as the Shannon entropy of the probability distribution, which measures the uncertainty or randomness of the time series:
H P ( m ) = j = 1 k P j ln ( P j )
To facilitate comparison, the PE is often normalized by dividing it by the maximum value it can achieve, denoted as H P ( m ) max = ln ( m ! ) . This normalization yields a normalized PE H P :
H P = H P ( m ) ln ( m ! )
where H P 0 , 1 . The value of H P reflects the randomness or irregularity of the time series, and it is positively correlated with the degree of randomness. The variation in H P amplifies the subtle details of the time series.

2.3. SVD Algorithm

SVD is a matrix decomposition method that decomposes a matrix into the product of three matrices: an orthogonal matrix, a diagonal matrix, and the transpose of another orthogonal matrix. SVD is commonly used in signal processing for tasks such as data dimensionality reduction, signal denoising, and feature extraction [25,26,27,28,29,30,31,32]. In signal denoising, SVD aims to eliminate noise by constructing a matrix that contains the signal information and then decomposing this matrix into a series of singular values and corresponding singular vectors representing the time–frequency subspaces. This approach helps retain useful information related to the faulty signal and finds extensive applications in signal analysis and processing [43].
Let X ( t ) be a time series of length N: { X ( t ) , t = 1 , 2 , 3 . , N } . The specific steps of SVD are as follows:
Step 1. Phase-space reconstruction is carried out by constructing an m × n order Hankel matrix ( H ) for the signal to be decomposed.
H = x ( 1 ) x ( 2 ) x ( n ) x ( 2 ) x ( 3 ) x ( n + 1 ) x ( m ) x ( m + 1 ) x ( N )
where in Equation (14),  N = m + n 1 N > m n > 1 , and H R m × n . The number of rows,  m , and the number of columns,  n , can be determined based on the following principles: when the length N is even, we set m = N 2 + 1 and n = N 2 ; when the length N is odd, we have n = m = N + 1 2 .
Step 2. Perform SVD decomposition on the matrix H .
H = U S V T = U Λ 0 0 0 V T
In Equation (15),  U R m × m and V R n × n are orthogonal matrices, representing the left and right singular vectors, respectively. S R m × n is a diagonal matrix with its diagonal elements being the singular values of H . Let Λ denote the diagonal matrix formed by arranging the singular values λ i in descending order, given by Λ = diag ( λ 1 , λ 2 , , λ q ) , where q = min ( m , n ) and λ 1 λ 2 λ q .
Step 3. Compute the differences between adjacent singular values, resulting in a sequence of differential singular values b i . Choosing the appropriate order of singular values is crucial for signal reconstruction, and the differential singular value spectrum effectively captures the singular value variations in noisy signals.
b i = λ i λ i + 1
In Equation (16),  i = 1 , 2 , 3 , , q and λ q + 1 = 0 .
k exists such that b k is the peak point and there is no abrupt change in the singular value after the values of k and b k gradually approach 0. In the case of the singular value sequence, the maximum peak represents the boundary between the noise and useful signal. When k = 1 , it indicates the presence of a DC component in the signal; we assign the second largest peak to k . On the other hand, when k > 1 , it suggests that the signal does not contain a DC component or that the DC component is smaller than the AC component [46]. In a noisy signal segment, the useful signal only contributes to the first singular values of the reconstructed matrix, while the noise uniformly contributes to the singular values [29].
Step 4. After determining the value of k , the singular values beyond are set to zero, resulting in a diagonal matrix where Λ 0 = diag ( λ 1 , λ 2 , , λ k , , λ q ) λ k + 1 = λ k + 2 = = λ q = 0 , and S 0 = Λ 0 0 0 0 , which is used to reconstruct the signal through SVD [28]. By analyzing the noise test results, we can obtain the denoised signal represented by the matrix H 0 .
H 0 = U S 0 V T = U Λ 0 0 0 0 V T

3. ICEEMDAN–PE–SVD-Based Vibration Signal Denoising for Hydropower Units

The proposed denoising method, ICEEMDAN–PE–SVD, combines the advantages of the aforementioned three algorithms to achieve dual denoising of turbine vibration signals under strong noise and complex electromagnetic interference backgrounds. The flowchart of the denoising method based on ICEEMDAN–PE–SVD for hydroelectric turbine vibration signals is shown in Figure 3. The detailed steps are shown in Algorithm 1.
Algorithm 1. Denoising method based on ICEEMDAN–PE–SVD for hydroelectric turbine vibration signals.
Input: Original signal x ( t ) .
Output: Denoised signal x ( t ) .
1: According to Equations (1)–(8), the ICEEMDAN decomposition of the original signal x ( t ) is calculated to obtain a set of intrinsic mode functions ( I M F k ( t ) ) and a residue component ( r k ( t ) ). X ( t ) k = 1 K I M F k ( t ) + r k ( t )
2: The PE of the resulting IMF components is calculated according to Equations (9)–(13). In [34], Brandt et al., after a lot of experiments and projections, recommended that the statistical results have high reasonableness when the embedding dimension m is taken from 3 to 7, and the delay time has less influence on the calculation of the PE. Therefore, in this paper, we chose the number of embedding bits of m = 5 , and the delay time of τ = 1 . I M F k H P ( I M F k )
3: The normalized PE threshold is set to 0.3 according to the results of multiple simulation experiments and combined with the PE calculation principle.
for  I M F k in I M F K do
if  H P ( I M F k ) 0.3 I M F k is selected as the valid IMF component
return all the valid IMF components.
4: The signal S i g i c e e m d a n p e is obtained by reconstructing it according to all the valid IMF components. I M F S i g i c e e m d a n p e
5: The SVD decomposition of S i g i c e e m d a n p e occurs according to Equations (14) and (15).
6: Calculate the difference spectrum of singular values obtained from the decomposition according to Equation (16).
7: Combined with the variation trend of the difference spectrum, an appropriate singular value order k is selected to reconstruct the characteristic matrix H 0 , and then the matrix H 0 is converted to the denoised signal x ( t ) .
8: return  x ( t ) . This achieves a double denoising effect and forms the final denoised signal, finishing the denoising process of the vibration signal of the hydropower unit.

4. Simulation Analysis

4.1. Construction of the Simulation Signal

To verify the effectiveness of the vibration signal denoising method based on ICEEMDAN–PE–SVD, the oscillation signal of hydropower units was selected for simulation analysis. The pendulum of the hydropower unit is mainly affected by mechanical excitation and hydraulic excitation, and the mechanical excitation is generally dominated by a medium frequency (1, 2, or 3 times the rotational frequency), while the hydraulic excitation is generally dominated by a low frequency (0.2~0.45 times the rotational frequency). Choosing the rotational frequency as 1.25 Hz , we used a frequency that is 1, 2, 3, 4, 0.2, 0.3, or 0.45 times the rotational frequency of the signal superposition and simulated the mechanical vibration and hydraulic vibration at the same time, so that the simulated signal is constructed as follows [47]:
f ( t ) = i = 1 7 A i sin 2 π f i t
In Equation (18),  A 1 ~ A 7 represent amplitudes of 20   μ m 4.5   μ m 2.55   μ m 1.5   μ m 0.4   μ m 0.3   μ m , and 0.2   μ m , respectively. f 1 ~ f 7 represent frequencies of 1.25   Hz 2 × 1.25 Hz 3 × 1.25 Hz 4 × 1.25 Hz 0.2 × 1.25 Hz 0.3 × 1.25 Hz , and 0.45 × 1.25 Hz . The sampling frequency was set to 1000 Hz . The simulated original signal without noise is shown in Figure 4. The simulated noisy signal was obtained by adding Gaussian white noise with a signal-to-noise ratio of 20 dB , as shown in Figure 5.

4.2. Noise Reduction in Simulated Signals by ICEEMDAN–PE–SVD

When performing ICEEMDAN decomposition on the simulated signal, it is necessary to use reasonable settings including the noise amplitude coefficient ε 0 that is added to the original signal, the average number of times I for measuring the signal, the total number of times M for performing ICEEMDAN and the maximum allowed number of iterations M a x I t e r .
Combined with the findings from [23], after several experiments; setting the parameters as ε 0 = 0.2 I = 100 M = 1000 , and M a x I t e r = 1000 ; and then performing ICEEMDAN decomposition on the synthesized simulated signal, 10 IMF components and 1 trend component R were obtained, as shown in Figure 6. The PE of these components was calculated and the results are shown in Table 1.
According to the principle of PE combined with the results of several simulation experiments, the final threshold value of PE was set to 0.3, and the IMF components (IMF6~IMF10) with a threshold value less than 0.3 were selected and reconstructed with the residual R. To make it more convenient to see the effect of reconstructing the signal by different number of components after ICEEMDAN decomposition during simulation analysis, the signals that were reconstructed with different numbers of components after ICEEMDAN decomposition are shown in Figure 7.
From Figure 7, it can be seen that the reconstructed signals containing IMF1, IMF2, IMF3, IMF4, and IMF5 components have obvious noise interference, while the reconstructed signals with IMF6~IMF10 and R components are smoother and closer to the original simulated signals, which are consistent with the structure of the components selected by the PE threshold. Therefore, the IMF6~IMF10 components are classified as valid components, and the reconstructed waveforms and spectrograms are shown in Figure 8.
From Figure 8, we can see that the waveform and spectrum after one reconstruction of ICEEMDAN–PE show a certain trend compared with the original noise-free waveform and spectrum, but the figure does not show the simulated pendulum signal clearly, while the spectrum can also show the characteristic frequency but there is still residual noise.
Then, the reconstructed signal was processed with SVD decomposition, and the decomposed singular value difference spectrum was obtained, as shown in Figure 9 below.
As shown in Figure 9, when the singular value order k = 12 , the singular value difference spectrum and the singular value energy are relatively large; thus, a singular value of k > 12 does not produce abrupt changes. According to the principle of the singular value difference spectrum, and considering the changes in the singular value and its energy value, the final order of reconstruction chosen in this paper is k = 12 , i.e., the components with a singular value less than 12 are selected for reconstruction, as shown in Figure 10, which is the waveform and spectrum after ICEEMDAN–PE–SVD noise reduction.
As shown in Figure 10, after the multiple noise reductions by ICEEMDAN–PE–SVD, the waveforms and spectrograms can be compared with the original noise-free waveforms and spectrograms. The waveforms and spectrograms are basically the same, and the main characteristic frequencies are completely extracted. The bottom noise is completely removed, indicating that the method can effectively remove the noise and maintain the integrity of the original signal.

4.3. Comparative Analysis of Related Denoising Method Indices

To facilitate a quantitative comparison of the denoising performance of different methods, the signal-to-noise ratio (SNR), root-mean-square error (RMSE), and mean absolute error (MAE) are defined below [32,36], and the larger the SNR and the smaller the RMSE and MAE, the better the denoising effect.
(1)
Signal-to-noise ratio (SNR):
SNR = 10 lg i = 1 N X i 2 i = 1 N [ Y i X i ] 2
(2)
Root-mean-square error (RMSE):
RMSE = 1 N i = 1 N [ Y i X i ] 2
(3)
Mean absolute error (MAE):
MAE = 1 N i = 1 N Y i X i
In Equation (19) to Equation (21): N is the number of sampling points; X i is the original un-noised simulated signal; and Y i is the denoised signal.
To verify the effectiveness of the ICEEMDAN–PE–SVD denoising method proposed in this paper in the denoising of the oscillation signal of hydropower units, the wavelet threshold denoising method, the CEEMDAN–PE [19] denoising method, the ICEEMDAN denoising method, and the ICEEMDAN–PE–SVD denoising method proposed in this paper were used to denoise the simulated additive noise signal. The denoising effect of these different methods was compared and analyzed by using the denoising waveforms and denoising performance indices.
The selection of parameters in the wavelet threshold denoising process has a great influence on the denoising effect of the original signal. In this paper, according to the method of determining wavelet threshold denoising parameters in [18], the wavelet threshold denoising parameters of the simulated noise-added signal are set as follows: ‘sym10′ as the wavelet basis function, a decomposition layer of 4, ’heursure’ as the threshold guidelines, and a hard threshold function. The waveforms after denoising are shown in Figure 11 and Figure 12, and the performance indices of the denoising effect are shown in Table 2.
By comparing Figure 11 and Figure 12 with the simulated un-noised signal in Figure 4, it is found that wavelet thresholding and CEEMDAN denoising can filter out most of the noise in the simulated denoised signal, but the burr phenomenon in the denoised signal waveform is more obvious and the denoising effect is less satisfactory.
Comparing Figure 8 with the waveform in Figure 4, it can be found that the ICEEMDAN–PE denoising signal is achieved by denoising the simulated additive noise signal. As can be seen from Figure 10, ICEEMDAN–PE–SVD denoising is a process of noise reduction in the mixed signal-noise components discarded by the ICEEMDAN–PE filtering and denoising method, so that some effective features of the original signal can be retained while keeping the waveform smooth.
According to the denoising performance evaluation criteria and by comparing the data in Table 2, it can be seen that the ICEEMDAN–PE–SVD denoising method has the best denoising effect, and compared with the wavelet threshold, SVD, CEEMDAN–PE, and the ICEEMDAN–PE denoising method, the signal-to-noise ratio (SNR) is improved by 6.9870 dB, 4.6789 dB, 8.9871 dB, and 4.3762 dB, respectively, and the root-mean-square error (RMSE) is reduced by 0.1425, 0.0824, 0.2092 and 0.0755, respectively.
In summary, by comparing and analyzing the waveform characteristics and denoising performance evaluation indices of denoised signals of different denoising methods, it was found that the combined ICEEMDAN–PE–SVD denoising method is significantly better than the other four methods.

5. Case Analysis

To verify the effectiveness of the proposed ICEEMDAN–PE–SVD denoising method in the process of denoising the vibration signal of a hydropower unit, an example of the pendulum monitoring of unit 1 of a hydropower station was further selected for analysis. The turbine model is HLA773a-LJ-200, with a rated head of 68 m. The Online monitoring device for unit vibration and swing is shown in Figure 13, which has an integrated NI high-speed acquisition card and an eddy current displacement sensor, CWY-DO-20Q08-50V.
The speed of the unit in the process of hydro-generator pendulum signal acquisition was 250 r / min , the sampling frequency f s of the test was 2048 Hz , the water head was 65.8m and 6000 points of the sampled data were intercepted for analysis to ensure that the characteristic parameters in the signal processing analysis can fully and truly reflect the actual working conditions.
The measured signal waveform of the X-direction swing of the upper guide of the hydroelectric unit is shown in Figure 14, which contains a large number of burrs. The ICEEMDAN decomposition was first performed on this measured signal to obtain 10 IMF components and 1 trend component R. The PE value of each component was calculated, and the results are shown in Table 3. Similarly, the six components (IMF5-IM10) with a PE value less than 0.3 were selected as the effective information components for reconstruction, and then the SVD algorithm (singular value order k = 7 ) performs the secondary noise reduction to complete the denoising of the measured signal of the upper conductivity of the hydropower unit. It can be seen from Figure 15 that the noise in the measured signal of the upper guidance swing degree has been effectively filtered.
Since the original, pure pendulum signal value cannot be obtained from the measured signal and the SNR and RMSE cannot be calculated directly, the noise rejection ratio (NRR) before and after signal denoising is defined below to characterize the prominence of the effective signal after denoising. The larger the value, the more prominent the effective signal after denoising [15].
N R R = 10 lg 10 ( σ 1 2 σ 2 2 )
In Equation (22),  σ 1 2 and σ 2 2 denote the variance in the signal before and after denoising, respectively. The noise suppression examples of the four denoising methods after processing the measured signal of the pendulum degree are shown in Table 4. From Table 4, we can see that ICEEMDAN–PE–SVD denoising has the highest NRR value and the best filtering effect, whereas the denoising with ICEEMDAN–PE and CEEMDAN–PE is in second, and wavelet threshold denoising is the worst.
The measured signal and denoised signal waveforms of the Y-direction swing of the lower guide and the X-direction swing of the water guide of the hydropower unit are shown in Figure 16, Figure 17, Figure 18 and Figure 19, respectively. ICEEMDAN decomposition was performed on these two signals to obtain the IMF component and trend component. The entropy value of each component was calculated, and the results are shown in Table 3. According to the PE threshold, the division of components was selected to reconstruct, and then the secondary denoising was performed on the reconstructed signal using the SVD algorithm, which finally completes the denoising of the measured signals of the Y-direction of the lower guide and the X-direction of the water guide of the hydropower unit, as shown in Figure 17 and Figure 19, respectively.
Table 3. PE of each component for case analysis by ICEEMDAN.
Table 3. PE of each component for case analysis by ICEEMDAN.
BSSF1F2F3F4F5F6F7F8F9F10R
UG0.9720.7480.5350.4020.2870.2250.1630.1480.1480.1390.017
LG0.9640.7840.550.3290.2690.1890.1630.150.150.135
WG0.9640.750.5360.3680.2820.2080.1770.1630.1510.1490.146
In Table 3 and Table 4, BSS denotes the bearing swing signal; UG denotes the upper guide bearing; LG denotes the lower guide bearing; and WG denotes the water guide bearing.
Table 4. The noise rejection ratio of de-noised swing signals by different methods.
Table 4. The noise rejection ratio of de-noised swing signals by different methods.
Denoising MethodUG NRR/dBLG NRR/dBWG NRR/dB
Wavelet Threshold2.15693.23874.7788
SVD11.186420.299516.5618
CEEMDAN–PE10.135612.186315.7854
ICEEMDAN–PE6.652516.869214.1723
ICEEMDAN–PE–SVD11.728620.31116.6323
From the above comparison of the figures, it can be seen that the method based on ICEEMDAN–PE–SVD can effectively remove a large amount of background noise contained in the measured signals of the lower guide and water guide pendulums of hydroelectric units. The noise rejection ratios of five different denoising methods for the down–conductor and hydro-conductor pendulum signals were also calculated. From Table 4, it can be seen that the ICEEMDAN–PE–SVD denoising method has the best denoising effect.

6. Conclusions

In this paper, an ICEEMDAN–PE–SVD-based denoising method for the vibration signals of hydroelectric units was proposed. The method can be used to denoise the simulated and three measured signals of hydropower unit swing, and the following conclusions are obtained after comparing the waveform differences and the magnitude of denoising performance indices with the methods of wavelet threshold denoising, CEEMDAN–PE denoising, and ICCEMDAN–PE denoising in comparison tests:
(1)
Through simulation tests, the ICEEMDAN–PE–SVD method proposed in this paper, after the double-noise reduction process, obtains a root-mean-square error as low as 0.1152 and the signal-to-noise ratio is improved to 42.0941, which maximizes the noise elimination while retaining the useful information within the fault signal. The method has a good noise reduction and pulse effect, and avoids modal mixing in the EMD decomposition process and the pseudo-modal problem of CEEMDAN decomposition.
(2)
Through the case analysis of the oscillation data of the measured hydro-generator unit’s upper guide in the X-direction, the lower guide in the Y-direction, and the water guide in the X-direction, it was found that the method can effectively reduce the noise of the measured unit data and extract the characteristic frequency of the vibration signal more accurately so that the cause of the unit vibration can be judged by the frequency. The denoising effect of the measured signal was better than that of the traditional method, as it can effectively filter out the noise components and provide a powerful tool for the online monitoring of equipment vibration signals.
(3)
The research results of this paper can also be widely applied to signal denoising and feature extraction of high-safety equipment in nuclear power, power grids, the petrochemical industry, and other industries.

Author Contributions

Conceptualization, F.Z., J.G., and F.Y.; methodology, F.Z. and J.G.; software, F.Z.; validation, F.Y., Y.S., and Z.L.; formal analysis, F.Z.; investigation, F.Y., Y.S., and Z.L.; resources, J.G.; data curation, F.Z. and Z.L.; writing—original draft preparation, F.Z.; writing—review and editing, J.G. and F.Y.; visualization, F.Z., Y.S., and Z.L.; supervision, J.G. and F.Y.; project administration, F.Y.; funding acquisition, J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Key Research and Development Program of China under Grant No. 2020YFB1709700, entitled “Key Technologies Research and Development of Intelligent Operation and Maintenance for High-Security Equipment”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw/processed data cannot be shared at this time. Due to the nature of this research, participants of this study did not agree for their data to be shared publicly.

Acknowledgments

The authors gratefully acknowledge the support of the National Key Research and Development Program of China (Grant No. 2020YFB1709700).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, X.; Lin, H.; Liu, Y.; Hu, B. Fusion of PCA and adaptive K-Means clustering for online fault detection of hydroelectric units. J. Electron. Meas. Instrum. 2022, 36, 260–267. [Google Scholar] [CrossRef]
  2. He, K.; Wang, W.; Jin, Y.; Li, C.; Liu, W.; Chen, Q. Research on intelligent fault diagnosis method of hydroelectric units based on CNN-SVM. Hydroelectr. Energy Sci. 2023, 41, 207–210+215. [Google Scholar] [CrossRef]
  3. Chen, X.; Xiao, Y.; Xiao, F.; Wang, Y.; Xiao, Z. Early warning method for fault analysis of hydroelectric units based on EEMD-SD vibration signal analysis. Hydroelectr. New Energy 2023, 37, 1–5. [Google Scholar] [CrossRef]
  4. Lu, N.; Zhang, G.; Liu, F.; Zhou, T.; Bian, H. Vibration fault diagnosis method of hydroelectric units based on LTSA and spectral clustering. J. Wuhan Univ. Eng. Ed. 2021, 54, 1064–1069. [Google Scholar] [CrossRef]
  5. Huang, G.; Zhuang, X.; Xie, L.; Zeng, X. MEMS gyroscope denoising algorithm based on CEEMDAN-WP-SG. J. Electron. Meas. Instrum. 2022, 36, 106–113. [Google Scholar] [CrossRef]
  6. Xu, L.; Yang, J.; Zhou, L.; Long, B. Interference fiber joint denoising method based on PE-VMD and wavelet threshold. Foreign Electron. Meas. Technol. 2022, 41, 39–46. [Google Scholar] [CrossRef]
  7. Wang, W.; He, K.; Jin, Y.; Mo, F.; Zhao, X.; Pi, J.; Xiao, Z. Vibration signal feature extraction of hydroelectric units based on CEEMDAN sample entropy and PSO-SVM. J. Wuhan Univ. Eng. Ed. 2022, 55, 1167–1175. [Google Scholar] [CrossRef]
  8. Chen, Q.; Huang, N.; Riemenschneider, S.; Xu, Y. A B-spline approach for empirical mode decompositions. Adv. Comput. Math. 2006, 24, 171–195. [Google Scholar] [CrossRef]
  9. Li, Y.; Li, Y.; Chen, X.; Yu, J.; Yang, H.; Wang, L. A new underwater acoustic signal denoising technique based on CEEMDAN, mutual information, permutation entropy, and wavelet threshold denoising. Entropy 2018, 20, 563. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, Z.; Huang, N.E. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2011, 1, 1–41. [Google Scholar] [CrossRef]
  11. Wang, Q.; Li, D. Analysis of Francis turbine vibration signal feature extraction based on EEMD. Renmin Chang. 2017, 48, 96–100. [Google Scholar] [CrossRef]
  12. Fu, X.; Yu, J.; Cui, X.; Dai, L.; Huang, J.; Ren, W. Comparative study on three empirical mode decomposition methods for vibration signal analysis of blasting. Eng. Blasting 2021, 27, 21–28. [Google Scholar] [CrossRef]
  13. Yeh, J.R.; Shieh, J.S.; Huang, N.E. Complementary ensemble empirical mode decomposition: A novel noise enhanced data analysis method. Adv. Adapt. Data Anal. 2010, 2, 135–156. [Google Scholar] [CrossRef]
  14. Dragomiretskiy, K.; Zosso, D. Variational mode decomposition. IEEE Trans. Signal Process. 2013, 62, 531–544. [Google Scholar] [CrossRef]
  15. Xiao, S.; Chen, B.; Shen, D.; Chen, H. Application of improved VMD and threshold algorithm in partial discharge denoising. J. Electron. Meas. Instrum. 2021, 35, 206–214. [Google Scholar] [CrossRef]
  16. Li, J.; Chen, Y.; Lu, C. Application of an improved variational mode decomposition algorithm in leakage location detection of water supply pipeline. Measurement 2021, 173, 108587. [Google Scholar] [CrossRef]
  17. Torres, M.E.; Colominas, M.A.; Schlotthauer, G.; Flandrin, P. A complete ensemble empirical mode decomposition with adaptive noise. In Proceedings of the IEEE International Conference on Acoustics, Prague, Czech Republic, 22–27 May 2011; pp. 4144–4147. [Google Scholar] [CrossRef]
  18. Cai, G.; Zhao, X.; Hu, X.; Huang, X.; Chen, H. Ball mill cylinder vibration signal denoising method based on CEEMDAN and wavelet threshold. Mech. Sci. Technol. 2020, 39, 1077–1085. [Google Scholar]
  19. Geng, D.; Wang, C.; Zhao, J.; Ning, Q.; Jiang, X. Research on denoising method of heart impulse signal based on CEEMDAN-PE. J. Instrum. Meas. 2019, 40, 155–161. [Google Scholar] [CrossRef]
  20. Zhao, Y.; Yan, S.; He, J.; Li, J.; Zou, X.; Qian, S. Joint denoising algorithm of HIFU echo signal based on ICEEMDAN combined with MMSVC and WT. J. Meas. Sci. Instrum. 2023, 14, 35–44. [Google Scholar]
  21. Zhai, Y.; Yang, X. Denoising algorithm for welding signal based on ICEEMDAN-ICA. Therm. Manuf. Process. 2022, 51, 96–102+106. [Google Scholar] [CrossRef]
  22. Li, C.; Yang, J.; Lian, H.; Zheng, D.; Lai, Y.; Liu, H. Combined prediction method of silicone oil dissolved gas concentration based on ICEEMDAN-IPSO-ELM. High Volt. Eng. 2023, 1–12. [Google Scholar] [CrossRef]
  23. Colominas, M.A.; Schlotthauer, G.; Torres, M.E. Improved complete ensemble EMD: A suitable tool for biomedical signal processing. Biomed. Signal Process. Control 2014, 14, 19–29. [Google Scholar] [CrossRef]
  24. Zhao, X.; Chang, X.; Li, M. Short-term wind speed prediction method based on hybrid multi-step decomposition with ICEEMDAN-PE/FE-IGWO-SVR. Mod. Electron. Tech. 2023, 46, 124–130. [Google Scholar] [CrossRef]
  25. Lei, B.; Yi, P.; Xiang, J.; Xu, W. A SVD-based signal denoising method with fitting threshold for EMAT. IEEE Access 2021, 9, 21123–21131. [Google Scholar] [CrossRef]
  26. Zhang, J.; Li, Z.; Huang, J.; Cheng, M.; Li, H. Study on Vibration-Transmission-Path Identification Method for Hydropower Houses Based on CEEMDAN-SVD-TE. Appl. Sci. 2022, 12, 7455. [Google Scholar] [CrossRef]
  27. Li, C.; Peng, T.; Zhu, Y.; Lu, S. Noise reduction method of shearer’s cutting sound signal under strong background noise. Meas. Control. 2022, 55, 783–794. [Google Scholar] [CrossRef]
  28. Wang, J.; Li, J.; Liu, Y. An improved method for determining the effective rank order of singular value decomposition denoising. J. Vib. Shock. 2014, 33, 176–180. [Google Scholar] [CrossRef]
  29. Li, H.; Chu, L.; Liu, Q.; Lu, J.; Li, F. Signal denoising method based on SG-VMD-SVD. J. Jilin Univ. Inf. Sci. Ed. 2021, 39, 158–165. [Google Scholar] [CrossRef]
  30. Ren, Y.; Liu, P.; Hu, L.; Huang, J.; Qiao, R.; Chen, H.; Li, X.; Huang, S. Research on noise reduction method of pressure pulsation signal of draft tube of hydropower unit based on ALIF-SVD. Shock. Vib. 2021, 2021, 5580319. [Google Scholar] [CrossRef]
  31. Li, H.; Liu, T.; Wu, X.; Li, S. Application of relative singular value ratio SVD in bearing fault diagnosis. J. Mech. Eng. 2021, 57, 138–149. [Google Scholar] [CrossRef]
  32. Liu, D.; Wang, X.; Huang, J.; Peng, L.; Xiao, Z. Research on vibration signal denoising method of hydropower units based on SVD. Hydropower New Energy 2018, 32, 38–43. [Google Scholar] [CrossRef]
  33. Li, C. Research on Data-Driven Cutting Mode Recognition Method for Shearer Sound Signal. Master’s Thesis, Anhui University of Science and Technology, Huainan, China, 2022. [Google Scholar] [CrossRef]
  34. Zhao, W. Research on Denoising Algorithm for Dual-Domain Gas Signal Based on CEEMDAN Fusion SVD Ratio Method. Master’s Thesis, Shandong University of Science and Technology, Qingdao, China, 2019. [Google Scholar] [CrossRef]
  35. Zhang, J.; Hou, G.; Bao, Z.; Zhang, Y.; Ma, Y. Vibration signal denoising method of spillway structure based on CEEMDAN and SVD. J. Vib. Shock. 2017, 36, 138–143. [Google Scholar] [CrossRef]
  36. Zhang, J.; Jiang, Y.; Li, X.; Luo, H.; Yin, S.; Kaynak, O. Remaining useful life prediction of lithium-ion battery with adaptive noise estimation and capacity regeneration detection. IEEE/ASME Trans. Mechatron. 2022, 28, 632–643. [Google Scholar] [CrossRef]
  37. Chang, Y.; Yang, Z.; Pan, F.; Tang, Y.; Huang, W. Ultra-short-term wind power prediction method based on CEEMDAN-PE-WPD and multi-objective optimization. Power Syst. Technol. 2022. [Google Scholar] [CrossRef]
  38. Zhang, J.; Zhang, K.; An, Y.; Luo, H.; Yin, S. An integrated multitasking intelligent bearing fault diagnosis scheme based on representation learning under imbalanced sample condition. IEEE Trans. Neural Netw. Learn. Syst. 2023, 1–12. [Google Scholar] [CrossRef]
  39. Bai, L.; Han, Z.; Ren, J. Application of CEEMDAN-PE-TFPF denoising method in gear fault diagnosis. Mach. Des. Manuf. 2020, 1, 80–83+88. [Google Scholar] [CrossRef]
  40. Zheng, J.; Bian, J. Rolling bearing fault diagnosis method based on integrated CEEMDAN-SVD and inverse frequency spectrum. J. Taiyuan Univ. Technol. 2021, 52, 495–501. [Google Scholar] [CrossRef]
  41. Gu, Y.; Zeng, L.; Zhang, M.; Li, W. Feature extraction of localized faults in gearbox based on CEEMDAN-SQI-SVD. Chin. J. Sci. Instrum. 2019, 40, 78–88. [Google Scholar] [CrossRef]
  42. Jiang, F.; Lu, J.; Liu, M.; Feng, Y. Rolling bearing fault diagnosis based on CEEMDAN and CNN-LSTM. Electron. Meas. Technol. 2023, 46, 72–77. [Google Scholar] [CrossRef]
  43. Feng, F.; Rao, G.; Si, A.; Wu, G. Research on permutation entropy algorithm and its application in vibration signal mutation detection. J. Vib. Eng. 2012, 25, 221–224. [Google Scholar] [CrossRef]
  44. Bandt, C.; Pompe, B. Permutation entropy: A natural complexity measure for time series. Phys. Rev. Lett. 2002, 88, 174102. [Google Scholar] [CrossRef] [PubMed]
  45. Luo, J.; Zhang, S.; Li, Y. Application of multi-resolution singular value decomposition in demodulation analysis of rolling bearing vibration signal. J. Vib. Eng. 2019, 32, 1114–1120. [Google Scholar] [CrossRef]
  46. Su, Y.; Wang, X.; Jin, X. Analysis of subway shunting test noise based on SVD differential spectrum denoising method. China Meas. 2019, 45, 42–45+65. [Google Scholar] [CrossRef]
  47. Liu, Y.; He, S.; Yang, Q.; Gao, B.; Liu, P.; Lei, Y. A new pendulum signal denoising method and its application. J. Vib. Test Diagn. 2019, 39, 1053–1060+1136. [Google Scholar] [CrossRef]
Figure 1. ICEEMDAN algorithm flow chart.
Figure 1. ICEEMDAN algorithm flow chart.
Sensors 23 06368 g001
Figure 2. PE algorithm flow chart.
Figure 2. PE algorithm flow chart.
Sensors 23 06368 g002
Figure 3. Process of signal denoising method based on ICEEMDAN–PE–SVD.
Figure 3. Process of signal denoising method based on ICEEMDAN–PE–SVD.
Sensors 23 06368 g003
Figure 4. Waveform and spectrum chat of simulation signal without noise.
Figure 4. Waveform and spectrum chat of simulation signal without noise.
Sensors 23 06368 g004
Figure 5. Waveform and spectrum chat of simulation signal with noise.
Figure 5. Waveform and spectrum chat of simulation signal with noise.
Sensors 23 06368 g005
Figure 6. ICEEMDAN decomposition of the simulation signals.
Figure 6. ICEEMDAN decomposition of the simulation signals.
Sensors 23 06368 g006
Figure 7. Reconstructed signals with different numbers of IMF components.
Figure 7. Reconstructed signals with different numbers of IMF components.
Sensors 23 06368 g007
Figure 8. Waveform and spectrum chat by ICEEMDAN–PE.
Figure 8. Waveform and spectrum chat by ICEEMDAN–PE.
Sensors 23 06368 g008
Figure 9. Singular value difference spectral graph of SVD. Red line means singular value; blue line means singular value difference.
Figure 9. Singular value difference spectral graph of SVD. Red line means singular value; blue line means singular value difference.
Sensors 23 06368 g009
Figure 10. Waveform and spectrum chart of ICEEMDAN–PE–SVD noise reduction.
Figure 10. Waveform and spectrum chart of ICEEMDAN–PE–SVD noise reduction.
Sensors 23 06368 g010
Figure 11. Waveform and spectrum chart after wavelet threshold noise reduction.
Figure 11. Waveform and spectrum chart after wavelet threshold noise reduction.
Sensors 23 06368 g011
Figure 12. Waveform and spectrum chart by CEEMDAN–PE.
Figure 12. Waveform and spectrum chart by CEEMDAN–PE.
Sensors 23 06368 g012
Figure 13. Online monitoring device for unit vibration and swing.
Figure 13. Online monitoring device for unit vibration and swing.
Sensors 23 06368 g013
Figure 14. Waveform and spectrum chart of real signals of upper–guide–bearing swing.
Figure 14. Waveform and spectrum chart of real signals of upper–guide–bearing swing.
Sensors 23 06368 g014
Figure 15. Waveform and spectrum chart of de-noised upper–guide–bearing swing signals by ICEEMDAN–PE–SVD (singular value order k = 7 .
Figure 15. Waveform and spectrum chart of de-noised upper–guide–bearing swing signals by ICEEMDAN–PE–SVD (singular value order k = 7 .
Sensors 23 06368 g015
Figure 16. Waveform and spectrum chart of real signals of lower–guide–bearing swing.
Figure 16. Waveform and spectrum chart of real signals of lower–guide–bearing swing.
Sensors 23 06368 g016
Figure 17. Waveform and spectrum chart of de-noised lower–guide–bearing swing signals by ICEEMDAN–PE–SVD (singular value order k = 15 ).
Figure 17. Waveform and spectrum chart of de-noised lower–guide–bearing swing signals by ICEEMDAN–PE–SVD (singular value order k = 15 ).
Sensors 23 06368 g017
Figure 18. Waveform and spectrum chart of real signals of water–guide–bearing swing.
Figure 18. Waveform and spectrum chart of real signals of water–guide–bearing swing.
Sensors 23 06368 g018
Figure 19. Waveform and spectrum chart of de-noised water–guide–bearing swing signals by ICEEMDAN–PE–SVD (singular value order k = 20 ).
Figure 19. Waveform and spectrum chart of de-noised water–guide–bearing swing signals by ICEEMDAN–PE–SVD (singular value order k = 20 ).
Sensors 23 06368 g019
Table 1. PE of each component for simulation analysis by ICEEMDAN.
Table 1. PE of each component for simulation analysis by ICEEMDAN.
IMF1IMF2IMF3IMF4IMF5IMF6IMF7IMF8IMF9IMF10R
0.93750.77520.56790.40930.3090.2120.13620.14890.14980.14360.0019
Table 2. Denoising performance of different methods.
Table 2. Denoising performance of different methods.
Denoising MethodSNR/dBRMSEMAE
Wavelet Threshold35.10710.25780.2027
SVD37.41520.19760.1424
CEEMDAN–PE33.10700.32450.2504
ICEEMDAN–PE37.71790.19080.1390
ICEEMDAN–PE–SVD42.09410.11520.0909
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, F.; Guo, J.; Yuan, F.; Shi, Y.; Li, Z. Research on Denoising Method for Hydroelectric Unit Vibration Signal Based on ICEEMDAN–PE–SVD. Sensors 2023, 23, 6368. https://doi.org/10.3390/s23146368

AMA Style

Zhang F, Guo J, Yuan F, Shi Y, Li Z. Research on Denoising Method for Hydroelectric Unit Vibration Signal Based on ICEEMDAN–PE–SVD. Sensors. 2023; 23(14):6368. https://doi.org/10.3390/s23146368

Chicago/Turabian Style

Zhang, Fangqing, Jiang Guo, Fang Yuan, Yongjie Shi, and Zhaoyang Li. 2023. "Research on Denoising Method for Hydroelectric Unit Vibration Signal Based on ICEEMDAN–PE–SVD" Sensors 23, no. 14: 6368. https://doi.org/10.3390/s23146368

APA Style

Zhang, F., Guo, J., Yuan, F., Shi, Y., & Li, Z. (2023). Research on Denoising Method for Hydroelectric Unit Vibration Signal Based on ICEEMDAN–PE–SVD. Sensors, 23(14), 6368. https://doi.org/10.3390/s23146368

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop