In this article, we apply the operational matrix to find the numerical solution of two- dimensional nonlinear Volterra integro-differential equation (2DNVIDE). Form this prospect, two-dimensional shifted Legendre functions (2DSLFs) has been presented for integration, product as well as differentiation. This method converts 2DNVIDE to an algebraic system of equations, so the numerical solution of 2DNVIDE is computable. The effectiveness and accuracy of the method were examined with some examples as well. The results and comparison with other methods have shown a remarkable performance.