[发明专利]双相各向异性介质交错网格有限差分模拟方法和装置

说明书

技术领域

本发明涉及计算数据处理技术领域，具体而言，涉及双相各向异性介质交错网格有限差分模拟方法和装置。

背景技术

目前，为了研究地球内部双相各向异性介质的特点，对储层介质加深认识，需对双相各向异性介质做一系列的正演模拟。经常采用交错网格对地震波在双相各向异性介质中的传播进行模拟。

相关技术中，采用交错网格对地震波在双相各向异性介质中的传播进行模拟的过程中，将网格点分为整网格点和半网格点，将质点振动速度和所受应力分别定义在两个不同且相邻的时间层上，通常将所有正应力放置于整网格点处，而切应力放置于半网格点处。利用相邻的整网格点和半网格点的场量值来分别计算速度分量和应力分量沿水平分量和竖直分量的差分，用差分代替微分，同时相邻两个时间层上的物理量在空间分布上恰好交错半个网格，以达到时间和空间交错的目的。

在通过交错网格对地震波在双相各向异性介质中的传播进行模拟的过程中，必须通过对切应力分量进行插值，才可以对地震波在双相各向异性介质中的传播进行模拟，插值计算本身并不精确，所以会增大模拟过程的计算误差，降低了得到的双相各向异性介质的介质质点的速度-应力数据的精度。

发明内容

有鉴于此，本发明实施例的目的在于提供双相各向异性介质交错网格有限差分模拟方法和装置，在通过交错网格对地震波在双相各向异性介质中的传播进行模拟的过程中，无需对切应力使用插值计算，以提高得到的双相各向异性介质质点的速度-应力数据的精度。

第一方面，本发明实施例提供了一种双相各向异性介质交错网格有限差分模拟方法，包括：

根据所述双相各向异性介质的运动方式，将所述双相各向异性介质的速度和应力分量分别设置在交错网格的网格点处，其中，将双相各向异性介质的介质质点的切应力分量设置在交错网格中的整网格点处；

根据确定的介质质点在所述交错网格中的速度和应力分量位置，计算速度和应力分量的离散化数据；

通过得到的所述交错网格对地震波在双相各向异性介质中的传播进行模拟，得到所述双相各向异性介质的介质质点的速度-应力数据。

结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，将双相各向异性介质的介质质点的切应力分量设置在交错网格中的整网格点处，包括：

通过公式