摘要: 为研究 三峡工程变动回水区的泥沙淤积.建造了长达800m的全沙试验模型.模型范围包括长江和嘉陵江约200km的天然河段.进行了清水、浑水验证以及蓄水位高程为175、180和156m的长系列模型试验.明确了各蓄水位方案下变动回水区河段的冲淤规律 、泥沙淤积对该区航道和沿江港口的影响 .为三峡工程的技术验证提供了科学 依据。
关键词: 变动回水区 全沙模型 泥沙淤积 冲淤平衡
1 引 言
长江三峡工程于1994年正式动工兴建。在此之前.对工程可行性进行过全面深入的论证。在论证工作中.直接影响可行性的一个关键问题 ,是变动回水区的泥沙淤积及其对该区航运的影响。由于这个原因,对变动回水区的泥沙淤积进行了大量的模型试验和一、二维数学模型计算 工作。由于问题的复杂性和重要性.需要建立一个变动回水区的长泥沙模型.以期对整个变动回水区的泥沙淤积及其对航运的影响作出全面深入的研究。变动回水区内上、下游河段之间有着内在的联系。下游河段的淤积将影响其上游河段的水位.从而影响上游河段的淤积量;而上游河段的淤积又将影响进入其下游河段的泥沙数量,从而影响其下游河段的淤积量。当进行河道整治试验研究时.这种上、下游之间的相互影响将更为强烈。进行变动回水区全河段长模型试验.就可较好地研究并解决这个问题。长模型的进口可置于变动回水区之上,不受囤水影响.其来沙量与天然情况下相同。模型的出口可做到变动回水区以下.位于常年水位之中。在常年回水区中,由于水面比降小.对河床糙率不敏感,因而,可由数学模型提供准确的模型出口水位。
由于泥沙运动的复杂性,在整个变动回水区长泥沙模型中准确地复演泥沙运动及冲淤变化是很困难的。我国于70年代围绕着长江葛洲坝工程泥沙的研究,开展了大规模的泥沙模型试验工作,使泥沙模型得到了迅速发展 和完善。并能在一个模型中,同时复演悬沙和底沙(包括卵石在内)的运动,从而掌握了进行全沙模型的试验技术。然而,葛洲坝的泥沙模型(包括全沙模蟹)仅限于复演较短的局部河段中的泥沙运动和冲淤变化。故对于进行整个变动回水区的长河段泥沙模型试验是否可行,必然有不少疑虑。进行长河段泥沙模型在技术上的主要困难是对模型的相似条件要求非常严格。只有各种相似条件能相应得到满足时,才有可能达到全河段各个部位的冲淤相似。因此,需要进一步提高泥沙模型试验和操作技术,以便更好地研究三峡工程变动回水区全河段的泥沙淤积情况及对该区航运的影响。
模型范围上起江津附近的青草背(航行里程①725km),下至涪陵附近的剪刀峡(航行里程550km),并包括嘉陵江18km(见图1)。自1985年按受三峡工程变动回水区全河段泥沙模型试验任务以来,完成了近800m长的模型制做、水流和泥沙冲淤验证、三峡大坝蓄水175m方案80年长系列淤积试验、水库运行100年后重庆洪水位抬高问题、蓄水180m方案80年长系列淤积试验以及175m水位方案中前期按156m水位运行30年等试验工作,为长江三峡工程的技术论证工作提供了可靠的科学依据。
①航行里程系指距离宜昌港的距离。
图1 三峡工程变动回水区河势
1、青草背 2、大中坝 3、大猫峡 4、渔洞溪 5、茄子溪 6、九龙滩 7、猪儿碛 8、重庆
9、寸滩10、铜锣峡11、明月峡 12、上洛碛13、下洛碛 14、长寿 15、黄草峡
16、金川碛 17、牛屎碛 18、剪刀峡
Fig. 1 Varying backwater zone of Three Gorges Project
2 模型设计和验证
长江水量大,沙量也大。河道迂回多弯,河床宽窄相间.坡陡流急,岸边石嘴、石梁众多,地形和流态均十分复杂。据寸滩水文站实测资料统计,泥沙年输移量约4.6亿t.泥沙粒径分布很广,从0.005mm以下直至200mm以上[1]。各种粒径的泥沙.其运动形式不同,淤积部位也不同。只有在一个模型中同时复演各种粒径泥沙的输移,才能更好地反映建库后河道各部位泥沙淤积的实际情况,因此,采用全沙模型相似理论 [2]设计泥沙运动的相似比尺。在模型设计中,除水流处于阻力平方区和满足重力相似、阻力相似,悬沙满足沉降、扬动和挟沙能力相似,底沙(包括卵石)满足起动、沉降和输移量相似外,还着重研究了悬沙和底沙级配相似。
计算表明,当模型的平面比尺λL选用250和垂直比尺λH选用100时,并采用比γs=1.46t/m3的电木粉作为模型沙,各相似比尺要求能得到较好的满足(见表1)。
表1 模型比尺汇总表
注:冲淤时间比尺λt采用120;模型沙γs=1.46t/m3;悬沙d50=0.026mm;底沙d50=0.12mm;卵石d50=3mm。
Tab. 1 Summary of Similarity scales for model
悬沙和底沙级配相似是全沙模型试验的关键。为保证原型沙和模型沙级配相似,在设计模型沙级配时,采用了文献 [3]中的统一沉降公式。将原型沙分为若干组,第i粒径组的直径为dp,i,相应沉降速度ωp,i,可由下列公式计算
其中,γs—泥沙比重;γ—水的比重;g—重力加速度;Rei—沉降雷诺数。由沉降相似要求可得到第i粒径组原型沙相对应的模型沙的沉速ωm,i,并由(1)、(2)和(3)式进一步计算得到模型沙粒径dm,i 。因此,模型沙的粒径比尺λd与原型沙的粒径有关。当原型沙的粒径范围为0.005~1.0mm时,粒径比尺的范围为1.06~2.15。
在制模中,对于关键的局部微地形亦进行了精细的塑造,保证了几何相似条件。原型河床糙率约为0.03~0.10.模型的糙率为0.022~0.074。模型河床采用梅花形排列的橡皮加糙。水面线验证试验表明,在寸滩流量为3150~21810m3/s范围内.水位误差一般在10cm以内(已换算成原型水位)个别站最大误差不超过20cm。为了验证边壁糙率,施放了寸滩站85700m3/s洪水流量,模型水位误差小于22cm。三峡建岸以后,河床将发生累积性泥沙淤积,动床糙率能否保证相似也是一个至关重要的问题。动床糙率一般由沙粒糙率和沙坡糙率所组成。R.J.Garde[4]在大量试验和原观的基础上,给出动床糙率系数Frs在0.1~1.0范围:
其中,H—水深;V—流速;d50—床沙中值粒径。在满足重力相似和采用满足相似要求的电木粉作为模型沙的条件下,由(4)式可得到λn=1.31~1.38,与阻力相似所要求的糙率比尺1.36基本一致,因此,动床的阻力相似是满足的。
在江津至剪刀峡近200km长江河道中,有大中坝、中堆、九龙滩、上洛碛、王家滩、金川碛和牛屎碛等著名宽浅滩,支流嘉陵江上还有石门和金沙碛两个宽浅段(参见图1)。在这10个宽浅段上,在一个水文年内,实测了3~9月的泥沙淤积量和9~10月的泥沙冲刷量。在模型中,模拟了实测水文年的来水来沙过程,并在模型相应10个宽浅河段,实测了3~9月的泥沙淤积量和9~10月的冲刷量。试验结果表明,模型淤积量和冲刷量与原型基本一致,误差一般在30%以内。原型河道宽浅汛期淤积、汛后冲刷,并在一个水文年内基本平衡的冲刷规律在模型中得到了较好的模拟,模型设计能反映原型河道的河床演变过程。
3 三峡工程各蓄水水位方案的泥沙淤积
三峡工程的重点论证方案为一级开发、一次建成、分期蓄水和连续移民的建设方案。在长模型中重点论证的有3种库水位运行方案(见表2)。各方案能否成立的关键,在于三峡水库按3种水位方案长期运行时变动回水区河段的泥沙淤积情况。
表2 长模型论证方案基本参数
注:(1)方案Ⅲ为方案Ⅰ初期运行方案,仅运行10年,其后按方案Ⅰ运行,但在试验中运行了30年。
(2)尾门水位由长江科学 院一维数模计算 提供。
Tab.2 Control parameters for at different selections
三峡水库蓄水后,变动回水区河道水位沿程抬高,流速普遍减小,河道的输沙能力随之降低。变动回水区河道的宽浅河段(往往是浅滩所在地)在蓄水前汛期是淤积的,汛后冲刷走沙,但蓄水后水库蓄水缩短了汛后冲刷时间,汛期淤积的泥沙在汛末不能全部冲走,变动回水区河段发生累积性淤积。泥沙淤积的数量与三峡水库水位运行方案和运行时间有关。各方案在变动回水区河段内泥沙淤积数量见表3。可见,三峡水岸按175m方案运行80年,变动回本区河段共淤积7.68亿m3;按180m方案运行80年共淤积9.76亿m3,增加2.08亿m3的淤积量。三峡水库按156m低水位运行时,30年共淤积4.36亿m3。各水位方案变动回水区河段的淤形态基率一致,差别在于淤积数量不同。各蓄水水位方案的变动回水区河段的泥沙冲淤规律 如下:
(1)三峡水库不论按何种水位运行,变动回水区均发生累积性淤积,淤积速率随水库运
用年限的增长而减缓,并在淤积过程中河道向单一、规顺、微弯和高滩深槽发展 ,并最终达到新的平衡。从图2的重庆河段在三峡水库按175m水位方案运行80年后的主要淤积部位图可见,主槽淤积较少,而边滩及回水沱或副汊则淤积较多。
(2)变动回水区淤积数量的分布与河床平面形态密切相关,宽浅河段(包括分汊河段)淤积较多,而窄深峡谷段淤积少。如三峡水库按175m水位运行80年,铜锣峡、明月峡和黄草峡每公里的淤积量仅为邻近宽浅段的8%~30%,其原因是汛期的累积性淤积主要发生在宽浅河段上。
表3 各水位运行方案变动回水区淤积汇总
注:(1)方案Ⅰ和方案Ⅱ长江段为红花碛至黄草峡,方案Ⅲ为铜锣峡至牛屎碛;
(2)嘉陵江段为入汇口至滋器口;
(3)重庆河段为李家沱至铜锣峡。
Tab.3 Summary of deposition quantities on varying backwater zone at different elevation selections
图2 重庆河段淤积形态(175m方案,80年)
Fig.2 Sketch deposition pattern on Chongqing reach(HRE 175m,80 years)
(3)淤沙粒径沿程分布的总趋势是上游河段粒径粗.越向下游粒径越细。最粗的卵石主要淤积在变动区的上端。因此,变动回水区河段的水力分选作用明显。变动回水区上端淤积相对较少,而下端淤积较多,主槽淤积较少,而边滩淤积较多。以175m水位方案为例,在7.68亿m3的总淤积量中,30%淤积在主槽中,70%淤在边滩。
(4)随着泥沙的累积性淤积,变动回水区原卵石河床逐渐为泥沙覆盖,河床糙率随之降低,水面比降也随之减小。以175m水位方案为例,建库前寸滩流量30400m3/s时,重庆至长寿河段的水面平均比降为2.0×10-4,水库运用30年、50年和80年后,其水面计算比降分别为建库前的69.0%、61.5%和61.0%;水库运用80年后,重庆以上河段的河床糙率系数相当于建库前的85%,重庆以下河段为75%。
(5)细泥沙在变动回水区河段中的造床作用不可忽略。淤沙的粒径分析 表明,各种颗粒的泥沙都参与了变动回水区的累积性淤积。以175m水位方案为例,在7.65亿m3的全部淤沙中,小于0.05mm的细沙为2.46亿m3,占总量的32.0%。在180m水位运行80年的试验中,细泥沙占更大的比例,在9.76亿m3的总淤积量中小于0.05mm的细沙为5.37亿m3,占55%。这说明,三峡水库运行水位越高,越不能忽视细泥沙的造床作用。
4 各水位方案对变动回水区河段航运的影响
三峡水库建成后,万吨级船队能否到达重庆九龙坡码头,也是三峡工程蓄水水位方案需要论证的问题 之一。试验表明,三峡水库按175m水位方案运用80年后,在水库消落期3.5m水深的最小航宽不小于150m,航道曲率半径一般均大于1000m,水流流速也较建库前大幅度降低,一般均小于2.5m/s。特别是窄深河段,如铜锣峡、明月峡和黄草峡,建库前的急流状况大大缓解,寸滩流量30400m3/s时,流速均小于2.5m/s。九龙坡码头位于变动回水区中段,九龙坡以下河道形成了一条比较稳定的深水航道,基本上满足万吨船队对航道尺寸的要求。试验过程中也发现,个别浅滩段(如洛碛)在个别枯水年的水位消落后期,3.5m水深航道宽度最小仅80m,需疏浚扩宽。某些浅滩段如九龙坡、金沙碛、金川碛的主航道在水库运用过程中发生倒槽,新航槽中的一些礁石需事先清除,以策航行安全。按180m水位方案运行80年后,九龙坡以下航道3.5m水深的最小航宽均在300m以上,航道曲率半径均大于1100m,水流流速一般均小于2.5m/s,其航道条件较175m水位方案优越,完全满足万吨船队到达重庆九龙坡码头的要求。
175m方案和180m方案都存在较严重的码头边滩淤积问题,除佛耳岩港和长寿港外,几乎所有重庆港码头、厂矿专用码头以及地方码头的前沿均出现大片边滩,将严重影响码头作业。例如在175m水位方案中,九龙坡码头前沿出现了宽约50~100m边滩(滩面高程约170~175m),原九龙坡码头作业区被淤废需要新建。由于嘉陵江入汇口的主流左摆,重庆朝天门港区嘉陵江沿岸1#~4#码头出现大片三角形边滩(最大宽度达300m,高程约170m),原码头作业区基本被淤废亦需重建。
5 175m水位方案的重庆洪水位
三峡水库长期运用后,重庆市洪水位抬高值是由一维数学模型提供的。考虑到数值中变动回水区河段糙率不易确定,加之,河道淤积数量及淤积部位对洪水位影响较大,数模成果宜在长模型中进行验证。长模型在复演重庆1981年大洪水时(寸滩流量85700m3/s),水位最大误差为0.22m,模型沙又严格遵守了各项相似比尺要求,特别是级配相似。河床淤积后,动床阻力也能满足相似要求,因此,用长模型预报三峡水库长期运用后重庆市洪水位,具有较高的精度。
由于变动回水区的下端位于长寿,模型尾门放在长寿并按一维数模的计算水位控制,并在175m方案80年淤积地形基础上进行水库运用100年的淤积试验。在100年淤积地形上分别观测了洪水频率为1%、5%和20%(洪水流量分别为88700、75300和61400m3/s)的沿江水位,相应频率的重庆洪水位分别为200.85m、197.65m和194.04m。相应频率建库前重庆洪水位分别为194.30、190.18和185.90m,即分别抬高6.55、7.47和8.14m。考虑到长模型试验的精度,洪水位的误差为±0.5m。
数学模型计算的重庆1%频率洪水位为199.09m,比长模型试验的结果偏低1.76m,5%频率洪水位偏低1.51m,20%频率洪水位偏低1.43m。数学模型所采用的长寿以上河道的综合糙率系数比长模型实际值偏小约8%~10%。为检验糙率对水位计算值的影响,在数学模型上进行了糙率敏感性分析[5]。结果表明,增、减糙率10%对常年回水区的水位和淤积量影响很小,而对变动回水区的影响较大。当糙率值增大10%时,百年一遇的重庆洪水位为201.21m,与长模型的预报上限值201.39m很接近。
6 结 语
(1)建立三峡工程变动回水区长泥沙模型不仅是必要的,也是完全可能的。只要认真把握住全沙模型相似律的基本点,就可以较好地复演近200km河道中水流和泥沙的运动规律及河床的冲淤变化,从而为全面研究 长河段泥沙问题提供新的手段。
(2)长江中各种粒径的泥沙均参与变动回水区的累积性淤积,在长河段上道行全沙试验能较好地反映河床淤积形态,从而能较好地明确泥沙淤积对变动回水区航运的影响。
(3)三峡工程不论何种水位运行方案,其变动回水区将发生累积性淤积,其淤积速率随着水库运用年限的增长而减缓,并最终达到新的平衡。在淤积发展的过程中,河道向单一、规顺、微弯、高滩、深槽演变。
(4)在变动回水区中淤积沿程分布的总趋势,是愈往上游淤积越少,但在靠近回水末端一段是粗沙卵石淤积区,淤积比较严重。在横向分布上,总的情况是边滩淤得多、主槽淤得少,但在发生倒槽河段原主槽将发生严重淤积。淤沙粒径的分布规律是上游粗、下游细,细颗粒主要淤在高滩上。
(5)随青河床的淤积,水位不断升高,但因淤沙覆盖原沙卵石河床的程度增大,河床糙率减小,水面比降也随之减小。蓄水位愈高,淤积量愈大,水位壅高愈多,水面比降亦愈小。
(6)蓄水位175m方案运用80年变动回水区河床接近平衡,运用100年基本平衡,长江九龙坡以下形成一条较好的航道,基本满足万吨级船队直达重庆九龙坡码头的要求,但现有沿江大部分码头将受到严重影响。对港口淤积造成碍航的问题应通过优化水库调度、港口改造、航道整治和疏浚等措施加以解决。
(7)蓄水位175m方案运用100年后发生百年一退洪水时,重庆水位约为200.85±0.5m,较建库前抬高约6~7m。
(8)在改善航道条件方面,180m方案优于175m方案,在长江九龙坡以下可形成一条良好航道,完全满足万吨级船队直达九龙坡码头的要求,且大大增加万吨级船队驶抵九龙坡的天数。但由于在180m方案中增加的淤积量都是小于0.04mm的细颗粒,使边滩淤高,对现有沿江码头的影响较175m方案更为严重,需要结合港口改造和整治来解决。
(9)在175m方案中前期按低水位156m运行,其前10年长江铜锣峡以下航道较建库前有一定改善。如果运用30年,则某些关键河段的航道条件已接近建库前的严重情况。因此,低水位运行阶段不宜太长。
参考 文献
1窦国仁,万声淦,陆长石.长江江津至涪陵河段水沙条件和河床演变分析.南京:南京水利 科学研究院,1989
2 窦国仁,全沙河工模型试验的研究. 科学道报,1981;(14)
3 窦国仁,紊流力学,北京:人民教育 出版社,1987
4 Garde BJ,Raju KGR. Resistance relationship fro alluvial channel flow .J of the Proc of the Am Soc of Civ Eng .Hydraulic Division ,HY4 ,1966;(6)
5长江科学院. 三峡工程水库泥沙淤积计算敏感性分析(长江三峡工程泥沙研究文集).水利 电力部科学技术司,长江科学院,北京:中国 科学技术出版社,1990:56-72
Investigation on sedimentation in varying
backwater zone of Three Gorges project
Abstract
A long river model,whose total length is about 800m,is used to simulate the varing backwater zone of the Three Gorges Project,covering 200km long of natural river reaches of the Yangtze River and Jialing River.The model is designed based of the total load similarity theory developed by first author.Verification tests show that the model has a basic similarity agreement with the prototype in the aspects of resistance characteristic,flow condition and laws governing scouring and deposition.By using this model,a series of tests with norma1 reservoir elevations at 175m,180m,and 156m with reservoir operation of 80 years and 30 years have been conducted to study the laws governing scouring and deposition in the varying backwater zone,the siltation configuration in channels and harbours for the feasibility studies of the Three Gorges Project.
Key words:varying backwater zone,total sediment transport model,siltation,fluvial equilibrium