【正文】 注入有功，節(jié)點(diǎn)注入無功 各節(jié)點(diǎn)的支路首端注入有功率，支路首端注入無功率，支路末端注入有功率，支路末端注入無功率，支路有功損耗，支路無功損耗 。 0 50+ 0 0 2。 0 40+ 10 0 2。 B2=[0 0 231 0 1。] B2 矩陣： 該節(jié)點(diǎn)發(fā)電機(jī)功率； 該節(jié)點(diǎn)負(fù)荷功率； PQ 節(jié)點(diǎn)電壓初始值，或平衡節(jié)點(diǎn)及 PV節(jié)點(diǎn)電壓的給定值。 2 3 + 0 220/ 1 1。 1 4 + 1 0 0。并且，將節(jié)點(diǎn) 1 設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn) 6設(shè)為 PV 節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為 PQ 節(jié)點(diǎn)。 ⑹計(jì)算各節(jié)點(diǎn)電壓的新值，即修正后值 ? ? ? ? ? ? ? ? ? ? ? ?001001 。 ⑵給這各節(jié)點(diǎn)電壓初始值 )0()0( ,fe 。潮流計(jì)算時(shí)一般只設(shè)一個(gè)平衡節(jié)點(diǎn)。等值負(fù)荷和等值電源的有功功率 LiP 、 GiP 是給定的，從而注入有功功率 iP 是給定的。0 ?? GiGi QP 對(duì)狀態(tài)變量 iU 的約束條件則是 m axm in iii UUU ?? 對(duì)某些狀態(tài)變量 i? 還有如下的約束條件 m axjiji ???? ??? 節(jié)點(diǎn)的分類： ⑴ 第一類稱 PQ節(jié)點(diǎn)。其中 , 1U 、 2U 主要受 1GQ 、 2GQ 的控制 , 1? 、 2? 主要受 1GP 、 2GP 的控制。常規(guī)的潮流計(jì)算一般給定其中的二個(gè)變量： PQ 節(jié)點(diǎn)（注入有功功率及無功功率），PV節(jié)點(diǎn)（注入有功功率及電壓的大?。胶夤?jié)點(diǎn)（電壓的大小及相角）。潮流計(jì)算是實(shí)現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)供電的必要手段和重要工作環(huán)節(jié)，因此潮流計(jì)算在電力系統(tǒng)的規(guī)劃設(shè)計(jì)、生產(chǎn)運(yùn)行、調(diào) 度管理及科學(xué)研究中都有著廣泛的應(yīng)用。為了保證電力系統(tǒng)的穩(wěn)定運(yùn)行，要進(jìn)行潮流調(diào)節(jié)。 4. 對(duì)不同的負(fù)荷變化，進(jìn)行潮流的調(diào)節(jié)控制，并說明調(diào)節(jié)控制的方法，并列表表示調(diào)節(jié)控制的參數(shù)變化。（變電所低壓母線電壓 10KV要求調(diào)整范圍在 之間；電壓 35KV要求調(diào)整范圍在 3536 之間） 5. 輪流斷開 支路雙回線中的一條 ，分析潮流的分布。 二、 課程設(shè)計(jì)基本內(nèi)容： 1. 對(duì)給定 的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖 。 課 程 設(shè) 計(jì) 報(bào) 告 學(xué)生姓名 : 學(xué)號(hào)： 學(xué) 院 : 班 級(jí) : 題 目 : 電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì) 指導(dǎo)教師： 年 月 日 課設(shè)題目及要求 一 .題目原始資料 系統(tǒng)圖： 兩個(gè)發(fā)電廠分別通過變壓器和輸電線路與四個(gè)變電所相連。 線路長(zhǎng)為 80km 。 2. 輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行分析。 （幾條支路斷幾次） 6. 利用 DDRTS 軟件，進(jìn)行繪制系統(tǒng)圖進(jìn)行上述各種情況潮流的分析，并進(jìn)行結(jié)果的比較。 5. 打印利用 DDRTS 進(jìn)行潮流分析繪制的系統(tǒng)圖，以及潮流分布圖。 隨著電力系統(tǒng)及在線應(yīng)用的發(fā)展，計(jì)算機(jī)網(wǎng)絡(luò)已經(jīng)形成，為電力系統(tǒng)的潮流計(jì)算提供了物質(zhì)基礎(chǔ)。它的發(fā)展主要圍繞這樣幾個(gè)方面：計(jì)算方法的收斂性、可靠性；計(jì)算速度的快速性；對(duì)計(jì)算機(jī)存儲(chǔ)容量的要求以及計(jì)算的方便、靈活等。 變量的分類： 負(fù)荷消耗的有功、無功功率 —— 1LP 、 1LQ 、 2LP 、 2LQ 電源發(fā)出的有功、無功功率 —— 1GP 、 1GQ 、 2GP 、 2GQ 母線或節(jié)點(diǎn)的電壓大小和相位 —— 1U 、 2U 、 1? 、 2? 在這十二個(gè)變量中，負(fù)荷消耗的有功和無功功率無法控制，因它們?nèi)Q于用戶，它們就稱為不可控變量或是擾動(dòng)變量。這四個(gè)變量就是簡(jiǎn)單系統(tǒng)的狀態(tài)變量。等值負(fù)荷功率 LiP 、 LiQ 和等值電源功率 GiP 、 GiQ 是給定的，從而注入功率 iP 、 iQ 是給定的，待求的則是節(jié)點(diǎn)電壓的大小 iU 和相位角 i? 。等值負(fù)荷的無功功率 LiQ 和節(jié)點(diǎn)電壓的大小 iU 也是給定的。等值負(fù)荷功率 LsP 、 LsQ是給定的，節(jié)點(diǎn)電壓的大小和相位也是給定的。 ⑶將以上電壓初始值代入式（ 4— 38a） ~式（ 4— 38c）或式（ 4— 45c）、（ 4— 45a），求出修正方程式中的不平衡量 )0(2)0()0( )(, iii UQP ??? 以及 。 iiiiii fffeee ?????? ⑺運(yùn)用各節(jié)點(diǎn)電壓的新值自第三步開始進(jìn)入下一次迭代。 二 .參數(shù)求取 設(shè)定基準(zhǔn)值 MVASB 100? ， KVUB 220? ，所以 4842 ?? BBB USZ? 根據(jù)題目原始資料，計(jì)算發(fā)電廠、變壓器及線路的參數(shù)。 1 7 + 1 0 0。 4 5 + 0 220/11 1 1。 節(jié)點(diǎn)所接無功補(bǔ)償并聯(lián)電容（感）的電納。 0 0 220 0 2。 60 0 231 0 3 0 0 220 0 2。] (1) 編寫程序見附錄 1 (2) 運(yùn)行并調(diào)試 表 2：調(diào)試變電所電壓要求 變電所 1 2 3 4 電壓范圍 35~36 ~ 35~36 ~ 初始矩陣運(yùn)行 結(jié)果 ： : 表 3：各節(jié)點(diǎn)電壓 節(jié)點(diǎn) 1 2 3 4 5 電壓 231 節(jié)點(diǎn) 6 7 8 9 10 電壓 231 表 4：各節(jié)點(diǎn)電壓相角 節(jié)點(diǎn) 1 2 3 4 5 相角 0 節(jié)點(diǎn) 6 7 8 9 10 相角 可以觀察到節(jié)點(diǎn) 10 的電壓都與題目給定的范圍相比偏高，因此調(diào)節(jié)發(fā)電 廠電壓和變壓器分接頭，通過對(duì)系統(tǒng)進(jìn)行多次調(diào)整，最終得到合理結(jié)果 。 DDRTS 仿真潮流計(jì)算結(jié)果圖（系統(tǒng)圖見附錄 2） DDRTS 圖中的 顯示的數(shù)據(jù)為各支路潮流分布，節(jié)點(diǎn)電 壓和節(jié)點(diǎn) 電壓 相角，與MATLAB 對(duì)比潮流計(jì)算結(jié)果相同。 0 0 220 0 2。 0 0 220 0 2。 0 0 220 0 2。 60 0 231 0 3。 0 49+ 10 0 2] 表 13：調(diào)節(jié)方式 調(diào)節(jié)方法 廠 1 電壓 廠 2 電壓 分接頭 1 分接頭 2 分接頭 3 分接頭 4 未調(diào)整 231 231 220/ 220/11 220/ 220/11 調(diào)整后 223 225 表 14： MATLAB 與 DDRTS 仿真的 調(diào)節(jié)結(jié)果 節(jié)點(diǎn)號(hào) 1 2 3 4 5 未調(diào)整 調(diào)整后 DDRTS 節(jié)點(diǎn)號(hào) 6 7 8 9 10 未調(diào)整 調(diào)整后 DDRTS 表 15：各 支路功率損耗 功率損耗 未調(diào)整 調(diào)整后 1~2 支路 1~4 支路 1~7 支路 6~9 支路 7~9 支路 2~3 支路 + + 4~5 支路 + + 7~8 支路 + + 9~10 支路 + + 總損耗 由表可以觀察到負(fù)荷上升后功率損耗有所減少 表 16：調(diào)節(jié)后潮流分布 各條支路的首端功率 Si 各條支路的末端功率 Sj 1~2 支路 + 1~4 支路 + 1~7 支路 6~9 支路 7~9 支路 2~3 支路 + 4~5 支路 + 7~8 支路 + 9~10 支路 + DDRTS 仿真結(jié)果圖 （ 3） 4 變電所的負(fù)荷同時(shí)以 2%的比例 下降， 2， 3 變電所的負(fù)荷同時(shí)以 2%比例上升 。 0 0 220 0 2。 0 + 35 0 2。 1 7 + 1 0 0。 4 5 + 0 220/11 1 1。 1 4 + 1 0 0。 2 3 + 0 220/ 1 1。] 下面進(jìn)行調(diào)節(jié)，使電壓達(dá)到要求范圍： 表 25：調(diào)節(jié)方式 調(diào)節(jié)方法 廠 1 電壓 廠 2 電壓 分接頭 1 分接頭 2 分接頭 3 分接頭 4 未調(diào)整 231 231 220/ 220/11 220/ 220/11 調(diào)整后 223 225 表 26： MATLAB 與 DDRTS 仿真的 調(diào)節(jié)結(jié)果 節(jié)點(diǎn)號(hào) 1 2 3 4 5 未調(diào)整 調(diào)整后 DDRTS 節(jié)點(diǎn)號(hào) 6 7 8 9 10 未調(diào)整 調(diào)整后 DDRTS 表 27：各支路功 率損耗 功率損耗 未調(diào)整 調(diào)整后 1~2 支路 1~4 支路 1~7 支路 6~9 支路 7~9 支路 2~3 支路 + + 4~5 支路 + + 7~8 支路 + + 9~10 支路 + + 總損耗 表 28：調(diào)節(jié)后潮流分布 各條支路的首端功率 Si 各條支路的末端功率 Sj 1~2 支路 + 1~4 支路 + 1~7 支路 6~9 支路 7~9 支路 2~3 支路 + 4~5 支路 + 7~8 支路 + 9~10 支路 + DDRTS 仿真結(jié)果圖 （ 6） 斷開 1， 7支路 雙回線中一條 B1矩陣中對(duì)應(yīng)的阻抗和導(dǎo)納發(fā)生了變化 ： B1=[1 2 + 1 0 0。 7 9 + 1 0 0。 9 10 + 0 220/11 1 1。 6 9 17+ 1 0 0。 7 8 + 0 220/ 1 1。 1 7 + 1 0 0。 4 5 + 0 220/11 1 1。 本次電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)中，通過計(jì)算在不同的負(fù)荷下各支路節(jié)點(diǎn)潮流分布的變化，分析負(fù)荷變動(dòng)時(shí)，對(duì)電壓，功率的影響，我們掌握了常用的調(diào)整電壓的措施。 這次課程設(shè)計(jì)主要是給定網(wǎng)絡(luò)接線圖，和各個(gè)變電所的負(fù)荷以及對(duì)個(gè)別節(jié)點(diǎn)的電壓和功率的要求，需要我們畫出等效電路圖，在計(jì)算出各元件參 數(shù)的基礎(chǔ)上，運(yùn)用牛頓 拉夫遜法，通過調(diào)節(jié)機(jī)端電壓，變壓器變比和無功補(bǔ)償來調(diào)節(jié)電壓到要求范圍內(nèi)。 實(shí)際的情況遠(yuǎn)比我們計(jì)算的情況復(fù)雜，這讓我深刻了解了潮流計(jì)算的重要性。只有扎實(shí)的基礎(chǔ)才能在將來的工作當(dāng)中不出現(xiàn)問題 。 附錄一 附錄二 clear。)。)。)。)。 6 9 + 1 0 0。 7 8 + 0 220/ 1 1。)。 0 0 220 0 2。 0 0 220 0 2。V=zeros(1,n)。q=B1(i,2)。 %非對(duì)角元 Y(q,p)=Y(p,q)。q=B1(i,2)。 %對(duì)角元 j 側(cè) +線路電納的一半 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./。disp(Y) %給定各節(jié)點(diǎn)初始電壓及給定各節(jié)點(diǎn)注入功率 G=real(Y)。 V(i)=abs(B2(i,3))。Q=imag(S)。N1=N0+1。 %求取各個(gè)節(jié)點(diǎn)的功率及功率偏差及 PV節(jié)點(diǎn)的電壓偏差 for i=1:n %n 個(gè)節(jié)點(diǎn) 2n 行 (每節(jié)點(diǎn)兩個(gè)方程 P 和 Q 或 U) p=2*i1。 for j1=1:n %第 i 行共 n 列 (n個(gè)節(jié)點(diǎn)間互導(dǎo)納及節(jié)點(diǎn)電壓相乘即電流 ) C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。%節(jié)點(diǎn)功率 Q 計(jì)算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej)