A surface integral equation (SIE) formulation is applied to the analysis of electromagnetic problems involving three-dimensional (3D) piecewise homogenized left-handed metamaterials (LHM). The resulting integral equations are discretized by the well-known method of moments (MoM) and solved via an iterative process. The unknowns are defined only on the interfaces between different media, avoiding the discretization of volumes and surrounding space, which entails a drastic reduction in the number of unknowns arising in the numerical discretization of the equations. Besides, the SIE-MoM formulation inherently includes the radiation condition at infinity, so it is not necessary to artificially include termination absorbing boundary conditions. Some 3D numerical examples are presented to confirm the validity and versatility of this approach on dealing with LHM, also providing some intuitive verifications of the singular properties of these amazing materials.