1.一种有限差分数值模拟方法,其特征在于,所述有限差分数值模拟方法包括:采用十字和菱形交错网格有限差分模板逼近时间偏导数;
采用隐式交错网格有限差分法逼近空间偏导数;
频散关系推导;
差分系数求取;
得到的时间高阶差分系数和隐式空间差分系数分别代入公式进行弹性波方程数值模拟,实现基于纵横波分离的P‑SV波高精度模拟;
所述十字和菱形交错网格有限差分模板逼近时间偏导数方法包括:二维弹性波速度应力方程表示为:
其中,(vx,vz)为偏振速度矢量,λ和μ为拉梅系数,ρ为密度;采用十字和菱形交错网格有限差分模板逼近时间偏导数,公式如下:其中,
其中,dm,0和dm,n为时间差分系数,Δt为时间步长,h为空间采样间隔,M和N为差分算子长度参数。
2.如权利要求1所述的有限差分数值模拟方法,其特征在于,所述隐式交错网格有限差分法逼近空间偏导数方法包括:隐式交错网格有限差分逼近空间偏导数的公式为:其中,cm和b为空间差分系数,
隐式交错网格有限差分中,空间偏导数通过对三对角线性方程组求解得到。
3.如权利要求2所述的有限差分数值模拟方法,其特征在于,所述频散关系推导方法包括:平面波理论下,波场表示公式为:
其中, 为虚数单位,k为波数,ω为角频率,θ为传播角度,kx=kcos(θ),kz=ksin(θ),将波场表示公式代入公式 并化简得到:其中,rP=vPΔt/h和rS=vSΔt/h为纵波和横波的库朗数,vP和vS分别为纵波和横波的传播速度,公式 和 为纵波
和横波的时间频散关系。
4.如权利要求3所述的有限差分数值模拟方法,其特征在于,所述差分系数求取方法包括:将公式 中的三角函数进行Taylor级数展开化简可得:
C+D=E,
其中,
其中,
2l 2l
比较方程C+D=E两端(kxh) 或(kzh) 的系数得:其中1≤l≤M比较方程C+D=E两端 的系数得,其中2≤l≤N,1≤ξ≤l‑1其中,int(x)为取整函数,通过公式:和
得到基于
Taylor级数展开的纵波时间高阶差分系数,将vP换为vS得到横波时间高阶差分系数;
隐式空间交错网格有限差分系数求取公式为:
5.如权利要求4所述的有限差分数值模拟方法,其特征在于,所述基于纵横波分离的P‑SV波高精度模拟方法包括:将第四步中得到的时间高阶差分系数和隐式空间差分系数分别代入公式 和进行弹性波方程数值模拟,采用波数域分离方法得到高精度的P‑SV转换波;其中;
模拟步骤为:
(1)采用纵波时间高阶差分系数和隐式空间差分系数求解弹性波方程,纵横波分离,只提取纵波场;
(2)采用横波时间高阶差分系数和隐式空间差分系数求解弹性波方程,纵横波分离,只提取横波场;
(3)将(1)中的纵波场和(2)中的横波场求和,得到最终的弹性波场;
(4)重复(1)‑(3),直到最大时刻。
6.一种计算机设备,其特征在于,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行权利要求1~
5任意一项所述的有限差分数值模拟方法。
7.一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,使得所述处理器执行权利要求1~5任意一项所述的有限差分数值模拟方法。
8.一种实施权利要求1~5任意一项所述有限差分数值模拟方法的有限差分数值模拟系统,其特征在于,所述有限差分数值模拟系统包括:时间偏导数获取模块,用于采用十字和菱形交错网格有限差分模板逼近时间偏导数;
空间偏导数获取模块,用于采用隐式交错网格有限差分法逼近空间偏导数;
频散关系获取模块,用于实现频散关系推导;
差分系数求取模块,用于差分系数求取;
数值模拟模块,用于基于纵横波分离实现P‑SV波高精度数值模拟。
9.一种如权利要求1~5任意一项所述有限差分数值模拟方法在弹性介质地震波正演、逆时偏移和全波形反演中的应用。