Grundlagen und Anwendung der seismischen Tomographie
Die seismische Tomographie ist eine geophysikalische Methode, mit der das Erdinnere mithilfe seismischer Wellen „fotografiert“ wird. Ähnlich wie die Computertomographie (CT) in der Medizin, die aus Röntgenstrahlen ein dreidimensionales Bild des Körpers erzeugt, modelliert die seismische Tomographie unterirdische Strukturen anhand von Laufzeit-, Amplituden- oder Wellenformdaten, die von einem Netzwerk von Seismometern aufgezeichnet werden. Diese Methode ist von entscheidender Bedeutung, da das Erdinnere nicht direkt beobachtet werden kann, während dynamische Prozesse wie Plattentektonik , Vulkanismus und Erdbeben stark von den Veränderungen der physikalischen Eigenschaften von Gesteinen in der Tiefe beeinflusst werden.
1. Grundbegriffe seismischer Wellen
Seismische Wellen sind elastische Wellen, die sich durch Gestein ausbreiten. Seismische Wellen werden im Allgemeinen unterteilt in:
1. Körperwellen
– P-Wellen (Primär-/Kompressionswellen): breiten sich am schnellsten aus, können durch Festkörper und Flüssigkeiten hindurchgehen und reagieren empfindlich auf Änderungen der Kompressionsgeschwindigkeit und der Dichte.
– S-Wellen (Sekundär-/Scherwellen): langsamer als P-Wellen, können sich nicht in Flüssigkeiten ausbreiten, reagieren empfindlich auf die Steifigkeit (Schermodul) des Gesteins.
2. Oberflächenwellen
Beispielsweise weisen Rayleigh- und Love-Wellen, die in Aufzeichnungen von Fernbeben im Allgemeinen dominant sind, eine Dispersion auf und liefern sehr viele Informationen über die Struktur der Lithosphäre bis hinunter zur flachen Asthenosphäre.
Die seismische Wellengeschwindigkeit hängt von elastischen Parametern (Modul) und der Dichte ab. Schwankungen der Temperatur, der Mineralzusammensetzung, des Drucks, der Porosität, von Klüften sowie das Vorhandensein von Fluiden oder Teilschmelzen beeinflussen die Geschwindigkeit. Dies ist die physikalische Grundlage für die Kartierung von Geschwindigkeitsvariationen, um Hinweise auf die geologischen Verhältnisse im Untergrund zu erhalten.
2. Prinzipien der Tomographie: Vorwärtsproblem und Rückwärtsproblem
Die seismische Tomographie basiert auf zwei Berechnungskonzepten:
a) Vorwärtsproblem (Vorwärtsproblem)
Ausgehend von einem Modell der Erde (z. B. der Wellengeschwindigkeitsverteilung) können wir Vorhersagen berechnen. seismische DatenDie Laufzeit einer Welle von der Quelle (Erdbeben oder künstliche Quelle) zur Station. In einem einfachen Ansatz der Strahlentheorie wird angenommen, dass der Wellenweg einem Strahl folgt, der dem Fermatschen Prinzip gehorcht: Die Welle wählt den Weg mit der kürzesten Laufzeit.
Mathematisch lässt sich die Reisezeit \(T\) wie folgt schreiben:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
wobei \(v(\mathbf{x})\) die Wellengeschwindigkeit an der Position \(\mathbf{x}\) ist und \(ds\) das Weglängenelement ist.
b) Inverses Problem
Das Umgekehrte gilt: Wir verfügen über Beobachtungsdaten (Laufzeiten, Laufzeitresiduen, Oberflächenwellendispersionen oder Wellenformen) und möchten das am besten passende Geschwindigkeitsmodell schätzen. Das inverse Problem ist im Allgemeinen nicht eindeutig und schlecht gestellt: Mehrere Modelle können gleich gut passen, und die Daten enthalten Unsicherheiten/Rauschen. Daher ist eine Regularisierung erforderlich, beispielsweise durch Glättung, Dämpfung oder geologische Randbedingungen.
In der Praxis wird die Inversion oft linearisiert durchgeführt: Das Ausgangsmodell wird leicht gestört, und dann wird die Laufzeitdifferenz wie folgt berechnet:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
wobei \(\delta s = \delta(1/v)\) die Langsamkeitsstörung ist. Diese Gleichungen werden zu einem großen linearen System \( \mathbf{d} = \mathbf{Gm} \) angeordnet und anschließend mit der Methode der gedämpften kleinsten Quadrate oder einer ihrer Varianten gelöst.
3. Arten der seismischen Tomographie
1) Laufzeittomographie
Die gängigste Methode nutzt die P- und S-Phasen-Ankünfte zahlreicher Erdbeben. Die Daten stellen Laufzeitresiduen gegenüber einem Referenzmodell dar (z. B. IASP91 oder ak135 auf globaler Ebene). Sie eignet sich für die dreidimensionale Geschwindigkeitsmodellierung in der Erdkruste und im oberen Erdmantel, insbesondere bei einem dichten Stationsnetz.
2) Oberflächenwellen-Tomographie
Die Tomographie nutzt die Dispersion (Änderung der Phasen-/Gruppengeschwindigkeit mit der Periode) von Rayleigh-/Love-Wellen. Da sie empfindlich auf flache bis mitteltiefe Strukturen reagiert, ist sie sehr nützlich zur Kartierung der Lithosphärendicke, von Zonen niedriger Wellengeschwindigkeit (LVZs) und von Temperaturvariationen.
3) Teleseismische Tomographie (teleseismische Tomographie)
Man nutzt teleseismische Fernbeben, deren Wellen sich durch den Erdmantel ausbreiten und von einem lokalen Netzwerk aufgezeichnet werden. Der Vorteil: Sie haben mehrere Ursprungsorte und kommen aus verschiedenen Richtungen, wodurch das Untergrundvolumen des Untersuchungsgebiets, beispielsweise unter Vulkanen oder Subduktionszonen, besser sichtbar wird.
4) Wellenformtomographie / vollständige Wellenforminversion (FWI)
Die Verwendung der kompletten Wellenform, nicht nur der Ankunftszeit, führt theoretisch zu einer hohen Auflösung, da Amplituden- und Phaseninformationen vollständig genutzt werden. Allerdings sind dafür große Rechenanforderungen und ein gutes Ausgangsmodell erforderlich, um lokale Minima zu vermeiden.
4. Allgemeine Verarbeitungs- und Inversionsphasen
1. Akuisisi-Daten
Sammlung von Aufzeichnungen von permanenten oder temporären Seismometern. Im explorativen Maßstab kann die Quelle eine Explosion oder ein Vibroseismus sein; im regionalen bis globalen Maßstab ist die primäre Quelle ein Erdbeben.
2. Phasenidentifizierung und -auswahl
Die Ankunftszeit von P-, S- oder Oberflächenwellen wird bestimmt. Die Qualität der Bestimmung bestimmt die Qualität des Modells.
3. Erste Korrektur und Modellierung
Zeitkorrektur (Uhrabweichung), Höhenkorrektur und Auswahl des Referenzmodells. Das Ausgangsmodell kann 1D oder ein einfaches 3D-Modell sein.
4. Strahlverfolgung / Wellensimulation
Berechne den Verlauf synthetischer Strahlen oder Wellenformen, um die Sensitivitätsmatrix zu erstellen.
5. Inversion und Regularisierung
Das System wird gelöst, um ein Geschwindigkeitsmodell zu erhalten. Die Regularisierung wird je nach Zielsetzung gewählt: ob scharfe Anomalien hervorgehoben oder Trends geglättet werden sollen.
6. Auflösungsbewertung und Zuverlässigkeitsprüfung
Zum Beispiel mit dem Schachbrettmuster-Test, dem Spike-Test oder der Kovarianz-/Punktspreizfunktionsanalyse, um zu sehen, welche Teile durch die Daten tatsächlich aufgelöst werden.
5. Interpretation des Geschwindigkeitsmodells
Das Ergebnis der Tomographie ist üblicherweise eine Karte relativer Geschwindigkeitsanomalien: Zonen mit hoher Geschwindigkeit werden oft als kühlere, dichtere oder steifere Gesteine interpretiert (z. B. subduzierte Platten). Zonen mit niedriger Geschwindigkeit können auf hohe Temperaturen, veränderte Gesteine, mit Fluiden gefüllte Bruchzonen oder partielles Aufschmelzen hinweisen – oft im Zusammenhang mit aktiven Vulkanen oder heißer Asthenosphäre.
Die Interpretation erfordert jedoch Vorsicht, da die Geschwindigkeit von vielen Faktoren beeinflusst wird. Idealerweise wird die seismische Tomographie mit anderen Daten wie Gravimetrie, Magnetotellurik (MT), Geodäsie, Oberflächengeologie und Petrologie kombiniert.
6. Anwendungen der seismischen Tomographie
a) Untersuchung von Subduktionszonen und Plattendynamik
Globale und regionale Tomographie ermöglicht die Kartierung der Mantel-Subduktionsplatte, einschließlich ihrer Geometrie, Tiefe und Segmentierung. Diese Informationen sind entscheidend für das Verständnis der Ursachen großer Erdbeben, der Mechanismen der Plattenkopplung und der tektonischen Entwicklung einer Region.
b) Vulkansysteme und Katastrophenvorsorge
In vulkanischen Gebieten kann die Tomographie Zonen niedriger Fließgeschwindigkeit identifizieren, die mit Magmakammern, aufsteigenden Fluidwegen oder verändertem Gestein in Verbindung stehen. Durch regelmäßige Überwachung können Geschwindigkeitsänderungen (zeitrafferbasierte Tomographie) Indikatoren für Fluid-/Magmabewegungen liefern, die für die Frühwarnung vor Vulkanausbrüchen relevant sind.
c) Energie- und Rohstoffexploration
Im Bereich der oberflächennahen Erdkruste wird die Tomographie in der Kohlenwasserstoff-, Geothermie- und Bergbauexploration eingesetzt. Geschwindigkeitsvariationen helfen bei der Kartierung von Lithologie, Verwerfungsstrukturen, Bruchzonen und Lagerstätten. In der Geothermie ist die Kombination von Tomographie und Magnetometer (MT) oft effektiv: Die Tomographie liefert Informationen zur Elastizität, das MT-Verfahren Informationen zur Fluidleitfähigkeit.
d) Charakterisierung aktiver Verwerfungen und Erdbebengefahren
Die lokale Tomographie kann Schwächezonen, Schadenszonen und Heterogenitäten im Bereich von Verwerfungen aufzeigen. Dies trägt zum Verständnis der Verwerfungssegmentierung, potenzieller Verriegelungen und Geschwindigkeitsvariationen bei, die die Wellenverstärkung beeinflussen (Standorteffekte).
e) Regionale Krusten- und Lithosphärenstruktur
Mithilfe der Oberflächenwellentomographie können Forscher die Dicke der Lithosphäre, die Moho-Grenze und Anomalien der Asthenosphäre kartieren. Die Ergebnisse bilden die Grundlage für geodynamische Modelle, darunter Modelle zur Beckenbildung, Orogenese und Kontinentalentwicklung.
7. Einschränkungen und Herausforderungen
Die seismische Tomographie ist stark von der Verteilung der Quellen und Stationen abhängig. Gebiete mit wenigen Erdbeben oder einer geringen Stationsdichte weisen eine schlechte Ausleuchtung und damit eine niedrige Auflösung auf. Zudem können die Annahmen der Strahlentheorie bei starker Heterogenität oder bestimmten Frequenzen ungenauer sein. Rauschen, Fehler bei der Erdbebenortung und Unsicherheiten bei der Hypozentrumsbestimmung können die Inversionsergebnisse ebenfalls beeinflussen. Daher sind die Bewertung der Auflösung und die Integration verschiedener Methoden entscheidend, um eine Überinterpretation zu vermeiden.
Penutup
Die seismische Tomographie ist ein Schlüsselinstrument zur Erforschung der inneren Struktur der Erde von lokalen bis hin zu globalen Maßstäben. Ihre theoretische Grundlage beruht auf dem Zusammenhang zwischen den elastischen Eigenschaften von Gestein und der Wellenausbreitung sowie auf der Lösung inverser Probleme, die Regularisierung und Auflösungstests erfordern. Ihre Anwendungsgebiete sind vielfältig: von der Kartierung subduzierter Platten und vulkanischer Magmasysteme bis hin zur Energieexploration und Erdbebengefährdungsanalyse. Dank Fortschritten bei seismischen Netzwerken, der Datenverarbeitung und Inversionsmethoden (einschließlich der Vollwellenforminversion ) liefert die seismische Tomographie zunehmend schärfere und aussagekräftigere Bilder für wissenschaftliche Zwecke und die Katastrophenvorsorge.
Auf Wunsch kann ich diesen Artikel in eine wissenschaftlichere Version (mit Zitaten und Bibliographie) umwandeln oder mich auf eine der Anwendungen konzentrieren (z. B. Tomographie für Vulkane oder für Geothermie).
