“…For example, the initial value ray tracing is achieved by fitting the rolling surface with cubic spline interpolation function (Yue et al, 2007), the fast seismic ray tracing is achieved by building complex geological models with block structures (Xu at al., 2004;Wu and Jia., 2008), travel time calculations in complex rolling-layered medium is achieved by using irregular grids (trapezoid grids and triangular grids) to fit the rolling surface as well as partitioned multistep computation technique (Huang and Bai, 2010;Zhao and Bai, 2010), and so on. Most researches on seismic ray tracing in anisotropic media are still based on simple layered models (Qin and Schuster, 1993;Zhou and Greenhalgh, 2006;Jiang and Zhou, 2011), and only a few researches (Li at al., 2010) take the rolling surface into consideration.…”