Next Article in Journal
Adaptive Ant Colony Optimization with Sub-Population and Fuzzy Logic for 3D Laser Scanning Path Planning
Next Article in Special Issue
Decoupling and Parameter Extraction Methods for Conical Micro-Motion Object Based on FMCW Lidar
Previous Article in Journal
Enhancing Stress Detection: A Comprehensive Approach through rPPG Analysis and Deep Learning Techniques
Previous Article in Special Issue
Design of Optical System for Ultra-Large Range Line-Sweep Spectral Confocal Displacement Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CEEMDAN-LWT De-Noising Method for Pipe-Jacking Inertial Guidance System Based on Fiber Optic Gyroscope

1
School of Engineering and Technology, China University of Geosciences, Beijing 100083, China
2
Technology Innovation Center for Directional Drilling Engineering, Ministry of Natural Resources, Langfang 065099, China
3
Innovation Base for Directional Drilling Engineering, Geological Society of China, Langfang 065099, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(4), 1097; https://doi.org/10.3390/s24041097
Submission received: 22 January 2024 / Revised: 5 February 2024 / Accepted: 6 February 2024 / Published: 7 February 2024
(This article belongs to the Special Issue Important Achievements in Optical Measurements in China 2022–2023)

Abstract

:
An inertial guidance system based on a fiber optic gyroscope (FOG) is an effective way to guide long-distance curved pipe jacking. However, environmental disturbances such as vibration, electromagnetism, and temperature will cause the FOG signal to generate significant random noise. The random noise will overwhelm the effective signal. Therefore, it is necessary to eliminate the random noise. This study proposes a hybrid de-noising method, namely complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN)—lifting wavelet transform (LWT). Firstly, the FOG signal is extracted using a sliding window and decomposed by CEEMDAN to obtain the intrinsic modal function (IMF) with N different scales and one residual. Subsequently, the effective IMF components are selected according to the correlation coefficient between the IMF components and the FOG signal. Due to the low resolution of the CEEMDAN method for high-frequency components, the selected high-frequency IMF components are decomposed with lifting wavelet transform to increase the resolution of the signal. The detailed signals of the LWT decomposition are de-noised using the soft threshold de-noising method, and then the signal is reconstructed. Finally, pipe-jacking dynamic and environmental interference experiments were conducted to verify the effectiveness of the CEEMDAN-LWT de-noising method. The de-noising effect of the proposed method was evaluated by SNR, RMSE, and Deviation and compared with the CEEMDAN and LWT de-noising methods. The results show that the CEEMDAN-LWT de-noising method has the best de-noising effect with good adaptivity and high accuracy. The navigation results of the pipe-jacking attitude before and after de-noising were compared and analyzed in the environmental interference experiment. The results show that the absolute error of the pipe-jacking pitch, roll, and heading angles is reduced by 39.86%, 59.45%, and 14.29% after de-noising. The maximum relative error of the pitch angle is improved from −0.74% to −0.44%, the roll angle is improved from 2.07% to 0.79%, and the heading angle is improved from −0.07% to −0.06%. Therefore, the CEEMDAN-LWT method can effectively suppress the random errors of the FOG signal caused by the environment and improve the measurement accuracy of the pipe-jacking attitude.

1. Introduction

Long-distance curved pipe-jacking has significant engineering and economic benefits in tunneling projects crossing densely built-up urban areas, heavily trafficked road sections, and oversized cross-sections [1,2]. However, the traditional pipe-jacking guidance method relies on the visibility conditions of the environment and does not apply to long-distance curved pipe-jacking. The fiber optic gyroscope (FOG) is an all-solid-state angular velocity-sensitive device, and it can be sensitive to the angular velocity of the Earth’s rotation, which has the advantage of autonomous navigation [3]. Therefore, the pipe-jacking inertial guidance system uses FOG as a primary sensor to provide real-time attitude measurements for long-distance curved pipe-jacking [4]. But in practical applications, large motor and transformer operations will produce vibrations, electromagnetic disturbances, and other environmental interference due to the rotation of the pipe-jacking cutter. The environment will affect components such as the fiber optic ring in the FOG, and the FOG signal will contain a large amount of random noise, resulting in insufficient long-term accuracy [5]. Therefore, it is necessary to reduce noise to obtain more accurate FOG signals and improve the accuracy of the pipe-jacking attitude [6].
Commonly used FOG random noise filtering methods include digital low-pass filters and time series forecasting methods [7,8]. The time series forecasting method is based on Kalman filtering by establishing an autoregressive (AR) or autoregressive moving average (ARMA) model for the FOG drifting signal, and optimal estimation is performed by a strong tracking Kalman filtering method [9,10]. The Kalman filter cannot strictly distinguish between the useful signal and interference noise in the high frequency part. Therefore, if the low-pass filter is too narrow, it will cause a loss of useful signal, and if the low-pass filter is too wide, it will lead to a poor filtering effect. However, the Kalman filtering method based on AR and ARMA models is proposed for smooth signals. Non-smooth FOG signals need to be smoothed before using these filtering methods.
The empirical mode decomposition (EMD) proposed by Huang E. et al. has achieved effective results in non-smooth and non-linear signal processing [11]. The EMD can adaptively decompose a signal into a series of intrinsic mode functions (IMFs) according to its own characteristics. The IMF is a description of the signal in different scales. Compared with wavelet transform, EMD has desirable adaptivity and decomposition effects in signal processing for non-linear and non-smooth signals. However, the EMD method has an issue that the IMFs may interact or alias with each other during the decomposition process, i.e., the signals of different scales and frequencies appear in the same IMF component, or signals of the same scale and frequency are decomposed into multiple IMF components [12]. Wu et al. [13] proposed an ensemble empirical mode decomposition (EEMD) method to solve the modal aliasing problem. The EEMD method is to decompose the original signal by adding white noise several times, and then averaging the results of the multiple decompositions to obtain the final IMF. Since the EEMD method adds white noise to the original signal several times, it will cause reconstruction errors. Yeh et al. [14] proposed a complementary ensemble empirical mode decomposition (CEEMD) method based on EEMD. The CEEMD works by adding two opposite white noise signals to the original signal for EMD decomposition. The residual white noise in the reconstructed signal is effectively eliminated under the premise of guaranteeing the decomposition effect is comparable to EEMD. However, not only do both EEMD and CEEMD have a problem with computation and storage complexity, but also the decomposition is too dependent on the amplitude of the added white noise and the number of accumulated averages. Thus, Torres et al. [15] proposed complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN), which obtains the IMF by calculating the residual signal by adding adaptive white noise at each stage of the EMD decomposition. The CEEMDAN method has complete decomposition, and the reconstruction error is very small each time, which solves the problem of the low computational efficiency of EEMD and CEEMD, and has a strong adaptability to non-smooth and non-linear signals [16]. Therefore, for the non-linear and non-smooth signals of the FOG in the pipe-jacking inertial guidance system, the CEEMDAN method is suitable for the noise reduction. However, the traditional CEEMDAN method directly removes the high-frequency IMF components for noise reduction, which can easily cause the loss of useful high-frequency signals, and random noise and spurious IMF components may also be contained in the mid-frequency and low-frequency IMF components.
Zhao et al. [17] propose a dynamic error compensation method for FOG based on CEEMDAN that properly separates the effective signal and dynamic error. There is no further algorithm that is added to filter out the noise signal. Wang et al. [18] propose a method for constructing an FOG temperature drift compensation model based on CEEMDAN, and use an adaptive Kalman filter (AKF) to filter mixed noise, which effectively reduces temperature errors. However, the CEEMDAN has low resolution for high-frequency signals, and the Kalman filter has no ability to improve the resolution of high-frequency signals. Therefore, a CEEMDAN-LWT denoising method combining CEEMDAN and lifting wavelet transform (LWT) is proposed in this study to preprocess the noisy signals of the FOG. The useful IMF components are screened by calculating the correlation coefficient between the IMF components and the original signal. Considering the low resolution of the CEEMDAN method for high-frequency components, the high-frequency IMF components are screened for the lifting wavelet transform to increase the resolution [19]. Then, wavelet soft threshold denoising is performed, and finally, the screened IMF components are reconstructed.
The proposed de-noising method has the complete features of CEEMDAN decomposition. At the same time, the method has the fast de-noising speed and high-resolution features of LWT decomposition. Furthermore, the method achieves the real-time processing of the FOG signal by using a sliding pane. By comparing the results of the denoising effects and Allan variance, the proposed method effectively reduces a majority of the random noise in the FOG signal within the pipe-jacking inertial guidance system. In addition, the proposed method also improves the accuracy of the pipe-jacking attitude results.

2. Denoising Principle of CEEMDAN-LWT

2.1. CEEMDAN

The CEEMDAN decomposition is achieved in the following steps.
Adding a specific white noise N ( t ) to the original signal s ( t ) as the signal to be decomposed, the i-th signal can be expressed as follows.
s i ( t ) = s ( t ) + N i ( t )
N i ( t ) = σ k ( t ) E i ( ε ( t ) )
In Equations (1) and (2), i represents the number of times Gaussian white noise is added, i = 1 , 2 , n ; σ k ( t ) is a parameter for adaptive adjustment, k = 1 , 2 , m ; ε ( t ) is the Gaussian white noise; E i ( ε ( t ) ) is the i th EMD integration average white noise.
First, the first-order IMF component is obtained by EMD decomposition on s i ( t ) , i.e., IMF1.
I M F 1 = s ( t ) 1 n i = 1 n M ( s i ( t ) )
where M(·) represents the local average of the envelope satisfying the IMF filtering threshold. The obtained decomposition residue is r 1 ( t ) = s ( t ) I M F 1 .
Next, the decomposition of the residual r 1 ( t ) is continued by adding white noise N ( t ) and the second-order IMF component is obtained.
I M F 2 = r 1 ( t ) 1 n i = 1 n M ( r 1 ( t ) + N i ( t ) )
The above steps are repeated until the residuals cannot be discretized, and the k-th order IMF components are as follows.
I M F k = r k ( t ) 1 n i = 1 n M ( r k ( t ) + N i ( t ) ) = r k 1 ( t ) 1 n i = 1 n M ( r k 1 ( t ) + σ k 1 ( t ) E i ( ε ( t ) ) )
The n is the number of white noises added, which determines the accuracy of the IMF and is usually set between 100 and 1000. Another important factor affecting the decomposition results is σ k 1 ( t ) . The value of σ k 1 ( t ) is usually set between 0.1 and 0.3.

2.2. Correlation Analysis of IMF Components

Two signals, x(i) and y(i), are set and the correlation coefficient R is an indicator of the degree of correlation between the two signals.
R = x ( i ) E ( x ) y ( i ) E ( y ) x ( i ) E ( x ) 2 y ( i ) E ( y ) 2
When R is close to 1, the two signals have a strong positive correlation; when R is close to −1, the two signals have a strong negative correlation; and when R is closer to 0, the two signals have a nearly weak correlation. After CEEMDAN decomposition, the effective IMF components are screened, and pseudo IMF components are removed by calculating the correlation coefficient between two signals. This process aims to enhance the accuracy of the signal [20,21].

2.3. LWT

The screened IMF components are de-noised by lifting wavelet transform, and the lifting algorithm is divided into split, prediction, and updating [22,23].
Step 1: Split. The input signal s ( k ) , k Z is divided into two wavelet subsets by parity.
s o ( k ) = s ( 2 k + 1 ) , k Z s e ( k ) = s ( 2 k ) , k Z
where s o ( k ) is the sequence of odd samples, s e ( k ) is the sequence of even samples, and k is the number of signal samples.
Step 2: Prediction. Set P represents the prediction operator, using s e ( k ) to predict the dual sequence s o ( k ) , and the predicted signal s o ( k ) is represented as P [ s e ( k ) ] . The prediction error of the even sample sequence is defined as the detail signal d ( k ) , which is the high-frequency component obtained after transformation.
d ( k ) = s o ( k ) P [ s e ( k ) ]
Step 3: Update. Set U represents the update operator that updates d ( k ) based on the detail signal s e ( k ) to approximate the overall characteristics of the original signal more closely while maintaining the scale properties of the original signal.
c ( k ) = s e ( k ) U [ d ( k ) ]
where c ( k ) is the approximation signal and is the low-frequency component obtained from the transform.

2.4. Soft Threshold De-Noising and Signal Reconstruction

Soft threshold de-noising is performed on the high-frequency components obtained from the LWT algorithm.
W ^ = sgn ( W ) ( W λ ) 0 W > λ W λ
where W represents the initial wavelet coefficients corresponding to the high-frequency component, W ^ denotes the wavelet coefficients of the high-frequency component after threshold processing, and λ represents the threshold.
λ = σ ln N
where σ is the standardized variance of the noise and N is the signal length.
After soft threshold de-noising, the de-noised signal is obtained by reconstructing each high-frequency component d ( k ) and low-frequency component c ( k ) . The flowchart of the CEEMDAN-LWT algorithm is shown in Figure 1.

3. Simulation Signal Analysis

Matlab is used to simulate the FOG signal and verify the effectiveness of the CEEMDAN-LWT de-noising algorithm. The simulated signal is a sinusoidal signal with frequency of 0.1 Hz, amplitude of 1 °/h, sampling frequency of 100 Hz, and sampling time of 50 s. Then, the white noise is added; the mean value of the white noise is 0, and the variance of the white noise is 1 °/h. The noisy signal as shown in Figure 2.
The noisy signal is subsequently conducted to CEEMDAN decomposition, resulting in 11 IMF components and a residue, shown in Figure 3.
In Figure 3, the frequencies of the IMF components gradually decrease from top to bottom, and the last item, Res, is the remaining residuals. It shows that the noise mainly exists in the IMF components at high frequencies, the IMF components at low frequencies contain most of the effective signal, and IMF8 and IMF9 are nearly close to the simulated signal. Then, the correlation between each of the obtained IMF components with the noise signal is computed, and the results of the correlation analysis are shown in Table 1.
According to the results of correlation analysis, the high-frequency IMF components: IMF1, IMF2, IMF3, and IMF4 with a correlation greater than 0.2 are extracted to conduct LWT soft threshold de-noising. The IMF components with correlation coefficients less than 0.1 are removed and the other IMF components are retained. The signal is reconstructed with the de-noised high-frequency IMF components and the remaining IMF components. The obtained de-noised signal is compared with the noise signal, and the simulated signal. The compared results are shown in Figure 4.
In Figure 4, the blue curve is the noise-containing signal, the red curve is the signal obtained after de-noising using the CEEMDAN-LWT method, and the yellow curve is the signal without noise in the simulation. The noise of the red curve is significantly reduced compared with the blue curve, and the trend of the red curve is consistent with the yellow curve. The red curve fluctuates around the yellow curve and without signal distortion proves that the algorithm proposed in this study can filter out the noise in the signals well.
The frequency analysis is conducted on the signals before and after de-noising in Figure 5.
In Figure 5, the frequency of the signal after noise reduction is mainly distributed at 0~0.2 Hz, and the signal frequency set by the simulation signal is 0.1 Hz, so the effective information of the signal is better retained. The signal is better decomposed to different frequencies by the CEEMDAN-LWT and is more effective in removing the high-frequency noise. The signal after noise reduction is close to the real signal, and the features of the effective signal are retained to the maximum extent.
This study applies the signal-to-noise ratio (SNR), standard deviation (SD), deviation (D) and root mean square error (RMSE) to identify the effectiveness of the signal denoising. The equations for them are shown in below [24].
S N R = 10 log 10 i = 1 N s 2 ( i ) i = 1 N s ( i ) s t ( i ) 2
R M S E = 1 N i = 1 N s t ( i ) s t ( i ) 2
D = 1 N i = 1 N s t ( i ) s ^ t ( i )
S D = 1 N i = 1 N s t ( i ) s ^ t ( i ) 2
where s ( i ) is the simulated signal, s t ( i ) is the noise-containing signal, s t ( i ) is the signal after noise reduction, s ^ t ( i ) is the average value of the signal after noise reduction, and N is the number of sampling points. The larger the SNR, the smaller the RMSE, SD and deviation, the better the de-noising effect. The SNR, RMSE, and deviation represent the de-noising effect of the simulated noise-containing signal. The de-noising effect of the signal is shown in Table 2.
Table 2 compares CEEMDAN, LWT, and CEEMDAN-LWT methods with the de-noising effects. The CEEMDAN method removes the first two high-frequency IMF components directly, and the LWT method conducts a 5-level decomposition and soft thresholds de-noising with the signal. According to SNR, SD, RESM, and D results, CEEMDAN-LWT has the best evaluation effect, proving that the algorithm proposed in this study has a better de-noising effect.

4. FOG Signal Analysis

4.1. Dynamic Pipe-Jacking Experiment

To verify the de-noising effect of the CEEMDAN-LWT method in the pipe-jacking working environment, a PA-GS300 FOG was installed on the running pipe-jacking machine to sample the FOG signal, as shown in Figure 6.
The sampling frequency was 100 Hz. In the dynamic experiment period, the pipe-jacking attitude is adjusted by the correction cylinders. The adjustment steps include setting the pipe-jacking head up first, then adjusting the head to turn from left to right. Finally, the pipe-jacking head is turned back and put down, as shown in Figure 7.
Sliding window data for the CEEMDAN-LWT de-noising method is chosen. The x-axis gyroscope data is taken as an example, as shown in Figure 8.
The original gyro signal is decomposed using the CEEMDAN method and 12 IMF components and a residue are obtained, as shown in Figure 9.
The correlation analysis of each IMF component with the FOG signal is shown in Table 3.
The high-frequency IMF1, IMF2, and IMF3 whose correlation coefficients greater than 0.2 are screened and processed with LWT and soft threshold de-noising. The level of decomposition is 5. The IMF components with correlation coefficients of less than 0.1 are removed, and the signal is reconstructed. The de-noised signals are obtained, as shown in Figure 10.
Figure 10 shows that the de-noised signal better retains the characteristics of the original signal and filters out most of the noise, which proves that the de-noising method proposed in this paper is suitable for the FOG signals. The frequency analysis of the noise-containing signal and the de-noised signal are shown in Figure 11.
Figure 11a shows that the frequency of the FOG signal is distributed from 0–50 Hz, and the main frequency is distributed near 0 Hz, 21 Hz, 25 Hz, and 37 Hz. During the dynamic experiment, the attitude adjustment is slow, and the effective signal is mainly distributed near 0 Hz. Figure 11b shows that the frequency of the de-noised signal is mainly distributed between 0–1 Hz, and the algorithm filters out most of the high-frequency noise signals, which proves the effectiveness of the algorithm proposed in this paper. The de-noising effect is analyzed and compared in Table 4.
Table 4 compares the de-noising effect of the FOG signal processed using the CEEMDAN, LWT, and CEEMDAN-LWT methods. The CEEMDAN method removes the first two high-frequency IMF components of the FOG signal directly, and the LWT method performs a 5-level decomposition and soft thresholding of the FOG signal. According to the SNR and SD results, CEEMDAN-LWT shows the best de-noising effect. It proves that the algorithm proposed in this study effectively improves the de-noising effect. The de-noised signal is solved with inertial navigation to obtain the pipe-jacking attitude, and the results are as follows.
In Figure 12, the blue curve is the result of the inertial navigation solving of the FOG signal, and the red curve is the result of the inertial navigation solving of the FOG signal after de-noising using the CEEMDAN-LWT algorithm. The pitch angle of the FOG signal changes from −2.4° to about −0.5°and then back to −2.4°, which is consistent with the trend of the attitude adjustment, proving that the de-noised signal is not distorted. The detailed graph shows that the fluctuation in the pitch angle solved using the CEEMDAN-LWT algorithm is minor and smoother. It proves that the algorithm proposed in this paper can filter out random noise better.
In Figure 13, the roll angle is set with no change, and the roll angle after the CEEMDAN-LWT algorithm fluctuates around 0.05°. The trend of the roll angle change before and after de-noising is consistent, which proves that there is no distortion of the FOG signal after de-noising. In the detailed graph, the fluctuation of the red curve is reduced compared with that of the blue curve, which proves that the algorithm proposed in this study is effective.
In Figure 14, the heading angle changes from −0.2° to −0.7° and then to 0.7° and finally to near 0°, which is consistent with the trend of the attitude adjustment from the left to the right, which proves that there is no distortion after de-noising. In the detailed graph, the heading angle after de-noising is smoother compared with the fluctuation before de-noising, and the algorithm effectively filters out part of the random noise. The attitude de-noising effect of the pipe-jacking dynamic experiment is evaluated as follows.
Table 5 compares the SD of the attitude before and after de-noising; the SD of the attitude after de-noising is reduced, which proves that the proposed algorithm improves the accuracy of attitude results. The SNR results of the de-noised attitude are improved, which proves that the proposed algorithm effectively filters out the random noise.

4.2. Pipe-Jacking Environmental Interference Experiment

A pipe-jacking environmental interference experiment was conducted to verify the adaptive de-noising effect of the CEEMDAN-LWT algorithm on the FOG during pipe-jacking operations. The sampling frequency was 100 Hz. The 10,000 sample points were intercepted and processed by the algorithm. The frequency converter and the motor were turned on to rotate the pipe jacking cutter at 4650 sampling points. The three-axis data of the FOG after the de-noising by the proposed algorithm are shown in Figure 15.
The blue curve in Figure 15 is the original signal measured by the FOG, and the noise is enhanced during the pipe-jacking operation at 4650 sampling points. The red curve is the signal after the de-noising algorithm, and the signal amplitude is consistent after de-noising in different environments, proving that the algorithm proposed in this paper has a certain degree of adaptivity. Meanwhile, the proposed algorithm filters out most of the random noise. The attitude results are compared before and after de-noising in Figure 16.
The blue curve in Figure 16 is the attitude result before de-noising. The pipe-jacking operation starts at 4650 sampling points, and the fluctuation of the attitude results is obviously increased. The red curve is the result of attitude solving using the algorithm proposed in this study, and the noise of attitude is obviously reduced, which proves that the de-noising algorithm proposed in this study can effectively filter out the environmental noise. Table 6 is a statistical analysis of the attitude results.
In Table 6, the initial attitude is the result of the attitude at the 4650th sampling point. The ultimate attitude is the result of the attitude at the 10,000th sampling point. The average attitude is the average value of between 4650 and 10,000 sampling points. The pipe-jacking attitude remains static throughout the environmental interference experiment. In theory, the pipe-jacking attitude should not change during the experiment, and the initial attitude value represents the true value. However, due to random noise and gyro drift error, the attitude will drift. By comparing the initial and ultimate attitudes before and after denoising, it is evident that the proposed method reduces a part of the drift error caused by random noise. Table 6 also shows that the average attitude of the CEEMDAN-LWT method is closer to the initial attitude, and the absolute error of the pitch angle is reduced by 39.86% relative to the original signal, and the roll and heading angle are reduced by 59.45% and 14.29%. The maximum relative error of the pitch angle is improved from −0.74% to −0.44%, the roll angle is improved from 2.07% to 0.79%, and the heading angle is improved from −0.07% to −0.06%. It is proven that the CEEMDAN-LWT algorithm proposed in this study filters out the random noise effectively and improves the accuracy of the attitude results. The de-noising effect is evaluated and compared in Table 7.
In Table 7, the CEEMDAN-LWT algorithm has the best de-noising effect compared to the CEEMDAN and LWT algorithms, but the heading angle de-noising effect is weak. Finally, the performance of gyro data before and after de-noising was analyzed using Allan variance as shown in Figure 17.
The blue curve in Figure 17 is the Allan variance curve of the original signal, and the red curve is the Allan variance curve of the CEEMDAN-LWT de-noised signal. Figure 17 shows that the red curves are all in the lower left of the blue curves, which indicates that the noise of the FOG signal de-noised using CEEMDAN-LWT is reduced. Therefore, the proposed de-noising method is effective. The specific error indicators are shown in Table 8.
In Table 8, Q is quantization noise, B is bias instability, N is random walk, and K is the rate of the random walk. The specific noise indicators of the CEEMDAN-LWT denoised signal become smaller, where the value of Q is very small and almost negligible. It is proven that the CEEMDAN-LWT method proposed in this study reduces the noise of the FOG signal effectively.

5. Conclusions

According to the characteristics of the FOG signals in the pipe-jacking inertial guidance system, this study proposes a hybrid de-noising method based on CEEMDAN-LWT. The FOG signal is decomposed by CEEMDAN, the main IMF components are screened by correlation coefficients, the high-frequency IMF components in the screening are de-noised using LWT and soft threshold de-noising, and finally, the signal is reconstructed. Simulation, dynamic, and environmental interference experiments were conducted, which compared the analysis of the frequency and the evaluation of the de-noising effect using the CEEMDAN and LWT de-noising methods. The results show that the method proposed in this study can effectively reduce the influence of random noise on the signal and accurately reflect the changing characteristics of the signal. The de-noised pipe-jacking attitude in the environmental interference experiment was closer to the theoretical value, the absolute error of the pitch angle was reduced by 39.86%, the roll angle was reduced by 59.45%, and the heading angle was reduced by 14.29%. The maximum relative error of the pitch angle was improved from −0.74% to −0.44%, the roll angle was improved from 2.07% to 0.79%, and the heading angle was improved from −0.07% to −0.06%. In conclusion, the CEEMDAN-LWT de-noising method proposed in this study combines the CEEMDAN and LWT decomposition features, which have the advantages of complete decomposition, fast speed, and high-resolution. At the same time, the method achieves real-time processing by using a sliding window. The study results show that the proposed method could filter out most of the random noise of the FOG signal in the pipe-jacking inertial guidance system and improve the accuracy of the pipe-jacking attitude results, which has better engineering application value.

Author Contributions

Conceptualization, L.W.; Formal analysis, Y.Z.; Resources, G.Y.; Data curation, Y.Z.; Writing—original draft, Y.Z.; Writing—review & editing, L.W.; Funding acquisition, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Technology Innovation Center for Directional Drilling Engineering, Ministry of Natural Resources, Innovation Base for Directional Driling Engineering Geoloeical Society of China (KF202301) and by National Natural Science Foundation of China (U22A20184 and 41902320).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ersoy, N.; Yavuz, E. Geodetic works in micro tunneling projects excavated by pipe-jacking method. In Proceedings of the SGEM International Multidisciplinary Scientific Geo-Conference and Expo, Albena, Bulgaria, 18–24 June 2011. [Google Scholar]
  2. Pan, G.R.; Fan, W. Automatic guidance system for long-distance curved pipe-jacking. KSCE J. Civ. Eng. 2020, 24, 2505–2518. [Google Scholar] [CrossRef]
  3. Liu, W.; Liu, J.; Shen, Q. Integrated modeling and filtering of fiber optic gyroscope’s random errors. Opto-Electron. Eng. 2018, 45, 180082. [Google Scholar]
  4. Xu, H.; Wang, L.; Zu, Y.; Gou, W.; Hu, Y. Application and Development of Fiber Optic Gyroscope Inertial Navigation System in Underground Space. Sensors 2023, 23, 5627. [Google Scholar] [CrossRef] [PubMed]
  5. Cao, Y.; Xu, W.; Lin, B. A method for temperature error compensation in fiber-optic gyroscope based on machine learning. Optik 2022, 256, 168765. [Google Scholar] [CrossRef]
  6. Honório, B.; Drummond, R.; Vidal, A. Well log denoising and geological enhancement based on discrete wavelet transform and hybrid thresholding. Energy Explor. Exploit. 2012, 30, 417–433. [Google Scholar] [CrossRef]
  7. Dang, S.; Tian, W.; Qian, F. EMD- and LWT-based stochastic noise eliminating method for fiber optic gyro. Measurement 2011, 44, 2190–2193. [Google Scholar] [CrossRef]
  8. Li, Z.; Meng, Z.; Liu, T. A novel method for determining and improving the quality of a quadrupolar fiber gyro coil under temperature variations. Opt. Express 2013, 21, 2521–2530. [Google Scholar] [CrossRef] [PubMed]
  9. Lv, P.; Liu, J.; Lai, J. Research on threshold algorithm in real-time wavelet de-noising of FOG. J. Proj. Rocket. Missiles Guid. 2009, 29, 18–22. [Google Scholar]
  10. Wang, P.; Li, G.; Gao, Y. A compensation method for gyroscope random drift based on unscented Kalman filter and support vector regression optimized by adaptive beetle antennae search algorithm. Appl. Intell. 2023, 53, 4350–4365. [Google Scholar] [CrossRef]
  11. Huang, E.; Shen, Z.; Long, S. The empirical mode de-composition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  12. Tang, B.; Dong, S.; Song, T. Method for eliminating mode mixing of empirical mode decomposition based on the revised blind source separation. Signal Process. 2012, 92, 248–258. [Google Scholar] [CrossRef]
  13. Wu, Z.; Huang, N. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009, 1, 1–41. [Google Scholar] [CrossRef]
  14. Yeh, J.R.; Shieh, J.R.; Huang, N. Complementary ensemble empirical mode decomposition: A novel noise enhanced data analysis method. Adv. Adapt. Data Anal. 2010, 2, 135–156. [Google Scholar] [CrossRef]
  15. Torres, M.; Colominas, M.; Schlotthauer, G. A complete ensemble empirical mode decomposition with adaptive noise. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing, Prague, Czech Republic, 22–27 May 2011; pp. 4144–4147. [Google Scholar]
  16. Belazi, A.; El-Latif, A.; Diaconu, A. Chaos-based partial image encryption scheme based on linear fractional and lifting wavelet transforms. Opt. Lasers Eng. 2017, 88, 37–50. [Google Scholar] [CrossRef]
  17. Zhao, S.; Zhou, Y.; Shu, X. Analysis of fiber optic gyroscope dynamic error based on CEEMDAN. Opt. Fiber Technol. 2022, 69, 102835. [Google Scholar] [CrossRef]
  18. Wang, X.; Cui, Y.; Cao, H. Temperature Drift Compensation of Fiber Optic Gyroscopes Based on an Improved Method. Micromachines 2023, 14, 1712. [Google Scholar] [CrossRef]
  19. Li, M.; Mao, H.; Zhang, Y. Fusion algorithm for multi–sensor images based on PCA and lifting wavelet transformation. N. Z. J. Agric. Res. 2007, 50, 667–671. [Google Scholar]
  20. Komaty, A.; Boudraa, A.; Augier, B. EMD-based filtering using similarity measure between probability density functions of IMFs. IEEE Trans. Instrum. Meas. 2014, 63, 27–34. [Google Scholar] [CrossRef]
  21. Boudraa, A.-O.; Cexus, J.-C. EMD-based signal filtering. IEEE Trans. Instrum. Meas. 2007, 56, 2196–2202. [Google Scholar] [CrossRef]
  22. Ding, W.; Wu, X.; Li, S. Adaptive directional lifting based wavelet transform for image coding. IEEE Trans. Image Process 2007, 16, 416–427. [Google Scholar] [CrossRef]
  23. Ding, W.; Wu, F.; Li, S. Lifting-based wavelet transform with directionally spatial prediction. In Proceedings of the Process Picture Coding Symposium, San Francisco, CA, USA, 15 December 2004; pp. 483–488. [Google Scholar]
  24. Zhang, J.; Guo, Y.; Shen, Y.; Zhao, D.; Li, M. Improved CEEMDAN–wavelet transform de-noising method and its application in well logging noise reduction. J. Geophys. Eng. 2018, 15, 775–787. [Google Scholar] [CrossRef]
Figure 1. The flowchart of the CEEMDAN-LWT de-noising algorithm.
Figure 1. The flowchart of the CEEMDAN-LWT de-noising algorithm.
Sensors 24 01097 g001
Figure 2. Schematic of the simulated noise signal.
Figure 2. Schematic of the simulated noise signal.
Sensors 24 01097 g002
Figure 3. The IMF components and a residue by CEEMDAN decomposition of the simulated signal.
Figure 3. The IMF components and a residue by CEEMDAN decomposition of the simulated signal.
Sensors 24 01097 g003
Figure 4. De-noising results of the simulated signal based on the CEEMDAN-LWT method.
Figure 4. De-noising results of the simulated signal based on the CEEMDAN-LWT method.
Sensors 24 01097 g004
Figure 5. Frequency analysis of noise-containing and de-noised signals: (a) noise-containing signal (b) denoised signal.
Figure 5. Frequency analysis of noise-containing and de-noised signals: (a) noise-containing signal (b) denoised signal.
Sensors 24 01097 g005
Figure 6. The FOG in pipe-jacking inertial guidance system.
Figure 6. The FOG in pipe-jacking inertial guidance system.
Sensors 24 01097 g006
Figure 7. Schematic diagram of pipe-jacking attitude adjustment.
Figure 7. Schematic diagram of pipe-jacking attitude adjustment.
Sensors 24 01097 g007
Figure 8. The x-axis gyroscope signal in the dynamic pipe-jacking experiment.
Figure 8. The x-axis gyroscope signal in the dynamic pipe-jacking experiment.
Sensors 24 01097 g008
Figure 9. The CEEMDAN decomposed results of the FOG signals for the dynamic pipe-jacking experiment.
Figure 9. The CEEMDAN decomposed results of the FOG signals for the dynamic pipe-jacking experiment.
Sensors 24 01097 g009
Figure 10. Denoising results of the FOG signal based on the CEEMDAN-LWT method in the dynamic pipe-jacking experiment.
Figure 10. Denoising results of the FOG signal based on the CEEMDAN-LWT method in the dynamic pipe-jacking experiment.
Sensors 24 01097 g010
Figure 11. Frequency analysis of the noise-containing and de-noised signals: (a) noise-containing signal (b) de-noised signal.
Figure 11. Frequency analysis of the noise-containing and de-noised signals: (a) noise-containing signal (b) de-noised signal.
Sensors 24 01097 g011
Figure 12. Pitch angle solving results of the dynamic pipe-jacking experiment.
Figure 12. Pitch angle solving results of the dynamic pipe-jacking experiment.
Sensors 24 01097 g012
Figure 13. Roll angle solving results of the dynamic pipe-jacking experiment.
Figure 13. Roll angle solving results of the dynamic pipe-jacking experiment.
Sensors 24 01097 g013
Figure 14. Heading angle solving results of the dynamic pipe-jacking experiment.
Figure 14. Heading angle solving results of the dynamic pipe-jacking experiment.
Sensors 24 01097 g014
Figure 15. Denoising results of the FOG signal based on the CEEMDAN-LWT method in pipe-jacking environmental interference experiment.
Figure 15. Denoising results of the FOG signal based on the CEEMDAN-LWT method in pipe-jacking environmental interference experiment.
Sensors 24 01097 g015
Figure 16. Attitude solving results of the pipe-jacking environmental interference experiment.
Figure 16. Attitude solving results of the pipe-jacking environmental interference experiment.
Sensors 24 01097 g016
Figure 17. Allan variance curve comparison of the pipe-jacking environmental interference experiment: (a) The Allan variance of Gryo X, (b) the Allan variance of Gryo Y, (c) the Allan variance of Gryo Z.
Figure 17. Allan variance curve comparison of the pipe-jacking environmental interference experiment: (a) The Allan variance of Gryo X, (b) the Allan variance of Gryo Y, (c) the Allan variance of Gryo Z.
Sensors 24 01097 g017
Table 1. Results of correlation analysis.
Table 1. Results of correlation analysis.
IMF NumberCorrelationIMF NumberCorrelation
IMF10.64IMF70.13
IMF20.47IMF80.22
IMF30.35IMF90.57
IMF40.27IMF100.54
IMF50.19IMF110.10
IMF60.14Res0.06
Table 2. Compared results of denoising effect evaluation.
Table 2. Compared results of denoising effect evaluation.
MethodSNR (dB)RMSE (deg/h)D (deg/h)SD (deg/h)
CEEMDAN8.520.2870.0120.85
LWT7.680.3110.0160.87
CEEDAN-LWT10.550.2100.0060.75
Table 3. Correlation analysis results of each IMF component.
Table 3. Correlation analysis results of each IMF component.
IMF NumberCorrelationIMF NumberCorrelation
IMF10.53IMF80.24
IMF20.43IMF90.26
IMF30.26IMF100.19
IMF40.18IMF110.50
IMF50.13IMF120.61
IMF60.15Res0.36
IMF70.21
Table 4. Comparison of denoising effect in the dynamic pipe-jacking experiment.
Table 4. Comparison of denoising effect in the dynamic pipe-jacking experiment.
MethodSNR (dB)SD (deg/s)
CEEMDAN13.011.28 × 10−6
LWT11.941.30 × 10−6
CEEDAN-LWT18.181.19 × 10−6
Table 5. Evaluation of the attitude denoising effect in the dynamic pipe-jacking experiment.
Table 5. Evaluation of the attitude denoising effect in the dynamic pipe-jacking experiment.
AttitudeSD (deg/s)SNR (dB)
Original signalPitch0.4469-
Roll0.0077-
Heading0.0036-
Denoised signalPitch0.446128.3121
Roll0.007633.5005
Heading0.003510.0399
Table 6. Statistics of attitude results in the pipe-jacking environmental interference experiment.
Table 6. Statistics of attitude results in the pipe-jacking environmental interference experiment.
SignalAttitude AngleInitial Attitude/degUltimate Attitude/degAverage Attitude/degAbsolute
Error/deg
Relative Error/%
Original signalPitch−0.9873−0.9896−0.9880−0.0070~0.0073−0.74~0.71
Roll1.79921.80571.8024−0.0140~0.0373−0.78~2.07
Heading−20.5249−20.5111−20.51880~0.0138−0.07~0
CEEMDAN-LWT denoised signalPitch−0.9873−0.9874−0.9875−0.0043~0.0043−0.44~0.44
Roll1.79921.80051.8021−0.0068~0.0140−0.38~0.79
Heading−20.5249−20.5118−20.51860~0.0131−0.06~0
Table 7. Evaluation of the attitude denoising effect in the pipe-jacking environmental interference experiment.
Table 7. Evaluation of the attitude denoising effect in the pipe-jacking environmental interference experiment.
MethodAttitudeSNR (dB)SD (deg)D (deg)
CEEMDANPitch13.45200.0010−0.0008
Roll13.14350.00230.0024
Heading21.65130.00440.0063
LWTPitch13.36760.0011−0.0009
Roll12.67150.00250.0029
Heading21.75310.00450.0065
CEEDAN-LWTPitch13.59260.0009−0.0002
Roll13.22950.00180.0012
Heading22.34700.00420.0063
Table 8. Comparison of Allan variance between the original signal and CEEMDAN-LWT denoised signal.
Table 8. Comparison of Allan variance between the original signal and CEEMDAN-LWT denoised signal.
Original SignalCEEMDAN-LWT Denoised Signal
Gyro axialGyro xGyro yGyro zGyro xGyro yGyro z
Q/(″)0.0030.0310.003---
N/deg/h1/23.436 × 10−37.809 × 10−38.740 × 10−41.118 × 10−32.703 × 10−32.760 × 10−4
B/deg/h1.762 × 10−27.964 × 10−47.382 × 10−41.694 × 10−25.446 × 10−46.182 × 10−4
K/deg/h3/28.359 × 10−29.730 × 10−24.484 × 10−27.461 × 10−28.655 × 10−24.441 × 10−2
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

Zu, Y.; Wang, L.; Hu, Y.; Yang, G. CEEMDAN-LWT De-Noising Method for Pipe-Jacking Inertial Guidance System Based on Fiber Optic Gyroscope. Sensors 2024, 24, 1097. https://doi.org/10.3390/s24041097

AMA Style

Zu Y, Wang L, Hu Y, Yang G. CEEMDAN-LWT De-Noising Method for Pipe-Jacking Inertial Guidance System Based on Fiber Optic Gyroscope. Sensors. 2024; 24(4):1097. https://doi.org/10.3390/s24041097

Chicago/Turabian Style

Zu, Yutong, Lu Wang, Yuanbiao Hu, and Gansheng Yang. 2024. "CEEMDAN-LWT De-Noising Method for Pipe-Jacking Inertial Guidance System Based on Fiber Optic Gyroscope" Sensors 24, no. 4: 1097. https://doi.org/10.3390/s24041097

APA Style

Zu, Y., Wang, L., Hu, Y., & Yang, G. (2024). CEEMDAN-LWT De-Noising Method for Pipe-Jacking Inertial Guidance System Based on Fiber Optic Gyroscope. Sensors, 24(4), 1097. https://doi.org/10.3390/s24041097

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