CN110601198B - 计及谐波和电压不平衡约束的混合微电网优化运行方法 - Google Patents

计及谐波和电压不平衡约束的混合微电网优化运行方法 Download PDF

Info

Publication number
CN110601198B
CN110601198B CN201911046585.XA CN201911046585A CN110601198B CN 110601198 B CN110601198 B CN 110601198B CN 201911046585 A CN201911046585 A CN 201911046585A CN 110601198 B CN110601198 B CN 110601198B
Authority
CN
China
Prior art keywords
harmonic
node
voltage
power
phase
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
CN201911046585.XA
Other languages
English (en)
Other versions
CN110601198A (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.)
Ningbo Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Ningbo Power Supply Co of State Grid Zhejiang 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 Ningbo Power Supply Co of State Grid Zhejiang Electric Power Co Ltd filed Critical Ningbo Power Supply Co of State Grid Zhejiang Electric Power Co Ltd
Priority to CN201911046585.XA priority Critical patent/CN110601198B/zh
Publication of CN110601198A publication Critical patent/CN110601198A/zh
Application granted granted Critical
Publication of CN110601198B publication Critical patent/CN110601198B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01Arrangements for reducing harmonics or ripples
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/26Arrangements for eliminating or reducing asymmetry in polyphase networks
    • 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
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics
    • 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
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/50Arrangements for eliminating or reducing asymmetry in polyphase networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种计及谐波和电压不平衡约束的混合微电网优化运行方法,包括以下步骤:步骤1:建立混合微电网优化运行模型,以运行成本最低为目标函数,并包含以下约束条件:功率平衡约束、节点电压约束、储能系统运行约束、谐波约束以及三相不平衡约束;步骤2:采用迭代式微电网三相潮流计算方法求解混合微电网优化运行模型,从而获得满足目标函数的各节点进行优化运行的基波电压、基波电流、谐波电流与谐波电压;迭代式微电网三相潮流计算方法是在三相基波潮流的基础上,通过谐波电流与基波电压的线性耦合方程来计算微电网的三相谐波潮流。本发明在考虑安全约束的前提下,满足负荷需求的同时实现综合成本的最优化运行。

Description

计及谐波和电压不平衡约束的混合微电网优化运行方法
技术领域
本发明涉及混合微电网优化运行技术领域。
背景技术
微电网作为消纳可再生能源的一种可行的技术手段备受关注,微电网的容量规模一般较小,接入的电压等级一般也不超过35kV。低压微网具有控制简单、运行灵活和成本较低等优点,已经成为微电网发展的主要方向。然而,低压微网中存在大量以电力电子装置为基波单元的单相负荷和单相微源,由此引起的谐波和三相不平衡问题日趋严重。因此,通过微电网的优化运行改善低压微电网中的电能质量具有重要的意义。微电网经济运行是微电网集成控制及能量管理研究中的重要组成部分。
发明内容
针对上述现有技术的不足,本发明提供一种计及谐波和电压不平衡约束的混合微电网优化运行方法,在考虑安全约束的前提下,满足负荷需求的同时实现综合成本的最优化运行。
为了解决上述技术问题,本发明采用了如下的技术方案:一种计及谐波和电压不平衡约束的混合微电网优化运行方法,包括以下步骤:
步骤1:建立混合微电网优化运行模型,以运行成本最低为目标函数,并包含以下约束条件:功率平衡约束、节点电压约束、储能系统运行约束、谐波约束以及三相不平衡约束;
步骤2:采用迭代式微电网三相潮流计算方法求解混合微电网优化运行模型,从而获得满足目标函数的各节点进行优化运行的基波电压、基波电流、谐波电流与谐波电压;迭代式微电网三相潮流计算方法是在三相基波潮流的基础上,通过谐波电流与基波电压的线性耦合方程来计算微电网的三相谐波潮流。
进一步的,迭代式微电网谐波潮流计算方法包括以下步骤:
步骤201:设第λ次迭代过程中节点i的h次谐波有功功率、谐波无功功率分别为
Figure BDA0002254280340000011
Figure BDA0002254280340000012
并且
Figure BDA0002254280340000013
节点i的基波有功功率、基波无功功率分别为
Figure BDA0002254280340000014
其中,Pi、Qi分别表示节点i的总有功功率、总无功功率;当λ=1时,各节点的初始谐波有功功率率与谐波无功功率均设为0;
步骤202:基于前推回代法计算微电网的三相基波潮流,从而得到λ+1次迭代时微电网中各节点的基波电压;其中,节点i的基波电压为
Figure BDA0002254280340000021
步骤203:将各节点的三相基波电压代入谐波电流与基波电压的线性耦合方程,从而得到λ+1次迭代时微电网中各节点的谐波电流;其中,节点i的h次谐波电流为
Figure BDA0002254280340000022
步骤204:根据节点电压方程得到λ+1次迭代时微电网中各节点的h次谐波电压
Figure BDA0002254280340000023
步骤205:计算λ+1次迭代时各节点的谐波有功功率、谐波无功功率、基波有功功率与基波无功功率;其中,λ+1次迭代时节点i的基波有功功率、基波无功功率分别为
Figure BDA0002254280340000024
分别表示节点i的h次谐波有功功率、h次谐波无功功率;
步骤206:判断相邻两次迭代对应的基波有功功率之差以及基波无功功率之差是否满足收敛条件;若是,则停止计算;若否,令λ=λ+2,并回到步骤202。
进一步的,节点谐波电流与基波电压的线性耦合方程如下:
Figure BDA0002254280340000025
式中,E是规定的单位阵;H表示最高谐波次数,H为奇数,h∈{3,5,...,H};
Figure BDA0002254280340000026
表示微电网的h次谐波电流向量:
Figure BDA0002254280340000027
其中,m表示节点总数;
Figure BDA0002254280340000028
表示微电网的基波电压向量:
Figure BDA0002254280340000029
Y1表示基波节点导纳矩阵:
Figure BDA0002254280340000031
Figure BDA0002254280340000032
B中的元素Zh表示谐波节点阻抗矩阵,B中的元素Yh,k表示h次谐波电压与k次谐波电压的耦合矩阵,k∈{3,5,...,H};
Figure BDA0002254280340000033
Figure BDA00022542803400000311
表示节点iφ相的h次谐波电压与k次谐波电压的耦合关系。
进一步的,混合微电网优化运行模型的目标函数如下:
minf=CESS+CGrid+CIL
式中,CESS表示储能系统运行成本:
Figure BDA0002254280340000034
CGrid表示微电网与配电网购售电交易成本:
Figure BDA0002254280340000035
Figure BDA0002254280340000036
CIL表示可中断负荷补偿成本:
Figure BDA0002254280340000037
式中,T为调度时段总数;NESS表示储能系统总数;
Figure BDA0002254280340000038
为储能系统ness的单位电量运行维护成本系数;
Figure BDA0002254280340000039
为储能系统nessφ相的充放电功率,充电时为正,放电时为负;PGrid(t)为微电网与配电网交互功率,配电网向微电网送电时为正,微电网向配电网送电时为负;pbuy(t)和psell(t)分别为t时刻从电网购电和向电网售电的价格;
Figure BDA00022542803400000310
为负荷nIL的可中断时段;NIL为可中断负荷总数;
Figure BDA0002254280340000041
为可中断负荷nIL在调度周期内单位电量中断成本;
Figure BDA0002254280340000042
为0-1二进制变量,表示可中断负荷nIL的开关状态;
Figure BDA0002254280340000043
表示负荷处于连接状态,
Figure BDA0002254280340000044
表示负荷处于中断状态;
Figure BDA0002254280340000045
为可中断负荷nIL中断功率。
进一步的,混合微电网优化运行模型的约束条件如下:
1)功率平衡约束
Figure BDA0002254280340000046
式中,NWT为风电机组总数;
Figure BDA0002254280340000047
为风机nwtφ相的出力;NPV为光伏电源总数;
Figure BDA0002254280340000048
为光伏电源npvφ相的出力;
Figure BDA0002254280340000049
为微电网与电网φ相的交互功率;
Figure BDA00022542803400000410
为可中断负荷nILφ相的功率;
Figure BDA00022542803400000411
为所有固定负荷φ相的功率;
Figure BDA00022542803400000412
为φ相的网损;
2)节点电压约束
Figure BDA00022542803400000413
Figure BDA00022542803400000414
式中,
Figure BDA00022542803400000415
为节点iφ相的电压;U i,φ(t)和
Figure BDA00022542803400000416
分别为节点iφ相的电压上下限;
Figure BDA00022542803400000417
为节点iφ相的h次谐波电压;
3)储能系统运行约束
Figure BDA00022542803400000418
Figure BDA00022542803400000419
Figure BDA00022542803400000420
式中,
Figure BDA00022542803400000421
Figure BDA00022542803400000422
分别为储能系统nessφ相的的充放电功率上下限;
Figure BDA00022542803400000423
为储能系统ness的电池荷电状态;
Figure BDA00022542803400000424
储能系统ness的容量;
Figure BDA00022542803400000425
Figure BDA00022542803400000426
分别为储能系统ness的荷电状态上下限;
4)谐波约束
Figure BDA0002254280340000051
式中,THDi,φ(t)为节点iφ相的电压畸变率;THDlim为电压畸变率的限值;
5)三相不平衡约束
根据基波潮流计算结果,运用对称分量法获得母线f基波电压正序分量
Figure BDA0002254280340000052
和负序分量
Figure BDA0002254280340000053
母线f的负序电压不平衡度
Figure BDA0002254280340000054
Figure BDA0002254280340000055
式中,
Figure BDA0002254280340000056
表示负序电压不平衡度限值。
与现有技术相比,本发明具有以下有益效果:
1、本发明的混合微电网优化运行模型以运行成本最低为目标函数,确立经济运行的目标;功率平衡约束、节点电压约束与储能系统运行约束是安全运行的保障;谐波约束以及三相不平衡约束则是电能质量的保障。
2、本发明对于混合微电网优化运行模型的求解,考虑到了微电网中负荷和电源一般情况下处于三相不平衡状态,因此放弃采用对称模型计算微电网中的潮流分布,而是采用三相潮流计算方法求解。
3、本发明采用前推回代法来计算三相基波潮流,并建立了谐波电流与基波电压的线性耦合方程,从而能够通过三相基波潮流来求解三相谐波潮流。
附图说明
图1是基于前推回代发的微电网三相潮流原理图。
具体实施方式
下面结合附图和优选实施方式对本发明作进一步的详细说明。
一种计及谐波和电压不平衡约束的混合微电网优化运行方法,包括以下步骤:
步骤1:建立混合微电网优化运行模型,以运行成本最低为目标函数,并包含以下约束条件:功率平衡约束、节点电压约束、储能系统运行约束、谐波约束以及三相不平衡约束;
步骤2采用迭代式微电网三相潮流计算方法求解混合微电网优化运行模型,从而获得满足目标函数的各节点进行优化运行的基波电压、基波电流、谐波电流与谐波电压;迭代式微电网三相潮流计算方法是在三相基波潮流的基础上,通过谐波电流与基波电压的线性耦合方程来计算微电网的三相谐波潮流。
下面分别对步骤1与步骤2进行详细说明。
一、建立混合微电网优化运行模型
①目标函数
将一天分为48个时段,根据新能源发电、不可控负荷预测曲线以及可控负荷的设定,利用能量管理系统对可控谐波源和储能电池运行时间或者功率进行调整,使得微电网在满足电能质量约束的前提下以最低成本运行。微电网运行成本主要由风电和光伏发电成本、储能充放电成本、与电网交易成本和可控负荷中断补偿成本组成。由于微网中风电、光伏以最大功率跟踪方式运行,其运行成本属于不可控的常数,因此目标函数中不考虑其运行成本。目标函数可表示为:
minf=CESS+CGrid+CIL
式中,CESS表示储能系统运行成本:
Figure BDA0002254280340000061
CGrid表示微电网与配电网购售电交易成本:
Figure BDA0002254280340000062
Figure BDA0002254280340000063
CIL表示可中断负荷补偿成本:
Figure BDA0002254280340000064
式中:T为调度时段总数,本具体实施方式中T=48;NESS表示储能系统总数;
Figure BDA0002254280340000065
为储能系统ness的单位电量运行维护成本系数;
Figure BDA0002254280340000066
为储能系统nessφ相的充放电功率,充电时为正,放电时为负;PGrid(t)为微电网与配电网交互功率,配电网向微电网送电时为正,微电网向配电网送电时为负;pbuy(t)和psell(t)分别为t时刻从电网购电和向电网售电的价格;
Figure BDA0002254280340000071
为负荷nIL的可中断时段;NIL为可中断负荷总数;
Figure BDA0002254280340000072
为可中断负荷nIL在调度周期内单位电量中断成本;
Figure BDA0002254280340000073
为0-1二进制变量,表示可中断负荷nIL的开关状态;
Figure BDA0002254280340000074
表示负荷处于连接状态,
Figure BDA0002254280340000075
表示负荷处于中断状态;
Figure BDA0002254280340000076
为可中断负荷nIL中断功率。
微电网优化运行模型的约束条件如下:
1)功率平衡约束
Figure BDA0002254280340000077
式中,NWT为风电机组总数;
Figure BDA0002254280340000078
为风机nwtφ相的出力;NPV为光伏电源总数;
Figure BDA0002254280340000079
为光伏电源npvφ相的出力;
Figure BDA00022542803400000710
为微电网与电网φ相的交互功率;
Figure BDA00022542803400000711
为可中断负荷nILφ相的功率;
Figure BDA00022542803400000712
为所有固定负荷φ相的功率;
Figure BDA00022542803400000713
为φ相的网损;
2)节点电压约束
Figure BDA00022542803400000714
Figure BDA00022542803400000715
式中:
Figure BDA00022542803400000716
为节点iφ相的电压;U i,φ(t)和
Figure BDA00022542803400000717
分别为节点iφ相的电压上下限;
Figure BDA00022542803400000718
为节点iφ相的h次谐波电压。
3)储能系统运行约束
Figure BDA00022542803400000719
Figure BDA00022542803400000720
Figure BDA00022542803400000721
式中:
Figure BDA00022542803400000722
Figure BDA00022542803400000723
分别为储能系统nessφ相的的充放电功率上下限;
Figure BDA00022542803400000724
为储能系统ness的电池荷电状态;
Figure BDA00022542803400000725
储能系统ness的容量;
Figure BDA00022542803400000726
Figure BDA00022542803400000727
分别为储能系统ness的荷电状态上下限。
4)谐波约束
Figure BDA0002254280340000081
式中,THDi,φ(t)为节点iφ相的电压畸变率;THDlim为电压畸变率的限值。由于公共连接点的电能质量直接关系到配电变压器的运行安全和上一级电网的电能质量,因此本具体实施方式仅考虑公共连接点的谐波电压畸变率约束。国标规定,公共连接点处的奇次谐波电压畸变率应不超过4%。
5)三相不平衡约束
根据基波潮流计算结果,运用对称分量法获得母线f基波电压正序分量
Figure BDA0002254280340000082
和负序分量
Figure BDA0002254280340000083
低压电网三相不平衡主要指负序电压不平衡,对零序电压不平衡不予考虑。母线f的负序电压不平衡度
Figure BDA0002254280340000084
Figure BDA0002254280340000085
与谐波电压畸变率类似,本具体实施方式仅考虑公共耦合点的电压不平衡约束。国标规定电力系统公共连接点负序电压不平衡度限值
Figure BDA0002254280340000086
为1.3%。
二、采用迭代式微电网谐波潮流计算方法求解混合微电网优化运行模型
首先建立谐波电流和基波电压的线性耦合方程,基于谐波电流和基波电压的线性耦合方程才能根据基波分量求解谐波分量,对于基波分量的求解采用前推回代法。
2.1立谐波电流和基波电压的线性耦合方程
①不控和相控型谐波源模型
不控和相控型整流装置广泛应用于各种用电设备,是一种典型的非线性谐波源。微电网中绝大多数负荷如空调、计算机、节能灯等均可视为该类型的整流设备。基于调制理论,典型的整流装置交流侧各次谐波电流与谐波电压的关系可表示为
Figure BDA0002254280340000087
上述方程可简写为
Figure BDA0002254280340000091
式中,Yi,+和Yi,-分别为谐波源i谐波耦合导纳矩阵,
Figure BDA0002254280340000092
表示谐波源h次谐波电压与k次谐波电流的耦合关系,而
Figure BDA0002254280340000093
表示整流装置h次谐波电压共轭与h次谐波电流的耦合关系,因为谐波电压共轭对谐波电流的影响极小,因此可将上式
Figure BDA0002254280340000094
忽略。
②全控型谐波源模型
由于不控和相控型谐波源模型不适用于分布式电源变流器此类全控型器件,对于此类的谐波源本文采用交叉频率导纳矩阵模型,如下式所示
Figure BDA0002254280340000095
通过逐次叠加谐波电压,求解谐波特性矩阵参数,该方法仅适用于稳态负载。而分布式电源在优化周期内功率水平处于动态变化,不同功率水平下谐波特性矩阵存在明显的差异,因此本文采用复数域最小二乘法求解全控型谐波源从10%额定功率到100%额定功率水平下,注入电流谐波与电压谐波之间的耦合导纳矩阵,然后利用插值函数获得实际功率水平下谐波特性矩阵。
除了不控、相控整流负荷以及全控型谐波源外,网络中还存在线性负荷、公共连接点等非谐波源节点。对于线性负荷,谐波耦合导纳矩阵为对角阵,矩阵参数为各次谐波对应的负荷等值导纳,而剩余其他节点的谐波特性导纳矩阵参数均设为0。负荷节点一般是包含线性负荷、非线性负荷多种负荷类型的综合节点,由于各种负荷是并联关系,因此该负荷节点的谐波特性矩阵即为各类负荷谐波耦合导纳矩阵之和。
③谐波电压和基波电流线性耦合方程
将网络中m个节点谐波特性矩阵方程组合,由此得到下式
Figure BDA0002254280340000096
式中:
Figure BDA0002254280340000097
和Yi分别为节点i的三相谐波电流相量、谐波电压相量和谐波特性导纳矩阵。
对矩阵方程等式两边进行初等变换,将谐波次数相等元素放在一起,即
Figure BDA0002254280340000101
式中:
Figure BDA0002254280340000102
分别是节点iφ相的h次谐波电流和电压;
Figure BDA0002254280340000103
节点iφ相的h次谐波电压与k次谐波电压的耦合关系。
结合谐波节点阻抗矩阵Zh,谐波电压可表示为
Figure BDA0002254280340000104
故式(5)可变换为
Figure BDA0002254280340000105
将基波电流提取出,余下的谐波电流方程为
Figure BDA0002254280340000106
式中:
Figure BDA0002254280340000111
对式(8)等式左右两端进行初等变换,可得
Figure BDA0002254280340000112
由此,获得微电网各节点谐波电流与基波电压的线性耦合关系,基波电流和谐波电压可根据谐波节点电压方程计算得到。
(2.2)迭代式微电网三相潮流计算:
考虑谐波功率时,微电网节点i的功率可用如下的功率方程表示
Figure BDA0002254280340000113
式中,Pi、Qi分别为节点i总有功和无功功率;Pi 1
Figure BDA0002254280340000114
分别节点i的基波有功和无功功率;Pi h
Figure BDA0002254280340000115
分别为节点i的h次谐波有功和无功功率,并且
Figure BDA0002254280340000116
由此可见,微电网节点功率由基波功率和各次谐波功率组成,而谐波功率可由谐波电流和谐波电压求解。结合式(10)可知,求解谐波参数之前,必须首先求解基波分量。
迭代式微电网三相潮流计算步骤如下:
步骤201:设第λ次迭代过程中节点i的h次谐波有功功率、谐波无功功率分别为
Figure BDA0002254280340000117
Figure BDA0002254280340000118
并且
Figure BDA0002254280340000119
节点i的基波有功功率、基波无功功率分别为
Figure BDA00022542803400001110
其中,Pi、Qi分别表示节点i的总有功功率、总无功功率;当λ=1时,各节点的初始谐波有功功率率与谐波无功功率均设为0;
步骤202:基于前推回代法计算微电网的三相基波潮流,从而得到λ+1次迭代时微电网中各节点的基波电压;其中,节点i的基波电压为
Figure BDA00022542803400001111
步骤203:将各节点的三相基波电压代入谐波电流与基波电压的线性耦合方程,从而得到λ+1次迭代时微电网中各节点的谐波电流;其中,节点i的h次谐波电流为
Figure BDA00022542803400001112
步骤204:根据节点电压方程得到λ+1次迭代时微电网中各节点的h次谐波电压
Figure BDA0002254280340000121
步骤205:计算λ+1次迭代时各节点的谐波有功功率、谐波无功功率、基波有功功率与基波无功功率;其中,λ+1次迭代时节点i的基波有功功率、基波无功功率分别为
Figure BDA0002254280340000122
Figure BDA0002254280340000126
分别表示节点i的h次谐波有功功率、h次谐波无功功率;
步骤206:判断相邻两次迭代对应的基波有功功率之差以及基波无功功率之差是否满足收敛条件;若是,则停止计算;若否,令λ=λ+2,并回到步骤202。
(2.3)基于前推回代法计算微电网的三相基波潮流
微电网中负荷和电源一般情况下处于三相不平衡状态,因此不能采用对称模型计算微电网中的潮流分布,有必要采用三相潮流计算方法求解。微电网通常为辐射型网络,而前推回代法在计算辐射型网络潮流具有计算效率高、收敛效果好等优点。
步骤2021:确定节点类型:由于微电网并网时,配电网可视为无穷大电源,因而将配电网等值电动势作为平衡节点电压,也是潮流计算的根节点。分布式电源节点根据DG模型已确定,而负荷节点统一视为PQ节点。
步骤2022:回代过程,更新支路电流:从最后一条支路的节点开始向根节点回代。将该节点的负荷电流与该节点下游支路电流相加,即得该节点上游支路电流。
图1为辐射型微电网潮流计算原理图,假设母线f+1的φ相(φ=A,B,C)负载为Sf+1,φ。对于母线f和f+1间的φ相支路电流
Figure BDA0002254280340000123
可表示为
Figure BDA0002254280340000124
步骤2023:前推过程,更新节点电压:从根节点开始向最后一层推进,然后继续向下一节点推进。处于母线f+1的三相电压与相邻母线f的三相电压关系可表示为
Figure BDA0002254280340000125
式中:Z为母线f和f+1间支路lf的阻抗矩阵。
一般情况下,相邻母线间的阻抗矩阵为线路阻抗矩阵,而根节点与公共连接点之间的阻抗矩阵为配电网系统等值阻抗矩阵。
步骤2024:判断是否满足收敛条件:计算相邻两次迭代节点电压差的模是否满足收敛条件,如不符合,重复步骤2022与步骤2023,直到满足要求为止。

Claims (3)

1.一种计及谐波和电压不平衡约束的混合微电网优化运行方法,其特征在于,包括以下步骤:
步骤1:建立混合微电网优化运行模型,以运行成本最低为目标函数,并包含以下约束条件:功率平衡约束、节点电压约束、储能系统运行约束、谐波约束以及三相不平衡约束;
步骤2:采用迭代式微电网三相潮流计算方法求解混合微电网优化运行模型,从而获得满足目标函数的各节点进行优化运行的基波电压、基波电流、谐波电流与谐波电压;迭代式微电网三相潮流计算方法是在三相基波潮流的基础上,通过谐波电流与基波电压的线性耦合方程来计算微电网的三相谐波潮流;
迭代式微电网谐波潮流计算方法包括以下步骤:
步骤201:设第λ次迭代过程中节点i的h次谐波有功功率、谐波无功功率分别为
Figure FDA00029054149400000111
Figure FDA0002905414940000011
并且
Figure FDA0002905414940000012
节点i的基波有功功率、基波无功功率分别为
Figure FDA0002905414940000013
其中,Pi、Qi分别表示节点i的总有功功率、总无功功率;当λ=1时,各节点的初始谐波有功功率与谐波无功功率均设为0;H表示最高谐波次数,H为奇数,h∈{3,5,...,H};
Figure FDA0002905414940000014
表示节点i的h次谐波电压,
Figure FDA0002905414940000015
表示节点i的h次谐波电流;
步骤202:基于前推回代法计算微电网的三相基波潮流,从而得到λ+1次迭代时微电网中各节点的基波电压;其中,节点i的基波电压为
Figure FDA0002905414940000016
步骤203:将各节点的三相基波电压代入谐波电流与基波电压的线性耦合方程,从而得到λ+1次迭代时微电网中各节点的谐波电流;其中,节点i的h次谐波电流为
Figure FDA0002905414940000017
步骤204:根据节点电压方程得到λ+1次迭代时微电网中各节点的h次谐波电压
Figure FDA0002905414940000018
步骤205:计算λ+1次迭代时各节点的谐波有功功率、谐波无功功率、基波有功功率与基波无功功率;其中,λ+1次迭代时节点i的基波有功功率、基波无功功率分别为
Figure FDA0002905414940000019
Figure FDA00029054149400000110
分别表示节点i的h次谐波有功功率、h次谐波无功功率;
步骤206:判断相邻两次迭代对应的基波有功功率之差以及基波无功功率之差是否满足收敛条件;若是,则停止计算;若否,令λ=λ+2,并回到步骤202;
节点谐波电流与基波电压的线性耦合方程如下:
Figure FDA0002905414940000021
式中,E是规定的单位阵;H表示最高谐波次数,H为奇数,h∈{3,5,...,H};
Figure FDA0002905414940000022
表示微电网的h次谐波电流向量:
Figure FDA0002905414940000023
其中,m表示节点总数;
Figure FDA0002905414940000024
表示微电网的基波电压向量:
Figure FDA0002905414940000025
Y1表示基波节点导纳矩阵:
Figure FDA0002905414940000026
Figure FDA0002905414940000027
B中的元素Zh表示谐波节点阻抗矩阵,B中的元素Yh,k表示h次谐波电压与k次谐波电压的耦合矩阵,k∈{3,5,...,H};
Figure FDA0002905414940000028
Figure FDA0002905414940000029
表示节点i的φ相的h次谐波电压与k次谐波电压的耦合关系;
基于前推回代法计算微电网的三相基波潮流,按如下步骤:
步骤2021:确定节点类型:将配电网等值电动势作为平衡节点电压,负荷节点统一视为PQ节点,分布式电源节点作为潮流计算的根节点;
步骤2022:回代过程,更新支路电流:从最后一条支路的节点开始向根节点回代:将该节点的负荷电流与该节点下游支路电流相加,即得该节点上游支路电流;
步骤2023:前推过程,更新节点电压:从根节点开始向最后一层推进,然后继续向下一节点推进;
步骤2024:计算相邻两次迭代节点电压差的模是否满足收敛条件;若是,则停止计算;若否,则重复步骤2022与步骤2023,直到满足收敛条件。
2.根据权利要求1所述的计及谐波和电压不平衡约束的混合微电网优化运行方法,其特征在于,混合微电网优化运行模型的目标函数如下:
minf=CESS+CGrid+CIL
式中,CESS表示储能系统运行成本:
Figure FDA0002905414940000031
CGrid表示微电网与配电网购售电交易成本:
Figure FDA0002905414940000032
Figure FDA0002905414940000033
CIL表示可中断负荷补偿成本:
Figure FDA0002905414940000034
式中,T为调度时段总数;NESS表示储能系统总数;
Figure FDA0002905414940000035
为储能系统ness的单位电量运行维护成本系数;
Figure FDA0002905414940000036
为储能系统ness的φ相的充放电功率,充电时为正,放电时为负;PGrid(t)为微电网与配电网交互功率,配电网向微电网送电时为正,微电网向配电网送电时为负;pbuy(t)和psell(t)分别为t时刻从电网购电和向电网售电的价格;
Figure FDA0002905414940000037
为负荷nIL的可中断时段;NIL为可中断负荷总数;
Figure FDA0002905414940000038
为可中断负荷nIL在调度周期内单位电量中断成本;
Figure FDA0002905414940000041
为0-1二进制变量,表示可中断负荷nIL的开关状态;
Figure FDA0002905414940000042
表示负荷处于连接状态,
Figure FDA0002905414940000043
表示负荷处于中断状态;
Figure FDA0002905414940000044
为可中断负荷nIL中断功率。
3.根据权利要求1所述的计及谐波和电压不平衡约束的混合微电网优化运行方法,其特征在于,混合微电网优化运行模型的约束条件如下:
1)功率平衡约束
Figure FDA0002905414940000045
式中,NWT为风电机组总数;NESS表示储能系统总数;
Figure FDA0002905414940000046
为风机nwt的φ相的出力;NPV为光伏电源总数;
Figure FDA0002905414940000047
为光伏电源npv的φ相的出力;
Figure FDA0002905414940000048
为微电网与电网的φ相的交互功率;
Figure FDA0002905414940000049
为可中断负荷nIL的φ相的功率;
Figure FDA00029054149400000410
为所有固定负荷的φ相的功率;
Figure FDA00029054149400000411
为φ相的网损;
2)节点电压约束
Figure FDA00029054149400000412
Figure FDA00029054149400000413
式中,
Figure FDA00029054149400000414
为节点i的φ相的电压;U i,φ(t)和
Figure FDA00029054149400000415
分别为节点i的φ相的电压下限、上限;
Figure FDA00029054149400000416
为节点i的φ相的h次谐波电压;
3)储能系统运行约束
Figure FDA00029054149400000417
Figure FDA00029054149400000418
Figure FDA00029054149400000419
式中,
Figure FDA00029054149400000420
Figure FDA00029054149400000421
分别为储能系统ness的φ相的充放电功率下限、上限;
Figure FDA00029054149400000422
为储能系统ness的电池荷电状态;
Figure FDA00029054149400000423
储能系统ness的容量;
Figure FDA00029054149400000424
Figure FDA00029054149400000425
分别为储能系统ness的荷电状态下限、上限;
4)谐波约束
Figure FDA0002905414940000051
式中,THDi,φ(t)为节点i的φ相的电压畸变率;THDlim为电压畸变率的限值;
5)三相不平衡约束
根据基波潮流计算结果,运用对称分量法获得母线f基波电压正序分量
Figure FDA0002905414940000052
和负序分量
Figure FDA0002905414940000053
母线f的负序电压不平衡度
Figure FDA0002905414940000054
Figure FDA0002905414940000055
式中,
Figure FDA0002905414940000056
表示负序电压不平衡度限值。
CN201911046585.XA 2019-10-30 2019-10-30 计及谐波和电压不平衡约束的混合微电网优化运行方法 Active CN110601198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911046585.XA CN110601198B (zh) 2019-10-30 2019-10-30 计及谐波和电压不平衡约束的混合微电网优化运行方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911046585.XA CN110601198B (zh) 2019-10-30 2019-10-30 计及谐波和电压不平衡约束的混合微电网优化运行方法

Publications (2)

Publication Number Publication Date
CN110601198A CN110601198A (zh) 2019-12-20
CN110601198B true CN110601198B (zh) 2021-02-26

Family

ID=68851932

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911046585.XA Active CN110601198B (zh) 2019-10-30 2019-10-30 计及谐波和电压不平衡约束的混合微电网优化运行方法

Country Status (1)

Country Link
CN (1) CN110601198B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112109583B (zh) * 2020-09-11 2021-09-10 国网山东省电力公司临沂供电公司 一种电网谐波畸变处理方法
CN112667962B (zh) * 2020-11-11 2023-02-14 福州大学 一种电网谐波潮流相量矩阵计算方法
CN112444675B (zh) * 2020-12-14 2022-06-14 南方电网科学研究院有限责任公司 一种输电网节点的谐波超标预警方法、装置、设备及介质
CN112670987B (zh) * 2020-12-30 2023-06-23 国网福建省电力有限公司 一种电网三相谐波潮流相量矩阵计算方法
CN112834820B (zh) * 2021-04-09 2024-01-23 杭州万高科技股份有限公司 一种电能表及其计量装置
CN113852104B (zh) * 2021-09-22 2023-10-10 国网山东省电力公司电力科学研究院 一种三相不对称配电网潮流计算方法
CN113964831B (zh) * 2021-10-29 2023-11-24 国网北京市电力公司 负载调节方法、装置、电子设备及计算机可读存储介质

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102611099A (zh) * 2012-02-09 2012-07-25 河海大学 一种微电网降低网损的方法
CN103972893B (zh) * 2014-04-18 2016-08-31 国家电网公司 一种含有分布式电源的配电网滤波器的优化配置方法
CN103956735B (zh) * 2014-05-12 2016-08-17 河海大学 一种分布式发电系统的谐波潮流分析方法
US10404064B2 (en) * 2015-08-18 2019-09-03 Virginia Tech Intellectual Properties, Inc. Modular multilevel converter capacitor voltage ripple reduction
CN106532710B (zh) * 2016-11-04 2019-04-09 国网山东省电力公司 计及电压稳定约束的微电网潮流优化方法
CN109309385A (zh) * 2018-12-05 2019-02-05 中南大学 一种有源配电网中混合有源滤波器优化配置方法
CN109728602A (zh) * 2018-12-21 2019-05-07 燕山大学 一种基于多功能并网逆变器容量分配的微电网谐波治理方法
CN109995034B (zh) * 2019-01-28 2021-01-19 西安交通大学 一种考虑谐波耦合特性的区域性电网谐波综合治理的方法
CN110071515A (zh) * 2019-05-14 2019-07-30 三峡大学 一种基于二级控制的多微网电压不平衡补偿方法

Also Published As

Publication number Publication date
CN110601198A (zh) 2019-12-20

Similar Documents

Publication Publication Date Title
CN110601198B (zh) 计及谐波和电压不平衡约束的混合微电网优化运行方法
Salazar et al. Energy management of islanded nanogrids through nonlinear optimization using stochastic dynamic programming
CN105914758A (zh) 考虑分布式电源无功补偿成本的电网无功优化方法
Karagiannopoulos et al. A centralised control method for tackling unbalances in active distribution grids
CN105761161A (zh) 一种交直流配电网供电模式评价方法
Mohamed et al. Microgrid online management and balancing using multiobjective optimization
CN110601224B (zh) 基于家庭温控负荷的低压台区电压三相不平衡调节方法
Zaree et al. An MILP formulation for centralized energy management strategy of microgrids
CN114709848A (zh) 一种考虑电能质量补偿的多相配电网储能优化调度方法
CN110994697B (zh) 含光储联合体的交直流配电网优化运行控制方法及系统
Xiao et al. Optimal power quality compensation of energy storage system in distribution networks based on unified multi-phase OPF model
Alanazi et al. Optimal energy storage sizing and siting in hybrid AC/DC microgrids
Garces et al. Optimal operation of distributed energy storage units for minimizing energy losses
Meliala et al. Application of off-grid solar panels system for household electricity consumptions in facing electric energy crisis
Capuder et al. Multi-energy microgrid ability to provide flexibility services to the system operator and security of supply to end-users
CN116404665A (zh) 基于柔性软开关的主动配电网电压三相不平衡抑制方法
Saini et al. Optimal siting and sizing of battery in varying PV generation by utilizing genetic algorithm in distribution system
Zhang et al. Multi-Resource Collaborative Service Restoration of a Distribution Network with Decentralized Hierarchical Droop Control
CN111614113B (zh) 一种在配电网储能系统优化运行的混合整数线性规划模型
He et al. Topology Evolution of AC-DC Distribution Network
Long et al. Multi-objective optimization method for distributed energy storage configuration under distribution network operation constraints
Lujano-Rojas et al. Massive Integration of Wind Power at Distribution Level Supported by Battery Energy Storage Systems
CN108336731A (zh) 一种配电网分布式电源准入容量的计算方法
CN117353304B (zh) 一种用于微电网群的合作运行全分布式控制方法及系统
Yang et al. Research on user energy storage optimization system considering a grid-on new energy

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