分汊河流江心洲洲头冲淤概化模型

李志威^1,2，李艳富³，王兆印²，胡世雄¹

（1.长沙理工大学水利工程学院，湖南 长沙 410004；2.清华大学水沙科学与水利水电工程国家重点实验室，北京 100084；3.南京水利科学研究院，江苏 南京 210029）

摘要：分汊河流的江心洲洲头是分汊河段分水分沙和泥沙冲淤的动态区域，其冲淤过程受来水来沙和地形条件变化呈现冲淤交替。通过几何概化洲头的浅滩逆坡和侧向顺坡，考虑逆坡促淤和顺坡分沙作用，建立洲头浅滩的推移质冲淤的概化模型，推导洲头由淤积转向冲刷的临界流量表达式。上游来流量增加，洲头淤积速率先增加到最大值后逐渐减少，直至达到临界流量，再转向冲刷。洲头淤积速率与上游水面纵比降和床沙粒径成正比，与浅滩逆坡、侧向水面比降与顺坡成反比。若其他条件相同，洲头浅滩形状越不对称和上游来沙量越小，越有利于冲刷。模型预测武汉河段天兴洲的临界流量范围与水文观测值较为接近，对其出现超过某个大临界流量仍出现淤积的特殊性，给出的新解释是天兴洲洲头的整治工程使得在大流量条件下洲头抵抗冲刷和阻挡作用促进泥沙淤积。

关键词：分汊河流；江心洲；淤积；冲刷；概化模型；临界流量
小尺度径流空间差异性及成因探索

——以青海省黄南州为例

段水强^1,2

（1.青海省水文水资源勘测局，青海 西宁 810001；2.青海省自然地理与环境过程重点实验室，青海 西宁 810008）

摘要：为探索资料稀疏地区小尺度径流的空间差异，2013年8月14—23日，对青海省黄南州50多条10~100 km²级别的河流进行了一次流量监测，并进行了时间升尺度和成因的探索。结果发现：①小流域日流量模数差异最大可达20多倍，降水条件相似的邻近河流，最大差异可达4倍以上；②研究区3处中等尺度水文站平水期存在良好的径流同步性，借助这一特性，可以将小流域的日尺度径流特征推绎到年尺度；③小尺度上，降水、坡度是径流空间差异的主导因素，构建的降水-降水损失-坡度径流模型，为实现当地的径流空间尺度转换提供了一座 “桥梁”。

关键词：径流；同步性；尺度；空间差异；青海黄南
坡面不同侵蚀沟断面特征及水力几何形态

赵春红^1,2，高建恩^1,3

（1.西北农林科技大学，陕西 杨凌 712100；2.清华大学水沙科学与水利水电工程国家重点实验室，北京 100084；3.中国科学院水利部水土保持研究所，陕西 杨凌 712100）

摘要：通过整理和分析国内外大量野外调查数据，系统研究了裸露坡耕地壤土土质不同侵蚀沟的断面特征和水力几何形态，并给出了相应的判别标准。结果表明：随坡面侵蚀沟逐步发展，宽深比参数ζ逐渐减小，细沟、浅沟、切沟和冲沟适宜的ζ范围分别为5~17、2~5、0.4~1.8和0.1~0.6；同一侵蚀沟条件下，壤土类型对ζ值也有影响。土壤质地越细，ζ值越小。不同侵蚀沟水力几何形态指数也不相同，细沟的形态指数β₁、β₂和β₃分别为0.3、0.4和0.3，浅沟和切沟为0.4、0.4和0.2，冲沟为0.5、0.4和0.1。验证结果表明给出的判别标准可以用来界定不同侵蚀沟类型。

关键词：坡耕地；侵蚀沟；宽深比；水力几何形态；判别标准
采用间歇灌溉进行土壤盐分淋洗的适用性

彭振阳，伍靖伟，黄介生

（武汉大学水资源与水电工程科学国家重点实验室，湖北 武汉 430072）

摘要：为比较一次性灌溉和间歇灌溉在盐分淋洗效率方面的优劣，采用HYDRUS-1D两区模型模拟分析了不同频次灌溉条件下的土壤水盐运移过程。结果表明存在一个间歇灌溉具有更高淋洗效率的最大深度，并将其称为临界深度。对砂土的模拟显示在灌水总量为500px时，随着潜在蒸发强度从0mm/d升高至6mm/d，临界深度从1500px以上降为0；而在潜在蒸发强度为2mm/d时，随着灌水总量从250px增加至1000px，临界深度相应从500px增加至接近2000px。比较而言，透水性好、不动水体占比高、两部分孔隙之间水分和溶质交换能力差的土壤，临界深度更大；对黏土、黏壤土、壤土、砂壤土和砂土分别进行模拟显示偏砂性土壤具有更大的临界深度，在灌水总量和潜在蒸发强度分别为750px和2mm/d时，黏土和砂土的临界深度分别为1400px和大于2500px。总的来看，对于间歇灌溉是否能提高盐分淋洗效率不能一概而论，引入临界深度可以在一定程度上解释研究者们对于间歇灌溉是否能提高盐分淋洗效率的不同认识。

关键词：间歇灌溉；盐分；HYDRUS-1D；非均匀流；土壤质地
土壤添加物对土-水曲线和土体收缩的影响

邢旭光^1,2，柳 烨^1,2，马孝义¹

（1. 西北农林科技大学水利与建筑工程学院，陕西 杨凌 712100；2. 西北农林科技大学旱区农业水土工程教育部重点实验室，陕西 杨凌 712100）

摘要：为探索不同添加物（小麦秸秆、麦壳和植物油渣）对土壤持水性能和土壤收缩特性的影响，定量掺混3种粉末，采用离心机法测定土-水曲线，并对土壤持水能力和土壤收缩特征进行评价。结果发现：各处理土-水曲线形态无差异，van Genuchten模型可用于此条件土-水曲线拟合；3种添加物均可提高土壤容重和土壤含水率并利于保水，且土壤持水能力表现为油渣>麦秆>麦壳>纯土；随着土壤脱水，土体以轴向收缩为主，三直线模型可用于土壤收缩曲线拟合，土体线缩率与吸力和含水率分别满足对数和Logistic函数关系；土体几何收缩量表现为麦壳<油渣<麦秆<纯土，土壤开裂程度表现为油渣<麦秆<麦壳<纯土。

关键词：水分特征曲线；土壤收缩；土壤裂缝；土壤添加物
横向振动格栅紊流紊动特性试验

李文杰¹，刘 华²，杨胜发¹，王 涛³，张 鹏²

（1．重庆交通大学国家内河航道整治工程技术研究中心，重庆 400074；2.重庆交通大学水利水运工程教育部重点实验室，重庆 400074；3.长江重庆航运工程勘察设计院，重庆 401147）

摘要：建立了一套大尺度格栅紊流试验系统，格栅进行有别于传统垂向模式的横向振动，采用粒子图像测速技术（Particle Image Velocimetry，PIV）测量流速。对瞬时流速的检验表明，该系统产生的紊流场具有较强的随机性和统计规律性。均方根流速在格栅片附近变化较大，在两片格栅中间处趋于稳定，纵向均方根流速明显大于垂向均方根流速，二者比值在1.5~2.0之间，接近天然明渠紊流。雷诺应力在距格栅越近处波动越大，随着距格栅距离的增加而减小，至两片格栅中间处雷诺应力基本为0。时间和长度积分尺度在格栅片处最小，随着距格栅距离的增加而线性增加，至两片格栅中间处达到最大值。流速能谱呈现Kolmogorov理论的-5/3次方规律。本系统生成的紊流场的统计规律与传统的垂向振动格栅紊流较为一致，但纵向和垂向的紊动强度更接近实际，为后续紊流中泥沙和污染物等运动机理的研究奠定了基础。

关键词：格栅紊流；横向振动；紊动特性；粒子图像测速
长江口枯季北港淡水向北支扩展的动力机制

李 路^1,2，朱建荣³

（1. 上海市水利工程设计研究院有限公司，上海 200061；2. 上海围海工程技术研究中心，上海 200061；3. 华东师范大学河口海岸学国家重点实验室，上海 200062）

摘要：2011年12月至2012年1月长江口的现场观测资料表明，大潮和大潮后中潮期间，北港的淡水向北支口门和下段扩展，减小了北支下段的盐度，而小潮及其后的中潮期间无此现象。为了确证淡水的来源，采用验证良好的三维数值模式再现了上述观测现象，并通过数值试验和淡水通量机制分解方法分析其动力成因。数值试验结果表明，大潮及其后的中潮期间，北港主槽流经北港北汊及其周围浅滩的水体是北支淡水的主要来源，拉格朗日余流输送是其向北支扩展的主要动力机制。北港北汊的增深增加了进入北支的淡水通量，减弱了北支的盐水入侵。偏北风风速的增加抑制了北港淡水向北支的扩展，增加了北支的盐水入侵。

关键词：盐水入侵；淡水通量；数值模式；长江口
非饱和带有限分析数值模拟的误差分析

张在勇，王文科，陈 立，王周锋，段 磊，安可栋

（长安大学旱区地下水文与生态效应教育部重点实验室，陕西 西安 710054）

摘要：由于描述包气带水分运移的Richard方程具有高度非线性，常用数值方法求解。有限分析法可以较好地保持原有问题的物理特性，近年来被广泛应用于环境及农业工程领域。基于不同的原理，目前有两种有限分析法模拟非饱和带水分运移问题，即基于局部线性化有限分析法和基于Kirchhoff变换有限分析法。为探讨两者在求解非饱和带水分运移问题时数值表现的差异性，采用已有解析解对两种方法进行初步验证；使用两个数值实验比较两种方法的计算精度；应用实测数据验证Kirchhoff变换有限分析法。相比于局部线性化有限分析法，Kirchhoff变换有限分析法能够更好地控制误差及获得更高精度的数值解。该研究对于完善非饱和带水分运移有限分析数值模拟的理论与方法有重要意义。

关键词：非饱和流；Richard方程；有限分析法；Kirchhoff变换
海湾水库蓄水初期的单次往复密度流

高增文，薛青美，李宇浩，隋玉柱

（青岛大学环境科学与工程学院，山东 青岛 266071）

摘要：为了探讨海湾水库蓄水初期单次往复水平密度流的产生与效应，进行砂质斜坡沉积物与水体之间盐分交换的水槽实验。通过沉积物孔隙水与其上覆水的多剖面电导率动态监测，分析盐分的时空分布特征，并计算深水区、浅水区和中心区的边界层单位面积含盐总量。实验结果显示，进水在浅水区沉积物表面产生明显的渗入-溢出现象，浅水区沉积物孔隙水盐分浓度显著低于中心区与深水区同一高度的盐分浓度，深水区初始底边界层含盐量与初始边界层下边界盐分浓度均为最高；这表明进水过程在倾斜沉积物表面产生了前进密度流，水流携带的盐分在坡底累积。水槽进满水后深水区边界层含盐总量与边界层下边界盐分浓度快速降低，并且浅水区表层沉积物孔隙水与深水区同一深度水体之间的盐分快速达到平衡；这表明蓄水初期在两者之间形成了返回密度流，从浅水区表层沉积物冲刷出的盐分在密度流作用下再次进入沉积物。为避免再次进入沉积物的盐分在后期继续影响水库泛咸，建议在蓄水结束后尽快实施坡底咸水排除方案。

关键词：海湾水库；密度流；往复流；斜坡沉积物；咸化
基于势流理论和OpenFOAM的耦合模型对多向不规则波浪的模拟

季新然^1,2，柳淑学¹， Bo T. PAULSEN^2,3，Harry B. BINGHAM²

（1.大连理工大学海岸和近海工程国家重点实验室，辽宁 大连 116024；2. 丹麦技术大学机械学院，灵比 2800；3.代尔夫特三角洲研究中心，代尔夫特 2629HD）

摘要：基于势流理论的数值水池可以快速计算波浪的传播及其对建筑物的作用，但是势流理论是基于波浪的无黏性假设的，而在工程中，通常需要在固体边界及波浪破碎的区域考虑黏性效应。针对基于求解Navier-Stokes方程的黏性水池计算量较大、速度较慢的缺点，采用耦合的方法模拟多向不规则波浪的传播，即在外域通过基于势流理论的数值水池产生多向不规则波浪，内域采用求解Navier-Stokes方程和流体体积法（VOF方法）对自由表面进行追踪，通过外域所提供的边界波浪，内域计算可以在较小区域进行计算，从而达到减少计算量、提高计算效率的目的。

关键词：势流理论；OpenFOAM；耦合；多向不规则波浪
基于随机理论的地下水环境风险评价

张 博¹，李国秀¹，程 品¹，孙法圣²，王洪义^1,3

（1. 吉林大学地下水资源与环境教育部重点实验室，吉林 长春 130021； 2. 山东省水利勘测设计院，山东 济南 250013；3.吉林省墒情监测中心，吉林 长春 130033）

摘要：针对确定性模型难以描述含水层非均质空间分布的问题，提出基于随机理论的地下水环境风险评价方法。以矩形场地地下水污染风险评价为例，采用蒙特卡罗法生成大量渗透系数随机场，模拟含水层参数各种可能的非均质空间分布，在此基础上建立场地地下水流模型与溶质运移模型，分别计算污染物在地下水中的迁移转化情况。统计大量随机模拟中污染事故发生的频率，当模拟次数足够多时，污染频率收敛于污染概率，污染风险即通过污染概率体现出来。该方法将模型参数设为满足一定分布特征的随机变量，避免了确定性方法得出的武断的评价结果，可为工厂的选址、水源地的选址等工作提供科学指导。

关键词：非确定性模型；环境风险评价；数值模拟；渗透系数随机场
论水库工程之管运洪水

丁 晶¹，何清燕¹，覃光华^1,2，李深奇¹，张泽慧¹

（1.四川大学水利水电学院，四川 成都 610065；2.四川大学水力学与山区河流开发保护国家重点实验室， 四川 成都 610065）

摘要：鉴于水库工程管理运用的重要性，提出管运洪水的新概念，并和设计洪水做了全面对比，指出二者异同。在此基础上提出了评判管运洪水是否合理的两个通用准则：条件合乎要求和防洪安全达标。将推求管运洪水的方法分为两大类：一类以分期内洪水最大值为依据，另一类以年最大值为依据。这两类方法本质不同，使用时各有优缺点。第一类方法的结果不完全符合现行标准的要求，而第二类方法的结果符合要求，但对资料条件要求较高。

关键词：管运洪水；设计洪水；水库防洪调度
水库群水沙调控的单-多目标调度模型及其应用

白 涛，阚艳彬，畅建霞，袁 梦

（西安理工大学西北旱区生态水利工程国家重点实验室培育基地，陕西 西安 710048）

摘要：为缓解内蒙古河段“二级悬河”形势，以黄河上游沙漠宽谷河段为研究对象，以龙羊峡水库、刘家峡水库为调控主体，开展黄河上游水沙调控研究。建立了输沙量、发电量最大的单目标模型以及多目标模型；分别采用自迭代模拟算法、逐次逼近动态规划算法（DPSA）和改进的非支配排序遗传优化算法（NSGA-II）求解模型；设置了初始、常规、优化和联合优化4种方案。通过实例计算，联合优化调度方案的区间总冲刷量达到了0.38亿t，梯级发电量148.22亿kW·h。 该方案以较小的电量损失换来了输沙量的大幅度增加，水沙调控效果显著，推荐为最优方案。研究成果量化了水沙调控效果和各目标间的转化规律，为开展黄河上游水沙调控提供了决策依据，具有重要的应用价值和实际指导意义。

关键词：水沙调控；多目标模型；改进的NSGA-II；联合优化运行
水库群供调水系统实时调度研究

孙万光¹，李成振¹，姜 彪¹，刘盈斐²

（1. 中水东北勘测设计研究有限责任公司，吉林 长春 130061；2. 中国水利水电科学研究院，北京 100048）

摘要：以大连市复杂水库群供、调水系统为背景，在对系统特性分析的基础上，提出了系统实时调度框架及二层耦合结构模式，设计了以“总量控制、耦合嵌套、多维决策、滚动修正”为核心的实时调度流程。在径流预报方面，提出了基于超越概率的水库群供、调水系统长期入库径流预报方法，结合中期GFS（Global Forecasting System）数值预报技术对不同时段入库径流进行滚动预报。在调度模型方面，采用基于动态规划的建模求解新方法和常规调度方法，为调度决策提供参考。实例研究表明，建立的实时调度系统实现了滚动预报和滚动调度，具备了实时性；实现了宏观总控与局部调整相结合的调度目标，具备了调度决策不同时间尺度耦合嵌套特性，证明了系统的适用性。研究成果对于同类系统实时调度具有重要参考价值。

关键词：水库群；供、调水系统；入库径流预报； 需水预报；实时调度；调度流程；大连市
海洋重要水文参数的卫星遥感反演研究综述

吴桂平^1,2，刘元波^1,2

（1.中国科学院南京地理与湖泊研究所，江苏 南京 210008；2.中国科学院流域地理学重点实验室，江苏 南京 210008 ）

摘要：经过近50年的发展，海洋卫星遥感日益成熟，各种反演算法和数据集应运而生，大大推动了海洋水文过程的研究。围绕与海洋水文相关的重要物理参数，简要回顾了海面高度、海水深度、海表洋流及海水质量变化的卫星遥感反演基本原理、主要算法以及存在的难点和前沿性研究问题，介绍了应用遥感反演算法制成的全球海洋水文参数主要数据集，并结合目前存在的问题探讨了卫星遥感反演海洋水文参数的发展趋势及所面临的挑战。研究发现，运用多卫星、多通道、多模式的联合反演方法可以显著地提升人们监测海洋水文参数时空动态变化的能力，必将成为未来全球海洋水文遥感的主要研究方向之一。

关键词：海洋水文参数；反演算法； 光学遥感； 微波遥感； 多传感器联合反演；全球数据集
计算河流动力学理论体系框架探讨

黑鹏飞¹，假冬冬²，冶运涛³，周 刚⁴，胡德超⁵

（1. 中央民族大学生命与环境科学学院,北京 100081；2. 南京水利科学研究院水文水资源与水利工程科学国家重点实验室, 江苏 南京 210029；3. 中国水利水电科学研究院流域水循环模拟与调控国家重点实验室, 北京 100038；4. 中国环境科学研究院国家环境保护河口与海岸带环境重点实验室, 北京 100012； 5. 长江科学院河流研究所, 湖北 武汉 430010）

摘要：计算河流动力学重在研究河流动力的时空变化，实质上是计算流体力学、计算土力学在河流动力学中的应用，但并非简单组合，主要难点源自相间作用机理复杂性及研究对象多尺度性所导致的应用可行性问题。问题解决的过程中产生大量新的适于河流动力计算的数学模型和数值模型，但当前理论体系还缺乏系统性。综合河流动力学、计算流体力学、计算土力学理论，给出计算河流动力学的理论体系框架，明确计算河流动力学各部分内容的基本联系。基于此理论体系，可以方便地判断模型准确性、计算量和复杂性，构建最优的数学模型，并判断数学模型的发展方向。

关键词：计算河流动力学；河流动力学；数学模型；河流；泥沙