地震叠前逆时偏移中的去噪与存储
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
地震叠前逆时偏移是当前公认的地震成像的有效途径,然而它面临着计算量甚巨,低频成像噪音以及存储量大等问题,因此,业内科研工作者对其研究乐此不疲.借助GPU/CPU协同计算可以有效解决计算量的难点,笔者已在另文中阐述,本文着重探讨成像噪音抑制以及存储问题.文中分析了叠前逆时偏移产生成像噪音的机制,据此提出在叠前地震资料中先对数据进行相位与振幅校正,进而在成像后运用拉普拉斯算子滤波法消除成像噪音,从而有效去除成像所产生的低频噪音;针对存储量,采用随机边界,用计算换存储,并借助GPU实现,节省了GPU与CPU之间的数据通讯,数值实验结果表明,采用随机边界方法的逆时偏移结果与直接存储波场的方法得到的结果差别甚小.
Pre-stack reverse time migration (RTM) is a very useful tool for seismic imaging, however, it has some problems such as highly intensive computation cost, imaging noise, and massy memory demand, which the researchers in the field have been striving to solve.The problem of intensive computation cost could be solved by GPU/CPU collaborative computing and has been discussed in the last paper and the other two problems will be emphasized in this paper. At the beginning, we analyze the generation principle of imaging noise and propose that the phase and frequency spectra of seismic data should be modified before migration and a Laplacian filter could be used to remove the low wave-number noise efficiently.Aiming at the problem of massy memory demand, we adopt the random boundary condition which reduces the memory demand but sacrifice the computation cost.The implementation could be performed on GPU and the communication between CPU and GPU could be saved which could reduce the computation cost in another way.The tests on synthetic data examples illustrate that the difference between the migration result from the method of random boundary condition and that from the traditional method of saving the footprint of the wave field could be overlooked.
引文
[1] Hemon C.Equations d'onde et modeles.Geophysical Prospecting,1978,26:790~821
    [2] Edip Baysal,Dan D Kosloff ,John W C Sherwood.Reverse time migration.Geophysics,1983,48(11) :1514~1524
    [3] McMechan G A.Migration by extrapolation of time-dependent boundary values.Geophysical Prospecting,1983,31 :413~420
    [4] Dan Loewenthal ,Irshad R Mufti.Reversed time migration in spatial frequency domain.Geophysics,1983,48(5) :627~635
    [5] Whitmore N D.Iterative depth migration by backward time propagation.53rd Annual International Meeting,SEG,Expanded Abstracts,1983. 827~830
    [6] McMechan G A,Chang W F.3D acoustic prestack reversetime migration.Geophysical Prospecting,1990,38(7) :737~756
    [7] Chang Wen-Fong,George A.McMechan.3-D elastic prestack,reverse-time depth migration.Geophysics,1994,59:597~609
    [8] Zhang Houzhu(James),Zhang Yu.Reverse time migration in 3D heterogeneous TTI media.78th Annual International Meeting,SEG Expanded Abstracts,2008,27. 2196~2200
    [9] Robin P Fletcher,Du Xiang,Paul J.FowlerReverse time migration in tilted transversely isotropic(TTI)media.Geophysics,2009,74:WCA179
    [10] Huang Tony,Zhang Yu,Zhang Houzhu,et al.Subsalt imaging using TTI reverse time migration.The Leading Edge,2009,28:448~452
    [11] Zhang Yu,Zhang Houzhu.A stable TTI reverse time migration and its implementation.79th Annual International Meeting,SEG Expanded Abstracts,2009,28:2794~2798
    [12] 邓玉琼,戴霆范,郭宗汾等.弹性波叠前有限元反时偏移.石油物探,1 990,29(3) :22~33 Teng Y C,Dai T F,John T.Kuo.Pre stack finite element reverse time migration for elastic waves.Geophysical Prospecting for Petroleum(in Chinese),1990,29(3) :22~33
    [13] 张美根,王妙月.各向异性弹性波有限元叠前逆时偏移.地球物理学报,2001,44(5) :711~719 Zhang M G,Wang M Y.Prestack finite element reverse-time migration for anisotropic elastic waves.Chinese J.Geophys.(in Chinese),2001,44(5) :711~719
    [14] 底青云,王妙月.弹性波有限元逆时偏移技术研究.地球物理学报,1997,40(4) :570~579 Di Q Y,Wang M Y.The study of finite element reverse-time migration for elastic wave.Chinese J.Geophys.(Acta Geophysica Sinica)(in Chinese),1997,40(4) :570~579
    [15] 王童奎,付兴深,朱德献等.谱元法叠前逆时偏移研究.地球物理学进展,2008,23(3) :681~685 Wang T K,Fu X S,Zhu D X,et al.Spectral-element method for prestack reverse-time migration.Progress in Geophysics(in Chinese),2008,23(3) :681~685
    [16] 陈可洋.基于高阶有限差分的波动方程叠前逆时偏移方法.石油物探,2009,48(5) :475~478 Chen K Y.High order finite difference prestack reverse time migration.Geophysical Prospecting for Petroleum(in Chinese),2009,48(5) :475~478
    [17] 王山山,聂勋碧,邓玉琼.混合法双程无反射波动方程偏移.石油地球物理勘探,1993,28(5) :537~542 Wang S S,Nie X B,Deng Y Q.Hybrid migration using two way nonreflecting wave equation.Oil Geophysical Prospecting(in Chinese),1993,28(5) :537~542
    [18] 薛东川,王尚旭.波动方程有限元叠前逆时偏移.石油地球物理勘探,2008,43(1) :17~21 Xue D C,Wang S X.Wave-equation finite-element prestack reverse time migration.Oil Geophysical Prospecting(in Chinese),2008,43(1) :17~21
    [19] 杜启振,秦童.横向各向同性介质弹性波多分量叠前逆时偏移.地球物理学报,2009,52(3) :801~807 Du Q Z,Qin T.Multicomponent prestack reverse-time migration of elastic waves in transverse isotropic medium.ChineseJ.Geophys.(in Chinese),2009,52(3) :801~807
    [20] 刘红伟,李博,刘洪等.地震叠前逆时偏移高阶有限差分算法及GPU实现.地球物理学报,2010,53(7) :1725~1733 Liu H W,Li B,Liu H,et al.The algorithm of high order finite difference pre-stack reverse time migration and GPU implementation.Chinese J.Geophys.(in Chinese),2010,53(7) :1725~1733
    [21] Baysal E,Kosloff D D,Sherwood J W C.A two-way nonreflecting wave equation.Geophysics,1984,49:132~141
    [22] Loewenthal D,Stoffa P A,Faria E L.Suppressing the unwanted reflections of the full wave equation.Geophysics,1987,52:1007~1012
    [23] Robin F Fletcher,Paul Fowler,Phil Kitchenside,et al.Suppressing artifacts in prestack reverse time migration.75th Annual International Meeting,SEG,Expanded Abstracts,2005,24(1) :2049~2051
    [24] Kwangjin Yoon,Kurt J Marfurt,William Starr.Challenges in reverse-time migration.74th Annual International Meeting,SEG Expanded Abstracts,2004,23:1057~1060
    [25] Bruno Kaelin and Antoine Guitton.Imaging condition for reverse time migration.76th Annual International Meeting,SEG,Expanded Abstracts,2006,25(1) :2594~2598
    [26] Liu Faqi,Zhang Guanquan,Scott A Morton,et al.Reverse time migration using one-way wavefield imaging condition.77th Annual International Meeting,SEG,Expanded Abstracts,2007,26(1) :2170~2174
    [27] André Bulcāo,Djalma Manoel Soares Filho and Webe Joāo Mansur.Improved quality of depth images using reverse time migration.77th Annual International Meeting,SEG,Expanded Abstracts,2007,26(1) :2407~2411
    [28] Antoine Guitton,Bruno Kaelin,and Biondo Biondi.Least square attenuation of reverse time migration artifacts.76th Annual International Meeting,SEG Expanded Abstracts,2006,25:2348~2352
    [29] Zhang Yu,James Sun.Practical issues of reverse time migration:true amplitude gathers,noise removal and harmonic-source encoding.CPS/SEG Beijing 2009 International Geophysical Conference & Exposition,2009,ID:5.
    [30] Matthew H Karazincir,Clive M Gerrard.Explicit high-order reverse time pre-stack depth migration.76th Annual International Meeting,SEG,Expanded Abstracts,2006,25(1) :2353~2357
    [31] Ruey-Chyuan Shih,Alan R Levander.Multigrid reverse-time migration.58th Annual International Meeting,SEG ,Expanded Abstracts,1988,7(1) :988~990
    [32] Christopher E Harris,George A.McMechan.Using downward continuation to reduce memory requirements in reverse-time migration.Geophysics,1992,57(6) :848~853
    [33] Guan H,Kim Y C,Ji J,et al.Multistep reverse time migration.The Leading Edge,2009,28 :442~447
    [34] William W Symes.Reverse time migration with optimal checkpointing.Geophysics,2007,72(5) :SM213-SM221
    [35] Robert G Clapp.Reverse time migration with random boundaries.79th Annual International Meeting,SEG Expanded Abstracts,2009,28:2809~2813
    [36] Isabelle Lecomte.Resolution and illumination analyses in PSDM:A ray-based approach.The Leading Edge,2008,27(5) :650~663
    [37] Eric Dussaud,William W Symes,Paul Williamson,et al.Computational strategies for reverse-time migration.78th Annual International Meeting,SEG Expanded Abstracts,2008,27:2267~2271

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