Gaussian packets in the computation of seismic wavefields

Ludek Klimes


A method of computing synthetic seismograms by means of the three-parametric summation of Gaussian packets over space is proposed. The expansion of the wavefield specified in terms of a complex-valued amplitude and a phase function into Gaussian packets is very similar to the expansion into Gaussian beams, with the additional summation over points along the rays. The initial conditions for the wavefield, for which the initial directions of the rays are time-independent (e.g. the point source), are assumed. The summation is performed in the frequency domain. The synthetic seismograms are then evaluated by the fast Fourier transform for any source time function. The method should not require a substantial increase of the computational time with respect to the ray and Gaussian beam methods.


Seismic waves, summation of Gaussian packets, space-time Gaussian beams, high-frequency asymptotics, paraxial approximation.

Whole paper

The image of the reprint in GIF 150dpi (5456 kB !) is designed for an instant screen preview.

The scanned image of the reprint in PDF 300dpi (1155 kB !) is designed for printing.

Geophys. J. int., 99 (1989), 421-433.
SW3D - main page of consortium Seismic Waves in Complex 3-D Structures .