Effect of Yield Power Law Fluid Rheological Properties on Cuttings Transport in Eccentric Horizontal Narrow Annulus

Narrow annular drilling such as casing-while-drilling technique is gaining popularity due to its ability to mitigate nonproductive time during oil and gas drilling operations. However, very little is known about the flow dynamics in narrow annular drilling. In this study, the Eulerian-Eulerian two-fluid model was used to examine the influence of Yield Power Law fluid rheological properties on cuttings transport in eccentric horizontal narrow annulus. The flow was assumed as fully developed, laminar, and transient state. The present simulation model was validated against experimental data, where a mean percent error of −1.2% was recorded. Results revealed an increase in the radial distribution of cuttings transport velocity in the wide annular region as the consistency index, K, and the flow behavior index, n, increase. Nonetheless, increasing the yield stress, τ o , had insignificant effect on the cuttings transport velocity. Three-dimensional profiles showed how cuttings preferred to travel in less resistant flow area, whereas cuttings concentration builds up in the narrow annular region. Furthermore, annular frictional pressure losses also increased as K, n, and


Introduction
Drilling fluid rheological properties are important parameters which contribute to effective hole cleaning.Adari et al. [1] indicated that drilling fluid rheological properties are highly influential on cuttings transport; hence, care must be taken in predicting the optimum parameters to enhance better hole cleaning.The Herschel-Bulkley (HB) viscosity model, also called Yield Power Law (YPL) model, is known to correlate better to the drilling fluid viscosity curve than most other rheological models [2,3].
In literature, most experimental studies on drilling fluids observed that the drilling fluid rheological curves (rheogram) conformed best to that of YPL fluid model.Ahmed and Miska [4] conducted extensive experimental study with polymerbased fluids in concentric and eccentric horizontal annuli without the presence of cuttings.The fluid rheological properties were accurately characterized using YPL model.The effects of inner pipe rotation, eccentricity, and flow rates on frictional pressure losses were analyzed.Other experimental studies [5][6][7][8] also accurately fitted their fluid rheological parameters using the Yield Power Law model.The authors proposed pressure loss equations for YPL fluid flow in pipes and annuli which is an integration of analytical, semianalytical, and empirical equations for laminar, transitional, and turbulent flows.Taghipour et al. [9] were among the first researchers to use YPL fluids to perform cuttings transport experimental study in inclined and horizontal annulus of diameter ratio,  = 0.50.They measured annular pressure losses, solids bed height, and drill string torque.
Several experimental and numerical studies have also evaluated the effects of drilling fluid rheological properties on cuttings transport.Cho et al. [10] observed that a decrease in flow behavior index, , resulted in a decrease of stationary bed, whereas moving bed layer increases.Adari et al. [1] also observed that the increase in the ratio of flow behavior index to consistency index (/) reduced cuttings bed height.In addition, the drilling fluid yield stress (YS) and plastic

Materials and Methods
The inhomogeneous (Eulerian-Eulerian) two-fluid model in ANSYS CFX 14.0, where both the liquid and solid phases are considered interpenetrating continua, is adopted in this study.Other models such as Eulerian-Lagrangian model, however, simulate the solid phase as a discrete phase and allows for particle tracking.The Eulerian-Eulerian model is preferred to the Eulerian-Lagrangian model due to its ability to handle high solid volume fractions.Furthermore, it accounts for particle-particle interaction and includes turbulence automatically.A drawback of this model is, however, the need for complex closure relations.The commercial software package ANSYS CFX 14.0 consists of the following five (5) specialized components: (a) DesignModeler for building the geometries, (b) CFX Mesh for mesh generation (c) CFX-Pre for flow model definition, (d) CFX-Solver for solving the governing equations, and (e) CFX-Post for analyzing the results.The following continuity and momentum equations representing the two-phase flow model are described for the sake of brevity.

Governing Equations.
The solid-liquid flow is assumed as (a) isothermal and (b) laminar and transient state.

Continuity Equations.
The governing continuity equations for both liquid and solid phases could be expressed, respectively, as [29] where the solid and liquid phase volume fractions sum up to unity as 2.1.2.Momentum Equations.The forces acting on each phase and interphase momentum transfer term that models the interaction between each phase are given below [29].
For liquid phase, Similarly, for solid phase,  For densely distributed solid particles, where the solid volume fraction k  < 0.2, the drag coefficient model by Wen-Yu [30] may be utilized.This model is modified and implemented in ANSYS CFX 14.0 to ensure the correct limiting behavior in the inertial regime as where   Re  = k   Re  and  Re  =   |  −   |  /  .For large solid volume fraction, k  > 0.2, the drag model by Gidaspow [31] may be used with the interphase drag force per unit volume defined as 2.2.2.Lift Force Model.For spherical solid particles, ANSYS CFX employs the lift force model by Saffman [32,33] as The authors correlated the lift force for low Reynolds number past a spherical solid particle, where    = 6.46 and 0 ≤  Re  ≤  Re  ≤ 1.For higher range of solid particle Reynolds number, Saffman's correlation was generalized by Mei and Klausner [34] as where  = 0.5( Re  / Re  ), and The suspension viscosity is expressed in terms of relative viscosity,   , as The relative viscosity of more concentrated suspension with particle-particle interactions is given by [35]   = 1 + 2.5k  + 10.05k 2  + 0.00273 exp (16.6k  ) .(13) The first two terms of (13) represent Einstein [36] equation.
Equations (11) to (13) are coded into CFX library as expression language to compute the solid viscosity.

Carrier Fluid Viscosity.
There is singularity problem associated with the classical YPL viscosity model at vanishing shear rate.To alleviate this, the proposed YPL viscosity function by Mendes and Dutra [37] is implemented in the present CFD study as This viscosity function is numerically stable and devoid of discontinuity.Equation ( 14) is currently not available in the library of variables of ANSYS CFX 14.0; hence, it is modeled as expression language and coupled with the momentum equations.

Boundary Conditions and Meshing.
At the inlet, a mass flow rate is specified, while a zero-gauge pressure is specified at the outlet boundary.At the pipe walls, different boundary conditions were used for both liquid and solids.The usual no-slip condition was imposed at the walls for the liquid phase, while for the solid phase, the free-slip condition was assumed at the walls to prevent the solid phase from adhering to the walls.This is consistent with real flow behavior of solid particles flowing near a solid boundary.The solid volume fraction in the domain was specified at the beginning of each simulation to correspond to the desired solid loading.
The annular 3D geometry of diameter ratio  = 0.90 was meshed into structured hexahedral grids of approximately 2.4 × 10 5 elements.Grid independent study of numerical results of pressure losses were carried out until variations in results were insignificant.In this study, the optimum grid sizes used for the radial, circumferential, and axial directions are 20() × 120() × 100().Figure 1 shows a section of the 2D and 3D annular mesh, while Figure 2 depicts the results from the grid independence study.

Description of Simulation Study.
The simulation of the two-phase solid-liquid flow was set up in three dimensions using ANSYS CFX 14.0 with the transport equations solved using CFX-Solver.The geometry dimensions, fluid rheological properties, solid properties, and operating parameters are presented in Table 1.To achieve a fully developed flow, the annular length, , must be longer than the hydrodynamic entrance length  ℎ of the flow.In single-phase Newtonian fluids flowing in pipes, the hydrodynamic entrance length is presented as [38]  ℎ = 0.05 Re () .
However, for a two-phase flow in annular gap with a non-Newtonian fluid, such expression as in (15) does not exist in literature.As a rule of thumb, the author adopted (15) by replacing the pipe diameter  with a hydraulic diameter The CFX-Solver is based on a finite volume method in which the flow equations are integrated over each control volume.The advection scheme is set to high resolution to satisfy both accuracy and boundedness, where the blend factor, Γ, is computed locally to be close to 1 without resulting in nonphysical values.The algorithm by Rhie and Chow [39] is used to solve the pressure-velocity coupling, since it can overcome pressure-velocity oscillations.
It should be noted that the complexity of the transport equations could not permit numerical convergence under steady state.However, all simulations were run in transient state.It is usually recommended that simulation of such steady state nature should first be run under transient state when it is difficult to attain convergence [23].A very small time step of 2.0 × 10 −5 s was used in all simulations to help the solution converge; hence, the solutions finally converged when the convergence criterion was met at a root mean square (RMS) value of 10 −4 .An average of 200 time steps was required to achieve convergence with 1-10 loop iterations at each time step.

Validation of Simulation Model.
There is scarcity of data for cuttings transport study using YPL fluids.The most recent and only experimental cuttings transport study is conducted by Taghipour et al., where annular pressure losses were measured for cuttings-YPL fluid flow in inclined (30 ∘ ) eccentric wellbore and are used to validate the two-phase CFD model.To validate the present model, the author simulated the experimental condition as presented by Taghipour et al. [9].Table 2 presents the experimental setup with operational and rheological parameters.It is worth noting that only a section of the experimental flow loop length of 12 m was simulated.A simulated length of 100 × (  −   ) was chosen based on the calculated hydrodynamic length in (15) to ensure a fully developed flow as well as saving computational time.The flow was assumed as laminar and isothermal, while the drag and  lift forces were modelled using the models by Wen-Yu and Saffman as presented by ( 5) and ( 8), respectively.Figure 3 depicts the comparison between experimental and simulated pressure loss data.The comparison shows good agreement with a mean percentage error (MPE) between the calculated and measured pressure loss data of −1.2%, thus confirming the validity of the current model setup.

Results and Discussion
The simulation results reported here include the effects of yield stress, consistency index, and flow behavior index on cuttings transport velocity in YPL fluid in eccentric horizontal narrow annulus.Measurements of the radial distributions of the cuttings velocity profiles were taken along sectors A-A and B-B or at the positions 0 ∘ and 180 ∘ (see Figure 1), representing the wide and narrow regions of the eccentric wellbore, respectively.To generate the cuttings velocity profiles in CFX-Post, a plane was first created parallel to the symmetry of the wellbore.Secondly, a polyline was also created with a boundary intersection method, where the boundary list of wellbore outlet was intersected with the plane.Finally, a chart was chosen where the velocity profiles were plotted using the polyline.
Studies have shown that YPL fluids used in the field have a wide range of values of rheological properties.These rheological properties could include the following range of values [13]: yield stress,   = 0.048-50.06Pa, consistency index,  = 0.043-10.27Pa⋅s n , and flow behavior index,  = 0.314-0.978.The range of rheological parameters used in this study is within the aforementioned range of values usually encountered during drilling operation using YPL fluids.yield stress in the wide annular region.This is an indication of better carrying capacity of high yield stress fluids.Nonetheless, this improvement in cuttings transport velocity is less significant especially in the core region as yield stress increases.In contrast, the narrow annular region provides greater flow resistance to high yield stress fluids, thus the least to transport cuttings, as shown in Figure 4(b).In the wide annular region (sector A-A), increasing   from 2 to 8 Pa led to ∼0.8% increment in the maximum cuttings velocity, whereas ∼5.8% reduction in cuttings velocity occurred as   increases from 2 to 8 Pa in the narrow annular region (sector B-B).

Consistency Index, 𝐾, Effect on Cuttings Transport Velocity.
The consistency index of a drilling fluid is a rheological property related to the cohesion of the individual particles of the fluid, its ability to deform, and its resistance to flow.Figures 5(a) and 5(b) illustrate the consistency index's effect on cuttings transport velocity in both wide and narrow annular regions, respectively.It is observed in Figure 5(a), that is, in the wide annular region, that when all other parameters are kept constant, increasing  results in an increase in cuttings velocity.
It should be noted that  is directly proportional to the fluid's effective viscosity.Therefore, cuttings will travel less in low  fluids with low effective viscosity as the cuttings tend to settle faster at the bottom of the annulus due to gravity.Whereas in high  fluids, having high effective viscosity, cuttings will suspend for a longer period and, hence, travel faster and farther.
On the contrary, due to the flow restriction induced by the narrow annular region, cuttings travelling in high  fluids had the least velocity as opposed to those travelling in low  fluids (see Figure 5(b)).Investigations on the range of  values presented in this study show that increasing  from 1.7 to 6.3 Pa⋅s n leads to ∼11.0% increment in the maximum cuttings velocity in the wide annular region (sector A-A), whereas ∼37.7% reduction in cuttings velocity was recorded as  increases from 1.7 to 6.3 Pa⋅s n in the narrow annular region (sector B-B).

Flow Behavior Index, 𝑛, Effect on Cuttings Transport
Velocity.The flow behavior index, , is a measure of a fluid's shear-thinning performance in both pseudoplastic and yield stress fluids; that is, apparent viscosity decreases with increasing shear rate.Fluids with higher shear-thinning properties,  < 1, exhibit a plug (wider flat) flow region in the central core.Figures 6(a) and 6(b) demonstrate the effect of  on cuttings transport velocity in both wide and narrow annular regions, respectively.In the wide annular region, as shown in Figure 6(a), cuttings transport velocity in the fluid domain decreases as the fluid becomes more shear-thinning, especially in the core region.In the vicinity of the outer and inner pipe boundaries, high shear-thinning fluids exhibit a flat radial cuttings velocity profile.This increases the high velocity zone towards the boundaries with improved cuttings transport.Within the study range, increasing  from 0.31 to 0.75 led to increment of ∼36.3% in the maximum cuttings velocity within the wide annular region (sector A-A) as shown in Figure 6(a).In the narrow annular region (sector B-B), cuttings travel faster in the radial distribution in low  fluids compared to high  fluids as shown in Figure 6(b).It is also observed that increasing  values from 0.31 to 0.75 led to a significant reduction in maximum cutting velocity by ∼84.6%.

Three-Dimensional Cuttings Velocity and Concentration
Profiles.Figure 7 presents the three-dimensional distributions of cuttings travelling in the annular geometry.Skewed cuttings velocity and concentration profiles were observed in the eccentric annulus as shown in Figures 7(a) and 7(b), respectively.These irregular distributions were due to the pipe-hole eccentricity which modified the flow by creating a high velocity zone in the wide annular region and a low or no velocity zone in the narrow annular region.The narrow region restricted the movement of the trapped cuttings and hence showed a very low velocity (see Figure 7(a)).Furthermore, high accumulation of cuttings bed can be observed in the narrow region of the annulus due to the effects of gravity and eccentricity as shown in Figure 7  3.5.Annular Pressure Losses.Accurate estimation of annular pressure losses is very vital when designing drilling hydraulic programs, particularly the equivalent circulating densities (ECD) required for efficient transport of drilled cuttings from the wellbore to the surface.Figure 8 examines the influence of , , and   on annular pressure loss.Figure 8(a) indicates that as the fluid becomes more shear-thinning, that is,  < 1, the pressure loss decreases.There is a gradual rise in pressure as loss from  = 0.31 to 0.51; however, a dramatic percentage increase of 316.4% is computed as  increased from 0.51 to 0.75.Furthermore, Figure 8(b) also reveals an increase in pressure loss as  increases from 1.7 to 6.3 Pa⋅s n with approximately 132.6% increase in pressure loss.Last but not the least, increase in   from 2 to 8 Pa had the least influence on pressure with nearly 13.6% increase.

Conclusions
A study on the effect of rheological parameters on cuttings transport velocity in YPL fluid flowing in eccentric narrow horizontal annulus is analyzed using Eulerian-Eulerian twofluid CFD model.The proposed viscosity model for YPL fluid  by Mendes and Dutra [37] was adopted since it is capable of avoiding numerical singularity difficulties at vanishing shear rate.For two-phase cuttings-YPL fluid flow, CFD calculated pressure loss data were validated using recent but rare experimental data from literature.The present model was used to analyze the effects of yield stress, consistency index, and flow behavior index of the carrier YPL fluid on cuttings transport velocity.
Within the context of this study, the radial distributions of cuttings transport velocity were observed in both wide and narrow annular regions.The increase in yield stress of the carrier fluid from 2 to 8 Pa did not have much influence on the cuttings transport velocity in the wide region especially in the core region.However, there was much improvement in the cuttings transport velocity at the vicinity of the walls, indicating less cuttings bed accumulation.In the narrow gap, cuttings travelled much faster in low yield stress fluids as a result of their low resistance to flow compared to high yield stress fluids.
The study also revealed that carrier fluids with high consistency index value enhanced more cuttings transport especially in the wide annular region due to the fluid's cuttings lifting ability.At the vicinity of the walls, cuttings travelled faster in low consistency index fluids due to the fluid's tendency to high shearing and, hence, resulted in less cuttings bed formation.In the narrow region, however, cuttings travelled faster in low consistency index fluids due to their less resistance to flow.Increasing the flow behavior index of the carrier fluid also showed much improvement in the cuttings transport velocity, especially in the core region of the wide annular gap.Meanwhile, carrier fluid with low flow behavior index transported cuttings better at the vicinity of the walls, an indication of less formation of cuttings bed.A reverse trend was, however, observed in the narrow annular gap, where there was significant cuttings transport in the carrier fluid with low flow behavior index.Three-dimensional flow distribution profiles have shown the actual dynamics of cuttings travelling in the eccentric annulus, where most of the cuttings were inclined to travel in the wide margin with less stress.
The YPL fluid model has been shown to fit much better rheological data of drilling fluids in the oil and gas industry compared to the Bingham plastic and power law fluid models.Furthermore, YPL fluids with high yield stress values reduce convective heat loss in most especially high temperature wellbores, thus possessing the very rheological properties which are important for an insulating fluid to perform well.It is noteworthy that YPL fluids have viscosities that increase significantly as shear-strain rate diminishes; hence, by increasing the viscosity of the fluid, the drilling engineer can gain partial control over convection heat loss.More importantly, YPL fluids tend to have relatively low viscosity at high shear rates (shear-thinning), making them easier to place initially, to bleed off pressure that may build up in the annulus equipped with venting capability, and to displace the drilling fluid in the event of well intervention.This study provides a guide to the drilling engineer on the selection of YPL fluid rheological properties which would enhance efficient transport of drilled cuttings in narrow horizontal wellbores.
It was further observed that the rheology of YPL fluids has significant effect on the annular pressure losses.Proper caution must be taken in selecting the fluid rheology to ensure efficient cuttings transport while maintaining a bottomhole pressure which will not fracture the formation.

Figure 3 :
Figure 3: Experimental and simulation comparison of pressure loss data.

Figure 4 :
Figure 4: Effect of yield stress on cuttings transport velocity: (a) wide annular region and (b) narrow annular region.

Figure 6 :
Figure 6: Effect of flow behavior index on cuttings transport velocity: (a) wide region and (b) narrow region.
ℎ =   −   .It should be noted that a much longer annular length would only result in a longer simulation run.