基于波动方程三维表面多次波预测方法研究
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
与传统的二维表面多次波预测算法相比,基于波动方程的全三维表面多次波预测方法无需对地下介质做简单近似,其更符合地震波在地下介质中传播的真实状况,是地震资料处理中解决多次波预测问题的强有力工具.本文从三维多次波预测的基本理论出发,给出了全三维多次波预测算法的预测矩阵表示、计算方法以及实现条件,采用GPU(图形处理器)加速全三维表面多次波预测,较传统的CPU串行计算,GPU并行预测表面多次波的计算效率约提高165倍.文中分别利用二维和三维表面多次波预测算法对理论模拟的含表面多次波的三维地震数据进行多次波预测计算,对比分析结果表明,相比于二维算法,文中所述的基于波动方程的全三维表面多次波预测效果明显改善,其计算精度更高,辅以合理有效的自适应相减算法,可获得高精度的地震勘探资料表面多次波压制数据.
Compare with the traditional 2D surface-related multiple prediction algorithm, 3D surface-related multiple prediction approach based on wave equation needs no assumption for underground medium. It is a powerful tool to deal with multiple prediction problem in seismic data processing, which is more conform to the actual situation that seismic wave propagate through the underground medium. Based on the fundamentals of multiple prediction, the matrix representation, calculation approach and realization condition for full 3D surface-related multiple prediction are proposed in this paper. Also using GPU (Graphic Processing Unit) to accelerate full 3D surface-related multiple prediction algorithm, compare to traditional CPU computation, the speed is 165 times faster. Multiple prediction tests are done on synthetic 3D seismic data with surface-related multiples, 2D and 3D surface-related multiple prediction effect from 2D and 3D algorithms are compared and analyzed, and the result demonstrates that full 3D surface-related multiple prediction approach based on wave equation has higher computation accuracy and better prediction effect than 2D algorithm. With the help of reasonable adaptive subtraction, we can obtain seismic prospecting data with high precision after surface-related multiple suppression.
引文
[1]郭书娟,李振春,仝兆岐等.表层多次波成像方法技术研究.地球物理学进展,2012,27(6):2570-2576.Guo S J,Li Z C,Tong Z Q,et al.Method and technique forimaging of surface-related multiples.Progress in Geophysics(in Chinese),2012,27(6):2570-2576.
    [2]赵昌垒,叶月明,姚根顺等.线性拉东域预测反褶积在海洋多次波去除中的应用.地球物理学进展,2013,28(2):1026-1032.Zhao C L,Ye Y M,Yao G S,et al.Prediction deconvolutionin linear radon domain on the application of ocean multiplesattenuation.Progress in Geophys.(in Chinese),2013,28(2):1026-1032.
    [3]胡天跃.地震资料叠前去噪技术的现状与未来.地球物理学进展,2002,17(2):218-223.Hu T Y.The current situation and future of seismic dataprestack noise attenuation techniques.Progress in Geophys.(in Chinese),2002,17(2):218-223.
    [4]王维红,崔宝文,刘洪.表面多次波衰减的研究现状与进展.地球物理学进展,2007,22(1):156-164.Wang W H,Cui B W,Liu H.Research progress in surface-related multiple attenuation.Progress in Geophys.(in Chinese),2007,22(1):156-164.
    [5]石颖,邢小林.表面多次波压制的研究进展:回顾与展望.地球物理学进展,2011,26(6):2046-2054.Shi Y,Xing X L.Investigation progress on surface-relatedmultiple suppression:review and outlook.Progress inGeophys.(in Chinese),2011,26(6):2046-2054.
    [6]吴迪,杨长春.数据驱动型多次波衰减方法的研究.地球物理学进展,2008,23(1):98-103.Wu D,Yang C C.Research on a fully-data driven multipleattenuation approach.Progress in Geophys.(in Chinese),2008,23(1):98-103.
    [7]赵保宗,孙永清,李学聪.基于波动方程的多次波压制方法应用研究.地球物理学进展,2010,25(1):272-281.Zhao B,Sun Y Q,Li X C.The progress of the multipleattenuation method based on the wave-equation.Progress inGeophys.(in Chinese),2010,25(1):272-281.
    [8]李鹏,刘伊克,常旭等.多次波问题的研究进展.地球物理学进展,2006,21(3):888-897.Li P,Liu Y K,Chang X,et al.Progress on the multipleproblems.Progress in Geophys.(in Chinese),2006,21(3):888-897.
    [9]陈小宏,刘华锋.预测多次波的逆散射级数方法与SRME方法及比较.地球物理学进展,2012,27(3):1040-1050.Chen X H,Liu H F.Comparison between inverse scatteringseries method and SRME method in free surface relatedmultiple prediction.Progress in Geophys.(in Chinese),2012,27(3):1040-1050.
    [10]石颖,井洪亮,李莹.反馈迭代法压制表面多次波效果分析.地球物理学进展,2012,27(4):1493-1500.Shi Y,Jing H L,Li Y.Surface-related multiple suppressioneffect analysis by feedback iteration approach.Progress inGeophys.(in Chinese),2012,27(4):1493-1500.
    [11]肖二莲,陈瑜,万欢等.SRME多次波衰减方法在海洋地震资料中的应用.地球物理学进展,2010,25(3):1057-1064.Xiao E L,Chen Y,Wan H,et al.Surface-related multipleelimination on marine seismic data.Progress in Geophys.(inChinese),2010,25(3):1057-1064.
    [12]Donno D,Chauris H,Noble M.Curvelet-based multipleprediction.Geophysics,2010,75(6):WB255-WB263.
    [13]石颖,王建民,井洪亮等.多道自适应匹配滤波方法压制表面多次波.地球物理学进展,2013,28(2):785-792.Shi Y,Wang J M,Jing H L,et al.Suppressing surface-related multiple by multi-trace adaptive matching filterapproach.Progress in Geophys.(in Chinese),2013,28(2):785-792.
    [14]Berkhout A J.Seismic Migration,Practical aspects Pt.B:Imaging of Acoustic Energy by Wave Field Extrapolation.Amsterdam:Elsevier Science Publishing Company,Inc,1982.20-43.
    [15]Verschuur D J,Berkhout A J,Wapenaar C P A.Adaptivesurface-related multiple elimination.Geophysics,1992,57(9):1166-1177.
    [16]Berkhout A J,Verschuur D J.Estimation of multiplescattering by iterative inversion,partⅠ:Theoreticalconsiderations.Geophysics,1997,62(5):1586-1595.
    [17]Verschuur D J,Berkhout A J.Estimation of multiplescattering by iterative inversion,partⅡ:Practical aspectsand examples.Geophysics,1997,62(5):1596-1611.
    [18]Kelamis P G,Verschuur D J.Surface-related multipleelimination on land seismic data—Strategies via case studies.Geophysics,2000,65(3):719-734.
    [19]Wang Y H.Multiple prediction through inversion:a fullydata-driven concept for surface-related multiple attenuation.Geophysics,2004,69(2):547-553.
    [20]Wang Y H.Multiple prediction through inversion:Theoreticaladvancements and real data application.Geophysics,2007,72(2):V33-V39.
    [21]Van Groenestjin G J A,Verschuur D J.Estimating primariesby sparse inversion and application to near-offset datareconstruction.Geophysics,2009,74(3):A23-A28.
    [22]Van Groenestjin G J A,Verschuur D J.Estimation ofprimaries and near-offset reconstruction by sparse inversion:Marine data applications.Geophysics,2009,74(6):R119-R128.
    [23]Lin D C,Young J,Lin W J,et al.3DSRME prediction andsubtraction practice for better imaging.75th Ann.Internat.Mtg.,Soc.of Expi.Geophys.,Expanded Abstracts,2005.2088-2091.
    [24]Van Dedem E J,Verschuur D J.Analysis of surface-relatedmultiples in 3-D media.67th Annual International Meeting,SEG,Expanded Abstracts,1997.1180-1183.
    [25]Levin S A.Prestack poststack 3D multiple prediction.72ndAnnual International Meeting,SEG,Expanded Abstracts,2002.2110-2113.
    [26]Van Dedem E J,Verschuur D J.3Dsurface-related multipleprediction using sparse inversion:Experience with field data.72nd Annual International Meeting,SEG,Expanded Abstracts,2002.2094-2097.
    [27]Kleemeyer G,Pettersson S E,Eppenga R,et al.It’s amagic:Industry first 3Dsurface multiple elimination and pre-stack depth migration on Ormen Lange.65th AnnualInternational Meeting,EAGE,Extended Abstracts,2003.
    [28]Lin D C,Young J,Huang Y,et al.3DSRME application inthe Gulf of Mexico.74th Annual International Meeting,SEG,Expanded Abstracts,2004.1257-1260.
    [29]Matson K H,Michell S,Abma R,et al.Advanced subsaltimaging and 3Dsurface multiple attenuation in Atlantis:Acase study.74th Annual International Meeting,SEG,Expanded Abstracts,2004.1269-1272.
    [30]Moore I,Dragoset W H.Practical,3D surface-relatedmultiple prediction(SMP).74th Annual International Meeting,SEG,Expanded Abstracts,2004.1249-1252.
    [31]Baumstein A,Hadidi M T.3D surface-related multipleelimination:Data reconstruction and application to field data.Geophysics,2006,71(3):E25-E33.
    [32]Dragoset B,Verschuur E,Moore L,et al.A perspective on3Dsurface-related multiple elimination.Geophysics,2010,75(5):75A245-75A261.
    [33]Zhang J H,Wang S Q,Yao Z X.Accelerating 3DFouriermigration with graphics processing units.Geophysics,2009,74(6):WCA129-WCA139.
    [34]Clapp R G,Fu H H,Lindtjorn O.Selecting the righthardware for reverse time migration.The Leading Edge,2010,29(1):48-58.
    [35]李博,刘国峰,刘洪.地震叠前时间偏移的一种图形处理器提速实现方法.地球物理学报,2009,52(1):245-252.Li B,Liu G F,Liu H.A method of using GPU to accelerateseismic pre-stack time migration.Chinese J.Geophys.(inChinese),2009,52(1):245-252.
    [36]刘国峰,刘洪,李博等.山地地震资料叠前时间偏移方法及其GPU实现.地球物理学报,2009,52(12):3101-3108.Liu G F,Liu H,Li B,et al.Method of prestack timemigration of seismic data of mountainous regions and its GPUimplementation.Chinese J.Geophys.(in Chinese),2009,52(12):3101-3108.
    [37]刘红伟,李博,刘洪等.地震叠前逆时偏移高阶有限差分算法及GPU实现.地球物理学报,2010,53(7):1725-1733.Liu H W,Li B,Liu H,et al.The algorithm of high orderfinite difference pre-stack reverse time migration and GPUimplemenation.Chinese J.Geophys.(in Chinese),2010,53(7):1725-1733.
    [38]Baumstein A,Hadidi M T,Hinkley D L,et al.A practicalprocedure for application of 3DSRME to conventional marinedata.The Leading Edge,2005,24(3):254-258.
    [39]石颖,刘洪,邹振.基于波动方程表面多次波预测与自适应相减方法研究.地球物理学报,2010,53(7):1716-1724.Shi Y,Liu H,Zou Z.Surface-related multiples predictionbased on wave equation and adaptive subtraction investigation.Chinese J.Geophys.(in Chinese),2010,53(7):1716-1724.
    [40]Xu S,Zhang Y,LambaréG.Antileakage Fourier transformfor seismic data regularization in higher dimensions.Geophysics,2010,75(6):WB113-WB120.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心