A Spectral Scheme for Viscoelastic Seismic Modeling,
Abstract
The pseudospectral method is especially valuable for seismic modeling because of its high accuracy compared to other numerical techniques. The method can be regarded as a limit of finite difference of increasing orders, and a process of trigonometric interpolation, thus it exhibits high accuracy. Stability of the method is also favorable. Fourier polynomials are especially efficient but have the disadvantage of forcing periodicity, and Chebyshev polynomials are somewhat less efficient but are more flexible in application of boundary conditions. We have used a Fourier pseudospectral method in the horizontal direction and Chebychev polynomials in the vertical direction. Curved grids conforming to the surface topography and major interfaces are made possible by coordinate transformations. A full viscoelastic formulation permits convenient implementation of attenuating layers to reduce wrap-around in the horizontal direction. The result is an efficient method for 2- and 3-D linear viscoelastic wave propagation.
Document Details
- Document Type
- Technical Report
- Publication Date
- Aug 14, 1995
- Accession Number
- ADP204464
Entities
People
- James C. Schatzman
- Zhaobo Meng
Organizations
- University of Wyoming