Potential der einfachen Schicht


Potential der einfachen Schicht

Als Potential der einfachen Schicht wird eine von Karl-Rudolf Koch in den 1970er-Jahren entwickelte Methode bezeichnet, mit der durch Einführung von Flächenbelegungen die Berechnung des Erdschwerefeldes vereinfacht oder beliebig verfeinert werden kann. Die Methode hat zwar – im Gegensatz etwa zu den harmonischen Modellen mittels Legendre- und Kugelfunktionen – bei der Approximation des Schwerepotentials eine unstetige Charakteristik an den Modellrändern, benötigt aber nur einen Bruchteil der Rechenzeit und ein etwa halb so großes Gleichungssystem wie die klassische Neumann'sche Methode.

Inhaltsverzeichnis

Geoid, Störpotential und Flächenbelegungen

Das Schwerepotential und seine Funktionale – die wichtigsten sind das Geoid, die Lotabweichungen und die Schwereanomalien) – werden großteils durch die Gravitation des Erdkörpers und die Fliehkraft der Erdrotation verursacht. Jedoch bewirken alle Unregelmäßigkeiten der Erdoberfläche (Gelände, unterschiedliche Dichten) und der geologische Aufbau der Erdkruste zahlreiche Abweichungen vom mathematischen „Normalfeld“, das dem theoretischen Schwerefeld eines mittleren Erdellipoids gleichgesetzt wird.

Jede zusätzliche oder gegenüber dem Erdellipsoid fehlende Masse – die in ihrer Gesamtheit als Massenverteilung bezeichnet wird – verändert das Schwerepotential der Erde, und zwar umso mehr, je näher die Anomalie am jeweiligen Messpunkt liegt. Die Potentialänderung wird als Störpotential bezeichnet und kann durch ein Skalar (eine lokale Differenz in der potentiellen Energie) ausgedrückt werden.

Ein Gebirgsmassiv erhöht z. B. das Erdpotential um einige Zehntel Promille, was sich in erhöhter Schwerkraft und einer geringfügigen Auswölbung der Niveauflächen und des Geoids auswirkt. Besser vorstellbar ist dieser Effekt, wenn man die Lotrichtung betrachtet (siehe nebenstehende Skizze):

Das Gebirge zieht eine frei hängende Lotschnur etwas zu sich, wodurch sich auch die darauf senkrechte örtliche Horizontale (eine Parallele zum Geoid) etwas verändert. Die resultierende Aufwölbung wird Geoidundulation genannt. Auf der anderen Seite des Gebirges verläuft der Effekt entgegengesetzt, sodass in größerer Entfernung davon nur mehr wenig überbleibt. Bei den entsprechenden Berechnungen, die mit Gesteinswürfeln oder Prismen durchgeführt werden, muss zwar die nächste Umgebung genau berücksichtigt werden, für das fernere Gelände genügt aber ein relativ grobes Modell. Die Lotabweichungen betragen in den Alpen maximal 50″ (0,015°), im Hügelland etwa ein Zehntel davon, sie sind aber bei jeder genauen Vermessung zu berücksichtigen (siehe topografische Reduktion).

Die Methode der Flächenbelegung kann diese Berechnungen vereinfachen, indem die Berge durch Massenbelegungen modelliert werden. Es sind dies unendlich dünne, übereinander legbare Platten, deren (fiktive) Masse dem Gesteinskörper entsprechen. Ein Massendefizit (z. B. ein tiefes Tal) wird durch negative Massen modelliert.

Das Störpotential und seine Funktionale

Zur Berechnung des Störpotentials T an einem bestimmten Punkt sind theoretisch alle Störmassen (d. h. die Abweichungen der Erdfigur von einem idealen Ellipsoid) ins Kalkül zu ziehen. De facto ist aber nur für die nächste Umgebung ein genaues digitales Geländemodell erforderlich, während für Entfernungen über 50 km ein sehr grobmaschiges Modell ausreicht.

Das gesamte Schwerepotential W an einem Punkt P der Erdoberfläche (der sog. Aufpunkt) mit den kartesischen Koordinaten X, Y, Z lässt sich formal als Volumenintegral über die gesamte Erdmasse anschreiben, indem das auf P wirkende Potential aller Massenpunkte der Erde summiert wird. Diese Massenpunkte mit dem Volumenelement dV haben die individuelle Dichte ρ (Gesteine 2.5–3.3 g/cm³, Erdmantel 4–6 g/cm, Erdkern ~10 g/cm³), die von ihrer Lage im Erdkörper (x, y, z) abhängt. Im Nenner steht der Vektor r zwischen dem Aufpunkt P und dem jeweiligen Massenpunkt ρ.dV.

Numerisch ist eine solche Lösung natürlich unmöglich, weil die Erde schon bei Zerlegung in 1 km³ große „Punktmassen“ in einer Billion Teile modelliert werden müsste. Außerdem kennt man den Verlauf der Dichte im tieferen Untergrund nicht genau genug. Die praktische Berechnung solcher Potentiale muss sich daher mit Näherungslösungen und gebietsweisen Abschätzungen begnügen.

Eine wesentliche Vereinfachung ergibt sich, wenn die obige formale Gleichung nur aus die sog. Störmassen beschränkt wird, welche die Abweichung vom Erdellipsoid repräsentieren. Das Erdinnere wird dabei mit Theorien wie der Gleichgewichtsfiguren erfasst, was heute auf wenige Millionstel genau möglich ist. Das zugehörige theoretische Potential wird mit U bezeichnet

Das Störpotential T ergibt sich somit formal zu … …

was bereits einer annähernd praktikablen Modellierung zugänglich ist.

Die weiteren interessierenden Größen des Schwerefeldes sind Funktionale dieses Störpotentials T, wobei γ die elliposoidische Normalschwere am Aufpunkt P ist und R der mittlere Erdradius. Das kartesische Koordinatensystem x,y,z wird durch ein lokales System u,v,w ersetzt, worin w in Richtung der örtlichen Vertikale weist, u nach Norden und v nach Osten:

Geoidundulation ζ = T/γ
Schwereanomalie Δg = \frac{\partial T}{\partial w} - T.\frac{2}{R}
Komponenten der ξ = - \frac{1}{γ} \frac{\partial T}{\partial u}
Lotabweichung η = - \frac{1}{γ} \frac{\partial T}{\partial v}

Approximation des Störpotentials mit Flächenbelegungen

Die Modellierung des Erdschwerefeldes erfolgt mittels Potential der einfachen Schicht (potential of a simple layer). Mehrere dieser dünnen, mit konkreten Massen behafteten Schichten werden fiktiv auf der Erdoberfläche ausgebreitet und können sich allenfalls überlagern. Ihre Dichten werden als Unbekannte angesetzt und mittels der gegebenen Schweredaten durch Ausgleichsrechnung nach kleinsten Quadraten ermittelt.

In der ersten Version der Methode (Koch 1970) wurden 192 Oberflächenelemente definiert, deren Massen einem harmonischen Potentialmodell (Kugelfunktionen bis zu Grad und Ordnung 15) der Satellitengeodäsie angepasst wurden, sowie einem großen Datensatz von terrestrischen Schwereanomalien.

Die Dichte dieser Kugelkalotten bezogen sich auf ein Referenzellipsoid, das dieselbe Abplattung hat wie ein Erdkörper im hydrostatischen Gleichgewicht. Die derart ermittelten Ergebnisse haben daher auch Bezug zu geophysikalischen Fragestellungen.

In späteren Anwendungen (Koch 1975f) wurden zusätzliche Kombinationslösungen zwischen Schwereanomalien und Satellitenaltimetrie durchgeführt. Diese Höhenmessungen waren bereits Mitte der siebziger Jahre so genau, dass eine gemeinsame Modellierung mit den Bahnstörungen versucht werden konnte.

Die zu bestimmenden Dichtewerte der Oberflächenelemente werden zwar als konstant angesetzt, doch lässt sich das Auflösungsvermögen des Modells durch eine Änderung ihrer Anzahl an die Qualität der Schweredaten anpassen.

Wenn die Altimeter- oder Schweredaten die Erde sehr dicht überdecken bzw. von hoher Genauigkeit sind, kann das Verfahren durch einen feineren Raster von Potentialschichten sehr flexibel gestaltet werden.

Bei der Berechnung dieser Oberflächenelemente kann die Ausgleichung durch zusätzliche Übergangs- oder Pufferzonen ergänzt werden und lässt sich zur Lösung sehr großer Gleichungssysteme in kleinere, unabhängige Subsysteme aufspalten.
Auch ein Übergang auf Kollokationsmethoden ist möglich.

Literatur und Quellen


Wikimedia Foundation.

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Paradoxon der schwachen jungen Sonne — Das Paradoxon der schwachen jungen Sonne bezeichnet den Widerspruch zwischen der in den vergangenen 3,8 Mrd. Jahren deutlich geringeren Strahlungsleistung der jungen Sonne und eindeutigen Hinweisen auf ein wärmeres Klima in der frühen… …   Deutsch Wikipedia

  • Potenzialtheorie — Die Potentialtheorie behandelt die mathematisch physikalischen Grundlagen von konservativen Kraftfeldern. Man nennt sie auch die Theorie der wirbelfreien Vektorfelder. Wichtige Anwendungen sind das elektrische Feld, das Magnetfeld und das… …   Deutsch Wikipedia

  • Potentialtheorie — Die Potentialtheorie oder Theorie der wirbelfreien Vektorfelder behandelt die mathematisch physikalischen Grundlagen konservativer Kraftfelder. Wichtige Anwendungen sind das elektrische Feld, das Magnetfeld und das Gravitations bzw. Schwerefeld.… …   Deutsch Wikipedia

  • Karl Rudolf Koch — (* 30. Juli 1935 in Hilchenbach, Kreis Siegen) ist ein deutscher Geodät und Professor an der Universität Bonn. Seit Juli 2000 ist er emeritiert; als sein Nachfolger wurde Wolf Dieter Schuh von der TU Graz berufen. In Fachkreisen ist er vor allem… …   Deutsch Wikipedia

  • Geländereduktion — Als topografische Reduktion wird in der Geodäsie und in der Geophysik die rechnerische Beseitigung der Topografie (des Geländes) bezeichnet. Diese Reduktion ist bei verschiedensten Aufgaben notwendig, um die vom Gelände verursachten Störungen im… …   Deutsch Wikipedia

  • Flächenbelegung — Als Flächenbelegung wird in der Potentialtheorie die theoretische oder rechnerische Applikation von Masse oder Ladungen an Modelle von Festkörpern bezeichnet, um die Wirkung von kleinen Änderungen mit einfachen Mitteln zu berücksichtigen oder… …   Deutsch Wikipedia

  • Karl-Rudolf Koch — (* 30. Juli 1935 in Hilchenbach, Kreis Siegen) ist ein deutscher Geodät und Professor an der Universität Bonn. Seit Juli 2000 ist er emeritiert; als sein Nachfolger wurde Wolf Dieter Schuh von der TU Graz berufen. In Fachkreisen ist er vor allem… …   Deutsch Wikipedia

  • Geoidbestimmung — Mit dem Begriff Geoidbestimmung bezeichnet die Geodäsie Verfahren, die zur Definition des Geoids führen. Diese Aufgabe ist eng mit der Bestimmung von Details des Erdschwerefeldes verbunden. Die Methoden zur Messung der benötigten physikalischen… …   Deutsch Wikipedia

  • Christian Tscherning — Carl Christian Tscherning (* 1942 in Kopenhagen) ist ein dänischer Geodät und Hochschullehrer. Er gilt in der Geodäsie als Wegbereiter für moderne mathematische Methoden, insbesondere für die genaue Bestimmung des irdischen Schwerefeldes.… …   Deutsch Wikipedia

  • Tesseral — Als tesseral werden mathematische Funktionen bezeichnet, mit denen man Phänomene auf einer Fläche beschreiben kann, die einen schachbrett artigen Verlauf besitzen. Kugelflächenfunktionen: oben zonal, Mitte sektoriell, unten tesseral. In den… …   Deutsch Wikipedia


Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.