This study uses iso-geometric investigation, which is based on the non-uniform rational B-splines (NURBS) basis function, to investigate natural oscillation of bi-directional functionally graded porous (BFGP) doubly-curved shallow microshells placed on Pasternak foundations with any boundary conditions. The characteristics of the present material vary in both thickness and axial directions along the x-axis. To be more specific, a material length-scale coefficient of the microshell varies in both thickness and length directions as the material’s mechanical properties. One is able to develop a differential equation system with varying coefficients that regulate the motion of BFGP double-curved shallow microshells by using Hamilton principle, Kirchhoff–Love hypothesis, and modified couple stress theory. The numerical findings are reported for thin microshells that are spherical, cylindrical, and hyperbolic paraboloidal, with a variety of planforms, including rectangles and circles. The validity and effectiveness of the established model are shown by comparing the numerical results given by the proposed formulations with previously published findings in many specific circumstances. In addition, influences of length scale parameters, power-law indexes, thickness-to-side ratio, and radius ratio on natural oscillation responses of BFGP microshells are investigated in detail.