Le probleme d'ecoulement biphasique est modelise par un systeme d'equations de conservation de masse en saturation-pression-vitesse et le probleme de transport reactif par un systeme d'equations de conservation de concentration de type advection- diffusion-reaction. On applique la methode des elements finis mixtes etendus pour approcher simultanement les pressions, les saturations et les champs de vitesse du modele biphasique. Le schema utilise est totalement implicite et non lineaire. Il est resolu par la methode de Newton inexacte. On utilise les espaces de Raviart-Thomas et des formules de quadratures pour obtenir un schema equivalent aux differences finies centrees. Pour le transport reactif, on applique la methode de splitting d'operateurs avec une discretisation temporelle semi-implicite, un schema de Godunov pour approcher le terme convectif, sous la condition de stabilite CFL, et les elements finis mixtes etendus pour le terme diffusif. Pour les reactions chimiques, on utilise la methode de Runge-Kutta si le type des reactions est cinetique, et la methode de minimisation point-interieur, si le type des reactions est equilibre."