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

Wenn die Erde bebt, trägt selbst ein winziges Zittern eine verborgene Geschichte darüber in sich, wie der Boden gebrochen ist. Seismologen sind schon lange in der Lage, den Ort eines Erdbebens und die Menge der freigesetzten Energie zu bestimmen, aber genau zu verstehen, wie diese Energie im Zeitverlauf freigesetzt wurde, blieb ein schwieriges Rätsel. Diese Geschichte wird durch etwas namens Quellzeitfunktion erzählt, einer Aufzeichnung des Herzschlags der Ruptur von ihrem ersten Riss bis zu ihrem letzten Stopp. Das Wissen um diesen Rhythmus ist entscheidend, da er die Geschwindigkeit der Ruptur, die Richtung, in die sie wanderte, und die Spannungsänderungen offenbart, die tief unter der Erde stattfanden. Es ist jedoch notorisch schwierig, diesen Rhythmus aus den unruhigen Wellen zu extrahieren, die an der Oberfläche aufgezeichnet werden. Traditionelle Methoden verlassen sich entweder darauf, ein kleineres, ähnliches Erdbeben als Referenz zu finden – ein Prozess, der subjektiv und unpräzise sein kann – oder sie versuchen, massiv rechenintensive Simulationen durchzuführen, die für den Routinegebrauch zu langsam sind.

In einer neuen Studie haben die Forscher Wojciech Dębski und Kamil Waśkiewicz vom Institut für Geophysik der Polnischen Akademie der Wissenschaften einen anderen Weg vorgeschlagen, um diesen verborgenen Rhythmen zu lauschen. Sie entwickelten eine Technik, die auf der einfachen physikalischen Idee basiert, dass Schall und seismische Wellen genauso rückwärts in der Zeit reisen können, wie sie vorwärts reisen. Indem sie die Aufzeichnungen an der Oberfläche nehmen, sie in der Zeit umkehren und sie mithilfe eines Computermodells zurück in den Boden senden, fokussieren sich die Wellen natürlich auf die Stelle, an der das Erdbeben begann. Wenn diese umgekehrten Wellen konvergieren, rekonstruieren sie das ursprüngliche Muster der Energiefreisetzung. Das Team testete diese Idee anhand von Computersimulationen kleiner, menschengemachter Erdbeben, wie sie beispielsweise durch Bergbau ausgelöst werden, unter Verwendung von rauschfreien synthetischen Daten. Ihre Ergebnisse zeigen, dass diese Methode die detaillierte zeitliche Abfolge der Ruptur erfolgreich wiederherstellen kann, ohne komplexe mathematische Inversionen oder die Suche nach einem Referenzerdbebens zu benötigen. Sie bietet einen schnelleren, robusteren Weg, um die Mechanik kleiner seismischer Ereignisse zu verstehen, indem sie ein schwieriges mathematisches Problem in einen geradlinigen physikalischen Prozess verwandelt.

Der Kern dieses neuen Ansatzes beruht auf einer grundlegenden Eigenschaft der Wellenphysik: der Zeitumkehrbarkeit. In einem perfekt elastischen Medium, in dem Energie nicht in Wärme oder Reibung umgewandelt wird, funktionieren die Gleichungen, die die Bewegung von Wellen beschreiben, gleich, ob die Zeit vorwärts oder rückwärts läuft. Stellen Sie sich das Video eines Steins vor, der in einen Teich fällt; wenn Sie das Video rückwärts abspielen, wandern die Wellen zurück nach innen und der Stein springt aus dem Wasser hervor. Obwohl man die Zeit in der realen Welt nicht tatsächlich umkehren kann, erlaubt die Mathematik es Wissenschaftlern, diesen Prozess zu simulieren. Die Forscher nutzten dieses Konzept, um einen Algorithmus zu erstellen, den sie TRSTF nennen. Der Prozess beginnt damit, dass die an verschiedenen Stationen an der Oberfläche aufgezeichneten seismischen Wellen genommen werden. Diese Wellen werden dann in der Zeit umgekehrt und in ein Computermodell der Erdkruste zurückgespeist. Während die umgekehrten Wellen durch das Modell wandern, prallen sie von Grenzen ab und interagieren mit der komplexen Geologie, bis sie schließlich wieder am ursprünglichen Quellort konvergieren.

Wenn diese umgekehrten Wellen am Ursprung zusammentreffen, interferieren sie miteinander. Am exakten Ort des Erdbebens richten sich die Wellen perfekt aus und erzeugen ein starkes, fokussiertes Signal. An jedem anderen Ort kommen die Wellen außertakt an und löschen sich gegenseitig aus, sodass nur ein schwaches Rauschen bleibt. Dieser Fokussierungseffekt ist der Schlüssel zur Methode. Das Signal, das am Fokuspunkt entsteht, ist eine direkte Rekonstruktion der Quellzeitfunktion, des eigentlichen Rhythmus der Energiefreisetzung des Erdbebens. Die Forscher fanden heraus, dass diese Technik in ihren rauschfreien synthetischen Experimenten bemerkenswert gut funktioniert, selbst wenn das Modell des Erdinneren nicht perfekt bekannt ist. In ihren Computerexperimenten simulierten sie drei verschiedene Szenarien, um reale Überwachungssituationen nachzubilden: eines, bei dem die Sensoren in der Nähe der Bergbauaktivität platziert wurden, eines, bei dem die Sensoren an der Oberfläche über dem Ereignis lagen, und ein drittes mit einer anderen Oberflächenanordnung. In jedem Fall stellte die Methode das ursprüngliche Rupturmuster erfolgreich wieder her.

Die Studie untersuchte auch, wie Fehler im Modell der Schallgeschwindigkeit der Erde die Ergebnisse beeinflussen. In der realen Welt kennen Geologen nicht die exakte Geschwindigkeit seismischer Wellen an jedem Punkt unter der Erde; sie verlassen sich auf Schätzungen. Die Forscher testeten ihre Methode, indem sie zufällige Fehler und blockartige Variationen in ihre Computermodelle einführten, um eine Situation zu simulieren, in der das Geschwindigkeitsmodell nur etwa 90 Prozent genau war. Die Ergebnisse zeigten, dass die Verwendung ungenauer Geschwindigkeitsmodelle systematische Fehler in die rekonstruierte Quellzeitfunktion einführt. Die Methode erwies sich als empfindlicher gegenüber organisierten, blockartigen Fehlern als gegenüber zufälligen, gestreuten Fehlern. Trotz dieser Fehler blieb das rekonstruierte Signal jedoch erkennbar und lieferte nützliche Informationen, was zeigt, dass die Methode auch bei unvollkommenen Modellen ihren Nutzen behält, obwohl die Genauigkeit gegenüber Modellungenauigkeiten nicht immun ist.

Eine der eindrucksvollsten Erkenntnisse war, wie die Anordnung der Sensoren die Qualität des Ergebnisses beeinflusste. Die Forscher entdeckten, dass es die Rekonstruktion tatsächlich verbesserte, wenn Sensoren näher an den Grenzen des untersuchten Gebiets platziert wurden. Dies liegt daran, dass die Reflexionen von diesen Grenzen wie zusätzliche virtuelle Sensoren wirken und dem System mehr Informationen liefern. In einem ihrer Testfälle, bei dem die Sensoren nahe der Oberkante des Modells platziert waren, war das rekonstruierte Signal schärfer und genauer als in den anderen Setups, obwohl die Rohdaten aufgrund von Interferenzen stärker verzerrt erschienen. Dieses kontraintuitive Ergebnis deutet darauf hin, dass die Methode von der Komplexität des Wellenfeldes profitiert und Reflexionen nutzt, anstatt sie als Rauschen zu behandない.

Die Forscher maßen auch spezifische Merkmale der rekonstruierten Signale, wie etwa wie schnell die Ruptur begann und wie lange sie dauerte. Sie fanden heraus, dass die rekonstruierten Signale zwar etwas breiter waren als die Originale – wahrscheinlich aufgrund der begrenzten Anzahl der in der Simulation verwendeten Sensoren –, aber die allgemeine Form und der zeitliche Ablauf bewahrt wurden. Die Fähigkeit, die Dauer der Ruptur abzuschätzen, ist entscheidend, da sie Wissenschaftlern hilft, den dynamischen Spannungsabfall zu berechnen, ein Maß dafür, wie viel Spannung während des Ereignisses freigesetzt wurde. Diese Information ist essenziell für das Verständnis der Mechanik induzierter Seismizität, wie etwa Erdbeben, die durch Bergbau oder Flüssigkeitseinpressung verursacht werden. Die Studie zeigte, dass die neue Methode diese Parameter mit einer Genauigkeit schätzen kann, die mit bestehenden Techniken vergleichbar oder sogar besser ist, jedoch mit wesentlich geringerem Rechenaufwand.

Im Gegensatz zur vollen Wellenforminversion (Full-Waveform Inversion), die Tausende von Simulationen und komplexe Optimierungsschleifen erfordert, um die beste Übereinstimmung zu finden, benötigt diese neue Methode nur einen einzigen Durchgang der umgekehrten Wellen durch das Modell. Dies macht sie recheneffizient und in der Lage, auf Standard-Desktop-Computern für kleine Modelle oder auf großen Supercomputern für regionale Skalen zu laufen. Sie vermeidet zudem die Notwendigkeit, ein Referenzerdbeben auszuwählen – ein Schritt, der bei traditionellen Methoden menschliche Voreingenommenheit und Unsicherheit einführen kann. Die Forscher betonen, dass ihr Ansatz besonders gut für kleine seismische Ereignisse geeignet ist, bei denen die Quelle klein genug ist, um als einzelner Punkt behandelt werden zu können. Sie räumen ein, dass die Methode davon ausgeht, dass der Boden sich wie ein perfektes elastisches Material verhält, was bedeutet, dass er keine Energie absorbiert, und dass die Quelle klein im Vergleich zur Wellenlänge der Wellen ist. Dies sind vernünftige Annahmen für viele bergbaurelevante Ereignisse, doch die Forscher merken an, dass weitere Studien notwendig sind, um zu verstehen, wie die Methode in komplexeren, realen Umgebungen mit signifikantem Energieverlust abschneidet.

Das Paper schließt mit dem Hinweis, dass diese Technik das Potenzial hat, ein Standardwerkzeug in der Seismologie zu werden. Indem sie einen schnellen, zuverlässigen und physikalisch intuitiven Weg bietet, die Quellzeitfunktion zurückzugewinnen, öffnet sie die Tür zur routinemäßigen Analyse kleiner seismischer Ereignisse, die zuvor zu schwer detailliert zu charakterisieren waren. Die Forscher schlagen vor, dass diese Methode in bestehende Überwachungssysteme integriert werden könnte, um sofortige Einblicke in die Natur der induzierten Seismizität zu geben. Während die aktuelle Studie auf Computersimulationen beschränkt war, ist das theoretische Fundament solide und die Ergebnisse sind vielversprechend. Der nächste Schritt, so die Autoren, wird darin bestehen, die Methode an realen seismischen Daten zu testen und das Verständnis ihrer Grenzen und Unsicherheiten zu verfeinern. Für den Moment steht die Studie als Beweis der Konzeptstärke dafür, dass die Zeitumkehrierung eine kraftvolle Linse sein kann, um in die seismische Vergangenheit der Erde zu blicken und die verborgene Geschichte darüber zu enthüllen, wie der Boden bricht.