Assessing Accelerator-Based HPC Reverse Time Migration

• Published in: IEEE transactions on parallel and distributed systems 2011-01, Vol.22 (1), p.147-162
• Main Authors: Araya-Polo, M, Cabezas, J, Hanzich, M, Pericas, M, Rubio, F, Gelado, I, Shafiq, M, Morancho, E, Navarro, N, Ayguade, E, Cela, J M, Valero, M
• Format: Article
• Language: English
• Subjects: Accelerators; Cell/B.E; Computation; Computer architecture; Economics; FPGA; geophysics; GPU; Graphics processing unit; Hardware; Heat transfer; Instruction sets; Mathematical model; Mathematical models; Microprocessors; Migration; Natural gas utilities; Platforms; Propagation; Resin transfer molding; Reverse time migration; Seismic phenomena; Studies
• ISSN: 1045-9219; 1558-2183
• DOI: 10.1109/TPDS.2010.144

Oil and gas companies trust Reverse Time Migration (RTM), the most advanced seismic imaging technique, with crucial decisions on drilling investments. The economic value of the oil reserves that require RTM to be localized is in the order of 10 13 dollars. But RTM requires vast computational power, which somewhat hindered its practical success. Although, accelerator-based architectures deliver enormous computational power, little attention has been devoted to assess the RTM implementations effort. The aim of this paper is to identify the major limitations imposed by different accelerators during RTM implementations, and potential bottlenecks regarding architecture features. Moreover, we suggest a wish list, that from our experience, should be included as features in the next generation of accelerators, to cope with the requirements of applications like RTM. We present an RTM algorithm mapping to the IBM Cell/B.E., NVIDIA Tesla and an FPGA platform modeled after the Convey HC-1. All three implementations outperform a traditional processor (Intel Harpertown) in terms of performance (10x), but at the cost of huge development effort, mainly due to immature development frameworks and lack of well-suited programming models. These results show that accelerators are well positioned platforms for this kind of workload. Due to the fact that our RTM implementation is based on an explicit high order finite difference scheme, some of the conclusions of this work can be extrapolated to applications with similar numerical scheme, for instance, magneto-hydrodynamics or atmospheric flow simulations.