GeoInv3D: A Scalable Forward Modeling Framework for Full Waveform Inversion Problem

Full Waveform Inversion is a high-resolution imaging method that has raised considerable interest in the oil industry since a decade. It has been mainly used as a P-wave velocity modeling building tool, while extension to multi-parameter elastic anisotropic reconstruction is now an active field of research. In this context, designing computationally-efficient and versatile FWI softwares, in which different numerical schemes for seismic modeling, wave physics and optimization algorithms can be interfaced easily, is of crucial interest. In this study we describe an object oriented framework based on the definition of abstract interfaces of components involved in frequency or time domain FWI workflow. We point out the numerical modeling of the full seismic wavefield as the central component and that the proposed design is suitable to interface different modeling approaches involving different discretisations and physics with the optimization kernel. Then we demonstrate the capability of the framework to preserve parallel scalability and efficiency of kernels even in an object oriented programming context. Lastly we present a concrete realisation of this abstract framework via an application of an acoustic 3D time-domain FWI on the Valhall field using a staggered grid finite difference scheme.