Finite-difference(FD)method is the most extensively employed numerical modeling technique.Nevertheless,when using the FD method to simulate the seismic wave propagation,the large spatial or temporal sampling interval ...Finite-difference(FD)method is the most extensively employed numerical modeling technique.Nevertheless,when using the FD method to simulate the seismic wave propagation,the large spatial or temporal sampling interval can lead to dispersion errors and numerical instability.In the FD scheme,the key factor in determining both dispersion errors and stability is the selection of the FD weights.Thus,How to obtain appropriate FD weights to guarantee a stable numerical modeling process with minimum dispersion error is critical.The FD weights computation strategies can be classified into three types based on different computational ideologies,window function strategy,optimization strategy,and Taylor expansion strategy.In this paper,we provide a comprehensive overview of these three strategies by presenting their fundamental theories.We conduct a set of comparative analyses of their strengths and weaknesses through various analysis tests and numerical modelings.According to these comparisons,we provide two potential research directions of this field:Firstly,the development of a computational strategy for FD weights that enhances stability;Secondly,obtaining FD weights that exhibit a wide bandwidth while minimizing dispersion errors.展开更多
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology (Qingdao) (No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group (No.SNKJ2022A06-R23)+2 种基金the Funds of Creative Research Groups of China (No.41821002)National Natural Science Foundation of China Outstanding Youth Science Fund Project (Overseas) (No.ZX20230152)the Major Scientific and Technological Projects of CNPC (No.ZD2019-183-003)。
文摘Finite-difference(FD)method is the most extensively employed numerical modeling technique.Nevertheless,when using the FD method to simulate the seismic wave propagation,the large spatial or temporal sampling interval can lead to dispersion errors and numerical instability.In the FD scheme,the key factor in determining both dispersion errors and stability is the selection of the FD weights.Thus,How to obtain appropriate FD weights to guarantee a stable numerical modeling process with minimum dispersion error is critical.The FD weights computation strategies can be classified into three types based on different computational ideologies,window function strategy,optimization strategy,and Taylor expansion strategy.In this paper,we provide a comprehensive overview of these three strategies by presenting their fundamental theories.We conduct a set of comparative analyses of their strengths and weaknesses through various analysis tests and numerical modelings.According to these comparisons,we provide two potential research directions of this field:Firstly,the development of a computational strategy for FD weights that enhances stability;Secondly,obtaining FD weights that exhibit a wide bandwidth while minimizing dispersion errors.
