Accuracy of the staggered-grid finite-difference method of the acoustic wave equation for marine seismic reflection modeling  

作　　者：钱进 吴时国 崔若飞 

机构地区：[1]Key Laboratory of Marine Geology and Environment,Institute of Oceanology,Chinese Academy of Sciences [2]School of Resource and Geoscience,China University of Mining and Technology

出　　处：《Chinese Journal of Oceanology and Limnology》2013年第1期169-177,共9页中国海洋湖沼学报（英文版）

基　　金：Supported by the National Natural Science Foundation of China(Nos. 41206043, 40930845);the Open Foundation of Key Laboratory of Marine Geology and Environment of Chinese Academy of Sciences(No. MGE2011KG07);the Knowledge Innovation Program of Chinese Academy of Sciences (No. KZCX2-YW-229);the National Basic Research Program of China (973 Program) (No. 2009CB219505)

摘　　要：Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient (SC) and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient （SC） and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.

关 键 词：marine seismic reflection modeling stability condition dispersion relation staggered grid finite-difference 

分 类 号：P714[天文地球—海洋科学]