Numerical solution of surface integral equations based on spline quasi-interpolation