Nous développons des algorithmes parallèles pour la résolution de problèmes non-linéaires de grande taille. Les cadres d'application sont la simulation de matériaux hyperélastiques incompressibles en grandes déformations et l'étude des milieux poreux, dont les modélisations choisies font apparaître des inconnues en déplacement et en pression.<br /><br />Nous retenons une stratégie éléments-finis associée à un solveur Newton-Raphson et une décomposition de domaine sans recouvrement combinée à un solveur de Krylov. <br /><br />Nous proposons des améliorations pour adapter ces approches à nos problèmes, puis pour des cas plus exigeants nous définissons une nouvelle approche de décomposition de domaine, appelée approche hybride, permettant de mieux respecter la physique des phénomènes et unifiant les approches classiques. Nous proposons également des stratégies d'accélération du processus non-linéaire. Enfin un cadre orienté objet est exposé pour la mise en oeuvre de l'ensemble des méthodes proposées.

Cette thèse s'intéresse à la conception de méthodes de discrétisation non-conforme pour un modèle de poromécanique. Le but de ce travail est de simplifier les couplages liant la géomécanique d'un milieu poreux à l'écoulement polyphasique compositionnel ayant cours en son sein tels qu'ils sont réalisés actuellement dans l'industrie pétrolière, en discrétisant sur un même maillage, typiquement non-conforme car à l'image de la lithologie, la mécanique et l'écoulement. La nouveauté consiste donc à traiter la mécanique par une méthode d'approximation non-conforme sur maillages généraux. Dans cette thèse, nous nous concentrons sur un modèle d'élasticité linéaire. Les difficultés inhérentes à son approximation non-conforme sont son manque de coercivité (se traduisant par la nécessité de satisfaire une inégalité de Korn sur un espace discret discontinu), ainsi que le phénomène de verrouillage numérique lorsque le matériau tend à devenir incompressible. Dans une première partie, nous construisons un espace d'approximation sur maillages généraux, s'apparentant à une extension de l'espace de Crouzeix-Raviart. Nous explicitons ses propriétés d'approximation et de conformité, et montrons que ce dernier est adapté à une discrétisation primale coercive et robuste au locking du modèle d'élasticité sur maillages généraux. La méthode proposée est moins coûteuse que son équivalent éléments finis (en termes de propriétés) P2. Nous nous intéressons dans une deuxième partie à l'approximation non-conforme d'un modèle couplé de poroélasticité. Nous étudions la convergence d'une famille de schémas numériques dont la discrétisation en espace utilise le formalisme des schémas Gradient, auquel appartient la méthode développée pour la mécanique. Nous prouvons la convergence de telles approximations vers la solution de régularité minimale du problème continu, indépendamment des paramètres physiques du système / This manuscript focuses on the conception of nonconforming discretization methods for a poromechanical model. The aim of this work is to ease the coupling between the geomechanics and the multiphase compositional Darcy flow in porous media by discretizing mechanics and flow on the same mesh, typically nonconforming as it represents the lithology. Hence, the novelty hinges on a nonconforming treatment of mechanics on general meshes. In this work, we focus on a linear elasticity model. The nonconforming approximation of such a model is not straightforward owing to its lack of coercivity (meaning that a discrete Korn's inequality must hold on a discontinuous discrete space) and to the numerical locking phenomenon occurring as the material becomes incompressible. In a first part, we design an approximation space on general meshes, which can be viewed as an extension of the so-called Crouzeix-Raviart space. We study its approximation and conformity properties, and prove that this latter is well-adapted to the design of a primal, coercive, and locking-free discretization of the elasticity model on general meshes. The proposed method is less costly than its finite element equivalent (in terms of properties) P2. In a second part, we tackle the nonconforming approximation of a coupled poroelasticity model. We study the convergence of a family of numerical schemes whose space discretization relies on the Gradient schemes framework, to which belongs the method developed for mechanics. We prove the convergence of such approximations toward the minimal regularity solution of the continuous problem, and independently of the choice of physical parameters

La propagation des ondes sismiques dans les milieux poreux multiphasiques présente des enjeux nombreux, tant sur le plan environnemental (risques naturels, géotechnique, pollutions de nappes...) que pour les réservoirs (aquifères, hydrocarbures, stockages de CO2...). L'utilisation des ondes sismiques pour étudier ces milieux se justifie par le fait qu'en se propageant, les ondes sont déformées par le milieu qu'elles traversent et contiennent ainsi des informations aux capteurs sur les phases fluides et solides et sur le squelette poreux. Ce travail de thèse s'intéresse aux caractéristiques des ondes sismiques dans les milieux multiphasiques (plusieurs phases fluides et solides), depuis la description physique jusqu'à la caractérisation des paramètres constitutifs par inversion, en passant par la modélisation numérique 2D de la propagation. La première partie du travail a consisté à décrire la physique des milieux multiphasiques (phase par phase et leurs intéractions dynamiques) en utilisant des méthodes d'homogénéisation pour se ramener à un milieu équivalent défini par sept paramètres. Ainsi, dans des milieux simple porosité saturés et dans des milieux plus complexes (double porosité, partiellement saturés ou visco-poroélastiques), je peux calculer la propagation des ondes sismiques sans approximation. En effet, j'utilise une méthode numérique dans le domaine fréquence-espace qui permet de prendre en compte tous les termes qui dépendent de la fréquence sans approximation. La discrétisation spatiale utilise une méthode d'éléments finis discontinus (Galerkin discontinu) qui permet de considérer des milieux hétérogènes.Je montre notamment que les attributs sismiques (vitesses et atténuations) des milieux poreux complexes sont fortement dispersifs et les formes d'ondes complètes, calculées sans approximation, sont fortement dépendantes de la description physique du milieu. La caractérisation des paramètres poroélastiques s'effectue par inversion. Une méthode en deux étapes a été proposée : la première consiste en une inversion ``classique`` (tomographie, inversion des formes d'ondes complètes) des données (sismogrammes) pour obtenir des paramètres macro-échelles (attributs sismiques). La seconde étape permet de reconstruire, à partir des paramètres macro-échelles, les paramètres poroélastiques micro-échelles. Cette étape d'inversion utilise une méthode d'optimisation semi-globale (algorithme de voisinage). Une analyse de sensibilité montre qu'en connaissant a-priori certains paramètres, on peut inverser avec précision les paramètres du squelette poroélastique ou retrouver la nature du fluide saturant, à partir des vitesses de propagation. En revanche, pour retrouver la saturation en fluide, il est préférable de connaître les atténuations. Deux applications réalistes (monitoring de réservoir et hydrogéophysique) mettent en oeuvre ce type d'inversion en deux étapes et démontrent qu'à partir de données estimées par des méthodes classiques d'imagerie, on peut remonter à certains paramètres poroélastiques constitutifs. / Seismic wave propagation in multiphasic porous media have various environmental (natural risks, geotechnics, groundwater pollutions...) and ressources (aquifers, oil and gas, CO2 storage...) issues. When seismic waves are crossing a given material, they are distorted and thus contain information on fluid and solid phases. This work focuses on the characteristics of seismic waves propagating in multiphasic media, from the physical complex description to the parameter characterisation by inversion, including 2D numerical modelling of the wave propagation. The first part consists in the description of the physics of multiphasic media (each phase and their interactions), using several upscaling methods, in order to obtain an equivalent mesoscale medium defined by seven parameters. Thus, in simple porosity saturated media and in complex media (double porosity, patchy saturation, visco-poroelasticity), I can compute seismic wave propagation without any approximation. Indeed, I use a frequency-space domain for the numerical method, which allows to consider all the frequency dependent terms. The spatial discretisation employs a discontinuous finite elements method (discontinuous Galerkin), which allows to take into account complex interfaces.The computation of the seismic attributes (velocities and attenuations) of complex porous media shows strong variations in respect with the frequency. Waveforms, computed without approximation, are strongly different if we take into account the full description of the medium or an homogenisation by averages. The last part of this work deals with the poroelastic parameters characterisation by inversion. For this, I develop a two-steps method: the first one consists in a classical inversion (tomography, full waveform inversion) of seismograms data to obtain macro-scale parameters (seismic attributes). The second step allows to recover, from the macroscale parameters, the poroelastic micro-scale properties. This downscaling step uses a semi-global optimisation method (neighbourhood algorithm), which allows the sampling of the full model space (thanks to the low numerical cost of the analytic direct model). With the a-priori knowledge of some parameters, a sensibility analysis shows that I can invert precisely skeleton parameters or the saturating fluid type, from the velocities only. Nevertheless, to recover the fluid saturation, it is preferable to use the attenuations. This two-steps procedure is tested on two realistic applications (reservoir monitoring and subsurface hydrogeophysics) and show that we can recover some constituve poroelastic parameters.

Dans la nature les plantes sont sans cesse soumises à des sollicitations mécaniques qui affectent et modifient leur croissance. Un aspect remarquable de cette réponse est qu’elle n’est pas seulement locale mais non-locale : la flexion d’une tige ou d’une branche inhibe rapidement la croissance loin de la zone sollicitée. Cette observation suggère l'existence d'un signal pouvant se propager à travers toute la plante. Parmi les différentes hypothèses, il a été suggéré que ce signal pouvait être purement mécanique, et provenir d’un couplage hydro/mécanique entre la déformation du tissu et la pression de l’eau contenue dans le système vasculaire de la plante. L’objectif de cette thèse est de comprendre l’origine physique de ce couplage par une approche biomimétique. Pour cela, nous avons développé des branches artificielles micro-fluidiques possédant des caractéristiques mécaniques et hydrauliques similaires à celles d'une branche d'arbre. Nous avons montré que la flexion de ces branches génère une surpression globale non-nulle dans le système, qui varie comme le carré de la déformation longitudinale. Un modèle simple basé sur un mécanisme analogue à l’ovalisation des tubes permet de prédire cette réponse poroélastique non-linéaire et d’identifier le paramètre physique clé pilotant cette réponse en pression : le module de compressibilité de la branche. A la lumière de ces résultats, des expériences sur des branches d'arbre ont ensuite été conduites et des signaux similaires sont obtenus et comparés au modèle théorique. La similitude suggère le caractère générique du mécanisme physique identifié pour la génération de signaux hydraulique dans les plantes. / Plants are constantly subjected to external mechanical loads such as wind or touch and respond to these stimuli by modifying their growth and development. A fascinating feature of this mechanical-induced-growth response is that it is not only local, but also non-local: bending locally a stem or a branch can induce a very rapid modification of the growth far away from the stimulated area, suggesting the existence of a signal that propagates across the whole plant. The nature and origin of this signal is still not understood, but it has been suggested recently that it could be purely mechanical and originate from the coupling between the local deformation of the tissues and the water pressure in the vascular system. The objective of this work is to understand the origin of this hydro/mechanical coupling using a biomimetic approach. Artificial microfluidic branches have been developed, that incorporate the mechanical and hydraulic key features of natural ones. We show that the bending of these branches generates a steady overpressure in the whole system, which varies quadratically with the bending deformation. A simple model based on a mechanism analogue to tube ovalization enables us to predict this non-linear poroelastic response, and identify the key physical parameter at play, namely the elastic bulk modulus of the branch. Further experiments conducted on natural tree branches reveal the same phenomenology. Once rescaled by the model prediction, both the biomimetic and natural branches falls on the same master curve, showing the universality of the identified mechanism for the generation of hydraulic signals in plants.

Les matériaux absorbants ont une efficacité très faible aux basses fréquences (<500Hz). Les mousses dites adaptatives ont été développées pour combler ce manque d'efficacité. Une mousse adaptative est la réunion d'un matériau absorbant et d'un actionneur commandé par un système de contrôle actif. L'objectif de cette thèse est de mener une étude numérique et expérimentale approfondie du concept de mousse adaptative et d'en dégager les mécanismes physiques et les limitations technologiques pour le contrôle de l'absorption acoustique. Le matériau absorbant est ici une mousse de mélamine et l'actionneur est un film piézoélectrique de PVDF. Le contrôle actif de l'absorption acoustique est réalisé en incidence normale avec l'hypothèse d'onde plane sur la plage de fréquence [0-1500Hz]. Les résultats font apparaître la possibilité d'absorber une pression incidente de 1Pa à 100Hz avec 100V et un bruit large bande de 94dB avec une centaine de Vrms à partir de 250Hz. Ces résultats ont été obtenus avec un prototype de mousse adaptative ayant une épaisseur moyenne de 4cm. Un frein important au contrôle large bande provient du fort niveau de distorsion des prototypes dans le bas (<500Hz) et le haut (<1500Hz) du spectre de fréquences. Un modèle éléments finis 3D couplant les domaines poroélastiques, acoustiques, élastiques et piézoélectriques a été proposé et validé expérimentalement. L'utilisation de ce modèle, complété par une étude analytique, a permis de mettre en lumière le mode d'action et les mécanismes de dissipation dans la mousse adaptative.

Cette thèse propose une étude du comportement des tissus mous intra-orbitaires dans le cadre du traitement des exophtalmies, caractérisées par une protrusion du globe oculaire. La partie chirurgicale de ce traitement consiste en une décompression de l'orbite par le biais d'une ostéotomie des parois orbitaires pour réduire la protrusion du globe. Nos travaux visent à prédire les relations entre le volume des tissus décompressés, la surface d'ostéotomie et le recul oculaire résultant pour apporter une aide au planning chirurgical. Pour cela, deux modèles ont été développés. Le premier est un modèle analytique simple de l'orbite assimilant les parois orbitaires à un cône et le globe à une sphère. Il permet de déterminer de façon satisfaisante le volume décompressé en fonction du recul souhaité. Le deuxième est un modèle biomécanique des tissus mous intra-orbitaires et de leurs interactions avec les parois osseuses et le globe. Il s'agit d'un maillage Eléments Finis utilisant un matériau poroélastique et prenant en compte la morphologie de l'orbite du patient, les propriétés mécaniques des tissus mous. Il permet de quantifier le recul oculaire et le volume de tissus décompressés en fonction de l'effort (ou du déplacement) imposé par le chirurgien et de l'ostéotomie (surface et position). Ses résultats sont relativement intéressants et permettent d'évaluer le comportement des tissus intra-orbitaires. La méthode de génération automatique de maillage éléments finis, développée dans cette thèse, a permis d'effectuer différentes simulations d'ostéotomies et de conclure que la morphologie de l'orbite patient a un impact sur le recul et le volume décompressé et que l'influence de la surface d'ostéotomie est modérée. Une analyse rhéologique des tissus mous orbitaires a été effectuée. Les tests in vitro ont permis de poser les bases de futures mesures. Les tests in vivo, réalisés avec un capteur ad hoc, ont déterminé la raideur des tissus mous. Bien que les modèles présentés dans cette thèse ne soient pas encore utilisables dans le cadre d'un planning chirurgical, ils fournissent des résultats satisfaisants et une bonne estimation des phénomènes observés lors d'une décompression orbitaire.

La propagation des ondes sismiques dans les milieux poreux multiphasiques présente des enjeux nombreux, tant sur le plan environnemental (risques naturels, géotechnique, pollutions de nappes...) que pour les réservoirs (aquifères, hydrocarbures, stockages de CO2...). L'utilisation des ondes sismiques pour étudier ces milieux se justifie par le fait qu'en se propageant, les ondes sont déformées par le milieu qu'elles traversent et contiennent ainsi des informations aux capteurs sur les phases fluides et solides et sur le squelette poreux. Ce travail de thèse s'intéresse aux caractéristiques des ondes sismiques dans les milieux multiphasiques (plusieurs phases fluides et solides), depuis la description physique jusqu'à la caractérisation des paramètres constitutifs par inversion, en passant par la modélisation numérique 2D de la propagation. La première partie du travail a consisté à décrire la physique des milieux multiphasiques (phase par phase et leurs intéractions dynamiques) en utilisant des méthodes d'homogénéisation pour se ramener à un milieu équivalent défini par sept paramètres. Ainsi, dans des milieux simple porosité saturés et dans des milieux plus complexes (double porosité, partiellement saturés ou visco-poroélastiques), je peux calculer la propagation des ondes sismiques sans approximation. En effet, j'utilise une méthode numérique dans le domaine fréquence-espace qui permet de prendre en compte tous les termes qui dépendent de la fréquence sans approximation. La discrétisation spatiale utilise une méthode d'éléments finis discontinus (Galerkin discontinu) qui permet de considérer des milieux hétérogènes.Je montre notamment que les attributs sismiques (vitesses et atténuations) des milieux poreux complexes sont fortement dispersifs et les formes d'ondes complètes, calculées sans approximation, sont fortement dépendantes de la description physique du milieu. La caractérisation des paramètres poroélastiques s'effectue par inversion. Une méthode en deux étapes a été proposée : la première consiste en une inversion ''classique'' (tomographie, inversion des formes d'ondes complètes) des données (sismogrammes) pour obtenir des paramètres macro-échelles (attributs sismiques). La seconde étape permet de reconstruire, à partir des paramètres macro-échelles, les paramètres poroélastiques micro-échelles. Cette étape d'inversion utilise une méthode d'optimisation semi-globale (algorithme de voisinage). Une analyse de sensibilité montre qu'en connaissant a-priori certains paramètres, on peut inverser avec précision les paramètres du squelette poroélastique ou retrouver la nature du fluide saturant, à partir des vitesses de propagation. En revanche, pour retrouver la saturation en fluide, il est préférable de connaître les atténuations. Deux applications réalistes (monitoring de réservoir et hydrogéophysique) mettent en oeuvre ce type d'inversion en deux étapes et démontrent qu'à partir de données estimées par des méthodes classiques d'imagerie, on peut remonter à certains paramètres poroélastiques constitutifs.

Les roches argileuses profondes de très faible perméabilité (10-20 m2), telles que l’argilite du Callovo-Oxfordien (COx) en France ou l’Argile à Opalinus en Suisse, sont des roches hôtes potentielles pour le stockage géologique des déchets radioactifs. Lors des différentes phases du stockage, ces roches seront soumises à des sollicitations thermo-hydro-mécaniques (THM) couplées. La détermination de leurs caractéristiques THM reste à compléter pour une meilleure compréhension de la réponse du champ proche des galeries de stockage.L’étude expérimentale des matériaux de faible perméabilité est délicate et plusieurs questions sur leur comportement THM restent posées. Dans ce contexte, deux systèmes expérimentaux originaux avec mesures précises des déformations locales et chemin de drainage (H) réduit ont été développés et utilisés pour la réalisation d’essais saturés drainés : une cellule de compression isotrope (H = 10 mm) et une cellule triaxiale standard dont le chemin de drainage a été réduit à l’aide d’un géotextile placé autour de l’échantillon (H = 19 mm) .Une caractérisation poroélastique détaillée de l’argilite du COx en cellule isotrope a permis, par des approches directes et indirectes compatibles, une détermination fiable des paramètres poroélastiques isotropes transverses du matériau, identifiés dans un cadre théorique permettant de déterminer les composantes du tenseur de Biot (coefficients b1 et b2, respectivement perpendiculaire et parallèle au litage). Un paramètre clé peu documenté à ce jour est le module de compression des grains solides, déterminé à l’aide d’un essai de compression sans membrane (Ks = 21.7 GPa).Un contrôle de température et un étalonnage précis de ces effets parasites ont permis l’étude de la réponse volumique thermique de l’argilite du COx lors d’un essai de chauffage drainé sous confinement isotrope constant proche de l’état in-situ. Une expansion thermoélastique suivie d’une contraction thermoplastique a été observée, avec une température de transition à 48°C, proche de la plus grande température supportée pendant l’histoire géologique de l’argilite. On confirme ainsi le fait que les argilites gardent en mémoire la température maximale supportée. La précision des mesures de déformations a aussi permis d’identifier un fluage volumique, dont l’amplitude est accentuée à 80°C.Des calculs poroélastiques avec les conditions aux limites du système de drainage amélioré ont permis de déterminer l’ordre de grandeur des taux de déformations axiales permettant un bon drainage lors du cisaillement triaxial (6.6×10-8 s-1). Les essais drainés sur l’argilite du COx présentent une bonne compatibilité avec un critère déjà publié. Des résultats cohérents ont aussi été obtenus sur le critère de rupture de l’Argile à Opalinus, par rapport auquel ont été comparés ceux d’essais de compression simple à différentes vitesses sur des échantillons équilibrés à 94% d’humidité relative qui a mis en évidence une dépendance de la résistance au pic vis-à-vis du taux de déformation axial.L’ensemble de ces résultats permet de réduire l’incertitude concernant les propriétés THM des argilites et devrait permettre une meilleure estimation de la réponse du champ proche des galeries au cours des différentes phases de leur période service.Mot clés : argilites, perméabilité, drainage, saturation, poroélasticité, coefficients de Biot, isotropie transverse, cisaillement triaxial, chauffage, fluage / Deep low permeability claystones (10-20 m2), such as Callovo-Oxfordian claystone (COx) in France or the Opalinus Clay in Switzerland, are potential host rocks for deep geological radioactive waste disposal. During the various phases of the storage, these rocks will be subjected to thermo-hydro-mechanical (THM) coupled effects. The determination of their THM parameters remains to be completed for a better understanding of the near-field response of the storage galleries.The experimental study of low permeability geomaterials is difficult and several questions about their THM behavior still remain. In this context, two original experimental systems with high precision local strain measurements and reduced drainage lengths (H) were developed and used for saturated drained tests: an isotropic compression cell (H = 10 mm) and a standard triaxial cell with a reduced drainage length using a geotextile placed around the sample (H = 19 mm).A detailed poroelastic characterization of the COx argillite in the isotropic cell provided a set of compatible transverse isotropic poroelastic parameters of the material, identified in a theoretical framework allowing to determine the Biot tensor components (coefficients b1 and b2, perpendicular and parallel to bedding plans, respectively). A key parameter not well documented to date is the unjacketed modulus that was determined by means of an unjacketed compression test (Ks = 21.7 GPa).A temperature control and an accurate calibration of thermal parasite effects allowed the investigation of the thermal volumetric response of the COx argillite during a drainage test under constant isotropic confining stress close in-situ state conditions. A thermoelastic expansion followed by a thermoplastic contraction was observed, with a transition at a temperature of 48 °C, close to the highest temperature supported during the geological history of the claystone. It is thus confirmed that such claystones keep in memory the maximum supported geological temperature. The precision of the deformation measurements also made it possible to identify a volumetric creep that is enhanced at 80 °C.Poroelastic calculations with the boundary conditions of the improved drainage system allowed to determine the magnitude of axial strain rates allowing good drainage during triaxial drained shearing (6.6×10-8 s-1). The drained tests carried out on the COx claystone showed a good compatibility with a criterion already published. Coherent data were also obtained on the Opalinus Clay failure criterion, that were compared to those of uniaxial compression tests at different speeds on samples equilibrated at 94% relative humidity, that exhibited a dependence of the peak strength on the shear rate.These results make it possible to reduce the uncertainties concerning the THM properties of claystones and should allow a better estimation of the response of the near field close to the galleries during the different phases of their service period.Key words: claystone, permeability, drainage, saturation, poroelasticity, Biot coefficient, transverse isotropy, triaxial testing, heating, creep

Une connaissance plus approfondie du comportement des roches argileuses, et plus particulièrement de leur tenue, s'avère aujourd'hui indispensable dans le cadre de plusieurs thèmes de recherche d'importance croissante : comme couverture des réservoirs pétroliers pour la production d’hydrocarbures et la séquestration de CO2, ou en tant que roche hôte pour le stockage des déchets nucléaires en couches géologiques profondes. Ces roches sont particulièrement difficiles à caractériser expérimentalement du fait de leur très faible perméabilité et de leur sensibilité à l'eau. L'objectif de cette thèse est d'améliorer la compréhension des couplages hydromécaniques intervenant dans ces matériaux. Tout d’abord, une contribution à la caractérisation de l’évolution de la structure des roches argileuses soumises à un chargement mécanique est apportée. A cet effet, un essai de compression uniaxiale sur l’argilite de Tournemire associant des observations au MEB et une méthode de corrélation d’images numériques est proposé. Par la suite, des simulations numériques des essais œdométriques ont été réalisées pour appréhender l'influence des propriétés poromécaniques, en particulier de la faible perméabilité sur le comportement de l'échantillon. L’étude des équations de diffusion et les simulations numériques effectuées nous permettent de proposer une méthode de détermination de la perméabilité intrinsèque de ces matériaux. Dans la dernière partie, une étude du comportement d’une roche couverture d’un réservoir pétrolier à grande échelle est réalisée sur un cas synthétique 3D. L’objectif de ce chapitre est d’étudier les variations du champ de contrainte au sein de la couverture au cours de la production du réservoir / Knowledge of shale poromechanical behavior is essential for various environmental issues such as CO2 deep geological storage, high level radioactive waste storage, oil field abandonment… Furthermore, due to their good confinement properties, there are considered as potential host rocks for high level radioactive waste storage. However, their experimental characterization is very complex because of their very low permeability and their saturating fluid sensitivity. The main objective of this study is to improve knowledge of hydromechanical coupling of such a material. To study the structure of argillites under stress, SEM in-situ uniaxial compression test is performed on parallelepiped sample; compression orientation is perpendicular to bedding planes direction. Correlation techniques are used to estimate strain distributions. Second part underlines the key points of shale experimental characterization using the framework of Biot’s mechanics of fluid saturated porous solids. Shales are well known as a more or less transverse isotropy material. The complete methodology to conduct œdometric tests on such sensitive and weakly permeable material is described. Measurements realized on Tournemire argillite are proposed and a comparison with the poroelastic parameters of Meuse / Haute-Marne is given. Furthermore, to investigate hydromechanical coupling, a transverse isotropic poroelastic model is used to estimate the influence of anisotropy on diffusion characteristic time. The experimental data are used in numerical simulations to estimate the intrinsic permeability of the samples. In the last part, the caprock behavior of a petroleum field is studied at large scale with a 3D synthetic model (finite element analysis). The purpose of this last chapter is to study changes in stress field during the oil production

Description analytique de la houle, effet de la déformabilité et de la perméabilité du fond marin sous la propagation de la houle, valeur de la pression cyclique induite sur le fond marin. Développement analytique pour l'effet de porosité et de perméabilité, propagation sur un fond incliné. Définition des contraintes induites et analyse de la réponse du sol sous ces sollicitations. Modélisation de l'action d'une houle sinusoïdale sur un fond marin élastique. Expansions analytiques des contraintes par la transformation de MELLIN. Calcul quasi statique de l'évolution d'un milieu poreux. Traitement numérique du problème par éléments finis. Etude paramétrique. La connaissance du comportement des massifs sous marins est devenue nécessaire au développement des travaux en mer.