CN108054789B - 一种内嵌无功和电压的安全约束经济调度方法 - Google Patents
一种内嵌无功和电压的安全约束经济调度方法 Download PDFInfo
- Publication number
- CN108054789B CN108054789B CN201711401494.4A CN201711401494A CN108054789B CN 108054789 B CN108054789 B CN 108054789B CN 201711401494 A CN201711401494 A CN 201711401494A CN 108054789 B CN108054789 B CN 108054789B
- Authority
- CN
- China
- Prior art keywords
- node
- indicate
- constraint
- power
- iteration
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000005457 optimization Methods 0.000 claims description 10
- 238000005553 drilling Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/26—Pc applications
- G05B2219/2639—Energy management, use maximum of cheap power, keep peak load low
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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/30—Reactive power compensation
-
- 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
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- 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
-
- 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
- Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
- Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
本发明提出一种内嵌无功和电压的安全约束经济调度方法,属于电力***安全约束经济调度领域,使得基于电网中各节点及线路潮流在安全约束范围内,达到电网运行费用最小。该方法首先构建由目标函数和约束条件构成的安全约束经济调度的最优潮流模型;对模型迭代求解后,对所有支路将检查标幺化误差:若存在线路不满足标幺化误差,则更新潮流方程约束,进行下一次迭代,直至所有线路均满足标幺化误差,模型收敛,迭代结束,得到包含节点电压幅值、电压相角、机组有功出力、无功出力的一整套电力***安全约束经济调度方案。本发明能最大程度满足电力***经济运行要求,降低电力***运行成本;精度高,收敛性好,具有很强的实际应用价值。
Description
技术领域
本发明属于电力***安全约束经济调度领域,特别涉及一种内嵌无功和电压的安全约束经济调度方法。
背景技术
为实现电网安全稳定运行,电网中的节点必须遵循一些规范制度,安全约束经济调度就是基于电网中各节点及线路潮流在安全约束范围内,使电网运行费用最小的调控行为。安全约束经济调度需要求解交流最优潮流问题,而涉及的交流潮流方程的非线性特征导致其为非凸优化问题,求解难度很大。
目前,仅考虑有功的直流潮流网络模型在电力工业界广泛应用。直流潮流网络模型忽略了无功和网损,并假设各节点电压幅值均为标准值(1.0p.u.)。直流潮流方程中有功和电压相角为线性关系。工业界在直流最优潮流解的基础上执行潮流计算以得到交流可行的调度方案。然而,由于直流潮流方程采用了大量近似,直流最优潮流模型可能得到不经济、甚至不安全的调度方案。在安全裕度较低、有功无功紧密耦合的电力***中,该现象愈加明显。目前,电力***运行逐渐逼近其稳定极限,有功调度频繁被***无功和电压约束制约。因此,电力工业界更迫切地需要能够考虑无功和电压的线性网络模型。
目前现有的考虑无功和电压的交流最优潮流模型为非线性模型,求解的方法有凸化松弛方法和内点法等。凸化松弛包括半正定规划(SemidifinateProgramming,SDP)松弛,二阶锥规划(SecondOrderConeProgramming,SOCP)松弛和二次规划(QuadraticProgramming,QP)松弛,但这类方法带来的误差较大,难以在实际电网中投入应用;而内点法的收敛性无法得到保证。
综上所述,现有的安全约束经济调度方法难以在考虑无功和电压的同时保证高计算精度、高计算效率和良好收敛性。
发明内容
本发明的目的在于为克服已有技术的不足之处,提出一种内嵌无功和电压的安全约束经济调度方法。本发明提出了内嵌无功和电压的最优潮流模型,对无功和有功实现了整体性优化,基于该模型得到的调度方案能最大程度满足电力***经济运行要求,降低电力***运行成本;精度高,收敛性好,具有很强的实际应用价值。
本发明提出的一种内嵌无功和电压的安全约束经济调度方法,其特征在于,该方法包括以下步骤:
1)构建安全约束经济调度的最优潮流模型,该模型由目标函数和约束条件构成;具体步骤如下:
1-1)确定模型的目标函数,表达式如下:
其中,Pg为机组g的有功出力,Qg为机组g的无功出力,为机组集合;和分别为Pg和Qg的成本函数;该目标函数表示机组有功出力和无功出力的总成本最小化;
1-2)确定模型的约束条件;具体如下:
1-2-1)基于降阶网络的潮流方程约束;
假设模型第k次迭代时的初值为(vi,k-1,θi,k-1),k=1,2,3…,i∈N,表示所有节点的集合;其中,当k>1时,vi,k-1为第k-1次迭代优化得到的节点i的电压幅值,θi,k-1为第k-1次迭代优化得到的节点i的相角;k=1时,初值(vi,0,θi,0)通过两种启动方式得到:冷启动和热启动,启动方式任意选择其中一种;
冷启动时,
热启动时,
vi,0=1p.u.,θi,0取直流最优潮流的解,
则潮流方程约束如下:
其中,Pij和Qij分别为线路(i,j)的有功潮流和无功潮流;vi表示节点i的电压幅值,θi表示节点i的相角;vij为线路(i,j)首末节点的电压幅值差,满足vij=vi-vj;θij为线路(i,j)首末节点的相角差,满足θij=θi-θj;和为第k次迭代有功平衡系数,计算表达式分别如(12)和(13)所示;和为第k次迭代无功平衡系数,计算表达式分别如(14)和(15)所示;
表示电压幅值对网损影响,满足如下约束:
k>1时,
k=1时,
其中,vij,est表示线路(i,j)首末节点的电压幅值差的估计值;vi,est表示节点i的电压幅值的估计值;vij,k-1为第k-1次迭代优化得到的线路(i,j)的电压幅值差;
为电压相角差对网损影响,满足如下约束:
k=1且在冷启动方式下,
其中,θij,est线路(i,j)首末节点的相角差的估计值,gij表示线路(i,j)的电导;
k>1且为任一种启动方式,或k=1且在热启动方式下,
有功平衡系数和的计算表达式分别如下:
无功平衡系数和的计算表达式分别如下:
其中,bij表示线路(i,j)的电纳;和为正弦函数泰勒展开表达式系数,计算表达式如(16)所示;和为余弦函数泰勒展开表达式系数,计算表达式如(17)所示:
1-2-2)节点平衡方程约束,表达式如下:
其中,表示与节点i通过线路相连的所有节点的集合;gii表示节点i的接地电导,bii表示节点i的接地电纳;
1-2-3)节点注入约束,表达式如下:
其中,表示与节点i相连的机组集合;Pd,i表示节点i的有功负荷,Qd,i表示节点i的无功负荷;
1-2-4)线路潮流复功率线性化约束,表达式如下:
其中,表示所有线路集合;
计算线路潮流复功率线性化约束采用多条线段近似圆弧,具体步骤如下:
1-2-4-1)选取参数为弧度值,取值范围[0,π],对和对应的圆弧进行分段线性近似;
1-2-4-2)选取参数M,N;M,N分别为上下半平面的分段线性数;将定义的圆弧均匀分成N份,用线段连接相邻的端点得到N条线段,即将定义的圆弧均匀分成M份,用线段连接相邻的端点得到M条线段,,即
1-2-4-3)用步骤1-2-4-2)得到的(M+N)条线段对应的线性约束组合近似二次复功率约束,得到复功率线性化约束,如式(22)所示;
1-2-5)机组有功出力、无功出力及电压幅值的上下限约束,表达式如下:
其中,Pg,max表示机组g有功出力上限,Pg,min表示机组g有功出力下限;Qg,max表示机组g无功出力上限,Qg,min表示机组g无功出力下限;vi,max表示节点i电压幅值上限,vi,mim表示节点i电压幅值下限;
2)对步骤1)建立的模型求解,得到电力***安全约束经济调度方案;
对步骤1)建立的模型求解,令当前迭代次数记为k,迭代完成后,得到第k次迭代的最优解vi,θi,Pg,Qg,并根据潮流方程约束式(4)和(5)分别得到模型第k次迭代后线路(i,j)的有功潮流和无功潮流
将第k次迭代的最优解vi,θi,代入线路交流潮流方程:
得到交流有功潮流Pij和交流无功潮流Qij;
根据PijQij,的值,对所有支路将检查标幺化误差,表达式如下:
其中,Δij表示线路(i,j)的潮流方程近似误差;
若所有线路均满足Δij<Δtol,Δtol表示线性近似误差上限,则模型收敛,迭代结束,输出vi,θi,Pg,Qg,得到了包含节点电压幅值、电压相角、机组有功出力、无功出力的电力***安全约束经济调度方案;否则,令k=k+1,将第k次迭代的最优解vi,θi作为初值vi,k-1和θi,k-1,更新潮流方程约束,重新返回步骤2),进行下一次迭代。
本发明的特点及有益效果在于:
1本发明方法不但考虑了有功和相角,而且内嵌无功和电压,得到的调度方案是同时保证了经济性和安全性,在安全裕度较低、有功无功紧密耦合的电力***中同样可以使用,有很高的应用价值。
2本发明提出的线性化方法大幅降低了最优潮流问题的非凸性和非线性度,从而确保其求解算法满足电力工业的应用要求;在迭代初值不佳的情况下,由于模型保留较高精度,仅需2~3次迭代即可大幅减小潮流方程线性化误差,达到收敛条件;每次迭代需要求解仅为线性模型,可方便利用现有的商业软件高效求解。
附图说明
图1是本发明一种内嵌无功和电压的安全约束经济调度方法的整体流程图。
图2是本发明中二次复功率约束的线性化方法示意图。
图3是本发明实施例的交流可行性恢复流程图。
具体实施方式
本发明提出的一种内嵌无功和电压的安全约束经济调度方法,下面结合附图和具体实施例进一步详细说明如下。
本发明提出的一种内嵌无功和电压的安全约束经济调度方法,其核心是采用逐次线性方法求解最优潮流模型,初次迭代可通过冷启动或热启动方式构建潮流方程约束,得到对应的最优潮流模型;之后的每次迭代中,基于上次迭代最优解更新潮流方程约束,使得模型精度进一步提高。本方法将线性近似潮流方程的误差作为收敛判据:如果最优解满足收敛判据,则认为线性近似的精度满足要求。此时,通过恢复交流可行性模块对优化结果做少量修正,得到的优化结果交流可行且接近交流最优。为修正由此引起的约束的潜在越限,方法将对是否存在越限进行检验;若存在,则收紧电网运行中对应的约束限值,重新建立最优潮流模型,并重新求解。最终调度方法得到的操作参数应当是安全且经济的。本方法整体流程如图1所示,包括以下步骤:
1)构建安全约束经济调度的最优潮流模型,该模型由目标函数和约束条件构成;具体步骤如下:
1-1)确定模型的目标函数,表达式如下:
其中,Pg为机组g的有功出力,Qg为机组g的无功出力,为机组集合;和分别为Pg和Qg的成本函数;该目标函数要求机组有功出力和无功出力的总成本最小化。
1-2)确定模型的约束条件;具体如下:
1-2-1)基于降阶网络的潮流方程约束;
假设模型第k次(k=1,2,3…)迭代时的初值为(vi,k-1,θi,k-1),i∈N(N表示所有节点的集合);其中,当k>1时,vi,k-1为第k-1次迭代后节点i的电压幅值,θi,k-1为第k-1次迭代后节点i的相角,(vi,k-1,θi,k-1)取第k-1次迭代后迭代得到的最优解;k=1时,初值(vi,0,θi,0)可通过两种启动方式得到:冷启动和热启动,启动方式可任意选择其中一种。
冷启动时,
热启动时,
vi,0=1p.u.,θi,0取直流最优潮流的解,
则潮流方程约束如下:
其中,Pij和Qij分别为线路(i,j)的有功潮流和无功潮流;vi表示节点i的电压幅值,θi表示节点i的相角,均为待求量;vij为线路(i,j)首末节点的电压幅值差,满足vij=vi-vj;θij为线路(i,j)首末节点的相角差,满足θij=θi-θj;和为第k次迭代有功平衡系数,计算表达式如(12)和(13)所示;和为第k次迭代无功平衡系数,计算表达式如(14)和(15)所示。
表示电压幅值对网损影响,需满足如下约束:
k>1时,
k=1时,
其中,vij,est表示线路(i,j)首末节点的电压幅值差的估计值,通常通过历史数据得到。vi,est表示节点i的电压幅值差的估计值。vij,k-1为第k-1次迭代优化得到的线路(i,j)的电压幅值差。
为电压相角差对网损影响,需满足如下约束:
k=1且在冷启动方式下
其中,θij,est线路(i,j)首末节点的相角差的估计值,通常通过历史数据得到;gij表示线路(i,j)的电导;
k>1且为任一种启动方式,或k=1且在热启动方式下,
有功平衡系数和的计算表达式分别如下:
无功平衡系数和的计算表达式分别如下:
其中,gij表示线路(i,j)的电导,bij表示线路(i,j)的电纳。和为正弦函数泰勒展开表达式系数,计算表达式如(16)所示;和为余弦函数泰勒展开表达式系数,计算表达式如(17)所示:
1-2-2)节点平衡方程约束,表达式如下:
其中,表示与节点i通过线路相连的所有节点的集合;gii表示节点i的接地电导,bii表示节点i的接地电纳。将v2视为独立变量,则式(18)和(19)均为线性方程。
1-2-3)节点注入约束,表达式如下:
其中,表示与节点i相连的机组集合;Pd,i表示节点i的有功负荷,Qd,i表示节点i的无功负荷。
1-2-4)线路潮流复功率线性化约束,表达式如下:
其中,表示所有线路集合。
计算线路潮流复功率线性化约束的核心思想即用多条线段近似圆弧,如图2所示,具体步骤如下:
1-2-4-1)选取参数为弧度值,取值范围[0,π]。基于线路潮流有功和无功的比例,通常只对和对应的圆弧做分段线性近似;
1-2-4-2)选取参数M,N;M,N分别为上下半平面的分段线性数,取值范围[10,100]。将定义的圆弧均匀分成N份。用线段连接相邻的端点(如图2中线段所示),由此可得到N条线段,即将定义的圆弧均匀分成M份,用线段连接相邻的端点得到M条线段,,即
1-2-4-3)用步骤1-2-4-2)得到的(M+N)条线段对应的线性约束组合近似二次复功率约束,得到复功率线性化约束,如式(22)所示。
1-2-5)机组有功出力、无功出力及电压幅值的上下限约束,表达式如下:
其中,Pg,max表示机组g有功出力上限,Pg,min表示机组g有功出力下限;Qg,max表示机组g无功出力上限,Qg,min表示机组g无功出力下限;vi,max表示节点i电压幅值上限,vi,mim表示节点i电压幅值下限。将θ和v2视为独立变量,最优潮流模型的约束为线性。
2)对步骤1)建立的模型求解,得到电力***安全约束经济调度方案;
利用现有的商业优化求解器(如Cplex,Gurobi)求解步骤1)建立的线性模型。令当前迭代次数记为k,迭代求解最优潮流模型后,在得到第k次迭代的最优解 的同时,根据潮流方程约束式(4)和(5)得到模型第k次迭代后线路(i,j)的有功潮流和无功潮流这是基于降阶网路近似得到的。
将第k次迭代的模型优解vi,θi,代入线路交流潮流方程:
得到交流有功潮流Pij和交流无功潮流Qij。
根据PijQij,的值,对所有支路将检查标幺化误差,表达式如下:
其中,Δij表示线路(i,j)的潮流方程近似误差;
若所有线路均满足Δij<Δtol,则模型收敛,若所有线路均满足Δij<Δtol,则模型收敛,迭代结束,输出vi,θi,Pg,Qg,得到了包含节点电压幅值、电压相角、机组有功出力、无功出力的电力***安全约束经济调度方案;否则,令k=k+1,将第k次迭代的最优解vi,θi作为初值vi,k-1和θi,k-1,更新潮流方程约束,重新返回步骤2),进行下一次迭代;其中,Δtol表示线性近似误差上限,取值范围[1e-4,1e-2]。
为进一步提高结果的准确度,本发明在步骤2)迭代结束后,还可以包括如下步骤:
3)对步骤2)得到的结果进行交流可行性恢复;
本发明在模型达到收敛后,采用常规的潮流计算恢复解的交流可行性。流程如图3所示。
根据步骤2)得到的最优潮流解(模型最优解),首先固定特定变量取值以准备潮流计算的输入数据。需要固定的特定变量有:参考节点的v和θ;PQ节点的P和Q;PV节点的P和v。基于输入条件,执行交流潮流计算。在交流潮流计算得到的解中,如果存在机组无功越限,则固定该机组的无功出力,将节点类型转化为PQ节点,并重新执行交流潮流计算;最终,得到交流可行的最优潮流解。
4)检查步骤3)得到的结果是否存在越限;
首先设置越限系数上限γmax(取值为小于1的正数),约束越限率超出该值才判定为越限。对步骤3)得到的交流可行的最优潮流解分别检查:
所有线路潮流复功率越限;
根据式(26)和(27)计算得到校正后的线路有功、无功潮流,分别平方再求和后再开根,得到校正之后准确地线路潮流复功率Sij,若其超出额定值Smax,ij的比率超过γmax,则判定其为越限,即如下情况:
Sij=Smax,ij(1+γij),γij≥γmax (29)
其中,γij为线路(i,j)潮流复功率的越限率。
松弛节点上机组的有功出力越限(根据式(23)和γmax判定),以及PQ的电压幅值越限(根据式(25)和γmax判定)类似可得。
若存在上述任一种越限,则进入步骤5);若无越限,则进入步骤6);
对步骤3)得到的交流可行的最优潮流解分别检查如下约束:所有线路潮流复功率约束,松弛节点上机组的有功出力约束,以及PQ的电压约束;
若上述任一种约束越限,即该约束超过设置越限系数上限γmax(取值为小于1的正数)时,则进入步骤5);若以上约束均无越限,则进入步骤6)。
5)越限约束校正;
对步骤4)中判定为线路潮流复功率越限的线路(i,j),
在下一次迭代中,线路(i,j)的线路潮流复功率上限将被修改为:
Sij,max=Sij,max(1-L·γij) (30)
对松弛节点上机组的有功出力越限,类似地在下一次迭代中调整机组有功界限;对PQ的电压幅值越限的情况,也类似地调整电压幅值界限。下一次迭代开始时,利用本次迭代模型优解vi,θi,作为降阶网络的潮流方程约束中的初值(vi,k-1,θi,k-1),更新潮流方程约束,重新返回步骤2),进行下一次迭代;
6)输出电力***安全约束经济调度方案;
迭代结束,步骤2)输出的vi,θi,Pg,Qg,为模型的最优解,得到了包含节点电压幅值、电压相角、机组有功出力、无功出力在内的一整套电力***安全约束经济调度方案。
附录)算例
此处将在若干IEEE标准测试***和规模较大的波兰测试***中验证本发明调度方法的有效性。数据来自开源仿真工具Matpower4.1。
算例采用Gurobi求解。计算机处理器为Intel(R)Core(TM)[email protected]。算例中采用的参数如下表所示:
表1算例参数设置表
a)使用冷启动的计算结果;
测试调度方法在冷启动(即以v=1.0p.u.,θ=0为初值)下的计算效果。大部分测试***均存在阻塞。方法计算结果如表2所:
表2冷启动下的计算结果列表
b)使用热启动的计算结果;
在基于直流最优潮流的热启动(v=1.0p.u.,θ取自直流最优潮流)下,本发明方法的计算结果如表3所示:
表3热启动下的计算结果列表
上述实验结果说明本专利给出的内嵌无功和电压的安全约束经济调度方法计算效果好,实际应用性强。
Claims (1)
1.一种内嵌无功和电压的安全约束经济调度方法,其特征在于,该方法包括以下步骤:
1)构建安全约束经济调度的最优潮流模型,该模型由目标函数和约束条件构成;具体步骤如下:
1-1)确定模型的目标函数,表达式如下:
其中,Pg为机组g的有功出力,Qg为机组g的无功出力,为机组集合;和分别为Pg和Qg的成本函数;该目标函数表示机组有功出力和无功出力的总成本最小化;
1-2)确定模型的约束条件;具体如下:
1-2-1)基于降阶网络的潮流方程约束;
假设模型第k次迭代时的初值为(vi,k-1,θi,k-1),k=1,2,3...,i∈N,表示所有节点的集合;其中,当k>1时,vi,k-1为第k-1次迭代优化得到的节点i的电压幅值,θi,k-1为第k-1次迭代优化得到的节点i的相角;k=1时,初值(vi,0,θi,0)通过两种启动方式得到:冷启动和热启动,启动方式任意选择其中一种;
冷启动时,
热启动时,
vi,0=1p.u.,θi,0取直流最优潮流的解,
则潮流方程约束如下:
其中,Pij和Qij分别为线路(i,j)的有功潮流和无功潮流;vi表示节点i的电压幅值,θi表示节点i的相角;vij为线路(i,j)首末节点的电压幅值差,满足vij=vi-vj;θij为线路(i,j)首末节点的相角差,满足θij=θi-θj;和为第k次迭代有功平衡系数,计算表达式分别如(12)和(13)所示;和为第k次迭代无功平衡系数,计算表达式分别如(14)和(15)所示;
表示电压幅值对网损影响,满足如下约束:
k>1时,
k=1时,
其中,vij,est表示线路(i,j)首末节点的电压幅值差的估计值;vi,est表示节点i的电压幅值的估计值;vij,k-1为第k-1次迭代优化得到的线路(i,j)的电压幅值差;
为电压相角差对网损影响,满足如下约束:
k=1且在冷启动方式下,
其中,θij,est线路(i,j)首末节点的相角差的估计值,gij表示线路(i,j)的电导;
k>1且为任一种启动方式,或k=1且在热启动方式下,
有功平衡系数和的计算表达式分别如下:
无功平衡系数和的计算表达式分别如下:
其中,bij表示线路(i,j)的电纳;和为正弦函数泰勒展开表达式系数,计算表达式如(16)所示;和为余弦0数泰勒展开表达式系数,计算表达式如(17)所示:
1-2-2)节点平衡方程约束,表达式如下:
其中,表示与节点i通过线路相连的所有节点的集合;gii表示节点i的接地电导,bii表示节点i的接地电纳;
1-2-3)节点注入约束,表达式如下:
其中,表示与节点i相连的机组集合;Pd,i表示节点i的有功负荷,Qd,i表示节点i的无功负荷;
1-2-4)线路潮流复功率线性化约束,表达式如下:
其中,表示所有线路集合;
计算线路潮流复功率线性化约束采用多条线段近似圆弧,具体步骤如下:
1-2-4-1)选取参数为弧度值,取值范围[0,π],对和对应的圆弧进行分段线性近似;
1-2-4-2)选取参数M,N;M,N分别为上下半平面的分段线性数;将定义的圆弧均匀分成N份,用线段连接相邻的端点得到N条线段,即将定义的圆弧均匀分成M份,用线段连接相邻的端点得到M条线段,,即
1-2-4-3)用步骤1-2-4-2)得到的(M+N)条线段对应的线性约束组合近似二次复功率约束,得到复功率线性化约束,如式(22)所示;
1-2-5)机组有功出力、无功出力及电压幅值的上下限约束,表达式如下:
其中,Pg,max表示机组g有功出力上限,Pg,min表示机组g有功出力下限;Qg,max表示机组g无功出力上限,Qg,min表示机组g无功出力下限;vi,max表示节点i电压幅值上限,vi,mim表示节点i电压幅值下限;
2)对步骤1)建立的模型求解,得到电力***安全约束经济调度方案;
对步骤1)建立的模型求解,令当前迭代次数记为k,迭代完成后,得到第k次迭代的最优解vi,θi,Pg,Qg,并根据潮流方程约束式(4)和(5)分别得到模型第k次迭代后线路(i,j)的有功潮流阳无功潮流
将第k次迭代的最优解vi,θi,代入线路交流潮流方程:
得到交流有功潮流Pij和交流无功潮流Qij;
根据PijQij,的值,对所有支路将检查标幺化误差,表达式如下:
其中,Δij表示线路(i,j)的潮流方程近似误差;
若所有线路均满足Δij<Δtol,Δtol表示线性近似误差上限,则模型收敛,迭代结束,输出vi,θi,Pg,Qg,得到了包含节点电压幅值、电压相角、机组有功出力、无功出力的电力***安全约束经济调度方案;否则,令k=k+1,将第k次迭代的最优解vi,θi作为初值vi,k-1和θi,k-1,更新潮流方程约束,重新返回步骤2),进行下一次迭代。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711401494.4A CN108054789B (zh) | 2017-12-22 | 2017-12-22 | 一种内嵌无功和电压的安全约束经济调度方法 |
PCT/CN2018/094555 WO2019119775A1 (zh) | 2017-12-22 | 2018-07-04 | 一种内嵌无功和电压的安全约束经济调度方法 |
US16/751,915 US11152786B2 (en) | 2017-12-22 | 2020-01-24 | Dispatch method and device for power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711401494.4A CN108054789B (zh) | 2017-12-22 | 2017-12-22 | 一种内嵌无功和电压的安全约束经济调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108054789A CN108054789A (zh) | 2018-05-18 |
CN108054789B true CN108054789B (zh) | 2019-08-16 |
Family
ID=62131474
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711401494.4A Active CN108054789B (zh) | 2017-12-22 | 2017-12-22 | 一种内嵌无功和电压的安全约束经济调度方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11152786B2 (zh) |
CN (1) | CN108054789B (zh) |
WO (1) | WO2019119775A1 (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108054789B (zh) * | 2017-12-22 | 2019-08-16 | 清华大学 | 一种内嵌无功和电压的安全约束经济调度方法 |
CN110991727A (zh) * | 2019-11-28 | 2020-04-10 | 海南电网有限责任公司 | 一种基于潮流网损模型和线路约束模型的电网规划方法 |
CN111130118B (zh) * | 2020-01-09 | 2021-02-02 | 清华大学 | 一种基于分段线性化的电力***最优潮流计算方法 |
CN111313425B (zh) * | 2020-01-15 | 2023-11-14 | 国网重庆市电力公司 | 基于变量空间最优选择的潮流模型线性化误差最小化方法 |
JP7254041B2 (ja) * | 2020-02-21 | 2023-04-07 | 株式会社日立製作所 | 系統計画支援システムおよび系統計画支援方法 |
CN111242389A (zh) * | 2020-03-02 | 2020-06-05 | 南方电网科学研究院有限责任公司 | 一种智能储能软开关规划方法、***、设备及介质 |
CN113765151B (zh) * | 2020-06-04 | 2023-09-22 | 株洲变流技术国家工程研究中心有限公司 | 一种能量优化调度方法及装置 |
CN111952956B (zh) * | 2020-07-02 | 2024-07-09 | 清华大学 | 一种考虑电压敏感负荷备用的电力***调度方法 |
CN111969609B (zh) * | 2020-07-06 | 2021-12-14 | 南方电网科学研究院有限责任公司 | 交直流输电网的二阶锥最优潮流模型构建方法和装置 |
CN111898811B (zh) * | 2020-07-17 | 2024-06-04 | 国网重庆市电力公司 | Dc-opf模型中修正二次目标函数分段线性化误差的方法 |
CN111799802B (zh) * | 2020-07-22 | 2023-12-05 | 国网重庆市电力公司 | 基于线性组合的潮流方程线性化方法 |
CN111756043B (zh) * | 2020-07-30 | 2021-10-19 | 广东电网有限责任公司电力科学研究院 | 一种考虑无功和电压幅值的最优潮流优化方法及*** |
CN112529614B (zh) * | 2020-11-27 | 2023-10-24 | 国网湖北省电力有限公司 | 电力市场出清流程的建立方法及*** |
CN112994014B (zh) * | 2021-02-19 | 2024-04-05 | 山东大学 | 一种考虑风电影响的电网最优运行点控制方法及*** |
CN112838588B (zh) * | 2021-03-19 | 2021-09-28 | 四川大学 | 一种考虑机组启停和线路投切的短路电流抑制方法 |
CN113191549B (zh) * | 2021-04-29 | 2023-02-24 | 国网山东省电力公司德州供电公司 | 考虑网损与简化交流潮流的配电网源网联合规划优化方法 |
CN113421004B (zh) * | 2021-06-30 | 2023-05-26 | 国网山东省电力公司潍坊供电公司 | 输配协同的主动配电网分布式鲁棒扩展规划***及方法 |
CN113627021A (zh) * | 2021-08-11 | 2021-11-09 | 东南大学 | 一种基于序列凸规划的电气互联***最优能流计算方法 |
CN113609699B (zh) * | 2021-08-17 | 2023-07-21 | 上海交通大学 | 辐射状配电网交流潮流模型的计算方法和*** |
DE102022204250A1 (de) | 2022-04-29 | 2023-11-02 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Steuerung eines Stromnetzes |
CN115360762A (zh) * | 2022-09-16 | 2022-11-18 | 东南大学溧阳研究院 | 交直流***有功功率最大化消纳方法 |
CN115603381A (zh) * | 2022-10-31 | 2023-01-13 | 华能湖北新能源有限责任公司(Cn) | 一种新能源发电并网管理方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105552915A (zh) * | 2016-02-25 | 2016-05-04 | 国网福建省电力有限公司 | 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 |
CN105790265A (zh) * | 2016-04-21 | 2016-07-20 | 三峡大学 | 一种考虑交流潮流约束的不确定性机组组合模型及求解方法 |
CN105870912A (zh) * | 2016-03-03 | 2016-08-17 | 广州华工弈高科技有限公司 | 主动配电网孤岛划分混合整数规划模型的建模方法 |
CN106684857A (zh) * | 2016-09-05 | 2017-05-17 | 国网江苏省电力公司电力科学研究院 | 一种含统一潮流控制器的线性化最优潮流模型 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101752903B (zh) * | 2009-11-27 | 2011-06-01 | 清华大学 | 一种时序递进式电力调度方法 |
CN103927590B (zh) * | 2014-03-20 | 2017-01-25 | 清华大学 | 信息不对称情况下的电力***优化控制及费用计算方法 |
CN106682808A (zh) * | 2016-09-20 | 2017-05-17 | 北京恒泰实达科技股份有限公司 | 在线滚动优化调度模型 |
CN107069706B (zh) * | 2017-02-17 | 2019-08-16 | 清华大学 | 一种基于多参数规划的输配电网协调的动态经济调度方法 |
CN108054789B (zh) * | 2017-12-22 | 2019-08-16 | 清华大学 | 一种内嵌无功和电压的安全约束经济调度方法 |
CN108054757B (zh) * | 2017-12-22 | 2019-08-16 | 清华大学 | 一种内嵌无功和电压的n-1闭环安全校核方法 |
-
2017
- 2017-12-22 CN CN201711401494.4A patent/CN108054789B/zh active Active
-
2018
- 2018-07-04 WO PCT/CN2018/094555 patent/WO2019119775A1/zh active Application Filing
-
2020
- 2020-01-24 US US16/751,915 patent/US11152786B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105552915A (zh) * | 2016-02-25 | 2016-05-04 | 国网福建省电力有限公司 | 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 |
CN105870912A (zh) * | 2016-03-03 | 2016-08-17 | 广州华工弈高科技有限公司 | 主动配电网孤岛划分混合整数规划模型的建模方法 |
CN105790265A (zh) * | 2016-04-21 | 2016-07-20 | 三峡大学 | 一种考虑交流潮流约束的不确定性机组组合模型及求解方法 |
CN106684857A (zh) * | 2016-09-05 | 2017-05-17 | 国网江苏省电力公司电力科学研究院 | 一种含统一潮流控制器的线性化最优潮流模型 |
Non-Patent Citations (1)
Title |
---|
Optimal Reactive Power Dispatch With Accurately;Zhifang Yang et al.;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20170531;第32卷(第3期);第2435-2444页 |
Also Published As
Publication number | Publication date |
---|---|
CN108054789A (zh) | 2018-05-18 |
US20200161863A1 (en) | 2020-05-21 |
WO2019119775A1 (zh) | 2019-06-27 |
US11152786B2 (en) | 2021-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108054789B (zh) | 一种内嵌无功和电压的安全约束经济调度方法 | |
He et al. | Decentralized optimization of multi-area electricity-natural gas flows based on cone reformulation | |
Contreras et al. | Improved assessment of the flexibility range of distribution grids using linear optimization | |
US9954362B2 (en) | Systems and methods for optimal power flow on a radial network | |
ES2816027T3 (es) | Procedimiento para controlar la distribución de energía | |
US10158229B2 (en) | Distributed gradient descent for solving optimal power flow in radial networks | |
CN105322546B (zh) | 交直流解耦混合潮流算法 | |
Sun et al. | Adaptive decoupled power control method for inverter connected DG | |
Shchetinin et al. | On the construction of linear approximations of line flow constraints for AC optimal power flow | |
CN111355241B (zh) | 基于潮流可行域的交直流混合微网下垂系数选取方法 | |
CN105470955A (zh) | 一种含分布式电源的输配电网牛拉-前推协调潮流算法 | |
CN105138799A (zh) | 适用于模块化多电平换流器的直流电抗器参数设计方法 | |
CN109038641A (zh) | 高压直流输电***快速提升功率能力的计算方法及*** | |
CN109560547A (zh) | 一种考虑输配协同的主动配电网n-1安全评估方法 | |
CN107994583B (zh) | 直流联络线在协调层建模的交直流***分布式潮流算法 | |
Karagiannopoulos et al. | Data-driven control design schemes in active distribution grids: Capabilities and challenges | |
CN110556853A (zh) | 一种为电磁暂态仿真提供初值的计算方法及*** | |
JP6782181B2 (ja) | 太陽光発電量推定装置および方法 | |
Nian et al. | Transient modeling and analysis of VSC based DC microgrid during short circuit fault | |
CN113555891B (zh) | 基于平方和的含不确定参数交直流***稳定域估计方法 | |
Baghsorkhi et al. | Embedding AC power flow in the complex plane part II: a reliable framework for voltage collapse analysis | |
Colot et al. | Optimal power flow pursuit via feedback-based safe gradient flow | |
Tavares et al. | Practical method for computing the maximum loading point using a load flow with step size optimisation | |
CN109217384A (zh) | 基于有功相角下垂控制的多分布式电源功率分配方法 | |
Sepulveda et al. | Sequential convex optimization for the dynamic optimal power flow of active distribution networks |
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 |