2018 Movehedi Et Al. - Modeling and Simulation of Barite Deposition in An Annulus Space of A Well
2018 Movehedi Et Al. - Modeling and Simulation of Barite Deposition in An Annulus Space of A Well
A R T I C L E I N F O A B S T R A C T
Keywords: In drilling industry, barite particles settling and barite sag as a major problem can potentially impose significant
Barite settling operational issues. Static conditions, in which well undergoes an extended shut-in period, could occur during
Eulerian-eulerian model different drilling and completion operations such fishing operation, tripping, and logging. Despite its importance,
Single particle model such phenomenon is not well understood yet. To avoid issues related to barite settlement and barite sag, a good
CFD simulation understanding of the impact of different drilling parameters on barite settlement and sag phenomenon is required.
Recently, the mathematical formulation and modeling of settlement and sag processes have gained more atten-
tion. In order to better design the drilling fluids and to optimize the operational conditions such as flow rate and
drill pipe rotation speed, as well as drilling well trajectory design to avoid barite settlement, the mathematical
formulation and modeling of such physical phenomena is required. This study is aimed at mathematically
describing and analyzing the settling of barite particles in annular concentric space without axial flow by the
means of computational fluid dynamic (CFD) simulation. The modeling approach consists of two main steps: a)
simple vertical annular domain and b) complex three dimensional barite settlement process in a deviated well
with rotational drill pipe. The current work is based on Eulerian two phase model. The settling of the single
particle and the hindering effect of other particles on settling velocity have been investigated. Based on single
particle model, the rate of barite settlement and fluid pressure reduction have been modeled. However, the single
particle model cannot describe the particle to particle interactions and their impacts on settlement such as particle
shearing viscosity, kinetic particle viscosity, and collisional effects. In this study, the accurate and comprehensive
modeling of barite settlement has been done by applying the Eulerian-Eulerian approach. The models have been
validated against available published experimental data in the literature. In addition, in order to further validate
the model presented here, drill pipe rotation, fluid rheology, and inclination angle impacts have been experi-
mental investigation and compared with results of the proposed model. The results indicate that the model based
on using Eulerian-Eulerian approach accurately predicts the experimental data results. Afterward, the sensitivity
analyses were carried out to show the impact of inclination angle, inner pipe rotation, fluid rheology, and barite
content on the rate of settlement. This study showed that drag forces induced resuspension plays a vital role in
barite settlement inside a deviated well. The results reveal that the settling rate is decreased with increasing the
inner pipe rotation and the vortex effect created by inner pipe rotational, play a major role in the reduction of
barite settlement. The results also show that the maximum settlement rate accrues at the deviation angles between
of 45–60 .
1. Introduction such as barite and Hematite play a significant role providing the required
hydrostatic pressure. Barite settling can be defined as a phenomenon in
In petroleum industry, drilling fluids are used to not only remove which weighing materials settle due to the impact of gravitational forces
drilling cuttings but also to lubricate and cool the drilling bit and drill applied on barite particles. Well trajectories as well as properties of mud
string. The other functionality of drilling fluids is to control the formation and operational conditions are among factors that influence the degree of
pressure by providing enough hydrostatics pressure. Weighting agents barite settling occurrence. For instant, in directional wells, addition of
* Corresponding author.
E-mail addresses: [email protected], [email protected] (S. Shad).
https://doi.org/10.1016/j.petrol.2017.12.014
Received 19 March 2017; Received in revised form 24 November 2017; Accepted 5 December 2017
Available online 7 December 2017
0920-4105/© 2017 Elsevier B.V. All rights reserved.
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
drag reducing material to the oil based mud could potentially lead to involved in this process.
significant barite settlement. Operational issues and challenges such as Drilling management and optimization requires an accurate modeling
tripping, logging, waiting on materials, reaming and fishing increases the and simulation of settlement process. Despite its importance, modeling
chance and magnitude barite settlement inside a well. This is due to the and CFD based simulation of barite settlement has not received much
fact that under such conditions the drilling process will stop and well will attention. Kaushal et al. (2012) studied and modeled the flow of
be shut down for days or weeks allowing barite particles to settle. Due to mono-dispersed fine particles slurry in a pipe at high concentration.
the settlement of barite particles, well will experience a reduction of Recently, Hashemian et al. (2014) conducted a simulation study on barite
density with depth. The density of the fluids is directly related to barite sag phenomenon. However, the impact of vortex created by rotational
concentration. Non-uniform density and barite concentration could lead shaft and spiral boundary layer on particle velocity and settlement has
to the hydrostatic pressures smaller than that required at the zone of not been fully analyzed.
interest. The higher reservoir pressure could then lead to well blow out Lately, the application of CFD methods in two phase flow and
and loss of well control. The settlement of barite can lead to creation of determination of velocity distribution in annular space during drilling
different barite concentration profile. Fig. 1 shows the different segments has gained more attention. Hajidavalloo et al. (2013) have conducted
of separated fluids at the heel of a well including slump, sag, and sus- simulation study on gas solid transport in annular space using CFD
pension and clarified zones (Omland, 2009). Despite the differences methods that are applicable in air drilling process. Their simulation study
between barite sag and settlement, the barite sag terminology has often is based on a Euler-Euler two phase model. They studied the drill pipe
been used to explain settlement process. rotation and well direction effects on volume fraction of particle in the
Despite its importance, barite settlement is not well understood yet. gas stream. Mme and Skalle (2012) conducted a simulation study on
Most of the researchers have focused on experimental analysis of such cutting transport phenomenon in concentric annular space using CFD
process and the modeling and simulation of this phenomenon has method. This study is focused on cuttings which are much larger than
received rather less attention. In order to achieve a good understanding barite particles. The impacts of particle size distribution, different
of barite settling process, it is vital to study the dynamics of fluid flow in annular velocity and angle of well on transport of drilling cuttings have
pipes under stationary and dynamic conditions. been investigated in their work. Akhshik et al. (2016) have performed
Hanson et al. (1990) and Saasen et al. (1995) are among pioneers who simulation on hole cleaning process in a deviated well using CFD-DEM
experimentally evaluated the barite sag or settlement in inclined tubes. simulation. They applied a two-way coupled CFD–DEM model to simu-
Their study suggested that barite sag is similar to that introduced by late cutting transport including fluidparticle interactions. The parti-
Boycott (1920) for blood corpuscles settlement. Later, in 1998, Bern et al. cleparticle, particlewall, and particledrill pipe interactions were
(1998) introduced an experimental model to predict barite settlement. modeled using the HertzMindlin model. The effects of particle shape on
Even though that their model predicts the barite settlement for stationary the deposition phenomenon also have been investigated at different inlet
pipe system, it fails to provide an accurate estimate of such process for velocities and drill pipe inclination angles. Other studies in the field of
conditions in which the pipe is rotating. cutting transport including experimental-modeling and simulation has
Skalle et al. (1999) experiments have shown that barite sag occurs been well documented (Dritselis, 2017; Elgaddafi et al., 2012; Kamyab
under both static and dynamics conditions. Skalle et al. (1999) also and Rasouli, 2016; Mishra and Ein-Mozaffari, 2017; Mondal et al., 2016;
showed that settling process will increase during laminar flow condition. Movahedi et al., 2017; Shahi and Kuru, 2016; Wang and Long, 2010).
Nguyen et al. (2009) conducted several experiments inside an annular Understanding the behavior of flow through the annulus gap in a
flow loop and showed that the particle settling can be mitigated at lower vertical or deviated position is important in modeling of fine size barite
rotation speeds. Nguyen et al. (2014) applied Taguchi and ANOVA particle settlement. In static conditions with no axial flow with an inner
methods to design the experiments and to study the effects of drilling shaft rotation, particles tend to follow a circular path due to the lift force
parameters on barite settlement under dynamic conditions. They showed that exerted on such particles. This force is resulted from rotation of the
that, in comparison with different drilling parameters, the fluid velocity inner shaft or drill pipe. The pipe rotation will change the particle settling
had the most dominant effect on preventing settlement. Parvizinia et al. regime such that a simple gravity dominated settlement will change into
(2011) conducted an experimental study on barite sag behaviors of a complex process. Barite particles are significantly smaller than drilling
oil-based fluids. The results indicated that fluids subjected to high shear cuttings and therefore vortices created by pipe rotation will have a
rates and elevated temperature experience significant barite sag. Amighi greater impact on them.
and Shahbazi (2010) suggested that using a cesium formate based dril- The complexity involved in barite settlement in an inclined annulus
ling fluid or polymer-coated ultra-fine barite for drilling HPHT wells will under the impact of a rotating inner pipe, motivated this study. This work
prevent barite particles from sagging. Zamora (2009) investigated the is aimed at modeling and simulation of barite settlement in the annulus
mechanisms, measurement methods used, and barite sag mitigation space. The effects of inner pipe rotation and pipe inclination on the barite
under static and dynamic conditions showing the degree of complexity settling, particle volume fraction and pressure profile of the fluid have
been studied. The model presented in this work enables us to better
describe barite settlement. Inaddition, particle settlement in annular
space with rotational inner shaft has been modeled and simulated. To
achieve this, a Eulerian two phase model for different pipe angle as well
as different rotational speeds has been used.
2.1. Barite settlement in a vertical well with stationary inner shaft and fluid
at rest condition
477
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
to represent the well and the drill pipe (see Fig. 2). The outer cylinder in which
representing the well has an inner diameter (ID) equal to 0.1 m. The inner
cylinder with an outer diameter (OD) equal to 0.05 m has been consid- ΣF ¼ Fbody þ Fsurface (5)
ered to represent the drill pipe. It is assumed that during shut-in period, In here the body force is gravitational force and surface forces consist
both the drill pipe and fluid are at stationary condition. of normal and shear forces.
By applying mass conservation to a control volume of an annulus The body force can simply be shown by the Buoyant force:
system, we will have
!
πdp3
_ z mΔt
mΔt _ zþΔz ¼ MtþΔt Mt (1) FBuoyant ¼ ρp ρf (6a)
6
h
Δt: ρp vset A z ρp vset A zþdz ¼ ρp Adz tþΔt ρp Adz t (1a) The normal force acting on a particle is equal to the magnitude of
forced caused by pressure and can be determined using the
following equation:
∂ðCvset AÞ ∂ðCAÞ
¼ (2)
∂Z ∂t Fp ¼ ∫ ps !
n dS (6b)
cs
∂C ∂C ∂vset
¼ vset þ C (3) The shear stress force is resulted from fluid acting on a settling par-
∂t ∂Z ∂z
ticle and is equal to:
Fτ ¼ ∫ τ:!
where C is the volumetric concentration.
n dS (6c)
As it is shown in Eq. (3), to find the barite concentration profile, it is cs
critical to find the isovelocity of the particles. To simplify the problem,
first we have analyzed only a single particle system. The forces applied on Under the assumption of the stationary condition and by assuming a
such particle include gravitational, buoyancy and drag forces (see Fig. 3). constant settlement velocity, the Stokes drag force FD will represents the
Therefore, the settling velocity can be found by considering a simple magnitude of the shear force applied on the particle surface. The Stokes
force balance and their impact on particle velocity. Such force balance drag force cab be found using the following equation:
can be written as follows: 1
FD ¼ ρf CD Ap u2 (7)
dv 2
ΣF ¼ m (4)
dt where vset is the particle settling velocity, CD is drag force coefficient, Ap is
the particle project area, ρf is the particle density, dp is the particle
diameter and (u) is slip velocity of particle influid.
Assuming constant settling velocity dv
dt ¼ 0 and by combining
the buoyant force and shear force equations, the terminal velocity can be
determined using the following equation:
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u4 ðρ ρ Þ
u3 s f gd
t ρf p
vt ¼ (8)
CD
Due to the impact of flow regime on drag coefficient, the solid particle
drag coefficient or CD will be a function of particle Reynolds number such
that for low Reynolds number the CD can be estimated by using the
following equation:
Fig. 3. Force balance on a solid particle submerged in a quiescent liquid.
478
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
24 Using the same argument, the Reynolds number for Bingham model is
CD ¼ (9)
Rep generated as follows:
Conservation equations
In order to model a complex flow of solid particles in side an annular
domain, volume fraction, conservation of mass and momentum have
been applied.
Table 1
Values for the Richardson and Zaki index m.
Fig. 4. Drag coefficient as a function of particle Reynolds number (Duan et al., 2015). Rep 500 2.4
479
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
dq vq dp vp
Fvm ¼ 0:5ρq Cp (22)
dt dt
The interphase force for a mixture with n phases is shown below:
X
n
! X
n
R pq ¼ Kpq !vp!
vq (23)
p¼1 p¼1
Cs ρs f
Ksl ¼ (24)
τs
Fig. 5. Barite bed characteristics (Nguyen et al., 2009).
τq ¼ Cq μq ∇!
vqþ!
¼
Flift ¼ 0:5ρq Cp !vp!
T
v q ∇v_q (21) vq (29)
In multiphase flow systems, when the discontinuous phase (particles) For the solid phase we have:
is accelerated through a continuous phase (fluid), the fluid experiences
T 2 ¼
τ s ¼ Cs μs ∇! v s þ Cs λs μs ∇:!
vsþ!
acceleration due to the work that is done on particles. In another word, ¼
vs I (30)
the difference in the acceleration of continuous and discontinues phases 3
leads to a virtual mass force. For further detail on this subject, please refer ¼
to Crowe et al. (2011) and Kolev (2007). Such force is defined by the where superscript T indicates the transpose and I is the identity tensor. λs
following equation: is the bulk viscosity of the solids phase that can be obtained by
the following:
480
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Cs
g0 ¼ 1 (32)
CsMax Fluid-fluid momentum equations
The conservation of momentum for a continuous fluid phase l can be
In this study CsMax is equal to 0.67. In addition, ess is the restitution written as follows:
coefficient and is equal to 0.9.
∂
θs is the granular temperature and it is proportional to kinetic energy ðCl ρl !
v l Þ þ ∇:ðCl ρl !
v l Þ ¼ Cl ∇p þ ∇:τl þ Cl ρl ! g þ ðKsl ð! vs!
¼
v lÞ
of the fluctuating particle motion. ∂t
! ! !
μq is the viscosity of liquid phase and μs is the shear viscosity of solids þ m_ sl v sl Þ þ Cl ρl F lift;l þ F vm;l
phase that contains collisional, kinetic and the optional fric-
(39)
tional elements:
where ! g is the gravitational acceleration.
μs ¼ μskin þ μsfr þ μscol (33)
Similar to the approach followed for the continuous phase, the con-
Fig. 6 shows the three terms of solid shear viscosity including colli- servation of momentum for the discontinuous solid phases can be
sional, kinetic, and frictional in the inclined bed such as inclined pipe determined by:
or well.
∂
ðCs ρs !
v s Þ þ ∇:ðCs ρs !
v s!
v s Þ ¼ Cs ∇p ∇ps þ ∇:τs þ Cs ρs !
¼
The collisional viscosity is modeled as: g
∂t
12 ! ! !
þ ðKls ð v l v s Þ þ m_ ls v ls Þ
4 θs
μs;col ¼ Cs ρs ds g0 ð1 þ ess Þ (34) ! !
5 π þ Cs ρs F lift;s þ F vm;s (40)
Following Syamlal et al. (1993) kinetic viscosity term can be
where Ps is the Sth solids pressure, Kls ¼ Ksl is the momentum exchange
modeled by:
coefficient between phases and n is the total number of phases.
pffiffiffiffiffiffiffi
Coupling the above partial differential equations leads to a nonlinear
Cs ds ρs θs π 2
μs;kin ¼ 1 þ ð1 þ ess Þð3ess 1ÞCs g0 (35) system of equations. Such system can only be solved by means of nu-
6ð3 ess Þ 5
merical methods.
In addition, 1987 In additionthe frictional viscosity term is modeled
using model proposed by Schaeffer (1987) 3. Numerical solution
481
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 8. The Pål Skalle experimental set up (a) and model developed by CFD (b). Figure b shows the volume concentration profile along new pipe test.
482
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
In order to further validate the CFD model, the effects of pipe rotation,
inclination angle, and rheological properties of base fluid have been
experimentally examined and compared with the results of the numerical
model. The experimental approach and validation process is presented in
next section.
Table 3
The base fluid contents.
Contents Volume
Weight ,%
Water 40%vol.
Oil 45%vol.
Primary emulsifier 5%vol.
Secondary emulsifier 10%vol.
Lime 3%wt.
Fig. 9. The comparison of between result of Pål Skalle experiment and CFD model. CMC 5%wt.
483
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Table 4
The fluids and barite particles properties.
Table 5
Properties of fluid and particle used in simulation.
Fig. 11. a. Schematic of designed annular flow loop. b. Laboratory apparatus used to create the flow loop system.
484
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 13. Comparison of the results of simulation with experiments on the fluid pressure change at the bottom of annular section during static test.
in the annulus section was recorded continuously with time passes by and then changed to linear trend. This behavior is based on this fact that
using multiple pressure transducers. This practice is based on this idea the particles settling velocity at early time is an accelerated movement,
that changes in the fluid density as a result of particles settling will cause and it is not possible to reach the constant velocity at the initial times.
pressure changes in the fluid column. Using the pressure transducers at However as shown, there is generally good agreement between experi-
the different segments of the annular section, it was possible to monitor mental data and CFD results for different inclination angles Fig. 15 shows
the fluid density and rate of barite settlement with time passes. A high the effects of rotational speed on barite settlement rate in the annular
accuracy pressure transducer was used in different segments of the section under no axial flow conditions and the inclination angle of
annular section of the flow loop in order to record the fluid pressure 30 and 45 . This figure obviously shows that the low rotational speed of
change. Moreover, as indicated in Fig. 11a it was possible to get the 70 rpm has not significant impact on barite settlement rate.
sample from the fluid in different segments of the annular section using
the sample cap. Therefore, it was possible to analyze the fluid density or 5. Modeling results and sensitivity analyses
particles concentration in the fluid sample.
As stated before, the fluid rheology of tests conducted in this work This section provides the results of sensitivity analyses conducted on
differ from fluid used by Skalle et al. (1999). Fig. 13 shows the experi- settlement in both vertical and inclined wells. The tests conducted on a
mental results of pressure reduction percent for different inclination single particle assume a static drill pipe (shut-in period only). However,
angle. Comparison between Figs. 9 and 13 indicate that the rate of barite under multi-particle system, both static and dynamic drill pipes were
settlement decreases as fluid yield point increases. Under no axial flow modeled numerically. All simulations were performed by using a model
conditions, the barite settlement was analyzed through measuring the geometry indicated in Fig. 2.
pressure reduction in the different segment of annular section. Fig. 14
shows the results of CFD simulation and experimental data on pressure 5.1. Single particle in a vertical well during shut-in period
reduction at the location of bottom pressure transducer port. As shown, at
early times, the rate of pressure change was similar to exponential trend Governing Eq. (3) in combination with Eq. (17) is used to model a
Fig. 14. Experimental data for pressure reduction percent under static condition.
485
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 15. Pressure reduction vs. time in annular section with rotational speed of 70 rpm, a: inclination angle ¼ 30 and b: inclination angle ¼ 45 .
simple vertical well. In order to numerically model the barite settling where ρavg is the average density of fluid in at each time step and Cb;v is
phenomenon, finite deference method is used. An element of a well is the average barite volume concentration:
simulated with a length of 1 m. The simulation domain has been divided PN t
into ten equal size grids (see Fig. 16). i¼k Cb;i
C b;v ¼ (42)
The simulation results depicted in Fig. 17 shows the variation of N ðk 1Þ
concentration profile in each grid at different time. The results suggest
that concentration variation with respect to time follows an S-shape In this equation, N is total number of grids and k is the grid number in
profile. This simply means that rate of concentration variation with time which pressure is measured at. The pressure at each grid can be calcu-
for each cell is a function of neighboring cells above and below. Fig. 18 lated using the following formula:
shows the variation of barite concentration along pipe for different time
1
durations. As shown the barite concentration at the bottom increases Pti ¼ ρavg g N i dL ði ¼ 1; 2; … NÞ (43)
2
with time due to barite settlement.
For each grid in the model and at each time step, the average fluid
density can be calculated by determining the average concentration: P0i
PN ¼ (44)
Pti
ρavg ¼ ρfluid 1 C b;v þ Cb;v ρbarite (41)
here Pti is pressure in the grid i at time t > t0 , PN and P0i are the normalized
pressure and initial pressure in the grid i, respectively.
Figs. 19 and 20 show the pressure and normalized pressure variations
at the different levels of fluid column in the pipe. As shown, the base fluid
hydrostatic pressure has been observed at the top of the column sooner
than that experienced in the bottom of the column. This is due to the fact
that the top element is acting just as a source and therefore barite con-
centration in such element reduces with time. However, other elements
experience a balance between particles entering such elements and
particles that leave them.
An example of such calculation is provided for grid 2 of the vertical
pipe model:
if t ¼ t0 ⇒C b;v ¼ C0
Similar to the initial condition, the model is bounded by the following
boundary condition:
if x ¼ L⇒ C b;v ¼ 0
For any other grid and at time t > t0 , the average concentration can be
determined as follows:
486
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 17. Variation of concentration profile in each grid (after complete settling).
1
P12 ¼ ρavg g 10 2 0:1 ¼ ρavg g 0:85
2
P02
PN ¼
P12
5.2. Eulerian two phase model on an inclined pipe with rotational inner
shaft
487
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 21. Effect of inner pipe rotation on the particle volume fraction at the bottom of the pipe (z ¼ 0), Pipe angle ¼ 45 and t ¼ 600sec.
Fig. 25 has been generated for an annulus with inclination angle of 45 in magnitude and direction. The force exerted by vortex flow on settled
with stationary inner shaft. Since pipe angle is considered as an influ- particles will push them away from the settlement bed and will lead to
ential parameter, contours of barite volume fraction along annulus for an resuspension of such fine particles. As a result, the thickness of particle
inclination angle of 60 with static inner shaft has been analyzed as well. bed on the lower side of the annulus will reduce. Fig. 29 shows the
Such results are shown in Fig. 26. vectors of velocity and the created vortices due to inner pipe rotation. As
Considering the importance of inner shaft rotation speed on barite sag shown, there is an upward-down ward motion of fluids element due to
settlement, contours of barite concentration along inclined annulus, the impact of barite settlement. When the barite particles starts to settling in
results of for two different rotation speeds under similar inclination angle a same direction, the elements will move in opposite direction and
have been compared in Fig. 27. In this figure, barite concentration pro- upward-downward motion is occurred.
files for a case with inner shaft rotation of 100 RPM have been compared To better understand how influential the vortex flow is, velocity of the
with those of a pipe with stationary inner shaft. As shown by increasing particles in the settlement layer next to the annulus wall has been stud-
the inner pipe rotation to 100 rpm the particle settlement rate is reduced. ied. Fig. 30 depicts the axial velocity of particles of such layer at two
The contours of barite axial velocity (called z velocity here) are also different inner shaft rotation speeds. In other words, this figure shows
shown in Fig. 28. As it is shown in this figure and due to the rotation of that the volume of particles settled on the lower wall of the pipe can be
the inner shaft, a counter current flow of particles inside the pipe can be reduced by increasing the pipe rotation speed.
observed. Considering the size of barite particles, the vortexes created by In order to characterize the barite settling under different conditions
the rotating central pipe significantly impacts their settling velocity both such as different inner pipe rotations and different pipe angles, a barite
488
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 22. The contours of barite volume fraction at the bottom of the pipe for different time duration. Pipe angle ¼ 75 and inner pipe rotation ¼ 50 rpm.
sag coefficient needs to be introduced to better understand the settlement from the upper elements and quantity of barite left such elements to the
phenomenon inside the annulus of a well under no axial flow condition. lower elements and therefore have not been used to define this index.
The “top” and “bottom” sections are used as reference locations due to One may use this concept to define a dynamic index for the whole length
the fact that maximum difference in barite concentration can be found of the pipe by slightly adjusting the sag coefficient defined above. In
between such locations. In another word these locations are used as petroleum industry both approaches have been used.
reference locations to define a. The dimensionless index sag coefficient is The top and bottom sections in simulation study are defined by using
defined as is follows: the cells in such locations. This means that the first set of cells represent
the “top” reference location and same applies for “bottom” reference
Cb;bottom Cb;top Cb;top location. Let's look at the top as an example, the top section is consisted of
Sag Coefficient ¼ ¼ 1 (45)
Cb;bottom Cb;bottom number of cells with known specific volume and barite concentrations.
Using such information the average barite concentration for this section
where Cb is volume fraction of barite in the fluid. has been calculated. The top and bottom sets of cells are illustrate
The index will result in a value between 0 and 1 in which represents in Fig. 31.
the degree of sagging. If the system is homogeneous with regards to the Therefore, by increasing the settlement of particles and because the
barite concentration, the index will be equal to zero and under extreme barite concentration at the top of the inclined bed approaches to zero, the
case in which there is no barite left in the top section, the index will be Sag coefficient will approach to 1. Similarly, when particles are well
equal to one. In addition it should be noted that barite concentration in suspended and settling does not occur, Ctop is equal to Cbottom and sag
other locations depends on the amount of barite entered such elements coefficient becomes zero.
489
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 23. Effect of rotational speed on particle volume fraction at the bottom of annular section under no axial flow conditions.
Fig. 24. The contours of barite volume fraction at the bottom of inclined annular section under no axial flow condition and rotation speed of 70 rpm.
Fig. 32 shows the variation of sag coefficient with time for different pronounced. At certain time, the cohesive force becomes the dominant
pipe rotation. As shown by increasing the pipe rotation, sag coefficient is factor and the settlement occurs at a reduced and constant rate. The time
decreased which means the settling rate is decreased. Also, the results span in which the first settlement pattern can be observed strongly de-
reveal that settlement rate will decrease as inner shaft rotation increases. pends on the inner shaft's speed. At higher rotations, the vortex flow will
However, at high rotation speeds (higher than 60 rpm), the difference in not allow the particles to create a bed. Therefore the particles will be
sag coefficient and rate of settlement is negligible. In addition, this figure resuspended, leading to lower settlement rates.
shows that two major settlement regimes can be observed. At early time Fig. 33 shows the volumetric concentration profile along pipe test at
the drag force and gravity force are dominant forces which control the pipe angle of 45 at two different rotations (stationary inner shaft and
magnitude of settlement and its rate. However, by time the impact of inner pipe rotation of 100 rpm). As shown the volumetric concentration
rolling mechanism and as a result the cohesive force becomes more is measured at the both side of pipe test in two different locations above
490
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 25. The front of constant barite concentration along pipe test with angle of 45 and no inner rotation shaft.
Fig. 26. The contours of barite volume fraction along the 1 m concentric annuli for pipe angle ¼ 60 and inner pipe rotation speed ¼ 0 rpm.
and below the pipe center line (y ¼ 0.049 and y ¼ 0.049 m). The results Fig. 34 shows the volumetric concentration along a vertical pipe after one
show that for each rotation speed, the volumetric barite concentration hour. It can be seen that in a system with stationary inner pipe, only a
below the centerline is more than its concentration above the centerline. small change in concentration has been observed. However, at high
491
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 27. The contours of barite volume fraction along the 1 m concentric annuli during 1 h with Pipe angle ¼ 45 . Case a with inner pipe rotation ¼ 100 rpm and case b without inner pipe
rotation.
Fig. 28. Contours of barite axial (z) velocity (m/s), t ¼ 600sec, pipe angle ¼ 45 and inner pipe rotation ¼ 100 rpm indicating counter current particle flow.
rotation speed particle settlement did not occur. 6. Summary and conclusion
Figs. 35 and 36 show the variation of sag coefficient with time for
different pipe angles at inner pipe rotations of 0 rpm and 50 rpm. As In this paper, static barite particle settling phenomenon in an annular
shown the maximum sag coefficient was observed at pipe angle of 45 . By space similar to drilling well was studied using single particle settling
decreasing or increasing the inclination angle the sag coefficient will modeling approach and Eulerian two-phase flow model. The results of
reduce. This result is consistent with other researcher's experimental theoretical modeling and CFD simulation of barite settling phenomenon
works (Dye et al., 1999, 2003, 2006; Skalle et al., 1999). is presented here. The impact of different forces and mechanisms on
Figs. 37 and 38 show the effect of base fluid viscosity on barite set- settling process have been investigated in detail. The particle-particle and
tlement phenomenon for a Newtonian fluid. As shown, the rate of barite particle-fluid interactions have been taken into account while modeling
settlement decreases by increasing the fluid viscosity. Figs. 39 and 40 this process based on three-dimensional Eulerian modeling. In order to
show the effect of initial barite content on rate of settlement for a base validate the numerical model, the experimental analyses have been
fluid with viscosity of 15 cp. As shown by increasing the barite content in conducted. The numerical model predictions indicated a good agreement
the base fluid, the sag coefficient slightly decreases due to the hindering with both available data in the literature and experimental results. Next,
effect of particles (eq. (17)). the sensitivity analyses were performed to evaluate the impact of inner
492
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 29. Particle velocity vectors. Effect of pipe rotation on particle velocity direction.
1.2
30 rpm 0 rpm
1
0.8
Pipe length(m)
0.6
0.4
0.2
0
0 -0.005 -0.01 -0.015 -0.02 -0.025 -0.03
Particle Axial Velocity (m/s)
Fig. 30. Effect of pipe rotation on axial velocity of particles in first settlement layer on the lower pipe wall (at y ¼ 0.04994 m). Pipe angle ¼ 45 .
Fig. 31. Reference locations on numerical grids: a) The upper section cells b) location of both top and bottom reference locations.
pipe rotation, deviation angle, barite content, and base fluid rheology on study showed that barite settlement phenomenon strongly depends on
the rate of barite settlement. operational parameters such as drill string rotation speed, well inclina-
In this study a dimensionless sag coefficient has been defined to tion angle as well as fluid rheology. The results showed that the settling
characterize the degree of sagging under different conditions with time. rate of particles decreases by increasing the inner pipe rotation. However,
Combination of such parameter with detailed numerical and experi- after a critical rotation speed the rate of settlement is constant and
mental analysis enable us to better characterize and predict the barite increasing the rotation speed has minimal impact on settlement. The
settlement behavior under various situations during drilling process. This critical rotation speed depends on fluid viscosity and rheology type and
493
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 34. Volume concentration profile along pipe test at pipe angle ¼ 90 during 1 h.
Fig. 32. Variation of sag coefficient with time for different inner pipe rotation. Pipe
angle ¼ 45 .
fluid properties. The simulation results also showed that the vortices
created by inner pipe rotation along the well significantly impact the
sagging process through resuspension of particles. In addition, the effect
of well orientation and inclination on barite settlement has been
analyzed and investigated. The results showed that the barite particles
settling rate maximizes at the pipe angle range of 45 –60 and then
decreased with increasing pipe angle to 90 . In vertical orientation (90
pipe) and under no axial flow condition settlement of multiple particles
affected by the amount of the drag force caused by a microflow between
particles as well as collision between them which have a hindering effect,
gravitational forces and buoyancy force. This hindering effect increases
by increasing particles or increasing the concentration of particles. In
addition, in order for these particles to settle, they have to displace same
volume of water. In vertical orientation, combination of these effects
reduces the settling velocity in comparison with that of a deviated well.
Fig. 35. The comparison of Sag coefficient variations for different pipe angles and no
When it comes to deviated pipes, the resultant force applied to each
inner pipe rotation condition.
particle leads to a vertical settlement of particles on the side of the pipe.
Unlike the settlement of particles in vertical pipe, the resistance faced by
water placed between particle and pipe body is reduced. The particles
settled on the side of the pipe is simply “rolling down” on the pipe body decreases as the barite content increases which indicate the hindering
and therefore the hindering effect due to presence of other particles will effect on the rate of settlement. The effect of fluid rheology on barite
significantly reduce. Such movement mechanism is known as rolling settlement was also investigated. The results show that the rate of barite
mechanism. The mechanism of settling in a deviated pipe is referred to as settlement decreases as fluid yield point increases. The results also show
boycott settling. that increasing the base fluid viscosity decrease the rate of barite
Analysis of sag coefficient showed that the sag coefficient slightly settlement.
0.5
0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2
Length(m)
Fig. 33. Volumetric concentration profile along pipe test at pipe angle ¼ 45 .
494
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
Fig. 39. The effect of initial barite content on barite concentration distribution along the
pipe test during 30 min.
Fig. 36. The comparison of Sag coefficient variations for different pipe angles and inner
pipe rotation ¼ 50 rpm.
Fig. 40. The comparison of Sag coefficient variations with initial barite content during
30 min.
Fig. 37. The effect of base fluid viscosity on barite concentration distribution along the
pipe test during 30 min.
Nomenclature
495
H. Movahedi et al. Journal of Petroleum Science and Engineering 161 (2018) 476–496
References Madlener, K., Ciezki, H., 2005. Theoretical Investigation of the Flow Behavior of Gelled
Fuels of the Extended Herschel Bulkley Type.
Mishra, P., Ein-Mozaffari, F., 2017. Using computational fluid dynamics to analyze the
Akhshik, S., Behzad, M., Rajabi, M., 2016. CFD-DEM simulation of the hole cleaning
performance of the Maxblend impeller in solid-liquid mixing operations. Int. J.
process in a deviated well drilling: the effects of particle shape. Particuology 25,
Multiph. Flow 91, 194–207.
72–82.
Mme, U., Skalle, P., 2012. CFD calculations of cuttings transport through drilling annuli at
Amighi, M., Shahbazi, K., 2010. Effective Ways to Avoid Barite Sag and Technologies to
various angles. Int. J. Petrol. Sci. Technol. 6 (2), 129–141.
Predict Sag in HPHT and Deviated Wells, SPE Deep Gas Conference and Exhibition.
Mondal, S., Wu, C.-H., Sharma, M.M., 2016. Coupled CFD-DEM simulation of
Society of Petroleum Engineers.
hydrodynamic bridging at constrictions. Int. J. Multiph. Flow 84, 245–263.
Bern, P., et al., 1998. Barite Sag: Measurement, Modelling and Management, IADC/SPE
Movahedi, H., Farahani, M.V., Jamshidi, S., 2017. Application of Hydrated Basil Seeds
Asia Pacific Drilling Technology. Society of Petroleum Engineers.
(HBS) as the herbal fiber on hole cleaning and filtration control. J. Petrol. Sci. Eng.
Boycott, A., 1920. Sedimentation of blood corpuscles. Nature 104, 532.
152, 212–228.
Crowe, C.T., Schwarzkopf, J.D., Sommerfeld, M., Tsuji, Y., 2011. Multiphase Flows with
Nguyen, T.C., Yu, M., Takach, N., 2009. Predicting Dynamic Barite Sag in Newtonian Oil-
Droplets and Particles. CRC press.
based Drilling Fluids, SPE Annual Technical Conference and Exhibition. Society of
Dritselis, C., 2017. Numerical study of particle deposition in a turbulent channel flow
Petroleum Engineers.
with transverse roughness elements on one wall. Int. J. Multiph. Flow 91, 1–18.
Nguyen, T.C., Miska, S., Saasen, A., Maxey, J., 2014. Using Taguchi and ANOVA methods
Duan, Z., He, B., Duan, Y., 2015. Sphere drag and heat transfer. Scientific reports 5.
to study the combined effects of drilling parameters on dynamic barite sag. J. Petrol.
Dye, W., Hemphill, T., Gusler, W., Mullen, G., 1999. Correlation of Ultra-low Shear Rate
Sci. Eng. 121, 126–133 .
Viscosity and Dynamic Barite Sag in Invert-emulsion Drilling Fluids, SPE Annual
Omland, T.H., 2009. Particle Settling in Non-newtonian Drilling Fluids.
Technical Conference, pp. 543–553.
Parvizinia, A., Ahmed, R.M., Osisanya, S.O., 2011. Experimental study on the
Dye, W., Mullen, G., Gusler, W., 2003. Drilling Processes: the Other Half of the Barite Sag
phenomenon of barite sag. In: International Petroleum Technology Conference.
Equation, SPE Asia Pacific Oil and Gas Conference and Exhibition. Society of
International Petroleum Technology Conference.
Petroleum Engineers.
Richardson, J.F., Zaki, W.N., 1954. Sedimentation and fluidisation: Part I. Trans. Inst.
Dye, W.M., Mullen, G.A., Gusler, W.J., 2006. Field-proven Technology to Manage
Chem. Eng. 32, 35–53.
Dynamic Barite Sag, IADC/SPE Drilling Conference. Society of Petroleum Engineers.
Saasen, A., Liu, D., Marken, C., 1995. Prediction of Barite Sag Potential of Drilling Fluids
Elgaddafi, R., Ahmed, R., George, M., Growcock, F., 2012. Settling behavior of spherical
from Rheological Measurements, SPE/IADC Drilling Conference. Society of
particles in fiber-containing drilling fluids. J. Petrol. Sci. Eng. 84, 20–28.
Petroleum Engineers.
Hajidavalloo, E., Sadeghi-Behbahani-Zadeh, M., Shekari, Y., 2013. Simulation of
Schaeffer, D.G., 1987. Instability in the evolution equations describing incompressible
gas–solid two-phase flow in the annulus of drilling well. Chem. Eng. Res. Des. 91 (3),
granular flow. J. Differ. Equ. 66 (1), 19–50.
477–484.
Shahi, S., Kuru, E., 2016 Aug 10. Experimental Investigation of the settling velocity of
Hanson, P., Trigg Jr., T., Rachal, G., Zamora, M., 1990. Investigation of Barite" Sag" in
spherical particles in Power-law fluids using particle image shadowgraph technique.
Weighted Drilling Fluids in Highly Deviated Wells, SPE Annual Technical Conference
Int. J. Miner. Process. 153, 60–65.
and Exhibition. Society of Petroleum Engineers.
Skalle, P., Backe, K., Lyomov, S., Sveen, J., 1999. Barite segregation in inclined boreholes.
Hashemian, Y., et al., 2014. Experimental study and modelling of barite sag in annular
J. Can. Petrol. Technol. 38 (13).
flow. J. Can. Petrol. Technol. 53 (06), 365–374.
Syamlal, M., O'Brien, T.J., 1989. Computer simulation of bubbles in a fluidized bed. In:
Fluent Inc, 2006. Fluent User's Guide. Fluent Inc., Lebanon, NH, USA, Version 6.3.26. .
AIChE Symp. Ser, pp. 22–31.
Kamyab, M., Rasouli, V., 2016. Experimental and numerical simulation of cuttings
Syamlal, M., Rogers, W., O'Brien, T.J., 1993. MFIX Documentation: Theory Guide.
transportation in coiled tubing drilling. J. Nat. Gas Sci. Eng. 29, 284–302.
Technical Note DOE/METC-95/1013 and NTIS/DE95000031. National Energy
Kaushal, D., Thinglas, T., Tomita, Y., Kuchii, S., Tsukamoto, H., 2012. CFD modeling for
Technology Laboratory, Department of Energy.
pipeline flow of fine particles at high concentration. Int. J. Multiph. Flow 43, 85–100.
Turton, R., Levenspiel, O., 1986. A short note on the drag correlation for spheres. Powder
Kolev, N.I., 2007. Multiphase Flow Dynamics 3: Turbulence, Gas Absorption and Release,
Technol. 47 (1), 83–86.
Diesel Fuel Properties, v. 3. .
Wang, Z.-m., Long, Z.-h., 2010. Study on three-layer unsteady model of cuttings transport
Lun, C., Savage, S.B., Jeffrey, D., Chepurniy, N., 1984. Kinetic theories for granular flow:
for extended-reach well. J. Petrol. Sci. Eng. 73 (1), 171–180.
inelastic particles in Couette flow and slightly inelastic particles in a general flow
Zamora, M., 2009. Mechanisms, measurement and mitigation of barite sag. In: Offshore
field. J. fluid Mech. 140, 223–256.
Mediterranean Conference and Exhibition. Offshore Mediterranean Conference.
496