Abstract Based on the linear traveltime interpolation (LTI) algorithm, we propose an improved LTI (referred as ILTI) algorithm for multiphase seismic ray tracing, which uses a velocity node in model parameterization and introduces secondary nodes between adjacent nodes. To better fit undulating surface model, irregular is used near interface, regular are still region far away from interface. We...