Friedrich Robert Helmert definierte die Aufgaben der Geodäsie als die Lehre von der Ausmessung und Abbildung der Erdoberfläche. Diese Definition trifft auch heute noch auf einen Teilbereich wie bspw. die Landesvermessung zu, grenzt aber gleichzeitig auch Bereiche aus. So finden sich einige Aufgaben, die von Vermessungsingenieuren ausgeübt werden, in dieser Definition nicht wieder. Beispielhaft seien hierbei Tätigkeitsfelder der präzisen Ingenieurvermessung und Messtechnik wie die Qualitätsüberwachung und -sicherung, Konvergenzmessungen oder die Deformationsanalyse genannt. Neben der Sicherheit für Mensch und Maschinen machen auch finanzielle Aspekte diese Aufgaben so wichtig. Aus diesem Grund soll die Deformationsanalyse in diesem Beitrag näher betrachtet und an einem Beispiel mit dem Ausgleichungsprogramm JAG3D, welches seit der Version v3.1 über eine entsprechende Option verfügt, durchgeführt werden. In einem zweiten Teil dieses Artikels sollen weitere Beispiele aus der Fachliteratur näher betrachtet werden.

Grundsätzlich existieren verschiedene Berechnungsmöglichkeiten bei der Durchführung einer Deformationsanalyse. Allen gemein ist, dass sie signifikante Veränderungen in den zu überprüfenden Objektpunkten X_O aufzeigen und instabile Referenz- bzw. Stabilpunkte X_R lokalisieren sollen. Eine häufig anzutreffende Variante ist die Koordinatenmethode (Jäger et al., 2005; Niemeier, 2008). Hierbei liegen die Koordinaten der Punkte des Deformationsnetzes inklusive der vollbesetzten Varianz-Kovarianz-Matrix C aus zwei (oder mehr) zu vergleichenden Epochen vor, die jeweils aus einer freien Netzausgleichung mit identischen Datum gewonnen wurden. Das funktionale Modell der gemeinsamen Ausgleichung lautet für zwei Epochen unter Verwendung von Einheitsmatrizen I:

Koordinatenmethode

Die als stabil angenommenen Referenzpunkte werden aus beiden Epochen zusammengefasst, während die Objektpunkte pro Epoche einmal eingeführt werden. Das stochastische Modell setzt sich aus den beiden Varianz-Kovarianz-Matrizen C₁₁ bzw. C₂₂ der betrachteten Epochen zusammen:

VKM zweier unkorrelierter Epochen

Mathematische Abhängigkeiten zwischen den betrachteten Epochen werden hierbei i.d.R. nicht berücksichtig, sodass die Nebendiagonalelemente C₁₂ bzw. C₂₁ Null sind. Veränderungen in den Objektpunkten können unter der Voraussetzung, dass die Referenzpunkte stabil geblieben sind, direkt abgeleitet werden.

Verschiebungsvektor

Mittels Hypothesentest lässt sich abschließend klären, ob die ermittelten Punktveränderungen signifikant oder eher dem Messrauschen zuzuordnen sind. Die auf die Fischer-Verteilung bezogene Testgröße T ergibt sich aus:

Fischer-verteilte Testgröße

mit

Gewichtsmatrix des Verschiebungsvektors

Veränderungen der als stabil angenommenen Referenzpunkte lassen sich durch einen Globaltest aufdecken. Da jeder Referenzpunkt für das ablehnen der Nullhypothese verantwortlich sein kann, sind die Referenzpunkte sukzessiv auf Veränderung zu prüfen (Niemeier, 2008). Neben dieser impliziten Suche nach instabilen Referenzpunkten schlagen Jäger et al. (2005) eine effektivere explizite Methode vor, die der Vorgehensweise des verallgemeinerten Datasnopping entspricht. Für nähere Details sei der interessierte Leser auf die zitierten Quellen verwiesen.

Neben der koordinatenbasierenden Vorgehensweise gibt es noch die Beobachtungsmethode, die auch im Überwachungssystem GOCA Anwendung findet. Statt die Koordinaten und Varianz-Kovarianz-Matrizen separater freier Netzausgleichungen zu betrachten, werden hierbei die originären Messwerte zweier (oder mehrerer) Epochen in einer gemeinsamen Ausgleichung ausgewertet. Das funktionale Modell für k Epochen lautet hierbei:

Funktionales Modell der Beobachtungsmethode

Anhand der Struktur der Designmatrix A lässt sich erkennen, dass in allen Epochen identische Referenzpunkte verwendet werden und somit ein identisches Datum zugrunde liegt. Die identischen Objektpunkte werden hingegen einmal pro Epoche als Neupunkte in die Ausgleichung eingeführt. Dies entspricht auch der Vorgehensweise der o.g. gemeinsamen Ausgleichung nach der koordinatenbasierten Methode. Eine Prüfung der Stabilpunkte kann wiederum implizit und explizit erfolgen (Jäger et al., 2005). Die vollbesetzte Varianz-Kovarianz-Matrix der Netzausgleichung ist abschließend wiederum für den Hypothesentest zum Aufdecken signifikanter Verschiebungen heranzuziehen. Im Gegensatz zur Koordinatenmethode sind hierbei die Nebendiagonalelemente C₁₂ bzw. C₂₁ nicht Null. Es existieren somit Korrelationen zwischen den Epochen. Da die Beobachtungsmethode auf einer reinen Netzausgleichung aufbaut, bietet sie sich für eine Integration in bestehende Ausgleichungsprogramme an und wurde in der freien Netzausgleichungssoftware JAG3D implementiert.

Ein fiktives 2D-Beispiel einer Überwachungsmessung an einer Staumauer soll eine mögliche Vorgehensweise mit JAG3D demonstrieren. Hierzu liegen die terrestrischen Beobachtungen der Null- und einer Folgeepoche vor. Das Netz besteht aus elf Objektpunkten, die sich auf der Mauerkrone befinden. Auf der Luftseite wurden sechs Vermessungspfeiler geschaffen, von denen die Punkte auf der Staumauer aus beobachtet werden können. Diese sechs Pfeiler bilden die Referenzpunkte des Deformationsnetzes.

Skizze des Deformationsnetzes

Das Beobachtungskonzept sieht vor, dass alle Pfeiler einmal besetzt werden. Von jedem Pfeiler aus werden Richtungs- und Streckenmessungen zu den übrigen durchgeführt. Von den Pfeilern 1001, 1002 und 1006 werden darüber hinaus auch die Objektpunkte mitbestimmt.

Das in der jeweiligen Messepoche erhobenen Datenmaterial ist vor der eigentlichen Deformationsanalyse auf fehlerbehaftete Beobachtungen zu prüfen. Ferner ist bei der epochenweisen Auswertung das stochastische Modell der Beobachtungen abzustimmen, da sich dieses Aufgrund von unterschiedlichen Instrumenten zwischen den Epochen ändern könnte. Aus diesen Gründen erfolgt die getrennte Auswertung der beiden Einzelepochen jeweils durch eine freie Netzausgleichung. Durch die freie Netzausgleichung ist gewährleistet, dass kein zusätzlicher Zwang auf das Netz ausgeübt wird und die Netzstruktur allein aus dem erhobenen Beobachtungsmaterial resultiert. Das geodätische Datum der freien Netzausgleichung wird durch die als stabil angenommenen Messpfeiler 1001-1006 definiert. Die Objektpunkte 100-110 werden als Neupunkte in die Ausgleichung eingeführt.

Mittels der Varianzkomponentenschätzung kann die Abstimmung des stochastischen Modells erfolgen. Die Gruppenvarianzfaktoren sollten in etwa zwischen 0.6 und 1.4 liegen (Benning, 2010). Es zeigt sich, dass sich bei den horizontalen Strecken ein Genauigkeitsansatz von 2mm und 2ppm anbietet. Bei den Richtungen ergeben sich ein absoluter Anteil von 0.2mgon und eine Zentriergenauigkeit von 0.3mm für beide Epochen. Grobe Fehler lassen sich in beiden Epochen nicht nachweisen.

Werden die ausgeglichenen Koordinaten der Nullepoche als erste Näherung in den Folgeepochen eingeführt, so kann eine erste Abschätzung möglicher Veränderungen anhand der ermittelten Verbesserungen abgeleitet werden, da das Netz der Folgeepoche so direkt auf die Koordinaten der Nullepoche transformiert wird.

Die Analyse der Restklassen erlaubt erste Rückschlüsse auf mögliche Deformationen. Die Objektpunkte weisen fast alle Veränderungen in der Größenordnung von 1-2cm auf. Auch bei den als stabil angenommenen Punkten sind Auffälligkeiten zu erkennen. Mindestens die Punkte 1004 und 1005 sollten in der sich anschließenden Deformationsanalyse kritisch betrachtet werden.

Um eine Deformationsanalyse innerhalb von JAG3D durchführen zu können, sind die Daten beider Epochen in einem Projekt zu vereinen. Wie oben bereits beschrieben, müssen die Objektpunkte pro Epoche separat eingeführt werden. Da die Punktnummer innerhalb von JAG3D einzigartig sein muss, sind die Punktenummern somit zwischen den Epochen umzubenennen und das Beobachtungsmaterial entsprechend der Punktnummern anzupassen. In dem Beispiel hier haben die Objektpunkte in der Nullepoche 100er und in der Folgeepoche 200er Nummern. Ein einheitliches Datum beider Epochen wird dadurch erreicht, dass die Referenzpunkte in beiden Epochen gleich sind, wie aus dem funktionalen Modell der Beobachtungsmethode ersichtlich wird.

Um die als unverändert angenommenen Referenzpunkte überhaupt auf signifikante Verschiebungen testen zu können, muss dem Programmsystem noch mitgeteilt werden, welche Beobachtungsgruppen zur Null- und welche zur Folgeepoche gehören. Hierzu bietet JAG3D im Eigenschaftendialog der jeweiligen Beobachtungsgruppe eine entsprechende Option an.

JAG3D - Gruppendialog der Beobachtungen

Beobachtungsgruppen, die zur Nullepoche gehören, sind der Referenzepoche zugeordnet, und die Beobachtungen der Kontrollmessung der Folgeepoche.

Weiterhin sind die auf Veränderung hin zu untersuchenden Punkte bzw. Vektoren d im Projekt zu definieren, indem die Punkte der Referenzepoche denen der Folgeepoche zugeordnet werden. Abschließend ist in den Projekteinstellungen lediglich die Option Deformationsanalyse durchführen auszuwählen und die Berechnung zu starten.

JAG3D - Definition der Verschiebungsvektoren

Bevor die Objektpunkte auf mögliche Änderungen untersucht werden, ist sicherzustellen, dass die Referenzpunkte stabil geblieben sind. Die Voranalyse durch Vergleichen der Ergebnisse der Einzelepochen ergab bereits, dass möglicherweise einige Referenzpunkte zwischen den Epochen nicht unverändert geblieben sind. Diese Vermutung wird durch den durchgeführten Stabilpunkttest innerhalb der Deformationsanalyse bestätigt. Eine Verschiebung des Punktes 1004 zwischen den Epochen ist hochsignifikant. Aufgrund der bekannten Verschmierungseffekte innerhalb der Ausgleichungsrechnung werden auch die anderen Referenzpunkte als verändert ausgewiesen. Ob weitere Punkte instabil sind, muss sukzessiv geprüft werden. In der Annahme, dass sich nur ein Punkt verändert hat, wird zunächst nur der Punkt 1004 als verschoben verdächtigt, da dieser mit Abstand die größten Testgrößen erhalten hat.

JAG3D - Ergebnis des Stabilpunkttests

Grundsätzlich könnte der Punkt 1004 aus der Ausgleichung gestrichen werden, indem er deaktiviert oder gelöscht wird. Da aber die Einzelepochenauswertung gezeigt hat, dass keine Unstimmigkeiten innerhalb des Beobachtungsmaterials vorliegen, hätte ein Ausschluss des Punktes 1004 den unnötigen Verzicht aller Beobachtungen, an denen der Punkt beteiligt ist, zur Folge. Sinnvoller erscheint es daher, den Punkt 1004 in die Gruppe der instabilen Objektpunkte zu schieben und einen Verschiebungsvektor zu definieren. Der Punkt 1004 wird in der Folgeepoche mit der Punktnummer 2004 eingeführt und die Berechnung wiederholt. Es zeigt sich, dass anschließend keiner der verbleibenden fünf Referenzpunkte als verschoben ausgewiesen wird, sodass nun eine Betrachtung der Ergebnisse der Objektpunkte erfolgen kann. Hier zeigt sich, dass bis auf die Randpunkte 100 und 110 der Staumauer alle Punkte signifikant verschoben sind. Durch die geschickte Festlegung des Koordinatensystems, bei dem die X-Achse parallel zur Staumauer verläuft, lässt sich auch die Verschiebungsrichtung leicht ableiten. Während in X-Richtung quasi keine Änderungen stattgefunden haben, sind deutliche Deformationen der Mauer in Y-Richtung zu erkennen. Diese sind zudem einseitig und verlaufen von der Wasserseite in Richtung der Luftseite.

Skizze des Deformationsnetzes mit Verschiebungsvektoren und Konfidenzellipsen

Auch der Deformationsvektor 1004-2004 weist deutliche Verschiebungsbeträge auf. Insgesamt ergibt sich eine Verschiebung von über 4cm für den Messpfeiler, der damit als Referenzpunkt in weiteren Epochen nicht mehr genutzt werden kann.

JAG3D - Veränderungen der Objektpunkte

Deformationsmessungen und -analysen liefern Hinweise auf Verformungen bzw. Veränderungen eines Körpers und beantwortet Fragen z.B. zur Standsicherheit. Verschiedene Modelle zum Durchführen dieser Aufgabe sind bekannt und in der geodätischen Literatur beschrieben, sodass sie in diesem Beitrag nur kurz angeschnitten wurden. Besonders die in Jäger et al. (2005) publizierte Beobachtungsmethode lässt sich leicht in bestehende Netzausgleichungsprogramme integrieren und wurde im freien Programmpaket JAG3D implementiert. Am Beispiel einer fiktiven Staumauer wurde eine Vorgehensweise bei der Auswertung von Deformationsmessungen kurz umrissen. Alle in diesem Testnetz eingebrachten Deformationen konnten dabei von JAG3D aufgedeckt werden. Auch bei dem in Jäger et al. (2005) abgedruckten Beispiel eines Höhendeformationsnetzes konnten übereinstimmende Ergebnisse abgeleitet werden.

Die für dieses Beispiel simulierten Daten stehen zum Download zur Verfügung, sodass die Möglichkeit besteht, die abgeleiteten Ergebnisse mit anderen Produkten zu vergleichen.

Beispieldatensatz zur Deformationsanalyse

22.09.2010 von Michael Lösler