Performance assessment of two-dimensional hydraulic models for generation of flood inundation maps in mountain river basins

Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.

HIGH RESOLUTION PARAMETRIC MODELLING FOR TWO-DIMENSIONAL RADAR TARGET USING PRONY ALGORITHM

On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high resolution estimation algorithm to extract the model parameters is also developed by building the relation of the scattering model and Prony model. The analysis of Cramer-Rao bound and simulation show that the method here has better statistical performance. The simulated analysis also indicates that the accurate extraction of the diffraction coefficient of scattering center is restricted by signal to noise ratio, radar center frequency and radar bandwidth.

AN EXPLORATION TO THE MODELLING OF TWO-DIMENSIONAL COMPLEX TURBULENT SHEAR FLOWS

The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.

基于MIKE FLOOD模型的平原水网地区内涝治理效果模拟研究

为科学评估平原水网地区内涝治理效果,以昆山市淀山湖镇为研究对象,基于MIKE FLOOD模型,耦合一维河网、二维地形和一维管网模拟了平原水网地区城镇内涝情况。采用实测降雨数据对模型进行了率定和验证,并模拟了2年一遇、5年一遇、10年一遇和20年一遇4种重现期降雨条件下现状和治理后的最大淹没水深。模拟结果表明:20年一遇降雨条件下内涝范围变化明显,其余3种重现期降雨条件下内涝范围变化不明显;对于降雨重现期较小条件下的城市内涝治理,应优先考虑管网改造,而针对短历时的强暴雨事件,应优先考虑海绵城市设施建设改造;为了应对不同重现期、短历时降雨造成的城市内涝风险,应在管网改造的基础上配合海绵城市设施建设改造。

Development and application of a two-dimensional water quality model for the Daqinghe River Mouth of the Dianchi Lake

Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.

Simulation of two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Full 2D Hydrodynamic Modelling of Rainfall-induced Flash Floods

Mountain catchments are prone to flash flooding due to heavy rainfall. Enhanced understanding of the generation and evolution processes of flash floods is essential for effective flood risk management. However, traditional distributed hydrological models based on kinematic and diffusion wave approximations ignore certain physical mechanisms of flash floods and thus bear excessive uncertainty. Here a hydrodynamic model is presented for flash floods based on the full two-dimensional shallow water equations incorporating rainfall and infiltration. Laboratory experiments of overland flows were modelled to illustrate the capability of the model. Then the model was applied to resolve two observed flash floods of distinct magnitudes in the Lengkou catchment in Shanxi Province, China. The present model is shown to be able to reproduce the flood flows fairly well compared to the observed data. The spatial distribution of rainfall is shown to be crucial for the modelling of flash floods. Sensitivity analyses of the model parameters reveal that the stage and discharge hydrographs are more sensitive to the Manning roughness and initial water content in the catchment than to the Green-Ampt head. Most notably, as the flash flood augments due to heavier rainfall, the modelling results agree with observed data better, which clearly characterizes the paramount role of rainfall in dictating the floods. From practical perspectives, the proposed model is more appropriate for modelling large flash floods.

Modelling Chorabari Lake outburst flood, Kedarnath, India

In this study, the Glacier Lake Outburst Flood(GLOF) that occurred over Kedarnath in June 2013 was modeled using integrated observations from the field and Remote Sensing(RS). The lake breach parameters such as area, depth, breach, and height have been estimated from the field observations and Remote Sensing(RS) data. A number of modelling approaches, including Snow Melt Runoff Model(SRM), Modified Single Flow model(MSF), Watershed Management System(WMS), Simplified Dam Breach Model(SMPDBK) and BREACH were used to model the GLOF. Estimations from SRM produced a runoff of about 22.7 m3 during 16–17, June 2013 over Chorabari Lake. Bathymetry data reported that the lake got filled to its maximum capacity(3822.7 m3) due to excess discharge. Hydrograph obtained from the BREACH model revealed a peak discharge of about 1699 m3/s during an intense water flow episode that lasted for 10–15 minutes on 17 th June 2013. Excess discharge from heavy rainfall and snowmelt into the lake increased its hydrostatic pressure and the lake breached cataclysmically.

Analytical model of 2D leakoff in waterflood-induced fractures

Waterflood-induced fractures,also known as self-induced fractures,spontaneously form at injection wells during waterflooding.These fractures propagate long distances through rock,allowing injected fluids to travel far away from a well,both within and outside the flooding layer.Essentially,the me-chanics of waterflood-induced fracture propagation is similar to that of hydraulic fractures,which are intentionally created for reservoir stimulation.Fracturing models developed for hydraulic fractures can also be applied to waterflood-induced fractures.However,waterflood-induced fractures are typically pumped with much larger volumes of water or brine and grow much longer in time.As a result,fluid leakoff from waterflood fractures into the formation is more extensive and two-dimensional(2D),a characteristic that is often ignored in a majority of modern fracturing simulators,making their appli-cation to waterflood fractures unreliable.In this work,we revisit the problem of leakoff for long-growing waterflood-induced fractures and develop a new analytical model for fluid leakoff that provides improved predictions of fracture geometry and can be easily implemented in fracturing simulators.We incorporate the developed solution into the classical Perkins-Kern-Nordgren(PKN)model of fracture growth,which shows that the choice of the Carter or a 2D leakoff model greatly impacts fracture ge-ometry at large time.The conducted parametric study shows while a toughness-dominated regime af-fects fracture evolution,most of fracture lifetime occurs in a viscosity-and-leakoff-dominated regime.We also develop an asymptotic solution for a leakoff profile in the limiting case of 2D leakoff domination(~~M and~~K).Finally,we study 3D fracture growth and out-of-zone injection with three layers and a complex structure of zones.The study shows that ignoring the 2D leakoff during simulation results in a significant overestimation of fracture geometry predictions.The present work,thus,plays an important role in improving waterflood fracture modelling,as it highlights the significance of 2D leakoff in waterflood-induced fractures and provides a reliable analytical model for fluid leakoff that can be incorporated into modern fracture simulators.

Flood Forecasting and Warning System: A Survey of Models and Their Applications in West Africa

Flood events occurrences and frequencies in the world are of immense worry for the stability of the economy and life safety. Africa continent is the third continent the most negatively affected by the flood events after Asia and Europe. Eastern Africa is the most hit in Africa. However, Africa continent is at the early stage in term of flood forecasting models development and implementation. Very few hydrological models for flood forecasting are available and implemented in Africa for the flood mitigation. And for the majority of the cases, they need to be improved because of the time evolution. Flash flood in Bamako (Mali) has been putting both human life and the economy in jeopardy. Studying this phenomenon, as to propose applicable solutions for its alleviation in Bamako is a great concern. Therefore, it is of upmost importance to know the existing scientific works related to this situation in Mali and elsewhere. The main aim was to point out the various solutions implemented by various local and international institutions, in order to fight against the flood events. Two types of methods are used for the flood events adaptation: the structural and non-structural methods. The structural methods are essentially based on the implementation of the structures like the dams, dykes, levees, etc. The problem of these methods is that they may reduce the volume of water that will inundate the area but are not efficient for the prediction of the coming floods and cannot alert the population with any lead time in advance. The non-structural methods are the one allowing to perform the prediction with acceptable lead time. They used the hydrological rainfall-runoff models and are the widely methods used for the flood adaptation. This review is more accentuated on the various types non-structural methods and their application in African countries in general and West African countries in particular with their strengths and weaknesses. Hydrologiska Byråns Vattenbalansavdelning (HBV), Hydrologic Engineer Center Hydrologic Model System (HEC-HMS) and Soil and Water Assessment Tool (SWAT) are the hydrological models that are the most widely used in West Africa for the purpose of flood forecasting. The easily way of calibration and the weak number of input data make these models appropriate for the West Africa region where the data are scarce and often with bad quality. These models when implemented and applied, can predict the coming floods, allow the population to adapt and mitigate the flood events and reduce considerably the impacts of floods especially in terms of loss of life.

Establishment of Unstable Flow Model and Well Testing Analysis for Viscoelastic Polymer Flooding

At present, the polymer solution is usually assumed to be Newtonian fluid or pseudoplastic fluid, and its elasticity is not considered on the study of polymer flooding well testing model. A large number of experiments have shown that polymer solutions have viscoelasticity, and disregarding the elasticity will cause certain errors in the analysis of polymer solution seepage law. Based on the percolation theory, this paper describes the polymer flooding mechanism from the two aspects of viscous effect and elastic effect, the mathematical model of oil water two-phase three components unsteady flow in viscoelastic polymer flooding was established, and solved by finite difference method, and the well-test curve was drawn to analyze the rule of well test curve in polymer flooding. The results show that, the degree of upward warping in the radial flow section of the pressure recovery curve when considering polymer elasticity is greater than the curve which not considering polymer elasticity. The relaxation time, power-law index, polymer injection concentration mainly affect the radial flow stage of the well testing curve. The relaxation time, power-law index, polymer injection concentration and other polymer flooding parameters mainly affect the radial flow stage of the well testing curve. The larger the polymer flooding parameters, the greater the degree of upwarping of the radial flow derivative curve. This model has important reference significance for well-testing research in polymer flooding oilfields.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

Global Two-Dimensional Chemistry Model and Simulation of Atmospheric Chemical Composition

A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.

基于FloodArea模型的呼和浩特市城区大青山南麓山洪致灾临界面雨量研究

以呼和浩特市城区北部大青山南麓山洪沟流域为研究区域,利用流域内气象观测站逐小时降水数据、数字高程模型数据、土地利用数据等资料,针对1998年7月12日和2020年8月1日暴雨时段,采用FloodArea淹没模型进行洪水动态淹没模拟。结果表明:1998年7月12日暴雨过程的FloodArea水动力模型模拟结果与实际情况较为吻合,且9条山洪沟的洪水淹没深度与累计9 h滑动面雨量的相关程度最高,构建了面雨量与淹没深度的回归关系方程,基于隐患点3个风险等级,确定了不同等级下的临界面雨量,以红山口沟为例,分别为69.9 mm(1级)、39.3 mm(2级)、17.2 mm(3级)。2020年8月1日暴雨过程的模型模拟结果基本与山洪过程一致。

Comprehensive two-dimensional river ice model based on boundary-fitted coordinate transformation method

River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal processes, and they are also influenced by weather and hydrologic conditions. Because natural rivers are serpentine, with bends, narrows, and straight reaches, the commonly-used one-dimensional river ice models and two-dimensional models based on the rectangular Cartesian coordinates are incapable of simulating the physical phenomena accurately. In order to accurately simulate the complicated river geometry and overcome the difficulties of numerical simulation resulting from both complex boundaries and differences between length and width scales, a two-dimensional river ice numerical model based on a boundary-fitted coordinate transformation method was developed. The presented model considers the influence of the frazil ice accumulation under ice cover and the shape of the leading edge of ice cover during the freezing process. The model is capable of determining the velocity field, the distribution of water temperature, the concentration distribution of frazil ice, the transport of floating ice, the progression, stability, and thawing of ice cover, and the transport, accumulation, and erosion of ice under ice cover. A MacCormack scheme was used to solve the equations numerically. The model was validated with field observations from the Hequ Reach of the Yellow River. Comparison of simulation results with field data indicates that the model is capable of simulating the river ice process with high accuracy.

Spatial diversion and coordination of flood water for an urban flood control project in Suzhou, China

Suzhou City,located in the Yangtze River Delta in China,is prone to flooding due to a complex combination of natural factors,including its monsoon climate,low elevation,and tidally influenced position,as well as intensive human activities.The Large Encirclement Flood Control Project(LEFCP)was launched to cope with serious floods in the urban area.This project changed the spatiotemporal pattern of flood processes and caused spatial diversion of floods from the urban area to the outskirts of the city.Therefore,this study developed a distributed flood simulation model in order to understand this transition of flood processes.The results revealed that the LEFCP effectively protected the urban areas from floods,but the present scheduling schemes resulted in the spatial diversion of floods to the outskirts of the city.With rainstorm frequencies of 10.0%to 0.5%,the water level differences between two representative water level stations(Miduqiao(MDQ)and Fengqiao(FQ))located inside and outside the LEFCP area,ranged from 0.75 m to 0.24 m and from 1.80 m to 1.58 m,respectively.In addition,the flood safety margin at MDQ and the duration with the water level exceeding the warning water level at FQ ranged from 0.95 m to 0.43 m and from 4 h to 22 h,respectively.Rational scheduling schemes for the hydraulic facilities of the LEFCP in extreme precipitation cases were developed ac-cording to food simulations under seven scheduling scenarios.This helps to regulate the spatial flood diversion caused by the LEFCP during extreme precipitation.

Flood Modeller在河流洪水风险管理中的应用

本文构建了某河流的洪水建模器(Flood Modeller)一维、二模模型,通过与实际洪水数据对比,校准模型的曼宁系数n值以提高模拟准确性。使用50年、100年和200年一遇的洪水重现期模拟洪水事件,评估潜在淹没区域,为洪水风险评估提供科学的可视化工具。

Numerical modeling of blast-induced rock fragmentation in deep mining with 3D and 2D FEM method approaches

To optimize the excavation of rock using underground blasting techniques,a reliable and simplified approach for modeling rock fragmentation is desired.This paper presents a multistep experimentalnumerical methodology for simplifying the three-dimensional(3D)to two-dimensional(2D)quasiplane-strain problem and reducing computational costs by more than 100-fold.First,in situ tests were conducted involving single-hole and free-face blasting of a dolomite rock mass in a 1050-m-deep mine.The results were validated by laser scanning.The craters were then compared with four analytical models to calculate the radius of the crushing zone.Next,a full 3D model for single-hole blasting was prepared and validated by simulating the crack length and the radius of the crushing zone.Based on the stable crack propagation zones observed in the 3D model and experiments,a 2D model was prepared.The properties of the high explosive(HE)were slightly reduced to match the shape and number of radial cracks and crushing zone radius between the 3D and 2D models.The final methodology was used to reproduce various cut-hole blasting scenarios and observe the effects of residual cracks in the rock mass on further fragmentation.The presence of preexisting cracks was found to be crucial for fragmentation,particularly when the borehole was situated near a free rock face.Finally,an optimization study was performed to determine the possibility of losing rock continuity at different positions within the well in relation to the free rock face.

DETERMINATION OF PARAMETERS IN A TWO-DIMENSIONAL MODEL FOR FIXED-BED ADSORBERS

A new method is proposed for the determination of the parameters in a two-dimensionalmodel which characterizes the properties of axial and radial mixing and mass transport in afixed-bed adsorber.Parameter estimation for the model is carried out with methane-air-5A molecularsieve in a bed under the condition of step injection of tracer from a point on the main axis of thebed by the curve fitting method in the time domain.

Two-dimensional hydrodynamic robust numerical model of soil erosion based on slopes and river basins

Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.