Extended Gassmann Equation with Dynamic Volumetric Strain: Modeling Wave Dispersion and Attenuation of Heterogenous Porous Rocks

Geophysics ◽  
2021 ◽  
pp. 1-97
Author(s):  
Luanxiao Zhao ◽  
Yirong Wang ◽  
Qiuliang Yao ◽  
Jianhua Geng ◽  
Hui Li ◽  
...  

Sedimentary rocks are often heterogeneous porous media inherently containing complex distributions of heterogeneities (e.g., fluid patches, cracks). Understanding and modeling their frequency-dependent elastic and adsorption behaviors is of great interest for subsurface rock characterization from multi-scale geophysical measurements. The physical parameter of dynamic volumetric strain (DVS) associated with wave-induced fluid flow is proposed to understand the common physics and connections behind known poroelastic models for modeling dispersion behaviors of heterogeneous rocks. We derive the theoretical formulations of DVS for patchy saturated rock at mesoscopic scale and cracked porous rock at microscopic grain scales, essentially embodying the wave-induced fluid pressure relaxation process. By incorporating the DVS into the classical Gassmann equation, a simple but practical “dynamic equivalent” modeling approach, extended Gassmann equation, is developed to characterize the dispersion and attenuation of complex heterogeneous rocks at non-zero frequencies. Using the extended Gassmann equation, the effect of microscopic or mesoscopic heterogeneities with complex distributions on the wave dispersion and attenuation signatures can be captured. The proposed theoretical framework provides a simple and straightforward analytical methodology to calculate wave dispersion and attenuation in porous rocks with multiple sets of heterogeneities exhibiting complex characteristics. We also demonstrate that, with the appropriate consideration of multiple crack sets and complex fluids patches distribution, the modeling results can better interpret the experimental data sets of dispersion and attenuation for heterogeneous porous rocks.

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA135-WA145 ◽  
Author(s):  
Fabian Krzikalla ◽  
Tobias M. Müller

Elastic upscaling of thinly layered rocks typically is performed using the established Backus averaging technique. Its poroelastic extension applies to thinly layered fluid-saturated porous rocks and enables the use of anisotropic effective medium models that are valid in the low- and high-frequency limits for relaxed and unrelaxed pore-fluid pressures, respectively. At intermediate frequencies, wave-induced interlayer flow causes attenuation and dispersion beyond that described by Biot’s global flow and microscopic squirt flow. Several models quantify frequency-dependent, normal-incidence P-wave propagation in layered poroelastic media but yield no prediction for arbitrary angles of incidence, or for S-wave-induced interlayer flow. It is shown that generalized models for P-SV-wave attenuation and dispersion as a result of interlayer flow can be constructed by unifying the anisotropic Backus limits with existing P-wave frequency-dependent interlayer flow models. The construction principle is exact and is based on the symmetry properties of the effective elastic relaxation tensor governing the pore-fluid pressure diffusion. These new theories quantify anisotropic P- and SV-wave attenuation and velocity dispersion. The maximum SV-wave attenuation is of the same order of magnitude as the maximum P-wave attenuation and occurs prominently around an angle of incidence of [Formula: see text]. For the particular case of a periodically layered medium, the theoretical predictions are confirmed through numerical simulations.


Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. WA141-WA151 ◽  
Author(s):  
Lin Zhang ◽  
Jing Ba ◽  
José M. Carcione ◽  
Weitao Sun

Understanding acoustic wave dispersion and attenuation induced by local (squirt) fluid flow between pores and cracks (compliant pores) is fundamental for better characterization of the porous rocks. To describe this phenomenon, some squirt-flow models have been developed based on the conservation of the fluid mass in the fluid mechanics. By assuming that the cracks are represented by isotropically distributed (i.e., randomly oriented) penny-shaped inclusions, this study applies the periodically oscillating squirt flow through inclusions based on the Biot-Rayleigh theory, so that the local squirt flow and global wave oscillation of rock are analyzed in the same theoretical framework of Hamilton’s principle. The governing wave-propagation equations are derived by incorporating all of the crack characteristics (such as the crack radius, crack density, and aspect ratio). In comparison with the previous squirt models, our model predicts the similar characteristics of wave velocity dispersion and attenuation, and our results are in agreement with Gassmann equations at the low-frequency limit. In addition, we find that the fluid viscosity and crack radius only affect the relaxation frequency of the squirt-flow attenuation peak, whereas the crack density and aspect ratio also affect the magnitudes of dispersion and attenuation. The application of this study to experimental data demonstrates that when the differential pressure (the difference between confining pressure and pore pressure) increases, the closure of cracks can lead to a decrease of attenuation. The results confirm that our model can be used to analyze and interpret the observed wave dispersion and attenuation of real rocks.


Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. D169-D181 ◽  
Author(s):  
Marco Milani ◽  
J. Germán Rubino ◽  
Tobias M. Müller ◽  
Beatriz Quintal ◽  
Eva Caspari ◽  
...  

Understanding and quantifying seismic energy dissipation in fluid-saturated porous rocks is of considerable interest because it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. An important, if not dominant, attenuation mechanism prevailing in the seismic frequency band is wave-induced fluid pressure diffusion in response to the contrasts in elastic stiffness in the mesoscopic-scale range. An effective way to estimate seismic velocity dispersion and attenuation related to this phenomenon is through the application of numerical upscaling procedures to synthetic rock samples of interest. However, the estimated seismic properties are meaningful only if the underlying sample volume is at least of the size of a representative elementary volume (REV). In the given context, the definition of an REV and the corresponding implications for the estimation of the effective seismic properties remain largely unexplored. To alleviate this problem, we have studied the characteristics of REVs for a set of idealized rock samples sharing high levels of velocity dispersion and attenuation. For periodically heterogeneous poroelastic media, the REV size was driven by boundary condition effects. Our results determined that boundary condition effects were absent for layered media and negligible in the presence of patchy saturation. Conversely, strong boundary condition effects arose in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. The results thus point to the importance of carefully determining the REV sizes of heterogeneous porous rocks for computing effective seismic properties, especially in the presence of strong dry frame stiffness contrasts.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. L45-L51 ◽  
Author(s):  
Boris Gurevich ◽  
Dina Makarynska

The Hashin-Shtrikman (HS) bounds define the range of bulk and shear moduli of an elastic composite, given the moduli of the constituents and their volume fractions. Recently, the HS bounds have been extended to the quasi-static moduli of composite viscoelastic media. Because viscoelastic moduli are complex, the viscoelastic bounds form a closed curve on the complex plane. We analyze these general viscoelastic bounds for a particular case of a porous solid saturated with a Newtonian fluid. In our analysis, for poroelastic media, the viscoelastic bounds for the bulk modulus are represented by a semicircle and a segment of the real axis, connecting formal HS bounds that are computed for an inviscid fluid. Importantly, viscoelastic bounds for poroelastic media turn out to be independent of frequency. However, because the bounds are quasi-static, the frequency must be much lower than Biot’s characteristic frequency. Furthermore, we find that the bounds for the bulk modulus are attainable (realizable). We also find that these viscoelastic bounds account for viscous shear relaxation and squirt-flow dispersion, but do not account for Biot’s global flow dispersion, because the latter strongly depends on inertial forces.


Sensors ◽  
2014 ◽  
Vol 14 (8) ◽  
pp. 15067-15083 ◽  
Author(s):  
Maria Strantza ◽  
Olivia Louis ◽  
Demosthenes Polyzos ◽  
Frans Boulpaep ◽  
Danny van Hemelrijck ◽  
...  

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. A51-A56 ◽  
Author(s):  
J. Germán Rubino ◽  
Gabriel A. Castromán ◽  
Tobias M. Müller ◽  
Leonardo B. Monachesi ◽  
Fabio I. Zyserman ◽  
...  

Numerical simulations of seismic wave propagation in fractured media are often performed in the framework of the linear slip theory (LST). Therein, fractures are represented as interfaces and their mechanical properties are characterized through a compliance matrix. This theory has been extended to account for energy dissipation due to viscous friction within fluid-filled fractures by using complex-valued frequency-dependent compliances. This is, however, not fully adequate for fractured porous rocks in which wave-induced fluid flow (WIFF) between fractures and host rock constitutes a predominant seismic attenuation mechanism. In this letter, we develop an approach to incorporate WIFF effects directly into the LST for a 1D system via a complex-valued, frequency-dependent fracture compliance. The methodology is validated for a medium permeated by regularly distributed planar fractures, for which an analytical expression for the complex-valued normal compliance is determined in the framework of quasistatic poroelasticity. There is good agreement between synthetic seismograms generated using the proposed recipe and those obtained from comprehensive, but computationally demanding, poroelastic simulations.


Sign in / Sign up

Export Citation Format

Share Document