算法调论文、.docVIP

计科1303班 徐莹 201309010322 变状态空间动态规划法在水库优化调度中的应用 摘要: 通过分析三峡梯级水电系统日优化调度的数学模型,提出了一种变状态空间动态规划法,并给出了计算实现步骤。该方法有效地克服了常规动态规划法中存在的“维数灾”问题,从而为解决此类问题提供了一种现实可行的途径。针对三峡梯级水电系统进行仿真计算,结果表明方法是有效的。 关键词: 优化调度; 变状态空间; 动态规划法 三峡日发电计划不仅是三峡梯级水电系统优化调度中的一个重要课题,而且也为三峡电力生产部门形成决策提供有效的技术支持,对其进行研究具有十分重要的理论与实用价值。由于电网负荷瞬息万变,同时,负荷预测不可能完全准确, 加之梯级系统中发电设备的事故与故障,使得三峡梯级电站各自的负荷图(或整个梯级电站负荷图)经常发生变化,于是形成了在负荷预测基础上的所谓负荷带宽。因此,制定日发电计划必须满足电力系统的负荷约束并符合梯级电站的实际情况。另外,梯级水电站存在水力联系与电力联系, 优化调度又必须考虑水力平衡约束、离散运行区间和离散状态约束以及非线性水头影响等诸因素,而使用线性规划、非线性规划、网络技术和大系统分解协调等算法求解比较困难,因此,采用动态规划法较为合适。通过分析日发电计划 模型,本文提出了一种变状态空间动态规划法,有效地克服了常规动态规划法中存在的计算“维数 灾”问题。 1 数学模型 1. 1 优化准则与目标函数 优化准则: 给定三峡电站日用水量,在满足水电机组运行约束和电力系统负荷约束等前提下, 如何选择用水(或出力)策略,使梯级系统在一个日调度周期内总电能达到最大。根据优化准则和最优性原理,目标函数为: (1)式中,表示梯级水库从时段初水库蓄水量出发,从时段1到时段t 时的最优发电量; V(s,t)、V(g ,t)表示t时三峡与葛洲坝水库各自的蓄水量; Qs, t、Qg, t表示t 时三峡与葛洲坝水库各自的发电流量; T 为计算时段。 1. 2 约束条件 系统负荷约束 式中, Pi , t、Pi , t分别表示t 时i 水库的系统负荷下限与上限。 出力约束 式中,Ni , t、Ni, t分别表示t 时i 水库的出力下限与上限。 发电流量约束 式中,Qi、Qi 分别表示i水库的发电流量下限与上限。 库容约束 式中,Vi、Vi 分别表示i水库的最小与最大蓄水库容。 水量平衡方程 式中, T为单位转换系数; Q入表示入库径流量; Q弃表示弃水流量;f表示时滞。三峡电站是一个 具有年调节能力的大型水库,而葛洲坝属于径流式水库,因此,优化调度暂不考虑三峡弃水流量, 但要考虑葛洲坝存在的可能弃水量,这一点在式( 6)、式( 7)中反映出来。 1. 3 初始条件 已知两水库期初水位Zs,0、Zg,0 , 梯级电站提供的可投入运行的机组型号和台数,以及三峡电 站日用水量Ws。 2 变状态空间动态规划法及实现步骤应用常规动态规划法求解水库优化调度问题时,状态空间范围的设置在整个调度期内一般是固定不变的,这种处理方式对于梯级水库优化调度的求解会导致计算困难,即出现所谓“维数灾” 问题。为此,针对梯级水库日发电计划模型,本文提出了一种变状态空间动态规划法,有效地解决了上述问题。具体思路是: 由梯级系统负荷约束式(2)通过厂内经济运行方式和水头迭代得到梯级两水库发电流量廊道范围,继而分别确定三峡与葛洲坝的状态空间范围(即水库库容廊道) ,最后在库容廊道范围内用二维动态规划法进行求解。整个优化算法包括如下几个模块: 2. 1 数据准备模块 数据准备模块包括: 基本参数数据的录入; 利用线性插值得到两水库的库容特性和下游水位流 量特性的子程序; 三峡梯级电站厂内最优运行子程序,该程序是整个梯级最优运行方式的基础。 2. 2 预评估和校核模块 根据三峡梯级提供可投入运行的机组台数和预设工作水头,调用厂内最优运行子程序,对机组是否能承担给定的系统负荷带宽进行评估并反馈; 对初始条件如初始库容、用水量Ws 等的设置是否合理进行检验。该模块主要对调度决策的可行性作简单的推理和判断。 2. 3 优化调度启动模块 该模块主要确定变状态空间— 库容廊道范围。三峡库容廊道范围的确定: 由初始水位Zs,0以及三峡电站在t 时段的系统负荷约束上限值Ps, t 与下限值Ps , t ,可以求出t时段的发电流量上限值Qs, t与下限值Qs , t ,然后再以发电流量的上限(下限)为基准逐时段构造出库容廊道的下限(上限) , 在时段t= 95时收敛于期末库容,这样就构成了所谓三峡水库库容廊道。显而易见,以这种方式构造的库容廊道范围随时段变化而发生改变,其范围远比正常蓄水位相应的库容和死库容之间的范围要窄得多,却包含了比较全面且合