【正文】 0 1 5 + 1 0 0 1 7 + 1 0 0 2 9 + 1 0 0 3 4 + 0 231/ 1 1 5 6 + 0 231/11 1 1 7 8 + 0 231/ 1 1 7 9 + 1 0 0 9 10 + 0 231/11 1 1 用修改后的 B2 矩陣帶入程序，用 MATLAB 進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整， 對(duì)應(yīng)的調(diào)節(jié)方法和結(jié)果見下表： 調(diào)節(jié)方法： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭四 未調(diào)整 235 235 231 231 231 231 調(diào)整后 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 調(diào)整后 由上表觀察到， 調(diào)節(jié)后 的電壓均在題目允許的范圍內(nèi)，對(duì)線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表： 功率損耗 未調(diào)整 調(diào)整后 (1,3)支路 B1= (1,5)支路 (1,7)支路 (2,9)支路 (3,4)支路 + + (5,6)支路 + + (7,8)支路 + + (7,9)支路 (9,10)支路 + + 總損耗 由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。 將 PSAT仿真得到的潮流圖中個(gè)點(diǎn)電壓的有名值與 MATLAB計(jì)算得到的各節(jié)點(diǎn)電壓進(jìn)行比較： 電壓 1 2 3 4 5 6 7 8 9 10 DDTRS MATLAB 由上述比較可觀察到 : 應(yīng)用 DDTRS仿真和 MATLAB計(jì)算的到的結(jié)果大致相同，可以認(rèn)為此潮流計(jì)算是正確的。具體的潮流分布如下： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)= S(9,2)=+ S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)=+ S(9,7)= S(9,10)=+ S(10,9)= 各節(jié)點(diǎn)的功率 S為 (節(jié)點(diǎn)號(hào)從小到大排列 )： + + + + 4 號(hào)變電所負(fù)荷成 2%下降， 3 號(hào)變電負(fù)荷成 2%上升 此時(shí)的各負(fù)荷變?yōu)椋? 變電所 1 LS =+ 變電所 2 LS = 51+ 變電所 3 LS =+ 變電所 4 LS =+ 相應(yīng)的， B2 矩陣變?yōu)椋? 0 0 235 0 1 200 0 235 0 3 0 0 220 0 2 0 + 35 0 2 0 0 220 0 2 0 51+ 10 0 2 0 0 220 0 2 0 + 35 0 2 0 0 220 0 2 0 + 10 0 2 將程序中的 B2 矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，此系統(tǒng)經(jīng)過調(diào)節(jié)后達(dá)到要求，對(duì)應(yīng)的調(diào)節(jié)方法和結(jié)果見下表： 調(diào)節(jié)方法： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭四 未調(diào)整 235 235 231 231 231 231 調(diào)整后 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 調(diào)整后 由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對(duì)線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表： 功率損耗 未調(diào)整 調(diào)整后 B2= (1,3)支路 (1,5)支路 (1,7)支路 (2,9)支路 (3,4)支路 + + (5,6)支路 + + (7,8)支路 + + (7,9)支路 (9,10)支路 + + 總損耗 由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。 變電所負(fù)荷變化時(shí)的潮流計(jì)算 4 個(gè)變電所負(fù)荷同時(shí)以 2%的比例增大 由于各個(gè)變電所的負(fù)荷同時(shí)以 2%的比例增大，所以此時(shí)的各負(fù)荷變?yōu)椋? 變電所 1 LS =+ 變電所 2 LS = 51+ 變電所 3 LS =+ 變電所 4 LS =+ 相應(yīng)的， B2 矩陣變?yōu)椋? 0 0 235 0 1 200 0 235 0 3 0 0 220 0 2 0 + 35 0 2 0 0 220 0 2 0 51+ 10 0 2 0 0 220 0 2 0 + 35 0 2 0 0 220 0 2 0 + 10 0 2 將原程序中的 B2 矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，經(jīng)過調(diào)節(jié)后達(dá)到要求，對(duì)應(yīng)的調(diào)節(jié)方法和結(jié)果見下表： 調(diào)節(jié)方法： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭四 未調(diào)整 235 235 231 231 231 231 調(diào)整后 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 調(diào)整后 由上表觀察到， 調(diào)節(jié)后 的電壓均在題目允許的范圍內(nèi)，對(duì)線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表： 功率損耗 未調(diào)整 調(diào)整后 (1,3)支路 (1,5)支路 (1,7)支路 (2,9)支路 (3,4)支路 + + (5,6)支路 + + (7,8)支路 + + (7,9)支路 (9,10)支路 + + 總損耗 由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。 且得到各節(jié)點(diǎn)的功率 (節(jié)點(diǎn)號(hào)從小到大排列 )： + + + + 仿真比較 由 MATLAB 得到電壓、功角和具體的潮流分布圖，調(diào)整好以后的數(shù)據(jù)如圖： 在 DDTRS 中繪制出電路圖，仿真并進(jìn)行調(diào)節(jié)以后可以得到潮流分布圖： 在用 DDTRS 繪制的系統(tǒng)潮流分布圖中，各個(gè)節(jié)點(diǎn)電壓在調(diào)試之后均能滿足要求，均達(dá)到系統(tǒng)電壓等級(jí)標(biāo)準(zhǔn)。 調(diào)節(jié)前后各支路損耗對(duì)比如下： 功率損耗 未調(diào)整 調(diào)整后 (1,3)支路 (1,5)支路 (1,7)支路 (2,9)支路 (3,4)支路 + + (5,6)支路 + + (7,8)支路 + + (7,9)支路 (9,10)支路 + + 總損耗 由此可以看出，隨著 發(fā)電機(jī)端電壓的降低，網(wǎng)絡(luò) 有 功損耗逐漸增大，符合實(shí)際情況，而最終調(diào)節(jié)使得電壓在規(guī)定范圍內(nèi)。 根據(jù)題目的不同要求對(duì)參數(shù)進(jìn)行調(diào)整，通過調(diào)節(jié)變壓器變比和發(fā)電廠的電壓，求解出合理的潮流分布，最后用 DDTRS 進(jìn)行潮流分析，將兩者進(jìn)行比較。 %)51(2201 ????k 1123111 %)51(2201 ????k 2 %)(2 2 02 ????k 11 2 511 %)(2 2 02 ????k 12203 ???k 1122020 12203 ???k 1 %)(2 2 04 ????k 11 1 411 %)(2 2 04 ????k %)51(2205 ????k 1120911 %)51(2205 ????k 根據(jù)題目給定的參數(shù)和要求，手繪 形成了 網(wǎng)絡(luò)的 等值電路圖 ,并 且在圖中標(biāo)注出了各參數(shù)和節(jié)點(diǎn)編號(hào)。 參數(shù)計(jì)算 根 據(jù)題目原始資料，計(jì)算發(fā)電廠、變壓器及線路的 有名值 參數(shù)。并且，將節(jié)點(diǎn) 1設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn) 2設(shè)為 PV 節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為 PQ 節(jié)點(diǎn)。 7. 最終形成課程設(shè)計(jì)成品說明書。（變電所低壓母線電壓 10KV 要求調(diào)整范圍在 之間；電壓 35KV 要求調(diào)整范圍在 3536 之間） 5. 輪流斷開支路雙回線中的一條，分析潮流的分布。 3. 跟隨變電所負(fù)荷按一定比例發(fā)生變化，進(jìn)行潮流計(jì)算分析。 1. 對(duì)給定 的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖。 發(fā)電廠資料： 母線 1 和 2 為發(fā)電廠高壓母線，發(fā)電 廠一總裝機(jī)容量為（ 300MW ）， 母線 3為 機(jī)壓母線，機(jī)壓母線上裝機(jī)容量（ 100MW），最大負(fù)荷和最小負(fù)荷為 50MW和 20MW；發(fā)電廠二總裝機(jī)容量（ 200MW）。如果不收斂，轉(zhuǎn)回②進(jìn)行下次迭代計(jì)算，直到收斂為止。 2． 牛頓法的框圖及求解過程 （ 1）用牛頓法計(jì)算潮流時(shí)，有以下的步驟： ①給這各節(jié)點(diǎn)電壓初始值 )0()0( ,fe ； ② 將 以 上 電 壓 初 始 值 代 入 公 式 ， 求 修 正 方 程 的 常 數(shù) 項(xiàng) 向 量 )0(2)0()0( )(, VQP ??? ； ③將電壓初始值在帶入上述公式，求出修正方程中系數(shù)矩陣的各元素。按上述得到的修正向量 nXXX ?????? , 21 ? 后，得到如下關(guān)系 nnn XXX ?????? 這比 nXXX ??? , 21 ? 更接近真實(shí)值。 用同樣的方法考慮，給出 n 個(gè)變量的 n 個(gè)方程： ??????????0),(0),(0),(21212211nnnnXXXfXXXfXXXf???????? 對(duì)其近似解 1X? 得修正量 1X?? 可以通過解下邊的方程來確定： ????????????????????????????????????????????????????????????????????????????????nnnnnnnnnnnXXXxfxfxfxfxfxfxfxfxfXXXfXXXfXXXf?????????2121222121211121212211),(),(),( 式中等號(hào)右邊的矩陣nnxf?? 都是對(duì)于 nXXX ??? , 21 ? 的值。 這一方法還可以做下面的解釋，設(shè)第 n 次迭代得到的解語真值之差，即 )(nX的誤差為 ? 時(shí)，則： 0)( )( ???nXf 把 )( )( ??nXf 在 )(nX 附近對(duì) ? 用泰勒級(jí)數(shù)展開 0. . . . . .)(!2)()()( )(2)()()( ????????? nnnn XfXfXfXf ??? 上式省略去 2? 以后部分 0)()( )()( ??? nn XfXf ? )(nX 的誤差可以近似由上式計(jì)算出來。 — 拉夫遜法概要 牛頓 — 拉夫遜法迭代原理 已知一個(gè)變量 X函數(shù)為： 0)( ?Xf ， 由適當(dāng)?shù)慕浦?)0(X 出發(fā)，根據(jù)： , . . . . . . )2,1()( )( )( )()()1( ????? nXf XfXX nnnn 反復(fù)進(jìn)行計(jì)算，當(dāng) )(nX 滿足適當(dāng)?shù)氖諗織l件就是上面方程的根。等值負(fù)荷功率