Variational Specific Mode Extraction: A Novel Method for Defect Signal Detection of Ferromagnetic Pipeline

The non-contact detection of buried ferromagnetic pipeline is a long-standing problem in the field of inspection of outside pipelines, and the extraction of magnetic anomaly signal is a prerequisite for accurate detection. Pipeline defects can cause the fluctuation of magnetic signals, which are easily submerged in wide-band background noise without external excitation sources. Previously, Variational Mode Decomposition (VMD) was used to separate modal components; however, VMD is based on narrow-band signal processing algorithm and the calculation is complex. In this article, a method of pipeline defect signal based on Variational Specific Mode Extraction (VSME) is employed to extract the signal of a specific central frequency by signal modal decomposition, i.e., the specific mode is weak magnetic anomaly signal of pipeline defects. VSME is based on the fact that a wide-band signal can be converted into a narrow-band signal by demodulation method. Furthermore, the problem of wide-band signal decomposition is expressed as an optimal demodulation problem, which can be solved by alternating direction method of multipliers. The proposed algorithm is verified by artificially synthesized signals, and its performance is better than that of VMD. The results showed that the VSME method can extract the magnetic anomaly signal of pipeline damage using experimental data, while obtaining a better accuracy.


Introduction
Pipeline transportation is one of the most economical and effective means to transport oil and gas to different regions [1]. However, the ferromagnetic pipeline is extremely vulnerable to damage due to internal and external defects, third-party damage, and manufacturing defects when it is buried underground, resulting in oil and gas leakage. To prevent the danger of oil and gas leakage caused by pipeline damage, we need to check the pipeline regularly to ensure the safety and health of the pipeline. The routine method of long-distance pipelines inspection is the In-Line Inspection (ILI) [2], such as Ultrasonic Testing [3] and Magnetic Flux Leakage testing [4,5]. However, most of the small diameter pipes and urban underground pipelines have characteristics that restrict the piping, therefore non-contact detection of pipeline defects is an urgent problem to be solved. The use of non-contact approaches for inspection and damage detection is of great relevance. Not only this is indispensable for buried pipelines and other structural elements otherwise impossible to approach physically, but it also allows performing Structural Health Monitoring (SHM) without altering the structural element under investigation. This is, for instance, the case with video-based SHM and computer vision [6,7]. While not directly applicable to buried pipelines, this is doable for external estimation based on the modal Instantaneous Frequency (IF) and Instantaneous Amplitude (IA) spectrum. For the extraction of pipeline damage signals, the VMD algorithm still seems too redundant, because it needs to analyze all signal modal components before stopping calculation. Moreover, EMD and VMD can only separate modal components, but cannot directly extract mode characteristics such as IF and IA [32]. Variational Mode Extraction (VME) [35] is a low-frequency modal signal extraction algorithm, which has a similar mathematical basis as VMD. However, the algorithm is based on the assumption that the residual signal after extracting a specific mode should not contain energy at the mode center frequency.
We tried to find a more effective technique for MAD to make up for the shortcomings of existing detection methods. Inspired by VMD, NVCMD, and VME algorithms, we propose a novel Variational Specific Mode Extraction (VSME) method to analyze the pipeline defects signals which can easily be submerged in geomagnetic background noise. The proposed method has a reliable mathematical foundation as VMD, which is based on Wiener filtering, variational decomposition, and Hilbert transform [33]. In addition, VSME has another theoretical basis that a wide-band signal can be converted into a narrow-band signal by matching demodulation transform technique [36]. Furthermore, the problem of wide-band signal decomposition is expressed as an optimal demodulation problem, which can be solved by Alternative Direction Method of Multipliers (ADMM) [37]. Therefore, we convert the problem of pipeline defect signals extraction into a demodulation problem. Actually, the proposed algorithm can be regarded as extracting a Time-Frequency filter with a specific center frequency, and its characteristics eliminate the need to determine the number of signal modes. The optimization method can be briefly described as: for a specific mode, two demodulated orthogonal signals can be iteratively updated through two filters with the estimated instantaneous frequency; then, the estimated instantaneous frequency is further updated by the arc-tangent demodulation method which obtains phase of the two orthogonal signals. Firstly, we use the novel method on three artificial signals to evaluate its performance. Then, applying it to pipeline defect detection, the results show that has significantly improved calculation efficiency and accuracy compared with the EMD, VMD, and EWT method. The methods mentioned above can be taken as different filters, and have their respective characteristics, as summarized in Table 1. The article is organized as follows. In Section 2, the signal model is defined and related concepts are introduced. We propose the VSME method, including the main ideas and algorithm flow in Section 3. In Section 4, the experiment design is presented, and this method of reliability was verified by simulation and laboratory experiments. Finally, Section 5 concludes the paper.

Signal Model
The purpose of VSME is to extract pipeline defect signal mode (PDSM) from magnetic anomaly detection signal. Definition 1. The PDSM is essentially an Amplitude Modulated and Frequency Modulated (AM-FM) function, which can be expressed as where A k (t) > 0 and f (t) > 0 are the instantaneous amplitude (IA) and instantaneous frequency (IF) of signal, respectively. φ is the initial phase.
In this paper, the analytic signal z k (t) can be expressed by Hilbert transform as follows: where H represent Hilbert transform. The analytical signal z k (t) expressed by amplitude and frequency can be defined as: where ϕ(t) denotes the modulation source. Original signal can be represented by g k (t) = z k (t) , which means the component signal g k (t) can be expressed as functions of amplitude A k (t) and frequency f (t).

Bandwidth of the Specific Mode
As is known, there is no strict definition for the bandwidth of AM-FM signal. In this subsection, it is assumed that the signal function has a bandwidth limitation and follows Carson's rule. Thus, the bandwidth of nonlinear intrinsic mode function can be defined as [33]. Definition 2. The bandwidth of a mode component can be expressed as: where BW AM is the maximum frequency of IA resulting from the Amplitude Modulated (AM), and BW FM stands for the bandwidth caused by Frequency Modulated (FM). Theoretically, the IA has no bandwidth limitation, but we assume that we can find an BW AM , which includes most energy of IA A k (t). BW FM increases both with the maximum deviation ∆ f of the IF and with the rate of excursion f FM , which can be defined as:

Matching Demodulation Transform
The Hilbert transform cannot provide the orthogonal transform for the signal g k (t). Thus, we can get the approximate results in Equation (2) under some restrictive assumptions [38]. We can then define two operators, which are called Modulation Operator (MO) and Demodulation Operator (DO), as follows: Algorithms 2020, 13, 105 5 of 20 where f d (s) stands for frequency of the operators, and carrier frequency f c is a constant value [36]. Demodulation technology can be used to eliminate the FM of component g k (t), thus generating narrow-band demodulation mode component. The demodulated signal can be expressed as follows: where the frequency of demodulated signal g d k (t) is around the carrier frequency. The original signal g k (t) can be obtained by modulation with the MO, and the equation is as follow: It should be noted that we get mixed multi-component signals, and each component has its own DO. These DOs demodulate the signal into different forms, and their Time-Frequency (TF) representations are concentrated around the corresponding IF trace. Then, the demodulated signal is represented in the TF plane, i.e., a two-dimensional demodulated signal is converted into another two-dimensional function in TF domain. In this paper, we aim to extract specific mode components, and the extraction strategy is to extract from the corresponding TF plane based on its trajectory. The TF upper boundary and TF lower boundary of the TF subregion of the component g d Then, the TF subregion of k − th component:

Main Idea
The proposed algorithm is based on the premise that the wide-band signal can be converted into narrow-band signal, and matching demodulation transform technique is used to modulate each component mode to its corresponding position. The IF of signal z k (t) in Equation (3) is Modulation source ϕ(t) in Equation (3) can be expressed as In the process of signal demodulation, a bivariate function of time t and time-shift u is introduced as a backward MO to transform the mode component z k (t) into a bivariate function: Algorithms 2020, 13, 105 6 of 20 The one-dimensional signal, which only takes time as a function, is transformed into a two-dimensional bivariate function of time variable and time-shift variable. The IF of bivariate function can be expressed as follows: where f i,g k (t) is bivariate MO, and it suppresses the high order terms of ϕ(t). The instantaneous frequency of signal is the same as that of the corresponding time instant u, i.e., f i,g k (t) = f i,g (u). The demodulation process can be performed in two steps, as shown in Figure 1. In the first step, the original signal is demodulated to carrier signal A k (t)e i2π f c t by using forward DO e −iϕ(t) , i.e., the nonlinear instantaneous frequency is converted to the constant frequency. Therefore, the spectrum of demodulated signal is more concentrated than that of original signal. However, there is a problem in the demodulated spectrum, i.e., it cannot reflect the actual time-frequency characteristics of the original signal. Therefore, the demodulated signal needs to be further modulated into another pure carrier signal by the backward MO e −i[ϕ(t)−ϕ(u)−ϕ (u)(t−u)] to get a bivariate function about time t and time-shift u in the second step. The demodulation process can be performed in two steps, as shown in Figure 1. In the first step, the original signal is demodulated to carrier signal , i.e., the nonlinear instantaneous frequency is converted to the constant frequency. Therefore, the spectrum of demodulated signal is more concentrated than that of original signal. However, there is a problem in the demodulated spectrum, i.e., it cannot reflect the actual time-frequency characteristics of the original signal. Therefore, the demodulated signal needs to be further modulated into another pure carrier signal by the backward MO   can be obtained as: where the instantaneous amplitude can be given as Furthermore, the evaluation of signal bandwidth is expressed as: where the square root of 2-Norm of second-derivative is employed to evaluate the mode component bandwidth. For the highly nonlinear FM signal, we adopt an update strategy of real-time adjusting IF in the iterative process [39]. This method is different from the traditional recursive scheme [33,40], which can better separate the similar modes. As a whole, the procedure can estimate instantaneous frequency of mode component, with which the demodulated baseband signal has the narrowest frequency band. To evaluate the baseband signal bandwidth, we calculate the 2-Norm of signal derivative, and express the mixed signal g(t) as follows: (18) wheref k (t) stands for the frequency functions of demodulation operators, and u k (t) and v k (t) can be obtained as: where the instantaneous amplitude can be given as A k (t) = u k 2 (t) + v k 2 (t). Furthermore, the evaluation of signal bandwidth is expressed as: where the square root of 2-Norm of second-derivative is employed to evaluate the mode component bandwidth. For the highly nonlinear FM signal, we adopt an update strategy of real-time adjusting IF in the iterative process [39]. This method is different from the traditional recursive scheme [33,40], which can better separate the similar modes.

Algorithm
In this section, we introduce the procedure of variational specific mode extraction algorithm in detail. To enhance the robustness of this algorithm, we introduce auxiliary variables ξ, i.e., signals with noise. Then, we transform the problem in Equation (21) into: where we explicitly consider the influence of noise by adding a noise function ψ(ξ) in Equation (22), ensuring that this method can still be effective in noisy environment. To solve the above constrained optimization problem, augmented Lagrangian can be transformed into: where λ is a Lagrangian multiplier, and α is a quadratic penalty parameter. Then, we use alternating direction method of multipliers (ADMM), which is a distributed computing frame [37], to solve this problem. This method is suitable for solving distributed convex optimization problems, and it integrates the advantages of dual ascent method and method of multipliers. The purpose of ADMM is to decompose the original function and the extended function. When the objective function is complex, ADMM can divide the objective function into several sub-problems. The general expression of ADMM problem solving is: where b m is the sub-function decomposed by ADMM, and x m is its corresponding variable. For the improved ADMM to solve optimization details, we will describe in the following content. Firstly, VSME algorithm is to update the auxiliary variables, thus, the sub-problem to be solved is as follows: The second step is to update two orthogonal demodulation signals u k (t) and v k (t), as follow: Algorithms 2020, 13, 105 where Λ is the second-order difference operator. We use Equations (26) and (27) to update the corresponding mode components. Next, the instantaneous frequency increment is obtained by arc-tangent demodulation method, and it is updated by: Since the instantaneous frequency fluctuation is a smooth function, we can update the instantaneous frequency by: where 0 < χ < 1 is scale factor for stabilizing algorithm. We let χ = 0.5, the algorithm can provide good results. Finally, Lagrange multiplier λ is updated by the dual ascent algorithm as follows: Therefore, the final algorithm (Algorithm 1) can be summarized as follows: In this algorithm, we employed the same technique as [33] to solve boundary problem, i.e., the mirror extension of each side of the mode component is extended by half length. If the initial IFs are far away from the real IFs, the algorithm is easy to diverge. Therefore, we use restart method [41], which has been widely used in accelerating the process and stabilizing the numerical algorithm.

Performance Analysis
In this subsection, a theoretical analysis of VSME performance is provided, including quantitative analysis of instantaneous frequency estimation error and convergence conditions.
Firstly, the problem of instantaneous frequency quantization is discussed, ϕ (m) (t) is the phase function of the original signal g(t) in step m of iteration, and its corresponding demodulation operator Algorithms 2020, 13, 105 where phase ∆ϕ (m),u (t) is the remainder of the first-order Taylor expansion of estimated phase ϕ (m) (t) in the vicinity of u, as follows: The time-frequency energy density of a signal can be defined as: where T is the sampling interval, nT stands for the discrete sequence, s represents the Gaussian function, and ω is the angular frequency.
In the time-frequency dimension, the value of instantaneous angular frequency can be estimated as: General expression of instantaneous angular frequency estimation error: where δ ∆φ,∆ϕ is the error caused by the residual error of deviation, and ∂P (m) (t, ω)/∂ω 0 δ ε is caused by the noise. For the problem of convergence condition, the function expression of original signal phase ϕ(t) can be established by the determined parameter α, i.e., the phase function is ϕ(t, α). The m − th phase function r (m) (t) is recorded as ϕ (m) (t, α), and the deviation ∆ω (M) (t) of IF estimation error of the M − th component is smaller than other ∆ω(t) values on the spectrum. Thus, the estimated instantaneous frequency r (M+1) (t) = ϕ t, α (M+1) is closer to the true IF ϕ(t, α), i.e., the estimated parameter α (M+1) of the IF function ϕ t, α (M+1) is closer to the actual parameter α of the true IF function ϕ (t, α). In short, the implementation of iterative procedure will gradually match the true IF, and eventually converge.
In addition, increasing the variance of Gaussian function would reduce the variance of the IF estimation error, but would increase its bias. To solve this problem, the variance of Gaussian function is solved by the optimization problem of minimizing the mean square error of the IF estimation: where the variance σ(t) varies with the third-derivative of the phase function ϕ(t) and the (m) (t) is much smaller than the third-derivative ϕ (3) (t) of phase function, thus a time invariant variance of the window function can minimize the estimation error of the IF.

Results and Analysis
In this section, we concentrate on the application of VSME algorithm in pipeline defect signal extraction. As discussed above, variational mode theory provides a satisfactory solution to MAD problem. Therefore, both VMD and VSME algorithms can be applied to solve this problem. However, in the MAD problem, VSME is a better choice than VMD because central frequency of the expected mode component can be determined.

Method Validation
In this subsection, we used VSME to analyze a artificial signal, which is a intrawave amplitude modulated and frequency modulated signal. Artificial data generated from mathematical expression were calculated to evaluate the performance of VSME. All the calculation work was running in MATLAB R2016b.
The artificial signal contains a signal whose frequency increases linearly with time, a signal whose frequency decreases linearly with time, and a signal with time-varying frequency and amplitude, as shown in Equation (36).
For signal frequency extraction, we used VSME and RPRG (Ridge Path Regrouping) [42] algorithm to compare the extraction effect. RPRG method is a novel non-parametric algorithm to extract IF from signal with overlapping components in time-frequency representation. Figure 2a is an image of a mixed signal f (t), and Figure 2b shows the IF of original component signal f sig1 (t), f sig2 (t) and f sig3 (t). The IF of component signal extracted by RPRG is shown in Figure 2e. It is observed that at the intersection of IF of three component signals, an abnormal shape is generated. Using VSME to extract the IF of mode component, we can get a better result than RPRG, as shown in Figure 2c. To quantify the separation effect of VSME and RPRG on mixed signals f (t), Figure 2d,f show the absolute error between the estimated IF extracted by the two methods and the true IF. In this subsection, we used VSME to analyze a artificial signal, which is a intrawave amplitude modulated and frequency modulated signal. Artificial data generated from mathematical expression were calculated to evaluate the performance of VSME. All the calculation work was running in MATLAB R2016b.
The artificial signal contains a signal whose frequency increases linearly with time, a signal whose frequency decreases linearly with time, and a signal with time-varying frequency and amplitude, as shown in Equation (36).
For signal frequency extraction, we used VSME and RPRG (Ridge Path Regrouping) [42] algorithm to compare the extraction effect. RPRG method is a novel non-parametric algorithm to extract IF from signal with overlapping components in time-frequency representation. Figure 2a  To compare with different methods, we employed VSME and VMD to extract the mode components of original signal ( ) f t . The red line shown in Figure 3a represents the estimated three mode components, and the blue line represents the absolute error between the estimated modes and To compare with different methods, we employed VSME and VMD to extract the mode components of original signal f (t). The red line shown in Figure 3a represents the estimated three mode components, and the blue line represents the absolute error between the estimated modes and the original modes in the vicinity of plane y = 0. The absolute error range is less than 0.05, and the relative error is not more than 1%. By comparing Figure 3a,b, we can conclude that the VMD algorithm has a poor separation result on mixed signals f (t).
Algorithms 2020, 13, x FOR PEER REVIEW 11 of 20 the original modes in the vicinity of plane 0 y = . The absolute error range is less than 0.05, and the relative error is not more than 1%. By comparing Figure 3a,b, we can conclude that the VMD algorithm has a poor separation result on mixed signals ( ) In this paper, the performance of VSME is studied by adding white noise to the artificial signal ( ) f t . The signal-to-noise ratio (SNR) is defined as: ( ) 10 10 log where x P represents the energy of noiseless signal, and P ε is the energy of the noise signal. The mean value of ( ) t ε is 0, and its variance 2 =1 ε σ .
Figures 4 and 5 analyze the mode extraction performance of VSME algorithm for mixed signal ( ) f t , when the SNR is 12 dB and 18 dB respectively. The results show that the algorithm is effective when the SNR is 18 dB; when the SNR is lower than 12 dB, the IF has obvious fluctuation effect in the time domain, as shown in Figure 4c. We will continue to optimize the algorithm in the future work, aiming at extract mode components more accurately in low SNR condition. In this paper, the performance of VSME is studied by adding white noise to the artificial signal f (t). The signal-to-noise ratio (SNR) is defined as: SNR = 10 log 10 (P x /P ε ) (40) where P x represents the energy of noiseless signal, and P ε is the energy of the noise signal. The mean value of ε(t) is 0, and its variance σ 2 ε = 1. Figures 4 and 5 analyze the mode extraction performance of VSME algorithm for mixed signal f (t), when the SNR is 12 dB and 18 dB respectively. The results show that the algorithm is effective when the SNR is 18 dB; when the SNR is lower than 12 dB, the IF has obvious fluctuation effect in the time domain, as shown in Figure 4c. We will continue to optimize the algorithm in the future work, aiming at extract mode components more accurately in low SNR condition.
The root mean squared error (RMSE) is used to measure the filtered signal error under different SNR levels, which is defined as: whereĝ k (t) is k − th estimated signal component, and g k (t) is the true signal. The 2-Norm error expression is: where · is the 2-Norm of the function, andĝ k (t) is the k − th estimated component signal.
ε Figures 4 and 5 analyze the mode extraction performance of VSME algorithm for mixed signal ( ) f t , when the SNR is 12 dB and 18 dB respectively. The results show that the algorithm is effective when the SNR is 18 dB; when the SNR is lower than 12 dB, the IF has obvious fluctuation effect in the time domain, as shown in Figure 4c. We will continue to optimize the algorithm in the future work, aiming at extract mode components more accurately in low SNR condition.   The root mean squared error (RMSE) is used to measure the filtered signal error under different SNR levels, which is defined as: where  is the 2-Norm of the function, and ( ) k g t is the -k th estimated component signal. Figure 6 shows the calculated values of RMSE and 2-Norm of the estimated component signal under different SNR conditions. When the SNR is 12 dB, the values of RMSE is less than 0.02, and 2-  Figure 6 shows the calculated values of RMSE and 2-Norm of the estimated component signal under different SNR conditions. When the SNR is 12 dB, the values of RMSE is less than 0.02, and 2-Norm does not exceed 0.2; when the SNR is more than 20 dB, the values of RMSE and 2-Norm tend to decline steadily and are less than 0.05.
In this section, we simplify the problem by extracting the interested component of a simulated signal, and verify the reliability of VSME algorithm. It uses the same mode definition as VMD, especially using Wiener filter to detect the signal mode. The difference is that the number of modal components is determined before VMD calculation, and VSME only extracts the mode of interest through its center frequency. VMD uses iterative method to determine all modal components, and the calculation can be stopped until all modes reach the optimal response. Therefore, this property of VMD reduces the convergence rate of calculation. In addition, the performance of VMD is greatly impacted by the number of modal components, which is given at the beginning of calculation. However, VSME only searches for the component of interest, thus its computational efficiency is not affected by modal number. Without setting the property of specific mode in advance, VSME can extract the interested mode component from a mixed signal, such as the defect signal detection of buried pipeline. The results showed that VSME can achieve better performance and higher calculation efficiency.
frequency extracted by VSME.
The root mean squared error (RMSE) is used to measure the filtered signal error under different SNR levels, which is defined as: where  is the 2-Norm of the function, and ( ) k g t is the -k th estimated component signal. Figure 6 shows the calculated values of RMSE and 2-Norm of the estimated component signal under different SNR conditions. When the SNR is 12 dB, the values of RMSE is less than 0.02, and 2-Norm does not exceed 0.2; when the SNR is more than 20 dB, the values of RMSE and 2-Norm tend to decline steadily and are less than 0.05.

Experimental Setup
To verify the effectiveness of proposed VSME method for pipeline detection, an experiment was carried out using the experimental equipment shown in Figure 7 The experimental device consists of cross array with five magnetic sensor elements, a data acquisition card, a GPS positioning system, and a data acquisition system. Sampling frequency of the experimental instrument is 1 kHz, and each sensor is three-axis differential output. Testing equipment can locate these lines and accurately measure the size and specific location of all defects.
Algorithms 2020, 13, x FOR PEER REVIEW 13 of 20 In this section, we simplify the problem by extracting the interested component of a simulated signal, and verify the reliability of VSME algorithm. It uses the same mode definition as VMD, especially using Wiener filter to detect the signal mode. The difference is that the number of modal components is determined before VMD calculation, and VSME only extracts the mode of interest through its center frequency. VMD uses iterative method to determine all modal components, and the calculation can be stopped until all modes reach the optimal response. Therefore, this property of VMD reduces the convergence rate of calculation. In addition, the performance of VMD is greatly impacted by the number of modal components, which is given at the beginning of calculation. However, VSME only searches for the component of interest, thus its computational efficiency is not affected by modal number. Without setting the property of specific mode in advance, VSME can extract the interested mode component from a mixed signal, such as the defect signal detection of buried pipeline. The results showed that VSME can achieve better performance and higher calculation efficiency.

Experimental Setup
To verify the effectiveness of proposed VSME method for pipeline detection, an experiment was carried out using the experimental equipment shown in Figure 7 The experimental device consists of cross array with five magnetic sensor elements, a data acquisition card, a GPS positioning system, and a data acquisition system. Sampling frequency of the experimental instrument is 1 kHz, and each sensor is three-axis differential output. Testing equipment can locate these lines and accurately measure the size and specific location of all defects.
Array element adopts the triaxial sensor of Tunnel Magnetoresistance (TMR) produced by MultiDimension Technology, and the model is TMR2905. To prevent the acquisition system from interfering with the array, we use digital signal communication to transmit signals by two Universal Serial Bus cables. Data acquisition system is responsible for receiving signals for display, storage, and related control functions. Global Positioning System (GPS) is integrated into the experimental instrument for precise positioning and trajectory tracking. Its P3 mobile terminal and I70 carry out differential conversion of satellite signals through digital radio communication with accuracy up to centimeter level, which uses the GPS Fix Data in the current international NMEA0183 standard format for conversion and preservation. GPS data can be directly converted into Keyhole Markup Language format file through Google map API, which can directly display the actual map, and can mark the location of pipe body damage abnormally. The laboratory experimental pipeline is made of Q235 steel material with a wall thickness of 3 mm, and a diameter of 75 mm. We fabricated six defects at different locations of the pipe to study the magnetic anomaly signals of different defects, as shown in Figure 8.  Array element adopts the triaxial sensor of Tunnel Magnetoresistance (TMR) produced by MultiDimension Technology, and the model is TMR2905. To prevent the acquisition system from interfering with the array, we use digital signal communication to transmit signals by two Universal Serial Bus cables. Data acquisition system is responsible for receiving signals for display, storage, and related control functions. Global Positioning System (GPS) is integrated into the experimental instrument for precise positioning and trajectory tracking. Its P3 mobile terminal and I70 carry out differential conversion of satellite signals through digital radio communication with accuracy up to centimeter level, which uses the GPS Fix Data in the current international NMEA0183 standard format for conversion and preservation. GPS data can be directly converted into Keyhole Markup Language format file through Google map API, which can directly display the actual map, and can mark the location of pipe body damage abnormally.
The laboratory experimental pipeline is made of Q235 steel material with a wall thickness of 3 mm, and a diameter of 75 mm. We fabricated six defects at different locations of the pipe to study the magnetic anomaly signals of different defects, as shown in Figure 8. The laboratory experimental pipeline is made of Q235 steel material with a wall thickness of 3 mm, and a diameter of 75 mm. We fabricated six defects at different locations of the pipe to study the magnetic anomaly signals of different defects, as shown in Figure 8. The shapes of these six defects are axial crack, transverse crack, square groove, through hole, blind hole, and 45-degree crack in turn, as shown in Figure 9. There were four different lift-off distance detection experiments, which are 0.5 m, 1.0 m, 1.5 m and 2.0 m, respectively. The shapes of these six defects are axial crack, transverse crack, square groove, through hole, blind hole, and 45-degree crack in turn, as shown in Figure 9. There were four different lift-off distance detection experiments, which are 0.5 m, 1.0 m, 1.5 m and 2.0 m, respectively. The shapes of these six defects are axial crack, transverse crack, square groove, through hole, blind hole, and 45-degree crack in turn, as shown in Figure 9. There were four different lift-off distance detection experiments, which are 0.5 m, 1.0 m, 1.5 m and 2.0 m, respectively.

Results Analysis
Buried pipelines are affected by external factors, thus the algorithm for extracting defect signals should have a certain anti-noise performance, and it should also have the ability to extract abnormal signal of subtle change. Therefore, we compare the performance of band-pass filter (BPF), EWT, VMD, and VSME for pipeline defect signals.
Before extracting the pipeline defect signal in original signal, as shown in Figure 10a, we first use STFT to carry out time-frequency spectrum analysis of the signal. Figure 10b shows that the main energy of acquired signal is concentrated in the low-frequency component, and the frequency spectrum changes abnormally in the middle of detected pipeline. According to the actual pipeline being tested, as shown in Figure 9, the position of through hole defect in pipeline is just detected at about 10th second of constant speed detection, where pipeline damage is relatively serious. This situation is usually not effective for the BPF method, and its response contains some additional fluctuations, which increase with the increase of noise level. Figure 10d-f show that analysis results of EWT, VMD, and VSME are satisfactory, and it is difficult to identify the subtle differences manually. In the following analysis, we use quantitative indicators to measure the filtering effect.

Results Analysis
Buried pipelines are affected by external factors, thus the algorithm for extracting defect signals should have a certain anti-noise performance, and it should also have the ability to extract abnormal signal of subtle change. Therefore, we compare the performance of band-pass filter (BPF), EWT, VMD, and VSME for pipeline defect signals.
Before extracting the pipeline defect signal in original signal, as shown in Figure 10a, we first use STFT to carry out time-frequency spectrum analysis of the signal. Figure 10b shows that the main energy of acquired signal is concentrated in the low-frequency component, and the frequency spectrum changes abnormally in the middle of detected pipeline. According to the actual pipeline being tested, as shown in Figure 9, the position of through hole defect in pipeline is just detected at about 10th second of constant speed detection, where pipeline damage is relatively serious. This situation is usually not effective for the BPF method, and its response contains some additional fluctuations, which increase with the increase of noise level. Figure 10d-f show that analysis results of EWT, VMD, and VSME are satisfactory, and it is difficult to identify the subtle differences manually. In the following analysis, we use quantitative indicators to measure the filtering effect.
Considering the length of article, only one group data of analysis figures was listed, and other groups of data gave the calculation results. To analyze the accuracy of mentioned algorithms in extracting MAD signal, frequency band of magnetic anomaly signal is obtained from the defect signal of experimental pipeline, and then detection results are analyzed statistically. We use SNR, normalized cross correlation coefficient (NCCC), and root mean square error (RMSE) to evaluate the performance difference between VME detectors and other techniques. The calculation equations of the two indexes are in Equations (40) and (41), and the NCCC is as follows: Algorithms 2020, 13, 105 16 of 20 wherex(n) is filtered output signal, x(n) is input signal, and N is signal dimension. The three mentioned indexes were calculated for Figure 11, and the average value is shown in Table 2.
use STFT to carry out time-frequency spectrum analysis of the signal. Figure 10b shows that the main energy of acquired signal is concentrated in the low-frequency component, and the frequency spectrum changes abnormally in the middle of detected pipeline. According to the actual pipeline being tested, as shown in Figure 9, the position of through hole defect in pipeline is just detected at about 10th second of constant speed detection, where pipeline damage is relatively serious. This situation is usually not effective for the BPF method, and its response contains some additional fluctuations, which increase with the increase of noise level. Figure 10d-f show that analysis results of EWT, VMD, and VSME are satisfactory, and it is difficult to identify the subtle differences manually. In the following analysis, we use quantitative indicators to measure the filtering effect. Considering the length of article, only one group data of analysis figures was listed, and other groups of data gave the calculation results. To analyze the accuracy of mentioned algorithms in extracting MAD signal, frequency band of magnetic anomaly signal is obtained from the defect signal of experimental pipeline, and then detection results are analyzed statistically. We use SNR, normalized cross correlation coefficient (NCCC), and root mean square error (RMSE) to evaluate the performance difference between VME detectors and other techniques. The calculation equations of the two indexes are in Equations (40) and (41), and the NCCC is as follows: where ( ) x n is filtered output signal, ( ) x n is input signal, and N is signal dimension. The three mentioned indexes were calculated for Figure 11, and the average value is shown in Table 2.  mentioned three filtering indicators shows that the VSME algorithm has good performance in the extraction of pipeline defect signals.  Calculation time of the same data set is an important index to measure algorithm. To compare the computational complexity of different algorithms, average time required to extract the MAD signal from the 100,000 sample points of experimental data is given in Table 3. According to the table, the calculation time of VSME is less than one-tenth that of VMD and EWT, and it can be compared with the computational complexity of BPF. Results showed that VSME has a good performance by evaluating the index of filtered experimental data. Table 3. Mean and standard deviation (SD) of calculation time.  To further evaluate the detection effect of pipeline defects at different lift-off distances, the SNR, RMSE, and NCCC of defect signals were computed using these algorithms. Figure 11 plots the evaluation indicators for the six sets of data, which are obtained at different lift-off distances. SNR of extracted defect signal gradually decreases as the lift-off distance increases, as shown in Figure 11a. The performance of VSME is slightly higher than VMD and EWT, but significantly higher than BPF. RMSE of the VSME-processed signal is generally lower than the other three methods, indicating that the proposed algorithm is more reliable, as shown in Figure 11b. When the lift-off distance increased from 0.5 m to 2 m, the value of RMSE showed a trend of increasing first and then decreasing. Through experimental analysis, when the lift-off distance is 1 m-1.5 m, the excessive RMSE fluctuation is due to the joint effect of the pipeline magnetic field signal and the surrounding interference, which leads to the instability of the magnetic field. As the lift-off height increases, the pipeline magnetic field signal weakens, thus the value of RMSE decreases. Figure 11c shows the NCCC, the performance of BPF is the worst, and VSME is slightly better than the other two methods. The analysis of the mentioned three filtering indicators shows that the VSME algorithm has good performance in the extraction of pipeline defect signals.
Calculation time of the same data set is an important index to measure algorithm. To compare the computational complexity of different algorithms, average time required to extract the MAD signal from the 100,000 sample points of experimental data is given in Table 3. According to the table, the calculation time of VSME is less than one-tenth that of VMD and EWT, and it can be compared with the computational complexity of BPF. Results showed that VSME has a good performance by evaluating the index of filtered experimental data. Computational time for processing a segment of 100,000 samples (100 s) of acquired signal (in second). Computing platform adopts DELL-EQVS5M4, the hardware configuration is Intel (R) core (TM) i5-4430 CPU @ 3.0GHz and the memory is 8.0 GB. The system version is Windows 10 professional, and we use Matlab 2016b as signal processing software.
Comparing defect signals at different lift-off distances, we transformed the extracted MAD signal into a normalized gradient energy operator. By analyzing the test data of pipeline in laboratory, it is found that with the increase of the lift-off distance, the signal of the less damaged location becomes weaker. Six defects are easily detected when the lift-off distance is 0.5 m, as shown in Figure 7b. When lift-off distance is 2.0 m, six defects can still be detected. However, the signal with larger damage degree has a certain submergence effect on the position with smaller damage degree. Figure 12d shows that the signals are weak at the axial crack and 45-degree crack positions. In this section, we concentrate on the performance of VSME in extracting MAD signals of pipeline defects. As mentioned, variational theory can provide effective solutions for pipeline defect signal extraction. Herein, we used and VSME to extract pipeline defect signals, and compared with BPF, EWT, and VMD. To ensure the performance of MAD signal extraction, we randomly selected six sets of experimental data, and extracted the MAD of defect signals by four algorithms on the same computing platform for several times. Then, the acquired MAD signals were evaluated by calculating the cross correlation coefficients with the reference pipeline defect signals. The results showed that VSME has certain advantages.

Conclusions
In the paper, we present a novel algorithm of extracting a specific mode component for noncontact inspection of pipeline defects. The proposed VSME is based on the fact that a wide-band signal can be converted into a narrow-band signal by demodulation method. Then, the problem of wide-band signal decomposition is expressed as an optimal demodulation problem, which can be solved by an improved ADMM. This algorithm does not need a parametric model, which makes it more suitable for practical problems. To verify the performance of VSME, we used this algorithm to analyze a complex simulation signal. Simulation results showed that VSME is more effective than VMD and EWT in extracting a specific mode signal, even in the case of modes overlap severely. Finally, we applied this method to the actual pipeline defect detection, and compared the detection results at different lift-off distances, which showed that it has a better anti-noise performance. Further work will be focused on improving the robustness of VSME in separating severely crossed component signals and on the applications of pipeline detection in more complex environment.   In this section, we concentrate on the performance of VSME in extracting MAD signals of pipeline defects. As mentioned, variational theory can provide effective solutions for pipeline defect signal extraction. Herein, we used and VSME to extract pipeline defect signals, and compared with BPF, EWT, and VMD. To ensure the performance of MAD signal extraction, we randomly selected six sets of experimental data, and extracted the MAD of defect signals by four algorithms on the same computing platform for several times. Then, the acquired MAD signals were evaluated by calculating the cross correlation coefficients with the reference pipeline defect signals. The results showed that VSME has certain advantages.

Conclusions
In the paper, we present a novel algorithm of extracting a specific mode component for non-contact inspection of pipeline defects. The proposed VSME is based on the fact that a wide-band signal can be converted into a narrow-band signal by demodulation method. Then, the problem of wide-band signal decomposition is expressed as an optimal demodulation problem, which can be solved by an improved ADMM. This algorithm does not need a parametric model, which makes it more suitable for practical problems. To verify the performance of VSME, we used this algorithm to analyze a complex simulation signal. Simulation results showed that VSME is more effective than VMD and EWT in extracting a specific mode signal, even in the case of modes overlap severely. Finally, we applied this method to the actual pipeline defect detection, and compared the detection results at different lift-off distances, which showed that it has a better anti-noise performance. Further work will be focused on improving the robustness of VSME in separating severely crossed component signals and on the applications of pipeline detection in more complex environment.