A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

A finite difference method for earthquake sequences in poroelastic solids

Induced seismicity (earthquakes caused by injection or extraction of fluids in Earth’s subsurface) is a major, new hazard in the USA, the Netherlands, and other countries, with vast economic consequences if not properly managed. Addressing this problem requires development of predictive simulations of how fluid-saturated solids containing frictional faults respond to fluid injection/extraction. Here, we present a finite difference method for 2D linear poroelasticity with rate-and-state friction faults, accounting for spatially variable properties. Semi-discrete stability and accuracy are proven using the summation-by-parts, simultaneous-approximation-term (SBP-SAT) framework for discretization and boundary condition enforcement. Convergence rates are verified using the method of manufactured solutions and comparison to the analytical solution to Mandel’s problem. The method is then applied to study fault slip triggered by fluid injection and diffusion through high-permeability fault damage zones. We demonstrate that in response to the same, gradual forcing, fault slip can occur in either an unstable manner, as short-duration earthquakes that radiate seismic waves, or as stable, aseismic, slow slip that accumulates over much longer time scales. Finally, we use these simulation results to discuss the role of frictional and elastic properties in determining the stability and nature of slip.     

Dynamic response of a soft soil layer to flow and periodical disturbance

The dynamic response of a soft soil layer of finite thickness under the mutual effects of flow and periodical disturbance at the free surface is discussed in this work. The homogeneous water is governed by potential theory and the soil layer obeys Biot's theory of poroelasticity. The boundary‐value problem is solved by an analytical algorithm, in which the wave number is found first. Secondly, the closed form solutions are found by a two‐parameter perturbation method with the boundary‐layer correction. The results are also compared with those of the poroelastic soil layer of infinite thickness to show the impermeable rigid boundary effect. Copyright © 2003 John Wiley & Sons, Ltd.     

A direct boundary element method for plane strain poroelasticity

This paper presents a direct boundary element method (BEM), formulated in the Laplace transform space, for plane strain poroelasticity. The paper expands on work by Cheng and Liggett by recasting the theoretical foundation of BEM within the framework of Rice and Cleary's formulation of the Biot theory of poroelasticity. Furthermore, the numerical algorithm is generalized to deal with both interior and exterior domain problems, and a method for indirectly calculating the Cauchy principal value of the singular integrals is presented. Formulae for the stress and flux inside the domain are also derived. Finally, the algorithm is validated by comparing the numerical results with the analytical solution of a borehole subject to a far-field deviatoric stress (exterior domain) and with the solution of Mandel's problem (interior domain). These two examples provide a critical test of the algorithm.     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.     

Pore‐scale modeling of fluid‐particles interaction and emerging poromechanical effects

A micro‐hydromechanical model for granular materials is presented. It combines the discrete element method for the modeling of the solid phase and a pore‐scale finite volume formulation for the flow of an incompressible pore fluid. The coupling equations are derived and contrasted against the equations of conventional poroelasticity. An analogy is found between the discrete element method pore‐scale finite volume coupling and Biot's theory in the limit case of incompressible phases. The simulation of an oedometer test validates the coupling scheme and demonstrates the ability of the model to capture strong poromechanical effects. A detailed analysis of microscale strain and stress confirms the analogy with poroelasticity. An immersed deposition problem is finally simulated and shows the potential of the method to handle phase transitions. Copyright © 2013 John Wiley & Sons, Ltd.     

FINITE ELEMENT ANALYSES OF ANISOTROPIC POROELASTICITY: A GENERALIZED MANDEL'S PROBLEM AND AN INCLINED BOREHOLE PROBLEM

The finite element equations for non-linear, anisotropic poroelasticity are cast in the form of measurable engineering constants. Two problems of importance to the rock and petroleum industry are analysed by the FEM. First, the classical Mandel's problem with an extension to transversely isotropic case is investigated. Second, the problem of an inclined borehole is explored. In particular, the effect of material anisotropy on stress concentration near the wall with implication to borehole instability is examined in detail.     

Explanation of locking of four‐node plane element by considering it as elastic Dirichlet‐type boundary value problem

The aim of this study is to arrive at a better understanding of the phenomenon of locking of low‐order compatible displacement type of finite elements in particular for the hour‐glass mode of the plane four‐node element and dilative materials. To this end the properties of finite elements are investigated in an analytical way, where a finite element is considered as a plane boundary value problem with prescribed boundary displacement (Dirichlet problem). In this paper for the sake of simplicity the simplest possible linear comparison solid, namely isotropic linear elasticity, is applied, although recognizing fully that for a dilative material elasto‐plasticity would be more realistic. From the study described in this paper it is concluded that locking of the four‐node element is not due to any particular numerical formulation of this compatible finite element since, even the analytical solution suffers from this problem. The locking of this element is not related to incompressibility of the material either as the analytical solution shows locking to occur at a parameter set which differs significantly from the one in case of incompressibility. It is shown that locking is a consequence of the combination of the dilative material behaviour and the compatible displacement type of boundary conditions, which leads to infinite isotropic stresses in the element. These infinite isotropic stresses occur at the limit of uniqueness of the solution, which for this element is shown to occur outside the parameter range of the sufficiency of uniqueness. Copyright © 2000 John Wiley & Sons, Ltd.     

A stabilized assumed deformation gradient finite element formulation for strongly coupled poromechanical simulations at finite strain

An adaptively stabilized finite element scheme is proposed for a strongly coupled hydro‐mechanical problem in fluid‐infiltrating porous solids at finite strain. We first present the derivation of the poromechanics model via mixture theory in large deformation. By exploiting assumed deformation gradient techniques, we develop a numerical procedure capable of simultaneously curing the multiple‐locking phenomena related to shear failure, incompressibility imposed by pore fluid, and/or incompressible solid skeleton and produce solutions that satisfy the inf‐sup condition. The template‐based generic programming and automatic differentiation (AD) techniques used to implement the stabilized model are also highlighted. Finally, numerical examples are given to show the versatility and efficiency of this model. Copyright © 2013 John Wiley & Sons, Ltd.     

A supercoarsening multigrid method for poroelasticity in 3D coupled flow and geomechanics modeling

The Galerkin finite-element discretization of the force balance equation typically leads to large linear systems for geomechanical problems with realistic dimensions. In iteratively coupled flow and geomechanics modeling, a large linear system is solved at every timestep often multiple times during coupling iterations. The iterative solution of the linear system stemming from the poroelasticity equations constitutes the most time-consuming and memory-intensive component of coupled modeling. Block Jacobi, LSOR, and Incomplete LU factorization are popular preconditioning techniques used for accelerating the iterative solution of the poroelasticity linear systems. However, the need for more effective, efficient, and robust iterative solution techniques still remains especially for large coupled modeling problems requiring the solution of the poroelasticity system for a large number of timesteps. We developed a supercoarsening multigrid method (SCMG) which can be multiplicatively combined with commonly used preconditioning techniques. SCMG has been tested on a variety of coupled flow and geomechanics problems involving single-phase depletion and multiphase displacement of in-situ hydrocarbons, CO₂ injection, and extreme material property contrasts. Our analysis indicates that the SCMG consistently improves the convergence properties of the linear systems arising from the poroelasticity equations, and thus, accelerates the coupled simulations for all cases subject to investigation. The joint utilization of the two-level SCMG with the ILU1 preconditioner emerges as the most optimal preconditioning/iterative solution strategy in a great majority of the problems evaluated in this work. The BiCGSTAB iterative solver converges more rapidly compared to PCG in a number of test cases, in which various SCMG-accelerated preconditioning strategies are applied to both iterators.     

A stabilized Smoothed Particle Hydrodynamics,Taylor–Galerkin algorithm for soil dynamics problems

Modelling of failure under dynamic conditions in geomaterials with finite elements presents a series of complex problems, among which we can mention those of (i) volumetric locking, which results on higher failure loads, (ii) influence of mesh alignment, resulting to unrealistic failure surfaces, (iii) diffusion of the shear band over some element widths, (iv) nonoptimal propagation properties (numerical diffusion and dispersion), (v) fulfilling Babuska–Brezzi conditions when using the same order of interpolation for displacement and pressures in coupled problems and (vi) large deformation analysis. This paper is based on previous work done by the authors, who developed a mixed approximation based on (i) casting the dynamic problem in the form of a system of first order PDEs and (ii) using stresses and velocities as nodal variables. The equations were discretized following a Taylor–Galerkin algorithm, first in time using a Taylor expansion and then in space using Galerkin method. The model was limited to small deformations. The purpose of this paper is to show how Taylor–Galerkin method can be extended to meshless formulations, such as the SPH method. The algorithm consists of (i) discretizing in time using a Taylor series expansion complemented with integration of source terms using a Runge–Kutta scheme and then (ii) discretizing in space using the SPH method. It is shown how the proposed method keeps the advantages of the Taylor–Galerkin method in Finite Elements (good propagation properties and capturing of shear bands) and avoids the tensile instability. A set of test problems ranging from elastic propagation of a wave in a bar to failure of a slope on a cohesive softening material are used to assess the performance of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

A note on the consolidation settlement of a rigid circular foundation on a poroelastic halfspace

This paper uses Biot's poroelasticity approach to examine the consolidation behaviour of a rigid foundation with a frictionless base in contact with a poroelastic halfspace. The mathematical development of the mixed boundary value problem involves a set of dual integral equations in the Laplace transform domain which cannot be conveniently solved by employing conventional procedures. In this paper, a numerical solution is developed using a scheme where the contact normal stress is approximated by a discretized equivalent. The influence of limiting drainage boundary conditions at the entire surface of the halfspace on the degree of consolidation of the rigid circular foundation is investigated. The results obtained in this study are compared with the corresponding results given in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Plain strain problem of poroelasticity using eigenvalue approach

A plain strain problem of an isotropic elastic liquid-saturated porous medium in poroelasticity has been studied. The eigenvalue approach using the Laplace and Fourier transforms has been employed and these transforms have been inverted by using a numerical technique. An application of infinite space with concentrated force at the origin has been presented to illustrate the utility of the approach. The displacement and stress components in the physical domain are obtained numerically. The results are shown graphically and can be used for a broad class of problems related to liquid-saturated porous media.     

尾矿连续沉积过程中超孔隙水压力演化的温度-渗流-力学-化学耦合模型

尾矿胶结充填技术可避免矿渣的地表堆存、改善采场围岩稳定性、提高矿石回采率,因此在国内外的矿山开采中得到了广泛应用。充填系统的稳定性取决于尾矿连续充填过程中多场耦合作用造成的孔隙水压力演化。本研究基于经典Biot孔隙弹性理论建立尾矿的温度-渗流-力学-化学耦合模型框架,进而提出连续沉积过程中尾矿热-化学固结的一维超压模型,并推导化学反应造成水压变化的临界温度闭合公式。通过分析不同沉积速率、不同初始和边界温度条件下的充填过程,揭示了多场耦合作用对尾矿超孔隙水压力的影响机理。     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

A high‐resolution local RBF collocation method for steady‐state poroelasticity and hydromechanical damage analysis

In this work, we describe a meshless numerical method based on local collocation with RBFs for the solution of the poroelasticity equation. The RBF finite collocation approach forms a series of overlapping nodal stencils, over which an RBF collocation is performed. These local collocation systems enforce the governing PDE operator throughout their interior, with the intersystem communication occurring via the collocation of field variables at the stencil periphery. The method does not rely on a generalised finite differencing approach, whereby the governing partial differential operator is reconstructed at the global level to drive the solution of the PDE. Instead, the PDE governing and boundary operators are enforced directly within the local RBF collocation systems, and the sparse global assembly is formed by reconstructing the value of the field variables at the centrepoint of the local stencils. In this way, the solution of the PDE is driven entirely by the local RBF collocation, and the method more closely resembles the approach of the full‐domain RBF collocation method. By formulating the problem in this fashion, high rates of convergence may be attained without the computational cost and numerical ill‐conditioning issues that are associated with the full‐domain RBF collocation approach. An analytical solution is formulated for a 2D poroelastic fluid injection scenario and is used to verify the proposed implementation of the method. Highly accurate solutions are produced, and convergence rates in excess of sixth order are observed for each field variable (i.e. pressure and displacement) and field‐variable derivative (i.e. pressure gradients and stresses). The stress and displacement fields resulting from the solution of the poroelasticity equation are then used to describe the formation and propagation of microfractures and microfissures, which may form in the presence of large shear strain, in terms of a continuous damage variable which modifies the mechanical and hydraulic properties of the porous medium. The formation of such hydromechanical damage, and the resulting increase in hydraulic conductivity, is investigated for a pressurised injection into sandstone. Copyright © 2014 John Wiley & Sons, Ltd.     

Extension of the Griffith's fracture criteria to saturated clays

Inglis [1] has solved the problem of distribution of stress in an elastic plate around an elliptical hole. His works clarify the role of cracks in the failure of an elastic material. However, his solution cannot be applied to saturated clay because he considers only total stresses, while, in saturated clay, the criterion of rupture should be expressed in terms of effective and not total stresses. The solution of Atkinson and Craster [2] using Biot's poroelasticity theory, shows that there is no high pore pressure in the vicinity of the crack tips for saturated clay. The major difference between this approach and the Biot's theory of is that, in saturated clay, strain is a function of the variation of the effective stress [3], while, in poroelastic media, strain is only a function of the variation of the total stress [4, Equation 2.2]. Also in their solution there is continuity between the pore fluid and the inner fluid in the crack. Their solution is valid for poroelastic media involving a movement of the pore fluid. In our solution there is no movement of the pore fluid (Undrained condition). In this paper we have solved the same problem as Inglis [1], but for the particular case of saturated clay obeying elastic law. By solving this problem we obtained the expressions for pore pressure, effective stress, total stress and displacements. The results show that not only the total stress but also the pore pressure and the effective stress are also high in the vicinity of the crack tips. A new failure criterion, based on Griffith's strain energy principle [5] and maximum tensile stress [6], valid for saturated clay is developed in this paper. Copyright © 2003 John Wiley & Sons, Ltd.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

A multiscale finite element formulation for axisymmetric elastoplasticity with volumetric locking

Similar to plane strain, axisymmetric stress problem is also highly kinematics constrained. Standard displacement‐based finite element exhibits volumetric locking issue in simulating nearly/fully incompressible material or isochoric plasticity under axisymmetric loading conditions, which severely underestimates the deformation and overestimates the bearing capacity for structural/geotechnical engineering problems. The aim of this paper is to apply variational multiscale method to produce a stabilized mixed displacement–pressure formulation, which can effectively alleviate the volumetric locking issue for axisymmetric stress problem. Both nearly incompressible elasticity and isochoric J₂ elastoplasticity are investigated. First‐order 3‐node triangular and 4‐node quadrilateral elements are tested for locking issues. Several representative simulations are provided to demonstrate the performance of the linear elements, which include the convergence study and comparison with closed‐form solutions. A comparative study with pressure Laplacian stabilized formulation is also presented. Copyright © 2009 John Wiley & Sons, Ltd.     

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.