Next Article in Journal
Wind Resource Assessment and Economic Viability of Conventional and Unconventional Small Wind Turbines: A Case Study of Maryland
Next Article in Special Issue
An Analytical Model for Production Analysis of Hydraulically Fractured Shale Gas Reservoirs Considering Irregular Stimulated Regions
Previous Article in Journal
Design of High-Speed Permanent Magnet Motor Considering Rotor Radial Force and Motor Losses
Previous Article in Special Issue
Numerical Investigation on Proppant–Water Mixture Transport in Slot under High Reynolds Number Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Processing for EM Telemetry Signal Recovery: Field Data from Sichuan Province

1
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
Frontier technology and Equipment Development Center for Deep Resources Exploration, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
3
Innovation Academy for Earth Science, Chinese Academy of Sciences, Beijing 100029, China
4
College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(22), 5873; https://doi.org/10.3390/en13225873
Submission received: 15 September 2020 / Revised: 5 November 2020 / Accepted: 6 November 2020 / Published: 10 November 2020
(This article belongs to the Special Issue Development of Unconventional Reservoirs 2020)

Abstract

:
This paper deals with the study of multi-channel adaptive noise cancellation with a focus on its application in electromagnetic (EM) telemetry. We presented new variable step-size least mean square (LMS) techniques: regularized variable step-size LMS and regularized sigmoid variable size LMS, for electromagnetic telemetry data processing. Considering the complexity and spatial distribution of environmental noise, algorithms with multiple reference signals were used to retrieve transmitted EM signals. The feasibility of the regularized variable step size LMS algorithms with numerical simulation was analyzed and presented. The adaptive processing techniques were applied in the recovery of frequency and binary phase shift key modulated signal. The proposed multi-channel adaptive technique achieves fast convergence speed, low mean squared error and is shown to have good convergence characteristics compared to conventional methods. In addition to attaining good results from the multi-channel adaptive filter and performing the signal analysis in real-time, we implemented combined fast effective impulse noise removal techniques. The combination of median and mean filters was effective in removing a wide range of impulsive noises without distorting any other data points. Further, electromagnetic telemetry data were acquired during a drilling operation in Sichuan province, China, for real field application. Data processing workflow was designed for EM telemetry data processing based on the noise characteristics, simulation results and expected result for demodulation. To establish a comprehensive overview, a performance comparison of the acquisition array system is also provided. Conclusively, the introduced multichannel adaptive noise canceling techniques are very effective in recovering transmitted EM telemetry signals.

1. Introduction

Over the past decade, bi-directional transmission of data from Bottom Hole Assembly (BHA) to the rig floor has played a vital role in reducing the risk involved in geosteering, well landing and saving considerable drilling time. In electromagnetic telemetry (EMT), downhole measured data are transmitted as coded electromagnetic waves by installing an electrical insulator sub along the drill string. Surface receivers are used in measuring the voltage difference between the drilling string and an earth antenna fixed into the ground at tens to hundreds of meters away from the rig. One of the fundamental challenges affecting the performance of EMT is the impact of noise on the measured signal. When the transmitted signals are from a transceiver at great depths, EM telemetry signal strength measured at the surface is very low due to attenuation, and noisy due to majorly surface interference. The noise could either distort or mask the transmitted signal and are caused by system nonlinearity, noise interference from adjacent environment and more.
Geophysical signal processing is either carried out post-data recovery or pre-data recovery from receiver unit, where a trade-off between time and data processing quality is considered. Digital signal processing (DSP) techniques use ranges from simple data filtering and smoothening by moving the average method to a more complex workflow in the neural network method [1,2,3,4,5,6,7,8,9]. However, recent practices have relied on high-resolution information to be applied or accessed in real-time. Conventional communication systems use filters to cancel noise in a received signal. Based on the previous study [5,6] and confirmation from our findings, noise sources on or near the drill rig usually have large power spectral densities over a broad frequency band including the transmitting frequency. Very low signal to noise ratio data, with noise and main signals operating in the same frequency, require highly effective data processing techniques. This explains why there are many difficulties associated with the retrieval of transmitted current signals from recorded EMT signals. Therefore, to maximize the acquired information and optimize the data recovery procedure, a significant and efficient data processing system is essential for accurate analysis and signal interpretation.
Two classes of noise affect the performance of electromagnetic telemetry systems: coherent harmonic noises from mainly industrial power supply equipment, which are simple to filter out either by spectral analysis or notch filtering and the random noises from atmospheric activities (e.g., thunderstorms), power grid switching, telluric, pipe scraping, etc. Therefore, EMT noise can be characterized as a combination of both random and impulsive noise. The random impulsive noises are mostly created by the power grid switching and by thunderstorms. In comparison to related fields, limited study relating to EMT (EM MWD) data processing has been published. Therefore, part of the review will reference other fields with similar signal-noise characteristics.
Noise suppression techniques have been introduced in several signal processing studies with similar noise characteristics, where the magnitude of noise recorded is higher than recorded signals and or data with very low signal to noise ratio is recovered. This includes adaptive filtering techniques such as least mean square, recursive least square, normalized least mean square, multichannel acoustic echo canceler, nonlinear adaptive filtering and other noise reduction techniques, spectral subtraction method and neural network techniques [5,6,10,11,12,13,14].
However, most of these methods achieve better noise cancellation only under certain conditions or by combining several denoising strategies; for example, spectral subtraction method used by reference [6] in EM MWD data processing resulted in improvement of the signal to noise ratio in the presence of typical rig surface noise. However, the improvement was limited to effective suppression of receiver-filtering related noise, with limitations faced in the retrieving of the coded signal. Reference [8], in their study of the effectiveness of Wiener filtering and adaptive noise canceling techniques in multi-channel magnetic resonance sounding (MRS) signal processing, showed/implied that different filters must be used for effective noise cancellation of different noise events as the two techniques performed less efficiently with signals containing spikes. In addition, neither methods reacted to the random phase shift or the difference in sampling frequencies, which led to a reduction in the noise canceling efficiency.
Furthermore, the study on least mean square (LMS) and normalized LMS (NLMS) techniques in noise cancellation shows that these techniques work well for medium to high signal to noise ratio signals but results in either noisy recovered signals which are post-processed or over-processed signal with information lost [12,15,16] when the signal to noise ratio ranges from low to very low. Reference [17] demonstrated that data collected from different sites show the line noise increases as the site falls within the power grid. In contrast, the extremely low frequency (ELF) noise power changes slowly as the site location is changed. This phenomenon is highly favorable to the application of an adaptive method. However, this study is focused on the application of single-channel LMS filter. With an improvement in adaptive noise cancellation algorithms, we present a comparison of the noise cancellation effect from 4 different adaptive methods with a multi-channel approach.
In general, the performance of nonlinear adaptive filtering algorithms is better than linear adaptive filtering algorithms [18,19,20]. However, the computation cost and the number of hardware required for these techniques are higher. Therefore, for practical engineering application, linear adaptive filtering algorithms such as LMS algorithm, whose calculations are smaller and hardware requirements are lower, are frequently used [21]. However, this paper aims at presenting the study on the application of nonlinear multichannel adaptive noise attenuation techniques and novel fast processing workflow in EMT to detect EM signals at the well surface.
In addition, to attain the best results from the multi-channel adaptive filter and perform the signal analysis in real-time, we introduced a combined fast and effective impulse noise removal technique. An ideal filter for impulsive noise removal is one that removes only the impulses without distorting any other data points.
When impulsive noise is present in digital data, median filters are often used as a remedy. The median filters run a sliding window over the data and at each step replace the point in the middle of the window by the median of the points inside the window [22,23,24,25]. However, study has shown the median filter technique to introduce distortion in the recovered data [26]. A more recent instantaneous approach involves applying a modified Z score to locate spikes and a simple moving average filter to remove the located spikes [26]. The Z scores are calculated using the detrend differentiated series. However, the technique is only effective at annihilating localized linear and slow-moving curve linear trends. Other techniques, such as generalized and modified wavelet methods [24], and spectral method [27,28] have been used effectively, but at the expense of obtaining timely results. Likewise, several techniques have been effectively applied in the attenuation of salt and pepper noise, which include new switching based median filter (NSMF) [29], adaptive weight median filter (AWMF) [30,31] and iterative adaptive fuzzy filter using alpha-trimmed mean [32]. References [33,34] demonstrated the application of maximum correntropy criterion-based sparse adaptive filtering algorithms in eliminating large impulsive noise. [33] indicates that the exponential term in the MCC updating equation is responsible for the elimination of large errors, which turns the algorithm robust against impulsive-like noises as observed in the pattern of the MSD training curves. Reference [35] introduced a technique that combines multilevel weighted graphs technique with induced generalized order weighted average (IGOWA) to remove salt and pepper noise from images. In addition, these techniques are highly effective in removing salt and pepper noise but are less time-efficient as they employ an iterative process in filter weight estimation.
As part of the developmental stage in geophysical equipment design, signal retrieval using effective noise cancellation tools are considered in this manuscript, with focus on the adaptive noise removal and inclusion of combined impulse noise removal technique. In this study, we compare results obtained from different adaptive noise cancellation techniques (NLMS, SVSLMS) with multi-reference array techniques and two new forms of sigmoidal variable step-size adaptive processing techniques. Signals generated from frequency shifting and phase-shifting form of signal modulation techniques (FSK and PSK) were considered. We considered the improvements in the algorithms for adaptive noise canceling techniques. Several reference dipoles, physically displaced from the primary dipole, are used to describe the noise characteristics in the drilling area. The high degree of noise correlation between distant locations makes adaptive noise canceling a good choice for noise removal. In addition, an adaptive filtering system can modify its parameters to adjust to the different noise conditions.

2. General Theory for Adaptive Noise Canceling

The adaptive noise-canceling technique has been an essential tool in the field of signal processing. Signal retrieval using effective noise cancellation tools, with a focus on the nonlinear adaptive noise removal, has a very wide application, including borehole studies. For EM telemetry, the data recorded at the surface, as shown in (1), are described as a combination of input current/voltage source Sv with specified waveform, earth filter Se and noise such as industrial noise, cultural noise, telluric and external/extraterrestrial noise Sn (receiver function excluded).
S d = S v + S e + S n ,
where Sd is the signal recorded at the receiver. In additive characteristics, the earth response and other noise characteristics are combined, and the simplified noise elimination expression is given as
S v = S d S n
However, we cannot achieve optimal filtering signal with fixed coefficient filter when random, unpredictable signal and noise characteristics are recorded. Therefore, nonlinear adaptive noise cancellation technique is essential. This technique requires the measurement of at least a single reference datum point that will adequately define the noise characteristics Sn for the specified dataset Sd. To collect the environmental noise without the transmitted signal being present, the reference dipole is placed colinear but away from the primary dipole and close to noise sources as much as possible. However, a single reference channel circumstance can only be achieved if there is only one and uniformly distributed noise source. This could be less effective in the situation of complexity and non-homogeneous spatial noise distribution. Therefore, our study will focus on using multiple reference signals.
In general, the ANC system has two parts: a primary channel with the main signal, and reference signals referred to as the noise signal (Figure 1). The recorded primary signal is represented as a time-harmonic signal by
S d ( n ) = A s · cos ( ω t + n 1 )
n =   [ 1 , 2 , 3 , , m ] , while the reference signals are given as
S ( n ) = [ S 1 ( n ) , S 2 ( n ) , S 3 ( n ) , , S k ( n ) ]
S k ( n ) = A r · cos ( ω t + n k )
and their respective filter weights are given as
W i ( n ) =   [ W 1 ( n ) , W 2 ( n ) , W 3 ( n ) , , W k ( n ) ]
such that the average reference signal is given as
S a ( n ) = i = 1 K [ W i ( n ) S i ( n ) ]
Based on (2), the estimated error signal can be written as follows:
E [ e ] =   E [ S v + S n S a ]
To adjust the parameters in a linear adaptive system, the least mean square (LMS) algorithm is used as a fast search algorithm using a gradient estimate [15,36]. The basic idea of the LMS algorithm is that it makes mean square error between the derived noise signal and desired signal to be minimum by adjusting the weight coefficients. This is expressed as
E | e ( n ) | 2 =   E | S v ( n ) + S n ( n ) S a ( n ) | 2
Consequently, the general expression for calculation of weight coefficients for the different LMS technique is given as,
W i ( n ) =   W i ( n 1 ) + 2 μ ( n ) e ( n ) S i ( n )
where μ(n) is the step size and e(n) is the estimated error signal. For a nonlinear adaptive system, the filtering process requires continuous adjustment of the average reference signal with the filter weight coefficient. Its application efficiency is owned to the ability to estimate the noise characteristics from provided reference data.

2.1. Fixed and Step-Size LMS Algorithms

Four classes of adaptive algorithms are considered in this paper. The adaptive algorithms are: least mean square (LMS), normalized LMS (NLMS), sigmoid variable size LMS (SVSLMS) [37,38], and regularized variable step size (VSS) techniques: regularized variable step size LMS (RVSSLMS) and regularized sigmoid variable size LMS (RSVSLMS). The expression for the calculation of the weight coefficients for the different LMS techniques to be considered is given as,
LMS:
W ( n ) =   W ( n 1 ) + 2 μ e ( n ) S d ( n )
that is, μ is constant.
NLMS:
W i ( n ) =   W i ( n 1 ) + 2 ( α β + S d T ( n ) S d ( n ) ) e ( n ) S i ( n )
where 0 < β < 1 , and 1 < α < 10 .
β in (12) is a small positive constant used to avoid division by zero when the input vector Sd is zero and α is the adaptation positive constant.
SVSLMS:
W i ( n ) =   W i ( n 1 ) + 2 ( β · min ( a b s ( e ( n ) ) · [ 1 exp ( α | e ( n ) | 2 ] , [ 1 exp ( α | e ( n ) | 2 ] ) ) e ( n ) S i ( n )
where 0 < β < 1 λ max , and 1 < α < 10 .
α is a parameter controlling the shape of the sigmoid function, while β controls the range of the sigmoid function.
RVSSLMS:
W i ( n ) =   W i ( n 1 ) + 2 ( β · min ( a b s ( e ( n ) ) · [ 1 exp ( α | e ( n ) | 2 ] + ( λ [ | μ ( n 1 ) | 2 ] ) ) e ( n ) S i ( n )
where 0 < β < 1 λ max , 1 < α < 10 , and 0.1 < λ < 1 .
λ is a parameter controlling the trade-off between the shape of sigmoid function and range of the sigmoid function.
RSVSLMS:
W i ( n ) =   W i ( n 1 ) + 2 ( ( β λ ) · min ( a b s ( e ( n ) ) · [ 1 exp ( α | e ( n ) | 2 ] + ( ( β / λ ) · [ 1 exp ( α | e ( n ) | 2 ] ) ) e ( n ) S i ( n )
where 0 < β < 1 λ max , 1 < α < 10 , and 0.1 < λ < 1 .
β in (13) to (15) controls the convergence rate, and both α and λ control the error when the algorithm was almost steady. In addition, λ controls the error at the initial state for the RSVSLMS technique. This ensures a possible shift in the function characteristics to ensure a faster convergence at the early stage.
LMS is the most commonly used technique due to its simplicity. For each iteration, the LMS algorithm requires 2N additions and 2N+1 multiplications. However, using a constant step size, the filter is modified by an amount relative to the immediate estimate of the gradient of the error. This could lead to very slow convergence and eventual settling for a less accurate result, especially when dealing with random noise. The normalized least mean square algorithm (NLMS) technique is an extension of the conventional LMS technique which solves the above-mentioned issue by selecting a different step size value during the iteration process. In the NLMS algorithm, the step size varies as a function of the input energy and hence gives enhanced convergence properties [10,14]. The value of β in (12) is a small positive constant, and it is introduced in the algorithm to avoid division by zero when the values of the input vector are zero. The robustness against low SNR could be improved by reusing coefficient schemes [39,40], particularly the case for non-stationary (time-varying) noise signals, and where a sudden rise in the noise power may occur. SVSLMS algorithm, on the other hand, is a variable step-size LMS algorithm based on sigmoid function and are controlled by two constants, α and β. The two constants control the shape and size of the step size function, respectively, and in-turn the error in the steady state and the convergence rate. The condition for the convergence of this algorithm is 0 < μ ( n ) < 1 / ε max , where ε max is the maximum eigenvalue of the autocorrelation matrix of input signals [10]. Since β controls the value range and initial value of the step size, its maximum limit is the same as the inverse of the maximum eigenvalue of the autocorrelation matrix of input signals. The algorithm has a faster convergence rate, but μ(n) varies too fast when e(n) approaches zero.
Similarly, in regularized VSS techniques, RVSSLMS and RSVSLMS, α affects the shape, while λ acts to regularize both the shape and the initial step size, and in turn the rate of convergence. This guides against instability in SVSLMS, which is related to the choice of β (too high or too low) with respect to the given initial value.

2.2. Analysis of Regularized VSS Techniques’ Efficiency

For the analysis of the efficiency of the RSVSLMS technique and selection of appropriate value for λ, we presented plots of the change of step size with a change in error function at various values of α and λ (Figure 2a–d). From the plot, we deduce that the value of λ within the range of 0.1–1 gives an effective convergence rate with the value of α ranging from 1 to 10 and effective value of 5.
Figure 2 indicates the shape of μ(n) changes a lot with a change in the values of λ. A rapid convergence rate is expected with the obtained function curve at the initial adaptive phase, meaning big value of μ(n). However, as e(n) closes to 0, there is a small detuning at the later phase, meaning μ(n) changes slowly. In general, the following description describes the influence of λ on the shape of the function plot, vis-a-vis the convergence rate;
λ ≡ 0.1: The rate of change of the step size at the initial phase is slow, keeping the step size high and the convergence rate fast.
λ ≡ 1: The rate of change of the step size at the initial phase is higher, having a quick transition from the value of high step size to a low value and the convergence rate reduced. This is expected to be beneficial in case of the moderately noisy signal.
In extreme cases,
λ < 0.1: The step size at the initial phase is constant and very high, keeping the step size high and the convergence rate fast.
λ > 1: The rate of change of the step size is fast and almost linear, and there is smaller detuning at the steady-state phase that means μ(n) changes slowly when e(n) closes to 0. The function curve is similar to the RVSSLMS curve (Figure 3). Based on its characteristics at the steady-state phase, it is expected to be less efficient in wide range noise canceling when α is set higher than 5.
Therefore, for the comparison of the function curve from other SVSLMS algorithm, we set α and λ to be 5 and 0.5, respectively, considering the demands of larger initial step size and a steady change when e(n) is close to 0. Figure 3 shows the curves of the step size, μ(n) changing with the error function, e(n), adopting methods from the modified SVSLMS algorithm (blue dot/scattered plot) for comparison [38]. The proposed RSVSLMS algorithm has a large but slowly varying step size in the initial phase to ensure rapid convergence rate, and after convergence, the step size varies slowly in a steady-state phase. The update parameter, λ, has a stable value for the initial step size, which is determined by β, within the specified effective range, while the shape changes with change in λ, thereby increasing or reducing the rate of change of the step size to accommodate the difference in noise level, convergence rate and time. When large, the update parameter λ controls both the amplitude and shape of the step size if its range is wider. Therefore, β and α are kept constant while only the update parameter is changed for utmost performance, ensuring a perfect trade-off between the function curve shape and magnitude of the step size.
RVSSLMS function curve, on the other hand, has a constant shape, characterized by a relatively linear decay of step size. However, in comparison with the variable step size algorithm presented by Reference [41], there is a small detuning at the steady-state phase as e(n) closes to 0, meaning μ(n) changes slowly within the steady-state phase, thereby resulting in steady convergence.

3. Implementation and Simulation

In this section, synthetic data with defined noise characteristics, running as an infinitely continuous signal are denoised by implementing adaptive filtering. The results are presented in the section for a period of 0.1 s. The design for the adaptive filter simulation in developed Matlab Simulink for real-time noise suppression using adaptive filters; NLMS, RVSSLMS and RSVSLMS, are shown in Appendix A. This generates synthetic data with random noise. The interpolated (digital FIR) low pass filter was introduced to generate the variation in noise characteristics by varying the cut of frequency between 50 Hz and 350 Hz and introducing delay time to the signal to represent reference noise (Figure 4d). Continuous moving average was used to generate the variation in noise characteristics to represent reference noise acquisition (Figure 4c). The design allows an operator to change parameters at any time in the design to achieve optimization filter, and observe and compare the results from designed filters for the adaptive techniques with LMS or NLMS filter using a built-in Matlab function.
We illustrate the performance of the proposed algorithms by comparing the output of the computer simulations in a system with an identical scenario and technique parameters given in Table 1. The adaptive filter and the unknown channel are assumed to have the same number of taps. The input signals are obtained by filtering the coded signal with: (1) random sequence and (2) random sequence and Gaussian noise (Figure 4d). Two reference noises with a time lag of 0.1 s and 0.25 s were added to S(n), which are independent of the input signal.
Firstly, the frequency shift key (FSK) modulation was selected in data coding with a sample rate of 2400 Hz and a cycle period of 0.05 s. The coded signal with a carrier frequency of 60 Hz is shown in Figure 4b. The input signals are obtained by filtering the coded signal with; (1) random sequence and (2) random sequence and Gaussian noise. Two reference noise with a time lag of 0.1 s and 0.25 s were added to S(n), which are independent of the input signal.
The left-hand side plots in Figure 5 shows the comparison of the overlay of the error signal on the input signal for all the techniques. These show the convergence speeds of the newly introduced algorithms (RVSSLMS and RSVSLMS) are higher than those of the old ones (NLMS and SVSLMS). Conversely, considering a measured signal, with random noise having the same sample time as the recorded signal, the proposed algorithms are faster, more accurate and efficient in noise cancellation than the regular varying step size NLMS algorithm. After a short period, this algorithm convergence enters a steady-state, and its detuning is smaller, which means that the fluctuation/step size changed within a smaller range. It also can be seen that the performance of SVSLMS is limited as the rate of convergence of the technique is slower, and it becomes unstable if the initial step size is reduced by using a value of β = 0.000001.
Furthermore, we carry out a comparison of the performance between the algorithms in the presence of random sequence with a sample time of 0.001 s and Gaussian noise with a sample time of 0.1 s. Both primary channel and reference channels contain a certain ratio of Gaussian noise running at a sample time different from that of the main signal. For a better representation of observed field scenarios, we added impulsive noise to the reference signal (Figure 6) and included the impulse noise reduction technique in the noise reduction workflow.
Later in this paper, we applied a two-stage method to suppress the impulsive noise in the field data but right now, only the second stage is applied to the simulated data. However, the complete two-stage method is introduced. Thus, in the first phase, we used the mean filter to remove generalized more repetitive impulsive signals. This is done by following the given steps.
Let Y = Y1, …, Yj, …, Yn, represent the values of the data within a given time window. Y j and Y n represent the mid-time and the end-tine data for the given time window, respectively. The mean M n t mean(Y) of the series is calculated.
The obtained mean M n t is compared with previous mean M n t 1 , and the lowest of both is selected as M n t .
The difference between the mean value ( M n t ) and mid-time value ( Y j ) is then calculated; Z m = M n t Y j . Z m value is then compared against a set of upper and lower thresholds, which represents the thresholds for background noise in the area which are acquired prior to the start of the field data.
Interpolated value Y j is then obtained for the mid-time data by calculating the mean of its immediate neighbors, using the following mathematical expression
Y j = 1 n i = 1 n Y i × I ( | Z m | > τ l < τ h )
where I ( u ) is an indicator function taking value one if the condition u is satisfied and 0 otherwise, and Y j remains unchanged. The width of the moving average neighborhood is controlled by the parameter i; which was set in our applications to 8.
In the second phase, we implemented the modified Z median filter [35]. However, the update for the interpolated data considers both the upper and lower thresholds using (8) above. This technique is effective at annihilating localized linear and slow-moving curve linear trends. The method proposed here requires low computation time and still reduces the total error considerably. This algorithm is computationally efficient and inexpensive, easy to execute and should be of great utility in a wide range of impulse noise removal. In this section, the criterion | Z m | > 1 was used as a guideline based on an outlier-labeling rule with the objective of screening large datasets. In practice, the scientist will have immediate control over this threshold parameter.
Figure 6 shows the synthesized main signal with added impulsive noise at a time of 2.25 s, the combined reference signal with added noise at time 3.75 s, and their respective denoised signals. The frequency spectrum for the signals (Figure 6a) indicates that the largest noise amplitude is concentrated within the first 50 Hz, which is within the range of the transmitted signal and that combined reference signal record higher amplitude at higher frequencies. The most notable difference between noisy and filtered data is observed in the combined references at high frequencies. In addition, no notable spike is recorded within the range of frequencies displayed. However, the effectiveness of the modified Z mean technique in removing localized impulsive noise and the effect of this processing step is established in this section.
Figure 7 demonstrates the processing results obtained from the different adaptive techniques. By considering the recovered signal from the impulsive noisy signal, Figure 7c,d shows the introduced algorithms (RVSSLMS and RSVSLMS) at a given time range performs well and better than the traditional varying step size NLMS algorithm and SVSLMS algorithm. The signal waveform extracted by the RSVSLMS algorithm is the closest to the ideal signal waveform with less noise. Although the result obtained with the NLMS technique reflects a much lesser distortion from the introduced impulsive noise, the recovered signal, in general, is very noisy and shows no good representation of the transmitted data at the given time frame.
The effectiveness of the modified Z median filter is observed by comparing the results obtained from noisy data with those obtained after the impulse noise has been filtered out (Figure 8a,b), using the RSVSLMS technique. The effect of impulse noise is observed as a distortion in the recovered transmitted signal at times 2.25 s, 3.75 s and 3.84 s. However, these distortions were effectively corrected, and the true transmitted waveform was restored in the recovered transmitted signal from median filtered data. In addition, the impulsive noise was more pronounced in Figure 8a where the noise was added to the main signal than in Figure 8b where the impulsive noise was added to only one of the reference signals at 3.75 s. The application of multi-channel reference signals resulted in the reduction of the effect of the impulsive noise until it became relatively insignificant. However, when the impulsive signal exists in all the reference signals at 3.84 s, the effect is notable, so also is the effectiveness of the median filter.
Furthermore, the binary phase shift key (BPSK) modulation was selected in data coding with a sample rate of 2400 Hz and a cycle period of 1 s. The coded signal with a central frequency of 10 Hz is shown in Figure 9b. The selected modulation technique and parameter is similar to the field transmitted current/voltage signal from the bottom hole assembly to the surface but with periodic cycle and larger sample rate. We illustrate the performance of the new algorithms by comparing the output of the computer simulations in a given system, using the parameters given in Table 2. The adaptive filter and the unknown channel are assumed to have the same number of taps.
Figure 10 shows comparisons of the denoised signal obtained by using LMS, SVSLMS, RVSSLMS, RSVSLMS techniques over a time range of 4 s. Figure 10b–d shows the introduced algorithms and the traditional varied step size LMS algorithm performed very well and better than the traditional LMS algorithm. The signal waveform extracted by these algorithms were close to the ideal signal waveform in Figure 9b. Figure 11 shows the comparison of the mean squared error plot for the various adaptive techniques. The mean squared error plots show that the result obtained from the variable step size techniques were close, and are more accurate than the conventional LMS technique. A high convergence rate characterizes the LMS technique result, but with low accuracy and a steady mean squared error value of 2.
Amongst the variable step size techniques, the RSVSLMS technique had the highest convergence rate with a mean squared error value of about 0.2. The traditional SVSLMS technique, on the other hand, had the largest convergence rate with a significant initial drop in MSE error for a period of about 10 s. However, the SVSLMS shows a continuous steady decrease in MSE error over a long period at the steady-state phase. This results in the lowest MSE value of about 0.01 for this study. However, in EM telemetry, a trade-off between convergence rate and MSE value is required to ensure real-time processing of the transmitted signal.
Similarly, we simulated synthetic signal using a combination of random and white noise but with a frequency difference of range of 10, that is, 10 and 20 Hz, respectively. Similar to previous results (Figure 10), the LMS technique had the highest convergence rate but with a very low MSE value, which is indicative of noisy retrieved signal as observed in Figure 12. The MSE plot (Figure 13) also shows a close similarity in the accuracy of the retrieved signal from the RSVSLMS and RVSSLMS techniques. In comparison with results from less complex noise characteristics, the retrieved signal quality for LMS and SVSLMS techniques decreased significantly, and their convergence rate increased by a large margin. Thus, we can infer that with more complex noise data from field measurement, the above techniques would be less effective in real-time EM telemetry signal recovery than introduced RVSSLMS and RSVSLMS techniques.
In general, the obtained results show that the newly introduced techniques are very effective in complex noise suppression. The proposed algorithms are faster, more accurate and efficient than the regular varying step size NLMS algorithm within the short time frame expected for EM telemetry signal recovery. After a short period, these algorithms convergence enters a steady-state, and its detuning is smaller, which means that the fluctuation changes within a smaller range. The rapid and then steady change in step size and changes in function curve for RSVSLMS algorithm ensures the stability of the introduced techniques over a wide range of applications. Although the computation cost and the number of hardware required for these techniques are higher, considering the convergence rate and accuracy of the results obtained within the first 5 to 10 s, the newly introduced RVSSLMS technique is expected to be most effective in adaptive signal retrieval for EM telemetry. The fast convergence rate makes this technique applicable in real-time engineering application with high output data quality irrespective of the complexity of the noise characteristics. Lastly, the multi-channel adaptive noise canceling performed above on the synthetic data gave an insight into the characteristics of the retrieved phase shift signal.

4. Field and Simulation Study

Electromagnetic telemetry data were acquired during a drilling operation in Sichuan province, China. The starting depth of the downhole transceiver was about 2000 m. Phase shift coded signal with a frequency of 10 Hz was transmitted from the insulating gap along the drill string, which was outside the borehole casing. The electrode array used in EM telemetry signal measurement is shown in Figure 14. The voltage difference between dipoles (L1 to L10), which varies in length from 30 m to 80 m was recorded using the DRU 45 system developed by Institute of Geology and Geophysics, Chinese Academy of Sciences at a sampling frequency of 2400 Hz.
Firstly, we considered the performance of the techniques if the frequency shift signal had to be the modulation type by convolving the field noise signal with the generated frequency shift signal. The noise spectrum from dipoles L2, L4, L5 and L7, were used in this study (Figure 15a). The spectrum shows noise signals with an amplitude range of the order of 10-1 to 10-5 mv for the frequency range of interest, with the largest magnitude recorded by L2. The processed main signal is obtained by convolving the noise signal from the dipole L5 (blue) with simulated frequency shift key modulated signal (Figure 15b). For this study, lines L2, L4 and L7 were used as reference signals. The frequency shift signal was generated using a sample rate of 2400 Hz, signal frequencies of 10 and 5 Hz and a period of 0.5 sec. The transmitted signal amplitude was varied between 0.001 V and 0.0002 V to allow for an assessment of the processing technique under different signal to noise ratio (SNR). The plot of the amplitude spectrum of the generated data and combined reference noise within the frequency range of interest is presented in Figure 15c. The frequency spectrum (Figure 15c) indicates that the combined reference signal amplitude is generally higher, especially at the low frequencies. In addition, a significantly high amplitude spectrum was recorded around frequencies 200 Hz, 400 Hz and 600 Hz. The performance of the variable step size techniques was compared by analyzing the output (obtained error signal) signals over specified time ranges, the mean squared error plots and bar chart of SNR and the power of the output signal in dB.
Figure 16 and Figure 17 demonstrates the processing results obtained by the variable step size techniques. The obtained results show that the introduced algorithms (RVSSLMS and RSVSLMS) performed better than the traditional varying step size LMS algorithm at both early time and late time in this study. The inclusion of the impulsive noise filter in the workflow resulted in the removal of major impulsive noise signals as recorded between the time band of 135 s and 139 s (Figure 15b and Figure 16a). This is also evident in the MSE plot.
The signal waveform extracted by the SVSLMS algorithm still contained the footprint of high and low-frequency noise, even though the signal recovered still indicate the alternating form of the input signal in Figure 16b is close to the ideal signal waveform. The recovered signal from the regularized techniques, considering RSVSLMS technique as an example, show a decrease in quality with a decrease in signal to noise ratio from −55.5 dB to 69.4 dB. At the current state, a good correlation with the original signal was observed when the SNR of the input/main signal is higher than −69 dB. The power of the output signal, on the other hand, remain mostly constant and unhinged by the change in the main signal’s SNR (Figure 18). The recovered signal from RSVSLMS technique has the highest signal power with a value of −25 dB. From the results, it is evident that the proposed regularized LMS techniques achieve better noise tracking by providing an overall quality improvement of the processed synthetic EMT signal over the regular SVSLMS method.
Furthermore, the mean squared error plots show a very fast initial convergence for the regularized LMS techniques and a much slower convergence for the SVSLMS technique (Figure 19). However, after the first 170 s, the SVSLMS technique records the lowest MSE value until all the techniques converge at a later time. The smoothness of the regularized LMS technique’s MSE values represents a consistent, smooth and well denoised output signal. On the other hand, the MSE plot for the SVSLMS technique is more erratic, and it represents a noisy output signal. In all, the RSVSLMS technique with very fast convergence, the lowest long-range MSE value and highest recovered signal power, is expected to perform the best in EMT data processing with the combined impulse noise filtering included in the processing workflow.
Lastly, the acquired EMT data was processed using the newly introduced adaptive techniques. The noise characteristics for the area and individual dipole were considered in selecting the reference dipole, while the main signal was selected based on expected maximum telemetry signal from previous simulation studies. The noise spectrum from dipoles L1 to L7, are shown in Appendix B. The spectrum shows general noise signals with an amplitude range of the order of 10−1 to 10−5 v.
For the frequency range of interest (Appendix B), the measured voltage ranges from 10−1 to 10−3 v with the largest magnitude recorded by L3 for electrodes set colinear to the main dipole L1, and L7 for electrodes orthogonal to the main dipole. This shows that the magnitude of data noise recorded increases as the offset electrode moves closer to the drilling rig, especially with electrodes set orthogonal to the main dipole. The ELF noise power changes slowly as the site location is changed. This phenomenon is highly favorable to the application of an adaptive method. Further, several high amplitude signals were recorded over the broad frequency spectrum, which includes 50 Hz industrial/power noise. The main input signal is obtained from dipole L1 with a dipole length of 80 m. Several sets of reference noise were used to identify the most effective dipole array for the estimation of overall main data noise. However, only the results from the sets of reference data, set 1: L2, L8 and L9, and set 2: L4, L5 and L7, will be presented in this publication.
Figure 20 demonstrates the data processing workflow adopted in this study considering the noise characteristics, simulation results and expected result for phase shift demodulation. The primary input to the signal processor is the signal received from the well by the antenna(s) connected to the wellhead. The primary input contains the desired EMT signals plus all the noise that is to be filtered out. The secondary inputs are used as the references for the adaptive filter and contain only the surface noise.
To achieve the best results and still perform the signal processing in real-time, the number of data points for a given time frame is reduced by setting the new sampling rate between 512 Hz and 200 Hz, which is larger than the anti-aliasing filter frequency. These settings provide the system with enough information to process the data with acceptable accuracy in real-time. Equiripple (digital FIR) low pass filter was introduced before the impulsive noise canceling to reduce the excessive frequency bandwidth to a reasonable range with a crossover frequency of 80 Hz.
When impulsive noise is present in digital data, combined impulsive noise filter runs a 4 point sliding window over the data and at each step replaces the data with new values. The data are scanned and only changed if impulses are detected. The impulse data are replaced by a piecewise-linear estimated value calculated from previous and current data points. The used techniques meet the requirement of an ideal filter for impulsive noise removal, which is one that removes the impulses without distorting any other data points. In addition, the adaptive noise counselling result and time-continuous sine wave, with the same frequency as the transmitted signal (10 Hz), were combined to generate the final processing output to enhance the phase shifts for easy demodulation process.
The performance of the regularized variable step-size LMS techniques were compared by analyzing the output signals over specified time ranges and the mean squared deviation plots. Figure 21, Figure 22, Figure 23 and Figure 24 demonstrate the processing results obtained by the selected techniques. The recovered signal between 50 to 60 s and 80 to 90 s of the data is presented for comparison. Figure 21 and Figure 22 show that the introduced techniques are highly effective in EMT data processing even in the presence of complex field noise with wide frequency band. The results obtained from the regularized techniques are very similar to each other. However, the signal waveform extracted by the RVSSLMS algorithm still contained the footprint of low-frequency noise. This may result in masking of the phase shift signal in the final processing output even though the signal recovered is close to the ideal signal waveform.
Furthermore, the results obtained from single and multiple reference signals using the RSVSLMS technique are compared in Figure 23. Though the two results are similar in most cases, few differences such as the high amplitude distortion between 137 s and 138 s, and distortions between the time 132 s and 134 s caused by high-frequency signal interference shows the single reference technique could present less accurate result over time. In another study, the effect of change in the reference signals’ orientation was considered. The result obtained from the different sets of reference signals is presented in Figure 24. Result obtained from the RSVSLMS technique shows a significant decrease in accuracy when the orthogonal dipole signals were used for data processing. Therefore, the selection of reference signal array in EM telemetry data processing when phase shift modulation is used should be limited to colinear dipoles and dipoles within 300 of the main signal.
Lastly, the mean squared deviation (MSD) plots show a significant correlation between the processing output and the MSD value (Figure 25). The two techniques, RSVSLMS in green and RVSSLMS in blue, converged significantly differently from the conventional SVSLMS and NLMS techniques. The MSD of NLMS displays fast convergence but has the highest MSD value at an average value above 10−4. The SVSLMS technique, on the other hand, had the slowest convergence rate with lowest MSD value about 10-5. The MSD plot for RSVSLMS and RVSSLMS results display the fastest convergence and the lowest MSD value. The MSD plot becomes steady at an average value of about 10−7. The characteristics of the MSD plot is in agreement with the level of accuracy of signal retrieved from these techniques. In all, the RSVSLMS technique shows more consistent good performance in noise-canceling and MSD value at the steady-state. The correlation between the performance of the regularized LMS techniques, both with synthetic and real field data, signifies its effectiveness in EM telemetry signal processing even when the recorded data have a low signal to noise ratio. The obtained results also show that the newly introduced technique can adapt effectively, even in the presence of several changes in noise characteristics. In comparison to the phase modulated signal, the adaptive processed results from the frequency modulated signal, in general, show a readily definite representation of the input data.

5. Conclusions

Novel variable step-size algorithms, regularized SVSLMS and VSSLMS, were proposed in this study for EM telemetry signal processing. Analysis of the simulation results and MSE error show an increase in convergence rate and steady-state detuning at the later phase. The function curve and error curve illustration show that the regularized variable step-size LMS algorithms combine the virtues of a trade-off between the shape of the step size function curve and the step size initial value, which leads to fast convergence rate and lower final misalignment. The regularized SVSLMS and VSSLMS algorithms show more advantages in convergence rate during the initial phase and steady-state detuning than a traditional variable step NLMS algorithm and SVSLMS algorithm. Through comparative analysis between the proposed algorithms with multi-channel reference noise signal and previous LMS techniques, it is concluded that the new algorithms perform very well on signal waveform retrieval and data fitting error. The electromagnetic telemetry test data demonstrate that a good denoising effect can be obtained with multiple reference channels and also, the dipole orientation has a significant effect on the shift signal recovery.
The effectiveness of the designed regularized LMS algorithm and reference dipole arrays to extract the EM telemetry signal from field data with complex noise estimated based on the progressive weight calculation from multiple colinear and near colinear reference channels was established. The regularized variable step size approach significantly improved the separation performance level and could increase the reach of EM telemetry. The proposed processing workflow with combined impulsive noise filter and the multi-channel adaptive technique method can be used in developing timely computerized automating EMT signal retrieval. Based on our findings, further study on the application of neural network for further processing and effective demodulation is proposed.

Author Contributions

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

Funding

This work was supported in part by the Strategic Priority Research Program of the Chinese Academy of Sciences (NO. XDA140501000).

Acknowledgments

The Author thank will like to thank Du from the Institute of Geology and Geophysics, Chinese Academy of Sciences for his helpful suggestion.

Conflicts of Interest

The authors declare no conflict of interest. In addition, the funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Figure A1. Simulink diagram for noise cancellation using normalized least mean square technique (NLMS).
Figure A1. Simulink diagram for noise cancellation using normalized least mean square technique (NLMS).
Energies 13 05873 g0a1
Figure A2. Simulink diagram for multi-channel noise cancellation using the regularized variable step-size LMS (RVSSLMS).
Figure A2. Simulink diagram for multi-channel noise cancellation using the regularized variable step-size LMS (RVSSLMS).
Energies 13 05873 g0a2
Figure A3. Simulink diagram for multi-channel noise cancellation using the regularized sigmoid variable step-size LMS (RSVSLMS).
Figure A3. Simulink diagram for multi-channel noise cancellation using the regularized sigmoid variable step-size LMS (RSVSLMS).
Energies 13 05873 g0a3

Appendix B

Figure A4. EM telemetry noise characteristics for dipoles L1, L2 and L3: (a) frequency spectrum of measured noisy data, and (b) frequency spectrum of recorded data for the frequency range of interest.
Figure A4. EM telemetry noise characteristics for dipoles L1, L2 and L3: (a) frequency spectrum of measured noisy data, and (b) frequency spectrum of recorded data for the frequency range of interest.
Energies 13 05873 g0a4
Figure A5. EM telemetry noise characteristics for dipoles L4, L5 and L7 orthogonal to the main signal: (a) frequency spectrum of measured noisy data and (b) frequency spectrum of recorded data for the frequency range of interest.
Figure A5. EM telemetry noise characteristics for dipoles L4, L5 and L7 orthogonal to the main signal: (a) frequency spectrum of measured noisy data and (b) frequency spectrum of recorded data for the frequency range of interest.
Energies 13 05873 g0a5

References

  1. IEEE. Fifty Years of Signal Processing: The IEEE Signal Processing Society and its Technologies 1948–1998; IEEE: Piscataway, NJ, USA, 1998. [Google Scholar]
  2. Jiang, C.; Lin, J.; Duan, Q.; Sun, S.; Tian, B. Statistical stacking and adaptive notch filter to remove high-level electromagnetic noise from MRS measurements. Near Surf. Geophys. 2011, 9, 459–468. [Google Scholar] [CrossRef]
  3. Engelsen, S.B.; Van Den Berg, F.W. Quantitative Analysis of Time Domain NMR Relaxation Data. Mod. Magn. Reson. 2017, 1–19. [Google Scholar] [CrossRef]
  4. Jayant, A.; Singh, Y.; Kaur, M. Different Techniques to Remove Baseline Wander from ECG Signal. Int. J. Emerg. Res. Manag. Technol. 2013, 2, 16–19. [Google Scholar]
  5. Suh, A.E. Noise Cancellation in Electromagnetic Measurement While Drilling Using Spectral Subtraction. Master’s Thesis, California Polytechnic State University, San Luis Obispo, CA, USA, 2004. [Google Scholar]
  6. Whitacre, T.P. A Neural Network Receiver for EM-MWD Communication. Master’s Thesis, California Polytechnic State University, San Luis Obispo, CA, USA, 2016. [Google Scholar]
  7. Dalgaard, E.; Auken, E.; Larsen, J.J. Adaptive noise cancelling of multichannel magnetic resonance sounding signals. Geophys. J. Int. 2012, 191, 88–100. [Google Scholar] [CrossRef]
  8. Dalgaard, E.; Christiansen, A.V.; Larsen, J.J.; Auken, E. A temporal and spatial analysis of anthropogenic noise sources affecting SNMR. J. Appl. Geophys. 2014, 110, 34–42. [Google Scholar] [CrossRef]
  9. Cannistraci, C.V.; Abbas, A.; Gao, X. Median Modified Wiener Filter for nonlinear adaptive spatial denoising of protein NMR multidimensional spectra. Sci. Rep. 2015, 5, 8017. [Google Scholar] [CrossRef] [Green Version]
  10. Paul, T.K.; Ogunfunmi, T. A Kernel Adaptive Algorithm for Quaternion-Valued Inputs. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 2422–2439. [Google Scholar] [CrossRef]
  11. Islam, T.; Shahnaz, C.; Fattah, S. Speech enhancement based on a modified spectral subtraction method. In Proceedings of the 2014 IEEE 57th International Midwest Symposium on Circuits and Systems (MWSCAS), College Station, TX, USA, 3–6 August 2014; pp. 1085–1088. [Google Scholar]
  12. Kumar, A.; Agarwal, R.P. Adaptive Signal Processing: A Comparative Approach. Int. J. Contemp. Res. Eng. Technol. 2019, 9, 90–93. [Google Scholar]
  13. Kumar, A. Mamta Comparision of Different Types of IIR Filters. Int. J. Adv. Res. Electron. Commun. Eng. 2016, 5, 393–402. [Google Scholar]
  14. Yasukawa, H.; Shimada, S.; Furukawa, I. Acoustic echo canceller with high speech quality. In Proceedings of the 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Dallas, TX, USA, 6–9 April 1987; pp. 2125–2128. [Google Scholar]
  15. Kumar, A.; Goel, P.; Gupta, V.K.; Chandra, M. Comparative Research on Various Adaptive Algorithms for NoiseCancellation in Speech Signals. In Proceedings of the IEEE International Conference on Control, Computing, Communication and Materials, UIT, Allahabad, India, 21–22 October 2016; pp. 1–5. [Google Scholar]
  16. Jebastine, J.; Rani, B.S. Implementation of Block Least Mean Square Adaptive Algorithm for effective Noise Cancellation in Speech Signal. Int. J. Electr. Electron. Eng. Res. 2011, 1, 1–11. [Google Scholar]
  17. Kher, R. Signal Processing Techniques for Removing Noise from ECG Signals. J. Biomed. Eng. 2019, 1, 1–9. [Google Scholar]
  18. Ao, W.; Xiang, W.-Q.; Zhang, Y.-P.; Wang, L.; Lv, C.-Y.; Wang, Z.-H. A New Variable Step Size LMS Adaptive Filtering Algorithm. In Proceedings of the 2012 International Conference on Computer Science and Electronics Engineering, Hangzhou, China, 23–25 March 2012; pp. 265–268. [Google Scholar]
  19. Sharma, N.; Singh, M.K.; Low, S.Y.; Kumar, A. Weighted Sigmoid-Based Frequency-Selective Noise Filtering for Speech Denoising. Circuits Syst. Signal Process. 2020, 1–20. [Google Scholar] [CrossRef]
  20. Akhtar, M.T. On Active Impulsive Noise Control (AINC) Systems. Circuits Syst. Signal Process. 2020, 39, 4354–4377. [Google Scholar] [CrossRef]
  21. Yin, L.; Yang, R.; Gabbouj, M.; Neuvo, Y. Weighted median filters: A tutorial. IEEE Trans. Circuits Syst. II Express Briefs 1996, 43, 157–192. [Google Scholar] [CrossRef]
  22. Chan, R.H.; Wa, C.; Nikolova, M. Salt-and-pepper noise removal by median-type noise detectors and detail-preserving regularization. IEEE Trans. Image Process. 2005, 14, 1479–1485. [Google Scholar] [CrossRef]
  23. Gupta, G. Algorithm for image processing using improved median filter and comparison of mean, median and improved median filter. Int. J. Soft Comput. Eng. (IJSCE) 2011, 1, 2231–2307. [Google Scholar]
  24. Patel, A.X.; Kundu, P.; Rubinov, M.; Jones, P.S.; Vértes, P.E.; Ersche, K.D.; Suckling, J.; Bullmore, E.T. A wavelet method for modeling and despiking motion artifacts from resting-state fMRI time series. NeuroImage 2014, 95, 287–304. [Google Scholar] [CrossRef] [Green Version]
  25. David, C. Stone Application of median filtering to noisy data. Can. J. Chem. 1995, 73, 1573–1581. [Google Scholar]
  26. Whitaker, D.; Hayes, K. A simple algorithm for despiking Raman spectra. Chemom. Intell. Lab. Syst. 2018, 179, 82–84. [Google Scholar] [CrossRef]
  27. Yu, G.; Niu, S.; Ma, J. A hybrid spectral gradient method for removing salt-and-pepper impulse noise. In Proceedings of the 2011 4th International Congress on Image and Signal Processing, Shanghai, China, 15–17 October 2011; pp. 765–768. [Google Scholar]
  28. Liu, J.; Li, S. Spectral gradient method for impulse noise removal. Optim. Lett. 2015, 9, 1341–1351. [Google Scholar] [CrossRef]
  29. Patanavijit, V.; Thakulsukanant, K. Simulated evaluation of new switching based median filter for suppressing SPN and RVIN. Indones. J. Electr. Eng. Comput. Sci. 2019, 15, 688–696. [Google Scholar] [CrossRef]
  30. Khan, S.; Lee, D.-H. An adaptive dynamically weighted median filter for impulse noise removal. EURASIP J. Adv. Signal Process. 2017, 2017, 67. [Google Scholar] [CrossRef] [Green Version]
  31. Zhang, P.; Li, F. A New Adaptive Weighted Mean Filter for Removing Salt-and-Pepper Noise. IEEE Signal Process. Lett. 2014, 21, 1280–1283. [Google Scholar] [CrossRef]
  32. Ahmed, F.; Das, S. Removal of High-Density Salt-and-Pepper Noise in Images with an Iterative Adaptive Fuzzy Filter Using Alpha-Trimmed Mean. IEEE Trans. Fuzzy Syst. 2013, 22, 1352–1358. [Google Scholar] [CrossRef]
  33. Haddad, D.B.; Petraglia, M.R.; Petraglia, A. A unified approach for sparsity-aware and maximum correntropy adaptive filters. In Proceedings of the 2016 24th European Signal Processing Conference (EUSIPCO), Budapest, Hungary, 29 August–2 September 2016; pp. 170–174. [Google Scholar]
  34. Ma, W.; Qu, H.; Gui, G.; Xu, L.; Zhao, J.; Chen, B. Maximum correntropy criterion based sparse adaptive filtering algorithms for robust channel estimation under non-Gaussian environments. J. Frankl. Inst. 2015, 352, 2708–2727. [Google Scholar] [CrossRef] [Green Version]
  35. Xu, Q.; Zhang, Q.; Hu, D.; Liu, J. Removal of Salt and Pepper Noise in Corrupted Image Based on Multilevel Weighted Graphs and IGOWA Operator. Math. Probl. Eng. 2018, 2018, 7975248. [Google Scholar] [CrossRef]
  36. Akhtar, M.; Abe, M.; Kawamata, M. A new variable step size LMS algorithm-based method for improved online secondary path modeling in active noise control systems. IEEE Trans. Audio Speech Lang. Process. 2006, 14, 720–726. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, L.; Wang, Q.; Bai, H.P. An Improved Variable Step Size LMS Algorithm Based on Sigmoid Function. Appl. Mech. Mater. 2014, 602, 3593–3596. [Google Scholar] [CrossRef]
  38. Tian, B.-F.; Zhou, Y.-Y.; Zhu, H.; Jiang, C.; Yi, X.-F. Noise cancellation of a multi-reference full-wave magnetic resonance sounding signal based on a modified sigmoid variable step size least mean square algorithm. J. Cent. South Univ. 2017, 24, 900–911. [Google Scholar] [CrossRef]
  39. Kim, S.-E.; Lee, J.; Song, W.-J. Steady-state analysis of the NLMS algorithm with reusing coefficient vector and a method for improving its performance. In Proceedings of the 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Prague, Czech Republic, 22–27 May 2011; pp. 4120–4123. [Google Scholar]
  40. Resende, L.C.; Haddad, D.B.; Petraglia, M.R. A Variable Step-Size NLMS Algorithm with Adaptive Coefficient Vector Reusing. In Proceedings of the 2018 IEEE International Conference on Electro/Information Technology (EIT), Rochester, MI, USA, 3–5 May 2018; pp. 181–186. [Google Scholar]
  41. Liu, K. Adaptive noise cancellation for electromagnetic-while drilling system. In Proceedings of the 2016 3rd International Conference on Information Science and Control Engineering (ICISCE), Beijing, China, 8–10 July 2016; pp. 1253–1256. [Google Scholar]
Figure 1. A diagram of the layout of multi-channel adaptive noise canceling with reference signals having three noise sources. The primary channel receives the main signal, Sd, recorded with the second electrode positioned tens to hundreds of meters away from the borehole. Reference channels record signals Sn1, Sn2 and Sn3 from different locations colinear with the main signal, close to the borehole and has a higher intensity of environmental and/or industrial noise interference. Weight coefficients w1(n), w2(n) and w3(n) are adjusted by adaptive filters whose purpose is to produce an output, Sa(n), which is the sum of the corresponding outputs of each adaptive filter. If Sa(n) is as similar as possible to the noise characteristics in Sn(n) by the minimum mean square error, by subtraction of Sa(n) from the primary channel, the noise can be canceled to produce the initially transmitted signal Si(n).
Figure 1. A diagram of the layout of multi-channel adaptive noise canceling with reference signals having three noise sources. The primary channel receives the main signal, Sd, recorded with the second electrode positioned tens to hundreds of meters away from the borehole. Reference channels record signals Sn1, Sn2 and Sn3 from different locations colinear with the main signal, close to the borehole and has a higher intensity of environmental and/or industrial noise interference. Weight coefficients w1(n), w2(n) and w3(n) are adjusted by adaptive filters whose purpose is to produce an output, Sa(n), which is the sum of the corresponding outputs of each adaptive filter. If Sa(n) is as similar as possible to the noise characteristics in Sn(n) by the minimum mean square error, by subtraction of Sa(n) from the primary channel, the noise can be canceled to produce the initially transmitted signal Si(n).
Energies 13 05873 g001
Figure 2. Function curves of μ(n) changing with e(n) at different values of λ: (a) 0.01, (b) 0.5, (c) 1 and (d) 2.
Figure 2. Function curves of μ(n) changing with e(n) at different values of λ: (a) 0.01, (b) 0.5, (c) 1 and (d) 2.
Energies 13 05873 g002
Figure 3. Function curves of step size ‘μ(n)’ changing with error value ‘e(n)’ for three variable step size algorithms. The value of the parameters λ and α was set as 1 and 5, respectively, for regularized variable step size lean mean square (RVSSLMS) algorithm (green dot plot). For RSVSLMS algorithms, α was set at a constant value of 5, while λ was varied between 0.5 and 1 for RSVSLMS 2 and 1, respectively.
Figure 3. Function curves of step size ‘μ(n)’ changing with error value ‘e(n)’ for three variable step size algorithms. The value of the parameters λ and α was set as 1 and 5, respectively, for regularized variable step size lean mean square (RVSSLMS) algorithm (green dot plot). For RSVSLMS algorithms, α was set at a constant value of 5, while λ was varied between 0.5 and 1 for RSVSLMS 2 and 1, respectively.
Energies 13 05873 g003
Figure 4. (a) Plots of signal modulator waveform; (b) frequency shift key (FSK) generated signal for data transmission; (c) generated input signal representing recorded main data; and (d) generated weighted reference signal representing the noise characteristics.
Figure 4. (a) Plots of signal modulator waveform; (b) frequency shift key (FSK) generated signal for data transmission; (c) generated input signal representing recorded main data; and (d) generated weighted reference signal representing the noise characteristics.
Energies 13 05873 g004
Figure 5. Comparison plots of generated synthetic data 1 with random noise processed by multi-channel noise cancellation techniques: (a) normalized LMS; (b) sigmoid variable step size LMS (SVSLMS); (c) regularized sigmoid variable step size LMS (RVSSLMS); (d) regularized sigmoid variable step size LMS (RSVSLMS). On the left side: The input data in blue and error output in pink are juxtaposed for effective comparison of the convergence rate. On the right side: The last 10 ms time band (area highlighted by black rectangle on the right side) for the error output from each technique is shown.
Figure 5. Comparison plots of generated synthetic data 1 with random noise processed by multi-channel noise cancellation techniques: (a) normalized LMS; (b) sigmoid variable step size LMS (SVSLMS); (c) regularized sigmoid variable step size LMS (RVSSLMS); (d) regularized sigmoid variable step size LMS (RSVSLMS). On the left side: The input data in blue and error output in pink are juxtaposed for effective comparison of the convergence rate. On the right side: The last 10 ms time band (area highlighted by black rectangle on the right side) for the error output from each technique is shown.
Energies 13 05873 g005
Figure 6. Comparison plots of generated synthetic signals with random sequence, Gaussian noise and impulse noise and median filter processed signals: (a) frequency spectrum of displayed signals; (b) The main signal with added impulsive noise at a time of 2.25 s; (c) filtered main signal; (d) the combined reference signal with added noise at time 3.75 s and (e) filtered reference signal.
Figure 6. Comparison plots of generated synthetic signals with random sequence, Gaussian noise and impulse noise and median filter processed signals: (a) frequency spectrum of displayed signals; (b) The main signal with added impulsive noise at a time of 2.25 s; (c) filtered main signal; (d) the combined reference signal with added noise at time 3.75 s and (e) filtered reference signal.
Energies 13 05873 g006
Figure 7. Comparison of retrieved signal from generated synthetic data with random sequence, Gaussian noise and impulse noise, using multi-channel noise cancellation techniques: (a) input FSK signal (b) NLMS; (c) sigmoid variable step-size LMS (SVSLMS); (d) regularized variable step-size LMS (RVSSLMS); (e) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 2 to 2.5 s is extracted from the first 5 s time band data.
Figure 7. Comparison of retrieved signal from generated synthetic data with random sequence, Gaussian noise and impulse noise, using multi-channel noise cancellation techniques: (a) input FSK signal (b) NLMS; (c) sigmoid variable step-size LMS (SVSLMS); (d) regularized variable step-size LMS (RVSSLMS); (e) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 2 to 2.5 s is extracted from the first 5 s time band data.
Energies 13 05873 g007
Figure 8. Comparison of retrieved signal from generated synthetic data with random sequence, Gaussian noise and impulse noise, using regularized sigmoid variable step-size LMS (RSVSLMS) technique. (a) The signal recorded over a time frame of 2 to 2.5 s and (b) the signal recorded over a time frame of 3.5 to 4 s. Each time frame contains a plot of the transmitted signal, recovered signal form data without median filter and recovered signal after median filter was applied to the data.
Figure 8. Comparison of retrieved signal from generated synthetic data with random sequence, Gaussian noise and impulse noise, using regularized sigmoid variable step-size LMS (RSVSLMS) technique. (a) The signal recorded over a time frame of 2 to 2.5 s and (b) the signal recorded over a time frame of 3.5 to 4 s. Each time frame contains a plot of the transmitted signal, recovered signal form data without median filter and recovered signal after median filter was applied to the data.
Energies 13 05873 g008
Figure 9. Plots of signal modulator waveform (a), binary phase shift key (BPSK) generated signal for data transmission (b), generated input signal representing recorded main data (c) and generated weighted reference signal representing the noise characteristics (d).
Figure 9. Plots of signal modulator waveform (a), binary phase shift key (BPSK) generated signal for data transmission (b), generated input signal representing recorded main data (c) and generated weighted reference signal representing the noise characteristics (d).
Energies 13 05873 g009
Figure 10. Comparison of retrieved signal from random sequence and Gaussian noise using multi-channel noise cancellation techniques: (a) LMS; (b) sigmoid variable step-size LMS (SVSLMS); (c) regularized variable step-size LMS (RVSSLMS); (d) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 40 to 44 s is extracted from the first 100 s time band data.
Figure 10. Comparison of retrieved signal from random sequence and Gaussian noise using multi-channel noise cancellation techniques: (a) LMS; (b) sigmoid variable step-size LMS (SVSLMS); (c) regularized variable step-size LMS (RVSSLMS); (d) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 40 to 44 s is extracted from the first 100 s time band data.
Energies 13 05873 g010
Figure 11. Comparison of MSE for existing variable step-size LMS and traditional LMS algorithms with new regularized sigmoid variable step-size LMS techniques.
Figure 11. Comparison of MSE for existing variable step-size LMS and traditional LMS algorithms with new regularized sigmoid variable step-size LMS techniques.
Energies 13 05873 g011
Figure 12. Comparison of retrieved signal from random sequence and Gaussian noise generated at different center frequencies, using multi-channel noise cancellation techniques: (a) LMS; (b) sigmoid variable step-size LMS (SVSLMS); (c) regularized variable step-size LMS (RVSSLMS); (d) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 40 to 44 s is extracted from the first 100 s time band data.
Figure 12. Comparison of retrieved signal from random sequence and Gaussian noise generated at different center frequencies, using multi-channel noise cancellation techniques: (a) LMS; (b) sigmoid variable step-size LMS (SVSLMS); (c) regularized variable step-size LMS (RVSSLMS); (d) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 40 to 44 s is extracted from the first 100 s time band data.
Energies 13 05873 g012
Figure 13. Comparison of MSE for existing variable step-size LMS and traditional LMS algorithms with new regularized sigmoid variable step-size LMS techniques.
Figure 13. Comparison of MSE for existing variable step-size LMS and traditional LMS algorithms with new regularized sigmoid variable step-size LMS techniques.
Energies 13 05873 g013
Figure 14. Data acquisition array for multi-channel adaptive noise canceling in EM telemetry.
Figure 14. Data acquisition array for multi-channel adaptive noise canceling in EM telemetry.
Energies 13 05873 g014
Figure 15. Electromagnetic telemetry noise characteristics and generated synthetic signal: (a) frequency spectrum of measured noise data, (b) synthetic EM telemetry data and (c) normalized amplitude spectrum of the main and the reference signals. The synthetic EMT data are used as the main signal (input data) in the multi-channel adaptive noise canceling system.
Figure 15. Electromagnetic telemetry noise characteristics and generated synthetic signal: (a) frequency spectrum of measured noise data, (b) synthetic EM telemetry data and (c) normalized amplitude spectrum of the main and the reference signals. The synthetic EMT data are used as the main signal (input data) in the multi-channel adaptive noise canceling system.
Energies 13 05873 g015
Figure 16. Comparison of pseudo-synthetic telemetry processing results for time range: (a) 133 to 139 s, and (b) 336 to 344 s. The results were obtained by using multi-channel noise cancellation techniques. Top: Sigmoid variable step-size LMS (SVSLMS). Middle: Regularized variable step-size LMS (RVSSLMS). Below: Regularized sigmoid variable step-size LMS (RSVSLMS).
Figure 16. Comparison of pseudo-synthetic telemetry processing results for time range: (a) 133 to 139 s, and (b) 336 to 344 s. The results were obtained by using multi-channel noise cancellation techniques. Top: Sigmoid variable step-size LMS (SVSLMS). Middle: Regularized variable step-size LMS (RVSSLMS). Below: Regularized sigmoid variable step-size LMS (RSVSLMS).
Energies 13 05873 g016
Figure 17. Comparison of pseudo-synthetic telemetry processing results for time range 336 to 346 s. The main signal SNR was varied between −55 dB and 69.4 dB. The results were obtained by using multi-channel noise cancellation techniques: (a) regularized variable step-size LMS (RVSSLMS) technique; and (b) regularized sigmoid variable step-size LMS (RSVSLMS) technique.
Figure 17. Comparison of pseudo-synthetic telemetry processing results for time range 336 to 346 s. The main signal SNR was varied between −55 dB and 69.4 dB. The results were obtained by using multi-channel noise cancellation techniques: (a) regularized variable step-size LMS (RVSSLMS) technique; and (b) regularized sigmoid variable step-size LMS (RSVSLMS) technique.
Energies 13 05873 g017
Figure 18. Bar chart showing a plot of the output signal power for the different LMS techniques and the SNR of the main signal. The signal to noise ratios for the three datasets are 55.5 dB, 61.5 dB and 69.4 dB for data 1, 2 and 3, respectively. The RSVSLMS technique returns the retrieved signals with the highest power, and its followed closely by those obtained from the RVSSLMS technique.
Figure 18. Bar chart showing a plot of the output signal power for the different LMS techniques and the SNR of the main signal. The signal to noise ratios for the three datasets are 55.5 dB, 61.5 dB and 69.4 dB for data 1, 2 and 3, respectively. The RSVSLMS technique returns the retrieved signals with the highest power, and its followed closely by those obtained from the RVSSLMS technique.
Energies 13 05873 g018
Figure 19. Mean squared error plot for results obtained from the variable step LMS techniques when the main signal has an SNR of −55.5 dB. The regularized LMS techniques converged faster than the conventional. This quality makes the techniques more suitable for EMT data processing. The performance of the different techniques was assessed by considering the characteristics of the MSE for the early time and the whole periods.
Figure 19. Mean squared error plot for results obtained from the variable step LMS techniques when the main signal has an SNR of −55.5 dB. The regularized LMS techniques converged faster than the conventional. This quality makes the techniques more suitable for EMT data processing. The performance of the different techniques was assessed by considering the characteristics of the MSE for the early time and the whole periods.
Energies 13 05873 g019
Figure 20. Adopted EM telemetry data processing workflow and corresponding data processing result from each step.
Figure 20. Adopted EM telemetry data processing workflow and corresponding data processing result from each step.
Energies 13 05873 g020
Figure 21. Comparison of retrieved EMT signal using multi-channel noise cancellation techniques: (a) regularized variable step-size LMS (RVSSLMS); and (b) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 50 to 60 s is extracted from the first 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Figure 21. Comparison of retrieved EMT signal using multi-channel noise cancellation techniques: (a) regularized variable step-size LMS (RVSSLMS); and (b) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 50 to 60 s is extracted from the first 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Energies 13 05873 g021
Figure 22. Comparison of retrieved EMT signal using multi-channel noise cancellation techniques: (a) regularized variable step-size LMS (RVSSLMS); and (b) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 80 to 90 s is extracted from the first 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Figure 22. Comparison of retrieved EMT signal using multi-channel noise cancellation techniques: (a) regularized variable step-size LMS (RVSSLMS); and (b) regularized sigmoid variable step-size LMS (RSVSLMS). The displayed result with a time frame of 80 to 90 s is extracted from the first 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Energies 13 05873 g022
Figure 23. Comparison of retrieved EMT signal using different number of reference channels: (a) multiple reference signals were used, and (b) single reference signal was used. The displayed result with a time frame of 132 to 138 s is extracted from the other 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Figure 23. Comparison of retrieved EMT signal using different number of reference channels: (a) multiple reference signals were used, and (b) single reference signal was used. The displayed result with a time frame of 132 to 138 s is extracted from the other 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Energies 13 05873 g023
Figure 24. Comparison of retrieved EMT signal using different reference dipole arrays: (a) dipoles L2, L3 and L8 were used as reference signals; and (b) dipoles L4, L5 and L7 were used as the reference signal. The displayed result with a time frame of 43 to 50 s is extracted from the first 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Figure 24. Comparison of retrieved EMT signal using different reference dipole arrays: (a) dipoles L2, L3 and L8 were used as reference signals; and (b) dipoles L4, L5 and L7 were used as the reference signal. The displayed result with a time frame of 43 to 50 s is extracted from the first 100 s time band data. The images are the processing result after combining the adaptive noise canceling output with 10 Hz sine wave signal.
Energies 13 05873 g024
Figure 25. Comparison of mean squared deviation for the NLMS technique, existing variable step-size LMS and new regularized sigmoid variable step-size LMS techniques. The data processing was carried out in a single set with the time range limited to 200 s. The plot shows a possible similarity in the data processing quality between the two regularized processing techniques, with the RSVSLMS technique having a better consistency/continuity in the MSD value after convergence.
Figure 25. Comparison of mean squared deviation for the NLMS technique, existing variable step-size LMS and new regularized sigmoid variable step-size LMS techniques. The data processing was carried out in a single set with the time range limited to 200 s. The plot shows a possible similarity in the data processing quality between the two regularized processing techniques, with the RSVSLMS technique having a better consistency/continuity in the MSD value after convergence.
Energies 13 05873 g025
Table 1. Adaptive noise-canceling technique parameters.
Table 1. Adaptive noise-canceling technique parameters.
TechniqueβαλσμConvergence Time (*0.05 s)
NLMSNANANA0.010.001550
SVSLMS0.000015NANANA500
RVSSLMS0.0000150.5NANA120
RSVSLMS0.0000150.5NANA80
Table 2. Adaptive noise-canceling technique parameters.
Table 2. Adaptive noise-canceling technique parameters.
Techniqueβαλμ
LMSNANANA0.00001
SVSLMS0.000055NANA
RVSSLMS0.0000550.5NA
RSVSLMS0.0000550.5NA
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fayemi, O.; Di, Q.; Zhen, Q.; Wang, Y.L. Adaptive Processing for EM Telemetry Signal Recovery: Field Data from Sichuan Province. Energies 2020, 13, 5873. https://doi.org/10.3390/en13225873

AMA Style

Fayemi O, Di Q, Zhen Q, Wang YL. Adaptive Processing for EM Telemetry Signal Recovery: Field Data from Sichuan Province. Energies. 2020; 13(22):5873. https://doi.org/10.3390/en13225873

Chicago/Turabian Style

Fayemi, Olalekan, Qingyun Di, Qihui Zhen, and Yu L. Wang. 2020. "Adaptive Processing for EM Telemetry Signal Recovery: Field Data from Sichuan Province" Energies 13, no. 22: 5873. https://doi.org/10.3390/en13225873

APA Style

Fayemi, O., Di, Q., Zhen, Q., & Wang, Y. L. (2020). Adaptive Processing for EM Telemetry Signal Recovery: Field Data from Sichuan Province. Energies, 13(22), 5873. https://doi.org/10.3390/en13225873

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