【正文】 3) 1 和 4 號變電所的負荷同時以 2%的比例下降，而 2 和 3 號變電所的負荷同時以 2%的比例上升； 4. 在不同的負荷情況下，分析潮流計算的結(jié)果，如果各母線電壓不滿足要求，進行電壓的調(diào)整。 （ 1）運用下列公式計算變壓器 低壓側(cè) 參數(shù)： SUPR NNkT 221000? SUUX NNKT 100 2%? XRZ TTT ?? （ 2）計算線路參數(shù) LjxrjXRZ )( ???? 將參數(shù)整理成如下表格： 首端號 末端號 阻抗有名值 電納有名值 1 3 + 1 5 + 1 7 + 2 9 + 3 4 + 0 5 6 + 0 7 8 + 0 7 9 + 9 10 + 0 （ 3）變電所負荷分別為： 變電所 1 LS =60+ 變電所 2 LS =50+ 變電所 3 LS =40+ 變電所 4 LS =60+ (4）計算變壓器分接頭變比 變壓器有 5個抽頭，電壓調(diào)節(jié)范圍為 NU ? 2*%， NU 對應(yīng)的分接頭開始時設(shè)變壓器高壓側(cè)接主接頭 ,降壓變壓器 5個分接頭時變比 K以 備調(diào)壓時選用 。 將 PSAT仿真得到的潮流圖中個點電壓的有名值與 MATLAB計算得到的各節(jié)點電壓進行比較： 電壓 1 2 3 4 5 6 7 8 9 10 DDTRS MATLAB 由上述比較可觀察到 : 應(yīng)用 DDTRS仿真和 MATLAB計算的到的結(jié)果大致相同，可以認為此潮流計算是正確的。具體的潮流分布如下： 各條支路的首端功率 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(1,3)=+ S(10,9)= 各節(jié)點的功率 S為 (節(jié)點號從小到大排列 )： + + + + + 斷開 7 支路的一條線 當(dāng)斷開 7節(jié)點間的雙回線中的一條時 時， B1矩陣 的改變， B1 矩陣變 1 3 + 1 0 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 B1= 7 9 + 1 0 0 9 10 + 0 231/11 1 1 用修改后的 B2 矩陣帶入程序，用 MATLAB 進行潮流計算并進行調(diào)整， 對應(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)，對線路損耗進行分析，統(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)支路 + + 總損耗 由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認為是合理的。對規(guī)劃中的電力系統(tǒng)，通過潮流計算可以檢驗所提出的電力系統(tǒng)規(guī)劃方案能否滿足各種運行方式的要求；對運行中的電力系統(tǒng)，通過潮流計算可以預(yù)知各種負荷變化和網(wǎng)絡(luò)結(jié)構(gòu)的改變會不會危及系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許的范圍以內(nèi)，系統(tǒng)中各種元件 (線路、變壓器等 )是否會出現(xiàn)過負荷，以及可能出現(xiàn)過負荷時應(yīng)事先采取哪些預(yù)防措施等。由于該調(diào)壓措施的調(diào)壓范圍比較窄及上述數(shù)據(jù)分析可知，各變電所的電壓均偏高，若降低機端電壓，在此時的負荷條件下可滿足個別變電所的電壓要求，但出廠母線電壓均降到 220KV 以下，是系統(tǒng)運行在 低于額定電壓等級狀況下 ,在系統(tǒng)的負荷發(fā)生改變時 ,在某些幾點處的電壓會驟然下降 ,引起系統(tǒng)的癱瘓 ,造成大面積停電 . 現(xiàn)代電力系統(tǒng)網(wǎng)絡(luò)中所常用的變壓器多為有載變壓器，能在有載條件下進行調(diào)節(jié)，不需退出運行就能改變變比。 使我在潮流計算、 MATLAB 的使用方面均有所提高，但也暴露出了一些問題；理論知識的儲備還是不足，對MATLAB 的性能和特點還不能有一個全面的把握，對 MATLAB 的 使用不夠熟悉，相信通過以后的學(xué)習(xí)能彌補這些不足，從而達到一個新的層次。 參考文獻： 《電力系統(tǒng)分析基礎(chǔ) 》 ，韋剛， 中國電力出版社 。 nl=9。 pr=。 3 4 + 0 231/ 1 1。 0 60+ 35 0 2。)。 SB=100。 %UB=input(39。 end disp(39。 end disp(39。 %非對角元 Y(q,p)=Y(p,q)。 %分解出導(dǎo)納陣的實部和虛部 for i=1:n %給定各節(jié)點初始電壓的實部和虛部 e(i)=real(B2(i,3))。IT2=1。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。 %節(jié)點有功功率差 DQ=Q(i)Q1。q=2*j11。J(m,q)=X2。q=2*j11。% 擴展列△ P J(m,q)=X2。X6=0。q=q+1。J(p,N)=DV。j1~=i %非平衡節(jié)點 amp。 J(m,q)=X1。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 J(m,q)=X1。 % dP/de=dQ/df X2=B(i,j1)*e(i)G(i,j1)*f(i)。% 節(jié)點功率 Q 計算 fi Σ(Gij*ejBij*fj)eiΣ (Gij*fj+Bij*ej) %求 P39。 while IT2~=0 IT2=0。 %i 節(jié)點注入功率 SGSL B(i,i)=B(i,i)+B2(i,5)。導(dǎo)納矩陣 Y=39。q=B1(i,2)。 B2(i,2)=B2(i,2)./SB。 BB1=B1。)。V=zeros(1,n)。 0 40+ 35 0 2。 9 10 + 0 231/11 1 1。 B1=[1 3 + 1 0 0。 isb=1。 n=10。但是由于考慮的不夠周全，沒有使兩個發(fā)電廠均衡的承擔(dān)負荷，使發(fā)電廠一的功率過小。對于一個網(wǎng)絡(luò)的電壓調(diào)整，常在變電所添加調(diào)相機、電容器等設(shè)備。當(dāng)系統(tǒng)的電壓偏離允許值時，電力系統(tǒng)必須應(yīng)用電壓調(diào)節(jié)技術(shù)調(diào)節(jié)系統(tǒng)電壓的大小，使其維持在允許值范圍內(nèi)。 將 PSAT仿真得到的潮流圖中個點電壓的有名值與 MATLAB計算得到的各節(jié)點電壓進行比較 。具體的潮流分布如下： 各條支路的首端功率 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é)點的功率 S為 (節(jié)點號從小到大排列 )： + + + + 仿真比較 負荷每次發(fā)生變化時相應(yīng)的修改 DDTRS 中電路的參數(shù)，仿真并進行調(diào)節(jié)，可以得到相應(yīng)的潮流分布圖，現(xiàn)于由 MATLAB 得到的潮流分布數(shù)據(jù)圖比較（以 4 個變電所負荷增加 2%為例）： 由 MATLAB 得到的數(shù)據(jù)圖： 由 DDTRS 得到的潮流分布圖： 在用 DDTRS 繪制的系統(tǒng)潮流分布圖中，各個節(jié)點電壓在調(diào)試之后均能滿足要求，均達到系統(tǒng)電壓等級標準。 五、潮流計算 給定負荷下的潮流計算 B1/B2 矩陣的形成 根據(jù)所求參數(shù)，以及 B B2 矩陣的含義，列寫 B B2 矩陣如下： 1 3 + 1 0 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 0 0 235 0 1 200 0 235 0 3 0 0 220 0 2 0 60+ 35 0 2 0 0 220 0 2 0 50+ 10 0 2 0 0 220 0 2 0 40+ 35 0 2 0 0 220 0 2 0 60+ 10 0 2 潮流調(diào)整及分析 將編好的的程序 (程序見附錄 )輸入 MATLAB 中運行，可以得到如下運行結(jié)果 : 節(jié)點 1 2 3 4 5 電壓 節(jié)點 6 7 8 9 10 電壓 由表觀察到節(jié)點 4,6,8,10 的電壓值都較正常范 圍偏高，因此調(diào)節(jié)發(fā)電廠電壓B1= B2= 和變壓器分接頭，經(jīng)過多次調(diào)整，使節(jié)點電壓均滿足要求，調(diào)節(jié)結(jié)果如下： 調(diào)整過程： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭四 未調(diào)整 235 235 231 231 231 231 第一次 230 230 231 231 231 231 第二次 225 225 231 231 231 231 第三次 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 第一次 第二次 第三次 由上表觀察到，先后進行發(fā)電機機端電壓調(diào)節(jié)和變壓器分接頭的調(diào)節(jié)使得節(jié)點 4,6,8,10 的節(jié)點電壓均在允許范圍內(nèi)。 四、初步分析 根據(jù)系統(tǒng)圖可知此系統(tǒng)為兩端供電網(wǎng)路， 為了潮流計算的完整性， 將母線 1，2設(shè)為節(jié)點 1， 2，將變 電所 4的高低壓側(cè)分別設(shè)為節(jié)點 10。 電廠一 電廠二 輸電線路資料： 發(fā)電廠和變電所之間的輸電線路的電壓等級及長度標于圖中，單位長度的電阻為 ? ，單位長度的電抗為 ? ，單位長度的電納為 S610* 。 ④解修正方程式 )0()0( , fe ?? ； ⑤修正各節(jié)點電壓 )0()0()1( eee ??? ， )0()0()1( fff ??? ； ⑥將 )1(e ， )1(f 在帶入方程式，求出 )1(2)1()1( )(, VQP ??? ； ⑦檢驗是否收斂，即 ? ? ???? )()( ,m a x kiki QP 如果收斂，迭代到此結(jié)束，進一步計算各線路潮流和平衡節(jié)點功率，并打印 輸出結(jié)果。這樣的方法就是所謂的牛頓 — 拉夫遜法。 第二類稱 PV節(jié)點：等值負荷和等值電源的有功功率 PLi、 PGi是給定的，從而注入有功功率 Pi是給定的。它的發(fā)展主要圍繞這樣幾個方面：計算方法的收斂性、可靠性；計算速度的快速性；對計算機存儲容量的要求以及計算的方便、靈活等。此外，電力系統(tǒng)的潮流計算也是計算機系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)，所以