Accelerating f inite difference wavef ield-continuation depth migration by GPU  

用GPU提速地震资料单程波有限差分叠前深度偏移(英文)

作　　者：刘国峰[1,2] 孟小红[1] 刘洪[2] 

机构地区：[1]中国地质大学地球探测与信息技术教育部重点实验室,北京100083 [2]中国科学院地质与地球物理研究所,北京100029

出　　处：《Applied Geophysics》2012年第1期41-48,115,共9页应用地球物理（英文版）

基　　金：supported by the National Natural Science Foundation of China (Nos. 41104083 and 40804024); Fundamental Research Funds for the Central Universities (No, 2011YYL022)

摘　　要：The most popular hardware used for parallel depth migration is the PC-Cluster but its application is limited due to large space occupation and high power consumption. In this paper, we introduce a new hardware architecture, based on which the finite difference （FD） wavefield-continuation depth migration can be conducted using the Graphics Processing Unit （GPU） as a CPU coprocessor. We demonstrate the program module and three key optimization steps for implementing FD depth migration： memory, thread structure, and instruction optimizations and consider evaluation methods for the amount of optimization. 2D and 3D models are used to test depth migration on the GPU. The tested results show that the depth migration computational efficiency greatly increased using the general-purpose GPU, increasing by at least 25 times compared to the AMD 2.5 GHz CPU.复杂介质情况下,地震波延多路径传播,此时基于波动方程延拓的深度成像方法,相对于Kirchhoff方法能够获得更为精确的成像效果,但是,该深度偏移方法由于高昂的计算消耗阻碍了它在生产中的应用。譬如,叠前深度偏移计算需要大规模的计算机集群,占地面积和电能消耗大。本文介绍了应用一种新的GPU计算架构来辅助CPU进行偏移计算。基于新架构的波动方程深度偏移提高了计算效率,而且机器占地面积和电能消耗也大幅度减少。本文以有限差分波动方程深度偏移为例,介绍了其编程模型和程序优化环节,提高了深度偏移计算效率。2D和3D测试表明,与相同单位个AMD2.5GHz CPU计算相比,该架构下的有限差分波动方程叠前深度偏移计算效率提高至少35倍。

关 键 词：Wavefield-continuation depth migration finite difference Graphic Processing Unit EFFICIENCY 

分 类 号：P631.44[天文地球—地质矿产勘探]