CN109193754B - 一种分布式光伏规划方法及系统 - Google Patents

一种分布式光伏规划方法及系统 Download PDF

Info

Publication number
CN109193754B
CN109193754B CN201811019913.2A CN201811019913A CN109193754B CN 109193754 B CN109193754 B CN 109193754B CN 201811019913 A CN201811019913 A CN 201811019913A CN 109193754 B CN109193754 B CN 109193754B
Authority
CN
China
Prior art keywords
scene
population
power supply
planning model
photovoltaic power
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
Application number
CN201811019913.2A
Other languages
English (en)
Other versions
CN109193754A (zh
Inventor
曾伟
孙旻
陈波
何伟
黄扬琪
熊俊杰
余侃胜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811019913.2A priority Critical patent/CN109193754B/zh
Publication of CN109193754A publication Critical patent/CN109193754A/zh
Application granted granted Critical
Publication of CN109193754B publication Critical patent/CN109193754B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H02J3/383
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers

Landscapes

  • Photovoltaic Devices (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

一种分布式光伏规划方法及系统,方法包括获取配电网系统参数;将所述电网系统参数代入预先建立的双层规划模型,采用分子微分进化算法进行求解,以求解结果作为规划方案;其中所述双层规划模型包括外层规划模型和内层规划模型,所述外层规划模型为以投资者对分布式光伏电源的年净收益最大为目标,所述内层规划模型以分布式光伏电源的年消纳率最大为目标。系统包括数据获取模块和计算模块。本发明运用了分子微分进化算法来求解双层多场景规划模型,具有收敛速度快,局部搜索能力强的优点,又能提高种群的多样性,增强其全局寻优能力。本发明分布式光伏规划方法可为配电网中分布式光伏电源的合理规划提供科学指导。

Description

一种分布式光伏规划方法及系统
技术领域
本发明涉及一种分布式光伏规划方法及系统,属电力新能源规划技术领域。
背景技术
分布式光伏作为一种清洁能源,接入配电网中不仅可以有效缓解当地的电力紧张格局,还能为光伏投资者带来收益。然而由于传统光伏电源规划时未计及光伏的消纳率,外加国家补贴政策的激励以及不断下降的光伏电源安装成本等因素,常会使得配电网中光伏电源的容量规划值过高。在实际运行中光伏电源渗透率过高往往会使得光伏上网电量过多,导致配电网中出现节点电压越限、系统潮流倒送等电能质量问题。为了确保配电网系统能够安全稳定地运行,光伏电源常被要求弃光限电来减少其上网电量,结果既造成了能源的浪费,又会导致投资者的收益受损和投资的积极性下降。现有规划模型默认了分布式电源所发出的电量能被全额消纳,实际上在对分布式电源进行规划时,如果忽略了其消纳率,则会导致规划结果存在较大的偏差。
发明内容
本发明的目的是,为了克服现有技术存在的上述缺陷,本发明提出了一种分布式光伏规划方法及系统,使得分布式光伏电源容量规划更为合理准确,减少限电弃光现象的发生,同时提高了投资者的收益。
实现本发明目的的技术方案是,一种分布式光伏规划方法,包括获取配电网系统参数;将所述电网系统参数代入预先建立的双层规划模型,采用分子微分进化算法进行求解,以求解结果作为规划方案;其中所述双层规划模型包括外层规划模型和内层规划模型,所述外层规划模型为以投资者对分布式光伏电源的年净收益最大为目标,所述内层规划模型以分布式光伏电源的年消纳率最大为目标。
所述分子微分进化算法对双层规划模型求解,其步骤如下:
(1)取配电网系统线路参数,以及全年的负荷和光照强度数据;采用所述改进K-均值聚类法对全年的“负荷-光照强度”场景进行缩减,得到缩减后指定数量场景的负荷和光照强度数据及其对应的概率;
(2)设置外层分子微分进化算法的基本参数,其中包括种群的规模NP、维度D,最大迭代次数Gmax、分子间的临界距离r0、变异因子F以及交叉因子CR,令外层进化代数G=0;对外层规划模型中的种群变量即光伏电源的安装容量进行初始化,初始种群为S=[X1,X2,…,XNP],其中的第i个种群个体表示为Xi=[xi,1,xi,2,…,xi,D],第i个个体各维的数据随机产生:
Figure BDA0001787048450000021
其中,d=1,2,…,D;randi(·)为产生均匀分布的随机整数函数;而L则为分布式光伏电源的单位额定容量;xd,max为第d维数据的最大值,xd,min为第d维数据的最小值;
(3)将外层规划模型中的光伏电源安装容量值传递到内层,设置内层分子微分进化算法的基本参数,参数类型与外层相同;求解内层规划模型中内层目标函数的最大值;最后计算出种群个体每个场景对应的消纳量然后传递回外层规划模型;
(4)根据外层光伏电源的初始容量并结合内层中求解的光伏消纳量计算外层规划模型的初始目标函数值;
(5)计算G=G+1,对外层种群进行变异、交叉以及校核操作来更新种群;
(6)种群变异,根据构造的变异机制对每个目标向量Xi,G进行变异,产生变异向量Yi,G+1
Yi,G+1=Xbest,G+F(Xr1,G-Xr2,G+Vr1,r2,G);
其中,Xbest,G为进化到第G代时种群的最优个体向量;Xr1,G、Xr2,G分别为第G代时的第r1、r2个个体向量,且i≠best≠r1≠r2;F为变异因子;Vr1,r2,G为第G代时种群个体之间的作用势能;
Figure BDA0001787048450000031
其中,εr1,r2,d为进化到第G代时种群个体向量之间的维差距;ε0为分子间的临界距离;vr1,r2,d为种群个体在d维度时的作用势能;εd为进化到第G代时种群个体向量之间的维差距;
(7)种群交叉,将外层种群的变异向量Yi,G+1与目标向量Xi,G的各维分量进行随机重组,产生新向量Zi,G+1
Figure BDA0001787048450000032
其中,R为[0,1]之间的随机数;CR为交叉因子;zi,d为新向量Zi,G+1第i个个体d维数据;xi,d为第i个种群个体d维数据;yi,d为变异向量第i个个体d维数据;
(8)种群校核,变异交叉操作后部分种群个体可能会出现越界的情况,按下式进行校核:
Figure BDA0001787048450000041
(9)种群选择,计算校核后的外层种群目标函数值,此时需要再次进入内层规划模型中进行寻优,方法参考步骤(3),最后按照下式做出贪婪选择:
Figure BDA0001787048450000042
(10)若G<Gmax,则返回到步骤(5),进入下一次循环;否则终止寻优过程,输出外层目标的最优值以及相应的光伏电源最佳规划方案。
所述采用改进K-均值聚类法对全年的场景进行缩减,采用如下步骤:
(1)假设用ωi(i=1,2,…,8760)表示缩减前全年不同时段的“负荷-光照强度”场景,ωi=[ωi,1i,2],ωi,1、ωi,2分别表示为i时刻的负荷和光照强度值,K则为指定的缩减场景数;分别计算全年场景集中每个场景与剩余场景的距离di,j、每个场景的平均差异度Mi,ave和全年场景集的总体平均差异度Mave
Figure BDA0001787048450000043
Figure BDA0001787048450000044
Figure BDA0001787048450000045
其中,j=1,2,....,8760,且i与j不相等;
(2)将平均差异度最大的场景为第一个聚类中心,并从全年场景集中删掉该场景;
(3)继续寻找剩余场景集中平均差异度最大的场景,作为待选聚类中心,计算其与已选取的聚类中心的距离,并将该场景从全年场景集中删除;
(4)如果步骤(3)中待选为聚类中心的场景与各个已选的聚类中心的距离均大于全年场景集的总体平均差异度,则此场景可作为聚类中心,否则返回到步骤(3);
(5)如果聚类中心数达到K,则将K个聚类中心存入聚类中心集合Tc={ωk}(k=1,2,…,K),否则返回到步骤(3);
(6)在得到聚类中心集合Tc后,接下来确定剩余场景的集合Tm={ωs}(s=1,2,…,8760-K),然后分别计算剩余各场景到各个聚类中心的距离,形成距离矩阵D:
D=D(ωsk)=||ωsk||;
(7)根据距离的大小将剩余各场景重新归类到距离自己最近的聚类中心,此时可以得到K个同类的场景聚类集合Ck(k=1,2,…,K);
(8)依次计算K个同类的场景聚类集合Ck中的每个场景与剩余场景的距离之和,然后以最小距离之和值对应的场景ωk为新的聚类中心;
(9)重复步骤(6)-(8),直到K个聚类中心不再变化为止,结束场景缩减过程;此时聚类后的每个场景概率为同类场景的概率之和。
所述外层规划模型,其目标函数表示为:
Figure BDA0001787048450000051
其中,pi为场景i发生的概率;N为安装光伏电源的位置总数;rs、rz、rb分别表示为火电机组脱硫电价、用户购电电价以及国家补贴电价;Psj,i、Pzj,i分别表示为i场景下位置j处上网的光伏电量和供给用户的光伏电量;Sj为光伏电源在位置j处安装的额定容量值;C0为单位安装容量成本;γ0为折现率;n为光伏电源使用寿命;μ0为运行维护率;
外层规划模型的约束条件为分布式光伏电源安装容量约束为:
Sj,min≤Sj≤Sj,max
其中,Sj,min、Sj,max分别表示为位置j处所安装光伏电源容量的下限和上限。
所述内层规划模型,在外层光伏电源容量确定的基础上,考虑配电网运行约束条件,以光伏年消纳率最大化为目标对光伏电源上网有功出力进行优化削减,由此确定各光伏电源的最优出力值,并传递回外层规划模型用于计算光伏电源年净收益目标值。
所述内层规划模型的目标函数表示为:
Figure BDA0001787048450000061
其中,Pxj,i表示为i场景下位置j处的光伏电源上网的有功出力削减量;PGj,i为i场景下位置j处的光伏电源的有功功率,采用如下计算方法计算PGj,i
Figure BDA0001787048450000062
其中,Ei为i场景下的光照强度,Es为额定光照强度;Sj为光伏电源在位置j处安装的额定容量值。
所述内层规划模型的约束条件表示为:
(1)配电网系统潮流方程等式约束:
Figure BDA0001787048450000063
其中,Pj,i、Qj,i分别为i场景下位置j处注入的有功出力和无功出力;Uj,i为i场景下位置j处的电压值;Uk,i分别为i场景下位置k处的电压值;NL为配电网节点总数;Gj,k表示为配电网位置j和位置k之间的电导值;Bj,k表示为配电网位置j和位置k之间的电纳值;δj,k为位置j和位置k之间的电压相角差。
(2)配电网系统电压约束:
Uj,min≤Uj,i≤Uj,max
其中,Uj,min、Uj,max分别为位置j处电压允许值的下限和上限。
(3)配电网导线支路传输电流约束:
Il,i≤Il,max,l∈Ω
其中,Il,i为i场景下第l条支路流过的电流;Ω为配电网系统线路的集合;Il,max为第l条支路允许流过电流的最大值。
(4)光伏电源上网的有功出力削减量约束:
0≤Pxj,i≤PGj,i-PLj,i
其中,PLj,i分别为i场景下位置j处的有功负荷。
(5)配电网反向潮流约束:
Figure BDA0001787048450000071
其中,PD,i、QD,i分别表示为i场景下位置j之后的有功功率和无功功率;QGj,i表示为i场景下位置j处的光伏电源的无功功率;Qxj,i表示为i场景下位置j处的光伏电源的无功出力削减量;QLj,i分别为i场景下位置j处的无功负荷。
一种分布式光伏规划系统,包括:
数据获取模块,用于获取配电网系统参数。
计算模块,与所述数据获取模块连接,将所述电网系统参数代入预先建立的双层规划模型采用分子微分进化算法进行求解,以求解结果作为规划方案。
其中,所述双层规划模型包括外层规划模型和内层规划模型,所述外层规划模型为以投资者对分布式光伏电源的年净收益最大为目标,所述内层规划模型以分布式光伏电源的年消纳率最大为目标。
本发明的有益效果是,本发明在对光伏电源的安装容量进行规划时,充分考虑光伏的消纳率,可以减少实际运行中光伏电源的弃光限电问题,显著提高投资者的年净收益。本发明采用改进K-均值聚类法对全年的场景进行缩减,有效改善了传统K-均值聚类法存在随机选取初始聚类中心而导致聚类结果不精确的缺点,提高了聚类场景的拟合精度,减少了双层规划模型的计算量。在求解算法上,本发明运用了分子微分进化算法来求解双层多场景规划模型,与DEA/best和DEA/rand相比,MDEA既具有收敛速度快,局部搜索能力强的优点,同时又能提高种群的多样性,增强其全局寻优能力,兼顾了种群个体寻优的速度和深度。
本发明一种分布式光伏规划方法及系统可为配电网中分布式光伏电源的合理规划提供科学指导。
附图说明
图1为本发明分布式光伏规划方法的流程框图。
具体实施方式
本发明的具体实施方式如图1所示。
本实施例一种分布式光伏规划方法实施如以下计算步骤:
1、采用改进K-均值聚类法对全年的场景进行缩减,采用如下步骤:
(1)假设用ωi(i=1,2,…,8760)表示缩减前全年不同时段的“负荷-光照强度”场景,ωi=[ωi,1i,2],ωi,1、ωi,2分别表示为i时刻的负荷和光照强度值,K则为指定的缩减场景数。分别计算全年场景集中每个场景与剩余场景的距离di,j以及每个场景的平均差异度Mi,ave和全年场景集的总体平均差异度Mave
Figure BDA0001787048450000091
Figure BDA0001787048450000092
Figure BDA0001787048450000093
其中,j=1,2,…,8760,且i与j不相等。
(2)将平均差异度最大的场景为第一个聚类中心,并从全年场景集中删掉该场景。
(3)继续寻找剩余场景集中平均差异度最大的场景,作为待选聚类中心,计算其与已选取的聚类中心的距离,并将该场景从全年场景集中删除。
(4)如果步骤(3)中待选为聚类中心的场景与各个已选的聚类中心的距离均大于全年场景集的总体平均差异度,则此场景可作为聚类中心,否则返回到步骤(3)。
(5)如果聚类中心数达到K,则将K个聚类中心存入聚类中心集合Tc={ωk}(k=1,2,…,K),否则返回到步骤(3)。
(6)在得到聚类中心集合Tc后,接下来确定剩余场景的集合Tm={ωs}(s=1,2,…,8760-K),然后分别计算剩余各场景到各个聚类中心的距离,形成距离矩阵D。
D=D(ωsk)=||ωsk||
(7)根据距离的大小将剩余各场景重新归类到距离自己最近的聚类中心,此时可以得到K个同类的场景聚类集合Ck(k=1,2,…,K)。
(8)依次计算K个同类的场景聚类集合Ck中的每个场景与剩余场景的距离之和,然后以最小距离之和值对应的场景ωk为新的聚类中心。
(9)重复步骤(6)-(8),直到K个聚类中心不再变化为止,结束场景缩减过程;此时聚类后的每个场景概率为同类场景的概率之和。
2、外层规划模型的目标函数表示为:
Figure BDA0001787048450000101
其中,pi为场景i发生的概率;N为安装光伏电源的位置总数;rs、rz、rb分别表示为火电机组脱硫电价、用户购电电价以及国家补贴电价;Psj,i、Pzj,i分别表示为i场景下位置j处上网的光伏电量和供给用户的光伏电量;Sj为光伏电源在位置j处安装的额定容量值;C0为单位安装容量成本;γ0为折现率;n为光伏电源使用寿命;μ0为运行维护率。
外层规划模型的约束条件为分布式光伏电源安装容量约束:
Sj,min≤Sj≤Sj,max
其中,Sj,min、Sj,max分别表示为位置j处所安装光伏电源容量的下限和上限。
3、内层规划模型的目标函数及其约束条件
内层规划模型在外层光伏电源容量确定的基础上,考虑配电网运行约束条件,以光伏年消纳率最大化为目标对光伏电源上网有功出力进行优化削减,由此确定各光伏电源的最优出力值,并传递回外层规划模型用于计算光伏电源年净收益目标值。
内层规划模型的目标函数表示为:
Figure BDA0001787048450000111
Pxj,i表示为i场景下位置j处的光伏电源上网的有功出力削减量;PGj,i为i场景下位置j处的光伏电源的有功功率,采用如下计算方法计算PGj,i
Figure BDA0001787048450000112
其中,Ei为i场景下的光照强度,ES为额定光照强度。
内层规划模型的约束条件表示为:
(1)配电网系统潮流方程等式约束:
Figure BDA0001787048450000113
其中,Pj,i、Qj,i分别为i场景下位置j处注入的有功出力和无功出力;Uj,i、Uk,i分别为i场景下位置j处和位置k处的电压值;NL为配电网节点总数;Gj,k、Bj,k分别表示为配电网位置j和位置k之间的电导和电纳值;δj,k为位置j和位置k之间的电压相角差。
(2)配电网系统电压约束:
Uj,min≤Uj,i≤Uj,max
其中,Uj,min、Uj,max分别为位置j处电压允许值的下限和上限。
(3)配电网导线支路传输电流约束:
Il,i≤Il,max,l∈Ω
其中,Il,i为i场景下第l条支路流过的电流;Ω为配电网系统线路的集合;Il,max为第l条支路允许流过电流的最大值。
(4)光伏电源上网的有功出力削减量约束:
0≤Pxj,i≤PGj,i-PLj,i
其中,PLj,i分别为i场景下位置j处的有功负荷。
(5)配电网反向潮流约束:
Figure BDA0001787048450000121
其中,PD,i、QD,i分别表示为i场景下位置j之后的有功功率和无功功率;QGj,i表示为i场景下位置j处的光伏电源的无功功率;Qxj,i表示为i场景下位置j处的光伏电源的无功出力削减量;QLj,i分别为i场景下位置j处的无功负荷。
4、采用分子微分进化算法对双层规划模型进行求解,其步骤如下:
(1)获取配电网系统线路参数,以及全年的负荷和光照强度数据。采用改进K-均值聚类法对全年的“负荷-光照强度”场景进行缩减,得到缩减后指定数量场景的负荷和光照强度数据及其对应的概率。
(2)设置外层分子微分进化算法的基本参数,其中包括种群的规模NP、维度D,最大迭代次数Gmax、分子间的临界距离r0、变异因子F以及交叉因子CR等,令外层进化代数G=0。对外层规划模型中的种群变量即光伏电源的安装容量进行初始化,初始种群为S=[X1,X2,…,XNP],其中的第i个种群个体表示为Xi=[xi,1,xi,2,…,xi,D],第i个个体各维的数据随机产生。
Figure BDA0001787048450000122
其中,d=1,2,…,D;randi(·)为产生均匀分布的随机整数函数;而L则为分布式光伏电源的单位额定容量;xd,max为第d维数据的最大值,xd,min为第d维数据的最小值。
(3)将外层规划模型中的光伏电源安装容量值传递到内层,设置内层分子微分进化算法的基本参数,参数类型与外层相同。求解内层目标函数的最大值。最后计算出种群个体每个场景对应的消纳量然后传递回外层规划模型。
(4)根据外层光伏电源的初始容量并结合内层中求解的光伏消纳量计算外层规划模型的初始目标函数值。
(5)计算G=G+1,对外层种群进行变异、交叉以及校核操作来更新种群。
(6)种群变异。根据构造的变异机制对每个目标向量Xi,G进行变异,产生变异向量Yi,G+1
Yi,G+1=Xbest,G+F(Xr1,G-Xr2,G+Vr1,r2,G)
其中,Xbest,G为进化到第G代时种群的最优个体向量;Xr1,G、Xr2,G分别为第G代时的第r1、r2个个体向量,且i≠best≠r1≠r2;F为变异因子;Vr1,r2,G为第G代时种群个体之间的作用势能。
Figure BDA0001787048450000131
其中,εd为进化到第G代时种群个体向量之间的维差距;ε0为分子间的临界距离。
(7)种群交叉。将外层种群的变异向量Yi,G+1与目标向量Xi,G的各维分量进行随机重组,产生新向量Zi,G+1
Figure BDA0001787048450000132
其中,R为[0,1]之间的随机数;CR为交叉因子;zi,d为新向量Zi,G+1第i个个体d维数据;xi,d为第i个种群个体d维数据;yi,d为变异向量第i个个体d维数据。
(8)种群校核。变异交叉操作后部分种群个体可能会出现越界的情况,按下式进行校核。
Figure BDA0001787048450000141
(9)种群选择。计算校核后的外层种群目标函数值,此时需要再次进入内层规划模型中进行寻优,方法参考步骤(3),最后按照下式做出贪婪选择。
Figure BDA0001787048450000142
(10)若G<Gmax,则返回到步骤(5),进入下一次循环;否则终止寻优过程,输出外层目标的最优值以及相应的光伏电源最佳规划方案。
本实施例一种分布式光伏规划系统,包括:
数据获取模块,用于获取配电网系统参数;
计算模块,与所述数据获取模块连接,将所述电网系统参数代入预先建立的双层规划模型采用分子微分进化算法进行求解,以求解结果作为规划方案。
其中,所述双层规划模型包括外层规划模型和内层规划模型,所述外层规划模型为以投资者对分布式光伏电源的年净收益最大为目标,所述内层规划模型以分布式光伏电源的年消纳率最大为目标。
本实施例中计算模块采用分子微分进化算法进行求解,包括种群变量初始化、内外层数据转换、种群变异子、种群交叉、种群校核与选择;所述种群变量初始化对外层规划模型中的种群变量即光伏电源的安装容量进行初始化;所述内外层数据转换将外层规划模型中的光伏电源安装容量值传递到内层,计算出种群个体每个场景对应的消纳量然后传递回外层规划模型;所述种群变异对每个目标向量进行变异,产生变异向量;所述种群交叉将外层种群的变异向量目标向量的各维分量进行随机重组,产生新向量;所述种群校核与选择用于变异交叉后的校核和外层种群目标函数值的寻优。内层规划型块通过对光伏电源上网有功出力的优化削减,计算确定各光伏电源的最优出力值;外层规划模型用于计算光伏电源年净收益目标值。
本实施例中数据获取模块获取配电网系统的参数。

Claims (8)

1.一种分布式光伏规划方法,其特征在于,所述方法获取配电网系统参数;将所述电网系统参数代入预先建立的双层规划模型,采用分子微分进化算法进行求解,以求解结果作为规划方案;其中所述双层规划模型包括外层规划模型和内层规划模型,所述外层规划模型为以投资者对分布式光伏电源的年净收益最大为目标,所述内层规划模型以分布式光伏电源的年消纳率最大为目标;
采用分子微分进化算法对所述双层规划模型求解的步骤如下:
(1)取配电网系统线路参数,以及预设周期的负荷和光照强度数据;采用改进K-均值聚类法对全年的“负荷-光照强度”场景进行缩减,得到缩减后指定数量场景的负荷和光照强度数据及其对应的概率;
(2)设置外层分子微分进化算法的基本参数,其中包括种群的规模NP、维度D,最大迭代次数Gmax、分子间的临界距离r0、变异因子F以及交叉因子CR,令外层进化代数G=0;对外层规划模型中的种群变量即光伏电源的安装容量进行初始化,初始种群为
Figure FDA0003143153970000012
其中的第i个种群个体表示为Xi=[xi,1,xi,2,…,xi,D],第i个个体各维的数据随机产生:
Figure FDA0003143153970000011
其中,d=1,2,…,D;randi(·)为产生均匀分布的随机整数函数;而L则为分布式光伏电源的单位额定容量;xd,max为第d维数据的最大值,xd,min为第d维数据的最小值;
(3)将外层规划模型中的光伏电源安装容量值传递到内层,设置内层分子微分进化算法的基本参数,参数类型与外层相同;求解内层规划模型中内层目标函数的最大值;最后计算出种群个体每个场景对应的消纳量然后传递回外层规划模型;
(4)根据外层光伏电源的初始容量并结合内层中求解的光伏消纳量计算外层规划模型的初始目标函数值;
(5)计算G=G+1,对外层种群进行变异、交叉以及校核操作更新种群;
(6)若G<Gmax,则返回进入下一次循环;否则终止寻优过程,输出外层目标的最优值以及相应的光伏电源最佳规划方案。
2.根据权利要求1所述的一种分布式光伏规划方法,其特征在于,所述对外层种群进行变异、交叉以及校核操作更新种群的步骤如下:
(1)种群变异,根据构造的变异机制对每个目标向量Xi,G进行变异,产生变异向量Yi,G+1
Yi,G+1=Xbest,G+F(Xr1,G-Xr2,G+Vr1,r2,G);
其中,Xbest,G为进化到第G代时种群的最优个体向量;Xr1,G、Xr2,G分别为第G代时的第r1、r2个个体向量,且i≠best≠r1≠r2;F为变异因子;Vr1,r2,G为第G代时种群个体之间的作用势能;
Figure FDA0003143153970000021
其中,εr1,r2,d为进化到第G代时种群个体向量之间的维差距;ε0为分子间的临界距离;vr1,r2,d为种群个体在d维度时的作用势能;
(2)种群交叉,将外层种群的变异向量Yi,G+1与目标向量Xi,G的各维分量进行随机重组,产生新向量Zi,G+1
Figure FDA0003143153970000022
其中,R为[0,1]之间的随机数;CR为交叉因子;zi,d为新向量Zi,G+1第i个个体d维数据;xi,d为第i个种群个体d维数据;yi,d为变异向量第i个个体d维数据;
(3)种群校核,变异交叉操作后部分种群个体可能会出现越界的情况,按下式进行校核:
Figure FDA0003143153970000031
(4)种群选择,计算校核后的外层种群目标函数值,此时需要再次进入内层规划模型中进行寻优,方法参考步骤(3),最后按照下式做出贪婪选择:
Figure FDA0003143153970000032
3.根据权利要求1所述的一种分布式光伏规划方法,其特征在于,所述采用改进K-均值聚类法对全年的“负荷-光照强度”场景进行缩减,采用如下步骤:
(1)假设用ωi(i=1,2,…,8760)表示缩减前全年不同时段的“负荷-光照强度”场景,ωi=[ωi,1i,2],ωi,1、ωi,2分别表示为i时刻的负荷和光照强度值,K则为指定的缩减场景数;分别计算全年场景集中每个场景与剩余场景的距离di,j、每个场景的平均差异度Mi,ave和全年场景集的总体平均差异度Mave
Figure FDA0003143153970000033
Figure FDA0003143153970000034
Figure FDA0003143153970000035
其中,j=1,2,....,8760,且i与j不相等;
(2)将平均差异度最大的场景为第一个聚类中心,并从全年场景集中删掉该场景;
(3)继续寻找剩余场景集中平均差异度最大的场景,作为待选聚类中心,计算其与已选取的聚类中心的距离,并将该场景从全年场景集中删除;
(4)如果步骤(3)中待选为聚类中心的场景与各个已选的聚类中心的距离均大于全年场景集的总体平均差异度,则将此场景作为聚类中心,否则返回到步骤(3);
(5)如果聚类中心数达到K,则将K个聚类中心存入聚类中心集合Tc={ωk}(k=1,2,…,K),否则返回到步骤(3);
(6)在得到聚类中心集合Tc后,接下来确定剩余场景的集合Tm={ωs}(s=1,2,…,8760-K),然后分别计算剩余各场景到各个聚类中心的距离,形成距离矩阵D:
D=D(ωsk)=||ωsk||;
(7)根据距离的大小将剩余各场景重新归类到距离自己最近的聚类中心,此时可以得到K个同类的场景聚类集合Ck(k=1,2,…,K);
(8)依次计算K个同类的场景聚类集合Ck中的每个场景与剩余场景的距离之和,然后以最小距离之和值对应的场景ωk为新的聚类中心;
(9)重复步骤(6)-(8),直到K个聚类中心不再变化为止,结束场景缩减过程;此时聚类后的每个场景概率为同类场景的概率之和。
4.根据权利要求1所述的一种分布式光伏规划方法,其特征在于,所述外层规划模型,其目标函数表示为:
Figure FDA0003143153970000041
其中,pi为场景i发生的概率;N为安装光伏电源的位置总数;rs、rz、rb分别表示为火电机组脱硫电价、用户购电电价以及国家补贴电价;Psj,i、Pzj,i分别表示为i场景下位置j处上网的光伏电量和供给用户的光伏电量;Sj为光伏电源在位置j处安装的额定容量值;C0为单位安装容量成本;γ0为折现率;n为光伏电源使用寿命;μ0为运行维护率;
外层规划模型的约束条件为分布式光伏电源安装容量约束为:
Sj,min≤Sj≤Sj,max
其中,Sj,min、Sj,max分别表示为位置j处所安装光伏电源容量的下限和上限。
5.根据权利要求1所述的一种分布式光伏规划方法,其特征在于,所述内层规划模型,在外层光伏电源容量确定的基础上,考虑配电网运行约束条件,以光伏年消纳率最大化为目标对光伏电源上网有功出力进行优化削减,由此确定各光伏电源的最优出力值,并传递回外层规划模型用于计算光伏电源年净收益目标值;
所述内层规划模型的目标函数表示为:
Figure FDA0003143153970000051
其中,Pxj,i表示为i场景下位置j处的光伏电源上网的有功出力削减量;PGj,i为i场景下位置j处的光伏电源的有功功率,采用如下计算方法计算PGj,i
Figure FDA0003143153970000052
其中,Ei为i场景下的光照强度,ES为额定光照强度;Sj为光伏电源在位置j处安装的额定容量值;
所述内层规划模型的约束条件表示为:
(1)配电网系统潮流方程等式约束:
Figure FDA0003143153970000061
其中,Pj,i、Qj,i分别为i场景下位置j处注入的有功出力和无功出力;Uj,i为i场景下位置j处的电压值;Uk,i分别为i场景下位置k处的电压值;NL为配电网节点总数;Gj,k表示为配电网位置j和位置k之间的电导值;Bj,k表示为配电网位置j和位置k之间的电纳值;δj,k为位置j和位置k之间的电压相角差;
(2)配电网系统电压约束:
Uj,min≤Uj,i≤Uj,max
其中,Uj,min、Uj,max分别为位置j处电压允许值的下限和上限;
(3)配电网导线支路传输电流约束:
Il,i≤Il,max,l∈Ω
其中,Il,i为i场景下第l条支路流过的电流;Ω为配电网系统线路的集合;Il,max为第l条支路允许流过电流的最大值;
(4)光伏电源上网的有功出力削减量约束:
0≤Pxj,i≤PGj,i-PLj,i
其中,PLj,i分别为i场景下位置j处的有功负荷;
(5)配电网反向潮流约束:
Figure FDA0003143153970000062
其中,PD,i、QD,i分别表示为i场景下位置j之后的有功功率和无功功率;QGj,i表示为i场景下位置j处的光伏电源的无功功率;Qxj,i表示为i场景下位置j处的光伏电源的无功出力削减量;QLj,i分别为i场景下位置j处的无功负荷。
6.实现如权利要求1~5所述一种分布式光伏规划方法的一种分布式光伏规划系统,其特征在于,所述系统包括:
数据获取模块,用于获取配电网系统参数;
计算模块,与所述数据获取模块连接,将所述电网系统参数代入预先建立的双层规划模型采用分子微分进化算法进行求解,以求解结果作为规划方案;
其中,所述双层规划模型包括外层规划模型和内层规划模型,所述外层规划模型为以投资者对分布式光伏电源的年净收益最大为目标,所述内层规划模型以分布式光伏电源的年消纳率最大为目标。
7.根据权利要求6所述的一种分布式光伏规划系统,其特征在于,所述计算模块采用分子微分进化算法进行求解,包括种群变量初始化、内外层数据转换、种群变异、种群交叉、种群校核与选择;所述种群变量初始化对外层规划模型中的种群变量即光伏电源的安装容量进行初始化;所述内外层数据转换将外层规划模型中的光伏电源安装容量值传递到内层,计算出种群个体每个场景对应的消纳量然后传递回外层规划模型;所述种群变异对每个目标向量进行变异,产生变异向量;所述种群交叉将外层种群的变异向量目标向量的各维分量进行随机重组,产生新向量;所述种群校核与选择用于变异交叉后的校核和外层种群目标函数值的寻优。
8.根据权利要求6所述的一种分布式光伏规划系统,其特征在于,所述内层规划模型通过对光伏电源上网有功出力的优化削减,计算确定各光伏电源的最优出力值;所述外层规划模型用于计算光伏电源年净收益目标值。
CN201811019913.2A 2018-09-03 2018-09-03 一种分布式光伏规划方法及系统 Active CN109193754B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811019913.2A CN109193754B (zh) 2018-09-03 2018-09-03 一种分布式光伏规划方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811019913.2A CN109193754B (zh) 2018-09-03 2018-09-03 一种分布式光伏规划方法及系统

Publications (2)

Publication Number Publication Date
CN109193754A CN109193754A (zh) 2019-01-11
CN109193754B true CN109193754B (zh) 2021-10-19

Family

ID=64917876

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811019913.2A Active CN109193754B (zh) 2018-09-03 2018-09-03 一种分布式光伏规划方法及系统

Country Status (1)

Country Link
CN (1) CN109193754B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110266010B (zh) * 2019-05-30 2023-01-17 武汉大学 直流配电网网架结构与分布式光伏多目标协同优化方法
CN112234615B (zh) * 2020-10-09 2022-07-26 国网青海省电力公司 新能源发展协同效率的计算方法和系统
CN112884276A (zh) * 2021-01-14 2021-06-01 华南理工大学 一种分布式光伏容量多目标优化规划方法
CN113690916A (zh) * 2021-09-24 2021-11-23 广东电网有限责任公司 基于时序生产模拟的直流电网储能配置方法、设备及介质
CN118157231B (zh) * 2024-05-11 2024-07-05 国网浙江省电力有限公司 基于光伏集群划分的提高光伏集群功率预测精度的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107301470A (zh) * 2017-05-24 2017-10-27 天津大学 一种配电网扩展规划与光储选址定容的双层优化方法
CN107834540A (zh) * 2017-10-17 2018-03-23 国网宁夏电力公司固原供电公司 一种基于概率约束的分布式光伏接入容量的确定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107301470A (zh) * 2017-05-24 2017-10-27 天津大学 一种配电网扩展规划与光储选址定容的双层优化方法
CN107834540A (zh) * 2017-10-17 2018-03-23 国网宁夏电力公司固原供电公司 一种基于概率约束的分布式光伏接入容量的确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
主动配电网中分布式光伏发电的优化配置研究;邱碧丹 等;《水电能源科学》;20160831;第34卷(第8期);全文 *
基于需求响应的光伏微网储能系统多目标容量优化配置;周楠 等;《电网技术》;20160630;第40卷(第6期);第1709-1716页 *

Also Published As

Publication number Publication date
CN109193754A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109193754B (zh) 一种分布式光伏规划方法及系统
CN107301470B (zh) 一种配电网扩展规划与光储选址定容的双层优化方法
CN112132427B (zh) 一种考虑用户侧多种资源接入的电网多层规划方法
CN107301472B (zh) 基于场景分析法和电压调节策略的分布式光伏规划方法
CN109871989A (zh) 一种含分布式电源的配电网分层规划方法
WO2017000853A1 (zh) 主动配电网多时间尺度协调优化调度方法和存储介质
CN109325608A (zh) 考虑储能并计及光伏随机性的分布式电源优化配置方法
CN109948849B (zh) 一种计及储能接入的配电网网架规划方法
CN107169631B (zh) 基于改进加权Voronoi图的有源配电网变电站规划方法
WO2018077016A1 (zh) 一种基于价值网络的储能系统调度优化方法
CN107732949A (zh) 一种综合风电全年多季节特性的储能布点定容方法
CN108551176B (zh) 一种结合储能均衡技术的储能电池系统容量配置方法
CN112381375B (zh) 一种基于潮流分配矩阵的电网经济运行域快速生成方法
CN111476475B (zh) 多约束条件下梯级水电站短期优化调度方法
CN115640963A (zh) 一种考虑投资运营模式的海上风电接入系统鲁棒规划方法
CN116760068A (zh) 储能参与新型电力市场的自适应调频辅助服务方法
CN114970986A (zh) 一种基于纳什均衡的分布式电源及储能的协同规划方法
CN115833105A (zh) 一种基于集群划分的配电网规划方法
CN117134409A (zh) 考虑电-氢-热互补的微网系统及其多目标优化配置方法
CN116307562A (zh) 一种轨道交通自洽能源系统的光储规划配置方法及装置
CN110970939B (zh) 一种分布式能源集群优化方法及系统
CN111293687B (zh) 一种基于三维粒子群算法的分布式电源选址定容的方法
CN116995740A (zh) 一种配电网分布式风电及储能优化配置方法和系统
CN115459320B (zh) 一种多点分布式储能系统聚合控制智能决策方法及装置
CN116865270A (zh) 一种含嵌入式直流的柔性互联配电网优化调度方法及系统

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