基于GPU加速策略的稀疏反演表面多次波压制方法
|
摘要
基于稀疏反演的表面多次波压制方法(EPSI)能够避免传统表面多次波压制方法对有效波的损伤,但是计算成本较高。研究基于新GPU加速策略的稀疏反演表面多次波压制方法,将所有频率成分整体导入设备端,并将中间结果存储于设备端,实现所有频率成分的同步运算,并减少主机端和设备端的数据传输时间。SMAART模型测试结果表明:基于新GPU加速策略的EPSI方法的计算效率较常规EPSI方法的提升11倍,较基于传统GPU加速策略的EPSI方法的提升4倍。基于新的GPU加速策略的稀疏反演表面多次波压制方法正确有效,计算效率更高。
Suppression of surface-related multiple by sparse inversion can avoid the damage to primary in traditional multiple suppression method.However,the computational cost of this method is high.In order to reduce the computational cost of the method,a sparse inversion surface multiple suppression method based on GPU acceleration was studied and a new GPU acceleration strategy was adopted.All frequency components were integrated into the device and the intermediate results were stored on the device.It realizes the synchronous operation of all frequency components and reduces the data transmission time between the host and the device.According to the SMAART model test,the computing efficiency of the EPSI method based on the new GPU acceleration strategy is 11 times higher than that of the conventional EPSI method,and the computation efficiency of the EPSI method based on the traditional GPU acceleration strategy is increased by 4 times.The sparse inversion surface multiple suppression method based on the new GPU acceleration technology is correct and effective,and the calculation efficiency is higher.
Suppression of surface-related multiple by sparse inversion can avoid the damage to primary in traditional multiple suppression method.However,the computational cost of this method is high.In order to reduce the computational cost of the method,a sparse inversion surface multiple suppression method based on GPU acceleration was studied and a new GPU acceleration strategy was adopted.All frequency components were integrated into the device and the intermediate results were stored on the device.It realizes the synchronous operation of all frequency components and reduces the data transmission time between the host and the device.According to the SMAART model test,the computing efficiency of the EPSI method based on the new GPU acceleration strategy is 11 times higher than that of the conventional EPSI method,and the computation efficiency of the EPSI method based on the traditional GPU acceleration strategy is increased by 4 times.The sparse inversion surface multiple suppression method based on the new GPU acceleration technology is correct and effective,and the calculation efficiency is higher.
引文
[1]Weglein A B.Multiple attenuation:An overview of recent advances and the road ahead[J].The Leading Edge,1999,18(1):40-44.
[2]Robison E A.Predictive decomposition of time series with applications to seismic exploration[D].Cambridge:MIT,2009.
[3]Foster D J,Mosher C C.Suppression of multiple reflections using the Radon transform[J].Geophysics,1992,57(3):386-395.
[4]Kabir M M N,Marfurt K J.Toward true amplitude multiple removal[J].The Leading Edge,1999,18(1):66-73.
[5]王维红,张振,石颖,等.λ-f域抛物Radon变换多次波压制方法[J].东北石油大学学报,2015,39(1):17-22.Wang Weihong,Zhang Zhen,Shi Ying,et al.Investigation of multiple suppression by parabolic Radon transform inλ-fdomain[J].Journal of Northeast Petroleum University,2015,39(1):17-22.
[6]Taner M T.Long period sea floor multiples and their suppression[J].Geophysics Prospecting,1980,28(1):30-48.
[7]Berryhill J R,Kim Y C.Deep-water peg legs and multiples:Emulation and suppression[J].Geophysics,1986,51(12):2177-2184.
[8]Wiggins J W.Attenuation of complex water-bottom multiples by wave-equation based prediction and subtraction[J].Geophysics,1988,53(12):1527-1539.
[9]Berkhout A J.Seismic migration,imaging of acoustic energy by wave field extrapolation,A:Theoretical aspects[M].Amsterdam:Elsevier,1982.
[10]Verschuur D J.Surface-related multiple elimination,an inversion approach[D].Delft:Delft University of Technology,1991.
[11]Berkhout A J,Verschuur D J.Estimation of multiples scattering by iterative inversion,part I:Theoretical considerations[J].Geophysics,1997,62(5):1586-1595.
[12]Weglein A B,Gasparotto F A,Garvalho P M,et al.An inverse-scatting series method for attenuating multiples in seismic reflection data[J].Geophysics,1997,62(6):1975-1989.
[13]李鹏.复杂介质多次波处理方法研究[D].北京:中国科学院地质与地球物理研究所,2007.Li Peng.The Study on surface related multiple elimination[D].Beijing:Institute of Geology and Geophysics,2007.
[14]Groenestijn G J A V,Verschuur D J.Estimating primaries by sparse inversion and application to near-offset data reconstruction[J].Geophysics,2009,74(3):A23-A28.
[15]Vanborselen R G,Fokkema J T,Vandenberg P M.Removal of surface-related wave phenomena:The marine case[J].Geophysics,1996,61(1):202-210.
[16]Biersteker J.MAGIC:Shell's surface multiple attenuation technique[J].Seg Technical Program Expanded Abstracts,2001,20(1):1301-1304.
[17]Amundsen L.Elimination of free-surface related multiples without need of the source wavelet[J].Geophysics,2001,66(1):327-341.
[18]Wang Y.Multiple prediction through inversion:A fully data-driven concept for surface-related multiple attenuation[J].Geophysics,2004,69(2):547-553.
[19]Berkhout A J.Seismic processing in the inverse data space[J].Geophysics,2006,71(4):A29-A33.
[20]马继涛.SRME压制多次波方法研究[D].北京:中国石油大学(北京),2009.Ma Jitao.Study on inultiple suppression method with SRME[D].Beijing:China University of Petroleum(Beijing),2009.
[21]Herrmann F J,B9niger U,Verschuur D J.Non-linear primary multiple separation with directional curvelet frames[J].Geophysical Journal International,2007,170(2):781-799.
[22]宋家文,Verschuur D J,陈小宏.多次波压制的研究现状与进展[J].地球物理学进展,2014,29(1):240-247.Song Jiawen,Verschuur D J,Chen Xiaohong.Research status and progress in multiple elimination[J].Progress in Geophysics,2014,29(1):240-247.
[23]李博,刘国峰,刘洪.地震叠前时间偏移的一种图形处理器提速实现方法[J].地球物理学报,2009,52(1):245-252.Li Bo,Liu Guofeng,Liu Hong.A method of using GPU to accelerate seismic pre-stack time migration[J].Chinese Journal of Geophysics,2009,52(1):245-252.
[24]石颖,刘洪.基于GPU的表面多次波预测技术[J].石油地球物理勘探,2010,45(4):540-544.Shi Ying,Liu Hong.GPU-based surface multiple prediction technique[J].Oil Geophysical Prospecting,2010,45(4):540-544.
[25]王天野,白军辉,王维红,等.基于反泄漏Fourier变换的地震数据重建方法[J].东北石油大学学报,2016,40(5):102-107.Wang Tianye,Bai Junhui,Wang Weihong,et al.Seismic data reconstruction method based on anti leakage Fourier transform[J].Journal of Northeast Petroleum University,2016,40(5):102-107.
[26]宋家文,赵波,李合群.应用稀疏反演方法压制自由表面表次波[J].石油地球物理勘探,2015,50(4):613-618.Song Jiawen,Zhao Bo,Li Hequn.Surface-related multiple elimination with sparse inversion[J].Oil Geophysical Prospecting,2015,50(4):613-618.
[27]Cheng J,Grossman M,Mckercher T.Professional CUDA C programming[M].[S.l.]:Wrox Press Ltd,2014.
[2]Robison E A.Predictive decomposition of time series with applications to seismic exploration[D].Cambridge:MIT,2009.
[3]Foster D J,Mosher C C.Suppression of multiple reflections using the Radon transform[J].Geophysics,1992,57(3):386-395.
[4]Kabir M M N,Marfurt K J.Toward true amplitude multiple removal[J].The Leading Edge,1999,18(1):66-73.
[5]王维红,张振,石颖,等.λ-f域抛物Radon变换多次波压制方法[J].东北石油大学学报,2015,39(1):17-22.Wang Weihong,Zhang Zhen,Shi Ying,et al.Investigation of multiple suppression by parabolic Radon transform inλ-fdomain[J].Journal of Northeast Petroleum University,2015,39(1):17-22.
[6]Taner M T.Long period sea floor multiples and their suppression[J].Geophysics Prospecting,1980,28(1):30-48.
[7]Berryhill J R,Kim Y C.Deep-water peg legs and multiples:Emulation and suppression[J].Geophysics,1986,51(12):2177-2184.
[8]Wiggins J W.Attenuation of complex water-bottom multiples by wave-equation based prediction and subtraction[J].Geophysics,1988,53(12):1527-1539.
[9]Berkhout A J.Seismic migration,imaging of acoustic energy by wave field extrapolation,A:Theoretical aspects[M].Amsterdam:Elsevier,1982.
[10]Verschuur D J.Surface-related multiple elimination,an inversion approach[D].Delft:Delft University of Technology,1991.
[11]Berkhout A J,Verschuur D J.Estimation of multiples scattering by iterative inversion,part I:Theoretical considerations[J].Geophysics,1997,62(5):1586-1595.
[12]Weglein A B,Gasparotto F A,Garvalho P M,et al.An inverse-scatting series method for attenuating multiples in seismic reflection data[J].Geophysics,1997,62(6):1975-1989.
[13]李鹏.复杂介质多次波处理方法研究[D].北京:中国科学院地质与地球物理研究所,2007.Li Peng.The Study on surface related multiple elimination[D].Beijing:Institute of Geology and Geophysics,2007.
[14]Groenestijn G J A V,Verschuur D J.Estimating primaries by sparse inversion and application to near-offset data reconstruction[J].Geophysics,2009,74(3):A23-A28.
[15]Vanborselen R G,Fokkema J T,Vandenberg P M.Removal of surface-related wave phenomena:The marine case[J].Geophysics,1996,61(1):202-210.
[16]Biersteker J.MAGIC:Shell's surface multiple attenuation technique[J].Seg Technical Program Expanded Abstracts,2001,20(1):1301-1304.
[17]Amundsen L.Elimination of free-surface related multiples without need of the source wavelet[J].Geophysics,2001,66(1):327-341.
[18]Wang Y.Multiple prediction through inversion:A fully data-driven concept for surface-related multiple attenuation[J].Geophysics,2004,69(2):547-553.
[19]Berkhout A J.Seismic processing in the inverse data space[J].Geophysics,2006,71(4):A29-A33.
[20]马继涛.SRME压制多次波方法研究[D].北京:中国石油大学(北京),2009.Ma Jitao.Study on inultiple suppression method with SRME[D].Beijing:China University of Petroleum(Beijing),2009.
[21]Herrmann F J,B9niger U,Verschuur D J.Non-linear primary multiple separation with directional curvelet frames[J].Geophysical Journal International,2007,170(2):781-799.
[22]宋家文,Verschuur D J,陈小宏.多次波压制的研究现状与进展[J].地球物理学进展,2014,29(1):240-247.Song Jiawen,Verschuur D J,Chen Xiaohong.Research status and progress in multiple elimination[J].Progress in Geophysics,2014,29(1):240-247.
[23]李博,刘国峰,刘洪.地震叠前时间偏移的一种图形处理器提速实现方法[J].地球物理学报,2009,52(1):245-252.Li Bo,Liu Guofeng,Liu Hong.A method of using GPU to accelerate seismic pre-stack time migration[J].Chinese Journal of Geophysics,2009,52(1):245-252.
[24]石颖,刘洪.基于GPU的表面多次波预测技术[J].石油地球物理勘探,2010,45(4):540-544.Shi Ying,Liu Hong.GPU-based surface multiple prediction technique[J].Oil Geophysical Prospecting,2010,45(4):540-544.
[25]王天野,白军辉,王维红,等.基于反泄漏Fourier变换的地震数据重建方法[J].东北石油大学学报,2016,40(5):102-107.Wang Tianye,Bai Junhui,Wang Weihong,et al.Seismic data reconstruction method based on anti leakage Fourier transform[J].Journal of Northeast Petroleum University,2016,40(5):102-107.
[26]宋家文,赵波,李合群.应用稀疏反演方法压制自由表面表次波[J].石油地球物理勘探,2015,50(4):613-618.Song Jiawen,Zhao Bo,Li Hequn.Surface-related multiple elimination with sparse inversion[J].Oil Geophysical Prospecting,2015,50(4):613-618.
[27]Cheng J,Grossman M,Mckercher T.Professional CUDA C programming[M].[S.l.]:Wrox Press Ltd,2014.