CN105406459B - 利用基于一致性分布式处理估计电网中的最优潮流的方法 - Google Patents

利用基于一致性分布式处理估计电网中的最优潮流的方法 Download PDF

Info

Publication number
CN105406459B
CN105406459B CN201510570334.7A CN201510570334A CN105406459B CN 105406459 B CN105406459 B CN 105406459B CN 201510570334 A CN201510570334 A CN 201510570334A CN 105406459 B CN105406459 B CN 105406459B
Authority
CN
China
Prior art keywords
bus
virtual
spirte
variable
agency
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
CN201510570334.7A
Other languages
English (en)
Other versions
CN105406459A (zh
Inventor
M·本诺斯曼
刘杰
A·拉格胡娜汉
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Publication of CN105406459A publication Critical patent/CN105406459A/zh
Application granted granted Critical
Publication of CN105406459B publication Critical patent/CN105406459B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
    • H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F13/00—Interconnection of, or transfer of information or other signals between, memories, input/output devices or central processing units
    • 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
    • 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
    • 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
    • Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
    • 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
    • Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20—End-user application control systems
    • Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
    • 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)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)

Abstract

利用基于一致性分布式处理估计电网中的最优潮流的方法。一种估计电网中的最优潮流OPF的方法，所述电网被表示成图形，所述图形被分割为虚拟子图形，各个虚拟子图形包括至少一个母线并且与代理关联，所述代理测量局部变量并且更新一致性变量CV。利用代理交换并更新相邻虚拟子图形的一致性变量。利用代理基于CV和局部变量针对虚拟子图形求解OPF问题。迭代所述交换和所述求解步骤，直至满足终止条件为止，此时针对各个虚拟子图形输出最优OPF。

Description

利用基于一致性分布式处理估计电网中的最优潮流的方法
技术领域
本发明总体上涉及电网，更具体地讲，涉及利用分布式处理估计电网中的潮流。
背景技术
电力网经由传输线和变电站从发电机向消费者(负载)分配电力。最优潮流(OPF)是在电网运行中要求解的关键问题。在发电限制、电压限制和传输线路热限制的约束下，OPF问题的解估计有功和无功发电以及电网中的母线的电压以使发电成本最小化。OPF对电网的可靠运行很关键。
OPF难以求解，因为该问题是非凸和非线性的。通常，由于相邻母线的电压之间的二次关系，非凸性难以确保全局最优解。OPF是在发电、所有母线中的电压以及传输线路热限制方面具有大量决策变量和约束的大规模优化问题。对于实际实现，需要降低求解该问题的计算复杂性。
实践中，通过用DC潮流问题(线性规划问题)近似AC潮流问题来避免求解OPF的困难。这种近似具有输电网可接受的精度水平，并且被许多独立系统运营商(ISO)使用。然而，DC潮流解对于配电网而言不够精确，无法满足智能电网的要求(包括可再生能源、分布式发电和存储)。
另一方法基于半正定规划(SDP)松弛使用AC OPF问题的全局最优解。即使在多个情况下可应对该问题的非凸性，但是集中式方法仍无法满足现代电网的要求。由于连接到电网的不同的、时变的并且易失性的负载和存储，集中式OPF问题的解应当反映电网的任何显著改变。
为了满足电网运行的实时要求，通常诸如区域输电运营商(RTO)和ISO的发电商每五分钟求解一次OPF问题。结果，大规模OPF问题需要集中式方法无法提供的效率和精度。
因此，大规模OPF问题的可扩展、快速求解使用分布式方法，其中OPF问题被分割为多个小规模子问题。各个子问题通过具有代理至代理通信能力的单个代理作为计算实体来求解。代理根据通信协议交换数据。因此，所有代理合作地以分布式方式参与求解OPF问题。
在一个分布式方法中，电网被分割为极大团，参见Lam等人的“Distributedalgorithms for optimal power flow problem”(IEEE第51届决策和控制(CDC)年会，第430–437页，2012年)。该方法利用松弛的半正定规划(SDP)求解OPF问题，以凸化(即，近似)该问题。在问题的SDP松弛之后，该方法使用交替方向乘子法(ADMM)来以分布式方式求解凸化问题。凸化方法在收敛时无法确保提供最优潮流问题的可行点。
在另一方法中，各个母线被当作单独的代理，其维持相邻母线的电压的估计，并且通过使发电的局部成本和估计误差最小化来构造约化OPF问题，参见Dall’Anese等人的“Distributed optimal power flow for smart microgrids”(IEEE智能电网会刊，第4卷，第3期，第1464–1475页，2013年)。对各个代理的局部等式约束是各个母线的功率平衡方程，其中其它母线的电压变量由估计取代。当方法达到可行解时，所述估计收敛于真值。该方法还使用SDP松弛，然后使用对偶算法来进行分布式优化，其也无法确保OPF问题的可行点。
另外，所述方法仅考虑电网作为树建模的情况，并且那两种方法均使用极大团分解。在这两种方法中，用于表示电网的图形(graph)限于树结构，所述图形被分解为极大团。早前的图形和团中的节点和边缘仅表示电网的真实物理组件(例如，实际母线和传输线路)。
发明内容
本发明的实施方式提供了一种利用基于一致性(consensus)的分布式处理来估计电网中的最优潮流(OPF)的方法。该方法由代理执行，其中各个代理与电网的至少一条母线关联。各个代理与局部变量和一致性变量关联，所述一致性变量是各个局部变量的估计。
局部变量表示：其电压和功率变量；以及其邻近代理的电压和功率变量的估计。局部变量由代理利用局部优化问题来计算，所述局部优化问题使母线处的发电成本以及局部变量相对于其一致性变量的偏差最小化。所述优化问题受到局部功率平衡约束，其表示连接代理与其邻居的线路上的潮流。
各个代理向其邻近代理通信：(i)其电压和功率以及(ii)其邻近代理的电压和功率的一致性变量。各个代理应用一致性滤波器以更新其一致性变量。优化和通信的处理被重复，直至各个代理的局部变量和一致性变量收敛于容限内。
所述方法将电网表示成图形，所述图形被分割为虚拟子图形。与现有技术不同，子图形的顶点(也被称作节点)无需表示真实物理组件(即，与真实物理发电机关联的母线和真实负载)。也就是说，所述分割步骤可从真实组件生成虚拟组件，即，虚拟发电机或虚拟负载。这有利于求解例如单个发电机连接到多个负载等时的OPF问题(其无法通过其它方法来处理)。
所述方法可使用两个不同的实施方式。在一个实施方式中，各个虚拟子图形仅包括与一个代理关联的一个母线。所述代理共享电压并且估计通过传输线路连接的邻近母线的电压。
在另一实施方式中，虚拟子图形包括至少两个母线，各个虚拟子图形与一个代理关联。仅针对连接虚拟子图形的母线估计潮流。
附图说明
图1A是根据本发明的实施方式的电网中的母线和关联的代理的框图；
图1B是根据本发明的实施方式的利用基于一致性的分布式处理来估计电网中的最优潮流(OPF)的方法的框图；
图1C是根据本发明的实施方式的具有三个发电机并且分割为两个虚拟子图形的网络的框图；
图2是根据本发明的实施方式的具有两个发电机和一个负载网络并且利用虚拟发电机和虚拟负载分割为两个虚拟子图形的网络的框图；
图3是根据本发明的实施方式的具有两个发电机和一个负载网络并且利用两个虚拟负载分割为两个虚拟子图形的网络的框图；
图4是根据本发明的实施方式的具有一个发电机和两个负载并且利用虚拟负载和虚拟发电机分割为两个虚拟子图形的网络的框图；
图5是根据本发明的实施方式的求解分散式OPF问题的一般情况的框图，其利用虚拟子图形方法将各个虚拟子图形与局部子OPF问题关联；以及
图6是根据本发明的实施方式的求解分散式OPF问题的具体情况的框图，其利用直接方法将各个母线与局部子OPF问题关联。
具体实施方式
本发明的实施方式提供了一种利用基于一致性的分布式处理用于电网中的最优潮流(OPF)的方法。
电网表示
如图1B所示，我们将电网建模成图形50，其中顶点集包括电网中的所有母线，顶点之间的链路ε是连接母线的高架或地下传输线路。
各个母线可连接到发电机或负载。连接到发电机的母线集合被表示成为了方便，连接到发电机的母线被称为发电机母线。没有任何发电机的母线被称为负载母线。我们使用i～j来表示母线i连接到母线j。
母线i处的有功和无功发电功率被表示成和母线i处的功率需求被表示成母线i处的复数电压被表示成Vi＝ei+jfi，其中ei是电压的实部(Re)，fi是虚部(Im)，为了简单，我们将ei和fi统称为Vi。我们将母线i的接地导纳表示成yii。我们将母线i和j之间的线路导纳表示成yij。导纳是复数yij＝gij+jbij。如果则yij＝0。
最优潮流(OPF)问题
所有母线的OPF问题可被表示成
使得
以及 (4)
其中，成本函数Fi可被表示成二次方程
其中c1，i、c2，i是预定正常数。
方程(2)和(3)是母线i的真实无功功率平衡方程。方程(4)中的不等式是电压大小的上下界。方程(5)中的不等式是真实无功发电功率的界。
OPF的非凸性来自功率平衡方程(2)和(3)中的电压之间的二次关系。这使得难以求解OPF问题。另外，电网中的母线的数量对于以集中方式求解OPF而言常常过大。
因此，我们在方程(1)中描述求解OPF问题的分布式方法。我们的方法具有两个变体。在一个实施方式中，各个母线与一个代理关联。在另一实施方式中，电网的图形表示被分割为虚拟子图形，各个虚拟子图形具有一个母线并且与一个代理关联。然而，我们需要指出的是，与现有技术不同，我们仍直接求解非凸问题。另外，代替使用ADMM来求解该问题，我们使用包括惩罚的规划。
基于一致性的分布式潮流估计的直接方法
如图1A所示，对于一个实施方式，电网中的各个母线10与一个代理15关联。在下述的另一实施方式中，代理与两个或更多母线关联。
各个代理包括计算和通信能力。例如，代理可包括连接到存储器21、输入/输出(I/O)接口22和收发器23的处理器20。该处理器可以是多核处理器、微处理器、并行处理器等。
图1B示出利用基于一致性的分布式处理来估计电网中的最优潮流(OPF)的一般方法的框图。
由图形50表示的电网被分割(55)成虚拟子图形60。各个虚拟子图形至少包括母线10。代理15与虚拟子图形60关联(70)。图形和子图形可被存储在存储器中。
代理测量(80)局部变量82并且利用处理器和I/O接口来更新(85)虚拟子图形中的一致性变量。
局部变量包括电压和功率变量。如果代理母线的母线以及与所述母线相邻的那些母线，则局部变量的一致性变量也由代理存储。局部变量受到虚拟子图形中的功率平衡约束。利用惩罚函数客观地对局部变量相对于一致性变量的偏差进行惩罚。所述变量可被存储在存储器中。
代理利用收发器在虚拟子图形60之间交换(85)局部变量。通过应用一致性滤波器86来更新一致性变量。
利用代理基于局部变量和更新的一致性变量求解(90)子图形的OPF问题。交换、更新和求解迭代，直至满足终止条件，输出虚拟子图形的最优OPF 95。
为了描述简单，我们不在代理15与母线10之间进行区分，除非特别需要如此做。
我们将OPF问题规划为多个子问题，各个子问题由对应代理求解。为了简单，我们假设所有母线均连接到发电机，即，对于各个母线i，我们规划OPF的局部子问题。母线i处的局部决策变量为需要指出的是，如方程(2)和(3)中所表示的，母线i与相邻母线j之间存在功率平衡方程中的电压的约束耦合。因此，在求解最优问题时，母线i估计相邻母线j的电压，其被表示为ej(i)，fj(i)。电压Vj(i)＝ej(i)+jfj(i)表示由母线i估计的母线j～i的复数电压。所述估计用于取代母线i的功率平衡方程中的母线j的真实电压。那些估计也是母线i的决策变量的部分。
局部OPF问题
母线i的局部OPF问题可被如下表示。
以及 (11)
其中方程(12)中的约束确保在母线i处形成的母线j的电压的估计与真值Vj一致。
为了求解方程(6)中的局部OPF问题，各个母线i使发电成本最小化，并且估计的ej(i),fj(i)遵循真实电压ej，fj，即，满足方程(12)中的约束。发电成本被包括在方程(6)中的局部OPF问题中。
估计的一致性方法
对于各对真实局部电压ej(i)、fj(i)，母线i维持一对一致性变量并且通过应用以下一致性滤波器来利用从母线j获得的ej，fj的真实局部值更新这些变量：
并且
其中，0＜γ＜1是一致性增益。
即，在(k+1)迭代更新的一致性变量的值的估计是(k)迭代的先前值的估计与(先前真实值与先前估计值之间)的差乘以一致性增益之和。
在我们的方法中，估计(例如，ej(i)和fj(i))未直接被来自其它母线的真值(例如，ej，fj)取代。相反，所述值通过一致性滤波器。一致性变量在该滤波器中可被视为中间变量，并且用于在局部优化中更新估计。
母线i和j在各个迭代k交换真实电压及其一致性变量。换言之，母线i接收真实电压ej(k)，fj(k)和一致性变量并且与相邻母线交换电压ei(k)，fi(k)和一致性变量由于交换仅发生在相邻母线之间，所以与需要获取来自每一个母线的所有数据并且由单个处理中心处理的集中式方法相比，通信开销的额外负担很小。
局部优化问题
为了使发电成本和估计误差二者最小化，重新规划方程(6)中的局部优化问题。对于母线局部优化问题被表示为
以及 (17)
其中ρi＞0是以较高优先级确保估计误差被最小化的大的正惩罚因子。这是必要的，因为在没有这种惩罚的情况下，各个母线可基于有偏估计集合来自私地使它自己的发电成本最小化。
方程(14)中的问题的决策变量仅与母线i和相邻母线有关，因此该问题是具有相对较小规模的局部非线性优化问题。可利用非线性规划求解器来有效地求解该问题，例如内点优化器(IPOPT)或者“寻找约束非线性多元函数的最小值”(fmincon)。
我们可将母线i的分布式估计方法总结如下。
(1)初始化变量：
ei(0)，fi(0)，ej(i)(0)，fj(i)(0)，
(2)在迭代k+1，如果满足终止条件，则结束。否则，各个代理与相邻母线j～i交换变量ei(k)，fi(k)，
(3)如方程(13)中那样更新一致性变量
(4)通过求解方程(14)中的优化问题来更新决策变量，即，
然后，进行到步骤(2)。
在另一实施方式中，局部优化问题的目标函数中的惩罚可利用1范数来规划成
此方法非常适用于求解所有母线均为发电机母线的情况。然而，对于具有负载母线的电网，该方法可对初始值敏感，并且常常对负载母线具有稳态估计误差。此问题是由于负载母线中缺少自由度导致的。为了解决此问题，我们对虚拟子图形使用一致性滤波。
基于虚拟子图形的类一致性分布式优化方法
此实施方式求解电网中存在负载母线时的问题。我们将电网分割为虚拟子图形。在一个实施方式中，各个虚拟子图形与两个母线之间的潮流有关。虚拟子图形与代理关联，该代理按照与上述直接方法相似的方式执行分布式方法。
虚拟子图形
图1C示出图形101，其具有分别通过线105和106连接的母线i 111、母线j和母线j′。分割55将图形101分成虚拟子图形102和103。即，真实母线i被分成两个虚拟母线，例如，母线i1和母线i2121。虚拟组件被示出为虚线。如果母线i111是发电机母线，母线j是负载母线，则虚拟子图形与从母线i至母线j的潮流有关，其中母线i与母线j之间的潮流以及母线i与母线j′之间的潮流根据两个虚拟子图形102和103进行连接。
存在三种类型的虚拟子图形：
情况1：一个发电机母线和一个负载母线；
情况2：两个发电机母线；以及
情况3：两个负载母线。
(1)一个发电机母线和一个负载母线
以及 (22)
其中Ni是母线i的相邻母线的数量。需要指出的是，功率项和电压Vi|(i，j)，Vj|(i，j)是由虚拟子图形维持的局部变量。具体地讲，电压Vi|(i，j)是由虚拟子图形(i，j)估计的母线i的电压。该电压仅是估计，因为另一虚拟子图形也可包括母线i并且确定它自己的估计。为了确保相同变量的一致性估计，对变量强加以下条件。
以及 (25)
功率的局部变量(例如，)可被解释为在虚拟子图形(i，j)的分支流中所涉及的那部分功率。换言之，我们将真实发电或负载分割为多个部分，并且将各个部分指派给不同的虚拟子图形。因此，如方程(24)中的条件所强制的，那些局部变量(例如，发电功率)之和必须等于母线i处的发电的真实值。
(2)两个发电机母线
发电机母线也可具有直接附接的负载。因此，灵活的是将发电机母线建模成虚拟子图形中的负载母线或发电机母线。建模可针对各种情况而不同。在具有彼此连接的两个发电机母线的虚拟子图形的情况下，考虑两种情况。
在第一种情况中，发电机母线没有足够的功率来支持连接的负载或者相邻虚拟子图形中的负载母线。在这种情况下，该发电机母线可被建模成具有供应所需功率的另一发电机母线(其也被建模成包括所述负载的虚拟子图形中的发电机母线)的虚拟子图形中的负载母线。
对于两个发电机母线i1和i2，我们可假设潮流从母线i1至母线i2。因此，我们将母线i2建模成负载母线。流方程可被表示成
需要指出的是，母线i2的负载可被视为消耗来自母线i2的潮流的虚拟负载。这意味着的值可超过真实负载这发生在电网超过母线j的部分可能需要通过母线i2供应的母线i1的功率时。为了一致性，母线i2还必须在其它虚拟子图形中作为发电机母线出现。
图2示出图形210，其表示包括发电机母线i1201、发电机母线i2202和负载母线j203的电网。我们使用圆圈来表示发电机母线，正方形表示负载母线。图形210被分割(55)成两个虚拟子图形211和212。然后，在虚拟子图形212中，则看到母线i2连接到虚拟负载221，真实无功功率被定义为在虚拟子图形212中，则看到母线i2连接到虚拟发电机222。
因此，在虚拟子图形(i2，j)212中，存在与虚拟发电机222关联的虚拟发电功率项使得并且
图2示出具有两个发电机母线的虚拟子图形的示例。在第一种情况下，虚拟子图形(i2，j)的功率平衡方程可被表示成
图3上所示的第二种情况是两个发电机母线共享母线之一的负载，使得母线i1301和母线i2302向母线j 303处的负载供应功率。在这种情况下，虚拟子图形可通过将母线i2302视为虚拟子图形(i1，i2)311中的负载母线来形成，其中母线i2321处的负载部分被定义为使得剩余负载可出现在其它虚拟子图形中，其中母线i2被视为发电机母线324。
(3)两个负载母线
对于包含两个负载母线j1402和j2403的虚拟子图形，潮流从母线j1至母线j2。我们假设存在连接到母线j1422的虚拟发电机，其供应功率我们可将潮流方程写为
需要指出的是虚拟发电项作为虚拟负载402出现在连接到母线j1的发电机母线401中。这可维持电网的功率平衡。
图4示出连接到两个负载母线j1402和j2403的发电机母线401的图形。该图形被分割为虚拟子图形411和412，其中负载母线j1402现在连接到发电机母线i1，具有真实无功功率和的虚拟发电机母线422连接到负载母线i2。因此，在虚拟子图形412中，虚拟发电功率项被增加到母线j1使得并且
具体地讲，虚拟子图形(i，j1)的功率平衡方程可被表示成
流方程之和等同于功率平衡方程
在电网中，对于各个母线i，假定方程(24-26)中对局部变量的条件有效，则包含母线i的所有虚拟子图形(例如，(i，j)和(j，i))中的流方程之和等同于母线i的功率平衡方程。
考虑简单的示例。假定所有虚拟子图形包含发电机母线和负载母线。假定所有发电机母线充当虚拟子图形中的发电机母线。对于各个发电机母线i，我们可在所有相邻母线j上对发电机母线的流方程求和。
通过方程(24-26)中的条件，我们知道与母线i有关的所有虚拟子图形上的发电功率之和等于母线i的原始发电功率，所有其它局部变量也是如此。因此，我们可恢复功率平衡方程如下。
以及 (45)
对于虚拟子图形中的负载母线j，可通过对出现在所有虚拟子图形(i，j)中的流方程求和，使得i～j来恢复功率平衡方程。方程之和被表示为
对于具有两个负载母线的虚拟子图形，这两个负载母线中的一个具有额外虚拟发电功率，此额外功率被增加到具有真实发电机母线的虚拟子图形的项抵消。考虑(j1，j2)均为负载母线的虚拟子图形中的母线j1。母线j1也在虚拟子图形(i，j1)中，其中母线i是发电机母线。然后，虚拟子图形(j1，j2)中的母线j1的流方程被表示成
因此，在虚拟子图形(i，j1)中，母线j1的流方程为
由于我们强制所以当关于母线j1的所有流方程被加起来时，项被抵消。这确保了流方程之和等同于负载的功率平衡方程。
虚拟子图形方法中的局部估计的一致性方法
在分布式优化中各个虚拟子图形是独立的代理。如上所述利用一致性滤波器来执行一致性变量的交换和更新。需要指出的是，在此实施方式中，各个虚拟子图形(i，j)具有母线i和j的电压值。然而，由于母线i可与不止一个虚拟子图形关联，所以电压vi|(i，j)的值必须与其它母线处的电压相容。也对发电功率和负载功率的值强加相同的一致性要求。那些关系由方程(24-26)来强加，其对获得全局问题的可行解而言是必要的。
为了满足这些要求，各个虚拟子图形(i，j)维持母线i的电压、发电和其它虚拟子图形(i，j′)，j′≠j，j′～i的负载的估计。因此，由虚拟子图形(i，j)估计的虚拟子图形(i，j′)中的i的发电功率被表示为负载估计为
对于与虚拟子图形(i，j)共享母线i的虚拟子图形，各个估计的一致性变量为根据以下规则来更新一致性变量
以及
其中常数0＜γ＜1是一致性增益，其可被适当地选择以确保一致性算法的收敛。一致性处理与其它虚拟子图形通信以交换一致性变量。
局部优化问题
类似于直接方法，针对各个虚拟子图形规划局部优化问题以使发电成本和估计误差最小化。在不损失一般性的情况下，我们考虑虚拟子图形(i，j)，其具有还出现在其它虚拟子图形中的母线i和j。局部优化问题可被表示成以下目标函数：
各个虚拟子图形内的局部优化问题针对两个母线。各个虚拟子图形的相邻虚拟子图形的数量等于虚拟子图形内的相邻母线i和j的数量。因此，这是小规模OPF问题，可通过类似fmincon和IPOPT的可用的非线性优化求解器来有效地求解。
(1)初始化fk|(i，j)(0)(k＝i，j)以及估计变量 (k＝i，j)和一致性变量的值。
(2)在迭代k+1，如果达到终止条件(例如，最大迭代次数或者估计的误差容限)，则方法终止。如果没有，则各个代理将局部变量发送给相邻虚拟子图形。
(3)根据(53)更新一致性变量
(4)通过求解优化问题(54)来更新决策变量，然后进行到步骤(2)。
图5是根据本发明的实施方式的利用直接方法求解将各个母线与局部OPF子问题关联的分散式OPF问题的一般情况的框图。OPF问题600被分割(605)成子OPF问题610，每个虚拟子图形1、2、...、n-1、n一个。然后，对于各个虚拟子图形的各个子问题610，以下步骤迭代，直至达到终止条件。在各个虚拟子图形处应用一致性滤波器620以更新一致性变量和针对各个虚拟子图形求解(630)子OPF问题以更新决策变量然后，在终止处，输出(640)各个虚拟子图形的最优电压和功率。
图6是根据本发明的实施方式的利用直接方法求解将各个母线与局部OPF子问题关联的分散式OPF问题的具体情况的框图。OPF问题700被分割(705)为子OPF问题710，每个母线1、2、...、n-1、n一个。然后，对于各个子问题，以下步骤迭代，直至达到终止条件。应用一致性滤波器720以更新一致性变量和求解(730)子OPF问题以更新局部决策变量然后，在终止处，输出(740)各个母线的最优电压和功率。
本发明的效果
实施方式提供了一种基于一致性的分布式优化方法。该方法可用于求解最优潮流问题。该方法以分布式方式在彼此通信的多个代理中执行。描述两个实现方式。在直接方法中，各个母线与代理功率。各个代理在关联的母线的功率平衡方程的约束下求解局部OPF问题。代理估计相邻代理的电压并且与其它代理交换一致性变量以达成一致性。使用一致性滤波方法来确保估计的收敛。在虚拟子图形方法中，各个虚拟子图形与相邻母线之间的分支流有关。此方法在多个负载连接到发电机母线的情况下更有效。针对所有变量达成一致性。
尽管通过优选实施方式的示例描述了本发明，将理解，在本发明的精神和范围内可进行各种其它更改和修改。因此，所附权利要求书的目的在于涵盖落入本发明的真实精神和范围内的所有这些变化和修改。

Claims (14)

1.一种用于估计电网中的最优潮流OPF的方法，其中，所述电网包括通过母线连接的发电机和负载，所述方法包括以下步骤：
将所述电网表示成图形，该图形被分割为虚拟子图形，其中，各个虚拟子图形包括至少一个发电机、负载和一个母线；
将代理与各个虚拟子图形关联，其中，所述代理包括计算和通信能力；
由所述代理在各个虚拟子图形中测量局部变量并且获得一致性变量CV，其中，所述局部变量和所述一致性变量包括受到所述虚拟子图形中的功率平衡约束的电压变量和功率变量；
利用所述代理交换和更新相邻虚拟子图形的所述CV；
利用所述代理，基于所测量的局部变量和所交换的一致性变量来求解各个虚拟子图形的OPF问题；
迭代进行所述交换步骤和所述求解步骤，直至满足终止条件为止；
输出各个虚拟子图形的最优OPF；以及
利用惩罚函数来惩罚所述局部变量与所述一致性变量之间的偏差，所述惩罚函数利用1范数来构造。
2.根据权利要求1所述的方法，其中，各个代理包括计算和通信能力以及存储器。
3.根据权利要求2所述的方法，该方法还包括以下步骤：
将所述图形和子图形存储在所述存储器中。
4.根据权利要求1所述的方法，其中，所述子图形中的顶点表示虚拟发电机和虚拟负载。
5.根据权利要求1所述的方法，其中，所述OPF问题是
其中，母线i处的有功发电功率和无功发电功率为和母线i的复数电压为Vi＝ei+jfi，其中ei是所述电压的实部，fi为虚部，并且ei，fi是局部变量。
6.根据权利要求1所述的方法，其中，所述OPF问题是由相邻母线的电压之间的二次关系导致的非凸。
7.根据权利要求1所述的方法，其中，所述更新步骤使用一致性滤波器。
8.根据权利要求1所述的方法，其中，所述交换步骤仅发生在相邻母线之间。
9.根据权利要求1所述的方法，其中，各个虚拟子图形仅包括一个母线。
10.根据权利要求1所述的方法，其中，各个虚拟子图形包括两个母线。
11.根据权利要求10所述的方法，其中，所述两个母线是发电机母线和负载母线。
12.根据权利要求10所述的方法，其中，所述两个母线是发电机母线。
13.根据权利要求10所述的方法，其中，所述两个母线是负载母线。
14.根据权利要求13所述的方法，其中，所述负载母线中的一个在所述子图形中被表示成真实负载母线和虚拟发电机母线。
CN201510570334.7A 2014-09-10 2015-09-09 利用基于一致性分布式处理估计电网中的最优潮流的方法 Active CN105406459B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US14/482,003 2014-09-10
US14/482,003 US20160071013A1 (en) 2014-09-10 2014-09-10 Method for Estimating Optimal Power Flows in Power Grids using Consensus-Based Distributed Processing

Publications (2)

Publication Number Publication Date
CN105406459A CN105406459A (zh) 2016-03-16
CN105406459B true CN105406459B (zh) 2018-03-16

Family

ID=55437806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510570334.7A Active CN105406459B (zh) 2014-09-10 2015-09-09 利用基于一致性分布式处理估计电网中的最优潮流的方法

Country Status (3)

Country Link
US (1) US20160071013A1 (zh)
JP (1) JP6584214B2 (zh)
CN (1) CN105406459B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9853690B1 (en) 2016-06-13 2017-12-26 International Business Machines Corporation Generating high resolution inferences in electrical networks
WO2018209476A1 (zh) * 2017-05-15 2018-11-22 深圳大学 获取直流电力网功率传输系数的等量电导补偿型偏心方法
CN107508284B (zh) * 2017-08-15 2020-05-19 华北电力大学 计及电气互联的微电网分布式优化调度方法
CN107947152B (zh) * 2017-11-02 2020-04-17 广西电网有限责任公司电力科学研究院 一种基于虚拟电网的中长期动态仿真安稳策略建模方法
JP7310191B2 (ja) * 2019-03-19 2023-07-19 富士電機株式会社 運転計画方法および運転計画装置
US11689618B2 (en) * 2019-12-18 2023-06-27 Hitachi Energy Switzerland Ag Data exchange and processing synchronization in distributed systems
CN111555370B (zh) * 2020-05-20 2023-08-11 云南电网有限责任公司电力科学研究院 一种基于云边协同的配电网分层协调调度方法及装置
CN114818784B (zh) * 2022-04-01 2023-08-04 武汉工程大学 一种改进的联合协方差矩阵和admm算法的稳健波束形成方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101141064A (zh) * 2007-09-14 2008-03-12 清华大学 通过交换边界节点状态和网损信息完成分布式潮流分析的方法
CN101795018A (zh) * 2009-12-31 2010-08-04 华北电力大学 基于可视化的电网智能调度技术支持系统
CN102033993A (zh) * 2010-12-07 2011-04-27 中国电力科学研究院 大规模电力系统动态仿真中继电保护模型的构建方法
CN103530736A (zh) * 2013-10-23 2014-01-22 华北电力大学 基于web-svg的分布式潮流计算系统图形平台的构建方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3399778B2 (ja) * 1997-04-25 2003-04-21 株式会社日立製作所 電力系統評価装置及び電力系統潮流最適化方法と装置並びに電力系統計画支援方法と装置
JP3677668B2 (ja) * 1997-12-12 2005-08-03 株式会社日立製作所 電力系統安定化システム及び電力系統安定度監視システム
US8655499B2 (en) * 2010-02-19 2014-02-18 The Boeing Company Controlling virtual power circuits
US8626353B2 (en) * 2011-01-14 2014-01-07 International Business Machines Corporation Integration of demand response and renewable resources for power generation management
US20130253898A1 (en) * 2012-03-23 2013-09-26 Power Analytics Corporation Systems and methods for model-driven demand response
WO2013187975A1 (en) * 2012-06-15 2013-12-19 Abb Research Ltd. Parallel computation of dynamic state estimation for power system
US9093842B2 (en) * 2012-08-16 2015-07-28 Mitsubishi Electric Research Laboratories, Inc. Method for globally optimizing power flows in electric networks
US9184589B2 (en) * 2013-02-27 2015-11-10 Mitsubishi Electric Research Laboratories, Inc. Method for optimizing power flows in electric power networks
WO2015028840A1 (en) * 2013-08-26 2015-03-05 Ecole Polytechnique Federale De Lausanne (Epfl) Composable method for explicit power flow control in electrical grids
US9419437B2 (en) * 2013-12-19 2016-08-16 Mitsubishi Electric Research Laboratories, Inc. Finite time power control for smart-grid distributed system
US9785130B2 (en) * 2014-04-10 2017-10-10 Nec Corporation Decentralized energy management platform
US9954362B2 (en) * 2014-05-23 2018-04-24 California Institute Of Technology Systems and methods for optimal power flow on a radial network

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101141064A (zh) * 2007-09-14 2008-03-12 清华大学 通过交换边界节点状态和网损信息完成分布式潮流分析的方法
CN101795018A (zh) * 2009-12-31 2010-08-04 华北电力大学 基于可视化的电网智能调度技术支持系统
CN102033993A (zh) * 2010-12-07 2011-04-27 中国电力科学研究院 大规模电力系统动态仿真中继电保护模型的构建方法
CN103530736A (zh) * 2013-10-23 2014-01-22 华北电力大学 基于web-svg的分布式潮流计算系统图形平台的构建方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Distributed Consensus-Based Economic Dispatch With Transmission Losses;Giulio Binetti,et al.;《IEEE transactions on power systems》;20140731;第29卷(第4期);第1711-1720页 *

Also Published As

Publication number Publication date
CN105406459A (zh) 2016-03-16
US20160071013A1 (en) 2016-03-10
JP6584214B2 (ja) 2019-10-02
JP2016059261A (ja) 2016-04-21

Similar Documents

Publication Publication Date Title
CN105406459B (zh) 利用基于一致性分布式处理估计电网中的最优潮流的方法
Wang et al. Peer-to-peer control for networked microgrids: Multi-layer and multi-agent architecture design
Liu et al. Energy management of cooperative microgrids: A distributed optimization approach
Wu et al. Distributed optimal coordination for distributed energy resources in power systems
Chen et al. Distributed optimal active power control of multiple generation systems
Guo et al. Distributed economic dispatch for smart grids with random wind power
Nair et al. Multi-agent systems for resource allocation and scheduling in a smart grid
Li et al. Decentralized multi-area dynamic economic dispatch using modified generalized benders decomposition
Xu et al. Novel multiagent based load restoration algorithm for microgrids
Min et al. Total transfer capability computation for multi-area power systems
Duan et al. Robust consensus-based distributed energy management for microgrids with packet losses tolerance
Zhao et al. Consensus-based distributed optimal energy management with less communication in a microgrid
Chang et al. A distributed robust optimization approach for the economic dispatch of flexible resources
Liu et al. Consensus-based distributed optimal power flow algorithm
Chen et al. Distributed finite-step iterative algorithm for economic dispatch of generation
Zhou et al. Hierarchical distributed voltage regulation in networked autonomous grids
Mohammadi et al. Asynchronous distributed approach for DC optimal power flow
Mangiatordi et al. Multi agent system for cooperative energy management in microgrids
Zhang et al. Distributed DCOPF with flexible transmission
Du et al. Efficient forecasting scheme and optimal delivery approach of energy for the energy Internet
Wang et al. A distributed proximal primal–dual algorithm for energy management with transmission losses in smart grid
Cheng et al. Distributed fixed-time secondary control for voltage restoration and economic dispatch of DC microgrids
Hao et al. Distributed piecewise approximation economic dispatch for regional power systems under non-ideal communication
Chen et al. A framework of decentralized electricity market based on the collaborative mechanism of blockchain and edge computing
Srikantha et al. Hierarchical signal processing for tractable power flow management in electric grid networks

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