In-situ stress calculation and application on the basis of anisotropic rock physics model

2016 ◽  
Author(s):  
Wang Changsheng* ◽  
Shi Yujiang ◽  
Wang Daxing ◽  
Zhang Haitao
Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. MR167-MR185 ◽  
Author(s):  
Romain Prioul ◽  
Richard Nolen-Hoeksema ◽  
MaryEllen Loan ◽  
Michael Herron ◽  
Ridvan Akkurt ◽  
...  

We have developed a method using measurements on drill cuttings as well as calibrated models to estimate anisotropic mechanical properties and stresses in unconventional reservoirs, when logs are not available in lateral wells. We measured mineralogy and organic matter on cuttings using diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). We described the methodology and illustrated it using two vertical control wells in the Vaca Muerta Formation, Argentina, and one lateral well drilled in the low-maturity oil-bearing reservoir. The method has two steps. First, using a vertical control well containing measurements from cuttings, a comprehensive logging suite, cores, and in situ stress tests, we define and calibrate four models: petrophysical, rock physics, dynamic-static elastic, and geomechanical. The petrophysical model provides petrophysical constituent volumes (mineralogy, organic matter, and fluids) from logs or DRIFTS inputs to the rock-physics model for calculating the dynamic anisotropic elastic moduli. The dynamic-static elastic and geomechanics models provide the relationships for computing static elastic properties and the minimum stress. Second, we acquire DRIFTS data on cuttings in the target lateral well and apply the four models for calculating stresses. We find that the method is successful for two reasons. First, the sonic-log-derived elastic moduli could be reconstructed accurately from the rock-physics model using input from petrophysical volumes from logs and DRIFTS data. A striking observation is that the elastic-property heterogeneity in those wells is explainable almost solely by compositional variations. Second, petrophysical volumes can be reconstructed by the petrophysical model and DRIFTS data. In the lateral well, we observed horizontal variations of mineralogy and organic matter, which controlled variations of elastic moduli and its anisotropy, and, in turn, affected partitioning of the gravitational and tectonic components in the minimum stress. This methodology promises accurate in situ stress estimates using cutting-based measurements and assessments of unconventional-reservoir heterogeneity.


2019 ◽  
Vol 38 (5) ◽  
pp. 358-365 ◽  
Author(s):  
Colin M. Sayers ◽  
Sagnik Dasgupta

This paper presents a predictive rock-physics model for unconventional shale reservoirs based on an extended Maxwell scheme. This model accounts for intrinsic anisotropy of rock matrix and heterogeneities and shape-induced anisotropy arising because the dimensions of kerogen inclusions and pores are larger parallel to the bedding plane than perpendicular to this plane. The model relates the results of seismic amplitude variation with offset inversion, such as P- and S-impedance, to the composition of the rock and enables identification of rock classes such as calcareous, argillaceous, siliceous, and mixed shales. This allows the choice of locations with the best potential for economic production of hydrocarbons. While this can be done using well data, prestack inversion of seismic P-wave data allows identification of the best locations before the wells are drilled. The results clearly show the ambiguity in rock classification obtained using poststack inversion of P-wave seismic data and demonstrate the need for prestack seismic inversion. The model provides estimates of formation anisotropy, as required for accurate determination of P- and S-impedance, and shows that anisotropy is a function not only of clay content but also other components of the rock as well as the aspect ratio of kerogen and pores. Estimates of minimum horizontal stress based on the model demonstrate the need to identify rock class and estimate anisotropy to determine the location of any stress barriers that may inhibit hydraulic fracture growth.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. D143-D150 ◽  
Author(s):  
Xinding Fang ◽  
Michael Fehler ◽  
Zhenya Zhu ◽  
Tianrun Chen ◽  
Stephen Brown ◽  
...  

Formation elastic properties near a borehole may be altered from their original state due to the stress concentration around the borehole. This could result in a biased estimation of formation properties but could provide a means to estimate in situ stress from sonic logging data. To properly account for the formation property alteration, we propose an iterative numerical approach to calculate stress-induced anisotropy around a borehole by combining a rock physics model and a finite-element method. We tested the validity and accuracy of our approach by comparing numerical results to laboratory measurements of the stress-strain relation of a sample of Berea sandstone, which contains a borehole and is subjected to a uniaxial stress loading. Our iterative approach converges very fast and can be applied to calculate the spatially varying stiffness tensor of the formation around a borehole for any given stress state.


2019 ◽  
Vol 17 (1) ◽  
pp. 70-85
Author(s):  
Ke-Ran Qian ◽  
Tao Liu ◽  
Jun-Zhou Liu ◽  
Xi-Wu Liu ◽  
Zhi-Liang He ◽  
...  

Abstract The brittleness prediction of shale formations is of interest to researchers nowadays. Conventional methods of brittleness prediction are usually based on isotropic models while shale is anisotropic. In order to obtain a better prediction of shale brittleness, our study firstly proposed a novel brittleness index equation based on the Voigt–Reuss–Hill average, which combines two classical isotropic methods. The proposed method introduces upper and lower brittleness bounds, which take the uncertainty of brittleness prediction into consideration. In addition, this method can give us acceptable predictions by using limited input values. Secondly, an anisotropic rock physics model was constructed. Two parameters were introduced into our model, which can be used to simulate the lamination of clay minerals and the dip angle of formation. In addition, rock physics templates have been built to analyze the sensitivity of brittleness parameters. Finally, the effects of kerogen, pore structure, clay lamination and shale formation dip have been investigated in terms of anisotropy. The prediction shows that the vertical/horizontal Young’s modulus is always below one while the vertical/horizontal Poisson’s ratio (PR) can be either greater or less than 1. Our study finds different degrees of shale lamination may be the explanation for the random distribution of Vani (the ratio of vertical PR to horizontal PR).


2017 ◽  
Vol 14 (4) ◽  
pp. 463-479 ◽  
Author(s):  
Ke-Ran Qian ◽  
Zhi-Liang He ◽  
Ye-Quan Chen ◽  
Xi-Wu Liu ◽  
Xiang-Yang Li

Geophysics ◽  
2021 ◽  
pp. 1-43
Author(s):  
Dario Grana

Rock physics models are physical equations that map petrophysical properties into geophysical variables, such as elastic properties and density. These equations are generally used in quantitative log and seismic interpretation to estimate the properties of interest from measured well logs and seismic data. Such models are generally calibrated using core samples and well log data and result in accurate predictions of the unknown properties. Because the input data are often affected by measurement errors, the model predictions are often uncertain. Instead of applying rock physics models to deterministic measurements, I propose to apply the models to the probability density function of the measurements. This approach has been previously adopted in literature using Gaussian distributions, but for petrophysical properties of porous rocks, such as volumetric fractions of solid and fluid components, the standard probabilistic formulation based on Gaussian assumptions is not applicable due to the bounded nature of the properties, the multimodality, and the non-symmetric behavior. The proposed approach is based on the Kumaraswamy probability density function for continuous random variables, which allows modeling double bounded non-symmetric distributions and is analytically tractable, unlike the Beta or Dirichtlet distributions. I present a probabilistic rock physics model applied to double bounded continuous random variables distributed according to a Kumaraswamy distribution and derive the analytical solution of the posterior distribution of the rock physics model predictions. The method is illustrated for three rock physics models: Raymer’s equation, Dvorkin’s stiff sand model, and Kuster-Toksoz inclusion model.


Sign in / Sign up

Export Citation Format

Share Document