La modélisation des hydrosystèmes souterrains est devenue un outil crucial dans la gestion des ressources d'eaux souterraines ainsi que pour la surveillance des sites contaminés. L'objectif de ce travail de thèse est de modéliser l'écoulement d'une phase non-aqueuse liquide dense (DNAPL) dans les sols, et de développer un code pour la simulation des écoulements triphasiques de fluides incompressibles en milieux poreux. Le modèle mathématique pour les écoulements multiphasiques de fluides en milieux poreux est généralement constitué d'une équation de pression et de deux équations de saturation. Notre approche est fondée sur une formulation en pression globale : elle permet un découplage partiel des équations de pression et de saturations, et elle est plus efficace en termes de résolution numérique. Le nouveau modèle est discrétisé selon un schéma IMPES. Dans ce travail de thèse, la méthode des Eléments Finis Discontinus de Galerkine est combinée à un schéma de Godunov généralisé pour la résolution de la partie convective de l'équation de saturation. Un écoulement diphasique non miscible avec effets gravitaires importants et sans capillarité a été simulé. Une analyse fonctionnelle démontre que le profil de saturation de la phase non-aqueuse entrante dépend essentiellement du rapport de la vitesse totale sur la différence de densité entre la phase non aqueuse et l'eau. Le code développé a permis de simuler un drainage gravitaire à grande échelle. Enfin, un modèle numérique non-linéaire d'écoulement triphasique de fluides utilisant la Méthode des Lignes est introduit. Les causes d'oscillations instables du système en zone elliptique sont examinées. Des méthodes de relaxation et la construction d'un modèle de perméabilités relatives vérifiant la condition de Différentielle Totale sont envisagées. / Numerical simulation has become a crucial tool in addressing water-resource management and other environmental problems such as polluted sites monitoring. The aim of this work is to model the flow of a dense non aqueous phase liquid (DNAPL) in the subsurface by developing a numerical code to simulate three-phase (DNAPL, water, and gas), incompressible flow in porous media. The mathematical model for multiphase flow in porous media is generally composed of a system of one pressure and two saturation equations. Our approach is based on a global pressure model : it leads to a partial decoupling of the pressure and the saturation equation and is more efficient from the computational point of view. The new model is discretized by a Mixed Hybrid Finite Elements Method (MHFEM), Discontinuous Galerkin Finite Elements (DGFEM), IMPES resolution method. In this work, the DGFEM scheme is combined with a generalised Godunov scheme to solve the convective part of the saturation equation. An immiscible two-phase flow with predominant gravity effects and whithout capillary effects has been modelled. It has been shown that the saturation profile of a displacing non-aqueous phase liquid (NAPL) in an initially water-saturated porous medium depends strongly on the ratio between the total specific discharge and the density difference between the NAPL and water. The 2D-code enabled a simulation of a large-scale gravity drainage. Finally, a non-linear three-phase 1D flow formulation using Method of Lines (MOL) has been introduced. Unstable oscillatory behavior of the system when the initial state are in the elliptic region of the ternary diagram is examined. Non-equilibrium formulation and construction of relative permeability model satisfying the Total-Differential are foreseen.