Optimally accurate second-order time-domain finite difference scheme for the elastic equation of motion: one-dimensional case

SUMMARY We previously derived a general criterion for optimally accurate numerical operators for the calculation of synthetic seismograms in the frequency domain (Geller & Takeuchi 1995). We then derived modi¢ed operators for the Direct Solution Method (DSM) (Geller & Ohminato 1994) which satisfy this general criterion, thereby yielding signi¢cantly more accurate synthetics (for any given numerical grid spacing) without increasing the computational requirements (Cummins et al .1 994; Takeuchi, Geller & Cummins 1996; Cummins, Takeuchi & Geller 1997). In this paper, we derive optimally accurate time-domain ¢nite diierence (FD) operators which are second order in space and time using a similar approach. As our FD operators are local, our algorithm is well suited to massively parallel computers. Our approach can be extended to other methods (e.g. pseudo-spectral) for solving the elastic equation of motion. It might also be possible to extend this approach to equations other than the elastic equation of motion, including non-linear equations.

[1]  Robert J. Geller,et al.  DSM complete synthetic seismograms : SH, spherically symmetric, case , 1994 .

[2]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[3]  R. M. Alford,et al.  ACCURACY OF FINITE‐DIFFERENCE MODELING OF THE ACOUSTIC WAVE EQUATION , 1974 .

[4]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[5]  Gerard T. Schuster,et al.  Parsimonious staggered grid finite‐differencing of the wave equation , 1990 .

[6]  Robert J. Geller,et al.  Computation of complete synthetic seismograms for laterally heterogeneous models using the Direct Solution Method , 1997 .

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

[8]  Tatsuhiko Hara,et al.  Two efficient algorithms for iterative linearized inversion of seismic waveform data , 1993 .

[9]  Robert J. Geller,et al.  Computation of synthetic seismograms and their partial derivatives for heterogeneous media with arbitrary natural boundary conditions using the Direct Solution Method , 1994 .

[10]  Bertrand Meyer,et al.  Crustal thickening in Gansu‐Qinghai, lithospheric mantle subduction, and oblique, strike‐slip controlled growth of the Tibet plateau , 1998 .

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

[12]  R. Geller,et al.  A new method for computing highly accurate DSM synthetic seismograms , 1995 .

[13]  A. Tarantola.  Inversion of seismic reflection data in the acoustic approximation , 1984 .

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

[15]  Robert J. Geller,et al.  Highly accurate P-SV complete synthetic seismograms using modified DSM operators , 1996 .

[16]  James Hardy Wilkinson,et al.  Eigenvalues of Ax = λBx with band symmetric A and B , 1969, Comput. J..

[17]  Gene H. Golub,et al.  Matrix computations , 1983 .

[18]  Guoming Xu,et al.  The symmetries of the system matrix and propagator matrix for anisotropic media and of the system matrix forperiodically layered media , 1998 .