In this paper, we study convergence of a fully discrete scheme for the two-dimensional nonstationary Kelvin-Voigt model. This scheme is based on a finite element approximation for space discretization and the Crank-Nicolson-type scheme for time discretization, which is a two step method. Moreover, we obtain error estimates of velocity and pressure. At last, the applicability and effectiveness of the present algorithm are illustrated by numerical experiments