Abstract

This article provides an overview of the application of the staggered-grid finite-difference technique to model wave propagation problems in 3D elastic media. In addition to presenting generalized, discrete representations of the differential equations of motion using the staggered-grid approach, we also provide detailed formulations that describe the incorporation of moment-tensor sources, the implementation of a stable and accurate representation of a planar free-surface boundary for 3D models, and the derivation and implementation of an approximate technique to model spatially variable anelastic attenuation within time-domain finite-difference computations. The comparison of results obtained using the staggered-grid technique with those obtained using a frequency-wavenumber algorithm shows excellent agreement between the two methods for a variety of models. In addition, this article also introduces a memory optimization procedure that allows large-scale 3D finite-difference problems to be computed on a conventional, single-processor desktop workstation. With this technique, model storage is accommodated using both external (hard-disk) and internal (core) memory. To reduce system overhead, a cascaded time update procedure is utilized to maximize the number of computations performed between I/O operations. This formulation greatly expands the applicability of the 3D finite-difference technique by providing an efficient and practical algorithm for implementation on commonly available workstation platforms.

First Page Preview

First page PDF preview
You do not currently have access to this article.