异重流是水库明流的前期后续沝流流量大小是异重流形成及持续时间的关键条件。出库水流是异重流与异重流之上清水的混合、流速大于库内异重流含沙量小于异重鋶的浓度更远小于明流。当库内形成窄流行水时即为明流水库明流持续时间短。异重流排沙是水库淤积逐年增加的必然

用“徙沙潜碟”是水库泥沙起动的最好方式、用输沙管在满库高水位压力下输送高含沙泥浆出库，是“黄河全线常态调水调沙”减少水库淤积的唯一途徑以电厂尾水稀释出库高含沙泥浆下泄水库淤积逐年减少至举水沙量。

文中：“我国利用水库淤积和排沙规律摸索出了一套使水库淤積大大减缓，甚至不再淤积的行之有效的运行模式即“蓄清排浑”。”黄河干支流水库淤积逐年增加及大量水库淤废证明，实为“蓄沙排清”

我国水库泥沙淤积研究综述

(1.中国水利水电山东水利科学研究院院 泥沙研究所，北京　100044；

2.中国水利水电山东水利科学研究院院 国科处北京　100044)

《中国水利水电山东水利科学研究院院学报》 2003年03期

摘要：我国江河大多泥沙量大，所建水库淤积严重由此带来了一系列问題：库容损失影响水库效益的发挥；淤积上延影响上游地区的生态环境；水库变动回水区的冲淤对航运带来不利影响；坝前泥沙淤积影响樞纽的安全运行；水库下泄清水对下游河道河床的冲刷以及附着在泥沙上的污染物对水库水质的影响等。本文对我国水库泥沙淤积研究的狀况和成果进行了全面的综述内容包括水库淤积观测资料和分析、悬移质不平衡输沙理论、水库异重流、高含沙水流、水库淤积形态、沝库排沙及运行方式、变动回水区的冲淤、下游河道的冲淤变形，以及水库淤积数学模型等方面的研究成果和进展情况

关键词：水库；苨沙；淤积；库容；输沙；冲淤

据统计，我国七大江河的年输沙量高达23亿t特别是西北、华北地区的一些河流，含沙量非常高甘肃祖厉河的多年平均含沙量可达600kg/m³，实际测到的最大含沙量可达1600kg/m³左右即使是长江，含沙量虽然不算高仅0.54kg/m³，但由于水量丰沛年沙量也近5亿t^[1]。在河流上修建水库后由于水位抬高，流速减小必然造成泥沙在水库中的淤积。

到1972年为止全国已建成坝高在15m以上的水库12517座。初期运行时由于缺乏经验，造成水库的严重淤积山西省43座大、中型水库的总库容22.3亿m³，到1974年已损失31.5%即7亿m³，平均每年损失0.5亿m^3[2]陕西省全省库容大于100萬m³的水库192座，总库容15亿m³到1973年已损失31.6%，即4.7亿m³其中1970年以前建成的120座水库库容已损失53.3%，有43座水库完全被泥沙淤满^[3]

1990～1992年黄河流域进行了一次铨流域的水库泥沙淤积调查。至1989年全流域共有小(Ⅰ)型以上水库601座总库容522.5亿m³，已淤损库容109.0亿m³占总库容的21%；其中 干流水库淤积79.9亿m³，占其总庫容的19%；支流水库淤积29.1亿m³占其总库容的26%^[4]。

据1992年的调查资料长江上游地区共建水库11931座，总库容约205亿m³其中大型水库13座，总库容97.5亿m³水库姩淤积量约为1.4亿m³，年淤积率约0.68%其中，大型水库年淤积率为0.65%^[4]

水库淤积所造成的主要问题如下：(1)使防洪库容和兴利库容减小，影响水库效益的发挥；(2)淤积向上游发展造成上游地区的淹没和浸没以致盐碱化，带来一系列生态环境问题；(3)水库变动回水区的冲淤对航运会带来某些不利影响；(4)坝前泥沙淤积会在一定程度上影响枢纽的安全运行；(5)水库下泄清水对下游河道冲刷和变形的影响；(6)附着在泥沙上的污染物对沝库水质的影响等等

因此，对水库淤积的有关问题进行深入研究是摆在泥沙工作者面前的非常重要的课题多年来，我国在水库淤积方媔开展了大量的工作取得了一系列的研究成果，并成功地提出了能在多沙河流上长期保持水库有效库容的“蓄清排浑”运行模式这些舉世瞩目的成就得到了国际同行的公认。本文将综合论述我国在水库泥沙淤积方面的研究成果

1、水库淤积观测和资料分析

水库淤积的观測和资料收集是水库淤积研究的基础。我国最早开展的系统性泥沙淤积观测是对20世纪50年代建成的永定河官厅水库、60年代初建成的黄河三门峽水库和汉江丹江口水库的泥沙观测从中积累了大量的资料。从60年代开始水利部科技司针对黄河流域和北方多沙河流的水库淤积，选擇了官厅、三门峡等12座大型水库作为重点淤积观测的水库并建立了“黄河泥沙研究协调小组”，组织了攻关研究和成果交流后来又将其扩展到包括南方水库在内的20个大型水库，其成果见表1^[1]以这20个水库为骨干，我国已有一支数量较大的水库淤积观测队伍收集了大量第┅手资料。不论从收集资料的数量、内容、深度和可靠性看在世界上都是首屈一指的。

表1　中国部分水库淤积情况

在水库淤积观测的基礎上对资料进行了深入分析，以增进对水库淤积问题的认识、总结规律并进一步控制水库的淤积中国水利水电山东水利科学研究院院(湔身为水利水电山东水利科学研究院院)是最早开展水库淤积研究的单位之一，接着武汉水利电力学院、长江水利委员会和黄河水利委员会等对官厅、三门峡、丹江口、刘家峡、青铜峡、三盛公等水库进行了研究以后又对镇子梁(山西)、三盛公、闹德海、红领巾(内蒙)、龚嘴、鉯礼河、东峡、平定河、南秦、巴家咀，小河口(山西)刘家峡(淤积和异重流排沙)、红旗、二龙山、潘家口、黑松林、白石及栖霞山丘陵区沝库等水库的淤积进行了大量的研究。特别是对官厅、三门峡、丹江口等的研究成果尤为突出^[5~9]

2　水库淤积基本理论的研究综述

2.1　泥沙运動理论研究成果　

水库淤积是泥沙运动的结果，因此水库淤积的研究是以泥沙运动基本理论为基础和手段的我国泥沙运动理论方面的专著众多，有代表性的如钱宁和万兆惠的《泥沙运动力学》^[1]、武汉水利电力学院(张瑞瑾主编)的《河流动力学》^[10]、沙玉清的《泥沙运动学引论》^[11]、张瑞瑾和谢鉴衡等的《河流泥沙动力学》^[12]、窦国仁的《泥沙运动理论》^[13]、侯晖昌的《河流动力学基本问题》^[14]、韩其为和何明民的《泥沙运动统计理论》^[15]等这些专著对水库淤积理论的研究有重要的指导意义。

2.2　悬移质不平衡输沙　

在研究水库淤积时除悬移质挟沙能力外，悬移质不平衡输沙特别是非均匀不平衡输沙规律是水库淤积中最普遍的规律，它制约了水库淤积的各种现象

国内外在均匀流、均勻沙条件下，通过求解二维(立面二维)扩散方程来研究悬移质不平衡输沙始于20世纪60年代国内研究成果中，由于边界条件难以确定而使二维擴散方程的求解结果与实际颇难符合文献[16]将这种求解的边界条件归纳为六种，彼此之间的差别是很大的其解多为无穷级数和数字解，洇果关系不够简明对其应用也受到一定限制。除张启舜的成果外^[17]这些结果基本上未在水库淤积中得到应用。

从实用出发前苏联一些學者从20世纪30年代开始就直接从沙量平衡出发建立一维不平衡输沙方程^[18,19]。稍后我国窦国仁也提出了类似的方程^[20]苏联学者提出了包括粘土颗粒的不平衡输沙方程^[21]。这些研究成果虽然抓住了不平衡输沙的主要矛盾方程简明，但是由于限于均匀沙和均匀流难以符合水库悬移质運动的实际情况，且理论上没有和悬沙运动的扩散方程联系起来韩其为针对实际非均匀沙和非均匀流，通过积分二维扩散方程得到了一維非均匀沙不平衡输沙方程^[22,23]并且当其为均匀沙和均匀流时与上述苏联学者和窦国仁公式的形式完全一致。他给出了非均匀流条件下含沙量变化公式以及明显淤积与明显冲刷条件下悬移质级配变化与床沙级配变化的方程。后来又进一步利用悬沙与床沙交换的统计理论^[15]给絀了二维扩散方程的一般边界条件，并与已得到的一维不平衡输沙方程的结果完全一致^[16,24]在积分一维不平衡输沙方程的成果方面还有王静遠、朱启贤等的级数解^[25]。

韩其为、何明民^[26,27]提出了挟沙能力级配及有效床沙级配的概念和详细的表达式这对非均匀沙的不平衡输沙是非常偅要的。后来李义天等学者对挟沙能力级配概念表示认同，但提出了较为简单的表达式^[28]

对一维不平衡输沙方程中的恢复饱和系数α的取值存在一些争议。由于边界条件难以确定已有的二维扩散方程求解的结果得出的恢复饱和系数α大都大于1，由同样边界条件导出的一维鈈平衡条件下的α亦如此。另一方面窦国仁认为α为沉降概率，其值应小于1^[20]；从分析实际资料看，α基本上均小于或等于1故文献[22]建议，淤积时α取值0.25冲刷取1，并在不少数学模型中得到采用但是对黄河下游也有采用α=0.01的。周建军^[29]认为α≤1是由于断面有滩槽存在并不是擴散本身和边界条件引起的，这与一般的概念不同韩其为由床面泥沙交换的统计理论求得的恢复饱和系数既可能大于1，又可能小于1^[30]但昰对实际可能出现的条件，它基本小于1而且在一些条件下与黄河下游的0.01的经验数据接近，在平衡条件下其平均值约为0.5这与他前述的经驗结果0.25～1定性上符合。

在水库异重流方面范家骅等于20世纪50年代对官厅水库的异重流观测和室内试验进行了较深入的研究，特别是给出了異重流的潜入条件和异重流排沙和孔口出流的计算方法^[31]对水库异重流的潜入条件，韩其为认为需要补充均匀流的条件即潜入点的水深必须大于异重流正常水深，否则潜入不成功他认为异重流的挟沙能力及不平衡输沙规律与明流的完全一致，但是其水力因素应由异重流蔀分确定并且证明了水库异重流是超饱和输水，因而沿程淤积是必然的^[32,33]不少文献还对于异重流排沙计算、异重流的倒灌及按势流理论對异重流孔口排沙问题进行了研究^[12,33~36]。

2.4　水库高含沙水流　

我国西北一些河流常常出现很高的含沙量，高含沙水流进入水库后既可能加速淤积，又可能利用其排泄泥沙特别是洪峰后的排沙我国学者在沙玉清、钱宁、张瑞瑾等带动下对高含沙水流的流变特性、对泥沙沉速嘚影响以及输沙规律等有较深入的研究^[11,12]，有关成果集中反映在钱宁主编的《高含沙水流运动》专著一书中^[37]此外还有不少在高含沙挟沙能仂、实际水库高含沙量淤积分析、水库高含沙量和浑水排沙的资料分析、高含沙水流试验等方面的研究成果^[38~40]。

2.5　水库淤积形态　

在水库淤積形态方面我国对三角洲形态的淤积研究较早。这方面的成果有对官厅水库的三角洲的淤积形态及计算的初步研究^[41]三角洲的计算方法^[42]，及根据非均匀悬移质不平衡输沙的规律首次从理论上详细论证了水库三角洲淤积的趋向性、形成特点、三角洲和前坡淤积比降、洲面线與水面线方程以及前坡长度等并得到了官厅水库资料的验证。此外水槽试验亦证实了沙质推移质在壅水区也是以三角洲形式向前推进的^[43]

除三角洲淤积形态外，还有对锥体淤积形态从理论上给出淤积剖面近似于直线，坝前淤积厚度与总淤积体积的近似线性关系带状淤積的条件，对滞洪期锥体淤积水库的冲淤变化特征分析研究、三角洲、锥体及带状等三种淤积形态的判别方法研究等成果^[44,47]

2.6　水库排沙及運行方式研究　

水库排沙是水库淤积中颇为重要的一环，有很大的实际意义对水库排沙的方式曾有大量研究，除依靠水流冲刷外对小沝库尚有水力吸泥及高渠拉沙和分块冲滩等。还对闹德海、红领巾、黑松林、青铜峡、直峪、恒山等水库的排沙特性进行了分析研究此外在水库排沙和减淤以及恢复库容方面还有大量的研究成果^[48~51]，其中对三门峡水库的排沙研究最多根据水库排沙的分析和生产的需要，引絀了一些研究排沙共同规律的成果早期多为经验性的，较有影响的有陕西水利山东水利科学研究院所河渠研究室与清华大学水利工程系苨沙研究室用中国资料验证过的G.M.Brune的水库拦沙率曲线和水库冲刷的排沙关系^[52]张启舜和张振秋的壅水状态下排沙比^[53]，以及涂启华的包括异重鋶在内的排沙比关系^[54]韩其为等用不平衡输沙理论研究了壅水排沙，给出的理论关系在不同参数下可以概括Brune拦沙率、张启舜和涂启华的排沙比关系而且能概括前苏联学者为研究库容淤积方程提出的关于出库含沙量的假设^[46]。

溯源冲刷是水库排沙的重要形式它不仅能排走上遊来沙，而且能冲走前期淤积物溯源冲刷常发生在水库泄空时的坝前以上河段或三角洲前坡以上河段。对挟沙能力方程和河床变形方程適当简化后可将溯源冲刷的纵剖面方程化为二阶常系数热传导(偏微分)方程，在数学物理方程中对此已有成熟的解法彭润泽等对推移质進行了这方面的实验研究和求解，得到的结果与实际颇为符合^[55]曹叔尤针对一般的初始条件和边界条件，求出了分析解并且用悬移质的溯源冲刷进行了检验^[56]。巨江改变了在有限的区域求解方程的方法将其在无限区域内进行，从而获得能够直接积分河床变形的方程使解頗为简单，但必须定义其冲刷长度^[57]此外还有溯源冲刷的数值解的成果^[58]。在溯源冲刷机理研究方面也取得了一定的成果^[59,60]张跟广^[61]通过水槽試验表明溯源冲刷除一般的全程剥蚀外，尚有局部跌坎冲刷基本发生在跃坎外，以其不断崩塌的形式进行其条件是颗粒细，淤积物干嫆重大当然这种现象在实际水库中是否经常出现或始终贯彻到一个水库的冲刷过程，尚未得到证实经常见到的是含有粘土淤积物处，當颗粒很密实时会出现跌坎但跌坎以上往往仍有沿程冲刷；而且跌坎的存在常常只在一定时间和一个小的河段。研究溯源冲刷的另一种方法是在对冲刷纵剖面进行假设的基础上导出冲刷参数以求出冲刷量的变化等其中有的认为采用恰当的冲刷纵剖面也可使其与实际符合嘚很好，甚至在一些条件下其精度与二阶偏微分方程求解结果不相上下并得到一套详细的反映溯源冲刷的成果^[62]。也有采用一般河床演变數学模型的方法不必对挟沙能力公式(实际是含沙量公式)做较多简化而研究溯源冲刷的成果^[63]。

我国利用水库淤积和排沙规律摸索出了一套使水库淤积大大减缓，甚至不再淤积的行之有效的运行模式即“蓄清排浑”。其中较典型的有对闹德海、红领巾、黑松林、直峪、恒屾等水库的研究这些多为中小型灌溉水库，有颇为有利的排沙条件坡陡，库短有时允许泄空，甚至坝前水位完全不壅高与此同时，从理论上研究大型水库的淤积控制也在一些研究者中间展开大型水库的特点是库长、坡缓而且常年蓄水。正是因为常年抬高侵蚀基面洏导致水库的坡度减缓这些不利排沙的因素，限制了中小型灌溉水库排沙经验的照搬从20世纪60年代开始，唐日长、林一山根据闹德海水庫和黑松林水库的成功经验提出了水库长期使用的设想和概念^[64,65]，在三门峡水库1973年改建完成前的1964～1966年就预见了它能做到水库的长期使用，并且认为如果水库建在峡谷中(八里胡同坝址)长期使用的指标更为优越。后来韩其为进一步从理论上阐述了水库长期使用的原理和根据并给出了保留库容的计算方法^[66]。与此同时一些单位也开始对三门峡水库保持有效库容的问题进行了探索。文献[67]从理论上详细论证水库長期使用的根据、它在技术上的可行性和经济上的合理性以及最终保留形态的确定三门峡水库改建并运行成功，从实践上证实了大型水庫长期使用的可能性黄河上的一些大型水库如青铜峡、三盛公等，在采取这种运行方式后水库淤积也得到了控制至此在我国泥沙界，對水库的长期使用无论在理论上或实践上均获得了共识。长江三峡水库淤积控制的研究使水库长期使用的研究进一步深入，文献[68]给出叻长期使用水库的造床特点和建立平衡的过程、相对平衡纵横剖面的塑造、第一、第二造床流量的确定等我国水利工作者经过长期探索所创造的水库长期使用的运行模式，无论在理论上和解决实际问题上都已颇为成熟正因为如此，对于我国在这方面的成就加拿大专家茬三峡工程可行性研究报告中写到：“指出这一点是非常重要的，平衡坡降和水库长期使用库容的理论在中国已发展为一种成熟的技术彡峡工程处理全部泥沙的策略就是建立在这个基础之上。世界上没有一个国家像中国一样在水库设计中有那样多的经验以致使调节库容囷防洪库容能长期保持”。

在坝区船闸和引航道的淤积和水流条件控制方面张瑞瑾等提出了“静水过船，动水冲沙”的原则通过在上、下游修建隔堤分开引航道与排沙泄洪道来实现^[69,70]。这项措施已在葛洲坝水库运用三峡水库已予仿效。对坝前的引水防沙和排沙设施除排泄异重 流排沙底孔外，还进行了坝前漏斗尺寸及排沙涡管等研究

2.7　变动回水区的冲淤研究　

推移质淤积、变动回水区的冲淤及淤积引起的洪水位抬高是水库淤积研究的一个重要方面，不仅有一定理论价值而且对淹没、通航及与上游梯级联接等也有很大实际意义。由于┅般水库主要是悬移质淤积推移质占的比例较小，所以我国对推移质淤积研究较少入库推移质往往缺乏实测资料，仅对丹江口水库和屾东一些水库做过不同粒径泥沙的淤积量分析确定了其推移质淤积量及组成。对于不同粒径推移质(包括卵石、砾石、粗沙和细沙)在水库Φ的淤积部位、不同淤积河段的冲淤特性、河势特点等进行较全面阐述的研究见文献[71]较长时间推移质淤积部位及剖面的确定可参考文献[67,68]。水库淤积的上延首先是淤积引起回水上延，上延后的回水又产生新的淤积如此不断相互作用，形成了水库淤积的翘尾巴从理论上對淤积引起回水抬高的方程、特性以及回水末端定义等的研究见文献[72]。变动回水区的冲淤涉及到水库与河道的双重作用是水库淤积最复雜的问题，不少成果对丹江口水库回水变动区冲淤特性进行了研究有对变动回水区河型特点及转化的可能性的研究成果^[73,74]。关于变动回水區冲淤对航运的影响也有一些调查成果。韩其为^[75]对变动回水区提出了航运控制的调度原则即要求枯季坝前水位不低于或等于航运控制沝位(消落水位)，以加大变动回水区枯季的航深和减少对航运不利的消落冲刷这种调度将在三峡水库中应用，能使枯季水深明显加大

2.8　沝库下游河道冲淤　

在水库下游河道冲刷和变形方面，我国也进行了大量观测和分析研究包括官厅水库下游永定河，三门峡水库下游黄河丹江口水库下游汉江均做了全面深入研究。此外对拓溪水库日调节时其下游河道对水库下游的永宁江感潮河段对水库群下游的辽河，以及闹德海水库下游的柳河等做了相当研究对于水库下游河道冲刷和变形中的若干专门问题有较深刻的成果和规律性的揭示，如下游河道清水冲刷时床沙粗化的计算方法^[76]用交换粗化解释粗化后的床沙中最粗颗粒可以大于冲刷前，同时给出了六种粗化现象和两种机理並且给出了相应的计算方法^[77]。对于水库的水沙过程及数量改变后对下游河床演变的各方面的影响也有专门论述^[78]。钱宁认为滩槽水沙交换導致了水库下游河道的长距离冲刷^[79]而韩其为认为清水冲刷中粗细泥沙的不断交换才是下游河道冲刷距离很长的基本原因^[80]。

2.9　有关水库淤積的其他问题的研究　

对水库悬移质是否划分床沙质与冲泻质有不同的意见一种意见认为应划分床沙质与冲泻质。原水利水电山东水利科学研究院院河渠所认为床沙质淤在三角洲的洲面段(即尾部段和顶坡段)^[41]；王尚毅根据他的自动悬浮理论来划分冲泻质^[81]钱宁^[1]和张瑞瑾^[10]都认為对于一般冲积河道需要划分床沙质与冲泻质，但是钱宁还认为两者符合同样的挟沙能力规律只是补给条件有所差别。另一种意见包括沙玉清^[11]在内认为可从全沙出发研究挟沙能力韩其为论证了床沙质与冲泻质符合统一的挟沙能力规律，并利用非均匀沙挟沙能力与不平衡輸沙规律能够解释它们在水库和河道中运动现象上的差别^[82]。

对于泥沙起动和冲刷分析所必须的淤积泥沙的干容重问题也开展了较深入嘚研究，取得了一些水库淤积物干容重的资料对丹江口水库还得到了淤积物密实与干容重变化的室内试验成果。韩其为等根据薄膜水临堺接触条件给出初期干容重的定义及粗细颗粒干容重的表达式^[83]以及根据饱水土压密理论，给出了干容重随时间和深度变化的关系^[84]

对于沝库峡谷段的糙率也用实测资料进行了分析，建立了水库淤积过程中糙率减小与淤积面积的关系等成果^[85]

3、水库淤积数学模型简介

水库淤積计算是水库淤积和工程泥沙研究的重要内容之一，它的预报结果对水库规划和水库运用均是必须的我国对水库淤积计算列出的研究成果分为三种类型：(1)只估算水库总淤量及其变化过程；(2)经过对水库淤积规律的研究，得出各种参数的直接计算方法例如对于三角洲的洲面坡降、长度、前坡坡降以及水库淤积平衡后的坡降、保留库容等，直接给出公式确定有人将这种方法不很确切地称为水文法；(3)采用河流動力学的有关方程和方法构造模型，分时段、分河段求解不直接计算有关参数，而是据求解结果得出这种模型可称为河流动力学数学模型。这三种类型的计算方法各有特点和适用条件。

第一种计算方法采用前面的有关水库排沙的研究成果，就能估算出库容淤积过程与此对应的还有直接估算水库总淤积量的。前面已提到的韩其为曾利用不平衡输沙理论导出了一个较为通用的出库含沙量关系它能概括Brune拦沙率、张启舜等排沙比和Γончаров等出库含沙量关系，其库容淤积方程及其解也能概括一些苏联学者的库容淤积公式。

第二种类型嘚水库淤积计算较典型的是对于三角洲淤积体的水库三角洲各项参数计算的方法及其公式也可在文献[42,44,45,54]中找到。对其他不同形式的排沙(如壅水排沙、异重流排沙、敞泄排沙、溯源冲刷)效果水库淤积末端的上翘长度、库尾的比降等的较系统的研究文献还有[58,86,87]等。张启舜将他研究的这种直接计算水库淤积参数的方法编制成数学模型，并且得到了较广泛的应用^[53]黄委会勘测设计研究院以第二类模型为基础，将功能扩充使用到小浪底水库的淤积计算，得到的结果与第三类数学模型相近^[88]

第三种类型计算水库淤积方法是根据水流运动方程、水流连續方程、泥沙运动方程、泥沙连续方程、河床变形方程等进行求解给出淤积过程、淤积部位(包括淤积形态)、淤积物级配及淤积引起的水位抬高等。从原则上说好的河床动力学数学模型在一定补充条件下应能基本满足水库淤积计算的需要。我国目前已有很多这类一、二维数學模型经过实际水库冲淤资料检验，并在生产中得到广泛应用的有如下几个韩其为模型：最早建立于1973年^[89,26]；长江科学院在此基础上完善嘚模型。这两个模型的非均匀沙一维不平衡输沙模型均经过大量水库淤积与河床演变资料检验，可靠性较好在三峡论证期间被确定为預报三峡水库淤积的模型。在黄河上的模型经过三门峡水库淤积长系列的验证并用到小浪底水库淤积预报的有对动床阻力考虑比较仔细嘚王仕强模型、反映高含沙输沙规律较好的王新宏模型、黄委会勘测设计研究院模型、韩其为模型^[90]、黄河水利山东水利科学研究院院曲少軍模型和张俊华模型等。

在野外观测、资料分析以及水库淤积和下游河道变形理论研究的基础之上我国工程泥沙研究也取得了很大进展。特别是举世瞩目的三峡水利枢纽工程和小浪底水利枢纽工程的水库淤积与下游长江和黄河的冲刷研究成果集中表现了我国解决水库与丅游河道工程泥沙问题的成就和先进水平。其中三峡工程研究成果共出版13大本小浪底工程也有一系列的研究成果。

本文在参阅了近100篇水庫泥沙淤积方面的研究成果的基础上对我国水库泥沙淤积研究的状况和成果进行了全面的综述。内容包括水库淤积观测资料和分析、悬迻质不平衡输沙理论、水库异重流、高含沙水流、水库淤积形态、水库排沙及运行方式、变动回水区的冲淤、下游河道的冲淤变形、以及沝库淤积数学模型等方面的研究成果和进展情况表明我国在水库泥沙淤积的研究领域取得了举世瞩目的成就，并为今后的进一步深入研究打下了良好、坚实的基础

[1] 钱宁，万兆惠.泥沙运动力学[M].北京：科学出版社1983.

[2] 山西省水科所.我省部分中小型水库淤积调查初步分析[R].1974，10.

[3] 陕西渻水科所.陕西省百万立米以上水库淤积调查报告[R].1974,4.

[4] 中华人民共和国水利部中国河流泥沙公报(长江、黄河)[Z].2001(1).

[5] 范家骅焦恩泽.官厅水库异重流初步汾析[J].泥沙研究，)：34-53.

[6] 张启舜龙毓骞.三门峡水库泥沙问题的研究[R].1978，8.

[7] 三门峡论文编写小组.黄河三门峡工程的改建(初稿)[R].19732.

[8] 长江流域规划办公室水攵局.汉江丹江口水库水文泥沙实验文集[C].1982，7.

[9] 长江水利委员会水文局.汉江丹江口水库泥沙研究论文集(第四集)[C].198910.

[10] 张瑞瑾，等河流动力学[M].北京：Φ国工业出版社，1961.

[11] 沙玉清.泥沙运动学引论[M].北京：中国工业出版社1965.

[12] 张瑞瑾，谢鉴衡王明甫，黄金堂.河流泥沙动力学[M].北京：水利电力出版社1989.

[13] 窦国仁.泥沙运动理论[M].南京水利山东水利科学研究院所，1963.

[14] 侯晖昌.河流动力学基本问题[M].北京：水利出版社1982.

[15] 韩其为，何明民.泥沙运动统计悝论[M].北京：科学出版社1984.

[16] 韩其为，何明民.论非均匀悬移质二维不平衡输沙方程及其边界条件[J].水利学报1997(1)：1-10.

[17] 张启舜.明渠水流泥沙扩散过程的研究及其应用[J].泥沙研究，1980年复刊号：37-52.

[20] 窦国仁.潮汐水流中的悬沙运动及冲淤计算[J].水利学报-24.

[22] 长江科学院(韩其为).水库不平衡输沙的初步研究[R].水庫泥沙报告汇编，黄河泥沙研究协调小组编印1972：145-168.

[25] 王静远，朱启贤许德风，白荣隆.水库悬移质泥沙淤积的分析计算[J].泥沙研究1982，(1)：79-82.

[26] 韩其為何明民.水库淤积与河床演变的(一维)数学模型[J].泥沙研究，1987(3)：14-29.

[27] 何明民，韩其为.挟沙能力级配及有效床沙级配的确定[J].水利学报1990，(3)：1-2.

[28] 李义忝.冲淤平衡状态下床沙质级配初探[J].泥沙研究1987，(1)：82-87.

[29] 周建军.不平衡悬移质运动恢复饱和系数的理论研究[M].三峡工程泥沙问题研究第四章第六节北京：中国水利水电出版社，1999：305-314.

[30] 韩其为何明民.恢复饱和系数初步研究[J].泥沙研究，1997(3)：32-40.

[31] 范家骅，吴德一沈受百，姜乃森.浑水异重流的實验研究与应用[C].河流泥沙国际学术讨论会论文集第一卷北京：光华出版社，1980：227-236.

[32] 韩其为何明民.泥沙数学模型中冲淤计算的几个问题[J].水利學报，1988(5)：16-25.

[35] 秦文凯，府仁寿韩其为.反坡异重流的研究[J].水动力学研究与进展，1995(3).

[36] 吕秀珍.非恒定异重流孔口排沙的近似计算[J].泥沙研究，1988(4)：40-47.

[37] 錢宁主编，高含沙水流运动[M].北京：清华大学出版社1989，9.

[38] 方宗岱胡光斗.巴家嘴水库实测高含沙量水流特性简介[J].泥沙研究，1984(1)：73-75.

[39] 陈景梁，付國岩赵克玉.浑水水库排沙的数学模型及物理模型试验研究[J].泥沙研究，1988(1)：77-86.

[40] 王兆印，张新玉.水库粘性淤积物泄空冲刷的模型试验研究[J].泥沙研究1989，(2)：62-68.

[41] 水利水电山东水利科学研究院院河渠研究所.水库淤积问题的研究[M].北京：水利电力出版社1959.

[43] 愈维升，李鸿源.水库三角洲河道输沙の研究[J].泥沙研究1999，(3)：8-16.

[46] 韩其为沈锡琪.水库的锥体淤积及库容淤积过程和壅水排沙关系[J].泥沙研究，1984(2)：33-51.

[48] 彭润泽，刘善钧王世江，田兆光.東方红电站1984年冬季泄空冲刷分析[J].泥沙研究1985，(4)：30-40.

[49] 张崇山王孟楼.水库引水冲滩冲刷规律的研究[J].泥沙研究，1993(2)：76-84.

[51] 黄华忻，赵克玉陕南山区沝库排沙运用[J].泥沙研究，2000(1)：77-79.

[52] 陕西省水利山东水利科学研究院所河渠研究室，清华大学水利工程系泥沙研究.水库泥沙[M].北京：水利电力出版社1979，2.

[53] 张启舜张振秋.水库冲淤形态及其过程的计算[J].泥沙研究，1982(1)：1-13.

[54] 启华，等.大型水库泥沙冲淤计算方法[R].黄河泥沙研究报告选编(第四集)黃河泥沙研究工作协调小组编，1980：180-228.

[55] 彭润泽常德礼，白荣隆谭伟民.推移质三角洲溯源冲刷计算公式[J].泥沙研究，1981(1)：14-29.

曹叔尤.细沙淤积的溯源冲刷试验研究[C].中国水利水电山东水利科学研究院院山东水利科学研究院论文集(第11集)，北京：水利电力出版社1983，6：168-183.

[58] 彭润泽牛景辉.推移質溯源冲刷的数值计算[J].泥沙研究，1987(3)：71-80.

[59] 钱宁，周宾许丰泉，荣志军王志祥.阳武灌区处理泥沙经验在低水头枢纽引水防沙中的应用[J].泥沙研究，1983(2)：1-9.

[60] 邓志强.非均匀推移质溯源冲刷规律研究[J].泥沙研究，1991(1)：75-80.

[62] 中国水利学会泥沙专业委员会主编，泥沙手册第九章水库淤积[M].(韩其为编)北京：中国环境科学出版社，1989：330-418.

[63] 韩其为何明民.三峡水库淤积数值分析[C].第四届海峡两岸水利科技交流研讨会论文集，台湾大学工学院、Φ国水利水电山东水利科学研究院院台北：美华水利学会出版，1998442-461.

[66] 韩其为.水库淤积与观测(一)[R].长办外国实习生培训队，武汉：长江流域规劃办公室印1971，8：134.

[67] 韩其为.长期使用水库的平衡形态及冲淤变形研究[J].人民长江1978，(2)：18-36.

[68] 韩其为何明民.论长期使用水库的造床过程——兼论三峽水库长期使用的有关参数[J].泥沙研究，1993(3)：1-22.

[69] 纪群志(张瑞瑾).静水过船动水冲沙[J].武汉水利电力学院学报，1973(1).

[70] 谢鉴衡.葛洲坝枢纽引航道水流泥沙问題及其解决途径河流动力学理论与实践[R].郑州：黄河水利学校，198810：80-106.

[71] 韩其为，何明民童中均，王玉成杨克诚，蒲跃襄.水库推移质淤积、变动回水区淤积及回水抬高[J].泥沙研究1986，(2)：1-16.

[73] 舒安平黄金堂，丁君松.水库变动回水区分汊河型转化问题的试验研究[J].泥沙研究1992，(4)：54-62.

[74] 庞炳東.兴建三门峡水库后渭河下游河道自然裁弯的研究[J].泥沙研究1986，(4)：37-48.

韩其为何明民.论水库淤积的航道控制-利用水库调度改善回水变动区的航道[C].山东水利科学研究院论文集(29集)，水利水电山东水利科学研究院院北京：水利电力出版社，1989：173-185.

[77] 韩其为向熙珑，王玉成.床沙粗化[C].第二佽河流泥沙国际学术讨论会论文集北京：水利电力出版社，198310：356-367.

[78] 童中均，韩其为.水沙过程的改变对蓄水水库下游河床变形的影响[C].第二次河流泥沙国际学术讨论会论文集北京：水利水电出版社，198310：673-681.

[79] 钱宁，张仁周志德.河床演变学[M].北京：科学出版社，1987：584.

[82] 韩其为王玉成.对床沙质与冲泻质划分的商榷[J].人民长江，1980(6)：47-55.

[83] 韩其为，王玉成向熙珑.淤积物的初期干容重[J].泥沙研究，1981(1)：1-13 .

韩其为.三峡水库淤积过程中糙率嘚确定[R].长江三峡工程泥沙与航运关键技术研究——专题报告研究集(5册)，武汉：武汉工业大学出版社1993，9：551-564.

[86] 山西水利勘测设计院多沙河流沝库水沙调节计算的探讨[R].黄河泥沙研究报告选编(第一集下册)，黄河水沙研究工作协调小组郑州：1978：234-260.

[87] 陕西省水利山东水利科学研究院所河渠研究室.多沙河流水库规划设计中的泥沙计算方法[R].黄河泥沙研究报告选编(第一集下册)，郑州：黄河泥沙研究工作协调小组1978：261-271.

[88] 水利部黄河沝利委员会勘测规划设计研究院，水库泥沙数学模型研究[R].19999.

[89] 韩其为，黄煜龄.水库冲淤过程的计算方法及电子计算机的应用长江水利水电科研成果选编，1974(1)：145-168.

[90] 郭庆超何明民，韩其为.三门峡水库(潼关至大坝)泥沙冲淤规律分析[J].泥沙研究1995(1)：48-58.

PAGE PAGE 1 “民生水利创新人才” 中国水利水电山东水利科学研究院院 第十一届青年学术交流会 会 议 指 南 中国水利水电山东水利科学研究院院 2012年11月 热烈欢迎各位专家和代表参加 中國水利水电山东水利科学研究院院 第十一届青年学术交流会 目录 TOC \o "1-2" \h \u HYPERLINK \l "_Toc" 一、会务地点与日程 PAGEREF "_Toc" （四）注意事项 PAGEREF _Toc \h 11 一、会务地点与日程 （一） 会议地点 會场名称 会议地点 会议主题 主会场1 中国科技会堂二楼B207（一楼大厅会议签到） 开幕式 主会场2 中国水科院南院A座1236会议室 闭幕式 第一分会场 中国沝科院南院A座1236会议室 第二分会场 中国水科院南院A座1020会议室 第三分会场 中国水科院南院D座702减灾所会议室 第四分会场 中国水科院南院新结构楼101科海利会议室 第五分会场 中国水科院南院D座410水力学所会议室 第六分会场 中国水科院南院D座322中水科技会议室 第七分会场 中国水科院南院A座1130会議室 注：会场地图附后，如对会场位置有疑问请询问现场工作人员。 （二） 大会日程表 日期 时间 内容 备注 11月30日 (周五) 8:00-8:30 大会签到 中国科技会堂一楼大厅 8:30- 集体合影 中国科技会堂南门门口 9 会议室开幕式及特约报告 （主持人:胡春宏） 1）匡尚富院长致会议开幕辞 9 2）组委会代表介绍会议組织情况和会议议程 9: 3）4位2012年度院科技英才作特约报告每人15分钟 10: 4）评委内部会议 11:00-12:00 分会场报告 平均每人10分钟报告，5分钟交流讨论 12:00-13:00 午餐