1.一种基于初至波走时约束的动态波场传播模拟方法，其特征在于，包括以下步骤：A、将给定的速度模型进行网格化，设置模型设置横向网格数nx，纵向网格数nz，以及相应的网格间距Δx和Δz，得到离散速度模型v(x,z)，x∈[0,nx],z∈[0,nz]；

B、设置差分阶数M和差分系数

C、将离散的速度模型v(x,z)代入公式(2)中的程函方程，计算初至波走时T(x,z)，初至波走时是一系列等时线；

D、将初至波走时T(x,z)中每一条等时间线上的点对应的空间坐标x和z通过time-bin技术按时间增加方向重新排序，封装到每一时刻的动态结构体数组D[c(t)]中，其中，c(t)代表对应t时刻的等时线上的点数；

E、定义基于初至波走时约束的动态模拟范围Ω

F、对波动方程进行离散，构建波动方程高阶有限差分格式，如公式(5)，其中，G、利用离散波动方程的高阶有限差分格式，逐一时刻在动态的空间窗口范围ΩH、随着传播时刻增加，当波场传播至模型边界时，加载PML边界条件进行处理，消除由于有限的模型边界产生的反射。

2.根据权利要求1所述的一种基于初至波走时约束的动态波场传播模拟方法，其特征在于：步骤E，定义基于初至波走时约束的动态模拟范围，具体包括以下步骤：对每个时刻对应的D[c(t)]中的空间坐标x和z分别进行比较，确定每一个时刻数值模拟的动态窗口边界。

3.根据权利要求2所述的一种基于初至波走时约束的动态波场传播模拟方法，其特征在于：对于任意时刻t，最小的x坐标赋值给模拟窗口的左边界