Wir beschreiben numerische Verfahren zur Lösung stationärer nichtlinearer Probleme mit stochastischen Unsicherheiten im Operator, in den Randbedingungen und in den Lasten. Dabei vergleichen wir hochdimensionale (Smolyak-) Quadraturverfahren mit Monte Carlo-Integrationsverfahren zur direkten Berechnung von Statistiken der Lösung. Zusätzlich berechnen wir die Lösung mit einem Galerkin-Verfahren in einem Raum stochastischer Ansatzfunktionen und untersuchen die Auswertung der im Residuum auftretenden hochdimensionalen Integrale. Für das entstehende große gekoppelte nichtlineare Gleichungssystem stellen wir einen effizienten Löser vor. Schließlich berechnen wir die Lösung durch direkte Projektion auf den orthogonalen stochastischen Ansatzraum und vergleichen die beschriebenen Lösungsverfahren anhand von Modellproblemen. Die Verfahren wenden wir auf ein prototypisches Grundwasserfluss-Problem an.