2.5-D modelling of elastic waves using the pseudospectral method

SUMMARY The pseudospectral method is a high-accuracy numerical modelling technique which requires less memory and CPU time than traditional numerical modelling techniques, such as finite difference. However, full 3-D modelling is still computationally intensive. As we must compromise between realism and computational efficiency, we have developed a 2.5-D pseudospectral scheme for calculating 3-D elastic wavefields in media varying in two dimensions. To demonstrate the feasibility of the scheme described in this paper, we use data from a refraction survey carried out in 1984 by the Explosion Seismology Group in the southeastern foothills of the Hidaka Mountains, Hokkaido, Japan.

[1]  Bengt Fornberg,et al.  The pseudospectral method; accurate representation of interfaces in elastic wave calculations , 1988 .

[2]  Jr. John W. Stockwell.  2.5-D wave equations and high-frequency asymptotics , 1995 .

[3]  Ekkehart Tessmer,et al.  3-D elastic modeling with surface topography by a Chebychev spectral method , 1994 .

[4]  Moshe Reshef,et al.  Three-dimensional elastic modeling by the Fourier method , 1988 .

[5]  M. Korn,et al.  Incorporation of attenuation into time-domain computations of seismic wave fields , 1987 .

[6]  L. Eskola,et al.  THE SOLUTION OF THE STATIONARY ELECTRIC FIELD STRENGTH AND POTENTIAL OF A POINT CURRENT SOURCE IN A 2 1/2‐ DIMENSIONAL ENVIRONMENT* , 1981 .

[7]  R. Kosloff,et al.  Absorbing boundaries for wave propagation problems , 1986 .

[8]  Takashi Furumura,et al.  Three-Dimensional Q Structure in and Around the Hidaka Mountains, Hokkaido, Japan , 1990 .

[9]  Jack K. Cohen,et al.  Two and one‐half dimensional Born inversion with an arbitrary reference , 1987 .

[10]  Edip Baysal,et al.  Forward modeling by a Fourier method , 1982 .

[11]  George A. McMechan,et al.  3-D physical modeling and pseudospectral simulation of seismic common-source data volumes , 1993 .

[12]  Ove Skovgaard,et al.  A spectral method for seismic wave propagation in elastic media , 1990 .

[13]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[14]  Taro Okamoto.  Effects of Sedimentary Structure and Bathymetry Near the Source On Teleseismic P Waveforms From Shallow Subduction Zone Earthquakes , 1993 .

[15]  Keiiti Aki,et al.  Three-dimensional seismic velocity anomalies and their relation to local seismicity , 1979 .

[16]  P. Williamson,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part 1, 2.5-D modeling method , 1995 .

[17]  George A. McMechan,et al.  Computation of multi-attribute seismic wavefields by solution of the elastodynamic equations , 1992 .

[18]  Dean Clifford Witte.  The Pseudospectral Method for Simulating Wave Propagation. , 1989 .

[19]  Pierre-Yves Bard,et al.  The seismic response of sediment-filled valleys. Part 2. The case of incident P and SV waves , 1980 .

[20]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[21]  J. E. Luco,et al.  Three-dimensional response of a cylindrical canyon in a layered half-space , 1990 .

[22]  George A. McMechan,et al.  Two-dimensional elastic pseudo-spectral modeling of wide-aperture seismic array data with application to the Wichita Uplift-Anadarko Basin region of southwestern Oklahoma , 1990, Bulletin of the Seismological Society of America.

[23]  J. Carcione.  The wave equation in generalized coordinates , 1994 .

[24]  Takeo Moriya,et al.  VELOCITY STRUCTURE AND AFTERSHOCK DISTRIBUTION OF THE 1982 URAKAWA-OKI EARTHQUAKE , 1987 .

[25]  N. Bleistein.  Two-and-One-Half Dimensional In-Plane Wave Propagation. , 1984 .

[26]  A. R. Gregory.  Aspects of Rock Physics From Laboratory and Log Data that are Important to Seismic Interpretation: Section 1. Fundamentals of Stratigraphic Interpretation of Seismic Data , 1977 .

[27]  Hiroshi Takenaka,et al.  Effect of 2-D topography on the 3-D seismic wavefield using a 2.5-D discrete wavenumber-boundary integral equation method , 1996 .

[28]  Keiiti Aki,et al.  Surface motion of a layered medium having an irregular interface due to incident plane SH waves , 1970 .

[29]  L. W. Braile,et al.  Reply to J. Vidale's “Comment on ‘A comparison of finite-difference and fourier method calculations of synthetic seismograms’” , 1989, Bulletin of the Seismological Society of America.

[30]  Hiroshi Takenaka,et al.  A wraparound elimination technique for the pseudospectral wave synthesis using an antiperiodic extension of the wavefield , 1995 .

[31]  Robert B. Herrmann,et al.  SH-wave generation by dislocation sources—A numerical study , 1979, Bulletin of the Seismological Society of America.

[32]  Takeo Moriya,et al.  VELOCITY STRUCTURE BENEATH THE HIDAKA MOUNTAINS IN HOKKAIDO, JAPAN , 1984 .

[33]  東京大学地震研究所爆破地震動研究グル-プ.  Explosion Seismic Observations on the Southwestern Foot of the Hidaka Mountains, Hokkaido, Japan : Niikappu-Samani Profile , 1988 .

[34]  Cengiz Esmersoy,et al.  Reverse-time wave-field extrapolation, imaging, and inversion , 1988 .

[35]  B. Fornberg.  The pseudospectral method: Comparisons with finite differences for the elastic wave equation , 1987 .