On numerical cubatures of singular surface integrals in boundary element methods