Next Article in Journal
Multi-Aspect SAR Target Recognition Based on Non-Local and Contrastive Learning
Next Article in Special Issue
Numerical Simulation of Hole Cleaning of a Horizontal Wellbore Model with Breakout Enlargement Section
Previous Article in Journal
On Recovering Sturm–Liouville-Type Operators with Global Delay on Graphs from Two Spectra
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating the Critical Velocity of the Incipient Motion of Particles on the Cuttings Bed Surface: An Experimental and Theoretical Analysis

1
Guangzhou Marine Geological Survey, China Geological Survey, Guangzhou 510075, China
2
National Engineering Research Center of Gas Hydrate Exploration and Development, Guangzhou 510075, China
3
College of Petroleum Engineering, Northeast Petroleum University, Daqing 163318, China
4
Department of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, NF A1B 3X5, Canada
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2689; https://doi.org/10.3390/math11122689
Submission received: 4 May 2023 / Revised: 4 June 2023 / Accepted: 10 June 2023 / Published: 13 June 2023

Abstract

:
During the drilling of highly deviated and horizontal wells, a pump shutdown causes drill cuttings to settle and form a cuttings bed in the annulus. This study investigated the incipient motion law of the particles on the cuttings bed surface when the drilling fluid starts circulating again. This work could help field engineers to determine a reasonable incipient pump displacement to improve hole-cleaning efficiency. In this study, the effects of the well inclination angle, cuttings size, and different cuttings densities on the critical velocity of particle motion are analyzed experimentally, using a large-scale flow loop. Next, based on a stress analysis of the particles on the cutting bed surface and on the boundary layer flow around the particles, an analytical formula for the surface shear force of the drilling fluid on particles is derived and a critical velocity model for incipient motion is established. Verification is then carried out and combined with the experimental results. This study has important implications for the design of drilling operations and for the management of cuttings transport in oil and gas wells. It can guide the setting and prediction of pump discharge to improve hole-cleaning efficiency.

1. Introduction

A significant amount of drill cuttings are generated during the drilling of extended-reach and horizontal wells. These cuttings are removed from the wellbore by circulating drilling fluid to ensure that the wellbore is properly cleaned [1,2,3]. However, during the construction process, circulation must sometimes be stopped, for example, when changing a dull bit or making a connection, and then drill cuttings are deposited to form a cuttings bed. As the height of the cuttings bed increases, the drill pipe becomes buried. When circulation is restarted, the cutting particles on the cuttings bed surface are mobilized and transported to the wellhead by the drilling fluid. The inadequate removal of the cuttings bed can lead to hole cleaning problems such as high rotary torque, a stuck pipe, formation breakdown, a slow rate of penetration, and a loss of circulation [4,5,6,7]. Therefore, studying the incipient motion law of cuttings beds after the drilling fluid is recycled can guide field engineering practices in determining a reasonable pump displacement to improve hole-cleaning efficiency.
Currently, two types of theoretical models are used to describe cuttings transportation. One is the solid–liquid stratification model, which is based on theoretical assumptions and experimental observations, such as the two-layer model and the three-layer model [8,9,10]. This type of model can be used to calculate parameters such as the cuttings bed height and transportation efficiency. The other type of model calculates the critical velocity required for drilling fluid to carry and transport drill cuttings [11,12,13]. In order to inhibit the formation of a cuttings bed and to keep the wellbore’s hole clean, the drilling fluid velocity needs to be greater than the critical velocity.
In terms of the stratification model, the two-layer model was initially established. However, the disadvantage of the two-layer model was that it did not consider the diffusion of particles from the bed surface to the liquid suspension layer. In their experiment, Nguyen and Rahman [14] observed that, at high flow rates, a mobile diffusion layer forms between the cuttings bed and the flow area of the drilling fluid. Therefore, based on the two-layer model, the mass and momentum balance equations of the mobile diffusion layer were added and a three-layer cuttings transport model was established. Kamp and Rivero [15] improved the solid phase concentration in the diffusion layer of the three-layer model by introducing the turbulence pulsation term, taking into account the promotion of turbulence in cuttings transportation. Cheng and Wang [16] divided the extended-reach well into three parts, where the horizontal and extended-reach well sections conformed to the characteristics of the three-layer model, and the moderately inclined well sections conformed to the characteristics of the two-layer model. The cuttings transportation in the vertical or near-vertical well sections only included the suspended layer. In addition to analyzing the control volume of the three-layer model, Ozbayoglu et al. [17] and Zhang et al. [18] also considered the influence of the force on the particles at the surface of the cuttings bed, making the three-layer model more complete. Wang et al. [19] solved the three-layer model using the SETS method and found that the cutting bed height and the drilling fluid displacement mutually affect each other, and that the well depth and displacement have a positive correlation with the cutting bed height. When the well depth reaches approximately 4000 m, the cutting bed height is no longer affected by the displacement and drill pipe rotation.
Regarding the incipient velocity model, when the velocity of the drilling fluid exceeds the critical velocity, the cuttings bed becomes unstable and continues to move forward. Compared to the two-layer and three-layer models, the incipient velocity model yields lower calculation values, providing a more pronounced advantage for parameter optimization. To ensure safe and efficient drilling, it is crucial to determine the incipient velocity accurately and reasonably. Clark and Bickham [20] conducted a force analysis on cutting particles at the bed surface, establishing the minimum transportation velocity (MTV) model. They identified three possible particle transportation states, including sliding, rolling, and lifting. Building on the MTV model, Ramadan et al. [21] considered the drag reduction effect of shear diluent liquid and introduced a downward plastic force, perpendicular to the flow direction of the force on the particles on the bed surface. The model results revealed that the rheological characteristics of the liquid phase also hindered particle start-up. Duan [11] incorporated the van der Waals force and used the near-bed velocity profile to establish an incipient velocity model for particles on the bed surface, accounting for factors such as particle size, rheology, and the inclination angle. Compared to other models, the incipient velocity model is easier to understand, providing a more direct guide for on-site engineers to optimize the discharge process. In recent years, in the field of fluid mechanics, machine learning methods have been used to process and analyze large amounts of experimental or simulated data to extract patterns or regularities related to the behavior of solid-liquid two-phase flows [22,23,24]. This combination can help researchers better understand and predict the characteristics of solid-liquid two-phase flows. In particular, GPU (Graphics Processing Unit) acceleration of machine learning has become a common method, which can significantly improve the speed of the training and inference processes. [25,26,27]
Currently, most incipient velocity models calculate the drag force based on the drag coefficient. However, there are numerous drag coefficient models available, which makes it challenging for field engineers to determine which model to use. This study analyzed the force applied to the particles on the surface of the cutting bed and the mechanism of the boundary layer flow to derive an analytical formula for the surface shear force of the drilling fluid’s flow around the particles on the cutting bed. Thus, this study established a critical velocity model that considered the effects of the flow around the particles. The effects of the well inclination angle, the cuttings size, and the different cuttings densities on the critical velocity of particle motion were analyzed experimentally, using a large-scale flow loop. The findings in this study can offer guidance in setting and predicting pump discharge, which, ultimately, could improve the performance of hole cleaning.

2. Experimental Study

2.1. Experimental Setup

The experiment was based on the fluid–solid coupling experimental equipment for drilling fluid and drill cuttings. This equipment is crucial for developing a better understanding of the behavior of drilling fluid in a wellbore, which is essential for achieving efficient and effective drilling operations. The experimental setup consisted of a transparent acrylic pipe, which simulated the annulus of a horizontal well section, with a total length of 5.5 m. The inner wellbore diameter of the acrylic pipe was 120 mm, which is a standard size for many oil and gas wells. The outer diameter of the drill pipe was 73 mm, which, again, is a common size for the drilling equipment used in the industry (Figure 1).

2.2. Test Materials

This study conducted experiments on the incipient velocity of drill cuttings using tap water and a poly aluminum chloride (PAC) solution. To prepare the PAC solution, this study mixed PAC powder with water and allowed the mixture to hydrate for a sufficient period. Prior to and after the experiments, the rheological properties of the test fluids were measured using an Anton Paar MCR92 advanced rheometer. The rheology of the PAC solution is best described by a power-law fluid model (k = 0.28 Pa.sn, n = 0.65). The rheology of the tap water is best described by a Newtonian fluid model (μ = 0.0013 Pa.s). The solid phase of experiment consisted of glass and ceramic with different densities and particle sizes.

2.3. Experimental Setup

To accurately model incipient velocity, a stable cuttings bed with particles at a specific degree of subsidence and incipient motion is required. For deviated or horizontal well sections, the cuttings bed must be prepared in advance. Therefore, the incipient velocity measurement experiment was conducted as follows:
Step 1. Fluid and Particle Preparation: in this study, the shale shaker and recovery device were cleaned, then a specific concentration of experimental drilling fluid mother liquor was prepared before the experiment and it was stored it in the liquid storage tank for stirring and dilution. This study also added a sufficient quantity of experimental particles to the sand injection tank.
Step 2. Cuttings Bed Deposition: in this study, the sand-adding tank was opened to allow the experimental particles to enter the annulus using higher displacement. This created enough particles to form a stable cuttings bed. This study then reduced the displacement, kept the cuttings bed stationary and distributed the particles uniformly on the existing cuttings bed.
Step 3. Cuttings Bed Height Calculation: in this study, the displacement was reduced which made the particles static after sinking. Once a stable cuttings bed was formed, this study used protruding particles on the bed surface as the observation objects. Next, we took fixed-point photos and measured the bed height. Since the outer wall of the annulus is circular, the arc length was measured to obtain the corresponding cuttings bed height.
Step 4. Incipient Velocity Observation: in this study, we increased the flow rate of the drilling fluid by 0.004 m/s, until the observed particles on the bed surface started rolling and moving forward. To study particle motion on a bed surface, this study recorded the discharge required to initiate movement, reduced the discharge until the particles were stationary, and repeated the process.
During the incipient motion of a particle on the cuttings bed, a high-speed camera (Revealer 2F04C) was used to capture the incipient motion process at a frame rate of 100 frames/sec, as shown in Figure 2. The camera tracks particles using Lagrangian frames of reference. The observed particles on the bed surface first slowly rolled forward, then rolled slightly backwards, and, finally, continued to roll forward under the continuous drag of drilling fluid. These observations are depicted in Figure 2.

3. Results and Analysis

3.1. Effect of Different Inclination Angles

The influence of the inclination angle on incipient velocity is illustrated in Figure 3. It is apparent from the figure that the incipient velocity was the highest in the wellbore sections with an inclination angle of 60–80°, making the incipient motion of drill cuttings particularly arduous.
The analysis of the forces acting on the drill cuttings indicated that only the net weight force in the direction of the particle lifting was affected by the inclination angle, whereas all the other forces were independent of the wellbore inclination. Furthermore, the net weight force in the direction of the particle lifting increased with an increase in the wellbore inclination angle. In the horizontal wellbore sections, the drill cuttings were the least likely to be lifted and were primarily transported by a rolling mechanism. In the wellbore sections with a small inclination angle, the resistance force in the direction of the particle lifting was minimal, making the drill cuttings more prone to being lifted into the mainstream of the drilling fluid. In fact, in near-vertical wellbore sections, the gravitational force acting on the drill cuttings in the direction of settling towards the wellbore wall was entirely nullified, thereby fully entraining rock particles into the main flow.

3.2. Effect of Different Cuttings Sizes

The variations in the incipient velocity under different cutting size conditions are presented in Figure 4. The incipient velocity increased with the increase in the cutting sizes under two different drilling fluid conditions. The velocity required for the drill cuttings particles to move in the PAC experimental fluid was lower than that in the tap water experimental fluid. The specific pattern of the cuttings movement tended to remain constant within a smaller range of particle sizes, e.g., within the range of 1.0–2.0, and increased rapidly when the particle size exceeded a certain value.
This phenomenon can be explained by analyzing the changes in the inertial forces concerning particle size. On the surface of the cuttings bed, the drag force of the drilling fluid caused the particles to start moving. When the particle size was small, the inertial force was greater than the shear stress applied to the particles and the movement of the particles on the surface of the cuttings bed tended to be stable. However, when the particle size exceeded a certain threshold, the inertial force experienced by the moving particles was smaller than the drag force of the drilling fluid, which made the effect of the particle size on the movement velocity more pronounced.

3.3. Effect of Different Cuttings Densities

The influence of the cuttings density on the incipient velocity is illustrated in Figure 5. When the particles on the surface of the cuttings bed were mobilized, the required flow rate of the drilling fluid increased with the density of the particles. The velocity required for the mobilization of the cuttings particles in the PAC experimental fluid was lower than that in the tap water experimental fluid. Increasing the density of the cuttings had a negative impact on their mobilization, leading to an increase in the height of the cuttings bed and a decrease in the effectiveness of the cleaning in the annular space.
The incipient velocity of drill cuttings increased with the increase in their density. The reason for this was that, as the density of the cuttings increased, the mass of the individual drill cuttings also increased, even though the volume remained the same. As a result, the energy required for the drill cuttings to start moving from a static state also increased, along with the incipient velocity. Due to the increase in the cuttings density, the decrease in the particle buoyancy, and the increase in the submerged weight, the average flow rate for the drilling fluid required for incipient velocity increased.

4. Model Development and Analysis

In this study, the following assumptions were made before establishing the incipient velocity model: (1) the drilling fluid is an incompressible medium; (2) spheres replace the actual irregular cutting particles with the same diameter; (3) the average bed surface replaces the actual cutting bed surface; and (4) the annulus is concentric, and the drill pipe does not rotate. Based on these assumptions, the incipient motion of the cutting particles on the bed surface was analyzed and the model was derived.

4.1. Shear Force of Flow around Cutting Particles on Cuttings Bed Surface

Due to the wellbore and the drill pipe forming a double-wall annular space, there was a velocity gradient distribution present in both the drill pipe and the wellbore. The effect of this velocity gradient distribution on particles can be simplified as the effect of different velocity layers on the particle.
The main force that moves particles on the cuttings bed surface is the drag force of the drilling fluid’s flow along the flow direction. Since the annular space formed by the wellbore and the drill pipe had two walls, the annular drilling fluid had dual velocity gradient distribution characteristics. As shown in Figure 6, there was a maximum velocity on the annular surface, with a certain radius of the annulus.
Assuming that the radius, Rc, corresponds to the maximum velocity of the annular surface, it can be calculated using the following expression:
R c = ( R p + R w ( R p R w ) 7 / 20 ) / ( 1 + ( R p R w ) 7 / 20 )
where Rp is the drill pipe radius and Rw is the wellbore radius.
The exposure level, ζ, was introduced to define the degree to which a particle was exposed to the flow:
ζ = ζ r s
where ζ’ is the distance from the lowest point of the particles to the average bed surface; rs is the average radius of the drill cuttings.
As shown in Figure 7, the calculation point was selected along the circumference as the origin to establish a Cartesian coordinate system. The potential flow velocity distribution around the circumference of the volume element was as follows:
U ( α ) = 2 u ( r 1 ) sin α
where u(r1) is the corresponding mainstream velocity of each layer’s volume element; r1 is the distance from the particle center to the central axis.
The maximum flow velocity, umax, in the annulus can be expressed as follows:
u max = u ¯ + R p R c u 1 ( r ) d r + R c R w u 2 ( r ) d r
where umax is the maximum flow velocity in the annulus; u ¯ is the average flow velocity above the bed surface, which can be calculated as the ratio of the volume flow rate of the flow area above the bed surface to the flow-section area; u1(r) is the velocity distribution from the drill pipe to the maximum velocity annular surface; and u2(r) is the velocity distribution from the maximum velocity annular surface to the wellbore.
The expression for the velocity distribution in the annulus can be summarized as follows:
u i ( r 1 ) = u max + u * i 2 κ ( R c r 1 Δ R ) + u * i 2 κ ( r s Δ R ) cos β
where κ is the Karman constant, with κ = 0.41 when the velocity is distributed between Rp and Rc (i.e., i = 1, ΔR = RcRp), and when the velocity distribution is between Rc and Rw (i.e., i = 2, ΔR = RwRc); u* is the friction velocity, u * = τ / ρ ; u*1 is the friction velocity on the outer wall of the drill pipe; and u*2 is the friction velocity on the inner wall of the wellbore.
The surface shear force can be expressed as follows:
τ = 0.5 f f ρ u ¯ 2
where ff is the average friction coefficient of the liquid flow [28], which can be calculated as follows:
1 f f = 2 lg ( ε 3.72 D h + 2.51 Re f f )
where ε is the surface average roughness; Dh is the effective hydraulic diameter; and Re is the Reynolds number.
Let us suppose that the velocity distribution at the calculation point and the vertical height, y, from the particle surface meets u(y) = ay2 + by + c, and that the boundary conditions are met as follows:
y = 0 ,   u ( y ) | 0 = 0 y = δ ,   u ( y ) | δ = U ( α ) ,   u ( y ) y | δ = 0
According to the simplification of the integral equation, in this study we obtained the following:
u ( r 1 ) 5 ν r β ( 3 δ 2 cos α + 2 3 δ δ sin α ) = 1
where rβ is the radius of the corresponding area element of β; v is the fluid kinematic viscosity; δ is the boundary layer thickness for the flow around the particle; and δ´ represents the derivative of the boundary layer thickness, with respect to radian α.
The fluid kinematic viscosity was calculated according to the following formula:
ν = K ρ [ ( 2 n + 1 ) ( R w R p ) 18 n u ¯ ] 1 n
where n is the flow index and K is the consistency coefficient.
The boundary layer thickness, δ(α), can be calculated as follows:
δ ( α ) = b C 2 cos α H y p e r g e o m e t r i c 2 F 1 [ 1 2 , 7 2 ; 3 2 ; cos 2 α ] b sin 9 α
where H y p e r g e o m e t r i c 2 F 1 [ 1 2 , 7 2 ; 3 2 ; cos 2 α ] is a Gaussian hypergeometric function, which is expanded by series. When determining the constant C, it is assumed that the initial boundary layer thickness, δ0, starts to develop from the front stagnation point, α0, as follows:
δ 0 = 0.00655864 b sin 9 α 0
Based on this assumption, in this study we obtained the following:
δ = 0.859 ( 2 cos α 7 3 cos 3 α + 7 4 cos 5 α 5 8 cos 7 α ) b sin 9 α
In the above formula, b = 2 u ( r 1 ) 15 v r β .
The shear force acting on any point around the circumference of the element is as follows:
τ 0 = 2 ν ρ U ( α ) δ = ( 2 ν ρ 2 u ( r 1 ) sin α ) 2 b sin 9 α 0.859 ( 2 cos α 7 3 cos 3 α + 7 4 cos 5 α 5 8 cos 7 α )
The total shear force, Fτa, is obtained by integrating the shear stress along the undivided circumference and then along the axis, l.
F τ α = 16 ρ v r s 3 15 ( 0 π sin β u 3 2 ( r 1 ) 2 F ( α ) | 0 9 π 20 d β β 1 β 2 sin β u 3 2 ( r 1 ) F ( α ) | π 2 θ 9 π 20 d β )
The function, F(α), is obtained by the regression of the integral upper limit function and the integral value, as follows:
F ( a ) = 0.1359 a 3 + 1.0201 a 2 0.5179 a + 0.0426
where θ is the included angle between the line from the particle exposure point to the particle center and the vertical direction.
β1 and β2 meet the following requirements:
β 1 = π r arccos r s ζ r s β 2 = π γ + arccos r s ζ r s
where γ is the position angle of the particles on the bed surface, as shown in Figure 2.
The particles on the bed surface are affected by both the annulus velocity gradient and also by the internal shear force of the viscous sublayer, which is expressed as Fτb. As shown in Figure 8, the volume differential element, C, was integrated from the outer boundary of the viscous sublayer to the inner boundary, and the boundary layer thickness can be expressed as follows:
δ p = 5 u * n 2 n ( K ρ ) 1 n
τ b = 0.5 f b ρ u ¯ 2
u * = τ b / ρ
where fb is the average friction coefficient of the bed surface; τb is the average shear stress of bed surface; and δp is the bed boundary layer thickness.
At the same time:
{ β out = arccos ( 1 ζ δ p / r s ) β in = arccos ( 1 ζ )
where βout is the outer boundary angle of the viscous sublayer integral, which refers to the angle from the top to the outer boundary of the viscous sublayer; and βin is the integral internal boundary angle of the viscous sublayer, which is the angle from the top of the particle to the contact point with the average bed surface.
The velocity distribution inside the viscous sublayer is expressed as follows:
u in = u * 2 v r s ( ( 1 ζ ) + cos β )
The integral expression can be calculated as follows:
F τ b = 8 ρ 2 v r s 3 15 β out β in ( u i 3 sin β ) 0.5 F ( a ) | 0 9 π 20 d β
The total surface shear force along the flow direction, Fτ, was obtained from the superposition of the shear force of the viscous sublayer’s internal bed boundary layer and the external annulus’s surface shear force, as follows:
F τ = F τ a + F τ b

4.2. Lifting Force of Bed Particles

The particles are also subject to a lifting force when the drilling fluid flows through the surface of the cuttings bed. The method for calculating the velocity, ud, varies with the deposition degree of the cuttings bed.
When the deposition of the drill cuttings was dense and the flow was laminar, the linear Darcy’s law was used for the calculation, as follows.
u d = Δ p K p μ L
where Kp is the permeability coefficient of the cutting bed; Δp/L is the pressure drop gradient in the L measuring section; and μ is the dynamic viscosity of the drilling fluid.
When the drill cuttings deposition is non-dense, non-Darcy’s law applies, as follows:
Δ p L = μ K p u d + β D ρ u d 2
where βD is the non-Darcy factor, which can be calculated as follows:
β D = 1.59 × 10 5 q 5.5 K p
where q is the average porosity of the cuttings bed and Kp is the porous media permeability, which can be expressed as follows:
K p = c 0 q 3 ( 1 q ) 2 M s 2
where C0 is the coefficient (taking C0 = 0.2); Ms is the specific surface area of the particles (for spherical particles, Ms = 3/rs).
According to Bernoulli’s principle, the lifting force and the average flow velocity on the bed surface, u ¯ , as well as the flow velocity inside the bed, ud, are related as follows:
F L = 0.5 ρ A s ( u ¯ 2 u d 2 )
where Δps is the difference between the upper and lower pressures of the particles; and As is the projected area of the particles that are perpendicular to the liquid flow direction.

4.3. Interaction Force between Bed Particles

There was also an adhesive force between particles, which is expressed as follows:
F p = 2 ( 5 3 ) m ξ s j m r s
where ξ s = 4.744 × 10 6   kg / m , m = 10.
With respect to the ratio, j, of the dry volume weight and the volume weight, its expression is as follows:
j = { 0.52 ( r s r s + 2 δ s ) 3 r s 0.5   mm 0.70 0.18 exp ( 0.095 r s r 0 r 0 ) r s > 0.5   mm
where δs is the thickness of the film water on the particles’ surface, δs = 0.4 μm; and r0 is the reference particles’ radius, r0 = 0.5 mm.

4.4. Incipient Velocity Calculation Process

First, this study determined the surface friction, Fτ, between the particles and the bed surface; the lifting force, FL; and the interaction force, Fp, between the particles. Then, this study determined the incipient velocity required for the particles to start rolling, according to the condition that the total positive torque was greater than or equal to the total negative torque.
As shown in Figure 9, for highly deviated and horizontal well sections, the conditions for the particle rolling to start are as follows:
( F τ a + F τ b ) L 1 + ( F L F P ) L 2 F g L 3 > 0
L 1 = r s ( 1 ζ )
L 2 = r s sin [ arccos ( 1 ζ ) ]
L 3 = r s cos [ φ arccos ( 1 ζ ) ]
where φ is the well deviation angle; Fg is the floating weight of the particles in the drilling fluid, as follows:
F g = 4 3 π r s 3 ( ρ s ρ ) g
where ρs is the average density of the drill cuttings.
The calculation steps for determining the incipient velocity are depicted in Figure 10. The surface friction was computed based on the corresponding integral value. Next, the calculations of the lifting force, inter-particle interaction force, and other relevant factors were conducted. The iterative calculation of the average velocity continued until the critical rolling motion conditions were met, after which the velocity was deemed the incipient velocity, uc.

4.5. Measured and Predicted Initial Velocity

In this paper, in addition to predicting the critical velocity of the particles on the cuttings bed surface through theoretical derivation, it was also necessary to compare and analyze the incipient velocity predicted by the model with the experimental measurement data. The comparison between the model’s prediction results and the experimental measurement results under three operating conditions is shown in Figure 11. As shown in Figure 11, the maximum relative deviation between the simulation and experimental results did not exceed 20%. This proves that the model established was relatively accurate and will be able to further simulate the incipient velocity under variable operating conditions. In addition, the relative error with the small velocity was relatively large, while the error after that was small. This was because the fluid field in the annulus underwent significant changes during the initial motion with small-circulated velocity. This error will decrease when the circulated velocity stabilizes.

5. Conclusions

In this paper, the effects of the well’s inclination angle, the cuttings size, and the different cuttings densities on the critical velocity of particle motion were analyzed experimentally, using a large-scale flow loop. Next, based on a stress analysis of the particles on the cuttings bed surface and of the boundary layer flow around the particles, an analytical formula for the surface shear force of the drilling fluid on the particles was derived and a critical velocity model was established. A verification was carried out and combined with the experimental results, as follows:
(1)
The surface shear force can be calculated analytically from the angle of flow around the surface of the particles on the bed. The drag force can be approximately replaced by the surface shear force. The influence of the differential pressure resistance should be taken into account when calculating the drag force.
(2)
By using the critical velocity, the critical displacement of the effective drilling fluid needed to maintain the erosion of the cuttings bed can be determined. The sinking of the cuttings bed has a greater impact on small particle sizes, while the effect of the pressure difference resistance on them is smaller.
(3)
When the drilling fluid flows over the surface of the particles, the surface shear force plays a major driving role. If the particle size increases, the negative pressure formed by the vortex behind the particles also increases, which cannot be ignored.

Author Contributions

Methodology, K.S.; software, G.J. and Y.W.; validation, Z.X.; writing—original draft, K.S., Z.X. and G.J.; writing—review and editing, G.J. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Marine Economy Development Foundation of Guangdong Province (GDNRC [2022]44) and the National Natural Science Foundation of China (Major Program, 51991361).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yan, T.; Qu, J.Y.; Sun, X.F.; Li, W.; Chen, Y.; Hu, Q.B. A novel predictive model of drag coefficient and settling velocity of drill cuttings in non-Newtonian drilling fluids. Pet. Sci. 2021, 18, 1729–1738. [Google Scholar] [CrossRef]
  2. Heydari, O.; Sahraei, E.; Skalle, P. Investigating the impact of drillpipe’s rotation and eccentricity on cuttings transport phenomenon in various horizontal annuluses using computational fluid dynamics (CFD). J. Pet. Sci. Eng. 2017, 156, 801–813. [Google Scholar] [CrossRef]
  3. Hu, Q.; Li, Y.; Sun, X.; Chen, M.; Bu, Q.; Gong, B. Integrating test device and method for creep failure and ultrasonic response of methane hydrate-bearing sediments. Rev. Sci. Instrum. 2023, 94, 025105. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, F.; Miska, S.; Yu, M.; Ozbayoglu, E.; Takach, N.; Osgouei, R.E. Is well clean enough? A fast approach to estimate hole cleaning for directional drilling. In Proceedings of the Society of Petroleum Engineers—Coiled Tubing and Well Intervention Conference and Exhibition 2015, The Woodlands, TX, USA, 24–25 March 2015; pp. 603–618. [Google Scholar]
  5. Mahmoud, H.; Hamza, A.; Nasser, M.S.; Hussein, I.A.; Ahmed, R.; Karami, H. Hole cleaning and drilling fluid sweeps in horizontal and deviated wells: Comprehensive review. J. Pet. Sci. Eng. 2020, 186, 106748. [Google Scholar] [CrossRef]
  6. Alawami, M.; Bassam, M.; Gharbi, S.; Al Rubaii, M. A real-time indicator for the evaluation of hole cleaning efficiency. In Proceedings of the Society of Petroleum Engineers—SPE/IATMI Asia Pacific Oil and Gas Conference and Exhibition 2019, APOG 2019, Bali, Indonesia, 29–31 October 2019; pp. 29–31. [Google Scholar]
  7. Naganawa, S.; Kudo, H.; Matsubuchi, H. Simulation study on influences of wellbore tortuosity on hole cleaning in extended-reach drilling. In Proceedings of the Society of Petroleum Engineers—Abu Dhabi International Petroleum Exhibition and Conference 2016, Denver, CO, USA, 5–6 May 2016. [Google Scholar]
  8. Naganawa, S.; Sato, R.; Ishikawa, M. Cuttings-Transport Simulation Combined With Large-Scale-Flow-Loop Experimental Results and Logging-While-Drilling Data for Hole-Cleaning Evaluation in Directional Drilling. SPE Drill. Complet. 2017, 32, 194–207. [Google Scholar] [CrossRef]
  9. Naganawa, S.; Nomura, T. Simulating transient behavior of cuttings transport over whole trajectory of extended-reach well. In Proceedings of the IADC/SPE Asia Pacific Drilling Technology Conference and Exhibition, Bangkok, Thailand, 13–15 November 2006; pp. 480–488. [Google Scholar] [CrossRef]
  10. Yan, X.; Bu, Y.; Cheng, R. Research on numerical simulation of transient cuttings transport with foam fluid in horizontal section. Drill. Fluid Complet. Fluid 2011, 28, 24–28. [Google Scholar] [CrossRef]
  11. Duan, M.; Miska, S.Z.; Yu, M.; Takach, N.E.; Zettner, R.M.A.C.M. Critical Conditions for Effective Sand-Sized Solids Transport in Horizontal and High-Angle Wells. SPE Drill. Complet. 2009, 24, 10. [Google Scholar] [CrossRef]
  12. Huo, H.J.; Wang, R.H.; Ni, H.J.; Li, Y.X.; Tan, C.Y.; Xue, S. Study of critical annulus up-returning velocity of cuttings carried by supercritical CO2 in deviated well. J. CO2 Util. 2017, 20, 105–112. [Google Scholar] [CrossRef]
  13. Qu, J.; Sun, X.; Yan, T.; Chen, Y.; Hu, Q. A novel prediction model for cuttings bed height in horizontal and deviated wells. J. Nat. Gas Sci. Eng. 2021, 95, 104192. [Google Scholar] [CrossRef]
  14. Nguyen, D.; Rahman, S.S. Three-layer hydraulic program for effective cuttings transport and hole cleaning in highly deviated and horizontal wells. In Proceedings of the SPE/IADC Asia Pacific Drilling Technology, Kuala Lumpur, Malaysia, 9–11 September 1996; pp. 163–173. [Google Scholar] [CrossRef]
  15. Kamp, A.M.; Rivero, M. Layer modeling for cuttings transport in highly inclined wellbores. In Proceedings of the Latin American and Caribbean Petroleum Engineering Conference, Caracas, Venezuela, 21–23 April 1999. [Google Scholar] [CrossRef]
  16. Cheng, R.; Wang, R. A Three-Segment Hydraulic Model for Annular Cuttings Transport with Foam in Horizontal Drilling. J. Hydrodyn. 2008, 20, 67–73. [Google Scholar] [CrossRef]
  17. Ozbayoglu, M.E.; Miska, S.Z.; Reed, T.; Takach, N. Using foam in horizontal well drilling: A cuttings transport modeling approach. J. Pet. Sci. Eng. 2005, 46, 267–282. [Google Scholar] [CrossRef]
  18. Zhang, F.; Miska, S.; Yu, M.; Ozbayoglu, E.M.; Takach, N. Pressure Profile in Annulus: Solids Play a Significant Role. J. Energy Resour. Technol. Trans. ASME 2015, 137, 064502. [Google Scholar] [CrossRef]
  19. Wang, Z.; Zhai, Y.; Hao, X.; Guo, X.; Sun, L. Numerical Simulation on Three Layer Dynamic Cutting Transport Model and Its Application on Extended Well Drilling. In Proceedings of the IADC/SPE Asia Pacific Drilling Technology Conference and Exhibition, Ho Chi Minh City, Vietnam, 1–3 November 2010; pp. 1–10. [Google Scholar] [CrossRef]
  20. Clark, R.K.; Bickham, K.L. Mechanistic model for cuttings transport. In Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, LA, USA, 25–28 September 1994; pp. 139–153. [Google Scholar] [CrossRef]
  21. Ramadan, A.; Saasen, A.; Skalle, P. Application of the minimum transport velocity model for drag-reducing polymers. J. Pet. Sci. Eng. 2004, 44, 303–316. [Google Scholar] [CrossRef]
  22. Mehmood, Y.; Shafqat, R.; Sarris, I.E.; Bilal, M.; Sajid, T.; Akhtar, T. Numerical Investigation of MWCNT and SWCNT Fluid Flow along with the Activation Energy Effects over Quartic Auto Catalytic Endothermic and Exothermic Chemical Reactions. Mathematics 2022, 10, 4636. [Google Scholar] [CrossRef]
  23. Barkalov, K.; Lebedev, I.; Usova, M.; Romanova, D.; Ryazanov, D.; Strijhak, S. Optimization of Turbulence Model Parameters Using the Global Search Method Combined with Machine Learning. Mathematics 2022, 10, 2708. [Google Scholar] [CrossRef]
  24. Li, S.; Chen, Z.; Li, W.; Yan, T.; Bi, F.; Tong, Y. An FE Simulation of the Fracture Characteristics of Blunt Rock Indenter Under Static and Harmonic Dynamic Loadings Using Cohesive Elements. Rock Mech. Rock Eng. 2023, 56, 2935–2947. [Google Scholar] [CrossRef]
  25. Zhao, M.; Jha, A.; Liu, Q.; Millis, B.A.; Mahadevan-Jansen, A.; Lu, L.; Landman, B.A.; Tyska, M.J.; Huo, Y. Faster Mean-shift: GPU-accelerated clustering for cosine embedding-based cell segmentation and tracking. Med. Image Anal. 2021, 71, 102048. [Google Scholar] [CrossRef]
  26. Li, Y.; Tang, B.; Geng, B.; Jiao, S. Fractional Order Fuzzy Dispersion Entropy and Its Application in Bearing Fault Diagnosis. Fractal Fract. 2022, 6, 544. [Google Scholar] [CrossRef]
  27. Zhao, M.; Liu, Q.; Jha, A.; Deng, R.; Yao, T.; Mahadevan-Jansen, A.; Tyska, M.J.; Millis, B.A.; Huo, Y. VoxelEmbed: 3D Instance Segmentation and Tracking with Voxel Embedding Based Deep Learning; Springer International Publishing: Berlin/Heidelberg, Germany, 2021; Volume 12966, ISBN 9783030875886. [Google Scholar] [CrossRef]
  28. Dodge, D.W.; Metzner, A.B. Turbulent flow of non-newtonian systems. AIChE J. 1959, 5, 189–204. [Google Scholar] [CrossRef]
Figure 1. Photograph of the critical velocity of particle motion experiment.
Figure 1. Photograph of the critical velocity of particle motion experiment.
Mathematics 11 02689 g001
Figure 2. Schematic diagram of incipient motion of bed particles. (Red circle indicates observed particles).
Figure 2. Schematic diagram of incipient motion of bed particles. (Red circle indicates observed particles).
Mathematics 11 02689 g002
Figure 3. Influence of inclination angles on the incipient velocity of drill cuttings.
Figure 3. Influence of inclination angles on the incipient velocity of drill cuttings.
Mathematics 11 02689 g003
Figure 4. Influence of cuttings sizes on the incipient velocity of drill cuttings.
Figure 4. Influence of cuttings sizes on the incipient velocity of drill cuttings.
Mathematics 11 02689 g004
Figure 5. Influence of cuttings densities on the incipient velocity of drill cuttings.
Figure 5. Influence of cuttings densities on the incipient velocity of drill cuttings.
Mathematics 11 02689 g005
Figure 6. Geometric relationship of the control volume of drill cuttings bed surface in the annulus.
Figure 6. Geometric relationship of the control volume of drill cuttings bed surface in the annulus.
Mathematics 11 02689 g006
Figure 7. Cuttings particles incipient motion at different positions on the cuttings bed surface.
Figure 7. Cuttings particles incipient motion at different positions on the cuttings bed surface.
Mathematics 11 02689 g007
Figure 8. Geometric relationship for boundary layer of flow around the particles on the bed surface.
Figure 8. Geometric relationship for boundary layer of flow around the particles on the bed surface.
Mathematics 11 02689 g008
Figure 9. Schematic diagram of bed particles in the incipient state.
Figure 9. Schematic diagram of bed particles in the incipient state.
Mathematics 11 02689 g009
Figure 10. Calculation flow chart for incipient velocity of bed particles.
Figure 10. Calculation flow chart for incipient velocity of bed particles.
Mathematics 11 02689 g010
Figure 11. Calculation flow chart of incipient velocity for bed particles.
Figure 11. Calculation flow chart of incipient velocity for bed particles.
Mathematics 11 02689 g011
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

Shen, K.; Xu, Z.; Ju, G.; Wang, Y.; Li, Z. Estimating the Critical Velocity of the Incipient Motion of Particles on the Cuttings Bed Surface: An Experimental and Theoretical Analysis. Mathematics 2023, 11, 2689. https://doi.org/10.3390/math11122689

AMA Style

Shen K, Xu Z, Ju G, Wang Y, Li Z. Estimating the Critical Velocity of the Incipient Motion of Particles on the Cuttings Bed Surface: An Experimental and Theoretical Analysis. Mathematics. 2023; 11(12):2689. https://doi.org/10.3390/math11122689

Chicago/Turabian Style

Shen, Kaixiang, Zhenqiang Xu, Guoshuai Ju, Yingsheng Wang, and Zijian Li. 2023. "Estimating the Critical Velocity of the Incipient Motion of Particles on the Cuttings Bed Surface: An Experimental and Theoretical Analysis" Mathematics 11, no. 12: 2689. https://doi.org/10.3390/math11122689

APA Style

Shen, K., Xu, Z., Ju, G., Wang, Y., & Li, Z. (2023). Estimating the Critical Velocity of the Incipient Motion of Particles on the Cuttings Bed Surface: An Experimental and Theoretical Analysis. Mathematics, 11(12), 2689. https://doi.org/10.3390/math11122689

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