Discrete fracture model for simulating waterflooding processes under fracturing conditions

Gallyamov, Emil; Garipov, Timur; Voskov, Denis; Hoek, Paul van den

doi:10.1002/nag.2797

Cited by 22 publications

(18 citation statements)

References 61 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Boundary stresses can significantly impact the onset, type, and magnitude of failure in rocks subjected to hydraulic stimulation because they determine the solution of the equilibrium equation, Equation (35). We focus on the effect of anisotropy in the boundary stresses on permeability growth and poroelastic damage.…”

Section: Case V: Boundary Stress Ratiomentioning

confidence: 99%

“…In cohesive zone models, the crack path is often defined a priori such that the locations of the cohesive elements are known. The discrete fracture model approach is similar in terms of its requirement to know the fracture locations and path a priori . These approaches lead to mesh‐dependent fracture path solutions because the fractures can propagate only along element edges.…”

Section: Introductionmentioning

confidence: 99%

“…The discrete fracture model approach is similar in terms of its requirement to know the fracture locations and path a priori. 35 These approaches lead to mesh-dependent fracture path solutions because the fractures can propagate only along element edges. Simulation of flow and deformation processes in fractured reservoirs using a mesh that explicitly represents multiple fractures and conform to individual fracture geometry is computationally prohibitive.…”

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Coupled poromechanics‐damage mechanics modeling of fracturing during injection in brittle rocks

Bubshait¹,

Jha²

2019

Numerical Meth Engineering

View full text Add to dashboard Cite

Summary Subsurface injection of fluids in a stress‐sensitive naturally fractured rock faces the problem of near‐wellbore fracture evolution and associated changes in rock properties. Numerical modeling of the changes in permeability and poroelastic properties of the near‐wellbore region is challenging due to the coupling between fracture dynamics and poromechanics across multiple length scales of fractures and the host rock. We present a numerical framework to model anisotropic and dynamic evolution in rock properties based on a coupled formulation of fluid flow, rock mechanics, and fracturing, where fracturing‐induced damage is used to update the rock properties. A generalized fixed‐stress method, which accounts for damage‐induced anisotropy in flow and deformation processes, is developed to sequentially solve the equations of flow, mechanics, and fracture evolution. We demonstrate the usefulness of our framework in quantifying the effects of injection rate variability, initial fracture distribution, and in‐situ stress state on the evolution in permeability, elastic stiffness, and the Biot parameters. Our framework does not require the computational mesh to conform to existing or future fractures, allows simultaneous growth of multiple randomly distributed fractures, and can be implemented relatively easily in existing coupled flow‐geomechanics simulators to extend them to model fracturing at the reservoir scale.

show abstract

Section: Case V: Boundary Stress Ratiomentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Coupled poromechanics‐damage mechanics modeling of fracturing during injection in brittle rocks

Bubshait¹,

Jha²

2019

Numerical Meth Engineering

View full text Add to dashboard Cite

show abstract

“…Considering fractures that have a negligible aperture compared to the modelled domain, this representation can be combined with an approach in which fractures are modelled as lower-dimensional structures (Martin et al 2005;Karimi-Fard et al 2003) and discretised with elements of zero thickness (Király 1988;Boon et al 2018;Flemisch et al 2018;Berre et al 2020). For poroelastic media with fractures, this allows for the application of standard finite element (Salimzadeh et al 2017), finite volume (Ucar et al 2018;Berge et al 2020) and combined finite element/finite volume schemes (Garipov et al 2016;Settgast et al 2017;Garipov and Hui 2019), more recently also including fracture contact mechanics (Garipov et al 2016;Gallyamov et al 2018;Franceschini et al 2020) and tensile fracture propagation (Settgast et al 2017;Salimzadeh et al 2017). Thermal effects on fracture deformation and propagation in poroelastic media are less studied, although some recent studies model deformation of existing fractures in thermo-poroelastic media (Salimzadeh et al 2016(Salimzadeh et al , 2018aStefansson et al 2021;Garipov and Hui 2019).…”

Section: Introductionmentioning

confidence: 99%

Numerical Modelling of Convection-Driven Cooling, Deformation and Fracturing of Thermo-Poroelastic Media

et al. 2021

View full text Add to dashboard Cite

Convection-driven cooling in porous media influences thermo-poro-mechanical stresses, thereby causing deformation. These processes are strongly influenced by the presence of fractures, which dominate flow and heat transfer. At the same time, the fractures deform and propagate in response to changes in the stress state. Mathematically, the model governing the physics is tightly coupled and must account for the strong discontinuities introduced by the fractures. Over the last decade, and motivated by a number of porous media applications, research into such coupled models has advanced modelling of processes in porous media substantially. Building on this effort, this work presents a novel model that couples fracture flow and heat transfer and deformation and propagation of fractures with flow, heat transfer and thermo-poroelasticity in the matrix. The model is based on explicit representation of fractures in the porous medium and discretised using multi-point finite volume methods. Frictional contact and non-penetration conditions for the fractures are handled through active set methods, while a propagation criterion based on stress intensity factors governs fracture extension. Considering both forced and natural convection processes, numerical results show the intricate nature of thermo-poromechanical fracture deformation and propagation.

show abstract

“…To overcome the drawbacks of DFM, several upscaling approaches were proposed (Karimi-Fard et al, 2006;Gong, 2007;Karimi-Fard et al, 2012;Karimi-Fard et al, 2016). In addition, the DFM approach was extended to account for geomechanics effects (Garipov et al, 2016) and fracture propagation (Gallyamov et al, 2018).…”

Section: Introductionmentioning

confidence: 99%

Multi-Level Discrete Fracture Model For Carbonate Reservoirs

Voskov

2018

Proceedings

View full text Add to dashboard Cite

The main challenge for predictive simulation of carbonate reservoirs is associated with large uncertainties in the geological characterization of the reservoir with multiple features including fractures and cavities. This type of reservoirs requires robust and efficient forward-simulation capabilities to apply data assimilation or optimization technique under uncertainties. The interaction between reservoir matrix and various features introduces a complex multi-scale flow response driven by global boundary conditions. The Discrete Fracture Models (DFM), which represent fractures explicitly, is capable to accurately depict all important features of flow behavior. However, these models are constrained by many degrees of freedom and corresponding computational efficiency when the fracture network becomes complicated. The Embedded DFM, which represents the interaction between matrix and fractures analytically, is an efficient approximation. However, it cannot accurately reproduce the effect of local flow conditions, especially when the secondary fractures are present. In this study, we applied a numerical upscaling of DFM models to a triple continuum model where large features are represented explicitly using the numerical EDFM and small features are upscaled as a third continuum. In this approach, we discretize the original geo-model with unstructured grid based on DFM and associate the mesh geometry with large features in the model. Using the global solution, we generate local boundary conditions for the model capturing the response of primary features to the flow. Applying local boundary conditions, we resolve all secondary features using a fine scale solution and update the local boundary conditions. This procedure is applied iteratively using the local-global-upscaling formalism. To demonstrate the accuracy of the Multi-Level Discrete Fracture Model, several realistic cases have been tested. By comparing with fine scale DFM solution and the traditional EDFM technique, we demonstrate that the proposed model is accurate enough to capture the flow behavior in complex fractured systems with advanced computational efficiency.

show abstract

Discrete fracture model for simulating waterflooding processes under fracturing conditions

Cited by 22 publications

References 61 publications

Coupled poromechanics‐damage mechanics modeling of fracturing during injection in brittle rocks

Coupled poromechanics‐damage mechanics modeling of fracturing during injection in brittle rocks

Numerical Modelling of Convection-Driven Cooling, Deformation and Fracturing of Thermo-Poroelastic Media

Multi-Level Discrete Fracture Model For Carbonate Reservoirs

Contact Info

Product

Resources

About