610 likes | 773 Views
Daniel Köhn Kiel, den 17. Januar 2005. Einführung in die Inversionstheorie und Regularisierung. Was ist Inversion ?. Was ist Inversion ?. Physiker. Was ist Inversion ?. Physikalische Messung. Was ist Inversion ?. Physikalisches Modell. Test des Modells. Was ist Inversion ?.
E N D
Daniel Köhn Kiel, den 17. Januar 2005 Einführung in die Inversionstheorie und Regularisierung
Was ist Inversion ? Physiker
Was ist Inversion ? Physikalische Messung
Was ist Inversion ? Physikalisches Modell
Was ist Inversion ? Modellparameter Meßdaten
Vorwärtsmodellierung Rate Modellparameter
Vorwärtsmodellierung Lange Iterationszeit beigroßem Parameterraum
Inversion Modell- parameter: m Physikalisches Modell: g Meßdaten: b Inversion: m = g-1(bobs) Vorwärtsmodellierung: bmod = g(m)
Beispiel: Bestimmung des Temperaturverlaufes in einem Bohrloch Physikalisches Modell:
Beispiel: Bestimmung des Temperaturverlaufes in einem Bohrloch Setze Messwerte in Modell ein: Messdaten Modellparameter
Beispiel: Bestimmung des Temperaturverlaufes in einem Bohrloch Lösung:
Lineare Inversionsprobleme Exakt bestimmtes Problem: Es existieren exakt soviele Messungen, wie unbekannte Modellparameter =>Quadratische Koeffizientenmatrix
Lineare Inversionsprobleme Überbestimmtes Problem: Es existieren mehr Messungen, als unbekannte Modellparameter =>Koeffizientenmatrix ist nicht quadratisch
Lineare Inversionsprobleme Lösung eines Überbestimmten Problems: Residuum e = Abweichung zwischen gemessenen und modellierten Daten Gauss: “Minimiere Summe der Quadrate des Residuums” Objektfunktion E(m)
Lineare Inversionsprobleme Damit folgen die optimalen Lösungsparameter x zu: Gauss-Newton Verfahren
Gauss-Newton: Beispiel 2 1D Love-Wellen Inversion
Gauss-Newton: 1D Love-Welleninversion Ausbreitung von Love-Wellen
Gauss-Newton: 1D Love-Welleninversion Entstehung von Love Wellen
Gauss-Newton: 1D Love-Welleninversion Bestimmung von vs(z) aus den gemessenen Phasengeschwindigkeiten vph(T) A vs = vph Phasengeschwindigkeits-residuen Untergrundmodell S-Wellengeschwindigkeits-residuen
Gauss-Newton: 1D Love-Welleninversion MOHO Asthenosphäre Oberer Mantel Startmodell vs(z)
Vorwärtsmodellierung Avs = vph
Gauss-Newton: 1D Love-Welleninversion Dvph Startmodell vph(z)
Inversion mit Gauss-Newton Dvs = (ATA)-1ATDvph
Gauss-Newton: 1D Love-Welleninversion Lösung vs(z) nach einem Iterationsschritt
Gauss-Newton: 1D Love-Welleninversion Problem: Schlecht konditioniertes Gleichungssystem.
Regularisierung Regularisierung: Die Untersuchung des Lösungsverhaltens und die anschließende Lösung eines schlecht konditionierten Problems.
Regularisierung Singulärwertzerlegung (SVD): V = Matrix aus den Eigenvektoren von AATU = Matrix aus den Eigenvektoren von ATA l = Eigenwerte von A