Nghiên cứu giải pháp giảm tổn thất điện năng trên lưới điện phân phối huyện sốp cộp

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI KỸ THUẬT CÔNG NGHIỆP NGUYỄN DUY TUỆ NGHIÊN CỨU GIẢI PHÁP GIẢM TỔN THẤT ĐIỆN NĂNG TRÊN LƯỚI ĐIỆN PHÂN PHỐI HUYỆN SỐP CỘP LUẬN VĂN THẠC SĨ KHOA HỌC KỸ THUẬT ĐIỆN THÁI NGUYÊN - NĂM 2020 ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI KỸ THUẬT CÔNG NGHIỆP NGUYỄN DUY TUỆ NGHIÊN CỨU GIẢI PHÁP GIẢM TỔN THẤT ĐIỆN NĂNG TRÊN LƯỚI ĐIỆN PHÂN PHỐI HUYỆN SỐP CỘP Chuyên ngành: Kỹ thuật điện Mã ngành: 8.52.02.01 LUẬN VĂN THẠC SĨ KHOA HỌC KỸ THUẬT ĐIỆN NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Vũ Văn Thắng THÁI NGUYÊN - NĂM 2020 CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc BẢN XÁC NHẬN CHỈNH SỬA LUẬN VĂN THẠC SĨ Họ tên tác giả luận văn: Nguyễn Duy Tuệ Đề tài luận văn: Nghiên cứu giải pháp giảm tổn thất điện lưới điện phân phối huyện Sốp Cộp Chuyên ngành: Kỹ thuật điện Mã số: : 8.52.02.01 Tác giả, Cán hướng dẫn khoa học Hội đồng chấm luận văn xác nhận tác giả sửa chữa, bổ sung luận văn theo biên họp Hội đồng ngày 10/10/2020 với nội dung sau: - Sửa sai sót thuật ngữ, lỗi tả, format, in ấn - Đã chỉnh sửa lỗi câu chữ, thuật ngữ luận văn cho phù hợp với chuyên ngành theo theo góp ý Hội đồng Thái Nguyên, ngày 26 tháng 10 năm 2020 Cán hướng dẫn Tác giả luận văn TS Vũ Văn Thắng Nguyễn Duy Tuệ CHỦ TỊCH HỘI ĐỒNG PGS.TS Nguyễn Như Hiển 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 đƣợc trình bày luận văn trung thực chƣa đƣợc công bố luận văn trƣớc Tác giả luận văn Nguyễn Duy Tuệ 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 đƣợc ý kiến đóng góp để luận văn đƣợc 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 Duy Tuệ 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Ị ix MỞ ĐẦU x I Lý chọn đề tài x II Mục đích nghiên cứu xi III Đối tƣợng phạm vi nghiên cứu xi VI Ý nghĩa khoa học thực tiễn xi 1.1 Cấu trúc đặc điểm lƣới điện phân phối trung áp 1.1.1 Định nghĩa lƣới điện phân phối trung áp 1.1.2 Đặc điểm lƣới điện 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 10 1.2 Chất lƣợng điện LĐPPTA 11 1.2.1 Điện áp 11 1.2.2 Hệ số công suất 12 1.2.3 Tần số 12 1.2.4 Sóng hài 13 1.2.5 Sự nhấp nháy điện áp 13 1.2.6 Dòng ngắn mạch thời gian loại trừ cố 14 1.3 Đặc tính tải LĐPP trung áp 15 1.3.1 Phụ tải đồ thị phụ tải 15 1.3.2 Tính ngẫu nhiên phụ tải điện 16 iv 1.4 Kết luận chƣơng 17 CHƢƠNG PHƢƠNG PHÁP GIẢM TỔN THẤT ĐIỆN NĂNG TRONG LƢỚI ĐIỆN PHÂN PHỐI TRUNG ÁP 18 2.1 Các nguyên nhân gây tổn thất LĐPP 18 2.1.1 Tổn thất kỹ thuật 18 2.1.2 Tổn thất phi kỹ thuật 19 2.2 Các yếu tố ảnh hƣởng đến tổn thất khả giảm thiểu tổn thất 20 2.2.1 Điện áp làm việc trang thiết bị 20 2.2.2 Truyền tải CSPK 20 2.2.3 Các biện pháp giảm thiểu tổn thất LĐPPTA 22 2.3 Giải pháp quản lý nhu cầu điện 23 2.3.1 Khái niệm quản lý nhu cầu điện 23 2.3.2 Chiến lƣợc DSM hệ thống điện 24 2.4 Phƣơng pháp bù CSPK hiệu biện pháp bù CSPK giảm tổn thất LĐPPTA 27 2.4.1 Khái niệm CSPK 27 2.4.2 Bù CSPK hệ thống điện 29 2.4.3 Hệ số công suất quan hệ với bù CSPK 30 2.4.4 Nhận xét 32 2.5 Kết luận chƣơng 32 CHƢƠNG PHƢƠNG PHÁP TÍNH TỐN BÙ TRONG LƢỚI ĐIỆN PHÂN PHỐI TRUNG ÁP XÉT ĐẾN XÁC SUẤT CỦA TẢI 33 3.1 Đặt vấn đề 33 3.2 Thiết bị bù CSPK 33 3.2.1 Máy phát máy bù đồng 34 3.2.2 Tụ bù tĩnh 34 3.2.2.1 Tụ bù tĩnh cố định 34 3.2.2.2 Tụ bù tĩnh điều chỉnh có cấp 35 3.2.3 Thiết bị bù điều chỉnh vô cấp SVC (Static Var Compensater) 35 3.2.4 Động điện 36 3.3 Phƣơng thức bù LĐPP 37 3.4 Các phƣơng pháp tính tốn bù LĐPP 38 v 3.4.1 Bù CSPK nâng cao hệ số cos  38 3.4.2 Cực tiểu tổn thất công suất 41 3.4.3 Theo điều kiện chỉnh điện áp 42 3.4.4 Phƣơng pháp bù kinh tế 43 3.4.4.1 Cực đại hóa lợi nhuận 44 3.4.4.2 Cực tiểu chi phí tính tốn Zmin 45 3.4.4.3 Cực đại lợi nhuận 46 3.4.5 Nhận xét 47 3.5 Xây dựng mơ hình tốn xác định vị trí dung lƣợng bù tối ƣu tụ điện LĐPPTA xét đến tải ngẫu nhiên 47 3.5.1 Hàm mục tiêu 47 3.5.2 Các ràng buộc 48 3.6 Cơng cụ tính tốn 50 3.6.1 Đặt vấn đề 50 3.6.2 Giới thiệu phần mềm PSS/Adept 50 3.6.2.1 Chức PSS/Adept 50 3.6.2.2 Các bước thực 51 3.6.3 Lập chƣơng trình tính tốn vị trí dung lƣợng bù tối ƣu ngơn ngữ lập trình GAMS 52 3.6.3.1 Giới thiệu ngơn ngữ lập trình GAMS 52 3.6.3.2 Thuật tốn solver BONMIN chương trình GAMS 54 3.7 Ví dụ minh họa 54 3.8 Kết luận chƣơng 57 CHƢƠNG NGHIÊN CỨU ÁP DỤNG BÙ TỐI ƢU CHO LƢỚI ĐIỆN PHÂN PHỐI TRUNG ÁP HUYỆN SỐP CỘP 55 4.1 Hiện trạng LĐPPTA huyện Huyện Sốp Cộp 55 4.1.1 Hiện trạng nguồn cung cấp điện 55 4.1.2 Hiện trạng LĐPPTA trạm biến áp 55 4.1.3 Hiện trạng bù LĐPPTA 59 4.1.4 Sơ đồ thông số lƣới điện 59 4.2 Hiện trạng tổn thất thông số chế độ LĐTA Sốp Cộp 67 4.2.1 Kết tính tốn tổn thất với thơng số trạng 67 4.2.2 Nhận xét 69 vi 4.3 Tính tốn vị trí dung lƣợng bù tối ƣu cho LĐPPTA Sốp Cộp xét đến tính ngẫu nhiên phụ tải 69 4.3.1 Sơ đồ thông số LĐPPTA Huyện Sốp Cộp 69 4.3.2 Thông số tải tụ bù 70 4.3.3 Kết tính tốn, lựa chọn vị trí dung lƣợng bù tối ƣu 71 4.3.4 Đánh giá tổn thất điện độ lệch điện áp 74 4.3.5 Nhận xét 78 4.4 Kết luận chƣơng 78 KẾT LUẬN VÀ KIẾN NGHỊ 80 TÀI LIỆU THAM KHẢO 82 PHỤ LỤC 84 vii DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT CSTD CSPK DSM management) ĐD HTĐ GAMS Công suất tác dụng Công suất phản kháng Quản lý nhu cầu sử dụng điện (Demand Side Đƣờng dây Hệ thống điện 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 78 38.39 AC50 12.73 2.29 24.89 AC50 12.73 0.08 39.40 AC50 12.73 2.05 89.90 AC50 12.73 0.04 40.41 AC50 12.73 1.94 25.91 AC70 16.06 0.51 41.42 AC50 12.73 1.8 91.92 AC70 16.06 0.45 42.43 AC50 12.73 1.72 92.93 AC70 16.06 0.35 43.44 AC50 12.73 1.54 93.94 AC70 16.06 0.31 44.45 AC50 12.73 1.38 94.95 AC70 16.06 0.26 45.46 AC50 12.73 0.96 95.96 AC70 16.06 0.15 46.47 AC50 12.73 0.75 96.97 AC70 16.06 0.03 47.48 AC50 12.73 0.68 91.98 AC70 16.06 0.08 48.49 AC50 12.73 0.64 98.99 AC70 16.06 0.07 4.3.5 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ĐPPTA Sốp Cộp 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 đồng thời lựa chọn, vị trí đƣợc lựa chọn nút xa nguồn đảm bảo đƣợc độ lệch điện áp nút theo yêu cầu nhƣ đƣờng dây thiết bị khác không bị tải Mặc dù công suất bù tăng so với trạng nhƣng tổn thất công suất tổn thất điện suốt giai đoạn tính toán giảm, điện áp tất nút đƣợc đảm bảo nâng cao so với trạng nâng cao chất lƣợng điện áp LĐPPTA Sốp Cộp 4.4 Kết luận chƣơng Kết tính tốn áp dụng cho LĐPPTA Sốp Cộp – lộ 371-E17.30 cho thấy, 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 phụ tải theo mơ hình xác suất đề xuất có hiệu cao tính tốn nhƣ thơng số hệ thống 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 đảm bảo cực tiểu suốt thời gian tính tốn Chế độ 78 79 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 tốn tƣơng ứng với tất trạng thái xuất phụ tải đƣợc đảm bảo theo yêu cầu Phƣơng pháp thƣờng đƣợc sử dụng với vị trí đƣợc lựa chọn cảm tính, cơng suất xác định theo hệ số công suất cos, thông số chế độ cần đƣợc tính tốn kiểm tra nên khó lựa chọn đƣợc vị trí bù tối ƣu thời gian tính tốn dài 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 lựa chọn đƣợc vị trí tối ƣu đồng thời đảm bảo cực tiểu tổn thất 79 80 KẾT LUẬN VÀ KIẾN NGHỊ * Kết luận Một số vấn đề cấu trúc, đặc điểm LĐPP nhƣ tốn tính tốn tổn thất, phƣơng pháp giảm tổn thất LĐPP thiết bị bù đƣợc giới thiệu Từ đó, phân tích lựa chọn phƣơng pháp giảm tổn thất công suất tổn thất điện LĐPP nói chung LĐPPTA Sốp Cộp nói riêng 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 đƣợc phân tích lựa chọn 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 đƣợc sử dụng Từ đó, lập chƣơng trình tính tốn ngơng ngữ lập trình GAMS, tính tốn 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 Sốp Cộp – lộ 371E17.30 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 đồng thời, thời gian tính tốn nhỏ Hơn nữa, tổn thất điện suốt thời gian tính tốn ln đảm bảo cực tiểu Thông số chế độ hệ thống đƣợc đảm bảo với công suất truyền tải thiết bị đƣợc đảm bảo giới hạn cho phép, độ lệch điện áp tất vị trí trạng thái vận hành tất năm tính tốn đƣợc đảm bảo theo yêu cầu Hơn nữa, thông số tụ bù đƣợc lựa chọn theo dung lƣợng bù định mức, rời rạc theo thiết bị bù thực tiễn kết tính tốn đảm bảo tối ƣu đồng thời dễ dàng áp dụng thực tế 80 81 * Hướng nghiên cứu - Nghiên cứu xét đến chi phí tụ điện nhằm nâng cao hiệu kinh tế phƣơng án bù - Nghiên cứu ứng dụng thiết bị bù có điều khiển vơ cấp (SVC) nhằm nâng cao chất lƣợng điện áp hệ thống 81 82 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] Bộ Công Thƣơng (2018), Về việc phê duyệt qui hoạch phát triển điện lực tỉnh Sơn La giai đoạn 2019-2025 có xét đến năm 2035 – Qui hoạch phát triển hệ thống điện 110kV, Hà Nội [4] Báo cáo kỹ thuật, Điện Lực Sông Mã, Sơn La 2019 [5] Đỗ Xuân 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 [6] 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 [7] 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 [8] 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 [9] 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 [10] 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  Tiếng Anh [11] Anthony J Pansini (2007), Electrical distribution engineering, The Fairmont Press, Inc [12] Pierre Bonami and Jon Lee (2007), BONMIN Users’ Manual, GAMS Development Corporation, Washington, DC, USA [13] Richard E Rosenthal (2010), GAMS - A User's Guide, GAMS Development Corporation, Washington, DC, USA [14] PSS/Adept™ 5, Users Manual, Shaw Power Technologies, Inc 2004 [15] M.Jayalakshmi, K.Balasubramanian, Simple Capacitors Supercapacitors-An Overview, International Journal Electrochemical Science, Vol.3, 2008 to of 82 83 [16] 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 [17] M.Dixit, P.Kundu, H R.Jariwala, Optimal Allocation and Sizing of Shunt Capacitor in Distribution System for Power Loss Minimization, 2016SCEECS, India, 2016 [18] 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 [19] A.A.Eajal, M.E.El-Hawary, Optimal capacitor placement and sizing in distorted radial distribution systems part III: Numerical results, ICHQP2010, Italy, 2010 [20] K.R Devabalaji, A.M.Imranb, T.Yuvaraj, K.Ravi, Power Loss Minimization in Radial Distribution System, Energy Procedia 79 (2015), 917-923 [21] 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 [22] 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 [23] Solutions for power factor correction at medium voltage, CIRCUTOR, S.A 2013 83 84 PHỤ LỤC PL1 Các bƣớc lập chƣơng trình tính PSS/Adept * Bƣớc 1: Thiết lập thông số lƣới điện Trong bƣớc này, thực khai báo thông số lƣới điện cần tính tốn để mơ PSS/ADEPT gồm nội dung: - Xác định thƣ viện dây dẫn nhƣ Hình P1 Hình P1: Giao diện xác định thƣ viện dây dẫn - Xác định thơng số thuộc tính lƣới điện: Bƣớc này, nhằm khai báo cho phần mềm PSS/Adept thiết lập từ đầu thuộc tính lƣới điện nhƣ: Điện áp qui ƣớc điện áp pha hay điện áp dây trị số, tần số, công suất biểu kiến bản… 84 85 Hình P2: Giao diện xác định thuộc tính lƣới điện * Bƣớc: Thiết lập sơ đồ Thiết lập sơ đồ lƣới điện giao diện PSS/Adept, thiết lập nhật số liệu đầu vào cho phần tử nhƣ nút, cái, nguồn điện, đƣờng dây, trạm biến áp, tải, cầu dao… Hình P3: Giao diện thiết lập thơng số phần tử lƣới điện * Bƣớc: Tính tốn Có phân hệ tính tốn phần mềm PSS/Adept 5.0 hộp thoại tùy chọn chƣơng trình tính tốn nhƣ Hình P4 Trong nghiên cứu sử dụng phân hệ tính tốn trào lƣu cơng suất hay luồng tải (Load Flow) Hình P4: Giao diện hộp tùy chọn chƣơng trình tính tốn * Bƣớc 4: Hiển thị kết Kết tính tốn chƣơng trình đƣợc phân tích vị trí nhƣ sau: - Xem hiển thị kết phân tích sơ đồ nhƣ Hình P5 85 86 - Xem kết tính tốn số progress view nhƣ Hình P6 - Xem kết tính tốn chi tiết từ phần report phần mềm PSS/Adept nhƣ Hình P7 Hình P6: Hiển thị kết tính tốn sơ đồ Hình P7: Hiển thị kết tính tốn số progress view Hình P9: Hiển thị kết tính tốn cửa sổ report 86 87 PL2 Thơng số ví dụ * Thông số tải Nút tải Pi.0, kW 240 348 384 192 432 360 360 192 10 264 11 174 12 192 13 552 14 264 15 672 16 312 Tổng Qi.0, kVAr Nút tải 192 17 288 18 300 19 156 20 372 21 360 22 360 23 144 24 192 25 132 26 162 27 522 28 216 29 552 30 240 31 12,810 Pi.0, kW Qi.0,kVAr 432 384 588 528 228 168 348 264 228 168 468 408 468 420 504 420 264 240 792 630 672 630 432 372 504 444 360 240 660 564 10,040 * Thông số đường dây 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 Sij*,t , MVA Rij,  26 26 26 26 26 10 10 10 10 10 10 8 8 0.15 0.31 0.08 0.23 0.69 1.89 0.46 0.63 0.50 0.84 1.26 1.07 0.83 1.12 1.84 0.95 Xij,  0.38 0.75 0.19 0.57 1.70 1.84 0.45 0.61 0.49 0.82 1.22 0.75 0.59 0.79 1.30 0.67 Nút ij 17.18 2.19 19.20 20.21 21.22 3.23 23.24 24.25 6.26 26.27 27.28 28.29 29.30 30.31 31.32 32.33 Sij*,t , MVA Rij,  8 8 8 8 10 10 10 10 8 8 2.37 1.78 1.18 1.48 1.18 1.60 1.18 1.48 1.18 1.39 1.47 1.89 1.18 2.07 2.37 2.13 Xij,  1.67 1.25 0.84 1.05 0.84 1.13 0.84 1.05 1.14 1.35 1.43 1.84 0.84 1.46 1.67 1.50 87 88 PL3 Chƣơng trình tính vị trí dung lƣợng bù tối ƣu xét đến tính ngẫu nhiên tải lập GAMS ****** KHAI BAO MANG DU LIEU VA DU LIEU DAU VAO ****** set i So nut /1*97/; alias (i,j); set t Thoi gian tinh toan /1*3/; set k So trang thai /1*16/; 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*97/ cp(i) Nut dau tu DG /1*97/ 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 du lieu CP /Qmin, Qmax/ Head6 He so cua tai va xac suat /kk0, xs0/ ; 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.12/; 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); PmnG(gen,t) = Generat(gen,"Pmin")/(SBase); QmxG(gen,t) = Generat(gen,"Qmax")/(SBase); QmnG(gen,t) = Generat(gen,"Qmin")/(SBase); TABLE Generat1(cp,Head5) Du lieu CP Qmin Qmax ; Parameter QmxG1(cp,t), QmnG1(cp,t); 88 89 QmxG1(cp,t) = 1*Generat1(cp,"Qmax")/(SBase); QmnG1(cp,t) = 1*Generat1(cp,"Qmin")/(SBase); TABLE Demand1(k,Head10) He so tai va xac suat kk0 xs0 ; Parameter kk(k), xs(k); kk(k) = Demand1(k,"kk0"); xs(k) = 1*Demand1(k,"xs0"); TABLE Demand(i, Head3) Cong suat phu tai PD0 QD0 ; Parameter PD(i,t), QD(i,t); PD(i,t)= 1.0*Demand(i,"PD0")*(1+GFD*(ord(t)-1))/(1000*SBase); QD(i,t)= 1.0*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")); Parameter G(i,j,t) , B(i,j,t) ; 89 90 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,k) Dien ap nut Delta(i,t,k) Goc pha dien ap nut P(i,t,k) Cong suat TD nhan tu nguon Q(i,t,k) Cong suat PK nhan tu nguon QG(i,t) Cong suat PK QG1(i,t) Cong suat PK Loss Tong chi phi U(n,i,t) Bien nhi phan ; Positive variable P, Q, QG; Binary variable U; Parameter VLevel(Gen) /1 1.05/; V.l(i,t,k) = 1.0; V.Fx(Gen,t,k) = VLevel(Gen); Delta.l(i,t,k) = 0.0; 90 91 Delta.fx("1",t,k) = 0.0; ************ HAM MUC TIEU VA RANG BUOC *************** Equations CostEq Equn1(i,t,k) Equn2(i,t,k) Equn3(i,t) Equn4(cp) Equn5 Equn6(i,t) ; ************** HAM MUC TIEU ************** CostEq Loss =e= 8760*0.5*Sum((i,j,t,k),xs(k)*G(i,j,t)*(V(i,t,k)*V(i,t,k) + V(j,t,k)*V(j,t,k) 2*V(i,t,k)*V(j,t,k)*cos(Delta(j,t,k)-Delta(i,t,k))) + QG(i,t)*0.0002*0)*Sbase*1000; ************** CAN BANG CONG SUAT NUT ************** Equn1(i,t,k)$(ord(i) gt 1) -QG(i,t)*0.0002- PD(i,t)*kk(k) =e= Sum(j, Y(i,j,t)*V(i,t,k)*V(j,t,k)*Cos(theta(i,j,t)+Delta(j,t,k) - Delta(i,t,k))); Equn2(i,t,k)$(ord(i) gt 1) QG1(i,t)- QD(i,t)*kk(k) =e= -Sum(j, Y(i,j,t)*V(i,t,k)*V(j,t,k)*Sin(theta(i,j,t)+Delta(j,t,k) - Delta(i,t,k))); ************** RANG BUOC NANG CAP THIET BI ************** Equn3(i,t)$(ord(t) gt 0) QG(i,t) =e= 1*(0.1*U("1",i,t)+0.15*U("2",i,t)+0.2*U("3",i,t)+0.3*U("4",i,t)+0.4 5*U("5",i,t)+0.05*U("6",i,t) +0.5*U("7",i,t)+0.6*U("8",i,t)+0.0*U("9",i,t)+1.2*U("10",i,t) +1.3*U("11",i,t)+1.0*U("12",i,t)); Equn4(cp) sum((n,t),U(n,cp,t)) =l= 1; Equn5 sum((n,cp,t),U(n,cp,t)) =l= 3; 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); ************ RANG BUOC DIEN AP NUT ********************** 91 92 V.Up(load,t,k) = 1.1; V.Lo(load,t,k) = 0.9; ************** GOI SLOVER TINH TOAN ********************* Option MINLP=BONMIN; Model OPF /all/; Solve OPF using MINLP Minimizing Loss; *************** TINH TOAN THONG SO VA HIEN THI ********** Display V.l, U.l, PDtt, Vnmin, Vnmax, DeltaP, Qg.l, loss.l; 92 ... NGUYỄN DUY TUỆ NGHIÊN CỨU GIẢI PHÁP GIẢM TỔN THẤT ĐIỆN NĂNG TRÊN LƯỚI ĐIỆN PHÂN PHỐI HUYỆN SỐP CỘP Chuyên ngành: Kỹ thuật điện Mã ngành: 8.52.02.01 LUẬN VĂN THẠC SĨ KHOA HỌC KỸ THUẬT ĐIỆN NGƯỜI HƯỚNG... tác giả luận văn: Nguyễn Duy Tuệ Đề tài luận văn: Nghiên cứu giải pháp giảm tổn thất điện lưới điện phân phối huyện Sốp Cộp Chuyên ngành: Kỹ thuật điện Mã số: : 8.52.02.01 Tác giả, Cán hướng dẫn... LĐPPTA Huyện Sốp Cộp Từ đó, đánh giá hiệu phƣơng pháp bù đồng thời đề xuất giải pháp bù cho LĐPPTA Huyện Sốp Cộp III Đối tƣợng phạm vi nghiên cứu - Đối tƣợng nghiên cứu: LĐPPTA Huyện Sốp Cộp -