Seismic Random Noise Attenuation Method Based on Variational Mode Decomposition and Correlation Coefficients

Seismic data is easily affected by random noise during field data acquisition. Therefore, random noise attenuation plays an important role in seismic data processing and interpretation. According to decomposition characteristics of seismic signals by using variational mode decomposition (VMD) and the constraint conditions of correlation coefficients, this paper puts forward a method for random noise attenuation in seismic data, which is called variational mode decomposition correlation coefficients VMDC. Firstly, the original signals were decomposed into intrinsic mode functions (IMFs) with different characteristics by VMD. Then, the correlation coefficients between each IMF and the original signal were calculated. Next, based on the differences among correlation coefficients of effective signals and random noise as well as the original signals, the corresponding treatment was carried out, and the effective signals were reconstructed. Finally, the random noise attenuation was realized. After adding random noise to simple sine signals and the synthetic seismic record, the improved complementary ensemble empirical mode decomposition (ICEEMD) and VMDC were used for testing. The testing results indicate that the proposed VMDC has better random noise attenuation effects. It was also used in real-world seismic data noise attenuation. The results also show that it could effectively improve the signal-to-noise ratio (SNR) of seismic data and could provide high-quality basic data for further interpretation of seismic data.


Introduction
The seismic signal is a typical nonlinear and nonstationary signal.The seismic exploration process is affected by various factors.There are effective waves and large amounts of random noise in seismic data.Therefore, effective treatment for random noise attenuation could not only improve the signal-to-noise ratio (SNR) and quality of seismic data, but also provide benefit for further interpretation of seismic data, lithology parameter inversion and seismic attributes analyses [1].At present, there are the median filter method, f-x prediction filter method, polynomial fitting method, wavelet transform and empirical mode decomposition (EMD) method for random noise attenuation in seismic exploration.They have their own advantages and disadvantages.The median filter is a smoothing-based method.By this method, the basic frequency-domain signals tend to shift to low-frequency signals, and the high-frequency signals may be damaged [2,3].In case of relatively low SNR in high frequency, the f-x prediction filter method may easily cause severe distortion of high-frequency signals and reduce the fidelity of signals and SNR of the seismic profile [4].The polynomial fitting method requires the original seismic signals to have good continuity.The false seismic events may occur after data processing [5].In the application of a wavelet transform, the selection of generating functions and de-noising thresholds has significant impacts on the effects of random noise attenuation [6].EMD has the problems of end effect and mode mixing, which may lead to unsatisfactory effects of random noise attenuation [7].
Variational mode decomposition (VMD) is an adaptive signal processing method put forward by Dragomiretskiy.Compared with EMD, it has stronger noise-resistance ability.Moreover, it could successfully separate two harmonics with very similar frequencies, and the separating effects are not affected by the sampling frequency [8][9][10][11].Li et al. introduced the principles of VMD and proposed a lateral consistency preserved VMD method [12].Liu et al. studied the seismic time-frequency representation based on VMD [13].Li et al. proposed a hybrid de-noising method based on thresholding variational mode decomposition [14].Li et al. used VMD to analyze the depositional sequence characterization [15].Jia et al. proposed a method to improve the resolution by using generalized S-transform based on VMD [16].Zhao et al. extracted intrinsic mode functions (IMFs) based on VMD from seismic amplitudes to constrain self-organizing map facies analysis [17].Lyu et al. analyzed the discontinuities with VMD-based coherence [18].
Combining VMD with correlation coefficients, this paper developed a new method for seismic data random noise attenuation.VMD was firstly used to decompose the original signals into IMFs with different characteristics.Then, the correlation coefficients between each IMF component and the original signal were calculated.The corresponding treatment was carried out based on differences among correlation coefficients of effective signals, random noise and original signals.The effective signals were reconstructed.Finally, the random noise attenuation was achieved.The results show that the VMDC method performs well in seismic random noise attenuation.

The Improved Complementary Ensemble Empirical Mode Decomposition ICEEMD
EMD was developed by Huang et al., and is a powerful analytical tool for nonlinear nonstationary signals [19].However, it has the problems of end effect and mode mixing.ICEEMD was put forward by Tary et al., and this method can solve the above problems to some extent [20][21][22][23].The calculation steps of this method are as follows: (1) Use EMD to calculate the local mean of the i-th iteration x i = x + ε 0 w i , to get the first residual error.
(3) Calculate the second residuals and the second IMF.
In the above formula, E k (•) represents the operator generating the k-th IMF, M(•) represents the operator generating local mean of signal, x is the input signal, w i is the decomposition of white noise with zero mean unit variance, ε k is a constant greater than zero, r i is the i-th residuals and I is the number of iterations.

VMD
In order to avoid the frequency mixture issue of the EMD [19], Dragomiretskiy et al. proposed a signal decomposition method with varying scales, which is the VMD method [8].Compared with EMD, it has a solid mathematical basis and could be used to effectively solve the mode mixing problem.By VMD, the original signals could be decomposed into k band-limited signals u k with the center frequency of ω k , where k is the default decomposition scale.It is assumed that each mode function u k is a limited bandwidth near its center frequency.The adaptive decomposition of the signal is realized by searching the optimal solution of the constrained variational model.The center frequency and bandwidth of each IMF are constantly updated in the iterative solution of the variational model.According to frequency-domain characteristics of actual signals, the adaptive decomposition of the signal band could be completed and some narrow-band IMFs could be obtained.
The steps of estimating the bandwidth of u k are in the following [8]: (1) To calculate the analytic function of each u k and obtain the corresponding one-sided frequency spectrum by Hilbert transform.
(2) To adjust the estimated central spectrum by adding exponential terms and modulate the frequency spectrum of each mode into the corresponding basic frequency band.
(3) To estimate the bandwidth by Gaussian smoothness of the demodulated signal and gradient energy criterion.
By following the aforementioned steps, the obtained constrained variational problem is as follows: , where {u k } = {u 1 , u 2 , . . . ,u k } is the function of each mode.{ω k } = {ω 1 , ω 2 , . . . ,ω k } is the center frequency and is the sum of each mode.
(4) To transform the above constrained variational problem into an unconstrained variational problem by introducing the Lagrange multiplier λ(t) and two-penalty factor, the formula of the augmented Lagrange multiplier could be obtained, as follows: The alternating direction method of multipliers is applied to solve the above variational problems.The iterative optimization of u k+1 , ω k+1 k and λ k+1 could get the saddle points of the augmented Lagrange multiplier.The iteration steps are as follows: (1) To initialize u 1 , ω 1 , λ 1 , n = 0.
(3) To execute the first inner loop, update u k according to and repeat step (3) until the completion of the first loop when k = K.
(5) To execute the second inner loop and update λ according to (6) To repeat step 2 and step 5 until meeting 2 ) < ε.The whole looping will end and k IMFs could be obtained.

Correlation Coefficient Method
The Pearson correlation coefficient (PCC) is a statistical method to quantitatively measure correlations between two random variables.One of its important mathematical characteristics is that the variations of positions and scales will not cause the changes of correlation coefficients, so it is suitable for correlation evaluation of geophysical data [24][25][26].The Pearson correlation coefficient could be expressed as follows [27]: where cov(X, Y) refers to covariance of X and Y. σ X and σ Y are the standard deviations.x = E(X) and y = E(Y) are the expected values of X and Y, respectively.The bigger the absolute values of the correlation coefficients, the stronger the correlation.The closer the correlation coefficient is to 1 or −1, the stronger the correlation between X and Y.The closer the correlation coefficient is to 0, the weaker the correlation between X and Y. Generally, the correlation intensity among variables could be judged according to Table 1, as follows:  1 shows that when the correlation coefficient is greater than 0.4, X and Y have good correlation.If the correlation coefficient is less than 0.2, the correlation coefficient of X and Y is poor.When the correlation coefficient is 0.2-0.4,the correlation between X and Y is general.

Random Noise Attenuation Method Based on VMD and Correlation Coefficients
It was assumed that a random noise signal was made of a noise-free signal and random noise.The VMD algorithm was firstly used to decompose the target signals into various IMFs which may include effective signals, effective signals with partial noise, and noise signals.Formula (9) was used to calculate the correlation coefficients between each IMF and the original signal.Then, on this basis, the effective signals were reconstructed.The reconstruction principles are as follows: (1) The correlation coefficient of less than 0.2 represented that the IMF component had no correlation with the original signal, and had only the random noise, so the IMF component did not participate in signal reconstruction.
(2) The correlation coefficient of bigger than 0.4 represented that there was good correlation between the IMF component and the original signal, so the IMF component participated in signal reconstruction.
(3) The correlation coefficient between 0.2 and 0.4 represented that there was weak correlation between the IMF component and the original signal.The IMF component contained the effective signal and random noise.It was decomposed into s n = s k + n k , where s k denotes the effective signal and n k denotes the residual random noise.Then, the VMD was used for IMF treatment to get s k and n k .s k participated in the signal reconstruction.The process of random noise attenuation based on VMDC is shown in Figure 1.
(2) The correlation coefficient of bigger than 0.4 represented that there was good correlation between the IMF component and the original signal, so the IMF component participated in signal reconstruction.

Simple Signal Model Test
The simple signal was expressed as z(t) = x1(t) + x2(t), where x1(t) = sin(12 × π × t) and x2(t) = sin(32 × π × t).The number of sampling points was 1000.ICEEMD and VMD were used for signal decomposition.The results are shown in Figures 2 and 3. Developed by Colominas et al., ICEEMD has developed the EMD and could solve the problems of modal mixing in the application of EMD [20,28].As shown in Figures 2 and 3    Electronics 2018, 7, x FOR PEER REVIEW 6 of 13 The simple signal was expressed as ( )     The above results show that both ICEEMD and VMD have good decomposition effects, and VMD has higher calculation efficiency.
Electronics 2018, 7, 280 7 of 13 20% random noise was added to the above simple signal using z(t) = x1(t)+x2(t)+0.2*rand(t).Figures 4 and 5 show the decomposition results by ICEEMD and VMD, respectively.In Figure 4, some IMF components contained modal mixing, which affected the reconstructed signals after superposition of the IMFs.By VMD, the correlation coefficients between each IMF and the original signal were calculated to be 0.73, 0.70 and 0.03.According to the principles of signal reconstruction in Section 2.3, IMF3 did not participate in signal reconstruction.Figure 6 shows the reconstruction results of IMF1 and IMF2.The blue line denotes the signals without noise, and the red line denotes the reconstruction results.The correlation coefficient between the decomposition signal and the original signal was 0.99 and the root mean square error (RMSE) was only 0.0985.In terms of calculation efficiency, ICEEMD took 22.90 s while VMD only took 0.51 s.
The above results indicate that VMD still has good decomposition effects with random noise in signals.Moreover, it is more beneficial to suppress random noise with higher operating efficiency.
Electronics 2018, 7, x FOR PEER REVIEW 7 of 13 The above results show that both ICEEMD and VMD have good decomposition effects, and VMD has higher calculation efficiency.
20% random noise was added to the above simple signal using z(t) = x1(t)+x2(t)+0.2*rand(t).Figures 4 and 5 show the decomposition results by ICEEMD and VMD, respectively.In Figure 4, some IMF components contained modal mixing, which affected the reconstructed signals after superposition of the IMFs.By VMD, the correlation coefficients between each IMF and the original signal were calculated to be 0.73, 0.70 and 0.03.According to the principles of signal reconstruction in Section 2.3, IMF3 did not participate in signal reconstruction.Figure 6 shows the reconstruction results of IMF1 and IMF2.The blue line denotes the signals without noise, and the red line denotes the reconstruction results.The correlation coefficient between the decomposition signal and the original signal was 0.99 and the root mean square error (RMSE) was only 0.0985.In terms of calculation efficiency, ICEEMD took 22.90 s while VMD only took 0.51 s.
The above results indicate that VMD still has good decomposition effects with random noise in signals.Moreover, it is more beneficial to suppress random noise with higher operating efficiency.

Synthetic Seismic Records Model Test
To further test the new method's application effects in seismic data processing, the synthetic seismic records model test was conducted after adding random noise.This paper established synthetic seismic records with the sampling interval of 0.1 ms and the dominant frequency of the wavelet of 45 HZ, and the synthetic seismic records with 20% random noise, as shown in Figure 7.The comparison in Figure 7 shows that after adding noise, the SNR of the synthetic seismic records was reduced, the events became blurred and some information of seismic horizons was masked by the random noise.
The decomposition of synthetic seismic records with 20% random noise by ICEEMD could obtain the seismic records and noise profiles, as shown in Figure 8.It could be seen that there was strong random noise in the profile after suppressing the random noise.Figure 9 shows the seismic

Synthetic Seismic Records Model Test
To further test the new method's application effects in seismic data processing, the synthetic seismic records model test was conducted after adding random noise.This paper established synthetic seismic records with the sampling interval of 0.1 ms and the dominant frequency of the wavelet of 45 HZ, and the synthetic seismic records with 20% random noise, as shown in Figure 7.The comparison in Figure 7 shows that after adding noise, the SNR of the synthetic seismic records was reduced, the events became blurred and some information of seismic horizons was masked by the random noise.
The decomposition of synthetic seismic records with 20% random noise by ICEEMD could obtain the seismic records and noise profiles, as shown in Figure 8.It could be seen that there was strong random noise in the profile after suppressing the random noise.Figure 9 shows the seismic

Synthetic Seismic Records Model Test
To further test the new method's application effects in seismic data processing, the synthetic seismic records model test was conducted after adding random noise.This paper established synthetic seismic records with the sampling interval of 0.1 ms and the dominant frequency of the wavelet of 45 HZ, and the synthetic seismic records with 20% random noise, as shown in Figure 7.The comparison in Figure 7 shows that after adding noise, the SNR of the synthetic seismic records was reduced, the events became blurred and some information of seismic horizons was masked by the random noise.
The decomposition of synthetic seismic records with 20% random noise by ICEEMD could obtain the seismic records and noise profiles, as shown in Figure 8.It could be seen that there was strong random noise in the profile after suppressing the random noise.Figure 9 shows the seismic records after suppressing the random noise by the VMDC.Compared with the noise attenuation results by ICEEMD and VMDC, the de-noising effects by the VMDC are superior to the effects by ICEEMD.In Figure 9, the seismic events were better restored and the SNR was greatly improved.In addition, the noise in Figure 8 contained a small amount of effective waves, while noise was dominant in Figure 9.
The above results show that, under the synthetic seismic records with the random noise, the reconstructed events by the VMDC are more obvious and continuous, and the random noise reduction effects are better.
records after suppressing the random noise by the VMDC.Compared with the noise attenuation results by ICEEMD and VMDC, the de-noising effects by the VMDC are superior to the effects by ICEEMD.In Figure 9, the seismic events were better restored and the SNR was greatly improved.In addition, the noise in Figure 8 contained a small amount of effective waves, while noise was dominant in Figure 9.
The above results show that, under the synthetic seismic records with the random noise, the reconstructed events by the VMDC are more obvious and continuous, and the random noise reduction effects are better.Electronics 2018, 7, x FOR PEER REVIEW 9 of 13 records after suppressing the random noise by the VMDC.Compared with the noise attenuation results by ICEEMD and VMDC, the de-noising effects by the VMDC are superior to the effects by ICEEMD.In Figure 9, the seismic events were better restored and the SNR was greatly improved.In addition, the noise in Figure 8 contained a small amount of effective waves, while noise was dominant in Figure 9.
The above results show that, under the synthetic seismic records with the random noise, the reconstructed events by the VMDC are more obvious and continuous, and the random noise reduction effects are better.records after suppressing the random noise by the VMDC.Compared with the noise attenuation results by ICEEMD and VMDC, the de-noising effects by the VMDC are superior to the effects by ICEEMD.In Figure 9, the seismic events were better restored and the SNR was greatly improved.In addition, the noise in Figure 8 contained a small amount of effective waves, while noise was dominant in Figure 9.
The above results show that, under the synthetic seismic records with the random noise, the reconstructed events by the VMDC are more obvious and continuous, and the random noise reduction effects are better.

Case Study
To fully verify the application effects of the proposed VMDC in real-world seismic data, this paper selected the 3-D seismic data in Inner Mongolia in China to carry out the test.The acquisition of seismic data was in the winter, so the wind was blowing very hard.At the same time, the gangue field and air shafts were under construction.What is worse, the random noise interference was more serious in this area because there were many vehicles in adjacent industrial areas.Figure 10 shows the actual seismic profiles with random noise.It could be found that the existing random noise reduced the SNR of the seismic data and influenced the continuity of events in seismic records.
The ICEEMD and the VMDC were used for noise attenuation of post-stack seismic data.The signal reconstruction process was introduced by taking the 50th channel as an example.Firstly, ICEEMD and VMDC were used to decompose seismic signals.Then, by ICEEMD, the correlation coefficients between each IMF and the original signal were calculated as 0.1432, 0.5058 and 0.7329.By VMD, the correlation coefficients between each IMF and the original signal were calculated as 0.7221, 0.6611 and 0.1156.Therefore, the IMF2 and IMF3 decomposed by ICEEMD, and IMF1 and IMF2 decomposed by VMDC, participated in the signal constructions.The signals in other seismic traces were reconstructed by similar steps, as shown in Figures 11 and 12. Through the comparative analysis of Figures 10-12, it is found that both ICEEMD and the VMDC could suppress the random noise to a certain extent, and the latter could significantly enhance the continuity of events.

Case Study
To fully verify the application effects of the proposed VMDC in real-world seismic data, this paper selected the 3-D seismic data in Inner Mongolia in China to carry out the test.The acquisition of seismic data was in the winter, so the wind was blowing very hard.At the same time, the gangue field and air shafts were under construction.What is worse, the random noise interference was more serious in this area because there were many vehicles in adjacent industrial areas.Figure 10 shows the actual seismic profiles with random noise.It could be found that the existing random noise reduced the SNR of the seismic data and influenced the continuity of events in seismic records.
The ICEEMD and the VMDC were used for noise attenuation of post-stack seismic data.The signal reconstruction process was introduced by taking the 50th channel as an example.Firstly, ICEEMD and VMDC were used to decompose seismic signals.Then, by ICEEMD, the correlation coefficients between each IMF and the original signal were calculated as 0.1432, 0.5058 and 0.7329.By VMD, the correlation coefficients between each IMF and the original signal were calculated as 0.7221, 0.6611 and 0.1156.Therefore, the IMF2 and IMF3 decomposed by ICEEMD, and IMF1 and IMF2 decomposed by VMDC, participated in the signal constructions.The signals in other seismic traces were reconstructed by similar steps, as shown in Figures 11 and 12. Through the comparative analysis of Figure 10, Figure 11 and Figure 12, it is found that both ICEEMD and the VMDC could suppress the random noise to a certain extent, and the latter could significantly enhance the continuity of events.

Case Study
To fully verify the application effects of the proposed VMDC in real-world seismic data, this paper selected the 3-D seismic data in Inner Mongolia in China to carry out the test.The acquisition of seismic data was in the winter, so the wind was blowing very hard.At the same time, the gangue field and air shafts were under construction.What is worse, the random noise interference was more serious in this area because there were many vehicles in adjacent industrial areas.Figure 10 shows the actual seismic profiles with random noise.It could be found that the existing random noise reduced the SNR of the seismic data and influenced the continuity of events in seismic records.
The ICEEMD and the VMDC were used for noise attenuation of post-stack seismic data.The signal reconstruction process was introduced by taking the 50th channel as an example.Firstly, ICEEMD and VMDC were used to decompose seismic signals.Then, by ICEEMD, the correlation coefficients between each IMF and the original signal were calculated as 0.1432, 0.5058 and 0.7329.By VMD, the correlation coefficients between each IMF and the original signal were calculated as 0.7221, 0.6611 and 0.1156.Therefore, the IMF2 and IMF3 decomposed by ICEEMD, and IMF1 and IMF2 decomposed by VMDC, participated in the signal constructions.The signals in other seismic traces were reconstructed by similar steps, as shown in Figures 11 and 12. Through the comparative analysis of Figure 10, Figure 11 and Figure 12, it is found that both ICEEMD and the VMDC could suppress the random noise to a certain extent, and the latter could significantly enhance the continuity of events.Above all, the VMDC method proposed in this paper has obvious noise attenuation effects and could make the events more clear and continuous.What is more, it could improve the SNR of seismic data and the smoothing of each channel of the seismic record.It also could better reflect shapes of strata.It is shown that the method could achieve good effects in random noise attenuation in realworld seismic data.

Conclusions
This paper proposed a new method for seismic random noise attenuation based on VMD and correlation coefficients, called VMDC.Under the situation of simple sine signals without noise, both ICEEMD and VMDC have better decomposition results, and VMDC has better calculation efficiency.After adding random noise to the simple sine signals and the synthetic seismogram, the testing results show that the VMDC has better noise attenuation effects.The application results in real-world seismic data indicate that the new method could significantly improve the SNR of seismic data, enhance the continuity of events as well as provide reliable basic data for further seismic data interpretation.Above all, the VMDC method proposed in this paper has obvious noise attenuation effects and could make the events more clear and continuous.What is more, it could improve the SNR of seismic data and the smoothing of each channel of the seismic record.It also could better reflect shapes of strata.It is shown that the method could achieve good effects in random noise attenuation in real-world seismic data.

Conclusions
This paper proposed a new method for seismic random noise attenuation based on VMD and correlation coefficients, called VMDC.Under the situation of simple sine signals without noise, both ICEEMD and VMDC have better decomposition results, and VMDC has better calculation efficiency.After adding random noise to the simple sine signals and the synthetic seismogram, the testing results show that the VMDC has better noise attenuation effects.The application results in real-world seismic data indicate that the new method could significantly improve the SNR of seismic data, enhance the continuity of events as well as provide reliable basic data for further seismic data interpretation.
, ICCEMD decomposes signals according to the frequency from high to low, while the VMD was the opposite.When the signal did not contain noise, the correlation coefficients between the decomposed signals by ICEEMD and VMD and the original signals could reach more than 0.97.In terms of calculation efficiency, ICEEMD took 21.96 s while VMD took only 2.15 s.The simple signal was expressed as ( ) 1( ) 2( ) z t x t x t   , where 1( ) sin(12 ) The number of sampling points was 1000.ICEEMD and VMD were used for signal decomposition.The results are shown in Figures 2 and 3. Developed by Colominas et al., ICEEMD has developed the EMD and could solve the problems of modal mixing in the application of EMD [20,28].As shown in Figures 2 and 3, ICCEMD decomposes signals according to the frequency from high to low, while the VMD was the opposite.When the signal did not contain noise, the correlation coefficients between the decomposed signals by ICEEMD and VMD and the original signals could reach more than 0.97.In terms of calculation efficiency, ICEEMD took 21.96 s while VMD took only 2.15 s.

Figure 2 .
Figure 2. Decomposition of the artificially mixed signal by improved complementary ensemble empirical mode decomposition (ICEEMD): (a) the mixed signal; (b) IMF1 (intrinsic mode function 1) of the signal; (c) IMF2 of the signal.

Figure 3 .
Figure 3. Decomposition of the artificially mixed signal by variational mode decomposition (VMD): (a) the mixed signal; (b) IMF1 of the signal; (c) IMF2 of the signal.

Figure 2 .
Figure 2. Decomposition of the artificially mixed signal by improved complementary ensemble empirical mode decomposition (ICEEMD): (a) the mixed signal; (b) IMF1 (intrinsic mode function 1) of the signal; (c) IMF2 of the signal.
The number of sampling points was 1000.ICEEMD and VMD were used for signal decomposition.The results are shown in Figures2 and 3. Developed by Colominas et al., ICEEMD has developed the EMD and could solve the problems of modal mixing in the application of EMD[20,28].As shown in Figures2 and 3, ICCEMD decomposes signals according to the frequency from high to low, while the VMD was the opposite.When the signal did not contain noise, the correlation coefficients between the decomposed signals by ICEEMD and VMD and the original signals could reach more than 0.97.In terms of calculation efficiency, ICEEMD took 21.96 s while VMD took only 2.15 s.

Figure 2 .
Figure 2. Decomposition of the artificially mixed signal by improved complementary ensemble empirical mode decomposition (ICEEMD): (a) the mixed signal; (b) IMF1 (intrinsic mode function 1) of the signal; (c) IMF2 of the signal.

Figure 3 .
Figure 3. Decomposition of the artificially mixed signal by variational mode decomposition (VMD): (a) the mixed signal; (b) IMF1 of the signal; (c) IMF2 of the signal.

Figure 3 .
Figure 3. Decomposition of the artificially mixed signal by variational mode decomposition (VMD): (a) the mixed signal; (b) IMF1 of the signal; (c) IMF2 of the signal.

Figure 4 .
Figure 4. Decomposition of the artificially mixed signal with 20% additive noise by ICEEMD: (a) signal; (b) IMF1 of the signal; (c) IMF2 of the signal; (d) IMF3 of the signal; (e) IMF4 of the signal; and (f) IMF5 of the signal.

Figure 4 .
Figure 4. Decomposition of the artificially mixed signal with 20% additive noise by ICEEMD: (a) the signal; (b) IMF1 of the signal; (c) IMF2 of the signal; (d) IMF3 of the signal; (e) IMF4 of the signal; and (f) IMF5 of the signal.

Figure 5 .
Figure 5. Decomposition of the artificially mixed signal with 20% additive noise by VMD: (a) the signal; (b) IMF1 of the signal; (c) IMF2 of the signal; (d) IMF3 of the signal.

Figure 6 .
Figure 6.Reconstruction results by VMDC (red) and the original signals without noise (blue).

Figure 5 .
Figure 5. Decomposition of the artificially mixed signal with 20% additive noise by VMD: (a) the signal; (b) IMF1 of the signal; (c) IMF2 of the signal; (d) IMF3 of the signal.

Figure 5 .
Figure 5. Decomposition of the artificially mixed signal with 20% additive noise by VMD: (a) the signal; (b) IMF1 of the signal; (c) IMF2 of the signal; (d) IMF3 of the signal.

Figure 6 .
Figure 6.Reconstruction results by VMDC (red) and the original signals without noise (blue).

Figure 6 .
Figure 6.Reconstruction results by VMDC (red) and the original signals without noise (blue).

Figure 8 .
Figure 8. Random noise attenuation results (left) and the noise (right) by ICEEMD.

Figure 9 .
Figure 9. Random noise attenuation results (left) and the noise (right) by VMDC.

Figure 8 .
Figure 8. Random noise attenuation results (left) and the noise (right) by ICEEMD.

Figure 9 .
Figure 9. Random noise attenuation results (left) and the noise (right) by VMDC.

Figure 8 .
Figure 8. Random noise attenuation results (left) and the noise (right) by ICEEMD.

Figure 8 .
Figure 8. Random noise attenuation results (left) and the noise (right) by ICEEMD.

Figure 9 .
Figure 9. Random noise attenuation results (left) and the noise (right) by VMDC.Figure 9. Random noise attenuation results (left) and the noise (right) by VMDC.

Figure 9 .
Figure 9. Random noise attenuation results (left) and the noise (right) by VMDC.Figure 9. Random noise attenuation results (left) and the noise (right) by VMDC.

Table 1 .
Person correlation coefficients and correlation intensity.
4 represented that there was weak correlation between the IMF component and the original signal.The IMF component contained the effective signal and random noise.It was decomposed into n