基于共散射点道集的多次波压制方法研究
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
多次波压制是地震数据处理中的重要环节。一种常用的多次波压制方法是根据一次波和多次波在成像空间的曲率差异并借助Radon变换进行分离和压制。成像空间既可以是叠前时间域也可以是叠前深度域。考虑到对复杂介质的适应性以及叠前处理对计算效率的要求,提出基于共散射点(CSP)道集的高分辨率Radon变换多次波压制方法。该方法只需要一个简单的初始速度场,就可将常规共中心点道集(CMP)映射到覆盖次数更高、炮检距覆盖范围更广的共散射点道集,然后在该道集上应用高分辨率双曲Radon变换,可较好地分离一次波和多次波。相对于以水平层状介质为假设条件的CMP道集,CSP道集更适应复杂地质构造,且时距关系满足双曲规律。模型和实际资料测试结果表明,该方法可以较好地实现速度谱能量团的聚焦,有利于较复杂地质条件下的多次波压制。
Multiple elimination is a key step in seismic data processing.One of common methods is using Radon transform to attenuate multiples based on curvature difference of primaries and multiples in the image domain.The image domain can be prestack time or depth domain after migration.Considering the capability to process complex medium and efficiency required by processing approach, we propose high-resolution Radon transform based on the CSP(common scattering point)gathers to eliminate offset range.This method maps CMP gathers to CSP gathers,which have higher folds and wider offset range only using rough velocity field.In order to better eliminate multiples,we apply sparse constraints to hyperbolic Radon transform in CSP gathers,which can provide better focus of velocity spectrum.Compared with multiple elimination method in the CMP gather,the method based on CSP gather can avoid the assumption of layered media and better process complex media even under the hyperbolic trajectory.Synthetic and real data examples demonstrate that the proposed method can better focus and separate of primary and multiple waves in the velocity spectrum.
引文
[1]Morely L.Predictive deconvolution for marine multiple suppression.Dissertation,Stanford University1982.
    [2]Yilmaz O.Seismic Data Analysis.Tulsa:Society of Exploration Geophysicists,2001.
    [3]Zhou B,Greenhalgh S A.Wave-equation extrapolation-based multiple attenuation:2-D filtering in the f-k domain.Geophysics,1994,59(9):1377-1391.
    [4]胡天跃,王润秋,White R E.地震资料处理中的聚束滤波方法.地球物理学报,2000,43(1):105-115.Hu Tianyue,Wang Runqiu,White R E.Beam forming in seismic data processing.Chinese Journal Geophysics,2000,43(1):105-115.
    [5]Verschuur D J,Berkhout A J,Wapenaar C P A.Adaptive surface-related multiple elimination.Geophysics,1992,57(9):1166-1177.
    [6]Sava P,Guitton A.Multiple attenuation in the image space.Geophysics,2005,70(1):V10-V20.
    [7]Bancroft J C.The equivalent offset method of prestack time migration.Geophysics,1998,63(6):2042-2053.
    [8]Mayne W H.Common reflection point horizontal data stacking techniques.Geophysics,1962,27(6):927-938.
    [9]Zhang Kai,Cheng Jiubing and Ma Zaitian.Pre-stack time migration and velocity analysis methods with common scatter-point gathers.Journal of Geophysics and Engineering,2006,3(3):283-289.
    [10]Stoffa P L,Buhl P,Diebold J et al.Direct mapping of seismic data to the domain of intercept time and ray parameter:A plane-wave decomposition.Geophysics,1981,46(3):255-267.
    [11]Cary P W.The simplest discrete Radon transform.SEG Technical Program Expanded Abstracts,1998,17:1999-2002.
    [12]Herrmann P,Mojesky T,Magesan M et al.Dealiased,high-resolution Radon transforms.SEG Technical Program Expanded Abstracts,2000,19:1953-1956.
    [13]Beylkin G.Discrete Radon transform:Acoustics,speech and signal processing.IEEE Transactionson,1987,35:162-172.
    [14]Sacchi M D and Ulrych T J.High-resolution velocity gathers and offset space reconstruction.Geophysics,1995,60(4):1169-1177.
    [15]Trad D T U and Sacchi M.Latest views of the sparse radon transform.Geophysics,2003,68(1):386-399.
    [16]Hu J,Fomel S,Demanet L et al.A fast butterfly algorithm for the hyperbolic Radon transform.SEG Technical Program Expanded Abstracts,2012,31:1-5.
    [17]Hampson D.Inverse velocity stacking for multiple elimination.SEG Technical Program Expanded Abstracts,1986,5:422-424.
    [18]Trad D O,Ulrych T J and Sacchi M D.Accurate interpolation with high-resolution time-variant Radon transforms.Geophysics,2002,67(2):644-656.
    [19]熊登,赵伟,张剑锋.混合域高分辨率抛物Radon变换及其在衰减多次波中的应用.地球物理学报,2009,52(4):1068-1077.Xiong Deng,Zhao Wei,Zhang Jianfeng.Hybrid-domain high-resolution parabolic Radon transform and its application to demultiple.Chinese Journal Geophysics,2009,52(4):1068-1077.
    [20]Liu Yujin,Peng Z,Symes W W et al.Sparse Radon transform with dual gradient ascent method.SEG Technical Program Expanded Abstracts,2013,32:4650-4655.
    [21]张凯.叠前偏移速度分析方法研究[博士学位论文].上海:同济大学,2008,22-27.Zhang Kai.Prestack migration velocity analysis[D].Tongji University,Shanghai,2008,22-27.
    [22]Schleicher J,Tygel M and Hubral P.3-D true-amplitude finite offset migration.Geophysics,1993,58(8):1112-1126.
    [23]Hale D.Dip move-out by Fourier transform.Geophysics,1984,49(6):741-757.
    [24]刘喜武,刘洪,李幼铭.高分辨率Radon变换方法及其在地震信号处理中的应用.地球物理学进展,2004,19(1):8-15.Liu Xiwu,Liu Hong,Li Youming.High resolution radon transform and its application in seismic signal processing.Progress in Geophysics,2004,19(1):8-15.
    [25]Ji J.CGG method for robust inversion and its application to velocity-stack inversion.Geophysics,2006,71(4):R59-R67.

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