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

申请号 CN201510570334.7 申请日 2015-09-09 公开(公告)号 CN105406459A 公开(公告)日 2016-03-16
申请人 三菱电机株式会社; 发明人 M·本诺斯曼; 刘杰; A·拉格胡娜汉;
摘要 利用基于一致性分布式处理估计 电网 中的最优潮流的方法。一种估计电网中的最优潮流OPF的方法,所述电网被表示成图形,所述图形被分割为虚拟子图形,各个虚拟子图形包括至少一个 母线 并且与代理关联,所述代理测量局部变量并且更新一致性变量CV。利用代理交换并更新相邻虚拟子图形的一致性变量。利用代理基于CV和局部变量针对虚拟子图形求解OPF问题。 迭代 所述交换和所述求解步骤,直至满足终止条件为止,此 时针 对各个虚拟子图形输出最优OPF。
权利要求

1.一种用于估计电网中的最优潮流OPF的方法,其中,所述电网包括通过母线连接的发电机和负载,所述方法包括以下步骤:
将所述电网表示成图形,该图形被分割为虚拟子图形,其中,各个虚拟子图形包括至少一个发电机、负载和一个母线;
将代理与各个虚拟子图形关联,其中,所述代理包括计算和通信能
由所述代理在各个虚拟子图形中测量局部变量并且获得一致性变量CV,其中,所述局部变量和所述一致性变量包括受到所述虚拟子图形中的功率平衡约束的电压变量和功率变量;
利用所述代理交换和更新相邻虚拟子图形的所述CV;
利用所述代理,基于所测量的局部变量和所交换的一致性变量来求解各个虚拟子图形的OPF问题;
迭代进行所述交换步骤和所述求解步骤,直至满足终止条件为止;以及输出各个虚拟子图形的最优OPF。
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.根据权利要求7所述的方法,其中,所述惩罚函数利用1范数来构造。
11.根据权利要求1所述的方法,其中,各个虚拟子图形仅包括一个母线。
12.根据权利要求1所述的方法,其中,各个虚拟子图形包括两个母线。
13.根据权利要求12所述的方法,其中,所述两个母线是发电机母线和负载母线。
14.根据权利要求12所述的方法,其中,所述两个母线是发电机母线。
15.根据权利要求12所述的方法,其中,所述两个母线是负载母线。
16.根据权利要求15所述的方法,其中,所述负载母线中的一个在所述子图形中被表示成真实负载母线和虚拟发电机母线。

说明书全文

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

技术领域

[0001] 本发明总体上涉及电网,更具体地讲,涉及利用分布式处理估计电网中的潮流。

背景技术

[0002] 电网经由传输线和变电站从发电机向消费者(负载)分配电力。最优潮流(OPF)是在电网运行中要求解的关键问题。在发电限制、电压限制和传输线路热限制的约束下,OPF问题的解估计有功和无功发电以及电网中的母线的电压以使发电成本最小化。OPF对电网的可靠运行很关键。
[0003] OPF难以求解,因为该问题是非凸和非线性的。通常,由于相邻母线的电压之间的二次关系,非凸性难以确保全局最优解。OPF是在发电、所有母线中的电压以及传输线路热限制方面具有大量决策变量和约束的大规模优化问题。对于实际实现,需要降低求解该问题的计算复杂性。
[0004] 实践中,通过用DC潮流问题(线性规划问题)近似AC潮流问题来避免求解OPF的困难。这种近似具有输电网可接受的精度平,并且被许多独立系统运营商(ISO)使用。然而,DC潮流解对于配电网而言不够精确,无法满足智能电网的要求(包括可再生能源、分布式发电和存储)。
[0005] 另一方法基于半正定规划(SDP)松弛使用AC OPF问题的全局最优解。即使在多个情况下可应对该问题的非凸性,但是集中式方法仍无法满足现代电网的要求。由于连接到电网的不同的、时变的并且易失性的负载和存储,集中式OPF问题的解应当反映电网的任何显著改变。
[0006] 为了满足电网运行的实时要求,通常诸如区域输电运营商(RTO)和ISO的发电商每五分钟求解一次OPF问题。结果,大规模OPF问题需要集中式方法无法提供的效率和精度。
[0007] 因此,大规模OPF问题的可扩展、快速求解使用分布式方法,其中OPF问题被分割为多个小规模子问题。各个子问题通过具有代理至代理通信能力的单个代理作为计算实体来求解。代理根据通信协议交换数据。因此,所有代理合作地以分布式方式参与求解OPF问题。
[0008] 在一个分布式方法中,电网被分割为极大团,参见Lam等人的“Distributed algorithms for optimal power flow problem”(IEEE第51届决策和控制(CDC)年会,第430–437页,2012年)。该方法利用松弛的半正定规划(SDP)求解OPF问题,以凸化(即,近似)该问题。在问题的SDP松弛之后,该方法使用交替方向乘子法(ADMM)来以分布式方式求解凸化问题。凸化方法在收敛时无法确保提供最优潮流问题的可行点。
[0009] 在另一方法中,各个母线被当作单独的代理,其维持相邻母线的电压的估计,并且通过使发电的局部成本和估计误差最小化来构造约化OPF问题,参见Dall’Anese等人的“Distributed optimal power flow for smart microgrids”(IEEE智能电网会刊,第4卷,第3期,第1464–1475页,2013年)。对各个代理的局部等式约束是各个母线的功率平衡方程,其中其它母线的电压变量由估计取代。当方法达到可行解时,所述估计收敛于真值。该方法还使用SDP松弛,然后使用对偶算法来进行分布式优化,其也无法确保OPF问题的可行点。
[0010] 另外,所述方法仅考虑电网作为树建模的情况,并且那两种方法均使用极大团分解。在这两种方法中,用于表示电网的图形(graph)限于树结构,所述图形被分解为极大团。早前的图形和团中的节点和边缘仅表示电网的真实物理组件(例如,实际母线和传输线路)。

发明内容

[0011] 本发明的实施方式提供了一种利用基于一致性(consensus)的分布式处理来估计电网中的最优潮流(OPF)的方法。该方法由代理执行,其中各个代理与电网的至少一条母线关联。各个代理与局部变量和一致性变量关联,所述一致性变量是各个局部变量的估计。
[0012] 局部变量表示:其电压和功率变量;以及其邻近代理的电压和功率变量的估计。局部变量由代理利用局部优化问题来计算,所述局部优化问题使母线处的发电成本以及局部变量相对于其一致性变量的偏差最小化。所述优化问题受到局部功率平衡约束,其表示连接代理与其邻居的线路上的潮流。
[0013] 各个代理向其邻近代理通信:(i)其电压和功率以及(ii)其邻近代理的电压和功率的一致性变量。各个代理应用一致性滤波器以更新其一致性变量。优化和通信的处理被重复,直至各个代理的局部变量和一致性变量收敛于容限内。
[0014] 所述方法将电网表示成图形,所述图形被分割为虚拟子图形。与现有技术不同,子图形的顶点(也被称作节点)无需表示真实物理组件(即,与真实物理发电机关联的母线和真实负载)。也就是说,所述分割步骤可从真实组件生成虚拟组件,即,虚拟发电机或虚拟负载。这有利于求解例如单个发电机连接到多个负载等时的OPF问题(其无法通过其它方法来处理)。
[0015] 所述方法可使用两个不同的实施方式。在一个实施方式中,各个虚拟子图形仅包括与一个代理关联的一个母线。所述代理共享电压并且估计通过传输线路连接的邻近母线的电压。
[0016] 在另一实施方式中,虚拟子图形包括至少两个母线,各个虚拟子图形与一个代理关联。仅针对连接虚拟子图形的母线估计潮流。附图说明
[0017] 图1A是根据本发明的实施方式的电网中的母线和关联的代理的框图
[0018] 图1B是根据本发明的实施方式的利用基于一致性的分布式处理来估计电网中的最优潮流(OPF)的方法的框图;
[0019] 图1C是根据本发明的实施方式的具有三个发电机并且分割为两个虚拟子图形的网络的框图;
[0020] 图2是根据本发明的实施方式的具有两个发电机和一个负载网络并且利用虚拟发电机和虚拟负载分割为两个虚拟子图形的网络的框图;
[0021] 图3是根据本发明的实施方式的具有两个发电机和一个负载网络并且利用两个虚拟负载分割为两个虚拟子图形的网络的框图;
[0022] 图4是根据本发明的实施方式的具有一个发电机和两个负载并且利用虚拟负载和虚拟发电机分割为两个虚拟子图形的网络的框图;
[0023] 图5是根据本发明的实施方式的求解分散式OPF问题的一般情况的框图,其利用虚拟子图形方法将各个虚拟子图形与局部子OPF问题关联;以及
[0024] 图6是根据本发明的实施方式的求解分散式OPF问题的具体情况的框图,其利用直接方法将各个母线与局部子OPF问题关联。

具体实施方式

[0025] 本发明的实施方式提供了一种利用基于一致性的分布式处理用于电网中的最优潮流(OPF)的方法。
[0026] 电网表示
[0027] 如图1B所示,我们将电网建模成图形 50,其中顶点集 包括电网中的所有母线,顶点之间的链路ε是连接母线的高架或地下传输线路。
[0028] 各个母线 可连接到发电机或负载。连接到发电机的母线集合被表示成为了方便,连接到发电机的母线被称为发电机母线。没有任何发电机的母线被
称为负载母线。我们使用i~j来表示母线i连接到母线j。
[0029] 母线i处的有功和无功发电功率被表示成 和 母线i处的功率需求被表示成 母线i处的复数电压被表示成Vi=ei+jfi,其中ei是电压的实部(Re),fi是虚
部(Im), 为了简单,我们将ei和fi统称为Vi。我们将母线i的接地导纳表示成
yii。我们将母线i和j之间的线路导纳表示成yij。导纳是复数yij=gij+jbij。如果则yij=0。
[0030] 最优潮流(OPF)问题
[0031] 所有母线的OPF问题可被表示成
[0032]
[0033] 使得
[0034]
[0035] 以及 (4)
[0036]
[0037] 其中,成本函数Fi可被表示成二次方程
[0038]
[0039] 其中c1,i、c2,i是预定正常数。
[0040] 方程(2)和(3)是母线i的真实无功功率平衡方程。方程(4)中的不等式是电压大小的上下界。方程(5)中的不等式是真实无功发电功率的界。
[0041] OPF的非凸性来自功率平衡方程(2)和(3)中的电压之间的二次关系。这使得难以求解OPF问题。另外,电网中的母线的数量对于以集中方式求解OPF而言常常过大。
[0042] 因此,我们在方程(1)中描述求解OPF问题的分布式方法。我们的方法具有两个变体。在一个实施方式中,各个母线与一个代理关联。在另一实施方式中,电网的图形表示被分割为虚拟子图形,各个虚拟子图形具有一个母线并且与一个代理关联。然而,我们需要指出的是,与现有技术不同,我们仍直接求解非凸问题。另外,代替使用ADMM来求解该问题,我们使用包括惩罚的规划。
[0043] 基于一致性的分布式潮流估计的直接方法
[0044] 如图1A所示,对于一个实施方式,电网中的各个母线10与一个代理15关联。在下述的另一实施方式中,代理与两个或更多母线关联。
[0045] 各个代理包括计算和通信能力。例如,代理可包括连接到存储器21、输入/输出(I/O)接口22和收发器23的处理器20。该处理器可以是多核处理器微处理器、并行处理器等。
[0046] 图1B示出利用基于一致性的分布式处理来估计电网中的最优潮流(OPF)的一般方法的框图。
[0047] 由图形 50表示的电网被分割(55)成虚拟子图形60。各个虚拟子图形至少包括母线10。代理15与虚拟子图形60关联(70)。图形和子图形可被存储在存储器中。
[0048] 代理测量(80)局部变量82并且利用处理器和I/O接口来更新(85)虚拟子图形中的一致性变量。
[0049] 局部变量包括电压和功率变量。如果代理母线的母线以及与所述母线相邻的那些母线,则局部变量的一致性变量也由代理存储。局部变量受到虚拟子图形中的功率平衡约束。利用惩罚函数客观地对局部变量相对于一致性变量的偏差进行惩罚。所述变量可被存储在存储器中。
[0050] 代理利用收发器在虚拟子图形60之间交换(85)局部变量。通过应用一致性滤波器86来更新一致性变量。
[0051] 利用代理基于局部变量和更新的一致性变量求解(90)子图形的OPF问题。交换、更新和求解迭代,直至满足终止条件,输出虚拟子图形的最优OPF 95。
[0052] 为了描述简单,我们不在代理15与母线10之间进行区分,除非特别需要如此做。
[0053] 我们将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的决策变量的部分。
[0054] 局部OPF问题
[0055] 母线i的局部OPF问题可被如下表示。
[0056]
[0057]
[0058]
[0059]
[0060]
[0061] 以及 (11)
[0062]
[0063] 其中方程(12)中的约束确保在母线i处形成的母线j的电压的估计与真值Vj一致。
[0064] 为了求解方程(6)中的局部OPF问题,各个母线i使发电成本最小化,并且估计的ej(i),fj(i)遵循真实电压ej,fj,即,满足方程(12)中的约束。发电成本被包括在方程(6)中的局部OPF问题中。
[0065] 估计的一致性方法
[0066] 对于各对真实局部电压ej(i)、fj(i),母线i维持一对一致性变量 并且通过应用以下一致性滤波器来利用从母线j获得的ej,fj的真实局部值更新这些变量:
[0067]
[0068] 并且
[0069] 其中,0<γ<1是一致性增益。
[0070] 即,在(k+1)迭代更新的一致性变量的值的估计是(k)迭代的先前值的估计与(先前真实值与先前估计值之间)的差乘以一致性增益之和。
[0071] 在我们的方法中,估计(例如,ej(i)和fj(i))未直接被来自其它母线的真值(例如,ej,fj)取代。相反,所述值通过一致性滤波器。一致性变量在该滤波器中可被视为中间变量,并且用于在局部优化中更新估计。
[0072] 母线i和j在各个迭代k交换真实电压及其一致性变量。换言之,母线i接收真实电压ej(k),fj(k)和一致性变量 并且与相邻母线交换电压ei(k),fi(k)和一致性变量 由于交换仅发生在相邻母线之间,所以与需要获取来自每一个母线的
所有数据并且由单个处理中心处理的集中式方法相比,通信开销的额外负担很小。
[0073] 局部优化问题
[0074] 为了使发电成本和估计误差二者最小化,重新规划方程(6)中的局部优化问题。对于母线 局部优化问题被表示为
[0075]
[0076]
[0077]
[0078] 以及 (17)
[0079]
[0080] 其中ρi>0是以较高优先级确保估计误差被最小化的大的正惩罚因子。这是必要的,因为在没有这种惩罚的情况下,各个母线可基于有偏估计集合来自私地使它自己的发电成本最小化。
[0081] 方程(14)中的问题的决策变量仅与母线i和相邻母线有关,因此该问题是具有相对较小规模的局部非线性优化问题。可利用非线性规划求解器来有效地求解该问题,例如内点优化器(IPOPT)或者“寻找约束非线性多元函数的最小值”(fmincon)。
[0082] 我们可将母线i的分布式估计方法总结如下。
[0083] (1)初始化变量:
[0084] ei(0),fi(0),ej(i)(0),fj(i)(0),
[0085] (2)在迭代k+1,如果满足终止条件,则结束。否则,各个代理与相邻母线j~i交换变量ei(k),fi(k),
[0086] (3)如方程(13)中那样更新一致性变量
[0087] (4)通过求解方程(14)中的优化问题来更新决策变量,即,
[0088]
[0089] 然后,进行到步骤(2)。
[0090] 在另一实施方式中,局部优化问题的目标函数中的惩罚可利用1范数来规划成[0091]
[0092] 此方法非常适用于求解所有母线均为发电机母线的情况。然而,对于具有负载母线的电网,该方法可对初始值敏感,并且常常对负载母线具有稳态估计误差。此问题是由于负载母线中缺少自由度导致的。为了解决此问题,我们对虚拟子图形使用一致性滤波。
[0093] 基于虚拟子图形的类一致性分布式优化方法
[0094] 此实施方式求解电网中存在负载母线时的问题。我们将电网分割为虚拟子图形。在一个实施方式中,各个虚拟子图形与两个母线之间的潮流有关。虚拟子图形与代理关联,该代理按照与上述直接方法相似的方式执行分布式方法。
[0095] 虚拟子图形
[0096] 图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进行连接。
[0097] 存在三种类型的虚拟子图形:
[0098] 情况1:一个发电机母线和一个负载母线;
[0099] 情况2:两个发电机母线;以及
[0100] 情况3:两个负载母线。
[0101] (1)一个发电机母线和一个负载母线
[0102]
[0103]
[0104] 以及 (22)
[0105]
[0106] 其中Ni是母线i的相邻母线的数量。需要指出的是,功率项 和电压Vi|(i,j),Vj|(i,j)是由虚拟子图形维持的局部变量。具体地讲,电压Vi|(i,j)是由虚拟子图形(i,j)估计的母线i的电压。该电压仅是估计,因为另一虚拟子图形也可包括母线i并且确定它自己的估计。为了确保相同变量的一致性估计,对变量强加以下条件。
[0107]
[0108] 以及 (25)
[0109]
[0110] 功率的局部变量(例如, )可被解释为在虚拟子图形(i,j)的分支流中所涉及的那部分功率。换言之,我们将真实发电或负载分割为多个部分,并且将各个部分指派给不同的虚拟子图形。因此,如方程(24)中的条件所强制的,那些局部变量(例如,发电功率 )之和必须等于母线i处的发电的真实值。
[0111] (2)两个发电机母线
[0112] 发电机母线也可具有直接附接的负载。因此,灵活的是将发电机母线建模成虚拟子图形中的负载母线或发电机母线。建模可针对各种情况而不同。在具有彼此连接的两个发电机母线的虚拟子图形的情况下,考虑两种情况。
[0113] 在第一种情况中,发电机母线没有足够的功率来支持连接的负载或者相邻虚拟子图形中的负载母线。在这种情况下,该发电机母线可被建模成具有供应所需功率的另一发电机母线(其也被建模成包括所述负载的虚拟子图形中的发电机母线)的虚拟子图形中的负载母线。
[0114] 对于两个发电机母线i1和i2,我们可假设潮流从母线i1至母线i2。因此,我们将母线i2建模成负载母线。流方程可被表示成
[0115]
[0116]
[0117]
[0118]
[0119] 需要指出的是,母线i2的负载 可被视为消耗来自母线i2的潮流的虚拟负载。这意味着 的值可超过真实负载 这发生在电网超过母线j的部分可能
需要通过母线i2供应的母线i1的功率时。为了一致性,母线i2还必须在其它虚拟子图形中作为发电机母线出现。
[0120] 图2示出图形210,其表示包括发电机母线i1201、发电机母线i2202和负载母线j203的电网。我们使用圆圈来表示发电机母线,正方形表示负载母线。图形210被分割(55)成两个虚拟子图形211和212。然后,在虚拟子图形212中,则看到母线i2连接到虚拟负载
221,真实无功功率被定义为 在虚拟子图形212中,则看到母线i2连
接到虚拟发电机222。
[0121] 因此,在虚拟子图形(i2,j)212中,存在与虚拟发电机222关联的虚拟发电功率项使得 并且
[0122] 图2示出具有两个发电机母线的虚拟子图形的示例。在第一种情况下,虚拟子图形(i2,j)的功率平衡方程可被表示成
[0123]
[0124]
[0125]
[0126]
[0127] 图3上所示的第二种情况是两个发电机母线共享母线之一的负载,使得母线i1301和母线i2302向母线j 303处的负载供应功率。在这种情况下,虚拟子图形可通过将母线i2302视为虚拟子图形(i1,i2)311中的负载母线来形成,其中母线i2321处的负载部分被定义为使得 剩余负载可出现在其它虚拟子图形中,其中母线i2被视为发电机母线324。
[0128] (3)两个负载母线
[0129] 对于包含两个负载母线j1402和j2403的虚拟子图形,潮流从母线j1至母线j2。我们假设存在连接到母线j1422的虚拟发电机,其供应功率 我们可将潮流方程写为
[0130]
[0131]
[0132]
[0133]
[0134] 需要指出的是虚拟发电项 作为虚拟负载402出现在连接到母线j1的发电机母线401中。这可维持电网的功率平衡。
[0135] 图4示出连接到两个负载母线j1402和j2403的发电机母线401的图形。该图形被分割为虚拟子图形411和412,其中负载母线j1402现在连接到发电机母线i1,具有真实无功功率 和 的虚拟发电机母线422连接到负载母线i2。因
此,在虚拟子图形412中,虚拟发电功率项 被增加到母线j1使得
并且
[0136] 具体地讲,虚拟子图形(i,j1)的功率平衡方程可被表示成
[0137]
[0138]
[0139]
[0140] 流方程之和等同于功率平衡方程
[0141] 在电网中,对于各个母线i,假定方程(24-26)中对局部变量的条件有效,则包含母线i的所有虚拟子图形(例如,(i,j)和(j,i))中的流方程之和等同于母线i的功率平衡方程。
[0142] 考虑简单的示例。假定所有虚拟子图形包含发电机母线和负载母线。假定所有发电机母线充当虚拟子图形中的发电机母线。对于各个发电机母线i,我们可在所有相邻母线j上对发电机母线的流方程求和。
[0143]
[0144]
[0145] 通过方程(24-26)中的条件,我们知道与母线i有关的所有虚拟子图形上的发电功率 之和等于母线i的原始发电功率,所有其它局部变量也是如此。因此,我们可恢复功率平衡方程如下。
[0146] 以及 (45)
[0147]
[0148] 对于虚拟子图形中的负载母线j,可通过对出现在所有虚拟子图形(i,j)中的流方程求和,使得i~j来恢复功率平衡方程。方程之和被表示为
[0149]
[0150]
[0151] 对于具有两个负载母线的虚拟子图形,这两个负载母线中的一个具有额外虚拟发电功率,此额外功率被增加到具有真实发电机母线的虚拟子图形的项抵消。考虑(j1,j2)均为负载母线的虚拟子图形中的母线j1。母线j1也在虚拟子图形(i,j1)中,其中母线i是发电机母线。然后,虚拟子图形(j1,j2)中的母线j1的流方程被表示成
[0152]
[0153]
[0154] 因此,在虚拟子图形(i,j1)中,母线j1的流方程为
[0155]
[0156]
[0157] 由于我们强制 所以当关于母线j1的所有流方程被加起来时,项被抵消。这确保了流方程之和等同于负载的功率平衡方程。
[0158] 虚拟子图形方法中的局部估计的一致性方法
[0159] 在分布式优化中各个虚拟子图形是独立的代理。如上所述利用一致性滤波器来执行一致性变量的交换和更新。需要指出的是,在此实施方式中,各个虚拟子图形(i,j)具有母线i和j的电压值。然而,由于母线i可与不止一个虚拟子图形关联,所以电压vi|(i,j)的值必须与其它母线处的电压相容。也对发电功率和负载功率的值强加相同的一致性要求。那些关系由方程(24-26)来强加,其对获得全局问题的可行解而言是必要的。
[0160] 为了满足这些要求,各个虚拟子图形(i,j)维持母线i的电压、发电和其它虚拟子图形(i,j′),j′≠j,j′~i的负载的估计。因此,由虚拟子图形(i,j)估计的虚拟子图形(i,j′)中的i的发电功率被表示为 负载估计为
[0161] 对于与虚拟子图形(i,j)共享母线i的虚拟子图形,各个估计的一致性变量为根据以下规则来更新一致性变量
[0162]
[0163]
[0164]
[0165]
[0166] 以及
[0167]
[0168] 其中常数0<γ<1是一致性增益,其可被适当地选择以确保一致性算法的收敛。一致性处理与其它虚拟子图形通信以交换一致性变量。
[0169] 局部优化问题
[0170] 类似于直接方法,针对各个虚拟子图形规划局部优化问题以使发电成本和估计误差最小化。在不损失一般性的情况下,我们考虑虚拟子图形(i,j),其具有还出现在其它虚拟子图形中的母线i和j。局部优化问题可被表示成以下目标函数:
[0171]
[0172]
[0173]
[0174]
[0175]
[0176]
[0177]
[0178]
[0179]
[0180]
[0181] 各个虚拟子图形内的局部优化问题针对两个母线。各个虚拟子图形的相邻虚拟子图形的数量等于虚拟子图形内的相邻母线i和j的数量。因此,这是小规模OPF问题,可通过类似fmincon和IPOPT的可用的非线性优化求解器来有效地求解。
[0182] (1)初始化fk|(i,j)(0)(k=i,j)以及估计变量
(k=i,j)和一致性变量的值。
[0183] (2)在迭代k+1,如果达到终止条件(例如,最大迭代次数或者估计的误差容限),则方法终止。如果没有,则各个代理将局部变量发送给相邻虚拟子图形。
[0184] (3)根据(53)更新一致性变量
[0185] (4)通过求解优化问题(54)来更新决策变量,然后进行到步骤(2)。
[0186] 图5是根据本发明的实施方式的利用直接方法求解将各个母线与局部OPF子问题关联的分散式OPF问题的一般情况的框图。OPF问题600被分割(605)成子OPF问题610,每个虚拟子图形1、2、...、n-1、n一个。然后,对于各个虚拟子图形的各个子问题610,以下步骤迭代,直至达到终止条件。在各个虚拟子图形处应用一致性滤波器620以更新一致性变量 和 针对各个虚拟子图形求解(630)子OPF问题以更新决策变量然后,在终止处,输出(640)各个虚拟子图形的最优电压和功率。
[0187] 图6是根据本发明的实施方式的利用直接方法求解将各个母线与局部OPF子问题关联的分散式OPF问题的具体情况的框图。OPF问题700被分割(705)为子OPF问题710,每个母线1、2、...、n-1、n一个。然后,对于各个子问题,以下步骤迭代,直至达到终止条件。应用一致性滤波器720以更新一致性变量 和 求解(730)子OPF问题以更新局部
决策变量 然后,在终止处,输出(740)各个母线的最优电压和功率。
[0188] 本发明的效果
[0189] 实施方式提供了一种基于一致性的分布式优化方法。该方法可用于求解最优潮流问题。该方法以分布式方式在彼此通信的多个代理中执行。描述两个实现方式。在直接方法中,各个母线与代理功率。各个代理在关联的母线的功率平衡方程的约束下求解局部OPF问题。代理估计相邻代理的电压并且与其它代理交换一致性变量以达成一致性。使用一致性滤波方法来确保估计的收敛。在虚拟子图形方法中,各个虚拟子图形与相邻母线之间的分支流有关。此方法在多个负载连接到发电机母线的情况下更有效。针对所有变量达成一致性。
[0190] 尽管通过优选实施方式的示例描述了本发明,将理解,在本发明的精神和范围内可进行各种其它更改和修改。因此,所附权利要求书的目的在于涵盖落入本发明的真实精神和范围内的所有这些变化和修改。
QQ群二维码
意见反馈