Modélisation directe de la propagation d'ondes dans les milieux visco-élastiques : approche 2D par différences finies dans le domaine temporel / Amine Dhemaied ; sous la direction de Roger Guérin et Fayçal Rejiba

Date :

Editeur / Publisher : [S.l.] : [s.n.] , 2011

Type : Livre / Book

Type : Thèse / Thesis

Langue / Language : français / French

Ondes sismiques

Élastodynamique

Sédiments géologiques -- Tunis (Tunisie)

Guérin, Roger (1963-.... ; géophysicien) (Directeur de thèse / thesis advisor)

Rejiba, Fayçal (1974-.... ; enseignant-chercheur en géophysique) (Directeur de thèse / thesis advisor)

Université Pierre et Marie Curie (Paris ; 1971-2017) (Organisme de soutenance / degree-grantor)

Relation : Modélisation directe de la propagation d'ondes dans les milieux visco-élastiques : approche 2D par différences finies dans le domaine temporel / Amine Dhemaied ; sous la direction de Roger Guérin et Fayçal Rejiba / Lille : Atelier national de reproduction des thèses , 2011

Résumé / Abstract : La modélisation de la propagation d'ondes pour des géométries bidimensionnelles avec des structures complexes associée à un comportement intrinsèque non-élastique quelconque, constitue un préalable essentiel à la compréhension et à l'interprétation des données sismiques, quelle que soit leur mise en œuvre, dans les milieux naturels. La résolution des équations de l'élastodynamique par la méthode des différences finies offre une flexibilité appréciable pour la représentation de modèles complexes, mais nécessite des développements numériques très spécifiques pour 1) inclure un modèle de comportement anélastique quelconque, 2) définir des conditions absorbantes efficaces afin de simuler un domaine ouvert, 3) définir des conditions de surfaces optimales afin de permettre l'évaluation des champs en présence de forts contrastes d'impédance, particulièrement en présence d'une topographie à l'interface air-sol. Dans le cadre de cette thèse plusieurs approches algorithmiques ont été développées et testées afin de répondre aux trois impératifs précédents. La discrétisation des équations de l'élastodynamique par différences finies dans le domaine temporel a été effectuée au second ordre en temps et quatrième ordre en espace (O(2,4)). Les conditions absorbantes sont de type 'Conditional Perfectly Matched Layer' afin de maximiser leur efficacité pour un grand nombre de configurations (caractéristiques du milieu, angles d'incidence). Le maillage est proposé soit sous sa forme classique dite orthogonale, soit avec une rotation appliquée (rotated grid) pour une évaluation des champs de vitesses et de contraintes plus précise à l'interface air-sol en présence de topographie. Concernant le modèle de comportement visco-élastique, une implémentation classique de type récursive est proposée, ainsi qu'une version originale de la méthode 'Auxiliary Differential Equation' pour la généralisation à des modèles visqueux quelconques (linéaire ou non-linéaire). Les validations numériques de rigueur et les comparaisons dans divers cas canoniques pour une validation expérimentale sur des assemblages aluminium/polyméthacrylate, puis un milieu granulaire disposé dans une cuve ont permis 1) d'asseoir la validité des modèles numériques vis à vis des temps de propagation et 2) d'établir les sources d'incertitudes sur les amplitudes essentiellement dues au caractère bidimensionnel des simulations. Par ailleurs, une étude du bassin de Tunis (Tunisie) initiée par une campagne de sismique réfraction et la modélisation de formes schématiques de bassin a été abordée. Cette dernière a permis de mettre en évidence les effets de site prévisibles liés à la forme et au type de remplissage d'un bassin sédimentaire.