An Efficient Second-Order Finite-Difference Scheme for 3D Time-Space Domain Acoustic Wave Equation

Author(s):  
N. Amini ◽  
C. Shin
Geophysics ◽  
2021 ◽  
pp. 1-83
Author(s):  
Navid Amini ◽  
Changsoo Shin ◽  
Jaejoon Lee

We propose compact implicit finite-difference (FD) schemes in time-space domain based on second-order FD approximation for accurate solution of the acoustic wave equation in 1D, 2D, and 3D. Our method is based on weighted linear combination of the second-order FD operators with different spatial orientations to mitigate numerical error anisotropy and weighted averaging of the mass acceleration term over the grid-points of the second-order FD stencil to reduce the overall numerical dispersion error. We present derivation of the schemes for 1D, 2D, and 3D cases and obtain their corresponding dispersion equations, then we find optimum weights by optimization of the time-space domain dispersion function and finally tabulate the optimized weights for each case. We analyze the numerical dispersion, stability and convergence rates of the proposed schemes and compare their numerical dispersion characteristics with the standard high-order ones. We also discuss efficient solution of the system of equations associated with the proposed implicit schemes using conjugate gradient method. The comparison of dispersion curves and the numerical solutions with the analytical and the pseudo-spectral solutions reveals that the proposed schemes have better performance than the standard spatial high-order schemes and remain stable for relatively large time-steps.


2018 ◽  
Vol 49 (6) ◽  
pp. 898-905 ◽  
Author(s):  
Wenquan Liang ◽  
Xiu Wu ◽  
Yanfei Wang ◽  
Jingjie Cao ◽  
Chaofan Wu ◽  
...  

2013 ◽  
Vol 56 (6) ◽  
pp. 840-850 ◽  
Author(s):  
LIANG Wen-Quan ◽  
YANG Chang-Chun ◽  
WANG Yan-Fei ◽  
LIU Hong-Wei

2016 ◽  
Vol 34 (4) ◽  
Author(s):  
Laura Lara Ortiz ◽  
Reynam C. Pestana

ABSTRACT. In this work we show that the solution of the first order differential wave equation for an analytical wavefield, using a finite-difference scheme in time, follows exactly the same recursion of modified Chebyshev polynomials. Based on this, we proposed a numerical...Keywords: seismic modeling, acoustic wave equation, analytical wavefield, Chebyshev polinomials. RESUMO. Neste trabalho, mostra-se que a solução da equação de onda de primeira ordem com um campo de onda analítico usando um esquema de diferenças finitas no tempo segue exatamente a relação de recorrência dos polinômios modificados de Chebyshev. O algoritmo...Palavras-chave: modelagem sísmica, equação da onda acústica, campo analítico, polinômios de Chebyshev.


Sign in / Sign up

Export Citation Format

Share Document