Ứng dụng facts trong vấn đề nâng cao ổn định điện áp

Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến PHẦN I MỞ ĐẦU I ĐẶT VẤN ĐỀ : Hiện nay,vấn đề ổn định điện áp khơng cịn vấn đề lạ tất Tuy nhiên, đóng vai trị quan trọng việc nghiên cứu đánh giá hệ thống điện, gần vấn đề ổn định điện áp nguyên nhân dẫn đến tượng sụp đổ điện áp Vì việc khảo sát, phân tích ổn định điện áp để xác định thông tin giới hạn ổn định chế gây ổn định ổn định điện áp để từ đưa dự báo, tìm biện pháp cải thiện, khắc phục kịp thời tượng sụp đổ điện áp xảy Nếu khơng có dự báo tương đối xác ổn định điện áp để đưa biện pháp cải thiện kịp thời gây hậu nghiêm trọng, gây ảnh hưởng đến phát triển kinh tế an ninh hệ thống điện Việc dự báo sụp đổ điện áp hệ thống điện tốn quan trọng q trình phân tích ổn định điện áp, đặc biệt hệ thống điện lớn ,đường dây dài phức tạp Để phần hiểu rõ vấn đề này, luận văn tìm hiểu phương pháp phân tích ổn định điện áp hệ thống điện với phương pháp phân tích modal kết hợp việc thành lập đường cong V-P, Q-V thành lập từ toán phân bố cơng suất Khi ấy, ta kết luận hệ thống điện ổn định, ổn định hay sụp đổ Từ đó, ta dựa vào giá trị riêng Lamda bé hệ số tham gia để đánh giá xem nút có độ ổn định hay gần với điểm tới hạn ranh giới ổn định nhất, biết nhánh máy phát quan trọng việc tham gia giữ ổn địnhđiện áp để có biện pháp cải thiện kịp thời Trong luận văn dự kiến sử dụng thiết bị bù tĩnh tụ điện thiết bị bù cơng suất phản kháng có điều khiển ứng dụng công nghệ FACTS mà cụ thể thiết bị SVC để bù vào nút khảo sát ổn định hệ thống để cải thiện nâng cao ổn định điện áp hệ thống Vì lí chọn tên đề tài “ ỨNG DỤNG FACTS TRONG VẤN ĐỀ NÂNG CAO ỔN ĐỊNH ĐIỆN ÁP” Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Ở nước ta đà kinh tế phát triển gia nhập vào WTO nhu cầu phụ tải “nóng” cơng suất nguồn phát cịn hạn chế Do vấn đề an ninh hệ thống mối quan tâm lớn ngành lượng nói chung ngành điện nói riêng, việc khảo sát, phân tích đánh giá ổn định điện áp đặt lên mối quan tâm hàng đầu để đưa giải pháp ngắn hạn chiến lược lâu dài để phát triển nguồn điện nhằm đáp ứng nhu cầu phụ tải ngày lớn II MỤC TIÊU VÀ NHIỆM VỤ: Khảo sát ổn định điện áp, phương pháp nghiên cứu đánh giá ổn định điện áp Giới thiệu ứng dụng thiết bị bù công suất phản kháng sử dụng công nghệ FACTS SVC nâng cao ổn định điện áp Áp dụng khảo sát , đánh giá ổn định điện áp cho lưới truyền tải 220kV hệ thống điện miền Nam chưa bù có bù SVC Một số giải pháp nhằm nâng cao ổn định điện áp hệ thống III PHẠM VI NGHIÊN CỨU Đề tài nghiên cứu đánh giá ổn định điện áp cho mạng mẫu IEEE 14 nút mạng truyền tải 220kV hệ thống điện miền Nam Đưa biện pháp cải thiện ổn định ứng dụng FACTS IV GIÁ TRỊ THỰC TIỄN CỦA ĐỀ TÀI Để đảm bảo ổn định cho hệ thống 220kV , nâng cao vấn đề an ninh lượng hệ thống điện Mở rộng cho mạng phức tạp nhiều nút V NỘI DUNG DỰ KIẾN: Phần 1: MỞ ĐẦU Phần 2: KHẢO SÁT ỔN ĐỊNH ĐIỆN ÁP VÀ TÌM HỂU VỀ FACTS ĐẾN ỔN ĐỊNH ĐIỆN ÁP Chương Giới thiệu tổng quan Chương2 Các khái niệm liên quan đến ổn định điện áp Chương3 Facts Chương4 Phương pháp phân tích modal Q-V đánh giá ổn định điện áp hệ thống điện Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Phần 3: KHẢO SÁT, PHÂN TÍCH VÀ ĐÁNH GIÁ ỔN ĐỊNH ĐIỆN ÁP, BÙ NGANG NÂNG CAO ỔN ĐỊNH ĐIỆN ÁP CHO LƯỚI TRUYỀN TẢI 220KV HỆ THỐNG ĐIỆN MIỀN NAM Chương Phân tích, đánh giá ổn định điện áp phương pháp phân tích modal Q-V , bù ngang nâng cao ổn định tụ SVC cho mạng mẫu IEEE 14 nút Chương Áp dụng phân tích đánh giá ổn định điện áp, bù nâng cao ổn điện áp định cho lưới 220Kv hệ thống điện miền Nam Phần 4: KẾT LUẬN Phần 5: PHỤ LỤC DỮ LIỆU MẠNG IEEE 14 NÚT &CHƯƠNG TRÌNH PHÂN TÍCH MODAL&TÀI LIỆU THAM KHẢO Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến PHẦN II KHẢO SÁT ỔN ĐỊNH ĐIỆN ÁP VÀ TÌM HIỂU VỀ FACTS ĐẾN ỔN ĐỊNH ĐIỆN ÁP CHƯƠNG GIỚI THIỆU TỔNG QUAN I KHÁI NIỆM CHUNG VỀ ỔN ĐỊNH HỆ THỐNG ĐIỆN: 1.1 Khái niệm ổn định hệ thống điện: Hệ thống điện tập hợp phần tử phát, dẫn, phân phối có mối quan hệ tương tác lẫn phức tạp Tồn vô số nhiễu tác động lên hệ thống Hệ thống phải đảm bảo tính ổn định có tác động nhiễu Ổn định hệ thống điện khả trở lại vận hành bình thường ổn định sau chịu tác động nhiễu Đây điều kiện thiết yếu để hệ thống tồn vận hành Như chế độ xác lập chẳng hạn, để tồn cần phải có cân cơng suất hệ (khi thơng số hệ giữ khơng đổi) đồng thời phải trì độ lệch nhỏ thơng số định mức kích động ngẫu nhiên nhỏ (làm thông số lệch khỏi giá trị điểm cân bằng) Hoặc tác động thao tác đóng cắt, hệ thống điện cần phải chuyển từ chế độ xác lập sang chế độ xác lập khác Hậu việc ổn định: Khi hệ thống ổn định phải cắt hàng loạt tổ máy, phụ tải, làm tan rã hệ thống gây thiệt hại nghiêm trọng cho kinh tế Do cần nghiên cứu ổn định thiết kế vận hành hệ thống nhằm đảm bảo [5]: - Ổn định tình vận hành bình thường sau cố - Có thể vận hành bình thường tình thao tác vận hành kích động cố 1.2 Ổn định động ổn định tĩnh: Ổn định tĩnh khả hệ thống sau kích động nhỏ (nhiễu nhỏ) phục hồi chế độ ban đầu Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Nếu liên tưởng đến ổn định hai hệ sau ta hiểu thêm khái niệm ổn định tĩnh Ở hệ (H.1.1a) vị trí cân lắc ổn định B A a) b) Hình 1.1: Minh họa ổn định hai hệ học Nghĩa giả sử có nhiễu nhỏ (thường xuyên hữu khơng gian) lắc dao động Tuy nhiên lực cản khơng khí, dao đơng tắt dần lắc trở vị trí ban đầu Tương tự vậy, vị trí A hệ (H.1.1b) ổn định Tuy nhiên vị trí B vị trí giả sử cân khơng ổn định Lý cần nhiễu nhỏ (gió nhẹ) bi rời khỏi điểm B Hệ thống điện, nói, chế độ xác lập có điều kiện cân cơng suất có thông số không thay đổi biến thiên nhỏ xung quanh giá trị ban đầu Tuy nhiên, hệ thống vận hành có nhiều tác động nhỏ, ngẫu nhiên, ví dụ thay đổi công suất phụ tải Hệ thống phải trì độ lệch nhỏ trở trị ban đầu thơng số chế độ Tính chất tính ổn định tĩnh hệ Ổn định động hệ thống khả hệ phục hồi trạng thái ban đầu gần với trạng thái ban đầu sau kích động lớn (nhiễu lớn) Như nói chuyển sang chế độ xác lập hệ có tính ổn định động Các kích động lớn hiểu như: - Ngắn mạch phần tử lưới điện - Đóng cắt phần tử lưới điện - Tăng giảm tải đột ngột Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Giả sử có máy phát điện phát công suất P0 nhu cầu cần phải tăng công suất lên P1 Ở thời điểm có tăng đột ngột cơng suất P = P1 – P0 nên máy phát quay nhanh lên Nếu tăng tải đảm bảo cho máy phát theo thời gian hãm tốc lại trở vị trí cân ổn định ta nói hệ có tính ổn định động Ngược lại, máy phát liên tục tăng tốc, rời khỏi đồng bắt buộc phải cắt khỏi lưới, hệ ổn định động 1.3 Mô hình hệ thống điện – Khơng gian trạng thái Mơ hình hệ thống điện – Khơng gian trạng thái, trạng thái hệ thống động hệ thống điện mơ tả tập hợp N phương trình vi phân bậc dạng phi tuyến sau : X i  f i ( x1 , xn ; u1 , u r ; t ) Trong : N :là số biến hệ thống R :số biến điều khiển Hệ phương trình viết lại dạng ma trận sau: X  f i ( x , u , t ) Trong đó: X= (x1,………,xn)T vectơ trạng thái hệ thống U= (u 1,………,u r)T vectơ đầu vào (biến điều khiển )của hệ thống f= (f1,………,fn)T Trong trường hợp đạo hàm biến trạng thái khơng hàm explicit hệ thống gọi hệ thống autonomus Khi phương trình viết đơn giản lại X  f i ( x , u ) Ngòai , cần phải quan tâm biến đầu có mối quan hệ đại số với biến trạng thái biến điều khiển diễn tả dạng : Y=g(x,u) Trong : y= (y1,………,ym)T vec tơ biến đầu (output ) hệ thống g= (g1,………,gm)T Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Đối với hệ thống máy điện đồng kích từ dạng IEEE lọai , phương trình mơ tả hệ thống viết lại sau ; Các phương trình vi phân kích từ điều tốc : d i  i   s dt , , , , d i TMi [ E qi  X di I di ]I qi [ E di  X qi I qi ]I di Di ( i   s )     dt Mi Mi Mi Mi dEqi' dt  E qi' ' Tdoi  [ X di  X di, ]I di ] ' Tdoi  E ,fdi ' Tdoi I qi dE di' E di'   '  ' ( X qi  X qi' ) dt Tqoi Tqoi dE 'fdi dt  K Ei  S E E fdi TEi E fdi  V Ri TEi dVRi V K K K K   Ri  Ai R fi  Ai Fi E fdi  Ai (Vrefi  Vi ) dt T Ai TAi T AiTFi T Ai dRFi R K Fi   Fi  E Fdi dt TFi (TFi ) Các phương trình đại số stator: E di'  Vi sin( i   i )  Rsi I di  X qi' I qi  E qi'  Vi cos( i   i )  R si I qi  X di' I di  Các phương trình hệ thống: n I diVi sin( i   i )  I qiVi cos( i   i )  PLi (Vi )   ViVk Yik cos( i   k   ik )  k 1 n I diVi cos( i   i )  I qiVi sin( i   i )  Q Li (Vi )   ViVk Yik sin( i   k   ik )  k 1 n PLi (Vi )   ViVk Yik cos( i   k   ik )  k 1 n QLi (Vi )   ViVk Yik sin( i   k   ik )  k 1 1.4 Tuyến tính hóa mơ hình hệ thống Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Cho x0 vectơ trạng thái cân , u0 vectơ đầu vào tương ứng với trạng thái cân với dao động nhỏ Do x0, u0 biến trạng thái cân nên thỏa mãn phương tình sau : x  f ( x , u )  Cho dao động nhỏ trạng thái X = x0 +  x u = u0 +  x Trạng thái hệ thống thoải mãn với phương trình :   x  x  x i  f i  x0  x , ( x0   u )     f i ( x0 , u )  Vì x i  f i f f f xi    i x n  i ui   i u r xi x ni ui u r x i  f i ( x0 , u ) , nên : f i f f f xi    i xn  i u i   i u r xi xni ui u r với i=1,2,3 n Với cách tương tự ta có : y i  g i g g g xi    i xm  i u i   i u r xi xm ui u r Với i=1,2, m x  Ax  Bu y  Cx  Du Trong : f1   f i  x x  n    A =     f n f n   x1 x n  f1   f i  u u  r   B=      f n f n   u1 u n  Học viên: Huỳnh Vũ Quốc Khánh  g i g1   x x  n   C=      g n g m   x1 x n  CBHD: TS Hồ Văn Hiến  g i g1   u u  r   D=      g n g n   u1 u r  Trong : x : vectơ trạng thái bậc n y vectơ đầu bậc m u :vectơ đầu vào( vectơ điều khiển ) bậc r A =ma trận trạng thái (ma trận Jacobi )bậc nxn B= ma trận điều khiển bậc nxr C= ma trận đầu bậc mxn D= ma trận hồi tiếp, tỉ lệ tín hiệu đầu vào ảnh hưởng trực tiếp đến tín hiệu đầu với ma trận A ma trận Jacobi mơ hình động hệ thống, định tính chất động học hệ thống 1.5 Điểm cân (Equibibrium or singular points) Điểm cân điểm mà tất đạo hàm x1 x n đồng thời không, hệ thống xem dừng tất biến số, không thay đổi theo thời gian Điểm cân thỏa mãn điều kiện sau : Fi(x0,u)=0 Nếu phương trình fi(i=1,2, n)là tuyến tính, hệ thống gọi hệ tuyến tính có điểm cân Ngươc lại, hệ thống phi tuyến thơng thường khảo sát có nhiều điểm cân đặc trưng khảo sát hệ thống động đó, đồng thời rút tính chất ổn định hệ thống Sau khảo sát điểm cân biết trạng thái ổn định như: + Ổn định cục +Ổn định giới hạn +Ổn định tòan cục Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến I KHÁI NIỆM VỀ ỔN ĐỊNH ĐIỆN ÁP: 2.1 Khái niệm ổn định điện áp Theo [4] ổn định điện áp hiểu khả hệ thống điện giữ điện áp chấp nhận tất nút hệ thống điện vận hành chế độ bình thường vận hành chế độ bất bình thường Ổn định điện áp thường thường chia thành hai loại : - Ổn định nhiễu nhỏ :là khả điều khiển điện áp hệ thống xuất tượng nhiễu nhỏ : ví dụ gia tăng dần phụ tải - Ổn định nhiễu lớn :là khả diều khiển điện áp hệ thống xuất tượng nhiễu lớn Ví dụ ngắn mạch, phụ tải đột ngột ( có giá trị cơng suất đáng kể ), máy phát Việc ổn định điện áp chia thành nhiều cấp độ khác theo thứ tự thời gian - Mất ổn định độ : thời gian xảy khỏang 10 giây - Mất ổn định lâu dài :thời gian tồn đến 10 phút 2.2 Các nguyên nhân đặc trưng tượng sụp đổ điện áp : Như biết, có nhiều nguyên nhân gây tượng ổn định điện áp, tượng dẫn đến tượng sụp đổ điện áp Biểu đặc trưng sụp đổ điện áp hệ thống có gia tăng đột ngột cơng suất kháng, nguồn công kháng thêm vào phải bổ sung máy phát thiết bị bù Nếu độ dự trữ cơng suất phản kháng đủ hệ thống trở lại ổn định, ngược lại công suất phản kháng khơng đủ đáp ứng hệ thống ổn định dẫn đến tượng sụp đổ điện áp Khi tượng sụp đổ điện áp xảy ảnh hưởng hưởng đến trình cung cấp điện vùng hay toàn vùng Hậu việc ổn định điện áp phải cắt hàng loạt tổ máy, phụ tải, nghiêm trọng gây tượng rã lưới, dẫn đến thiệt hại kinh tế đáng kể [6] Vì cần phải nghiên cứu vấn đề cach cụ thể để vận dụng vào thiết kế, vận hành hệ thống cho đảm bảo tình trạng ổn định hệ thống tình vận hành bình thường sau cố xảy Sụp đổ điện áp thường xảy khi: 10 Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến if kb(n) == P(n)=Pk; Q(n) = Qk; end % Swing bus P if kb(n) == Q(n)=Qk; if Qmax(n) ~= Qgc = Q(n)*basemva + Qd(n) - Qsh(n); if iter % the Mvar of generator buses are if Qgc < Qmin(n), % tested If not within limits Vm(n) Vm(n) = Vm(n) + 0.01; % is changed in steps of 0.01 pu to elseif Qgc > Qmax(n), % bring the generator Mvar within Vm(n) = Vm(n) - 0.01;end % the specified limits else, end else,end else,end end if kb(n) ~= A(nn,nn) = J11; %diagonal elements of J1 DC(nn) = P(n)-Pk; end if kb(n) == A(nn,lm) = 2*Vm(n)*Ym(n,n)*cos(t(n,n))+J22; %diagonal elements of J2 A(lm,nn)= J33; %diagonal elements of J3 A(lm,lm) =-2*Vm(n)*Ym(n,n)*sin(t(n,n))-J44; %diagonal of elements of J4 DC(lm) = Q(n)-Qk; end end % ket thuc vong lap for dem=1:nbus DX=A\DC'; for n=1:nbus nn=n-nss(n); lm=nbus+n-ngs(n)-nss(n)-ns; if kb(n) ~= delta(n) = delta(n)+DX(nn); end if kb(n) == Vm(n)=Vm(n)+DX(lm); end end maxerror=max(abs(DC)); if iter == maxiter & maxerror > accuracy fprintf('\nWARNING: Iterative solution did not converged after ') fprintf('%g', iter), fprintf(' iterations.\n\n') fprintf('Press Enter to terminate the iterations and print the results \n') converge = 0; pause, else, end end % ket thuc vong lap while if converge ~= tech= (' ITERATIVE SOLUTION DID NOT CONVERGE'); else, tech=(' Power Flow Solution by Newton-Raphson Method'); end V = Vm.*cos(delta)+j*Vm.*sin(delta); deltad=180/pi*delta; i=sqrt(-1); k=0; for n = 1:nbus if kb(n) == k=k+1; S(n)= P(n)+j*Q(n); Pg(n) = P(n)*basemva + Pd(n); Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n); Pgg(k)=Pg(n); Qgg(k)=Qg(n); %june 97 elseif kb(n) ==2 k=k+1; S(n)=P(n)+j*Q(n); 144 Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n); Pgg(k)=Pg(n); Qgg(k)=Qg(n); % June 1997 end yload(n) = (Pd(n)- j*Qd(n)+j*Qsh(n))/(basemva*Vm(n)^2); end busdata(:,3)=Vm'; busdata(:,4)=deltad'; Pgt = sum(Pg); Qgt = sum(Qg); Pdt = sum(Pd); Qdt = sum(Qd); Qsht = sum(Qsh); End 2.3.5 Tập tin jacobireduce.m tính matrận jacobi A1 ma trận -1 Jpd,Jpv,Jqd,Jqv , sau tính matrân Jacobi rút gọn Jr %tinh ma tran Jacobi A1 tinh den cac nut phu tai (P,Q) va cac nut dieu %khien dien ap (P,V), c1=2*nbus-2*ns; m1=nbus-ns; A1=zeros(c1,c2); for n=1:nbus nn=n-nss(n); lm=nbus+n-nss(n)-ns; J11=0; J22=0; J33=0; J44=0; for chay=1:nbr if nl(chay) == n | nr(chay) == n if nl(chay) == n, l = nr(chay); end if nr(chay) == n, l = nl(chay); end J11=J11+ Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l)); J33=J33+ Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n) + delta(l)); if kb(n)~=1 J22=J22+ Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n) + delta(l)); J44=J44+ Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l)); else, end if kb(n) ~= & kb(l) ~=1 lk = nbus+l-nss(l)-ns; ll = l -nss(l); % off diagonalelements of J1 A1(nn, ll) =-Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l)); % off diagonal elements of J2 A1(nn, lk) =Vm(n)*Ym(n,l)*cos(t(n,l)- delta(n) + delta(l)); % off diagonal elements of J3 A1(lm, ll) =-Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n)+delta(l)); % off diagonal elements of J4 A1(lm, lk) =-Vm(n)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l)); else end else , end end if kb(n) ~= A1(nn,nn) = J11; %diagonal elements of J1 A1(nn,lm) = 2*Vm(n)*Ym(n,n)*cos(t(n,n))+J22; %diagonal elements of J2 A1(lm,nn)= J33; %diagonal elements of J3 A1(lm,lm) =-2*Vm(n)*Ym(n,n)*sin(t(n,n))-J44; %diagonal of elements of J4 end end %ket thuc vong lap for i=1:nbus % tinh cac ma tran Jacobi rut gon Jpd,Jpv,Jqd,Jqv Jpd=zeros(m1,m1); Jpv=zeros(m1,m1); Jqd=zeros(m1,m1); Jqv=zeros(m1,m1); for couter1=1:m1 for couter2=1:m1 Jpd(couter1,couter2)=A1(couter1,couter2); 145 Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến end, end for couter1=1:m1 for couter2=1:m1 Jpv(couter1,couter2)=A1(couter1,m1+couter2); end, end for couter1=1:m1 for couter2=1:m1 Jqd(couter1,couter2)=A1(m1+couter1,couter2); end, end for couter1=1:m1 for couter2=1:m1 Jqv(couter1,couter2)=A1(m1+couter1,m1+couter2); end, end % tinh ma tran Jacobi rut gon Jr va Jr^-1 Jr=Jqv-(Jqd*(Jpd^-1)*Jpv); Jr1=Jr^-1; clear A A1 DC DX J11 J22 J33 J44 Qk delta lk ll lm 2.3.6 Tập tin in kết điện áp nút busout.m disp(tech) fprintf(' Sai so cong suat = %g \n', maxerror) fprintf(' So lan lap = %g \n\n', iter) head =[' Bus Voltage Angle Load -Generation - Injected' ' No Mag Degree MW Mvar MW Mvar Mvar ' ' ']; disp(head) for n=1:nbus fprintf(' %5g', n), fprintf(' %7.3f', Vm(n)), fprintf(' %8.3f', deltad(n)), fprintf(' %9.3f', Pd(n)), fprintf(' %9.3f', Qd(n)), fprintf(' %9.3f', Pg(n)), fprintf(' %9.3f ', Qg(n)), fprintf(' %8.3f\n', Qsh(n)) end fprintf(' \n'), fprintf(' Total ') fprintf(' %9.3f', Pdt), fprintf(' %9.3f', Qdt), fprintf(' %9.3f', Pgt), fprintf(' %9.3f', Qgt), fprintf(' %9.3f\n\n', Qsht) 2.3.7 Tập tin Modal.m phân tích ổn định điện áp từ ma trận Jr-1 %Chuong trinh tinh ma tran Modal va cac he so tham gia Pki, Pji,Pmi % tinh ma tran Modal [Ri,E]=eig(Jr1); % E cac vector tri rieng cua ma tran Jacobi Jr1 Le=Ri^-1; % Ri: vecto rieng phai; Le:vecto rieng trai lamda=diag(E); % cac gia tri Lamda sap theo cot for k=1:length(lamda) lamda(k)=lamda(k)^-1;%Cac gia tri Lamda nghich dao end; min_i=find(lamda==min(lamda));% vi tri gia tri lamda min(sap)theo cot minlamda=min(lamda);% gia tri lamdamin for c1=1:(length(Jr1)-1)%sap xep gia tri lamda theo thu tu tang dan for c2=1:(length(Jr1)-c1) if lamda(c2)>lamda(c2+1) tam=lamda(c2); lamda(c2)=lamda(c2+1); 146 Học viên: Huỳnh Vũ Quốc Khánh CBHD: TS Hồ Văn Hiến lamda(c2+1)=tam; end end end % tinh cac he so tham gia cua nut pki ung voi tri so lamda nho nhat kk=0 for n=1:nbus kk1=n-nss(n); if kk1~=kk kk=kk1; Bus_Pki(kk,1)=n; else,end; end for k=1:length(Jr1) Bus_Pki(k,2)=Ri(k,min_i)*Le(min_i,k); end for c1=1:(length(Jr1)-1) for c2=1:(length(Jr1)-c1) if Bus_Pki(c2,2)