Design and Operation of Decentralized Reservoirs in Urban Drainage Systems

Poor drainage of urban storm water can lead to urban inundation which presents a risk to people and property. Previous research has presented various measures to prevent and reduce urban flooding and these measures can be classified into costly but effective structural measures, and economical but less effective non-structural measures. This study suggests a new approach to reduce urban flooding by combining structural and non-structural measures in a target watershed in Seoul, South Korea. Inlet design modification in a detention reservoir (Decentralized Reservoir, DR) is examined in conjunction with combined inlet/outlet management for the DR. Monitoring nodes used to control DR inlet/outlet operations are selected by locating the first flooding node, maximum flooding node and DR inlet node. This new approach demonstrates outstanding flood volume reduction for historical flooding events that occurred in Seoul during 2010 and 2011. Flood volumes during the 2010 event using the combined inlet/outlet operation in the DR were between 1656 m 3 and 1815 m 3 compared to a flood volume of 6617 m 3 using current DR operation. Finally, the suggested operating level for the DR based on the best hydraulic section, system resilience index, and local regulations is 1.2 m.


Introduction
The frequency of extreme rainfall events has increased with recent environmental change.Furthermore, the extent of impervious surfaces in urban areas has risen rapidly.Consequently, the incidence of urban inundation has increased markedly.To mitigate the risk of urban inundation, many solutions for preventing urban flooding have been suggested [1].These solutions can be categorized into structural and non-structural measures.Structural measures refer to infrastructure changes including the construction of new drainage facilities and rehabilitation or replacement of sewer networks.Non-structural measures refer to technological or behavioral innovations which do not require significant infrastructure investment, including flood prediction, improved rainfall forecasting, and revised operating procedures for drainage facilities.Non-structural measures are significantly cheaper to apply as compared to structural measures however, they are also somewhat limited in their ability to prevent or reduce flooding.For example, improving the operation of drainage facilities can do little to reduce flooding if the facility was designed to withstand a 30-year maximum rainfall event and a 100-year maximum rainfall event occurs.
Table 1.Classification of the measures employed by previous studies to reduce inundation in urban drainage systems.reducing urban inundation.This study focus on combined measure for detention reservoir in urban drainage system such as alternative inlet design and inlet/outlet operation of detention reservoir.The operation of detention reservoirs as a non-structural measure has a structural limitation such as small capacity of drainage pumps because the design of drainage pumps in detention reservoir were not designed considering inflow to detention reservoir in Korea.For overcoming this limitation, inlet and outlet operation in detention reservoir are required to obtain additional space in detention reservoir.A detention reservoir has no ability to mitigate urban inundation if it is filled with received water.Hence, this study examines a new operation method for inlet and outlet controls in a detention reservoir (Decentralized Reservoir, DR) combined with a new inlet design for the DR.This study consists of four parts: (1) selection of the monitoring nodes used to operate the inlets and outlets in the DR; (2) an investigation of inlet design and replacement of the inlet weir with an inlet gate in the DR; (3) an investigation of the impacts of the new inlet/outlet operating procedures in the DR; and (4) a determination of the operating level in the DR using a system resilience index and the best hydraulic section for urban drainage systems.

Monitoring Node Selection
Careful selection of monitoring nodes is essential to effectively operate inlet gates and outlet drainage pumps in the DR.Almost all detention reservoirs in the urban drainage systems of South Korea have inlets that are designed as weirs.Consequently, managed inlet operation is impossible because the inlet weir has a fixed height and all flow over the weir enters the detention reservoir.Managed inlet operation is made possible by replacing inlet weirs with gates.The effects of managed inlet gate operations, which are controlled by the level of the monitoring nodes, are compared to the current inlet weirs in DR.
The DR outlet is operated based on monitoring nodes within the urban drainage system and the effects of managed outlet operation were verified in an earlier study [1].The inlet gate is open and outlet drainage pumps are not operated when water levels at the monitoring nodes are high.Conversely, the inlet gate is closed and outlet drainage pumps are operated when water levels at the monitoring nodes are low.The effects of combined operation, when the inlet gate and outlet drainage pumps are both operated according to the level of the monitoring nodes, is compared with the current operation in DR.The operating levels in the DR are determined by the resilience index and best hydraulic section in the urban drainage system.
A flowchart for this study is shown in Figure 1.The first part of this study is the selection of monitoring nodes for inlet/outlet operations in DR.Two methods were used to select the monitoring nodes.Firstly, the first flooding node was detected using a rainfall-runoff simulation because the first flooding node is one of the most susceptible nodes for inundation in an urban drainage system.To maintain a safe urban drainage system, identification of the first flooding node is critical for effective DR operation.Secondly, the maximum flooding node was identified.This node is important because it illustrates the most at risk region of the urban drainage system.inlet gate operations, which are controlled by the level of the monitoring nodes, are compared to the current inlet weirs in DR.
The DR outlet is operated based on monitoring nodes within the urban drainage system and the effects of managed outlet operation were verified in an earlier study [1].The inlet gate is open and outlet drainage pumps are not operated when water levels at the monitoring nodes are high.Conversely, the inlet gate is closed and outlet drainage pumps are operated when water levels at the monitoring nodes are low.The effects of combined operation, when the inlet gate and outlet drainage pumps are both operated according to the level of the monitoring nodes, is compared with the current operation in DR.The operating levels in the DR are determined by the resilience index and best hydraulic section in the urban drainage system.
A flowchart for this study is shown in Figure 1.The first part of this study is the selection of monitoring nodes for inlet/outlet operations in DR.Two methods were used to select the monitoring nodes.Firstly, the first flooding node was detected using a rainfall-runoff simulation because the first flooding node is one of the most susceptible nodes for inundation in an urban drainage system.To maintain a safe urban drainage system, identification of the first flooding node is critical for effective DR operation.Secondly, the maximum flooding node was identified.This node is important because it illustrates the most at risk region of the urban drainage system.Thorough evaluation of monitoring nodes is required because the design frequency of conduits in urban drainage systems frequently differ.For example, the design frequency of branch conduits tends to be smaller than that of main conduits which means that flooding can easily occur in nodes Thorough evaluation of monitoring nodes is required because the design frequency of conduits in urban drainage systems frequently differ.For example, the design frequency of branch conduits tends to be smaller than that of main conduits which means that flooding can easily occur in nodes between branch conduits.However, these nodes should not be selected for monitoring purposes because they do not represent the entire drainage system.Hence, rainfall-runoff simulations should be conducted to facilitate the selection of monitoring nodes.Synthetic rainfall data by Huff (1967) [36] were used in this study because the DR in the target watershed was designed using the Huff distribution.
Selection of first flooding nodes starts with the determination of rainfall distribution.Then, several durations are selected considering the time of concentration in target watershed.The initial value of total rainfall volume is 1 mm and the total rainfall is distributed by selected distribution.The distributed rainfall is applied to rainfall-runoff model and flooding results are obtained.If flooding occurs, the flooding node is selected as the first flooding node and the flooding volume is also checked.If no flooding occurs, the total rainfall volume is assumed 2 mm with 1 mm increment.This process is repeated until first flooding nodes are obtained.
Selection of maximum flooding nodes also begins with the determination of rainfall distribution.Various frequencies and durations are selected and generated by selected distribution.The distributed rainfall is applied to rainfall-runoff model and flooding results are obtained.The maximum flooding node among flooding nodes is examined in each rainfall data, respectively.If it is hard to select the monitoring node because maximum flooding nodes is different from each other, historical rainfall data are applied to rainfall-runoff model for obtaining the maximum flooding node.

Generation of Synthetic Rainfall Data
Synthetic rainfall data are necessary for monitoring node selection and to determine the operating level of the DR.The design of almost all hydraulic facilities in South Korea is based on the Huff distribution [36].This distribution consists of four quartiles and the difference between each quartile represents the timing of peak rainfall.For example, the first quartile has the peak value of rainfall in the first quarter of rainfall duration.The regression equations for the first, second, third, and fourth quartiles of the Huff distribution for Seoul, South Korea are shown in Equations ( 1)-(4), respectively (Korea Precipitation Frequency Data Server) [37], y = 29.289x 6 where y represents the cumulative rainfall ratio and x is the cumulative time ratio in total duration.In Seoul, the third-quartile of the Huff distribution was determined to be an appropriate rainfall distribution based on an examination of rainfall events and calculating discharge in the urban watershed (Yoon et al., 2013) [38].

Alteration of Inlet Type in the DR
In South Korea, almost all DR in urban drainage systems have an inlet weir and the design frequency of the weir is generally smaller than the design frequency of the DR itself.The present operation in the DR is based on outlet management using drainage pumps.Reserved water in the DR is released when one of two conditions is satisfied; firstly, water reaches the limit level of the DR, and secondly, operators of the DR determine that the rainfall event has concluded.The capacity of the drainage pumps in the DR is determined at the design phase of the reservoir using the designed drainage time of the DR.The general structure of the DR is shown in Figure 2. Inlet management of the DR using a fixed inlet weir is impossible and the storage capacity of the DR decreases during continuous rainfall.The DR cannot reserve capacity for follow-up inflow because a DR with a fixed inlet weir receives initial low frequency inflow during continuous rainfall events.To address this issue, the inlet weir is replaced with an inlet gate, which is operated based on the water level at the monitoring and inlet nodes.
Water 2017, 9, 246 5 of 14 the inlet weir is replaced with an inlet gate, which is operated based on the water level at the monitoring and inlet nodes.

Determination of DR Operational Level
The resilience index of an urban drainage system can be defined as the ability of the system to recover from a failure such as the malfunction of drainage facilities or inundation [1].Resilience index is defined using Equation ( 5), where u(T)i is the performance evaluation function at time i, Fi is the flooding volume (m 3 ) at time i, tc is time of concentration (min) in the target watershed [39], Ri is the rainfall at the present time (i), Rd is accumulated rainfall during the present time (i) and i-tc (10 −3 m), and Au is the unit area (1 ha = 10 −4 m 2 = 0.01 km 2 ).The time of concentration is a concept to be used in hydrology and it is the estimated time while rainfall in the farthest region of watershed flows to watershed outlet.Equation ( 5) cannot be calculated if flooding occurs when there is no rainfall because the denominator of Equation ( 5) is 0 if there is no rainfall.The negative value of Equation ( 5) does not exist because the value of rainfall, flooding volume, and area is not the negative value.The unit area instead of the catchment area is used because the big value of the catmint area could enormously reduce failure depth.A high performance evaluation value indicates that the drainage system has a high ability to recover after flooding or system failure.The resilience of the urban drainage system by performance evaluation function per minute is calculated as shown in Equation ( 6), where Rs is the resilience of the urban drainage system in target watershed, T0 is the start time of the rainfall event and Tn is the end time of the rainfall event.In this study, 48 rainfall datasets with four frequencies (30, 50, 80, and 100 years) and three durations (30, 60, and 90 min) from the Huff Distribution [36] are used.The operational level of the DR is calculated at the point where the system begins to flood.

Determination of DR Operational Level
The resilience index of an urban drainage system can be defined as the ability of the system to recover from a failure such as the malfunction of drainage facilities or inundation [1].Resilience index is defined using Equation ( 5), where u(T) i is the performance evaluation function at time i, F i is the flooding volume (m 3 ) at time i, t c is time of concentration (min) in the target watershed [39], R i is the rainfall at the present time (i), R d is accumulated rainfall during the present time (i) and i-tc (10 −3 m), and A u is the unit area (1 ha = 10 −4 m 2 = 0.01 km 2 ).The time of concentration is a concept to be used in hydrology and it is the estimated time while rainfall in the farthest region of watershed flows to watershed outlet.Equation ( 5) cannot be calculated if flooding occurs when there is no rainfall because the denominator of Equation ( 5) is 0 if there is no rainfall.The negative value of Equation ( 5) does not exist because the value of rainfall, flooding volume, and area is not the negative value.The unit area instead of the catchment area is used because the big value of the catmint area could enormously reduce failure depth.A high performance evaluation value indicates that the drainage system has a high ability to recover after flooding or system failure.The resilience of the urban drainage system by performance evaluation function per minute is calculated as shown in Equation (6), where R s is the resilience of the urban drainage system in target watershed, T 0 is the start time of the rainfall event and T n is the end time of the rainfall event.In this study, 48 rainfall datasets with four frequencies (30, 50, 80, and 100 years) and three durations (30, 60, and 90 min) from the Huff Distribution [36] are used.The operational level of the DR is calculated at the point where the system begins to flood.

Target Area and Rainfall Runoff Model
Seoul, South Korea is traversed by the Han River (Figure 3).In the target watershed, flooding occurred twice, in 2010 and 2011, after the DR was completed [40,41].The selected rainfall-runoff model is the Storm Water Management Model (SWMM), developed by the US Environmental Protection Agency [42].The Seoul Metropolitan Government supplied urban drainage system information and Geographic Information System (GIS) data for Seoul.The sewer network for the Daerim 3 pump station is shown in Figure 4.In the target watershed, flooding occurred twice, in 2010 and 2011, after the DR was completed [40,41].The selected rainfall-runoff model is the Storm Water Management Model (SWMM), developed by the US Environmental Protection Agency [42].The Seoul Metropolitan Government supplied urban drainage system information and Geographic Information System (GIS) data for Seoul.The sewer network for the Daerim 3 pump station is shown in Figure 4.

Application Results and Discussion
To select monitoring nodes for the DR, we examined the first flooding nodes in the target watershed.Considering the time of concentration in the target watershed (30 min), three durations of 30, 60, and 90 min were selected for the rainfall-runoff simulation.To be a candidate for monitoring, nodes needed to be located between the main conduits and have a subcatchment area greater than 0.12 km 2 .First flooding nodes located in main conduits of the target watershed are shown in Figure 4.The first flooding node for 30-and 60-min durations is node 560, when the total rainfall is 79 mm and 104 mm, respectively.The first flooding node for the 90-min duration simulation is node 575, with total rainfall of 154 mm.Considering the time of concentration in the target watershed, node 560 was selected as a monitoring node.
To identify the maximum flooding nodes, nine rainfall datasets (frequency: 30, 50 and 100 years; duration: 30, 60 and 90 min) were generated using the Huff distribution [36].However, these results were inconclusive, with nine different maximum flooding nodes identified by the nine simulations.Consequently, historical rainfall data, from 23 September 2010 and 27 July 2011, when flooding occurred in the watershed, were used to identify the maximum flooding nodes.Based on the historical data, node 550 was found to be the maximum flooding node and was selected as a monitoring node.To summarize, nodes 550 and 560, and the inlet node of the Daerim DR were selected as monitoring nodes.
The depth of the monitoring nodes is 1.5 m.For the combined inlet/outlet operation of the detention reservoir, operating levels between 0.8 m (approximately half of node depth) and 1.5 m (node depth) were selected because the simulated operating levels from 0.0 m to 0.7 m showed similar or worse results as compared to the current operation of the DR.Low operating levels mean that detention reservoir receives water and reserved water (not discharge) when water levels in the sewer network are low.
The depth of the monitoring nodes means the pipe diameter.Small diameter is chosen as the

Application Results and Discussion
To select monitoring nodes for the DR, we examined the first flooding nodes in the target watershed.Considering the time of concentration in the target watershed (30 min), three durations of 30, 60, and 90 min were selected for the rainfall-runoff simulation.To be a candidate for monitoring, nodes needed to be located between the main conduits and have a subcatchment area greater than 0.12 km 2 .First flooding nodes located in main conduits of the target watershed are shown in Figure 4.The first flooding node for 30-and 60-min durations is node 560, when the total rainfall is 79 mm and 104 mm, respectively.The first flooding node for the 90-min duration simulation is node 575, with total rainfall of 154 mm.Considering the time of concentration in the target watershed, node 560 was selected as a monitoring node.
To identify the maximum flooding nodes, nine rainfall datasets (frequency: 30, 50 and 100 years; duration: 30, 60 and 90 min) were generated using the Huff distribution [36].However, these results were inconclusive, with nine different maximum flooding nodes identified by the nine simulations.Consequently, historical rainfall data, from 23 September 2010 and 27 July 2011, when flooding occurred in the watershed, were used to identify the maximum flooding nodes.Based on the historical data, node 550 was found to be the maximum flooding node and was selected as a monitoring node.To summarize, nodes 550 and 560, and the inlet node of the Daerim DR were selected as monitoring nodes.
The depth of the monitoring nodes is 1.5 m.For the combined inlet/outlet operation of the detention reservoir, operating levels between 0.8 m (approximately half of node depth) and 1.5 m (node depth) were selected because the simulated operating levels from 0.0 m to 0.7 m showed similar or worse results as compared to the current operation of the DR.Low operating levels mean that detention reservoir receives water and reserved water (not discharge) when water levels in the sewer network are low.
The depth of the monitoring nodes means the pipe diameter.Small diameter is chosen as the depth of the monitoring nodes if the diameter of the inlet pipe is different from it of the outlet pipe.If the depth of the monitoring nodes is different from each other, each depth is applied to the range from 0.0 D to 1.0 D. For example, 1.5 m is same with 1.0 D if the depth of the monitoring node is 1.5 m and 1.2 m is same with 1.0 D if the depth of the monitoring node is 1.2 m.
The inlet weir was replaced with and inlet gate in the DR and the impact of combined inlet/outlet operation with operating levels from 0.8 m to 1.5 m were assessed.Figure 5 shows the comparison between the current operation, outlet operation, and combined inlet/outlet operation in the DR. Figure 5a shows the results of the three operations in the 2010 flooding event and Figure 5b shows the results of the three operations in the 2011 event.The inlet weir was replaced with and inlet gate in the DR and the impact of combined inlet/outlet operation with operating levels from 0.8 m to 1.5 m were assessed.Figure 5 shows the comparison between the current operation, outlet operation, and combined inlet/outlet operation in the DR. Figure 5a shows the results of the three operations in the 2010 flooding event and Figure 5b shows the results of the three operations in the 2011 event.For the 2010 event, the result of current operation is 6617 m 3 regardless of the level of monitoring nodes because current operation has no relation with the level of monitoring nodes.The results of outlet operation in the DR are distributed from 4466 m 3 to 6617 m 3 .There are no differences between the current operation and outlet operation in the DR when the water level of the monitoring nodes is low, that is, when drainage pumps in the DR were stopped.The flooding volume using combined inlet/outlet operation in the DR range between 1656 and 1815 m 3 .The results of combined inlet/outlet operation in the DR show the best performance of the three operating methods for all levels of monitoring nodes.
For the 2011 event, flood volume using current operation is 664 m 3 for all monitoring node levels.Results for outlet operation in the DR are between 552 m 3 and 703 m 3 .Poor results are apparent for low water levels at monitoring nodes between 0.8 m and 0.9 m. Results are similar to the current operation when the water levels at monitoring nodes are between 1.0 m and 1.1 m.When the monitoring node water levels are high (1.2 m, 1.3 m, 1.4 m, and 1.5 m), outlet operation shows better results as compared to current operation.The flooding volume for inlet/outlet operation in the DR is distributed from 551 m 3 to 556 m 3 .The results of combined inlet/outlet operation in the DR show the best performance of the three operations.It is notable that combined inlet/outlet operation results are superior for all monitoring node levels.Combined inlet/outlet operation shows the smallest variance of flooding volume and total flooding volume for both the 2010 and 2011 events.These results imply that that combined inlet/outlet operations of the DR can substantially reduce flooding volume.
Outlet operation performs better than current operation except when the water level at the monitoring nodes is low.Outlet operation in the DR is limited by the capacity of the drainage pumps (0.3 m 3 /s).The simulated maximum inflow to the DR was approximately 1.33 m 3 /s when rainfall frequency distribution was 100 years.This means that operators cannot empty or fill the DR.Inlet/outlet operation shows outstanding performance compared to both current operation and outlet operation.The results of system resilience based on Huff distribution synthetic rainfall data are shown in Table 2. Values of system resilience in Table 2 are calculated by values of the performance evaluation function ranging from 0 to 1.The system resilience of urban drainage system is 1 if there are no flood.System resilience is high in low frequency rainfall events such as 30-year and it is low in high frequency rainfall events such as 100-year.Values of system resilience in Table 2 are different from each other according to levels of monitoring nodes.For example, an average system resilience in all frequency and duration is 0.946162 when levels of monitoring nodes are 0.8 m and it is 0.949053 when levels of monitoring nodes are 1.5 m.System resilience in high levels of monitoring nodes from 1.2 m to 1.5 m is stably high.
Determining the operating level in the DR is difficult because it is directly linked with urban inundation in the target watershed.In this section, the operating level in the DR was determined using a resilience index and the concept of best hydraulic section.For calculating the system resilience in the target watershed, the flooding volume is required.Performance evaluation of the target watershed using Equation ( 5) is calculated using the flooding volume, and system resilience is calculated using Equation (6).A total of 12 rainfall events with four frequencies (30, 50, 80, and 100 years) and three durations (60, 120 and 180 min) were generated using the Huff distribution for rainfall-runoff simulation.Frequencies of rainfall scenarios were all over 30 years because there was no inundation when rainfall frequency distributions under 30 years were applied to the rainfall-runoff simulations.Furthermore, durations of rainfall data are all more than one hour because there is no difference between DR operating methods when rainfall durations are short.The results of system resilience are shown in Table 2.
The average of total system resilience in Table 2 was 0.948891 and the average system resilience for each operating level was considered to determine the operating level in the DR.Average system resilience for each operating level and the total average system resilience are shown in Figure 6.The average system resilience of operating levels from 1.2 m to 1.5 m is higher than the total average system resilience.
The ratio between partial flow Q and full flow QF is shown in Equation (10).
Equation ( 8) is differentiated about θ to obtain the maximum discharge (Qmax) of a conduit.Qmax is calculated when θ is 5.28 radian (θ = 303°).The ratio of d and D is 0.938 when θ is 303°.Hence, the depth of maximum discharge in the monitoring nodes (d) is 1.407 m because the total depth of monitoring nodes (D) is 1.5 m.The range of the ratio of d and D when Q is larger than Qmax is from 0.81 to 1.00 and the range of monitoring nodes is approximately from 1.2 m to 1.5 m.The level of conduits also gratifies the qualification of best hydraulic section if the level of the monitoring nodes satisfies the condition of best hydraulic section.The best hydraulic section in conduits is considered.In closed conduits, water flows with gravity if it has a free water surface.The best hydraulic section in rectangular conduits cannot be calculated though the best hydraulic section in circular conduits because the velocity and flow rate in rectangular conduits linearly increase.However, the best hydraulic section in circular conduits is suggested for determination of operating level in detention reservoir because more than 80% of all conduits are circular conduits in target watershed and the majority of conduits in Korea are circular conduits.Circular conduits in urban drainage systems have different relation between discharge and velocity when whole conduit sections are filled with water compared to when partial section of conduits are filled with water.This relation can be described as a hydraulic characteristic curve.Figure 7 shows the cross section of circular conduit.Discharge area, wetted perimeter and hydraulic mean radius in conduits is calculated by Equation ( 7), where A is the discharge area in circular conduit (m 2 ), D is the diameter of conduit (m), θ is the angle between the center of conduit and water surface (radian), P is the wetted perimeter (m), and R h is the hydraulic mean radius (m).Discharge in the conduit can be described by Manning's equation, Equation (8), where Q is the discharge (m 3 /s), n is Manning's roughness coefficient, S 0 is the slope of conduit, and K is the constant ratio between S 0 and n. θ is determined from 0 to 2π.If the conduit is filled with water, then the discharge of conduit can be explained as Q F in Equation ( 9).
Water 2017, 9, 246 11 of 14 Engineers should consider a safety factor including differences between real discharge and design discharge when they design urban drainage conduits in South Korea.Safety factors for urban drainage conduits are based on the diameter of the conduit.For example, regulations stipulate that large conduits with diameters over 1500 mm should have an excess capacity of between 25% and 50% (Korea Water and Wastewater Association, 2011) [44].Safety factors in monitoring nodes range from 1.25 to 1.5.The operating level including a safety factor is from 1.0 m (safety factor = 1.5 m/1.5 = 1.0 m) to 1.2 m (safety factor = 1.5 m/1.25 = 1.2 m).The possible operating level in the DR based on the resilience index, best hydraulic section, and local regulation is shown in Table 3.

Conclusions
This study suggested combined changes to inlet design and inlet/outlet operation of DR in urban drainage systems to reduce the incidence of urban inundation.This study was comprised of four steps.Firstly, monitoring nodes for operating DR inlet and outlet were selected.Secondly, the DR inlet weir was replaced with an inlet gate after an investigation of inlet design.Thirdly, the impact of a combined inlet/outlet operating procedure for the DR was investigated using the rainfall-runoff model SWMM and compared to the current operation and outlet only operation.Finally, the operating level in the DR was determined using a system resilience index and best hydraulic section for an urban drainage system in the target watershed.
Synthetic rainfall data were generated using the Huff distribution, and the most suitable quartile of synthetic rainfall data for the target watershed was selected.Rainfall durations of 30, 60, and 90 min were selected for a rainfall-runoff simulation to find the first flooding nodes considering the time of concentration in the target network.Total rainfall was increased in 1 mm increments until first flooding occurred.Node 560 was the first flooding node in 30-and 60-min duration simulations and this node was selected as a monitoring node.Maximum flooding nodes were identified using historical rainfall data.In two historical rainfall events, node 550 was the maximum flooding node; hence, this node was selected as a monitoring node.The DR inlet node was selected as the final monitoring node.
For combined inlet/outlet operation in the DR, the inlet weir was replaced with an inlet gate.The The ratio between partial flow Q and full flow Q F is shown in Equation (10).
Equation ( 8) is differentiated about θ to obtain the maximum discharge (Q max ) of a conduit.Q max is calculated when θ is 5.28 radian (θ = 303 • ).The ratio of d and D is 0.938 when θ is 303 • .Hence, the depth of maximum discharge in the monitoring nodes (d) is 1.407 m because the total depth of monitoring nodes (D) is 1.5 m.The range of the ratio of d and D when Q is larger than Q max is from 0.81 to 1.00 and the range of monitoring nodes is approximately from 1.2 m to 1.5 m.The level of conduits also gratifies the qualification of best hydraulic section if the level of the monitoring nodes satisfies the condition of best hydraulic section.
Engineers should consider a safety factor including differences between real discharge and design discharge when they design urban drainage conduits in South Korea.Safety factors for urban drainage conduits are based on the diameter of the conduit.For example, regulations stipulate that large conduits with diameters over 1500 mm should have an excess capacity of between 25% and 50% (Korea Water and Wastewater Association, 2011) [44].Safety factors in monitoring nodes range from 1.25 to 1.5.The operating level including a safety factor is from 1.0 m (safety factor = 1.5 m/1.5 = 1.0 m) to 1.2 m (safety factor = 1.5 m/1.25 = 1.2 m).The possible operating level in the DR based on the resilience index, best hydraulic section, and local regulation is shown in Table 3.

Conclusions
This study suggested combined changes to inlet design and inlet/outlet operation of DR in urban drainage systems to reduce the incidence of urban inundation.This study was comprised of four steps.Firstly, monitoring nodes for operating DR inlet and outlet were selected.Secondly, the DR inlet weir was replaced with an inlet gate after an investigation of inlet design.Thirdly, the impact of a combined inlet/outlet operating procedure for the DR was investigated using the rainfall-runoff model SWMM and compared to the current operation and outlet only operation.Finally, the operating level in the DR was determined using a system resilience index and best hydraulic section for an urban drainage system in the target watershed.
Synthetic rainfall data were generated using the Huff distribution, and the most suitable quartile of synthetic rainfall data for the target watershed was selected.Rainfall durations of 30, 60, and 90 min were selected for a rainfall-runoff simulation to find the first flooding nodes considering the time of concentration in the target network.Total rainfall was increased in 1 mm increments until first flooding occurred.Node 560 was the first flooding node in 30-and 60-min duration simulations and this node was selected as a monitoring node.Maximum flooding nodes were identified using historical rainfall data.In two historical rainfall events, node 550 was the maximum flooding node; hence, this node was selected as a monitoring node.The DR inlet node was selected as the final monitoring node.
For combined inlet/outlet operation in the DR, the inlet weir was replaced with an inlet gate.The inlet gate is open and outlet drainage pumps stop when the water level of monitoring nodes is high.Conversely, the inlet gate is closed and outlet drainage pumps are operational when the water level of monitoring nodes is low.Rainfall-runoff simulations for the current operation, outlet operation, and combined inlet/outlet operation in the DR were conducted.The flooding volume of outlet operation is generally smaller than that of current operation.The exceptional cases arise based on two reasons; the frequency of applied rainfall data and reserve gradient of some sewer conduits.The frequency distributions of the 2010 and 2011 rainfall events are approximately 100 years.Extreme rainfall with a 100-year return period produces exceptional results because the maximum design frequency of urban drainage system in South Korea is 30 years.Flooding volume using the combined inlet/outlet operation is considerably smaller than the flooding volume using the current operation and outlet operation methods.
DR operating level was determined using the system resilience index and best hydraulic section.Average system resilience was higher than total average system resilience when the operating level was between 1.2 m and 1.5 m.The best hydraulic section in circular conduits is when the ratio (d/D) between water level (d) and total level (D) is 0.938.The range of water levels that produce a discharge greater than Q F with full flow was 1.2 m to 1.5 m.Hence, the operating level of the Daerim detention reservoir was determined to be 1.2 m, accounting for local regulations.This study highlighted the benefits of inlet/outlet operation in the DR and determined the operating level for the DR.Follow up research based on these findings is planned including the operation of multiple detention reservoirs, determination of monitoring nodes using entropy, and integrated operations between river and urban drainage systems.

Figure 1 .
Figure 1.Flowchart for new detention reservoir (DR) operation and inlet design.

Figure 1 .
Figure 1.Flowchart for new detention reservoir (DR) operation and inlet design.

Figure 2 .
Figure 2. Structure of a DR in an urban drainage system.

Figure 2 .
Figure 2. Structure of a DR in an urban drainage system.
One of tributaries of the Han River is the Anyang Stream, which also has several tributaries including the Dorim Stream.The Dorim Stream has two detention reservoirs (Daerim detention reservoir and Gwanak detention reservoir), ten pump stations (Guro 1, Guro 2, Guro 3, Guro 4, Mullae, Dorim 1, Daerim 2, Daerim 3, Sinlim, and Sinlim 5), and two tributaries (Daebang Stream and Bongchun Stream).The Daerim detention reservoir within the urban drainage system of the Daerim 3 pump station was selected for this study.The drainage system of the Daerim 3 pump station is adjacent to the drainage system of Daerim 2 pump station.In the 1990s, two drainage areas including Daerim 3 and Daerim 2 pump stations were combined into one drainage area and Daerim 3 pump station drained all the flow in target watershed.The drainage area was divided into two drainage areas and Daerim 2 pump station was established because it is an area where flooding is frequent.The drainage area of Daerim 3 pump station including Daerim detention reservoir is 248.5 ha.The Daerim detention reservoir was constructed upstream of the drainage network in 2009 [17].Water 2017, 9, 246 6 of 14 reservoir within the urban drainage system of the Daerim 3 pump station was selected for this study.The drainage system of the Daerim 3 pump station is adjacent to the drainage system of Daerim 2 pump station.In the 1990s, two drainage areas including Daerim 3 and Daerim 2 pump stations were combined into one drainage area and Daerim 3 pump station drained all the flow in target watershed.The drainage area was divided into two drainage areas and Daerim 2 pump station was established because it is an area where flooding is frequent.The drainage area of Daerim 3 pump station including Daerim detention reservoir is 248.5 ha.The Daerim detention reservoir was constructed upstream of the drainage network in 2009 [17].

Figure 3 .
Figure 3. Location and features of the target watershed in Seoul (from Google Earth Imagery © Google Inc., 2015).

Figure 3 .
Figure 3. Location and features of the target watershed in Seoul (from Google Earth Imagery © Google Inc., 2015).

Figure 4 .
Figure 4. Location of the first and maximum flooding nodes in the main conduits.

Figure 4 .
Figure 4. Location of the first and maximum flooding nodes in the main conduits.
0 D to 1.0 D. For example, 1.5 m is same with 1.0 D if the depth of the monitoring node is 1.5 m and 1.2 m is same with 1.0 D if the depth of the monitoring node is 1.2 m.

Figure 5 .
Figure 5. Flooding volumes compared to the level of monitoring nodes for three DR operation methods: current operation (black); outlet operation (grey); and combined inlet/outlet operation (white).

Figure 5 .
Figure 5. Flooding volumes compared to the level of monitoring nodes for three DR operation methods: current operation (black); outlet operation (grey); and combined inlet/outlet operation (white).

Figure 7 .
Figure 7. Schematic of uniform flow in a circular conduit.

Figure 7 .
Figure 7. Schematic of uniform flow in a circular conduit.

Table 2 .
System resilience based on Huff distribution synthetic rainfall data.

Table 3 .
DR operating level based on the system resilience index, best hydraulic section and local regulations.

Table 3 .
DR operating level based on the system resilience index, best hydraulic section and local regulations.