高阶交错网格有限差分弹性波场模拟的精度分析  

Accuracy analysis of elastic wave field simulation based on high-order staggered grid finite difference scheme

作　　者：岳晓鹏[1,2] 白超英[1] 岳崇旺[1] YUE Xiaopeng;BAI Chaoying;YUE Chongwang(Department of Geophysics, College of Geology Engineering and Geomatics, Chang’an University, Xi’an 710064,China;School of Mathematics and Statistics, Xuchang University, Xuchang 461000, China)

机构地区：[1]长安大学地质工程与测绘学院地球物理系,陕西西安710064 [2]许昌学院数学与统计学院,河南许昌461000

出　　处：《煤田地质与勘探》2017年第1期125-130,共6页Coal Geology & Exploration

基　　金：国家自然科学基金青年科学基金项目(41504090);中央高校基金基础研究项目(310826161021)~~

摘　　要：交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明:时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。Seismic wavefield numerical simulation based on the staggered grids is widely used in seismic forward atpresent,there will be some differences in efficiency and accuracy by using different order difference scheme to simulate.In order to compare these differences,four kinds of difference formulas and coefficients of4order temporal and2N order spatial of elastic wave equation are deduced,and the stability conditions of high order staggered grid underdifferent difference schemes are calculated.The elastic wave field is simulated by using the four difference formulasand coefficients,and the results of wave field snapshot,synthetic seismogram and CPU time are compared and analyzed.The results shows that the staggered grid with accuracy of4order temporal and6+6orders spatial has highcalculation precision and computational efficiency in numerical simulation of seismic wave field.

关 键 词：精度分析 弹性波场 交错网格 高阶有限差分 

分 类 号：P315.69[天文地球—地震学]