基于辛算法及FCT的地震波场模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
为了解决在长时间和复杂结构的地震波场数值模拟情况时的数值频散问题,现在辛算法的基础上,主要结合通量校正传输(FCT)、褶积微分算子、完全匹配层(PML)等数值模拟技术,寻找一种更为优秀的地震波场数值模拟方法.地震波场的数值模拟结果表明,辛算法不仅具有保持体系结构的特性,并且具有长时间跟踪能力,具有很强的数值模拟稳定性;FCT方法基于通量守恒原理,压制网格频散效果明显;褶积微分算子突出了空间微分的局部属性.因此,通过合理应用各种技术,可以对地震波场特征进行更精确的数值模拟.
To find out a solution of numerical dispersion in long-time tracing and complicated structure modeling,based upon the symplectical algorithm,this paper introduced a 2D seismic wavefield modeling technique which combins with the Flux-corrected Transport method(FCT) and convolutional differentiator,and uses the Perfectly Matched Layer(PML) as a absorbing boundaries.The numerical modeling shows that the Symplectical algorithm not only has a structure-preserving feature but also has a long-time tracing ability,it has a good stability during numerical modeling.Due to adopting convolutional differentiator and FCT techniques,the spatial feature is enhanced,and the numerical dispersion caused by the gridding is suppressed to produce a more accurate seismic wavefield.So we can get a better result with combins these technique reasonable.
引文
Aki K,Larner K L,1970.Surface motion of a layered medium having an irregular interface due to incident plane SH waves[J].J.Geophys.Res,75(5):933-954.
    Alterman Z,Karal F C Jr.1968.Propagation of ela.stic waves in layered media by finite difference methods[J].Bull.Seism.Soc.Am.,58(1):367-398.
    Boris J P,Book D L.1973.Flux-corrected transport.I.SHASTA,A fluid transport algorithm that works[J].J.Comp.Phys.,11(1):38-69.
    Feng K.Qin M Z.2003.Symplectlc Geometric Algorithm for Hamiltonian Systems(in Chinese)[M].Hangzhou;Zhejiang Science and Technology Press.
    Huang C,Dong L G.2009.Staggered-grid high-order finitedifference method in elastic wave simulation with variable grids and local time-steps[J].Chinese Journal of Geophysics(in Chinese),52(11):2870-2878,doi:10.3969/j.issa 0001-5733.2009.11.022.
    Iwatsu R 2009.Two new solutions to the third-order symplectic integration method[J].Phys.Lett.A,373(34):3056-3060.
    Li S X.2001.The High-Frequency Approximation of Wave Equations and Symplectic Algorithm(in Chinese)[M].Beijing:Science Press.
    Li Y Q,Li X F,Zhu T.2011.The seismic scalar wave field modeling by symplectic scheme and singular kernel convolutional differentiator[J].Chinese J.Geophys.(in Chinese),54(7):1827-1834,doi:10.3969/j.issa 0001-5733.2011.07.016.
    Li Y Q,Li X F,Zhang M G.2012.The elastic wave fields modeling by symplectic discrete singular convolution differentiator method[J].Chinese J.Geophys.(in Chinese),55(5):1725-1731.doi:10.6038/j.issn.0001-5733.2012.05.029.
    Long G H,Li X F,Zhang M G.2009.The application of convolutional differentiator in seismic modeling based on Shannon singular kernel theory[J].Chinese J.Geophys.(in Chinese),52(4):1014-1024.
    Luo M Q,Liu H,Li Y M 2001.Hamiltonian description and symplectic method of seismic wave propagation[J].Chinese Journal of Geophysics(in Chinese),44(1):120-128.
    Qin M Z,Chen J B.2000.Maslov asymptotic theory and symplectic algorithm[J].Chinese Journal of Geophysics(in Chinese),43(4):522-533.
    Smith W D.1975.The Application of finiteelement analysis to body wave propagation problems[J].Geophysics J.Roy Astr.Soc.,42(2):747-768.
    Sun L J,Yin X Y.2011.A finite-difference scheme based on PML boundary condition with high power grid step variation[J].Chinese Journal of Geophysics(in Chinese),54(6):1614-1623,doi:10.3969/j.issn.0001-5733.2011.06.021.
    Suzuki M.1992.General theory of higher-order decomposition of exponential operators and symplectic integrators[J].Phys.Lett.A,165(5-6):387-395.
    Sun Y H,Zhou Y C,Li S G,et al.2006.A windowed Fourier pseudospectral method for hyperbolic conservation laws[J].J.Comp.Appl.Math.,214(2):466-490.
    Wang W S,Li X F,Lu M W,et al.2012.Structure-preserving modeling for seismic wavefields based upon a multisymplectic spectral element method[J].Chinese Journal of Geophysics(in Chinese),55(10):3427-3439,doi:10.6038/j.issn.0001-5733.2012.10.026.
    Wei G W.1999.Discrete singular convolution for the solution of the Fokker-Planck equation[J].J.Chem.Phys.,110(18):8930-8942.
    Yang K D,Song G J,Li J S.2011.FCT compact difference simulation of wave propagation based on the Biot and the squirtflow coupling interaction[J].Chinese J.Geophys.(in Chinese),54(5):1348-1357,doi:10.3969/j.issn.0001-5733.2011.05.024.
    Yoshida H.1990.Construction of higher order symplectic integrators[J].Phys.Lett.A,150(5-7):262-269.
    Zhou B,Greenhalgh S A.1992.Seismic scalar wave equation modeling by a convolutional differentiator[J].Bull.Seism.Soc.Am,82(1):289-303.
    冯康,秦孟兆.2003.哈密尔顿系统的辛几何算法[M].杭州:浙江科学技术出版社.
    黄超,董良国.2009.可变网格与局部时间步长的交错网格高阶差分弹性波模拟[J].地球物理学报,52(11):2870-2878,doi:10.3969/j.issn.0001-5733.2009.11.022.
    李世雄.2001.波动方程的高频近似与辛几何[M].北京:科学出版社.
    李一琼,李小凡,朱童.2011.基于辛格式奇异核褶积微分算子的地震标量波场模拟[J].地球物理学报,54(7):1827-1834,doi:10.3969/j.issn.0001-5733.2011.07.016.
    李一琼,李小凡,张美根.2012.基于辛格式离散奇异褶积微分算子的弹性波场模拟EJ3.地球物理学报,55(5):1725-1731,doi:10.6038/j.issn.0001-5733.2012.05.029.
    龙桂华,李小凡,张美根.2009.基于Shannon奇异核理论的褶积微分算子在地震波场模拟中的应用[J].地球物理学报,52(4):1014-1024,doi:10.3969/j.issn.0001-5733.2009.04.018.
    罗明秋,刘洪,李幼铭.2001.地震波传播的哈密顿表述及辛几何算法[J].地球物理学报,44(1):120-128.
    秦孟兆,陈景波.2000.Maslov渐近理论与辛几何算法[J].地球物理学报,43(4):522-533.
    孙林洁,印兴耀.2011.基于PML边界条件的高倍可变网格有限差分数值模拟方法[J].地球物理学报,54(6):1614-1623,doi:10.3969/j.issn.0001-5733.2011.06.021.
    汪文帅,李小凡,鲁明文,等.2012.基于多辛结构谱元法的保结构地震波场模拟[J].地球物理学报,55(10):3427-3439,doi:10.6038/j.issn.0001-5733.2012.10.026.
    杨宽德,宋国杰,李静爽.2011.Biot流动和喷射流动耦合作用下波传播的FCT紧致差分模拟[J].地球物理学报,54(5):1348-1357,doi:10.3969/j.issn.0001-5733.2011.05.024.

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