劉暢　張志　康麗

摘 要：電力系統(tǒng)潮流計算是指電力系統(tǒng)運行方式處于穩(wěn)定狀態(tài)時，根據(jù)已知運行條件和系統(tǒng)結(jié)構(gòu)確定其運行狀態(tài)的計算方法。由于其涉及算法迭代的思想和繁瑣的矩陣變換，潮流計算是《電力系統(tǒng)分析》課程教學(xué)過程中的難點和重點。本文針對潮流計算中牛頓-拉夫遜算法，以IEEE6-bus為實例，采用Matlab軟件進(jìn)行潮流計算，取得可視化的收斂結(jié)果，為潮流計算的教學(xué)過程提供了新的思路。

關(guān)鍵詞：電力系統(tǒng)；潮流計算；Matlab仿真

中圖分類號：TM744 文獻(xiàn)標(biāo)志碼：B 文章編號：2095-9214（2016）01-0135-02

引言

潮流計算是《電力系統(tǒng)分析》教學(xué)過程中的重難點，該部分內(nèi)容涉及迭代算法和矩陣變化的應(yīng)用，其繁瑣的數(shù)學(xué)推導(dǎo)過程往往降低了學(xué)生的學(xué)習(xí)興趣與積極性[1]。潮流計算教學(xué)通常包含牛頓-拉夫遜潮流算法和P-Q分解潮流算法，兩種方法的核心都是將對非線性功率方程組轉(zhuǎn)變?yōu)閷€性微分方程組的求解；其先假定一個近似值，與真值的誤差為修正量，構(gòu)建修正方程，它的系數(shù)矩陣為雅克比矩陣，不斷地求解修正方程，用求得的修正變量反復(fù)修正近似解，使其越來越接近真實解，直至滿足收斂精度 [2]。

針對潮流算法的教學(xué)的繁瑣性，近年來，市面上推出了眾多潮流計算的教學(xué)軟件，其模塊化與可視化程度較高，并可以靈活的通過軟件中的各種電氣組件，自由搭建不同的電力系統(tǒng)，最終獲得系統(tǒng)準(zhǔn)確的潮流分布[3]。然而，此類軟件的最大缺陷是其忽略了潮流算法的整個流程，只通過組件的搭建獲得結(jié)果，雖然方便了教學(xué)過程，但影響了學(xué)生對潮流算法的深入理解與主動的探索。

按學(xué)校的教學(xué)安排，在潮流計算中引入Matlab仿真教學(xué)，引導(dǎo)學(xué)生通過此種仿真軟件，靈活的搭建系統(tǒng)并設(shè)計算法過程。本文在此思路下，針對典型電力系統(tǒng)的仿真實例，獲得可視化的收斂圖形結(jié)果，進(jìn)一步為潮流計算的深化教學(xué)奠定基礎(chǔ)。

一、節(jié)點分類與基本方程

1.節(jié)點分類。電力系統(tǒng)中存在許多節(jié)點，按照節(jié)點已知條件的不同，一般可將節(jié)點分為三類：

（1）PQ節(jié)點。這類節(jié)點已知節(jié)點的有功功率Pi和無功功率Qi，待求量是節(jié)點電壓幅值i和相位δi。變電所通常都是屬于PQ節(jié)點，PQ節(jié)點還包括在某些時段發(fā)出穩(wěn)定功率的發(fā)電廠，系統(tǒng)中的負(fù)荷節(jié)點以及有功和無功都為零的浮游節(jié)點。

（2）PU節(jié)點。這類節(jié)點已知節(jié)點的有功功率Pi和電壓幅值i，待求節(jié)點的無功功率Qi和電壓相位δi。這類節(jié)點一般要求選擇在有一定的無功補償設(shè)備或者無功功率儲備的母線，這樣可以維持電壓的恒定。

（3） 平衡節(jié)點。潮流計算之中，只用一個節(jié)點來平衡網(wǎng)絡(luò)中的有功功率，它就是平衡節(jié)點。在一個普通網(wǎng)絡(luò)中，各支路的功率損耗一般是未知的，這就需要至少有一個承擔(dān)著平衡系統(tǒng)有功功率的節(jié)點，其電壓幅值i和電壓相位δi已知，待求節(jié)點的有功功率Pi和無功功率Qi。在實際系統(tǒng)中主調(diào)頻發(fā)電廠或出線數(shù)最多的發(fā)電廠都可作為平衡節(jié)點。

2.基本方程。若電氣網(wǎng)絡(luò)的節(jié)點數(shù)為n，其n個節(jié)點的電力系統(tǒng)潮流方程的一般形式可表示為：

將上式的實部和虛部展開，每個節(jié)點可得到2個實數(shù)方程；根據(jù)1.1中節(jié)點的節(jié)點分類，每個節(jié)點的2個實數(shù)方程對應(yīng)2個待求變量，其總共的方程與待求變量的個數(shù)均為2n；采用直角坐標(biāo)表示節(jié)點電壓和導(dǎo)納，即i=ei+jfi、ij=Gij+jBij，則展開式為：

上式為直接坐標(biāo)系下的功率基本方程。

二、修正方程

若系統(tǒng)總節(jié)點數(shù)為n，其中PQ節(jié)點個數(shù)為m，PU節(jié)點個數(shù)為n-m-1，平衡節(jié)點個數(shù)為1；則對于PQ節(jié)點 （i=1，2，…，m-1），根據(jù)式（2）可得功率不平衡量為：

根據(jù)以上分析，對于節(jié)點數(shù)為n的系統(tǒng)，可列出2n-2個非線性方程，其中有功功率方程n-1個，無功功率方程m-1個，電壓方程n-m。將式（3）、（4）聯(lián)立形成n節(jié)點系統(tǒng)的非線性方程組，并按泰勒級數(shù)在電壓初始值e（0）i和f（0）i （i=1，2，…，n，i≠m）處展開，忽略高次項后，可得修改方程的矩陣形式為：

照上述方法反復(fù)迭代，直到達(dá)到預(yù)設(shè)的收斂精度后迭代停止。收斂條件一般取所有節(jié)點的ΔP≤ε和ΔQ≤ε或ΔU2≤ε，ε為很小的數(shù)[4]。

三、仿真結(jié)果

如圖1所示，潮流計算的仿真采用IEEE-6BUS標(biāo)準(zhǔn)試驗系統(tǒng)，網(wǎng)絡(luò)中的各類節(jié)點如下：（1）PQ節(jié)點，網(wǎng)絡(luò)中節(jié)點標(biāo)號為前m個的節(jié)點為PQ節(jié)點，現(xiàn)在系統(tǒng)中以節(jié)點1至節(jié)點4為PQ節(jié)點；（2）PU節(jié)點，記為5號節(jié)點，這類節(jié)點的有功功率和電壓幅值是已給定；（3）平衡節(jié)點，它的電壓幅值和電壓相位是常數(shù)，記為6號節(jié)點。

直角坐標(biāo)程序仿真時，收斂精度為ε=10-5，圖2至圖5分別表示4個PQ節(jié)點迭代過程有功和無功的不平衡量的變化過程；圖6表示PU節(jié)點迭代過程電壓和無功的不平衡量的變化過程；圖7表示6個節(jié)點在迭代過程中電壓幅值的變化過程，其中平衡節(jié)點的電壓幅值一直保持固定值；從以上各節(jié)點收斂曲線中可看出，經(jīng)過四至五次算法迭代后，各節(jié)點的有功功率、無功功率或電壓幅值均達(dá)到收斂要求。

四、結(jié)束語

電力系統(tǒng)潮流計算是《電力系統(tǒng)分析》課程教學(xué)內(nèi)容中的重要組部分，是電力系統(tǒng)安全運行和經(jīng)濟(jì)規(guī)劃的基礎(chǔ)。本文分析了基于直角坐標(biāo)系的牛頓-拉夫遜電力系統(tǒng)潮流計算方法，以IEEE6-Bus為實例，采用Matlab軟件對其進(jìn)行的仿真，充分利用了Matlab的運算功能，豐富了教學(xué)手段的多樣性，同時也適合學(xué)生自己獨立設(shè)計個性化的系統(tǒng)和完成實驗項目，符合當(dāng)前電力系統(tǒng)潮流分析數(shù)值計算的發(fā)展趨勢。

（作者單位：東莞理工學(xué)院電子工程學(xué)院）

基金項目：東莞理工學(xué)院教學(xué)改革基金資助項目。

參考文獻(xiàn)：

[1] 何仰贊，溫增銀.電力系統(tǒng)分析[M]. 第三版. 武漢： 華中科技大學(xué)出版社，2002.

[2] 陳帥，王勇，楊恒. 基于MATLAB和PSASP的電力系統(tǒng)潮流分析與計算[J].上海電力學(xué)院學(xué)報. 2012，28（1）：19-22.

[3] 廖青華，趙斌. 電力系統(tǒng)可視化仿真軟件PowerWorld的應(yīng)用[J].上海電力學(xué)院學(xué)報. 2008，16（6）：18-19.

[4] 胡浩，閏英敏，陳永利. 基于MATLAB的電力系統(tǒng)潮流計算[J].應(yīng)用天地. 2012， 31（12） ：55-60.