The frequency and group velocity dispersion behaviors, and characteristic surfaces of waves in a hybrid multilayered piezoelectric circular cylinder are investigated. The associated frequency dispersion equation is developed using an analytical-numerical method. In this method, the cylinder is modeled using the three-nodal-line layer element; the coupling between the elastic field and the electric field is considered in each element. A system of governing differential equations of each layer element is obtained following the Hamilton Principle. The phase velocity and slowness as well as the group velocity and slowness are established in terms of the Rayleigh quotient. Six characteristic wave surfaces, e.g. the phase velocity, slowness and wave surfaces as well as the group velocity, slowness and wave surfaces, are introduced to visualize the effects of anisotropy and piezoelectricity on wave propagation. A corresponding program code is developed and numerical examples are presented for hybrid multilayered piezoelectric circular cylinders with two ratios of radius to thickness.