Due to some technical issues that can appear during the manufacturing process of Functionally Graded Materials (FGMs), it can be extremely difficult to produce perfect materials. Indeed, one of the biggest problems is the presence of porosities. For this purpose, the vibrational behavior of doubly-curved shells made of FGM including porosities is investigated in this paper. With respect to previous research, the porosity has been added to the mechanical model that characterizes the through-the-thickness distribution of the graded constituents and applied to doubly-curved shell structures. Few papers have been published on this topic. In fact, it is easier to find works related to one-dimensional structures and beam models that take account the effect of porosities. The First-order Shear Deformation Theory (FSDT) is considered as the theoretical framework. In addition, the mechanical properties of the constituents vary along the thickness direction. For this purpose, two power-law distributions are employed to characterize their volume fraction. Strain components are established in an orthogonal curvilinear coordinate system and the governing equations are derived according to the Hamilton's principle. Finally, Navier's solution method is used and the numerical results concerning three different types of shell structures are presented.