A method for approximating the solution of a class of two-dimensional nonlinear Volterra integral equations is discussed. The properties of two-dimensional block-pulse functions are presented. The operational matrices of integration and product together with the collocation points are utilized to reduce the solution of the integral equation to the solution of a system of nonlinear algebraic equations. A numerical example is presented to show the efficiency of the proposed method