A walkaway VSP survey is carried out using a receiver array. First arrivals to a plurality of receivers are picked and used to estimate a normal-moveout (NMO) velocity. Using the NMO velocity and vertical velocities estimated from the VSP data, two anisotropy parameters are estimated for each of the layers. The anisotropy parameters may then be used to process surface seismic data to give a stacked image in true depth and for the interpretation purposes. For multi-azimuthal walkaway or 3D VSP data, we determine two VTI parameters ε and δ for multi-azimuth vertical planes. Then we determine five anisotropic interval parameters that describe P-wave kinematics for orthorhombic layers. These orthorhombic parameters may then be used to process surface seismic data to give a stacked image in true depth and for the interpretation purposes.