A Compact High-Order Finite-Difference Method with Optimized Coefficients for 2D Acoustic Wave Equation

被引:4
作者
Chen, Liang [1 ]
Huang, Jianping [1 ]
Fu, Li-Yun [1 ]
Peng, Weiting [1 ]
Song, Cheng [1 ]
Han, Jiale [1 ]
机构
[1] China Univ Petr East China, Sch Earth Sci, Qingdao 266000, Peoples R China
基金
国家重点研发计划;
关键词
wavefield simulation; compact finite difference; acoustic wave equation; numerical dispersion; optimization schemes; SPECTRAL-ELEMENT SIMULATIONS; PSEUDOSPECTRAL METHOD; NUMERICAL-SIMULATION; HETEROGENEOUS MEDIA; CONVOLUTIONAL PML; P-SV; PROPAGATION; SCHEMES; ACCURACY; DISPERSION;
D O I
10.3390/rs15030604
中图分类号
X [环境科学、安全科学];
学科分类号
08 ; 0830 ;
摘要
High-precision finite difference (FD) wavefield simulation is one of the key steps for the successful implementation of full-waveform inversion and reverse time migration. Most explicit FD schemes for solving seismic wave equations are not compact, which leads to difficulty and low efficiency in boundary condition treatment. Firstly, we review a family of tridiagonal compact FD (CFD) schemes of various orders and derive the corresponding optimization schemes by minimizing the error between the true and numerical wavenumber. Then, the optimized CFD (OCFD) schemes and a second-order central FD scheme are used to approximate the spatial and temporal derivatives of the 2D acoustic wave equation, respectively. The accuracy curves display that the CFD schemes are superior to the central FD schemes of the same order, and the OCFD schemes outperform the CFD schemes in certain wavenumber ranges. The dispersion analysis and a homogeneous model test indicate that increasing the upper limit of the integral function helps to reduce the spatial error but is not conducive to ensuring temporal accuracy. Furthermore, we examine the accuracy of the OCFD schemes in the wavefield modeling of complex structures using a Marmousi model. The results demonstrate that the OCFD4 schemes are capable of providing a more accurate wavefield than the CFD4 scheme when the upper limit of the integral function is 0.5 pi and 0.75 pi.
引用
收藏
页数:19
相关论文
共 59 条
[1]   Three-dimensional seismic full-waveform inversion using the finite-difference contrast source inversion method [J].
Abubakar, A. ;
Pan, G. ;
Li, M. ;
Zhang, L. ;
Habashy, T. M. ;
van den Berg, P. M. .
GEOPHYSICAL PROSPECTING, 2011, 59 (05) :874-888
[2]  
Adam Y., 1975, Computers & Mathematics with Applications, V1, P393, DOI 10.1016/0898-1221(75)90041-3
[3]   ACCURACY OF FINITE-DIFFERENCE MODELING OF ACOUSTIC-WAVE EQUATION [J].
ALFORD, RM ;
KELLY, KR ;
BOORE, DM .
GEOPHYSICS, 1974, 39 (06) :834-842
[4]   REVERSE TIME MIGRATION [J].
BAYSAL, E ;
KOSLOFF, DD ;
SHERWOOD, JWC .
GEOPHYSICS, 1983, 48 (11) :1514-1524
[5]  
Bouchon M, 1995, B SEISMOL SOC AM, V85, P1679
[6]   A three-point combined compact difference scheme [J].
Chu, PC ;
Fan, CW .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 140 (02) :370-399
[7]   An efficient fourth-order low dispersive finite difference scheme for a 2-D acoustic wave equation [J].
Das, Sambit ;
Liao, Wenyuan ;
Gupta, Anirudh .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2014, 258 :151-167
[8]   NUMERICAL-SIMULATION OF ELASTIC-WAVE PROPAGATION USING A FINITE-VOLUME METHOD [J].
DORMY, E ;
TARANTOLA, A .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 1995, 100 (B2) :2123-2133
[9]   Numerical modeling of seismic wavefields in transversely isotropic media with a compact staggered-grid finite difference scheme [J].
Du Qizhen ;
Bin, Li ;
Bo, Hou .
APPLIED GEOPHYSICS, 2009, 6 (01) :42-49
[10]   2D and 3D elastic wave propagation by a pseudo-spectral domain decomposition method [J].
Faccioli, E ;
Maggio, F ;
Paolucci, R ;
Quarteroni, A .
JOURNAL OF SEISMOLOGY, 1997, 1 (03) :237-251