A numerical laboratory for simulation and visualization of seismic wavefields1

Computational seismic modelling (CSM) plays an important role in the geophysical industry as an established aid to seismic interpreters. Numerical solution of the elastic wave equations has proved to be a very important tool for geophysicists in both forward modelling and migration. Among the techniques generally used in CSM, we consider the finite-element method (FEM) and investigate its computational and visualization requirements. The CSMFEM program, designed for this purpose and developed on an IBM 3090 computer with vector facility, is described in detail. It constitutes a numerical laboratory for performing computer experiments. Two Newmark type algorithms for time integration are compared with other time integration schemes, and both direct and iterative methods for solving the corresponding large sparse system of linear algebraic equations are analysed. Several numerical experiments to simulate seismic energy propagation through heterogeneous media are performed. Synthetics in the form of common shot gathers, vertical seismic profiles and snapshots are suitably displayed, since with the large amounts of data obtained from CSM research, methods for visualization of the computed results must be developed. The FEM is compared with other numerical tools, such as finite-difference and pseudo-spectral methods.