Luận văn thạc sĩ nghiên cứu bù tối ưu cho lưới điện phân phối trung áp xét đến xác suất của phụ tải

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP NGUYỄN SỸ TÙNG NGHIÊN CỨU BÙ TỐI ƯU CHO LƯỚI ĐIỆN PHÂN PHỐI TRUNG ÁP XÉT ĐẾN XÁC SUẤT CỦA PHỤ TẢI LUẬN VĂN THẠC SỸ KỸ THUẬT Chuyên ngành: Kỹ thuật điện THÁI NGUYÊN – 2020 ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP NGUYỄN SỸ TÙNG NGHIÊN CỨU BÙ TỐI ƯU CHO LƯỚI ĐIỆN PHÂN PHỐI TRUNG ÁP XÉT ĐẾN XÁC SUẤT CỦA PHỤ TẢI Chuyên ngành: Kỹ thuật điện Mã số: 852 02 01 LUẬN VĂN THẠC SĨ KỸ THUẬT KHOA CHUYÊN MÔN NGƯỜI HƯỚNG DẪN KHOA HỌC TRƯỞNG KHOA TS VŨ VĂN THẮNG PHÒNG ĐÀO TẠO THÁI NGUYÊN – 2020 i LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng tôi, nghiên cứu kết trình bày luận văn trung thực chưa công bố luận văn trước Tác giả luận văn Nguyễn Sỹ Tùng ii LỜI CẢM ƠN Tôi xin chân thành cảm ơn giúp đỡ, hướng dẫn bảo tận tình thầy giáo TS Vũ Văn Thắng thầy giáo, cô giáo môn Hệ thống điện, Khoa điện, trường Đại học Kỹ thuật Công nghiệp, giúp đỡ chân tình bạn đồng nghiệp, gia đình tạo điều kiện giúp tơi hồn thành luận văn Trong trình thực hiện, thời gian hạn hẹp nên luận văn có thiếu sót Tơi mong muốn nhận ý kiến đóng góp để luận văn hồn thiện thêm kết nghiên cứu thực có ý nghĩa góp phần nâng cao chất lượng điện hệ thống điện Việt Nam Tôi xin trân trọng cảm ơn! Tác giả luận văn Nguyễn Sỹ Tùng iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT vii DANH MỤC CÁC BẢNG viii DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ viii MỞ ĐẦU x I Lý chọn đề tài x II Mục đích nghiên cứu x III Đối tƣợng phạm vi nghiên cứu xi VI Ý nghĩa khoa học thực tiễn xi CHƢƠNG LƢỚI ĐIỆN PHÂN PHỐI VÀ BÙ CÔNG SUẤT PHẢN KHÁNG TRONG LĐPPTA 1.1 Lƣới điện phân phối trung áp 1.1.1 Giới thiệu lƣới điện phân phối trung áp 1.1.2 Đặc điểm lƣới điện phân phối trung áp 1.1.2.1 1.1.2.2 1.1.2.3 1.1.2.4 Phân loại lưới điện trung áp Vai trò lưới điện trung áp Các phần tử lưới điện trung áp Cấu trúc lưới điện trung áp 1.1.3 Hiện trạng lƣới điện trung áp Việt Nam 1.1.3.1 Tình hình phát triển lưới điện trung áp nước ta 1.1.3.2 Tình hình phát triển phụ tải điện 1.2 Đặc tính tải LĐPP 1.2.1 Đồ thị phụ tải ngày 1.2.2 Tính ngẫu nhiên phụ tải điện 11 1.3 Chất lƣợng điện LĐPPTA 11 1.3.1 Điện áp 11 1.3.2 Hệ số công suất 12 1.3.3 Tần số 13 1.3.4 Sóng hài 13 iv 1.3.5 Sự nhấp nháy điện áp 13 1.3.6 Dòng ngắn mạch thời gian loại trừ cố 14 1.4 Tổn thất vấn đề giảm tổn thất LĐPPTA 14 1.4.1 Các nguyên nhân gây tổn thất LĐPP .14 1.4.1.1 Tổn thất kỹ thuật 15 1.4.1.2 Tổn thất phi kỹ thuật 15 1.4.2 Các yếu tố ảnh hƣởng đến tổn thất khả giảm thiểu tổn thất 16 1.4.2.1 Điện áp làm việc trang thiết bị .16 1.4.2.2 Truyền tải CSPK 16 1.4.2.3 Các biện pháp giảm thiểu tổn thất LĐPPTA .18 1.5 Hiệu biện pháp bù CSPK giảm tổn thất LĐPPTA18 1.5.1 Khái niệm CSPK 18 1.5.2 Bù CSPK hệ thống điện 19 1.5.3 Hệ số công suất quan hệ với bù CSPK .20 1.6 Kết luận chƣơng 22 CHƢƠNG PHƢƠNG PHÁP TÍNH TỐN BÙ TRONG LĐPPTA .23 2.1 Đặt vấn đề 23 2.2 Thiết bị bù CSPK 23 2.2.1 Máy phát máy bù đồng 23 2.2.2 Tụ bù tĩnh 24 2.2.2.1 Tụ bù tĩnh cố định 24 2.2.2.2 Tụ bù tĩnh điều chỉnh có cấp 24 2.2.3 Thiết bị bù điều chỉnh vô cấp SVC (Static Var Compensater) .25 2.2.4 Động điện 25 2.2.5 Nhận xét 26 2.3 Phƣơng thức bù LĐPP 26 2.4 Các phƣơng pháp tính tốn bù LĐPP 27 2.4.1 Bù CSPK nâng cao hệ số cos  27 2.4.2 Cực tiểu tổn thất công suất 29 2.4.3 Theo điều kiện chỉnh điện áp 30 2.4.4 Phƣơng pháp bù kinh tế 31 2.4.4.1 Cực đại hóa lợi nhuận 32 2.4.4.2 Cực tiểu chi phí tính tốn Zmin 33 v 2.4.4.3 Cực đại lợi nhuận 2.4.5 Nhận xét 2.5Kết luận chƣơng CHƢƠNG XÂY DỰNG MƠ HÌNH VÀ CHƢƠNG TRÌNH TÍNH TỐN BÙ TRONG LĐPPTA XÉT ĐẾN XÁC SUẤT CỦA PHỤ TẢI 3.1Đặt vấn đề 3.2Mơ hình tốn 3.2.1 Hàm mục tiêu 3.2.2 Các ràng buộc 3.3Cơng cụ tính tốn 3.3.1 Đặt vấn đề 3.3.2 Giới thiệu phần mềm PSS/Adept 3.3.2.1 Chức PSS/Adept 3.3.2.2 Các bước thực 3.3.3 Lập chƣơng trình tính tốn vị trí dung 42 3.3.3.1 Giới thiệu ngơn ngữ lập trình GAMS 3.3.3.2 Cấu trúc chương trình 3.3.3.3 Thuật tốn solver BONMIN chươn 3.4Ví dụ 3.5Kết luận chƣơng CHƢƠNG NGHIÊN CỨU ÁP DỤNG BÙ TỐI ƢU CHO LĐPPTA HUYỆN TIÊN DU 4.1Giới thiệu LĐPPTA huyện Tiên Du 4.1.1 Hiện trạng nguồn cung cấp điện 4.1.2 Hiện trạng LĐPPTA trạm biến áp c 4.1.3 Hiện trạng bù LĐPPTA 4.1.4 Đồ thị phụ tải 4.2Hiện trạng tổn thất thông số chế độ lộ 4.2.1 Sơ đồ thông số lộ 478-E27.1 Tiên D 4.2.2 Kết tính toán 4.2.3 Nhận xét vi 4.3 Tính tốn vị trí dung lƣợng bù tối ƣu lộ 478-E27.1 xét đến xác xuất phụ tải 67 4.3.1 Sơ đồ thông số LĐPPTA lộ 478-E27.1 .67 4.3.2 Kết tính tốn vị trí dung lƣợng bù tối ƣu 67 4.3.3 Đánh giá tổn thất điện chất lƣợng điện áp .68 4.3.4 Nhận xét 73 4.4 Kết luận chƣơng 73 KẾT LUẬN VÀ KIẾN NGHỊ 79 TÀI LIỆU THAM KHẢO 80 PHỤ LỤC 82 vii DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT CSTD Công suất tác dụng CSPK Công suất phản kháng ĐD Đƣờng dây HTĐ Hệ thống điện GAMS Ngơn ngữ lập trình (The General Algebraic Modeling System) MC Máy cắt MBA Máy biến áp LĐPP Lƣới điện phân phối LĐPPTA Lƣới điện phân phối trung áp SCADA Điều khiển giám sát thu thập liệu (Supervisory Control And Data Acquisition) PSS/Adept Phần mềm (Power System Simulator/Avancer Distribution Enginering Productivity tool) TBPĐ Thiết bị phân đoạn viii DANH MỤC CÁC BẢNG Bảng 1.1: Độ biến dạng sóng hài điện áp 13 Bảng 1.2: Giới hạn độ nhấp nháy điện áp 14 Bảng 1.3: Dòng ngắn mạch lớn cho phép thời gian loại trừ cố 14 Bảng 2.1: Giá trị kkt theo phƣơng thức cấp điện 28 Bảng 3.1: Modul thuật toán giải GAMS 44 Bảng 4.1: Hiện trạng nguồn cấp 55 Bảng 4.2: Bảng thông số trạng tải MBA phân phối 56 Bảng 4.3: Bảng thống kê dung lƣợng bù tụ điện 58 Bảng 4.4: Thông số phụ tải 59 Bảng 4.5: Thông số đƣờng dây lộ 478-E27.1 63 Bảng 4.6: Xác suất tải 67 Bảng 4.7: Thông số tụ bù 67 Bảng 4.8: Vị trí dung lƣợng bù tối ƣu 68 Bảng 4.9: Tổn thất điện thời gian tính tốn 68 Bảng 4.10: Công suất cực đại đƣờng dây thời gian tính tốn .72 DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1.1: Sơ đồ lƣới phân phối hình tia Hình 1.2: Sơ đồ lƣới phân phối hình tia có phân đoạn Hình 1.3: Sơ đồ lƣới kín vận hành hở nguồn cung cấp Hình 1.4: Sơ đồ lƣới kín vận hành hở nguồn cung cấp độc lập Hình 1.5: Sơ đồ lƣới điện kiểu đƣờng trục Hình 1.6: Sơ đồ lƣới điện có đƣờng dây dự phịng chung .7 Hình 1.7: Sơ đồ hệ thống phân phối điện Hình 1.8: Đồ thị phụ tải ngày 10 Hình 1.9: Đồ thị phụ tải ngày điển hình theo mùa 10 Hình 1.10: Biến thiên điện áp lƣới điện 12 Hình 1.11: Sơ đồ tham số mạch điện 18 72 tính tốn u cầu đƣợc đảm bảo chƣơng trình kết tính tốn kiểm tra cơng suất truyền tải đoạn đƣờng dây đƣợc trình bày Bảng 4.10 Kết cho thấy, công suất truyền tải tất đƣờng dây đảm bảo yêu cầu, công suất lớn đạt 68.1% công suất giới hạn đƣờng dây đoạn 1.2 công suất truyền tải đoạn 37.38 0.87% Bảng 4.10: Công suất cực đại đƣờng dây thời gian tính tốn Nút i-j S (MVA max 1.2 2.3 3.4 4.5 5.6 6.7 7.8 8.9 8.29 9.10 10.11 11.12 11.30 12.13 12.33 13.14 14.15 15.16 15.36 16.17 17.18 18.19 19.20 19.39 20.21 21.22 22.23 73 23.24 24.25 24.44 25.26 26.27 4.3.4 Nhận xét Kết tính tốn lựa chọn vị trí dung lƣợng bù tối ƣu xét đến thay đổi tải theo mơ hình xác suất lộ 478-E27.1 cho thấy, chƣơng trình tính đƣợc lập trình có khả tính tốn phù hợp cho LĐPPTA thực tế Vị trí dung lƣợng bù tối ƣu đƣợc lựa chọn đồng thời đảm bảo đƣợc độ lệch điện áp nút theo yêu cầu nhƣ đảm bảo vận hành đƣờng dây thiết bị khác không bị tải Tổng công suất bù đƣợc lựa chọn 1800kVAr với tụ bù, gam dung lƣợng 3x300kVAr/bộ đƣợc đặt vị trí xa nguồn giảm đƣợc tổn thất điện Trong giai đoạn tính tốn năm, tổn thất điện giảm đƣợc 0.24% so với bù trạng 1.24 so với phƣơng án không sử dụng tụ bù Hơn nữa, phƣơng pháp đảm bảo đƣợc độ lệch điện áp cho phép nhƣ giới hạn công suất truyền tải đƣờng dây Do đó, ln đảm bảo đƣợc yêu cầu chất lƣợng điện áp trạng thái vận hành nhƣ đƣờng dây khơng bị q tải suốt giai đoạn tính tốn 4.4 Kết luận chƣơng Từ kết tính toán áp dụng cho LĐPPTA huyện Tiên Du – lộ 478-E27.1 cho thấy hiệu thiết bị bù nhƣ phƣơng pháp chƣơng trình tính tốn lựa chọn vị trí dung lƣợng bù tối ƣu xét đến thay đổi tải theo mơ hình xác suất đề xuất Vị trí, dung lƣợng thời gian lắp đặt tụ điện đƣợc lựa chọn với tổn thất điện suốt thời gian tính tốn đảm bảo cực tiểu Chế độ làm việc thiết bị đƣợc đảm bảo không bị tải đồng thời độ lệch điện áp tất vị trí chế độ vận hành suốt giai đoạn tính tốn đƣợc đảm bảo theo u cầu Phƣơng pháp chƣơng trình tính đề xuất khắc phục đƣợc hạn chế phƣơng pháp tính tốn vị trí, dung lƣợng bù sử dụng cho LĐPPTA lựa chọn đƣợc vị trí xa nguồn dẫn đến đảm bảo giảm tổn thất nhỏ phƣơng pháp sử dụng 74 KẾT LUẬN VÀ KIẾN NGHỊ * Kết luận Nghiên cứu tổng hợp, nghiên cứu số vấn đề tổn thất, phƣơng pháp giảm tổn thất LĐPP thiết bị bù Lựa chọn phƣơng pháp tính tốn thơng số bù xét đến thay đổi phụ tải theo mơ hình xác suất nhằm nâng cao tính xác kết tính tốn đồng thời đảm bảo độ lệch điện áp hầu hết chế độ vận hành hệ thống Giới thiệu phân tích mơ hình tốn, chƣơng trình tính lựa chọn vị trí dung lƣợng bù tối ƣu xét đến ảnh hƣởng thay đổi phụ tải theo mơ hình xác suất phát triển tải giai đoạn tính tốn Chƣơng trình đƣợc sử dụng để tính tốn cho ví dụ đơn giản sau kiểm chứng kết phần mềm PSS/Adept cho thấy phƣơng thức bù chƣơng trình tính tốn phù hợp với LĐPPTA, đảm bảo yêu cầu hệ thống, công suất vị trí bù đƣợc lựa chọn tối ƣu, chƣơng trình tính tốn đơn giản, dễ áp dụng Kết tính tốn áp dụng cho LĐPPTA huyện Tiên Du – lộ 478-E27.1 cho thấy hiệu thiết bị bù nhƣ phƣơng pháp chƣơng trình tính tốn lựa chọn vị trí dung lƣợng bù tối ƣu theo mơ hình xác suất tải Vị trí, dung lƣợng thời gian lắp đặt tụ điện đƣợc lựa chọn với tổn thất điện suốt thời gian tính tốn đảm bảo cực tiểu Chế độ làm việc thiết bị đƣợc đảm bảo không bị tải đồng thời độ lệch điện áp tất vị trí vận hành, mùa tất năm tính toán đƣợc đảm bảo theo yêu cầu * Hướng nghiên cứu Nghiên cứu phƣơng pháp giải mơ hình tốn với số lƣợng trạng thái tính tốn lớn nhằm nâng cao độ xác kết tính tốn Nghiên cứu xét đến chi phí tụ điện nhƣ giá điện nhằm nâng cao hiệu kinh tế phƣơng án bù 75  TÀI LIỆU THAM KHẢO Tiếng Việt [1] Bộ Công Thƣơng (2015), Quy định hệ thống điện phân phối, Thông tƣ 39/2015/TT-BCT ngày 18 tháng 11 năm 2015, Hà Nội [2] Bộ Công Thƣơng (2015), Quy chuẩn ký thuật quốc gia kỹ thuật điện, QCVN 621: 2015/BCT ngày 20 tháng 09 năm 2015, Hà Nội [3] Đỗ Xn Khơi (1998), Tính tốn phân tích hệ thống điện, Nhà xuất Khoa học Kỹ thuật, Hà Nội [4] Nguyễn Xuân Phú, Nguyễn Công Hiền, Nguyễn Bội Khuê (2003), Cung cấp điện, NXB Khoa học kỹ thuật, Hà Nội [5] Trần Quang Khánh (2012), Cung cấp điện theo tiêu chuẩn IEC, Nhà xuất Khoa học Kỹ thuật, Hà Nội [6] Viện Năng lƣợng (2011), Qui hoạch phát triển điện lực quốc gia giai đoạn 2011-2020 có xét đến năm 2030, Hà Nội [7] Viện Năng lƣợng (2016), Điều chỉnh qui hoạch phát triển điện lực quốc gia giai đoạn 2011-2020 có xét đến năm 2030, Hà Nội [8] Vũ Văn Thắng, Bạch Quốc Khánh (2017), Vị trí cơng suất tối ưu tụ điện qui hoạch cải tạo hệ thống phân phối, Tạp chí KHCN – Đại học Đà nẵng, Số 3(112).2017  [9] Báo cáo kỹ thuật, Điện lực Tiên Du, Bắc Ninh 2019 Tiếng Anh [10] Anthony J Pansini (2007), Electrical distribution engineering, The Fairmont Press, Inc [11] Pierre Bonami and Jon Lee (2007), BONMIN Users’ Manual, GAMS Development Corporation, Washington, DC, USA [12] Richard E Rosenthal (2010), GAMS - A User's Guide, GAMS Development Corporation, Washington, DC, USA [13] PSS/Adept™ 5, Users Manual, Shaw Power Technologies, Inc 2004 [14] M.Jayalakshmi, K.Balasubramanian, Simple Capacitors to Supercapacitors-An Overview, International Journal of Electrochemical Science, Vol.3, 2008 [15] Trần Vinh Tịnh, T.V Chƣơng, Bù tối ƣu công suất phản kháng LPP, Tạp chí KH&CN Đại học Đà Nẵng, số 2, 2008 M.Dixit, P.Kundu, H R.Jariwala, Optimal Allocation and Sizing of Shunt Capacitor in Distribution System for Power Loss Minimization, 2016SCEECS, India, 2016 [16] 76 A.A.A.El-Ela, A.M.Kinawy, M.T.Mouwafi, R.A.El-Sehiemy, Optimal sitting and sizing of capacitors for voltage enhancement of distribution systems, 2015UPEC, UK, 2015 [17] A.A.Eajal, M.E.El-Hawary, Optimal capacitor placement and sizing in distorted radial distribution systems part III: Numerical results, ICHQP2010, Italy, 2010 [18] [19] K.R Devabalaji, A.M.Imranb, T.Yuvaraj, K.Ravi, Power Loss Minimization in Radial Distribution System, Energy Procedia 79 (2015), 917-923 N.Rugthaicharoencheep, S.Nedphograw, W.Wanaratwijit, Distribution system operation for power loss minimization and improved voltage profile with distributed generation and capacitor placements, 2011DRPT, China, 2011 [20] [21] Y M Atwa, E F El-Saadany, M M A Salama, and R Seethapathy, Optimal Renewable Resources Mix for Distribution systems Energy Loss Minimization, IEEE Tran on Power Sytems, Vol.25, No.1, 2010 [22] Solutions for power factor correction at medium voltage, CIRCUTOR, S.A 2013 77 PHỤ LỤC PL1 Thơng số ví dụ * Thông số tải Nút tải 10 11 12 13 14 15 16 Tổng * Thông số đường dây Sij* ,t , MVA Nút ij 1.2 2.3 3.4 4.5 5.6 6.7 7.8 8.9 9.10 10.11 11.12 12.13 13.14 14.15 15.16 16.17 7 78 PL2 Chƣơng trình tính ngơn ngữ lập trình GAMS OPTION LIMROW = 30; OPTION LIMCOL = 30; Option ITERLIM = 50000; Option Solprint = Off; Option Sysout = Off; Option Reslim = 5000000; OPTION OPTCA = 5.35; OPTION OPTCR = 5.35; Option decimals = 2; ******* KHAI BAO MANG DU LIEU VA DU LIEU DAU VAO********** set i So nut /1*40/; alias (i,j); set t Thoi gian tinh toan /1*3/; set h So gio /1*24/; set s mua nam /1*4/; set n so bien tuong ung cong suat roi rac cua Tu /1*12/; Set Gen(i) Nut nguon /1/ Load(i) Nut tai /2*40/ cp(i) Nut dau tu Tu /1*40/ Head1 Bang du lieu duong day /Rf, Xf, Chf, sf/ Head2 Bang du lieu nguon /Pmin, Pmax, Qmin, Qmax/ Head3 Bang du lieu tai /PD0, QD0/ Head4 Bang cong suat gioi han TBA /Ss0/ Head5 Bang cong suat gioi han duong day /Sf0/ Head6 Bang du lieu CP /Qmin, Qmax/ Head10 mua nam /1*4/ ; Scalar phi/3.141592654/; Scalar SBase Cong suat co so MVA /1/; Scalar VBase Dien ap co so kV /35/; Scalar GFD He so phat trien cua phu tai /0.2/; TABLE Generat(Gen,Head2) Du lieu nguon Pmin Pmax Qmin Qmax ; Parameter PmxG(gen,t), PmnG(gen,t), QmxG(gen,t), QmnG(gen,t); PmxG(gen,t) = Generat(gen,"Pmax")/(SBase); 79 PmnG(gen,t) = Generat(gen,"Pmin")/(SBase); QmxG(gen,t) = Generat(gen,"Qmax")/(SBase); QmnG(gen,t) = Generat(gen,"Qmin")/(SBase); TABLE Generat1(cp,Head6) Du lieu CP ; Parameter QmxG1(cp,t), QmnG1(cp,t); QmxG1(cp,t) = 10*Generat1(cp,"Qmax")/(SBase); QmnG1(cp,t) = 1*Generat1(cp,"Qmin")/(SBase); TABLE Demand1(s,h,Head10) He so tai va xac suat ; Parameter kt(s,h); kt(s,h) = Demand1(s,h,"s1"); TABLE Demand(i, Head3) Cong suat phu tai PD0 QD0 ; Parameter PD(i,t), QD(i,t); PD(i,t)= 1.2*Demand(i,"PD0")*(1+GFD*(ord(t)-1))/(1000*SBase); QD(i,t)= 1.2*Demand(i,"QD0")*(1+GFD*(ord(t)-1))/(1000*SBase); Table LineData(i,j,head1) Rf Xf Chf sf ; LineData(j,i,head1) = LineData(i,j,head1) ; Parameter Se(i,j,t); Se(i,j,t) = LineData(i,j,"Sf"); Se(i,j,t) = Se(j,i,t); Parameter Re(i,j,t,head1), Xe(i,j,t,head1), Che(i,j,t,head1); Re(i,j,t,"Rf") = Linedata(i,j,"Rf")*(SBase/(VBase*VBase)) ; Xe(i,j,t,"Xf") = Linedata(i,j,"Xf")*(SBase/(VBase*VBase)) ; Che(i,j,t,"Chf")=Linedata(i,j,"Chf"); Re(j,i,t,"Rf")$ (Re(i,j,t,"Rf") gt 0) = Re(i,j,t,"Rf") ; Xe(j,i,t,"Xf")$ (Xe(i,j,t,"Xf") gt 0) = Xe(i,j,t,"Xf") ; Che(j,i,t,"Chf")$ (Che(i,j,t,"Chf") gt 0) = Che(i,j,t,"Chf") ; Parameter Z(i,j,t), GG(i,j,t), BB(i,j,t); Z(i,j,t) = (Re(i,j,t,"Rf")*Re(i,j,t,"Rf")) + (Xe(i,j,t,"Xf")*Xe(i,j,t,"Xf")) ; GG(i,j,t)$(z(i,j,t) ne 0.00) = Re(i,j,t,"Rf")/Z(i,j,t) ; BB(i,j,t)$(z(i,j,t) ne 0.00) = -Xe(i,j,t,"Xf")/Z(i,j,t); Parameter YCL(i,t); YCL(i,t) = sum(j, Che(i,j,t,"Chf")); 80 Parameter G(i,j,t) , B(i,j,t) ; B(i,i,t) = sum(j,BB(i,j,t)); G(i,i,t) = sum(j,GG(i,j,t)); G(i,j,t)$(ord(i) ne ord(j)) = -GG(i,j,t); B(i,j,t)$(ord(i) ne ord(j)) = -BB(i,j,t); Parameter Y(i,j,t); Y(i,j,t) = sqrt(G(i,j,t)*G(i,j,t) + B(i,j,t)*B(i,j,t)); Parameter ZI(i,j,t); ZI(i,j,t)$(G(i,j,t) ne 0.00) = abs(B(i,j,t))/abs(G(i,j,t)) ; Parameter theta(i,j,t); theta(i,j,t) = arctan(ZI(i,j,t)); theta(i,j,t)$((b(i,j,t) eq 0) and (g(i,j,t) gt 0)) = 0.0 ; theta(i,j,t)$ ((b(i,j,t) eq 0) and (g(i,j,t) lt 0)) = phi ; theta(i,j,t)$((b(i,j,t) gt 0) and (g(i,j,t) gt 0)) = theta(i,j,t) ; theta(i,j,t)$((b(i,j,t) lt 0) and (g(i,j,t) gt 0)) = 2*phi - theta(i,j,t) ; theta(i,j,t)$((b(i,j,t) gt 0) and (g(i,j,t) lt 0)) = phi - theta(i,j,t); theta(i,j,t)$((b(i,j,t) lt 0) and (g(i,j,t) lt 0)) = phi + theta(i,j,t); theta(i,j,t)$((b(i,j,t) gt 0) and (g(i,j,t) eq 0)) = 0.5*phi; theta(i,j,t)$((b(i,j,t) lt 0) and (g(i,j,t) eq 0)) = -0.5*phi; theta(i,j,t)$((b(i,j,t) eq 0) and (g(i,j,t) eq 0)) = 0.0 ; Parameter G(i,j,t); G(i,j,t) = -Y(i,j,t)*cos(Theta(i,j,t)); *********** KHAI BAO CAC BIEN VA XAP XI DAU ***************** VARIABLES V(i,t,s,h) Dien ap nut Delta(i,t,s,h) Goc pha dien ap nut P(i,t,s,h) Cong suat TD nhan tu nguon Q(i,t,s,h) Cong suat PK nhan tu nguon QG(i,t) Cong suat PK Tu QG1(i,t) Cong suat PK Tu1 Loss U(n,i,t) Tong ton that Bien nhi phan ; Positive variable P, Q, QG; Binary variable U; Parameter VLevel(Gen) /1 V.l(i,t,s,h) V.Fx(Gen,t,s,h) = VLevel(Gen); Delta.l(i,t,s,h) Delta.fx("1",t,s,h) = 0.0; 81 ************** HAM MUC TIEU VA RANG BUOC ******************* Equations CostEq Equn1(i,t,s,h) Equn2(i,t,s,h) Equn3(i,t) Equn4(cp) Equn5(i,j,t) Equn6(i,t) ; ************** HAM MUC TIEU ************** CostEq Loss =e= (365/2)*0.5*Sum((i,j,t,s,h),G(i,j,t)*(V(i,t,s,h)*V(i,t,s,h) + V(j,t,s,h)*V(j,t,s,h) 2*V(i,t,s,h)*V(j,t,s,h)*cos(Delta(j,t,s,h)-Delta(i,t,s,h))) + QG(i,t)*0.0002)*Sbase*1000; ************** CAN BANG CONG SUAT NUT ************** Equn1(i,t,s,h)$(ord(i) gt 1) -QG(i,t)*0.0002- PD(i,t)*kt(s,h) =e= Sum(j, Y(i,j,t)*V(i,t,s,h)*V(j,t,s,h)*Cos(theta(i,j,t)+Delta(j,t,s,h) Delta(i,t,s,h))); Equn2(i,t,s,h)$(ord(i) gt 1) (QG(i,t-4)+QG(i,t-3)+QG(i,t-2)+QG(i,t-1)+QG(i,t))QD(i,t)*kt(s,h) =e= -Sum(j, Y(i,j,t)*V(i,t,s,h)*V(j,t,s,h)*Sin(theta(i,j,t)+Delta(j,t,s,h) Delta(i,t,s,h))); ************** RANG BUOC NANG CAP THIET BI ************** Equn2(i,t)$(ord(t) gt 0) QG(i,t) =e= (0.05*U("12",i,t)+0.1*U("11",i,t)+0.15*U("10",i,t)+0.21*U("9",i,t)+0 24*U("8",i,t)+0.3*U("7",i,t) +0.33*U("6",i,t)+0.39*U("5",i,t)+0.42*U("4",i,t)+0.45*U("3",i,t) +0.5*U("2",i,t)+0.55*U("1",i,t) ); Equn3(cp) sum((n,t),U(n,cp,t)) =l= 1; Equn4 sum((n,cp,t),U(n,cp,t)) =l= 5; Equn5(i,j,t)$(ord(i) lt ord(j)) (Smax((s,h),abs(Sbase*V(i,t,s,h)*Y(i,j,t)*(V(j,t,s,h)V(i,t,s,h)))))=l= Se(i,j,t); Equn6(i,t)$(ord(t) gt 0) QG1(i,t) =e= QG1(i,t-1)+QG(i,t); ********************** GIOI HAN CAC BIEN ********************** QG.Up(cp,t) = 1*QmxG1(cp,t); QG.Lo(cp,t) = 1*QmnG1(cp,t); 82 ***************** RANG BUOC DIEN AP NUT ********************** V.Up(load,t,s,h) V.Lo(load,t,s,h) ******************* GOI SLOVER TINH TOAN ********************** Option MINLP=BONMIN; Model OPF /all/; Solve OPF using MINLP Minimizing Loss; ************ TINH TOAN THONG SO VA HIEN THI ****************** Display Loss.l, V.l; ... CÔNG NGHIỆP NGUYỄN SỸ TÙNG NGHIÊN CỨU BÙ TỐI ƯU CHO LƯỚI ĐIỆN PHÂN PHỐI TRUNG ÁP XÉT ĐẾN XÁC SUẤT CỦA PHỤ TẢI Chuyên ngành: Kỹ thuật điện Mã số: 852 02 01 LUẬN VĂN THẠC SĨ KỸ THUẬT KHOA CHUYÊN MÔN... Section 1.1 Lƣới điện phân phối trung áp 1.1.1 Giới thiệu lưới điện phân phối trung áp Lƣới điện phân phối trung áp (LĐPPTA) phần hệ thống điện, làm nhiệm vụ phân phối điện từ trạm trung gian, trạm... Lƣới điện phân phối trung áp 1.1.1 Giới thiệu lƣới điện phân phối trung áp 1.1.2 Đặc điểm lƣới điện phân phối trung áp 1.1.2.1 1.1.2.2 1.1.2.3 1.1.2.4 Phân loại lưới điện trung áp