兩機五節(jié)點網(wǎng)絡潮流計算牛拉法

1、畢 業(yè) 設 計 (論 文) 信息工程學院 系 電氣工程及其自動化 專業(yè)畢業(yè)設計（論文）題目 兩機五節(jié)點網(wǎng)絡潮流計算-牛拉法 學 生 姓 名 田玖婷 班 級 10電氣一班 學 號 指 導 教 師 基于matlab的電力系統(tǒng)潮流仿真計算the power flow simulation of power system based on matlab總計 畢業(yè)設計（論文） 頁 表 格 個插 圖 幅摘 要潮流計算是電力系統(tǒng)的一項重要分析功能，是進行故障計算，繼電保護整定，安全分析的必要工具。傳統(tǒng)的潮流計算程序缺乏圖形用戶界面，結果顯示不直觀，難于與其他分析功能集成。網(wǎng)絡原始數(shù)據(jù)輸入工作量大且易于出錯。

2、隨著計算機技術的飛速發(fā)展，microsoft windows操作系統(tǒng)早已被大家所熟悉，其友好的圖形用戶界面已成為pc機的標準，而dos操作系統(tǒng)下的應用程序因其界面不夠友好，開發(fā)具有windows風格界面的電力系統(tǒng)分析軟件已成為當前的主流趨勢。另外，傳統(tǒng)的程序設計方法是結構化程序設計方法，該方法基于功能分解，把整個軟件工程看作是一個個對象的組合，由于對某個特定問題域來說，該對象組成基本不變，因此，這種基于對象分解方法設計的軟件結構上比較穩(wěn)定，易于維護和擴充。本文介紹了圖形化潮流計算軟件的開發(fā)設計思想和總體結構，闡述了該軟件所具備的功能和特點。結合電力系統(tǒng)的特點，軟件采用 matlab語言運行于w

3、indows操作系統(tǒng)的圖形化潮流計算軟件。本系統(tǒng)的主要特點是操作簡單，圖形界面直觀，運行穩(wěn)定.計算準確。計算中，算法做了一些改進，提高了計算速度，各個類的有效封裝又使程序具有很好的模塊性.可維護性和可重用性。關鍵詞：電力系統(tǒng)潮流計算；牛頓拉夫遜法潮流計算； matlab目錄摘要第一章 電力系統(tǒng)潮流計算概述 1.1電力系統(tǒng)簡介1.2潮流計算簡介1.3潮流計算的意義及其發(fā)展 第二章 潮流計算的數(shù)學模型2.1導納矩陣的原理及計算方法 2.1.1自導納和互導納的確定方法2.1.2節(jié)點導納矩陣的性質及意義 2.1.3非標準變比變壓器等值電路2.2潮流計算的基本方程 2.3電力系統(tǒng)節(jié)點分類2.4潮流計算的

4、約束條件(1第三章 牛頓拉夫遜法概述3.1牛頓-拉夫遜法基本原理 3.3牛頓-拉夫遜法求解過程3.2牛頓-拉夫遜法程序框圖 第四章 matlab概述4.1matlab簡介 4.2矩陣的生成4.3矩陣的運算4.4牛頓拉夫遜法潮流計算程序 總結參考文獻 (34)第一章 電力系統(tǒng)潮流計算概述 1.1 電力系統(tǒng)敘述電力工業(yè)發(fā)展初期，電能是直接在用戶附近的發(fā)電站(或稱發(fā)電廠)中生產(chǎn)的，各發(fā)電站孤立運行。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源(如煤田)和水能資源豐富的地區(qū)又往往遠離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠距離

5、輸送給電力用戶。同時，為了提高供電可靠性以及資源利用的綜合經(jīng)濟性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。這種由發(fā)電機、升壓和降壓變電所，送電線路以及用電設備有機連接起來的整體，即稱為電力系統(tǒng)。電力系統(tǒng)加上發(fā)電機的原動機(如汽輪機、水輪機)，原動機的力能部分(如熱力鍋爐、水庫、原子能電站的反應堆)、供熱和用熱設備，則稱為動力系統(tǒng)?，F(xiàn)代電力系統(tǒng)提出了“靈活交流輸電與新型直流輸電”的概念。靈活交流輸電技術是指運用固態(tài)電子器件與現(xiàn)代自動控制技術對交流電網(wǎng)的電壓、相位角、阻抗、功率以及電路的通斷進行實時閉環(huán)控制，從而提高高壓輸電線路的輸送能力和電力系統(tǒng)的穩(wěn)定水平。新型直流輸電技術

6、是指應用現(xiàn)電力電子技術的最新成果，改善和簡化變流站的造價等。運行方式管理中，潮流是確定電網(wǎng)運行方式的基本出發(fā)點；在規(guī)劃領域，需要進行潮流分析驗證規(guī)劃方案的合理性；在實時運行環(huán)境，調度員潮流提供了電網(wǎng)在預想操作情況下電網(wǎng)的潮流分布以校驗運行可靠性。在電力系統(tǒng)調度運行的多個領域都涉及到電網(wǎng)潮流計算。潮流是確定電力網(wǎng)絡運行狀態(tài)的基本因素，潮流問題是研究電力系統(tǒng)穩(wěn)態(tài)問題的基礎和前提。待設計電氣設備系統(tǒng)圖系統(tǒng)接線圖（其中節(jié)點1為平衡節(jié)點，節(jié)點2、3、4、5為pq節(jié)點。）1.2潮流計算簡介電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀

7、態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設計和現(xiàn)有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性?？煽啃院徒?jīng)濟性。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎。所以潮流計算是研究電力系統(tǒng)的一種很重要和基礎的計算。 電力系統(tǒng)潮流計算也分為離線計算和在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經(jīng)常監(jiān)視及實時控制。 利用電子數(shù)字計算機進行電力系統(tǒng)潮流計算從50年代中期就已經(jīng)開始。在這20年內(nèi)，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要圍繞著對潮流計算的一些基本要求進

8、行的。對潮流計算的要求可以歸納為下面幾點：（1）計算方法的可靠性或收斂性；（2）對計算機內(nèi)存量的要求；（3）計算速度；（4）計算的方便性和靈活性。電力系統(tǒng)潮流計算問題在數(shù)學上是一組多元非線性方程式求解問題，其解法都離不開迭代。因此，對潮流計算方法，首先要求它能可靠地收斂，并給出正確答案。由于電力系統(tǒng)結構及參數(shù)的一些特點，并且隨著電力系統(tǒng)不斷擴大，潮流計算的方程式階數(shù)也越來越高，對這樣的方程式并不是任何數(shù)學方法都能保證給出正確答案的。這種情況成為促使電力系統(tǒng)計算人員不斷尋求新的更可靠方法的重要因素。1.3 潮流計算的意義及其發(fā)展電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)

9、正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。實際電力系統(tǒng)的潮流技術那主要采用牛頓-拉夫遜法。在運行方式管理中，潮流是確定電網(wǎng)運行方式的基本出發(fā)點；在規(guī)劃領域，需要進行潮流分析驗證規(guī)劃方案的合理性；在實時運行環(huán)境，調度員潮流提供了多個在預想操作

10、情況下電網(wǎng)的潮流分布以校驗運行可靠性。在電力系統(tǒng)調度運行的多個領域都涉及到電網(wǎng)潮流計算。潮流是確定電力網(wǎng)絡運行狀態(tài)的基本因素，潮流問題是研究電力系統(tǒng)穩(wěn)態(tài)問題的基礎和前提。在用數(shù)字計算機解電力系統(tǒng)潮流問題的開始階段，普遍采取以節(jié)點導納矩陣為基礎的逐次代入法。這個方法的原理比較簡單，要求的數(shù)字計算機內(nèi)存量比較下，適應50年代電子計算機制造水平和當時電力系統(tǒng)理論水平。但它的收斂性較差，當系統(tǒng)規(guī)模變大時，迭代次數(shù)急劇上升，在計算中往往出現(xiàn)迭代不收斂的情況。這就迫使電力系統(tǒng)計算人員轉向以阻抗矩陣為基礎的逐次代入法。阻抗法改善了系統(tǒng)潮流計算問題的收斂性，解決了導納法無法求解的一些系統(tǒng)的潮流計算，在60年代

11、獲得了廣泛的應用。阻抗法的主要缺點是占用計算機內(nèi)存大，每次迭代的計算量大。當系統(tǒng)不斷擴大時，這些缺點就更加突出。為了克服阻抗法在內(nèi)存和速度方面的缺點，60年代中期發(fā)展了以阻抗矩陣為基礎的分塊阻抗法。這個方法把一個大系統(tǒng)分割為幾個小的地區(qū)系統(tǒng)，在計算機內(nèi)只需要存儲各個地區(qū)系統(tǒng)的阻抗矩陣及它們之間聯(lián)絡線的阻抗，這樣不僅大幅度地節(jié)省了內(nèi)存容量，同時也提高了計算速度?？朔杩狗ㄈ秉c的另一途徑是采用牛頓-拉夫遜法。這是數(shù)學中解決非線性方程式的典型方法，有較好的收斂性。在解決電力系統(tǒng)潮流計算問題時，是以導納矩陣為基礎的，因此，只要我們能在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓法潮流

12、程序的效率。自從60年代中期，在牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。內(nèi)存要求。速度方面都超過了阻抗法，成為60年代末期以后廣泛采用的優(yōu)秀方法。第二章 潮流計算的數(shù)學模型2.1導納矩陣的原理及計算方法2.1.1自導納和互導納的確定方法電力網(wǎng)絡的節(jié)點電壓方程： (2-1)為節(jié)點注入電流列向量，注入電流有正有負，注入網(wǎng)絡的電流為正，流出網(wǎng)絡的電流為負。根據(jù)這一規(guī)定，電源節(jié)點的注入電流為正，負荷節(jié)點為負。既無電源又無負荷的聯(lián)絡節(jié)點為零，帶有地方負荷的電源節(jié)點為二者代數(shù)之和。為節(jié)點電壓列向量，由于節(jié)點電壓是對稱于參考節(jié)點而言的，因而需先選定參考節(jié)點。在電力系統(tǒng)中一般以地為參考節(jié)點。如整個網(wǎng)

13、絡無接地支路，則需要選定某一節(jié)點為參考。設網(wǎng)絡中節(jié)點數(shù)為（不含參考節(jié)點），則，均為n*n列向量。為n*n階節(jié)點導納矩陣。節(jié)電導納矩陣的節(jié)點電壓方程： 展開為： ： (2-2)是一個n*n階節(jié)點導納矩陣，其階數(shù)就等于網(wǎng)絡中除參考節(jié)點外的節(jié)點數(shù)。 節(jié)點導納矩陣的對角元素 (i=1，2，n)成為自導納。自導納數(shù)值上就等于在i節(jié)點施加單位電壓，其他節(jié)點全部接地時，經(jīng)節(jié)點i注入網(wǎng)絡的電流，因此，它可以定義為： (2-3)節(jié)點i的自導納數(shù)值上就等于與節(jié)點直接連接的所有支路導納的總和。節(jié)點導納矩陣的非對角元素 (j=1，2，n;i=1，2，。，n;j=i)稱互導納，由此可得互導納數(shù)值上就等于在節(jié)點i施加單位

14、電壓，其他節(jié)點全部接地時，經(jīng)節(jié)點j注入網(wǎng)絡的電流，因此可定義為： (2-4)節(jié)點j，i之間的互導納數(shù)值上就等于連接節(jié)點j，i支路到導納的負值。顯然，恒等于?；Ъ{的這些性質決定了節(jié)點導納矩陣是一個對稱稀疏矩陣。而且，由于每個節(jié)點所連接的支路數(shù)總有一個限度，隨著網(wǎng)絡中節(jié)點數(shù)的增加非零元素相對愈來愈少，節(jié)點導納矩陣的稀疏度，即零元素數(shù)與總元素的比值就愈來愈高。2.1.2節(jié)點導納矩陣的性質及意義節(jié)點導納矩陣的性質：（1）為對稱矩陣，=。如網(wǎng)絡中含有源元件，如移相變壓器，則對稱性不再成立。（2）對無接地支路的節(jié)點，其所在行列的元素之和均為零，即 。對于有接地支路的節(jié)點，其所在行列的元素之和等于該點接地

15、支路的導納。利用這一性質，可以檢驗所形成節(jié)點導納矩陣的正確性。（3）具有強對角性：對角元素的值不小于同一行或同一列中任一元素（4）為稀疏矩陣，因節(jié)點i ，j 之間無支路直接相連時=0，這種情況在實際電力系統(tǒng)中非常普遍。矩陣的稀疏性用稀疏度表示，其定義為矩陣中的零元素與全部元素之比，即 ， 式中z 為中的零元素。s 隨節(jié)點數(shù)n 的增加而增加：n=50，s可達92%；n=100，s 可達90%；n=500，s可達99%，充分利用節(jié)點導納矩陣的稀疏性可節(jié)省計算機內(nèi)存，加快計算速度，這種技巧稱為稀疏技術。節(jié)點導納矩陣的意義：是n*n階方陣，其對角元素 (i=1，2，-n)稱為自導納，非對角元素(i，j

16、=1，2，n， )稱為互導納。將節(jié)點電壓方程展開為可見， (2-5)表明，自導納在數(shù)值上等于僅在節(jié)點i施加單位電壓而其余節(jié)點電壓均為零（即其余節(jié)點全部接地）時，經(jīng)節(jié)點i注入網(wǎng)絡的電流。其顯然等于與節(jié)點i直接相連的所有支路的導納之和。同時可見。表明，互導納在數(shù)值上等于僅在節(jié)點j施加單位電壓而其余節(jié)點電壓均為零時，經(jīng)節(jié)點i注入網(wǎng)絡的電流，其顯然等于()即=。為支路的導納，負號表示該電流流出網(wǎng)絡。如節(jié)點ij之間無支路直接相連，則該電流為0，從而=0。注意字母幾種不寫法的不同意義：粗體黑字表示導納矩陣，大寫字母代矩陣中的第i行第j列元素，即節(jié)點i和節(jié)點j之間的互導納。小寫字母i，j支路的導納等于支路阻

17、抗的倒數(shù)數(shù)，。根據(jù)定義直接求取節(jié)點導納矩陣時，注意以下幾點：1)。節(jié)點導納矩陣是方陣，其階數(shù)就等于網(wǎng)絡中除去參考節(jié)點外的節(jié)點數(shù)。參考節(jié)點一般取大地，編號為零。2)。節(jié)點導納矩陣是稀疏矩陣，其各行非零非對角元素就等于與該行相對應節(jié)點所連接的不接地支路數(shù)。3)。節(jié)點導納矩陣的對角元素就等于各該節(jié)點所連接導納的總和。因此，與沒有接地支路的節(jié)點對應的行或列中，對角元素為非對角元素之和的負值。4)。節(jié)點導納矩陣的非對角元素等于連接節(jié)點i，j支路導納的負值。因此，一般情況下，節(jié)點導納矩陣的對角元素往往大于非對角元素的負值。5)。節(jié)點導納矩陣一般是對稱矩陣，這是網(wǎng)絡的互易特性所決定的。從而，一般只要求求取這

18、個矩陣的上三角或下三角部分。2.1.3非標準變比變壓器等值電路變壓器型等值電路更便于計算機反復計算,更適宜于復雜網(wǎng)絡的潮流計算.雙繞組變壓器可用阻抗與一個理想變壓器串聯(lián)的電路表示.理想變壓器只是一個參數(shù),那就是變比?，F(xiàn)在變壓器阻抗按實際變比歸算到低壓側為例,推導出變壓器型等值電路. a 雙繞組變壓器原理圖b 變壓器阻抗歸算到低壓側等值模型流入和流出理想變壓器的功率相等 (2-6)式中, 是理想變壓器的變比,和 分別為變壓器高,低繞組的實際電壓.從圖b直接可得： （2-7） 從而可得: （2-8）式中,又因節(jié)點電流方程應具有如下形式: （2-9）將式（1-8）與（1-9）比較，得： 因此可得各支

19、路導納為: （2-10） 由此可得用導納表示的變壓器型等值電路:圖 c2.2潮流計算的基本方程在潮流問題中，任何復雜的電力系統(tǒng)都可以歸納為以下元件（參數(shù)）組成。（1）發(fā)電機（注入電流或功率）（2）負荷（注入負的電流或功率）（3）輸電線支路（電阻，電抗）（4）變壓器支路（電阻，電抗，變比）（5）母線上的對地支路（阻抗和導納）（6）線路上的對地支路（一般為線路充電點容導納）集中了以上各類型的元件的簡單網(wǎng)絡如圖 (a) 潮流計算用的電網(wǎng)結構圖(b) 潮流計算等值網(wǎng)絡采用導納矩陣時，節(jié)點注入電流和節(jié)點電壓構成以下線性方程組 (2-11)其中 可展開如下形式 (2-12)由于實際電網(wǎng)中測量的節(jié)點注入量一

20、般不是電流而是功率，因此必須將式中的注入電流用節(jié)點注入功率來表示。節(jié)點功率與節(jié)點電流之間的關系為 (2-13)式中，因此用導納矩陣時，pq節(jié)點可以表示為把這個關系代入式中 ，得（2-14）式（3-4 ）就是電力系統(tǒng)潮流計算的數(shù)學模型-潮流方程。它具有如下特點：（1）它是一組代數(shù)方程，因而表征的是電力系統(tǒng)的穩(wěn)定運行特性。（2）它是一組非線性方程，因而只能用迭代方法求其數(shù)值解。（3）由于方程中的電壓和導納既可以表為直角坐標，又可表為極坐標，因而潮流方程有多種表達形式-極坐標形式，直角坐標形式和混合坐標形式。a。取 ，得到潮流方程的極坐標形式： (2-15)b。 取 ， ，得到潮流方程的直角坐標形式

21、： (2-16)c。取， ，得到潮流方程的混合坐標形式： (2-17)不同坐標形式的潮流方程適用于不同的迭代解法。例如：利用牛頓-拉夫遜迭代法求解，以直角坐標和混合坐標形式的潮流方程為方便；而p-q解耦法是在混合坐標形式的基礎上發(fā)展而成，故當然采用混合坐標形式。（4）它是一組n個復數(shù)方程，因而實數(shù)方程數(shù)為2n個但方程中共含4n個變量：p，q，u和，i=1，2，n，故必須先指定2n個變量才能求解。2.3電力系統(tǒng)節(jié)點分類用一般的電路理論求解網(wǎng)絡方程，目的是給出電壓源(或電流源)研究網(wǎng)絡內(nèi)的電流(或電壓)分布，作為基礎的方程式，一般用線性代數(shù)方程式表示。然而在電力系統(tǒng)中，給出發(fā)電機或負荷連接母線上電

22、壓或電流(都是向量)的情況是很少的，一般是給出發(fā)電機母線上發(fā)電機的有功功率(p)和母線電壓的幅值(u)，給出負荷母線上負荷消耗的有功功率(p)和無功功率(q)。主要目的是由這些已知量去求電力系統(tǒng)內(nèi)的各種電氣量。所以，根據(jù)電力系統(tǒng)中各節(jié)點性質的不同，很自然地把節(jié)點分成三類： pq節(jié)點對這一類點，事先給定的是節(jié)點功率(p，q)，待求的未知量是節(jié)點電壓向量(u，)，所以叫pq節(jié)點。通常變電所母線都是pq節(jié)點，當某些發(fā)電機的輸出功率p。q給定時，也作為pq節(jié)點。pq節(jié)點上的發(fā)電機稱之為pq機(或pq給定型發(fā)電機)。在潮流計算中，系統(tǒng)大部分節(jié)點屬于pq節(jié)點。 pu節(jié)點這類節(jié)點給出的參數(shù)是該節(jié)點的有功功率

23、p及電壓幅值u，待求量為該節(jié)點的無功功率q及電壓向量的相角。這類節(jié)點在運行中往往要有一定可調節(jié)的無功電源。用以維持給定的電壓值。通常選擇有一定無功功率儲備的發(fā)電機母線或者變電所有無功補償設備的母線做pu節(jié)點處理。pu節(jié)點上的發(fā)電機稱為pu機(或pu給定型發(fā)電機) 平衡節(jié)點在潮流計算中，這類節(jié)點一般只設一個。對該節(jié)點，給定其電壓值，并在計算中取該節(jié)點電壓向量的方向作為參考軸，相當于給定該點電壓向量的角度為零。也就是說，對平衡節(jié)點給定的運行參數(shù)是u和，因此有城為u節(jié)點，而待求量是該節(jié)點的p。q，整個系統(tǒng)的功率平衡由這一節(jié)點承擔。關于平衡節(jié)點的選擇，一般選擇系統(tǒng)中擔任調頻調壓的某一發(fā)電廠(或發(fā)電機)

24、，有時也可能按其他原則選擇，例如，為提高計算的收斂性?？梢赃x擇出線數(shù)多或者靠近電網(wǎng)中心的發(fā)電廠母線作平衡節(jié)點。以上三類節(jié)點4個運行參數(shù)p。q。u。中，已知量都是兩個，待求量也是兩個，只是類型不同而已。2.4潮流計算的約束條件 電力系統(tǒng)運行必須滿足一定技術和經(jīng)濟上的要求。這些要求夠成了潮流問題中某些變量的約束條件，常用的約束條件如下：1. 節(jié)點電壓應滿足2. (2-18)從保證電能質量和供電安全的要求來看，電力系統(tǒng)的所有電氣設備都必須運行在額定電壓附近。pu節(jié)點電壓幅值必須按上述條件給定。因此，這一約束條件對pq節(jié)點而言。3. 節(jié)點的有功功率和無功功率應滿足4. (2-19) pq節(jié)點的有功功率

25、和無功功率，以及pu節(jié)點的有功功率，在給定是就必須滿足上述條件，因此，對平衡節(jié)點的p和q以及pu節(jié)點的q應按上述條件進行檢驗。5. 節(jié)點之間電壓的相位差應滿足 (2-30) 為了保證系統(tǒng)運行的穩(wěn)定性，要求某些輸電線路兩端的電壓相位不超過一定的數(shù)值。這一約束的主要意義就在于此。 因此，潮流計算可以歸結為求解一組非線性方程組，并使其解答滿足一定的約束條件。常用的方法是迭代法和牛頓法，在計算過程中，或得出結果之后用約束條件進行檢驗。如果不能滿足要求，則應修改某些變量的給定值，甚至修改系統(tǒng)的運行方式，重新進行計算。第三章 牛頓拉夫遜法概述3.1牛頓-拉夫遜法基本原理電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一

26、種最基本的計算，是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。實際電力系統(tǒng)的潮流技術那主要采用牛頓-拉夫遜法。牛頓-拉夫遜法(簡稱牛頓法)在數(shù)學上是求解非線性代數(shù)方程式的有效方法。其要點是把非線性方程式的求解過程變成反復地對相應

27、的線性方程式進行求解的過程。即通常所稱的逐次線性化過程。對于非線性代數(shù)方程組： 即 (3-1-1)在待求量x的某一個初始估計值附近，將上式展開成泰勒級數(shù)并略去二階及以上的高階項，得到如下的經(jīng)線性化的方程組： (3-1-2)上式稱之為牛頓法的修正方程式。由此可以求得第一次迭代的修正量 (3-1-3)將和相加，得到變量的第一次改進值。接著就從出發(fā)，重復上述計算過程。因此從一定的初值出發(fā)，應用牛頓法求解的迭代格式為： (3-1-4) (3-1-5)上兩式中：是函數(shù)對于變量x的一階偏導數(shù)矩陣，即雅可比矩陣j;k為迭代次數(shù)。有上式可見，牛頓法的核心便是反復形式并求解修正方程式。牛頓法當初始估計值和方程的

28、精確解足夠接近時，收斂速度非常快，具有平方收斂特性。牛頓潮流算法突出的優(yōu)點是收斂速度快，若選擇到一個較好的初值，算法將具有平方收斂特性，一般迭代45次便可以收斂到一個非常精確的解。而且其迭代次數(shù)與所計算網(wǎng)絡的規(guī)?；緹o關。牛頓法也具有良好的收斂可靠性，對于對以節(jié)點導納矩陣為基礎的高斯法呈病態(tài)的系統(tǒng)，牛頓法也能可靠收斂。牛頓法所需的內(nèi)存量及每次迭代所需時間均較高斯法多。牛頓法的可靠收斂取決于有一個良好的啟動初值。如果初值選擇不當，算法有可能根本不收斂或收斂到一個無法運行的節(jié)點上。對于正常運行的系統(tǒng)，各節(jié)點電壓一般均在額定值附近，偏移不會太大，并且各節(jié)點間的相位角差也不大，所以對各節(jié)點可以采用統(tǒng)一

29、的電壓初值(也稱為平直電壓)，如假定： 或 (3-1-6) 這樣一般能得到滿意的結果。但若系統(tǒng)因無功緊張或其它原因導致電壓質量很差或有重載線路而節(jié)點間角差很大時，仍用上述初始電壓就有可能出現(xiàn)問題。解決這個問題的辦法可以用高斯法迭代12次，以此迭代結果作為牛頓法的初值。也可以先用直流法潮流求解一次以求得一個較好的角度初值，然后轉入牛頓法迭代。3.2牛頓-拉夫遜法潮流求解過程以下討論的是用直角坐標形式的牛頓拉夫遜法潮流的求解過程。當采用直角坐標時，潮流問題的待求量為各節(jié)點電壓的實部和虛部兩個分量由于平衡節(jié)點的電壓向量是給定的，因此待求兩共需要2(n-1)個方程式。事實上，除了平衡節(jié)點的功率方程式在

30、迭代過程中沒有約束作用以外，其余每個節(jié)點都可以列出兩個方程式。對pq節(jié)點來說，是給定的，因而可以寫出 （3-2-1）對pv節(jié)點來說，給定量是，因此可以列出 （3-2-2）求解過程大致可以分為以下步驟：（1）形成節(jié)點導納矩陣（2）將各節(jié)點電壓設初值u，（3）將節(jié)點初值代入相關求式，求出修正方程式的常數(shù)項向量（4）將節(jié)點電壓初值代入求式，求出雅可比矩陣元素（5）求解修正方程，求修正向量（6）求取節(jié)點電壓的新值（7）檢查是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉入下一步（8）計算支路功率分布，pv節(jié)點無功功率和平衡節(jié)點柱入功率。 以直角坐標系形式表示. 迭

31、代推算式 采用直角坐標時,節(jié)點電壓相量及復數(shù)導納可表示為: (3-2-3)將以上二關系式代入上式中,展開并分開實部和虛部;假定系統(tǒng)中的第1,2,m號為pq節(jié)點,第m+1,m+2,n-1為pv節(jié)點,根據(jù)節(jié)點性質的不同,得到如下迭代推算式: 對于pq節(jié)點 (3-2-4)對于pv節(jié)點 (3-2-5)對于平衡節(jié)點 平衡節(jié)點只設一個,電壓為已知,不參見迭代,其電壓為: (3-2-6). 修正方程式(2-3-5)和(2-3-6)兩組迭代式工包括2(n-1)個方程.選定電壓初值及變量修正量符號之后代入式(2-3-5)和(2-3-6),并將其按泰勒級數(shù)展開,略去二次方程及以后各項,得到修正方程如下: (3-2

32、-7) (3-2-8).雅可比矩陣各元素的算式式(3-2-8)中, 雅可比矩陣中的各元素可通過對式(3-2-4)和(3-2-5)進行偏導而求得.當時, 雅可比矩陣中非對角元素為 (3-2-9)當時,雅可比矩陣中對角元素為: (3-2-10)由式(3-2-9和(3-2-10)看出,雅可比矩陣的特點:矩陣中各元素是節(jié)點電壓的函數(shù),在迭代過程中,這些元素隨著節(jié)點電壓的變化而變化;導納矩陣中的某些非對角元素為零時,雅可比矩陣中對應的元素也是為零.若,則必有;雅可比矩陣不是對稱矩陣;雅可比矩陣各元素的表示如下: 2.5牛頓拉夫遜法的程序框圖第四章 matlab概述4.1matlab簡介目前電子計算機已廣

33、泛應用于電力系統(tǒng)的分析計算，潮流計算是其基本應用軟件之一?，F(xiàn)有很多潮流計算方法。對潮流計算方法有五方面的要求：（1）計算速度快（2）內(nèi)存需要少（3）計算結果有良好的可靠性和可信性（4）適應性好，亦即能處理變壓器變比調整、系統(tǒng)元件的不同描述和與其它程序配合的能力強（5）簡單。 matlab是一種交互式、面向對象的程序設計語言，廣泛應用于工業(yè)界與學術界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。matlab程序設計語言結構完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計算問題，特別是關于矩陣和矢量的計算

34、。matlab與c語言和fortran語言相比更容易被掌握。通過m語言，可以用類似數(shù)學公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構思而不是編程上。另外，matlab提供了一種特殊的工具：工具箱（toolboxes）.這些工具箱主要包括：信號處理（signal processing）、控制系統(tǒng)（control systems）、神經(jīng)網(wǎng)絡（neural networks）、模糊邏輯(fuzzy logic)、小波(wavelets)和模擬（simulation）等等。不同領域、不同層次的用戶通過相應工具的學習和應用，可以方便地進行計算、分析及設計工作

35、。matlab設計中，原始數(shù)據(jù)的填寫格式是很關鍵的一個環(huán)節(jié)，它與程序使用的方便性和靈活性有著直接的關系。原始數(shù)據(jù)輸入格式的設計，主要應從使用的角度出發(fā)，原則是簡單明了，便于修改。4.2矩陣的運算矩陣是matlab數(shù)據(jù)存儲的基本單元，而矩陣的運算是matlab語言的核心，在matlab語言系統(tǒng)中幾乎一切運算均是以對矩陣的操作為基礎的。矩陣的基本數(shù)學運算包括矩陣的四則運算、與常數(shù)的運算、逆運算、行列式運算、秩運算、特征值運算等基本函數(shù)運算，這里進行簡單介紹。四則運算矩陣的加、減、乘運算符分別為“+，*” ，用法與數(shù)字運算幾乎相同，但計算時要滿足其數(shù)學要求 在matlab中矩陣的除法有兩種形式：左除

36、“”和右除“/”。在傳統(tǒng)的matlab算法中，右除是先計算矩陣的逆再相乘，而左除則不需要計算逆矩陣直接進行除運算。通常右除要快一點，但左除可避免被除矩陣的奇異性所帶來的麻煩。在matlab6中兩者的區(qū)別不太大。與常數(shù)的運算 常數(shù)與矩陣的運算即是同該矩陣的每一元素進行運算。但需注意進行數(shù)除時，常數(shù)通常只能做除數(shù)。基本函數(shù)運算矩陣的函數(shù)運算是矩陣運算中最實用的部分，常用的主要有以下幾個：det(a) 求矩陣a的行列式eig(a) 求矩陣a的特征值inv(a)或a (-1) 求矩陣a的逆矩陣rank(a) 求矩陣a的秩trace(a) 求矩陣a的跡（對角線元素之和）我們在進行工程計算時常常遇到矩陣對

37、應元素之間的運算。這種運算不同于前面講的數(shù)學運算，為有所區(qū)別，我們稱之為數(shù)組運算?；緮?shù)學運算數(shù)組的加、減與矩陣的加、減運算完全相同。而乘除法運算有相當大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對應元素之間的乘除法，它們的運算符為“.*”和“./”或“.”。前面講過常數(shù)與矩陣的除法運算中常數(shù)只能做除數(shù)。在數(shù)組運算中有了“對應關系”的規(guī)定，數(shù)組與常數(shù)之間的除法運算沒有任何限制。另外，矩陣的數(shù)組運算中還有冪運算（運算符為 . ）、指數(shù)運算（exp）、對數(shù)運算（log）、和開方運算（sqrt）等。有了“對應元素”的規(guī)定，數(shù)組的運算實質上就是針對數(shù)組內(nèi)部的每個元素進行的。矩陣的冪運算與數(shù)組的冪運算有很大的區(qū)

38、別。邏輯關系運算 邏輯運算是matlab中數(shù)組運算所特有的一種運算形式，也是幾乎所有的高級語言普遍適用的一種運算。 4.3牛頓拉夫遜法潮流計算程序clccleare=1.06 1 1 1 1;%給功率和電壓初值。f=0 0 0 0 0;p2=0.2;q2=0.2;p3=-0.45;q3=-0.15;p4=-0.4;q4=-0.05;p5=-0.6;q5=-0.1;y =6.25-18.75i -5+15i -1.25+3.75i 0 0 -5+15i 10.834-32.5i -1.667+5i -1.667+5i -2.5+7.5i -1.25+3.75i -1.667+5i 12.917-

39、38.75i -10+30i 0 0 -1.667+5i -10+30i 12.917-38.75i -1.25+3.75i 0 -2.5+7.5i 0 -1.25+3.75i 3.75-11.25i;disp(節(jié)點導納矩陣y為:);disp(y);g=real(y);b=imag(y);k=0;%迭代次數(shù)。for v=1:10i=0,0;0,0;0,0;0,0;0,0;for m=1:5for n=1:5i(m,1)=i(m,1)+g(m,n)*e(n)-b(m,n)*f(n);%計算。 i(m,2)=i(m,2)+g(m,n)*f(n)+b(m,n)*e(n);%計算 。endendh=;

40、n=;m=;l=;j=;p2=0.2-e(2)*i(2,1)-f(2)*i(2,2);%計算功率不平衡量。q2=0.2-f(2)*i(2,1)+e(2)*i(2,2);p3=-0.45-e(3)*i(3,1)-f(3)*i(3,2);q3=-0.15-f(3)*i(3,1)+e(3)*i(3,2);p4=-0.4-e(4)*i(4,1)-f(4)*i(4,2);q4=-0.05-f(4)*i(4,1)+e(4)*i(4,2);p5=-0.6-e(5)*i(5,1)-f(5)*i(5,2);q5=-0.1-f(5)*i(5,1)+e(5)*i(5,2);for m=2:5 for n=2:5 i

41、f(m=n) h(m,m)=-b(m,m)*e(m)+ g(m,m)*f(m)+i(m,2); n(m,m)=g(m,m)*e(m)+b(m,m)*f(m)+i(m,1); m(m,m)=-g(m,m)*e(m)-b(m,m)*f(m)+i(m,1); l(m,m)=-b(m,m)*e(m)+g(m,m)*f(m)-i(m,2); else h(m,n)=-b(m,n)*e(m)+g(m,n)*f(m); n(m,n)= g(m,n)*e(m)+b(m,n)*f(m); m(m,n)=-n(m,n); l(m,n)=h(m,n);%計算雅可比矩陣的元素。 end endendj=h(2,2),

42、n(2,2),h(2,3),n(2,3),h(2,4),n(2,4),h(2,5),n(2,5);m(2,2),l(2,2),m(2,3),l(2,3),m(2,4),l(2,4),m(2,5),l(2,5);h(3,2),n(3,2),h(3,3),n(3,3),h(3,4),n(3,4),h(3,5),n(3,5);m(3,2),l(3,2),m(3,3),l(3,3),m(3,4),l(3,4),m(3,5),l(3,5);h(4,2),n(4,2),h(4,3),n(4,3),h(4,4),n(4,4),h(4,5),n(4,5);m(4,2),l(4,2),m(4,3),l(4,3)

43、,m(4,4),l(4,4),m(4,5),l(4,5);h(5,2),n(5,2),h(5,3),n(5,3),h(5,4),n(5,4),h(5,5),n(5,5);m(5,2),l(5,2),m(5,3),l(5,3),m(5,4),l(5,4),m(5,5),l(5,5);%形成雅可比矩陣。disp(雅克比矩陣j:); disp(j);c=p2;q2;p3;q3;p4;q4;p5;q5%形成功率不平衡量矩陣。a=inv(j)*c%解修正方程式。wucha=max(abs(a);if (wucha0.00001) for n=2:5 f(n)=f(n) +a(2*n-3); e(n)=e

44、(n) +a(2*n-2); %計算新值 endu=e+f*i; disp(節(jié)點電壓的第c(k)次近似值:); disp(u); disp(各點的電壓實部di(單位:v)為(節(jié)點號從小到大排列):); disp(e); disp(各點的電壓虛部fi(單位:v)為(節(jié)點號從小到大排列):);disp(f);k=k+1; disp(迭代次數(shù):); disp(k);elsebreak;endenddisp(最后電壓誤差:);disp(wucha);s1=(1.06+0*i)*(g(1,1)-i*b(1,1)*(1.06-0*i)+(g(1,2)-i*b(1,2)*(e(2)-i*f(2)+(g(1,3)-i*b(1,3)*(e(3)-i*f(3) %計算平衡節(jié)點的功率s12=(1.06+0*i)*(1.06-e(2)+i*f(2)*(-g(1,2)+i*b(1,2) %計算各支路功率s21=(e(2)+i*f(2)*e(2)-i*f(2)-1.06*(-g(2,1)+i*b(2,1)s13=(1.06+0*i)*(1.06-e(3)+i*f(3)*(-g(1,3)+i*b(1,3)s31=(e(3)+i*f(3)*e(3)-i*f(3)-1.06*(-g(3,1)+i