【正文】 障時發(fā)電機暫態(tài)過程中的搖擺曲線基本相同。對于單臺發(fā)電機而言，歸算后的 TJ=， H=， 22( ) 0 . 7 2 2GN Bd T d B G NV SX X X VS? ? ? ? ? ?，22( ) 0 . 1 2 2GN Bd T d B G NV SX X X VS???? ? ? ? ?， 22( ) 0 . 6 1 9GN Bq T q B G NV SX X X VS? ? ? ? ? ? 具體參數(shù)設(shè)計如圖 1 1 13 所示： 圖 11. 平衡機 模型 的參數(shù)設(shè)置 30 圖 12. 發(fā)電機模型 的參數(shù)設(shè)置 圖 13. 發(fā)電機的勵磁機模型參數(shù)設(shè)置 線路 24 中點三相短路時的故障分析 1）在事件窗口添加： 1s 時 線路 24 中點 發(fā)生三相對稱短路，并在 后切除故障，得出的搖擺曲線如 圖 14，此時系統(tǒng)暫態(tài)穩(wěn)定。r39。t2。 x04=x3(end,:)39。,[tc,],x02,h_opt,a,b,c)。simb39。tc=。,t,x(:,2),39。simb39。%短路后的 功率 [t1,x1]=ode45(39。%定義轉(zhuǎn)子運動方程組 end 主程序： w0=1?；疚呛?。m=90176。和發(fā)電機的電磁功率方程如表 15 所示。) 穩(wěn)態(tài)無故障 + + + 三相短路 + + + 切除故障線路 + + + 單相接地故障 + + 功角特性計算 由系統(tǒng)潮流計算可得到平衡節(jié)點 1 的標幺值為： V1*=， P1*=，Q1*=。 四種運行方式經(jīng)網(wǎng)絡(luò)變換后，均可變換為如圖 7 所示的網(wǎng)架結(jié)構(gòu)。7 地 。*+XT* 經(jīng)等效后簡化模型如圖 。、內(nèi)阻為 Xd39。則歸算到全網(wǎng)基準后的發(fā)電機參數(shù)可按式 9 來計算： 2()* 2()GN Bdd G N BV SXX SV??? ? ? （ 9） 帶入?yún)?shù)計算可得： Xd*=， Xq*=， Xd39。說明兩種計算方法都是可行的。 1）各節(jié)點的電壓如表 6 所示。具體參數(shù)設(shè)置如圖 。而且由于其可交互性和可繪 17 圖性，它也可以用于向非專業(yè)用戶解釋電力系統(tǒng)的運行操作。由結(jié)果可知，共用 6 次迭代便能滿足精度要求，共用時間 僅 為 ； PQ 分解法潮流計算的結(jié)果分析 利用 節(jié)的 PQ 分解法的 matlab 程序計算的電網(wǎng)潮流結(jié)果如附錄中附 2 圖所示。 至此，形成節(jié)點導(dǎo)納矩陣 Y。 網(wǎng)絡(luò)信息處理 形成電網(wǎng)線路的節(jié)點導(dǎo)納矩陣是進行潮流計算前必須要做的準備。 系主任 批準 年 月 日 教員組主任 審核 年 月 日 指導(dǎo)老師 簽 發(fā) 20xx 年 12 月 23 日 11 線路參數(shù)的計算和標幺化 220kV 分網(wǎng)結(jié)構(gòu)的 π 型等值電路圖 1 如下圖所示： 1 2 3 4 5 6G~ 圖 待求解電網(wǎng)的基準功率取為 100MVA，基準電壓 220kV，則有 100220484BBBBS MVAV kVZI kA??? ??? ???? ?? （ 1） 線路的正序參數(shù)： j 0 . 0 5 4 j0 . 3 0 8 / k m0 . 0 1 1 6μ F / k mz r xC ? ? ? ? ??? ?? （ 2） 線路阻抗標幺值計算公式： 2BBSZ zl V? ?? （ 3） 線路等值電納標幺值計算公式： 22 BN BVB j f Cl S?? ?? （ 4） 12 當線路為雙回線路時，線路阻抗為單回線線路的一半，導(dǎo)納為單回線路的兩倍。 函數(shù)形式： [t, y]=ode45(odefun,[t1,tf],x0,options)。 用 PowerWorld 軟件的分析問題 4，并與編程計算結(jié)果進行比較校核。穩(wěn)定仿真中不考慮發(fā)電廠的調(diào)速器和原動機模型。F/km； 40186。 作者簽名： 日期： 年 月 日 導(dǎo)師簽名： 日期： 年 月 日 4 注 意 事 項 （論文）的內(nèi)容包括： 1）封面（按教務(wù)處制定的標準封面格式制作） 2）原創(chuàng)性聲明 3）中文摘要（ 300 字左右）、關(guān)鍵詞 4）外文摘要、關(guān)鍵詞 5）目次頁（附件不統(tǒng)一編入） 6）論文主體部分：引言（或緒論）、正文、結(jié)論 7）參考文獻 8）致謝 9）附錄（對論文支持必要時） ：理工類設(shè)計（論文）正文字數(shù)不少于 1 萬字（不包括圖紙、程序清單等），文科類論文正文字數(shù)不少于 萬字。 作 者 簽 名： 日 期： 指導(dǎo)教師簽名： 日 期： 使用授權(quán)說明 本人完全了解 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計（論文）的印刷本和電子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學(xué)?？梢圆捎糜坝?、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目的前提下，學(xué)?？梢怨颊撐牡牟糠只蛉績?nèi)容。 作者簽名： 日 期： 3 學(xué)位論文原創(chuàng)性聲明 本人鄭重聲明：所呈交的論文是本人在導(dǎo)師的指導(dǎo)下獨立進行研究所取得的研究成果。 ：任務(wù)書、開題報告、外文譯文、譯文原文（復(fù)印件）。C 長期運行允許的最大電流： 1190A。負荷采用恒阻抗模型，負序阻抗標幺制取。 發(fā)電機采用三階模型，用 PowerWorld 作為分析工具，對問題 4 的兩種故障方案進行穩(wěn)定計算，給出搖擺曲線，并計算故障的極限切除時間。 說明： odefun：列向量 1*n，通過函數(shù)計算柄輸出的微分方程的右端項； t1,tf: 分別制定積分的時間起點和終點； x0：列向量 1*n，狀態(tài)變量初值 options：：微分優(yōu)化參數(shù)，是一個結(jié)構(gòu)體，使用 odeset可以設(shè)置其具體參數(shù)，詳細內(nèi)容查看幫助。各線路參數(shù)的標幺化計算如下： 線路 L12 長度 15km，則： Z12= 15 （ +） =（ +） Ω B12=2 2π 50 15 106= 104S 標幺值： Z12*=（ +） 100/2202= 104+ B12*= 104 2202/100= 線路 L13 長度 18km，則： Z13= 18 （ +） =（ +） Ω B13=2 2π 50 18 106= 104S 標幺值： Z13*=（ +） 100/2202= 103+ B12*= 104 2202/100= 線路 L24 長度 12km，則： Z24= 12 （ +） =（ +） Ω B24=2 2π 50 12 106= 105S 標幺值： Z24*=（ +） 100/2202= 104+ B24*= 105 2202/100= 線路 L45 長度 10km，則： Z45= 10 （ +） =（ +） Ω B45=2 2π 50 10 106= 105S 標幺值： Z45*=（ +） 100/2202= 104+ B45*= 105 2202/100= 線路 L46 長度 12km，則： Z46=12 （ +） =（ +） Ω B46=2π 50 12 106= 105S 標幺值： Z46*=（ +） 100/2202= 103+ B46*= 105 2202/100= 線路 L3G長度 16km，則： Z3G= 16 （ +） =（ +） Ω B3G=2 2π 50 16 106= 104S 標幺值： Z3G*=（ +） 100/2202= 104+ B3G*= 104 2202/100= 線路 L5G長度 25km，則： Z5G= 25 （ +） =（ +） Ω B5G=2 2π 50 25 106= 104S 標幺值： Z5G*=（ +） 100/2202= 103+ B5G*= 104 2202/100= 13 線路 L6G長度 12km，則： Z6G= 12 （ +） =（ +） Ω B6G=2 2π 50 12 106= 105S 標幺值： Z6G*=（ +） 100/2202= 104+ B6G*= 105 2202/100= 綜上，電網(wǎng)的線路參數(shù)的有名值和標幺值如下表所示： 表 線路名 長度 /km 有名值 標幺值 阻抗 Z/Ω 導(dǎo)納 B/S 阻抗 Z*/Ω 導(dǎo)納 B*/S L12 2 15 + 104 104+ L13 2 18 + 104 103+ L24 2 12 + 105 104+ L45 2 10 + 105 104+ L46 12 + 105 103+ L3G 2 16 + 104 104+ L5G 2 25 + 104 103+ L6G 2 12 + 105 104+ 線路零序參數(shù)的標幺值計算與正序相似，如表 。首先利用線路的標幺值參數(shù)整合成矩陣，使其包含線路的所有信息，再形成節(jié)點導(dǎo)納矩陣。 5） 初始化各節(jié)點電壓（幅值和相角）的存儲矩陣和各節(jié)點注入功率（實部和虛部）的存儲矩陣，以存放迭代過程中的中間值和輸出最終值； 6） 構(gòu)造各節(jié)點的潮流方程，進入潮流計算程序 ， 計算不平衡功率 ΔP， 判斷 ΔP 是否小于可接受誤差。將 matlab 程序計算的結(jié)果整理如下： 1）各節(jié)點的電壓和注入功率如表 。 該仿真器是一個集成的產(chǎn)品，其核心是一個全面、強大的潮流計算程序。 圖 2. 平衡機的參數(shù)設(shè)置 節(jié)點 6 為 PV 節(jié)點， P=837MW， V=，每臺發(fā)電機的輸出功率 P=297MW。 19 表 各節(jié)點的電壓 節(jié)點編號 基準電壓 kV 標幺電壓 實際電壓 (kV) 相角 (度 ) 1 220 0 2 220 3 220 4 220 5 220 6 220 7 220 2）發(fā)電機的狀態(tài)如表 。 2）各支路狀態(tài)的對比如表 10 所示。=。的電壓源來等效原來的三臺發(fā)電機組。 23 1 2 3 4 5 6E ’ 7Z 1 *Z 2 *Z 3 *Z 4 *Z 5 *Z 6 *Z 7 *Z 1 3 *Z 1 2 *Z 2 4 *Z 4 6 *Z 4 5 * Z 6 7 *Z 3 7 *Z 5 7 *X Σ * 圖 5. 等效后簡化模型 網(wǎng)絡(luò)變換（穩(wěn)態(tài)時） 消去節(jié)點 6 利用 YΔ變換，可將節(jié)點 6 消去。將變換后各節(jié)點線路對地的電阻并聯(lián)，即*4* 4 * 4 *1Z= 1 1 1Z Z Z? ??? ??4 地 地，*7* **1Z= 1 1 1Z Z Z? ??? ??7 7 地 7 地，4 7 * *1Z= 11ZZ?? ??47* 47，處理后的網(wǎng)架圖如圖 6 所 示。 25 圖 。由公式（ 11） （ 14）可得到機端電壓 E39。 26 表 運行方式 ΔV δV PE39。(176。 5. 故障仿真 采用 Matlab 軟件分析故障 采用 Matlab 軟件分析故障時，發(fā)電機采用二階模型，即 E39。tc=。simb39。,[tc,40