The production cross-sections of $Upsilon(1S)$, $Upsilon(2S)$ and $Upsilon(3S)$ mesons in proton-proton collisions at $sqrt{s}$= 13 TeV are measured with a data sample corresponding to an integrated luminosity of $277 pm 11$ $rm pb^{-1}$ recorded by the LHCb experiment in 2015. The $Upsilon$ mesons are reconstructed in the decay mode $Upsilontomu^{+}mu^{-}$. The differential production cross-sections times the dimuon branching fractions are measured as a function of the $Upsilon$ transverse momentum, $p_{rm T}$, and rapidity, $y$, over the range $0 < p_{rm T}< 30$ GeV/c and $2.0 < y < 4.5$. The ratios of the cross-sections with respect to the LHCb measurement at $sqrt{s}$= 8 TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.