Workstation computation of synthetic seismograms for vertical and horizontal profiles: A full wavefield response for a two-dimensional layered half-space

FORTRAN code for generation of full wavefield synthetic seismograms is presented for two-dimensional horizontally layered models bounded by a free surface and a half space. Model layers are user defined by compressional and shear velocities, Q factors, densities and thicknesses. The algorithm is based on the reflectivity method and uses the propagator matrix approach. Explosion (point) and double couple (fault) sources are generated with a moment tensor representation. As evaluation of the slowness integrals involves time consuming numerical Hankel transforms, these computations are made with a generalized Filon method that saves computational time. The architecture of the program is unusual because the outermost loop is over temporal frequency and the innermost loop is over slowness. This permits the use of frequency-dependent seismic velocities, necessary for causality, while giving a factor of seven speed-up from vectorization. The codes are applicable for both vector computers and workstations. Two test cases demonstrate successful applications of the codes for both horizontal seismic profiles (receivers at one depth at successively larger offsets) and for vertical seismic profiles (receivers arranged in a vertical array at any offset). Receivers and source may be positioned within any layer. The seismograms display direct, refracted, reflected, and head-wave arrivals and their multiples. Mode converted events of compressional and shear propagation are generated and identified. The code generates seismograms for pressure, vertical and horizontal displacement sensors and for models combining acoustic and elastic layers.