One-way wave equation migration of common-offset vector gathers: parallel multi CPU/GPU implementation

We present an original algorithm for seismic imaging, based on the depth wavefield extrapolation by one-way wave equation. Parallel implementation of the algorithm is based on the several levels of parallelism. Parallelism of the input data parallelism allows processing full coverage for some area (up to one square km). Mathematical approach allows dealing with each frequency independently and treating solution layer-by-layer, thus a set of 2D instead of initial 3D common-offset vector gathers are processed simultaneously. Next, each common-offset vector image can be stacked and stored independently. As the result, we designed and implemented parallel algorithm which allows computing common-offset vector images using one-way wave equation-based amplitude preserving migration.