浅水方程模拟中确定糙率的回归系数法

基于水流模拟本身就是模拟模型回归求真(即实体原型)过程的认识,从模拟模型和实体原型之间的相关关系出发,应用统计学中的线性回归分析方法,通过回归系数建立了模型与原型糙率系数之间的关系式,提出了二维浅水方程模拟中糙率确定的回归系数法。数值试验表明,该方法并不依赖于糙率系数的初值,具有稳定性好、收敛速度快的特点;对磨刀门水道水流的数值模拟应用显示,该方法能有效地改善模拟计算的精度,具有较好的适应能力,简便实用,为二维浅水流动模拟计算时率定糙率系数提供了一种有效途径。     

河道洪水实时概率预报模型与应用

通过数据同化方法合理地将实时水文观测数据融入到洪水预报模型中,可提高洪水预报模型的实时性和精确度。选取沿程断面流量、水位和糙率系数作为代表水流状态的基本粒子,以监测断面实测水位数据作为观测信息,建立了基于粒子滤波数据同化算法的河道洪水实时概率预报模型。模型应用于黄河中下游河道洪水预报计算的结果表明,采用粒子滤波方法同化观测水位后,不仅可以直接校正水位,同时也可以有效地校正流量和糙率,为未来时刻模型预报计算提供更准确的水流初始条件和糙率取值区间,进而有效地提高模型预报结果的精度,给出合理的概率预报区间。不同预报期的预报结果表明,随着预报期的增长,同化效果减弱,模型预报结果的精度会有所降低,水位概率预报结果受粒子间糙率不同的影响不确定性增加,而流量概率预报结果受给定模型边界条件的影响不确定性降低。所提出模型可以有效同化真实水位观测数据,适合应用于实际的洪水预报工作中。     

Numerical modeling of the Sakuma Dam reservoir sedimentation

The present study attempts to predict the reservoir sedimentation in 32 km region of the Tenryu River between the Hiraoka and Sakuma Dams in Japan. For numerical simulations of the reservoir sedimentation, the one-dimensional model of the Hydrologic Engineering Centre-River Analysis System (HEC-RAS) is used together with the inclusion of channel geometry, bed gradation curve, Exner-5 bed sorting mechanisms, fall velocity of the particle, and flow and sediment boundary conditions pertaining to modeling region. The modeling region of the Tenryu River is divided into 48 river stations with 47 reaches in the numerical simulations. The numerical model is calibrated using the available data for 48 years from 1957 to 2004. The formulae of sediment transport function, Manning’s roughness coefficient, computational increment and fall velocity have been identified for getting the best estimation of the Sakuma Dam reservoir sedimentation. Combination of obtained sensitive parameters and erodible limits of 2 m gave the best comparison with the measured bed profile. The computed results follow the trend of measured data with a small underestimation. Although Manning’s roughness coefficient has an effect on the sedimentation, no direct relation is found between the Manning’s roughness coefficient and reservoir sedimentation. It is found that the temperature of water has no effect on the reservoir sedimentation.     

黄河下游动床阻力变化及其计算方法

动床阻力在冲积河道洪水演进与河床冲淤计算中具有十分重要的作用。黄河下游不同水沙条件下的床面形态变化较大,动床阻力变化规律十分复杂,因此需要研究动床阻力的计算方法。利用黄河下游花园口、高村、利津等7个水文站1958—1990年的686组实测数据,确定了影响动床阻力变化的关键水沙因子——水流弗劳德数（Fr）与相对水深（h/D₅₀）,前者表示水流强度,后者表示床面相对粗糙度;建立了基于水流能态分区的动床阻力计算公式,并采用这些实测数据率定了公式中的相关参数;利用黄河下游各水文站1991—2016年的2 288组实测资料,进一步验证了公式的计算精度。计算结果表明:动床阻力的大小随弗劳德数或相对水深的增加而减小;基于水流能态分区动床阻力公式的计算精度明显高于未分区的公式及其他4个动床阻力公式,且决定系数（R2）总体接近0.80,说明水流强度与床面相对粗糙度对动床阻力影响十分显著。     

冰湖溃决洪水的二维水动力学数值模拟

为分析冰湖溃决洪水过程,建立了一个二维非恒定流数值模型,可以定量的开展冰湖溃决洪水预测.数值模型的基本控制方程采用二维浅水水流方程,数值方法采用可以捕捉激波的WAF TVD二阶精度格式,并用非界面追踪的方法处理干湿边界.通过理想情况下一维解析解的验证和二维复杂地形情况的数值试验,证明数值模型具有较好的精度和可靠性,并且可以适应复杂地形.将数值模型应用于一个典型的冰湖溃决事件,定量地预估了该冰湖溃决洪水过程,结果表明:水流流速很高,床面剪切应力很大,可能诱发活跃的泥沙运动,形成泥石流灾害.     

Paleohydraulic interpretation of Quaternary alluvium near Golden,Colorado

Quaternary terraces and pediments along Ralston Creek and Clear Creek, near Golden, Colorado, are associated with Verdos, Slocum, Louviers, and Broadway Alluviums. Terrace deposits can be locally correlated on the basis of elevation and relict paleosols. The terrace sediments probably represent aggradation by braided streams flowing from glaciated drainage basins. Engineering hydraulic calculation procedures suggest that flood flows were 2–3 m deep on steep gradients (0.008–0.01). Discharges were as great as 1400 m³/sec, nearly an order of magnitude greater than modern flood discharges. The most useful paleohydraulic calculation techniques were found to be the dimensionless shear approach applied to stream competence and bedload function theory applied to stream capacity.     

泥石流对桥墩冲击力的试验研究

泥石流冲毁桥墩是桥梁在遭受泥石流冲击时的常见破坏形式。为了研究泥石流对桥墩的冲击力大小,通过调整黏土、沙、石子、水的不同含量,配置不同流变特性、不同密度的泥石流,使用所配置的原料在泥石流槽内对两种形状（圆形、方形）的桥墩缩尺模型进行冲击,综合考察了流变特性、流速、桥墩形状以及冲击力的关系。试验表明：试验配置的泥石流原料流变特性差异明显,且可以用简单的选择流变仪测得,用牛顿流体或宾汉体描述。泥石流的流速可用曼宁公式求得,而公式中的糙率系数与泥石流黏度满足幂函数关系。相同工况下,不同形状桥墩所受的冲击力差异明显,方形桥墩阻力系数普遍大于圆形桥墩。使用非牛顿流体雷诺数（Re）可以综合反映流变特性和流速,因此,圆墩的阻力系数可表达为Re的函数,而方墩则没有明显关系。为方便工程应用,可根据黏性泥石流、稀性泥石流对圆墩的阻力系数分别为2.3、0.9,对方墩分别为2.6、1.9进行选用。     

GIS支持下的小流域设计洪水过程线推求

浙江省衢州庙源溪为典型的山区河流,庙源溪雨量丰富,河流比降大,源短流急,洪水陡涨陡落,汇流速度快,洪峰流量大。基于GIS的不规则三角网和空间分析功能,实现水文特征值等值线的内插。根据GIS获取的水文特征值,计算不同历时平均点雨量,考虑点面折算系数,得到平均面雨量。通过频率分析,获得不同频率下的设计雨量。采用瞬时单位线法,进行汇流计算,推求设计洪水过程线及洪峰流量,为制作山洪风险图,建立小流域防洪避洪保障体系提供依据。     

Reconstruction and post-event analysis of a flash flood in a small ungauged basin: a case study in Slovak territory

Flash floods are one of the major natural hazards occurring in small streams with a negative effect on the country as well as on human lives. Heavy rainfall occurred on July 20, 2014 and July 21, 2014 and caused severe surface water flooding and a flash flood in the Malá Fatra National Park (Slovakia). The most affected was Vrátna Valley with the Varínka stream. This study presents a reconstruction and post-event analysis of a flash flood on small ungauged basin located in this protected area of Slovakia. The reconstruction included hydraulic terrain measurements on estimating the flood’s culmination and documenting the flood’s development. The measurements were taken at three cross sections of the Varínka stream. This paper is focused mainly on post-event analysis of the Varínka stream in two profiles: Strá?a (gauged profile) and Tiesňavy (ungauged cross section). Subsequently, the extremeness of the flash flood was preliminary evaluated. Results of the post-event analysis showed that the July 2014 flood was not the highest flood in this area despite its catastrophic consequences. By studying historical materials, we came to the conclusion that in the past (e.g. in 1848 or 1939) some devastating floods in this area had occurred, which had disastrous consequences for the population. The second part of the study is focused on comparing this flash flood with three major floods which have occurred in Slovak territory since 1998. The first flood occurred on the 20th of July, 1998 on the Malá Svinka stream, and the two others are floods which occurred on the 7th of June, 2011 in the Small Carpathian Mountains: on the Gidra stream in Píla village and on the Parná stream in Horné Ore?any village. Such comparison of flash floods from different geographical regions and different rainfall events can provide comprehensive information about their regimes, threats and disastrous effects.     

北京市朝阳区沙子营垃圾污染物在潜水含水层中迁移规律的现场模拟实验

选择北京平原区水文地质、环境地质等方面都比较典型的沙子营垃圾堆放场,建立了其水文地质模型。在充分收集资料、分析得出部分计算所需参数后,采用现场弥散实验、勘查取样测试等方法,求得了该含水层的弥散系数等参数;用二维非稳定流溶质运移方程对污染物在此含水层中的迁移扩散规律、速度和污染范围等进行了模拟计算;采用现场钻探、取样测试分析等方法,评价了该含水层的实际污染状况。实验模拟计算和现场调查结果表明:污染物在潜水含水层中的运移规律遵循二维非稳定流场中的溶质运移方程,污染物在潜水含水层中的运移速度约为86.25m/a,迁移扩散主要发生在地下水流向上,侧向扩散宽度极小,是地下水流向上的1/17。     

Ice stream sticky spots: A review of their identification and influence beneath contemporary and palaeo-ice streams

Rapidly-flowing ice streams are the arterial drainage routes in continental ice sheets and exert a major influence on ice sheet mass balance. Recent observations have revealed that ice stream flow exhibits considerable variability, with relatively rapid changes taking place in speed and direction. This spatial and temporal variability is intimately linked to the conditions at the base of the ice streams and the distribution of localised patches of basal friction, known as ‘sticky spots’. In this paper, we provide a detailed review of sticky spot observations from both contemporary and palaeo-ice stream beds in order to better understand their nature and influence. Observations and theoretical considerations reveal four primary causes of ‘stickiness’: (i), bedrock bumps; (ii), till-free areas; (iii), areas of ‘strong’ (well drained) till; and (iv), freeze-on of subglacial meltwater. These may act together in one location, or in isolation; and a progressive increase in their distribution could lead to ice stream shut-down. Bedrock bumps are influential under active ice streams, where they provide form drag and can create thinner ice which increases the likelihood of basal freeze-on. Increased bed roughness may prevent the lateral migration of some ice streams but bedrock bumps are unlikely to cause ice stream shut-down because, over long time-scales, ice stream erosion might be expected to reduce their amplitude. The influence of till-free areas beneath an ice stream will depend critically on the amount of water that might be drawn out of the surrounding till to lubricate such areas. They are likely to be most important in ice stream onset zones but their identification has proved difficult beneath active ice streams. If an ice stream operates solely by till deformation, it is conceivable that a progressive increase in the exposure of till-free areas could lead to shut-down through a process of sediment exhaustion. Areas of strong, well drained till have been identified beneath both active and ancient ice streams and are most likely to result from the reorganisation of subglacial meltwater. The collapse of an inefficient ‘cannalised’ system to a more efficient ‘channelised’ system can occur rapidly and this mechanism has been hypothesised as a candidate for ice stream shut-down in both contemporary and palaeo-settings. Basal freeze-on has also been observed and inferred from beneath modern and palaeo-ice streams, and a reduction in basal meltwater supply coupled with ice stream drawdown and the advection of cold ice increases the likelihood of switching off an ice stream. A paucity of data from ice stream sticky spots limits a better understanding of their nature, distribution and evolution beneath ice streams. Future technological advances are likely to improve the resolution of the data collected from the beds of modern ice streams but well-preserved palaeo-ice stream beds also hold potential for investigating their influence on ice stream flow and we present simple landsystems models to aid their identification. Such data will considerably enhance the basal boundary condition in ice stream models which will, ultimately, refine our predictions of the response of contemporary ice sheets to future changes in climate.     

水箱模型在大尺度流域实时洪水预报模型研制中的应用

水箱模型用于实时洪水作业预报的具体成果尚不多，本文研究的目的在于将水箱模型用于大尺度流域实时洪水预报。本文介绍了应用水箱模型建立实时洪水预报模型的方法。按照河段流量传播时间将寸滩以上干、支流划分为若干子河段，各子河段按照计算时段长分成若干单元河段，各单元河段区间降雨径流预报采用水箱模型．河道流量演算采用连续马斯京根法。河系预报模型精度在85％以上，能够满足实时洪水预报的要求。     

黄河下游非恒定输沙数学模型——Ⅰ模型方程与数值方法

构建起具有通用性的黄河下游一维非恒定输沙数学模型.该模型建立了新的泥沙连续性方程与河床变形方程,克服了以往数学模型计算中取饱和恢复系数小于1等缺陷,引入了符合黄河下游河道水沙特点的水流挟沙力和河床糙率计算等公式,给出了悬移质含沙量以及悬移质泥沙平均粒径沿横向分布的计算方法,以及阐明了河槽在冲淤过程中河宽变化规律的模拟技术.运用Preissmann四点差分格式离散水流方程,并与泥沙连续性方程进行非耦合求解.     

Modified governing equation and numerical simulation of seepage flow in a single fracture with three-dimensional roughness

Roughness and tortuosity influence groundwater flow through a fracture. Steady flow through a single fracture can be described primitively by the well-known Cubic Law and Reynolds equation with the assumption that the fracture is made of smooth parallel plates. However, ignoring the roughness and tortuosity of the fracture will lead to inaccurate estimations of the flow rate. To obtain a more accurate flow rate through a rough fracture, this paper has derived a modified governing equation, taking into account the three-dimensional effect of the roughness. The equation modifies the Reynolds equation by adding correction coefficients to the terms of the flow rates, which are relative to the roughness angles in both the longitudinal and transverse directions. Experiments of steady seepage flow through sawtooth fractures were conducted. The accuracy of the modified equation has been verified by comparing the experimental data and the theoretical computational data. Furthermore, three-dimensional numerical models were established to simulate the steady flow in rough fractures with the triangular, sinusoidal surfaces and the typical joint roughness coefficient (JRC) profiles. The simulation results were compared with the calculation results of the modified equation and the current equations. The comparison indicates that the flow rate calculated by the modified equation is the closest to the numerical result.     

山区阶梯式河道水流阻力研究

山区阶梯式河道是自然界来水来沙和山区大比降河床相互作用所形成的,其水力特性、河道阻力较为复杂。为了深入研究其阻力机理,从河床结构阻水的角度提出阶梯式结构糙率尺度的概念,得到了适用于阶梯式河道跌落水流的形体阻力系数表达式,并由不同来源的阻力分量构成了新的山区阶梯式河道水流总阻力系数公式。新公式结果表明,通过不同类型河段的野外试验数据,表面泥沙、阶梯形状以及漂石产生的阻力系数分量之和与试验测量总阻力系数数据的误差较小。随着雷诺数和阶梯单元尺寸的增加,河床表面的泥沙肤面阻力系数基本不变,阶梯式结构形成的跌落水流过渡到滑行水流,阶梯形体阻力系数逐渐变小,而较大粒径松散漂石产生的水流阻力系数只随着雷诺数的增大而减小。     

长江流域“2017·07”暴雨洪水分析

2017年6月下旬至7月初,受持续强降雨影响,长江发生中游区域性大洪水。以实时报汛数据为基础,分析长江"2017·07"暴雨洪水特性,依据洪峰水位判断,强降雨导致洞庭湖水系湘江发生超历史最高水位特大洪水,资水、沅江发生超保证水位大洪水,洞庭湖超过保证水位;鄱阳湖水系乐安河上游发生超历史最高水位特大洪水,昌江、乐安河中下游、修水发生10a一遇较大洪水,鄱阳湖超过警戒水位;长江干流莲花塘以下江段全线超过警戒水位。在应对此次洪水过程中,长江上中游重点水库防洪效益十分明显,有效避免中游干流莲花塘至螺山江段超保,缩短洞庭湖城陵矶站超保时间6d左右。     

冰封河道综合糙率计算方法比较与分析

综合糙率是冰封河道水力计算的重要输入参数。根据水流连续性和力的平衡原理,分别推导得到冰封河道综合糙率的通用公式,并选取抛物线形、矩形和梯形作为典型断面,系统地总结Pavlovskiy、Einstein、Lotter、Sabaneev、Larsen、改进的Larsen等常见公式由通用公式基于何种假定（分区的湿周、流速或水力半径相等）简化得到。采用水槽实测数据和野外观测资料,对比上述公式的计算精度,结果表明改进的Larsen公式的计算结果与实测值吻合较好,推荐在工程设计中使用。前5种公式误差较大的原因在于分区的流速、水力半径相等对于冰封河道不具有普适性,尤其是Lotter公式和Pavlovskiy公式在冰盖和河床糙率差异较大时会明显低估冰封河道的综合糙率,在工程设计中应慎重采用。     

GIS applications in morphometric analysis of Koshalya-Jhajhara watershed in northwestern India

The present paper deals with morphometric analysis of Koshalya-Jhajhara (K-J) watershed, tributaries of Ghaggar in northwestern India. The area forms a rugged topography having elevation range from 399 m to 1810 m MSL. Aster DEM and SOI toposheets have been used to calculate various parameters using Geographic Information System (GIS). Linear, areal and relief aspects for morphometric analysis were calculated using GIS. Total number of streams in the watershed are 991 out of which 543 are of 1^st order, 259 are of 2^nd order, 124 are of 3^rd order, 58 are of 4^th order and only 7 are of 5^th order. The order of stream determines the chances of flood in the stream. Higher the stream order more the probabilities of flood. The present paper indicates that K-J watershed is of 5^th order and less elongated in shape having lower peak flows of longer duration with dendritic pattern and having fine drainage texture. Bifurcation ratio varies from 3 to 5 which indicate that geological structures don’t have dominant influence on drainage pattern. High slope is witnessed in NE part and low slope in SW part with very low gradient ratio. Variables like stream frequency and drainage density determines the volume of the water discharge and its speed of flow in the river channels. Higher the stream frequency and drainage density more will be the probability of floods. These variables also effects temporal variations in the speed at which the water flows in the stream when flood reaches its peak.     

山区卵砾石河道水流阻力计算

卵砾石河道广泛存在于山区河流中,在卵砾石河床近底层的水流流速低于上方流速,速度剖面出现拐点,类似于自由剪切流动,传统的指数型和对数型水流阻力公式计算误差偏大。将自由剪切流理论类比到山区卵砾石河道,并考虑山区河流特殊流态、流场和剪切力的影响,引入雷诺数(Re)、弗劳德数(Fr)和摩阻流速利用水槽数据进一步修正了该类水流阻力公式。筛选位于意大利南部的亚平宁山脉140条山区河流野外数据对新公式进行验证,并检验对比了已发表的多个山区河流阻力公式。结果表明:修正后的水流阻力公式Nash-Sutcliffe效率指数最接近1、均方根误差最小、相对误差最小。自由剪切流动的类比是从湍流结构角度推导河道水流阻力,能更好地阐释特殊水流结构,进而提高计算精度。     

入湖河道生态修复工程对水流流速的影响研究

在河流中开辟生化反应区、直接布设仿生填料净化河流微污染水体可以达到改善入湖水质、减少投资、降低能源消耗的目的。但在河道内直接布设填料势必会对河流的水动力条件及其河流原有的使用功能造成一定的影响。从填料给水流带来的阻力、流速及糙率变化等方面进行河道内直接布设生物填料对河流水动力条件的影响试验分析。结果表明,由于河道流速不大(属于缓流水体),9根/m2的填料布设密度对水流的扰动作用程度不高,填料段的糙率仅为无填料段糙率的1.08倍,不会对河道的泄水能力造成影响。