Basisprincipes en toepassingen van seismische tomografie
Seismische tomografie is een geofysische methode die wordt gebruikt om het binnenste van de aarde te "fotograferen" met behulp van seismische golven. Net zoals een CT-scan in de geneeskunde een driedimensionaal beeld van het lichaam creëert op basis van röntgenstralen, modelleert seismische tomografie structuren in de ondergrond met behulp van gegevens over reistijd, amplitude of golfvormen die worden geregistreerd door een netwerk van seismometers. Deze methode is cruciaal omdat het binnenste van de aarde niet direct kan worden waargenomen, terwijl dynamische processen zoals platentektoniek, vulkanische activiteit en aardbevingen sterk worden beïnvloed door variaties in de fysische eigenschappen van gesteenten op grote diepte.
1. Basisbegrip van seismische golven
Seismische golven zijn elastische golven die zich door gesteente voortplanten. Seismische golven worden over het algemeen onderverdeeld in:
1. Lichaamsgolven
– P-golven (primaire/compressiegolven): planten zich het snelst voort, kunnen door vaste stoffen en vloeistoffen heen gaan en zijn gevoelig voor veranderingen in compressiesnelheid en dichtheid.
– S-golven (secundaire/schuifgolven): langzamer dan P-golven, kunnen zich niet voortplanten in vloeistoffen, gevoelig voor de stijfheid (schuifmodulus) van het gesteente.
2. Oppervlaktegolven
Rayleigh- en Love-golven, die doorgaans dominant zijn in opnames van aardbevingen op grote afstand, vertonen bijvoorbeeld dispersie en zijn zeer informatief voor de structuur van de lithosfeer tot aan de ondiepe asthenosfeer.
De snelheid van seismische golven is afhankelijk van elastische parameters (modulus) en dichtheid. Variaties in temperatuur, mineraalsamenstelling, druk, porositeit, breuken en de aanwezigheid van vloeistoffen of gedeeltelijke smelt beïnvloeden de snelheid. Dit vormt de fysieke basis voor het in kaart brengen van snelheidsvariaties, wat inzicht geeft in de geologische omstandigheden in de ondergrond.
2. Principes van tomografie: Voorwaarts probleem en achterwaarts probleem
Seismische tomografie is gebaseerd op twee rekenconcepten:
a) Voorwaarts probleem (voorwaarts probleem)
Aan de hand van een model van de aarde (bijvoorbeeld een golfsnelheidsverdeling) kunnen we een voorspelling doen van seismische gegevens: de reistijd van een golf van een bron (aardbeving of kunstmatige bron) naar een station. In een eenvoudige benadering met behulp van de straaltheorie wordt aangenomen dat het golfpad een straal is die het principe van Fermat volgt: de golf kiest het pad met de kortste reistijd.
De reistijd \(T\) kan wiskundig als volgt worden geschreven:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
waarbij \(v(\mathbf{x})\) de golfsnelheid is op positie \(\mathbf{x}\), en \(ds\) het padlengte-element is.
b) Inverse probleem
Het omgekeerde is ook waar: we beschikken over observatiegegevens (reistijden, reistijdresidualen, oppervlaktegolfdispersie of golfvormen) en willen het best passende snelheidsmodel schatten. Het inverse probleem is over het algemeen niet uniek en slecht geformuleerd: meerdere modellen kunnen even goed passen en de gegevens bevatten onzekerheid/ruis. Daarom is regularisatie nodig, zoals gladmaken, dempen of geologische beperkingen.
In de praktijk wordt inversie vaak op een gelineariseerde manier uitgevoerd: het oorspronkelijke model wordt lichtjes verstoord, waarna het reistijdverschil als volgt wordt berekend:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
waarbij \(\delta s = \delta(1/v)\) de traagheidsverstoring is. Deze vergelijkingen worden samengevoegd tot een groot lineair systeem \( \mathbf{d} = \mathbf{Gm} \), dat vervolgens wordt opgelost met de gedempte kleinste-kwadratenmethode of varianten daarvan.
3. Soorten seismische tomografie
1) Reistijdtomografie
De meest gebruikte methode maakt gebruik van P- en S-fase aankomsten van vele aardbevingen. De data zijn reistijdresidualen ten opzichte van een referentiemodel (bijvoorbeeld IASP91 of ak135 op wereldschaal). Geschikt voor 3D-snelheidsmodellering in de aardkorst en de bovenste mantel, vooral als het stationsnetwerk dicht is.
2) Oppervlaktegolftomografie
Door gebruik te maken van de dispersie (verandering in fase-/groepsnelheid met de periode) van Rayleigh-/Love-golven. Deze tomografie is gevoelig voor ondiepe tot middeldiepe structuren en is zeer nuttig voor het in kaart brengen van de dikte van de lithosfeer, zones met lage snelheden (LVZ's) en temperatuurvariaties.
3) Teleseismische tomografie (teleseismische tomografie)
Door gebruik te maken van verre aardbevingen (teleseismische aardbevingen) waarvan de golven door de mantel reizen en vervolgens worden geregistreerd door een lokaal netwerk. Het voordeel: ze hebben meerdere bronnen en komen uit verschillende richtingen, waardoor het ondergrondse volume van het onderzoeksgebied, bijvoorbeeld onder vulkanen of subductiezones, beter in kaart wordt gebracht.
4) Golfvormtomografie / volledige golfvorminversie (FWI)
Het gebruik van de volledige golfvorm, en niet alleen de aankomsttijd, levert theoretisch een hoge resolutie op omdat het de amplitude- en fase-informatie volledig benut, maar vereist wel veel rekenkracht en een goed initieel model om valkuilen van lokale minima te vermijden.
4. Algemene fasen van verwerking en inversie
1. Gegevensverzameling
Het verzamelen van opnames van permanente of tijdelijke seismometers. Op exploratieschaal kan de bron een explosie of vibroseisme zijn; op regionale tot mondiale schaal is de primaire bron een aardbeving.
2. Fase-identificatie en picking
Het bepalen van de aankomsttijd van P-, S- of oppervlaktegolven. De kwaliteit van de meting bepaalt de kwaliteit van het model.
3. Initiële correctie en modellering
Tijdcorrectie (klokafwijking), hoogtecorrectie en selectie van het referentiemodel. Het initiële model kan 1D of eenvoudig 3D zijn.
4. Straaltracering / golfsimulatie
Bereken het pad van synthetische stralen of golfvormen om de gevoeligheidsmatrix samen te stellen.
5. Inversie en regularisatie
Het systeem wordt opgelost om een snelheidsmodel te verkrijgen. De regularisatie wordt gekozen op basis van het doel: scherpe anomalieën of vloeiende trends benadrukken.
6. Resolutie-evaluatie en betrouwbaarheidstesten
Bijvoorbeeld met behulp van de schaakbordtest, de piektest of de covariantie-/puntspreidingsfunctieanalyse om te zien welke delen daadwerkelijk door de data worden opgelost.
5. Interpretatie van het snelheidsmodel
Het resultaat van tomografie is meestal een kaart van relatieve snelheidsanomalieën: zones met hoge snelheden worden vaak geïnterpreteerd als koelere, dichtere of stijvere gesteenten (bijvoorbeeld subducerende platen). Zones met lage snelheden kunnen wijzen op hoge temperaturen, veranderde gesteenten, met vloeistof gevulde breukzones of gedeeltelijke smelting – vaak geassocieerd met actieve vulkanen of een hete asthenosfeer.
De interpretatie moet echter met de nodige voorzichtigheid gebeuren, omdat de snelheid door vele factoren wordt beïnvloed. Idealiter wordt seismische tomografie gecombineerd met andere gegevens, zoals zwaartekrachtmetingen, magnetotellurische metingen (MT), geodetische metingen, oppervlaktegeologie en petrologie.
6. Toepassingen van seismische tomografie
a) Studie van subductiezones en plaatdynamiek
Globale en regionale tomografie kan de subductiezone van de mantel in kaart brengen, inclusief de geometrie, diepte en segmentatie ervan. Deze informatie is cruciaal voor het begrijpen van de oorzaken van grote aardbevingen, de koppelingsmechanismen tussen platen en de tektonische evolutie van een regio.
b) Vulkanische systemen en rampenbestrijding
In vulkanische gebieden kan tomografie zones met lage seismische snelheden identificeren die verband houden met magmakamers, opstijgende vloeistofstromen of veranderde gesteenten. Door regelmatige monitoring kunnen snelheidsveranderingen (time-lapse tomografie) indicatoren leveren voor vloeistof-/magmabewegingen die relevant zijn voor een vroegtijdig waarschuwingssysteem voor erupties.
c) Energie- en grondstoffenexploratie
Op de ondiepe schaal van de aardkorst wordt tomografie gebruikt bij de exploratie van koolwaterstoffen, geothermie en mijnbouw. Variaties in stroomsnelheid helpen bij het in kaart brengen van lithologie, breukstructuren, breukzones en reservoirs. Bij geothermie is de integratie van tomografie met MT vaak effectief: tomografie levert informatie over de elasticiteit, MT levert informatie over de vloeistofgeleidbaarheid.
d) Karakterisering van actieve breuklijnen en aardbevingsrisico's
Lokale tomografie kan zwakke zones, beschadigingszones en heterogeniteit rondom breuken aan het licht brengen. Dit helpt bij het begrijpen van breuksegmentatie, mogelijke blokkering en snelheidsvariaties die de golfversterking beïnvloeden (locatie-effecten).
e) Regionale korst- en lithosfeerstructuur
Met behulp van oppervlaktegolftomografie kunnen onderzoekers de dikte van de lithosfeer, de Moho-grens en asthenosferische anomalieën in kaart brengen. De resultaten vormen de basis voor geodynamische modellen, waaronder bekkenvorming, orogenese en continentale evolutie.
7. Beperkingen en uitdagingen
Seismische tomografie is sterk afhankelijk van de verdeling van bronnen en stations. Gebieden met weinig aardbevingen of een lage dichtheid aan stations hebben een slechte belichting, wat resulteert in een lage resolutie. Bovendien kunnen de aannames van de straaltheorie minder nauwkeurig zijn bij sterke heterogeniteit of bepaalde frequenties. Ruis, selectiefouten en onzekerheid over de hypocentrumlocatie kunnen ook de inversieresultaten beïnvloeden. Daarom zijn resolutie-evaluatie en integratie van meerdere methoden cruciaal om overinterpretatie te voorkomen.
Sluitend
Seismische tomografie is een essentieel instrument voor het onderzoeken van de interne structuur van de aarde, van lokale tot mondiale schaal. De theoretische basis ervan berust op de relatie tussen de elastische eigenschappen van gesteente en golfvoortplanting, evenals op het oplossen van inverse problemen die regularisatie en resolutietesten vereisen. De toepassingen zijn breed: van het in kaart brengen van subducerende platen en vulkanische magmasystemen tot energie-exploratie en het beoordelen van aardbevingsrisico's. Dankzij de vooruitgang in seismische netwerken, computertechnologie en inversiemethoden (waaronder volledige golfvorminversie) is seismische tomografie steeds beter in staat om scherpere en meer informatieve beelden te produceren, zowel voor wetenschappelijke doeleinden als voor rampenbestrijding.
Indien gewenst kan ik dit artikel bewerken tot een meer academische versie (met bronvermelding en een bibliografie), of me richten op een van de toepassingen (bijvoorbeeld tomografie voor vulkanen of voor geothermische energie).
