Bernd Kummer, Alfred Behle; Second-order finite-difference modeling of SH-wave propagation in laterally inhomogeneous media. Bulletin of the Seismological Society of America ; 72 (3): 793–808. doi:
Download citation file:
Following a homogeneous formulation, a new finite-difference (FD) representation of the SH-wave field at first-order lateral discontinuities is introduced. The new scheme, in contrast to previous schemes, represents a consistent second-order approximation of the truncation error. By means of a new, generalized sufficient criterion of stability, the time increment can always be chosen so that the numerical errors of the new scheme remain bounded. Consequently, this FD approximation necessarily converges to the exact solution of the corresponding initial value problem. A numerical example of a Ricker wavelet, vertically incident on a basin structure, illustrates that the new FD scheme converges faster than previous schemes.