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

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Ĩ NGUYỄN THỊ PHƯƠNG QUỲNH ỔN ĐỊNH HỆ THỐNG ĐIỆN CHO HỆ NHIỀU MÁY S K C 0 9 NGÀNH: THIẾT BỊ MẠNG VÀ NHÀ MÁY ĐIỆN - 605250 S KC 0 Tp Hồ Chí Minh, 2012 Luan van 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Ĩ NGUYỄN THỊ PHƯƠNG QUỲNH ỔN ĐỊNH HỆ THỐNG ĐIỆN CHO HỆ NHIỀU MÁY NGÀNH: THIẾT BỊ, MẠNG VÀ NHÀ MÁY ĐIỆN - 605250 Tp.Hồ Chí Minh, tháng 1/2013 Luan van 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Ĩ NGUYỄN THỊ PHƯƠNG QUỲNH ỔN ĐỊNH HỆ THỐNG ĐIỆN CHO HỆ NHIỀU MÁY NGÀNH: THIẾT BỊ, MẠNG VÀ NHÀ MÁY ĐIỆN - 605250 Hướng dẫn khoa học: PGS.TS QUYỀN HUY ÁNH Tp.Hồ Chí Minh, tháng 1/2013 Luan van LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƢỢC: Họ & tên: NGUYỄN THỊ PHƢƠNG QUỲNH Giới tính: Nữ Ngày, tháng, năm sinh: 14-12-1985 Nơi sinh: TP Hồ Chí Minh Quê quán: TP Hồ Chí Minh Dân tộc: Kinh Chỗ riêng địa liên lạc: 99F, Đƣờng 379, P Tăng Nhơn Phú A, Quận Điện thoại quan: Điện thoại nhà riêng: 0973432949 Fax: E-mail: phuongquynhcdn@gmail.com II QUÁ TRÌNH ĐÀO TẠO: Trung học chuyên nghiệp: Hệ đào tạo: Thời gian đào tạo từ ……/…… đến ……/ …… Nơi học (trƣờng, thành phố): Ngành học: Đại học: Hệ đào tạo: Chính quy Thời gian đào tạo từ: 09/2004 đến 05/2009 Nơi học (trƣờng, thành phố): Đại Học Sƣ Phạm Kỹ Thuật TP Hồ Chí Minh Ngành học: Kỹ Thuật Điện-Điện Tử Tên đồ án, luận án mơn thi tốt nghiệp: Điều khiển lập trình nâng cao, Thiết kế hệ thống điện, Quản lý dự án Ngày & nơi bảo vệ đồ án, luận án thi tốt nghiệp: 03/2009, Đại Học Sƣ Phạm Kỹ Thuật TP Hồ Chí Minh Ngƣời hƣớng dẫn: III Q TRÌNH CƠNG TÁC CHUN MƠN KỂ TỪ KHI TỐT NGHIỆP ĐẠI HỌC: Nơi công tác Thời gian Từ 5/2009 đến Công việc đảm nhiệm Trƣờng Cao Đẳng Nghề Kỹ Thuật Công Nghệ TPHCM i Luan van Giảng viên 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 Luan van 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 cô bạn bè, anh chị để luận văn đƣợc hồ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 Luan van 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 tồ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 Luan van 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 Luan van 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 Luan van 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 Luan van  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 Luan van 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 Luan van [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 Luan van 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 Luan van 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 Luan van 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 Luan van 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 Luan van 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 Luan van 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 q độ 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ố qn 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 Luan van 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 Luan van 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 Luan van 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 Luan van 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 Luan van 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 Luan van Luan van ... 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ó phức tạp Thơng thường phân 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... đá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... lệch cho phép mà hệ thống có ổn định Đó miền giới hạn ổn định hệ thống với kích động lớn Ổn định động hệ thống điện nghiên cứu sở khái niệm Lyapunov Xét lại ví dụ ổn định động cho hệ thống điện