CN105552915A - 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 - Google Patents
一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 Download PDFInfo
- Publication number
- CN105552915A CN105552915A CN201610103626.4A CN201610103626A CN105552915A CN 105552915 A CN105552915 A CN 105552915A CN 201610103626 A CN201610103626 A CN 201610103626A CN 105552915 A CN105552915 A CN 105552915A
- Authority
- CN
- China
- Prior art keywords
- direct current
- delta
- power
- constraint
- lag
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/02—Circuit arrangements for ac mains or ac distribution networks using a single network for simultaneous distribution of power at different frequencies; using a single network for simultaneous distribution of ac power and of dc power
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,该方法以系统网损最小和N-1故障后线路功率、节点电压不越限为目标,设计一个交直流电网最优潮流与静态安全校核迭代的计算模型,通过求解电网中发电机组以及柔性直流系统的有功功率及无功功率设定值,在满足N-1故障后不越限的前提下获得交直流混合电网的最优运行状态,以解决含模块化多电平换流器柔性直流交直流混合电网的优化运行调度问题。
Description
技术领域
本发明涉及电力系统优化技术领域,特别是一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法。
背景技术
随着IGBT等全控型电力电子器件的不断发展和日趋成熟,基于模块化多电平换流器(ModularMultilevelConverter,MMC)的柔性高压直流输电(HighVoltageDirectCurrent,HVDC)已成为下一代直流输电的发展方向。模块化多电平换流器由于采用了多个子模块串联构成换流阀,加上最近电平逼近调制策略的应用,使得阀组可以在较低的单元件开关频率下获得很高的等效开关频率,显著改善了输出电压的波形,无需另外装设滤波器。
随着柔性直流输电技术在电力系统的逐步推广应用,含柔性直流的交直流混合系统运行控制技术已得到研究人员的广泛关注。相比传统的相控换流器(PhaseControlConverter,PCC)直流输电,柔性直流输电系统每个换流站可以独立调整有功功率和无功功率,提高了电网控制运行灵活度,在实际运行中,如何在保障系统静态安全的前提下调整电网中发电机组和柔性直流系统的有功及无功设定值,实现交直流混合电网的安全、优化运行,是摆在电网运行和方式人员面前的一道亟待解决的问题。
发明内容
有鉴于此,本发明的目的是提出一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,该方法以系统网损最小和N-1故障后线路功率、节点电压不越限为目标,设计一个交直流电网最优潮流与静态安全校核迭代修正的计算模型,通过求解电网中发电机组以及柔性直流系统的有功功率及无功功率设定值,在满足N-1故障后线路功率、节点电压不越限的前提下获得交直流混合电网的最优运行状态。
本发明采用以下方案实现:一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,具体包括以下步骤:
步骤S1:获取含柔性直流的交直流混合电网相关参数,所述电网相关参数包括交流侧电网各设备参数以及柔性直流换流站的参数;
步骤S2:根据所述电网相关参数,初始化含MMC-HVDC交直流电网最优潮流模型的各项参数:模型中自变量x、约束下限松弛变量l、约束上限松弛变量u赋初值,拉格朗日乘子z、w、y赋初值,设定优化迭代收敛精度10-6;静态安全校核的故障集SC包含交流电网各线路、变压器、发电机、负荷单一开断故障以及柔性直流系统闭锁故障,初始化修正次数k=0;
步骤S3:应用原-对偶内点法求解步骤S2中的含MMC-HVDC交直流电网最优潮流模型,获得电网中发电机组及柔性直流系统有功功率及无功功率最优设定值;
步骤S4:基于所述故障集SC,对步骤S3的计算结果进行静态安全校核,判断最优潮流计算结果是否存在N-1越限;
步骤S5:若步骤S4中静态安全校核的结果不存在N-1越限,则优化计算结果满足静态安全要求,获得满足电网安全要求的最优解,结束计算;若存在N-1越限,则进入步骤S6;
步骤S6:结合静态安全校核结果,得到优化模型约束条件的修正策略,同时修正次数k=k+1;
步骤S7:若k已达到最大次数限制kmax,则无法获得满足静态安全约束的最优解;否则根据步骤S6得到的修正策略,对优化模型的相应约束条件进行修正,并返回步骤S3进行下一次优化计算。
进一步地,在步骤S1中,获取的交流侧电网相关参数还包括:母线节点名称编号、母线电压上下限约束、负荷有功及无功、各节点无功补偿容量、发电机组有功出力以及经济参数、发电机有功无功上下限约束、线路以及主变首末端节点编号、线路及主变等值参数、线路及主变传输功率约束;柔性直流换流站的参数包括:换流变电阻、换流变电抗、MMC换流阀的有功损耗特性系数、MMC换流阀的有功及无功上下限约束、MMC换流阀的电压调制比和直流电压上下限约束、换流站间直流线路电阻。
进一步地,在步骤S2中,最优潮流模型为:
obj.minf(x)
s.t.h(x)=0(1)
gm≤g(x)≤gM
其中x为最优潮流数学模型的优化变量,维度为n,f(x)为最优潮流模型的目标函数,h(x)为等式约束合集,维度为m,g(x)为不等式约束合集,维度为s,gm、gM分别为不等式约束的上下限值;
其中,(1)式中目标函数f(x)的表达式为:
其中PGi为交直流电网中第i台发电机组的有功出力,Ng为交直流电网中发电机组的数目;
其中(1)式中优化变量x具体内容为:
x=[V,θ,Pg,Qg,Pc,Qc,Vdc]T(3)
其中V、θ分别为所有交流节点的电压幅值及相角,包括MMC换流阀交流侧节点;Pg、Qg分别为电网中发电机有功及无功出力;Pc、Qc以及Vdc分别代表MMC换流阀的交流侧注入有功、无功与直流侧节点电压;对一个含p个交流节点、q台发电机组、r个MMC换流站的系统,优化变量x的维度是n=2p+2q+3r;
其中(1)式中等式约束主要包括节点功率平衡方程,由于引入MMC换流阀注入功率作为优化变量,节点功率方程可以针对交流节点和柔性直流节点分别列写,其中交流节点功率平衡方程如下:
其中P、Q、V分别代表交流节点的注入有功功率、注入无功功率及电压,下标gi、di以及ci分别代表交流节点i连接的发电机、负荷以及MMC换流阀;Gij、Bij为交流节点i,j之间线路导纳的实部、虚部,θi、θj为交流节点i、j的相角;
MMC换流阀直流侧节点的有功平衡方程如下:
其中NDC为MMC换流阀个数,VDCj、VDCk为换流阀j、k直流侧节点电压,YDCjk为换流阀j、k之间直流线路电导,Pcj、Qcj、Vcj分别为换流阀j的交流侧注入有功功率、无功功率及交流侧节点电压,K0、K1、K2为换流阀有功损耗系数;
为系统指定一个相角参考点,设定参考点相角为0,即:
θslack=0(7)
(1)式中不等式约束包含交流侧电网及直流侧电网相关约束,其中交流侧电网约束条件包括:节点电压幅值及相角约束,线路/变压器输送功率约束,发电机有功及无功出力约束;直流侧电网约束条件包括:直流母线电压约束,换流阀容量约束,换流阀电压调制比约束,直流线路输送功率约束。
进一步地,所述步骤S3中所述应用原-对偶内点法求解步骤S2中的含MMC-HVDC交直流电网最优潮流模型,步骤如下:
步骤S31:构造式(1)的拉格朗日函数为:
其中,l=[l1,…,ls]T、u=[u1,…,us]T为(1)式优化模型不等式约束的松弛变量,y=[y1,…,ym]T和z=[z1,…,zs]T、w=[w1,…,ws]T是拉格朗日乘子;
根据Karush-Kuhn-Tucker定理,可得到最优解的KKT条件为:
其中L=diag(l1,…,ls),U=diag(u1,…,us),Z=diag(z1,…,zs),W=diag(w1,…,ws),I为单位矩阵,Lagx、Lagy、Lagz、Lagw、分别为式(8)拉格朗日函数关于x、y、l、u、z、w的偏导数;
步骤S32:计算原-对偶互补间隙Gap,并判断Gap是否小于迭代收敛精度ε,是则迭代收敛,输出结果并结束计算,否则转入步骤S33;原-对偶互补间隙Gap的表达式如下:
Gap=lTz-uTw(10)
步骤S33:计算本次迭代的障碍参数μ;障碍参数μ的表达式如下:
步骤S34:求解修正方程组,得到各状态变量及拉格朗日乘子的修正量,所述修正方程组根据(9)式KKT条件得到,表达式为:
上式中:
求解(12)式得到各状态变量及拉格朗日乘子的修正量Δx、Δl、Δu、Δz、Δw、Δy;
步骤S35:计算原始变量和对偶变量的迭代步长,并修正优化数学模型各状态变量及拉格朗日乘子;原始变量的迭代步长stepP为:
对偶变量的迭代步长stepD为:
根据式(15)修正最优潮流数学模型各状态变量及拉格朗日乘子:
进一步的,步骤S4针对线路/变压器单一故障以及发电机、负荷单一故障对步骤S3的最优潮流计算结果进行静态安全校核,所采用的方法是目前已公开并普遍使用的补偿法。
进一步地,所述步骤S6中根据静态安全校核结果得到优化模型约束条件的修正策略的方法如下:已公开的相关文献指出实际物理电网具有同调特性,即在不改变电网拓扑结构及参数的前提下,对故障前某一元件状态变量Xi0(包括线路/变压器有功功率、节点电压)的调整将获得故障后该状态量Xik同一趋势的改变。对某一元件状态变量Xi0以及获得故障后该状态量Xik的关系进行线性化处理:
X′ik-Xik=εik(X′i0-Xi0)(16)
其中Xi、X'i分别代表调整前及调整后的状态变量,下标0以及k分别代表故障前及故障k之后的值,εik为与故障k相应的线性化系数;
设静态安全校核结果中故障k引起的状态变量Xi越限量为Δik,Xi包括线路/变压器输送有功功率及节点电压幅值,以越上限为例,为使该状态变量回到约束范围以内,需有下式成立:
X′ik-Xik≤-Δik(17)
结合(16)、(17)式,可得:
对故障集SC中各故障,均有上式成立,则静态安全校核后,下一次优化计算中Xi的约束上限应满足:
同理,对于静态安全校核结果Xi越下限情况,下一次优化计算中Xi的约束下限应满足:
根据(19)、(20)式,结合静态安全校核结果得到下一次优化计算时约束条件的修正策略。
与现有技术相比,本发明有以下有益效果:本方法以系统网损最小和N-1故障后线路功率、节点电压不越限为目标,设计一个交直流电网最优潮流与静态安全校核迭代修正的计算模型,通过求解电网中发电机组以及柔性直流系统的有功功率及无功功率设定值,在满足N-1故障后线路功率、节点电压不越限的前提下获得交直流混合电网的最优运行状态。该方法弥补了含MMC-HVDC交直流混合电网静态安全约束优化潮流计算方面的不足,具有很强的实用性和广阔的应用前景。
附图说明
图1为本发明的流程示意图。
图2是本发明所涉及的MMC-HVDC换流站模型示意图。
图3为本发明实施例中采用的修改后的IEEE39节点模型图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
如图1所示,本实施例提供了一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,具体包括以下步骤:
步骤S1:获取含柔性直流的交直流混合电网相关参数,所述电网相关参数包括交流侧电网各设备参数以及柔性直流换流站的参数;
步骤S2:根据所述电网相关参数,初始化含MMC-HVDC交直流电网最优潮流模型的各项参数:模型中自变量x、约束下限松弛变量l、约束上限松弛变量u赋初值,拉格朗日乘子z、w、y赋初值,设定优化迭代收敛精度10-6;静态安全校核的故障集SC包含交流电网各线路、变压器、发电机、负荷单一开断故障以及柔性直流系统闭锁故障,初始化修正次数k=0;
步骤S3:应用原-对偶内点法求解步骤S2中的含MMC-HVDC交直流电网最优潮流模型,获得电网中发电机组及柔性直流系统有功功率及无功功率最优设定值;
步骤S4:基于所述故障集SC,对步骤S3的计算结果进行静态安全校核,判断最优潮流计算结果是否存在N-1越限;
步骤S5:若步骤S4中静态安全校核的结果不存在N-1越限,则优化计算结果满足静态安全要求,获得满足电网安全要求的最优解,结束计算;若存在N-1越限,则进入步骤S6;
步骤S6:结合静态安全校核结果,得到优化模型约束条件的修正策略,同时修正次数k=k+1;
步骤S7:若k已达到最大次数限制kmax,则无法获得满足静态安全约束的最优解;否则根据步骤S6得到的修正策略,对优化模型的相应约束条件进行修正,并返回步骤S3进行下一次优化计算。
在本实施例中,在步骤S1中,获取的交流侧电网相关参数还包括:母线节点名称编号、母线电压上下限约束、负荷有功及无功、各节点无功补偿容量、发电机组有功出力以及经济参数、发电机有功无功上下限约束、线路以及主变首末端节点编号、线路及主变等值参数、线路及主变传输功率约束;柔性直流换流站的参数包括:换流变电阻、换流变电抗、MMC换流阀的有功损耗特性系数、MMC换流阀的有功及无功上下限约束、MMC换流阀的电压调制比和直流电压上下限约束、换流站间直流线路电阻。
在本实施例中,在步骤S2中,最优潮流模型为:
obj.minf(x)
s.t.h(x)=0(1)
gm≤g(x)≤gM
其中x为最优潮流数学模型的优化变量,维度为n,f(x)为最优潮流模型的目标函数,h(x)为等式约束合集,维度为m,g(x)为不等式约束合集,维度为s,gm、gM分别为不等式约束的上下限值;
其中,(1)式中目标函数f(x)的表达式为:
其中PGi为交直流电网中第i台发电机组的有功出力,Ng为交直流电网中发电机组的数目;
其中(1)式中优化变量x具体内容为:
x=[V,θ,Pg,Qg,Pc,Qc,Vdc]T(3)
其中V、θ分别为所有交流节点的电压幅值及相角,包括MMC换流阀交流侧节点;Pg、Qg分别为电网中发电机有功及无功出力;Pc、Qc以及Vdc分别代表MMC换流阀的交流侧注入有功、无功与直流侧节点电压;对一个含p个交流节点、q台发电机组、r个MMC换流站的系统,优化变量x的维度是n=2p+2q+3r;
其中(1)式中等式约束主要包括节点功率平衡方程,由于引入MMC换流阀注入功率作为优化变量,节点功率方程可以针对交流节点和柔性直流节点分别列写,其中交流节点功率平衡方程如下:
其中P、Q、V分别代表交流节点的注入有功功率、注入无功功率及电压,下标gi、di以及ci分别代表交流节点i连接的发电机、负荷以及MMC换流阀;Gij、Bij为交流节点i,j之间线路导纳的实部、虚部,θi、θj为交流节点i、j的相角;
MMC换流阀直流侧节点的有功平衡方程如下:
其中NDC为MMC换流阀个数,VDCj、VDCk为换流阀j、k直流侧节点电压,YDCjk为换流阀j、k之间直流线路电导,Pcj、Qcj、Vcj分别为换流阀j的交流侧注入有功功率、无功功率及交流侧节点电压,K0、K1、K2为换流阀有功损耗系数;
为系统指定一个相角参考点,设定参考点相角为0,即:
θslack=0(7)
(2)式中不等式约束包含交流侧电网及直流侧电网相关约束,其中交流侧电网约束条件包括:节点电压幅值及相角约束,线路/变压器输送功率约束,发电机有功及无功出力约束;直流侧电网约束条件包括:直流母线电压约束,换流阀容量约束,换流阀电压调制比约束,直流线路输送功率约束。
在本实施例中,所述步骤S3中所述应用原-对偶内点法求解步骤S2中的含MMC-HVDC交直流电网最优潮流模型,步骤如下:
步骤S31:构造式(1)的拉格朗日函数为:
其中,l=[l1,…,ls]T、u=[u1,…,us]T为(1)式优化模型不等式约束的松弛变量,y=[y1,…,ym]T和z=[z1,…,zs]T、w=[w1,…,ws]T是拉格朗日乘子;
根据Karush-Kuhn-Tucker定理,可得到最优解的KKT条件为:
其中L=diag(l1,…,ls),U=diag(u1,…,us),Z=diag(z1,…,zs),W=diag(w1,…,ws),I为单位矩阵,Lagx、Lagy、Lagz、Lagw、分别为式(8)拉格朗日函数关于x、y、l、u、z、w的偏导数;
步骤S32:计算原-对偶互补间隙Gap,并判断Gap是否小于迭代收敛精度ε,是则迭代收敛,输出结果并结束计算,否则转入步骤S33;原-对偶互补间隙Gap的表达式如下:
Gap=lTz-uTw(10)
步骤S33:计算本次迭代的障碍参数μ;障碍参数μ的表达式如下:
步骤S34:求解修正方程组,得到各状态变量及拉格朗日乘子的修正量,所述修正方程组根据(9)式KKT条件得到,表达式为:
上式中:
求解(12)式得到各状态变量及拉格朗日乘子的修正量Δx、Δl、Δu、Δz、Δw、Δy;
步骤S35:计算原始变量和对偶变量的迭代步长,并修正优化数学模型各状态变量及拉格朗日乘子;原始变量的迭代步长stepP为:
对偶变量的迭代步长stepD为:
根据式(15)修正最优潮流数学模型各状态变量及拉格朗日乘子:
在本实施例中,步骤S4针对线路/变压器单一故障以及发电机、负荷单一故障对步骤S3的最优潮流计算结果进行静态安全校核,所采用的方法是目前已公开并普遍使用的补偿法。
在本实施例中,所述步骤S6中根据静态安全校核结果得到优化模型约束条件的修正策略的方法如下:已公开的相关文献指出实际物理电网具有同调特性,即在不改变电网拓扑结构及参数的前提下,对故障前某一元件状态变量Xi0(包括线路/变压器有功功率、节点电压)的调整将获得故障后该状态量Xik同一趋势的改变。对某一元件状态变量Xi0以及获得故障后该状态量Xik的关系进行线性化处理:
X′ik-Xik=εik(X′i0-Xi0)(16)
其中Xi、X'i分别代表调整前及调整后的状态变量,下标0以及k分别代表故障前及故障k之后的值,εik为与故障k相应的线性化系数;
设静态安全校核结果中故障k引起的状态变量Xi越限量为Δik,Xi包括线路/变压器输送有功功率及节点电压幅值,以越上限为例,为使该状态变量回到约束范围以内,需有下式成立:
X′ik-Xik≤-Δik(17)
结合(16)、(17)式,可得:
对故障集SC中各故障,均有上式成立,则静态安全校核后,下一次优化计算中Xi的约束上限应满足:
同理,对于静态安全校核结果Xi越下限情况,下一次优化计算中Xi的约束下限应满足:
根据(19)、(20)式,结合静态安全校核结果得到下一次优化计算时约束条件的修正策略。
图2为本发明所涉及的MMC-HVDC换流站模型。图2代表本发明所述交直流电网中的一个MMC换流站,其中ACSystem指换流站所连接的交流侧电网,Rtr、Xtr分别为换流变压器的电阻、电抗,Vs、θs、Vc、θc分别为换流变压器交流电网侧节点的电压幅值、相角以及换流阀侧的电压幅值、相角,Ps、Qs、Pc、Qc分别为交流电网流入换流变压器的有功功率、无功功率以及换流变压器流入换流阀的有功功率、无功功率,VDC、IDC分别为换流阀直流侧的电压、流入直流系统的电流,Ploss为换流阀有功损耗。Rf为启动电阻,稳态运行过程中退出,无需考虑。
图3是本发明方法一实施例中采用的修改后的IEEE39节点模型图。图3中,1-39分别为IEEE39节点电网中各个母线编号,为电网中各发电机组,↓为电网中节点上的负荷,MMC1及MMC2分别为电网中组成双端柔性直流输电系统的两个换流站。选取修改后的IEEE39节点测试系统的最优潮流问题对本发明提出的方法进行测试,观察本发明提出的静态安全约束最优潮流模型的计算效果。
修改后的IEEE39节点电网中MMC1及MMC2柔性直流换流站相关参数如表1所示。
表1换流站相关参数
表中Rtr以及Xtr分别为换流变压器的电阻和电抗值;为换流阀有功最大值;Mmax及Mmin分别为电压调制比上下限;K2~K0为换流阀损耗公式系数,根据实际工程设备参数拟合计算得到;及分别为直流母线电压的上下限。表中各参数均已按标幺值处理,直流侧电压基值取VdcBase=320kV,功率基值与交流侧一致,均为100MVA。
本实施例包含如下四个计算方式:
Case1:不考虑静态安全约束,仅考虑正常方式下的节点电压/线路功率约束;
Case2:在Case1的基础上增加N-1故障后的各节点电压约束;
Case3:在Case1的基础上增加N-1故障后的各支路功率约束。
Case4:综合考虑N-1故障后的节点电压和支路功率约束。
节点电压和支路功率约束修正的线性化系数分别取值:εV=2.0,εP=2.5。
本实施例四个计算方式的主要结果如表2所示。
表2修改后的IEEE39节点系统Case1~Case4计算结果
Case2在Case1基础上增加N-1节点电压约束,通过相应调整发电机及换流站无功,对不满足N-1约束节点的电压上下限进行修正,在损失一定最优性的前提下重新获得安全约束可行解。交流节点电压约束取Vac∈[0.9,1.1],包括正常及故障后状态,计算过程对节点电压约束的修正情况如表3所示。
表3Case2节点电压修正情况
Case3考虑N-1支路功率约束,本实施例在Case1结果基础上设定支路功率限额,换流站近区部分线路及部分断面出现N-1越限,计算过程中修正情况如表4。
表4Case3支路功率修正情况
表中各支路功率限额均取301MW,对于N-1越限的支路,在修正其功率约束上限后,下一次OPF计算时,发电机及柔性直流有功相应调整,使得支路或断面功率下降,数次迭代修正后满足静态安全要求。从表4结果还可发现,在第1次支路功率约束修正后,由于柔性直流功率调整,导致支路(14-15)出现了新的N-1越限。
Case4由于综合考虑了支路功率和节点电压N-1约束,迭代修正次数和系统网损进一步增加。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (5)
1.一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,其特征在于:包括以下步骤:
步骤S1:获取含柔性直流的交直流混合电网相关参数,所述电网相关参数包括交流侧电网各设备参数以及柔性直流换流站的参数;
步骤S2:根据所述电网相关参数,初始化含MMC-HVDC交直流电网最优潮流模型的各项参数:模型中自变量x、约束下限松弛变量l、约束上限松弛变量u赋初值,拉格朗日乘子z、w、y赋初值,设定优化迭代收敛精度10-6;静态安全校核的故障集SC包含交流电网各线路、变压器、发电机、负荷单一开断故障以及柔性直流系统闭锁故障,初始化修正次数k=0;
步骤S3:应用原-对偶内点法求解步骤S2中的含MMC-HVDC交直流电网最优潮流模型,获得电网中发电机组及柔性直流系统有功功率及无功功率最优设定值;
步骤S4:基于所述故障集SC,对步骤S3的计算结果进行静态安全校核,判断最优潮流计算结果是否存在N-1越限;
步骤S5:若步骤S4中静态安全校核的结果不存在N-1越限,则优化计算结果满足静态安全要求,获得满足电网安全要求的最优解,结束计算;若存在N-1越限,则进入步骤S6;
步骤S6:结合静态安全校核结果,得到优化模型约束条件的修正策略,同时修正次数k=k+1;
步骤S7:若k已达到最大次数限制kmax,则无法获得满足静态安全约束的最优解;否则根据步骤S6得到的修正策略,对优化模型的相应约束条件进行修正,并返回步骤S3进行下一次优化计算。
2.根据权利要求1所述的一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,其特征在于:在步骤S1中,获取的交流侧电网相关参数还包括:母线节点名称编号、母线电压上下限约束、负荷有功及无功、各节点无功补偿容量、发电机组有功出力以及经济参数、发电机有功无功上下限约束、线路以及主变首末端节点编号、线路及主变等值参数、线路及主变传输功率约束;柔性直流换流站的参数包括:换流变电阻、换流变电抗、MMC换流阀的有功损耗特性系数、MMC换流阀的有功及无功上下限约束、MMC换流阀的电压调制比和直流电压上下限约束、换流站间直流线路电阻。
3.根据权利要求1所述的一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,其特征在于:
在步骤S2中,最优潮流模型为:
obj.minf(x)
s.t.h(x)=0(1)
gm≤g(x)≤gM
其中x为最优潮流数学模型的优化变量,维度为n,f(x)为最优潮流模型的目标函数,h(x)为等式约束合集,维度为m,g(x)为不等式约束合集,维度为s,gm、gM分别为不等式约束的上下限值;
其中,(1)式中目标函数f(x)的表达式为:
其中PGi为交直流电网中第i台发电机组的有功出力,Ng为交直流电网中发电机组的数目;
其中(1)式中优化变量x具体内容为:
x=[V,θ,Pg,Qg,Pc,Qc,Vdc]T(3)
其中V、分别为所有交流节点的电压幅值及相角,包括MMC换流阀交流侧节点;Pg、Qg分别为电网中发电机有功及无功出力;Pc、Qc以及Vdc分别代表MMC换流阀的交流侧注入有功、无功与直流侧节点电压;对一个含p个交流节点、q台发电机组、r个MMC换流站的系统,优化变量x的维度是n=2p+2q+3r;
其中(1)式中等式约束主要包括节点功率平衡方程,由于引入MMC换流阀注入功率作为优化变量,节点功率方程可以针对交流节点和柔性直流节点分别列写,其中交流节点功率平衡方程如下:
其中P、Q、V分别代表交流节点的注入有功功率、注入无功功率及电压,下标gi、di以及ci分别代表交流节点i连接的发电机、负荷以及MMC换流阀;Gij、Bij为交流节点i,j之间线路导纳的实部、虚部,为交流节点i、j的相角;
MMC换流阀直流侧节点的有功平衡方程如下:
其中NDC为MMC换流阀个数,VDCj、VDCk为换流阀j、k直流侧节点电压,YDCjk为换流阀j、k之间直流线路电导,Pcj、Qcj、Vcj分别为换流阀j的交流侧注入有功功率、无功功率及交流侧节点电压,K0、K1、K2为换流阀有功损耗系数;
为系统指定一个相角参考点,设定参考点相角为0,即:
θslack=0(7)
(1)式中不等式约束包含交流侧电网及直流侧电网相关约束,其中交流侧电网约束条件包括:节点电压幅值及相角约束,线路/变压器输送功率约束,发电机有功及无功出力约束;直流侧电网约束条件包括:直流母线电压约束,换流阀容量约束,换流阀电压调制比约束,直流线路输送功率约束。
4.根据权利要求1所述的一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,其特征在于:
所述步骤S3中所述应用原-对偶内点法求解步骤S2中的含MMC-HVDC交直流电网最优潮流模型,步骤如下:
步骤S31:构造式(1)的拉格朗日函数为:
其中,l=[l1,…,ls]T、u=[u1,…,us]T为(1)式优化模型不等式约束的松弛变量,y=[y1,…,ym]T和z=[z1,…,zs]T、w=[w1,…,ws]T是拉格朗日乘子;
根据Karush-Kuhn-Tucker定理,可得到最优解的KKT条件为:
其中L=diag(l1,…,ls),U=diag(u1,…,us),Z=diag(z1,…,zs),W=diag(w1,…,ws),I为单位矩阵,Lagx、Lagy、Lagz、Lagw、分别为式(8)拉格朗日函数关于x、y、l、u、z、w的偏导数;
步骤S32:计算原-对偶互补间隙Gap,并判断Gap是否小于迭代收敛精度ε,是则迭代收敛,输出结果并结束计算,否则转入步骤S33;原-对偶互补间隙Gap的表达式如下:
Gap=lTz-uTw(10)
步骤S33:计算本次迭代的障碍参数μ;障碍参数μ的表达式如下:
步骤S34:求解修正方程组,得到各状态变量及拉格朗日乘子的修正量,所述修正方程组根据(9)式KKT条件得到,表达式为:
上式中:
求解(12)式得到各状态变量及拉格朗日乘子的修正量Δx、Δl、Δu、Δz、Δw、Δy;
步骤S35:计算原始变量和对偶变量的迭代步长,并修正优化数学模型各状态变量及拉格朗日乘子;原始变量的迭代步长stepP为:
对偶变量的迭代步长stepD为:
根据式(15)修正最优潮流数学模型各状态变量及拉格朗日乘子:
5.根据权利要求1所述的一种考虑静态安全约束的含MMC-HVDC交直流混合电网优化潮流计算方法,其特征在于:
所述步骤S6中根据静态安全校核结果得到优化模型约束条件的修正策略的方法如下:对某一元件状态变量Xi0以及获得故障后该状态量Xik的关系进行线性化处理:
X′ik-Xik=εik(X′i0-Xi0)(16)
其中Xi、X'i分别代表调整前及调整后的状态变量,下标0以及k分别代表故障前及故障k之后的值,εik为与故障k相应的线性化系数;
设静态安全校核结果中故障k引起的状态变量Xi越限量为Δik,Xi包括线路/变压器输送有功功率及节点电压幅值,以越上限为例,为使该状态变量回到约束范围以内,需有下式成立:
X′ik-Xik≤-Δik(17)
结合(16)、(17)式,可得:
对故障集SC中各故障,均有上式成立,则静态安全校核后,下一次优化计算中Xi的约束上限应满足:
同理,对于静态安全校核结果Xi越下限情况,下一次优化计算中Xi的约束下限应满足:
根据(19)、(20)式,结合静态安全校核结果得到下一次优化计算时约束条件的修正策略。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610103626.4A CN105552915B (zh) | 2016-02-25 | 2016-02-25 | 一种考虑静态安全约束的含mmc‑hvdc交直流混合电网优化潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610103626.4A CN105552915B (zh) | 2016-02-25 | 2016-02-25 | 一种考虑静态安全约束的含mmc‑hvdc交直流混合电网优化潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105552915A true CN105552915A (zh) | 2016-05-04 |
CN105552915B CN105552915B (zh) | 2018-02-06 |
Family
ID=55831929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610103626.4A Active CN105552915B (zh) | 2016-02-25 | 2016-02-25 | 一种考虑静态安全约束的含mmc‑hvdc交直流混合电网优化潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105552915B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021754A (zh) * | 2016-05-26 | 2016-10-12 | 重庆大学 | 考虑vsc无功越限调整策略的混联电网概率潮流算法 |
CN107171352A (zh) * | 2017-05-17 | 2017-09-15 | 华北电力大学 | 并列混合直流系统中柔性直流稳态运行区间的计算方法 |
CN107221935A (zh) * | 2017-06-12 | 2017-09-29 | 国网江西省电力公司经济技术研究院 | 基于置信度变换的交直流混联输电系统区间潮流方法 |
CN107342594A (zh) * | 2017-07-25 | 2017-11-10 | 中国电力科学研究院 | 一种柔性直流输送功率优化方法及装置 |
CN108808715A (zh) * | 2018-06-07 | 2018-11-13 | 国电南瑞科技股份有限公司 | 考虑直流网络故障功率的多端柔直系统静态安全分析方法 |
CN109066692A (zh) * | 2018-09-14 | 2018-12-21 | 吉林大学 | 一种分布式能源接入的电力网络潮流优化方法 |
CN108054789B (zh) * | 2017-12-22 | 2019-08-16 | 清华大学 | 一种内嵌无功和电压的安全约束经济调度方法 |
CN110690724A (zh) * | 2019-10-29 | 2020-01-14 | 国网新疆电力有限公司 | 考虑mmc内部动态约束的换流站安全稳定控制方法 |
CN108695905B (zh) * | 2018-04-17 | 2021-04-06 | 华南理工大学 | 一种含110kV等级电网安全约束优化调度方法 |
CN114268120A (zh) * | 2021-12-03 | 2022-04-01 | 国网福建省电力有限公司经济技术研究院 | 一种mmc交流侧近端不对称故障短路电流计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102522746A (zh) * | 2011-12-13 | 2012-06-27 | 河海大学 | 基于原对偶内点法的vsc-hvdc交直流最优潮流方法 |
CN102545207A (zh) * | 2011-12-22 | 2012-07-04 | 河海大学 | 基于预测-校正内点法的vsc-hvdc交直流最优潮流方法 |
CN104269861A (zh) * | 2014-10-13 | 2015-01-07 | 国家电网公司 | 基于柔性环网控制器电磁环网无功功率环流优化控制方法 |
-
2016
- 2016-02-25 CN CN201610103626.4A patent/CN105552915B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102522746A (zh) * | 2011-12-13 | 2012-06-27 | 河海大学 | 基于原对偶内点法的vsc-hvdc交直流最优潮流方法 |
CN102545207A (zh) * | 2011-12-22 | 2012-07-04 | 河海大学 | 基于预测-校正内点法的vsc-hvdc交直流最优潮流方法 |
CN104269861A (zh) * | 2014-10-13 | 2015-01-07 | 国家电网公司 | 基于柔性环网控制器电磁环网无功功率环流优化控制方法 |
Non-Patent Citations (1)
Title |
---|
卫志农等: "《含VSC-HVDC的交直流系统内点法最优潮流计算》", 《中国电机工程学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021754A (zh) * | 2016-05-26 | 2016-10-12 | 重庆大学 | 考虑vsc无功越限调整策略的混联电网概率潮流算法 |
CN106021754B (zh) * | 2016-05-26 | 2019-04-16 | 重庆大学 | 考虑vsc无功越限调整策略的混联电网概率潮流算法 |
CN107171352B (zh) * | 2017-05-17 | 2020-07-10 | 华北电力大学 | 并列混合直流系统中柔性直流稳态运行区间的计算方法 |
CN107171352A (zh) * | 2017-05-17 | 2017-09-15 | 华北电力大学 | 并列混合直流系统中柔性直流稳态运行区间的计算方法 |
CN107221935A (zh) * | 2017-06-12 | 2017-09-29 | 国网江西省电力公司经济技术研究院 | 基于置信度变换的交直流混联输电系统区间潮流方法 |
CN107221935B (zh) * | 2017-06-12 | 2019-12-31 | 国网江西省电力公司经济技术研究院 | 基于置信度变换的交直流混联输电系统区间潮流方法 |
CN107342594A (zh) * | 2017-07-25 | 2017-11-10 | 中国电力科学研究院 | 一种柔性直流输送功率优化方法及装置 |
CN107342594B (zh) * | 2017-07-25 | 2022-02-01 | 中国电力科学研究院 | 一种柔性直流输送功率优化方法及装置 |
CN108054789B (zh) * | 2017-12-22 | 2019-08-16 | 清华大学 | 一种内嵌无功和电压的安全约束经济调度方法 |
CN108695905B (zh) * | 2018-04-17 | 2021-04-06 | 华南理工大学 | 一种含110kV等级电网安全约束优化调度方法 |
CN108808715A (zh) * | 2018-06-07 | 2018-11-13 | 国电南瑞科技股份有限公司 | 考虑直流网络故障功率的多端柔直系统静态安全分析方法 |
CN108808715B (zh) * | 2018-06-07 | 2021-07-27 | 国电南瑞科技股份有限公司 | 考虑直流网络故障功率的多端柔直系统静态安全分析方法 |
CN109066692A (zh) * | 2018-09-14 | 2018-12-21 | 吉林大学 | 一种分布式能源接入的电力网络潮流优化方法 |
CN110690724A (zh) * | 2019-10-29 | 2020-01-14 | 国网新疆电力有限公司 | 考虑mmc内部动态约束的换流站安全稳定控制方法 |
CN110690724B (zh) * | 2019-10-29 | 2021-11-12 | 国网新疆电力有限公司 | 考虑mmc内部动态约束的换流站安全稳定控制方法 |
CN114268120A (zh) * | 2021-12-03 | 2022-04-01 | 国网福建省电力有限公司经济技术研究院 | 一种mmc交流侧近端不对称故障短路电流计算方法 |
CN114268120B (zh) * | 2021-12-03 | 2024-02-27 | 国网福建省电力有限公司经济技术研究院 | 一种mmc交流侧近端不对称故障短路电流计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105552915B (zh) | 2018-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105552915A (zh) | 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 | |
US20060282239A1 (en) | Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system | |
CN104868500B (zh) | 一种适用于微电网逆变器并联运行控制方法 | |
Fuerte-Esquivel et al. | Unified power flow controller: a critical comparison of Newton–Raphson UPFC algorithms in power flow studies | |
CN107171341B (zh) | 一种基于分布式计算的输配电网一体化无功优化方法 | |
CN105140907B (zh) | 直流微网多智能体自适应下垂一致性协调控制方法及装置 | |
CN110504691A (zh) | 一种计及vsc控制方式的交直流配电网最优潮流计算方法 | |
CN108134401A (zh) | 交直流混合系统多目标潮流优化及控制方法 | |
CN105119275B (zh) | 一种计及统一潮流控制器的电力系统动态最优潮流的算法 | |
CN105322546B (zh) | 交直流解耦混合潮流算法 | |
CN102570477B (zh) | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 | |
CN110718918B (zh) | 一种基于hvac和lcc-vsc hvdc混合系统的潮流计算方法 | |
CN102545207B (zh) | 基于预测-校正内点法的vsc-hvdc交直流最优潮流方法 | |
CN103793755B (zh) | 基于原对偶内点法的混合直流输电系统最优潮流方法 | |
CN108493940B (zh) | 一种含柔性直流的交直流混合系统潮流交替迭代计算方法 | |
CN105514971A (zh) | 一种适用于各种运行模式微电网的潮流计算方法 | |
CN109802381A (zh) | 一种基于模糊控制的直流微电网多源动态协调控制方法 | |
CN106960394A (zh) | 一种基于蒙特卡罗的交直流混联电网输电能力评估方法 | |
CN107171328A (zh) | 一种基于adpss的分布式潮流控制器建模及仿真方法 | |
CN109617092A (zh) | 一种交直流混合电网动态无功优化方法及系统 | |
CN107093901A (zh) | 一种分布式潮流控制器的机电暂态模型与仿真方法 | |
CN106208099B (zh) | 一种基于二层规划的电力系统无功优化方法及其应用 | |
CN106786599A (zh) | 交直流混合微电网双向dc‑ac互联装置智能控制方法 | |
CN107482665A (zh) | 一种含柔性直流的交直流混合电网越限校正控制方法 | |
CN105958496A (zh) | 基于连续潮流的含vsc系统最大输电能力计算方法 |
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 |