La technologie la plus utilisée pour séparer les principaux composants de l'air est aujourd'hui ladistillation en colonne à garnissage structuré. Ce procédé se caractérise par un écoulement gazliquide contre-courant au sein d'une structure constituée de plaques corruguées placées parallèlement et agencées en packs. La description d'un tel procédé est rendue difficile par les très grandes dimensions du système et par la complexité des phénomènes à la petite échelle. La méthode de la prise de moyenne volumique, utilisée pour les problèmes de changement d'échelleen milieu poreux, est utilisée pour décrire le système à une échelle qui permet une résolution. Le travail est organisé en trois étapes. Dans un premier temps, pour les débits modérés, une méthode est proposée pour évaluer la perte de charge au sein de la structure en prenant en compte des rugosités de structure ou créées par des instabilités du film liquide. A ce stade, l'effet de la surface rugueuse est caractérisé par une condition limite effective. Le problème aux limites effectif pour la phase gaz est ensuite moyenné en volume pour obtenir un système d'équations à grande échelle. Le bilan de quantité de mouvement à grande échelle est une loi de Darcy généralisée aux écoulements inertiels, dans laquelle les paramètres effectifs contiennent les effets des instabilités de surfaces de la petite échelle. La seconde étape est dédiée à l'interaction entre les deux phases à plus hauts débits. On montre que des modèles qui incluent explicitement des termes croisés à grande échelle permettent de décrire l'écoulement au sein du garnissage à grands nombres de Reynolds. Plus généralement, ces modèles, peu utilisés dans la littérature sur les milieux poreux, s'avèrent adaptés pour les écoulements dans les milieux très perméables, pour lesquels des variations importantes de la perte de charge et des saturations sont observées. Enfin, on s'intéresse à la description de la distribution de la phase liquide au sein de la structurede garnissage. Une approche multiphasique, où la phase liquide est séparée en plusieurs pseudophases, est adaptée pour modéliser l'anisotropie de l'écoulement. Deux méthodes impliquant une approche à deux pseudo-phases et une approche à quatre pseudo-phases pour la phase liquide sont comparées. Cette dernière méthode est notamment utile pour décrire des régimes d'écoulement très différents, et permet de capturer à grande échelle les chemins préférentielssuivis par le film liquide au sein du garnissage. / Distillation in columns equipped with structured packings is today the most used technology for separating air in its primary components. This process is characterized by a counter-current gasliquid flow in a structure made of parallel corrugated sheets arranged in packs. The description of such system is constrained by the large dimensions of the columns and by the complexity of the local-scale phenomena. This leads to consider a strategy of upscaling, based on the volume averaging method, to describe the system at a scale at which a resolution is possible. The work is organized in three steps. As a first step, considering moderate flow rates, a methodology ofupscaling is developed to predict the pressure drop in the flow of the gas phase taking into account small scale roughnesses due to the structure itself or perturbations of the liquid film. At this stage, the effect of this rough surface is characterized by an effective boundary condition. The boundary value problem for the flow of the gas phase is volume averaged in order to derive a system of equations at large scale. The resulting momentum balance is a generalized Darcy's law for inertial flows, involving effective parameters accounting for the roughness at the microscale. The second step of this work focuses on the interaction between the two phases at higher flow rates. It is shown that models involving non-standard macroscopic cross-terms are more prone to describe the flow in packings at high Reynolds numbers than the models usually used in porousmedia sciences. More generally, these models are shown to characterize accurately processes in highly permeable media, where drastic changes of pressure drop and retention are observed. We finally study the distribution of the liquid phase in the structured packing. It is shown that a specific approach involving a multiphase model with liquid decomposition is required to capture the anisotropy generated in the flow of the liquid phase. Two methods involving two pseudo-phases and four pseudo-phases for the liquid phase are compared. This last method captures a number of very different distribution regimes in the column and offers additional flexibility to describepreferential paths of the liquid.

Pour la future génération d’aimants supraconducteurs à hauts champs magnétiques refroidis avec de l’hélium superfluide comme les prochains aimants développés pour l’amélioration des performances du Large Hadron Collider (LHC), les ingénieurs prennent en considération l’utilisation de supraconducteurs en niobium-étain. Dans leur environnement, ces aimants supraconducteurs vont être soumis à des pertes thermiques bien plus importantes que celles des aimants actuels des accélérateurs de particules. Alors que la résistance thermique due aux isolations électriques des câbles supraconducteurs constitue la principale résistance thermique au refroidissement, un nouveau type d’isolation, à base de céramique poreuse, est considéré. Afin de comprendre la capacité de refroidissement et la stabilité thermique de ces aimants, il devient alors nécessaire d’étudier les transferts thermiques en hélium superfluide dans les milieux poreux. Dans ce but, plusieurs dispositifs expérimentaux ont été mis au point afin de tester des milieux poreux modèles dont on connaît les principales caractéristiques (taille de pore, porosité, épaisseur). Par ailleurs, une étude théorique utilisant une méthode de changement d’échelle de type prise de moyenne volumique a permis d’obtenir une modélisation macroscopique des phénomènes étudiés et a montré la validité de la loi de Darcy en hélium superfluide sous certaines hypothèses. Ceci a été vérifié par des simulations numériques directes. Les résultats expérimentaux ont permis de mettre en exergue différents résultats ; tout d’abord, le régime d’écoulement est fonction de la taille des pores. Pour les pores les plus petits (10-7 et 10-6 m), le régime d’écoulement est le régime de Landau et ceci a permis de calculer la perméabilité en He II. Pour des tailles de pores de 10-5 m, le régime d’écoulement est le régime de Gorter-Mellink et les transferts de chaleur ont pu être modélisés à partir de la loi de Gorter-Mellink en introduisant la notion de tortuosité / For the next generation of high field magnets cooled with superfluid helium like magnets for the next particles accelerator or the upgrade of the Large Hadron Collider (LHC), engineers are taking into consideration the use of Nb3Sn superconductors. In their environment, these superconducting magnets will undergo much higher heat losses than in current particles accelerator. Since the thermal resistance due to the electrical insulations of the superconducting cables constitutes the main thermal resistance for cooling, new type of insulation, based on ceramic materials, are considered. In order to understand the cooling capacity and the thermal stability of these magnets, it is necessary to investigate the heat transfer through porous media in superfluid helium. In this purpose, several experiments were set up in order to study porous medium models for which we know the main characteristics (average pore size, porosity, thickness). Also, a theoretical study using an up-scaling method of type volume averaging allowed a macroscopic modeling of the studied phenomena and showed that Darcy’s law is valid in superfluid helium under given hypothesis. This has been verified by direct numerical simulations. Experimental results allowed to highlight several results; first, the flow regime depends on the pore size diameter. For the smallest pore (10-7 and 10-6 m), the flow regime is the Landau regime and this permitted to calculate the permeability in He II. For an average pore size of 10-5 m, the flow regime is the Gorter-Mellink regime and heat transfer were modeled using the Gorter-Mellink’s law by introducing the notion of tortuosity

Ce travail de thèse, réalisé en collaboration avec l’IMFT et TOTAL, traite de la modélisation de la combustion in-situ appliquée à une huile lourde Vénézuéliène. Il a été initié suite à une observation simple : même si le procédé est étudié depuis plusieurs décénies, on ne peut pas encore le modéliser correctement. Des résultats expérimentaux, issus d’expérience à l’échelle du laboratoire (tubes à combustions), ne sont pas reproductibles avec des outils numériques commerciaux de types simulateurs réservoirs thermiques. Par conséquent, et face à ce constat, nous avons été contraint d’explorer plusieurs pistes pour améliorer la modélisation du procédé : – La chimie et les méthodes de détermination de mécanismes réactionnels. – La description thermodynamique d’une huile lourde et le calcul d’équilibre triphasique. – Le transport de masse et de chaleur dans un milieu poreux, en situation multiphasique, réactive et miscible. – La conception d’un modèle mathématique et numérique d’un modèle complet. Nous pensons que le problème pluridisciplinaire et fortement complexe peut trouver une réponse si l’ensemble des mécanismes et leurs liens sont traités de façon adéquate. Une campagne expérimentale (expériences de cellules cinétiques), portant sur l’étude des effets de l’eau sur les réactions chimiques de l’huile, a permis de mettre en évidence des effets inattendus et nouveaux. Ces données, complétées par des expériences de types tubes à combustion, fournissent une importante base de données expérimentale. Pour modéliser les expériences de cellules cinétiques, nous avons tout d’abord développer un nouvel outil de simulation directe, reposant sur une description compositionnelle de l’huile où les comportements de toutes les phases sont prédits par les équations d’états. Le calcul d’équilibre est fait grâce à un flash diphasique. Afin de déterminer un mécanisme réactionnel paramétré, nous avons couplé ce dernier outil à un algorithme génétique. Finalement, dans le but de simuler les expériences de tubes à combustion, un nouveau simulateur compositionel, triphasique, thermique et réactif a été développé. Il est spécialement adapté à la simulation de ce genre d’expérience. Le calcul d’équilibre de phase est réalisé grâce à un nouvel outil développé pour l’occasion. Ce dernier repose sur l’hypothèse free water et repose sur une formulation originale et novatrice. / The study of this PhD, realized jointly with IMFT and TOTAL, deals with modeling of in-situ combustion applied to a Venezuelan heavy oil. It has begun with a relatively simple observation: even if the process has been extensively studied since some decades, we cannot correctly model it. Experiment data provided by lab scale experiments (combustion tubes) mismatches numerical results obtained from commercial thermal simulator, especially for wet experiments. The need to better understand the process related to this issue forced us to explore multiple tracks for various scientific fields. Thus, one can cite: • The chemistry and methods of reduction of reactive mechanisms. • The thermodynamic description of the heavy oil and the calculations of three-phase equilibrium. • Heat and mass transport in multiphase, reactive and miscible porous medium. • Mathematical and numerical design of a full model. The problem exceedingly complex can find a complete and consistent answer if one takes into account the whole mechanisms and links between them. We have followed this way in order to determine a robust reactive scheme using both theoretical numerical and experimental developments. A whole set of kinetic cell manipulations was conducted to better understand and discriminate the effects of water on chemistry on a certain type of heavy oils. New interactions and effects on steam on heavy oil combustion have been discovered and studied. These manipulations, supplemented by a set of some combustion tubes provide a large set of experimental data. This will compose our base case that we will try to match later using some new tools devised during this study. To model kinetic experiments, we firstly developed a new simulation tool based on a compositional description and a full equation of state formulation. Equilibrium calculation is made by a two-phase flash. To determine consistent kinetic parameters, we used a genetic algorithm coupled with the new tool. Finally, in order to validate the kinetic model and simulate combustion tube experiment, a new threephase compositional simulator has been developed. It is especially fitted to take into account characteristic of the experimental device. Three-phase equilibrium calculation is computed by a new free-water

Des phénomènes de colmatage de la structure granulaire des réservoirs pétroliers sont observés lors de la réinjection d'eau de production dans le réservoir (endommagement hydrique). Le dépôt des particules solides en suspension dans l'eau de production peut engendrer une réduction de la perméabilité du réservoir notamment autour des puits d'injection. Il en résulte une diminution significative de l'injectivité des puits. L'objectif de cette thèse est d'étudier le processus de transport et de dépôt des particules solides dans un réservoir de sable et son influence sur la perméabilité du milieu. Nous étudions également l'influence de certains paramètres tels que : la concentration en particules, le débit d'injection, la taille de particules injectées, la taille des pores du milieu et la rugosité de la surface des grains. La finalité de ce travail est de proposer un modèle prédictif de transport-dépôt-endommagement afin de pouvoir estimer la pérennité du système. Dans ce but, une campagne d'essais expérimentaux est menée sur le processus de transport et de dépôt des particules siliceuses de taille de quelques microns dans une éprouvette de sable dense. Ces essais ont été réalisés sur trois dispositif différents ; deux dispositifs d'injection unidimensionnelle (colonne d'injection et cellule d'injection) et un dispositif d'injection radiale en chambre d'étalonnage en vue de simuler les conditions d'injection en puits. Dans ces essais, les conditions hydrauliques sont contrôlées et la concentration en particules dans la suspension est imposée. Au cours de l'injection, le suivi de l'évolution de la pression interstitielle est réalisé en utilisant des capteurs de pression installés à différents endroits dans le milieu. La concentration des particules dans le fluide sortant est mesurée en utilisant un turbidimètre installé à la sortie du milieu. L'étude expérimentale réalisée sous différentes conditions a permis de mettre en évidence une forte hétérogénéité du dépôt des particules dans le milieu poreux. Cette hétérogénéité se traduit par la formation d'un cake interne et/ou externe au voisinage du point d'injection. L'analyse de ce cake montre un milieu à double porosité avec deux familles de pores correspondant à la porosité interne des agglomérats formés par les particules déposées et la porosité du sable réduite progressivement par le dépôt de ces derniers. Une étude paramétrique a permis de mettre en évidence une distinction entre l'effet des faibles et des fortes concentrations, ainsi que l'importance de l'effet du débit d'injection, de la taille des pores, de la taille des particules de la suspension ainsi que de la rugosité de la surface des grains. Ce travail expérimental s'accompagne d'une modélisation des phénomènes de transport et de dépôt des particules dans un milieu poreux. Dans un premier temps, nous avons utilisé un modèle de convection dans lequel la loi décrivant le dépôt des particules prend en compte la densité des particules déposées ainsi que le relargage des particules. Les paramètres de la loi constitutive de dépôt, ainsi qu'une loi empirique porosité-perméabilité ont été calibrés en simulant les résultats expérimentaux d'essais d'injection en cellule sous différentes conditions. Un nouveau modèle semi analytique basé sur le modèle d'Eylander (1988) a été développé dans lequel nous prenons en compte une distribution spatiale des particules déposées et une double porosité du milieu (cake et milieu granulaire). Cette distribution est inspirée du modèle convection étudié auparavant. La simulation d'essais d'injection en cellule a permis de montrer les capacités du modèle pour la prédiction des variations de la pression en cours de l'essai. Ce modèle permet une prédiction de l'évolution de la perméabilité pendant l'injection de façon plus simple et directement applicable dans l'ingénierie de réservoirs / Hydraulic damage and plugging of the granular structure of oil reservoirs are commonly observed during the re-injection of produced water. These phenomena can lead to a significant decrease of the reservoir permeability around the injectors and affect the injectivity of these wells. The objective of this thesis is to study the transport and deposition process of solid particles in a sand reservoir and its influence on the permeability of the medium. We also study the influence of certain parameters such as particle concentration, injection flow rate, the size of injected particles, the pore size of the medium and the surface roughness of the grains. The aim of this work is to propose a predictive model of transport-deposition-damage in order to estimate the durability of the system. For this purpose, an experimental program was conducted on the transport and deposition process of the siliceous particles of size of a few microns in dense sand specimens. These tests were carried out on three different devices; two one-dimensional injection devices (injection column and injection cell) and a radial injection device in “chambre d'étalonnage” in order to simulate the injection wells conditions. In these tests, the flow conditions are controlled and the concentration of particles in the suspension is imposed. During injection, the monitoring of the evolution of pore pressure is carried out using pressure sensors located in different places in the medium. The particle concentration in the outgoing fluid is measured using a turbidimeter installed at the outlet. The experimental study carried out under different conditions allowed to highlight a highly heterogeneous particle deposition in the porous medium. This heterogeneity results in the formation of an internal and/or external cake in the vicinity of the injection point. The analysis of this cake shows a double porosity medium with two pore families corresponding to the internal porosity of the agglomerates formed by the deposited particles and porosity of the sand gradually reduced by the particles deposition. A parametric study permitted to show a distinction between the effect of low and high particles concentrations, as well as the importance of the effect of injection flow rate, the pore size, the particle size of the suspension and the roughness of the surface of the grains. This experimental work is associated with modelling of transport and deposition of particle in a porous medium. We used first the convection model in which a constitutive law describing the particle deposition takes into account the effects of the deposited particles concentration and the release of particles. The parameters of the particle deposition constitutive law as well as an empirical porosity-permeability law have been calibrated by simulation of the results of injection tests performed in the injection cell under different conditions. A semi-analytical model based on the Eylander (1988) model was developed in which we take into account a spatial distribution of the deposited particles and a double porosity of the medium (cake and granular medium). This distribution is inspired from the results if convection model studied before. The simulation of the results of injection experiments permitted to show the performances of this model for prediction of the pore pressure evolution during the fluid injection. This model allows prediction of the permeability evolution in a simpler and directly applicable manner in the reservoir engineering

Ce travail a pour objectif d'identifier et de modéliser les mécanismes physico-chimiques régissant les écoulements d'eau et les transferts de solutés dans un milieu poreux consolidé fissuré.<br />Un dispositif expérimental original a été mis en place, où un milieu fissuré modèle est reproduit en empilant des cubes poreux consolidés de 5 cm de côté : le « cube ». En parallèle, des études sont menées sur des colonnes du milieu poreux consolidé homogène. La même technique de traçage anionique est utilisée dans les deux cas. Selon une approche de dynamique des systèmes, on sollicite le dispositif avec des créneaux de concentration en traceur, pour obtenir des courbes de percée. Après avoir identifié bilan de masse et temps de séjour, nous avons calé les modèles CD et MIM sur les données expérimentales. <br />Le modèle MIM permet de reproduire les courbes expérimentales sur le milieu poreux consolidé homogène de façon plus satisfaisante que le modèle CD. On obtient une fraction d'eau mobile en cohérence avec la géométrie du milieu poreux. L'étude de l'influence de la vitesse de l'écoulement met en évidence un régime de dispersion d'interférence. L'influence de la longueur d'observation n'a pas pu être mise en évidence dans ce cas.<br />Nous avons par contre mis en évidence un effet de l'échelle d'observation sur le milieu poreux et fissuré en comparant les résultats obtenus sur un petit « cube » et un grand « cube ». Le modèle CD n'est pas satisfaisant dans ce cas. Si le modèle MIM permet de reproduire les courbes de percée expérimentales, il n'a pas été possible de caler des paramètres uniques pour toutes les expériences.

La thèse vise à quantifier les effets de la variabilité spatiale d'un milieux poreux sur l'évolution d'un système géochimique. Les réactions de dissolution ou de précipitation des minéraux modifient la structure microscopique du milieu, et par suite les caractéristiques hydrodynamiques du système, notamment la perméabilité. La variabilité spatiale du milieu initial peut-elle alors être à l'origine de la formation de digitations ou de chenaux ? La première partie traite du changement d'échelle, pour le passage d'une simulation géostatistique sur grille fine au calcul de transport sur un maillage plus grossier. Dans le cas du code Hytec, qui utilise un schéma aux volumes finis fondé sur une discrétisation en polygones de Voronoï, plusieurs méthodes de calcul de la perméabilité équivalente sont comparées, suivant différents critères. La deuxième partie aborde les calculs de transport réactif sur une famille de simulations géostatistiques du milieu. L'influence de la variabilité spatiale initiale sur l'évolution des systèmes est quantifiée à l'aide d'observables adaptées. Deux réactions distinctes sont étudiées, un cas de dissolution, de façon approfondie, et, plus rapidement, un cas de précipitation tendant au rééquilibrage du système.

La mouillabilité mixte ou partielle dans les milieux poreux influe considérablement sur le comportement des écoulements diphasiques. Dans le monde, la plupart des réserves pétrolières se situent dans des réservoirs possédant une mouillabilité hétérogène. Des études détaillées de l'influence de la mouillabilité partielle sur le comportement des écoulements dans les milieux poreux de type modèle micro fluidique sont par conséquent d'un intérêt majeur, tant pour la recherche fondamentale que pour les études appliquées. Nous avons développé une nouvelle technique permettant pour la première fois la modification et la structuration de la mouillabilité de surface dans les systèmes micro fluidique de type complexe, système d'une grande homogénéité et de haute qualité fabriqué en poly(dimethylsiloxane) (PDMS). Ces progrès obtenus dans le traitement des surfaces, nous ont permis d'étudier l'influence de la mouillabilité partielle sur le comportement des écoulements diphasiques et d'observer de nouveaux régimes d'écoulements. L'observation des régimes d'écoulement particuliers en fonction de la mouillabilité dans des carottes de roche nous a permis de développer une nouvelle technique de détection basée sur une visualisation optique de celle-ci. Ce concept de mesure d'écoulements par visualisation optique à l'échelle du pore fut étendu et un nouvel outil d'imagerie microscopique a été développé.

L'objectif de cette thèse est l'étude du problème de Cauchy pour les solutions faibles de trois problèmes (systèmes paraboliques dégénérés et fortement couplés) modélisant des écoulements diphasiques et compressibles en milieu poreux. La motivation de ce travail est un "benchmark" du GNR MoMaS pour l'étude de l'impact de l'écoulement du gaz d\^{u} à la corrosion des matériaux ferreux dans un site de stockage de déchets radioactifs. Cette thèse est divisée en trois chapitres indépendants. Premièrement, on s'intéresse à l'analyse mathématique d'un problème modélisant l'écoulement de deux phases immiscibles et en considérant qu'une phase est compressible et l'autre est incompressible (eau/gaz). Deuxièmement, on traite le cas général du déplacement de deux fluides compressibles et immiscibles dans un milieu poreux. Enfin, le dernier chapitre est consacré à la construction et à la convergence de la méthode des volumes finis pour le système eau-gaz sous l'hypothèse que la densité du gaz est une fonction de la pression globale.

Cette thèse porte sur l'étude du comportement hydromécanique des déchets biodégradables à deux échelles : au laboratoire et dans les sites de stockage. Après une introduction sur les Installations de Stockage et le traitement final des déchets(Chapt.1), les paramètres physiques, mécaniques et hydrauliques sont définis par analogie avec les milieux granulaires(Chapt.2).Les recherches présentées se rapportent d'une part aux propriétés de transfert et d'autre part aux propriétés mécaniques d'un déchet biodégradable. Des essais originaux de compression-perméabilité au gaz ont été réalisés sur des déchets jeunes et des déchets biodégradés au laboratoire dans le but de mesurer l'évolution de la conductivité au gaz en fonction du tassement et du degré de saturation en lixiviat (Chapt.3). Ces essais permettent de valider le modèle de double porosité proposé par Stoltz (Thèse 2009) (Chapt.4). Dans un deuxième temps, on s'intéresse aux propriétés mécaniques : une nouvelle version du logiciel ISPM est proposée, il s'agit d'un outil de modélisation destiné à la prévision du tassement des déchets -Olivier (Thèse 2003). Ce nouveau modèle appliqué à l'étude du tassement de deux alvéoles en service , a permis par méthode inverse de déterminer le coefficient de compressibilité secondaire de ces déchets (Chapt.5). Enfin, la résistance au cisaillement de déchets biodégradables issus de forages est caractérisée en fonction de différents paramètres d'état (Chapt.6) et ces valeurs sont utilisées pour un calcul de stabilité sur site.

Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark « Ecoulement diphasique » proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé / The two-phase flow in porous media is a complex phenomenon and which relate to many industrial problems. EDF works on the feasibility and the safety of a storage in deep geologic layer of nuclear waste. In this domain the simulation of the two-phase flow in porous media is particularly important in at least three domains : first of all during the phase of ventilation of the galleries of the storage which could desaturate the rock and so modify its properties, but also during the phase of resaturation of the materials and finally during the arrival of the water on the metal parts contained in the storage which will then involve phenomena of corrosion and a hydrogen release. In this context, EDF wishes to obtain robust numerical methods without restrictive condition on the mesh. This work is dedicated at first to the development of the finite volume scheme SUSHI (Scheme Using Stabilization and Hybrid Interfaces) in the code of mechanics of EDF, Code_Aster in order to simulate the two-phase flow in porous media. This scheme was developed in 2D and in 3D. At the same time a new formulation which allows to simulate in a uniform way the flows in saturated and unsaturated porous media for miscible and immiscible problems is proposed. Various studies simulating difficulties related to the problems of the storage of nuclear waste in deep geological layers were study. We can quote the study of a bi-material which advances the capillary rebalancing of a material by an other one possessing properties and initial very heterogeneous conditions in saturation. We will also quote the study of the injection of hydrogen in an porous media initially saturated in pure water which is proposed by the benchmark "two-phase Flow " proposed by the GNR MOMAS. This study had for objective to bring to light the good treatment of the appearance of a phase in a saturated porous media and thus the relevance of our new formulation to study with a way unified a problem of saturated flow and a problem of unsaturated flow