任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟

被引：69

作者：

裴正林 ^{[1
]}

机构：

[1] 石油大学(北京)CNPC物探重点实验室

来源：

石油地球物理勘探 | 2004年 / 06期

关键词：

任意起伏地表; 弹性波方程; 交错网格; 高阶有限差分;

D O I：

10.13810/j.cnki.issn.1000-7210.2004.06.002

中图分类号：

P631.4 [地震勘探];

学科分类号：

0818 ; 081801 ; 081802 ;

摘要：

任意起伏地表弹性波数值模拟主要涉及两个问题 ,其一是如何求解弹性波方程 ;其二是如何处理自由边界条件。本文首先从泰勒级数展开式出发 ,推导出交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型应力—速度弹性波方程交错网格任意偶数阶精度差分格式求解方程 ;然后采用将零速度法和广义虚像法相结合的方法来处理自由边界 ,并在自由边界上采用四阶精度差分格式 ;运用上述方法对光滑起伏模型和任意起伏模型进行了数值模拟试验。结果表明 ,本文所述方法稳定性好、模拟精度高 ,且适合任意起伏地表弹性波波场模拟。

引用

页码：629 / 634 +621-751

页数：8

共 14 条

[1] 任意偶数阶精度有限差分法数值模拟
刘洋
李承楚
牟永光
[J]. 石油地球物理勘探, 1998, (01) : 1 - 10
[2] High-order (space and time )finite-difference modeling of elastic wave equation. Crase E. Expanded Abstracts of SEG 60th Ann Internat Mtg . 1990
[3] 2D finite-difference elastic wave modeling including surface topography. Hestholm S O and Ruud B O. Geophysical Prospecting . 1994
[4] Spatial finite difference approximations for wave-type equation. Fornberg B and Ghrist M. SIAM Journal on Numerical Analysis . 1999
[5] Simulating seismic wave propagation in 3-D elastic media using staggered-grid finite difference. Graves R W. Bulletin of the Seismological Society of America . 1996
[6] A free-surface boundary conditions for including 3-D topography in the finite-difference method. Ohminato T and Chouet B A. Bulletin of the Seismological Society of America . 1996
[7] 3-D elastic modeling with surface topography by Chebychev spectral method. Tessmer E and Kosloff D. Geophysics . 1994
[8] Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Collino F and Tsogka C. Geophysics . 2001
[9] A numerical free-surface condition for elastic viscoelastic finite-difference modeling in the presence of topography. Robertsson J O A. Geophysics . 1996
[10] 3-D finite difference elastic wave modeling including surface topography. Hestholm S O and Ruud B O. Geophysics . 1998

← 1 2 →