Zr-L2,3 XANESs of tetravalent zirconium oxides with different coordination numbers and local symmetries are systematically investigated by ab initio multiplet calculations using fully relativistic molecular spinors for model clusters. Experimental Zr-L2,3 XANESs are obtained for SrZrO3, m-ZrO2 (monoclinic) and t-ZrO2 (tetragonal). The theoretical spectra are in good agreement with the experimental data. The multiplet effects are found to play essential roles in determining the peak shape. The shapes of L3- and L2-edges are systematically different. The intensity ratios of the doublet peaks at both L3- and L2-edges are found to be sensitive to the coordination number of Zr. The ratio can therefore be used to estimate the coordination number of Zr in such oxides.