Practical aspects of reflectivity modeling

This paper is intended to help those not familiar with the “lore” of layered earth modeling to avoid some common problems. In the computation of the reflectivity function, an easily incorporated phase‐integral approximation is used away from turning points when the velocity gradient is smaller than the frequency. Hanning windows, or segments thereof, work well for both the slowness integral and the frequency integral. For the quadrature of the slowness integral the Filon method of Frazer is easily coded and vectorizes well; Levin’s Filon method and the Clenshaw‐Curtis‐Filon method of Xu and Mal are more difficult to vectorize, but more powerful because they require fewer evaluations of the reflectivity function. A modification of Strick’s power law is a convenient way to calculate complex frequency‐dependent seismic velocities. The complex frequency technique for avoiding time aliasing is explained by use of the Poisson sum formula. In writing code for vector computers, such as the CRAY, if frequency‐inde...

[1]  H. Takeuchi,et al.  Seismic Surface Waves , 1972 .

[2]  Don L. Anderson,et al.  Velocity dispersion due to anelasticity; implications for seismology and mantle composition , 1976 .

[3]  C. Chapman,et al.  Diffracted Seismic Signals and Their Numerical Solution , 1972 .

[4]  Robert A. Phinney,et al.  Rapid generation of synthetic seismograms in layered media by vectorization of the algorithm , 1987 .

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

[6]  L. Neil Frazer,et al.  On a generalization of Filon's method and the computation of the oscillatory integrals of seismology , 1984 .

[7]  Lawrence W. Braile,et al.  Guide to the Interpretation of Crustal Refraction Profiles , 1975 .

[8]  M. E. O'Neill,et al.  Causal absorption: Its effect on synthetic seismograms computed by the reflectivity method , 1979, Bulletin of the Seismological Society of America.

[9]  J. Gettrust,et al.  Ophiolites, synthetic seismograms, and ocean crustal structure: 1. Comparison of ocean bottom seismometer data and synthetic seismograms for the Bay of Islands ophiolite , 1982 .

[10]  Robert A. Phinney,et al.  Theoretical calculation of the spectrum of first arrivals in layered elastic mediums , 1965 .

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

[12]  Ajit K. Mal,et al.  An adaptive integration scheme for irregularly oscillatory functions , 1985 .

[13]  David Levin,et al.  Procedures for computing one- and two-dimensional integrals of functions with rapid irregular oscillations , 1982 .

[14]  C. H. Chapman,et al.  Diffraction of P Waves by the Core and an Inhomogeneous Mantle , 1970 .

[15]  D. L. Anderson,et al.  Importance of Physical Dispersion in Surface Wave and Free Oscillation Problems: Review (Paper 6R0680) , 1977 .

[16]  A. Papoulis,et al.  The Fourier Integral and Its Applications , 1963 .

[17]  G. Franssens.  Calculation of the elasto-dynamic Green's function in layered media by means of a modified propagator matrix method , 1983 .

[18]  B. Kennett,et al.  Seismic waves in a stratified half-space — IV: P—SV wave decoupling and surface wave dispersion , 1983 .

[19]  N. Kerry.  Synthesis of seismic surface waves , 1981 .

[20]  U. Ascher,et al.  Calculation of complete theoretical seismograms in vertically varying media using collocation methods , 1983 .

[21]  B. Kennett,et al.  Seismic waves in a stratified half space—III. Piecewise smooth models , 1981 .

[22]  R. C. Y. Chin,et al.  Matrix methods in synthetic seismograms , 1984 .

[23]  E. Strick,et al.  The Determination of Q, Dynamic Viscosity and Transient Creep Curves from Wave Propagation Measurements* , 1967 .

[24]  Bernard Budiansky,et al.  Measures of dissipation in viscoelastic media , 1978 .

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

[26]  Gerard J. Fryer,et al.  A slowness approach to the reflectivity method of seismogram synthesis , 1980 .

[27]  G. J. Tango,et al.  Efficient global matrix approach to the computation of synthetic seismograms , 1986 .

[28]  V. Cormier.  The synthesis of complete seismograms in an earth model specified by radially inhomogeneous layers , 1980 .

[29]  G. L. Choy.  Theoretical seismograms of core phases calculated by frequency‐dependent full wave theory, and their interpretation , 1977 .