Large-scale highly-accurate extended full waveform inversion using convergent Born series

Full-waveform inversion (FWI) is a seismic imaging method that provides quantitative inference about subsurface properties with a wavelength-scale resolution. Its frequency-domain formulation is computationally efficient when processing only a few discrete frequencies. However, classical FWI, which is formulated on the reduced-parameter space, requires starting the inversion with a sufficiently-accurate initial model and low frequency to prevent being stuck in local minima due to cycle skipping. FWI with extended search space has been proposed to mitigate this issue. It contains two main steps: first, data-assimilated (DA) wavefields are computed by allowing for wave-equation errors to match the data at receivers closely. Then, subsurface parameters are estimated from these wavefields by minimizing the wave-equation errors. The DA wavefields are the least-squares solution of an overdetermined system gathering the wave and observation equations. The numerical bandwidth of the resulting normal-equation system is two times that of the wave-equation system, which can be a limiting factor for 3D large-scale applications. Therefore, computing highly accurate DA wavefields at a reasonable computational cost is an issue in extended FWI. This issue is addressed here by rewriting the normal system such that its solution can be computed by solving the time-harmonic wave equation several times in sequence. Moreover, the computational burden of multi-right-hand side (RHS) simulations is mitigated with a sketching method. Finally, we solve the time-harmonic wave equation with the convergent Born series method, which conciliates accuracy and computational efficiency. Application of the new extended FWI algorithm on the salt benchmark shows that it reconstructs at a reasonable cost subsurface models that are similar to those obtained with the classical extended FWI.