A predictive anisotropic rock-physics model for estimating elastic rock properties of unconventional shale reservoirs

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.

2016 ◽  
Author(s):  
Wang Changsheng* ◽  
Shi Yujiang ◽  
Wang Daxing ◽  
Zhang Haitao

2006 ◽  
Author(s):  
Kyle Spikes ◽  
Jack Dvorkin ◽  
Gary Mavko

2020 ◽  
Vol 223 (1) ◽  
pp. 622-631
Author(s):  
Lin Zhang ◽  
Jing Ba ◽  
José M Carcione

SUMMARY Determining rock microstructure remains challenging, since a proper rock-physics model is needed to establish the relation between pore microstructure and elastic and transport properties. We present a model to estimate pore microstructure based on porosity, ultrasonic velocities and permeability, assuming that the microstructure consists on randomly oriented stiff equant pores and penny-shaped cracks. The stiff pore and crack porosity varying with differential pressure is estimated from the measured total porosity on the basis of a dual porosity model. The aspect ratio of pores and cracks and the crack density as a function of differential pressure are obtained from dry-rock P- and S-wave velocities, by using a differential effective medium model. These results are used to invert the pore radius from the matrix permeability by using a circular pore model. Above a crack density of 0.13, the crack radius can be estimated from permeability, and below that threshold, the radius is estimated from P-wave velocities, taking into account the wave dispersion induced by local fluid flow between pores and cracks. The approach is applied to experimental data for dry and saturated Fontainebleau sandstone and Chelmsford Granite.


2020 ◽  
Vol 117 (32) ◽  
pp. 18991-18997 ◽  
Author(s):  
Xin Gu ◽  
Gary Mavko ◽  
Lisa Ma ◽  
David Oakley ◽  
Natalie Accardo ◽  
...  

In weathered bedrock aquifers, groundwater is stored in pores and fractures that open as rocks are exhumed and minerals interact with meteoric fluids. Little is known about this storage because geochemical and geophysical observations are limited to pits, boreholes, or outcrops or to inferences based on indirect measurements between these sites. We trained a rock physics model to borehole observations in a well-constrained ridge and valley landscape and then interpreted spatial variations in seismic refraction velocities. We discovered that P-wave velocities track where a porosity-generating reaction initiates in shale in three boreholes across the landscape. Specifically, velocities of 2.7 ± 0.2 km/s correspond with growth of porosity from dissolution of chlorite, the most reactive of the abundant minerals in the shale. In addition, sonic velocities are consistent with the presence of gas bubbles beneath the water table under valley and ridge. We attribute this gas largely to CO2produced by 1) microbial respiration in soils as meteoric waters recharge into the subsurface and 2) the coupled carbonate dissolution and pyrite oxidation at depth in the rock. Bubbles may nucleate below the water table because waters depressurize as they flow from ridge to valley and because pores have dilated as the deep rock has been exhumed by erosion. Many of these observations are likely to also describe the weathering and flow path patterns in other headwater landscapes. Such combined geophysical and geochemical observations will help constrain models predicting flow, storage, and reaction of groundwater in bedrock systems.


2021 ◽  
Author(s):  
Vagif Suleymanov ◽  
Abdulhamid Almumtin ◽  
Guenther Glatz ◽  
Jack Dvorkin

Abstract Generated by the propagation of sound waves, seismic reflections are essentially the reflections at the interface between various subsurface formations. Traditionally, these reflections are interpreted in a qualitative way by mapping subsurface geology without quantifying the rock properties inside the strata, namely the porosity, mineralogy, and pore fluid. This study aims to conduct the needed quantitative interpretation by the means of rock physics to establish the relation between rock elastic and petrophysical properties for reservoir characterization. We conduct rock physics diagnostics to find a theoretical rock physics model relevant to the data by examining the wireline data from a clastic depositional environment associated with a tight gas sandstone in the Continental US. First, we conduct the rock physics diagnostics by using theoretical fluid substitution to establish the relevant rock physics models. Once these models are determined, we theoretically vary the thickness of the intervals, the pore fluid, as well as the porosity and mineralogy to generate geologically plausible pseudo-scenarios. Finally, Zoeppritz (1919) equations are exploited to obtain the expected amplitude versus offset (AVO) and the gradient versus intercept curves of these scenarios. The relationship between elastic and petrophysical properties was established using forward seismic modeling. Several theoretical rock physics models, namely Raymer-Dvorkin, soft-sand, stiff-sand, and constant-cement models were applied to the wireline data under examination. The modeling assumes that only two minerals are present: quartz and clay. The appropriate rock physics model appears to be constant-cement model with a high coordination number. The result is a seismic reflection catalogue that can serve as a field guide for interpreting real seismic reflections, as well as to determine the seismic visibility of the variations in the reservoir geometry, the pore fluid, and the porosity. The obtained reservoir properties may be extrapolated to prospects away from the well control to consider certain what-if scenarios like plausible lithology or fluid variations. This enables building of a catalogue of synthetic seismic reflections of rock properties to be used by the interpreter as a field guide relating seismic data to volumetric reservoir properties.


2018 ◽  
Vol 6 (4) ◽  
pp. SM1-SM8 ◽  
Author(s):  
Tingting Zhang ◽  
Yuefeng Sun

Fractured zones in deeply buried carbonate hills are important because they often have better permeability resulting in prolific production than similar low-porosity rocks. Nevertheless, their detection poses great challenge to conventional seismic inversion methods because they are mostly low in acoustic impedance and bulk modulus, hardly distinguishable from high-porosity zones or mudstones. A proxy parameter of pore structure defined in a rock-physics model, the so-called Sun model, has been used for delineating fractured zones in which the pore structure parameter is relatively high, whereas the porosity is low in general. Simultaneous seismic inversion of the pore structure parameter and porosity proves to be difficult and nontrivial in practice. Although the pore structure parameter is well-defined at locations where density, P-, and S-velocity are known from logs, estimation of P- and S-velocity information, especially density information from prestack seismic data is rather challenging. A three-step iterative inversion method, which uses acoustic, gradient, and elastic impedance from angle-stacked seismic data as input to the rock-physics model for calculating porosity and bulk and shear pore structure parameters simultaneously, is proposed and implemented to solve this problem. The methodology is successfully tested with well logs and seismic data from a deeply buried carbonate hill in the Bohai Bay Basin, China.


Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. MR1-MR13 ◽  
Author(s):  
Humberto S. Arévalo-López ◽  
Jack P. Dvorkin

Interpreting seismic data for petrophysical rock properties requires a rock-physics model that links the petrophysical rock properties to the elastic properties, such as velocity and impedance. Such a model can only be established from controlled experiments in which both groups of rock properties are measured on the same samples. A prolific source of such data is wellbore measurements. We use data from four wells drilled through a clastic offshore oil reservoir to perform rock-physics diagnostics, i.e., to find a theoretical rock-physics model that quantitatively explains the measurements. Using the model, we correct questionable well curves. Moreover, a crucial purpose of rock-physics diagnostics is to go beyond the settings represented in the wells and understand the seismic signatures of rock properties varying in a wider range via forward seismic modeling. With this goal in mind, we use our model to generate synthetic seismic gathers from perturbational modeling to address “what-if” scenarios not present in the wells.


2021 ◽  
Vol 9 ◽  
Author(s):  
Zhiqi Guo ◽  
Xueying Wang ◽  
Jian Jiao ◽  
Haifeng Chen

A rock physics model was established to calculate the P-wave velocity dispersion and attenuation caused by the squirt flow of fluids in gas hydrate-bearing sediments. The critical hydrate saturation parameter was introduced to describe different ways of hydrate concentration, including the mode of pore filling and the co-existence mode of pore filling and particle cementation. Rock physical modeling results indicate that the P-wave velocity is insensitive to the increase in gas hydrate saturation for the mode of pore filling, while it increases rapidly with increasing gas hydrate saturation for the co-existence mode of pore filling and particle cementation. Meanwhile, seismic modeling results show that both the PP and mode-converted PS reflections are insensitive to the gas hydrate saturation that is lower than the critical value, while they tend to change obviously for the hydrate saturation that is higher than the critical value. These can be interpreted that only when gas hydrate begins to be part of solid matrix at high gas hydrate saturation, it represents observable impact on elastic properties of the gas hydrate-bearing sediments. Synthetic seismograms are calculated for a 2D heterogeneous model where the gas hydrate saturation varies vertically and layer thickness of the gas hydrate-bearing sediment varies laterally. Modeling results show that larger thickness of the gas hydrate-bearing layer generally corresponds to stronger reflection amplitudes from the bottom simulating reflector.


Sign in / Sign up

Export Citation Format

Share Document