CN105262107A - 一种风电并网的无功优化方法 - Google Patents
一种风电并网的无功优化方法 Download PDFInfo
- Publication number
- CN105262107A CN105262107A CN201510641997.3A CN201510641997A CN105262107A CN 105262107 A CN105262107 A CN 105262107A CN 201510641997 A CN201510641997 A CN 201510641997A CN 105262107 A CN105262107 A CN 105262107A
- Authority
- CN
- China
- Prior art keywords
- target function
- wind
- disaggregation
- idle work
- work optimization
- 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.)
- Pending
Links
Classifications
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
-
- 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/30—Reactive power compensation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种风电并网的无功优化方法,包括:步骤一:选取有功网损作为目标函数,确定目标函数表达式;步骤二:对目标函数表达式进行处理,以目标函数的均值作为电力***的经济性指标,以目标函数的方差作为电力***的安全性指标;步骤三:创建以所述经济性指标和安全性指标为多目标的无功优化数学模型,对所述无功优化数学模型进行求解,获得决策变量的优化解集,继而进行相应的优化操作。本发明通过一个目标函数的均值和方差即可同时获得风电并网的经济性和安全性指标,运算简单,易于控制无功优化的操作。
Description
技术领域
本发明涉及一种风电并网的无功优化方法,应用于风力发电***。
背景技术
随着世界能源消耗的持续增加,能源短缺矛盾以及气候环境问题的不断凸显,深入开发、利用可持续的清洁能源发电成为各个国家的共识和重要能源战略。风能由于具有可再生、无污染、分布广的特点,风力发电逐渐显示出其光明的应用前景,在电力***中发挥着越来越重要的作用。风电场的并网运行给电力***的规划、调度和运行带来了许多不同于其他发电场的影响,主要表现在以下方面:电压、***稳定性改变;电力***的有功、无功潮流变化等。同时,风力发电机多为异步发电机,它的运行需要无功电源的支持,这也要增加电网对无功的需求,使得电力***的有功潮流、无功潮流分布更加复杂。因此,在现有的电力***中,研究风电机组的电力***潮流、风电并网对电力***电压稳定性和潮流分布的影响具有深远的现实意义。
无功优化是指电力***在一定运行方式,满足各种约束条件的前提下,使***的一个或多个性能指标(如***有功网损最小化、节点电压偏移量最小、静态电压稳定裕度最大、***的总无功补偿量最小等)达到最优的无功调节手段。无功优化是电力***潮流优化理论的重要组成部分,它涉及无功补偿装置投入地点选择、投入容量确定、变压器分接头的调节和发电机机端电压的配合等。通过电力***无功优化,可以实现两个目的:第一,通过改变各个节点之间的电压差可以有效降低有功损耗。第二,通过改变电网的无功分布可以改变整个电力***的潮流分布,从而实现电压运行在安全范围之内。
随着电力***的发展,往往需要同时考虑电力***的经济性和安全性,所以出现了同时考虑有功网损和节点电压偏离最小及电压稳定裕度最大等的多目标无功优化模型。现有技术的多目标无功优化的处理方法大多采用将多目标问题转化为单目标问题进行求解,如建立有功网损最小化的目标函数,往往只能从经济性的角度分析,如果同时考虑安全性,需要结合其他的的目标模型,即将常用的目标函数组合起来实现某些目标,这种方法处理过程繁琐,运算较为复杂,无法通过一个目标函数的处理同时指示电力***的经济性和安全性问题。
发明内容
本发明所要解决的技术问题是提供一种通过一个目标函数的处理同时获得风电并网的经济性和安全性指标的无功优化方法。
为了实现上述目的,本发明采用如下技术方案予以实现:一种风电并网的无功优化方法,包括:
步骤一:选取有功网损作为目标函数,确定目标函数表达式;
步骤二:对目标函数表达式进行处理,以目标函数的均值作为电力***的经济性指标,以目标函数的方差作为电力***的安全性指标;
步骤三:创建以所述经济性指标和安全性指标为多目标的无功优化数学模型,对所述无功优化数学模型进行求解,获得决策变量的优化解集,继而进行相应的优化操作。
优选的,步骤一中所述目标函数表达式为:
其中,f为目标函数;n为电力***的总节点数;Ui、Uj分别为节点i、j的电压,单位为kV;θij为节点i、j的电压相位角的差值;Gij为节点i、j之间支路的电导,单位为S。
优选的,步骤二中采用公式(1)、(2)、(3)对目标函数表达式进行处理,获得风速与节点电压的关系:
公式(1)为Pm=0.5ρAV3Cp,
其中,Pm为风机在风力作用下产生的机械功率,单位为MW,ρ为空气密度,单位为kg/m3,A为风机叶片的扫掠面积,单位为m2,V为风速,单位为m/s,Cp为风机的风能利用系数;
公式(2)为
其中,P为风机的有功功率,单位为MW,Q为风机的无功功率,单位为Mvar,为恒定的功率因数;
公式(3)为
其中,Qij为节点i、j流动的无功功率,单位为Mvar,Ui、Uj分别为节点i、j的电压,单位为kV;c为线路电抗,单位为Ω,δ为节点i与j的电压相位角的差值;
在上述公式中,Pm=P,Qij=Q;然后,根据风速与节点电压的关系建立目标函数的均值-方差模型:
其中,Ploss为目标函数均值,Vloss为目标函数方差。
进一步的,步骤二中公式(2)适用于以恒功率因数运行方式运行的双馈风场。
优选的,步骤三中所述无功优化数学模型可描述为:
{y=F(x)=(f1(x),f2(x),...,fk(x))s.t.:g(x)=(g1(x),g2(x),...,gl(x))≤0h(x)=(hl+1(x),hl+2(x),...,hm(x))=0,其中 X={(x1,x2,...,xn)|si≤xi≤ui},s=(s1,s2,...,sn),u=(u1,u2,...,un)},
其中,g(x)≤0和h(x)=0是F(x)的约束条件,x为决策变量,y为总目标函数,X为决策空间,Rn为n维实数空间,u和s分别为x的上界和下界,m为约束条件总数,k为目标函数总数,所述g(x)=(g1(x),g2(x),...,gl(x))≤0定义了l个不等式约束,所述h(x)=(hl+1(x),hl+2(x),...,hm(x))=0定义了m-l个等式约束。
进一步的,决策变量x为无功补偿电源的容量、变压器分接头的位置和发电机端电压的幅值中的一个或多个。
进一步的,步骤三还包括:对所述无功优化数学模型的约束条件进行处理,对候选解集进行划分。
进一步的,所述无功优化数学模型的约束条件处理方法包括:
步骤1:将约束条件转化为一个函数 其中,Gi(x)=max{0,gi(x)},1≤i≤l,Hj(x)=|hj(x)|,l≤j≤m;
步骤2:根据函数fk+1(x)的大小将候选解集划分为可行解集和不可行解集,若fk+1(x)=0,则x满足约束条件为可行解,若fk+1(x)>0,则x不满足约束条件为不可行解;
步骤3:根据帕累托占优方法将可行解集划分为帕累托占优解集和非帕累托占优解集,根据fk+1(x)的值将不可行解集划分为有益解集和非有益解集;
步骤4:保存有益解集和帕累托占优解集作为优化解集。
与现有技术相比,本发明通过对有功网损目标函数表达式进行处理,以目标函数的均值作为电力***的经济性指标,以目标函数的方差作为电力***的安全性指标,既考虑了经济性又考虑了安全性,然后,创建以所述经济性指标和安全性指标为多目标的无功优化数学模型,对所述无功优化数学模型进行求解,获得决策变量的优化解集,继而进行相应的优化操作。本发明通过一个目标函数的均值和方差即可同时获得风电并网的经济性和安全性指标,运算简单,易于控制无功优化的操作。
附图说明
图1是本发明的无功优化方法的步骤示意图;
图2是本发明一实施例中无功优化数学模型约束条件处理方法的步骤示意图。
具体实施方式
下面结合附图对本发明作详细描述:
如图1和图2所示,本发明的风电并网的无功优化方法,包括:
步骤一:选取有功网损作为目标函数,确定目标函数表达式。其中,所述目标函数表达式为:
其中,f为目标函数;n为电力***的总节点数;Ui、Uj分别为节点i、j的电压,单位为kV;θij为节点i、j的电压相位角的差值;Gij为节点i、j之间支路的电导,单位为S。
步骤二:对目标函数表达式进行处理,以目标函数的均值作为电力***的经济性指标,以目标函数的方差作为电力***的安全性指标。
以有功网损的均值Ploss作为经济性指标,可以具有代表性地体现电力***的经济性,以有功网损的方差Vloss来作为电力***的安全性指标,当***的有功网损方差值比较大时,表明***实际有功网损偏离期望值很远,说明该无功控制策略不能很好地处理风速不确定性的波动场景,则该***在风速波动情况下表现出不稳定的特性,所以可以利用有功网损的方差Vloss来作为***安全性的指标。通过一个目标函数的均值和方差即可同时获得电力***的经济性和安全性指标,运算简单,易于控制无功优化的操作。
较佳的,采用公式(1)、(2)、(3)对目标函数表达式进行处理,获得风速与节点电压的关系:
公式(1)为Pm=0.5ρAV3Cp,
其中,Pm为风机在风力作用下产生的机械功率,单位为MW,ρ为空气密度,单位为kg/m3,A为风机叶片的扫掠面积,单位为m2,V为风速,单位为m/s,Cp为风机的风能利用系数;所述Cp为叶尖速比和叶片桨距角的函数,风机的Cp特性曲线由实验数据获得。
公式(2)为
其中,P为风机的有功功率,单位为MW,Q为风机的无功功率,单位为Mvar,为恒定的功率因数;公式(2)适用于以恒功率因数运行方式运行的双馈风场。
公式(3)为
其中,Qij为节点i、j流动的无功功率,单位为Mvar,Ui、Uj分别为节点i、j的电压,单位为kV;c为线路电抗,单位为Ω,δ为节点i与j的电压相位角的差值;
在上述公式中,Pm=P,Qij=Q;然后,根据风速与节点电压的关系建立目标函数的均值-方差模型:
其中,Ploss为目标函数均值,Vloss为目标函数方差。
通过公式(1)、(2)、(3)的处理,获得风速与节点电压的关系,其中,风电场的风速V作为影响因子,决定风电场的节点电压,通过风电场的节点电压最终可以影响目标函数的均值Ploss和方差Vloss。根据目标函数的均值Ploss和方差Vloss的要求,通过无功优化数学模型的求解,可以获得风电场的节点电压的优化措施,继而对风速V进行调节。也可以根据风速V的实际情况,通过无功优化数学模型的求解,获得电力***的无功优化。
步骤三:创建以所述经济性指标和安全性指标为多目标的无功优化数学模型,对所述无功优化数学模型进行求解,获得决策变量的优化解集,继而进行相应的优化操作。
其中,所述无功优化数学模型可描述为:
{y=F(x)=(f1(x),f2(x),...,fk(x))s.t.:g(x)=(g1(x),g2(x),...,gl(x))≤0h(x)=(hl+1(x),hl+2(x),...,hm(x))=0,其中 X={(x1,x2,...,xn)|si≤xi≤ui},s=(s1,s2,...,sn),u=(u1,u2,...,un)},
其中,g(x)≤0和h(x)=0是F(x)的约束条件,x为决策变量,y为总目标函数,X为决策空间,Rn为n维实数空间,u和s分别为x的上界和下界,m为约束条件总数,k为目标函数总数,所述g(x)=(g1(x),g2(x),...,gl(x))≤0定义了l个不等式约束,所述h(x)=(hl+1(x),hl+2(x),...,hm(x))=0定义了m-l个等式约束。
其中,决策变量x为无功补偿电源的容量、变压器分接头的位置和发电机端电压的幅值中的一个或多个。
根据经济性指标和安全性指标的要求,通过对无功优化数学模型的求解,获得决策变量即无功优化控制措施的优化解集,继而进行相应的优化操作,得到合理的无功补偿位置和无功补偿容量以使电力***安全稳定运行。
步骤三还包括:对所述无功优化数学模型的约束条件进行处理,对候选解集进行划分。
较佳的,所述无功优化数学模型的约束条件处理方法包括:
步骤1:将约束条件转化为一个函数 其中,Gi(x)=max{0,gi(x)},1≤i≤l,Hj(x)=|hj(x)|,l≤j≤m;
步骤2:根据函数fk+1(x)的大小将候选解集划分为可行解集和不可行解集,若fk+1(x)=0,则x满足约束条件为可行解,若fk+1(x)>0,则x不满足约束条件为不可行解;
步骤3:根据Pareto占优方法将可行解集划分为帕累托(Pareto)占优解集和非帕累托(Pareto)占优解集,根据fk+1(x)的值将不可行解集划分为有益解集和非有益解集;
步骤4:保存有益解集和帕累托(Pareto)占优解集作为优化解集。
Claims (8)
1.一种风电并网的无功优化方法,其特征在于,包括:
步骤一:选取有功网损作为目标函数,确定目标函数表达式;
步骤二:对目标函数表达式进行处理,以目标函数的均值作为电力***的经济性指标,以目标函数的方差作为电力***的安全性指标;
步骤三:创建以所述经济性指标和安全性指标为多目标的无功优化数学模型,对所述无功优化数学模型进行求解,获得决策变量的优化解集,继而进行相应的优化操作。
2.根据权利要求1所述的风电并网的无功优化方法,其特征在于,步骤一中所述目标函数表达式为:
其中,f为目标函数;n为电力***的总节点数;Ui、Uj分别为节点i、j的电压,单位为KV;θij为节点i、j的电压相位角的差值;Gij为节点i、j之间支路的电导,单位为S。
3.根据权利要求1所述的风电并网的无功优化方法,其特征在于,步骤二中采用公式(1)、(2)、(3)对目标函数表达式进行处理,获得风速与节点电压的关系:
公式(1)为Pm=0.5ρAV3Cp,
其中,Pm为风机在风力作用下产生的机械功率,单位为MW,ρ为空气密度,单位为kg/m3,A为风机叶片的扫掠面积,单位为m2,V为风速,单位为m/s,Cp为风机的风能利用系数;
公式(2)为
其中,P为风机的有功功率,单位为MW,Q为风机的无功功率,单位为为恒定的功率因数;
公式(3)为
其中,Qij为节点i、j流动的无功功率,单位为Mvar,Ui、Uj分别为节点i、j的电压,单位为KV;c为线路电抗,单位为Ω,δ为节点i与j的电压相位角的差值;
在上述公式中,Pm=P,Qij=Q;
然后,根据风速与节点电压的关系建立目标函数的均值-方差模型:
其中,Ploss为目标函数均值,Vloss为目标函数方差。
4.根据权利要求3所述的风电并网的无功优化方法,其特征在于,步骤二中公式(2)适用于以恒功率因数运行方式运行的双馈风场。
5.根据权利要求1所述的风电并网的无功优化方法,其特征在于,步骤三中所述无功优化数学模型可描述为:
{y=F(x)=(f1(x),f2(x),...,fk(x))s.t.:g(x)=(g1(x),g2(x),...,gl(x))≤0h(x)=(hl+1(x),hl+2(x),...,hm(x))=0,其中 X={(x1,x2,...,xn)|si≤xi≤ui},s=(s1,s2,...,sn),u=(u1,u2,...,un)},
其中,g(x)≤0和h(x)=0是F(x)的约束条件,x为决策变量,y为总目标函数,X为决策空间,Rn为n维实数空间,u和s分别为x的上界和下界,m为约束条件总数,K为目标函数总数,所述g(x)=(g1(x),g2(x),...,gl(x))≤0定义了l个不等式约束,所述h(x)=(hl+1(x),hl+2(x),...,hm(x))=0定义了m-l个等式约束。
6.根据权利要求5所述的风电并网的无功优化方法,其特征在于,决策变量x为无功补偿电源的容量、变压器分接头的位置和发电机端电压的幅值中的一个或多个。
7.根据权利要求5所述的风电并网的无功优化方法,其特征在于,步骤三还包括:对所述无功优化数学模型的约束条件进行处理,对候选解集进行划分。
8.根据权利要求7所述的风电并网的无功优化方法,其特征在于,所述无功优化数学模型的约束条件处理方法包括:
步骤1:将约束条件转化为一个函数 其中,Gi(x)=max{0,gi(x)},1≤i≤l,Hj(x)=|hj(x)|,l≤j≤m;
步骤2:根据函数fk+1(x)的大小将候选解集划分为可行解集和不可行解集,若fk+1(x)=0,则x满足约束条件为可行解,若fk+1(x)>0,则x不满足约束条件为不可行解;
步骤3:根据帕累托占优方法将可行解集划分为帕累托占优解集和非帕累托占优解集,根据fk+1(x)的值将不可行解集划分为有益解集和非有益解集;
步骤4:保存有益解集和帕累托占优解集作为优化解集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510641997.3A CN105262107A (zh) | 2015-09-30 | 2015-09-30 | 一种风电并网的无功优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510641997.3A CN105262107A (zh) | 2015-09-30 | 2015-09-30 | 一种风电并网的无功优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105262107A true CN105262107A (zh) | 2016-01-20 |
Family
ID=55101663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510641997.3A Pending CN105262107A (zh) | 2015-09-30 | 2015-09-30 | 一种风电并网的无功优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105262107A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106410819A (zh) * | 2016-11-29 | 2017-02-15 | 重庆大学 | 含电压不可行节点的受端电网动态无功优化方法及装置 |
CN109861297A (zh) * | 2019-04-11 | 2019-06-07 | 上海电机学院 | 一种基于灰狼优化算法的电力***黑启动方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102611119A (zh) * | 2012-03-14 | 2012-07-25 | 华北电力大学 | 电力***多目标无功优化方法 |
WO2015134187A1 (en) * | 2014-03-06 | 2015-09-11 | General Electric Company | System and method for robust wind turbine operation |
-
2015
- 2015-09-30 CN CN201510641997.3A patent/CN105262107A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102611119A (zh) * | 2012-03-14 | 2012-07-25 | 华北电力大学 | 电力***多目标无功优化方法 |
WO2015134187A1 (en) * | 2014-03-06 | 2015-09-11 | General Electric Company | System and method for robust wind turbine operation |
Non-Patent Citations (2)
Title |
---|
张雪霞: "智能优化算法及其在电力***无功优化中的应用研究", 《万方学位论文》 * |
王瑞琪: "分布式发电与微网***多目标优化设计与协调控制研究", 《万方学位论文》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106410819A (zh) * | 2016-11-29 | 2017-02-15 | 重庆大学 | 含电压不可行节点的受端电网动态无功优化方法及装置 |
CN106410819B (zh) * | 2016-11-29 | 2018-10-09 | 重庆大学 | 含电压不可行节点的受端电网动态无功优化方法及装置 |
CN109861297A (zh) * | 2019-04-11 | 2019-06-07 | 上海电机学院 | 一种基于灰狼优化算法的电力***黑启动方法 |
CN109861297B (zh) * | 2019-04-11 | 2022-11-01 | 上海电机学院 | 一种基于灰狼优化算法的电力***黑启动方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Armghan et al. | Nonlinear integral backstepping based control of a DC microgrid with renewable generation and energy storage systems | |
Venkatesan et al. | Optimal power flow control of hybrid renewable energy system with energy storage: A WOANN strategy | |
CN110247438B (zh) | 基于天牛须算法的主动配电网资源优化配置 | |
CN104835080B (zh) | 微电网日内调度计划混合整数规划模型的建模方法 | |
CN103903073A (zh) | 一种含分布式电源及储能的微电网优化规划方法及*** | |
CN107666155B (zh) | 基于Markov模型的多能互补***随机稳定性分析方法 | |
Olulope et al. | Modeling and simulation of hybrid distributed generation and its impact on transient stability of power system | |
CN105243516A (zh) | 基于主动配电网的分布式光伏发电最大消纳能力计算*** | |
CN102856899B (zh) | 一种微电网降低网损的方法 | |
CN105322535A (zh) | 含统一潮流控制器的电力***两阶段最优潮流计算方法 | |
CN103368186A (zh) | 一种风电***无功优化方法 | |
CN110206686A (zh) | 一种用于风力发电机组的自适应最大功率跟踪控制方法 | |
CN105375513A (zh) | 一种基于实时在线等值的110千伏风电场自动电压控制方法 | |
CN108306334A (zh) | 基于粒子群优化算法的风电场内部无功优化策略 | |
Singh et al. | Frequency regulation of micro-grid connected hybrid power system with SMES | |
CN107732945A (zh) | 一种基于模拟退火粒子群算法的储能单元优化方法 | |
CN112132363A (zh) | 一种增强***运行鲁棒性的储能选址定容方法 | |
CN104578091A (zh) | 一种含多源电网的无延时最优无功协调控制***及方法 | |
CN104638654A (zh) | 兼顾风电场和网络节点电压调节的statcom控制方法 | |
CN106780119A (zh) | 基于多主动管理的配电网可再生能源消纳方法 | |
CN104362639B (zh) | 基于改进微分进化方法的电网全网无功优化方法 | |
CN105262107A (zh) | 一种风电并网的无功优化方法 | |
Altimania | Modeling of doubly-fed induction generators connected to distribution system based on eMEGASim® real-time digital simulator | |
Paleba et al. | Optimal placement and sizing distributed wind generation using particle swarm optimization in distribution system | |
CN104037806B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160120 |