Parabolic Radon transform parallel algorithm for CPU-GPU heterogeneous platform
ZHANG Quan1,2, LIN Baiyue1, YANG Bo3, PENG Bo1, ZHANG Wei1, TU Ran1
1. School of Computer Science, Southwest Petroleum University, Chengdu, Sichuan 610500, China; 2. School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu, Sichuan 611731, China; 3. State Grid Chongqing Electric Power Research Institute, Chongqing, 404100, China
Abstract：Parabolic Radon transform is widely used for suppressing and removing multiple waves in pre-stack seismic data.Although parabolic Radon transform in mixed domain can well suppress multiples,it still needs a long time when processing seismic gathers larger and larger.In this paper, first we use the GPU to optimize the parabolic Radon transform algorithm in parallel,and use the CUDA library and other optimizing techniques to accelerate the most time-consuming Fourier transform and algebraic operations to nearly 13x speedup over serial execution in the calculating process;then based on the CPU-GPU heterogeneous platform,we propose a CPU-GPU parallel scheme,which makes full use of computer hardware,to realize the parabolic Radon transform parallel algorithm through CPU multi-threading and multiple GPUs in parallel.The speedup can be nearly 30x over serial execution.
李志娜.多次波分离与成像方法研究[D].山东青岛:中国石油大学(华东),2015.LI Zhina.Study on Methods of Multiple Separation and Multiple Imaging[D].China University of Petroleum(East China),Qingdao,Shandong,2015.
[2]
胡浩.Radon变换及其在地震资料去噪中的应用研究[D].四川成都:成都理工大学,2014.HU Hao.Research on Radon Transform and Its Application in Seismic Data Noise Elimination[D].Chengdu University of Technology,Chengdu,Sichuan,2014.
[3]
侯维丽,刘洪林,刘云鑫,等.多次波压制方法及优化组合压制技术研究[J].内陆地震,2015,29(2):122-128.HOU Weili,LIU Honglin,LIU Yunxin,et al.Study on multiple suppression methods and optimum combination suppression technique[J].Inland Earthquake,2015,29(2):122-128.
Deans S R.The Radon Transform and Some of Its Applications[M].John Wiley and Sons Inc,New York,1983.
[6]
Claerbout J F,Johnson A G.Extrapolation of time- dependent waveforms along their path of propagation[J].Geophysical Journal International,1971,26(1):285-293.
[7]
Thorson R,Claerbout J.Velocity stack and slant stack- stochastic inversion[J].Geophysics,1985,50(12):2727-2741.
[8]
Scales J,Gersztenkorn A and Treitel S.Fast lp solution of large,sparse,linear systems:Application to seismic travel time tomography[J].Journal of Computational Physics,1988,75(2):314-333.
[9]
Sacchi M,Ulrych T.High-resolution velocity gathers and offset space reconstruction[J].Geophysics, 1995, 60(4):1169-1177.
[10]
Trad D,Ulrych T and Sacchi M.Latest views of the sparse Radon transform[J].Geophysics,2003,68(1):386-399.
[11]
刘喜武,刘洪,李幼铭,等.高分辨率Radon变换方法及其在地震信号处理中的应用[J].地球物理学进展,2004,19(1):8-15.LIU Xiwu,LIU Hong,LI Youming,et al.High resolution radon transform and its application in seismic signal processing[J].Progress in Geophysics,2004,19(1):8-15.
[12]
Schonewille M A and Aaron P A.Applications of time-domain high-resolution Radon demultiple[C].SEG Technical Program Expanded Abstracts,2007,26:2565-2569.
[13]
Lu W K.An accelerated sparse time-invariant Radon transform in the mixed frequency-time domain based on iterative 2D model shrinkage[J].Geophysics,2013,78(4):V147-V155.
[14]
熊登,赵伟,张剑锋.混合域高分辨率抛物Radon变换及在衰减多次波中的应用[J].地球物理学报,2009,52(4):1068-1077.XIONG Deng,ZHAO Wei,ZHANG Jianfeng.Hybrid-domain high-resolution parabolic Radon transform and its application to demultiple[J].Chinese Journal of Geophysics,2009,52(4):1068-1077.
[15]
张军华,臧胜涛,单联瑜,等.高性能计算的发展现状及趋势[J].石油地球物理勘探,2010,45(6):918-925.ZHANG Junhua,ZANG Shengtao,SHAN Lianyu,et al.Development status and trends for high perfor-mance computing[J].Oil Geophysical Prospecting,2010,45(6):918-925.
[16]
吴吉忠,左虎.叠前衰减补偿时间偏移及GPU实现[J].石油地球物理勘探,2019,54(1):84-92.WU Jizhong,ZUO Hu.Attenuation compensation in prestack time migration and its GPU implementation[J].Oil Geophysical Prospecting,2019,54(1):84-92.
[17]
刘帅,季晓慧,芦俊,等.CPU与GPU协同并行的多分量地震数据各向异性叠前时间偏移[J].石油地球物理勘探,2019,54(1):65-72.LIU Shuai,JI Xiaohui,LU Jun,et al.Multi-component anisotropy prestack time migration based on collaborative parallel computing with CPU and GPU[J].Oil Geophysical Prospecting,2019,54(1):65-72.
[18]
赵虎,武泗海,尹成,等.基于OpenACC编程模型的逆时偏移多级并行的设计与优化[J].石油地球物理勘探,2018,53(6):1307-1313,1325.ZHAO Hu,WU Sihai,YIN Cheng,et al.Multi-level parallel design and optimization for reverse time migration based on OpenACC programming model[J].Oil Geophysical Prospecting,2018,53(6):1307-1313.
[19]
王一达,赵长海,李超,等.异构计算环境下的三维Kirchhoff叠前深度偏移混合域并行算法[J].石油地球物理勘探,2018,53(3):478-486.WANG Yida,ZHAO Changhai,LI Chao,et al.Hybrid domain parallel algorithm for 3D Kirchhoff prestack depth migration in the heterogeneous environment[J].Oil Geophysical Prospecting, 2018,53(3):478-486.
[20]
张全,林柏栎,彭博,等.基于CUDA的地震相干体并行算法[J].地质与勘探,2020,56(1):147-153.ZHANG Quan,LIN Baiyue,PENG Bo,et al.Seismic coherence parallel algorithm based on CUDA[J].Geology and Exploration,2020,56(1):147-153.
[21]
Beck A,Teboulle M.A fast iterative shrinkage-thresholding algorithm for linear inverse problems[J].SIAM Journal on Imaging Sciences,2009,2(1):183-202.
[22]
贾连凯.南海深水区多次波压制技术研究[D].山东青岛:中国科学院研究生院(海洋研究所),2014.JIA Liankai.Study of Multiple Elimination Techniques in the Deepwater Area of South China Sea[D].Graduate University of the Chinese Academy of Sciences(The Institute of Oceanology),Qingdao,Shandong,2014.
[23]
Whitehead N, Fit-Florea A.Precision & performance:Floating point and IEEE 754 compliance for NVIDIA GPUs[J].rn (A+B), 2011, 21(1):18749-19424.
