(Luận văn thạc sĩ) ổ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 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 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 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 Giảng viên LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu tơi 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 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 Chƣơng TỔNG QUAN 1.1 Tổng quan lĩnh vực nghiên cứu Lịch sử phát minh sử dụng điện đánh dấu kiện đáng ghi nhớ  Năm 1883, phát minh hệ thống điện xoay chiều pha  Năm 1884, Pháp lần điện tải xa dòng điện xoay chiều  Năm 1885, đánh dấu việc phát minh máy biến áp điện lực  Năm 1891, lần đường dây tải điện pha vận hành thử nghiệm khoảng cách 175km Từ lịch sử phát triển điện gắn liền với lịch sử hình thành phát triển hệ thống điện xoay chiều pha công suất lớn: khoảng cách truyền tải điện ngày tăng, công suất đơn vị tổ máy phát điện đồng pha ngày lớn, số tổ máy phát làm việc song song hệ thống ngày nhiều Chính đặc điểm cuối dẫn đến yêu cầu phát triển lý thuyết ổn định hệ thống điện Năm 1928, nhà bác học người Mỹ R Park lần đặt móng nghiên cứu ổn định hệ thống điện dựa sở thiết lập hệ phương trình vi phân trình độ điện máy điện đồng hệ tọa độ quay Gần đồng thời với Park, loạt công trình cơng bố độc lập nhà bác học người Nga A.A Goriev năm 1930 – 1935 mơ hình q trình q độ máy điện quay làm phát triển lý thuyết nghiên cứu ổn định hệ thống điện thêm bước, sau hệ phương trình gọi tên Park – Goriev Mơ hình q trình q độ hệ thống điện hệ tọa độ quay làm đơn giản đáng kể hệ phương trình vi phân mơ tả trạng thái độ hệ thống điện Dựa sở phương pháp tốn ổn định hệ thống nói chung nghiên cứu áp dụng cho hệ thống điện Trước hết phải kể đến phương pháp dựa khái niệm cân lượng, gọi tiêu chuẩn lượng Phương pháp đơn giản, nhận kết dễ áp dụng nhiều trường hợp Tuy nhiên, phương pháp khác sau rằng, khái niệm ổn định theo ý nghĩa cân lượng không đầy đủ, không phát hiện tượng ổn định dao động quán tính Hơn nữa, phương pháp cân lượng khơng có sở chặt chẽ phương pháp để áp dụng hệ thống điện phức tạp Tiếp đến phương pháp toán nghiên cứu ổn định hệ vật lý nói chung A M Lyapunov phát triển áp dụng cho hệ thống điện, đặc biệt phương pháp dao động bé Chính Goriev cơng trình cách nghiên cứu ổn định hệ thống điện theo Lyapunov, đặc biệt hệ thống điện đơn giản, phức tạp Ông chứng minh sử dụng tiêu chuẩn dấu dương số hạng tự phương trình đặc trưng hệ phương trình vi phân trình q độ phát hầu hết trường hợp ổn định hệ thống điện, gọi tiêu chuẩn ổn định phi chu kỳ Dựa tiêu chuẩn cách tính tốn phân tích ổn định cho hệ thống điện phức tạp dạng chung xây dựng áp dụng phổ biến chương trình phân tích hệ thống điện Tiêu chuẩn đặc biệt tiện lợi tìm thơng số giới hạn chế độ theo điều kiện ổn định Dựa sở tiêu chuẩn ổn định phi chu kỳ I M Markovits chứng minh chất tiêu chuẩn lượng, chúng trường hợp riêng đảm bảo ổn định phi chu kỳ Nhờ tiêu chuẩn riêng đánh giá nhanh mức độ ổn định hệ thống điện cụ thể vận hành Vấn đề chỗ, điều kiện cụ thể tồn hay vài tiêu chuẩn (dạng đạo hàm dP dQ dE ' ' ) dễ bị vi phạm Trong vận hành d dU dU  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 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ố 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