Method of using CPU-GPU platform for seismic wave reverse-time migration imaging

Abstract

The invention provides a method of using a CPU-GPU platform for seismic wave reverse-time migration imaging, which can be applied to the technical field of reflected wave seismic data processing. The method comprises steps: a main control node generates a task pool according to shot gather data; the main control node carries out reverse-time migration on single shot seismic data according to the GPU for calculating the needed memory size, a CPU-GPU parallel processing strategy is determined, and the CPU-GPU parallel processing strategy is sent to each slave node; the main control node starts each slave node; each slave node receives a tack from the task pool, and according to the CPU-GPU parallel processing strategy, shot seismic data and a speed model included in the current task are used for reverse-time migration calculation, and a single shot reverse-time migration result according to the current task is obtained; the main control node carries out stacked processing on single shot reverse-time migration results according to various tasks so as to obtain a reverse-time migration imaging profile corresponding to the shot gather data. The method makes full use of the CPU-GPU platform parallel calculation ability, and has the advantages of high parallel degree and easy realization, and a shot point wave field simulation process is saved.

Description

technical field [0001] The present invention relates to the technical field of reflected wave seismic data processing, in particular to a method for imaging reverse time migration of seismic waves using a CPU-GPU platform. Background technique [0002] The reverse time migration method based on the full-acoustic wave equation directly solves the full-acoustic partial differential equation with an explicit high-order finite difference algorithm in the time-space domain, and truly simulates the wave propagation phenomenon. This method completely complies with the wave equation and has no inclination limitation. It is applicable to sharp changes in the velocity field and has obvious advantages in imaging complex three-dimensional structures with high imaging accuracy. [0003] Due to the use of the finite difference algorithm and the need to save a large amount of source wave field data, compared with other migration algorithms, the computational cost of the reverse time migrat...

AI Technical Summary

Problems solved by technology

Although this scheme does not have the problem of wave field access, it has the following problems: the wave field of the shot point is calculated twice; in order to reconstruct the wave field of the shot point in reverse time, the energy of the wave field must be kept in the calculation area during the whole process Therefore, it is not suitable to use the absorbing boundary condition; in order to avoid the coherent calculation domain boundary reflection, it is necessary to add a random velocity band on each boundary surface of the calculation domain, and replace the coherent artificial reflection with random artificial reflection, but then The speed band will bring reflection noise, which will affect the calculation results

Images