本文在对地震波场进行模拟时,采用辛差分格式对波动方程进行时间离散,采用奇异核褶积微分算子对波动方程进行空间离散.该方法尽管增加了一些计算量,但提高了计算精度和稳定性;相对于其他非辛算法,它是全局保结构的,并且具有较强的长时间跟踪能力.该方法为解决大尺度、长时程地震波场的高精度模拟问题提供了一种新的、有效的选择.
The paper adopts both symplectic difference scheme and singular kernel convolutional differentiator for modeling of seismic scalar wave field. For this method, the computational accuracy and stability have been greatly improved with a slight increase of calculation amount. Compared with other non-symplectic numerical methods, major advantages of the method presented in this paper are the structure-preserving property and long-term simulations. These appealing characters of the algorithm would make it effective to model the large-scale and long- term seismic wave propagation.