Optimized Model Parameterization using Compact Full Waveform Inversion

Abstract

Summary Seismic Full Waveform Inversion (FWI) has the potential to produce high-resolution subsurface images, but the computational resources required for realistically sized problems can be prohibitively large. In terms of computational costs, Gauss-Newton algorithms are more attractive than the commonly employed conjugate gradient methods, because the former have favorable convergence properties. However, efficient implementations of Gauss-Newton algorithms require an excessive amount of computer memory for larger problems. To address this issue, we introduce Compact Full Waveform Inversion (CFWI). Here, a suitable inverse model parameterization is sought that allows representing all subsurface features, potentially resolvable by a particular source-receiver deployment, but using only a minimum number of model parameters. In principle, an inverse model parameterization, based on the Eigenvalue decomposition, would be optimal, but this is computationally not feasible for realistic problems. Instead, we present two alternative parameter transformations, namely the Haar and the Hartley transformations, with which similarly good results can be obtained. By means of a suite of numerical experiments, we demonstrate that these transformations allow the number of model parameters to be reduced to only a few percent of the original parameterization without any significant loss of spatial resolution. This facilitates efficient solutions of large-scale FWI problems with explicit Gauss-Newton algorithms.