Luận văn thạc sĩ 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 Phân loại lưới điện trung áp 1.1.2.2 Vai trò lưới điện trung áp 1.1.2.3 Các phần tử lưới điện trung áp 1.1.2.4 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 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 toá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 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 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] Đỗ 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 [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 to Supercapacitors-An Overview, International Journal of Electrochemical Science, Vol.3, 2008 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 M.Dixit, P.Kundu, H R.Jariwala, Optimal Allocation and Sizing of Shunt Capacitor in Distribution System for Power Loss Minimization, 2016SCEECS, India, 2016 [17] 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 [18] A.A.Eajal, M.E.El-Hawary, Optimal capacitor placement and sizing in distorted radial distribution systems part III: Numerical results, ICHQP2010, Italy, 2010 [19] [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 toá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 10 11 12 13 14 15 16 Tổng * Thông số đường dây Nút ij Sij* ,t , MVA 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 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 V.l(i,t,k) V.Fx(Gen,t,k) = VLevel(Gen); Delta.l(i,t,k) 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... 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... 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 Phạm