Đây là luận văn do tôi thực hiện. Các số liệu, kết quả nêu trong luận văn là trung thực và chưa từng được ai công bố trong bất kì công trình nào khác. Dữ liệu tham khảo trong luận văn này đều được dẫn nguồn và có độ chính xác cao nhất trong phạm vi hiểu biết của tôi.Đây là luận văn do tôi thực hiện. Các số liệu, kết quả nêu trong luận văn là trung thực và chưa từng được ai công bố trong bất kì công trình nào khác. Dữ liệu tham khảo trong luận văn này đều được dẫn nguồn và có độ chính xác cao nhất trong phạm vi hiểu biết của tôi.
ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA ĐIỆN-ĐIỆN TỬ BỘ MÔN HỆ THỐNG ĐIỆN - - LUẬN VĂN TỐT NGHIỆP ĐẠI HỌC ỨNG DỤNG SVC ĐỂ ĐIỀU KHIỂN ĐIỆN ÁP VÀ DỊNG CƠNG SUẤT PHẢN KHÁNG TRONG HỆ THỐNG ĐIỆN GVHD: TS.NGUYỄN VĂN LIÊM SVTH: HOÀNG MẠNH CƯỜNG MSSV: 41100450 TP.HỒ CHÍ MINH, THÁNG 12 NĂM 2015 LỜI CAM KẾT Tơi cam kết: • Đây luận văn tơi thực • Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác • Dữ liệu tham khảo luận văn dẫn nguồn có độ xác cao phạm vi hiểu biết LỜI CÁM ƠN Lời đầu tiên, tơi xin tỏ lòng biết ơn sâu sắt tới TS Nguyễn Văn Liêm, người tận tình giúp đỡ, bảo tạo điều kiện tốt để tơi hồn thành luận văn Tôi xin chân thành cảm ơn quý thầy, cô Bộ môn Hệ thống điện, Trường Đại học Bách Khoa Tp Hồ Chí Minh tận tình truyền đạt kiến thức năm học tập Vốn kiến thức tiếp thu trình học tản cho trình nghiên cứu làm đề tài luận văn tơi, chúng hành trang q báu để bước vào đời cách vững tự tin Tơi xin cảm ơn gia đình, anh, em bạn bè bên cạnh giúp đỡ suốt năm học, đặc biệt trình thực luận văn Với vốn kiến thức, kinh nghiệm thời gian hạn chế nên khơng tránh khỏi sai sót q trình thực đề tài luận văn Tơi kính mong nhận ý kiến phê bình, đóng góp q thầy, để luận văn tơi hồn thiện Tp Hồ Chí Minh, ngày….tháng….năm 2015 Sinh viên Hồng Mạnh Cường MỤC LỤC TÓM TẮT LUẬN VĂN Ngày nay, với phát triển nhanh chóng kinh tế, hệ thống điện ngày phát triển với quy mô độ phức tạp ngày lớn Do vấn đề đặt phải đảm bảo cho hệ thống vận hành an toàn hiệu Với đường dây truyền tải dài, việc đảm bảo điện áp giữ ổn định ranh giới cho phép yêu cầu bắt buộc để giữ ổn định hệ thống Nhằm đáp ứng yêu cầu này, thiết bị điều khiển mới, hoạt động dựa tản thiết bị điện tử công suất áp dụng ngày rộng rãi vào việc điều khiển hệ thống điện Các thiết bị gọi chung thiết bị điều khiển xoay chiều linh hoạt FACTS Có nhiều loại thiết bị FACTS thường dùng hệ thống, điển SVC, STACOM, UPFC, TCSC…Trong đó, thiết bị SVC ( Static VAR Compensor) sử dụng công cụ giúp cân công suất nút, qua đảm bảo điện áp giữ ổn định giới hạn cho phép Nhiệm vụ luận văn tìm hiểu cấu tạo, nguyên lý hoạt động SVC từ xây dựng mơ hình SVC chế độ xác lập Thực mô mạng điện thực tế thấy lợi ích vượt trội sử dụng SVC hệ thống điện Những kết thu từ luận văn gồm: - Xây dựng mơ hình SVC chế độ xác lập - Xác định nguyên lý điều khiển điện áp bù công suất phản kháng SVC - Kiểm chứng khả điều khiển điện áp SVC tất nút tải hệ thống điện đơn giản nút giản đồ điện áp - Kiểm chứng khả bù công suất nút mạng điện đơn giản nút giản đồ điện áp - Đã khảo sát mức độ tác động yếu tố độ dốc, giới hạn điện nạp, giá trị đặt, lên q trình điều khiển điện áp bù cơng suất SVC Từ đó, đưa nhận xét để điều chỉnh yếu tố cách phù hợp, để tối ưu hiệu hoạt động SVC DANH MỤC CÁC HÌNH ẢNH DANH SÁCH CÁC BẢNG BIỂU CHÚ GIẢI CÁC TỪ VIẾT TẮT Ký hiệu Định nghĩa Đơn vị FACTS Hệ thống truyền tải điện xoay chiều linh hoạt - SVC Thiết bị bù tĩnh điều khiển Thyristor - STACOM Thiết bị bù ngang điều khiển Thyristor - TCSC Thiết bị bù dọc điều khiển Thyristor - UPFC Thiết bị điều khiển dòng cơng suất hợp - TCR Cuộn kháng điều khiển Thyristor - TSR Tụ điện đóng mở Thyristor - nSVC Số SVC lắp hệ thống - Scal ( k ) Công suất nút thứ k absV ( k ) Biên độ điện áp nút thứ k kV δ (k ) Góc pha điện áp nút thứ k radian kload (k ) Mảng chứa thống số nút tải thứ k - kgen(k ) Mảng chứa thông số nút máy phát thứ k - kSVC (k ) Mảng chứa thông số nút SVC - hSVC (k ) Nút điện áp cao SVC thứ k - lSVC (k ) Nút điện áp thấp SVC thứ k - apha(k ) Độ dốc SVC thứ k % B Điện nạp SVC - Yij Biên độ tổng dẫn nhánh - θ ij Góc pha tổng dẫn nhánh radian Qinj (k ) Công suất phản kháng bù nút thứ k MVAr Qgen ( k ) Công suất phản kháng máy phát thứ k MVAr Qgen ( k ) Công suất phản kháng cực tiểu máy phát thứ k MVAr Qgen max ( k ) Công suất phản kháng cực đại máy phát thứ k MVAr XB(k ) 10 Điện kháng máy biến áp SVC thứ k MVA - PHỤ LỤC B Jagen22(k,n) = -imag(Scal(n))-absV(n)*absV(n)*imag(ybus(n,n)); else Jagen11(k,n) = 0; Jagen21(k,n) = abs((absV(kgen(k)) +j*theta(kgen(k)))*ybus(kgen(k),n))*cos(angle(ybus(kgen(k),n))+ theta(n)-theta(kgen(k))); Jagen22(k,n)=-abs((absV(n)+j*theta(n))*(absV(kgen(k)) +j*theta(kgen(k)))*ybus(kgen(k),n))*sin(angle(ybus(kgen(k),n))+ theta(n)-theta(kgen(k))); end end end else for k=1:G Qgen(k) = imag(Scal(kgen(k))-Sg_load(k))-Qinj(kgen(k)); if round(Qgen(k)*1e3)/1e3 =Qgenmax(k) fgen1(k)=imag(Scal(kgen(k)))-imag(Sg_load(k))-Qgenmax(k); fgen2(k) = real(Scal(kgen(k))) -Pgen(k)-real(Sg_load(k)); kt=1; else fgen1(k) = absV(kgen(k))-Vgen(k); fgen2(k) = real(Scal(kgen(k))) -Pgen(k)-real(Sg_load(k)); kt=0; end for n = 1:N if kt==1 if n == kgen(k) Jagen11(k,n)=imag(Scal(n))/absV(n)-(absV(n))*imag(ybus(n,n)); Jagen12(k,n)=real(Scal(n))-absV(n)*absV(n)*real(ybus(n,n)); Jagen21(k,n) = real(Scal(n))/absV(n)+absV(n)*real(ybus(n,n)); Jagen22(k,n) = -imag(Scal(n))-(absV(n))^2*imag(ybus(n,n)); else Jagen11(k,n)=-absV(kgen(k))*abs(ybus(kgen(k),n))*sin(angle(ybus(kgen(k),n))+ theta(n)theta(kgen(k))); Jagen12(k,n)=-absV(n)*absV(kgen(k))*abs(ybus(kgen(k),n))*cos(angle(ybus(kgen(k),n))+ theta(n)-theta(kgen(k))); Hoàng Mạnh Cường 122 PHỤ LỤC B Jagen21(k,n) = absV(kgen(k))*abs(ybus(kgen(k),n))*cos(angle(ybus(kgen(k),n))+ theta(n)-theta(kgen(k))); Jagen22(k,n)=-absV(n)*absV(kgen(k))*abs(ybus(kgen(k),n))*sin(angle(ybus(kgen(k),n))+ theta(n)-theta(kgen(k))); end elseif kt==0 Jagen12(k,n) = 0; if n == kgen(k) Jagen11(k,n) = Jagen21(k,n) = Jagen22(k,n) = else Jagen11(k,n) = Jagen21(k,n) = theta(kgen(k))); 1; real(Scal(n))/absV(n)+absV(n)*real(ybus(n,n)); -imag(Scal(n))-(absV(n))^2*imag(ybus(n,n)); 0; abs(V(kgen(k))*ybus(kgen(k),n))*cos(angle(ybus(kgen(k),n))+ theta(n)- Jagen22(k,n)=-abs(V(kgen(k))*V(n)*ybus(kgen(k),n))*sin(angle(ybus(kgen(k),n))+ theta(n)-theta(kgen(k))); end end end end end %% process_SVC.m %%su ly nut SVC %% if A(1,3) > for k=1:nSVC B(k)=-(absV(hSVC(k))-absV(lSVC(k))*cos(theta(lSVC(k))-theta(hSVC(k))))/(XB(k)*absV(hSVC(k))); if (round(B(k)*10^3)/10^3 = Bmax(k) Bcal(k)=Bmax(k); fSVC1(k)= (absV(hSVC(k))-absV(lSVC(k))*cos(theta(lSVC(k))-theta(hSVC(k))))/XB(k) +Bmax(k)*absV(hSVC(k)); fSVC2(k)=(absV(hSVC(k))*absV(lSVC(k))*sin(theta(lSVC(k))-theta(hSVC(k))))/XB(k); kt=1; else Bcal(k)=B(k); if typeSVC(k)==0 fSVC1(k)=absV(hSVC(k))-Vref(k)-apha(k)*(absV(hSVC(k))-absV(lSVC(k))*cos(theta(lSVC(k))theta(hSVC(k))))/XB(k); else % fSVC1(k)=imag(V(hSVC(k))*conj((V(lSVC(k))-V(hSVC(k)))/(XB(k)*j)))-Qref(k); % end fSVC2(k)=(absV(hSVC(k))*absV(lSVC(k))*sin(theta(lSVC(k))-theta(hSVC(k))))/XB(k); kt=0; end for n=1:N if n==lSVC(k) if kt==0 if typeSVC(k)==0 JaSVC11(k,n)=apha(k)*cos(theta(lSVC(k))-theta(hSVC(k)))/XB(k); JaSVC12(k,n)=-absV(lSVC(k))*apha(k)*sin(theta(lSVC(k))-theta(hSVC(k)))/XB(k); else % -JaSVC11(k,n)= (absV(hSVC(k))*cos(theta(hSVC(k))-theta(lSVC(k))))/XB(k); JaSVC12(k,n)=absV(hSVC(k))*absV(lSVC(k))*sin(theta(hSVC(k))-theta(lSVC(k)))/XB(k); % end elseif kt==1 JaSVC11(k,n)=-cos(theta(lSVC(k))-theta(hSVC(k)))/XB(k); JaSVC12(k,n)=absV(lSVC(k))*sin(theta(lSVC(k))-theta(hSVC(k)))/XB(k); Hoàng Mạnh Cường 124 PHỤ LỤC B end JaSVC21(k,n)=absV(hSVC(k))/XB(k)*sin(theta(lSVC(k))-theta(hSVC(k))); JaSVC22(k,n)=absV(hSVC(k))*absV(lSVC(k))/XB(k)*cos(theta(lSVC(k))-theta(hSVC(k))); elseif n==hSVC(k) if kt==0; if typeSVC(k)==0 JaSVC11(k,n)=1-apha(k)*absV(lSVC(k))/XB(k); JaSVC12(k,n)=absV(lSVC(k))*apha(k)*sin(theta(lSVC(k))-theta(hSVC(k)))/XB(k); else % -JaSVC11(k,n)= (absV(lSVC(k))*cos(theta(hSVC(k))-theta(lSVC(k)))2*absV(hSVC(k)))/XB(k); JaSVC12(k,n)= -absV(lSVC(k))*absV(hSVC(k))*sin(theta(hSVC(k))theta(lSVC(k)))/XB(k); % end elseif kt==1; JaSVC11(k,n)=1/XB(k)+Bcal(k); JaSVC12(k,n)=-absV(lSVC(k))*sin(theta(lSVC(k))-theta(hSVC(k)))/XB(k); end JaSVC21(k,n)=absV(lSVC(k))/XB(k)*sin(theta(lSVC(k))-theta(hSVC(k))); JaSVC22(k,n)=-absV(hSVC(k))*absV(lSVC(k))/XB(k)*cos(theta(lSVC(k))-theta(hSVC(k))); else if typeSVC(k)==0 JaSVC11(k,n)=0; JaSVC12(k,n)=0; else % JaSVC11(k,n)=0; JaSVC12(k,n)=0; % end Hoàng Mạnh Cường 125 PHỤ LỤC B JaSVC21(k,n)=0; JaSVC22(k,n)=0; end end %% end end caculate_current_and_power_loss.m %% tinh dong dien va dong cong suat tren cac nhanh nhanh=[nl nr]; for k=1:length(nhanh) pp=nhanh(k,1); qq=nhanh(k,2); In(k)=(V(pp)-V(qq))*(y(k)/a(k))+V(pp)*((1/a(k))*(1/a(k)-1))*y(k)+V(pp)*Bc(k); Ik(k)=(V(qq)-V(pp))*(y(k)/a(k))+V(qq)*(1-1/a(k))*y(k)+V(qq)*Bc(k); Iloss(k)=In(k)+Ik(k); Snk(k)=V(pp)*conj(In(k))*baseVA; Skn(k)=V(qq)*conj(Ik(k))*baseVA; SL(k)=Snk(k)+Skn(k); end store_data.m if A(1,3) > BB=[BB Bcal]; end VV=[VV absV]; ttheta=[ttheta theta]; SScal=[SScal Scal]; PLL=[PLL sum(real(SL))]; QLL=[QLL sum(imag(SL))]; plot_V_bus.m x=[1:1:5]; load('PV0_data.mat'); hold on; grid on; Hoàng Mạnh Cường 126 PHỤ LỤC B xlabel('bus'); ylabel('Voltage'); if module==1 if A(1,3)==0 plot(x,VV0','c','LineWidth',2); legend('dien ap he thong khong co SVC'); else absVV=absV; absVV(6,:)=[]; plot(x,absVV','c','LineWidth',2); legend('dien ap he thong co SVC'); end elseif module==2 if A(1,3)>0 VV(6,:)=[]; title('khao sat anh huong cua SVC len dien ap'); plot(x,VV(:,1)','r','LineWidth',1.5); plot(x,VV(:,2)','b','LineWidth',1.5); plot(x,VV(:,3)','g','LineWidth',1.5); plot(x,VV0','y','LineWidth',1.5); legend('SVC dat tai nut so 3','SVC dat tai nut so 4','SVC dat tai nut so 5','khong co SVC'); else plot(x,VV0','c','LineWidth',2); legend('khong co SVC'); end elseif module==3 if A(1,3)>0 VV(6,:)=[]; title('khao sat anh huong cua gioi han len dien ap'); plot(x,VV(:,1)','r','LineWidth',1.5); plot(x,VV(:,2)','b','LineWidth',1.5); plot(x,VV(:,3)','g','LineWidth',1.5); plot(x,VV(:,4)','y','LineWidth',1.5); plot(x,VV(:,5)','k','LineWidth',1.5); legend('Bmax=0.05','Bmax=0.1','Bmax=0.15','Bmax=0.2','Bmax=0.25'); else plot(x,VV0','c','LineWidth',2); legend('khong co SVC'); Hoàng Mạnh Cường 127 PHỤ LỤC B end elseif module==4 if A(1,3)>0 VV(6,:)=[]; title('khao sat anh huong cua doc apha len dien ap'); plot(x,VV(:,1)','r','LineWidth',1.5); plot(x,VV(:,2)','b','LineWidth',1.5); plot(x,VV(:,3)','g','LineWidth',1.5); plot(x,VV(:,4)','y','LineWidth',1.5); plot(x,VV(:,5)','k','LineWidth',1.5); legend('apha=0','pha=2','apha=4','apha=6','apha=8'); else plot(x,VV0','c','LineWidth',2); legend('khong co SVC'); end elseif module==5 if A(1,3)>0 xlabel('Ploss'); ylabel('Vref'); VVref=[0.9:0.02:1.04]; title('anh huong Vref len dien ap'); plot(PLL,VVref,'r','LineWidth',1.5); legend('Ploss-Volt'); else plot(x,VV0','c','LineWidth',2); legend('khong co SVC'); end % elseif module==6 if A(1,3)>0 VV(6,:)=[]; title('anh huong bu bang SVC len dien ap'); plot(x,VV(:,1)','r','LineWidth',1.5); Hoàng Mạnh Cường 128 PHỤ LỤC B plot(x,VV(:,2)','b','LineWidth',1.5); plot(x,VV(:,3)','g','LineWidth',1.5); plot(x,VV(:,4)','y','LineWidth',1.5); plot(x,VV(:,5)','k','LineWidth',1.5); legend('Qref=0','Qref=-0.1','Qref=-0.2','Qref=-0.3','Qref=-0.4'); else title('bieu dien ap he thong'); plot(x,VV0','c','LineWidth',2); legend('khong co SVC'); end else end if A(1,3)>0 x=[1 5]; load('PV0_data.mat'); hold on; grid on; xlabel('bus'); ylabel('Voltage'); VV(6,:)=[]; title('khao sat tac dong cua SVC len dien ap o che bu cong suat'); legend('SVC dat tai nut so 3','SVC dat tai nut so 4','SVC dat tai nut so 5','khong co SVC'); plot(x,VV(:,1)','r','LineWidth',1.5); plot(x,VV(:,2)','b','LineWidth',1.5); plot(x,VV(:,3)','g','LineWidth',1.5); plot(x,VV0','y','LineWidth',1.5); else legend('khong co SVC'); plot(x,VV0','c','LineWidth',2); end display_result.m %% DISPLAY RESULT %% header=[ ' ' Hoàng Mạnh Cường 129 PHỤ LỤC B |' |' |' |' ' | ' |De tai: ung dung SVC dieu khien dien ap va bu cong suat ' |GVHD : TS Nguyen Van Liem ' |SVTH : Hoang Manh Cuong ' |MSSV : 41100450 |' Luan an tot nghiep ' | |']; disp(header); fprintf(' \n'); fprintf(' \n'); fprintf(' \n'); fprintf(' DATA SUMMARY\n'); fprintf(' \n'); fprintf('- Base VA =%5.0f [MVA] \n', baseVA); fprintf('- Number Of Bus =%5.0f bus \n', N); fprintf('- Number Of Branch =%5.0f branch \n', nbr); if A(1,3) > fprintf('- Number Of SVC =%5.0f device \n', nSVC); fprintf('- Bus SVC connect =%5.0f \n', hSVC); fprintf('- Bus SVC =%5.0f \n', lSVC); end fprintf('- Bus is Slack bus.\n'); fprintf('- Number Of PQ bus =%5.0f bus.\n', L ); fprintf('- Number Of PV bus =%5.0f bus.\n', G ); fprintf(' \n'); fprintf(' \n'); fprintf(' \n'); % -header1=[ ' -BUS DATA -'Bus type notation: ' [0]: PQ bus Hoàng Mạnh Cường ' ' ' 130 PHỤ LỤC B ' [1]: Slack bus ' ' [2]: PV bus ' ' ']; header2=[ ' ' ' ' ' Bus Bus Usp -Load -Generator Static Mvar' ' No Type (PU.) MW Mvar MW Mvar Qmin Qmax +Qc/-Ql ' ' ' ' ']; disp(header1); disp(header2); for n=1:Nbd fprintf(' %3.0f', busdata(n,1)), fprintf(' %5.0f', busdata(n,2)), fprintf(' %10.3f', busdata(n,3)), fprintf(' %9.3f', busdata(n,4)), fprintf(' %9.3f', busdata(n,5)), fprintf(' %9.3f', busdata(n,6)) fprintf(' %9.3f ', busdata(n,7)), fprintf(' %9.3f', busdata(n,8)), fprintf(' %9.3f ', busdata(n,9)), fprintf(' %9.3f \n', busdata(n,10)); end fprintf(' _\n\n') % header3=[ ' ' ' -BRANCH DATA - ' ' Bus Bus R X B Tr.Tap Setting ' ' nl nr PU PU PU PU ' ' - ']; fprintf('\n'); disp(header3); fprintf('\n'); for n=1:nbr fprintf(' %4.0f', linedata(n,1)), fprintf(' %7.0f', linedata(n,2)), fprintf(' %9.3f', linedata(n,3)), fprintf(' %9.3f', linedata(n,4)), fprintf(' %9.3f', linedata(n,5)), fprintf(' %7.0f \n', linedata(n,6)), end fprintf(' '); fprintf('\n'); Hoàng Mạnh Cường 131 PHỤ LỤC B fprintf('\n') % fprintf(' \n'); fprintf(' \n'); header3=[ ' SVC DATA -' 'stypeSVC=0: control volt bus ' 'stypeSVC=1: control power bus ' ' _' 'lowbus highbus stypeSVC apha Bmax Bmin Vref Qref ' ' -']; if A(1,3) > disp(header3); for k=1:nSVC ISVC(k)=(absV(hSVC(k))-absV(lSVC(k))*cos(theta(lSVC(k))-theta(hSVC(k))))/XB(k); if typeSVC(k)==0 fprintf('%5.0f %5.0f %8.0f %9.2f %8.3f %8.3f %6.3f %8.3f \n',lSVC(k),hSVC(k),typeSVC(k),apha(k) ,Bmax(k),Bmin(k),Vref(k),''); else fprintf('%5.0f %5.0f %8.0f %9.3f %8.3f %8.3f %14.3f %14.3f \n',lSVC(k),hSVC(k),typeSVC(k),apha(k) ,Bmax(k),Bmin(k),'',Qref(k)); end end fprintf(' _'); fprintf(' \n'); fprintf(' \n'); end % -fprintf('\n\n'); fprintf(' -\n'); fprintf(' POWER FLOW SOLUTION BY NEWTON-RAPHSON ALGORITHM\n'); fprintf(' -\n'); fprintf('- Maximum Power Mismatch : %g \n', err); fprintf('- No of Iterations : %g \n', ite); fprintf(' -\n'); Hoàng Mạnh Cường 132 PHỤ LỤC B fprintf('\n\n'); % header4=[ ' -BUS VOLTAGE ' ' _' ' BUS absV theta ' ' -']; disp(header4); theta = theta*180/pi; for i=1:N if A(1,3) > % for k=1:nSVC % if i==hSVC(k) % fprintf(' -\n'); % fprintf('%5.0f %9.3f %9.3f \n',i,absV(i),theta(i)); % fprintf(' -\n'); % end % % end fprintf('%5.0f %9.3f %9.3f \n',i,absV(i),theta(i)); else fprintf('%5.0f %9.3f %9.3f \n',i,absV(i),theta(i)); end end fprintf(' \n') % -fprintf('\n') fprintf(' \n'); header5=[ ' CURRENT FLOW - ' ' AT EACH TRANSMISSION LINE ' ' _' ' Bus Bus Current flow ' ' nl nr nl to nr ' ' -']; disp(header5); Hoàng Mạnh Cường 133 PHỤ LỤC B for L=[1:nbr] fprintf(' %5.0f %4.0f %12.3f %7.3fj',nl(L),nr(L),real(In(L)),imag(In(L))); fprintf(' \n'); end fprintf(' _'); fprintf(' \n'); % -fprintf(' \n'); fprintf(' \n'); header6 =[ ' ' ' TRANSMISSION LINE FLOW AND LOSS ' ' ' ' Bus Bus Power Flow Power Flow -Power Loss - ' ' nl nr nl to nr nr to nl MW Mvar ' ' ']; disp(header6); for L = [1:nbr] fprintf(' %5.0f', nl(L)), fprintf(' %5.0f', nr(L)), fprintf(' %9.3f %7.3fj', real(Snk(L)),imag(Snk(L))), fprintf(' %9.3f %7.3fj', real(Skn(L)),imag(Skn(L))), fprintf(' %9.3f %7.3fj', real(SL(L)),imag(SL(L))), fprintf(' \n'); end fprintf(' \n'); fprintf(' Total ') fprintf(' %9.3f %7.3fj', sum(real(Snk)),sum(imag(Snk))), fprintf(' %9.3f %7.3fj', sum(real(Skn)),sum(imag(Skn))), fprintf(' %9.3f', sum(real(SL))), fprintf(' %7.3fj\n', sum(imag(SL))) fprintf(' '); fprintf(' \n'); fprintf(' \n'); fprintf(' \n'); fprintf(' \n'); %% header7=[ ' SVC RESULT ' Hoàng Mạnh Cường 134 PHỤ LỤC B 'lSVC hSVC Bmin Bcal Bmax QSVC ISVC delVmin delVmax ' ' -']; if A(1,3) > disp(header7); limit_voltage_control; for k=1:nSVC ISVC(k)=(absV(hSVC(k))-absV(lSVC(k))*cos(theta(lSVC(k))-theta(hSVC(k))))/XB(k); fprintf('%4.0f %5.0f %9.3f %9.3f %9.3f %9.3f %9.3f %9.3f %12.3f \n',lSVC(k),hSVC(k),Bmin(k),Bcal(k),Bmax(k) ,imag(Snk(nbr-k+1)),ISVC(k),delVmin(k),delVmax(k)); end fprintf(' _'); fprintf(' \n'); fprintf(' \n'); end %% Hoàng Mạnh Cường 135 PHỤ LỤC B TÀI LIỆU THAM KHẢO [1] E S ABB-AB, "SVC Static Var Compensator -An insurance for improved grid system stability and reliability," ABB – the pioneer 2010-11 [2] H a L Gyugyi, Understanding FACTS: concepts and technology of flexible ac transmission systems IEEE Press, 1999 [3] H Tyll, R K Varma, H Bilodeau, and C Horwill, "Tutorial on Static Var Compensators," San Francisco June 12, 2005 [4] SamitShah, "static var control (SVC) project at baymon t.c.inpureto rico," IEEE PES PSCE, MARCH 2009 2009 [5] B t s o F Didactic, "Electricity and New Energy Static Var Compensator (SVC," vol 1, 01/2015 [6] M A U Rahman and M S u Islam (2014, Voltage Control and Dynamic Performance of Power Transmission Using Static VAR Compensator Vol [7] P Kundur, N J Balu, and M G Lauby, Power system stability and control vol 7: McGraw-hill New York, 1994 [8] K Prabha, Power system stability and control, 1994 [9] C R F.-E Enrique Acha, Hugo AMbriz-perez,Cesar Andeles-Camacho, FACTS modelling and Simulation in Power Networks, 2004 [10] L Sahu, "Modelling of STACOM and SVC for power sytem steady state operation and anhancement of transient stability of a multi-machine power sytem by STACOM " Hoàng Mạnh Cường 136 ... nguyên lý điều khiển điện áp bù công suất phản kháng SVC - Kiểm chứng khả điều khiển điện áp SVC tất nút tải hệ thống điện đơn giản nút giản đồ điện áp - Kiểm chứng khả bù công suất nút mạng điện. .. đến điện áp nút hệ thống Công suất phản kháng khu vực q thừa xảy tượng điện áp (điện áp cao) Ngược lại, thiếu công suất phản kháng điện áp bị sụt xuống thấp Do đó, giống cơng suất tác dụng, công. .. động SVC sau (Bảng 2.1): Bảng 2.1 Tóm tắt mơ hình điều khiển SVC Chế độ điều khiển Điều khiển điện áp Phương trình thứ f 1svc = Vh − Vref − α Điều khiển dòng cống suất phản kháng Ngồi vùng điều khiển