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