研究生學(xué)位課電力系統(tǒng)潮流計(jì)算1

1、north china electric power universitydepartment of electrical engineeringbaoding2009.11-2009.01 目目 錄錄 一電力系統(tǒng)潮流計(jì)算一電力系統(tǒng)潮流計(jì)算二電力系統(tǒng)狀態(tài)估計(jì)二電力系統(tǒng)狀態(tài)估計(jì) 四電力系統(tǒng)復(fù)雜故障分析四電力系統(tǒng)復(fù)雜故障分析 三三電力系統(tǒng)靜態(tài)安全分析電力系統(tǒng)靜態(tài)安全分析參考書參考書 1. 現(xiàn)代電力系統(tǒng)分析現(xiàn)代電力系統(tǒng)分析王錫凡王錫凡 主編主編2. 高等電力網(wǎng)分析高等電力網(wǎng)分析 張伯明張伯明 3. 電力系統(tǒng)狀態(tài)估計(jì)電力系統(tǒng)狀態(tài)估計(jì)于爾鏗于爾鏗 主編主編 4 電力系統(tǒng)故障分析電力系統(tǒng)故障分析劉萬順劉萬

2、順第一章第一章 電力系統(tǒng)潮流計(jì)算電力系統(tǒng)潮流計(jì)算一概述一概述 二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 三三潮流計(jì)算的幾種基本方法潮流計(jì)算的幾種基本方法 四四保留非線性潮流算法保留非線性潮流算法 五最小化潮流算法五最小化潮流算法六六潮流計(jì)算中的自動(dòng)調(diào)整潮流計(jì)算中的自動(dòng)調(diào)整 七最優(yōu)潮流問題七最優(yōu)潮流問題八八交直流電力系統(tǒng)的潮流計(jì)算交直流電力系統(tǒng)的潮流計(jì)算 九九幾種特殊性質(zhì)的潮流計(jì)算問題簡(jiǎn)介幾種特殊性質(zhì)的潮流計(jì)算問題簡(jiǎn)介 電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的基本電氣計(jì)算，電力系穩(wěn)態(tài)運(yùn)行情況的基本電氣計(jì)算，電力系統(tǒng)潮流計(jì)算的任務(wù)是根據(jù)給定的網(wǎng)絡(luò)結(jié)統(tǒng)潮流計(jì)

3、算的任務(wù)是根據(jù)給定的網(wǎng)絡(luò)結(jié)構(gòu)及運(yùn)行條件，求出電網(wǎng)的運(yùn)行狀態(tài)，構(gòu)及運(yùn)行條件，求出電網(wǎng)的運(yùn)行狀態(tài)，其中包括各母線的電壓、各支路的功率其中包括各母線的電壓、各支路的功率分布以及功率損耗等。分布以及功率損耗等。一概述一概述離線計(jì)算離線計(jì)算: :規(guī)劃設(shè)計(jì)規(guī)劃設(shè)計(jì); ;運(yùn)行方式分析運(yùn)行方式分析; ;其他其他計(jì)算的配合計(jì)算的配合在線計(jì)算在線計(jì)算: :安全監(jiān)控和安全分析安全監(jiān)控和安全分析潮流計(jì)算是電力系統(tǒng)中應(yīng)用最為廣泛、潮流計(jì)算是電力系統(tǒng)中應(yīng)用最為廣泛、最基本和最重要的一種電氣計(jì)算。最基本和最重要的一種電氣計(jì)算。一概述一概述 常用的潮流計(jì)算方法歸納到數(shù)學(xué)上屬常用的潮流計(jì)算方法歸納到數(shù)學(xué)上屬于多元非線性代數(shù)方程

4、組的求解問題，于多元非線性代數(shù)方程組的求解問題，一般需采用迭代計(jì)算方法進(jìn)行求解計(jì)算。一般需采用迭代計(jì)算方法進(jìn)行求解計(jì)算。 2020世紀(jì)世紀(jì)5050年代中期起，電力系統(tǒng)潮流年代中期起，電力系統(tǒng)潮流計(jì)算的研究就是如何使用電子計(jì)算機(jī)計(jì)計(jì)算的研究就是如何使用電子計(jì)算機(jī)計(jì)算電力系統(tǒng)的潮流問題。算電力系統(tǒng)的潮流問題。 一概述一概述 對(duì)于潮流算法，其基本要求可歸納成對(duì)于潮流算法，其基本要求可歸納成以下四個(gè)方面：以下四個(gè)方面： （1 1）計(jì)算速度；）計(jì)算速度； （2 2）計(jì)算機(jī)內(nèi)存占用量；）計(jì)算機(jī)內(nèi)存占用量； （3 3）算法的收斂可靠性；）算法的收斂可靠性； （4 4）程序設(shè)計(jì)的方便性以及算法擴(kuò)充）程序設(shè)計(jì)的

5、方便性以及算法擴(kuò)充移植等的靈活通用性。移植等的靈活通用性。 此外，程序使用的方便性及良好的人此外，程序使用的方便性及良好的人- -機(jī)界面也越來越受到人們的關(guān)注。機(jī)界面也越來越受到人們的關(guān)注。一概述一概述 本章安排本章安排: : 潮流計(jì)算問題的數(shù)學(xué)模型潮流計(jì)算問題的數(shù)學(xué)模型 三種最基本的潮流算法三種最基本的潮流算法 保留非線性潮流法保留非線性潮流法 最小化潮流計(jì)算最小化潮流計(jì)算 自動(dòng)調(diào)整計(jì)算功能自動(dòng)調(diào)整計(jì)算功能 最優(yōu)潮流最優(yōu)潮流 特殊用途的潮流計(jì)算問題特殊用途的潮流計(jì)算問題 一概述一概述 電力系統(tǒng)由發(fā)電機(jī)、變壓器、輸配電電力系統(tǒng)由發(fā)電機(jī)、變壓器、輸配電線路及負(fù)荷等組成。線路及負(fù)荷等組成。 進(jìn)行潮

6、流計(jì)算時(shí)，發(fā)電機(jī)和負(fù)荷一般進(jìn)行潮流計(jì)算時(shí)，發(fā)電機(jī)和負(fù)荷一般可用接在相應(yīng)節(jié)點(diǎn)上的一個(gè)電流注入量可用接在相應(yīng)節(jié)點(diǎn)上的一個(gè)電流注入量表示。表示。 電力網(wǎng)絡(luò)中的變壓器、線路、電容器、電力網(wǎng)絡(luò)中的變壓器、線路、電容器、電抗器等元件可用集中參數(shù)表示的由線電抗器等元件可用集中參數(shù)表示的由線性電阻、電抗構(gòu)成的等值電路模擬。性電阻、電抗構(gòu)成的等值電路模擬。 二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 對(duì)這樣的線性網(wǎng)絡(luò)一般采用節(jié)點(diǎn)電壓對(duì)這樣的線性網(wǎng)絡(luò)一般采用節(jié)點(diǎn)電壓法進(jìn)行分析。節(jié)點(diǎn)電壓與節(jié)點(diǎn)注入電流法進(jìn)行分析。節(jié)點(diǎn)電壓與節(jié)點(diǎn)注入電流之間的關(guān)系為之間的關(guān)系為: :或或二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)

7、學(xué)模型 iuyizu 展開為展開為或或二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 niiuyinjjij, 2 , 11niizunjjiji, 2 , 11 在實(shí)際中，已知的節(jié)點(diǎn)注入量往往不在實(shí)際中，已知的節(jié)點(diǎn)注入量往往不是節(jié)點(diǎn)電流而是節(jié)點(diǎn)功率，為此用節(jié)點(diǎn)是節(jié)點(diǎn)電流而是節(jié)點(diǎn)功率，為此用節(jié)點(diǎn)功率代替節(jié)點(diǎn)電流功率代替節(jié)點(diǎn)電流, ,得得 (1-6) 或或 (1-7) 二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 niujqpuyiiinjjij, 2 , 11niujqpzunjjjjiji, 2 , 11 上兩上兩式式是潮流計(jì)算問題的基本方程式，是潮流計(jì)算問題的基本方程式，是一個(gè)以節(jié)點(diǎn)

8、電壓為變量的非線性代數(shù)是一個(gè)以節(jié)點(diǎn)電壓為變量的非線性代數(shù)方程組。而采用節(jié)點(diǎn)功率作為節(jié)點(diǎn)注入方程組。而采用節(jié)點(diǎn)功率作為節(jié)點(diǎn)注入量是造成方程組呈非線性的根本原因。量是造成方程組呈非線性的根本原因。由于方程組為非線性的，因此必須采用由于方程組為非線性的，因此必須采用迭代方法進(jìn)行數(shù)值求解。迭代方法進(jìn)行數(shù)值求解。 根據(jù)對(duì)方程組的不同處理方式，形成根據(jù)對(duì)方程組的不同處理方式，形成了不同的潮流算法。了不同的潮流算法。二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 對(duì)于電力系統(tǒng)中的每個(gè)節(jié)點(diǎn)，需要對(duì)于電力系統(tǒng)中的每個(gè)節(jié)點(diǎn)，需要p p、q q 、u u和相角四個(gè)變量才能確定其運(yùn)行和相角四個(gè)變量才能確定其運(yùn)行狀態(tài)

9、。狀態(tài)。n n個(gè)節(jié)點(diǎn)總共有個(gè)節(jié)點(diǎn)總共有4n4n個(gè)運(yùn)行變量。個(gè)運(yùn)行變量。而基本方程式只有而基本方程式只有n n個(gè)個(gè), ,將實(shí)部與虛部分將實(shí)部與虛部分開，則形成開，則形成2n2n個(gè)實(shí)數(shù)方程式，僅可解得個(gè)實(shí)數(shù)方程式，僅可解得2n2n個(gè)未知運(yùn)行變量。必須將另外個(gè)未知運(yùn)行變量。必須將另外2n2n個(gè)變個(gè)變量作為已知量而預(yù)先給定。也即對(duì)每個(gè)量作為已知量而預(yù)先給定。也即對(duì)每個(gè)節(jié)點(diǎn)，要給定兩個(gè)變量為已知條件，而節(jié)點(diǎn)，要給定兩個(gè)變量為已知條件，而另兩個(gè)變量作為待求量。另兩個(gè)變量作為待求量。二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 根據(jù)電力系統(tǒng)的實(shí)際運(yùn)行條件，按照根據(jù)電力系統(tǒng)的實(shí)際運(yùn)行條件，按照預(yù)先給定的變

10、量的不同，電力系統(tǒng)的節(jié)預(yù)先給定的變量的不同，電力系統(tǒng)的節(jié)點(diǎn)可分成點(diǎn)可分成pqpq節(jié)點(diǎn)、節(jié)點(diǎn)、pvpv節(jié)點(diǎn)及平衡節(jié)點(diǎn)三節(jié)點(diǎn)及平衡節(jié)點(diǎn)三種類型。種類型。 對(duì)平衡節(jié)點(diǎn)來說，其電壓相角一般作對(duì)平衡節(jié)點(diǎn)來說，其電壓相角一般作為系統(tǒng)電壓相角的基準(zhǔn)。為系統(tǒng)電壓相角的基準(zhǔn)。 二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型 交流電力系統(tǒng)中的復(fù)數(shù)電壓變量可交流電力系統(tǒng)中的復(fù)數(shù)電壓變量可以用兩種坐標(biāo)形式表示以用兩種坐標(biāo)形式表示或或 而復(fù)數(shù)導(dǎo)納為而復(fù)數(shù)導(dǎo)納為二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型jiieuuiiijfeuijijijjbgy 將以上三式代入以導(dǎo)納矩陣為基礎(chǔ)的將以上三式代入以導(dǎo)納矩陣為基礎(chǔ)的

11、式式(1-6)(1-6)，并將實(shí)部與虛部分開，可得，并將實(shí)部與虛部分開，可得到兩種形式的潮流方程。到兩種形式的潮流方程。二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型直角坐標(biāo)形式直角坐標(biāo)形式 (1-11) (1-12)極坐標(biāo)形式極坐標(biāo)形式 (1-13) (1-14) 二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型niebfgffbegepijijjijjijijijjijii, 2 , 1)()(niebfgefbegfqijijjijjijijijjijii, 2 , 1)()(nibguupijijijijijjii, 2 , 1)sincos(sincos)1,2,iijijijij

12、ijj iquugbin 若以若以p p、u u、x x分別表示擾動(dòng)變量、控分別表示擾動(dòng)變量、控制變量、狀態(tài)變量，則潮流方程可以用制變量、狀態(tài)變量，則潮流方程可以用更簡(jiǎn)潔的方式表示為更簡(jiǎn)潔的方式表示為 (1-15)(1-15) 根據(jù)式根據(jù)式(1-15)(1-15)，潮流計(jì)算的含義就是，潮流計(jì)算的含義就是針對(duì)某個(gè)擾動(dòng)變量針對(duì)某個(gè)擾動(dòng)變量p p，根據(jù)給定的控制，根據(jù)給定的控制變量變量u u，求出相應(yīng)的狀態(tài)變量，求出相應(yīng)的狀態(tài)變量x x。二潮流計(jì)算問題的數(shù)學(xué)模型二潮流計(jì)算問題的數(shù)學(xué)模型0),(puxf一一 高斯高斯- -塞德爾法塞德爾法 以導(dǎo)納矩陣為基礎(chǔ)，并應(yīng)用高斯以導(dǎo)納矩陣為基礎(chǔ)，并應(yīng)用高斯- -

13、塞塞德爾迭代的算法是電力系統(tǒng)應(yīng)用最早的德爾迭代的算法是電力系統(tǒng)應(yīng)用最早的潮流計(jì)算方法。潮流計(jì)算方法。 三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法 討論電力系統(tǒng)中除討論電力系統(tǒng)中除1 1個(gè)平衡節(jié)點(diǎn)外，其個(gè)平衡節(jié)點(diǎn)外，其余都是余都是pqpq節(jié)點(diǎn)的情況。節(jié)點(diǎn)的情況。 由式由式(1-6)(1-6)可得可得 (1-16)(1-16) 式中：式中： 、 為已知的節(jié)點(diǎn)注入有功、無為已知的節(jié)點(diǎn)注入有功、無功功率。功功率。三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法niuyujqpyunijjjijiiiiii, 2 , 111ipiq 假定節(jié)點(diǎn)假定節(jié)點(diǎn)l l為平衡節(jié)點(diǎn)，其給定電壓為平衡節(jié)點(diǎn)，其給定電

14、壓為為 。平衡節(jié)點(diǎn)不參加迭代。于是對(duì)應(yīng)。平衡節(jié)點(diǎn)不參加迭代。于是對(duì)應(yīng)于這種情況的高斯于這種情況的高斯- -塞德爾迭代格式為塞德爾迭代格式為 (1-17)(1-17) 上式是該算法最基本的迭代計(jì)算公式。上式是該算法最基本的迭代計(jì)算公式。 其迭代收斂的判據(jù)是其迭代收斂的判據(jù)是 三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法su111(1)( )1121( )1(2,3,inkskkiiiiijiijijj ikiiipjquy uy uy uinyu kikiiuu1max 本算法的突出優(yōu)點(diǎn)是原理簡(jiǎn)單，程序設(shè)本算法的突出優(yōu)點(diǎn)是原理簡(jiǎn)單，程序設(shè)計(jì)容易。導(dǎo)納矩陣對(duì)稱且高度稀疏，因計(jì)容易。導(dǎo)納矩陣對(duì)稱

15、且高度稀疏，因此占用內(nèi)存非常節(jié)省。此占用內(nèi)存非常節(jié)省。 該算法的主要缺點(diǎn)是收斂速度慢。該算法的主要缺點(diǎn)是收斂速度慢。由于各節(jié)點(diǎn)電壓在數(shù)學(xué)上松散耦合，所由于各節(jié)點(diǎn)電壓在數(shù)學(xué)上松散耦合，所以節(jié)點(diǎn)電壓向精確值的接近非常緩慢。以節(jié)點(diǎn)電壓向精確值的接近非常緩慢。另外，算法的迭代次數(shù)隨著網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)另外，算法的迭代次數(shù)隨著網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)的增加而上升，因此在用于較大規(guī)模電的增加而上升，因此在用于較大規(guī)模電力系統(tǒng)的潮流計(jì)算時(shí)，速度顯得非常緩力系統(tǒng)的潮流計(jì)算時(shí)，速度顯得非常緩慢。慢。三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法 為提高算法收斂速度，常用的方法是為提高算法收斂速度，常用的方法是在迭代過程中加入加速因

16、子在迭代過程中加入加速因子 ，即取，即取 式中：式中： 是通過式是通過式(1-17)(1-17)求得的節(jié)點(diǎn)求得的節(jié)點(diǎn)i i電壓的第電壓的第k+1k+1次迭代值；次迭代值； 是修正后是修正后節(jié)點(diǎn)節(jié)點(diǎn)i i電壓的第電壓的第k+1k+1次迭代值；次迭代值； 為加速為加速因子，一般取因子，一般取 。三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法 kikikikiuuuu111 kiu1kiu21 對(duì)于具有下述所謂病態(tài)條件的系統(tǒng)，高對(duì)于具有下述所謂病態(tài)條件的系統(tǒng)，高斯斯- -塞德爾迭代法往往會(huì)發(fā)生收斂困難：塞德爾迭代法往往會(huì)發(fā)生收斂困難： (l)(l)節(jié)點(diǎn)間相位角差很大的重負(fù)荷系統(tǒng)；節(jié)點(diǎn)間相位角差很大

17、的重負(fù)荷系統(tǒng)； (2)(2)包含有負(fù)電抗支路（如某些三繞組變包含有負(fù)電抗支路（如某些三繞組變壓器或線路串聯(lián)電容等）的系統(tǒng)；壓器或線路串聯(lián)電容等）的系統(tǒng)； (3)(3)具有較長(zhǎng)的輻射形線路的系統(tǒng)；具有較長(zhǎng)的輻射形線路的系統(tǒng)； (4)(4)長(zhǎng)線路與短線路接在同一節(jié)點(diǎn)上，而長(zhǎng)線路與短線路接在同一節(jié)點(diǎn)上，而且長(zhǎng)短線路的長(zhǎng)度比值又很大的系統(tǒng)。且長(zhǎng)短線路的長(zhǎng)度比值又很大的系統(tǒng)。 此外，選擇不同的節(jié)點(diǎn)為平衡節(jié)點(diǎn)，也此外，選擇不同的節(jié)點(diǎn)為平衡節(jié)點(diǎn)，也會(huì)影響到收斂性能。會(huì)影響到收斂性能。三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法 為克服基于節(jié)點(diǎn)導(dǎo)納矩陣的高斯為克服基于節(jié)點(diǎn)導(dǎo)納矩陣的高斯- -塞塞德爾迭代法

18、的這些缺點(diǎn)，德爾迭代法的這些缺點(diǎn)，2020世紀(jì)世紀(jì)6060年代年代初提出了基于節(jié)點(diǎn)阻抗矩陣的高斯初提出了基于節(jié)點(diǎn)阻抗矩陣的高斯- -塞塞德爾迭代法。但在牛頓法潮流出現(xiàn)后，德爾迭代法。但在牛頓法潮流出現(xiàn)后，即很少再被便用。即很少再被便用。 目前基于節(jié)點(diǎn)導(dǎo)納矩陣的高斯目前基于節(jié)點(diǎn)導(dǎo)納矩陣的高斯- -塞德塞德爾法主要為牛頓法等對(duì)于待求量的迭代爾法主要為牛頓法等對(duì)于待求量的迭代初值要求比較高的算法提供初值，一般初值要求比較高的算法提供初值，一般只需迭代只需迭代1 12 2次就可以滿足要求。次就可以滿足要求。 三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法二二 牛頓牛頓- -拉夫遜法拉夫遜法（一）牛頓

19、（一）牛頓- -拉夫遜法的一般概念拉夫遜法的一般概念牛頓牛頓- -拉夫遜法在數(shù)學(xué)上是求解非線性代拉夫遜法在數(shù)學(xué)上是求解非線性代數(shù)方程式的有效方法。其要點(diǎn)是把非線數(shù)方程式的有效方法。其要點(diǎn)是把非線性方程式的求解過程變成反復(fù)地對(duì)相應(yīng)性方程式的求解過程變成反復(fù)地對(duì)相應(yīng)的線性方程式進(jìn)行求解的過程，即通常的線性方程式進(jìn)行求解的過程，即通常所稱的逐次線性化過程。所稱的逐次線性化過程。三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法牛頓法解非線性方程牛頓法解非線性方程原理：原理：將非線性方程線性化將非線性方程線性化 taylor 展開展開取取 x0 x*，將將 f (x)在在 x0 做一階做一階taylor

20、展開展開:20000)(! 2)()()()(xxfxxxfxfxf ， 在在 x0 和和 x* 之間。之間。將將 (x* x0)2 看成高階小量，則有：看成高階小量，則有：)*)()(*)(0000 xxxfxfxf )()(*000 xfxfxx 線性線性xyx*x0 x1)()(1kkkkxfxfxx 迭代公式迭代公式： 將非線性代數(shù)方程組將非線性代數(shù)方程組 (1-22)(1-22) 在待求量在待求量 的某一個(gè)初始估計(jì)值的某一個(gè)初始估計(jì)值 附附近，展開成泰勒級(jí)數(shù)并略去二階及以上近，展開成泰勒級(jí)數(shù)并略去二階及以上的高階項(xiàng)，得到線性化方程組的高階項(xiàng)，得到線性化方程組 (1-24)(1-24)

21、 稱為牛頓法的修正方程式。稱為牛頓法的修正方程式。 牛頓牛頓-拉夫遜法拉夫遜法0)(xfx)0(x(0)(0)(0)()()0f xfxx 由上式根據(jù)初值由上式根據(jù)初值 可求得第一次迭可求得第一次迭代的修正量代的修正量 (1-25)(1-25) 將將 和和 相加，得到變量的第一次相加，得到變量的第一次改進(jìn)值改進(jìn)值 。 牛頓牛頓-拉夫遜法拉夫遜法(0)(0)1(0)()()xfxf x )0(x)0(x)0(x)1(x 因此，應(yīng)用牛頓法求解的迭代格式為因此，應(yīng)用牛頓法求解的迭代格式為 (1-26)(1-26) (1-27) (1-27) 上兩式中：上兩式中： 是函數(shù)是函數(shù) 對(duì)于對(duì)于 的的一階偏導(dǎo)

22、數(shù)矩陣，即雅可比矩陣一階偏導(dǎo)數(shù)矩陣，即雅可比矩陣 ，為，為迭代次數(shù)。迭代次數(shù)。 牛頓法當(dāng)初值牛頓法當(dāng)初值 和方程的精確解足和方程的精確解足夠接近時(shí)，具有平方收斂特性。夠接近時(shí)，具有平方收斂特性。 牛頓牛頓-拉夫遜法拉夫遜法( )( )( )()()kkkfxxf x )()()1(kkkxxx( )fx)(xfxjk)0(x（二）牛頓潮流算法的修正方程式（二）牛頓潮流算法的修正方程式 將牛頓法用于求解電力系統(tǒng)潮流計(jì)算問將牛頓法用于求解電力系統(tǒng)潮流計(jì)算問題時(shí)，由于所采用的數(shù)學(xué)表達(dá)式以及復(fù)題時(shí)，由于所采用的數(shù)學(xué)表達(dá)式以及復(fù)電壓變量采用的坐標(biāo)形式的不同，可以電壓變量采用的坐標(biāo)形式的不同，可以形成牛頓

23、潮流算法的不同形式。形成牛頓潮流算法的不同形式。 以下討論用得最廣泛的以下討論用得最廣泛的 采用功率采用功率方程式模型，而電壓變量分別采用極坐方程式模型，而電壓變量分別采用極坐標(biāo)和直角坐標(biāo)的兩種形式。標(biāo)和直角坐標(biāo)的兩種形式。牛頓牛頓-拉夫遜法拉夫遜法)(xf 1 1 極坐標(biāo)形式極坐標(biāo)形式 令令 ，對(duì)每個(gè)，對(duì)每個(gè) 節(jié)點(diǎn)及節(jié)點(diǎn)及 節(jié)點(diǎn)，節(jié)點(diǎn)，根據(jù)式根據(jù)式(1-13)(1-13)，有，有 (1-28)(1-28) 對(duì)每個(gè)對(duì)每個(gè) 節(jié)點(diǎn)，根據(jù)式節(jié)點(diǎn)，根據(jù)式(1-14)(1-14)，有，有 (1-29)(1-29)牛頓牛頓-拉夫遜法拉夫遜法iiiuupqpv0)sincos(iijijijijijjiip

24、bguuppq0)cossin(iijijijijijjiiqbguuq 將上述方程式在某個(gè)近似解附近用泰將上述方程式在某個(gè)近似解附近用泰勒級(jí)數(shù)展開，略去二階及以上的高階項(xiàng)勒級(jí)數(shù)展開，略去二階及以上的高階項(xiàng)后，得到以矩陣形式表示的修正方程式后，得到以矩陣形式表示的修正方程式 (1-30)(1-30) 式中：式中： 為節(jié)點(diǎn)個(gè)數(shù)，為節(jié)點(diǎn)個(gè)數(shù)， 為為 節(jié)點(diǎn)數(shù)，節(jié)點(diǎn)數(shù)，雅可比矩陣是雅可比矩陣是 階非奇異方陣。階非奇異方陣。牛頓牛頓-拉夫遜法拉夫遜法1111pnnhnmlqu unmnm nmpv22mn 雅可比矩陣各元素的表示式如下：雅可比矩陣各元素的表示式如下：牛頓牛頓-拉夫遜法拉夫遜法2(sinc

25、os) () (1-31) () (1-32) ijijijijjiiijiiiijuu gbj iphu bqj i2(cossin ) (j i) (1-33) (j i) (1-34)ijijijijijiijjiiiijuu gbpnuugpu2(cossin ) (j i) (1-35) (j i) (1-36)ijijijijijiijiiiijuu gbqmu gp2(sincos ) (j i) (1-37) (j i) (1-38)ijijijijijiijjiiiijuu gbqluu bqu 2 2 直角坐標(biāo)形式直角坐標(biāo)形式 令令 ，此時(shí)每個(gè)節(jié)點(diǎn)，都有兩，此時(shí)每個(gè)節(jié)點(diǎn)，都

26、有兩個(gè)方程式。因此共有個(gè)方程式。因此共有 個(gè)方程式。個(gè)方程式。 對(duì)每個(gè)對(duì)每個(gè)pqpq 節(jié)點(diǎn)，根據(jù)式節(jié)點(diǎn)，根據(jù)式(1-11)(1-11)和式和式(1-12)(1-12)有：有： (1-39)(1-39) (1-40) (1-40)牛頓牛頓-拉夫遜法拉夫遜法iiijfeu2(1)npq()()0iiijjijjiijjijjij ipe g eb ff g fb ep ()()0iiijjijjiijjij iij iqf g eb fe g fb eq 對(duì)每個(gè)對(duì)每個(gè) 節(jié)點(diǎn)，除了有與式節(jié)點(diǎn)，除了有與式(1-39)(1-39)相同的有功功率方程式之外，還有相同的有功功率方程式之外，還有 (1-41)

27、(1-41) 采用直角坐標(biāo)形式的修正方程式為采用直角坐標(biāo)形式的修正方程式為 (1-42)(1-42)牛頓牛頓-拉夫遜法拉夫遜法pv2222()0iiiiuefu 21e11 f1nphnnnmqmlnrsmu 雅可比矩陣各元素的表示式如下：雅可比矩陣各元素的表示式如下：牛頓牛頓-拉夫遜法拉夫遜法j i() (ji) (1-43)() (ji) (1-44)ijiijiiijijjijjii iiiijg eb fphg eb fg eb fe (ji) (1-45)() (ji) (1-46)ijiijiiijijjijjii iiiijj ib eg fpng fb eb eg ff (ji

28、) (1-47)() (ji) (1-48)ijiijiiijijjijjii iiiijj ib eg fqmg fb eb eg fej i (ji) (1-49)() (ji) (1-50)ijiijiiijijjijjii iiiijg eb fqlg eb fg eb ff牛頓牛頓-拉夫遜法拉夫遜法20 (ji) (1-51)2 (ji) (1-52)iijijuree20 (ji) (1-53)2 () (1-54)iijijusfjif 分析以上兩種類型的修正方程式，分析以上兩種類型的修正方程式，可以看出兩者具有以下的共同特點(diǎn)?？梢钥闯鰞烧呔哂幸韵碌墓餐攸c(diǎn)。 (1)(1) 修正

29、方程式的數(shù)目分別為修正方程式的數(shù)目分別為 及及 個(gè)，在個(gè)，在 節(jié)點(diǎn)比例不大時(shí)，兩節(jié)點(diǎn)比例不大時(shí)，兩者的方程式數(shù)目基本接近者的方程式數(shù)目基本接近 個(gè)。個(gè)。 (2)(2) 雅可比矩陣的元素都是節(jié)點(diǎn)電壓雅可比矩陣的元素都是節(jié)點(diǎn)電壓的函數(shù)，每次迭代，雅可比矩陣都需要的函數(shù)，每次迭代，雅可比矩陣都需要重新形成。重新形成。牛頓牛頓-拉夫遜法拉夫遜法mn1212npv12n (3)(3) 從雅可比陣非對(duì)角元素的表示式從雅可比陣非對(duì)角元素的表示式可見，某個(gè)非對(duì)角元素是否為零決定于可見，某個(gè)非對(duì)角元素是否為零決定于相應(yīng)的節(jié)點(diǎn)導(dǎo)納矩陣元素相應(yīng)的節(jié)點(diǎn)導(dǎo)納矩陣元素 是否為零。是否為零。如將修正方程式按節(jié)點(diǎn)號(hào)的次序排列

30、，如將修正方程式按節(jié)點(diǎn)號(hào)的次序排列，并將雅可比矩陣分塊，把每個(gè)并將雅可比矩陣分塊，把每個(gè) 階子階子陣陣 作為分塊矩陣的作為分塊矩陣的 元素，則按節(jié)點(diǎn)號(hào)順序而構(gòu)成的分塊雅元素，則按節(jié)點(diǎn)號(hào)順序而構(gòu)成的分塊雅可比矩陣將和節(jié)點(diǎn)導(dǎo)納矩陣具有同樣的可比矩陣將和節(jié)點(diǎn)導(dǎo)納矩陣具有同樣的稀疏結(jié)構(gòu)，是一個(gè)高度稀疏的矩陣。稀疏結(jié)構(gòu)，是一個(gè)高度稀疏的矩陣。牛頓牛頓-拉夫遜法拉夫遜法ijy22等如ijijijijijijijijsrnhlmnh (4)(4) 和節(jié)點(diǎn)導(dǎo)納矩陣具有相同稀疏結(jié)和節(jié)點(diǎn)導(dǎo)納矩陣具有相同稀疏結(jié)構(gòu)的分塊雅可比矩陣在位置上對(duì)稱，但構(gòu)的分塊雅可比矩陣在位置上對(duì)稱，但由于由于 ，所以雅可比矩陣所以雅可比矩

31、陣不是對(duì)稱陣不是對(duì)稱陣。 修正方程式的這些特點(diǎn)決定了牛頓法修正方程式的這些特點(diǎn)決定了牛頓法潮流程序特點(diǎn)，在設(shè)計(jì)算法時(shí)應(yīng)重點(diǎn)考潮流程序特點(diǎn)，在設(shè)計(jì)算法時(shí)應(yīng)重點(diǎn)考慮。慮。 牛頓牛頓-拉夫遜法拉夫遜法jiijjiijjiijjiijllmmnnhh,（三）修正方程式的處理和求解（三）修正方程式的處理和求解 有效地處理修正方程式是提高牛頓法潮流程有效地處理修正方程式是提高牛頓法潮流程序計(jì)算速度并降低內(nèi)存需求量的關(guān)鍵。序計(jì)算速度并降低內(nèi)存需求量的關(guān)鍵。 結(jié)合修正方程式的求解，目前實(shí)用的牛頓結(jié)合修正方程式的求解，目前實(shí)用的牛頓法潮流程序的程序特點(diǎn)主要有以下三個(gè)方面，法潮流程序的程序特點(diǎn)主要有以下三個(gè)方面，

32、這些程序特點(diǎn)對(duì)牛頓法潮流程序性能的提高起這些程序特點(diǎn)對(duì)牛頓法潮流程序性能的提高起著決定性的作用。著決定性的作用。 1 1 對(duì)于稀疏矩陣，在計(jì)算機(jī)中只儲(chǔ)存其非對(duì)于稀疏矩陣，在計(jì)算機(jī)中只儲(chǔ)存其非零元素，且只有非零元素才參加運(yùn)算。零元素，且只有非零元素才參加運(yùn)算。 牛頓牛頓-拉夫遜法拉夫遜法 牛頓牛頓-拉夫遜法拉夫遜法3 3 節(jié)點(diǎn)編號(hào)優(yōu)化節(jié)點(diǎn)編號(hào)優(yōu)化。經(jīng)過消元運(yùn)算得到。經(jīng)過消元運(yùn)算得到的上三角矩陣一般仍為稀疏矩陣，但由的上三角矩陣一般仍為稀疏矩陣，但由于消元過程中有新的非零元素注入，使于消元過程中有新的非零元素注入，使得它的稀疏度比原雅可比矩陣有所降低。得它的稀疏度比原雅可比矩陣有所降低。分析表明，

33、新增非零元素的多少和消元分析表明，新增非零元素的多少和消元的順序或節(jié)點(diǎn)編號(hào)有關(guān)。的順序或節(jié)點(diǎn)編號(hào)有關(guān)。牛頓牛頓-拉夫遜法拉夫遜法 節(jié)點(diǎn)編號(hào)優(yōu)化的作用即在于找到一種節(jié)點(diǎn)編號(hào)優(yōu)化的作用即在于找到一種網(wǎng)絡(luò)節(jié)點(diǎn)的重新編號(hào)方案，使得按此構(gòu)網(wǎng)絡(luò)節(jié)點(diǎn)的重新編號(hào)方案，使得按此構(gòu)成的節(jié)點(diǎn)導(dǎo)納矩陣以及和它相應(yīng)的雅可成的節(jié)點(diǎn)導(dǎo)納矩陣以及和它相應(yīng)的雅可比矩陣在高斯消元或三角分解過程中新比矩陣在高斯消元或三角分解過程中新增的非零元素?cái)?shù)目能盡量減少。增的非零元素?cái)?shù)目能盡量減少。 牛頓牛頓-拉夫遜法拉夫遜法 節(jié)點(diǎn)編號(hào)優(yōu)化通常有三種方法：節(jié)點(diǎn)編號(hào)優(yōu)化通常有三種方法：(1) (1) 靜態(tài)法靜態(tài)法按各節(jié)點(diǎn)靜態(tài)連接支路數(shù)按各節(jié)點(diǎn)靜

34、態(tài)連接支路數(shù)的多少順序編號(hào)。由少到多編號(hào)；的多少順序編號(hào)。由少到多編號(hào)；(2) (2) 半動(dòng)態(tài)法半動(dòng)態(tài)法一按各節(jié)點(diǎn)動(dòng)態(tài)連接支路一按各節(jié)點(diǎn)動(dòng)態(tài)連接支路數(shù)的多少順序編號(hào)；數(shù)的多少順序編號(hào)；(3) (3) 動(dòng)態(tài)法動(dòng)態(tài)法一按各節(jié)點(diǎn)動(dòng)態(tài)增加支路數(shù)一按各節(jié)點(diǎn)動(dòng)態(tài)增加支路數(shù)的多少順序編號(hào)。的多少順序編號(hào)。 消去節(jié)點(diǎn)后出現(xiàn)新支路數(shù)最少的節(jié)點(diǎn)。消去節(jié)點(diǎn)后出現(xiàn)新支路數(shù)最少的節(jié)點(diǎn)。牛頓牛頓-拉夫遜法拉夫遜法 三種節(jié)點(diǎn)編號(hào)優(yōu)化方法中動(dòng)態(tài)法效果三種節(jié)點(diǎn)編號(hào)優(yōu)化方法中動(dòng)態(tài)法效果最好，但優(yōu)化本身所需計(jì)算量也最多，最好，但優(yōu)化本身所需計(jì)算量也最多，而靜態(tài)法則反之。對(duì)于牛頓法潮流計(jì)算而靜態(tài)法則反之。對(duì)于牛頓法潮流計(jì)算來說，一般

35、認(rèn)為，采用半動(dòng)態(tài)法似乎是來說，一般認(rèn)為，采用半動(dòng)態(tài)法似乎是較好的選擇。較好的選擇。牛頓牛頓-拉夫遜法拉夫遜法牛頓牛頓-拉夫遜法拉夫遜法牛頓牛頓-拉夫遜法拉夫遜法牛頓牛頓-拉夫遜法拉夫遜法 三潮流計(jì)算的幾種基本方法三潮流計(jì)算的幾種基本方法 p-q分解法分解法p-q分解法原理分解法原理p-q分解法分解法rx p-q分解法分解法hp)/(uulqmn121n1mnp-q分解法分解法2010ijijbgijijijijbgsin1cos；2/iiuqiibiiiibuq2 p-q分解法分解法ijjiijbuuhijjiijbuulhub ulub u/bb 及1n1mnp-q分解法分解法p-q分解法分

36、解法p-qp-q分解法的修正方程式分解法的修正方程式p-q分解法分解法/()p ub u /q ubu /bb 及p-q分解法分解法/bbp-q分解法分解法b p-q分解法分解法/p ub/q ububbp-q分解法分解法222211;ijijiiijijj ij ij ij iijijijijiiioiij iijijijijj ibxbbxxxbbbbbrxrx ijijijbbb ijijxrp-q分解法分解法1n1mnmn12圖圖1-3 1-3 牛頓法和牛頓法和p-qp-q分解法的典型收斂特性分解法的典型收斂特性nrnr牛頓法；牛頓法；fdlffdlfp-qp-q分解法分解法 p-q分

37、解法分解法 p-q分解法分解法bbp-q分解法分解法jbbp-q分解法分解法p-q分解法分解法xr xrp-q分解法分解法xrxrxr 解決這個(gè)問題的途徑主要有以下兩種。解決這個(gè)問題的途徑主要有以下兩種。 對(duì)大對(duì)大 比值支路的參數(shù)加以補(bǔ)償比值支路的參數(shù)加以補(bǔ)償 1 1 對(duì)大對(duì)大 比值支路的參數(shù)加以補(bǔ)償比值支路的參數(shù)加以補(bǔ)償 對(duì)大對(duì)大 比值支路的參數(shù)加以補(bǔ)償，又比值支路的參數(shù)加以補(bǔ)償，又分成分成串聯(lián)補(bǔ)償法串聯(lián)補(bǔ)償法及及并聯(lián)補(bǔ)償法并聯(lián)補(bǔ)償法兩種。兩種。p-q分解法分解法xrxrxr (1) (1) 串聯(lián)補(bǔ)償法串聯(lián)補(bǔ)償法 這種方法的原理這種方法的原理見見圖圖1-61-6，其中，其中 為增為增加的虛構(gòu)

38、節(jié)點(diǎn)，加的虛構(gòu)節(jié)點(diǎn)， 為新增的補(bǔ)償電容。為新增的補(bǔ)償電容。 數(shù)值的選擇應(yīng)滿足數(shù)值的選擇應(yīng)滿足 支路支路 的條件。的條件。p-q分解法分解法mcjxcxmi rxxc)( p-q分解法分解法 這種方法的缺點(diǎn)是如果原來支路的這種方法的缺點(diǎn)是如果原來支路的 比值非常大，從而使比值非常大，從而使 的值選得過大時(shí)，的值選得過大時(shí)，新增節(jié)點(diǎn)新增節(jié)點(diǎn) 的電壓值有可能偏離節(jié)點(diǎn)的電壓值有可能偏離節(jié)點(diǎn) 及及 的電壓很多，這種不正常的電壓將的電壓很多，這種不正常的電壓將導(dǎo)致潮流計(jì)算收斂緩慢，甚至不收斂。導(dǎo)致潮流計(jì)算收斂緩慢，甚至不收斂。p-q分解法分解法xrcxmij p-q分解法分解法jbgjbjbbbjgyff

39、fij)2121(1)(ji ji mmufbji 圖圖1-7 1-7 對(duì)大對(duì)大r/xr/x比值支路的井聯(lián)補(bǔ)償比值支路的井聯(lián)補(bǔ)償(a) (a) 原支路；原支路；(b) (b) 補(bǔ)償后支路補(bǔ)償后支路 p-q分解法分解法p-q分解法分解法xrxrp-q分解法分解法p-q分解法分解法bxbbxbbxbbxxxrbbbx 方案采用的是嚴(yán)格的方案采用的是嚴(yán)格的 , , 交替迭代方案，這也是該算法和現(xiàn)在通交替迭代方案，這也是該算法和現(xiàn)在通行的行的 方案的標(biāo)準(zhǔn)型方案的標(biāo)準(zhǔn)型p-qp-q分解法的第二分解法的第二個(gè)差別。新算法若仍采用老的迭代方案，個(gè)差別。新算法若仍采用老的迭代方案，將會(huì)出現(xiàn)周期性的使功率偏差不

40、再下降將會(huì)出現(xiàn)周期性的使功率偏差不再下降的的 , , ， ， 循環(huán)迭代過程。循環(huán)迭代過程。p-q分解法分解法bxpuq xbpqqpqqpqqp-q分解法分解法xrxrxr四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法 四保留非線性潮流算法四保留非線性潮流算法 四保留非線性潮流算法四保留非線性潮流算法222()()iijijijijijijijijj iiijijijijijijijijj iiiipg eeb e fg f fb f eqg f eb f fg e fb eeuef 四保留非線性潮流算法四保留非線性潮流算法 為推導(dǎo)方便，將上述潮流方程寫成更普

41、為推導(dǎo)方便，將上述潮流方程寫成更普遍的齊次二次方程的形式。這里先定義：遍的齊次二次方程的形式。這里先定義： n n維未知變量向量維未知變量向量 n n維函數(shù)向量維函數(shù)向量 n n維函數(shù)給定值向量維函數(shù)給定值向量 四保留非線性潮流算法四保留非線性潮流算法tnxxxx,21tnxyxyxyxy)(,),(),()(21tsnsssyyyy,21 一個(gè)具有一個(gè)具有n n個(gè)變量的齊次二次代數(shù)方個(gè)變量的齊次二次代數(shù)方程式的普遍形式為程式的普遍形式為 四保留非線性潮流算法四保留非線性潮流算法111 1121 211212 1222 2221122( ) ()()() ()()() ()()()iiin i

42、niin inninninnn in ny xaxxaxxaxxax xax xax xax xax xax x(1-69)(1-69)于是潮流方程組可以寫成如下的矩陣形式：于是潮流方程組可以寫成如下的矩陣形式： (1-70)(1-70)或或 (1-71)(1-71)四保留非線性潮流算法四保留非線性潮流算法12( )snx xx xyy xax x( )( )0sf xy xy 式式(1-70)(1-70)中，系數(shù)矩陣為：中，系數(shù)矩陣為： (1-72)(1-72)四保留非線性潮流算法四保留非線性潮流算法11 112 11 121 122 12 11 12 1111 212 21 221 222

43、 2221 22 22111212122212( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )nnnnnnnnnnnnnnn nnnn nn nnnnn naaaaaaaaaaaaaaaaaaaaaaaaaaaa（二）泰勒級(jí)數(shù)展開式（二）泰勒級(jí)數(shù)展開式 對(duì)式對(duì)式(1-69)(1-69)在初值附近展開，可得到在初值附近展開，可得到?jīng)]沒有截?cái)嗾`差有截?cái)嗾`差的精確展開式為：的精確展開式為： (1-73)(1-73)即即 (1-74)(1-74) 四保留非線性潮流算法四保留非線性潮流算法2

44、(0)(0)(0)1111()2nnnsiijjkjjkjjkyyyy xxxxxxxxxxx ！12(0)1()2snxxxxyy xj xhxx 式中：式中： 為修正量向量。為修正量向量。 為雅可比矩陣。為雅可比矩陣。四保留非線性潮流算法四保留非線性潮流算法tnxxxxxx,21)0(1111222212(0)12 nn nnnnnnyyyxxxyyyxxxjjrxxyyyxxx 是一個(gè)常數(shù)矩陣，其階數(shù)很高，但高度是一個(gè)常數(shù)矩陣，其階數(shù)很高，但高度稀疏。稀疏。 四保留非線性潮流算法四保留非線性潮流算法nnnnnnnnnnnnnnnnnnnnnnnnnxxyxxyxxyxxyxxyxxyxx

45、yxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyxxyh222122222212212212112222221222222222122221222122112221221211221222121212112211211122nnrh 四保留非線性潮流算法四保留非線性潮流算法)()()0(xyxjxyysixiiixxx)0(0)(0)(0)(0)(0)(0)()() ijiijjijijjiijx xxxxxxxxxxxxx 四保留非線性潮流算法四保留非線性潮流算法)0(xnnjinnjinnjinnjisxxxxxxxxax

46、xxxxxxxaxxxxxxxxaxxxxxxxxay2111)0()0()0(21)0(11)0()0(2)0(11)0(1)0()0()0()0()0(2)0(1)0(1)0(1四保留非線性潮流算法四保留非線性潮流算法)()0(xy 四保留非線性潮流算法四保留非線性潮流算法( )y x)0(x)0(221122221121221111xxnnnnnnnnnxxyxxyxxyxxyxxyxxyxxyxxyxxyxj四保留非線性潮流算法四保留非線性潮流算法jiijjjjiiijijixxxxxxxxxxxxxxd)()()( 四保留非線性潮流算法四保留非線性潮流算法)0(x)0(2)0()0(

47、)0()0(212)0(1)0(111)0(1xxxxxxxxxxxxxxxxannnjiji 四保留非線性潮流算法四保留非線性潮流算法)( xy x)0(x四保留非線性潮流算法四保留非線性潮流算法x)0(xx)()()0(1xyyxyjxs)()()()0(1) 1(kskxyyxyjx式中：為迭代次數(shù)；式中：為迭代次數(shù)； 按按 估計(jì)而得。估計(jì)而得。 進(jìn)行第一次迭代時(shí)，進(jìn)行第一次迭代時(shí)， ，令，令 , ,同牛頓法的第一次迭代計(jì)算完全相同。同牛頓法的第一次迭代計(jì)算完全相同。 算法的收斂判據(jù)為算法的收斂判據(jù)為 (1-84)(1-84)也可以采用也可以采用 (1-85)(1-85)作為收斂判據(jù)作為

48、收斂判據(jù)。 式式(1-85)(1-85)是比式是比式(1-84)(1-84)更合理的收斂更合理的收斂判據(jù)。判據(jù)。四保留非線性潮流算法四保留非線性潮流算法kj)0(xx 0k0)()0(xy)()1(maxkikiixx)()1(maxkiikiiixyxy四保留非線性潮流算法四保留非線性潮流算法啟動(dòng)輸入原始數(shù)據(jù)賦初值形成雅可比矩陣j形成j 因子表計(jì)算二階項(xiàng)形成節(jié)點(diǎn)導(dǎo)納矩陣()()kyx是k=0(0)0 x用式(1-83)求解(1)kx(1)()max?kkiiixx(1)(0)(1)kkxxx計(jì)算支路潮流輸出結(jié)果停機(jī)k=k+1否四保留非線性潮流算法四保留非線性潮流算法( )( )0sf xy

49、xy)()() 1()(1)()()()(kkkskkkxxxyxyxjx)1()0()1()()0(1)0()1()()()(kkkskxxxxyyxyxjx四保留非線性潮流算法四保留非線性潮流算法 四保留非線性潮流算法四保留非線性潮流算法 0 x四保留非線性潮流算法四保留非線性潮流算法 x)(kx)(kx)(kx)0(x四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法 xf四保留非線性潮流算法四保留非線性潮流算法 0 xf(0)(0)()()()0f xfxxhx )( xh (0)(1)(0)( )()()()kkfxxf xhx四保留非線性潮流算法四保留

50、非線性潮流算法0k0)0(x0)()0(xh(0)(1)(0)()()fxxf x四保留非線性潮流算法四保留非線性潮流算法1k(0)(2)(0)(1)()()()fxxf xhx(0)(1)(0)(0)(1)(1)()()()()f xxf xf xxhx)()()1()0()1(xxfxh(0)(2)(0)(0)(1)()()()fxxf xf xx四保留非線性潮流算法四保留非線性潮流算法 k(0)( )(0)(1)()()()kkfxxf xhx11)()0()1()()(kiikxxfxh, 1 , 0)()()(1)() 0 () 0 () 1() 0 (kxxfxfxxfkiik四保

51、留非線性潮流算法四保留非線性潮流算法 四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法 四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法xr四保留非線性潮流算法四保留非線性潮流算法()()iiijjijjiijjijjjipe g eb ff g fb e ()()iiijjijjiijjijijiqfgebfegfbe 222()iiiuef 四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法四保留非線性潮流算法iiijj ij iiiijj ij iggbb 四保留非線性潮流算

52、法四保留非線性潮流算法220() ()()iiiiiijjijjiijjijjj ipgefe g eb ff g fb e220()()()iiiiiijjijjiijjijjj iqbeff g eb fe g fb e00iibg 、iiiiibg 及四保留非線性潮流算法四保留非線性潮流算法/iiqp 及)(220/iiiiifegpp)(220/iiiiifebqq四保留非線性潮流算法四保留非線性潮流算法)()(feqfepii，、，ijijjijjiiispffpeepfepfep)()()0()0(，ijijjijjiiisqffqeeqfeqfeq),()()0()0(，iisq

53、sp 及四保留非線性潮流算法四保留非線性潮流算法0)0()0(jejfesii0),(),()0()0()0()0(feqfepii四保留非線性潮流算法四保留非線性潮流算法)0()0(iiqp及sqspeflmnhqp/四保留非線性潮流算法四保留非線性潮流算法ijsjiijiisiiiiijsjiijiisiiiiijsjiijiisiiiiijsjiijiisiiiibeeqlbeeqlgefqmgefqmgeepngeepnbefphbefph,四保留非線性潮流算法四保留非線性潮流算法sqspefbggbeqps/qpsqsprqrpssrp ebgfrq egbe四保留非線性潮流算法四保留非線性潮流算法effe 、ijjijjijijijjijiiijjijjijijijjijiiebfgefbegfsqebfgffbegesp)()()()(四保留非線性潮流算法四保留非線性潮流算法()()ijjijjiij iijjijjisj igebfrp egfberq eiiisiisiiisiisspe rp efrq esqfrp ee rq e 四保留非線性潮流算法四保留非線性潮流算法(1)( )( )( )( )(1)( )( )( )( )kkkkkiiisiiskkkkkiiisiissperpefr