Master thesis. Major Disaster Management: Research on large flood management of Tra Bong River basin, Quang Ngai province

This study aimed at proposing solutions for managing large flood in Tra Bong riverbasin, Quang Ngai province, The research includes 2 main parts: create a floodinundation map based on hi

Firstly, I would like to express my sincere gratitude to my advisors Assoc.Prof.

Hoang Thanh Tung and Dr Pham Thanh Hai for the continuous support of my M.Sc.

study and related research, for their patience, motivation, and immense knowledge Their guidance helped me in all the time of research and writing of this thesis I could not have imagined having a better advisor and mentor for my M.Sc study.

Besides, I am especially grateful to lecturers in the Department of Hydrology and Water resources, Thuy Loi University and foreigner lectures from NICHE project who supported me for all the lectures and useful advices throughout my course.

My sincere thanks also goes to my colleagues in National

Hydro-Meteorological Service and Hydro-Hydro-Meteorological Department of Middle Centre

region who supported me for data collection and analysis Without they precious support

it would not be possible to conduct this research.

Last but not the least, I would like to thank my family: my parents, my wife and

children for supporting me spiritually throughout the course and my life in general.

Hanoi, Nov 08" 2016

Nguyen Phu Luan

This study aimed at proposing solutions for managing large flood in Tra Bong riverbasin, Quang Ngai province, The research includes 2 main parts: create a floodinundation map based on historical flood marks collected and flood forecastingexperiment for the downstream of river by simulating flood flow with the combination

of hydrological and hydrodynamic models At the first part, 42 flood marks of the

2009 historical flood had been collected with their exactly coordinate and maximumwater level The distribution and elevation of flood marks have been used to calculatewith topography data (DEM) using GIS tools and the result was the 2009 floodinundation map At the second part, a rainfall ~ runoff hydrological model (NAM) hasbeen firstly used to simulate flow from upstream of basin to the section at thebeginning of the main river The input data were collected from 3 rain gauge stationsfor simulation After that, the computed water discharge got from rainfall-runoff

‘model have been used as the upstream boundary for simulating flood flow in the mainriver The model has been used was I-dimensional hydrodynamic model Aftersimulating the flood flow in the river, the study tried flood forecasting for thedownstream part Based on the flood inundation map and the result of floodforecasting experiment, the research has proposed solutions to manage large floods onTra Bong river basin, These solutions can be used to develop large flood managementplans for local authorities, in order to enhance efficiency in water resources using andreduce losses caused by large flood,

Digital Elevation ModelDisaster Management CycleMinistry of Natural Resources and EnvironmentMinistry of Construction

National Hydro-Meteorological Sevice

‘Nedbor-Afstromnings-Model (Danish, meaning rainfall-runoltmodel)

Flood ForeeastingHydrodynamic1-dimensiona

BLE OF CONTENTS

INTRODUCTION

1 General Introduction

2 Description of the Study Area

3 Problems and Need of Study

4 Objectives of Study

5 Scope of Study

CHAPTER I: LITERATURE REVIEW

1.1 Related researches about study site

1.2 Flood inundation mapping,

1.3 Geographical Information Systems in Hydrology and Water Resources

1.4 Flood Forecasting

(CHAPTER I: APPROACH AND METHODOLODY:

2.1 Approach of study

2.2 Flood inundation map

2.3 Using GIS tools to develop flood inundation map

24 Mike 11 general description

2.5 Theoretical Foundation of rainfall - runoff hydrological model (NAM),

25.1 The ic parameters of NAM model

2.5.2 Basic modelling components

2.5.3 Initial conditions ofthe model

2.5.4 Model calibration

2.6 Flood Forecasting (Mike 11 FF): Updating procedure

2.6.1 Two unique features MIKE 11 FPS updating procedure

2.6.2 The calibration updating parameters

27 Flood Forecasting Error

(CHAPTER III: RESULTS AND DISCUSSIONS

3.1 Analysing flood features of Tra Bong river basin

3.1.1 Rain features

3.1.2 Flood features

303242

4650sỉsỉ

5537

”s5961

3.2 Develop flood inundation mapping for Tra Bong river basin.

NAM3.3 Flow forming simulation using Mi

3.31 Caleulation layout

3.3.2 Data analysis

3.3.3 NAM model calibration and verifieation

3.4 Flood flow in the downstream using Mike 11 hydraulic model

3.4.1 Boundary condition

3.42 Calibration and verification of flood flow simulation model

3.5 Flood forecasting experiment for Tra Bong river system

3.6 Propose solutions of large flood managing in Tra Bong river basin

n

$284

95

9

tô

101102

103

LIST OF FIGURES

Figure 1: Administration map of Quang Ngai Province 9Figure 2.1: Structure of literature review 13Figure 2.2: Conceptual framework for flood hazard and risk calculations 19Figure 2.3: Different flood map types 20Figure 2.4: Process for developing a flood forecasting model 25Figure 3.1: Conceptual Framework 2Figure 3.2: Overview of study 28Figure 3.3: Digital Elevation Model with square grid 31Figure 3.4: Error based on the topography data in flood inundation mapping 31Figure 3.5: The structure of NAM model 38Figure 3.6: Channel section with computational grid “Figure 37: The shape of the computational grid around a node which has three

branches 45

Figure 3.8: The shape ofthe grid points and the nodes in the complete modkl 45Figure 3.9: Branch matrix before reducing 4Figure 3.10: Branch matrix after reducing 47Figure 3.11: Three-branch node with limit for continuity equation 48

Figure 3.12: River branch with discharge boundary 49

Figure 3.13: Illustration of amplitude and phase error 33Figure 3.14: The updating results of simulations 35Figure 3.15: Example of measured and simulated discharge at an update location 57Figure 3.16: Example of updating parameters 31Figure 4.1: Chart of the possibility of flooding which reach 2nd alarm level or higher

in flood season - Tra Bong river basin, 6Figure 4.2: 6 hour rainfall chart, from 19:00 27 Sep 2009 to 01:00 30 Sep 2009 - Tra

Bong station 65Figure 4.3: 6 hour rainfall chart, from 19:00 27 Sep 2009 to 01:00 30 Sep 2009 ~

Chau 0 station 66Figure 4.4: Hourly water level process of the flood from 28 Sep to 01 Oct, 2009 ~

Chau 0 station, 67

Figure 4.5: 2009 Flood marks map of Tra Bong river downstream 70Figure 4.6: 2009 Flood inundation map of Tra Bong river basin T73Figure 4.7: Map of hydro-meteorological stations network 4Figure 4.8: Calculation layout T5

Figure 4.9: Weight factor distribution layout of rain gauge stations in Quang Ngai

province T6

Figure 4.10: Basin parameters declaring dialog T¡Figure 4.11: Parameters calibration for rainfall-runoff model (NAM), 7Figure 4.12: Calculated flow process at Binh Minh station 8Figure 4.13: Rainfall — runoff model (NAM) calibration, compare observed and

ated flood discharge at Bình Minh ~ Tra Bong river, from 16 to 20

2008 $0Figure 4.14: Rainfall - runoff model (NAM) calibration, compare observed and

simulated flood discharge at Binh Minh - Tra Bong river, from 28 Sep to

5 Oct, 2009 $0Figure 4.15: Rainfall = runoff model (NAM) calibration, compare observed and

simulated flood discharge at Binh Minh ~ Tra Bong river, from 13 to 19Nov, 2010 81Figure 4.16: Rainfall — runoff model (NAM) verification, compare observed and

simulated flood discharge at Binh Minh ~ Tra Bong river, from 14 to 20Oct, 2011 81Figure 4.17: Rainfall ~ runoff model (NAM) verification, compare observed and

simulated flood discharge at Binh Minh ~ Tra Bong river, from 5 to 9 Nov,

201 82

Figure 4.18: Rainfall — runoff model (NAM) veri

simulated flood discharge at Binh Minh —

jon, compare observed and

a Bong river, from 25 to 29Nov, 2011 9Figure 4.19: Hydraulic routing layout of Tra Bong river downstream 83

Figure 420: Initial conditions ealibration dialog $6

Figure 421: Bed resistance calibration dialog $6Figure 4.22: Bed resistance calibration for cross-sections 87Figure 4.23: MIKE 11 (HD) calibration, compare observed and simulated water level

at Chau O station — Tra Bong river, from 16 to 20 Oct, 2008 88Figure 4.24: MIKE 11 (HD) calibration, compare observed and simulated water level

‘at Chau O station ~ Tra Bong river, from 28 Sep to 5 Oct, 2009 88

Figure 4.25: MIKE 11 (HD) verification, compare observed and simulated water level

at Chau O station — Tra Bong river, from 13 to 19 Nov, 2010 90Figure 4.26: MIKE 11 (HD) verification, compare observed and simulated water level

at Chau O station — Tra Bong river, from 14 to 20 Oct, 2011 90

Figure 4.27: Flood flow forecasting results from rainfall data at Binh Minh — The rains

from 5 to 9 Nov, 2011 92

Figure 4.28: Observed and 5 hour predicted flood water level process in Chau O

station - from 5 to 9 Nov, 2011 93Figure 4.29: Flood flow forecasting results from rainfall data at Binh Minh — The rains

from 25 to 29 Nov, 2011 93Figure 4.30: Observed and 5 hour predicted flood water level process in Chau O

station - from 25 to 29 Nov, 2011 9Figure 5.1: Synthesis of large flood management solutions in Tra Bong river basin 97Figure 5.2: Disaster Management Cycle 99

corresponding scales 32

‘Table 3.2: Quality of forecasting classification 59

‘Table 4.1: Rainfall causes flood rising in Tra Bong river basin 61Table 4.2: Rainfall in the history flood on 28-30 Sep 2009 61

‘Table 4.3: Maximum daily rainfall in Tra Bong station, from 2006 to 2011 61

‘Table 4.4: Annual peak flood of Chau O station (2006-2011) 2Table 4.5: Typical flood intensity and amplitude in Chau O station (with the peak

flood higher than first alarm level) @

‘Table 4.6: The typical floods in 2006-2011 period of Chau O station 6Table 4.7: 6 hour rainfall data, from 19:00 27 Sep 2009 to 01:00 30 Sep 2009 ~ Tra

n 65acteristics of the flood from 28 to 30 Sep 2009 66Table 4.9: The data of 2009 flood marks in Tra Bong river basin 67

‘Table 4.10: The rain gauges used for hydrological calculating 75Table 4.11: Results of parameters calibration for rainfall — runoff model 79

‘Table 4.12: Results of NAM model calibration and verification at Bình Minh, 79Table 4.13: Position of nod xin hydraulic calculation layout of Tra Bong river 83

‘Table 4.14: Analysis of model calibration efficieney and error 87Table 4.15: Analysis of model verification efficiency and error 89Table 4.16: Results synthesis of flood forecasting of Tra Bong river in Chau O

station 95

1 General Introduetion

Large flood is a kind of disaster occur regularly and seriously annually around theworld Flood has some benefits such as bring fertile soil to replace nutrient-poor soils,but it also impacts and causes enormous damage, constantly threatening people li ngand the economic and social development, People use many methods to prevent andreduce the impact of flood including management, structures and non-structures

‘methods Therein, large flood management is always an important objective, which

requires detailed and specific research for each region,

The research “Research on Large flood management of Tra Bong river basin, QuangNgai province” aims at approaching a new point of view in flood management, inorder to achieve high efficiency in reducing the harmful effects of flood to theeconomic development and environmental protection of study area

2 Description of the Study Area

Quang Ngai province is located at the latitudes 1432: 15°25' North, longitudes108°06' - 109°04" East, lean on Truong Son mountain range, overlooking the Eastem

ca Quang Ngai abuts Quang

South, Kon Tum province in the West, and Eastem sea in the East Located in the

fam province in the North, Binh Dinh province in the

‘middle of the country, it is 883km from Hanoi capital to Quang Ngai and 838km from

Ho Chi Mink City

Due to the steep topography and poor vegetation, so the transfer speed of flood is very

high, flood’s damage is huge Every year floods have caused extensive damage andlosses of life and property 600,000 people affected by flooding, especially somecommunes which suffered flood inundation elevation more than 3 m with thefrequency of floods are 20% and 10%

‘Tra Bong river is one of the 4 biggest rivers of Quang Ngai province including Tra

Bong, Tra Khue, Ve and Tra Cau Tra Bong starts from the western mountains of Tra

Bong District, running through Binh Son district to the sea at the mouth of Son Tra

‘The length is 45 km, which flows from southwest to northeast The majority of the

length of the river is mountainous terrain with the elevation of 200-1300m, the rest of

river flowing in alternating hills and barren plains and sandy beaches

Upstream of Tra Bong river has many tributaries including rivers and streams, as Nun

stream, Ca Du stream and Tra Boi river in Tra Thuy, Tra Giang communes Indownstream, the eastern of Binh Son district is relatively high, so Tra Bong river has

‘no longer flow rapidly as upstream, Tra Bong river has four I" level branches, At the

downstream, there are some small rivers and streams enter the main river before

flowing into the sea,

Tra Bong river basin covers most of Tra Bong district and Binh Son district The

catchment area is about 697 km”, The average elevation is 196m, average slope is

10.5%, river density is 0.43km/km”

More than half of the catchment area in the upstream is high mountain jungle, forest

revival, the rest in downstream is barren hills interspersed with agricultural land.

Sandy areas are along the estuaries and coastal account for a small part

Table 1.1: The morphological features of Tra Bong river and major tributaries

AverageName of rivers or | Lengthof | Lengthof | Area of

+ | width of

tributaries rivers (km) | basins (km) | basins (km?)

basins (km)Tra Bong 45 36 697 124

‘SaThin 2 8 30 62

Tra Boi 3 8 29

Nau River 1g | T8 8a

Ban Dien 4 Boa 33

‘Tributary Nos bì Bo ane 86

= Water level measuring station network

Currently, there are water level stations on all main rivers of the province However,

only Tra Khuc river and Ve river measure water level throughout the year, Tra Bongand Tra Cau river measure only during flood season (September to December)Specific:

‘The station on Tra Bong river is Chau O, in the town of Chau O, Binh Son district

- Rain gauge station network:

‘There are three rain gauge stations that have been collected data for use in this studyare Tay Tra, Tra Bong and Chau O station

“The meteorological and hydrological situations of study area including rain and floodfeatures are described and analysed in detail in Chapter IIL

3 Problems and Need of Study

Every year, Quang Ngai province suffered so many disasters, especially large floods

In the year 2009, an extreme flood caused severe flooding on a large scale of

downstream areas of Quang Ngai Some river nd tributary basins have peak floodelevation even higher than history flood in 1999 Through the statistic data, we see the

‘annual flood damage caused hundreds of billions: VND lost, hundreds of people werekilled and injured in Quang Ngai province: especially the 2009 flood caused over fourthousand billions VND lost and up to 51 people dead, 506 people injured The

government at central and local levels took immediate action during and after the

floods to assist the local communities to cope with the crisis and to help restore foodsupplies, essential services, and essential infrastructure Further more, An IntegratedNatural Disaster Mitigation Policy for Central Viemam has been prepared whichinvolves both structural and non-structural measures However, the non-structuralmeasures are more concentrated and preferable because of the current condition offund limitation in Vietnam, In these non-structure measures, improvement of floodmanagement measures is those of highest priority Therefore, I ehose this study inorder not only to learn myself knowledge of flood management, but also paving the

‘way for the further researches, which are useful for local province

4 Objectives of Study

‘The overall objective of this study is to research large flood management measures for

‘Quang Ngai province of Vietnam,

To obtain this objective, some sub-objectives that need to be achieved are as follows

‘Analyse flood features of Tra Bong river basin, flood season and the typical flood,

* Develop flood inundation mapping for Tra Bong river basin based on flood mark'sdistribution and characteristics of the major flood in 2009

« Using rainfall - runoff hydrological model and I-dimensional hydrodynamicmodel to simulate flood of Tra Bong river

+ Flood forecasting experiment for Tra Bong river system,

‘© Propose solutions for large flood management

socio-+ Based on data collected, analysis flood mode of Tra Bong River basin include:flood season, flood characteristics (peak flood, flood intensity, floods hydrographand typical flood process)

‘+ Topographic data of the basin collected as a digital elevation model (DEM) Based

on the data of flood marks coordinates and elevation of 2009 major flood, usingGIS tools to take position of flood marks on the map, thereby creating flood.inundation maps for the downstream area of the basin,

‘© Using rainfall — runoff hydrological model to simulate the process of flow formingfrom data of rainfall gauges in the basin, The result is the hydrograph at theposition of last section in the upstream part of Tra Bong river

+ Using hydrograph calculated from rainfall-runoff model as the upstream boundary

to simulate flood flow in the downstream, The model used is I-dimensionalhydrodynamic model

‘+ Applying the hydrological — hydrodynamic combine model to forecast flood for

‘Tra Bong river downstream, Using the measured data of several floods to verify theforecasting reliability of the model

‘© Propose solutions for large flood management of the basin,

(CHAPTER I: LITERATURE REVIEW

‘There are many researches on flood management using GIS tools in Vietnam and

around the world, This part of the study summaries the review of related literatures instructured way shown below:

— mapping

Figure 2.1: Structure of literature review

1-1 Related researches about study site

‘The study area is the Tra Bong River Basin has some characteristics that can causemajor flooding as the total rainfall in the basin is big, the basin is short and steepCharacteristics of the terrain and climate are also common features of the centralprovinces of Vietnam in general and the river systems of Quang Ngai province in

ts andparticular Thus, the study of floods in this area is relatively respected by scient

the local government, There were many studies of floods and disasters in Centralprovinces, Quang Ngai or even in Tra Bong river basin in recent years, which are very

‘200d references for the research objectives of the thesis,

In 2000-2001, the Vietnam Institute of Meteorology, Hydrology and Climate Changehhas performed a national independent reseach: "Investigations, studies and flood

‘warning support disaster prevention in the Central Basin” (Du, 2001) This tesearch

played an important role in hydro ~ meteorological research activities supportingdisaster prevention for the region, as the basis for subsequent detailed studies,

Within the scope of Quang Ngai province, there were some studies on flooding inrecent years as follow:

Disaster Mitigation Project in Quang Ngai province was supported by the AustralianGovernment and implemented from 2003 to 2007 One of the contents ofthis project

‘was constructing hydraulic model for floodplain management, Hydraulic model was a

tool to verify the management plans include floodplain basins of Tra Bong, Tra Khuc,

Ve and Tra Cau rivers Constructed hydraulic model was a detail 2D model, usingSOBEK software from Delft Hydraulics and some other supporting software products

In addition, the hydraulic model was also used to examine the impact of futuredevelopment to flooding,

In 2010-2011, Institute of Geography - Vietnam Academy of Science and Technology

hhas implemented the research: Research on flood prevention and drainage planning of

‘Tra Khue river and Ve river in Quang Ngai province (Huong, 2011) This study hasbeen taken follow the objectives:

+ Setting up the database of floods and flood evolution supporting flood preventionand drainage planning of Tra Khuc river and Ve river in Quang Ngai province

+ Determination of flood drainage way for estuary deltas (from the Thach Nham dam

to Dai estuary in Tra Khuc River; from Hanh Tin Dong to Lo estuary in Ve river)corresponding to a frequency of 1%, 5% , 10%

+ Proposed planning options of flood prevention and drainage planning of Tra Khucriver and Ve river in Quang Ngai province

In particular, research and management agency on hydrometeorology in the region isthe Hydrometeorology Department of Mid-Central region has taken some specificresearches on floods related to Tra Bong river basin, as follows:

In 2000-2001, implemented the research: Develop flood zoning maps and forecastwarning flood plans for rivers in Quang Ngai province (Lien, 2001) The study has

ed the flood regime of Tra Bong, Tra Khuc and Ve river basins, survey floodsituation in 1999, thereby mapped downstream inundation on the rivers for the 1999large flood and also flood risk maps corresponding to the design frequency This studyhas also built the first 32 flood marks in Quang Ngai province

In 2006-2007, implemented the research: Additional investigation and flood markbuilding in floodplains of main rivers downstream in Quang Ngai province (Lien,2007), Results of study were the new flood alert levels on the three main rivers of theprovince included Tra Bong, Tra Khuc and Ve rivers The flood alert levels proposed

by research had been officially used due to Decision No 632/QĐ-TTg on 10/5/2010 ofthe Prime Minister and Decision No 1142/QD-UBND on 25/8/2010 of the President

of the Committee of Quang Ngai province In addition, the study continued to buildadditional 32 flood marks, the flood risk maps base on the flood alert levels Most ofthe results of study were used for flood prevention in recent years,

Expanded and continued previous studies, (Chien & Thiem, 2012) had researchedabout flood in Tra Bong river in order to correcting flood risk map and flood warningmark following the new alarming levels In this research, the authors using GIStechnology for mapping and the research aim at correcting the warning marks, whichare often used in hydrological stations system, and they are very useful for farmerthemself see the warning level of flood, The research had a good overview about flood

in Quang Ngai province, but it did not analyze much about flood management,especially large flood However, the research helps us to have an overview on thesituation of floods in the area of research and is the basis for us to compare, evaluatethe accuracy of the data collected

In another research, (Hanh, 2014) has simulated flooding in downstream area of TraBong river The major research area is Binh Son district, where the Highway 1Arunning through, The main content ofthe research is using some 1D and 2D hydraulic

‘models to simulate flood flows and inundation calculations for the study area, whichproposed solutions in terms of both construction and non-construction to prevent floodinundation for the 1A national Highway section running through Quang NgaiResearch has simulated very specific and potentially apply to the local, However, thecontent of the study did not focus on flood management issues and the main objective

is to develop solutions to prevent flooding for the highway

Research topies related to flood management in Tra Bong river basin and Quang Ngaiprovince were mainly carried out by the local authorities as Provincial People’sCommittee, Regional Hydro-meteorological Department, Center for DisasterPrevention and Reduction, In flood management solutions had been proposed,solutions of downstream flood forecasting from measurable rainfall data anddownstream inundation mapping got high efficiency There were two main methodsused in studi in the region, including:

+ Traditional methods: survey, collect marks of the historic floods, based on the

distribution and elevation of flood marks to build inundation maps Construct thecorrelation between rainfall data in the upstream and downstream flood levels,thereby develop the flood alert levels to warn local residents, The researches usingthis method include a number of researches implemented by HydrometeorologyDepartment of Mid-Central region (Chien & Thiem, 2012; Lien, 2001, 2007),

+ Modelling methods: Apply the combination hydrological rainfall ~ runoff modelsand hydrodynamic models to simulate the process of flow forming in the basinusing rainfall data from rain gauges and flood flow progression in the riverdownstream Apply GIS tools to build flood map for each flood alert level to

provide spatial view of the flood risk to the user Typical researches using this

007; Hanh, 2014:method is (Disaster Mitigation Project in Quang Ngai province,

Huong, 2011),

Each methodology had advantages and disadvantages However, despite using anymethod, the studies of flood in this area had great impacts on socioeconomicdevelopment in localities, especially in flood warning and forecasting on river basins.

“The studies also provided an overview on the situation of flood and natural conditions,topography, climate and other features for the author to perform this study

1.2 Flood inundation mapping

In studies on flood management in Vietnam and around the world, the construction offlood inundation maps plays a very important role Flood inundation maps provide aspatial view of the possibility of flooding for each region, thereby plan for preventingand mitigating the damage caused by floods The construction of flood inundationmaps using many methods: Investigation and survey method, statistical method,system analysis method, modelling method etc or the combination of two oF moremethods, depending on the requirements and objectives of the research,

Studied on large floods and large flood management solutions, especially extremefloods, (Buchele et al., 2006) aim to enhance existing methods for hazard and risk

assessment for extreme flood events In this study, the authors analysed the process to

develop flood risk maps in three major steps:

‘+ Regional estimation of flood discharges (basin-, site-specific hydrological loads)+ Estimation of flow characteristics in potential inundation areas (local hydraulic

impacts)

+ Estimation of the resulting damages (area- or object-specific risk assessment)

Furthermore, the authors also analysed two types of requirements in developing floodrisk map The first are minimum requirements on data and methods for a standardquality of hazard and risk assessment on local scale, The second are moresophisticated approaches which require more spatial information and more complexcalculations up to fully dynamic simulations of unobserved extreme flood situations

‘The overview of suitable approaches for these steps of hazard and vulnerabilityassessments is given in Table 1

Table 2,1: Overview of methods and data for high-resolution flood-risk mapping in

Germany (Búchele et al., 2006)

Brac apa ‘Mor deed approach Qasr sta =—— ‘remain, icon znc

pec)

—— "" ` s.Ả

for gauged basis) —¬

coer dor engaged as) vay bi probabil

‘one elimate end)

Estimation of bơi y=

——— "—

— fetes of Sow Seber

— Cakealtion of wate Ieee ana Sesectow)

{undo zones depths based — waeaty bywĐymmắc saul

sa tydraae modes sả đại ten ef ceeee food seman

cet)

| Simpliied approaches (onl

large-scale 1/0000) Daygc chai

"`

fed quái a eg ATKIS vidual bulla or BmôUMG

mm (eg ALK dao)

— consideration of further -đansgeddetmlang — fader

(ean tow veo prea

flonary mesures Waraing

‘ie

Inthe study of (Duy Kieu, 2012), the flood risk maps is developed by a combination ofStatistical method and system analysis method The map was developed through thefollowing steps: 1 Zoning basin into smaller areas base on the location of monitoring

stations, 2 Develop criteria to assess the risk of large floods occurred in each region

(the criteria are mainly based on the measured meteorological and hydrological data),

3 Scoring for each region based on the above criteria 4 Using GIS to establish floodrisk maps This method has the advantages of simple, easy to calculate, easy t0 carry,

‘mainly based on the systems analysis and synthesis However, the study area of theresearch should be relatively large and the number of monitoring stations is largeenough to be able to partition flood risks in detail and with high accuracy In my

research, I referred the method of statistical data analysing on floods and hydrological

characteristics of (Duy Kieu, 2012) to analyse flood features of Tra Bong river basin

(Tủ & Tingsanchali, 2010) had made a research to create flood risk map for HoangLong river basin Researchers have used the method of flood flow simulation byhydraulic model combined with the use of GIS tools to develop flood risk maps, Thehydrologic data collected were good and the topographic data using DEM 90x90 The

research used ID Model (Mike 11) to run the simulation of ereating flood flow To

‘modeling 2D flood process, there is an addition method that called FP4 of DHL hadbeen used

In an aggregate and systematic research, (de Moel, van Alphen, & Aerts, 2009) gave

an overview of existing flood mapping practices in 29 countries in Europe and shows

‘what maps are already available and how such maps are used The research indicatedthat flood maps existed in many different forms, but in general it was possible todistinguish between flood hazard and flood risk maps Flood hazard maps containedinformation about the probability and/or magnitude of an event whereas flood risk

‘maps contained additional information about the consequences (e.g economic

damage, number of people affected) Within these two general types, however, there

‘were different methods available to quantify hazards and risks, resulting in different

types of flood maps (Fig 22)

Figure 2.2: Conceptual framework for flood hazard and risk calculations (de Moel et

al, 2009)

To represent flood situation, there were many ways to ereate flood maps (de Moel et4l, 2009) had listed a total of six types of flood maps The maps selection depended onthe requirements and objectives of the research, because each type of map represented

4 different characteristic of flood and different use value The examples for all map

types were shown in Figure 2.3

risk (damage) maps (de Moel et al., 2009)

To assess flood damage to the study site takes into account the vulnerability of theprotected objects, researches often build the flood inundation maps (flood hazardmaps) and the flood risk maps The investigation and gathering socio-economic data todevelop flood risk maps often need more time Meanwhile in order to focus on thework of flood waming, forecasting and large flood management, the flood inundationmaps can relatively meet the requirements of the topic Thus in our research, theauthor chose to build flood elevation maps base on the distribution and elevation data

‘of 2009 historic flood marks This flood inundation map will help readers have spatial

‘view on the possibility of flooding in downstream areas of Tra Bong river basin,

1.3 Geographical Information Systems in Hydrology and Water Resources

In the past few decades, the development in computing capacity of the hydrology hascreated significant development about hydrology research and operationalGeographical Information Systems (GIS) introduction has created a powerful tool tocalculate and represent the result of hydrological studies, particularly floodmanagement From the first years of application of GIS in hydrology and waterresourrees management, (Singh & Fiorentino, 1996) in his esearch proposed a conceptfof GIS: GIS is a general-purpose computer-based technology for handling

‘geographical data in digital form It is designed to capture, store, manipulate, analyseand display divers sets of spatial or geo-referenced data,

In hydrology and water resources studies, concer for resource management, andenvironmental quality requires application of distributed models One characteristic ofdistributed-parameter models is that they are data intensive, Multiple data types such

as hydrometeorology, topography, land use, soils, geology, streamflows, ete are

commonly required Many of these data are often used to identify a hydrologic unit,because application of these models requires partitioning of the watershed intohomogeneous units This often proves to be cumbersome Therefore, a spatial dataanalysis and manipulation tool is desirable and that tool is GIS Furthermore, the use

‘of computers in hydrologic modelling has become so widespread that the marriagebetween hydrology and GIS should be a logical step (Singh & Fiorentino, 1996),

‘The book “Hydrological Applications of GIS” (Gumell & Montgomery, 2003) isbased on the art ‘was published in the journal "Hydrological Processes", focusing

‘on the main theme is the development and application of geographic informationsystem (GIS) to address hydrological problems, During the period before this book

‘was published, the application of GIS in environmental modelling has flourished,utilizing the advantages of linking the data types of GIS with models,

2

In the book, the article shows that the application of GIS in hydrology is very diverse,while scientists have simulated the hydrological processes through models, fromcommon characteristics model to semi-disuibuted and distributed model, Themanagement of water resources also followed a similar path to increase the spatialresolution, especially on infrastructure With both objects, GIS is useful because it canprovide data with higher resolution, Furthermore, various applications of hydrologyoften require controlled by multiple users by access to common databases, so thedatabase must support GIS, data quality and how management database is the "heart"for the development of GIS applications With the development of digital elevation

‘models (DEMs), the most widely GIS application in hydrology defining catchmentflow and overflow areas based on terrain types and combinations of them withhydrological models

In recent years, the development of GIS technologies and the availability of a widerange of geodatasets provide great potential for modeling physiographic objects, in

particular, river systems An adequate automatic method of river network digitizing

was developed and tested using the digital elevation model based on applying the

Complex Energy Index (CEI) in a research of (Gartsman & Shekman, 2016) The

index form includes the independent parameterization of main geomorphologic andclimatic factors of the first-order stream generation The method testing proved that it

is efficient and the results of ils use can be reliably interpreted, The software wasdeveloped on the basis of ArcGIS tools It implements the full algorithm of theautomatic digitizing of river network using the digital elevation model It also allowsadvancing studies on this issue, and could be used for solving applied problems,

Studies have shown that for each individual an effective tool to classify ma

characte ic according to the needs of the user Thereby, the user ean use the mathoperations to calculate between the attributes of each layer such as addition,subtraction, filtering, or set conditions (0 ereate a new map layer as required by thestudy These features of GIS are useful in creating flood maps of the study area, and

can also display maps in accordance with the standards for presentation and

publication,

1.4 Flood Forecasting

“The main objective of this research is to find a suitable flood forecasting method for

‘Tra Bong river basin, in terms number of monitoring stations is limited ‘Thus, thereview part first examined flood forecasting methods that used in studies in thecountry and around the world to get the most appropriate selection In a report of 4project implemented by Cooperative Research Centre for Catchment Hydrology’ ~

‘Australia (Srikanthan, Elliott, & Adams, 1994) had determined that flood forecastingmethods were considered under two separate heading: rainfall ~ runoff methods andfood routing methods

+ Real-time forecasts of discharge, obtained by rainfall ~ runoff modelling, aregenerally less accurate than those obtained by channel routing of hydrographobserved at an upstream gauging site, However, real-time forecasting methodsbased on rainfal-runoff modelling are necessary because in head water basins there

is no alternative since upstream stations do not exist and in some circumstances itmay yield forecasts with greater lead time

+ The second heading, flood routing methods offers a satisfactory means of floodforecasting for long river systems For this approach to be successful, the travel time

of flood peaks from upstream to the downstream site needs to be long enough toallow adequate period of warning Flood routing can be classified under four

‘methods: experience methods, statistical methods, hydrologic routing methods andhydraulic routing methods

Using three of four methods mentioned above: statistical, hydrologic and hydraulicrouting methods, (Duy Kieu, 2012) in his engineering doctoral thesis had researched afully procedure of flood forecasting, in order to develop large flood managementsolutions for Lam river basin By statistical methods, he used the rainfall data,discharge and water level data of all monitoring stations in the basin to analyse thecharacteristic of climate and flood in study area Thereby, he analysed and partitioned

ng GIS tools,

he could develop the flood hazard map for the basin And applied hydrologic andflood tisk by give marks for each sub-area of study site, By the mark

2B

hydraulic routing methods, he simulated the flood flow in the river using somehhydraulic modern as Mike 11, Mike Flood and Mike 21 In this study, the combination

Of flood problems in the river system is very complex But the author stated that

crucial to develop the scientific basis for flood management, Research of thecombination of large flood on Lam river system in this thesis focused on the issues:The total amount of flood water, the meeting of the flood peak flow from upstream todownstream, the origin of flood water from tributaries

Referring to the role of flood forecasting in flood management, (WMO, 2011)indicated that flood forecasting is a necessary part of flood management, given that no

preventative or defence measures can be completly effective, The reality of economic

limits to the provision of defences, together with the possibility that the capacity ofdefence systems may be exceeded or they may be fal require that other measures are

in place Provision of flood forecasting will also form part of flood managementplanning and development strategies, which recognize that there are occupied flood

plain areas where noi structural measures cạn be effective

This manual book provided the basic knowledge and guidance to develop or to set up

an appropriate and tailored system for any case in which a flood forecasting and

forecasting and waming system, which were:

+ Design of a flood forecasting system:

+ Implementation and operation ofa flod forecasting sytem;

+ Flood warnings:

+ Training

The Manual does not set out a step-by-step process for the design of a flood

forecasting and warning system along the lines of a particular template or practice inany one country Rather, presented in all chapters are a number of examples ofdifferent practices and technologi „ whíeh may reflect different levels of

development, ranges of needs and also capacities in a number of different situations

Figure 2.4: Process for developing a flood forecasting model (WMO, 2011)

An important part of flood forecasting is selecting a model or a combination of models

to appropriately simulate the process of flood flow forming in the river Hydrological

tainfall-ranoff model was selected to simulate the process of flow forming from rainfall data when the data series of flow, water levels and the number of monitoring

stations is limited The structured and methodology of this model are relatively simple,

0 it is not too complicated to select

To simulate the flood flow process in downstream river, a 1-dimensional (1D) or

2-dimentional (2D) hydraulic model can be used 2D model gives more detailed results,

provide ä full range of information on flood flows, such as flood depth, flood velocity,

flow direction, and easily combined with GIS tools to build the flood maps However,

25

2D model has high requirements on the input data, especially the topography data(usually high resolution DEM): time to run the model is long and requires a verytrong computer, In contrast, ID model gives less detailed results and more difficultycombined with GIS tools for flood mapping, but it requires less data over terrain (justthe data in the cross sections) and model running time is shorter In addition, to buildflood maps from the results of ID model, an additional tool (such as FP4 of DHI) can

be used to interpolate and make a map with obtained results of cross sections

To evaluate the efficiency of 1D and 2D numerical models for predicting river floodinundation, (Horritt & Bates, 2002) had made a test for 1D and 2D models of floodhydraulics (HEC-RAS, LISFLOOD-FP and TELEMAC-2D) on a 60 km reach of theriver Severn, UK The three models are calibrated, using floodplain and channel friction

as free parameters, against both the observed inundated area and records of downstreamdischarge About the data, while HEC-RAS only needed 19 cross

channel, LISFLOOD-FP and TELEMAC-2D, the two 2D models need 50 m resolution

ections of the

DEM for entire channel and floodplain The results of testing are shown in Table 2

Table 2.2: Predictive performance of the 3 models using independent

calibradon/validation data (Horritt & Bates, 2002)

mo ae oma Sasi Niwa nae

In terms of simulations where the model is calibrated against the data it is trying topredict, the TELEMAC-2D marginally outperforms the other models apart from forthe 2000 inundation data In terms of predictive performance (Table 3), the HEC-RASmodel performs marginally better (on 4 out of 6 measures) than TELEMAC-2D,Which shows the best performance on 2 measures and which in tum is better thanLISFLOOD-FP LISFLOOD-FP does not produce the best predictive performance forany of the measures and data sets examined It is interesting to note that TELEMAC-2D is better when predicting inundation for the 1998 than HEC-RAS, the situation

being reversed for the 2000 event, Finally, study had shown that HEC-RAS was thebest model in flood forecasting for the study area, with the terms and requirements,stated in the study.

mm

CHAPTER II: APPROACH AND METHODOLODY:

2.1 Approach of study

Study area

2009 flood marks Data collection for

data collection flood simulation

GS tools Hydrotogical model

DEM data Hydrodynamic

model

2009 flood Flood forecasting

inundation map experiment

Assessing the large floodwe the large

situation of the basin |<

rivers of Quang Ngai province in particular, he characteristis are short and steep, fast

"water concentration time, so these maps are fastest sources of information about thescale and level of flooding, helpful for the evacuation if flood occurs

Flood maps give the information about flooding and non-flooding areas, and also flooddepth Bes les information on flooding, the maps also show the other basicinformation, such as landmarks, transportation, resident, river Flood maps give alook from general to detail of flood level due caused by a certain rain,

‘The use of flood maps built with some certain scenarios will help local authoritiesactively choose quickly responding plan when there are information of flood warnings,fore: isting in downstream,

Flood maps can be classified into two basic types:

‘+ Type 1: Map of flooding caused by historic floods, such as the floods in 1999, 2009

‘These maps show the real information of flood so they support the most objectiveforecasting and warning for residents These maps are also useful in verifying theother flood forecasting maps,

flood on this map is the designing floods with the designing frequency of rainfall ordischarge in one or more gauging stations on rivers, In other words, these maps arebuilt according to some rainfall - flood scenarios

‘Most of the flood maps are built on the basis of the depth of flooding Flooding depth

‘was determined from the flood water elevation and the ground level, by the following

1)

‘iy: Flood water depth;

H„: Flood water elevation;

?

Z,: Ground level

Hy, and Z, have to be in the same elevation system,

Thus, it is necessary to determine the flood water elevation and ground level at aposition to calculate the depth of flooding there,

Ground level is developed from topographic maps or field measurements, It is a very

important factor for the flood simulating results, The detail level of the topographydata will affect the accuracy of flood risk maps TẾ the point density is not enough, notfully reflect the transformation of the terrain, the floodplain simulations will not bedetail Some non-flooded areas with be considered as flooded areas, and vice versaFlood water elevation: For a map of the historic flood, flood water elevation can beinvestigated and actual measured, combined with computational methods For floodmaps corresponding to a designing flood, the flood water elevation is determined bythe hydrological and hydraulic model, and calibrated, verified with observed data

2.3 Using GIS tools to develop flood inundation map

‘The inundation map will be created using Are GIS software (Are Map 10, ESRI Inc hup:/wwwesti.com) Are GIS includes very powerful tools to manage, update,analyse and publish information to create a complete Geographic Information System

-G S), allows

+ Create and edit integrated data (integrating spatial data with attribute data), allowsuse a variety of different data formats and even the data obtained from the Internet,+ Query spatial data and attribute data from many sources and in many different ways;

+ Display, query, and analyse spatial data in combination with atribute date

+ Create specific map and print with professional presentation quality

‘Topographic data used in Arc Map as contours or elevation grid - often presentedsquare or triangle grid (this type of data is often called the Digital Elevation Model ~DEM),

Figure 3.2: Digital Elevation Model with square grid (Source: www.p-gis.com)

Geographic factor determines the classification of submerged depth, For example, if

the terrain with contour different is 5 m, the flooding boundary determining will have

significant error because some small but high zones between 2 contours can be

considered to be flooded, but in fact, itis not, (Figure 3.4)

Zone considered+ foods

Figure 3.3: Error based on the topography data in flood inundation mapping.

In general, to reduce errors due to topography data, elevation difference between thecontours should be equal to or less than the elevation difference of submerged depth

lassification, If submerged depth classification is Im, the corresponding elevation

difference between the contours on the map is also Im or less According to Vietnam

31

Construction Standard 309:2004: Surveying work in building construction - Generalrequirements (MOC, 2005), the scale of map must be 1/2000 or more detail (Table3.1), The average point density for creating @ 1/2000 seale topography map is from 12

to 18 points per a 1 square km

‘Table 3.1: Technical requirements of elevation different between the contours with

24 Mike 11 general description

The research focuses on using rainfall ~ runoff hydrological model (NAM) and dimensional hydrodynamic model (Mike 11 HD) to simulate the process of flow

1-forming in the basin, thereby experiment flood forecasting for Tra Bong river system.

‘The Theoretical Foundations of rainfall ~ runoff hydrological model (NAM) and dimensional hydrodynamic model (Mike 11 HD) are mainly based on the book Mike

1-11 A Modelling System for Rivers and Channels Reference Manual (DHI, 2007).MIKE 11 is a professional engineering software package for the simulation of flows,

‘water quality and sediment transport in estuaries, rivers, irigation systems, channels

and other water bodies MIKE 11 is a user-friendly, fully dynamic, one-dimensional

‘modelling tool for the detailed analysis, design, management and operation of bothsimple and complex river and channel systems, With its exceptional flexibility, speedand user friendly environment, MIKE 11 provides a complete and effective designenvironment for engineering, water resources, water quality management and planningapplications (DHI, 2007)

“The Hydrodynamic (HD) module is the nucleus of the MIKE 11 modelling system andforms the basis for most modules including Flood Forecasting, Advection-Dispersion,

Water Quality and Non-cohesive sediment transport modules The MIKE 11 HD

module solves the vertically integrated equations for the conservation of continuityand momentum, ie the Saint Venant equations

Applications related to the MIKE 11 HD module include:

+ Flood forecasting and reservoir operation

+ Simulation of flood control measures

* Operation of irrigation and surface drainage systems

‘sign of channel systems

+ Tidal and storm surge studies in rivers and estuaries

‘The primary feature of the MIKE 11 modelling system is the integrated modularstructure with a variety of add-on modules each simulating phenomenon related to

river systems

In addition to the HD module described above, MIKE 11 includes add-on modules for:+ Hydrology

* Advection-Dispersion

+ Models for vatious aspects of Water Quality

* Cohesive sediment transport

+ Non-cohesive sediment transport

MIKE 11 HD applied with the dynamic wave description solves the verticallyintegrated equations of conservation of continuity and momentum (the Saint Venantequations), based on the following assumptions:

+ The water incompressible and homogeneous, ic negligible variation in density+ The bottom-slope is small, thus the cosine of the angle it makes with the horizontal

may be taken as 1

3ã

+ The wave lengths are large compared to the water depth This ensures that the floweverywhere can be regarded as having a direction parallel to the bottom, ie verticalaccelerations can be neglected and a hydrostatic pressure variation along the verticalcan be assumed

+ The flow is subcritical

Continuity equation (Conservation of mass, in differential form):

aoe Xa)» = panar= pCa at G8)

‘Momentum equation (Conservation of momentum, in difference form):

Weight:

‘The solution of the equations of continuity and momentum is based on an implicitfinite difference scheme developed by Abbott and Ionescu (1967) In this scheme, thewater level and flow along the tributaries are included in a system of alternating gridpoints as shown in Figure 3.6 below

Figure 3.6: Channel section with computational grid

‘MIKE-11 can handle multiple branches, and in the tributaries where branches meet Anode will be created in which the water level is calculated The shape of the

‘computational grid around a node which has three branches is shown in Figure 3.7,

Figure 3.7: The shape of the computational grid around a node which has three

branches

35

‘The shape of the grid points and the nodes of a complete model is shown in Figure 3.8

It should be noted that at the border point, we create a node that we can calculate thewater level

In a grid point, the relationship between the variables Z, (including water level h, orgrid poi p 1s

flow discharge Q,) at the point and at the neighboring point is represented by using a

linear equation as follows:

“27 + BZ) + 72h G.10)

‘The lower number in the equation (3.10) indicates the position along the branch andthe upper one indicates time step The coefficient ơ, 8, 8 and y in the equation at thepoint h is calculated by implicit finite difference of continuity equation and at the point

Q by implicit finite difference of momentum equation

At every grid point along the brane the equation (3.10) is created Suppose a branchhave m grid points; 1 is an odd number, so the beginning and end grid points of abranch is always water level points (h) This creates n linear equations with n + 2unknowns, Two additional unknowns are from the equations created at the beginningand the end points, at which Z, and Z2, respectively become the water level values at

the node, whereby the upstream and downstream of the river are connected The

following section describes the linear equation:

aft + BH + QM = 8 ate + BQ + yA" = 3, a,Q3" + pny" + +” = 8 ans + B,QI + yhệ” 84 a,Qy" + Bh" + G5" = 8,

AINE Tie h Gan

ty 205s + Đa¿hg2 + Yn2Qht = Sn2

GD + Base + Tạ” = By

qUO” + Bh + HE = 8

iin the frst equation and Hạ, in the final equation are respectively water levelk inthe upstream node and the downstream nodes In MIKE-11 compatibility water level isapplied at nodes, that means the water level atthe frst point of @ branch is the waterlevel at the node, whereby the last section of upstream of the branch is connected Inother words, hy = Hyy This means a=, BI=1, 71=0 and 61=0 Similarly, the final{rd point at which hy = Hds so on=0, Bi

using standard reduction algorithm, However, because MIKE-11 can handle multiple

branches so we apply another approach To explain this problem, the above equation

‘will be presented in the matrix in Figure 3.9

Figure 3.9: Branch matrix before reducing

By standard reduetion method, we can change matrix 3.9 to matrix 3.10

By this matrix we can see at any grid point, variable Z (water levels and flows) isrepresented as a function of the water level and flow upstream and downstream:

Zp" =e,~a,HẸT -b,HS" 6.12)

@ weesse

— Riverbank:

hụ Water level in poi)

trachX Tre Comertine

Figure 3.11: Three-branch node with limit for continuity equation

At the intersection, the continuity equation at the node of intersection is set as follows:

In equation (3.14) Q2 Qui and Qc; at time level n+1 can be substituted according

to equation (3.13), and the following equation is achieved