In this paper, we propose a method to deal with numerical integral
by using two kinds of $C^2$ quasi-interpolation operators on the
bivariate spline space, and also discuss the convergence properties
and error estimates. Moreover, the proposed method is applied to the
numerical evaluation of 2-D singular integrals. Numerical
experiments will be carried out and the results will be compared
with some previously published results.