Accelerating Kirchhoff Migration on GPU Using Directives

Author(s):  
Rengan Xu ◽  
Maxime Hugues ◽  
Henri Calandra ◽  
Sunita Chandrasekaran ◽  
Barbara Chapman
Keyword(s):  
2000 ◽  
Vol 43 (2-4) ◽  
pp. 281-295 ◽  
Author(s):  
Mark L. Moran ◽  
Roy J. Greenfield ◽  
Steven A. Arcone ◽  
Allan J. Delaney

Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1592-1603 ◽  
Author(s):  
Yonghe Sun ◽  
Fuhao Qin ◽  
Steve Checkles ◽  
Jacques P. Leveille

A beam implementation is presented for efficient full‐volume 3-D prestack Kirchhoff depth migration of seismic data. Unlike conventional Kirchhoff migration in which the input seismic traces in time are migrated one trace at a time into the 3-D image volume for the earth’s subsurface, the beam migration processes a group of input traces (a supergather) together. The requirement for a supergather is that the source and receiver coordinates of the traces fall into two small surface patches. The patches are small enough that a single set of time maps pertaining to the centers of the patches can be used to migrate all the traces within the supergather by Taylor expansion or interpolation. The migration of a supergather consists of two major steps: stacking the traces into a τ-P beam volume, and mapping the beams into the image volume. Since the beam volume is much smaller than the image volume, the beam migration cost is roughly proportional to the number of input supergathers. The computational speedup of beam migration over conventional Kirchhoff migration is roughly proportional to [Formula: see text], the average number of traces per supergather, resulting a theoretical speedup up to two orders of magnitudes. The beam migration was successfully implemented and has been in production use for several years. A factor of 5–25 speedup has been achieved in our in‐house depth migrations. The implementation made 3-D prestack full‐volume depth imaging feasible in a parallel distributed environment.


Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1270-1274 ◽  
Author(s):  
Le‐Wei Mo ◽  
Jerry M. Harris

Traveltimes of direct arrivals are obtained by solving the eikonal equation using finite differences. A uniform square grid represents both the velocity model and the traveltime table. Wavefront discontinuities across a velocity interface at postcritical incidence and some insights in direct‐arrival ray tracing are incorporated into the traveltime computation so that the procedure is stable at precritical, critical, and postcritical incidence angles. The traveltimes can be used in Kirchhoff migration, tomography, and NMO corrections that require traveltimes of direct arrivals on a uniform grid.


Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. S211-S219 ◽  
Author(s):  
Siwei Li ◽  
Sergey Fomel

The computational efficiency of Kirchhoff-type migration can be enhanced by using accurate traveltime interpolation algorithms. We addressed the problem of interpolating between a sparse source sampling by using the derivative of traveltime with respect to the source location. We adopted a first-order partial differential equation that originates from differentiating the eikonal equation to compute the traveltime source derivatives efficiently and conveniently. Unlike methods that rely on finite-difference estimations, the accuracy of the eikonal-based derivative did not depend on input source sampling. For smooth velocity models, the first-order traveltime source derivatives enabled a cubic Hermite traveltime interpolation that took into consideration the curvatures of local wavefronts and can be straightforwardly incorporated into Kirchhoff antialiasing schemes. We provided an implementation of the proposed method to first-arrival traveltimes by modifying the fast-marching eikonal solver. Several simple synthetic models and a semirecursive Kirchhoff migration of the Marmousi model demonstrated the applicability of the proposed method.


2021 ◽  
Author(s):  
Vadym Plakhtii ◽  
Oleksandr Dumin ◽  
Oleksandr Pryshchenko

2021 ◽  
Author(s):  
Yu Pu ◽  
Gang Liu ◽  
Diancheng Wang ◽  
Hui Huang ◽  
Ping Wang

Geophysics ◽  
1999 ◽  
Vol 64 (6) ◽  
pp. 1793-1805 ◽  
Author(s):  
Herman H. Jaramillo ◽  
Norman Bleistein

The Kirchhoff approximation provides a representation of seismic data as a summation of imaged data along isochron surfaces (demigration). The asymptotic inversion of this representation provides a migration as a summation of seismic data along diffraction surfaces. We replace Born inversion techniques with Kirchhoff inversion techniques and further show the link between the Kirchhoff and Born representations after the Born linearized reflection coefficient is replaced by the Kirchhoff reflection coefficient.


Sign in / Sign up

Export Citation Format

Share Document