電力系統(tǒng)課程設(shè)計(jì)報(bào)告

1、東南大學(xué)電氣工程學(xué)院2015年電力系統(tǒng)綜合課程設(shè)計(jì)報(bào)告班 級 160124姓 名 齊 濟(jì)學(xué) 號 16012429指導(dǎo)老師 陳 中小組成員 王武森 石應(yīng)隆目錄線路數(shù)據(jù)2一、網(wǎng)絡(luò)拓?fù)?二、發(fā)電機(jī)穩(wěn)態(tài)數(shù)據(jù)2三、變電站數(shù)據(jù)2四、線路數(shù)據(jù)3五、發(fā)電機(jī)暫態(tài)數(shù)據(jù)3典型方式電力系統(tǒng)潮流計(jì)算4一、發(fā)電廠變壓器選擇及參數(shù)計(jì)算4二、形成導(dǎo)納矩陣4三、求出系數(shù)矩陣B和B5四、設(shè)置各節(jié)點(diǎn)電壓初值5五、求功率不平衡量5六、解修正方程5七、求新的初值6八、判斷是否已收斂6九、計(jì)算結(jié)果6變電站變壓器選擇和主接線設(shè)計(jì)7一、變壓器選擇71.變電站主變壓器72.發(fā)電廠變壓器113.變壓器等效電路11二、主接線設(shè)計(jì)121.發(fā)電廠主接

2、線122.變電站主接線13三、總結(jié)與提升16線路斷路器選擇18仿真曲線計(jì)算20感想與建議23附錄一：PQ 解耦法潮流計(jì)算 Matlab 程序24附錄二：改進(jìn)歐拉法計(jì)算功角特性 Matlab 程序25線路數(shù)據(jù)一、 網(wǎng)絡(luò)拓?fù)鋱D 1二、 發(fā)電機(jī)穩(wěn)態(tài)數(shù)據(jù)名稱額定電壓機(jī)端電壓典型方式輸出有功G110KV1.05平衡機(jī)G210KV1.051.60G310KV1.050.90三、 變電站數(shù)據(jù)負(fù)荷額定電壓為10KV變電站編號負(fù)荷1負(fù)荷2負(fù)荷3最大負(fù)荷（100MVA）1.451.11.2典型方式負(fù)荷（100MVA）1.250.91.0最小負(fù)荷（100MVA）1.050.60.7Tmax （小時）50004000

3、4500 功率因數(shù) cos0.850.8 0.85低壓側(cè)電壓（kV）101010 對備用要求30%40%30% 對調(diào)壓要求逆調(diào)壓順調(diào)壓常調(diào)壓負(fù)荷重要性重要重要一般四、 線路數(shù)據(jù)名稱電壓電阻電抗對地導(dǎo)納（B）7-8110KV0.010.0850.0448-9110KV0.020.0920.0407-5110KV0.030.1020.0765-4110KV0.010.1620.0899-6110KV0.020.1750.0376-4110KV0.030.1830.042負(fù)序和正序相等、零序?yàn)檎虻?倍五、 發(fā)電機(jī)暫態(tài)數(shù)據(jù)名稱D軸電抗D軸暫態(tài)電抗阻尼系數(shù)Tj（慣性時間常數(shù)）G10.140.072.0

4、47.2G20.900.221.912.8G21.320.132.16.02典型方式電力系統(tǒng)潮流計(jì)算一、 發(fā)電廠變壓器選擇及參數(shù)計(jì)算三個發(fā)電廠的變壓器選擇如下：T1：單臺 SFP7-120000/110T2：雙臺 SFP7-90000/110T3：單臺 SFP7-90000/110經(jīng)計(jì)算得到三臺變壓器的阻抗、導(dǎo)納參數(shù)為：名稱電阻 R電抗 X對地電導(dǎo) G對地電納 BT10.00290.08750.001060.006T20.00210.05830.00170.0108T30.00420.1170.000850.0054二、 形成導(dǎo)納矩陣根據(jù)線路數(shù)據(jù)和發(fā)電機(jī)變壓器數(shù)據(jù)，可以算出電力系統(tǒng)網(wǎng)絡(luò)中各節(jié)點(diǎn)

5、的 自導(dǎo)納和各節(jié)點(diǎn)之間的互導(dǎo)納。y14=-8.54iy27=-11.423iy39=-17.13iy45=0.38-6.15iy46=0.872-5.32iy57=2.654-9.02iy69=0.645-5.64iy78=1.365-11.6iy89=2.256-10.38iy11=-8.5454iy22=-11.429iy33=-17.16iy44=1.56-19.945iy55=3.034-15.09iy66=1.52-10.92iy89=4.4-31.983iy88=3.621-21.938iy99=3.52-33.13i由此可得電力系統(tǒng)網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣為:-8.545i008.54

6、i00000-11.429i000011.423i000-17.16i000008.54i00-19.945i6.15i5.32i000006.15i-15.09i09.02i00005.32i0-10.92i00011.423i009.02i0-31.983i11.6i00000011.6i-21.983i0017.15i005.64i010.38i三、 求出系數(shù)矩陣B和B節(jié)點(diǎn) 2 作為平衡節(jié)點(diǎn)，節(jié)點(diǎn) 1、3 為 PV 節(jié)點(diǎn)，其余節(jié)點(diǎn)為 PQ 節(jié)點(diǎn)，由系 數(shù)矩陣的定義可得B =-8.54508.54000000-17.160000017.158.540-19.9456.155.3200000

7、6.15-15.0909.0200005.320-10.92005.640009.020-31.98311.600000011.6-21.98310.83017.15005.64010.38-33.13B =-19.9456.155.320006.15-15.0909.02005.320-10.92005.6409.020-31.98311.6000011.6-21.98310.83005.64010.38-33.13四、 設(shè)置各節(jié)點(diǎn)電壓初值Uii=1,2,3=1.05 ii=1,2,3=0Ujj=4,5,6,7,8,9=1 jj=4,5,6,7,8,9=0五、 求功率不平衡量由公式：Pi=P

8、is-Uij=1j=nUi(Gijcosij+Bijsinij)Qi=Qis-Uij=1j=nUi(Gijsinij-Bijcosij)可求得有功不平衡量Pi(0)和無功不平衡量Qi(0)，從而得到Pi(0)/Ui(0)和Qi0/Ui0(i=1,3,4,5,6,7,8,9)六、 解修正方程求解矩陣方程P/U=-BUQ/U=-BU得出各節(jié)點(diǎn)電壓相角的修正量i(0)和各節(jié)點(diǎn)電壓模值的修正量Ui(0)。七、 求新的初值求出修正后的值i(1)=i(0)+i(0)Ui(1)=Ui(0)+Ui(0)八、 判斷是否已收斂由收斂判斷依據(jù)i0 和 Ui0dd)|(abs(dtheta(9)dd)|(abs(dt

9、heta(3)dd)|( abs(dtheta(4)dd)|(abs(dtheta(5)dd)|(abs(dtheta(6)dd)|(abs(dtheta(7)dd)|(abs(dtheta(8)dd)|(abs(du(1)dd)|(abs(du(9)dd)|(abs(du(3)dd)|(abs(du(4)dd)|(abs(du(5)dd)|(abs(du(6)dd)|(abs (du(7)dd)|(abs(du(8)dd) dp1 = 0.9 - 1.05 * u(4) * 8.54 * sind(theta(1) - theta(4); dp3 = 1.6 - 1.05 * 17.15 *

10、 u(9) * sind(theta(3) - theta(9); dp4 = 0 - u(4) * (1.05 * 8.54 * sind(theta(4) - theta(1) + u(5) * (-0.38 * cosd(theta(4) - theta(5) + 6.15 * sind(theta(4) - theta(5) + u(6) * (-0.872 * cosd(theta(4) - theta(6) + 5.32 * sind(theta(4) - theta(6) + u(4) * 1.56); dp5 = -1.0625 - u(5) * (u(4) * (-0.38

11、* cosd(theta(5) - theta(4) + 6.15 * sind(theta(5) - theta(4) + u(7) * (-2.654 * cosd(theta(5) - theta(7) + 9.02 * sind(theta(5) - theta(7) + u(5) * 3.034); dp6 = -0.85 - u(6) * (u(4) * (-0.872 * cosd(theta(6) - theta(4) + 5.32 * sind(theta(6) - theta(4) + u(9) * (-0.645 * cosd(theta(6) - theta(9) +

12、5.64 * sind(theta(6) - theta(9) + u(6) * 1.52); dp7 = 0 - u(7) * (11.423 * 1.05 * sind(theta(7) - theta(2) + u(5) * (-2.654 * cosd(theta(7) - theta(5) + 9.02 * sind(theta(7) - theta(5) + u(8) * (-1.365 * cosd(theta(7) - theta(8) + 11.6 * sind(theta(7) - theta(8) + u(7) * 4.4); dp8 = -0.72 - u(8) * (

13、u(7) * (-1.365 * cosd(theta(8) - theta(7) + 11.6 * sind(theta(8) - theta(7) + u(9) * (-2.256 * cosd(theta(8) - theta(9) + 10.38 * sind(theta(8) - theta(9) + u(8) * 3.621); dp9 = 0 - u(9) * (17.15 * 1.05 * sind(theta(9) - theta(3) + u(6) * (-0.645 * cosd(theta(9) - theta(6) + 5.64 * sind(theta(9) - t

14、heta(6) + u(8) * (-2.256 * cosd(theta(9) - theta(8) + 10.38 * sind(theta(9) - theta(8) + u(9) * 3.52); dp_u = dp1/u(1) dp3/u(3) dp4/u(4) dp5/u(5) dp6/u(6) dp7/u(7) dp8/u(8) dp9/u(9); udtheta = -B1_ * dp_u; dtheta(1) = udtheta(1) / u(1); dtheta(3) = udtheta(2) / u(3); dtheta(4) = udtheta(3) / u(4); d

15、theta(5) = udtheta(4) / u(5); dtheta(6) = udtheta(5) / u(6); dtheta(7) = udtheta(6) / u(7); dtheta(8) = udtheta(7) / u(8); dtheta(9) = udtheta(8) / u(9); theta = theta + dtheta; dq4 = 0 - u(4) * (-1.05 * 8.54 * cosd(theta(4) - theta(1) + u(5) * (-0.38 * sind(theta(4) - theta(5) - 6.15 * cosd(theta(4

16、) - theta(5) + u(6) * (-0.872 * sind(theta(4) - theta(6) - 5.32 * cosd(theta(4) - theta(6) + 19.945 * u(4); dq5 = -0.1875 - u(5) * (u(4) * (-0.38 * sind(theta(5) - theta(4) - 6.15 * cosd(theta(5) - theta(4) + u(7) * (-2.654 * sind(theta(5) - theta(7) - 9.02 * cosd(theta(5) - theta(7) + 15.09 * u(5);

17、 dq6 = -0.15 - u(6) * (u(4) * (-0.872 * sind(theta(6) - theta(4) - 5.32 * cosd(theta(6) - theta(4) + u(9) * (-0.645 * sind(theta(6) - theta(9) - 5.64 * cosd(theta(6) - theta(9) + 10.92 * u(6); dq7 = 0 - u(7) * (-1.05 * 11.423 * cosd(theta(7) - theta(2) + u(5) * (-2.654 * sind(theta(7) - theta(5) - 9

18、.02 * cosd(theta(7) - theta(5) + u(8) * (-1.365 * sind(theta(7) - theta(8) - 11.6 * cosd(theta(7) - theta(8) + 31.983 * u(7); dq8 = -0.18 - u(8) * (u(7) * (-1.365 * sind(theta(8) - theta(7) - 11.6 * cosd(theta(8) - theta(7) + u(9) * (-2.256 * sind(theta(8) - theta(9) - 10.38 * cosd(theta(8) - theta(

19、9) + 21.938 * u(8); dq9 = 0 - u(9) * (-1.05 * 17.15 * cosd(theta(9) - theta(3) + u(6) * (-0.645 * sind(theta(9) - theta(6) - 5.64 * cosd(theta(9) - theta(6) + u(8) * (-2.256 * sind(theta(9) - theta(8) - 10.38 * cosd(theta(9) - theta(8) + 33.13 * u(9); dq_u = dq4/u(4) dq5/u(5) dq6/u(6) dq7/u(7) dq8/u

20、(8) dq9/u(9); dutemp = -B2_ * dq_u; du(4) = dutemp(1); du(5) = dutemp(2); du(6) = dutemp(3); du(7) = dutemp(4); du(8) = dutemp(5); du(9) = dutemp(6); u = u + du; i = i + 1; enduthetai附錄二：改進(jìn)歐拉法計(jì)算功角特性 Matlab 程序clearclcload(Yn.mat)G11 = real(Yn(1,1);G12 = real(Yn(1,2);G13 = real(Yn(1,3);G21 = real(Yn(2

21、,1);G22 = real(Yn(2,2);G23 = real(Yn(2,3);G31 = real(Yn(3,1);G32 = real(Yn(3,2);G33 = real(Yn(3,3);B12 = imag(Yn(1,2);B13 = imag(Yn(1,3);B21 = imag(Yn(2,1);B23 = imag(Yn(2,3);B31 = imag(Yn(3,1);B32 = imag(Yn(3,2);Tj1 = 6.02;Tj2 = 47.2;Tj3 = 12.8;E1 = abs(1.082 + 0.04i);E2 = abs(1.069 + 0.105i);E3 =

22、abs(1.098 + 0.298i);PT1 = 0.9;PT2 = 1.57;PT3 = 1.6;delta1 = 5.9; %1delta2 = 5.61; %2delta3 = 17.34; %3omega1 = 1; %1omega2 = 1; %2omega3 = 1; %3t = 0;h = 0.01;omg1 = omega1; %1omg2 = omega2; %2omg3 = omega3; %3dlt1 = delta1; %1dlt2 = delta2; %2dlt3 = delta3; %3PE1 = E1 * E1 * G11 + E1 * (E2 * (abs(Yn(1,2) * sind(delta1 - delta2 + atand(G12/B12) + E3 * (abs(Yn(1,3) * sind(delta1 - delta3 + atand(G13