基于GPU加速的地震图像重建技术

Seismic Image Reconstruction Based on GPU Acceleration

摘要:

针对目前地层层析成像算法中正演算法存在计算量大、计算速度慢的问题，以图像处理器(GPU)为核心，研究并实现了一种基于GPU平台的时域有限差分(FDTD)正演算法。CUDA是一种由NVIDIA推出的GPU通用并行计算架构，也是目前较为成熟的GPU并行运算架构。而FDTD正演算法本身在算法特性上满足并行的要求，二者的结合将极大地加速程序的计算速度。在基于标准Marmousi速度模型的正演模拟中，程序速度提升30倍，而GPU正演图像与CPU正演结果误差小于千分之一。算例表明CUDA可以大大加速目前的FDTD正演算法，并且随着GPU硬件自身的发展和计算架构的不断改进，加速效果还将进一步提升，这将有利于后续波形反演工作的进展。

Abstract:

In order to solve the problems of high computation and low speed in traditional waveform tomography algorithm, FDTD algorithm was proposed based on CUDA platform with GPU as the core processor. CUDA is a general parallel computing architecture introduced by NVIDIA and also one of the most popular architectures in GPU calculation. Combined the advantage of strong calculation ability of GPU with the natural parallel characteristic of FDTD algorithm, the performance of the program can be greatly increased. Forward modeling based on Marmousi model shows that the speed is increased by 30 times but the error of the result between GPU and CPU is less than millesimal. Simulation results indicate CUDA can be used in acceleration of FDTD algorithm. As the development of hardware and computing architecture, the accelerating effect can be increased, which will contribute a lot for the research of wave inversion.