Numerical formulations based on surface integral equations (SIEs) provide an accurate and efficient framework for the solution of the electromagnetic scattering problem by three-dimensional plasmonic nanostructures in the frequency domain. In this paper, we present a unified description of SIE formulations with both singular and nonsingular kernel and we study their accuracy in solving the scattering problem by metallic nanoparticles with spherical and nonspherical shape. In fact, the accuracy of the numerical solution, especially in the near zone, is of great importance in the analysis and design of plasmonic nanostructures, whose operation critically depends on the manipulation of electromagnetic hot spots. Four formulation types are considered: the N-combined region integral equations, the T-combined region integral equations, the combined field integral equations and the null field integral equations. A detailed comparison between their numerical solutions obtained for several nanoparticle shapes is performed by examining convergence rate and accuracy in both the far and near zone of the scatterer as a function of the number of degrees of freedom. A rigorous analysis of SIE formulations and their limitations can have a high impact on the engineering of numerous nano-scale optical devices such as plasmon-enhanced light emitters, biosensors, photodetectors, and nanoantennas.