[发明专利]一种在复杂各向异性介质中模拟地震波的方法及装置

摘要

本发明提供了一种在复杂各向异性介质中模拟地震波的方法及装置，方法包括：给数值化的三维地质模型的每个网格点的弹性常数和密度赋值；获取三维地质模型中网格点沿三个坐标轴的速度分量的一阶空间导数；将速度分量的一阶空间导数代入到弹性介质中弹性波方程组中的第一方程式，并进行求解获取应力分量；获取三维地质模型中网格点沿三个坐标轴的应力分量的一阶空间导数；将应力分量的一阶空间导数代入到弹性介质中弹性波方程组中的第二方程式，并进行求解获取速度分量；根据通过第二方程式获取的速度分量得到当前时间第一时刻的速度分量的值，作为当前时刻的地震波的记录；获取下一时刻的地震波的记录，直至获取时间长度t的地震波的波场值。

主权项

一种在复杂各向异性介质中模拟地震波的方法，其特征在于，所述方法包括：给数值化的地质模型的每个网格点的弹性常数c和密度ρ赋值；通过获取三维地质模型中在网格点坐标为（i，j，k）的沿三个坐标轴的速度分量的一阶空间导数；其中，当网格点的坐标i+j+k=偶数，Ti,j,k=v；当网格点的坐标i+j+k=奇数，Ti,j,k的速度分量为网格点（i，j，k）的上、下、左、右、前、后的六个坐标值为偶数的网格点的速度分量的平均值，我们称这种网格为菱形交错网格；将获取的所述速度分量v的一阶空间导数代入到弹性介质中弹性波方程组中的第一方程式，并对第一方程式进行求解获取应力分量σ；其中，应力分量σ为速度分量v的表达式；利用所述应力分量σ，根据Ti,j,k=σor c orρ,(i+j+k)%2=1v,(i+j+k)%2=0,]]>通过获取三维地质模型中在网格点坐标为（i，j，k）的沿三个坐标轴的应力分量σ的一阶空间导数；将获取的所述应力分量σ的一阶空间导数代入到弹性介质中弹性波方程组中的第二方程式，并对第二方程式进行求解获取速度分量v；根据通过第二方程式获取的速度分量得到当前时刻的速度分量的值，作为当前时刻的地震波的记录；获取下一时刻的地震波的记录，直至获取时间长度t的地震波的波场值。