考虑历史轨迹知识的梯级水库群生态调度优化方法技术

本发明专利技术涉及梯级水库群生态调度技术领域，公开了一种考虑历史轨迹知识的梯级水库群生态调度优化方法，在给定初始调度过程后，以离散微分动态规划方法(DDDP)为基础框架，首先在每次迭代前，通过挖掘历史轨迹过程携带的知识信息，对当前轨迹进行变尺度重心反向学习，得到结果更好的优化轨迹；然后在优化轨迹的基础上进行离散迭代求解；最后通过迭代寻优逐次逼近全局最优解。本发明专利技术考虑历史轨迹知识的梯级水库群生态调度优化方法，充分利用了历史轨迹知识信息，通过变尺度重心反向学习来挖掘生成多样性的搜索信息，提高了全局搜索能力，有效克服了传统算法的局部收敛问题。

【技术实现步骤摘要】
考虑历史轨迹知识的梯级水库群生态调度优化方法
本专利技术涉及梯级水库群生态调度
，具体涉及一种考虑历史轨迹知识的梯级水库群生态调度优化方法。
技术介绍
伴随着我国水利事业飞速迅猛发展，各个流域大规模梯级水库群不断建成运营，发挥了巨大的防洪、发电、航运等社会效益和经济效益，但也急剧改变了原始的天然径流水文情势，对河流生态系统的结构和工程也产生了诸多不利的影响，包括河道萎缩、泥沙淤积、生物多样性减少。为减缓水库运行对河道生态环境的负面影响，水库群生态调度作为一种典型的非工程措施，通过调整水库群的运行方式，获取生态系统稳定及保持物种多样性最为合适的径流过程。可见，开展水库群生态调度非常有必要，受到国内外的广泛关注。梯级水库群生态调度问题具有大规模、高维度、多阶段、强约束、非线性等特点，传统的动态规划方法(DP)应用于求解梯级水库群联合优化调度问题时，无法避免地面临严重的维数灾问题。离散微分动态规划方法(DDDP)是一种以逐次渐进逼近理论为核心的动态规划改进方法。其基本思路是：首先，根据经验或其他方法获得满足各项复杂约束条件的初始轨迹；然后，在该轨迹的邻域内对各电站不同时段的状态变量进行离散形成廊道；其次，基于常规动态规划方法在各时段离散状态组合间进行优化，寻找出新的最优轨迹作为下次迭代的试验轨迹，反复迭代直至满足收敛条件。但DDDP在处理大规模水库群优化调度问题时，仍然存在维数灾、早熟收敛、搜索能力不足等缺陷。因此，亟待对DDDP的计算机理进行有效改进，提高水库群联合调度问题的计算效率和求解精度。重心反向学习是近年来依托反向学习提出的新智能计算技术，其基本思想是充分利用群体的搜索信息，评估当前状态及其重心反向状态，择优使用，从而加速搜索进程。由于DDDP是逐阶段计算过程，没有群体的概念，本次将各阶段的状态视为一个个体，这样若干个不同阶段的状态则可视为一个群体。为了有效利用历史轨迹信息，本次提出了考虑不同历史轨迹的变尺度重心反向学习，即选取不同尺度的阶段个数，然后组成一个群体并计算重心反向状态，来挖掘生成多样性的搜索信息，提高全局搜索能力。本专利技术成果考虑历史轨迹所携带的知识信息，采用变尺度重心反向学习与DDDP相结合，通过变尺度重心反向学习寻找优化轨迹，来有效克服DDDP的局部收敛问题，提高方法的求解效率和计算精度，对大规模水库群优化调度问题具有良好的支撑应用价值。
技术实现思路
本专利技术的目的就是针对上述技术的不足，提供一种考虑历史轨迹知识的梯级水库群生态调度优化方法，解决DDDP在复杂水电系统优化调度求解时存在的局部收敛问题。为实现上述目的，本专利技术所设计的考虑历史轨迹知识的梯级水库群生态调度优化方法，包括如下步骤：1)确定初始计算条件，包括梯级水库群生态调度的优化目标函数、约束条件和决策变量；2)设定计算参数，包括最大迭代次数M、最大变尺度重心反向学习次数K、各水库初始离散步长、收敛精度ε、水库数目N和调度期时段数目T；3)采用常规动态规划方法生成满足各约束条件的各个水库的初始轨迹Ω1，并将其存入历史轨迹库4)初始化迭代次数m＝1；5)当历史轨迹库数目达到指定数量Y时，历史轨迹库实施考虑历史轨迹知识的变尺度重心反向学习，否则直接转步骤6)，变尺度重心反向学习包括如下步骤：①初始化变尺度重心反向次数k＝1；②随机生成在1到Y之间的正整数y，有1≤y≤Y；③从历史轨迹库中，选取从y次开始依次存入的Y+1-y条历史轨迹{Ωy,Ωy+1,Ωy+2,...,ΩY}，计算被选取所有历史轨迹的各个水库在各个时段的状态平均值其中，为第y条历史轨迹的水库i在时段j的水位值；④由重心反向原理，计算当前轨迹的重心反向轨迹⑤比较重心反向轨迹和当前轨迹，如果优于则以重心反向轨迹替换当前轨迹，否则不作处理；⑥令k＝k+1，如果k≤K，转步骤②，否则转步骤(6)；6)实施常规DDDP算法，在当前轨迹的可行范围内形成一个搜索廊道，在搜索廊道内使用常规动态规划方法寻求当前最优轨迹，并将当前最优轨迹按递进顺序，存入并扩充历史轨迹库7)算相邻两次最优轨迹各时段的水位差值，若则收缩所有水库离散步长，转至步骤8)，否则转至步骤5)；8)令m＝m+1，若m＞M，则转至步骤9)，否则转至步骤5)；9)停止计算，输出最终的最优轨迹。本专利技术与现有技术相比，具有以下优点：1、以DDDP为基础框架，通过挖掘历史轨迹过程携带的知识信息，对当前轨迹进行变尺度重心反向学习，得到结果更好的优化轨迹；2、充分利用了历史轨迹知识信息，通过变尺度重心反向学习来挖掘生成多样性的搜索信息，提高了全局搜索能力，有效克服了传统算法的局部收敛问题。附图说明图1为本专利技术考虑历史轨迹知识的梯级水库群生态调度优化方法的流程图。具体实施方式下面结合附图和具体实施例对本专利技术作进一步的详细说明。如图1所示，如图所示，一种考虑历史轨迹知识的梯级水库群生态调度优化方法，包括如下步骤：1)确定初始计算条件，包括梯级水库群生态调度的优化目标函数、约束条件和决策变量，其中，梯级水库群生态调度模型可以描述为：已知各水库调度期初始水位、末水位和入库径流过程以及区间径流过程，在满足各水库相应的水位、流量等复杂约束的情况下，以梯级水库群的生态溢水量和生态缺水量之和最小为优化目标，梯级水库群生态调度目标函数的数学表达式为：minf＝minVeco＝min(VecoOVER+VecoLACK)式中：Veco为梯级水库群总生态溢缺水量(m3)，VecoOVER为梯级水库群总生态溢水量(m3)，VecoLACK为梯级水库群总生态缺水量(m3)；Oi,j为水库i在时段j的出库流量(m3/s)，为水库i在时段j的适宜生态流量上限值(m3/s)，为水库i在时段j的适宜生态流量下限值(m3/s)，N为水库数目，i为水库序号，且i＝1,2,…,N，T是调度期时段数目，j为时段序号，且j＝1,2,…,T，Δj为时段长度(h)。需要满足的约束条件，主要包括：(1)水量平衡约束：式中，Vi,j为水库i在j时段末的库容值(m3)，Ii,j为水库i在j时段的入库流量(m3/s)；(2)水力联系约束：式中，Ri,j为水库i在j时段的区间流量(m3/s)，Ui为水库i的上游水库数目，Ou,j为上游水库i在j时段的出库流量(m3/s)；(3)时段水位约束：式中，Zi,j为水库i在j时段的坝前水位值(m)，和分别为水库i在j时段的坝前水位最低值和坝前水位最高值(m)；(4)出库流量约束：式中，和分别为水库i在j时段的最小出库流量(m3/s)和最大出库流量(m3/s)；(5)初末水位约束：Zi,0＝Zi,start,Zi,T＝Zi,end式中，本文档来自技高网...

【技术保护点】
1.一种考虑历史轨迹知识的梯级水库群生态调度优化方法，其特征在于：包括如下步骤：/n1)确定初始计算条件，包括梯级水库群生态调度的优化目标函数、约束条件和决策变量；/n2)设定计算参数，包括最大迭代次数M、最大变尺度重心反向学习次数K、各水库初始离散步长、收敛精度ε、水库数目N和调度期时段数目T；/n3)采用常规动态规划方法生成满足各约束条件的各个水库的初始轨迹Ω

【技术特征摘要】
1.一种考虑历史轨迹知识的梯级水库群生态调度优化方法，其特征在于：包括如下步骤：
1)确定初始计算条件，包括梯级水库群生态调度的优化目标函数、约束条件和决策变量；
2)设定计算参数，包括最大迭代次数M、最大变尺度重心反向学习次数K、各水库初始离散步长、收敛精度ε、水库数目N和调度期时段数目T；
3)采用常规动态规划方法生成满足各约束条件的各个水库的初始轨迹Ω1，并将其存入历史轨迹库
4)初始化迭代次数m＝1；
5)当历史轨迹库数目达到指定数量Y时，历史轨迹库实施考虑历史轨迹知识的变尺度重心反向学习，否则直接转步骤6)，变尺度重心反向学习包括如下步骤：
①初始化变尺度重心反向次数k＝1；
②随机生成在1到Y之间的正整数y，有1≤y≤Y；
③从历史轨迹库中，选取从y次开始依次存入的Y+1-y条历史轨...

【专利技术属性】
技术研发人员：邹强，鲁军，王学敏，洪兴骏，李肖男，严凌志，
申请(专利权)人：长江勘测规划设计研究有限责任公司，
类型：发明
国别省市：湖北;42