Global analysis of the Sivers functions at NLO+NNLL in QCD