High‐resolution fracture characterization from azimuthal P‐wave seismic surveys using rock physics model‐based Bayesian inversion — Pinedale field, Wyoming

Author(s):  
Mita Sengupta ◽  
Ran Bachrach
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.


2018 ◽  
Vol 71 (0) ◽  
pp. 43-55 ◽  
Author(s):  
Yusuke Ohta ◽  
Tada-nori Goto ◽  
Katsuaki Koike ◽  
Koki Kashiwaya ◽  
Weiren Lin ◽  
...  

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.


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