Mathematical model for flood routing based on cellular automaton

Increasing frequency and severity of flooding have caused tremendous damage in China, requiring more essential countermeasures to alleviate the damage. In this study, the dynamic simulation property of a cellular automaton was used to make further progress in flood routing. In consideration of terrain＇s influence on flood routing, we regarded the terrain elevation as an auxiliary attribute of a two-dimensional cellular automaton in path selection for flood routing and developed a mathematical model based on a cellular automaton. A numerical case of propagation of an outburst flood in an area of the lower Yangtze River was analyzed with both the fixed-step and variable-step models. The results show that the flood does not spread simultaneously in all directions, but flows into the lower place first, and that the submerged area grows quickly at the beginning, but slowly later on. The final submerged areas obtained from the two different models are consistent, and the flood volume balance test shows that the flood volume meets the requirement of the total volume balance. The analysis of the case shows that the proposed model can be a valuable tool for flood routing.

Mathematical model for flood routing in Jingjiang River and Dongting Lake network

The main stream of the Yangtze River, Dongting Lake, and the river network in the Jingjiang reach of the Yangtze River constitute a complex water system. This paper develops a one-dimensional （l-D） mathematical model for flood routing in the river network Of the Jingjiang River and Dongting Lake using the explicit finite volume method. Based on observed data during the flood periods in 1996 and 1998, the model was calibrated and validated, and the results show that the model is effective and has high accuracy. In addition, the one-dimensional mathematical model for the river network and the horizontal two-dimensional （2-D） mathematical model for the Jingjiang flood diversion area were coupled to simulate the flood process in the Jingjiang River, Dongting Lake, and the Jingjiang flood diversion area. The calculated results of the coupled model are consistent with the practical processes. Meanwhile, the results show that the flood diversion has significant effects on the decrease of the peak water level at the Shashi and Chenjiawan hydrological stations near the flood diversion gates, and the effect is more obvious in the downstream than in the upstream.

Analysis of dynamic wave model for flood routing in natural rivers

Flooding is a common natural disaster that causes enormous economic, social, and human losses. Of various flood routing methods, the dynamic wave model is one of the best approaches for the prediction of the characteristics of floods during their propagations in natural rivers because all of the terms of the momentum equation are considered in the model. However, no significant research has been conducted on how the model sensitivity affects the accuracy of the downstream hydrograph. In this study, a comprehensive analysis of the input parameters 9f the dynamic wave model was performed through field applications in natural rivers and routing experiments in artificial channels using the graphical multi-parametric sensitivity analysis （GMPSA）. The results indicate that the effects of input parameter errors on the output results are more significant in special situations, such as lower values of Manning＇s roughness coefficient and/or a steeper bed slope on the characteristics of a design hydrograph, larger values of the skewness factor and/or time to peak on the channel characteristics, larger values of Manning＇s roughness coefficient and/or the bed slope on the space step, and lower values of Manning＇s roughness coefficient and/or a steeper bed slope on the time step and weighting factor.

Google Earth-based dynamic visualization system for storm surge flood routing

To describe the dynamic process of flood routing intuitively and realistically when storm surge disaster occurs,a method for ArcGIS data and Google Earth(GE) data integration is proposed,which realizes the importing and integrating of basic geographic information into GE. Based on SketchUp and AutoCAD software,threedimension(3D) visualization of seawall and other tidal defense structures is achieved. By employing Microsoft Foundation Class Library(MFC),the related system modules and storm surge flood routing dynamic visualization system are developed. Therefore,dynamic visualization of flood routing process and interactive query of submerged area and inundated depth are implemented. A practical application case study of Tianjin Binhai New Area provides decision-making support for coastal seawall planning and storm surge disaster prevention and reduction.

Flood routing model incorporating intensive streambed infiltration

Flood routing models are critical to flood forecasting and confluence calculations. In the streams that dry up and disconnect from groundwater, the streambed infiltration is intensive and has a significant effect on flood wave movement. Streambed infiltration should be considered in flood routing. A flood routing model incorporating intensive streambed infiltration is proposed. In the model a streambed infiltration simulation method based on soil infiltration theory is developed. In this method the Horton equation is used to calculate infiltration capacity. A trial-and-error method is developed to calculate infiltration rate and determine whether the flood wave can travel downstream. A formula is derived to calculate infiltration flow per unit length. The Muskingum-Cunge method with streambed infiltration flow as lateral outflow is used for flood routing. The proposed model is applied to the stream from the downstream of the Yuecheng Reservoir to the Caixiaozhuang Hydrometric Station in the Zhangwei River of the Haihe River Basin. Simulation results show that the accuracy of the model is high, and the infiltration simulation method can represent infiltration processes well. The proposed model is simple and practical for flood simulation and forecasting, and can be used in river confluence calculations in a rainfall-runoff model for arid and semiarid regions.

A NEW GENETIC SIMULATED ANNEALING ALGORITHM FOR FLOOD ROUTING MODEL

In this paper, a new approach, the Genetic Simulated Annealing (GSA), wasproposed for optimizing the parameters in the Muskingum routing model. By integrating the simulatedannealing method into the genetic algorithm, the hybrid method could avoid some troubles oftraditional methods, such as arduous trial-and-error procedure, premature convergence in geneticalgorithm and search blindness in simulated annealing. The principle and implementing procedure ofthis algorithm were described. Numerical experiments show that the GSA can adjust the optimizationpopulation, prevent premature convergence and seek the global optimal result. Applications to theNanyunhe River and Qingjiang River show that the proposed approach is of higher forecast accuracyand practicability.

POINTWISE AND UPWIND DISCRETIZATIONS OF SOURCE TERMS IN OPEN-CHANNEL FLOOD ROUTING

Upwind algorithms are becoming progressively popular for river flood routing due to their capability of resolving trans-critical flow regimes. For consistency, these algorithms suggest natural upwind discretization of the source term, which may be essential for natural channels with irregular geometry. Yet applications of these upwind algorithms to natural river flows are rare, and in such applications the traditional and simpler pointwise, rather than upwind discretization of the source term is used. Within the framework of a first-order upwind algorithm, this paper presents a comparison of upwind and pointwise discretizations of the source term. Numerical simulations were carried out for a selected irregular channel comprising a pool-riffle sequence Jn the River Lune, England with observed data. It is Shown that the impact of pointwise discretization, compared to the upwind, is appreciable mainly in flow zones with the Froude number closer to or larger than unity. The discrepancy due to pointwise and upwind discretizations of the source term is negligible in flow depth and hence in water surface elevation, but well manifested in mean velocity and derived flow quantities. Also the occurrence of flow reversal and equalisation over the pool-riffle sequence in response to increasing discharges is demonstrated.

Application of Computer Simulation in the Reservoir Flood Control

[Objective] To discuss the application of simulation technique in preliminary solution of reservoir flood control of.[Method] Based on the study of river system simulation theory,the InfoWorks RS system was adopted to build a digital analog model for flood routing in certain river system and the flood discharge was analyzed with the check criterion.[Result] The flood information system based on InfoWorks RS simulated the inundation of downstream during flood discharge and drew mutation pictures to indicate how hydraulic parameters of maximum runoff at each control section,runoff of each control section and flow velocity varied with time,how hydraulic elements mutated with time as well as the most likely inundated area analysis of lower reservoir.[Conclusion] To provide references for reducing serious loss caused by flood and making out preliminary solution to flood control.

Application of Muskingum routing method with variable parameters in ungauged basin

This paper describes a flood routing method applied in an ungauged basin, utilizing the Muskingum model with variable parameters of wave travel time K and weight coefficient of discharge x based on the physical characteristics of the river reach and flood, including the reach slope, length, width, and flood discharge. Three formulas for estimating parameters of wide rectangular, triangular, and parabolic cross sections are proposed. The influence of the flood on channel flow routing parameters is taken into account. The HEC-HMS hydrological model and the geospatial hydrologic analysis module HEC-GeoHMS were used to extract channel or watershed characteristics and to divide sub-basins. In addition, the initial and constant-rate method, user synthetic unit hydrograph method, and exponential recession method were used to estimate runoff volumes, the direct runoff hydrograph, and the baseflow hydrograph, respectively. The Muskingum model with variable parameters was then applied in the Louzigou Basin in Henan Province of China, and of the results, the percentages of flood events with a relative error of peak discharge less than 20% and runoff volume less than 10% are both 100%. They also show that the percentages of flood events with coefficients of determination greater than 0.8 are 83.33%, 91.67%, and 87.5%, respectively, for rectangular, triangular, and parabolic cross sections in 24 flood events. Therefore, this method is applicable to ungauged basins.

Hydraulic model with roughness coefficient updating method based on Kalman filter for channel flood forecast

A real-time channel flood forecast model was developed to simulate channel flow in plain rivers based on the dynamic wave theory. Taking into consideration channel shape differences along the channel, a roughness updating technique was developed using the Kalman filter method to update Manning＇s roughness coefficient at each time step of the calculation processes. Channel shapes were simplified as rectangles, triangles, and parabolas, and the relationships between hydraulic radius and water depth were developed for plain rivers. Based on the relationship between the Froude number and the inertia terms of the momentum equation in the Saint-Venant equations, the relationship between Manning＇s roughness coefficient and water depth was obtained. Using the channel of the Huaihe River from Wangjiaba to Lutaizi stations as a case, to test the performance and rationality of the present flood routing model, the original hydraulic model was compared with the developed model. Results show that the stage hydrographs calculated by the developed flood routing model with the updated Manning＇s roughness coefficient have a good agreement with the observed stage hydrographs. This model performs better than the original hydraulic model.

The characteristics of dike-break flood in different scenarios of the lower Yellow River based on numerical simulations

The lower Yellow River still faces the threat of flood due to the unusual precipitation caused by global environmental change, river channel sedimentation, hidden danger in the dike and unfavorable river regime of ＂hanging river＂. According to the characteristics of the dike-break flood of the Yellow River, this paper has simulated, in six different scenarios, the dike-break flood routing by inputting the terrain data, typical historical flood data and land use data of study area to two-dimensional unsteady flow model. The results show that： firstly, the routing process of flood will occupy other rivers on the way and return to the rivers after reaching the lower reaches; secondly, in the same river reach, flood inundating area of north band is bigger than that at corresponding location of south bank under the same historical flood; thirdly, it is different in the degree of flood inundation in different regions due to different geographical locations in flood plain; fourthly, the area of mainstream where flood is deep and flow velocity is quick is relatively smaller, but the area of non-mainstream, where flood is shallow and flow velocity is slow, is relatively big; and finally, the possible influenced area of the dike-break flood is 141,948 km^2.

High-efficiency calculation method for watershed rainfallrunoff routing using one-dimensional dynamic wave equations

Efficiency in solving the Saint-Venant equations for watershed rainfall-runoff routing is important in flood hydrology. This paper presents a high-efficiency numerical solution of one-dimensional dynamic wave equations(HEDWE) for watershed rainfall-runoff routing, in which the full momentum equation is written as a quadratic equation with only one unknown variable Q, water depth is derived from the continuity equation using the two-step predictorcorrector method, and the discrete scheme is the explicit upwind scheme. The results of numerical tests showed the HEDWE approach has several major advantages. 1) It is a stable numerical method, even for an initially dry area. 2) Its computational efficiency is higher than 4.76E+05 times/s. 3) It can be used for overland flow, river flow, and combinations thereof. The primary disadvantages of the HEDWE approach are its unsuitability for rapidly varying flow, such as dam-break floods.

A framework of numerical simulation on moraine-dammed glacial lake outburst floods

Glacial outburst floods（GLOFs） in alpine regions tend to be relatively complicated, multi-stage catastrophes, capable of causing significant geomorphologic changes in channel surroundings and posing severe threats to infrastructure and the safety and livelihoods of human communities. GLOF disasters have been observed and potential hazards can be foreseen due to the newly formed glacial lakes or the expansion of existing ones in the Poiqu River Basin in Tibet, China. Here we presented a synthesis of GLOF-related studies including triggering mechanism（s）, dam breach modeling, and flood routing simulation that have been employed to reconstruct or forecast GLOF hydrographs. We provided a framework for probability-based GLOFs simulation and hazard mapping in the Poiqu River Basin according to available knowledge. We also discussed the uncertainties and challenges in the model chains, which may form the basis for further research.

Directional Routing Algorithm for Deep Space Optical Network

With the development of science, economy and society, the needs for research and exploration of deep space have entered a rapid and stable development stage. Deep Space Optical Network(DSON) is expected to become an important foundation and inevitable development trend of future deepspace communication. In this paper, we design a deep space node model which is capable of combining the space division multiplexing with frequency division multiplexing. Furthermore, we propose the directional flooding routing algorithm(DFRA) for DSON based on our node model. This scheme selectively forwards the data packets in the routing, so that the energy consumption can be reduced effectively because only a portion of nodes will participate the flooding routing. Simulation results show that, compared with traditional flooding routing algorithm(TFRA), the DFRA can avoid the non-directional and blind transmission. Therefore, the energy consumption in message routing will be reduced and the lifespan of DSON can also be prolonged effectively. Although the complexity of routing implementation is slightly increased compared with TFRA, the energy of nodes can be saved and the transmission rate is obviously improved in DFRA. Thus the overall performance of DSON can be significantly improved.