CN105631528B - 一种基于nsga-ii和近似动态规划的多目标动态最优潮流求解方法 - Google Patents
一种基于nsga-ii和近似动态规划的多目标动态最优潮流求解方法 Download PDFInfo
- Publication number
- CN105631528B CN105631528B CN201510605157.1A CN201510605157A CN105631528B CN 105631528 B CN105631528 B CN 105631528B CN 201510605157 A CN201510605157 A CN 201510605157A CN 105631528 B CN105631528 B CN 105631528B
- Authority
- CN
- China
- Prior art keywords
- dynamic
- dynamic programming
- power flow
- target
- solution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及基于NSGA‑II和近似动态规划的多目标动态最优潮流求解方法,属于电力***日前调度计划研究领域,包括以下步骤:获取下一调度周期内***相关数据;以发电耗费、污染物排放量与网损最小为目标,以节点安全电压和无功出力为约束,建立多目标动态最优潮流模型;采用NSGA‑II求解各单时段Pareto解集,将单时段连续决策空间离散化,以非支配解的模糊满意度作为值函数,基于时段耦合的动态约束将模型转化为多时段动态规划问题;采用策略迭代动态规划法求解模型得到累加满意度最大的路径,即多目标动态最优潮流的最优解。本发明结合NSGA‑II和近似动态规划法求解多目标动态最优潮流问题,综合处理其中包含连续变量的多目标优化问题与传统动态规划的“维数灾”问题。
Description
技术领域
本发明属于电力***日前调度计划研究领域,涉及一种基于 NSGA-II和近似动态规划的多目标动态最优潮流求解方法。
背景技术
电力***日前调度计划的研究和制定对于整个调度周期内的电力***优化运行与控制具有重要意义。
多目标动态最优潮流是日前调度计划研究的一个分支,旨在机组组合确定的情况下,基于对电力***下一调度周期负荷等情况的合理预测,制定***中可调节手段如机组出力、机端电压、无功补偿投入和需求响应资源等的调度计划,以实现电力***在满足机组爬坡及电压安全等约束下的经济、环保、节能等方面多目标动态优化运行。这类问题是具有大量混合变量和约束条件的复杂多目标非线性优化问题,其求解方法的有效性和快速性一直是研究重点和难点。
传统的电力***动态最优调度模型的求解方法主要是内点法,随着动态问题和优化技术研究的深入和发展,现有的动态优化算法大体可划分为如下三类:一类是数学规划算法,如内点法、梯度投影法、动态规划等;二类是人工智能算法,如遗传算法,差分进化算法,粒子群算法等;三类是混合算法,将前述算法中的两种或以上融合起来,协调发挥各自优势实现更好的优化性能,如粒子群算法与序列二次规划法的结合等,这方面研究尚有长足发展空间。基于贝尔曼最优性原理的传统动态规划法广泛应用于各领域多阶段决策问题的建模和求解,然而传统动态规划基本用于求解优化单目标问题,并且随着优化问题的复杂化与扩大化,求解过程中的“维数灾”问题日益凸显,严重影响优化计算效率。近年来近似动态规划的研究取得进步,通过合理的方法近似逼近动态规划中的值函数可解决“维数灾”,受到广泛关注。最优潮流问题是包含有大量连续变量的非线性优化问题,基于遗传思想的快速非支配排序遗传算法(non-dominated sorting genetic algorithm-II,NSGA-II)具有很好的非线性优化能力和鲁棒性,并且可获取多目标Pareto最优解集,保证最优个体多样性从而为决策者提供不同偏好选择,成为求解多目标优化问题的优秀智能算法之一。综上所述,研究多目标动态最优潮流求解方法,以有效解决多目标动态最优潮流求解中存在的包含连续变量的多目标优化问题与动态规划的“维数灾”问题,对于提升电力***日前调度水平具有积极意义。
发明内容
针对现有技术的不足,本发明“一种基于NSGA-II和近似动态规划的多目标动态最优潮流求解方法”,提出将NSGA-II与近似动态规划结合,利用近似动态规划对多目标动态最优潮流问题进行建模和转换,结合NSGA-II求取Pareto最优解集和最大满意度决策的思想,从而发挥近似动态规划在解决多阶段决策问题的同时避免“维数灾”的优势和快速非支配排序遗传算法在获取多目标Pareto最优解集的优势,获得该问题全局和整个调度周期多目标最优解。
本发明采用如下技术方案:一种基于近似动态规划和快速非支配排序遗传算法的多目标动态最优潮流求解方法,该方法包括如下步骤:
步骤1:获取电力***在下一个完整调度周期的数据,并进行负荷预测。
步骤2:以发电耗费最小、污染物排放量最小、网损最小为目标,考虑电压稳定约束和无功出力约束,对电力***多目标动态最优潮流问题建模。
步骤3:利用NSGA-II求解的多目标静态Pareto最优解集中的非支配解的模糊满意度表示值函数,将单时段变量的连续决策空间离散化,在此基础上计及时段耦合的动态约束,从而将多目标动态最优潮流问题数学模型表示为传统动态规划模型。
步骤4:将传统动态规划问题进一步转化成近似动态规划问题,采用策略迭代近似动态规划求解,用近似值函数代替精确计算值函数,通过迭代过程更新各状态的近似值函数,最终获取累加满意度最大的路径,即为原多目标动态最优潮流问题的解。
附图说明
图1:本发明基于近似动态规划和快速非支配排序遗传算法的多目标动态最优潮流求解方法的整体实施流程图;
图2:本发明的NSGA-II获取多目标最优潮流Pareto解集流程图;
图3:本发明的偏小型满意度函数曲线;
图4:IEEE30节点***的接线示意图;
图5:IEEE30节点***的典型负荷曲线图。
具体实施方式
下面结合附图及具体实施例,对本发明做进一步详述。
本发明提出的基于近似动态规划和快速非支配排序遗传算法的多目标动态最优潮流求解方法,其整体实施流程见图1,下面以 IEEE30节点***为具体实施例对其进行详细说明,其接线情况见图 4。实施例用于说明但不限于本发明。
步骤1:获取电力***在下一个完整调度周期的数据,并进行负荷预测。
对于本实施例直接输入IEEE30节点***的数据即可。负荷预测曲线采用图5典型双峰曲线为例。
步骤2:以发电耗费最小、污染物排放量最小、网损最小为目标,考虑电压稳定约束和无功出力约束,对电力***多目标动态最优潮流问题建模。
本专利建立的多目标动态最优潮流模型如下:
决策变量:
PGi(t)是火电机组i的有功出力(i=1,2,…nG),UGi(t)是火电机组i 所在节点的电压,Bk(t)(k=1,2,…nSC)是无功补偿器k投入量。本实例中,传统火电机组6台,无功补偿器2台,平衡机组有功出力不纳入决策变量,则有13个决策变量。
目标函数:
约束条件:
1.静态约束
(1)功率平衡约束
对于节点bi,需满足有功和无功功率平衡
(2)发电约束
(3)节点电压约束
(4)***旋转备用约束
取μ=5%
2.动态约束
火电机组的爬坡约束
考虑爬坡约束后的火电机组出力上下限由下式决定
步骤3:利用NSGA-II求解的多目标静态Pareto最优解集中的非支配解的模糊满意度表示值函数,将单时段变量的连续决策空间离散化,在此基础上计及时段耦合的动态约束,从而将多目标动态最优潮流问题数学模型表示为传统动态规划模型。
电力***多目标最优潮流问题的多个目标通常彼此矛盾,难以同时得到优化,本发明利用NSGA-II处理该问题,以获取Pareto最优解集,具体求解流程如图2所示。在NSGA-II算法中,种群规模设为100,Pareto分数设为0.35,在优化目标函数中加入罚函数以实现约束。
采用偏小型模糊满意度计算公式,其示意图如图3。对于Pareto 解集中的每个非支配解,计算其每个目标值的满意度
其中fm是目标m的值(m=1,2,…,M),本专利中M=3;fmax,fmin是设定的最大和最小的目标值,取为单时段单目标最优解的目标值。
本专利中N=35,即在Pareto解集中有35个非支配解。第n个非支配解的满意度计算公式为
这样,三个目标综合最优体现为满意度最大,可将动态多目标优化问题转化为累计满意度最大问题。传统动态规划问题多用于处理决策空间是离散化的问题,本专利通过NSGA-II求解和满意度计算的方法,将多目标动态最优潮流这一多阶段决策问题的连续决策空间离散化,根据动态规划理论,多目标动态最优潮流可表述为如下贝尔曼方程:
对于每个{t,St}
V(St)=max[μ(St,at)+V(St+1)] (15)
t表示时段,本专利中对一个调度周期(一天)按小时划分为24 时段;St是第t时段非支配解的信息状态向量;at是第t时段采取的决策,即处于St时选择了当下可行的35个解中的哪一个;St+1是St执行了at后下一时段的决策解信息状态。V(St+1)为St+1的值函数(valuefunction),表征第t时段的调度计划对未来时段满意度产生的影响,取为当t+1时段***状态为St+1时,第t+1到T时段的最大满意度。根据贝尔曼最优性原理,当调度计划{S0,…,ST}使第0到T时段的累加满意度最大时,其第t到T时段优化问题的负的累加满意度一定也是最大的。如此,问题转化为简化的子问题集合,可用递归方式求解。
步骤4:采用策略迭代近似动态规划求解,用近似值函数代替精确计算值函数,通过迭代过程更新各状态的近似值函数,最终获取累加满意度最大的路径,即为原多目标动态最优潮流问题的解。
前文提出的动态多目标调度问题算法若以传统动态规划的逆序求解法求解,具体到实例提出的问题,对需对所有可行的Pareto非支配解进行计算,穷举则存在3524个解路径,即存在“维数灾”问题。本专利采用策略迭代近似动态规划的方法来避免穷举的求解,与传统动态规划逆序求解法不同,策略迭代近似动态规划方法采用正向求解,用近似值函数的方法来指导决策。
对前述动态规划模型作进一步定义:
状态空间:所有可能的Pareto非支配解构成的集合;
采用近似值函数为:
式中αk-1∈[0,1]是步长(或称平滑系数、学习因数)。
选择使得观测满意度达到最大值的动作为预决策状态的最优动作,即决策确定原则如下
通过策略迭代动态规划嵌套NSGA-II寻优,最终选出累加满意度最大的路径,即构成多目标动态最优潮流的解。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的保护范畴。
Claims (4)
1.一种基于近似动态规划和快速非支配排序遗传算法的动态多目标最优潮流求解方法,其特征是,该方法包括如下步骤:
步骤1:获取电力***在下一个完整调度周期的数据,并进行负荷预测;
步骤2:以发电耗费最小、污染物排放量最小、网损最小为目标,考虑电压稳定约束和无功出力约束,对电力***多目标动态最优潮流问题建模;
步骤3:利用NSGA-II求解的多目标静态Pareto最优解集中的非支配解的模糊满意度表示值函数,将单时段变量的连续决策空间离散化,在此基础上计及时段耦合的动态约束,从而将多目标动态最优潮流问题数学模型表示为传统动态规划模型;
步骤4:将传统动态规划问题进一步转化成近似动态规划问题,采用策略迭代方法对近似动态规划问题求解,用近似值函数代替精确计算值函数,通过迭代过程更新各状态的近似值函数,最终获取累加满意度最大的路径,即为原多目标动态最优潮流问题的解。
2.根据权利要求l所述的基于近似动态规划和快速非支配排序遗传算法的动态多目标最优潮流求解方法,其特征是,步骤3中的多目标动态最优潮流问题数学模型转化为传统动态规划模型,是通过单时段NSGA-II寻优获取Pareto解集,将每时段决策空间离散化,从而将多目标动态最优潮流问题数学模型表示为传统动态规划模型,传统动态规划基于贝尔曼最优原理,贝尔曼方程表示为:V(St)=max[μ(St,at)+V(St+1)],式中t表示时段,St是第t时段非支配解的信息状态向量,at是第t时段采取的决策,即处于St时选择了当下可行的35个解中的哪一个,St+1是St执行了at后下一时段的决策解信息状态,V(St)为St的值函数,V(St +1)为下一阶段状态下的值函数。
3.根据权利要求1至2中任意一项权利要求所述的基于近似动态规划和快速非支配排序遗传算法的动态多目标最优潮流求解方法,其特征是,步骤3中的将多目标动态最优潮流问题数学模型表示为传统动态规划模型,是以非支配解的模糊满意度为动态规划问题的值函数,通过获取多时段累加满意度最大的路径来求取多目标动态最优潮流的最优解,在公式V(St)=max[μ(St,at)+V(St+1)]中,当调度计划{S0,…,ST}使第0到T时段的累加满意度最大时,其第t到T时段优化问题的负的累加满意度一定也是最大的,也就是将问题转化为简化的子问题集合,进一步用递归方式求解,式中t表示时段,St是第t时段非支配解的信息状态向量,at是第t时段采取的决策,即处于St时选择了当下可行的35个解中的哪一个,St+1是St执行了at后下一时段的决策解信息状态,V(St)为St的值函数,V(St+1)为下一阶段状态下的值函数。
4.根据权利要求1至2中任意一项权利要求所述的基于近似动态规划和快速非支配排序遗传算法的动态多目标最优潮流求解方法,其特征是,步骤4中通过策略迭代方法,将传统动态规划问题进一步转化成近似动态规划问题并求解,考虑到多目标动态最优潮流问题若以传统动态规划的逆序求解法求解,具体到实例提出的问题,需对所有可行的Pareto非支配解进行计算,穷举则存在3524个解路径,即存在“维数灾”问题,因此步骤4中采用策略迭代近似动态规划的方法来避免“维数灾”问题,与传统动态规划逆序求解法不同,策略迭代近似动态规划方法采用正向求解法,用近似值函数来指导决策,也就是通过策略迭代动态规划嵌套NSGA-II寻优,最终选出累加满意度最大的路径,即构成多目标动态最优潮流问题的解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510605157.1A CN105631528B (zh) | 2015-09-22 | 2015-09-22 | 一种基于nsga-ii和近似动态规划的多目标动态最优潮流求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510605157.1A CN105631528B (zh) | 2015-09-22 | 2015-09-22 | 一种基于nsga-ii和近似动态规划的多目标动态最优潮流求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631528A CN105631528A (zh) | 2016-06-01 |
CN105631528B true CN105631528B (zh) | 2020-01-10 |
Family
ID=56046436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510605157.1A Active CN105631528B (zh) | 2015-09-22 | 2015-09-22 | 一种基于nsga-ii和近似动态规划的多目标动态最优潮流求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631528B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529733A (zh) * | 2016-11-22 | 2017-03-22 | 浙江工业大学 | 基于Gamma Test和NSGA‑II的分布式光伏出力预测输入变量降维方法 |
CN108985532B (zh) * | 2017-06-02 | 2021-05-11 | 上海交通大学 | 基于碳排放的网源荷调度评估***及方法 |
CN107196343B (zh) * | 2017-07-19 | 2020-01-10 | 国家电网公司 | 一种多端柔性直流孤岛电网送端的日前电压计划生成方法 |
CN107832542B (zh) * | 2017-11-22 | 2020-09-11 | 国网河南省电力公司电力科学研究院 | 一种基于时空尺度消纳风光的机组组合优化调度方法 |
CN107944634A (zh) * | 2017-12-11 | 2018-04-20 | 华中科技大学 | 一种时段耦合嵌套的水电站群发电优化调度方法 |
CN108255780B (zh) * | 2017-12-20 | 2021-11-16 | 中国人民解放军海军大连舰艇学院 | 一种基于优化目标的舰队防空资源控制参量的求解计算方法 |
CN108363303B (zh) * | 2018-03-09 | 2020-09-04 | 重庆科技学院 | 基于ar偏好信息的差分进化铝电解多目标优化方法 |
CN109085752B (zh) * | 2018-03-09 | 2020-09-29 | 重庆科技学院 | 基于角度支配关系的铝电解偏好多目标优化算法 |
CN110533263A (zh) * | 2019-09-11 | 2019-12-03 | 绍兴建元电力集团有限公司 | 一种基于改进nsga-ii算法的电-气-热综合***多目标优化方法 |
CN116756808B (zh) * | 2023-05-23 | 2024-04-23 | 济青高速铁路有限公司 | 一种铁路线形优化方法、装置、电子设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867083A (zh) * | 2012-08-30 | 2013-01-09 | 浙江大学 | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8015127B2 (en) * | 2006-09-12 | 2011-09-06 | New York University | System, method, and computer-accessible medium for providing a multi-objective evolutionary optimization of agent-based models |
-
2015
- 2015-09-22 CN CN201510605157.1A patent/CN105631528B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867083A (zh) * | 2012-08-30 | 2013-01-09 | 浙江大学 | 考虑不确定性的压力机滑块机构高刚度轻量化设计方法 |
Non-Patent Citations (1)
Title |
---|
基于NSGA-II的含风电场电力***多目标调度计划研究;罗斌;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20131215(第S2期);摘要,第10-50页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105631528A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105631528B (zh) | 一种基于nsga-ii和近似动态规划的多目标动态最优潮流求解方法 | |
CN112186743B (zh) | 一种基于深度强化学习的动态电力***经济调度方法 | |
Zhang et al. | Multi-objective optimal reactive power dispatch of power systems by combining classification-based multi-objective evolutionary algorithm and integrated decision making | |
Yuan et al. | An extended NSGA-III for solution multi-objective hydro-thermal-wind scheduling considering wind power cost | |
Polimeni et al. | Development and experimental validation of hierarchical energy management system based on stochastic model predictive control for Off-grid Microgrids | |
CN109242188B (zh) | 一种钢铁煤气***长期区间预测及其结构学习方法 | |
CN113572157A (zh) | 一种基于近端策略优化的用户实时自治能量管理优化方法 | |
CN111367349A (zh) | 一种基于预测模型的光伏mppt控制方法及*** | |
He et al. | Stochastic model predictive control of hybrid energy storage for improving AGC performance of thermal generators | |
Gao et al. | An iterative optimization and learning-based IoT system for energy management of connected buildings | |
Chatziagorakis et al. | Enhancement of hybrid renewable energy systems control with neural networks applied to weather forecasting: the case of Olvio | |
CN116914751A (zh) | 一种智能配电控制*** | |
CN115411728A (zh) | 一种融合q学习与势博弈的多微电网***协调控制方法 | |
Gassi et al. | Improving real-time energy decision-making model with an actor-critic agent in modern microgrids with energy storage devices | |
Pacaud et al. | Stochastic optimal control of a domestic microgrid equipped with solar panel and battery | |
Zhou et al. | Optimization of building demand flexibility using reinforcement learning and rule-based expert systems | |
Liu et al. | A novel integral reinforcement learning-based control method assisted by twin delayed deep deterministic policy gradient for solid oxide fuel cell in DC microgrid | |
CN111342471B (zh) | 一种基于机器学习的家庭产消者功率优化管理方法 | |
CN115459320A (zh) | 一种多点分布式储能***聚合控制智能决策方法及装置 | |
CN115115276A (zh) | 考虑不确定性和隐私保护的虚拟电厂调度方法及*** | |
Kousounadis-Knousen et al. | A New Co-Optimized Hybrid Model Based on Multi-Objective Optimization for Probabilistic Wind Power Forecasting in a Spatiotemporal Framework | |
Casagrande et al. | Learning-based MPC using differentiable optimisation layers for microgrid energy management | |
CN114091767A (zh) | 数据驱动的综合能源***鲁棒随机规划的方法及*** | |
Cheng et al. | Reinforcement learning-based scheduling of multi-battery energy storage system | |
Omondi et al. | Smart grid energy management system for industrial applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |