(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy

(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy(Luận văn thạc sĩ) Ổn định hệ thống điện cho hệ nhiều máy

LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu Các số liệu, kết nêu luận văn trung thực chƣa đƣợc cơng bố cơng trình khác Tp Hồ Chí Minh, ngày … tháng … năm 201… Họ tên Nguyễn Thị Phƣơng Quỳnh ii LỜI CẢM TẠ Lời đầu tiên, tác giả luận văn xin bày tỏ lòng biết ơn chân thành đến tất thầy cô giáo Khoa Điện – Điện Tử , Trƣờng Đại Học Sƣ Phạm Kỹ Thuật TP HCM truyền đạt cho tác giả kinh nghiệm kiến thức thật quý báu thời gian học tập trƣờng Đặc biệt xin chân thành cám ơn thầy PGS TS Quyền Huy Ánh, ngƣời trực tiếp hƣớng dẫn tác giả hoàn thành luận văn thạc sĩ Xin đƣợc gởi lời cám ơn chân thành đến gia đình, bạn bè ngƣời thân, ngƣời giúp đỡ động viên tác giả tinh thần nhƣ vật chất suốt thời gian thực luận văn Trong trình thực luận văn, thời gian kinh nghiệm có hạn nên chắn khơng tránh đƣợc thiếu xót, mong đƣợc quan tâm đóng góp thầy bạn bè, anh chị để luận văn đƣợc hoàn chỉnh Mọi đóng góp luận văn xin gởi qua địa email phuongquynhcdn@gmail.com Xin chân thành cảm ơn ! Tp Hồ Chí Minh, tháng 01/2013 Học viên thực Nguyễn Thị Phƣơng Quỳnh iii TÓM TẮT Hệ thống điện Việt Nam năm qua có bƣớc phát triển nhảy vọt công suất quy mô lãnh thổ Ngành công nghiệp hệ thống điện lĩnh vực có thay đổi liên tục ngày trở nên phức tạp chúng liên kết với Do đó, vấn đề đặt phải đảm bảo ổn định cho toàn hệ thống điện Chính vậy, việc nghiên cứu để tìm phƣơng pháp đánh giá ổn định nhanh xác thực trở nên cần thiết hệ thống điện Việt Nam Mục đích đề tài khảo sát đánh giá tính ổn định hệ thống điện cho hệ nhiều máy phát phƣơng pháp mô Để đạt đƣợc mục đích ngƣời nghiên cứu tiến hành công việc: - Xây dựng phƣơng pháp đánh giá ổn định hệ thống điện - Xây dựng chƣơng trình tính tốn phân bố cơng suất hệ thống điện cho hệ nhiều máy phần mềm Matlab - Xây dựng mơ hình mơ ngơn ngữ lập trình m.file có phần mềm Matlab - Mơ đánh giá kết đạt đƣợc từ việc đánh giá tính ổn định hệ thống điện nút máy phát Những kết mà đề tài đạt đƣợc là: - Phƣơng pháp đánh giá tính ổn định hệ thống điện - Cơng cụ đánh giá ổn định hệ thống nút máy phát Có thể kết luận rằng, phƣơng pháp cơng cụ đánh giá tính ổn định hệ thống điện cho hệ nhiều máy phƣơng pháp mô góp phần vào việc nghiên cứu tìm phƣơng pháp đánh giá xác, nhanh chóng trực quan để ứng dụng vào thực tiễn hệ thống điện Việt Nam iv ABSTRACT Vietnam's power system in recent years has seen dramatic development in power and territorial size Industrial power system is a field with constantly changing and becoming more complex as they relate to each other Therefore, the problem posed is to ensure the stability of the entire electrical system Therefore, the research to find methods to evaluate stability quickly and accurately becomes really necessary for Vietnam's power system The purpose of thesis is set up the method and tools to evaluate the dynamic stability of power system by simulation To get this objective, the researcher have to carry out following steps: - Set up the method to evaluate the dynamic stability of power system - Set up the calculated equations to distribute the power of power system with many generators by Matlab software - Set up the simulation model by programming language m.file in Matlab software - Simulation and accessing the result from bus, generators The result of thesis reached: - The new method to evaluate dynamic stability of power system - Tools to evaluate the dynamic stability of bus, generators system It could be said that the algorithm and tools to evaluate the dynamic stability of power system through simulation method is contributed to find the method which has accurately, quickly and visually This new method is able to apply in the power system in Vietnam v MỤC LỤC Trang tựa Quyết định giao đề tài Lý lịch khoa học i Lời cam đoan ii Lời cảm tạ iii Tóm tắt iv Abstract………………………………………………………….…………….v Mục lục vi Danh sách hình ix Danh sách bảng xii Chƣơng 1: TỔNG QUAN 1.1 Tổng quan lĩnh vực nghiên cứu 1.1.1 Kết nghiên cứu nƣớc 1.1.2 Kết nghiên cứu nƣớc 1.2 Mục đích đề tài 1.3 Nhiệm vụ đề tài giới hạn đề tài 1.3.1 Nhiệm vụ đề tài 1.3.2 Giới hạn đề tài 1.4 Phƣơng pháp nghiên cứu Chƣơng 2: CƠ SỞ LÝ THUYẾT 2.1 Các chế độ hệ thống điện………………………………………………… 2.2 Các khái niệm ổn định hệ thống điện………………………… 2.2.1 Ổn định tĩnh hệ thống điện………………………………… 2.2.2 Ổn định động hệ thống điện 11 2.3 Sơ lƣợc ứng dụng phƣơng pháp toán học nghiên cứu ổn định hệ thống điện………………………………………………………………… .13 vi 2.3.1 Phƣơng pháp phân tích đánh giá ổn định theo tiêu chuẩn lƣợng …………………………………………………….………… 14 2.3.2 Phƣơng pháp phân tích đánh giá ổn định theo Lyapunov …… 16 2.3.2.1 Phƣơng pháp trực tiếp (hay gọi phƣơng pháp thứ hai Lyapunov)…… ……………………………………………… 20 2.3.2.2 Phƣơng pháp xấp xỉ bậc Lyapunov (hay gọi phƣơng pháp thứ Lyapunov)………………… …………22 2.4 Hậu ổn định yêu cầu đảm bảo ổn định cho hệ thống điện …………………………… …………………………………………30 2.5 Kết luận ……………………… ………………………………………… 31 Chƣơng 3: MƠ HÌNH TỐN 3.1 Các giả thuyết………………… ………………………………………… 33 3.2 Giải thuật phân tích đánh giá ổn định hệ nhiều máy phƣơng pháp mô phỏng……………… …………………………………………………….34 3.2.1 Phƣơng pháp phân bố công suất phép lặp Gauss – Seidel… 36 3.2.2 Phƣơng pháp phân bố công suất phép lặp Newton – Grapson 40 3.3 Thời gian cắt tới hạn 50 3.4 Kết luận 51 Chƣơng 4: XÂY DỰNG MƠ HÌNH TỐN 4.1 Cơng suất điện máy phát……………………….……… 52 4.1.1 Công suất điện máy phát thời điểm xảy cố………………………………………………………………….……… 52 4.1.2 Công suất điện máy phát thời điểm sau cố xảy ra……………………………………………………………….…… 53 4.2 Trạng thái máy phát………………………….……………………53 4.2.1 Trạng thái máy phát thời điểm cố xảy ra………………………………………………………………… 53 vii 4.2.2 Trạng thái máy phát sau cố xảy ra………………………………………………………………… 56 4.2.3 Trạng thái máy phát tồn q trình có cố xảy ra…………………………………………………….………… 57 4.3 Kết luận 57 Chƣơng 5: KẾT QUẢ MÔ PHỎNG 5.1 Hệ thống điện máy – nút………………………………………………58 5.1.1 Tính tốn thơng số đầu vào phƣơng pháp giải tích…… …… 60 5.1.2 Khảo sát đánh giá ổn định hệ thống điện máy – nút………….63 Chƣơng 6: KẾT LUẬN VÀ HƢỚNG PHÁT TRIỂN ĐỀ TÀI 6.1 Tổng kết …………………… ……………………………………………83 6.2 Những mặt tồn hƣớng phát triển cho đề tài……………………84 TÀI LIỆU THAM KHẢO PHỤ LỤC viii DANH SÁCH CÁC HÌNH HÌNH TRANG Hình 2.1: Sơ đồ hệ thống điện đơn giản Hình 2.2: Đặc tính cơng suất điện từ máy phát đặc tính cơng suất tuabin Hình 2.3: Sơ đồ hệ thống điện xét nút phụ tải tƣơng quan cân công suất phản kháng 10 Hình 2.4: Đặc tính Q _ V 11 Hình 2.5: Hệ thống ổn định động 12 Hình 2.6: Hệ thống ổn định động 12 Hình 2.7: Ổn định dao động bé theo Lyapunov 18 Hình 2.8: Tiêu chuẩn ổn định Mikhailov 28 Hình 2.9a: Quỹ đạo ổn định hệ thống 30 Hình 2.9 b, c: Quỹ đạo không ổn định hệ thống 30 Hình 3.1: Lƣu đồ giải thuật phân tích ổn định hệ nhiều máy 35 Hình 3.2: Sơ đồ thay hình  39 Hình 3.3: Sơ đồ thay tƣơng đƣơng hệ thống điện nút 40 Hình 3.4: Sơ đồ thay hệ thống điện nhiều máy 46 Hình 3.5: Lƣu đồ thuật tốn xác định thời gian cắt tới hạn 51 Hình 5.1: Sơ đồ hệ thống điện nút – máy 58 Hình 5.2: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tc = 0.1s) 64 Hình 5.3: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tc = 0.16s) 64 Hình 5.4: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tc = 0.17s) 65 Hình 5.5: Vị trí góc pha máy 65 ix Hình 5.6: Mối quan hệ vận tốc góc w21 w31 (tc = 0.1s) 66 Hình 5.7: Mối quan hệ vận tốc góc w21 w31 (tc = 0.16s) 66 Hình 5.8: Mối quan hệ vận tốc góc w21 w31 (tc = 0.17s) 67 Hình 5.9: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 10%) 69 Hình 5.10: Vị trí góc pha máy (tăng tải 10%) 69 Hình 5.11: Mối quan hệ vận tốc góc w21 w31 (tăng tải 10%) 70 Hình 5.12: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 20%) 70 Hình 5.13: Vị trí góc pha máy (tăng tải 20%) 71 Hình 5.14: Mối quan hệ vận tốc góc w21 w31 (tăng tải 20%) 71 Hình 5.15: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 30%) 72 Hình 5.16: Vị trí góc pha máy (tăng tải 30%) 72 Hình 5.17: Mối quan hệ vận tốc góc w21 w31 (tăng tải 30%) 73 Hình 5.18: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 40%) 73 Hình 5.19: Vị trí góc pha máy (tăng tải 40%) 74 Hình 5.20: Mối quan hệ vận tốc góc w21 w31 (tăng tải 40%) 74 Hình 5.21: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 50%) 75 Hình 5.22: Vị trí góc pha máy (tăng tải 50%) 75 Hình 5.23: Mối quan hệ vận tốc góc w21 w31 (tăng tải 50%) 76 Hình 5.24: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 60%) 76 Hình 5.25: Vị trí góc pha máy (tăng tải 60%) 77 Hình 5.26: Mối quan hệ vận tốc góc w21 w31 (tăng tải 60%) 77 Hình 5.27: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 70%) 78 x Hình 5.28: Vị trí góc pha máy (tăng tải 70%) 78 Hình 5.29: Mối quan hệ vận tốc góc w21 w31 (tăng tải 70%) 79 Hình 5.30: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 80%) 79 Hình 5.31: Vị trí góc pha máy (tăng tải 80%) 80 Hình 5.32: Mối quan hệ vận tốc góc w21 w31 (tăng tải 80%) 80 Hình 5.33: Góc lệch pha máy phát máy phát so với máy phát chuẩn (tăng tải 90%) 81 Hình 5.34: Vị trí góc pha máy (tăng tải 90%) 81 Hình 5.35: Mối quan hệ vận tốc góc w21 w31 (tăng tải 90%) 82 xi  Tiếp tục xây dựng phát triển chương trình tính tốn Matlab, tạo giao diện thân thiện gần gũi với người sử dụng  Để nâng cao tính ổn định hệ thống người ta sử dụng AVR, PSS hay kết hợp AVR PSS 85 TÀI LIỆU THAM KHẢO Tài liệu tham khảo tiếng việt: [1] Nguyễn Hoàng Việt, Ngắn mạch ổn định hệ thống điện, Nhà xuất Đại học Quốc Gia TP.HCM, 2010 [2] PGS.TS Quyền Huy Ánh, Mơ hình hố mơ phỏng, Đại học Sư Phạm Kỹ Thuật TP Hồ Chí Minh [3] GS.TS Nguyễn Cơng Hiền, Mơ hình hóa hệ thống mô phỏng, Đại Học Bách Khoa Hà Nội [4] Lã Văn Út, Phân tích điều khiển ổn định hệ thống điện , Nhà xuất Khoa Học Kỹ thuật, 2000 [5] Nguyễn Phùng Quang, Matlab and Simulink, Nhà xuất KhoaHọc Kỹ Thuật, 2005 Tài liệu tham khảo tiếng anh: [6] D.P.Kothari, T.S Bhatti, Ramnarayan Patel, Matlab/Simulink – based transient stability analysis of a multimachine power system, Centre for Energy Studies, India, 1999 [7] P.Kundur, Power System Stability And Control, Mc.Graw – Hill International Editions, 1999 [8] D.P.Kothari, I.J.Nagrath, Modern Power System Analysis, The McGraw – Hill Companies, 2003 [9] Simulink, Dynamic System Simulation for Matlab, MathWorks Inc.,1998 [10] M.A.Pai and Peter W.Sauer, Stability Analysis of Power System by Lyapunov ‘s Direct Method, IEEE Control Systems, 1989 86 [11] Tharam S Dillon, Dagmar Biebur, Multimachine Eigenvalue Sensitivities of Power Systems Parameters, IEEE Transaction on Power System, 2000 [12] Anthony N Michel, A A Fouad and Vijay Vittal, Power System Transient Stability Using Individual Machine Energy Functions, IEEE Transaction on Power System, 1983 [13] Hadi Saadat, Power System Analysis, Mc.Graw – Hill International Editions, 1999 87 PHỤ LỤC % PHÂN TÍCH ĐỘ ỔN ĐỊNH QUÁ ĐỘ CỦA HỆ NHIỀU MÁY TRONG HỆ THỐNG ĐIỆN global Pm f H E Y th ngg f=60; % zdd=gendata(:,2)+j*gendata(:,3); ngr=gendata(:,1); % H=gendata(:,4); ngg=length(gendata(:,1)); %% for k=1:ngg zdd(ngr(k))=gendata(k, 2)+j*gendata(k,3); % H(ngr(k))=gendata(k, 4); H(k)=gendata(k,4); % new end %% for k=1:ngg I=conj(S(ngr(k)))/conj(V(ngr(k))); % Ep(ngr(k)) = V(ngr(k))+zdd(ngr(k))*I; % Pm(ngr(k))=real(S(ngr(k))); Ep(k) = V(ngr(k))+zdd(ngr(k))*I; % new Pm(k)=real(S(ngr(k))); % new end E=abs(Ep); d0=angle(Ep); for k=1:ngg nl(nbr+k) = nbus+k; nr(nbr+k) = gendata(k, 1); %R(nbr+k) = gendata(k, 2); %X(nbr+k) = gendata(k, 3); R(nbr+k) = real(zdd(ngr(k))); X(nbr+k) = imag(zdd(ngr(k))); Bc(nbr+k) = 0; a(nbr+k) = 1.0; yload(nbus+k)=0; end nbr1=nbr; nbus1=nbus; nbrt=nbr+ngg; nbust=nbus+ngg; 88 linedata=[nl, nr, R, X, -j*Bc, a]; [Ybus, Ybf]=ybusbf(linedata, yload, nbus1,nbust); fprintf('\Ma trận tổng dẫn nút rút gọn trước cố\n') Ybf Y=abs(Ybf); th=angle(Ybf); Pm=zeros(1, ngg); disp([' G(i) E''(i) d0(i) Pm(i)']) for ii = 1:ngg for jj = 1:ngg Pm(ii) = Pm(ii) + E(ii)*E(jj)*Y(ii, jj)*cos(th(ii, jj)-d0(ii)+d0(jj)); end, fprintf(' %g', ngr(ii)), fprintf(' %8.4f',E(ii)), fprintf(' %8.4f', 180/pi*d0(ii)) fprintf(' %8.4f \n',Pm(ii)) end respfl='y'; while respfl =='y' | respfl=='Y' nf=input('Nhập vào nút bị cố -> '); fprintf('\Ma trận tổng dẫn nút rút gọn xảy cố\n') Ydf=ybusdf(Ybus, nbus1, nbust, nf) %Fault cleared [Yaf]=ybusaf(linedata, yload, nbus1,nbust, nbrt); fprintf('\nMa trận tổng dẫn nút rút gọn sau cố\n') Yaf resptc='y'; while resptc =='y' | resptc=='Y' tc=input('Nhập vào thời gian cắt cố tính giây, tc = '); tf=input('Nhập vào thời gian mơ cuối tính giây, tf = '); clear t x del t0 = 0; w0=zeros(1, length(d0)); x0 = [d0, w0]; tol=0.0001; Y=abs(Ydf); th=angle(Ydf); %[t1, xf] =ode23('dfpek', t0, tc, x0, tol); % Trong trình xảy cố (sử dụng MATLAB 4) tspan=[t0, tc]; %sử dụng MATLAB [t1, xf] =ode23('dfpek', tspan, x0); % Trong trình xảy cố, sử dụng MATLAB 5) x0c =xf(length(xf), :); Y=abs(Yaf); th=angle(Yaf); %[t2,xc] =ode23('afpek', tc, tf, x0c, tol); % Sau cố(sử dụng MATLAB 4) tspan = [tc, tf]; % sử dụng MATLAB [t2,xc] =ode23('afpek', tspan, x0c); % Quá trình sau cố (sử dụng MATLAB 5) 89 t =[t1; t2]; x = [xf; xc]; fprintf('\n cố cắt %4.3f Sec \n', tc) for k=1:nbus if kb(k)==1 ms=k; else, end end fprintf('\Góc lệch pha máy phát\n') fprintf(‘độ máy phát chuẩn.\n') fprintf(' t - sec') kk=0; for k=1:ngg if k~=ms kk=kk+1; del(:,kk)=180/pi*(x(:,k)-x(:,ms)); fprintf(' d(%g,',ngr(k)), fprintf('%g)', ngr(ms)) else, end end fprintf(' \n') disp([t, del]) h=figure; figure(h) plot(t, del) title([‘Góc lệch pha máy (cắt cố tại', num2str(tc),'s)']) xlabel('t, sec'), ylabel('Delta, degree'), grid resp=0; while strcmp(resp, 'n')~=1 & strcmp(resp, 'N')~=1 & strcmp(resp, 'y')~=1 & strcmp(resp, 'Y')~=1 resp=input(‘ nhập vào thời gian cắt cố? Nhập ''y'' or ''n'' ngoặc kép -> '); if strcmp(resp, 'n')~=1 & strcmp(resp, 'N')~=1 & strcmp(resp, 'y')~=1 & strcmp(resp, 'Y')~=1 fprintf('\n Khơng cố gắng thử lại \n\n'), end end resptc=resp; end resp2=0; while strcmp(resp2, 'n')~=1 & strcmp(resp2, 'N')~=1 & strcmp(resp2, 'y')~=1 & strcmp(resp2, 'Y')~=1 resp2=input('Vị trí xảy cố khác: Nhập vào ''y'' or ''n'' dấu ngoặc kép -> '); if strcmp(resp2, 'n')~=1 & strcmp(resp2, 'N')~=1 & strcmp(resp2, 'y')~=1 & strcmp(resp2, 'Y')~=1 fprintf('\n Nếu không đúng, cố gắng thử lại \n\n'), end respf1=resp2; 90 end if respf1=='n' | respf1=='N', return, else, end end THÔNG SỐ ĐẦU VÀO basemva = 100; accuracy = 0.0001; maxiter = 10; busdata=[1 1.03 0.0 00.000 0.00 71.62 00.00 0 2 1.02 0.0 00.000 0.00 162.98 00.00 150 1.02 0.0 00.000 0.00 85.04 00.00 100 1.02 0.0 0.0000 0.00 00.00 00.00 0 1.02 0.0 125.00 50.00 00.00 00.00 0 1.02 0.0 90.000 30.00 00.00 00.00 0 1.02 0.0 0.0000 0.000 00.00 00.00 0 1.02 0.0 100.00 35.00 00.00 00.00 0 1.02 0.0 0.000 0.000 00.00 00.00 0 0; linedata=[1 0.000 0.0576 0.00 1.0 0.000 0.0625 0.00 1.0 0.000 0.0586 0.0 1.0 0.010 0.085 0.088 1.0 0.017 0.092 0.079 1.0 0.032 0.161 0.153 1.0 0.039 0.170 0.179 1.0 0.0085 0.072 0.0745 1.0 0.0119 0.1008 0.1045 1.0]; lfybus lfnewton busout % Gen gendata=[ trstab % ma trận tổng dẫn nút % Tính dịng cơng suất phương pháp Newton-Raphson % In dịng cơng suất Ra Xd' 0.0608 0.1198 0.1813 H 23.64 6.4 3.01]; 91 VẬN TỐC GÓC for n=1:(sl-1) tocdo(n,:)=(del(n+1,:)-del(n,:))/(t(n+1)-t(n)); tocdo1(n,:)= (tocdo (n,:)*pi)/180; end tocdo1(sl,:)=tocdo1(sl-1,:); plot(t,tocdo1) VỊ TRÍ GĨC PHA CỦA TỪNG MÁY syms tam1 tam1 = [x(:,1), x(:,2), x(:,3)]*60; plot (t, tam1) 92 NGƢNG LAI KHƠNG IN ………………………………………… PHỤ LỤC % PHÂN TÍCH ĐỘ ỔN ĐỊNH QUÁ ĐỘ CỦA HỆ NHIỀU MÁY TRONG HỆ THỐNG ĐIỆN % Chương trình ``trstab'' sử dụng kết hợp với chương trình dịng cơng suất Bất kỳ dịng công suất ``lfgauss'', ``lfnewton'', ``decoupled (ngắt liên kết)'' ``perturb (nhiễu loạn)'' sử dụng ưu tiên chương trình ``trstab'' Chương trình dịng cơng suất gồm cơng suất, độ lớn điện áp góc pha cho nút Vì vậy, tổng dẫn tải tính tốn lại ma trận có tên ``yload'' Thêm vào đó, liệu cung cấp thêm dịng cơng suất, điện kháng độ số quán tính máy xác lập Điều xác định ma trận tên ``gendata'', Mỗi hàng chứa số nút kết nối với máy phát, điện trở phần ứng, điện kháng độ, số quán tính máy % Chương trình ``trstab'' chứa ma trận tổng dẫn nút trước cố bao gồm tổng dẫn tải Điện áp phía sau điện kháng độ % Ma trận tổng dẫn trước, sau cố tìm thấy `Phương trình máy phát phương trình biến trạng thái ode23 Matlab sử dụng để giải phương trình hệ nhiều máy Các góc lệch pha máy phát so với nút chuẩn hiển thị hình vẽ % Việc mơ lặp lại với nhiều thời gian cắt cố khác vị trí cắt khác % global Pm f H E Y th ngg f=60; % zdd=gendata(:,2)+j*gendata(:,3); ngr=gendata(:,1); % H=gendata(:,4); ngg=length(gendata(:,1)); %% for k=1:ngg zdd(ngr(k))=gendata(k, 2)+j*gendata(k,3); % H(ngr(k))=gendata(k, 4); H(k)=gendata(k,4); % new end %% for k=1:ngg 93 I=conj(S(ngr(k)))/conj(V(ngr(k))); % Ep(ngr(k)) = V(ngr(k))+zdd(ngr(k))*I; % Pm(ngr(k))=real(S(ngr(k))); Ep(k) = V(ngr(k))+zdd(ngr(k))*I; % new Pm(k)=real(S(ngr(k))); % new end E=abs(Ep); d0=angle(Ep); for k=1:ngg nl(nbr+k) = nbus+k; nr(nbr+k) = gendata(k, 1); %R(nbr+k) = gendata(k, 2); %X(nbr+k) = gendata(k, 3); R(nbr+k) = real(zdd(ngr(k))); X(nbr+k) = imag(zdd(ngr(k))); Bc(nbr+k) = 0; a(nbr+k) = 1.0; yload(nbus+k)=0; end nbr1=nbr; nbus1=nbus; nbrt=nbr+ngg; nbust=nbus+ngg; linedata=[nl, nr, R, X, -j*Bc, a]; [Ybus, Ybf]=ybusbf(linedata, yload, nbus1,nbust); fprintf('\Ma trận tổng dẫn nút rút gọn trước cố\n') Ybf Y=abs(Ybf); th=angle(Ybf); Pm=zeros(1, ngg); disp([' G(i) E''(i) d0(i) Pm(i)']) for ii = 1:ngg for jj = 1:ngg Pm(ii) = Pm(ii) + E(ii)*E(jj)*Y(ii, jj)*cos(th(ii, jj)-d0(ii)+d0(jj)); end, fprintf(' %g', ngr(ii)), fprintf(' %8.4f',E(ii)), fprintf(' %8.4f', 180/pi*d0(ii)) fprintf(' %8.4f \n',Pm(ii)) end respfl='y'; while respfl =='y' | respfl=='Y' nf=input('Nhập vào nút bị cố -> '); 94 fprintf('\Ma trận tổng dẫn nút rút gọn xảy cố\n') Ydf=ybusdf(Ybus, nbus1, nbust, nf) %Fault cleared [Yaf]=ybusaf(linedata, yload, nbus1,nbust, nbrt); fprintf('\nMa trận tổng dẫn nút rút gọn sau cố\n') Yaf resptc='y'; while resptc =='y' | resptc=='Y' tc=input('Nhập vào thời gian cắt cố tính giây, tc = '); tf=input('Nhập vào thời gian mô cuối tính giây, tf = '); clear t x del t0 = 0; w0=zeros(1, length(d0)); x0 = [d0, w0]; tol=0.0001; Y=abs(Ydf); th=angle(Ydf); %[t1, xf] =ode23('dfpek', t0, tc, x0, tol); % Trong trình xảy cố (sử dụng MATLAB 4) tspan=[t0, tc]; %sử dụng MATLAB [t1, xf] =ode23('dfpek', tspan, x0); % Trong trình xảy cố, sử dụng MATLAB 5) x0c =xf(length(xf), :); Y=abs(Yaf); th=angle(Yaf); %[t2,xc] =ode23('afpek', tc, tf, x0c, tol); % Sau cố(sử dụng MATLAB 4) tspan = [tc, tf]; % sử dụng MATLAB [t2,xc] =ode23('afpek', tspan, x0c); % Quá trình sau cố (sử dụng MATLAB 5) t =[t1; t2]; x = [xf; xc]; fprintf('\n cố cắt %4.3f Sec \n', tc) for k=1:nbus if kb(k)==1 ms=k; else, end end fprintf('\Góc lệch pha máy phát\n') fprintf(‘độ máy phát chuẩn.\n') fprintf(' t - sec') kk=0; for k=1:ngg if k~=ms kk=kk+1; del(:,kk)=180/pi*(x(:,k)-x(:,ms)); fprintf(' d(%g,',ngr(k)), fprintf('%g)', ngr(ms)) else, end 95 end fprintf(' \n') disp([t, del]) h=figure; figure(h) plot(t, del) title([‘Góc lệch pha máy (cắt cố tại', num2str(tc),'s)']) xlabel('t, sec'), ylabel('Delta, degree'), grid resp=0; while strcmp(resp, 'n')~=1 & strcmp(resp, 'N')~=1 & strcmp(resp, 'y')~=1 & strcmp(resp, 'Y')~=1 resp=input(‘ nhập vào thời gian cắt cố? Nhập ''y'' or ''n'' ngoặc kép -> '); if strcmp(resp, 'n')~=1 & strcmp(resp, 'N')~=1 & strcmp(resp, 'y')~=1 & strcmp(resp, 'Y')~=1 fprintf('\n Không cố gắng thử lại \n\n'), end end resptc=resp; end resp2=0; while strcmp(resp2, 'n')~=1 & strcmp(resp2, 'N')~=1 & strcmp(resp2, 'y')~=1 & strcmp(resp2, 'Y')~=1 resp2=input('Vị trí xảy cố khác: Nhập vào ''y'' or ''n'' dấu ngoặc kép -> '); if strcmp(resp2, 'n')~=1 & strcmp(resp2, 'N')~=1 & strcmp(resp2, 'y')~=1 & strcmp(resp2, 'Y')~=1 fprintf('\n Nếu không đúng, cố gắng thử lại \n\n'), end respf1=resp2; end if respf1=='n' | respf1=='N', return, else, end end 96 THÔNG SỐ ĐẦU VÀO basemva = 100; accuracy = 0.0001; maxiter = 10; busdata=[1 1.03 0.0 00.000 0.00 71.62 00.00 0 2 1.02 0.0 00.000 0.00 162.98 00.00 150 1.02 0.0 00.000 0.00 85.04 00.00 100 1.02 0.0 0.0000 0.00 00.00 00.00 0 1.02 0.0 125.00 50.00 00.00 00.00 0 1.02 0.0 90.000 30.00 00.00 00.00 0 1.02 0.0 0.0000 0.000 00.00 00.00 0 1.02 0.0 100.00 35.00 00.00 00.00 0 1.02 0.0 0.000 0.000 00.00 00.00 0 0; linedata=[1 0.000 0.0576 0.00 1.0 0.000 0.0625 0.00 1.0 0.000 0.0586 0.0 1.0 0.010 0.085 0.088 1.0 0.017 0.092 0.079 1.0 0.032 0.161 0.153 1.0 0.039 0.170 0.179 1.0 0.0085 0.072 0.0745 1.0 0.0119 0.1008 0.1045 1.0]; lfybus lfnewton busout % Gen gendata=[ trstab % ma trận tổng dẫn nút % Tính dịng cơng suất phương pháp Newton-Raphson % In dịng công suất Ra Xd' 0.0608 0.1198 0.1813 H 23.64 6.4 3.01]; 97 VẬN TỐC GÓC for n=1:(sl-1) tocdo(n,:)=(del(n+1,:)-del(n,:))/(t(n+1)-t(n)); tocdo1(n,:)= (tocdo (n,:)*pi)/180; end tocdo1(sl,:)=tocdo1(sl-1,:); plot(t,tocdo1) VỊ TRÍ GĨC PHA CỦA TỪNG MÁY syms tam1 tam1 = [x(:,1), x(:,2), x(:,3)]*60; plot (t, tam1) 98 ... đánh giá ổn định hệ thống điện nhiều máy phát  Xây dựng phương pháp đánh giá ổn định hệ thống điện nhiều máy phương pháp mô  Xây dựng công cụ đánh giá ổn định hệ thống điện nhiều máy phương... đến tính chất ổn định động hệ thống điện, đặc biệt hệ thống điện nhiều máy phương pháp mơ 32 Chƣơng MƠ HÌNH TỐN Phân tích ổn định động cho hệ thống điện đặc biệt hệ thống điện nhiều máy tốn khó... tích ổn định cho hệ thống điện nhiều máy hệ thống điện nhiều máy biến đổi tương đương thành hệ thống có máy nối với hệ thống vô lớn Trong chương tác giả trình bày phương pháp phân tích ổn định cho