The Preferences menu allows to define the general configuration which is saved :
in the .thyrsis.ini file of the project folder where the computation database is saved,
in the thyrsis.ini file of the $HOME/.thyrsis folder.
At the first use of the plugin, the $HOME/.thyrsis folder is created and the thyrsis.ini file is set with default values. The Preferences dialog is opened and the user is asked to set his configuration (cf. infra).
When a project is re-opened, the configuration loaded is the .thyrsis.ini file of the database folder. Therefore it is possible to define a configuration for each project.
The variables preferences provide the opportunity to define which variables should be saved :
potential,
Darcy’s velocity,
concentration,
activity,
saturation,
and the display configuration, by defining:
the unit,
the scale,
the number of classes.
The scale defines the ratio between the minimal value to display and the maximal value obtained from simulation. By default the scale is set to 0.001, but a lower value can be chosen to get a better view of the plume extent.
The number of classes is used to display the plume with a graduated scale.
This preferences window enables also to define the mass unit for the mass balance and to check a box avoiding to display units, for communication purposes.
This window allows the user to connect to an experimental database by defining a connection service.
The push button Synchronize experimental database with site database enables the loading of boreholes, points of interest, hydraulic heads and volumic concentrations from the experimental base to the sites.sqlite database.
The creation of a new project can be made through the menu, by selecting an existing site or by creating a mesh for a new model (Figure 40) :
Figure 40 New menu allowing to select a site or build a new mesh
When selecting an existing site, a new drop-down menu is opened, which shows the available models and the available saved studies (defined by the simulations table in the sites.sqlite database) (Figure 41).
When selecting the entry Mesh…, a mesh project is created and a new model can be created (see Model building).
Figure 41 Open menu for selecting a study or a model
The “Open” item in the THYRSIS menu enables to open an existing spatialite computation database calcul.sqlite by selecting it in a file dialog.
If a calcul.qgs project exists in the same folder as calcul.sqlite, an alert dialog is opened, asking the user if the project has to be overwrite :
if the answer is “Yes”, a new calcul.qgs project file is built using the calcul.sqlite database and the existing project is overwritten;
if the answer is “No”, no action is done.
To open an existing calcul.qgs file - including the calcul.sqlite file - the Project>Open menu of QGIS should be used.
When a THYRSIS project is opened, several actions are added to the THYRSIS menu as shown in Figure 42. These actions are detailed below.
After computation, toolbars are updated (dates, variables, units), the “results” layer is added to the QGIS canvas and set on the last simulation date, and the vertical columns of the unsaturated zone are displayed in the ZNS window.
A simulation creates a new subfolder calcul_tmp in the computation database calcul.sqlite folder. This folder contains all input/output files of the hydrogeological computations. These files are managed in subfolders as follows :
samplexxxxx folders, where xxxxx is the id number of the parameters set (cf. table parametres_simulation),
in each of these folders :
insatyyyyy folders, where yyyyyy is the injection area number. These folders contains inpout/output files of each computation in unsaturated area,
a sature folder, which contains input/output files for computation in saturated zone.
Based on these results, THYRSIS generates numpy files in the computation database folder and these numpy files are then displayed in the interface. These files are :
calcul.variable.type.npy : values of the variable “variable” for each requested date (saturated zone), at each node if type=node or at each center of mesh if type=elem,
calcul.variable.insat.npy : values of the variable “variable” at each node of the ZNS column, for each requested date (unsaturated zone).
The variable is one of the following (according to the parameters chosen in Preferences) :
in saturated zone :
concentration,
potential,
Darcy’s velocity,
flow rate
in unsaturated zone :
saturation,
Darcy’s velocity,
concentration,
activity
The calcul_tmp folder can be deleted without consequence for Qgis interface use, exception made for probability map computation, which needs to access to the whole results for each parameters set.
This folder is overwritten at each “Simulate” action.
All input files needed by the hydrogeological software are automatically generated. But some particular cases need to use specific files, called “ user files ”. These files have a “u_” prefix and are placed at the same level that the computation database. They are described in the User files appendix.
The menu Thyrsis>Create provides the following actions (Figure 44) :
Injection zones
This action removes the existing injections and automatically creates new injection zones covering the whole surface of the area study. THYRSIS groups ZNS columns by thickness and permeability, according to the pas_echant_zns and rapport_max_permeability values, defined in the simulations table (Diffuse injection on the whole domain). Cf. example Figure 45.
Figure 45 Example of a study area covered by injection zones.
Latin Hypercube
This action allows to create several parameter sets with the Latin Hypercube method. Requested informations are (Figure 46) :
the number of simulations - defining the number of parameter sets to create,
and for each parameter (from the “parametres” table) :
the min value,
the max value,
the probability law to use, chosen betwen “constant”, “uniform”, “loguniform”, “normal”, “lognormal”.
Parameter sets are created in the parametres_simulation table.
The simulate action will generate as many simulations as there are parameter sets. These simulations are parallelized on the available processors.
After computation, THYRSIS creates two result files with the mean and standard deviation for each node or element and each date. Probability maps are available through the menu Thyrsis>Create>Probability map
Probability map action creates a probability layer based on several simulations. Threshold and unit should be specified in a pop-up window (Figure 47).
The result is added as a new result layer, as concentration or potential. In consequence the probability map is available through the combobox of the THYRSIS menu (Menu bar).
Isovalues action generates isovalues polylines on the active layer which can be the result layer or any of the raster layers describing the model properties (cf. 2D Panel) as long as the layer is connected to a legend.
Isovalues definition is specified in a pop-up dialog (Figure 48). By default, the values are extracted from the legend. These values can be manually modified/removed/added. Values can also be generated automatically, by providing the unit, the number of isovalues, the min and the max values, and the scale type (linear, log). The creation of values is triggered by the following button .
Isovalues (gdal) action generates isovalues polylines on the active layer, as the previous action, but this action uses the gdal_contour tools which needs to create an intermediate raster. Resolution can be provided in pixel/meter (default value : 0.1) (Figure 49). This method avoids some problems encountered with the previous action.
This action opens the 1D panel which displays evolution graphs of the selected variable for each selected point. This action can also be triggered using the button of the Menu bar.
Second medium
This action duplicates the current computation base in a new one, adapted to the Dual porosity model :
the tables noeuds_second_milieu, mailles_second_milieu, simulations_second_milieu and injections_second_milieu are copied from the tables noeuds, mailles, simulations et injections;
ASHAPE, BSHAPE, CSHAPE and PL entries are added to the parameters table
the parametres_simulation table is reset.
To exploit this new dual-porosity computation base, it is necessary to open a new projet (Thyrsis>Open) based on this database.
Chemical element
This action allows to create or remove chemical elements.
The menu Thyrsis>Export provides the following actions :
Raster
The Raster action exports a PNG or TIF image of the current 2D layer, with the current date and the legend.
Animation
The Animation action generates several PNG images on a simulation period, including the 2D display (Qgis canvas) but also the ZNS and the 1D panel. The record starts by clicking on the button of the Menu bar, and stops by clicking on the same button, became .
Video
The Video action generates an AVI video including the 2D display (Qgis canvas) but also the ZNS and the 1D panel. The video is created by concatenating the PNG images. The record starts at the current date and stops by clicking on the button.
Template
The Template action saves the current study in the sites.sqlite database, adding an entry to the following table :
simulations,
dates_resultats,
parametres,
metis_options,
injections.
If a study with the same name already exists, the old one is overwritten.
opens a panel to display the lithology of the boreholes contained in the database,
displays the results of the inversion process, like the iteration dependent mean error value and the correlation between measured and calculated piezometric values,
opens the inversion process dialog, where user can choose between pilot points or permeability field method,
launches the simulation (cf. THYRSIS panel to define injection parameters),
the first combobox contains the simulation dates. Selecting a date displays the result layer at the selected date,
the cursor offers a faster way to access to a specific date,
starts an animation, displaying result at each date,
displays the first simulation date results,
displays the previous simulation date results,
displays the next simulation date results,
displays the last simulation date results,
the checkbox “Adjust scale” adjusts the color scale for each date displayed. If uncheck, the color scale is calculated on the whole simulation,
the two next comboboxes are used to select the variable to display and its unit,
allows to request the value of the result layer on the selected location in the QGIS canvas,
opens a matplotlib panel displaying evolution graphs at interest points,
opens a matplotlib panel displaying the mass balance, which shows the time dependent values of the injected masses, distributed between the unsaturated and saturated zones, and leaving the domain.
Figure 53 Example of 2D panel, with result layer displayed after a simulation
The 2D panel is the main QGIS interface, more precisely the map canvas. It is connected with the “Layers” panel and displays on various background maps the specific parameters of the hydrogeological model :
mesh nodes,
mesh,
points of interest,
boreholes,
contours,
potential (piezometric level),
permeability,
topographic elevation,
aquifer substratum altitude,
Darcy’s velocity,
unsatured zone thickness,
groundwater thickness,
reference potential, which is obtained by interpolation of the piezometric data, and could have been used to build the hydrogeological model, through inversion computation.
Injection zones, defined in the THYRSIS panel, are displayed in the 2D panel.
At the end of a computation, a result layer is also displayed in the 2D panel. The dates and the variable to display are controled by the Menu bar.
Figure 54 Example of a Thyrsis panel, used to define injection zones
The THYRSIS panel is used to define injection zones.
The first values to define are :
the simulation name,
the chemical to use, selected in the drop-down menu,
the partition coefficient and the solubility limit in case of mass injection,
the starting date of the simulation, with format DD/MM/YYYY [hh:mm:ss],
a set of tuples (duration, step), that define the simulation output dates.
Then the type of injection can be chosen by clicking on :
flux,
concentration,
mass.
Figure 55 Example of THYRSIS panel to define a flux injection zone
According to the injection type selected, relevant fields may differ, as mention in the injection part. A non-null water volume is a leak that induces a transient flow simulation.
It is possible to add as many injections as wanted by clicking many times on the selected button, but it is impossible to mix several injection types in the same simulation.
The center coordinates of the injection zone can be defined manually or by clicking on \({\odot}\) and then clicking on the wanted location in the map canvas.
Similary the surface can be defines manually in the Area field, or automatically calculated by clicking on and then drawing the shape on the map canvas.
Below the THYRSIS panel, three buttons allows to cancel or save the modifications, and also to launch computation :
“Cancel” will cancel current modifications and reload the saved informations that are stored in the database.
“Save” will save the modifications (data and geometry) in the computation database.
will launch the simulation.
Note
If you made modifications, it is important to save it by clicking on “Save” before launching the computation by clicking on .
Figure 56 Example of stratigraphic log diplayed in the borehole panel
When stratigraphic and lithologic informations are provided in the sites.sqlite database, the stratigraphic log can be displayed in the borehole panel. Click on in the bar menu and select the borehole directly in the map canvas.
If the fracturated rate is known, it is displayed in the stratigraphic log.
This graphic representation is based on differents codifications presented in the Stratigraphic logs annexe :
codification on the geological time scale (BRGM),
RGB colorimetric codification of the Commission for the geological map of the world (CGMW),