Computation of the spatial derivatives with non-local differential operators, such as the Fourier pseudospectral method, may cause strong numerical artifacts in the form of noncausal ringing. This situation occurs when regular grids are used. The problem is attacked by using a staggered pseudospectral technique. The nature and the causes of acausal ringing in regular grid methods and the reasons why staggered grid methods eliminate this problem are explained in the papers by Fornberg [7] and Öz envar and McMechan [8]. Thus, the objective here is not to propose a new method, but to develop the algorithm for the poroviscoelastic wave equation [3]. Özdenvar and McMechan [9] developed a pseudospectral staggered-grid algorithm for the poroelastic differential equations expressed in the displacement formulation. Here, the equations are expressed in the velocity–stress formulation, including attenuation mechanisms due to fluid/matrix interactions. As is well known, one of the advantages of using the velocity–stress differential equations is the avoidance of the differentiation of the material properties, which eliminates numerical artifacts and makes the algorithm more stable [11]. The algorithm is illustrated for a physical situation that requires very high accuracy, such as the gas–water contact in a natural gas reservoir, where large contrasts in material properties occur. The staggered-grid solution is noise-free in the dynamic range where regular grids generate artifacts that may have amplitudes similar to those of physical arrivals.
[1]
George A. McMechan,et al.
Algorithms for staggered‐grid computations for poroelastic, elastic, acoustic, and scalar wave equations
,
1997
.
[2]
B. Fornberg.
High-order finite differences and the pseudospectral method on staggered grids
,
1990
.
[3]
N. Dutta,et al.
Seismic reflections from a gas‐water contact
,
1983
.
[4]
José M. Carcione,et al.
SOME ASPECTS OF THE PHYSICS AND NUMERICAL MODELING OF BIOT COMPRESSIONAL WAVES
,
1995
.
[5]
G. McMechan,et al.
Causes and reduction of numerical artefacts in pseudo-spectral wavefield extrapolation
,
1996
.
[6]
José M. Carcione,et al.
Wave propagation in anisotropic, saturated porous media: Plane‐wave theory and numerical simulation
,
1996
.
[7]
J. Carcione.
Viscoelastic effective rheologies for modelling wave propagation in porous media
,
1998
.
[8]
J. Virieux.
P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method
,
1986
.
[9]
M. Biot,et al.
THE ELASTIC COEFFICIENTS OF THE THEORY OF CONSOLIDATION
,
1957
.
[10]
M. Biot.
MECHANICS OF DEFORMATION AND ACOUSTIC PROPAGATION IN POROUS MEDIA
,
1962
.