郭峰，趙逸嵐，楊歡

（1．國網(wǎng)陜西省電力公司電力科學(xué)研究院，陜西西安 710100；2．杭州市電力設(shè)計院有限公司，浙江 杭州 310062；3．浙江大學(xué)，浙江杭州 310027）

通信延時環(huán)境下虛擬發(fā)電廠分布式控制收斂速度分析

郭峰1，趙逸嵐2，3，楊歡3

（1．國網(wǎng)陜西省電力公司電力科學(xué)研究院，陜西西安 710100；2．杭州市電力設(shè)計院有限公司，浙江 杭州 310062；3．浙江大學(xué)，浙江杭州 310027）

針對一種考慮通信延時環(huán)境下的應(yīng)用于虛擬發(fā)電廠的分布式一致性控制策略進(jìn)行研究，分析了控制策略收斂速度的影響因素，并據(jù)此提出了收斂速度衡量指標(biāo)。在仿真算例中，分析了線路權(quán)重、領(lǐng)導(dǎo)節(jié)點位置以及算法收斂系數(shù)大小對收斂速度的影響，并對比了有/無延時兩種情況下一致性算法的收斂速度差異，為通信網(wǎng)的優(yōu)化設(shè)計提供了參考。

分布式控制；虛擬發(fā)電廠；通信延時；收斂速度

隨著當(dāng)今全球能源形勢的逐漸緊張以及環(huán)境污染問題的日益凸顯，可再生清潔能源的大力開發(fā)利用已成為電力行業(yè)的革新趨勢，這也促進(jìn)了分布式電源的使用和發(fā)展。分布式電源是指位于負(fù)荷側(cè)，就近接入中低壓配電網(wǎng)的小型發(fā)電系統(tǒng)[1]，主要包含光伏發(fā)電、風(fēng)力發(fā)電、生物質(zhì)能發(fā)電等可再生能源發(fā)電單元，以及燃?xì)廨啓C等清潔能源發(fā)電單元[2]。

虛擬發(fā)電廠（virtual power plant，VPP）是一種電源管理技術(shù)，該技術(shù)將多個分散小電源聚合起來，使其整體的出力外特性與傳統(tǒng)發(fā)電機類似。應(yīng)用虛擬發(fā)電廠技術(shù)，能有效地對地理位置分散，單機容量小的分布式電源進(jìn)行控制管理。在虛擬發(fā)電廠的控制模式選取上，傳統(tǒng)的集中式控制模式在管理規(guī)模較大且地理分布較為廣泛的受控電源時，存在通信帶寬較窄、控制中心計算負(fù)擔(dān)較重等問題，并且某一信道的通信故障有可能會對整個系統(tǒng)造成嚴(yán)重的影響[3]。相比之下，分布式控制模式則無需通過一個控制中心來對所有電源下發(fā)控制指令；各電源在與少數(shù)鄰居電源通信的基礎(chǔ)上即可計算得出各自的控制參數(shù)值。因此，分布式控制模式具有經(jīng)濟性好、魯棒性強以及通信容錯程度高等特點，更適用于對虛擬發(fā)電廠內(nèi)部分布式電源的管理。

一致性算法是一種廣泛應(yīng)用于多智能體協(xié)同控制的分布式算法。多智體協(xié)同控制的本質(zhì)在于，通過多個價格低廉且結(jié)構(gòu)簡單的智能體的聯(lián)合工作，代替單個復(fù)雜且昂貴的個體[4]。目前，一致性算法在控制、計算機領(lǐng)域已有大量的研究成果，可應(yīng)用于解決智能體編隊控制[5]、無線傳感網(wǎng)絡(luò)[6]、多航天器控制[4]等問題。除此之外，一致性算法也被引入到了電力系統(tǒng)領(lǐng)域，用以設(shè)計分布式控制策略以對發(fā)電設(shè)備輸出的有功、無功功率進(jìn)行控制[7-8]。文獻(xiàn)[9]提出了一種用于光伏系統(tǒng)有功、無功功率自適應(yīng)控制的分布式策略，實現(xiàn)了系統(tǒng)有功功率的自平衡，并維持了電壓穩(wěn)定。文獻(xiàn)[10]結(jié)合多智體一致性協(xié)同控制理論，設(shè)計了一種分散自治的自動發(fā)電控制功率分配框架。文獻(xiàn)[11]基于一致性算法，提出了一種適用于直流孤島微網(wǎng)的控制策略，實現(xiàn)了變虛擬電阻的自適應(yīng)下垂控制。

分布式控制需要通過多次迭代計算才能實現(xiàn)控制目標(biāo)，因此收斂速度成為了評判控制效果優(yōu)劣的一個重要標(biāo)準(zhǔn)。已有研究表明，相較于無通信延時的情況，一致性算法在考慮通信延時下的收斂性較差[12]。在此基礎(chǔ)上，文獻(xiàn)[12]結(jié)合代數(shù)圖論、矩陣?yán)碚?、控制理論分析了連續(xù)一致性算法的收斂性與延時長短的關(guān)系；文獻(xiàn)[13]將一致性算法應(yīng)用于兩階離散系統(tǒng)，考慮不均勻通信延時以及時變通信結(jié)構(gòu)，推導(dǎo)出了系統(tǒng)實現(xiàn)一致性收斂所需滿足的充分條件。

本文將一階離散一致性算法應(yīng)用于虛擬發(fā)電廠內(nèi)部分布式電源的有功功率控制之中，分析了在考慮均勻通信延時的情況下，控制方法收斂速度的影響因素，并提出了收斂速度衡量指標(biāo)。在仿真算例中，分析了線路權(quán)重、領(lǐng)導(dǎo)節(jié)點位置以及算法收斂系數(shù)大小對收斂速度的影響，并對比了有/無通信延時兩種情況下一致性算法的收斂速度差異，為通信網(wǎng)優(yōu)化設(shè)計提供了參考。

1 虛擬發(fā)電廠分布式控制策略

1.1 虛擬發(fā)電廠分布式經(jīng)濟調(diào)度

虛擬發(fā)電廠中分布式電源的調(diào)度控制結(jié)構(gòu)為：能量管理系統(tǒng)（energy management system，EMS）向虛擬發(fā)電廠下發(fā)總功率指令，在滿足這一總出力目標(biāo)的前提下，各分布式電源的功率分配將由分布式控制方法實現(xiàn)，即分布式電源不直接接受上級的調(diào)度。如圖1所示，根據(jù)地理位置的不同將配電網(wǎng)中的部分分布式電源劃分為一個虛擬發(fā)電廠。虛擬發(fā)電廠中設(shè)有一個領(lǐng)導(dǎo)節(jié)點，用以直接接收上級EMS下發(fā)的功率調(diào)度指令，并將功率偏差信息通過通信網(wǎng)傳遞給所在虛擬發(fā)電廠中的鄰居電源；其余分布式電源作為跟隨節(jié)點，相互之間可以進(jìn)行通信，并根據(jù)鄰域信息決定發(fā)電量。需要注意的是，為確?？刂频挠行?，虛擬發(fā)電廠內(nèi)部的通信網(wǎng)絡(luò)須保持連通[14]。

圖1 虛擬發(fā)電廠分布式控制結(jié)構(gòu)Fig.1 Distributed control architecture of VPP

在該結(jié)構(gòu)中，控制目標(biāo)選取為整體運行經(jīng)濟性最優(yōu)，據(jù)此在分布式電源間進(jìn)行有功出力的分配：

式中：n為該虛擬發(fā)電廠中分布式電源的總數(shù)；Pi為第i個分布式電源的有功出力；Ci（Pi）為第i個分布式電源的發(fā)電成本函數(shù)；Pref為由EMS下發(fā)的整個虛擬發(fā)電廠有功總出力的參考值。

分布式電源的成本函數(shù)可用如下的二次函數(shù)近似表示[15]：

式中，αi、βi、γi為發(fā)電機i的發(fā)電成本系數(shù)。

為實現(xiàn)經(jīng)濟性最優(yōu)目標(biāo)，各分布式電源需滿足等微增率原則[16]，并結(jié)合公式可以得到：

式中，xi為DG i的成本微增率。

在虛擬發(fā)電廠的經(jīng)濟調(diào)度中，分布式控制需要實現(xiàn)以下目標(biāo)：

1）僅使用自己和相鄰單元的信息，無集中控制信號；

2）最終每個分布式電源的成本微增率xi及有功出力Pi需滿足式（3）。

1.2 考慮均勻通信延時的分布式控制算法

假設(shè)分布式電源之間的通信網(wǎng)絡(luò)為雙向通信，對應(yīng)的無向圖G?（V，S）。其中，V={1，…，n}為節(jié)點（分布式電源）集合，S?V×V為邊（通信線路）集合。當(dāng)Vi、Vj間無通信時，邊Sij對應(yīng)的權(quán)重eij=0；當(dāng)Vi、Vj間有通信時，邊Sij對應(yīng)的權(quán)重eij＞0。

本文采用離散一致性算法。為便于描述，暫不考慮來自EMS的總功率指令，僅先實現(xiàn)各分布式電源成本微增率xi相等。

設(shè)迭代次數(shù)為k，信號采樣周期為T，信號由Vi傳輸?shù)絍j的通信延時為τij。在（k+1）T時刻，DGi的微增率xi由下式?jīng)Q定：

當(dāng)采樣周期T大于通信網(wǎng)最大延時τmax時，式（5）等價于由式（6）所表示的無延時情況下的一致性算法：

盡管增長采樣周期可以有效抑制通信延時對控制過程帶來的不利影響，但同時也會導(dǎo)致算法收斂時間的延長。為選取合適的采樣周期，現(xiàn)針對網(wǎng)絡(luò)通信延時較均勻的情況開展研究，以使φ（τij）=1。

公式（6）可簡化為

設(shè)矩陣D為通信網(wǎng)的權(quán)重矩陣，D的對角元素為零，非對角元素Dij等于邊權(quán)eij。有eij=eji，故D為對稱矩陣。將式（7）表示為如下的矩陣形式（x∈R1×1，A∈Rn×n，D∈Rn×n）：

下一步，在公式（8）的基礎(chǔ)上引入總功率約束。即把總功率偏差量信息反饋給領(lǐng)導(dǎo)節(jié)點，構(gòu)成主-從一致性算法，作為整個虛擬發(fā)電廠有功控制策略的核心：

式中：ε為收斂系數(shù)常數(shù)；列向量dn×1的每行對應(yīng)一個分布式電源，領(lǐng)導(dǎo)節(jié)點所在行為1，跟隨節(jié)點所在行為0；矩陣a、b中的元素可由各臺發(fā)電機成本函數(shù)系數(shù)確定。當(dāng)對于任意兩個分布式電源DGi、DGj均滿足時，系統(tǒng)微增率狀態(tài)值完成收斂，從而在完成虛擬電廠內(nèi)部的有功出力分配的同時實現(xiàn)運行經(jīng)濟性最優(yōu)的目標(biāo)。

式（9）所示的一致性算法由兩部分組成：

1）均攤部分：算法x（k+1）=Ax（k）+x（k-1）能通過多次迭代，削減向量x每行元素之間的差別，最終使得每個單元的微增率相等。

2 收斂速度分析

2.1 均攤算法

首先暫不計及領(lǐng)導(dǎo)節(jié)點接收到的功率偏差反饋信號，僅考慮算法的均攤部分。在通信延時均勻、延時時長為τ時，令采樣周期T的取值略大于τ/2，使得對于所有通信線路Sij均存在φ（τij）=1。

現(xiàn)將公式（8）改寫成如下形式：

式中：I為n×n的單位矩陣；W為非負(fù)矩陣，滿足W·1n×1=1n×1，因此W為行隨機矩陣，譜半徑為1[3]。

定理1[3]：定義矩陣A=[aij]∈Rn×n對應(yīng)的拓?fù)鋱D為Γ（A），則當(dāng)且僅當(dāng)Γ（A）為強連通時，矩陣A不可約。

定理2[3]：若矩陣A=[Aij]∈Rn×n為非負(fù)矩陣且不可約，則其譜半徑ρ（A）為其代數(shù)重度為1的特征值。

定理3[3]：設(shè)矩陣A=[aij]∈Rn×n為行隨機矩陣，若A有代數(shù)重度為1的特征值λ=1，且其他特征值均滿足|λ|＜1，則A為SIA（Indecomposable and Aperiodic），有向量v滿足vTA=vT且

由1.2節(jié)可知，矩陣A+D對應(yīng)的拓?fù)鋱D為G，可以證明，對于連通的無向圖G，矩陣W對應(yīng)的有向圖為強連通；由定理1和定理2可知，矩陣W不可約，且有代數(shù)重度為1的特征值λ=1，其他特征值均滿足|λ|＜1。由定理3可得，W為SIA，故當(dāng)k→∞時，列向量中的所有元素會收斂至相等，即能夠?qū)崿F(xiàn)使所有分布式電源的成本微增率相同。

因為W僅有一個模值為1的特征值，而其余特征值的模值均小于1，故收斂過程的收斂速度由矩陣W中模值第二大的特征值決定，記該特征值的模值為μ（W）。本文將μ（W）作為均攤部分收斂速度衡量指標(biāo)，即μ（W）越小，公式式收斂速度越快。

由上述分析知，矩陣A、D為非負(fù)且通信拓?fù)浔３诌B通，是公式成功收斂的充分條件。

2.2 主-從一致性算法

當(dāng)計及接收功率反饋信息的領(lǐng)導(dǎo)節(jié)點后，在通信延時均勻，延時時長為τ的情況下，令采樣周期T的取值略大于τ/2，使得對于所有通信線路Sij均存在φ（τij）。主-從一致性算法（9）可進(jìn)一步改寫為式（11）所示的矩陣形式：

定義誤差ξ（k+1）為：

將式（12）與式（13）相減，有：

當(dāng)且僅當(dāng)M為收斂矩陣，即譜半徑ρ（M）＜1時，ξ（k）會在k→∞的過程中逐漸消減為0[17]。ρ（M）能夠反映Mk的收斂速度，故可作為整個主-從一致性算法收斂速度的衡量指標(biāo)，即ρ（M）越小，算法的收斂速度越快。

通過觀察矩陣M可以發(fā)現(xiàn)，ρ（M）會受到矩陣W、向量dn×1、系數(shù)ε以及對角陣a的影響。這些因素各自的含義為：通信網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)及通信線路權(quán)重、領(lǐng)導(dǎo)節(jié)點位置、功率反饋收斂系數(shù)ε以及各分布式電源發(fā)電成本函數(shù)。因此，在設(shè)計分布式控制的通信網(wǎng)絡(luò)時，可以ρ（M）最小化為目標(biāo)，對上述因素進(jìn)行優(yōu)化。

3 仿真算例

虛擬發(fā)電廠中的分布式電源群組均采用如下形式的發(fā)電成本函數(shù)[7]：

設(shè)虛擬發(fā)電廠中分布式電源數(shù)量為n=15，通信線路數(shù)目為nl=22。能量管理系統(tǒng)下發(fā)的有功出力參考值為Pref=6 MW。

將分布式電源分為5組，每個群組對應(yīng)的發(fā)電成本函數(shù)系數(shù)如表1所示。

表1 分布式電源發(fā)電成本函數(shù)系數(shù)Tab.1 Cost function coefficients of DGs

各分布式電源之間的通信網(wǎng)拓?fù)浣Y(jié)構(gòu)如圖2所示，記該通信網(wǎng)絡(luò)對應(yīng)的拓?fù)鋱D為圖G（圖中所示線路權(quán)重為代數(shù)連通度最大權(quán)重D3，詳見后文）。

采取如下3種不同的權(quán)重選擇方式[18]，應(yīng)用策略式（9）分別進(jìn)行仿真實驗，并記錄相應(yīng)的收斂所需的迭代次數(shù)：

1）最大度權(quán)重D1

圖2 通信網(wǎng)拓?fù)浣Y(jié)構(gòu)Fig.2 Topology structure of a communication network

其中，di為節(jié)點Vi的度，即與節(jié)點Vi直接相連的節(jié)點的數(shù)量。定義此時的通信網(wǎng)為G-D1，在此通信網(wǎng)下，以分布式電源DGk為領(lǐng)導(dǎo)節(jié)點時的矩陣M的譜半徑為ρ（M1-DGk）。

2）均勻權(quán)重D2

算例中，取α=5，定義此時的通信網(wǎng)為G-D2，在此通信網(wǎng)下，以分布式電源DGk為領(lǐng)導(dǎo)節(jié)點時的矩陣M的譜半徑為ρ（M2-DGk）。

3）代數(shù)連通度最大權(quán)重D3

對于通信拓?fù)浯_定的虛擬發(fā)電廠，令l為圖G對應(yīng)的權(quán)重矩陣的集合，則有當(dāng)Sij?G且i≠j時，當(dāng)Sij∈G時

需求解如下問題：

依照本算例的優(yōu)化方法計算出相應(yīng)的權(quán)重結(jié)果如圖2所示，定義此時的通信網(wǎng)為G-D3，在此通信網(wǎng)下，以分布式電源DG k為領(lǐng)導(dǎo)節(jié)點時的矩陣M的譜半徑為ρ（M3-DG k）。

4 仿真結(jié)果分析

4.1. 線路權(quán)重及領(lǐng)導(dǎo)節(jié)點選取對收斂速度的影響

對于本節(jié)的仿真算例，設(shè)收斂變量初始值為x（0）=[1 1 … 1]T。當(dāng)同時滿足條件：

1）總功率差額ΔP＜0.001 MW。

2）對于任意一個分布式電源均有|xi（k+1）-xi（k）|＜0.001，迭代過程結(jié)束。

在通信網(wǎng)G-D1、G-D2和G-D3中，依次選擇1～n號分布式電源作為領(lǐng)導(dǎo)節(jié)點，取收斂系數(shù)ε=0.1，分別計算矩陣M對應(yīng)的譜半徑ρ（M），結(jié)果如圖3所示。按照同樣的步驟，應(yīng)用公式（9）所示策略進(jìn)行仿真，并分別記錄收斂所需的迭代次數(shù)，如圖4所示。需要注意的是，當(dāng)采用通信網(wǎng)G-D3進(jìn)行仿真時，無論領(lǐng)導(dǎo)節(jié)點如何選擇，系統(tǒng)均無法收斂，故未展示于圖中。

圖3 3種通信網(wǎng)的M矩陣譜半徑Fig.3 Spectral radius of matrix M in three kinds of communication networks

由圖3可知，對于通信網(wǎng)G-D3，無論領(lǐng)導(dǎo)節(jié)點如何選擇，均有ρ（M）＞1成立。前文2.2節(jié)的分析表明當(dāng)ρ（M）＞1時，一致性算法無法收斂，仿真結(jié)果也印證了這一結(jié)論。通信網(wǎng)G-D3的線路權(quán)重能夠使得圖G代數(shù)連通度達(dá)到最大，即能使通信拓?fù)浯_定的無通信延時的一致性算法的收斂達(dá)到最快。但需要注意的是，求解過程中所采用的公式（19）所示的半定規(guī)劃方法的線路權(quán)重計算結(jié)果，不一定能夠滿足為非負(fù)矩陣的要求，即在含有通信延時的情況下，求解出不考慮通信延時情況下的最優(yōu)權(quán)重，可能導(dǎo)致算法的不收斂。在本仿真算例中，矩陣A=diag（1-中的元素A2，2、A6，6、A7，7、A8，8、A11，11和A13，13均為負(fù)，因此矩陣W不是非負(fù)的行隨機矩陣，且譜半徑為1.017 5，則在所示算例條件下無領(lǐng)導(dǎo)一致性算法式（8）不收斂，即在通信網(wǎng)G-D3的權(quán)重設(shè)置下，控制策略無法完成成本微增率的均攤；當(dāng)引入能量管理系統(tǒng)的總功率指令后，更無法使系統(tǒng)收斂至穩(wěn)定點。

圖4 通信網(wǎng)G-D1、G-D2對應(yīng)的收斂所需迭代次數(shù)Fig.4 Convergence iterations in communicationnetworks of G-D1and G-D2

此外，對比通信網(wǎng)G-D1及G-D2中矩陣M對應(yīng)的譜半徑曲線可以發(fā)現(xiàn)，通信網(wǎng)G-D2的ρ（M）普遍比通信網(wǎng)G-D1的小，僅在選擇12號節(jié)點為領(lǐng)導(dǎo)節(jié)點時，二者相近。同時，圖4的仿真結(jié)果表明，當(dāng)采用通信網(wǎng)G-D2時，該情況下控制策略式（9）的收斂速度普遍比采用通信網(wǎng)G-D1時快，且僅在選擇12號節(jié)點為領(lǐng)導(dǎo)節(jié)點時兩組通信網(wǎng)迭代次數(shù)結(jié)果相差不大，這一結(jié)論與譜半徑計算結(jié)果一致。此外，由圖3及圖4可以發(fā)現(xiàn)，一致性算法的迭代次數(shù)曲線與相應(yīng)通信網(wǎng)的M矩陣譜半徑曲線走勢相吻合，這一現(xiàn)象也進(jìn)一步說明了采用ρ（M）作為一致性算法收斂速度衡量指標(biāo)的合理性。

因此，對于拓?fù)浣Y(jié)構(gòu)相同的通信網(wǎng)，通信線路權(quán)重以及領(lǐng)導(dǎo)節(jié)點位置的選擇均會對算法的收斂速度產(chǎn)生影響。在優(yōu)化通信網(wǎng)線路權(quán)重，選擇領(lǐng)導(dǎo)節(jié)點位置時，可以以指標(biāo)ρ（M）最小化為目標(biāo)進(jìn)行。

4.2 收斂系數(shù)對收斂速度影響

現(xiàn)選取通信網(wǎng)G-D2為研究對象。算例中，收斂系數(shù)ε在[0.002，1]之間均勻地取500個值；分別選取通信網(wǎng)G-D2中能夠使得指標(biāo)ρ（M）為最小和最大的分布式電源作為領(lǐng)導(dǎo)節(jié)點（分別為DG1和DG14），并記錄收斂所需的迭代次數(shù)（見圖5、圖6）。其中，灰色背景區(qū)域所對應(yīng)的部分表示在相應(yīng)的取值下，算法無法收斂。另外，為凸顯曲線變化較大的區(qū)域，將圖中曲線變化平緩的區(qū)域進(jìn)行折疊（在圖5中折疊0.2＜ε＜0.8區(qū)域部分，在圖6中折疊0.2＜ε＜0.4區(qū)域部分）。

圖5 領(lǐng)導(dǎo)節(jié)點為DG1：迭代次數(shù)、與收斂系數(shù)的關(guān)系曲線Fig.5 Graph of the relationships between iterations/and with leader DG1

圖6 領(lǐng)導(dǎo)節(jié)點為DG14：迭代次數(shù)、與收斂系數(shù)的關(guān)系曲線Fig.6 Graph of the relationships between iterations/and with leader DG14

由圖5及圖6可知，在收斂系數(shù)ε由0.002逐漸增加到1的過程中，收斂所需迭代次數(shù)曲線是一個U型曲線。當(dāng)ε增大到一定程度后（如圖6中ε＞0.504時），M矩陣譜半徑可能會大于1，此時算法不再收斂，系統(tǒng)失穩(wěn)。圖6中的兩種情況的譜半徑曲線同樣都是U型曲線，與各自的迭代次數(shù)曲線的走勢一致，并且在ρ（M）＞1時系統(tǒng)均失穩(wěn)，這也再次印證本文所提出的將ρ（M）作為收斂速度衡量因子的方法是切實有效的。

在實際確定收斂系數(shù)ε的數(shù)值時，可以采取在一定區(qū)間內(nèi)多次取值、計算相應(yīng)的ρ（M）的方法，最終采用U型曲線最低點對應(yīng)的數(shù)值為收斂系數(shù)ε的取值。

4.3 有/無通信延時仿真結(jié)果對比

分別將考慮通信延時的分布式算法式（9）以及如下不考慮通信延時的分布式算法式（21）應(yīng)用于通信網(wǎng)G-D1及G-D2中：

并依次選擇1～n號分布式電源作為領(lǐng)導(dǎo)節(jié)點，取收斂系數(shù)ε=0.1，分別記錄收斂所需的迭代次數(shù)，如圖7及圖8所示。其中，紅色圓形曲線對應(yīng)有通信延時一致性算法；藍(lán)色正方形曲線對應(yīng)無通信延時一致性算法；藍(lán)色三角形曲線對應(yīng)無通信延時一致性算法迭代所需收斂次數(shù)乘2。

圖7 通信網(wǎng)G-D1：有/無通信延時情況下算法收斂速度Fig.7 Convergence speed of algorithm with/without time-delay in communication network G-D1

圖8 通信網(wǎng)G-D2：有/無通信延時情況下算法收斂速度Fig.8 Convergence speed of algorithm with/without time-delay in communication network G-D2

在通信延時均勻、延時時長為τ的條件下，可以通過合理設(shè)置采樣時間T1的大小以消除延時因素對一致性算法的影響，此時應(yīng)滿足關(guān)系T1＞τ；若采用含通信延時的一致性算法，則采樣時間T2只需滿足關(guān)系T2＞τ/2。當(dāng)在滿足上述條件的前提下，采樣時間T1及T2取盡可能小的值時，近似地有關(guān)系T1=2T2。定義在有/無通信延時情況下的算法迭代所需收斂次數(shù)分別為k1和k2，則最終收斂所用時間t1=k1·T1=k1·2T2，因此比較2k1與k2的數(shù)值大小即可比較出最終收斂所用時間t1與t2的大小。

由圖7及圖8可知，對于通信網(wǎng)G-D1及G-D2，無通信延時的一致性算法收斂所需迭代次數(shù)均少于有通信延時的情況，這說明通信延時確實會導(dǎo)致算法迭代次數(shù)的增加，對于通信網(wǎng)G-D1而言，這樣的差異非常明顯（如圖7所示，有延時情況下的迭代次數(shù)基本在200以上，最高達(dá)1 159次，而無延時情況下的迭代次數(shù)全部在100次以下）。但需要注意的是，當(dāng)考慮采樣時間因素后，無通信延時的一致性算法反而可能比有通信延時的一致性算法收斂所需時間更長。例如，在圖8中，當(dāng)選取1～7號以及11～15號節(jié)點為領(lǐng)導(dǎo)節(jié)點時，均存在k1＜k2＜2k1，說明在這幾種情況下，采用有延時的一致性算法收斂所需時間反而更短。

綜上所述，對于以一致性算法為核心的分布式控制策略，對于有/無通信延時算法的選取，應(yīng)當(dāng)綜合考慮通信延時時間、通信網(wǎng)拓?fù)洹⑼ㄐ啪€路權(quán)重以及領(lǐng)導(dǎo)節(jié)點位置等因素。

5 結(jié)論

針對虛擬發(fā)電廠一致性控制策略中通信延時對收斂速度的影響進(jìn)行了研究，提出了均勻通信延時環(huán)境下的一致性算法收斂速度衡量因子。結(jié)果表明：

1）所提出的收斂速度衡量因子能夠較為準(zhǔn)確地反映通信延時環(huán)境下一致性算法的收斂性及收斂速度，并可以被用作領(lǐng)導(dǎo)節(jié)點選擇的指標(biāo)。

2）在設(shè)計通信網(wǎng)時，將收斂速度衡量因子最小化作為目標(biāo)對通信網(wǎng)線路結(jié)構(gòu)及權(quán)重進(jìn)行優(yōu)化，可以有效提升一致性算法的收斂速度。除此之外，本文還對線路權(quán)重、領(lǐng)導(dǎo)節(jié)點位置以及算法收斂系數(shù)大小對收斂速度的影響進(jìn)行了分析，并對比了有/無通信延時兩種情況下的收斂速度，為一致性算法通信網(wǎng)優(yōu)化設(shè)計提供了參考。

[1]許佳佳.基于虛擬發(fā)電廠的分布式電源調(diào)度管理模式的研究[D].北京：華北電力大學(xué)，2013.

[2]鮑薇，胡學(xué)浩，何國慶，等.分布式電源并網(wǎng)標(biāo)準(zhǔn)研究[J].電網(wǎng)技術(shù)，2012，36（11）：46-52.BAO Wei，HU Xuehao，HE Guoqing，et al.Study on stan －dard for grid-integration of distributed resources[J].Power System Technology，2012，36（11）：46-52（in Chinese）.

[3]REN W，BEARD R W.Distributed consensus in multivehicle cooperative control[M].London：Springer-Verlag，2008：1-13

[4]孟子陽.基于一致性算法的多航天器姿態(tài)協(xié)同控制研究[D].北京：清華大學(xué)，2010.

[5]REN W.Consensus strategies for cooperative control of vehicle formations[J].Control Theory&Applications，IET，2007，1（2）：505-512.

[6]何建平.基于一致性的無線傳感器網(wǎng)絡(luò)時鐘同步算法研究[D].杭州：浙江大學(xué)，2013.

[7]XIN H，GAN D，LI N，et al.Virtual power plant-based distributed control strategy for multiple distributed generators[J].Control Theory&Applications，IET，2013，7（1）：90-98.

[8]曹倩，宋永端，王磊，等.基于比率一致性算法的飛輪儲能矩陣系統(tǒng)分布式雙層控制[J].電網(wǎng)技術(shù)，2014，38（11）：3024-3029.CAO Qian，SONG Yongduan，WANG Lei，et al.Hierarchical distributed control for flywheel energy storage matrix system based on ratio consensus algorithm[J].Power System Technology，2014，38（11）：3024-3029（in Chinese）.

[9]XIN H，QU Z，SEUSS J，et al.A self-organizing strategy for power flow control of photovoltaic generators in a distribution network[J].Power Systems，IEEE Transactions on，2011，26（3）：1462-1473.

[10]張孝順，余濤.互聯(lián)電網(wǎng)AGC功率動態(tài)分配的虛擬發(fā)電部落協(xié)同一致性算法[J].中國電機工程學(xué)報，2015，35（15）：3750-3759.ZHANG Xiaoshun，YU Tao.Virtual generation tribe based collaborative consensus algorithm for dynamic generation dispath of AGC in interconnected power grids[J].Proceedings of the CSEE，2015，35（15）：3750-3759（in Chinese）.

[11]呂振宇，吳在軍，竇曉波，等.基于離散一致性的孤立直流微網(wǎng)自適應(yīng)下垂控制.中國電機工程學(xué)報，2015，35（17）：4397-4407.Lü Zhenyu，WU Zaijun，DOU Xiaobo，et al.An adaptive droop control for the islanded DC microgrid based on discrete consensus algorithm[J].Proceedings of the CSEE，2015，35（17）：4397-4407（in Chinese）.

[12]OLFATI SABER R，MURRAY R M.Consensus problems in networks of agents with switching topology and time-delays[J].Automatic Control，IEEE Transactions on， 2004，49（9）：1520-1533.

[13]LIN P，JIA Y.Consensus of second-order discrete-time multi-agent systems with nonuniform time-delays and dynamically changing topologies[J].Automatica，2009，45（9）：2154-2158.

[14]趙亦嵐.虛擬發(fā)電廠分布式控制收斂速度分析及優(yōu)化[D].杭州：浙江大學(xué)，2016.

[15]NUTKANI I U，POH CHIANG L，BLAABJERG F.Costbased droop scheme with lower generation costs for microgrids[J].Power Electronics，IET，2014，7（5）：1171-1180.[16]WOOD A J，WOLLENBERG B F.Power generation，operation，and control[M].John Wiley&Sons，2012.

[17]VARGA R S.Matrix iterative analysis[M].Springer Science&Business Media，2009.

[18]LIN X，BOYD S.Fast linear iterations for distributed averaging[J].Systems&Control Letters，2004，53（1）：65-78.

Convergence Speed Analysis for Distributed Control of Virtual Power Plant with Time-Delays

GUO Feng1，ZHAO Yilan2，3，YANG Huan3
（1.State Grid Shaanxi Electric Power Research Institute，Xi’an 710100，Shaanxi，China;2.Hangzhou Electric Power Design Institute，Hangzhou 310062，Zhejiang，China;3.Zhejiang University，Hangzhou 310027，Zhejiang，China）

In this paper，a control strategy with time-delays is applied to the virtual power plant control，and mathematical analysis is presented to reveal the factors affecting convergence speed of the control strategy.Moreover，a convergence speed factor is proposed.To offer advice for communication network design，the influence of communication line weight，leader's position and feedback coefficient on convergence speed is analyzed through simulation，convergence speed comparison with and without time-delays is also presented.

distributed control;virtual power plant;timedelays;convergence speed

1674-3814（2017）07-0077-08

TM721

A

國家高技術(shù)研究發(fā)展計劃（863計劃）（2015-AA050202）；浙江省自然科學(xué)基金（LY15E070003）。

Project Supported by the National High Technology Research and Development of China （863 Programme）（2015AA050202）；the Zhejiang Provincial Natural Science Foundation of China（LY15E070003）．

2017-04-21。

郭 峰（1964—），男，碩士，高級工程師，研究方向為電網(wǎng)運行與維護(hù)管理。

（編輯 徐花榮）