Geophysics , 81 ( 5 ) , T 197-T 209 Wavefield reconstruction by interpolating significantly decimated boundaries

Pengliang Yang,Romain Brossier,Jean Virieux

2016-01-01

Abstract:Many practical seismic applications such as reverse time migration (RTM) and full waveform inversion (FWI) are usually computation and memory intensive. To perform cross-correlation in RTM or build the gradient for FWI, it is mandatory to access both forward wavefield and adjoint wavefield simultaneously. To do this, there are three methods: one is to read the stored forward wavefield from the disk, second using final snapshot and the stored boundaries via reverse propagation, and third is remodeling using checkpointing from stored state to another state. Among these techniques, wavefield reconstruction by reverse propagation appears to be quite straightforward approach, however, suffering a stringent memory bottleneck for 3D large scale imaging applications. The Courant-Friedrichs-Lewy (CFL) condition is a fundamental criterion to determine temporal sampling to achieve stable wavefield extrapolation. The injection of the boundary sequence in time is essentially determined by Nyquist sampling principle, rather than the time interval given by CFL, which is much smaller than the Nyquist requirement. Based on this recognition, in this paper we present three boundary interpolation techniques, namely the discrete Fourier transform (DFT) interpolation, Kaiser windowed sinc interpolation and Lagrange polynomial interpolation, for wavefield reconstruction to move from CFL to Nyquist limit. The wavefield reconstruction via DFT interpolation can be implemented by folding and unfolding steps in the forward simulation and backward reconstruction on the fly. Compared with the DFT interpolation, the wavefield reconstruction methods using Kaiser windowed sinc interpolation and Lagrange polynomial interpolation have better efficiency while remaining a competitive accuracy. These methods allow us to dramatically decimate the boundary without significant loss of information, and nicely reconstruct the boundary elements in between the samples, making the in-core memory saving of the boundaries feasible in 3D large scale imaging applications.