Abstract

The landing string is an important component of deepwater riserless drilling systems. Determination of the dynamic characteristics of the landing string plays an essential role in its design for ensuring its safe operation. In this paper, a dynamic model is developed to investigate the dynamic response characteristics of a landing string, where a landing string in a marine environment is modeled as a flexible slender tube undergoing coupled transverse and axial motions. The heaving motion of the drilling platform is taken as the upper boundary condition and the motion of the drilling bit caused by the interaction between the rock and the bit as the lower boundary condition. A semiempirical Morison equation is used to simulate the effect of the load imposed by the marine environment. The dynamic model, which is nonlinearly coupled and multibody, is discretized by a finite element method and solved by the Newmark technique. Using the proposed model, the dynamic responses of the displacement, axial force, and moment in the landing string are investigated in detail to find out the influences of driving depth of surface catheter, platform motion, bit movement, and marine environment on the dynamical characteristics of the landing string.

1. Introduction

Riserless drilling is a new deepwater drilling technology in which the drilling platform no longer bears the interaction load induced by the ocean waves and the ocean currents which are imposed on the riser and the mass buoyancy block. Compared with the traditional deepwater drilling method with riser, this technique has many remarkable advantages, such as reducing the requirement of storage capacity of drilling platform, greatly shortening the construction period, saving the cost of drilling, and improving the safety of drilling operation [1]. However, there are still some thorny problems in the technology. One salient issue of them is the dynamic load response in the landing string used in a jet technique for lowering and lifting the surface catheter and the BHA (bottom hole assembly) (Figure 1). Due to the lacking of protection from riser as well as the exposure to the ocean environment directly, the landing string bears very complex patterns of external loads, such as working load, marine environmental load, drilling platform movement, and their combinations. Therefore, the internal force and deformation of the landing string are often complex and changeable, making the pipe string prone to be instability, fracture, and fatigue failure and causing serious accidents and economic loss [2]. Thus, in order to ensure the safety and reliability of deepwater riserless drilling operation, it is very necessary to find out the dynamic characteristics of the landing string, including vibration modes, critical speed of the bit, deformation, and internal load response under the actions of environmental load and platform movement [35].

Many works, using theoretical analyses [6, 7], numerical simulations [8], and experimental investigations [9, 10], have been conducted to investigate the dynamics behaviour of on-land drilling string to determine the change regulations of the friction, torsion, and axial force under various working conditions. These research results guided the operation of on-land drilling operation. However, there is a clear difference between offshore drilling and on-land drilling. For the former, the traditional drilling string is in contact with the riser in motion. For the latter, the drilling string is in direct contact with the fixed casing. Therefore, the dynamic behaviour of offshore drilling pipe string is more complex than that of on-land drilling pipe string [11]. Most of the research studies on the mechanics analysis of marine pipe string focused on the dynamic behaviour of the riser [1215]. Some works [16, 17] investigated the interaction between the completion pipe string (drilling pipe or test pipe) and the riser.

In comparison with the conventional deepwater drilling with a riser, the riserless landing pipe system is obviously different in structure and load [18]. At first, the upper end of the pipe string is fixed directly to the drilling ship, so the bending deformation of the pipe string cannot be reduced by applying the top tensioning force as that used in conventional deepwater riser. Secondly, since the landing string is not constrained by the riser in the radial direction, it is prone to buckle under the disturbance of the transverse load as the axial tension in the pipe is small. Therefore, the landing string is more prone to damage induced by large deformation and complex vibration [19] and has attracted the attentions of many scholars [2022].

An approach to the design of deepwater riserless landing strings was proposed by Adams [23], which focuses on the selection of pipe joints and tool joints. The design, manufacture, testing, and operation of landing strings have also been discussed by Breihan et al. [24], Jellison et al. [25], and Cantrell et al. [26]. However, in these studies, dynamic analysis did not get enough attention. Sathuvalli et al. [27] pointed out that deepwater landing strings can usually withstand great axial loads, so slip collapse is the main factor to restrict the operational capability of such strings. They also presented a new method for calculating the slip collapse load. Bradford et al. [28] carried out five different types of experiments on slip collapse. Everage et al. [29] and Zheng et al. [30] investigated methods for calculating the longitudinal vibration load on a landing string induced by platform heave motion during lowering of the casing. Samuel et al. [31] presented a simulation approach for hook load, torque, and drag calculations and investigated various parameters’ influence on these load, such as wellhead offset from the rig center, wellbore inclination, curvature, wellbore torsion, angle of entry into the wellhead, wave forces, and ocean loop currents. The dynamic problems involved in these studies are either for the longitudinal vibration of the string or for the transverse vibration. Han and Benaroya [32] pointed out that the effects of longitudinal and transverse coupling vibration (LTCV) are not negligible for slender columns and investigated the LTCV behaviour of a tension leg platform. The LTCV behaviour of a beam simply supported at both ends was examined by Xing and Liang [33]; Hu and Yang [34]; Ghayesh and Farokhi [35]; and Zhao and Wu [36]. Liu et al. [37] established a LTCV model to investigate the effects of the main parameters such as marine load, platform motion, and suspending length on the dynamic behaviour of a deepwater hard suspension riser under emergency evacuation conditions.

Besides the effect caused by the nonlinear LTCV of a pipe string or a beam discussed by Han and Benaroya [32]; Zhao and Wu [36]; and Liu et al. [37], the landing string in riserless drilling is affected simultaneously by the combined actions of the platform, drilling bit, ocean environment, driving depth, and so on. Therefore, its vibrational behaviour is strongly nonlinear. However, up to now, there has been very limited research on this issue. This paper presents a nonlinear dynamic model of a deepwater landing string under the combined action of platform movement, drill movement, and marine load, taking into account the LTCV and driving depth effects. Based on the model, the nonlinear vibration modes and mechanical response mechanism of a landing string are investigated in detail, aimed at providing a theoretical basis for the design, construction, and operation of deepwater landing strings in riserless drilling.

2. The Proposed Nonlinear Dynamic Model

In this section, the Hamilton variational principle is used to establish the nonlinear vibration control equation of the landing string, based on the energy analysis of system.

2.1. Basic Assumptions

According to the structure and working conditions of the deepwater drilling string system shown in Figure 1, the following assumptions are made to derive the dynamic governing equations of the system:(1)The material property is homogeneous and isotropic, and the elasticity remains linear during movement and deformation.(2)The influence of the joints on the bending rigidity of the landing string is not taken into account.(3)The effect of structural damping on the model of the riser is neglected.(4)The directions of ocean currents and the deformation of the drilling string are assumed to lie in the same plane.(5)The landing string is connected rigidly with the platform and the surface catheter. The transverse displacement of the surface catheter is restricted.(6)The deformation of the surface catheter and BHA is ignored.(7)The surface catheter keeps vertical during driving.

2.2. Energy Equation of the System

The total length of the landing string is denoted by . The elements of the string are labeled by their locations in the undeformed configuration . The reference and current configurations of the midline of the string are shown in Figure 2. The reference configuration is indicated by the vertical straight dashed line and the current configuration by the curved dashed line. Similar coordinate system can also be found in Han and Liu’s work [32, 37].

According to Han and Benaroya [32], the Green strain based on Kirchhoff’s hypothesis can be written aswhere and are the deflections of the midline in the and directions, respectively; is the transverse distance from the midline to the point of interest on the cross section in the reference frame. Assuming a symmetrical cross section, is also the average deflection of the beam element.

According to the definition of strain energy and (1), the strain energy can be given as follows and its detailed derivation is presented in Appendix A:where the primes indicate derivatives with respect to .

The total kinetic energy of the system consists of two parts: the kinetic energy of the landing string and the kinetic energy of the combined catheter plus BHA.

The kinetic energy of the landing string is given bywhere the dots indicate derivatives with respect to time.

The total kinetic energy of the landing system can be expressed as follows and its detailed derivation can be found in Appendix B:where is the Rayleigh rotational term, which is the kinetic energy due to the rotation of the cross section and is usually small compared with the kinetic energy due to translation, .

2.3. Control Equation and Boundary Conditions via Hamilton’s Principle

Hamilton’s principle is used to formulate the governing dynamical equations of the landing string system. The Lagrangian integrated over time, , is given byand the virtual work done by the transverse forces and axial forces is given by

Using Hamilton’s principle, control equation and boundary condition can be formulated asand their detailed derivation can be found in Appendix C.

Boundary conditions:where is the heaving displacement of the drilling platform.

The hydrodynamic force is simulated by the Morison equation whose concise form is given by Han and Benaroya [32] as

3. Model Solution and Verification

In this section, the control equation is discretized by the finite element method and solved by the Newmark technique. The validity of the model is tested using the ANSYS software.

3.1. Finite Element Scheme for Nonlinear Governing Equations

We use a linear Lagrange function and a cubic Hermitian function to approximate the transverse displacement and longitudinal displacement of a beam element:where

Using (10) and the Galerkin weighted residual method, the finite element scheme for the nonlinear governing equation (7) is given as follows:where is the generalized coordinate vector; , , and are, respectively, the stiffness, mass, and load matrices, which, respectively, consist of element stiffness, mass, and load matrices, , , and shown in Appendix D.

Using Newmark’s integral method in Appendix E, (12) yields the following equation:

The system responses, including the displacements, velocities, and accelerations at the instant , can be obtained by substituting the solution of (13) into (E.1) and (E.2) in Appendix E.

3.2. Verification of the Nonlinear Dynamic Model

A longitudinal-transverse coupling dynamic analysis code (LTND) based on the above model is developed. The validity of the model and the code is tested using the ANSYS software, which does not take into account the LTCV effect (Figure 3). The basic parameters of the example are shown in Table 1.

The longitudinal and transverse displacement responses of the middle section of the standpipe under two kinds of working conditions are simulated, respectively, by the coupling method and ANSYS. ANSYS has adopted BEAM188 element, and case 1 and case 2 have divided 10 elements and 100 elements, respectively. The greatest difference between the two cases is the slenderness ratio, which is much larger in case 2. It is worth noting that as the slenderness ratio increases, the LTCV effect becomes more obvious [32].

Figures 4(a) and 5(a) show that, for longitudinal vibration, the results obtained with the coupling method are close to those obtained using ANSYS. For the larger slenderness ratio, the vibrational amplitude obtained with the coupling method is greater than that from ANSYS. Figures 4(b) and 5(b) show that, for the smaller slenderness ratio, the transverse vibrational response obtained with the coupling method is close to that from ANSYS. However, when the slenderness ratio is increased, the difference between the results becomes very large. These comparisons show that LTCV behaviour depends significantly on the slenderness ratio, with transverse vibrations being more strongly affected by this ratio than longitudinal vibrations.

4. LTCV Dynamic Characteristic Analysis

Using the developed LTND code, the LTCV characteristics of a practical deepwater riserless landing string system are investigated. The structure and environment parameters are presented in Table 2. The motion of a bit is simulated by the following formula [38]:wherewhere ; is the initial displacement, in this paper .

The lateral displacement of the surface casing is constrained, and the longitudinal friction can be calculated by the following formula:

The calculation formulas for formation pressure and friction factor can be found in literature [39, 40], respectively, and are listed as follows:

4.1. Vibration Modal Analysis

From the submatrix of the stiffness matrix of the nonlinear system, , , and , it can be seen that the stiffness matrix is time-variant. For an arbitrarily small-time quantum, the frequencies and modes are assumed not to change with time and their calculation equations are given aswhich shows that the frequencies and modes in the nonlinear system in this paper are also time-variant, where is the vector of mode shape. Using (18), the first four order frequencies and modes of the deepwater landing string have been calculated and, respectively, are shown in Figures 6 and 7.

Figure 6 shows that the frequency of transverse vibration varies greatly with time and increases with the increase of order. The mode diagram in Figure 7 reflects the same trend. At the initial stage of the vibration, the frequency of each order rises sharply and then changes up and down at a certain value. The main reason is that at the beginning of the vibration, the longitudinal tension as well as the stiffness of the landing string increases sharply due to the string’s flexes. Later, the landing string vibrates in a certain bending shape, and the axial force and stiffness change periodically.

To further illustrate the effect of LTCV on the vibration frequency of a landing string, the vibration frequency without considering LTCV is also shown in Figure 6. It is demonstrated that the frequencies taking into account the LTCV are much larger than those that do not consider coupling. The increase and changeability of the vibration frequency make the landing string more inclined to resonate with the platform and bit.

4.2. Dynamic Response Analysis

In this section, the influence of LTCV, shear flow velocity, platform motion, drilling bit speed, and driving depth on the dynamic response analysis characteristic of a landing string is analyzed.

4.2.1. LTCV Effect

To investigate the effect of LTCV on the dynamical behaviour of the landing string, a 1000 m long landing string in a shear flow with 1.2 m/s maximum velocity at sea level and 0.036 m/s minimum velocity at the seabed is considered. The velocity of the shear flow is approximately considered to change linearly. Other parameters are shown in Table 2. As can be seen from Figure 8, longitudinal-transverse coupling has little effect on the amplitude of longitudinal vibrations compared with the uncoupled case. However, when coupling is taken into account, the first-order frequency of longitudinal vibrations obviously decreases, and, as can be seen from Figure 9, the amplitude and the first-order frequency of transverse vibrations are both significantly affected. Therefore, in the dynamic analysis of a deepwater landing string, the influence of nonlinear coupling of vertical and horizontal motions cannot be ignored.

4.2.2. Effect of Shear Flow Velocity

The dynamic behaviour of a landing string of length 1500 m is analyzed to investigate the influence of shear velocity on vibration of the pipe string. Figure 10 (where is the flow velocity at sea level and that at the seabed) shows that the shear velocity has little influence on the longitudinal displacement and frequency of the tube, and the longitudinal vibrational displacement and frequency are controlled mainly by the motion of the platform. Figure 11 shows that the shear velocity has a strong influence on the transverse displacement of the tube: the higher the shear velocity, the greater is the transverse displacement at the center of the tube.

As shown in Figure 12, the maximum dynamic tensile force appears in the upper and middle parts of the landing string, whereas Figure 13 shows that the maximum dynamic compressive force appears in the lower part. The two axial forces decrease with increasing shear velocity. Figure 14 shows that the maximum dynamic moment increases with increasing water depth. The main reason for this is that the tension in the tube gradually decreases with increasing water depth owing to the influence of the pipe weight, and thus, bending gradually becomes more and more significant.

4.2.3. Effect of Platform Motion

The influence of platform motion on the dynamic behaviour of the landing string is analyzed under the condition that other parameters remain unchanged. As can be seen from Figures 15 and 16, the longitudinal vibration at the center of the pipe string is controlled mainly by the longitudinal motion of the platform. The longitudinal vibration of the platform has no obvious influence on the transverse displacement but does have a significant influence on the higher-order vibrational frequencies of transverse vibration. Figures 1719 show that longitudinal vibration of the platform will increase the maximum dynamic axial force and bending moment of the whole landing string, especially its lower part. Therefore, platform motion has an important influence on the safety of the landing string and the working efficiency of the drilling bit.

4.2.4. Effect of Drilling Bit Speed

The influence of drilling bit speed on the dynamic behaviour of the landing string is investigated under the condition that other parameters remain unchanged. Figures 20 and 21 show that the bit speed has no significant effect on the displacement response of the landing string. Figures 22 and 23 show that the higher the drill speed, the greater is the dynamic axial force on the pipe string. Figure 24 shows that the bit speed has no significant influence on the maximum dynamic bending moment of the whole pipe string.

4.2.5. Effect of Driving Depth

The influence of driving depth on the dynamic behaviour of the landing string is analyzed under the condition that other parameters remain unchanged. The depth from the platform in the section contains the driving depth and the length of the landing string. Figures 25 and 26 show that the increase of the driving depth has no significant influence on the vibration displacement in the middle of the landing string, whereas Figures 2729 show that the maximum dynamic axial force and moment are influenced significantly by the driving depth. The maximum dynamic axial force and moment increase with the increase of driving depth; the main reason is that the increasing driving depth increases the longitudinal friction force on the pipe section in the formation section. Figures 2729 also show that the driving depth has the greater influence on the dynamic internal forces in the middle-upper part of the landing string than those in the lower part.

5. Conclusions

(1)A longitudinal-transverse coupling dynamic analysis code (LTND) for deepwater landing string is developed and tested by comparing its simulation result with that from finite element software ANSYS. The comparison and analysis show that the code has good calculation accuracy, and the LTCV vibration is obvious for the pipe with large slenderness ratio. Deepwater landing strings are typical members with large slenderness ratio. The vibration model and displacement response analysis of an actual landing string show that the LTCV effect is significant and cannot be ignored to a landing string, especially for its transverse vibration.(2)The effects of shear flow, platform motion, bit moment, and driving depth on the dynamic response of the actual landing string are investigated, by examining the displacement, dynamic axial force, and dynamic moment. It is shown that the axial force in the pipe string is significantly influenced by platform motion, shear flow velocity, and bit movement. The greater the amplitude of platform motion or the drill speed, the greater is the maximum dynamic axial force. By contrast, the higher the shear flow velocity, the smaller is the maximum dynamic axial force. It is also found that the driving depth has the greater influence on the dynamic internal forces in the middle-upper part of the landing string than those in the lower part.(3)In actual drilling, a special emergency disconnection joint is present in the landing string near the underwater wellhead. When a threat such as a typhoon is imminent, this joint is automatically disconnected by moving the platform to increase the axial tension near the joint to a certain value. However, it is important that this emergency joint remains connected under normal drilling conditions. Therefore, under normal conditions, the axial tension must be less than the set value, which requires particular measures to be taken, such as using a heave compensation device to limit the movement of the platform and controlling the drilling speed of the bit.

Appendix

A. Derivation of

The second Kirchhoff stress, ignoring the Poisson effect, is given by

The general expression for the strain energy is

Using the expressions for the second Kirchhoff stress and the Green strain, we obtain

Expanding the integrand, we obtain

On integrating over the area, the terms that are odd functions of (the second and fifth terms) disappear because of the symmetry of the cross section, and becomes the area moment of inertia about the neutral axis. The remaining terms are functions of and only, that is, (2) can be obtained.

B. Derivation of

According to Kirchhoff’s hypothesis, we can writewhere , and are the displacement field functions corresponding to the coordinate systems , and .

Replacing the axial displacements and in (B.1) using (3), we can writewhere, on integration over the cross section, the term with factor has disappeared and the term with factor has become as before. The kinetic energy of the combined catheter plus BHA is given by

C. Formulation of Hamilton’s Principle

The variation of (5) is given by

Integrating by parts and noting that the variations in time at the two endpoints are assumed to be zero, we obtain

According to Hamilton’s principle, , and thus, the nonlinear coupled control equations of the landing string system and its boundary condition are obtained.

D. , , and in (12)

where

E. Newmark’s Integral Method

Newmark’s integral method is used to solve (12), in which the displacement and velocity at the instant can be expressed, respectively, aswhere is the time step length and and are parameters related to the precision and stability and under normal circumstances take the values and .

Using (D.2), the velocity and acceleration at the instant can be expressed approximately by the responses at the instant and the displacement at the instant :

Nomenclature

:Young’s modulus
:Cross-sectional area of the drilling string
:Volume of the deformed landing string
:Area moment of inertia of the cross section around the neutral axis
:Length of drilling string
:Material density
:Added mass coefficient
:Drag coefficient for an oscillating cylinder in still water
:Density of the surrounding fluid
:Cross section of the displaced volume
:Length of one finite element
:Mass of drilling string with unit length
:Mass of the combined catheter plus BHA
:Outer radius of the drill string
:Transverse force of unit length on the drill string
:Axial force of unit length on the drill string
:Heave displacement of the drilling platform
:Angular frequency of the system
:Natural frequency of the drill string
:Natural frequency of damped system
:Equivalent stiffness of drill string
:Equivalent mass of drill string
:Damping ratio of the system
:Drilling speed
:System viscosity coefficient
:Formation pressure acting on surface casing
:Area of the surface casing subjecting formation pressure
:Casing formation friction factor
:Overlying formation pressure
:Vertical depth of a sea
:Formation density
:Total vertical depth including sea section and formation section
:Pressure gradient of overlying formation
:Static fluid gradient of sea section.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Jun Liu was responsible for conception and design of the study, Hongliang Zhao was involved in acquisition of data, Simon X. Yang performed analysis and interpretation of data, Guorong Wang drafted the manuscript, and Qingyou Liu revised the manuscript critically for important intellectual content. Jun Liu, Hongliang Zhao, Simon X. Yang, Qingyou Liu, and Guorong Wang approved the version of the manuscript to be published.

Acknowledgments

This work has been supported in part by the National Natural Science Foundation of China (Grant No. 51875489): “Study on axial-lateral-torsion nonlinear coupling vibration characteristics of deepwater riserless drillstring,” and the Open Fund (OGE201701-01) of Key Laboratory of Oil & Gas Equipment, Ministry of Education (Southwest Petroleum University).