一种基于灰色系统预测的离散微分动态规划方法
【专利摘要】本发明属于水电优化调度技术领域，公开了一种基于灰色系统预测的离散微分动态规划方法(GDDDP)。该发明的技术方案为：以离散微分动态规划方法(DDDP)为基础，每次迭代前把历史迭代中得到的最优轨迹和预测轨迹组成的数据序列作为灰色系统预测方法的输入进行预测；然后在各个水库预测轨迹的基础上进行离散迭代求解，以此来提高算法的求解精度及计算效率；其中，针对灰色系统预测方法对于振荡序列有较大误差的缺点，采用等差数列递推式改进了灰色系统预测方法。本发明通过改进的灰色系统预测使算法更好地向全局最优解收敛，增加了算法的全局搜索能力；同时通过预测加快了收敛速度，提高了计算效率。
【专利说明】
一种基于灰色系统预测的离散微分动态规划方法
技术领域
[0001] 本发明属于水电优化调度技术领域，涉及到一种基于灰色系统预测的离散微分动 态规划方法(GDDDP)。
【背景技术】
[0002] 水库优化调度模型经过多年的不断发展，已有大量的优化算法相继被提出，这些 优化算法中大多存在"维数灾"问题。离散微分动态规划(DDDP)是一种经典的水库群"降维" 优化算法，但该算法受初始轨迹影响较大。为了提高初始解的质量，路志宏等(路志宏，施保 昌，周晓阳.梯级电站优化调度模型中全局寻优策略[J].人民长江，2007，38(8): 72-74， 84.)把遗传算法求得的近似解作为DDDP的初始解，并对梯级水电站的优化调度模型进行求 解；白小勇等（白小勇，王晨华，李允军，等.人工鱼群算法与离散微分动态规划结合在水库 优化调度中的应用[J].水电自动化与大坝监测，2008，32(6): 66-69.)把人工鱼群算法求得 的解作为DDDP的初始解进行优化调度。冯仲恺等（冯仲恺，廖胜利，牛文静，等.梯级水电站 群中长期优化调度的正交离散微分动态规划方法[J].中国电机工程学报，2015,35(18): 4635-4644.)用正交试验设计选取部分DDDP状态组合对水库中长期发电调度模型进行迭代 求解。然而，这些研究通常只是将两种算法进行简单的串联式耦合，对DDDP算法没有实质性 的改进，甚至牺牲了求解效率;或者为了提高算法的求解效率，而降低最终求解结果的全局 最优性。因此，在以DDDP法为基础的水库群优化调度中，对如何提高解的全局最优性和算法 效率问题进行研究具有重要意义。
[0003] 1982年邓聚龙教授建立的灰色系统理论，是一种研究少数据、贫信息不确定性问 题的新方法(刘思峰，谢乃明.灰色系统理论及其应用[M].第六版.北京:科学出版社，2013: 97-100,153-156.)。经过多年的发展该方法已经在工农业、能源交通等领域得到广泛地应 用。在水库调度中也逐渐有学者尝试把灰色模型与优化调度模型进行耦合，如马志鹏和陈 守伦(马忠鹏，陈守伦.水库预报调度的灰色动态规划模型[J].水力发电学报，2007,26(5): 7-9.)把灰色理论中的灰数理论方法耦合到动态规划算法中对水库进行调度；马志鹏和李 杰等(马志鹏，李杰，董延军，等.灰色动态规划方法在水库调度中的应用研究[J].中国农村 水利水电，2009,（6) :56-58.)利用区间比较的可能度公式及可能度矩阵解决区间灰数的排 序问题并与动态规划算法耦合进行水库调度;马志鹏和袁建国等(马志鹏，袁建国，石赞赞. 基于随机赋权法的梯级水库群灰色决策模型[J].水电能源科学，2009，27(1):77-80.)将灰 色理论中的灰色关联分析方法与随机赋权法结合进行水库调度;武新宇等（武新宇，范祥 莉，程春田，等.基于灰色关联度与理想点法的梯级水电站多目标优化调度方法[J].水利学 报，2012，43(4):422-428.)将灰色理论中的灰色关联分析方法与熵权理想点法相结合对梯 级水电站进行优化调度。但灰色系统预测方法作为灰色系统理论中最活跃的研究领域，却 很少有人将该预测方法应用到优化调度中。本发明利用灰色系统预测方法使用较少数据即 能预测的特点，使用该方法求出水库的预测轨迹，实现灰色预测方法与DDDP算法的结合，用 于提高梯级水库群发电优化调度的求解效率及精度。
【发明内容】
[0004] 针对现有技术的不足，本发明提供一种基于灰色系统预测的离散微分动态规划方 法(GDDDP)。本发明将改进后的灰色系统预测方法与离散微分动态规划方法(DDDP)相结合， 建立基于灰色系统预测的离散微分动态规划算法(GDDDP)。
[0005] 本发明的技术方案为：
[0006] 一种基于灰色系统预测的离散微分动态规划方法(GDDDP)，具体包括以下步骤： [0007]第一步，确定初始计算条件，包括梯级水库群优化调度目标函数、约束条件和决策 变量；
[0008] 第二步，设定状态离散数目K、最大迭代次数c和终止精度ε等计算参数；
[0009] 第三步，使用常规动态规划(DP)算法生成满足第一步中约束条件的各个水库水位 的最优轨迹；
[0010]第四步，初始化迭代次数1=〇,并以第三步的最优轨迹作为此次迭代的试验轨迹；
[0011] 第五步，在试验轨迹每一点的可行范围附近用公式（1)形成一个搜索廊道;在搜索 廊道中，使用常规DP算法寻求最优轨迹;令迭代次数I = 1+1;
[0012] {Zpt+ΔΖ · k},k=[-(K-l)/2,(K-l)/2] (1)
[0013] 其中，Δ Z为离散水位步长;Zpt为初始试验轨迹的水位值。
[0014] 第六步，计算相邻两次最优轨迹各时段的水位差值，若相邻两次最优轨迹各时段 水位差值满足精度要求或迭代次数超过设置的最大迭代次数，则算法搜索终止，输出最优 轨迹;否则转至第五步。
[0015] 其中，当I = 1时，以I = 〇时的最优轨迹为试验轨迹；当I = 2时，以第三步的最优轨 迹和1=〇、1 = 1时的最优轨迹作为原始数据序列，采用改进的灰色系统预测方法，求出各个 水库的预测轨迹，把预测轨迹作为DDDP算法的试验轨迹；当1>2时，把历史迭代求得的最优 轨迹和预测轨迹作为原始数据序列，采用改进的灰色系统预测方法，求出各个水库的预测 轨迹，把预测轨迹作为DDDP算法迭代的试验轨迹。
[0016] 灰色系统预测方法对振荡序列有较大误差，而改进的灰色系统预测方法对振荡序 列误差较小;所述的第六步中改进的灰色系统预测方法，包括以下子步骤：
[0017] 6.1)确定水库在迭代中的最优轨迹和预测轨迹在第h时段末水位值组成的原始数 据序列X (())= lx(())(l)，x(<))(2)，···，x((3)(m)};其中，m代表作为灰色系统预测方法输入序列的 轨迹个数，x ((3)(m)代表第m条轨迹在第h时段末的水位值;所述的h时段指的是调度的各个时 段。
[0018] 6.2)输入原始数据序列乂(())，判断乂(())是否为振荡序列，，得出水库预测轨迹第1 1时 段末的水位值；
[0019] 6.2.1)若Xw不是振荡序列，则根据原始差分GM(1，1)模型对该序列进行预测，得 出各个水库的预测轨迹第h时段末的水位值i(a)饮)；
[0020] 所述的得到原始差分GM(1，1)模型的步骤如下：
[0021] ①对该原始数据序列Xw进行一次累加生成新的序列:Χω = {Χω(1)，χ(υ(2)，…， x(1)(m)}。其中，
[0022]
( 1 )
[0023] ②建立原始差分GM(1，1)方程
[0024] X ⑶（k)+ax ⑴（k)=b (2)
[0025] ③令3 = [α,/?；Γ，运用最小二乘法估计求解待估向量J
[0029 ] ④求得a和b后，代入公式（2)解方程式可求出[0030]
(5)
[0026]
[0027]
[0028]
[0031] 其中，i((M(/c)为x((3)(k)的水位预测值。
[0032] 6.2.2)若乂(())是振荡序列，由于原始差分61(1，1)模型对振荡序列预测的误差较 大，则去掉振荡序列中趋势变化前的旧序列值，保留振荡序列中趋势变化后最新的两个序 列值，对保留序列进行线性延伸预测，得到该水库预测轨迹第h时段末的水位值为
[0033]本发明的有益效果为：本发明以离散微分动态规划DDDP为基础，采用等差数列递 推式改进了灰色系统预测方法并用于水库预测轨迹，在梯级水库群中各个水库预测轨迹的 基础上进行离散迭代求解，有效提高了算法的求解精度及计算效率。
【附图说明】
[0034] 图1是基于灰色系统预测的离散微分动态规划算法计算流程图。
【具体实施方式】
[0035] 下面通过实施例，并结合图1，对本发明做进一步说明。
[0036]白山水库、丰满水库是第二松花江流域两座以发电为主的水利枢纽。白山水库下 距丰满坝址250km，地处吉林省东部山区桦甸与靖宇两县交界处，汛后允许最高蓄水位 416m，正常蓄水位413m，汛期限制水位413m，死水位380m，电站最大引用流量1500m 3/s，保证 出力16 · 7 X 104kW，电站最大出力155 X 104kW。丰满水库是位于第二松花江距吉林市城区东 南24km处，控制流域面积4.25 X 104km2，占第二松花江流域面积的57.9%，为多年调节水库， 汛后允许最高蓄水位263.5m，正常蓄水位26lm，汛期限制水位26 lm，死水位242m，电站最大 引用流量1126 · 5m3/s，保证出力16 · 6 X 104kW，电站最大出力60 · 25 X 104kW。为保证水库安 全，要求两水库在7,8月份水位不超过汛限水位，到9月初以后才允许超蓄，直到汛后允许最 高蓄水位。下面以白山-丰满梯级水电站水库为例，采用1987-1988年(枯水年）的径流资料 分别采用GDDDP与DDDP进行水电优化调度计算。
[0037]步骤1，采用⑶DDP进行对该梯级水库优化调度计算。
[0038] 第一步，确定初始计算条件，包括梯级水库群优化调度目标函数、约束条件、决策 变量；
[0039] 该梯级水库群优化调度是根据入库的流量过程，以发电引用流量为决策变量，以 水库水电站群在调度周期内总发电量最大为目标进行调度。其目标函数及约束条件如下：
[0040] (1)目标函数：
[0041]
[0042] 其中，E为梯级电站总发电量;Qpt为电站p在时段t的发电流量;Mt为时段t的小时 数;R Pt为电站P在一定水头条件下时段t的平均耗水率，单位:立方米/千瓦时;η为水库个数 (这里η = 2) ;Τ为计算总时段(这里η = 12);
[0043] (2)约束条件
[0044] 水量平衡约束 Vp,t+i = Vpt+(Ipt-Qpt-SPt)
[0045] 关联方程 Ip+i,t = Qpt+Spt+IZpt;
[0046] 水库蓄水量约束 Vpt min彡Vpt彡Vpt max
[0047] 水轮机过水流量约束Qpt min彡Qpt彡Qpt max
[0048] 溢洪道下泄水量约束0彡spt彡sptmax
[0049] 边界约束 Vpi = Vpic;Vp,t+i = Vp,t+i,c
[0050] 水库泄流量约束 qPt min彡Qpt+Spt彡qPt max
[0051] 水电站出力约束 NP min彡3600 · Qpt/Rpt彡
[0052] 水位库容关系 Vpt = UP(Zpt)
[0053] 其中，VP,t+1为水库p在t时段末的库容;Vpt为水库p在t时段初的库容；I pt为水库p在 时段t的平均入库流量；IZpt为水库p在时段t的区间入流;Spt为水库p在时段t的弃水流量； Spt max为水库p的溢洪道在时段t的泄流能力；△ t为第t时段长;Vpt min为水库p在t时段初应 保证的最小库容;Vpt max为水库p在t时段初允许的最大库容;Qpt min为电站p在时段t允许的 最小引用流量;Qpt max为电站p在时段t允许的最大引用流量，由于水库群的发电用水与下游 灌概及城市用水相结合，因此在Qpt〈Qpt max时，Spt = 0，而^Qpt = Qpt max时，Spt多0;Vplc，Vp,T+l,c 分别为给定的P水库在优化调度期的初、末库容;qPt min为水库p在时段t的最小下泄流量； qPt max为水库p在时段t的下游河道安全泄流量;Np min为电站p的最小出力;Np max为电站p的 最大出力限制;Zpt为水库p在t时段初的水位;UP为水库p的水位-库容关系函数。
[0054] 本梯级水电站优化调度数学模型，目标函数是求η · T个决策变量Qn，Q12，…，QnT的 总发电量最大问题，由于发电流量Q Pt是水库水位Zpt的隐函数，问题可以转化为求η · T个水 位Zll，Zl2，…，ΖηΤ的总发电量最大问题。
[0055] 第二步，设定计算参数，其中状态离散数目Κ = 7、最大迭代次数c = 150、终止精度ε =0.005m；
[0056] 第三步，使用常规DP算法，对梯级水库水位进行离散，其中白山水库水位离散区间 为380m-413m，丰满水库水位离散区间242m-261m，并生成满足各约束条件的各个水库的初 始试验轨迹；
[0057]第四步，初始化迭代次数1 = 0,并以第三步的结果作为此次迭代的试验轨迹；
[0058]第五步，在试验轨迹每一点的可行范围附近用{Zpt+AZ · k}，k=[_4,4](其中ΔΖ =0.01m)形成一个搜索廊道。在廊道中，使用常规DP算法寻求最优轨迹，完成一次迭代，令I = 1+1;
[0059] 第六步，计算相邻两次最优轨迹各时段的水位差值，若相邻两次最优轨迹各时段 水位差值满足精度要求或迭代次数超过设置的最大迭代次数，则算法搜索终止，输出最优 轨迹;否则转至第五步。
[0060] 其中，当1 = 1，以1 = 0时的最优轨迹为试验轨迹；当I = 2，以第三步的最优轨迹和I =0、1 = 1时的最优轨迹作为原始数据序列，采用改进的灰色系统预测方法，求出各个水库 的预测轨迹，把预测轨迹作为DDDP算法的试验轨迹；当1>2,把历史迭代求得的最优轨迹和 预测轨迹作为原始数据序列，采用改进的灰色系统预测方法，求出各个水库的预测轨迹，把 预测轨迹作为DDDP算法迭代的试验轨迹。
[0061] 改进的灰色系统预测方法详见
【发明内容】
。
[0062] 步骤2,采用DDDP进行对该梯级水库优化调度计算，DDDP为已知算法，不再赘述，其 状态离散数目K、最大迭代次数c、终止精度ε以及离散水位步长△ Z均和⑶DDP相同。
[0063] 步骤3，比较两种算法优化结果。两算法的运行结果见表1，由结果可见，用⑶DDP算 法可以求得比DDDP算法更优的解，发电量增加556 X 104kW · h，同时计算机运行时间还可以 大幅降低，为DDDP的68%左右。该结果表明⑶DDP与DDDP相比，不仅具有更高的全局搜索能 力，而且有更高的计算效率。主要原因是GDDDP算法在每次迭代前采用灰色系统预测方法预 测的轨迹作为本次迭代的试验轨迹。因此，通过灰色预测可以使算法更好地向全局最优解 收敛，增加了算法的全局搜索能力；同时通过预测加快了收敛速度，提高了计算效率。
[0064] 表1 1987-1988年来水条件下的运行结果对比
[0065]
【主权项】
1. 一种基于灰色系统预测的离散微分动态规划方法，其特征在于，包括W下步骤： 第一步，确定初始计算条件，包括梯级水库群优化调度目标函数、约束条件、决策变量； 第二步，设定状态离散数目K、最大迭代次数C和终止精度ε等计算参数； 第Ξ步，使用常规动态规划算法生成满足各约束条件的各个水库的初始试验轨迹； 第四步，初始化迭代次数I = 0，并W第Ξ步的结果作为此次迭代的试验轨迹； 第五步，在试验轨迹每一点的可行范围附近用公式（1)形成一个捜索廊道;在捜索廊道 中，使用常规动态规划算法寻求最优轨迹;令迭代次数I = 1+1; (Zpt+ΔΖ · k}，k=[-化-1)/2,化-1)/2] (1) 其中，A Z为离散水位步长;Zpt为初始试验轨迹的水位值； 第六步，计算相邻两次最优轨迹各时段的水位差值，若相邻两次最优轨迹各时段水位 差值满足精度要求或迭代次数超过设置的最大迭代次数，则算法捜索终止，输出最优轨迹； 否则转至第五步； 其中，当1 = 1时，Wl=0时的最优轨迹为试验轨迹；当I = 2时，W第Ξ步的最优轨迹和I =0、1 = 1时的最优轨迹作为原始数据序列，采用改进的灰色系统预测方法，求出各个水库 的预测轨迹，把预测轨迹作为DDDP算法的试验轨迹；当1〉2时，把历史迭代求得的最优轨迹 和预测轨迹作为原始数据序列，采用改进的灰色系统预测方法，求出各个水库的预测轨迹， 把预测轨迹作为DDDP算法迭代的试验轨迹。2. 根据权利要求1所述的一种基于灰色系统预测的离散微分动态规划方法，其特征在 于，所述的第六步中改进的灰色系统预测方法W下子步骤： 6.1) 确定水库在迭代中的最优轨迹和预测轨迹在第h时段末水位值组成的原始数据序 列xW = {xW(l)，xW(2)，…，其中，m代表作为灰色系统预测方法输入序列的轨迹 个数，xW(m)代表第m条轨迹在第h时段末的水位值;所述的h时段指的是调度的各个时段； 6.2) 输入原始数据序列XW，判断xW是否为振荡序列，得出水库预测轨迹第h时段末的 水位值； 6.2.1)若XW不是振荡序列，则根据原始差分GM(1，1)模型对该序列进行预测，得出各 个水库的预测轨迹第h时段末的水位值公反;I; 所述的得到原始差分GM( 1，1)模型的步骤如下： ① 对该原始数据序列xW进行一次累加生成新的序列:χω = Ιχ^α)，χω(2)，…，xW (m)};其中，(1 ) ② 建立原始差分GM( 1，1)方程 x(〇)(k)+ax(i)(k)=b (2)③令a = [a剧T，运用最小二乘法估计求解待估向量為(3) 其中，(4) ④求得a和b后，代入公式(2)解方程式可求出(5) 其中，批)为X W (k)的水位预测值. 6.2.2)若XW是振荡序列，由于原始差分GM( 1，1)模型对振荡序列预测的误差较大，贝U 去掉振荡序列中趋势变化前的旧序列值，保留振荡序列中趋势变化后最新的两个序列值， 对保留序列进行线性延伸预测，得到该水库预测轨迹第h时段末的水位值为 -(巧? +1) = 2_\：(('')(/巧)-Λ'")'(巧;-1)。
【文档编号】G06Q50/06GK105976058SQ201610290117
【公开日】2016年9月28日
【申请日】2016年5月3日
【发明人】彭勇, 张弛, 周惠成, 史亚军, 李昱, 丁伟
【申请人】大连理工大学