A general optimal method for a 2D frequency-domain finite-difference solution of scalar wave equation
A general optimal method for a 2D frequency-domain finite-difference solution of scalar wave equation
复制标题
DOI:
10.1190/geo2016-0457.1
复制
发表时间:
2017-04
期刊:
影响因子:
3.3
通讯作者:
N. Fan;Lian‐Feng Zhao;Xiao-Bi Xie;Xin-Gong Tang;Z. Yao
中科院分区:
文献类型:
--
作者:
N. Fan;Lian‐Feng Zhao;Xiao-Bi Xie;Xin-Gong Tang;Z. Yao
ABSTRACTWe have developed a general optimal method for 2D frequency-domain finite-difference simulation of the scalar wave equation. For a given finite-difference stencil, this method can generate the dispersion equation and optimize the expansion coefficients. Many commonly used frequency-domain finite-difference schemes (e.g., grids with different numbers of points, rotated grids, and grid spaces with different aspect ratios) can be derived as special cases under this framework. The possibility of expanding this method to 3D does exist. Based on the 2D scalar wave equation, the optimized coefficients of 25-point, 9-point, 17-point, and 15-point schemes have been worked out. The dispersion analysis indicates that our 25-point scheme has much higher accuracy than the average-derivative method 25-point scheme. The number of grid points per the smallest wavelength is reduced from 2.78 to 2.13 for a maximum phase velocity errors of 1%. The synthetic seismograms and the wavefield snapshots calculated using ou...