In this paper, two efficient two-grid finite element algorithms are proposed for solving two-dimensional nonlinear pseudo-parabolic integro-differential equations. Firstly, we obtain optimal error estimates of the fully discrete finite element method using a temporal-spatial error splitting technique in H1 and Lp norms. Then the two-grid technique to improve computation efficiency of the proposed finite element method. Error estimates in H1 and Lp norms of two-grid solutions are presented. Theoretical analysis shows that the two-grid algorithms maintain asymptotically optimal accuracy. Finally, numerical examples are provided to support our theoretical results and demonstrate the effectiveness of these methods.