BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ HUỲNH HỒNG THANH PHONG XÁC ĐỊNH VỊ TRÍ VÀ DUNG LƯỢNG TCSC DỰA TRÊN MẶT CẮT TỐI THIỂU CẢI TIẾN NHẰM NÂNG CAO CÔNG SUẤT LƯỚI ĐIỆN TRUYỀN TẢI NGÀNH: KỸ THUẬT ĐIỆN - 8520201 SKC008018 Tp Hồ Chí Minh, tháng 2/2023 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ HUỲNH HỒNG THANH PHONG XÁC ĐỊNH VỊ TRÍ VÀ DUNG LƯỢNG TCSC DỰA TRÊN MẶT CẮT TỐI THIỂU CẢI TIẾN NHẰM NÂNG CAO CÔNG SUẤT LƯỚI ĐIỆN TRUYỀN TẢI NGÀNH: KỸ THUẬT ĐIỆN - 8520201 Hướng dẫn khoa học: PGS.TS TRƯƠNG VIỆT ANH Tp Hồ Chí Minh, tháng …/2023 17 17 18 18 19 19 20 20 21 18 22 21 21 20 20 23 23 22 0.0018 0.0135 0.0033 0.0033 0.0051 0.0051 0.0028 0.0028 0.0087 0.0144 0.1053 0.0259 0.0259 0.0396 0.0396 0.0216 0.0216 0.0678 0.0303 0.2212 0.0545 0.0545 0.0833 0.0833 0.0455 0.0455 0.1424 500 600 625 500 600 625 500 600 625 500 600 625 500 600 625 500 600 625 500 600 625 500 600 625 500 600 625 0 0 0 0 0 1 1 1 1 -360 -360 -360 -360 -360 -360 -360 -360 -360 360; 360; 360; 360; 360; 360; 360; 360; 360; ]; %% - OPF Data -%% %% generator cost data % startup shutdown n x1 y1 xn yn % startup shutdown n c(n-1) c0 mpc.gencost = [ % bus Pmin Pmax Qmin Qmax Unit Code 1500 130 400.6849; % 16 20 10 U20 1500 130 400.6849; % 16 20 10 U20 1500 0.014142 16.0811 212.3076; % 15.2 76 -25 30 U76 1500 0.014142 16.0811 212.3076; % 15.2 76 -25 30 U76 1500 130 400.6849; % 16 20 10 U20 1500 130 400.6849; % 16 20 10 U20 1500 0.014142 16.0811 212.3076; % 15.2 76 -25 30 U76 1500 0.014142 16.0811 212.3076; % 15.2 76 -25 30 U76 1500 0.052672 43.6615 781.521; % 25 100 60 U100 1500 0.052672 43.6615 781.521; % 25 100 60 U100 1500 0.052672 43.6615 781.521; % 25 100 60 U100 1500 0.00717 48.5804 832.7575; % 13 69 197 80 U197 1500 0.00717 48.5804 832.7575; % 13 69 197 80 U197 1500 0.00717 48.5804 832.7575; % 13 69 197 80 U197 1500 0 0; % 14 SynCond 1500 0.328412 56.564 86.3852; % 15 2.4 12 U12 1500 0.328412 56.564 86.3852; % 15 2.4 12 U12 1500 0.328412 56.564 86.3852; % 15 2.4 12 U12 1500 0.328412 56.564 86.3852; % 15 2.4 12 U12 1500 0.328412 56.564 86.3852; % 15 2.4 12 U12 1500 0.008342 12.3883 382.2391; % 15 54.3 155 -50 80 U155 1500 0.008342 12.3883 382.2391; % 16 54.3 155 -50 80 U155 69 1500 1500 1500 1500 1500 1500 1500 1500 1500 U155 1500 U155 1500 ]; 0 0 0 0 3 3 3 3 0.000213 0.000213 0.001 0.001 0.001 0.001 0.001 0.001 0.008342 4.4231 395.3749; % 18 100 400 -50 200 U400 4.4231 395.3749; % 21 100 400 -50 200 U400 0.001; % 22 10 50 -10 16 U50 0.001; % 22 10 50 -10 16 U50 0.001; % 22 10 50 -10 16 U50 0.001; % 22 10 50 -10 16 U50 0.001; % 22 10 50 -10 16 U50 0.001; % 22 10 50 -10 16 U50 12.3883 382.2391; % 23 54.3 155 -50 80 0.008342 12.3883 382.2391; % 23 54.3 155 -50 80 0.004895 11.8495 665.1094; % 23 140 350 -25 150 U350 Chương trình chạy phân bố công suất cho lưới điện với 140% công suất clear clc close all mpc = loadcase('case24_ieee_rts'); HeSoQuaTai=1.4; mpc.bus(:,3)=mpc.bus(:,3)*HeSoQuaTai; mpc.bus(:,4)=mpc.bus(:,4)*HeSoQuaTai; result=runpf(mpc); figure plot(result.bus(:,1),result.bus(:,8),'-*k') grid on grid minor % legend('Vol') title('Cong suat tai la 140%') xlabel('nut') ylabel('Dien Ap (p.u.)') figure nBranch=size(result.branch,1); Temp=sqrt(result.branch(:,14).^2+result.branch(:,15).^2); plot(1:nBranch,100*abs(Temp./result.branch(:,6)),'-*k') 70 % plot(1:nBranch,100*abs(result.branch(:,14)./result.branch(:,6)),'-*k') grid on grid minor % legend('Vol') title('Cong suat tai la 140%') xlabel('nhanh') ylabel('Cong suat (%)') chương trình xác định trạng thái tải hay thấp áp lưới điện 24 nút clear clc close all mpc = loadcase('case24_ieee_rts'); HeSoQuaTai=1.35; mpc.bus(:,3)=mpc.bus(:,3)*HeSoQuaTai; mpc.bus(:,4)=mpc.bus(:,4)*HeSoQuaTai; i=1; for step1=1:size(mpc.branch,1) j=1; for step2=-0.8:0.1:0.2 data(i,1)=step1; % nhanh data(i,2)=step2; % Xtcsc data(i,9)=j; i=i+1; j=j+1; end end for step=1:size(data,1) % step=1; mpc1=mpc; mpc1.branch( data(step,1),4)=mpc1.branch( data(step,1),4)*(1+data(step,2)); result=runpf(mpc1); Temp=sqrt(result.branch(:,14).^2+result.branch(:,15).^2); Temp= 100*abs(Temp./result.branch(:,6)); [a, b]=max(Temp); data(step,3)=b; 71 data(step,4)=a; if a>100 data(step,7)=1; end [c, d]=min(result.bus(:,8)); data(step,5)=d; data(step,6)=c; if c