时-空局域化地震波传播方法:Dreamlet叠前深度偏移
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度.
We propose a complete time-space localized seismic wavefield decomposition and propagation method─dreamlet migration.The time-space wavefield is projected to the time-frequency-space-wavenumber domain by the local cosine transform,and the wavefield is extrapolated to next depth by the local phase space dreamlet one-way propagator.The wavefield and propagator after decomposition are very sparse and most importantly,they stay in the real data domain,which means the process of propagation and imaging complete in the real data domain.In this paper,we derive the formulas of the wavefield decomposition and phase-space dreamlet one-way propagator.The dreamlet domain wave propagation is sparse matrix multiplication,and the valid time length of the receiver wavefield is shrinking as the depth increasing.Numerical examples and tests on the synthetic two dimensional SEG/EAGE salt model and SIGSBEE model demonstrate the validity of this method.Migration results show that this method has good imaging quality in laterally high contrast heterogeneous media.
引文
[1]Steinberg B Z.Evolution of local spectra in smoothly varyingnonhomogeneous environments local canonization and marchingalgorithms.J.Acoust.Soc.Am.,1993,93(5):2566-2580.
    [2]Steinberg B Z,Birman R.Phase-space marching algorithm inthe presence of a planar wave velocity discontinuity-Aqualitative study.J.Acoust.Soc.Am.,1995,98(1):484-494.
    [3]Wu R S,Wang Y Z,Gao J H.Beamlet migration based onlocal perturbation theory.∥70th Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,2000:1008-1011.
    [4]Wu R S,Chen L.Wave propagation and imaging usingGabor-Daubechies beamlets.∥Theory and ComputationalAcoustics,World Scientific,New Jersey,2002:661-670.
    [5]Chen L,Wu R S,Chen Y.Target-oriented beamletmigration based on Gabor-Daubechies frame decomposition.Geophysics,2006,71(2):s37-s52.
    [6]Wu R S,Wang Y Z,Luo M Q.Beamlet migration using localcosine basis.Geophysics,2008,73(5):207-217.
    [7]毛剑,吴如山,高静怀.利用局部谐和基小波束的高精度叠前深度域偏移成像方法研究.地球物理学报,2010,53(10):2442-2451.Mao J,Wu R S,Gao J H.High-accuracy prestack depthmigration and imaging using local harmonic basis beamlet.Chinese J.Geophys.(in Chinese),2010,53(10):2442-2451.
    [8]高静怀,周艳辉,毛剑等.局部角度域波传播步进算法研究.地球物理学报,2007,50(1):248-259.Gao J H,Zhou Y H,Mao J,et al.A wave propagationmethod in local angle domain.Chinese J.Geophys.(inChinese),2007,50(1):248-259.
    [9]Ma Y W,Margrave G F.Seismic depth migration with theGabor transform.Geophysics,2008,73(3):s91-s97.
    [10]Pedersen,Brandsberg-Dahl S,Ursin B.Seismic imagingusing lateral adaptive windows.Geophysics,2010,75(2):s73-s79.
    [11]Candès E J,Demanet L,Donoho D,Li Y.Fast discretecurvelet transforms.SIAM Multiscale Modeling and Simulation,2006,5(3):861-899.
    [12]Herrmann F.Optimal seismic imaging with curvelets.∥73th Ann.Internat.Mtg.,Soc.Expl.Geophys.,ExpandedAbstracts,2003:997-1000.
    [13]Candès E J,Demanet L.The curvelet representation of wavepropagators is optimally sparse.Communications on Pure andApplied Mathematics,2004,58:1472-1528.
    [14]Douma H,de Hoop M V.Leading-order seismic imagingusing curvelets.Geophysics,2007,72(6):s231-s248.
    [15]Wu R S,Wu B Y,Geng Y.Seismic wave propagation andimaging using time-space wavelets.∥78th Ann.Internat.Mtg.,Soc.Expl.Geophys.,Expanded Abstracts,2008:2983-2987.
    [16]Wu R S,Wu B Y,Geng Y.Imaging in compressed domainusing dreamlets.∥CPS/SEG Beijing′2009,InternationalGeophysical Conference,Expanded Abstracts,ID:57.
    [17]Mao J,Wu R S,Gao J H.Directional illumination analysisusing the local exponential frame.Geophysics,2010,75(4):s163-s174.
    [18]Püschel M,Moura J.The algebraic approach to the discretecosine and sine transform and their fast algorithms.SIAMJournal of Computing,2003,32(5):1280-1316.
    [19]Coifman R R,Meyer Y.Remaxques sur 1′analyse de Fouriera fengtre.Comptes Rendus de l′Acad6 mie des Sciences,1991,312:259-261.

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