Iteratively coupled solution strategies for a four‐field mixed finite element method for poroelasticity

In this paper, we consider numerical algorithms for modeling of the time‐dependent coupling between the fluid flow and deformation in elastic porous media. Here, we employ a four‐field formulation which uses the total stress, displacement, flux, and pressure as its primary variables and satisfies Darcy's law and linear elasticity in mixed weak form. We present four different iteratively coupled methods, known as drained, undrained, fixed‐strain, and fixed‐stress splits, in which the diffusion operator is separated from the elasticity operator and the two subproblems are solved in a staggered way while ensuring convergence of the solution at each time step. A‐priori convergence results for each iterative coupling which differs from those found when using a traditional two‐field or three‐field formulation are presented. We also present some numerical results to support the convergence estimates and to show the accuracy and efficiency of the algorithms. Copyright © 2016 John Wiley & Sons, Ltd.     

Convergence analysis of fixed stress split iterative scheme for anisotropic poroelasticity with tensor Biot parameter

We perform a convergence analysis of the fixed stress split iterative scheme for the Biot system modeling coupled flow and deformation in anisotropic poroelastic media with tensor Biot parameter. The fixed stress split iterative scheme solves the flow subproblem with all components of the stress tensor frozen using a multipoint flux mixed finite element method, followed by the poromechanics subproblem using a conforming Galerkin method in every coupling iteration at each time step. The coupling iterations are repeated until convergence and Backward Euler is employed for time marching. The convergence analysis is based on studying the equations satisfied by the difference of iterates to show that the fixed stress split iterative scheme for anisotropic poroelasticity with Biot tensor is contractive. We also demonstrate that the scheme is numerically convergent using the classical Mandel’s problem solution for transverse isotropy.     

Convergence of iterative coupling for coupled flow and geomechanics

In this paper, we study solving iteratively the coupling of flow and mechanics. We demonstrate the stability and convergence of two widely used schemes: the undrained split method and the fixed stress split method. To our knowledge, this is the first time that such results have been rigorously obtained and published in the scientific literature. In addition, we propose a new stress split method, with faster convergence rate than known schemes. These results are specially important today due to the interest in hydraulic fracturing (Dean and Schmidt SPE J. 14:707–714, 2009; Ji et al. SPE J. 14:423–430, 2009; Samier and De Gennaro 2007; Settari and Maurits SPE J. 3:219–226, 1998), in oil and gas shale reservoirs.     

A three-phase XFEM model for hydraulic fracturing with cohesive crack propagation

A three-phase hydro-mechanical model for hydraulic fracturing is proposed. Three phases include: porous solid, fracturing fluid and host fluid. Discontinuity is handled using extended finite element method (XFEM) while cohesive crack model is used as fracturing criterion. Flow through fracture is defined as one-dimensional laminar flow, and flow through porous medium (host reservoir) is defined as two-dimensional Darcy flow. Coupling between two fluids in each space, fracture and pore, is captured through capillary pressure–saturation relationship, while the identical fluids in fracture and pore are coupled through a so-called leak-off mass transfer term. Coupling between fluids and deformation is captured through compatibility of volumetric strain of fluids within fracture and pore, and volumetric strain of the matrix. Spatial and temporal discretisation is achieved using the standard Galerkin method and the finite difference technique, respectively. The model is verified against analytical solutions available from literature. The leaking of fracturing fluid into the medium and suction of porous fluid into the fracture around the tip, are investigated. Sensitivity analyses are carried out for cases with slow and fast injection rates. It is shown that the results by single-phase flow may underestimate the leak-off.     

Robust iterative schemes for non-linear poromechanics

We consider a non-linear extension of Biot’s model for poromechanics, wherein both the fluid flow and mechanical deformation are allowed to be non-linear. Specifically, we study the case when the volumetric stress and the fluid density are non-linear functions satisfying certain assumptions. We perform an implicit discretization in time (backward Euler) and propose two iterative schemes for solving the non-linear problems appearing within each time step: a splitting algorithm extending the undrained split and fixed stress methods to non-linear problems, and a monolithic L-scheme. The convergence of both schemes are shown rigorously. Illustrative numerical examples are presented to confirm the applicability of the schemes and validate the theoretical results.     

Performance studies of the fixed stress split algorithm for immiscible two-phase flow coupled with linear poromechanics

In this work, we measure the performance of the fixed stress split algorithm for the immiscible water-oil flow coupled with linear poromechanics. The two-phase flow equations are solved on general hexahedral elements using the multipoint flux mixed finite element method whereas the poromechanics equations are discretized using the conforming Galerkin method. We introduce a rigorous calculation of the update in poroelastic properties during the iterative solution of the coupled system equations. The effects of the coupling parameter on the performance of the fixed stress algorithm is demonstrated in two field studies: the Frio oil reservoir and the Cranfield injection site.

     

A priori error estimates for the numerical solution of a coupled geomechanics and reservoir flow model with stress-dependent permeability

In this paper we consider the numerical solution of a coupled geomechanics and a stress-sensitive porous media reservoir flow model. We combine mixed finite elements for Darcy flow and Galerkin finite elements for elasticity. This work focuses on deriving convergence results for the numerical solution of this nonlinear partial differential system. We establish convergence with respect to the L ²-norm for the pressure and for the average fluid velocity and with respect to the H ¹-norm for the deformation. Estimates with respect to the L ²-norm for mean stress, which is of special importance since it is used in the computation of permeability for poro-elasticity, can be derived using the estimates in the H ¹-norm for the deformation. We start by deriving error estimates in a continuous-in-time setting. A cut-off operator is introduced in the numerical scheme in order to derive convergence. The spatial grids for the discrete approximations of the pressure and deformation do not need be the same. Theoretical convergence error estimates in a discrete-in-time setting are also derived in the scope of this investigation. A numerical example supports the convergence results.     

Multiscale simulations of porous media flows in flow-based coordinate system

In this paper, we propose a multiscale technique for the simulation of porous media flows in a flow-based coordinate system. A flow-based coordinate system allows us to simplify the scale interaction and derive the upscaled equations for purely hyperbolic transport equations. We discuss the applications of the method to two-phase flows in heterogeneous porous media. For two-phase flow simulations, the use of a flow-based coordinate system requires limited global information, such as the solution of single-phase flow. Numerical results show that one can achieve accurate upscaling results using a flow-based coordinate system.     

Homogenization of a Darcy–Stokes system modeling vuggy porous media

We derive a macroscopic model for single-phase, incompressible, viscous fluid flow in a porous medium with small cavities called vugs. We model the vuggy medium on the microscopic scale using Stokes equations within the vugular inclusions, Darcy's law within the porous rock, and a Beavers–Joseph–Saffman boundary condition on the interface between the two regions. We assume periodicity of the medium and obtain uniform energy estimates independent of the period. Through a two-scale homogenization limit as the period tends to zero, we obtain a macroscopic Darcy's law governing the medium on larger scales. We also develop some needed generalizations of the two-scale convergence theory needed for our bimodal medium, including a two-scale convergence result on the Darcy–Stokes interface. The macroscopic Darcy permeability is computable from the solution of a cell problem. An analytic solution to this problem in a simple geometry suggests that: (1) flow along vug channels is primarily Poiseuille with a small perturbation related to the Beavers–Joseph slip, and (2) flow that alternates from vug to matrix behaves as if the vugs have infinite permeability.     

Computation of the filtration laws through porous media for a non-Newtonian fluid obeying the power law

We consider a stationary flow of an incompressible non-Newtonian flow through a porous medium, induced by an injection velocity when inertial effects are negligible. At the pore scale, the governing equations are based on a nonlinear relation between the stress and the rate of deformation. In such a situation, the limit problem obtained when the pore size tends to zero, is called the homogenized problem that leads to the filtration law. This filtration law is given by a non-linear system coupling a local problem on a typical cell of the porous medium to a global problem at the scale of the whole porous medium. We propose, in this work, a numerical method to solve this homogenized problem and apply this method when the velocity dependent viscosity is given by the power law. Finally, we propose some numerical experiments to illustrate our approach.     

煤层水力压裂应力与裂隙演化的细观规律

水力压裂作为煤层强化增透技术的一种,其应力演化特征及裂隙形态与扩展范围的判断尤为重要。采用离散元数值方法,以导向压裂为背景,建立水力压裂流固耦合模型;通过应力路径、裂纹热点图等手段,探究水力压裂过程中压裂排量、泊松比、天然裂隙密度对应力演化和裂隙演化的影响及其细观规律。结果表明:不同压裂排量下的应力演化方向及最终应力路径曲线形状有着明显的不同,低排量下裂隙附近的应力比值逐渐增大,而在高排量下先增大后减小;煤层泊松比越大,平均压裂半径越低,但对起裂时间及裂隙的扩展形态影响不明显;天然裂隙的发育情况对水力裂隙的扩展起着关键性作用,高裂隙发育煤层水力裂隙扩展的方向性无法预测,应力演化方向会出现反转现象;压裂过程中不同区域的应力演化特征能够反映出裂隙的扩展状态,现场可通过监测压裂区域附近应力变化,判断水力压裂缝网的扩展范围。      

The representer method for data assimilation in single-phase Darcy flow in porous media

We apply the representer method, a data assimilation algorithm, to single-phase Darcy flow in porous media. The measurement array that yields the assimilated data can be expressed as a vector of linear functionals of pressure. The a priori discretization errors in the representer method are analyzed in terms of the convergence properties of the underlying numerical schemes used in each part of the algorithm. We formulate some proof-of-concept numerical experiments that illustrate the error analysis.     

Fully coupled hydro-mechanical numerical manifold modeling of porous rock with dominant fractures

Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.     

超临界CO2分段多簇压裂井间干扰规律研究

水平井分段多簇压裂是非常规油气藏开发的关键技术,在合理利用压裂诱导应力增大储层改造体积的同时,避免井间干扰所导致的裂缝砂堵和压裂窜扰,是压裂工艺优化中的关键科学问题。文章针对超临界CO₂分段多簇压裂的缝间干扰和井间干扰问题,采用流固耦合的扩展有限元方法研究单井及多井裂缝诱导应力演化特征,充分考虑裂缝内超临界CO₂流动和滤失,从非常规油气储层岩性特征、地应力场分布及施工工艺等多方面对压裂扰动应力进行系统研究,揭示单井分段多簇压裂缝扩展机制及应力扰动特征,在此基础上研究多井井间压裂缝干扰规律。结果表明：高水平应力差、高弹性模量的储层中压裂干扰界限较大,低水平应力差、低弹性模量地层需适度增大簇间距,减小簇间干扰;老井压裂后,其邻井压裂缝非对称系数随井间距呈先增大、后减小的趋势;当井间距等于压裂干扰界限时,非对称系数λ达到最大,且井周改造范围最大,但裂缝两翼的非对称性可能导致储层动用不充分。本研究为水平井细分切割压裂和立体式井网设计优化提供理论基础,在“双碳”战略背景下对非常规油气资源高效开发具有重要意义。     

The effect of stress boundary conditions on fluid-driven fracture propagation in porous media using a phase-field modeling approach

A phase-field approach for fluid-driven fracture propagation in porous media with varying constant compatible stress boundary conditions is discussed and implemented. Since crack opening displacement, fracture path, and stress values near the fracture are highly dependent on the given boundary conditions, it is crucial to take into account the impact of in situ stresses on fracturing propagation for realistic applications. We illustrate several numerical examples that include the effects of different boundary conditions on the fracture propagation. In addition, an example using realistic boundary conditions from a reservoir simulator is included to show the capabilities of our computational framework.     

Numerical analysis of thermal fracturing in subsurface cold water injection by finite element methods

Integration of poromechanics and fracture mechanics plays an important role in understanding a series of thermal fracturing phenomena in subsurface porous media such as cold water flooding for enhanced oil recovery, produced‐water reinjection for waste disposal, cold water injection for geothermal energy extraction, and CO₂ injection for geosequestration. Thermal fracturing modeling is important to prevent the potential risks when fractures propagate into undesired zones, and it involves the coupling of heat transfer, mass transport, and stress change as well as the fracture propagation. Analytical method, finite element method, and finite difference method as well as boundary element method have been used to perform the thermal fracturing modeling considering different degrees and combinations of coupling. In this paper, extended finite element method is employed for the thermal fracturing modeling in a fully coupled fashion with remeshing avoided, and the stabilized finite element method is employed to account for the convection‐dominated heat transfer in the fracturing process with numerical oscillation circumvented. With the thermal fracturing model, a hypothetical numerical experiment on cold water injection into a deep warm aquifer is conducted. Results show that parameters such as injection rate, injection temperature, aquifer stiffness, and permeability can affect the fracture development in different ways and extended finite element method and stabilized finite element method provide effective tools for thermal fracturing simulation. Copyright © 2012 John Wiley & Sons, Ltd.     

Modeling fractures as interfaces with nonmatching grids

We consider a model for fluid flow in a porous medium with a fracture. In this model, the fracture is treated as an interface between subdomains, on which specific equations have to be solved. In this article, we analyze the discrete problem, assuming that the fracture mesh and the subdomain meshes are completely independent, but that the geometry of the fracture is respected. We show that despite this nonconformity, first-order convergence is preserved with the lowest-order Raviart–Thomas(-Nedelec) mixed finite elements. Numerical simulations confirm this result.     

A discontinuous Galerkin method for two-phase flow in a porous medium enforcing H(div) velocityand continuous capillary pressure

We consider the slightly compressible two-phase flow problem in a porous medium with capillary pressure. The problem is solved using the implicit pressure, explicit saturation (IMPES) method, and the convergence is accelerated with iterative coupling of the equations. We use discontinuous Galerkin to discretize both the pressure and saturation equations. We apply two improvements, which are projecting the flux to the mass conservative H(div)-space and penalizing the jump in capillary pressure in the saturation equation. We also discuss the need and use of slope limiters and the choice of primary variables in discretization. The methods are verified with two- and three-dimensional numerical examples. The results show that the modifications stabilize the method and improve the solution.     

Optimal design of hydraulic fracturing in porous media using the phase field fracture model coupled with genetic algorithm

We present a framework for the coupling of fluid-filled fracture propagation and a genetic inverse algorithm for optimizing hydraulic fracturing scenarios in porous media. Fracture propagations are described by employing a phase field approach, which treats fracture surfaces as diffusive zones rather than of interfaces. Performance of the coupled approach is provided with applications to numerical experiments related to maximizing production or reservoir history matching for emphasizing the capability of the framework.     

Accelerating the convergence of coupled geomechanical‐reservoir simulations

The pressure variations during the production of petroleum reservoir induce stress changes in and around the reservoir. Such changes of the stress state can induce marked deformation of geological structures for stress sensitive reservoirs as chalk or unconsolidated sand reservoirs. The compaction of those reservoirs during depletion affects the pressure field and so the reservoir productivity. Therefore, the evaluation of the geomechanical effects requires to solve in a coupling way the geomechanical problem and the reservoir multiphase fluid flow problem. In this paper, we formulate the coupled geomechanical‐reservoir problem as a non‐linear fixed point problem and improve the resolution of the coupling problem by comparing in terms of robustness and convergence different algorithms. We study two accelerated algorithms which are much more robust and faster than the conventional staggered algorithm and we conclude that they should be used for the iterative resolution of coupled reservoir‐geomechanical problem. Copyright © 2006 John Wiley & Sons, Ltd.