A three-dimensional phase-field model was developed to study the effects of surface energy and diffusivity anisotropy on void growth behavior in irradiated Zr. The gamma surface energy function used in the phase-field model was developed with the surface energy anisotropy calculated from the molecular dynamics simulations. The model assumes that vacancies have larger mobility in the c-axis than in the a-axis and b-axis, whereas interstitials have larger mobility in the basal plane than in the c-axis. The equilibrium void morphology and the effect of defect concentrations and mobility anisotropy on void growth behavior were simulated using the developed model. The simulations demonstrated that (1) the developed phase-field model correctly reproduces the faceted void morphology predicted by the Wulff construction; (2) with isotropic diffusivity, the void prefers to grow on the basal plane; and (3) when the vacancy has large mobility along the c-axis and interstitial has a large mobility on the basal plane of hexagonal closed-packed Zr alloys, a platelet void grows in the c-direction and shrinks on the basal plane. These findings are in agreement with the observation of void morphology and growth behavior in irradiated Zr, and reveal that the strong mobility anisotropy of defects results in the anisotropic growth of voids in Zr.