Quasi-interpolation based on the ZP-element for the numerical solution of integral equations on surfaces in $mathbb R^3$