Next Article in Journal
A CMOS-Compatible Carrier-Injection Plasmonic Micro-Ring Modulator (CIPMRM) with Stable Performance as Temperature Varying around 60 K
Next Article in Special Issue
Seismic Performance Evaluation of Highway Bridges under Scour and Chloride Ion Corrosion
Previous Article in Journal
Energy Evolution Analysis of Coal Fracture Damage Process Based on Digital Image Processing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Recoil Control of Deepwater-Drilling Riser with Optimal Guaranteed Cost H Control

1
College of Automation and Electronic Engineering, Qingdao University of Science & Technology, Qingdao 266061, China
2
College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(8), 3945; https://doi.org/10.3390/app12083945
Submission received: 26 February 2022 / Revised: 9 April 2022 / Accepted: 11 April 2022 / Published: 13 April 2022
(This article belongs to the Special Issue Multi-Hazard Analysis of Structures under Wind and Seismic Loads)

Abstract

:
In deepwater-drilling engineering, it is necessary to disconnect the bottom equipment of the lower marine-riser package from the blowout preventer when encountering multi-hazard environmental factors. In order to reduce the impact of recoil on the drilling platform after the sudden disconnection of the riser, in this paper, an optimal guaranteed cost H recoil control problem is considered for the drilling riser. First, a three-element mass-damper-spring deepwater-drilling riser model subject to fluid discharge and heave motion of offshore platform is given. Then, an optimal guaranteed cost H controller (OGCHC) is designed to suppress the recoil response of the drilling riser, and the sufficient conditions for the asymptotic stability of the closed-loop system are derived. Third, it is found through simulation results that the designed OGCHC can reduce the recoil response effectively. In order to further analyze the advantages of the OGCHC, the performance indices of the riser without active-recoil control and with optimal control (OC) and OGCHC are compared. It is shown that the average response amplitudes of three mass blocks of the riser are almost the same, while the control cost by the OGCHC is less than that by the OC. Further, under the designed recoil control, no riser compression occurs, thereby ensuring the safety of the riser system.

1. Introduction

Typhoon, surge, wave, and earthquake generally put forward strict requirements for the safety of engineering structures [1,2,3]. For example, ocean current, wind, earthquake, and wave act on the deep-sea-drilling equipment for a long time, potentially leading to economic losses and pipe damage [4,5]; therefore, the drilling riser is a weak connecting part for deepwater oil and gas exploration between the submarine wellhead and the drilling ship [6,7,8]. The components of a deepwater-drilling riser system mainly include the drilling platform, tensioner, split joint, single riser, buoyancy joint, filling valve, lower marine-riser package (LMRP), blowout preventer (BOP), etc. In case of extreme harsh environments or dynamic-positioning failure, emergency disconnection of the LMRP and BOP is required to prevent damage to the drilling vessel [9,10].
During the disconnection operation, the deepwater-drilling riser will move upward rapidly under the tension of the tensioner and the elastic potential energy, resulting in recoil movements and endangering the safety of the drilling platform; therefore, many scholars have conducted recoil analysis on the riser, see [11,12,13,14,15,16,17], and the references therein. It is known that during riser emergency disconnection, the external load in the recoil process of the deepwater-drilling riser mainly comes from viscous resistance caused by discharging drilling fluid, and more and more research results concerning recoil response have been reported. In [18], the auxiliary program is used to simulate the mud column discharging process without the filling valve. Simulation results show that the auxiliary program and the RESTEN program achieve good agreement in mud kinematic applications. In [11], an one-dimensional finite volume model is proposed to calculate the fluid pressure and velocity. In [19], a slug fluid model is established for seawater injection by the refilled valve. Based on this model, the drilling mud discharge model with a whole fluid column model is developed, where seawater and mud are considered as a whole for force analysis [20]. In [21], the friction force calculation method is introduced for slug fluid model and whole fluid column model. In particular, based on the whole fluid column model and the action force of the fluid on the riser inner wall, the drilling fluid discharge is analyzed in [15], where the drill pipe is also considered. In [22], by comparing the existing drilling mud discharge models, a more accurate three-dimension computational fluid dynamics model is developed. This shows that, to reflect practical response of drilling riser, a more accurate model of the drilling mud discharge process is required.
In recent decades, much more research attention mainly focuses on riser recoil analysis, while there are relatively few studies with regard to the riser recoil control of the drilling riser. In [23], a recoil control model of the drilling riser is proposed first, and then the linear quadratic optimal control method is proposed to control the riser recoil displacement. It is found that, under the designed optimal controller, the recoil response of the drilling riser has been reduced effectively. Note that, during the design of the optimal controller for the drilling riser, the effects of the mud discharge on the riser are not considered specifically. In fact, the classic linear quadratic optimal controller design is generally based on the known exact dynamic model, and it is not always effective to cope with uncertain systems subject to external disturbances. In other words, there still exists room for further reduction in the recoil response of the riser, which motivates our study in this paper.
Due to the existence of external disturbances and parameter perturbations, the quadratic performance index of the system with an optimal controller may generally not be minimized. In this case, one of the feasible ways is to resort to the guaranteed cost optimal control, by which a minimum upper bound of the performance index of the system may be possible. As a result, the guaranteed cost optimal control H has attracted much attention and applied in many practical applications, see references [24,25,26,27,28,29,30], and therein. Inspired by the above studies, in this paper, an optimal guaranteed cost H control scheme is designed to control the recoil displacement of the riser subject to friction force of fluid discharge and heave motion of the offshore platform.
The main contributions of this study are given as follows:
(i)
An optimal guaranteed cost H recoil control scheme is proposed for the drilling riser system in the event of emergency disconnection. Based on this scheme, the minimum upper bound of the performance index of the riser subject to friction force of drilling fluid discharge as well as heave motion of platform can be achieved.
(ii)
The existence conditions of the optimal guaranteed cost H recoil controller are derived, and the design algorithm of the recoil controller of the riser is presented. The effectiveness and advantages of the proposed active-recoil control scheme are verified through simulation results.

2. Problem Formulation

In this section, the whole fluid column model [20,23] of fluid discharge is introduced to compute the friction force of the fluid discharge on the riser. In the two cases of constant and time-varying friction coefficients, the friction force of the fluid discharge is analyzed briefly. Then, a three-element mass-damper-spring recoil control system of the drilling riser system subject to platform heave motion and fluid discharge is given. Third, an optimal guaranteed cost H recoil control problem of the drilling riser is formulated.

2.1. Whole Fluid Column Model of Fluid Discharge and Heave Motion Model of Platform

It is known that the viscous resistance of fluid discharge and the heave motion of platform have obvious effects on the recoil response of the riser. In this paper, the whole fluid column model of fluid discharge is utilized to analyze the influence on the recoil response of the drilling riser. As stated in [20,23], the characteristics of the fluid discharge process can be described by
ρ w g A s ( L r L m ) + ρ m g A s L m ρ w g A s L r F m F w F d = [ ρ w A s ( L r L m ) + ρ m A s L m ] a
where
F m = f m A s ρ m v 2 L m 2 D s , F d = ρ m A s v 2 2 F w = f w A s ρ w v 2 ( L r L m ) 2 D s
where ρ w and ρ m are the densities of seawater and drilling mud, respectively, kg/m 3 ; L m is the length of drilling mud, m; F m is the friction force of drilling mud, N; f m is the friction coefficient of drilling mud; A s is the flow area of the riser, m 2 ; D s and L r are the hydraulic diameter and the length of the riser, respectively, m; F d is the frontal force, N; F w is the friction force of seawater, N; f w is the friction coefficient of seawater; v is the discharge velocity of fluid at time t, m/s; a is the acceleration of fluid at time t, m/s 2 .
The drilling platform is connected to the deepwater-drilling riser by the tensioners; the heave motion of the drilling platform can be approximated as the superposition of the sine functions:
w ( t ) = i = 1 n α i sin ( ω i t + φ i )
where α i is the random amplitude of heave motion, ω i is the random frequency of heave motion, φ i is the random phase of heave motion, and n is a given positive integer.

2.2. Dynamic Model of a Riser-Tension System

Notice that the tension system uses the hydraulic cylinder to drive the riser directly, and provides the required tension force for the riser. The tension force generated by a single tensioner can be calculated by [23]
N t ( t ) = f t k 1 x p u ( t )
where
f t = P H 0 A r P L 0 A p , k 1 = ϵ P H 0 A r 2 V H 0 + ϵ P L 0 A p 2 V L 0
with f t representing the top tension force of the initial disconnection moment, N; P H 0 and P L 0 are the pressure of the high pressure air vessel and low pressure nitrogen vessel at the initial time, respectively, Pa; A r and A p are the area of piston rod side and piston rod-less side, respectively, m 2 ; k 1 is the spring stiffness of tensioner, N/m; ϵ is the gas constant, ϵ [ 1.0 , 1.4 ] ; V H 0 and V L 0 are the volume of the high pressure gas and low pressure nitrogen vessel at the initial time, respectively, m 3 ; x p is the displacement of piston relative to hydraulic cylinder at time t, m; u ( t ) is the pressure drop.
Now, we aim to present a dynamic model of the riser-tension recoil control system. Similar to [15,23], the system can be simplified as a three-element mass-damper-spring model as seen in Figure 1, where m 1 , m 2 , and m 3 are the masses of three blocks of the whole riser, x 1 , x 2 , and x 3 are displacements of three mass blocks of the riser, c 1 , c 2 , and c 3 are the damping on the three mass blocks, k 2 and k 3 are equivalent spring stiffness of riser system with k 2 = k 3 = E A e / L . E is the elastic modulus of riser, A e is the cross-sectional area of riser, and L is the effective length of the riser.
To present dynamic equations of the riser, the following assumptions are required [15,23].
Assumption A1.
The upward direction of force and displacement is positive, and the coordinate origin of riser system is located at the position of respective initial disconnection moment.
Assumption A2.
The influence of buoyancy modules on the riser stiffness is ignored, the drilling pipe and the riser are supposed to be concentric. The viscous resistance of drilling mud discharge is evenly distributed along the axial direction of the riser.
Based on the Assumptions aforesaid, by Newton’s second law, the dynamic equation of riser-tension control system can be expressed as
m 1 x ¨ 1 = k 2 ( x 1 x 2 ) c 1 x ˙ 1 c 2 ( x ˙ 1 x ˙ 2 ) + N t F 1 m 2 x ¨ 2 = k 2 ( x 1 x 2 ) + c 2 ( x ˙ 1 x ˙ 2 ) c 3 ( x ˙ 2 x ˙ 3 ) k 3 ( x 2 + x 0 x 3 ) F 2 m 3 x ¨ 3 = k 3 ( x 2 + x 0 x 3 ) + c 3 ( x ˙ 2 x ˙ 3 ) F 3
where x 0 = ( f t m 1 g + F b 1 ) / k 2 , and
F 1 = F w 1 + F m 1 F b 1 + m 1 g F 2 = F w 2 + F m 2 F b 2 + m 2 g F 3 = F w 3 + F m 3 F b 3 + m 3 g
with F w i = F w / 3 , F m i = F m / 3 , i = 1 , 2 , 3 , and F b 1 , F b 2 , F b 3 representing the buoyancy during recoil of the riser.
The mass of the three blocks can be computed as
m 1 = M riser / 3 + M slip m 2 = M riser / 3 m 3 = M riser / 3 + M LMRP
where M riser = N b M buoyancy is the dry weight of the riser, kg; N b is the number of buoyancy joint; M buoyancy is the dry weight of the buoyancy joint, kg; M slip is the dry weight of slip joint outer cylinder, kg; M LMRP is the dry weight of LMRP, kg.
The damping of the three blocks can be calculated by
c 1 = 0.5 k 1 m 1 , c 2 = 2 ζ k 2 m 2 , c 3 = 2 ζ k 3 m 3
where the damping ratio ζ is 0.01.

2.3. State Space Model of the Riser System and Active-Recoil Control Problem

Let
z 1 ( t ) = x 1 ( t ) , z 2 ( t ) = x ˙ 1 ( t ) , z 3 ( t ) = x 2 ( t ) z 4 ( t ) = x ˙ 2 ( t ) , z 5 ( t ) = x 3 ( t ) , z 6 ( t ) = x ˙ 3 ( t )
and denote z ( t ) : = z 1 ( t ) z 2 ( t ) z 3 ( t ) z 4 ( t ) z 5 ( t ) z 6 ( t ) T . Then, the dynamic Equation (6) can be written as
z ˙ ( t ) = A z ( t ) + B u ( t ) + D 1 w ( t ) + D ˜ 2 h 1 ( t ) , z ( 0 ) = z 0
where z 0 is the system initial value as
z 0 = [ 0 0 x 0 0 x ı 0 ] T
with x ı = ( f t m 1 g m 2 g + F b 1 + F b 2 ) / k 3 , and
A = 0 1 0 0 0 0 k 1 + k 2 m 1 c 1 + c 2 m 1 k 2 m 1 c 2 m 1 0 0 0 0 0 1 0 0 k 2 m 2 c 2 m 2 k 2 + k 3 m 2 c 2 + c 3 m 2 k 3 m 2 c 3 m 2 0 0 0 0 0 1 0 0 k 3 m 3 c 3 m 3 k 3 m 3 c 3 m 3 , B = 0 1 m 1 0 0 0 0 D 1 = 0 k 1 m 1 0 0 0 0 , D ˜ 2 = 0 0 0 1 m 1 0 0 0 0 0 0 1 m 2 0 0 0 0 0 0 1 m 3 , h 1 ( t ) = F 1 f t F 2 + k 3 x 0 F 3 k 3 x 0
Note that the equilibrium state of the system is given by
z e 1 = κ k 1 z e 3 = κ k 1 m 2 g + m 3 g F b 2 F b 3 k 2 z e 5 = κ k 1 + κ k 2 m 3 g F b 3 k 3 z e 2 = z e 4 = z e 6 = 0
where κ = f t m 1 g m 2 g m 3 g + F b 1 + F b 2 + F b 3 .
Let
z ˜ i ( t ) = z i ( t ) z e i , i = 1 , 2 , , 6 .
Then from (10)–(14), one yields an incremental equation of the system as
z ˜ ˙ ( t ) = A z ˜ ( t ) + B u ( t ) + D 1 w ( t ) + D 2 h ( t ) , z ˜ ( 0 ) = z ˜ 0
where
D 2 = 0 1 m 1 0 1 m 2 0 1 m 3 T , h ( t ) = F w + F m 3
Let
h ^ ( t ) = [ w T ( t ) , h T ( t ) ] T , D = [ D 1 D 2 ]
Then, the system Equation (15) can be written as
z ˜ ˙ ( t ) = A z ˜ ( t ) + B u ( t ) + D h ^ ( t ) , z ˜ ( 0 ) = z ˜ 0
The output equation of the system is given as
η ( t ) = C 1 z ˜ ( t ) + E 1 h ^ ( t )
where C 1 and E 1 = [ E w E h ] are given matrices with appropriate dimensions.
Introducing the following quadratic performance index as
J = 0 z ˜ T ( t ) Q z ˜ ( t ) + u T ( t ) R u ( t ) d t
where the weight matrices Q = Q T 0 and R = R T > 0 .
A state feedback control law is designed as
u ( t ) = K z ˜ ( t )
where K is a 1 × 6 gain matrix to be designed.
From (17) and (20), one yields the closed-loop system as
z ˜ ˙ ( t ) = ( A + B K ) z ˜ ( t ) + D h ^ ( t ) , z ˜ ( 0 ) = z ˜ 0
The aim of this paper is to design the optimal guaranteed cost H recoil controller (20) such that
  • The closed-loop system (21) with w ( t ) = 0 and h ( t ) = 0 is asymptotically stable, and J J * , J * is the upper bounds of the performance;
  • The H performance index
    η ( t ) γ h ^ ( t )
    of the closed-loop system can be guaranteed for a given γ > 0 .

3. Design of Optimal Guaranteed Cost H Recoil Controller

In this section, to control the recoil responses of the riser system, an optimal guaranteed cost H control scheme is designed for the system. The existence condition and design algorithm of an optimal guaranteed cost H controller are presented for the riser recoil control system.

3.1. The Existence Condition of the Optimal Guaranteed Cost H Recoil Controller

To analyze the asymptotic stability of the closed-loop system (21), consider the following Lyapunov function as
V ( z ˜ ( t ) ) = z ˜ T ( t ) P z ˜ ( t )
where P is a 6 × 6 positive definite matrix. Then, the time derivative of V ( z ˜ ( t ) ) along system (21) is given as
V ˙ ( z ˜ ( t ) ) = z ˜ T ( t ) Π z ˜ ( t ) + 2 z ˜ T ( t ) P D h ^ ( t )
where
Π = A T P + K T B T P + P A + P B K
To prove asymptotic stability of the closed-loop system (21), set h ^ ( t ) = 0 . Then from (24), one obtains
V ˙ ( z ˜ ( t ) ) = z ˜ T ( t ) ( Y Q K T R K ) z ˜ ( t )
where Y = Π + Q + K T R K .
Note that
V ˙ ( z ˜ ( t ) ) z ˜ T ( t ) Y z ˜ ( t ) λ m i n ( Q + K T R K ) z ˜ ( t ) 2
where λ m i n represents the minimum eigenvalue of the matrix.
To guarantee the asymptotic stability of the system, we require Y < 0 . In fact, by the Schur Complement, it is equivalent to the following matrix inequality as
Π K T R Q * R 0 * * Q < 0
If the matrix inequality (28) holds, then there exists a sufficiently small positive scalar c such that V ˙ ( z ˜ ( t ) ) < c z ˜ ( t ) 2 , z ˜ ( t ) 0 , which indicates that riser system (17) is asymptotically stable.
Integrating both sides of (26) from 0 to , one obtains
0 z ˜ T ( t ) ( Q + K T R K ) z ˜ ( t ) d t < V ( 0 )
That is to say, the quadratic performance index J satisfies J < z ˜ T ( 0 ) P z ˜ ( 0 ) .
Then, we try to find the condition guaranteeing the H performance index (22) of the riser system. Denote
β ( t ) = [ z ˜ T ( t ) h ^ T ( t ) ] T
Then, it follows from (24) and (18) that
V ˙ ( z ˜ ( t ) ) + η T ( t ) η ( t ) γ 2 h ^ T ( t ) h ^ ( t ) β T ( t ) Ψ β ( t )
where
Ψ = Y + C 1 T C 1 P D + C 1 T E 1 * γ 2 I + E 1 T E 1
Note that Ψ < 0 is equivalent to the following matrix inequality as
Π P D K T R Q C 1 T * γ 2 I 0 0 E 1 T * * R 0 0 * * * Q 0 * * * * I < 0
If the above inequality holds, then from (31), one yields
V ˙ ( z ˜ ( t ) ) + η T ( t ) η ( t ) γ 2 h ^ T ( t ) h ^ ( t ) < 0
Further, integrating both sides of (34) and noticing the zero initial state condition of the system yields the H performance index inequality (22) directly.
Note that the inequality (33) implies that (28) holds, which indicates if the inequality (33) holds, then the closed-loop system (21) with h ^ ( t ) = 0 is asymptotically stable, the quadratic performance index (19) is bounded, and the H performance (22) is guaranteed for the friction force h ^ ( t ) .
To obtain the minimum upper bound of the quadratic performance index J, introducing an upper bound μ > 0 , i.e., z ˜ T ( 0 ) P z ˜ ( 0 ) < μ , which is equivalent to the one as
μ z ˜ T ( 0 ) * P 1 < 0
Now, we state a proposition which provides the existence of the optimal guaranteed cost H recoil controller for the riser system.
Proposition 1.
For a given scalar γ > 0 , matrices Q 0 , R > 0 , and the initial state value z ˜ ( 0 ) , if there exist 6 × 6 matrix P > 0 and 1 × 6 matrix K such that the inequalities (33) and (35) hold, then the optimal guaranteed cost H recoil controller (20) for the riser system (6) can be designed, and the gain matrix K of the controller is solvable.

3.2. Computation of the Gain Matrix K of the Controller

To solve the gain matrix K of the optimal guaranteed cost H recoil controller, let
Φ = diag P 1 , I , I , I , I
and pre- and post-multiply both sides of (33) by Φ and its transpose, respectively; denote P ¯ = P 1 and K ¯ = K P 1 . Then, from (33) and (35), one obtains the following linear matrix inequalities as:
Π ¯ D K ¯ T R P ¯ Q P ¯ C 1 T * γ 2 I 0 0 E 1 T * * R 0 0 * * * Q 0 * * * * I < 0
μ z ˜ T ( 0 ) * P ¯ < 0
where Π ¯ = P ¯ A T + K ¯ T B T + A P ¯ + B K ¯ .
Further, the computation of the gain matrix K can be formulated as the following optimization problem:
Minimize μ s . t . ( 36 ) a n d ( 37 ) .
Now, we have the following proposition.
Proposition 2.
For a given scalar γ > 0 , matrices Q 0 , R > 0 , and the initial state value z ˜ ( 0 ) , if there exist 6 × 6 matrix P ¯ > 0 and 1 × 6 matrix K ¯ such that there exists the feasible solution to the optimization problem (38), then the gain matrix K of the optimal guaranteed cost H recoil controller can be solved by K = K ¯ P ¯ 1 .
Remark 1.
Proposition 2 provides a method to solve the optimal guaranteed cost H recoil controller for the riser system. Compared with the classic linear quadratic optimal controller designed in [23], where the external disturbance, i.e., the characteristics of friction force of discharging fluid on the riser are not used to design the controller, while in this paper, based on the introduced H performance index, the characteristics the friction force are used to design the recoil controller. In this sense, the designed optimal guaranteed cost H controller may be better than the optimal controller to reduce the recoil response of the riser, which was verified by simulation results, and presented in the next section.
Remark 2.
Note that there always exists physical limitations of the controller from the application point of view. By adjusting the opening λ ( t ) of the recoil control valve, the flow of hydraulic oil in the tensioner can be controlled. The opening of recoil control valve λ ( t ) can be written as [23]:
λ ( t ) = ρ d q ¯ A r 3 x ˙ p 2 ( t ) 2 S 2 u ( t ) , u ( t ) 0
where the density of the hydraulic oil is represented by ρ d , kg/m 3 ; q ¯ is the flow resistance coefficient; S is the maximum flow area of recoil control valve port, m 2 . In the next section, during the recoil controller design and application, the above relation of the control valve opening λ ( t ) and the control input u ( t ) is considered.

4. Simulation Results

In this section, in the two cases of constant and time-varying friction coefficients, the friction force of the fluid discharge is computed, respectively, and the corresponding discharge length, velocity and total friction of fluid are analyzed. Then, an optimal guaranteed cost H controller and an optimal controller adopted in [23] are designed and used to control the recoil response of the riser.

4.1. Parameters of Deepwater-Drilling Riser System

A drilling riser with a length of 1000 m and six tensioners is applied to simulation study in what follows. The system parameters of the riser are taken from [23], and the main values are given by Table 1. Based on the settings, one yields the system matrices in (17) as
A = 0 1 0 0 0 0 47.7777 0.5794 46.9756 0.1315 0 0 0 0 0 1 0 0 51.0082 0.1428 102.0164 0.3116 51.0082 0.1688 0 0 0 0 0 1 0 0 36.5424 0.1209 36.5424 0.1209 B = 10 6 × 0 2.8153 0 0 0 0 , D 1 = 0 0.8021 0 0 0 0 , D 2 = 10 6 × 0 2.8153 0 3.0569 0 2.1900
In the output Equation (18), the matrices C 1 and E 1 are set as
C 1 = 0.01 0 0 0 0 0 0 0 0.01 0 0 0 0 0 0 0 0.01 0 , E 1 = 0.01 0.01 0 0 0.01 0 0.01 0 0.01 0 0 0.01
In (3), set n = 6 . Then, the response of heave motion of the platform is depicted in Figure 2.

4.2. Analyses of the Fluid Discharge Model

Now, we compare the whole fluid column model from the length, velocity, and total friction. The two cases of constant and time-varying friction coefficients are considered, respectively.
Case I: The friction coefficients, f w and f m of the seawater and drilling mud, are constant. In (2), set f w = 0.002 and f m = 0.004 ; the values are taken from [23].
Case II: The fluid friction coefficients f w and f m are time-varying, and satisfy the Haaland formula as [21]:
1 f w = 1.8 log 6.9 R e w + ε 3.7 D s 1.11
1 f m = 1.8 log 6.9 R e m + ε 3.7 D s 1.11
where ε is the roughness parameter of the riser; R e w and R e m are the Reynolds numbers of the seawater and drilling mud, respectively.
In the above two cases, the friction force of seawater and mud discharge is computed, and the corresponding discharging time, velocity, and the total friction force of the fluid given by Table 2. It can be observed that the discharge velocity of the fluid for the time-varying friction coefficients case is slower than that for the constant case. Consequently, the overall discharge time of the fluid becomes large.
In next subsection, based on the obtained data regarding the different friction force, we investigate the anti-recoil effectiveness of the proposed optimal guaranteed cost H control scheme.

4.3. Effectiveness of the Optimal Guaranteed Cost H Controller

To design an optimal guaranteed cost H controller (OGCHC), the H performance index level γ is set as 0.55, and the weight matrices Q and R in (19) are set as
Q = diag 0.01 , 0.25 , 0.01 , 0.25 , 0.01 , 0.25 , R = 2.5 × 10 12
Then, by Proposition 2, solving the optimization problem (38) yields an OGCHC with gain matrix K as
K OGCHC = 10 5 × 1.0510 4.0249 1.1883 1.1332 0.1114 1.4368
For comparison purposes, a classic optimal controller (OC) developed in [23] was designed. Set weight matrices Q and R of the linear quadratic performance index as:
Q = 10 4 × diag 1 , 1 , 1 , 1 , 1 , 1 , R = 10 7
Then the gain matrix of an OC is computed as
K OC = 10 5 × 3.6435 4.4226 0.9550 1.8504 0.6362 2.6558
As the OGCHC and OC are applied to the riser system, and the controller constraint (39) is considered at the same time, the curves of the displacement responses of the three mass blocks of the riser and control force are depicted—see Figure 3, Figure 4, Figure 5 and Figure 6 for the case of constant friction coefficients, and Figure 7, Figure 8, Figure 9 and Figure 10 for the case of time-varying friction coefficients, respectively.
To further compare the two controllers quantitatively, the quadratic performance index of the riser system with OGCHC and OC are computed. Note that the performance index (19) are weight-matrix-dependent, and the designed OGCHC and OC are based on a different weight matrix pair, ( Q , R ) . In this situation, the performance index (19) is modified as two weight-matrix-independent performance indices:
J z ˜ = 1 t f 0 t f z ˜ T ( t ) z ˜ ( t ) d t
J u = 1 t f 0 t f u 2 ( t ) d t
Clearly, J z ˜ and J u are state-related and control-related root mean square values of the riser system, respectively. Denote the average displacement peak value of three mass blocks of the riser by d ^ . Then, the above indices of the riser system without control, with OC and OGCHC in the two cases of constant and time-varying friction coefficients are computed and listed in Table 3.
It can be found from Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10, and Table 3, that the controllers OGCHC and OC can reduce the recoil response of the riser system efficiently; however, the OGCHC is more economical than the OC. In fact, the average displacement amplitudes of the three mass blocks of the riser are almost the same, while the control cost by the OGCHC is less than that by the OC. In addition, the reduction in the average displacement peak value of the riser under OGCHC is larger than that under OC.
Note that if the riser is compressed in the recoil process, the riser system may be collapsed. To further verify the safety of the riser under the designed recoil controller in this paper, the elongations of springs of the riser are computed. The elongation curves of the springs under the above two cases are shown in Figure 11 and Figure 12, respectively. It can be found that the spring elongations are both greater than zero, which means that the riser is not compressed during the recoil process thereby ensuring the safety of the riser system.

5. Conclusions

The recoil control issue of the deepwater-drilling riser system has been studied in this paper. First, an incremental model of the riser system has been established by considering the heave motion, linear tension force, and drilling mud discharge model. Then, an optimal guaranteed cost H recoil control scheme has been presented for the riser system, and the sufficient conditions for the asymptotic stability of the closed-loop system have been derived based on Lyapunov stability theory. In addition, under the designed recoil controller, the H performance index and the minimum upper bound of the quadratic performance index of the riser can be guaranteed.
Simulation results show that the designed optimal guaranteed cost H recoil controller is more efficient than the exiting quadratic optimal recoil controller to refrain the recoil movements of the riser. In fact, the quadratic performance index of the riser with the former is smaller than that with the latter. It is also found that under the designed recoil controller, the compression of the riser has been avoided, and the safety of the riser has been ensured.
This paper provides an active recoil control scheme for the riser system in the event of emergency disconnection; however, there are still important issues regarding recoil analysis and control topics to be investigated in the future. It is important to further explore and understand the recoil mechanism of the riser, to provide an in-depth development potential and reasonable active-recoil control schemes.

Author Contributions

Conceptualization and methodology, B.-L.Z.; software, B.-L.Z. and Y.-T.S.; validation, Y.-D.Z., W.Z. and H.S.; formal analysis and investigation, Y.-T.S.; resources, Y.-D.Z., W.Z. and H.S.; writing—original draft preparation, Y.-T.S.; writing—review and editing, B.-L.Z.; visualization, W.Z. and H.S.; supervision, B.-L.Z. and W.Z.; project administration and funding acquisition, B.-L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant 61773356.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data presented in this paper are freely available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, C.; Li, H.-N.; Hao, H.; Bi, K.; Chen, B.-K. Seismic fragility analyses of sea-crossing cable-stayed bridges subjected to multi-support ground motions on offshore sites. Eng. Struct. 2018, 165, 441–456. [Google Scholar] [CrossRef]
  2. Li, C.; Liu, Y.; Li, H.-N. Fragility assessment and optimum design of a steel-concrete frame structure with hybrid energy-dissipated devices under multi-hazards of earthquake and wind. Eng. Struct. 2021, 245, 112878. [Google Scholar] [CrossRef]
  3. Pan, H.-Y.; Li, H.-N.; Li, C. Seismic behaviors of free-spanning submarine pipelines subjected to multi-support earthquake motions within offshore sites. Ocean Eng. 2021, 237, 109606. [Google Scholar] [CrossRef]
  4. Wu, Z.; Liu, J.; Liu, Z.; Lu, Z. Parametrically excited vibrations of marine riser under random wave forces and earthquake. Adv. Struct. Eng. 2016, 19, 449–462. [Google Scholar] [CrossRef]
  5. Wang, J.; Xu, L.; Cao, J.; Sheng, L.; Li, C. Numerical and experimental investigation for extreme storm-safe drilling riser. Ship. Offshore Struc. 2021, 1–13. [Google Scholar] [CrossRef]
  6. Chang, Y.; Chen, G.; Sun, Y.; Xu, L.; Peng, P. Nonlinear dynamic analysis of deepwater drilling risers subjected to random loads. China Ocean Eng. 2008, 22, 683–691. [Google Scholar]
  7. Rocha, L.A.S.; Junqueira, P.; Roque, J.L. Overcoming deep and ultra deepwater drilling challenges. In Proceedings of the Offshore Technolog Conference, Houston, TX, USA, 5–8 May 2003. OTC-15233-MS. [Google Scholar]
  8. Su, K.; Butt, S.; Yang, J.; Qiu, H. Coupled dynamic analysis for the riser-conductor of deepwater surface BOP drilling system. Shock. Vib. 2018, 2018, 6568537. [Google Scholar] [CrossRef] [Green Version]
  9. Nguyen, D.H.; Nguyen, D.T.; Quek, S.T.; Sorensen, A.J. Position-moored drilling vessel in level ice by control of riser end angles. Cold Reg. Sci. Technol. 2011, 66, 65–74. [Google Scholar] [CrossRef]
  10. Yamamoto, M.; Morooka, C.K. Feedback control system for blow-out preventer positioning. Appl. Ocean Res. 2019, 82, 362–369. [Google Scholar] [CrossRef]
  11. Lang, D.W.; Real, J.; Lane, M. Recent developments in drilling riser disconnect and recoil analysis for deepwater applications. In Proceedings of the ASME 2009 28th International Conference on Ocean, Offshore and Arctic Engineering, Honolulu, HI, USA, 31 May–5 June 2009; pp. 305–318. [Google Scholar]
  12. Liu, X.; Li, Y.; Zhang, N.; Sun, H.; Chang, Y.; Chen, G.; Xu, L.; Sheng, L. Improved axial dynamic analysis of risers based on finite element method and data-driven models. Ocean Eng. 2020, 214, 107782. [Google Scholar] [CrossRef]
  13. Pestana, R.G.; Roveri, F.E.; Franciss, R.; Ellwanger, G.B. Marine riser emergency disconnection analysis using scalar elements for tensioner modelling. Appl. Ocean Res. 2016, 59, 83–92. [Google Scholar] [CrossRef]
  14. Wang, Y.; Gao, D.; Fang, J. Axial dynamic analysis of marine riser in installation. J. Nat. Gas Sci. Eng. 2014, 21, 112–117. [Google Scholar] [CrossRef]
  15. Wang, Y.; Gao, D. Recoil analysis of deepwater drilling riser after emergency disconnection. Ocean Eng. 2019, 198, 106406. [Google Scholar] [CrossRef]
  16. Wang, X.; Liu, X.; Liu, Z.; Qu, N.; Hu, P.; Chang, Y.; Chen, G.; Li, C. Dynamic recoil response of tensioner and riser coupled in an emergency disconnection scenario. Ocean Eng. 2022, 247, 110730. [Google Scholar] [CrossRef]
  17. Liu, X.; Liu, Z.; Wang, X.; Qu, N.; Zhang, N.; Chang, Y.; Chen, G. An intelligent recoil controller for riser system based on fuzzy control theory. Int. J. Nav. Archit. Ocean Eng. 2022, 14, 100439. [Google Scholar] [CrossRef]
  18. Young, R.D.; Hock, C.J.; Karlsen, G.; Miller, J.E. Analysis and design of anti-recoil system for emergency disconnect of a deepwater riser: Case study. In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 4–7 May 1992. OTC-6891-MS. [Google Scholar]
  19. Grønevik, A. Simulation of Drilling Riser Deiconnection-Recoil Analysis. Master’s Thesis, Norwegian Uniersity of Science and Technology, Trondheim, Norway, 2013. [Google Scholar]
  20. Li, C.; Fan, H.; Wang, Z.; Ji, R.; Ren, W.; Feng, X. Two methods for simulating mud discharge after emergency disconnection of a drilling riser. J. Nat. Gas Sci. Eng. 2016, 28, 142–152. [Google Scholar] [CrossRef] [Green Version]
  21. Meng, S.; Che, C.-D.; Zhang, W.-J. Discharging flow effect on the recoil response of a deep-water drilling riser after an emergency disconnect. Ocean Eng. 2018, 151, 199–205. [Google Scholar] [CrossRef]
  22. Wang, X.; Liu, X.; Zhang, S.; Chen, G.; Chang, Y. Study on mud discharge after emergency disconnection of deepwater drilling risers. J. Pet. Sci. Eng. 2020, 190, 107105. [Google Scholar] [CrossRef]
  23. Liu, X.; Liu, Z.; Wang, X.; Zhang, N.; Qiu, N.; Chang, Y.; Chen, G. Recoil control of deepwater drilling riser system based on optimal control theory. Ocean Eng. 2021, 220, 108473. [Google Scholar] [CrossRef]
  24. Li, H.; Wang, J.; Wu, L.; Lam, H.K.; Gao, Y. Optimal guaranteed cost sliding-mode control of interval Type-2 fuzzy time-delay systems. IEEE Trans. Fuzzy Syst. 2018, 26, 246–257. [Google Scholar] [CrossRef] [Green Version]
  25. Ye, D.; Song, T. Decentralized reliable guaranteed cost control for large-scale nonlinear systems using actor-critic network. Neurocomputing 2018, 320, 121–128. [Google Scholar] [CrossRef]
  26. Wu, J.; Peng, C.; Zhang, J.; Yang, M.; Zhang, B.-L. Guaranteed cost control of hybrid-triggered networked systems with stochastic cyber-attacks. ISA Trans. 2020, 104, 84–92. [Google Scholar] [CrossRef]
  27. Xie, X.; Lam, J. Guaranteed cost control of periodic piecewise linear time-delay systems. Automatica 2018, 94, 247–282. [Google Scholar] [CrossRef]
  28. Zhao, Y.; Guo, G.; Ding, L. Guaranteed cost control of mobile sensor networks with Markov switching topologies. ISA Trans. 2015, 58, 206–213. [Google Scholar] [CrossRef]
  29. Wu, Z.-G.; Dong, S.; Shi, P.; Su, H.; Huang, T.; Lu, R. Fuzzy-model-based nonfragile guaranteed cost control of nonlinear markov jump systems. IEEE Trans. Syst. Man Cybern. Syst. 2017, 47, 2388–2397. [Google Scholar] [CrossRef]
  30. Zhang, J.; Peng, C. Guaranteed cost control of uncertain networked control systems with a hybrid communication scheme. IEEE Trans. Syst. Man Cybern. Syst. 2018, 50, 3126–3135. [Google Scholar] [CrossRef]
  31. Yang, C.; Du, J.; Cheng, Z.; Wu, Y.; Li, C. Flexibility investigation of a marine riser system based on an accurate and efficient modelling and flexible multibody dynamics. Ocean Eng. 2020, 207, 107407. [Google Scholar] [CrossRef]
Figure 1. A riser-tension recoil control system model. Adapted with permission from [23,31]. 2020, 2021, Elsevier.
Figure 1. A riser-tension recoil control system model. Adapted with permission from [23,31]. 2020, 2021, Elsevier.
Applsci 12 03945 g001
Figure 2. Response of heave motion of platform.
Figure 2. Response of heave motion of platform.
Applsci 12 03945 g002
Figure 3. Displacement 1 of the riser without control, with OC and OGCHC under case I.
Figure 3. Displacement 1 of the riser without control, with OC and OGCHC under case I.
Applsci 12 03945 g003
Figure 4. Displacement 2 of the riser without control, with OC and OGCHC under case I.
Figure 4. Displacement 2 of the riser without control, with OC and OGCHC under case I.
Applsci 12 03945 g004
Figure 5. Displacement 3 of the riser without control, with OC and OGCHC under case I.
Figure 5. Displacement 3 of the riser without control, with OC and OGCHC under case I.
Applsci 12 03945 g005
Figure 6. Control force with OC and OGCHC under case I.
Figure 6. Control force with OC and OGCHC under case I.
Applsci 12 03945 g006
Figure 7. Displacement 1 of the riser without control, with OC and OGCHC under case II.
Figure 7. Displacement 1 of the riser without control, with OC and OGCHC under case II.
Applsci 12 03945 g007
Figure 8. Displacement 2 of the riser without control, with OC and OGCHC under case II.
Figure 8. Displacement 2 of the riser without control, with OC and OGCHC under case II.
Applsci 12 03945 g008
Figure 9. Displacement 3 of the riser without control, with OC and OGCHC under case II.
Figure 9. Displacement 3 of the riser without control, with OC and OGCHC under case II.
Applsci 12 03945 g009
Figure 10. Control force with OC and OGCHC under case II.
Figure 10. Control force with OC and OGCHC under case II.
Applsci 12 03945 g010
Figure 11. The elongation of springs under case I.
Figure 11. The elongation of springs under case I.
Applsci 12 03945 g011
Figure 12. The elongation of springs under case II.
Figure 12. The elongation of springs under case II.
Applsci 12 03945 g012
Table 1. Parameters of riser system.
Table 1. Parameters of riser system.
ParametersValueParametersValueParametersValue
P L 0 150,000 V L 0 2.25 N b 42
P H 0 3,050,000 V H 0 4.28 D s 0.4826
M slip 28,082 A p 0.2463 ε 5 × 10 5
M LMRP 129,496 A r 0.2048 ν w 1 10 4
M buoyancy 23,366 A e 0.0405 ν m 2 1.15 × 10 6
1νw is the viscosity of seawater. 2 νm is the viscosity of mud.
Table 2. Discharging time, velocity, and friction force of the fluid discharge under two different cases.
Table 2. Discharging time, velocity, and friction force of the fluid discharge under two different cases.
CaseDischarging Time (s)Velocity (m/s 2 )Friction Force ( 10 5 N)
Case I51.1324.30945.6475
Case II114.2911.98158.0178
Table 3. Performance indices of the riser without control, with OC and OGCHC.
Table 3. Performance indices of the riser without control, with OC and OGCHC.
CasesControllers J z ˜ (m) d ^ (m) J u ( 10 4 N)
INo control0.23626.2985
OC [23]0.18884.30661.0859
OGCHC0.19123.80651.0625
IINo control0.33024.9055
OC [23]0.30013.67410.9758
OGCHC0.30492.87220.8956
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, Y.-T.; Zhao, Y.-D.; Zhang, B.-L.; Zhang, W.; Su, H. Recoil Control of Deepwater-Drilling Riser with Optimal Guaranteed Cost H Control. Appl. Sci. 2022, 12, 3945. https://doi.org/10.3390/app12083945

AMA Style

Sun Y-T, Zhao Y-D, Zhang B-L, Zhang W, Su H. Recoil Control of Deepwater-Drilling Riser with Optimal Guaranteed Cost H Control. Applied Sciences. 2022; 12(8):3945. https://doi.org/10.3390/app12083945

Chicago/Turabian Style

Sun, Yue-Ting, Yan-Dong Zhao, Bao-Lin Zhang, Wei Zhang, and Hao Su. 2022. "Recoil Control of Deepwater-Drilling Riser with Optimal Guaranteed Cost H Control" Applied Sciences 12, no. 8: 3945. https://doi.org/10.3390/app12083945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop