厂内经济运行动态规划.ppt

水电站厂内经济运行的动态规划方法 水电站厂内经济运行的动态规划方法 阶段变量k——以机组编号为序，机组台数k为阶段变量； 状态变量 ——k台机组的总负荷； 决策变量 ——第k台机组的负荷； 状态转移方程： 最优值函数 ——k台机组共同带负荷 时的最优流量； 递推方程： 初始(边界)条件： 动态规划递推示意图 1 k n 例：已知某水电站4台机组的流量特性如下表所示，求电站的综合流量特性。 k 1 2 3 4 N 0 0 0 0 0 1 30 31 29 28 2 34 34 35 35 3 38 37 38 37 4 40 k 1 2 3 4 J 0 0 0 0 0 0 0 0 0 1 1 30 0 30 1 29 1 28 2 2 34 0\2 34 0 34 0 34 3 3 38 3 37 0 37 0\3 37 4 4 40 0 40 0 40 5 4 70 1 69 1 68 6 4 74 0 74 0 74 7 4 78 0\3 78 3 77 k 1 2 3 4 J 8 1 107 1\3 106 9 3 112 3 111 10 3 116 3 115 11 1\3 144 12 3 149 13 3 153 1 2 3 4 1 2 3 4 0 0 0 0 0 0 7 0 4 0 3 77 1 0 0 0 1 28 8 3\0 4 0\1 1\3 106 2 0\2 2\0 0 0 34 9 2 4 0 3 111 3 0 3\0 0 0\3 37 10 3\0 4 0\3 3 115 4 0 4 0 0 40 11 3 4 1\3 3\1 144 5 0 4 0 1 68 12 2 4 3 3 149 6 2 4 0 0 74 13 3 4 3 3 153 动态规划方法计算程序框图 变量说明 k：机组台号及台数; N(k)：第k号机组的装机容量; Nk：1~k号机组的总容量; Nk1： 1~k-1号机组的总容量; j：状态变量，总出力的取值; i：决策变量， k号机组的出力; Q2(j)：1~k号机组的总工作流量（最优）; Q1(j)：1~k-1号机组的总工作流量（最优）; N0(j,k)：存放优化负荷分配结果，当负荷j由k台机组承担时，k号机组所承担的负荷（最优决策）。 动态规划计算框图 变量初始化，数据输入， k=1 k=1? Nk=N(k) J=0 Q1(J)=Q(k,J), N0(J,k)=J J=J+1 JN(k)? Nk1=Nk,Nk=Nk+N(k) J=0 Y=106,i=0 J-iNk1? iJ? X=Q(k,i)+Q1(J-i) XY? Y=X,I0=i i=i+1 iN(k)? Q2(J)=Y,N0(J,k)=I0 J=J+1 JNk? Q1(J)=Q2(J), Q2(J)=0,J=0~Nk k=k+1 kn? y y n n n y n n y y n y n y n y (1) 负荷分配表 (1) J=Nk M=J k=n N0(J,k)=N0(M,k) M=M-N0(M,k) k=k-1 k1? J=J-1 J0? stop y n n y 动态规划方法的特点 直接利用流量特性，对流量特性没有特殊要求（凸、可微） 可同时确定开机台数、机组组合和负荷最优分配 能保证所得结果是最优解 最优化原理的一个推广定理 定理 设 为电站负荷为N在n台机组之间的最优分配， 是任一机组组合状态向量，其中 则 是对应于这一机组组合状态，负荷为Na的最优分配。 设第l,m号机组停机，若电站负荷为J0，现在的问题是如何把负荷在不包含第l,m号机的机组之间进行分配。步骤： (1)选定J=J0行； (2)计算 ； (3)若J1=J0 ，则得最优分配； 若J1J0 ，则J=J+1 ，J2=J1 ，转(2)； 若J1J0 ，则用线性插值进行负荷分配，即 汽蚀振动区的处理方法 1.罚函数法 为避免机组在汽蚀和振动区[ N1 , N2 ] 运行,可用惩罚函数对机组流量曲线进行修正Qk ( Nk) = Qk ( Nk) +αA , 式中, A 为惩罚量,α= 1 时在汽蚀和振动区,α= 0 为其他工况 2.虚拟机组法 机组流量特性曲线中可运行区间[ N min , N1 ] 和[ N2 ,Nmax] 的部分分别当成虚拟机组K1和K2的流量特性曲线。 对有多个汽蚀振动区的情况也可类似处理 时间最优问题模型 若电力系统在某时段内给电厂的发电任务已确定,于是电厂的运行调度应以耗水量最