Imagerie Haute Résolution Des Structures Lithosphériques Par Inversion De Formes D'ondes P Télésismiques Courte Période

Yi Wang

2017-01-01

Abstract:La tomographie sismique permet d'imager l'interieur de la Terre a partir de l'observation des ondes sismiques faite a la surface. L'inversion de forme d'ondes completes est une methode tomographique qui permet d'imager les structures lithospheriques de petite echelle. Cette approche demande des methodes numeriques efficaces et precises pour resoudre l'equation des ondes dans des milieux heterogenes complexes. En theorie, la limite de resolution que l'on peut atteindre avec cette technique est de l'ordre de grandeur de la plus petite longueur d'onde presente dans le champ d'onde utilise. Du fait de son cout eleve en temps de calcul, l'inversion de formes d'ondes completes constituait encore recemment un formidable defi pour le sismologue. Cependant, cette situation est en train d'evoluer rapidement du fait des progres recents a la fois des moyens de calcul haute performance ainsi que des methodes numeriques, mais aussi des deploiements de reseaux denses a l'echelle regionale. Dans cette these, nous nous sommes interesses a l'imagerie haute resolution des structures lithospheriques sous la chaine des Pyrenees par inversion de formes d'ondes P telesismiques courte periode. L'objectif principal est d'apporter des contraintes nouvelles sur le taux de raccourcissement subi par cette chaine de montagnes pendant la convergence alpine. Nous utilisons une methode de modelisation hybride qui couple une methode de propagation d'onde globale 1D a une methode d'elements spectraux 3D a l'echelle regionale. Cette methode hybride permet de coupler les champs globaux et regionaux sur les bords du domaine regional 3D. Elle limite les calculs 3D qui sont extremement couteux a l'interieur du domaine regional, ce qui permet de reduire considerablement le temps de calcul. La methode hybride permet ainsi de modeliser des sismogrammes synthetiques jusqu'a des periodes de l'ordre de la seconde, en prenant en compte toutes les complexites qui peuvent affecter la propagation des ondes dans le domaine regional 3D. A l'aide de cette methode, il est egalement possible de calculer par la methode de l'adjoint les derivees de Frechet qui relient les perturbations des formes d'onde observees aux perturbations des parametres elastiques et de la densite dans le milieu. Ces noyaux de sensibilite sont utilises pour formuler un probleme inverse resolu par un algorithme iteratif de type L-BFGS. Nous inversons les donnees de 5 sources telesismiques enregistrees par deux transects denses deployes au niveau des Pyrenees occidentales et centrales pendant l'experience PYROPE. Nous avons ainsi obtenu les premieres sections haute resolution de vitesses des ondes P et S au travers d'une chaine de montagnes. Les modeles tomographiques apportent des evidences claires en faveur d'un sous charriage de la plaque iberique sous la plaque europeenne. Ils montrent egalement l'importance de l'heritage et en particulier des structures liees a l'episode de rifting cretace dans la structuration de la chaine.