電力系統(tǒng)潮流上機(jī)設(shè)計(jì)教學(xué)指導(dǎo)書061210(存儲(chǔ)版)
【正文】 用simulink搭建仿真模型等方法，進(jìn)行在變電所的某一負(fù)荷情況下的潮流計(jì)算及仿真，并對(duì)計(jì)算結(jié)果進(jìn)行分析。 進(jìn)行潮流計(jì)算圖2為仿真模型圖。由表6的電壓調(diào)節(jié)前后功率損耗對(duì)比，可以看出有功功率隨著變壓器分接頭變比的增大而逐漸增大，使得變壓器的低壓側(cè)的電壓處于允許范圍內(nèi)，符合其要求。e=zeros（1，n）；Y(p,q)=Y(p,q)。V(i)=B2(i3)。whilemax(abs(De))ipamp。V=abs(v)。type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%該子程序是用來求取Df for i=1:nIfi=phC(i)=0。Df(p)=P(i)P1。endfunction [De,Df]=hxf(J,Df,ph,n,no)%該子函數(shù)是為求取De Df DX=JDf。[x,y]=size(DX)。case 3 P=2*i1。X2=G(i,j)*e(i)+B(i,j)*f(i)。q=2*j1。J(m,q)=X4。X3= C(i)+G(i,j)*e(i)B(i,j)*f(i)。end J(m,q)=X3。End2實(shí)例驗(yàn)證 【例題】設(shè)有一系統(tǒng)網(wǎng)絡(luò)結(jié)線見圖1，各支路阻抗和各節(jié)點(diǎn)功率均已以標(biāo)幺值標(biāo)示于圖1中，其中節(jié)點(diǎn)2連接的是發(fā)電廠，設(shè)節(jié)點(diǎn)1電壓保持U1＝，試計(jì)算其中的潮流分布，請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=5 請(qǐng)輸入支路數(shù)：m=7 請(qǐng)輸入平衡母線的節(jié)點(diǎn)號(hào)：ph=l 請(qǐng)輸入支路信息：BI=[ l 2 0．02+0．06i O l 1；1 3 0．08+0．24i 0 1 2；2 3 0．06+0．18i 0 l 3： 2 4 0．06+0．18i O l 4： 2 5 0．04+0．12i 0 l 5： 3 4 0．01+0．03i 0 l 6： 4 5 0．08+0．24i O 1 7] 請(qǐng)輸入節(jié)點(diǎn)信息：B2=[ 0 0 1．06 3；0．2+0．20i 0 1 1；一O．45一O．15i 0 l l；一0．40．05i 0 l 1；一0．6—0．1i 0 1 l] 請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 請(qǐng)輸入修正值：ip=0．000 0l參考文獻(xiàn)[1]陳珩．電力系統(tǒng)穩(wěn)定分析[M]．北京：中國(guó)電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實(shí)用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯(lián)系統(tǒng)全過 程動(dòng)態(tài)電壓研究[J】．中國(guó)電力，2008，4l(10)：l4． SHU Hong—ch吼，SUN Shiyun，et a1．Research on fun prc39。J(:,no)=[]。end J(p,q)=X1。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。J(p,q)=X2。p=2*i1。X1=G(i)*f(i)B(i,j)*e(i)。%該子程序是用來求取jacci矩陣for i=1:n X3=0。End ElseDX=[DX1,0,0]。end end end Df=Df’。V1=e(i)^2+f(i)^2。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。end S(ph)=sum(hh)*v(ph)。[De,Di]=hxf(J,D,F,ph,n,no)。For i=1:m S(i)=B2(i1)B2(i2)。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個(gè)端點(diǎn) %第二列存貯支路的另一個(gè)端點(diǎn) %第三列存貯支路阻抗%第四列存貯支路的對(duì)地導(dǎo)納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號(hào)B2=input(‘n請(qǐng)輸入節(jié)點(diǎn)信息：B2=’)； %第一列為電源側(cè)的功率 %第二列為負(fù)荷側(cè)的功率 %第三列為該點(diǎn)的電壓值%第四列為該點(diǎn)的類型：1為PQ，2為PV節(jié)點(diǎn)，3為平衡節(jié)點(diǎn) A=input（‘n請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗：A=’）； ip=input（‘n請(qǐng)輸入修正值：ip=’）； %ip為修正值)。因此我們可以看出，經(jīng)過調(diào)節(jié)后，節(jié)點(diǎn)9點(diǎn)電壓已經(jīng)滿足了系統(tǒng)的要求。并且，將節(jié)點(diǎn)1設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn)10設(shè)為pv節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為pq節(jié)點(diǎn)。④輸電線路資料：發(fā)電廠和變電所之間的輸電線路的電壓等級(jí)及長(zhǎng)度標(biāo)于圖中，*106s。實(shí)驗(yàn)數(shù)據(jù)如圖所示三母線電力系統(tǒng)中，支路電抗和節(jié)點(diǎn)注入的攻入如圖所示，編寫程序，求個(gè)各條支路的有功潮流分布。三、直流潮流算法計(jì)算步驟1）選擇平衡節(jié)點(diǎn)。237。i,j185。相當(dāng)于注入的有功功率。在一些應(yīng)用場(chǎng)所，如輸電網(wǎng)絡(luò)中，只要考慮的是電力系統(tǒng)中有功功率的分布，而不需要計(jì)算各個(gè)節(jié)點(diǎn)電壓幅值，且要計(jì)算速度要快，這勢(shì)必要對(duì)潮流計(jì)算進(jìn)行簡(jiǎn)化處理，本節(jié)實(shí)驗(yàn)就是研究直流潮流計(jì)算，編程與調(diào)試，獲得電力系統(tǒng)中各支路的有功分布，為進(jìn)一步進(jìn)行電力系統(tǒng)分析作準(zhǔn)備。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。endk=k+1。L(i,i)=L(i,i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))endendN(i,i)=N(i,i)2*(U(i))^2*G(i,i)。endx2=[oP,oQ]39。M(i,j)=0。% 賦初值，U為節(jié)點(diǎn)電壓的幅值，a為節(jié)點(diǎn)電壓的相位角U=ones(1,5)。B=imag(Y)。y(2,3)=1/(+)。2 基本步驟和設(shè)計(jì)流程圖形成了雅克比矩陣并建立了修正方程式，運(yùn)用牛頓拉夫遜法計(jì)算潮流的核心問題已經(jīng)解決，已有可能列出基本計(jì)算步驟并編制流程圖。同時(shí)基于MATLAB的計(jì)算機(jī)算法以雙精度類型進(jìn)行數(shù)據(jù)的存儲(chǔ)和運(yùn)算, 數(shù)據(jù)精確度高，能進(jìn)行潮流計(jì)算中的各種矩陣運(yùn)算，使得傳統(tǒng)潮流計(jì)算方法更加優(yōu)化。for(j=k+1。float d,t,x[NBUS]。在潮流計(jì)算中要把高斯消去法子程序放到潮流計(jì)算主程序后面。 n”)。sL[i].X=d2。sB[i].LoadP=d5,sB[i].LoadQ=d6。e=。int i1,i2,i3,kp,kq。float G,B。5,7,0 1,2 2,0 3,0 4,0 5,0 1 1 3 2 3 4 3 4 5 4 2 5 電力系統(tǒng)潮流上機(jī)課程設(shè)計(jì)教學(xué)指導(dǎo)書 1 2 6 2 3 7 2 4 4． 節(jié)點(diǎn)導(dǎo)納矩陣的形成 struct Line { int Num,NumI,NumJ。8． 建議用《電力系統(tǒng)穩(wěn)態(tài)分析》書中的例43 來校驗(yàn)結(jié)果是否正確。學(xué)會(huì)看幫助（F1），在調(diào)試過程中結(jié)束程序運(yùn)行（Ctrl+F2）5． 《電力系統(tǒng)穩(wěn)態(tài)分析》教材中第四章涉及到的基本知識(shí)。各種條件判斷語句，循環(huán)語句的使用。3． Tc數(shù)組或結(jié)構(gòu)體從0開始計(jì)數(shù)。2． 將我提供的程序拷貝到d:lx目錄下，并將Options――Directories中的輸出目錄改為d:lx，就可以運(yùn)行。struct Bus { int Num。float X[NBUS]。struct Bus sB[NBUS]。isB[i].Num=i1。for(i=0。} for(i=0。
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號(hào)-1