An Efficient Numerical Method for Computing Synthetic Seismograms for a Layered Half-space with Sources and Receivers at Close or Same Depths

Abstract — It is difficult to compute synthetic seismograms for a layered half-space with sources and receivers at close to or the same depths using the generalized R/T coefficient method (Kennett, 1983; Luco and Apsel, 1983; Yao and Harkrider, 1983; Chen, 1993), because the wavenumber integration converges very slowly. A semi-analytic method for accelerating the convergence, in which part of the integration is implemented analytically, was adopted by some authors (Apsel and Luco, 1983; Hisada, 1994, 1995). In this study, based on the principle of the Repeated Averaging Method (Dahlquist and Björck, 1974; Chang, 1988), we propose an alternative, efficient, numerical method, the peak-trough averaging method (PTAM), to overcome the difficulty mentioned above. Compared with the semi-analytic method, PTAM is not only much simpler mathematically and easier to implement in practice, but also more efficient. Using numerical examples, we illustrate the validity, accuracy and efficiency of the new method.

[1]  Horace Lamb,et al.  On the propagation of tremors over the surface of an elastic solid , 1904, Proceedings of the Royal Society of London.

[2]  Xiaofei Chen,et al.  A systematic and efficient method of computing normal modes for multilayered half-space , 1993 .

[3]  D. Helmberger.  Generalized ray theory for shear dislocations , 1974, Bulletin of the Seismological Society of America.

[4]  B. Kennett,et al.  Seismic Wave Propagation in Stratified Media , 1983 .

[5]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[6]  Marijan Dravinski,et al.  On evaluation of the green functions for harmonic line loads in a viscoelastic half space , 1988 .

[7]  R. J. Apsel,et al.  On the Green's functions for a layered half-space. Part II , 1983 .

[8]  Keiiti Aki,et al.  Discrete wave-number representation of seismic-source wave fields , 1977, Bulletin of the Seismological Society of America.

[9]  M. Bouchon.  Discrete wave number representation of elastic wave fields in three-space dimensions , 1979 .

[10]  Robert B. Herrmann,et al.  A numerical study of P-, SV-, and SH-wave generation in a plane layered medium , 1980 .

[11]  B. Kennett,et al.  Seismic waves in a stratified half space. , 1979 .

[12]  M. Bouchon.  The complete synthesis of seismic crustal phases at regional distances , 1982 .

[13]  N. Ricker.  Transient waves in visco-elastic media , 1977 .

[14]  D. Harkrider,et al.  A generalized reflection-transmission coefficient matrix and discrete wavenumber method for synthetic seismograms , 1983 .

[15]  Yoshiaki Hisada,et al.  An efficient method for computing Green's functions for a layered half-space with sources and receivers at close depths (part 2) , 1995, Bulletin of the Seismological Society of America.

[16]  K. Fuchs,et al.  Computation of Synthetic Seismograms with the Reflectivity Method and Comparison with Observations , 1971 .