In this paper we utilize affine biquadratic elements and a two-step temporal discretization to develop a finite volume element method for parabolic problems on quadrilateral meshes. The method is proved to have an optimal order convergence rate in L2(0, T; H1(Ω)) under the ‘asymptotically parallelogram’ mesh assumption. Numerical experiments that corroborate the theoretical analysis are also presented.