基于matlab的電力系統(tǒng)潮流計算課程設(shè)計論文-資料下載頁

　　

【正文】 變變壓器的變比、無功補償設(shè)備的添加等。 對于電力系統(tǒng)，調(diào)節(jié)電壓首先要考慮改變機端電壓，其是最經(jīng)濟最直接的調(diào)壓措施。適用于出廠母線附近沒有其他的用電設(shè)備網(wǎng)絡(luò)，調(diào)壓范圍在 — 之間。由于該調(diào)壓措施的調(diào)壓范圍比較窄及上述數(shù)據(jù)分析可知，各變電所的電壓均偏高，若降低機端電壓，在此時的負荷條件下可滿足個別變電所的電壓要求，但出廠母線電壓均降到 220KV 以下，是系統(tǒng)運行在 低于額定電壓等級狀況下 ,在系統(tǒng)的負荷發(fā)生改變時 ,在某些幾點處的電壓會驟然 下降 ,引起系統(tǒng)的癱瘓 ,造成大面積停電 . 現(xiàn)代電力系統(tǒng)網(wǎng)絡(luò)中所常用的變壓器多為有載變 壓器， 能在有載條件下進行調(diào)節(jié)， 不需 退出運行 就能改變變比。用變壓器調(diào)壓的前提是 必須保證系統(tǒng)的無功功率充足。有上述的數(shù)據(jù)分析， 當調(diào)節(jié)機端電壓不能滿足各個變電所的電壓要求時， 若 個別 變電所電壓偏高， 則升高變壓器的變比；若 個別 變電所電壓偏低，則降低變壓器的變比， 通常可以達到調(diào)壓要求。 無功補償設(shè)備的添加常見的有：串聯(lián)電抗器、并聯(lián)電容器、調(diào)相機等。對于一個網(wǎng)絡(luò)的電壓調(diào)整，常在變電所添加調(diào)相機、電容器等設(shè)備。雖然通過無功補償可以滿足要求，但其成本投入比較大，造價高。 七、心得體會 從本次電力系統(tǒng)潮流計算課程設(shè)計中，通過計算在不同的負荷下各支路節(jié)點潮流分布的變化，分析負荷變動時，對電壓，功率的影響，我們掌握了常用的調(diào)整電壓的措施 ， 鍛煉了自己實際操作分析能力 。 使我在潮流計算、 MATLAB 的使用方面均有所提高，但也暴露出了一些問題；理論知識的儲備還是不足，對MATLAB 的性能和特點還不能有一個全面的把握，對 MATLAB 的 使用不夠熟悉，相信通過以后的學(xué)習(xí)能彌補這些不足，從而達到一個新的層次。 在潮流計算的過程中，為了滿足調(diào)壓要求，需要一次一次的嘗試，首選的調(diào)壓措施是改變機端電壓，一般不能滿足所有的節(jié)點的電壓要求，這就需要調(diào)節(jié)個別的變壓器分接頭，這是兩種較為經(jīng)濟的調(diào)壓措施。 剛開始計算時，沒有充分考慮到負荷分配的問題，雖然電能質(zhì)量能夠滿足，但是導(dǎo)致網(wǎng)絡(luò)損耗過大，網(wǎng)絡(luò)運行不經(jīng)濟；然后重新審視了網(wǎng)絡(luò)結(jié)構(gòu)，對兩個發(fā)電廠所承擔的負荷合理分配以后，使網(wǎng)損較低很多，經(jīng)濟性得到了提高。這個方案充分利用了前 兩種調(diào)壓措施，使電壓都滿足了要求，保證了調(diào)壓的經(jīng)濟性。在 DDRTS仿真時，線路的連接和 參數(shù)的設(shè)置需要 都 細心，有 一點疏忽就會導(dǎo)致運行結(jié)果不合理， 一般 經(jīng)過調(diào)整 就可使得 兩種方法的計算結(jié)果相差甚小。 此次課程設(shè)計，讓我充分了解了解到 精準的潮流計算不僅可以節(jié)省能源并充分發(fā)揮能源的作用，還可以使電網(wǎng)處于穩(wěn)定且平穩(wěn)的狀態(tài)，這對電網(wǎng)的完全運行起到關(guān)鍵性的作用。電力系統(tǒng)的安全運行不是兒戲，一個很小的失誤就可能造成整個電網(wǎng)的崩潰，給國家和社會帶來極大的損失。我們應(yīng)該努力學(xué)好專業(yè)課，打好專業(yè)基礎(chǔ)，在以后的工作崗位上有所成就。 參考文獻： 《電力系統(tǒng)分析基礎(chǔ) 》 ，韋剛， 中國電力出版社 。 《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩，中國電力出版社 。 《 MATLAB 實用指南》，蘇金明，王永利，電子工業(yè)出版社。 附錄 1：潮流分布圖 附錄 2： MATLAB 程序 %本程序的功能是用牛頓 —— 拉夫遜法進行潮流計算 % B1 矩陣： 支路首端號； 末端號； 支路阻抗； 支路對地電納 % 支路的變比； 支路首端處于 K側(cè)為 1， 1側(cè)為 0 % B2 矩陣： 該節(jié)點發(fā)電機功率； 該節(jié)點負荷功率； 節(jié)點電壓初始值 % PV 節(jié)點電壓 V 的給定值； 節(jié)點所接的無功補償設(shè)備的容量 % 節(jié)點分類標號 clear。 n=10。%input(39。請輸入節(jié)點數(shù)： n=39。)。 nl=9。%input(39。請輸入支路數(shù)： nl=39。)。 isb=1。%input(39。請輸入平衡母線節(jié)點號： isb=39。)。 pr=。%input(39。請輸入誤差精度： pr=39。)。 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。 7 9 + 1 0 0。 9 10 + 0 231/11 1 1。 B2=[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] %input(39。請輸入各節(jié)點參數(shù)形成的矩陣： B2=39。)。 Y=zeros(n)。e=zeros(1,n)。f=zeros(1,n)。V=zeros(1,n)。sida=zeros(1,n)。S1=zeros(nl)。 %－－－－－－－修改部分－－－－－－－－－－－－ ym=0。 SB=100。UB=220。 %ym=input(39。您輸入的參數(shù)是標么值？（若不是則輸入一個不為零的數(shù)值）39。)。 if ym~=0 %SB=input(39。請輸入功率基準值 :SB=39。)。 %UB=input(39。請輸入電壓基準值 :UB=39。)。 YB=SB./UB./UB。 BB1=B1。 BB2=B2。 for i=1:nl B1(i,3)=B1(i,3)*YB。 B1(i,4)=B1(i,4)./YB。 end disp(39。B1 矩陣 B1=39。)。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。 B2(i,2)=B2(i,2)./SB。 B2(i,3)=B2(i,3)./UB。 B2(i,4)=B2(i,4)./UB。 B2(i,5)=B2(i,5)./SB。 end disp(39。B2 矩陣 B2=39。)。 disp(B2) end % % % for i=1:nl %支路數(shù) if B1(i,6)==0 %左節(jié)點處于低壓側(cè) p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 %非對角元 Y(q,p)=Y(p,q)。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 %對角元 K側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。 %對角元 1 側(cè) end %求導(dǎo)納矩陣 disp(39。導(dǎo)納矩陣 Y=39。)。 disp(Y) % G=real(Y)。B=imag(Y)。 %分解出導(dǎo)納陣的實部和虛部 for i=1:n %給定各節(jié)點初始電壓的實部和虛部 e(i)=real(B2(i,3))。 f(i)=imag(B2(i,3))。 V(i)=B2(i,4)。 %PV 節(jié)點電壓給定模值 end for i=1:n %給定各節(jié)點注入功率 S(i)=B2(i,1)B2(i,2)。 %i 節(jié)點注入功率 SGSL B(i,i)=B(i,i)+B2(i,5)。 %i 節(jié)點無功補償量 end %=================================================================== P=real(S)。Q=imag(S)。 ICT1=0。IT2=1。N0=2*n。N=N0+1。a=0。 while IT2~=0 IT2=0。a=a+1。 for i=1:n if i~=isb %非平衡節(jié)點 C(i)=0。D(i)=0。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。%Σ(Gij*ej Bij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。%Σ(Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i)。% 節(jié)點功率 P 計算 ei Σ(Gij*ejBij*fj)+fiΣ (Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。% 節(jié)點功率 Q 計算 fi Σ(Gij*ejBij*fj)eiΣ (Gij*fj+Bij*ej) %求 P39。,Q39。 V2=e(i)^2+f(i)^2。 %電壓模平方 %========= 以下針對非 PV 節(jié)點來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV 節(jié)點 DP=P(i)P1。 %節(jié)點有功功率差 DQ=Q(i)Q1。 %節(jié)點無功功率差 %=============== 以 上 為 除 平 衡 節(jié) 點 外 其 它 節(jié) 點 的 功 率 計 算 ================= %================= 求取 Jacobi 矩陣 =================== for j1=1:n if j1~=isbamp。j1~=i %非平衡節(jié)點 amp。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 % dP/de=dQ/df X2=B(i,j1)*e(i)G(i,j1)*f(i)。 % dP/df=dQ/de X3=X2。 % X2=dp/df X3=dQ/de X4=X1。 % X1=dP/de X4=dQ/df p=2*i1。q=2*j11。 J(p,q)=X3。J(p,N)=DQ。m=p+1。 J(m,q)=X1。J(m,N)=DP。q=q+1。 J(p,q)=X4。J(m,q)=X2。 elseif j1==iamp。j1~=isb %非平衡節(jié)點 amp。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 % dQ/de X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。% dQ/df p=2*i1。q=2*j11。J(p,q)=X3。J(p,N)=DQ。%擴展列△ Q m=p+1。 J(m,q)=X1。q=q+1。J(p,q)=X4。J(m,N)=DP。% 擴展列△ P J(m,q)=X2。 end end else %=============== 下面是針對 PV 節(jié)點來求取 Jacobi 矩陣的元素 =========== DP