Simulation numérique de la propagation d'ondes en milieu géologique complexe

Simulation numérique de la propagation d'ondes en milieu géologique complexe PDF Author: Elise Delavaud
Publisher:
ISBN:
Category :
Languages : fr
Pages : 174

Get Book Here

Book Description
Ce travail de thèse est consacré au développement d'un outil numérique capable de modéliser la propagation 3D d'ondes sismiques dans les milieux géologiques complexes caractérisés par des effets de site. La méthode des éléments spectraux (SEM) est particulièrement adaptée à la problématique de la réponse sismique dans de tels milieux: estimation précise des ondes de surface, capacité à prendre en compte des géométries compliquées (topographie des surfaces libres, interfaces) et possibilité d'ajuster la résolution des longueurs d'ondes sismiques dans des conditions de milieux hétérogènes par des maillages non-structurés. Des extensions numériques de la SEM liées à la réponse de bassin ont été développées: Perfectly Matched Layers filtrantes et introduction efficace d'un champ incident. La principale limitation actuelle des SEM reste cependant liée au manque de flexibilité des maillages hexaédriques malgré l'utilisation de ces maillages non-structurés dont la génération constitue une tâche complexe. Ce travail de thèse souligne ainsi les difficultés associées à la prise en compte de telles structures et décrit le processus complet qu'implique la modélisation de la réponse sismique dans des milieux géologiques complexes, de l'élaboration du modèle physique et numérique à l'analyse des résultats.Je présente des simulations 3D de la réponse sismique dans la région de Caracas (Venezuela), incluant la géométrie 3D du bassin et la chaîne de montagnes Ávila qui borde la ville. L'importance des effets liés à ces deux types de structures est évaluée, par différents scénarios d'excitation par une onde plane. Ces effets géométriques 3D révèlent des phénomènes complexes d'amplification associés à la réflexion et à la focalisation d'ondes se propageant dans toutes les directions, ainsi qu'à la génération d'ondes de surface au niveau des bords et des zones de rétrécissement du remplissage sédimentaire. Une comparaison avec des simulations 2D montre l'intérêt d'une modélisation 3D, de nettes différences en terme de temps de résidence et de niveau d'énergie et d'amplification ayant été constatées

Modélisation à base topologique de couches géologiques

Modélisation à base topologique de couches géologiques PDF Author: Allan Fousse
Publisher:
ISBN:
Category :
Languages : fr
Pages : 185

Get Book Here

Book Description
LA SIMULATION DE PROPAGATION D'ONDES DANS DES COUCHES GEOLOGIQUES POUVANT PRESENTER UNE TOPOGRAPHIE DE SURFACE ET DES INTERFACES COURBES, PASSE PAR LA MODELISATION DU MILIEU ET DE SON MAILLAGE. IL EST NECESSAIRE DE FOURNIR AU CODE NUMERIQUE DEVANT ASSURER LA PROPAGATION D'ONDES, UN TABLEAU TRIDIMENSIONNEL OU SONT INDEXEES LES REFERENCES DE CHAQUE POINT AINSI QUE LA NATURE DU MILIEU AUQUEL IL APPARTIENT. AINSI, LA METHODE DES DIFFERENCES FINIES NECESSITE LA CONNAISSANCE DES VALEURS CARACTERISTIQUES ET CALCULEES DU SOL SUR UN RESEAU STRUCTURE UNIFORME, C'EST-A-DIRE CARTESIEN. GRACE AUX NOUVEAUX DEVELOPPEMENTS SUR CETTE METHODE, CE RESEAU PEUT ETRE ENVISAGE DEFORME POUR SUIVRE LES DIFFERENTES INTERFACES. ON TRAVAILLE ALORS DANS UN ESPACE DE CALCUL OBTENU PAR UN CHANGEMENT DE COORDONNEES POUR PASSER D'UN SYSTEME CARTESIEN A DES COORDONNEES CURVILIGNES. L'AVANTAGE MAJEUR DE CETTE COMPLICATION DU CALCUL RESIDE DANS L'ELIMINATION DES DIFFRACTIONS PARASITES SUR LES INTERFACES COURBES MODELISEES EN MARCHES D'ESCALIER SUR UN MAILLAGE CARTESIEN. DE PLUS, LA PRISE EN COMPTE D'UNE TOPOGRAPHIE DE SURFACE N'EST FACILEMENT ACCESSIBLE QUE DANS UN SYSTEME DE COORDONNEES GENERALISEES. NOUS PROPOSONS UN ENSEMBLE DE METHODE ET D'OUTILS POUR SIMULER LA PROPAGATIONS D'ONDES DANS UN MILIEU TRIDIMENSIONNEL COMPLEXE. PLUS PARTICULIEREMENT NOUS TRAITONS LES TROIS PHASES PREALABLES A LA SIMULATION : LA MODELISATION 2-D OU 3-D DU DOMAINE OU DOIT S'EFFECTUER LA PROPAGATION. POUR CECI, IL S'AGIT D'UTILISER LES FONCTIONNALITES GENERIQUES DU LOGICIEL DE MODELISATION GEOMETRIQUE TOPOFIL, BASE SUR LES CARTES GENERALISEES, ET AUSSI DE NOUVELLES OPERATIONS REPONDANT AUX BESOINS SPECIFIQUES DE CETTE APPLICATION. - LA GENERATION D'UN MAILLAGE HEXAEDRIQUE REGULIER

RESOLUTION ET ETUDE NUMERIQUES DE QUELQUES PROBLEMES DE PROPAGATION D'ONDES ACOUSTIQUES EN GEOPHYSIQUE

RESOLUTION ET ETUDE NUMERIQUES DE QUELQUES PROBLEMES DE PROPAGATION D'ONDES ACOUSTIQUES EN GEOPHYSIQUE PDF Author: Quang Huy Trâǹ
Publisher:
ISBN:
Category :
Languages : fr
Pages : 160

Get Book Here

Book Description
CETTE THESE EST COMPOSEE DE QUATRE TRAVAUX DIFFERENTS QUI SE RAPPORTENT TOUS A UN MEME SUJET: LA SIMULATION NUMERIQUE DE LA PROPAGATION D'ONDE EN GEOPHYSIQUE. 1. PROPAGATION DES ONDES ACOUSTIQUES EN MILIEU ALEATOIRE: ETUDES NUMERIQUES EN 1-D ET 3-D. IL S'AGIT D'ETABLIR EXPERIMENTALEMENT LES PROPRIETES DE CETTE PROPAGATION, SACHANT QUE DE NOMBREUX TRAVAUX THEORIQUES, BASES SUR L'HYPOTHESE DE LINEARISATION, EXISTENT DEJA. 2. OPTIMISATION SOUS CONTRAINTE EN TOMOGRAPHIE DE REFLEXION: LA METHODE DU LAGRANGIEN AUGMENTE. ON RELATE ICI LA RECHERCHE D'UNE METHODE D'OPTIMISATION EFFICACE POUR RESOUDRE UN PROBLEME INVERSE CONCRET MAIS DIFFICILE. 3. DECOMPOSITION DE DOMAINES APPLIQUEE A L'EQUATION DES ONDES ACOUSTIQUES. ON PRESENTE UNE NOUVELLE METHODE DE SIMULATION QUI PERMET D'ADAPTER LE PAS DU MAILLAGE AUX CARACTERISTIQUES DE CHAQUE SOUS-DOMAINE. UNE ANALYSE DE CONVERGENCE EST FOURNIE. 4. ANALYSE DE DISPERSION ET DU COUT DES SCHEMAS AUX DIFFERENCES FINIES EN 3-D. A UNE PRECISION IMPOSEE SUR L'ERREUR DE GROUPE RELATIVE ET UNE TAILLE DONNEE DU MODELE GEOLOGIQUE, ON COMPARE LES SCHEMAS 2-4 ET 2-8 DU POINT DE VUE DU NOMBRE D'OPERATIONS FLOTTANTES NECESSAIRES

Modélisation numérique de la propagation des ondes dans les milieux poreux anélastiques

Modélisation numérique de la propagation des ondes dans les milieux poreux anélastiques PDF Author: Hormoz Modaressi
Publisher:
ISBN:
Category :
Languages : fr
Pages : 193

Get Book Here

Book Description
Influence des facteurs géomorphologiques et géologiques sur les réponses sismiques des sites constitués de sols dits déformables. Formulation physique des milieux biphasiques saturés et modélisation du comportement élastoplastique des sols. Problèmes liés à la dispersion et à la réflexion numérique des ondes dues à la discrétisation des équations différentielles en espace et en temps. Test de validation et quelques exemples complexes de vallée alluviale constituée de matériaux anélastiques et biphasiques

Contribution à la modélisation numérique de la propagation des ondes sismiques sur architectures multicœurs et hiérarchiques

Contribution à la modélisation numérique de la propagation des ondes sismiques sur architectures multicœurs et hiérarchiques PDF Author: Fabrice Dupros
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0

Get Book Here

Book Description
En termes de prévention du risque associé aux séismes, la prédiction quantitative des phénomènes de propagation et d'amplification des ondes sismiques dans des structures géologiques complexes devient essentielle. Dans ce domaine, la simulation numérique est prépondérante et l'exploitation efficace des techniques de calcul haute performance permet d'envisager les modélisations à grande échelle nécessaires dans le domaine du risque sismique.Plusieurs évolutions récentes au niveau de l'architecture des machines parallèles nécessitent l'adaptation des algorithmes classiques utilisées pour la modélisation sismique. En effet, l'augmentation de la puissance des processeurs se traduit maintenant principalement par un nombre croissant de cœurs de calcul et les puces multicœurs sont maintenant à la base de la majorité des architectures multiprocesseurs. Ce changement correspond également à une plus grande complexité au niveau de l'organisation physique de la mémoire qui s'articule généralement autour d'une architecture NUMA (Non Uniform Memory Access pour accès mémoire non uniforme) de profondeur importante.Les contributions de cette thèse se situent à la fois au niveau algorithmique et numérique mais abordent également l'articulation avec les supports d'exécution optimisés pour les architectures multicœurs. Les solutions retenues sont validées à grande échelle en considérant deux exemples de modélisation sismique. Le premier cas se situe dans la préfecture de Niigata-Chuetsu au Japon (événement du 16 juillet 2007) et repose sur la méthode des différences finies. Le deuxième exemple met en œuvre la méthode des éléments finis. Un séisme hypothétique dans la région de Nice est modélisé en tenant compte du comportement non linéaire du sol.

Modélisation de la propagation des ondes sismiques dans les milieux viscoélastiques

Modélisation de la propagation des ondes sismiques dans les milieux viscoélastiques PDF Author: Fateh Bouchaala
Publisher:
ISBN:
Category :
Languages : fr
Pages :

Get Book Here

Book Description
La modélisation de la propagation des ondes dans les milieux complexes 3d, inhomogènes, isotropes et viscoélastiques (milieux I.I.V.), est un sujet intéressant et largement d’actualité. cette modélisation aide pour la connaître la géométrie et les propriétés physiques des milieux géologiques. La viscosité du milieu de propagation est modélisée par la prise en compte de vitesses de propagation complexes. Le phénomène de propagation est modélisé par la théorie des rais avec la prise en compte de l’hypothèse de faible viscosité (hypothèse non restrictive dans le domaine de la géophysique). Lles expressions analytiques des coefficients de réflexion transmission entre deux milieux viscoélastiques ont été établies avec une étude numérique de l’effet de la faible viscosité sur ces coefficients a été menée. Le code anray est, a l’origine, un code de simulation de la propagation des ondes dans les milieux inhomogènes 3d purement élastiques. son développement repose sur la théorie des rais. Le code anray a ensuite été modifie (anrayq), afin de l’adapter a la simulation de la propagation des ondes dans les milieux I.I.V. les résultats numériques du code anrayq, ont été confrontés aux résultats numériques issus du code 5kb, code de simulation base sur la méthode de réflectivité. Enfin, des données sismiques 3d, acquises lors de la campagne hydratech (Ifremer 2002), dans une région caractérisée par la présence des hydrates de gaz et du gaz libre, ont été simulées a l’aide du code anrayq. le profil de l’atténuation viscoélastique a été aussi déterminé a partir de ces données sismiques, afin de localiser les hydrates de gaz et le gaz libre dans les couches sédimentaires.

Simulation numérique de la propagation d'ondes sismiques dans les milieux stratifiés à deux ou trois dimensions

Simulation numérique de la propagation d'ondes sismiques dans les milieux stratifiés à deux ou trois dimensions PDF Author: Laurence Martineau-Nicoletis
Publisher:
ISBN:
Category : Seismic waves
Languages : fr
Pages : 288

Get Book Here

Book Description


Correction et simplification de modèles géologiques par frontières

Correction et simplification de modèles géologiques par frontières PDF Author: Pierre Anquez
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0

Get Book Here

Book Description
Les modèles géologiques numériques 2D et 3D permettent de comprendre l'organisation spatiale des roches du sous-sol. Ils sont également conçus pour réaliser des simulations numériques afin d'étudier ou de prédire le comportement physique du sous-sol. Pour résoudre les équations qui gouvernent les phénomènes physiques, les structures internes des modèles géologiques peuvent être discrétisées spatialement à l'aide de maillages. Cependant, la qualité des maillages peut être considérablement altérée à cause de l'inadéquation entre, d'une part, la géométrie et la connectivité des objets géologiques à représenter et, d'autre part, les contraintes requises sur le nombre, la forme et la taille des éléments des maillages. Dans ce cas, il est souhaitable de modifier un modèle géologique afin de pouvoir générer des maillages de bonne qualité permettant la réalisation de simulations physiques fidèles en un temps raisonnable. Dans cette thèse, j'ai développé des stratégies de réparation et de simplification de modèles géologiques 2D dans le but de faciliter la génération de maillages et la simulation de processus physiques sur ces modèles. Je propose des outils permettant de détecter les éléments des modèles qui ne respectent pas le niveau de détail et les prérequis de validité spécifiés. Je présente une méthode pour réparer et simplifier des coupes géologiques de manière locale, limitant ainsi l'extension des modifications. Cette méthode fait appel à des opérations d'édition de la géométrie et de la connectivité des entités constitutives des modèles géologiques. Deux stratégies sont ainsi explorées : modifications géométriques (élargissements locaux de l'épaisseur des couches) et modifications topologiques (suppressions de petites composantes et fusions locales de couches fines). Ces opérations d'édition produisent un modèle sur lequel il est possible de générer un maillage et de réaliser des simulations numériques plus rapidement. Cependant, la simplification des modèles géologiques conduit inévitablement à la modification des résultats des simulations numériques. Afin de comparer les avantages et les inconvénients des simplifications de modèles sur la réalisation de simulations physiques, je présente trois exemples d'application de cette méthode : (1) la simulation de la propagation d'ondes sismiques sur une coupe au sein du bassin houiller lorrain, (2) l'évaluation des effets de site liés à l'amplification des ondes sismiques dans le bassin de la basse vallée du Var, et (3) la simulation d'écoulements fluides dans un milieu poreux fracturé. Je montre ainsi (1) qu'il est possible d'utiliser les paramètres physiques des simulations, la résolution sismique par exemple, pour contraindre la magnitude des simplifications et limiter leur impact sur les simulations numériques, (2) que ma méthode de simplification de modèles permet de réduire drastiquement le temps de calcul de simulations numériques (jusqu'à un facteur 55 sur une coupe 2D dans le cas de l'étude des effets de site) tout en conservant des réponses physiques équivalentes, et (3) que les résultats de simulations numériques peuvent être modifiés en fonction de la stratégie de simplification employée (en particulier, la modification de la connectivité d'un réseau de fractures peut modifier les écoulements fluides et ainsi surestimer ou sous-estimer la quantité des ressources produites).

SIMULATION NUMERIQUE DE LA PROPAGATION D'ONDES SISMIQUES EN MILIEU FRACTURE

SIMULATION NUMERIQUE DE LA PROPAGATION D'ONDES SISMIQUES EN MILIEU FRACTURE PDF Author: BRUNO.. RIOLLET
Publisher:
ISBN:
Category :
Languages : fr
Pages : 103

Get Book Here

Book Description
LA PRODUCTION D'UN RESERVOIR PETROLIER PEUT ETRE PLUS OU MOINS AFFECTEE LORSQUE DES PUITS SONT FORES DANS DES ROCHES FRACTUREES. LA DENSITE DE FRACTURATION ET L'ALIGNEMENT DES FRACTURES CONSTITUENT ALORS DES PARAMETRES DONT LA DETERMINATION EST IMPORTANTE POUR MAXIMISER CETTE PRODUCTION. CE MEMOIRE EST DEDIE A LA CARACTERISATION D'UN ENSEMBLE DE FRACTURES EN S'APPUYANT SUR LA MODELISATION NUMERIQUE DE LA DIFFUSION DES ONDES SISMIQUES. LA PROPAGATION DES ONDES SISMIQUES EST ICI SIMULEE A L'AIDE DE LA METHODE DES DIFFERENCES FINIES D'ORDRE ELEVE POUR UN MILIEU BIDIMENSIONNEL. CETTE METHODE EST CAPABLE DE CALCULER PRECISEMENT LE CHAMP DES DEPLACEMENTS DES ONDES SISMIQUES LORSQUE CELLES-CI SE PROPAGENT DANS DES MODELES DE VITESSE TRES FLUCTUANTS TELS QUE LES MILIEUX FRACTURES. D'UNE PART, NOUS MONTRONS, QU'UNE DISTRIBUTION ALEATOIRE DE FRACTURES PARALLELES, ET PETITES PAR RAPPORT A LA LONGUEUR D'ONDE, INDUIT UNE ANISOTROPIE DES VITESSES DES ONDES SISMIQUES EN ACCORD AVEC LA THEORIE DE MILIEU EQUIVALENT D'HUDSON (1981). PAR AILLEURS, CE MODELE THEORIQUE FOURNIT UNE BONNE APPROXIMATION DES MILIEUX FRACTURES, NOTAMMENT POUR DES DENSITES DE FRACTURATION ELEVEES. D'AUTRE PART, LORSQUE TOUTES LES FISSURES NE SONT PLUS PETITES DEVANT LA LONGUEUR D'ONDE, TEL QU'UN ENSEMBLE DE FRACTURES DE LONGUEURS FRACTALES, L'ATTENUATION APPARENTE DES ONDES DEVIENT UN PHENOMENE DE PREMIER ORDRE. LA DEPENDANCE FREQUENTIELLE DU FACTEUR D'ATTENUATION Q##1 PARAIT ALORS ETRE CONTROLEE DAVANTAGE PAR LA DIRECTION PREFERENTIELLE DES FRACTURES QUE PAR LE TYPE DE DISTRIBUTION REGISSANT LA LONGUEUR DES FRACTURES. LORSQUE LES FRACTURES NE SONT PLUS PARALLELES MAIS REPARTIES EN DEUX ENSEMBLES QUASI-ORTHOGONAUX, ON OBSERVE ALORS UN TRANSFERT PLUS IMPORTANT DE L'ENERGIE SISMIQUE DANS LES CODAS. CES DERNIERES OBSERVATIONS NOUS FONT PENSER QUE LES FRACTURES, QUELQUE SOIT LEUR NATURE, JOUENT UN ROLE DETERMINANT DANS LE CARACTERE FRACTAL DE LA CROUTE.

Integration of 3D Geological and Numerical Models Based on Tetrahedral Meshes for Hydrogeological Simulations in Fractured Porous Media

Integration of 3D Geological and Numerical Models Based on Tetrahedral Meshes for Hydrogeological Simulations in Fractured Porous Media PDF Author: Daniela Blessent
Publisher:
ISBN:
Category : Hydrogeology
Languages : en
Pages : 0

Get Book Here

Book Description
Une nouvelle approche de modélisation des milieux géologiques fracturés représentés par un modèle conceptuel de fractures discrètes et déterministes est présentée dans cette thèse. L'objectif principal de l'étude est de reproduire l'hétérogénéité et la complexité des milieux poreux fracturés dans un modèle géométrique tridimensionnel afin d'effectuer des simulations numériques dans le but d'améliorer les capacités de modélisation en hydrogéologie. Ceci est réalisé à travers le couplage entre une plateforme de modélisation géologique (GOCAD) et un code numérique (HydroGeoSphere). Les principaux défis à relever sont: la représentation géométrique du réseau de fractures, la sélection d'un maillage approprié pour la discrétisation spatiale du domaine de simulation et l'adaptation du code numérique à ce maillage. La nouvelle approche est basée sur une première phase de modélisation géologique 3D, suivie par la génération d'un maillage tétraédrique 3D et par la simulation numérique de l'écoulement souterrain en conditions saturées et du transport de solutés. En général, le maillage tétraédrique s'avère plus adéquat que les maillages de blocs ou de prismes pour discrétiser les geometries complexes telles que les milieux fracturés. De plus, une définition alternative du maillage "dual", qui est essentiel pour appliquer la méthode numérique élément finis - volume de contrôle utilisée par HydroGeoSphere, est analysée et intégrée dans le code numérique. Le code numérique proposé est d'abord vérifié par l'intermédiaire de simples scénarios de simulation dont les solutions, analytiques et numériques, sont déjà connues. La complexité des simulations est augmentée de façon graduelle. L'approche de modélisation est finalement appliquée au site Olkiluoto (Finlande) où un laboratoire de recherche souterrain est en construction afin d'évaluer la faisabilité du stockage géologique profond de déchets nucléaires à haute activité. Les techniques de modélisation géologique mises au point permettent de modéliser facilement la géométrie des fractures identifiées à travers la caractérisation géologique in situ. De plus, le modèle numérique s'avère adéquat pour la simulation de l'écoulement et du transport de solutés dans ce site complexe. Ce travail de recherche présente une contribution au développement des techniques de modélisation hydrogéologique des milieux fracturés.