三峡水库淤积计算预测与原型实测结果比较分析

三峡水库淤积计算预测与原型实测结果比较分析摘要在分析三峡工程蓄水运行10A以来水库淤积规律的基础上，对长江科学院在三峡工程初步设计、技术设计以及三峡后续工作规划等不同阶段的水库淤积预测成果与三峡工程蓄水运行以来库区淤积量、淤积分布以及排沙比的实测结果进行了对比分析。结果表明，各阶段预测的水库淤积规律与实测的相同，由于预测采用的入库水沙条件及水库运用方式与三峡工程运用以来实际发生的差别较大，因此预测的库区淤积量及排沙比较实测的大，如近似将预测采用的入库水沙条件换算至与实际发生的相近时，则预测成果与实测值接近，表明预测采用的数学模型及预测成果是可靠的。今后需进一步研究揭示水库泥沙运动机理，改进和完善预测模型，并随着实测资料的积累，对数学模型进行不断地率定与改进，以提高预测精度。关键词淤积量；淤积分布；排沙比；对比分析；三峡水库中图分类号TV145文献标志码A文章编号10015485201312000106RESEARCHADVANCESINTRANSLATIONALLANDSLIDESMECHANISMANDRISKMITIGATIONABSTRACTTRANSLATIONALLANDSLIDEISASPECIALMODEOFSLOPEINSTABILITYITSSLIPSURFACEISOFEXTREMELYGENTLEANGLEOREVENCOUNTERINCLINED，ANDITSINSTABILITYMECHANISMISINCONSISTENTWITHTRADITIONALLIMITEQUILIBRIUMTHEORYRESEARCHONTHETRANSLATIONALLANDSLIDEHAZARDANDTHERISKMITIGATIONMEASURESISOFGREATSIGNIFICANCEINTHISPAPER，THERECENTRESEARCHFINDINGSANDADVANCESINCHINAANDABROADAREREVIEWEDFROMFIVEASPECTSMAINFEATURESOFTRANSLATIONALLANDSLIDE，FORMATIONCONDITIONSOFTRANSLATIONALLANDSLIDE，FORMATIONMECHANISMOFTRANSLATIONALLANDSLIDE，THRESHOLDCRITERIONOFTRANSLATIONALLANDSLIDEANDRISKMITIGATIONASSOCIATEDWITHTRANSLATIONALLANDSLIDEONTHEBASISOFTHECRITICALREVIEW，PROBLEMSINTHECURRENTRESEARCHESAREEXAMINEDTHETRENDANDPROSPECTOFTHERESEARCHAREALSODISCUSSEDFROMTHEABOVEFIVEASPECTSKEYWORDSTRANSLATIONALLANDSLIDE；FORMATIONMECHANISM；RISKMITIGATIONMEASURES；THRESHOLDCRITERION三峡工程具有防洪、发电、航运、灌溉等综合效益，水库正常蓄水位175M，防洪库容2215亿M3。水库蓄水运用后库区将发生大量泥沙淤积，库容减少，变动回水区洪水位抬高，直接影响水库长期使用和综合效益的发挥。为此，在三峡工程可行性论证、初步设计、技术设计以及三峡工程后续工作规划等各个阶段，水库泥沙问题都被列为需要解决的关键技术问题之一。为此，长江科学院采用各种研究手段和方法对水库泥沙淤积问题进行了大量研究，为工程设计、建设与运行提供了科学依据。本文在分析三峡工程运行以来水库泥沙冲淤规律的基础上，重点对比分析了不同阶段库区泥沙淤积量、淤积分布及排沙比等数学模型计算预测成果与10A来实测结果的差异及原因，为数学模型的改进和预测精度的提高提供参考。1三峡水库调度方式及运行过程三峡工程于2003年6月进入围堰蓄水期，坝前水位按汛期135M、枯季139M运行4A。2006年汛后进入初期蓄水期，坝前水位按汛期144M、枯季156M运行2A。2008年汛末开始进入175M试验性蓄水期，2008年9月28日开始蓄水，起蓄水位14527M，至11月4日达到最高水位17280M，之后水位基本稳定在170M左右运行。2009年1月1日至6月20日，坝前水位由16917M逐渐消落至14531M，汛期坝前水位基本控制在14491465M之间；9月15日开始蓄水，起蓄水位14587M，至11月24日坝前水位达到17143M，为2009年试验性蓄水最高蓄水位；汛期8月初水库进行了一次防洪运用，坝前水位达到15353M8月9日，拦蓄洪量427亿M3。2010年1月1日至6月10日，坝前水位由16939M消落至1465M以下；汛期6月10日至9月9日水库进行了7次防洪调度，最高库水位16102M，累计拦蓄洪量2643亿M3，汛期平均库水位为15154M；9月10日水库开始蓄水，起蓄水位1602M，至10月26日坝前水位首次蓄水至175M，之后库水位维持在1745175M之间。2011年1月1日至6月14日，坝前水位从17466M逐步消落至14566M，汛期水库实施了4次中小洪水调度，拦蓄洪量24716亿M3，最高运行水位为15384M；9月10日水库开始蓄水，起蓄水位15224M，至10月31日库水位蓄至1750M。2012年1月份坝前水位即开始从17467M逐步消落，至6月14日坝前水位下降至最低14539M，汛期水库实施了4次洪水调度，拦蓄洪量200亿M3，最高运行水位为16311M，其中7月24日三峡水库迎来蓄水成库以来的最大洪峰71200M3/S，经水库削峰拦洪，最大下泄流量44100M3/S；9月10日，三峡水库开始蓄水，起蓄水位15892M，至10月30日蓄水位达到1750M。此外，三峡水库还开展了生态调度和水库减淤调度试验。图1为三峡工程蓄水运用以来坝前水位过程图。可以看出，蓄水后三峡水库实际运行方式与初步设计确定的运行方式有所不同一是初期蓄水期较短；二是试验性蓄水期的汛期实施了中小洪水调度，坝前水位抬高较多；三是试验性蓄水期汛后蓄水时间提前。2三峡工程运用以来库区泥沙淤积情况21三峡工程运用以来入库水沙条件2003年至2012年，三峡工程平均年入库水量3680亿M3，较蓄水前多年平均值3972亿M3减少735，其中2006年、2011年为小水年，2005年、2012年为丰水年，这10A来水总体偏枯；平均年入库沙量19亿T，与蓄水前多年平均值46亿T相比减少了59，其中2006年和2011年入库沙量最少，仅为10亿T，入库沙量最多的2005年也只有254亿T。总的来看，三峡工程蓄水以来的10A，入库水量较蓄水前有所减少，入库沙量减少较多，只有蓄水前多年平均值的411。三峡水库入出库水沙及库区泥沙淤积量见表1。22三峡工程运用以来库区总淤积量三峡水库2003年6月蓄水至2012年12月，入库悬移质泥沙19008亿T，出库悬移质泥沙464亿T黄陵庙站，不考虑三峡库区区间来沙，水库淤积泥沙14368亿T，近似年均淤积泥沙144亿T见表1。按地形法统计，三峡水库2003年3月至2012年10月，干流库区铜锣峡大坝累积淤积13666亿M3。为比较输沙量法和地形法的计算结果，利用20052010年三峡库区淤积物干密度观测成果对20032010年库区淤积量进行了换算比较2，结果表明，根据输沙量法计算20032010年寸滩大坝段淤积泥沙120亿T；根据断面法，铜锣峡大坝段淤积119亿M3，采用干密度实测值分段进行换算，得到淤积量为116亿T，与输沙量法计算结果基本吻合。由图1、表1看出，三峡库区泥沙淤积与入库水沙及坝前水位关系密切，相同入库水量条件下来沙量多，则库区淤积量相应增加；相同来水来沙条件下库区淤积随坝前水位抬高而增加。2006年汛后坝前水位按156M运行，汛期水位按144M运行，比围堰期水位高9M，库区淤积量也相应多些，如2007年的来沙量比2005年的少，而淤积量则比2005年的多；2008年9月水库按175M试验性蓄水，汛期水位、蓄水位都相应抬高，库区淤积量也随着增加，如2008，2010，2012年的来沙量比2005年的少，而淤积量则比2005年多1530。如将2003年6月至2006年12月、20072008年和20092012年分别近似地作为135M围堰蓄水期、156M初期蓄水期和175M试验性蓄水期，则3个时段的淤积量分别为4703亿，3551亿，6114亿T，分别占总淤积量的3273，2472和4255。近似年均淤积量分别为1176亿，1776亿，1529亿T。23干流库区淤积量分布三峡水库蓄水运用后，大部分泥沙淤积在常年回水区清溪场以下，水库运用至2012年末，干流库区朱沱大坝累计淤积14368亿T，其中清溪场以下淤积1321亿T，占干流库区淤积量92，清溪场以上淤积1158亿T，约占干流库区淤积量的81。具体淤积量及分布见表2。库区泥沙主要淤积在宽谷河段和深槽。宽谷段淤积量占总淤积量的90以上。随着水库蓄水位逐步抬高，回水末端逐渐上延，变动回水区相应出现泥沙淤积，但淤积量相对较少，占干流库区总淤积量比重较小。24出库排沙比三峡水库2003年6月蓄水运用至2012年末，年平均排沙比为244。实测结果表明，蓄水后年际间出库排沙比变化较大，除受来水来沙条件的影响外，蓄水方式变化及汛期水位、蓄水位抬高也有直接的影响。围堰蓄水期坝前水位维持在139135M运行，出库排沙比平均为370，其中2006年因入库水小沙少，特别是汛期基本无大流量，入库流量大于25000M3/S仅4D，最大流量只有29700M3/S，致使出库泥沙很少，排沙比仅为87，其它年份在3841之间；初期蓄水阶段，水库按156144M方式运行，汛期水位、蓄水位均相对围堰蓄水期抬高，排沙比大幅度减小，约1523，平均为188；2008年汛后水库按175M试验性蓄水，除2011年小水少沙外，出库排沙比在1421之间，由于汛期实施中小洪水拦洪，排沙比相对正常运行情况小，其中2010年拦洪次数较多，排沙比较小，试验性蓄水期平均排沙比为161。3三峡水库泥沙淤积预测成果自三峡工程论证以来的不同阶段，长江科学院利用自行研制开发的数学模型对三峡水库不同方案的泥沙淤积进行了大量计算预测，为工程设计、建设与运行提供了科学依据。本文选取三峡工程初步设计、技术设计和三峡后续工作规划阶段正常蓄水位175M、汛期限制水位145M和枯季消落低水位155M方案计算预测的代表性成果与蓄水以来实际淤积情况进行比较分析。31计算条件311计算范围干流自合江至三斗坪，全长779KM；嘉陵江自合川至入汇口，长927KM；乌江自武隆至入汇口，长686KM。312计算初始地形初步设计阶段计算初始地形采用1960年实测地形，技术设计阶段采用2003年3月实测地形，三峡后续工作规划阶段采用2007年11月实测地形。313计算水沙条件初步设计阶段选用19611970年干流寸滩站和乌江武隆站实测水文系列简称60系列，该系列包括了丰水丰沙、中水中沙、少水少沙等不同典型年。这10A的水、沙量平均值与多年19531984年平均值接近见表3。技术设计阶段除采用60系列计算外，还采用上游来沙减少的19912000年系列简称90系列进行了计算，该系列寸滩武隆站年均水量较60系列减少71，而沙量则减少了259见表3。鉴于近期长江上游来沙大幅度减少及金沙江溪洛渡、向家坝等大型水库即将建成运行，三峡后续工作规划阶段计算时前6A2007年10月至2012年12月入库水沙条件采用近期20012006年水沙系列，该系列寸滩武隆站年均水量较60系列减少12，而沙量则减少了571见表3；2013年以后入库水沙条件采用90年代水沙系列，并考虑三峡水库上游建库拦沙。314计算时段划分时段划分以基本能反映年内水沙变化过程为原则。每个水文年划分为80个时段，其中汛期流量大或变化大时取1D为1个时段，其它根据水位、流量变化情况划分时段，枯期最长的时段为10D。315水库调度方式三峡水库采取“蓄清排浑”的方式运用，水位按175，145，155M方案的调度方式为汛期69月份库水位维持在汛期限制水位145M，腾空库容，以满足防洪需要，同时进行排沙；汛后入库沙量、水量减少，水库采用高水位175M运用，以求得较大的发电和航运效益；汛前14月份的枯水期，为了满足航运要求，库水位维持在枯水期消落水位155M以上，保证航运的畅通，同时水库泄水以满足保证发电出力和坝下游航运要求；此后逐步降低，至6月10日降至汛期限制水位145M。各阶段水库运用初期计算采用的水库调度方式为初步设计阶段水库运用初期10A坝前水位按156，135，140M方式运用。技术设计阶段水库运用初期10A蓄水位逐步上升，2003年6月至2006年9月，坝前水位按139，135M方式运用；2006年10月至2013年9月，坝前水位按156，135，140M方式运用。三峡后续工作阶段，水库运用初期2007年10月至2008年9月按156，144M方式运行1年；2008年10月至2013年12月水库按正常蓄水位175，145，155M方式运用，入库水沙条件采用近期20012006年水沙系列年。32计算预测成果由于三峡水库运用后可用于比较的资料只有10A，因此，以下仅给出各阶段公布的水库运用初期的预测成果。321初步设计阶段预测成果初步设计阶段计算预测水库运行前30A水库淤积857亿M3，淤积绝大部分在涪陵以下的常年回水区，占总淤积量的953。计算预测水库运用初期10A平均排沙比为299，库区淤积量按入、出库沙量统计，年入库沙量509亿T寸滩武隆，19611970年平均值，见表3，年出库沙量按排沙比299计，约152亿T，库区淤积量357亿T。322技术设计阶段预测成果入库水沙按60系列计算，水库运用10A末，干流库区淤积289亿M3，其中涪陵以下常年回水区淤积量占干流库区总量的93，涪陵以上淤积量占干流库区总量的7；入库水沙90系列计算，水库运用10A末，干流库区淤积1677亿M3，其中涪陵以下淤积量占干流库区总量的945，涪陵以上淤积量占干流库区总量的55。预测成果详见表44。323三峡后续工作规划阶段预测成果计算结果表明，三峡水库运用6A2007年10月至2013年12月，库区总淤积量748亿M3，其中变动回水区涪陵以上淤积051亿M3，占库区总淤积量685。4三峡水库淤积计算预测成果与蓄水运用以来10A实测结果比较本文比较所采用的水库淤积计算预测成果均摘自长江科学院各阶段公开发表的研究报告24，三峡水库蓄水运用以来10A实测结果直接引用长江水利委员会水文局公布的监测成果1。41淤积量比较三峡工程初步设计阶段计算预测水库运用初期10A，入库沙量年均509亿T，出库排沙比为299，库区年淤积量357亿T，是实测值144亿T的248倍。技术设计阶段计算预测水库运用初期10A，入库沙量为60系列时，三峡水库淤积289亿M3，年均淤积289亿M3，是实测值137亿M3的211倍；入库沙量为90系列时，库区年淤积1677亿M3，较实测值大227。三峡后续工作规划阶段入库沙量为近期20012006年系列时，计算预测三峡水库运用至2012年末2007年汛末至2012年库区淤积量748亿M3，年均淤积1496亿M3，较实测大94。综上可见，不同阶段预测的库区淤积量较实测值大，主要受入库沙量变化的影响，淤积量与入库沙量呈正相关，当入库沙量与实际蓄水的来水来沙环境相近时，库区总淤积量与实测值接近。42淤积分布比较三峡水库蓄水运用至2012年，大部分泥沙淤积在常年回水区，按输沙量法统计，涪陵以下常年回水区淤积量占干流库区的92，涪陵以上变动回水区淤积量占干流库区的8。初步设计阶段60年代水沙系列预测，水库运用30A末，涪陵以下淤积量占干流库区的95，涪陵以上淤积量占干流库区的5。技术设计阶段60水沙系列预测，水库运用10A涪陵以下淤积量占干流库区的931，涪陵以上淤积量占干流库区的69；90水沙系列预测，水库运用10A涪陵以下淤积量占干流库区的945，涪陵以上淤积量占干流库区的55。三峡后续工作规划阶段近期水沙系列预测，涪陵以下淤积量占干流库区的93，涪陵以上淤积量占干流库区的7。预测与实测对比如表5所示。结果表明，三峡水库运用初期10A计算预测的水库淤积分布与实测的分布基本接近，变动回水区淤积量占总淤积量的比例预测较实测略小。43出库排沙比比较初步设计阶段预测水库运用10A年均出库排沙比为299，实际运行10A的排沙比均值为244，可见初步设计阶段预测出库排沙比较实际运行的大55个百分点。44影响水库泥沙淤积预测精度的主要因素综上比较可见，三峡水库库区总淤积量预测值比实测大009148倍，淤积分布及排沙比也存在一定差别。影响水库泥沙淤积预测精度的因素主要有对水库泥沙输移规律的认识、边界条件及水库运用方式等3个方面。河流上修建水库后，库区形成了与天然河流完全不同的水流条件，水深大幅增加，水流从库尾天然状态至坝前几乎静止，这种大水深强不平衡水流条件下的泥沙输移规律与天然河流有很大差别，如水流和泥沙传播规律变化、细颗粒泥沙絮凝作用、水流挟沙力、变动回水区推移质泥沙运动规律、坝前段异重流形成与运动规律等。这些基本问题虽然以往也进行过许多研究，但认识还不够深入，因此，以往水库淤积预测模型中尚未考虑细颗粒泥沙絮凝作用、异重流等对水库淤积的影响，水流挟沙力等多为基于天然条件获得计算公式，对参系数做一些改进、率定，这些都对预测精度带来影响，需要结合观测资料进一步揭示水库泥沙运动机理，改进和完善预测模型。预测采用的水沙条件、河床泥沙组成、库区地形等边界条件对水库淤积预测结果影响较大，尤其是入库水沙条件。三峡水库预测采用的入库水沙条件与近10A的实际条件有较大变化，尤其是入库沙量差别较大。预测采用20世纪60年代水沙系列年，入库沙量较三峡水库实际蓄水以来10A的来沙量大168倍，其库区淤积量比实测大148倍；预测采用20世纪90年代系列年，入库沙量较实测大098倍，库区年淤积量差别相应缩小，库区淤积量预测值较实测值大227；预测入库沙量采用近期水沙系列年，入库沙量较实测大149，库区年淤积量预测值较实测值大94。因此库区淤积量与入库沙量关系密切，入库沙量大，库区淤积量亦大；入库沙量与实测来量变化小，则库区淤积量预测值与实测值相差亦小，表明不同时期预测的库区淤积量基本上是合理的。以初步设计阶段预测成果为例，若预测的入库沙量按实测比例减少，库区淤积量将比实测值少20左右，再考虑入库水量的差别、试验蓄水期提前蓄水、汛期限制水位抬高15M左右及中小洪水调度等运行方式的差异，库区淤积量预测值与实测值差别将小于20。水库运用方式对水库淤积影响也比较大，尤其是对排沙比影响相对更大。三峡水库预测采用水库蓄水方式与实际运行蓄水方式存在较大差别。初步设计阶段水库运用初期10A坝前水位按156，135，140M方式运用；技术设计阶段水库蓄水前4A坝前水位按135M运用，后6A按156，135，140M方式运用；三峡后续工作规划阶段水库初期运用5A2007年汛末至2012年坝前按175，145，155M方式运用。而实际运行中坝前水位前4A按139，135M方式运用，中间2A按156，144M方式运用，后4A水库进行试验性蓄水，按175，145，155M方式运用。按水库淤积与蓄水方式的关系推算，库区淤积预测应比实测少，排沙比应较实测大，与实际有出入，这是预测的淤积量因入库沙量大而被掩盖，排沙比预测大于实测与上述分析吻合。这也间接说明水库蓄水方式对库区淤积量的影响相对较小，对排沙比的影响相对较大，三峡水库蓄水10A实测结果表1也说明这种趋势。水库排沙比同样也与水库运用方式和来水来沙条件密切相关。初设阶段预测计算时，水库初期运用10A所采用的调度方式为156，135，140M，与实测运行10A有较大差别，而且2010年、2012年汛期水位远高于145M，最高水位分别达到16102，16309M，较初步设计预测采用的145M抬高了1618M，是排沙比相差55个百分点的主要原因。其次是这10A中遇上2006年、2011年2个小水少沙水文年，特别是汛期流量小于30000M3/S2006年最大流量为29700M3/S，2011年最大流量为27400M3/S，汛期“排浑”几率甚少，与预测的条件差别较大，致使预测的排沙比大于实测值。5结语水库泥沙输移规律的研究水平、边界条件及水库运用方式等是影响水库淤积预测精度的主要因素。三峡工程不同阶段水库淤积计算预测成果与蓄水以来10A实测结果比较分析表明，预测的水库淤积规律与实测的相同，淤积量及排沙比较实测的大，除与水库泥沙输移规律认识不够深入外，主要是预测所采用的水沙条件与蓄水以来实际的入库水沙条件差别较大，以及水库运用方式的差别所致；库区淤积量与入库沙量关系密切，入库沙量大，库区淤积量亦大，入库沙量与实测来量差异小，则库区淤积量预测值与实测值相差亦小，表明以往不同阶段的预测成果是合理可靠的，所采用的数学模型是可靠的；今后还有待需进一步揭示水库泥沙运动机理、改进和完善预测模型，并随着实测资料的积累，对数学模型进行不断的率定、改进，以提高预测精度。参考文献1长江水利委员会水文局2012年度三峡水库进出库水沙特性、水库淤积及坝下游河道冲刷分析R武汉长江水利委员会，2013BUREAUOFHYDROLOGYOFCWRCANALYSISONFEATURESOFINCOMINGANDOUTFLOWINGWATERANDSEDIMENTDISCHARGES，RESERVOIRSEDIMENTAIONANDDOWNSTREAMSCOURINGOFTGPIN2012RWUHANCHANGJIANGWATERRESOURCESCOMMISSION，2013INCHINESE2长江水利委员会水文局2010年度三峡水库进出库水沙特性、水库淤积及坝下游河道冲刷分析R武汉长江水利委员会，2011BUREAUOFHYDROLOGYOFCWRCANA