Dans le domaine des barrages, la pratique classique de l'ingénierie considère une seule cote maximale de remplissage de la retenue, la cote des plus hautes eaux (PHE), qui fixe l'intensité maximale de la sollicitation hydraulique. Cet article présente une méthode probabiliste pour modéliser la cote de remplissage de la retenue. Elle met en ½uvre de multiples hydrogrammes de crue générés aléatoirement par simulations et elle prend en compte la variabilité des cotes de remplissage de la retenue. L'application de la méthode à trois barrages en service permet d'apprécier son intérêt dans le cadre d'une utilisation en analyse de sureté de fonctionnement.