Fast full-wavefield seismic inversion using encoded sources

Full-wavefield seismic inversion (FWI) estimates a subsurface elastic model by iteratively minimizing the difference between observed and simulated data. This process is extremely computationally intensive, with a cost comparable to at least hundreds of prestack reverse-time depth migrations. When FWI is applied using explicit time-domain or frequency-domain iterative-solver-based methods, the seismic simulations are performed for each seismic-source configuration individually. Therefore, the cost of FWI is proportional to the number of sources. We have found that the cost of FWI for fixed-spread data can be significantly reduced by applying it to data formed by encoding and summing data from individual sources. The encoding step forms a single gather from many input source gathers. This gather represents data that would have been acquired from a spatially distributed set of sources operating simultaneously with different source signatures. The computational cost of FWI using encoded simultaneous-source gathers is reduced by a factor roughly equal to the number of sources. Further, this efficiency is gained without significantly reducing the accuracy of the final inverted model. The efficiency gain depends on subsurface complexity and seismic-acquisition parameters. There is potential for even larger improvements of processing speed.