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

Wanneer de aarde schudt, draagt zelfs een kleine trilling een verborgen verhaal in zich over hoe de grond brak. Seismologen zijn al lang in staat om te bepalen waar een aardbeving begon en hoeveel energie deze vrijgaf, maar het begrijpen van de exacte wijze waarop die energie in de loop van de tijd werd vrijgegeven, bleef een moeilijke puzzel. Dit verhaal wordt verteld door iets dat de 'source time function' wordt genoemd, een registratie van de hartslag van de breuk, van de eerste scheur tot de laatste stop. Het kennen van dit ritme is cruciaal omdat het de snelheid van de breuk, de richting waarin deze zich voortbewoog en de spanningsveranderingen die diep onder de grond plaatsvonden, onthult. Het extraheren van dit ritme uit de chaotische golven die aan het oppervlak worden geregistreerd, is echter berucht moeilijk. Traditionele methoden vertrouwen ofwel op het vinden van een kleinere, vergelijkbare aardbeving om als referentie te gebruiken — een proces dat subjectief en onnauwkeurig kan zijn — of ze proberen massieve, computerintensieve simulaties uit te voeren die te traag zijn voor routinegebruik.

In een nieuwe studie hebben onderzoekers Wojciech Dębski en Kamil Waśkiewicz van het Instituut voor Geofysica van de Poolse Academie van Wetenschappen een andere manier voorgesteld om naar deze verborgen ritmes te luisteren. Ze ontwikkelden een techniek gebaseerd op het eenvoudige fysieke idee dat geluid en seismische golven net zo goed achteruit in de tijd kunnen reizen als vooruit. Door de opnames aan het oppervlak in de tijd om te keren en ze met behulp van een computermodel weer de grond in te sturen, focussen de golven zich van nature weer op de plek waar de aardbeving begon. Terwijl deze omgekeerde golven samenkomen, reconstrueren ze het oorspronkelijke patroon van energievrijgave. Het team testte dit idee met behulp van computersimulaties van kleine, door de mens veroorzaakte aardbevingen, zoals die getriggerd door mijnbouw, met behulp van ruisvrije synthetische data. Hun resultaten laten zien dat deze methode de gedetailleerde timing van de breuk succesvol kan herstellen zonder complexe wiskundige inversies of het zoeken naar een referentieaardbeving nodig te hebben. Het biedt een snellere, robuustere manier om de mechanica van kleine seismische gebeurtenissen te begrijpen, door een moeilijk wiskundig probleem te transformeren in een rechttoe eigenlijk fysiek proces.

De kern van deze nieuwe aanpak rust op een fundamentele eigenschap van de golfmechanica: tijdsomkeerbaarheid. In een perfect elastisch medium, waar energie niet verloren gaat aan warmte of wrijving, werken de vergelijkingen die beschrijven hoe golven bewegen op dezelfde manier of de tijd nu vooruit of achteruit beweegt. Stel je een video voor van een steen die in een vijver valt; als je de video achteruit afspeelt, reizen de rimpelingen weer naar binnen en springt de steen uit het water. Hoewel je de tijd in de echte wereld niet daadwerkelijk kunt terugdraaien, staat de wiskunde toe dat wetenschappers dit proces simuleren. De onderzoekers gebruikten dit concept om een algoritme te creëren dat zij TRSTF noemen. Het proces begint met het nemen van de seismische golven die op verschillende stations aan het oppervlak zijn geregistreerd. Deze golven worden vervolgens in de tijd omgekeerd en teruggevoerd in een computermodel van de aardkorst. Terwijl de omgekeerde golven door het model reizen, stuiteren ze tegen grenzen en interageren ze met de complexe geologie, om uiteindelijk weer samen te komen op de oorspronkelijke locatie van de bron.

Wanneer deze omgekeerde golven elkaar bij de bron ontmoeten, interfereren ze met elkaar. Op de exacte locatie van de aardbeving lijnen de golven zich perfect uit, wat een sterk, gefocust signaal creëert. Op elke andere locatie komen de golven uit de pas aan en heffen ze elkaar op, waardoor alleen een zwakke ruis overblijft. Dit focuseffect is de sleutel tot de methode. Het signaal dat bij het brandpunt tevoorschijn komt, is een directe reconstructie van de source time function, het eigenlijke ritme van de energievrijgave van de aardbeving. De onderzoekers ontdekten dat deze techniek opmerkelijk goed werkt in hun ruisvrije synthetische experimenten, zelfs wanneer het model van de binnenkant van de aarde niet perfect bekend is. In hun computerexperimenten simuleerden ze drie verschillende scenario's om monitoringopstellingen in de echte wereld na te bootsen: één waarbij sensoren nabij het niveau van de mijnbouwactiviteit waren geplaatst, één waarbij de sensoren aan het oppervlak boven het evenement waren geplaatst, en een derde met een andere oppervlakteopstelling. In elk geval slaagde de methode erin om de vorm van het oorspronkelijke breukpatroon te herstellen.

De studie onderzocht ook hoe fouten in het model van de geluidssnelheid in de aarde de resultaten beïnvloeden. In de echte wereld kennen geologen niet de exacte snelheid van seismische golven op elk punt onder de grond; ze vertrouwen op schattingen. De onderzoekers testten hun methode door willekeurige fouten en blokvormige variaties in hun computermodellen te introduceren, om een situatie te simuleren waarin het snelheidmodel slechts ongeveer 90 procent accuraat was. De resultaten gaven aan dat het gebruik van onnauwkeurige snelheidsmodellen systematische fouten introduceert in de gereconstrueerde source time function. De methode bleek gevoeliger voor georganiseerde, blokvormige fouten dan voor willekeurige, verspreide fouten. Ondanks deze fouten bleef het gereconstrueerde signaal herkenbaar en leverde het nuttige informatie, wat aantoont dat de methode bruikbaar blijft zelfs met imperfecte modellen, hoewel de nauwkeurigheid niet immuun is voor modelonnauwkeurigheden.

Een van de meest opvallende bevindingen was hoe de opstelling van de sensoren de kwaliteit van het resultaat beïnvloedde. De onderzoekers ontdekten dat het hebben van sensoren dichter bij de grenzen van het bestudeerde gebied de reconstructie feitelijk hielp. Dit komt omdat de reflecties van deze grenzen fungeren als aanvullende virtuele sensoren, die meer informatie aan het systeem leveren. In een van hun testgevallen, waarbij de sensoren nabij de bovenrand van het model waren geplaatst, was het gereconstrueerde signaal scherper en nauwkeuriger dan in de andere opstellingen, ook al zag de ruwe data er door interferentie meer verstoord uit. Dit contra-intuïtieve resultaat suggereert dat de methode profiteert van de complexiteit van het golfveld, door reflecties als een voordeel te gebruiken in plaats van ze als ruis te behandelen.

De onderzoekers maten ook specifieke kenmerken van de gereconstrueerde signalen, zoals hoe snel de breuk begon en hoe lang deze duurde. Ze vonden dat, hoewel de gereconstrueerde signalen iets breder waren dan de originelen — waarschijnlijk door het beperkte aantal sensoren in de simulatie — de algemene vorm en timing behouden bleven. Het vermogen om de duur van de breuk te schatten is essentieel omdat het wetenschappers helpt de dynamische spanningsval (dynamic stress drop) te berekenen, een maatstaf voor hoeveel spanning er tijdens het evenement werd vrijgegeven. Deze informatie is essentieel voor het begrijpen van de mechanica van geïnduceerde seismiciteit, zoals aardbevingen veroorzaakt door mijnbouw of vloeistofinjectie. De studie toonde aan dat de nieuwe methode deze parameters kan schatten met een nauwkeurigheid die vergelijkbaar met, of beter dan, bestaande technieken is, maar met veel minder computationele inspanning.

In tegen tegenstelling tot volledige golfvorminversie (full-waveform inversion), die duizenden simulaties en complexe optimalisatielussen vereist om de beste match te vinden, heeft deze nieuwe methode slechts één enkele passage van de omgekeerde golven door het model nodig. Dit maakt het computationeel efficiënt, in staat om te draaien op standaard desktopcomputers voor kleine modellen of op grote supercomputers voor regionale schalen. Het voorkomt ook de noodzaak om een referentieaardbeving te selecteren, een stap die in traditionele methoden voor menselijke bias en onzekerheid kan zorgen. De onderzoekers benadrukken dat hun aanpak bijzonder geschikt is voor kleine seismische gebeurtenissen, waarbij de bron klein genoeg is om als een enkel punt te worden behandeld. Ze erkennen dat de methode ervan uitgaat dat de grond zich gedraagt als een perfect elastisch materiaal, wat betekent dat het geen energie absorbeert, en dat de bron klein is in verhouding tot de golflengte van de golven. Dit zijn redelijke aannames voor veel mijnbouwgerelateerde gebeurtenissen, maar de onderzoekers merken op dat verder onderzoek nodig is om te begrijpen hoe de methode presteert in complexere, echtere omgevingen met aanzienlijk energieverlies.

Het artikel concludeert door te benadrukken dat deze techniek de potentie heeft om een standaardinstrument in de seismologie te worden. Door een snelle, betrouwbare en fysiek intuïtieve manier te bieden om de source time function terug te halen, opent het de deur naar routinematige analyse van kleine seismische gebeurtenissen die voorheen te moeilijk in detail te karakteriseren waren. De onderzoekers suggereren dat deze methode geïntegreerd kan worden in bestaande monitoringsystemen om direct inzicht te geven in de aard van geïnduceerde seismiciteit. Hoewel de huidige studie beperkt was tot computersimulaties, is de theoretische basis solide en de resultaten veelbelovend. De volgende stap, volgens de auteurs, is het testen van de methode op echte seismische data en het verfijnen van het begrip van de beperkingen en onzekerheden ervan. Voor nu staat de studie als een bewijs van concept dat tijdreversie een krachtige lens kan zijn om terug te kijken in het verleden van de aarde, om het verborgen verhaal te onthullen van hoe de grond breekt.