用于抑制数值频散的八阶NAD-RK方法、系统、设备及介质技术方案

本发明专利技术属于地震波数值模拟技术领域，公开了一种用于抑制数值频散的八阶NAD

【技术实现步骤摘要】
用于抑制数值频散的八阶NAD
‑
RK方法、系统、设备及介质


[0001]本专利技术属于地震波数值模拟领域，尤其涉及一种一种用于抑制数值频散的八阶NAD
‑
RK方法、系统、设备及介质。

技术介绍

[0002]目前，开发具有高效率和高精度的地震波数值模拟方法是当前研究地震波动方程反演的重要课题。传统的数值模拟方法由于存在严重的数值频散现象而不能满足当前大尺度地震波模拟的要求，杨顶辉教授团队2003年首次在地震波正演模拟中引入近似解析离散化算子(NAD)，目前已获得系列NAD算子的数值模拟方法，这些方法均有着良好的数值频散效果，但是这些方法仅为四阶的空间离散精度。
[0003]通过上述分析，现有技术存在的问题及缺陷为：
[0004]现有方法的空间离散精度较低。

技术实现思路

[0005]针对现有技术存在的问题，本专利技术提供了一种用于抑制数值频散的八阶NAD
‑
RK方法、系统、设备及介质。
[0006]本专利技术是这样实现的，一种用于抑制数值频散的八阶NAD本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种用于抑制数值频散的八阶NAD
‑
RK方法，其特征在于，包括：将二维声波传播方程转为常微分方程，利用八阶NAD算子对所述常微分方程的空间高阶偏导数进行离散化处理，再采用三阶Runge
‑
Kutta方法对所述常微分方程的时间高阶偏导数进行离散化处理，得到八阶NAD
‑
RK方法；并通过比较波形记录分析所述八阶NAD
‑
RK方法的数值频散关系。2.如权利要求1所述用于抑制数值频散的八阶NAD
‑
RK方法，其特征在于，所述二维声波传播方程为：式中，u表示位移，c表示声波速度；令则得到以下向量方程：式中，三阶空间微分算子A定义为其中I3×3是三阶单位矩阵；再令V＝(U,W)
T
，则向量方程转换为常微分方程，表达式为：式中，D表示空间偏微分算子矩阵，且3.如权利要求1所述用于抑制数值频散的八阶NAD
‑
RK方法，其特征在于，所述八阶NAD算子离散化处理的具体过程为：所述高阶偏微分算子矩阵D中包含位移u和粒子速度w的二阶和三阶偏导数，对所述常微分方程中的高阶偏微分算子矩阵D进行空间高阶偏导数离散化，将所述常微分方程转化为关于时间t的离散化常微分方程；所述三阶Runge
‑
Kutta方法对所述半离散化常微分方程进行时间高阶偏导数离散化，计算公式为：式中，Δt为时间步长，V
(1)
、V
(2)
为中间变量，V
n
＝V(nΔt)；消去公式中的中间变量V
(1)
、V
(2)
，得到最终计算公式为：
4.如权利要求1所述用于抑制数值频散的八阶NAD
‑
RK方法，其特征在于，所述八阶NAD
‑
RK方法的计算公式为：RK方法的计算公式为：这里A2＝A
·
A。5.如权利要求1所述用于抑制数值频散的八阶NAD
‑
RK方法，其特征在于，所述八阶NAD
‑
RK方法的数值频散关系为：在步长h＝Δx＝Δz的情况下，Δx＝Δz表示网格步长，所述三阶Runge
‑
Kutta方法的最终计算公式的谐波解为：式中：式中：将所述谐波解代入所述最终计算公式并结合高阶偏导数的离散表达式得到数值频散方程，表达式为：Det(e
iγ
I6‑
H)＝0式中，γ＝ω
num
Δt(ω
num
为数值角速度，...

【专利技术属性】
技术研发人员：陈丽，张朝元，
申请(专利权)人：大理大学，
类型：发明
国别省市：