地震叠前时间偏移的一种图形处理器提速实现方法
|
摘要
新近发展的图形处理器(GPU,Graphic Processing Unit)通用计算技术,现已日趋实用成型,并获得诸多应用领域的广泛关注.对油气勘探专项资料处理技术的运用而言,概因GPU与中央处理器(CPU)的计算性能的甚大差异,致使GPU这一通用计算技术在石油工业中的应用研究正在有效开展.本文仅借助于油气勘探中广泛使用的叠前时间偏移,旨在于扼要阐明其基于GPU应用的有效性;文中还提出一种利用GPU实现地震叠前时间偏移的软件构件方法,并针对非对称走时叠前时间偏移所拓展的应用软件提供一种具体实现架构.与以往用个人计算机(PC,Personal Computer)或者PC集群所用的叠前时间偏移相比,本文方法可甚大地提高计算效率,从而在石油物探资料处理中可显著地节约计算成本和维护费用.文中实际例证也表明,基于GPU进行高性能并行计算,当是适应目前石油工业中大规模计算需求的一个重要发展途径.
General Purpose GPU technology has being becoming mature,it has been applied in many industry area.However,due to the differences of computing feather between CPU and GPU,the study of GPU in petroleum industry application should be developed effectively.In this article,we introduce the General Purpose GPU technology and propose a method to realize pre-stack time migration software on GPU.Compared with traditional pre-stack time migration running on Personal Computer(PC)or PC-Cluster,the new programming method greatly improves computational efficiency,and then dramatically save power and maintenances cost.the actual tests in real seismic data illustrate that high performance computing based on General Purpose GPU technology(GPGPU)is an important direction of developments to meet the requirements of large scale computing in petroleum industry.
General Purpose GPU technology has being becoming mature,it has been applied in many industry area.However,due to the differences of computing feather between CPU and GPU,the study of GPU in petroleum industry application should be developed effectively.In this article,we introduce the General Purpose GPU technology and propose a method to realize pre-stack time migration software on GPU.Compared with traditional pre-stack time migration running on Personal Computer(PC)or PC-Cluster,the new programming method greatly improves computational efficiency,and then dramatically save power and maintenances cost.the actual tests in real seismic data illustrate that high performance computing based on General Purpose GPU technology(GPGPU)is an important direction of developments to meet the requirements of large scale computing in petroleum industry.
引文
[1]柳有权,刘学慧,吴恩华.基于GPU带有复杂边界的三维实时流体模拟.软件学报,2006,3:568～576Liu Y Q,Liu X H,Wu EN H.Real-time3Dfluid simulation on GPU with complex obstacles.Journal of Soft ware(in Chinese),2006,3:568～576
[2]Kenneth Moreland,Edward Angel.The FFTon a GPU.In:SIGGRAPH/Euro-graphics Workshop on Graphics Hardware2003Proceedings,2003July.112～119
[3]杨冰,蒋杰,应龙等.基于GPU的流体动力学模拟.计算机工程与应用,2007,11:7～10Yang B,Jiang J,Ying L,et al.Fluid dynamics simulation on GPU.Computer Engineering and Applications(in Chinese),2007,11:7～10
[4]Technical Brief:NVIDI A GeForce8800GPU Architecture Overview.Santa Clara,California,November2006
[5]NVIDI A CUDA Programming Guide Version1.1.Santa Clara,California,Jan,2008
[6]Luebke D,Humphreys G.How GPUs work.IEEE Computer,2007,2:96～100
[7]M.Turhan Taner,Fulton Koehler.Velocity spectra-digital computer derivation and application of velocity functions.Geophysics,1969,34:859～881
[8]刘洪,袁江华,陈景波等.大步长波场深度延拓的理论.地球物理学报,2006,49(6):1779～1793Liu H,Yuan J H,Chen J B,et al.Theory andimplementing scheme of large-step wavefield depth extrapolation.Chinese J.Geophys.(in Chinese),2006,49(6):1779～1793
[9]刘洪,刘国峰,武威等.三维波动方程逆散射的基础理论研究.石油物探,2007,46(6):569～581Liu H,Liu G F,Wu W,et al.Study on the fundamental theory of inverse scattering of multi-dimensional wave equation.Geophysical Prospecting for Petroleum(in Chinese),2007,46(6):569～581
[10]刘洪,王秀闽,曾锐等.单程波算子积分解的象征表示.地球物理学进展,2007,22(2):463～471Liu H,Wang X M,Zeng R,et al.Symbol description to integral solution of oneway wave operator.Progress in Geophysics(in Chinese),2007,22(2):463～471
[11]刘洪,刘国锋,李博等.基于横向导数的走时计算方法及其叠前时间偏移应用.石油物探,2009,48(1):Liu H,Liu G F,Li B,et al.The travel time calculation method via laterally derivative to velocity and its application in Prestack time migration.Geophysical Prospecting for Petroleum(in Chinese),2009,48(1):
[12]Lumley D,Claerbout J,Bevc D.Anti-aliased Kirchhoff3-D migration.64th Ann.Internet Mtg,Soc.Expl.Geophys.Expanded Abstracts,1994.1282～1285
[13]王秀闽.积分偏移的计算几何与并行实现[博士论文].北京:中国科学院研究生院,2008Wang X M.Computational geometry and parallel implementation of integral migration[Ph.D.thesis](in Chinese).Beijing:Graduate University of Chinese Academy of Sciences,2008
[14]张剑锋,卢宝坤,刘礼农.波动方程深度偏移的频率相关变步长延拓方法.地球物理学报,2008,51(1):221～227Zhang J F,Lu B K,Liu L N.Frequency-dependent varying-step depth extrapolation scheme for wave equation based migration.Chinese J.Geophys.(in Chinese),2008,51(1):221～227
[15]张金海,王卫民,赵连锋等.傅里叶有限差分法三维波动方程正演模拟.地球物理学报,2007,50(6):1854～1862Zhang J H,Wang W M,Zhao L F,et al.Modeling3-D scalar waves using the Fourier finite-difference method.Chinese J.Geophys.(in Chinese),2007,50(6):1854～1862
[2]Kenneth Moreland,Edward Angel.The FFTon a GPU.In:SIGGRAPH/Euro-graphics Workshop on Graphics Hardware2003Proceedings,2003July.112～119
[3]杨冰,蒋杰,应龙等.基于GPU的流体动力学模拟.计算机工程与应用,2007,11:7～10Yang B,Jiang J,Ying L,et al.Fluid dynamics simulation on GPU.Computer Engineering and Applications(in Chinese),2007,11:7～10
[4]Technical Brief:NVIDI A GeForce8800GPU Architecture Overview.Santa Clara,California,November2006
[5]NVIDI A CUDA Programming Guide Version1.1.Santa Clara,California,Jan,2008
[6]Luebke D,Humphreys G.How GPUs work.IEEE Computer,2007,2:96～100
[7]M.Turhan Taner,Fulton Koehler.Velocity spectra-digital computer derivation and application of velocity functions.Geophysics,1969,34:859～881
[8]刘洪,袁江华,陈景波等.大步长波场深度延拓的理论.地球物理学报,2006,49(6):1779～1793Liu H,Yuan J H,Chen J B,et al.Theory andimplementing scheme of large-step wavefield depth extrapolation.Chinese J.Geophys.(in Chinese),2006,49(6):1779～1793
[9]刘洪,刘国峰,武威等.三维波动方程逆散射的基础理论研究.石油物探,2007,46(6):569～581Liu H,Liu G F,Wu W,et al.Study on the fundamental theory of inverse scattering of multi-dimensional wave equation.Geophysical Prospecting for Petroleum(in Chinese),2007,46(6):569～581
[10]刘洪,王秀闽,曾锐等.单程波算子积分解的象征表示.地球物理学进展,2007,22(2):463～471Liu H,Wang X M,Zeng R,et al.Symbol description to integral solution of oneway wave operator.Progress in Geophysics(in Chinese),2007,22(2):463～471
[11]刘洪,刘国锋,李博等.基于横向导数的走时计算方法及其叠前时间偏移应用.石油物探,2009,48(1):Liu H,Liu G F,Li B,et al.The travel time calculation method via laterally derivative to velocity and its application in Prestack time migration.Geophysical Prospecting for Petroleum(in Chinese),2009,48(1):
[12]Lumley D,Claerbout J,Bevc D.Anti-aliased Kirchhoff3-D migration.64th Ann.Internet Mtg,Soc.Expl.Geophys.Expanded Abstracts,1994.1282～1285
[13]王秀闽.积分偏移的计算几何与并行实现[博士论文].北京:中国科学院研究生院,2008Wang X M.Computational geometry and parallel implementation of integral migration[Ph.D.thesis](in Chinese).Beijing:Graduate University of Chinese Academy of Sciences,2008
[14]张剑锋,卢宝坤,刘礼农.波动方程深度偏移的频率相关变步长延拓方法.地球物理学报,2008,51(1):221～227Zhang J F,Lu B K,Liu L N.Frequency-dependent varying-step depth extrapolation scheme for wave equation based migration.Chinese J.Geophys.(in Chinese),2008,51(1):221～227
[15]张金海,王卫民,赵连锋等.傅里叶有限差分法三维波动方程正演模拟.地球物理学报,2007,50(6):1854～1862Zhang J H,Wang W M,Zhao L F,et al.Modeling3-D scalar waves using the Fourier finite-difference method.Chinese J.Geophys.(in Chinese),2007,50(6):1854～1862
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心