水库淤积研究综述.doc

我国水库泥沙淤积研究综述韩其为1，杨小庆2(1.中国水利水电科学研究院 泥沙研究所，北京100044；2.中国水利水电科学研究院 国科处，北京100044)摘要：我国江河大多泥沙量大，所建水库淤积严重，由此带来了一系列问题：库容损失影响水库效益的发挥；淤积上延影响上游地区的生态环境；水库变动回水区的冲淤对航运带来不利影响；坝前泥沙淤积影响枢纽的安全运行；水库下泄清水对下游河道河床的冲刷以及附着在泥沙上的污染物对水库水质的影响等。本文对我国水库泥沙淤积研究的状况和成果进行了全面的综述。内容包括水库淤积观测资料和分析、悬移质不平衡输沙理论、水库异重流、高含沙水流、水库淤积形态、水库排沙及运行方式、变动回水区的冲淤、下游河道的冲淤变形，以及水库淤积数学模型等方面的研究成果和进展情况。关键词：水库；泥沙；淤积；库容；输沙；冲淤中图分类号：TV145文献标识码：A据统计，我国七大江河的年输沙量高达23亿t，特别是西北、华北地区的一些河流，含沙量非常高，甘肃祖厉河的多年平均含沙量可达600kg/m3，实际测到的最大含沙量可达1600kg/m3左右。即使是长江，含沙量虽然不算高，仅0.54kg/m3，但由于水量丰沛，年沙量也近5亿t1。在河流上修建水库后，由于水位抬高，流速减小，必然造成泥沙在水库中的淤积。到1972年为止，全国已建成坝高在15m以上的水库12517座。初期运行时，由于缺乏经验，造成水库的严重淤积。山西省43座大、中型水库的总库容22.3亿m3，到1974年已损失31.5%，即7亿m3，平均每年损失0.5亿m32。陕西省全省库容大于100万m3的水库192座，总库容15亿m3，到1973年已损失31.6%，即4.7亿m3，其中1970年以前建成的120座水库库容已损失53.3%，有43座水库完全被泥沙淤满3。19901992年黄河流域进行了一次全流域的水库泥沙淤积调查。至1989年全流域共有小()型以上水库601座，总库容522.5亿m3，已淤损库容109.0亿m3，占总库容的21%；其中 干流水库淤积79.9亿m3，占其总库容的19%；支流水库淤积29.1亿m3，占其总库容的26%4。据1992年的调查资料，长江上游地区共建水库11931座，总库容约205亿m3。其中大型水库13座，总库容97.5亿m3。水库年淤积量约为1.4亿m3，年淤积率约0.68%，其中，大型水库年淤积率为0.65%4。水库淤积所造成的主要问题如下：(1)使防洪库容和兴利库容减小，影响水库效益的发挥；(2)淤积向上游发展，造成上游地区的淹没和浸没以致盐碱化，带来一系列生态环境问题；(3)水库变动回水区的冲淤对航运会带来某些不利影响；(4)坝前泥沙淤积会在一定程度上影响枢纽的安全运行；(5)水库下泄清水对下游河道冲刷和变形的影响；(6)附着在泥沙上的污染物对水库水质的影响等等。因此，对水库淤积的有关问题进行深入研究是摆在泥沙工作者面前的非常重要的课题。多年来，我国在水库淤积方面开展了大量的工作，取得了一系列的研究成果，并成功地提出了能在多沙河流上长期保持水库有效库容的“蓄清排浑”运行模式，这些举世瞩目的成就得到了国际同行的公认。本文将综合论述我国在水库泥沙淤积方面的研究成果。1水库淤积观测和资料分析水库淤积的观测和资料收集是水库淤积研究的基础。我国最早开展的系统性泥沙淤积观测是对20世纪50年代建成的永定河官厅水库、60年代初建成的黄河三门峡水库和汉江丹江口水库的泥沙观测，从中积累了大量的资料。从60年代开始，水利部科技司针对黄河流域和北方多沙河流的水库淤积，选择了官厅、三门峡等12座大型水库作为重点淤积观测的水库，并建立了“黄河泥沙研究协调小组”，组织了攻关研究和成果交流。后来又将其扩展到包括南方水库在内的20个大型水库，其成果见表11。以这20个水库为骨干，我国已有一支数量较大的水库淤积观测队伍，收集了大量第一手资料。不论从收集资料的数量、内容、深度和可靠性看，在世界上都是首屈一指的。表1中国部分水库淤积情况序号水库名称河流控制面积/km2坝高/m设计库容/亿m3统计年限总淤积量/亿m3淤积量占库容/%1刘家峡黄河181700147.0 57.20196819785.8010.12盐锅峡黄河18280057.02.20196119781.6072.73八盘峡黄河20470043.00.490197519770.1835.74青铜峡黄河28500042.76.20196619774.8578.25三盛公黄河314000闸坝式0.800196119770.4050.06天桥黄河38800042.00.680197619780.07511.07三门峡黄河688400106.096.401960197837.6039.08巴家嘴蒲河352274.05.25196019781.9437.09冯家山千河323273.03.89197419780.235.910黑松林冶峪河37045.50.086196119770.03439.011汾河汾河526860.07.00195919772.6037.112官厅永定河4760045.022.70195319775.5224.313红山西辽河2449031.025.60196019774.7518.514闹德海柳河450141.51.9619420.3819.515冶原弥河78623.71.68195919720.127.216岗南滹沱河1590063.015.60196019762.3515.117龚嘴大渡河7640088.03.51196719781.3338.018碧口白龙江27600101.05.21197619780.285.419丹江口汉江95220110.0160.50196819746.253.920新桥红柳河132747.02.0014年1.5678.0在水库淤积观测的基础上，对资料进行了深入分析，以增进对水库淤积问题的认识、总结规律并进一步控制水库的淤积。中国水利水电科学研究院(前身为水利水电科学研究院)是最早开展水库淤积研究的单位之一，接着武汉水利电力学院、长江水利委员会和黄河水利委员会等对官厅、三门峡、丹江口、刘家峡、青铜峡、三盛公等水库进行了研究。以后又对镇子梁(山西)、三盛公、闹德海、红领巾(内蒙)、龚嘴、以礼河、东峡、平定河、南秦、巴家咀，小河口(山西)，刘家峡(淤积和异重流排沙)、红旗、二龙山、潘家口、黑松林、白石及栖霞山丘陵区水库等水库的淤积进行了大量的研究。特别是对官厅、三门峡、丹江口等的研究成果尤为突出59。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，由同样边界条件导出的一维不平衡条件下的亦如此。另一方面，窦国仁认为为沉降概率，其值应小于120；从分析实际资料看，基本上均小于或等于1。故文献22建议，淤积时取值0.25，冲刷取1，并在不少数学模型中得到采用。但是对黄河下游也有采用=0.01的。周建军29认为1是由于断面有滩槽存在，并不是扩散本身和边界条件引起的，这与一般的概念不同。韩其为由床面泥沙交换的统计理论求得的恢复饱和系数既可能大于1，又可能小于130。但是对实际可能出现的条件，它基本小于1，而且在一些条件下与黄河下游的0.01的经验数据接近，在平衡条件下其平均值约为0.5，这与他前述的经验结果0.251定性上符合。2.3水库异重流在水库异重流方面，范家骅等于20世纪50年代对官厅水库的异重流观测和室内试验进行了较深入的研究，特别是给出了异重流的潜入条件和异重流排沙和孔口出流的计算方法31。对水库异重流的潜入条件，韩其为认为需要补充均匀流的条件，即潜入点的水深必须大于异重流正常水深，否则潜入不成功。他认为异重流的挟沙能力及不平衡输沙规律与明流的完全一致，但是其水力因素应由异重流部分确定，并且证明了水库异重流是超饱和输水，因而沿程淤积是必然的32,33。不少文献还对于异重流排沙计算、异重流的倒灌及按势流理论对异重流孔口排沙问题进行了研究12,3336。2.4水库高含沙水流我国西北一些河流，常常出现很高的含沙量，高含沙水流进入水库后，既可能加速淤积，又可能利用其排泄泥沙特别是洪峰后的排沙。我国学者在沙玉清、钱宁、张瑞瑾等带动下对高含沙水流的流变特性、对泥沙沉速的影响以及输沙规律等有较深入的研究11,12，有关成果集中反映在钱宁主编的高含沙水流运动专著一书中37。此外还有不少在高含沙挟沙能力、实际水库高含沙量淤积分析、水库高含沙量和浑水排沙的资料分析、高含沙水流试验等方面的研究成果3840。2.5水库淤积形态在水库淤积形态方面，我国对三角洲形态的淤积研究较早。这方面的成果有对官厅水库的三角洲的淤积形态及计算的初步研究41，三角洲的计算方法42，及根据非均匀悬移质不平衡输沙的规律首次从理论上详细论证了水库三角洲淤积的趋向性、形成特点、三角洲和前坡淤积比降、洲面线与水面线方程以及前坡长度等，并得到了官厅水库资料的验证。此外水槽试验亦证实了沙质推移质在壅水区也是以三角洲形式向前推进的43。除三角洲淤积形态外，还有对锥体淤积形态，从理论上给出淤积剖面近似于直线，坝前淤积厚度与总淤积体积的近似线性关系，带状淤积的条件，对滞洪期锥体淤积水库的冲淤变化特征分析研究、三角洲、锥体及带状等三种淤积形态的判别方法研究等成果44,47。2.6水库排沙及运行方式研究水库排沙是水库淤积中颇为重要的一环，有很大的实际意义。对水库排沙的方式曾有大量研究，除依靠水流冲刷外，对小水库尚有水力吸泥及高渠拉沙和分块冲滩等。还对闹德海、红领巾、黑松林、青铜峡、直峪、恒山等水库的排沙特性进行了分析研究。此外在水库排沙和减淤以及恢复库容方面还有大量的研究成果4851，其中对三门峡水库的排沙研究最多。根据水库排沙的分析和生产的需要，引出了一些研究排沙共同规律的成果。早期多为经验性的，较有影响的有陕西水利科学研究所河渠研究室与清华大学水利工程系泥沙研究室用中国资料验证过的G.M.Brune的水库拦沙率曲线和水库冲刷的排沙关系52，张启舜和张振秋的壅水状态下排沙比53，以及涂启华的包括异重流在内的排沙比关系54。韩其为等用不平衡输沙理论研究了壅水排沙，给出的理论关系在不同参数下可以概括Brune拦沙率、张启舜和涂启华的排沙比关系，而且能概括前苏联学者为研究库容淤积方程提出的关于出库含沙量的假设46。溯源冲刷是水库排沙的重要形式，它不仅能排走上游来沙，而且能冲走前期淤积物。溯源冲刷常发生在水库泄空时的坝前以上河段或三角洲前坡以上河段。对挟沙能力方程和河床变形方程适当简化后，可将溯源冲刷的纵剖面方程化为二阶常系数热传导(偏微分)方程，在数学物理方程中对此已有成熟的解法。彭润泽等对推移质进行了这方面的实验研究和求解，得到的结果与实际颇为符合55。曹叔尤针对一般的初始条件和边界条件，求出了分析解，并且用悬移质的溯源冲刷进行了检验56。巨江改变了在有限的区域求解方程的方法，将其在无限区域内进行，从而获得能够直接积分河床变形的方程，使解颇为简单，但必须定义其冲刷长度57。此外还有溯源冲刷的数值解的成果58。在溯源冲刷机理研究方面也取得了一定的成果59,60。张跟广61通过水槽试验表明溯源冲刷除一般的全程剥蚀外，尚有局部跌坎，冲刷基本发生在跃坎外，以其不断崩塌的形式进行，其条件是颗粒细，淤积物干容重大。当然这种现象在实际水库中是否经常出现或始终贯彻到一个水库的冲刷过程，尚未得到证实。经常见到的是含有粘土淤积物处，当颗粒很密实时会出现跌坎，但跌坎以上往往仍有沿程冲刷；而且跌坎的存在常常只在一定时间和一个小的河段。研究溯源冲刷的另一种方法是在对冲刷纵剖面进行假设的基础上导出冲刷参数以求出冲刷量的变化等。其中有的认为采用恰当的冲刷纵剖面也可使其与实际符合得很好，甚至在一些条件下其精度与二阶偏微分方程求解结果不相上下，并得到一套详细的反映溯源冲刷的成果62。也有采用一般河床演变数学模型的方法，不必对挟沙能力公式(实际是含沙量公式)做较多简化而研究溯源冲刷的成果63。我国利用水库淤积和排沙规律，摸索出了一套使水库淤积大大减缓，甚至不再淤积的行之有效的运行模式，即“蓄清排浑”。其中较典型的有对闹德海、红领巾、黑松林、直峪、恒山等水库的研究。这些多为中小型灌溉水库，有颇为有利的排沙条件，坡陡，库短，有时允许泄空，甚至坝前水位完全不壅高。与此同时，从理论上研究大型水库的淤积控制也在一些研究者中间展开。大型水库的特点是库长、坡缓而且常年蓄水。正是因为常年抬高侵蚀基面而导致水库的坡度减缓。这些不利排沙的因素，限制了中小型灌溉水库排沙经验的照搬。从20世纪60年代开始，唐日长、林一山根据闹德海水库和黑松林水库的成功经验，提出了水库长期使用的设想和概念64,65，在三门峡水库1973年改建完成前的19641966年，就预见了它能做到水库的长期使用，并且认为如果水库建在峡谷中(八里胡同坝址)，长期使用的指标更为优越。后来韩其为进一步从理论上阐述了水库长期使用的原理和根据，并给出了保留库容的计算方法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、黄河水利科学研究院曲少军模型和张俊华模型等。4工程泥沙研究成果在野外观测、资料分析以及水库淤积和下游河道变形理论研究的基础之上，我国工程泥沙研究也取得了很大进展。特别是举世瞩目的三峡水利枢纽工程和小浪底水利枢纽工程的水库淤积与下游长江和黄河的冲刷研究成果，集中表现了我国解决水库与下游河道工程泥沙问题的成就和先进水平。其中三峡工程研究成果共出版13大本，小浪底工程也有一系列的研究成果。5结语本文在参阅了近100篇水库泥沙淤积方面的研究成果的基础上，对我国水库泥沙淤积研究的状况和成果进行了全面的综述。内容包括水库淤积观测资料和分析、悬移质不平衡输沙理论、水库异重流、高含沙水流、水库淤积形态、水库排沙及运行方式、变动回水区的冲淤、下游河道的冲淤变形、以及水库淤积数学模型等方面的研究成果和进展情况。表明我国在水库泥沙淤积的研究领域取得了举世瞩目的成就，并为今后的进一步深入研究打下了良好、坚实的基础。参考文献：1 钱宁，万兆惠.泥沙运动力学M.北京：科学出版社，1983.2 山西省水科所.我省部分中小型水库淤积调查初步分析R.1974，10.3 陕西省水科所.陕西省百万立米以上水库淤积调查报告R.1974,4.4 中华人民共和国水利部中国河流泥沙公报(长江、黄河)Z.2001(1).5 范家骅，焦恩泽.官厅水库异重流初步分析J.泥沙研究，1958.3(4)：34-53.6 张启舜，龙毓骞.三门峡水库泥沙问题的研究R.1978，8.7 三门峡论文编写小组.黄河三门峡工程的改建(初稿)R.1973，2.8 长江流域规划办公室水文局.汉江丹江口水库水文泥沙实验文集C.1982，7.9 长江水利委员会水文局.汉江丹江口水库泥沙研究论文集(第四集)C.1989，10.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.18 ,.M.，1959.19 ， M. ，1960.20 窦国仁.潮汐水流中的悬沙运动及冲淤计算J.水利学报，1963(4):13-24.21 . , ，124，1965.22 长江科学院(韩其为).水库不平衡输沙的初步研究R.水库泥沙报告汇编，黄河泥沙研究协调小组编印，1972：145-168.23 韩其为.非均匀悬移质不平衡输沙研究J.科学通报，1979.(17)，804-808.24 Han Qiwei, He Mingmin. A Mathematical Model for Reservoir Sedimentation and Fluvial Processes J. Inter. J. of Sediment Research,1990,5,(2):43-84.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.33 吴德一.关于水库异重流的计算方法J.泥沙研究，1983，(2)：54-63.34 金德春.浑水异重流的运动和淤积J.水利学报，1980，(3)：39-47.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.42 张威.水库三角洲淤积及其近似计算J.人民长江，1964，(2)：39-45.43 愈维升，李鸿源.水库三角洲河道输沙之研究J.泥沙研究，1999，(3)：8-16.44 韩其为.论水库的三角洲淤积(一)J.湖泊科学，1995，(2)：107-118.45 韩其为.论水库的三角洲淤积(二)J.湖泊科学，1995，(3)：213-225.46 韩其为，沈锡琪.水库的锥体淤积及库容淤积过程和壅水排沙关系J.泥沙研究，1984，(2)：33-51.47 杨克诚.滞洪水库淤积规律J.泥沙研究，1981，(3)：73-81.48 彭润泽，刘善钧，王世江，田兆光.东方红电站1984年冬季泄空冲刷分析J.泥沙研究，1985，(4)：30-40.49 张崇山，王孟楼.水库引水冲滩冲刷规律的研究J.泥沙研究，1993，(2)：76-84.50 焦恩泽.三门峡水库拦粗排细的分析J.泥沙研究，1997，(3)：62-65.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.56 曹叔尤.细沙淤积的溯源冲刷试验研究C.中国水利水电科学研究院科学研究论文集(第11集)，北京：水利电力出版社，1983，6：168-183.57 巨江.溯源冲刷的计算方法及其应用J.泥沙研究，1990，(1)：30-39.58 彭润泽，牛景辉.推移质溯源冲刷的数值计算J.泥沙研究，1987，(3)：71-80.59 钱宁，周宾，许丰泉，荣志军，王志祥.阳武灌区处理泥沙经验在低水头枢纽引水防沙中的应用J.泥沙研究，1983，(2)：1-9.60 邓志强.非均匀推移质溯源冲刷规律研究J.泥沙研究，1991，(1)：75-80.61 张跟广.水库溯源冲刷模式初探J.泥沙研究，1993，(3)：86-94.62 中国水利学会泥沙专业委员会主编，泥沙手册第九章水库淤积M.(韩其为编)，北京：中国环境科学出版社，1989：330-418.63 韩其为，何明民.三峡水库淤积数值分析C.第四届海峡两岸水利科技交流研讨会论文集，台湾大学工学院、中国水利水电科学研究院，台北：美华水利学会出版，1998，442-461.64 唐日长.水库淤积调查报告J.人民长江，1964，(3)：8-14.65 林一山.水库长期使用问题J.人民长江，1978，(2)：1-8.66 韩其为.水库淤积与观测(一)R.长办外国实习生培训队，武汉：长江流域规划办公室印，1971，8：134.67 韩其为.长期使用水库的平衡形态及冲淤变形研究J.人民长江，1978，(2)：18-36.68 韩其为，何明民.论长期使用水库的造床过程兼论三峡水库长期使用的有关参数J.泥沙研究，1993(3)：1-22.69 纪群志(张瑞瑾).静水过船，动水冲沙J.武汉水利电力学院学报，1973(1).70 谢鉴衡.葛洲坝枢纽引航道水流泥沙问题及其解决途径，河流动力学理论