diff --git a/docs/_quarto-tutorial.yml b/docs/_quarto-tutorial.yml
index 52af07c..7791000 100644
--- a/docs/_quarto-tutorial.yml
+++ b/docs/_quarto-tutorial.yml
@@ -6,7 +6,7 @@ book:
chapters:
- index.qmd
- tutorial.qmd
- - tutorial_QGIS_Tim.qmd
+ - tutorial_Rijsenhout.qmd
- tutorial_TheHague.qmd
format:
diff --git a/docs/_quarto-website.yml b/docs/_quarto-website.yml
index 07c4d7e..7ee48f7 100644
--- a/docs/_quarto-website.yml
+++ b/docs/_quarto-website.yml
@@ -31,7 +31,7 @@ website:
align: left
contents:
- href: tutorial.qmd
- - href: tutorial_QGIS_Tim.qmd
+ - href: tutorial_Rijsenhout.qmd
- href: tutorial_TheHague.qmd
format:
diff --git a/docs/tutorial.qmd b/docs/tutorial.qmd
index 2c51638..0af6447 100644
--- a/docs/tutorial.qmd
+++ b/docs/tutorial.qmd
@@ -3,13 +3,13 @@ title: "{{< fa solid book >}} Tutorials"
listing:
type: grid
contents:
- - tutorial_QGIS_Tim.qmd
+ - tutorial_Rijsenhout.qmd
- tutorial_TheHague.qmd
image: figures/logo/iMOD-tutorial.svg
description: 'Learn how to use QGIS Tim plugin'
index: "5"
---
-Modelling with analytic elements is not very common. The QGIS-Tim plugin makes it more easy to build your model in a graphical user interface starting from simple and making it as complex as you like.
+Modelling with analytic elements is not very common. The QGIS-Tim plugin makes it more easy to build your model in a graphical user interface starting from simple and making it as complex as you like.
Therefore, we have compiled a set of tutorials to help you get started and give you an overview of the capabilities of analytical elements.
diff --git a/docs/tutorial_QGIS_Tim.qmd b/docs/tutorial_Rijsenhout.qmd
similarity index 88%
rename from docs/tutorial_QGIS_Tim.qmd
rename to docs/tutorial_Rijsenhout.qmd
index d6fa6b9..5a04f3e 100644
--- a/docs/tutorial_QGIS_Tim.qmd
+++ b/docs/tutorial_Rijsenhout.qmd
@@ -1,13 +1,13 @@
---
-title: "Tutorial Building Pit"
+title: "Building Pit Rijsenhout"
---
## Requirements
-- Be sure QGIS version 3.22.00 or higher is installed.
-- Be sure the ``gistim`` Python package is installed (see [installation](install.qmd) for instructions).
+- Be sure QGIS version 3.22.00 or higher is installed.
+- Be sure the ``gistim`` Python package is installed (see [installation](install.qmd) for instructions).
- Download the tutorial material. [Follow this link.](https://deltares.thegood.cloud/s/zryWnfjT43Ryoea)
-- Installation of the QGIS-Tim plugin and the MOD plugin is part of this Tutorial. The necessary ZIP files are included in the tutorial material.
-- Internet connections is optional during this Tutorial. It is only required for installation of additional plugins and the use of an online topographic background map.
+- Installation of the QGIS-Tim plugin and the MOD plugin is part of this Tutorial. The necessary ZIP files are included in the tutorial material.
+- Internet connections is optional during this Tutorial. It is only required for installation of additional plugins and the use of an online topographic background map.
## Description
In this tutorial, you will learn how to:
@@ -23,41 +23,41 @@ Calculation of a pumping well extraction.
## Introduction case *Rijsenhout*
-In one of the fields west of Rijsenhout industrial activities are planned, the construction of green houses. From historical sources it is known that some bombs from World War 2 are still in the surface. Before building starts, these bombs are removed for safety reasons. The question is: what dewatering is necessary to be able to remove the bombs in dry conditions? The contractor plans to drill sheet piles to prevent the excavation from collapsing. The groundwater level inside the building pit is lowered with a set of wells just in the top of the aquifer within the building pit. All pumped water is infiltrated over a stretch of wells near the existing buildings.
-At the project location, the top layer is 12 m thick and consists of a low permeable combination of clay and peat. Below this top layer we find an aquifer of 50 m. The Westeinderplassen, a lake system with a depth of 3 m., is located East of the project area.
+In one of the fields west of Rijsenhout industrial activities are planned, the construction of green houses. From historical sources it is known that some bombs from World War 2 are still in the surface. Before building starts, these bombs are removed for safety reasons. The question is: what dewatering is necessary to be able to remove the bombs in dry conditions? The contractor plans to drill sheet piles to prevent the excavation from collapsing. The groundwater level inside the building pit is lowered with a set of wells just in the top of the aquifer within the building pit. All pumped water is infiltrated over a stretch of wells near the existing buildings.
+At the project location, the top layer is 12 m thick and consists of a low permeable combination of clay and peat. Below this top layer we find an aquifer of 50 m. The Westeinderplassen, a lake system with a depth of 3 m., is located East of the project area.
![Geological cross section near Rijssenhout \[source: BRO GeoTOP v1.4.1\]](figures/tutorial/GeoTOP-Verticaledoorsnede.png){width=100% #fig-GeoTOP-Verticaledoorsnede}
-## Getting Started
+## Getting Started
Let's first configure the gistim python installation again, to be sure that QGIS can find the gistim software.
-(@) Open the Deltaforge prompt (search in Windows Start for "Deltaforge Prompt"). A black window pops up.
+(@) Open the Deltaforge prompt (search in Windows Start for "Deltaforge Prompt"). A black window pops up.
(@) In this window type ``python -m gistim configure`` and press ENTER.
-(@) You can close this window now.
+(@) You can close this window now.
(@) Launch QGIS from your START menu, from your desktop or click on …\\QGIS3.28.0\\bin\\qgis-bin.exe.
> **Intermezzo:** *QGIS language settings*
->
-> Perhaps your QGIS was installed in another language than English. Because the Tutorial refers to the English version, let’s change to English.
->
-> a. From the main menu click on *Settings* and select *Options* (e.g. in Dutch *Extra* and *Opties*).
-> a. In the new window go to the *General section* (Dutch: *Algemeen*) on the left.
-> a. Check the box to allow *Override System Locale* (Dutch: *Landinstellingen negeren*) and expand this sub menu.
-> a. From the drop-down menu *“User interface translation”* (Dutch: *Vertaling gebruikers-interface*) select *American English* and click *OK*.
+>
+> Perhaps your QGIS was installed in another language than English. Because the Tutorial refers to the English version, let’s change to English.
+>
+> a. From the main menu click on *Settings* and select *Options* (e.g. in Dutch *Extra* and *Opties*).
+> a. In the new window go to the *General section* (Dutch: *Algemeen*) on the left.
+> a. Check the box to allow *Override System Locale* (Dutch: *Landinstellingen negeren*) and expand this sub menu.
+> a. From the drop-down menu *“User interface translation”* (Dutch: *Vertaling gebruikers-interface*) select *American English* and click *OK*.
> a. Close QGIS and open it again to activate your language change.
We start with the creation of a new QGIS project.
-(@) From the main menu click on *Project* and select *New*.
+(@) From the main menu click on *Project* and select *New*.
The case in this tutorial is located in The Netherlands, so next we select the appropriate projection.
-(@) From the main menu click on *Project* and select *Properties*.
+(@) From the main menu click on *Project* and select *Properties*.
(@) In the *Properties* window select the category *CRS*, search for “EPSG:28992” and you find “Amersfoort / RD New”. Select this option and click the *Apply* button, followed by the *OK* button to close the window.
> In case your work is mostly in The Netherlands and in the “Amersfoort / RD New” projection, consider making this your default projection.
->
+>
> - From the main menu click on *Settings* and select *Options...*.
> - In the section *CRS and Transforms* select *CRS (handling)*, pick the radio button *Use a default CRS* and select “EPSG:28992 -Amersfoort / RD New”.
> - Click OK.
@@ -68,39 +68,39 @@ This is the moment to download/import four plugins needed for this tutorial:
- the QGIS-Tim plugin. The development version imported from ZIP file.
- the iMOD plugin. The development version imported from ZIP file.
-- the Value Tool. The official version installed via QGIS (internet connection required).
+- the Value Tool. The official version installed via QGIS (internet connection required).
- the PDOK plugin. The official version installed via QGIS (internet connection required).
This plugin gives access to a large database from which we will load the topographic maps and use the navigation option.
-(@) Go to *Plugins* from the main menu and select *Manage and Install Plugins…* to open the plugin window.
+(@) Go to *Plugins* from the main menu and select *Manage and Install Plugins…* to open the plugin window.
(@) On the left section select *Install from ZIP*.
-(@) Click the Browse button ({width=4%}) and from the tutorial dataset select the ZIP file "QGIS-Tim_Tutorial\\QGIS-iMOD-plugin.zip".
+(@) Click the Browse button ({width=4%}) and from the tutorial dataset select the ZIP file "QGIS-Tim_Tutorial\\QGIS-iMOD-plugin.zip".
(@) Click *Install Plugin*.
-(@) In the same way, install the QGIS-Tim plugin using the ZIP file "QGIS-Tim_Tutorial\\QGIS-Tim-plugin.zip".
+(@) In the same way, install the QGIS-Tim plugin using the ZIP file "QGIS-Tim_Tutorial\\QGIS-Tim-plugin.zip".
If you have an internet connection continue with the installation of the next two plugins from the QGIS plugin library.
-(@) From the left section, select the group *All* to see all available plugins.
-(@) Search for “Value Tool” and install it.
-(@) Search for “PDOK services plugin” and install it.
+(@) From the left section, select the group *All* to see all available plugins.
+(@) Search for “Value Tool” and install it.
+(@) Search for “PDOK services plugin” and install it.
(@) Make sure that under *Plugins > Manage and Install Plugins > Installed* now the 4 added plugins are checked.
-(@) *Close* the Plugins window.
+(@) *Close* the Plugins window.
-See in the toolbar section of QGIS that the plugins are installed:
+See in the toolbar section of QGIS that the plugins are installed:
-- iMOD Toolbar {width=30%}
+- iMOD Toolbar {width=30%}
- QGIS-Tim {width=4%}
- Value Tool {width=4%}
- PDOK Services Plugin {width=4%}
-Further in this Tutorial we will use some default toolbars that might be hidden at the moment. Let's check that and unhide if necessary.
+Further in this Tutorial we will use some default toolbars that might be hidden at the moment. Let's check that and unhide if necessary.
(@) Select *View* from the main menu and choose *Toolbars*.
(@) Be sure the "Advanced Digitizing Toolbar", the "Snapping Toolbar" and the "Attributes Toolbar" are checked.
### Prepare your project
-For navigation purposes, let's load a topographic map for The Netherlands from the online PDOK database.
+For navigation purposes, let's load a topographic map for The Netherlands from the online PDOK database.
-> No internet connection? Follow the next steps to import a simple PNG file as a background.
+> No internet connection? Follow the next steps to import a simple PNG file as a background.
>
> - Go to *Layer* in the main menu, go to *Add layer* and select *Add Raster layer*.
> - Use the browse button ({width=4%}) and from the tutorial material select “…\\QGIS-Tim_Tutorial\\dbase\\TopographicMapRijssenhout.png”.
@@ -109,10 +109,10 @@ For navigation purposes, let's load a topographic map for The Netherlands from t
> - Continue after step @continuewithoutinternet.
-(@) If you do have an internet connection click on the PDOK plugin button ({width=4%}) to open the "PDOK Services Plugin" window.
-(@) From the tab *PDOK Services* search for "pastel" and you will find a WMTS type layer called "BTRM Achtergrondkaart WMTS".
+(@) If you do have an internet connection click on the PDOK plugin button ({width=4%}) to open the "PDOK Services Plugin" window.
+(@) From the tab *PDOK Services* search for "pastel" and you will find a WMTS type layer called "BTRM Achtergrondkaart WMTS".
(@) Select the layer.
-(@) In the section "laag toevoegen" click the button *Onder*.
+(@) In the section "laag toevoegen" click the button *Onder*.
(@) Close the PDOK window.
Our project area is near the town of Rijsenhout so let's navigate to that town using the PDOK plugin.
@@ -137,47 +137,47 @@ Let's save this project to be able to return to it later or in case of a crash o
(@) Go to *Project* in the main menu, select *Save As* and select a folder and a file name for your project, e.g. “…\\QGIS-Tim_Tutorial\\Rijsenhout.qgz”
-## Start your Tim model
-Now we are ready to activate the QGIS-Tim plugin.
+## Start your Tim model
+Now we are ready to activate the QGIS-Tim plugin.
(@) Click on the QGIS-Tim plugin button ({width=4%}) and the QGIS-Tim panel appears.
-(@) Go to the tab *GeoPackage*.
Here we will create an empty database (geopackage) to store all elements and parameters for the model.
+(@) Go to the tab *GeoPackage*.
Here we will create an empty database (geopackage) to store all elements and parameters for the model.
(@) Click the *New* button to create the GeoPackage and save it for instance in the folder with your tutorial data, e.g. "..\\QGIS-Tim_Tutorial\\dbase\\case-Rijsenhout.gpkg".
Your window looks like in @fig-Panel-QGIS-Tim.
{width=60% #fig-Panel-QGIS-Tim}
-(@) Check in the *Layers* panel on the left that your new geopackage is added as a group.
A sub group **timml** for the steady state model input, the sub group **ttim** for the transient model input and a series of output formats (vector/mesh/raster).
+(@) Check in the *Layers* panel on the left that your new geopackage is added as a group.
A sub group **timml** for the steady state model input, the sub group **ttim** for the transient model input and a series of output formats (vector/mesh/raster).
-If you had no introduction to the Tim plugin, read the Intermezzo below for a general explanation of the components.
+If you had no introduction to the Tim plugin, read the Intermezzo below for a general explanation of the components.
-> **Intermezzo:** *introduction Tabs on the Tim panel*
->
-> a. GeoPackage: an overview of the elements in your geopackage. In case you switch to transient modelling, an extra column with *ttim* elements is added.
-> a. Elements: a list of 14 Tim elements from which you can build your model.
-> a. Compute: here you can define your domain and cell size, decide if your model is transient or not and change the output name.
-> a. Extract: open an existing 3d geohydrological model (NC file) and extract the data for your project area.
+> **Intermezzo:** *introduction Tabs on the Tim panel*
+>
+> a. GeoPackage: an overview of the elements in your geopackage. In case you switch to transient modelling, an extra column with *ttim* elements is added.
+> a. Elements: a list of 14 Tim elements from which you can build your model.
+> a. Compute: here you can define your domain and cell size, decide if your model is transient or not and change the output name.
+> a. Extract: open an existing 3d geohydrological model (NC file) and extract the data for your project area.
Now we are ready to define our first steady state model by parameterizing our Aquifer.
## Model 1: single aquifer
We start with a very simple 'model', only the parameters of a single aquifer.
-Important message: all editing of model parameters and model elements you do in the *Layers* panel on the left!
+Important message: all editing of model parameters and model elements you do in the *Layers* panel on the left!
(@) So select the layer "timml Aquifer:Aquifer" on the left.
(@) Click your right mouse button and from the menu select *Attribute Table* to open the table in a new window.
**NB** Alternative is to press F6 or use the button *Open Attribute Table* ({width=3%}).
(@) Start the editing mode with a click on the *Toggle Editing Mode* button ({width=4%}).
-(@) Hover with your mouse over the buttons and find the *Add Feature* button ({width=20%}). Perhaps number 5 from left.
-(@) Add a new feature with this *Add Feature* button. In this case a feature is a hydrological layer.
-(@) Fill the feature with the values from the table below:
+(@) Hover with your mouse over the buttons and find the *Add Feature* button ({width=20%}). Perhaps number 5 from left.
+(@) Add a new feature with this *Add Feature* button. In this case a feature is a hydrological layer.
+(@) Fill the feature with the values from the table below:
{width=80%}
| parameter | value | unit | comment|
|--- |--- |--- |--- |
-| fid | Autogenerate | [-] | ID is autogenerated by QGIS |
-| Layer | 0 | [-] |(see remark below) |
+| fid | Autogenerate | [-] | ID is autogenerated by QGIS |
+| Layer | 0 | [-] |(see remark below) |
| aquifer_top | +2 | [m MSL] | |
| aquifer_bottom | -10 | [m MSL]| |
| aquifer_c | *NULL* | [d]| |
@@ -186,7 +186,7 @@ Important message: all editing of model parameters and model elements you do in
| semiconfined_head | *NULL* | [m MSL] | |
**NB!** In the real world counting starts with 1. However, Tim is programmed in Python and in Python counting starts with 0. You will get used to it.
-Also be aware that the *Aquifer* element you just edited is just a table. This is indicated with the icon ({width=3%}) just before the layer name in the left panel. For Tim this means that the properties in this table apply for the full model domain. We will introduce some inhomogeneities later as polygons within that domain.
+Also be aware that the *Aquifer* element you just edited is just a table. This is indicated with the icon ({width=3%}) just before the layer name in the left panel. For Tim this means that the properties in this table apply for the full model domain. We will introduce some inhomogeneities later as polygons within that domain.
(@) Save all changes with the *Save Edits* button ({width=4%}).
(@) Stop the *Editing Mode* with a click on the button ({width=4%}).
@@ -195,19 +195,19 @@ Also be aware that the *Aquifer* element you just edited is just a table. This i
(@) Click the button *Set to current extent*.
(@) Start the calculation with a click on *Compute*.
-A black Python.exe window pops up indicating that the TIM calculation started on the background. You can ignore this window but keep it open. Of course you van minimize it. If the calculation was completed successful, you will see this echo in QGIS.
-{width=80%}.
+A black Python.exe window pops up indicating that the TIM calculation started on the background. You can ignore this window but keep it open. Of course you van minimize it. If the calculation was completed successful, you will see this echo in QGIS.
+{width=80%}.
-After the calculation you see that the result is automatically added to a new geopackage, probably called "case-Rijsenhout output". Results are saved both as mesh and raster. For each format a separate sub group is created. Although these layers / groups are checked, the data is not visible. That is because the geopackage was added last, and QGIS adds layers at the end of the list. Let's move the layer "pastel" to the background.
+After the calculation you see that the result is automatically added to a new geopackage, probably called "case-Rijsenhout output". Results are saved both as mesh and raster. For each format a separate sub group is created. Although these layers / groups are checked, the data is not visible. That is because the geopackage was added last, and QGIS adds layers at the end of the list. Let's move the layer "pastel" to the background.
(@) Select the layer "pastel" and drag it with your left mouse button to the bottom of the list of layers.
The calculation result is now visible and we see a raster with just the value 0, not a very exciting result because no other elements are present yet. Let's now add a well element.
> **Intermezzo:** *Mesh and Raster format explained*
->
+>
> - Mesh: an unstructured grid usually with temporal and other components. Preferred format in QGIS-Tim to animate temporal data, to create cross sections or evaluate values at your mouse position.
-> - Raster: is made up of pixels (also referred to as grid cells). They are regularly spaced and square. Preferred format in QGIS-Tim to perform calculations with the *Raster Calculator* tool.
+> - Raster: is made up of pixels (also referred to as grid cells). They are regularly spaced and square. Preferred format in QGIS-Tim to perform calculations with the *Raster Calculator* tool.
## Model 2: single well
@@ -221,15 +221,15 @@ See that layer "timml Well:DewateringWell" is added to the timml sub group in yo
(@) Select the layer "timml Well:DewateringWell".
(@) Click your right mouse button and select the *Toggle Editing Mode* ({width=4%}).
(@) To add a new well (feature) to the layer click the *Add Point Feature* button ({width=4%}).
-(@) With your mouse, click on a location in the centre of the group of bombs.
-(@) Fill the feature with the values from the example below.
+(@) With your mouse, click on a location in the centre of the group of bombs.
+(@) Fill the feature with the values from the example below.
{width=40%}
(@) Click *OK*.
(@) Click the *Toggle Editing* button ({width=4%}) and if you are asked to save changes and select *Save*.
An alternative way to save your edits is to click the *Save Layer Edits* button ({width=4%}) and then stop editing with the *Toggle Editing* button ({width=4%}).
-Next step is to rerun the model including the new well.
+Next step is to rerun the model including the new well.
(@) In QGIS-Tim go to the tab *GeoPackage* and see that the Well is added.
(@) Go to the tab *Compute* and click on *Compute*.
@@ -241,8 +241,8 @@ The output in your ouput group "case-Rijsenhout output" is directly overwritten
Do you like to see the values of the calculated Head under you mouse?
(@) Deselect the output sub group *mesh*.
-(@) Select the *Value Tool* button ({width=4%}).
-(@) Hover over the area and in the “Value Tool” panel you see the value within the raster file at your mouse location.
+(@) Select the *Value Tool* button ({width=4%}).
+(@) Hover over the area and in the “Value Tool” panel you see the value within the raster file at your mouse location.
{width=40% #fig-figure_ValueTool-in-QGIS fig-align="left"}
@@ -254,11 +254,11 @@ In the empty cross-section we can add a selection of (geological) layers. For no
(@) On the iMOD Cross Section Plot click on the button *Select location* and draw your cross-section line from north to south (right mouse button to close the line).
(@) From the dropdown menu on the left of this toolbar, select the raster ({width=3%}) named "case-Rijsenhout-head_layer_0".
(@) Click the button *Add* to add this layer to the cross-section manager below.
-(@) Click the button *Plot* to draw this layer in the cross section.
+(@) Click the button *Plot* to draw this layer in the cross section.
-Your screen might look like @fig-figure_HeadsInCrosssection01.
+Your screen might look like @fig-figure_HeadsInCrosssection01.
-TIP: If you do not see any line, perhaps the axes are not defined well. To view all data, click you right mouse button in the figure and select the option “View All”. The alternative is to click on the small A symbol ({width=3%}) in the lower left of the chart.
+TIP: If you do not see any line, perhaps the axes are not defined well. To view all data, click you right mouse button in the figure and select the option “View All”. The alternative is to click on the small A symbol ({width=3%}) in the lower left of the chart.
{width=100% #fig-figure_HeadsInCrosssection01}
@@ -275,7 +275,7 @@ A single aquifer is of course not enough to describe the geology of our project
The Holocene top layer is not a distinguished layer but a resistance on top of layer 0.
-(@) Fill the features with the values from the table below:
+(@) Fill the features with the values from the table below:
|fid|layer|aquifer_top|aquifer_bottom|aquitard_c|aquifer_k|semiconf_top|semiconf_head|
|---|---|---|---|---|---|---|---|
@@ -292,34 +292,34 @@ Let's now check the calculated heads for the 4 layers, first with the *Value Too
(@) From the layers select the layer "case-Rijsenhout-output:raster" and deselect the "case-Rijsenhout-output:mesh".
(@) Activate the *Value Tool* ({width=4%}).
(@) In the *Value Tool* panel go to the tab *Options*.
-(@) For *Show Layers* choose "Visible layers" and for *Show bands* choose "Active bands".
-(@) Return to the tab *Table* and hoover over the calculated values, especially near your well.
-(@) Activate the *iMOD Cross section* widget ({width=3%}) to create your cross section with the four layers (Need instructions? Go back to step @imodcrosssection for support).
**NB** You don't need to 'add' all four layers separately. Just use the raster file "case-Rijsenhout-head_layer_0" and from the dropdown menu *Variable:* select the four bands/layers.
+(@) For *Show Layers* choose "Visible layers" and for *Show bands* choose "Active bands".
+(@) Return to the tab *Table* and hoover over the calculated values, especially near your well.
+(@) Activate the *iMOD Cross section* widget ({width=3%}) to create your cross section with the four layers (Need instructions? Go back to step @imodcrosssection for support).
**NB** You don't need to 'add' all four layers separately. Just use the raster file "case-Rijsenhout-head_layer_0" and from the dropdown menu *Variable:* select the four bands/layers.
## Model 4: add a set of sheet piles
-Now we introduce a new but frequently used element: the sheet pile.
+Now we introduce a new but frequently used element: the sheet pile.
(@) Return to the QGIS-Tim panel and go to the tab *Elements*.
(@) Click on the element "Leaky Line Doublet".
-(@) Give the layer a name, e.g. "sheetpile".
+(@) Give the layer a name, e.g. "sheetpile".
(@) In your geopackage select the layer "timml Leaky Line Doublet:sheetpile" and start editing using the *Toggle Editing Mode* button ({width=4%}).
-(@) In the editing mode, the *Add Line Feature* button ({width=4%}) is available. Click the button and draw a box around the 5 bombs (click left to start, click right to close).
+(@) In the editing mode, the *Add Line Feature* button ({width=4%}) is available. Click the button and draw a box around the 5 bombs (click left to start, click right to close).
(@) In the *Feature Attributes* window make resistance = 1500 and layer = 0.
(@) Click *OK*.
-You probably tried to close the building pit but either your lines cross or do not close actually when zooming in. To close the pit, let's snap both start and end point together.
+You probably tried to close the building pit but either your lines cross or do not close actually when zooming in. To close the pit, let's snap both start and end point together.
(@) First enable the snapping mode with a click on the button ({width=4%}).
(@) Then activate the *Vertex Tool* with a click on the button {width=4%}.
(@) Move your mouse to the start point location and when it shows a small circle, click your left mouse button.
-(@) Next, move your mouse towards the end point and you see the snapping active: the end point gets a small fuchsia box.
-(@) Click your left mouse and you see that both vertices now connect.
+(@) Next, move your mouse towards the end point and you see the snapping active: the end point gets a small fuchsia box.
+(@) Click your left mouse and you see that both vertices now connect.
-Perhaps you are not satisfied with your created building pit.
+Perhaps you are not satisfied with your created building pit.
-(@) You can play around with other nice options in QGIS to move ({width=4%}), rotate ({width=4%}) or scale ({width=4%}) your pit.
+(@) You can play around with other nice options in QGIS to move ({width=4%}), rotate ({width=4%}) or scale ({width=4%}) your pit.
(@) End the *Toggle Editing Mode* with {width=4%} and save your changes.
-(@) In the QGIS-Tim panel go to *Compute*.
+(@) In the QGIS-Tim panel go to *Compute*.
(@) Rerun your model and check the calculated heads.
Disappointed in the effect of the sheet piles on the calculated head? You are right if you expected no major effect because of the shallow penetration in this high permeable thick aquifer. If we would increase the aquitard_c between layers 0 and 1, perhaps we would see more effect of the sheet piles.
@@ -334,7 +334,7 @@ Disappointed in the effect of the sheet piles on the calculated head? You are ri
## {width=4%} How to reopen Tim after closing QGIS.
Before we introduce you another new Tim element we must explain you how to reopen your Tim project after a closure of QGIS. A crash might close QGIS or just a long day of work is a reason to shut down your laptop. Always remember two things:
-- your geopackage (*.gpkg) contains the complete set of model features and parameter values. Every time you change a model element, you will remember, you had to save the changes. So the GPKG file is always up to date.
+- your geopackage (*.gpkg) contains the complete set of model features and parameter values. Every time you change a model element, you will remember, you had to save the changes. So the GPKG file is always up to date.
- your QGIS project is saved in the *.qgz file and contains all the added layers and their format (legend, line color, labels etc.).
Now let's experience what it is to close and open QGIS including your QGIS-Tim project.
@@ -343,14 +343,14 @@ Now let's experience what it is to close and open QGIS including your QGIS-Tim p
(@) Select *Project* from the main menu and choose *New*.
(@) Again select *Project* but now choose *Open Recent* and find your own project in the list.
(@) Open the QGIS-Tim plugin with a click on {width=4%}.
-
+
*The Tim panel is empty but why?*
-It is not possible for the plugin to read your geopackage as it is loaded in the layer list. To load your Tim project you *must* open the geopackage from the Tim panel. This will add the geopackage to the Layers panel next to the existing one with the same name, which makes it at least confusing. We better remove the original Geopackage first before opening it with Tim. Disadvantage is that we also remove the calculated results including all your formatting efforts. The developers try to fix that in the next version.
-The following steps are the right way to go for now:
-
-(@) Select layer *case-Rijsenhout input* in the panel on the left in order to remove it.
-(@) Click your right mouse button and select *Remove Group...* and click *OK*.
-(@) In the same way remove layer *case-Rijsenhout output*.
+It is not possible for the plugin to read your geopackage as it is loaded in the layer list. To load your Tim project you *must* open the geopackage from the Tim panel. This will add the geopackage to the Layers panel next to the existing one with the same name, which makes it at least confusing. We better remove the original Geopackage first before opening it with Tim. Disadvantage is that we also remove the calculated results including all your formatting efforts. The developers try to fix that in the next version.
+The following steps are the right way to go for now:
+
+(@) Select layer *case-Rijsenhout input* in the panel on the left in order to remove it.
+(@) Click your right mouse button and select *Remove Group...* and click *OK*.
+(@) In the same way remove layer *case-Rijsenhout output*.
(@) In the QGIS-Tim panel on the tab *GeoPackage* click *Open* and open the GPKG file containing your model.
(@) Your geopackage is now added to the Tim panel and the Layers list. Don't forget to move down the background layer "pastel".
@@ -362,7 +362,7 @@ To prevent damage on the private properties caused by the dewatering, the author
(@) On the tab *Elements* in QGIS-Tim click on the element "Line Sink Ditch".
(@) Give the layer a name, e.g. "InfiltrationWell".
-See that layer "timml Line Sink Ditch:InfiltrationWell" is added to the timml sub group in your geopackage with 'line' as geometry. For the transient sum a table with the same name is added to the ttim sub group. Next step is to add the location of the well and its capacity.
+See that layer "timml Line Sink Ditch:InfiltrationWell" is added to the timml sub group in your geopackage with 'line' as geometry. For the transient sum a table with the same name is added to the ttim sub group. Next step is to add the location of the well and its capacity.
(@) Click the *Measure Line" button ({width=4%}) from the "Attributes Toolbar".
(@) With your left mouse button try to get an idea what a distance of 400 m. looks like. Close the window.
@@ -370,9 +370,9 @@ See that layer "timml Line Sink Ditch:InfiltrationWell" is added to the timml su
(@) Start editing with the *Toggle Editing Mode*.
(@) Use the *Add Line Feature* button ({width=4%}) to draw a 400 m line (click left to start, click right to close).
(@) Fill these feature within the table: discharge = -1000 (negative abstraction is into the model), resistance = 1, width = 1, layer = 0! not *NULL* ;-). By default order = 4.
-(@) Click *OK*.
-(@) Stop editing the layer, save your changes.
-(@) Rerun the model.
+(@) Click *OK*.
+(@) Stop editing the layer, save your changes.
+(@) Rerun the model.
(@) Analyse the results with the iMOD Cross Section widget. Your graph may look like @fig-figure_HeadsInCrosssection02.
{width=70% #fig-figure_HeadsInCrosssection02}
@@ -383,99 +383,99 @@ We have the Value Tool and the Cross section tool to check for calculated values
(@) On the tab *Elements* in QGIS-Tim click on the element "Observation".
(@) Give the layer a name, e.g. "Piezometers".
(@) Start editing with the *Toggle Editing Mode* and use *Add Point Feature* button ({width=4%}) to add 5 points in line with 1 really close to the dewatering well (see @fig-figure_ObservationsLocations).
-(@) Save your changes and rerun the model.
+(@) Save your changes and rerun the model.
{width=80% #fig-figure_ObservationsLocations}
-The observed values are saved in a new layer under the existing group "case-Rijsenhout output" and sub group "vector".
+The observed values are saved in a new layer under the existing group "case-Rijsenhout output" and sub group "vector".
(@) In this sub group select the layer "case-Rijsenhout-timml Observation:Piezometers" and check the calculated values within this layer by opening the Attribute Table (F6).
(@) Close the Attribute Table.
-One way of displaying the calculated results is still missing: contour lines.
+One way of displaying the calculated results is still missing: contour lines.
-(@) In QGIS-Tim go to the tab *Compute*.
+(@) In QGIS-Tim go to the tab *Compute*.
(@) In the section *Contour* select layer "case-Rijsenhout-head_layer_0".
(@) Define the contours from -9 to -5 with an increment of 0.10 and click *Export contours*.
-A layer with your contours is saved in the group "case-Rijsenhout-output:vector". Your map hopefully looks like @fig-figure_Contourlines. If not, try to change the contour settings and create a better contour map.
+A layer with your contours is saved in the group "case-Rijsenhout-output:vector". Your map hopefully looks like @fig-figure_Contourlines. If not, try to change the contour settings and create a better contour map.
{width=90% #fig-figure_Contourlines}
## Model 7: determine the influence of the nearby lake
East of the project area there is a large lake. What will be the influence of this lake? Let's find out.
-In Tim a lake is added as an inhomogeneity within the total model domain. This domain is described in table "timml Aquifer:Aquifer" we filled in with the first model.
+In Tim a lake is added as an inhomogeneity within the total model domain. This domain is described in table "timml Aquifer:Aquifer" we filled in with the first model.
(@) On the tab *Elements* in QGIS-Tim click on the element "Polygon Semi-Confined Top".
(@) Give the layer a name, e.g. "Lake" and find out that layer "timml Polygon Semi-Confined Top:Lake" is added to the geopackage.
-(@) Start editing with the *Toggle Editing Mode* and use *Add Polygon Feature* button ({width=4%}) to add the Lake contour.
**NB** Try to minimize the number of segments because each segment slows down the calculations. With QGIS it is tempting to use detail but the influence of detail on the outcome is small.
+(@) Start editing with the *Toggle Editing Mode* and use *Add Polygon Feature* button ({width=4%}) to add the Lake contour.
**NB** Try to minimize the number of segments because each segment slows down the calculations. With QGIS it is tempting to use detail but the influence of detail on the outcome is small.
(@) Fill the Lake feature with the following parameters: aquitard_c = 500, semiconf_top = -7 (bottom of the Lake) and semiconf_head = -2 (water level).
-(@) Save your changes and rerun the model.
-(@) Analyse the results and perhaps create contour lines from the tab "Compute".
+(@) Save your changes and rerun the model.
+(@) Analyse the results and perhaps create contour lines from the tab "Compute".
{width=100% #fig-figure_Contourlines02}
-You can not miss the effect of the Lake in this calculation. In many cases you like to analyse or even save this difference. Let's see how we can use Tim and QGIS to isolate the effect. Therefor we rerun the model without the Lake and rename the output.
+You can not miss the effect of the Lake in this calculation. In many cases you like to analyse or even save this difference. Let's see how we can use Tim and QGIS to isolate the effect. Therefor we rerun the model without the Lake and rename the output.
(@) In QGIS-Tim go to the tab *Geopackage*.
(@) In the list of Steady State elements switch off "timml Polygon Semi-Confined Top:Lake".
-(@) Go to the tab *Compute*.
-(@) Click the button *Set path as ...* and change the name of the output, e.g. "case-Rijsenhout-NoLake".
-(@) Click *Compute* to run the model without the Lake.
+(@) Go to the tab *Compute*.
+(@) Click the button *Set path as ...* and change the name of the output, e.g. "case-Rijsenhout-NoLake".
+(@) Click *Compute* to run the model without the Lake.
-The new results are added as new layer to the project. Check for instance sub group "raster" under the group "case-Rijsenhout output". Now we use QGIS to calculate the head difference for layer 0. For raster calculations we prefer *Raster Calculator* over the *Mesh Calculator* in QGIS.
+The new results are added as new layer to the project. Check for instance sub group "raster" under the group "case-Rijsenhout output". Now we use QGIS to calculate the head difference for layer 0. For raster calculations we prefer *Raster Calculator* over the *Mesh Calculator* in QGIS.
-(@) From the QGIS main menu *Raster* select the tool *Raster Calculator* and the tool window opens.
+(@) From the QGIS main menu *Raster* select the tool *Raster Calculator* and the tool window opens.
(@) In the section *Raster Calculator Expression* you can create an expression using the elements from *Raster Bands* and *Operators*. Create this expression: "case-Rijsenhout-head_layer_001 - case-Rijsenhout-NoLake-head_layer_001".
(@) In *Output layer* give a new name for this new GeoTIFF file, e.g. "head-effect-lake".
(@) Click *OK* to start the calculation and see that the GeoTIFF file is added as a layer.
(@) In the properties of layer "head-effect-lake" go to the section *Symbology* and set the *Render type* to "Singelband pseudocolor".
-(@) See that the difference in head is maximum 3 meter.
+(@) See that the difference in head is maximum 3 meter.
(@) In *Geopackage* check again the Lake element.
(@) In *Compute* reset the path to the old name.
## Model 8: abstraction effect over time
-Now you are ready for the last step: make your model transient. We keep it simple: we run the model for 30 days and see the increasing effect of the well.
+Now you are ready for the last step: make your model transient. We keep it simple: we run the model for 30 days and see the increasing effect of the well.
(@) In *Compute* within the *output* section change "steady state" into "transient".
-(@) Go to *Geopackage* and see that only 5 elements have a transient component.
+(@) Go to *Geopackage* and see that only 5 elements have a transient component.
-By selecting the "transient" option, columns containing extra parameters necessary for a transient model have been unhidden now. First of all we see it in the layer "timml Aquifer:Aquifer" where columns "storage" and "porosity are visible now.
+By selecting the "transient" option, columns containing extra parameters necessary for a transient model have been unhidden now. First of all we see it in the layer "timml Aquifer:Aquifer" where columns "storage" and "porosity are visible now.
(@) In your geopackage select the layer "timml Aquifer:Aquifer".
-(@) Open its table, start the editing mode and for every layer set both acquifer_s and acquitard_s to 0.001. Acquitard_s (layer 0) = 0.25. You can use copy and paste to do it quickly.
+(@) Open its table, start the editing mode and for every layer set both acquifer_s and acquitard_s to 0.001. Acquitard_s (layer 0) = 0.25. You can use copy and paste to do it quickly.
**NB!** In transient mode, Tim can not handle aquitard with zero thickness so let's set the aquitard thickness to 0.1 m
(@) Change the values in the column "aquifer_top" to -17.0, -22.1, -27.1 and -47.1 m.
(@) In your geopackage select the layer "ttim Temporal Settings:Aquifer" to define the temporal properties of the model.
-(@) Open its table, start the editing mode and add a new feature.
+(@) Open its table, start the editing mode and add a new feature.
(@) Make sure that tmin=0.01, tmax=30, tstart=0, reference_date=2023-03-23 00:00:00.
-(@) In your geopackage select the layer "ttim Computation Times:Domain" in order to define the moments in time for which raster output is saved.
+(@) In your geopackage select the layer "ttim Computation Times:Domain" in order to define the moments in time for which raster output is saved.
(@) Open its table, start the editing mode and add 7 features / periods (click 7 times "Add feature").
-(@) In the column "time" add the moments 1,2,5,10, 15, 20 and 30.
-(@) Save and Close the table.
+(@) In the column "time" add the moments 1,2,5,10, 15, 20 and 30.
+(@) Save and Close the table.
-For the observations we can define a different set of output moments.
+For the observations we can define a different set of output moments.
-(@) In your geopackage select the layer "ttim Observation:Piezometers".
+(@) In your geopackage select the layer "ttim Observation:Piezometers".
(@) Open its table, start the editing mode and add 10 features / periods (click 10 times "Add feature").
-(@) In the column "time" add the moments 1,2,3,4,5,10,15,20,25 and 30.
-(@) Save and Close the table.
+(@) In the column "time" add the moments 1,2,3,4,5,10,15,20,25 and 30.
+(@) Save and Close the table.
Special attention for the Well package while it has 2 options to make it transient:
- Simple: a well can be switched on and off only once. Parameters are set in the stationary part of the model, so in "timml Well:DewateringWell".
- Detailed: a well can be switched on and off multiple times within the modelled period. Parameters are set in the transient part of the model, so in "ttim Well:DewateringWell".
-In this tutorial we show you the simple version.
+In this tutorial we show you the simple version.
-(@) In your geopackage select the layer "timml Well:DewateringWell".
-(@) Open its table, start the editing mode and see that columns with time dependent parameters are visible now.
+(@) In your geopackage select the layer "timml Well:DewateringWell".
+(@) Open its table, start the editing mode and see that columns with time dependent parameters are visible now.
(@) Switch off the Steady State abstraction, make discharge = 0.
(@) Switch on the Transient abstraction: time_start = 2, time_end = 30, discharge_transient = 2000, caisson_radius = 1.
(@) Save your changes.
-(@) Compute your transient model in the tab *Compute*.
+(@) Compute your transient model in the tab *Compute*.
In case of a successful calculation a new layer is added to your Vector output: case-Rijsenhout-ttim Observation:Piezometers. This layer contains transient data which is indicated with the clock icon right from the layer name. There are 2 ways to visualize these calculated heads in the observations point a) animation over time and b) timeseries at location.
@@ -485,39 +485,39 @@ In case of a successful calculation a new layer is added to your Vector output:
(@) In the Temporal Controller panel click the green play button ({width=14%}). Navigation buttons appear.
(@) Increase the *Step* to 16 hours and start the animation with a click on the *Play* button ({width=4%}).
-Only the first few days you see most prominent drawdown. Let's now display the timeseries at the point of your mouse with data from the mesh.
+Only the first few days you see most prominent drawdown. Let's now display the timeseries at the point of your mouse with data from the mesh.
-(@) Go to the iMOD Time Series panel.
-(@) Be sure the layer "case-Rijsenhout-head_layer_0" is selected.
-(@) For "variable:" select *head* and for "layers:" select *0*.
-(@) Click the button *Select Points*, your mouse changes into a . Be sure *Update on Selection* is checked and hover over the mesh. The graph shows the timeseries at the location of your mouse.
-(@) Deselect the checkbox *Update on Selection*.
-(@) Click the button *Select Points* (your mouse becomes a  and with the left mouse button select 3 points ad random.
-(@) Finally click the *Plot* button and the 3 timeseries are added to the chart.
+(@) Go to the iMOD Time Series panel.
+(@) Be sure the layer "case-Rijsenhout-head_layer_0" is selected.
+(@) For "variable:" select *head* and for "layers:" select *0*.
+(@) Click the button *Select Points*, your mouse changes into a . Be sure *Update on Selection* is checked and hover over the mesh. The graph shows the timeseries at the location of your mouse.
+(@) Deselect the checkbox *Update on Selection*.
+(@) Click the button *Select Points* (your mouse becomes a  and with the left mouse button select 3 points ad random.
+(@) Finally click the *Plot* button and the 3 timeseries are added to the chart.
-In the same way you can display the time series from the Observations.
+In the same way you can display the time series from the Observations.
-(@) In the iMOD time series panel select the layer "case-Rijsenhout-ttim Observation:Piezometers".
-(@) For "ID column:" select *label* and for "Variable:" select *head_layer0*. Don't forget to **deselect *fid***.
-(@) Click the button *Select Points* and draw a box with your mouse to select one or more observation points.
-(@) Click the *Plot* button and the selected timeseries are added to the chart.
+(@) In the iMOD time series panel select the layer "case-Rijsenhout-ttim Observation:Piezometers".
+(@) For "ID column:" select *label* and for "Variable:" select *head_layer0*. Don't forget to **deselect *fid***.
+(@) Click the button *Select Points* and draw a box with your mouse to select one or more observation points.
+(@) Click the *Plot* button and the selected timeseries are added to the chart.
{width=100% #fig-figure_HeadTimeSeries01}
-Before closing the display we change the line properties and save the graph as a PNG file.
+Before closing the display we change the line properties and save the graph as a PNG file.
(@) Select the single time series for the observation point closest to your dewatering well. The color of the line is now visible in the field next to the button *Line Color*.
(@) Change the color to red.
(@) Select the timeseries with the most shallow timeseries and change the color to green.
(@) Check the box *Draw markers*.
-(@) Move your mouse to the display and click your right mouse button. Explore the different options.
+(@) Move your mouse to the display and click your right mouse button. Explore the different options.
(@) From this menu (or from the display) select the function *Export...*.
-(@) Set the export format to PNG and export the file and don't forget to close the *Export* window.
+(@) Set the export format to PNG and export the file and don't forget to close the *Export* window.
## Export your Tim model to Python
-For scenario calculation or sensitivity analysis you probably want to switch from QGIS to Python for efficiency reasons. In this tutorial we only show you the first step: export your model to a Python file (*.py). Your Tim model is just a set of elements (points, lines, polygons) and its parameters so the Python file is not very large.
+For scenario calculation or sensitivity analysis you probably want to switch from QGIS to Python for efficiency reasons. In this tutorial we only show you the first step: export your model to a Python file (*.py). Your Tim model is just a set of elements (points, lines, polygons) and its parameters so the Python file is not very large.
-(@) In QGIS-Tim go to the tab *GeoPackage*.
+(@) In QGIS-Tim go to the tab *GeoPackage*.
(@) Click the button *Convert GeoPackage to Python script* and you can save the *.py file wherever you like.
(@) Check the content of the file with you text editor (e.g. Notepad).
diff --git a/docs/tutorial_TheHague.qmd b/docs/tutorial_TheHague.qmd
index e4e124e..a2bc6fa 100644
--- a/docs/tutorial_TheHague.qmd
+++ b/docs/tutorial_TheHague.qmd
@@ -42,15 +42,15 @@ ol.w {list-style-type: inherit;}
-->
---
-title: "Tutorial The Hague building pit"
+title: "Building Pit The Hague "
---
## Requirements
-- Be sure QGIS version 3.22.00 or higher is installed.
-- Be sure the ``gistim`` Python package is installed (see [installation](install.qmd) for instructions).
+- Be sure QGIS version 3.22.00 or higher is installed.
+- Be sure the ``gistim`` Python package is installed (see [installation](install.qmd) for instructions).
- Download the tutorial material. Follow this link.
-- Installation of the QGIS-Tim plugin and the MOD plugin is part of this Tutorial. The necessary ZIP files are included in the tutorial material.
-- Internet connections is optional during this Tutorial. It is only required for installation of additional plugins and the use of an online topographic background map.
+- Installation of the QGIS-Tim plugin and the MOD plugin is part of this Tutorial. The necessary ZIP files are included in the tutorial material.
+- Internet connections is optional during this Tutorial. It is only required for installation of additional plugins and the use of an online topographic background map.
## Description
In this tutorial, you will learn how to:
@@ -68,12 +68,12 @@ Calculation of a pumping well extraction.
In the old city centre of The Hague a building plan with a parking basement is constructed.
The dimensions of the building pit are **34.2 m * 83.2 m.** The area is shown in @fig-TheHagueCityCentre.
-{width=50% #fig-TheHagueCityCentre fig-align="left"}
+{width=50% #fig-TheHagueCityCentre fig-align="left"}
-To carry out the construction in dry conditions, the groundwater level needs to be lowered by **3.5 m.** The required drainage time is 4 months.
+To carry out the construction in dry conditions, the groundwater level needs to be lowered by **3.5 m.** The required drainage time is 4 months.
The environment is vulnerable with valuable buildings and monuments (see @fig-TheMonumentalEnvironment) that have a shallow foundation. The subsoil consists of sandy layers which locally contain shallow peat/clay/silt layers and deeper layers of clay/silt. The builder contractor must use a sheet pile wall that needs to be installed in a water-retarding layer (aquitard).
-{width=50% #fig-TheMonumentalEnvironment fig-align="left"}
+{width=50% #fig-TheMonumentalEnvironment fig-align="left"}
**The questions to be elaborated in the hydrogeological advice is:**
*"How deep should a sheet pile wall be placed to prevent damage? Can risks arise from a leak?"*
@@ -82,13 +82,13 @@ Applications of a sheet pile wall (with a hydraulic resistance of 100 days) and
The soil profile has been investigated with borings and CPTs. An example of a CPT is shown in @fig-CPT01.
-{width=70% #fig-CPT01 fig-align="left"}
+{width=70% #fig-CPT01 fig-align="left"}
The interpretation of the CPT and borings gives the following soil layer scheme:
- The soil mainly consists of fine sand in Holocene layers (max depth NAP -15 m) and coarse sand in deeper Pleistocene layers.
-- Just below the groundwater level (NAP -0.5 m), there is an approximately 1.5 m thick peat/clay/silt layer (not in CPT).
-- The layer between NAP –5.5 m and –8.5 m consists of clay lenses. There is little certainty about the resistance of that layer.
+- Just below the groundwater level (NAP -0.5 m), there is an approximately 1.5 m thick peat/clay/silt layer (not in CPT).
+- The layer between NAP –5.5 m and –8.5 m consists of clay lenses. There is little certainty about the resistance of that layer.
- Between NAP –13 m and –15 m there is a clay/silt layer.
Based on experience in this region, the following geohydrological schematization is known (@tbl-GeohydrologicalModelSubsoil). The layer depth is indicated relative to the groundwater level of NAP-0.5m, because only that part is relevant for our model. In @tbl-GeohydrologicalModelSubsoil c = vertical hydraulic resistance, k= horizontal permeability, S’= specific storage: :
@@ -120,78 +120,78 @@ Tasks:
- What are the risks? Look at the groundwater drawdown outside the construction pit in situations with different sheet pile depth.
- Also check for point leaks. Model these leaks with a hole in the sheet pile in the southwest corner of the area.
-## Getting Started
+## Getting Started
Let's first configure the gistim python installation again, to be sure that QGIS can find the gistim software.
-(@) Open the Deltaforge prompt (search in Windows Start for "Deltaforge Prompt"). A black window pops up.
+(@) Open the Deltaforge prompt (search in Windows Start for "Deltaforge Prompt"). A black window pops up.
(@) In this window type ``python -m gistim configure`` and press ENTER.
-(@) You can close this window now.
+(@) You can close this window now.
(@) Launch QGIS from your START menu, from your desktop or click on …\\QGIS3.28.0\\bin\\qgis-bin.exe.
> **Intermezzo:** *QGIS language settings*
->
-> Perhaps your QGIS was installed in another language than English. Because the Tutorial refers to the English version, let’s change to English.
->
-> a. From the main menu click on *Settings* and select *Options* (e.g. in Dutch *Extra* and *Opties*).
-> a. In the new window go to the *General section* (Dutch: *Algemeen*) on the left.
-> a. Check the box to allow *Override System Locale* (Dutch: *Landinstellingen negeren*) and expand this sub menu.
-> a. From the drop-down menu *“User interface translation”* (Dutch: *Vertaling gebruikers-interface*) select *American English* and click *OK*.
+>
+> Perhaps your QGIS was installed in another language than English. Because the Tutorial refers to the English version, let’s change to English.
+>
+> a. From the main menu click on *Settings* and select *Options* (e.g. in Dutch *Extra* and *Opties*).
+> a. In the new window go to the *General section* (Dutch: *Algemeen*) on the left.
+> a. Check the box to allow *Override System Locale* (Dutch: *Landinstellingen negeren*) and expand this sub menu.
+> a. From the drop-down menu *“User interface translation”* (Dutch: *Vertaling gebruikers-interface*) select *American English* and click *OK*.
> a. Close QGIS and open it again to activate your language change.
We start with the creation of a new QGIS project.
-(@) From the main menu click on *Project* and select *New*.
+(@) From the main menu click on *Project* and select *New*.
The case in this tutorial is located in The Netherlands, so next we select the appropriate projection.
-(@) From the main menu click on *Project* and select *Properties*.
+(@) From the main menu click on *Project* and select *Properties*.
(@) In the *Properties* window select the category *CRS*, search for “EPSG:28992” and you find “Amersfoort / RD New”. Select this option and click the *Apply* button, followed by the *OK* button to close the window.
> In case your work is mostly in The Netherlands and in the “Amersfoort / RD New” projection, consider making this your default projection.
->
+>
> - From the main menu click on *Settings* and select *Options...*.
> - In the section *CRS and Transforms* select *CRS (handling)*, pick the radio button *Use a default CRS* and select “EPSG:28992 -Amersfoort / RD New”.
> - Click OK.
> - Close this window.
### Install plugins
-This is the moment to download/import four plugins needed for this tutorial. This is the list:
+This is the moment to download/import four plugins needed for this tutorial. This is the list:
- the QGIS-Tim plugin. The development version, imported from a ZIP file.
- the iMOD plugin. The development version, imported from a ZIP file.
-- the Value Tool. The official version, installed via the Plugin Manager of QGIS (internet connection required).
+- the Value Tool. The official version, installed via the Plugin Manager of QGIS (internet connection required).
- the PDOK plugin. The official version, installed via the Plugin Manager of QGIS (internet connection required). This plugin gives access to a large database from which we will load the topographic maps and use the navigation option.
-(@) Go to *Plugins* from the main menu and select *Manage and Install Plugins…* to open the plugin window.
+(@) Go to *Plugins* from the main menu and select *Manage and Install Plugins…* to open the plugin window.
(@) On the left section select *Install from ZIP*.
-(@) Click the *Browse* button ({width=4%}) and from the tutorial dataset select the ZIP file "QGIS-Tim_Tutorial-TheHague\\QGIS-iMOD-plugin.zip".
+(@) Click the *Browse* button ({width=4%}) and from the tutorial dataset select the ZIP file "QGIS-Tim_Tutorial-TheHague\\QGIS-iMOD-plugin.zip".
(@) Click *Install Plugin*.
-(@) In the same way, install the QGIS-Tim plugin using the ZIP file "QGIS-Tim_Tutorial-TheHague\\QGIS-Tim-plugin.zip".
+(@) In the same way, install the QGIS-Tim plugin using the ZIP file "QGIS-Tim_Tutorial-TheHague\\QGIS-Tim-plugin.zip".
If you have an internet connection continue with the installation of the next two plugins from the QGIS plugin library.
-(@) From the left section, select the group *All* to see all available plugins.
-(@) Search for “Value Tool” and install it.
-(@) Search for “PDOK services plugin” and install it.
+(@) From the left section, select the group *All* to see all available plugins.
+(@) Search for “Value Tool” and install it.
+(@) Search for “PDOK services plugin” and install it.
(@) Make sure that under *Plugins > Manage and Install Plugins > Installed* now the 4 added plugins are checked.
-(@) *Close* the Plugins window.
+(@) *Close* the Plugins window.
-See in the toolbar section of QGIS that the plugins are installed:
+See in the toolbar section of QGIS that the plugins are installed:
-- iMOD Toolbar {width=30%}
+- iMOD Toolbar {width=30%}
- QGIS-Tim {width=4%}
- Value Tool {width=4%}
- PDOK Services Plugin {width=4%}
-Further in this Tutorial we will use some default toolbars that might be hidden at the moment. Let's check that and unhide if necessary.
+Further in this Tutorial we will use some default toolbars that might be hidden at the moment. Let's check that and unhide if necessary.
(@) Select *View* from the main menu and choose *Panels* and be sure these two toolbars are checked:
- Layers
- Browser.
(@) Select *View* from the main menu and choose *Toolbars* and be sure these three toolbars are checked:
- Advanced Digitizing Toolbar
- Snapping Toolbar
- Attributes Toolbar.
### Prepare your project
-For navigation purposes, let's load a topographic map for The Netherlands from the online PDOK database.
+For navigation purposes, let's load a topographic map for The Netherlands from the online PDOK database.
-> No internet connection? Follow the next steps to import a simple PNG file as a background.
+> No internet connection? Follow the next steps to import a simple PNG file as a background.
>
> - Go to *Layer* in the main menu, go to *Add layer* and select *Add Raster layer*.
> - Use the *Browse* button ({width=4%}) and from the tutorial material select “…\\QGIS-Tim_Tutorial-TheHague\\TopographicMapTheHague.png”.
@@ -199,10 +199,10 @@ For navigation purposes, let's load a topographic map for The Netherlands from t
> - If you do not see the map, select the layer "TopographicMapTheHague", click your right mouse button and select "Zoom to Layer(s)".
> - Continue after step @continuewithoutinternet.
-(@) If you do have an internet connection click on the PDOK plugin button ({width=4%}) to open the "PDOK Services Plugin" window.
-(@) From the tab *PDOK Services* search for "pastel" and you will find a WMTS type layer called "BRT Achtergrondkaart WMTS".
+(@) If you do have an internet connection click on the PDOK plugin button ({width=4%}) to open the "PDOK Services Plugin" window.
+(@) From the tab *PDOK Services* search for "pastel" and you will find a WMTS type layer called "BRT Achtergrondkaart WMTS".
(@) Select the layer.
-(@) In the section "laag toevoegen" click the button *Onder*.
+(@) In the section "laag toevoegen" click the button *Onder*.
(@) Close the PDOK window.
Our project area is in the centre of the city of The Hague so let's navigate to that city using the PDOK plugin.
@@ -220,7 +220,7 @@ Let's now open a shape file containing the circumference of the building locatio
(@) In the *Layers* panel on the left, select the layer "building_pit".
(@) Click your right mouse button and from the menu select *Properties*.
-(@) In the new window go to the section *Symbology* on the left and try to pick a polygon style with only a contour color.
+(@) In the new window go to the section *Symbology* on the left and try to pick a polygon style with only a contour color.
(@) Click on *OK* a to save and close the window.
Let's save this project to be able to return to it later or in case of a crash of QGIS.
@@ -228,47 +228,47 @@ Let's save this project to be able to return to it later or in case of a crash o
(@) Go to *Project* in the main menu, select *Save As* and select a folder and a file name for your project, e.g. “…\\QGIS-Tim_Tutorial-TheHague\\TheHague.qgz”
### Open the QGIS-Tim panel
-Now we are ready to activate the QGIS-Tim plugin.
+Now we are ready to activate the QGIS-Tim plugin.
(@) Click on the QGIS-Tim plugin button ({width=4%}) and the QGIS-Tim panel appears.
-(@) Go to the tab *GeoPackage*.
Here we will create an empty database (geopackage) to store all elements and parameters for the model.
+(@) Go to the tab *GeoPackage*.
Here we will create an empty database (geopackage) to store all elements and parameters for the model.
(@) Click the *New* button to create the GeoPackage and save it for instance in the folder with your tutorial data, e.g. "..\\QGIS-Tim_Tutorial-TheHague\\case-TheHague.gpkg".
Your window looks like in @fig-Panel-QGIS-Tim_v02.
{width=50% #fig-Panel-QGIS-Tim_v02 fig-align="left"}
-(@) Check in the *Layers* panel on the left that your new geopackage is added as a group.
A sub group **timml** for the steady state model input, the sub group **ttim** for the transient model input and a series of output formats (vector/mesh/raster).
+(@) Check in the *Layers* panel on the left that your new geopackage is added as a group.
A sub group **timml** for the steady state model input, the sub group **ttim** for the transient model input and a series of output formats (vector/mesh/raster).
-If you had no introduction to the Tim plugin, read the Intermezzo below for a general explanation of the components.
+If you had no introduction to the Tim plugin, read the Intermezzo below for a general explanation of the components.
-> **Intermezzo:** *introduction Tabs on the Tim panel*
->
-> a. GeoPackage: an overview of the elements in your geopackage. In case you switch to transient modelling, an extra column with *ttim* elements is added.
-> a. Elements: a list of 14 Tim elements from which you can build your model.
-> a. Compute: here you can define your domain and cell size, decide if your model is transient or not and change the output name.
-> a. Extract: open an existing 3d geohydrological model (NC file) and extract the data for your project area.
+> **Intermezzo:** *introduction Tabs on the Tim panel*
+>
+> a. GeoPackage: an overview of the elements in your geopackage. In case you switch to transient modelling, an extra column with *ttim* elements is added.
+> a. Elements: a list of 14 Tim elements from which you can build your model.
+> a. Compute: here you can define your domain and cell size, decide if your model is transient or not and change the output name.
+> a. Extract: open an existing 3d geohydrological model (NC file) and extract the data for your project area.
Let's save this project to be able to return to it later or in case of a crash of QGIS.
(@) Go to *Project* in the main menu, select *Save As* and select a folder and a file name for your project, e.g. “…\\QGIS-Tim_Tutorial-TheHague\\TheHague.qgz”
## Checking available data from LHM
-We will now research what possible layer info is available from the LHM database (Landelijk hydrologisch model). Information from LHM is made available in the provided file “…\\QGIS-Tim_Tutorial-TheHague\\LHM4.1-ondergrondmodel.nc". The file contains rasterdata of both the 3D geological layering within The Netherlands and the corresponding geological parameterization (permeability, resistance, thickness).
+We will now research what possible layer info is available from the LHM database (Landelijk hydrologisch model). Information from LHM is made available in the provided file “…\\QGIS-Tim_Tutorial-TheHague\\LHM4.1-ondergrondmodel.nc". The file contains rasterdata of both the 3D geological layering within The Netherlands and the corresponding geological parameterization (permeability, resistance, thickness).
-(@) In the QGIS-Tim panel go to the *Extract* tab.
-(@) Open the file "..\\QGIS-Tim_Tutorial-TheHague\\LHM4.1-ondergrondmodel.nc" (this may take a minute because the file is 350 MB).
-(@) Click the button *Select by Polygon* and draw a large polygon (> 250x250m) to select the area around the building location between *Parkstraat* and *Oranjestraat* (see @fig-TheHagueCityCentre). Left-click several times and close the polygon with a right click.
-(@) Then click the *Extract* button and save the extracted subsoil data to a CSV file, e.g. "Ex2DH.csv".
+(@) In the QGIS-Tim panel go to the *Extract* tab.
+(@) Open the file "..\\QGIS-Tim_Tutorial-TheHague\\LHM4.1-ondergrondmodel.nc" (this may take a minute because the file is 350 MB).
+(@) Click the button *Select by Polygon* and draw a large polygon (> 250x250m) to select the area around the building location between *Parkstraat* and *Oranjestraat* (see @fig-TheHagueCityCentre). Left-click several times and close the polygon with a right click.
+(@) Then click the *Extract* button and save the extracted subsoil data to a CSV file, e.g. "Ex2DH.csv".
(@) A separate python window is opened and the selection is performed. You can follow the progress in your tool bar. There is no message indicating the process has ended. Just ignore the black screen and continue with the next step.
-> {width=90% #fig-figure_ExtractionError}
QGIS displays this error message in case your polygon is too small. In that case redraw your polygon for a larger area.
+> {width=90% #fig-figure_ExtractionError}
QGIS displays this error message in case your polygon is too small. In that case redraw your polygon for a larger area.
-Automatically 2 CSV files will be created. One with average values for bottom and top of layers and aquifer permeabilities and aquitard resistances. The other contains statistics on the data.
+Automatically 2 CSV files will be created. One with average values for bottom and top of layers and aquifer permeabilities and aquitard resistances. The other contains statistics on the data.
-(@) From a file manager on your laptop (e.g. Explorer) open the CSV files.
+(@) From a file manager on your laptop (e.g. Explorer) open the CSV files.
-An example of values from the first CSV are presented in @tbl-SubsoilCharacteristicsLHM. In each LHM layer the Aquitard information is positioned on top of Aquifer.
+An example of values from the first CSV are presented in @tbl-SubsoilCharacteristicsLHM. In each LHM layer the Aquitard information is positioned on top of Aquifer.
| fid
-| layer
-| aquitard
c| aquitard
npor| aquitard
s| aquifer
k| aquifer
npor| aquifer
s| aquifer
top| aquifer
bottom| semiconf
top| semiconf
head|
|---|---|:--:|:--:|:--:|:--:|:--:|:--:|:-:|:--:|:--:|:--:|
@@ -289,77 +289,77 @@ For marine deposits at -12 a value of 50 to 100 days per meter layer thickness i
## Start your Tim model
We are now ready to define our first steady state model by parameterizing our Aquifer.
-(@) From the *Layers* panel, select the layer *timml Aquifer:Aquifer*.
+(@) From the *Layers* panel, select the layer *timml Aquifer:Aquifer*.
(@) Click your right mouse button and from the menu select *Open Attribute table* ({width=3%}). Search for the same icon somewhere on the *Attributes Toolbar*.
(@) Switch to the editing mode in the table with a click on ({width=4%}) *Toggle Editing*.
-(@) Then by clicking 3 times on ({width=4%}) *Add feature* you will see new rows are added.
-(@) Now you are able to fill in the desired values shown in @fig-AquiferAttrTable_v01.
-**NB! Do not fill in column FID. 'Autogenerate' will take care.**
+(@) Then by clicking 3 times on ({width=4%}) *Add feature* you will see new rows are added.
+(@) Now you are able to fill in the desired values shown in @fig-AquiferAttrTable_v01.
+**NB! Do not fill in column FID. 'Autogenerate' will take care.**
-{width=70% #fig-AquiferAttrTable_v01}
+{width=70% #fig-AquiferAttrTable_v01}
-**NB!** Closing this table window with the {width=2.5%} will not save the filled in data or stop the editing mode!
+**NB!** Closing this table window with the {width=2.5%} will not save the filled in data or stop the editing mode!
-(@) Click on the *Save-Edits* button ({width=3%}) to save you data during the process or click on the *Toggle Editing Mode* button ({width=4%}) to stop edditing and QGIS askes you if your changes should be saved.
+(@) Click on the *Save-Edits* button ({width=3%}) to save you data during the process or click on the *Toggle Editing Mode* button ({width=4%}) to stop edditing and QGIS askes you if your changes should be saved.
(@) Now you can close the Attribute Table window.
-In the QGIS-Tim Window we now introduce the following elements for ground water modelling in the QGIS-Tim tab *Element*:
+In the QGIS-Tim Window we now introduce the following elements for ground water modelling in the QGIS-Tim tab *Element*:
- Leaky Line Doublet
- Well
### Adding a Leaky Line Doublet
-(@) In the QGIS-Tim window go to tab *Elements* and select the button *Leaky Line Doublet*.
-(@) Fill in the name of the layer in the pop-up panel, e.g. 'sheet_pile', and click *OK* and this new Layer is added to the *Layers* panel.
+(@) In the QGIS-Tim window go to tab *Elements* and select the button *Leaky Line Doublet*.
+(@) Fill in the name of the layer in the pop-up panel, e.g. 'sheet_pile', and click *OK* and this new Layer is added to the *Layers* panel.
-First, we will *draw* the location of the sheet pile before *adding* the parameter values.
+First, we will *draw* the location of the sheet pile before *adding* the parameter values.
-(@) In the *Layers* panel right click on the layer {width=4%} *timml Leaky Line Doublet:sheet_pile* and start the editing mode by clicking the *Toggle Editing* button ({width=4%}).
+(@) In the *Layers* panel right click on the layer {width=4%} *timml Leaky Line Doublet:sheet_pile* and start the editing mode by clicking the *Toggle Editing* button ({width=4%}).
-Two important remarks before drawing the sheet pile:
+Two important remarks before drawing the sheet pile:
-- Start drawing at the most southern corner of the building pit to make a variation for an extra purpose later on in this tutorial.
-- In order to make a closed area surrounded with the sheet piles, the first and last point must have the same coordinate and intersect. For this we must activate the Snapping option in QGIS.
+- Start drawing at the most southern corner of the building pit to make a variation for an extra purpose later on in this tutorial.
+- In order to make a closed area surrounded with the sheet piles, the first and last point must have the same coordinate and intersect. For this we must activate the Snapping option in QGIS.
-(@) In the main menu go to *Project* and select *Snapping Options..*.
-(@) From the new toolbar (see @fig-button_ToolbarSnappingSettings), click on the *Enable Snapping* button ({width=3%}).
+(@) In the main menu go to *Project* and select *Snapping Options..*.
+(@) From the new toolbar (see @fig-button_ToolbarSnappingSettings), click on the *Enable Snapping* button ({width=3%}).
-{width=90% #fig-button_ToolbarSnappingSettings}
+{width=90% #fig-button_ToolbarSnappingSettings}
-(@) Then select the *Add line feature* icon ({width=3%}) to start the drawing of the sheet pile.
-(@) Draw the sheet pile around the building pit (start south corner) with your left mouse button and close the feature on the first point.
+(@) Then select the *Add line feature* icon ({width=3%}) to start the drawing of the sheet pile.
+(@) Draw the sheet pile around the building pit (start south corner) with your left mouse button and close the feature on the first point.
(@) In the pop-up window fill in the resistance (1000 d, see \* below) and layer number (layer=0, see \*\* below) as in @fig-LeakylinedoubletFeatureAttributes and click OK to close the window.
-{width=60% #fig-LeakylinedoubletFeatureAttributes fig-align="left"}
+{width=60% #fig-LeakylinedoubletFeatureAttributes fig-align="left"}
> Remarks to the provided values:
> \* Because of a known issue in the TimML software, until this issue is resolved, the value of the resistance must be multiplied by the permeability. For normal sheet pile walls a resistance value of 100-250 d for interlock leakage is a good choice. In this case the chosen value needs to be increased by multiplying with the applicable aquifer permeability in the effected layer (R=100 d becomes R=10*100 d).
> \*\* In the real-world counting starts with 1. However, Tim is programmed in Python and in Python counting starts with 0. You will get used to it.
-(@) Close the editing mode with a click on ({width=4%}) and you are asked to save your changes.
-(@) Open the attribute table for *Leaky Line Doublet* (click {width=3%} or press F6) to check hydraulic resistance values for the sheet pile walls.
+(@) Close the editing mode with a click on ({width=4%}) and you are asked to save your changes.
+(@) Open the attribute table for *Leaky Line Doublet* (click {width=3%} or press F6) to check hydraulic resistance values for the sheet pile walls.
### Adding a Well
-(@) In GIS-Tim go to the tab *Elements* and select the *Well* element.
-(@) A name for the element can be given in the pop-up panel, e.g. *pumping_wells*.
-(@) In the *Layers* panel right click on the layer {width=3%} *timml Well:pumping_wells* and start the editing mode by clicking the *Toggle Editing* button ({width=4%}).
+(@) In GIS-Tim go to the tab *Elements* and select the *Well* element.
+(@) A name for the element can be given in the pop-up panel, e.g. *pumping_wells*.
+(@) In the *Layers* panel right click on the layer {width=3%} *timml Well:pumping_wells* and start the editing mode by clicking the *Toggle Editing* button ({width=4%}).
(@) Next click on *Add point feature* ({width=4%}).
-(@) With you left mouse button add the first well location similar to point 1 in @fig-SchematizationTheHague.
+(@) With you left mouse button add the first well location similar to point 1 in @fig-SchematizationTheHague.
(@) In the pop-up window, fill in the well parameters discharge (50 m3/d), radius (0.1 m), resistance (1 d) and layer (0).
-(@) Now add the other 7 wells locations in a fast way: do not import parameters with every single point. We will do that later. Just click *OK* on each *Feature Attribute* window.
+(@) Now add the other 7 wells locations in a fast way: do not import parameters with every single point. We will do that later. Just click *OK* on each *Feature Attribute* window.
-{width=60% #fig-SchematizationTheHague fig-align="left"}
+{width=60% #fig-SchematizationTheHague fig-align="left"}
-(@) In the *Layers* panel right select layer *timml Well:pumping_wells* and open the *Attribute Table (F6)*.
+(@) In the *Layers* panel right select layer *timml Well:pumping_wells* and open the *Attribute Table (F6)*.
(@) Start the editing mode and fill in the values shown in @fig-WellAttributeTable. Don't forget the column 'Label'.
The discharge value in the table is a first guess to be adjusted later to desired level of -3.5 m groundwater lowering in the building pit.
-(@) Stop the editing mode, save your work and close the window.
-(@) Select layer *timml Well:pumping_wells* again, click right and from the menu select *Show labels*.
+(@) Stop the editing mode, save your work and close the window.
+(@) Select layer *timml Well:pumping_wells* again, click right and from the menu select *Show labels*.
-{width=100% #fig-WellAttributeTable}
+{width=100% #fig-WellAttributeTable}
-When you finished the input of the wells and their parameters, you have to look back to the combination of wells and sheet pile wall. The reason is that near a well the flow is strong and therefore the flow through the wall is increased. This may lead to extremities in the calculation if wall segments are chosen too large. Tim uses control points but they are set at regular distances on the doublet (the number of control points is 1+order). It is necessary to divide the wall in smaller sections near well locations. The length of each section should be chosen equal to approximately the distance of the wall to the nearest well.
+When you finished the input of the wells and their parameters, you have to look back to the combination of wells and sheet pile wall. The reason is that near a well the flow is strong and therefore the flow through the wall is increased. This may lead to extremities in the calculation if wall segments are chosen too large. Tim uses control points but they are set at regular distances on the doublet (the number of control points is 1+order). It is necessary to divide the wall in smaller sections near well locations. The length of each section should be chosen equal to approximately the distance of the wall to the nearest well.
(@) Select the Leaky Line Doublet in the *Layers* panel and start the editing mode.
(@) Activate the Vertex Tool ({width=4%}) and then hoover along the sheet pile line in your drawing. The line element near your pointer lights up.
@@ -367,50 +367,50 @@ When you finished the input of the wells and their parameters, you have to look
@fig-figure_SheetPileExtraPoints shows the result of positioned extra points in the line.
-{width=60% #fig-figure_SheetPileExtraPoints fig-align="left"}
+{width=60% #fig-figure_SheetPileExtraPoints fig-align="left"}
## Computing the groundwater head drawdown
(@) Zoom in or out to desired domain for which you want to see the model results.
-(@) In the QGIS-Tim panel select the tab *Compute*.
-(@) Select the button *Set to current extent* to define the Domain.
-(@) Grid spacing will follow automatically but for now make the results mesh more dense by changing "Grid spacing" to 3.00 m.
-(@) In the "Output" section give the name of the file where you want to store the results.
-(@) For contouring, select the check box "Auto-generate contours".
-(@) Set the increment for contouring to a proper value: 0.5, 0.25 or 0.1 as applicable for your study.
-(@) Press the *Compute* button to have the program perform the calculations.
+(@) In the QGIS-Tim panel select the tab *Compute*.
+(@) Select the button *Set to current extent* to define the Domain.
+(@) Grid spacing will follow automatically but for now make the results mesh more dense by changing "Grid spacing" to 3.00 m.
+(@) In the "Output" section give the name of the file where you want to store the results.
+(@) For contouring, select the check box "Auto-generate contours".
+(@) Set the increment for contouring to a proper value: 0.5, 0.25 or 0.1 as applicable for your study.
+(@) Press the *Compute* button to have the program perform the calculations.
-A black Python.exe window pops up indicating that the TIM calculation started on the background. You can ignore this window but keep it open. Of course you van minimize it. If the calculation was completed successful, you will see this echo in QGIS.
-{width=80%}.
+A black Python.exe window pops up indicating that the TIM calculation started on the background. You can ignore this window but keep it open. Of course you van minimize it. If the calculation was completed successful, you will see this echo in QGIS.
+{width=80%}.
## Studying output results
-After the calculation you see that the result is automatically added to the *Layers* panel, probably called "case-TheHague output". Results are presented as Mesh, Raster and Vector data. Contours are saved as vector.
Although these layers / groups are checked, the data is not visible. That is because the geopackage was added last, and QGIS adds layers at the end of the list.
Let's move the layer "pastel" to the background.
+After the calculation you see that the result is automatically added to the *Layers* panel, probably called "case-TheHague output". Results are presented as Mesh, Raster and Vector data. Contours are saved as vector.
Although these layers / groups are checked, the data is not visible. That is because the geopackage was added last, and QGIS adds layers at the end of the list.
Let's move the layer "pastel" to the background.
(@) Select the layer "pastel" and drag it with your left mouse button to the bottom of the list of layers.
(@) Uncheck mesh and raster to only visualize contours on the base map.
(@) Uncheck contour lines of layer 1 and 2 to get only the result for layer 0 on screen as in @fig-DrawdownConstructionPit_v01.
-{width=70% #fig-DrawdownConstructionPit_v01 fig-align="left"}
+{width=70% #fig-DrawdownConstructionPit_v01 fig-align="left"}
-The lowering in the building pit is shown to be around -6.20 m. To get it around -3.50 m as demanded, the well discharge should be decreased accordingly. An extra calculation with half the flow per well will be sufficient.
+The lowering in the building pit is shown to be around -6.20 m. To get it around -3.50 m as demanded, the well discharge should be decreased accordingly. An extra calculation with half the flow per well will be sufficient.
(@) Go to the well attribute table, adjust the flow to 60% (30 m3/d per well) and compute the model again.
### Adding observations wells
The city authorities that perform quality checks on the effect of construction projects insisted on the installation of some piezometers to assure reduction of risks for surrounding old monumental structures.
-(@) In the QGIS-Tim panel go to the *Elements* tab and select the element *Observation*.
-(@) Give a name in the pop-up panel, e.g. "observations".
-(@) Then go to the *Layers* panel and select the layer {width=3%}*timml Observation:observations*.
-(@) Go to the toolbar with drawing options and activate Toggle editing ({width=4%}).
+(@) In the QGIS-Tim panel go to the *Elements* tab and select the element *Observation*.
+(@) Give a name in the pop-up panel, e.g. "observations".
+(@) Then go to the *Layers* panel and select the layer {width=3%}*timml Observation:observations*.
+(@) Go to the toolbar with drawing options and activate Toggle editing ({width=4%}).
(@) Go right to the Add point feature ({width=3%}) and drop some piezometers in the drawing.
We propose to use 5 points (see @fig-SchematizationTheHague): 1 point in the centre of the building pit, 1 near the lower corner outside, 1 outside near the middle of the eastern sheet pile section and just 2 near street corners.
(@) Compute the model again.
Results of calculations at the observation locations are presented as *Vector* data in the *Layers* panel.
-(@) Uncheck/check layers in order to display observation results only.
+(@) Uncheck/check layers in order to display observation results only.
-By default, the label at each observation location is the calculated head for layer 0. Perhaps you see a second number near the location. For your information: this is the "location number" label belonging to the model input in layer *timml Observation:observations*.
+By default, the label at each observation location is the calculated head for layer 0. Perhaps you see a second number near the location. For your information: this is the "location number" label belonging to the model input in layer *timml Observation:observations*.
We can observe that the lowering of groundwater around the building pit is quite high due to leakage of the sheet piles or leakage through the bottom clay layer with small resistance. Therefore, it is needed to improve the wall quality or the length. First we check the effect of the clay layer by studying a cross section.
@@ -424,22 +424,22 @@ We can observe that the lowering of groundwater around the building pit is quite
(@) Satisfied with your line? Click the button *Plot* to draw this layer in the cross section. Your screen might look like @fig-HeadsInCrosssection03.
(@) By using the *Export* button you can store results from the cross section in a CSV file.
-{width=100% #fig-HeadsInCrosssection03}
+{width=100% #fig-HeadsInCrosssection03}
## Making calculations with parameter variations or checking bandwidth
The authorities demand a drawdown effect of dewatering at a maximum of 0.10 m at surrounding buildings. This means that improvements for leakage control are needed but first we need to discover what parameter to focus on.
To check whether the wall resistance or the bottom resistance of the layer 1 (below the building pit) is more important we can make 2 variations; one with C-clay=200 d and one with R-wall=500 d.
-Of course you can change your model input, rerun the model and overwrite your model results. The next steps show you how to change the model input and save the results in separate *.gpkg and *.nc files.
+Of course you can change your model input, rerun the model and overwrite your model results. The next steps show you how to change the model input and save the results in separate *.gpkg and *.nc files.
(@) In the input group select layer *timml Aquifer:...*
-(@) Open the *Attrribute Table* (F6) and change the value for "aquitard_c" in layer 1 into 200 d.
+(@) Open the *Attrribute Table* (F6) and change the value for "aquitard_c" in layer 1 into 200 d.
(@) In the QGIS-Tim panel go to the tab *Compute* and change the name of the output, e.g. case-TheHague_v1.
(@) Click *Compute* to run variant 1.
-Check in the *Layers* panel and see that the results are not overwritten but added to the groups, e.g. layer *case-TheHague_v1-timml Observation:observations* is added to the group *Vector*.
+Check in the *Layers* panel and see that the results are not overwritten but added to the groups, e.g. layer *case-TheHague_v1-timml Observation:observations* is added to the group *Vector*.
-(@) Fill in your calculated heads at the observation locations in @tbl-CalculatedHeads2Variants or use Excel.
+(@) Fill in your calculated heads at the observation locations in @tbl-CalculatedHeads2Variants or use Excel.
|Observation
location |Default:
Cc=40d
Rw=100d*|{width=20%}
your value:|Variant 1:
Cc=200d
Rw=100d*|your value:|Variant 2:
Cc=40d
Rw=500d*|your value:|
|---|---|---|---|---|---|---|
@@ -454,11 +454,11 @@ Check in the *Layers* panel and see that the results are not overwritten but add
Let's now run Variant 2.
-(@) In layer *timml Aquifer:...* reset the value for "aquitard_c" in layer 1 to the default of 40 d.
-(@) In layer *timml Leaky Line Doublet:...* change the value for "resistance" into 5000 d (500x10).
+(@) In layer *timml Aquifer:...* reset the value for "aquitard_c" in layer 1 to the default of 40 d.
+(@) In layer *timml Leaky Line Doublet:...* change the value for "resistance" into 5000 d (500x10).
(@) In the QGIS-Tim panel go to the tab *Compute* and change the name of the output, e.g. case-TheHague_v2.
(@) Click *Compute* to run variant 2.
-(@) Fill in your calculated heads at the observation locations in the table above.
+(@) Fill in your calculated heads at the observation locations in the table above.
To get the same lowering in the building pit, in the second variation the well flow might be reduced to 33% (10m3/d per well). For the sheet pile wall, increasing the wall quality or decreasing interlock leakage doesn’t make a big difference.
We can conclude that the best investment during the phase of design would be to perform extra hydrogeological research, e.g. by making more cpt’s, borings or performing a pumping test.
@@ -466,39 +466,39 @@ We can conclude that the best investment during the phase of design would be to
## Sheet piles with extra depth
Suppose the best guess value of the clay layer resistance was right, then a mitigation measure for the effect of dewatering in the construction phase could be the installation of the wall to a deeper level where additional hydraulic resistance of 100d can be found at a depth of -13 m to -15 m NAP.
-If we want to create extra depth of the sheet pile we will have to introduce it in a deeper layer. There are 2 options to implement it in your model:
+If we want to create extra depth of the sheet pile we will have to introduce it in a deeper layer. There are 2 options to implement it in your model:
-- Add a copy of the geometry of the sheet pile wall to the existing *Leaky Line Doublet* shape and assign it to layer 1.
-- We recommend to create an extra *Leaky Line Doublet* element. In this case it is more easy to switch on/off this additional element in your sensitivity analysis.
+- Add a copy of the geometry of the sheet pile wall to the existing *Leaky Line Doublet* shape and assign it to layer 1.
+- We recommend to create an extra *Leaky Line Doublet* element. In this case it is more easy to switch on/off this additional element in your sensitivity analysis.
-How to copy the sheet pile wall to an extra *Leaky Line Doublet* element?
+How to copy the sheet pile wall to an extra *Leaky Line Doublet* element?
(@) In the QGIS-Tim panel go to the tab *Elements* and add a second *Leaky Line Doublet* and give it a name, e.g. "sheet_pile_L1"
-(@) Go to the tab *GeoPackage* and see that the element separately is added to the list. Here is can switch this element on / off for a calculation.
+(@) Go to the tab *GeoPackage* and see that the element separately is added to the list. Here is can switch this element on / off for a calculation.
(@) In the *Layers* panel select the new layer *timml Leaky Line Doublet:sheet_pile_l1*.
-(@) Open its *Attribute Table* (F6) and start the editing mode. The table is empty.
-(@) Also open the *Attribute Table* of the first *Leaky Line Doublet* and select the existing element.
+(@) Open its *Attribute Table* (F6) and start the editing mode. The table is empty.
+(@) Also open the *Attribute Table* of the first *Leaky Line Doublet* and select the existing element.
(@) Click on the *Copy* button ({width=3%}) in the source table to copy the selected row to the clipboard.
- {width=70%}
+ {width=70%}
(@) In the target table, paste it with the *Paste* button ({width=3%}) as a new layer.
(@) Assign this new sheet pile to layer=1.
-(@) Stop editing and save the new element.
-(@) Click *Compute* to start the computation again.
+(@) Stop editing and save the new element.
+(@) Click *Compute* to start the computation again.
The results are directly visible in the contours and cross-section again.
-We can conclude that the drawdown in the building pit increases with a factor of almost 2 (-7.48 m at the centre of the building pit). Therefore, we adjust the well flow to 3.95/7.48*30=15.54 m3/d per well.
+We can conclude that the drawdown in the building pit increases with a factor of almost 2 (-7.48 m at the centre of the building pit). Therefore, we adjust the well flow to 3.95/7.48*30=15.54 m3/d per well.
(@) Implement this change in the *timm Well* element and recalculate the model.
-{width=100% #fig-figure_HeadsInCrosssection04}
+{width=100% #fig-figure_HeadsInCrosssection04}
-We conclude that the drawdown of groundwater level around the building pit with deep sheet piles decreased significantly.
+We conclude that the drawdown of groundwater level around the building pit with deep sheet piles decreased significantly.
-Next also alternatives with a shallow concrete cut-off wall will be calculated for a shallow and a deep wall. In that case we have R=1000 d, but note that the input in the attribute than becomes k*R=10000 due to the error in Tim.
+Next also alternatives with a shallow concrete cut-off wall will be calculated for a shallow and a deep wall. In that case we have R=1000 d, but note that the input in the attribute than becomes k*R=10000 due to the error in Tim.
-(@) After changing the value in attribute tables of wall elements, we can compute again, switching off and on the element for the deep wall section (tab Geopackage on the QGIS-Tim panel).
+(@) After changing the value in attribute tables of wall elements, we can compute again, switching off and on the element for the deep wall section (tab Geopackage on the QGIS-Tim panel).
Again extra calculation is needed to adjust well extractions for drawdown in the building pit.
Results of calculations are gathered in the following table, showing extractions and head outside the wall at South East monitoring position.
@@ -512,19 +512,19 @@ Results of calculations are gathered in the following table, showing extractions
: Effect of 4 Wall alternatives on extraction rate and drawdown South East. {#tbl-CalculatedHeads2Variants}
Installation of 15 m deep sheet pile wall or cut-off wall can be elaborated in a geotechnical design.
-Still, probably some decrease of interlock leakage is needed when sheet piles are chosen. Interlock sealing or maybe irrigation of water in a shallow drain pipe around the building pit could lead to approval by authorities.
+Still, probably some decrease of interlock leakage is needed when sheet piles are chosen. Interlock sealing or maybe irrigation of water in a shallow drain pipe around the building pit could lead to approval by authorities.
-## Effect of a not closed wall (about 20 cm)
-Authorities are afraid that leakage incidents could be harmful for surrounding monuments.
-The effect of leakage can be studied by creating a fictitious little opening in the wall. This can be done e.g. by changing the values of first and last coordinate of the leaky line doublet.
+## Effect of a not closed wall (about 20 cm)
+Authorities are afraid that leakage incidents could be harmful for surrounding monuments.
+The effect of leakage can be studied by creating a fictitious little opening in the wall. This can be done e.g. by changing the values of first and last coordinate of the leaky line doublet.
(@) Select the *Leaky Line Doublet* in the *Layers* panel
-(@) In the editing toolbar go to the toggle editing option ({width=4%}) and check the vertex option ({width=4%}).
-(@) On the drawing panel select a point in the line element of the wall and right click.
+(@) In the editing toolbar go to the toggle editing option ({width=4%}) and check the vertex option ({width=4%}).
+(@) On the drawing panel select a point in the line element of the wall and right click.
-The list with coordinates appears in the vertex editor at the lower left corner of the screen. There you can edit the coordinates of all points in the line. Another option is to zoom in and select a point in the line element and drag it to a new position.
+The list with coordinates appears in the vertex editor at the lower left corner of the screen. There you can edit the coordinates of all points in the line. Another option is to zoom in and select a point in the line element and drag it to a new position.
-When the coordinates differ a gap results, in the studied case we created a 0.2 m wide gap. It might be elaborate to perform but the result is shown in the next figure. The result is calculated by using a small grid spacing. At this created gap a large flow results in the corner of the building pit, leading to an insufficient drawdown in the building pit and a 0.1 – 0.2 m larger lowering outside the building pit at that location.
+When the coordinates differ a gap results, in the studied case we created a 0.2 m wide gap. It might be elaborate to perform but the result is shown in the next figure. The result is calculated by using a small grid spacing. At this created gap a large flow results in the corner of the building pit, leading to an insufficient drawdown in the building pit and a 0.1 – 0.2 m larger lowering outside the building pit at that location.
{width=80% #fig-figure_LeakageEffectCorner}
@@ -535,8 +535,8 @@ QGIS-Tim offers the opportunity to export the geopackage of the created model to
- Calculation of model results in other Python environments, like Anaconda or Spyder or in a notebook.
- Use in other Python oriented programs, like the Probabilistic Toolkit.
-(@) If input of all elements is ready and the model has proved to run properly, go to the QGIS-Tim panel and the tab *Geopackage*.
-(@) At the bottom press the button *Convert GeoPackage to Python script*.
+(@) If input of all elements is ready and the model has proved to run properly, go to the QGIS-Tim panel and the tab *Geopackage*.
+(@) At the bottom press the button *Convert GeoPackage to Python script*.
(@) After a short period for translation in Python the explorer panel appears where you can enter the name you want to give for the python file, e.g. “case-TheHague.py” and store it in a directory you choose to save your work. The file looks like:
{width=90% #fig-figure_PythonExport01}
@@ -545,15 +545,15 @@ As can be seen the converted Python file start with calls (e.g. import timml) to
As can be seen the Python script for TimML is written in a very dedicated and condensed manner.
-To get the Python file running in a platform like Anaconda or Spyder, extra lines should be added at wish, to get the output that the user needs.
+To get the Python file running in a platform like Anaconda or Spyder, extra lines should be added at wish, to get the output that the user needs.
-Next this file can be used for geostatistical and scenario-analysis. There are several ways of handling this kind of study, like writing an additional Python program to perform repeated calculations and statistical analysis on results. But an easy way is to use the Probabilistic Toolkit (PTK), a platform for statistical analysis to be used together with geotechnical design programs, developed by Deltares. The PTK can be used for study of model sensitivity for variation of parameter values or reliability analysis.
+Next this file can be used for geostatistical and scenario-analysis. There are several ways of handling this kind of study, like writing an additional Python program to perform repeated calculations and statistical analysis on results. But an easy way is to use the Probabilistic Toolkit (PTK), a platform for statistical analysis to be used together with geotechnical design programs, developed by Deltares. The PTK can be used for study of model sensitivity for variation of parameter values or reliability analysis.
The PTK can be downloaded free of charge at Probabilistic Toolkit - download.
-(@) Open the **Probabilistic Toolkit** from your desktop ({width=4%}) or Windows menu (Deltares folder).
+(@) Open the **Probabilistic Toolkit** from your desktop ({width=4%}) or Windows menu (Deltares folder).
-The Toolkit opens at the first of 5 tabs: *Model*.
+The Toolkit opens at the first of 5 tabs: *Model*.
(@) In the section *Model Type* check if the dropdown menu *Type* is set to *Internal*.
(@) In the section *Model Type* check if the dropdown menu *Language* is set to *Python* and a the field *Version* appears.
@@ -561,7 +561,7 @@ The Toolkit opens at the first of 5 tabs: *Model*.
(@) Than you copy the content of the Python file we converted from QGIS-Tim into the window "Model code".
We handle the process in this way because we want to change some lines in the source code to get the program running in PTK.
-
+
Next the specific parameters must be selected that are expected to be probably most relevant to variations in results. In the source code used in the PTK, those parameters will not have input on a value but need to be mentioned with a name that the PTK can use for input selection in the calculations.
The parameters that seem to be important are:
@@ -570,20 +570,20 @@ The parameters that seem to be important are:
- The permeability of the sand layer k01.
- The resistance c2 of the second clay layer at -14 m NAP.
-{width=80% #fig-figure-PTK-VariableNameValue}
+{width=80% #fig-figure-PTK-VariableNameValue}
In PTK we have to set these parameters:
(@) In the PTK panel *Input* click 4 times on the *Add* button ({width=3%}), to add 4 variables and name them Rshp, czba, khol and cbasis (see @fig-figure-PTK-VariableNameValue)
-(@) In the tab *Variables* give these 4 variables the values 10000, 40, 10 and 100.
+(@) In the tab *Variables* give these 4 variables the values 10000, 40, 10 and 100.
(@) Return to the tab *Model* and assign these variable values to the parameters in the source code by copying this code block to the top of the Python source code (see @fig-figure-PTK02).
``Rwall = Rshp``
``c0 = czba``
``k0 = khol``
-``c2 = cbasis``
-(@) Finaly, replace the fixed value (**NB!** 2 times "Rwall", for each sheet pile definition) in the Python code for the parameter name (see red elements in @fig-figure-PTK02).
+``c2 = cbasis``
+(@) Finaly, replace the fixed value (**NB!** 2 times "Rwall", for each sheet pile definition) in the Python code for the parameter name (see red elements in @fig-figure-PTK02).
-{width=80% #fig-figure-PTK02}
+{width=80% #fig-figure-PTK02}
(@) At the end of the source code, eliminate the following lines from the converted file:
@@ -602,55 +602,55 @@ In PTK we have to set these parameters:
(@) In the panel *Output* of the PTK use the *Add* button ({width=3%}) the new defined variables "pb0" and "pb1".
(@) Check if the program works properly. Field *Analysis* must contain "Run model", field *Results* must contain "Single run"
-(@) Press the RUN button ({width=3%}) to calculate the model.
+(@) Press the RUN button ({width=3%}) to calculate the model.
In the tab *Run model* we find the results of our calculation.
(@) Check if it complies with your earlier calculations in QGIS-Tim.
-{width=40% #fig-figure-PTK03 fig-align="left"}
+{width=40% #fig-figure-PTK03 fig-align="left"}
-If results are as expected, we can step to a sensitivity analyses.
+If results are as expected, we can step to a sensitivity analyses.
(@) In the tab *Field*, set the Analysis option to "Sensitivity*.
(@) Go to the tab *Variables*.
For each selected parameter a distribution is defined with certain limits or characteristic values. Distribution formulas can be chosen based on knowhow of the user. Best guess of the parameter value distributions is given in the next example window.
-(@) Change the distributon types for each Variable in the column "Distrubution" and fill in the other values.
+(@) Change the distributon types for each Variable in the column "Distrubution" and fill in the other values.
-{width=100% #fig-figure-PTK04}
+{width=100% #fig-figure-PTK04}
-(@) Select the row with sheetpile resistance (Rshp) and in the panel below the distribution is shown (see @fig-figure-PTK-graph01).
+(@) Select the row with sheetpile resistance (Rshp) and in the panel below the distribution is shown (see @fig-figure-PTK-graph01).
Parameter value distributions of silt layer resistance (czba), basic peat layer resistance (cbasis) and permeability of Holocene sand (khol) are shown in graphs below.
-{width=100% #fig-figure-PTK-graph01}
+{width=100% #fig-figure-PTK-graph01}
-{width=100% #fig-figure-PTK-graph02}
+{width=100% #fig-figure-PTK-graph02}
-{width=100% #fig-figure-PTK-graph03}
+{width=100% #fig-figure-PTK-graph03}
+
+{width=100% #fig-figure-PTK-graph04}
-{width=100% #fig-figure-PTK-graph04}
-
Low and high value of resistivity are set in the tab *Calculation* to 10% and 90%.
-(@) Press RUN ({width=3%}) to recalculate the model in the Sensitivity mode.
+(@) Press RUN ({width=3%}) to recalculate the model in the Sensitivity mode.
(@) After the run is finished, a new tab *Sensitivity* is visible. Go to the tab.
Here we can find what extent parameter variations contribute to model results. In @fig-figure-PTK-Sensitivity01 it is shown what parameter variation means for drawdown in the building pit. We conclude that for a situation with a sheet pile wall in only the first sand layer the variation of the resistance of the loamy layer (czba) determines the drawdown, and with that this factor determines the amount of extraction in that situation for the largest part. Translated to practical considerations, it is worthwhile to spend extra budget on determining the homogeneity of that layer and the vertical permeability of the loamy layer in more detail.
-{width=100% #fig-figure-PTK-Sensitivity01}
+{width=100% #fig-figure-PTK-Sensitivity01}
-However, with this model we get a varying outcome for the level in the building pit and this is not realistic for a practical situation. In real projects, the dewatering contractor would like to know what the effect on dewatering demand and interaction with environment would be if the aim is to reach the designed groundwater head in the building pit while adjusting the rate of extraction in the wells.
+However, with this model we get a varying outcome for the level in the building pit and this is not realistic for a practical situation. In real projects, the dewatering contractor would like to know what the effect on dewatering demand and interaction with environment would be if the aim is to reach the designed groundwater head in the building pit while adjusting the rate of extraction in the wells.
We can do this quite simply with the aid of the Python file. To perform this analysis we follow the assumption that the relation between drawdown and well extraction is linear.
-(@) Go to the tab *Model* and in the panel *Model code* remove from the end of the script the lines for ``pb0`` and the line for ``pb1``.
-(@) On the top of the script, below the line ``k01 = khol`` add a value for the fixed total extraction flux and the fixed head within the building pit:
-``Qws = 100``
+(@) Go to the tab *Model* and in the panel *Model code* remove from the end of the script the lines for ``pb0`` and the line for ``pb1``.
+(@) On the top of the script, below the line ``k01 = khol`` add a value for the fixed total extraction flux and the fixed head within the building pit:
+``Qws = 100``
``pb0d = -3.5``
-(@) Now add the following lines to the code. **NB!** *Use the same observation variable name as is used in your script!*
+(@) Now add the following lines to the code. **NB!** *Use the same observation variable name as is used in your script!*
``pb0 = observation_observations_0[0]``
``fact = pb0d/pb0``
``Qtot = 8*Qws*fact``
@@ -662,36 +662,36 @@ We can do this quite simply with the aid of the Python file. To perform this ana
Again, we see that the parameter value for the hydraulic resistance of the loamy clay layer (czba) is most relevant to the variation of the outcome, not only for the drawdown effect in the area surrounding the building pit (pb1, see @fig-figure-PTK-Sensitivity02) but also for the total flow to the extraction wells (Qtot, see @fig-figure-PTK-Sensitivity03).
{width=100% #fig-figure-PTK-Sensitivity02}
-
+
{width=100% #fig-figure-PTK-Sensitivity03}
Based on the same statistic distributions for the parameter values we now are also able to calculate the uncertainty in output for this case, using a sheet pile wall in the first aquifer.
The distribution of the results for drawdown just outside the building pit is given in @fig-figure-PTK-Uncertainty01, which in the lower part shows the distribution for total flow rate from the dewatering in the building pit.
-
+
{width=100% #fig-figure-PTK-Uncertainty01}
{width=100% #fig-figure-PTK-Uncertainty02}
The probabilistic model gives a 50% probability on a flow of 297 m3/d, even with a 10% chance of an overrun to 659 m3/d. The groundwater drawdown at the standpipe outside the building pit amounts -0.44 m with 50% probability but with 10% value of -1.10 m.
-Moreover, it is possible to perform a calculation of reliability with the PTK.
+Moreover, it is possible to perform a calculation of reliability with the PTK.
-(@) Go to the tab *Analysis* and select the option *Reliability*.
-(@) Then we check the variable distributions again.
-(@) Next, we open the tab *Calculation*.
+(@) Go to the tab *Analysis* and select the option *Reliability*.
+(@) Then we check the variable distributions again.
+(@) Next, we open the tab *Calculation*.
-This gives the possibility to introduce a failure criterium for certain variable, like in our case pb1. We choose a less strict criterium because with a shallow sheetpile wall it is not possible to meet the requirements anyhow.
+This gives the possibility to introduce a failure criterium for certain variable, like in our case pb1. We choose a less strict criterium because with a shallow sheetpile wall it is not possible to meet the requirements anyhow.
-(@) Choose a comparison with undershooting a critical value of REF -0.5 m for the groundwater head in pb1 (see @fig-figure-PTK-Calculation01).
+(@) Choose a comparison with undershooting a critical value of REF -0.5 m for the groundwater head in pb1 (see @fig-figure-PTK-Calculation01).
At the right hand side in the settings window a calculation method needs to be selected. Many statistical methods are available, like Monte Carlo method etc. They are described in the help of the PTK in the main menu. We chose FORM (First Order Reliability Method) as a fast method because it is an automated intelligent iteration process to find a design point for reliability by using first calculated realisations (p.97 in the Help manual, see PTK menu).
-
+
{width=90% #fig-figure-PTK-Calculation01}
(@) Recalculate the model and go to the tab *Reliability*.
From the result (see @fig-figure-PTK-Reliability01) we conclude that the FORM method is able to find a solution in just a few iterations.
-
+
{width=100% #fig-figure-PTK-Reliability01}
It follows that the probability of failure for even a not very discriminating criterium as a drawdown outside the building pit of REF -0.5 m is 42%. The calculation method also returns the contributions of all relevant parameters to the outcome. It is obvious that the loamy layer at shallow depth below the building pit does not deliver enough hydraulic resistance to isolate the building activities from surrounding monuments.
@@ -707,16 +707,16 @@ At this moment, it is not possible to create a transient (time dependent) model
Nevertheless in this section we want to guide you on setting up a basic transient model.
(@) Return to QGIS and in the QGIS-Tim panel go to the tab *Compute*.
-(@) In the section *Output* turn on the option "Transient".
+(@) In the section *Output* turn on the option "Transient".
(@) In the panel *Layers* select *timml Aquifer:Aquifer*.
As can be seen in @fig-Panel-TimmlAquiferAttrTable_v02, the matrix of properties is extended in the transient mode with cells for transient parameters, like specific storage coefficients (aquitard_s, and aquifer_s) and porosities (aquitard_npor and aquifer_npor). Take care that specific storage coefficient is a storage per meter layer thickness. In older MWell course material a speadsheet is presented for estimation of coefficient values but there we calculated storage coefficients per whole layer. Those values should be divided by layer thickness.
{width=100% #fig-Panel-TimmlAquiferAttrTable_v02}
-(@) Fill in the Table with the corresponding values.
+(@) Fill in the Table with the corresponding values.
(@) In the TTim section in the *Layers panel* go to the list of elements for transient calculations.
-(@) Select *ttim Temporal Settings:Aquifer* and open its *Attribute Table* (F6). Here the length of the time series is specified.
+(@) Select *ttim Temporal Settings:Aquifer* and open its *Attribute Table* (F6). Here the length of the time series is specified.
(@) Set the starting time to 0.
But calculations take place for the range of time steps between "time_min" and "time_max". The time unit is according to central settings in the project, mostly time in days. "Time_min" should not be zero but given a slight offset. Here we also find a "Stehfest_M" number, related to the number of terms in the calculation method performing transformation of formulas in solving differential equations. A reference date can be set to relate the calculation to real time data.
@@ -727,7 +727,7 @@ But calculations take place for the range of time steps between "time_min" and "
**NB!**: In QGIS it is an obligation to use a "reference_date" to get the model running and present outcomes!
-(@) Select *ttim Computation Times:Domain* and set the time values for time steps where we want to get calculations of spatially distributed heads in mesh or raster points and contours.
+(@) Select *ttim Computation Times:Domain* and set the time values for time steps where we want to get calculations of spatially distributed heads in mesh or raster points and contours.
{width=30% #fig-Panel-TTim_DomainAttrTable_v01}
@@ -741,6 +741,6 @@ Don’t choose exactly the same time as when changing the flow. So if 60 is an e
In the *Attribute Table* of layer *ttim Well:pumping_wells* we can set the discharge amount. The value might change during time when excavation is considered or construction is performed in phases. For each change we have to add a line of input. In case there are a lot of switch on/off moments, this is a lot of work.
We just leave this option without any change because there is an other option in case each well has only a singel start and end time.
-(@) Therefor go to the layer *timml Well:pumping_wells* layer. This layer contains simple elements to make it behave as a simple Timm element.
+(@) Therefor go to the layer *timml Well:pumping_wells* layer. This layer contains simple elements to make it behave as a simple Timm element.
(@) In the *Attribute Table* make start_time=5, end_time=30 and discharge_transient=15.
(@) But fill in a value of 0 m3/d in column "Discharge" because this variable is part of TimML well element. We use the "discharge_transient" variable to define flow in TTim. Remember: Solution of TTim is superposed on solution of TimML.
\ No newline at end of file