美国科罗拉多河水沙变化分析_吴保生

美国科罗拉多河流域水资源开发与利用研究科罗拉多河是美国西部最重要的河流之一，流经科罗拉多州、犹他州、亚利桑那州以及其他州，被誉为“美国的生命线”。

这条河流不仅给数百万人提供生活用水，还支持着农业、工业和能源生产等多个经济部门。

然而，随着人口增长和气候变化等因素的影响，科罗拉多河流域的水资源面临着日益紧张的局势，因此对此进行深入的研究是十分必要的。

首先，科罗拉多河流域的水资源开发与利用面临的最大挑战之一是干旱。

科罗拉多州以及周边地区的干旱问题长期存在，并对水资源的供应和分配带来了极大的压力。

据统计，科罗拉多河流域地区每年平均只有不到15英寸的降水量，远低于全美平均水平。

因此，科罗拉多河流域的居民和农田主要依靠河水和地下水来满足他们的生活和农业用水需求。

其次，大规模水坝建设和水库管理也是科罗拉多河流域水资源开发与利用的重点。

科罗拉多河流域的水坝和水库被用于调节河流流量，提供优质的灌溉水，以及发电等目的。

其中，科罗拉多河上最知名的水坝是胡佛大坝，这是世界上最大的重力混凝土坝。

胡佛大坝不仅为该地区提供了丰富的水资源，还成为了一个重要的旅游景点，吸引了大量的游客。

尽管已经有了大规模的水坝和水库，但科罗拉多河流域的水资源依然十分有限。

因此，科学合理的水资源管理和节约措施是至关重要的。

一方面，需要改善现有水资源的利用效率，减少水的浪费。

例如，采用节水型灌溉技术和设备，合理安排农业灌溉的时间和量。

另一方面，通过水资源的再利用和回收，可以有效增加可用水量。

科罗拉多河流域的某些城市和农田已经开始采用回用废水来供应非饮用水的需求。

此外，生态环境保护和河流生态系统的恢复也应该成为科罗拉多河流域水资源开发与利用的关注点。

科罗拉多河是一个重要的生态系统，支持着众多物种的繁衍生息。

然而，过度的水资源开发和使用对河流生态系统造成了巨大的压力。

为了保护河流的生态功能，需要制定并执行相应的环境保护政策，限制盲目的水资源开发，并进行河道的生态修复工作。

第３２卷第３期２０２１年５月㊀㊀水科学进展ＡＤＶＡＮＣＥＳＩＮＷＡＴＥＲＳＣＩＥＮＣＥＶｏｌ.３２ꎬＮｏ.３Ｍａｙ２０２１ＤＯＩ:１０ １４０４２/ｊ ｃｎｋｉ ３２ １３０９ ２０２１ ０３ ００９三峡库区汛期洪峰和沙峰异步运动特性的三维数值模拟张帮稳ꎬ吴保生ꎬ章若茵(清华大学水沙科学与水利水电工程国家重点实验室ꎬ北京㊀１０００８４)摘要:三峡水库蓄水后水深的增加导致洪水传播过程中沙峰滞后于洪峰的时间增加ꎬ利用洪峰和沙峰异步运动特性进行沙峰排沙调度是减缓水库淤积的重要手段之一ꎮ采用三维数值模型ＳＣＨＩＳＭ(Ｓｅｍｉ￣ｉｍｐｌｉｃｉｔＣｒｏｓｓ￣ｓｃａｌｅＨｙｄｒｏ￣ｓｃｉｅｎｃｅＩｎｔｅｇｒａｔｅｄＳｙｓｔｅｍＭｏｄｅｌ)研究三峡库区洪峰和沙峰的异步运动规律ꎬ基于２０１３年的汛期水沙传播过程验证模型的可靠性和准确性ꎬ初步分析不同洪水位下洪峰和沙峰的异步运动特性ꎮ结果表明:数值模型能够准确地模拟三峡水库汛期洪峰和沙峰长距离的传播过程ꎻ坝前蓄水位对洪峰传播时间的影响不大ꎬ但是对沙峰传播时间的影响较为显著ꎻ坝前水位增加使水流的流速减小和挟沙能力降低ꎬ沙峰传播速度减慢且峰值沿程不断降低ꎬ最终导致库区中沙峰滞后于洪峰出现的时间沿程越来越长ꎮ研究结果可为进一步优化三峡水库的沙峰排沙调度等提供科学依据ꎮ关键词:沙峰ꎻ洪峰ꎻ滞后时间ꎻ含沙量ꎻ三峡水库中图分类号:ＴＶ１４５㊀㊀㊀文献标志码:Ａ㊀㊀㊀文章编号:１００１￣６７９１(２０２１)０３￣０４０８￣１０收稿日期:２０２０￣０４￣２０ꎻ网络出版日期:２０２０￣１２￣０３网络出版地址:ｈｔｔｐｓ:ʊｋｎｓ.ｃｎｋｉ.ｎｅｔ/ｋｃｍｓ/ｄｅｔａｉｌ/３２.１３０９.Ｐ.２０２０１２０３.１４５４.００２.ｈｔｍｌ基金项目:国家重点研发计划资助项目(２０１７ＹＦＣ０４０５２０２)ꎻ国家自然科学基金资助项目(５１６３９００５)作者简介:张帮稳(１９８７ )ꎬ男ꎬ山东菏泽人ꎬ助理研究员ꎬ博士后ꎬ主要从事河流动力学方面研究ꎮＥ￣ｍａｉｌ:ｚｈａｎｇｂａｎｇｗｅｎ＠ｍａｉｌ.ｔｓｉｎｇｈｕａ.ｅｄｕ.ｃｎ通信作者:吴保生ꎬＥ￣ｍａｉｌ:ｂａｏｓｈｅｎｇ＠ｔｓｉｎｇｈｕａ.ｅｄｕ.ｃｎ三峡水库泥沙淤积问题不仅关系到水库的使用寿命和功能发挥ꎬ而且对水库及下游的生态环境具有重要的影响[１]ꎮ近年来长江上游梯级水库的建设导致三峡入库径流量和泥沙量发生变化ꎬ严重影响了三峡水库综合效益的发挥ꎬ为此ꎬ三峡水库在 蓄清排浑 运用方式的基础上ꎬ开展了中小洪水调度和汛末提前蓄水调度ꎬ目的是提高发电和航运效益ꎬ减轻坝下游的防洪压力[２]ꎮ但汛期中小洪水调度和汛末提前蓄水会抬高库区平均水位ꎬ降低库区流速ꎬ加重库区泥沙淤积ꎮ为缓解中小洪水调度及汛末提前蓄水调度带来的泥沙淤积影响ꎬ三峡水库于２０１２年和２０１３年汛期利用洪峰和沙峰的异步运动特性开展了沙峰调度试验ꎬ取得了较好的排沙效果[３￣４]ꎮ研究者基于实测资料对三峡库区洪峰和沙峰的运动特性进行了分析ꎬ发现三峡水库蓄水后ꎬ库水位抬高ꎬ库区水流流速减慢ꎬ沙峰输移时间较蓄水前大幅增加ꎬ沙峰多滞后于洪峰ꎬ且沙峰滞后于洪峰的时间沿程逐渐增加ꎬ有利于进行水库的沙峰调度[５￣７]ꎮ但目前关于洪峰和沙峰异步运动特性多以实测资料进行分析ꎬ基于试验研究及数值模拟开展相对较少ꎬ对洪峰和沙峰异步运动特性缺乏较为深入的认识ꎮ利用数值模拟方法对不同情况下洪峰和沙峰运动特性进行系统分析ꎬ可为应对未来愈加复杂的水沙变化和地形变化情势下水库的沙峰排沙调度提供参考ꎮ水沙数学模型在研究水沙输移㊁泥沙淤积及河床演变规律中有着重要的作用ꎮ基于断面平均或垂向平均的一维和二维数值模型[８￣９]能够对工程水沙问题进行一定程度的研究ꎬ但是泥沙淤积㊁冲刷和输移是一个非常复杂的水动力现象ꎬ特别是在复杂的自然地形ꎬ例如河漫滩㊁弯曲河段和沙波地形等产生二次流和流体分离等三维水动力特性ꎬ严重影响了数值模拟结果的准确性ꎮ随着计算机科学技术和高性能计算机的发展ꎬ三维模型[１０￣１１]逐渐应用到工程实际当中ꎬ目前大多数采用三维数值模拟方法针对局部河段或库区的泥沙输运和淤积进行研究ꎬ但对洪水传播过程中洪峰和沙峰异步运动特性研究较少ꎮ主㊀第３期张帮稳ꎬ等:三峡库区汛期洪峰和沙峰异步运动特性的三维数值模拟４０９㊀要原因是洪峰和沙峰的运动特性与研究河流或水库的尺度(时间尺度和空间尺度)有关ꎬ较短的距离和较短的时间难以显示出洪峰和沙峰运动规律ꎮ三峡库区地形较为复杂ꎬ总长超６６０ｋｍꎬ对此进行长距离洪水演进的三维数值模拟具有较大的难度ꎮ掌握洪峰与沙峰异步运动特性㊁开展入库泥沙实时监测与预报ꎬ是进行沙峰排沙调度的基础ꎮ为进一步提高沙峰预报精度ꎬ掌握更精准的沙峰调度时机㊁制订更科学合理的沙峰排沙调度方案ꎬ本文采用三维水沙数值模型ＳＣＨＩＳＭ(Ｓｅｍｉ￣ｉｍｐｌｉｃｉｔＣｒｏｓｓ￣ｓｃａｌｅＨｙｄｒｏｓｃｉｅｎｃｅＩｎｔｅｇｒａｔｅｄＳｙｓｔｅｍＭｏｄｅｌ)对三峡库区汛期洪水传播过程中洪峰和沙峰运动进行模拟ꎬ初步分析三峡水库不同蓄水位下洪峰和沙峰在万县 三峡大坝库区的异步运动特性ꎬ可为三峡水库的沙峰调度提供参考依据ꎮ１㊀数值模拟方法１.１㊀水流的基本方程ＳＣＨＩＳＭ水动力学模型[１２￣１３]的控制方程采用基于Ｒｅｙｎｏｌｄｓ时均Ｎ￣Ｓ方程ꎬ满足静压假定和Ｂｏｕｓｓｉｎｅｓｑ涡黏性假定ꎮ在笛卡尔坐标系下ꎬ对于不可压缩流体Ｎ￣Ｓ方程的连续性方程为∂ｕ∂ｘ＋∂ｖ∂ｙ＋∂ｗ∂ｚ＝０(１)动量守恒方程:ｄｕｄｔ＝ｆｕ－ｇ∂η∂ｘ－１ρ０∂ｐａ∂ｘ－ｇρ０ʏηｚｂ∂ρ∂ｘｄｚ＋∂∂ｚｋｍｖ∂ｕ∂ｚ()＋ｋｍｈ∂２ｕ∂ｘ２＋∂２ｕ∂ｙ２()(２)ｄｖｄｔ＝ｆｖ－ｇ∂η∂ｙ－１ρ０∂ｐａ∂ｙ－ｇρ０ʏηｚｂ∂ρ∂ｙｄｚ＋∂∂ｚｋｍｖ∂ｖ∂ｚ()＋ｋｍｈ∂２ｖ∂ｘ２＋∂２ｖ∂ｙ２()(３)式中:ｘ和ｙ分别表示笛卡尔水平坐标ꎬｚ为垂向坐标ꎬ向上为正ꎻｕꎬｖꎬｗ分别表示３个方向的流速ꎻｔ为时间ꎻｆ为柯氏力系数ꎻη为自由水面ꎻｚｂ为河床底高程ꎻρ０和ρ分别表示参考密度和混合流体的密度ꎻｇ为重力加速度ꎻＫｍｈ和Ｋｍｖ分别为水平与垂直涡黏性系数ꎬ其中垂向涡黏性系数根据紊流模型进行封闭ꎬ水平涡黏性系数采用常数化处理ꎻｐａ为自由水面大气压强ꎮ自由水面采用水位函数法处理ꎬ对连续方程(１)沿水深方向积分ꎬ可得自由水面方程∂η∂ｔ＋∂∂ｘʏηｚｂｕｄｚ＋∂∂ｙʏηｚｂｖｄｚ＝０(４)模型在河床底部的动力学边界条件由底床摩擦剪应力和底层水体的雷诺应力平衡给出Ｋｍｖ∂ｕ∂ｚꎬ∂ｖ∂ｚ()＝τｂｘꎬτｂｙ()㊀㊀ｚ＝ｚｂ(５)式中:ｔｂｘ和ｔｂｙ分别为床面的ｘ㊁ｙ方向摩擦剪应力ꎮ对于垂向紊动涡黏性系数Ｋｍｖꎬ利用紊流闭合模型ＧＬＳ(ＧｅｎｅｒｉｃＬｅｎｇｔｈＳｃａｌｅ)[１４]进行求解ꎬ包括紊动能ｋ方程和通用紊动长度ψ方程ꎮＧＬＳ紊流模型的控制方程为:ｄｋｄｔ＝∂∂ｚ(υψｋ)＋ＫｍｖＭ２＋μＮ２－ε(６)ｄψｄｔ＝∂∂ｚυψ∂ψ∂ｚ()＋ψｋ(ｃψ１υＭ２＋ｃψ３μＮ２－ｃψ２Ｆｗε)(７)υψｋ＝Ｋｍｖσｋ㊀㊀υψ＝Ｋｍｖσψ(８)Ｍ２＝∂ｕ∂ｚ()２＋∂ｖ∂ｚ()２㊀㊀㊀Ｎ２＝ｇρ０∂ρ∂ｚ(９)式中:ｋ为紊动能ꎻψ为通用紊动长度参数ꎻμ为盐度㊁温度等物质的垂向扩散系数ꎻυｋ和υｙ分别表示紊动４１０㊀水科学进展第３２卷㊀能和通用紊动长度的垂向扩散系数ꎻε为紊动耗散项ꎻＦｗ为壁函数ꎬ在ｋ￣ε模式中为１ꎻσｋ㊁σψ㊁ｃψ１㊁ｃψ２和ｃψ３为模型系数ꎻＭ２和Ｎ２分别表示由于剪切变形和密度分层而引起的紊动能产生项ꎬ通用紊动长度ψ和紊动耗散项ε作为紊流模型中的关键参量ꎬ其表达式如下:ψ＝(ｃμ)ｍｋｎｌｐꎬε＝(ｃμ)ｍｋｎｌｐ(１０)式中:ｌ为紊动掺混长度ꎻｃμ为常数０.３ꎮ在ＧＬＳ模型ｋ￣ε模式双方程紊流模式的参数取值见表１ꎮ表１㊀ＧＬＳ模型ｋ￣ε模式双方程紊流参数值Ｔａｂｌｅ１Ｐａｒａｍｅｔｅｒｖａｌｕｅｓｏｆｔｈｅｔｕｒｂｕｌｅｎｃｅｉｎｋ￣εｅｑｕａｔｉｏｎｂａｓｅｄｏｎＧＬＳｍｏｄｅｌ参数ｍｎｐσｋσψｃψ１ｃψ２ｃψ３数值３１.５－１１１.３１.４４１.９２１１.２㊀悬移质泥沙输移基本方程在悬移质泥沙动力学中ꎬ对流扩散理论将泥沙视为单一的连续介质ꎬ并假设悬移质泥沙运动与水流运动在垂直方向上存在速度差ꎬ且等于泥沙颗粒的沉降速度ꎮ基于对流￣扩散理论的三维悬移质泥沙输运方程表达式为∂Ｃｉ∂ｔ＋ｕ∂Ｃｉ∂ｘ＋ｖ∂Ｃｉ∂ｙ＋(ｗ＋ωｓꎬｉ)∂Ｃｉ∂ｚ＝∂∂ｚＫｓｖ∂Ｃｉ∂ｚæèçöø÷＋Ｋｓｈ∂２Ｃｉ∂ｘ２＋∂２Ｃｉ∂ｙ２æèçöø÷－Ｄｑ＋Ｅｑ(１１)式中:Ｃｉ为ｉ组含沙量ꎻｗｓꎬｉ为ｉ组泥沙颗粒的沉降速度ꎻＫｓｖ为泥沙垂向扩散系数ꎬ通常假定与水流紊动黏性系数呈倍数关系ꎬ可通过紊流模型求解或采用经验关系估计ꎬＫｓｖ＝Ｋｍｖ/σｓꎬσｓ为Ｓｃｈｍｉｄｔ数ꎬ通常取值在０.６~１.２之间ꎻＫｓｈ为泥沙水平扩散系数ꎬ考虑到水平扩散的量级远小于垂向扩散ꎬ通常忽略不计ꎻＤｑ为泥沙的沉积通量ꎻＥｑ为泥沙的冲刷通量ꎮ三峡水库中细颗粒泥沙存在絮凝现象ꎬ在泥沙运动力学中ꎬ采用絮凝因数Ｆ反映絮凝对泥沙颗粒的沉速的影响[１５]ꎬ其表达式为Ｆ＝ωｓꎬｉ/ω０(１２)式中:ω０为泥沙颗粒静水中的沉速ꎮ广泛采用的絮凝因数Ｆ＝αｄβꎬ其中ꎬα取值为０.０１３[１６]或０.００５５[１７]ꎬβ取值为－１.９[１５]或－１.０２[１６]ꎮ张地继等[１６]采用α＝０.００５５ꎬβ＝－１.０２研究了万县 庙河库区沙峰的衰减规律ꎬ本研究中采用相同的絮凝因数系数ꎮ１.３㊀河床冲刷演变基本方程根据泥沙质量守恒方程ꎬ悬移质泥沙输移过程中可将床面边界条件视为床面附近泥沙通量的处理ꎬ包括床面泥沙的沉积通量Ｄｂ和冲刷通量Ｅｂꎬ其表达式分别为:Ｄｂ＝ωｓꎬｉｃ１ꎬｉ(１３)Ｅｂ＝Ｅ０ꎬｉ(１－ｑｉ)(τｆ/τｃｒꎬｉ－１)㊀㊀τｆ>τｃｒꎬｉ(１４)式中:ｃ１ꎬｉ为数值模拟中最底部一层网格的ｉ组泥沙浓度ꎻＥ０ꎬｉ为经验冲刷率系数ꎬ取决于局部床面泥沙颗粒条件ꎬ取值大小范围为１０－４~１０－２ｋｇ/(ｍ２ ｓ－１)ꎻｑｉ为床面层ｉ组泥沙体积分数ꎻτｃｒꎬｉ为ｉ组泥沙颗粒临界起动剪切应力ꎻτｆ为水流底部的剪切应力ꎮ河床泥沙的冲淤和悬移质泥沙的净输移导致河床表面的地形变化ꎬ其公式为Δｈ＝(Ｄｑ－Ｅｑ)Δｔρｓ(１－ｐ)(１５)式中:Δｈ为床面高程变化量ꎻρｓ为泥沙颗粒的密度ꎮ２㊀数值模型的建立及验证为了研究三峡库区汛期洪水传播过程中洪峰和沙峰异步运动特性ꎬ本文选择万县站 三峡大坝的库区段㊀第３期张帮稳ꎬ等:三峡库区汛期洪峰和沙峰异步运动特性的三维数值模拟４１１㊀作为研究对象(图１)ꎮ若把寸滩或清溪场作为入口边界条件ꎬ则河段过长ꎬ计算耗时过大ꎬ计算效率极低ꎮ万县至大坝长２８０ｋｍꎬ不仅长度适中ꎬ并且也是目前少有的长距离㊁大范围实际水库的三维数值模拟ꎬ难度较高ꎻ另一方面原因ꎬ万县 三峡大坝库区之间仅有靠近大坝庙河站可以测量流量和含沙量ꎬ长度不足以模拟洪峰和沙峰的传播特性ꎬ因此只有万县站满足本文的研究要求ꎬ可作为入口条件ꎮ模拟中采用的资料分别来源于万县站㊁奉节站㊁巫山站㊁庙河站和茅坪站ꎬ如图１(ａ)所示ꎬ其中水文站有万县站和庙河站ꎬ具有日平均水位㊁流量和含沙量信息ꎬ水位站有奉节站㊁巫山站和茅坪站ꎮ图１㊀研究河段的地形高程㊁水文站和网格设置Ｆｉｇ.１Ｔｏｐｏｇｒａｐｈｉｃｅｌｅｖａｔｉｏｎꎬｈｙｄｒｏｌｏｇｉｃａｌｓｔａｔｉｏｎａｎｄｇｒｉｄｓｅｔｔｉｎｇｏｆｔｈｅｓｔｕｄｙｒｅａｃｈ２.１㊀模型网格设置万县 三峡大坝的地形较为复杂ꎬ宽阔和狭窄的河段交替变化ꎬ水平方向上进行网格设置时主要混合使用２种网格类型ꎬ在宽阔河段混合采用四边形和三角形网格ꎬ主河道采用四边形网格ꎬ岸滩采用三角形网格ꎬ在狭窄河段采用三角形网格ꎬ网格尺度约为２０~２３ｍꎬ共得到４９５６６２个网格节点和８５６０５６个网格单元ꎮ计算初始地形条件由２０１１年实测地形插值得到ꎬ插值后的局部地形如图１(ｂ)所示ꎮ此外ꎬ为了更好地模拟河道至坝前水深变化幅度大的特点并且提高计算效率ꎬ垂向上采用分层ＬＳＣ２坐标[１７]ꎬ最大水深处可达３６层ꎬ最浅处为１６层ꎬ平均约为１９层ꎬ每层厚度约５~６ｍꎮ图１(ｃ)为局部河段横断面的垂向网格示意图ꎮ为了数值模拟的稳定性ꎬ时间步长设置３０ｓꎮ整个计算模型在清华大学高性能计算集群上采用２８０个核进行计算ꎬ模拟时间为８ｄꎮ２.２㊀模型边界条件本文以三峡库区２０１３年汛期７月１日 ８月１日作为模拟时间段ꎬ万县站的流量和含沙量过程作为模型的入口边界条件ꎬ茅坪站的水位过程作为出口边界条件(图２)ꎮ三峡库区泥沙主要以悬移质输移为主ꎬ万县站７月实测的泥沙粒径分别为０.００２ｍｍ㊁０.００４ｍｍ㊁０.００８ｍｍ㊁０.０１６ｍｍ㊁０.０３２ｍｍ和０.０６４ｍｍꎬ对应的泥沙级配分别为１０.３％㊁１２.３％㊁２０.９％㊁２６.０％㊁１８.５％和１２.０％ꎬ考虑到泥沙０.００２ｍｍ和０.００４ｍｍ的粒径很小ꎬ统一归为０.００３ｍｍ的泥沙进行计算ꎬ床沙依据大断面实测的泥沙级配进行插值设置ꎬ并且选择与悬移质同样的５组粒径泥沙ꎮ数值模拟泥沙输移过程中暂时没有考虑温度和盐度的影响ꎮ４１２㊀水科学进展第３２卷㊀图２㊀数值模型边界条件Ｆｉｇ.２Ｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎｓｉｎｎｕｍｅｒｉｃａｌｍｏｄｅｌ２.３㊀数值模型验证为验证数值模型模拟洪水传播和泥沙输移过程的正确性ꎬ分别对比分析了实测与模拟的奉节站和庙河站的水位ꎬ以及庙河站的流量㊁平均流速和含沙量ꎮ图３对比了奉节站和庙河站水位的模拟值和实测值ꎬ可以看出数值模拟的结果和实测值吻合较好ꎬ能够较好地反映洪水传播过程中水位的变化ꎮ图４对比了庙河站断面平均的流量㊁流速和含沙量的模拟值和实测值ꎬ可以看出洪水流量模拟值的变化和实测值的变化过程基本一致ꎬ局部的差别可能是由于万县 三峡大坝图３㊀水位模拟值和实测值的对比结果Ｆｉｇ.３Ｃｏｍｐａｒｉｓｏｎｏｆｗａｔｅｒｌｅｖｅｌｂｅｔｗｅｅｎｓｉｍｕｌａｔｅｄａｎｄｍｅａｓｕｒｅｄｒｅｓｕｌｔｓ图４㊀庙河站平均流量㊁流速和平均含沙量的模拟值和实测值对比Ｆｉｇ.４ＣｏｍｐａｒｉｓｏｎｏｆａｖｅｒａｇｅｆｌｏｗｄｉｓｃｈａｒｇｅꎬｖｅｌｏｃｉｔｙａｎｄｓｅｄｉｍｅｎｔｃｏｎｃｅｎｔｒａｔｉｏｎｂｅｔｗｅｅｎｓｉｍｕｌａｔｅｄａｎｄｍｅａｓｕｒｅｄｒｅｓｕｌｔｓａｔＭｉａｏｈｅｓｔａｔｉｏｎ㊀第３期张帮稳ꎬ等:三峡库区汛期洪峰和沙峰异步运动特性的三维数值模拟４１３㊀之间支流入汇和复杂的局部地形糙率的影响ꎬ经过分析２０１３年７月支流平均流量的总和占万县站流量的１％ꎬ本文中暂时没有考虑支流径流对主河道洪水传播的影响ꎮ平均含沙量的模拟值与庙河站实测时刻的含沙量点较为接近ꎬ和２０１３年水文年鉴中日均含沙量有局部的差别ꎮ这是因为水文年鉴中日均含沙量是基于实测某一时刻的含沙量回归插值得到的[１８]ꎬ本身存在一定的误差ꎬ也可能是限于万县 三峡大坝之间实测资料的限制ꎬ模拟区域局部的冲淤不能很好地反映出来ꎮ图５给出了庙河站断面流速和含沙量分布的瞬时模拟结果ꎮ从图中可以看出流速和含沙量分布与断面的几何形态有关ꎬ断面含沙量的分布不均匀ꎬ呈现出分层的现象ꎮ图６给出了庙河站垂向流速和含沙量的测量值和模拟值的对比结果ꎬ测点的位置见图５ꎬ测点之间间隔１２０ｍꎮ从图中可以看出垂向流速的数值模拟结果和测量值吻合较好ꎬ垂向含沙量的模拟结果和测量值存在一定的误差ꎬ数值模拟断面两侧的测点底部含沙量偏小ꎬ可能与数值模型没有考虑推移质的泥沙运动和断面几何形态有关ꎮ三维水沙数值模型能够进行断面垂向流速和含沙量的分析ꎬ相对于一维和二维数值模拟的精度更高ꎮ图５㊀庙河站断面流速和含沙量模拟结果的瞬时分布Ｆｉｇ.５ＩｎｓｔａｎｔａｎｅｏｕｓｄｉｓｔｒｉｂｕｔｉｏｎｏｆｓｉｍｕｌａｔｅｄｆｌｏｗｖｅｌｏｃｉｔｙａｎｄｓｅｄｉｍｅｎｔｃｏｎｃｅｎｔｒａｔｉｏｎａｔＭｉａｏｈｅｓｔａｔｉｏｎ图６㊀庙河站断面垂向流速和含沙量模拟值和实测值的对比结果Ｆｉｇ.６ＣｏｍｐａｒｉｓｏｎｏｆｖｅｒｔｉｃａｌｆｌｏｗｖｅｌｏｃｉｔｙａｎｄｓｅｄｉｍｅｎｔｃｏｎｃｅｎｔｒａｔｉｏｎｂｅｔｗｅｅｎｓｉｍｕｌａｔｅｄａｎｄｍｅａｓｕｒｅｄｒｅｓｕｌｔｓａｔＭｉａｏｈｅｓｔａｔｉｏｎ３　三峡库区坝前不同蓄水位下洪峰和沙峰异步运动特性三峡工程于２００３年６月进入围堰蓄水期ꎬ坝前水位汛期按１３５ｍ㊁枯季１３９ｍ运行ꎻ２００６年汛后初期蓄水后ꎬ坝前水位按汛期１４４ｍ㊁枯季１５６ｍ运行ꎻ自２００８年汛末三峡水库进行１７５ｍ试验性蓄水以来ꎬ工４１４㊀水科学进展第３２卷㊀程进入１７５ｍ试验性蓄水期ꎮ为了更好地分析三峡水库蓄水以来汛期不同蓄水位下洪峰和沙峰异步运动特性的规律ꎬ本文设置３个研究方案ꎮ方案２为２０１３年７月实测洪水期间坝前水位ꎻ方案１和方案３分别为２０１３年７月实测洪水期间坝前水位减去和加上１０ｍꎮ图７给出了３种研究方案的坝前水位变化ꎮ基于数值模拟的结果分别提取万县 三峡大坝之间重要水文站的流量和含沙量随时间的变化ꎬ用以分析洪峰和沙峰沿程的异步运动特性ꎮ图７㊀不同方案的坝前水位变化Ｆｉｇ.７Ｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｉｎｗａｔｅｒｌｅｖｅｌｉｎｆｒｏｎｔｏｆｔｈｅｄａｍｉｎｄｉｆｆｅｒｅｎｔｓｃｈｅｍｅｓ图８给出了整个研究区域三维流场瞬时的模拟结果ꎬ从图中可以看出上游库区的水流流速较大ꎬ坝前段水流流速较小ꎬ主要由于坝前段水深较大ꎮ图８㊀研究区域三维流场瞬时的模拟结果Ｆｉｇ.８Ｉｎｓｔａｎｔａｎｅｏｕｓｒｅｓｕｌｔｓｏｆｔｈｅｔｈｒｅｅ￣ｄｉｍｅｎｓｉｏｎａｌｓｉｍｕｌａｔｅｄｆｌｏｗｆｉｅｌｄｉｎｔｈｅｓｔｕｄｙｚｏｎｅ图９分别给出了不同方案下重要水文站的流量和含沙量随时间变化的模拟结果对比ꎮ从沿程各水文站流量的模拟结果可以看出ꎬ蓄水位的变化对洪峰传播时间的影响较小ꎻ从含沙量模拟结果可以看出ꎬ蓄水位的变化对沙峰传播时间和大小的影响较大ꎬ随着蓄水位增加ꎬ沙峰滞后洪峰的时间逐渐增加ꎬ其实质是洪水在河道型水库向下游传播过程中ꎬ随着水深的增加ꎬ流速降低导致沙峰传播越来越慢ꎬ沙峰滞后洪峰的时间也越来越大ꎻ同时由于水流流速减小ꎬ水流的挟沙能力降低ꎬ泥沙沿程不断地落淤导致造成沙峰坦化ꎬ沙峰的峰值逐渐减小ꎮ从奉节站㊁巫山站和巴东站的含沙量随时间变化的曲线形态可以看出在巫山站 巴东站库区间存在局部的泥沙冲刷ꎮ为了定量地反映坝前蓄水位的变化对三峡库区洪峰和沙峰异步运动特性的影响ꎬ表２分别给出了不同方案下沿程各水文站的洪峰和沙峰到达时间及沙峰滞后洪峰的时间ꎬ表明坝前蓄水位的变化对洪峰传播时间的影响较小ꎬ对沙峰传播时间的影响较大ꎻ并且沿程各水文站沙峰滞后于洪峰的时间随着坝前蓄水位增加越来越大ꎬ即沙峰滞后洪峰的时间随着水流流速的降低越来越大ꎮ方案２和方案３相对于方案１坝前沙峰滞后于洪峰的时间分别增加了６.４％和１６％ꎮ㊀第３期张帮稳ꎬ等:三峡库区汛期洪峰和沙峰异步运动特性的三维数值模拟４１５㊀图９㊀３种不同方案下沿程各水文站流量和含沙量的模拟结果对比Ｆｉｇ.９Ｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｆｌｏｗｄｉｓｃｈａｒｇｅａｎｄｓｅｄｉｍｅｎｔｃｏｎｃｅｎｔｒａｔｉｏｎａｔｅａｃｈｈｙｄｒｏｇｒａｐｈｉｃｓｔａｔｉｏｎｕｎｄｅｒｄｉｆｆｅｒｅｎｔｓｃｈｅｍｅｓ表２㊀不同方案下沿程各水文站洪峰和沙峰到达时间及滞后时间Ｔａｂｌｅ２Ａｒｒｉｖａｌｔｉｍｅａｎｄｌａｇｔｉｍｅｏｆｆｌｏｏｄｐｅａｋａｎｄｓｅｄｉｍｅｎｔｐｅａｋｒｅａｃｈｉｎｇａｔｅａｃｈｈｙｄｒｏｇｒａｐｈｉｃｓｔａｔｉｏｎｕｎｄｅｒｔｈｅｄｉｆｆｅｒｅｎｔｓｃｈｅｍｅｓ位置方案１方案２方案３洪峰时间沙峰时间滞后时间洪峰时间沙峰时间滞后时间洪峰时间沙峰时间滞后时间奉节１３日０５:００１７日２１:００㊀４ｄ１６ｈ１３日０３:００１８日０６:００５ｄ３ｈ１３日０３:００１８日１６:００５ｄ１３ｈ巫山１３日０６:００１８日１５:００㊀５ｄ９ｈ１３日０５:００１９日０４:００６ｄ０ｈ１３日０４:００１９日１７:００６ｄ１３ｈ巴东１３日０６:００１９日０９:００㊀６ｄ３ｈ１３日０５:００１９日２３:００６ｄ１８ｈ１３日０４:００２０日１３:００７ｄ９ｈ庙河１３日０６:００２０日１９:００㊀７ｄ１３ｈ１３日０５:００２１日０１:００７ｄ２０ｈ１３日０４:００２１日１９:００８ｄ１５ｈ茅坪１３日０６:００２１日０１:００㊀７ｄ１９ｈ１３日０５:００２１日１２:００８ｄ７ｈ１３日０４:００２２日０５:００９ｄ１ｈ４㊀结㊀㊀论本文采用三维水沙数值模型ＳＣＨＩＳＭ对万县 三峡大坝２８０ｋｍ的库区进行了洪水传播和泥沙输移的大尺度数值模拟ꎬ主要对比分析了不同坝前蓄水位下沿程各水文站洪峰和沙峰异步运动特性ꎬ结果表明:(１)三维数值模型ＳＣＨＩＳＭ能够较好地模拟三峡库区长距离的洪水传播和泥沙输移过程ꎬ与实测值验证结果较好ꎮ４１６㊀水科学进展第３２卷㊀(２)三峡水库坝前蓄水位的变化对洪峰传播时间的影响不明显ꎬ对沙峰传播时间的影响较为显著ꎻ坝前水位的增加导致水流流速减小ꎬ沙峰传播减慢ꎬ引起库区主要水文站沙峰滞后于洪峰的时间越来越大ꎮ(３)在万县 三峡大坝库区洪水传播过程中ꎬ随着水深的增加ꎬ水流流速减小ꎬ水流挟沙能力降低ꎬ泥沙不断的落淤导致沙峰峰值减小ꎮ致谢:本研究得到了长江水利委员会水文局许全喜教高㊁武汉大学张为副教授㊁郑珊副教授的帮助和支持ꎬ数值模型在清华大学探索１００集群上计算ꎬ特此致谢!参考文献:[１]陈进.三峡水库建成后长江中下游防洪战略思考[Ｊ].水科学进展ꎬ２０１４ꎬ２５(５):７４５￣７５１.(ＣＨＥＮＪ.ＡｎａｐｐｒｏａｃｈｏｎｆｌｏｏｄｃｏｎｔｒｏｌｓｔｒａｔｅｇｙｉｎｍｉｄｄｌｅａｎｄｌｏｗｅｒｒｅａｃｈｅｓｏｆＹａｎｇｔｚｅＲｉｖｅｒａｆｔｅｒｔｈｅｃｏｍｐｌｅｔｉｏｎｏｆｔｈｅＴｈｒｅｅＧｏｒｇｅｓＤａｍｐｒｏｊｅｃｔ[Ｊ].Ａｄ￣ｖａｎｃｅｓｉｎＷａｔｅｒＳｃｉｅｎｃｅꎬ２０１４ꎬ２５(５):７４５￣７５１.(ｉｎＣｈｉｎｅｓｅ))[２]ＺＨＥＮＧＳＲ.ＲｅｆｌｅｃｔｉｏｎｓｏｎｔｈｅＴｈｒｅｅＧｏｒｇｅｓＰｒｏｊｅｃｔｓｉｎｃｅｉｔｓｏｐｅｒａｔｉｏｎ[Ｊ].Ｅｎｇｉｎｅｅｒｉｎｇꎬ２０１６ꎬ２(４):３８９￣３９７. [３]董炳江ꎬ乔伟ꎬ许全喜.三峡水库汛期沙峰排沙调度研究与初步实践[Ｊ].人民长江ꎬ２０１４ꎬ４５(３):７￣１１.(ＤＯＮＧＢＪꎬＱＩＡＯＷꎬＸＵＱＸ.ＳｔｕｄｙｏｎｓｅｄｉｍｅｎｔｐｅａｋｒｅｇｕｌａｔｉｏｎｏｆＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒｉｎｆｌｏｏｄｓｅａｓｏｎａｎｄｐｒｅｌｉｍｉｎａｒｙｐｒａｃｔｉｃｅ[Ｊ].Ｙａｎ￣ｇｔｚｅＲｉｖｅｒꎬ２０１４ꎬ４５(３):７￣１１.(ｉｎＣｈｉｎｅｓｅ))[４]董炳江ꎬ陈显维ꎬ许全喜.三峡水库沙峰调度试验研究与思考[Ｊ].人民长江ꎬ２０１４ꎬ４５(１９):１￣５.(ＤＯＮＧＢＪꎬＣＨＥＮＸＷꎬＸＵＱＸ.ＩｎｖｅｓｔｉｇａｔｉｏｎｓａｎｄｃｏｎｓｉｄｅｒａｔｉｏｎｓｏｎｐｅａｋｓｅｄｉｍｅｎｔｒｅｇｕｌａｔｉｏｎｏｆＴｈｒｅｅＧｅｏｒｇｅｓＲｅｓｅｒｖｏｉｒ[Ｊ].ＹａｎｇｔｚｅＲｉｖｅｒꎬ２０１４ꎬ４５(１９):１￣５.(ｉｎＣｈｉｎｅｓｅ))[５]黄仁勇.长江上游梯级水库泥沙输移与泥沙调度研究[Ｄ].武汉:武汉大学ꎬ２０１６.(ＨＵＡＮＧＲＹ.ＳｔｕｄｙｏｎｔｈｅｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔａｎｄｓｅｄｉｍｅｎｔｏｐｅｒａｔｉｏｎｏｆｃａｓｃａｄｅｒｅｓｅｒｖｏｉｒｓｉｎｔｈｅＵｐｐｅｒＣｈａｎｇｊｉａｎｇＲｉｖｅｒ[Ｄ].Ｗｕｈａｎ:ＷｕｈａｎＵｎｉｖｅｒｓｉｔｙꎬ２０１６.(ｉｎＣｈｉｎｅｓｅ))[６]张地继ꎬ董炳江ꎬ杨霞ꎬ等.三峡水库库区沙峰输移特性研究[Ｊ].人民长江ꎬ２０１８ꎬ４９(２):２３￣２８ꎬ６８.(ＺＨＡＮＧＤＪꎬＤＯＮＧＢＪꎬＹＡＮＧＸꎬｅｔａｌ.ＲｅｓｅａｒｃｈｏｎｓｅｄｉｍｅｎｔｐｅａｋｍｏｖｅｍｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｉｎＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ[Ｊ].ＹａｎｇｔｚｅＲｉｖｅｒꎬ２０１８ꎬ４９(２):２３￣２８ꎬ６８.(ｉｎＣｈｉｎｅｓｅ))[７]ＲＥＮＪＱꎬＺＨＡＯＭＤꎬＺＨＡＮＧＷꎬｅｔａｌ.ＩｍｐａｃｔｏｆｔｈｅｃｏｎｓｔｒｕｃｔｉｏｎｏｆｃａｓｃａｄｅｒｅｓｅｒｖｏｉｒｓｏｎｓｕｓｐｅｎｄｅｄｓｅｄｉｍｅｎｔｐｅａｋｔｒａｎｓｐｏｒｔｖａｒｉａｔｉｏｎｄｕｒｉｎｇｆｌｏｏｄｅｖｅｎｔｓｉｎｔｈｅＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ[Ｊ].Ｃａｔｅｎａꎬ２０２０ꎬ１８８:１０４４０９.[８]张小峰ꎬ徐顺ꎬ卢新华ꎬ等.一维悬移质泥沙输移方程及其滞后和离散效应[Ｊ].泥沙研究ꎬ２０１８(４):８￣１５.(ＺＨＡＮＧＸＦꎬＸＵＳꎬＬＵＸＨꎬｅｔａｌ.Ｈｙｓｔｅｒｅｓｉｓａｎｄｄｉｓｐｅｒｓｉｏｎｉｎｏｎｅｄｉｍｅｎｓｉｏｎａｌｓｕｓｐｅｎｄｅｄｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔｅｑｕａｔｉｏｎ[Ｊ].ＪｏｕｒｎａｌｏｆＳｅｄ￣ｉｍｅｎｔＲｅｓｅａｒｃｈꎬ２０１８(４):８￣１５.(ｉｎＣｈｉｎｅｓｅ))[９]夏军强ꎬ张晓雷ꎬ邓珊珊ꎬ等.黄河下游高含沙洪水过程一维水沙耦合数学模型[Ｊ].水科学进展ꎬ２０１５ꎬ２６(５):６８６￣６９７.(ＸＩＡＪＱꎬＺＨＡＮＧＸＬꎬＤＥＮＧＳＳꎬｅｔａｌ.ＭｏｄｅｌｌｉｎｇｏｆｈｙｐｅｒｃｏｎｃｅｎｔｒａｔｅｄｆｌｏｏｄｓｉｎｔｈｅＬｏｗｅｒＹｅｌｌｏｗＲｉｖｅｒｕｓｉｎｇａｃｏｕｐｌｅｄａｐｐｒｏａｃｈ[Ｊ].ＡｄｖａｎｃｅｓｉｎＷａｔｅｒＳｃｉｅｎｃｅꎬ２０１５ꎬ２６(５):６８６￣６９７.(ｉｎＣｈｉｎｅｓｅ))[１０]ＬＵＹＪꎬＷＡＮＧＺＹ.３￣ＤｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｆｏｒｗａｔｅｒｆｌｏｗｓａｎｄｓｅｄｉｍｅｎｔｄｅｐｏｓｉｔｉｏｎｉｎｄａｍａｒｅａｓｏｆｔｈｅＴｈｒｅｅＧｏｒｇｅｓＰｒｏｊｅｃｔ[Ｊ].ＪｏｕｒｎａｌｏｆＨｙｄｒａｕｌｉｃＥｎｇｉｎｅｅｒｉｎｇꎬ２００９ꎬ１３５(９):７５５￣７６９.[１１]ＪＩＡＤＤꎬＳＨＡＯＸＪꎬＺＨＡＮＧＸＮꎬｅｔａｌ.Ｓｅｄｉｍｅｎｔａｔｉｏｎｐａｔｔｅｒｎｓｏｆｆｉｎｅ￣ｇｒａｉｎｅｄｐａｒｔｉｃｌｅｓｉｎｔｈｅｄａｍａｒｅａｏｆｔｈｅＴｈｒｅｅＧｏｒｇｅｓＰｒｏｊｅｃｔ:３￣Ｄｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ[Ｊ].ＪｏｕｒｎａｌｏｆＨｙｄｒａｕｌｉｃＥｎｇｉｎｅｅｒｉｎｇꎬ２０１３ꎬ１３９(６):６６９￣６７４.[１２]ＺＨＡＮＧＹＬꎬＢＡＰＴＩＳＴＡＡＭ.ＳＥＬＦＥ:ａｓｅｍｉ￣ｉｍｐｌｉｃｉｔＥｕｌｅｒｉａｎ￣Ｌａｇｒａｎｇｉａｎｆｉｎｉｔｅ￣ｅｌｅｍｅｎｔｍｏｄｅｌｆｏｒｃｒｏｓｓ￣ｓｃａｌｅｏｃｅａｎｃｉｒｃｕｌａ￣ｔｉｏｎ[Ｊ].ＯｃｅａｎＭｏｄｅｌｌｉｎｇꎬ２００８ꎬ２１(３/４):７１￣９６.[１３]ＺＨＡＮＧＹＪꎬＹＥＦꎬＳＴＡＮＥＶＥＶꎬｅｔａｌ.Ｓｅａｍｌｅｓｓｃｒｏｓｓ￣ｓｃａｌｅｍｏｄｅｌｉｎｇｗｉｔｈＳＣＨＩＳＭ[Ｊ].ＯｃｅａｎＭｏｄｅｌｌｉｎｇꎬ２０１６ꎬ１０２:６４￣８１.[１４]ＵＭＬＡＵＦＬꎬＢＵＲＣＨＡＲＤＨ.Ａｇｅｎｅｒｉｃｌｅｎｇｔｈ￣ｓｃａｌｅｅｑｕａｔｉｏｎｆｏｒｇｅｏｐｈｙｓｉｃａｌｔｕｒｂｕｌｅｎｃｅｍｏｄｅｌｓ[Ｊ].ＪｏｕｒｎａｌｏｆＭａｒｉｎｅＲｅｓｅａｒｃｈꎬ２００３ꎬ６１(２):２３５￣２６５.[１５]涂启华ꎬ杨赉斐.泥沙设计手册[Ｍ].北京:中国水利水电出版社ꎬ２００６.(ＴＵＱＨꎬＹＡＮＧＬＦ.Ｓｅｄｉｍｅｎｔｄｅｓｉｇｎｍａｎｕａｌ[Ｍ].Ｂｅｉｊｉｎｇ:ＣｈｉｎａＷａｔｅｒＰｏｗｅｒＰｒｅｓｓꎬ２００６.(ｉｎＣｈｉｎｅｓｅ))㊀第３期张帮稳ꎬ等:三峡库区汛期洪峰和沙峰异步运动特性的三维数值模拟４１７㊀[１６]张地继ꎬ董炳江ꎬ杨霞ꎬ等.三峡水库库区沙峰输移特性研究[Ｊ].人民长江ꎬ２０１８ꎬ４９(２):２３￣２８.(ＺＨＡＮＧＤＪꎬＤＯＮＧＢＪꎬＹＡＮＧＸꎬｅｔａｌ.ＲｅｓｅａｒｃｈｏｎｓｅｄｉｍｅｎｔｐｅａｋｍｏｖｅｍｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｉｎＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ[Ｊ].ＹａｎｇｔｚｅＲｉｖｅｒꎬ２０１８ꎬ４９(２):２３￣２８.(ｉｎＣｈｉｎｅｓｅ))[１７]ＺＨＡＮＧＹＪꎬＡＴＥＬＪＥＶＩＣＨＥꎬＹＵＨＣꎬｅｔａｌ.Ａｎｅｗｖｅｒｔｉｃａｌｃｏｏｒｄｉｎａｔｅｓｙｓｔｅｍｆｏｒａ３￣Ｄｕｎｓｔｒｕｃｔｕｒｅｄ￣ｇｒｉｄｍｏｄｅｌ[Ｊ].ＯｃｅａｎＭｏｄｅｌｌｉｎｇꎬ２０１５ꎬ８５:１６￣３１.[１８]周波ꎬ许全喜ꎬ李雨.三峡水库入库泥沙实时监测试验研究[Ｊ].水文ꎬ２０１６ꎬ３６(４):５３￣５７.(ＺＨＯＵＢꎬＸＵＱＸꎬＬＩＹ.Ｏｎｎｅｗｔｅｃｈｎｏｌｏｇｙｏｆｓｅｄｉｍｅｎｔｒｅａｌ￣ｔｉｍｅｍｏｎｉｔｏｒｉｎｇｉｎＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ[Ｊ].ＪｏｕｒｎａｌｏｆＣｈｉｎａＨｙｄｒｏｌｏｇｙꎬ２０１６ꎬ３６(４):５３￣５７.(ｉｎＣｈｉｎｅｓｅ))Ｔｈｒｅｅ￣ｄｉｍｅｎｓｉｏｎａｌｎｕｍｅｒｉｃａｌｏｆｔｈｅａｓｙｎｃｈｒｏｎｏｕｓｍｏｖｅｍｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｆｌｏｏｄａｎｄｓｅｄｉｍｅｎｔｐｅａｋｓｉｎｔｈｅＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ∗ＺＨＡＮＧＢａｎｇｗｅｎꎬＷＵＢａｏｓｈｅｎｇꎬＺＨＡＮＧＲｕｏｙｉｎ(ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙｏｆＨｙｄｒｏｓｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇꎬＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙꎬＢｅｉｊｉｎｇ１０００８４ꎬＣｈｉｎａ)Ａｂｓｔｒａｃｔ:ＡｆｔｅｒｔｈｅｉｍｐｏｕｎｄｍｅｎｔｏｆｔｈｅＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒꎬｔｈｅｉｎｃｒｅａｓｅｏｆｗａｔｅｒｄｅｐｔｈｃａｕｓｅｄｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｌａｇｔｉｍｅｏｆｔｈｅｓｅｄｉｍｅｎｔｐｅａｋｂｅｈｉｎｄｔｈｅｆｌｏｏｄｐｅａｋｉｎｔｈｅｆｌｏｏｄｐｒｏｐａｇａｔｉｏｎ.Ｔｈｅａｓｙｎｃｈｒｏｎｏｕｓｍｏｖｅｍｅｎｔｃｈａｒａｃ￣ｔｅｒｉｓｔｉｃｓｂｅｔｗｅｅｎｔｈｅｆｌｏｏｄａｎｄｓｅｄｉｍｅｎｔｐｅａｋｓｃａｎｂｅｕｔｉｌｉｚｅｄｔｏｃｏｎｄｕｃｔｔｈｅｓｅｄｉｍｅｎｔｐｅａｋｒｅｇｕｌａｔｉｏｎꎬｏｎｅｏｆｔｈｅｃｒｕ￣ｃｉａｌｍｅａｓｕｒｅｓｔｏｒｅｄｕｃｅｔｈｅｒｅｓｅｒｖｏｉｒｓｅｄｉｍｅｎｔａｔｉｏｎ.Ｔｈｅｔｈｒｅｅ￣ｄｉｍｅｎｓｉｏｎａｌｎｕｍｅｒｉｃａｌｍｏｄｅｌｎａｍｅｄＳＣＨＩＳＨ(Ｓｅｍｉ￣ｉｍｐｌｉｃｉｔＣｒｏｓｓ￣ｓｃａｌｅＨｙｄｒｏｓｃｉｅｎｃｅＩｎｔｅｇｒａｔｅｄＳｙｓｔｅｍＭｏｄｅｌ)ｗａｓａｄｏｐｔｅｄｔｏｓｔｕｄｙｔｈｅａｓｙｎｃｈｒｏｎｏｕｓｍｏｖｅｍｅｎｔｃｈａｒａｃ￣ｔｅｒｉｓｔｉｃｓｏｆｔｈｅｆｌｏｏｄａｎｄｓｅｄｉｍｅｎｔｐｅａｋｓｉｎｔｈｅＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ.Ｔｈｅｒｅｌｉａｂｉｌｉｔｙａｎｄａｃｃｕｒａｃｙｏｆｔｈｅｍｏｄｅｌｗｅｒｅｖｅｒｉｆｉｅｄｂｙｃｏｍｐａｒｉｎｇｔｈｅｗａｔｅｒａｎｄｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔｄｕｒｉｎｇｔｈｅｆｌｏｏｄｓｅａｓｏｎｉｎ２０１３ꎬｔｈｅｎｔｈｅａｓｙｎｃｈｒｏｎｏｕｓｍｏｖｅ￣ｍｅｎｔｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｔｈｅｆｌｏｏｄａｎｄｓｅｄｉｍｅｎｔｐｅａｋｓｕｎｄｅｒｄｉｆｆｅｒｅｎｔｉｍｐｏｕｎｄｅｄｌｅｖｅｌｓｗｅｒｅｐｒｅｌｉｍｉｎａｒｉｌｙａｎａｌｙｚｅｄｂａｓｅｄｏｎｔｈｅｍｏｄｅｌｒｅｓｕｌｔｓ.Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｅｎｕｍｅｒｉｃａｌｍｏｄｅｌｃｏｕｌｄａｃｃｕｒａｔｅｌｙｓｉｍｕｌａｔｅｔｈｅｌｏｎｇ￣ｄｉｓｔａｎｃｅｐｒｏｐａｇａｔｉｏｎｏｆｔｈｅｆｌｏｏｄａｎｄｓｅｄｉｍｅｎｔｐｅａｋｓ.Ｔｈｅｉｍｐｏｕｎｄｅｄｌｅｖｅｌｉｎｆｒｏｎｔｏｆｔｈｅｄａｍｈａｄｌｉｔｔｌｅｉｎｆｌｕｅｎｃｅｏｎｔｈｅｐｒｏｐ￣ａｇａｔｉｏｎｔｉｍｅｏｆｔｈｅｆｌｏｏｄｐｅａｋꎬｂｕｔｈａｄａｓｉｇｎｉｆｉｃａｎｔｅｆｆｅｃｔｏｎｔｈｅｐｒｏｐａｇａｔｉｏｎｔｉｍｅｏｆｔｈｅｓｅｄｉｍｅｎｔｐｅａｋ.Ｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｉｍｐｏｕｎｄｅｄｌｅｖｅｌｒｅｓｕｌｔｅｄｉｎｔｈｅｄｅｃｒｅａｓｅｏｆｔｈｅｆｌｏｗｖｅｌｏｃｉｔｙａｎｄｔｈｅｓａｎｄ￣ｃａｒｒｙｉｎｇｃａｐａｃｉｔｙꎬｓｏｔｈｅｐｒｏｐａｇａｔｉｏｎｓｐｅｅｄｏｆｔｈｅｓｅｄｉｍｅｎｔｐｅａｋｗａｓｓｌｏｗｅｒａｎｄｔｈｅｓｅｄｉｍｅｎｔｐｅａｋｖａｌｕｅｗａｓｓｍａｌｌｅｒꎬｗｈｉｃｈｌｅｄｔｏｔｈｅｌｏｎｇｅｒｌａｇｔｉｍｅｏｆｔｈｅｓｅｄｉｍｅｎｔｐｅａｋｂｅｈｉｎｄｔｈｅｆｌｏｏｄｐｅａｋａｌｏｎｇｔｈｅｒｅｓｅｒｖｏｉｒ.ＴｈｅｒｅｓｅａｒｃｈｒｅｓｕｌｔｓｃａｎｐｒｏｖｉｄｅｓｃｉｅｎｔｉｆｉｃｅｖｉｄｅｎｃｅｆｏｒｆｕｒｔｈｅｒｏｐｔｉｍｉｚａｔｉｏｎｏｆｔｈｅｏｐｅｒａｔｉｏｎｏｆｔｈｅｓｅｄｉｍｅｎｔｐｅａｋｒｅｇｕｌａｔｉｏｎｉｎｔｈｅＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ.Ｋｅｙｗｏｒｄｓ:ｓｅｄｉｍｅｎｔｐｅａｋꎻｆｌｏｏｄｐｅａｋꎻｌａｇｔｉｍｅꎻｓｅｄｉｍｅｎｔｃｏｎｃｅｎｔｒａｔｉｏｎꎻＴｈｒｅｅＧｏｒｇｅｓＲｅｓｅｒｖｏｉｒ∗ＴｈｅｓｔｕｄｙｉｓｆｉｎａｎｃｉａｌｌｙｓｕｐｐｏｒｔｅｄｂｙｔｈｅＮａｔｉｏｎａｌＫｅｙＲ＆ＤＰｒｏｇｒａｍｏｆＣｈｉｎａ(Ｎｏ.２０１７ＹＦＣ０４０５２０２)ａｎｄｔｈｅＮａｔｉｏｎａｌＮａｔｕｒａｌＳｃｉｅｎｃｅＦｏｕｎｄａｔｉｏｎｏｆＣｈｉｎａ(Ｎｏ.５１６３９００５).。

水资源配置中的生态环境用水问题浅析国家电力公司华东勘测设计研究院傅菁菁张祖生摘要：根据目前国内外生态环境用水的研究情况，对解决水资源配置中的生态环境用水问题提出了几点思路，并综述了各项生态环境用水量的计算方法。

关键词：生态环境用水计算方法生态环境是人类生存发展的基本自然条件，水对维持人类赖以生存的生态环境起着决定性作用，而长期以来，人类在水资源的开发利用中，只重视解决生产、生活用水的需求，对水资源保护利用与生态环境的相互依存关系缺乏认识，忽视生态环境用水，从而导致了一系列的生态环境问题。

海河流域曾是水域、湿地广泛分布的地区，生态环境状况良好，而目前其水资源开发利用率已达80％以上，严重超出世界公认的30％开发利用率上限，已造成流域内的许多河流断流、湖泊湿地干涸以及地面沉降。

新疆塔里木河及甘肃、内蒙古黑河下游，河水断流、湖泊干涸、湿地消失、林草植被大面积死亡；黄河下游河段长时间枯水、断流、黄河三角洲湿地退化、地下水位下降等，无一不是由于生态环境水量减少引起的。

美国中央河谷工程和加利福尼亚州水道工程都是以萨克拉门托河为取水水源，调水量较大，造成海湾水质恶化，影响海湾水生生物，海水倒灌使旧金山湾地区的土地盐渍化，给生态环境带来了较大不利影响[1]。

科罗拉多河的地表水，绝大部分用于灌溉美国7个州的含盐土地，每年引水量约为95亿m3，致使河水含盐量不断增加，被浇灌的土地盐碱化[1]。

因此，近年来，有许多研究者强力疾呼我国所面临的水资源问题，向全社会提出，“在水资源配置中，要从不重视生态环境用水，转变为在保证生态环境用水的前提下，合理规划和保障经济社会的用水”[16]。

“应该在水资源总量中专门划出一部分作为生态环境用水，另一部分为国民经济各部门的用水，包括工、农业及城市生活用水等”[17]。

1.生态环境用水的基本概念目前，对生态环境用水的概念尚无统一的定义。

生态用水量应该是特定区域内生态系统用水量的总称，包括生物体自身的用水量和生物体赖以生存的环境用水量，其实质就是维持生态系统生物群落及其栖息环境动态稳定所需的用水量。

人类活动影响下营养盐向河口/近海的输出和模型研究= Summary of human activities on global nutrient export from watersheds to estuaries and coastal water:biogeochemical cycles and modeling/晏维金地理研究.2006,25(5).825~ 835人类活动(人口增长和经济发展所带来的对粮食和能源的需求和消费等)已经大大改变了地表圈层(土壤圈-水圈-大气圈)营养盐(主要是氮、磷、硅)的生物地球化学循环,结果使许多陆地、淡水和海洋生态系统收到显著影响并导致了严重的环境质量退化.本文对2002-2005年由联合国教科文组织(UNESCO)和美国国家自然科学基金委员会资助的全球营养盐输出和模型!(Global NEWS)项目进行了总结.图3表2参52Be P7362007010105海南岛滨珊瑚18O对环境条件的响应=Response to envi ronmental conditions of coral18O of proites lutea from Shalao,Hainan Island/宋少华,周卫健海洋地质与第四纪地质.2006,26(4).23~28通过与同期观测记录的海洋表面温度、海洋表层盐度进行比较分析,提出珊瑚18O对环境条件变化的响应特征.该地区珊瑚18O的年际变化主要受海洋表面温度的控制,而以夏季为代表的季节性变化则主要受海洋表层盐度的影响图6参29(许桂红)N P731.232007010106辐射沙洲东大港潮流水道悬沙输移机制分析=Mechanism of suspended sediment transport in Dongdagang tidal channel of radial sand r idges/吴德安,张忍顺河海大学学报. 2006,34(2).216~222以东大港潮流水道4号站位的水、沙实测资料为例,利用悬沙通量机制分解法、功率谱分析法,对东大港潮流水道悬沙输移的动力机制进行分析探讨.图4表6参12(许桂红)N#湖####沼####学#湖##泊##学P343.32007010107 2002年前后博斯腾湖水位变化及其对中亚气候变化的响应=Water level change in Bo sten Lake under the climatic var ia tion background of central Asian around2002/王润,孙占东冰川冻土.2006,28(3).324~329通过对湖泊补给水量的分析,博斯腾湖2002年前后水位由高向低的转折变化是由开都河径流变化造成的.对2002年前后天山南北其它主要湖泊的水位和河流径流变化作了比较分析结论,由于山区河源径流补给组成不同,表现出在同一中亚气候背景下,即2002年前后该地区气温降低而降水东西各有差异,使得天山西段受降水补给为主的河流,2002年后径流仍有增加.图5表1参11(洪明)Bi X1412007010108乌梁素海沉积物中全磷的分布特征=Distr ibution of total phosphorus in sediments of Wuliangsuhai Lake/孙惠民,何江沉积学报.2006,24(4).579~584开展了浅水草型湖泊乌梁素海沉积物中TP的分布特征研究.结果表明,乌梁素海表层沉积物中TP含量在0.37 -0.99g/kg;水平分布上,TP含量呈现从湖区四周向湖心递减的变化趋势,大部分湖区已具有生态毒性效应,并对湖区底栖生物群落和生态环境构成较大威胁.图4表2参15 (许桂红)N P343.3,P512.322007010109内蒙古黄旗海不同粒级湖泊沉积物R、S组成与环境变化=f f f RS ff z f f Huangqihai Lake,Inner Mongolia/申洪源,贾玉连地理学报.2006,61(11).1208~1217内蒙古黄旗海湖泊沉积物剖面层次结构清晰完整,通过对该剖面全样、77m-20m和小于20m三个不同粒级Rb、Sr含量和磁化率的系统测定,结果表明,湖泊沉积物中小于20m粒级成分的Rb、Sr含量和Rb/Sr值变化阶段明显,Rb/Sr值与磁化率分布曲线十分吻合,因此,小于20m 粒级成分的Rb/Sr值与磁化率包含更加真实的古气候环境信息.R b/Sr值反映了风化的强度,因而可以作为指示夏季风环流强度的代用指标.黄旗海湖泊沉积物的地球化学Rb、Sr元素含量和Rb/Sr值以及磁化率值波动特征具有深刻的环境演变背景,由此揭示的环境特征在一定程度上反映了全新世气候千百年尺度的波动变化,与全球气候波动有关.图6表1参32(广兵)Be P534.632007010110青藏高原中部湖泊岩芯记录的第四纪湖平面变化及气候意义=Lake level change r ecorded by core of the Quaternar y la custrine sediment in t he centr al T ibetan Plateau and its climatic implicat ions/伊海生,张小青地学前缘.2006,13 (5).300~307通过对羌塘地区中部羌D1井长1816m第四系岩芯的研究,分析湖泊沉积记录的岩相旋回和Fe/Mn和Sr/Ba等微量元素古环境参数的变化,结合样品的热释光(TL)测年数据,讨论晚更新世该区湖泊扩张和湖面升降变化过程.结果表明,晚更新世以来该区湖泊环境的演化,经历了二次湖进过程,二次湖平面上升期之间出现了一次广泛的干化低水位时期其中5万年期间的湖泛事件可以在区域上进行追踪对比,它对应于深海氧同位素的第5阶段,可以作为高原中央气候转型期的标志图表参B5b rEnvironmental change in erred rom distribution o b and r in di erent grain si e ractions rom lacustrine sediments in .1-7.2220e1X 522007010111吉林松花湖叶绿素的垂直分布及其影响因素=T he vertical dynamic changes in chlorophylla and its relative factors I n Songhua Lake in Jilin/金香琴,盛连喜干旱环境监测.2006,20(1).23~26,封底图2表1参9(王彦)SDP641.62007010112岱海水盐变化原因及影响研究=Study on the causes and im pacts of water salt change of the Daihai Lake/孙占东,姜加虎干旱区研究.2006,23(2).264~268从湖泊水文要素变化出发,基于湖泊和流域的关系,利用相关分析方法,对近50a 来岱海水位、矿化度等水文要素变化序列进行了分析.并对变化原因进行了探讨,结果表明:相对封闭及较小尺度的流域面积等自身特性,决定了其水资源系统自我调节能力的脆弱性;而近年来流域上的人类活动,增大了降水径流滞留时间和蒸发作用:蒸发的水汽随环流流出流域,最终导致汇人流域的持水量和汇入湖泊的水量不断减少,矿化度在湖水量不断减少的情况下浓缩升高.最后就变化趋势及影响进行了讨论.图6表1参13(作者)SDP512.322007010113青藏高原东北部希门错湖岩心粒度特征及其环境意义=T he grain size characteristics of Ximencuo Lake core in the northeast Tibet Plateau and its environmental significance/类延斌,张成君海洋地质与第四纪地质.2006,26(3).31~38位于青藏高原东北部的希门错为一冰川深水湖,流域面积小,受人类活动影响较小,忠实地记录了湖泊的自然演化过程.结合希门错湖区特点,对沉积物的粒度特征进行了详细分析.图5参22(许桂红)NP597.22007010114长江中下游湖泊水体氧同位素组成=Character istics of oxy gen isotopic composition of the lakes along the mid lower reach es of the Yangtze River /吴敬禄,林森海洋地质与第四纪地质.2006,26(3).53~56图1表2参9(许桂红)NP426.22007010115鄱阳湖水面蒸发量的计算与变化趋势分析(1955-2004年)=Calculat ion of lake evaporation and tr end analysis of Lake Poyang:19552004/闵骞,刘影湖泊科学.2006,17(5).452~457利用器测折算法与气候模式法,分别计算鄱阳湖周围康山、棠荫、都昌、星子、湖口5站的单站水面蒸发量,以5站两种方法计算值的平均值代表鄱阳湖大湖面的水面蒸发量,求得鄱阳湖1955-2004年各月的水面蒸发量和蒸发水量.图1表6参8(许桂红)NP34.32007010116鄱阳湖湖口、外洲、梅港三站水沙变化及趋势分析(1955-2001年)=Analysi s on the runoff and sediment transpor tation ,W z M f L y 55郭鹏,陈晓玲湖泊科学6,(5)5~63对鄱阳湖流域三个主要控制站湖口、外洲和梅港多年(1955-2001年)水沙数据进行了统计分析,利用滑动平均法,Spearman 秩次相关检验、线性回归检验方法对三个测站的水沙变化趋势进行了分析检验.从赣江和信江水沙总体变化趋势来看,赣江径流量变化趋势不明显,而输沙量具有明显减少的趋势;信江径流量增加趋势明显,输沙量基本无明显趋势.图2表3参8(许桂红)NP343.32007010117洪湖分蓄洪区洪水淹没风险动态识别与可能损失评估=Dynamic inundation risk identification and est imation of the po tential loss in Honghu flood diversion area,China/Marco GEM MER ,王国杰湖泊科学.2006,17(5).464~469以洪湖分蓄洪区为案例,采用基于GIS 栅格数据整合于Ar cview3.x 的二维水文-水动力学模型进行洪水淹没风险动态识别,并且根据土地利用分类及其单位面积价值,建立洪水淹没损失函数,进行洪水淹没动态损失评估,建立了东洪湖分蓄洪区洪水淹没动态损失数据库,为东洪湖分蓄洪区的合理利用提供定量科学依据.图3表2参18(许桂红)NP463.252007010118气候背景下冰川在博斯腾湖水量平衡中的作用=Effect of glaciers change to water balance of Lake Bosten under climatic backgrounds/孙占东,王润湖泊科学.2006,17(5).484~489基于博斯腾湖水量平衡关系,对博斯腾湖水量变化做了分析,认为博斯腾湖近年水位的显著变化与占其入湖水量85%的开都河流量变化有直接关系.冰川作为开都河上游重要的水资源形势,造就了开都河稳定的基流.图5表1参16(许桂红)NP463.242007010119西太湖湖滨典型河网区与太湖水量的交换=Water quantity exchanges between typical river network area and western Lake Taihu/庄巍,逄勇湖泊科学.2006,17(5).490~494以西太湖湖滨典型河网区宜兴大浦河网区为研究区域,对于该典型河网区与太湖之间的水量交换情况进行了系统的计算研究.建立起适合该河网区域的非稳态河网水量模型,结合径流计算结果,对各年内典型河网区内水流的动态变化进行了计算模拟,对典型河网区与太湖水体的水量交换情况进行系统研究,分别对顺流与逆流情况下交换水量的年内分布与年际变化规律进行了分析.图6参13(许桂红)NP461.52007010120风浪扰动对太湖水体悬浮物重金属含量的影响=Effects of wind wave disturbance on heavy metal contents in suspended solids of lake Taihu/池俏俏,朱广伟湖泊科学.2006,17(5).495~498通过离心浓缩的方法,获取太湖梅梁湾口东岸处(即梅梁湾与贡湖湾的交界处)不同风浪条件下的悬浮颗粒物,冷冻干燥,微波消解,I S 的方法测定了其中、、、N 、、Z 等重金属元素及、、F 、M 等相关金属元6in the Hukou ai hou and eigang stations o ake Po ang dur ing 192001/.20017.484CP -AE Co Cr Cu i Pb n Al Ca e n 1素的含量.图4参14(许桂红)N P461.52007010121广州麓湖水-气界面多环芳烃的通量分析=Deposition and exchange of polycyclic aromatic hydrocarbons across air water interface in an ur ban lake:Lake Luhu,Guangzhou/陈灿云,张干湖泊科学.2006,17(5).499~502图2表1参11(许桂红)N X522007010122微生物作用对太湖沉积物磷释放影响的模拟实验研究= M icrobial effects on phosphorus release in Taihu Lake%s sedi ments/孙晓杭,张昱环境化学.2006,25(1).24~ 27图3参9(许桂红)N P512.22007010123青藏高原湖泊沉积研究及其进展=Researches on lacustr ine sediments in Tibetan Plateau and their pr ogress/于守兵,李世杰山地学报.2005,23(6).663~671青藏高原湖泊沉积研究主要围绕青藏高原隆升和全球变化响应开展.研究揭示了青藏高原隆升具有整体性、阶段性和后期加速性.中晚更新世以来和末次间冰期-冰期气候演化过程得到重建,并能够与冰芯和深海氧同位素记录对比;同时也存在MI S3阶段强烈暖湿和末次冰盛期冷湿等区域特征;新仙女木事件在湖泊沉积物也有明显记录.全新世研究表明青藏高原早期暖湿并经历冷事件,大暖期普遍出现高水位,后期气候向干冷化方向发展.湖泊沉积环境定量化重建也得到研究.青藏高原湖泊沉积应在高分辨率纹层沉积和环境指标定量分析基础理论方面加强研究.参70 (作者)Che K992007010124太湖流域主要湖泊的水域动态变化=Study on surface ar ea change of major lakes in Taihu basin based on topographical and r emotely sensed information/李新国,江南水资源保护.2006,22(3).20~23利用GIS和遥感技术对不同时期的地形图和遥感影像数据进行处理,获取3个时期的太湖流域主要湖泊的水域面积.结合多种文献资料,对太湖流域主要湖泊的水域面积变化进行了研究,并分析了引起水域面积变化的主要因素.结果表明,1971-2002年期间太湖流域主要湖泊的水域面积减少了188.87km2;湖泊水域面积减少的主要原因是由于人类活动的影响.在实地调查研究的基础上,对合理利用与保护湖泊水资源的对策进行了探讨.图1表1参10(盛春蕾)Cha#沼##泽##学P343.42007010125中国海岸湿地环境特点与开发利用=Character istics and ex ploitation of wetland of China/王颖,朱大奎长江流域资源与环境.2006,15(5).553~559中国海岸湿地环境受季风波浪、潮汐与大河影响的河海交互作用为特征,地跨3个纬度带,其面积约占中国湿地总面积的5按成因与表相可分为河口芦苇、草滩湿地,平原海岸草滩盐沼湿地及在华南隐蔽港湾顶部红树林沼泽3种类型.均具有潮上带、潮间带与潮下带分带特征.指出:海岸湿地是珍贵的新生空间资源,但由于滥伐垦殖、外来物种属入侵危害、陆源水沙减少及海平面上升招致海岸侵蚀、盐湖入侵与内涝频繁等灾害影响,海岸湿地环境质量下降,面临着解决海岸湿地生态保护与沿海人民生计这一新生矛盾问题.图6表2参11(宋金叶)HN P343.42007010126长江下游典型水域湿地近百年变化分析=Changes in typical water wetland in the lower reaches of the Yangtze R iver in t he past100years/赵广举,张鹏长江流域资源与环境. 2006,15(5).560~562长江中下游地区分布着以湖泊和河流湿地为主的大面积湿地,同时也是人口密集区,水域湿地在近百年发生了巨大变化.选择长江下游太湖平原上的典型湿地为例,根据历史史料记载和近代精确的地形图遥感资料,首次得出了研究区近百年水域湿地的精确面积和空间分布状况,分析结果显示:1916年建昌圩水域面积为11.30km2,沼泽面积较小;1957年水域面积减少到10.32km2,沼泽面积增加到0. 86;1966年水域湿地面积较少到9.18km2,到1981年水域湿地面积减至7km2,到2000年水域湿地面积减至5. 50km2.并对影响其变化的因素进行了分析.图2表2参9 (宋金叶)HN P343.42007010127长江中下游湿地保护与流域生态管理=Wetland protection and basin management in the middle and lower reaches of t he Yangtze River/王学雷,许厚泽长江流域资源与环境.2006,15(5).564~567长江中下游湿地保护必须从流域管理角度进行规划和保护,由过去主要是单块湿地保护向按流域系统保护转变.具体可采取以下措施:在流域内建立统一协调机制,对流域湿地进行保护与合理利用,合理布局,统一规划;按湖泊流域和物种分布整合现有保护区,建立新的湿地保护区,解决目前管理上的制约问题;在大力保护的前提下科学合理利用长江中下游湿地资源,开创新的生产力;同时加强湿地科学研究.表1参6(宋金叶)HN P941.272007010128江汉平原湿地水系统综合评价与水资源管理探讨=Com prehensive evaluat ion on wetland water system and water r e sour ces management in the Jianghan plain/王学雷,吕宪国地理科学.2006,26(3).311~315在探讨江汉平原湿地水系统结构和功能的基础上,提出了合理利用湿地水资源和实施科学水管理的途径,首先要对水资源和水环境进行积极的保护;其次积极实施退田还湖,增强江汉平原湿地的调蓄洪水能力,减轻洪涝灾害损失;其三是建立江汉平原水环境安全调控体系;最后将江汉平原湿地水系统管理纳入整个长江中游及长江流域的统一管理与规划之中,进行总体协调与管理图3参(盛春蕾)691/.:.7Cha P41200701012917。

美国科罗拉多河高速水流实验对下游水生生态系统的影响及启示袁文洁;潘明麒;于秀波【摘要】The adaptive management of the Colorado River in USA is well known worldwide. US government agencies and scientists have conducted three experiments on high-speed water flow to study the impacts of Glen Canyon Dam on the aquatic ecosystems in the lower reaches of the river. These experiments have not only generated significant ecological benefits, but accumulated a wealth of experience in this field. The research showed that the participation of stakeholders is the key to the success of these experiments institutionally, while monitoring and assessment provide strong technical support for the successful execution of such experiments, and dam adaptive regulation serves as an effective tool to protect and manage the river. China can learn a lot in terms of dam regulation from the execution processes and outcomes of the experiments on the Colorado River. Protecting the aquatic wildlife and wetland ecosystem in the lower reaches of a river should be integrated as one of the objectives of dam regulation. In addition, experiments and researches on dam adaptive regulation should be enhanced, and the relevant stakeholders should be more involved in the river protection and management.%科罗拉多河的适应性管理处于世界领先水平.美国为研究格伦峡谷大坝对下游水生生态系统的影响而进行了3次高速水流实验,取得了显著的生态效益,也积累了许多宝贵经验.分析显示,利益相关方参与是该实验成功的体制保障,跟踪监测与评估是实验成功的技术保障,闸坝适应性调度是河流保护与管理的有效手段.该实验的实施过程和效果对我国闸坝调度有许多有益的启示.例如,对下游水生生物和湿地生态系统的保护应作为闸坝调度的目标之一,应加强闸坝适应性调度的实验和研究,促进利益相关方参与河流保护与管理.【期刊名称】《湿地科学与管理》【年(卷),期】2012(008)002【总页数】5页(P49-53)【关键词】高速水流实验;科罗拉多河格伦峡谷大坝;水生态系统;适应性管理;利益相关方【作者】袁文洁;潘明麒;于秀波【作者单位】山西大学环境与资源学院,山西太原030006;中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室,北京100101;中国科学院科技政策与管理科学研究所,北京100190;中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室,北京100101【正文语种】中文美国科罗拉多河格伦峡谷大坝对下游生态系统的影响受到普遍关注。

从美国科罗拉多河看黄河朱尔明刘斌摘要美国科罗拉多河水资源开发利用与管理，可以说是美国乃至世界的一个典范，其河流特性与水资源状况与黄河有较大的相似性。

本文在比较两河异同，总结科罗拉多河水资源开发与管理的经验和不足后指出：应系统、全面、综合、辩证地看待黄河水资源问题，着眼于经济社会的大系统和水资源自身的小系统；逐步建立健全流域水资源法制体系，加强水资源全流域统一管理；利用政策法规和经济杠杆等引导和鼓励流域水资源的节约与保护；流域外调水要与黄河流城内挖潜相结合，适时合理开发；必须加强对黄河水资源的合理开发利用及节约与保护的宣传。

关键词黄河水资源管理利用科罗拉多河一、美国科罗拉多河与黄河的比较从河流特性及水资源开发的角度看，科罗拉多河与黄河有极其相似的一面，如水资源分配极不均匀，河流含沙量高，水资源相对短缺等等，从表1、表2我们会有更清晰的概念。

表1河流水文特征河流名称科罗拉多河黄河河长2333km5464km流域面积66.8×104km275.3×104km2年平均径流量185×108m3（立佛里）580×108m3（河口）年平均径流深27.7mm77.0mm 汛期3个月径流量70%60%百分比河流含沙量11.6kg/m335.7kg/m3年输沙量1.8×108t（立佛里）16×108t（花园口）年径流量丰枯比 4.0 3.4表2 水资源开发与利用河流名称科罗拉多河黄河年平均径流量185×108m3580×108m3供水区人口2500万人15300万人供水区人均水资源量740m3379m3供水区灌溉面积1800万亩11000万亩流域内总库容760×108m3（有效库容）574×108m3引水量165亿m3（1996年）395亿m3（1997年）水资源利用率89%68%农业用水比例80%80%供水区人均用水量660m3258m3从表中可以看出，科罗拉多河流域面积与黄河流域相近，也是一条多沙河流。

《泥沙研究》工程师职称文章【期刊简介】《泥沙研究》主要刊登泥沙与有关方面的学术文章和科研成果，交流技术经验，开展学术讨论，报道学术动态。

其专业范围包括：流域产沙、泥沙利用、水土保持、泥沙运动基本理论、河流动力学、河流地貌、河床演变、河流及航道整治、河口港湾泥沙、渠系泥沙、电站泥沙、水库泥沙、工业取水防沙、管路固体输送、泥石流、风沙、泥沙淤积与生态影响，模型试验及测试仪器与技术及其它工业建设中的泥沙问题。

【主要栏目】文章、短文、学术综述杂志优秀目录参考：塔里木河干流河道水沙数学模型及其应用郭庆超，陆琴，张玉江，邓安军，GUO Qing-chao，LU Qin，ZHANG Yu-jiang，DENG An-jun黄河内蒙古河道冲淤及同流量水位的变化特点吴保生，王永强，王彦君，周丽艳，WU Bao-sheng，WANG Yong-qiang，WANG Yan-jun，ZHOU Li-yan消落期三峡水库库尾泥沙冲淤特性研究袁晶，许全喜，董炳江，朱玲玲，YUAN Jing，XU Quan-xi，DONG Bing-jiang，ZHU Ling-ling长江口横沙通道演变对北槽深水航道上段回淤的影响郭兴杰，程和琴，计娜，胡浩，GUO Xing-jie，CHENG He-qin，JI Na，HU Hao下荆江二元结构河岸崩退过程模拟及影响因素分析张翼，夏军强，宗全利，许全喜，ZHANG Yi，XIA Jun-qiang，ZONG Quan-li，XU Quan-xi里下河地区沿海主要港道冲淤规律及影响要素探析朱建英，李春华，陆体成，ZHU Jian-ying，LI Chun-hua，LU Ti-cheng近十年灰鳖洋动力地貌演变特征研究刘光生，吴修广，LIU Guang-sheng，WU Xiu-guang长江口北槽洪季走航断面水沙分布特征及通量研究谢火艳，张国安，何青，李占海，XIE Huo-yan，ZHANG Guo-an，HE Qing，LI Zhan-hai人类活动驱动下伶仃洋洪季大潮水沙异变谢丽莉，刘霞，杨清书，黄广灵，XIE Li-li，LIU Xia，YANG Qing-shu，HUANG Guang-ling变坡陡比降河道强输沙下泥沙淤积与水位激增的试验研究李彬，郭志学，陈日东，张斌，LI Bin，GUO Zhi-xue，CHEN Ri-dong，ZHANG Bin珠江河口大型涉水工程方案优化研究——以港珠澳大桥工程为例何用，何贞俊，徐峰俊，王华，HE Yong，HE Zhen-jun，XU Feng-jun，WANG Hua中级职称文章发表：整体叶盘机器人自动化抛光技术【摘要】整体叶盘技术是一种在航天航空领域应用较为广泛的技术，其技术的核心在于结构与技术的综合。

美国科罗拉多河三角洲人造洪水再研究2014年5月，来自科罗拉多河的人造脉冲洪水注进加利福尼亚海湾，这场历时8个星期的短期洪水，具有持续的生态和社会效应。

早在2014年3月23日，位于美国―墨西哥边境亚利桑那州尤马附近的莫雷洛斯大坝就打开了。

在接下来的8个星期里，大坝里的水持续流入干旱的科罗拉多河河床，润泽着该河的三角洲地区，这场人造的洪水，像在大坝截流之前的天然洪水一样冲刷着科罗拉多下游的河槽，业内人士称之为人造脉冲洪水。

科罗拉多河的部分河段，由于常年干旱，一些青少年从来没有见过这样的情形，脉冲洪水的到来，科学家终于可以用数据记录这段河流了。

在2016年的10月19日，一个国际性的研究小?M公布了一项研究成果，其研究的主要内容是脉冲洪水对三角洲生态系统的影响效果。

研究人员发现，脉冲洪水对三角洲生态的影响是迅速的。

研究人员惊奇地发现，在脉冲洪水期间，许多干旱期多年没有发芽的棉木和柳树重新发芽，引起了许多生物学家的研究兴趣。

脉冲洪水带来许多土壤中的无机盐，有助于本地树种的生长。

植物的繁盛又吸引了很多鸟类成群结队前来筑巢，例如一些当地特有的啄木鸟、棕色头的牛鹂、灰喉捕蝇鹊等，还有一种黄胸鸟，也十分罕见。

河口三角洲植被的恢复使当地19种关键鸟类的数量增加了49%。

人造脉冲洪水的水文效应也相继显现，有近95%的水渗入地下含水层，补充了地下水的储量。

水流到达加利福尼亚海湾，距离莫雷洛斯大坝有124公里，是自1997年河道断流以来的第一次。

人造脉冲洪水有着明显的社会效益和政治效益，人造脉冲水流是美国和墨西哥政府之间的协议，目前一个新的协议正在酝酿中，新一轮关于环境流量的讨论将会提上日程。

研究人员希望，脉冲洪水试验能够促使新的分水协议分配更多的水到科罗拉多河三角洲。

新的研究，注重2015年12月以来科罗拉多河口三角洲地区水生态的变化，这期间包括两个生长季节，这些研究成果对科罗拉多河三角洲水生态系统的恢复意义重大，因为这里是一些国际性鸟类和几种濒危物种重要的栖息地。

少量相差21%。

造成这个差值的部分原因可能是只包含松散的泥沙。

表1　实测与计算库容的比较高程/m库容/103m3(1973)库容(1989)/103m3库容减少/103m3实测模拟实测模拟10184912161027633822202554824300237981248157025515304946749026206325044　结语为估计八堂大坝建成后水库内的河床演变及随之而引起的库容减少,开发了三维输沙模型和拟三维水工模型。

被模拟河床的总演变图与大坝建成后16a中实测结果十分吻合。

八堂水库的地形变化导致了泥沙净淤积。

模拟结果的精度表明本模型可以用来作为估计浅水河床演变的一个工具。

俞瑞堂　译自《第十九届国际大坝会议论文集》 第Ⅲ卷,Q74,R19.1997.5欧阳常恒　校(收稿日期1998-03-10)水资源拉斯维加斯盆地的水资源保护[美]　R.莫里斯　D.德维特　A.克赖茨等 摘　要　内华达州南部是美国西部发展最快的城市集中地区之一。

到2010年,水资源将近告罄。

拉斯维加斯建市初期,人均用水量高于干旱的西部其他大城市。

其原因是人们以为水资源取之不尽,政府鼓励打井,满足城市发展需求,吸引游客;加之水资源执法不严,40年代,第一次水资源危机时,采取过度开采地下水的办法应付,地下水过度开采带来地质问题。

70年代,第二次水资源危机,是通过增加从科罗拉多河取水解决。

政府采取了用水配给、制订水利法规、加强研究和教育等措施,解决水资源保护问题。

主题词　水资源保护　水法　环境教育 自由词　拉斯维加斯　美国1　引言1.1　地理条件拉斯维加斯盆地位于内华达州与犹他州大盆地南端的莫哈韦沙漠地区。

当地居民使用地下水和按计划分配的科罗拉多河水。

盆地东缘离科罗拉多及米德湖以西8.1km,东西延伸24.2km,面积1550km2。

四面环山,主峰海拔3632.6m。

盆地西北高,东南低。

在都市化前,地表水来自夏季山洪暴发、冬季降雨以及喷泉水,形成湿地,流入米德湖。

美国科罗拉多河流域的水资源分配与治理美国科罗拉多河流域是一个重要的水资源来源，在这个地区，水资源分配与治理成为关键问题。

本文将探讨科罗拉多河流域的水资源分配与治理的挑战和解决方案。

科罗拉多河流域是一个重要的水资源来源，它涵盖了美国西部八个州的大片土地。

这个地区的干旱气候和不断增长的人口使得水资源分配和治理变得尤为重要。

然而，由于不同州和利益相关者之间的竞争，水资源的分配和管理成为一个困难的问题。

首先，科罗拉多河流域的水资源分配面临着地理和地缘政治的挑战。

由于科罗拉多河流经的州份众多，每个州都希望获得足够的水供应来满足其发展需求。

然而，由于上游州的水资源供应相对较多，下游州往往面临供水不足的问题。

这导致了水资源的分配不均衡和利益冲突。

其次，科罗拉多河流域的水资源治理面临着技术和管理挑战。

水资源的存储、分配和管理需要先进的技术和有效的管理措施。

然而，由于科罗拉多河流域的地理复杂性和多个州之间的合作问题，水资源的管理变得更加困难。

加上政府部门和私人利益相关者之间的不一致，水资源治理的效率和公平性都受到了挑战。

面对这些挑战，科罗拉多河流域的水资源分配与治理需要综合解决方案。

首先，各州和利益相关者之间应加强合作与协调。

通过开展对话和谈判，各方可以找到互利共赢的解决方案，以平衡不同州的需求和实现水资源的公平分配。

其次，推动科技创新和可持续发展也是解决水资源分配和治理问题的重要途径。

技术的进步可以提高水资源的利用效率和节约用水。

例如，利用先进的灌溉系统和水资源循环利用技术可以减少用水量和节约水资源。

此外，推动可持续发展也有助于减少对水资源的需求，例如鼓励水资源的保护和生态系统的恢复。

另外，加强政府的监管和管理是确保水资源分配和治理的关键。

政府可以制定更加严格的水资源管理政策和规定，确保水资源的公平分配和合理利用。

此外，加强监测和监控系统也可以及时发现和解决水资源的问题，以保护科罗拉多河流域的水资源。

最后，公众的参与和教育也是解决水资源分配和治理问题的重要因素。

科罗拉多河水资源管理
刘斌;姚文广
【期刊名称】《海河水利》
【年(卷),期】2000(000)006
【总页数】2页(P43-44)
【作者】刘斌;姚文广
【作者单位】水利部,北京;水利部,北京
【正文语种】中文
【中图分类】TV213.4
【相关文献】
1.科罗拉多河水权分配的启示 [J], 杨彦明
2.从科罗拉多河河水之争看美河水资源分配 [J],
3.严重持续干旱对美国科罗拉多河水资源的影响及政策选择评价 [J], 杨宏伟;孙桢;黄微波
4.美国科罗拉多河水沙变化分析 [J], 吴保生;褚明华;府仁寿
5.科罗拉多河水权分配历程及其启示 [J], 周婷;郑航
因版权原因，仅展示原文概要，查看原文内容请购买。

Liverpool湾水沙输运及长期地形演变数值模拟中期报告本研究旨在利用数值模拟方法研究Liverpool湾的水沙输运及长期地形演变。

本文报告了中期研究成果，包括模型建立、参数设定和模拟结果。

一、模型建立本研究采用了基于三维流体动力学（CFD）的水沙输运和地形演变模型。

模型采用了隐式有限体积方法（IFVM）求解粘性不可压缩流动方程组。

在模拟过程中，考虑了浮力、湍流、颗粒间相互作用等因素。

二、参数设定为了准确模拟实际情况，本研究对模型的参数进行了设定。

其中，涉及水流和 suspended-load sediment的参数包括粘度、密度、粗糙度、浮力等。

沉积物参数包括颗粒大小、密度、分布等。

三、模拟结果本研究的模拟结果包括水沙输运和长期地形演变。

首先，我们对模型进行了验证，并与实测数据进行了比较。

结果表明，模型的仿真结果与实测数据吻合，具有较高的可信度。

接着，我们进行了水沙输运模拟。

通过模拟，我们可以得到Liverpool湾中不同水深和底质条件下的水流和泥沙输运特征。

结果显示，在不同条件下，水流和泥沙分选、输运规律各异。

最后，我们通过模拟研究了Liverpool湾的长期地形演变规律。

结果显示，在不同条件下，湾内的地形变化迥异，有些地区出现了淤积，而有些地区则出现了冲刷。

四、总结与展望本研究基于CFD模拟方法，对Liverpool湾的水沙输运和长期地形演变进行了模拟研究。

模拟结果显示，在不同条件下，湾内的水沙输运和地形演变规律各异。

未来，我们将继续完善模型和参数设定，并结合实测数据对模型进行进一步验证，在模拟研究中更进一步提高精度和可信度。