Article Denoising of hydropower unit vibration signal based on variational mode decomposition and approximate entropy Transactions of the Institute of Measurement and Control 016, Vol. 38(3) 8 9 Ó The Author(s) 015 Reprints and permissions: sagepub.co.u/journalspermissions.nav DOI: 10.1177/0143311559064 tim.sagepub.com Xueli An 1 and Junjie Yang Abstract A denoising method for a hydropower unit vibration signal based on variational mode decomposition (VMD) and approximate entropy is proposed. The signal was decomposed by VMD into discrete numbers of modes, then the approximate entropy of each mode was computed. These modes were reconstructed according to a preset threshold of the approximate entropy. Finally, the denoising of the hydropower unit vibration signal can be achieved. A simulation signal and real signals of hydropower unit vibration were used to verify the proposed method. The results showed that the proposed method had a good denoising performance and was better than the wavelet transform method in the signal-to-noise ratio (SNR), root mean square error (RMSE) and partial correlation index. It was ideally suited for the online denoising of the hydropower unit vibration signal. Keywords Approximate entropy, denoising, hydropower unit, variational mode decomposition, vibration signal. Introduction The vibration and pressure fluctuation of a hydropower unit are important indexes for measuring the stability of the running unit (Zhang et al., 013). In the process of vibration fault diagnosis of the hydropower units, signal processing is the ey factor. During the vibration signal acquisition, due to a variety of noise interference, the collected signal may not truly reflect the operating status of the unit (Merisz and Waligo rsi, 014; Su et al., 01; Wang et al., 013). Therefore, fast and accurate extraction of the vibration signal from the noise is significant for the timely detection of equipment failure. Chiementin et al. (01) highlighted two aspects of denoising in vibration analysis, which included revealing singularities and noise elimination to eep the useful signal. They provided a reflection on the use of discrete wavelet transform through the various parameters used during processing. Wang et al. (014) proposed a sub-band adaptive denoising method for detective impulsive signatures based on the minimum description length principle with improved normalized maximum lielihood density model. In the proposed method, the threshold is determined automatically, and the noise variance does not need to be estimated. Simulation and practical signals are used to verify the effectiveness of the method. Another effective method based on a second-generation wavelet transform for the denoising of electrocardiogram signals has been introduced by Li and Wang (013). In this method, an improved half-soft threshold based on the lifting wavelet is utilized to overcome the drawbacs of thresholds applied in the classic wavelet. Ji et al. (015) proposed a method based on identification, which uses adaptive noise cancellation to process the signal from the electromechanical inertial sensors. They built a semi-physical simulation based on the identification of sensor ontology to obtain the noise that is used as the bacground in the test. A normalized least-mean-squares method is used to abstract the sensor output that is used as an adaptive noise canceller. Cheng et al. (014) reduced the strong noise of vibration signals of mining machinery gears by using a wavelet transform correlation filter, which is used to calculate the Shannon entropy of denoised wavelet coefficients. The entropy can reflect the vibration signal complexity. Experimental results showed that this method is effective for gear vibration signal processing in strong noise. Wavelet transform (An et al., 011; Gaied, 015; Li and Wang, 013; Li et al., 014; Modarresi et al., 014; Mondal et al., 014; Su et al., 013; Zhang et al., 015) is a very common signal analysis technique. This method can analyse the time frequency characteristics of the signal and has a variable resolution. The wavelet transform can decompose the signal into 1 China Institute of Water Resources and Hydropower Research, Haidian District, Beijing, China School of Information Science and Technology, Lingnan Normal University, Zhanjiang, Guangdong Province, China Corresponding author: Xueli An, China Institute of Water Resources and Hydropower Research, Haidian District, Beijing 100038, China. Email: an_xueli@163.com
An and Yang 83 many signal components with different frequency bands (scale). Based on this feature, the signal can be filtered using this method. This is a basic premise: the wavelet transform is used to extract and reconstruct signal. However, mixing will arise between the different bands, as the window width of the wavelet is just a mean square equivalent concept. This means that the energy of a band is not entirely focused on this band, and it may lea into the adjacent bands. If the signal energy of different frequencies is similar, the method of setting some band signals to zero can be used to filter in order to ignore the energy leaage. However, if the difference of energy of different frequency signals is large, the lea of the large energy signal will overcome the adjacent band small energy signal. Empirical mode decomposition (EMD) (Luo et al., 01; Wang et al., 01) has a high frequency resolution and adaptive decomposition characteristics. However, it has drawbacs of mode mixing and end effect. For the mode mixing problem of the EMD, Dragomiretsiy and Zosso (014) proposed a new variational mode decomposition (VMD), which can analyse non-linear and non-stationary signals, and solve the defects of the EMD. Pincus proposed an approximate entropy theory (Kang et al., 013; Pincus, 1991; Zhang and Liang, 013), which uses a non-negative number to represent the complexity of a time series. The more complex the time series, the bigger the approximate entropy. Conversely, the more regular the time series, the smaller the approximate entropy. This theory has been successfully applied to the non-linear analysis of complex signal (Kang et al., 013; Pincus, 1991; Zhang and Liang, 013). In this paper, a denoising method of a hydropower unit vibration signal based on the VMD and the approximate entropy is presented. Firstly, the VMD is used to decompose a noisy vibration signal into different modes of separate spectral bands. Then the approximate entropy of each mode is computed. The modes are reconstructed according to a preset threshold of the approximate entropy, then the denoising of the hydropower unit vibration signal can be achieved. Variational mode decomposition The VMD method is a novel variation method (Dragomiretsiy and Zosso, 014), which can decompose a complex signal into discrete numbers of modes u and is centred around v. It can overcome some limits of the EMD method. The decomposition is given as follows: ( X min t d(t)+ j ) u (t) e jvt u, v pt where the P u = f, u ( =1,,., K) andv ( =1,,., K) are the set of all nodes and their centre pulsation. Using () to set the signal s absolute integarable property: L(u, v, l)= a X t d(t)+ j u (t) pt + f X u D + l, f X E u e jvt ð1þ ðþ The alternate direction method of multipliers (ADMM) is used in the VMD to produce different decomposed modes and the centre frequency during each shifting operation. The ADMM is given as: (a) ADMM concept. Initialize u 1, v1, l1, n 0 Repeat n n + 1 for = 1 : K do Update u : u n + 1 arg min v Lu n 1 + 1,, u n 1 + 1, u, u n + 1,, un K, vn 1,, vn K, ð3þ ln end for for = 1 : K do Update v : v n + 1 arg min v Lu1 n + 1,, u n K + 1, v n 1 + 1,, v n 1 + 1, v, v n + 1,, vn K, ln ð4þ end for Dual ascent: l n + 1 l n + t f X until convergence P. u n + 1 u n u n + 1! u n \e. (b) Minimization w.r.t. u. The (3) is rewritten to update the modes u, as follow: u n + 1 = arg min v R ( a t d(t)+ j pt + f X u i + l g u (t) e jvt The Parseval/Plancherel Fourier isometry under the L norm is used to find the solution of this quadratic optimization problem, as follows: ð5þ ð6þ ^u n + 1 = ^f X ^u i + ^l! 1 i6¼ 1 + a(v v ) ð7þ (c) Minimization w.r.t. v. The minimization w.r.t. v can be done by using (8) (10). ( v n + 1 = arg min t d(t)+ j ) u (t) e jvt v pt v n + 1 = arg min v Z 0 (v v ) j^u (v) j dv ð8þ ð9þ
84 Transactions of the Institute of Measurement and Control 38(3) v n + 1 = R 0 v j ^u (v) j dv R 0 j^u (v) j dv ð10þ where the new v is put at the centre of gravity of the corresponding mode s power spectrum. The obtained steps for the complete optimization of the VMD are given as: Initialize ^u 1, v1, ^l 1, n 0 Repeat for = 1 : K do Update ^u for all v 0: ^u n + 1 Update v : n n + 1 ^f Pi\ ^un i + 1 P i. ^un i + ^l n 1 + a(v v ) ð11þ v n + 1 end for Dual ascent for all v 0: R 0 v ^un + 1 (v) dv R 0 ^u n + 1 (v) dv ^l n + 1 ^l n + t ^f X until convergence: P. ^u n + 1 ^u n Approximate entropy ^u n + 1 ^u n! \e. ð1þ ð13þ The approximate entropy (Kang et al., 013; Pincus, 1991; Zhang and Liang, 013) is defined as an aggregation degree of phase vector in high-dimensional space. The complexity of time series is calculated from a multi-dimensional perspective, which includes the information of time pattern. For the time series of N points {u(i)}, the approximate entropy can be obtained through the following steps: Step 1. The series {u(i)} consists of m-dimensional vectors X(i), i.e. XðÞ= i ½ui ðþ, ui+1 ð Þ,..., ui+ ð m 1ÞŠ ð14þ where the i =1,,., N-m +1. Step. For each i value, the distance between vector X(i) and the rest vectors X(j) is calculated. dx(i), ½ X (j) Š= max = 0, 1,, m 1 j u(i + ) u(j + ) j ð15þ where the j =1,,., N m +1. Step 3. Giving a threshold valve r (r. 0), for each i value, the number of d[x(i), X(j)] \ r is counted, and the ratio C m i (r) of this number to the total number of vectors is denoted as Ci m thenumberofd½x(i), X(j)Š\r (r)= f g (N m + 1) ð16þ Step 4. The logarithm of C m i (r) is computed, and its average for all i is calculated as follows: F m (r)= 1 N m + 1 N m X+ 1 i = 1 ln C m i (r) ð17þ Step 5. Repeat the processes of Step 1 to Step 4 for m + 1; the F m + 1 (r) can be obtained. Step 6. Finally, the approximate entropy is defined as Ae(m, r, N)=F m (r) F m + 1 (r) ð18þ where the m is a preselected pattern dimension, r is the preselected similar tolerance. In the actual calculation, m = and r = (0.1 0.)SD x are usually used, where the SD x is the standard deviation of the time series {x(i), i =1,,., N}. Denoising of hydropower unit vibration signal based on the VMD and the approximate entropy The VMD can decompose a signal into an ensemble of bandlimited intrinsic mode functions, and the noise is mainly concentrated in the high-frequency components. The approximate entropy can be used to describe the complexity of the non-linear dynamics parameters of the time series quantitatively. The procedure of the denoising of hydropower unit vibration signal based on the VMD and the approximate entropy contains the following three steps: 1) Use the VMD to decompose the noisy vibration signal of the hydropower unit into a discrete number of subsignals (modes) u i. ) Calculate the approximate entropy of each mode u i, and determine the threshold A y of the approximate entropy. 3) These modes, where the approximate entropy is less than a given threshold A y, are reconstructed to complete the denoising of hydropower unit vibration signal. The signal-to-noise ratio (SNR) (Su et al., 01), the root mean square error (RMSE), and the partial correlation index (Li et al., 006) are defined to assess the denoising performance of the proposed method based on the VMD and the approximate entropy. Suppose x() is an ideal signal, ^x() is the estimated signal of x() andm is the number of samples. 1) SNR: 0 R s = 10 lgb @ P M = 1 P M = 1 1 ½x()Š C ½^x() x() Š A ð19þ
An and Yang 85 Figure 1. Simulation signal of the hydropower unit vibration. ) RMSE: 3) Partial correlation index: vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u 1 X M E r = t j^x() x() j M = 1 P M ^x()x() = 1 C loc = sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi P M ½^x()Š PM ½x()Š = 1 Simulation signal analysis = 1 ð0þ ð1þ The vibration of hydropower unit is mainly affected by mechanical and hydraulic excitation. The mechanical excitations are mainly intermediate frequencies (one times, two times and three times the unit rotational speed). The hydraulic excitations are mainly low frequencies (0. 0.45 times the unit rotational speed). In this paper, the following expression is used to simulate the vibration signal of the hydropower unit s main shaft at guide bearing: s(t)= X6 i = 1 A i sin pf i t ðþ where A 1 A 6 are 0, 4.5,.55, 1.5, 0.4 and 0.3 mm, respectively; f 1 f 6 are 1.5, 3 1.5, 3 3 1.5, 4 3 1.5, 0. 3 1.5 and 0.3 3 1.5 Hz, respectively. The sampling frequency is 1000 Hz. The simulation vibration signal of the hydropower unit is shown in Figure 1(a). A shoc signal with a period of 1.5 Hz was added to the signal and each shoc sustained for 3 ms, then a Gaussian white noise with SNR of 0 db was also added to the signal. The signal with shoc and noise is shown Figure 1(b). The VMD is used to decompose the noisy signal into 10 components (Figure ). The approximate entropy of these components is computed (Table 1). It can be seen from the table that the approximate entropies of the u 3 u 10 are significantly greater than the others. On the basis of the rule that the approximate entropy threshold is less than 0.4, the u 1 u are summed to complete the denoising of the hydropower unit vibration signal. Figure 3(a) shows the results of the denoising method based on the VMD. The wavelet transform is widely used in signal denoising. In this paper, the db4 wavelet is used to process the noisy signal as shown in Figure 1(b). For purposes of comparison, the decomposition level is also 10 layers. The common sqtwolog
86 Transactions of the Institute of Measurement and Control 38(3) Figure. Variational mode decomposition (VMD) decomposition of the vibration signal.
An and Yang 87 Figure 3. The results of simulation signals after denoising. Table 1. Approximate entropy of each component. u 1 u u 3 u 4 u 5 u 6 u 7 u 8 u 9 u 10 0.043 0.043 0.546 0.66 0.590 0.553 0.538 0.563 0.694 0.833 Table. Comparison of denoising performance of the two methods. Algorithm Signal-to-noise ratio Root mean square error Correlation coefficient VMD 3.8300 0.3346 0.9997 db4 wavelet 4.9446 0.896 0.9985 VMD, variational mode decomposition. Table 3. Approximate entropy of shaft vibration signal s components at generator upper guide bearing. u 1 u u 3 u 4 u 5 u 6 u 7 u 8 u 9 u 10 0.04 0.454 0.584 0.606 0.51 0.619 0.604 0.598 0.578 0.608 is used as the threshold in processing. The results using the db4 wavelet are shown in Figure 3(b). A comparison of the denoising performance of the two methods is shown in Table. It can be seen from the table
88 Transactions of the Institute of Measurement and Control 38(3) Figure 4. Real signals of shaft vibration signal at generator upper guide bearing of the hydropower unit. Figure 5. Denoised shaft vibration signal at generator upper guide bearing. that three indicators of the SNR, the RMSE and the correlation coefficient in the VMD-based method are better than that in the db4 wavelet method, i.e. the VMD-based method has a better denoising performance. It can be seen from Figure 3 that db4 wavelet is not great for the removal of periodic shoc noise. Case study The field test data from a large hydropower unit were used to verify denoise performance of the proposed method. The data included the main shaft vibration signal of hydropower unit at the generator upper guide bearing, the generator lower guide bearing and the turbine guide bearing. The speed of the unit was 75 rpm, the sampling frequency of the test was 1000 Hz and the number of sampling points is 6000. The real signal of shaft vibration at the generator upper guide bearing of the hydropower unit is shown in Figure 4. It can be seen from the figure that the vibration waveform contains a lot of noise and periodic shocs. The period of the shocs was the rotational frequency of the unit. The shocs
An and Yang 89 Table 4. Approximate entropy of shaft vibration signal s components at generator lower guide bearing. u 1 u u 3 u 4 u 5 u 6 u 7 u 8 u 9 u 10 0.048 0.043 0.54 0.5 0.591 0.604 0.611 0.578 0.590 0.653 Table 5. Approximate entropy of shaft vibration signal s components at turbine guide bearing. u 1 u u 3 u 4 u 5 u 6 u 7 u 8 u 9 u 10 0.06 0.115 0.609 0.597 0.64 0.598 0.581 0.576 0.640 0.710 Figure 6. Real signals of shaft vibration signal at generator lower guide bearing of the hydropower unit. resulted from the coarse main shaft of the unit. The VMD was used to decompose this signal into 10 components. The approximate entropy of each component was calculated and the results are shown in Table 3. It can be seen from the table that the approximate entropy of u 1 is significantly less than the rest. On the basis of the rule that the approximate entropy threshold is less than 0.4, the u 1 itself was the denoised signal of shaft vibration at the generator upper guide bearing of the hydropower unit. Figure 5(a) shows the results of the VMDbased denoising method. It can be seen from the figure that the VMD method can eliminate a lot of noise and periodic shocs of shaft vibration signal at the generator upper guide bearing. Figure 5(b) shows the results based on db4-wavelet denoising method. For the purposes of comparison, the decomposition level was also 10 layers. The common sqtwolog was used as the threshold in processing. It can be seen from the figure that the db4-wavelet method cannot remove the periodic shocs noise and the bacground noise of shaft vibration signal at the generator upper guide bearing well. The real shaft vibration signal of the hydropower unit at the generator lower guide bearing and the turbine guide bearing are shown in Figures 6 and 8, respectively. It can be seen from the figures that the two vibration signals contain considerable noise. The VMD was used to decompose the two signals and 10 components for each signal can be obtained. The approximate entropy of each component was calculated, and the results are shown in Tables 4 and 5. On the basis of the rule that the approximate entropy threshold is less than 0.4, the components u 1 u 3 of shaft vibration signal at generator lower guide bearing were summed to achieve the denoising of this signal (Figure 7a). The components u 1 u of the shaft vibration signal at the turbine guide bearing were summed to achieve the denoising of this signal (Figure 9a). It can be seen from the Figures 7(a) and 9(a) that the VMD-based method can remove a lot of noise from shaft vibration signal of the hydropower unit at the generator lower guide bearing and the turbine guide bearing. The results of the shaft vibration signal at the generator lower guide bearing and the turbine guide bearing using the db4-wavelet method are shown in Figures 7(b) and 9(b), respectively. For purposes of comparison, the decomposition level of the two signals was also 10 layers. The sqtwolog was used as the processing threshold. It can be seen from the two figures that the db4-wavelet method cannot remove the periodic shocs noise and the bacground noise of the shaft vibration signal at the generator lower guide bearing and the turbine guide bearing well. Conclusions This paper presents a denoising method for the hydropower unit vibration signal based on the VMD and the approximate entropy. The VMD was used to decompose a noisy vibration
90 Transactions of the Institute of Measurement and Control 38(3) Figure 7. Denoised shaft vibration signal at generator lower guide bearing. Figure 8. Real signals of shaft vibration signal at turbine guide bearing of the hydropower unit. signal into a discrete number of sub-signals (modes). The criterion of the approximate entropy threshold was used to determine the number of reconstructed modes, then the denoising of the hydropower unit vibration signal was completed. The simulation signal and the real signals of the hydropower unit vibration were used to verify the proposed method. The results showed that the VMD-based denoising method can eliminate a lot of bacground noise and periodic shocs of the main shaft vibration signal of the hydropower unit at the generator upper guide bearing, the generator lower guide bearing and the turbine guide bearing. The db4-wavelet denoising method cannot remove the periodic shocs noise and bacground noise of these three vibration signals well. This means that the denoising performance of traditional signal processing method (e.g. wavelet transform) for the unit vibration signal is unfavourable, whereas the VMD-based
An and Yang 91 Figure 9. Denoised shaft vibration signal at turbine guide bearing. method has very good denoising performance. Therefore, the proposed method is very suitable for the real-time noise cancellation of the hydropower unit vibration signals, and the real operating conditions of the hydropower unit can be obtained quicly and accurately. Declaration of conflicting interest The authors declare that there is no conflict of interest. Funding This wor was supported by the National Natural Science Foundation of China (grant number 5130958). References An X, Jiang D, Liu C, et al. (011) Wind farm power prediction based on wavelet decomposition and chaotic time series. Expert Systems with Applications 38(9): 1180 1185. Cheng G, Chen X, Shan X, et al. (014) A new method of gear fault diagnosis in strong noise based on multi-sensor information fusion. Journal of Vibration and Control doi: 10.1177/ 107754631454187. Chiementin X, Kilundu B, Rasolofondraibe L, et al. (01) Performance of wavelet denoising in vibration analysis: highlighting. Journal of Vibration and Control 18(6): 850 858. Dragomiretsiy K and Zosso D (014) Variational mode decomposition. IEEE Transactions on Signal Processing 6(3): 531 544. Gaied K (015) Wavelet-based prognosis for fault-tolerant control of induction motor with stator and speed sensor faults. Transactions of the Institute of Measurement and Control 37(1): 100 113. Ji Y, Li X and Zhao J (015) Electromechanical inertial sensors test and analysis method based on identification using adaptive noise cancellation. Transactions of the Institute of Measurement and Control 37(4): 475 484. Kang Y, Cai H and Song S (013) Study and application of complexity model for hydrological system. Journal of Hydroelectric Engineering 3(1): 5 10. Li H and Wang X (013) Detection of electrocardiogram characteristic points using lifting wavelet transform and Hilbert transform. Transactions of the Institute of Measurement and Control 35(5): 574 58. Li T, Gao L, Nie Y, et al. (006) A new adaptive direct-threshold algorithm to partial discharge data processing based on empirical mode decomposition. Proceedings of the CSEE 6(15): 9 34. Li T, Zhao D, Huang Z, et al. (014) A wavelet-based grey particle filter for self-estimating the trajectory of manoeuvring autonomous underwater vehicle. Transactions of the Institute of Measurement and Control 36(3): 31 335. Luo X, Xue Y and Wang H (01) Singular data restoration for hydraulic generator fault based on EMD-wavelet analysis. Journal of Hydroelectric Engineering 31(): 45 50. Merisz J and Waligórsi M (014) Digital vibroacoustic signal processing for combustions in Heavy-Duty Diesel engine for operational and environmental machine efficiency. Vibroengineering Procedia 3: 148 153. Modarresi J, Hooshmand R and Nosratabadi S (014) Enhancement of power measurement using a modified method based on wavelet
9 Transactions of the Institute of Measurement and Control 38(3) with preprocessing: electric arc furnace case study. Transactions of the Institute of Measurement and Control 01433114556798, first published on November 10. Mondal U, Sengupta A and Patha R (014) Servomechanism for periodic reference input: discrete wavelet transform-based repetitive controller. Transactions of the Institute of Measurement and Control 0143311456657, first published on December 19. Pincus S (1991) Approximate entropy as a measure of system complexity? Proceeding of the National Academy Sciences USA 88(6): 97 301. Su L, Nan H, Yu X, et al. (01) Application of wavelet denoising of improved threshold function to vibration signal analysis of hydroelectric units. Journal of Hydroelectric Engineering 31(3): 46 51. Su W, Zhu Y, Liu F, et al. (013) An online outlier detection method based on wavelet technique and robust RBF networ. Transactions of the Institute of Measurement and Control 35(8) 1046 1057. Wang H, Luo X, Xui Y, et al. (01) EMD and index energy-based extraction of the draft tube dynamic characteristic of hydraulic turbine. Journal of Hydroelectric Engineering 31(5): 86 91. Wang J, Pan L and Cao S (013) Wavelet transforms applied to cavitation noise analysis for hydro-turbine. Journal of Hydroelectric Engineering 3(4): 15 0. Wang Y, Xiang J, Jiang Z, et al. (014) Vibration signals denoising using minimum description length principle for detecting impulsive signatures. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 8(10): 1818 188. Zhang F, Pan L and An X (013) Statistical characteristics and maintenance alarm strategy research on stability parameter of hydraulic turbine generator unit. Journal of Hydroelectric Engineering 3(5): 69 7. Zhang X and Liang J (013) Chaotic time series prediction model of wind power based on ensemble empirical mode decompositionapproximate entropy and reservoir. Acta Physica Sinica 6(5): 1 10. Zhang X, Feng N, Wang Y, et al. (015) Acoustic emission detection of rail defect based on wavelet transform and Shannon entropy. Journal of Sound and Vibration 339: 419 43.