In this paper, a novel size-dependent functionally graded (FG) cylindrical shell model is developed based on the nonlocal strain gradient theory in conjunction with the Gurtin-Murdoch surface elasticity theory. The new model containing a nonlocal parameter, a material length scale parameter, and several surface elastic constants can capture three typical types of size effects simultaneously, which are the nonlocal stress effect, the strain gradient effect, and the surface energy effects. With the help of Hamilton’s principle and first-order shear deformation theory, the non-classical governing equations and related boundary conditions are derived. By using the proposed model, the free vibration problem of FG cylindrical nanoshells with material properties varying continuously through the thickness according to a power-law distribution is analytically solved, and the closed-form solutions for natural frequencies under various boundary conditions are obtained. After verifying the reliability of the proposed model and analytical method by comparing the degenerated results with those available in the literature, the influences of nonlocal parameter, material length scale parameter, power-law index, radius-to-thickness ratio, length-to-radius ratio, and surface effects on the vibration characteristic of functionally graded cylindrical nanoshells are examined in detail.