模拟降雨条件下坡面流水动力学特性研究

为探明降雨条件下沙黄土坡面水流水动力学特性,以流体力学和泥沙运动力学理论为依据,通过5个坡度和5个雨强组合条件下室内模拟降雨试验,系统研究了坡面水流水力参数的变化规律。结果表明,降雨条件下坡面薄层水流雷诺数均小于580,处于层流失稳区;水流流型随雨强和坡度的变化而发生转捩,当坡度较缓、雨强较小时,床面形态处于低能态区和过渡区,床面出现沙纹现象,水流宏观上呈缓流,反之,坡度较陡、雨强较大时,床面由沙纹和沙垄向动平床过渡,宏观上多呈急流;并根据薄层水流阻力组成特点,推导出沙黄土坡面薄层水流阻力计算公式。验证结果表明,该式误差较小,可为坡面侵蚀预报模型的构建提供参考。     

模拟植被覆盖条件下坡面流水动力学特性

利用人工模拟植被试验,系统研究了6个坡度和10个流量、5种覆盖度条件下坡面流水力参数水力要素关系及阻力的变化特征,以期揭示坡面植被水流阻力的内在规律性.结果表明:流态指数随植被覆盖度的增加而显著增加,相同植被覆盖度下,试验坡度对其影响并不显著.各试验工况下,水流流态主要分布在"虚拟层流区"和过渡流区,未达到紊流区.随试验坡度的增加,水流流态由缓流向急流区域延伸,而随覆盖度的增加,流态向相反方向延伸.植被覆盖条件下坡面水流阻力与雷诺数的变化规律并非呈单调增加或者单调递减的趋势,而与植被覆盖度紧密相关,当植被覆盖度较低时,阻力系数随雷诺数的增加而减小,而覆盖度较高时,阻力系数随雷诺数增加呈线性增加的趋势.     

Study of surface runoff using physical models

. The profiles of overland water flow as a function of space and time obtained by applying kinematic wave approximations combined with the Darcy-Weisbach resistance formula to laminar flow are presented. Rainfall-excess is assumed to remain constant during a certain period of time. Runoff from surfaces of constant slope, with uniform surface texture, and the effects of different parameters on overland flow have been studied. Comparisons of runoff using Darcy-Weisbach, Manning, and Chezy resistance formulas have been made. It was found that the lower the rainfall-excess rates, the longer the surface runoff starting time, peak time, and smaller the peak runoff value at any distance. It was also found that the overland flow increases rapidly to the peak value, followed by a rapid decline which begins at the moment the rainfall-excess ceases, and then approaches zero slowly. Comparison of the theoretical calculations of runoff with the measured data of the Los Angeles field tests on concrete-paved surface shows good agreement.     

输沙能力模式在坡面非恒定输沙过程中的应用

在坡面土壤侵蚀输沙计算中,可选用的输沙能力模式较多,但这些模式多运用于恒定输沙,在非恒定输沙中运用不多。建立了坡面土壤侵蚀非恒定输沙数学模型,模型中的输沙能力运用水流切应力、水流功率、单位水流功率3种模式,对模型过程采用有限差分格式离散求解。根据实测水沙资料进行模型参数率定,运用3种输沙能力模式于不同降雨强度、不同坡度的非恒定坡面输沙过程中。结果表明:在坡面非恒定输沙计算中,不同输沙能力模式对计算结果有明显影响,在降雨强度较小时,单位水流功率模式结果较其他两个模式为好,而在雨强较大时,切应力模式计算结果较好。     

人工降雨条件下冲沟型泥石流起动试验研究

下垫面以位于贡嘎山东坡的熊家沟为模型,开展了不同降雨强度条件下冲沟型泥石流起动的模拟试验,初步研究了冲沟型泥石流的形成机理和演化特征.试验研究表明:(1)在强降雨条件下,水体入渗速度、不同深度土体含水量变化与降雨强度呈反比例关系,降雨强度越大,越不利于水体入渗,而有利于坡面汇流、冲沟径流和下切侵蚀; (2)在强降雨和径流条件下,土体破坏方式、破坏程度以及泥石流形成机理表现出差异性.相对较小雨强降雨条件下,土体破坏方式以滑坡为主,泥石流形成模式表现为滑坡液化与转化起动,雨强较大降雨条件下,土体破坏方式以侵蚀垮塌为主,泥石流形成模式为洪流席卷垮塌体和沟床揭底; (3)起动试验中泥石流阵性特征明显.在强降雨条件下,雨强与泥石流的规模、黏度之间没有正相关性,雨强越大,泥石流黏度越小,试验中多出现的是高含砂洪流,而相对较小雨强作用下由土体液化转化形成的泥石流黏度较大.试验现象和结果与熊家沟泥石流起动、发生过程具有较高的一致性.     

不同坡度草地含沙水流水力学特性及其拦沙机理

参照黄土区侵蚀降雨和坡面片蚀产沙特征,采用恒流泥沙输送装置模拟坡面上方来水来沙,探讨不同坡度草地含沙水流的水力学特性及其对上方来沙的拦蓄机理。结果表明,草地坡面的水流弗劳德数随坡度增大而增加,而阻力系数与坡度呈反势。按明渠水流的一般标准,不同坡度草地水流均为层状缓流。草地坡面拦沙效应随坡度增大而减小,且径流前期的减沙作用较后期更为显著。不同坡度草地坡面的出流泥沙平均直径和大颗粒(>10μm)泥沙含量均显著小于上方来沙,这说明草地的拦沙效应主要体现在对大粒径泥沙的拦蓄上。     

坡地氮磷流失过程模拟

根据农田坡面氮、磷流失的主要过程,建立了基于次降雨事件的坡面氮、磷迁移模型。模型采用改进的Green-Ampt方程和运动波方程计算下渗与坡面流,土壤侵蚀采用修改的欧洲土壤侵蚀模型计算,用考虑了侵蚀影响与扩散作用的迁移模型计算坡面氮、磷迁移过程,采用一维对流扩散方程计算氮、磷在土壤中的迁移过程。利用室内人工降雨资料对模型进行了率定及验证。结果表明,模型可以较好地模拟地表径流中氮、磷浓度,率定期与验证期模型效率系数均在0.89以上。根据坡地氮、磷流失机理,分析了降雨强度与地表坡度对坡面径流中氮、磷浓度的影响。     

Experimental study of free-surface flow instability and bedforms in shallow flows

To test erosion and development of bed sculpture in unstable, pulsating shallow flow, three series of experiments were performed in a 10-m flume with flows of Reynolds Number (Re) up to 10³, Froude Number (Fr) up to 3.0, and slope up to 2%. In the first series of experiments (59 runs), the bed was rigid (glass flume-bottom), in the second (50 runs) and third (10 runs) series, planar beds were made from well-sorted and poorly sorted sand, respectively. Development of flow surface instability above the rigid bed showed the following stages: smooth surface; small wavelets; well-developed waves; and rough flow. Development of bedform morphology in the second series of tests showed the following sequence: smooth bed; flow-aligned ridge and trough pattern; combined ridge and rhomboid pattern; and rhomboid configuration. In presence of flow instability, ridge patterns form in flows with mean values indicative of laminar-subcritical regime, combined ridge and rhomboid patterns in laminar-supercritical regime, and rhomboid patterns in flows with mean values indicative of transitional turbulent-supercritical regime. The third series of experiments demonstrated transverse segregation of the bed material across the ridges.The results, in line with stability theory, show that laminar and transitional flows already become unstable when $\text{Fr}\text{=}\text{1}\text{2}$ and suggest that pulsations may appear in sheet-floods and overland flow, intensifying erosion and transport of sediment. The precise mechanism of origin of the parent kinematic flow configurations instrumental in generation of the ridge and the rhomboid patterns is not yet clear. It may be related to the mechanism operating in transition and generation of turbulence in parallel flows, or flows with a superimposed fluctuation, but other types of hydrodynamic instability are also possible.     

Depositional processes,bedform development and hybrid bed formation in rapidly decelerated cohesive (mud–sand) sediment flows

Flows with high suspended sediment concentrations are common in many sedimentary environments, and their flow properties may show a transitional behaviour between fully turbulent and quasi‐laminar plug flows. The characteristics of these transitional flows are known to be a function of both clay concentration and type, as well as the applied fluid stress, but so far the interaction of these transitional flows with a loose sediment bed has received little attention. Information on this type of interaction is essential for the recognition and prediction of sedimentary structures formed by cohesive transitional flows in, for example, fluvial, estuarine and deep‐marine deposits. This paper investigates the behaviour of rapidly decelerated to steady flows that contain a mixture of sand, silt and clay, and explores the effect of different clay (kaolin) concentrations on the dynamics of flow over a mobile bed, and the bedforms and stratification produced. Experiments were conducted in a recirculating slurry flume capable of transporting high clay concentrations. Ultrasonic Doppler velocity profiling was used to measure the flow velocity within these concentrated suspension flows. The development of current ripples under decelerated flows of differing kaolin concentration was documented and evolution of their height, wavelength and migration rate quantified. This work confirms past work over smooth, fixed beds which showed that, as clay concentration rises, a distinct sequence of flow types is generated: turbulent flow, turbulence‐enhanced transitional flow, lower transitional plug flow, upper transitional plug flow and a quasi‐laminar plug flow. Each of these flow types produces an initial flat bed upon rapid flow deceleration, followed by reworking of these deposits through the development of current ripples during the subsequent steady flow in turbulent flow, turbulence‐enhanced transitional flow and lower transitional plug flow. The initial flat beds are structureless, but have diagnostic textural properties, caused by differential settling of sand, silt and cohesive mud, which forms characteristic bipartite beds that initially consist of sand overlain by silt or clay. As clay concentration in the formative flow increases, ripples first increase in mean height and wavelength under turbulence‐enhanced transitional flow and lower transitional plug‐flow regimes, which is attributed to the additional turbulence generated under these flows that subsequently causes greater lee side erosion. As clay concentration increases further from a lower transitional plug flow, ripples cease to exist under the upper transitional plug flow and quasi‐laminar plug flow conditions investigated herein. This disappearance of ripples appears due to both turbulence suppression at higher clay concentrations, as well as the increasing shear strength of the bed sediment that becomes more difficult to erode as clay concentration increases. The stratification within the ripples formed after rapid deceleration of the transitional flows reflects the availability of sediment from the bipartite bed. The exact nature of the ripple cross‐stratification in these flows is a direct function of the duration of the formative flow and the texture of the initial flat bed, and ripples do not form in cohesive flows with a Reynolds number smaller than ca 12 000. Examples are given of how the unique properties of the current ripples and plane beds, developing below decelerated transitional flows, could aid in the interpretation of depositional processes in modern and ancient sediments. This interpretation includes a new model for hybrid beds that explains their formation in terms of a combination of vertical grain‐size segregation and longitudinal flow transformation.     

Application of physical modelling of debris flow triggering to field conditions: Limitations posed by boundary conditions

Debris flows are often triggered by Hortonian overland flow during high-intensity rainstorms. Data derived from debris flow trigger zones in the southern French Alps were fed into a physical model of debris flow triggering based on Takahashi. Using a Monte Carlo approach with 1000 runs, the results show a wide distribution of safety factor values, indicating that physical modelling based on actual field measurements may not always be practical.As all safety factor values obtained are well below 1 even though debris flows only occur during very high-intensity rainstorms, the model used must be inappropriate. Apparently, the composition of the overland flow plays an important role: during high-intensity rainstorms it usually has a very high sediment content and contains stones. This prevents it from flowing through the pores of coarse debris accumulations in the central gully of a trigger zone; it will rather run over the debris. This situation is more stable than with the fluid flowing through the pores. The behaviour switch of the fluid above a certain sediment and stone content thus drastically changes the triggering conditions for debris flows and it is concluded that debris flow triggering in the area requires the occurrence of both overland flow and landsliding.     

Influence of hydrate exploitation on stability of submarine slopes

Urban pluvial flash floods have become a matter of widespread concern, as they severely impact people’s lives in urban areas. Hydrological and hydraulic models have been widely used for urban flood management and urban planning. Traditionally, to reduce the complexity of urban flood modelling and simulations, simplification or generalization methods have been used; for example, some models focus on the simulation of overland water flow, and some models focus on the simulation of the water flow in sewer systems. However, the water flow of urban floods includes both overland flow and sewer system flow. The overland flow processes are impacted by many different geographical features in what is an extremely spatially heterogeneous environment. Therefore, this article is based on two widely used models (SWMM and ANUGA) that are coupled to develop a bi-directional method of simulating water flow processes in urban areas. The open source overland flow model uses the unstructured triangular as the spatial discretization scheme. The unstructured triangular-based hydraulic model can be better used to capture the spatial heterogeneity of the urban surfaces. So, the unstructured triangular-based model is an essential condition for heterogeneous feature-based urban flood simulation. The experiments indicate that the proposed coupled model in this article can accurately depict surface waterlogged areas and that the heterogeneous feature-based urban flood model can be used to determine different types of urban flow processes.

     

对坡面径流挟沙力研究的几点认识

坡面径流挟沙力是土壤侵蚀过程模型的核心输入参数。在系统界定坡面径流挟沙力概念的基础上,深入分析了坡面径流水动力学特性和泥沙特性对坡面径流挟沙力的影响及其机制,受降雨、地形、土壤、植被、土地利用、结皮形成、细沟发育和砾石出露等多种因素时空变异的综合影响,坡面径流挟沙力具有强烈的时空变化特征,导致坡面径流挟沙力的测定较为困难。将确定挟沙力常用的方法划分为河流挟沙力方程修正法、推理法、模型模拟法和测量法,比较分析了不同方法的优缺点。未来的研究应该重点关注挟沙力测定方法与技术、坡面径流输沙过程的动力机制、典型土壤的挟沙力、挟沙力方程构建与验证及变化条件下坡面径流挟沙力等方面。     

Geochemistry and Genesis of the Pampean and Post-Pampean Formations (Late Pleistocene-Holocene), Central and Northern Argentina Pampas and Its Adjacent Region-An Approaching Case Study

Littoral siliciclastic shallow marine horizontal conformable beds (a heterolithic succession), Middle to Late Miocene, outcropping in northeastern Argentina (Ituzaingó Formation) and overlying transitional conformable horizontal regolithic mantle-rock bed derived from them (the Pampean and Post-Pampean Formations), were geochemically analyzed. The focus of this study is placed on the application of geochemical parameter and signature analyses related with the aforementioned geological units, which are of subcontinental extension into South America. The encountered results show an outstandingly similar geochemical behaviour between them. The main conclusion is that regolithic mantle-rock beds were derived from the littoral shallow marine mudstone (silty-argillaceous) beds. This is in oposition to previous aeolian processes proposed early in the 50's and later. These mudstone beds constitute important sections of the littoral shallow marine sequence beds (outcropping HST parasequence). Such regolitization proposed for the Pampean and Post-Pampean Formations predominantly developed in-situ during the Late Pleistocene and Holocene. The main erosional and mobilized agents were the surface free water (pluvial, fluvial and laminar water sheets) and vadose water. So, the wind flows and/or the immense air flow hurricanes are of insignificant sedimentological influence. Otherwise, there is not a proved appropiate sandy-silty reservoir for a reasonable support of the "aeolian hypothesis", as well as clear aeolian structures settled in the regolithic mantle-rock bed. In spite of these lines of equality textural-structural evidence, the geochemical values for both major elements and trace elements, demonstrate that the Miocene heterolithic marine succession was the mother rock of the overlying regolithic mantle-rock bed, which was formed as an in-situ mantle-rock bed.     

黑河上游高寒山区降雨-径流形成过程的同位素示踪

为揭示中、低纬度高寒山区降雨-径流的形成过程,指导水资源的合理开发利用,选择黑河上游红泥沟小流域为研究区,基于河道径流量及雨水和河水稳定同位素的观测数据,构建二元混合模型,计算了2013年7月23日及8月21日两次典型降雨-径流事件中事件水（降雨）和事件前水（流域前期储水）对河道径流的贡献及其动态变化.结果显示:两次降雨事件中事件前水的贡献比例分别为68.69%和54.46%;事件前水的贡献比例在涨水阶段减小,在退水阶段增大.结合河水电导率的观测结果,进一步分析了降雨-径流的形成过程:河道径流的形成主要受饱和区蓄满产流、河岸带地下径流和山坡地下径流3种产流机制控制;事件水主要源于蓄满产流,事件前水主要源于河岸带和山坡地下径流;事件初期和末期以河岸带地下水补给为主,涨水阶段后期和退水阶段前期转为以蓄满产流和山坡地下水补给为主,洪峰期间蓄满产流的贡献达到最大.两次事件的对比表明,事件前的湿度条件和降雨强度对降雨-径流的形成过程有着重要影响:前期越湿润,流域储水能力越弱,导水能力越强,事件水的贡献越大,河道径流对降雨的响应越迅速;降雨强度越大,蓄满产流及其中的事件水比例越高,河道径流中事件水的比例也越高.      

A coupled surface/subsurface flow model accounting for air entrapment and air pressure counterflow

This work introduces the soil air system into integrated hydrology by simulating the flow processes and interactions of surface runoff, soil moisture and air in the shallow subsurface. The numerical model is formulated as a coupled system of partial differential equations for hydrostatic (diffusive wave) shallow flow and two-phase flow in a porous medium. The simultaneous mass transfer between the soil, overland, and atmosphere compartments is achieved by upgrading a fully established leakance concept for overland-soil liquid exchange to an air exchange flux between soil and atmosphere. In a new algorithm, leakances operate as a valve for gas pressure in a liquid-covered porous medium facilitating the simulation of air out-break events through the land surface. General criteria are stated to guarantee stability in a sequential iterative coupling algorithm and, in addition, for leakances to control the mass exchange between compartments. A benchmark test, which is based on a classic experimental data set on infiltration excess (Horton) overland flow, identified a feedback mechanism between surface runoff and soil air pressures. Our study suggests that air compression in soils amplifies surface runoff during high precipitation at specific sites, particularly in near-stream areas.     

坡面流层流区动力学特性

为完善坡面水流的基础理论研究,基于定床水槽试验,以甘油溶液为试验流体,采用超声波测量技术,进行了4种粗糙度、5种坡度和13种单宽流量条件下的组合试验,研究了坡面流层流区的动力学特性及滚波特征。研究结果表明:在雷诺数为8~160的范围内,流态指数在理论值0.33附近波动,且随粗糙度的增大呈现先增大后减小的变化趋势,在粗糙度为0.10 mm附近达到峰值;阻力系数受坡度和粗糙度影响显著,可较好地由坡度、粗糙度和雷诺数的关系式表示;随着雷诺数的增大,滚波波速和波峰均呈幂函数形式递增,滚波周期大小无明显变化;坡度的增大会使层流失稳临界单宽流量减小,粗糙度的增大会使临界弗劳德数的均值减小。     

流域产沙模型的研究

由坡面泥沙颗粒的动力平衡条件,推导出坡面上细颗粒泥沙的起动切应力,由坡面径流的剩余输沙能力,导得坡面径流的侵蚀量计算公式,与文献[4]的雨滴溅蚀量计算公式一起,构成了概念清楚、过程明确、结构完整的流域产沙模型。产沙量取决于水流挟沙能力与可供沙量的对比关系,若水流挟沙能力小于供沙量,则产沙量等于水流挟沙能力;反之,水流将进一步冲刷表土,形成径流侵蚀,产沙量就等于供沙量与径流侵蚀量之和。据此,可计算出坡面和流域产沙量。模型在岔巴沟流域上应用,结果令人满意。     

The stratigraphic record and processes of turbidity current transformation across deep‐marine lobes

Sedimentary facies in the distal parts of deep‐marine lobes can diverge significantly from those predicted by classical turbidite models, and sedimentological processes in these environments are poorly understood. This gap may be bridged using outcrop studies and theoretical models. In the Skoorsteenberg Formation (South Africa), a downstream transition from thickly bedded turbidite sandstones to argillaceous, internally layered hybrid beds, is observed. The hybrid beds have a characteristic stratigraphic and spatial distribution, being associated with bed successions which generally coarsen and thicken‐upward reflecting deposition on the fringes of lobes in a dominantly progradational system. Using a detailed characterization of bed types, including grain size, grain‐fabric and mineralogical analyses, a process‐model for flow evolution is developed. This is explored using a numerical suspension capacity model for radially spreading and decelerating turbidity currents. The new model shows how decelerating sediment suspensions can reach a critical suspension capacity threshold beyond which grains are not supported by fluid turbulence. Sand and silt particles, settling together with flocculated clay, may form low yield strength cohesive flows; development of these higher concentration lower boundary layer flows inhibits transfer of turbulent kinetic energy into the upper parts of the flow ultimately resulting in catastrophic loss of turbulence and collapse of the upper part of the flow. Advection distances of the now transitional to laminar flow are relatively long (several kilometres) suggesting relatively slow dewatering (several hours) of the low yield strength flows. The catastrophic loss of turbulence accounts for the presence of such beds in other fine‐grained systems without invoking external controls or large‐scale flow partitioning and also explains the abrupt pinch‐out of all divisions of these sandstones. Estimation of the point of flow transformation is a useful tool in the prediction of heterogeneity distribution in subsurface systems.     

A semi‐analytical solution of surface flow and subsurface flow on a slope land under a uniform rainfall excess

A new approach is proposed to analyze the surface flow and subsurface flow passing over a pervious ground under a uniform rainfall excess. The flow field is divided into two regions that are called water layer and soil layer. To figure out the hydraulic behavior of overland flow on an inclined plane under a rainfall event, the simplified Navier–Stokes equations are employed for the surface water flow, and the flow inside the soil layer is porous media flow, which is governed by Biot's (1956, 1962) theory of poroelasticity. The velocity distribution of overland flow is nonzero at the ground surface. The relation between water depth and slope length was developed first. The profile of surface water flow was then found backwards from the downstream end of the flow section by the Runge–Kutta method. After that, the flow velocity and flow discharge of each layer could also be obtained via the water depth. Finally, the variation of fluid shear stress inside the soil layer is also discussed. Copyright © 2011 John Wiley & Sons, Ltd.