《畢業(yè)設計論文 電力系統(tǒng)潮流計算》由會員分享，可在線閱讀，更多相關《畢業(yè)設計論文 電力系統(tǒng)潮流計算（40頁珍藏版）》請在裝配圖網上搜索。

1、摘要 電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種重要的分析計算，它根據給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各局部的運行狀態(tài)。在電力系統(tǒng)規(guī)劃設計和現有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性，可靠性和經濟性。MATLAB使用方便，有著其他高級語言無法比較的強大的矩陣處理功能。這樣使MATLAB成為電力系統(tǒng)潮流計算的首選計算機語言。 牛頓-拉夫遜法師電力系統(tǒng)潮流計算的常用算法之一，它收斂性好，迭代次數少。介紹了電力系統(tǒng)潮流計算機輔分析的根本知識及潮流計算牛頓-拉普遜法，最后介紹了利用matlab GUI 制作潮流計算軟件的過程。

2、 關鍵詞：電力系統(tǒng)潮流計算；牛頓-拉普遜法；matlab GUI Abstract Power Flow Calculation of Power System is an important analysis and calculation of power system steady-state operation, which according to the given operating conditions and system wiring to determine the various parts

3、of the power system running state. In the study of power system design and the current operation mode are required Power Flow Calculation to quantitatively analyzed and compared to the program or run mode power supply reasonable, reliability and economy or not. MATLAB is easy to use, the powerful m

4、atrix processing is the other high-level language can not be compared with. This allows MATLAB to become the preferred computer language of power flow calculation. Newton Raphson power flow calculation is one of the most commonly used algorithms, which has good convergence and fewer iterations .Thi

5、s article describes the power flow computer assisted analysis of the basic knowledge and power flow Newton - Raphson method, Finally it describe how to use matlab GUI to make the power flow calculation software. Keywords: power flow calculation; Newton - Raphson method; matlab GUI

6、 目錄 第1章 緒論 1 課題背景 1 選題意義 1 潮流計算及其現狀及其開展趨勢 2 本畢業(yè)設計主要工作 3 第2章 電力系統(tǒng)潮流計算根本原理 4 電力網絡的數學模型 4 電力網絡的根本方程式 4 2.1.2 自導納和互導納確實定方法 5 2.1.3 節(jié)點導納矩陣的性質及意義 7 2.1.4 非標準變比變壓器等值電路 8 潮流計算的數學模型 10 2.2.1 潮流計算的節(jié)點類型 10 2.2.2 潮流計算根本方程 10 潮流計算的約束條件 12 2.4 潮流計算方法 13 2.4.1

7、 牛頓——拉夫遜法 13 2.4.2 高斯——賽德爾法 13 2.4.3 PQ分解法 14 2.4.4 擬牛頓算法 16 Matlab簡介 16 2.5.1 Matlab概述 16 2.5.2 matlab GUI 簡介 16 2.5.3 GUI 設計模板及設計窗口 17 2.5.4 GUI 設計的根本操作 17 第3章 牛頓拉夫遜潮流計算理論分析 18 概述 18 牛頓法根本原理 18 牛頓法潮流計算方程 22 節(jié)點功率方程 22 3.3.2 修正方程 23 牛頓法潮流計算主要流程 26 第4章 基于matlab潮流計算軟件的實現 28 登

8、陸界面的設計實現 28 潮流計算主界面設計實現 28 4.2.1 主界面介紹 28 4.2.2 數據初始化 29 4.2.3 潮流計算 30 4.2.4 數據處理 32 數據的傳遞問題 32 第5章 實例仿真與分析 33 實例仿真 33 運行結果分析 34 第6章 小結 35 第1章 緒論 1.1 課題背景 電力是衡量一個國家經濟開展的主要指標，也是反映人民生活水平的重要標志，它已成為現代工農業(yè)生產、交通運輸以及城鄉(xiāng)生活等許多方面不可或缺的能源和動力。電力系統(tǒng)是由發(fā)電、輸電、變電、配電和用電等環(huán)節(jié)組成的電能生產與消費系統(tǒng)。它的功能是將自然界的一次能

9、源通過發(fā)電動力裝置轉化成電能，再經輸電、變電和配電將電能供給到各用戶。為實現這一功能，電力系統(tǒng)在各個環(huán)節(jié)和不同層次還具有相應的信息與控制系統(tǒng)，對電能的生產過程進行測量、調節(jié)、控制、保護、通信和調度，以保證用戶獲得平安、經濟、優(yōu)質的電能。 電力系統(tǒng)的出現，使電能得到廣泛應用，推動了社會生產各個領域的變化，開創(chuàng)了電力時代，出現了近代史上的第二次技術革命。20世紀以來，電力系統(tǒng)的開展使動力資源得到更充分的開發(fā)，工業(yè)布局也更為合理，使電能的應用不僅深刻地影響著社會物質生產的各個側面，也越來越廣地滲透到人類日常生活的各個層面。電力系統(tǒng)的開展程度和技術水準已成為各國經濟開展水平的標志之一。 潮流

10、計算是在給定電力系統(tǒng)網絡結構、參數和決定系統(tǒng)運行狀態(tài)的邊界條件的情況下確定系統(tǒng)穩(wěn)態(tài)運行狀態(tài)的一種根本方法，是電力系統(tǒng)規(guī)劃和運營中不可缺少的一個重要組成局部?？梢哉f，它是電力系統(tǒng)分析中最根本、最重要的計算，是系統(tǒng)平安、經濟分析和實時控制與調度的根底。是電力系統(tǒng)研究人員長期研究的一個課題。MATLAB自1980年問世以來，它的強大的矩陣處理功能給電力系統(tǒng)的分析、計算帶來許多方便。在處理潮流計算時，其計算機軟件的速度已無法滿足大電網模擬和實時控制的仿真要求，而高效的潮流問題相關軟件的研究已成為大規(guī)模電力系統(tǒng)仿真計算的關鍵。隨著計算機技術的不斷開展和成熟，對MATLAB潮流計算的研究為快速、詳細地解決

11、大電網的計算問題開辟了新思路。 1.2 選題意義 電力系統(tǒng)已經與我們的生活息息相關，不可分割。進行電力系統(tǒng)潮流計算是保證電力系統(tǒng)正常運行的必要計算。具體來講電力系統(tǒng)潮流計算具有以下意義： (1) 在電網規(guī)劃階段,通過潮流計算,合理規(guī)劃電源容量及接入點,合理規(guī)劃網架,選擇無功補償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調峰、調相、調壓的要求。 　　 (2) 在編制年運行方式時,在預計負荷增長及新設備投運根底上,選擇典型方式進行潮流計算,發(fā)現電網中薄弱環(huán)節(jié),供調度員日常調度控制參考,并對規(guī)劃、基建部門提出改良網架結構,加快基建進度的建議。 　 (3) 正常檢修及特殊運行方式下的潮

12、流計算,用于日運行方式的編制,指導發(fā)電廠開機方式,有功、無功調整方案及負荷調整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質量要求。 　 (4) 預想事故、設備退出運行對靜態(tài)平安的影響分析及作出預想的運行方式調整方案?？偨Y為在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經濟性。同時，為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài)，也需要進行大量而快速的潮流計算。 因此，潮流計算是電力系統(tǒng)中應用最廣泛、最根本和最重要的一種電氣運算。 基于電力系統(tǒng)計算對保證電力系統(tǒng)正常運行具有如此正要的意義，這就要求我們能夠快速準確的進行潮流計算，計算機技術的開展使

13、電力系統(tǒng)機輔分析成為可能，各種潮流計算軟件也相繼出現。MATLAB使用方便，有著其他高級語言無法比較的強大的矩陣處理功能。MATLAB擁有600多個工程數學運算函數，可實現潮流計算的矩陣求積、求逆、稀疏矩陣形成、復數運算以及初等數學運算。同時MATLAB語言允許用戶以數學形式的語言編寫程序，這樣編程的工作量就大為減少。要到達較高的計算精度，且兼顧矩陣程序設計的難易程度，使MATLAB成為首選潮流計算的計算機語言。 因此本次設計提出了基于MATLAB潮流計算軟件的分析與設計。 該軟件能快速準確的對電力系統(tǒng)潮流進行計算，并具有一定的輔助分析功能。 1.3 潮流計算及其現狀及其開展趨勢

14、 利用電子計算機進行潮流計算從20世紀50年代中期就已經開始。此后，潮流計算曾采用了各種不同的方法，這些方法的開展主要是圍繞著對潮流計算的一些根本要求進行的。電力系統(tǒng)潮流計算屬于穩(wěn)態(tài)分析范疇，不涉及系統(tǒng)元件的動態(tài)特性和過渡過程。因此其數學模型不包含微分方程，是一組高階非線性方程。非線性代數方程組的解法離不開迭代，因此，潮流計算方法首先要求它是能可靠的收斂，并給出正確答案。隨著電力系統(tǒng)規(guī)模的不斷擴大，潮流問題的方程式階數越來越高，目前已到達幾千階甚至上萬階，對這樣規(guī)模的方程式并不是采用任何數學方法都能保證給出正確答案的。這種情況促使電力系統(tǒng)的研究人員不斷尋求新的更可靠的計算方法。知道現在潮流算法

15、的研究仍然非?；顫?，但是大多數研究都是圍繞改良牛頓法和P-Q分解法進行的。此外，隨著人工智能理論的開展，遺傳算法、人工神經網絡、模糊算法也逐漸被引入潮流計算。但是，到目前為止這些新的模型和算法還不能取代牛頓法和P-Q分解法的地位。由于電力系統(tǒng)規(guī)模的不斷擴大，對計算速度的要求不斷提高，計算機的并行計算技術也將在潮流計算中得到廣泛的應用，成為重要的研究領域。 通過幾十年的開展，潮流算法日趨成熟。近幾年，對潮流算法的研究仍然是如何改善傳統(tǒng)的潮流算法，即高斯-塞德爾法、牛頓法和快速解耦法。牛頓法，由于其在求解非線性潮流方程時采用的是逐次線性化的方法，為了進一步提高算法的收斂性和計算速度，人們考

16、慮采用將泰勒級數的高階項或非線性項也考慮進來，于是產生了二階潮流算法。后來又提出了根據直角坐標形式的潮流方程是一個二次代數方程的特點，提出了采用直角坐標的保存非線性快速潮流算法。巖本伸一等提出了一種保存非線性的快速潮流計算法,但用的是指教坐標系，因而沒法利用P-Q解耦。為了更有利于大電網的潮流計算,將此原理推廣用于P-Q解耦。這樣,既利用了保存非線性的快速算法,在迭代中使用常數雅克比矩陣，又保存了P-Q解耦的優(yōu)點。對于一些病態(tài)系統(tǒng)，應用非線性潮流計算方法往往會造成計算過程的振蕩或者不收斂，從數學上講，非線性的潮流計算方程組本來就是無解的。這樣，人們提出來了將潮流方程構造成一個函數，求此函數的最

17、小值問題，稱之為非線性規(guī)劃潮流的計算方法。優(yōu)點是原理上保證了計算過程永遠不會發(fā)散。如果將數學規(guī)劃原理和牛頓潮流算法有機結合一起就是最優(yōu)乘子法。另外，為了優(yōu)化系統(tǒng)的運行，從所有以上的可行潮流解中挑選出滿足一定指標要求的一個最正確方案就是最優(yōu)潮流問題。最優(yōu)潮流是一種同時考慮經濟性和平安性的電力網絡分析優(yōu)化問題。OPF 在電力系統(tǒng)的平安運行、經濟調度、可靠性分析、能量管理以及電力定價等方面得到了廣泛的應用。可信域和線性搜索方法是保證最優(yōu)化算法全局收斂性能的兩類技術，將內點法和可信域、線性搜索方法有機結合，構造新的優(yōu)化算法，是數學規(guī)劃領域的研究熱點。對于一些特殊性質的潮流計算問題有直流潮流計算方法、隨

18、機潮流計算方法和三相潮流計算方法。 1.4 本畢業(yè)設計主要工作 本文致力于研究分析電力網絡的運行情況。結合電力系統(tǒng)潮流計算的特點，設計一款基于matlab的潮流計算軟件，該軟件能夠進行電力系統(tǒng)潮流計算并且具有一定的輔助分析功能。具體來講要完成如下工作： 1：研究電力系統(tǒng)潮流計算的根本原理和根本方法。 2：完成電力系統(tǒng)網絡的數學建模。 3：利用matlab的M語言進行編程實現電力系統(tǒng)的潮流計算。 4：利用matlab GUI 完成軟件的登陸界面及主界面的制作 5：利用該軟件進行某電力系統(tǒng)的潮流計算，并對計算結果進行分析以驗證該軟件的可用性。

19、 第2章 電力系統(tǒng)潮流計算根本原理 2.1 電力網絡的數學模型 電力網絡的根本方程式 電力網絡可以用結點方程式或回路方程式表示出來。在結點方程式中表示網絡狀態(tài)的變量是各節(jié)點的電壓，在回路方程式中是各回路中的回路電流。 一般假設給出網絡的支路數b，結點數n，那么回路方程式數m為 m=b-n+1 結點方程式數為 =n-1 因此，回路方程式數比結點方程式數多 d=m-=b-2n+2 在一般電力系統(tǒng)中，各結點(母線)和大地間有發(fā)電機、負荷 、線路電容等對地支路，還有結點和結點之間也有輸電線路和變壓器之路，一般b

20、>2n，用結點方程式表示比用回路方程式表示方程式數目要少。而且如以下所示，用結點方程式表示容易建立直觀的方程式，輸電線的連接狀態(tài)等變化時也很容易變更網絡方程式?；谏鲜隼碛?，電力系統(tǒng)的根底網絡方程式一般都用結點方程式表示。 如圖2-1所示, 2 1 Net k N 圖2-1 把電力系統(tǒng)的發(fā)電機端子和負荷端子〔同步調相機等的端子也作為發(fā)電機端來處理〕抽出來，剩下的輸電線路及其它輸電

21、系統(tǒng)概括為網絡Ｎet表示 。在發(fā)電機結點和負荷結點上標出任意順序的記號：1,2,…,I,…,n.在輸電系統(tǒng)Net的內部不包含電源，并且各節(jié)點和大地間連接的線路對地電容、電力電容器等都作為負荷來處理。 令端子1,2……,n的對地電壓分別為，由各端子流向輸電系統(tǒng)Net 的電流相應為，那么此網絡方程組可以表示為 (2-1) (2-1)式可以簡單寫成 (I=1,2,…,n) (2-2) 或者寫成 I＝YV

22、 (2-3) 其中 (2-4) (2-4)的Y稱為節(jié)點導納矩陣。因輸電系統(tǒng)Net只是由無源元件構成的，而導納矩陣是對稱矩陣，于是有以下關系 〔2-5〕 電壓V和電流Ｉ的關系用式(2-1)～(2-5) 表示時稱為節(jié)點導納方程式。這里電壓Ｖ用電流Ｉ的方程式表示時，那么(2-3)式化為 V＝ZI (2-6) 其中 (2-6)式稱為結點阻抗

23、方程式，當然，阻抗矩陣也是對稱矩陣。 自導納和互導納確實定方法 電力網絡的節(jié)點電壓方程： (2-7) 式(2-7)為節(jié)點注入電流列向量，注入電流有正有負，注入網絡的電流為正，流出網絡的電流為負。根據這一規(guī)定，電源節(jié)點的注入電流為正，負荷節(jié)點為負。既無電源又無負荷的聯絡節(jié)點為零，帶有地方負荷的電源節(jié)點為二者代數之和。 式(2-7)為節(jié)點電壓列向量，由于節(jié)點電壓是對稱于參考節(jié)點而言的，因而需先選定參考節(jié)點。在電力系統(tǒng)中一般以地為參考節(jié)點。如整個網絡無接地支路，那么需要選定某一節(jié)點為參考。設網絡

24、中節(jié)點數為〔不含參考節(jié)點〕，那么，均為n*n列向量。為n*n階節(jié)點導納矩陣。 節(jié)電導納矩陣的節(jié)點電壓方程：，展開為： (2-8) 是一個n*n階節(jié)點導納矩陣，其階數就等于網絡中除參考節(jié)點外的節(jié)點數。 節(jié)點導納矩陣的對角元素 (i=1，2，n)成為自導納。自導納數值上就等于在i節(jié)點施加單位電壓，其他節(jié)點全部接地時，經節(jié)點i注入網絡的電流，因此，它可以定義為： 〔2-9〕 節(jié)點i的自導納數值上就等于與節(jié)點直接連接的所有支路導納的總和。 節(jié)點導納矩陣的非對角元素 (j=1，2，…，n;i=1，2，…

25、。，n;j=i)稱互導納，由此可得互導納數值上就等于在節(jié)點i施加單位電壓，其他節(jié)點全部接地時，經節(jié)點j注入網絡的電流，因此可定義為： (2-10) 節(jié)點j，i之間的互導納數值上就等于連接節(jié)點j，i支路到導納的負值。顯然，恒等于。互導納的這些性質決定了節(jié)點導納矩陣是一個對稱稀疏矩陣。而且，由于每個節(jié)點所連接的支路數總有一個限度，隨著網絡中節(jié)點數的增加非零元素相對愈來愈少，節(jié)點導納矩陣的稀疏度，即零元素數與總元素的比值就愈來愈高。 節(jié)點導納矩陣的性質及意義 節(jié)點導納矩陣的性質： 〔1〕為對稱矩陣，=。如網絡中含有源元件，如移相變壓器，那么對

26、稱性不再成立。 〔2〕對無接地支路的節(jié)點，其所在行列的元素之和均為零，即 。對于有接地支路的節(jié)點，其所在行列的元素之和等于該點接地支路的導納。利用這一性質，可以檢驗所形成節(jié)點導納矩陣的正確性。 〔3〕具有強對角性：對角元素的值不小于同一行或同一列中任一元素。 〔4〕為稀疏矩陣，因節(jié)點i ，j 之間無支路直接相連時=0，這種情況在實際電力系統(tǒng)中非常普遍。矩陣的稀疏性用稀疏度表示，其定義為矩陣中的零元素與全部元素之比，即 ， 式中Z 為中的零元素。S 隨節(jié)點數n 的增加而增加：n=50，S可達92%；n=100，S 可達90%；n=500，S可達99%，充分利用節(jié)點導納矩陣的稀疏性可節(jié)省計

27、算機內存，加快計算速度，這種技巧稱為稀疏技術。 節(jié)點導納矩陣的意義： 是n*n階方陣，其對角元素 (i=1，2，----n)稱為自導納，非對角元素(i，j=1，2，n， )稱為互導納。將節(jié)點電壓方程展開為： 〔2-11〕 可見 (2-12) 說明，自導納在數值上等于僅在節(jié)點i施加單位電壓而其余節(jié)點電壓均為零〔即其余節(jié)點全部接地〕時，經節(jié)點i注入網絡的電流。其顯然等于與節(jié)點i直接相連的所有支路的導納之和。同時可見。說明，互導納在數值上等于僅在節(jié)點j施加單位電壓而其余節(jié)點電壓均為零時，經節(jié)點i注入網絡的電流，其顯然等于()即=。為支路的導納，負號

28、表示該電流流出網絡。如節(jié)點ij之間無支路直接相連，那么該電流為0，從而=0。 注意字母幾種不寫法的不同意義：粗體黑字表示導納矩陣，大寫字母代矩陣中的第i行第j列元素，即節(jié)點i和節(jié)點j之間的互導納。小寫字母i，j支路的導納等于支路阻抗的倒數數，。 根據定義直接求取節(jié)點導納矩陣時，注意以下幾點： 1) 節(jié)點導納矩陣是方陣，其階數就等于網絡中除去參考節(jié)點外的節(jié)點數。參考節(jié)點一般取大地，編號為零。 2) 節(jié)點導納矩陣是稀疏矩陣，其各行非零非對角元素就等于與該行相對應節(jié)點所連接的不接地支路數。 3) 節(jié)點導納矩陣的對角元素就等于各該節(jié)點所連接導納的總和。因此，與沒有接地支路的節(jié)點對應的行或列中

29、，對角元素為非對角元素之和的負值。 4) 節(jié)點導納矩陣的非對角元素等于連接節(jié)點i，j支路導納的負值。因此，一般情況下，節(jié)點導納矩陣的對角元素往往大于非對角元素的負值。 5) 節(jié)點導納矩陣一般是對稱矩陣，這是網絡的互易特性所決定的。從而，一般只要求求取這個矩陣的上三角或下三角局部。 非標準變比變壓器等值電路 變壓器型等值電路更便于計算機反復計算,更適宜于復雜網絡的潮流計算.雙繞組變壓器可用阻抗與一個理想變壓器串聯的電路表示.理想變壓器只是一個參數,那就是變比。現在變壓器阻抗按實際變比歸算到低壓側為例,推導出變壓器型等值電路。 圖2-2雙繞組變壓器原理圖 圖2-3

30、變壓器阻抗歸算到低壓側等值模型 流入和流出理想變壓器的功率相等 (2-13) 式(2-13)中, 是理想變壓器的變比,和 分別為變壓器高,低繞組的實際電壓.從圖2-3直接可得： 〔2-14〕 從而可得: 〔2-15〕 式〔2-14〕中,又因節(jié)點電流方程應具有如下形式: 〔2-16〕 將式〔2-14〕與

31、〔2-15〕比較，得：，； ，。 因此可得各支路導納為： 〔2-17〕 由此可得用導納表示的變壓器型等值電路： 圖2-4變壓器型等值電路 2.2 潮流計算的數學模型 潮流計算的節(jié)點類型 用一般的電路理論求解網絡方程，目的是給出電壓源(或電流源)研究網絡內的電流(或電壓)分布，作為根底的方程式，一般用線性代數方程式表示。然而在電力系統(tǒng)中，給出發(fā)電機或負荷連接母線上電壓或電流(都是向量)的情況是很少的，一般是給出發(fā)電機母線上發(fā)電機的有功功率(P)和母線電壓的幅值(U)，給出負荷母線上負荷消耗的有功功率(P)和無功功率(Q)。主要目的是由這些量

32、去求電力系統(tǒng)內的各種電氣量。所以，根據電力系統(tǒng)中各節(jié)點性質的不同，很自然地把節(jié)點分成三類： （1） 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é)點。 （2） PU節(jié)點 這類節(jié)點給出的參數是該節(jié)點的有功功率P及電壓幅值U，待求量為該節(jié)點的無功功率Q及電壓向量的相角。這類節(jié)點在運行中往往要有一定可調節(jié)的無功電源。用以維持給定的電壓值。通常選擇有一定無功功率儲藏

33、的發(fā)電機母線或者變電所有無功補償設備的母線做PU節(jié)點處理。PU節(jié)點上的發(fā)電機稱為PU機(或PU給定型發(fā)電機) （3） 平衡節(jié)點 在潮流計算中，這類節(jié)點一般只設一個。對該節(jié)點，給定其電壓值，并在計算中取該節(jié)點電壓向量的方向作為參考軸，相當于給定該點電壓向量的角度為零。也就是說，對平衡節(jié)點給定的運行參數是U和，因此有城為U節(jié)點，而待求量是該節(jié)點的P。Q，整個系統(tǒng)的功率平衡由這一節(jié)點承當。 關于平衡節(jié)點的選擇，一般選擇系統(tǒng)中擔任調頻調壓的某一發(fā)電廠(或發(fā)電機)，有時也可能按其他原那么選擇，例如，為提高計算的收斂性?？梢赃x擇出線數多或者靠近電網中心的發(fā)電廠母線作平衡節(jié)點。 以上三類節(jié)點4個運行

34、參數P、Q、U、中，量都是兩個，待求量也是兩個，只是類型不同而已。 潮流計算根本方程 電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最根本的計算，是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為根底。 采用導納矩陣時，節(jié)點注入電流和節(jié)點電壓構成如式(2-7)所示線性方程組可展開如下形式： (2

35、-18) 由于實際電網中測量的節(jié)點注入量一般不是電流而是功率，因此必須將式中的注入電流用節(jié)點注入功率來表示。 節(jié)點功率與節(jié)點電流之間的關系為： (2-19) 式中， 因此用導納矩陣時，PQ節(jié)點可以表示為 把這個關系代入式中 ，得 〔2-20〕 式〔2-20〕就是電力系統(tǒng)潮流計算的數學模型-----潮流方程。它具有如下特點： 1：它是一組代數方程，因而表征的是電力系統(tǒng)的穩(wěn)定運行特性。 2：它是一組非線性方程，因而只能用迭代方法求其數值解。 3：由于方程中的電壓和導納

36、既可以表為直角坐標，又可表為極坐標，因而潮流方程有多種表達形式---極坐標形式，直角坐標形式和混合坐標形式。 〔1〕取 ，，得到潮流方程的極坐標形式： (2-21) (2) 取 ， ，得到潮流方程的直角坐標形式： (2-22) (3) 取 ，得到潮流方程的混合坐標形式： (2-23) 不同坐標形式的潮流方程適用于不同的迭代解法。例如：利用牛頓---拉夫遜迭代法求解，以直角坐標和混合坐標形式的潮流方程為方便；而P-Q解耦法是在混合坐標形式的根底上開展而成，故當然采用混合坐標形式。

37、 4: 它是一組n個復數方程，因而實數方程數為2n個但方程中共含4n個變量：P，Q，U和，i=1，2，，n，故必須先指定2n個變量才能求解。 2.3 潮流計算的約束條件 電力系統(tǒng)運行必須滿足一定的技術和經濟上的要求。這些要求構成了潮流問題中某些變量的約束條件，常用的約束條件如下： ①節(jié)點電壓應滿足小于節(jié)點最大額定電壓并大于最小額定電壓，即： (2-24) 從保證電能質量和供電平安的要求來看，電力系統(tǒng)的所有電氣設備都必須運行在額定電壓附近。PV節(jié)點電壓幅值必須按上述條件給定。因此，這一約束條件對PQ節(jié)點而言。 ②節(jié)點的有功功率和

38、無功功率應滿足小于節(jié)點最大額定功率并大于最小額定功率，即： (2-25) PQ節(jié)點的有功功率和無功功率，以及PV節(jié)點的有功功率，在給定時就必須滿足上述條件，因此，對平衡節(jié)點的P和Q以及PV節(jié)點的Q應按上述條件進行檢驗。 ③節(jié)點之間電壓的相位差應滿足小于最小額定相角差，即： (2-26) 為了保證系統(tǒng)運行的穩(wěn)定性，要求某些輸電線路兩端的電壓相位不超過一定的數值。這一約束的主要意義就在于此。 因此，潮流計算可以歸結為求解一組非線性方程組，并使其解答滿足一定的約束條件。常用的方法是迭代法和牛頓法，在計算過

39、程中，或得出結果之后用約束條件進行檢驗。如果不能滿足要求，那么應修改某些變量的給定值，甚至修改系統(tǒng)的運行方式，重新進行計算。 2.4 潮流計算方法 牛頓——拉夫遜法 牛頓法是數學中求解非線性方程式的典型方法，它是通過泰勒級數展開，忽略二階以上高階項，原理是逐次將非線性方程組線性，在屢次形成和求解修正方程，直至滿足要求，具體的內容參照第三章。 高斯——賽德爾法 高斯-塞德爾法原理比較簡單，主要以節(jié)點導納矩陣為根底。下面簡單介紹下其原理和潮流計算過程。 〔1〕高斯-塞德爾法的根本原理 設有n個聯立的非線性方程 〔2-

40、27〕 解此方程組可得 〔2-28〕 假設已經求得各變量的第k此迭代值，那么第〔k+1〕次迭代值為 〔2-29〕 只要給定變量的初值就可以按式〔2-10〕迭代計算，一直進行到所有變量都滿足收斂條件：即可。 〔2〕高斯-塞德爾潮流計算過程 假設有n個節(jié)點的電力系統(tǒng)，沒有PV節(jié)點，平衡節(jié)點編號為s，功率方程可寫成以下復數方程式： 〔2-30〕 對每一個PQ節(jié)點都可列出一個方程式，因而有n-1個方程式。在這些方程式中，注入功率和都是給定的，平衡節(jié)點電壓也是的，因而只有n-

41、1個節(jié)點的電壓為未知量，從而有可能求得唯一解。 將上式寫成高斯-塞德爾法的迭代形式 (2-31) 如系統(tǒng)內存在PV節(jié)點，假設節(jié)點p為PV節(jié)點，設定的節(jié)點電壓為Up0。假定高斯-塞德爾迭代法已完成第k次迭代，接著要做第k+1次迭代前，先按下式求出節(jié)點p的注入無功功率： (2-32) 然后代入下式，求出p點電壓 (2-33) 在迭代過程中，按上式求得的節(jié)點p的電壓大小不一定等于設定的節(jié)點電壓Up0，所有在下一次的迭代中，應以設定的Up0對電壓進

42、行修正，但其相角仍保持上式所求得的值，使得 (2-34) 如果所求得PV節(jié)點的無功功率越限，那么無功功率在限，該 PV節(jié)點轉化為PQ節(jié)點。 歸納起來，高斯-塞德爾迭代法計算潮流的步驟為： 1.設定各節(jié)點電壓的初值，并給定迭代誤差判據； 2.對每一個PQ節(jié)點，以前一次迭代的節(jié)點電壓值代入功率迭代方程式求出新值； 3．對于PV節(jié)點，求出其無功功率，并判斷是否越限，如越限那么將PV節(jié)點轉化為PQ節(jié)點； 4.判別各節(jié)點電壓前后二次迭代值相量差的模是否小于給定誤差，如不小于，那么回到第2步，繼續(xù)進行計算，否那么轉到第5步； 5.根

43、據功率方程求出平衡節(jié)點注入功率； 6求支路功率分布和支路功率損耗。 PQ分解法 PQ分解法是牛頓法的一種簡化方法，它利用了電力系統(tǒng)特有的運行特性，改良和提高了運行速度。由牛頓法的修正方程進行展開可得： (2-35) 根據電力系統(tǒng)的運行特性進行簡化： 1. 考慮到電力系統(tǒng)中有功功率分布主要受節(jié)點電壓相角的影響，無功功率分布主要受節(jié)點電壓幅值的影響，所以可以近似的忽略電壓幅值變化對有功功率和電壓相位變化對無功功率分布的影響，即 (2-36) 2. 根據

44、電力系統(tǒng)的正常運行條件還可作以下假設： 1) 電力系統(tǒng)正常運行時線路兩端的電壓相位角一般變化不大〔不超過10~20度〕； 2) 電力系統(tǒng)中一般架空線路的電抗遠大于電阻； 3) 節(jié)點無功功率相應的導納Q/U*U遠小于該節(jié)點的自導納的虛部。 用算式表示如下： 由以上假設，可得到雅克比矩陣的表達式： 〔2-37〕 修正方程式為 (2-38) U為節(jié)點電壓有效值的對角矩陣，B為電納矩陣〔由節(jié)點導納矩陣中各元素的虛部構成〕. 根據不同的節(jié)點還要做一些改變

45、： 1. 在有功功率局部，要除去與有功功率和電壓相位關系較小的因素，如不包含各輸電線路和變壓器支路等值Π型電路的對地電納。 2. 在無功功率局部，PV節(jié)點要做相應的處理。 那么修正方程表示為： (2-39) 一般，由于以上原因，B’和B’’是不相同的，但都是對稱的常數矩陣 。 PQ分解法的特點： 1. 以一個n-1階和一個n-m-1階線性方程組代替原有的2n-m-1階線性方程組； ’和B〞為對稱常數矩陣，且在迭代過程中保持不變； 3.P-Q分解法具有線性收斂特性，與牛頓-拉夫遜法相比，當收斂到同樣的精度時需要的迭代次

46、數較多； 4.P-Q分解法一般只適用于110KV及以上電網的計算。因為35KV及以下電壓等級的線路r/x比值很大，不滿足上述簡化條件，可能出現迭代計算不收斂的情況 擬牛頓算法 擬牛頓法是從牛頓法派生出來的新的算法，它一出現就引起廣泛的重視。近年來，擬牛頓法的研究十分活潑，它成為解非線性方程組及優(yōu)化問題的重要方法。它能在計算電力系統(tǒng)的潮流分布中，成功地減少每步迭代的計算量，并保持著超線性收斂速度。 2.5 Matlab簡介 2.5.1 Matlab概述 MATLAB (Matrix Laboratory)為美國Mathworks公司1983年首次推出的一套高性能的數值分析和

47、計算軟件，其功能不斷擴充，版本不斷升級。 MATLAB將矩陣運算、數值分析、圖形處理、編程技術結合在一起，為用戶提供了一個強有力的科學及工程問題的分析計算和程序設計工具，它還提供了專業(yè)水平的符號計算、文字處理、可視化建模仿真和實時控制等功能，是具有全部語言功能和特征的新一代軟件開發(fā)平臺。MATLAB具有編程效率高、用戶使用方便、擴充能力強、語句簡單，內涵豐富、高效方便的矩陣和數組運算、方便的繪圖功能等特點，給用戶帶來了極大的方便。 MATLAB 已開展成為適合眾多學科，多種工作平臺、功能強大的大型軟件。在歐美等國家的高校，MATLAB已成為線性代數、自動控制理論、數理統(tǒng)計、數字信號處理、時間

48、序列分析、動態(tài)系統(tǒng)仿真等高級課程的根本教學工具。成為攻讀學位的本科、碩士、博士生必須掌握的根本技能。在設計研究單位和工業(yè)開發(fā)部門，MATLAB被廣泛的應用于研究和解決各種具體問題。在中國，MATLAB也已日益受到重視，短時間內就將盛行起來，因為無論哪個學科或工程領域都可以從MATLAB中找到適宜的功能。 2.5.2 matlab GUI 簡介 圖形用戶界面〔GUI〕是用戶與計算機程序之間的交互方式，是用戶與計算機進行信息交流的方式。計算機在屏幕顯示圖形和文本，假設有揚聲器還可產生 聲音。用戶通過輸入設備，如：鍵盤、鼠標、跟蹤球、繪制板或麥克風，與計算機通訊。用戶界面設定了如何觀看和如何感

49、知計算機、操作系統(tǒng)或應用程序。通常， 多是根據悅目的結構和用戶界面功能的有效性來選擇計算機或程序。圖形用戶界面或GUI是包含圖形對象，如：窗口、圖標、菜單和文本的用戶界面。以某種方式 選擇或激活這些對象，通常引起動作或發(fā)生變化。最常見的激活方法是用鼠標或其它點擊設備去控制屏幕上的鼠標指針的運動。按下鼠標按鈕，標志著對象的選擇或 其它動作。 Matlab作為強大的數學計算軟件，同樣也提供了圖像用戶界面設計的功能。在matlab中，根本的圖形用戶界面對象包含3類：用戶控件對象〔uicontrol〕、下拉式菜單對象〔uimenu〕、和快捷菜單對象〔uicontexmenu〕。根據這些對象可以設計出

50、界面友好、操作方便的圖形用戶界面。 GUI 設計模板及設計窗口 Matlab為GUI設計準本了四個模板，分別是Blank GUI(默認)、GUI with Uicontronl〔帶控件對象的GUI〕、GUI with Axes and Menu〔帶坐標軸與菜單的GUI〕、Modal Question Dialog〔帶模式問話對話框的GUI模板〕，GUI 設計模板如圖3-1所示。 當用戶選擇不同模板時，在GUI設計模板界面的右邊就會顯示與該模板對應的GUI圖形。 圖2-5 GUI設計模板 選擇設計模板后就進如GUI 設計窗口，GUI設計窗口由菜單欄、工具欄、控件工具欄以及圖形對

51、象設計區(qū)組成。 在GUI設計窗口的工具欄上有位置調整器、菜單編輯器、tab順序編輯器、屬性查看器等可視化設計工具?？丶ぞ邫诎≒ush Button、Check Box、Edit Box、Popup Menu、Axes、table等控件對象，他們是構成GUI 的根本元素。 GUI 設計的根本操作 為了添加對象控件，可以從GUI 設計窗口的控件工具欄中選擇一個對象，然后以拖曳方式在對象設計區(qū)建立該對象，其對象創(chuàng)立方式方便、簡單。在GUI 設計窗口創(chuàng)立對象后，通過雙擊該對象，就會顯示該對象的屬性查看器，通過它可以設計該對象的屬性值。 在選中對象的前提下，單擊鼠標右鍵，會彈出一個快捷菜單

52、，可以從中某個子菜單進行相應的操作。在對象設計區(qū)右擊鼠標，會顯示與圖形窗口有關的快捷菜單。 第3章 牛頓拉夫遜潮流計算理論分析 3.1 概述 牛頓法收斂性好，迭代次數少，在潮流計算方法中得到廣泛的應用，目前為止還沒有更好的方法能夠完全取代它。 牛頓拉夫遜法〔下面簡稱牛頓法〕是數學中求解非線性方程的典型方法，能快速求出其他方法求不出或者難以求出的解。本章將主要針對牛頓法的理論進行具體介紹。 3.2 牛頓法根本原理 牛頓-拉夫遜法是解非線性方程式的有效方法。牛頓拉夫遜法潮流計算是目前最為廣泛、效果最好的一種潮流計算方法。這種把非線性方程式的求解過程變成反復對相應的線性方程式的求解過程，即

53、逐次線性化過程，這就是牛頓法的核心。我們以如下非線性方程式的求解過程為例來說明： 〔3-1〕 設為該方程式的初值。而真正解x在它的近旁： 〔3-2〕 式中：為初始值的修正量。如果求得，那么由式〔3-2〕就可以得到真正解x。 為此將式 〔3-3〕 按泰勒級數展開 〔3-4〕 當我們選擇的初始值比較好，即很小時，式〔3-4〕

54、中包含的和更高階次項可以略去不計。因此，式〔3-4〕可以簡化為 〔3-5〕 這是對于變量的形式方程式，用它可以求出修正量。 由于式〔3-5〕是式〔3-4〕的簡化結果，所以由式〔3-5〕解出后，還不能得到方程式〔3-1〕的真正解。實際上，用對修正后得到的： 〔3-6〕 只是向真正解更逼近一些?，F在如果再以作為初值，解式〔3-5〕 就能得到更趨近真正解的：

55、 〔3-7〕 這樣反復下去，就構成了不斷求解非線性方程式的逐次線性化過程。第t次迭代時的參數方程為 〔3-8〕 或 〔3-9〕 上式左端可以看成是近似解引起的誤差，當時，就滿足了原方程式〔3-1〕，因而就成為該方程的解。式中是函數 在點的一次導數，也就是曲線在點的斜率，如圖〔3-1〕所示，修正量那么是由點的切線與橫軸的交點來確定，由圖〔3-1〕可以直觀的看出牛頓法的求解過程。 圖3-1 牛頓-拉夫遜法幾何解釋 現在把牛頓法推廣到

56、多變量非線性方程組的情況。設有變量的非線性聯立方程組： 〔3-10〕 給定各變量初值，假設為其修正量，并使其滿足 〔3-11〕 對以上n個方程式分別按泰勒級數展開，當忽略所組成的二次項和高次項時，可以得到 〔3-12〕 式中：為函數對自變量的偏導數在點〔〕處的值。把上式寫成矩陣形式： 〔3-13〕 這是變量的線性方程組，稱為牛頓法的修正方程，通過它可以解出 ，并可以進一步求得

57、 〔3-14〕 式中向真正解逼近了一步，如果再以它們作為初值重復解式〔3-13〕修正方程式，等到更接近真解的，如此迭代下去，并按式〔3-14〕進行修正，直到滿足收斂要求為止并停止迭代計算，這就構成了牛頓法的迭代過程。 一般第t次迭代式的修正方程為 〔3-15〕 上式可以簡寫為 (3-16) 其中 ， 其中的為第t次迭代時的雅克比矩陣； 同理可以得到第t次迭代時的修正量： 〔3-17〕 同樣，也可以寫出類似〔3-14〕的算

58、式 〔3-18〕 這樣反復交替的解式〔3-16〕及式〔3-18〕就可以使逐步趨近方程式的真正解。當滿足人為收斂條件時，即 或 〔3-19〕 迭代結束，式中為預先給定的小正數。 3.3 牛頓法潮流計算方程 節(jié)點功率方程 電力系統(tǒng)的負荷習慣用功率表示，對于有n個節(jié)點的電力系統(tǒng)，系統(tǒng)中各節(jié)點注入電流與注入功率以標幺值表示的關系為 i=1，2，……，n 〔3-20〕 式中表示其共軛復數。將此關系式代入節(jié)點電壓方程的通式，可得到以節(jié)點注入功率表示的節(jié)點電壓方程:

59、 (3-21) 上述的方程式，通常稱為功率方程。根據方程中的節(jié)點電壓向量表示的不同，可以得到不同形式的功率方程。 假設節(jié)點電壓向量以直角坐標表示，即以復數平面上實軸與虛軸上的投影表示可寫成 〔3-22〕 其共軛值為 〔3-23〕 導納表示為

60、 〔3-24〕 把這兩關系式代回式〔3-21〕的功率方程中，展開后再將功率方程的實部和虛局部別寫成有功、無功功率別離的節(jié)點方功率方程： 〔3-25〕 式中：i=1，2，……，n為各節(jié)點的編號。 假設節(jié)點電壓以極坐標表示，那么 或寫成 〔3-26〕 將其同導納的復數表達式一起代入式〔3-21〕的功率方程，進整理可以得到 〔3-27〕 式中：——i與j節(jié)點電壓的相角差。 由式〔3-25〕和〔3-27〕給出的功率方程表示方法防止了復數運算，因此，在潮流計算

61、中普遍采用。 修正方程 采用牛頓法計算潮流時，需要對功率方程進行修改。下面將根據在不同坐標內的修改良行討論： 〔1〕在直角坐標系內時，由PQ節(jié)點功率方程〔3-25〕可知：節(jié)點i的注入功率是各點電壓的函數，設節(jié)點的電壓，代入式〔3-25〕，可以求出節(jié)點i的有功及無功功率，它們與給定的PQ 節(jié)點的注入功率的差值應滿足以下方程 〔3-28〕 對于PV 節(jié)點，節(jié)點的注入有功功率及節(jié)點電壓大小，記作，其節(jié)點的有功功率應滿方程： 〔3-29〕 對于平衡節(jié)點，因為其電壓給定，故不需要迭代求解。 通過以上分析可見，式〔3-28〕和式〔3-29〕共2〔n-1〕個方程，待求量共2〔n-

62、1〕個。將上述2〔n-1〕個方程按泰勒級數展開，并略去修正量的高次方項后得到修正方程如下： 〔3-30〕 其中雅克比矩陣的各元素可以對式〔3-28〕和式〔3-29〕求偏導數獲得。 對于非對角元素〔〕有 〔3-31〕 對于對角元素〔有 〔3-32〕 由上述表達式可以看到，雅克比矩陣具有以下特點： 1） 各元素是各節(jié)點電壓的函數，迭代過程中每迭代一次各節(jié)點電壓都要變化，因而各元素每次也變化； 2） 雅克比矩陣不具有對稱

63、性； 3） 互導納，與之對應的非對角元素亦為零，此外因非對角元素，故雅克比矩陣是稀疏矩。 （1） 當在極坐標系內時，由功率方程〔3-27〕可知節(jié)點i的注入功率是各節(jié)點電壓幅值和相角的函數。代入式〔3-27〕可以求出節(jié)點i的有功功率和無功功率，它們與給定的PQ節(jié)點的注入功率的差值滿足下面方程： 〔3-33〕 式中：——i與j節(jié)點電壓的相角差。 在有n個節(jié)點的系統(tǒng)中，假定第號節(jié)點為PQ節(jié)點，第m+1~n-1號節(jié)點為PV節(jié)點，第n號節(jié)點為平衡節(jié)點。和是給定的，PV節(jié)點的電壓幅值也是給定的，因此，只剩下n-1個節(jié)點的電壓相角和m個節(jié)點的電壓幅值

64、是未知量。由〔3-33〕可知一共包含了n-1+m方程式，正好同未知量的數目相等，而直角坐標形式的方程少了n-1-m個。 由方程〔3-33〕可以寫出修正方程 (3-34) 式中 〔3-35〕 其中：H是階方陣，其元素為；N是階矩陣，其元素為；K是階矩陣，其元素為；L是階矩陣，其元素為。 對式〔3-33〕求偏導數，可得雅克比矩陣元素的表達式如下： 非對角元素〔〕 〔3-36〕 對角元素〔〕

65、 〔3-37〕 3.4 牛頓法潮流計算主要流程 1) 形成節(jié)點導納矩陣； 2) 給各節(jié)點電壓設初值； 3) 將節(jié)點電壓初值代入〔3-28〕〔3-29〕，求出修正方程式的常數項向量； 4) 將節(jié)點電壓初值代入〔3-31〕，〔3-32〕，求出雅可比矩陣元素； 5) 求解修正方程式〔3-30〕，求出變量的修正向量； 6) 求出節(jié)點電壓的新值； 7) 如有PV節(jié)點，那么檢查該類節(jié)點的無功功率是否越限； 8) 檢查是否收斂，由式〔3-19〕可知，假設電壓趨近于真解時，功率偏移量將趨于零。如不收斂，那么以各節(jié)點電壓的新值作為初值自第3步重新開始下一次迭代，否那么轉入下一步。 9)

66、 計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點注入功率，最后輸出結果，并結束。 牛頓-拉夫遜潮流計算程序框圖如圖3-2所示 第4章 基于matlab潮流計算軟件的實現 4.1 登陸界面的設計實現 登陸界面是用戶對軟件進行操作的門戶，也對用戶身份進行驗證。本設計的登陸界面設計如圖4-1所示： 圖 4-1登陸界面 當用戶輸入正確的用戶信息是進入潮流計算主界面，否那么彈出錯誤提示對話框，錯誤對話框如圖4-2所示： 圖4-2 錯誤提示對話框 4.2 潮流計算主界面設計實現 主界面介紹 主界面是進行潮流計算的的操作界面，能完成數據的初始化、進行潮流計算、并顯示潮流計算的相關結果以供用戶分析。具體來講就是點擊第一個彈出框選擇潮流計算數據組，并在主界面上顯示所選擇數據組的電路接線圖，單擊計算按鈕就能完成潮流計算并將潮流計算結果保存下來，單擊第二個彈出框可以查看具體的結果，同時還將畫出潮流計算電壓迭代曲線。潮流計算主界面如圖4-3所示： 圖4-3 主界面 數據初始化 本設計將數據保存在excel表格里，在主界面上單擊彈出框選擇數據組從而從e