Application of a spectral method to numerical modeling of the propagation of seismic waves in porous media with energy dissipation

Abstract

This paper presents an algorithm based on the spectral Laguerre method for approximation of time derivatives as applied to a problem of seismic wave propagation in porous media with energy dissipation. The initial system of equations is written as a first-order hyperbolic system in terms of velocities, stresses, and pore pressure. To numerically solve the problem, a combination of an analytical Laguerre transform and a finite-difference method is used. The method proposed in the paper is an analog of a well-known spectral method based on the Fourier transform. However, unlike the Fourier transform, the integral Laguerre transform with respect to time reduces the initial problem to a system of equations in which the expansion parameter is present only in the right-hand side of the equations as a recurrence relation. As compared to the finite-difference method, with an analytical transform in the spectral method it is possible to reduce the original problem to a system of differential equations having only derivatives with respect to the spatial coordinates. This allows using the known stable difference scheme for the recurrence solutions to similar systems. Such an approach is effective when solving dynamic problems for porous media. Because of the presence of a second longitudinal wave with low velocity, the use of difference schemes in all the coordinates to obtain stable solutions requires a small step consistent both with respect to time and space, which inevitably increases the execution time.