Next Article in Journal
Simulation on Buffet Response and Mitigation of Variant-Tailed Aircraft in Maneuver State
Previous Article in Journal
Analytical Study of Nonlinear Flexural Vibration of a Beam with Geometric, Material and Combined Nonlinearities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stick–Slip Suppression in Drill String Systems Using a Novel Adaptive Sliding Mode Control Approach

1
Department of Mechanical Engineering, Texas A&M University at Qatar, Doha 23874, Qatar
2
Laboratory of Energy Applications and Renewable Energy Efficiency (LAPER), Faculty of Sciences of Tunis, El Manar University, Tunis 1068, Tunisia
3
Department of Engineering Technology and Industrial Distribution, Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Vibration 2024, 7(2), 479-502; https://doi.org/10.3390/vibration7020026
Submission received: 22 March 2024 / Revised: 3 May 2024 / Accepted: 3 May 2024 / Published: 23 May 2024

Abstract

:
A novel control technique is presented in this paper, which is based on a first-order adaptive sliding mode that ensures convergence in a finite time without any prior information on the upper limits of the parametric uncertainties and/or external disturbances. Based on an exponent reaching law, this controller uses two dynamically adaptive control gains. Once the sliding mode is reached, the dynamic gains decrease in order to loosen the system’s constraints, which guarantees minimal control effort. The proof of convergence was demonstrated according to Lyapunov’s criterion. The proposed algorithm was applied to a drill string system to evaluate its performance because such systems present variable operating conditions caused by, for example, the type of rock. The effectiveness of the proposed controller was evaluated by conducting a comparative study that involved comparing it against a commonly used sliding mode controller, as well as other recent adaptive sliding mode control techniques. The different mathematical performance measures included energy consumption. The proposed algorithm had the best performance measures with the lowest energy consumption and it was able to significantly improve the functioning of the drill string system. The results indicated that the proposed controller had 20% less chattering than the classic SM controller. Finally, the proposed controller was the most robust to uncertainties in system parameters and external disturbances, thus demonstrating the auto-adjustable features of the controller.

1. Introduction

Sliding mode controllers (SMCs) are advanced control methods that are widely used in uncertain systems with external disturbances that satisfy a matching condition [1]. As a variable structure control method, SMCs utilize a discontinuous control law that enables the closed-loop system to converge in a finite time. SMCs have gained significant attention from researchers because of their robustness to parameter variations and external perturbations [2,3], which makes them an attractive scheme for controlling different kinds of systems [4]. It has been used in different research fields such as power electronics [5], power systems [6], electromechanical systems [3,7], biological economics [8], the space industry [9], robotics [10], and drilling systems [11,12]. Despite the successful use of SMCs, it has some drawbacks [2]. In fact, prior knowledge of the upper bounds of the uncertainties and external disturbances is necessary for designing a sliding mode (SM) controller. These upper bounds are poorly known; generally, their values are overestimated in order to compensate for their effects. However, overestimation of these bounds can lead to excessive gains. In contrast, frequent switching of the control gain can lead to the chattering phenomenon, especially for first-order SMCs because of the discontinuous switching law. This harmful phenomenon is undesirable because it causes energy losses and damage to the system’s components, and reduces its lifetime. Various methods have been studied to address the issue of chattering. Several artificial intelligence-based techniques have been applied to the SM technique [13,14], such as fuzzy rule-based algorithms that result in smooth dynamics when the system is near the switching region. Other techniques include the use of high-order SM controllers to reduce chattering. The authors of [15] proposed a second-order SMC for the control of robot manipulators that showed good tracking performance with minimal chattering. However, prior awareness of the upper limits of the uncertainties was necessary for implementing the proposed controller and it also needed information about the successive derivatives of the sliding surface. In [16], the authors implemented a super-twisting algorithm for controlling a four-rotor helicopter. The studied control algorithm managed to reduce the chattering phenomenon taking into account modelling imprecisions and external factors. The authors in [17] combined the fractional-order nonsingular SMC, the super-twisting SMC, and a double-loop fuzzy neural network algorithm for controlling a micro-gyroscope with unknown uncertainties. The proposed controller had the advantage of the fractional order controller, which improved the flexibility of the control and the super-twisting algorithm, and, thus, effectively solved the chattering problem. In [18], the authors presented a sliding mode control (SMC) method that utilizes a double-hidden-layer recurrent neural network with fuzzy logic to regulate a single-phase active power filter. The proposed algorithm guaranteed the control of the active power filter with great precision. However, the definition of the fuzzy rules was very complex, and the controller’s performance depended directly on the number of fuzzy rules implemented. An experimental investigation of a recurrent neural network fractional-order SMC applied on an active power filter was provided in [19]. The experimental outcomes confirmed the efficiency of the proposed control algorithm regarding the robustness and compensation performance compared with the standard neural SMC. However, the main disadvantage of the proposed algorithm was that it required much more data than traditional machine learning algorithms. With the exception of the super-twisting algorithm, higher-order sliding mode (SM) algorithms are only suitable for systems with a relative degree above one. One way to effectively reduce chattering and eliminate the need for a priori knowledge of uncertainty bounds is by defining dynamic SM controller gains with respect to the sliding surface.
Different adaptive SMC schemes have been proposed in the literature. To address the issue of chattering, the authors in [20] developed an adaptive SMC strategy for controlling a micro-AUV in the presence of water currents and parametric uncertainties. In [21], the authors designed an embedded flight controller based on the adaptive SM technique. Both of these previous studies investigated similar algorithms for adaptive gains based on index-reaching laws. Depending on the sliding surface, the dynamic gain was defined either by a constant that set the minimal value of the gain or by an expression containing the sign function of the sliding variable multiplied by the rate of adaptation. Both control strategies demonstrated the ability to effectively track the desired velocity within an acceptable time frame. Nevertheless, the proposed controllers lacked experimental validation to confirm their efficacy. The authors of [22] proposed an adaptive form of the controller proposed in [18]. The proposed adaptive sliding mode control demonstrated superior compensation performance and tracking accuracy when compared to the basic SMC. The authors introduced two first-order adaptive sliding mode controllers that employ constant reaching laws. The gain dynamics of the first controller included a low-pass filter to obtain the average of the sign function when the sliding region had been reached, and a constant reaching law to ensure convergence to the SM. If this filter is not properly adjusted, it may distort the signal, for example, through loss of information by attenuating magnitude. However, a phase shift is inevitable with a low-pass filter. To avoid the use of such a filter, the second control law is defined by a dynamic gain involving the sliding surface’s absolute value. The authors of [23] used an adaptive gain similar to [1]. However, the discontinuous control law included the sliding surface term without any functions, including the sign function. The authors of [24] proposed a second-order SM controller in which the gain gradually increased until the SM had been reached. However, the gain of the proposed algorithm never decreases, which can result in an overestimation of the gain and more chattering. The authors of [25] proposed a higher-order adaptive SM controller. This approach featured a unique dual-layer structure that relied on the application of equivalent injection principles and necessitated the existence of the first and second derivatives of the uncertainty. In [26], the authors investigated an adaptive SMC based on an exponential-reaching law. A simple adaptation scheme based on a single adaptive gain was used as an adaptive law for the controller gain. The controller showed great robustness to parameter uncertainties and external disturbances. Yet, a considerable amount of fluctuation remained at the bit velocity when applied with measurement noise.
It is important to note that inadequate estimation of either the uncertainties or the external disturbances results in more chattering and minimizes the robustness of the controller. For complex systems, such as drill strings, the uncertainties and external perturbations are significantly hard to predict since they are directly related to several factors such as the operating mode and condition. However, adapting the SM controller gains makes it possible to ensure a SM that is independent of these upper bounds. In this paper, a new adaptive SMC design is proposed to tackle the uncertainties and disturbances in a system that needs to be controlled, without prior knowledge. Various simulations were carried out to test its ability to reduce torsional vibrations while maintaining good tracking performance, and also to test its robustness regarding some parametric uncertainties. We also evaluated its performance under different operating conditions and external perturbations. The novel adaptation law proposed in this paper is distinguished by the utilization of two dynamically adapted gains that are inversely proportional to the sliding variables. The new adaptive scheme can make the system converge rapidly in a finite time and can reduce the chattering phenomenon in the traditional SMC. The proposed reaching law guarantees that the control constraints on the system will be minimal since both gains decrease once the SM has been reached. However, the proposed controller can be effectively implemented, ensuring good performance even in the presence of external perturbations and across diverse operating conditions.
This paper makes the following contributions to the literature:
  • Design of a novel adaptive SMC with first-order dynamics, without requiring advance information of the uncertainties and external perturbations’ upper bounds, thus significantly decreasing the chattering phenomenon while enhancing the controller component’s lifespan;
  • Application of the proposed controller to a drill string system under various operating regimes.
  • Comparative study of the first-order SM adaptive and proposed controllers. It includes a data analysis and a discussion to evaluate the performance of the proposed controller in comparison to the ones presented in the literature.
The structure of this paper can be summarized as follows. Section 2 provides the design strategy of the proposed controller along with a demonstration of its finite-time convergence via Lyapunov’s direct method. Section 3 describes the application of the proposed controller to a drill string system. In Section 4, the simulation results and an evaluation of the closed-loop system’s performance are discussed. Finally, the conclusions are presented in Section 5.

2. Control Strategy

This section provides a detailed description of the design process for the proposed adaptive sliding mode controller. The goal of this paper was to design an adaptive SMC with an interesting exponent-based convergence law. The novel scheme ensures dynamically adaptive gains, resulting in fast and precise convergence.

2.1. Problem Statement

Let us suppose a nonlinear system that includes external disturbances, which can be presented as follows:
x ˙ ( t ) = A x ( t ) + B u ( t ) + d ( t ) y ( t ) = C x ( t )
where A and B are input matrices, C is defined to have y ( t ) . u ( t ) is the control input, x ( t ) ( n ) is the state vector, and d ( t ) presents the bounded external perturbations satisfying | | d ( t ) | | < L where L > 0 is the upper limit of the disturbance. Note that the system mentioned in Equation (1) has a relative degree that is equal to one. The main target of the controller is to force the output of the system y ( t ) to follow a desired value y d ( t ) , which means that the error presenting the deviation from the desired trajectory can be written as follows:
e ( t ) = y ( t ) y d ( t )
If e ( t ) reaches a small value (i.e., e→0), the controller’s objective is achieved. The sliding surface term can be presented as follows:
s ( t ) = e ˙ ( t ) + λ e ( t )
Substitution of Equations (1) and (2) into Equation (3) results in:
s ( t ) = C A x ( t ) + C d ( t ) + λ ( y ( t ) y d ( t ) ) y ˙ d ( t ) + C B u ( t )
As previously stated, let us suppose that the sliding variable s has a relative degree of 1 with respect to the control input. The expression for the first derivative of the sliding surface described in Equation (4) can be expressed as follows:
s ˙ = ψ ( x ( t ) ) + Γ ( x ( t ) ) u ( t )
It is assumed that the functions ψ ( x ( t ) ) and γ ( x ( t ) ) , which belong to the set , are uncertain and can be expressed as follows:
ψ ( x ( t ) ) = ψ n o m ( x ( t ) ) + Δ ψ ( x ( t ) ) Γ ( x ( t ) ) = Γ n o m ( x ( t ) ) + Δ Γ ( x ( t ) )
where the nominal known functions are ψ n o m ( x ( t ) ) and Γ ( x ( t ) ) and Δ ψ ( x ( t ) ) and Δ Γ ( x ( t ) ) present the uncertainties terms. We also assume that ψ n o m ( x ( t ) ) and Γ ( x ( t ) ) are limited in the following manner:
ψ ( x ( t ) ) ψ M 0 < Γ m Γ ( x ( t ) ) Γ M
where ψ M , Γ m and Γ M are the bounds of ψ ( x ( t ) ) and Γ ( x ( t ) ) , respectively. These bounds are supposed to exist but are not known. According to [25], the sliding surfaces for a given system must be defined as the following set with ρ > 0 being the upper bound of s:
s = { x χ , | s ( x , t ) | < ρ }

2.2. Design of the Proposed Controller

The control law must be chosen so that the sliding surface, s, approaches zero in finite time. The control input u is expressed as follows:
u = u e q + u d i s c
The equivalent control u e q can be obtained by solving s ˙ = 0 , while u d i s c represents the controller’s discontinuous component. The adaptive SMC proposed in this study can be expressed as follows:
u d i s c = k 1 ^ ( t ) | s | 1 / 2 s i g n ( s ) k 2 ^ ( t ) s
where
k 2 ^ ( t ) = a k 1 ^ ( t )
The two equations defining the dynamic gains of the controller depend on the switching conditions given below:
k 1 ˙ ^ ( t ) = μ 1 | s | b s i g n ( | s | δ ) if k 1 ^ ( t ) > α μ 2 if k 1 ^ ( t ) α
The positive constants b, δ , μ 1 , and μ 2 , and the small positive constants α < < 1 and a < < 1 are used in the following equations. At the initial time, t = 0 , the value of k 1 ^ ( 0 ) should be greater than α to ensure that k 1 ^ ( t ) remains positive.
The dynamics of the gains ensure that the adaptive gains decrease once the SM has been reached. The control law described by the system of Equation (12) imposes few constraints on the controlled system. A general diagram of the proposed control law is shown in Figure 1.
Remark 1. 
In the work of [26], the authors investigated an adaptive SMC based on an exponential-reaching law. A simple adaptation scheme based on a single adaptive gain was used as an adaptive law for the controller gain. The controller showed great robustness to parameter uncertainties and external disturbances. Yet, a considerable amount of fluctuation remained at the bit velocity when applied with measurement noise. The proposed control schedule in this paper includes a real-time update of both gains ( k 1 , k 2 ) of u d i s c . Therefore, it proposes an adaptive switching control based on the exponent reaching law. The adaptive update of the controller gains allows the controller to maintain good performance even when there are external disturbances or uncertainties in the system’s parameters. This update also eliminates the need to reset the controller gains every time the operating conditions change. In this paper, the use of the absolute value of the sliding variable in the adaptive behavior, which was added to the classic scheme, improves the controller performance even for one operating condition without considering any perturbations/uncertainties. Compared to the work in [26], both adaptive gains in the convergence region increase quickly to a large value to counteract any disturbances. Then, the dynamic gains decrease until small values at the sliding region. Hence, it is possible to further reduce the problem related to the chattering phenomenon.

Proof of a General Convergence

This section employs Lyapunov’s direct method to prove the convergence of the closed-loop system.
Lemma 1. 
The proposed controller, given by Equation (11), and implemented on system (1) with dynamic gains defined by Equations (12) and (13), ensures that the gain k 1 ^ ( t ) is upper-bounded by a positive constant k 1 ¯ , as demonstrated below.
0 < k 1 ^ ( t ) k 1 ¯ , t > 0
Proof. 
If we take the first derivative of the sliding surface presented in Equation (5) and the controller described by Equations (11)–(13), the candidate Lyapunov function is chosen to be [26]:
V ( t ) = 1 2 s 2 + 1 2 γ ( k 1 ^ ( t ) k 1 ¯ ) 2
The derivative along the trajectory of the function V ( t ) is given by the following:
V ( t ) ˙ = s s ˙ + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 | s | b s i g n ( | s | δ )
Substituting the expression of s, as defined in Equations (6) and (15), can be rewritten as follows:
V ( t ) ˙ = s ψ ( x ( t ) ) + s Γ ( x ( t ) ) u ( t ) + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ )
Equation (11) can now be substituted into Equation (16), yielding the following:
V ( t ) ˙ = s ψ ( x ( t ) ) + s Γ ( x ( t ) ) ( k 1 ^ ( t ) s 1 / 2 s i g n ( s ) a k 1 ^ ( t ) s ) + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ )
V ( t ) ˙ = s ψ ( x ( t ) ) s Γ ( x ( t ) ) k 1 ^ ( t ) s 1 / 2 s i g n ( s ) a Γ k 1 ^ ( t ) s 2 + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ )
V ( t ) ˙ s ψ M s Γ m k 1 ^ ( t ) s 1 / 2 a Γ M k 1 ^ ( t ) s 2 + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ )
Adding and substituting the term Γ m k 1 ¯ ( | s | 1 / 2 a | s | ) gives the following:
V ( t ) ˙ s ψ M s Γ m k 1 ^ ( t ) ( s 1 / 2 a s ) + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ ) + s Γ M k 1 ¯ ( s 1 / 2 a s ) s Γ m k 1 ¯ ( s 1 / 2 a s )
After simplification, this gives the following:
V ( t ) ˙ s [ ψ M Γ m k 1 ¯ ( s 1 / 2 a s ) ] + ( k 1 ^ ( t ) k 1 ¯ ) [ s Γ m k 1 ^ ( t ) ( s 1 / 2 a s ) + 1 γ μ 1 s b s i g n ( s δ ) ]
Let us define a positive constant β . In this case, we have the following:
V ( t ) ˙ s [ ψ M Γ m k 1 ¯ ( s 1 / 2 a s ) ] + ( k 1 ^ ( t ) k 1 ¯ ) [ s Γ m k 1 ^ ( t ) ( s 1 / 2 a s ) + 1 γ μ 1 s b s i g n ( s δ ) ] + β k 1 ^ ( t ) k 1 ¯ β k 1 ^ ( t ) k 1 ¯
Therefore,
V ( t ) ˙ s ( ψ M + Γ m k 1 ¯ ( s 1 / 2 a s ) ) β k 1 ^ ( t ) k 1 ¯ k 1 ^ ( t ) k 1 ¯ s 3 / 2 ( Γ m k 1 ^ ( t ) ( 1 a s 1 / 2 ) ) k 1 ^ ( t ) k 1 ¯ ( 1 γ μ 1 s b s i g n ( s δ ) + β )
If the following two conditions are satisfied, then the derivative of the Lyapunov function V ( t ) is negative, as follows:
1 a | s | 1 / 2 > 0 1 γ μ 1 s i g n ( | s | δ ) + β > 0
The first assumption, a < 1 / | s | 1 / 2 , is a direct consequence of the selected value for the constant a. As for the second inequality in Equation (24), two scenarios can be examined:
i 
For the case | s | > δ , V ( t ) ˙ is negative-definite if μ 1 > γ β .
ii 
For the case | s | δ , V ( t ) ˙ is negative-definite if μ 1 > γ β
Thus, the condition can be written as follows:
μ 1 > s i g n ( γ β ) γ β
where μ 1 and μ 2 are the main parameters of the proposed controller. Tuning them is carried out as follows: To ensure that k 1 ^ ( t ) stays positive, it is important to choose a small value for μ 1 , such that the initial value of the proposed gain at t = 0 , k 1 ^ ( 0 ) , is greater than α . A small value of μ 2 must be chosen, as this represents the value of the dynamic gain when the system stabilizes at the sliding region. The selection of the control parameters was done through a trial-and-error process to achieve a balance between the rate of convergence and the magnitude of chattering in the system. Moreover, a small value of b must be chosen to reduce the instant overshoot caused by the gain dynamics. □

3. Drill String Application

A drill string system is composed of various components including a rotary table, drill pipes, drill collars, and a drill bit (see Figure 2). The rotary table generates the driving torque, which is then transmitted through the drill pipes to the drill bit for rock-cutting purposes [27]. To understand the structure of a drill string system, it can be divided into three parts: the surface, middle, and bottom [12]. The surface part includes the rotary table that is controlled by a regulator to maintain precise and robust drilling speed. The middle section contains the drill collars and drill pipes, which can span over several hundred meters to several kilometers. The bottom section consists of the drill bit that is tasked with overcoming the complicated and nonlinear bit-rock interaction, which can sometimes involve time delays.
The rotary motion generated by the rotary table is transmitted to the drill bit through the drill string. In addition to transmitting the rotary motion, drill strings also provide the necessary weight on bit (WOB), which aids in the drilling process [28]. The system’s interaction with the borehole results in various unwanted oscillations, including the stick–slip phenomenon, which is a severe form of torsional vibration that causes the drill string to stop rotating and then periodically rotate freely [29]. The main cause of this particular vibration is the highly nonlinear torque applied on the bit. This nonuniform rotation can reduce the drilling efficiency or even damage the drill string components [27]. Moreover, drill strings’ parameters vary depending on the operating modes and conditions [30]. For example, it is very challenging to know the correct stiffness and damping of the drill pipes. Moreover, the WOB varies, depending on the type of rock being cut. This always makes the drill string system susceptible to external perturbations [12]. The objective of this study is to evaluate and confirm the effectiveness of the controller in eliminating the stick–slip occurrence, regardless of any external disturbances or uncertainties in the system’s parameters.

3.1. Mathematical Model of the Drill String

There are various modeling techniques available to simulate the behavior of a drill string system numerically, such as those discussed in [31,32]. For pure torsional studies in particular, lumped parameter modeling has been widely used in the literature to assess stick–slip behavior [33,34]. In this study, a model with four degrees of freedom (4DOF) has been considered that includes the four subsystems described in Figure 2 [35]. The model used in this study is based on three main assumptions: (i) Only torsional vibrations are taken into account, (ii) the drill string is assumed to be vertical, and (iii) all drill pipes have the same inertia.
The dynamics of the torsional vibrations are described by the following equation [35]:
J d t φ ¨ d t + c d t φ ˙ d t + c p i ( φ ˙ d t φ ˙ p i ) + k p i ( φ d t φ p i ) = u J p i φ ¨ p i + c p i ( φ ˙ p i φ ˙ d t ) + c c o ( φ ˙ p i φ ˙ c o ) + k p i ( φ p i φ d t ) + k c o ( ψ p i ψ c o ) = 0 J c o φ ¨ c o + c c o ( φ ˙ c o φ p i ˙ ) + c b i ( φ ˙ c o φ ˙ b i ) + k c o ( φ c o φ p i ) + k b ( φ c o φ b i ) = 0 J b i φ ¨ b i + c b i t φ ˙ b i + c b i ( φ ˙ b i φ ˙ c o ) + k b i ( φ b i φ c o ) + T b i t = 0
where J d t , J p i , J c o , and J b i are the moments of inertia of the top drive, drill pipes, drill collars, and drill bit; φ d t , φ p i , φ c o , and φ b i are the angular positions of the top drive, drill pipes, drill collars, and drill bit; k p i , k c o , and k b i are the torsional stiffness of the drill pipes, drill collars, and drill bit; and c p i , c c o , and c b i represent the torsional damping of the drill pipes, drill collars, and drill bits. We can now define a state vector x ( t ) as follows:
x ( t ) = ( φ ˙ d t , φ d t φ p i , φ ˙ p i , φ p i φ c o , φ ˙ c o , φ c o φ b , φ ˙ b i ) T x ( t ) = ( x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 ) T
Using Equation (27), the system can be written in the following form [34]:
x 1 ˙ = 1 J d t [ ( c d t + c p i ) x 1 k p i x 2 + c p i x 3 + u ] x 2 ˙ = x 1 x 3 x 3 ˙ = 1 J p i [ c p i x 1 + k p i x 2 ( c p i + c c o x 3 ) k c o x 4 + c c o x 5 ] x 4 ˙ = x 3 x 5 x 5 ˙ = 1 J c o [ c c o x 3 + k c o x 4 ( c c o + c b i ) x 5 k b i x 6 + c b i x 7 ] x 6 ˙ = x 5 x 7 x 7 ˙ = 1 J b i [ c b i x 5 + k b i x 6 ( c b i + c b i t ) x 7 T b i t ( x ( t ) ) ]
The system presented in Equation (28) can be expressed in state space form as follows [34]:
x ˙ ( t ) = A ( x ( t ) ) + B u ( t ) + D T b i t ( x ( t ) ) y = C x ( t )
Here, the matrices A, B, C, and D have dimensions of 7 × 7 , 7 × 1 , 1 × 7 , and 7 × 1 , respectively. The matrices can be expressed as follows:
A = ( c d t + c p i ) J d t k p i J d t c p i J d t 0 0 0 0 1 0 1 0 0 0 0 c p i J p i k p i J p i ( c p i + c c o ) J p i k c o J p i c c o J p i 0 0 0 0 1 0 1 0 0 0 0 c c o J c o k c o J c o ( c c o + c b i ) J c o k b i J c o c b i J c o 0 0 0 0 1 0 1 0 0 0 0 c b i J b i k b i J b i ( c b i + c b i t ) J b i B = 1 J d t 0 0 0 0 0 0 D = 0 0 0 0 0 0 1 J b i C = 1 0 0 0 0 0
The torque on the bit must accurately reproduce the reaction of the borehole to the bit. The well-known models of Stribeck and Karnopp [36,37] are generally used for presenting bit–rock interactions. Both models have shown good results, but Karnopp’s model has an advantage in reproducing friction at the bit level due to its continuity at zero velocity, unlike Stribeck’s friction model. The Karnopp’s friction model can be expressed as follows:
T b i t = T if φ b i ˙ Δ w , T T a T a s i g n ( φ b i ˙ ) if φ b i ˙ Δ w , T > T a [ T 0 + ( T a T 0 ) e ξ φ b i ˙ Δ w ] s i g n ( φ b i ˙ ) if φ b i ˙ > Δ w
The maximum friction torque is denoted as T a = μ a W O B R , and the sliding friction torque is represented by T 0 = μ 0 W O B R . The parameter Δ w is used as a threshold value between static and sliding friction, while ξ denotes the rate of decline in friction torque.
In the field, drill strings are generally driven by motors with a nominal velocity ( v m ) that can reach 1500 RPM. However, the maximal velocity of the drill bit ( v b ) is generally limited to 200 RPM; n = v b / v m = 0.13 < < 1 . We can conclude that the actuator’s dynamics do not significantly interfere with the drill string’s dynamics. In many systems, disturbances and uncertainness can significantly reduce the controller’s performance. However, the SM technique has proven to be one of the most robust control techniques because it can estimate parameter variations and external disturbances that satisfy a matching condition (the disturbances are applied to the system through the same channels as the input signals, without affecting the order of the system). For the drill string system under study, the disturbances under investigation are those that act directly on the axial motion of the drill string (variation of the WOB) and directly affect the control input. It should be noted that in cases where disturbances cannot be calculated directly, they can be estimated based on the available measurements. However, a suitable control strategy needs to be developed either using these estimations or based on the measurements of the sensors to effectively counteract the perturbations.

3.2. Proof of Closed Loop Stability

The primary goal of the proposed control strategy is to regulate the bit velocity. To accomplish this, a sliding function is defined as per [37]:
s = ( x 1 x d ) + λ 0 t ( x 1 x d ) d t + λ 0 t ( x 1 x 7 ) d t
In the proposed control strategy, the desired drill string velocity is denoted by x d ( t ) . As the system reaches the sliding region s = 0 , it reaches the desired equilibrium point where the bit velocity approaches the desired velocity. The first-order derivative of the variable defined in the above equation can be given as follows:
s ˙ = x 1 ˙ + λ ( x 1 x d ) + λ ( x 1 x 7 )
Replacing x 1 ˙ in the previous equation results in the following equation:
s ˙ = ( 1 / J d t ) [ u ( c d t + c p i ) x 1 k p i x 2 + c p i x 3 ] + λ ( x 1 x d ) + λ ( x 1 x 7 )
In our case, Γ ( x ( t ) ) is a positive constant defined as follows:
Γ ( x ( t ) ) = 1 / J d t
For this particular case, Γ ( x ( t ) ) is restricted to a single constant value, and the function ψ ( x ( t ) ) can be expressed as follows:
ψ ( x ( t ) ) = 1 / J d t [ ( c d t + c p i ) x 1 k p i x 2 + c p i x 3 ] + λ ( x 1 x d ) + λ ( x 1 x 7 )
The values of c d t , c p i , and k p i are uncertain and can be expressed as follows:
c d t = c d t n o m + c d t c p i = c p i n o m + c p i k p i = k p i n o m + k p i
It is also assumed that c d t , c p i , and k p i are bounded as follows:
c d t c d t M c p i c p i M k p i k p i M
The parameters c d t M , c p i M , and k p i M represent the upper limits of c d t , c p i , and k p i , respectively. Although these upper limits exist, their values are unknown. As a result, we can derive the following inequality:
ψ ( x ( t ) ) ψ m a = 1 / J d t [ ( c d t M + c p i M ) x 1 k p i M x 2 + c p i M x 3 ] + λ ( x 1 x d ) + λ ( x 1 x 7 )
Since Γ ( x ( t ) ) is a positive constant equal 1 / J d t , Γ m can be replaced with 1 / J d t and this could not affect the convergence proof. Substituting Equation (37) into Equation (19) gives the following:
V ( t ) ˙ s ψ m a s ( 1 / J d t ) k 1 ^ ( t ) s 1 / 2 a ( 1 / J d t ) k 1 ^ ( t ) s 2 + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ )
Adding and substituting the term 1 / J d t k 1 ¯ ( | s | 1 / 2 a | s | ) gives the following:
V ( t ) ˙ s ψ m a s ( 1 / J d t ) k 1 ^ ( t ) ( s 1 / 2 a s ) + 1 γ ( k 1 ^ ( t ) k 1 ¯ ) μ 1 s b s i g n ( s δ ) + s ( 1 / J d t ) k 1 ¯ ( s 1 / 2 a s ) s ( 1 / J d t ) k 1 ¯ ( s 1 / 2 a s )
After simplification, this gives the following:
V ( t ) ˙ s [ ψ m a ( 1 / J d t ) k 1 ¯ ( s 1 / 2 a s ) ] + ( k 1 ^ ( t ) k 1 ¯ ) [ s ( 1 / J d t ) k 1 ^ ( t ) ( s 1 / 2 a s ) + 1 γ μ 1 s b s i g n ( s δ ) ]
Let us define a positive constant β ; in this case:
V ( t ) ˙ s [ ψ m a ( 1 / J d t ) k 1 ¯ ( s 1 / 2 a s ) ] + ( k 1 ^ ( t ) k 1 ¯ ) [ s ( 1 / J d t ) k 1 ^ ( t ) ( s 1 / 2 a s ) + 1 γ μ 1 s b s i g n ( s δ ) ] + β k 1 ^ ( t ) k 1 ¯ β k 1 ^ ( t ) k 1 ¯
Therefore,
V ( t ) ˙ s ( ψ m a + ( 1 / J d t ) k 1 ¯ ( s 1 / 2 a s ) ) β k 1 ^ ( t ) k 1 ¯ k 1 ^ ( t ) k 1 ¯ s 3 / 2 ( ( 1 / J d t ) k 1 ^ ( t ) ( 1 a s 1 / 2 ) ) k 1 ^ ( t ) k 1 ¯ ( 1 γ μ 1 s b s i g n ( s δ ) + β )
If we consider these two assumptions, then the time derivative of V ( t ) , denoted as V ˙ ( t ) , will be negative.
1 a | s | 1 / 2 > 0 1 γ μ 1 s i g n ( | s | δ ) + β > 0
Therefore, following the same steps as the general convergence section, the closed-loop convergence can be proved. The equivalent control is defined as the solution of s ˙ = 0 , which can be expressed as follows:
u e q = c p i ( x 1 x 3 ) + k p i x 2 + c d t x 1 J d t λ ( x 1 x d ) J d t λ ( x 1 x 7 )
In Equation (30), the sliding surface depends on the velocities of both the drill bit and the rotary table following a desired reference velocity. Generally, the drill string can be driven with various reference velocities depending on the length and diameter of the drill string and the type of rock of the borehole. These velocities are chosen in order to have an adequate rate of penetration, providing fast drilling while maintaining the integrity of the drilling components. All variables existing in the control command are supposed to be bounded as follows:
x 1 [ x 1 m i n , x 1 m a x ] R P M x 7 [ x 7 m i n , x 7 m a x ] R P M
The desired velocity x d is bounded as well. This implies that both u e q and u d i s c are limited signals, which ensures that the control input u is a bounded signal with the following:
u [ u m i n , u m a x ] k N . m

4. Validation and Results

To evaluate the effectiveness of the suggested control approach, simulations were carried out using a drill string system. The aim of the controller was to follow a desired reference velocity and to counteract external disturbances in various operational conditions and with varying parametric uncertainties. The proposed controller (Proposed) was compared with various other controllers, including classic first-order SMC (CSMC), modified first-order SMC (MSMC), and the well-known adaptive SMC (ASMC). The switching law of CSMC, as described in [34], is as follows:
μ s w 1 = k 1 s i g n ( s )
According to [21], the MSMC is described as an exponent reaching law as follows:
μ s w 2 = k 2 | s | 1 / 2 s i g n ( s ) k 3 s
The well-known adaptive controller proposed in [38] is described by the following switching law and gain dynamics:
μ d i s c = k 4 ^ ( t ) s i g n ( s )
where the dynamics of the controller gain are as follows:
k 4 ^ ˙ ( t ) = μ 11 | s | s i g n ( | s | δ 2 ) if k 4 ^ ( t ) > μ 22 μ 22 if k 4 ^ ( t ) μ 22
where k 1 is large enough to counteract the external perturbations, k 2 and k 3 are positive constants, and u 11 > 0 and u 22 > 0 are very small. Table 1 gives the different parameter values for each controller with optimal adjustments.
The values of the parameters were selected to achieve optimal performance while maintaining a balance between the overshoot and settling time of the controllers, based on a given reference velocity and WOB.

4.1. Performance Criteria

To assess the performance of the various controllers, several performance metrics were used. The first measure was the integral square error (ISE), which evaluates the energy associated with the error by integrating the square of the error over time, and it emphasizes large errors. The ISE can be expressed as follows:
I S E = e 2 d t
where:
e = V r e f V b i t
The integral absolute error (IAE) measures the integral of the absolute error over time and provides an idea of the cumulative errors. It is useful for evaluating the performance of a controller in terms of the overall deviation from the desired trajectory. The IAE can be defined as follows:
I A E = | e | d t
In order to evaluate the steady-state error, the root mean square error (RMSE) was employed by calculating the square root of the average of squared errors over a specified time period, the RMSE provides an effective measure of the discrepancy between desired and actual values, taking into account the entire sample set. This metric assumes a Gaussian distribution of errors and is widely adopted as a reliable method for quantifying the accuracy of a system’s response.
R M S E = e 2 n d t
where n is the number of samples over time. The energy consumption of each controller can be quantified by the root mean square deviation (RMSD), which is defined as follows:
R M S D = T 2 n d t
where T is the torque input of the drive table. This criterion is very important since a controller with minimal energy consumption is greatly desired in order to save money in terms of energy, maintenance, etc. Overshoot (O%) and settling time (ST) criteria were also used to evaluate the controllers. These are calculated as follows:
O % = p e a k v a l u e f i n a l v a l u e f i n a l v a l u e × 100 %
The variable “ST” signifies the time required to achieve and sustain stability within a defined range, typically set at 5% of the final value. This measure captures the duration it takes for the system’s response to settle and remain within an acceptable margin of variation.

4.2. Simulation Results

Multiple experiments were performed using the MATLAB/Simulink platform, employing a sampling frequency of f s = 1 kHz. These tests aimed to evaluate the resilience of each controller when subjected to diverse operating conditions, including measurement noise and parametric uncertainties related to the stiffness and damping of the drill pipes, as well as the damping of the drive table. The selected parameters listed in Table 2 were extracted from an actual drill string system documented in [19]. These parameters were chosen based on their relevance and applicability to the specific drill string system under investigation.
In the remainder of the paper, the results of the various tests mentioned above will be presented alongside an analysis of the results obtained via the evaluation criteria already defined (Equations (50)–(55)).

4.2.1. Robustness to Measurement Noise

The robustness of the measurement noise section was tested to approximate experimental tests as closely as possible. Various sources can be implemented to represent noise. They all share the ability to add random variation to the signal under study. In the simulation studies, random white noise was selected as the noise source. This type of noise has a stochastic nature and is characterized by having equal power across all frequencies, making it suitable for imitating measurement noise [26]. By studying the behavior of controllers under such conditions, we can gain valuable insights into their robustness and adaptability. The system under control was subjected to a controlled level of white noise with a noise amplitude of 10 4 dB. The bit velocities of the controllers were then analyzed, and the outcomes were recorded in Figure 3. As shown in Figure 4a,b, the proposed controller exhibited the least response to measurement noise. Quantitative measurements revealed that the proposed controller achieved an average chattering magnitude reduction of 20% compared to the CSMC, 14% compared to the MSMC, and 18% compared to the ASMC. These numbers highlight the substantial advantage of the proposed controller in mitigating chattering.
The proposed controller showed a remarkable reduction in chattering magnitude compared to the other controllers while consuming the least energy. This proves that the proposed controller was the least affected by measurement noise. In addition to its success in chattering reduction, the proposed adaptive methodology has also demonstrated superior performance in terms of energy consumption compared to a previously proposed adaptive sliding mode controller [26]. Figure 5 shows the RMSD values of the proposed controller versus CSMS, MSMC, and ASMC controllers. The developed controller in [26] exhibited an RMSD value of 13,016, indicating a higher level of energy consumption. However, with the introduction of the new proposed controller, an improvement has been achieved, resulting in a smaller RMSD value of 13,012. The reduction signifies a noteworthy enhancement in energy efficiency. The reduction is due to the | s | b term in the proposed gain dynamics, which presents the key parameter that helped decrease the chattering effect. This improvement indicates that the proposed adaptive methodology optimizes the energy consumption of the system, resulting in reduced waste and increased overall energy efficiency.

4.2.2. Robustness to Parametric Uncertainties

In the literature, different stochastic and deterministic approaches have been applied to drill strings to investigate the problem of uncertainties. The research in [39,40] presented some examples of stochastic models. The authors in [41] proposed a stochastic model for the bit-rock interaction torque. This was achieved by introducing uncertain parameters into the dry friction model and subsequently conducting a Monte Carlo simulation. Adaptive techniques for the deterministic approach were also explored, such as the direct adaptive controller proposed in [42], which proved to be robust against changes in friction characteristics and induced noise. It is worth mentioning that the convergence of the controller can be affected by an inappropriate setup of the adaptive algorithms. In the subsequent simulations, parametric uncertainties were introduced for the drill pipes’ stiffness and damping, as well as the drive table’s damping. According to Equation (35), the equivalent control incorporating parametric uncertainties can be expressed as follows:
u e q = Δ 1 c p i ( x 1 x 3 ) + Δ 2 k p i x 2 + Δ 3 c d t x 1 J d t λ ( x 1 x d ) J d t λ ( x 1 x 7 )
To evaluate the performance of various controllers under diverse operating conditions, an investigation was conducted to analyze the influence of different parametric uncertainties on the system’s performance. The purpose here is to gain insights into the controllers’ effectiveness and assess their robustness in handling varying parametric uncertainties. These uncertainties are categorized into three conditions: condition 1, condition 2, and condition 3. Condition 1 refers to small parametric uncertainties, where the variations in the parameters are relatively minor. Condition 2 corresponds to medium parametric uncertainties, where the variations in the parameters are of moderate magnitude. Finally, condition 3 represents large parametric uncertainties, characterized by significant variations in the parameters. By examining the system’s behavior under these distinct conditions, we aim to assess its robustness and performance across a range of parametric uncertainties, as shown in Table 3.
To clearly demonstrate the responses of each controller to different variations, separate figures were used. Figure 6 depicts that the CSMC remained robust for small parametric uncertainties. However, when subjected to medium parametric uncertainties, the system exhibited convergence but with a significant static error of 4 rad/s from the desired reference velocity. In addition, the system lost its robustness and experienced the stick–slip phenomenon for large parametric uncertainties. In Condition 1, the controller was robust with a smaller overshoot value, and no parametric adjustment was required. However, for Conditions 2 and 3, the responses became mediocre and asymptotically stable, respectively. The controller’s parameters need to be adjusted in these cases.
Figure 7 demonstrates that the system response diverges completely, independently from the magnitudes of the parametric uncertainties. Therefore, it can be concluded that the MSMC is not robust even to minimal parametric uncertainties. In order to ensure stability and acceptable performance in the presence of any parametric variation, parametric adjustment is necessary.
For all cases of parametric uncertainty, the ASMC system’s trajectory converged to the reference velocity (see Figure 8). It should be noted that the settling time of the system response increases considerably by nearly 20% for large parametric uncertainties.
Figure 9 shows that the proposed controller maintains its robustness regardless of variation in the parameters. This proves that the proposed controller does not need any parameter adjustments for any of the cases of parametric uncertainty, which demonstrates the fast adaptability of the gains.
Figure 10, Figure 11 and Figure 12 show that the proposed controller had the smallest integral square error whereas the MSMC had the greatest value, with a difference of 35%. This indicates that the proposed algorithm had the least amount of error energy, while the MSMC had the highest amount of errors. The proposed algorithm had a smaller IAE value than the MSMC controller by almost 20%, indicating a closer response to the reference velocity. The RMSD values were used to compare the energy consumption of the different controllers, both with and without added measurement noise. As can be seen, the maximum energy consumption was found for the ASMC, whereas the minimum was found for the proposed controller. This criterion is very important since it affects the lifetime of the component. The proposed controller had the smallest RMSE value, which proves that it had the least steady-state error. This shows the good precision of the regulator and results in good static performance, which is required for such an application.
Table 4 presents a summary of the results comparing four controllers based on several performance metrics. The metrics considered are EC (energy consumption), RPU (robustness to parametric uncertainty), RMN (robustness to measurement noise), and P (static precision). These metrics provide insights into the controllers’ energy efficiency, their ability to handle uncertainties in system parameters and measurement noise, and their resilience to variations in weight on bit and reference velocity. The table offers a comprehensive evaluation of the controllers, facilitating a thorough assessment of their performance in different operational aspects. The controller with the best performance for each category is indicated by the “High” or “Very High” rating, while the worst performance is indicated by the “Low” rating. Table 4 shows that the proposed controller performed better in terms of precision, convergence time, and overshoot compared with both classical and adaptive SMCs. Given its robustness to parametric uncertainty and measurement noise, the proposed controller guaranteed minimal chattering, resulting in a smooth input signal. It is worth mentioning that a previous article presented an adaptive controller that yielded favorable results [43]. However, it is important to note that this controller did not undergo tests specifically examining its performance in the presence of measurement noise.

5. Conclusions

In summary, this paper presents a novel adaptive sliding mode control (SMC) approach aimed at suppressing stick–slip vibrations in drilling string systems. A new switching adaptive SM controller based on an exponent reaching law is introduced, which dynamically updates controller gains in real time. Through comprehensive simulations and analysis, the effectiveness of the proposed control strategy in improving drilling string performance while minimizing stick–slip phenomena is demonstrated. The proposed controller offers several advantages over traditional control methods, including robustness to measurement noise, resilience to parametric uncertainties, and convergence in finite time without requiring prior knowledge of uncertainties or disturbances. By adjusting the controller gains in real-time and using dynamic adjustments, the adaptive SMC approach achieves superior performance in suppressing stick–slip vibrations compared to conventional SMC techniques. Furthermore, to help other researchers effectively address the stick–slip phenomenon in drill string systems, the following engineering suggestions are recommended:
  • Select appropriate modeling approach: choose the most suitable modeling approach based on the specific requirements of the drilling operation. Torsional models are particularly effective for addressing stick–slip vibrations, while coupled models are better suited for mitigating bit bounce or whirling phenomena.
  • Consider certain assumptions for drill string modeling: this includes considering stiffness and damping factors and the relationship between weight on bit (WOB) and torque on bit to refine the modeling approach.
  • Account for nonlinear effects in torque on bit calculations: account for nonlinear effects by choosing adequate models such as the Karnopp friction model. Understanding and accounting for these nonlinearities are crucial for developing robust control strategies.
  • Prioritize the integration of adaptive robust control strategies into drill string systems: such control mechanisms can enhance system robustness and performance by dynamically adjusting to varying operating conditions and disturbances.

Author Contributions

Conceptualization, M.G. and L.S.; methodology, F.Z. and L.S.; software, F.Z.; validation, F.Z.; formal analysis, F.Z. and L.S.; investigation, F.Z. and L.S.; resources, M.G.; data curation, M.G.; writing—original draft preparation, F.Z.; writing—review and editing, F.Z. and M.G.; visualization, M.G.; supervision, M.G.; project administration, M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Qatar National Research Fund grant number NPRP10-0101-170081.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

This publication was made possible by NPRP grant [NPRP10-0101-170081] from the Qatar National Research Fund (a member of Qatar Foundation).

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AcronymDefinition
ASMCAdaptive Sliding Mode Controller
CPDECoupled Partial Differential Equation
CSMCClassic Sliding Mode Controller
DOFDegrees Of Freedom
LQGLinear–Quadratic–Gaussian
ODEsOrdinary Differential Equations
PDProportional–Derivative
PIProportional–Integral
PIDProportional–Integral–Derivative
ROPRate Of Penetration
SMSliding Mode
SMCSliding Mode Controller
SMCmModified Sliding Mode Controller
WOBWeight On Bit
SymbolParameter
sSliding surface
, ψ , Γ Uncertain functions belonging to the set
uControl input
λ Positive constant
s ˙ The time derivative of the sliding surface
VThe candidate Lyapunov function
δ , u 1 , u 2 , aPositive constants
k 1 ¯ Adaptive gain upper bound
k 1 , k 2 , k 3 Positive gains
k 4 ^ Adaptive gain of ASMC
V r e f Required reference velocity
V b i t Bit velocity
Δ 2 , u 1 1 , u 2 2 Positive constants

References

  1. Taleb, M.; Plestan, F.; Bououlid, B. Higher order sliding mode control based on adaptive first order sliding mode controller. IFAC Proc. Vol. 2014, 47, 1380–1385. [Google Scholar] [CrossRef]
  2. Sira-Ramirez, H. On the dynamical sliding mode control of nonlinear systems. Int. J. Control 1993, 57, 1039–1061. [Google Scholar] [CrossRef]
  3. Shao, K.; Zheng, J.; Wang, H.; Xu, F.; Wang, X.; Liang, B. Recursive sliding mode control with adaptive disturbance observer for a linear motor positioner. Mech. Syst. Signal Process. 2021, 146, 107014. [Google Scholar] [CrossRef]
  4. Li, H.; Maghareh, A.; Montoya, H.; Uribe, J.W.C.; Dyke, S.J.; Xu, Z. Sliding mode control design for the benchmark problem in real-time hybrid simulation. Mech. Syst. Signal Process. 2021, 151, 107364. [Google Scholar] [CrossRef]
  5. Tan, S.C.; Lai, Y.; Tse, C.K.; Cheung, M.K. Adaptive feedforward and feedback control schemes for sliding mode controlled power converters. IEEE Trans. Power Electron. 2006, 21, 182–192. [Google Scholar]
  6. Liu, X.; Su, X.; Shi, P.; Shen, C.; Peng, Y. Event-triggered sliding mode control of nonlinear dynamic systems. Automatica 2020, 112, 108738. [Google Scholar] [CrossRef]
  7. Zou, M.; Fang, P.; Hou, Y.; Wang, Y.; Hou, D.; Peng, H. Synchronization analysis of two eccentric rotors with double-frequency excitation considering sliding mode control. Commun. Nonlinear Sci. Numer. Simul. 2021, 92, 105458. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Shi, Z. Sliding mode control for uncertain TS fuzzy singular biological economic system. IEEE Access 2019, 7, 14387–14395. [Google Scholar] [CrossRef]
  9. Tian, B.; Cui, J.; Lu, H.; Zong, Q. Reentry attitude control for RLV based on adaptive event-triggered sliding mode. IEEE Access 2019, 7, 68429–68435. [Google Scholar] [CrossRef]
  10. Ahmed, S.; Ahmed, A.; Mansoor, I.; Junejo, F.; Saeed, A. Output feedback adaptive fractional-order super-twisting sliding mode control of robotic manipulator. Iran. J. Sci. Technol. Trans. Electr. Eng. 2021, 45, 335–347. [Google Scholar] [CrossRef]
  11. Nobahar Sadeghi, A.; Arıkan, K.B.; Özbek, M.E.; Baranoğlu, B. Robust and adaptive control design of a drilling rig during the operating modes. Meas. Control 2019, 52, 702–719. [Google Scholar] [CrossRef]
  12. Sadeghimehr, R.; Nikoofard, A.; Khaki Sedigh, A. Predictive-based sliding mode control for mitigating torsional vibration of drill string in the presence of input delay and external disturbance. J. Vib. Control 2020, 27, 2432–2450. [Google Scholar] [CrossRef]
  13. Ahamed, N.U.; Yusof, Z.; Hamedon, Z.; Rabbi, M.F.; Sikandar, T.; Palaniappan, R.; Ali, M.A.; Rahman, S.M.; Sundaraj, K. Fuzzy logic controller design for intelligent drilling system. In Proceedings of the 2016 IEEE International Conference on Automatic Control and Intelligent Systems (I2CACIS), Selangor, Malaysia, 22 October 2016; pp. 208–213. [Google Scholar]
  14. Raviraj, V.; Sen, P.C. Comparative study of proportional-integral, sliding mode, and fuzzy logic controllers for power converters. IEEE Trans. Ind. Appl. 1997, 33, 518–524. [Google Scholar] [CrossRef]
  15. Capisani, L.M.; Ferrara, A.; Magnani, L. Design and experimental validation of a second-order sliding-mode motion controller for robot manipulators. Int. J. Control 2009, 82, 365–377. [Google Scholar] [CrossRef]
  16. Derafa, L.; Fridman, L.; Benallegue, A.; Ouldali, A. Super twisting control algorithm for the four rotors helicopter attitude tracking problem. In Proceedings of the 2010 11th International Workshop on Variable Structure Systems (VSS), Mexico City, Mexico, 26–28 June 2010; pp. 62–67. [Google Scholar]
  17. Fei, J.; Feng, Z. Fractional-order finite-time super-twisting sliding mode control of micro gyroscope based on double-loop fuzzy neural network. IEEE Trans. Syst. Man Cybern. Syst. 2020, 51, 7692–7706. [Google Scholar] [CrossRef]
  18. Fei, J.; Chen, Y. Fuzzy double hidden layer recurrent neural terminal sliding mode control of single-phase active power filter. IEEE Trans. Fuzzy Syst. 2020, 29, 3067–3081. [Google Scholar] [CrossRef]
  19. Fei, J.; Wang, H. Experimental investigation of recurrent neural network fractional-order sliding mode control of active power filter. IEEE Trans. Circuits Syst. II Express Briefs 2019, 67, 2522–2526. [Google Scholar] [CrossRef]
  20. Rodriguez, J.; Castañeda, H.; Gordillo, J. Design of an Adaptive Sliding Mode Control for a Micro-AUV Subject to Water Currents and Parametric Uncertainties. J. Mar. Sci. Eng. 2019, 7, 445. [Google Scholar] [CrossRef]
  21. Castañeda, H.; Gordillo, J.L. Embedded Flight Control Based on Adaptive Sliding Mode Strategy for a Quadrotor Micro Air Vehicle. Electronics 2019, 8, 793. [Google Scholar] [CrossRef]
  22. Fei, J.; Chen, Y. Dynamic terminal sliding-mode control for single-phase active power filter using new feedback recurrent neural network. IEEE Trans. Power Electron. 2020, 35, 9904–9922. [Google Scholar] [CrossRef]
  23. Huang, Y.J.; Kuo, T.C.; Chang, S.H. Adaptive sliding-mode control for nonlinearsystems with uncertain parameters. IEEE Trans. Syst. Man, Cybern. Part B (Cybern.) 2008, 38, 534–539. [Google Scholar] [CrossRef]
  24. Li, P.; Zheng, Z.Q. Robust adaptive second-order sliding-mode control with fast transient performance. IET Control Theory Appl. 2012, 6, 305–312. [Google Scholar] [CrossRef]
  25. Edwards, C.; Shtessel, Y.B. Adaptive continuous higher order sliding mode control. Automatica 2016, 65, 183–190. [Google Scholar] [CrossRef]
  26. Zribi, F.; Sidhom, L.; Gharib, M.; Refaat, S.S. New algorithm based active method to eliminate stick-slip vibrations in drill string systems. Syst. Sci. Control Eng. 2022, 10, 468–487. [Google Scholar] [CrossRef]
  27. Cheng, J.; Wu, M.; Wu, F.; Lu, C.; Chen, X.; Cao, W. Modeling and Control of Drill-String System with Stick-Slip Vibrations Using LPV Technique. IEEE Trans. Control. Syst. Technol. 2020, 29, 718–730. [Google Scholar] [CrossRef]
  28. Dao, N.H.; Menand, S.; Isbell, M. Mitigating and understanding stick-slip in unconventional wells. In Proceedings of the SPE/IADC International Drilling Conference and Exhibition, The Hague, The Netherlands, 5–7 March 2019; p. SPE-194117-MS. [Google Scholar]
  29. Yang, L.; Tian, J.; Liu, Q.; Dai, L.; Hu, Z.; Li, J. The multidirectional vibration and coupling dynamics of drill string and its influence on the wellbore trajectory. J. Mech. Sci. Technol. 2020, 34, 2681–2692. [Google Scholar] [CrossRef]
  30. Jansen, J.D.; Van Den Steen, L. Active damping of self-excited torsional vibrations in oil well drillstrings. J. Sound Vib. 1995, 179, 647–668. [Google Scholar] [CrossRef]
  31. Dunayevsky, V.; Abbassian, F. Application of stability approach to bit dynamics. SPE Drill. Complet. 1998, 13, 99–107. [Google Scholar]
  32. Navarro-López, E.M.; Suárez, R. Practical approach to modelling and controlling stick-slip oscillations in oilwell drillstrings. In Proceedings of the 2004 IEEE International Conference on Control Applications, Taipei, Taiwan, 2–4 September 2004; Volume 2, pp. 1454–1460. [Google Scholar]
  33. Puebla, H.; Alvarez-Ramirez, J. Suppression of stick-slip in drillstrings: A control approach based on modeling error compensation. J. Sound Vib. 2008, 310, 881–901. [Google Scholar] [CrossRef]
  34. Navarro-López, E.M.; Licéaga-Castro, E. Non-desired transitions and sliding-mode control of a multi-DOF mechanical system with stick-slip oscillations. Chaos Solitons Fractals 2009, 41, 2035–2044. [Google Scholar] [CrossRef]
  35. Abdulgalil, F.; Siguerdidjane, H. Backstepping design for controlling rotary drilling system. In Proceedings of the 2005 IEEE Conference on Control Applications, Toronto, ON, Canada, 28–31 August 2005; pp. 120–124. [Google Scholar]
  36. Fu, M.; Zhang, P.; Li, J.; Wu, Y. Observer and reference governor based control strategy to suppress stick-slip vibrations in oil well drill-string. J. Sound Vib. 2019, 457, 37–50. [Google Scholar] [CrossRef]
  37. Liu, Y. Control of a class of multibody underactuated mechanical systems with discontinuous friction using sliding-mode. Trans. Inst. Meas. Control 2018, 40, 514–527. [Google Scholar] [CrossRef]
  38. Plestan, F.; Shtessel, Y.; Bregeault, V.; Poznyak, A. New methodologies for adaptive sliding mode control. Int. J. Control 2010, 83, 1907–1919. [Google Scholar] [CrossRef]
  39. Lobo, D.; Ritto, T.; Castello, D. A novel stochastic process to model the variation of rock strength in bit-rock interaction for the analysis of drill-string vibration. Mech. Syst. Signal Process. 2020, 141, 106451. [Google Scholar] [CrossRef]
  40. Chen, J.; Liao, H.; Zhang, Y.; Liang, H.; Liu, C.; Qi, D. A torsional-axial vibration analysis of drill string endowed with kinematic coupling and stochastic approach. J. Pet. Sci. Eng. 2021, 198, 108157. [Google Scholar] [CrossRef]
  41. Trindade, M.A. Robust evaluation of stability regions of oil-well drilling systems with uncertain bit-rock nonlinear interaction. J. Sound Vib. 2020, 483, 115481. [Google Scholar] [CrossRef]
  42. Wasilewski, M.; Pisarski, D.; Konowrocki, R.; Bajer, C.I. A new efficient adaptive control of torsional vibrations induced by switched nonlinear disturbances. Int. J. Appl. Math. Comput. Sci. 2019, 29, 285–303. [Google Scholar] [CrossRef]
  43. Zribi, F.; Sidhom, L.; Krama, A.; Gharib, M. Enhancement of drill string system operations with adaptive robust controller and hardware in-the-loop validation. J. Sound Vib. 2023, 556, 117716. [Google Scholar] [CrossRef]
Figure 1. Overall architecture of the proposed control scheme.
Figure 1. Overall architecture of the proposed control scheme.
Vibration 07 00026 g001
Figure 2. Schematic for the main parts of a drill string.
Figure 2. Schematic for the main parts of a drill string.
Vibration 07 00026 g002
Figure 3. Effects of measurement noise on drill string bit velocities.
Figure 3. Effects of measurement noise on drill string bit velocities.
Vibration 07 00026 g003
Figure 4. Effects of measurement noise on input signals; (a) global scale; and (b) local scale.
Figure 4. Effects of measurement noise on input signals; (a) global scale; and (b) local scale.
Vibration 07 00026 g004
Figure 5. RMSD values of different controllers.
Figure 5. RMSD values of different controllers.
Vibration 07 00026 g005
Figure 6. System responses of the CSMC with parametric uncertainty.
Figure 6. System responses of the CSMC with parametric uncertainty.
Vibration 07 00026 g006
Figure 7. System responses of the MSMC with parametric uncertainty.
Figure 7. System responses of the MSMC with parametric uncertainty.
Vibration 07 00026 g007
Figure 8. System responses of the ASMC with parametric uncertainty.
Figure 8. System responses of the ASMC with parametric uncertainty.
Vibration 07 00026 g008
Figure 9. System responses of the proposed controller with parametric uncertainty.
Figure 9. System responses of the proposed controller with parametric uncertainty.
Vibration 07 00026 g009
Figure 10. Root mean square error (RMSE) of different controllers.
Figure 10. Root mean square error (RMSE) of different controllers.
Vibration 07 00026 g010
Figure 11. The ISE (integral square error) and IAE (integral absolute error) values for all controllers.
Figure 11. The ISE (integral square error) and IAE (integral absolute error) values for all controllers.
Vibration 07 00026 g011
Figure 12. RMSD values with and without measurement noise for all controllers.
Figure 12. RMSD values with and without measurement noise for all controllers.
Vibration 07 00026 g012
Table 1. Parameter adjustments for each controller.
Table 1. Parameter adjustments for each controller.
ControllerParametersValue
CSMC k 1 3
SMCm k 2 3
k 3 0.001
ASMC μ 11 0.1
μ 22 0.01
Proposed μ 1 3
μ 2 0.001
α 0.1
a0.01
b0.125
Table 2. Numerical values for the simulation parameters.
Table 2. Numerical values for the simulation parameters.
ParameterSymbolValueUnits
Moment of inertia of the drive table J d t 930kg·m2
Moment of inertia of the drill bit J b i 471.96kg·m2
Moment of inertia of the drill pipes J p i 2782.25kg·m2
Moment of inertia of the drill collar J c o 750kg·m2
Torsional damping between drive table and drill pipe C d t 425N·m·s/rad
Torsional damping between the drill pipe and BHA C b i t 50N·m·s/rad
Torsional damping between BHA and drill bit C c o 190N·m·s/rad
Torsional damping between BHA and drill pipe C p i 193.61N·m·s/rad
Torsional stiffness between drive table and drill pipe k p i 698.06N·m/rad
Torsional stiffness between drill collar and BHA k c o 1080N·m/rad
Torsional stiffness between BHA and drill bit k b i 907.48N·m/rad
Decline rate of friction torque ξ 0.5s/rad
Static and sliding friction threshold Δ w 0.001rad/s
Radius of the drill bitR0.15m
Coulomb friction coefficient μ a 0.8
Static friction coefficient μ 0 0.5
Table 3. Conditions of parametric uncertainties.
Table 3. Conditions of parametric uncertainties.
Parametric Uncertainties Δ 1 (%) Δ 2 (%) Δ 3 (%)
Nominal condition000
Condition 1403030
Condition 2605040
Condition 3757575
Table 4. Comparative analysis of controllers based on performance metrics.
Table 4. Comparative analysis of controllers based on performance metrics.
ControllerECRPURMNP
CSMCHighLowLowHigh
MSMCLowLowLowLow
ASMCLowHighHighHigh
Presented in [43]UnknownVery HighUnknownVery High
ProposedVery HighVery HighVery HighVery High
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zribi, F.; Sidhom, L.; Gharib, M. Stick–Slip Suppression in Drill String Systems Using a Novel Adaptive Sliding Mode Control Approach. Vibration 2024, 7, 479-502. https://doi.org/10.3390/vibration7020026

AMA Style

Zribi F, Sidhom L, Gharib M. Stick–Slip Suppression in Drill String Systems Using a Novel Adaptive Sliding Mode Control Approach. Vibration. 2024; 7(2):479-502. https://doi.org/10.3390/vibration7020026

Chicago/Turabian Style

Zribi, Fourat, Lilia Sidhom, and Mohamed Gharib. 2024. "Stick–Slip Suppression in Drill String Systems Using a Novel Adaptive Sliding Mode Control Approach" Vibration 7, no. 2: 479-502. https://doi.org/10.3390/vibration7020026

APA Style

Zribi, F., Sidhom, L., & Gharib, M. (2024). Stick–Slip Suppression in Drill String Systems Using a Novel Adaptive Sliding Mode Control Approach. Vibration, 7(2), 479-502. https://doi.org/10.3390/vibration7020026

Article Metrics

Back to TopTop