Le problčme d'écoulement biphasique est modélisé par un systčme d'équations de conservation de masse en saturation-pression-vitesse et le problčme de transport réactif par un systčme d'équations de conservation de concentration de type advection- diffusion-réaction. On applique la méthode des éléments finis mixtes étendus pour approcher simultanément les pressions, les saturations et les champs de vitesse du modčle biphasique. Le schéma utilisé est totalement implicite et non linéaire. Il est résolu par la méthode de Newton inexacte. On utilise les espaces de Raviart-Thomas et des formules de quadratures pour obtenir un schéma équivalent aux différences finies centrées. Pour le transport réactif, on applique la méthode de splitting d'opérateurs avec une discrétisation temporelle semi-implicite, un schéma de Godunov pour approcher le terme convectif, sous la condition de stabilité CFL, et les éléments finis mixtes étendus pour le terme diffusif. Pour les réactions chimiques, on utilise la méthode de Runge-Kutta si le type des réactions est cinétique, et la méthode de minimisation point-intérieur, si le type des réactions est équilibre.