一种基于全局优化的交错网格有限差分法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
在地震波场数值模拟中,交错网格有限差分技术得到了广泛的应用,但是在弹性模量变化较大时,通常会因插值而导致模拟误差增大.旋转交错网格可以很好地克服这个缺点,因而适合于各向异性介质正演模拟.但是对于同样大小的网格单元,旋转交错网格需要的步长比常规交错网格要大,这会使梯度和散度算子的误差增大因而更易产生空间数值频散.针对这些问题,本文提出了旋转交错网格与紧致有限差分相结合的方法,并基于模拟退火算法进行全局优化,压制数值频散,拓宽波数范围.数值模拟结果表明,此方法可以有效地压制数值频散,且具有较高的模拟精度.
Staggered grid finite difference techniques have been widely used in numerical simulation of seismic wave field,but the interpolation might enlarge the simulation errors when the elastic moduli change dramatically.Rotated staggered grid can overcome this drawback thus be more applicable for modeling the anisotropic media.However,for the grid cells with the same size,rotated staggered grid needs larger step than conventional staggered grid does,consequently,the gradient operator and divergence operator generate larger errors which are more likely to cause spatial numerical dispersion.Aiming at the above problems,a combination of rotated staggered grid and compact finite difference method is proposed,simultaneously,global optimization is preformed based on simulated annealing algorithm to suppress numerical dispersion and to broaden the range of wavenumber.Numerical simulation results show that this method can effectively suppress the numerical dispersion and exhibits high simulation accuracy at the same time.
引文
白晓寅.2008.基于地震波衰减理论的地层吸收参数提取方法研究[D].东营:中国石油大学(华东)地球科学与技术学院:24--54.Bai X Y.2008.The Study on Method of Layer Absorbing Parameters Extraction Based on the Theory of Seismic WaveAttenuation[D].Dongying:School of Geosciences,China University of Petroleum:24--54(in Chinese).
    陈浩,王秀明,赵海波.2006.旋转交错网格有限差分及其完全匹配层吸收边界条件[J].科学通报,51(17):1985--1994.Chen H,Wang X M,Zhao H B.2006.Rotated staggered grid finite-difference and the PML boundary[J].Chinese Science Bulletin,51(17):1985--1994(in Chinese).
    董良国.2003.地震波数值模拟与反演中几个关键问题研究[D].上海:同济大学海洋与地球科学学院:1--2.Dong L G.2003.Several Key Problems in Seismic Wave Numerical Simulation and Inversion[D].Shanghai:School ofOcean and Earth Science,Tongji University:1--2(in Chinese).
    何燕.2008.正交各向异性弹性波高阶有限差分正演模拟研究[D].东营:中国石油大学(华东)地球科学与技术学院:45--80.He Y.2008.High-Order Finite-Difference Forward Modeling of Elastic-Wave in Orthorhombic Anisotropic Media[D].Dongying:School of Geosciences,China University of Petroleum:45--80(in Chinese).
    李敏,刘洋.2012.高阶旋转交错网格有限差分方法模拟TTI介质中横波分裂[J].物探与化探,36(6):934--940.Li M,Liu Y.2012.Modeling of the S-wave splitting in TTI media using high-order rotated staggered grid scheme[J].Geophysical and Geochemical Exploration,36(6):934--940(in Chinese).
    裴正林,王尚旭.2005.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟[J].地震学报,27(4):441--451.Pei Z L,Wang S X.2005.A staggered-grid high-order finite-difference modeling for elastic wave field in arbitrary tilt anisotropic media[J].Acta Seismologica Sinica,27(4):441--451(in Chinese).
    孙成禹,肖云飞,印兴耀,彭洪超.2010.黏弹介质波动方程有限差分解的稳定性研究[J].地震学报,32(2):147--156.Sun C Y,Xiao Y F,Yin X Y,Peng H C.2010.Study on the stability of finite difference solution of visco-elastic waveequations[J].Acta Seismologica Sinica,32(2):147--156(in Chinese).
    孙林洁.2012.非均匀孔隙介质有限差分正演模拟方法研究[D].青岛:中国石油大学(华东)地球科学与技术学院:64--108.Sun L J.2012.Finite-Difference Modeling Research in Heterogeneous Porous Media[D].Qingdao:School of Geosciences,China University of Petroleum:64--108(in Chinese).
    王守东.2003.声波方程完全匹配层吸收边界[J].石油地球物理勘探,38(1):31--34.Wang S D.2003.Absorbing boundary condition for acoustic wave equation by perfectly matched layer[J].Oil Geophysi-cal Prospecting,38(1):31--34(in Chinese).
    王书强,杨顶辉,杨宽德.2002.弹性波方程的紧致差分方法[J].清华大学学报:自然科学版,42(8):1128--1131.Wang S Q,Yang D H,Yang K D.2002.Compact finite difference scheme for elastic equations[J].Journal of TsinghuaUniversity:Science and Technology,42(8):1128--1131(in Chinese).
    严红勇,刘洋.2012.黏弹TTI介质中旋转交错网格高阶有限差分数值模拟[J].地球物理学报,55(4):1354--1365.Yan H Y,Liu Y.2012.Rotated staggered grid high-order finite-difference numerical modeling for wave propagation inviscoelastic TTI media[J].Chinese Journal of Geophysics,55(4):1354--1365(in Chinese).
    殷文,印兴耀,吴国忱,梁锴.2006.高精度频率域弹性波方程有限差分方法及波场模拟[J].地球物理学报,49(2):561--568.Yin W,Yin X Y,Wu G C,Liang K.2006.The method of finite difference of high precision elastic wave equations in thefrequency domain and wave field simulation[J].Chinese Journal of Geophysics,49(2):561--568(in Chinese).
    Boersma B J.2011.A 6th order staggered compact finite difference method for the incompressible Navier-Stokes and scalar transport equations[J].J Comput Phys,230(12):4940--4954.
    Collino F,Tsogka C.2001.Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,66(1):294--307.
    Dai N,Vafidis A,Kanasewich E R.1995.Wave propagation in heterogeneous,porous media:A velocity-stress,finitedifference method[J].Geophysics,60(2):327--340.
    Du Q Z,Li B,Hou B.2009.Numerical modeling of seismic wavefields in transversely isotropic media with a compactstaggered-grid finite difference scheme[J].Appl Geophys,6(1):42--49.
    Gold N,Shapiro S A,Burr E.1997.Modelling of high contrasts in elastic media using a modified finite difference scheme[C]∥SEG Technical Program Expanded Abstracts.Dallas,Texas:SEG:1850--1853.
    Holberg O.1987.Computational aspects of the choice of operator and sampling interval for numerical differentiation inlarge-scale simulation of wave phenomena[J].Geophys Prospect,35(6):629--655.
    Igel H,Mora P,Riollet B.1995.Anisotropic wave propagation through finite-difference grids[J].Geophysics,60(4):1203--1216.
    Kim J W,Lee D J.1996.Optimized compact finite difference schemes with maximum resolution[J].AIAA J,34(5):887--893.
    Kirkpatrick S,Gelatt C D Jr,Vecchi M P.1983.Optimization by simulated annealing[J].Science,220(4598):671--680.
    Kosloff D,Pestana R C,Tal-Ezer H.2010.Acoustic and elastic numerical wave simulations by recursive spatial derivative operators[J].Geophysics,75(6):T167--T174.
    Lele S K.1992.Compact finite difference schemes with spectral-like resolution[J].J Comput Phys,103(1):16--42.
    Liu Y,Sen M K.2009.An implicit staggered-grid finite-difference method for seismic modelling[J].Geophys J Int,179(1):459--474.
    Liu Y.2013.Globally optimal finite-difference schemes based on least squares[J].Geophysics,78(4):T113--T132.
    Saenger E H,Bohlen T.2004.Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid[J].Geophysics,69(2):583--591.
    Saenger E H,Gold N,Shapiro S A.2000.Modeling the propagation of elastic waves using a modified finite-differencegrid[J].Wave Motion,31(1):77--92.
    Shan G.2009.Optimized implicit finite-difference and Fourier finite-difference migration for VTI media[J].Geophysics,74(6):WCA189--WCA197.
    Virieux J.1984.SH-wave propagation in heterogeneous media:Velocity-stress finite-difference method[J].Geophysics,49(11):1933--1957.
    Zeng Y Q,He J Q,Liu Q H.2001.The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media[J].Geophysics,66(4):1258--1266.
    Zhang J H,Yao Z X.2012.Optimized finite-difference operator for broadband seismic wave modeling[J].Geophysics,78(1):A13--A18.
    Zhou H B,Zhang G Q.2011.Prefactored optimized compact finite-difference schemes for second spatial derivatives[J].Geophysics,76(5):WB87--WB95.

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