Tomographic Full Waveform Inversion Research Articles

Simultaneous inversion of full data bandwidth by tomographic full-waveform inversion

The convergence of full-waveform inversion can be improved by extending the velocity model along either the subsurface-offset axis or the time-lag axis. The extension of the velocity model along the time-lag axis enables us to linearly model large time shifts caused by velocity perturbations. This linear modeling was based on a new linearization of the scalar wave equation in which perturbation of the extended slowness squared was convolved in time with the second time derivative of the background wavefield. The linearization was accurate for reflected events and transmitted events. We determined that it can effectively model conventional reflection data as well as modern long-offset data containing diving waves. It also enabled the simultaneous inversion of reflections and diving waves, even when the starting velocity model was far from being accurate. We solved the optimization problem related to the inversion with a nested algorithm. The inner iterations were based on the proposed linearization and on a mixing of scales between the short- and long-wavelength components of the velocity model. Numerical tests performed on synthetic data modeled on the Marmousi model and on the “Caspian Sea” portion of the well-known BP model demonstrated the global-convergence properties as well as the high-resolution potential of the proposed method.

Tomographic full-waveform inversion (TFWI) by combining FWI and wave-equation migration velocity analysis

Convergence of full waveform inversion can be improved by extending the velocity model along either the subsurface-offset axis or the time-lag axis. The extension of the velocity model along the time-lag axis enables us to linearly model large time shifts caused by velocity perturbations. The extension is based on a new linearization of the scalar wave equation where the extended-velocity perturbation is convolved in time with the Laplacian of the background wavefield. This linearization is accurate for both reflected events and transmitted events and, in particular, for diving waves recorded at large offsets. The modeling capabilities of the proposed linearization enable the simultaneous inversion of reflections and diving waves even when the starting velocity model is far from being accurate. We solve the resulting optimization problem with a nested algorithm. The inner iterations are based on the proposed linearization and on a mixing of scales between the short- and long-wavelength components of the velocity model. Numerical tests performed on two synthetic data sets modeled on the Marmousi model and on the Caspian Sea portion of the well-known BP model show the global-convergence properties as well as the high-resolution potential of the method.