发明内容
本发明的目的在于克服现有技术的不足,提出考虑风光无功调节的配电网分布式电源双层规划方法,更有利于源网协同绿色经济运行。在增加DG的控制策略后,可以有效降低支路的网络损耗,减少网损费用;此外,新能源的接入量增加,减少了对上级电网购电需求,降低了年综合成本。
本发明解决其技术问题是采取以下技术方案实现的:
考虑风光无功调节的配电网分布式电源双层规划方法,包括以下步骤:
步骤1、构建以经济性作为主要目标函数,对DG进行选址定容的上层规划模型;
步骤2、构建在得到DG选址定容规划结构下,以每个场景下节点电压波动最小为目标函数,计算最优DG控制策略的下层运行模型;
步骤3、根据步骤1的上层规划模型和步骤2的下层运行模型,构建考虑DG无功控制方案的双层规划模型;
步骤4、采用粒子群算法求解步骤3得到的考虑DG无功控制方案的双层规划模型,得到DG的规划方案。
而且,所述步骤1中上层规划模型的目标函数为:年综合成本最小,年综合成包括配电网运行成本、DG年投资费用、DG年运行维护费用、政府补贴收益和电网年网损费用;决策变量为DG接入位置及容量;约束条件为DG容量约束。
而且,所述上层规划模型中目标函数C为:
min C=CBuy+CDG-cons+CDG-oper-CDG-subs+CLoss
其中,CBuy为配电网向上级电网购电成本,CDG-cons和CDG-oper分别为DG折算到每年的建成成本以及每年运行维修成本,CDG-subs为政府对新能源出力补贴收益,CLoss为有功网损费用;配电网向上级电网购电成本CBuy的计算方法为:
其中,Ns表示典型场景数目;ts表示典型场景内的时段数;Psub,t,s为在典型场景s中时刻t 配电网从上级电网的购电功率;ξt为时刻t的购电电价;
DG折算到每年的建成成本CDG-cons的计算方法为:
其中,r为贴现率;n代表规划总年限;Npv和Nwt分别表示安装光伏和风电的数目;cpv和cwt分别为光伏和风电的单位容量投资费用;PPVi和PWTi为在i节点安装的光伏和风电的装机容量; DG每年运行维修成本CDG-oper的计算方法为:
其中,Nb表示负荷节点总数目;ts为每个典型场景的日内时段总数,dt为一年内每个典型场景的天数,PPVi,t,s为在典型场景s中时刻t时节点i处光伏的出力,PWTi,t,s为在典型场景s中时刻t时节点i处风电的出力,τpv和τwt分别为光伏和风电的单位发电运行成本;
政府对新能源出力补贴收益CDG-subs的计算方法为:
其中,cb,pv和cb,wt分别光伏和风电的单位功率补贴费用;
有功网损费用CLoss的计算方法为:
其中,PLoss为典型场景s中时刻t时系统的总有功网损。
而且,所述上层规划模型中DG容量约束包括:待选节点的DG最大安装容量约束:
其中,PPVi,max和PWTi,max分别为节点i处允许安装的光伏和风机的最大容量,PPVi和PWTi为在i 节点安装的光伏和风电的装机容量。
而且,所述步骤2中下层运行模型的目标函数为:节点电压波动偏差;决策变量为DG 协调优化运行策略;约束条件为光伏/风电运行约束、潮流约束、电压幅值约束、热极限约束和反向潮流约束。
而且,所述下层运行模型的目标函数为:
其中,fu为所有场景下各段时刻下节点电压偏差总和,Vi,t,s为在典型场景s中时刻t时节点 i处电压值,Vp.u.为电压额定幅值,Ns表示典型场景数目;ts表示典型场景内的时段数;Nb表示负荷节点总数目;
而且,所述下层运行模型的约束条件中光伏运行约束为:
QPVi,t,s=PPVi,t,s·tan(cos-1(pfi,t,s))
其中,PPVi,t,s为在典型场景s中时刻t时节点i处光伏的出力,PPVi,max为节点i处允许安装的光伏的最大容量,SPVi为节点i处逆变器容量,QPVi,t,s为光伏无功出力值;pfi,t,s为在典型场景s 中时刻t时节点i处的恒功率因数;i∈{1,...,Npv}代表系统中接入光伏数目;
风电运行约束为:
其中,Q
WTi,min和Q
WTi,max分别为无功出力的极限值,V
WTi,s,t为风电接入节点的电压值,
为恒电压模式下的电压恒定值。
和
为节点电压与
之间的偏差;当无功出力未达到极限值时,
和
均为0,此时
为恒电压控制;若当无功出力达到上限时,此时
为非负数,V
WTi,s,t将低于电压设定值
不再保持恒电压控制;
潮流约束:
其中,PG,sub,i和QG,sub,i为场景s时t时段变电站传输的有功功率和无功功率;PLi,t,s和QLi,t,s为场景s时t时段节点i的有功负荷和无功负荷功率;Gij和Bij为节点i和j之间线路的等效电导和电纳,θij,t,s为节点i和j之间的相角差;
电压幅值约束:
Vi min≤Vi,t,s≤Vi max i∈{1,2,...,NB}
其中,Vi min和Vi max为节点i电压的上下限值;
热极限约束包括变电站功率约束和线路容量约束,变电站功率约束为:
0<PG,sub,t,s<Psub,max
其中,Psub,max为变电站允许的最大有功功率;
线路容量约束为:
其中,S
l,t,s为场景s时t时段线路l的视在功率;
为线路l的容量上限;
反向潮流约束:
其中,P
l-inv,t,s为场景s时t时段线路l的逆向有功功率;
为支路的逆向潮流最大值。
而且,所述步骤3的具体实现方法为:通过上层规划模型计算得到的DG接入位置和容量信息用于求解下层运行模型目标函数下的DG协调优化运行策略,并将下层运行模型结果中各时段的DG运行参数以及网损返回上层规划模型,用于准确计算上层规划模型目标函数,并更新上层规划模型,最终得到DG的最优规划方案和最优运行策略。
而且,所述步骤4包括以下步骤:
步骤4.1:输入所规划配电网网络数据,初始粒子数和迭代最大次数;
步骤4.2:根据上层规划模型,输入初始粒子群,并设置迭代次数T=1;
步骤4.3:将上层规划模型中计算得到DG接入位置和容量传递下层运行模型,针对上层规划模型中的每个场景,使用前推回代法进行潮流计算,通过粒子群算法计算出下层运行模型目标函数的最优解,得到DG的最优控制策略;
步骤4.4:将下层运行模型的DG最优控制策略返回上层规划模型;
步骤4.5:结合上层规划模型DG接入位置和容量以及下层规划模型返回的DG控制策略,求解每个粒子的上层目标函数,获得上层规划模型目标函数的最优适应值及最优解;
步骤4.6:设置T=T+1,判断若T大于迭代最大次数,则计算结束,输出结果,否则进行下一步;
步骤4.7:判断是否符合收敛条件,若不满足条件则更新粒子的速度和位置,返回步骤 4.3,否则输出DG的规划方案。
本发明的优点和积极效果是:
1、本发明建立了以经济性和稳定性为目标的含DG控制模型的双层配电网规划模型。在上层中以最小年综合费用为目标对DG位置容量进行决策,下层计算电压偏差最小化时的DG 控制策略,使用粒子群算法对该模型进行分层求解,求得年费用和电压波动最优的DG安装位置和容量。本发明提出了考虑风光无功调节的配电网分布式电源双层规划模型,该模型将分布式电源逆变器的控制模型加入配电网规划模型中,提升了配电网规划模型的准确性和可靠性。在该规划模型中对配电网的年综合费用和电压偏差两个重要指标进行优化,在保证经济性基础上保证配电网用户电能质量。
2、在进行分布式电源选址定容时,DG的无功控制会影响接入位置和容量的选择。本发明在考虑DG的控制策略后,DG的无功功率可以与有功出力协同配合,提高了电网接纳新能源的能力,增加了新能源的接入容量。
3、本发明所提规划模型和方法更有利于源网协同绿色经济运行。在增加DG的控制策略后,可以有效降低支路的网络损耗,减少网损费用;此外,新能源的接入量增加,减少了对上级电网购电需求,降低了年综合成本。
4、本发明规划阶段的DG控制模型对配电网电能质量有较大改善作用。对比不考虑DG 的无功控制得到的规划结果,在本文所得到的规划结果下运行时,节点电压更加平稳,其中相比于传统规划结果,电压波动可降低45%,对比未加DG的基础潮流电压波动可降低67.6%。因此在规划阶段考虑DG的控制策略对于解决新能源波动性给电网电压带来负面影响是十分有必要的。
具体实施方式
以下结合附图对本发明做进一步详述。
由于新能源出力具有较强的不确定性,其高比例的接入电网会导致电压波动、影响了潮流的走向及线路的有功网损,对系统的安全稳定运行造成很大冲击。而随着分布式电源和逆变器控制技术的成熟,光伏和风电等分布式电源也可以作为控制元件为电网提供无功补偿,参与配电网的运行调度,实现源网协同控制。然而,分布式电源的容量和选址会影响到分布式电源运行策略的优化;反之,与配电网相互配合的分布式电源控制策略也会为系统提供无功补偿维持电压稳定,提高规划的合理性和经济型。因此,有必要在规划阶段考虑源网协同控制模型。
考虑风光无功调节的配电网分布式电源双层规划方法,包括以下步骤:
步骤1、构建以经济性作为主要目标函数,对DG进行选址定容的上层规划模型。
上层规划模型的目标函数为:年综合成本最小,年综合成包括配电网运行成本、DG年投资费用、DG年运行维护费用、政府补贴收益和电网年网损费用;决策变量为DG接入位置及容量,属于离散变量优化问题;约束条件为DG容量约束。
上层规划模型中目标函数C为:
min C=CBuy+CDG-cons+CDG-oper-CDG-subs+CLoss
其中,CBuy为配电网向上级电网购电成本,CDG-cons和CDG-oper分别为DG折算到每年的建成成本以及每年运行维修成本,CDG-subs为政府对新能源出力补贴收益,CLoss为有功网损费用;
配电网向上级电网购电成本CBuy的计算方法为:
其中,Ns表示典型场景数目;ts表示典型场景内的时段数;Psub,t,s为在典型场景s中时刻t 配电网从上级电网的购电功率;ξt为时刻t的购电电价;
DG折算到每年的建成成本CDG-cons的计算方法为:
其中,r为贴现率(固定利率);n代表规划总年限;Npv和Nwt分别表示安装光伏和风电的数目;cpv和cwt分别为光伏和风电的单位容量投资费用;PPVi和PWTi为在i节点安装的光伏和风电的装机容量;
DG每年运行维修成本CDG-oper的计算方法为:
其中,Nb表示负荷节点总数目;ts为每个典型场景的日内时段总数,dt为一年内每个典型场景的天数,PPVi,t,s为在典型场景s中时刻t时节点i处光伏的出力,PWTi,t,s为在典型场景s中时刻t时节点i处风电的出力,τpv和τwt分别为光伏和风电的单位发电运行成本;
政府对新能源出力补贴收益CDG-subs的计算方法为:
其中,cb,pv和cb,wt分别光伏和风电的单位功率补贴费用;
有功网损费用CLoss的计算方法为:
其中,PLoss为典型场景s中时刻t时系统的总有功网损。
上层规划模型中DG容量约束包括:待选节点的DG最大安装容量约束:
其中,PPVi,max和PWTi,max分别为节点i处允许安装的光伏和风机的最大容量,PPVi和PWTi为在i 节点安装的光伏和风电的装机容量。
步骤2、构建在得到DG选址定容规划结构下,以每个场景下节点电压波动最小为目标函数,计算最优DG控制策略的下层运行模型。
下层运行模型的目标函数为:节点电压波动偏差;决策变量为DG协调优化运行策略;约束条件为光伏/风电运行约束、潮流约束、电压幅值约束、热极限约束和反向潮流约束。
下层运行模型的目标函数为:
其中,fu为所有场景下各段时刻下节点电压偏差总和,Vi,t,s为在典型场景s中时刻t时节点 i处电压值,Vp.u.为电压额定幅值,Ns表示典型场景数目;ts表示典型场景内的时段数;Nb表示负荷节点总数目;
在并网过程中,光伏通过逆变器接入电网,而逆变器的剩余容量可为系统提供无功支撑,可以减少系统无功补偿设备的投入,但其无功输出能力与系统容量和有功功率值有关。图2 为光伏的有功/无功曲线,其功率的理论范围应在图中所示的半圆内。下层运行模型的约束条件中光伏运行约束为:
QPVi,t,s=PPVi,t,s·tan(cos-1(pfi,t,s))
其中,PPVi,t,s为在典型场景s中时刻t时节点i处光伏的出力,PPVi,max为节点i处允许安装的光伏的最大容量,SPVi为节点i处逆变器容量,通常为额定功率PPVi,max的1.1倍,,QPVi,t,s为光伏无功出力值;此外在本发明中考虑逆变器采用恒功率因数控制,即t时刻的无功出力值由该时段内恒功率因数pfi,t,s和有功值所决定。一般情况下,逆变器功率因数应控制在超前0.9 到滞后0.9范围内。pfi,t,s为在典型场景s中时刻t时节点i处的恒功率因数;i∈{1,...,Npv}代表系统中接入光伏数目;
风电并网时采用可以提供无功支撑的双馈异步风力发电机,可提供动态无功补偿保持母线电压恒定。在下层规划中,目标函数时电压波动,因此当无功充足时,风电通过调节无功来保持母线电压恒定,是恒电压模型;但当无功补偿达到极限值时,无法满足恒电压控制需转化为恒无功控制,此时无功出力维持在极限值。风电运行约束为:
其中,Q
WTi,min和Q
WTi,max分别为无功出力的极限值,V
WTi,s,t为风电接入节点的电压值,
为恒电压模式下的电压恒定值。
和
为节点电压与
之间的偏差;当无功出力未达到极限值时,
和
均为0,此时
为恒电压控制;若当无功出力达到上限时,此时
为非负数,V
WTi,s,t将低于电压设定值
不再保持恒电压控制;
潮流约束:
其中,PG,sub,i和QG,sub,i为场景s时t时段变电站传输的有功功率和无功功率;PLi,t,s和QLi,t,s为场景s时t时段节点i的有功负荷和无功负荷功率;Gij和Bij为节点i和j之间线路的等效电导和电纳,θij,t,s为节点i和j之间的相角差;
电压幅值约束:
Vi min≤Vi,t,s≤Vi max i∈{1,2,...,NB}
其中,Vi min和Vi max为节点i电压的上下限值;
热极限约束包括变电站功率约束和线路容量约束,变电站功率约束为:
0<PG,sub,t,s<Psub,max
其中,Psub,max为变电站允许的最大有功功率,为了避免反向潮流传输最小值设为0;
线路容量约束为:
其中,S
l,t,s为场景s时t时段线路l的视在功率;
为线路l的容量上限;
反向潮流约束:
其中,P
l-inv,t,s为场景s时t时段线路l的逆向有功功率;
为支路的逆向潮流最大值。
步骤3、如图1所示,根据步骤1的上层规划模型和步骤2的下层运行模型,构建考虑DG无功控制方案的双层规划模型。
通过上层规划模型计算得到的DG接入位置和容量信息用于求解下层运行模型目标函数下的DG协调优化运行策略,并将下层运行模型结果中各时段的DG运行参数以及网损返回上层规划模型,用于准确计算上层规划模型目标函数,并更新上层规划模型,最终得到DG的最优规划方案和最优运行策略。
步骤4、采用粒子群算法求解步骤3得到的考虑DG无功控制方案的双层规划模型,得到DG的规划方案。
针对上述考虑风光无功调节的双层规划模型,本文采用粒子群算法进行求解。粒子群算法是一种模拟鸟群运动机制,追随当前最优解不断更新位置寻求全局最优解的算法,具有以下优点:收敛速度快,有多种措施避免求解陷入局部最优;原理相简单,可扩展性强在很多领域被广泛应用;可调参数少,有成熟的参数选择方案。使用粒子群计算时,在每次迭代过程中需要对每个粒子的位置和方向进行更新:
其中,
和
表示在第T次迭代中粒子i的速度和位置,
是第i个粒子在第T次迭代后的最佳位置,
是第T次迭代后的所有粒子中最佳位置。w是惯性系数,为保证算法在初始阶段的全局搜索能力,w取较大值;而在迭代后期,w值可以减少以提高迭代效率。r
1,r
2是随机系数,round()是一个离散化粒子位置的函数。粒子群算法的求解流程如图3所示。
步骤4.1:输入所规划配电网网络数据,初始粒子数和迭代最大次数;
步骤4.2:根据上层规划模型,输入初始粒子群,并设置迭代次数T=1;
步骤4.3:将上层规划模型中计算得到DG接入位置和容量传递下层运行模型,针对上层规划模型中的每个场景,使用前推回代法进行潮流计算,通过粒子群算法计算出下层运行模型目标函数的最优解,得到DG的最优控制策略;
步骤4.4:将下层运行模型的DG最优控制策略返回上层规划模型;
步骤4.5:结合上层规划模型DG接入位置和容量以及下层规划模型返回的DG控制策略,求解每个粒子的上层目标函数,获得上层规划模型目标函数的最优适应值及最优解;
步骤4.6:设置T=T+1,判断若T大于迭代最大次数,则计算结束,输出结果,否则进行下一步;
步骤4.7:判断是否符合收敛条件,若不满足条件则更新粒子的速度和位置,返回步骤 4.3,否则输出DG的规划方案。
根据上述考虑风光无功调节的配电网分布式电源双层规划方法,使用Matlab软件以IEEE 33节点配电网为例进行验证。网络中总有功负荷为3715kW、总无功负荷为2300kVar,电压等级为12.66kV。除节点0以外,其他节点均可接入分布式电源,且每个节点只能安装一种类型分布式电源,待选安装容量步长为0.1MW,技术参数如表1所示。风电模型采用双参数威布尔分布,形状参数k取2.3,尺度参数c取8.92;光伏模型使用Beta概率分布来模拟,参数α取0.85,β取0.85。使用仿射传播聚类方法对一年8760h的风光荷数据和场景进行缩减,场景个数取30。其余规划运行的参数如表2所示。使用粒子群算法求解双层规划中,设定最大粒子群数目为50,最大迭代次数为100。
表1分布式电源参数
表2配电网相关参数
其中,峰时段:10:00-12:00,17:00-20:00;平时段:6:00-9:00,13:00-16:00,21:00-24:00;谷时段:1:00-5:00。
为了验证本文所提出的考虑风光无功控制的双层规划模型效果,本发明设计了两种规划方案进行对比。方案1:考虑DG控制模型的双层优化,使用本文提出的双层粒子群算法求解;方案2:不考虑下层中DG的控制模型,使用简单的DG有功出力上下限约束:光伏出力约束和风电出力约束。使用粒子群算法求解得到的两种方案结果如表3所示。
表3考虑源网协同控制和不考虑源网协同控制的规划结果对比
如表3中所示,相比于不考虑分布式电源控制模型优化的双层规划,增加DG控制后得到规划方案的新能源接入比例明显提升。这是因为增加分布式电源的无功优化调控后,节点的电压可以维持在较稳定的波动范围内,使系统可以接纳更高比例的新能源接入。分布式电源总安装容量增加后,提高了新能源出力,减少向上级电网购电需求,主网购电成本降低;同时DG的协同控制使得潮流分布得到改善,降低网损成本,因此方案一碳排放量更少,符合双碳目标需求。
为了进一步说明方案1对系统各支路网损的改善情况,选择典型场景的12h和24h时,将两种方案下的支路网损与IEEE33系统的基础支路网损进行对比,如图4所示。方案1和方案2时系统网损均小于基础潮流网损,验证了合理的DG接入可以降低系统网损。在运行过程中,方案1的DG无功根据下层规划运行的结果进行调节,如图5为在12h时,方案1 以15min为间隔的六个分布式电源的无功出力情况。在方案1的风光无功调节下相比于方案 2,无论时在12h还是24h,方案1的网损均小于方案2的网损。在方案1的12h时,总的有功网损相比基础潮流网损从202.67kw降低到94.48kw,降低了53.38%,在24h时,有功网损降低到85.39kw,降低了57.86%。而在方案2中,12h总网损为131.29kw,24h总网损为98.53kw,虽然相比于基础潮流网损均有降低但仍高于方案1的网损值。上述证实了方案1中,DG的协同控制更有利于系统网损降低。
为研究双层规划结果对系统电压的影响,如图6所示为方案1和方案2在某一场景中所选时刻的节点电压曲线,图中黑色折线为IEEE33算例在基础潮流下的电压曲线。其中基础潮流时的节点电压偏差较大,多个节点电压值偏低,而接入分布式电源后系统的电压均有了明显抬升。在图6中,相同颜色表示同一时刻时电压波动,其中实线为方案1的电压值,虚线为方案2的电压值,对比同种颜色的实线和虚线曲线,考虑DG优化调控后,节点的电压得到了有效的改善。以24点时为例,相比于基础潮流,总的节点与标幺值的差降低了67.6%,相比于未考虑优化调控时,电压波动降低了45%。这是因为在本文所提的双层规划模型中,光伏和风电的无功出力均会根据有功值的变化而进行调整,因此可以维持电压在较稳定的状态内变化,减少新能源波动对系统电压的影响。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。