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