王川-任意偶数阶精度交错网格粘声波方程数值模拟.docVIP

任意偶数阶精度的交错网格粘声波方程数值模拟

PAGE5

任意偶数阶精度交错网格粘声波方程数值模拟

摘要

如何提高正演模拟的精度，使之更符合实际地震波的传播规律，是地震波正演模拟的目标。本文基于粘声波方程，详细推导了其空间任意偶数阶精度交错网格差分格式。综合分析各参数的其模拟效果，应用空间八阶精度差分格式，结合完全匹配层（PML）吸收边界条件来模拟地震波在各种模型中的传播规律。常见模型及Marmousi模型实验结果显示，粘声波方程的高阶交错网格有限差分方法可精确描述多种地震波的传播规律，且能反映出地震波在传播过程中的频率耗散变化，其可应用于复杂构造模型的高精度地震波模拟。

关键字：交错网格、完全匹配层、粘声波方程

Baseonstaggered-gridanyeven-orderaccuracyseismicmodelingofviscoelasticmedia

ABSTRACT

Thegoalofseismicmodelingishowtoimprovetheaccuracyofseimicmodelingtomakeitmorerealisticofseismicwavesinpropagation.Basedontheviscoelasticwaveequationandstaggered-gridtechnical,thispaperderivedanyeven-orderaccuratedifferenceschemeinspace.Basetheanalysisofthemodelingresults,choosingeight-orderaccuracydifferenceschemeinspace,combinedwithperfectlymatchedlayer(PML)absorbingboundarytomodelingthepropagationofseismicwaveinvariousmodels.TheresultsofmodelinginMarmousimodelandsomeothercomplexmodelshowthatstaggered-stickhigh-orderaccuracyfinite-differenceequationwhichbaseonviscoelastictheorycanaccuratelydescribeavarietyofseismicwavepropagation,andcanreflectthefrequencydissipationofseismicwave.Aboveall,thetectoniccanbeappliedtomodelingofthehigh-precisionseismicwavepropagationincomplexmedia.

Keywords:Staggered-grid,viscoelastic,PML

目录

TOC\o1-3\u1前言 2

2一阶应力-速度方程组波动方程模拟 3

2.1一阶应力-速度方程组 3

2.2一阶应力-速度声波方程组交错网格任意偶数阶精度差分 4

2.2.1一阶应力-速度声波方程组交错网格计算实现 4

2.2.2任意偶数阶交错网格计算差分实现 6

2.2.3一阶应力-速度方程组的任意偶数阶精度交错网格下的离散形式实现 7

2.3地震子波选取 8

2.4完全匹配层（PML）吸收边界 9

2.5声波方程稳定性 14

2.6声波方程的正演数值模拟 14

2.6.1不同精度的声波方程模拟效果对比 14

2.6.2完全匹配层吸收边界效果对比 18

2.6.3完全弹性声波方程模拟效果 21

3粘弹性介质中波动方程 25

3.1粘弹性介质的基本情况 25

3.2标准线性模型的波动方程 25

3.3交错网格实现 26

3.4吸收边界与子波 27

3.5粘声波方程的正演数值模拟示例 28

4复杂介质中粘声波方程模拟 30

4.1断层模型 30

4.2倾斜界面地质模型 32

4.3盆地模型 34

4.4基于Marmousi复杂模型的模拟 36

5总结与建议 40

6致谢 41

参考文献 42

前言

随着地震勘探需求的增大，地震勘探技术的发展也十分迅速。作为地震勘探技术重要分支的地震勘探数值模拟技术，在地震波传播研究、偏移成像、反演等各方面有着关