基于辛格式奇异核褶积微分算子的地震标量波场模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
本文在对地震波场进行模拟时,采用辛差分格式对波动方程进行时间离散,采用奇异核褶积微分算子对波动方程进行空间离散.该方法尽管增加了一些计算量,但提高了计算精度和稳定性;相对于其他非辛算法,它是全局保结构的,并且具有较强的长时间跟踪能力.该方法为解决大尺度、长时程地震波场的高精度模拟问题提供了一种新的、有效的选择.
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.
引文
[1] 李世雄.波动方程的高频近似与辛几何.北京:科学出版社,2001Li S X. The High-Frequency Approximation of Wave Equations and Symplectic Algorithm (in Chinese). Beijing:Science Press, 2001
    [2] 陈景波,秦孟兆.辛几何算法在射线追踪中的应用.数值计算与计算机应用,2000,21(4) :255~265Chen J B, Qin M Z. Ray tracing by symplectic algorithm. Journal of Numerical Methods and Computer Applications (in Chinese), 2000, 21(4) : 255-265
    [3] 高亮,李幼铭,陈旭荣等.地震射线辛几何算法初探.地球物理学报,2000,43(3) :402~410Gao L, Li Y M, Chen X R, et al. An attempt to seismic ray tracing with symplectic algorithm. Chinese J . Geophys. (inChinese), 2000, 43(3) : 402-410
    [4] Chen J B. Lax-Wendroff and Nystrom methods for seismic modeling. Geophysical Prospecting, 2009, 57(6) : 931-941
    [5] Holberg O. Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena. Geophysical Prospecting,2006, 35(6) : 629-655
    [6] Zhou B, Greenhalgh S. Seismic scalar wave equation modeling by a convolutional differentiator. Bull. Seism.Am. , 1992, 82(1) : 289-303
    [7] 张中杰,滕吉文,杨顶辉.声波与弹性波场数值模拟中的褶积微分算子法.地震学报,1996,18(1) :63~69Zhang Z J, Teng J W, Yang D H. The convolutional differentiator method for numerical modeling of acoustic and elastic wave-field. Acta Seismologica Sinica (in Chinese),1996, 18(1) : 63-69
    [8] 戴志阳,孙建国,查显杰.地震波场模拟中的褶积微分算子法.吉林大学学报(地球科学版),2005,35(4) :520~524Dai Z Y, Sun J G, Cha X J. Seismic wave field modeling with convolutional differentiator algorithm. Journal of Jilin University (Earth Science Edition) (in Chinese), 2005, 35(4) : 520-524
    [9] 龙桂华,李小凡,张美根.基于Shannon奇异核理论的褶积微分算子在地震波场模拟中的应用.地球物理学报,2009,52(4) :1014~1024Long G H, Li X F, Zhang M G. The application of convolutional differentiator in seismic modeling based on Shannon singular kernel theory. Chinese J. Geophys. (in Chinese), 2009, 52(4) : 1014-1024
    [10] Suzuki M. General theory of higher-order decomposition of exponential operators and symplectic integrators. Phys.Lett. A, 1992, 165(5-6) : 387-395
    [11] Yoshida H. Construction of higher order symplectic integrators. Phys. Lett. A, 1990, 150(5-7) : 262-269
    [12] Ruth R D. A canonical integration technique. IEEE Transactions on Nuclear Science , 1983, 30(4) : 2669-2671
    [13] Iwatsu R. Two new solutions to the third-order symplectic integration method. Phys. Lett. A, 2009, 373(34) : 3056-3060
    [14] Wei G W. Discrete singular convolution for the solution of the Fokker-Planck equation. J. Chem. Phys. , 1999, 110 (18) ; 8930-8942
    [15] 黄志祥,沙威,吴先良等.高阶辛算法的稳定性与数值色散性分析.计算物理,2010,27(1) :82~88Huang Z X, Sha W, Wu X L, et al. Stability and numerical dispersion of high order symplectic schemes. Chinese Journal of Computational Physics (in Chinese), 2010, 27(1) : 82-88
    [16] Hirono T, Lui W, Seki S, et al. A three-dimensional fourth-order finite-difference time-domain scheme using a symplectic integrator propagator. IEEE Transactions on Microwave Theory and Techniques, 2001, 49(9) . 1640-1648
    [17] Sha W, Huang Z X, Wu X L, et al. Total field and scattered field technique for fourth-order symplectic finite-difference time-domain method. Chinese Physics Letters , 2006, 23(1) : 103-105
    [18] Van de Vyver H. A symplectic Runge-Kutta-Nystrom method with minimal phase-lag. Physics Letters A , 2007,367(1-2) : 16-24
    [19] Dopico F M, Johnson C R. Complementary bases in symplectic matrices and a proof that their determinant is one.Linear Algebra and Its Applications , 2006, 419(2-3) : 772-778
    [20] 董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究.地球物理学报,2000,43(6) :856~864Dong L G, Ma Z T, Cao J Z. A staggeredgrid high-order difference method of one-order elastic wave equation. Chinese J. Geophys. (in Chinese), 2000, 43(6) : 856-864
    [21] Kong L H, Liu R, Zheng X H. A survey on symplectic and multi-symplectic algorithms. Applied Mathematics and Computation, 2007, 186(1) : 670-684

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