Next Article in Journal
Effect of Direct-Contact Ultrasonic and Far Infrared Combined Drying on the Drying Characteristics and Quality of Ginger
Next Article in Special Issue
Numerical Simulation of Proppant Transport in Transverse Fractures of Horizontal Wells
Previous Article in Journal
Bentonite Modified with Surfactants—Efficient Adsorbents for the Removal of Non-Steroidal Anti-Inflammatory Drugs
Previous Article in Special Issue
A New Bottom-Hole Assembly Design Method to Maintain Verticality and Reduce Lateral Vibration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wellbore Pressure Modeling for Pumping and Tripping Simultaneously to Avoid Severe Pressure Swab

1
School of Petroleum Engineering, Yangtze University, Wuhan 430100, China
2
National Engineering Research Center for Oil & Gas Drilling and Completion Technology, Yangtze University, Wuhan 430100, China
3
Hubei Key Laboratory of Oil and Gas Drilling and Production Engineering, Yangtze University, Wuhan 430100, China
4
PetroChina Jidong Oilfield Company, Tangshan 063200, China
5
Bohai Drilling Engineering Technology Institute, Tianjin 300457, China
*
Author to whom correspondence should be addressed.
Processes 2024, 12(1), 97; https://doi.org/10.3390/pr12010097
Submission received: 28 November 2023 / Revised: 23 December 2023 / Accepted: 28 December 2023 / Published: 31 December 2023
(This article belongs to the Special Issue Study of Multiphase Flow and Its Application in Petroleum Engineering)

Abstract

:
A pumping-while-tripping method is proposed to mitigate pressure swabs during tripping out in wells with a narrow mud density window and extended reach. In the proposed tripping-out process, the fluid circulation is started by using a special pump from a customized circulation line before tripping is initiated. During the tripping out, drilling fluid is circulated in the wellbore simultaneously while the drilling string is moving. A model to simulate the dynamic pressure changes in this process is developed based on the Navier–Stokes (N-S) equations and a damped free vibration system. The model was initially developed for Herschel–Bulkley (H-B) fluid; however, it can be applied to other fluid models by eliminating the non-existing terms. An analysis was conducted to investigate the effect of tripping velocity and circulation pumping rate on the pressure changes. The results show that pumping-while-tripping is effective in mitigating the pressure swab during tripping out, which is especially useful for extended-reach wells. It can also help to increase tripping out velocity and save tripping time for drilling operations.

1. Introduction

As more and more extended-reach drilling (ERD) is being conducted, the mud density window left for drilling operations becomes much narrower. The pressure surge/swab caused by tripping becomes significant as the measured depth of the wells increases. In drilling operations, even short durations of drilling outside of the safe mud window has led to costly well complications [1,2]. To maintain the wellbore pressure within the drilling window while tripping drill string, running casing, and cementing, the surge/swab pressures must be controlled precisely [3,4,5].
Existing approaches to determine the pressure swab and surge include steady-flow models and dynamic models. The steady-flow model focuses on the enduring equilibrium state at the wellbore bottom over an extended period, while the dynamic model takes into account the system’s instantaneous variations occurring within a brief timeframe. The primary distinction between the two models lies in the temporal scale and the approach adopted for handling pressure fluctuations. Burkhardt, Fontenot and Clark, and Schuh present good representations of the steady-flow model. Burkhardt first proposed a theoretical description of surge pressure and analyzed the influence of different parameters (such as the rheological model of drilling fluid, mud properties, inner-pipe end types, geometry, and tripping velocity) on surge/swab pressure during tripping [5]. Schuh provided an approximate numerical model for power-law fluids [6]. Fontenot and Clark improved the equations proposed by presenting a comprehensive and general approach, including both power-law and Bingham plastic fluids, to calculate bottom hole pressure while tripping and implemented them into a computer program [5,6,7]. Bing et al. studied the steady-state surge/swab pressure of the laminar flow of an H-B fluid [8]. In their models, the drilling fluid is assumed to be incompressible so that the drilling fluid is perfectly displaced by the drill string motion. The pressure losses caused by a moving inner pipe through a concentric annulus are used for predicting swab/surge pressures. An experimental study by Hussain and Sharif shows that the fluctuation pressure decreases with the increase in eccentricity [9]. Srivastav verified this rule later and developed a generic model to simulate the fluctuating pressure reduction caused by eccentricity via a concentric annulus modeling method [10]. As the research continues, more influence factors on surge/swab pressure have been taken into account, such as fluid compressibility, the elasticity of formation and drill string, axial vibration of the string, the effect of temperature on fluid rheology, and inflow of reservoir fluid, which improves the accuracy of model calculation [11,12]. Freddy presented a new steady-state model that can account for fluid and formation compressibility and pipe elasticity in H-B fluid [13]. Comparisons of the model predictions with the measurements showed a satisfactory agreement. Vadim’s model accounts for the mud compressibility, the elasticity of the drill string and the wellbore walls, and the string’s longitudinal vibrations [14], which shows that the amplitude and wave velocity of pressure decrease if elasticity effects are taken into account. Robello provided a coupled model not only considering the effects of fluid inertia, compressibility, wellbore elasticity, axial elasticity of the pipe and temperature-dependent fluid properties but also included a reservoir fluid influx [15], which shows the effect of the downhole pressures owing to mixed fluids under fluid influx results not only in lower wellbore pressures at the bottom of the well but also at the shoe. To investigate the effects of pipe speed (i.e., tripping speed), fluid properties, and borehole geometry on surge and swab pressures under laboratory conditions, Freddy performed a series of tests in varying the tripping speed and developed a model to calculate surge and swab pressures under steady-state flow conditions [16]. Experimental results and model predictions confirm that the trip speed, fluid rheology, annular clearance and pipe eccentricity significantly affect the surge and swab pressures. Tang presented a new model that takes the effects of the velocity of the drilling string on the boundary conditions of surge/swab pressure into consideration [17]. Meng established an elastodynamics model and a poroelastodynamics model, provided a numerical solution by the implicit finite difference method, and verified the analytical solution achieved from the Laplace transform and Crump inversion [18], which filled the gap in calculating dynamic wellbore stress under tripping conditions.
The first dynamic model was developed by Lubinski [19]. The model emphasized the importance of drilling fluid compressibility in swab–surge calculations. Lubinsiki’s dynamic approach is based on Bergeron’s method in the 1930s for solving complicated engineering problems of plane wave propagation. Different from traditional methods, this method is used to find the relationship between pressure and flow rate by imaginary observers moving upstream and downstream with the velocity of propagation. Lal corrected a number of deficiencies in Lubinski’s model and investigated the effects of pipe acceleration as well as the formation elasticity and borehole expansion [20]. Mitchell added the effect of pipe axial elasticity to the dynamic surge analysis [21]. The results of the dynamic models were verified by field data, and the verifications showed that a dynamic model is capable of accurately simulating the swab–surge during drilling. Gjerstadd and Time developed a set of simplified flow equations for H-B fluid in laminar Couette Poiseuille flow for a dynamic surge/swab model, replacing the complex part of the solution with a simpler approximation [22]. Zhang presented a new transient swab/surge model, which considers drill string components, wellbore structure, formation elasticity, pipe elasticity, fluid compressibility, fluid rheology, and the flow between wellbore and formation [23].
Typically, there is no pumping during the tripping operations. For tripping in process, tripping should be stopped approximately every 5000 feet, and “break circulation” needs to be run to fill the empty volume inside the drill string. Moreover, pumping is useful in reducing the pressure swab during tripping out. With fluid circulation during tripping out, the tripping speed can be increased without increasing the pressure swab; thus, significant tripping time can be saved. Furthermore, in choosing the mud weight, the hydra static pressure can be designed closer to the pore pressure without worrying about the pressure swab, which leads to a larger margin in the drilling window [24]. Thus, pumping is required for the tripping out process in back reaming, tight hole conditions and situations where swab pressure is concerned. For tripping out without circulation, there would be a pressure swab; for tripping out with circulation, there would be either a pressure surge or pressure swab depending on the tripping speed and the pumping flow rate [25]. If fluid is circulating during tripping out, the traditional models for surge/swab calculations are no longer applicable, and new methods are needed for the pressure prediction of the situations of pumping and tripping simultaneously [26,27,28,29,30].
This work focuses on the effects of fluid circulation on pressure change during tripping out. The model is developed based on the N-S equations for H-B fluids and can be applied to the other three major fluid models (Bingham plastic, power law, and Newtonian) by eliminating the relevant terms. The H-B model is primarily characterized by three essential parameters: the flow behavior index, yield point, and consistency index. The model can be simplified to the Bingham plastic model by setting the “flow behavior index” to 1 and the “yield point” to 0. The main idea of this approach is to calculate the pressure gradient through the study of the fluid flow velocity profile in the wellbore by moving the drill string. The relationships between flow rate, pressure, and drill string tripping velocity are obtained through N-S equations for H-B fluid. Numerical methods are used in solving the model. The effects of tripping speed, tripping velocity profile, and fluid circulation flow rate on the wellbore pressure during tripping and circulating operations are analyzed.

2. Model

During the process of tripping, it is imperative to account for all potential scenarios within the model, encompassing both the internal flow within the drill string and the circulation within the annulus formed between the wellbore and the moving drill string. The approach involves the incorporation of the H-B model within the framework of the N-S equations. This is followed by the establishment of the N-S equations based on appropriate assumptions and simplified geometrical properties. Subsequently, the application of boundary conditions enables the derivation of an expression delineating frictional loss as a function of both tripping velocity and circulation rate. The absence of an analytical solution for the N-S equations pertaining to the H-B model necessitates the utilization of numerical integration techniques for computation purposes.
To simplify the problem, the following assumptions are made in the formulation of the model.
  • The fluid is incompressible;
  • Only the axial direction flow is considered;
  • The flow is fully developed;
  • No slip condition of the fluid flow at the wall;
  • Drilling string is assumed to be concentric with the wellbore.

2.1. Flow Inside Drill String

The velocity distribution of fluid flow inside the drill string in the context of the H-B model is depicted in Figure 1. This distribution exhibits two plausible scenarios: one where the fluid and the drill string are moving in different directions (Figure 1a), and the other where both the fluid and the drill string move in the same direction (Figure 1b). In the context of drilling applications, the occurrence of scenario (b) is infrequent. This study specifically investigates scenario (a), which manifests during the tripping process involving fluid circulation.
The velocity profile for case (a) has two distinct regions: I. the sheared region flows near the pipe wall, II. the plug zone, where the shear rate is zero and the flow velocity is constant.
The expression of viscosity for the H-B model is:
μ a p p = τ y γ + K γ m 1
where μ a p p is the apparent viscosity, cP, τ y is the yield point, psi, γ is the shear rate, s−1, K is the fluid consistency index, lb-sn/ft2, and m is the fluid behavior index.
The model derivation commences with the N-S equations. Initially, Equation (1) is substituted into the N-S equations in cylindrical coordinates. Subsequently, the equations are simplified based on assumptions, and boundary conditions are applied to determine the fluid flow velocity distribution within the pipe. Finally, integrating the velocity distribution across the cross-sectional area of the drill string yields the pressure gradient as a function of the net flow rate and the drill string’s movement velocity. Detailed steps of the derivation are provided in Appendix A.
After some derivations, the fluid flow velocity profile for region (I) and (II) can be expressed as Equation (2):
v I = A 1 r + A 2 r A 3 1 m d r + C 2
where
A 1 = 1 2 K Δ P Δ L
A 2 = C 1 r K
A 3 = τ y K
where C 1 and C 2 are constants from the integration; they can be calculated by applying the boundary conditions.
The radius of the plug flow region is:
a = 2 τ y Δ P / Δ L
The boundary conditions for the velocity profile are:
At   r = a ,   d v d r = 0
At   r = R p ,   v = V p
where a is the radius of the plug flow region, R p is the inner radius of the pipe, ft, and V p is the pipe movement velocity, ft/s.
The net flow rate inside the pipe can be expressed as:
Q p = 2 π 0 R p v r r d r = 2 π 0 a v I I r r d r + a R p v I r r d r
where v I I is the plug region flow velocity, v I I = v I ( a ) .
The analytical solution for the velocity profile in region I is notably intricate and impractical for real-world applications. Utilizing numerical integration methods becomes necessary to attain the velocity profile. By employing the Simpson’s rule, the integration presented in Equation (2) can be reformulated as:
r 0 r 2 A 1 r + A 2 r A 3 d r h 3 F r 0 + 4 F r 1 + F r 2
where h is the numerical integration step, and
F r = A 1 r + A 2 r A 3
The step-by-step procedure to calculate the pressure gradient by using this approach is shown below:
  • Assuming a pressure gradient, Δ P / Δ L ;
  • Calculate A 1 in Equation (2);
  • Calculate the plug flow region radius by using Equation (6);
  • Calculate C 1 in Equation (2) by using the boundary condition: At r = a , d v d r = 0 ;
  • Assuming the plug flow region velocity v I I ;
  • Integrate the velocity profile from r = a to r = R p by using Equation (10), to obtain the velocity at the inner wall of the pipe V p _ a s s u m p t i o n ;
  • Compare V p with V p _ a s s u m p t i o n , if the difference is larger than tolerance, then go back to step 5 and change the guessed flow region velocity v I I and repeat step 6 and 7. Otherwise, continue to step 8;
  • Calculate the net flow rate in the pipe, Q a s s u m p t i o n , by integrating Equation (9) numerically;
  • Compare Q a s s u m p t i o n with the real flow rate Q r e a l , if Q a s s u m p t i o n > Q r e a l , go back to step 1 and decrease Δ P / Δ L ; if Q a s s u m p t i o n < Q r e a l , go back to step 1 and increase Δ P / Δ L , until convergence is obtained.

2.2. Flow in the Annulus

The fluid flow in the annulus is more complex than inside the drill string. The three potential scenarios depend on the tripping out velocity and the pumping flow rate: (1) Q p u m p i n g > Q d i s p l a c e d , the direction of the average flow velocity is the same as the tripping direction; (2) Q p u m p i n g = Q d i s p l a c e d , the average flow velocity is zero; (3) Q p u m p i n g < Q d i s p l a c e d , the direction of the average flow velocity is opposite with the tripping direction. The flow velocity profiles for different scenarios are shown in Figure 2.
In Figure 2, the drill string is moving out. Assume there is no slip at the wall, so the fluid velocity is the same as the drill string velocity at the wall of the drill string, and the fluid velocity is zero at the wall of the wellbore. In this figure, profile (a) is the case that Q p u m p i n g < Q d i s p l a c e d , profile (b) is the case that Q p u m p i n g = Q d i s p l a c e d , and profile (c) is the case that Q p u m p i n g > Q d i s p l a c e d .
The approach to calculate the velocity profile for H-B fluid in the annulus is similar to the approach used inside the drill string. Because of the complexity of the geometry and the boundary conditions, the concentric annulus is expressed by an equivalent narrow slot.
Start from the N-S equations, after a series of derivations, detailed steps of the derivation are provided in Appendix B. In Figure 2a,b, the velocity profile for each region can be expressed as:
Region I:
v ~ I = π 1 y ~ 1 y ~ b y ~ 1 b ;               0 y ~ y ~ 1
Region III:
v ~ I I I = 1 π 1 1 y ~ 1 b y ~ y ~ 2 b ;               y ~ 2 y ~ 1
where
π I = m m + 1 H V p Δ P Δ L K H 1 m
and
b = n + 1 n ;   v ~ 1 = v 1 V p   ;   v ~ 2 = v 2 V p   ;   y ~ 1 = y 1 H   ;    y ~ 2 = y 2 H   ;
In Figure 2c, the velocity profile for each region can be expressed as:
Region I:
v ~ I = π 1 y ~ 1 b y ~ 1 y ~ b ;               0 y ~ y ~ 1
Region III:
v ~ I I I = 1 + π 1 1 y ~ 1 b y ~ y ~ 2 b ;               y ~ 2 y ~ 1
The thickness of the plug region can be expressed as:
y ~ 2 y ~ 1 = π 2 = 2 τ 0 / H Δ P / Δ L
The fluid flow velocities at y = y 1 and y = y 2 should be the same. By combining Equations (15)–(17), the following can be obtained:
1 y ~ 1 π 2 b y ~ 1 b 1 π 1 = 0
Similarly, In Figure 2a,b, combine Equations (9), (10), and (18) to obtain
1 y ~ 1 π 2 b y ~ 1 b + 1 π 1 = 0
The total flow rate is the sum of the flow rate in each region:
Q t = 2 π W 0 H v y d y = 2 π W 0 y 1 v I y d y + y 1 y 2 v I I y d y + y 2 H v I I I y d y
The total flow rate is the sum of the pumping flow rate, Q p u m p i n g , and the fluid volume displaced by drill string in unit time, Q d i s p l a c e d . The Q d i s p l a c e d is positive for tripping in and negative for tripping out.
The algorithm for this method is shown as below:
  • Calculate the total flow rate in the annulus, Q t , based on the pumping flow rate and the tripping velocity, then, obtain the direction of the average flow rate;
  • Choose the velocity profile from Figure 2, and determine the equations that can represent the velocity profile;
  • Guess a pressure gradient based on the direction of the average flow velocity;
  • Transfer all the parameters dimensionless by using Equations (12) and (13);
  • Calculate the thickness of the plug region by using Equation (18);
  • Calculate y ~ 1 by solving Equations (19) or (20) (these two equations are for different velocity profiles) numerically (Newton’s method is a good choice);
  • Calculate the velocity profile for each region (Equations (9) and (10) or Equations (15) and (16));
  • Integrate the velocity profile and calculate the assumptive total flow rate in the annulus, Q t _ a s s u m p t i o n ;
  • Compare the real flow rate, Q t , and assumptive flow rate, Q t _ a s s u m p t i o n , if the difference is larger than tolerance, go back to step 3 and change the pressure gradient, and then repeat steps 3–8. Otherwise, the system becomes converged and outputs the pressure gradient.

2.3. Non-Converging Region

It needs to be mentioned that for cases when the direction of the average velocity is the same as the direction of pipe movement, there is a possibility that the velocity profile does not have a plug flow region even if the fluid has yield power law properties. The fluid flow is mainly caused by the movement of the drill pipe rather than pressure gradient. In such circumstances, if the velocity profile depicted in Figure 2 is continuously applied, convergence might not be achieved during the iterative process. Therefore, alternative functions are required to handle this specific region.
The first step is to detect the non-converging region and find the boundaries for this region. The non-converging happens when the pressure gradient is close to zero. Therefore, the way to detect the non-converging region is to start from a relatively large pressure gradient to calculate the according flow rate. Decrease the pressure gradient gradually until the function cannot find a converging flow rate, and then, record the first non-converging value as the boundary.
The pressure gradient can be negative or positive in the wellbore. Therefore, when detecting non-convergent regions, both positive and negative aspects need to be considered to identify two boundaries.
Once the boundaries are detected, a square function is plugged into the non-converging region to represent the relationship between the net flow rate and the pressure gradient.

2.4. Turbulent Model

The calculation of pressure loss in the turbulent zone during tripping is predicated upon the regular pressure loss characteristics exhibited in turbulent flow. Addressing the impact of drill string movement on pressure loss within the turbulent regime involves considering the relative flow velocity between the fluid flow and the moving drill string. For example, the absolute velocity of the fluid flow is U, and the drill string is moving at velocity V; thus, the relative velocity between the fluid flow and drill string is U + V. The equation to calculate the friction loss inside a drill string for a stationary wall is:
D P D L = ρ f 2 g d Q A Q A
where d is the drill pipe inner diameter.
If the drill pipe is moving at velocity V , the following equation can be used:
D P D L = ρ f 2 g d Q A + V Q A + V
The friction loss in annulus with both walls stationary can be expressed by:
D P D L = ρ 2 g d o 2 d i 2 d i f i Q A Q A + d o f o Q A Q A
where d i is the inner diameter of the annulus, and d o is the outer diameter.
In the event of drill string movement, the above equation can be modified for the annulus with a moving inner wall and a stationary outer wall, and the frictional losses can be expressed as:
D P D L = ρ 2 g d o 2 d i 2 d i f i Q A + V Q A + V + d o f o Q A Q A

2.5. Effect of the Inertia

The surge/swab component of inertial pressure fluctuations arises due to the tendency of the mud column to resist changes in motion. It can be expressed by the following:
For the closed-end strings:
Δ P Δ L = ρ a p D p 2 g D w 2 D p 2
For the open-end strings:
Δ P Δ L = ρ a p D p 2 D i 2 g D w 2 D p 2 + D i 2
where ρ is the density of the drilling fluid, ppg, a p is the acceleration of the string, ft/s², D p is the outer diameter of the string, ft, D i is the inner diameter of the string, ft, and D w is the diameter of the well, ft.

2.6. Effect of the Gel

The formulations to calculate the pressure required to break the mud gel and initiate fluid circulation have been given by Melrose [31].
Inside the drill string:
Δ P Δ L = 4 ζ D i
In the annulus:
Δ P Δ L = 4 ζ D w D p
where ζ is the gel of drilling fluid, lb/100 ft². Once the gel is broken, the effect of the gel would be insignificant; the pressure gradient would be mainly determined by viscous and inertia.

2.7. Model Validation

Crespo conducted an experiment to determine the fluctuating pressure of H-B fluid in a concentric annulus. The fundamental parameters used in the experiment are presented in Table 1, and Figure 3 illustrates the comparison between the results of our model and the experimental finding [17].
Based on the comparative results depicted in Figure 3, it is evident that the model’s calculated outcomes align well with the experimental results in the non-pumping condition. However, after initiating pumping, the model’s calculated results are comparatively lower than the experimental findings conducted by Crespo. This discrepancy is primarily attributed to the reduction in surge pressure caused by the pumping operation.

3. Sensitivity Analysis

The analysis below was run for a practical case in which there is a drill bit at the end of the drill string. If there is no pumping, it is treated as a closed-end scenario since the nozzles on the drill bit are extremely small, and the flow through the nozzles can be neglected. If there is pumping, the fluid will flow through the nozzle and then flow back to the surface through the annulus. All the scenarios described below are conducted on a vertical well of 5000 feet depth with a diameter of 8.5 inches by 4.5 inches. The fluid parameters and base geometry parameters are shown in Table 2 and Table 3.

3.1. Effect of Tripping Velocity

Equivalent circulating density (ECD) is a parameter that evaluates the impact of fluid circulation on the bottom hole static pressure during drilling operations. The effects of tripping out speed on ECD at different circulation flow rates for the well given in Table 3 are shown in Figure 4. The ECD decreases as the tripping velocity increases, which is caused by the pressure swab. At lower tripping speeds, there is a relatively higher ECD, whereas at higher speeds, the ECD decreases due to the influence of pressure swabbing. This relationship requires precise control during drilling operations to ensure that the ECD remains within a safe range, simultaneously enhancing drilling efficiency. When there is no fluid circulation (0 gpm) during tripping, the ECD is always lower than the fluid density at static conditions, and there is a significant decrease in ECD as the tripping velocity increases to 200 feet/min. When the pumping flow rate is 160 gpm, at a tripping velocity lower than 100 ft/min, the ECD is higher than the equivalent static density (ESD) because of the pressure surge caused by the over-pumping; at a tripping speed larger than 100 ft/min, the pumping flow rate is not enough to balance the pressure swab caused by tripping, so the ECD is lower than the ESD. The term ESD refers to the static density of the fluid in the wellbore under specified conditions, excluding the effects of fluid circulation or drilling operations.
When the pumping flow rate reaches 320 gpm, it consistently results in an ECD exceeding the ESD. As the rate of tripping the drill string out of the borehole increases, the amplitude of pressure surges gradually diminishes. It is readily discernible that, in this particular scenario, the flow rate of 320 gpm exceeds the optimal requirement for the tripping-out procedure. In typical field operations, the conventional tripping velocity stands at approximately 90 feet per minute. Hence, in situations where the drilling safety density window is narrower than the suction pressure without fluid circulation, it is recommended to reduce the fluid pumping flow rate to 160 gpm during the tripping-out phase.

3.2. Effect of Pumping Flow Rate

The effect of the pumping flow rate on ECD is shown in Figure 5. The augmentation in pumping flow rate results in an escalated fluid velocity and pressure at the wellbore bottom, potentially leading to an increase in ECD. Under a high pumping flow rate, a pressure swab may occur due to the phenomenon of rapid fluid accumulation at the wellbore bottom, causing transient elevations in ECD. It can be observed from this figure that if the pumping flow rate is very low (which is not practical in field operations), ECD remains nearly unchanged. After a certain point, the ECD increases with increasing the pumping flow rate, regardless of the tripping velocity. When increasing the pump flow rate to make the ECD value equal to the ESD value, this pumping flow rate is referred to as the “balanced pumping flow rate”. The balanced pumping flow rate increases as the tripping velocity increases. Once the circulation flow rate is larger than the balanced pumping flow rate, there would be a pressure surge during tripping out, and the ECD is larger than the ESD. The ECD keeps increasing as the pumping flow rate increases. In the cases of over-pumping, the pressure surge is more severe for a low tripping speed than a high tripping speed. At high tripping speed, a little bit of over-pumping has a slight effect on the ECD.

3.3. Effect of Inertia

The effect of inertia on the pressure change during tripping is not very significant [32]. The results are shown in Figure 6 for open-end cases and Figure 7 for closed-end cases. Normally, the acceleration is no more than 0.5 ft/s2, which means the drill string velocity reaches 90 ft/min in 3 s. In this situation, the pressure change caused by the drill string acceleration is 5 psi for open-end cases and about 20 psi for the closed-end cases for given examples in Table 2 and Table 3, which is insignificant for a 5000 ft depth well.
In scenarios involving the continuous circulation of drilling fluid during the process of tripping the drill string out of the borehole, the fluid is inherently in motion, thus mitigating the influence of the acceleration of the drill string on the fluid. In succinct summation, it can be ascertained that the impact of inertial forces on the alteration of pressure dynamics during the tripping-out procedure remains negligible.

3.4. Effect of the Fluid Rheology

Alterations in fluid rheological parameters, such as yield point, fluid behavior index, and fluid consistency index, have the potential to modify the flow characteristics of the fluid within the wellbore, consequently impacting the ECD. The effects of fluid rheology parameters on the ECD during tripping out with pumping are shown in Figure 8, Figure 9 and Figure 10.
Figure 8 shows the effect of yield point on ECD at different tripping velocities. The pressure surge/swab is more severe for fluids with high yield points than fluids with low yield points. For fluids with high yield points, it takes more pressure to move it in the annulus, which leads to a lower ECD for inefficient pumping and a higher ECD for over-pumping. Therefore, accurately controlling the pumping flow rate is required for a high yield point fluid.
The effects of the fluid consistency index and fluid behavior index are shown in Figure 9 and Figure 10, separately. Increasing either of these two indexes leads to a higher effective viscosity, which causes a more severe pressure surge/swab.

4. Tripping out Case

In this section, the preceding research findings are practically applied to drilling operations to demonstrate the impact of pump initiation on bottom-hole pressure during the drilling process. The well information pertinent to the case study is visually represented in Figure 11. In this instance, the kick-off depth is approximately 10,269 ft, where the inclination angle reaches 86.1° at 13,794 ft and, subsequently, remains constant up to 20,080 ft. The azimuth angle remains unchanged at 10.97°.
After the drilling finished, the 4 ½ inch drill strings were pulled out of the hole at the tripping velocity profile shown in Figure 12.
The bottom hole pressure during tripping out of the first stands for this case is shown in Figure 13. At time zero, the drill bit is at the bottom hole. The blue plot is the bottom hole pressure during tripping without pumping, and the red plot is the bottom hole pressure with pumping at 150 gpm.
Without fluid circulation, there is a significant pressure swab after the tripping starts. After the tripping velocity becomes steady, the bottom hole pressure also becomes steady (in this case, a steady model was used and fluid compressibility was not considered). At the end of the tripping of one stand, the tripping velocity gradually decreases to zero, and the bottom hole pressure returns back to the hydra static pressure. It keeps constant at the hydra static pressure during the time to disconnect the drill string until the tripping of the second stand starts.
The pressure changes for tripping with fluid circulation (red plot in Figure 13) are different with no fluid circulation. In this case, the operation procedure is to turn on the fluid circulation before the tripping starts and turn off the fluid circulation after the tripping of one stand finishes. Therefore, there are pressure surges at the beginning and end of the tripping, which are caused by the fluid circulation. In fact, in this particular scenario, there is no need for concern regarding the pressure surge: the fluid circulation flow rate during tripping is less than the fluid circulation rate during drilling, which means the pressure surge here is much less than the pressure increase caused by the fluid circulation during drilling.
The pressure swab during the tripping out process with fluid circulation is significantly lower compared to without fluid circulation. When designing the well, it is possible to bring the hydrostatic pressure closer to the pore pressure. Consequently, a larger safety margin can be preserved within the drilling window, enabling further drilling in extended-reach drilling operations. Furthermore, with fluid circulation during tripping, the tripping speed can be increased without increasing the risk of fluid influx, which means the tripping time can be saved.
In Figure 14, the change in bottom hole pressure during the whole tripping out process under different operation conditions was plotted. The horizontal axis of this figure is the position of the drill bit during tripping. Assuming that the drill bit was positioned at the bottom hole at the start of the tripping operation, as the tripping goes on, the drill bit moves closer to the wellhead, and the total pressure swab caused by tripping decreases as the length of the drill string in the well decreases. From Figure 14, the following remarks can be concluded: (1) the pressure swab is significant only in deep wells; (2) higher tripping speed leads to a more severe pressure swab if other operation parameters are the same; (3) fluid circulation during tripping out can reduce the pressure swab; (4) insufficient pumping during tripping out at high speed can still lead to significant pressure swab (150 gpm, 200 ft/min); (5) over-pumping may lead to a pressure surge in the wellbore (260 gpm, 100 ft/min); (6) the severe pressure swab is mainly caused by the tight hole conditions (the pressure changing ratio is the highest in an 83/4 inch open hole).
Due to the relations between the pressure swab, drill bit position in the well, and the well geometries summarized above, tripping operational parameters can be adjusted during various phases of the tripping process to optimize the cost-effective retrieval of the drill string. For the case shown above, two tripping plans are presented herein (Figure 15 and Figure 16).
In Figure 15, it is assumed that the drill bit was positioned at the bottom hole at the start of the tripping operation, and the draw works have enough power to pull out the drill string at the given velocity. This plan is designed to trip out the drill strings at a constant high tripping velocity (200 ft/min). As shown in Figure 15, the whole tripping process has three periods: high pumping flow rate, low pumping flow rate, and no pumping. In this case, the pore pressure (green curve in figure) is about 890 psi lower than the hydrostatic pressure (red curve in figure). To keep the bottom hole pressure above the pore pressure during tripping, the pressure swab must be controlled within 890 psi. At the beginning, the clearance between the drill string and the wall of the well bore is narrow, and the length of the drill string inside the wellbore is very long. To avoid a severe pressure swab, the drilling fluid is circulated at 400 gpm. When the drill bit reaches the liner section of the well, the clearance becomes larger and the drill string length is shorter, the pumping flow rate can be lowered to 150 gpm. At last, the pumping is no longer needed at the last stage of the tripping.
The second tripping plan is shown in Figure 16. In this case, it is assumed that the winch lacks sufficient horsepower to pull the drill string out of the wellbore at a rate of 200 ft/min from the bottom. Therefore, at the beginning, the tripping velocity is 100 ft/min, and the pumping flow rate is 150 gpm. Increase the tripping velocity to 200 ft/min when the load of the drill string acting on the draw works is low enough. At last, stop the pumping when the drill bit position is very shallow.

5. Conclusions

From the studies above, the following conclusions can be obtained:
  • In deep well or tight well conditions, there is a significant pressure swab during tripping out without fluid circulation. Pumping fluid properly during tripping out is able to decrease or diminish the pressure swab;
  • With fluid circulation during tripping out, the tripping speed can be increased without increasing the pressure swab; thus, significant time can be saved during tripping;
  • In choosing the mud weight, the hydra static pressure can be designed to be closer to the pore pressure without worrying about the pressure swab by circulating fluid flow during tripping out, which leads to a larger margin in the drilling window for drilling deeper or saves some casing work;
  • By using the modeling results, tripping out operation parameters can be optimized based on the formation properties, well properties, and rig facilities.
Although this study has achieved some results, it is not without limitations. In the next phase of research, we intend to validate our findings through additional instances and consider other potential influencing factors. For instance, we plan to explore the impact of temperature on the ECD.

Author Contributions

C.S. conceived and designed the program, and F.Z. supervised the study and debugged the program. Y.L. and X.L. assisted in the development of the algorithm. Y.T. analyzed the computed results and wrote the paper. 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 (No. 52374003) (Z.F.), Project of Scientific Research Program of Hubei Provincial Department of Education (Grant No: T2021004) (Z.F.), Hubei Province key research and development project: Research on digital twin technology and equipment for intelligent drilling of long horizontal Wells of shale gas in western Hubei Province (No: 2023BCB111) (Z.F.).

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgments

The authors also gratefully acknowledge the financial support of the CNPC Bohai Drilling Corporation.

Conflicts of Interest

Author Yaoyao Tang was employed by the PetroChina Jidong Oilfield Company. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

The velocity profiles for H-B fluid flow under different situations are derived from the N-S equations. In this case, only the axial flow of fluid needs to be considered. The original N-S equation for the z direction under cylindrical coordinates is:
ρ V z t + V r V z r + V θ r V z θ + V z V z z = P z 1 r r r τ r z + 1 r θ τ θ z + θ τ z z + ρ g z
To simplify the equation, the following assumptions are made:
  • The fluid is incompressible;
  • The flow is in steady state and only happens in the axial direction, which means V z t = 0 ,   V r = 0 ,   V θ = 0 ;
  • The flow regime is laminar;
  • The flow is fully developed, so V z z = 0 .
Therefore, all the terms on the left side of the equation are neglected. Here, the influence of gravity is not currently taken into consideration, so Equation (A1) can be simplified to:
P z = 1 r r r τ r z + 1 r θ τ θ z + θ τ z z
From the definition of shear stress:
τ θ z = μ a p p 1 r V z θ + V θ z τ r z = μ a p p V z r + V r z τ z z = μ a p p 2 V z z
where
μ a p p = τ y γ + K γ m 1
From the assumptions, it is easy to see that τ θ z and τ z z are zero. Substitute Equation (A4) into τ r z :
τ r z = τ y + K d V z d r m
Substitute Equation (A4) into Equation (A2):
P z = 1 r r r τ y + K d V z d r m
Integrate Equation (A6) twice intothe expression for V z :
V z = A 1 r + A 2 r A 3 1 m d r + C 2
where
A 1 = 1 2 K Δ P Δ L A 2 = C 1 r K A 3 = τ y K

Appendix B

The annular velocity profile (Figure 2) in the wellbore during tripping operations has three distinct regions. (1) Outer sheared region (Region I) within the boundary limits 0 y y 1 ; (2) plug zone (Region II) within the boundary limits y 1 y y 2 and, (3) inner sheared region (Region I) within the boundary limits y 2 y H . For the sheared regions I and II, the momentum balance is expressed as:
Region I:
τ 1 = Δ P Δ L y + τ w
Region III:
τ 3 = Δ P Δ L y Δ P Δ L y 2 + τ 0
where τ w is wall shear stress, which is the average of the wall shear stresses acting on the inner and outer walls. For the H-B fluid, the local shear stresses in Region I and III are related to the local shear rates using the constitutive equation as:
τ 1 = τ 0 + k d v d y n
and
τ 3 = τ 0 + k d v d y n
Respectively, where k and n are the consistency index and fluid behavior index. Combining Equations (A9)–(A12), the dimensionless velocity distributions can be obtained.

References

  1. Yin, Q.; Yang, J.; Li, Z. A Field Case Study of Managed Pressure Drilling in Offshore Ultra High-Pressure High-Temperature Exploration Well in the South China Sea. SPE Drill. Complet. 2020, 35, 503–524. [Google Scholar] [CrossRef]
  2. Nauduri, S.; Parker, M.; Shimi, A.; Guefack, G. A Review of Conventionally Undrillable Wells with Ultra-Narrow Pressure Windows Successfully Drilled Utilizing MPD Over the Last Five Years in USA and Canada. In Proceedings of the SPE/IADC Managed Pressure Drilling & Underbalanced Operations Conference & Exhibition, Kuala Lumpur, Malaysia, 27–28 September 2022. [Google Scholar]
  3. Lou, W.; Wang, Z.; Zhang, J. Flow regime evolution mechanism and falling flux prediction model for bypass injection of viscous liquid in vertical T-junction. Chem. Eng. J. 2023, 476, 146601. [Google Scholar] [CrossRef]
  4. Carpenter, C. Drilling-System Optimization Leads to Expedited Wells Delivery. J Pet. Technol. 2022, 74, 69–71. [Google Scholar] [CrossRef]
  5. Damski; Carlos; Sameh, E. The Hard Case for Whole Drilling Optimization. In Proceedings of the SPE/IADC Middle East Drilling Technology Conference and Exhibition, Abu Dhabi, United Arab Emirates, 25–27 May 2021. [Google Scholar]
  6. Burkhardt, J. Wellbore Pressure Surges Produced by Pipe Movement. J. Pet. Technol. 1961, 13, 595–605. [Google Scholar] [CrossRef]
  7. Schuh, F.J. Computer Makes Surge-Pressure Calculations Useful. Oil Gas J. 1964, 62, 96–104. [Google Scholar]
  8. Fontenot, J.E.; Clark, R.K. An Improved Method for Calculating Swab and Surge Pressures and Circulating Pressures in a Drilling Well. Soc. Pet. Eng. J. 1974, 14, 451–462. [Google Scholar] [CrossRef]
  9. Bing, Z.; Kaiji, Z.; Qiji, Y. Equations Help Calculate Surge and Swab Pressures in Inclined Wells. Oil Gas J. 1995, 93, 204–215. [Google Scholar]
  10. Hussain, Q.E.; Sharif, M.A.R. Viscoplastic Fluid Flow in Irregular Eccentric Annuli Due to Axial Motion of the Inner Pipe. Can. J. Chem. Eng. 1997, 75, 1038–1045. [Google Scholar] [CrossRef]
  11. Srivastav, R.; Enfis, M.S.; Crespo, F.E. Surge and Swab Pressures in Horizontal and Inclined Wells. In Proceedings of the SPE Latin America and Caribbean Petroleum Engineering Conference, Mexico City, Mexico, 16–18 April 2012. [Google Scholar]
  12. Shwetank, K.; Syahrir, R.; Scott, C. Experimental evaluation of surge/swab pressure in varying annular eccentricities using non-Newtonian fluid under Couette-Poiseuille flow for drilling applications. J. Pet. Sci. Eng. 2021, 206, 108982. [Google Scholar]
  13. Erge, O.; Akin, S. Numerical and Surrogate Modeling of Surge and Swab Pressure in Concentric and Eccentric Annuli. SSRN. 2023. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=4416311 (accessed on 8 September 2023).
  14. Freddy, C.; Ramadan, A.; Arild, S. Surge and Swab Pressure Predictions for Yield-Power-Law Drilling Fluids. In Proceedings of the SPE Latin American and Caribbean Petroleum Engineering Conference, Lima, Peru, 1–3 December 2010. [Google Scholar]
  15. Vadim, T.; Olga, B.; Valery, L. Development of Model and Software for Simulation of Surge-Swab Process at Drilling. In Proceedings of the SPE Russian Petroleum Technology Conference and Exhibition, Moscow, Russia, 24–26 October 2016. [Google Scholar]
  16. Robello, S.; Randy, L. Swab and Surge Pressures with Reservoir Fluid Influx Condition during MPD. In Proceedings of the IADC/SPE Drilling Conference and Exhibition, Fort Worth, TX, USA, 1–3 March 2016. [Google Scholar]
  17. Freddy, C.; Ramadan, A. A simplified surge and swab pressure model for yield power law fluids. J. Pet. Sci. Eng. 2013, 101, 12–20. [Google Scholar]
  18. Tang, M.; Xiong, J.; He, S. A new model for computing surge/swab pressure in horizontal wells and analysis of influencing factors. J. Nat. Gas Sci. Eng. 2014, 19, 337–343. [Google Scholar] [CrossRef]
  19. Meng, M.; Zahra, Z.; Stefan, M. Dynamic stress distribution around the wellbore influenced by surge/swab pressure. J. Pet. Sci. Eng. 2019, 172, 1077–1091. [Google Scholar] [CrossRef]
  20. Lubinski, A.; Hsu, F.H.; Nolte, K.G. Transient Pressure Surges Due to Pipe Movement in an Oil Well. Oil Gas Sci. Technol. 1977, 32, 307–348. [Google Scholar] [CrossRef]
  21. Lal, M. Surge and Swab Modeling for Dynamic Pressures and Safe Trip Velocities. In Proceedings of the IADC/SPE Drilling Conference, New Orleans, LA, USA, 20 February 1983. [Google Scholar]
  22. Mitchell, R.F. Dynamic Surge/Swab Pressure Predictions. SPE Drill. Eng. 1988, 3, 325–333. [Google Scholar] [CrossRef]
  23. Gjerstadd, K.; Time, W.R. Simplified Explicit Flow Equations for Herschel-Bulkley Fluids in Couette-Poiseuille Flow-For Real-Time Surge and Swab Modeling in Drilling. SPE J. 2015, 20, 610–627. [Google Scholar] [CrossRef]
  24. Zhang, F.; Kang, Y.; Wang, Z. Transient Coupling of Swab/Surge Pressure and In-Situ Stress for Wellbore-Stability Evaluation During Tripping. SPE J. 2018, 23, 1019–1038. [Google Scholar] [CrossRef]
  25. Xie, Y.; Min, O.; Zhao, H. Wellbore pressure management of offshore wells during casing running in narrow safety pressure window formations. J. Pet. Sci. Eng. 2021, 199, 108284. [Google Scholar] [CrossRef]
  26. Bacon, R. Optimised Tripping and Drilling Performance Through Drilling Control System Automation to Control Axial Velocity. In Proceedings of the International Petroleum Technology Conference, Bangkok, Thailand, 1–3 March 2023. [Google Scholar]
  27. Tengesdal, N.K.; Holden, C. Identification and optimal control for surge and swab pressure reduction while performing offshore drilling operations. Model. Identif. Control. 2020, 41, 165–184. [Google Scholar] [CrossRef]
  28. Van Thang, N.; Thang Vinh, P.; Rogachev, M. A comprehensive method for determining the dewaxing interval period in gas lift wells. J. Pet. Explor. Prod. Technol. 2023, 13, 1163–1179. [Google Scholar]
  29. Islamov, S.R.; Bondarenko, A.V.; Gabibov, A.F. Polymer compositions for well killing operation in fractured reservoirs. In Advances in Raw Material Industries for Sustainable Development Goals; Taylor & Francis: London, UK, 2021; pp. 343–351. [Google Scholar]
  30. Islamov, S.R.; Bondarenko, A.V.; Mardashov, D.V. A selection of emulsifiers for preparation of invert emulsion drilling fluids. In Topical Issues of Rational Use of Natural Resources; Taylor & Francis: London, UK, 2020; pp. 487–494. [Google Scholar]
  31. Cioklosz, R.; Gunn, P.; Palacios, J. Calibrating ECD, Surge, and Swab Models Using Downhole Pressure Measurements to Avoid Wellbore Instability and Allow Safe Tripping. In Proceedings of the SPE/IADC Managed Pressure Drilling & Underbalanced Operations Conference & Exhibition, Kuala Lumpur, Malaysia, 27–28 September 2022. [Google Scholar]
  32. Thi, T.L.; Thorsten, K. Effect of inertia force on the interface stability of a tangential-velocity discontinuity in porous media. Int. J. Multiph. Flow 2023, 169, 104612. [Google Scholar]
Figure 1. Velocity profile inside the drill string for H-B fluid. (a) the fluid and the drill string are moving in different directions. (b) both the fluid and the drill string move in the same direction. I the sheared region, II the plug zone.
Figure 1. Velocity profile inside the drill string for H-B fluid. (a) the fluid and the drill string are moving in different directions. (b) both the fluid and the drill string move in the same direction. I the sheared region, II the plug zone.
Processes 12 00097 g001
Figure 2. Velocity profile in the annulus for H-B fluid during tripping out. (a) the direction of the average flow velocity is opposite with the tripping direction. (b) the average flow velocity is zero. (c) the direction of the average flow velocity is the same as the tripping direction.
Figure 2. Velocity profile in the annulus for H-B fluid during tripping out. (a) the direction of the average flow velocity is opposite with the tripping direction. (b) the average flow velocity is zero. (c) the direction of the average flow velocity is the same as the tripping direction.
Processes 12 00097 g002
Figure 3. Comparison between calculation results of the presented model and experimental results in reference [17].
Figure 3. Comparison between calculation results of the presented model and experimental results in reference [17].
Processes 12 00097 g003
Figure 4. Effect of tripping speed on ECD.
Figure 4. Effect of tripping speed on ECD.
Processes 12 00097 g004
Figure 5. Effect of pumping flow rate on ECD.
Figure 5. Effect of pumping flow rate on ECD.
Processes 12 00097 g005
Figure 6. Effect of inertia for open-end cases.
Figure 6. Effect of inertia for open-end cases.
Processes 12 00097 g006
Figure 7. Effect of inertia for closed-end cases.
Figure 7. Effect of inertia for closed-end cases.
Processes 12 00097 g007
Figure 8. Effect of yield point on ECD during tripping out.
Figure 8. Effect of yield point on ECD during tripping out.
Processes 12 00097 g008
Figure 9. Effect of fluid consistency index on ECD during tripping out.
Figure 9. Effect of fluid consistency index on ECD during tripping out.
Processes 12 00097 g009
Figure 10. Effect of fluid behavior index on ECD during tripping out.
Figure 10. Effect of fluid behavior index on ECD during tripping out.
Processes 12 00097 g010
Figure 11. Well schematic.
Figure 11. Well schematic.
Processes 12 00097 g011
Figure 12. Tripping velocity profile.
Figure 12. Tripping velocity profile.
Processes 12 00097 g012
Figure 13. Bottom hole pressure during drilling.
Figure 13. Bottom hole pressure during drilling.
Processes 12 00097 g013
Figure 14. Bottom hole pressure drilling process.
Figure 14. Bottom hole pressure drilling process.
Processes 12 00097 g014
Figure 15. BHP for tripping out at 200 ft/min.
Figure 15. BHP for tripping out at 200 ft/min.
Processes 12 00097 g015
Figure 16. Optimized tripping out plan.
Figure 16. Optimized tripping out plan.
Processes 12 00097 g016
Table 1. Rheology of test fluids.
Table 1. Rheology of test fluids.
Test FluidsK (lbf/100 ft2 × sn)n τ 0 (lbf/100 ft2)Pipe OD (Inch)Well (Inch)
0.44% xanthan gum0.00750.520.0720.4170.893
Pipe OD means pipe outer diameter (OD).
Table 2. Fluid rheology parameters.
Table 2. Fluid rheology parameters.
Tal_y (lbf/100 ft2)K (lbf/100 ft2 × sn)mDensity (ppg)Gel Strength (lbf/100ft2)
10 s10 min30 min
16.890.70912.712.7273350
Table 3. Geometry parameters.
Table 3. Geometry parameters.
Pipe ID (Inch)Pipe OD (Inch)Well (Inch)Length (Feet)
3.8264.58.55000
Pipe ID means pipe inner diameter (OD).
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

Sheng, C.; Zhang, F.; Tang, Y.; Li, Y.; Liu, X. Wellbore Pressure Modeling for Pumping and Tripping Simultaneously to Avoid Severe Pressure Swab. Processes 2024, 12, 97. https://doi.org/10.3390/pr12010097

AMA Style

Sheng C, Zhang F, Tang Y, Li Y, Liu X. Wellbore Pressure Modeling for Pumping and Tripping Simultaneously to Avoid Severe Pressure Swab. Processes. 2024; 12(1):97. https://doi.org/10.3390/pr12010097

Chicago/Turabian Style

Sheng, Cancheng, Feifei Zhang, Yaoyao Tang, Yafeng Li, and Xuesong Liu. 2024. "Wellbore Pressure Modeling for Pumping and Tripping Simultaneously to Avoid Severe Pressure Swab" Processes 12, no. 1: 97. https://doi.org/10.3390/pr12010097

APA Style

Sheng, C., Zhang, F., Tang, Y., Li, Y., & Liu, X. (2024). Wellbore Pressure Modeling for Pumping and Tripping Simultaneously to Avoid Severe Pressure Swab. Processes, 12(1), 97. https://doi.org/10.3390/pr12010097

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