[pt] A partir da amostragem de blocos de travertino cisalhados contendo materiais representativos de zona de dano e/ou núcleo de falhas geológicas, corpos-de-prova cilíndricos constituídos tanto por rocha intacta como por rocha de falha, foram submetidos a ensaios triaxiais multi-estágios com medição de permeabilidade na direção axial, com confinamentos efetivos variáveis de 2, 5, 11, 19 e 29 MPa. O presente trabalho avalia a adequação do modelo Mohr-Coulomb com endurecimento em representar o comportamento tensão-deformação de materiais de zona de falha em travertino observado nos ensaios de laboratório e determina os parâmetros que melhor representam o comportamento experimental, além de analisar a variação da permeabilidade. Para isso, é utilizado o método de elementos finitos, através do software ABAQUS, para modelar os corpos-de-prova e simular o comportamento tensão-deformação de dois dos ensaios triaxiais realizados. Ao final do estudo são determinados os parâmetros de rocha de falha que melhor ajustam aos dados experimentais. Além disso, são determinados fatores multiplicadores de permeabilidade que representam as variações de permeabilidade esperadas a partir das deformações volumétricas nesses materiais. Esse conhecimento tem importante contribuição na previsibilidade do comportamento hidromecânico de falhas geológicas, cujo equilíbrio é alterado com a explotação de reservatórios de petróleo. / [en] From sheared travertine blocks, containing both damage zone and/or gouge representative materials, cylindrical plugs were obtained containing both intact rock and fault materials, and were submitted to triaxial tests with axial permeability measurements in effective confinement pressures of 2, 5, 11, 19 e 29 MPa. This work assess the suitability of Mohr-Coulomb hardening model to represent stress-strain behaviour of travertine fault zone materials observed in laboratory measurements and determines the parameters that best fit the experimental behavior. Furthermore, a permeability variation analysis is conducted. For this, finite element method is used through ABAQUS software, to model the plugs and to simulate the stress-strain behaviour of three triaxial tests. It is possible after all, to determine rock parameters that adjust to experimental results. Besides that, permeability multipliers are determined to adjust permeability changes due to volumetric deformations in these materials. This knowledge is an important contribution to forecast hidromechanical behavior of geological faults that may have stability altered by hydrocarbon exploitation in petroleum reservoirs.

L'importance croissante de l'évaluation de la performance des structures soumis au chargement sismique souligne la nécessité d'estimer le risque de liquéfaction. Dans ce scénario extrême de la liquéfaction du sol, des conséquences dévastatrices sont observées, par exemple des tassements excessifs et des instabilités de pentes. Dans le cadre de cette thèse, la réponse dynamique et l'interaction d'un système ouvrage en terre-fondation sont étudiées, afin de déterminer quantitativement le mécanisme de ruine dû à la liquéfaction du sol de la fondation. Par ailleurs, les chargements sismiques peuvent induire dans les ouvrages en terre un mode de rupture générant des bandes de cisaillement. Une étude de sensibilité aux maillages a donc été engagée pour quantifier la dépendance des résultats de l'analyse dynamique. Par conséquent, l'utilisation d'une méthode de régularisation est évaluée au cours des analyses dynamiques. Le logiciel open-source Code_Aster, basé sur la méthode des Eléments Finis et développé par EDF R&D, est utilisé pour les simulations numériques, tandis que le comportement du sol est représenté par le modèle de comportement de l'ECP, développé à CentraleSupélec. En premier lieu, un modèle simplifié de propagation 1D des ondes SH dans une colonne de sol avec comportement hydromécanique couplé non linéaire a été simulé. L'effet des caractéristiques du signal sismique et de la perméabilité du sol sur la liquéfaction est évalué. Le signal sismique d'entrée est un élément important pour l'apparition de la liquéfaction, puisque la durée du choc principal peut conduire à de fortes non linéarités et à un état de liquéfaction étendu. En outre, quand une variation de perméabilité en fonction de l'état de liquéfaction est considérée, des changements significatifs sont observés pendant la phase de dissipation de la surpression interstitielle de l'eau et au comportement du matériau. En revanche, ces changements ne suivent pas une tendance unique. Puis, l'effet d'une méthode de régularisation avec cinématique enrichie, appelée premier gradient de dilatation, sur la propagation des ondes SH est étudié au travers d'une solution analytique. Des problèmes à la réponse dynamique du sol sont observés et discutés quand cette méthode de régularisation est appliquée. Ensuite, un modèle 2D d'un déblai est simulé et sa réponse dynamique est évaluée en conditions sèches, complètement drainées et hydromécanique couplées. Deux critères sont utilisés pour définir le début de la rupture de la structure. Le travail du second ordre est utilisé pour décrire l'instabilité locale à des instants spécifiques du mouvement sismique, tandis que l'estimation d'un facteur de sécurité locale est proposée prenant en compte la résistance résiduelle du sol. En ce qui concerne le mode de ruine, l'effet de la surpression interstitielle de l'eau est de grande importance, puisqu'un déblai stable en conditions sèches et complètement drainées, devient instable lors de l'analyse couplée à cause de la liquéfaction de la fondation. Enfin, un système digue-fondation est simulé et l'influence de la perméabilité du sol, la profondeur de la couche liquéfiable, ainsi que, les caractéristiques du séisme sur la ruine induite par la liquéfaction du sol est évaluée. Pour ce modèle de digue, le niveau de dommages est fortement lié à la fois à l'apparition de la liquéfaction dans la fondation et la dissipation de la surpression d'eau. Une surface d'effondrement circulaire est générée à l'intérieur de la couche du sol liquéfié et se propage vers la crête dans les deux côtés de la digue. Pourtant, lorsque la couche liquéfiée est située en profondeur, la digue n'est pas affectée par la liquéfaction de la fondation pour ce cas particulier de chargement. Ce travail de recherche se concentre sur une étude de cas de référence pour l'évaluation sismique des ouvrages en terre soumis à un séisme et fournit des méthodes et outils de calculs numériques performants accessibles aux ingénieurs. / The increasing importance of performance-based earthquake engineering analysis points out the necessity to assess quantitatively the risk of liquefaction. In this extreme scenario of soil liquefaction, devastating consequences are observed, e.g. excessive settlements, lateral spreading and slope instability. The present PhD thesis discusses the global dynamic response and interaction of an earth structure-foundation system, so as to determine quantitatively the collapse mechanism due to foundation’s soil liquefaction. As shear band generation is a potential earthquake-induced failure mode in such structures, the FE mesh dependency of results of dynamic analyses is thoroughly investigated and an existing regularization method is evaluated. The open-source FE software developed by EDF R&D, called Code_Aster is used for the numerical simulations, while soil behavior is represented by the ECP constitutive model, developed at CentraleSupélec. Starting from a simplified model of 1D SH wave propagation in a soil column with coupled hydromechanical nonlinear behavior, the effect of seismic hazard and soil’s permeability on liquefaction is assessed. Input ground motion is a key component for soil liquefaction apparition, as long duration of mainshock can lead to important nonlinearity and extended soil liquefaction. Moreover, when a variation of permeability as function of liquefaction state is considered, changes in the dissipation phase of excess pore water pressure and material behavior are observed, which do not follow a single trend. The effect of a regularization method with enhanced kinematics approach, called first gradient of dilation model, on 1D SH wave propagation is studied through an analytical solution. Deficiencies of the use of this regularization method are observed and discussed, e.g. spurious waves apparition in the soil’s seismic response. Next, a 2D embankment-type model is simulated and its dynamic response is evaluated in dry, fully drained and coupled hydromechanical conditions. Two criteria are used to define the onset of the structure’s collapse. The second order work is used to describe the local instability at specific instants of the ground motion, while the estimation of a local safety factor is proposed by calculating soil’s residual strength. Concerning the failure mode, the effect of excess pore water pressure is of great importance, as an otherwise stable structure-foundation system in dry and fully drained conditions becomes unstable during coupled analysis. Finally, a levee- foundation system is simulated and the influence of soil’s permeability, depth of the liquefiable layer, as well as, characteristics of input ground motion on the liquefaction-induced failure is evaluated. For the current levee model, its induced damage level (i.e. settlements and deformations) is strongly related to both liquefaction apparition and dissipation of excess pore water pressure on the foundation. A circular collapse surface is generated inside the liquefied region and extends towards the crest in both sides of the levee. Even so, when the liquefied layer is situated in depth, no effect on the levee response is found. This research work can be considered as a reference case study for seismic assessment of embankment-type structures subjected to earthquake and provides a high-performance computational framework accessible to engineers.

Lors du dimensionnement des ouvrages en terre : remblais, digues, on observe que la plupart des matériaux sont compactés à l’optimum Proctor ou coté humide. En général, ce compactage implique que le sol se trouve dans un état où le degré de saturation est très élevé. Cruz et al (1985) ont montré qu'à un degré de saturation élevé (supérieur à 85%, voire 90% dans le cas de certains sols), la phase liquide est continue alors que l’air présent sous forme de bulles est occlus ; ce qui rend le comportement du sol complexe. L’élaboration d’un modèle de comportement pour ce type de sols nécessite une compréhension approfondie des phénomènes physico-mécaniques intervenant au sein de l’air occlus, de l'eau liquide contenant de l'air dissous et du squelette solide. Dans ce sens, un nouveau modèle hydromécanique a été développé. Ce modèle prend en compte le comportement physico-mécanique et la cinématique propre de chacun des constituants du milieu polyphasé (eau liquide, air dissous, air sous forme gazeuse et matrice solide). En particulier, dans ce modèle, nous tenons compte de la tension de surface, de la migration des phases gazeuse et liquide qui ont des impacts importants sur le comportement mécanique des sols. Le développement du modèle conduit à un système d’équations aux dérivées partielles fortement non linéaire qui peut être résolu numériquement en utilisant la méthode des éléments finis. Ce nouveau modèle a été implémenté dans un code de calcul écrit en C++ « Hydromech », développé à l'origine par Pereira (2005), qui permet de simuler les essais oedométriques suivant différents trajets de chargement hydromécanique. En particulier, ce code de calcul permet de simuler de façon cohérente la transition entre différents régimes de saturation, aussi bien dans l'espace (translation progressive d'une frontière entre deux régimes voisins) que dans le temps (passage d'un régime à l'autre en un point donné) ; ce qui constitue un problème de modélisation difficile. Les études numériques réalisées montrent que ce modèle donne des résultats cohérents et mettent en évidence sa capacité à simuler avec précision le comportement hydromécanique des sols quasi-saturés comportant de l'air occlus. / The behaviour of quasi-saturated materials is an important factor to be considered when designing cuttings and embankments in which earthwork materials are compacted to the optimum proctor density. Typically, soil compaction is performed at the optimum Proctor or on the wet side of the optimum, which means that the soil is in a highly saturated state. Cruz et al (1985) have shown that at a high degree of saturation (greater than 85% or even 90% in the case of certain soils), the liquid phase is continuous whereas the gas phase in the form of entrapped air bubbles is discontinuous. It is the presence of the entrapped air bubbles which makes the soil behaviour complex. The construction of a theoretical model for this type of soils requires the consideration of various physical-mechanical phenomena and their couplings occurring within the tri-phasic medium consisting of the solid grains, liquid water containing dissolved air and the entrapped air bubbles. In this sense, a new hydromechanical model has been developed that takes into account the physical-mechanical interactions between different phases as well as the kinematics of each constituent (liquid water, dissolved air, gaseous air and solid grains). In particular, the model accounts for the interfacial tension, migration of gaseous and liquid phases, which have important impacts on the mechanical behaviour. The development leads to a system of highly non-linear partial differential equations which can be solved numerically using the finite element method. This new model has been implemented in a numerical code “Hydromech” written in C++, developed originally by Pereira (2005) that has been used to simulate oedometer tests with different hydromechanical loading paths. In particular, this code allows to simulate consistently the transition across different regimes of saturation, both with respect to space (progressive translation of a boundary between two neighbouring regimes) and to time (transition of one regime to another at a fixed material point); which constituted a difficult modelling problem at the start. Numerical studies carried out show that this model gives consistent results providing a clear demonstration of its ability to simulate with precision the hydro-mechanical behaviour of quasi-saturated soils containing entrapped air.

[pt] O desenvolvimento e exploração das formações shale gas a nível mundial é relativamente recente, tendo seu início nos Estados Unidos no final da década 1990. A partir dos resultados obtidos com a aplicação do procedimento de fraturamento hidráulico, como método de estimulação, o estudo e avaliação de outros prospectos shale gas em outras regiões do mundo foi estimulado. No entanto, a análise, estudo e caracterização deste tipo de reservatórios são dificultados uma vez que devem ser levados em conta múltiplos fatores relacionados a geologia, mineralogia, petrofísica, geoquímica, entre outros. O fraturamento hidráulico é um processo hidromecânico acoplado com alto grau de dificuldade especialmente em reservatórios shale gas, onde existe a presença de fraturas. Um estudo numérico é conduzido neste trabalho com o objetivo de investigar o comportamento hidromecânico de uma fratura natural durante a injeção de fluido. O software UDEC (Universal Distinct Element Code) baseado no método dos elementos discretos foi empregado no desenvolvimento da modelagem numérica. O UDEC tem a capacidade de modelar o comportamento hidromecânico de uma fratura incluindo os fenômenos de abertura/fechamento das fraturas, cisalhamento e dilatação sob condições de contato ou separação. Nesta investigação numérica foram criados modelos de reservatório com fraturas naturais aleatórias por meio de uma distribuição voronoi visando representar a complexidade da rede de fraturas naturais presentes nas formações shale gas. Os resultados do estudo numérico mostram que o comportamento hidromecânico do sistema é fortemente dependente da variação de parâmetros como tensões in situ, viscosidade do fluido de fraturamento e taxa de injeção de fluido. Portanto, estes resultados proporcionam um melhor entendimento dos mecanismos de fraturamento e resposta da pressão de um tratamento de fraturamento hidráulico em um reservatório de gás não-convencional naturalmente fraturado. / [en] The development and deployment of shale gas formations around the world are relatively recent, starting in the United Stated in the late 1990. From the results obtained with the application of hydraulic fracturing as a method of stimulation, the study and evaluation of other prospects of shale gas in others places in the world was encouraged. However, the analysis, study and characterization of this type of reservoirs are difficult, because it must be taken into account several factors such as geology, mineralogy, petrophysics, geochemistry among others. Hydraulic fracturing is a complicated hydro-mechanical coupled process, with high difficulty degree especially in shale gas reservoir, where natural fractures exist. A numerical study is conducted to investigate the hydromechanical behavior of a natural fracture during fluid injection. UDEC (Universal Distinct Element Code) software based on discrete elements method was employed to numerical modeling development. UDEC has the ability to model the hydro-mechanical behavior of a fracture including phenomena like fracture enlargement, closure, slippage, and dilation under contact or separation condition. In this numerical investigation, numerical reservoir models, with random natural fractures through a distribution voronoi were created aiming to represent the network complexity of natural fractures present in shale gas formations. The numerical study results show that the hydromechanical system behavior is strongly dependent on the parameters variation such as in situ stress, fluid fracturing viscosity and fluid injection rate. Therefore, these results provide a better understanding of fracturing mechanisms and pressure response of a hydraulic fracturing treatment in a non- conventional naturally fractured reservoir.