Ce travail concerne l'utilisation de la technique de Fracturation Hydraulique (FH) pour exploiter l'énergie géothermique des réservoirs profonds de roches sèches chaudes (HDR). La fracturation hydraulique est réalisée par injection de fluides géothermiques dans des réservoirs partiellement fracturés de faible perméabilité. Les fluides à haute pression sont destinés à faire évoluer les fissures et leur connectivité. Les valeurs de débit/pression auxquelles les fluides géothermiques doivent être pompés, ainsi que le calendrier de pompage pour initier la fracturation hydraulique, dépendent principalement des conditions géostatiques existantes (contraintes géostatiques, pression fluide et température initiales de l'HDR) ainsi que des propriétés des fissures de l'HDR (longueur, épaisseur, densité et distribution directionnelle initiales moyennes de fissures). Tous ces éléments, en sus de leurs effets sur la stabilité des forages, sont analysés dans cette recherche. Des modèles de fracturation, qui sont capables de suivre l'évolution des fissures dans toutes les orientations spatiales possibles, sont utilisés pour obtenir le tenseur anisotrope de perméabilité. Ces modèles sont intégrés dans un code domestique d'éléments finis qui est développé pour résoudre des problèmes aux limites thermo-poroélastiques. Pour supprimer/diminuer les oscillations qui accompagnent les solutions paraboliques et/ou hyperboliques lors de la convection forcée, plusieurs techniques de stabilisation ont dû être implémentées. / The application of the Hydraulic Fracturing (HF) technology to exploit geothermal energy from Hot Dry Rocks (HDR) reservoirs is addressed. HF is achieved by extensively pumping geothermal fluids to already existing fractured HDR reservoirs of low permeability. High fluid pressures are expected to drive cracks to evolve and connect. The newly created burgeoning hydraulic conduits should supposedly enhance the permeability of the existing HDR reservoirs. The flow rate/pressure values at which geothermal fluids should be pumped, as well as the pumping schedule to initiate HF, depend primarily on the existing geostatic conditions (geostatic stresses, initial HDR pressure and temperature) as well as on HDR fracture properties (initial mean fracture length, mean fracture aperture, density and orientational distribution of fractures). While these components, in addition to their effects on borehole stability, are scrutinized in this research, focus is on the evolution during circulation processes of the fracture properties. A fracturing model that is capable of tracking fracture evolution in all possible spatial orientations is used to obtain the time course of the anisotropic permeability tensor. This evolving property is integrated into a domestic finite element code which is developed to solve thermo–poroelastic BVPs: emphasis is laid on the efficiency of the doublet flow technique where a fluid gains thermal energy by circulating through the HDR reservoir from the injection well to the production well. The spurious oscillations in the hyperbolic solutions of the approximated finite element approach that are commensal with the phenomenon of forced heat convection are healed/mitigated through several stabilization approaches.

The mechanical and acoustic anisotropy of media is a governing factor in the behaviour of multilayered systems including such media. The understanding of the mechanisms conditioning the dynamic behaviour of multilayered systems is of paramount importance. In particular, the intrinsic anisotropy of poroelastic media presents a potential for the optimal design of systems for multifunctional performances. Indeed, these multilayered systems are bound by stiffness, thermal and acoustic performance constraints in simultaneously. A plane wave method is presented to study the influence of material orientation in the dynamic behaviour of multilayered systems composed of anisotropic poroelastic media. The method is applied to a system composed of an anisotropic open-celled melamine foam core in between two metal sheets. This particular multilayered configuration allows to shed light on phenomena intrinsic to layers composed of anisotropic poroelastic materials, such as the frequency shift of the fundamental resonance of the panel, or the compression-shear coupling effects taking place in the poroelastic core layers. The latter phenomena is of particular importance, as it is evidenced on the unconventional polarisation of waves in anisotropic poroelastic media. Finally, the method is adapted to the optimisation of multilayered systems for acoustic performance. the design variables are consequently the core material orientations with respect to the global coordinate system. The solutions to the optimisation problem are analysed in terms of dynamic behaviour, thus allowing to correlate acoustic performance of the overall structure, and the response of each individual layer. / <p>QC 20161110</p>

Au cours de la dernière décennie, la méthode d’analyse isogéométrique (AIG) a attiré l’attention des chercheurs grâce à ses capacités supérieures à la méthode standard des éléments finis (MEF). Le concept AIG utilise les mêmes fonctions de base que celles utilisées dans la conception assistée par ordinateur (CAO) pour l’approximation des champs inconnus tels que les déplacements, pression interstitielle ou la température dans la solution des éléments finis d’un problème thermo–hydro–mécanique (éventuellement couplé). Parmi les caractéristiques les plus importantes d’AIG, la régularité, le taux de convergence et surtout sa continuité intrinsèque d’ordre supérieur représentent une nette amélioration par rapport à la méthode standard des éléments finis, permettant d’obtenir des avantages computationnels significatifs en termes de précision de la solution et de efficacité.Ce travail tente d’exploiter les caractéristiques d’AIG pour la résolution numérique des problèmes hydromécaniques (HM) couplés dans les géomatériaux de second gradient de type poro–élastoplastiques partiellement saturés. D’une part, le modèle second gradient appartenant à la théorie des milieux continus avec microstructure assure l’objectivité des résultats en présence de phénomènes de localisation de la déformation en termes d’indépendance de maillage de la solution numérique, ce qui ne peut être réalisé avec des modèles constitutifs classiques qui n’implique pas l’intervention d’une longueur interne. D’autre part, la continuité C1 réalisable au moyen de fonctions de base AIG permet une implémentation directe de tels modèles constitutifs d’ordre supérieur, dans une formulation HM dérivée de l’approche de mélange classique. De plus, la régularité des fonctions de base AIG s’est révélée très efficace dans la modélisation de processus couplés caractérisés par de forts gradients hydrauliques – comme la simulation de la propagation d’un front de saturation dans une pente partiellement saturée. Dernier point, mais non des moindres, il convient de noter que, par rapport aux approches existantes basées sur les multiplicateurs de Lagrange, la méthode AIG pour résoudre les problèmes hydromécaniques (HM) couplés dans les matériaux du second gradient saturé et partiellement saturé permet une réduction considérable du nombre de degrés de libertés requis pour atteindre le même niveau de précision. Cela entraîne non seulement une augmentation significative de l’efficacité de calcul, mais permet également d’étendre la formulation du second gradient à l’analyse de problèmes réalistes en 3D, dont la solution a été présentée pour la première fois dans ce travail.La formulation poro–élastoplastique du second gradient développée dans ce travail est mise en œuvre dans le code orienté vers la recherche GeoPDEs, un code IAG–MEF open source écrit en Matlab et développé à l’Université de Pavia. Sur la base des résultats obtenus dans une large série de problèmes aux limites en 2D et 3D analysées dans ce travail, on peut conclure que la combinaison de AIG et d’élastoplasticité du second gradient représente un outil puissant pour la simulation numérique de problèmes géotechniques caractérisés par de forts couplages multiphysiques, un comportement fortement non linéaire du sol, et des gradients de déplacement et de pression interstitielle fortement localisés. / During the last decade, Isogeometric Analysis (IGA) has drawn the attention of the Finite Element community to its superior capabilities over the standard Finite Element Method (FEM). The IGA concept uses the same basis functions used in Computed Aided Design (CAD) for the approximation of the unknown fields such as displacements, pore pressure or temperature in the Finite Element solution of a (possibly coupled) thermo– hydro–mechanical problem. Among the most relevant features of IGA, its smoothness, its convergence rate and particularly its intrinsic higher–order continuity between elements represent a definite improvement over the standard FEM, which allow to obtain significant computational advantages in terms of accuracy of the solution and computa- tional efficiency.This work attempts to exploit the characteristics of IGA for the numerical solution of coupled hydro–mechanical (HM) problems in saturated and partially saturated second gradient poro–elastoplastic geomaterials. On one hand, the second gradient model belonging to the theory of continua with microstructure ensures the objectivity of the results in presence of strain localization phenomena in terms of mesh independence of the numerical solution, which cannot be achieved with classical constitutive models without an internal length scale. On the other hand, the C1–continuity achievable by means of IGA basis functions allows a straightforward implementation of such higher order constitutive models, within a HM formulation derived from the classical mixture approach. In addition, the smoothness of the IGA basis functions proved to be very efficient in the modeling of coupled processes characterized by strong hydraulic gradients – such as the simulation of the downward propagation of a saturation front in a partially saturated slope subject to rainfall infiltration. Last but not least, it is worth noting that, as compared to the existing approaches based on Lagrange multipliers, the IGA approach to the solution of coupled hydro-mechanical (HM) problems in saturated and partially saturated second gradient materials allows a dramatic reduction in the number of degrees of freedoms required to achieve the same level of accuracy. This not only results in a significant increase of the computational efficiency, but also allows to extend the complete second gradient formulation to the analysis of realistic 3D problems, the solution of which has been presented in this work for the first time.The local second gradient poro–elastoplastic formulation developed in this work is implemented in the research-oriented code GeoPDEs, a Matlab open source IGA–FEM code developed at the University of Pavia. Based on the results obtained in a large series of representative 2D and 3D initial–boundary value problems analyzed in this work, it can be concluded that the combination of IGA and the second gradient elastoplasticity represents a powerful tool for the numerical simulation of geotechnical problems characterized by strong multiphysics couplings, highly nonlinear behavior of the soil, and strongly localized displacement and pore pressure gradients.

Dans un puits pétrolier, une gaine de ciment est coulée entre la roche réservoir et le cuvelage d'acier. Cette gaine de ciment a un rôle primordial dans la stabilité et l'étanchéité des puits pétroliers. Durant les différentes étapes de la vie d'un puits, depuis le forage jusqu'à l'exploitation et finalement l'abandon du puits, le ciment est soumis à différentes sollicitations mécaniques et thermiques. Ces sollicitations peuvent endommager le ciment et dégrader ses propriétés mécaniques et sa perméabilité et par conséquent modifier sa contribution dans la stabilité et l'étanchéité du puits. La connaissance du comportement de ce ciment dans des conditions de fond de puits (hautes températures, fortes pressions) et de son évolution au cours du temps est essentielle pour l'analyse de fonctionnement du puits pendant l'exploitation et aussi pour garantir son étanchéité pour le stockage et la séquestration des gaz à effet de serre dans des réservoirs pétroliers abandonnés. C'est dans ce but qu'est développée cette étude sur la caractérisation du comportement thermo-poro-mécanique des ciments pétroliers sous fortes contraintes pour des températures entre 20 et 90°C. Les paramètres poro-mécaniques du ciment et leurs variations avec l'état de contraintes et la température ont été évalués dans un programme d'essais de compression isotrope drainés, non drainés et sans membrane. Un phénomène de dégradation des propriétés élastiques a été observé dans les essais. Les observations microscopiques ont montré que cette dégradation est provoquée par la microfissuration du matériau y compris sous chargement isotrope. La bonne compatibilité et la cohérence des paramètres poro-mécaniques obtenus montrent que le comportement de ciment peut être décrit dans le cadre de la théorie des milieux poreux. Les essais poro-mécaniques ont été couplés avec une évaluation de la perméabilité du ciment sous chargement en régime transitoire. Des essais déviatoriques ont été également réalisés pour étudier l'endommagement du matériau sous chargement déviatorique. Lorsqu'un puits pétrolier est soumis à de rapides variations de températures, les déformations et la surpression interstitielle qui en résultent peuvent porter atteinte à l'intégrité du puits. Ce phénomène de pressurisation thermique induit lors d'un chauffage en conditions non-drainées est dû à la différence entre le coefficient d'expansion thermique du fluide interstitiel et celui de la matrice solide. Une étude préliminaire, théorique et expérimentale a été menée sur un grès. Les effets des déformations mécaniques et thermiques du système de drainage lors d'un essai de chauffage non-drainé ont été analysés et une méthode simple a été développée pour apporter une correction aux valeurs de pressions interstitielles mesurées. Le comportement thermique du ciment a été étudié expérimentalement dans les essais de chauffage drainés et non-drainés et le coefficient de pressurisation thermique ainsi que les coefficients d'expansion thermique drainé et non-drainé du ciment ont été évalués. Le coefficient de pressurisation thermique a été évalué à 0.6 MPa/°C, ce qui montre l'importance des changements de contraintes effectives que peut occasionner un brusque changement de température. L'anomalie du comportement thermique du fluide interstitiel du ciment a été attribuée à son confinement dans les pores de très petite taille et cette anomalie a été étudiée à partir des résultats expérimentaux. En complément de travail sur le comportement du ciment, une étude expérimentale a été menée sur la loi de contrainte effective qui contrôle la perméabilité d'un calcaire et un modèle conceptuel basé sur la micro-texture du matériau a été proposé.

Ce travail concerne l'utilisation de la technique de Fracturation Hydraulique (FH) pour exploiter l'énergie géothermique des réservoirs profonds de roches sèches chaudes (HDR). La fracturation hydraulique est réalisée par injection de fluides géothermiques dans des réservoirs partiellement fracturés de faible perméabilité. Les fluides à haute pression sont destinés à faire évoluer les fissures et leur connectivité. Les valeurs de débit/pression auxquelles les fluides géothermiques doivent être pompés, ainsi que le calendrier de pompage pour initier la fracturation hydraulique, dépendent principalement des conditions géostatiques existantes (contraintes géostatiques, pression fluide et température initiales de l'HDR) ainsi que des propriétés des fissures de l'HDR (longueur, épaisseur, densité et distribution directionnelle initiales moyennes de fissures). Tous ces éléments, en sus de leurs effets sur la stabilité des forages, sont analysés dans cette recherche. Des modèles de fracturation, qui sont capables de suivre l'évolution des fissures dans toutes les orientations spatiales possibles, sont utilisés pour obtenir le tenseur anisotrope de perméabilité. Ces modèles sont intégrés dans un code domestique d'éléments finis qui est développé pour résoudre des problèmes aux limites thermo-poroélastiques. Pour supprimer/diminuer les oscillations qui accompagnent les solutions paraboliques et/ou hyperboliques lors de la convection forcée, plusieurs techniques de stabilisation ont dû être implémentées.

Abstract In Finland there is a food monitoring program which has found elevated dioxin and polychlorinated biphenyl concentrations in the muscle of semi-domesticated reindeer (Rangifer tarandus tarandus) calves. This led to further research on the concentrations of persistent organic pollutants in reindeer muscle, liver, and other internal organs. The research was further expanded on wild moose (Alces alces) muscle and liver. The main objective of this thesis is to increase knowledge of polychlorinated dibenzo-p-dioxin (PCDD), polychlorinated dibenzofuran (PCDF), polychlorinated biphenyl (PCB), and polybrominated diphenyl ether (PBDE) pollution levels in the Finnish terrestrial environment, and in semi-domesticated reindeer and wild moose in particular. The research gives information of exposure conditions in the reindeer’s food chain, as well as species differences and individual variation in accumulation and distribution of pollutants in reindeer and moose. Local differences between the contaminant concentrations were explored. Toxic equivalencies (TEQs) set by the World Health Organization (WHO) of PCDD/Fs and PCBs were calculated in order to assess the validity of selling reindeer and moose tissue. It was observed that there is a species-, individual-, and tissue-specific accumulation of dioxins, dioxin-like PCB, and PBDE compounds in reindeer and moose. Varying exposure conditions mainly explain the differences, although taking into account the age of an individual animal, its metabolic patterns have a role, too. Reindeer placenta and milk proved to be important factors in the transporting of compounds from hind to calf. The highest PCDD/F and PCB concentrations (as WHO-TEQs) were observed in reindeer calves in the study area in which animals are fed in natural pastures only. Despite the findings, it was concluded that it is safe to eat reindeer and moose meat since the concentrations of dioxins and dioxin-like PCBs in the muscle are relatively low, and because of the low fat content in these animals. Reindeer liver, for its part, had quite a lot of dioxin-like compounds that may compromise its safety as food, at least on a regular basis. The study shows that Finnish semi-domesticated reindeer and wild moose are good indicator species of POP contamination in a terrestrial environment, reindeer describing the situation in northern parts of the country especially. / Tiivistelmä Suomalaisessa ruokamonitorointitutkimuksessa löydettiin kohonneita dioksiinien ja polykloorattujen bifenyylien pitoisuuksia puolivillin poron vasojen lihaksista. Se johti lisätutkimuksiin pysyvien orgaanisten yhdisteiden pitoisuuksista poron lihaksessa, maksassa ja muissa sisäelimissä. Tutkimus laajeni koskemaan myös hirven (Alces alces) lihasta ja maksaa. Väitöskirjatyö lisää tietämystä polykloorattujen dibentso-p-dioksiinien (PCDD), polykloorattujen dibentsofuraanien (PCDF), polykloorattujen bifenyylien (PCB) ja polybromattujen difenyylieettereiden (PBDE) pitoisuuksista suomalaisessa maaympäristössä ja erityisesti porossa ja hirvessä. Tutkimus antaa tietoa yhdisteille altistumisesta sekä viitteitä lajien- ja yksilöiden välisistä eroavai¬suuksista yhdisteiden kertymisessä ja niiden jakaantumisessa eri kudosten välillä. Alueellista vaihtelua yhdisteiden pitoisuuksissa selvitettiin myös. Maailman terveysjärjestön (WHO) asettamia PCDD/F- ja PCB-yhdisteiden toksisuusekvivalenttiarvoja (TEQ) tarkasteltaessa pystyttiin arvioimaan kemiallista elintarvikekelpoisuutta suhteessa EU:n antamiin sallittuihin pitoisuuksiin. Työssä havaittiin laji-, yksilö-, ikä- ja kudosspesifistä dioksiinien, dioksiinienkaltaisten PCB- ja PBDE-yhdisteiden kerääntymistä porossa ja hirvessä. Tämä on todennäköisimmin seurausta vaihtelevasta altistumisesta yhdisteille, mutta myös lajien metabolisissa toiminnoissa yksilön eri ikäkausina voi olla eroavaisuuksia. Poron istukan ja maidon havaittiin olevan tärkeitä yhdisteiden kulkeutumisessa emolta sikiölle ja vasalle. Korkeimmat dioksiinien ja dioksiinien kaltaisten PCB-yhdisteiden konsentraatiot (WHO-TEQ-pitoisuuksina) havaittiin poron vasoilla tutkimusalueella, jossa eläimet olivat laiduntaneet ainoastaan luonnonlaitumilla. Huolimatta havaituista haitta-ainepitoisuuksista todettiin, että sekä luonnostaan vähärasvaisen poron että hirvenlihan syönti on turvallista suhteellisen alhaisten dioksiinien ja dioksiininkaltaisten PCB-yhdisteiden pitoisuuksien perusteella. Poron maksa puolestaan sisälsi melko korkeita dioksiinien kaltaisten yhdisteiden pitoisuuksia, mikä voi vaikuttaa sen turvalliseen käyttöön elintarvikkeena ainakin usein syötynä. Suomalainen puolivilli poro ja villi hirvi sopivat hyvin POP-kontaminoitumisen indikaattoreiksi maaympäristössä; poron erityisesti kuvaten tilannetta maan pohjoisosissa.

The aim of this diploma thesis is design of multifunctional building in the town Zábřeh in dis-trict Šumperk. The house is no basement. It is divided into two operating parts – commercial space and residential space. The building has four floors. Vertical supporting structure are formed by ceramic system Porotherm, horizontal supporting structures are formed by ceramic inserts Porotherm. The building is covered with a flat roof.

Des milieux polyphasiques aux milieux complexes

Dans cette thèse nous développons des estimations d'erreur a posteriori par équilibrage de flux pour la poro-mécanique et la poro-plasticité.En se basant sur ces estimations, nous proposons des algorithmes adaptatifs pour la résolution numérique de problèmes en mécanique des sols.Le premier chapitre traite des problèmes en poro-élasticité linéaire.Nous obtenons une borne garantie sur l'erreur en utilisant des reconstructions équilibrées et $H({rm div})$-conformes de la vitesse de Darcy et du tenseur de contraintes mécaniques.Nous appliquons cette estimation dans un algorithme adaptif pour équilibrer les composantes de l'erreur provenant de la discrétisation en espace et en temps pour des simulations en deux dimensions.La contribution principale du chapitre porte sur la reconstruction symétrique du tenseur de contraintes.Dans le deuxième chapitre nous proposons une deuxième technique de reconstruction du tenseur de contraintes dans le cadre de l'élasticité nonlinéaire.En imposant la symétrie faiblement, cette technique améliore les temps de calcul et facilite l'implémentation.Nous démontrons l'éfficacité locale et globale des estimateurs obtenus avec cette reconstruction pour une grande classe de lois en hyperélasticité.En ajoutant un estimateur de l'erreur de linéarisation, nous introduisons des critères d'arrêt adaptatifs pour le solveur de linéarisation.Le troisième chapitre est consacré à l'application industrielle des résultats obtenus. Nous appliquons un algorithme adaptatif à des problèmes poro-mécaniques en trois dimensions avec des lois de comportement mécanique élasto-plastiques. / In this Ph.D. thesis we develop equilibrated flux a posteriori error estimates for poro-mechanical and poro-plasticity problems.Based on these estimations we propose adaptive algorithms for the numerical solution of problems in soil mechanics.The first chapter deals with linear poro-elasticity problems.Using equilibrated $H({rm div})$-conforming flux reconstructions of the Darcy velocity and the mechanical stress tensor, we obtain a guaranteed upper bound on the error.We apply this estimate in an adaptive algorithm balancing the space and time discretisation error components in simulations in two space dimensions.The main contribution of this chapter is the symmetric reconstruction of the stress tensor.In the second chapter we propose another reconstruction technique for the stress tensor, while considering nonlinear elasticity problems.By imposing the symmetry of the tensor only weakly, we reduce computation time and simplify the implementation.We prove that the estimate obtained using this stress reconstuction is locally and globally efficient for a wide range of hyperelasticity problems.We add a linearization error estimator, enabling us to introduce adaptive stopping criteria for the linearization solver.The third chapter adresses the industrial application of the obtained results.We apply an adaptive algorithm to three-dimensional poro-mechanical problems involving elasto-plastic mechanical behavior laws.