GPU implementation of minimal dispersion recursive operators for reverse time migration

Summary The implementation of new recursive operators for computation of numerical derivatives results in minimal dispersion in Reverse Time Migration Prestack Depth Migration (i.e. RTM PSDM). Compared to the more commonly used finite difference operators, the presented new method enables imaging of higher frequencies in RTM PSDM. Since RTM PSDM is the modern method of choice for imaging many exploration targets, computer optimization is an integral part of code development. We present here details of the GPU implementation for newly developed spatial derivative operators which enable routine use in industrial settings.