Flood Modeller Tutorial Module02 Provisional: Difference between revisions
Content deleted Content added
(26 intermediate revisions by 3 users not shown) | |||
Line 1:
= Introduction =
In this module,
This will include:
* Modification of the floodplain topography through the creation of a 3D TIN surface.
* Revision of the land use.
* Addition of pipes and pits in ESTRY to represent the underground network.
* Linking of the pipe network in ESTRY with the Flood Modeller network.
* Introduction of an inflow into the pipe system.
* Addition of a river reach represented in ESTRY downstream of the Flood Modeller network.
=GIS and Model Inputs=
The steps required to modify each of the GIS inputs are demonstrated in QGIS using SHP and GPKG formats. Instructions for completing the module in ArcGIS or MapInfo are available on the archive page for [[Flood_Modeller_Tutorial_Module02_Archive |Tutorial Module 02]].
===Define Elevations (Building a TIN)===
The GIS layers necessary to modify the ground elevations to represent the proposed development are provided. This part of the tutorial demonstrates how a TIN is created from these GIS layers. The GIS defining the road crest level is also updated. Instructions are available below for the preferred GIS format:
* [[FM Tutorial M02_QGIS_SHP_Define_Elevations | QGIS - SHP]]
* [[FM Tutorial M02_QGIS_GPKG_Define_Elevations | QGIS - GPKG]]
===Define Surface Roughness===
The GIS layers necessary to modify the land use areas affected by the proposed development are provided. This part of the tutorial demonstrates how to populate the layer attributes to assign Manning’s n roughness values to each land use. Instructions are available below for the preferred GIS format:
* [[FM Tutorial M02_QGIS_SHP_Define_Roughness | QGIS - SHP]]
* [[FM Tutorial M02_QGIS_GPKG_Define_Roughness | QGIS - GPKG]]
===Define Pipe Network===
This part of the module creates the GIS layers that define the sub-surface pipe network. The inlets and pits of the pipe network are linked to the 2D domain. The pit inlet database is also created, linking the GIS layers to depth-discharge curves. Instructions are available below for the preferred GIS format:
* [[FM Tutorial M02_QGIS_SHP_Pipe_Network | QGIS - SHP]]
* [[FM Tutorial M02_QGIS_GPKG_Pipe_Network | QGIS - GPKG]]
===Define Boundary Conditions===
This part of the module demonstrates how an inflow can be applied directly to the pits of the pipe network. A GIS layer of the inflow boundary is provided. The existing Boundary Conditions Database is also modified to include these new inflows. Instructions are available below for the preferred GIS format:
* [[FM Tutorial M02_QGIS_SHP_Boundary_Conditions | QGIS - SHP]]
* [[FM Tutorial M02_QGIS_GPKG_Boundary_Conditions | QGIS - GPKG]]
==Flood Modeller 1D/ESTRY 1D Link==
This part of the module demonstrates how TUFLOW 1D (ESTRY) domains can be dynamically linked with Flood Modeller using "X1DH" and "X1DQ" links. Instructions are available below for the preferred GIS format:
* [[FM Tutorial M02_QGIS_SHP_Flood Modeller1D/ESTRY 1D Link | QGIS - SHP]]
* [[FM Tutorial M02_QGIS_GPKG_Flood Modeller1D/ESTRY 1D Link | QGIS - GPKG]]
The Flood Modeller 1D/ESTRY 1D link can be employed for a number of reasons, including:
* Inclusion of the pipe network and manhole modelling capabilities of ESTRY within a Flood Modeller – TUFLOW linked model (see <u>[https://docs.tuflow.com/classic-hpc/manual/2025.1/DomainLinking-1.html#OneD2DLinkingFM-2 Section 10.5]</u> of the TUFLOW Manual for more details).
* Extension of a Flood Modeller network within ESTRY to overcome Flood Modeller node licence limits.
* Representation of a steeper tributary in ESTRY which can then be connected to the main river represented in Flood Modeller.
ESTRY nodes are automatically generated at the upstream and downstream extremities of an ESTRY link, so manual generation of a node is not mandatory. If no node is manually added, the Flood Modeller–ESTRY link is assumed to be an “X1DH” link. If an ESTRY node is manually generated, the ESTRY node can have a 1d_nwk layer Conn_1D_2D attribute of either “X1DH” or “X1DQ”.
An "
=Modify Simulation Control Files=
With the input GIS layers
== TUFLOW Geometry Control File (TGC) ==
At this stage, the following changes will be made to the geometry:
* A 3D TIN is created to represent changes to the ground elevations.
<ol>
<li>Open <b>FMT_M01_001.
<li>Open <b>FMT_M02_001.tgc</b>.</li>
<li>Add the commands to modify the topography to represent the proposed development. These commands should be placed after the <tt>Read GIS Z Shape</tt> line:
<br><u>'''QGIS - SHP'''</u>
<br><tt><font color="blue">Create TIN Zpts WRITE TIN</font> <font color="red">==</font> gis\2d_ztin_fmt_m02_development_001_R.shp | gis\2d_ztin_fmt_m02_development_001_P.shp</tt>
<br><u>'''QGIS - GPKG'''</u>
<br><tt><font color="blue">Create TIN Zpts WRITE TIN</font> <font color="red">==</font> 2d_ztin_fmt_m02_development_001_R | 2d_ztin_fmt_m02_development_001_P</tt>
<br>The <tt><font color="blue">Create TIN Zpts WRITE TIN</font></tt> command creates and writes a .2dm mesh file to the same location as the GIS layer (in this case the TUFLOW\model\gis folder). The .2dm TIN can be viewed, checked, and modified in QGIS. It can then be read into the model directly using the <tt>Read TIN Zpts</tt> command for subsequent model simulations.
The 2d_mat layers created in this module build upon the existing commands that modify roughness. The new layers overwrite the existing layers at the location of the proposed development. This process of layering provides a powerful tool in TUFLOW that minimises data duplication and offers a means of quality control. The commands reading in the new 2d_mat layers must be placed after the existing commands.
<br><u>'''QGIS - SHP'''</u>
<br><tt><font color="blue">Read GIS Mat</font> <font color="red">==</font> gis\2d_mat_fmt_m02_dev_001_R.shp</tt>
<br><tt><font color="blue">Read GIS Mat</font> <font color="red">==</font> gis\2d_mat_fmt_m02_dev_buildings_001_R.shp</tt>
<br><u>'''QGIS - GPKG'''</u>
<br><tt><font color="blue">Read GIS Mat</font> <font color="red">==</font> 2d_mat_fmt_m02_dev_001_R</tt>
<br><tt><font color="blue">Read GIS Mat</font> <font color="red">==</font> 2d_mat_fmt_m02_dev_buildings_001_R</tt>
</li>
<br>
<li>Save the file. The TGC is now ready to be used.</li>
</ol>
== ESTRY Control File (ECF) ==
At this stage, the following changes are made to the ECF file:
* A 1d_nwk layer is created to represent the culverts of the proposed pipe network.
* A 1d_nwk layer is created to represent the pits of the proposed pipe network.
* A pit inlet database is created to link depth-discharge curves to the pit inlet type.
<ol>
<li>Open <b>FMT_M01_001.ecf</b> in your text editor. Save the file as <b>FMT_M02_001.ecf</b>.</li>
<li>Open <b>FMT_M02_001.ecf</b><br></li>
<li>Add the following commands at the bottom of the file as follows:
<div>
<u>'''QGIS - SHP'''</u><br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> gis\1d_nwk_FMT_M02_Pipes_001_L.shp<br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> gis\1d_nwk_FMT_M02_Pits_001_P.shp<br>
<font color="blue"><tt>Pit Inlet Database</tt></font> <font color="red"><tt>==</tt></font> ..\pit_dbase\pit_inlet_dbase.csv<br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> gis\1d_nwk_FMT_M02_Channels_001_L.shp<br>
<font color="blue"><tt>Read GIS Table Links</tt></font> <font color="red"><tt>==</tt></font> xs\1d_xs_FMT_M02_Creek_001_L.shp<br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> gis\1d_nwke_X1DQ_P.shp<br>
<font color="blue"><tt>Read GIS BC</tt></font> <font color="red"><tt>==</tt></font> gis\1d_BC_FMT_M02_001_P.shp<br>
<u>'''QGIS - GPKG'''</u><br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> 1d_nwk_FMT_M02_Pipes_001_L<br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> 1d_nwk_FMT_M02_Pits_001_P<br>
<font color="blue"><tt>Pit Inlet Database</tt></font> <font color="red"><tt>==</tt></font> ..\pit_dbase\pit_inlet_dbase.csv<br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> 1d_nwk_FMT_M02_Channels_001_L<br>
<font color="blue"><tt>Read GIS Table Links</tt></font> <font color="red"><tt>==</tt></font> 1d_xs_FMT_M02_Creek_001_L<br>
<font color="blue"><tt>Read GIS Network</tt></font> <font color="red"><tt>==</tt></font> 1d_nwke_X1DQ_P<br>
<font color="blue"><tt>Read GIS BC</tt></font> <font color="red"><tt>==</tt></font> 1d_BC_FMT_M02_001_P
</div>
</li>
<li value="4">Save the file. The ECF file is now ready to be used.</li>
</ol>
== TUFLOW Boundary Control File (TBC) ==
<ol>
<li>Open <b>FMT_M01_001.tbc</b> in
<li>Insert the following commands after the existing <font color="blue"><tt>Read GIS SA</tt></font><font color="red"><tt> ==</tt></font> mi\2d_sa_M01_002_R.
<u>'''QGIS - SHP'''</u><br>
<font color="blue"><tt>Read GIS SA PITS</tt></font><font color="red"><tt> ==</tt></font> gis\2d_sa_FMT_M02_001_R.shp<br>
<u>'''QGIS - GPKG'''</u><br>
<font color="blue"><tt>Read GIS SA PITS</tt></font><font color="red"><tt> ==</tt></font> 2d_sa_FMT_M02_001_R<br>
<li>
</ol>
==TUFLOW Control File (TCF)==
We will need to create a new tcf file that references the new tgc, ecf and tbc files.
Line 137 ⟶ 149:
:<font color="blue"><tt>BC Control File</tt></font><font color="red"><tt> ==</tt></font> ..\model\FMT_M02_001.tbc
<li>We have also created a new bc_dbase in this module which will need to be referenced. Update the command as follows:
:<font color="blue"><tt>BC Database</tt></font><font color="red"><tt> ==</tt></font> ..\bc_dbase\
<li>We are using an updated GIS layer for the Flood Modeller Nodes. Update the command as follows:
<br><u>'''QGIS - SHP'''</u><br><font color="blue"><tt>Read GIS X1D Nodes</tt></font> <font color="red"><tt>== </tt></font>..\model\gis\1d_x1d_FMT_M01_nodes_001_P.shp <font color="green"><tt>! GIS layer referencing node IDs from Flood Modeller </tt></font><br><u>'''QGIS - GPKG'''</u><br><font color="blue"><tt>Read GIS X1D Nodes</tt></font> <font color="red"><tt>== </tt></font> 1d_x1d_FMT_M01_nodes_001_P <font color="green"><tt>! GIS layer referencing node IDs from Flood Modeller</tt></font>
<li>Lastly, update the following command to specify a new output folder for the results of this module:
:<font color="blue"><tt>Output Folder</tt></font><font color="red"><tt> ==</tt></font> ..\results\FMT_M02\2d
Line 145 ⟶ 159:
<ol>
<li>Open the Flood Modeller Pro model as per [[Flood_Modeller_Tutorial_Module01_Provisional | Flood Modeller Tutorial Module 1]]
<li>Create a copy of '''FMT_M01_001.ief''' and save as '''FMT_M02_001.ief'''.
<li>In the 'Links' tab with the 2D scheme set as TUFLOW, change the full path of the 2D control file to the <b>FMT_M02_DEV_001.tcf</b> from the <b>\FMT_Tutorial\FMT_M02\TUFLOW\runs</b> folder
<li>Save the Scenario Data.</li></ol>
Use
=Review Check Files=
Once the model has compiled and the simulation started, we can review the check files to ensure the changes have been correctly applied. The following section of this module outlines how the generated check files can be used to review each of the key changes we have made to the model. Note that there is often more than one check file that can be used to review each component of the model. The below steps outline just how some of these check files can be used.
From the <b>TUFLOW\check\2d\</b> folder open within QGIS:
*<b>FMT_M02_001_DEM_Z.TIF</b>
Line 160 ⟶ 174:
The FMT_M02_001_DEM_Z.TIF is a grid of the final zpts used by TUFLOW after processing of each of the layers within the TGC. It can be read by most mainstream GIS software and also visualised in SMS.
From the <b>TUFLOW\check\2d\</b> folder open within QGIS:
*FMT_M02_001_grd_check
Line 170 ⟶ 183:
*By right clicking on the grd_check layer within the table of contents and click Properties. Once the properties dialogue window is open select the Symbology tab. Choose the Unique values option under Categories in the left-hand list and map the Value Field of interest.
From the <b>TUFLOW\check\1d\</b> folder open within QGIS:
*FMT_M02_001_nwk_C_check
Line 178 ⟶ 190:
These two check files provide information on all 1D elements within the model. The _nwk_C_check layer provides information on all 1D channels (including structures) in the model, whilst the _nwk_N_check layer provides information on the nodes. Notice how the _nwk_C check layer shows a series of dashed lines at the locations where we have specified pits in the model. These lines represent 'pit channels', zero length channels that connect the 1D pipe network to the 2D floodplain. Viewing in conjunction with the _nwk_N check layer shows two nodes are created at the upstream (ground level) and downstream (invert of the pipe network) ends of the pit channel.
From the <b>TUFLOW\check\2d\</b> folder open within QGIS:
*FMT_M02_001_1d_to_2d_check
The _1d_to_2d_check layer highlights the location of all 1D/2D boundary links within the model. Where we have digitised the pipe network, this check file shows the SX boundaries that have been created at the location of the pits. This check file also shows the ZC elevation of the SX boundary which can be compared to surrounding ZC elevations by viewing alongside the grd_check layer. The check layer also shows the HX boundaries digitised along the river banks in Module 01.
=== Review Flood Modeller/ESTRY 1D/1D Links ===
The successful connection of the Flood Modeller and ESTRY networks can be checked by the presence of CHECK 1393 messages highlighting the presence of "X1DH" and "X1DQ" links to external nodes within the TUFLOW Log File. Alternatively, the success of the connectivity of the "X1DQ" and "X1DH" links can be assessed through reference to the _messages.SHP/GPKG GIS layer which contains CHECK 1393 at each ESTRY node linked to a Flood Modeller node.
=Review the Results=
Open the 2D results in
[[File:M07 2d Results.PNG|800px]]<br>
<br>
The 1D ESTRY results can be viewed by opening the time series results using the approach described in the following page: [[TUFLOW Viewer - Load Results - Time Series]]. You can also open results from Flood Modeller within the TUFLOW Viewer using the instructions here: [[TUFLOW Viewer - Load Results - Time Series FM]]
It can be seen that there is negative flow in Pipe16 and Pipe18 when the water levels in the Flood Modeller node and the connected ESTRY nodes are such that backwater effects are present.
[[File:X1dh.png|alt=X1DH Flows from Estry to Flood Modeller at Pipe 16 and 18|none|thumb|500x500px|X1DH Flows from Estry to Flood Modeller at Pipe 16 and 18]]
A sense check of the flows being transferred across the "X1DQ" link can also be performed by plotting and comparing the flow time series for the first ESTRY node (ds3) and the last Flood Modeller node (ds2).
[[File:X1dq.png|alt=X1DQ Flow from Flood Modeller to Estry at DS3|none|thumb|500x500px|X1DQ Flow from Flood Modeller to Estry at DS3]]
We have now completed the tutorial and you should now be familiar with the approaches of linking Estry to Flood Modeller and vice versa.
{{Tips Navigation
|uplink=[[Flood_Modeller_Tutorial_Model_Provisional| Back to Tutorial Introduction Main Page]]
|