Abstract: Using the rotated staggered-grid of finite difference form, the numerical modeling of the first-order stress-velocity equation of quasi-P wave in VTI (vertical axis of symmetry of the transversely isotropic) medium is achieved. Under the PML boundary conditions and the stability condition, high precision wave field snapshots and synthetic record of complex models such as Marmousi are obtained. Besides, we analyze the effect of the anisotropy on seismic wave. The results show that the rotated staggered-grid of finite difference can get high precision seismic simulation data, and the PML boundary has better absorption effect.