In this study, the Bernoulli polynomials are used to ob tain an approximate solution of a class of nonlinear two-dimensional integral equations. To this aim, the operational matrices of inte- gration and the product for Bernoulli polynomials are derived and utilized to reduce the considered problem to a system of nonlinear algebraic equations. Some examples are presented to illustrate the eciency and accuracy of the method.