Next Article in Journal
Studying Unimodal, Bimodal, PDI and Bimodal-PDI Variants of Multiple Soil Water Retention Models: I. Direct Model Fit Using the Extended Evaporation and Dewpoint Methods
Next Article in Special Issue
Attitude and Actual Behaviour towards Water-Related Green Infrastructures and Sustainable Drainage Systems in Four North-Western Mediterranean Regions of Italy and France
Previous Article in Journal
Hydraulic Conductivity Estimation Using Low-Flow Purging Data Elaboration in Contaminated Sites
Previous Article in Special Issue
Analysis of Hyetographs for Drainage System Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Urban Flood Prediction Using Deep Neural Network with Data Augmentation

Department of Civil Engineering, Kyungpook National University, 80 Daehak-ro, Buk-gu Daegu 41566, Korea
*
Author to whom correspondence should be addressed.
Water 2020, 12(3), 899; https://doi.org/10.3390/w12030899
Submission received: 11 February 2020 / Revised: 18 March 2020 / Accepted: 19 March 2020 / Published: 22 March 2020

Abstract

:
Data-driven models using an artificial neural network (ANN), deep learning (DL) and numerical models are applied in flood analysis of the urban watershed, which has a complex drainage system. In particular, data-driven models using neural networks can quickly present the results and be used for flood forecasting. However, not a lot of data with actual flood history and heavy rainfalls are available, it is difficult to conduct a preliminary analysis of flood in urban areas. In this study, a deep neural network (DNN) was used to predict the total accumulative overflow, and because of the insufficiency of observed rainfall data, 6 h of rainfall were surveyed nationwide in Korea. Statistical characteristics of each rainfall event were used as input data for the DNN. The target value of the DNN was the total accumulative overflow calculated from Storm Water Management Model (SWMM) simulations, and the methodology of data augmentation was applied to increase the input data. The SWMM is one-dimensional model for rainfall-runoff analysis. The data augmentation allowed enrichment of the training data for DNN. The data augmentation was applied ten times for each input combination, and the practicality of the data augmentation was determined by predicting the total accumulative overflow over the testing data and the observed rainfall. The prediction result of DNN was compared with the simulated result obtained using the SWMM model, and it was confirmed that the predictive performance was improved on applying data augmentation.

1. Introduction

It is very important to predict the possibility of flooding in urban basins and ascertain the degree of flooding in advance. Hazardous flooding events that cause residential and commercial property damage can occur in many forms and with frequent frequencies. Numerical analysis, trend analysis, and flood prediction for various rainfall scenarios can be utilized as important basic data for urban planning or flood response. One- and two-dimensional flood analysis programs can be used for this purpose; however, in the case of a numerical-based model, it may take some time to adjust parameters, collect data, and perform post-processing of output data.
For this reason, data-driven models using machine learning have recently been applied to the prediction of urban runoff [1]. Granata et al. [2] predicted runoff in urban regions via support vector regression (SVR) and compared it with the results of EPA-SWMM simulations. Although there was a lack of accuracy in estimating the peak runoff, the possibility of improvement was confirmed through verification. If more data are applied to the learning of SVR, the prediction power of peak flow will be further enhanced. The study of input data for a machine learning method is also important to find a reliable prediction model. Talei et al. [3] suggested the evaluation result of rainfall and discharge inputs with adaptive network fuzzy inference systems (ANFIS). In this study, the criteria for determining the optimal rainfall-runoff analysis model among sequential, pruned sequential, and non-sequential time-series input data were presented. In addition to a general artificial neural network (ANN), the deep learning (DL) technique has begun to be applied to the water resources field to improve the predictive power of hydrologic data and to include more ideas in the model. According to Shen [4], the application of the DL technique to the water resources and hydrology fields is becoming more common, and it is likely to yield good predictions and analysis results in the water resources field as well as in most scientific fields. Hu et al. [5] used a long short-term memory (LSTM)-based DL model for rainfall-runoff simulation and used data on 86 rainfall-runoff patterns. This research was possible because the database for LSTM could be collected by hydrologic observations in the target river basin.
Although many studies have been conducted to predict rainfall, hydrologic data, and flood events using ANN and DL techniques, research related to the supplementation of input data for DL is deficient. Indeed, some studies have raised the problem of insufficiency data for application to a data-driven model. Li et al. [6] raised the problem of the lack of input data and attempted to predict the possibility of flooding in urban regions. For flood forecasting, a conceptual model considering the drainage system was combined with the logistic regression technique. The hydraulic model was compared with the prediction results, and the verification results that could be developed were shown. Nikhil et al. [7] analyzed a dam basin that had been recently affected by a disaster, which was caused by dam outflow. The dam outflow could not be predicted due to the lack of data. Rainfall, flow discharge, and water level data were used to predict the flood downstream, and as observation data were added to the ANN, the predictive power was shown to decrease. It appears that the performance of a neural network becomes degraded when using uncertain observational data when data are insufficient. However, this study does not provide a sure solution for a lack of data.
In this study, urban flood prediction was performed using a deep neural network (DNN), and the problem of a lack of data was resolved by data augmentation and observed rainfall data that caused the urban flooding. The total amount of overflow was predicted via the trained DNN for the drainage basin of Samseong-dong, Seoul, which was damaged by heavy flooding during the years of 2010–2011. The input data for DNN learning consisted of 10 statistical characteristics, including total rainfall, rainfall intensity, hourly maximum rainfall, standard deviation, skewness, and kurtosis. The total accumulative overflow calculated by SWMM simulations was selected as the target value. To apply abundant data to the predictive model, data augmentation was applied. The predictive power according to the amount of augmented data was evaluated to estimate the applicability of the suggested methodology. In this paper, the effectiveness of data augmentation for the urban flood prediction is discussed. In addition, this paper attempts to confirm the optimal input combination by using statistical characteristics of rainfall distribution. The proposed methodology can supplement the lack of data on actual rainfall that has caused flooding and establish an optimal database that can be applied to DL techniques that require various data patterns. Finally, a prediction model is selected to predict the exact total accumulative overflow for the study area.

2. Methodology

2.1. Observed Rainfall

Heavy rainfall that caused severe flooding damage to cities (Seoul, Incheon, Cheonan, Chungju, Busan, Ulsan, Changwon city area) was investigated. For collecting rainfall data, the Automated Synoptic Observing System (ASOS) and observed data from the Automatic Weather System (AWS) were obtained from the Korea Meteorological Administration (KMA) website [8]. The actual rainfall from 2009 to 2017 was investigated and a total of 70 heavy rainfall events were collected. The duration of rainfall was considered to be 6 hours, since most urban floods occur within 6 hours [9]. The total rainfall ranged from a minimum of 7.5 mm to a maximum of 284.4 mm , and the rainfall that did not cause flooding was also considered for input data of the DNN. By collecting data on observed cases of rainfall that induced urban flood, it is possible to compensate for the lack of rainfall data in the study area. Table 1 shows the names and codes of observatories, the dates of observation, and the amounts of rainfall employed in this study.

2.2. Urban Runoff Simulation

EPA-SWMM, which is used to calculate the total accumulative overflow in the study area, is an urban runoff model that can be used for both rainfall–runoff and river routing, such as the analysis of surface and underground flows, and flows in the drainage pipe network in urban basins with drainage systems [10]. Prior to analysis of the scope of flooding for urban basins, the runoff for the rainfall was calculated using SWMM (version 5.1, United States Environmental Protection Agency, Washington, DC, USA), and operation of the pump station was applied with reference to the report of the Seoul Comprehensive Planning System for Flood Damage Reduction (2015). RUNOFF and EXTRAN were used among the execution blocks of the SWMM model, and the RUNOFF block performed initial calculations for the outflow, water quality, and inflow hydrographs for drainage basins based on the rainfall scenario [11]. The EXTRAN block was used to calculate the flow rate and depth of the drainage pipe system using the output data of the runoff block, and it was possible to analyze the backflow and overflow amount in the pipe based on the flow rate and the water level in the drainage network for each calculation [12]. For analysis of the overflow according to various rainfall events in the study area, the Saint–Venant equations (Equations (1) and (2)) were used.
Q = W × 1 n ( d d p ) 5 3 × S 1 2
Q t + g A S f 2 V A t V 2 A x + g A H x = 0
Here, Q is runoff ( m 3 / s ), W is the subwatershed width ( m ), n is the Manning’s roughness coefficient, d is the depth ( m ), d p is the ground reservoir lost depth ( m ), S is the subwatershed slope, A is the surface flow cross-sectional area of sub-watershed ( m 2 ), and V is the surface flow velocity ( m / s ). One-dimensional runoff analysis results of SWMM were used as the target value of the DNN model, although they had some limitations [13]. Because there is no observed flow data in the study area, simulated SWMM results were considered as a reasonable runoff from rainfall events.

2.3. Deep Neural Network

The DNN is the most basic DL technique that can be used for DL by using two or more hidden layers in a general ANN. The basic structure of a DNN is shown as Figure 1a, where x 1 , x 2 , …, x n represent the input data with n attributes, and y ^ represents the value predicted through the DNN. When the neural network layer is deepened, it can have a high level of expressive power and is strong in terms of elaborate learning and expression efficiency. The DNN can require large amounts of data and must be supported by high-performance computing technology and data storage capabilities. DNNs are also composed of an input layer, a hidden layer, and an output layer, and various activation functions can be applied to ensure that input data is completely reflected the hidden layer [14].
The rectified linear units (ReLU) function as shown in Figure 1b was used as the activation function of hidden layers. As a result, if x is greater than 0, the derivative value is 1, and even if input data go through several hidden layers, the characteristics of the data remain without disappearing to the output layer [15]. Adam used for optimization for the learning process of DNN, which is effective for dealing with nonlinear problems, including outliers [16]. In this study, these parameters could improve the accuracy and reliability of the learning process.
For the error function, the mean absolute error (MAE) shown in Equation (3) was used to estimate the loss rate for learning. Here, n is the number of target data and, in this study, it represents the total number of total accumulative overflows to be predicted. y is the amount of total accumulative overflow calculated using EPA-SWMM, and y ^ is the total amount of accumulative overflow predicted by the DNN. In order to make up for the shortcomings of MAE’s size-dependent error analysis, the mean absolute percentage error (MAPE) was also used. The MAE and MAPE were applied in the error analysis to verify the performance of the DNN.
MAE = 1 n | y ^ y |
MAPE = 1 n | y ^ y y | × 100 ( % )

2.4. Data Augmentation

Data augmentation is a method of increasing the quantitative amount of data by finely tuning the initial input data. It is effective for neural network training to make invariant predictions [17]. In other words, some changes are made to specific data x to obtain x and this is treated as new input data. In general, it is also an effective method to uniformly apply random noise along a Gaussian distribution [18]. In this study, it is necessary to build a DL model with only observed rainfall scenario and so, it was determined that the application of data augmentation is effective. Therefore, data augmentation was attempted on the input data. The standard deviation for the properties (total rainfall, maximum rainfall by duration, rainfall intensity, etc.) of various samples (rainfall scenarios) was multiplied by C (constant number) and N ( 0 ,   1 ) (the noise of a Gaussian distribution). Then, the new set of input data was constructed by adding multiplication results to the initial input data. In Equation (5), x represents the original data, σ i ( i = 1 ~ 9 ) is the standard deviation of each 9 data feature, noise is a random number between 0 and 1 generated in the Gaussian distribution, x represents newly generated data, and this was used as new input data. The constant number of C is user-defined hyper-parameter and it was empirically specified as a 0.05 value. The reason for using 0.05 in C is that when the data augmentation technique is applied to the observed rainfall, it changed less than 10 mm . There are not many examples of the application of data augmentation in the water resources field, and the focus should be on confirming the validity and reliability of this technique. The final goal was to increase the accuracy of the prediction of the total overflow rate by linking data augmentation and DL.
x = x + C × σ i × N ( 0 ,   1 )
The flowchart for applying data augmentation is shown in Figure 2. The database of rainfall and the total accumulative overflow was constructed by simulation results of SWMM which validated by the usage of FLO-2D (two-dimensional flood analysis program) and flood marks. The DNN was trained using the initial database, and the prediction of total accumulative overflow was performed on the new input data, which came out of the data augmentation. An outlier check for the predicted results was performed by comparing with SWMM simulation results. Excel was used to check whether the negative prediction result was shown, and the unacceptable value was predicted compared to the simulated results. Once the data were confirmed to be sound, the processed input data and the predicted total accumulative overflow data were added to the database for DL. The proposed series of processes was repeated, and this study attempted to identify how many rounds of data augmentation would be best according to the input data combination.

3. Application

3.1. Study Area

The subject area was set to the Samseong-dong district including Samseong 1 and Samseong 2 drainage basins in Seoul, Korea. Figure 3a reveals the drainage basin borders, drainpipe network, and manhole locations. Figure 3b shows the nearby AWS stations around the study area. The Gangnam AWS station in the Samsung-dong basin. Seocho and Songpa AWS stations are 4.1 km , 4.6 km from Gangnam Observatory, respectively. The drainpipe network and manhole information were used in performing SWMM-based one-dimensional urban runoff analysis. The rainfall on the urban surface flows into the Tancheon stream through a conduit network. The subject area experienced extreme flood damage because of torrential rain on 21 September 2010 and 27 July 2011. It appears that an urban flood occurred because of the relatively low land, complicated drainage system, and inflow of rain exceeding the conduit capacity. According to Seoul City’s Storm and Flood Damage Reduction Master Plan [19], urban flood-prone areas (red section in Figure 4) in Samseong-dong were Samseong-Seonreung Stations (around Teheran street, Figure 4a) and Bongeunsa-Samseong Stations (Yeongdong street, Figure 4b). It is likely that in the event of torrential rain exceeding sewage conduit capacity, a flood might have occurred in the two regions described above. According to the inundation trace map in 2010 and 2011, a flood took place along the street between Samseong Station and Seonreung Station.

3.2. SWMM Validation with Flood Trace Mark

To decide the adequacy of the one-dimensional urban runoff analysis results, this study performed two-dimensional inundation analysis, using overflow at each manhole point and attempted to compare the results with those for flood-prone areas and flood marks. It is because there are only records of flood marks and data on the water level and discharge in the conduit are absent. To validate the adequacy of the SWMM model results, the two-dimensional model adopted the use of the finite difference method (FDM)-based hydraulic analysis program ’FLO-2D’. FLO-2D is a grid-based hydraulic model designed to analyze flood wave propagation approved by the Federal Emergency Management Agency (FEMA) [20]. If the overflows of each manhole were calculated from SWMM, these were entered into the input file of FLO-2D. The exact location of a manhole in the two-dimensional space was found by using the GIS program (spatial join tool in ArcGIS). The scope of calculation on the target basin was examined through an expected flood map and flood history records. After completing the construction of two-dimensional grids, flood waves were analyzed, using the continuity and momentum equations.
Sufficient time was assigned for simulation, delivering flood waves between buildings and to the road as much as possible. To properly reflect the influence of buildings and roads on the flood waves, two-dimensional inundation analysis was performed on 5 m2 grids. Using Equations (6) and (7), the synthetic roughness coefficient was calculated, and the result (0.025) was applied. Here, n is the synthetic roughness coefficient and n 0 is the bed roughness coefficient. In addition, θ refers to building coverage (%) while n 1 , n 2 , and n 3 are 0.06 (farmland), 0.047 (road), and 0.05 (others), respectively. Furthermore, A 1 , A 2 , and A 3 denote the agricultural area, road area, and other land-use areas, respectively, while h (m) refers to the depth of water [21,22]. Considering that the study area is a metropolitan watershed, only the roughness coefficient of the roads ( n 2 ) was used in Equation (7). In addition, this study applied 0.014 as the roughness coefficient of the roads by referring to the road design manual [23].
n 2 = n o 2 + 0.02 × θ 100 θ × h 4 / 3
n o 2 = n 1 2 A 1 + n 2 2 A 2 + n 3 2 A 3 A 1 + A 2 + A 3
As shown in Figure 5, the flood mark of 2010 and 2011 was compared to the two-dimensional flood analysis results applying observed rainfall for 21 September 2010 and 27 July 2011. Flood marks were prepared by Seoul city. The hydrograph of the SWMM result (10 min interval) for each observed rainfall was indicated in Figure 5b,d. The ratio of inundation was calculated in 5 × 5 grid units. The percentage of simulated inundation grid within the flood mark was 75.8% and 71.2% for 2010 and 2011 rainfall events, respectively. In addition, 67% and 63.4% of the two-dimensional simulation results for 2010 and 2011 events occurred within the flood mark from Samseong Station to Seonreung Station. Because a flood trace map is prepared at 0.3 m or deeper, and Figure 5b,d represents the flood depth through surface flow analysis results reflecting roads and buildings, it is concluded that two-dimensional flood analysis and SWMM simulation results in the study area are acceptable.

3.3. Input Data

Because the observed heavy rainfall is diverse in terms of the shape of distribution and location of the peak, it is hard to predict urban flood in real-time. It appears that runoff patterns in the urban basins are influenced by rainfall distribution characteristics. Based on 70 observed rainfall events, this study calculated the total rainfall, maximum 1-hour rainfall, maximum 2-hour rainfall, maximum 3-hour rainfall, rainfall intensity, peak rainfall position(%), standard variation, skewness, inter-event time, and kurtosis. Using the same statistical features of rainfall events as those of the DNN input data, this study attempted to build a total accumulative overflow prediction model in which spatiotemporal characteristics of rainfall were learned. The maximum rainfall for a duration of 1–3 h represented the maximum value of ordinate with 1-, 2-, and 3-hour durations while rainfall intensity was estimated by dividing the 6-hour duration of the total rainfall. The peak rainfall location represents the location of the peak value of the rainfall distribution in percentage. The inter-event time is obtained by estimating the time, which represents 0 rainfall events in minutes. The results of statistical analysis on a part of the collected rainfall event data are stated in Table 2 and Table 3.
The statistical characteristics on rainfall event data were entered into the DNN. In this study, four different cases of input data sets were applied (Table 4). Then, a correlation analysis of the total accumulative overflow of statistical characteristics was performed. It is estimated that the correlation between total rainfall and total accumulative overflow is relatively higher than other characteristics. In terms of a peak rainfall position (%), the correlation converged to 0, not being used in combination with input data. In addition, a relatively high positive correlation with total accumulative overflow was found at maximum 3-hour rainfall and rainfall intensity as well. In terms of skewness, kurtosis, and inter-event time on observed rainfall, a negative correlation with the target value was observed. Case 1 reveals the use of all statistical analysis values excluding a peak rainfall position while Case 2 refers to a combination of input data without two sub-factors among characteristics in negative correlation. Case 3 uses five data sets in a high positive correlation. Lastly, Case 4 uses the highest positive correlation data and three of the negative correlation data. Regardless of the combination of input data, the target values were the same as the total accumulative overflow.

3.4. Prediction Model and Data Augmentation

A total of 69 rainfall events were used for DNN learning except for the prediction target rainfall among a total of 70 observed rainfall events. The DNN was composed of one input layer, eight hidden layers, and one output layer, and three hidden layers were set to have 18, 18, 12, 12, 12, 10, 10, and 8 nodes, respectively. The number of hidden layers and the number of nodes of each hidden layer were determined through trial and error. The final structure of DNN was determined by the empirical process and the MAE of the testing process. The summarized process of finding the optimal structure of the DNN is shown in Table 5. The activation function of the hidden layer uses the ReLU function. In the case of the activation function for output nodes, the linear function was used to derive the results in units of floating decimal points. The epoch for maximum learning was set to 5000. In DNN learning, 70% of input data was used for training, 20% was used for validation, and 10% was used for testing. The data for training, validation, and testing were randomly selected in 69 rainfall events. This is to avoid over-fitting a particular data set. As another way to avoid over-fitting, the early stopping function was used. If the MAE value for validation does not decrease even after repeating the DNN training, the learning process was stopped early.
The MAE with training and validation process is shown in Table 6. Table 7 shows the predictive results for the testing data through MAE and MAPE analysis. In order to show the overall error variation according to the application of the data augmentation, the number of data augmentation was applied up to 10 times. The reason for doing data augmentation up to 10 times is that too much data augmentation can make the original data meaningless. The dataset for DNN learning is increased by applying noise to the original data. In addition, gaps in the original data can be filled with noise-added data through the data augmentation technique. For these reasons, the DNN used in this study learns a more abundant rainfall-runoff pattern by data augmentation more than using only simulated results. Considering that abundant data is a fundamental condition of deep learning models, it is expected that the predictive power of the suggested model would be improved. As data augmentation was applied, the training, validation, and testing error were decreased (Table 6 and Table 7). In the case of the MAE, the overall decrease was observed without significant variation. However, in the case of the MAPE, which shows the error ratio of the observed value and the variability appeared to be large depending on the number of times data augmentation was applied.
Figure 6 shows the results of DL training for each combination of the input data after applying the initial input data and data augmentation. The horizontal axis of each graph represents an epoch, and the vertical axis represents the mean absolute error (MAE). The red dotted line shows the graph for the validation data set, and the blue solid line shows the error graph for the training data set. Each graph shows that after applying data augmentation, predictive power is improved not only for the training data set for the DL model, but also for the validation data set as well. In other words, the data augmentation technique improves the analysis result of validation data, rather than simply over-fitting with the training data.

3.5. Prediction Results

The result of SWMM simulations with consideration of the observed rainfall on 27 July 2011, was used as the target data for total accumulative overflow prediction (Figure 7). The DNN trained with 69 rainfall scenarios was used to perform the prediction. The total rainfall for the target rainfall was 184.5 mm , the rainfall intensity was 30.75   mm / h , the 1-, 2-, and 3-hour maximum rainfall was 71.0, 105.5, and 143.0 mm , respectively, the standard deviation was 4.38 mm , the skewness was 1.27, the kurtosis was 1.23, and the inter-event time was 10 min. As a result of EPA-SWMM simulations, the total accumulated overflow was calculated as 144.7 m 3 / s . Table 6 shows the predictive results before applying data augmentation and the predictive results of the DNN according to the number of data augmentations applied. When the data augmentation was applied four times for the input data of CASE 1, seven times for CASE 2, ten times for CASE 3, and once for CASE 4, predictions were the closest to the results of SWMM. The absolute average was calculated for the correlation values between the target value and the statistical characteristics in each input data (refer to Table 4). They were calculated as 0.4395, 0.5366, 0.5715, and 0.3877 in CASE 1, CASE 2, CASE 3, and CAS 4, respectively, and the higher the mean of the correlation coefficients was, the lower was the standard deviation of the predictive results according to data augmentation. The absolute average correlation is shown in Table 8.
In the case of CASE 4, which had the lowest absolute mean value of the correlation coefficient, the difference with the results of SWMM was confirmed to increase as data augmentation was applied. In contrast, in CASE 3, which had the highest mean value of the correlation coefficient, it was shown that predictive power could be increased by attempting data augmentation. In Table 8, * (asterisk) indicates that the predictive results of the DNN for each input data combination was the closest to the results of SWMM.

4. Conclusions

In general, there is a lack of historical data on flooding for a single urban drainage basin and basic data for constructing a data-driven model. In this study, a DNN was used to predict the total accumulative overflow in an urban drainage basin. Cases in which significant damage was caused by flooding in the target area were limited to the events of 21 September 2010, and 27 July 2011, and therefore an additional 6-hour rainfall observation on the national scale was applied. SWMM was used to calculate the total accumulative overflow in 10-minute units. Data augmentation was used to increase the quantitative amount of input data, and the variation in the predictive performance was examined for each of the input data conditions. The main findings are as follows.
(1)
Flood analysis was performed on the drainage basin of Samseong-dong, Seoul. The total accumulative overflow results for each heavy rainfall event were calculated using 70 observed rainfall events and EPA-SWMM. The characteristics of each rainfall event were analyzed and the correlation with the total accumulative overflow was calculated. As a result, rainfall characteristics that had a high correlation with urban outflow phenomena were identified. Although the highest correlation was found for the total rainfall and rainfall intensity, it was also found that maximum rainfall in 1- to 3-hour units during 6 hours of heavy rainfall was also highly correlated. It was found that the peak rainfall location did not have significant influence during rainfall events and that kurtosis, skewness, and rainfall time for rainfall events had a negative correlation with urban runoff.
(2)
As data augmentation was applied, it was found that the mean absolute error (MAE) and mean absolute percentage error (MAPE) values of the predictive results for testing data decreased for all input data combinations, which verified that prediction performance for data that were not applied to the training can be sufficiently improved using data augmentation. The largest difference of error analysis between the initial input condition (predicted with original DNN, Table 6) and after applying the 10th data augmentation was indicated at the CASE 4-based DNN model. Although input data with low correlation was used in CASE 4, it was judged that data augmentation could be helpful to partially overcome the poor predictive power.
(3)
The total accumulative overflow for the rainfall event on 27 July 2011, was predicted using the DNN constructed according to input data combinations. Predictions were made according to the data augmentation method, and the predictions of CASE 3 using a highly correlated input data combination were the closest to the results of SWMM. In the case of CASE 4, which used the least correlated input data, as more data augmentations were applied, the poorer the predictive results were. Although data augmentation can be used to make up for the lack of input data and reduce errors in learning, it is necessary to conduct proper correlation analysis between the input data and the target value data beforehand. It took 14 minutes for the one-dimensional urban runoff analysis of the SWMM model to be completed; however, predictions using the DNN took 2–3 seconds. In an event of heavy rainfall causing actual flooding, it is likely to save a lot of time in estimating the degree of urban flooding.
(4)
By successfully performing predictions using observed rainfall data and applying data augmentation, basic research on data supplementation techniques in the data-based analysis could be performed. The result of the proposed DNN model is expected to be used as basic data for the real-time flood response in urban areas. If a predictive model is constructed not only for the drainage basin in Samseong-dong but also for all drainage basins in Seoul, it seems that it would be possible to use practically for the entire flood forecasting and warning system in Seoul. Furthermore, if the predicted total accumulative overflow is linked with the expected inundation map, the rapid-simulation of a two-dimensional flood map could be possible. This methodology would be helpful to identify the flood risk area in an urban watershed.

Author Contributions

Conceptualization, H.I.K.; Methodology, H.I.K.; Writing original draft, H.I.K.; Supervision, K.Y.H.; Writing-review & editing, K.Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Korea Environment Industry & Technology Institute(KEITI) though Water Management Research Program, funded by Korea Ministry of Environment(MOE)(79609).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mosavi, A.; Ozturk, P.; Chau, K.W. Flood Prediction Using Machine Learning Models: Literature Review. Water 2018, 10, 1536. [Google Scholar] [CrossRef] [Green Version]
  2. Granata, F.; Gargano, R.; Marinis, G. Support Vector Regression for Rainfall-Runoff Modeling in Urban Drainage: A Comparison with the EPA’s Storm Water Management Model. Water 2016, 8, 69. [Google Scholar] [CrossRef]
  3. Talei, A.; Chua, L.H.C.; Wong, T.S.W. Evaluation of Rainfall and Discharge Inputs used by Adaptive Network-based Fuzzy Inference Systems (ANFIS) in rainfall-runoff modeling. J. Hydrol. 2010, 391, 248–262. [Google Scholar] [CrossRef]
  4. Shen, C. A Transdisciplinary Review of Deep Learning Research and Its Relevance for Water Resources Scientists. Water Resour. Res. 2018, 54, 8558–8593. [Google Scholar] [CrossRef]
  5. Hu, C.; Wu, Q.; Li, H.; Jian, S.; Li, N.; Lou, Z. Deep Learning with a Long Short-Term Memory Networks Approach for Rainfall-Runoff Simulation. Water 2018, 10, 1543. [Google Scholar] [CrossRef] [Green Version]
  6. Li, X.; Willems, P. A Data-Driven Hybrid Urban Flood Modeling Approach. In EPiC Series in Engineering, HIC 2018, Proceedings of the 13th International Conference on Hydroinformatics, Palermo, Italy, 1–6 July 2018; Loggia, G.L., Freni, G., Puleo, V., Marchis, M.D., Eds.; EasyChair: Manchester, UK, 2018; Volume 3, pp. 1193–1200. [Google Scholar]
  7. Nikhil, B.C.; Arjun, N.; Keerthi, C.; Sreerag, S.; Ashwin, H.N. Flood Prediction using Flow and Depth Measurement with Artificial Neural Network in Canals. In Proceedings of the Third International Conference on Computing Methodologies and Communication (ICCMC 2019), Erode, India, 27–29 March 2019; pp. 798–801. [Google Scholar]
  8. Korea Meteorological Agency Meteorological Database. 2020. Available online: https://data.kma.go.kr (accessed on 2 September 2019).
  9. Yoon, S.S.; Bae, D.H.; Choi, Y.J. Urban Inundation Forecasting Using Predicted Radar Rainfall: Case Study. J. Korean Soc. Hazard Mitig. 2014, 14, 117–126. [Google Scholar] [CrossRef] [Green Version]
  10. Huber, W.C.; Dickson, R.E. Storm Water Management Model. User’s Manual Version 4; Environmental Protection Agency: Washington, DC, USA, 1988.
  11. United States Environmental Protection Agency (EPA). Storm Water Management Model User’s Manual Version 5.0; Environmental Protection Agency: Washington, DC, USA, 2010.
  12. Park, J.H.; Kim, S.H.; Bae, D.H. Evaluating Appropriateness of the Design Methodology for Urban Sewer System. J. Korea Water Resour. Assoc. 2019, 52, 411–420. [Google Scholar]
  13. Mark, O.; Weesakul, S.; Apirumanekul, C.; Arronnet, S.B.; Djordjevic, S. Potential and Limitations of 1D Modeling of Urban Flooding. J. Hydrol. 2004, 299, 284–299. [Google Scholar] [CrossRef]
  14. Izumi, T.; Miyoshi, M.; Kobayashi, N. Runoff Analysis Using a Deep Neural Network. In Proceedings of the 12th International Conference on Hydroscience & Engineering, Hydro-Science & Engineering for Environmental Resilience, Taiwan, China, 6–10 November 2016. [Google Scholar]
  15. Agarap, A.F. Deep Learning using Rectified Linear Units (ReLU). arXiv 2019, arXiv:1803.08375v2. [Google Scholar]
  16. Kingma, D.P.; Ba, J.L. ADAM: A Method for Stochastic Optimization. arXiv 2015, arXiv:1412.6980. [Google Scholar]
  17. Cui, X.; Goel, V.; Kingsbury, B. Data Augmentation for Deep Neural Network Acoustic Modeling. In Proceedings of the 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 4–9 May 2014; pp. 1469–1477. [Google Scholar]
  18. Xu, Y.; Jia, R.; Mou, L.; Li, G.; Chen, Y.; Lu, Y.; Jin, Z. Improved Relation Classification by Deep Recurrent Neural Networks with Data Augmentation. arXiv 2016, arXiv:1601.03651v2. [Google Scholar]
  19. Seoul Metropolitan City. Comprehensive Plan for Storm and Flood Damage Reduction; Seoul Metropolitan City: Seoul, Korea, 2015; Volume 1, pp. 374–375. [Google Scholar]
  20. Risi, R.D.; Jalayer, F.; Paola, F.D. Meso-scale hazard zoning of potentially flood prone areas. J. Hydrol. 2015, 527, 316–325. [Google Scholar] [CrossRef]
  21. Moore, M.R. Development of a High-Resolution 1D/2D Coupled Flood Simulation of Charles City, Iowa. Master’s Thesis, University of Iowa, Iowa City, IA, USA, 2011. [Google Scholar]
  22. Son, A.L.; Kim, B.H.; Han, K.Y. A study on prediction of inundation area considering road network in urban area. J. Korean Soc. Civ. Eng. 2015, 35, 307–318. [Google Scholar] [CrossRef] [Green Version]
  23. Korea Institute of Civil Engineering and Building Technology. Road Design Manual; KICT: Korea, Seoul, 2001. [Google Scholar]
Figure 1. (a) Deep neural network (DNN) and (b) rectified linear units (ReLU) activation function.
Figure 1. (a) Deep neural network (DNN) and (b) rectified linear units (ReLU) activation function.
Water 12 00899 g001
Figure 2. Flowchart for data augmentation and training the DNN.
Figure 2. Flowchart for data augmentation and training the DNN.
Water 12 00899 g002
Figure 3. Study area and nearby Automatic Weather System (AWS) stations. (a) shows the drainage system, and (b) shows rainfall observatory.
Figure 3. Study area and nearby Automatic Weather System (AWS) stations. (a) shows the drainage system, and (b) shows rainfall observatory.
Water 12 00899 g003
Figure 4. Flood hazard zone in Samsung-dong. (a) shows Samsung–Seonreung Station, and (b) shows Bongeunsa–Samsung Station.
Figure 4. Flood hazard zone in Samsung-dong. (a) shows Samsung–Seonreung Station, and (b) shows Bongeunsa–Samsung Station.
Water 12 00899 g004
Figure 5. Verification of 2D flood analysis results. (a) shows flood mark (2010), (b) shows simulation result (21 September 2010), (c) shows flood mark (2011), and (d) shows simulation result (27 July 2011).
Figure 5. Verification of 2D flood analysis results. (a) shows flood mark (2010), (b) shows simulation result (21 September 2010), (c) shows flood mark (2011), and (d) shows simulation result (27 July 2011).
Water 12 00899 g005
Figure 6. The trend of mean absolute error (MAE) with validation and testing data. (a) Shows Case 1: original data, (b) shows Case 1: after the 10th data augmentation, (c) shows Case 2: original data, (d) shows Case 2: after the 10th data augmentation, (e) shows Case 3: original data (f) shows Case 3: after the 10th data augmentation, (g) shows Case 4: original data and (h) shows Case 4: after the 10th data augmentation.
Figure 6. The trend of mean absolute error (MAE) with validation and testing data. (a) Shows Case 1: original data, (b) shows Case 1: after the 10th data augmentation, (c) shows Case 2: original data, (d) shows Case 2: after the 10th data augmentation, (e) shows Case 3: original data (f) shows Case 3: after the 10th data augmentation, (g) shows Case 4: original data and (h) shows Case 4: after the 10th data augmentation.
Water 12 00899 g006aWater 12 00899 g006b
Figure 7. Observed rainfall on 27 July 2011.
Figure 7. Observed rainfall on 27 July 2011.
Water 12 00899 g007
Table 1. Survey of observed rainfall nationwide.
Table 1. Survey of observed rainfall nationwide.
Rainfall StationDays of Observation Total   Rainfall   ( m m )
Seoul AreaSeoul (108) ASOS21 September 2010,
27 July 2011
267.5, 116.5
Gangnam (400) AWS21 September 2010,
27 July 2011,
15 August 2012,
22 July 2013
253.5, 184.5,
131.5, 140.5
Kwanak Mountain (116) Weather Radar21 September 2010,
27 July 2011
110.0, 184.5
Seocho (401) AWS21 September 2010,
27 July 2011,
15 August 2012,
22 July 2013
258.0, 201.0,
129.5, 128.0
Gangdong (402) AWS21 September 2010,
27 July 2011
275.5, 111.5
Songpa (403) AWS21 September 2010,
27 July 2011
252.0, 194.5
Gangseo (404) AWS21 September 2010,
27 July 2011
70.0, 178.5
Yangchoen (405) AWS21 September 2010,
27 July 2011
52.0, 180.5
Dobong (406) AWS21 September 2010166.5
Nowon (407) AWS21 September 2010,
27 July 2011
57.5, 88.5
Dongdaemun (408) AWS21 September 2010,
27 July 2011
232.5, 127.5
Jungrang (409) AWS21 September 2010,
27 July 2011
267.5, 112.0
KMA (410) AWS21 September 2010,
27 July 2011
251.5, 164.5
Mapo (411) AWS21 September 2010,
27 July 2011
241.5, 177.5
Seodaemun (412) AWS21 September 2010,
27 July 2011
154.5, 121.5
Gwangjin (413) AWS21 September 2010,
27 July 2011
249.0, 117.0
Seongbuk (414) AWS21 September 2010,
27 July 2011
160.0, 106.5
Yongsan (415) AWS21 September 2010,
27 July 2011
180.0, 136.5
Eunpyeong (416) AWS21 September 2010,
27 July 2011
238.5, 25.5
Geumcheon (417) AWS21 September 2010,
27 July 2011
160.0, 47.0
Hangang (418) AWS21 September 2010,
27 July 2011
240.5, 150.0
Seongdong (421) AWS21 September 2010,
27 July 2011
216.5, 118.0
Bukak Mountain (422) AWS27 July 2011125.0
Gangbuk (424) AWS21 September 2010,
27 July 2011
119.5, 81.0
Namhyeon (425) AWS27 July 2011248.0
Kwanak (509) AWS21 September 2010,
27 July 2011
91.0, 245.0
Yeongdeungpo (510) AWS21 September 2010,
27 July 2011
260.0, 181.0
Gwacheon (590) AWS21 September 2010,
27 July 2011
63.5, 140.0
Other AreasIncheon (112) ASOS23 July 201769.3
Cheonan (232) ASOS16 July 2017223.6
Chungju (127) ASOS16 July 2017284.4
Busan (159) ASOS16 July 2009,
27 July 2011,
15 July 2012,
25 August 2014,
11 September 2017
243.0, 229.5,
114.0, 91.0,
254.1,
Ulsan (152) ASOS5 October 2016,
11 September 2017
91.2
Changwon (155) ASOS25 August 2014,
5 October 2016
209.5, 116.7
Table 2. Observed rainfall and statistical features (four observation stations in the Seoul area).
Table 2. Observed rainfall and statistical features (four observation stations in the Seoul area).
Statistical CharacteristicGangnam StationGangnam StationDobong StationGangbuk Station
Observation Date21 September 201027 July 201121 September 201027 July 2011
Total Rainfall (mm)253.50184.50166.5081.00
Max. Rainfall in 1 h (mm)78.0071.0057.0025.00
Max. Rainfall in 2 h (mm)136.50105.50106.0038.00
Max. Rainfall in 3 h (mm)181.50143.00136.0041.00
Rainfall Intensity (mm/h)42.2530.7527.7513.35
Peak Rainfall Location (%)25.0055.5666.6777.78
Standard Deviation (mm)5.174.384.682.49
Skewness0.391.280.401.05
Inter-event Time (min)6010100130
Kurtosis−0.49731.2289−1.47150.3522
Table 3. Observed rainfall and statistical features (four observation stations in other areas).
Table 3. Observed rainfall and statistical features (four observation stations in other areas).
Statistical CharacteristicBusan StationUlsan StationChangwon StationCheonan Station
Observation Date11 September 20175 October 20165 October 201616 July 2017
Total Rainfall (mm)254.10233.80116.70223.60
Max. Rainfall in 1 h (mm)84.90103.1074.9069.30
Max. Rainfall in 2 h (mm)144.90166.4097.10118.50
Max. Rainfall in 3 h (mm)176.20203.30107.20142.50
Rainfall Intensity (mm/h)42.3538.9719.4537.27
Peak Rainfall Location (%)36.1177.7875.0075.00
Standard Deviation (mm)5.156.064.964.05
Skewness1.050.892.250.33
Inter-event Time (min)00200
Kurtosis0.6922−0.45084.7215−0.8722
Table 4. Combination of input data.
Table 4. Combination of input data.
CASE 1
(9 Inputs)
CASE 2
(6 Inputs)
CASE 3
(5 Inputs)
CASE 4
(4 Inputs)
R 2   with
Total Accumulative Overflow
Total RainfallTotal RainfallTotal RainfallTotal Rainfall0.8144
Max. Rainfall in 1 h Max. Rainfall in 1 h Max. Rainfall in 1 h -0.5698
Max. Rainfall in 2 h Max. Rainfall in 2 h Max. Rainfall in 2 h -0.5114
Max. Rainfall in 3 h Max. Rainfall in 3 h Max. Rainfall in 3 h -0.3957
Rainfall Intensity Rainfall Intensity Rainfall Intensity -0.566
Standard DeviationStandard Deviation--0.3623
Skewness--Skewness(−) 0.3778
Kurtosis--Kurtosis(−) 0.2499
Inter-event Time--Inter-event Time(−) 0.1087
Table 5. Empirical process for finding the DNN structure.
Table 5. Empirical process for finding the DNN structure.
No.MAE (Testing)Number of BatchesStructure
119.6344/9/18/12/10
221.0244/9/18/18/12/12/10
323.244/9/18/18/12/10
420.8444/18/18/12/12/12/10/10
519.8188/18/18/12/12/12/10/10
619.6488/18/18/12/12/12/10/10
7(Selected)17.51818/18/12/12/12/10/10/8
818.281018/18/12/12/12/10/10/8
919.05818/18/18/12/12/12/10/10/8
Table 6. Training and validation error with data augmentation.
Table 6. Training and validation error with data augmentation.
Training and Validation ResultInput Criterion
CASE 1 (MAE)CASE 2 (MAE)CASE 3 (MAE)CASE 4 (MAE)
Train.Val.Train.ValTrain.ValTrain.Val
Original DNN 17.1325.2410.7914.808.4211.2825.1818.32
Number of Data Augmentations Applied110.4015.044.0110.276.689.839.8021.71
27.4511.175.097.693.815.6510.3713.52
35.886.815.435.293.844.708.6311.79
46.697.254.305.333.984.487.9811.25
55.335.014.914.233.823.257.246.08
65.064.713.163.283.113.087.407.77
75.17.743.064.623.024.886.438.83
84.923.943.242.422.952.707.156.04
95.126.152.823.093.62.945.605.80
104.174.713.433.743.313.976.045.78
Table 7. Testing error with data augmentation.
Table 7. Testing error with data augmentation.
Testing ResultInput Criterion
CASE 1CASE 2CASE 3CASE 4
MAEMAPEMAEMAPEMAEMAPEMAEMAPE
Original DNN 23.8761.7513.0237.8311.1136.8923.7749.30
Number of Data Augmentations Applied113.9146.947.9248.755.9021.3127.3271.77
210.0611.045.416.184.756.4612.7514.28
310.69.487.127.135.995.4614.7414.72
48.3312.624.206.733.926.7611.8022.24
59.1716.294.4418.995.4635.138.4027.83
66.2357.273.396.966.2114.347.9612.97
75.1315.914.1211.192.778.076.548.93
81.9611.563.7727.993.777.247.2742.32
96.1413.563.0721.743.1522.836.4432.0
104.4634.423.3811.494.227.834.157.02
Table 8. Results of the simulation and prediction.
Table 8. Results of the simulation and prediction.
Total   Accumulative   Overflow   ( m 3 / s ) Input Criterion
CASE 1CASE 2CASE 3CASE 4
Average of Absolute R-square
(refer to R 2 in Table 4)
0.43950.53660.57150.3877
Simulated with SWMM144.7
Prediction with Original DNN168.66159.10142.15149.90
Predict with Data Augmentation1174.92147.67143.36145.61 *
2149.27160.80152.75140.37
3130.88166.44149.01134.75
4146.63 *148.31150.18132.25
5152.48147.14146.77119.99
6154.49148.82160.14142.34
7162.85143.57 *147.71136.49
8129.85160.21160.21136.32
9151.82149.65162.69137.04
10147.39142.69144.51 *138.47
Standard Deviation13.217.586.326.47

Share and Cite

MDPI and ACS Style

Kim, H.I.; Han, K.Y. Urban Flood Prediction Using Deep Neural Network with Data Augmentation. Water 2020, 12, 899. https://doi.org/10.3390/w12030899

AMA Style

Kim HI, Han KY. Urban Flood Prediction Using Deep Neural Network with Data Augmentation. Water. 2020; 12(3):899. https://doi.org/10.3390/w12030899

Chicago/Turabian Style

Kim, Hyun Il, and Kun Yeun Han. 2020. "Urban Flood Prediction Using Deep Neural Network with Data Augmentation" Water 12, no. 3: 899. https://doi.org/10.3390/w12030899

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop