1. 前言
黄土高原严重的水土流失不仅造成当地生态环境脆弱,生活贫困,而且致使黄河下游河道淤积,严重威胁着黄河下游的河道稳定和防洪安全。根据三门峡水库修建以前多年平均,每年有4亿t泥沙淤积在下游河道中,形成著名的“地上悬河” [1] 。20世纪90年代进入下游河道泥沙量(以花园口站为代表)仅为1950~1960年间平均值的46.6%。然而,由于径流量大幅度减少,黄河下游河道淤积仍较严重。1986~1999年间,下游河道年均淤积量为2.23亿t,下游各测站3000 m3/s水位年均抬升0.10~0.15 m [2] 。为此,黄河水利委员会提出了“维持黄河健康生命”的治河理念,其中安全通畅的水沙通道是黄河健康的重要标志之一 [3] 。黄河下游河道成为一条地上悬河,是水流携带泥沙长期淤积的结果,洪水携带大量泥沙进入下游河道,因水沙条件的差异在下游河道发生冲刷或淤积。国内外学者对黄河下游河道水沙关系进行了大量研究,如严军和胡春宏 [4] 根据1950~2000年黄河下游河道实测水沙资料,分析得到黄河下游河道输沙水量与单位输沙水量的计算公式;郭少磊等 [5] 以1960~1999年黄河下游水沙资料为基础,用统计方法对黄河下游不同量级洪水与不同含沙量洪水的输沙规律进行了分析,提出了实现均衡输沙的水沙条件;袁吉栋等 [6] 根据1950~2000年黄河下游河道汛期实测径流量和输沙量资料,计算得到下游河段处于冲淤平衡状态时的输沙量、输沙水量与单位输沙水量。黄河下游泥沙冲淤特征也是学者的研究热点,如胡春宏等 [7] 根据1974~1979年实测水沙资料,得到下游冲淤率为0的临界含沙量为21 kg/m3;龙毓骞等 [8] 根据1950年以来实测水沙资料研究得到黄河下游临界平衡的含沙量也接近21 kg/m3;韩其为 [9] 从黄河下游河型沿程变化的实际出发,建立河相系数与流量的关系,研究得出了黄河下游泥沙冲淤的临界流量。上述研究主要侧重黄河下游不同量级流量的输沙量及其关系,而对于黄河下游不同水文年型的冲淤平衡临界输沙量研究较少。并且,黄河上游河段的来水来沙量、不同泥沙颗粒级配以及不同河段河槽条件等均影响着河流输沙能力,如何在上述约束条件下求出河流达到冲淤平衡时的临界输沙量,对于预测下游河道的发展趋势,以及制定合理的河道治理规划和水库调度运行方式具有重要的现实意义和理论价值。因此,笔者在总结分析前人有关黄河下游冲淤平衡时的不同含沙量的临界流量成果的基础上,利用花园口、高村水文站逐日水沙资料,推算出不同水文年黄河下游水沙冲淤平衡的临界输沙量,以期为黄河下游河道治理提供科学依据。
2. 研究资料与方法
黄河下游河道的冲淤平衡与河道来水来沙量密切相关,河道水沙的平衡是一个相对状态,因此,从平均意义上来探讨黄河下游冲淤平衡。利用花园口和高村水文站逐日水沙资料(1962~1988年),通过年径流量大小将下游河道划分为不同水文年(丰水年、平水年和枯水年),选择不同水文年典型代表年分别为平水年(1976~1978年)、丰水年(1982~1985年)、枯水年(1986~1988年)。根据黄河流域降雨特性,将1年划分为汛期(7~10月)和非汛期(11月至翌年6月)。利用不同水文年逐日水沙资料,划分不同流量级,统计不同流量级出现频次及占汛期和非汛期总时间(d)的比例(Pi),以Pi作为黄河下游冲淤平衡临街输沙量的权重系数,结合前人 [9] [10] 对黄河下游冲淤平衡时不同含沙量的临界流量值,总结出不同流量级临界含沙量。在此基础上,利用加权平均法,计算出汛期和非汛期不同流量级流量的临界含沙量,再乘以汛期和非汛期径流量,即得到汛期和非汛期的临界输沙量。临界含沙量计算公式为 [11] :
(1)
式中:SC为汛期或非汛期的临界含沙量,kg/m3;Sci为i流量级的临界含沙量,kg/m3;Pi为i流量级出现频次占某时段(汛期或非汛期)的比例,%。
3. 黄河下游冲淤特征
3.1. 三门峡水库不同运用时期的冲淤特征
1960年三门峡水库开始蓄水运用,黄河下游来水来沙的变化使得河道的冲淤情形也随之发生改变 [12] 。
(1) 蓄水运用期。至1964年10月期间共拦沙44.10亿t,下游河道平均含沙量10.28 kg/m3,为1960~1996年平均含沙量的38.0%。在此期间下游河道共冲刷泥沙27.80亿t,花园口至利津河段冲刷量占全下游冲刷量的66.0%。
(2) 滞洪排沙期。1964年11月至1973年10月,水库采用滞洪排沙运行方式,平均含沙量达到38.33 kg/m3,为1960~1996年含沙量平均值的1.43倍,期间河道共淤积泥沙39.42亿t,花园口至利津河段淤积量占全下游的77.6%。
(3) 蓄清排浑运用。1973年11月至1996年6月,下游河段平均含沙量25.92 kg/m3,为1960~1996年平均值的97%;下游淤积有所减弱,淤积泥沙24.70亿t,其中花园口以上冲刷0.61亿t,淤积主要发生在花园口至艾山河段。
由上述分析可知,除三门峡水库运用初期蓄水运用下游河道发生冲刷外,其余两个运用期,花园口至艾山河段均是主要淤积河段,淤积量分别占全下游的66.8%和85.0%。
3.2. 小浪底水库运用初期下游冲淤特征
黄河小浪底水库于1999年10月开始蓄水,由于水库下泄清水使得下游河道全线冲刷。如表1所示:2000~2007年黄河下游(小浪底至利津河段)汛期来水691.5亿m3,仅占总来水量的38.1%;同期下游河道汛期来沙4.97亿t,占全部来沙的93.6%,下游年均来沙量不到1986~1999年平均来沙量的9.0%。河道共冲刷15.55亿t,年均冲刷达1.94亿t,其中汛期冲刷10.35亿t,占总冲刷量的66.6%。期间小浪底水库库区年均淤积泥沙约4.20亿t。
由此可见,黄河下游河道的冲淤与来水来沙条件密切相关。从三门峡水库历年冲淤资料来看,下游淤积主要集中在花园口至艾山河段,然而1999年小浪底水库开始运行后,水库下泄清水,下游河道普遍发生冲刷。胡一三和张晓华 [13] 根据黄河下游多年断面资料研究表明,花园口至高村是二级悬河最为发育的河段,各断面滩槽高差均为负值,大部分达到−1.0 m以上。这与黄河下游泥沙淤积主要发生在花园口至高村河段密切相关。通过以上分析,选择黄河下游花园口至高村河段作为研究的代表河段,研究花园口至高村河段达到冲淤平衡时的临界输沙量显得尤为重要。
4. 临界输沙量研究
以前人 [9] [10] 对黄河下游河道冲淤临界含沙量和流量关系的研究结果为基础,总结出不同流量级的临界含沙量,即流量级为<1000、1000~2000、2000~3000和3000~4000 m3/s的临界含沙量分别为10、20、25和35 kg/m3。以花园口至高村河段为代表,依据孙高虎等 [14] 根据花园口水文站逐年径流量与多年平均值(1961~1999年)比值的变异程度得出的水文年型,即平水年(1976~1978年)、丰水年(1982~1984年)和枯水年(1986~1988年)平均径流量分别为411.1亿m3、525.7亿m3和292.3亿m3,并统计各流量级流量出现频次,见表2。
根据式(1)和表2计算出不同水文年汛期和非汛期达到冲淤平衡的临界含沙量,由于流量大于4000 m3/s以后下游河道易发生冲刷和漫滩现象 [12] ,故忽略此流量级的临界含沙量计算。利用典型代表年汛期和非汛期来水量,求得不同水文年典型代表年份汛期和非汛期临界输沙量(表3)。花园口断面不同水文年典型代表年份汛期和非汛期临界输沙量分别为平水年4.60和1.94亿t、丰水年5.95和2.88亿t、枯水年2.43和1.52亿t,不同水文年临界输沙量分别为平水年6.54亿t、丰水年8.83亿t和枯水年3.95亿t,平均值为6.44亿t。也就是说在平均状态下,进入下游河道泥沙在6.44亿t左右时,下游河道将处于临界的冲淤平衡状态;当进入下游河道泥沙大于6.44亿t时,下游河道将处于淤积状态;相反,当来沙量小于6.44亿t时,下游河道将出现冲刷。
而高村断面不同水文年典型代表年份汛期和非汛期临界输沙量分别为平水年4.25和1.64亿t、丰水年5.88和2.47亿t、枯水年2.00和1.22亿t,不同水文年临界输沙量分别为平水年5.89亿t、丰水年8.35亿t、枯水年
![](Images/Table_Tmp.jpg)
Table 1. Condition of water and sediment with annual distribution of scour and deposition in the lower reaches of the Yellow River
表1. 2000~2007年黄河下游水沙条件及冲淤的年内分布
注:数据引自《小浪底水库拦沙初期运用分析评估报告》(黄河水利委员会,2007年)。“-”表示冲刷。
![](Images/Table_Tmp.jpg)
Table 2. Typical characteristics of annual flood and non-flood seasons daily flow characteristics in Huayuankou Hydrological station
表2. 花园口水文站典型代表年汛期和非汛期逐日流量特征
![](Images/Table_Tmp.jpg)
Table 3. The critical sediment discharge in the balance of erosion and deposition of the Huayuankou and Gaocun section
表3. 花园口和高村断面冲淤平衡时的临界输沙量
3.22亿t,平均值为5.82亿t,相比花园口断面高村断面计算结果偏小0.62亿t,偏小比例为9.6%。这主要是因为洪水的沿程衰减和人类引水引沙有关。黄河下游(花园口至利津)河段是沿黄引水灌溉工程主要河段之一,年均引沙达1.2亿吨 [15] 。
5. 结论与讨论
根据前人研究资料,分析了三门峡水库和小浪底水库运行后黄河下游河道冲淤特征,从三门峡水库历年冲淤资料来看,下游淤积主要集中在花园口至艾山河段,然而1999年小浪底水库运行后下游河道普遍发生冲刷。利用黄河下游花园口和高村水文站逐日水沙资料,计算得到下游花园口至高村河段水沙平衡时的临界输沙量为平水年6.54亿t、丰水年8.83亿t和枯水年3.95亿t,平均值为6.44亿t。胡春宏等 [7] 根据1974~1979年实测水沙资料,得到下游冲淤率为0的临界含沙量为21 kg/m3;龙毓骞等 [8] 根据1950年以来实测水沙资料研究得到黄河下游临界平衡的含沙量也接近21 kg/m3。依据以上研究结果,根据花园口站多年平均(1974~2007年)径流量320.7亿m3,得到下游河道冲淤平衡时的临界输沙量为6.7亿t。许炯心等 [1] 以1962~1985年水沙资料为依据,得出下游河道泥沙不淤的年全沙来沙量应小于7.0亿t。袁吉栋等 [6] 根据1950~2000年黄河下游河道汛期实测径流量和输沙量资料,分析得到汛期花园口和高村河段冲淤平衡的输沙量平均值分别为8.24亿t和7.82亿t,花园口至高村河段平均值为8.03亿t,由于没有考虑非汛期情况使得结果有偏差。因此,本研究得出的6.44亿t作为平衡状态下黄河下游冲淤平衡时的临界输沙量具有一定的合理性。但需要指出的是,本研究没有考虑高含沙洪水、下游不同河段的河床及边界条件的差异、洪水沿程衰减及人类引水引沙的影响,通过高村水文站的计算,所得结果略小于花园口站,偏小9.6%。对于大于4000 m3/s的洪水,可能会形成漫滩,或主槽冲刷滩区淤积现象 [12] ,在计算临界含沙量时没有考虑此流量级流量的影响,可能对结果产生一定的偏差;此外,结果是在全沙条件下得出的,对于不同粒径组泥沙,下游的输沙能力是不同的,尤其是粒径 > 0.025 mm的粗泥沙是黄河下游淤积主体 [1] 。另外,黄河下游水沙条件与黄土高原不同侵蚀类型区产沙量、水土保持综合治理效果及河道泥沙输移比密切相关。因此,冲淤平衡条件下的临界输沙量有待进一步深入研究,以得出更合理的结果。