一种基于异距21点拉伸型网格的频率域正演模拟方法

本发明公开了一种基于异距21点拉伸型网格的频率域正演模拟方法，涉及石油地球物理勘探技术领域，输入速度场、频域震源子波、观测系统参数；推导DD21点差分格式；求取系数优化，并进行频散分析；施加PML边界条件，构建DD21点格式PML波动方程；基于DD21点格式PML波动方程构造阻抗矩阵；求解阻抗矩阵，计算频率‑空间域下的波场值；将频率‑空间域下的波场变换到时间‑空间域下；抽取检波点波场，输出频率域正演模拟结果。本方法更适用于矩形网格采样的异距21点差分格式，当网格间距比大于1.5时，每个波长仅需2.13个网格点即可将相速度误差控制在1％以内，提高了频率域正演模拟的运算效率和精度，为后续的反演及成像处理提供了一种更高效的算法。