Waveform Imaging Based on Linear Forward Representations for Scalar Wave Seismic Data

The current reverse-time migration, which is based on wave equations for imaging wavefields, employs an imaging formula derived from Claerbout's imaging principle. This imaging formula is only valid for plane waves with small incident angles on the perfectly flat reflecting surface. However, the complexity of seismic wave propagation may lead to situations that do not meet this requirement. Therefore, this paper divides the subsurface into local scattering and reflecting bodies. It proposes linear forward representations for scattering and reflection data based on perturbations in the physical parameters and wave impedance, respectively. To further describe the effect on the reflecting body boundary, the local reflection coefficient is defined and the linear forward representation for the reflection data based on it is obtained. After that, the proposed linear forward representations are used as the forward equations for the linear inverse of the seismic data, and the seismic data waveform imaging method is developed based on linear inversion theory. At the same time, the specific waveform imaging calculation formulas for scalar wave scattering data and scalar wave reflection data are provided and validated via numerical experiments. Compared with the current reverse-time migration, waveform migration not only has the correct phase and higher resolution in theory but also does not increase the computational complexity. To some extent, it improves the deficiencies of the current structural imaging and provides a basis for subsurface lithological imaging.