Ce travail cherche à élaborer un cadre théorique unifié et cohérent, fondé sur les principes rigoureux de la thermodynamique et permettant de modéliser une classe étendue de phénomènes physiques intervenant au sein des milieux poreux non saturés, ainsi que leur interaction avec l'état mécanique des structures. Cette classe est celle des phénomènes réactifs, parmi lesquels le changement de phase de l'eau interstitielle (dessiccation), modélisé dans l'ensemble ses manifestations (fluage, mais surtout retrait) sera abordé comme exemple particulier d'application. Le premier chapitre rappelle les fondements de la description adoptée du milieu poreux, ainsi que le cadre thermodynamique qui sous-tend l'ensemble de la modélisation. Le chapitre II s'attache à la formulation -essentiellement nouvelle- et à l'identification d'une loi de comportement élastique non linéaire du milieu. Différents éléments de réflexion liés aux comportements microscopiques des constituants et à l'expérience, orientent le modèle vers certaines hypothèses fortes qui permettent la détermination complète et explicite de la loi de comportement macroscopique. Le chapitre III expose les différentes expériences à disposition pour l'identification de la loi de comportement, et permet d'intégrer ces expériences, en tant que caractéristiques comportementales, dans le modèle théorique. Chapitre IV et chapitre V sont des exemples d'application : le chapitre IV étudie le problème du retrait et du fluage de dessiccation en comportement couplé élastique linéaire ; le chapitre V se limite au cas du retrait de dessiccation, par application du modèle poroélastique non linéaire complet. Les résultats obtenus présentent une bonne adéquation avec l'expérience.

Cette thèse porte sur la modélisation et la discrétisation d’écoulements Darcy dans les milieux poreux fracturés. Nous suivons l’approche des modèles, dits à dimensions hybrides, qui représentent les réseaux de fractures comme des surfaces de codimension 1 immergées dans la matrice. Les modèles considérés prennent en compte les interactions entre matrice et fractures et permettent de traiter des fractures agissant comme conduites ou comme barrières, ce que nécessite de prendre en compte les sauts de pression aux interfaces matrice-fracture. Dans le cas des écoulements diphasiques, nous proposons des modèles, qui prennent en compte les sauts de saturations aux interfaces matrice-fracture, dû à la capillarité. L’analyse numérique est menée dans le cadre général de la méthode de discrétisations gradients, qui est étendue aux modèles considérés. Deux familles de schémas numériques, le schéma Vertex Approximate Gradient et le schéma Volumes Finis Hybrides sont adaptées aux modèles à dimensions hybrides. On prouve via des résultats de densité que ce sont des schémas gradients, pour lesquels la convergence est établie. En diphasique, l’existence d’une solution est obtenue en passant. Plusieurs cas tests sont présentés. En monophasique, on observe la convergence sur des différents types de mailles pour une famille de solutions dans un milieux fracturé hétérogène et anisotrope. En diphasique, nous présentons une série de cas tests afin de comparer les modèles à dimensions hybrides au modèle de référence, dans lequel les fractures ont la même dimension que la matrice. A part quantifier le gain en performance de calcul, ces tests montrent la qualité des différents modèles réduits. / This thesis investigates the modelling of Darcy flow through fractured porous media and its discretization on general polyhedral meshes. We follow the approach of hybrid dimensional models, invoking a complex network of planar fractures. The models account for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrix-fracture interfaces. In the case of two phase flow, we present two models, which permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The numerical analysis is performed in the general framework of the Gradient Discretisation Method, which is extended to the models under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to fit in the gradient scheme framework, which yields, in particular, convergence. For single phase flow, we obtain convergence of order 1 via density results. For two phase flow, the existence of a solution is obtained as a byproduct of the convergence analysis. Several test cases are presented. For single phase flow, we study the convergence on different types of meshes for a family of solutions. For two phase flow, we compare the hybrid-dimensional models to the reference equidimensional model, in which fractures have the same dimension as the matrix. This does not only provide quantitative evidence about computational gain, but also leads to deep insight about the quality of the proposed reduced models.

Dans un contexte de réduction internationale des émissions de gaz à effet de serre, les techniques de Captage Transport et Stockage de CO2 (CTSC) apparaissent comme une solution à moyen terme particulièrement efficace. En effet, les capacités de stockage géologique pourraient s'élever jusqu'à plusieurs millions de tonnes de CO2 injectées par an, soit une réduction substantielle des émissions atmosphériques de ce gaz. Une des cibles privilégiées pour la mise en place de cette solution sont les aquifères salins profonds. Ces aquifères sont des formations géologiques contenant une saumure dont la salinité est souvent supérieure à celle de la mer la rendant impropre à la consommation. Cependant, cette technique fait face à de nombreux défis technologiques ; en particulier la précipitation des sels, dissous dans l'eau présente initialement dans l'aquifère cible, suite à son évaporation par le CO2 injecté. Les conséquences de cette précipitation sont multiples, mais la plus importante est une modification de l'injectivité, c'est-à-dire des capacités d'injection. La connaissance de l'influence de la précipitation sur l'injectivité est particulièrement importante tant au niveau de l'efficacité du stockage et de l'injection qu'au niveau de la sécurité et de la durabilité du stockage. Le but de ces travaux de thèse est de comparer l'importance relative des phénomènes négatif (colmatage) et positif (fracturation) consécutifs à l'injection de CO2 et à la précipitation des sels. Au vu des nombreux résultats de simulations et de modélisation dans la littérature décrivant le colmatage de la porosité, il a été décidé de porter l'accent sur les effets mécaniques de la cristallisation des sels et la possible déformation de la roche mère. Une modélisation macroscopique et microscopique, tenant compte de deux modes possibles d'évaporation induits par la distribution spatiale de l'eau résiduelle a donc été développée afin de prédire le comportement mécanique d'un matériau poreux soumis à un assèchement par injection de CO2. Les résultats montrent que la pression de cristallisation consécutive à la croissance d'un cristal en milieu confiné peut atteindre des valeurs susceptibles localement de dépasser la résistance mécanique du matériau, soulignant ainsi l'importance de ces phénomènes dans le comportement mécanique global de l'aquifère. Sur le plan expérimental, les travaux ont porté sur l'utilisation d'un nouveau prototype de percolation réactive afin de reproduire le comportement d'une carotte de roche soumise à l'injection et ainsi obtenir l'évolution des perméabilités dans des conditions similaires à celle d'un aquifère.