The Effect of Geological and Geomechanical Parameters on Reservoir Stress Path and Its Importance in Studying Permeability Anisotropy

2000 ◽  
Vol 3 (05) ◽  
pp. 394-400 ◽  
Author(s):  
M. Khan ◽  
L.W. Teufel

Summary Reservoir stress path is defined as the ratio of change in effective horizontal stress to the change in effective vertical stress from initial reservoir conditions during pore-pressure drawdown. Measured stress paths of carbonate and sandstone reservoirs are always less than the total stress boundary condition (isotropic loading) and are either greater or less than the stress path predicted by the uniaxial strain boundary condition. Clearly, these two boundary-condition models that are commonly used by the petroleum industry to calculate changes in effective stresses in a reservoir and to measure reservoir properties in the laboratory are inaccurate and can be misleading if applied to reservoir management problems. A geomechanical model that incorporates geologic and geomechanical parameters was developed to more accurately predict the reservoir stress path. Numerical results show that reservoir stress path is dependent on the size and geometry of the reservoir and on elastic properties of the reservoir rock and bounding formations. In general, stress paths become lower as the aspect ratio of reservoir length to thickness increases. Lenticular sandstone reservoirs have a higher stress path than blanket sandstone reservoirs that are continuous across a basin. This effect is enhanced when the bounding formations have a lower elastic modulus than the reservoir and when the reservoir is transversely isotropic. In addition, laboratory experiments simulating reservoir depletion for different stress path conditions demonstrate that stress-induced permeability anisotropy evolves during pore-pressure drawdown. The maximum permeability direction is parallel to the maximum principal stress and the magnitude of permeability anisotropy increases at lower stress paths. Introduction Matrix permeability and pore volume compressibility are fundamentally important characteristics of hydrocarbon reservoirs because they provide measures of reservoir volume and reservoir producibility. Laboratory studies have shown that these properties are stress sensitive and are usually measured under hydrostatic (isotropic) loads that do not truly reflect the anisotropic stress state that exists in most reservoirs and do not adequately simulate the evolution of deviatoric stresses in a reservoir as the reservoir is produced. Recent laboratory studies1–3 have shown that permeability and compressibility are dependent on the deviatoric stress and change significantly with reservoir stress path. In-situ stress measurements in carbonate and clastic reservoirs indicate that the reservoir stress path is not isotropic loading (equal to 1.0) and can range from 0.14 to 0.76. 4 The measured reservoir stress paths are also inconsistent with the elastic uniaxial strain model5 commonly used to calculate horizontal stress and changes in horizontal stress with pore-pressure drawdown. The calculated uniaxial strain stress path can be significantly less or greater than the measured stress path.4 Knowledge of the stress path that reservoir rock will follow during production and how this stress path will affect reservoir properties is critical for reservoir management decisions necessary to increase reservoir producibility. However, in-situ stress measurements needed to determine reservoir stress path are difficult and expensive to conduct, and may take several years to collect. Various analytical models have been proposed to calculate in-situ horizontal stresses and they could be applied to the prediction of reservoir stress path during pore-pressure drawdown.5–9 However, none of these models addresses all of the essential geological and geomechanical factors that influence reservoir stress path, such as reservoir size and geometry or the coupled mechanical interaction between the reservoir and the bounding formations. Accordingly, a geomechanical model was developed to more accurately predict reservoir stress path. The model incorporates essential geological and geomechanical factors that may control reservoir stress path during production. In addition, laboratory results showing the effect of reservoir stress path on permeability and permeability anisotropy in a low-permeability sandstone are also presented. These experiments clearly demonstrate that during pore-pressure drawdown permeability decreases and that permeability parallel and perpendicular to the maximum stress direction decreases at different rates. The smallest reduction in permeability is parallel to the maximum principal stress. Consequently, stress-induced permeability anisotropy evolves with pore-pressure drawdown and the magnitude of permeability anisotropy increases at lower stress paths. Field Measurements of Stress Path in Lenticular Sandstone Reservoirs Salz10 presented hydraulic fracture stress data and pore-pressure measurements from reservoir pressure build-up tests in low-permeability, lenticular, gas sandstones of the Vicksburg formation in the McAllen Ranch field, Texas (Table 1). This work was one of the first studies to clearly show that the total minimum horizontal stress is dependent on the pore pressure. Hydraulic fractures were completed in underpressured and overpressured sandstone intervals from approximately 3100 to 3800 m. Some of the sandstones (9A, 10A, 11A, 12A, 13A, and 14A) were later hydraulically fractured a second time to improve oil productivity after several years of production. For initial reservoir conditions before production, the total minimum horizontal stress shows a decrease with decreasing pore pressure for different sandstone reservoirs. The effective stress can also be determined from these data. Following Rice and Cleary11 effective stress is defined by σ = S − α P , ( 1 ) where ? is the effective stress, S is the total stress, ? is a poroelastic parameter, and P is the pore pressure. For this study ? is assumed to equal unity. A linear regression analysis of the minimum horizontal and vertical effective stress data shows that at initial reservoir conditions the ratio of change in minimum effective horizontal stress to the change in effective vertical stress with increasing depth and pore pressure is 0.50.

1999 ◽  
Vol 2 (03) ◽  
pp. 266-272 ◽  
Author(s):  
H. Ruistuen ◽  
L.W. Teufel ◽  
D. Rhett

Summary The influence of production-induced changes in reservoir pore pressure on compressibility and permeability of weakly cemented sandstones has been analyzed. Laboratory experiments simulating reservoir depletion have been conducted over a range of stress paths that a reservoir may follow. The results suggest that compressibility of weakly cemented sandstones is stress path dependent. Compressibility measured under uniaxial strain conditions, or a stress path defined by a lower ratio of the rate at which the effective horizontal to effective vertical stress were increased than the one associated with uniaxial strain, is more than twice the corresponding value found from the hydrostatic loading experiment. In contrast, matrix permeability measured in the maximum stress direction show no significant stress path dependence. Experimental results suggest that a better understanding of the stress-sensitive behavior of weakly cemented sandstones can only be gained by dealing more directly with the microstructure of the rock. The stress-path-dependent nonlinear behavior of weakly cemented sandstones is related to effects of shear-enhanced compaction. Increasing cementation has been experimentally shown to reduce stress sensitivity. The observed nonlinearity is attributed to dilatancy rather than shear-enhanced compaction, also reflected by permeability measurements made in the maximum stress direction. Introduction Reliable data on rock compressibility and matrix permeability are essential in reservoir engineering due to the significant impact these parameters have on reserves and productivity estimations. Laboratory measurements of rock compressibility are applied to production forecasts, reservoir pressure maintenance evaluations, as well as reservoir compaction and subsidence studies,1–4 while matrix permeability heavily influences reservoir productivity and injectivity and is essential in performance forecasting.4 Formation compressibility is defined as the in situ bulk volume strain that results from changes in reservoir pore pressure: c = − 1 V i d V d P . ( 1 ) By adopting this definition, formation compressibility is not related to specific stress conditions. Formation compressibility is simply defined as the bulk response of the reservoir rock to production-induced changes in pore pressure. The stress changes that result from changes in pore pressure are uniquely defined by reservoir characteristics such as boundary conditions, reservoir geometry, and the mechanical properties of the reservoir rocks and bounding formations. A common procedure within the oil industry has been to use the so-called uniaxial correction factor to correct the results obtained from the hydrostatic compressibility test (cb) to "formation compressibility:"5 c = 1 + μ 3 ( 1 − μ ) c b . ( 2 ) An inherent assumption in this expression is that the rock is elastic throughout its production-induced deformation history, which may not be the case for weakly cemented reservoir rocks. The validity of the procedure also relies on the assumption that the uniaxial strain model adequately simulates reservoir conditions during depletion. Recent in situ stress measurements have demonstrated that this assumption is not necessarily valid. Since the early 1950's a number of researchers have investigated the relationships between rock matrix permeability and applied external pressure. Early observations suggested that permeability declines approximately exponentially with increasing confining pressure6 and that a relatively greater permeability reduction should be expected for a lower permeability matrix.7 These results were obtained from tests conducted under hydrostatic loading conditions. More recently, permeability measurements have also been performed under triaxial stress conditions.8–10 Matrix permeability has been related to compressibility and thus to the fabric and mineralogy of rocks. Bruno, Bovberg, and Nakagawa9 have shown that mineralogy may play a significant role in high-porosity rocks. Both increasing clay content and decreasing cementation resulted in a larger reduction in permeability with increasing stress. Holt8 reported experimental results on stress sensitivity of matrix permeability of a Jurassic sandstone. Samples were loaded under both triaxial compression and extension. No major differences in permeability were found between deviatoric and hydrostatic loading prior to yielding. At the yield stress, a sharp decline in permeability was observed. Most of the permeability reduction took place in the range of 60% of 90% of the peak shear stress. Teufel and Rhett3 introduced the term "stress path" to quantify the actual stress changes that take place in the reservoir during pressure depletion. (In this work, stress path is denoted K (not K0) to avoid confusion with uniaxial strain conditions, which is commonly denoted K0 test conditions. Also note that stress path here describes a constant ratio of change in stress state, which implies that different stress paths do not approach a common point in stress space.) The term describes the constant ratio of change in effective minimum (horizontal) stress to effective maximum (vertical) stress from initial reservoir conditions: K = Δ σ m i n Δ σ m a x . ( 3 ) The changes in the reservoir stress state resulting from depletion along stress paths of K=0, 0.5, and 1 are illustrated in Fig. 1. The importance of the reservoir stress path is that the shear stress has a larger increase for a lower stress path.


2021 ◽  
Author(s):  
Ahmed E. Radwan ◽  
Souvik Sen

Abstract The purpose of this study is to evaluate the reservoir geomechanics and stress path values of the depleted Miocene sandstone reservoirs of the Badri field, Gulf of Suez Basin, in order to understand the production-induced normal faulting potential in these depleted reservoirs. We interpreted the magnitudes of pore pressure (PP), vertical stress (Sv), and minimum horizontal stress (Shmin) of the syn-rift and post-rift sedimentary sequences encountered in the studied field, as well as we validated the geomechanical characteristics with subsurface measurements (i.e. leak-off test (LOT), and modular dynamic tests) (MDT). Stress path (ΔPP/ΔShmin) was modeled considering a pore pressure-horizontal stress coupling in an uniaxial compaction environment. Due to prolonged production, The Middle Miocene Hammam Faraun (HF) and Kareem reservoirs have been depleted by 950-1000 PSI and 1070-1200 PSI, respectively, with current 0.27-0.30 PSI/feet PP gradients as interpreted from initial and latest downhole measurements. Following the poroelastic approach, reduction in Shmin is assessed and reservoir stress paths values of 0.54 and 0.59 are inferred in the HF and Kareem sandstones, respectively. As a result, the current rate of depletion for both Miocene reservoirs indicates that reservoir conditions are stable in terms of production-induced normal faulting. Although future production years should be paid more attention. Accelerated depletion rate could have compelled the reservoirs stress path values to the critical level, resulting in depletion-induced reservoir instability. The operator could benefit from stress path analysis in future planning of infill well drilling and production rate optimization without causing reservoir damage or instability.


2021 ◽  
Author(s):  
Jianguo Zhang ◽  
Karthik Mahadev ◽  
Stephen Edwards ◽  
Alan Rodgerson

Abstract Maximum horizontal stress (SH) and stress path (change of SH and minimum horizontal stress with depletion) are the two most difficult parameters to define for an oilfield geomechanical model. Understanding these in-situ stresses is critical to the success of operations and development, especially when production is underway, and the reservoir depletion begins. This paper introduces a method to define them through the analysis of actual minifrac data. Field examples of applications on minifrac failure analysis and operational pressure prediction are also presented. It is commonly accepted that one of the best methods to determine the minimum horizontal stress (Sh) is the use of pressure fall-off analysis of a minifrac test. Unlike Sh, the magnitude of SH cannot be measured directly. Instead it is back calculated by using fracture initiation pressure (FIP) and Sh derived from minifrac data. After non-depleted Sh and SH are defined, their apparent Poisson's Ratios (APR) are calculated using the Eaton equation. These APRs define Sh and SH in virgin sand to encapsulate all other factors that influence in-situ stresses such as tectonic, thermal, osmotic and poro-elastic effects. These values can then be used to estimate stress path through interpretation of additional minifrac data derived from a depleted sand. A geomechanical model is developed based on APRs and stress paths to predict minifrac operation pressures. Three cases are included to show that the margin of error for FIP and fracture closure pressure (FCP) is less than 2%, fracture breakdown pressure (FBP) less than 4%. Two field cases in deep-water wells in the Gulf of Mexico show that the reduction of SH with depletion is lower than that for Sh.


2015 ◽  
Vol 2015 ◽  
pp. 1-7 ◽  
Author(s):  
Changqing Qi ◽  
Wei Lu ◽  
Jimin Wu ◽  
Xing Liu

Earthquake-induced liquefaction is one of the major causes of catastrophic earth dam failure. In order to assess the liquefaction potential and analyze the seismic performance of an earth dam in Fujian, Southeastern China, the in situ shear wave velocity test was firstly carried out. Results indicate that the gravelly filling is a type of liquefiable soil at present seismic setting. Then the effective stress model was adopted to thoroughly simulate the response of the soil to a proposed earthquake. Numerical result generally coincides with that of the empirical judgment based on in situ test. Negative excess pore pressure developed in the upper part of the saturated gravelly filling and positive excess pore pressure developed in the lower part. The excess pore pressure ratio increases with depth until it reaches a maximum value of 0.45. The displacement of the saturated gravelly soil is relatively small and tolerable. Results show that the saturated gravelly filling cannot reach a fully liquefied state. The dam is overall stable under the proposed earthquake.


1991 ◽  
Vol 28 (5) ◽  
pp. 650-659 ◽  
Author(s):  
Vinod K. Garga ◽  
Mahbubul A. Khan

Most of the laboratory testing methods available for the evaluation of in situ horizontal stresses are applicable to normally consolidated or lightly overconsolidated clays. This paper describes a new laboratory method for the determination of in situ horizontal stresses of heavily overconsolidated clays using a stress-path triaxial apparatus. The proposed method is based on the concept that if the radial stress exceeds the in situ horizontal stress, while maintaining the axial stress constant and equal to the in situ vertical effective stress, only then will the sample experience significant axial strain. The results obtained for undisturbed samples of an overconsolidated clay crust are found to be in agreement with some available methods. For verification of the applicability of the proposed method, K0 was determined for artificially prepared samples that had been subjected to known stress paths simulating field stress history. Key words: K0, overconsolidation, in situ stress, in situ test, clay crust, laboratory test.


Géotechnique ◽  
2022 ◽  
pp. 1-35
Author(s):  
S. L. Chen ◽  
Y. N. Abousleiman

A novel graphical analysis-based method is proposed for analysing the responses of a cylindrical cavity expanding under undrained conditions in modified Cam Clay soil. The essence of developing such an approach is to decompose and represent the strain increment/rate of a material point graphically into the elastic and plastic components in the deviatoric strain plane. It allows the effective stress path in the deviatoric plane to be readily determined by solving a first-order differential equation with the Lode angle being the single variable. The desired limiting cavity pressure and pore pressure can be equally conveniently evaluated, through basic numerical integrations with respect to the mean effective stress. Some ambiguity is clarified between the generalized (work conjugacy-based) shear strain increments and the corresponding deviatoric invariants of incremental strains. The present graph-based approach is also applicable for the determination of the stress and pore pressure distributions around the cavity. When used for predicting the ultimate cavity/pore pressures, it is computationally advantageous over the existing semi-analytical solutions that involve solving a system of coupled governing differential equations for the effective stress components. It thus may serve potentially as a useful and accurate interpretation of the results of in-situ pressuremeter tests on clay soils.


Author(s):  
Mojtaba P. Shahri ◽  
Stefan Z. Miska

There has been an increasing consciousness regarding stress changes associated with reservoir depletion as the industry moves towards more challenging jobs in deep-water or depleted reservoirs. These stress changes play a significant role in the design of wells in this condition. Therefore, accurate prediction of reservoir stress path, i.e., change in horizontal stresses with pore pressure, is of vital importance. In this study, the current stress path formulation is investigated using a Tri-axial Rock Mechanics Testing Facility. The reservoir depletion scenario is simulated through experiments and provides a better perspective on the currently used formulation and how it’s applicable during production and injection periods. The effect of fluid re-injection into reservoirs on the horizontal stress is also analyzed using core samples. According to the results, formation fracture pressure would not be equal to its initial value if pressure builds up using re-injection. The irrecoverable formation fracture pressure has a power law relation with pore pressure drawdown range. In order to avoid higher permanent fracture pressure reduction, it’s recommended to start the injection process as soon as possible during the production life of reservoirs. According to the experimental results, rocks behave differently during production and injection periods. Poisson’s ratio is greater during pressure build-up as compared to the depletion period. According to the current industry standards, Poisson’s ratio is usually obtained using fracturing data; i.e., leak-off test or mini-fracture test, or well logging methods. However, we are not able to use the same Poisson’s ratio for both pressure drawdown and build-up scenarios according to the experimental data. Corresponding to Poisson’s ratio values, the change in horizontal stress with pore pressure during drawdown (production) is higher than during build-up (injection) period. The outcomes of this study can significantly contribute to well planning and design of challenging wells over the life of reservoirs.


1997 ◽  
Vol 37 (1) ◽  
pp. 536
Author(s):  
R.R. Hillis ◽  
D.G. Crosby ◽  
A.K. Khurana

Theoretical fracture gradient relations are generally based on the assumption that the sedimentary sequence behaves elastically under conditions of lateral constraint. Hence the minimum horizontal stress (σhmin) is given by: where V is Poisson's ratio, σv is overburden stress, pp is pore pressure, and at is far -field tectonic stress. In driling practice, fracture initiation, or leak -off pressures, which are related to σhmin are most commonly predicted by the application of empirical stress /depth relations such as that proposed for offshore Western Australia by Vuckovic (1989): Leak -off pressure (psi) = 0.197D1145, where D is depth in feet. A modified form of the uniaxial elastic relation for the prediction of σhmin is proposed, such that: where the constants c and d are straight line regression constants derived from cross -plotting effective minimum horizontal stress and effective vertical stress. This relation, as opposed to previous empirical approaches to fracture gradient /σhmin determination, yields regression coefficients of physical significance: c represents the average Poisson's ratio term, v /(1 -v), and d represents an estimate of the tectonic (and inelastic) component of the minimum horizontal stress. This application of the modified fracture gradient relation, termed the effective stress cross -plot method, is tested successfully against published data from experimental wells in the East Texas Basin where independent estimates of Poisson's ratio are available. Leak -off pressures have been compiled from 61 wells in the Timor Sea. Leak -off pressures in the Timor Sea are somewhat lower than predicted by Vuckovic's (1989) stress /depth relation for offshore Western Australia, and a new, empirical stress /depth relation, which better fits the Timor Sea data is proposed: The effective stress cross -plot method is also applied to the Timor Sea data, yielding: Detailed pore pressure data were not available for the Timor Sea data -set and the effective stress cross -plot method does not fit the observed data any better than the new empirical stress /depth relation. However, the regression constants suggest an average Poisson's ratio of 0.26 and a relatively insignificant tectonic stress of 1 MPa for the Timor Sea.


Sign in / Sign up

Export Citation Format

Share Document