Extraction of the tomography mode with nonstationary smoothing for full-waveform inversion

Full-waveform inversion (FWI) includes migration and tomography modes. The tomographic component of the gradient from reflection data is usually much weaker than the migration component. To use the tomography mode to fix background velocity errors, it is necessary to extract the tomographic component from the gradient. Otherwise, the inversion will be dominated by the migration mode. We have developed a method based on nonstationary smoothing to extract the tomographic component from the raw gradient. By analyzing the characteristics of the scattering angle filtering, the wavenumber of the tomographic component at a given frequency is seen to be smaller than that of the migration component. Therefore, low-wavenumber-pass filtering can be applied to extract the tomographic component. The low-wavenumber-pass smoothing filters are designed with Gaussian filters that are determined by the frequency of inversion, the model velocity, and the minimum scattering angle. Thus, this filtering is nonstationary smoothing in the space domain. Because this filtering is carried out frequency by frequency, it works naturally and efficiently for FWI based on frequency-domain modeling. Furthermore, because the maximum opening angle of the reflections in a typical acquisition geometry is much smaller than the minimum scattering angle for the tomographic component, which is generally set at 160 degrees, there is a relatively large gap between the wavenumbers of the tomographic and migration components. In other words, the nonstationary smoothing can be applied once to a group of frequencies for time-domain FWI without leaking the migration component into the tomographic component. Analyses and numerical tests indicate that two frequency groups are generally sufficient to extract the tomographic component for the typical frequency range of time-domain FWI. The numerical tests also demonstrate that the nonstationary smoothing method is effective and efficient at extracting the tomographic component for reflection waveform inversion.