CN109713728A - 微电网控制方法及系统 - Google Patents
微电网控制方法及系统 Download PDFInfo
- Publication number
- CN109713728A CN109713728A CN201910172375.9A CN201910172375A CN109713728A CN 109713728 A CN109713728 A CN 109713728A CN 201910172375 A CN201910172375 A CN 201910172375A CN 109713728 A CN109713728 A CN 109713728A
- Authority
- CN
- China
- Prior art keywords
- film
- grid
- micro
- population
- moment
- 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.)
- Granted
Links
Landscapes
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种微电网控制方法,所述微电网包括分布式电源、储能装置和负荷,所述方法应用于微电网控制系统,所述方法包括建立多目标微电网能量管理模型;通过多目标粒子群膜优化算法得到Pareto最优解集;通过模糊决策选择Pareto最优解集中的折中解;根据所述折中解对微电网的分布式电源和储能装置进行调度。本发明通过建立多目标微电网能量管理模型,并采用多目标粒子群膜优化算法,可以得到优化且合理的控制数据对微电网进行调度,可靠性好,稳定性高。
Description
技术领域
本发明涉及微电网领域,具体涉及微电网控制方法及系统。
背景技术
随着分布式发电和有源分布式网络的不断增加,出现了微电网,微电网是负荷、分布式电源和储能等装置共同组成的供电网络,与传统发电厂相比,微电网具有特定的具体操作和约束条件,到目前为止,在一些发达国家,如日本、加拿大和美国,建立了许多微电网实验平台,近年来,我国很多科研院所也在研究微电网,由于微电网结构复杂,很难获得优化的控制数据对微电网系统进行合理调度。
发明内容
有鉴于此,本发明提供一种微电网控制方法,通过建立微电网多目标能量管理模型,并采用多目标粒子群膜优化算法,可以得到优化且合理的控制数据对微电网进行调度,可靠性好,稳定性高。为解决以上技术问题,本发明提供的技术方案如下:
微电网控制方法,所述微电网包括分布式电源、储能装置和负荷,所述方法应用于微电网控制系统,所述方法包括:
步骤1:建立微电网多目标能量管理模型,包括:
根据公式(1)-(5)建立包括最小化运行成本和污染物排放的微电网多目标能量管理函数minf2obj(x),根据公式(6)-(9)建立微电网运行成本目标函数f1(x),根据公式(10)-(12)建立微电网污染物排放目标函数f2(x),其中:
minf2obj(x)={f1(x),f2(x)} (1)
x=[Pg,Ug]1×2nT (2)
n=NDG+NS+1 (5)
costGrid(t)=uGrid(t)PGrid(t)CGrid(t)(9)
emissionGrid(t)=uGrid(t)PGrid(t)EGrid(t) (12)
其中,x为决策变量,i为分布式电源的编号,j为储能装置的编号,n为决策变量的个数,NDG为分布式电源的数量,NS为储能装置的数量,0≤i≤NDG,0≤j≤NS,Pg为分布式电源、储能和大电网的功率向量,Ug为分布式电源、储能和大电网的状态向量,为第i个分布式电源的功率向量,为第j个储能装置的功率向量,PGrid为大电网的功率向量,ui为第i个分布式电源的状态向量,uj为第j个储存装置的状态向量,uGrid为大电网的状态向量,T为时间长度,costDG(t)为t时刻各种分布式电源发电运行成本总和,costS(t)为t时刻储能装置运行成本总和,costGrid(t)为t时刻微电网与大电网的交互费用总和;ui(t)为第i种分布式电源在t时刻的启停状态,为第i种分布式电源在t时刻的发电功率,为第i种分布式电源在t时刻的发电运行成本费用,为第i种分布式电源启停费用;uj(t)为第j种储能装置在t时刻的启停状态,为第j种储能装置在t时刻的发电功率,为第j种储能装置在t时刻的储能运行成本费用,为第j种储能装置启停费用;uGrid(t)表示大电网在t时刻处于买电状态或卖电状态,PGrid(t)为微电网在t时刻向大电网购买或卖出的电能功率,CGrid(t)为微电网在t时刻向大电网买电或卖电时的电价;
其中,emissionDG(t)为t时刻各种分布式电源发电运行时产生的污染物.排放量,emissionS(t)为t时刻储能装置运行时产生的污染物排放量,emissionGrid(t)为t时刻微电网与大电网交互时产生的污染物排放量,为t时刻第i种分布式电源发电运行产生每兆瓦的功率时的污染物排放量,为t时刻第j个储能装置运行产生每兆瓦的功率时的污染物排放量,EGrid(t)为t时刻微电网与大电网交互产生每兆瓦的功率时的污染物排放量;
步骤2:通过多目标粒子群膜优化算法得到Pareto最优解集;
步骤3:通过模糊决策选择Pareto最优解集中的折中解;
步骤4:根据所述折中解对微电网的分布式电源和储能装置进行调度。
进一步地,建立微电网运行成本目标函数和污染物排放目标函数时满足公式(13)-(15),其中:
其中,Wmin≤Wj(t)≤Wmax,Wj(t)表示储能装置在t时刻的电量状态,和分别表示充电功率和放电功率,η和μ分别表示充电效率和放电效率。
进一步地,多目标粒子群膜优化算法中包括带转运规则的组织型膜系统,根据公式(16)建立组织型膜系统的膜结构,
Π=(O,ω1,,,ωi,,,,ωm,R1,,,Ri,,,Rm,i0) (16)
其中,0是对象的字母表,1≤i≤m,m为细胞个数,ωi表示在初始状态下存在于区域1,2,...,m中的种群,Ri={r1,r2}表示转运规则的有限集合,i0表示输出膜的区域。
进一步地,所述转运规则包括:
膜间进行通信的同向转运规则r1:(i,a,j),其中1≤i,j≤m,a是膜i中当前存在的对象,r1主要负责细胞膜间的信息共享,和
膜与环境进行通信的同向转运规则r2:(i,a,0),其中1≤i≤m,0表示环境,在满足迭代要求时,r2负责将膜内所有对象传输到环境中。
进一步地,步骤2包括:
步骤21:建立多目标粒子群膜优化算法的膜结构为[[]1[]2...[]m]0,其中,初始化参数包括种群大小nPop、外部归档集大小nRep、迭代次数MaxIt、惯性权重ω、个体学习系数c1、全局学习系数c2、膜数量m,0表示环境,并且种群nPop中的每个种群Pop为所述微电网运行成本目标函数的决策变量;
步骤22:在环境0内根据公式(16)初始化种群,
Pop.Vel=0 (16)
其中,Pop.Vel为种群速度;
步骤23:根据公式(17)将每个粒子的历史最优位置设为粒子本身,
Pbest(i)=Pop(i) (17)
其中,i=1,2,3,...,N,N是种群大小;
步骤24:根据所述微电网的运行成本目标函数和污染物排放目标函数计算每个粒子的适应度值;
步骤25:将环境0内已完成初始化的粒子随机分配到各个膜中,每个膜中至少包含一个粒子,环境0内不含粒子,构成并保存为每个膜内的父代种群;
步骤26:将每个膜中的父代种群按公式(18)进行速度位置更新,获得每个膜中的子种群,
Pop(i).Vel=ω×Pop(i).Vel+k1×(Pbest(i)-Pop(i))
+k2×(REP(h)-Pop(i))
Pop(i)=Pop(i)+Pop(i).Vel (18)
其中,ω是惯性权重,k1、k2是[0,1]之间的随机数;
步骤27:执行膜间进行通信的同向转运规则r1,将每层膜内的子种群传递到下一层膜内,与下一层膜内的子种群合并并进行速度位置更新,得到每个膜内的新子种群;
步骤28:将每层膜内的新子种群与其父代种群进行合并,再次计算每个粒子的适应度值,保留所有非支配粒子,得到备选最优解集;
步骤29:判断是否满足终止条件,若是,执行步骤210,否则,执行步骤26;
步骤210:执行膜与环境进行通信的同向转运规则r2,将膜1~m更新后得到的备选最优解集传送到环境0中,将环境0中所有非支配粒子存入外部归档集中,消除外部归档集中所有支配的粒子,输出Pareto最优解集。
进一步地,步骤3包括:采用具有隶属度函数的模糊决策从Pareto最优解集中得到折中解,其中:
隶属度函数为公式(19):
其中,表示目标函数β中的第α个解的隶属度,和是非支配解集中的最大值和最小值,每个解α的标准化隶属度函数可表示为公式(20):
其中,Nobj为目标函数个数,Np为种群数目,1≤β≤Nobj,1≤α≤Np,依次计算出值最大的解,则为选出的折中解。
微电网控制系统,包括:
微电网多目标能量管理模型建立模块:用于设置微电网多目标能量管理关于成本和污染物排放的控制目标函数;
多目标粒子群膜优化模块:用于通过多目标粒子群膜优化算法得到Pareto最优解集;
模糊决策模块:用于通过模糊决策方法选择Pareto最优解集中的折中解;
调度模块:用于根据所述折中解对微电网的分布式电源和储能装置进行调度。
进一步地,所述多目标粒子群膜优化模块包括:
参数初始化模块:用于接收微电网能量管理适应度值参数,微电网能量管理适应度值参数包括:24小时光伏、风机预测发电量、重要负荷24小时预测用电量、24小时预测实时电价和微电网各个微源的功率上下限、发电价格及污染物排放量;
初始化模块:用于初始化算法的相关参数及膜结构,算法和膜结构参数包括:种群大小nPop、外部归档集大小nRep、迭代次数MaxIt、惯性权重ω、个体学习系数c1、全局学习系数c2、膜数量m;并初始化种群的速度、位置与适应度值,将已完成初始化的粒子随机分配到各个膜中,每个膜中至少保证含有一个粒子,环境中不含粒子以获得并保存为每个膜中的父代种群;
粒子更新模块:用于更新每个粒子的位置与速度以获得每个膜中的子种群;
膜间进行通信的同向转运规则模块:用于将每层膜内的子种群传递到下一层膜内,与下一层膜内的子种群合并,并进行速度位置更新,得到并保存每个膜内的新子种群;
多目标粒子群膜优化算法计算模块:用于将每个膜中的新子种群和父代种群进行合并,计算种群中每个粒子的适应度值,保留所有非支配粒子,得到备选最优解集;
终止条件模块:判断是否达到终止条件,如果满足,则触发膜与环境进行通信的同向转运规则模块;否则,触发粒子更新模块;
膜与环境进行通信的同向转运规则模块:整合所有膜中的备选最优解集并传送到环境中,在环境中保留所有非支配解,得到Pareto最优解集。
本发明提供的微电网控制方法,通过建立多目标微电网能量管理模型,并采用多目标粒子群优化算法和模糊决策算法,可以得到优化且合理的控制数据对微电网进行调度。
附图说明
图1为根据一示例性实施例示出的一种微电网控制方法流程图。
图2为根据一示例性实施例示出的一种组织型膜系统的膜结构图。
图3为根据一示例性实施例示出的一种通过多目标粒子群优化算法得到Pareto最优解集的方法流程图。
图4为根据一示例性实施例示出的一种微电网控制系统结构框图。
具体实施方式
为了使本领域的技术人员更好地理解本发明的技术方案,下面结合附图和具体实施例对本发明作进一步的详细说明。
实施例1
如图1所示,本实施例提供一种微电网控制方法,所述微电网包括分布式电源、储能装置和负荷,所述方法应用于微电网控制系统,所述方法包括:
步骤1:建立微电网多目标能量管理模型,包括:
根据公式(1)-(5)建立包括最小化运行成本和污染物排放的微电网多目标能量管理函数minf2obj(x),根据公式(6)-(9)建立微电网运行成本目标函数f1(x),根据公式(10)-(12)建立微电网污染物排放目标函数f2(x),其中:
minf2obj(x)={f1(x),f2(x)} (1)
x=[Pg,Ug]1×2nT (2)
n=NDG+NS+1 (5)
costGrid(t)=uGrid(t)PGrid(t)CGrid(t)(9)
emissionGrid(t)=uGrid(t)PGrid(t)EGrid(t) (12)
其中,x为决策变量,i为分布式电源的编号,j为储能装置的编号,n为决策变量的个数,NDG为分布式电源的数量,NS为储能装置的数量,0≤i≤NDG,0≤j≤NS,Pg为分布式电源、储能和大电网的功率向量,Ug为分布式电源、储能和大电网的状态向量,为第i个分布式电源的功率向量,为第j个储能装置的功率向量,PGrid为大电网的功率向量,ui为第i个分布式电源的状态向量,uj为第j个储存装置的状态向量,uGrid为大电网的状态向量,T为时间长度,costDG(t)为t时刻各种分布式电源发电运行成本总和,costS(t)为t时刻储能装置运行成本总和,costGrid(t)为t时刻微电网与大电网的交互费用总和;ui(t)为第i种分布式电源在t时刻的启停状态,为第i种分布式电源在t时刻的发电功率,为第i种分布式电源在t时刻的发电运行成本费用,为第i种分布式电源启停费用;uj(t)为第j种储能装置在t时刻的启停状态,为第j种储能装置在t时刻的发电功率,为第j种储能装置在t时刻的储能运行成本费用,为第j种储能装置启停费用;uGrid(t)表示大电网在t时刻处于买电状态或卖电状态,PGrid(t)为微电网在t时刻向大电网购买或卖出的电能功率,CGrid(t)为微电网在t时刻向大电网买电或卖电时的电价;
其中,emissionDG(t)为t时刻各种分布式电源发电运行时产生的污染物排放量,emissions(t)为t时刻储能装置运行时产生的污染物排放量,emissionGrid(t)为t时刻微电网与大电网交互时产生的污染物排放量,为t时刻第i种分布式电源发电运行产生每兆瓦的功率时的污染物排放量,为t时刻第j个储能装置运行产生每兆瓦的功率时的污染物排放量,EGrid(t)为t时刻微电网与大电网交互产生每兆瓦的功率时的污染物排放量;实施本实施例时主要考虑三种主要的污染物:二氧化碳、二氧化硫和氮氧化物;
步骤2:通过多目标粒子群膜优化算法得到Pareto最优解集;
步骤3:通过模糊决策选择Pareto最优解集中的折中解;
步骤4:根据所述折中解对微电网的分布式电源和储能装置进行调度。
作为优选,本实施例在建立微电网运行成本目标函数和污染物排放目标函数时满足公式(13)-(15),其中:
其中,Wmin≤Wj(t)≤Wmax,Wj(t)表示储能装置在t时刻的电量状态,和分别表示充电功率和放电功率,η和μ分别表示充电效率和放电效率。
作为优选,如图2所示,多目标粒子群膜优化算法中包括带转运规则的组织型膜系统,根据公式(16)建立组织型膜系统的膜结构,
Π=(O,ω1,,,ωi,,,ωm,R1,,,Ri,,,Rm,i0) (16)
其中,0是对象的字母表,1≤i≤m,m为细胞个数,ωi表示在初始状态下存在于区域1,2,...,m中的种群,Ri={r1,r2}表示转运规则的有限集合,i0表示输出膜的区域。
作为优选,所述转运规则包括:
膜间进行通信的同向转运规则r1:(i,a,j),其中1≤i,j≤m,a是膜i中当前存在的对象,r1主要负责细胞膜间的信息共享,和
膜与环境进行通信的同向转运规则r2:(i,a,0),其中1≤i≤m,0表示环境,在满足迭代要求时,r2负责将膜内所有对象传输到环境中。
作为优选,如图3所示,步骤2包括:
步骤21:建立多目标粒子群膜优化算法的膜结构为[[]1[]2...[]m]0,其中,初始化参数包括种群大小nPop、外部归档集大小nRep、迭代次数MaxIt、惯性权重ω、个体学习系数c1、全局学习系数c2、膜数量m,0表示环境,并且种群nPop中的每个种群Pop为所述微电网运行成本目标函数的决策变量;
步骤22:在环境0内根据公式(16)初始化种群,
Pop.Vel=0 (16)
其中,Pop.Vel为种群速度;
步骤23:根据公式(17)将每个粒子的历史最优位置设为粒子本身,
Pbest(i)=Pop(i) (17)
其中,i=1,2,3,...,N,N是种群大小;
步骤24:根据所述微电网的运行成本目标函数和污染物排放目标函数计算每个粒子的适应度值;
步骤25:将环境0内已完成初始化的粒子随机分配到各个膜中,每个膜中至少包含一个粒子,环境0内不含粒子,构成并保存为每个膜内的父代种群;
步骤26:将每个膜中的父代种群按公式(18)进行速度位置更新,获得每个膜中的子种群,
Pop(i).Vel=ω×Pop(i).Vel+k1×(Pbest(i)-Pop(i))
+k2×(REP(h)-Pop(i))
Pop(i)=Pop(i)+Pop(i).Vel (18)
其中,ω是惯性权重,k1、k2是[0,1]之间的随机数;
步骤27:执行膜间进行通信的同向转运规则r1,将每层膜内的子种群传递到下一层膜内,与下一层膜内的子种群合并并进行速度位置更新,得到每个膜内的新子种群;
步骤28:将每层膜内的新子种群与其父代种群进行合并,再次计算每个粒子的适应度值,保留所有非支配粒子,得到备选最优解集;
步骤29:判断是否满足终止条件,若是,执行步骤210,否则,执行步骤26;
步骤210:执行膜与环境进行通信的同向转运规则r2,将膜1~m更新后得到的备选最优解集传送到环境0中,将环境0中所有非支配粒子存入外部归档集中,消除外部归档集中所有支配的粒子,输出Pareto最优解集。
作为优选,步骤3包括:采用具有隶属度函数的模糊决策从Pareto最优解集中得到折中解,其中:
隶属度函数为公式(19):
其中,表示目标函数β中的第α个解的隶属度,和是非支配解集中的最大值和最小值,每个解α的标准化隶属度函数可表示为公式(20):
其中,Nobj为目标函数个数,Np为种群数目,1≤β≤Nobj,1≤α≤Np,依次计算出值最大的解,则为选出的折中解。
实施例2
如图4所示,微电网控制系统,包括:
微电网多目标能量管理模型建立模块:用于设置微电网多目标能量管理关于成本和污染物排放的控制目标函数;
多目标粒子群膜优化模块:用于通过多目标粒子群膜优化算法得到Pareto最优解集;
模糊决策模块:用于通过模糊决策方法选择Pareto最优解集中的折中解;
调度模块:用于根据所述折中解对微电网的分布式电源和储能装置进行调度。
作为优选,所述多目标粒子群膜优化模块包括:
参数初始化模块:用于接收微电网能量管理适应度值参数,微电网能量管理适应度值参数包括:24小时光伏、风机预测发电量、重要负荷24小时预测用电量、24小时预测实时电价和微电网各个微源的功率上下限、发电价格及污染物排放量;
初始化模块:用于初始化算法的相关参数及膜结构,算法和膜结构参数包括:种群大小nPop、外部归档集大小nRep、迭代次数MaxIt、惯性权重ω、个体学习系数c1、全局学习系数c2、膜数量m;并初始化种群的速度、位置与适应度值,将已完成初始化的粒子随机分配到各个膜中,每个膜中至少保证含有一个粒子,环境中不含粒子以获得并保存为每个膜中的父代种群;
粒子更新模块:用于更新每个粒子的位置与速度以获得每个膜中的子种群;
膜间进行通信的同向转运规则模块:用于将每层膜内的子种群传递到下一层膜内,与下一层膜内的子种群合并,并进行速度位置更新,得到并保存每个膜内的新子种群;
多目标粒子群膜优化算法计算模块:用于将每个膜中的新子种群和父代种群进行合并,计算种群中每个粒子的适应度值,保留所有非支配粒子,得到备选最优解集;
终止条件模块:判断是否达到终止条件,如果满足,则触发膜与环境进行通信的同向转运规则模块;否则,触发粒子更新模块;
膜与环境进行通信的同向转运规则模块:整合所有膜中的备选最优解集并传送到环境中,在环境中保留所有非支配解,得到Pareto最优解集。
以上仅是本发明的优选实施方式,应当指出的是,上述优选实施方式不应视为对本发明的限制,本发明的保护范围应当以权利要求所限定的范围为准。对于本技术领域的普通技术人员来说,在不脱离本发明的精神和范围内,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.微电网控制方法,其特征在于,所述微电网包括分布式电源、储能装置和负荷,所述方法应用于微电网控制系统,所述方法包括:
步骤1:建立微电网多目标能量管理模型,包括:
根据公式(1)-(5)建立包括最小化运行成本和污染物排放的微电网多目标能量管理函数minf2obj(x),根据公式(6)-(9)建立微电网运行成本目标函数f1(x),根据公式(10)-(12)建立微电网污染物排放目标函数f2(x),其中:
minf2obj(x)={f1(x),f2(x)} (1)
x=[Pg,Ug]1×2nT (2)
n=NDG+NS+1 (5)
cosrGrid(t)=uGrid(t)PGrid(t)CGrid(t) (9)
emissionGrid(t)=uGrid(t)PGrid(t)EGrid(t) (12)
其中,x为决策变量,i为分布式电源的编号,j为储能装置的编号,n为决策变量的个数,NDG为分布式电源的数量,NS为储能装置的数量,0≤i≤NDG,0≤j≤NS,Pg为分布式电源、储能和大电网的功率向量,Ug为分布式电源、储能和大电网的状态向量,为第i个分布式电源的功率向量,为第j个储能装置的功率向量,PGrid为大电网的功率向量,ui为第i个分布式电源的状态向量,uj为第j个储存装置的状态向量,uGrid为大电网的状态向量,T为时间长度,costDG(t)为t时刻各种分布式电源发电运行成本总和,costS(t)为t时刻储能装置运行成本总和,costGrid(t)为t时刻微电网与大电网的交互费用总和;ui(t)为第i种分布式电源在t时刻的启停状态,为第i种分布式电源在t时刻的发电功率,为第i种分布式电源在t时刻的发电运行成本费用,为第i种分布式电源启停费用;uj(t)为第j种储能装置在t时刻的启停状态,为第j种储能装置在t时刻的发电功率,为第j种储能装置在t时刻的储能运行成本费用,为第j种储能装置启停费用;uGrid(t)表示大电网在t时刻处于买电状态或卖电状态,PGrid(t)为微电网在t时刻向大电网购买或卖出的电能功率,CGrid(t)为微电网在t时刻向大电网买电或卖电时的电价;
其中,emissionDG(t)为t时刻各种分布式电源发电运行时产生的污染物排放量,emissionS(t)为t时刻储能装置运行时产生的污染物排放量,emissionGrid(t)为t时刻微电网与大电网交互时产生的污染物排放量,为t时刻第i种分布式电源发电运行产生每兆瓦的功率时的污染物排放量,为t时刻第j个储能装置运行产生每兆瓦的功率时的污染物排放量,EGrid(t)为t时刻微电网与大电网交互产生每兆瓦的功率时的污染物排放量;
步骤2:通过多目标粒子群膜优化算法得到Pareto最优解集;
步骤3:通过模糊决策选择Pareto最优解集中的折中解;
步骤4:根据所述折中解对微电网的分布式电源和储能装置进行调度。
2.根据权利要求1所述的微电网控制方法,其特征在于,建立微电网运行成本目标函数和污染物排放目标函数时满足公式(13)—(15),其中:
其中,PL(t)表示t时刻负荷的需求功率,Wmin≤Wj(t)≤Wmax,Wj(t)表示储能装置在t时刻的电量状态,和分别表示充电功率和放电功率,η和μ分别表示充电效率和放电效率。
3.根据权利要求2所述的微电网控制方法,其特征在于,多目标粒子群膜优化算法中包括带转运规则的组织型膜系统,根据公式(16)建立组织型膜系统的膜结构,
∏=(O,ω1,,,ωi,,,ωm,R1,,,Ri,,,Rm,i0) (16)
其中,O是对象的字母表,1≤i≤m,m为细胞个数,ωi表示在初始状态下存在于区域1,2,...,m中的种群,Ri={r1,r2}表示转运规则的有限集合,i0表示输出膜的区域。
4.根据权利要求3所述的微电网控制方法,其特征在于,所述转运规则包括:
膜间进行通信的同向转运规则r1:(i,a,j),其中1≤i,j≤m,a是膜i中当前存在的对象,r1主要负责细胞膜间的信息共享,和
膜与环境进行通信的同向转运规则r2:(i,a,0),其中1≤i≤m,0表示环境,在满足迭代要求时,r2负责将膜内所有对象传输到环境中。
5.根据权利要求4所述的微电网控制方法,其特征在于,步骤2包括:
步骤21:建立多目标粒子群膜优化算法的膜结构为[[]1[]2...[]m]0,其中,初始化参数包括种群大小nPop、外部归档集大小nRep、迭代次数MaxIt、惯性权重ω、个体学习系数c1、全局学习系数c2、膜数量m,0表示环境,并且种群nPop中的每个种群Pop为所述微电网运行成本目标函数的决策变量;
步骤22:在环境0内根据公式(16)初始化种群,
Pop.Vel=0 (16)
其中,Pop.Vel为种群速度;
步骤23:根据公式(17)将每个粒子的历史最优位置设为粒子本身,
Pbest(i)=Pop(i) (17)
其中,i=1,2,3,...,N,N是种群大小;
步骤24:根据所述微电网的运行成本目标函数和污染物排放目标函数计算每个粒子的适应度值;
步骤25:将环境0内已完成初始化的粒子随机分配到各个膜中,每个膜中至少包含一个粒子,环境0内不含粒子,构成并保存为每个膜内的父代种群;
步骤26:将每个膜中的父代种群按公式(18)进行速度位置更新,获得每个膜中的子种群,
Pop(i).Vel=ω×Pop(i).Vel+k1×(Pbest(i)-Pop(i))+k2×(REP(h)-Pop(i))
nPop(i)=Pop(i)+Pop(i).Vel (18)
其中,ω是惯性权重,k1、k2是[0,1]之间的随机数;
步骤27:执行膜间进行通信的同向转运规则r1,将每层膜内的子种群传递到下一层膜内,与下一层膜内的子种群合并并进行速度位置更新,得到每个膜内的新子种群;
步骤28:将每层膜内的新子种群与其父代种群进行合并,再次计算每个粒子的适应度值,保留所有非支配粒子,得到备选最优解集;
步骤29:判断是否满足终止条件,若是,执行步骤210,否则,执行步骤26;
步骤210:执行膜与环境进行通信的同向转运规则r2,将膜1~m更新后得到的备选最优解集传送到环境0中,将环境0中所有非支配粒子存入外部归档集中,消除外部归档集中所有支配的粒子,输出Pareto最优解集。
6.根据权利要求5所述的微电网控制方法,其特征在于,步骤3包括:采用具有隶属度函数的模糊决策从Pareto最优解集中得到折中解,其中:
隶属度函数为公式(19):
其中,表示目标函数β中的第α个解的隶属度,和是非支配解集中的最大值和最小值,每个解α的标准化隶属度函数可表示为公式(20):
其中,Nobj为目标函数个数,Np为种群数目,1≤β≤Nobj,1≤α≤Np,依次计算出值最大的解,则为选出的折中解。
7.微电网控制系统,其特征在于,包括:
微电网多目标能量管理模型建立模块:用于设置微电网多目标能量管理关于成本和污染物排放的控制目标函数;
多目标粒子群膜优化模块:用于通过多目标粒子群膜优化算法得到Pareto最优解集;
模糊决策模块:用于通过模糊决策方法选择Pareto最优解集中的折中解;
调度模块:用于根据所述折中解对微电网的分布式电源和储能装置进行调度。
8.根据权利要求7所述的微电网控制系统,其特征在于,所述多目标粒子群膜优化模块包括:
参数初始化模块:用于接收微电网能量管理适应度值参数,微电网能量管理适应度值参数包括:24小时光伏、风机预测发电量、重要负荷24小时预测用电量、24小时预测实时电价和微电网各个微源的功率上下限、发电价格及污染物排放量;
初始化模块:用于初始化算法的相关参数及膜结构,算法和膜结构参数包括:种群大小nPop、外部归档集大小nRep、迭代次数MaxIt、惯性权重ω、个体学习系数c1、全局学习系数c2、膜数量m;并初始化种群的速度、位置与适应度值,将已完成初始化的粒子随机分配到各个膜中,每个膜中至少保证含有一个粒子,环境中不含粒子以获得并保存为每个膜中的父代种群;
粒子更新模块:用于更新每个粒子的位置与速度以获得每个膜中的子种群;
膜间进行通信的同向转运规则模块:用于将每层膜内的子种群传递到下一层膜内,与下一层膜内的子种群合并,并进行速度位置更新,得到并保存每个膜内的新子种群;
多目标粒子群膜优化算法计算模块:用于将每个膜中的新子种群和父代种群进行合并,计算种群中每个粒子的适应度值,保留所有非支配粒子,得到备选最优解集;
终止条件模块:判断是否达到终止条件,如果满足,则触发膜与环境进行通信的同向转运规则模块;否则,触发粒子更新模块;
膜与环境进行通信的同向转运规则模块:整合所有膜中的备选最优解集并传送到环境中,在环境中保留所有非支配解,得到Pareto最优解集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910172375.9A CN109713728B (zh) | 2019-03-07 | 2019-03-07 | 微电网控制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910172375.9A CN109713728B (zh) | 2019-03-07 | 2019-03-07 | 微电网控制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109713728A true CN109713728A (zh) | 2019-05-03 |
CN109713728B CN109713728B (zh) | 2021-11-05 |
Family
ID=66266461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910172375.9A Expired - Fee Related CN109713728B (zh) | 2019-03-07 | 2019-03-07 | 微电网控制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109713728B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111985811A (zh) * | 2020-08-20 | 2020-11-24 | 成都工业学院 | 一种面向电动汽车的多微网协调控制方法 |
CN113962828A (zh) * | 2021-10-26 | 2022-01-21 | 长春工程学院 | 考虑碳消耗的综合能源系统协调调度方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100891513B1 (ko) * | 2008-08-18 | 2009-04-06 | 주식회사 케이디파워 | 태양광 및 배터리 시스템을 이용한 계통 연계형 하이브리드발전 시스템 및 이를 이용한 발전 방법 |
CN107609693A (zh) * | 2017-08-31 | 2018-01-19 | 安徽大学 | 基于Pareto档案粒子群算法的微电网多目标优化方法 |
CN108347062A (zh) * | 2018-01-12 | 2018-07-31 | 华南理工大学 | 基于势博弈的微电网能量管理分布式多目标协同优化算法 |
CN109217351A (zh) * | 2017-07-06 | 2019-01-15 | 华北电力大学(保定) | 兼顾不同主体利益的主动配电网优化调度方法 |
CN109327042A (zh) * | 2018-09-27 | 2019-02-12 | 南京邮电大学 | 一种微电网多能源联合优化调度方法 |
-
2019
- 2019-03-07 CN CN201910172375.9A patent/CN109713728B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100891513B1 (ko) * | 2008-08-18 | 2009-04-06 | 주식회사 케이디파워 | 태양광 및 배터리 시스템을 이용한 계통 연계형 하이브리드발전 시스템 및 이를 이용한 발전 방법 |
CN109217351A (zh) * | 2017-07-06 | 2019-01-15 | 华北电力大学(保定) | 兼顾不同主体利益的主动配电网优化调度方法 |
CN107609693A (zh) * | 2017-08-31 | 2018-01-19 | 安徽大学 | 基于Pareto档案粒子群算法的微电网多目标优化方法 |
CN108347062A (zh) * | 2018-01-12 | 2018-07-31 | 华南理工大学 | 基于势博弈的微电网能量管理分布式多目标协同优化算法 |
CN109327042A (zh) * | 2018-09-27 | 2019-02-12 | 南京邮电大学 | 一种微电网多能源联合优化调度方法 |
Non-Patent Citations (2)
Title |
---|
MOHAMMAD H. MORADI等: ""Operational Strategy Optimization in an Optimal Sized Smart Microgrid"", 《IEEE TRANSACTIONS ON SMART GRID》 * |
曾嶒等: ""基于鸟群算法的微电网多目标运行优化"", 《电力系统保护与控制》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111985811A (zh) * | 2020-08-20 | 2020-11-24 | 成都工业学院 | 一种面向电动汽车的多微网协调控制方法 |
CN111985811B (zh) * | 2020-08-20 | 2024-02-02 | 成都工业学院 | 一种面向电动汽车的多微网协调控制方法 |
CN113962828A (zh) * | 2021-10-26 | 2022-01-21 | 长春工程学院 | 考虑碳消耗的综合能源系统协调调度方法 |
CN113962828B (zh) * | 2021-10-26 | 2024-05-10 | 长春工程学院 | 考虑碳消耗的综合能源系统协调调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109713728B (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Boaro et al. | Adaptive dynamic programming algorithm for renewable energy scheduling and battery management | |
CN105591406B (zh) | 一种基于非合作博弈的微电网能量管理系统的优化算法 | |
Kyriakarakos et al. | A fuzzy cognitive maps–petri nets energy management system for autonomous polygeneration microgrids | |
Niknam et al. | Impact of thermal recovery and hydrogen production of fuel cell power plants on distribution feeder reconfiguration | |
CN106058855A (zh) | 协调储能与柔性负荷的主动配电网多目标优化调度方法 | |
Parouha et al. | DPD: An intelligent parallel hybrid algorithm for economic load dispatch problems with various practical constraints | |
CN106329568B (zh) | 户商型光伏发电经济调度控制系统 | |
CN110138006A (zh) | 考虑含有新能源电动汽车的多微电网协调优化调度方法 | |
CN109713728A (zh) | 微电网控制方法及系统 | |
CN113887858A (zh) | 一种基于cnn-lstm负荷预测的充电站微电网系统优化调度方法 | |
Gbadega et al. | A hybrid constrained Particle Swarm Optimization-Model Predictive Control (CPSO-MPC) algorithm for storage energy management optimization problem in micro-grid | |
Abedinia et al. | Multiobjective environmental and economic dispatch using imperialist competitive algorithm | |
Bandopadhyay et al. | Application of hybrid multi-objective moth flame optimization technique for optimal performance of hybrid micro-grid system | |
CN107273968A (zh) | 一种基于动态模糊混沌粒子群的多目标调度方法及装置 | |
Shojaeiyan et al. | A novel bio-inspired stochastic framework to solve energy management problem in hybrid AC-DC microgrids with uncertainty | |
Maruf et al. | An enhancement of agent-based power supply-demand matching by using ann-based forecaster | |
CN116667325A (zh) | 一种基于改进布谷鸟算法的微电网并网运行优化调度方法 | |
Khorram-Nia et al. | Optimal switching in reconfigurable microgrids considering electric vehicles and renewable energy sources | |
Rahman et al. | Novel metaheuristic optimization strategies for plug-in hybrid electric vehicles: A holistic review | |
Mellouk et al. | Overview of mathematical methods for energy management optimization in smart grids | |
Zhang et al. | Dynamic control of wind/photovoltaic hybrid power systems based on an advanced particle swarm optimization | |
CN104092209A (zh) | 基于实时反馈的互动型微电网能量控制处理方法 | |
CN111178642A (zh) | 一种基于模拟退火粒子群算法的微电网优化方法 | |
Zhang et al. | Hierarchical energy management scheme for residential communities under grid outage event | |
Al-Messabi et al. | Heuristically enhanced dynamic neural networks for structurally improving photovoltaic power forecasting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211105 |
|
CF01 | Termination of patent right due to non-payment of annual fee |