【文章內(nèi)容簡(jiǎn)介】 ====上面是用線性變換方式將 Jacobi矩陣化成單位矩陣 ===== for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N)。 修改節(jié)點(diǎn) 電壓實(shí)部 k1=k+1。 f(L)=f(L)J(k1,N)。 修改節(jié)點(diǎn)電壓虛部 end 修改節(jié)點(diǎn)電壓 for k=3:N0 DET=abs(J(k,N))。 if DET=pr 電壓偏差量是否滿足要求 IT2=IT2+1。 不滿足要求的節(jié)點(diǎn)數(shù)加 1 end end ICT2(a)=IT2。 ICT1=ICT1+1。 end 16 用高斯消去法解 w=J*V disp(39。迭代次數(shù)： 39。)。 disp(ICT1)。 disp(39。沒有達(dá)到精度要求的個(gè)數(shù)： 39。)。 disp(ICT2)。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 sida(k)=atan(f(k)./e(k))*180./pi。 E(k)=e(k)+f(k)*j。 end %=============== 計(jì)算各輸出量 =========================== disp(39。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值 E為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。)。 disp(E)。 EE=E*UB。 disp(EE)。 disp(39。39。)。 disp(39。各節(jié)點(diǎn)的電壓大小 V為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。)。 disp(V)。 VV=V*UB。 disp(VV)。 disp(39。39。)。 disp(39。各節(jié)點(diǎn)的電壓相角 sida為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。)。 disp(sida)。 for p=1:n C(p)=0。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 end S(p)=E(p)*C(p)。 end 17 disp(39。各節(jié)點(diǎn)的功率 S為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。)。 disp(S)。 disp(39。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。)。 SS=S*SB。 disp(SS)。 disp(39。39。)。 disp(39。各條支路的首端功率 Si為 (順序同您輸入 B1時(shí)一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 if B1(i,6)==0 Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 Siz(i)=Si(p,q)。 else Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 Siz(i)=Si(p,q)。 end disp(Si(p,q))。 SSi(p,q)=Si(p,q)*SB。 ZF=[39。S(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(SSi(p,q))]。 disp(ZF)。 %disp(SSi(p,q))。 disp(39。39。)。 end disp(39。各條支路的末端功率 Sj為 (順序同您輸入 B1時(shí)一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 18 if B1(i,6)==0 Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 Sjy(i)=Sj(q,p)。 else Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 Sjy(i)=Sj(q,p)。 end disp(Sj(q,p))。 SSj(q,p)=Sj(q,p)*SB。 ZF=[39。S(39。,num2str(q),39。,39。,num2str(p),39。)=39。,num2str(SSj(q,p))]。 disp(ZF)。 %disp(SSj(q,p))。 disp(39。39。)。 end disp(39。各條支路的功率損耗 DS為 (順序同您輸入 B1時(shí)一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 DS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。 DDS(i)=DS(i)*SB。 ZF=[39。DS(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(DDS(i))]。 disp(ZF)。 %disp(DDS(i))。 disp(39。39。)。 end figure(1)。 19 subplot(2,2,1)。 plot(V)。 xlabel(39。節(jié)點(diǎn)號(hào) 39。)。ylabel(39。電壓標(biāo)幺值 39。)。 grid on。 subplot(2,2,2)。 plot(sida)。 xlabel(39。節(jié)點(diǎn)號(hào) 39。)。ylabel(39。電壓角度 39。)。 grid on。 subplot(2,2,3)。 bar(S)。 xlabel(39。節(jié)點(diǎn)注入有功 39。)。ylabel(39。節(jié)點(diǎn)注入無功 39。)。 grid on。 subplot(2,2,4)。 bar(Siz)。 xlabel(39。支路首端有功 39。)。ylabel(39。支路首端無功 39。)。 grid on。 —————————————— 程 —— 序 ———————————————— 20 變電所 1 變電所 2 變電所 3 變電所 4 10kV母線 10kV母線 35kV母線 35kV母線 一次側(cè)電壓 220kV 一次側(cè)電壓 220kV 母線 1 母線 3 母線 2 線路長(zhǎng)為 80km 線路長(zhǎng)為 80km 線路長(zhǎng)為 70km 線路長(zhǎng)為 90km 線路長(zhǎng)為 100km 2*QFQ502 3*QFs502 TQN1002 2*TQN1002 。。 線路長(zhǎng)為 90km 三 . 算例 電力系統(tǒng)潮流計(jì)算 課程設(shè)計(jì) 任務(wù)書 一、 題目原始資料 系統(tǒng)圖： 兩個(gè)發(fā)電廠分別通過變壓器和輸電線路與四個(gè)變電所相連。 發(fā)電廠資料： 母線 1和 2為發(fā)電廠高壓母線，發(fā)電 廠一總裝機(jī)容量為 350MW， 母線 3 為 機(jī)壓母線，機(jī)壓母線上裝機(jī)容量為 100MMW， 最大負(fù)荷和最小負(fù)荷分別為 40MW 和20MW；發(fā) 電廠二總裝機(jī)容量為 200WM。 變電所資料： (一 ) 變電所 4低壓母線的電壓等級(jí)分別為： 10KV 10KV 35KV 35KV (二 ) 變電所的負(fù)荷分別為： （ 1） 50MW 50MW 70MW 60MW 電廠一 電廠二 21 （ 2） 60MW 80MW 50MW 45MW （ 3） 50MW 40MW 55MW 70MW （ 4） 70MW 50MW 60MW 60MW （ 5） 70MW 60MW 50MW 45MW （ 6） 40MW 50MW 60MW 70MW （ 7） 70MW 40MW 50MW 60MW （ 8） 30MW 40MW 50MW 60MW （ 9） 40MW 40MW 50MW 60MW （ 10） 50MW 40MW 40MW 60MW (三 ) 每個(gè)變電所的功率因數(shù)均為 cosφ =； (四 ) 變電所 1 和變電所 2 分別配有兩臺(tái)容量為 75MVA 的變壓器，短路損耗 414KW，短路電壓（ %） =；變電所 3 和變電所 4 分別配有兩臺(tái)容量為 63MVA 的變壓器，短路損耗為 245KW，短路電壓（ %） =； 輸電線路資料： 發(fā)電廠和變電所之間的輸電線路的電壓等級(jí)及長(zhǎng)度標(biāo)于圖中，單位長(zhǎng)度的電阻為 ? ，單位長(zhǎng)度的電抗為 ? ，單位長(zhǎng)度的電納為 S610* 。 二、 課程設(shè)計(jì)基本內(nèi)容： 1. 對(duì)給定 的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖 2. 輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行分析。 3. 跟隨變電所負(fù)荷按一定比例發(fā)生變化，進(jìn)行潮流計(jì)算分析。 1) 4 個(gè)變電所的負(fù)荷同時(shí)以 2%的比例增大； 2) 4 個(gè)變電所的負(fù)荷同時(shí)以 2%的比例下降 3) 1 和 4 號(hào)變電所的負(fù)荷同時(shí)以 2%的比例下降，而 2 和 3 號(hào)變電所的負(fù)荷同時(shí)以 2%的比例上升； 4. 在不同的負(fù)荷情況下，分析潮流計(jì)算的結(jié)果，如果各母線電壓不滿足要求，進(jìn)行電壓的調(diào)整。 （變電所低壓母線電壓 10KV 要求調(diào)整范圍在 之間；電壓 35KV 要求調(diào)整范圍在 3536之間） 5. 輪流斷開環(huán)網(wǎng)一回線，分析潮流的分布。 22 6. 利用 DDRTS 軟件，進(jìn)行繪制系統(tǒng)圖進(jìn)行上述各種情況潮流的分析，并進(jìn)行結(jié)果的比較。 7. 最終形成課程設(shè)計(jì)成品說明書。 2. 等值電路的計(jì)算過程 電壓是衡量電力系統(tǒng)電能質(zhì)量的標(biāo)準(zhǔn)之一。電壓過高或過低，都將對(duì)人身及其用電設(shè)備產(chǎn)生重大的影響。保證用戶的電壓接近額定值是電力系統(tǒng)調(diào)度的基本任務(wù)之一。當(dāng)系統(tǒng)的電壓偏離允許值時(shí)，電力系統(tǒng)必須應(yīng)用電壓調(diào)節(jié)技術(shù)調(diào)節(jié)系統(tǒng)電壓的大小，使其維持 在允許值范圍內(nèi)。本文經(jīng)過手算形成了等值電路圖，并編寫好了程序得出節(jié)點(diǎn)電壓標(biāo)幺值，使其滿足所要求的調(diào)整范圍。 ⑴ 系統(tǒng)的等值電路圖 23 我們首先對(duì)給定的程序輸入部分作了簡(jiǎn)要的分析，程序開始需要我們確定輸入節(jié)點(diǎn)數(shù)、支路數(shù)、平衡母線號(hào)、支路參數(shù)矩陣、節(jié)點(diǎn)參數(shù)矩陣。 ⑵ 為了保證整個(gè)系統(tǒng)潮流計(jì)算的完整性，我們把凡具有母線處均選作節(jié)點(diǎn)，這樣，共有 10 條母線即選定 10 個(gè)節(jié)點(diǎn)，我們確定發(fā)電廠一母線為平衡節(jié)點(diǎn)，節(jié)點(diǎn)號(hào)為 ① ，發(fā)電廠二母線為 PV節(jié)點(diǎn)，節(jié)點(diǎn)號(hào)為 ② ，其余節(jié)點(diǎn)均為 PQ節(jié)點(diǎn)，節(jié)點(diǎn)號(hào)見等值電路圖。 ⑶ 確定完節(jié)點(diǎn)及編號(hào)后，各條支路也相應(yīng)確定了，網(wǎng)絡(luò)中總計(jì)有 10 條支路，我們對(duì)各支路參數(shù)進(jìn)行了計(jì)算 。 根據(jù)所給實(shí)際電路圖和題中的已知條件，有以下公式計(jì)算各輸電線路的阻抗和對(duì)地支路電容的標(biāo)幺值和變壓器的阻抗標(biāo)幺值。選擇 電 壓 基 準(zhǔn) 值 為 KVUB 220? 和 功 率 基 準(zhǔn) 值 MVASB 100?，所以SUSYBBBBBB SUZ 322 ??????? ， 。 ⑷ 計(jì)算各線路參數(shù)： 支路阻抗 BZZZjXRZLXLR ????? ?，， 支路對(duì)地電容 BYYYjBYLB ???? ?? ， 1 3 支路 80km 輸電線路： ????? jjXRZ SjLjjBY 46 1022 ?? ????? jjZZZB?????? B??? 24 1 4 支路 80km 輸電線路： ,4 jYSjYjZjZ???????????， 3 4 支路 90km 輸電線路： 12 07 4 jYSjYjZjZ???????????， 4 5 支路 70km 同桿架設(shè)雙回線： 1 8 8