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

: 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 ﬁlter and performing the signal analysis in real-time, we implemented combined fast e ﬀ ective impulse noise removal techniques. The combination of median and mean ﬁlters was e ﬀ ective 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 ﬁeld application. Data processing workﬂow 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 e ﬀ ective in recovering transmitted EM telemetry signals. variable step-size LMS techniques.


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 where S d 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 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 S n for the specified dataset S d . 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 ) S k (n) = A r · cos(ωt+n k ) (5) and their respective filter weights are given as W i (n) = [W 1 (n), W 2 (n), W 3 (n), . . . . . . , W k (n)] (6) such that the average reference signal is given as Based on (2), the estimated error signal can be written as follows: 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 (9) Energies 2020, 13, 5873 5 of 32 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) (10) 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.
Energies 2020, 13, x FOR PEER REVIEW 5 of 31 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).

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: (11) that is, μ is constant. NLMS: where max 1 0 λ β < < , and α is a parameter controlling the shape of the sigmoid function, while β controls the range of the  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, S d , recorded with the second electrode positioned tens to hundreds of meters away from the borehole. Reference channels record signals S n1 , S n2 and S n3 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 w 1 (n), w 2 (n) and w 3 (n) are adjusted by adaptive filters whose purpose is to produce an output, S a (n), which is the sum of the corresponding outputs of each adaptive filter. If S a (n) is as similar as possible to the noise characteristics in S n (n) by the minimum mean square error, by subtraction of S a (n) from the primary channel, the noise can be canceled to produce the initially transmitted signal S i (n).

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) (11) that is, µ is constant. NLMS: where 0 < β < 1, and 1 < α < 10. β in (12) is a small positive constant used to avoid division by zero when the input vector S d is zero and α is the adaptation positive constant. SVSLMS: W i (n) = W i (n − 1)+ 2 β· min abs(e(n))· 1 − exp(−α e(n) 2 , 1 − exp(−α e(n) 2 e(n)S i (n) (13) Energies 2020, 13, 5873 6 of 32 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: 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: 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.

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 α Energies 2020, 13, 5873 7 of 32 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. 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 x 10 -7 Step size parameter VS Error value e(N) Step size parameter VS Error value e(N) Step size parameter VS Error value Step size parameter VS Error value  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.
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. 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.

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 builtin 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. Step size parameter VS Error value e(N) RSVSLMS algorithm 1 SVSLMS algorithm RVSSLMS algorithm RSVSLMS algorithm 2 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.
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.

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.   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. 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 Energies 2020, 13, 5873 10 of 32 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.  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 = Y 1 , . . . , Y j , . . . , Y n , 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 nt mean(Y) of the series is calculated.
The obtained mean M nt is compared with previous mean M nt−1 , and the lowest of both is selected as M nt .
The difference between the mean value (M nt ) and mid-time value (Y j ) is then calculated; Z m = M nt − 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 (16) 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. 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.  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 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.    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. 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.  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 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.   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-

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.
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 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.  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. Figures 16 and 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 (Figures 15b and 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). 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.  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). 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. 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 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. Figures 21-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. Figures 21 and 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.
Energies 2020, 13, x FOR PEER REVIEW 24 of 31 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 piecewiselinear 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 timecontinuous 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. Figures 21-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. Figures 21 and 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 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 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. 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. Figure 25. Comparison of mean squared deviation for the NLMS technique, existing variable stepsize 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.

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 steadystate 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  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.

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.

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 B
Energies 2020, 13, x FOR PEER REVIEW 29 of 31 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.