【正文】 的存儲(chǔ)矩陣和各節(jié)點(diǎn)注入功率（實(shí)部和虛部）的存儲(chǔ)矩陣，以存放迭代過(guò)程中的中間值和輸出最終值；4）構(gòu)造各節(jié)點(diǎn)的潮流方程，進(jìn)入潮流計(jì)算程序；5）基于MATLAB平臺(tái)，得出計(jì)算結(jié)果，采用誤差精度為10e5。 PV節(jié)點(diǎn)電壓V的給定值；節(jié)點(diǎn)所接的無(wú)功補(bǔ)償設(shè)備的容量?；贛ATLAB平臺(tái)編寫(xiě)程序，我們對(duì)待求解電網(wǎng)進(jìn)行潮流計(jì)算，從而得出各PQ節(jié)點(diǎn)的節(jié)點(diǎn)電壓以及各線(xiàn)路的線(xiàn)路損耗。此外，節(jié)點(diǎn)6為PV節(jié)點(diǎn)，其中總發(fā)出的有功功率（5）。L12的線(xiàn)路長(zhǎng)度為27km，則Z12=*（+）*27/484Ω=(+) ΩB12=2*2**50**106*27*484S=L23的線(xiàn)路長(zhǎng)度為17km，則Z23=*（+）*17/484Ω=(+) ΩB23=2*2**50**106*17*484S=L24的線(xiàn)路長(zhǎng)度為6km，則Z24=*（+）*6/484Ω=(+) ΩB24=2*2**50**106*6*484S=L45的線(xiàn)路長(zhǎng)度為5km，則Z45=*（+）*5/484Ω=(+) ΩB45=2*2**50**106*5*484S=L36的線(xiàn)路長(zhǎng)度為18km，則Z36=*（+）*18/484Ω=(+) ΩB36=2*2**50**106*18*484S=L56的線(xiàn)路長(zhǎng)度為24km，則Z56=*（+）*24/484Ω=(+) ΩB56=2*2**50**106*24*484S=綜上，電網(wǎng)線(xiàn)路參數(shù)的標(biāo)幺值如表1所示：表線(xiàn)路參數(shù)的標(biāo)幺值統(tǒng)計(jì)線(xiàn)路L12L23L24L45L36L56阻抗Z/Ω++++++導(dǎo)納B/S 節(jié)點(diǎn)信息統(tǒng)計(jì)表節(jié)點(diǎn)給定參數(shù)統(tǒng)計(jì)編號(hào)123456節(jié)點(diǎn)類(lèi)型平衡節(jié)點(diǎn)PQ節(jié)點(diǎn)PQ節(jié)點(diǎn)PQ節(jié)點(diǎn)PQ節(jié)點(diǎn)PV節(jié)點(diǎn)P(MW)無(wú)250380280420Q(MVar)無(wú)40908595無(wú)V(kV)230無(wú)無(wú)無(wú)無(wú)由圖1可知，網(wǎng)絡(luò)中共有6個(gè)節(jié)點(diǎn)，其中節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，視為無(wú)窮大母線(xiàn)，電壓穩(wěn)定在230kV。 網(wǎng)絡(luò)變換法求解轉(zhuǎn)移阻抗的變換過(guò)程圖。 潮流計(jì)算結(jié)果（潮流圖）178。 以上設(shè)計(jì)任務(wù)每一部分的計(jì)算過(guò)程和結(jié)果分析；178。比較兩種模型的仿真結(jié)果，分析發(fā)電機(jī)模型選擇對(duì)于穩(wěn)定計(jì)算結(jié)果的影響，并且分析勵(lì)磁調(diào)節(jié)系統(tǒng)參數(shù)變化對(duì)于穩(wěn)定計(jì)算結(jié)果的影響；編制課程設(shè)計(jì)報(bào)告。與PowerWorld軟件的分析結(jié)果進(jìn)行比較校核。二、 設(shè)計(jì)的主要內(nèi)容：進(jìn)行參數(shù)計(jì)算和標(biāo)幺化，形成潮流計(jì)算參數(shù)；用Matlab編制潮流計(jì)算程序，要求采用PQ分解法潮流計(jì)算方法。穩(wěn)定仿真中不考慮發(fā)電廠的調(diào)速器和原動(dòng)機(jī)模型?？紤]該電廠開(kāi)機(jī)三臺(tái)，所有發(fā)電機(jī)保留10%的功率裕度。不計(jì)內(nèi)阻和空載損耗。發(fā)電機(jī)組主要參數(shù)如下表：機(jī)組臺(tái)數(shù)單臺(tái)容量（MW）額定電壓（EV）功率因數(shù)升壓變?nèi)萘縈VAXdXd’XqTd0’TJ=2H330035087當(dāng)發(fā)電機(jī)采用三階模型時(shí)，勵(lì)磁環(huán)節(jié)（含勵(lì)磁機(jī)和勵(lì)磁調(diào)節(jié)器）模型如下:圖中參數(shù)如下：TR=0，KA=20，TA=，Te=，KE=1，SE=0，KF=，TF=不考慮PSS，即圖中的附加信號(hào)Vs=0。燃煤發(fā)電廠G有三臺(tái)機(jī)組，均采用單元接線(xiàn)。F/km；40186。所有線(xiàn)路型號(hào)均為L(zhǎng)GJ2*300，基本電氣參數(shù)為：正序參數(shù)： r = , x = , C = 181。電力系統(tǒng)穩(wěn)定分析和計(jì)算設(shè)計(jì)報(bào)告電力系統(tǒng)穩(wěn)定分析和計(jì)算設(shè)計(jì)報(bào)告專(zhuān)業(yè)：電氣工程及其自動(dòng)化班級(jí)： 09電氣（4）班學(xué)號(hào)：32號(hào)、34號(hào)姓名：禤培正、陳子程華南理工大學(xué)電力學(xué)院2012229目錄課程設(shè)計(jì)題目：電力系統(tǒng)穩(wěn)定分析和計(jì)算 11 潮流的參數(shù)計(jì)算和標(biāo)幺化 4 線(xiàn)路參數(shù)計(jì)算與標(biāo)幺化 4 節(jié)點(diǎn)信息統(tǒng)計(jì) 52 采用MATLAB計(jì)算電網(wǎng)潮流 6 PQ分解法潮流計(jì)算 6 MATLAB平臺(tái)潮流計(jì)算結(jié)果分析 73 采用powerworld仿真軟件計(jì)算潮流 7 PowerWorld仿真軟件簡(jiǎn)介 7 PowerWorld軟件建模與潮流計(jì)算 8 PowerWorld仿真結(jié)果與Matlab計(jì)算結(jié)果對(duì)比分析 114 網(wǎng)絡(luò)變換法求解輸入阻抗與轉(zhuǎn)移阻抗 12 負(fù)荷、發(fā)電機(jī)、變壓器、線(xiàn)路等效電容的處理 12 采用網(wǎng)絡(luò)變換法的具體步驟 15 求解輸入阻抗和轉(zhuǎn)移阻抗 19 求解有功功率傳輸特性 195 故障方案仿真 20 選定故障1：節(jié)點(diǎn)3發(fā)生三相對(duì)稱(chēng)短路 20 選定故障2：線(xiàn)路L23中點(diǎn)處發(fā)生三相對(duì)稱(chēng)短路 266 采用等面積法求解極限切除角 31 故障方案一 31 故障方案二 327 基于PowerWorld，分析三階模型的發(fā)電機(jī)的暫態(tài)過(guò)程 33 故障方案一：節(jié)點(diǎn)3發(fā)生三相對(duì)稱(chēng)短路 33 故障方案二：線(xiàn)路L23中點(diǎn)處發(fā)生三相對(duì)稱(chēng)短路 358 發(fā)電機(jī)模型及勵(lì)磁調(diào)節(jié)參數(shù)對(duì)穩(wěn)定計(jì)算結(jié)果的影響 36 發(fā)電機(jī)模型對(duì)穩(wěn)定計(jì)算結(jié)果的影響 36 TJ對(duì)穩(wěn)定計(jì)算結(jié)果的影響 38 勵(lì)磁參數(shù)對(duì)穩(wěn)定計(jì)算結(jié)果的影響 39參考文獻(xiàn) 40附錄： 40[鍵入文字]課程設(shè)計(jì)題目：電力系統(tǒng)穩(wěn)定分析和計(jì)算姓名：指導(dǎo)教師：一、 一個(gè)220kV分網(wǎng)結(jié)構(gòu)和參數(shù)如下：500kV站（1）的220kV母線(xiàn)視為無(wú)窮大母線(xiàn)，電壓恒定在230kV。圖中，各變電站參數(shù)如下表：編號(hào)類(lèi)型220kV最大負(fù)荷，MVA1500kV站平衡節(jié)點(diǎn)2220kV站250+j403220kV站380+j904220kV站280+j855220kV站420+j95各線(xiàn)路長(zhǎng)度如圖所示。F/km；零序參數(shù)： r0 = , x0 = , C0 = 181。C長(zhǎng)期運(yùn)行允許的最大電流：1190A。電廠220kV側(cè)采用雙母接線(xiàn)。發(fā)電廠升壓變參數(shù)均為Vs%=%。發(fā)電廠按PV方式運(yùn)行。發(fā)電廠廠用電均按出力的7%考慮。負(fù)荷采用恒阻抗模型。用PowerWorld軟件進(jìn)行潮流計(jì)算并與自己編制的軟件計(jì)算結(jié)果進(jìn)行校核和分析；用Matlab編制穩(wěn)定計(jì)算程序（三臺(tái)機(jī)可并聯(lián)等值成一臺(tái)機(jī)），發(fā)電機(jī)采用二階經(jīng)典模型，要求給出網(wǎng)絡(luò)變換法求解輸入阻抗和轉(zhuǎn)移阻抗的變換過(guò)程圖；自行選擇23種故障方案，給出搖擺曲線(xiàn)，并計(jì)算故障的極限切除時(shí)間和極限切除角。對(duì)上面23種故障方案，采用等面積定則計(jì)算故障的極限切除角，并分析比較與時(shí)域仿真法計(jì)算結(jié)果的差別；用Matlab編制穩(wěn)定計(jì)算程序，發(fā)電機(jī)采用三階模型，并對(duì)第6步的23種故障方案進(jìn)行穩(wěn)定計(jì)算，給出搖擺曲線(xiàn)，并計(jì)算故障的極限切除時(shí)間；與PowerWorld軟件的分析結(jié)果進(jìn)行比較校核。三、 設(shè)計(jì)要求和設(shè)計(jì)成果：2位同學(xué)為一組，自行分工，但任務(wù)不能重復(fù)；每位同學(xué)對(duì)自己的設(shè)計(jì)任務(wù)編寫(xiě)課程設(shè)計(jì)說(shuō)明書(shū)一份；一組同學(xué)共同完成一份完整的設(shè)計(jì)報(bào)告；設(shè)計(jì)說(shuō)明和報(bào)告應(yīng)包含：178。 所編制的潮流和穩(wěn)定計(jì)算源程序（主要語(yǔ)句應(yīng)加注釋?zhuān)?78。 穩(wěn)定計(jì)算的功角曲線(xiàn)等；178。1 潮流的參數(shù)計(jì)算和標(biāo)幺化 線(xiàn)路參數(shù)計(jì)算與標(biāo)幺化圖電網(wǎng)的等效線(xiàn)路圖圖1是待求解電網(wǎng)的等效線(xiàn)路圖，由已知取基準(zhǔn)功率100MW，基準(zhǔn)電壓220kV：則有 線(xiàn)路正序參數(shù)： 線(xiàn)路阻抗標(biāo)幺值計(jì)算公式： 線(xiàn)路導(dǎo)納標(biāo)幺值計(jì)算公式： 由于電網(wǎng)線(xiàn)路為雙回路線(xiàn)路，即阻抗值為單回路的一半，導(dǎo)納值為單回路的兩倍。節(jié)點(diǎn)5均為PQ節(jié)點(diǎn)，帶有一定的負(fù)荷。綜上，各節(jié)點(diǎn)的參數(shù)如表2所示：2 采用MATLAB計(jì)算電網(wǎng)潮流通過(guò)對(duì)待求解電網(wǎng)的等效線(xiàn)路參數(shù)的計(jì)算與標(biāo)幺化，得出待求解電網(wǎng)的各PQ節(jié)點(diǎn)的有功功率和無(wú)功功率，各PV節(jié)點(diǎn)的有功功率和節(jié)點(diǎn)電壓。我們采用PQ分解法對(duì)待求解電網(wǎng)進(jìn)行潮流計(jì)算，大致流程如下：1）編程思想是將所有支路的狀態(tài)信息整合成一個(gè)矩陣進(jìn)行計(jì)算、將所有節(jié)點(diǎn)的狀態(tài)信息整合成另一個(gè)矩陣進(jìn)行計(jì)算，具體如下：構(gòu)造支路信息矩陣B1，每一行表示一條支路，一行中的各元素表示該支路的信息：支路首端號(hào)；末端號(hào)；支路阻抗；支路對(duì)地電納 支路的變比；支路首端處于K側(cè)為1，首端處于1側(cè)則為0；構(gòu)造節(jié)點(diǎn)信息矩陣B2，每一行表示一個(gè)節(jié)點(diǎn)，一行中的各元素表示該節(jié)點(diǎn)的信息：該節(jié)點(diǎn)發(fā)電機(jī)功率；該節(jié)點(diǎn)負(fù)荷功率；節(jié)點(diǎn)電壓初始值。 節(jié)點(diǎn)分類(lèi)標(biāo)號(hào)：1為平衡節(jié)點(diǎn)（應(yīng)為1號(hào)節(jié)點(diǎn)）；2為PQ節(jié)點(diǎn)；3為PV節(jié)點(diǎn)；根據(jù)前文對(duì)待求解電網(wǎng)的線(xiàn)路參數(shù)的歸算和標(biāo)幺化，所構(gòu)造的支路信息矩陣B1和構(gòu)造節(jié)點(diǎn)信息矩陣B2如下所示（MATLAB界面截圖表示），圖MATLAB界面截圖2）接下來(lái)，根據(jù)支路信息矩陣B1和節(jié)點(diǎn)信息矩陣B2，程序先形成待求解電網(wǎng)的節(jié)點(diǎn)導(dǎo)納矩陣：先對(duì)空矩陣加入對(duì)角線(xiàn)元素，即先加入每個(gè)節(jié)點(diǎn)的對(duì)地導(dǎo)納y0；根據(jù)支路信息矩陣B1，對(duì)導(dǎo)納矩陣再加入互導(dǎo)納ypq；再將互導(dǎo)納ypq加入到導(dǎo)納矩陣的對(duì)角線(xiàn)元素，形成自導(dǎo)納。由結(jié)果可知，共用7次迭代便能滿(mǎn)足精度要求，； MATLAB平臺(tái)潮流計(jì)算結(jié)果分析采用MATLAB平臺(tái)計(jì)算電網(wǎng)潮流的輸出結(jié)果如下所示，（1）各節(jié)點(diǎn)電壓如表3所示：表各節(jié)點(diǎn)電壓節(jié)點(diǎn)編號(hào)基準(zhǔn)電壓kV標(biāo)幺電壓實(shí)際電壓(kV)相角 (度)1220022203220