祝洪博，劉新偉，邵萌萌，賈建夫，楊龑亮

(1.東北電力大學，吉林省吉林市 132012;2.吉林省電力公司培訓中心，吉林長春 130022）

基于注入電流不平衡量的配電網(wǎng)改進潮流算法

祝洪博1，劉新偉1，邵萌萌1，賈建夫2，楊龑亮1

(1.東北電力大學，吉林省吉林市 132012;2.吉林省電力公司培訓中心，吉林長春 130022）

本文在節(jié)點注入電流模型配電網(wǎng)潮流計算的基礎上，采用了基于注入電流的不平衡量，對算法進行了改進和簡化，將Jacobian矩陣作為一個固定矩陣，從而使計算量大大減少，提高了計算速度;通過引入最優(yōu)乘子，從算法上保證了潮流計算的不發(fā)散，有效地解決了潮流計算對初值的敏感性以及一些病態(tài)系統(tǒng)的潮流計算問題。用30節(jié)點、141節(jié)點配電網(wǎng)算例和幾種典型的病態(tài)系統(tǒng)的計算驗證了算法的有效性和可行性。

配電網(wǎng);潮流計算;最優(yōu)乘子;病態(tài)

引言

配電網(wǎng)潮流計算是配電網(wǎng)經(jīng)濟運行、系統(tǒng)分析等的重要基礎［1］。目前，輸電系統(tǒng)潮流計算方法已較為成熟而且獲得了廣泛的實際應用。但由于配電網(wǎng)與輸電網(wǎng)有著明顯的差異:配電網(wǎng)通常呈輻射狀，支路比值較大，分支線較多;傳統(tǒng)的牛頓法和快速分解法在應用于配電網(wǎng)潮流計算時容易形成病態(tài)而無法收斂，因此，研究適合于配電網(wǎng)的潮流算法是至關重要的。

許多學者結(jié)合配電網(wǎng)特殊的網(wǎng)絡結(jié)構(gòu)，研究開發(fā)出一些適合于低壓配電網(wǎng)的母線類和支路類潮流算法，例如母線類的 Zbus［2］法和 Ybus［3］法，支路類的前推回推法、基于回路方程的算法、改進的牛頓-拉夫遜法等?；诠?jié)點注入電流模型的配電網(wǎng)潮流計算［4］是以節(jié)點注人電流為變量列出潮流方程，它同牛頓法相似，在每次迭代時都需要重新計算Jacobian矩陣，使計算速度大大降低，同時，配電網(wǎng)有時會因初值選取不當，無功緊張或其它原因?qū)е码妷嘿|(zhì)量很差，在進行潮流計算時會出現(xiàn)病態(tài)的問題。針對上述的問題，本文采用基于注入電流的不平衡量，對算法進行了改進和簡化，將Jacobian矩陣作為一個固定矩陣，從而使計算量大大減少，提高了計算速度;通過引入最優(yōu)乘子，有效地解決了配電網(wǎng)中病態(tài)潮流的計算問題。

1 節(jié)點注入電流模型的改進

對于單電源的配電網(wǎng)，取該電源點為平衡節(jié)點，其余節(jié)點為PQ節(jié)點;而對于多電源的配電網(wǎng)，選取其中的一個節(jié)點為平衡節(jié)點，另外的電源節(jié)點為PV節(jié)點，剩余的節(jié)點為PQ節(jié)點。這里節(jié)點電壓和注入功率以及線路導納都由直角坐標來表示。

(1）對PQ節(jié)點的處理

在文獻［4］的基礎上，引入了注入電流的不平衡量 ΔIi，設 ΔIi= Δai+jΔbi，將實部，虛部分開，則可得基于節(jié)點注入電流模型的非線性潮流方程。(式中Δai和Δbi為注入電流的不平衡量）。

再對 Δai，Δbi的非線性方程(即式1，2）運用梯度法求解，則可得潮流的修正方程為:

其中，Δei，Δfi為除平衡節(jié)點外的節(jié)點電壓增量的列向量。

式(3）中Jacobian矩陣各元素為:

當i=j時量級為10－1～10－2，+的值的數(shù)量級為10－1～100。因此可以認為式(4）、(5）、(6）、(7）中等號右邊后半部分的絕對值遠大于前面部分的絕對值，因此，可以對以上四個式子進行簡化:當i=j時，

在實際配電網(wǎng)潮流計算中，式(4）、(5）、(6）、(7）中等號右邊的兩部分的數(shù)值是相差很大的，通過許多例子可以發(fā)現(xiàn)二者的差距為1-3個數(shù)量級。一般來說，系統(tǒng)的線路電阻和電抗標幺值的數(shù)量級為10－3～10－1，從而 Gii和 Bii的值的標幺值的數(shù)量級為101～103。而負荷節(jié)點的功率值Pi和Qi的數(shù)

(2）對PV節(jié)點的處理

對于多電源的配電網(wǎng)，選取其中的一個節(jié)點為平衡節(jié)點，另外的電源節(jié)點為PV節(jié)點。對PV節(jié)點來說，由于無功注入量未知，無法得到ΔIi，因此，這里仍采用常規(guī)的牛頓法進行處理。

2 病態(tài)潮流問題的解決

在實際配電網(wǎng)潮流計算中，經(jīng)常會因初值設得不合理，無功緊張或其它原因?qū)е码妷嘿|(zhì)量很差、有重載線路而導致節(jié)點間相角差很大或其它原因，往往會出現(xiàn)計算過程的發(fā)散和振蕩的現(xiàn)象［5-6］。因此，引入了最優(yōu)乘子［7］，就是在第k次迭代求得狀態(tài)變量的修正量Δx(k）之后，不直接用Δx(k）去修正x(k），而是乘以一個標量乘子μ后再去修正，即用μ修正后的值作為k+1次迭代的初值。同時μ的選取要滿足:使潮流的誤差方程的平方和取最小值。這在某種程度上限制了對x(k）的修正，避免產(chǎn)生過修正或者欠修正，保證了潮流計算的不發(fā)散。因此，它能較好地處理配電網(wǎng)中病態(tài)潮流的計算問題。

采用功率型直角坐標的潮流方程的泰勒展開式可以精確地表示為

對于采用注入電流不平衡量的PQ節(jié)點［8］來說，此時上式中等號右端第三項中的各元素將為零，因此該項中僅對應于PV節(jié)點的元素有值，其余元素均為零。這樣構(gòu)造的目標函數(shù)為

通過對F(x）對μ求導，并令其等于零;

將上式展開，可得

式(16）可用牛頓法或卡丹(Cardan）公式來求解，所得解μ值即為最優(yōu)乘子μ*。

改進算法的計算流程如圖1所示。

其中各變量的定義如下:

3 算例分析

圖1 改進算法的流程圖Fig.1 Flow chart of improved algorithm

利用本文方法和文獻［4］算法分別對單電源的30節(jié)點配電網(wǎng)算例進行了潮流計算，計算結(jié)果比較如表1，原始數(shù)據(jù)參考文獻［9］。對多電源的配電網(wǎng)來說，選取雙電源的141節(jié)點配電網(wǎng)進行了潮流計算，計算結(jié)果比較如表2，原始數(shù)據(jù)參考文獻［4］。本文對30節(jié)點的初值選取收斂性也進行了比較如表3。這里的算法程序通過VB6.0來實現(xiàn)。還對一般方法難以收斂的病態(tài)系統(tǒng)(如13節(jié)點系統(tǒng)和43節(jié)點系統(tǒng)）進行了潮流計算，計算結(jié)果如表4。

表1 單電源30節(jié)點配電網(wǎng)的計算結(jié)果Tab.1 30-node 1-source distribution system

表2 雙電源141節(jié)點配電網(wǎng)的計算結(jié)果Tab.2 141-node 2-source distribution system

表3 對30節(jié)點初值選取的收斂性比較Tab.3 Convergence comparison for 30-node selected initial

表4 對處理病態(tài)系統(tǒng)的計算結(jié)果Tab.4 Calculation results for ill-conditioned system processing

由表1可以看出，對單電源的30節(jié)點系統(tǒng)來說，本文算法同牛頓-拉夫遜法相比，大大節(jié)省了內(nèi)存，而同文獻［4］相比，也節(jié)省了內(nèi)存，雖然收斂次數(shù)增加，但每次收斂的時間卻大大減少，從而總的收斂時間也減少了。

由表2的計算結(jié)果表明:對雙電源的141節(jié)點的系統(tǒng)來說，由于存在PV節(jié)點，對它的處理采用常規(guī)的牛頓法處理，每次迭代時只需更新與PV節(jié)點有關的Jacobian矩陣元素即可，再加上在實際配電網(wǎng)中PV節(jié)點較少，因此，雖然收斂次數(shù)增加，但它總的收斂時間比牛頓法和文獻［4］都要短，而且也節(jié)省了內(nèi)存。

由表3表明:當牛頓-拉夫遜法的節(jié)點電壓初值低于0.5時，潮流開始發(fā)散(或振蕩）;而文獻［4］算法的節(jié)點電壓初值低于0.3時，潮流也開始發(fā)散(或振蕩）;本文的算法由于引入了最優(yōu)乘子［10］，算法穩(wěn)定性比牛頓法和文獻［4］算法都好，當節(jié)點電壓初值低于0.1時，仍能可靠地收斂。

由表4可以看出:對一些病態(tài)系統(tǒng)(如13和43節(jié)點系統(tǒng)），牛頓法和文獻［4］算法均難以收斂，而本文的算法卻能保證潮流可靠地收斂。

4 結(jié)論

本文充分利用了配電網(wǎng)不同于輸電網(wǎng)的特點，在基于電流模型的配電網(wǎng)潮流計算的基礎上進行了簡化和改進，不僅提高了計算速度，也節(jié)省了內(nèi)存，而且通過引入最優(yōu)乘子，保證了潮流計算的不發(fā)散，有效地解決了病態(tài)潮流的計算問題，可有效地運用于復雜單電源和多電源配電網(wǎng)的潮流計算。同牛頓-拉夫遜法和文獻［4］算法相比，節(jié)省了內(nèi)存，雖然收斂次數(shù)增加，但每次收斂的時間卻大大減少，從而總的收斂時間也減少了，從而提高了計算速度;并且解決了潮流計算對初值的敏感性以及一些病態(tài)系統(tǒng)的潮流計算問題，通過算例證明了該算法的有效性和可行性。

References）:

［1］張學松，柳焯(Zhang Xuesong，Liu Zhuo）.配電網(wǎng)潮流算法比較研究 (Comparative study on power flow algorithm of distribution network）［J］.電網(wǎng)技術 (Power System Tech.），1998，22(4）:96-100.

［2］Chen T S，Chen M S.Distribution system power flow analysis a rigid approach［J］.IEEE Trans.PWD，1991，6(3）:1146-1152.

［3］Tripathy S C，Prasad G D，Malik O P.Load-flow solutions for ill-conditioned power system by a Newton-like method［J］.IEEE Trans.on PAS，1982，101(10）:3648-3657.

［4］鄢長春，張焰，陳章潮 (Yan Changchun，Zhang Yan，Chen Zhangchao）.基于節(jié)點注入電流模型的配電網(wǎng)潮流算法(Distribution system power flow calculation based on the node current-injection model）［J］.電力系統(tǒng)自動化 (Automation of Elec.Power System），1999，23(7）:55-58.

［5］諸駿偉(Zhu Junwei）.電力系統(tǒng)分析(Power system analysis）［M］.北京:中國電力出版社 (Beijing:China Elec.Power Press），1995.

［6］劉廣一，胡錫龍，于爾鏗 (Liu Guangyi，Hu Xilong，Yu Erkeng）.電力系統(tǒng)病態(tài)潮流計算新算法(A new calculation algorithm of power system ill condition power flow）［J］.中國電機工程學報 (Proc.CSEE），1991，11(S1）:27-36.

［7］Iwamoto S，Tamura Y.A fast load flow method retaining nonlinearity［J］.IEEE Trans.on PAS，1978，100(4）:1736-1737.

［8］王承民，蔣傳文，侯志儉(Wang Chengmin，Jiang Chuanwen，Hou Zhijian）.基于節(jié)點不平衡功率的病態(tài)潮流算法(Ill-conditioned load flow algorithm based on node lopsided powers）［J］.上海交通大學學報 (J Shanghai Jiaotong Univ.），2004，38(8）:24-28.

［9］K Cespedes.New method for analysis of distribution networks［J］.IEEE Trans.on Power Delivery，1990，1(22）:1736-1737.

［10］石飛，趙晉全，王毅 (Shi Fei，Zhao Jinquan，Wang Yi）.計及發(fā)電機無功約束的最優(yōu)乘子潮流計算方法的比較(A comparison of Newton-raphson optimal multipliers load flow methods condisidering generator reactive power limitation）［J］.電力系統(tǒng)保護與控制 (Power System Protection＆ Control），2008，37(2）:54-59.

［11］盧明富，梅生偉(Lu Mingfu，Mei Shengwei）.小世界電網(wǎng)生長演化模型及其潮流特性分析 (Evolutive model of small-world power grid and its power flow analysis）［J］.電工電能新技術 (Adv.Tech.of Elec.Eng.＆Energy），2010，29(1）:25-29.

Improved power flow calculation of distribution network based on injection current imbalance

ZHU Hong-bo1，LIU Xin-wei1，SHAO Meng-meng1，JIA Jian-fu2，YANG Yan-liang1
(1.Northeast Dian Li University，Jilin 132012，China;2.Training Center in Jilin Province Power Company，Changchun 130022，China）

Based on the power flow calculation of distribution network of the node current-injection model，the injection current imbalance is adopted to improve and simplify the algorithm，to make Jacobian matrix a fixed matrix，and computation is greatly reduced and the computing speed is improved.An optimal factor is introduced to ensure the convergence of the computing process，and the sensitivity of the initial value about power flow calculation and some ill-conditioned systems problems are also resolved effectively.Using the distribution network of 30-node，141-node and several typical examples of ill-conditioned systems verify the effectiveness and feasibility of the method.

distribution network;power flow calculation;optimal factor;ill-conditioned

TM744

A

1003-3076(2010）03-0022-04

2009-08-24

祝洪博(1942-），男，吉林籍，教授，研究方向為配電自動化、人工智能等;

劉新偉(1983-），男，山東籍，碩士研究生，主要研究方向為配電自動化。