This method is applied for determination of the interval velocities from arrival times and relative arrival times of compression (P) and vertically (SV) and horizontally (SH) polarized shear (S) seismic waves which are generated near surface and travel down to an array of vertically installed seismic sensors