基于面平均雨量误差修正的实时洪水预报修正方法

空间集总式水文模型的洪水预报精度会受到面平均雨量估计误差的严重影响.点雨量监测值的误差类型、误差大小以及流域的雨量站点密度和站点的空间分布都会影响到面平均雨量的计算.为提高实时洪水预报精度,本文提出了一种基于降雨系统响应曲线洪水预报误差修正方法.通过此方法估计降雨输入项的误差,从而提高洪水预报精度.此方法将水文模型做为输入和输出之间的响应系统,用实测流量和计算流量之间的差值做为信息,通过降雨系统响应曲线,使用最小二乘估计原理,对面平均雨量进行修正,再用修正后的面平均雨量重新计算出流过程.将此修正方法结合新安江模型使用理想案例进行检验,并应用于王家坝流域的16场历史洪水以及此流域不同雨量站密度的情况下,结果证明均有明显修正效果,且在雨量站密度较低时修正效果更加明显.该方法是一种结构简单且不增加模型参数和复杂度的实时洪水修正的新方法.     

黄土高原大理河流域水沙耦合模型应用研究

水沙模型是定量描述水沙关系及水沙规律的重要工具,现阶段国内外对于水沙模型的研究大都为基于某个典型流域的经验统计模型或基于流域大量基础资料的物理模型,极大限制了其使用范围及模拟精度.本文建立了结构与参数均具有物理意义的流域水沙耦合物理概念模型,其优点是物理概念清楚,模拟精度高,实用性强,易于深入研究泥沙基本规律.该模型将概念性水文模型和泥沙模型耦合,提出水流挟沙能力和土壤抗侵蚀能力概念,用对数曲线近似描述流域土壤抗侵蚀能力的空间变异性,在拜格诺河道水流悬移质泥沙公式基础上建立概念性沟蚀产沙公式,参照水流汇集相似性建立泥沙汇集演算公式.选取黄河中游大理河流域4个不同流域尺度的实际流域对模型进行应用检验,模拟结果表明,该模型的水流泥沙两部分均有很高的模拟精度,可以很好地模拟黄河中游地区不同流域尺度和年际尺度上的洪水过程和泥沙产生及输移过程,充分证明了该模型结构、参数和计算方法的合理性,可用于定量分析流域内各项水保措施的减水减沙效益及水沙关系变化趋势,对该模型的推广应用可做进一步分析研究.     

开都河流域融雪径流模拟研究

高山融雪是塔里木河源流区重要的产流方式,4个山区流域具有面积大、测站稀少、降雨与融雪混合补给径流和显著局部降雨等特征.以开都河流域为研究区,分析流域特征对SRM融雪径流模型参变量的影响,确定相应选取策略以提高融雪径流模拟预报精度,为相似流域融雪径流模拟提供参考.研究结果表明(i)气温输入控制模拟径流的整体趋势,对模拟精度起决定性作用.但测站日均气温数据通常不能代表流域同高程的平均水平,直接作为输入会导致很大误差.基于遥感积雪图和模拟结果分析,开都河流域选择0.5倍巴音布鲁克站日最大气温作为流域平均气温较为合理.(ii)由于雨量站稀少和局部降雨特征显著,计算各高程分带平均降水并不现实.将测站降雨乘以放大系数,并借助参数"降雨径流系数"进行校正,可以满足模型对降雨输入的需求.(iii)根据融雪和降雨位置变化,调整参数"滞时"取值对提高局部洪峰过程的模拟精度非常重要.(iv)随气温升高,降雨增多,未能被有限测站完全监测,导致模拟精度逐步降低.     

开都河流域融雪径流模拟研究

高山融雪是塔里木河源流区重要的产流方式,4个山区流域具有面积大、测站稀少、降雨与融雪混合补给径流和显著局部降雨等特征.以开都河流域为研究区,分析流域特征对SRM融雪径流模型参变量的影响,确定相应选取策略以提高融雪径流模拟预报精度,为相似流域融雪径流模拟提供参考.研究结果表明:(i)气温输入控制模拟径流的整体趋势,对模拟精度起决定性作用.但测站日均气温数据通常不能代表流域同高程的平均水平,直接作为输入会导致很大误差.基于遥感积雪图和模拟结果分析,开都河流域选择0.5倍巴音布鲁克站日最大气温作为流域平均气温较为合理.(ii)由于雨量站稀少和局部降雨特征显著,计算各高程分带平均降水并不现实.将测站降雨乘以放大系数,并借助参数"降雨径流系数"进行校正,可以满足模型对降雨输入的需求.(iii)根据融雪和降雨位置变化,调整参数"滞时"取值对提高局部洪峰过程的模拟精度非常重要.(iv)随气温升高,降雨增多,未能被有限测站完全监测,导致模拟精度逐步降低.     

基于土壤墒情抗差的降雨径流预报模型——设计与实验流域应用

本文旨在将实时监测得到的土壤墒情转化为流域水文模型可以直接使用的土壤含水量,论证将实时土壤墒情资料用于实时预报的可行性;利用实时监测土壤墒情,改进传统的模型结构,设计基于实测土壤墒情的降雨径流水文预报模型.采用土壤含水量误差抗差估计技术以抵御观测资料粗差的影响,提高系统的稳定性;并在此基础上提出了土壤含水量系统响应修正方法,以提高模型计算精度.将该模型应用于实验流域——宝盖洞流域进行应用检验,洪水模拟合格率达到92.3%,整体模拟精度达到甲级.     

基于微分响应的流域产流分单元修正方法

为考虑洪水预报误差的空间变化,提出一种基于微分响应的流域产流分单元修正方法.该方法建立了各单元流域产流与流域出口流量之间的微分响应关系,采用正则化最小二乘法结合逐步迫近进行反演求解,将产流误差估计量分配给相应单元流域实现流域产流分单元修正.将构建的方法应用于大坡岭流域和七里街流域进行新安江模型产流修正,比较分析了流域产流分单元修正、流域面平均产流修正和自回归修正的效果.结果表明：流域产流分单元修正效果优于流域面平均产流修正;随着预见期的增大,产流微分响应修正效果优于自回归修正.该方法通过汇流系统将流域出口断面流量信息进行分解用于修正各单元流域产流,有利于提高实时洪水预报精度.     

云南洱海桃溪河口净化工程的设计思路及初步净化效果

以太湖上游西苕溪流域为研究对象,通过分布式水文模型HEC-HMS模拟次降雨洪水过程:采用可视化数据存储系统HEC-DSS建立水文气象数据库,利用Geodatabase地理数据库技术集成流域自然属性数据库,通过距离平方倒数法对雨量数据进行空间插值,SCS曲线数法计算水文损失,运动波法计算直接径流与河道洪水演进,选用基流指数退水法模拟流域基流,并对模型中水库模拟部分进行适当修正.经模型校验,模拟结果表明,计算流量与观测流量拟合较好,效率系数大于0.8,洪峰流量误差低于4%,峰现时间误差低于2 h,该模型在土地利用变化对洪水水文要素的影响研究方面有较好的应用前景.     

降雨输入对实时洪水预报精度与预见期的影响

针对降雨输入不确定性对实时洪水预报影响的问题,本文采用不考虑未来预报降雨、考虑未来预报降雨、考虑预报降雨的降雨量误差和降雨时间误差4种方法,以陕西省两个半湿润流域（陈河流域和大河坝流域）为研究区域,分析不同预见期和不同降雨输入情况下洪水预报的精度.研究表明：相对于不考虑未来降雨情况,考虑未来降雨后在预报预见期较长时对预报结果精度提升较大,在预见期较短时对预报结果精度提升不显著;暴雨中心位置不同对预报精度影响也不同,当暴雨中心位于流域下游时降雨量误差对流量预报误差影响更大;降雨量误差主要影响洪量相对误差和洪峰相对误差,且这种影响是线性的,对确定性系数的影响是非线性的二次函数,降雨时间误差主要影响峰现时间误差.     

大尺度水文模型TOPX构建及其与区域环境系统集成模式RIEMS的耦合

基于改进型SIMTOP参数化径流方案和新安江模型的三层土壤水量平衡计算方法,本文构建了一个输入数据和率定参数较少、同时具有地形指数尺度转换机制、较好描述二维水文过程的简单高效的大尺度水文模型TOPX,并将其与区域环境系统集成模式RIEMS紧密耦合,以增强区域气候模式对大尺度流域径流量的定量数值模拟能力．TOPX模型在酉水河流域和泾河流域的离线测试表明：该模型对小尺度流域的径流量模拟精度较高,能够较好地描述流域水文变化过程;同时,该模型在大尺度上具有较强的分布式模拟能力,能够捕捉陆面水文过程的主要特征和时空演变特点．TOPX与RIEMS的耦合模式在泾河流域进行了在线测试,借助TOPX模型中的地形指数降尺度转换和水文过程产汇流机制,耦合模式实现了利用区域气候模式模拟的气象资料来驱动水文模型进行大尺度流域日径流量的模拟．进一步分析还表明：区域气候模式RIEMS模拟的降水时空分布数据的精度是影响耦合模式对径流量模拟效果的关键因素．     

基于四维变分同化法的巢湖流域南淝河水质模拟

针对河流模拟中未知不确定性源对模拟精度的影响,以巢湖流域南淝河为研究对象,建立了基于四维变分同化方法的南淝河干流水质模型,研究了含未知污染源的南淝河水质过程模拟.模型以未知污染负荷的动态变化过程为控制变量,通过同化沿河不同断面的逐日水质监测数据,识别不同河段的逐日入河污染负荷过程来实现水质过程的模拟,改变了常规模型模拟需提前预知并输入污染负荷的应用前提.模拟结果表明,采用四维变分同化方法的水质模拟结果有明显改进,重点河段水质模拟的纳什效率系数从小于0提高到0.5以上.识别的入河污染过程与降雨过程波动总体一致,证实南淝河的入河污染与降雨过程密切;同时,模型也可识别异常的入河负荷,提高模型对水环境问题的诊断分析能力.该方法可推广应用于复杂河流系统,为巢湖等流域污染来源定量解析、水质预测预警及污染管控提供支持.     

A semi-physical sediment yield model for estimation of suspended sediment in loess region

Sediment yield is a complex function of many environmental factors including climate,hydrology,vegetation,basin topography,soil types,and land cover.We present a new semi-physical watershed sediment yield model for the estimation of suspended sediment in loess region.This model is composed by three modules in slope,gully,and stream phases.For slope sediment yield,a balance equation is established based on the concept of hydraulic erosion capacity and soil erosion resistance capacity.According to the statistical analysis of watershed characteristics,we use an exponential curve to approximately describe the spatial variability of watershed soil erosion resistance capacity.In gully phase,the relationship between gully sediment concentration and flow velocity is established based on the Bagnold'stream power function.In the stream phase,we assume a linear dependence of the sediment volume in the reach on the weighted sediment input and output.The proposed sediment yield model is operated in conjunction with a conceptual hydrologic model,and is tested over 16 regions including testing grounds,and small,medium and large watersheds in the loess plateau region in the mid-reach of Yellow River.Our results indicate that the model is reasonable in structure and is able to provide a good simulation of sediment generation and transportation processes at both flood event scale and inter-annual time scale.The proposed model is generally applicable to the watersheds with soil texture similar to that of the loess plateau region in the Yellow River basin in China.     

Evaluation of an interrill soil erosion model using laboratory catchment data

Physically based soil erosion simulation models require input parameters of soil detachment and sediment transport owing to the action and interactions of both raindrops and overland flow. A simple interrill soil water transport model is applied to a laboratory catchment to investigate the application of raindrop detachment and transport in interrill areas explicitly. A controlled laboratory rainfall simulation study with slope length simulation by flow addition was used to assess the raindrop detachment and transport of detached soil by overland flow in interrill areas. Artificial rainfall of moderate to high intensity was used to simulate intense rain storms. However, experiments were restricted to conditions where rilling and channelling did not occur and where overland flow covered most of the surface. A simple equation with a rainfall intensity term for raindrop detachment, and a simple sediment transport equation with unit discharge and a slope term were found to be applicable to the situation where clear water is added at the upper end of a small plot to simulate increased slope length. The proposed generic relationships can be used to predict raindrop detachment and the sediment transport capacity of interrill flow and can therefore contribute to the development of physically‐based erosion models. Copyright © 1999 John Wiley & Sons, Ltd.     

Uncertainty based modeling of rainfall-runoff: Combined differential evolution adaptive Metropolis (DREAM) and K-means clustering

Simulation of rainfall-runoff process in urban areas is of great importance considering the consequences and damages of extreme runoff events and floods. The first issue in flood hazard analysis is rainfall simulation. Large scale climate signals have been proved to be effective in rainfall simulation and prediction. In this study, an integrated scheme is developed for rainfall-runoff modeling considering different sources of uncertainty. This scheme includes three main steps of rainfall forecasting, rainfall-runoff simulation and future runoff prediction. In the first step, data driven models are developed and used to forecast rainfall using large scale climate signals as rainfall predictors. Due to high effect of different sources of uncertainty on the output of hydrologic models, in the second step uncertainty associated with input data, model parameters and model structure is incorporated in rainfall-runoff modeling and simulation. Three rainfall-runoff simulation models are developed for consideration of model conceptual (structural) uncertainty in real time runoff forecasting. To analyze the uncertainty of the model structure, streamflows generated by alternative rainfall-runoff models are combined, through developing a weighting method based on K-means clustering. Model parameters and input uncertainty are investigated using an adaptive Markov Chain Monte Carlo method. Finally, calibrated rainfall-runoff models are driven using the forecasted rainfall to predict future runoff for the watershed. The proposed scheme is employed in the case study of the Bronx River watershed, New York City. Results of uncertainty analysis of rainfall-runoff modeling reveal that simultaneous estimation of model parameters and input uncertainty significantly changes the probability distribution of the model parameters. It is also observed that by combining the outputs of the hydrological models using the proposed clustering scheme, the accuracy of runoff simulation in the watershed is remarkably improved up to 50% in comparison to the simulations by the individual models. Results indicate that the developed methodology not only provides reliable tools for rainfall and runoff modeling, but also adequate time for incorporating required mitigation measures in dealing with potentially extreme runoff events and flood hazard. Results of this study can be used in identification of the main factors affecting flood hazard analysis.     

A three dimensional model for water-sediment movement

A three-dimensional k-ε-Ap two-fluid turbulence model is proposed to study liquid-particle two-phase flow and bed deformation.By solving coupled liquid-phase and solid-phase governing equations in a finite-volume method,the model can calculate the movement of both water and sediment.The model was validated by water-sediment transport in a 180° channel bend with a movable bed.The validation concerns two-phase time-averaged velocities,bed deformation,water depth,depth-averaged streamwise velocity,cross-stream bed profiles,and two-phase secondary flow velocity vectors.The agreement between numerical results and experimental results was generally good.The comparisons of the numerical results of different models show that the three-dimensional k-ε-Ap two-fluid turbulence model has a relatively higher accuracy than one-fluid model.     

Base flow separation for soil erosion simulation in a granitic forested headwater catchment using a process-based model,GeoWEPP

Distributed erosion models, which simulate the physical processes of water flow and soil erosion, are effective for predicting soil erosion in forested catchments. Although subsurface flow through multiple pathways is dominant for runoff generation in forested headwater catchments, the process-based erosion model, Geo-spatial interface for Water Erosion Prediction Project(Geo WEPP), does not have an adequate subsurface component for the simulation of hillslope water flow. In the current study, t...     

The effect of soil type,meteorological forcing and slope gradient on the simulation of internal erosion processes at the local scale

Numerical simulation experiments of water erosion at the local scale (20 × 5 m) using a process‐based model [Plot Soil Erosion Model_2D (PSEM_2D)] were carried out to test the effects of various environmental factors (soil type, meteorological forcing and slope gradient) on the runoff and erosion response and to determine the dominant processes that control the sediment yield at various slope lengths. The selected environmental factors corresponded to conditions for which the model had been fully tested beforehand. The use of a Green and Ampt model for infiltration explained the dominant role played by rainfall intensity in the runoff response. Sediment yield at the outlet of the simulated area was correlated positively with rainfall intensity and slope gradient, but was less sensitive to soil type. The relationship between sediment yield (soil loss per unit area) and slope length was greatly influenced by all environmental factors, but there was a general tendency towards higher sediment yield when the slope was longer. Contribution of rainfall erosion to gross erosion was dominant for all surfaces with slope lengths ranging from 4 to 20 m. The highest sediment yields corresponded to cases where flow erosion was activated. An increase in slope gradient resulted in flow detachment starting upstream. Sediment exported at the outlet of the simulated area came predominantly from the zone located near the outlet. The microrelief helped in the development of a rill network that controlled both the ratio between rainfall and flow erosion and the relationship between sediment yield and slope length. Copyright © 2010 John Wiley & Sons, Ltd.     

Mineral composition and particle size distribution of river sediment and loess in the middle and lower Yellow River

Hyperconcentrated flows often occur in the middle and lower Yellow River(MLYR)and its tributaries,within which the main sediment source originates from the Loess Plateau of China due to serious water erosion.Little is known about the properties of river sediment that is transported by hyperconcentrated flows,particularly with respect to the mineral composition and size distribution.Samples of sediment and loess were collected in the northern,middle,and southern Loess Plateau and the mainstream and tributaries of the MLYR.A total of 18 loess samples and 24 river sediment samples were analyzed to determine their sediment size distribution and mineral composition.The bottom loess samples reflected the original sedimentary features of the Loess Plateau,and the median particle size reduced,and the clay content increased from the north to the south of the study region.The surface loess has been weathered under the action of wind and rainfall,and the clay particle content in the surface loess samples was higher than that in the undisturbed bottom loess.Erosion of the surface soil due to rainfall and surface runoff means that fine particles(mostly clay)have been washed away.The median diameter of surface loess particles was a little larger than that of the bottom loess particles where water erosion dominates.The particle size became coarser with increasing distance from the estuary in the MLYR,which reflects depositional sorting in the river channel.Significant logarithmic relations were found between the median diameter of the sediment particles and the i)non-clay mineral content and ii)clay mineral content.Thus,clay and non-clay mineral compositions can be conveniently estimated from the particle size distribution.     

Piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events on the Loess Plateau,China

Establishing a universal watershed‐scale erosion and sediment yield prediction model represents a frontier field in erosion and soil/water conservation. The research presented here was conducted on the Chabagou watershed, which is located in the first sub‐region of the hill‐gully area of the Loess Plateau, China. A back‐propagation artificial neural model for watershed‐scale erosion and sediment yield was established, with the accuracy of the model, then compared with that of multiple linear regression. The sensitivity degree of various factors to erosion and sediment yield was quantitatively analysed using the default factor test. On the basis of the sensitive factors and the fractal information dimension, the piecewise prediction model for erosion and sediment yield of individual rainfall events was established and further verified. The results revealed the back‐propagation artificial neural network model to perform better than the multiple linear regression model in terms of predicting the erosion modulus, with the former able to effectively characterize dynamic changes in sediment yield under comprehensive factor conditions. The sensitivity of runoff erosion power and runoff depth to the erosion and sediment yield associated with individual rainfall events was found to be related to the complexity of surface topography. The characteristics of such a hydrological response are thus closely related to topography. When the fractal information dimension is greater than the topographic threshold, the accuracy of prediction using runoff erosion power is higher than that of using runoff depth. In contrast, when the fractal information dimension is smaller than the topographic threshold, the accuracy of prediction using runoff depth is higher than that of using runoff erosion power. The developed piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events, which introduces runoff erosion power and runoff depth using the fractal information dimension as a boundary, can be considered feasible and reliable and has a high prediction accuracy. Copyright © 2013 John Wiley & Sons, Ltd.     

Inverse Modeling Approach to Allogenic Karst System Characterization

Allogenic karst systems function in a particular way that is influenced by the type of water infiltrating through river water losses, by karstification processes, and by water quality. Management of this system requires a good knowledge of its structure and functioning, for which a new methodology based on an inverse modeling approach appears to be well suited. This approach requires both spring and river inflow discharge measurements and a continuous record of chemical parameters in the river and at the spring. The inverse model calculates unit hydrographs and the impulse responses of fluxes from rainfall hydraulic head at the spring or rainfall flux data, the purpose of which is hydrograph separation. Hydrograph reconstruction is done using rainfall and river inflow data as model input and enables definition at each time step of the ratio of each component. Using chemical data, representing event and pre-event water, as input, it is possible to determine the origin of spring water (either fast flow through the epikarstic zone or slow flow through the saturated zone). This study made it possible to improve a conceptual model of allogenic karst system functioning. The methodology is used to study the Bas-Agly and the Cent Font karst systems, two allogenic karst systems in Southern France.