A fully implicit and consistent finite element framework for modeling reservoir compaction with large deformation and nonlinear flow model. Part II: verification and numerical example

A fully implicit, fully coupled, and fully consistent finite element framework has been formulated in part I of this work for modeling reservoir compaction through linearizing coupled solid and flow field equations and constructing a local material integrator. In part II of this work, we focus on verification and performance analysis of our numerical formulation and computer implementation using several numerical examples. First, we design a cube problem in triaxial compression to verify our numerical formulation and computer code implementation especially for rock formation in compaction using cap plasticity models. The finite element prediction on stresses is compared with the analytical solution. The second problem we select is a strip footing problem popular in the geotechnical area where the evolution of soil consolidation degrees following the diffusion of pore pressure is the main interest. In this example, we demonstrate a good performance of the proposed numerical formulation on solving different shear and compaction-dominated deformation behaviors by varying the footing length. Importantly, an extremely sharp cap model based on real experimental data for Leda clays, a challenging cap model, is successfully applied in this footing problem. Our focus in this work is to model field reservoirs undergoing serious compaction. A reservoir with complex payzone geometries, multiple horizontal wells, and cap plasticity models with sharp cap surfaces has been successfully solved using our fully implicit formulation. The last example is to model a horizontal wellbore damage problem. Finally, the sensitivity of predicted subsidence to nonlinear flow model, cap hardening parameters, and Lode angles have been systemically investigated and documented in detail, which can provide a constructive guidance on how to successfully model field reservoir compaction problems with cap plasticity models.     

Discrete fracture model for coupled flow and geomechanics

We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well.     

Parameters controlling pressure and fracture behaviors in field injectivity tests: A numerical investigation using coupled flow and geomechanics model

Field injectivity tests are widely used in the oil and gas industry to obtain key formation characteristics. The prevailing approaches for injectivity test interpretation rely on traditional analytical models. A number of parameters may affect the test results and lead to interpretation difficulties. Understanding their impacts on pressure response and fracture geometry of the test is essential for accurate test interpretation. In this work, a coupled flow and geomechanics model is developed for numerical simulation of field injectivity tests. The coupled model combines a cohesive zone model for simulating fluid-driven fracture and a poro-elastic/plastic model for simulating formation behavior. The model can capture fracture propagation, fluid flow within the fracture and formation, deformation of the formation, and evolution of pore pressure and stress around the wellbore and fracture during the tests. Numerical simulations are carried out to investigate the impacts of a multitude of parameters on test behaviors. The parameters include rock permeability, the leak-off coefficient of the fracture, rock stiffness, rock toughness, rock strength, plasticity deformation, and injection rate. The sensitivity of pressure response and fracture geometry on each parameter is reported and discussed. The coupled flow and geomechanics model provides additional advantages in the understanding of the fundamental mechanisms of field injectivity tests.     

A finite element solution of a fully coupled implicit formulation for reservoir simulation

An iterative method is presented for solving a fully coupled and implicit formulation of fluid flow in a porous medium. The mathematical model describes a set of fully coupled three-phase flow of compressible and immiscible fluids in a saturated oil reservoir. The finite element method is applied to obtain the simultaneous solution (SS) for the resulting highly non-linear partial differential equations where fluid pressures are the primary unknowns. The final discretized equations are solved iteratively by using a fully implicit numerical scheme. Several examples, illustrating the use of the present model, are described. The increased stability achieved with this scheme has permitted the use of larger time steps with smaller material balance errors.     

A fully coupled thermo-hydro-mechanical nonlinear model for a frozen medium

This paper describes a nonlinear elasto-plastic simulation of freezing and thawing of rock. A mathematical formulation is described in which deformation, fluid flow and heat flow are fully coupled. A non-linear elasto-plastic constitutive relationship is presented and a two dimensional (plane stress) numerical modeling is performed based on the finite element method applied to thermo-poro-elastoplasticity. It is assumed that the Mohr-Coulomb's failure criterion is valid for yield locus and plastic potential. The numerical scheme employed in the code accommodates phase change of pore-water from liquid to solid (ice). The primary aim of this paper is to compare the temperature transfer and deformation prediction obtained from the numerical code with those obtained from the freezing and thawing experiments. It is found from the numerical simulation that a relatively good prediction can be made of temperature transfer and deformation behavior. The numerical code has also been applied to a hypothetical cavern problem to demonstrate its applicability.     

An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation

A numerical modeling framework is described that is able to calculate the coupled processes of fluid flow, geomechanics, and rock failure for application to general engineering problems related to reservoir stimulation, including hydraulic fracturing and shear stimulation. The numerical formulation employs the use of an embedded fracture modeling approach, which provides several advantages over more traditional methods in terms of computational complexity and efficiency. Specifically, the embedded fracture modeling strategy avoids the usual requirement that the discretization of the fracture domain conforms to the discretization of the rock volume surrounding the fractures. As fluid is exchanged between the two domains, conservation of mass is guaranteed through a coupling term that appears as a simple source term in the governing mass balance equations. In this manner, as new tensile fractures nucleate and propagate subject to mechanical effects, numerical complexities associated with the introduction of new fracture control volumes are largely negated. In addition, the ability to discretize the fractures and surrounding rock volume independently provides the freedom to choose an acceptable level of discretization for each domain separately. Three numerical examples were performed to demonstrate the utility of the embedded fracture model for application to problems involving fluid flow, mechanical deformation, and rock failure. The results of the numerical examples confirm that the embedded fracture model was able to capture accurately the complex and nonlinear evolution of reservoir permeability as new fractures propagate through the reservoir and as fractures fail in shear.     

Finite deformation analysis of geomaterials

The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd.     

Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

用流固耦合方法研究油藏压裂后应力应变和孔渗特性变化

油藏压裂后将引起地应力场发生变化，使岩石变形，导致孔隙度和渗透率变化，进而影响产量，为研究这一问题，作者建立了油藏压裂后流－固耦合渗流模型，考虑了以下因素：油藏岩石变形，地应力，孔隙度和渗透率变化，人工裂缝，流体渗流与岩石应变耦合，储藏渗流与裂缝渗流耦合，非达西效应等。较详细地给出了耦合方程及推导过程，控制方程包括的未知变量有压力，饱和度及位移，11个变量，和11个方程，用有限差分方法将流体渗流和岩石应变方程离散成主对角占优的七对角矩阵，可在修改已有三维二相渗流和三维固体力学程序的基础上，采用隐式迭代方法求解，示例分析表明，用此模型可以研究储层应力变变，孔隙度和渗透率随时间和空间变化规律，为开发方案制定，整体压裂设计，压后生产管理等方面提供定量分析技术。     

A fully coupled numerical modeling to investigate the role of rock thermo-mechanical properties on reservoir uplifting in steam assisted gravity drainage

One of the crucial consequences of steam assisted gravity drainage (SAGD) process is abnormal reservoir uplifting under thermal steam injection, which can significantly influence the reservoir rock deformation, specifically thin bed reservoirs and causes intensive failures and fractures into the cap rock formations. A thorough understanding of the influences of rock thermo-mechanical properties on reservoir uplifting plays an important role in preventing those aforementioned failures within design and optimization process in SAGD. In addition, coupling of reservoir porous medium and flowing of specific fluid with temperature as an additional degree of freedom with initial pore pressure and in-situ stress condition, are also very challenging parts of geomechanical coupled simulation which would be clearly explained. Thus, a fully coupled thermo-poro-elastic geomechanical model with finite element codes was performed in ABAQUS to investigate the role of rock thermo-mechanical parameters on reservoir vertical uplift during steam injection. It is clearly observed that, any increase in rock thermo-mechanical properties specifically rock’s thermal properties such as specific heat, thermal expansion, and formation’s thermal conductivity, have significant influences on reservoir uplift. So by coupling the temperature as an additional degree of freedom with the coupled pore-fluid stress and diffusion finite element model of SAGD process, the more realistic simulation will be conducted; hence, the errors related to not having heat as an additional degree of freedom will be diminished. In addition, Young’s modulus and specific heat are the rock thermo-mechanical parameters which have the maximum and minimum effects on the reservoir uplift, respectively.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

Cell to node projections: An assessment of error

Reservoir simulators typically use cell‐centered finite volume schemes and do not model directly the coupling of the flow processes with the geomechanics. Coupling of geomechanics with fluid flow can be important in many cases, but introducing fully coupled geomechanical effects in those simulators is not a trivial issue, because the geomechanics is better done by using the Galerkin vertex‐centered finite element methods by which the solid displacements are computed at the vertices of the cells. This creates difficulties in interfacing cell variables with nodal variables. Uncoupled or loosely coupled models are used by many researchers/practitioners by which a reservoir model is coupled to a geomechanical model by staggering in‐time flow and deformation via a sophisticated interface that repeatedly calls first flow and then mechanics. The method therefore requires projection of the reservoir cell variables onto the nodes of the geomechanics Galerkin finite element mesh. In this note, we attempt to quantify the errors associated with cell to node projection operations. For that purpose, we use a simple model of the pressure equation for a heterogeneous medium in one dimension. We are able to derive the exact analytical solution for this problem for both nodal and cell pressures. This allows us to compute the errors due to projection analytically, function of meshing refinement and permeability field variations. We compute upper and lower bounds for the errors, and analyze their magnitude for a variety of cases. We conclude that, in general, cell to node projection operations lead to substantial errors. Copyright © 2010 John Wiley & Sons, Ltd.     

A locally conservative finite element framework for the simulation of coupled flow and reservoir geomechanics

In this paper, we present a computational framework for the simulation of coupled flow and reservoir geomechanics. The physical model is restricted to Biot’s theory of single-phase flow and linear poroelasticity, but is sufficiently general to be extended to multiphase flow problems and inelastic behavior. The distinctive technical aspects of our approach are: (1) the space discretization of the equations. The unknown variables are the pressure, the fluid velocity, and the rock displacements. We recognize that these variables are of very different nature, and need to be discretized differently. We propose a mixed finite element space discretization, which is stable, convergent, locally mass conservative, and employs a single computational grid. To ensure stability and robustness, we perform an implicit time integration of the fluid flow equations. (2) The strategies for the solution of the coupled system. We compare different solution strategies, including the fully coupled approach, the usual (conditionally stable) iteratively coupled approach, and a less common unconditionally stable sequential scheme. We show that the latter scheme corresponds to a modified block Jacobi method, which also enjoys improved convergence properties. This computational model has been implemented in an object-oriented reservoir simulator, whose modular design allows for further extensions and enhancements. We show several representative numerical simulations that illustrate the effectiveness of the approach.     

A new locally conservative numerical method for two-phase flow in heterogeneous poroelastic media

We construct a new class of locally conservative numerical methods for two-phase immiscible flow in heterogeneous poroelastic media. Within the framework of the so-called iteratively coupled methods and fixed-stress split algorithm we develop mixed finite element methods for the flow and geomechanics subsystems which furnish locally conservative Darcy velocity and transient porosity input fields for the transport problem for the water saturation. Such hyperbolic equation is decomposed within an operator splitting technique based on a predictor–corrector scheme with the predictor step discretized by a higher-order non-oscillatory finite volume central scheme. The proposed scheme adopts an inhomogeneous dual mesh with variable cell size ruled by the local wave speed of propagation to compute numerical fluxes at cell edges. In the limit of small time steps the central scheme gives rise to a semidiscrete formulation for the water saturation capable of incorporating heterogeneous porosity fields and generalized flux functions including the water transport due to the solid phase velocity. Numerical simulations of a water-flooding problem in secondary oil recovery are presented for different realizations of the input random fields (permeability, Young modulus and initial porosity). Comparison between the accuracies of the proposed approach and the traditional one-way coupled hydro-geomechanical formulation are presented. The effects of the cross-correlation between the input random fields and compaction drive mechanism upon finger growth and breakthrough curves are also analyzed. A notable feature of the formulation proposed herein is the accurate prediction of the influence of geomechanical effects upon the unstable movement of the water front, whose evolution is dictated by rock heterogeneity and unfavorable viscosity ratio, without deteriorating the local conservative character of the numerical schemes.     

Response of saturated and nearly saturated porous media: Different formulations and their applicability

The response of saturated porous medium is of significant interest in many fields ranging from geomechanics to biomechanics. Biot was the first to formulate the basic equations governing the process of coupled flow and deformation in porous media. Depending on the nature of loading vis‐à‐vis the characteristics of the media, different formulations (fully dynamic, partly dynamic, quasi‐static) are possible. In this study, analytical solutions are developed for the response of saturated and nearly saturated porous media under plane strain condition. The solutions for different formulations are developed in terms of non‐dimensional parameters. The response is studied for various conditions and the regions of validity for various formulations are identified in a parametric space. An assessment of the needed formulation for few important problems is also presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Two‐dimensional finite element analysis of gravitational and lateral‐driven deformation in sedimentary basins

The paper deals with the modeling of some aspects, such as the formulation of constitutive equations for sediment material or finite element approach for basin analysis, related to mechanical compaction in sedimentary basins. In addition to compaction due to gravity forces and pore‐pressure dissipation, particular emphasis is given to the study of deformation induced by tectonic sequences. The numerical model relies upon the implementation of a comprehensive constitutive model for the sediment material formulated within the framework of finite poroplasticity. The theoretical model accounts for both hydromechanical and elasticity–plasticity coupling due to the effects of irreversible large strains. From the numerical viewpoint, a finite element procedure specifically devised for dealing with sedimentary basins as open systems allows to simulate within a two‐dimensional setting the process of sediment accretion or erosion. Several basin simulations are presented. The main objective is to analyze the behavior of a sedimentary basin during the different phases of its life cycle: accretion phase, pore‐pressure dissipation phase and compressive/extensional tectonic motions. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling mechanical behaviour of limestone under reservoir conditions

High porosity and low permeability limestone has presented pore collapse. As fluid is withdrawn from these reservoirs, the effective stresses acting on the rock increase. If the strength of the rock is overcome, pore collapse may occur, leading to irreversible compaction of porous media with permeability and porosity reduction. It impacts on fluid withdrawal. Most of reservoirs have been discovered in weak formations, which are susceptible to this phenomenon. This work presents a study on the mechanical behaviour of a porous limestone from a reservoir located in Campos Basin, offshore Brazil. An experimental program was undergone in order to define its elastic plastic behaviour. The tests reproduced the loading path conditions expected in a reservoir under production. Parameters of the cap model were fitted to these tests and numerical simulations were run. The numerical simulations presented a good agreement with the experimental tests. Copyright © 2006 John Wiley & Sons, Ltd.     

Bifurcation analysis of deep boreholes: I. Surface instabilities

In this paper a bifurcation analysis of boreholes in deep rock formations under uniform stress at infinity is presented. Rock is described by the constitutive equations of a deformation theory of plasticity for rigidplastic, incompressible and cohesive-frictional material. The corresponding bifurcation problem is solved numerically by the finite element method. The evidence gained from the numerical solution is utilized to establish a simplified borehole stability analysis that combines Biot's hodograph method with a surface instability analysis.     

Implicit integration algorithm for Hoek-Brown elastic-plastic model

A realistic strength criterion often used to describe the yielding behaviour of a jointed rock mass at a continuum level is the well-known Hoek and Brown criterion. This paper is concerned with a 3-D stress generalization of the Hoek-Brown failure criterion by means of an elliptical functional which leads to a smooth deviatoric trace in the stress space. For its incorporation into a finite element analysis involving plasticity calculations, the formulation of an implicit stress integration algorithm is presented. The key computational methodology alludes to the notion of consistent tangent modulus and implicit return mapping schemes (radial and closest point return) for stress integration in a finite element analysis. Within the context of non-linear elastoplastic analysis, it is found that formulation of such consistent modulus and success into achieving numerical efficiency are closely intertwined. Indeed, the procedure results into accurate and rapid convergence of the displacement finite element scheme during the search for equilibrium. This means that considerable savings in computational time can be achieved for large scale problems. Numerical examples which focus on the Hoek-Brown plasticity model are presented in order to fully appreciate the robustness of the algorithm, and hence the viability of such method to practical problems.