长江中下游河湖洪水演进的数值模拟

以长江中下游(宜昌-大通)防洪系统为对象,在水文、河道及湖盆地形现状条件下,建立了一个能适应各种复杂条件的一二维非恒定流模型,来进行长江干流、河网、湖泊、分蓄洪区垸及水库的洪水演进和调度仿真.所建模型的洞庭湖部分采用无结构网格二维非恒定有限体积格式,以适应湖区复杂的边界形状和保持水量平衡.河网区部分采用一维非恒定流显隐结合的分块三级算法,以准确实现河网汊点流量的自动分配和往复流动.为了提高模拟精度和扩展模拟功能,在水流数值模拟的范围内侧重讨论了内外动边界处理、分蓄洪运用及阻力项计算等环节,提出了合理可行的数值处理方法.采用20世纪80年代至90年代共6年汛期洪水资料对所建模型进行了严格的率定和检验,高精度的模拟结果证实了模型的合理性和有效性.     

不稳定河床非恒定流的水位流量关系校正计算

本文根据非恒定流运动的基本原理,并进行一系列简化假设,推导出不稳定河床并受多种水力因素影响下的水位～流量关系校正计算公式,用以改进水位～流量关系曲线,提高资料整编成果精度,效果较好,现整理成文,供参考。一、基本原理与公式推导明渠水流运动可用圣维南方程组描述,即     

漫滩水流摩阻因子与涡粘性系数的研究

对复式河道漫滩水流进行了研究,建立了摩阻因子与一维明渠流摩阻因子的关系,揭示出滩地的量纲为一涡粘性系数随着摩阻因子的增加而增加,并将二次流的影响表达成与雷诺切应力相同的形式,从而将其归入表观切应力中,由量纲为一表观涡粘性系数确定,后者可通过量纲为一表观涡粘性系数与摩阻因子之间的一个关系加以计算,该关系类似于量纲为一涡粘性系数与摩阻因子之间的关系。对光滑边界上对称的漫滩水流进行了计算,计算结果与实验资料吻合良好。     

明渠非恒定流流速分布及推移质运动研究进展

明渠非恒定流流速分布及非恒定流作用下推移质运动分别是水力学和泥沙运动力学研究中的重要问题之一。近几十年来,国内外许多学者在这些领域做过很多探索,取得了一定进展,也存在较多问题。在总结相关研究成果的基础上,简述了明渠非恒定流试验方法、流速分布、床面切应力、非恒定流输沙与水流运动不同步现象、产生原因、以及非恒定流推移质输沙公式等方面的成果,指出了现有研究中存在的不足,提出了需要进一步研究的关键问题:应用先进的量测仪器,研究非恒定流作用下非均匀沙运动输移规律。     

长江中游洞庭湖防洪系统水流模拟——Ⅰ．建模思路和基本算法

针对长江中游洞庭湖防洪系统规模庞大、水流复杂、资料短缺和预报时限紧迫的实际条件,提出了具有层次和模块结构特点、一维与二维水流模拟、水力学与水文学方法、理论模型与补充信息相结合的建模途径.所建模型的湖泊部分采用无结构网格二维非恒定流高性能有限体积格式,以适应湖区复杂的边界形状和保持水量平衡;河网部分采用一维非恒定流守恒型显格式,避免隐格式矩阵算法的复杂性,同时有利于与二维模型的耦合及与各种复杂连通关系的显式连接.这种一、二维混合非恒定流模型可用于长江干流、洞庭湖河网及湖泊、堤垸区的洪水演进和防洪调度的水流仿真.     

计算明渠非恒定流的一个Lattice BGK模型

分析了Lattice Boltzmann(LB)方法的建模原理,通过比较LB的多尺度方程和明渠St.Venant方程,建立了计算一维明渠非恒定流的一个Lattice BGK模型。对几个典型应用实例进行了模拟分析,与Preissmann隐格式及特征线法作了比较,证明将该方法引入实际明渠非恒定流计算是可行的。     

非恒定流动能与加速比降对流量计算的影响分析

非恒定流动能与加速比降对流量计算的影响分析罗光兰，章应耕（湖南省水文总站）一、非恒定流基本方程分析天然河道非恒定渐变流基本微分方程可表示为：且假定非恒定流中沿程水头损失与恒定流中沿程水头损失具有相同的规律，将（１）式对流程Ｌ的偏导数合并在一在同一时刻...     

城市暴雨内涝数学模型的研究与应用

以城市地表与明渠、河道水流运动为主要模拟对象,研制了模拟城市暴雨内涝积水的数学模型。模型以平面二维非恒定流的基本方程和无结构不规则网格划分技术为骨架,同时,针对小于离散网格尺度的河道或明渠,应用了一维非恒定流方程的算法。采用分类简化处理的方法,将通道分为河道型、路面型、特殊通道型(城市内的二级河道),根据不同类型简化动量方程,求任意网格各个通道上的单宽流量。采用一维非恒定流方程模拟地下排水管网内的水流,并给出泵站、闸门、淹没出流管道等排水系统的处理方法。根据无结构不规则网格的设计思路,按照天津、南京、南昌三市的地形地貌特征分别设计多边形的计算网格。介绍了城市面雨量的计算方法以及数学模型在天津市、南京市、南昌市的应用情况和误差分析。     

黄河下游高含沙洪水过程一维水沙耦合数学模型

采用浑水控制方程,建立了基于耦合解法的一维非恒定非均匀沙数学模型,用于模拟高含沙洪水演进时的河床冲淤过程.然后采用黄河下游游荡段1977年7—8月实测高含沙洪水资料对该模型进行率定,基于水沙耦合解法的各水文断面流量、总含沙量及分组含沙量的计算过程与实测过程符合更好,计算的沿程最高水位及累计河段冲淤量与实测值也较为符合.最后还采用2004年8月高含沙洪水资料对该模型进行了验证.模型率定及验证计算结果表明,采用一维水沙耦合模型计算高含沙洪水过程,能取得较高的精度.     

基于优化DTW算法的水文要素时间序列数据相似性分析

本文在分析基于日平均水文要素时间序列数据相似性分析的基础上,提出基于水文要素摘录数据的相似性分析优点和存在问题,引入动态时间弯曲距离算法(DTW)解决水文要素摘录数据时序间距不一致,时常不一致的问题,并根据水文时间序列数据周期性起伏变化频繁的特征,基于最长公共子串,对DTW算法进行优化改进,从而提高相似性匹配精度.以汉...     

Hydraulic conductivity prediction based on grain-size distribution using M5 model tree

The hydraulic conductivity, K_s, is one of the most important hydraulic properties which controls the water and solute movement into the soil. It is measured on soil specimens in the laboratory. On the other hand, sometimes it is obtained by tests carried out in the field by a number of researchers. Therefore, several experimental formulas have developed to predict it. Recently, soft computing tools have been used to evaluate the hydraulic conductivity. However, these tools are not as transparent as empirical formulas. In this study, another soft computing approach, i.e. model trees, have been used for predicting the hydraulic conductivity. The main advantage of model trees is that, unlike the other data learning tools, they are easier to use and represent understandable mathematical rules more clearly. In this paper, a new formula that includes some parameters is derived to estimate the hydraulic conductivity. To develop the new formulas, experimental data sets of hydraulic conductivity were used. A comparison is made between the estimated hydraulic conductivity by this new formula and formulas given by other’s researches.     

改进的无作物潜水蒸发经验公式

根据安徽省五道沟水文水资源实验站砂姜黑土和砂壤土2种原状土样1992-1994年的潜水蒸发实验观测数据,在分析无作物条件下潜水蒸发与各种影响因素关系、潜水蒸发系数与水面蒸发关系的基础上,采用非线性最小二乘法,以决定性系数、残差均值和方差为性能指标,对已有的潜水蒸发经验公式-阿维里扬诺夫公式和叶水庭公式进行了非线性回归拟合优度比较,并提出了2个公式的改进形式。结果表明,采用改进后的经验公式进行潜水蒸发的拟合和预测,较原有公式精度有明显提高。     

What can be learned from sequential multi-well pumping tests in fracture-karst media? A case study in Zhangji, China

A fracture-karst aquifer is a karst aquifer with a fractured rock matrix, and its parameters are difficult to determine. Two sequential pumping tests in a fracture-karst aquifer system at the Zhangji well field of China are considered, one carried out before (in 2000) and one after (in 2005) the operation of a pumping station in the well field (2003–2005). The sequential tests serve multiple purposes. First, they provide a cross check of the parameters obtained. Second, they can be used to assess the effect of long-term groundwater exploitation of the aquifer. A three-dimensional finite-element transient flow model has been developed to simulate groundwater flow at the site. Generally good agreement has been found between the simulated and observed hydraulic heads for both tests. The hydraulic parameters obtained from the 2005 test are generally consistent with their counterparts from the 2000 test. However, a small but steady increase of hydraulic conductivities from 2000 to 2005 at the site has been observed. A 10-year prediction of groundwater resources has been made and indicates that the well field can accommodate the proposed 8.0?×?10⁴ m³/day exploitation rate under relative drought conditions without causing a steady decline of groundwater levels.     

Estimating hydraulic conductivity from specific capacity for Hawaii aquifers,USA

Site-specific relationships between specific capacity and hydraulic parameters (transmissivity and hydraulic conductivity) were investigated for volcanic rocks in Maui, Hawaii, USA. Details about well construction were commonly ignored in previous studies. To improve on such efforts, specific-capacity values were normalized by the open interval of the well. Correcting specific capacity for turbulent head losses using step-drawdown tests and including aquifer penetration length improved the correlation between specific capacity and hydraulic conductivity and reduced uncertainty in the prediction of hydraulic parameters. The relationships provide estimates of aquifer parameters with correlation coefficients between 0.81 and 0.99. The relationships for Maui can probably be extended to other Hawaii islands, given the similarity of aquifer formations and a reasonable fit to step-drawdown data from Oahu. Hydraulic conductivity was estimated from 1,257 specific-capacity values in the Hawaii’s well database. Hydraulic-conductivity estimates for dike-free volcanic rocks are consistent on different islands. For all islands, the estimates range from 3 to 8,200 m/d, with a geometric-mean and median value of 272 and 291 m/d, respectively. A geostatistical approach was applied to Maui and Oahu to generate island-wide hydraulic-conductivity maps to facilitate groundwater management efforts.     

关于用稳定流抽水降落曲线q—s关系建立求解水文地质参数(Rs、K、R)新模型的探讨

本文以q/s曲线长度A和单位涌水量q、水位降深s之间关系,用达西渗透公式诱导和假设建立了求地下水向井q/s曲线及q/s曲线段间渗流速度公式,并用q/s曲线段间渗流速度Rs求得含水层渗透系数K,又用渗透系数K建立了反求影响半径R的假设公式。所建立的公式只适用于稳定流抽水试验中计算含水层的K和R的水文地质参数。     

饱和砂层驱水渗透注浆分析

运用渗流力学理论,对饱和砂层在毛细力作用下的平面径向流、球形向心流两种渗透性注浆模型进行了推导,得到了球形和柱形渗透性注浆扩散公式。以球形扩散公式为例,结合工程实例并与常用的Maag球形扩散公式进行比较,结果表明,由Maag公式计算注浆所需时间偏小,文中公式更接近工程实际。最后探讨了注浆参数对扩散半径的影响,表明扩散半径对注浆参数的变化均较敏感,文中公式与Maag公式所得到的变化趋势基本一致。     

均质非饱和土一维稳态流水头垂直分布的解析解

渗流场水头分布计算是进行渗流量和渗流水力坡降计算的基础,准确、有效地求取渗流场水头分布是渗流计算的关键环节。对均质非饱和土体一维稳态流的流动方程进行分析,考虑到渗透系数是与基质吸力相关的函数,通过数学变换,给出了稳定渗流场的解析通式,并基于渗透性函数中的Gardner模型,给出了非饱和土一维稳态流水头垂直分布的解析解。该解析通式表明,均质非饱和土一维稳态流水头垂直分布主要受地表水头、深度和流动率3个因素控制。分别计算了一维稳态蒸发条件下粉土和黏土两种典型土类水头沿垂直方向的分布。计算结果表明：稳态蒸发条件下粉土层和黏土层内的水头分布表现出相似的变化规律,即自地表至地下水位处随着土层深度的增加,水头分布呈现出加速递减的趋势;在相同的蒸发条件下,对于相同深度处的黏土和粉土而言,黏土层内水头更高些;对同一种土类而言,在较大的蒸发状态下同一深度处土层内水头更高。反之,则较低。     

Structural and statistical characterization of joints and multi-scale faults in an alternating sandstone and shale turbidite sequence at the Santa Susana Field Laboratory: Implications for their effects on groundwater flow and contaminant transport

This paper describes the properties of faults and fractures in the Upper Cretaceous Chatsworth Formation exposed at Santa Susana Field Laboratory and its surroundings (Simi Hills, California), where groundwater flow and contamination have been studied for over three decades.The complex depositional architecture of this turbidite consisting of alternating sandstones and shales, interacting with formative stress conditions are responsible for multi-scale fault hierarchies and permeable fractures in which nearly all groundwater flow occurs.Intensity and distribution of background fractures and their relation to bedding thickness are established for sandstones, the dominant lithology. The architecture of faults with increasing displacement is described, and relationships among fault dimensional parameters captured.Data from ∼400 boreholes and piezometers reveal the effect of faults and fractures on groundwater flow. Large hydraulic head differences, observed across fault zones with shale-rich cores, indicate these structures as cross-flow barriers. Moreover, hydraulic head profiles under ambient conditions, and pumping tests suggest strong hydraulic connectivity in all directions to depth of hundreds of meters.This outcrop-based structural characterization relates the horizontal hydraulic conductivity to the observed well-connected fracture network, and explains the strong vertical connectivity across low-hydraulic conductivity shales as faults and sheared fractures provide flow pathways.     

应用地下水匀坡与加坡原理计算渗透系数

应用匀坡及加坡原理,推导出了单井（孔）稳定流承压完整井渗透系数的计算公式,该公式已考虑到无效降深对渗透系数计算的影响,同时在计算过程中充分利用了水位恢复资料,而水位恢复的过程也排除了一些干扰因素的影响,因此,参数计算结果应更接近实际。     

一维稳态渗流时非饱和土中的基质吸力分布

首先对非饱和土中的土水势和干土吸水过程中土水势的变化进行了定性分析。然后利用水相连续性条件、广义Darcy定律及Gardner关于非饱和土导水系数和基质吸力的经验公式,建立了非饱和土一维稳态渗流情况下基质吸力沿深度变化的关系表达式,并利用差分方法对稳态入渗和稳态蒸发两种情况进行了求解。根据提出的方法,不仅可以方便地得到基质吸力沿深度的分布,而且可以求出水头和导水系数沿深度的变化。