Es wird Zeit, einmal etwas von der Methode zu berichten, mit der ich mich hauptsächlich in meiner Arbeit befasse. Zunächst nur als grundsätzliche Erklärung, aber da wir letzte Woche eine Feldmessung hatten, kann ich auch noch ein angewandtes Beispiel zeigen in einem Folgepost.
Die Problematik ist diese: Wir haben einen Boden bzw eine Oberfläche
in die wir hineinsehen wollen. Beispiele gibt es viele: Die
geophysikalischen Methoden kommen weitgehend aus der Suche nach
Ölvorkommen, aber das ist eine andere Skala und erfordert andere
Methoden als oberflächennahe Methoden mit denen ich mich befasse.
Außerdem würde ich als Umweltphysiker nie nach Öl suchen!
Bodennah gibt es auch Anwendungen, die meistens mit einer
Strukturerkundung des Untergrundes befasst sind. Das kann z.B. die
Suche nach vergrabenen Gegenständen sein wie Rohren oder
archäologischen Funden, oder nach Hohlräumen oder Höhlen, oder einfach
eine Erkundung der Schichten im Untergrund. Meine Interesse liegt bei
Flüssigkeiten im Boden. Da kann man ein Interesse haben,
Schadstoffausbreitung im Boden zu beobachten, oder Untersuchungen an
einem Grundwasserleiter durchführen. Oder einfach den Wassergehalt des
Bodens ermitteln.
Messmethode
Das grundsätzliche Prinzip der Electrical Resistivity Tomography (ERT) oder einfach Geoelektrik ist einfach zu beschreiben:
Man steckt vier Elektroden in den Boden. Dann speist man an zwei
Elektroden, die man im Allgemeinen mit A und B bezeichnet, einen Strom
in den Boden ein. Die roten Linien deuten die Stromlinien an, entlang
denen der Strom fließt. Man arbeitet hier mit Gleichstrom (es gibt auch
die Variante der Induzierten Polarisation, die auch frequenzabhängig
misst). An den weiteren zwei Elektroden misst die Spannung und kann
daraus den Wert des elektrischen Widerstandes ermitteln.
Im
Bild sieht man einen homogenen Untergrund, aber natürlich ist ein Boden
nicht homogen. Daher reicht nicht eine Messung dieser Art. Den
gemessenen elektrischen Widerstand ordnet man hier einer (Pseudo)tiefe
zu, die ungefähr 1/5tel der Distanz zwischen den äußeren Elektroden
entspricht. Dann misst man mit anderem Abstand zwischen den Elektroden
nochmal und erreicht eine andere Tiefe. So setzt man sich nach und nach
ein Tomogramm des Untergrundes zusammen.
Das geht aber heutzutage
einfacher als die Elektroden umzustecken und immer einzeln zu messen.
Stattdessen steckt man eine ganze Reihe von Elektroden, und misst für
alle möglichen Kombinationen von Elektroden auf diese Art. Dazu gibt es
Multiplexer-Kabel, die eine große Anzahl Elektroden verbinden (z.B. 24
oder 48) und eine Messbox, die automatisch die richtigen Elektroden
ansteuert.
Man misst nicht alle Kombinationen, das würde zu viel
werden bei großen Arrays. Es gibt einige etablierte Arrays, am
häufigsten trifft man diese an:
Wenner-Schlumberger befinden sich die Einspeiseelektroden A und B
außen. Die Potentialelektroden M und N sitzen in der Mitte, und der
Abstand von A nach M und B nach N ist gleich, und jeweils ein ganzes
Vielfaches des Abstandes a zwischen M und N.
Dipol-Dipol sieht
anders aus, hier ist der Abstand jeweils zwischen A und B und M und N
klein, also man hat jeweils einen Einspeise-Dipol und einen
Potential-Dipol mit einer Entfernung n*a.
Eine typische Messung mit 24 Elektroden hätte z.B. zwischen 300 und 400 dieser Quadrupole.
Inversion
Aber
damit ist man (leider) nicht fertig, denn man hat dann kein Bild der
Widerstandsverteilung im Untergrund. Stellt euch vor, jetzt haben wir
für jede unserer z.B. 300 Messungen zwar diese Messung einem Punkt im
Untergrund zugewiesen, aber wir haben angenommen dass die Stromlinien
so schön verlaufen wie beim homogenen Untergrund. Das stimmt natürlich
nicht, weil genau die Widerstandsvertielung, bzw die elektrischen
Leitfähgikeitsstrukturen im Boden, die uns interessieren, die
Stromlinien beeinflussen.
Daher müssen wir auf ein mathematisches
Verfahren zurückgreifen, um ein Modell der Widerstandsverteilung im
Untergrund zu erzeugen. Stellen wir uns dafür vor, dass wir den Boden
in kleine Einheiten aufteilen, die je einen einheitlichen elektrischen
Widerstand haben, und versuchen diesen Wert anzupassen. Wir können
Dreiecke wählen oder einfach rechteckige Kistchen.
Fangen wir mit
einem einfachen Versuch an: Wählen wir den gleichen Widerstand für alle
Kisten. Jetzt kann man einen Schritt machen, den man Vorwärtsmodellierung
nennt. Die Daten, das ist unsere Messung. Das Modell, das ist die Liste
der Widerstände unserer Blöcke. Wenn wir jetzt ausrechnen (durch ein
Computermodell), was wir messen würden, wenn wir tatsächlich auf einer
solchen Widerstandsverteilung messen würden, dann erhalten wir
modellierte Daten.
Die Aufgabe der Inversion ist es, den
Unterschied zwischen modellierten und gemessenen Daten zu minimieren.
Dazu macht sie immer wieder solche Vorwärtsmodellierungen und ändert
nach jeder das Modell, also die Widerstände pro Kiste, passend um,
sodass der Abstand zwischen modellierten und gemessenen Daten kleiner
wird.
Eine klassische least-squares-Minimierung setzt z.B. auf
diese Art an (ich erspare mir mal die Formel): Man verwendet sogenannte
Jacobi-Matrizen, die partielle Ableitungen der Daten nach den
Modellparametern enthalten. Das heißt nichts anderes als: Ich suche die
Richtung in die ich besonders an dem Modellparametern (den Widerständen
der Blöcke in meinem Modell des Untergrundes) drehen muss, um eine
bessere Anpassung zu erreichen.
Leider ist das Problem nicht
eindeutig, im Gegenteil, es existieren unendlich viele Lösungen. Manch
ein Block mag die Messung an vielen Quadrupolen beeinflusst haben,
manch einer gar keine oder zu wenige für eine eindeutige Lösung. Das
nennt man ein gemischt-bestimmtes Problem. Man muss zusätzliche
Annahmen hineinstecken, um den Prozess der Inversion, also die Rechnung
die der Computer durchführt, numerisch stabil zu halten. Eine beliebte
Annahme ist die der Glattheit, dass also keine großen Sprünge zwischen
benachbarten Blöcken auftauchen. Diese Annahme ist beliebig schlecht,
für qualitative Anwendungen aber völlig ausreichend gut. Für eine
Quantifizierung der Widerstände, die man beispielsweise braucht wenn
man Wassergehalte messen will, kann man damit aber Probleme bekommen.
Weswegen ich Forschungsarbeit habe.
Im nächsten Teil zeige ich dann Beispiele dazu!
Kommentare (6)