Dans cet exposé nous considérons les équations quasi-geostrophiques surfaciques et nous proposons des schémas multi-pas en temps pour ce modèle dans le cadre de la formulation stochastique sous incertitude de position (LU, location uncertainty). La modélisation stochastique LU est basée sur une décomposition de la vitesse lagrangienne en deux composantes : une composante lisse de grande échelle et une stochastique à petite échelle. Dans ce contexte, la généralisation du schéma d’Euler a une convergence forte d’ordre 0.5 et faible d’ordre 1. Ces convergences de bas ordre entraînent des contraintes assez restrictives sur le pas de temps à utiliser. Nous proposons un schéma de type Milstein, qui améliore la convergence forte. Ce schéma est ensuite utilisé dans un schéma multi-pas.