1. 引言
随着国民经济和社会发展水平的快速提高,能源需求量不断攀升 [1] [2] 。传统化石能源不仅总量有限,而且还会对环境造成有害影响,开发利用可循环再生的清洁能源符合社会发展的趋势。截至2017年底,我国已建各类水库的总库容为9035亿m3,年末全国水电装机容量34,168万kW,全年发电量11967亿kW·h [3] ,水电已经成为我国重要的清洁能源。
梯级水电站单一调度存在水能利用率不高、补偿效益较弱等典型问题,而联合调度从整体效益出发,对梯级水库统一调度和管理、提高综合利用效率具有重要意义。同时梯级水电站联合调度具有大规模、多维度、非线性等特征,约束条件和水力耦合联系复杂 [4] ,高效求解一直是研究的热点。传统的动态规划算法能遍历整个可行域,但随着水库群规模的增大,会出现“维数灾”问题,直接影响调度模型的有效求解。为此,人们依据动态规划算法进行改进,产生了多种优化算法。逐次优化算法(POA)具有计算效率高、耗费时间短、收敛性好等优点,逐次渐进动态规划算法(DPSA)将多维问题转为多个一维问题,使原问题得到简化。
近年来,众多学者对上述算法在水库调度中的应用做了大量研究。陈炯宏 [5] 利用POA算法求解清江梯级水库汛限水位动态控制模型,显著增加了汛期发电量。为避免POA求解时造成的“维数灾”现象,冯仲恺 [6] 提出了正交逐步优化算法,较于传统算法,求解时间大幅降低。张诚 [7] 基于逐步优化算法的寻优机制,采用逐步差分和变阶段优化的改进策略,算法具有较好的全局搜索能力。赵志鹏 [8] 等利用离散梯度下降法对逐步优化算法子优化问题求解,该方法能够确定最优搜索方向,优化计算规模。Zhang [9] 等将POA算法和DPSA算法结合起来,并且将黄金分割法应用到单纯形法中,能有效计算水库蓄水轨迹。
本文以雅砻江流域为研究对象,以“两河口~锦屏一级~二滩”梯级水电站为调控主体,根据雅砻江梯级水电站水库群实际情况,建立梯级水库联合调度模型,以DPSA-POA混合算法求解,探求雅砻江梯级水库最大发电效益潜力。
2. 研究区域及数据资料
2.1. 雅砻江流域概况
雅砻江是金沙江第一大支流,发源于青海省,经甘孜、新龙、雅江、木里、冕宁、盐边等县,在攀枝花市汇入金沙江。雅砻江干流全长约1571 km,天然落差3830 m,年径流量约600亿m3。流域呈北西—南东走向,平均宽度128 km,略呈“柳叶”状,面积达13.6万km2。
雅砻江流域水能资源丰富且非常集中,干流技术可开发量约3000万kW,占到四川全省的24.93%,调节补偿性能和技术经济指标优良,开发成本低于其他流域,是我国优质的水电能源基地 [10] 。其中两河口是具有多年调节能力的龙头水库,锦屏一级具有年调节能力,二滩水库为季调节水库。
2.2. 梯级水电站资料
两河口水电站,位于四川省甘孜州雅江县境内,为雅砻江中下游梯级电站的控制性水库工程,主要功能为发电,同时具有蓄水蓄能、改善航道枯水期航运条件的作用。锦屏一级水库水电站位于四川省盐源县和术里县境内,雅砻江下游大河湾处,距河口358 km主要开发任务是发电,具有良好的调节性能,是雅砻江下游的控制性工程,在雅砻江开发中起着承上启下的关键作用。二滩水库电站位于攀枝花市,距与金沙江的交汇口33 km,是一座具有季调节能力,以发电为主的大型水力发电枢纽。
Table 1. The main parameters and characteristics water level of cascade reservoirs
表1. 梯级电站主要参数和特征水位
2.3. 径流资料
雅砻江流域的径流资料中,两河口、锦屏一级和二滩水库采用1956年6月~2013年5月的旬入库流量资料(见表1)。下一级水库的入库径流由上一级水库的出库流量和同期两水库之间的区间入流组成。
3. 联合水库优化调度
3.1. 优化调度模型的建立
1) 目标函数
在满足各约束条件下,以梯级发电量最大为目标函数:
,
(1)
式中:T为总时数;n表示电站个数;
为第t时段梯级电站的出力;
为时段长;
为第i个电站第t时段的出力。
2) 约束条件
①水量平衡方程:
(2)
式中:
和
分别表示第i个电站第t时段初、末水库的蓄水容积;
表示第i个电站第t时段平均入库流量;
表示第i个电站第t时段平均出库流量。
②水库的蓄水限制:
(3)
式中:
表示第i个电站第t时段的最小蓄水,一般为死库容;
表示第i个电站第t时段的最大蓄水容积,在汛期为汛限水位对应的库容,在非汛期为正常高水位对应的库容。
③水库出库流量限制:
(4)
式中:
表示第i个电站第t时段电站放水量下限,
表示第i个电站第t时段电站放水量上限。
④水库之间的水量平衡方程:
(5)
式中:
表示第t时段第i个水库至第(
)个水库区间流量。
⑤电站出力限制:
(6)
式中:
、
分别表示第i个电站第t时段的出力下限和上限。
⑥水库边界条件:
(7)
式中:
表示第i个调度期初水库蓄水位;
表示第i个水库调度期末水库蓄水位。本研究均取为死水位。
3.2. 模型求解方法
针对该梯级电站联合调度问题,采用逐次逼近动态规划—逐次优化算法(DPSA-POA)求解。逐次逼近动态规划算法(DPSA)将原来n维状态变量的动态规划问题分解成只有一个状态变量的n个子问题 [11] ,并对n个子问题逐步求解,使这些子问题的优化序列收敛于原问题的解,可大大降低问题的维数,也可减少所需的存储空间。
逐次优化算法(POA)将一个多阶段的决策问题分解为多个两阶段问题,使原问题得到简化 [12] 。固定其他阶段的变量,只对所选两阶段的决策变量寻优,在解决该两阶段问题后再考虑下一两阶段问题,有效地降低了问题的维度。
如图1所示,DPSA- POA算法的步骤为:
1) 输入基本资料。
2) 采用单库最优的方法计算初始状态变量和决策变量,拟定初始调度轨迹。
3) 先对水库1进行优化,其余n-1个水库的运行决策和状态变量保持不变。利用一维优化算法POA对水库1优化调度进行求解,从而可得水库1的运行策略。
4) 按照水库顺序依次对后续水库进行优化,优化第n个水库时,第1~n-1个水库保持更新后的策略,其余水库仍然保持初始策略,同样利用POA算法求解第n个水库的最有运行策略。
5) 遍历所有水库得到所求的优化调度线,并以该调度线为初始调度线,重复上述步骤不断寻优,若优化调度线与初始调度线差异满足精度要求则输出最优调度轨迹;否则,重复步骤(3)~(4),直到满足精度要求为止。
4. 联合调度结果分析
采用DPSA-POA算法进行优化计算,得到梯级联合优化调度结果如表2所示,可知:
1) 通过梯级补偿,各电站不仅可满足各自的其他用水要求,而且梯级发电保证率、发电量均有较大程度提高,其中梯级年均发电量510.77亿kW·h,发电保证率达到97.61%。
2) 各年发电量相差较大,最大年份(2012)达到598.88亿kW·h,超过均值约17.25%;最小年份(1956)仅341.91亿kW·h,仅为均值的57.09%。雅砻江梯级年发电量~保证率关系如图2所示,可知在理想情况下,雅砻江梯级年发电量超过469.77亿kW·h的经验频率约为80%。
3) 各年的发电保证率也有一定差异,最小的年份为1956年,发电保证率为77.78%。雅砻江梯级出力~保证率关系如图2所示,可知通过优化调度,大概有97.56%的时段满足保证出力要求,有5%的时段保证满发。
Figure 1. Flow chart of optimal reservoir operation based on DPSA-POA
图1. DPSA-POA法求解水库优化调度流程图
Figure 2. Relationship between optimized power generation and output with guarantee probability
图2. 梯级电站优化的年发电量和出力与保证率关系图
对雅砻江梯级水库1956~2012年长系列历史资料进行优化调度,得到理论优化值如表3所示。保证出力取3136 MW时,梯级电站年均发电量为510.77亿kW·h,较常规调度方案可年均增加发电量44.03亿kW·h,提高发电量约9.43%;较常规调度减少弃水量57.49%。从对社会效益的贡献看,所增发的电量相当于节约标准煤178万吨,减少CO2排放量440万吨,对国家节能减排政策以及环境保护具有重要意义。可见,雅砻江梯级水库群调度存在优化的空间,得到的最优调度轨迹过程线可作为分析水库调度规则的基础。
Table 2. Joint optimization results of cascade power generation
表2. 梯级电站联合优化调度历年结果
Table 3. Results of cascade hydropower joint optimization scheduling (1956~2012)
表3. 雅砻江梯级水库电站优化调度结果(1956~2012年)
5. 结论和展望
以雅砻江流域两河口,锦屏一级和二滩水库为研究对象,建立联合优化调度模型,采用DPSA-POA法求解,得出如下结论:
1) 采用1956~2012年的资料进行分析计算,梯级水电站联合优化调度的年均发电量达到510.77亿kW·h,较常规调度增加44.03亿kW·h (9.43%),减少弃水量65.91亿m3 (57.49%)。
2) 结合雅砻江梯级水库实际情况,对梯级水库优化调度研究的重要性进行了探讨,优化调度能显著提高水能综合利用效率。下一步将考虑梯级调度生态等问题,深入开展梯级水库长系列多目标优化调度研究。
基金项目
国家重点研发计划(2016YFC0402206),国家自然科学基金(51879192)资助项目。