Full text loading...
-
Implementation of Large-Scale Integral Operators with Modern HPC Solutions
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, EAGE 2020 Annual Conference & Exhibition Online, Dec 2020, Volume 2020, p.1 - 5
Abstract
Numerical integral operators of convolution type lie at the foundations of most wave-equation-based methods for processing and imaging of seismic data. Several of such methods require the solution of an inverse problem, which in turn calls for multiple forward and adjoint passes of the modelling operator. In this abstract we provide some insights into the numerical aspects of solving such systems of integral equations and present a framework that leverages open-source libraries for distributed storage and computing as well as for high-level symbolic representation of linear operators. To validate the effectiveness of our implementation, we evaluate the scalability of the forward and adjoint operations of the well-known time-domain multi-dimensional convolution (MDC) operator with respect to increasing size of the input data and number of computational resources. Finally, we use this implementation to solve the Marchenko equations by means of least-squares inversion for a 3D synthetic dataset composed of up to 9801 sources and receivers.