Efficient Parallel Implementation of Multi-Arrival 3D Prestack Seismic Depth Migration
The goal of seismic migration is to reconstruct the image of Earth's depth inhomogeneities on the base of seismic data. Seismic data is obtained using shots in shallow wells that are located in a dense grid points. Those shots could be considered as special point sources. A reflected and scattered seismic waves from the depth inhomogeneities are received by geophones located also in a dense grid points on a surface. A seismic image of depth inhomogeneities can be constructed based on these waves. The implementation of 3-D seismic migration implies the solution of about 104÷5 3-D direct problems of wave propagation. Hence efficient asymptotic methods are of a great practical importance. The multi-arrival 3-D seismic migration program is implemented based on a new asymptotic method. It takes into account multi-pass wave propagation and caustics. The program uses parallel calculations in an MPI environment on hundreds and thousands of processor cores. The program was successfully tested on an international synthetic "SEG salt" data set and on real data. A seismic image cube for Timan-Pechora region is given as an example.
Ivanov, A., Khilkov, S.: Aiwlib library as the instrument for creating numerical modeling applications. Scientific Visualization 10(1), 110–127 (2018), DOI: 10.26583/sv.10.1.09
Pleshkevich, A., Ivanov, A., Khilkov, S.: Asymptotic solution of wavefield continuation problem in the ray parametric coordinates. In: SEG Technical Program Expanded Abstracts 2017. pp. 5551–5555 (2017), DOI: 10.1190/segam2017-17633541.1
Pleshkevich, A., Ivanov, A., Levchenko, V., Khilkov, S.: Multiarrival amplitude–preserving prestack 3D depth migration. Russian Geophysics (S), 76–84 (2017)
Semtchenok, N., Popov, P., Verdel, A.: Depth migration by the Gaussian beam summation method. Geophysics (75), 81–93 (2010), DOI: 10.1190/1.3361651
Vainberg, B.: Asymptotic methods in the equations of mathematical physics. Moscow: MSU (1982)