ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP MAI VIỆT DŨNG NGHIÊN CỨU VẬN HÀNH TỐI ƯU LƯỚI ĐIỆN PHÂN PHỐI KHI SỬ DỤNG CÁC NGUỒN ĐIỆN PHÂN TÁN VÀ ỨNG DỤNG TẠI MẠNG ĐIỆN TP CẨM PHẢ Chuyên ngành: Kỹ thuật điện Mã số: 60520202 LUẬN VĂN THẠC SĨ KỸ THUẬT THÁI NGUYÊN - 2016 i download by : skknchat@gmail.com LỜI CAM ĐOAN Tôi xin cam đoan vấn đề trình bày luận văn nghiên cứu riêng cá nhân tôi, hướng dẫn khoa học TS Nguyễn Minh Ý có tham khảo số tài liệu Các số liệu kết nghiên cứu nêu luận văn trung thực chưa công bố luận văn khác Tơi xin hồn tồn chịu trách nhiệm có sử dụng lại kết người khác Thái Nguyên, Ngày tháng năm 2016 Học viên Mai Việt Dũng ii download by : skknchat@gmail.com LỜI CẢM ƠN Tôi xin chân thành cảm ơn hướng dẫn nhiệt tình TS Nguyễn Minh Ý, giảng viên Bộ môn Hệ thống điện, Trường đại học kỹ thuật công nghiệp Thái Nguyên Người chịu trách nhiệm hướng dẫn tơi hồn thành luận văn Từ đáy lịng mình, tơi xin gửi lời cảm ơn sâu sắc tới thầy cô tham gia giảng dạy khóa học chuyên ngành Kỹ thuật điện, tạo điều kiện thuận lợi giúp tơi hồn thành khóa học Tơi xin gửi lời cảm ơn tới cán hành khoa Điện Phòng Đào tạo Sau đại học giúp đỡ chúng tơi q trình học tập trường Lời cuối cùng, chân thành cảm ơn động viên gia đình, bạn bè, đồng nghiệp người tạo điều kiện nhiều cho suốt chặng đường học tập qua Thái Nguyên, Ngày tháng năm 2016 Học viên Mai Việt Dũng iii download by : skknchat@gmail.com MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN iii MỤC LỤC iv DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT .vii DANH MỤC BẢNG viii DANH MỤC HÌNH VẼ ix MỞ ĐẦU 1 Mục đích nghiên cứu lý chọn đề tài Đối tượng nghiên cứu lý chọn đề tài Ý nghĩa khoa học tính thực tiễn đề tài Phương pháp nghiên cứu Chương TỔNG QUAN VỀ LƯỚI ĐIỆN PHÂN PHỐI TP CẨM PHẢ VÀ NGUỒN ĐIỆN PHÂN TÁN 1.1 Tổng quan lưới điện phân phối 1.1.1 Định nghĩa Lưới điện phân phối 1.1.2 Phân loại Lưới điện phân phối 1.1.3 Vai trò Lưới điện phân phối 1.2 Tổng quan lưới điện phân phối TP Cẩm Phả 1.2.1 Tổng quan TP Cẩm Phả 1.2.3 Nguồn điện cấp cho Cẩm Phả 1.2.4 Kết luận 1.3 Tổng quan nguồn điện phân tán (DG) 1.3.1 Định nghĩa nguồn điện phân tán 1.3.2 Đặc điểm công nghệ nguồn phát điện phân tán 1.3.3 Mạng điện Microgid 1.3.4 Một số loại nguồn phát điện phân tán 11 1.3.5 Xu hướng phát triển nguồn phân tán Việt Nam 28 1.3.6 Kết luận 31 iv download by : skknchat@gmail.com Chương TỐI ƯU HÓA VẬN HÀNH TRONG HỆ THỐNG ĐIỆN 33 2.1 Đặc tính tiêu hao nhiên liệu loại nguồn phát 33 2.1.1 Nguồn phát thủy điện 33 2.1.2 Nhà máy nhiệt điện 34 2.1.3 Đặc tính tiêu hao nhiên liệu 35 2.2 Tối ưu hóa vận hành truyền thống 36 2.3 Thị trường điện giá bán điện 37 2.3.1 Thị trường điện 37 2.3.2 Hiện trạng giá điện 42 2.4 Tối ưu hóa vận hành lưới điện thị trường điện 44 2.5 Mơ tả tốn nghiên cứu 45 2.6 Phương pháp giải toán tối ưu phối hợp tổ máy phát điện thị trường điện 49 2.7 Kết luận 52 Chương THUẬT TOÁN UNIT DECOMMITTMENT TRONG VẬN HÀNH TỐI ƯU CÁC TỔ MÁY PHÁT ĐIỆN 54 3.1 Đặt toán 54 3.1.1 Mô hình tốn 55 3.1.2 Tính tốn lợi nhuận thị trường điện cạnh tranh 57 3.2 Giải bái toán 61 3.3 Thuật toán Unit Decommittment 65 3.3.1 Tính tốn cơng suất dựa vào suất tăng công suất tổ máy 65 3.3.2 Tính tốn cơng suất dựa vào hàm chi phí tổ máy 69 3.3.3 Sơ đồ khối mô tả thuật toán 72 3.4 Kết luận 73 Chương ỨNG DỤNG THUẬT TOÁN TRONG BÀI TOÁN MINH HOẠ 74 4.1 Bài toán minh họa 74 4.1.1 Trường hợp thứ 76 4.1.2 Trường hợp thứ (TH2) 82 4.1.3 Trường hợp thứ (TH3) 88 v download by : skknchat@gmail.com 4.2 Kết luận 93 Chương KẾT LUẬN 94 5.1 Kết luận 94 5.2 Kiến nghị 95 TÀI LIỆU THAM KHẢO 96 PHỤ LỤC 98 vi download by : skknchat@gmail.com DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT CCĐ: Cung cấp điện CG: Nguồn phát điện trung tâm CN: Phụ tải công nghiệp DG: Nguồn điện phân tán ĐD: Đường dây ĐMT: Điện mặt trời LL: Tổng tổn thất công suất tác dụng đường dây lưới điện LPP: Lưới điện phân phối NLSK: Năng lượng sinh khối NLSK: Năng lượng sinh khối PCC: Điểm kết nối chung PQ: Nút phụ tải PV: Nút nguồn phát SH: Phụ tải sinh hoạt TĐ: Thủy Điện TĐN: Thuỷ điện nhỏ TM 1: Tổ máy nguồn phân tán số TM 2: Tổ máy nguồn phân tán số TM 3: Tổ máy nguồn phân tán số TM 4: Tổ máy nguồn phân tán số TM 5: Tổ máy nguồn phân tán số VI: Đại lượng đặc trưng cho chất lượng điện áp xuất tuyến VP: Hệ số cải thiện chất lượng điện áp lưới điện vii download by : skknchat@gmail.com DANH MỤC BẢNG Bảng 1.1 Dải công suất tương ứng công nghệ DG Bảng 1.2 Cường độ xạ trung bình tháng số khu vực phía Bắc 29 Bảng 1.3 Tổng xạ mặt trời Móng Cái 29 Bảng 1.4 Kế hoạch phát triển nguồn điện sử dụng lượng tái tạo 31 giai đoạn 2011 - 2020 có xét đến năm 2030 31 Phụ lục A: Thông số tổ máy phát điện 98 Phụ lục B: Thông số phụ tải 98 Phụ lục C: Thông số giá điện 99 Phục lục D: Công suất phát tổ máy sau tối ưu thuật toán 101 Phụ lục E: Mã lệnh thuật toán 104 viii download by : skknchat@gmail.com DANH MỤC HÌNH VẼ Hình 1.1 Sơ đồ phân loại nguồn phần tán Hình 1.2 Công nghệ điều khiển nối cho lưới điện siêu nhỏ Singapore 10 Hình 1.3 Cơng nghệ điều khiển nối cho lưới điện siêu nhỏ Hoa Kỳ 10 Hình 1.4 Sơ đồ nguyên lý hệ thống điện mặt trời 11 Hình 1.5 Hệ thống điện mặt trời lớn Việt Nam (tại TP Hồ Chí Minh) 12 Hình 1.6 Hình ảnh nhà máy điện gió 13 Hình 1.7 Nguyên lý cấu tạo tổ hợp tua-bin - máy phát điện gió 14 Hình 1.8 Nhà máy điện nhiệt kết hợp 17 Hình 1.9 Nhà máy điện nhiệt kết hợp công suất lớn 18 Hình 1.10 Sơ đồ nguyên lý máy phát điện tua-bin khí 19 Hình 1.11 Sơ đồ cấu tạo nguyên lý hoạt động pin nhiên liệu 20 Hình 1.12 Cơng trình xây dựng nhà máy thủy điện nhỏ 21 Hình 1.13 Nhà máy điện thủy triều kiểu đập cửa sông Rance (Pháp) 23 Hình 1.14 Hệ thống máy phát tua-bin thủy triều 23 Hình 1.15 Mơ hình phát điện sử dụng khí Biogass 25 Hình 1.16 Nhà máy điện sử dụng dạng lượng sinh khối 25 Hình 1.17 Nguyên lý sản xuất điện từ lượng địa nhiệt 27 Hình 1.18 Nhà máy địa nhiệt điện 27 Hình 1.19 Dự báo công suất nguồn phân tán Việt Nam đến năm 2030 30 Hình 2.1 Sơ đồ thị trường phát điện cạnh tranh 40 Hình 2.2 Sơ đồ thị trường bán bn điện cạnh tranh 40 Hình 2.3 Thị trường bán lẻ điện cạnh tranh 41 Hình 2.4 Biểu đồ giá điện 43 Hình 2.5 Lợi nhuận chi phí 47 Hình 2.6 Đường cong nhiên liệu đầu vào- công suất phát 48 Hình 3.1 Phối hợp tổ máy phát điện sở thị trường điện 56 Hình 3.2 Quan hệ Công suất Giá điện 58 Hình 3.3 Sơ đồ trạng thái tổ máy với Tjup = 3(h) Tịdown = 4(h) 61 Hình 3.4 Đồ thị giá điện 65 ix download by : skknchat@gmail.com Hình 3.5 Cơng suất phát theo giá điện λ 66 Hình 3.6 Công suất phát theo khả thực tế máy phát 67 Hình 3.7 Cơng suất phát thực tế máy phát 68 Hình 3.8 Quan hệ cơng suất lợi nhuận thực tế khả máy phát 69 Hình 3.9 Quan hệ giá điện công suất 70 Hình 3.10 Quan hệ cơng suất lợi nhuận 71 Hình 3.11 Quan hệ công suất lợi nhuận thực tế khả máy phát 72 Hình 4.1 Sơ đồ nguyên lý vị trí đặt nguồn điện 75 Hình 4.2 Đồ thị giá điện ngày 76 Hình 4.3 Đồ thị cơng suất tải ngày 77 Hình 4.4 Đồ thị cơng suất phát chưa tối ưu hóa 78 Hình 4.5 Đồ thị cơng suất phát tối ưu 79 Hình 4.6 Đồ thị chi phí vận hành tổ máy 80 Hình 4.7 Doanh thu 24h 81 Hình 4.8 Giá điện 24h cho TH2 82 Hình 4.9 Đồ thị cơng suất tải cho TH2 83 Hình 4.10 Đồ thị cơng suất phát chưa tối ưu cho TH2 84 Hình 4.11 Công suất phát tổ máy sau tối ưu cho TH2 85 Hình 4.12 Chi phí vận hành tổ máy TH2 86 Hình 4.13 Doanh thu 24h tổ máy 87 Hình 4.14 Giá điện 24h TH3 88 Hình 4.15 Cơng suất tải 24h TH3 89 Hình 4.16 Công suất phát chưa tối ưu TH3 90 Hình 4.17 Cơng suất phát tổ máy sau tối ưu TH3 91 Hình 4.18 Chi phí tổ máy TH3 92 Hình 4.19 Doanh thu ngày TH3 93 x download by : skknchat@gmail.com phát công suất lớn nhất, tổ máy cịn lại trì mức phát công suất 85% công suất định mức Tại thời điểm 23h 24h tổ máy 4, tắt, tổ máy 1, 2, phát với mức công suất gần nhỏ Chi phí vận hành tổ máy kVNĐ 18000 16000 14000 12000 10000 8000 6000 4000 2000 h TM TM 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 TM TM TM Nguồn lưới Hình 4.18 Chi phí tổ máy TH3 Hình 4-18 thể chi phí phát tổ máy trường hợp 3, hình 4-18 có đường đặc tính gần giống hình 4-17 chi phí phát phụ thuộc vào công suất phát tổ máy Chỉ thời điểm khởi động máy lúc 9h 10h chi phí phát lúc cịn kèm theo chi phí khởi động tổ máy 92 download by : skknchat@gmail.com Doanh thu 24h tổ máy kVNĐ 25000 20000 15000 10000 5000 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 h Hình 4.19 Doanh thu ngày TH3 Doanh thu ngày là: 322167 kVNĐ Lợi nhuận thu ngày: 7120 kVNĐ Nhận xét: Trong trường hợp từ thời điểm 1h đến 8h tổ máy tắt, thời điểm 9h có tổ máy thứ bật phát với công suất gần nhỏ Tại thời điểm từ 1h đến 9h giá điện gần mức thấp Từ 10h đến 20h tổ máy 1, 2, 3, bật phát với mức cơng suất trung bình, tổ máy chế độ tắt lúc giá điện mức trung bình Tại thời điểm 21h giá thành điện mức cao tất tổ máy bật phát với công suất gần lớn Nhưng sang thời điểm 22h giá điện mức trung bình tổ máy thứ khơng tắt có ràng buộc bật tắt tổ máy Vào 23h 24h tổ máy 1, bật công suất gần nhỏ tổ máy 4, chế độ tắt 4.2 Kết luận Với việc áp dụng vào tính tốn cho lưới điện khu vực trung tâm thành phố Cẩm Phả, số liệu đưa vào phần mềm MATLAB xử lý, kết ta thu chí phí vận hành giảm, lợi nhuận thu từ bán điện tăng lên Công suất phát tối ưu nhắm đảm bảo lợi nhuận ln lớn nhất, từ ta xây dựng kế hoạch phát điện nguồn phân tán ngày 93 download by : skknchat@gmail.com Chương KẾT LUẬN 5.1 Kết luận Với việc phân tích nguồn phân tán, xu hướng thị trường điện giá bán điện tương lai, luận văn xây dựng mơ hình tốn học toán nguồn phân tán vận hành thị trường điện với biểu giá điện 24h, tính tốn phối hợp tổ máy để chi phí đạt cực tiểu, nâng cao lợi nhuận phát điện lớn Sử dụng mơ hình tốn học để giải toán phân bố tối ưu tổ máy phát điện Đưa định kế hoạch vận hành tổ máy phát điện dựa vào điều kiện ràng buộc: - Về thơng số chi phí tổ máy phát điện - Về chi phí khởi động tổ máy phát điện với cấp trạng thái tổ máy khác giá khởi động tổ máy khác - Về thời gian bật tắt máy - Về ảnh hưởng tổn thật phát nguồn - Về giới hạn phát công suất tổ máy Sau lập kế hoạch vận hành, ta có trạng thái, cơng suất phát, chi phí lợi nhuận thu cụ thể tổ máy thời điểm khảo sát Kết tính tốn, giúp ta xác định được: - Các ảnh hưởng DG lưới làm thay đổi tổng tổn thất cơng suất lưới - Bên cạnh ảnh hưởng tới mơi trường tính kinh tế lưới điện - Phương pháp cho phép người vận hành nghiên cứu mạng phân phối bất kỳ, sử dụng thơng tin có sẵn để lập kế hoạch cho kết nối DG nhằm đạt mục tiêu tối thiểu hóa chi phí vận hành, cải thiện điện áp Do khả thời gian có hạn, đâu lĩnh vực tương đối nên nội dung luận văn chi tập trung nêu vấn đề nguồn phân tán lưới điện phân phối - Trạng thái cụ thể tổ máy thời điểm t (giờ) suốt chu kỳ khảo sát - Thể giá trị công suất phát tổ máy thời điểm t 94 download by : skknchat@gmail.com - Chi phí tổ máy đạt trình vận hành - Lợi nhuận tổ máy đạt trình vận hành - Tổng lợi nhuận hệ thống vận hành theo kế hoạch tính tốn - Từ kết ta đưa phương pháp vận hành tổ máy có nguồn phân tán ngày cho đạt lợi nhuận cao tổ máy - Ngoài tốn cịn đưa tới thị trường điện tương lại thị trường cạnh tranh giá điện cạnh tranh 5.2 Kiến nghị - Luận văn đưa đưa giá điện tham khảo trường điện Mỹ chưa tham khảo giá điện Việt Nam biểu giả điện Việt Nam có mức, nêu tiến tới thị trường điện cạnh tranh giá điện có biểu giá 24h - Luận văn chưa xét tới ảnh hưởng điện áp, ảnh hưởng tổn thất, ảnh hưởng vị trí nguồn phân tán kết cấu lưới điện vận hành nguồn phân tán Cần có nghiên cứu sâu ảnh hưởng nguồn phân tán tới chế độ vận hành lưới điện - Chưa xét tới việc ảnh hưởng nguồn phân tán tới độ tin cậy cung cấp điện lưới điện Đây hướng nghiên cứu tương lai - Luận văn nghiên cứu vận hành nguồn phân tán có khả điều khiển (là nguồn phân tán có tiêu thụ nhiên liệu) mà chưa xét đến nguồn lượng tái tạo phụ thuộc vào điều kiện tư nhiên 95 download by : skknchat@gmail.com TÀI LIỆU THAM KHẢO [1] Trần Bách (2008), Lưới điện Hệ thống điện, tập 1, 2, 3, NXB Khoa học Kỹ thuật, Hà Nội [2] Trần Bách (1999), Tối ưu hoá chế độ Hệ thống điện, Khoa Đại học Tại chức - Trường Đại học Bách khoa Hà Nội, Hà Nội [3] Bộ Công thương (2011), Tiêu chuẩn kỹ thuật điện Việt Nam 2011, Nhà xuất Lao động [4] Bộ Công thương (2015), Thông tư số 39/2015/TT-BCT ngày 18 tháng 11 năm 2015 Bộ Công thương “Quy định Hệ thống điện phân phối”, Hà Nội [5] Nguyễn Văn Đạm (2005), Mạng lưới điện 1, Nhà xuất Khoa học Kỹ thuật, Hà Nội [6] Lê Kim Hùng, Lê Thái Thanh, Tối ưu hố vị trí đặt cơng suất nguồn phân tán mơ hình lưới phân phối 22kV, Tạp chí Khoa học Cơng nghệ, Đại học Đà Nẵng 2008 [7] Nguyễn Công Hiền, Nguyễn Mạch Hoạch (2007), Hệ thống cung cấp điện xí nghiệp cơng nghiệp nhà cao tầng, NXB KHKT, Hà Nội [8] Nguyễn Xuân Phú (1998), Cung cấp điện, NXB KHKT, Hà Nội [9] Nguyễn Quang Khải, Những vấn đề phát triển lượng sinh khối Việt Nam, Báo cáo Hội thảo phát triển lượng bền vững Việt Nam [10] Đặng Đình Thống, Lê Danh Liên, Cơ sở lượng tái tạo, Đại học Bách khoa Hà Nội, NXB Khoa học Kỹ thuật, Hà Nội 2006 [11] Viện Năng lượng, Quy hoạch phát triển điện lực Quốc Gia giao đoạn 2006 2015, có xét triển vọng đến năm 2025 [12] Allen J Wood, Bruce F Wollenberg, “Power Generation, Operation and Control”, John Wiley and Sons, New York, New York, 1984 [13] Ioannis G Damousis, Anastasios G Bakirtzis, Petros S Dokopoulos, “A Solution to the Unit-Commitment Problem Using Integer-Coded Genetic Algorithm”, IEEE Transactions on Power Systems, Vol 19, No 2, May 2004 [14] Sayeed Salam, “Unit Commitment Solution Methods”, WorldAcademyaA of Science, Engineering and Technology 35 2007 96 download by : skknchat@gmail.com [15] Farid Benhamida, E N Abdallah and A H Rashed, “Thermal Unit Commitment Solution Using an ImprovedLagrangian Relaxation”, Department of Electrical Engineering, Alexandria University [16] Weerakorn Ongsakul, Nit Petcharaks, “Unit Commitment by Enhanced Adaptive Lagrangian Relaxation”, IEEE Transactions on Power Systems, Vol 19, No 1, February 2004 [17] T Senjyu, A.Y Saber, T Miyagi, K Shimabukuro, N Urasaki and T.Funabashi, “Fast Technique for Unit Commitment by Genetic Algorithm Based on Unit Clustering”, IEE Proc.-Gener Transm Distrib., Vol 152, No.5, September 2005 [18] V.N Dieu and W Ongsakul, “Improved Merit Order and Augmented Lagrange Hopfield Network for Unit Commitment”, IET Gener Transm Distrib., 2007, 1, (4), pp 548-556 [19] Joao Catalao, Silvio Mariano, Victor Mendes, Luís Ferreira, “Unit Commitment with Environmental Considerations: A Practical Approach”, Session 18, Paper 3, Page 1, 15th PSCC, Liege, 22-26 August 2005 [20] Calvin Jin San Chan, “Development of Profit Maximisation Unit Commitment Program”, MSc Thesis, The University of Manchester Institute of Science and Technology, 2000 [21] C K Pang, H C Chen “Optimal Short-Term Thermal Unit Commitment”, IEEE Transactions on Apparatus and Power Systems, Vol PAS-95, No.4, pp.1336-1346, Jul./Aug 1976 [22] Chen, H Wong, K.P Wang, X Chung, C.Y ”A Coevolutionary Approach to Modeling Oligopolistic Electricity Markets” In Proceedings of 2005 IEEE Power Engineering Society General Meeting, pp.230 - 236 Vol [23] Phan Thanh Tú, “Luận văn thạc sĩ”, Trường Đại học SPKT Thành phố Hồ Chí Minh 97 download by : skknchat@gmail.com PHỤ LỤC Phụ lục A: Thông số tổ máy phát điện Tổ máy Tổ máy Khả phát tổ máy P_mix (kW) P_max (kW) TM TM TM TM TM Nguồn lưới 300 300 400 450 500 -5000 800 900 1000 1000 1100 5000 Thời gian t_on (h) 2 3 bật tắt t_off (h) 2 2 a 0.0005 0.0005 0.0003 0.0003 0.0003 b 1.4 1.5 1.5 1.6 1.7 c 40 50 50 60 60 95 100 100 105 105 0.02 0.02 0.03 0.04 0.04 Hằng số chi phí Chi phí Cst khởi động (kVNĐ) Tổn thất Lf (%) 0.005 Phụ lục B: Thông số phụ tải Thời gian (h) Công suất tải (kW) Thời gian (h) Công suất tải (kW) 4833 13 9155 4814 14 8093 4790 15 7183 4876 16 6469 5190 17 6560 5950 18 6714 9269 19 6671 9071 20 8179 9069 21 9017 10 9124 22 8336 11 9619 23 6214 12 9552 24 5640 98 download by : skknchat@gmail.com Phụ lục C: Thông số giá điện Trường hợp thứ nhất: Thời gian (h) Giá điện (kVNĐ Thời gian (h) Giá điện (kVNĐ 1.644 13 1.863 1.642 14 1.747 1.568 15 1.672 1.578 16 1.608 1.679 17 1.606 1.770 18 1.678 1.908 19 1.789 2.013 20 2.136 2.328 21 2.608 10 2.288 22 2.233 11 2.083 23 1.738 12 1.937 24 1.685 Thời gian (h) Giá điện (kVNĐ Thời gian (h) Giá điện (kVNĐ 1.557 13 2.519 1.513 14 2.396 1.423 15 2.216 1.443 16 2.111 1.518 17 2.175 1.705 18 2.166 2.325 19 2.028 2.350 20 2.044 2.476 21 2.323 10 2.569 22 2.115 11 2.628 23 1.771 12 2.588 24 1.640 Trường hợp thứ hai: 99 download by : skknchat@gmail.com Trường hợp thứ ba: Thời gian (h) Giá điện (kVNĐ Thời gian (h) Giá điện (kVNĐ 1.538 13 2.018 1.504 14 1.969 1.304 15 1.918 1.221 16 1.921 1.325 17 1.973 1.428 18 2.025 1.412 19 1.998 1.538 20 2.021 1.776 21 2.425 10 1.947 22 2.154 11 1.988 23 1.875 12 2.025 24 1.549 100 download by : skknchat@gmail.com Phục lục D: Công suất phát tổ máy sau tối ưu thuật toán Trường hợp thứ nhất: Thời gian TM1 TM2 TM3 TM4 TM5 Nguồn lưới (h) (kW) (kW) (kW) (kW) (kW) (kW) 0 0 4858 0 0 4838 0 0 4815 0 0 4901 0 0 5217 343 0 0 5642 480 380 601 0 7884 582 482 770 570 6769 800 793 1000 1000 911 4727 10 800 754 1000 1000 846 4883 11 651 551 884 682 515 6467 12 507 407 647 0 8069 13 435 335 528 0 7928 14 320 0 0 7818 15 300 0 0 6924 16 300 0 0 6206 17 300 0 0 6297 18 300 0 0 6453 19 362 407 0 5951 20 704 604 970 768 601 4665 21 800 900 1000 1000 1100 4387 22 800 700 1000 925 758 4303 23 312 300 400 0 5253 24 300 300 400 0 4688 101 download by : skknchat@gmail.com Trường hợp thứ hai: Thời gian TM1 TM2 TM3 TM4 TM5 Nguồn lưới (h) (kW) (kW) (kW) (kW) (kW) (kW) 0 0 4858 0 0 4838 0 0 4815 0 0 4901 0 0 5217 0 0 5980 800 790 1000 1000 905 4936 800 815 1000 1000 946 4675 800 900 1000 1000 1100 4439 10 800 900 1000 1000 1100 4494 11 800 900 1000 1000 1100 4992 12 800 900 1000 1000 1100 4925 13 800 900 1000 1000 1100 4525 14 800 860 1000 1000 1019 3576 15 782 682 1000 896 730 3233 16 679 579 930 728 561 3113 17 742 642 1000 831 664 2812 18 733 633 1000 816 649 3013 19 598 498 796 595 500 3794 20 613 513 821 620 500 5228 21 800 788 1000 1000 903 4687 22 683 583 936 734 568 4961 23 344 300 400 0 5221 24 300 300 400 0 4688 102 download by : skknchat@gmail.com Trường hợp thứ ba: Thời gian TM1 TM2 TM3 TM4 TM5 Nguồn lưới (h) (kW) (kW) (kW) (kW) (kW) (kW) 0 0 4858 0 0 4838 0 0 4815 0 0 4901 0 0 5217 0 0 5980 0 0 9316 0 0 9117 349 0 0 8771 10 517 417 663 464 7155 11 558 458 731 531 7442 12 594 494 790 590 7189 13 588 488 779 579 6823 14 539 439 699 500 6005 15 489 389 616 450 5321 16 492 392 621 450 4592 17 543 443 705 505 4447 18 594 494 790 590 4336 19 567 467 746 545 4433 20 590 490 783 583 5829 21 800 888 1000 1000 1066 4431 22 722 622 1000 797 631 4702 23 447 347 546 0 4931 24 300 300 400 0 4688 103 download by : skknchat@gmail.com Phụ lục E: Mã lệnh thuật toán % - Unit decommittment -N_stage = length(Price); N_g = length(P_min) - 1; t = [1:N_stage]; for k = 1:N_stage for i = 1:N_g u(k,i) = 1; end; end; for k = 1:N_stage Miu(k) = Price(k) / (1 - lf(N_g+1)); for i = 1:N_g P(k,i) = (Miu(k) * (1 - lf(i)) - b(i)) / (2 * a(i)); if P(k,i) > P_max(i) P(k,i) = P_max(i); elseif P(k,i) < P_min(i) P(k,i) = P_min(i); end; end; P(k,N_g+1) = (P_load(k) + P(k, 1:N_g) * (lf(1:N_g))'- sum(P(k, 1:N_g))) / (1 - lf(N_g+1)); end; for k = 1:N_stage for i = 1:N_g C_fuel(k,i) = FuelCost(P(k,i), a(i), b(i), c(i)); end; P(k,N_g+1) = (P_load(k) + P(k, 1:N_g) * (lf(1:N_g))'- sum(P(k, 1:N_g))) / (1 - lf(N_g+1)); C_feeder(k) = Price(k) * P(k,N_g+1); C_unit(k) = sum(C_fuel(k,:)) + C_feeder(k); end; C_unit_total = sum(C_unit); 104 download by : skknchat@gmail.com for k = 1:N_stage for i = 1:N_g miu(k,i) = (P(k,i) * * a(i) + b(i)) / (1 - lf(i)); end; end; miu_de = miu; u_de = u; P_de = P; for i = 1:N_g for k = 1:N_stage if miu_de(k,i) > Miu(k) count = 1; cont = 1; for j = k+1:N_stage if miu_de(j,i) > Miu(j) && cont == count = count + 1; else cont = 0; end; end; time_off = max(t_off(i), count); T_off = min(N_stage, k+time_off-1); C_save = 0; for j = k:T_off C_save = C_save + FuelCost(P_de(j,i), a(i), P_de(j,i) * (1 - lf(i)) * Miu(j); end; C_save = C_save - Cst(i); if C_save >= for j = k:T_off u_de(j,i) = 0; P_de(j,i) = 0; miu_de(j,i) = Miu(j); end; end; time_off = 0; end; 105 download by : skknchat@gmail.com b(i), c(i)) - end; end; for k = 1:N_stage P_de(k,N_g+1) = (P_load(k) + P_de(k, 1:N_g) * (lf(1:N_g))'- sum(P_de(k, 1:N_g))) / (1 - lf(N_g+1)); end; for i = 1:N_g u_pre(i) = 1; end; for k = 1:N_stage for i = 1:N_g if u_de(k,i) == C_fuel_de(k,i) = FuelCost(P_de(k,i), a(i), b(i), c(i)); if u_pre(i) == C_fuel_de(k,i) = C_fuel_de(k,i) + Cst(i); end; else C_fuel_de(k,i) = 0; end; u_pre(i) = u_de(k,i); end; C_feeder_de(k) = Price(k) * P_de(k,N_g+1); C_deunit(k) = sum(C_fuel_de(k,:)) + C_feeder_de(k); end; C_deunit_total = sum(C_deunit); 106 download by : skknchat@gmail.com ... Nguyên, Ngày tháng năm 2016 Học viên Mai Việt Dũng ii download by : skknchat@gmail.com LỜI CẢM ƠN Tôi xin chân thành cảm ơn hướng dẫn nhiệt tình TS Nguyễn Minh Ý, giảng viên Bộ môn Hệ thống điện,... năm 2016 Học viên Mai Việt Dũng iii download by : skknchat@gmail.com MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN iii MỤC LỤC iv DANH MỤC CÁC KÝ HIỆU, CHỮ VIẾT TẮT... xin cam đoan vấn đề trình bày luận văn nghiên cứu riêng cá nhân tôi, hướng dẫn khoa học TS Nguyễn Minh Ý có tham khảo số tài liệu Các số liệu kết nghiên cứu nêu luận văn trung thực chưa công