Multi-scale seismic envelope inversion using a direct envelope Frechet derivative for strong-nonlinear full waveform inversion

Traditional seismic envelope inversion takes use of a nonlinear misfit functional which relates the envelope of seismogram to the observed wavefield records, and then derive the sensitivity kernel of envelope to velocity through the use of waveform Frechet derivative by linearizing the nonlinear functional. We know that the waveform Frechet derivative is based on the Born approximation of wave scattering and can only work for the case of weak scattering. Therefore, the traditional envelope inversion using waveform Frechet derivative has severe limitation in the case of strong-scattering for large-scale, strong-contrast inclusions. In this paper, we derive a new direct envelope Frechet derivative (sensitivity operator) based on energy scattering physics without using the chain rule of differentiation. This new envelope Frechet derivative does not have the weak scattering assumption for the wavefield and can be applied to the case of strong-nonlinear inversion involving salt structures. We also extend the envelope data into multi-scale window-averaged envelope (WAE), which contains much rich low-frequency components corresponding to the long-wavelength velocity structure of the media. Then we develop a joint inversion which combines the new multi-scale direct envelope inversion (MS-DEI) with standard FWI to recover large-scale strong-contrast velocity structure of complex models. Finally, we show some numerical examples of its successful application to the inversion of a 1-D thick salt-layer model and the 2D SEG/EAGE salt model.