Full waveform inversion (FWI) is an emerging subsurface imaging technique, used to locate oil and gas reservoirs. The key chal- lenges that hinder its adoption by industry are both algorithmic and computational in nature, including storage, communication, and process- ing of large-scale data structures, which impose cardinal impediments upon computational scalability. In this work we will present a complete matrix-free algorithmic formulation of a 3D elastic time domain spectral element solver for both the forward and adjoint wave-fields as part of a greater cloud based FWI framework. We discuss computational opti- misation (SIMD vectorisation, use of Many Integrated Core architec- tures, etc.) and present scaling results for two HPC systems, namely an IBM Blue Gene/Q and an Intel based system equipped with Xeon Phi coprocessors.
['Stephen Moore', 'Devi Sudheer Chunduri', 'Sergiy Zhuk', 'Tigran T. Tchrakian', 'Ewout van den Berg', 'Albert Akhriev', 'Alberto Costa Nogueira', 'Andrew Rawlinson', 'Lior Horesh']
Semi-discrete Matrix-Free Formulation of 3D Elastic Full Waveform Inversion Modeling