伪深度域声波数值模拟方法及应用
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
地震勘探所面对的地下介质中常包括低速体、起伏地表、高陡构造、孔缝单元等,对这些区域进行地震波场模拟时,为保证模拟精度和计算稳定性,网格间距需要取得很小,从而导致基于笛卡尔坐标系的常规有限差分正演模拟方法的局部过采样问题。为此,本文在前人研究成果的基础上,将伪深度思想引入到正演模拟中,利用曲线坐标系下的梯度、散度公式推导出伪深度域一阶速度—应力方程,实现了伪深度域地震正演模拟。同时,鉴于横向过采样问题没有得到有效解决的实际情况,采用自适应变差分算子计算横向空间导数,提出了伪深度域自适应变差分算子算法,综合考虑了纵、横向的空间采样问题,进一步提高了计算效率。最后,将算子扩展到逆时偏移中,实现了伪深度域自适应逆时偏移算法。模型试算结果表明:相比常规正演模拟算法,伪深度域自适应变差分算子算法在保证精度的前提下,计算效率提高了约25%,存储占用降低了近30%。
Surface seismic exploration usually meets challenges like low-velocity bodies,surface topography,steep structures,fractures and so on.Seismic wavefield simulation with conventional Cartesian finite difference forward modeling method in these areas often needs very small grid spacing to ensure the accuracy and stability.So it leads local oversampling problem.To overcome this problem,we introduce the idea of pseudo-depth to the forward modeling.Using the gradient and divergence formula in curvilinear coordinates,we derive the first-order velocity-stress equation in pseudo-depth domain and implement the pseudo-depth domain seismic modeling algorithm.Meanwhile,considering uneven horizontal sampling,we introduce adaptive variable-length spatial operators to calculate the horizontal spatial derivatives and propose seismic forward modeling method with adaptive variable-length spatial operators in the pseudodepth domain.Finally,by extending the forward operator to reverse time migration,we implement pseudo-depth domain reverse time migration with adaptive variable-length spatial operators.Model test results show that our pseudo-depth domain method with adaptive variable-length spatial operators can reduce about 25%of computer time consuming and nearly 30% of the storage without influencing the accuracy of the results.
引文
[1]董良国.地震波数值模拟与反演中几个关键问题研究[博士学位论文].上海:同济大学,2003.
    [2]裴正林,牟永光.地震波传播数值模拟.地球物理学进展,2004,19(4):933-941.Pei Zhenglin,Mu Yongguang.Numerical simulation of seismic wave propagation.Progress in Geophysics,2004,19(4):933-941.
    [3]张剑锋.弹性波数值模拟的非规则网格差分法.地球物理学报,1998,41(增刊):357-366.Zhang Jianfeng.Non-orthogonal grid finite-difference method for numerical simulation of elastic wave propagation.Chinese Journal of Geophysic,1998,41(Sl):357-366.
    [4]Mcozo P.Finite-difference technique for SH waves in2-D media using irregular grids-application to the seismic response problem.Geophys J Int,1989,99(2):321-329.
    [5]Kristek J,Moczo P.Seismic wave propagation in viscoelastic media with material discontinuities-a 3D4thorder staggered-grid finite-difference modeling.Bull Seism Soc Am,2003,93(5):2273-2280.
    [6]朱生旺,曲寿利,魏修成等.变网格有限差分弹性波方程数值模拟方法.石油地球物理勘探,2007,42(6):634-639.Zhu Shengwang,Qu Shouli,Wei Xiucheng et al.Numeric simulation by grid-various finite-difference elastic wave equation.OGP,2007,42(6):634-639.
    [7]李振春,张慧,张华.一阶弹性波方程的变网格高阶有限差分数值模拟.石油地球物理勘探,2008,43(6):711-716.Li Zhenchun,Zhang Hui,Zhang Hua.Variable-grid high-order finite-difference numeric simulation of first-order elastic wave equation.OGP,2008,43(6):711-716.
    [8]黄超,董良国.可变网格与局部时间步长的交错网格高阶差分弹性波模拟.地球物理学报,2009,52(11):2870-2878.Huang Chao,Dong Liangguo.Staggered-grid high-order finite-difference method in elastic wave simulation with variable grids and local time-steps.Chinese Journal of Geophysic,2009,52(11):2870-2878.
    [9]张慧,李振春.基于双变网格算法的地震波正演模拟.地球物理学报,2011,54(1):77-86.Zhang Hui,Li Zhenchun.Seismic wave simulation method based on dual-variable grid.Chinese Journal of Geophysic,2011,54(1):77-86.
    [10]Mi Tieliang,Ma Jianwei,HervéChauris et al.Multilevel adaptive mesh modeling for wave propagation in layer media.Journal of Seismic Exploration,2010,19(2):121-139.
    [11]Liu Y and Sen M K.A new time-space domain highorder finite-difference method for the acoustic wave equation.Journal of Computational Physics,2009,228(23):8779-8806.
    [12]Alkhalifah T,Fomel S and Biondi B.The space-time domain:Theory and modelling for anisotropic media.Geophysical Journal International,2001,144(1):105-113.
    [13]Ma X and Alkhalifah T.Wavefield extrapolation in pseudo-depth domain.73rd EAGE Conference&Exhibition Extended Abstracts,2011,A014.
    [14]Ma X and Alkhalifah T.Wavefield extrapolation in pseudo-depth domain.74rd EAGE Conference&Exhibition Extended Abstracts,2012,P212.
    [15]Liu Y,Sen M K.Finite-difference modeling with adaptive variable-length spatial operators.Geophysics,2011,76(4):T79-T89.
    [16]Thompson J F,Warsi Z U A,Mastin C W.Numerical Grid Generation:Foundations and Applications.Amsterdam:North-Holland,1985.

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