CN105119280A - 基于锥优化的交直流混合结构有源配电网运行优化方法 - Google Patents
基于锥优化的交直流混合结构有源配电网运行优化方法 Download PDFInfo
- Publication number
- CN105119280A CN105119280A CN201510546775.3A CN201510546775A CN105119280A CN 105119280 A CN105119280 A CN 105119280A CN 201510546775 A CN201510546775 A CN 201510546775A CN 105119280 A CN105119280 A CN 105119280A
- Authority
- CN
- China
- Prior art keywords
- prime
- node
- voltage source
- active power
- source converter
- 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
- Supply And Distribution Of Alternating Current (AREA)
- Direct Current Feeding And Distribution (AREA)
Abstract
一种基于锥优化的交直流混合结构有源配电网运行优化方法:根据选定的交直流混合结构有源配电网;依据交直流混合结构有源配电网,同时考虑交流网络损耗、直流网络损耗及电压源换流器的运行损耗,建立交直流混合结构有源配电网运行优化问题的数学模型;将数学模型中目标函数和约束条件的绝对值项进行线性化;根据锥优化的标准形式对绝对值项线性化后的目标函数和约束条件进行锥模型转化;对转化得到的线性化的目标函数、线性约束、二阶锥约束和二阶旋转锥约束的约束条件,采用CLPEX求解器进行求解,并输出结果。本发明可以对交直流混合结构有源配电网和电压源换流器运行优化问题进行统一描述,避免了繁琐的迭代和大量的测试,在计算速度上有大幅地提升。
Description
技术领域
本发明涉及有源配电网运行优化方法。特别是涉及基于锥优化的交直流混合结构有源配电网运行优化方法。
背景技术
新能源和可再生能源的利用在用户侧主要是通过分布式的方式广泛、高密度地接入配电网中,在提供电能供应的同时,也带动了电力电子技术和装置在配电网中的大量配置和使用。分布式发电和先进电力电子技术的广泛应用促进了直流配电网的发展,主要由于:1)大量接入的分布式电源很多为直流型电源或可经过简单整流后变为直流电源;2)随着电力电子技术的发展,电压源换流器(VSC)等电力电子变换装置在用户侧得到了广泛的应用,直流电源和直流负载在配电网中不断增加,如果能在电源侧或负荷侧采用直流供电的方式将大大减少电力电子变流装置的投资。因此,含直流环节的配电网日益引起人们的关注。
未来的智能配电网将很可能是一个含有多种直流环节的交直流全面互联的复杂系统。该系统在整体上采用效率更高的交直流混合的供电形式,从而更好地与上层交直流混合输电体系相协调,并满足配网层分布式电源和用户侧直流负荷的接入。与交流配电网相比,交直流混合结构的配电网可以将大容量可控AC/DC变换器作为直流电源,采用适当的控制策略,有效提高直流系统的供电连续性和电能质量。因为直流环节只存在电阻损耗,整个供电系统的损耗将有望大幅度降低,而且便于分布式电源、储能装置等接入。
在交直流混合结构有源配电网中可以利用的资源包括可控型分布式电源、可控负载、储能装置等。其中,将部分分布式电源由“不可调度”变为“可调度”,并与配电网的传统优化手段相互配合,不仅可以实现分布式电源的充分利用,还可以通过优化电压源换流器VSC传输的有功功率,同时考虑电压源换流器VSC的有功损耗,来降低系统阻塞,提高系统的经济性和安全性。在交直流混合结构有源配电网中,不仅能够优化系统的有功功率,而且可以进行无功功率的优化,从而改善系统潮流分布,提高系统供电能力。在传统的无功优化问题中,用来进行控制的手段主要有发电机无功出力的调节、变压器分接头的调整和无功补偿装置的投切等。而交直流混合结构有源配电网可以将分布式电源和电压源换流器VSC作为连续可调的无功电源参与到接入配电网的电压控制中,为配电网提供无功辅助服务,可解决传统配电网无功调压手段调节速度慢、难以实现电压连续调节的缺陷,并能够减少大容量无功补偿装置的投入。
对于交直流混合结构有源配电网的运行优化问题,其数学本质是大规模非线性规划问题,对于这类决策变量为连续量的非线性数学优化问题,主要的求解方法包括内点法、解析法、连续消去法等。其中,锥优化(ConicProgramming,CP)方法是线性规划与非线性规划的推广,因凸锥所具有的优美的几何结构和特殊的处理方式,能够实现连续优化问题的快速收敛和准确求解,与其他常见算法相比,锥优化方法大大提升了计算速度,并且能够保证所得解的最优性。
发明内容
本发明所要解决的技术问题是,提供一种能够综合考虑交直流网络损耗和电压源换流器的运行损耗,确定合理的电压源换流器的基于锥优化的交直流混合结构有源配电网运行优化方法。
本发明所采用的技术方案是:基于锥优化的交直流混合结构有源配电网运行优化方法,包括如下步骤:
1)根据选定的交直流混合结构有源配电网,输入系统基本参数信息,包括:线路参数、负荷水平、网络拓扑连接关系,系统运行电压水平和各支路电流限制,分布式电源接入位置、类型和容量,电压源换流器接入位置与容量,系统基准电压和基准功率;
2)依据步骤1)给出的交直流混合结构有源配电网,同时考虑交流网络损耗、直流网络损耗及电压源换流器的运行损耗,建立交直流混合结构有源配电网运行优化问题的数学模型,包括:选取根节点为平衡节点,设定有源配电网运行的总损耗最小为目标函数,分别考虑交直流混合结构有源配电网的潮流约束、运行约束及电压源换流器的运行约束和电压源换流器控制方式约束;
3)将步骤2)所建立的交直流混合结构有源配电网运行优化问题的数学模型中目标函数和约束条件的绝对值项进行线性化;根据锥优化的标准形式对绝对值项线性化后的目标函数和约束条件进行锥模型转化;
4)对步骤3)转化得到的线性化的目标函数,以及转化为线性约束、二阶锥约束和二阶旋转锥约束的约束条件,采用CLPEX求解器进行求解;
5)输出步骤4)的求解结果,包括电压源换流器传输的最优功率值、网络潮流结果以及目标函数值。
步骤2)所述的有源配电网运行总损耗最小为目标函数表示为
minf=Ploss,ac+Ploss,dc+Ploss,VSC
式中,交流网络损耗Ploss,ac、直流网络损耗Ploss,dc与电压源换流器的运行损耗Ploss,VSC分别用下式表示
式中,NN,ac为系统中交流网络的节点数,NN,dc为系统中直流网络的节点数,NVSC为系统中接入的电压源换流器个数;Ω(i)为交流节点i的相邻节点的集合;Ω(i′)为直流节点i′的相邻节点的集合;交流节点i和交流节点j是交流支路ij两端的节点,直流节点i′和直流节点j′是直流支路i′j′两端的节点;rij,ac为交流支路ij的电阻;ri′j′,dc为直流支路i′j′的电阻;Iij,ac为流过交流支路ij的电流值;Ii′j′,dc为流过直流支路i′j′的电流值;Pm,ac为第m个电压源换流器交流侧输出的有功功率,Am为第m个电压源换流器的有功损耗系数。
步骤2)所述的交直流混合结构有源配电网的潮流约束中交流配电网络的潮流约束表示为
Pi,ac=PDG,i+PVSC_ac,i-PLOAD,i
Qi,ac=QDG,i+QVSC_ac,i-QLOAD,i
式中,Φ(i)为以交流节点i为末端节点的支路首端节点集合,Ψ(i)为以交流节点i为首端节点的支路末端节点集合;Ui,ac为交流节点i的电压值,xij,ac为交流支路ij的电抗;Pji,ac为交流节点j流向交流节点i的有功功率,Qji,ac为交流节点j流向交流节点i的无功功率;Pi,ac为交流节点i上注入的有功功率之和,PDG,i、PVSC_ac,i、PLOAD,i分别为交流节点i上分布式电源注入的有功功率、电压源换流器交流侧输出的有功功率、负荷消耗的有功功率,Qi,ac为交流节点i上注入的无功功率之和,QDG,i、QVSC_ac,i、QLOAD,i分别为交流节点i上分布式电源注入的无功功率、电压源换流器交流侧输出的无功功率、负荷消耗的无功功率。
步骤2)所述的交直流混合结构有源配电网的潮流约束中的直流配电网络的潮流约束表示为
Pi′,dc=PDG,i′+PVSC_dc,i′-PLOAD,i′
式中,Φ(i′)为以直流节点i′为末端节点的支路首端节点集合,Ψ(i′)为以直流节点i′为首端节点的支路末端节点集合;Ui′,dc为直流节点i′的电压值;Pi′j′,dc为直流节点i′流向直流节点j′的功率;Pi′,dc为直流节点i′上注入的功率之和,PDG,i′、PVSC_dc,i′、PLOAD,i′分别为直流节点i′上分布式电源注入的功率、电压源换流器直流侧输出的功率、负荷消耗的功率。
步骤2)所述的电压源换流器的运行约束表示为
Pm,ac+Pm,dc+Am|Pm,ac|=0
-Qm,max≤Qm,ac≤Qm,max
式中,Pm,ac为第m个电压源换流器交流侧输出的有功功率,Pm,dc为第m个电压源换流器直流侧输出的功率,Qm,ac为第m个电压源换流器交流侧输出的无功功率;Sm,max、Qm,max分别为第m个电压源换流器的视在功率和无功功率的输出上限,m为大于等于1的整数。
步骤2)所述的电压源换流器控制方式约束分别表示为:
当电压源换流器采用交流侧有功功率和无功功率控制时,所述的电压源换流器控制方式约束表示为
Pm,ac=Pm,set
Qm,ac=Qm,set
式中,Pm,set为第m个电压源换流器交流侧有功功率的设定值,Qm,set为第m个电压源换流器交流侧无功功率的设定值;
当电压源换流器采用直流侧电压控制和交流侧无功功率控制时,所述的电压源换流器控制方式约束表示为
Um,ac=Um,set
Qm,ac=Qm,set
式中,Um,ac为第m个电压源换流器直流侧节点电压,Um,set为第m个电压源换流器直流侧节点电压的设定值。
步骤3)所述的将目标函数和约束条件中的绝对值项进行线性化是:
引入辅助变量M1=|Pm,ac|=max{Pm,ac,-Pm,ac}来替换目标函数和约束条件中含有绝对值的项|Pm,ac|,并增加如下约束条件:
M1≥0
M1≥Pm,ac
M1≥-Pm,ac。
步骤3)所述的进行锥模型转化包括:
(1)对目标函数和约束条件中含有的二次项和分别对应采用vi,ac、vi′,dc、lij,ac和li′j′,dc进行替换,替换后得到的线性化目标函数如下,
替换后得到的线性化约束条件如下,
∑j∈Φ(i)(Pji,ac-rji,aclji,ac)+Pi,ac=∑k∈Ψ(i)Pik,ac
∑j∈Φ(i)(Qji,ac-xji,aclji,ac)+Qi,ac=∑k∈Ψ(i)Qik,ac
∑j′∈Φ(i′)(Pj′i′,dc-rj′i′,dclj′i′,dc)+Pi′,dc=∑k′∈Ψ(i′)Pi′k′,dc
再将替换后得到的如下约束条件松弛为二阶锥约束,
松弛得到的二阶锥约束为
||[2Pij,ac2Qij,aclij,ac-vi,ac]T||2≤lij,ac+vi,ac
||[2Pi′j′,dcli′j′,dc-vi′,dc]T||2≤li′j′,dc+vi′,dc;
(2)将电压源换流器的运行约束由非线性约束转换为旋转锥约束
本发明的基于锥优化的交直流混合结构有源配电网运行优化方法,立足于解决交直流混合结构有源配电网的运行优化问题,其数学本质是大规模非线性规划问题(NLP)。本发明依据锥优化方法的基本原理,对优化模型的目标函数与约束条件进行了锥转化,将原问题转化为二阶锥规划问题(SOCP),大大降低了求解难度,便于使用求解工具进行求解。本发明所采用的锥优化方法可以对交直流混合结构有源配电网和电压源换流器运行优化问题进行统一描述,使得复杂的非线性规划问题的求解得以实现,避免了繁琐的迭代和大量的测试,在计算速度上有大幅地提升。并且,因为凸锥所具有的优美的几何结构和特殊的处理方式,使其能够保证所求解问题的解的最优性,将其应用到交直流混合结构有源配电网的运行优化问题中,可以快速获得最优的系统运行方案。
附图说明
图1是交直流混合结构的IEEE33节点算例以及分布式电源、电压源换流器接入位置图;
图2是本发明的基于锥优化的交直流混合结构有源配电网运行优化方法流程图。
具体实施方式
下面结合实施例和附图对本发明的基于锥优化的交直流混合结构有源配电网运行优化方法做出详细说明。
本发明的基于锥优化的交直流混合结构有源配电网运行优化方法,用于有源配电系统的运行优化问题研究中,可以采用集成于MATLAB上的MOSEK、CPLEX、GUROBI等求解器进行求解。本发明采用CPLEX求解器,以图1所示的交直流混合结构的IEEE33节点测试系统为实施例。
如图2所示,本发明的基于锥优化的交直流混合结构有源配电网运行优化方法,包括如下步骤:
1)根据选定的交直流混合结构有源配电网,输入系统基本参数信息,包括:线路参数、负荷水平、网络拓扑连接关系,系统运行电压水平和各支路电流限制,分布式电源接入位置、类型和容量,电压源换流器(VSC)接入位置与容量,系统基准电压和基准功率;
对于本实施例,如图1所示,首先输入交直流混合结构的IEEE33节点系统中线路元件的参数,负荷元件的有功功率、无功功率,以及网络拓扑连接关系。其中节点31-33所在区域改为直流供电,记为直流环节1,经过电压源换流器VSC1与交流系统互联,并且该换流器作为直流系统的平衡节点,采用UdcQ控制,电压等级为10.0kV。节点15-18所在区域改为直流供电,记为直流环节2,经过电压源换流器VSC2与交流系统互联,采用PQ控制,电压等级为10.0kV,两个换流器的视在功率均为600kVA,无功功率输出上限均为200kVar,损耗系数为0.01,详细参数见表1~表4;然后设定节点33接入1台光伏系统PV,接入容量为200kVA,节点16接入1台风电机组WT,接入容量为300kVA,节点18接入1台储能系统ESS,接入容量为600kVA,作为直流环节2的平衡节点;最后设置交流系统的基准电压为12.66kV,直流系统的基准电压为10.0kV,系统基准功率为1MVA。
表1交直流混合结构的IEEE33节点算例中交流负荷接入位置及功率
表2交直流混合结构的IEEE33节点算例中直流负荷接入位置及功率
节点编号 | 功率(kW) | 节点编号 | 功率(kW) |
16 | 60 | 32 | 210 |
17 | 60 | 33 | 60 |
18 | 90 |
表3交直流混合结构的IEEE33节点算例中交流线路参数
表4交直流混合结构的IEEE33节点算例中直流线路参数
2)依据步骤1)给出的交直流混合结构有源配电网,同时考虑交流网络损耗、直流网络损耗及电压源换流器的运行损耗,建立交直流混合结构有源配电网运行优化问题的数学模型,包括:选取根节点为平衡节点,本实施例是选取图1中节点1为平衡节点,设定有源配电网运行的总损耗最小为目标函数,分别考虑交直流混合结构有源配电网的潮流约束、运行约束及电压源换流器的运行约束和电压源换流器控制方式约束;其中,
1、所述的有源配电网运行总损耗最小为目标函数表示为
minf=Ploss,ac+Ploss,dc+Ploss,VSC(1)
式中,交流网络损耗Ploss,ac、直流网络损耗Ploss,dc与电压源换流器的运行损耗Ploss,VSC分别用下式表示
式中,NN,ac为系统中交流网络的节点数,NN,dc为系统中直流网络的节点数,NVSC为系统中接入的电压源换流器个数;Ω(i)为交流节点i的相邻节点的集合;Ω(i′)为直流节点i′的相邻节点的集合;交流节点i和交流节点j是交流支路ij两端的节点,直流节点i′和直流节点j′是直流支路i′j′两端的节点;rij,ac为交流支路ij的电阻;ri′j′,dc为直流支路i′j′的电阻;Iij,ac为流过交流支路ij的电流值;Ii′j′,dc为流过直流支路i′j′的电流值;Pm,ac为第m个电压源换流器交流侧输出的有功功率,Am为第m个电压源换流器的有功损耗系数。
2、所述的交直流混合结构有源配电网的潮流约束中的交流配电网络的潮流约束表示为
Pi,ac=PDG,i+PVSC_ac,i-PLOAD,i(7)
Qi,ac=QDG,i+QVSC_ac,i-QLOAD,i(8)
式中,Φ(i)为以交流节点i为末端节点的支路首端节点集合,Ψ(i)为以交流节点i为首端节点的支路末端节点集合;Ui,ac为交流节点i的电压值,xij,ac为交流支路ij的电抗;Pji,ac为交流节点j流向交流节点i的有功功率,Qji,ac为交流节点j流向交流节点i的无功功率;Pi,ac为交流节点i上注入的有功功率之和,PDG,i、PVSC_ac,i、PLOAD,i分别为交流节点i上分布式电源注入的有功功率、电压源换流器交流侧输出的有功功率、负荷消耗的有功功率,Qi,ac为交流节点i上注入的无功功率之和,QDG,i、QVSC_ac,i、QLOAD,i分别为交流节点i上分布式电源注入的无功功率、电压源换流器交流侧输出的无功功率、负荷消耗的无功功率。
3、所述的交直流混合结构有源配电网的潮流约束中的直流配电网络的潮流约束表示为
Pi′,dc=PDG,i′+PVSC_dc,i′-PLOAD,i′(12)
式中,Φ(i′)为以直流节点i′为末端节点的支路首端节点集合,Ψ(i′)为以直流节点i′为首端节点的支路末端节点集合;Ui′,dc为直流节点i′的电压值;Pi′j′,dc为直流节点i′流向直流节点j′的功率;Pi′,dc为直流节点i′上注入的功率之和,PDG,i′、PVSC_dc,i′、PLOAD,i′分别为直流节点i′上分布式电源注入的功率、电压源换流器直流侧输出的功率、负荷消耗的功率。
4、所述的交流配电网络的运行约束表示为
式中,Ui,min_ac和Ui,max_ac分别为交流节点i的最小允许电压值和最大允许电压值;Iij,max_ac为交流支路ij的最大允许电流值。
5、所述的直流配电网络的运行约束表示为
式中,Ui′,min_dc和Ui′,max_dc分别为直流节点i′的最小允许电压值和最大允许电压值;Ii′j′,max_dc为直流支路i′j′的最大允许电流值。
6、所述的电压源换流器的运行约束表示为
Pm,ac+Pm,dc+Am|Pm,ac|=0(19)
-Qm,max≤Qm,ac≤Qm,max(21)
式中,Pm,ac为第m个电压源换流器交流侧输出的有功功率,Pm,dc为第m个电压源换流器直流侧输出的功率,Qm,ac为第m个电压源换流器交流侧输出的无功功率;Sm,max、Qm,max分别为第m个电压源换流器的视在功率和无功功率的输出上限,m为大于等于1的整数。
7、所述的电压源换流器控制方式约束分别表示为
当电压源换流器采用交流侧有功功率和无功功率控制,即PQ控制时,所述的电压源换流器控制方式约束表示为
Pm,ac=Pm,set(22)
Qm,ac=Qm,set(23)
式中,Pm,set为第m个电压源换流器交流侧有功功率的设定值,Qm,set为第m个电压源换流器交流侧无功功率的设定值;
当电压源换流器采用直流侧电压控制和交流侧无功功率控制,即UdcQ控制时,所述的电压源换流器控制方式约束表示为
Um,ac=Um,set(24)
Qm,ac=Qm,set(25)
式中,Um,ac为第m个电压源换流器直流侧节点电压,Um,set为第m个电压源换流器直流侧节点电压的设定值。
建立了由目标函数和约束条件构成的交直流混合结构有源配电网运行优化问题的数学模型,同时考虑了VSC换流器的运行损耗。
3)将步骤2)所建立的交直流混合结构有源配电网运行优化问题的数学模型中目标函数和约束条件的绝对值项进行线性化;根据锥优化的标准形式对绝对值项线性化后的目标函数和约束条件进行锥模型转化;其中,
1、所述的将目标函数公式(4)和约束条件公式(19)中的绝对值项进行线性化是:
引入辅助变量M1=|Pm,ac|=max{Pm,ac,-Pm,ac}来替换目标函数与约束条件中含有绝对值的项|Pm,ac|,并增加如下约束条件:
M1≥0(26)
M1≥Pm,ac(27)
M1≥-Pm,ac(28)。
2、所述的进行锥模型转化包括:
(1)对目标函数公式(2)和(3)与约束条件公式(5)、(6)、(9)~(11)和(13)~(18)中含有的二次项和分别对应采用vi,ac、vi′,dc、lij,ac和li′j′,dc进行替,替换后得到的线性化目标函数如下,
替换后得到的线性化约束条件如下,
∑j∈Φ(i)(Pji,ac-rji,aclji,ac)+Pi,ac=∑k∈Ψ(i)Pik,ac(30)
∑j∈Φ(i)(Qji,ac-xji,aclji,ac)+Qi,ac=∑k∈Ψ(i)Qik,ac(31)
∑j′∈Φ(i′)(Pj′i′,dc-rj′i′,dclj′i′,dc)+Pi′,dc=∑k′∈Ψ(i′)Pi′k′,dc(33)
再将替换后得到的如下约束条件松弛为二阶锥约束,
松弛得到的二阶锥约束为
||[2Pij,ac2Qij,aclij,ac-vi,ac]T||2≤lij,ac+vi,ac(41)
||[2Pi′j′,dcli′j′,dc-vi′,dc]T||2≤li′j′,dc+vi′,dc(42);
(2)将电压源换流器运行约束公式(20)由非线性约束转换为旋转锥约束
4)对步骤3)转化得到的线性化的目标函数,以及转化为线性约束、二阶锥约束和二阶旋转锥约束的约束条件,采用CLPEX求解器进行求解;
5)输出步骤4)的求解结果,包括电压源换流器传输的最优功率值、网络潮流结果以及目标函数值。
本发明基于锥优化方法实现了交直流混合结构有源配电网运行优化问题的求解。
执行优化计算的计算机硬件环境为Intel(R)Core(TM)i5-4570,主频为3.20GHz,内存为4GB;软件环境为Windows7操作系统。
采用交直流混合结构的IEEE33节点算例如图1所示,本实施例对不同控制方式的电压源换流器VSC进行优化,对采用UdcQ控制方式的电压源换流器VSC优化其传输的无功功率,对采用PQ控制方式的电压源换流器VSC同时优化其传输的有功功率和无功功率,电压源换流器VSC传输功率的优化结果见表5,直流配电网络的电压结果见表6,交直流混合结构有源配电网的损耗结果见表7。通过优化电压源换流器VSC的传输功率,可以降低交直流混合结构配电网的总损耗,改善系统的运行电压水平,从而提高电能质量和供电的可靠性。
表5电压源换流器VSC传输功率的优化结果
换流器 | 有功功率(MW) | 无功功率(MVar) |
VSC1 | 0.07090 | 0.20000 |
VSC2 | 0.52359 | 0.20000 |
表6直流配电网电压结果
节点编号 | 电压幅值(p.u.) | 节点编号 | 电压幅值(p.u.) |
15′ | 0.98965 | 31′ | 1.00000 |
16 | 0.99364 | 32 | 0.99978 |
17 | 0.99742 | 33 | 1.00026 |
18 | 1.00000 |
表7交直流混合结构有源配电网损耗情况
交直流混合结构有源配电网的运行优化问题,其数学本质是大规模非线性规划问题。本发明提出的一种基于锥优化的交直流混合结构有源配电网运行优化方法,能够快速、准确的求解此类问题,并能保证解的最优性。针对交直流混合结构有源配电网的运行优化问题,同时采用基于内点法的优化工具包IPOPT进行求解,并对解的最优性和计算性能进行比较,比较结果见表8。
表8不同求解方法计算性能比较
/ | 本发明方法 | IPOPT |
目标函数值(MW) | 0.09112 | 0.09112 |
求解时间(s) | 0.028 | 0.154 |
Claims (7)
1.基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,包括如下步骤:
1)根据选定的交直流混合结构有源配电网,输入系统基本参数信息,包括:线路参数、负荷水平、网络拓扑连接关系,系统运行电压水平和各支路电流限制,分布式电源接入位置、类型和容量,电压源换流器接入位置与容量,系统基准电压和基准功率;
2)依据步骤1)给出的交直流混合结构有源配电网,同时考虑交流网络损耗、直流网络损耗及电压源换流器的运行损耗,建立交直流混合结构有源配电网运行优化问题的数学模型,包括:选取根节点为平衡节点,设定有源配电网运行的总损耗最小为目标函数,分别考虑交直流混合结构有源配电网的潮流约束、运行约束及电压源换流器的运行约束和电压源换流器控制方式约束;
3)将步骤2)所建立的交直流混合结构有源配电网运行优化问题的数学模型中目标函数和约束条件的绝对值项进行线性化;根据锥优化的标准形式对绝对值项线性化后的目标函数和约束条件进行锥模型转化;
4)对步骤3)转化得到的线性化的目标函数,以及转化为线性约束、二阶锥约束和二阶旋转锥约束的约束条件,采用CLPEX求解器进行求解;
5)输出步骤4)的求解结果,包括电压源换流器传输的最优功率值、网络潮流结果以及目标函数值。
2.根据权利要求1所述的基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,步骤2)所述的有源配电网运行总损耗最小为目标函数表示为
minf=Ploss,ac+Ploss,dc+Ploss,VSC
式中,交流网络损耗Ploss,ac、直流网络损耗Ploss,dc与电压源换流器的运行损耗Ploss,VSC分别用下式表示
式中,NN,ac为系统中交流网络的节点数,NN,dc为系统中直流网络的节点数,NVSC为系统中接入的电压源换流器个数;Ω(i)为交流节点i的相邻节点的集合;Ω(i′)为直流节点i′的相邻节点的集合;交流节点i和交流节点j是交流支路ij两端的节点,直流节点i′和直流节点j′是直流支路i′j′两端的节点;rij,ac为交流支路ij的电阻;ri′j′,dc为直流支路i′j′的电阻;Iij,ac为流过交流支路ij的电流值;Ii′j′,dc为流过直流支路i′j′的电流值;Pm,ac为第m个电压源换流器交流侧输出的有功功率,Am为第m个电压源换流器的有功损耗系数。
3.根据权利要求1所述的基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,步骤2)所述的交直流混合结构有源配电网的潮流约束中交流配电网络的潮流约束表示为
Pi,ac=PDG,i+PVSC_ac,i-PLOAD,i
Qi,ac=QDG,i+QVSC_ac,i-QLOAD,i
式中,Φ(i)为以交流节点i为末端节点的支路首端节点集合,Ψ(i)为以交流节点i为首端节点的支路末端节点集合;Ui,ac为交流节点i的电压值,xij,ac为交流支路ij的电抗;Pji,ac为交流节点j流向交流节点i的有功功率,Qji,ac为交流节点j流向交流节点i的无功功率;Pi,ac为交流节点i上注入的有功功率之和,PDG,i、PVSC_ac,i、PLOAD,i分别为交流节点i上分布式电源注入的有功功率、电压源换流器交流侧输出的有功功率、负荷消耗的有功功率,Qi,ac为交流节点i上注入的无功功率之和,QDG,i、QVSC_ac,i、QLOAD,i分别为交流节点i上分布式电源注入的无功功率、电压源换流器交流侧输出的无功功率、负荷消耗的无功功率。
步骤2)所述的交直流混合结构有源配电网的潮流约束中的直流配电网络的潮流约束表示为
Pi′,dc=PDG,i′+PVSC_dc,i′-PLOAD,i′
式中,Φ(i′)为以直流节点i′为末端节点的支路首端节点集合,Ψ(i′)为以直流节点i′为首端节点的支路末端节点集合;Ui′,dc为直流节点i′的电压值;Pi′j′,dc为直流节点i′流向直流节点j′的功率;Pi′,dc为直流节点i′上注入的功率之和,PDG,i′、PVSC_dc,i′、PLOAD,i′分别为直流节点i′上分布式电源注入的功率、电压源换流器直流侧输出的功率、负荷消耗的功率。
4.根据权利要求1所述的基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,步骤2)所述的电压源换流器的运行约束表示为
Pm,ac+Pm,dc+Am|Pm,ac|=0
-Qm,max≤Qm,ac≤Qm,max
式中,Pm,ac为第m个电压源换流器交流侧输出的有功功率,Pm,dc为第m个电压源换流器直流侧输出的功率,Qm,ac为第m个电压源换流器交流侧输出的无功功率;Sm,max、Qm,max分别为第m个电压源换流器的视在功率和无功功率的输出上限,m为大于等于1的整数。
5.根据权利要求1所述的基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,步骤2)所述的电压源换流器控制方式约束分别表示为:
当电压源换流器采用交流侧有功功率和无功功率控制时,所述的电压源换流器控制方式约束表示为
Pm,ac=Pm,set
Qm,ac=Qm,set
式中,Pm,set为第m个电压源换流器交流侧有功功率的设定值,Qm,set为第m个电压源换流器交流侧无功功率的设定值;
当电压源换流器采用直流侧电压控制和交流侧无功功率控制时,所述的电压源换流器控制方式约束表示为
Um,ac=Um,set
Qm,ac=Qm,set
式中,Um,ac为第m个电压源换流器直流侧节点电压,Um,set为第m个电压源换流器直流侧节点电压的设定值。
6.根据权利要求1所述的基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,步骤3)所述的将目标函数和约束条件中的绝对值项进行线性化是:
引入辅助变量M1=|Pm,ac|=max{Pm,ac,-Pm,ac}来替换目标函数和约束条件中含有绝对值的项|Pm,ac|,并增加如下约束条件:
M1≥0
M1≥Pm,ac
M1≥-Pm,ac。
7.根据权利要求1所述的基于锥优化的交直流混合结构有源配电网运行优化方法,其特征在于,步骤3)所述的进行锥模型转化包括:
(1)对目标函数和约束条件中含有的二次项和分别对应采用vi,ac、vi′,dc、lij,ac和li′j′,dc进行替换,替换后得到的线性化目标函数如下,
替换后得到的线性化约束条件如下,
∑j∈Φ(i)(Pji,ac-rji,aclji,ac)+Pi,ac=∑k∈Ψ(i)Pik,ac
∑j∈Φ(i)(Qji,ac-xji,aclji,ac)+Qi,ac=∑k∈Ψ(i)Qik,ac
∑j′∈Φ(i′)(Pj′i′,dc-rj′i′,dclj′i′,dc)+Pi′,dc=∑k′∈Ψ(i′)Pi′k′,dc
再将替换后得到的如下约束条件松弛为二阶锥约束,
松弛得到的二阶锥约束为
‖[2Pij,ac2Qij,aclij,ac-vi,ac]T‖2≤lij,ac+vi,ac
‖[2Pi′j′,dcli′j′,dc-vi′,dc]T‖2≤li′j′,dc+vi′,dc;
(2)将电压源换流器的运行约束由非线性约束转换为旋转锥约束
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510546775.3A CN105119280B (zh) | 2015-08-31 | 2015-08-31 | 基于锥优化的交直流混合结构有源配电网运行优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510546775.3A CN105119280B (zh) | 2015-08-31 | 2015-08-31 | 基于锥优化的交直流混合结构有源配电网运行优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105119280A true CN105119280A (zh) | 2015-12-02 |
CN105119280B CN105119280B (zh) | 2017-08-08 |
Family
ID=54667206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510546775.3A Active CN105119280B (zh) | 2015-08-31 | 2015-08-31 | 基于锥优化的交直流混合结构有源配电网运行优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105119280B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105321003A (zh) * | 2015-12-04 | 2016-02-10 | 东北电力大学 | 一种含vsc-hvdc的交直流系统多目标潮流优化方法 |
CN106130028A (zh) * | 2016-08-16 | 2016-11-16 | 中国电力科学研究院 | 一种基于交直流主动配电网电压无功协调控制方法 |
CN106329523A (zh) * | 2016-11-19 | 2017-01-11 | 中国南方电网有限责任公司电网技术研究中心 | 考虑不确定性的有源配电网智能软开关鲁棒优化建模方法 |
CN106655227A (zh) * | 2017-01-18 | 2017-05-10 | 天津大学 | 一种基于智能软开关的有源配电网馈线负载平衡方法 |
CN106655226A (zh) * | 2016-11-18 | 2017-05-10 | 天津大学 | 基于智能软开关的有源配电网不对称运行优化方法 |
CN106655177A (zh) * | 2017-01-18 | 2017-05-10 | 中国南方电网有限责任公司电网技术研究中心 | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 |
CN106786599A (zh) * | 2017-01-13 | 2017-05-31 | 天津大学 | 交直流混合微电网双向dc‑ac互联装置智能控制方法 |
CN107123988A (zh) * | 2017-05-12 | 2017-09-01 | 南京理工大学 | 一种考虑恢复量不确定的停电电网负荷恢复二阶锥规划方法 |
CN107332252A (zh) * | 2017-07-18 | 2017-11-07 | 国网江西省电力公司电力科学研究院 | 一种计及广义无功源的配电网低电压综合治理方法 |
CN107947153A (zh) * | 2017-11-06 | 2018-04-20 | 国网上海市电力公司 | 一种最小化交直流混合输电系统中有功损耗的方法 |
CN108599157A (zh) * | 2018-05-17 | 2018-09-28 | 东北电力大学 | 一种考虑电价型需求响应的交直流配合电网优化调度方法 |
CN108599170A (zh) * | 2018-05-14 | 2018-09-28 | 国网青海省电力公司电力科学研究院 | 一种适用于交直流系统的潮流获取方法 |
CN108649580A (zh) * | 2018-05-21 | 2018-10-12 | 武汉大学 | 一种基于二阶锥的交直流混合电力系统安全校正方法 |
CN109038551A (zh) * | 2018-07-20 | 2018-12-18 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 交直流并联输电通道功率优化分配的混合整数凸规划方法 |
CN110137970A (zh) * | 2019-05-28 | 2019-08-16 | 华中科技大学 | 一种基于棱锥近似的无松弛潮流获取方法 |
CN110932259A (zh) * | 2019-11-21 | 2020-03-27 | 国网江苏省电力有限公司南通供电分公司 | 考虑多场景模型及需求侧响应策略的配网重构方法 |
CN113097995A (zh) * | 2021-03-17 | 2021-07-09 | 河北工业大学 | 一种考虑多时间尺度的交直流混合配电网优化运行方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537258A (zh) * | 2015-01-12 | 2015-04-22 | 天津大学 | 分布式储能参与有源配电网运行调节的锥优化建模方法 |
-
2015
- 2015-08-31 CN CN201510546775.3A patent/CN105119280B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537258A (zh) * | 2015-01-12 | 2015-04-22 | 天津大学 | 分布式储能参与有源配电网运行调节的锥优化建模方法 |
Non-Patent Citations (1)
Title |
---|
刘一兵等: "基于混合整数二阶锥规划的三相有源配电网无功优化", 《电力系统自动化》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105321003A (zh) * | 2015-12-04 | 2016-02-10 | 东北电力大学 | 一种含vsc-hvdc的交直流系统多目标潮流优化方法 |
CN105321003B (zh) * | 2015-12-04 | 2019-08-06 | 东北电力大学 | 一种含vsc-hvdc的交直流系统多目标潮流优化方法 |
CN106130028A (zh) * | 2016-08-16 | 2016-11-16 | 中国电力科学研究院 | 一种基于交直流主动配电网电压无功协调控制方法 |
CN106655226A (zh) * | 2016-11-18 | 2017-05-10 | 天津大学 | 基于智能软开关的有源配电网不对称运行优化方法 |
CN106655226B (zh) * | 2016-11-18 | 2019-07-26 | 天津大学 | 基于智能软开关的有源配电网不对称运行优化方法 |
CN106329523A (zh) * | 2016-11-19 | 2017-01-11 | 中国南方电网有限责任公司电网技术研究中心 | 考虑不确定性的有源配电网智能软开关鲁棒优化建模方法 |
CN106329523B (zh) * | 2016-11-19 | 2018-09-25 | 中国南方电网有限责任公司电网技术研究中心 | 考虑不确定性的有源配电网智能软开关鲁棒优化建模方法 |
CN106786599A (zh) * | 2017-01-13 | 2017-05-31 | 天津大学 | 交直流混合微电网双向dc‑ac互联装置智能控制方法 |
CN106786599B (zh) * | 2017-01-13 | 2019-02-19 | 天津大学 | 交直流混合微电网双向dc-ac互联装置智能控制方法 |
CN106655177B (zh) * | 2017-01-18 | 2018-12-28 | 中国南方电网有限责任公司电网技术研究中心 | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 |
CN106655177A (zh) * | 2017-01-18 | 2017-05-10 | 中国南方电网有限责任公司电网技术研究中心 | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 |
CN106655227A (zh) * | 2017-01-18 | 2017-05-10 | 天津大学 | 一种基于智能软开关的有源配电网馈线负载平衡方法 |
CN106655227B (zh) * | 2017-01-18 | 2019-02-15 | 天津大学 | 一种基于智能软开关的有源配电网馈线负载平衡方法 |
CN107123988B (zh) * | 2017-05-12 | 2019-09-27 | 南京理工大学 | 一种考虑恢复量不确定的停电电网负荷恢复二阶锥规划方法 |
CN107123988A (zh) * | 2017-05-12 | 2017-09-01 | 南京理工大学 | 一种考虑恢复量不确定的停电电网负荷恢复二阶锥规划方法 |
CN107332252A (zh) * | 2017-07-18 | 2017-11-07 | 国网江西省电力公司电力科学研究院 | 一种计及广义无功源的配电网低电压综合治理方法 |
CN107332252B (zh) * | 2017-07-18 | 2019-12-06 | 国网江西省电力公司电力科学研究院 | 一种计及广义无功源的配电网低电压综合治理方法 |
CN107947153A (zh) * | 2017-11-06 | 2018-04-20 | 国网上海市电力公司 | 一种最小化交直流混合输电系统中有功损耗的方法 |
CN107947153B (zh) * | 2017-11-06 | 2020-11-06 | 国网上海市电力公司 | 一种最小化交直流混合输电系统中有功损耗的方法 |
CN108599170A (zh) * | 2018-05-14 | 2018-09-28 | 国网青海省电力公司电力科学研究院 | 一种适用于交直流系统的潮流获取方法 |
CN108599157A (zh) * | 2018-05-17 | 2018-09-28 | 东北电力大学 | 一种考虑电价型需求响应的交直流配合电网优化调度方法 |
CN108649580A (zh) * | 2018-05-21 | 2018-10-12 | 武汉大学 | 一种基于二阶锥的交直流混合电力系统安全校正方法 |
CN108649580B (zh) * | 2018-05-21 | 2021-09-24 | 武汉大学 | 一种基于二阶锥的交直流混合电力系统安全校正方法 |
CN109038551B (zh) * | 2018-07-20 | 2020-12-08 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 交直流并联输电通道功率优化分配的混合整数凸规划方法 |
CN109038551A (zh) * | 2018-07-20 | 2018-12-18 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 交直流并联输电通道功率优化分配的混合整数凸规划方法 |
CN110137970A (zh) * | 2019-05-28 | 2019-08-16 | 华中科技大学 | 一种基于棱锥近似的无松弛潮流获取方法 |
CN110932259A (zh) * | 2019-11-21 | 2020-03-27 | 国网江苏省电力有限公司南通供电分公司 | 考虑多场景模型及需求侧响应策略的配网重构方法 |
CN113097995A (zh) * | 2021-03-17 | 2021-07-09 | 河北工业大学 | 一种考虑多时间尺度的交直流混合配电网优化运行方法 |
CN113097995B (zh) * | 2021-03-17 | 2023-08-29 | 河北工业大学 | 一种考虑多时间尺度的交直流混合配电网优化运行方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105119280B (zh) | 2017-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105119280A (zh) | 基于锥优化的交直流混合结构有源配电网运行优化方法 | |
CN107800155B (zh) | 考虑智能软开关的有源配电网多时段孤岛运行方法 | |
CN107171341B (zh) | 一种基于分布式计算的输配电网一体化无功优化方法 | |
CN103730908B (zh) | 一种规模化离网型微电网中储能换流器控制方法 | |
CN108134401B (zh) | 交直流混合系统多目标潮流优化及控制方法 | |
CN105140907B (zh) | 直流微网多智能体自适应下垂一致性协调控制方法及装置 | |
CN105740973A (zh) | 基于混合整数锥规划的智能配电网综合电压无功优化方法 | |
CN109149620B (zh) | 一种自储能多端柔直系统控制方法及系统 | |
CN106410835A (zh) | 一种多端柔性直流输电系统的电网仿真实现方法 | |
CN109390962B (zh) | 一种多端柔直系统的不平衡功率自适应优化分配方法 | |
CN102751720B (zh) | 一种Flexible HVDC潮流计算模型及其计算方法 | |
CN104868500A (zh) | 一种适用于微电网逆变器并联运行控制方法 | |
CN104269861A (zh) | 基于柔性环网控制器电磁环网无功功率环流优化控制方法 | |
CN103715701B (zh) | 考虑电容器操作次数限制的主动配电网无功功率控制方法 | |
CN103593711B (zh) | 一种分布式电源优化配置方法 | |
CN106786805B (zh) | 一种计及光伏集群接入的交直潮流断面协调最优控制方法 | |
Aziz et al. | Admm-based distributed optimization of hybrid mtdc-ac grid for determining smooth operation point | |
CN109217297A (zh) | 交直流主动配电网日前调度二阶锥模型 | |
Montilla-DJesus et al. | Optimal reactive power allocation in an offshore wind farms with LCC-HVdc link connection | |
Dawoud et al. | Enhancing the performance of multi-microgrid with high penetration of renewable energy using modified droop control | |
CN104779609A (zh) | 一种用于互联电网的潮流协同控制方法 | |
CN107332290B (zh) | 一种基于直流线路的区域负荷转供方法 | |
CN106953361B (zh) | 交流微电网运行控制方法及装置 | |
CN107579546A (zh) | 基于双馈风电机组风机拓扑结构的海上风电场无功优化配置方法 | |
CN106992522B (zh) | 电压源换流器潮流计算模型及有源配电网潮流计算方法 |
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 |