水库拦沙对长江水沙态势变化的影响
王延贵1，史红玲1，刘 茜2
（1. 国际泥沙研究培训中心，北京 100048；2. 北京交通大学，北京 100044）
摘要：为探讨水库建设对水沙变化态势的影响，通过引入来沙系数、输沙模数和流域水库调控系数，深入分析了长江水库拦沙对水沙态势变异的作用。长江干支流水沙参数与流域水库库容参数之间有重要关系，主要水文站年输沙量随流域累积库容增加而逐渐减少，来沙系数和输沙模数随流域水库调控系数成指数关系衰减，表明水库蓄水拦沙是长江水沙态势变异的主要因素。长江干支流输沙量衰减规律有一定差异，支流来沙系数和输沙模数衰减规律相对独立，而干流上游站衰减幅度大于下游水文站，且干流各站输沙模数的衰减规律比较接近。
关键词：长江水库拦沙；水沙态势变异；流域水库调控系数；来沙系数；输沙模数
江心洲平衡形态水动力条件的理论分析
刘晓芳，黄河清，邓彩云
（中国科学院地理科学与资源研究所陆地水循环及地表过程重点实验室，北京 100101）
摘要：江心洲在何种水动力条件下具有稳定的形态一直是河床演变学中的核心问题，采用河道形态变分分析方法，综合分析分汊型河流在水流连续方程、阻力方程、推移质输沙方程和总流能量方程控制下的自动调整机理。以长江中下游两汊不对称的江心洲为研究对象，将其形态概化为规则的等腰三角形，详细数理分析发现，主汊的分沙比大于分流比是这类江心洲达到平衡的必要条件。当分汊河道达到平衡状态时，江心洲的形态会随着汊道间水沙分配比例的调整而变化，且分流比与分沙比越接近，江心洲越窄长。另外，利用长江中游监利汊道的实测水沙资料对理论分析进行验证，结果表明，江心洲形态的理论计算值与相应的实测值符合程度较好，相对误差为9.0%左右。
关键词：分汊河流；分流比；分沙比；变分方法；江心洲
未来气候情景下长江上游区域积雪时空变化分析
——基于CMIP5多模式集合数据
陆桂华1,2，杨 烨1,2，吴志勇1,2，何 海1,2，肖 恒1,2
（1.河海大学水文水资源学院，江苏 南京 210098；2. 河海大学水问题研究所，江苏 南京210098）
摘要：基于21个CMIP5全球气候模式集合数据，耦合VIC模型，预估了未来30年（2011—2040年）RCP2.6、RCP4.5和RCP8.5 三个情景下长江上游区域积雪的时空变化。结果表明：与基准期1970—1999年相比，长江上游区域未来30年的多年平均气温和各月平均气温都将升高1~2℃，其中冬季和春季升温较大；平均年降水量将增加3%~4%，但秋、冬季降水有所减小。未来30年平均积雪深相对于基准期将减小37.8%左右，在积雪过程中达到最大积雪深的时间与基准期基本相同，而融雪开始的时间略有延后；从空间变化来看，冬季（1月份）长江上游区域大部分地区的积雪深都呈现减小趋势，部分地区积雪深减小超过了50%。
关键词：气候情景；耦合模式比较计划；长江上游区域；VIC模型；积雪；时空变化
1961—2010年中国400mm和800mm等雨量线时空变化
袁 喆，严登华，杨志勇，尹 军，袁 勇
（中国水利水电科学研究院流域水循环模拟与调控国家重点实验室，北京 100038）
摘要：以中国1961—2010年逐日降水数据为基础，利用等值线加权平均位置的方式对400mm和800mm等雨量线的空间位置进行定量化处理。采用Kendall 秩次相关检验法、Morlet小波分析法、滑动t-检验法和Yamamoto法时间序列分析方法，对全国50年来400mm和800mm等雨量线的空间位移特征进行了系统揭示：1961—2010年期间，400mm和800mm等雨量线加权平均位置坐标分别为（106°07′12″E，39°25′13″N）和（110°16′31″E，34°12′04″N）。趋势性结果表明，全国400 mm和800mm等雨量线有向西和向南发生移动的趋势，其中400mm等雨量线向西移动明显，800mm等雨量线向南移动明显；周期性结果表明，全国400mm等雨量线在经向和纬向上发生迁移的主周期分别为9年和12年，800mm等雨量线在经向和纬向变化均存在7年的主周期；突变性结果表明，400mm等雨量线空间位置的突变年为1995—1996年（纬度），800mm等雨量线空间位置的突变年为1975年、2002年（经度）和1980年、1982年和1987年（纬度）。
关键词：等雨量线；加权平均位置；时间序列分析；中国
淮河流域近500年洪旱事件演变特征分析
杨传国1,2，陈 喜1,2，张润润1,2，胡 琪3，余钟波1,2，郝振纯1,2，林朝晖4
（1. 河海大学水文水资源与水利工程科学国家重点实验室，江苏 南京 210098; 2. 河海大学水文水资源学院，江苏 南京 210098; 3. School of Natural Resources and Department of Earth and Atmospheric Sciences, University of Nebraska Lincoln, Lincoln, NE685880987, USA;4. 中国科学院大气物理研究所，北京 100029）
摘要：为了认识淮河流域过去500年洪旱事件发生规律并鉴别当前的洪旱情势，收集并对比分析了流域实测降雨资料、重建历史雨季降雨资料、历史旱涝等级资料、历史洪旱文献记录和历史调查洪水资料等多源洪旱灾害数据。以重建历史雨季降雨资料和历史旱涝等级资料为主要依据，通过滑动平均、频率计算、小波分析和突变检验等方法，分析流域过去500年洪水干旱时空分布特征和演变规律。结果表明，17世纪淮河流域洪旱灾害最严重，但20世纪极端洪旱事件发生频次最多。淮河流域洪旱事件存在40年左右的稳定长周期，主周期从18世纪的15~20年逐渐减少到19世纪的5年周期，近20年来出现2~3年的主周期，洪旱灾害事件呈增加趋势，流域社会经济发展面临着严峻的洪旱灾害威胁。
关键词：洪水；干旱；重建历史降雨；演变特征；淮河流域
塔里木河下游胡杨林耗水数值模型
苏里坦1，关东海2，王兴勇3，赵天宇1
（1．中国科学院新疆生态与地理研究所荒漠与绿洲生态国家重点实验室，新疆 乌鲁木齐830011；2．新疆水利水电科学研究院，新疆 乌鲁木齐 830049；3．中国水利水电科学研究院水力学研究所，北京 100038）
摘要：基于SHAW（Simultaneous Heat and Water）模型，以基本观测要素、植被参数和土壤剖面水热观测数据为模型的输入，对河岸胡杨林的耗水过程、土壤剖面水分变化和能通量进行了较小时间尺度上的模拟研究。结果表明，采用SHAW模型模拟的胡杨耗水量与观测值间存在较大偏差。因此，为了进一步提升水热耦合SHAW模型在干旱区的实用性，引入了地下水位因子GSI（Groundwater-Soil water Interaction），建立了改进的SHAW（GSI-SHAW）模型，解决干旱区荒漠河岸林耗水过程模拟的方法问题。采用SHAW模型和GSI-SHAW模型对胡杨耗水量的模拟进行了对比研究。结果显示，SHAW模型和GSI-SHAW模型模拟的胡杨耗水量与观测值的相关性系数分别为0.8533、0.9075，其平均相对误差分别为21.4%、16.9%，可见，改进的SHAW模型的模拟值更加接近试验观测值。地下水位的考虑一定程度上提升了传统SHAW模型的模拟精度，为干旱区自然植被耗水量的计算提供了新的方法和科学依据。
关键词：胡杨林；耗水；改进SHAW模型；塔里木河下游；中国西北干旱区
近期新疆哈密代表性冰川变化及对水资源影响
王璞玉1，李忠勤1,2，周 平1，金 爽1，陈 辉2
（1. 中国科学院寒区旱区环境与工程研究所冰冻圈科学国家重点实验室/天山冰川站，甘肃 兰州 730000；2. 西北师范大学地理与环境科学学院，甘肃 兰州 730070）
摘要：新疆哈密是资源性缺水地区，冰川是该区主要的供给水源。为了对该区冰川和水文水资源的现状和未来变化做出合理评价和预估，以榆树沟6号冰川和庙尔沟冰帽为代表，结合野外实测资料以及近年来在冰川变化领域所取得的相关研究成果，综合分析了在气候变暖背景下新疆哈密地区冰川近年来的变化过程及对水资源的影响。结果表明，1972—2011年，榆树沟6号冰川厚度平均减薄20m，减薄速率约为0.51m/a，冰川末端退缩254m，年均退缩约6.5m。由于强烈消融，冰川表面径流发育明显，且在末端有冰碛湖形成。相比之下，庙尔沟冰帽消融速率相对较小，对气候变暖的敏感度较低，这与冰帽类型和所处海拔较高有直接关系。综合分析发现，流域有无冰川覆盖及覆盖的比率大小，直接造成哈密地区不同流域近期河川径流变化的差异。
关键词：冰川变化；水资源；榆树沟6号冰川；庙尔沟冰帽；哈密地区
基于人工降雨模拟试验的坡面水文连通性
王盛萍1，姚安坤2，赵小婵1
（1. 华北电力大学资源与环境研究院区域能源环境系统优化教育部重点实验室，北京102206；2. 北京林业大学水土保持学院，北京 100083）
摘要：为揭示降水入渗时坡面水文连通性特征以及坡度和雨强等对水文连通性的影响，研究基于人工降雨模拟试验，分别采用结构性指标径流长度（Flowlength，lF）和功能性指标相对地表连接函数（Relative Surface Connection Function, FRC）分析不同坡度坡面（5°、15°、25°；面积为3750px×1250px）在不同雨强（25mm/h、50mm/h、70mm/h和86mm/h）条件下的水文连通性。研究结果表明：坡面水文连通结构主要由lF为0~100mm的径流路径组成，其栅格频率高达90%以上，而lF为100~500mm的栅格频率为3%~8%。水文连通性的迅速发展主要集中在降水事件的开始阶段，而一定阶段以后入渗速率与地表储水速率基本持平，水文连通性呈现出较为稳定的发展趋势。坡度增大整体上有利于水文连通性的发展。5°坡面对应水文连通面积比率整体小于其他坡度。由于雨强增大到一定范围时将增加地表糙率，从而阻碍水文连通性的发展，研究中仅25mm/h对应水文连通性明显区别于其他雨强，而其他雨强间则在多数情况下水文连通性差别较小。
关键词：水文连通性；降雨模拟；相对地表连接函数；径流长度
闸控河网水文-水动力-水质耦合数学模型
——I.理论
陈炼钢1，施 勇1，钱 新2，栾震宇1，金 秋1
（1. 南京水利科学研究院水文水资源与水利工程科学国家重点实验室，江苏 南京 210029；
2. 南京大学环境学院污染控制与资源化研究国家重点实验室，江苏 南京 210023）
摘要：针对流域级闸控大型河网水环境日常管理的应用需求与实际特点，将资料适应性强的水文学方法与数据要求相对严苛的水动力-水质数值模型相融合，构建了一维与二维嵌套、分块组合的闸控大型河网水文-水动力-水质耦合数学模型DHQM（Hydrology, Hydrodynamics, and Water Quality Model for Impounded Rivers）。模型由河道径流模拟、闸坝调度过程模拟、河道水质模拟、区间入流及入河污染负荷估算和水质预警实时校正等5个模块组成。模型可服务于水环境实时预警和调度，也可为闸坝水文环境效应的量化提供基础工具。
关键词：环境水利；水量水质耦合模拟；模型集成；多闸坝河网
双时间步法在二维浅水方程求解中的应用
喻海军1，黄国如1,2，武传号1
（1. 华南理工大学土木与交通学院，广东 广州 510640；2. 华南理工大学亚热带建筑科学国家重点实验室，广东 广州 510640）
摘要：针对传统算法效率低的问题，将隐式双时间步法应用于求解二维浅水方程，建立了非结构网格下高效的有限体积模型。在应用双时间步法时，虚拟时间层中的定常问题采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法进行迭代求解。通过模拟计算4个典型算例以及与传统显式算法进行比较，对模型精度、效率及处理实际问题能力进行检验，分析了时间步长、内迭代次数对模型性能的影响。结果表明，双时间步法放宽了稳定性对时间步长的限制，时间步长可取到显式格式10倍以上，计算耗时减少了50%以上，模型具有良好精度与适应性，具有较好的推广应用价值。
关键词：双时间步法；二维浅水方程；有限体积法；非结构网格
感潮河网二维水流-输运耦合数学模型
宋利祥1,2，杨 芳1,2，胡晓张1,2，杨莉玲1,2，周建中3
（1. 珠江水利科学研究院，广东 广州 510611; 2. 水利部珠江河口动力学及伴生过程调控重点实验室，广东 广州 510611;3. 华中科技大学水电与数字化工程学院，湖北 武汉 430074）
摘要：感潮河网水流不仅受径流、潮汐动力的双重影响，还常受水闸调度等人类活动的影响，水动力条件极为复杂。针对水闸调度影响下感潮河网复杂水动力及其伴随物质输移扩散模拟，提出了考虑水闸调度并集成输运对流项的水流-输运通量耦合求解器，建立了基于非结构网格Godunov格式的二维水流-输运耦合数学模型。采用具有时空二阶精度的MUSCL-Hancock预测-校正格式，结合变量重构限制器技术，在保证计算精度的同时避免了数值振荡。运用斜底三角单元网格表达水闸线状地形，并通过地形调整模拟了水闸启闭过程。算例研究表明，该模型计算精度高，可有效模拟水闸调度影响下感潮河网水流运动及污染物输运过程，具有较好的推广应用价值。
关键词：物质输移扩散；感潮河网；水闸调度；Godunov格式；非结构网格
主槽一维和滩地二维侧向耦合洪水演进模型
张防修1,2，韩龙喜1，王 明2，黄 瑞1
（1. 河海大学环境科学与工程学院，江苏 南京 210098；2. 黄河水利科学研究院，河南 郑州 450003）
摘要：为实现对复式河道大尺度漫滩洪水的快速模拟，结合洪水在滩槽行洪特点和河道一维、二维水动力学模型研究成果，在主槽与滩地分别划分为矩形和三角形网格，以左右滩唇为分界实现两网格系统嵌套；建立主槽一维水动力学数学模型，内插大断面间网格处的水力因子，在滩地采用简化的二维水动力学模型模拟滩槽间、滩地网格间水量交换，实现主槽一维与滩地二维模型侧向耦合。利用黄河下游“96·8”大洪水实测资料，对花园口至夹河滩河段洪水演进过程验证，结果表明该方法模拟精度较高，能够利用该模型模拟复式河道大漫滩洪水过程。
关键词：复式河道；漫滩洪水；数学模型；侧向耦合
荆江段河岸粘性土体抗冲特性试验
宗全利1,2，夏军强1，张 翼1，许全喜3
（1. 武汉大学水资源与水电工程科学国家重点实验室，湖北 武汉 430072；2. 石河子大学水利建筑工程学院，新疆 石河子 832000；3. 长江水利委员会水文局，湖北 武汉 430010）
摘要：粘性土体具有较强的抗冲能力，其抗冲性强弱决定了河岸的侵蚀速率。为了估算河岸的侵蚀速率，确定粘性土体的起动条件、侵蚀系数及两者之间的数量关系非常重要，为此开展了荆江段河岸粘性土体起动条件与冲刷特性的水槽试验研究。根据试验结果获得了粘性土体的起动流速与土体液限/自然含水率之间关系以及起动切应力与干密度、起动切应力与液性指数之间的定量表达式，综合反映了粘性土起动条件与其物理性质指标之间的关系；结合冲刷特性试验结果，还获得了土体侵蚀系数随起动切应力的变化规律，并与已有其他试验结果进行了比较。结果表明：荆江河岸粘性土侵蚀系数均比相同条件下其他公式计算值偏大，这与该试验土体中粘粒含量相对较低、土体结构受到扰动等因素有关；根据试验结果，拟合得到了荆江河岸粘性土侵蚀系数与起动切应力之间定量关系式，其相关系数R2=0.90，故该关系式能为荆江段崩岸过程的计算提供参考依据。
关键词：粘性土；起动切应力；侵蚀系数；抗冲特性；荆江河岸
梯级水库设计洪水最可能地区组成法计算通式
刘章君，郭生练，李天元，徐长江
（武汉大学水资源与水电工程科学国家重点实验室，水资源安全保障湖北省协同创新中心，湖北 武汉 430072）
摘要：利用Copula函数建立各分区洪水的联合分布，基于联合概率密度最大原则，推导得到最可能地区组成法的计算通式，并用来推求梯级水库下游断面的设计洪水。选择清江流域水布垭—隔河岩—高坝洲梯级水库为例，开展了验证和方法比较研究。结果表明：最可能地区组成法计算得到的设计洪水值位于同频率地区组成法多方案计算结果的区间之内；受清江梯级水库调洪的影响，宜都断面设计洪水的削峰率十分显著，最可能地区组成法推求100年一遇设计洪水的削峰率达到30.2%。该法具有较强的统计基础，组成方案唯一，结果合理可行，为复杂梯级水库设计洪水的计算提供了一种新途径。
关键词：设计洪水；梯级水库；最可能地区组成法；计算通式；Copula函数；清江流域
不同微尺度膜下滴灌棉田土壤水盐空间变异特性
史文娟1，马 媛1，徐 飞1，王全九1,2
（1. 西安理工大学西北水资源与环境生态教育部重点实验室，陕西 西安 710048；2. 黄土高原土壤侵蚀与旱地农业国家重点实验室，陕西 杨凌 712100）
摘要：为揭示微尺度膜下滴灌农田土壤水盐的空间变异特性，通过大田试验，采用经典统计和地质统计方法，研究了3种微尺度（0.25m、1m和4m）和不同土层深度棉田土壤水盐的空间变异性，并确定其合理取样数。结果发现，3种微尺度下，土壤含水量的变异强度为中等偏弱，其变异性随尺度增大而增强，随深度增加呈先增强后减弱趋势。土壤含盐量的变异强度为中等偏强，其变异性随尺度增大而增强，随深度增大呈先减弱后增强趋势；3种尺度和不同深度条件下，土壤含水量和含盐量的半方差函数大部分可采用高斯模型模拟，且精度较高；样品的合理取样数为367个。研究结果可为制定膜下滴灌土壤水盐的监测方案和调控措施提供理论指导。
关键词：微尺度；膜下滴灌；土壤水盐；空间变异
变化环境下城市水文学的发展与挑战
——I.城市水文效应
张建云1,2，宋晓猛1,2，王国庆1,2，贺瑞敏1,2，王小军1,2
（1. 南京水利科学研究院水文水资源与水利工程科学国家重点实验室，江苏 南京 210029；2. 水利部应对气候变化研究中心，江苏 南京 210029 ）
摘要：全球气候变化和快速城市化导致的城市水循环过程变化是当前城市水文学研究的热点问题。为综合理解城市水循环演变过程，结合城市水文学的发展历程，剖析水循环过程对快速城市化进程的响应机制。总结了国内外城市化水文效应的主要成果，包括城市化对水循环过程、洪涝灾害、水生态系统以及水资源的影响。系统归纳了城市化水文效应的评估方法和技术手段。针对现有研究中的不足，指出变化环境下城市化水文效应研究面临的主要挑战及关键技术难题，提出未来研究的重点方向，如城市化降水效应的机理、不透水面的分布及有效性评估、城市化与水生态系统的响应关系与综合城市水资源管理及需水预测等。
关键词：城市水文学；水循环；水环境；水生态；水资源；水文模型
中国泥石流起动物理模拟试验研究进展
倪化勇1,2，唐 川1
（1. 成都理工大学地质灾害防治与地质环境保护国家重点实验室，四川 成都 610059；2. 成都地质矿产研究所，四川 成都 610081）
摘要：为进一步提高泥石流起动物理模拟试验的科学性，完善试验体系，采用资料收集与分析方法，总结了近年来中国泥石流起动物理模拟试验开展的现状，从水流冲刷与泥石流起动试验以及人工降雨与泥石流起动试验两个方面论述了中国泥石流起动试验取得的主要进展和理论成果。在国内外泥石流起动物理模拟试验对照基础上，提出中国泥石流起动物理模拟试验研究的建议：提高水流浓度、降雨雨型与土体特征的相似率；加强降雨或水流作用下土体物理力学特征变化与泥石流起动响应研究；推动降雨与汇流共同作用下泥石流起动综合性物理模拟试验的开展。
关键词：泥石流；起动；物理模拟试验；水流冲刷；人工降雨
