Bài tập lớn giải tích mô phỏng hệ thống điện

TRƯỜNG ĐẠI HỌC QUY NHƠN KHOA KỸ THUẬT VÀ CÔNG NGHỆ BÀI TẬP LỚN MƠN HỌC : GIẢI TÍCH VÀ MÔ PHỎNG TRONG HỆ THỐNG ĐIỆN GVHD : TS TRẦN MẠNH HÀ Sinh viên thực : TRẦN NGỌC SƠN Lớp : Kỹ Thuật Điện- Điện tử ĐH54 Khoa : ĐIỆN Bình Định / Ngày .Tháng .Năm Đề bài: Cho sơ đồ hệ thống điện mẫu IEEE 30 nút hình vẽ sau Các liệu nút liệu nhánh hệ thống điện cho bảng tương ứng Các thông số cho hệ đơn vị tương đối ứng với lượng công suất 100MVA Yê u cầu sinh viên thực nội dung lớn sau đây: Sinh viên: Nguyễn Văn Lâm Lớp : Lớp KTĐ-ĐT K38A + Hệ số tải Kt (%) =108 + Mất đường dây: 10 - 20 Hình Sơ đồ mơt sợi hệ thống điện IEEE 30 nút I ỨNG DỤNG POWERWORLD: 1.Xây dựng sơ đồ mô 2.Xuất kết ma trận Ybus 3.Xuất kết điện áp nút 4.Xuất kết phân bổ công suất hệ thống điện II ỨNG DỤNG MATLAB TÍNH TỐN MA TRẬN Ybus 1.Viết hàm tính tốn 2.Xuất kết ma trậnYbus 3.So sánh với kết từ phần mềm Powerworld III ỨNG DỤNG MATLAB ( PHƯƠNG PHÁP GAUSS-SEIDEL) 1.Viết đoạn lệnh tính tốn chế độ xác lập 2.Xuất kết điện áp nút 3.Xuất kết phân bố công suất hệ thống điện IV ỨNG DỤNG MATLAB ( PHƯƠNG PHÁP NEWTON-RAPHSON) 1.Viết đoạn lệnh tính tốn chế độ xác lập 2.Xuất kết điện áp nút 3.Xuất kết phân bố công suất hệ thống điện Bảng Dữ liệu nút hệ thống điện IEEE 30 nút Nút Mã nút Điện áp (pu) Góc pha (độ) Pt Qt Pf Qf Qmin Qmax Qshunt (MW) (Mvar) (MW) (Mvar) (Mvar) (Mvar) (Mvar) 1 1.06 0 0 0 0 2 23,436 13,716 40 48.822 0 2,592 1,296 0 0 1.06 8,208 1,728 0 0 101,73 20,52 35.975 0 0 0 0 0 24,624 11,772 0 0 32,4 32,4 30.826 0 0 0 0 0 10 6,264 2,16 0 0 19 11 0 0 16.119 0 12 12,096 8,1 0 0 13 0 0 10.423 0 14 6,696 1,728 0 0 15 8,856 2,7 0 0 16 3,78 1,944 0 0 17 9,72 6,264 0 0 18 3,456 0,972 0 0 19 10,26 3,672 0 0 20 2,376 0,756 0 0 21 18,9 12,096 0 0 22 0 0 0 0 23 3,456 1,728 0 0 24 9,396 7,236 0 0 4.3 25 0 0 0 0 26 3,78 2,484 0 0 27 0 0 0 0 28 0 0 0 0 29 2,592 0,972 0 0 30 11,448 2,052 0 0 Bảng Dữ liệu nhánh hệ thống điện IEEE 30 nút Nút Nút đến R (pu) X (pu) B/2 (pu) k 0.0192 0.0575 0.0264 1 0.0452 0.1852 0.0204 0.057 0.1737 0.0184 0.0132 0.0379 0.0042 0.0472 0.1983 0.0209 0.0581 0.1763 0.0187 0.0119 0.0414 0.0045 0.046 0.116 0.0102 0.0267 0.082 0.0085 0.012 0.042 0.0045 0.208 0.978 10 0.556 0.969 11 0.208 10 0.11 12 0.256 0.932 12 13 0.14 12 14 0.1231 0.2559 12 15 0.0662 0.1304 12 16 0.0945 0.1987 14 15 0.221 0.1997 16 17 0.0824 0.1923 15 18 0.1073 0.2185 18 19 0.0639 0.1292 19 20 0.034 0.068 10 20 0.0936 0.209 10 17 0.0324 0.0845 10 21 0.0348 0.0749 10 22 0.0727 0.1499 21 22 0.0116 0.0236 15 23 0.1 0.202 22 24 0.115 0.179 23 24 0.132 0.27 24 25 0.1885 0.3292 25 26 0.2544 0.38 25 27 0.1093 0.2087 28 27 0.396 0.968 27 29 0.2198 0.4153 27 30 0.3202 0.6027 29 30 0.2399 0.4533 28 0.0636 0.2 0.0214 28 0.0169 0.0599 0.065 I ỨNG DỤNG POWERWORLD: Xây dựng sơ đồ mơ Hình 2: Xây dựng sơ đồ mô phần mềm Power world 20 Bảng 3: Dữ liệu nút sơ đồ xây dụng Case Information Bảng 4: Dữ liệu nhánh sơ đồ xây dụng Case Information Xuất kết ma trận Ybus j = sqrt(-1); Ybus = sparse(30); Ybus( 1, 1) = 6.4684+ j*( -20.6959); Ybus( 1, 2) = -5.2246+ j*( 15.6467); Ybus( 1, 3) = -1.2437+ j*( 5.0960); Ybus( 2, 1) = -5.2246+ j*( 15.6467); Ybus( 2, 2) = 9.7523+ j*( -30.6487); Ybus( 2, 4) = -1.7055+ j*( 5.1974); Ybus( 2, 5) = -1.1360+ j*( 4.7725); Ybus( 2, 6) = -1.6861+ j*( 5.1165); Ybus( 3, 1) = -1.2437+ j*( 5.0960); Ybus( 3, 3) = 9.4392+ j*( -28.6023); Ybus( 3, 4) = -8.1954+ j*( 23.5309); Ybus( 4, 2) = -1.7055+ j*( 5.1974); Ybus( 4, 3) = -8.1954+ j*( 23.5309); Ybus( 4, 4) = 16.3141+ j*( -54.9186); Ybus( 4, 6) = -6.4131+ j*( 22.3112); Ybus( 4, 12) = -0.0000+ j*( 3.9062); Ybus( 5, 2) = -1.1360+ j*( 4.7725); Ybus( 5, 5) = 4.0900+ j*( -12.1906); Ybus( 5, 7) = -2.9540+ j*( 7.4493); Ybus( 6, 2) = -1.6861+ j*( 5.1165); Ybus( 6, 4) = -6.4131+ j*( 22.3112); Ybus( 6, 6) = 22.3416+ j*( -82.4350); Ybus( 6, 7) = -3.5902+ j*( 11.0261); Ybus( 6, 8) = -6.2893+ j*( 22.0126); Ybus( 6, 9) = -0.0000+ j*( 4.8077); Ybus( 6, 10) = -0.0000+ j*( 1.7986); Ybus( 6, 28) = -4.3628+ j*( 15.4636); Ybus( 7, 5) = -2.9540+ j*( 7.4493); Ybus( 7, 6) = -3.5902+ j*( 11.0261); Ybus( 7, 7) = 6.5442+ j*( -18.4567); Ybus( 8, 6) = -6.2893+ j*( 22.0126); Ybus( 8, 8) = 7.7333+ j*( -26.5275); Ybus( 8, 28) = -1.4440+ j*( 4.5408); Ybus( 9, 6) = -0.0000+ j*( 4.8077); Ybus( 9, 9) = 0.0000+ j*( -18.7063); Ybus( 9, 10) = -0.0000+ j*( 9.0909); Ybus( 9, 11) = -0.0000+ j*( 4.8077); Ybus( 10, 6) = -0.0000+ j*( 1.7986); Ybus( 10, 9) = -0.0000+ j*( 9.0909); Ybus( 10, 10) = 16.7791+ j*( -48.3791); Ybus( 10, 17) = -3.9560+ j*( 10.3174); Ybus( 10, 21) = -10.2037+ j*( 21.9614); Ybus( 10, 22) = -2.6193+ j*( 5.4008); Ybus( 11, 9) = -0.0000+ j*( 4.8077); Ybus( 11, 11) = 0.0000+ j*( -4.8077); Ybus( 12, 4) = -0.0000+ j*( 3.9062); Ybus( 12, 12) = 6.5740+ j*( -24.4242); Ybus( 12, 13) = -0.0000+ j*( 7.1429); Ybus( 12, 14) = -1.5266+ j*( 3.1734); Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( Ybus( 12, 12, 13, 13, 14, 14, 14, 15, 15, 15, 15, 15, 16, 16, 16, 17, 17, 17, 18, 18, 18, 19, 19, 19, 20, 20, 21, 21, 21, 22, 22, 22, 22, 23, 23, 23, 24, 24, 24, 24, 25, 25, 25, 25, 26, 26, 27, 27, 27, 27, 27, 28, 28, 15) = 16) = 12) = 13) = 12) = 14) = 15) = 12) = 14) = 15) = 18) = 23) = 12) = 16) = 17) = 10) = 16) = 17) = 15) = 18) = 19) = 18) = 19) = 20) = 19) = 20) = 10) = 21) = 22) = 10) = 21) = 22) = 24) = 15) = 23) = 24) = 22) = 23) = 24) = 25) = 24) = 25) = 26) = 27) = 25) = 26) = 25) = 27) = 28) = 29) = 30) = 6) = 8) = -3.0954+ j*( 6.0973); -1.9520+ j*( 4.1044); -0.0000+ j*( 7.1429); 0.0000+ j*( -7.1429); -1.5266+ j*( 3.1734); 4.0175+ j*( -5.4243); -2.4910+ j*( 2.2509); -3.0954+ j*( 6.0973); -2.4910+ j*( 2.2509); 9.3655+ j*( -16.0116); -1.8108+ j*( 3.6874); -1.9683+ j*( 3.9761); -1.9520+ j*( 4.1044); 3.8346+ j*( -8.4979); -1.8826+ j*( 4.3935); -3.9560+ j*( 10.3174); -1.8826+ j*( 4.3935); 5.8386+ j*( -14.7110); -1.8108+ j*( 3.6874); 4.8865+ j*( -9.9062); -3.0757+ j*( 6.2188); -3.0757+ j*( 6.2188); 8.9580+ j*( -17.9835); -5.8824+ j*( 11.7647); -5.8824+ j*( 11.7647); 5.8824+ j*( -11.7647); -10.2037+ j*( 21.9614); 26.9783+ j*( -56.0891); -16.7746+ j*( 34.1277); -2.6193+ j*( 5.4008); -16.7746+ j*( 34.1277); 21.9345+ j*( -43.4829); -2.5405+ j*( 3.9544); -1.9683+ j*( 3.9761); 3.4298+ j*( -6.9653); -1.4614+ j*( 2.9892); -2.5405+ j*( 3.9544); -1.4614+ j*( 2.9892); 5.3118+ j*( -9.2313); -1.3099+ j*( 2.2876); -1.3099+ j*( 2.2876); 4.4957+ j*( -7.8650); -1.2165+ j*( 1.8171); -1.9693+ j*( 3.7602); -1.2165+ j*( 1.8171); 1.2165+ j*( -1.8171); -1.9693+ j*( 3.7602); 3.6523+ j*( -32.1877); -0.0001+ j*( 25.2525); -0.9955+ j*( 1.8810); -0.6875+ j*( 1.2940); -4.3628+ j*( 15.4636); -1.4440+ j*( 4.5408); Ybus( 28, 27) = -0.0001+ j*( 25.2525); Ybus( 28, 28) = 5.8069+ j*( -45.1705); Ybus( 29, 27) = -0.9955+ j*( 1.8810); Ybus( 29, 29) = 1.9076+ j*( -3.6044); Ybus( 29, 30) = -0.9121+ j*( 1.7234); Ybus( 30, 27) = -0.6875+ j*( 1.2940); Ybus( 30, 29) = -0.9121+ j*( 1.7234); Ybus( 30, 30) = 1.5995+ j*( -3.0173); V( 1) = 1.060000 + j*( 0.000000); V( 2) = 0.995182 + j*(-0.098041); V( 3) = 0.969361 + j*(-0.148009); V( 4) = 0.947499 + j*(-0.176061); V( 5) = 0.914801 + j*(-0.265331); V( 6) = 0.932580 + j*(-0.208049); V( 7) = 0.914207 + j*(-0.239650); V( 8) = 0.925352 + j*(-0.221626); V( 9) = 0.919149 + j*(-0.264609); V( 10) = 0.894842 + j*(-0.289569); V( 11) = 0.951679 + j*(-0.273974); V( 12) = 0.879909 + j*(-0.286524); V( 13) = 0.894654 + j*(-0.291325); V( 14) = 0.853239 + j*(-0.300140); V( 15) = 0.844490 + j*(-0.299521); V( 16) = 0.876134 + j*(-0.291349); V( 17) = 0.881698 + j*(-0.292585); V( 18) = 0.799505 + j*(-0.318362); V( 19) = 0.778231 + j*(-0.326775); V( 20) = 0.776154 + j*(-0.327650); V( 21) = 0.883733 + j*(-0.290959); V( 22) = 0.882822 + j*(-0.290266); V( 23) = 0.846208 + j*(-0.294144); V( 24) = 0.860788 + j*(-0.282906); V( 25) = 0.889692 + j*(-0.255032); V( 26) = 0.866980 + j*(-0.257574); V( 27) = 0.918473 + j*(-0.234502); V( 28) = 0.924169 + j*(-0.223749); V( 29) = 0.889217 + j*(-0.252942); V( 30) = 0.871108 + j*(-0.266416); Xuất kết điện áp nút Nom kV 132.00 132.00 132.00 132.00 132.00 132.00 132.00 132.00 1.00 Volt PU Volt (kV) 106.000 139.920 100.000 132.000 0.97301 128.437 0.95498 126.058 0.94490 124.727 0.94444 124.666 0.93490 123.406 0.93958 124.024 0.94488 0.945 dlieunut(:,3)=Vm'; dlieunut(:,4)=deltad'; dienap=dlieunut(:,3) gocpha=dlieunut(:,4) SLT = 0; fprintf('\n') fprintf(' Line Flow and Losses \n\n') fprintf(' Line Power at bus & line flow Line loss-Transformer\n') fprintf(' from to MW Mvar MVA MW Mvar tap\n') for n = 1:nbus busprt = 0; for L = 1:nbr; if busprt == fprintf(' \n'), fprintf('%6g', n), fprintf(' %9.3f', P(n)*Scb) fprintf('%9.3f', Q(n)*Scb), fprintf('%9.3f\n', abs(S(n)*Scb)) busprt = 1; else, end if nl(L)==n k = nr(L); In = (V(n) - a(L)*V(k))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(n); Ik = (V(k) - V(n)/a(L))*y(L) + Bc(L)*V(k); Snk = V(n)*conj(In)*Scb; Skn = V(k)*conj(Ik)*Scb; SL = Snk + Skn; SLT = SLT + SL; elseif nr(L)==n k = nl(L); In = (V(n) - V(k)/a(L))*y(L) + Bc(L)*V(n); Ik = (V(k) - a(L)*V(n))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(k); Snk = V(n)*conj(In)*Scb; Skn = V(k)*conj(Ik)*Scb; SL = Snk + Skn; SLT = SLT + SL; else, end if nl(L)==n | nr(L)==n fprintf('%12g', k), fprintf('%9.3f', real(Snk)), fprintf('%9.3f', imag(Snk)) fprintf('%9.3f', abs(Snk)), fprintf('%9.3f', real(SL)), if nl(L) ==n & a(L) ~= fprintf('%9.3f', imag(SL)), fprintf('%9.3f\n', a(L)) else, fprintf('%9.3f\n', imag(SL)) end else, end end end SLT = SLT/2; fprintf(' \n'), fprintf(' Total loss ') fprintf('%9.3f', real(SLT)), fprintf('%9.3f\n', imag(SLT)) 2.Xuất kết điện áp nút TTCDXACLAPPPG_S dienap = dienap = 1.0600 1.0000 0.9964 0.9829 1.0000 0.9874 0.9839 1.0000 1.0006 0.9975 1.0000 0.9956 1.0000 0.9757 0.9665 0.9884 0.9888 0.9341 0.9189 0.9175 0.9829 0.9833 0.9627 0.9672 0.9741 0.9541 0.9882 0.9872 0.9658 0.9529 gocpha = -5.5829 -8.8395 -10.7038 -16.3364 -12.7441 -14.8068 -13.8032 -16.0383 -17.8013 -16.0386 -17.7664 -17.7665 -19.0409 -19.2722 -18.1242 -18.1496 -21.1384 -22.0399 -22.1322 -18.4031 -18.4101 -19.4099 -19.1608 -18.7623 -19.2578 -18.2020 -13.5387 -19.6303 -20.6630 Xuất kết phân bố công suất hệ thống điện Line Flow and Losses Line Power at bus & line flow from to MW Mvar MVA Line loss Transformer MW Mvar tap 289.884 68.799 297.936 197.210 50.536 203.582 7.135 15.761 92.654 18.274 94.439 3.624 10.530 16.560 -67.783 69.777 -190.075 -34.774 193.230 7.135 15.761 49.187 -5.857 49.535 1.388 0.613 91.042 -14.904 92.254 3.990 12.582 66.400 -12.252 67.520 2.624 4.270 -2.590 -1.300 2.898 -89.031 -7.744 89.367 3.624 10.530 86.441 6.414 86.679 1.000 2.048 -8.210 -1.730 -47.799 6.470 -85.441 -4.367 74.609 -31.230 12 50.427 27.403 8.390 48.235 1.388 0.613 85.553 1.000 2.048 80.881 0.803 1.919 57.392 0.000 7.582 -101.740 46.519 111.871 -87.052 27.486 91.288 3.990 12.582 -14.691 19.024 24.036 0.284 -1.291 0.000 0.000 0.000 -63.775 16.521 65.881 -73.806 33.149 80.908 40.047 -8.801 41.002 32.472 -38.872 50.651 27.910 5.190 28.388 10 16.116 4.655 16.774 28 21.025 -11.809 24.115 2.624 0.803 0.457 0.312 0.000 0.000 0.082 0.932 4.270 1.919 -0.249 0.202 1.644 0.978 1.507 0.969 -12.383 -24.620 -11.770 27.289 14.976 -20.315 25.238 0.284 -1.291 -39.590 8.552 40.503 0.457 -0.249 -32.400 43.389 54.151 -32.161 39.074 50.607 0.312 0.202 28 -0.226 4.315 4.321 0.027 -4.142 0.000 0.000 -27.910 -3.545 11 0.002 0.289 10 27.915 3.249 10 -6.260 -16.116 -27.915 17 9.631 21 18.774 22 9.350 11 0.000 28.134 0.000 1.644 0.289 -0.000 0.000 28.104 0.000 0.868 16.840 17.966 -3.148 16.420 -2.381 28.017 6.615 11.684 10.818 21.668 4.945 10.578 0.000 0.000 0.044 0.164 0.082 0.000 -0.289 0.289 -0.002 -0.289 0.289 -0.000 1.507 0.868 0.116 0.353 0.169 0.000 12 -12.100 -8.100 14.561 -50.427 -19.821 54.183 0.000 7.582 13 0.002 -3.125 3.125 -0.000 0.014 14 9.922 3.075 10.388 0.134 0.279 15 24.487 10.036 26.464 0.468 0.921 16 3.920 1.739 4.289 0.018 0.037 13 0.000 3.139 3.139 12 -0.002 3.139 3.139 -0.000 0.014 14 -6.700 -1.730 6.920 12 -9.788 -2.796 10.180 0.134 0.279 15 3.086 1.065 3.265 0.025 0.022 15 -8.860 -2.700 9.262 12 -24.020 -9.114 25.691 0.468 0.921 14 -3.061 -1.043 3.234 0.025 0.022 18 16.605 6.415 17.802 0.364 0.741 23 1.619 1.036 1.922 0.004 0.008 16 -3.780 -1.940 4.249 12 -3.903 -1.703 4.258 0.018 0.037 17 0.124 -0.236 0.267 0.000 0.000 17 -9.720 -6.260 11.561 16 -0.124 0.236 0.267 0.000 0.000 10 -9.587 -6.499 11.582 0.044 0.116 18 -3.460 -0.970 3.593 15 -16.241 -5.674 17.204 0.364 0.741 19 12.782 4.710 13.622 0.136 0.275 19 -10.260 -3.670 10.897 18 -12.646 -4.436 13.401 0.136 0.275 20 2.380 0.764 2.500 0.003 0.005 20 -2.380 -0.760 19 -2.378 -0.759 2.498 2.496 0.003 0.005 21 -18.900 -12.100 22.441 10 -18.610 -10.465 21.350 0.164 0.353 22 -0.304 -1.634 1.662 0.000 0.001 22 0.000 0.000 0.000 10 -9.269 -4.777 10.427 0.082 0.169 21 0.304 1.635 1.663 0.000 0.001 24 8.978 3.140 9.512 0.108 0.167 23 -3.460 -1.730 3.868 15 -1.615 -1.028 1.915 0.004 0.008 24 -1.842 -0.701 1.971 0.006 0.011 24 -9.400 22 -8.871 23 1.848 25 -2.378 -2.940 9.849 -2.972 9.356 0.108 0.167 0.713 1.980 0.006 0.011 -0.677 2.473 0.012 0.022 25 0.000 24 2.390 26 3.835 27 -6.225 0.000 0.000 0.698 2.490 0.012 0.022 2.569 4.616 0.057 0.085 -3.273 7.033 0.057 0.109 26 -3.780 -2.480 25 -3.778 -2.483 27 25 28 29 30 0.000 6.282 -20.691 6.705 7.690 0.000 3.381 -7.067 1.846 1.847 4.521 4.521 0.057 0.085 0.000 7.134 0.057 0.109 21.864 0.000 1.939 6.955 0.109 0.206 7.909 0.205 0.386 28 0.000 0.000 27 20.691 9.005 0.252 -8.458 -20.944 -0.574 0.000 22.566 0.000 1.939 0.968 8.461 0.027 -4.142 20.952 0.082 -12.383 29 -2.590 -0.970 2.766 27 -6.597 -1.640 6.798 0.109 0.206 30 4.011 0.665 4.066 0.043 0.080 30 -11.450 -2.050 11.632 27 -7.485 -1.461 7.626 0.205 0.386 29 -3.969 -0.585 4.011 0.043 0.080 Total loss 23.761 47.405 IV ỨNG DỤNG MATLAB ( PHƯƠNG PHÁP NEWTON-RAPHSON) Viết đoạn lệnh tính tốn chế độ xác lập ns=0; ng=0; Vm=0; delta=0; yload=0; deltad=0; nbus = length(dlieunut(:,1)); for k=1:nbus n=dlieunut(k,1); kb(n)=dlieunut(k,2); Vm(n)=dlieunut(k,3); delta(n)=dlieunut(k, 4); Pd(n)=dlieunut(k,5); Qd(n)=dlieunut(k,6); Pg(n)=dlieunut(k,7); Qg(n) = dlieunut(k,8); Qmin(n)=dlieunut(k, 9); Qmax(n)=dlieunut(k, 10); Qsh(n)=dlieunut(k, 11); if Vm(n) = accuracy & iter Qmax(n), Vm(n) = Vm(n) - 0.01;end else, end else,end else,end end if kb(n) ~= A(nn,nn) = J11; DC(nn) = P(n)-Pk; end if kb(n) == A(nn,lm) = 2*Vm(n)*Ym(n,n)*cos(t(n,n))+J22; A(lm,nn)= J33; A(lm,lm) =-2*Vm(n)*Ym(n,n)*sin(t(n,n))-J44; DC(lm) = Q(n)-Qk; end end DX=A\DC';det(A) for n=1:nbus nn=n-nss(n); lm=nbus+n-ngs(n)-nss(n)-ns; if kb(n) ~= delta(n) = delta(n)+DX(nn); end if kb(n) == Vm(n)=Vm(n)+DX(lm); end end maxerror=max(abs(DC)); if iter == maxiter & maxerror > accuracy fprintf('\nWARNING: Iterative solution did not converged after ') fprintf('%g', iter), fprintf(' iterations.\n\n') fprintf('Press Enter to terminate the iterations and print the results \n') converge = 0; pause, else, end end if converge ~= tech= (' ITERATIVE SOLUTION DID NOT CONVERGE'); else, tech=(' Power Flow Solution by Newton-Raphson Method'); end V = Vm.*cos(delta)+j*Vm.*sin(delta); deltad=180/pi*delta; i=sqrt(-1); k=0; for n = 1:nbus if kb(n) == k=k+1; S(n)= P(n)+j*Q(n); Pg(n) = P(n)*Scb + Pd(n); Qg(n) = Q(n)*Scb + Qd(n) - Qsh(n); Pgg(k)=Pg(n); Qgg(k)=Qg(n); elseif kb(n) ==2 k=k+1; S(n)=P(n)+j*Q(n); Qg(n) = Q(n)*Scb + Qd(n) - Qsh(n); Pgg(k)=Pg(n); Qgg(k)=Qg(n); end yload(n) = (Pd(n)- j*Qd(n)+j*Qsh(n))/(Scb*Vm(n)^2); end dlieunut(:,3)=Vm'; dlieunut(:,4)=deltad'; dienap=dlieunut(:,3) gocpha=dlieunut(:,4) Pgt = sum(Pg); Qgt = sum(Qg); Pdt = sum(Pd); Qdt = sum(Qd); Qsht = sum(Qsh); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% dienap=dlieunut(:,3) gocpha=dlieunut(:,4) SLT = 0; fprintf('\n') fprintf(' Line Flow and Losses \n\n') fprintf(' Line Power at bus & line flow Line loss-Transformer\n') fprintf(' from to MW Mvar MVA MW Mvar tap\n') for n = 1:nbus busprt = 0; for L = 1:nbr; if busprt == fprintf(' \n'), fprintf('%6g', n), fprintf(' %9.3f', P(n)*Scb) fprintf('%9.3f', Q(n)*Scb), fprintf('%9.3f\n', abs(S(n)*Scb)) busprt = 1; else, end if nl(L)==n k = nr(L); In = (V(n) - a(L)*V(k))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(n); Ik = (V(k) - V(n)/a(L))*y(L) + Bc(L)*V(k); Snk = V(n)*conj(In)*Scb; Skn = V(k)*conj(Ik)*Scb; SL = Snk + Skn; SLT = SLT + SL; elseif nr(L)==n k = nl(L); In = (V(n) - V(k)/a(L))*y(L) + Bc(L)*V(n); Ik = (V(k) - a(L)*V(n))*y(L)/a(L)^2 + Bc(L)/a(L)^2*V(k); Snk = V(n)*conj(In)*Scb; Skn = V(k)*conj(Ik)*Scb; SL = Snk + Skn; SLT = SLT + SL; else, end if nl(L)==n | nr(L)==n fprintf('%12g', k), fprintf('%9.3f', real(Snk)), fprintf('%9.3f', imag(Snk)) fprintf('%9.3f', abs(Snk)), fprintf('%9.3f', real(SL)), if nl(L) ==n & a(L) ~= fprintf('%9.3f', imag(SL)), fprintf('%9.3f\n', a(L)) else, fprintf('%9.3f\n', imag(SL)) end else, end end end SLT = SLT/2; fprintf(' \n'), fprintf(' Total loss fprintf('%9.3f', real(SLT)), fprintf('%9.3f\n', imag(SLT)) TTCDXACLAPPPG_S dienap = 1.0600 1.0000 0.9964 0.9829 1.0000 0.9874 0.9839 1.0000 1.0006 0.9975 1.0000 0.9956 1.0000 0.9757 0.9665 0.9884 0.9888 0.9341 0.9189 0.9175 0.9829 0.9833 0.9627 0.9672 ') 0.9742 0.9542 0.9882 0.9873 0.9658 0.9529 gocpha = -5.5827 -8.8397 -10.7038 -16.3362 -12.7440 -14.8067 -13.8033 -16.0382 -17.8009 -16.0382 -17.7666 -17.7666 -19.0409 -19.2724 -18.1244 -18.1497 -21.1378 -22.0390 -22.1314 -18.4026 -18.4098 -19.4099 -19.1603 -18.7616 -19.2581 -18.2002 -13.5386 -19.6293 -20.6607 Xuất kết phân bổ cong suất hệ thống điện Line Flow and Losses Line Power at bus & line flow Line loss Transformer from to MW Mvar MVA MW Mvar tap 289.860 68.804 297.914 197.205 50.537 203.577 7.135 15.760 92.656 18.267 94.439 3.624 10.530 16.560 -67.791 69.785 -190.070 -34.777 193.225 7.135 15.760 49.189 -5.859 49.537 1.388 0.613 91.041 -14.904 92.253 3.990 12.582 66.400 -12.252 67.521 2.624 4.270 -2.590 -1.300 2.898 -89.032 -7.737 89.368 3.624 10.530 86.442 6.437 86.681 1.000 2.048 -8.210 -1.730 -47.800 6.472 -85.442 -4.389 74.605 -31.220 12 50.428 27.406 8.390 48.237 1.388 0.613 85.555 1.000 2.048 80.873 0.802 1.918 57.395 0.000 7.583 -101.740 46.512 111.868 -87.052 27.486 91.288 3.990 12.582 -14.688 19.027 24.037 0.284 -1.291 0.000 0.000 0.000 -63.776 16.521 65.881 -73.802 33.138 80.900 40.049 -8.797 41.004 32.483 -38.876 50.660 27.910 5.188 28.388 10 16.115 4.655 16.774 28 21.021 -11.830 24.122 2.624 4.270 0.802 1.918 0.457 -0.249 0.312 0.202 0.000 1.644 0.978 0.000 1.506 0.969 0.082 -12.383 -24.620 -11.770 27.289 14.972 -20.317 25.238 0.284 -1.291 -39.592 8.547 40.505 0.457 -0.249 -32.400 43.388 54.151 -32.171 39.078 50.617 0.312 0.202 28 -0.229 4.310 4.316 0.026 -4.142 0.000 0.000 -27.910 -3.544 11 -0.000 0.290 10 27.910 3.253 10 -6.260 -16.115 -27.910 17 9.641 21 18.772 22 9.351 11 0.000 28.134 0.000 1.644 0.290 0.000 0.000 28.099 0.000 0.867 16.840 17.966 -3.149 16.419 -2.386 28.012 6.613 11.691 10.818 21.666 4.945 10.578 0.932 0.000 0.000 0.045 0.164 0.082 0.000 -0.290 0.290 0.000 -0.290 0.290 0.000 1.506 0.867 0.116 0.353 0.169 0.000 12 -12.100 -8.100 14.561 -50.428 -19.824 54.185 0.000 7.583 13 -0.000 -3.127 3.127 0.000 0.014 14 9.921 3.074 10.386 0.134 0.278 15 24.487 10.036 26.463 0.468 0.921 16 3.921 1.741 4.290 0.018 0.037 13 0.000 12 0.000 14 3.140 3.140 3.140 3.140 0.000 0.014 -6.700 -1.730 6.920 12 -9.787 -2.796 10.178 0.134 0.278 15 3.087 1.066 3.266 0.025 0.022 15 12 14 18 23 -8.860 -24.019 -3.062 16.602 1.619 -2.700 -9.114 -1.043 6.421 1.037 9.262 25.690 0.468 0.921 3.235 0.025 0.022 17.801 0.364 0.741 1.922 0.004 0.008 16 -3.780 -1.940 4.249 12 -3.903 -1.704 4.259 0.018 0.037 17 0.123 -0.236 0.267 0.000 0.000 17 -9.720 -6.260 11.561 16 -0.123 0.236 0.267 0.000 0.000 10 -9.597 -6.496 11.589 0.045 0.116 18 -3.460 -0.970 3.593 15 -16.238 -5.680 17.203 0.364 0.741 19 12.778 4.710 13.619 0.136 0.275 19 -10.260 -3.670 10.897 18 -12.643 -4.435 13.398 0.136 0.275 20 2.383 0.765 2.502 0.003 0.005 20 -2.380 -0.760 19 -2.380 -0.760 2.498 2.498 0.003 0.005 21 -18.900 -12.100 22.441 10 -18.608 -10.464 21.348 0.164 0.353 22 -0.292 -1.636 1.662 0.000 0.001 22 0.000 10 -9.269 21 0.292 24 8.977 23 -3.460 -1.730 3.868 15 -1.615 -1.029 1.915 0.004 0.008 24 -1.845 -0.701 1.974 0.006 0.011 0.000 -4.776 1.636 3.140 0.000 10.428 0.082 0.169 1.662 0.000 0.001 9.510 0.108 0.167 24 -9.400 -2.940 9.849 22 -8.869 -2.972 9.354 0.108 0.167 23 1.851 0.713 1.983 0.006 0.011 25 -2.381 -0.680 2.476 0.012 0.022 25 0.000 24 2.393 26 3.837 27 -6.231 26 -3.780 -2.480 25 -3.780 -2.480 27 25 28 29 30 0.000 6.288 -20.684 6.707 7.689 0.000 0.000 0.702 2.494 0.012 0.022 2.565 4.616 0.057 0.085 -3.267 7.035 0.057 0.109 0.000 3.376 -7.068 1.844 1.848 4.521 4.521 0.057 0.085 0.000 7.137 0.057 0.109 21.858 0.000 1.938 6.956 0.109 0.206 7.908 0.205 0.386 28 0.000 0.000 27 20.684 9.006 0.256 -8.452 -20.940 -0.553 29 -2.590 -0.970 2.766 27 -6.598 -1.638 6.799 0.109 0.206 30 4.008 0.668 4.064 0.042 0.080 30 -11.450 -2.050 11.632 27 -7.484 -1.462 7.626 0.205 0.386 29 -3.966 -0.588 4.009 0.042 0.080 Total loss 0.000 22.559 0.000 1.938 0.968 8.456 0.026 -4.142 20.947 0.082 -12.383 23.760 47.403 THE END ... suất hệ thống điện Bảng 5: Kết phân bổ công suất lưới điện II ỨNG DỤNG MATLAB TÍNH TỐN MA TRẬN clc; clear all; disp(' _ BÀI TẬP LỚNN _ '); disp(' | GIẢI TÍCH VÀ MƠ PHỎNG HỆ THỐNG ĐIỆN...Đề bài: Cho sơ đồ hệ thống điện mẫu IEEE 30 nút hình vẽ sau Các liệu nút liệu nhánh hệ thống điện cho bảng tương ứng Các thông số cho hệ đơn vị tương đối ứng với lượng... độ xác lập 2.Xuất kết điện áp nút 3.Xuất kết phân bố công suất hệ thống điện Bảng Dữ liệu nút hệ thống điện IEEE 30 nút Nút Mã nút Điện áp (pu) Góc pha (độ) Pt Qt Pf Qf Qmin Qmax Qshunt (MW) (Mvar)