可变网格多尺度有限差分模拟方法

《可变网格多尺度有限差分模拟方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、可变网格多尺度有限差分模拟方法　　摘要：针对典型各项同性地质模型，地震波波动方程和有限差分法对交错网格地震波场进行模拟仿真研究。为了得到一阶速度应力弹性波方程，首先对地震波波动方程原始形式进行推演，然后用时空二元高阶交错差分网格的方法对其进行离散化，最后对地质模型模拟仿真实验，并获取了对应的波场快照。最后由仿真结果表明，算法计算精度准确，计算效率高，而且可变网格的变换对地质模型的研究起到了很重要的作用。Abstract：Aimingatthetypicalgaygeologicalmodel，theseismicwaveequa

2、tionandfinitedifferencemethodareusedtodosimulationresearchtothestaggeredgridseismicwavefield.Inordertogetafirstorderrateofseismicwaveequation，theoriginalformoftheseismicwaveequationisdeducedfirstly.Thenspace-timedualhigherorderdifferencestaggeredgridmethodisusedtodod

3、iscretization.Finallythegeologicalmodelsimulationexperimentisdoneandthecorrespondingwavefieldsnapshotsareobtained.Thesimulationresultsshowthatthealgorithmisaccuratewith7highcomputationalefficiency，andthevarisionofvariablegridisofimportantsignificanceforthestudygeolog

4、icalmodel.关键词：有限差分；地震模型；波动方程；正演模型Keywords：finitedifference；seismicmodel；waveequation；forwardmodel中图分类号：O241.82文献标识码：A文章编号：1006-4311（2014）02-0315-020引言随着地震勘探面临的勘探任务越来越复杂，地震勘探成效在很大程度上，取决于对基于符合实际介质模型的方法理论研究。本文就详细介绍了有限差分法对交错网格地震波场进行的模拟仿真研究。一般的有限差分地震模拟方法是规则的交错网格。但是当地质体强纵横

5、向变速介质时计算效率却很低。为此，人们研究了基于可变网格和不规则网格的地震波数值模拟方法。它对地质模型的离散化更为合理，同时，对于保持模型计算的灵活性也非常重要。从空间采样的角度考虑，最有效的提高模拟精度同时又降低计算机内存需求的方法，就是在模型的不同区域采用不同的网格步长，即可变网格（也称为不规则网格）[1]。1正演模拟方法原理1.1均匀各向同性完全弹性介质的弹性波动方程7由本质上讲，地震波就是在地下岩层中传播的机械弹性波，但反映弹性波传播的基本规律是由弹性波的波动方程[2]。由研究岩石的不同物理参数得到均匀各项同性完全弹性介

6、质的波动方程，即■+■+■+ρg■=ρ■■+■+■+ρg■=ρ■■+■+■+ρg■=ρ■（1）式中：t为时间分量；u，v，ω为介质位移在x，y，z三个方向上的分量；σxx、σyy、σzz、τxy、τyz、τzx为应力分量；gx、gy、gz为体力密度分量；ρ为介质密度。1.2各向同性介质中的一阶弹性波方程弹性波理论是弹性体受力和应变的关系，根据各向同性介质表示的的本构方程和表示的柯西方程，可以推导出一阶速度—应力弹性波方程。二维二分量各向同性介质中的一阶应力—速度弹性波方程（假定体力为零）[4]，即ρ■=■+■ρ■=■+■■=（λ

7、+2μ）■+λ■■=λ■+（λ+2μ）■■=μ（■+■）（2）式中λ、μ表示拉梅弹性常数，ρ表示弹性体密度，σ■、σ■表示正应力，σ■表示切应力，vx、vz分别表示x方向、z方向速度分量。2网格有限差分法2.1有限差分法基本原理有限差分法（Finite7DifferentialMethod）是基于差分原理的一种数值计算方法，基本思想是：将波场域离散为许多小网格，应用差分原理，用差商来近似代替微商，将求解的连续函数问题转换为求解网格节点上的离散的差分方程组的问题[5]。2.2高阶方程化为低阶方程的一般方法解决高阶方程化低阶方程的问

8、题，其基本思路是将高阶差分分裂成低阶方程（二阶）。在每一延拓步长中交替使用低阶方程进行延拓，即可达到高阶差分延拓的目的。当色散方程Muir展开式的项数取得比较多时，一般式为■=-■■■p并可（在zj，zj+1间隔内）分裂成下列方程组■=-■■p■p■（z■）=p