This paper presents the construction of a ray tracing system adapted for any kind of curvilinear coordinate and the derivation of its main formulas in details. The applications of the ray tracing system on spherical coordinate are given and the ray tracing results of theoretical and real 3D model are discussed. Comparing with the bending method, the 3D ray tracing on spherical coordinate can describe the velocity distribution of medium and the seismic wave propagation more precisely.