一种“M(2，4)FDTD+FDTD”的低色散低频地波传播时延预测方法

摘要

本发明公开了一种“M(2，4)FDTD+FDTD”的低色散低频地波传播时延预测方法，该方法通过对整个计算区域进行分层处理，将粗网格M(2，4)FDTD算法(用于上层区域)与细网格传统FDTD算法(用于下层区域)相结合，利用修正因子以及亚网格技术，进行低频地波传播时延的高精度快速预测，本发明解决了传统全区域FDTD算法在减小数值色散误差与降低计算机消耗之间的矛盾，提高了低频地波传播时延预测精度的同时，减少了计算机内存资源占用，提高了计算速度。

主权项

一种“M(2，4)FDTD+FDTD”的低色散低频地波传播时延预测方法，其特征在于，通过对低色散低频地波计算区域进行分层处理，将粗网格M(2，4)FDTD与细网格传统FDTD方法相结合，进行低频地波传播时延的预测，从而达到在保证预测精度的同时，提高计算速度、减少计算机内存占用，具体步骤如下：步骤1、设置并生成模型文件、输入模型文件；步骤2、参数设置及初始化；步骤3、添加场源；步骤4、更新整个计算区域的ρ向电场分量E_ρ；步骤5、更新整个计算区域的z向电场分量E_z；步骤6、更新整个计算区域的φ向磁场分量H_φ；步骤7、更新激励源；步骤8、判断结束条件，循环；步骤9、观测点观测量计算并输出。