Quasi-Newton full-waveform inversion with a projected Hessian matrix

ABSTRACTWe present a method, in realistic-size full-waveform inversion (FWI), to explicitly construct a projected Hessian matrix and its inverse matrix, which we subsequently used to solve FWI with a quasi-Newton method. Newton’s method is practically unfeasible in solving realistic-size FWI problems because of the prohibitive cost (computing time and memory consumption) of calculating the Hessian matrix and the inverse Hessian. Therefore, the Gauss-Newton method and various quasi-Newton methods are proposed to approximate a Hessian matrix. Particularly, current quasi-Newton FWI (QNFWI) commonly uses the limited-memory BFGS (L-BFGS) method, which, however, only implicitly approximates an inverse Hessian. We repose FWI as a sparse optimization problem in a sparse model space, which contains substantially fewer model parameters that are constrained by structures of the model. With respect to fewer parameters in the sparse model, we can avoid the “limited-memory” approximation and are able to explicitly comp...