• 中国计算机学会会刊
  • 中国科技核心期刊
  • 中文核心期刊

J4 ›› 2015, Vol. 37 ›› Issue (05): 907-913.

• 论文 • 上一篇    下一篇

基于Xeon Phi平台的波动方程叠前深度偏移

杨祥森1,金君2,王鹏1,马召贵1,亢永敢1   

  1. (1.中国石化石油物探技术研究院,江苏 南京 211113;2.Intel亚太研发有限公司,上海 200241)
  • 收稿日期:2014-08-16 修回日期:2014-10-20 出版日期:2015-05-25 发布日期:2015-05-25

Wave equation prestack depth migration
based on Xeon Phi platform  

YANG Xiangsen1,JIN Jun2,WANG Peng1,MA Zhaogui1,KANG Yonggan1   

  1. (1.Sinopec Geophysical Research Institute,Nanjing 211113;
    2.Intel AsiaPacific Research and Development Ltd,Shanghai 200241,China)
  • Received:2014-08-16 Revised:2014-10-20 Online:2015-05-25 Published:2015-05-25

摘要:

波动方程叠前深度偏移适用于强横向变速介质,是一种高精度成像方法,但其巨大的计算量阻碍了该技术的应用。Xeon Phi是一种全新的高性能计算设备,为波动方程叠前深度偏移方法的推广应用提供了新的技术支持。以裂步傅里叶算子为例,介绍了面向Xeon Phi平台的偏移算法移植和优化方法,即采用offload模式将计算核函数加载到Xeon Phi设备上,在Xeon Phi 协处理器上采用多线程方式,并且调整程序结构,充分利用SIMD矢量引擎提高向量化处理效率。扩展负载动态均衡的并行框架,形成了一套适用于大规模异构系统、基于Xeon Phi平台的波动方程叠前深度偏移软件。实际数据测试表明Xeon Phi平台可以极大地提高地震偏移处理效率,具有良好的可扩展性。

关键词: Xeon Phi平台;波动方程叠前深度偏移;裂步傅里叶算子;并行框架;offload 模式

Abstract:

Wave equation prestack depth migration is a high-precision seismic imaging method,suitable for complex media with sharp lateral velocity variation,however,its huge computation hinders its application.Xeon Phi is a new product for throughput computing which provides a better choice for these methods.Taking the Split Step Fourier algorithm for example, we introduce the code migration and performance optimization on Xeon Phi platform.We offload the kernel which is fully threaded to Xeon Phi by predefined compiler pragmas so as to take full advantage of SIMD engine to improve the efficiency of the parallelization processing; and we change the code structure to ensure the vectorization.Besides,by extending the parallel processing framework with the dynamic workloadbalance mechanism,we develop a seismic imaging system based on Xeon Phi platform which is suitable for largescale,full hybrid heterogeneous clusters.Finally, we conduct a real workload test and the results show that Xeon Phi platform can significantly improve the seismic imaging efficiency,and has good core scalability.

Key words: Xeon Phi platform;wave equation prestack depth migration;split step Fourier algorithm;parallel processing framework;offload pattern