Hydrologic models coupled with 2D hydrodynamic model for high-resolution urban flood simulation

Natural Hazards - Floods are the most frequent natural disaster and pose a very challenging threat to many cities worldwide. Understanding the flood dynamic is essential for developing strategies...     

一维-二维耦合的防洪保护区洪水演进数学模型

针对防洪保护区溃堤及漫堤洪水演进数值模拟面临复杂计算域、河道-防洪保护区洪水耦合作用等问题,建立了基于侧向联解的一维-二维耦合水动力学模型.通过构造并求解Riemann问题实现一维-二维模型耦合,有效克服了基于堰流公式的传统方法难以处理模型间动量交换的缺点,也避免了堰流公式中流量系数选取的不确定性;提出了时间步长自适应匹配方法,解决了一维模型和二维模型时间步长不一致问题.算例研究表明,该模型可有效模拟河道-防洪保护区耦合系统中漫堤洪水和具有任意溃口形状的溃堤洪水演进过程,具有较好的推广应用价值.     

基于WCA2D与SWMM模型的城市暴雨洪涝快速模拟

为探讨计算高效的元胞自动机模型(WCA2D)与传统一维管网模型耦合的机制以及计算效果,尝试将WCA2D与SWMM模型耦合(SWMM/WCA2D),以广州市长湴片区为例探究一种暴雨洪涝快速二维模拟技术,对比实测积水数据以及SWMM/LISFLOOD-FP模拟,结果表明:SWMM/WCA2D模拟结果与"20180607"实测积水数据相近,表明模型精度良好;根据多指标评估结果,综合考虑主干渠道淤积以及建筑物阻挡情景的R_TPR、R_PPV、F₁值分别达到0.8、0.6、0.7,模拟精度最高,最能反映区域实际情况;通过与SWMM耦合,WCA2D和LISFLOOD-FP的模拟结果差异小(最大水深差值基本低于0.1 m)、相关性强(相关系数基本超过0.7),但前者计算效率是后者的3～5倍,表明WCA2D能够耦合SWMM且计算效率更高,为复杂城市化地区暴雨洪涝快速模拟提供了一种新方法。     

Multiple scenario analyses of Huangpu River flooding using a 1D/2D coupled flood inundation model

Huangpu River floodplain is historically vulnerable to flooding due to its location in the path of tropical cyclones, low elevation, relatively flat topography, rapid changes in sea level and fast rate of land subsidence due to urbanization. This paper presents a scenario-based study that investigates the fluvial flood potentials in the Huangpu River floodplain. Flood scenarios with return periods of 50, 100, 200, 500 and 1,000 years were designed to cover the probable situations. Further, a flood inundation model (FloodMap) that tightly couples a river flow model with a 1D solution of the full form of the St. Venant equations and a 2D floodplain flow model was used to predict the river flow and inundation extents. Flood characteristics obtained from the simulations were used in the exposure analysis to determine the spatial distribution of susceptible land uses under different scenarios. Results suggest that overtopping inundation mainly occurs within 1–2 km of the banks of the Huangpu River, with larger inundation extent predicted in the upper and middle reaches of the channel, a result of varying protection levels from relatively rural upstream to high urbanized floodplain in the vicinity of the middle reaches.     

太湖流域应对特大洪水防洪工程效益模拟

防洪效益评估对防洪工程投资决策与减灾对策制定具有重要意义。建立集成了与太湖流域防洪效益评估相关的系列模型和方法,包括含降雨产流与平原净雨计算的水文分析方法、由河网水动力学模型和平原区域洪水分析模型组成的大尺度水力学模型、综合流域社会经济和淹没因素的洪灾损失评估模型。模拟了太湖流域遇特大洪水的灾害损失,开展了不同防洪工程应对流域性特大洪水减灾效益的预测分析。结果表明：1999年型200年一遇降雨将会给太湖流域造成高达568.29亿元的直接经济损失,外排动力增强30%至100%的防洪效益介于26.69亿元到45.70亿元之间,新建圩区、太浦河拓宽的防洪效益依次减小,而圩区泵排能力增加30%的防洪效益仅为0.65亿元。基于研究成果提出了增设外排泵站、加强圩区科学调度、通过保险分担风险等应对特大洪水的对策措施建议,为太湖流域特大洪水的防治提供支撑和参考。     

太湖流域应对特大洪水防洪工程效益模拟

防洪效益评估对防洪工程投资决策与减灾对策制定具有重要意义。建立集成了与太湖流域防洪效益评估相关的系列模型和方法,包括含降雨产流与平原净雨计算的水文分析方法、由河网水动力学模型和平原区域洪水分析模型组成的大尺度水力学模型、综合流域社会经济和淹没因素的洪灾损失评估模型。模拟了太湖流域遇特大洪水的灾害损失,开展了不同防洪工程应对流域性特大洪水减灾效益的预测分析。结果表明:1999年型200年一遇降雨将会给太湖流域造成高达568.29亿元的直接经济损失,外排动力增强30%至100%的防洪效益介于26.69亿元到45.70亿元之间,新建圩区、太浦河拓宽的防洪效益依次减小,而圩区泵排能力增加30%的防洪效益仅为0.65亿元。基于研究成果提出了增设外排泵站、加强圩区科学调度、通过保险分担风险等应对特大洪水的对策措施建议,为太湖流域特大洪水的防治提供支撑和参考。     

一维-二维耦合的河湖系统整体水动力模型

为对自然界中的复杂水流系统进行整体水动力数值模拟,建立了适用于河湖系统的一维-二维耦合模型.将河网计算中用于处理河段间耦合的汊点水位预测校正法应用到一维-二维耦合边界的处理中,即耦合边界水位预测校正法.该方法具有一维、二维模型计算完全独立,可有效利用现有模型的优点.与传统一维-二维耦合处理方法相比,耦合连接条件的满足程度较高且可通过对计算容差的调整进行有效控制.利用理想及实际例子进行了计算验证,结果表明该方法具有较好的实用性.     

A model for coupled mechanical and hydraulic behaviour of a rock joint

Constitutive laws for rock joints should be able to reproduce the fundamental mechanical behaviour of real joints, such as dilation under shear and strain softening due to surface asperity degradation. In this work, we extend the model of Plesha to include hydraulic behaviour. During shearing, the joint can experience dilation, leading to an initial increase in its permeability. Experiments have shown that the rate of increase of the permeability slows down as shearing proceeds, and, at later stages, the permeability could decrease again. The above behaviour is attributed to gouge production. The stress–strain relationship of the joint is formulated by appeal to classical theories of interface plasticity. It is shown that the parameters of the model can be estimated from the Barton–Bandis empirical coefficients; the Joint Roughness Coefficient (JRC) and the Joint Compresive strength (JSC). We further assume that gouge production is also related to the plastic work of the shear stresses, which enables the derivation of a relationship between the permeability of the joint and its mechanical aperture. The model is implemented in a finite element code (FRACON) developed by the authors for the simulation of the coupled thermal–hydraulic–mechanical behaviour of jointed rock masses. Typical laboratory experiments are simulated with the FRACON code in order to illustrate the trends predicted in the proposed model. © 1998 by John Wiley & Sons. Ltd.     

Development of roughness updating based on artificial neural network in a river hydraulic model for flash flood forecasting

Flood is the worst weather-related hazard in Taiwan because of steep terrain and storm. The tropical storm often results in disastrous flash flood. To provide reliable forecast of water stages in rivers is indispensable for proper actions in the emergency response during flood. The river hydraulic model based on dynamic wave theory using an implicit finite-difference method is developed with river roughness updating for flash flood forecast. The artificial neural network (ANN) is employed to update the roughness of rivers in accordance with the observed river stages at each time-step of the flood routing process. Several typhoon events at Tamsui River are utilized to evaluate the accuracy of flood forecasting. The results present the adaptive n-values of roughness for river hydraulic model that can provide a better flow state for subsequent forecasting at significant locations and longitudinal profiles along rivers.     

Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.     

Development of a full 3D numerical model to investigate the hydraulic fracture propagation under the impact of orthogonal natural fractures

Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.

     

Generation of 1D and 2D flood maps of Sakarya river passing through Geyve district of Sakarya city in Turkey

Natural Hazards - In recent years, changes in precipitation caused by climate change resulting from global warming lead to great floods in world. This climate change causes drought in some regions,...     

流域防洪效益计算宏观模型

针对现有防洪效益计算模型所需资料多、计算复杂等问题,本文在符合防洪效益特点的宏观假定基础上,提出了防洪效益计算的洪水频率面积系数和洪水频率严重系数,并用水文学和数理统计的有关原理导出了流域防洪系统效益计算的宏观模型。算例研究的结果表明:模型理论上合理,实际应用简单,可操作性强。     

某水力枢纽基坑开挖过程渗流耦合数值模拟

岩土工程中存在多种耦合作用,在水利水电工程中,渗流场与应力场的耦合是最重要的。以某水利枢纽厂房基坑开挖为工程背景,比奥固结理论为理论基础,莫尔-库仑准则为土体屈服判据,对基坑开挖过程中渗流耦合进行了数值模拟。计算结果显示：回填后的应力略大于开挖过程中的应力,沉井附近出现一定程度的应力集中。回填前后塑性区变化不大且范围较小,因而土体的屈服破坏不是该工程的重点。由于土体渗透系数大,厂房区地下水位高,基坑开挖过程中渗水量较大,回填后沉井及基坑底部防渗底板的防渗效果显著,因而在基坑开挖过程中防渗排水是该工程的重点问题。     

On the potential for CO2 mineral storage in continental flood basalts - PHREEQC batch- and 1D diffusion-reaction simulations

ABSTRACT: Continental flood basalts (CFB) are considered as potential CO2 storage sites because of their high reactivity and abundant divalent metal ions that can potentially trap carbon for geological timescales. Moreover, laterally extensive CFB are found in many place in the world within reasonable distances from major CO2 point emission sources.Based on the mineral and glass composition of the Columbia River Basalt (CRB) we estimated the potential of CFB to store CO2 in secondary carbonates. We simulated the system using kinetic dependent dissolution of primary basalt-minerals (pyroxene, feldspar and glass) and the local equilibrium assumption for secondary phases (weathering products). The simulations were divided into closed-system batch simulations at a constant CO2 pressure of 100?bar with sensitivity studies of temperature and reactive surface area, an evaluation of the reactivity of H2O in scCO2, and finally 1D reactive diffusion simulations giving reactivity at CO2 pressures varying from 0 to 100?bar.Although the uncertainty in reactive surface area and corresponding reaction rates are large, we have estimated the potential for CO2 mineral storage and identified factors that control the maximum extent of carbonation. The simulations showed that formation of carbonates from basalt at 40?C may be limited to the formation of siderite and possibly FeMg carbonates. Calcium was largely consumed by zeolite and oxide instead of forming carbonates. At higher temperatures (60 - 100?C), magnesite is suggested to form together with siderite and ankerite. The maximum potential of CO2 stored as solid carbonates, if CO2 is supplied to the reactions unlimited, is shown to depend on the availability of pore space as the hydration and carbonation reactions increase the solid volume and clog the pore space. For systems such as in the scCO2 phase with limited amount of water, the total carbonation potential is limited by the amount of water present for hydration of basalt.     

Assessment of flood inundation mapping of Surat city by coupled 1D/2D hydrodynamic modeling: a case application of the new HEC-RAS 5

Surat city of India, situated 100 km downstream of Ukai Dam and 19.4 km upstream from the mouth of River Tapi, has experienced the largest flood in 2006. The peak discharge of about 25,770 m³ s^?1 released from the Ukai Dam was responsible for a disaster. To assess the flood and find inundation in low-lying areas, simulation work is carried out under the 1D/2D couple hydrodynamic modeling. Two hundred ninety-nine cross sections, two hydraulic structures and five major bridges across the river are considered for 1D modeling, whereas a topographic map at 0.5 m contour interval was used to produce a 5 m grid, and SRTM (30 and 90 m) grid has been considered for Surat and the Lower Tapi Basin. The tidal level at the river mouth and the release from the Ukai Dam during 2006 flood are considered as the downstream and upstream boundaries, respectively. The model is simulated under the unsteady flow condition and validated for the year 2006. The simulated result shows that 9th August was the worst day in terms of flooding for Surat city and a maximum 75–77% area was under inundation. Out of seven zones, the west zone had the deepest flood and inundated under 4–5 m. Furthermore, inundation is simulated under the bank protection work (i.e., levees, retaining wall) constructed after the 2006 flood. The simulated results show that the major zones are safe against the inundation under 14,430 m³ s^?1 water releases from Ukai Dam except for the west zone. The study shows the 2D capability of new HEC-RAS 5 for flood inundation mapping and management studies.     

Application of 1D and 2D hydrodynamic modeling to study glacial lake outburst flood (GLOF) and its impact on a hydropower station in Central Himalaya

Natural Hazards - This study aims to examine a set of alternatives as remedial measures for the flood management purposes in order to protect the northern part of Nicosia which has been affected by...