1. Introduction
In the process of dynamic measurement while drilling (MWD), there are complex noises in the attitude measurement signal due to the strong vibration of bottom-hole drilling assemblies, which are produced by the bit–rock interaction and the collision between the drill string and the borehole wall. On the one hand, the frequency range of the noise presents the characteristics of wideness and randomness, which are shown in
Table 1; on the other hand, the lateral vibration acceleration of the near-bit during drilling is generally about 10 g or up to 30 g [
1], while the amplitude of the useful signal is generally less than 1 g. Therefore, the signal-to-noise ratio (SNR) of the attitude measurement signal is usually as low as −20 dB, or even lower. Thus, the useful signal is annihilated in the strong vibration background noise, which negatively interferes with the tool attitude measurement and makes the MWD invalid [
2,
3]. Multi-frequency, high amplitude noise interference and extremely low SNR are always the technical difficulties in the field of the dynamic measurement of steerable drilling tools.
At present, the commonly used detection method for underground weak SNR signals is a digital filter based on time-frequency analysis, the premise of which is that the spectrum of useful signal and noise does not overlap, the useful signal can be retained, and the irrelevant noise components can be filtered in the filtering process [
4,
5]. However, in the collected attitude measurement signals, due to the variety of interference sources, the frequency distribution of the noise signal is very complex, as shown in
Table 1 [
6], which is bound to be close to the frequency of a useful signal. Therefore, while suppressing the noise, the useful signal will inevitably be suppressed or damaged [
7]. Even because of the complex and changeable strong noise interference, it is difficult to detect the weak SNR signal which overlaps with the noise spectrum, which leads to the invalidity of MWD.
Because of the limitations of the traditional detection method of weak signals, it is necessary to explore a new detection method for weak SNR signals in the pit. In recent years, the intensive study of nonlinear science has provided a new idea to understand, analyze and solve the problem of weak signal detection. Chaos is a special form of nonlinear dynamic system. Since the 1990s, many scholars have applied chaotic oscillators to weak signal detection, and have achieved remarkable results [
8,
9,
10,
11]. Moreover, the emergence of memristor has made weak signal detection based on chaotic oscillator enter a new stage of research [
12]. As a kind of nanodevice with memory and nonlinear characteristics, memristor is conducive to the generation of a high-frequency chaotic oscillation signal with complex dynamic characteristics [
13]; as such, it has wider application prospects in image encryption [
14], secure communication [
15] and other fields. In addition, the nonlinear system will produce a stochastic resonance (SR) phenomenon under the appropriate parameter conditions [
16]. Scholars have deeply studied this phenomenon, and have applied it to underwater acoustic communication [
17], fault diagnosis [
18,
19] and Internet of Things communication [
20]. The experimental results show that this method can successfully detect the weak characteristic signal in strong noise or a complex noise background.
Through the above research, the SNR of the signal which was detected by chaos and SR methods is much lower than that of the traditional detection methods. As such, it has wider application prospects. In many nonlinear systems, the second-order Duffing system, which has both SR effects and chaotic phase transition (CPT) characteristics, has become a research hot spot in the field of weak signal detection.
SR refers to the phenomenon by which, under the synergistic effect of input signal and noise, the SNR of the output signal increases first and then decreases with the increase of noise by using a nonlinear system, and the peak value appears at a certain noise intensity, resulting in resonance output [
21,
22]. At present, there are few reports about the application of SR in weak periodic signal detection while drilling. The related research includes Zhang et al. [
23], who proposed an adaptive optimization method of SR parameters based on a genetic simulated annealing algorithm for mud pulse signal detection under the complex background of MWD, and realized signal detection under low SNR. Based on the noise energy collection mechanism of SR, Park et al. [
24] transformed the high-intensity and low-speed vibration noise generated by drilling tools in the process of drilling. The results show that this method can effectively identify the frequency of the signal to be measured. According to the Langevin equation, Chen [
25] established the SR detection model of a down-hole acoustic communication signal, selected appropriate parameters, and finally obtained the frequency value of the useful signal.
Through the above research, it is not difficult to see that using the SR effect produced by the nonlinear system under the condition of appropriate parameters can effectively improve the SNR of the system output signal, identify the frequency of the down-hole useful signal, to realize the detection of a weak signal. However, SR cannot measure the ‘size’ of the measured signal, that is, the amplitude and phase parameters of the signal, so the complete information of the down-hole attitude measurement signal cannot be obtained.
Using the sensitivity of the CPT to system parameter perturbation and its noise immunity, it can achieve extremely low SNR detection in the case of fewer measurement data points and arbitrary colored noise background [
26,
27]. Therefore, compared with the traditional detection method, this method can detect the signal with lower SNR, and is not limited by the statistical characteristics of the noise [
28,
29].
Wu et al. [
30] designed an improved Duffing chaotic system and applied it to the resistivity detection of a near-bit electromagnetic logging signal. The test results showed that the smaller the amplitude of the tested signal, the more sensitive and noise resistant the system is to the signal, and the lower the SNR of the detectable signal. At present, the research on the application of a second-order Duffing system to the weak signal detection of dynamic MWD is under-developed, and the current research results have certain limitations in their application: they are affected by the strong vibration near the bit, and the frequency of the downhole attitude measurement signal constantly fluctuates based on the rotation speed of the drilling tool, rather than being fixed. Therefore, it cannot be directly used as an input signal and linear superposition with the system driving term. Although the scale transformation method is used to solve this problem in [
31], the method requires the high accuracy of frequency matching, so it will reduce the accuracy of the frequency detection in a field-drilling application, and it is not easy to meet the real-time requirements of MWD.
In order to overcome the shortcomings and limitations of SR and chaos detection, a new detection scheme for weak SNR signals in dynamic MWD is investigated in this paper. The scheme takes the second-order Duffing system as the theoretical framework, and uses the synergy of SR and CPT to obtain the complete parameter information of the useful signal. It not only avoids the matching requirement of the frequency of the measured signal in chaos detection but also makes up for the deficiency that the signal “size” cannot be measured by SR. The simulation test was carried out on the vibration platform system and the field-drilling data of a well in northern Shaanxi, which verified the feasibility and efficiency of the proposed algorithm being applied to down-hole weak SNR signal detection.
2. Frequency Detection of a Weak Signal Based on the SR of a Variable-Scale Duffing System
2.1. Basic Principle of Frequency Detection Based on the SR of a Duffing System
Consider a second-order Duffing system driven by useful signal
s(
t) and noise signal
n(
t), as follows:
where
k is the damping ratio,
a and
b are the system parameters,
x is the system output of the Duffing oscillator, and
f(
x) = −
ax +
bx3 is the nonlinear restoring force.
When
s(
t) = 0 and
n(
t) = 0, the system potential function of the second-order Duffing system is as follows:
The relationship between the potential function of the system and the nonlinear restoring force is as follows:
The system potential function is shown as the blue curve in
Figure 1.
It can be seen that Equation (2) is a typical bistable potential field structure. Let
du(
x)/dx = 0; the equilibrium point of the system can be found as follows:
where
xm1 and
xm2 are the stable equilibrium points, and
xn is the unstable equilibrium point. Thus, two potential wells separated by the middle barrier are formed. The height difference between the barrier and the well is easy to find, as follows:
Under the above conditions, Brownian particles will eventually fall into one of the potential wells after a transient process. The specific potential well is determined by the initial conditions of the system. Therefore, the second-order Duffing system is also called a bistable system.
Let the useful signal be
s(
t) =
Acos(2π
f0t +
φ) and the noise signal be
n(
t) = 0; then, the system potential function is periodically modulated by the useful signal, and changes from
U(
x) to
V(
x), as follows:
In this way, the two potential wells of the potential function are periodically raised or deepened, as shown in the black curve in
Figure 1. At this time, there is a critical amplitude
Ac. When
A <
Ac, Brownian particles will cause local periodic motion near a potential well; when
A >
Ac, Brownian particles will cause large-scale periodic transition motion between two potential wells. The critical amplitude
Ac means that
V(
x) changes from a bistable structure to a monostable structure when
s(
t) is at its maximum or minimum, and its theoretical value is as follows:
When the useful signal
s(
t) and noise
n(
t) exist at the same time and the system parameters are appropriate, Brownian particles will cause a large-scale transition between the two potential wells with the help of noise, even if
A <
Ac. At this time, the signal, noise, and system achieve synergy, and noise will have a positive effect, such that part of the noise energy is transferred to the useful signal. The signal energy is enhanced, the output SNR of the system reaches its maximum, and the SR is realized, as shown in the red curve in
Figure 1.
When the second-order Duffing system produces the SR effect, there will be an obvious peak at the frequency of the signal s(t), which can be used to detect the frequency value of the weak SNR signal in the strong noise background. This is the basic principle of using the SR characteristics of the Duffing system to detect the frequency value of a weak signal.
Although the SR characteristics of the Duffing system can effectively identify the frequency value of a weak signal with noise interference, there are still limitations when it is applied to the detection of a down-hole attitude measurement signal: the SR of a Duffing system based on adiabatic approximation theory is only applicable to the condition of a small frequency parameter (near 0.01 Hz). However, in the process of MWD, the frequency of the useful signal is mainly determined by the rotation speed of the drill string, which varies from 1 to 3 Hz during normal drilling, which is beyond the limit of Duffing system SR characteristics on the frequency parameters. In this case, the SR phenomenon cannot be obtained by direct numerical calculation, and the useful signal in the noise cannot be identified.
2.2. Frequency Detection Based on the SR of a Variable-Scale Duffing System
In order to obtain the SR phenomenon through the Duffing system under the condition of large frequency parameters, this section uses the step transformation method to improve the standard Duffing model, which was shown in Equation (1). Furthermore, let the measured signal
s(
t) be as follows:
where
f0 is the actual frequency value of the measured signal.
Keeping the other parameters of the Duffing system unchanged, a variable step-size coefficient
R is introduced in order to enlarge
s(
t) by
R times on the time axis, namely
t′ =
Rt. In this case, if the variable step coefficient
R = 100
f, the measured signal can be expressed as
It can be seen from Equation (9) that the frequency of the measured signal is compressed from the actual value to 0.01 Hz on the frequency axis. Under the action of the Duffing system, the existence of signal s(t′) can be detected, and then the frequency value of the actual measured signal can be determined according to the relationship of R = 100f. However, when processing the actual measured signal, it is impossible to directly transform the collected sensor signal through a linear transformation, such as Equation (9). Therefore, the step size of the numerical calculation is transformed in order to realize frequency parameter reconstruction. The specific steps are as follows.
Step 1: Assuming that the sampling frequency is
fs, the step size of the numerical calculation should be as follows:
Step 2: By introducing the variable step-size coefficient
R, the numerical calculation step size becomes
At this time, the signal time interval is increased by R times, and the corresponding signal frequency is compressed by R times. In other words, a down-hole attitude measurement signal with a sampling frequency of fS and an actual frequency of f is transformed into a signal with the sampling frequency of fS/R and a characteristic frequency of f/R through step transformation.
Step 3: Input the signal after the step transformation into the Duffing system, and make the calculation step meet Equation (11); complete the parameter reconstruction of the MWD signal, and then identify the frequency value of the weak SNR signal through SR, which was produced by the Duffing detective model in Equation (1).
It can be seen from the above analysis that the step size transformation does not change the value of the measured signal, but only changes the frequency or time scale of the signal, compresses or amplifies it on the time axis, and reorders the values, which will not affect the final calculation results.
It can be seen from the above variable-scale Duffing system detection principle that no matter what the frequency parameter of the measured signal is, it is reduced to 0.01 Hz. The frequency of the attitude measurement signal is compressed by step transform coefficient
R to meet the parameter condition of SR frequency detection. The process of frequency detection based on the SR phenomenon of the variable-scale Duffing system for the attitude signal in MWD is shown in
Figure 2.
3. Parameter Estimation of a Weak Signal Based on the CPT of a Duffing System
In order to calculate the real-time attitude information of drilling tools, it is necessary to restore the complete signal waveform. Therefore, in addition to the frequency information extracted in the previous section, it is also necessary to obtain the amplitude and phase parameters of the useful signal. Another nonlinear property of the Duffing system, namely CPT, is used. However, the traditional chaos detection method based on the Duffing system cannot estimate multiple parameters of the signal at the same time, and the existing phase estimation method has low computational efficiency, which makes it difficult to meet the requirements of accuracy and real-time MWD. For this reason, this section proposes a method of amplitude and phase synchronization estimation that is suitable for down-hole attitude measurement signals.
For the following Holmes-type Duffing chaotic system:
where
Acos(
ωt +
α) is driving term;
γcos(
ωt +
φ) is input, that is, the measured signal
s(
t) in Equation (1); and
γ and
φ are the amplitude and phase of the measured signal, that is, the quantity to be calculated.
When the total amplitude of the signal at the right end of the equation changes, the state of the system goes through homoclinic orbit, chaos, and a large-scale periodic state, namely phase transition. Chaos detection obtains the amplitude and phase parameters of the measured signal through this phase transition.
In order to make the Duffing system produce the phase transition effect smoothly, the angular frequency of the driving term and input term should be set at 1 rad/s at the same time; otherwise, the dynamic characteristics of the system will deteriorate, which will bring great difficulties and errors to the detection. Therefore, based on the scale transformation theory mentioned in the previous section, the step size of the numerical calculation is transformed to reconstruct the angular frequency. Therefore, the Duffing chaotic system after the frequency transformation is obtained as follows:
By comparing Formulae (12) and (13), it can be seen that the angular frequency of the driving term and the input term simultaneously change from ω to 1, where the value of ω is provided by the SR frequency detection in the previous section.
By simplifying the first two terms at the right end of the equal sign in Equation (13), we can obtain
It can be seen that the input term and the driver term are combined into a trigonometric function term. By observing the position of
θ in Formula (14), we can see that the value of
θ hardly affects the amplitude parameters of the Duffing chaotic system, but only changes the initial position of the trajectory solution. It is the amplitude term of the trigonometric function that affects the detection result. Therefore, in order to turn the Duffing system from a chaotic state to a large-scale periodic state, the equation of the
M(
t) amplitude term is established as follows:
where 0.827 is the critical amplitude that can turn the system output from a chaotic state to a large-scale periodic state.
Let
Aa denote the amplitude of the driving term when the Duffing system changes from a chaotic state to a large-scale periodic state. As such, there will be
where
Aa can be obtained by observing the phase state of the Duffing system.
It can be seen from the equation above that when Aa is determined, (16) is a binary equation with the amplitude γ and phase φ of the measured signal as variables. Therefore, we consider the use of a similar method to establish a binary equation regarding γ and φ, and then estimate the amplitude and phase of the measured signal synchronously by solving the binary function equations.
For the parameter estimation problem to be solved in this section, first let the initial phase α of the driving signal take 0 and π respectively, and input the measured signal into two Duffing chaotic systems with different α; then, adjust the amplitude of the driving signal of the two Duffing equations step by step, observe the change of the phase diagram of the system, and record the amplitude when the phase transition occurs. If
Aa1 and
Aa2 are used to represent the amplitude of the driving signal of which the output phase state changes when
α is 0 and π respectively, then the binary quadratic equations regarding the amplitude
γ and phase
φ of the measured signal can be expressed as
Formula (19) is the estimation formula of the amplitude
γ and phase
φ of the measured signal. The principle of parameter synchronization estimation based on a double Duffing system is shown in
Figure 3.
Based on the above analysis, the steps of the frequency detection and parameter estimation can be summarized as follows:
Determine the parameters of the SR model of the Duffing system. For example, in Equation (1), k = 0.5, a = b = 1, and the initial value of state (x(0), x′(0)) = (0, 0), and the critical amplitude λC = (4a3/27b) ≈ 0.148 are usually taken.
Estimate the strength of the measured sensor signal, and normalize the collected signal; that is, the amplitude of the signal is linearly compressed or expanded so that the signal strength is within the appropriate processing range, i.e., less than the critical amplitude λC.
Introduce the variable-scale coefficient R to change the signal with the sampling frequency of fs and signal frequency of f0 into a signal with a sampling frequency of fs/R and an actual frequency of f0/R. Take the calculation step dt′ = R/fs to solve the Duffing SR model and obtain the system waveform and signal spectrum.
Adjust the variable-scale coefficient R and observe the spectrum of the output signal. If the signal has an obvious peak at frequency 0.01, it will produce SR effect, and the value of R is the frequency value of the measured signal.
Set the relevant parameters of the Duffing chaotic system, such as the amplitude of driving term A = 0.825 in Equation (12).
Set the initial phase of the driving term to 0 and π, and record it as detection systems 1 and 2. The measured signal is input into two detection systems respectively, and the change of the output phase diagram of the system is observed. If the output of detection system 1 is chaotic after adding the measurement signal, the amplitude of the driving signal will be increased gradually with a certain amplitude until the system output jumps to the periodic state; if the system output is periodic, the amplitude of the driving signal will be decreased gradually until the system jumps from the periodic state to the chaotic state. The amplitude of the system jumping into a large-scale periodic state is denoted as Aa1.
Use the same method to obtain the amplitude of detection system 2 when it jumps to a large-scale periodic state, which is recorded as Aa2. Finally, Aa1 and Aa2 are substituted into Equation (14) in order to obtain the estimated values of the amplitude and phase parameters of the MWD signal.
5. Conclusions
In the process of oil and gas well drilling, there are multi-frequency and high amplitude interference signals in attitude measurement due to severe vibration and rapid rotation. In order to solve this problem, this paper proposes a method of the frequency detection and parameter estimation of weak SNR signals while drilling based on the SR and CPT cooperative of a variable-scale Duffing system. The conclusions are summarized as follows:
(1) The proposed SR detection algorithm based on a variable-scale Duffing system can complete the parameter reconstruction of the frequency value of the measured signal without changing the discrete value, and can meet the small frequency parameter requirements of SR detection of the Duffing system. The experimental results show that the method is feasible and effective.
(2) Two Duffing oscillators with a different initial phase of driving signal are combined in order to estimate the amplitude and phase parameters of the measured signal. The experimental results show that the method can detect the amplitude and phase of the signal synchronously in severe background noise, and the estimation accuracy is high, which can effectively improve the accuracy of the attitude angle of the drilling tool.
(3) The proposed detection algorithm is used to identify the stimulation signal with an amplitude of 0.05, and the results show that the SNR of the method is as low as −19 dB; the experimental results show that with the decrease of the amplitude of the signal to be detected, the SNR of the signal to be detected is further reduced. The algorithm proposed in this paper provides a solution to the problems of weak SNR signals in MWD.
In addition, the theoretical analysis and simulation analysis show that the proposed method also has some limitations in its application, which are mainly reflected in the following two aspects: (1) due to the influence of parameter setting (the influence of noise intensity in SR) and operation state judgment, there are still some shortcomings, such as the low detection efficiency and easy misjudgment of the state transition; (2) each detection needs to complete the construction of two models (SR and CPT), and it is difficult to deal with the tested object with multiple frequency signals, which affects the applicability of the proposed method.
Therefore, future research studies will be carried out with regard to the following aspects: (1) we should not be limited to the Duffing oscillator as a nonlinear system, and should consider other systems that can produce stochastic resonance or chaos, such as Lorenz system or Chua’s circuit, and identify the frequency of the measured signal through its chaotic synchronization phenomenon; (2) we should combine the nonlinear system with other methods (such as neural networks) to further improve the detection efficiency and applicability.