A computer-implemented method for pseudo acoustic quasi-P wave propagation which remain stable in anisotropic media with variable tilt and is not limited to weak anisotropic conditions. The method includes acquiring a seismic exploration volume for a subsurface region of interest, and determining a modeling geometry for the seismic exploration volume. The method further includes propagating at least one wavefield through the seismic exploration volume utilizing the modeling geometry and initial conditions and preventing the accumulation of energy along the axis of symmetry of the seismic exploration volume and ensuring positive stiffness coefficients in the stress-strain relations through the use of finite quasi-S wave velocities thereby producing a stable wavefield. The method includes utilizing the stable wavefield to generate subsurface images of the subsurface region of interest.