交错网格有限差分正演

人工智能-机器学习-高阶交错网格有限差分地震波场计算.pdf

optimum_FD MATLAB脚本可为交错网格有限差分计算最佳有限差分系数 由Erik Koene撰写，2018年5月17日版权归各自的作者和出版商所有。

INT（E）NSE是IN插值的二维（D）-Navier-Stokes-Equations的缩写（双关语）。 这是我在业余时间开发的不稳定的Navier-Stokes解算器。 使用混合的二阶和一阶精确方案对面部速度进行插值，以便能够捕获涡旋。 使用SIMPLE算法实现了压力-速度耦合，并且随着向后有限差分离散化，时间步长是一阶精确的。

以瑞利波频散方程为出发点的Abo-Zena传递矩阵等方法,只能研究层状介质中瑞利波的传播特性,对于非层状介质,传递矩阵无能为力。因此,不得不考虑有限元、有限差分等方法。根据弹性动力学方程,采用交错网格有限差分方法对均匀弹性半空间介质进行全波场模拟,并在Matlab环境下实现编程计算,再现了瑞利波在近地表的传播状态,实现了地震剖面曲线的正振幅充填问题,从而获得了类似野外地震勘探的剖面记录。模拟结果表明,采用一阶差分格式所得到的地震记录存在较严重的数值频散,在采用有限差分法研究瑞利波"之"字形成因时,须尽量压制数值频散现象,才能获得更接近实际情况的频散曲线。