Phương pháp điểm nội tối ưu phân bố công suất trong hệ thống điện AC DC song song

Đại Học Quốc Gia Tp Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA TRẦN QUỐC TUẤN PHƯƠNG PHÁP ĐIỂM NỘI TỐI ƯU PHÂN BỐ CÔNG SUẤT TRONG HỆ THỐNG ĐIỆN AC/DC SONG SONG Chuyên ngành : Thiết bị, mạng nhà máy điện Mã Ngành: 60 52 50 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng năm 2008 CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học : TS Vũ Phan Tú (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét : (Ghi rõ họ, tên, học hàm, học vị chữ ký) Cán chấm nhận xét : (Ghi rõ họ, tên, học hàm, học vị chữ ký) Luận văn thạc sĩ bảo vệ tại: HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA Ngày tháng năm TRƯỜNG ĐẠI HỌC BÁCH KHOA PHÒNG ĐÀO TẠO SĐH CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM ĐỘC LẬP – TỰ DO – HẠNH PHÚC Tp HCM, ngày tháng năm 200 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: TRẦN QUỐC TUẤN Phái: Nam Ngày, tháng, năm sinh: 30-09-78 Nơi sinh: Quảng Ngãi Chuyên ngành: Thiết Bị Mạng Và Nhà Máy Điện MSHV: 01806513 I- TÊN ĐỀ TÀI: PHƯƠNG PHÁP ĐIỂM NỘI TỐI ƯU PHÂN BỐ CÔNG SUẤT TRONG HỆ THỐNG ĐIỆN AC/DC SONG SONG II- NHIỆM VỤ VÀ NỘI DUNG: Tính phân bố công suất mạng AC/DC với ma trận Jacobi đầy đủ Tính tối ưu phân bố cơng suất mạng AC/DC song song phương pháp điểm nội đối ngẫu dự đoán hiệu chỉnh Xây dựng giao diện cho phân tích mơ hệ thống phần mềm Matlab III- NGÀY GIAO NHIỆM VỤ : 21 – 01 – 2008 IV- NGÀY HOÀN THÀNH NHIỆM VỤ: 30 – 06 – 2008 V- CÁN BỘ HƯỚNG DẪN (Ghi rõ học hàm, học vị, họ, tên): TS VŨ PHAN TÚ CÁN BỘ HƯỚNG DẪN (Học hàm, học vị, họ tên chữ ký) CN BỘ MÔN QL CHUYÊN NGÀNH Nội dung đề cương luận văn thạc sĩ Hội đồng chun ngành thơng qua TRƯỞNG PHỊNG ĐT – SĐH Ngày tháng năm TRƯỞNG KHOA QL NGÀNH (Ghi chú: Học viên phải đóng tờ nhiệm vụ vào trang tập thuyết minh LV) LỜI CẢM ƠN Cảm ơn Ba, Mẹ sinh thành nuôi nấng dưỡng dục con, cho trí tuệ sáng suốt, ln ln bên ủng hộ đường học vấn Cảm ơn người bạn lớp Cao Học Thiết Bị Mạng Nhà Máy Điện K2006, em Công Quyền, Lê Phương, bạn Hữu Tồn đóng góp lời nhận xét q báu Xin cảm ơn quan tơi đặc biệt ban lãnh đạo Công ty thủy điện Đa Nhim Hàm Thuận Đa Mi, Phi Anh, anh Oánh tạo điều kiện cho trau dồi nâng cao kiến thức Xin cảm ơn người bạn đồng nghiệp giúp đỡ đặc biệt vợ chồng anh Minh Thắng Cảm ơn Pinkit! Em người cho anh nghị lực niềm tin, bên anh, ủng hộ giúp đỡ anh lúc anh gặp khó khăn Với tất lịng chân thành, xin cảm ơn q thầy trường Đại Học Bách Khoa Tp Hồ Chí Minh, đặc biệt TS Vũ Phan Tú người thầy tận tình bảo, hướng dẫn, đề xuất ý tưởng hay để em hoàn thành luận văn Tp Hồ Chí Minh, ngày 30/07/2008 Trần Quốc Tuấn TÓM TẮT Truyền tải HVDC xu phát triển tập đoàn điện lực toàn giới kỷ 21, nhằm liên kết vùng, lãnh thổ hay quốc gia lại với làm tăng hiệu sử dụng nguồn, tăng độ tin cậy Nhằm giải toán quản lý vận hành tối ưu hệ thống, Luận Án giải toán tối ưu phân bố công suất mạng AC/DC song song Trước tiên Luận Án trình bày sơ lược lịch sử phát triển mạng truyền tải HVDC, so sánh điểm mạnh yếu với hệ thống HVAC Luận Án giới thiệu cơng trình xây dựng tiêu biểu qua cho nhìn tổng quan hệ thống truyền tải HVDC Để tìm hiểu sâu hệ thống HVDC, Luận Án mô tả chi tiết phận cấu thành nên trạm truyền tải HVDC chuyển đổi, máy biến áp, lọc, đường dây…, trình bày mơ hình liên kết hệ thống truyền tải HVDC HVAC chủ yếu cho tính tốn phân bố cơng suất Đồng thời Luận Án đưa mơ hình điều khiển cho đường dây HVDC, thành lập tính ma trận Jacobi đầy đủ cho tính tốn phân tích phân bố cơng suất Tiếp theo lý thuyết toán tối ưu Luận Án trình bày sơ lược phương pháp giải, sau lập mơ hình tốn học tổng qt giải cho toán tối ưu phương pháp điểm nội đối ngẫu dự đoán hiệu chỉnh Áp dụng lý thuyết giải giải toán tối ưu phân bố công suất hệ thống điện với hàm đối tượng cực đại hố lợi nhuận Cuối Luận Án mơ tả chương trình phân tích hệ thống mơ hình mơ phần tử hệ thống điện, mơ hình liệu cho phân tích tối ưu mơ hình cho phần tử khác Phần mơ tả bước tính giải thuật Newton Raphson cho phân bố công suất giải thuật điểm nội cho tính tốn tối ưu xây dựng chương trình Chữ viết tắt CHỮ VIẾT TẮT Ký hiệu Tên Tiếng Anh AC Alternate Current ATC Available transfer capacity CB Circiuit Breaker CC Constant Current Control CCC Capacitor Commutated Converter CEA Constant Extintion Angle CIA Constant Ignition Angle DAS Dual Affine Scaling DC Direct Current DG Distributed Generation FACTS Flexible AC Transmission Systems GRTS Ground Return Transfer Switch GTO Gate Turn-Off Thyristor HPP Hydro power plant HSNBS High-Speed Neutral Bus Switch HSGS High-Speed Ground Switch HVDC High Voltage Direct Current IGBT Insulated Gate Bipolar Transistor IPM Interior Point Method KKT Karush – Kuhn – Tucker LP Linear Program LVDC Low voltage direct current MRTB Metallic Return Transfer Breaker MVDC Medium Voltage Direct Current NR Newton Raphson OH Over Head i Chữ viết tắt OPF Optimal Power Flow PDIPA Primal-Dual Interior Point Algorithm PCPDIPA Predictor-Corrector Primal-Dual Interior Point Algorithm PF Power Flow PI Proportional Integral PSAT Power System Analysis Toolbox PWM Pulse Width Modulation QP Quadratic Program RoW Right of Way STATCOM Static Compensator VDCOL Voltage Dependent Current Order Limit VSC Voltage Source Converter ii Mục Lục MỤC LỤC Chương I:   TỔNG QUAN   I.  ĐẶT VẤN ĐỀ 5  II.  SƠ LƯỢC MẠNG TRUYỀN TẢI CAO ÁP MỘT CHIỀU 5  II.1.  LỊCH SỬ PHÁT TRIỂN MẠNG TRUYỀN TẢI HVDC 5  II.2.  MỘT SỐ ĐẶC TRƯNG – THUẬN LỢI CỦA HVDC 7  II.3.  CÁC CƠNG TRÌNH XÂY DỰNG TIÊU BIỂU 14  III.  PHẠM VI NGHIÊN CỨU – ĐIỂM MỚI 19  III.1.  ĐỀ TÀI LIÊN QUAN 19  III.2.  PHẠM VI NGHIÊN CỨU 21 III.3.  ĐIỂM MỚI 21 Chương II: MƠ HÌNH TRUYỀN TẢI HVDC                                                                    LIÊN I.  KẾT VỚI HVAC     TRẠM VÀ ĐƯỜNG DÂY TRUYỀN TẢI HVDC 22  I.1.  CÁC CẤU HÌNH TRẠM TRUYỀN TẢI HVDC 22  I.1.1.  KẾT NỐI ĐƠN CỰC 22  I.1.2.  KẾT NỐI LƯỠNG CỰC 22  I.1.3.  KẾT NỐI ĐỒNG CỰC 24  I.1.4.  CÁC KIỂU LIÊN KẾT BỘ CHUYỂN ĐỔI 24  I.2.  CÁC THÀNH PHẦN CHÍNH CỦA MỘT HỆ THỐNG HVDC 26  I.2.1.  TRẠM TRUYỀN TẢI 26  I.2.2.  PHƯƠNG TIỆN TRUYỀN TẢI 29  I.2.3.  ĐIỆN CỰC 31  II.  MƠ HÌNH HỆ THỐNG HVDC 33  II.1.  II.1.1.  MƠ HÌNH BỘ BIẾN ĐỔI 33  NGUYÊN LÝ HOẠT ĐỘNG 33  Trang Mục Lục II.1.2.  II.2.  II.2.1.  QUAN HỆ GIỮA ĐẠI LƯỢNG A.C VÀ D.C 35  ĐẶC TÍNH ĐIỀU KHIỂN CỦA HỆ THỐNG HVDC 37  ĐẶC TÍNH ĐIỀU KHIỂN CƠ BẢN: 37  II.2.1.1.  Sự chọn lựa cơ bản cho điều khiển 37  II.2.1.2.  Đặc tính điều khiển 38  II.2.2.  ĐẶC TÍNH ĐIỀU KHIỂN BỔ SUNG 40  II.2.2.1.  Giới hạn dòng 40  II.2.2.2.  Giới hạn góc kích 41  II.2.2.3.  Tóm tắt nguyên tắc điều khiển 41  III.  LỜI GIẢI PHÂN BỐ CÔNG SUẤT 42  III.1.  CÁC PHƯƠNG TRÌNH LIÊN HỆ GIỮA HVDC VÀ HVAC 42  III.2.  BIỂU DIỄN LỜI GIẢI PHÂN BỐ CÔNG SUẤT 44  III.3.  BIỂU DIỄN CÁC ĐẠI LƯỢNG TRONG HỆ ĐƠN VỊ 46  Chương III:                       BÀI TOÁN PHÂN BỐ CƠNG SUẤT  I.  CƠ SỞ TỐN HỌC 47  II.  BÀI TỐN PHÂN BỐ CƠNG SUẤT TRONG MẠNG AC 50  II.1.  CÁC PHƯƠNG TRÌNH LIÊN HỆ 50  II.2.  GIẢI BÀI TOÁN PHÂN BỐ CÔNG SUẤT 51  II.3.  LƯU ĐỒ GIẢI THUẬT 53  III.  BÀI TỐN PHÂN BỐ CƠNG SUẤT TRONG MẠNG AC/DC 55  III.1.  MƠ HÌNH BỘ ĐIỀU KHIỂN VÀ ĐƯỜNG DÂY HVDC 55  III.2.  LỜI GIẢI BÀI TOÁN PHÂN BỐ CÔNG SUẤT 57  III.3.  GIẢI THUẬT 59  Chương IV:            BÀI TOÁN TỐI ƯU PHÂN BỐ CƠNG SUẤT                                                                  CỰC I.  ĐẠI HĨA LỢI NHUẬN  BÀI TỐN TỐI ƯU PHÂN BỐ CÔNG SUẤT 60  I.1.  LỊCH SỬ PHÁT TRIỂN VÀ LÝ THUYẾT CƠ BẢN 60  I.2.  CÁC PHƯƠNG PHÁP GIẢI BÀI TOÁN OPF 62  Trang Mục Lục I.2.1.  SƠ LƯỢC CÁC PHƯƠNG PHÁP 62  I.2.2.  PHƯƠNG PHÁP ĐIỂM NỘI 66  I.2.2.1.  Tỷ lệ affine (PAS) 67  I.2.2.2.  Phương pháp điểm nội đối ngẫu affine (DAS) .68  I.2.2.3.  Phương pháp điểm nội sử dụng bước tiến Newton .70  I.2.2.4.  Phương pháp điểm nội thừa số chặn logarithm (Barrier) .72  I.2.2.5.  Phương pháp điểm nội thừa số Lagrange .72  I.2.2.6.  Phương pháp điểm nội đối ngẫu dự đoán hiệu chỉnh 72  II.  HÀM ĐỐI TƯỢNG – CỰC ĐẠI HOÁ LỢI NHUẬN 79  II.1.  HÀM CHI PHÍ PHÁT VÀ GIAO DỊCH ĐIỆN NĂNG 79  II.2.  HÀM CỰC ĐẠI CÔNG SUẤT TRUYỀN TẢI 80  II.3.  HÀM ĐỐI TƯỢNG ĐA MỤC TIÊU 81  II.4.  HÀM LAGRANGE TỔNG QUÁT 82  Chương V:  I.             CHƯƠNG TRÌNH PHÂN TÍCH HỆ THỐNG  MƠ HÌNH CÁC PHẦN TỬ TRONG HỆ THỐNG 83  I.1.  MƠ HÌNH HỆ THỐNG HVDC 83  I.2.  DỮ LIỆU CHO BÀI TOÁN TỐI ƯU 85  I.2.1.  NGUỒN CUNG CẤP 86  I.2.2.  MÁY PHÁT DỰ TRỮ 87  I.2.3.  ĐỘ TĂNG - GIẢM CỦA CÔNG SUẤT NGUỒN 88  I.2.4.  PHỤ TẢI 88  I.2.5.  ĐỘ TĂNG - GIẢM CỦA PHỤ TẢI 89  I.2.6.  ĐỒ THỊ PHỤ TẢI 90  II.  CHƯƠNG TRÌNH PHÂN BỐ CƠNG SUẤT 92  II.1.  CÀI ĐẶT CHUNG 92  II.2.  TRÌNH TỰ XỬ LÝ 93  III.  CHƯƠNG TRÌNH TỐI ƯU PHÂN BỐ CƠNG SUẤT 95  III.1.  CÀI ĐẶT VÀ GIAO DIỆN NGƯỜI DÙNG 95  III.2.  TRÌNH TỰ XỬ LÝ 96  Trang Chương V Chương trình phân tích + Gọi hàm fm_wcall, viết hàm fm_call cho phần tử gọi + Lấy bus chuẩn (bus cân bằng) + Gọi hàm fm_dynlf định nghĩa biến trạng thái cho phương trình phân bố cơng suất + Kiểm tra phương pháp tính phân bố cơng suất NR hay FD + Kiểm tra slack bus hay nhiều slack bus + Kiểm tra phần tử liên kết lưới hay khơng + Set thời gian tính tốn vịng lặp While sai số > sai số cho phép số vòng lặp < vòng lặp max if NR Gọi hàm refreshJac xoá ma trận Jacobi if hệ thống nhiều nút cân • Gọi hàm fm_call ‘kgpf’ tính tốn ma trận J hàm f • Giải phương trình J∆x = -∆f • Cập nhật thơng số elseif setting.robust else phương pháp NR • Gọi hàm fm_call ‘l’ tính tốn ma trận J hàm f • Refresh ma trận Jacobi • Giải phương trình J∆x = -∆f • Cập nhật thông số end else FD + + + + + + + + end ++ tăng vòng lặp ++ cập nhật thông số ++ kiểm tra sai số lớn hiển thị thơng báo end Tính thời gian tính tốn Tính tổng cơng suất tiêu thụ (Pl,Ql) bơm vào hệ thống (Pg, Qg) Cấp nhớ cho biến trạng thái biến động Hiển thị kết thời gian tính tốn Nếu hiển thị kết phân bố cơng suất gọi hàm giao diện fm_stat Kích hoạt tất phương trình refresh ma trận Jacobi, gọi hàm fm_call ‘i’ để đặt điểm ban đầu Xây dựng cấu trúc “Snap Shot” Hiển thị thông báo lỗi hội tụ gọi hàm fm_status (‘pf’, ‘close’) B4: Xuất kết phân bố công suất dạng file.txt hay trực tiếp mơ hình simulink hiển thị biên độ điện áp, góc pha, giá trị (p.u) chiều dịng chảy cơng suất đường dây, máy phát, phụ tải từ giao diện Simulink Settings Trang 94 Chương V Chương trình phân tích Hiển thị khơng hiển thị giá trị biên độ điện áp góc pha Hiển thị không hiển thị giá trị chiều công suất Hiển thị không hiển thị tên máy phát III CHƯƠNG TRÌNH TỐI ƯU PHÂN BỐ CÔNG SUẤT III.1 CÀI ĐẶT VÀ GIAO DIỆN NGƯỜI DÙNG ¾ Giao diện phân tích tối ưu phân bố công suất Newton Direction Mehrotra Predictor/Corrector Currents Active Powers Apparents powers Single OPF Pareto Set Dayly Forecast ACT (by CPF) ACT (by sentivity analysis) Flat start Actual power flow solution ¾ Thơng số cài đặt Kí hiệu Giải thích Weighting Factor [0,1] omega (ω): Trọng số chọn hàm đối tượng Min Loading Parameter lambda (λmin): Hệ số tải nhỏ Max Loading Parameter lambda (λmax): Hệ số tải lớn Time Interval deltat (∆t): Thời gian tính tích phân OP Tolerance esp2 (ε2): Sai số cho phân tích tối ưu PF eqs Tolerance esp1 (ε1): Sai số tính tốn phân bố cơng suất mu Tolerance eps_mu (εµ): Sai số cho thừa số chặn Safety Factor gamma (γ): hệ số an tồn Trang 95 Chương V Chương trình phân tích Centering Parameter beta (β): hệ số trung tâm Solving Method: 1: Newton Direction: phương pháp bước tiến newton 2: Mehrotra Predictor/Corrector: phương pháp dự đoán hiệu chỉnh theo Mehrotra Flow Limit 1: Currents: giới hạn dòng điện (giới hạn) 2: Active Powers: giới hạn công suất tác dụng 3: Apparents powers: giới hạn công suất biểu kiến OPF type 1: Single OPF tốn tối ưu thơng thường 2: Pareto Set sử dụng tập Pareto (loại 3: Dayly Forecast dự báo ngày toán OPF) 4: ACT (by CPF): khả truyền tải dùng phân tích CPF 5: ACT (by sentivity analysis): ACT phân tích độ nhạy Initial Guess 1: Flat start : Chọn V = 1, góc pha zero (trị ban đầu) 2: Actual power flow solution: kế thừa giá trị điện áp, góc pha, cơng suất có sau thực tốn phân bố cơng suất III.2 TRÌNH TỰ XỬ LÝ B1: Sau load file liệu cần phân tích, từ cửa sổ giao diện menubar > chọn OPF settings, chương trình gọi hàm fm_opffig hiển thị giao diện cho phân tích tối ưu Settings thơng số cần cho phân tích B2: Nhấn nút Optimal PF, chương trình gọi hàm fm_set opf Nếu lambda > lambda max thơng báo lambda phải nhỏ lambda max switch OPF type Kiểm tra xem omega lambda số vô hướng vector case Sử dụng phần tử thứ omega lambda if omega == Gọi hàm fm_opfm: hàm tối ưu cực đại hoá lợi nhuận else Gọi hàm fm_opfsdr: hàm tối ưu đa mục tiêu ràng buộc an ninh điện áp end case if omega số vô hướng if omega ==0 Gọi hàm fm_opfm else Gọi hàm fm_opfsdr end else Trang 96 Chương V Chương trình phân tích gán hàm OPF.fun = ‘fm_opfsdr’ gọi hàm fm_pareto end case báo chưa nghiên cứu dự báo ngày ! case gọi hàm fm_atc case gọi hàm fm_atc end Hàm fm_opfm tối ưu cực đại hoá lợi nhuận bắt đầu đếm thời gian if initial guess ==1 gán điều kiện đầu cho điện áp công suất kháng máy phát else lấy liệu từ phân bố cơng suất end - Kiểm tra có chọn hệ máy phát, phụ tải không - Đọc thơng số nguồn, phụ tải, nguồn dự phịng - Kiểm tra giới hạn điện áp, công suất kháng máy phát, dịng điện đường dây - Tính ma trận Jacobi - Chọn điểm ban đầu While - Tính hàm f(y) (chi tiết xem phụ lục B) - Điều kiện KKT: Tính đạo hàm riêng hàm Lagrange theo biến s, µ, y - Tính ma trận Hessian Switch phương pháp Case NR giải phương trình tìm hướng tiến ∆y, ∆s, ∆µ Case Mehrotra + Xem µs thừa số chặn zero, ước lượng NR cho µ, thay vào phương trình giải tìm đại lượng ước lượng ∆y’, ∆s’, ∆µ’ + Dùng bước centering corector ước lượng tính µs theo ∆(y’, s’, µ’) + Tính đại lượng hiệu chỉnh gs thay đại lượng hiệu chỉnh gs vào phương trình (B.2) giải tìm hướng tiến ∆y, ∆s, ∆µ thực End - Tính độ dài bước alpha theo (4.40) - Cập nhật biến - Giảm µs theo (4.37) - Kiểm tra tiêu chuẩn hội tụ - Vòng lặp > vòng lặp max, break End Xuất kết Trang 97 Chương VI Kết - nhận xét CHƯƠNG VI: KẾT QUẢ - NHẬN XÉT Chương trình bày kết thu sau kiểm tra với mạng điện 24 nút 118 nút chuẩn IEEE Kết chi tiết xem phụ lục D đánh giá so sánh kết hai trường hợp phân bố công suất tối ưu phân bố công suất thường Kết so sánh với chương trình tính toán Matpower trung tâm nghiên cứu hệ thống điện (PSERC) I KẾT QUẢ KIỂM TRA MẠNG AC 118 NÚT Tối ưu phân bố công suất phương pháp IP Phân bố công suất thường Bus # V [p.u.] Bus 0.955 - Bus 0.9714 Bus Chương trình phân tích Ps Qg [p.u.] [p.u.] Matpower V [p.u.] Ps [MW] Qg [Mvar] V [p.u.] Ps [MW] Qg [Mvar] -0.03 1.0331 26.1 15 1.033 26.48 15 - - 1.0388 - - 1.039 - - 0.9677 - - 1.038 - - 1.038 - - Bus 0.998 - -0.128 1.06 - 70.653 1.06 71.52 Bus 1.0018 - 0.4015 1.0575 - - 1.057 - - Bus 0.99 - 0.1627 1.0528 - 33.332 1.053 0.02 33.28 Bus 0.9893 - - 1.0511 - - 1.051 - - Bus 1.015 - 0.5842 1.0413 - -78.35 1.041 -78.64 Bus 1.0429 - - 1.06 - - 1.06 - - Bus 10 1.05 4.5 -0.51 1.0527 401.65 -100.7 1.053 401.87 -99.67 Bus 11 0.985 - - 1.046 - - 1.046 - - Bus 12 0.99 0.85 0.9169 1.05 85.783 54.537 1.05 85.79 54.44 Bus 13 0.9683 - - 1.0345 - - 1.034 - - Bus 14 0.9836 - - 1.0489 - - 1.049 - - Bus 15 0.97 - 0.0821 1.0487 21.088 24.096 1.049 20.88 23.85 Bus 16 0.9838 - - 1.0462 - - 1.046 - - Bus 17 0.9947 - - 1.06 - - 1.06 - - Bus 18 0.973 - 0.2915 1.0504 13.444 30.592 1.05 13.22 30.47 Bus 19 0.962 - -0.142 1.0479 21.772 24 1.048 21.58 23.96 Bus 20 0.9568 - - 1.0362 - - 1.036 - - Bus 21 0.9576 - - 1.0313 - 1.031 - - Trang 98 Chương VI Kết - nhận xét Bus V Ps Qg V Ps Qg V Ps Qg # [p.u.] [p.u.] [p.u.] [p.u.] [MW] [Mvar] [p.u.] [MW] [Mvar] Bus 22 0.9688 - - 1.0344 - - 1.035 - - Bus 23 0.9994 - - 1.0485 - - 1.049 - - Bus 24 0.992 - -0.149 1.0459 - -7.679 1.046 -7.3 Bus 25 1.05 2.2 0.5185 1.06 193.63 -47 1.06 193.81 -47 Bus 26 1.015 3.14 0.0694 1.0278 279.7 -27.45 1.027 279.76 -27.58 Bus 27 0.968 - 0.0411 1.0411 9.6079 25.311 1.041 9.92 25.53 Bus 28 0.9616 - - 1.0346 - - 1.035 - - Bus 29 0.9632 - - 1.0352 - - 1.035 - - Bus 30 0.9867 - - 1.0325 - - 1.031 - - Bus 31 0.967 0.07 0.3281 1.0384 7.2455 27.224 1.039 7.25 27.19 Bus 32 0.963 - -0.161 1.0408 14.62 21.815 1.041 14.86 21.83 Bus 33 0.9707 - - 1.0462 - - 1.046 - - Bus 34 0.984 - -0.291 1.0561 6.0714 -7.054 1.056 4.88 -7.9 Bus 35 0.9804 - - 1.0543 - - 1.054 - - Bus 36 0.98 - 0.0855 1.0547 11.723 21.206 1.054 10.66 19.74 Bus 37 0.9902 - 0.2451 1.06 - - 1.06 - - Bus 38 0.964 - - 1.0179 - - 1.015 - - Bus 39 0.9698 - - 1.0427 - - 1.042 - - Bus 40 0.97 - 0.2891 1.0429 49.763 29.074 1.042 49.32 28.75 Bus 41 0.9668 - - 1.0358 - - 1.035 - - Bus 42 0.985 - 0.4132 1.041 40.946 22.312 1.041 40.99 22.1 Bus 43 0.977 - - 1.0401 - - 1.04 - - Bus 44 0.9842 - -0.097 1.0309 - - 1.031 - - Bus 45 0.9862 - -0.097 1.0253 - - 1.025 - - Bus 46 1.005 0.19 -0.149 1.0378 19.029 -4.039 1.038 19.04 -4.15 - - 1.045 - - 1.045 - - Bus 48 1.0206 - -0.156 1.0473 - - 1.047 - - Bus 49 1.025 2.04 1.1699 1.0496 193.21 19.613 1.05 193.33 19.69 Bus 50 1.0011 - - 1.0393 - - 1.039 - - Bus 51 0.9668 - - 1.0224 - - 1.022 - - Bus 52 0.9568 - - 1.0168 - - 1.017 - - Bus 53 0.946 - - 1.0165 - - 1.016 - - Bus 54 0.955 0.48 0.0395 1.0315 49.557 34.767 1.031 49.54 34.86 Bus 55 0.952 - 0.0466 1.0315 32.314 19.89 1.031 32.13 20.01 Bus 56 0.954 - -0.023 1.0312 32.713 15 1.031 32.56 14.89 Bus 57 0.9706 - - 1.032 - - 1.032 - - Bus 58 0.959 - - 1.0241 - - 1.024 - - Bus 59 0.985 1.55 0.8027 1.0473 149.85 95.431 1.047 149.7 95.67 Bus 60 0.9931 - - 1.0467 - - 1.046 - - Bus 47 1.017 Trang 99 Chương VI Kết - nhận xét Bus V Ps Qg V Ps Qg V Ps Qg # [p.u.] [p.u.] [p.u.] [p.u.] [MW] [Mvar] [p.u.] [MW] [Mvar] Bus 61 0.995 1.6 -0.419 1.0491 148.38 29.194 1.048 148.41 31.96 Bus 62 0.998 - 0.0127 1.0446 - 0.5116 1.044 0.89 Bus 63 0.9697 - - 1.0174 - - 1.016 - - Bus 64 0.9842 - - 1.0245 - - 1.023 - - Bus 65 1.005 3.91 0.641 1.0188 352.3 -67 1.016 352.24 -67 Bus 66 1.05 3.92 0.0685 1.06 348.49 -67 1.06 348.86 -67 Bus 67 1.0197 - - 1.0469 - - 1.047 - - Bus 68 1.0038 - - 1.0183 - - 1.015 - - Bus 69 1.035 5.14 -0.712 1.06 452.65 -112.5 1.06 453.67 -111.1 Bus 70 0.984 - 0.0998 1.0392 - 29.25 1.039 28.88 Bus 71 0.9868 - - 1.0394 - - 1.039 - - Bus 72 0.98 - -0.111 1.0396 - -5.206 1.04 -5.2 Bus 73 0.991 - 0.0966 1.0382 - -2.215 1.038 -2.22 Bus 74 0.958 - -0.166 1.0217 16.253 1.022 16.93 Bus 75 0.9673 - - 1.024 - - 1.024 - - Bus 76 0.943 - 0.0528 1.0116 22.514 23 1.012 22.85 23 Bus 77 1.006 - 0.1218 1.047 - 70 1.047 70 Bus 78 1.0034 - - 1.0423 - - 1.042 - - Bus 79 1.0092 - -0.204 1.0437 - - 1.044 - - Bus 80 1.04 4.77 1.1086 1.06 431.08 -10.48 1.06 430.84 -9.43 Bus 81 0.9984 - - 1.0142 - - 1.011 - - Bus 82 0.9885 - -0.195 1.0377 - - 1.038 - - Bus 83 0.9844 - -0.097 1.04 - - 1.04 - - Bus 84 0.9797 - - 1.0427 - - 1.043 - - Bus 85 0.985 - -0.056 1.0505 - 23 1.051 23 Bus 86 0.9867 - - 1.0443 - - 1.044 - - Bus 87 1.015 0.04 0.1102 1.0571 3.6261 3.5405 1.057 3.63 3.54 Bus 88 0.9875 - - 1.0481 - - 1.048 - - Bus 89 1.005 6.07 -0.059 1.06 501.97 -23.99 1.06 501.84 -24.01 Bus 90 0.985 - 0.5931 1.0418 - 47.2 1.042 47.21 Bus 91 0.98 - -0.131 1.0456 - 1.5637 1.046 1.56 Bus 92 0.99 - -0.14 1.05 - 1.05 Bus 93 0.9854 - - 1.04 - - 1.04 - - Bus 94 0.9898 - - 1.0379 - - 1.038 - - Bus 95 0.9803 - - 1.0267 - - 1.027 - - Bus 96 0.9923 - - 1.0354 - - 1.035 - - Bus 97 1.0112 - - 1.0431 - - 1.043 - - Bus 98 1.0235 - - 1.0519 - - 1.052 - - Trang 100 Chương VI Kết - nhận xét Bus V Ps Qg V Ps Qg V Ps Qg # [p.u.] [p.u.] [p.u.] [p.u.] [MW] [Mvar] [p.u.] [MW] [Mvar] Bus 99 1.01 - -0.175 1.054 - -3.17 1.054 -3.17 Bus 100 1.017 2.52 0.9556 1.0584 231.33 45.716 1.058 231.29 45.77 Bus 101 0.9914 - - 1.0428 - - 1.043 - - Bus 102 0.9891 - - 1.0468 - - 1.047 - - 38.25 11.72 22.89 5.16 9.29 - - 29.03 3.26 Bus 103 1.01 0.4 0.7542 1.0514 38.249 Bus 104 0.971 - 0.0239 1.0426 - Bus 105 0.965 - -0.37 1.0403 5.2003 Bus 106 0.9611 - Bus 107 0.952 - Bus 108 0.9662 - - 1.0393 - - 1.039 - - Bus 109 0.967 - - 1.039 - - 1.039 - - Bus 110 0.973 - -0.054 1.0414 7.0559 19.668 1.041 7.03 19.67 Bus 111 0.98 0.36 -0.018 1.0491 35.24 -0.246 1.049 35.24 -0.25 Bus 112 0.975 - 0.4151 1.0341 36.493 10.28 1.034 36.48 10.28 Bus 113 0.993 - 0.0797 1.0554 - -10.08 1.055 -10.23 Bus 114 0.9601 - - 1.0364 - - 1.037 - - Bus 115 0.96 - - 1.036 - - 1.036 - - Bus 116 1.005 - 12.029 1.015 6.41 Bus 117 0.9738 - - 1.0351 - - 1.035 - - Bus 118 0.9494 - - 1.0122 - - 1.012 - - 4319.69 389.87 4319.40 388.26 Tổng 43.7553 Hàm đối tượng [$/h] - 1.0343 0.0112 1.0342 0.3766 1.0185 7.167 29.047 - 129678.207 11.701 1.051 23 1.043 9.2065 1.04 - 1.034 3.2567 1.034 129660.69 Tối ưu phân bố công suất Dự đốn Newton Phương pháp IP hiệu chỉnh thơng thường Số vòng lặp 13 67 3.5437  6.6575  Thời gian tính tốn (s) Nhận Xét: ¾ Kết phân tích chương trình ¾ Khi tối ưu phân bố cơng suất thì: Điện áp nâng lên đáng kể (đều lớn p.u) Tổng công suất nguồn phát thấp trường hợp phân bố công suất thường điều có nghĩa tổn hao thấp ¾ Thời gian tính tốn theo giải thuật dự đốn hiệu chỉnh nhanh phương pháp Newton thông thường Trang 101 Chương VI Kết - nhận xét II KẾT QUẢ KIỂM TRA MẠNG 24 NÚT Phân bố công suất thường IEEE One Area RTS -96 P=4 V=1.05 P =4 V =1.05 Bus 18 |V| = 1.05 p.u