Numerical Methods in Seismic Wave Propagation

Matt A. D. McDonald

The numerical modelling of wave equations is a common theme in many seismic applications, and is an important tool in understanding how the physical systems of interest react in the process of a seismic experiment. In this thesis we apply state-of-the-art numerical methods based on domain-decomposition combined with local pseudospectral spatial discretization, to three physically realistic models of seismic waves, namely their propagation in acoustic, elastic, and viscoelastic media. The Galerkin formulation solves the weak form of the partial differential equation representing wave propagation and naturally includes boundary integral terms to represent free surface, rigid, and absorbing boundary effects. Stability, accuracy, and computation issues are discussed in this context along with direct comparison with finite difference methodologies.

This works is an effort to bridge the gap between the development of accurate physical models to represent the real world, as seen in seismic modelling, and the implementation of modern numerical techniques for the accurate solutions of partial differential equations.