Abstract
Seismic imaging techniques based on two-way wave equations are computationally and data-intensive activities in the oil and gas industry. For instance, Reverse Time Migration (RTM) migrates a set of SEG-Y format data from the disk called a seismogram. Besides, during execution, the RTM application needs to store the forward-propagated wavefield (or source wavefield) on disk to build the final seismic image. Storing the source wavefield for multiple RTMs executing in parallel is even more challenging because the storage capacity can reach tens of Terabytes of information. Aiming to mitigate the storage demand, we develop a 3-D RTM with source wavefield reconstruction. The source wavefield is reconstructed by introducing a new wave equation to the problem and adjusting the initial and boundary conditions to take advantage of random boundary conditions’ (RBC) properties. The RBC does not suppress unwanted waves coming from the artificial boundary enabling full wavefield recovery. We also develop a hybrid OpenACC/MPI implementation for the 3-D RTM on a multi-GPU machine. We test the RTM implementation on the 3-D HPC4E Seismic Test Suite. The numerical experiments show that the OpenACC/MPI 3-D RTM, which implements the wavefield reconstruction, presents the best execution times and hard disk demands.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
Barbosa, C.H., Coutinho, A.L.: Enhancing reverse time migration: hybrid parallelism plus data compression. In: Proceedings of the XLI Ibero-Latin-American Congress on Computational Methods in Engineering. ABMEC (2020)
Barbosa, C.H., Coutinho, A.L.: Seismic modeling and migration with random boundaries on the NEC SX-Aurora TSUBASA. arXiv preprint arXiv:2204.03380 (2022)
Barbosa, C.H., et al.: A workflow for seismic imaging with quantified uncertainty. Comput. Geosci. 145, 104615 (2020)
Cerjan, C., Kosloff, D., Kosloff, R., Reshef, M.: A nonreflecting boundary condition for discrete acoustic and elastic wave equations. Geophysics 50(4), 705–708 (1985)
Clapp, R.G.: Reverse time migration: saving the boundaries. Stanford Exploration Project 137 (2008)
Clapp, R.G.: Reverse time migration with random boundaries. In: SEG Technical Program Expanded Abstracts 2009, pp. 2809–2813. Society of Exploration Geophysicists (2009)
Faria, E.: Migração antes do empilhamento utilizando propagação reversa no tempo, January 1986. http://www.cpgg.ufba.br/pgeof/resumos/gfm/gfm0056a.html
Givoli, D.: Time reversal as a computational tool in acoustics and elastodynamics. J. Comput. Acoust. 22(03), 1430001 (2014)
Komatitsch, D., Martin, R.: An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation. Geophysics 72(5), SM155–SM167 (2007)
Kukreja, N., Hückelheim, J., Louboutin, M., Hovland, P., Gorman, G.: Combining checkpointing and data compression to accelerate adjoint-based optimization problems. In: Yahyapour, R. (ed.) Euro-Par 2019. LNCS, vol. 11725, pp. 87–100. Springer, Cham (2019). https://doi.org/10.1007/978-3-030-29400-7_7
Kukreja, N., Hückelheim, J., Louboutin, M., Washbourne, J., Kelly, P.H., Gorman, G.J.: Lossy checkpoint compression in full waveform inversion: a case study with ZFPv0. 5.5 and the overthrust model. Geosci. Model Dev. 15(9), 3815–3829 (2022)
Kushida, N., Lin, Y.-T., Nielsen, P., Le Bras, R.: Acceleration in acoustic wave propagation modelling using OpenACC/OpenMP and its hybrid for the global monitoring system. In: Wienke, S., Bhalachandra, S. (eds.) WACCPD 2019. LNCS, vol. 12017, pp. 25–46. Springer, Cham (2020). https://doi.org/10.1007/978-3-030-49943-3_2
Li, Q., Fu, L.Y., Wu, R.S., Du, Q.: Efficient acoustic reverse time migration with an attenuated and reversible random boundary. IEEE Access 8, 34598–34610 (2020)
Lindstrom, P., Chen, P., Lee, E.J.: Reducing disk storage of full-3D seismic waveform tomography (F3DT) through lossy online compression. Comput. Geosci. 93, 45–54 (2016)
Liu, H., et al.: The issues of prestack reverse time migration and solutions with graphic processing unit implementation. Geophys. Prospect. 60(5), 906–918 (2012)
Nguyen, B.D., McMechan, G.A.: Five ways to avoid storing source wavefield snapshots in 2D elastic prestack reverse time migration. Geophysics 80(1), S1–S18 (2015)
Pasalic, D., McGarry, R.: Convolutional perfectly matched layer for isotropic and anisotropic acoustic wave equations. In: SEG Technical Program Expanded Abstracts 2010, pp. 2925–2929. Society of Exploration Geophysicists (2010)
Qawasmeh, A., Hugues, M.R., Calandra, H., Chapman, B.M.: Performance portability in reverse time migration and seismic modelling via OpenACC. Int. J. High Perform. Comput. Appl. 31(5), 422–440 (2017)
Schuster, G.T.: Seismic Inversion. Society of Exploration Geophysicists, Tulsa (2017)
Serpa, M.S., et al.: Energy efficiency and portability of oil and gas simulations on multicore and graphics processing unit architectures. Concurr. Comput. Pract. Exp. 33(18), e6212 (2021)
Silva, K.C.: Modelagem, mrt e estudos de iluminação empregando o conceito de dados sísmicos blended, January 2012. http://www.coc.ufrj.br/pt/dissertacoes-de-mestrado/112-msc-pt-2012/2265-karen-carrilho-da-silva
Sun, W., Fu, L.Y.: Two effective approaches to reduce data storage in reverse time migration. Comput. Geosci. 56, 69–75 (2013)
Symes, W.W.: Reverse time migration with optimal checkpointing. Geophysics 72(5), SM213–SM221 (2007)
Wang, Y.: Frequencies of the Ricker wavelet. Geophysics 80(2), A31–A37 (2015)
Zand, T., Malcolm, A., Gholami, A., Richardson, A.: Compressed imaging to reduce storage in adjoint-state calculations. IEEE Trans. Geosci. Remote Sens. 57(11), 9236–9241 (2019)
Zhou, H.W., Hu, H., Zou, Z., Wo, Y., Youn, O.: Reverse time migration: a prospect of seismic imaging methodology. Earth Sci. Rev. 179, 207–227 (2018)
Acknowledgements
This study was financed in part by CAPES, Brazil Finance Code 001. This work is also partially supported by FAPERJ, CNPq, and Petrobras. Computer time on Santos Dumont machine at the National Scientific Computing Laboratory (LNCC - Petrópolis).
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2022 The Author(s), under exclusive license to Springer Nature Switzerland AG
About this paper
Cite this paper
Barbosa, C.H.S., Coutinho, A.L.G.A. (2022). Multi-GPU 3-D Reverse Time Migration with Minimum I/O. In: Navaux, P., Barrios H., C.J., Osthoff, C., Guerrero, G. (eds) High Performance Computing. CARLA 2022. Communications in Computer and Information Science, vol 1660. Springer, Cham. https://doi.org/10.1007/978-3-031-23821-5_12
Download citation
DOI: https://doi.org/10.1007/978-3-031-23821-5_12
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-031-23820-8
Online ISBN: 978-3-031-23821-5
eBook Packages: Computer ScienceComputer Science (R0)