Prediction of Z→l+l- production cross section (where l±=e±,μ±) in proton-proton collisions at s=14 TeV is estimated up to next-to-next-to-leading order (NNLO) in perturbative QCD including next-to-leading order (NLO) electroweak (EW) corrections. The total inclusive Z boson production cross section times leptonic branching ratio, within the invariant mass window 66<mll<116 GeV, is predicted using NNLO HERAPDF2.0 at NNLO QCD and NLO EW as σZTot=2111.69-26.92+26.31 (PDF) ±11 (αs) ±17 (scale) -30.98+57.41 (parameterization and model). Theoretical prediction of the fiducial cross section is further computed with the latest modern PDF models (CT14, MMHT2014, NNPDF3.0, HERAPDF2.0, and ABM12) at NNLO for QCD and NLO for EW. The central values of the predictions are based on DYNNLO 1.5 program and the uncertainties are extracted using FEWZ 3.1 program. In addition, the cross section is also calculated as functions of μR and μF scales. The choice of μR and μF for scale variation uncertainty is further discussed in detail.