Transient response of a circular cavity in a poroelastic medium

This paper considers the transient response of a pressurized long cylindrical cavity in an infinite poroelastic medium. To obtain transient solutions, Biot's equations for poroelastodynamics are specialized for this problem. A set of exact general solutions for radial displacement, stresses, pore pressure and discharge are derived in the Laplace transform space by using analytical techniques. Solutions are presented for three different types of prescribed transient radial pressures acting on the surface of a permeable as well as an impermeable cavity surface. Time domain solutions are obtained by inverting Laplace domain solutions using a reliable numerical scheme. A detailed parametric study is presented to illustrate the influence of poroelastic material parameters and hydraulic boundary conditions on the response of the medium. Comparisons are also presented with the corresponding ideal elastic solutions to portray the poroelastic effects. It is noted that the maximum radial displacement and hoop stress at the cavity surface are substantially higher than the classical static solutions and differ considerably from the transient elastic solutions. Time histories and radial variations of displacement, hoop stress, pore pressure and fluid discharge corresponding to a cavity in two representative poroelastic materials are also presented.     

Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface

By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Thermal effects on fluid flow and hydraulic fracturing from wellbores and cavities in low-permeability formations

The coupled heat-fluid-stress problem of circular wellbore or spherical cavity subjected to a constant temperature change and a constant fluid flow rate is considered. Transient analytical solutions for temperature, pore pressure and stress are developed by coupling conductive heat transfer with Darcy fluid flow in a poroelastic medium. They are applicable to low permeability porous media suitable for liquid-waste disposal and also simulating reservoir for enhanced oil recovery, where conduction dominates the heat transfer process. A full range of solutions is presented showing separately the effects of temperature and fluid flow on pore pressure and stress development. It is shown that injection of warm fluid can be used to restrict fracture development around wellbores and cavities and generally to optimize a fluid injection operation. Both the limitations of the solutions and the convective flow effect are addressed. Copyright © 1999 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Modelling poroelastic hollow cylinder experiments with realistic boundary conditions

A general poroelastic solution for axisymmetrical plane strain problems with time dependent boundary conditions is developed in Laplace domain. Time‐domain results are obtained using numerical inversion of the Laplace transform. Previously published solutions can be considered as special cases of the proposed solution. In particular, we could reproduce numerical results for solid and hollow poroelastic cylinders with suddenly applied load/pressure (Rice and Cleary, Rev. Geophys. Space Phys. 1976; 14 :227; Schmitt, Tait and Spann, Int. J. Rock Mech. Min. Sci. 1993; 30 :1057; Cui and Abousleiman, ASCE J. Eng. Mech. 2001; 127 :391). The new solution is used to model laboratory tests on thick‐walled hollow cylinders of Berea sandstone subjected to intensive pressure drawdown. In the experiments, pressure at the inner boundary of the hollow cylinder is observed to decline exponentially with a decay constant of 3–5 1/s. It is found that solutions with idealized step‐function type inner boundary conditions overestimate the induced tensile radial stresses considerably. Although basic poroelastic phenomena can be modelled properly at long time following a stepwise change in pressure, realistic time varying boundary conditions predict actual rock behaviour better at early time. Experimentally observed axial stresses can be matched but appear to require different values for α and ν than are measured at long time. The proposed solution can be used to calculate the stress and pore pressure distributions around boreholes under infinite/finite boundary conditions. Prospective applications include investigating the effect of gradually changing pore pressure, modelling open‐hole cavity completions, and describing the phenomenon of wellbore collapse (bridging) during oil or gas blowouts. Copyright © 2004 John Wiley & Sons, Ltd.     

An analytical approach for cased borehole stress calculation in general anisotropic formations

We present an explicit analytical approach for calculating the distribution of stress around a circular cased borehole in a general anisotropic formation. This approach is developed based on a pure elastic cased borehole model, which can include multiple cemented concentric casing strings that are surrounded by a homogeneous anisotropic medium. The wall of the innermost casing is under the compression of borehole internal fluid pressure, and the formation can be subjected to arbitrary far‐field stresses. Every casing–cement boundary is assumed to be welded, so as the cement–formation boundary. The derived analytical solution can be applied to an arbitrary wellbore trajectory in an arbitrary anisotropic formation and is also applicable to degenerate isotropic formations. Accuracy and robustness of the analytical solution is validated through comparing its results with those calculated from the classical Kirsch solution for an isotropic example that has identical elastic properties for all materials and from a finite element method for an anisotropic cased borehole example.     

Deformation of a poroelastic layer overlying an elastic half‐space due to dip‐slip faulting

An analytical solution of the plane strain problem of the deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half‐space due to two‐dimensional seismic sources buried in the elastic half‐space has been obtained. The integral expressions for the displacements, stresses and pore pressure have been obtained using the stress function approach by applying suitable boundary conditions at the free surface and the interface. The solution obtained is in the Laplace–Fourier transform domain. The case of a vertical dip‐slip line dislocation for the oceanic crust model of Earth is studied in detail. Schapery's formula is used for the Laplace inversion and the extended Simpson's formula for the Fourier inversion. Diffusion of pore pressure in the layer is studied numerically. Contour maps showing the pore pressure in the poroelastic layer have been plotted. The effect of the compressibility of the solid and fluid constituents on pore pressure has also been studied. Copyright © 2015 John Wiley & Sons, Ltd.     

Analysis of refracturing in horizontal wells: Insights from the poroelastic displacement discontinuity method

In this paper, a fully coupled 2‐dimensional poroelastic displacement discontinuity method is used to investigate the refracturing process in horizontal wells. One of the objectives of refracturing is to access new reserves by adding new hydraulic fractures in zones that were bypassed in the initial fracturing attempt. Pore pressure depletion in the vicinity of old fractures directly affects the state of stress and eventually the propagation of newly created hydraulic fractures. Thus, a poroelastic analysis is required to identify guidelines for the refracturing process, in particular to understand the extension of the pore pressure depletion, and eventually, the orientation of new as well as old fractures. We propose a fully coupled approach to model the whole process of child fracture propagation in a depleted area between 2 parent fractures in the same wellbore. This approach omits the need of using multistep workflow that is regularly used to model the process. The maximum tensile stress criterion (σ criterion) is used for hydraulic fracture propagation. The proposed method is verified using available analytical solutions for total stress and pore pressure loading modes on a line fracture in drained and undrained conditions. Then, test cases of multifractured horizontal wells are studied to calculate the time evolution of the stress and pore pressure fields around old fractures and to understand the effect of these fields on the propagation path of newly created fractures. Finally, the effect of the pore pressure depletion on the propagation path of the newly created fractures in the bypassed area of the same wellbore is studied. The results show that the depleted areas around old fractures are highly affected by the extent and severity of the stress redistribution and pore pressure depletion. It is observed that a successful creation of new fractures may only happen in certain time frames. The results of this study provide new insights on the behavior of newly created fractures in depleted zones. They also clarify the relationship between stress change and pore pressure depletion in horizontal wells.     

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamic response to fluid extraction from a poroelastic half‐space

In this study, the dynamic response of a poroelastic half‐space to a point fluid sink is investigated using Biot's dynamic theory of poroelasticity. Based on Biot's theory, the governing field equations are re‐formulated in frequency domain with solid displacement and pore pressure. In a cylindrical coordinate system, a method of displacement potentials for axisymmetric displacement field is proposed to decouple the Biot's field equations to three scalar Helmholtz equations, and then the general solution to axisymmetric problems are obtained. The full‐space fundamental singular solution for a point sink is also derived using potential methods. The mirror‐image method is finally applied to construct the fundamental solution for a point sink buried in a poroelastic half‐space. Furthermore, a numerical study is conducted for a rock, that is, Berea sandstone, as a representative example. Copyright © 2013 John Wiley & Sons, Ltd.     

Analysis of Pore Pressure and Stress Distribution around a Wellbore Drilled in Chemically Active Elastoplastic Formations

Drilling in low-permeable reactive shale formations with water-based drilling mud presents significant challenges, particularly in high-pressure and high-temperature environments. In previous studies, several models were proposed to describe the thermodynamic behaviour of shale. Most shale formations under high pressure are expected to undergo plastic deformation. An innovative algorithm including work hardening is proposed in the framework of thermo-chemo-poroelasticity to investigate the effect of plasticity on stresses around the wellbore. For this purpose a finite-element model of coupled thermo-chemo-poro-elastoplasticity is developed. The governing equations are based on the concept of thermodynamics of irreversible processes in discontinuous systems. In order to solve the plastic problem, a single-step backward Euler algorithm containing a yield surface-correction scheme is used to integrate the plastic stress–strain relation. An initial stress method is employed to solve the non-linearity of the plastic equation. In addition, super convergent patch recovery is used to accurately evaluate the time-dependent stress tensor from nodal displacement. The results of this study reveal that thermal and chemical osmosis can significantly affect the fluid flow in low-permeable shale formations. When the salinity of drilling mud is higher than that of pore fluid, fluid is pulled out of the formation by chemical osmotic back flow. Similar results are observed when the temperature of drilling mud is lower than that of the formation fluid. It is found that linear elastic approaches to wellbore stability analysis appear to overestimate the tangential stress around the wellbore and produce more conservative stresses compared to the results of field observation. Therefore, the drilling mud properties obtained from the elastoplastic wellbore stability in shales provide a safer mud weight window and reduce drilling cost.     

Finite element analysis of land subsidence above depleted reservoirs with pore pressure gradient and total stress formulations

The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd.     

Vertical vibrations of a rigid disk embedded in a poroelastic medium

This paper considers the steady-state vertical vibrations of a rigid circular disk embedded at a finite depth below the free surface of a poroelastic medium. Biot's elastodynamic theory for porous media is used in the analysis. General solutions for axisymmetric poroelastic fields are obtained by using Hankel integral transforms. Analytical solutions for influence functions corresponding to four types of buried axisymmetric excitations are derived. The embedded disk problem is fomulated in terms of a set of coupled integral equations for unknown traction and pore pressure jumps across the disk. The kernel functions of the integral equations are the influence functions corresponding to buried vertical, radial and pore pressure ring loads. The system of integral equations is solved numerically by discretizing the disk into several concentric annular rings. Selected numerical solutions for displacements, vertical stress and pore pressure due to a buried fully flexible disk (uniform pressure) are also presented. The vertical compliances of a rigid disk are examined for different depths of embedment, poroelastic materials and hydraulic boundary conditions. Solutions for traction and pore pressure jumps are also examined. The present results are useful in the study of dynamic response of embedded foundations and anchors in poroelastic soils. Copyright © 1999 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 Fully Coupled Chemo-Poroelastic Analysis of Pore Pressure and Stress Distribution around a Wellbore in Water Active Rocks

Water active rocks consist of minerals that hold water in their crystalline structure and in pore spaces. Free water from drilling fluid can be attracted by the formation depending on the potential differences between pore space and drilling fluid. The fluid movement into the formation or out of the formation can lead to a change in effective stress, thus causing wellbore failures. In all previous studies it is found that the solute transport from or to the formation is primarily controlled by diffusion process and the effect of advection on solute transfer is negligible for a range of very low permeable shale formations (>10⁻⁵ mD). In this study a range of permeable shale formations (10⁻⁵ to 10⁻³ mD) commonly encountered in drilling oil and gas wells are considered to investigate the solute transfer between drilling fluid and formation due to advection. For this purpose a finite element model of fully coupled chemo-hydro-mechanical processes was developed. Results of this study revealed that the solute transfer between the drilling fluid and the shale formation is controlled primarily by permeability of the shale formations. For the range of shale formations studied here, there exists a threshold permeability below which the solute transfer is dominated by diffusion process and above which by fluid in motion (fluid flow). Results from the numerical experiments have shown that when the permeability of shales is greater than this threshold permeability, the chemical potential gradient between the pore fluid and drilling fluid reaches equilibrium faster than that when the permeability of shales is below this threshold value. Also it has been found that when advection is taken into account, effective radial and tangential stresses decrease around the wellbore, particularly near the wellbore wall where the solute concentration has reached near equilibrium.     

横观各向同性饱和地基三维瞬态Lamb问题

研究了横观各向同性饱和土地基在地表动力荷载作用下的三维瞬态响应。基于饱和多孔介质的三维Biot波动理论,利用Laplace变换,建立圆柱坐标系下横观各向同性饱和土的波动方程;解耦波动方程后,根据算子理论,并借助Fourier展开和Hankel变换技术,得到瞬态荷载作用下,饱和土介质的土骨架位移和应力、孔隙水相对位移和孔隙水压力的一般解;利用一般解,给出横观各向同性饱和地基在地表集中荷载激励下的瞬态Lamb问题的解答。数值算例结果表明,采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的瞬态动力特性。     

Analytical solutions of a finite two‐dimensional fluid‐saturated poroelastic medium with compressible constituents

An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd.     

Borehole Stability Analysis of Horizontal Drilling in Shale Gas Reservoirs

Serious wellbore instability occurs frequently during horizontal drilling in shale gas reservoirs. The conventional forecast model of in situ stresses is not suitable for wellbore stability analysis in laminated shale gas formations because of the inhomogeneous mechanical properties of shale. In this study, a new prediction method is developed to calculate the in situ stresses in shale formations. The pore pressure near the borehole is heterogeneous along both the radial and tangential directions due to the inhomogeneity in the mechanical properties and permeability. Therefore, the stress state around the wellbore will vary with time after the formation is drained. Besides, based on the experimental results, a failure criterion is verified and applied to determine the strength of Silurian shale in the Sichuan Basin, including the long-term strength of gas shale. Based on this work, horizontal well borehole stability is analyzed by the new in situ stress prediction model. Finally, the results show that the collapse pressure will be underestimated if the conventional model is used in shale gas reservoirs improperly. The collapse pressure of a horizontal well is maximum at dip angle of 45°. The critical mud weight should be increased constantly to prevent borehole collapse if the borehole is exposed for some time.     

Bending of fluid‐saturated linear poroelastic beams with compressible constituents

Analytical solutions are presented for fluid‐saturated linear poroelastic beams under pure bending. The stress‐free boundary condition at the lateral surfaces is satisfied in the St Venant's sense and the Beltrami–Michell compatibility conditions are resolved rigorously, rendering the flexure of the beams analytically tractable. Two sets of formulations are derived based on the coupled and uncoupled diffusion equations respectively. The analytical solutions are compared with three‐dimensional finite element simulations. Both sets of analytical formulations are capable of capturing exactly both the initial (undrained) and the steady‐state (fully drained) deflection of the beams. However, the analytical solutions are found to be deficient during the transient phase. The cause for the deficiency of the transient analytical solutions is discussed. The accuracy of the analytical solutions improves as Poisson's ratio and the compressibility of the constituents of the porous beam increase, where the St Venant's edge effect at the lateral surfaces is mitigated. Copyright © 2008 John Wiley & Sons, Ltd.     

Poro‐elasto‐plastic response of an unconsolidated formation confined with stiff seal rocks under radial injection

In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.