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

可变网格多尺度有限差分模拟方法 　　摘要： 针对典型各项同性地质模型，地震波波动方程和有限差分法对交错网格地震波场进行模拟仿真研究。为了得到一阶速度应力弹性波方程，首先对地震波波动方程原始形式进行推演，然后用时空二元高阶交错差分网格的方法对其进行离散化，最后对地质模型模拟仿真实验，并获取了对应的波场快照。最后由仿真结果表明，算法计算精度准确，计算效率高，而且可变网格的变换对地质模型的研究起到了很重要的作用。 Abstract： Aiming at the typical gay geological model， the seismic wave equation and finite difference method are used to do simulation research to the staggered grid seismic wave field. In order to get a first order rate of seismic wave equation， the original form of the seismic wave equation is deduced firstly. Then space-time dual higher order difference staggered grid method is used to do discretization. Finally the geological model simulation experiment is done and the corresponding wave field snapshots are obtained. The simulation results show that the algorithm is accurate with high computational efficiency， and the varision of variable grid is of important significance for the study geological model. 关键词： 有限差分；地震模型；波动方程；正演模型 Key words： finite difference；seismic model；wave equation；forward model 中图分类号：O241.82 文献标识码：A 文章编号：1006-4311（2014）02-0315-02 0 引言 随着地震勘探面临的勘探任务越来越复杂，地震勘探成效在很大程度上，取决于对基于符合实际介质模型的方法理论研究。本文就详细介绍了有限差分法对交错网格地震波场进行的模拟仿真研究。一般的有限差分地震模拟方法是规则的交错网格。但是当地质体强纵横向变速介质时计算效率却很低。为此，人们研究了基于可变网格和不规则网格的地震波数值模拟方法。它对地质模型的离散化更为合理，同时，对于保持模型计算的灵活性也非常重要。从空间采样的角度考虑，最有效的提高模拟精度同时又降低计算机内存需求的方法，就是在模型的不同区域采用不同的网格步长，即可变网格（也称为不规则网格）[1]。 1 正演模拟方法原理 1.1 均匀各向同性完全弹性介质的弹性波动方程 由本质上讲，地震波就是在地下岩层中传播的机械弹性波，但反映弹性波传播的基本规律是由弹性波的波动方程[2]。由研究岩石的不同物理参数得到均匀各项同性完全弹性介质的波动方程，即 ■+■+■+ρg■=ρ■■+■+■+ρg■=ρ■■+■+■+ρg■=ρ■（1） 式中：t为时间分量；u，v，ω为介质位移在x，y，z三个方向上的分量；σxx、σyy、σzz、τxy、τyz、τzx为应力分量；gx、gy、gz为体力密度分量；ρ为介质密度。 1.2 各向同性介质中的一阶弹性波方程 弹性波理论是弹性体受力和应变的关系，根据各向同性介质表示的的本构方程和表示的柯西方程，可以推导出一阶速度—应力弹性波方程。二维二分量各向同性介质中的一阶应力—速度弹性波方程（假定体力为零）[4]，即 ρ■=■+■ρ■=■+■■=（λ+2μ）■+λ■■=λ■+（λ+2μ）■■=μ（■+■）（2） 式中λ、μ表示拉梅弹性常数，ρ表示弹性体密度，σ■、σ■表示正应力，σ■表示切应力，vx、vz分别表示x方向、z方向速度分量。 2 网格有限差分法 2.1 有限差分法基本原理 有限差分法（Finite Differential Method）是基于差分原理的一种数值计算方法，基本思想是：将波场域离散为许多小网格，应用差分原理，用差商来近似代替微商，将求解的连续函数问题