Les tremblements de terre sont des catastrophes naturelles spectaculaires, comme en témoignent les tremblements de terre de Sumatra en 2004 et de Tohoku-Oki en 2011. La prévision des tremblements de terre reste l'un des plus grands défis sociétaux dans le domaine des sciences naturelles. Ce projet de recherche tentera de répondre à la question suivante : Dans quelle mesure les tremblements de terre sont-ils prévisibles ? Nous proposons une approche multidisciplinaire articulée autour de trois axes principaux : (i) la prévisibilité déterministe des tremblements de terre dans des failles simples et homogènes, étudiée en reproduisant et en comprenant les phénomènes sismiques en laboratoire, (ii) la prévisibilité déterministe des tremblements de terre dans des failles complexes et hétérogènes, étudiée par des expériences en laboratoire produisant de multiples cycles sismiques sur des failles avec des hétérogénéités contrôlées et (iii) la prévisibilité statistique des tremblements de terre, étudiée en prévoyant la distribution spatiale de la sismicité expérimentale à l'aide de l'apprentissage automatique.
Sur ces 50 dernières années, les séismes ont coûté 800 milliards de dollars — l’essentiel dans les pays développés — et 1.3 millions de vies humaines — l’essentiel dans les pays en développement. Face à ces chiffres, qui ne montrent pas de signe d’inflexion, la conscientisation face au risque continue d’appliquer l’approche classique où les connaissances détenues par les scientifiques sont traduites « vers le bas » pour le public et les décideurs. Une approche inverse, « ascendante », où les citoyens collectent et partagent des informations sur les séismes, pourrait-elle être un modèle alternatif ?
Le projet WIND a pour objectif de développer des méthodes d'imagerie sismique plus particulièrement adaptées à des dispositifs d'acquisition mis en oeuvre avec des parcs de stations sismiques (terrestres ou sous-marines) autonomes. Le fait que le dispositif de réception est fixe par rapport à celui des sources donne la flexibilité nécessaire pour positionner les capteurs par rapport aux sources avec des déports variables de manière à ce que les ondes enregistrées interagissent avec les hétérogénéités du sous sol de manière variée par transmission, réflexion ou diffraction. Le riche éclairage angulaire du sous-sol qui en découle constitue un facteur favorable à une imagerie haute résolution multi-paramètre sous réserve que la méthode d'imagerie puisse assimiler la richesse de l'information contenue dans les mesures sismiques. Dans ce contexte, nous développons des méthodes d'imagerie par inversion de formes d'ondes complètes (FWI: Full Waveform Inversion) qui visent à reconstruire les propriétés constitutives du sous-sol en minimisant une distance entre les données enregistrées et leur analogues simulées numériquement. Trois difficultés majeures sont à surmonter: (1) la nonlinéarité du problème inverse nécessitant de concevoir de nouvelles distances ou stratégies d'optimisation visant à étendre le régime linéaire de l'inversion ou l'espace de recherche. L'approche que nous développons repose sur la méthode de reconstruction de champs d'onde qui introduit une relaxation de l'équation d'onde généré par un terme de rappel aux observations pour mieux ajuster ces dernières lors du problème direct. (2) La non unicité de la solution résultant de l'éclairage incomplet du sous-sol depuis la surface nécessitant d'incorporer des a priori dans le processus inverse sous forme de contraintes ou de régularisations. Nous développons une librairie de régularisation incluant des functions non lisses (localement non différentiables) implémentées avec des algorithmes proximaux. (3) Le coût calculatoire de la simulation numérique des ondes dans des modèles de Terre en trois dimensions nécessitant l'utilisation maitrisée du calcul scientifique haute performance. Nous développons des algorithmes de simulation des ondes dans les domaines temps-espace et fréquence-espace où le problème numérique à résoudre est de nature très différente (problème d'évolution versus problème de condition aux limites). Dans le domaine fréquentiel, nous utilisons à la fois des méthodes directes fondées sur des solveurs multifrontaux tels que MUMPS et des méthodes hybrides fondées sur des préconditionneurs par décomposition en domaine. Nous développons également des méthodes tomographiques où les pentes d'évènements localement cohérents sont utilisées conjointement au temps de trajet des ondes pour améliorer la résolution des modèles et mieux poser le problème inverse via l'information supplémentaire fournie par les pentes. Ces modèles peuvent être utilisés comme modèles de réference ou modèles initiaux pour la migration avant-sommation profondeur et la FWI. On s'intéresse également au problème du design de l'acquisition qui doit se formuler comme un problème inverse. C'est un enjeu important pour des dispositifs de "nodes" où le nombre d'instruments est nécessairement limité et où il s'agit de minimiser l'empreinte de l'acquisition dans l'imagerie.
Site web : https://www.geoazur.fr/WIND
Financement: consortium de compagnies pétrolières (Chevron, Shell, Total), 405k€ sur 3 ans
PI: Stéphane Operto (email: operto@geoazur.unice.fr; tel: 04 83 61 87 52); Hossein Aghamiry (email: aghamiry@geoazur.unice.fr)
Sponsors: Chevron, Shell, Total.
Collaborations académiques: LJAD (UCA), LJLL (Sorbonne University), ENSEEIHT (Université de Toulouse), MUMPS-consortium, University of Tehran.
Dates : Janvier 2020-Décembre 2023
Les Alpes constituent un laboratoire naturel privilégié pour étudier l'orogenèse en raison de affleurements exceptionnels. La structure profonde des Alpes et les processus tectoniques associés restent néanmoins controversés en raison de la complexité de la chaîne mais aussi de la faible résolution avec laquelle ces niveaux profonds ont été imagés. Pour palier ce manque, le consortium Européen AlpArray a déployé un réseau dense de stations large bande (inter-station < 52km) sur l'ensemble de l'arc Alpin ainsi que 30 stations sismiques sous-marines dans le bassin Ligure. Le projet LisAlps propose d'appliquer la « Full Waveform Inversion » (FWI) aux données télésismiques enregistrées durant AlpArray pour construire :
(1) un nouveau modèle de référence multi-paramétrique de la lithosphère et de l'asthénosphère à l'échelle de l'arc alpin (1500kmx700km) jusqu'à une profondeur de 700km à partir de l'ensemble du réseau et d'un catalogue de 300 séismes,
(2) un modèle lithosphérique haute-résolution des Alpes Occidentales centré sur la zone complexe du noeud Ligure.
UMR Géoazur
Campus Azur du CNRS
250 rue Albert Einstein
- CS 10269 - F 06905 SOPHIA ANTIPOLIS Cedex
+33 (0)483 618 500