scholarly journals Numerical upscaling in 2-D heterogeneous poroelastic rocks: Anisotropic attenuation and dispersion of seismic waves

2016 ◽  
Vol 121 (9) ◽  
pp. 6698-6721 ◽  
Author(s):  
J. Germán Rubino ◽  
Eva Caspari ◽  
Tobias M. Müller ◽  
Marco Milani ◽  
Nicolás D. Barbosa ◽  
...  
Keyword(s):  
Seismic Waves ◽  
Attenuation And Dispersion

The transversely isotropic poroelastic wave equation including the Biot and the squirt mechanisms: Theory and application

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 309-318 ◽  
Author(s):  
Jorge O. Parra
Keyword(s):  
Wave Equation ◽  
Fluid Flow ◽  
Analytical Solution ◽  
Seismic Waves ◽  
Transversely Isotropic ◽  
P Wave ◽  
Fluid Properties ◽  
Permeability Anisotropy ◽  
Attenuation And Dispersion ◽  
Maximum Permeability

The transversely isotropic poroelastic wave equation can be formulated to include the Biot and the squirt‐flow mechanisms to yield a new analytical solution in terms of the elements of the squirt‐flow tensor. The new model gives estimates of the vertical and the horizontal permeabilities, as well as other measurable rock and fluid properties. In particular, the model estimates phase velocity and attenuation of waves traveling at different angles of incidence with respect to the principal axis of anisotropy. The attenuation and dispersion of the fast quasi P‐wave and the quasi SV‐wave are related to the vertical and the horizontal permeabilities. Modeling suggests that the attenuation of both the quasi P‐wave and quasi SV‐wave depend on the direction of permeability. For frequencies from 500 to 4500 Hz, the quasi P‐wave attenuation will be of maximum permeability. To test the theory, interwell seismic waveforms, well logs, and hydraulic conductivity measurements (recorded in the fluvial Gypsy sandstone reservoir, Oklahoma) provide the material and fluid property parameters. For example, the analysis of petrophysical data suggests that the vertical permeability (1 md) is affected by the presence of mudstone and siltstone bodies, which are barriers to vertical fluid movement, and the horizontal permeability (1640 md) is controlled by cross‐bedded and planar‐laminated sandstones. The theoretical dispersion curves based on measurable rock and fluid properties, and the phase velocity curve obtained from seismic signatures, give the ingredients to evaluate the model. Theoretical predictions show the influence of the permeability anisotropy on the dispersion of seismic waves. These dispersion values derived from interwell seismic signatures are consistent with the theoretical model and with the direction of propagation of the seismic waves that travel parallel to the maximum permeability. This analysis with the new analytical solution is the first step toward a quantitative evaluation of the preferential directions of fluid flow in reservoir formation containing hydrocarbons. The results of the present work may lead to the development of algorithms to extract the permeability anisotropy from attenuation and dispersion data (derived from sonic logs and crosswell seismics) to map the fluid flow distribution in a reservoir.

Numerical simulation of seismic waves in porous media components

2021 ◽  
Vol 9 (1) ◽  
pp. 4-30
Author(s):  
M. Azeredo ◽  
◽  
V. Priimenko ◽  
Keyword(s):  
Numerical Simulation ◽  
Porous Medium ◽  
High Frequency ◽  
Seismic Waves ◽  
Computational Algorithm ◽  
Low Frequency ◽  
Physical Parameters ◽  
Mathematical Algorithm ◽  
Biot Equations ◽  
Attenuation And Dispersion

This work presents a mathematical algorithm for modeling the propagation of poroelastic waves. We have shown how the classical Biot equations can be put into Ursin’s form in a plane-layered 3D porous medium. Using this form, we have derived explicit for- mulas that can be used as the basis of an efficient computational algorithm. To validate the algorithm, numerical simulations were performed using both the poroelastic and equivalent elastic models. The results obtained confirmed the proposed algorithm’s reliability, identify- ing the main wave events in both low-frequency and high-frequency regimes in the reservoir and laboratory scales, respectively. We have also illustrated the influence of some physical parameters on the attenuation and dispersion of the slow wave.

Seismic wave attenuation due to fluid pressure diffusion at the mesoscopic scale: an experimental and numerical study

2021 ◽  
Author(s):  
Samuel Chapman ◽  
Jan V. M. Borgomano ◽  
Beatriz Quintal ◽  
Sally M. Benson ◽  
Jerome Fortin
Keyword(s):  
Seismic Wave ◽  
Seismic Waves ◽  
Wave Attenuation ◽  
Fluid Pressure ◽  
Fluid Distribution ◽  
Mesoscopic Scale ◽  
Seismic Wave Attenuation ◽  
X Ray ◽  
Pressure Diffusion ◽  
Attenuation And Dispersion

<p>Monitoring of the subsurface with seismic methods can be improved by better understanding the attenuation of seismic waves due to fluid pressure diffusion (FPD). In porous rocks saturated with multiple fluid phases the attenuation of seismic waves by FPD is sensitive to the mesoscopic scale distribution of the respective fluids. The relationship between fluid distribution and seismic wave attenuation could be used, for example, to assess the effectiveness of residual trapping of carbon dioxide (CO2) in the subsurface. Determining such relationships requires validating models of FPD with accurate laboratory measurements of seismic wave attenuation and modulus dispersion over a broad frequency range, and, in addition, characterising the fluid distribution during experiments. To address this challenge, experiments were performed on a Berea sandstone sample in which the exsolution of CO2 from water in the pore space of the sample was induced by a reduction in pore pressure. The fluid distribution was determined with X-ray computed tomography (CT) in a first set of experiments. The CO2 exosolved predominantly near the outlet, resulting in a heterogeneous fluid distribution along the sample length. In a second set of experiments, at similar pressure and temperature conditions, the forced oscillation method was used to measure the attenuation and modulus dispersion in the partially saturated sample over a broad frequency range (0.1 - 1000 Hz). Significant P-wave attenuation and dispersion was observed, while S-wave attenuation and dispersion were negligible. These observations suggest that the dominant mechanism of attenuation and dispersion was FPD. The attenuation and dispersion by FPD was subsequently modelled by solving Biot’s quasi-static equations of poroelasticity with the finite element method. The fluid saturation distribution determined from the X-ray CT was used in combination with a Reuss average to define a single phase effective fluid bulk modulus. The numerical solutions agree well with the attenuation and modulus dispersion measured in the laboratory, supporting the interpretation that attenuation and dispersion was due to FPD occurring in the heterogenous distribution of the coexisting fluids. The numerical simulations have the advantage that the models can easily be improved by including sub-core scale porosity and permeability distributions, which can also be determined using X-ray CT. In the future this could allow for conducting experiments on heterogenous samples.</p>

PROPOSED ATTENUATION‐DISPERSION PAIR FOR SEISMIC WAVES

Geophysics ◽  
1972 ◽  
Vol 37 (3) ◽  
pp. 456-461 ◽  
Author(s):  
J. E. White ◽  
D. J. Walsh
Keyword(s):  
Experimental Data ◽  
Seismic Waves ◽  
Numerical Application ◽  
Low Frequencies ◽  
One Dimensional ◽  
Lumped Element ◽  
Resistive Element ◽  
The Hilbert Transform ◽  
Attenuation And Dispersion ◽  
Mathematical Behavior

Several papers in recent years have dealt with the causality‐imposed relation between attenuation and dispersion for waves in lossy solids, with emphasis on seismic waves. While the published formulas for dispersion within a particular frequency band are supported by experimental evidence within that band, the mathematical behavior of these expressions outside the band, particularly at low frequencies, is physically unacceptable. In the present paper, one‐dimensional seismic waves are modeled as propagation along a simple lumped‐element transmission line, leading to expressions for attenuation and velocity as functions of frequency which not only satisfy the experimental data available, but exhibit no objectionable behavior outside the range of available data. This is achieved by introducing a resistive element whose value is inversely proportional to frequency. Numerical application of the Hilbert transform shows the condition of causality to be satisfied by this model quite accurately.

Attenuation and dispersion of predicted seismic waves in the simplified poroelastic theory

2018 ◽  
Author(s):  
Haixia Zhao ◽  
Xiaokai Wang ◽  
Bangyu Wu ◽  
Jingrui Luo ◽  
Hui Li
Keyword(s):  
Seismic Waves ◽  
Poroelastic Theory ◽  
Attenuation And Dispersion

On White’s model of attenuation in rocks with partial gas saturation

Geophysics ◽  
1979 ◽  
Vol 44 (11) ◽  
pp. 1806-1812 ◽  
Author(s):  
N. C. Dutta ◽  
A. J. Seriff
Keyword(s):  
Seismic Waves ◽  
Numerical Calculations ◽  
Porous Solids ◽  
Approximate Theory ◽  
Porous Rocks ◽  
Low Frequencies ◽  
Exact Calculations ◽  
Attenuation And Dispersion ◽  
Gas Compressibility ◽  
Good Agreement

In two important papers, J. E. White and coauthors (White, 1975; White et al, 1976) have given an approximate theory for the calculation of attenuation and dispersion of compressional seismic waves in porous rocks filled mostly with brine but containing gas‐filled regions. Modifications of White’s formulas for [Formula: see text] and Q in the case of gas‐filled spheres brings the results into good agreement with the more exact calculations of Dutta and Odé (1979a, b, this issue), who used Biot’s theory for porous solids. In particular, the modified formulas give the expected Gassmann‐Wood velocity at very low frequencies. Inclusion of the finite gas compressibility in numerical calculations for gas‐filled spheres shows an interesting maximum of the attenuation at low gas saturations which is not seen if the gas is ignored. A comparison of the attenuation calculated for the same rock and fluids but for three different geometries of the gas‐filled regions suggests that the configuration of the gas‐filled zones does not have an important effect on the magnitude of the attenuation.

Modeling Attenuation of Diffusive-Viscous Wave Using Reflectivity Method

2019 ◽  
Vol 27 (03) ◽  
pp. 1850030
Author(s):  
Haixia Zhao ◽  
Jinghuai Gao ◽  
Jigen Peng ◽  
Gulan Zhang
Keyword(s):  
Seismic Waves ◽  
Numerical Models ◽  
Pore Space ◽  
Low Frequency ◽  
Compressional Wave ◽  
Viscous Medium ◽  
Intrinsic Attenuation ◽  
Reflectivity Method ◽  
Strong Amplitude ◽  
Attenuation And Dispersion

Seismic waves in earth materials are subject to attenuation and dispersion in a broad range of frequencies. The commonly accepted mechanism of intrinsic attenuation and dispersion is the presence of fluids in the pore space of rocks. The diffusive-viscous model was proposed to explain low-frequency seismic anomalies related to hydrocarbon reservoirs. But, the model is only a description of compressional wave. In this work, we firstly discuss the extended elastic diffusive-viscous model. Then, we extend reflectivity method to the diffusive-viscous medium. Finally, we present two numerical models to simulate the attenuation of diffusive-viscous wave in horizontal and dip multi-layered media compared with the results of viscoelastic wave. The modeling results show that the diffusive-viscous wave has strong amplitude attenuation and phase shift when it propagates across absorptive layers.

Attenuation and dispersion of compressional waves in fluid‐filled porous rocks with partial gas saturation (White model)—Part I: Biot theory

Geophysics ◽  
1979 ◽  
Vol 44 (11) ◽  
pp. 1777-1788 ◽  
Author(s):  
N. C. Dutta ◽  
H. Odé
Keyword(s):  
Seismic Waves ◽  
Equations Of Motion ◽  
Fluid Pressure ◽  
Present Theory ◽  
Biot Theory ◽  
Type I ◽  
Porous Rocks ◽  
Water Contact ◽  
Coupled Equations ◽  
Attenuation And Dispersion

An exact theory of attenuation and dispersion of seismic waves in porous rocks containing spherical gas pockets (White model) is presented using the coupled equations of motion given by Biot. Assumptions made are (1) the acoustic wavelength is long with respect to the distance between gas pockets and their size, and (2) the gas pockets do not interact. Thus, the present theory essentially is quite similar to that proposed by White (1975), but the problem of the radially oscillating gas pocket is solved in a more rigorous manner by means of Biot’s theory (1962). The solid‐fluid coupling is automatically included, and the model is solved as a boundary value problem requiring all radial stresses and displacements to be continuous at the gas‐brine interface. Thus, we do not require any assumed fluid‐pressure discontinuity at the gas‐water contact, such as the one employed by White (1975). We have also presented an analysis of all of the field variables in terms of Biot’s type I (the classical compressional) wave and, type II (the diffusion) wave. Our quantitative results are presented in Dutta and Odé (1979, this issue).

Sign in / Sign up

Export Citation Format

Share Document