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

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.

Geophysics ◽  
1979 ◽  
Vol 44 (11) ◽  
pp. 1789-1805 ◽  
Author(s):  
N. C. Dutta ◽  
H. Odé

In this investigation, Biot’s (1962) theory for wave propagation in porous solids is applied to study the velocity and attenuation of compressional seismic waves in partially gas‐saturated porous rocks. The Physical model, proposed by White (1975), is solved rigorously by using Biot’s equations which describe the coupled solid‐fluid motion of a porous medium in a systematic way. The quantitative results presented here are based on the theory described in Dutta and Odé (1979, this issue). We removed several of White’s questioned approximations and examined their effects on the quantitative results. We studied the variation of the attenuation coefficient with frequency, gas saturation, and size of gas inclusions in an otherwise brine‐filled rock. Anomalously large absorption (as large as 8 dB/cycle) at the exploration seismic frequency band is predicted by this model for young, unconsolidated sandstones. For a given size of the gas pockets and their spacing, the attenuation coefficient (in dB/cycle) increases almost linearly with frequency f to a maximum value and then decreases approximately as 1/√f. A sizable velocity dispersion (of the order of 30 percent) is also predicted by this model. A low gas saturation (4–6 percent) is found to yield high absorption and dispersion. An analysis of all of the field variables (stresses and displacements) is presented in terms of Biot’s type I (the classical compressional) wave and type II (the diffusion) wave. It is pointed out that the dissipation of energy in this model is mainly due to the relative fluid flow from the type II wave. From our formulation, many of White’s equations can be derived as suitable approximations, and it is shown that the discontinuity in fluid pressure assumed by White at the gas‐water interface is the discontinuity in the fluid pressure contribution by the type II wave. Our quantitative results are in reasonably good agreement with White’s (1975) approximate theory. However, the phase velocities computed by White’s approximate treatment do not approach the correct zerofrequency limit (Gassmann‐Wood) when compared to the present theory. Most of these disagreements disappear if the corrections to White’s theory as suggested by Dutta and Seriff (1979, this issue) are incorporated.


Geophysics ◽  
1972 ◽  
Vol 37 (3) ◽  
pp. 456-461 ◽  
Author(s):  
J. E. White ◽  
D. J. Walsh

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.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. N41-N45 ◽  
Author(s):  
Gracjan Lambert ◽  
Boris Gurevich ◽  
Miroslav Brajanovski

To explore the validity and limitations of the theoretical model of wave propagation in porous rocks with periodic distribution of planar fractures, we perform numerical simulation using a poroelastic reflectivity algorithm. The numerical results are found to be in good agreement with the analytical model, not only for periodic fractures, but also for random distribution of constant thickness fractures.


Geophysics ◽  
1979 ◽  
Vol 44 (11) ◽  
pp. 1777-1788 ◽  
Author(s):  
N. C. Dutta ◽  
H. Odé

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).


Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 309-318 ◽  
Author(s):  
Jorge O. Parra

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.


1988 ◽  
Vol 110 (4) ◽  
pp. 545-551 ◽  
Author(s):  
A. Cummings ◽  
I.-J. Chang

A quasi one-dimensional analysis of sound transmission in a flow duct lined with an array of nonlinear resonators is described. The solution to the equations describing the sound field and the hydrodynamic flow in the neighborhood of the resonator orifices is performed numerically in the time domain, with the object of properly accounting for the nonlinear interaction between the acoustic field and the resonators. Experimental data are compared to numerical computations in the time domain and generally very good agreement is noted. The method described here may readily be extended for use in the design of exhaust mufflers for internal combustion engines.


1982 ◽  
Vol 26 (04) ◽  
pp. 229-245 ◽  
Author(s):  
J. B. Roberts

By a combination of averaging techniques with the theory of Markov processes, an approximate theory is developed for the rolling motion of a ship in beam waves. A simple expression is obtained for the distribution of the roll angle, and is tested by a comparison with a set of digital simulation estimates due to Dalzell. Good agreement is obtained over a realistic range of damping values.


Author(s):  
M. Azeredo ◽  
◽  
V. Priimenko ◽  

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.


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 ◽  
...  

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