Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Wave propagation in heterogeneous, porous media: A velocity-stress, finite-difference method

Journal Article · · Geophysics
DOI:https://doi.org/10.1190/1.1443769· OSTI ID:37120
;  [1];  [2]
  1. Univ. of Alberta, Edmonton, Alberta (Canada)
  2. Technical Univ. of Crete, Chania (Greece)

A particle velocity-stress, finite-difference method is developed for the simulation of wave propagation in 2-D heterogeneous poroelastic media. Instead of the prevailing second-order differential equations, the authors consider a first-order hyperbolic system that is equivalent to Biot`s equations. The vector of unknowns in this system consists of the solid and fluid particle velocity components, the solid stress components, and the fluid pressure. A MacCormack finite-difference scheme that is fourth-order accurate in space and second-order accurate in time forms the basis of the numerical solutions for Biot`s hyperbolic system. An original analytic solution for a P-wave line source in a uniform poroelastic medium is derived for the purposes of source implementation and algorithm testing. In simulations with a two-layer model, additional slow compressional incident, transmitted, and reflected phases are recorded when the damping coefficient is small. This slow compressional wave is highly attenuated in porous media saturated by a viscous fluid. From the simulation they also verified that the attenuation mechanism introduced in Biot`s theory is of secondary importance for fast compressional and rotational waves. The existence of seismically observable differences caused by the presence of pores has been examined through synthetic experiments that indicate that amplitude variation with offset may be observed on receivers and could be diagnostic of the matrix and fluid parameters. This method was applied in simulating seismic wave propagation over an expanded steam-heated zone in Cold Lake, Alberta in an area of enhanced oil recovery (EOR) processing. The results indicate that a seismic surface survey can be used to monitor thermal fronts.

OSTI ID:
37120
Journal Information:
Geophysics, Journal Name: Geophysics Journal Issue: 2 Vol. 60; ISSN GPYSA7; ISSN 0016-8033
Country of Publication:
United States
Language:
English

Similar Records

Poroelastic model to relate seismic wave attenuation and dispersion to permeability anisotropy
Journal Article · Mon Jan 31 23:00:00 EST 2000 · Geophysics · OSTI ID:20015462

Wave simulation in anisotropic, saturated porous media
Conference · Sat Dec 30 23:00:00 EST 1995 · OSTI ID:542952

Full frequency-range transient solution for compressional waves in a fluid-saturated viscoacoustic porous medium
Journal Article · Sun Dec 31 23:00:00 EST 1995 · Geophysical Prospecting · OSTI ID:207851