粘土心墙坝漫顶溃坝过程离心模型试验与数值模拟

利用作者研制成功的溃坝离心模型试验系统,对粘土心墙坝漫顶溃决过程进行了试验研究,结果发现粘土心墙坝与均质坝溃决机理与溃口发展规律明显不同,随着漫坝水流对下游坝壳冲蚀程度的增加,粘土心墙发生剪断破坏,溃口洪水流量迅速增大.基于上述试验结果,提出了一个描述粘土心墙坝漫顶溃坝过程的数学模型,并建议了相应的数值计算方法.该模型...     

喜马拉雅中段冰湖终碛坝体溃决过程实验研究

终碛坝广泛分布于世界各地的高山和极高山区。为了探究终碛坝的溃决过程,了解溃口的演变特征,文章以嘉龙错终碛坝的原位实验,模拟了终碛湖漫顶溃决过程。通过分析实验结果发现:（1）根据观察,将终碛坝溃决过程划分为坝体下游坡面冲刷、“溯源侵蚀”、出水口下切和溃口拓宽四个阶段。（2）上游湖区崩塌体激发的涌浪会造成溃口内的瞬时流量增加数倍,从而使得在有涌浪和无涌浪的条件下,“溯源侵蚀”过程出现陡坎和斜坡两种下切型。（3）通过分析溃口下切侵蚀过程,发现溃口的下切侵蚀发展过程主要受到坝体孔隙比和细粒含量的影响,并且溃口中点侵蚀率与水流剪应力存在一定的线性关系,符合线性侵蚀模型。通过分析发现,嘉龙错终碛坝的侵蚀系数为0.051,临界启动应力为237.64 Pa。与堰塞坝相比,可侵蚀系数比更小,而临界启动应力更大。     

无粘性土堤漫顶溃决水流运动试验研究

堤防溃决时溃口附近水流形态及溃决流量变化过程的正确描述是指导溃口堵复及时避险、预测溃堤洪水演进并进一步评估淹没损失的前提条件.通过弯道水槽中的概化模型试验,根据无粘性土堤防漫顶溃决的溃口变化和水流运动特征,将无粘性土堤防溃决过程分为漫流、冲槽、展宽及稳定4个阶段.试验系统研究了堤防溃决后河道以及溃口附近水位的变化过程,研究了河道洪水流量、洪水位和筑堤材料对溃口处水位的影响.基于大量的试验资料,进一步归纳总结出溃决流量与溃口口门宽度及溃口处水头变化之间的关系.     

大型天然土石坝的溃坝方式及环境效应分析

岷江上游分布有大量的古滑坡堵江形成的土石坝体。本文首先调查了扣山古滑坡堵江坝形成的地质环境和坝体的几何特征、溃口形态;然后采取坝体土石样品,进行了室内大型土工试验,考虑湖水、坝体渗流和岩土体的耦合作用,采用有限元进行了古滑坡堵江坝的稳定性分析,计算了溃坝洪水的最大洪流量及其洪水在下游的演化特征,认为该天然滑坡堵江坝的溃决方式为漫坝瞬时全溃,天然土石坝的溃决将对沿岸及下游地区的生态环境产生巨大的影响。     

堰塞坝漫顶溃决过程溃口纵向下切与溯源

为进一步了解堰塞坝溃坝过程,开展了9组水槽模型试验,对溃口纵向下切和溯源发展过程进行了系统分析,并讨论了上溯源点移动速度与溃口水深之间的关系。研究发现：非黏性堰塞坝溃坝过程中,冲刷面与底床的夹角时刻发生变化,上、下溯源点位置不固定但也不能完全发展到坝踵;上、下坡面坡度增大到最大值1：1.5时,下溯源点到下游坝趾的最大距离与坝体沿水流方向长度的比值（x_p^*/x_d^*,反映下溯源点最终相对位置）对应降低到最小值0.24和0.18;坝体相对尺寸从1减小到1/2时,x_p^*/x_d^*值从0.38增大到0.47。上溯源点的无量纲移动速度是不断变化的,在无量纲时刻为0.13时,其x,y分量分别达到峰值0.94和0.32;上溯源点处溃口水深出现时刻相对移动速度峰值点出现时刻有延迟,大概延迟0.04个无量纲时间。     

土石坝溃决模拟及水流计算研究进展

对土石坝溃决模拟技术的最新研究成果进行了总结和评价,在此基础上对该研究领域今后的研究工作提出了若干建议,包括:应研究提出不同坝型溃决可能性分析评价方法;针对不同坝型开展溃坝离心模型试验,揭示在不同致灾因子作用下,土石坝溃决机理和溃口发展过程;开展土石坝初始管涌形成以后发展过程的试验研究,揭示孔流转变为堰流的控制条件;针对溃坝水流的流线曲率较大、溃坝水流大多是非恒定超临界流以及筑坝材料粒径级配范围变化大等特点,开展大型溃坝水工水力学模型试验,揭示不同坝型的溃口流量过程、泥沙输移及下游河道洪水演进规律。     

喜马拉雅中段冰湖终碛坝体溃决过程实验研究——以嘉龙错冰湖终碛坝为例

终碛坝广泛分布于世界各地的高山和极高山区。为了探究终碛坝的溃决过程,了解溃口的演变特征,文章以嘉龙错终碛坝的原位实验,模拟了终碛湖漫顶溃决过程。通过分析实验结果发现:(1)根据观察,将终碛坝溃决过程划分为坝体下游坡面冲刷、"溯源侵蚀"、出水口下切和溃口拓宽四个阶段。(2)上游湖区崩塌体激发的涌浪会造成溃口内的瞬时流量增加数倍,从而使得在有涌浪和无涌浪的条件下,"溯源侵蚀"过程出现陡坎和斜坡两种下切型。(3)通过分析溃口下切侵蚀过程,发现溃口的下切侵蚀发展过程主要受到坝体孔隙比和细粒含量的影响,并且溃口中点侵蚀率与水流剪应力存在一定的线性关系,符合线性侵蚀模型。通过分析发现,嘉龙错终碛坝的侵蚀系数为0.051,临界启动应力为237.64 Pa。与堰塞坝相比,可侵蚀系数比更小,而临界启动应力更大。     

堰塞坝溃坝模型实验研究综述

堰塞坝是天然形成的坝体,结构比较松散、稳定性差、渗透作用强,发生溃决危险性大、概率高、突发性强,而且破坏可能性高及产生的洪水威胁人们的生命和财产安全,因此需要系统、全面的研究。作者从单坝溃坝、级联溃坝及堰塞坝处置的溃坝试验进行详细的回顾,总结及分析了国内外学者在堰塞坝溃坝模型实验取得的成果及局限性,进一步分析了单坝溃决的颗粒级配、密实度、含水率、沟床坡度等因素,最后讨论了溃坝因素与溃决模式、溃决特征、溃决流量、溃口演化的关系。基于模型实验相似理论及模型比尺、实验测试手段、堰塞坝处置三个方面,提出了今后的研究重点。     

考虑河床坡度和泄流槽横断面影响的堰塞坝溃决过程试验研究

针对缺乏地形条件和工程处置措施对堰塞坝溃决过程影响研究的现状,采用4种河床坡度(0°、1°、2°、3°)和3种泄流槽横断面型式(三角形、梯形、复合型),开展了堰塞坝溃决的模型试验。通过分析堰塞坝的溃决流量、溃决历时、溃口发展和坝体纵截面演变过程,研究了不同河床坡度和泄流槽横断面对堰塞坝溃决过程的影响规律。试验结果表明:(1)堰塞坝溃决过程可分为3个阶段。阶段Ⅰ:溃口形成阶段,溃决流量较小;阶段Ⅱ:溃口发展阶段,水流下蚀及侧蚀强烈,溃决流量到达峰值;阶段Ⅲ:衰减-平衡阶段,粗化层形成,溃口停止发展。(2)河床坡度增加意味着下游坝坡、坝顶及泄流槽的坡度增加,导致水流侵蚀能力增强,溃口下切迅猛,因此在0°～3°范围内河床坡度越大,峰值流量越大,峰现时间越早,溃决流量过程曲线越趋于“高瘦型”,且残留坝高越小。(3)泄流槽横断面型式不同导致其槽深、槽宽和侧坡坡度不同,进而影响溃口发展和溃决流量。三角形槽的水土作用面积小,溃口下切及展宽速率最高,峰值流量最大,峰现时间最早;梯形槽的槽底高程最高,水土作用面积最大,溃口下切速率最低,峰现时间最晚;而复合槽介于前两者之间。试验成果将为堰塞坝应急抢险和工...     

精细地形下的尾矿坝稳定性及溃坝模拟分析

目前在尾矿坝稳定性和溃坝模拟分析方面,对溃口位置及水砂的流动状态难以做出准确判断。将尾矿坝稳定性和溃坝模拟有机结合,采用FLAC3D计算正常水位、洪水位、漫顶水位三种工况下尾矿坝稳定性,并利用Rhino与Fluent建立尾矿库及下游精细地形,开展尾矿库溃坝水砂在不同时刻及不同地形下的流动状态研究分析。结果表明:（1）浸润线的埋深随尾矿库水位的升高而变小,由正常水位升高至洪水位时浸润线埋深下降5～8 m,漫顶水位时坝顶浸润线沿坡面向下运移约8 m;（2）库水位对剪切带及尾矿坝稳定性有显著影响,从正常水位到洪水位时,剪切带纵向上不断向坝体内部延伸,横向上不断向坝脚延伸,剪切应变率增大为5.78×10?5,尾矿坝稳定系数由1.80下降至1.32;（3）达到坝顶时剪切带急剧缩短,而剪切应变率进一步增大为3.32×10?4,尾矿坝稳定系数由1.32下降到1.18。溃坝水砂的流动状态受地形影响明显,在山谷中表现为范围减小、流速增大的汇聚流动,在平坦农田处表现为范围增大、流速减小的发散流动特点。     

瞬时溃坝最大流量计算新通式推导及验证

为了计算瞬时溃坝最大流量,基于堰流与波流量相等原理,建立了瞬时溃坝最大流量与堰流关系,理论推导得到同一公式可以计算大坝全溃、横向局部溃坝、垂向局部溃坝及横垂向局部溃坝的瞬时最大流量的新通式,给出了该通式的基本和上限、下限及一般等通式。研究发现:瞬时溃坝流量与溃口堰坎类型有关,可以把瞬时溃坝最大流量、逐渐溃坝洪水或漫坝洪水的计算公式和相应系数统一到不同堰坎类型的堰流量计算中;而这些新通式相应系数的取值是被实验确认和经典著作肯定的,保证了系数取值的可靠和准确。这些新通式容易理解、形式简单、计算简便,参数取值可靠,把溃坝洪水计算的复杂问题转变为处理经典水力学的简单堰流计算。同时通过对比分析,充分论证了一般通式的合理性,并用国内外溃坝实例验证了计算方法和成果较科学、可靠、合理。     

明渠湍流涡运动尺度分布特性

为了解明渠湍流涡运动尺度沿水深的分布特征,使用高精度粒子图像测速法水槽试验系统测量充分发展的明渠湍流,利用改进后的功率谱分析方法,对明渠湍流涡运动尺度沿水深分布、上升流、下扫流和平均涡量等进行了分析。结果表明:明渠湍流涡运动尺度大小沿水深增加呈先增大后减小;沿水深增加下扫流比上升流更容易发生;各自水深的顺向涡涡量大于逆向涡涡量,在近壁区平均涡量最大;明渠湍流近壁区产生涡量大、尺度小的涡结构,沿水深增加,涡尺度发展变大,但受水面区下扫流抑制,涡尺度又变小。     

淹没植物条件下明渠KH涡的垂向几何尺度和频率

含淹没植物明渠中存在的Kelvin-Helmholtz（KH）涡,对水流能量交换和物质输移具有直接的影响。基于淹没刚性植物条件下的水槽实验,利用粒子图像测速技术（PIV）得到瞬时流速场和涡量场,定量刻画KH涡的垂向几何尺度（上下边界和垂向作用范围）,结合三维声学多普勒流速仪（ADV）测量的流速时间序列和小波变换方法,分析KH涡的频率,研究植物淹没度对KH涡上下边界位置、作用范围和频率的影响。研究发现:本实验确定的KH涡上、下边界与根据能量分布计算得到的混合层上、下边界基本一致;在本实验条件下,KH涡上边界为0.7~0.8倍水深,下边界在植物密度一定时围绕一常数波动,KH涡的频率在0.14~0.20 Hz之间;随着淹没度的增大,KH涡的上边界升高,下边界变化不明显,作用范围和频率增大。     

Bed forms in bimodal sand–gravel sediments: laboratory and field analysis

Bed forms were studied in Goodwin Creek and a laboratory flume channel. The bed sediment of the field site and flume had median diameters of 8·3 (modes of 0·4, 22·6 mm) and 1·82 mm (modes of 0·5, 5·6 mm), respectively. The laboratory and field channels had similar values of bimodal parameters, ratios of flow depth to median bed material diameter, and ratios of shear stress to critical shear stress and were judged to be comparable in the transport of bed load sediment and the resulting bed forms. Three groupings of bed forms from the laboratory flume experiments (ripple-like bed forms, bed load sheets, low-relief bed waves) were identified using the height and period of the bed forms. For the range of flow depths and discharges investigated in the flume, bed forms became higher and longer with increasing bed shear stress. Bed forms from Goodwin Creek were similar to those from the flume with comparable ratios between bed form length, height, and flow depth. The bed forms in the flume provide a positive link between rate and size fluctuations measured in the field and the bed forms. The smaller bed forms identified were sediment starved and are not considered to be dunes, while the largest bed forms in which all of the bed material sizes were mobilized in the field and laboratory were judged to be dunes.     

Effects of fluid resistance in the mine waste dam-break problem

A characteristic common to most mine tailings dam failures is that the mine tailings tend to liquefy and flow over substantial distances, with the potential for extensive damage to property and life. In order to be able to assess the potential for damage in case of such a failure, it is necessary to be able to predict the characteristics of the flow and the possible extent of flood movement. A procedure that has some potential for such flow analyses is described in this paper. When fluid resistanceis neglected, the solution of the simple dam-break problem is readily found using the shallow water theory. However, near the head of the mine waste dam-break flood wave the laminar resistance resistance effects are appreciable. In this paper the retardation of the wave-front behind the position predicted by the inviscid shallow water theory is found.     

An experimental study of flow, bedload transport and bed topography under conditions of erosion and deposition and comparison with theoretical models

The nature of flow, sediment transport and bed texture and topography was studied in a laboratory flume using a mixed size-density sediment under equilibrium and non-equilibrium (aggradational, degradational) conditions and compared with theoretical models. During each experiment, water depth, bed and water surface elevation, flow velocity, bed shear stress, bedload transport and bed state were continuously monitored. Equilibrium, uniform flow was established with a discharge of about 0.05 m³ s^?1, a flow depth of about 0.01 m, a flow velocity of about 0.81–0.88 m s^?1, a spatially averaged bed shear stress of about 1.7–2.2 Pa and a sediment transport rate of about 0.005–0.013 kg m^?1 s^?1 (i.e. close to the threshold of sediment transport). Such equilibrium flow conditions were established prior to and at the end of each aggradation or degradation experiment. Pebble clusters, bedload sheets and low-lying bars were ubiquitous in the experiments. Heavy minerals were relatively immobile and occurred locally in high concentrations on the bed surface as lag deposits. Aggradation was induced by (1) increasing the downstream flow depth (flume tilting) and (2) sediment overloading. Tilt-induced aggradation resulted in rapid deposition in the downstream half of the flume of a cross-stratified deposit with downstream dipping pebbles (pseudo-imbricated). and caused a slight decrease in the equilibrium mean water surface slope and total bedload transport rate. These differences between pre- and post-aggradation equilibrium flow conditions are due to a decrease in the local grain roughness of the bed. Sediment overloading produced a downstream fining and thinning wedge of sediment with upstream dipping pebbles (imbricated), whereas the equilibrium flow and sediment transport conditions remained relatively unchanged. Degradation was induced by (1) decreasing the downstream flow depth (flume tilting) and (2) cutting off the sediment feed. Tilt-induced degradation produced rapid downstream erosion and upstream deposition due to flow convergence with little change to the equilibrium flow and sediment transport conditions. The cessation of sediment feed produced degradation and armour development, a reduction in the mean water surface slope and flow velocity, an increase in flow depth, and an exponential decrease in bedload transport rate as erosion proceeded. A bedload transport model predicted total and fractional transport rates extremely well when the coarse-grained (or bedform trough) areas of the bed are used to define the sediment available to be transported. A sediment routing model, MIDAS, also reproduced the equilibrium and non-equilibrium flow conditions, total and fractional bedload transport rates and changes in bed topography and texture very well.     

基于结构网格的溃坝水流数值模拟

针对溃坝水流数值模拟面临的复杂地形和不规则边界等问题,基于结构网格建立了适应复杂地形和不规则边界的溃坝水流数值模拟有限体积模型（HydroM2D）。模型基于具有守恒特性的二维浅水方程,利用HLLC格式的近似Riemann解计算网格界面通量,利用MUSCL-Hancock法不断向前积分,使模型在时空上具有二阶精度;对源项进行离散处理确保模型的稳定性;模型引入有效干湿边界和不规则地形边界处理方法,准确模拟了干湿单元的动态交替和复杂边界上的水流特性。最后分别利用水槽试验、物理模型和实际算例对模型进行验证。结果表明,该模型对不同情景下的溃坝洪水模拟结果和实测资料以及现有模型模拟结果具有较高的一致性,模拟精度较高,稳定性较好,具有推广应用价值。     

明渠均匀流Q结构分布及运动特性

为了解明渠均匀流湍流Q结构(Quadrant events)第一象限事件(Q1)、第二象限事件(Q2)、第三象限事件(Q3)和第四象限事件(Q4)的时均与瞬态运动特性,使用高频粒子成像测速系统对3种明渠湍流二维瞬时流速矢量场进行测量,统计分析明渠湍流脉动流速场中Q结构的雷诺应力、发生频率、占据测量区域面积率以及Q结构与涡结构运动关系。结果表明:Q1与Q3结构时均特性相似,具有发生频率低、面积占据量少、负瞬时脉动动量通量小等特点;而Q2与Q4结构时均特性与之相反,且Q4结构发生频率沿水深增加而逐渐增多至0.8H附近达到最大,Q2结构发生频率在壁面区域为最大,沿水深增加而减小。明渠瞬态运动特性主要表现为水面区造成Q4结构向下运动,挤压、推进、爬升床面区的Q2结构伴随其自身形变的瞬态发展过程。     

Flow behaviour of tailings paste for surface disposal

Surface disposition of mine tailings in paste form is a new disposal technique, and to achieve a desired depositional geometry, it is necessary to characterize the paste's flow properties. This study investigates how water content affects the flow behaviour and depositional geometry of tailings and kaolinite pastes, which are shear-thinning, high solids content, mineral pastes. A stress-controlled rheometer and a strain-controlled viscometer with vane fixtures were used to characterize the yield behaviour of the pastes, and three types of yield stress were determined. A flume apparatus was used to simulate paste deposition under laboratory conditions. The depositional angle, determined from the flume tests, and the yield stresses, determined from the rheometers, decreased as the water content increased. For each type of yield stress, a linear relationship was found between the depositional angle and the Sofra–Boger dimensionless group (τ_y′Fr/Re), with the linear coefficient depending on paste type.     

榆木沟尾矿库自流式辐射井排渗系统施工关键技术

浸润线过高导致的渗流状况不良已经成为许多尾矿库的主要安全隐患,是引起尾矿库溃坝事故的最直接原因之一。自流式辐射井排渗系统具有自流排水安全可靠、辐射式集渗管水平伸入坝体范围广、排水量大、水位下降快、成本低、使用寿命周期长、后期管理运行方便等优势。在榆木沟尾矿库2号大口辐射井工程中应用了自流式辐射井排渗系统,取得了良好的技术和经济效果。结合该工程实践,对自流式辐射井排渗系统在尾矿库排渗工程中的施工关键技术问题进行了研究和探讨,为今后类似工程提供了有益的经验。