Seismic pore space compressibility and Gassmann’s relation

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1743-1749 ◽  
Author(s):  
Gary Mavko ◽  
Tapan Mukerji

The pore space compressibility of a rock provides a robust, model‐independent descriptor of porosity and pore fluid effects on effective moduli. The pore space compressibility is also the direct physical link between the dry and fluid‐saturated moduli, and is therefore the basis of Gassmann’s equation for fluid substitution. For a fixed porosity, an increase in pore space compressibility increases the sensitivity of the modulus to fluid substitution. Two simple techniques, based on pore compressibility, are presented for graphically applying Gassmann’s relation for fluid substitution. In the first method, the pore compressibility is simply reweighted with a factor that depends only on the ratio of fluid to mineral bulk modulus. In the second technique, the rock moduli are rescaled using the Reuss average, which again depends only on the fluid and mineral moduli.

Geophysics ◽  
2003 ◽  
Vol 68 (2) ◽  
pp. 430-440 ◽  
Author(s):  
Tad M. Smith ◽  
Carl H. Sondergeld ◽  
Chandra S. Rai

Fluid substitution is an important part of seismic attribute work, because it provides the interpreter with a tool for modeling and quantifying the various fluid scenarios which might give rise to an observed amplitude variation with offset (AVO) or 4D response. The most commonly used technique for doing this involves the application of Gassmann's equations. Modeling the changes from one fluid type to another requires that the effects of the starting fluid first be removed prior to modeling the new fluid. In practice, the rock is drained of its initial pore fluid, and the moduli (bulk and shear) and bulk density of the porous frame are calculated. Once the porous frame properties are properly determined, the rock is saturated with the new pore fluid, and the new effective bulk modulus and density are calculated. A direct result of Gassmann's equations is that the shear modulus for an isotropic material is independent of pore fluid, and therefore remains constant during the fluid substitution process. In the case of disconnected or cracklike pores, however, this assumption may be violated. Once the values for the new effective bulk modulus and bulk density are calculated, it is possible to calculate the compressional and shear velocities for the new fluid conditions. There are other approaches to fluid substitution (empirical and heuristic) which avoid the porous frame calculations but, as described in this tutorial, often do not yield reliable results. This tutorial provides the reader with a recipe for performing fluid substitutions, as well as insight into why and when the approach may fail.


Geophysics ◽  
2013 ◽  
Vol 78 (5) ◽  
pp. L87-L99 ◽  
Author(s):  
Gary Mavko ◽  
Nishank Saxena

Fluid and solid substitution of bulk modulus are exact and unique for materials whose elastic bulk and/or shear moduli fall on the Hashin-Shtrikman bounds. For materials whose moduli lie between the bounds, solid and fluid substitution of bulk moduli can be computed exactly, but not uniquely. Every initial bulk modulus can be realized with an infinite number of microstructures and therefore transform to an infinite number of moduli upon substitution of the pore fill. This nonuniqueness arises when detailed information on the material pore geometry is not available. We evaluated four embedded-bound constructions for fluid and solid substitution that were based on realizable materials. In the limiting case of pore fluids, two of these constructions reduced to the bounds of Gibiansky and Torquato, which illustrated that those bounds were optimum. For solids, the first two constructions corresponded to a homogeneous pore stiffness and predicted the smallest change in modulus. The third construction prediction corresponded to a pore space with heterogeneous stiffness, and it predicted a much larger change in modulus.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. L21-L32 ◽  
Author(s):  
Nishank Saxena ◽  
Gary Mavko

We derived exact equations, elastic bulk and shear, for fluid and solid substitution in monomineralic isotropic rocks of arbitrary pore shape and suggested methods to obtain the required substitution parameters. We proved that the classical Gassmann’s bulk modulus equation for fluid-to-fluid substitution is exact for solid-to-solid substitution if compression-induced mean stresses (pressure) in initial and final pore solids are homogeneous and either the shear modulus of the substituted solid does not change or no shear stress is induced in pores. Moreover, when compression-induced mean stresses in initial and final pore solids are homogeneous, we evaluated exact generalizations of Gassmann’s bulk modulus equation, which depend on usually known parameters. For the effective shear modulus, we found general exactness conditions of Gassmann and other approximations. Using the new exact substitution equations, we interpreted that predicting solid-filled rock stiffness from a dry rock stiffness measurement requires more information (i.e., assumptions about the pore shape) compared to predicting the same from a fluid-saturated rock stiffness.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. MR201-MR212
Author(s):  
Zhi-Qiang Yang ◽  
Tao He ◽  
Chang-Chun Zou

Velocity dispersion is a common phenomenon for fluid-charged porous rocks and carries important information on the pore structure and fluid in reservoir rocks. Previous ultrasonic experiments had measured more significant non-Biot velocity dispersion on saturated reservoir sandstones with increasing pore-fluid viscosity. Although wave-induced local squirt-flow effect could in theory cause most of the non-Biot velocity dispersion, its quantitative prediction remains a challenge. Several popular models were tested to predict the measured velocities under undrained conditions, but they either underestimated the squirt-flow effect or failed to simultaneously satisfy P- and S-wave velocity dispersions (especially for higher viscosity fluids). Based on the classic double-porosity theory that pore space is comprised of mainly stiff/Biot’s porosity and minor compliant porosity, an effective “wet frame” was hypothesized to account for the squirt-flow effect, whose compliant pores are filled with a hypothesized fluid with dynamic modulus. A new dynamic elastic model was then introduced by extending Biot theory to include the squirt-flow effect, after replacing the dry-frame bulk/shear moduli with their wet-frame counterparts. In addition to yielding better velocity predictions for P- and S-wave measurements of different fluid viscosities, the new model is also more applicable because its two key tuning parameters (i.e., the effective aspect ratio and porosity of compliant pores) at in situ reservoir pressure could be constrained with laboratory velocity measurements associated with pore-fluid viscosities.


2019 ◽  
Vol 89 (12) ◽  
pp. 1231-1249 ◽  
Author(s):  
Kristin W. Porten ◽  
Michał J. Warchoł ◽  
Ian A. Kane

ABSTRACT Well-developed detrital clay grain coats are observed in deep-marine sandstones of the Upper Cretaceous Springar Formation of the Vøring Basin in the Norwegian Sea. The detrital clay coats form thin and compact rims on individual sand grains and meniscus-shaped bridges between grains. These well-developed coats are found in high-density turbidites and proximal hybrid event beds with common to pervasive dewatering structures deposited in proximity to the base of a syndepositionally active basin high. Here, in one exploration well, detrital clay grain coats are common throughout a sandstone package 100 m thick. High-density turbidites and proximal and distal hybrid event beds drilled in mid- to distal-fan settings unaffected by seismically resolved seafloor topography show common dewatering features, but have only scattered detrital clay coats confined to individual dewatering pipes or dish structures. Hence, we propose that intense sediment dewatering has the potential to form detrital clay coats in deep-marine sandstones by a combination of elutriation and reorganization of clays during fluid escape from sediment bodies with pore fluid pressures significantly higher than the hydrostatic pressure. In submarine fan systems, deposition of sediment with coeval trapping of large volumes of interstitial pore fluid is most likely to occur where gravity flows undergo rapid deceleration in response to an abrupt decrease in confinement or gradient. Such environments include the channel–lobe transition and settings in proximity to seabed topography. The investigated sandstones are quartz arenites and subarkoses, with minor to moderate volumes of quartz cement (up to 6%). However, strongly to completely quartz-cemented intergranular pore space is observed where detrital clay coats or matrix does not cover quartz grains in the deepest part of the studied formation. Modeling of quartz cementation predicts that most intergranular macroporosity in the lower part of the Springar Formation would be quartz cemented if the sandstones were free of detrital clays. Based on our observations and modeling results we propose that intense sediment dewatering has the potential to form detrital clay coats, which can be important for retaining porosity in deeply buried sandstones and in basins with high present or past heat flow.


Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1222-1236 ◽  
Author(s):  
Nabil Akbar ◽  
Gary Mavko ◽  
Amos Nur ◽  
Jack Dvorkin

We investigate the effects of permeability, frequency, and fluid distribution on the viscoelastic behavior of rock. The viscoelastic response of rock to seismic waves depends on the relative motion of pore fluid with respect to the solid phase. Fluid motion depends, in part, on the internal wave‐induced pore pressure distribution that relates to the pore micro‐structure of rock and the scales of saturation. We consider wave‐induced squirt fluid flow at two scales: (1) local microscopic flow at the smallest scale of saturation heterogeneity (e.g., within a single pore) and (2) macroscopic flow at a larger scale of fluid‐saturated and dry patches. We explore the circumstances under which each of these mechanisms prevails. We examine such flows under the conditions of uniform confining (bulk) compression and obtain the effective dynamic bulk modulus of rock. The solutions are formulated in terms of generalized frequencies that depend on frequency, saturation, fluid and gas properties, and on the macroscopic properties of rock such as permeability, porosity, and dry bulk modulus. The study includes the whole range of saturation and frequency; therefore, we provide the missing link between the low‐frequency limit (Gassmann’s formula) and the high‐frequency limit given by Mavko and Jizba. Further, we compare our model with Biot’s theory and introduce a geometrical factor whose numeric value gives an indication as to whether local fluid squirt or global (squirt and/or Biot’s) mechanisms dominate the viscoelastic properties of porous materials. The important results of our theoretical modeling are: (1) a hysteresis of acoustic velocity versus saturation resulting from variations in fluid distributions, and (2) two peaks of acoustic wave attenuation—one at low frequency (caused by global squirt‐flow) and another at higher frequency (caused by local flow). Both theoretical results are compared with experimental data.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. M19-M28 ◽  
Author(s):  
Gary Mavko

The interaction of pore stiffness with pore fluid moduli leads to shifts in viscoelastic relaxation times of the overall rock relative to those of the fluids alone. Crack-based and fluid substitution models indicate that stiff pores cause little shift, whereas thin, soft cracks can shift relaxation times by several orders of magnitude toward lower frequencies (longer relaxation times). Pore stiffness also causes a shift in apparent temperature dependence of rock viscoelasticity toward higher temperatures when cracks are present. As with more conventional fluid substitution problems, quantifying the effects of pore fluids on rock properties requires information about the crack and pore stiffness distributions in addition to the complex moduli and viscosity of the pure fluid.


The dynamic response of isotropic composites of two viscoelastic isotropic phases mixed in fixed proportions is considered in the frequency range where the acoustic wavelength is much larger than the inhomogeneities. The effective bulk-modulus bounds of Hashin-Shtrikman-Walpole are extended to viscoelasticity in this quasi-static régime, where the properties of the isotropic composite can be described by complex bulk and shear moduli. The effective bulk modulus is shown to be constrained to a lens-shaped region of the complex plane bounded by the outermost pair of four circular arcs (three circular arcs in the case of two-dimensional elasticity). This is proved using a new variational principle for viscoelasticity together with two established techniques for deriving bounds on effective moduli, namely the translation method and the Hashin-Shtrikman method. In this application the Hashin-Shtrikman method needs to be generalized to allow the reference tensor to have an associated quasiconvex energy. Microstructures are identified which have bulk-moduli that correspond to various points on each of the circular arcs. Thus these microstructures have extremal viscoelastic behaviour when the associated arc forms one of the outermost pair. The bounds and the extremal microstructures are similar to those obtained for the complex dielectric constant, but the methods used here are entirely different.


Sign in / Sign up

Export Citation Format

Share Document