Time-space domain scalar wave modeling by a novel hybrid staggered-grid finite-difference method with high temporal and spatial accuracies

被引:6
作者
Zhou, Hongyu [1 ,3 ]
Liu, Yang [1 ,2 ,3 ]
Wang, Jing [1 ,3 ]
机构
[1] China Univ Petr, State Key Lab Petr Resources Prospecting, Beijing 102249, Peoples R China
[2] China Univ Petr, Karamay Campus, Karamay 834000, Peoples R China
[3] China Univ Petr, CNPC Key Lab Geophys Prospecting, Beijing 102249, Peoples R China
基金
中国国家自然科学基金;
关键词
Finite-difference; Wave equation; Modeling; Staggered-grid; HETEROGENEOUS MEDIA; ELEMENT SCHEMES; ORDER ACCURACY; LEAST-SQUARES; P-WAVE; PROPAGATION; EQUATION; GALERKIN; RESPECT;
D O I
10.1016/j.jcp.2022.111004
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Staggered-grid finite-difference (SFD) scheme is favored in wave equation simulation due to its superior accuracy and stability to center-grid finite-difference (CFD) scheme. However, for the scalar wave equation (SWE) modeling, the conventional SFD (CSFD) scheme only reaches second-order accuracy in space and time for the discrete SWE; the recently developed modified SFD (MSFD) scheme improves the accuracy to (2N)th-order (N < 4) but is costly because the MSFD scheme is coined by adding many extra grid points to the CSFD scheme. To tackle these issues, we develop a cost-effective hybrid SFD (HSFD) scheme, which combines the features of the CSFD and MSFD schemes; we prove that the new HSFD scheme can simultaneously reach (2N)th-order accuracy in space and time for the discrete wave equation. In addition, to deal with the optimization difficulties due to the nonlinear dispersion relation of the SFD schemes, we propose a two-step linear optimization method to improve the accuracy of the new HSFD scheme. The analyses on dispersion, stability properties and numerical simulation examples demonstrate that the new HSFD scheme owns better accuracy and stability than the CSFD scheme. Computational cost analysis shows that the HSFD scheme can be more efficient than the MSFD scheme because it only requires half the additional grid points. (C) 2022 Elsevier Inc. All rights reserved.
引用
收藏
页数:29
相关论文
共 62 条
[1]   ACCURACY OF FINITE-DIFFERENCE MODELING OF ACOUSTIC-WAVE EQUATION [J].
ALFORD, RM ;
KELLY, KR ;
BOORE, DM .
GEOPHYSICS, 1974, 39 (06) :834-842
[2]   SEISMIC WAVES IN A QUARTER AND 3-QUARTER PLANE [J].
ALTERMAN, ZS ;
LOEWENTH.D .
GEOPHYSICAL JOURNAL OF THE ROYAL ASTRONOMICAL SOCIETY, 1970, 20 (02) :101-&
[3]  
[Anonymous], 2005, Computational electrodynamics: The finite-difference time-domain method
[4]  
Berkhout A.J., 1982, SEISMIC MIGRATION IM
[5]  
Boore D. M., 1972, Methods Comput Phys Adv Res Appl, V11, P1, DOI [10.1016/B978-0-12-460811-5.50006-4, DOI 10.1016/B978-0-12-460811-5.50006-4]
[6]  
Castagna J. P., 1993, ROCK PHYS LINK ROCK
[7]  
Chen H., 1995, SEG TECHNICAL PROGRA, P1289
[8]  
Chen H., 1996, SEG TECHNICAL PROGRA, P797
[9]   A k-space operator-based least-squares staggered-grid finite-difference method for modeling scalar wave propagation [J].
Chen, Hanming ;
Zhou, Hui ;
Zhang, Qingchen ;
Xia, Muming ;
Li, Qingqing .
GEOPHYSICS, 2016, 81 (02) :T45-T61
[10]  
Claerbout J, 1985, Imaging the Earths Interior