An algorithm to solve the inverse problem of detecting inclusion interfaces in a piezoelectric structure is proposed. The material interfaces are implicitly represented by level sets which are identified by applying regularization using total variation penalty terms. The inverse problem is solved iteratively and the extended finite element method is used for the analysis of the structure in each iteration. The formulation is presented for three-dimensional structures and inclusions made of different materials are detected by using multiple level sets. The results obtained prove that the iterative procedure proposed can determine the location and approximate shape of material sub-domains in the presence of higher noise levels.