電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)

1、課 程 設(shè) 計(jì) 報(bào) 告學(xué)生姓名:學(xué) 號：學(xué) 院:電氣工程學(xué)院班 級:題 目:電力系統(tǒng)潮流計(jì)算初 壯指導(dǎo)教師： 職稱: 副教授 指導(dǎo)教師： 李翠萍 職稱: 副教授 2014年 01月10日摘 要潮流計(jì)算是電力系統(tǒng)最基本最常用的計(jì)算。根據(jù)系統(tǒng)給定的運(yùn)行條件，網(wǎng)絡(luò)接線及元件參數(shù)，通過潮流計(jì)算可以確定各母線的電壓幅值和相角，各元件流過的功率，整個(gè)系統(tǒng)的功率損耗。潮流計(jì)算是實(shí)現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)供電的必要手段和重要工作環(huán)節(jié)。因此，潮流計(jì)算在電力系統(tǒng)的規(guī)劃計(jì)算，生產(chǎn)運(yùn)行，調(diào)度管理及科學(xué)計(jì)算中都有著廣泛的應(yīng)用。潮流計(jì)算在數(shù)學(xué)上是多元非線性方程組的求解問題，求解的方法有很多種，牛頓拉夫遜法是數(shù)學(xué)上解非線性方程

2、組的有效方法，有較好的收斂性。將牛頓拉夫遜法用于潮流計(jì)算是以導(dǎo)納矩陣為基礎(chǔ)的，由于利用了導(dǎo)納矩陣的對稱性，稀疏性及節(jié)點(diǎn)編號順序優(yōu)劃等技巧，使牛頓拉夫遜法在收斂性，占用內(nèi)存，計(jì)算速度等方面的優(yōu)點(diǎn)都超過了阻抗法。運(yùn)用電子計(jì)算機(jī)計(jì)算一般要完成以下幾個(gè)步驟：建立數(shù)學(xué)模型，確定解算方法，制訂計(jì)算流程，編制計(jì)算程序。首先，畫出系統(tǒng)的等效電路圖，在計(jì)算出各元件參數(shù)的基礎(chǔ)上，應(yīng)用牛頓拉夫遜法以及matlab軟件進(jìn)行計(jì)算對給定系統(tǒng)圖進(jìn)行了四種不同負(fù)荷下的潮流計(jì)算，經(jīng)過調(diào)節(jié)均得到符合電壓限制及功率限制的潮流分布。其次，牛頓拉夫遜法具有較好的收斂性，上述計(jì)算過程經(jīng)過四到五次迭代后均能收斂。根據(jù)運(yùn)算結(jié)果，分析各支路損

3、耗和系統(tǒng)總損耗。關(guān)鍵詞：牛頓-拉夫遜法 matlab 潮流計(jì)算 matpow一、 題目原始資料1系統(tǒng)圖：兩個(gè)發(fā)電廠分別通過變壓器和輸電線路與四個(gè)變電所相連。變電所1變電所2變電所3變電所435kv母線10kv母線35kv母線10kv母線一次側(cè)電壓220kv一次側(cè)電壓220kv母線1母線3母線2線路長為90km線路長為80km線路長為90km線路長為100km2*qfq-50-22*qfs-50-2tqn-100-22*tqn-100-2。線路長為80km。電廠一電廠二線路長為70km2、發(fā)電廠資料：母線1和2為發(fā)電廠高壓母線，發(fā)電廠一總裝機(jī)容量為（ 300mw ），母線3為機(jī)壓母線，機(jī)壓母線上

4、裝機(jī)容量為100mw，最大負(fù)荷和最小負(fù)荷分別為20mw和50mw,發(fā)電廠二裝機(jī)容量為200mw。3、變電所資料：(一) 變電所1、2、3、4低壓母線的電壓等級分別為：35kv 10kv 35kv 10kv(二) 變電所的負(fù)荷分別為：60mw 40mw 70mw 50mw(三) 每個(gè)變電所的功率因數(shù)均為cos=0.85；(四) 變電所1和變電所3分別配有兩臺(tái)容量為75mva的變壓器，短路損耗414kw，短路電壓（%）=16.7；變電所2和變電所4分別配有兩臺(tái)容量為63mva的變壓器，短路損耗為245kw，短路電壓（%）=10.5。4、輸電線路資料：發(fā)電廠和變電所之間的輸電線路的電壓等級及長度標(biāo)于

5、圖中，單位長度的電阻為，單位長度的電抗為，單位長度的電納為。二、 課程設(shè)計(jì)內(nèi)容及要求課設(shè)內(nèi)容：1. 對給定的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖2. 輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對計(jì)算結(jié)果進(jìn)行分析。3. 跟隨變電所負(fù)荷按一定比例發(fā)生變化，進(jìn)行潮流計(jì)算分析。1) 4個(gè)變電所的負(fù)荷同時(shí)以2%的比例增大；2) 4個(gè)變電所的負(fù)荷同時(shí)以2%的比例下降3) 1和4號變電所的負(fù)荷同時(shí)以2%的比例下降，而2和3號變電所的負(fù)荷同時(shí)以2%的比例上升；4. 在不同的負(fù)荷情況下，分析潮流計(jì)算的結(jié)果，如果各母線電壓不滿足要求，進(jìn)行電壓的調(diào)整。（變電所

6、低壓母線電壓10kv要求調(diào)整范圍在9.5-10.5之間；電壓35kv要求調(diào)整范圍在35-36之間）5. 輪流斷開支路雙回線中的一條，分析潮流的分布。（幾條支路斷幾次）6. 最終形成課程設(shè)計(jì)成品說明書。課設(shè)要求：1. 在讀懂程序的基礎(chǔ)上畫出潮流計(jì)算基本流程圖2. 通過輸入數(shù)據(jù)，進(jìn)行潮流計(jì)算輸出結(jié)果3. 對不同的負(fù)荷變化，分析潮流分布，寫出分析說明。4. 對不同的負(fù)荷變化，進(jìn)行潮流的調(diào)節(jié)控制，并說明調(diào)節(jié)控制的方法，并列表表示調(diào)節(jié)控制的參數(shù)變化。5. 打印利用ddrts進(jìn)行潮流分析繪制的系統(tǒng)圖，以及潮流分布圖。三、 題目分析3.1牛頓拉夫遜法概要首先對一般的牛頓拉夫遜法作一簡單的說明。已知一個(gè)變量x

7、函數(shù)為：到此方程時(shí)，由適當(dāng)?shù)慕浦党霭l(fā)，根據(jù)：反復(fù)進(jìn)行計(jì)算，當(dāng)滿足適當(dāng)?shù)氖諗織l件就是上面方程的根。這樣的方法就是所謂的牛頓拉夫遜法。這一方法還可以做下面的解釋，設(shè)第次迭代得到的解語真值之差，即的誤差為時(shí)，則：把在附近對用泰勒級數(shù)展開上式省略去以后部分的誤差可以近似由上式計(jì)算出來。比較兩式，可以看出牛頓拉夫遜法的休整量和的誤差的一次項(xiàng)相等。用同樣的方法考慮，給出個(gè)變量的個(gè)方程：對其近似解得修正量可以通過解下邊的方程來確定：式中等號右邊的矩陣都是對于的值。這一矩陣稱為雅可比（jacobi）矩陣。按上述得到的修正向量后，得到如下關(guān)系這比更接近真實(shí)值。這一步在收斂到希望的值以前重復(fù)進(jìn)行，一般要反復(fù)計(jì)算

8、滿足為預(yù)先規(guī)定的小正數(shù)，是第n次迭代的近似值。3.2牛頓法的框圖及求解過程3.2.1用牛頓法計(jì)算潮流時(shí)，有以下的步驟：a) 給這各節(jié)點(diǎn)電壓初始值；b) 將以上電壓初始值代入公式，求修正方程的常數(shù)項(xiàng)向量 ；c) 將電壓初始值在帶入上述公式，求出修正方程中系數(shù)矩陣的各元素。d) 解修正方程式；e) 修正各節(jié)點(diǎn)電壓，；f) 將，在帶入方程式，求出；g) 檢驗(yàn)是否收斂，即h) 如果收斂，迭代到此結(jié)束，進(jìn)一步計(jì)算各線路潮流和平衡節(jié)點(diǎn)功率，并打印輸出結(jié)果。如果不收斂，轉(zhuǎn)回（2）進(jìn)行下次迭代計(jì)算，直到收斂為止。3.2.2程序框圖如下3.3 matpower的功能與使用方法matpower所用的所有數(shù)據(jù)文件均

9、為matlab的m文件或者mat文件，用來定義和返回變量、basemva、bus、branch、gen等。 basemva變量是一個(gè)標(biāo)量，用來設(shè)置基準(zhǔn)容量。對于計(jì)算中采用有名值，可以根據(jù)需要設(shè)置，如ioomva；對于計(jì)算中采用標(biāo)么值，一般設(shè)置為1。 bus變量是一個(gè)矩陣，用來設(shè)置電網(wǎng)中各母線參數(shù)，其格式為：bus_i、type、pd、qd、gs、bs、83ea、vm、va、basekv、zone、vmax、vmin。bus-i用來設(shè)置母線編號，范圍為129997。type用來設(shè)置母線類型，1為pq節(jié)點(diǎn)母線，2為pv節(jié)點(diǎn)母線，3為平衡（參考)節(jié)點(diǎn)母線。pd和qd用來設(shè)置母線注入負(fù)荷的有功和無功功

10、率。gs、bs用來設(shè)置與母線并聯(lián)電導(dǎo)和電納。basekv用來設(shè)置該母線基準(zhǔn)電壓。vm和va用來設(shè)置母線電壓的幅值、相位初值。vmax和vmin用來設(shè)置工作時(shí)母線最高、最低電壓幅值。area和zone用來設(shè)置電網(wǎng)斷面號和省耗分區(qū)號，一般都設(shè)置為1，前者可設(shè)置范圍為l100，后者可設(shè)置范圍為1999。 branch變量是一個(gè)矩陣，用來設(shè)置電網(wǎng)中各支路參數(shù)，其格式為：fbus、tbus、r、x、b、ratea、mteb、mtec、ratio、ajlgle、status。fbus和tbus用來設(shè)置該支路由起始節(jié)點(diǎn)(母線)編號和終止節(jié)點(diǎn)(母線)編號。r、x和b用來設(shè)置該支路的電阻、電抗和充電電納。rat

11、ea、 rateb和rate(；分別用來設(shè)置該支路長期、短期和緊急允許功率。ratio用來設(shè)置該支路的變比，如果支路元件僅僅是導(dǎo)線為o；如果支路元件為變壓器，則該變比為fl，us側(cè)母線的基準(zhǔn)電壓與thus側(cè)母線的基準(zhǔn)電壓之比。angle用來設(shè)置支路的相位角度，如果支路元件為變壓器，就是變壓器的轉(zhuǎn)角；如果支路元件不是變壓器，相位角度為0度。status用來設(shè)置支路工作狀態(tài)，1表示投入運(yùn)行，0表示退出運(yùn)行。 gen變量也是一個(gè)矩陣，用來設(shè)置接入電網(wǎng)中的發(fā)電機(jī)(電源)參數(shù)，其格式為：bus、魄、qg、qmax、qmin、vg、mbase、status、pmax、pmin。bus用來設(shè)置接人發(fā)電機(jī)(電

12、源)的母線編號。pg和qg用來設(shè)置接人發(fā)電機(jī)(電源)的有功和無功功率。pmax和pmin用來設(shè)置接入發(fā)電機(jī)(電源)的有功功率最大、最小允許值。9max和qmin用來設(shè)置接入發(fā)電機(jī)(電源)的無功功率最大、最小允許值。vg用來設(shè)置接人發(fā)電機(jī)(電源)的工作電壓。mbase用來設(shè)置接入發(fā)電機(jī)(電源)的功率基準(zhǔn)，如為缺省值，就是baselva變量的值。status用來設(shè)置發(fā)電機(jī)(電源)工作狀態(tài)，l表示投入運(yùn)行，0表示退出運(yùn)行。3.4 問題求解電壓是衡量電力系統(tǒng)電能質(zhì)量的標(biāo)準(zhǔn)之一。電壓過高或過低，都將對人身及其用電設(shè)備產(chǎn)生重大的影響。保證用戶的電壓接近額定值是電力系統(tǒng)調(diào)度的基本任務(wù)之一。當(dāng)系統(tǒng)的電壓偏離允

13、許值時(shí)，電力系統(tǒng)必須應(yīng)用電壓調(diào)節(jié)技術(shù)調(diào)節(jié)系統(tǒng)電壓的大小，使其維持在允許值范圍內(nèi)。本文經(jīng)過手算形成了等值電路圖，并編寫好了程序得出節(jié)點(diǎn)電壓標(biāo)幺值，使其滿足所要求的調(diào)整范圍。常規(guī)的電力系統(tǒng)潮流計(jì)算中一般具有三種類型的節(jié)點(diǎn)：pq、pv及平衡節(jié)點(diǎn)。一個(gè)節(jié)點(diǎn)有四個(gè)變量，即注入有功功率、注入無功功率，電壓大小及相角。常規(guī)的潮流計(jì)算一般給定其中的二個(gè)變量：pq節(jié)點(diǎn)（注入有功功率及無功功率），pv節(jié)點(diǎn)（注入有功功率及電壓的大?。胶夤?jié)點(diǎn)（電壓的大小及相角）。3.4.1變量的分類：負(fù)荷消耗的有功、無功功率、電源發(fā)出的有功、無功功率、母線或節(jié)點(diǎn)的電壓大小和相位、 、在這十二個(gè)變量中，負(fù)荷消耗的有功和無功功率無法

14、控制，因它們?nèi)Q于用戶，它們就稱為不可控變量或是擾動(dòng)變量。電源發(fā)出的有功無功功率是可以控制的自變量，因此它們就稱為控制變量。母線或節(jié)點(diǎn)電壓的大小和相位角是受控制變量控制的因變量。其中, 、主要受、的控制, 、主要受、的控制。這四個(gè)變量就是簡單系統(tǒng)的狀態(tài)變量。為了保證系統(tǒng)的正常運(yùn)行必須滿足以下的約束條件：對控制變量 對沒有電源的節(jié)點(diǎn)則為對狀態(tài)變量的約束條件則是 對某些狀態(tài)變量還有如下的約束條件 3.4.2節(jié)點(diǎn)設(shè)置及分類 根據(jù)系統(tǒng)圖可知此系統(tǒng)為兩端供電網(wǎng)路，將母線1，2設(shè)為節(jié)點(diǎn)1，6，將變電所1、2、3、4的高低壓側(cè)分別設(shè)為節(jié)點(diǎn)2、3、4、5、7、8、9、10。并且，將節(jié)點(diǎn)1設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn)6

15、設(shè)為pv節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為pq節(jié)點(diǎn)。3.4.3參數(shù)求取設(shè)定基準(zhǔn)值，所以根據(jù)題目原始資料，計(jì)算發(fā)電廠、變壓器及線路的參數(shù)。（1）運(yùn)用下列公式計(jì)算變壓器參數(shù)：采用變壓器參數(shù)為折算至高壓側(cè)的數(shù)值，其變比k1，其中， 75mva的變壓器： 兩臺(tái)變壓器并聯(lián)，得阻抗假設(shè)我們所采用的變壓器有5個(gè)抽頭，電壓調(diào)節(jié)范圍為2*2.5%， 對應(yīng)的分接頭開始時(shí)設(shè)變壓器高壓側(cè)主分接頭,設(shè)四個(gè)變電所變壓器的非標(biāo)準(zhǔn)變比，降壓變壓器5個(gè)分接頭的非標(biāo)準(zhǔn)變比如下，以備調(diào)壓時(shí)選用變壓器非標(biāo)準(zhǔn)變比按照上面數(shù)據(jù)輸入后到得到各個(gè)接點(diǎn)電壓標(biāo)幺值，具體檢驗(yàn)過程如下：可得出變壓器二次側(cè)的實(shí)際電壓：形成節(jié)點(diǎn)參數(shù)矩陣時(shí)，還需要我們計(jì)算出各節(jié)點(diǎn)所接發(fā)

16、電機(jī)功率，節(jié)點(diǎn)負(fù)荷功率（可根據(jù)負(fù)荷有功功率，發(fā)電機(jī)有功功率及已知的功率因數(shù)計(jì)算）變電所1 變電所2 變電所3 變電所4 平衡節(jié)點(diǎn)為節(jié)點(diǎn)，所設(shè)節(jié)點(diǎn)電壓的初始值為，給定值為1.05pv節(jié)點(diǎn)為節(jié)點(diǎn)，給定電壓值為1.05，。（2）計(jì)算線路參數(shù)支路阻抗 支路對地電容 支路70km輸電線路：支路80km輸電線路：支路90km輸電線路：支路90km同桿架設(shè)雙回線：支路80km同桿架設(shè)雙回線：支路100km同桿架設(shè)雙回線：3.5 等值電路的計(jì)算3.6代入程序并運(yùn)行從結(jié)果中可以觀察到各個(gè)節(jié)點(diǎn)電壓標(biāo)幺值分別為:節(jié)點(diǎn)12345678910電壓1.0180.9620.9911.0021.0401.0181.0251.

17、0321.0501.050折算到有名值以后，可以觀察到節(jié)點(diǎn)1、4的電壓都與題目給定的范圍相比偏高，因此調(diào)節(jié)變壓器分接頭和發(fā)電廠電壓，通過對系統(tǒng)進(jìn)行多次調(diào)整，最終得到合理結(jié)果，調(diào)整結(jié)果如下: 調(diào)節(jié)方法：調(diào)節(jié)方法分接頭1分接頭2分接頭3分接頭4未調(diào)整1111調(diào)整后1.0251.01.01.05調(diào)節(jié)結(jié)果： 節(jié)點(diǎn)號12345678910未調(diào)整1.0180.9620.9911.0021.0401.0181.0251.0321.0501.050調(diào)整后0.9920.9620.9910.9511.0401.0181.0251.0321.0501.050 具體的潮流分布見下表：各條支路的首端功率p各條支路的首端

18、功率q各條支路的末端功率p各條支路的末端功率qs(5,1)0.800.53-0.80-0.50s(6,2)1.411.03-1.40-0.87s(7,3)1.200.83-1.20-0.74s(8,4)0.800.55-0.80-0.50s(5,9)-0.80-0.530.81-0.10s(6,9)-0.22-1.090.240.87s(6,7)-1.180.061.19-0.19s(7,9)0.62-1.32-0.601.14s(7,8)-3.020.683.08-1.03s(8,10)-3.880.484.00-0.84變壓器分接頭的調(diào)高，使最終調(diào)節(jié)使得電壓在規(guī)定范圍內(nèi)時(shí)，有功損耗為11.

19、95mw。3.6.1 4個(gè)變電所負(fù)荷同時(shí)以2%的比例增大由于變電所負(fù)荷以2%的比例增大，所以各個(gè)變電所負(fù)荷變?yōu)樽冸娝? =0.816+j0.5057變電所2 =1.428+j0.8850變電所3 =1.224+j0.7585變電所4 =0.816+j0.5057變電所4 =0.816+j0.5057將程序中相應(yīng)矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，此系統(tǒng)經(jīng)過調(diào)節(jié)后達(dá)到要求，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：調(diào)節(jié)方法分接頭1分接頭2分接頭3分接頭4調(diào)整前1111調(diào)整后1.0251.01.01.05 調(diào)節(jié)結(jié)果： 節(jié)點(diǎn)電壓12345678910調(diào)整前1.0170.9600.9891.0001.0391.0171.0

20、241.0311.0501.050調(diào)整后0.9910.9600.9890.9501.0391.0171.0241.0311.0501.050 由上表觀察到，進(jìn)行最后一次調(diào)節(jié)時(shí)節(jié)點(diǎn)1,4的電壓均在題目允許的范圍內(nèi)，變壓器分接頭的調(diào)高，使最終調(diào)節(jié)使得電壓在規(guī)定范圍內(nèi)時(shí)，有功損耗為12mw。由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下： 各條支路的首端功率p各條支路的首端功率q各條支路的末端功率p各條支路的末端功率qs(5,1)0.820.54-0.82-0.51s(6,2)1.431.06-1.43-0.88s(7,3)1.230.85-1

21、.22-0.76s(8,4)0.820.56-0.82-0.51s(5,9)-0.82-0.540.82-0.08s(6,9)-0.25-1.110.270.89s(6,7)-1.180.051.19-0.18s(7,9)0.58-1.33-0.571.15s(7,8)-3.000.663.06-1.02s(8,10)-3.880.454.00-0.813.7.2 4個(gè)變電所負(fù)荷同時(shí)以2%的比例減小由于變電所負(fù)荷以2%的比例減小，所以各個(gè)變電所負(fù)荷變?yōu)椋?變電所1 =0.784+j0.4859變電所2 =1.372+j0.8502變電所3 =1.1776+j0.7287變電所4 =0.784+

22、j0.4859將程序中相應(yīng)矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，此系統(tǒng)經(jīng)過調(diào)節(jié)后達(dá)到要求，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：調(diào)節(jié)方法分接頭1分接頭2分接頭3分接頭4調(diào)整前1111調(diào)整后1.0251.01.01.05 調(diào)節(jié)結(jié)果： 節(jié)點(diǎn)電壓12345678910調(diào)整前1.0190.9640.9921.0031.0401.0191.0261.0321.0501.050調(diào)整后0.9930.9640.9920.9521.0401.0191.0261.0321.0501.050 由上表觀察到，進(jìn)行最后一次調(diào)節(jié)時(shí)節(jié)點(diǎn)1,4的電壓均在題目允許的范圍內(nèi)，變壓器分接頭的調(diào)高，使最終調(diào)節(jié)使得電壓在規(guī)定范圍內(nèi)時(shí)，有功損耗為11.95

23、mw。由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：各條支路的首端功率p各條支路的首端功率q各條支路的末端功率p各條支路的末端功率qs(5,1)0.790.52-0.78-0.49s(6,2)1.381.01-1.37-0.85s(7,3)1.180.81-1.18-0.73s(8,4)0.790.54-0.78-0.49s(5,9)-0.79-0.520.79-0.11s(6,9)-0.19-1.080.200.86s(6,7)-1.190.071.20-0.20s(7,9)0.66-1.30-0.641.13s(7,8)-3.030

24、.703.10-1.05s(8,10)-3.880.514.00-0.873.8 斷開環(huán)網(wǎng)回線的潮流的分布（1）分析斷開6、9支路時(shí)的潮流分析當(dāng)斷開6、9支路時(shí)，導(dǎo)致支路數(shù)變?yōu)?程序矩陣改變，其他參數(shù)不變。用matlab軟件進(jìn)行潮流計(jì)算進(jìn)行調(diào)整。節(jié)點(diǎn)號12345678910調(diào)整前1.0180.9160.9710.9911.0400.9751.0061.0211.0501.050調(diào)整后0.9920.9160.9710.9651.0400.9751.0061.0211.0501.050 由上表觀察到，進(jìn)行最后一次調(diào)節(jié)后節(jié)點(diǎn)1,2,3,4的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后

25、各個(gè)支路的有功損耗，記錄于下表： 功率損耗pqs(5,1)0.0010.04s(6,2)0.0060.18s(7,3)0.0030.09s(8,4)0.002 0.05s(5,9)0.0050.01s(6,9)0.0000.00s(6,7)0.0250.06s(7,9)0.0380.09s(7,8)0.0590.14總損耗0.2560.93具體的潮流分布如下表：各條支路的首端功率p各條支路的首端功率q各條支路的末端功率p各條支路的末端功率qs(5,1)0.800.53-0.80-0.50s(6,2)1.411.05-1.40-0.87s(7,3)1.200.83 -1.20-0.74s(8,4

26、)0.800.55-0.80-0.50s(5,9)-0.80-0.530.81-0.10s(6,9)0.00 0.000.000.00s(6,7)-1.41-1.051.430.96s(7,9)0.39-1.92-0.351.79s(7,8)-3.020.123.08-0.47s(8,10)-3.88-0.084.00-0.27 分析斷開6、7支路時(shí)的潮流分析各條支路的首端功率p各條支路的首端功率q各條支路的末端功率p各條支路的末端功率qs(5,1) 0.800.53-0.80-0.50s(6,2)1.411.04-1.40-0.87s(7,3)1.200.83-1.20-0.74s(8,4)

27、0.80 0.55-0.80-0.50s(5,9)-0.80-0.530.81 -0.10s(6,9)-1.41-1.041.440.87s(6,7)0.000.00 0.000.0s(7,9)1.81 -1.62-1.751.54 s(7,8)-3.020.793.08-1.14s(8,10)-3.880.604.00-0.95分析斷開7、9支路時(shí)的潮流分析各條支路的首端功率p各條支路的首端功率q各條支路的末端功率p各條支路的末端功率qs(5,1)0.800.53-0.80-0.50s(6,2)1.411.04 -1.40-0.87s(7,3)1.200.83 -1.20-0.74s(8,4

28、)0.800.55-0.80 -0.50s(5,9)-0.80-0.53 0.81 -0.10s(6,9)0.38-1.81-0.341.65s(6,7)-1.790.771.82-0.85s(7,9)0.000.000.000.00s(7,8)-3.020.013.08-0.36s(8,10)-3.88-0.194.00-0.16 四、總結(jié)從本次電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)中，通過計(jì)算在不同的負(fù)荷下各支路節(jié)點(diǎn)潮流分布的變化，分析負(fù)荷變動(dòng)時(shí)，對功率的影響，我們掌握了常用的調(diào)整電壓的措施。在本次課程設(shè)計(jì)過程當(dāng)中，鍛煉了自己實(shí)際操作分析能力，理論聯(lián)系實(shí)際，對運(yùn)行中的電力系統(tǒng)，通過潮流計(jì)算可以預(yù)知各種負(fù)

29、荷變化和網(wǎng)絡(luò)結(jié)構(gòu)的改變會(huì)不會(huì)危及系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許的范圍以內(nèi)，系統(tǒng)中各種元件(線路、變壓器等)是否會(huì)出現(xiàn)過負(fù)荷，以及可能出現(xiàn)過負(fù)荷時(shí)應(yīng)事先采取哪些預(yù)防措施等。實(shí)際的情況遠(yuǎn)比我們計(jì)算的情況復(fù)雜，這讓我深刻了解了潮流計(jì)算的重要性。精準(zhǔn)的潮流計(jì)算不僅可以使電網(wǎng)處于穩(wěn)定且平衡的狀態(tài)，這對電網(wǎng)的安全運(yùn)行起到關(guān)鍵性的作用，還可以實(shí)現(xiàn)運(yùn)行的經(jīng)濟(jì)性。所以我認(rèn)為學(xué)好電力系統(tǒng)的關(guān)鍵在于學(xué)好潮流計(jì)算，電力系統(tǒng)的安全運(yùn)行不是兒戲，一個(gè)很小的失誤就可能造成整個(gè)電網(wǎng)的崩潰，這給國家和社會(huì)帶來的危害和損失之大是可想而知的。學(xué)好電力系統(tǒng)這門課、學(xué)好潮流計(jì)算是我們電專業(yè)學(xué)生的當(dāng)務(wù)之急。沒有堅(jiān)實(shí)的基礎(chǔ)就不

30、可能建造出高樓大廈，所以，要想在以后的工作崗位上有所成就，就必須打好基礎(chǔ)，努力提高自己的專業(yè)素養(yǎng)。 參考文獻(xiàn)1、電力系統(tǒng)分析課程設(shè)計(jì)綜合實(shí)驗(yàn)，祝淑萍 等，中國電力出版社;2、電力系統(tǒng)穩(wěn)態(tài)分析，陳珩，中國電力出版社;3、matlab實(shí)用指南，蘇金明，王永利，電子工業(yè)出版社。4、陳珩.電力系統(tǒng)穩(wěn)態(tài)分析.北京：中國電力出版社，20055、/matpower/附 錄function mpc=caseks5% matpower case format:version 2mpc.version=2;% - power flow data - % sys

31、tem mva basempc.basemva=1;1% bus data % bus_i type pd qd gs bs area vm va basekv zone vmax vminmpc.bus= 1 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94; 2 1 1.4 0.8676 0 0 1 1 0 1 1 1.1 0.94 ; 3 1 1.2 0.7436 0 0 1 1 0 1 1 1.1 0.94; 4 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94 ; 5 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 6 1 0

32、0 0 0 1 1 0 1 1 1.1 0.95; 7 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 8 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 9 3 0 0 0 0 1 1 0 1 1 1.1 0.95; 10 2 0 0 0 0 1 1 0 1 1 1.1 0.95; ; % generator data% bus pg qg qmax qmin vg mbase status pmax pminmpc.gen= 9 0 0 10 -10 1.05 1 1 7 0; 10 4 0 10 -10 1.05 1 1 4 0; ; % branch data

33、% fbus tbus r x b ratea rateb ratec ratio angle status angmin angmax mpc.branch= 5 1 1.543*10-3 4.165*10-2 0 0 0 0 1.025 0 1 -360 360; 6 2 1.84*10-3 5.565*10-2 0 0 0 0 1 0 1 -360 360; 7 3 1.543*10-3 4.165*10-2 0 0 0 0 1 0 1 -360 360; 8 4 1.84*10-3 5.565*10-2 0 0 0 0 1.025 0 1 -360 360; 5 9 0.8034*10

34、-2 1.8998*10-2 0.5882 0 0 0 0 0 1 -360 360; 6 9 1.2854*10-2 3.0397*10-2 0.2353 0 0 0 0 0 1 -360 360; 6 7 0.8034*10-2 1.8998*10-2 0.1471 0 0 0 0 0 1 -360 360; 7 9 1.1248*10-2 2.6597*10-2 0.2059 0 0 0 0 0 0 -360 360; 7 8 0.6427*10-2 1.5199*10-2 0.4706 0 0 0 0 0 1 -360 360; 8 10 0.8034*10-2 1.8998*10-2

35、 0.5882 0 0 0 0 0 1 -360 360; ; return; function mpc=caseks2% matpower case format:version 2mpc.version=2;% - power flow data - % system mva basempc.basemva=1;1% bus data % bus_i type pd qd gs bs area vm va basekv zone vmax vminmpc.bus= 1 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94; 2 1 1.4 0.8676 0 0 1 1 0

36、 1 1 1.1 0.94 ; 3 1 1.2 0.7436 0 0 1 1 0 1 1 1.1 0.94; 4 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94 ; 5 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 6 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 7 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 8 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 9 3 0 0 0 0 1 1 0 1 1 1.1 0.95; 10 2 0 0 0 0 1 1 0 1 1 1.1 0.95; ; % generato

37、r data% bus pg qg qmax qmin vg mbase status pmax pminmpc.gen= 9 0 0 10 -10 1.05 1 1 7 0; 10 4 0 10 -10 1.05 1 1 4 0; ; % branch data % fbus tbus r x b ratea rateb ratec ratio angle status angmin angmax mpc.branch= 5 1 1.543*10-3 4.165*10-2 0 0 0 0 1.025 0 1 -360 360; 6 2 1.84*10-3 5.565*10-2 0 0 0 0

38、 1 0 1 -360 360; 7 3 1.543*10-3 4.165*10-2 0 0 0 0 1 0 1 -360 360; 8 4 1.84*10-3 5.565*10-2 0 0 0 0 - 0 1 -360 360; 5 9 0.8034*10-2 1.8998*10-2 0.5882 0 0 0 0 0 1 -360 360; 6 9 1.2854*10-2 3.0397*10-2 0.2353 0 0 0 0 0 1 -360 360; 6 7 0.8034*10-2 1.8998*10-2 0.1471 0 0 0 0 0 1 -360 360; 7 9 1.1248*10

39、-2 2.6597*10-2 0.2059 0 0 0 0 0 1 -360 360; 7 8 0.6427*10-2 1.5199*10-2 0.4706 0 0 0 0 0 1 -360 360; 8 10 0.8034*10-2 1.8998*10-2 0.5882 0 0 0 0 0 1 -360 360; ; return; function mpc=caseks31% matpower case format:version 2mpc.version=2;% - power flow data - % system mva basempc.basemva=1; % bus data % bus_i type pd qd gs bs area vm va basekv zone vmax vminmpc.bus= 1 1 0.816 0.5057 0 0 1 1 0 1 1 1.1 0.94; 2 1 1.428 0.8850 0 0 1 1 0 1 1 1.1