Diese Arbeit stellt das Konzept der kompressionssensitiven Magnetresonanzelastographie vor. Kompressionssensitive MRE analysiert die Ausbreitung von Kompressionswellen und liefert dadurch Erkenntnisse über die Kompressionseigenschaften eines Mediums auf Grundlage eines poroelastischen Modells. Anomalien bei der Regulation des Gewebedrucks stehen in Zusammenhang mit verschiedenen Krankheitsbildern, wie Normaldruck-Hydrozephalus und Pfortader-Hypertonie. Statischer Druck spielt als Porendruck eine zentrale Rolle in den poroelastischen Wellengleichungen; die kompressionssensitive MRE könnte daher ein nichtinvasives Diagnoseinstrument darstellen, das die durch konventionelle Scherwellen-Elastographie gewonnenen Informationen um weitere Aspekte ergänzt. Diese Arbeit beschreibt die Entwicklung einer schnellen Singleshot-EPI-Bildgebungssequenz, mit deren Hilfe die durch propagierende Druckwellen hervorgerufene volumetrische Verzerrung quantifiziert werden kann. Die Validierung der kompressionssensitiven MRE erfolgte an verschiedenen Systemen: an porösen Gelphantomen, an der menschlichen Lunge in zwei Atemzuständen, in einer ex-vivo Schafsleber bei unterschiedlichen hydrostatischen Drücken und schließlich am menschlichen Gehirn. Die Ergebnisse belegen, dass die Stärke der induzierten volumetrischen Verzerrung sensitiv gegenüber Druckänderungen ist, wohingegen die Scherverzerrung keine derartige Abhängigkeit aufweist. In einer weiteren Studie wurde intrinsische Pulsation des menschlichen Hirns anstelle einer externen Vibrationsquelle ausgenutzt. Dabei erzeugte die arterielle Pulswelle eine kurze lokale Expansion des Hirnparenchyms; in der sich anschließenden diastolischen Phase erfolgte eine langsame Rückkehr zum Ausgangszustand. Aus den gemessenen volumetrischen Verzerrungen wurden durch Inversion der Druckwellengleichung numerische Werte für den Druckwellenmodul M berechnet; Rauschen wurde als primäre Ursache für die systematische Unterschätzung von M identifiziert. / This thesis introduces the concept of compression-sensitive Magnetic Resonance Elastography. Compression-sensitive MRE detects the propagation of pressure waves, providing insight into the compressibility of a material based on a poroelastic tissue model. Poroelastic models incorporate compressibility through interaction of compartments, even as each individual compartment remains incompressible. Hydrostatic tissue pressure abnormalities are associated with a number of diseases, such as normal pressure hydrocephalus or hepatic portal hypertension. Since pore pressure plays a central role in the poroelastic wave equations, compression-sensitive MRE could potentially serve as a diagnostic tool, providing information complimentary to shear-wave MRE data. This thesis describes the development of a fast single-shot EPI MR sequence capable of quantifying volumetric strain induced by external vibrations. Compression-sensitive MRE was validated in porous gel phantoms, in the human lung at two different respiratory states, in an ex vivo sheep liver at varying levels of hydrostatic pressure, and finally in human liver and brain. Results illustrate that compression-sensitive MRE is capable of quantifying volumetric strain in phantoms and in human organs. It was found that volumetric strain was sensitive toward pressure changes associated with different physiological states, whereas shear strain remained constant. In an additional study, pulsation of the human brain, driven by the heart cycle, was used as the actuation source instead of the external vibration generator. Results indicate local expansion of brain parenchyma upon the arrival of the arterial pulse wave, followed by a slow return to the initial state during the diastolic phase. Numerical values for the pressure wave modulus M were calculated from measured volumetric strain through inversion of the pressure wave equation. Measurement noise was identified as the primary effect causing a severe underestimation of M.

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.

Submitted by Maria Cristina (library@lncc.br) on 2017-12-12T10:49:35Z No. of bitstreams: 1 Thesis - Thiago Quinelato.pdf: 2369263 bytes, checksum: 6a1ac9e2d37bb0377981785cfa50683c (MD5) / Approved for entry into archive by Maria Cristina (library@lncc.br) on 2017-12-12T10:50:02Z (GMT) No. of bitstreams: 1 Thesis - Thiago Quinelato.pdf: 2369263 bytes, checksum: 6a1ac9e2d37bb0377981785cfa50683c (MD5) / Made available in DSpace on 2017-12-12T10:50:14Z (GMT). No. of bitstreams: 1 Thesis - Thiago Quinelato.pdf: 2369263 bytes, checksum: 6a1ac9e2d37bb0377981785cfa50683c (MD5) Previous issue date: 2017-03-01 / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) / Esta tese é focada no desenvolvimento e na análise de aproximações em dimensão finita das equações que descrevem problemas de elasticidade linear e poroelasticidade. A estratégia de aproximação é baseada em formulações de elementos finitos mistas hibridas desses problemas e a construção dos espaços de dimensão finita é guiada por várias propriedades desejadas: continuidade das trações (conservação do momento linear), simetria do tensor de tensão (conservação do momento angular), número reduzido de graus de liberdade globais e robustez sob distorção de malha. A principal dificuldade está relacionada com o atendimento simultâneo da condição inf-sup e da simetria do tensor de tensão. O ultimo requisito é relaxado, sendo satisfeito de maneira fraca pela introdução de um multiplicador de Lagrange. A maior contribuição é o desenvolvimento e a análise de espaços de dimensão finita estáveis para aproximação mista dos problemas de elasticidade linear e poroelasticidade em malhas quadrilaterais arbitrárias. Esses espaços são capazes de fornecer convergência com taxa ótima do campo de tensão na norma H(div) em malhas de quadriláteros arbitrários, o que é provado pela análise numérica e confirmado por experimentação. / This thesis is focused on the development and analysis of finite dimensional approximations of the equations describing linear elasticity and poroelasticity problems. The approximation strategy is based on mixed hybrid finite element formulations of those problems and the construction of the finite dimensional spaces is guided by several desired properties: continuity of the tractions (conservation of linear momentum), symmetry of the stress tensor (conservation of angular momentum), reduced number of global degrees of freedom, and robustness under mesh distortion. The main difficulty is related with the simultaneous fulfillment of the inf-sup condition and the symmetry of the stress tensor. The last requirement is relaxed, being enforced in the weak sense through the introduction of a Lagrange multiplier. The main contribution is the development and analysis of stable finite dimensional spaces for mixed approximation of linear elasticity and poroelasticity problems on arbitrary quadrilateral meshes. These spaces are capable of providing optimal order convergence of the stress field in the H(div)-norm on meshes of arbitrary quadrilaterals, which is proved by numerical analysis and confirmed by experimentation.

O conhecimento de zonas potencialmente instáveis no fundo do mar é de fundamental importância para o desenvolvimento das estruturas marinhas, pois permite posicionar estruturas offshore em áreas mais seguras, reduzindo-se possíveis danos, custos e eventual poluição ambiental. Nesse contexto, o objetivo do presente trabalho é investigar, através de uma abordagem analítico-numérica, a estabilidade de maciços submarinos submetidos ao carregamento cíclico da onda. O efeito das ondas de água sobre o leito submerso é descrito pela propagação de uma onda de pressão ao longo de sua superfície, empregando-se a teoria linear de Stokes. São considerados maciços com superfície superior horizontal e inclinada, constituídos por material coesivo (argilas) e material granular (areias). Em maciços constituídos por solos finos, a capacidade resistente do material é modelada pelo critério de Tresca não-homogêneo e a análise da estabilidade é desenvolvida em condição não drenada. Por outro lado, em leitos granulares, a resistência do meio depende explicitamente do valor da poropressão, sendo descrita classicamente pelo critério de Coulomb sem coesão. A análise de estabilidade é então desenvolvida em tensões efetivas e o gradiente de poropressão atua como uma força volumétrica sobre o esqueleto, caracterizando o modo de carregamento principal deste material. Em razão da tendência a se densificar quando submetido a um estado de tensões desviadoras cíclicas, ocorre, em geral, a acumulação de excesso de poropressão no maciço granular. Consequentemente, o acoplamento entre o comportamento do material e o carregamento cíclico tem fundamental importância sobre o cálculo do excesso de poropressão desenvolvido. Para determinação das forças de percolação, considera-se uma abordagem simplificada baseada na partição das deformações em contribuições reversível e irreversível, que permite desacoplar o cálculo da pressão intersticial induzida pela onda. Aplicando-se conceitos da teoria da Análise Limite é possível formular limites inferiores e superiores da máxima amplitude segura do carregamento da onda. Finalmente, os efeitos da declividade da superfície do leito e da espessura de camada de solo sobre a estabilidade são analisados. / The knowledge of potentially unstable areas on the seabed is of fundamental importance to the development of marine structures, because it allows to install offshore structures in safer areas, reducing possible damages, costs and eventual environmental pollution. In this context, the objective of the present work is to investigate the stability of submarine soil masses subjected to the wave cyclic loading through analytical-numerical approaches. The effect of water waves on the submerged bed is described by the propagation of a pressure wave along its surface, using the Stokes’s linear theory. Soil masses with horizontal and sloped upper surface, composed of cohesive material (clays) and granular materiais (sands) are considered in this study. In soil masses constituted of fine soil, the material strength capacity is modeled by the non-homogeneous Tresca criterion and the stability analysis is carried out in undrained condition. On the other hand, in granular beds, the strength explicitly depends on the pore pressure value, being classically described by the Coulomb criterion without cohesion. Then, the stability analysis is developed in effective stress and the pore pressure gradient acts as a volumetric force on the skeleton, characterizing the main charging mode of this material. Due to the tendency to densify when subjected to a cyclic deviatoric stress state occurs, in general, the build-up of pore pressure excess in the granular mass. Consequently, the coupling between the material behavior and the cyclic loading has fundamental importance in the calculation of the pore pressure excess generated. In order to define the seepage forces, a simplified approach based on the partition of deformations in reversible and irreversible contributions is considered, which allows to decouple the wave-induced pore pressure calculation. Applying the concepts of the limit analysis theory it is possible to formulate upper and lower boundaries of the maximum safe amplitude of the wave loading. Finally, the effects of the seabed surface steepness and of the soil layer thickness on the stability are analyzed.

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.

[pt] A grande maioria dos problemas de estabilidade de poços de petróleo ocorre em trechos de folhelhos, rochas nas quais, uma especificação eficiente da pressão do fluido de perfuração requer previamente uma especificação correta da concentração salina e da temperatura. Todavia, para um dimensionamento adequado das características do fluido de perfuração necessárias à estabilidade do poço, é necessário o uso de modelos matemáticos que considerem um acoplamento adequado entre efeitos poroelásticos, químicos e térmicos. Entretanto, a complexidade matemática das equações de modelos acoplados normalmente leva à adoção de soluções numéricas, que consomem um tempo computacional muito grande e, por isso, esses modelos não são atrativos à aplicação na análise da estabilidade de poços. Este trabalho apresenta um modelo acoplado termo-químico-poroelástico representado por duas soluções, uma numérica, que utiliza o método dos elementos finitos, e outra analítica, baseada no método das transformadas de Laplace. Ao comparar ambas as soluções é demonstrado que a solução analítica consegue representar tão bem quanto à solução numérica os principais processos acoplados presentes durante a perfuração de folhelhos e que influenciam na sua estabilidade e, por esta razão, pode ser utilizada na análise de estabilidade de poços em folhelhos. Através de um estudo de caso, é verificado que um controle eficiente da estabilidade do poço é obtido especificando a pressão do fluido de perfuração em função da sua temperatura e concentração salina. Estes resultados também indicam as razões de alguns problemas não previstos por modelos desacoplados, e que quase sempre ocorrem durante a perfuração em folhelhos. / [en] Wellbore stability problems are most common when drilling through shales. In order to avoid such problems in this kind of rocks the solute concentration and temperature must be properly defined in the drilling fluid composition, which requires considering poroelastic, thermal and chemical effects in a coupled way. The equations complexity of coupled models usually results in numerical solutions that are very time consuming, thus, unattractive for stability analysis. In an opposite way, it is very difficult to develop closed- form solutions for coupled models. This work presents a thermochemoporoelastic model represented by a numerical solution based upon the finite element method and an analytical solution based upon the Laplace transform method. A comparison between the results of the numerical solution and analytical solution shows that the later can reproduce the coupled processes involved in the wellbore stability problem in shales as well as the former, and for this reason the closed-form solution can be applied as a practical tool in wellbore stability analysis. The analysis of a typical wellbore drilled through shales showed that an efficient control of wellbore stability can be obtained through an adequate specification of the drilling fluid pressure when taking in account its solute concentration and temperature. The model was also able to explain some problems not predicted by uncoupled models, but almost always seen while drilling through shales.

O conhecimento de zonas potencialmente instáveis no fundo do mar é de fundamental importância para o desenvolvimento das estruturas marinhas, pois permite posicionar estruturas offshore em áreas mais seguras, reduzindo-se possíveis danos, custos e eventual poluição ambiental. Nesse contexto, o objetivo do presente trabalho é investigar, através de uma abordagem analítico-numérica, a estabilidade de maciços submarinos submetidos ao carregamento cíclico da onda. O efeito das ondas de água sobre o leito submerso é descrito pela propagação de uma onda de pressão ao longo de sua superfície, empregando-se a teoria linear de Stokes. São considerados maciços com superfície superior horizontal e inclinada, constituídos por material coesivo (argilas) e material granular (areias). Em maciços constituídos por solos finos, a capacidade resistente do material é modelada pelo critério de Tresca não-homogêneo e a análise da estabilidade é desenvolvida em condição não drenada. Por outro lado, em leitos granulares, a resistência do meio depende explicitamente do valor da poropressão, sendo descrita classicamente pelo critério de Coulomb sem coesão. A análise de estabilidade é então desenvolvida em tensões efetivas e o gradiente de poropressão atua como uma força volumétrica sobre o esqueleto, caracterizando o modo de carregamento principal deste material. Em razão da tendência a se densificar quando submetido a um estado de tensões desviadoras cíclicas, ocorre, em geral, a acumulação de excesso de poropressão no maciço granular. Consequentemente, o acoplamento entre o comportamento do material e o carregamento cíclico tem fundamental importância sobre o cálculo do excesso de poropressão desenvolvido. Para determinação das forças de percolação, considera-se uma abordagem simplificada baseada na partição das deformações em contribuições reversível e irreversível, que permite desacoplar o cálculo da pressão intersticial induzida pela onda. Aplicando-se conceitos da teoria da Análise Limite é possível formular limites inferiores e superiores da máxima amplitude segura do carregamento da onda. Finalmente, os efeitos da declividade da superfície do leito e da espessura de camada de solo sobre a estabilidade são analisados. / The knowledge of potentially unstable areas on the seabed is of fundamental importance to the development of marine structures, because it allows to install offshore structures in safer areas, reducing possible damages, costs and eventual environmental pollution. In this context, the objective of the present work is to investigate the stability of submarine soil masses subjected to the wave cyclic loading through analytical-numerical approaches. The effect of water waves on the submerged bed is described by the propagation of a pressure wave along its surface, using the Stokes’s linear theory. Soil masses with horizontal and sloped upper surface, composed of cohesive material (clays) and granular materiais (sands) are considered in this study. In soil masses constituted of fine soil, the material strength capacity is modeled by the non-homogeneous Tresca criterion and the stability analysis is carried out in undrained condition. On the other hand, in granular beds, the strength explicitly depends on the pore pressure value, being classically described by the Coulomb criterion without cohesion. Then, the stability analysis is developed in effective stress and the pore pressure gradient acts as a volumetric force on the skeleton, characterizing the main charging mode of this material. Due to the tendency to densify when subjected to a cyclic deviatoric stress state occurs, in general, the build-up of pore pressure excess in the granular mass. Consequently, the coupling between the material behavior and the cyclic loading has fundamental importance in the calculation of the pore pressure excess generated. In order to define the seepage forces, a simplified approach based on the partition of deformations in reversible and irreversible contributions is considered, which allows to decouple the wave-induced pore pressure calculation. Applying the concepts of the limit analysis theory it is possible to formulate upper and lower boundaries of the maximum safe amplitude of the wave loading. Finally, the effects of the seabed surface steepness and of the soil layer thickness on the stability are analyzed.

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

O conhecimento de zonas potencialmente instáveis no fundo do mar é de fundamental importância para o desenvolvimento das estruturas marinhas, pois permite posicionar estruturas offshore em áreas mais seguras, reduzindo-se possíveis danos, custos e eventual poluição ambiental. Nesse contexto, o objetivo do presente trabalho é investigar, através de uma abordagem analítico-numérica, a estabilidade de maciços submarinos submetidos ao carregamento cíclico da onda. O efeito das ondas de água sobre o leito submerso é descrito pela propagação de uma onda de pressão ao longo de sua superfície, empregando-se a teoria linear de Stokes. São considerados maciços com superfície superior horizontal e inclinada, constituídos por material coesivo (argilas) e material granular (areias). Em maciços constituídos por solos finos, a capacidade resistente do material é modelada pelo critério de Tresca não-homogêneo e a análise da estabilidade é desenvolvida em condição não drenada. Por outro lado, em leitos granulares, a resistência do meio depende explicitamente do valor da poropressão, sendo descrita classicamente pelo critério de Coulomb sem coesão. A análise de estabilidade é então desenvolvida em tensões efetivas e o gradiente de poropressão atua como uma força volumétrica sobre o esqueleto, caracterizando o modo de carregamento principal deste material. Em razão da tendência a se densificar quando submetido a um estado de tensões desviadoras cíclicas, ocorre, em geral, a acumulação de excesso de poropressão no maciço granular. Consequentemente, o acoplamento entre o comportamento do material e o carregamento cíclico tem fundamental importância sobre o cálculo do excesso de poropressão desenvolvido. Para determinação das forças de percolação, considera-se uma abordagem simplificada baseada na partição das deformações em contribuições reversível e irreversível, que permite desacoplar o cálculo da pressão intersticial induzida pela onda. Aplicando-se conceitos da teoria da Análise Limite é possível formular limites inferiores e superiores da máxima amplitude segura do carregamento da onda. Finalmente, os efeitos da declividade da superfície do leito e da espessura de camada de solo sobre a estabilidade são analisados. / The knowledge of potentially unstable areas on the seabed is of fundamental importance to the development of marine structures, because it allows to install offshore structures in safer areas, reducing possible damages, costs and eventual environmental pollution. In this context, the objective of the present work is to investigate the stability of submarine soil masses subjected to the wave cyclic loading through analytical-numerical approaches. The effect of water waves on the submerged bed is described by the propagation of a pressure wave along its surface, using the Stokes’s linear theory. Soil masses with horizontal and sloped upper surface, composed of cohesive material (clays) and granular materiais (sands) are considered in this study. In soil masses constituted of fine soil, the material strength capacity is modeled by the non-homogeneous Tresca criterion and the stability analysis is carried out in undrained condition. On the other hand, in granular beds, the strength explicitly depends on the pore pressure value, being classically described by the Coulomb criterion without cohesion. Then, the stability analysis is developed in effective stress and the pore pressure gradient acts as a volumetric force on the skeleton, characterizing the main charging mode of this material. Due to the tendency to densify when subjected to a cyclic deviatoric stress state occurs, in general, the build-up of pore pressure excess in the granular mass. Consequently, the coupling between the material behavior and the cyclic loading has fundamental importance in the calculation of the pore pressure excess generated. In order to define the seepage forces, a simplified approach based on the partition of deformations in reversible and irreversible contributions is considered, which allows to decouple the wave-induced pore pressure calculation. Applying the concepts of the limit analysis theory it is possible to formulate upper and lower boundaries of the maximum safe amplitude of the wave loading. Finally, the effects of the seabed surface steepness and of the soil layer thickness on the stability are analyzed.

Submitted by José Henrique Krähenbühl Ambiel (zeambiel@hotmail.com) on 2018-09-28T13:24:06Z No. of bitstreams: 1 Dissertação_JoséAmbiel.pdf: 8234688 bytes, checksum: e8355af378aacfa5b31cc3b2d4f77de7 (MD5) / Approved for entry into archive by Lucilene Cordeiro da Silva Messias null (lubiblio@bauru.unesp.br) on 2018-09-28T18:11:05Z (GMT) No. of bitstreams: 1 ambiel_jhk_me_bauru.pdf: 6828894 bytes, checksum: e1a58fe94084b497e2236056afe8f889 (MD5) / Made available in DSpace on 2018-09-28T18:11:05Z (GMT). No. of bitstreams: 1 ambiel_jhk_me_bauru.pdf: 6828894 bytes, checksum: e1a58fe94084b497e2236056afe8f889 (MD5) Previous issue date: 2018-07-24 / Materiais porosos constituem uma grande gama de materiais que podem ser encontrados na natureza ou em forma artificial. Rochas reservatório é um exemplo importante desse tipo de material, sendo o estudo delas a motivação principal desse trabalho. O estudo de rochas reservatório, de onde são extraídos gases e petróleo, consiste em um problema físico no qual os sistemas mecânico e hidráulico são acoplados. O acoplamento ocorre pois as deformações (no sistema mecânico) inuenciam as pressão (no sistema hidráulico), que por sua vez inuenciam as tensões (sistema mecânico). As equações governantes do sistema mecânico são mostradas e as do hidráulico deduzidas. Para a resolução do problema, o Método dos Elementos Finitos (MEF) foi utilizado para ambos os sistemas físicos, logo, as equações governantes são apresentadas em sua forma fraca e, então, aproximada pelo MEF. Numericamente, o acoplamento pode ser tratado de diferentes maneiras, seja considerando um dos sistemas de maneira bem pobre tal como fórmulas empíricas simplistas, seja considerado os sistemas de maneira individual, ou então de maneira completa. Essa última maneira de considerar um acoplamento, o acoplamento total, é formulada, programada e testada nesse trabalho. Para validar a implementação, dois problemas foram analisados: Problema de Terzaghi e Problema Mandel, ambos com solução analítica conhecidas. Os resultados obtidos numericamente comparados aos analíticos indicam que o método totalmente acoplado foi bem implementado, tanto em 2D quanto em 3D. Nesse trabalho também é mostrada a oscilação numérica que há em problemas de acoplamento hidromecânico e uma das formas de amenizá-la. / Porous materials constitute a wide range of materials that can be found in nature and arti cially. Reservoir rock is an important example of this kind of material, which is the main motivation of this work. The study of reservoir rocks, from which gases and oil are extracted, consists of a physical problem in which mechanical and hydraulic systems are coupled. The coupling occurs because the deformations (in the mechanical system) in uence the pressure (in the hydraulic system), which in turn in uence the stresses (mechanical system). The governing equations of the mechanical system are shown and those of the hydraulic system are deduced. To solve the problem, the Finite Element Method (FEM) is used for both physical systems, so the governing equations are presented in their weak form and then approximated according to the FEM. Numerically, the coupling can be handled in di erent ways, either by considering one of the systems in a very poor way by using simplistic empirical formulas, by considering the systems individually, or in a complete manner. The latter one, the fully-coupled treatment, is formulated, programmed and tested in this work. To validate the implementation, two problems has been analyzed: Terzaghi Problem and Mandel Problem, both with known analytical solutions. The comparison between the results obtained numerically and analytically indicates that the fully coupled method has been well implemented in both 2D and 3D cases. The numerical oscillation existing in hydrmechanical coupled problems is also shown and one of the ways to minimize it is presented.