Nous présentons les outils de la modélisation de fiabilité des remblais homogènes avec une représentation de la variabilité spatiale des propriétés des sols par des champs stochastiques spatiaux en introduisant l'autocorrélation horizontale suivant les deux axes x et y, ainsi que l'autocorrélation verticale suivant l'axe z de la cohésion non drainée. La répartition spatiale de la cohésion non drainée a été modélisée au moyen d'un champ stochastique. Cette approche nous a permis d'évaluer l'espérance et la variance de la valeur moyenne de la résistance au cisaillement le long d'une surface potentielle de rupture dans le cadre d'une analyse limite tridimensionnelle de stabilité d'un remblai homogène (Hovland, 1977). Cette approche nous a permis de montrer que la fiabilité des mécanismes 3D peut être plus faible que celle des mécanismes 2D bien que le facteur de sécurité déterministe soit plus fort, du fait que la variance de la cohésion moyenne sur une surface de rupture de faibles dimensions reste forte. L'application de cette approche pour le cas du barrage de Mirgenbach avec une surface de rupture comprenant une partie plane (bloc), bien qu'elle corresponde à un facteur de sécurité déterministe plus élevé, peut être moins fiable que celle d'une surface de rupture de type ellipsoïde de glissement, du fait que la variance de la cohésion moyenne sur la partie horizontale reste très forte en raison de la forte corrélation horizontale. Ceci a été montré par Auvinet et Rossa (1991) dans le cas 2D sur le même ouvrage.