摘要
<正>有限差分方法因其计算速度快、编程简便、并行性强等优势,在计算地球物理学中得到了极其广泛的应用。然而,长期以来,有限差分一直存在着两个亟待解决的问题,一是数值频散,二是人工边界条件。数值频散是一种由于连续的偏微分方程在差分算子的作用下,数值波速变得依赖于波数,因而与真实波速产生偏离的现象。数值频散扭曲了我们对真实波场和地震反演结果的认识,应当尽可能得去除。近似解析保辛分部Runge-Kutta(NSPRK)方法就是一种具有低数值频散的有限差分方法。它
引文
1.冯康.哈密尔顿系统的辛几何算法[M]//浙江科技出版社,2003.
2.马啸,杨顶辉,张锦华.求解声波方程的辛可分Runge-Kutta方法.地球物理学报,2010,53(008):1993-2003.
3.Komatitsch D,Martin R.An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation[J].Geophysics,2007,72(5):SM155.
4.Ma X,Yang D,Song G,et al.A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations:I.Scheme and Theoretical Analysis[J].Bulletin of the Seismological Society of America,2014,104.