Transient Stresses and Displacement Around a Wellbore Due to Fluid Flow in Transversely Isotropic, Porous Media: I. Infinite Reservoirs

1968 ◽  
Vol 8 (01) ◽  
pp. 63-78 ◽  
Author(s):  
M.S. Seth ◽  
K.E. Gray

Abstract Equations of elasticity for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure were formulated. Using an analogy between thermal and porous body stresses, it was shown that the solution for a transversely isotropic porous body may be obtained by incorporating body forces and the stresses due to a boundary load into the corresponding solution for the thermal stress problem. Equations of elasticity and the thermal analogy method were used to determine transient horizontal, tangential, and vertical s tresses and radial displacement in a semi-infinite cylindrical region when either a constant pressure or a constant rate of flow is maintained at the wellbore. The vertical and tangential displacements are zero from the conditions of the problem. A numerical analysis was made of the solutions obtained by using a digital computer to determine the relative influence of each system variable. Considering rock as a porous body with internal fluid pressure generally gives results significantly different than considering the rock to be nonporous; the directional character of rocks leads to significant differences as compared to results based upon the common assumption of isotropy. Stress gradients are high near the wellbore but die out away from the well. Radial stresses are compressive or neutral, whereas tangential stresses are tensile, neutral or compressive, depending upon the boundary conditions and physical properties of the system. Vertical stresses are compressive for an unbounded system. For constant wellbore injection rate, the vertical stress is proportional to the rate of fluid injection and decreases with time, whereas the radial and tangential stresses increase with time. At a given location, the radial displacement generally is very dependent upon time. INTRODUCTION A realistic appraisal of the state of stress in subsurface rock formations would be of considerable interest and use to the petroleum industry. For example, knowing the state of stress in proximity to a wellbore would be of fundamental importance in designing a fracturing operation or, more important, of clearly understanding the conditions necessary to produce rock failure of desired dimensions and geometry. Understanding conditions necessary for rock failure at the wellbore would also be of utility in a preventive sense. For example, borehole stability is an important consideration for many rock formations, and knowledge of the stress state at and near the wellbore under conditions of substantial pressure gradient due to fluid flow would be of great value. When fluid flows through a porous body which is initially at some uniform stress level, the following forces generate stresses at any point in the body.Forces due to nonuniform pressure distribution. With increasing pressure the elements of a body are compressed. Such compression cannot proceed freely in a continuum when the pressure is not uniform throughout, and thus, stresses due to flow of fluid are set up.Pore fluid pressure. This gives rise to normal stresses whose value at any point is the product of areal porosity and the fluid pressure. Although the fluid exerts uniform pressure, the stresses it creates in an anisotropic body may not be the same in all directions since the areal porosity in an anisotropic porous body is a direction-dependent quantity. This consideration leads to the concept of directional porosity.

1968 ◽  
Vol 8 (01) ◽  
pp. 79-86 ◽  
Author(s):  
M.S. Seth ◽  
K.E. Gray

Abstract In Part 1 of this work,1 equations of elasticity were formulated for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure. Equations of elasticity and the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant rate of pressure or a constant rate of flow is maintained at the wellbore. In this paper, the approach presented earlier is extended to finite reservoirs for the cases ofsteady-state flow,constant pressures at the well bore and outer boundary andconstant pressure at the wellbore and no flow at the outer boundary. Results of this work show that radial and tangential stress gradients are high near the wellbore but diminish rapidly away from the well; the vertical stress gradient behaves in the same way but is less severe. Radial stresses are compressive or neutral, whereas tangential and vertical stresses may be tensile, neutral or compressive, depending upon the boundary conditions, the physical properties of the system and the radial distance involved (vertical stresses are always compressive in an unbounded system1). For constant boundary pressures, both radial and tangential stresses increase with time whereas they both decrease for a closed outer boundary and constant pressure at the wellbore. The vertical stress decreases with time for both systems. For steady-state systems, radial displacement may be positive or negative, depending upon the dimensions of the system, the pressure differential and the porosity. Radial displacement may be positive or negative for a closed outer boundary but is positive for constant pressures at both boundaries. INTRODUCTION The importance, utility and complexity of a realistic appraisal of the stress state at and local to a wellbore were indicated in Part 1. In this paper the analytical approach presented earlier is extended to finite, cylindrical reservoir geometry for the cases ofsteady-state flow,constant pressures at wellbore and outer boundary andconstant pressure at the wellbore and no flow at the outer boundary. Other than the outer boundary of the reservoir being finite, the physical system and assumptions pertinent thereto are the same as before. The reader may wish to review the mathematical development through Eq. 49 of Part 1 before proceeding here.


2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Minh H. Tran ◽  
Younane N. Abousleiman

The porochemoelectroelastic analytical models have been used to describe the response of chemically active and electrically charged saturated porous media such as clay soils, shales, and biological tissues. However, existing studies have ignored the anisotropic nature commonly observed on these porous media. In this work, the anisotropic porochemoelectroelastic theory is presented. Then, the solution for an inclined wellbore drilled in transversely isotropic shale formations subjected to anisotropic far-field stresses with time-dependent down-hole fluid pressure and fluid activity is derived. Numerical examples illustrating the combined effects of porochemoelectroelastic behavior and anisotropy on wellbore responses are also included. The analysis shows that ignoring either the porochemoelectroelastic effects or the formation anisotropy leads to inaccurate prediction of the near-wellbore pore pressure and effective stress distributions. Finally, wellbore responses during a leak-off test conducted soon after drilling are analyzed to demonstrate the versatility of the solution in simulating complex down-hole conditions.


1999 ◽  
Author(s):  
Gerard A. Ateshian ◽  
Michael A. Soltz

Abstract Experimental studies have demonstrated that the moduli of articular cartilage in compression are one to two orders of magnitude smaller than in tension. However, only a few analyses of cartilage mechanics have been performed which account for this tension-compression nonlinearity (Soulhat et al., 1998; Ateshian and Soltz, 1999a,b). In order to understand the state of stress under loading conditions which simulate the physiologic environment of diarthrodial joints, and the possible implications for tissue failure and the pathomechanics of osteoarthritis, it is important to determine whether the tension-compression nonlinearity of cartilage significantly affects our current understanding of its response in contact mechanics. Most analyses of cartilage contact have employed linear isotropic or transversely isotropic models for cartilage, either within the context of elasticity theory or porous media theories. In this study, we present a finite element solution for the contact of a rigid spherical impermeable sphere against a cartilage layer supported on a rigid impermeable subchondral bone foundation, where cartilage is modeled using our recently proposed biphasic conewise linear elasticity model (Ateshian and Soltz, 1999a,b). A comparison is also provided with the more frequently used linear isotropic biphasic model, under similar conditions.


1974 ◽  
Vol 14 (03) ◽  
pp. 271-278 ◽  
Author(s):  
Milos Kojic ◽  
J.B. Cheatham

Introduction A number of problems occur in the fields of drilling and rock mechanics for which consideration must be given to the interaction of fluid flow and rock deformation. Such problems include those of borehole stability, chip removal from under a drill bit, drilling in the presence of a fluid pressure gradient between the drilling fluid and formation fluid, and drilling by use of hydraulic jets. We have recently developed a general theory of the influence of fluid pressure gradients and gravity on the plasticity of porous media. The solution of the problem considered here serves as an example of the application of that theory. The illustrative problem is to determine the load required on a flat problem is to determine the load required on a flat axially symmetric punch for incipient plasticity of the porous medium under the punch when fluid flows through the bottom face of the punch. The rock is assumed to behave as a Coulomb plastic material under the influence of body forces plastic material under the influence of body forces due to fluid pressure gradients and gravity. Numerical methods that have been used by Cox et al. for analyzing axially symmetric plastic deformation in soils with gravity force are applied to the problem considered here. Involved is an iterative process for determining the slip lines. The fluid flow field ‘used for calculating the fluid pressure gradient is based upon the work by Ham pressure gradient is based upon the work by Ham in his study of the potential distribution ahead of the bit in rotary drilling. The effective stresses in the porous rock and the punch force for incipient plasticity are computed in terms of the fluid plasticity are computed in terms of the fluid pressure and the cohesive strength and internal pressure and the cohesive strength and internal friction of the rock. PLASTICITY OF POROUS MEDIA PLASTICITY OF POROUS MEDIA A recently developed general theory of plasticity of porous media under the influence of fluid flow is summarized in this section. The equation of motion for the porous solid for the case of incipient plastic deformation reduces to the following equilibrium equation:(1) where Ts is the partial stress tensor of the solid; Fs is the body force acting on the solid per unit volume of the solid material; P is the interaction force between the solid and the fluid; and is the porosity, which is defined as the ratio of the pore porosity, which is defined as the ratio of the pore volume to the total volume of the solid-fluid mixture. The partial stress tensor Ts can be considered as the effective stress tensor that is used in sod mechanics. With the acceptance of the effective stress principle defined in Ref. 5, the yield function, f, in the following form is satisfied for plastic deformation of the porous medium. plastic deformation of the porous medium.(2) where EP is the plastic strain tensor and K and the work-hardening parameter. From the equation of motion for the fluid, the interaction force P can be expressed in the form(3) where is the inertial force of the fluid per unit volume of the mixture and F is the body force acting on the fluid per unit volume of fluid. For the case of incipient plastic deformation the solid can be considered static (velocities of the solid particles are zero), and the problem of determining particles are zero), and the problem of determining the fluid flow field is the one usually analyzed in petroleum engineering. petroleum engineering. Consider a flow of be fluid such that the inertial forces of the fluid can be neglected and assume that Darcy's law is applicable. SPEJ P. 271


2009 ◽  
Vol 131 (9) ◽  
Author(s):  
Lars Johansson ◽  
Ulf Edlund ◽  
Anna Fahlgren ◽  
Per Aspenberg

A model where bone resorption is driven by stimulus from fluid flow is developed and used as a basis for computer simulations, which are compared with experiments. Models for bone remodeling are usually based on the state of stress, strain, or energy density of the bone tissue as the stimulus for remodeling. We believe that there is experimental support for an additional pathway, where an increase in the amount of osteoclasts, and thus osteolysis, is caused by the time history of fluid flow velocity, fluid pressure, or other parameters related to fluid flow at the bone/soft tissue interface of the porosities in the bone.


2015 ◽  
Vol 2015 ◽  
pp. 1-10 ◽  
Author(s):  
Fei Liu ◽  
Li-qiang Zhao ◽  
Ping-li Liu ◽  
Zhi-feng Luo ◽  
Nian-yin Li ◽  
...  

This paper proposes a numerical model for the fluid flow in fractured porous media with the extended finite element method. The governing equations account for the fluid flow in the porous medium and the discrete natural fractures, as well as the fluid exchange between the fracture and the porous medium surrounding the fracture. The pore fluid pressure is continuous, while its derivatives are discontinuous on both sides of these high conductivity fractures. The pressure field is enriched by the absolute signed distance and appropriate asymptotic functions to capture the discontinuities in derivatives. The most important advantage of this method is that the domain can be partitioned as nonmatching grid without considering the presence of fractures. Arbitrarily multiple, kinking, branching, and intersecting fractures can be treated with the new approach. In particular, for propagating fractures, such as hydraulic fracturing or network volume fracturing in fissured reservoirs, this method can process the complex fluid leak-off behavior without remeshing. Numerical examples are presented to demonstrate the capability of the proposed method in saturated fractured porous media.


Sign in / Sign up

Export Citation Format

Share Document