Retrieving the Seismic Moment Release Rate of Small Seismic Events Using the Adjoint Method - Theory and Numerical Validation

Lorsque la Terre tremble, même un infime séisme porte en lui une histoire cachée sur la manière dont le sol s'est rompu. Les sismologues sont capables depuis longtemps de localiser l'origine d'un tremblement de terre et de déterminer la quantité d'énergie qu'il a libérée, mais comprendre exactement comment cette énergie a été libérée au fil du temps est resté un casse-tête difficile. Cette histoire est racontée par ce que l'on appelle la fonction de temps de la source, un enregistrement du rythme de la rupture, de sa première fissure à son arrêt final. Connaître ce rythme est crucial car il révèle la vitesse de la rupture, la direction de son trajet et les changements de contrainte qui se sont produits profondément sous terre. Cependant, extraire ce rythme des ondes désordonnées enregistrées à la surface est notoirement difficile. Les méthodes traditionnelles reposent soit sur la recherche d'un séisme plus petit et similaire pour l'utiliser comme référence — un processus qui peut être subjectif et imprécis — soit sur la tentative d'exécuter des simulations massives et gourmandes en ressources informatiques, trop lentes pour une utilisation de routine.

Dans une nouvelle étude, les chercheurs Wojciech Dębski et Kamil Waśkiewicz, de l'Institut de Géophysique de l'Académie Polonaise des Sciences, ont proposé une autre façon d'écouter ces rythmes cachés. Ils ont développé une technique basée sur l'idée physique simple que le son et les ondes sismiques peuvent voyager vers l'arrière dans le temps tout comme ils voyagent vers l'avant. En prenant les enregistrements de la surface, en les inversant dans le temps et en les renvoyant dans le sol à l'aide d'un modèle informatique, les ondes se concentrent naturellement sur l'endroit où le tremblement de terre a commencé. À mesure que ces ondes inversées convergent, elles reconstruisent le motif original de la libération d'énergie. L'équipe a testé cette idée en utilisant des simulations informatiques de petits séismes d'origine humaine, tels que ceux déclenchés par l'activité minière, en utilisant des données synthétiques sans bruit. Leurs résultats montrent que cette méthode peut récupérer avec succès la chronologie détaillée de la rupture sans nécessiter d'inversions mathématiques complexes ou la recherche d'un séisme de référence. Elle offre un moyen plus rapide et plus robuste de comprendre la mécanique des petits événements sismiques, transformant un problème mathématique difficile en un processus physique direct.

Le cœur de cette nouvelle approche repose sur une propriété fondamentale de la physique des ondes : la réversibilité temporelle. Dans un milieu parfaitement élastique, où l'énergie n'est pas perdue sous forme de chaleur ou de friction, les équations qui décrivent le mouvement des ondes fonctionnent de la même manière que le temps avance ou recule. Imaginez une vidéo d'une pierre tombant dans un étang ; si vous jouez la vidéo à l'envers, les ondulations remontent vers l'intérieur et la pierre jaillit hors de l'eau. Bien que vous ne puissiez pas réellement inverser le temps dans le monde réel, les mathématiques permettent de simuler ce processus. Les chercheurs ont utilisé ce concept pour créer un algorithme qu'ils appellent TRSTF. Le processus commence par la prise des ondes sismiques enregistrées par diverses stations à la surface. Ces ondes sont ensuite inversées dans le temps et réinjectées dans un modèle informatique de la croûte terrestre. À mesure que les ondes inversées traversent le modèle, elles rebondissent sur les limites et interagissent avec la géologie complexe, finissant par converger vers le lieu d'origine de la source.

Lorsque ces ondes inversées se rejoignent à la source, elles interfèrent entre elles. À l'emplacement exact du tremblement de terre, les ondes s'alignent parfaitement, créant un signal fort et focalisé. À tout autre endroit, les ondes arrivent décalées et s'annulent mutuellement, ne laissant qu'un faible bruit. Cet effet de focalisation est la clé de la méthode. Le signal qui émerge au point focal est une reconstruction directe de la fonction de temps de la source, le rythme même de la libération d'énergie du séisme. Les chercheurs ont constaté que cette technique fonctionne remarquablement bien dans leurs expériences synthétiques sans bruit, même lorsque le modèle de l'intérieur de la Terre n'est pas parfaitement connu. Dans leurs expériences informatiques, ils ont simulé trois scénarios différents pour imiter des configurations de surveillance du monde réel : un cas où les capteurs étaient placés près du niveau de l'activité minière, un autre où les capteurs étaient en surface au-dessus de l'événement, et un troisième avec une disposition de surface différente. Dans chaque cas, la méthode a réussi à récupérer la forme du motif de rupture original.

L'étude a également exploré comment les erreurs dans le modèle de la vitesse du son de la Terre affectent les résultats. Dans le monde réel, les géologues ne connaissent pas la vitesse exacte des ondes sismiques en chaque point du sous-sol ; ils s'appuient sur des estimations. Les chercheurs ont testé leur méthode en introduisant des erreurs aléatoires et des variations de type blocs dans leurs modèles informatiques, simulant une situation où le modèle de vitesse n'était précis qu'à environ 90 %. Les résultats ont indiqué que l'utilisation de modèles de vitesse inexacts introduit des erreurs systématiques dans la fonction de temps de la source reconstruite. La méthode s'est révélée plus sensible aux erreurs organisées de type bloc qu'aux erreurs aléatoires et dispersées. Cependant, malgré ces erreurs, le signal récupéré restait reconnaissable et fournissait des informations utiles, démontrant que la méthode conserve son utilité même avec des modèles imparfaits, bien que la précision ne soit pas à l'abri des inexactitudes du modèle.

L'une des découvertes les plus frappantes fut la manière dont la disposition des capteurs influençait la qualité du résultat. Les chercheurs ont découvert que le fait d'avoir des capteurs plus proches des limites de la zone étudiée aidait en réalité la reconstruction. Cela s'explique par le fait que les réflexions de ces limites agissent comme des capteurs virtuels supplémentaires, fournissant plus d'informations au système. Dans l'un de leurs cas de test, où les capteurs étaient placés près du bord supérieur du modèle, le signal reconstruit était plus net et plus précis que dans les autres configurations, même si les données brutes paraissaient plus déformées en raison des interférences. Ce résultat contre-intuitif suggère que la méthode tire parti de la complexité du champ d'ondes, utilisant les réflexions à son avantage plutôt que de les traiter comme du bruit.

Les chercheurs ont également mesuré des caractéristiques spécifiques des signaux reconstruits, telles que la rapidité avec laquelle la rupture a commencé et sa durée. Ils ont constaté que, bien que les signaux reconstruits soient légèrement plus larges que les originaux, probablement en raison du nombre limité de capteurs utilisés dans la simulation, la forme globale et le minutage étaient préservés. La capacité d'estimer la durée de la rupture est vitale car elle aide les scientifiques à calculer la chute de contrainte dynamique, une mesure de la quantité de contrainte libérée lors de l'événement. Cette information est essentielle pour comprendre la mécanique de la sismicité induite, comme les tremblements de terre causés par l'exploitation minière ou l'injection de fluides. L'étude a montré que la nouvelle méthode pouvait estimer ces paramètres avec une précision comparable, voire supérieure, aux techniques existantes, mais avec un effort de calcul bien moindre.

Contrairement à l'inversion de forme d'onde complète, qui nécessite l'exécution de milliers de simulations et de boucles d'optimisation complexes pour trouver la meilleure adéquation, cette nouvelle méthode ne nécessite qu'un seul passage des ondes inversées à travers le modèle. Cela la rend informatiquement efficace, capable de fonctionner sur des ordinateurs de bureau standards pour des modèles de petite taille ou sur de grands supercalculateurs pour des échelles régionales. Elle évite également la nécessité de sélectionner un séisme de référence, une étape qui peut introduire un biais humain et une incertitude dans les méthodes traditionnelles. Les chercheurs soulignent que leur approche est particulièrement adaptée aux petits événements sismiques, où la source est suffisamment petite pour être traitée comme un point unique. Ils reconnaissent que la méthode suppose que le sol se comporte comme un matériau parfaitement élastique, ce qui signifie qu'il n'absorbe pas d'énergie, et que la source est petite par rapport à la longueur d'onde des ondes. Ce sont des hypothèses raisonnables pour de nombreux événements liés aux mines, mais les chercheurs notent que des études supplémentaires sont nécessaires pour comprendre comment la méthode se comporte dans des environnements réels plus complexes présentant une perte d'énergie significative.

L'article conclut en soulignant le potentiel de cette technique à devenir un outil standard en sismologie. En fournissant un moyen rapide, fiable et physiquement intuitif pour récupérer la fonction de temps de la source, elle ouvre la voie à l'analyse de routine des petits événements sismiques qui étaient auparavant trop difficiles à caractériser en détail. Les chercheurs suggèrent que cette méthode pourrait être intégrée aux systèmes de surveillance existants pour fournir des informations immédiates sur la nature de la sismicité induite. Bien que l'étude actuelle ait été limitée à des simulations informatiques, le fondement théorique est solide et les résultats sont prometteurs. La prochaine étape, selon les auteurs, sera de tester la méthode sur des données sismiques réelles et d'affiner la compréhension de ses limites et de ses incertitudes. Pour l'instant, l'étude constitue une preuve de concept démontrant que la réversibilité temporelle peut être une lentille puissante pour regarder dans le passé de la Terre, révélant l'histoire cachée de la rupture du sol.