3D Frequency-Domain Seismic Modeling With a Block Low-Rank Algebraic Multifrontal Direct Solver.

SUMMARY

Three-dimensional frequency-domain full waveform inversion of fixed-spread data can be efficiently performed in the visco-acoustic approximation when seismic modeling is based on a sparse direct solver. We present a new algebraic Block Low-Rank (BLR) multifrontal solver which provides an approximate solution of the time-harmonic wave equation with a reduced operation count, memory demand and volume of communication relative to the full-rank solver. We show some preliminary simulations in the 3D SEG/EAGE overthrust model, that give some insights on the memory and time complexities of the low-rank solver for frequencies of interest in full-waveform inversion (FWI) applications.