【正文】 發(fā)電機照成的危害遠大于單相接地故障。x3。simb39。 [t3,x3]=ode45(39。c=()*pi/180。b=。由此試探法 可求得極限切除時間 =。x2]。b=。%給定微分方程組的初值 a=。此外，為了加快功角的收斂速度，添加了阻尼部分，阻尼系數(shù) D 設(shè)為 10。=，代入式 +(176。) 切除故障線路 +(176。 1 1 1 21 1 1 239。 表 運行方式 Z11 α11(176。 各種運行方式下轉(zhuǎn)移阻抗的求解 三相短路故障時，將線路 24 的一回看成是中點直接接地，其余處理方式與穩(wěn)態(tài)時相同。39。等效過程的具體做法入下： 例如對于節(jié)點 2，其負荷阻抗 ZL2=+，與節(jié)點 2 相連的支路對地電納B12=， B24=，則等效導(dǎo)納為： 2 * 1 2 * 2 4 *2*110 . 5 ( ) 0 . 5 ( 0 . 0 5 2 9 2 0 . 0 4 2 3 1 ) 7 . 3 9 5 0 . 9 5 2 40 . 4 2 4 5 0 . 2 3 5 8LY B B j j jZj? ? ? ? ? ? ? ? ? ??即 Z2*=1/Y2*=+。根據(jù)題目要求，發(fā)電機采用二階經(jīng)典模型，即 E 39。 21 4. 轉(zhuǎn)移阻抗的求解 負荷、發(fā)電機、變壓器、線路的處理 負荷的處理 根據(jù)題目的要求，負荷采用恒阻抗模型，其阻抗標幺值的計算公式為： 2* 2BL BLSUZ VS??? （ 8） 其中， U 是負荷所在節(jié)點的電壓， LS? 是負荷的共軛值。于是網(wǎng)損率 η=（ ） /=%，與 Matlab 計算結(jié)果形同。 圖 4. 單項接線圖 采用 power world 仿真軟件的潮流計算結(jié)果 運行模式下，點擊 tools 中的 slove 選項即可對電網(wǎng)進行潮流計算，在 case information中點擊 Power Flow List 可以查看各個節(jié)點的電壓和線路的潮流。仿真器廣泛地使用了圖形和動畫功能，大大地增強了用戶對系統(tǒng)特性、問題和約束的理解，以便于用戶對系統(tǒng)進行維護。表中的功率值均為標幺值，基準功率 SB=100MW。 7）當 和 都滿足要求時，計算結(jié)束。待求解網(wǎng)架中， 1 號節(jié)點為平衡節(jié)點，2~6 號節(jié)點為 PQ 節(jié)點， 7 號節(jié)點（即母線 G）為 PV 節(jié)點。 表 編號 1 2 3 4 5 6 G 節(jié)點類型 平衡節(jié)點 PQ 節(jié)點 PQ 節(jié)點 PQ 節(jié)點 PQ 節(jié)點 PQ 節(jié)點 PV 節(jié)點 有功功率 P/MW － 180 350 260 250 150 837 無功功率 Q/MVar － 100 180 120 120 70 － 電壓 V/kV 230 － － － － － 242 采用 Matlab 計算電網(wǎng)潮流 通過對待求解電網(wǎng)的等效線路參數(shù)的計算與標幺化， 形成節(jié)點導(dǎo)納矩陣。 =1e8。 附注： ODE 函數(shù)說明 Matlab提供了一階常微分方程組求解的系列函數(shù)： ode**。 時三相重合閘。不計內(nèi)阻和空載損耗。變電站參數(shù)和負荷如下表： 編號 類型 220kV 最大負荷， MVA 1 500kV 站 平衡節(jié)點 2 220kV 站 180+j110 3 220kV 站 350+j180 4 220kV 站 260+j120 5 220kV 站 250+j120 6 220kV 站 150+j70 各線路長度如圖所示。 作者簽名： 日期： 年 月 日 學(xué)位論文版權(quán)使用授權(quán)書 本學(xué)位論文作者完全了解學(xué)校有關(guān)保留、使用學(xué)位論文的規(guī)定，同意學(xué)校保留并向國家有關(guān)部門或機構(gòu)送交論文的復(fù)印件和電子版，允許論文被查閱和借閱。電力系統(tǒng)穩(wěn)定分析與計算 課程設(shè)計報告 電力系統(tǒng)穩(wěn)定分析與計算 2 畢業(yè)設(shè)計（論文）原創(chuàng)性聲明和使用授權(quán)說明 原創(chuàng)性聲明 本人鄭重承諾：所呈交的畢業(yè)設(shè)計（論文），是我個人在指導(dǎo)教師的指導(dǎo)下進行的研究工作及取得的成果。本人完全意識到本聲明的法律后果由本人承擔。圖中，各變電站負荷均等效接于 220kV 母線。發(fā)電機組主要參數(shù)如下表： 機組臺數(shù) 單臺容 量（ MW） 額定電壓（ EV） 功率因數(shù) 升壓變?nèi)萘縈VA 升壓變Vs% Xd Xd’ Xq Td0’ TJ=2H 3 300 350 4 11 當發(fā)電機采用三階模型時，勵磁環(huán)節(jié)（含勵磁機和勵磁調(diào)節(jié)器）模型如下（不考慮 PSS） : 18km 12km 12km 12km 25km 16km 15km 1 2 4 5 6 3 G 10km 8 11 RTs? 1AAKTs?Σ Σ Σ1eTs1FFKsTs?KE++VGV1Vr e fVFVRVf+ 上圖中參數(shù)如下： TR=0， KA=25， TA=， Te=， KE=1， KF=， TF= 發(fā)電廠升壓變均采用 Y/Y0 接線，變比 。 用 PowerWorld 軟件進行潮流計算并與自己編制的軟件計算結(jié)果進行校核和分析； 設(shè) 2 和 4 母線之間雙回線路中一回的中點分別發(fā)生以下 2 種故障： 1） 1s 時發(fā)生三相短路， 同時切除故障線路三相； 2） 1s 時單相接地短路， 時同時切除故障線路三相。 三、 設(shè)計要求和設(shè)計成果： 2 位同學(xué)為一組，自行分工，但任務(wù)不能重復(fù)； 每位同學(xué)對自己的設(shè)計任務(wù)編寫課程設(shè)計說明書一份； 一組同學(xué)共同完成一份完整的設(shè)計報告； 9 設(shè)計說明和報告應(yīng)包含： ? 以上設(shè)計任務(wù)每一部分的計算過程和結(jié)果分析； ? 所編制的潮流和穩(wěn)定計算源程序（主要語句應(yīng)加注釋）； ? 潮流計算結(jié)果（潮流圖） ? 穩(wěn)定計算的功角曲線等； ? 網(wǎng)絡(luò)變換法求解轉(zhuǎn)移阻抗的變換過程圖。 Ode45函數(shù)引用如下： x0= [ 1, 0] options=odeset。此外，節(jié)點 6 為 PV 節(jié)點，由于發(fā)電機并不總是滿載運行，在正常運行時，考慮到 7%的廠用電，因而發(fā)電機機組發(fā)出的總有功功率為： 滿載運行時， 3 (1 7 % ) 3 0 0 8 3 7P M W? ? ? ? ? （ 6） 于是 ，當發(fā)電機機組滿載時，每臺機組的出力為 837/3=279MW，發(fā)電機機端額定電壓為 ，出線側(cè)的高壓母線電壓穩(wěn)定在 ，各個節(jié)點的參數(shù)信息如表 3 所示。 15 構(gòu)造節(jié)點信息矩陣 B2，題目給定的網(wǎng)絡(luò)共有 7 個節(jié)點，本文編程中的 B2 矩陣為 7 6 的矩陣，每一行代表一個節(jié)點的信息，每一行中的各元素的含義： 該節(jié)點的發(fā)電機功率（注入）； 該節(jié)點的負荷功率（輸出）； 空行（保持矩陣格式）； PV 節(jié)點電壓給定值和 PQ 節(jié)點的初始值； 各節(jié)點電壓相角； 節(jié)點類型（其中， 1 表示平衡節(jié)點， 2 表示 PQ 節(jié)點， 3 表示 PV 節(jié)點）。若 ΔQ 大于可接受誤差，則進行電壓幅值修正。 2）各支路的狀態(tài)如表 。此外，系統(tǒng)模型可以通過使用仿真軟件的圖形編輯工具很容易地進行修改，用戶只需輕輕點擊幾下鼠標就可以在檢修期間切換線路、增加新的線路或發(fā)電機、確定新的交易容量。 18 圖 3. 發(fā)電機的參數(shù)設(shè)置 采用 power world 畫出的單項接線圖如圖 。即公式（ 7） = 1 0 0 %?DLD發(fā) 電 機 輸 出 功 率 P 負 荷 P網(wǎng) 損 率 發(fā) 電 機 輸 出 功 率 P （ 7） 其中， PD=發(fā)電機發(fā)出的有功功率 +平衡節(jié)點與外界交換的有功功率， PL=負荷的有功率代數(shù)和， PD=3 279+=，總負荷 PL=180+350+260+250+150=1190MW。說明兩種計算方法都是可行的。=。 輸電線路的處理 在輸電線路的 π 值等效中，線路兩端存在并聯(lián)電容，它們與負荷阻抗構(gòu)成并聯(lián)關(guān)系，在網(wǎng)絡(luò)變換時首先應(yīng)將這些并聯(lián)電容和負荷阻抗等效為一個等值阻抗。47， Z39。 E ’ V （ 即 節(jié) 點 1 ）Z AZ B Z C 圖 其中， ZA=+， ZB=+， ZC=+。由此計算得的各量的值如表 。 ( )E V V V?? ? ? ? （ 13） 239。) 三相短路故障 +(176。m=，穩(wěn)定時 PE39。 三相短路故障 二階模型的轉(zhuǎn)子運動方程： 0( ( ) ( 1 ) )NeJd wdtwdw P P D wd t T??? ??????? ? ? ? ??? （ 15） 由于 Pe 的表達式在不同的狀態(tài)下（故障前、故障后、切除后）是不一樣的，所以需要 Matlab 定義多個時間段進行仿真。w0]。%給第二個計算區(qū)間賦初值 a=。 x=[x1?？芍敼收锨谐龝r間為 tc=，時發(fā)電機暫態(tài)穩(wěn)定；當故障切除時間為 tc=，時發(fā)電機暫態(tài)不穩(wěn)定，為了更明顯的看出發(fā)電機暫態(tài)過程的失穩(wěn)現(xiàn)象， tc= 時仿真時間改為 15s。 a=。b=。c=()*pi/180。 [t4,x4]=ode45(39。x2。由此可以判斷發(fā)電機在線路24 經(jīng)歷 發(fā)生單相接地故障， 切除故障線路， 重合閘， 再次斷開等一系列操作后依然能夠保持暫態(tài)穩(wěn)定。由此可得出發(fā)電機采用二階模型時系統(tǒng)的極限切除時 間 =。三臺機組則定義為經(jīng)典 三 階模型，即 E39。 37 6. 發(fā)電機模型選擇對穩(wěn)定計算結(jié)果的影響 由 節(jié)和 節(jié)可知， 發(fā)電機分別采用二階模型和三階模型時， 極限切除時間 分別為 t2=， t3=，相差 Δ=，二者相比，極限切除時間相差并不大 。 % 4 6 + 1 0 0。 % 0 + 0 1 0 2。 %B2= % [0 0 0 0 1。 參考文獻 【 1】 《電力系統(tǒng)分析》（上冊）華中科技大學(xué)出版 ，何仰贊 溫增銀 等 【 2】 《電力系統(tǒng)分析》（上冊）華中科技大學(xué)出版 ，何仰贊 溫增銀 等 【 3】 《發(fā)電廠電氣部分》中國電力出版社 【 4】 Power World 中文版 使用手冊 【 5】 基于 Matlab/Simulink 的系統(tǒng)仿真技術(shù)與應(yīng)用 38 附錄： 附 1. PQ 分解法的 Matlab 程序 % PQ分解法進行潮流計算 % 單位默認電壓（ KV），容量（ MVA），有功（ MW），無功（ Mvar） ,電流（ kA）； % 數(shù)據(jù)輸入格式 % 支路參數(shù)矩陣 B1(列 )：變壓器支路化為阻抗串聯(lián)變比 1： K形式 ,且需化成阻抗側(cè)變比為 1側(cè) % 支路首端節(jié)點號； 末端節(jié)點號； 支路阻抗 (變壓器支路為等效阻抗）； % 支路對地電納（變壓器支路該參數(shù)為 0）； 支路的變比 K（無變壓器支路變比為 1）； % 支路首端處于 K側(cè)為 1，處于 1側(cè)為 0； 是否為變壓器支路（是為 1，不是為 0） % 節(jié)點參數(shù)矩陣 B2(列 )：輸入時第 n行即為節(jié)點編號 n % 該節(jié)點發(fā)電機功率（注入）； 該節(jié)點負荷功率（輸出）； N_R法各節(jié)點電壓初始值（直角坐標） /PQ分解法時為空行 % PV節(jié)點 V的幅值 /非 PV節(jié)點 (N_R法為 0， PQ法為初始值 )； PQ法各節(jié)點電壓相角 /N_R法為空行； % 節(jié)點分類標號： 1為平衡節(jié)點； 2為 PQ節(jié)點； 3為 PV節(jié)點； %本例中的矩陣為： %B1= % [1 2 + 1 0 0。 圖 24. 線路 發(fā)生三相短路故障， 切除時發(fā)電機的暫態(tài)圖像 3）用試探法可找出發(fā)電機的極限切除時間， 線路 24 中點發(fā)生三相對稱短路故障， 切除故障（圖 25 左）時系統(tǒng)穩(wěn)定， 切除故障（圖 25 右）時系統(tǒng)不穩(wěn)定，如圖 。 表 16. Power World 和 Matlab 仿真結(jié)果對比 極限切除時間 是否失穩(wěn) 曲線是否一致 Power World Matlab Power World Matlab Power World Matlab 三相短路 穩(wěn)定 穩(wěn)定 一致 一致 單相接地故障 — — 穩(wěn)定 穩(wěn)定 一致 一致 相對誤差（ s） — — — — 絕對誤差（ %） — — — — 由上表可以看出， Matlab 和 Power World 的仿真結(jié)果基本相同，短路故