ỨNG DỤNG MATLAB LẬP TRÌNH TÍNH TOÁN PHÂN BỐ CÔNG SUẤT CHO HỆ THỐNG TRUYỀN TẢI CAO ÁP ĐỒNG BẰNG SÔNG CỬU LONG

TRƯỜNG ĐẠI HỌC CẦN THƠ CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAMKHOA CÔNG NGHỆ ặ Độc lập Tự do Hạnh phúcBộ MỒN KỸ THUẬT ĐIỆNNHẬN XÉT VÀ ĐÁNH GIÁ CỦA CÁN BỘ HƯỚNG DẨN1. Cán bộ hướng dẫn:2. Đề tài: ứng dụng MATLAB lập trình tính toán phân bố công suất chohệ thống truyền tải cao áp Đồng Bằng Sông Cửu Long.3. Sinh viên thực hiện: Phạm Ngọc Minh MSSV: 10911294. Lớp: Kỹ Thuật ĐiệnTRƯỜNG ĐẠI HỌC CẦN THƠ CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAMKHOA CÔNG NGHỆ ặ Độc lập Tự do Hạnh phúcBộ MỒN KỸ THUẬT ĐIỆNNHẬN XÉT VÀ ĐÁNH GIÁ CỦA CÁN BỘ HƯỚNG DẨN1. Cán bộ hướng dẫn:2. Đề tài: ứng dụng MATLAB lập trình tính toán phân bố công suất chohệ thống truyền tải cao áp Đồng Bằng Sông Cửu Long.3. Sinh viên thực hiện: Phạm Ngọc Minh MSSV: 10911294. Lớp: Kỹ Thuật Điện

B ộ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC CẦN THƠ KHOA CƠNG NGHỆ•_ _ LUẬN VĂN TỐT NGHIỆP ĐẠI HỌC • • • • ỨNG DỤNG MATLAB LẬP TRÌNH TÍNH TỎÁN PHÂN BỐ CồNG SUẤT CHO HỆ THỐNG TRUYỀN TẢI CAO ÁP ĐỒNG BẰNG SÔNG c u LONG CÁN B ộ HƯỚNG DẪN: SINH VIÊN THựC HIỆN: K Trần Anh Nguyện Phạm Ngọc Minh (MSSV: 1091129) Ngành: Kỹ thuật điện - Khóa: 35 Tháng 05/2013 TRƯỜNG ĐẠI HỌC CẦN THƠ KHOA CÔNG NGHỆ B ộ MỒN KỸ THUẬT ĐĨỆN CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - T ự - H ạnh phức Cần Thơ, ngày 09 tháng 01 năm 2013 PHIẾU ĐỀ NGHỊ ĐỀ TÀI TỐT NGHIỆP CỦA SINH VIÊN HỌC KỲ II NĂM HỌC 2012 - 2013 Họ tên sinh viên: Ngành: Phạm Ngọc M inh Kỹ T huật Điện MSSV: 1091129 Khoá: 35 Tên đề tài: ứng dụng MATLAB lập trình tính tốn phân bố công suất cho hệ công cho hệ thống truyền tải cao áp Đồng Bằng Sông Cửu Long Đia điểm thưc hiên: Bô môn Kỹ T huât Điên - Khoa Công Nghê - T rường Đai Học Cần Thơ Họ tên người hướng dẫn khoa học (NHDKH): Kỹ sư.Trần Anh Nguyện Mục tiêu đề tài: Viết chương trình Matlab tính tốn phân bố cơng suất tần thất hệ thống truyền tải cao áp Đồng Bằng Sông Cửu Long dựa phương pháp lặp Fast Decoupled, đồng thời hướng đến khả áp dụng vào hệ thống điện khác thực tế Các nội dung giới hạn đề tài: A MỞ ĐÀU B NỘI DUNG CHƯƠNG 1: TỔNG QUAN MATLAB CHƯƠNG 2: CÁC PHƯƠNG PHÁP TÍNH PHÂN BỐ CƠNG SUẤT CHƯƠNG 3: VIẾT CHƯƠNG TRÌNH MATLAB TÍNH PHÂN BỐ CƠNG SUẤT CHƯƠNG 4: ỨNG DỤNG CHƯƠNG TRÌNH ĐÊ TÍNH PHÂN BỐ CÔNG SUẤT CHO HỆ THỐNG TRUYỀN TẢI CAO ÁP THựC TẾ CHƯƠNG 5: KIỂM TRA LẠI CHƯƠNG TRÌNH BẰNG PHẦN MỀM POWERWORLD c K Ế T LUẬN Các yêu cầu hỗ trợ cho việc thực đề tài: Do kiến thức em hạn chế cần có sử hổ trợ cán hướng dẫn thầy Bộ môn để em hồn thành đề tài SINH VIÊN ĐỀ NGHỊ (Ký tên ghi rõ họ tên) Ý KIẾN CỦA NHDKH Ý KIẾN CỦA B ộ MÔN Ý KIẾN CỦA HỘI ĐỒNG LV&TLTN TRƯỜNG ĐẠI HỌC CẦN THƠ KHOA CÔNG NGHỆ B ộ MỒN KỸ THUẬT ĐIỆN CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM ặ Độc lập - Tự - Hạnh phúc NHẬN XÉT VÀ ĐÁNH GIÁ CỦA CÁN BỘ HƯỚNG DẨN Cán hướng dẫn: Đề tài: ứng dụng MATLAB lập trình tính tốn phân bố cơng suất cho hệ thống truyền tải cao áp Đồng Bằng Sông Cửu Long Sinh viên thực hiện: Phạm Ngọc Minh MSSV: 1091129 Lớp: Kỹ Thuật Điện Khoá: 35 Nội dung nhận xét: a Nhận xét hình thức LVTN: b Nhận xét vẽ: c Nhận xét nội dung LVTN: ❖ Đánh giá nội dung thực đề tài: ❖ Những vấn đề hạn chế: d Kết luận đề nghị: Điểm đánh giá: Cần Thơ, ngày tháng ểểnăm 2013 Cán hướng dẫn TRƯỜNG ĐẠI HỌC CẦN THƠ KHOA CÔNG NGHỆ Bộ MỒN KỸ THUẬT ĐIỆN CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc NH Ậ N XÉT V À Đ Á N H G IÁ CỦA CÁN B ộ PH Ả N BIỆN • • • Cán phản biện: Tên đề tài: Khảo sát tính ổn định mạng điện Quận Cái Răng phần mềm Power-Word Sinh viên thực hiện: Lê Việt Tuấn (MSSV: 1091225) Lớp: Kỹ Thuật Điện - Khóa 35 Nội dung nhận xét: a Nhận xét hình thức tập thuyết minh: bệ Nhận xét vẽ: cệ Nhận xét nội dung tiểu luận: Các công việc đạt được: Những vấn đề hạn chế: d Nhận xét sinh viên thực đề tài: e Kết luận đề nghị: Điểm đánh giá: Cần Thơ, ngày tháng năm 2013 CÁN Bộ PHẢN BIỆN LỜI NÓI ĐÀU Cùng với đổi phát triển đất nước giai đoạn nay, kinh tế nước nhà tăng trưởng mạnh Đảng phủ quan tâm đến mặt xã hội, vấn đề truyền tải điện cững nhận quan tâm đặc biệt nhà nước tổ chức, thu hút nhiều dự án đầu tư cơng trình phát triển nhằm giải cách tốt vấn đề độ tin cậy hệ thống điện nói chung vấn đề cung cấp đầy đủ điện cho nhân dân chất lẫn lượng nói riêng Do đó, vấn đề cấp thiết phải vận hành hệ thống điện tốt điều kiện cân pha bình thường điều kiện ổn định, muốn đạt hệ thống điện phải đảm bảo yêu cầu sau: • Nguồn điện (generation) cung cấp thỏa mãn yêu cầu công suất phụ tải lượng công suất tổn thất hệ thống • Độ lớn điện áp nút phải trì gần với giá trị quy định • Vận hành máy phát (generator) cung cấp cơng suất tác dụng công suất phản kháng xác định tới hệ thống nằm giới hạn cho phép • Tất đường dây truyền tải máy biến áp điều khơng trạng thái q tảiễ Bài tốn khảo sát phân bố công suất áp dụng cho hệ thống điện ba pha cân dựa sơ đồ tương đương pha Thông thường sử dụng hệ đơn vị tương đối Khảo sát phân bố cơng suất địi hỏi số liệu hệ thống chi tiết giá trị góc lệch pha Giải tốn phân bố công suất cho giá trị điện áp (biên độ diện áp) góc pha (góc điện áp) điểm nút, dịng cơng suất nhánh tổn thất cơng suất mạng điện Bài tốn phục vụ cho công tác thiết kế vận hành hệ thống, điều chỉnh điện áp công suất, tối ưu hóa kinh tế vận hành hệ thống, v.v Sinh Viên TH: Phạm Ngọc Minh i LỜI CẢM ƠN cso&CŨIsosd Trong trình học tập luận văn tốt nghiệp tổng hợp kiến thức vốn có sinh viên q trình học tập, hội để sinh viên kiểm chứng lại kiến thức trước bước vào môi trường làm việc Dưới hướng dẫn quý thầy, mơn sinh viên có thêm kiến thức kỹ sinh viên Sau 14 tuần thực hiện, đến luận văn em hoàn thành Nhân em xin gửi lời cảm ơn đến thầy Trần Anh Nguyện, lời cảm ơn chân thành lòng biết ơn sâu sắc hướng dẫn chu đáo, nhiệt tình thầy cảm ơn quý thầy, cô dạy cho em kiến thức vô quý báu Em xin gửi lời cảm ơn đến tất quý thầy Bộ môn Kỹ Thuật Điện Khoa Công Nghệ - Trường Đại Học cần Thơ tạo điều kiện để em hoàn thành luận văn Nhân em xin gửi lời cảm ơn đến gia đình, bạn bè giúp đỡ động viên em suốt trình học Tuy nhiên, kiến thức cịn hạn hẹp, khơng có kinh nghiệm thực tế thời gian làm luận văn có hạn nên khơng thể tránh khỏi thiếu sót Rất mong thơng cảm nhận đóng góp ý kiến quý báu quý thầy, cô để luận văn em hoàn thiện Sinh Viên Thực Hiện Phạm Ngọc Minh Sinh Viên TH: Phạm Ngọc Minh ii MỤC LỤC CBHD: Trần Anh Nguyện MỤC LỤC LỜI NÓI Đ Ầ U i LỜI CÁM Ơ N ii MỤC LỤC iii MỤC LỤC HÌNH VÀ HÌNH VẼ vii MỤC LỤC BẢNG BIÊU ix A MỞ ĐẦU GIỚI THIỆU CHUNG .1 MỤC ĐÍCH NGHIÊN c ứ u VÀ PHẠM VI CỦA ĐỀ TÀI KẾT CẤU CỦA ĐỀ TÀI B NỘI DUNG CHƯƠNG TỔNG QUAN MATLAB 1.1 GIAO ĐIỆN CỦA CHƯƠNG TRÌNH 1.2ỄCÁC PHÉP TOÁN - TÊN BIẾN - CÁC HÀM c BẢN 1.2.1 Các phép toán .5 1.2.2 Cách đặt tên biến 1.2.3 Điều khiển vào r a 1.2.4 Một số hàm toán học 1.3 SỐ PHỨC TRONG MATLAB 1.3.1 Nhập số phức 1.3.2 Các phép toán với số phức 1.4 MA TRẬN VÀ ỨNG DỤNG 1.4.1 Ma trận 1.4.2 Các phép toán với ma trận Matlab 1.4.3 ứng dụng ma trận vào giải hệ phương trình .8 1.5 CẤU TRUC ĐIỀU KIỆN 10 1.5.1 Cấu trúc if-end 10 1.5.2 Cấu trúc if-elseif-else-end .11 1.6 CẤU TRÚC LẶP 12 1.6.1 Cấu trúc for-end 12 1.6.2 Cấu trúc while-end 12 Sinh Viên TH: Phạm Ngọc Minh iii MỤC LỤC CBHD: Trần Anh Nguyện 1.7 GIAO DIỆN ĐỒ HỌA GUIDE TRONG MATLAB .13 1.8 LIÊN KÉT GIỮA MATLAB VA MICROSOFT EXCEL 15 1.8.1 Các tiện ích liên kết 15 1.8.2 Chuyển liệu từ Excel sang Matlab 16 1.8.3 Chuyển liệu từ Matlab sang Excel 17 CHƯƠNG CÁC PHƯƠNG PHÁP TÍNH PHÂN BỐ CƠNG SUẤT 2.1 GIỚI THIỆU 19 2.2 MƠ HÌNH NÚT LƯỚI ĐIỆN 19 2.3 CÁC PHƯƠNG PHÁP GIẢI TÍCH LƯỚI ĐIỆN 24 2.3.1 Phương pháp lặp GAUSS-SEIDEL (GS) 24 2.3.1.1 Giải phương trình phương pháp lặp G S 25 2.3.1.2 Giải phương hệ trình phương pháp lặp GS 26 2.3.1.3 Áp dụng GS vào giải tích lưới điện 28 2.3.2 Phương pháp lặp NEWTON-RAPSHON (NR) .31 2.3.2.1 Sơ lược phương pháp lặp N R 31 2.3.2.2 Giải phương trình phương pháp lặp NR .31 2.3.2.3 Giải phương hệ trình phương pháp lặp N R 32 2.3.2.4 Áp dụng NR vào giải tích lưới điện 34 2.3.3 Phương pháp lặp FAST-DECOUPLED (FD) 38 2.3.3.1 Sơ lược phương pháp lặp FD 38 2.3.3.2 Áp dụng FD vào giải tích lưới điện 41 2.3.3.3 Mơ hình hóa phương pháp FD giải tích lưới điện 45 2.3.3.4 Ví dụ áp dụng FD vào giải tích lưới điện 46 CHƯƠNG VIẾT CHƯƠNG TRÌNH MATLAB TÍNH PHÂN BỐ CƠNG SUẤT 3.1 KHÁI QUÁT VỀ CHƯƠNG TRÌNH 49 3.1.1 Giới thiệu khái quát chương trình viết theo phương pháp F D 49 3.1.2 Các thông số đầu vào 50 3.2 CÁC CHƯƠNG TRÌNH CON TRONG CHƯƠNG TRÌNH TÍNH TỐN PHÂN BỐ CƠNG SUẤT VIẾT THEO PHƯƠNG PHÁP FD .52 3.2.1 Chương trình đổi sang đơn vị tương đối Busdata 52 3.2.2 Chương trình đổi sang đơn vị tương đối Linedata 53 3.2.3 Chương trình đổi sang đơn vị tương đối Tranformerdata 54 3.2.4 Chương trình tạo ma trận Ybus 55 Sinh Viên TH: Phạm Ngọc Minh iv MỤC LỤC CBHD: Trần Anh Nguyện 3.2.5 Chương trình tính tốn phân bố công suất hệ thống truyền tải cao áp Đồng Bằng Sông Cửu Long phương pháp Fast Decouple 57 CHƯƠNG ỨNG DỤNG CHƯƠNG TRÌNH ĐẺ TÍNH PHÂN BỐ CƠNG SUẤT CHO HỆ THỐNG TRUN TẢI CAO ÁP THựC TẾ 4.1 TỔNG QUAN VỀ HỆ THỐNG TRUYỀN TẢI CAO ÁP ĐỒNG BẰNG SÔNG CỬU LONG 'Ế 60 4.2 NHẬN XÉT KÉT QUẢ TÍNH TỐN 61 4.3 NHẬN XÉT PHƯƠNG PHÁP TÍNH TỐN 62 CHƯƠNG KIỂM TRA LẠI CHƯƠNG TRÌNH BẰNG PHẦN MÈM POWERWORLD 5.1 NHẬP CÁC THÔNG SỐ TỪ HỆ THỐNG TRUYỀN TẢI THựC TẾ VÀ XUẤT RA KẾT QUẢ Ể 64 5.1.1 Một số lưu ý q trình nhập thơng số vào PowerWorld 64 5.1.2 Kết mô sau chạy PowerWorld 67 5.1.3 Ma trận Ybus PowerWorld 68 5.1.4 Kết Bus 69 5.1.5 Kết đường dây MBA 70 5.1.6 Tổng hợp thông số máy phát, phụ tải, tổn thất toàn hệ thống „ằẾẾằ’ằẾẾ , Ếằ! Ếằ"ẾẾằằ' Ế 70 5.2 SO SÁNH KẾT QUẢ MÔ PHỎNG POWERWORLD VỚI KÉT QUẢ TỪ CHƯƠNG TRÌNH MATLAB 71 5.2.1 Một số thông số ma trận Ybus 71 5.2.2 Các thông số 10 Bus đầu 72 5.2.3 Thông số công suất truyền đường dây qua MBA 73 5.2.4 Công suất phát, phụ tải tổn thất công suất trông hệ thống 74 C.KẾT LUẬN TÀI LIỆU THAM KHẢO 77 PHU LUC • • Chương trình đổi sang đom vị tương đối Busdata 78 Chương trình đổi sang đơn vị tương đối Linedata 78 Sinh Viên TH: Phạm Ngọc Minh V CBHD: Trần Anh Nguyện TẠI LIỆU THAM KHẢO TÀI LIỆU THAM KHẢO ôer&r&s' Tiếng Việt [1] Võ Ngọc Điều, Bài giảng mơn học Giải Tích Hệ Thống Điện Nâng Cao, Trường Đại Học Bách Khoa [2] Lê Kim Hùng, Giải Tích Mạng, Trường Đại Học Bách Khoa Đà Nắng [3] TS Trần Trung Tính (9/2009), Bài giảng mơn học Hệ Thống Điện 1, Trường Đại Học Cần Thơ [4] TS Trần Trung Tính (2010), Bài giảng mơn học Ồn Định Hệ Thống Điện, Trường Đại Học cần Thơ [5] TS Nguyễn Đăng Toản, Bài giảng môn học Phần Mầm Tính Tốn Hệ Thống Điện, 32, Trường Đại Học Điện Lực [6] Nguyễn Trọng Toàn & Vũ Anh Mỹ, Đề cương giảng môn học Phương Pháp Tỉnh Tốn số, Bộ Mơn Tốn, Khoa Cơng Nghệ Thơng Tin [7] Nguyễn Hoài Sơn, Đỗ Thanh Việt, Bùi Xuân Lâm (2002), ửng Dụng Matỉab Trong Tỉnh Toán Kỹ Thuật (Tập I), Nhà Xuất Bản Đại Học Quốc Gia Thảnh Phố Hồ Chí Minh Tiếng Anh [8] Prof J Nanda, P.R Bijwe, J Henry, V Bapi Raju, (March 1992) “General purpose fast decoupled power flow”, IEE PROCEEDINGS-C, Vol 139, No 2, pp 91 [9] Ray D Zimmerman, Hsiao-Dong Chiang, “Fast Decoupled Power Flow for Unbalanced Radial Distribution”, Systems School of Electrical Engineering ,Cornell University, Ithaca, NY 14853 USA, pp [10] B Stott and o Alsaẹ, “Fast Decoupled Load Flow”, IEEE Transactions on Power Apparatus and Systems, Vol PAS-93, May/June 1974 Sinh Viên TH: Phạm Ngọc Minh Trang 77 CBHD: Trần Anh Nguyện PHỤLỤC PHỤ LỤC ■ Danh sách biến Míĩle “Dulieudauvao.m” Bảng 6: Ỷ nghĩa biến MFile “Duỉieudauvao.m ” STT Tên Biến busdata linedata transdata -donvituongdoi -Ybus -FD -phanbocongsuat -xuatketqua Giải thích Mơ tả Chứa ứiông số Bus Chứa thông số đường dây Chứa thông số máy biến áp Liên kết liệu file “Dulieudauvao.m” với file “donvituongdoi.m” Liên kết liệu file “Dulieudauvao.m” với file “Ybus.m” Liên kết liệu file “Dulieudauvao.m” với file “FD.m” Liên kết liệu file “Dulieudauvao.m” với file “phanbocongsuatm” Liên kết liệu file “Dulieudauvao.m” với file “xuatketqua.m” ìs, Chưoug trình đổi sang đon vị tưtmg đối Busdata: %BUSDATA [ g , h ] = s i z e ( b u s d a t a ) ; b a s e = 0 %g=24,h=ll pubus= zeros(g,h); pubus( :,1 ) =busdata(:,1 ); pubus(:,1 )=busdata(:,11); for i= l:g p u b u s ( i , 2)= p u b u s ( i , ) = b u s d a t a ( i , ) * p u b u s ( i , 2) fo r j=3:9 p u b u s ( i ,j )= b u s d a ta ( i,j) /base; end end ìs, Chưtmg trình đổi sang đơn vị tưong đối Linedata: %LINEDATA [ g ,h ] = s i z e ( b u s d a t a ) ; [ v , w ] = s i z e ( l i n e d a t a ) ; %g=24,h=ll,v=18,w=9 p u lin e = z e r o s ( v ,7); for i= l:v f o r j = l :7 p u l i n e ( i , j )= l i n e d a t a ( i , j ); Sinh Viên TH: Phạm Ngọc Minh Trang 78 CBHD: Trần Anh Nguyện PHỤLỤC end end for i= l:v fo r x=l:g i f l i n e d a t a ( i , ) = = b u s d a t a ( x , 1) p u l i n e ( i , ) = ( l i n e d a t a ( i , ) * l i n e d a t a ( i , ) ) / ( b u s d a t a ( x , ) A2 / 0 ) ; p u l i n e ( i , ) = ( l i n e d a t a ( i , ) * l i n e d a t a ( i , ) ) / ( b u s d a t a ( x , ) A2 / 0 ) ; p u l i n e ( i , 6)=- ( ( l i n e d a t a ( i , ) / l i n e d a t a ( i , ) ) *le6) / ( b u s d a t a ( x , ) A2 / 0 ) ; %x u a t d i e n " l e " l a d o n v i c u a R s h u n t l a Mohm/km p u l i n e ( i , ) = ( l i n e d a t a ( i , ) ) / (b a s e * 0 /( s q r t(3 )* b u s d a ta (x ,2)));% x u a t h i e n 1000 l a d on v i c u a S c b l l e c o n c u a Ucb l a l e d o n g i a n v o i end end end ìs, Chương trình đổi sang đơn vị tưong đối Tranformerdata: %MBA [ t , s ] = s i z e ( t r a n s d a t a ) ; % t=ll,s=7 p u t r a n s = z e r o s ( t ,4); p u tra n s(:,1 :3 )= tran sd ata(:,1 :3 ); for i= l : t fo r x=l:g i f t r a n s d a t a ( i , ) = = b u s d a t a ( x , 1) p u t r a n s ( i , ) = ( t r a n s d a t a ( i , ) / 0 ) * ( ( t r a n s d a t a ( i , ) A2 ) / t r a n s d a t a ( i , ) ) * ( 0 / b u s d a t a ( x , ) A2 ) ; end end end B Danh sách biến MfQe “donvituongdoi.m” Bảng 7: Ỷ nghĩa biến MFile “donvituongdoỉ.m ’’ STT Tên Biến g h V w t s pubus puline Giải thích Mơ tả Số hàng busdata Số cột busdata Số hàng linedata Số cột linedata Số hàng transdata Số cột transdata Ma trận busdata đổi sang đơn vị tương đối Ma trận linedata đổi sang đơn vị Sinh Viên TH: Phạm Ngọc Minh Trang 79 CBHD: Trần Anh Nguyện PHỤLỤC putrans 10 11 base luly tương đôi Ma trận transdata đổi sang đơn vị tương đối Có giá trị 100 Chuỗi thông tin muốn hiển thị ■53 Chương trình tạo ma trận Ybus: clc % g = , h = l l , v = , w = , x = l l , ỵ=7 [ g , h ] = s i z e ( b u s d a t a ) ; [ v , w]= s i z e ( l i n e d a t a ) ; [ x , y ] = s i z e ( t r a n s d a t a ) ; Z d ata= zero s(v ,3 ) ; Y e=zeros(x,6 ) ;Y = zeros(g,g); Z = zero s(v ,1 ) ; Z sh = zero s(v ,1); for i= l:v z ( i , 1)= co m p le x ( p u l i n e ( i , ) , p u l i n e ( i , ) ) ; Zs h ( i , 1) = c o m p l e x ( , p u l i n e ( i , 6) ) ,ẵ Z d a t a ( i , 1) = l i n e d a t a ( i , 2) ; Z d a t a ( i , 2) = l i n e d a t a ( i , 3) ; Z d a t a ( i , 3) =z ( i , 1) ,ẵ end nl= Z data(:,1 ); nr=Z data(:,2 ); for i= l:v for j= l:3 ỵ ( i , j ) = lin e d a ta ( i , j ); end end for i= l:v f o r j = l :6 y (i,4 )= l/Z (i,l); y s h ( i , 1) = / Z s h ( i , 1) ; y ( i , 5) = y s h ( i , ) * ; y (i , ) = y (i , ) ; end end d = z e ro s (x ,2); for i= l:x for j= l:g i f t r a n s d a t a ( i , ) = = b u s d a t a ( j ,1) d ( i , ) = t r a n s d a t a ( i , ) / b u s d a t a (j , 2) ; end i f t r a n s d a t a ( i , ) = = b u s d a t a ( j ,1) d ( i , ) = t r a n s d a t a ( i , ) / b u s d a t a (j , ) ; end end Sinh Viên TH: Phạm Ngọc Minh Trang 80 PHỤLỤC CBHD: Trần Anh Nguyện end Z t= z e ro s (x ,1 ) ;ỵe=zeros (x ,1); tio = z e ro s (x ,1); for i= l:x r a t i o ( i , 1)= d ( i , ) / d ( i , ) ; end for i= l:x for j= l:3 Y e( i , j ) = t r a n s d a t a ( i , j ); end end for i= l:x z t ( i , 1)= c o m p le x (ũ , p u t r a n s ( i , ) ) ; y e ( i , ) = l / z t ( i , 1) ; Y e ( i , 4)= y e ( i , ) / r a t i o ( i , 1); Y e ( i , )= (1 -(1 / r a t i o ( i , ) ) ) * y e ( i , 1); Ye ( i , 6) = ( ( / ( r a t i o ( i , 1) A2) ) - ( / r a t i o ( i , 1) ) ) *ye ( i , 1) ; end % Ma t r a n Ybus n g o a i d u o n g c h e o c h i n h for i= l:v i f n l (i) > & n r (i) > Y ( n l ( i ) , n r ( i ) ) = Y ( n l ( i ) , n l ( i ) ) - ỵ ( i , 4) ; Y ( n r ( i ) , n l ( i ) ) = Y ( n l ( i ) , n r ( i ) ); end end for j= l:x i f Y e (j , ) > & Y e (j , ) > Y (Y e(:,2) ,Y e ( j,3 ) )=Y(Ye(j,2) ,Y e ( j,3 ) ) - Y e ( j,4 ) ; Y (Y e(:,3) ,Y e ( :,2 ) )=Y(Ye(:,2) ,Y e ( :,3 ) ) ; end end % Ma t r a n Ybus t r e n d u o n g c h e o c h i n h for i= l:g for j= l:x i f Y e (j , ) = = i Y ( i , i ) = Y ( i , i ) + Y e (j , ) + Y e (j , 6) ; e l s e i f Y e(j,3)== i Y ( i , i ) = Y ( i , i ) + Y e (j , ) + Y e (j , 5) ; end end end for i= l:g for j= l:v i f y (j,2)= = i Y ( i , i ) =Y ( i , i ) +y (j , 4) +y (j , 5) ; e ls e if y (j,3 )= = i Y ( i , i ) =Y ( i , i ) +y (j , 4) +Y (j , 5) ; Sinh Viên TH: Phạm Ngọc Minh Trang 81 CBHD: Trần Anh Nguyện PHỤLỤC end end end ■ Danh sách biến Mfile “Ybus.m” Bảng 8: Ỷ nghĩa biến MFile “Ybus.m ” STT Tên Biến Zdata Ye Y z Zsh y ysh d zt 10 nl 11 nr 12 ratio Giải thích Mơ tả Ma trận chứa giá trị: số thứ tự, tò bus đầu đến bus cuối đường dây ma trận linedata Ma trận chứa giá trị: số thứ tự, từ bus đầu đến bus cuối, điện kháng, dung dẫn hai phía MBA Ma trận Ybus hệ thống điện Ma trận chứa số phức điện trở điện kháng ma trận puline Ma trận chứa số phức dung kháng ma trận puline Ma trận chứa giá trị: số thứ tự, từ bus đầu đến bus cuối, điện dẫn (1/Z), dung dẫn (ysh) đường dây Ma trận dung dẫn (1/Zsh) đường dây Ma trận chứa giá trị: tỷ số (điện áp sơ cấp MBA điện áp bus phía trước MBA) tỷ số (điện áp thứ cấp MBA điện áp bus phía sau MBA) Ma trận chứa giá trị số phức điện kháng MBA Ma trận chứa giá trị cột thứ ma trận Zdata Ma trận chứa giá trị cột thứ ma trận Zdata Ma trận chứa tỷ số cột ma trận d Sinh Viên TH: Phạm Ngọc Minh Trang 82 CBHD: Trần Anh Nguyện PHỤLỤC Chương trình tính tốn phân bố công suất tổn thất hệ thống điện phương pháp Fast Decoupled: Ì3, clear V anglel d e l t a P d e l t a Q deltaangle deltav p jj j c l e a r P s c h Q s c h Bp Bq i t e r B p l B q l Pk Qk Pg Pd Qg Qd DP PI Q1 Ps Qs nbus=length(busdata( : , ) ) ; n b r= le n g th (lin e d a ta ( : , ) ) ; [g ,h ]= size(b u sd ata); base=10ũ; in d ex = l; nb=0; n lo a d = ; format s h o rt g Ybus=abs(Y); t= a n g le ( Y ) ; %%%Vong l a p b a c dau%%%%%% P s c h = [ ] ; P g = [ ] ; P d = [ ] ; Q s c h = [ ] ; Qg=[ ] ; Qd=[ ] ; k b = []; for i= l:g i f p u b u s ( i , ) =acc%%%%%%Su l a p d i l a p l a i dau%%%%%% i t e r = i t e r + l ; %%%Tính P ( k ) v a Q(k)%%% chuan b i bac for i= l:g i f k b (i )==0 I k b (i )==2 n b = n b + l;end, i f k b ( i ) = = n l o a d = n l o a d + l ; end end B pl=inv(B p);B ql=inv(B q); j=0; for i= l:g i f k b (i )==0 I k b (i )==2 j= j+ l; Ps (j , ) = P s c h ( i , l ) ; e n d end j=0; for i= l:g i f kb(i)==0 j= j+l; Qs ( j , ) = Q s c h ( i , ) ; e n d end n n = l e n g t h ( P s ) ; P x = z e r o s ( g , ) ; %%%Tinh p ( k ) %%% for i= l:g for j= l:g P ( i , l ) = v ( i , l ) * v ( j , l ) * Y b u s ( i , j ) * c o s ( t ( i , j ) - a n g l e l ( i , ) + a n g l e l (j , ) ) ; Px ( i , 1) =Px ( i , 1) +p ( i , 1) ; end end P l = z e r o s ( n n ,1 ) ; j=0; for i= l:g i f kb ( i ) ==0 I kb ( i ) ==2 j= j+ l; P1 ( j , l ) = P x ( i , l ) ,Ệ V I(j,l)= v (i,l); end end f o r i = l : nn%%%Tinh d e l t a P ( k ) %%%% d e l t a P ( i , 1) =Ps ( i , l ) - P l ( i , l ) , ẵ d e c o y ( i , 1) = d e l t a P ( i , 1) / V I ( i , 1) ,ẵ end Sinh Viên TH: Phạm Ngọc Minh Trang 84 CBHD: Trần Anh Nguyện PHỤLỤC d e lta a n g le = - B p l* decoy j=0; for i= l:g i f k b (i )==0 I k b (i )==2 j= j+ l; a n g l e l ( i , ) = a n g l e l ( i , ) + d e l t a a n g l e (j , ) ; end end n u n = l e n g t h ( Q s ) ; Q x = z e r o s ( g , ) ; %%%Tinh Q ( k ) %%% for i= l:g for j= l:g Q( i , ) = - ( v ( i , )* v ( j , ) * Y b u s (i,j) * s i n ( t ( i , j ) a n g l e l ( i , 1) + a n g l e l (j , 1) ) ) ,ẵ Q x ( i , 1)= Q x (i,1 )+ Q ( i,1); end end j=0; for i= l:g i f kb(i)==0 j= j+l; Q1 ( j , 1) =Qx ( i , 1) ; V (j,l)= v (i,l) ; end end f o r i = l :nun%%%Tinh del taQ( k) %%% d e l t a Q ( i , 1) =Qs ( i , 1) -Q1 ( i , 1) ; d e c o y l ( i , 1)= d e l t a Q ( i ,1 ) /V 2( i , 1); end deltav=-B ql*decoyl; j=0; f o r i i = l : g%%%Thay t h e g i a t r i moi c u a " v " %%% i f k b ( i i ) ==0 j=j+l/ẳ V (ii,1)=v(ii,1)+ d elt a v (j , 1); end end m axerror=m ax(m ax(abs(deltaP)) ,m ax (ab s(d eltaQ ))); nb=0; n lo a d = ; end Q g en= zeros(g,1); for i= l:g i f k b (i )==2 Q g e n ( i , ) = Q x ( i , 1) + Q d ( i , 1) ; Sinh Viên TH: Phạm Ngọc Minh Trang 85 CBHD: Trần Anh Nguyện PHỤLỤC i f Q g e n ( i , ) < = p u b u s ( i , 6) Q g ( i , ) = p u b u s (i , ) ; kb ( i ) =0; index=l; end i f Q g e n ( i , ) > = p u b u s ( i , 7) Q g ( i , ) = p u b u s ( i , 7) ; kb ( i ) =0; index=l; end i f Q g e n ( i , ) > p u b u s ( i , 6)& Q g e n