Extra neutral gauge bosons (Z) are predicted in many extensions of the Standard Model (SM). In the minimal anomaly-free Z model (AFZ), the phenomenology is controlled by only three parameters beyond the SM ones, the Z mass and two effective coupling constants g_Y and g_{BL}. We study the Z 5-sigma discovery potential in e+e- collisions at 1.4 and 3 TeV CLIC. Assuming LHC discovers a Z of 5 TeV mass, the expected accuracies on the Zmu+mu- couplings are presented. We discuss also the requirements on detector performance and beam polarization.