Paraxial ray approximations in the computation of seismic wavefields in inhomogeneous media