1. Introduction
Fluid losses during drilling operations, a significant challenge in the field, are primarily influenced by mud filtration in permeable areas (at the wellbore and fracture walls) and flow through fractured zones. This phenomenon not only results in fluid losses but also affects other crucial aspects of drilling operations. For instance, changes in filtrate concentration can alter resistivity log readings, impacting data accuracy [1,2,3,4,5]. Furthermore, the rate of mud filtrate advance is known to cause formation damage and contribute to the growth of mud cake, complicating drilling operations [6,7,8,9,10].
The issues analyzed contribute to non-productive time, which costs hundreds of millions of dollars annually and is defined as the time during which the drilling rate is reduced or stopped. States that worldwide, around $20 billion is invested annually in drilling, of which 15% corresponds to NPT. On the other hand, some studies carried out in Colombia in the Castilla field, have shown that lost circulation is the main cause occupying 31% of the NPT [11,12].
The dynamics of mud filtration have been extensively studied both through laboratory experiments and mathematical modeling. Experimental studies focusing on fluid displacement in rock samples have shed light on the impact of mud filtration on formation damage, allowing us to determine the reduction in permeability due to the effect of fluid filtration invasion [13,14,15]. Complementing laboratories are numerical approaches, which have been instrumental in modeling the filtration process and mud cake growth, employing both linear (laboratory) and radial (physical models) coordinates [16,17,18,19]. These models enable the analysis of various interaction scenarios between the drilling fluid and the permeable zone. Proposals in the literature are based on mass balances to determine fluid filtration and fluid concentration changes, studying advective, concentration-gradient, and pressure-gradient transport. Filtrate invasion affects variables such as fluid cake growth, formation permeability, and fluid saturation at the wellbore. This research seeks to couple all of the mentioned phenomena and characterize filtration in a realistic manner.
In naturally fractured zones, the analysis of drilling fluid losses has primarily been based on fluid flow in channels, linking fluid losses to fracture geometry. Some researchers have concentrated on characterizing natural fractures concerning losses observed in the field [20,21], while others have estimated fluid losses based on fracture geometry [22,23,24,25,26,27,28,29,30,31]. Some authors have linked the permeable zone with the fracture zone without including the growth of the fluid cake [24,26,32,33] (Table 1). In addition, sometimes fractures connected by the drilling fluid can reactivate and propagate, which can increase losses, so it is important to strengthen fluid flow and geomechanically couplings (Solid-liquid) [34] and monitor rock properties and stresses over time [35]. The models that analyze fluid losses in fractures have in common that they use the momentum balance equation to propose an expression that relates shear stress to the pressure gradient and then equates it to the rheological model that represents the behavior of the drilling fluid under different shear rates.
In this context, our article introduces a coupled numerical model designed to widely represent the interaction between drilling fluid and the wellbore. This model encompasses the phenomena of mud filtration in permeable zones and fluid flow in natural fractures. For the permeable zone, the mud filtration rate, changes in filtrate concentration, permeability reduction, and mud cake growth were calculated based on the law of continuity and the theory of convection-dispersion. In natural fracture zones, the flow was evaluated, by determining the velocity profile, flow rate, and accumulated volume of losses, using the motion equation balance and incorporating both Newtonian and non-Newtonian rheological models.
This novel modeling approach enhances the precision and understanding of the interaction between drilling fluid and the wellbore. It provides crucial support for decision-making during well-planning and subsequent operational phases.
2. Developed Methodology
Our methodology to analyze the interaction between the drilling fluid and the wellbore through mathematical modeling is divided into two distinct parts, each focusing on different types of zones: permeable and fractured. This structured approach allows for a comprehensive and detailed analysis of each zone type. Figure 1 illustrates the workflow of our methodology.
A.. Analysis of Permeable Zones
In permeable zones, our focus is threefold: Growth of the Mud Cake: Here, we assess the rate and pattern of mud cake formation on the wellbore. This analysis helps in understanding the filtration dynamics and the associated impact on wellbore stability. Filtration Rate: We evaluate the filtration rate of drilling fluid, considering pressure changes at the wellbore. This involves analyzing the pressure differentials and their influence on fluid dynamics and filtration rates. Changes in Filtrate Concentration: This part of the analysis looks at how the filtrate concentration varies over time and its implications for wellbore integrity and drilling efficiency. B.. Analysis of Fractured Zones
The study of fractured zones involves examining the flow velocity of the drilling fluid, incorporating different rheological models. This segment is further subdivided into: Power Law Velocity: Analysis under this model focuses on the flow behavior of drilling fluids following a Power Law fluid model. Bingham Plastic Velocity: Here, we explore the flow characteristics of drilling fluids modeled as Bingham Plastic fluids, which is crucial for understanding non-Newtonian fluid behavior in fractured zones. Modified Power Law Velocity: This analysis involves a modified approach to the Power Law model, offering insights into complex fluid dynamics in fractured geological structures.
2.1. Numerical Modeling—Permeable Zone
The interaction between the drilling fluid and the permeable zone (Figure 2) depends on the movement of fluids due to the pressure gradient and diffusion due to the concentration gradient. These factors affect the growth of the cake and present the following assumptions.
The drilling mud is water-based, and the filtrate exists as a single phase.
The properties of the mud cake, such as density and porosity, remain unchanged over time, after the growth of the fluid cake stabilizes over time.
The formation is homogeneous, isotropic, and characterized by isothermal flow (This assumption applies to each node in the radial direction, as you move to the next node in the vertical direction, the layer characteristics may change.).
Formation properties, are considered constant. (This assumption applies to each node in the radial direction, as you move to the next node in the vertical direction, the layer characteristics may change.).
Formation fluid displacement leads to changes in saturation and relative permeabilities.
The bottom-hole pressure (pw) is maintained constant (Applies if the density of the drilling fluid remains constant, otherwise the value is updated).
The thickness of the cake varies in “r” as a function of time and in the “z” direction as a function of layer properties.
In Section 2.1.1 the proposed mathematical model is derived, based on the assumptions.
2.1.1. Growth of the Mud Cake
The growth of the mud cake is fundamentally driven by the dynamics of particle deposition and erosion in the drilling mud, as outlined in previous studies [17,18,36,37]. It is premised on the assumption that particles suspended in the drilling mud water-based do not displace within the porous medium but only deposit on the wellbore, as expressed in Equation (1).
(1)
the growth of the fluid cake is determined by the balance between the rates of particle deposition and erosion, as shown in Equation (2):(2)
the rate of deposition is defined by Equation (3), and the erosion rate is characterized by Equation (4):(3)
(4)
Literature suggests that the constant is derived from experiments involving latex suspensions (0.04 µm) filtered through glass beads (diameter 0.397 mm and porosity 0.35), using a specific velocity expression to model adhesion and suspension. This results in being defined as [15,38]. Likewise, the constant has been successfully implemented in previous research [36,37], which used core flooding data to validate your model.
To establish an expression for mass flow per unit of time, Equation (3) is multiplied by the cross-sectional area of filtration, and Equation (4) by the porous volume of the mud cake. These are then substituted into Equation (2), yielding Equation (5):
(5)
Further, by differentiating Equation (1) with respect to time for a radial system and equating it to Equation (5), we derive a new model to determine the growth of the mud cake at the wellbore, as shown in Equation (6). To apply the equation, it is assumed that the properties of the fluid cake such as density and porosity do not change over time
(6)
where is the initial concentration, is the fluid cake density, is the wellbore radius, is the fluid cake thickness and is the erosion constant and the terms and are the porosity and the filtration rate. The derivation of the equation is seen in Appendix A.2.1.2. Filtration Rate
In our model, although the bottom hole pressure () is considered constant, fluid invasion at the wellbore leads to a pressure increase. To represent this, the continuity equation was used (Equation (7)), which is solved using forward implicit finite differences. To apply the equation, it is assumed that the formation is homogeneous and isotropic, therefore the intrinsic petrophysical properties in the radial direction do not change and as it is assumed isothermal, the viscosity of the filtrate remains constant for the node of study.
(7)
where is the pressure variation between well and reservoir, is the compressibility of the rock, is the viscosity of the fluid, is the porosity, is the permeability of the rock, and is the pressure variation over time. To analyze the filtration rate in the mud cake, Darcy’s law is applied in radial coordinates, as shown in Equation (8):(8)
where is the permeability of the fluid cake, is the difference between the bottom hole pressure and the reservoir pressure, is the radius of the well, is the thickness of the fluid cake and is the radius of the reservoir affectation.2.1.3. Concentration of Filtrate
To examine changes in filtrate concentration, we adopt previous research [16,18] which employs the Crank-Nicolson method to approximate and solve the relevant differential equations. The transport of mud particles is analyzed through two primary mechanisms:
Advection Transport: This involves the movement of dissolved solids along with the flowing water, also known as advective transport. The quantity of solute transported depends on its concentration in the water and the volume of water flowing, as expressed in Equation (9) [39]:
(9)
Concentration Gradient Transport: Here, solutes in water migrate from areas of higher concentration to those of lower concentration, a process termed diffusion. This occurs even in static fluid and is governed by Fick’s first law. The mass of fluid diffusing is proportional to the concentration gradient, described by Equation (10) [39]:
(10)
Advective transport is represented by Equation (11) and dispersive transport by Equation (12). The total mass flow is the sum of these two mechanisms, indicated in Equation (13):
(11)
(12)
(13)
Including Equations (11)–(13) in a mass balance and considering fluid filtration leads to Equation (14):
(14)
where is the variation in filtrate concentration over time, is the diffusion constant, and the terms and are the porosity and the filtration rate defined previously. In scenarios with no filtration, movement is solely due to concentration changes, resulting in Equation (15):(15)
For correcting the diffusion coefficient in porous media, we use the correlation [40], with exponent values for consolidated formations varying between ϵ = 1.8 to ϵ = 2.
Section 2.1.1, Section 2.1.2 and Section 2.1.3 can also be developed in linear coordinates to validate the phenomenon using laboratory data, as shown in Equations (16) to (18):
(16)
(17)
(18)
For pressure variation, we implement the continuity equation in linear coordinates (Equation (19)):
(19)
Flow velocity is estimated using the sample pressure [9], as shown in Equation (20):
(20)
The concentration of solids at the sample face and for the diffusion coefficient, previous research is used [13,16,17,39], as shown in Equation (21):
(21)
Since it is assumed that the displacement of formation fluid leads to changes in saturation and relative permeabilities, it is necessary to calculate their variation in time. Permeability correction due to filtrate invasion is calculated for each time step using the mud saturation equation proposed in previous research [16], as modified in Equation (22):
(22)
where is the filtrate concentration at each distance , is the saturation that can change, the difference between the maximum water saturation and the irreducible water saturation, and represents the increase in filtrate saturation of the fluid that invades the formation. Relative permeability is then calculated using Equation (23) [41]:(23)
For studies lacking documented permeability of the cake, an estimation can be made using the Equation (24) [42]:
(24)
where is the effective particle radius, is the permeability decay constant, γ is the volume fraction of particles in the filter cake, and .2.2. Numerical Modeling—Fractured Zone
In this part of our study, we focus on the velocity profile of the drilling fluid within fractures (refer to Figure 3). The behavior of the fluid depends on its rheology and the pressure generated by the fluid column. This aspect is approached through the balance of the motion equation [43]. As a result of this balance, a relationship is established between the shear stress perpendicular to the “Z” axis in the “r” direction and the pressure gradient, as shown in Equation (25) (Appendix B).
(25)
By equating Equation (24) with the rheological model of the fluid, we derive an equation for calculating the velocity profile. This is summarized in Table 2.
Table 2Velocity Equation Based on Rheological Model.
Rheology | Equation | |
---|---|---|
Newtonian | (26) | |
Power Law | (27) | |
Bingham Plastic | (28) | |
(29) | ||
Plug Flow | and | |
Modified Power Law | (30) | |
(31) |
is the consistency index, is the fracture width, is the flow behavior index, is the pressure variation between the fracture and the formation, is the initial shear stress used in modified power law and refers to the fracture width variation.
3. Application of the Numerical Model and Results
The solution of the equations proposed is carried out by means of an algorithm developed in MATLAB-2022. Differential equations whose analytic solution is difficult to implement, are approximated with finite differences using the implicit Crank-Nicolson method (Equation (32)).
(32)
after obtaining the approximations, the matrix is built, for example, the structure for pore pressure would be as follows:(33)
where the constants would group the following properties(34)
As for the boundary conditions, there and and for the meshing, similar geometry to the literature models was constructed in order to compare and validate, having a constant radius and angle increment.
3.1. Interaction with Permeable Zone
For the permeable zone, Ref. [14] analyzed three types of fluids: fresh water, flocculated, and dispersed (Table 3 and Table 4).
To validate our model using data from Table 2 and Table 3, an algorithm was designed in MATLAB, incorporating Equations (16)–(24) in linear coordinates [14].
The accumulated volume of mud filtrate reported in the literature [14] and calculated with our proposed model, exhibits an excellent fit for the three types of fluids, as depicted in Figure 4. However, for cake growth, our model shows a divergence after the initial 50 min, stabilizing at a constant value. This variation appears because the proposed model incorporates the erosion term, unlike the continuous increase reported in the literature [14], which did not consider this factor. The effectiveness of the fluid filter cake growth model was verified by comparing the erosion trend of a physical model published in the literature [44] with the erosion trend generated by Equation (18). In the physical model, upon reaching the maximum thickness, erosion intensifies, decreasing its thickness to a constant thickness as shown by the proposed numerical model.
To validate the model for radial coordinates, the data in Table 5 were used [45]. An algorithm was developed in MATLAB, incorporating Equations (6)–(15) and (22)–(24) in radial coordinates. The scenario reported in the literature [45], focusing on the impact of drill string eccentricity and rotational speed on filtrate cake thickness, was used for validation. Their findings indicate that rotational speed does not significantly influence cake size, which corroborates the validity of our base case comparison.
The results for filtration rate and cake growth are shown in Figure 5. Our findings align with literature-reported mud filtration rates, but the cake growth is lower due to the inclusion of the erosion term in our model. By including the erosion term, the fluid cake will grow to a maximum value, and then erosion will intensify until it stabilizes at a constant value.
Moreover, the drilling fluid concentration and pore pressure at the wellbore, crucial inputs for Equation (6), were calculated as a function of well radius. This analysis was essential for estimating the invasion radius of the filtrate, which reaches up to 0.5 m, as illustrated in Figure 6. This information is vital for correcting resistivity logs and calculating treatment volumes to mitigate formation damage.
3.2. Interaction with Natural Fractures
To evaluate the interaction between drilling fluid and natural fractures, we applied Equations (26)–(31) suited for different rheological models. The fluid velocity within fractures was calculated based on fracture width and pressure differential. To validate the model with the modified Power Law (Herschel–Bulkley, Equations (30)–(31)) and for the Power Law (Equation (27)). The relevant data for validation is presented in Table 6, and the results of each case are represented in Figure 7.
The validation outcomes with Power Law and modified Power Law models provide credibility to other models based on similar theoretical principles. This suggests that the Bingham Plastic and Newtonian models could also be utilized, contingent upon the specific rheology of the fluid.
Additionally, to verify the model’s applicability, a comparison was made with a Computational Fluid Dynamics (CFD) study focusing on a horizontal fracture [46]. The data used for this comparison are detailed in Table 7, and the results for scenarios at 200, 500, 800, and 1000 psi are illustrated in Figure 8. The rheology applied in this comparative analysis was Herschel–Bulkley.
Finally, to confirm the effectiveness of the proposed methodology and before analyzing field data, the results of Equation (27) are compared with the data generated from a physical model [47], which consists of measuring the time it takes for a liter of drilling fluid to flow through a fracture with known dimensions. The fracture and the comparison of results are shown in Figure 9.
4. Field Application
After implementing the numerical model and validating it with data from the literature, the next step is to carry out an application in a Colombian field, estimating the rate of drilling fluid losses in fracture zones and filtrate in the permeable zone. The field of study is the Llanos Orientales basin. For the application, the GeoDrill-Flow tool was developed using MATLAB’s APPDESIGNER, which is an interactive development environment that facilitates the design of applications and the programming of their behavior. The information required for the application is loss events, fracture geometry (inclination and azimuth), geomechanical model, petrophysical model, mud weight, bit position, and mud rheology, as shown in Figure 10, Figure 11 and Figure 12. In addition to the information mentioned, the data in Table 8 is required [48].
With the field information and the model developed in Section 2, we proceed to estimate the growth of the fluid cake, the filtration rate, the velocity of fluid losses, the permeability of the fluid cake, the pressure changes, and the concentration changes as observed in Figure 13.
The variables calculated in Figure 13 allow us to determine the pressure variation within the fracture (Figure 14a) and total fluid loss rate (Figure 14b), including the flow in the fracture and the filtrate in the permeable zone, which is compared with field data.
Figure 14b is applied to a drill bit diameter 12 ¼ in to a depth of around 6500 ft. From that section on, the diameter changes, so another modeling would have to be performed. This figure shows an excellent fit between field data and modeled data, which highlights the effectiveness of the coupled model implemented.
5. Discussion
The developed model facilitates a comprehensive analysis of the interaction between drilling fluid and the permeable zone, accounting for variables such as petrophysics, rheology, and operational conditions. Additionally, the model assesses the impact of the mud’s exposure time on the wellbore, which is crucial for understanding cake growth dynamics. This aspect is particularly relevant in field operations, where the well remains exposed to the mud until casing placement. For prolonged exposure periods, our model determines the maximum cake growth value, as illustrated in Figure 15a. Furthermore, recording the filtrate invasion radius (Figure 15b) is essential to comprehend the extent of the invaded zone and to strategize cleaning treatment volumes, therefore minimizing formation damage.
In Figure 16, the maximum thickness observed in the mud cake is approximately 3 mm. This value is lower than the thickness reported in previous research [45], which can be attributed to the erosion factor included in our model. Regarding the invasion radius, it extends to 0.6 m, 20% deeper than indicated in Figure 6. This expansion is due to the longer study duration, suggesting that prolonged mud exposure leads to a larger affected radius, independent of pressure differential. The filtrate movement can occur due to concentration changes.
Another critical aspect is the size of the mud particles, which significantly influences permeability reduction, mud cake growth, and filtrate invasion. For instance, Figure 16 demonstrates the permeability reduction as a function of particle size. Smaller particles result in lower permeability and thinner mud cakes, which are advantageous in drilling operations.
As indicated in Figure 8, an increase in fracture radius leads to a decrease in fluid velocity. This observation can be combined with other variables to mitigate fluid losses. For example, a similar trend is evident when pressure decreases (Figure 17a) and with varying initial shear stress levels (Figure 17b). Higher initial stress slows down fluid movement, suggesting that fluids with high yield stress can help decrease losses. Similarly, smaller particles improve the mud cake quality, reducing filtrate invasion.
6. Conclusions
This study introduces a comprehensive methodology to improve the interaction between drilling fluid and the wellbore. It effectively estimates fluid velocity within fractures in naturally fractured reservoirs, and examines the influence of rheology, pressure, and operational conditions on this interaction.
Key contributions of the proposed model include: Analysis of Fluid-Wellbore Interactions: The model offers insights into the dynamics between drilling fluid and wellbore, analyzing the permeable and fracture zone, and factoring in petrophysical properties, fluid rheology, and operational conditions. Impact of Mud Exposure Time: It highlights the significance of mud exposure time on the wellbore, a critical aspect in field operations where the well is exposed until casing placement. This feature of the model is vital for understanding the implications of prolonged mud exposure. Quantification of Drilling Dynamics: Through the coupled model, we quantitatively evaluate various factors such as filtration velocity, filtrate concentration, filter cake formation, and invasion radius. This approach incorporates time and distance-based analysis and applies convection-dispersion theory to explore different mass transfer mechanisms in permeable zones. Informing Operational Decisions: The understanding gathered from studying mud behavior in both permeable and fractured zones equips drilling and mud personnel with valuable information. This knowledge is crucial for optimizing mud formulations and strategies, ultimately aiding in the reduction of non-productive time (NPT) associated with wellbore stability issues.
In summary, the proposed model enhances the understanding of complex fluid dynamics within wellbores, thereby contributing to more efficient and effective drilling operations, particularly in challenging geological settings.
Conceptualization, methodology, software, D.A.V.S.; validation, D.A.V.S., D.C.M. and Z.H.C.; formal analysis, D.A.V.S.; investigation, D.A.V.S. and Z.H.C.; resources, D.A.V.S.; data curation, D.C.M.; writing—original draft preparation, D.A.V.S.; writing—review and editing, Z.H.C. and G.E.R.; visualization, Z.H.C. and G.E.R.; supervision, Z.H.C. and D.C.M.; project administration, D.A.V.S.; funding acquisition, (no funding was received). All authors have read and agreed to the published version of the manuscript.
The data previously published in the master’s thesis
Thanks to the Well Stability Research Group (GIEP) for their advice.
The authors declare no conflicts of interest.
The following nomenclature is used in this paper, with units provided for clarity:
Term | Nomenclature | Units * |
Cake volume | ||
Cake mass | ||
Cake porosity | ||
Deposition constant | ||
Erosion constant | ||
Deposition rate | ||
Erosion rate | ||
Filtrate concentration | ||
Filtration velocity | ||
A | Cross-sectional area | |
Mass flow rate | ||
Cake thickness | ||
Rock compressibility | ||
Fluid viscosity | ||
Fluid permeability | ||
Formation porosity | ||
Bottom hole and reservoir pressures | ||
Cake permeability | ||
Well radius | ||
Cake thickness | ||
Diffusion coefficient | ||
Sample length | ||
Effective particle radius | mm | |
Permeability decay constant | Dimensionless | |
γ is the volume fraction of particles in the filter cake | Dimensionless | |
Consistency index | ||
Flow behavior index | Dimensionless | |
Initial yield stress | ||
Fracture width | ||
* Note: For equations in linear coordinates, the cgs system is used. |
Footnotes
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.
Figure 1 Workflow to Analyze the Interaction Between Drilling Fluid and Wellbore.
Figure 2 Radial Geometry with Side and Top View of the interaction between the drilling fluid and the permeable zone.
Figure 3 Velocity Profile of Drilling Fluid Within the Fracture.
Figure 4 Comparison of Filtration Rate and Cake Growth for Different Fluids; (a,b) for flocculated; (c,d) for dispersed and (e,f) for fresh water [
Figure 5 (a) Filtration rate through the fluid cake and permeable zone, (b) Fluid cake growth as a function of time [
Figure 6 (a) Pore Pressure Distribution in the wellbore, (b) Filtrate Concentration in the Wellbore.
Figure 7 Velocity Profile for Different Rheologies. (a) Herschel–Bulkley [
Figure 8 Fluid Velocity Within Fractures as a Function of Fracture Radius for different pressure scenarios (a) Pressure differential of 200 psi (b) Pressure differential of 500 psi (c) Pressure differential of 800 psi (d) Pressure differential of 1000 psi [
Figure 9 (a) Physical model (b) comparison of results between Mathematical model and Physical model.
Figure 10 Colombian basin study well events. (a) Drilling fluid loss events in the well (b) Cumulative volume of drilling fluid losses, (c) Drilling fluid with loss control material implemented in the well to control fluid losses (d) Time the event lasted.
Figure 11 Fracture geometry of a study well in the Colombian basin (a) Width of the fractures recorded with image registration, (b) Inclination of the fractures, (c) Azimuth of the fractures.
Figure 12 Geomechanical model of a study well in the Colombian basin (a) Mechanical properties of the formations in the study well (b) Internal friction angle (c) Pore pressure (d) Principal stresses in the study area.
Figure 13 Interaction between drilling fluid and the wellbore (a) Pore pressure in wellbore (b) Filtrate concentration in wellbore (c) Filtrate velocity (d) Fluid cake permeability (e) Flow velocity within the fracture (f) Cake thickness over time.
Figure 14 (a) Pressure distribution within the fracture. (b) Modeled data vs. field data.
Figure 15 Impact of Time on (a) Cake Growth, (b) Filtrate Invasion [
Figure 16 Influence of Particle Size on Mud Cake Permeability.
Figure 17 Fluid Velocity Profile Within the Fracture for (a) Different Pressure Values, (b) Different Shear Stress Values.
The table shows the most important aspects of each model.
Author | Rheological Model | Fracture Flow | Filtered | Fluid Cake Growth |
---|---|---|---|---|
[ | Newtonian | Yes | No | No |
[ | Bingham Plastic | Yes | No | No |
[ | Power law | Yes | No | No |
[ | Power law Bingham Plastic | Yes | No | No |
[ | Power law | Yes | Yes | No |
[ | Power law | Yes | Yes | No |
[ | Herschel–Bulkley | Yes | No | No |
[ | Herschel–Bulkley | Yes | No | No |
[ | Herschel–Bulkley | Yes | No | No |
[ | Truncated power-law | Yes | No | No |
[ | Power law-Herschel–Bulkley | Yes | Yes | No |
[ | Power law-Bingham Plastic | Yes | Yes | No |
Laboratory Properties to Validate the Proposed Model.
Variable | Value | Unit | Reference |
---|---|---|---|
Circulating Pressure | 0.6 | MP | [ |
Diameter | 2.54 | cm | [ |
Length | 20.32 | cm | [ |
Permeability | 0.102 | μm2 | [ |
Porosity | 21.4 | % | [ |
Confining Pressure | 1000 | Psi | [ |
Density | 1.03 | g/cm3 | [ |
Cake Permeability ranges between | 1 × 10−5–1 × 10−7 | μm2 | [ |
Dispersion Coefficient | 0.0063 | cm2/s | [ |
Erosion Constant | 4.35 × 10−4 | 1/s | [ |
Rheological Properties [
Fluid | K (Pa.sn) | n |
---|---|---|
Fresh Water | 0.8 | 0.319 |
Dispersed | 0.7 | 0.313 |
Flocculated | 1.7 | 0.235 |
Data for Model Validation in Radial Coordinates.
Variable | Value | Unit | Reference |
---|---|---|---|
Pressure | 200 | kPa | [ |
Radius of interest | 3 | m | |
Permeability | 9.8 × 10−13 | m2 | [ |
Porosity | 30 | % | |
Fluid Density | 1345 | kg/m3 | [ |
Initial Concentration assumed similar to density | 1000 | kg/m3 | [ |
Cake Permeability | correlation | | [ |
Dispersion Coefficient | 5 × 10−9 | m2/s | [ |
Erosion Constant | 4.35 × 10−4 | 1/s | [ |
String Radius | 0.127 | m | [ |
Well Radius | 0.254 | m | [ |
Particle Radius | 5 | | [ |
Effective Particle Radius in Cake | 0.5 | | [ |
Power Law Exponent (n) | 0.7 | [ | |
Consistency Index | 0.1238 | Pa.sn | [ |
Filtrate Viscosity | 0.0014 | Pa.s | |
Volume Fraction of Solids in Cake | 0.6 | [ | |
Decay Constant | 0.001 | 1/s | [ |
Particle Density | 4500 | kg/m3 | [ |
Data for Validating Velocity Profile for Different Rheologies.
Rheology | K (Pa.sn) | n | τo (Pa.s) | W (m) | ΔP | Reference |
---|---|---|---|---|---|---|
Herschel–Bulkley | 0.1 | 0.8 | 10 | 0.001 | 6 psi | [ |
Herschel–Bulkley | 0.2231 | 0.7 | 2.5 | 0.001 | 40 kPa/m | [ |
Power Law | 0.005 | 0.3 | N/A | 0.001 | 1-5-75 Pa/m | [ |
Data for CFD Validation [
Variable | Value | Unit |
---|---|---|
Well Radius (rw) | 0.11 | m |
Fracture Radius (rf) | 5 | m |
Fracture Width (w) | 8.8 × 10−4 | m |
Pressure (DP) | 200 | Psi |
Power Law Exponent (n) | 0.94 | |
Consistency Index (k) | 0.04 | Pa.sn |
Yield Stress (τo) | 4.022 | Pa |
Additional data required for field application [
Variable | Value | Unit |
---|---|---|
Porosity ( | 0.2 | |
Permeability (k) | 1 × 10−13 | m2 |
Internal Diameter (ID) | 12.5 | in |
Particle diameter (Dp) | 5 × 10−7 | m |
Power Law Exponent (n) | 0.411 | |
Consistency Index (k) | 0.2877 | Pa*sn |
Yield Stress (τo) | 5 | Pas |
Appendix A. Fluid Cake Growth
Starting from a mass balance we have
Solving and simplifying
Applying the derivative and factoring,
Remembering the Equation (5),
Equating the Equations (A6) and (A7), the following expression is obtained,
Knowing that the volume of the cake is
Simplifying
Appendix B. Balance of the Equation of Motion
Taking the radial component (Equation (A11))
Knowing that,
Integrating we have
Considering
1. Semmelbeck, M.E.; Holditch, S.A. Effects of mud-filtrate invasion on the interpretation of induction logs. SPE Form. Eval.; 1988; 3, pp. 386-392. [DOI: https://dx.doi.org/10.2118/14491-PA]
2. Lane, H.R. Numerical simulation of mud filtrate invasion and dissipation. Proceedings of the SPWLA 34th Annual Logging Symposium; Calgary, AB, Canada, 13–16 June 1993.
3. Wu, J.; Torres-Verdín, C.; Sepehrnoori, K.; Delshad, M. Numerical simulation of mud filtrate invasion in deviated wells. Proceedings of the SPE Annual Technical Conference and Exhibition; Orleans, LA, USA, 30 September–3 October 2001; pp. 3701-3715. Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2523/71739-MS]
4. Chowdhury, S.; Torres-Verdín, C.; Mezzatesta, A. Interpretation of borehole measurements acquired in laminated clastic sequences subject to mud-filtrate invasion. Proceedings of the SPE Annual Technical Conference and Exhibition; Houston, TX, USA, 26–29 September 2004; Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2118/90199-MS]
5. Jiang, Y.; Sun, J.; Gao, J.; Cui, J. Mud invasion and array laterolog responses in horizontal well based on numerical simulation. J. Geophys. Eng.; 2017; 14, pp. 15-25. [DOI: https://dx.doi.org/10.1088/1742-2132/14/1/15]
6. Dewan, J.T.; Chenevert, M.E. A model for filtration of water-base mud during drilling: Determination of mudcake parameters. Petrophysics; 2001; 42, pp. 237-250.
7. Amanullah, M. A novel laboratory method for assessing the erosional characteristics of mudcakes. Proceedings of the SPE/IADC Drilling Conference; Abu Dhabi, United Arab Emirates, 20–22 October 2003; pp. 411-420. Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2523/85340-MS]
8. Calçada, L.; Scheid, C.; de Araújo, C.; Waldmann, A.; Martins, A. Analysis of dynamic and static filtration and determination of mud cake parameters. Braz. J. Pet. Gas; 2011; 5, pp. 159-170. [DOI: https://dx.doi.org/10.5419/bjpg2011-0016]
9. Khalid, M. Investigating the Effect of Fluid and Formation Parameters on Mud Cake Thickness, Filtration Velocity and Invasion Depth. ResearchGate 2016; Department of Geoscience & Petroleum Engineering Universiti Teknologi PETRONAS, Malaysia. Available online: https://www.researchgate.net/publication/311304765_Investigating_the_Effect_of_Fluid_and_Formation_Parameters_on_Mud_Cake_Thickness_Filtration_Velocity_and_Invasion_Depth?channel=doi&linkId=5841210908ae8e63e6213cf2&showFulltext=true (accessed on 1 January 2023).
10. Adewole, J.K.; Kazeem, T.S.; Oyehan, T.A. Molecular transport of solutes and ions through filter-cake: Implications for wellbore instability and formation damage. Proceedings of the SPE Kingdom of Saudi Arabia Annual Technical Symposium and Exhibition 2018; Dammam, Saudi Arabia, 23–26 April 2018; Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2118/192408-MS]
11. Moazzeni, A.R.; Nabaei, M.; Ghadami, S. Nonproductive time reduction through a deep-rig time analysis: Case study. Proceedings of the SPE Digital Energy Conference and Exhibition; The Woodlands, TX, USA, 19–21 April 2011; SPE–142799.SPE [DOI: https://dx.doi.org/10.2118/142799-MS]
12. Emhana, S. Analysis of non-productive time (NPT) in drilling operations—A case study of the Ghadames Basin. Proceedings of the Second Scientific Conference of Oil and Gas; Ajdabiya, Libya, 5–7 February 2018; pp. 5-7.
13. Donaldson, E.C.; Chernoglazov, V. Characterization of drilling mud fluid invasion. J. Pet. Sci. Eng.; 1987; 1, pp. 3-13. [DOI: https://dx.doi.org/10.1016/0920-4105(87)90010-6]
14. Jiao, D.; Sharma, M. Mechanism of cake buildup in crossflow filtration of colloidal suspensions. J. Colloid Interface Sci.; 1993; 169, pp. 264-270. [DOI: https://dx.doi.org/10.1006/jcis.1994.1060]
15. Corapcioglu, M.Y.; Abboud, N.; Haridas, A. Governing equations for particle transport in porous media. Advances in Transport Phenomena in Porous Media; Springer: Berlin/Heidelberg, Germany, 1987; pp. 269-342. [DOI: https://dx.doi.org/10.1007/978-94-009-3625-6_8]
16. Civan, F.; Engler, T. Drilling mud filtrate invasion—Improved model and solution. J. Pet. Sci. Eng.; 1994; 11, pp. 183-193. [DOI: https://dx.doi.org/10.1016/0920-4105(94)90039-6]
17. Liu, X.; Civan, F. Formation damage and filter cake buildup in laboratory core tests: Modeling and model-assisted analysis. SPE Form. Eval.; 1996; 11, pp. 26-30. [DOI: https://dx.doi.org/10.2118/25215-PA]
18. Parn-Anurak, S.; Engler, T.W. Modeling of fluid filtration and near-wellbore damage along a horizontal well. J. Pet. Sci. Eng.; 2005; 46, pp. 149-160. [DOI: https://dx.doi.org/10.1016/j.petrol.2004.12.003]
19. Shad, S.; Salmanpour, S.; Zamani, H.; Zivar, D. Dynamic analysis of mud loss during overbalanced drilling operation: An experimental study. J. Pet. Sci. Eng.; 2021; 196, 107984. [DOI: https://dx.doi.org/10.1016/j.petrol.2020.107984]
20. Sanfillippo, F.; Brignoli, M.; Santarelli, F.; Bezzola, C. Characterization of conductive fractures while drilling. Proceedings of the SPE European Formation Damage Conference and Exhibition; The Hague, The Netherlands, 2–3 June 1997; SPE–38177.Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2118/38177-MS]
21. Liétard, O.; Unwin, T.; Guillo, D.J.; Hodder, M.H. Fracture width logging while drilling and drilling mud/loss-circulation-material selection guidelines in naturally fractured reservoirs. SPE Drill. Complet.; 1999; 14, pp. 168-177. [DOI: https://dx.doi.org/10.2118/57713-PA]
22. Lavrov, A.; Tronvoll, J. Modeling mud loss in fractured formations. Proceedings of the Abu Dhabi International Petroleum Exhibition and Conference; Abu Dhabi, United Arab Emirates, 10–13 October 2004; SPE–88700.Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2118/88700-MS]
23. Lavrov, A.; Tronvoll, J. Mechanics of borehole ballooning in naturally-fractured formations. Proceedings of the SPE Middle East Oil and Gas Show and Conference; Kingdom of Bahrain, 12–15 March 2005; SPE–93747.Society of Petroleum Engineers [DOI: https://dx.doi.org/10.2118/93747-MS]
24. Xia, Y.; Jin, Y.; Chen, M.; Chen, K.P.; Lin, B.T.; Hou, B. Hydrodynamic modeling of mud loss controlled by the coupling of discrete fracture and matrix. J. Pet. Sci. Eng.; 2015; 129, pp. 254-267. [DOI: https://dx.doi.org/10.1016/j.petrol.2014.07.026]
25. Xia, Y.; Jin, Y.; Chen, M. Comprehensive methodology for detecting fracture aperture in naturally fractured formations using mud loss data. J. Pet. Sci. Eng.; 2015; 135, pp. 515-530. [DOI: https://dx.doi.org/10.1016/j.petrol.2015.10.017]
26. Vargas, D.; Calderon, Z.; Mateus, D. Assessment of the drilling fluid losses phenomenon in naturally-fractured reservoir by numerical modeling that couples fracture plane and permeable zone. Proceedings of the ISRM 8th International Symposium Geomechanics; Bucaramanga, Colombia, 6–10 May 2019; pp. 978-981.
27. Majidi, R.; Miska, S.Z.; Ahmed, R.; Yu, M.; Thompson, L.G. Radial flow of yield-power-law fluids: Numerical analysis, experimental study, and the application for drilling fluid losses in fractured formations. J. Pet. Sci. Eng.; 2010; 70, pp. 334-343. [DOI: https://dx.doi.org/10.1016/j.petrol.2009.12.005]
28. Sun, Y.; Huang, H. Effect of rheology on drilling mud loss in a natural fracture. Proceedings of the 49th US Rock Mechanics/Geomechanics Symposium; San Francisco, CA, USA, 28 June–1 July 2015; Volume 4, pp. 2678-2684.
29. Shahri, P.; Zeyghami, M.; Majidi, R. Investigation of fracture ballooning and breathing in naturally fractured reservoirs: Effect of fracture deformation law. Proceedings of the Nigeria Annual International Conference and Exhibition; Abuja, Nigeria, 30 July–3 August 2011; OnePetro [DOI: https://dx.doi.org/10.2118/150817-MS]
30. Shahri, P.; Mehrabi, M. A new approach in modeling of fracture ballooning in naturally fractured reservoirs. Proceedings of the SPE Kuwait International Petroleum Conference and Exhibition; Kuwait City, Kuwait, 10–12 December 2012; OnePetro [DOI: https://dx.doi.org/10.2118/163382-MS]
31. Lavrov, A. Flow of truncated power-law fluid between parallel walls for hydraulic fracturing applications. J. Non-Newton. Fluid Mech.; 2015; 223, pp. 141-146. [DOI: https://dx.doi.org/10.1016/j.jnnfm.2015.06.005]
32. Bychina, M.; Thomas, G.M.; Khandelwal, R.; Samuel, R. A robust model to estimate the mud loss into naturally fractured formations. Proceedings of the SPE Annual Technical Conference and Exhibition; San Antonio, TX, USA, 9–11 October 2017; D011S005R008.SPE [DOI: https://dx.doi.org/10.2118/187219-MS]
33. Bao, K.; Lavrov, A.; Nilsen, H. Numerical modelling of non-Newtonian fluid flow in fractures and porous media. Proceedings of the ECMOR XV-15th European Conference on the Mathematics of Oil Recovery; Amsterdam, The Netherlands, 29 August–1 September 2016; cp–494.European Association of Geoscientists & Engineers [DOI: https://dx.doi.org/10.3997/2214-4609.201601773]
34. Li, Q.; Li, Q.; Cao, H.; Wu, J.; Wang, F.; Wang, Y. The crack propagation behaviour of CO2 fracturing fluid in unconventional low permeability reservoirs: Factor analysis and mechanism revelation. Processes; 2025; 13, 159. [DOI: https://dx.doi.org/10.3390/pr13010159]
35. Li, Q.; Li, Q.; Wu, J.; Li, X.; Li, H.; Cheng, Y. Wellhead stability during the development process of hydrate reservoir in the Northern South China Sea: Evolution and mechanism. Processes; 2024; 13, 40. [DOI: https://dx.doi.org/10.3390/pr13010040]
36. Windarto, G.; Sukarno, P.; Soewono, E. Modeling of mud filtrate invasion and damage zone formation. J. Pet. Sci. Eng.; 2011; 77, pp. 359-364. [DOI: https://dx.doi.org/10.1016/j.petrol.2011.04.011]
37. Windarto, G.; Sukarno, P.; Soewono, E. Modelling of formation damage due to mud filtrate invasion in a radial flow system. J. Pet. Sci. Eng.; 2012; 100, pp. 99-105. [DOI: https://dx.doi.org/10.1016/j.petrol.2012.11.003]
38. Ring, T. Bed filtration of colloidal particles. Prepr. Soc. Min. Eng. AIME; 1983; 276, pp. 1838-1842.
39. Fetter, C.; Boving, T.; Kreamer, D. Contaminant hydrogeology. Contaminant hydrogeology; Waveland Press: Long Grove, IL, USA, 2017; [DOI: https://dx.doi.org/10.1016/0016-7037(93)90162-p]
40. Coutelieris, F.; Delgado, J. Advanced Structured Materials: Transport Processes in Porous Media; Springer: Berlin/Heidelberg, Germany, 2012; Volume 20.
41. Proett, M.A.; Chin, W.; Wu, J.; Manohar, M.; Belanger, D. Sample quality prediction with integrated oil and water-based mud invasion modeling. Proceedings of the SPE—Asia Pacific Oil and Gas Conference; Melbourne, Australia, 8–10 October 2002; pp. 1081-1094. [DOI: https://dx.doi.org/10.2118/77964-ms]
42. Fisher, K.; Wakeman, R.; Chiu, T.; Meuric, O. Numerical modelling of cake formation and fluid loss from non-Newtonian muds during drilling using eccentric/concentric drill strings with/without rotation. Chem. Eng. Res. Des.; 2000; 78, pp. 707-714. [DOI: https://dx.doi.org/10.1205/026387600527888]
43. Byron, B. Fenómenos de Transporte; 2nd ed. Wikipedia: San Francisco, CA, USA, 2006; Volume 2.
44. van Opedal, N.; Cerasi, P.; Ytrehus, J.D. Dynamic fluid erosion on filter cakes. Proceedings of the SPE European Formation Damage Conference and Exhibition; Noordwijk, The Netherlands, 5–7 June 2013; SPE–165107.SPE [DOI: https://dx.doi.org/10.2118/165107-MS]
45. Hashemzadeh, S.M.; Hajidavalloo, E. Numerical investigation of filter cake formation during concentric/eccentric drilling. J. Pet. Sci. Eng.; 2016; 145, pp. 161-167. [DOI: https://dx.doi.org/10.1016/j.petrol.2016.03.024]
46. Ahammad, F.; Mahmud, S.; Islam, S. Computational fluid dynamics study of yield power law drilling fluid flow through smooth-walled fractures. J. Pet. Explor. Prod. Technol.; 2019; 9, pp. 2717-2727. [DOI: https://dx.doi.org/10.1007/s13202-019-0646-5]
47. Vargas, D.; Uribe, L.; Calderón, Z.; Bernal, R. Estimación de pérdidas de lodo a través de fracturas mediante modelamiento matemático y su validación mediante un modelo físico. Boletín De Geol.; 2024; 46, pp. 103-114. [DOI: https://dx.doi.org/10.18273/revbol.v46n1-2024006]
48. Vargas, D. Modelamiento Numérico de Pérdidas de Fluido de Perforación En Zonas Fracturadas; Universidad Industrial de Santander: Bucaramanga, Colombia, 2017; Available online: https://noesis.uis.edu.co/handle/20.500.14071/37566 (accessed on 1 January 2023).
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
In well drilling operations, interactions between drilling fluid water-based and the well-bore present significant challenges, often escalating project costs and timelines. Particularly, fractures (both induced and natural) and permeable zones at the wellbore can result in substantial mud loss or increased filtration. Addressing these challenges, our research introduces a novel coupled numerical model designed to precisely calculate fluid losses in fractured and permeable zones. For the permeable zone, fundamental variables such as filtration velocity, filtrate concentration variations, permeability reduction, and fluid cake growth are calculated, all based on the law of continuity and convection-dispersion theory. For the fracture zone, the fluid velocity profile is determined using the momentum balance equation and both Newtonian and non-Newtonian rheology. The model was validated against laboratory data and physical models, and adapted for field applications. Our findings emphasize that factors like mud particle size, shear stress, and pressure differential are pivotal. Effectively managing these factors can significantly reduce fluid loss and mitigate formation damage caused by fluid invasion. Furthermore, the understanding gathered from studying mud behavior in both permeable and fractured zones equips drilling personnel with valuable information related to the optimal rheological properties according to field conditions. This knowledge is crucial for optimizing mud formulations and strategies, ultimately aiding in the reduction of non-productive time (NPT) associated with wellbore stability issues.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details



1 Wellbore Stability Research Group, Petroleum Engineering School, Industrial University of Santander, Bucaramanga 680002, Colombia; [email protected]
2 ECOPETROL, Bogotá 681011, Colombia; [email protected]
3 Polymer Research Group, Chemistry Engineering School, Industrial University of Santander, Bucaramanga 680002, Colombia; [email protected]