Numerical-analytical algorithm of seismic wave propagation in inhomogeneous media Original Research Article

This paper is the latest version of the numerical-analytical algorithm for solving the forward seismic problem. The main concept of the algorithm is in splitting of 2D and 3D seismic problems to a series of 1D problems with the help of finite integral Fourier, Fourier-Bessel, or Legendre transforms and with the use of a 1D finite difference technique (see [1–3]). In this paper, instead of the finite difference method with respect to one spatial coordinate and time, we propose to use the finite difference approximation with respect to only the spatial coordinate with a subsequent analytical solution of the obtained system of ordinary differential equations.