1  Einführung

1.1 Gliederung Teil Geoelektrik

  1. Grundgleichungen und elektrische Felder
  2. Petrophysikalische Grundlagen
  3. Anordnungen, Profilmessungen und Sensitivitäten
  4. Widerstandstiefensondierung
  5. 2D-Geoelektrik, Modellierung und Inversion
  6. 3D-Geoelektrik und crosshole-ERT
  7. Prozessmonitoring und Timelapse-Verfahren
  8. Induzierte Polarisation im Zeit- und Frequenzbereich

1.2 Gliederung Teil Elektromagnetik

  1. Grundlagen der Geoelektromagnetik
  2. Magnetotellurik
  3. CSEM-Methoden im Frequenzbereich
  4. Transientelektromagnetik
  5. Gemeinsame Inversion
  6. Nuklearmagnetische Resonanz

1.3 Empfohlene Lehrbücher

Telford, Geldart, Sherriff (1990): Cambridge

Knoedel, Krummel, Lange (2005): Springer

Everett (2013): Cambridge

1.4 Literatur

1.5 Lernpyramide

  • Vorlesung und Scripte
  • Notebooks zum Spielen
  • Übungen
  • praktische Messungen
  • gemeinsame Auswertung

2 Einführung

2.1 Elektromagnetische Größen

  • magn. Feldstärke \(H\) [A/m]
  • magn. Flussdichte \(B\) [Vs/m²] \(B=\mu H\) (oft “Magnetfeld”)
  • magn. Permeabilität \(\mu\) [Vs/Am] \(\mu=\mu_0\mu_r\) , \(\mu_0\)=4\(\pi\) 1e-7Vs/Am
  • mag. Suszeptibilität \(\chi_m=\mu_r\)-1
  • magn. Fluss \(\Phi\) [Vs]
  • spez. el. Widerstand \(\rho\) [\(\Omega\)m]
  • elektr. Feldstärke \(E\) [V/m]
  • dielektr. Verschiebung \(D\) [As/m²] \(D=\epsilon E\) (Hilfsfeld)
  • elektr. Permittivität \(\epsilon\) [As/Vm] \(\epsilon=\epsilon_0\epsilon_r\) , \(\epsilon_0\)=8.854e-12As/Vm
  • el. Suszeptibilität \(\chi_e=\epsilon_r\)-1
  • elektr. Stromdichte \(j\) [A/m²]
  • el. Leitfähigkeit \(\sigma\) [S/m]

2.2 Maxwell-Gleichungen

\[ \curl\vb E = -\pdv{\vb B}{t} \]

\[ \curl \vb H = \pdv{\vb D}{t} + \vb j \]

2.3 Maxwell - Quellen und Material

Ladungsträger erzeugen eine elektrische Verschiebung \[\div \vb D = \varrho\]

Das Magnetfeld ist frei von Quellen und Senken \[\div \vb B = 0\]

Materialgleichungen: \(\vb B=\mu\vb H\), \(\vb D=\epsilon\vb E\), \(\vb j=\sigma\vb E\)

2.4 Amperesches Gesetz in der Geophysik

\[ \curl \vb H = \pdv{\vb D}{t} + \vb j \]

  • Wellenausbreitung \(\curl\vb H=\pdv*{\vb D}{t}\): Georadar
  • Induktion \(\curl\vb H=\vb j\), \(\curl\vb E=-\pdv*{\vb B}{t}\)
    Magnetotellurik, EMI in Zeit- und Frequenzbereich
  • Magnetostatik \(\curl\vb H=\vb j\), \(\curl\vb E=0\)

2.5 Näherung und Methoden

\[ \curl \vb H = \pdv{\vb D}{t} + \vb j \qquad \curl\vb E = -\pdv{\vb B}{t} \]

Vorgänge/Felder \(\curl\vb H\) \(\curl\vb E\) Methoden
Stationär \(\vb j\) 0 Geoelektrik
Quasistationär \(\vb j\) \(-\pdv*{\vb B}{t}\) Induktion
Wellenvorgänge \(\pdv*{\vb D}{t}+\vb j\) \(-\pdv*{\vb B}{t}\) Georadar

2.6 Elektrische Leitfähigkeit

Elektrischer Widerstand
  • Länge \(L\), Querschnitt \(A\)
  • Spannung \(U\) [V], Stromstärke \(I\) [A]
  • Ohmscher Widerstand \(R\)=\(U\)/\(I\) [\(\Omega\)]
  • Stromdichte \(j\)=\(I\)/\(A\) [A/m²]
  • Elektr. Feld \(E\)=\(U\)/\(L\) [V/m]
  • Ohmsches Gesetz
    \(j=\sigma E\) (Leitfähigkeit \(\sigma\) [A/Vm])
    \(\Rightarrow \rho=1/\sigma=E/j\) [\(\Omega\)m] spezifischer Widerstand

\(\Rightarrow \rho = \frac{U/L}{I/A}=\frac{U}{I}\frac{A}{L}=R\cdot k\) Konfigurations-Faktor \(k\) [m]

2.7 Spezifischer Widerstand (Knödel, 2005)

25 Dekaden:

  • Diamant: 1e18 \(\Omega\)m
  • Kupfer: 1e-7 \(\Omega\)m

abhängig von:

  • Wassergehalt
  • Fluidleitfähigkeit
  • Tonanteil

3 Anwendung

  • Ingenieurgeophysik (Festgestein, Hohlräume, Bäume/Wurzeln)
  • Hydrogeologie (Aquifer-Architektur, Kontaminationen)
  • Salzwasser-Intrusion in Grundwasser-Aquifere
  • Geologisch/tektonische Untersuchungen (Störungen)
  • Archäologie (Fundamente, Hütten und Schlacken)
  • Landwirtschaft (Feuchtigkeit und Tonanteil)
  • Georisiken (Hangrutschungen, Verkarstung)
  • Prozessmonitoring (Wasser- und Stoff-Transport)

3.1 Kleinskalig (Martin&Günther, 2013)

3.2 Mittelskalig: Geologie

Tanner et al. (2019): Störungszone Irland

3.3 Großskalig: Geologie

Nickschick et al. (2017): CO2-Aufstieg an Störung im Egergraben

3.4 Salzwasser-Intrusion von Süßwasseraquiferen

3.5 Hydrogeologie (Hübner et al., 2015)

Wasserspeicher in Hangdellen am Teststandort Mulda

3.6 Untersuchung von Halden (Günther et al., 2006)

3.7 Zielfunktion guter Leiter

  • 3-Schicht-Fall
  • Mächtigkeiten 20m
  • spez. Widerstände 100/1000/100 \(\Omega\)m

3.8 Sensitivitäten

3.9 Wenner-Anordnung

3.10 Schlumberger-Anordnung

3.11 Wenner-\(\beta\) Anordnung

3.12 Dipol-Dipol-Anordnung

3.13 Mess-Strategien

Kartierung

“Bewegung” einer festen Anordnung entlang eines Profils

Sondierung

Vergrößerung der Auslage bei festem Mittelpunkt

Sondierungskartierung/2D-Geoelektrik

Kombination aus Kartierung und Sondierung, meist bei festem Aufbau einer Elektrodenkette (Multielektroden-Geoelektrik): Electrical Resistivity Tomography (ERT)

3.14 Widerstands-Kartierung (Wenner)

3.15 Tiefensondierung (Schlumberger)

4 2D-Geoelektrik

  • Kombination aus Kartierung und Sondierung, automatisiert mit Relais
  • Entwicklung von Multi-Elektroden-Apparaturen (1990er Jahre)
  • Hersteller weltweit: ABEM (GeoGuideline, SE), IRIS Instruments (FR), SuperSting (US), Resecs, Lippmann, GeoTom (DE), ARES (CZ), Eigenbau-Systeme (OhmPi, FR)

4.1 Praktische Realisierung

  • Kombination aus positiver & negativer Phase mit Ausschaltpausen

Stromsignal (square wave with 50% duty cycle)

4.2 2D-Geoelektrik (Schlumberger)

4.3 2D-Geoelektrik (Wenner-\(\beta\))

4.4 Beispiel: Modell

4.5 Stromverteilung homogener Fall

4.6 Beispiel: Strom & Potential

4.7 Beispiel: Pseudosektion

4.8 Beispiel: Untergrundmodell

4.9 Pseudosektion Wenner-Array

4.10 Pseudosektion Schlumberger-Array

4.11 Dipol-Dipol mit leitfähiger Einlagerung

4.12 Dipol-Dipol mit resistiver Einlagerung

5 Das pyGIMLi-Modul ERT

  • Generierung von Datensätzen (Quadrupolanordnung)
  • Berechnung von Geometriefaktoren analytisch oder numerisch
  • Synthetische Simulation von Widerstandsverteilungen
  • Daten-Import, -Filterung, Fehlerbestimmung
  • 2D/3D-Inversion von Datensätzen
  • Auswertung von Frequenz- und Zeitbereichs-IP-Daten
  • Timelapse-Inversion von Monitoring-Daten

6 Das pyGIMLi-Modul VES

  • Beliebige Anordnungen, typisch Schlumberger (AB/2, MN/2)
  • Simulation beliebiger Schichtmodelle
  • Auswertung von n-Schicht-Modellen (Widerstände & Mächtigkeiten)
  • Auswertung glatter Verteilung (nur Widerstände)

6.1 Vergleich von Elektrodenanordnungen

WA WB SL DD MG
Signalstärke H M M L H
Horizontale Auflösung M M M H M
Vertikale Auflösung L L M L M
Erkundungstiefe H H M L L
Vertikale Sondierung L L H M M
Kartierung H H L H M

(Wenner-\(\alpha\)/\(\beta\), Schlumberger, Dipol-dipol, Multi-gradient)

6.2 Beispiel: Modell

6.3 Beispiel: Pseudosektion

6.4 Beispiel: Untergrundmodell

6.5 Inversion

6.6 Inversion - Methodik

Datenvektor \(\vb d=\{d_i\}=\{\rho^a_i\}\), Modellvektor \(\vb m=\{m_j\}=\{\rho_j\}\)

Datenmisfit \(\vb \delta \vb d = \vb d - \vb f(\vb m)\) soll klein sein (Quadrate)

Modellveränderung \(\delta\vb m\) über Sensitivitätsmatrix \(S_{i,j}=\pdv{f_i(\vb m)}{m_j}=\pdv{\rho_a(\rho)}{\rho}\)

\[\Rightarrow \vb S \vb \delta \vb m = \delta \vb d\]

Least-SquaresLösung \(\vb S^T \vb S \delta \vb m=\vb S^T \delta \vb d\) +Zusatzterme