A multi-GPU benchmark for 2D Marchenko Imaging

2021 
The Marchenko method allows estimating Green's functions with a virtual source in the subsurface from a reflection response on the surface. It is an inverse problem that can be solved directly or by an iterative scheme, with the latter being more feasible computationally. In this work we present a multi-GPU implementation of a well-established iterative Marchenko algorithm based on (the) Neumann series. The time convolution and space integration performed on each iteration, also referred to as synthesis, are here represented as a segmented dot product, which can be accelerated on modern GPUs through the usage of warp-shuffle instructions and CUDA libraries. The original CPU version is benchmarked on 36 CPU cores versus the implemented version on 4 GPUs, over three different reflection data sets, with sizes ranging from 3 GB to 250 GB.
    • Correction
    • Source
    • Cite
    • Save
    • Machine Reading By IdeaReader
    0
    References
    2
    Citations
    NaN
    KQI
    []