CN107171341B - 一种基于分布式计算的输配电网一体化无功优化方法 - Google Patents

一种基于分布式计算的输配电网一体化无功优化方法 Download PDF

Info

Publication number
CN107171341B
CN107171341B CN201710452090.1A CN201710452090A CN107171341B CN 107171341 B CN107171341 B CN 107171341B CN 201710452090 A CN201710452090 A CN 201710452090A CN 107171341 B CN107171341 B CN 107171341B
Authority
CN
China
Prior art keywords
distribution network
transmission
boundary
reactive power
network
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
CN201710452090.1A
Other languages
English (en)
Other versions
CN107171341A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201710452090.1A priority Critical patent/CN107171341B/zh
Publication of CN107171341A publication Critical patent/CN107171341A/zh
Application granted granted Critical
Publication of CN107171341B publication Critical patent/CN107171341B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
    • 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于分布式计算的输配电网一体化无功优化方法，采用广义主从法分解输配电网优化问题为输电网无功优化、配电网无功优化和边界信息交互三部分，采用二次罚函数‑原对偶内点法分别交替求解输网、配网无功优化子问题，通过交换边界影响因子(边界连接节点电压、等值功率和目标函数边界拉格朗日乘子信息)参与协调，实现全网无功优化的分布式计算。其中，无功优化普遍存在的离散变量使边界乘子的计算出现困难，引入二次罚函数规整离散变量，使得在输配网全局优化目标函数可微可分的基础上，实现离散变量的处理。本发明采用分布式计算方法求解输配电网一体化无功优化问题，保持了输、配电网现有的计算模式，且具有较高的计算精度。

Description

一种基于分布式计算的输配电网一体化无功优化方法
技术领域
本发明涉及一种输配电网一体化无功优化方法，特别是涉及一种基于分布式计算的输配电网一体化无功优化方法，属于电力系统优化运行技术领域。
背景技术
随着DG在配电网中渗透比例的持续增长，配电网将兼具电源接纳、电能分配、协同输电网运行的作用。配电网的运行特性和运营模式将变得日益复杂，与主网联系更加紧密，无功调节手段更加丰富。输电网无功优化不能再忽视配网的无功电压支撑作用，传统相对独立的输配电网无功优化经济性和安全性难以保证，输配电网一体化无功优化可在更大范围内实现无功资源的合理配置。基于输、配网分布自治、集中协调的管理模式，实现输配网一体化分布式无功电压控制具有重要的研究意义。
输配电网电压等级、网络结构等方面存在较大差异，管理调控分属不同控制中心。因此输配电网分布式计算需要分解协调计算。目前关于输配电网分布式无功优化有一定研究。例如：《大电网自动电压控制技术的研究与发展》(电力科学与技术学报，2007年第22卷第1期第7页)提出采用关口电压或无功功率作为协调变量的分层自动电压控制(AVC)，但是此种方法不适用于大量DG接入的输配网计算。《输—配电网联合静态电压安全评估与无功优化》(山东大学，2015年)和《输配电网潮流与优化的理论研究》(山东大学，2013年)采用主从分裂法在输配电网进行了考虑双方潮流影响的协同无功优化，模型中将输配网直接分解，难以保证最优解。《Coordinated transmission and distribution AC optimal powerflow》(IEEE Transactions on Smart Grid，2016年第PP卷第99期第1页)提出广义主从算法将输配电网分解并进行分布式有功优化。该方案利用交换边界变量和乘子的方法可以一定程度上精确的全局最优解，但是边界乘子的互换对输配电网的求解算法有一定要求，对离散变量的处理也带来了一定难度。《Decentralized reactive power optimizationmethod for transmission and distribution networks accommodating large-scaleDG integration》(IEEE Transactions on Sustainable Energy，2017年第1卷第8期第363页)采用广义benders分解法将输配电网无功优化问题进行分解并计算。
发明内容
本发明所要解决的技术问题是：提供一种基于分布式计算的输配电网一体化无功优化方法，结合KKT最优解条件、广义主从分裂法、异步迭代分布式计算和无功优化技术对输配电网组成的全局电力系统进行无功优化。
本发明为解决上述技术问题采用以下技术方案：
一种基于分布式计算的输配电网一体化无功优化方法，所述输配电网包括一个输电网和多个配电网，各配电网和输电网之间通过各自对应的边界节点联系；包括如下步骤：
步骤1，对第k个配电网边界电压幅值和相角赋初值，根据所赋初值对输配电网进行一体化分布式潮流计算，得到潮流结果并判断潮流结果是否收敛；若不收敛则停止计算，否则输出一体化潮流边界节点状态变量
Figure GDA0002354697940000021
Figure GDA0002354697940000022
k＝1,2,…,NDis，NDis为配电网的数量，开始输配电网一体化分布式无功优化计算；
步骤2，利用广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型；迭代次数赋初值m＝1，第k个配电网边界影响因子赋初值
Figure GDA0002354697940000023
对应的输电网边界影响因子赋初值
Figure GDA0002354697940000024
第k个配电网边界电压幅值、相角分别为
Figure GDA0002354697940000025
步骤3，根据步骤2分解得到的第k个配电网无功优化模型，采用二次罚函数-原对偶内点法进行无功优化计算，第k个配电网边界电压幅值、相角分别为
Figure GDA0002354697940000026
输电网边界影响因子取
Figure GDA0002354697940000027
优化计算得到边界功率
Figure GDA0002354697940000028
及配电网边界影响因子
Figure GDA0002354697940000029
并发送给输电网作为迭代数据；
步骤4，根据步骤2分解得到的输电网无功优化模型，采用二次罚函数-原对偶内点法进行无功优化计算，输电网中与各配电网对应的边界节点的负荷功率值取
Figure GDA0002354697940000031
配电网边界影响因子取
Figure GDA0002354697940000032
优化计算得到输电网中与各配电网对应的边界节点的电压幅值
Figure GDA0002354697940000033
相角
Figure GDA0002354697940000034
和输电网边界影响因子
Figure GDA0002354697940000035
并发送给配电网作为迭代数据；
步骤5，步骤4完成后，判断边界节点的状态变量和影响因子是否收敛，若收敛则计算成功，否则m＝m+1，将步骤4计算得到的电压幅值、相角作为下一次迭代时各个配电网边界节点的电压幅值、相角，即
Figure GDA0002354697940000036
并返回步骤3。
作为本发明的一种优选方案，步骤2所述利用广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型，其中输配电网分解前数学模型为：
输配电网全局无功优化目标函数为：
Figure GDA0002354697940000037
其中，NDis为配电网的数量，xM表示输电网内部变量向量，
Figure GDA0002354697940000038
xpcc,k分别表示第k个配电网内部变量向量、边界变量向量，fM、
Figure GDA0002354697940000039
分别表示输电网、第k个配电网优化的目标函数；
输电网内部等式与不等式约束为：
Figure GDA00023546979400000310
gM(xM)≤0
其中，hM表示输电网潮流约束，gM表示输电网的状态变量与控制变量的不等式约束；
第k个配电网内部等式与不等式约束为：
Figure GDA00023546979400000311
Figure GDA00023546979400000312
其中，
Figure GDA00023546979400000313
表示第k个配电网潮流约束，
Figure GDA00023546979400000314
表示状态变量与控制变量的不等式约束；
输配电网边界约束如下：
Figure GDA00023546979400000410
gpcc,k(xpcc,k)≤0k＝1，…，NDis
其中，hpcc,k为输电网与第k个配电网的边界等式约束条件，gpcc,k为输电网与第k个配电网的边界不等式约束条件。
作为本发明的一种优选方案，步骤2所述基于广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型，其中输电网分解后无功优化模型为：
Figure GDA0002354697940000041
其中，
Figure GDA0002354697940000042
表示分解后的输电网目标函数，ΩM表示输电网所有节点的集合，Vi、Vj分别表示节点i、j的电压幅值，Gij表示节点i、j之间的电导，θij表示节点i、j之间的相角差，NDis为配电网的数量，xpcc,k表示第k个配电网边界变量向量，xpcc,k＝[VM,pcc,k θM,pcc,k]T，VM,pcc,k、θM,pcc,k分别为输电网中与第k个配电网的边界节点的幅值和相角，～表示变量在上述输电网分解后无功优化模型求解过程中为定值，T表示转置，αk具体表示为：
Figure GDA0002354697940000043
其中，
Figure GDA0002354697940000044
表示第k个配电网优化的目标函数，
Figure GDA0002354697940000045
表示第k个配电网潮流约束，
Figure GDA0002354697940000046
表示状态变量与控制变量的不等式约束，
Figure GDA0002354697940000047
分别表示分解后的第k个配电网边界等式约束、不等式约束，
Figure GDA0002354697940000048
分别为第k个配电网边界电压幅值、相角，
Figure GDA0002354697940000049
分别表示第k个配电网的等式约束、不等式约束、边界等式约束、边界不等式约束对应的拉格朗日乘子向量；
输电网的等式与不等式约束为：
Figure GDA0002354697940000051
其中，Pdi、Qdi分别表示输电网节点i的有功、无功负荷，PGi、QGi分别表示节点i发电机的有功、无功出力，QCi表示节点i电容器的无功出力，Bij表示节点i、j之间的电纳；QGi,max、QGi，min分别表示节点i发电机的无功出力上、下限；Tl、Tl，max、Tl，min分别表示支路l变压器的可调变比及其上、下限；Vi,max、Vi，min分别表示节点i电压的上、下限；Ci、Ci，max Ci，min分别表示节点i电容器的档位及其上、下限。
作为本发明的一种优选方案，步骤2所述利用广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型，其中第k个配电网分解后无功优化模型为：
边界不等式约束包括有功、无功功率上下限：
Figure GDA0002354697940000052
Figure GDA0002354697940000053
其中，
Figure GDA0002354697940000054
分别表示边界有功功率及其上、下限，
Figure GDA0002354697940000055
Figure GDA0002354697940000056
分别表示边界无功功率及其上、下限；
配电网分解后的目标函数为：
Figure GDA0002354697940000057
其中，
Figure GDA0002354697940000058
表示分解后的配电网目标函数，
Figure GDA0002354697940000059
表示第k个配电网所有节点的集合，Gab表示节点a、b之间的电导，Va、Vb分别表示节点a、b的电压幅值，θab表示节点a、b之间的相角差，
Figure GDA00023546979400000510
为上次输电网迭代计算结果，表示边界优化计算结果对输网优化目标的影响，在此时为定值；
βk＝yM.pcc,k
Figure GDA0002354697940000066
其中，yM.pcc,k表示输电网中与第k个配网的边界等式约束对应的拉格朗日乘子，
Figure GDA0002354697940000061
分别为第k个配电网边界有功、无功功率；
第k个配电网的等式与不等式约束为：
Figure GDA0002354697940000062
其中，Pd,a、Qd,a分别表示节点a的有功、无功负荷，PDG,a、QDG,a分别表示节点a并网DG的有功、无功出力，QC,a表示节点a处电容器的无功出力，Bab表示节点a、b之间的电纳，QDG,a,max、QDG,a，min分别表示节点a处并网DG的无功出力上、下限，Va,max、Va，min分别表示节点a处电压上、下限；Ca、Ca，max、Ca，min分别表示节点a处电容器的档位及其上、下限。
作为本发明的一种优选方案，步骤5所述判断边界节点的状态变量和影响因子是否收敛，其中收敛的判据为：
Figure GDA0002354697940000063
其中，k＝1,2,…,NDis，NDis为配电网的数量，
Figure GDA0002354697940000064
分别为第k个配电网边界有功、无功功率，
Figure GDA0002354697940000065
分别为第k个配电网边界电压幅值、相角，αk为第k个配电网中的配网边界影响因子，βk为输电网中与第k个配电网对应的输网边界影响因子，m+1、m、m-1均表示迭代次数，ε＝10-5，|·|表示绝对值。
本发明采用以上技术方案与现有技术相比，具有以下技术效果：
1、本发明将全局优化模型分解，引入边界影响因子，保证了输、配电网优化目的的一致性和同步性。
2、本发明校正环节中输、配电网无功优化的分布式计算保持了输、配电网现有的计算模式，通过PCC点处电气信息和边界影响因子的交换即可实现全网无功优化的分布式计算。
3、本发明采用罚函数规整离散变量，使得在输配网全局优化目标函数可微可分的基础上，保证了边界影响因子的计算，保证了方法的可行性和适用度。
附图说明
图1是本发明输配全局无功优化问题分解示意图。
图2是本发明输配电网分布式无功优化算法流程图。
图3是独立优化、集中优化与本发明方法下的电压水平变化趋势图。
图4是独立优化、集中优化与本发明方法下的无功补偿设备投入量图。
具体实施方式
下面详细描述本发明的实施方式，所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
本发明将输配电网无功优化问题分解为输电网无功优化、配电网无功优化和边界信息交互三部分，引入代表输配电网的相互影响的边界拉格朗日乘子项(简称边界影响因子)，并将边界状态变量及边界影响因子作为协同无功优化的边界协调量，采用二次罚函数-原对偶内点法分别交替求解输、配电网无功优化子问题，通过分布式迭代实现输配电网一体化无功优化计算。针对离散变量导致边界乘子无法计算的问题，本发明引入二次罚函数规整离散变量，使得在输配网全局优化目标函数可微可分的基础上，实现离散变量的处理。
如图1所示的输配电网全局无功优化问题分解图，主要包括分解前和分解后两个过程，分解前包括一个输电网(M)连接多个配电网Dk(k＝1,2,…,NDis)，其边界节点简称PCC点。分解后为输电网无功优化、配电网无功优化和边界交互三个部分。具体分解过程为：
输配网分解前数学模型：
1)输配网全局无功优化目标函数为：
Figure GDA0002354697940000081
其中，NDis表示配电网的数量，xM、
Figure GDA0002354697940000082
xpcc.k分别表示输网内部变量向量、配网k内部变量向量及其边界变量向量，fM、
Figure GDA0002354697940000083
分别表示输网、第k个配网优化的目标函数。
2)输电网内部等式与不等式约束为：
Figure GDA0002354697940000084
gM(xM)≤0
其中，hM表示输网潮流约束，gM表示输网的状态变量与控制变量的不等式约束。将PCC点进行复制，规定xM,pcc,k与
Figure GDA0002354697940000085
为电气上的同一点，拥有相同的电气量，
Figure GDA0002354697940000086
分别表示输网和第k个配网的PCC点在输、配网中的状态变量向量。
3)配网k内部等式与不等式约束为：
Figure GDA0002354697940000087
Figure GDA0002354697940000088
其中，
Figure GDA0002354697940000089
为第k个配网潮流约束，
Figure GDA00023546979400000810
表示状态变量与控制变量的不等式约束。
4)输配电网边界约束如下：
Figure GDA00023546979400000811
gpcc,k(xpcc,k)≤0k＝1，…，NDis
其中，hpcc,k为输电网与配网k的边界等式约束条件，gpcc,k为输电网与配网k的边界不等式约束条件。
全网无功优化的KKT条件分解：
1)为将输配电网潮流解耦，引入变量向量Spcc.k，表示输配网边界传输功率向量。则边界等式约束hpcc,k变为：
hM,pcc.k(xM,xM,pcc.k)＝Spcc.k
Figure GDA0002354697940000091
其中，hM,pcc.k、
Figure GDA0002354697940000092
分别表示分解后输网、配网k边界等式约束。
2)输配网全局目标函数根据拉格朗日理论分解，下面进行分解。全网拉格朗日函数为：
Figure GDA0002354697940000093
其中，y、w分别表示等式、不等式约束的对偶乘子，上标T表示转置。全网取得最优解的KKT条件为：
Figure GDA0002354697940000094
设
Figure GDA0002354697940000095
则全网取得最优解的KKT条件也可以写为：
Figure GDA0002354697940000096
若要满足上述条件，进行一阶线性近似，输、配电网目标函数改进为：
Figure GDA0002354697940000097
Figure GDA0002354697940000098
Figure GDA0002354697940000101
βk＝yM.pcc,k
其中，
Figure GDA0002354697940000102
分别表示分解后的输配网目标函数，αk、βk分别称为配网边界影响因子、输网边界影响因子。可以看出，输电网目标函数的增加项αk Txpcc.k代表边界对配电网的影响，配电网的目标函数增加项βk TSpcc.k代表边界对输电网的影响。
因此通过以上处理，本发明将全局无功优化问题分解为输网无功优化子问题、配网无功优化子问题和边界一致性协调子问题。输配网可通过交换少量边界信息和有限次的迭代完成一体化分布式计算，从而达到全局的优化配置。
3)输电网分解后无功优化数学模型
输网通过调节发电机无功出力、电容器分接头及可调变压器的分接头，同时要满足无功电源出力、节点电压幅值的约束和电容器(C)、可调变压器(T)分接头位置的限制。目标函数为：
Figure GDA0002354697940000103
其中，Vi表示节点i的电压幅值；Gij、Bij表示节点i、j之间的电导、电纳；θij＝θi-θj表示节点i、j之间的相角差；ΩM表示输电网所有节点的集合；
Figure GDA0002354697940000104
来自配电网上一次迭代计算结果，此时为定值，具体表示为：
Figure GDA0002354697940000105
xpcc,k＝[VM,pcc,k θM,pcc,k]T
输网的等式与不等式约束为：
Figure GDA0002354697940000111
其中，Pdi、Qdi分别表示输电网节点i的有功、无功负荷，PGi、QGi分别表示节点i发电机的有功、无功出力，QCi表示节点i电容器的无功出力，Bij表示节点i、j之间的电纳；QGi,max、QGi，min分别表示节点i发电机的无功出力上、下限；Tl、Tl，max、Tl，min分别表示支路l变压器的可调变比及其上、下限；Vi,max、Vi，min分别表示节点i电压的上、下限；Ci、Ci，max Ci，min分别表示节点i电容器的档位及其上、下限。
输电网无功优化问题是个连续变量和离散变量共存的、非线性混合整数规划问题。离散变量的存在使边界影响因子的计算面临了困难，因此本发明采用在目标函数中加入二次罚函数项的方法迫使离散控制到达它的一个分级上，同时也可保证乘子项的计算。
4)配网k分解后无功优化数学模型
输配电网边界功率约束，配电网优化过程中，传输功率需符合以下约束，如果超过限制，采用边界值。边界不等式约束包括有功、无功功率上下限：
Figure GDA0002354697940000112
Figure GDA0002354697940000113
其中，
Figure GDA0002354697940000114
分别表示边界有功功率及其上、下限，
Figure GDA0002354697940000115
Figure GDA0002354697940000116
分别表示边界无功功率及其上、下限。配电网分解后的目标函数为：
Figure GDA0002354697940000117
其中，
Figure GDA0002354697940000118
表示分解后的配电网目标函数，
Figure GDA0002354697940000119
表示第k个配电网所有节点的集合，Gab表示节点a、b之间的电导，Va、Vb分别表示节点a、b的电压幅值，θab表示节点a、b之间的相角差，
Figure GDA0002354697940000121
为上次输电网迭代计算结果，表示边界优化计算结果对输网优化目标的影响，在此时为定值。
βk＝yM.pcc,k
Figure GDA0002354697940000122
其中，yM.pcc,k表示输电网中与第k个配网的边界等式约束对应的拉格朗日乘子，
Figure GDA0002354697940000123
分别为第k个配电网边界有功、无功功率；
配网无功优化控制变量包括可控DG的无功出力和电容器：
Figure GDA0002354697940000124
其中，Pd,a、Qd,a分别表示节点a的有功、无功负荷，PDG,a、QDG,a分别表示节点a并网DG的有功、无功出力，QC,a表示节点a处电容器的无功出力，Bab表示节点a、b之间的电纳，QDG,a,max、QDG,a，min分别表示节点a处并网DG的无功出力上、下限，Va,max、Va，min分别表示节点a处电压上、下限；Ca、Ca，max、Ca，min分别表示节点a处电容器的档位及其上、下限。
本发明所涉及的基于分布式计算的输配电网一体化无功优化算法流程图如图2所示，包括如下步骤：
1)配电网k的PCC点赋电压初值(此后默认配网k＝1，…，NDis)，输配电网进行一体化分布式潮流计算。若潮流不收敛，停止计算；若收敛得到一体化潮流PCC点状态变量
Figure GDA0002354697940000125
开始输配电网一体化分布式无功优化计算，同时赋初值：迭代次数m＝1，输、配电网边界影响因子
Figure GDA0002354697940000126
Figure GDA0002354697940000127
配网边界电压幅值和相角
Figure GDA0002354697940000128
2)配电网进行无功优化，配网PCC点电压初值为
Figure GDA0002354697940000129
输网边界影响因子取
Figure GDA00023546979400001210
根据配网k无功优化问题模型，采用的罚函数-原对偶内点法进行无功优化计算，优化计算得边界功率
Figure GDA00023546979400001211
及配网边界影响因子
Figure GDA00023546979400001212
发给输网作为迭代数据。
3)根据输网的无功优化模型，PCC点负荷功率值取
Figure GDA0002354697940000131
Figure GDA0002354697940000132
配电网边界影响因子取
Figure GDA0002354697940000133
采用罚函数-原对偶内点法进行无功优化计算。输网无功优化算得输网PCC点电压幅值和相角
Figure GDA0002354697940000134
以及下次迭代计算的输电网边界影响因子
Figure GDA0002354697940000135
发送给配网作为迭代数据，
Figure GDA0002354697940000136
Figure GDA0002354697940000137
4)判断边界收敛条件是否收敛。若满足则计算成功；否则m＝m+1，返回2)。
其中边界收敛判据为：
Figure GDA0002354697940000138
其中，
Figure GDA0002354697940000139
分别为第k个配电网边界有功、无功功率，
Figure GDA00023546979400001310
分别为第k个配电网边界电压幅值、相角，αk为第k个配电网中的配网边界影响因子，βk为输电网中与第k个配电网对应的输网边界影响因子，m+1、m、m-1均表示迭代次数，ε为很小的数，取ε＝10-5，|·|表示绝对值。
结果验证：为了测试本发明所提方法的有效性，应用本发明方法对IEEE-30和IEEE-33节点系统进行了仿真验证。为不改变原输网负荷分布，将输网PCC点的部分负荷值等效为配网系统，形成输配电网优化系统模型。输电网功率基准值为100MW，配电网电压基准值取为12.66kV，功率基准值取为10000kW，电压安全范围为0.95～1.05p.u.。由于输配电网采用的标幺值不同，因此边界乘子在进行交替迭代时也需要进行标幺值归算。
IEEE30节点系统电容器调节范围标幺值为0～0.3，设为6档，每档0.05，最大、最小离散惩罚系数为100,10；变压器分接头调节范围为0.95～1.05，档位间隙0.01，最大、最小系数为200,50。
IEEE33节点系统在节点4、8、15安装DG(称为DG4,DG8,DG15)，有功出力为300kw，无功出力在-300～300kvar范围内可调；29号节点接入SVC，可调范围为0～300kvar；在10、24号节点接入C，调节范围标幺值为0～0.03，设为6档，每档0.005，最大、最小系数为100,10。
为验证本发明算法正确性，构建输配网集中优化算例系统(输、配网传送数据到同一控制中心优化)，从而找到具有参照性的精确最优解，为接下来的对比分析做准备。将配网网络参数标幺值归算至输网侧，其中为不改变原输网负荷分布，将输网26号节点的部分负荷值等效为配网系统，形成输配集中优化系统算例，输配电网简称M，D26。
采用集中优化、独立优化、本文一体化分布式优化三种方式分别对算例进行仿真，计算结果及对比分析见表1。独立优化表示输配网采用一体化潮流计算的边界结果分别进行无功优化，配电网控制变量初值选为调节范围上限的1/3。
表1集中、独立、一体化分布式无功优化结果对比
Figure GDA0002354697940000141
由表1可知，本发明优化算法和集中优化结果非常接近，证明本发明算法合理有效。相比于独立优化，本发明算法更合理地分配了无功资源，明显减少网损并提高了PCC点电压幅值。输配网可以相互支撑，增加了电网运行的安全性和经济性。给出输配电网各节点电压幅值曲线如图3所示，无功补偿设备投入对比如图4所示。由图3、图4可知，本发明算法可以在降低系统网损的同时，合理分布无功资源。
算例结果表明，本发明可以充分调度输配网无功资源，提供输配网电压水平，配网可为输电网提供无功支撑；输网也能为配网的运行提供支持；输配网一体化分布式无功优化比传统主配网独立优化更加合理，将为大量DG接入后的主配网一体化运行分析提供基础理论和技术支撑。
以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (5)

1.一种基于分布式计算的输配电网一体化无功优化方法，所述输配电网包括一个输电网和多个配电网，各配电网和输电网之间通过各自对应的边界节点联系；其特征在于，包括如下步骤：
步骤1，对第k个配电网边界电压幅值和相角赋初值，根据所赋初值对输配电网进行一体化分布式潮流计算，得到潮流结果并判断潮流结果是否收敛；若不收敛则停止计算，否则输出一体化潮流边界节点状态变量
Figure FDA0002354697930000011
Figure FDA0002354697930000012
NDis为配电网的数量，开始输配电网一体化分布式无功优化计算；
步骤2，利用广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型；迭代次数赋初值m＝1，第k个配电网边界影响因子赋初值
Figure FDA0002354697930000013
对应的输电网边界影响因子赋初值
Figure FDA0002354697930000014
第k个配电网边界电压幅值、相角分别为
Figure FDA0002354697930000015
步骤3，根据步骤2分解得到的第k个配电网无功优化模型，采用二次罚函数-原对偶内点法进行无功优化计算，第k个配电网边界电压幅值、相角分别为
Figure FDA0002354697930000016
输电网边界影响因子取
Figure FDA0002354697930000017
优化计算得到边界功率
Figure FDA0002354697930000018
及配电网边界影响因子
Figure FDA0002354697930000019
并发送给输电网作为迭代数据；
步骤4，根据步骤2分解得到的输电网无功优化模型，采用二次罚函数-原对偶内点法进行无功优化计算，输电网中与各配电网对应的边界节点的负荷功率值取
Figure FDA00023546979300000110
配电网边界影响因子取
Figure FDA00023546979300000111
优化计算得到输电网中与各配电网对应的边界节点的电压幅值
Figure FDA00023546979300000112
相角
Figure FDA00023546979300000113
和输电网边界影响因子
Figure FDA00023546979300000114
并发送给配电网作为迭代数据；
步骤5，步骤4完成后，判断边界节点的状态变量和影响因子是否收敛，若收敛则计算成功，否则m＝m+1，将步骤4计算得到的电压幅值、相角作为下一次迭代时各个配电网边界节点的电压幅值、相角，即
Figure FDA00023546979300000115
并返回步骤3。
2.根据权利要求1所述基于分布式计算的输配电网一体化无功优化方法，其特征在于，步骤2所述利用广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型，其中输配电网分解前数学模型为：
输配电网全局无功优化目标函数为：
Figure FDA0002354697930000021
其中，NDis为配电网的数量，xM表示输电网内部变量向量，
Figure FDA0002354697930000029
xpcc,k分别表示第k个配电网内部变量向量、边界变量向量，fM、
Figure FDA0002354697930000022
分别表示输电网、第k个配电网优化的目标函数；
输电网内部等式与不等式约束为：
Figure FDA0002354697930000023
gM(xM)≤0
其中，hM表示输电网潮流约束，gM表示输电网的状态变量与控制变量的不等式约束；
第k个配电网内部等式与不等式约束为：
Figure FDA0002354697930000024
Figure FDA0002354697930000025
其中，
Figure FDA0002354697930000026
表示第k个配电网潮流约束，
Figure FDA0002354697930000027
表示状态变量与控制变量的不等式约束；
输配电网边界约束如下：
Figure FDA0002354697930000028
gpcc,k(xpcc,k)≤0 k＝1，…，NDis
其中，hpcc,k为输电网与第k个配电网的边界等式约束条件，gpcc,k为输电网与第k个配电网的边界不等式约束条件。
3.根据权利要求1所述基于分布式计算的输配电网一体化无功优化方法，其特征在于，步骤2所述基于广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型，其中输电网分解后无功优化模型为：
Figure FDA0002354697930000031
其中，
Figure FDA0002354697930000032
表示分解后的输电网目标函数，ΩM表示输电网所有节点的集合，Vi、Vj分别表示节点i、j的电压幅值，Gij表示节点i、j之间的电导，θij表示节点i、j之间的相角差，NDis为配电网的数量，xpcc,k表示第k个配电网边界变量向量，xpcc,k＝[VM,pcc,k θM,pcc,k]T，VM,pcc,k、θM,pcc,k分别为输电网中与第k个配电网的边界节点的幅值和相角，～表示变量在上述输电网分解后无功优化模型求解过程中为定值，T表示转置，αk具体表示为：
Figure FDA0002354697930000033
其中，
Figure FDA0002354697930000034
表示第k个配电网优化的目标函数，
Figure FDA0002354697930000035
表示第k个配电网潮流约束，
Figure FDA0002354697930000036
表示状态变量与控制变量的不等式约束，
Figure FDA0002354697930000037
分别表示分解后的第k个配电网边界等式约束、不等式约束，
Figure FDA0002354697930000038
分别为第k个配电网边界电压幅值、相角，
Figure FDA0002354697930000039
分别表示第k个配电网的等式约束、不等式约束、边界等式约束、边界不等式约束对应的拉格朗日乘子向量；
输电网的等式与不等式约束为：
Figure FDA00023546979300000310
其中，Pdi、Qdi分别表示输电网节点i的有功、无功负荷，PGi、QGi分别表示节点i发电机的有功、无功出力，QCi表示节点i电容器的无功出力，Bij表示节点i、j之间的电纳；QGi,max、QGi，min分别表示节点i发电机的无功出力上、下限；Tl、Tl，max、Tl，min分别表示支路l变压器的可调变比及其上、下限；Vi,max、Vi，min分别表示节点i电压的上、下限；Ci、Ci，max Ci，min分别表示节点i电容器的档位及其上、下限。
4.根据权利要求1所述基于分布式计算的输配电网一体化无功优化方法，其特征在于，步骤2所述利用广义主从分裂法将输配电网一体化分布式无功优化计算模型分解为输电网无功优化、第k个配电网无功优化和边界节点信息交互三个模型，其中第k个配电网分解后无功优化模型为：
边界不等式约束包括有功、无功功率上下限：
Figure FDA0002354697930000041
Figure FDA0002354697930000042
其中，分别表示边界有功功率及其上、下限，
Figure FDA0002354697930000044
Figure FDA0002354697930000045
分别表示边界无功功率及其上、下限；
配电网分解后的目标函数为：
Figure FDA0002354697930000046
其中，
Figure FDA0002354697930000047
表示分解后的配电网目标函数，
Figure FDA0002354697930000048
表示第k个配电网所有节点的集合，Gab表示节点a、b之间的电导，Va、Vb分别表示节点a、b的电压幅值，θab表示节点a、b之间的相角差，
Figure FDA0002354697930000049
为上次输电网迭代计算结果，表示边界优化计算结果对输网优化目标的影响，在此时为定值；
βk＝yM.pcc,k
Figure FDA00023546979300000410
其中，yM.pcc,k表示输电网中与第k个配网的边界等式约束对应的拉格朗日乘子，
Figure FDA00023546979300000411
分别为第k个配电网边界有功、无功功率；
第k个配电网的等式与不等式约束为：
Figure FDA0002354697930000051
其中，Pd,a、Qd,a分别表示节点a的有功、无功负荷，PDG,a、QDG,a分别表示节点a并网DG的有功、无功出力，QC,a表示节点a处电容器的无功出力，Bab表示节点a、b之间的电纳，QDG,a,max、QDG,a，min分别表示节点a处并网DG的无功出力上、下限，Va,max、Va，min分别表示节点a处电压上、下限；Ca、Ca，max、Ca，min分别表示节点a处电容器的档位及其上、下限。
5.根据权利要求1所述基于分布式计算的输配电网一体化无功优化方法，其特征在于，步骤5所述判断边界节点的状态变量和影响因子是否收敛，其中收敛的判据为：
Figure FDA0002354697930000052
其中，k＝1,2,…,NDis，NDis为配电网的数量，
Figure FDA0002354697930000053
分别为第k个配电网边界有功、无功功率，
Figure FDA0002354697930000054
分别为第k个配电网边界电压幅值、相角，αk为第k个配电网中的配网边界影响因子，βk为输电网中与第k个配电网对应的输网边界影响因子，m+1、m、m-1均表示迭代次数，ε＝10-5，|·|表示绝对值。
CN201710452090.1A 2017-06-15 2017-06-15 一种基于分布式计算的输配电网一体化无功优化方法 Active CN107171341B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710452090.1A CN107171341B (zh) 2017-06-15 2017-06-15 一种基于分布式计算的输配电网一体化无功优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710452090.1A CN107171341B (zh) 2017-06-15 2017-06-15 一种基于分布式计算的输配电网一体化无功优化方法

Publications (2)

Publication Number Publication Date
CN107171341A CN107171341A (zh) 2017-09-15
CN107171341B true CN107171341B (zh) 2020-03-17

Family

ID=59818580

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710452090.1A Active CN107171341B (zh) 2017-06-15 2017-06-15 一种基于分布式计算的输配电网一体化无功优化方法

Country Status (1)

Country Link
CN (1) CN107171341B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107947180B (zh) * 2017-12-12 2020-06-02 北京四方继保自动化股份有限公司 一种适用于发输配一体化潮流计算的拓扑分析方法
CN108536917A (zh) * 2018-03-15 2018-09-14 河海大学 一种输配电网全局电压稳定控制的分布式计算方法
CN109038532B (zh) * 2018-07-25 2022-07-26 深圳供电局有限公司 一种基于直流潮流等值的分布式输配协同机组组合方法
CN109586313B (zh) * 2018-11-08 2022-04-22 国电南瑞科技股份有限公司 无功优化的高效内点求解方法
CN109829525B (zh) * 2019-01-24 2020-04-07 清华大学 一种基于群智能的建筑控制方法及系统
CN110120673B (zh) * 2019-05-09 2020-10-02 山东大学 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统
CN110071503B (zh) * 2019-05-09 2020-10-02 山东大学 分布式输配协同无功优化的二次规划模型构建方法及系统
CN110611315B (zh) * 2019-10-14 2022-11-22 国网四川省电力公司经济技术研究院 一种考虑输配协同的分布式电源最大接入容量计算方法
CN110880758B (zh) * 2019-11-22 2021-01-01 清华大学 电力系统中输电网与配电网的分解协调最优潮流控制方法
CN111463794B (zh) * 2020-04-29 2021-07-30 山东大学 一种附加对地阻抗的输配电网协同潮流控制方法及系统
CN112308331A (zh) * 2020-11-09 2021-02-02 中国南方电网有限责任公司 输配网分离电力市场下的最优潮流模型及其迭代求解方法
CN115411746B (zh) * 2022-11-03 2023-02-14 江苏金智科技股份有限公司 一种大电网电压无功优化的分解协调系统及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101635456A (zh) * 2009-06-23 2010-01-27 上海申瑞电力科技股份有限公司 输配电网联合状态估计的方法
CN103001218A (zh) * 2012-11-29 2013-03-27 清华大学 基于输配电网综合效益的无功电压协调控制方法
CN104319780A (zh) * 2014-10-29 2015-01-28 国家电网公司 一种输配电网全局无功电压优化方法
CN104701867A (zh) * 2015-03-27 2015-06-10 河海大学 基于分支定界法以及原对偶内点法的日前无功优化方法
CN105048468A (zh) * 2015-07-27 2015-11-11 河海大学 基于分布式计算的输配电网一体化电压稳定评估方法
CN105098797A (zh) * 2015-08-25 2015-11-25 中国电力科学研究院 一种基于网络分解与区域协调的配电网无功优化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101635456A (zh) * 2009-06-23 2010-01-27 上海申瑞电力科技股份有限公司 输配电网联合状态估计的方法
CN103001218A (zh) * 2012-11-29 2013-03-27 清华大学 基于输配电网综合效益的无功电压协调控制方法
CN104319780A (zh) * 2014-10-29 2015-01-28 国家电网公司 一种输配电网全局无功电压优化方法
CN104701867A (zh) * 2015-03-27 2015-06-10 河海大学 基于分支定界法以及原对偶内点法的日前无功优化方法
CN105048468A (zh) * 2015-07-27 2015-11-11 河海大学 基于分布式计算的输配电网一体化电压稳定评估方法
CN105098797A (zh) * 2015-08-25 2015-11-25 中国电力科学研究院 一种基于网络分解与区域协调的配电网无功优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Master–Slave-Splitting Based Distributed Global Power Flow Method for Integrated Transmission and Distribution Analysis;Hongbin Sun et al.;《IEEE Transactions on Smart Grid》;20150531;第6卷(第3期);全文 *

Also Published As

Publication number Publication date
CN107171341A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN107171341B (zh) 一种基于分布式计算的输配电网一体化无功优化方法
Zhao et al. Coordinated restoration of transmission and distribution system using decentralized scheme
Pandey et al. A literature survey on load–frequency control for conventional and distribution generation power systems
Chai et al. Hierarchical distributed voltage optimization method for HV and MV distribution networks
CN110690732B (zh) 一种光伏无功分区计价的配电网无功优化方法
CN105119280A (zh) 基于锥优化的交直流混合结构有源配电网运行优化方法
CN102684201B (zh) 一种基于电压越限概率的含风电场电网无功优化方法
CN114362267B (zh) 考虑多目标优化的交直流混合配电网分散式协调优化方法
Xu et al. A hierarchically coordinated operation and control scheme for DC microgrid clusters under uncertainty
Aziz et al. Admm-based distributed optimization of hybrid mtdc-ac grid for determining smooth operation point
CN103490428B (zh) 微电网无功补偿容量配置方法及系统
Dong et al. A decentralized optimal operation of AC/DC hybrid microgrids equipped with power electronic transformer
Zhai et al. Distributed optimal power flow for VSC-MTDC meshed AC/DC grids using ALADIN
CN112134309A (zh) 一种适用于配电网分布式电压控制的新型分区方法
Mallikarjunaswamy et al. Implementation of an effective hybrid model for islanded microgrid energy management
CN115333110A (zh) 基于admm的配电网-微电网群协同分布式优化调度方法及系统
CN112467748A (zh) 三相不平衡主动配电网双时标分布式电压控制方法及系统
Zou et al. Optimized Robust Controller Design Based on CPSOGSA Optimization Algorithm and H 2/H∞ Weights Distribution Method for Load Frequency Control of Micro-Grid
Sivadanam et al. Dynamic performance enhancement of interconnected hybrid thermal power system in the presence of electric vehicles
Dong et al. A distributed power transfer limit calculation method for multi-area interconnection power networks
Wei et al. Advanced Control and Optimization for Complex Energy Systems.
Gündüz et al. The impact of electric vehicles aggregator on the stability region of micro-grid system with communication time delay
CN114709848A (zh) 一种考虑电能质量补偿的多相配电网储能优化调度方法
Shrivastwa et al. Frequency control using V2G and synchronous power controller based HVDC links in presence of wind and PV units
Shayeghi et al. Management of Autonomous Microgrids Using Multi-Agent Based Online Optimized NF-PID Controller

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