Integrating GIS methods for the analysis of GeoSystems

Abstract

Modelling of hydrosystems is of increasing demand in fields of water and environmental research and basically includes modelling of surface water and groundwater flow. Due to the limitation of available data and the complexity of geo-systems themselves in most modelling tasks, Geographical Information System (GIS) is taken as a useful tool due to the ability to store diverse kinds of hydrogeological information in a common understandable structure combined with its capabilities. The GIS data are served as the input parameters for the numerical modelling. Modern GIS , not yet based on a real 3D spatial representation, focus on 2 and 2.5 dimensions to represent surface and groundwater system. The construction of 3D model is challenging and of vital importance in order to accurately simulate hydrodynamic flow. This work presents the integration implementation of GIS methods within the numerical modelling software GeoSys/Rockflow (Kolditz et al. 2004) and demonstrates how GIS technology and its integration make the geometric construction of 2D and 3D hydrosystem models easier. The 3D structural model was built based on the implemented integration techniques using existing GIS data. The GIS project contains hydrogeological boundaries extracted from geological maps, Digital Elevation Model (DEM), raster maps for the aquifer layers, fault system, river network, lithology (cross section), boreholes, and springs. The data are kept in a GIS project and converted into geometric components of the modeling system : GEOLib, and then can be used as boundary conditions, initial conditions, source sink terms, material groups and model domain. Fundamental geometric objects are points, polylines, surfaces, volumes and domains. In addition, the concept of advanced layer models, such as layer polylines, layer surfaces with Triangulated Irregular Network (TIN), layer volumes are implemented in the finite element program, which are the methodology for the construction of a 3D structural model. In order to accomplish these aims described above, the following contributions were achieved: To design and implement the geometric objects library (GEOLib) within GeoSys/RockFlow; to create the interface between GIS data and GeoSys/RockFlow; to extend the existing Graphical User Interface (GUI) of GeoSys/RockFlow for the features mentioned. As a result, the GIS integrated concept was demonstrated in several case studies in Middle East, Europe, and in local area Tübingen, which was also used as an education example for the students.

Die Nachfrage nach Simulationen nimmt in der Hydrologie und Umweltforschung zu. Diese Simulationen beinhalten die Abbildung von Oberflächen- und Grundwasserströmungen. Daten über die betrachteten Systeme stehen nur in beschränktem Umfang zur Verfügung. Außerdem weisen die Systeme eine hohe Komplexität auf. Verschiedene hydrogeologische Daten können in GIS verwaltet und dargestellt werden und werden daher als nützliches Modellierwerkzeug herangezogen. Die Daten aus einem GIS Projekt dienen dann als Eingabeparameter zur numerischen Simulation. Das heutige GIS kann Strukturen noch nicht in ihrer eigentlichen 3-dimensionalen Gestalt bearbeiten. Stattdessen stellt es die Geometrie von Oberflächen und Grundwassersystemen in 2 und 2,5 Dimensionen dar. Die Generierung 3-dimensionaler Modelle ist jedoch erforderlich, um hinreichend genaue Ergebnisse bei der Modellierung von Hydrosystemen zu erhalten. Die vorliegende Arbeit behandelt die Einbindung des GIS in das numerische Softwarepaket GeoSys/Rockflow (Kolditz et al. 2004). Gezeigt wird auch, inwiefern somit die Erstellung 2- und 3-dimensionaler numerischer hydrologischer Modelle vereinfacht wird. Ein 3-dimensionales Modell wurde exemplarisch mithilfe einese GIS-Datensatzes erstellt. Der Datensatz beinhaltet hydrogeologische Ränder, digitale Höhenmodelle (DEM), Rasterkarten, Kluftsysteme, Flussnetze, Lithologien (Querschnitte), Bohrlöcher und Quellen. Diese GIS-Daten werden als geometrische Objekte in der geometrischen Programmbibliothek GEOLib gespeichert. Fundamentale geometrische Objekte sind dabei Punkte, Polylinien, Oberflächen, Volumen und Gebiete. Diese Objekte werden zur Definition der Modellgeometrie, Anfangs- und Randbedingungen, Quelltermen verwendet. Darüber hinaus wird ein Konzept für Schichtenmodellierung aus Schicht-Polylinien, Schicht-Oberflächen (TINs) und Schicht-Volumen ins Finite Elemente Programm implementiert. Um die genannten Ziele zu erreichen, wurde die geometrische Bibliothek GEOLib entwickelt und wurde das bestehende Graphical User Interface (GUI) erweitert und verbessert. Die Neuerungen wurden anhand einiger Fallbeispiele entwickelt und getestet. Diese Fallstudien betrachten Gebiete im Mittleren Osten, Europa und dem Tübinger Umland. Letztere wurden auch zur Lehre herangezogen.