Compressed Implicit Jacobian Scheme for Elastic Full-Waveform Inversion

A. Abubakar,M. Li,Y. Lin,T. M. Habashy

DOI: https://doi.org/10.1111/j.1365-246x.2012.05439.x

2011-01-01

Proceedings

Abstract:We present a regularized GaussNewton (GN) inversion method for solving the elastic full-waveform inversion problem in the frequency domain. The main bottleneck of this method is the Jacobian matrix storage and the computational cost of calculating the GN step (the inner-loop calculation). In this work, we managed to reduce the memory usage by calculating the Jacobian matrix on the fly in each inner-loop iteration. By doing so the computational cost of calculating the GN step increases; however, this overhead is mitigated by compressing the field matrices using the adaptive cross approximation scheme. For some cases, this compressed implicit Jacobian scheme may even speed-up the GN step calculation and further regularizes the GN method. As examples, we present inversion results of cross-well seismic and surface seismic data.