On s'intéresse à l'approximation numérique de la solution du couplage entre l'équation d'écoulement dans un milieu poreux aléatoire et l'équation d'advection-diffusion. Plus précisément, on présente et analyse une méthode numérique pour calculer la valeur moyenne de l'extension d'un soluté introduit au temps initial, et la valeur moyenne de la macro-dispersion, définie comme la dérivée temporelle de l'extension. On propose une méthode de Monte-Carlo pour tenir compte des incertitudes, c'est à dire du caractère aléatoire du champ de perméabilité. L'équation d'écoulement est alors résolue en utilisant des éléments finis. L'équation d'advection-diffusion est vue comme une équation de Fokker-Planck, et sa solution est donc approchée grâce à une méthode particulaire probabiliste. L'extension peut en effet être exprimée comme l'espérance d'une fonction de la solution de l'équation différentielle stochastique correspondante, et est calculée grâce à un schéma d'Euler pour l'équation différentielle stochastique et une méthode de Monte-Carlo. On donne des estimations d'erreur pour l'extension moyenne et la dispersion moyenne, sous différentes hypothèses, en particulier sur le champ de perméabilité.