(广东粤电南水发电有限责任公司,广东 韶关 512700)
流域径流不仅是水资源领域的重要组成部分,也是自然界水文循环的重要组成部分。径流资源的时空分布,尤其是年内径流的分布,是水资源评价的重要组成部分。
年内径流量分布不均匀,不仅直接影响水资源开发利用和水利工程规模,而且增加水文预报,水库调度,洪涝灾害控制难度。目前,常用径流非均匀性分布的定量研究指标主要表现为月或季平均年径流量占年径流量的相对比例和汛期、枯水季的相对数来表示。此外,还有连续4个月的最大年径流量占全年径流量的百分数来表示年度径流分布。
南水水库以上流域多年平均降水量2211.8mm。流域内的降水量年内分布不均匀,降雨主要集中在4—9月份,尤其以5月上旬至9月下旬为暴雨爆发的集中期,占全年降雨量的75%左右[1]。
南水水库库区,由于热带低压和南海季风云团趋于活跃,并将逐步加强北抬,对流域内冷空气的强度、途径,以及冷气环流的循环有很大的影响,进而使该流域内产生的降水量级也具有很大的差异性。受到副热带高压的阻挡效应而生成静止切变,此时常常可能会在流域内形成大暴雨。台风也是影响南水水库流域水文气象特征的重要因素,但是由于该水库处于粤北地区,偏内陆,且四周有群山环绕,因此很少有台风中心直接经过水库库区,通常南水水库库区降水是受台风外围波及而造成的。
年径流量的分布决定了水库运行的模式,径流的变化会导致水资源供需的变化,从而影响水资源的开发利用和管理。本文采用径流年内不均匀分配系数、年内完全调节系数等,并利用径流分布的集中度来衡量年内径流量分布的不均匀性;利用汛期的集中度和集中期来揭示汛期径流的集中程度[2]。
定义Cv为径流年内分配的不均匀系数,Cv值越大,表明径流年内分配越不均匀。
其计算公式如下:
(1)
Ri——第i月径流量,108m3;
σ——均方差。
定义Cr为径流年内分配的完全调节系数,Cr值越大,表明径流年内分配越集中。
其计算公式如下:
(2)
(3)
式中, 符号含义同式(1)。
集中期和集中度是利用矢量原理来确定时间分布特征的参数,将该原理应用于南水水库汛期(4月1日至10月31日),可以定量描述汛期的集中度和集中期,对于洪水灾害的发生具有很好的判断作用。该原理是将汛期期间的各旬径流量作为一个向量,将流量的大小视为向量的长度,所处的旬为向量的方向。将汛期的214d作为一个圆周,1d等于1.682°,以每旬的中间日期5日、15日、25日作为这个旬的向量的方向[2]。然后各矢量累加,其和占汛期径流量的百分数即集中度,用RCDx表示,其意义反映在汛期中径流的集中度。如果某年汛期中,该月的径流量完全集中在该月中的某一旬内,则表示集中度达到最大值,即为1,越接近1,表明径流越集中;如果该月每个旬的径流量都相差不大,则表示集中度达到最小值,即为0,越接近0,表明各旬的径流量处于比较平均的水平。对于集中度的计算则是合成该旬代表的向量的方位角,从而确定在全部旬期内降水量在合成后的整体效应,进而确定最大旬径流量出现在汛期的哪个时段[3-4]。
按以上方法合成后径流向量的方向可求得集中期D,即反映径流量最集中的部分。
计算公式如下:
(4)
(5)
(6)
式中Rx——21个旬汛期径流量的分量之和所构成的水平分量;
Ry——21个旬汛期径流量的分量之和所构成的垂直分量;
ri——第i旬的汛期径流量;
θi——第i旬径流代表的向量角度(见表1);
i——旬序号(i=1,2,3,…,21)。
表1 汛期各旬包含的角度及各旬代表的角度值
对于南水水库流域来说,短时间内大量的集中降雨是产生流域洪水的主要原因,也是其主要的库容来源,因此其年内径流量分配过程跟降雨量有一定的相似性,两者都随着季节变化,有着明显的降雨和径流集中现象。南水水库流域1970—2016年各月径流年内变化趋势如图1所示,各月径流在年内分配比例见表2。
图1 南水水库流域内径流量的年内变化
从图1可以看出,各年代径流量的年内分配变化基本保持一致,均呈“单峰型”分布,1月、2月处于低值,3月处于上升期,最大值出现在5月,6月以后进入下降期,冬季(10—12月)达到低值。1970—1999年,径流呈现逐渐减少的趋势;2010年以来,整体偏丰,月平均径流量高于前40年,月平均径流量明显偏大。
图2 南水水库汛期旬内平均入库流量与旬内径流均值量占全汛期总径流量的百分比
由图2可以看出,南水水库的降雨径流时空分布特征呈“三峰”状态:
a.汛期旬平均入库流量由4月上旬开始逐步增大,到5月下旬出现旬平均入库流量最高峰,6月为台风高发期,所以在6月下旬平均入库流量再创新高,7月上旬到10月下旬,入库流量逐步减少,但在8月下旬受降水的影响,入库流量依然还存在一个小峰值。
b.汛期径流中,5月下旬至8月下旬径流比例最大,约占汛期径流总量的60%,三峰结构在旬径流占总径流比例图中表现得很明显,这与当地的气候气象条件有关,尤其是台风降雨对径流的影响极大。
c.汛期中的旬平均入库流量在5—8月比较大且波动也大,约为30m3/s—75.57m3/s,最大的旬平均入库流量出现在5月下旬,为75.57m3/s,6月下旬次之,为72.61m3/s,此后平均入库流量呈下降趋势。
由表2可以看出南水水库流域径流量年内分配不均匀,径流主要集中在汛期(4—9月),汛期径流量占到全年径流量的70%以上,径流量连续最大的2个月(5—6月)就占到全年径流量的1/3以上,是连续最小2个月(12月和次年1月)的7倍左右。
表2 南水水库流域内径流量的年内分配
由此可以看出影响径流年内分配的主要因素是气候变化。因为流域地处亚热带气候湿润区,当副热带高压位于流域北部且位置稳定时,在南方或东南方将有大量的暖湿气流会沿副热带高压东侧北上,此时常常可能会在流域内形成大暴雨,由此造成流域降水年内分配极不均匀,这种气候特点决定了该流域夏秋径流量大而春冬径流量小的特性,即径流年内分配的不均匀性跟该流域气候有关,且气候特点占主要因素[5]。
表3为南水水库流域的Cv(径流年内分配不均匀系数)和Cr(完全调节系数)的计算成果。由此可知,该流域的Cv极小值是0.71,极大值是1.80,均值是1.25。Cv∈(0.2,0.8),即径流年内分配完全调节系数一般都处于0.2~0.8之间,整体水平较高,说明汛期径流年内分配较不均匀,因为该流域的补蓄水以降雨为主,降雨量随季节存在周期性变化,这就会导致径流量同样存在较为明显的季节性差异变化[6]。
表3 南水水库流域Cv(径流年内分配的不均匀度)和Cr(完全调节系数)
从图3可以看出,1970—2016年南水水库流域Cv总体上呈锯齿形波动,20世纪70、80年代的Cv比较低,主要是由于20世纪70、80年代径流偏枯;90年代以后,Cv呈缓慢增加的态势;进入21世纪后,Cv呈缓慢降低态势,且变幅增大,这表明南水水库流域年内径流量的波动程度在逐渐增大,其径流年内分配正在趋向不均匀,对水资源的充分利用处于负面影响状态。
图3 南水水库流域Cv的变化趋势
从汛期径流分配的集中度(见表4)来看,1970—2016年集中度多在30%~80%之间,最低出现在1980年、1987年和2009年,都是17%;最高则为1974年、2002年和2013年,都为79%;多年平均集中度比较高,为49.5%。对集中度的变化趋势进行分析,发现集中度先趋于降低后趋于增大(见图4),这说明在这47年中,年内径流分配集中度由高到低再到高,年内分配经历了先从不均匀趋于较均匀后又趋于不均匀。
汛期的集中度≥70%的有5年,其中1974年、2002年和2013年属于偏丰年,1980年、1987年和2009年属于偏枯年, 就集中期而言,汛期的径流量集中于旬序号为5、6、10、14、15,分别是5月1—21日、7月1—10日、8月11—31日;而目前的南水水库汛期为4月15日至9月15日,计算结果表明,汛期径流分配的集中期和集中度越大,表明该流域发生洪水的概率也就越高;汛期径流分配集中度呈增大趋势,表明汛期径流分配逐渐趋于不均匀。
表4 南水水库流域汛期径流分配的集中度和集中期
图4 南水水库汛期集中度变化趋势
以珠江流域北江水系的南水水库的1970—2016年的径流量资料为基础,在多年开展汛期水情趋势预报的基础上,利用南水水库流域47年的水库径流量资料,分析了南水流域的月径流、旬径流变化特征和径流的年内分配不均匀性以及汛期径流集中性等,研究得出南水水库流域年内,年际的径流变化一般特性,对南水水库汛期降水量惊醒了预测和分析,可得以下结论:
径流量年内分配的变化基本保持一致,分布形态均显示为“三峰型”,月径流一般会在5月达到极大值。径流主要集中在汛期,占到全年径流量的70%以上。径流年内分配很不均匀,径流年内分配不均匀系数处于逐渐增加的状态,且变幅也在逐渐加大,这表明南水水库流域年内径流量的波动程度在逐渐增大,汛期的集中度偏大,发生洪水的概率就大,因此需要高度重视未来汛期集中度较高的年份。
基于集中度的南水水库径流年内分配研究作用在于计算相对简单,精度较好,有利于判断汛期大小年,提前做好大汛大旱等极端年份的防汛抗旱工作,有利于领导者在发电效益的最优化分配过程中的决策判断,对水资源的充分利用和发电调度等提供数据性支持。建议在水库调度过程中抓住一个或两个大的降水提高水库储存量,以确保水库蓄水量可以满足向韶关供水等各方面的需要,实现水资源的合理利用。
[1]焦晓东,邵广哲.南水水库年末消落水位的确定[J].广东水利水电,2017(1):26-29.
[2]刘贤赵,李嘉竹,宿庆,等.基于集中度与集中期的径流年内分配研究[J].地理科学,2007,27(6):791-795.
[3]姜洋.基于集中度的大伙房水库汛期分期探讨[J].水利水电技术,2016,47(5):119-123.
[4]蒋丽娟,施能,周文志,等.桂林市汛期降水集中度和集中期与旱涝关系研究[J].气象研究与应用,2007,28(3):23-26.
[5]吴新新.基于改进集对分析的新疆黄水沟年径流预测[J].水资源开发与管理,2016(1):65-67.
[6]丁瑞,焦晓东.基于模糊集合分析法的南水水库汛期气象水文特征分析[J].人民珠江,2017(5):5-8.