Anisotropy generally exists in underground media. Nowadays, seismic wave ray tracing based on anisotropy is mostly carried out in weakly anisotropic media by using the approximate formula of group velocity, which will lead to large errors in strongly anisotropic media, and can not really simulate the propagation law of seismic wave. According to the anisotropy in underground media and the basic propagation law of seismic wave, this paper proposed the Newton iteration method for solving the group velocity. Based on the Paraview platform, a 3D geological model was automatically constructed, and the shortest path method was used to simulate and visualize the seismic ray tracing. Taking the northern Shanxi rift zone in the North China Craton as an example, the expression of velocity heterogeneity and anisotropy in complex 3D geology was realized. The experimental results showed that the proposed method can mitigate the influence of anisotropy on seismic wave propagation simulation, clearly express the geological structure and anisotropic characteristics of the study area, so it can be better applied in the interpretation of complex 3D geological structure.