99 2012年12月
现代地质 GEOSCIENCE
应用图形处理器快速计算逆时偏移
刘雅宁，刘国峰，张致付
（中国地质大学地球物理与信息技术学院，北京
100083）
VoL 26Na. 6 Dee.2012
摘要：逆时偏移是目前精度最高的地震数据叠前深度偏移方法，但高强度的计算需求限制了其在工业生产领域的大规模应用。可编程图形处理器的发展为逆时偏移的快速计算提供了一种新的计算选择。围绕如何在图形处理器上开展逆时偏移计算展开，总结了图形处理器计算的优化关键，并根据逆时偏移的特点着重介绍了两个优化环节：一个是应用随机边界条件，以计算换存储，减少数据在主机和图形处理器间的传输；二是应用共享存储器来存储正演计算的波场，相比全局存储器，提高了数据读取的带宽。应用Mammousi模型数据对经过上述优化后的程序进行了测试，结果表明，
图形处理器逆时偏移程序得到了很好的优化，提高了计算效率。关键词：图形处理器；逆时偏移；随机边界；共享存储器
中图分类号：P631.3
文献标志码：A
文章编号：10008527(2012)06128905
StrategiesofAcceleratingReverseTimeMigration
UsingGraphicProcessingUnit LIU Ya-ning, LIU Guo-feng, ZHANG Zhi-fu
(School of Gesphysies aad Irformation Technology, China Unitersity of Geosciences, Bejing100083, China)
eee computing needs of its high-intensity limit its large-scale computing applications. The development of program mable Graphics Processing Unit( GPU) provides an altermative method for rapid calculation of the reverse-time migration. In this paper, we focus on the main steps using GPU to calculate the reverse time migration. We sum-eessaonaeoeooae random boundary conditions.The data transmission between the host and GPU is reduced.The second is the ap-plication of shared memory to store the forward wave field. Compared to the global memory, data reading band widths increase greatly. Finally, we use the Marmousi data to test our code, and the results show that the reverse-
time migration program running on GPU has been well optimized to improve the computational efficiency Key words: Graphic Processing Unit; reverse time migration; random boundary condition; share memory
0引言
像效果更好」，同时，逆时偏移能对回折波、多次波等波场成像，被认为是目前精度最高的地震
在逆时偏移中，地震波能够沿各个方向传播，
资料成像方法。
它比基于高频近似的Kirchhoff类方法以及假设地
逆时偏移并不是最新发展的方法，早在1978
震波只是沿垂向方向单向传播的单程波方法的成年，Hemon就提出了逆时计算的基础理论[2)。在
收稿日期：2012-0610；改回日期：2012-07-18；责任编辑：孙义梅。
基金项目：国家白热科学基金项目（51074138，41104083，41074095）：国家专项项目（SinoProbe-01，201011035）；中央高校基
本科研业务费专项资金资助项目。
作者简介：刘雅宁，女，硕士研究生，1990年出生，地球物理学专业，从事复杂介质地震波传播成像研究。
Email;sakuda728@sinacom。
通讯作者：刘国峰，男，讲师，1981年出生，地球物理学专业，从事地震波成像、反演、高性能计算等研究工作。
Email; liugfe eugh, edu. cn。