床面附近泥沙交换率在悬移质输沙计算中的应用

总结并分析了关于悬移质运动方程的源/汇项或底部边界条件的已有的处理方法.认为在悬移质输沙计算中,这些方法因引入了平衡输沙概念、忽略了推移质颗粒对悬移质运动的影响或因缺乏床面附近泥沙质量交换率的理论表达式,均存在着某些不足.运用床面附近泥沙质量交换率理沦,对一般的悬移质输沙方程的底部边界条件进行了新的探讨,并提供了关于一维和垂向二维悬移质输沙的计算实例.     

基于B超成像的泥沙起动流速测量方法

床面上泥沙颗粒的起动测量一直是泥沙运动规律研究的难题。为此,提出了一种基于B超成像技术的泥沙起动流速测量方法。在模型试验水槽中,该方法利用B超仪获取水下地形及其附近粒子的运行图像,通过图像处理技术分析统计床面附近运动粒子的成像光斑个数,并分析其与流速之间的关系。结果表明,B超成像光斑个数在泥沙起动过程中存在一个突变过程,且该突变过程与泥沙起动运动相对应,可用来判定泥沙起动及其对应的起动流速,并利用该流速下的床面地形变化验证了该方法的正确性。该方法具有无扰动无干扰、适合于清水和浑水、易于实现自动化测量的特点。     

粘性泥沙分层运动特征的试验研究

应用水槽试验研究了粘性泥沙的运动特性及水流条件和泥沙浓度对粘性泥沙运动的影响。实验结果表明,在一定的水流和泥沙浓度条件下,粘性泥沙运动存在着特殊运动形式——"分层运动","分层运动"形成时水槽中泥沙浓度沿垂线分布成"阶梯"形状,水流流速沿垂线分布为中部最大,向水面和水槽底部方向流速减小。分层运动的形成及其特征与水流条件泥沙浓度等因素密切相关,不同的泥沙浓度下形成分层运动的水动力条件不同,泥沙浓度增加,形成分层运动的水流强度也相应增加。     

泥沙知识痤讲 (三)推移质运动

河流泥沙运动的主要形式有沿床面滚动、跳跃和悬浮。此外在个别条件下可能发生滑动(当床面光滑时)和床面成层移动(当水流强度很大时)。悬浮运动的泥沙称为悬移质;滚动、滑动、跳跃、层移的泥沙称为推移质。推移质运动是由底部水动力和颗粒碰撞所引起,只发生在床面附近的区域(底层)内。1.泥沙起动当水流速度由小增大到一定数值后,床面上原来静止的个别泥沙,或因其粒径较细,或由于其位置突出,再碰上较大的纵向瞬时流速,就开始滚动。此时这种泥沙处于起动状态;而相应的水流时均速度则称为该类泥沙的起动流速。泥沙起动流速关系到开始测验推移质的时机,不冲渠道的设计、河床冲刷的极限深度,以及截流和抢险抛石的尺寸等等,故有重要的实用意义。泥沙起动时所受的作用力有:水流正面推力 P=     

冰盖下冰花浓度分布的数值模拟

封冻河道内的冰花分布规律对冰期水文及环境分析有着重要的意义. 视冰花为连续介质, 基于多相流理论和欧拉-欧拉模型方法, 建立了河渠垂向二维标准k-ε紊流数值模型. 在不同的水流流速、 冰花颗粒尺寸、 上游来冰量条件下, 对冰盖下冰花浓度场进行了模拟分析, 得到了冰盖下冰花分布特征. 结果表明: 水流流速越大, 冰花沿垂向分布范围越大, 冰花浓度分布曲线斜率越大, 冰花浓度沿垂向的均匀性越大, 同时浓度最大值越小; 冰花颗粒尺寸越大, 冰花沿垂向分布范围越大, 冰花浓度分布曲线斜率越大; 上游来冰量越大, 冰花沿垂向分布范围越大, 在同一水深其浓度值也越大, 冰花浓度沿垂向分布的均匀性越大.     

水生植物对水沙运动影响的三维湍流模型

水生植物的存在改变了明渠内水流的流动结构,也影响着泥沙的输运。建立三维湍流模型,在水流控制方程中加入植被阻力项和植被密度项来考虑刚性植物对水动力特性和泥沙输运特性的影响。应用该三维数值模型计算了矩形水槽内淹没植被对水流水平时均流速垂向分布的影响、复式明渠边滩栽种挺水植被对水流深度、平均流速分布以及植被岛周围泥沙床面的冲淤变形的影响。数值计算结果与实测结果吻合良好,表明本模型可以有效地描述刚性水生植物对水流泥沙运动的影响。     

泥沙起动判别标准探讨

分析了现有起动标准的不足之处。在制定起动标准时,不但要考虑水流条件,还必须考虑泥沙颗粒在床面的相对位置。对同一粒径的泥沙颗粒,其起动条件不是一个常量,而是位于一个区间,它随颗粒在床面的相对暴露度而变化。根据力学原理,推导了泥沙颗粒起动时的临界无量纲切应力公式,并对其系数取值进行了理论计算。计算结果表明,对同一粒径、在同一起动状态下,泥沙颗粒的起动临界条件并非一个常量。并从理论上分析得到了临界值的范围,对个别及少量起动,其无量纲临界切应力为0.021～0.042及0.041～0.062。     

关于泥沙运动基本概率的研究

明渠挟沙水流运动具有两重性,既有确定性的一面,也有随机性的一面。分析了明渠水流运动、床面泥沙颗粒分布和水流中运动泥沙的随机性及其表达方式,指出泥沙运动是由静止、滚动、跃移和悬移４种状态组成,每种状态之间是不断转化的,可由６个转移概率来描述,并可进一步简化为由４个基本概率来表达,从理论上推导得到了泥沙运动基本概率计算方法;并采用实测资料进行了验证。这一研究成果在理论上和生产实际中都具有重要的意义。     

泥沙影响流速分布规律的试验研究

挟沙水流中泥沙颗粒的存在将影响水流的紊动结构,并进一步影响水流流速分布。采用超声流速仪MicroADV测量挟沙水流的时均流速沿垂线分布,分析了在不同水力条件、不同泥沙条件(包括泥沙浓度及颗粒粒径)下泥沙与水流相互作用对流速沿垂线分布的影响,在此基础上,考虑到现有流速分布公式在实际河道水流流速沿垂线分布中的应用,提出新的指数型流速分布公式形式;并根据实验资料,初步检验了该公式的精度,同时分析了指数型流速公式中系数k和指数m随泥沙浓度s、粒径d以及水力条件的变化规律。结果表明:系数k随泥沙浓度的增加而增大,同时水力条件以及泥沙颗粒粒径同样也影响着k的变化;而指数m随粒径d的增大而减小,随泥沙浓度s的增加而增大。     

非定常不稳定液固两相流动中旋涡对颗粒运动影响的数值研究

构建起双向耦合的液固两相流动旋涡动力学模型与数值方法;应用离散涡方法,计算非定常不稳定水流场;采用Lagrange方法模拟颗粒运动,颗粒对流体的反作用通过修正涡泡运动速度来实现。利用所建模型,计算了两种St数的泥沙粒子在圆柱绕流场中的运动。结果证明了液固两相流动中颗粒运动与旋涡存在着明确的相关结构:(1)当水沙混合物中的泥沙颗粒碰上旋涡时,泥沙颗粒被卷入旋涡中,被卷入旋涡中的泥沙颗粒在运动过程中始终分布于旋涡区;(2)均匀水沙混合物绕圆柱流动,由于流体流过圆柱时产生剧烈分离流动,使得在尾迹流内中等St数 (St～o (1))的泥沙颗粒从均匀水沙混合物中分离出来而往旋涡区聚集。     

Threshold for particle entrainment into suspension

ABSTRACT Laboratory observations regarding the limit conditions for particle entrainment into suspension are presented. A high‐speed video system was used to investigate conditions for the entrainment of sediment particles and glass beads lying over a smooth boundary as well as over a rough bed. The results extend experimental conditions of previous studies towards finer particle sizes. A criterion for the limit of entrainment into suspension is proposed which is a function of the ratio between the flow shear velocity and particle settling velocity. Observations indicate that particles totally immersed within the viscous sublayer can be entrained into suspension by the flow, which contradicts the conclusions of previous researchers. A theoretical analysis of the entrainment process within the viscous sublayer, based on force–balance considerations, is used to show that this phenomenon is related to turbulent flow events of high instantaneous values of the Reynolds stress, in agreement with previous observations. In the case of experiments with a rough bed, a hiding effect was observed, which tends to preclude the entrainment of particles finer than the roughness elements. This implies that, as the ratio between particle and roughness element sizes becomes smaller, progressively higher bed shear stresses are required to entrain particles into suspension. On the other hand, an overexposure effect was also observed, which indicates that a particle moving on a smooth bed is more prone to be entrained than the same particle moving on a bed formed by identical particles.     

Sensitivity of incipient particle motion to fluid flow penetration depth within a packed bed

A discrete element method is applied to a three‐dimensional analysis related to sediment entrainment on a micro‐scale. Sediment entrainment is the process by which a fluid medium accelerates particles from rest and advects them upward until they are either transported as bedload or suspended by the flow. Modelling of the entrainment process is a critically important aspect for studies of erosion, pollutant resuspension and transport, and formation of bedforms in environmental flows. Previous discrete element method studies of sediment entrainment have assumed the flow within the particle bed to be negligible and have only allowed for the motion of the topmost particles. At the same time, micro‐scale experimental studies indicate that there is a small slip of the fluid flow at the top of the bed, indicating the presence of non‐vanishing fluid velocity within the topmost bed layers. The current study demonstrates that the onset of particle incipient motion, which immediately precedes particle entrainment, is highly sensitive to this small fluid flow within the topmost bed layers. Using an exponential decay profile for the inner‐bed fluid flow, the discrete element method calculations are repeated with different fluid penetration depths within the bed for several small particle Reynolds numbers. For cases with slip velocity corresponding to that observed in previous experiments with natural sediment, the predicted particle velocity is found to be a few percent of the fluid velocity at the top of the viscous wall layer, which is a reasonable range of velocities for observation of incipient particle motion. This method for prescribing the fluid flow within the particle bed allows for the current discrete element method to be extended in future studies to the analysis of sediment entrainment under the influence of events such as turbulent bursting. Additionally, predictions for the slip velocities and fluid flow profile within the bed suggest the need for further experimental studies to provide the data necessary for additional improvement of the discrete element method models.     

Bed sediment entrainment by streamflow: State of the science

In fluvial sedimentology, bed sediment entrainment by streamflow has a decisive role in controlling several fluvial processes. Owing to its huge practical importance, the subject has been painstakingly explored for over a century. However, a detailed understanding of the mechanism of the bed sediment entrainment phenomenon achieved hitherto is far from complete. The central theme of bed sediment entrainment is occupied by the sediment entrainment threshold, which varies enormously in its qualitative definition, identification and quantification encompassing a broad range of spatiotemporal scales. This article presents the state of the science of the entrainment of non‐cohesive bed sediments under a steady‐unidirectional streamflow. It begins with the diverse definitions and representations of the entrainment threshold criterion from both qualitative and quantitative perspectives, scrutinising its suitability and ambit of applicability. Then, the effects of energetic factors that drive the entrainment threshold criterion are critically appraised. The indispensable mechanisms of bed sediment entrainment, including the theoretical background and modelling strategies, the role of turbulent bursting phenomenon and the phenomenological perspective into the origin of the scaling laws of sediment entrainment, are explained. Throughout the article, special emphasis is given to the strengths and weaknesses of the current state of the science. In addition, a deliberate attempt is made to invoke the thought‐provoking ideas on the multifarious features of bed sediment entrainment. Finally, the innovative perspectives on the bed sediment entrainment are provided and the concluding remarks are made, elucidating the major challenges and suggesting the prospective ways to resolve them as a future scope of research.     

The grain–fluid interaction as a self-stabilizing mechanism in fluvial bed load transport

A grain-scale model of fluvial bed load transport is described, with particular emphasis on the equilibrium between the saltating grains and the near bed flow, and its role in determining transport rate. The model calculates, explicitly, the modification of the velocity profile by the moving grains, together with the consequential reduction in surface fluid shear stress. As the surface fluid shear stress is reduced by the moving grains, so the entrainment rate decreases and the model reaches a steady state. The results provide insight into two important questions at a macroscopic level. First, they show that, in the absence of large static roughness, the dynamic roughness caused by the moving grains may be a significant contributor to flow resistance. Secondly, the model indicates the manner in which transport may be limited by a combination of the transport capacity of the flow and the availability of sediment for entrainment. Only in the case of high sediment availability does the fluid shear stress acting at the surface approach the critical entrainment value, reproducing the behaviour suggested by Bagnold (1956 ) and Owen (1964 ). This suggests that prediction formulae based on this assumption only describe the bed load transport system under particular conditions.     

泥沙起动概率的高速摄影测量方法

为了通过实验测量泥沙起动概率,基于封闭槽道内开展的泥沙起动实验,系统论述了基于高速摄影的起动概率测量方法所涉及的起动概率定义、相机采样频率、图像处理方法、起动比例与运动比例的关系等关键问题.根据泥沙在水流紊动猝发作用下起动的物理机制,定义实测起动概率为单个紊动猝发周期内,从一定面积床面上冲刷外移的泥沙比例.根据推移质运动概念模型,提出高速图片拍摄间隔应小于推移质运动的中间时间尺度,以捕获泥沙的每一次起动.在高速图像处理时,通过连续两帧图片之间的灰度差实现运动比例的无偏差提取.进一步分析表明,运动泥沙可处于起动、止动或滑翔状态,运动比例等于起动比例、止动比例和滑翔比例之和,起动比例与运动比例之比等于两张图片之间的时间间隔除以推移质运动的中间时间尺度.基于上述关键技术,将高速摄影测量方法用于泥沙起动概率的测量,结果表明,临界起动水流条件时,床沙处于少量起动现象,对应的临界起动概率与窦国仁所推导的2.28%接近.     

The physical character of subaqueous sedimentary density flows and their deposits

The complexity of flow and wide variety of depositional processes operating in subaqueous density flows, combined with post‐depositional consolidation and soft‐sediment deformation, often make it difficult to interpret the characteristics of the original flow from the sedimentary record. This has led to considerable confusion of nomenclature in the literature. This paper attempts to clarify this situation by presenting a simple classification of sedimentary density flows, based on physical flow properties and grain‐support mechanisms, and briefly discusses the likely characteristics of the deposited sediments. Cohesive flows are commonly referred to as debris flows and mud flows and defined on the basis of sediment characteristics. The boundary between cohesive and non‐cohesive density flows (frictional flows) is poorly constrained, but dimensionless numbers may be of use to define flow thresholds. Frictional flows include a continuous series from sediment slides to turbidity currents. Subdivision of these flows is made on the basis of the dominant particle‐support mechanisms, which include matrix strength (in cohesive flows), buoyancy, pore pressure, grain‐to‐grain interaction (causing dispersive pressure), Reynolds stresses (turbulence) and bed support (particles moved on the stationary bed). The dominant particle‐support mechanism depends upon flow conditions, particle concentration, grain‐size distribution and particle type. In hyperconcentrated density flows, very high sediment concentrations (>25 volume%) make particle interactions of major importance. The difference between hyperconcentrated density flows and cohesive flows is that the former are friction dominated. With decreasing sediment concentration, vertical particle sorting can result from differential settling, and flows in which this can occur are termed concentrated density flows. The boundary between hyperconcentrated and concentrated density flows is defined by a change in particle behaviour, such that denser or larger grains are no longer fully supported by grain interaction, thus allowing coarse‐grain tail (or dense‐grain tail) normal grading. The concentration at which this change occurs depends on particle size, sorting, composition and relative density, so that a single threshold concentration cannot be defined. Concentrated density flows may be highly erosive and subsequently deposit complete or incomplete Lowe and Bouma sequences. Conversely, hydroplaning at the base of debris flows, and possibly also in some hyperconcentrated flows, may reduce the fluid drag, thus allowing high flow velocities while preventing large‐scale erosion. Flows with concentrations <9% by volume are true turbidity flows (sensu 4 ), in which fluid turbulence is the main particle‐support mechanism. Turbidity flows and concentrated density flows can be subdivided on the basis of flow duration into instantaneous surges, longer duration surge‐like flows and quasi‐steady currents. Flow duration is shown to control the nature of the resulting deposits. Surge‐like turbidity currents tend to produce classical Bouma sequences, whose nature at any one site depends on factors such as flow size, sediment type and proximity to source. In contrast, quasi‐steady turbidity currents, generated by hyperpycnal river effluent, can deposit coarsening‐up units capped by fining‐up units (because of waxing and waning conditions respectively) and may also include thick units of uniform character (resulting from prolonged periods of near‐steady conditions). Any flow type may progressively change character along the transport path, with transformation primarily resulting from reductions in sediment concentration through progressive entrainment of surrounding fluid and/or sediment deposition. The rate of fluid entrainment, and consequently flow transformation, is dependent on factors including slope gradient, lateral confinement, bed roughness, flow thickness and water depth. Flows with high and low sediment concentrations may co‐exist in one transport event because of downflow transformations, flow stratification or shear layer development of the mixing interface with the overlying water (mixing cloud formation). Deposits of an individual flow event at one site may therefore form from a succession of different flow types, and this introduces considerable complexity into classifying the flow event or component flow types from the deposits.     

Effects of bed longitudinal inflexion and sediment porosity on basal entrainment mechanism: insights from laboratory debris flows

Channel morphology and bed sediment erodibility are two crucial factors that significantly affect debris flow entrainment processes. Current debris flow entrainment models mostly hypothesize the erodible beds are infinite with uniform slopes. In this study, a series of small-scale flume experiments were conducted to investigate the effects of bed longitudinal inflexion and sediment porosity on basal entrainment characteristics. Experimental observations revealed that sediment entrainment is negligible at early stages and accelerates rapidly as several erosion points appear. Continual evolution of flow-bed interfaces changes interactions between debris flows and bed sediments, rendering the interfacial shear action involved into a mixed shear and frontal collisional action. Lower bed sediment porosity will change the spatial arrangement and orientation of particle mixture, strengthen the interlocking and anti-slide forces of adjacent sediment particles, and promote the formation of particle clusters, all of which will increase bed sediment resistance to erosion. By examining the post-experimental bed morphology, the slope-cutting amounts and topographic reliefs are determined to positively correlate with longitudinal transition angles. These high topographic reliefs may indicate the propensity of triangular slab erosion, rather than strip-shaped slab erosion, in non-uniform channels with relatively steep erodible beds. Empirical formulas are obtained that denote the relationships among bed sediment strength, channel curvature radius, and sediment porosity through a multi-parameter regression analysis. This study may aid in clarifying the complex coupling effects of spatial variations in debris flow dynamics as well as sediment erodibility and bed morphology in non-uniform channels with abundant seismic loose material.

     

推移质翻越低坝输移特性的试验研究

对于山区河流低坝而言,平时淤积在坝前的推移质粗沙可能会在洪水期集中翻越坝顶,形成高强度输沙。本文开展水槽试验,研究推移质粗沙自上游起动、推进、再翻越坝顶后向下游输移的过程,分析了输沙参数的变化特性及数理规律,描述了翻坝输沙模式及运动特征,揭示了输沙规律与河床形态之间的自然联系。取得如下认识:①输沙量随时间大致以幂函数规律增长。②低坝附近区域河床形态终将趋于稳定,上游和下游均形成相对稳定的曲面斜坡淤积体。③在不同的水流强度下推移质翻坝输移模式存在差异。对于一般水流强度工况,上游淤积体曲面斜坡表面泥沙颗粒以滚动或滑动模式起动,推移至接近坝顶位置时再跃移翻坝,后向下游输移;对于更高水流强度工况,后期的翻坝输沙模式可能发生显著转变,周期性边壁漩涡成为翻坝输沙的主要动力来源。     

推移质翻越低坝输移特性的试验研究

对于山区河流低坝而言，平时淤积在坝前的推移质粗沙可能会在洪水期集中翻越坝顶，形成高强度输沙。本文开展水槽试验，研究推移质粗沙自上游起动、推进、再翻越坝顶后向下游输移的过程，分析了输沙参数的变化特性及数理规律，描述了翻坝输沙模式及运动特征，揭示了输沙规律与河床形态之间的自然联系。取得如下认识：①输沙量随时间大致以幂函数规律增长。②低坝附近区域河床形态终将趋于稳定，上游和下游均形成相对稳定的曲面斜坡淤积体。③在不同的水流强度下推移质翻坝输移模式存在差异。对于一般水流强度工况，上游淤积体曲面斜坡表面泥沙颗粒以滚动或滑动模式起动，推移至接近坝顶位置时再跃移翻坝，后向下游输移；对于更高水流强度工况，后期的翻坝输沙模式可能发生显著转变，周期性边壁漩涡成为翻坝输沙的主要动力来源。     

波流边界层水沙运动数值模拟——II.泥沙运动模拟

为解析波流边界层内泥沙运动,建立了基于水动力-泥沙-床面互馈过程的波流边界层1DV泥沙数学模型,可用于模拟不同床面形态下粉沙-沙的含沙量过程。床面形态模块提供床面形态类型和相应参数;给出了平底和沙波床面粗糙高度和泥沙扩散系数的确定方法;采用了适宜粉沙及沙的制约沉速、底部参考浓度和起动剪切应力等公式;引入含沙量层化效应和制约沉降反映水动力与泥沙之间的相互影响。水槽试验资料验证表明,建立的模型较好地模拟了不同床面不同波流组合条件下的含沙量剖面。在此基础上,讨论了不同床面含沙量剖面模拟方法的差异,指出床面形态是决定含沙量变化的重要因素之一,仅通过改变床面粗糙高度不足以反映漩涡沙波床面的含沙量剖面特征。该模型可为研究波流边界层内泥沙运动和物质输运提供工具。