摘要: 对于大多数速度场，地震波沿射线传播的初至波走时，可以用有限差分外推的方法在二维或三维数值网格上计算出来. 在保证精度的条件下，为提高计算效率和适应性，本文推导了基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法. 另外，该方法对首波和散射波做了合适的处理，而且不会碰到传统射线法存在的阴影区和焦散区等问题. 简单模型和复杂的Marmousi模型试算的结果表明，该方法精度较高并适用于强纵、横向变速的复杂介质. 基于该方法的Kirchhoff叠前深度偏移, 在主要构造和目的层位置的成像效果上基本达到了波动方程法叠前深度偏移的位置成像效果. 由于未考虑续至波等有效能量，在成像的保幅性上不如波动方程法叠前深度偏移的效果，但其计算效率则明显高于全格林函数法和波动方程法.