Imagerie sismique quantitative de la marge convergente d'Equateur-Colombie : Application des mèthodes tomographiques aux données de sismique réflexion multitrace et réfraction-réflexion grand-angle des campagnes SISTEUR et SALIERI

• Main Author: Agudelo, William
• Language: FRE
• Published: Université Pierre et Marie Curie - Paris VI 2005
• Subjects: [SDU:STU:GP] Sciences of the Universe/Earth Sciences/Geophysics; Géodynamique; marges convergents; subduction; zones de subduction; zone de subduction d'Equateur-Colombie; grandes séismes dans les zones de subduction; zone sismogénique; géophysique marine; sismique réfléxion multitrace; sismique grand-angle; mèthodes d'inversion; tomographie sismique; migration avant sommation; modélisation des amplitudes; modélisation des vitesses; chennal de subduction; décollement interplaque; Bottom Simulating Reflector; Asperités sismologiques; splay fault
• Online Access: http://tel.archives-ouvertes.fr/tel-00010415; http://tel.archives-ouvertes.fr/docs/00/04/85/43/PDF/tel-00010415.pdf; http://tel.archives-ouvertes.fr/docs/00/04/85/43/ANNEX/tel-00010415.ppt

Mon travail de thèse se propose d'étudier la structure, les propriétés physiques et les processus géodynamiques de la zone de subduction d'Equateur-Colombie grâce à l'adaptation et le développement d'outils d'imagerie sismique (inversion de formes d'ondes 'alias' tomographie en diffraction) et à leur application aux données de sismique marine multitrace (MCS) et grand-angle OBS (WA) acquises en Equateur-Colombie pendant les campagnes SISTEUR et SALIERI. Ces outils m'ont permis de réaliser une imagerie fine et quantitative à trois niveaux : l'imagerie superficielle (~ 0-3 km), l'imagerie à profondeur intermédiaire (~ 3-10 km) et l'imagerie profonde (~ 10-30 km). Dans le domaine superficiel, j'ai effectué une cartographie fine et quantitative des propriétés physiques des sédiments au voisinage du BSR (Bottom Simulating Reflector), interpreté comme la base de stabilité des hydrates de gaz. Sur le profil SIS-40 situé sur la marge sud de la Colombie, j'ai pu identifier la présence de failles qui perturbent localement le BSR. Les résultats présentés sous la forme d' une série de logs adjacents de l'image migrée en profondeur, montrent que certaines régions du BSR sont caractérisées par une augmentation de la vitesse (1470-1650 m/s), indiquant la présence d'une faible quantité d'hydrates de gaz au dessus du BSR; d'autres zones situées immédiatement sous le BSR sont caractérisées par une diminution de la vitesse (~1200 m/s), liée à la présence de gaz libres piégés sous la couche d'hydrate de gaz. A des profondeurs moyennes j'ai étudié la structure du chenal de subduction (profil SIS-72). Le chenal constitue la limite mécanique entre la plaque chevauchante et la plaque plongeante. Il est délimité à son toit par un fort réflecteur interprété comme le décollement interplaque et à sa base par le toit très réflectif de la croûte océanique en subduction. L'imagerie fine et quantitative des propriétés physiques du décollement interplaque permet de mieux comprendre le rôle de la circulation des fluides et des variations lithologiques et physiques, sur le couplage mécanique inter-plaque. En raison de la sensibilité de la méthode de tomographie en diffraction au macro-modèle de vitesse, un code de correction de ce modèle a été implémenté, afin d'obtenir des images tomographiques fiables (i.e. géométrie et amplitudes correctes). Du fait de la bande passante limitée de la source et de la longueur du dispositif d'acquisition limitée à 4.5 km, les images tomographiques ont une résolution spatiale limitée : l'image tomographique présente un déficit des petits et grands nombre d'onde (fréquences spatiales) limitant ainsi l' interprétation géologique des paramètres physiques cartographiés. Un traitement spécifique basé sur la modélisation des traces sismiques a été implémenté. L'image tomographique, traitée comme une série de traces verticales, constitue la donnée observée. L'espace des modèles est constitué par un ensemble de modèles impulsionnels et unidimensionnels de Terre construits aléatoirement. Ces modèles sont dégradés par convolution avec une estimation de l'ondelette source afin de fournir une représentation synthétique de l'image tomographique « observée ». La minimisation de la fonction coût entre les traces migrées et les traces synthétiques est effectuée dans le cadre d'une inversion globale par recuit simulé (VFSA= « Very Fast Simulated Annealing »). Le modèle moyen issu de cette procédure fournit un modèle 2D fin de vitesse, fonction de la profondeur et comparable à la limite de la résolution théorique de la source. A l'issue de ce traitement, des perturbations de vitesse positives sont mises en évidence au toit de la croûte, et d'autres négatives accompagnent certains segments du niveau du décollement. Ces dernières sont probablement associées à la présence de fluides. Le domaine plus profond a été étudié à partir des données MCS et WA dans le double but (1) d'améliorer la résolution spatiale des images sismiques du Moho et du contact interplaque en relation avec la zone sismogène, et (2) de détecter la présence d'anomalies crustales de vitesse et d'analyser leur relation avec les zones d'aspérité sismologiques. L'utilisation conjointe des données de sismique MCS et WA a été mise en oeuvre pour prolonger vers le bas les images de sismique verticale et tenter ainsi d'établir une relation entre les processus profonds et les manifestations en surface. L'application de la chaîne de traitement au profil SIS-44 a permis d'obtenir un modèle de vitesse bien contraint jusqu'à 25 km de profondeur. Ce modèle met en évidence des réflecteurs profonds (Moho et contact interplaque ) et des réflecteurs plus superficiels (splay fault), dont l'interprétation était initialement incertaine sur les images migrées en temps.