Tổng quan nghiên cứu

Trong cơ cấu nguồn điện hiện đại, nhiệt điện truyền thống chiếm tỷ trọng lớn từ 50% đến 60% tổng sản lượng điện năng cung cấp cho nền kinh tế quốc dân. Tuy nhiên, các nguồn nhiên liệu hóa thạch như than đá, dầu diesel và khí tự nhiên đang đối mặt với nguy cơ cạn kiệt nghiêm trọng, trong khi chi phí nhiên liệu chiếm khoảng 60% đến 70% tổng chi phí vận hành của toàn hệ thống điện. Bài toán điều độ kinh tế (Economic Dispatch - ED) đóng vai trò then chốt nhằm phân bổ tối ưu công suất phát giữa các nhà máy nhiệt điện sao cho tổng chi phí nhiên liệu đạt mức thấp nhất, đồng thời đáp ứng trọn vẹn nhu cầu phụ tải và thỏa mãn nghiêm ngặt các giới hạn kỹ thuật vận hành.

Mục tiêu cốt lõi của nghiên cứu là xây dựng, lập trình và kiểm chứng thuật toán tìm kiếm dựa trên lực hấp dẫn (Gravitational Search Algorithm - GSA) trên môi trường MATLAB để giải quyết bài toán điều độ kinh tế phức tạp. Đề tài tập trung xử lý các đặc tính phi tuyến thực tế của nhà máy điện như hiệu ứng điểm van công suất (valve-point loading effects), vận hành đa nhiên liệu, giới hạn tốc độ tăng giảm công suất (ramp rate) và các vùng cấm vận hành.

Nghiên cứu được triển khai thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia TP.HCM trong giai đoạn từ tháng 1 năm 2013 đến tháng 9 năm 2014. Phạm vi thử nghiệm bao gồm hệ thống 10 máy phát nhiệt điện không xét điểm van, hệ thống chuẩn IEEE 13 tổ máy và hệ thống quy mô lớn IEEE 40 tổ máy phát nhiệt điện có xét hiệu ứng điểm van. Kết quả nghiên cứu mang lại ý nghĩa kinh tế kỹ thuật to lớn, giúp giảm chi phí nhiên liệu từ 0,85% đến 1,5% so với các giải thuật truyền thống, rút ngắn hơn 20% thời gian tính toán và nâng cao độ an toàn vận hành hệ thống điện.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng lý thuyết tối ưu hóa hệ thống điện kết hợp với lý thuyết trí tuệ nhân tạo hiện đại để giải quyết bài toán điều độ kinh tế:

Thứ nhất, mô hình toán học bài toán điều độ kinh tế (ED): Hàm mục tiêu tổng quát hướng đến cực tiểu hóa tổng chi phí nhiên liệu của N tổ máy phát điện trong hệ thống. Trong điều kiện vận hành lý thuyết, hàm chi phí của từng tổ máy được mô tả bằng hàm bậc hai kinh điển. Tuy nhiên, khi xét đến hiệu ứng đóng mở của các van hơi nước, hàm chi phí xuất hiện các gợn sóng phi tuyến phức tạp và được mô hình hóa bằng việc cộng thêm thành phần dao động sin trị tuyệt đối. Hệ thống ràng buộc bao gồm cân bằng công suất tác dụng có tính đến tổn thất đường dây thông qua ma trận tổn thất B-matrix kích thước N x N, giới hạn công suất cực tiểu và cực đại, giới hạn ramp rate và các vùng cấm vận hành (Prohibited Operating Zones).

Thứ hai, lý thuyết thuật toán tìm kiếm dựa trên lực hấp dẫn (GSA): Được xây dựng dựa trên định luật vạn vật hấp dẫn và định luật II Newton. Trong không gian tìm kiếm N chiều, mỗi giải pháp điều độ công suất được xem như một vật thể có khối lượng quán tính tỷ lệ thuận với giá trị độ phù hợp (fitness). Khối lượng càng lớn tương ứng với phương án có chi phí nhiên liệu càng thấp, tạo ra lực hấp dẫn mạnh mẽ để hút các vật thể khác về phía mình, định hướng toàn bộ quần thể tiến nhanh về điểm tối ưu toàn cục.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu kỹ thuật từ các hệ thống điện chuẩn quốc tế được công bố bởi IEEE, bao gồm hệ thống 10 tổ máy không xét điểm van, hệ thống 13 tổ máy phát với 5 điểm van công suất và hệ thống quy mô lớn 40 tổ máy phát nhiệt điện phục vụ mức phụ tải cực đại lên tới 10.500 MW.

Cỡ mẫu mô phỏng được thiết lập gồm 100 lần chạy độc lập cho từng kịch bản thử nghiệm nhằm bảo đảm tính chính xác và độ tin cậy về mặt thống kê. Phương pháp chọn mẫu áp dụng kỹ thuật khởi tạo quần thể ngẫu nhiên có kiểm soát không gian tìm kiếm với kích thước bầy từ 30 đến 50 cá thể vật thể.

Lý do lựa chọn thuật toán GSA là vì các phương pháp toán học cổ điển như nguyên lý Lagrange hay Gradient rất dễ bị kẹt tại các điểm cực trị địa phương hoặc không thể hội tụ khi gặp hàm chi phí phi tuyến, không lồi và gián đoạn. Toàn bộ thuật toán được lập trình trên phần mềm MATLAB R2013, kết hợp kỹ thuật hàm phạt để xử lý các ràng buộc bất đẳng thức. Tiến trình nghiên cứu được triển khai xuyên suốt từ tháng 1 năm 2013 đến khi nghiệm thu vào tháng 9 năm 2014.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng thực nghiệm trên các hệ thống điện thử nghiệm đã ghi nhận bốn phát hiện khoa học quan trọng:

Thứ nhất, đối với hệ thống 10 máy phát nhiệt điện tiêu chuẩn với mức tải nhu cầu từ 2.000 MW đến 2.700 MW, thuật toán GSA đã tìm ra phương án phân bổ công suất tối ưu giúp giảm tổng chi phí nhiên liệu khoảng 0,85% so với phương pháp Gradient và tiết kiệm 0,92% so với thuật toán Di truyền (GA).

Thứ hai, trên hệ thống chuẩn IEEE 13 tổ máy phát có xét hiệu ứng điểm van phức tạp ở mức phụ tải 1.800 MW và 2.520 MW, GSA đạt mức chi phí vận hành tối ưu 17.960 $/h, thấp hơn khoảng 1,5% so với thuật toán Tối ưu hóa bầy đàn (PSO) cơ bản, đồng thời triệt tiêu hoàn toàn hiện tượng dao động không hội tụ do 5 điểm van gây ra.

Thứ ba, trên hệ thống mạng lưới quy mô lớn IEEE 40 tổ máy phát với đặc tính đa cực trị cực kỳ phức tạp (phụ tải 10.500 MW), GSA đạt mức chi phí nhiên liệu tối ưu xấp xỉ 121.412 $/h, vượt trội hơn các phương pháp ABC và GA khoảng 0,94%, đồng thời duy trì thời gian thực thi tính toán trung bình dưới 45 giây cho 500 thế hệ lặp.

Thứ tư, thuật toán GSA thể hiện độ ổn định vượt bậc với tỷ lệ tìm thấy nghiệm tối ưu thành công đạt 98% qua 100 lần chạy ngẫu nhiên, độ lệch chuẩn chi phí giữa các lần chạy chỉ dao động ở mức cực nhỏ từ 0,02% đến 0,05%.

Thảo luận kết quả

Hiệu năng vượt trội của GSA xuất phát từ cơ chế tương tác lực hấp dẫn giúp cân bằng xuất sắc giữa hai giai đoạn then chốt: thăm dò diện rộng (exploration) ở những bước lặp đầu và khai thác sâu (exploitation) ở những bước lặp cuối. Trong khi các giải thuật như GA dễ mất đa dạng di truyền và PSO dễ bị bẫy cục bộ bởi vận tốc hạt quá lớn, GSA tự động giảm dần hằng số hấp dẫn theo thời gian, cho phép các cá thể hội tụ êm thuận về cực tiểu toàn cầu.

Dữ liệu mô phỏng được trực quan hóa thông qua biểu đồ đường cong hội tụ (Convergence Characteristics) biểu diễn sự suy giảm của chi phí nhiên liệu theo số vòng lặp từ 1 đến 500, cùng bảng đối sánh chi tiết mức huy động công suất từng tổ máy từ P1 đến P40. Kết quả đối chiếu với các công trình nghiên cứu quốc tế sử dụng PSO, ACO hay Cuckoo Search (CS) cho thấy GSA giảm được khoảng 30% dung lượng bộ nhớ máy tính và có cấu trúc mã nguồn đơn giản hơn đáng kể.

Đề xuất và khuyến nghị

Nhằm chuyển hóa các kết quả nghiên cứu thành giá trị thực tiễn trong quản lý và vận hành hệ thống điện, bốn giải pháp cụ thể được đề xuất:

Thứ nhất, tích hợp thuật toán GSA vào hệ thống SCADA/EMS tại Trung tâm Điều độ Hệ thống điện Quốc gia và các trung tâm điều độ miền. Mục tiêu đặt ra là tối ưu hóa biểu đồ huy động công suất ngày tới (Day-ahead) và thời gian thực, giúp giảm từ 1,0% đến 2,5% chi phí mua nhiên liệu hàng ngày, thực hiện trong lộ trình 12 tháng.

Thứ hai, phát triển mô hình thuật toán lai ghép giữa GSA và PSO (PSOGSA) nhằm tăng tốc độ xử lý thêm 20% đến 35%, đáp ứng hoàn hảo chu kỳ điều độ tức thời 5 đến 15 phút. Nhóm nghiên cứu tại các trường đại học khối kỹ thuật chủ trì hoàn thiện trong vòng 6 đến 9 tháng.

Thứ ba, mở rộng mô hình bài toán sang điều độ kinh tế - môi trường (Combined Economic Emission Dispatch - CEED), tích hợp hàm mục tiêu cắt giảm 10% đến 15% lượng khí thải CO2 và NOx. Tập đoàn Điện lực Việt Nam chủ trì thực hiện giai đoạn 2024 - 2026.

Thứ tư, chuẩn hóa cơ sở dữ liệu ngõ vào - ngõ ra và đặc tính điểm van của 100% các tổ máy nhiệt điện than, tuabin khí trên toàn hệ thống điện quốc gia. Các Tổng công ty Phát điện (GENCO 1, 2, 3) chịu trách nhiệm hoàn thành trong thời hạn 18 tháng.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

Nhóm 1: Kỹ sư điều độ và chuyên gia lập phương thức vận hành hệ thống điện tại các Trung tâm Điều độ. Nghiên cứu cung cấp công cụ tính toán phân bổ công suất chính xác, giúp tiết kiệm hàng tỷ đồng chi phí nhiên liệu hàng năm và đảm bảo an ninh năng lượng.

Nhóm 2: Cán bộ quản lý kỹ thuật tại các nhà máy nhiệt điện và Tổng công ty Phát điện. Luận văn giúp các đơn vị nắm bắt chính xác đặc tính tiêu hao nhiên liệu, tối ưu hóa điểm vận hành của từng lò hơi - tuabin và nâng cao hiệu suất nhiệt của nhà máy thêm 1,5% đến 3,0%.

Nhóm 3: Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Hệ thống điện, Kỹ thuật điện và Tự động hóa. Tài liệu là nguồn tham khảo học thuật chất lượng cao về phương pháp mô hình hóa bài toán phi tuyến phức tạp và ứng dụng thuật toán tiến hóa trên MATLAB với hơn 90 trang phân tích chi tiết.

Nhóm 4: Các công ty công nghệ và đơn vị phát triển phần mềm năng lượng thông minh (Smart Grid). Đề tài cung cấp cấu trúc thuật toán tối ưu để nhúng trực tiếp vào các module phần mềm quản trị năng lượng thế hệ mới.

Câu hỏi thường gặp

Thuật toán GSA có nguyên lý hoạt động cơ bản như thế nào trong bài toán điều độ kinh tế? GSA mô phỏng định luật vạn vật hấp dẫn, xem mỗi phương án phân bổ công suất như một vật thể trong không gian N chiều. Khối lượng vật thể tỷ lệ với nghịch đảo chi phí nhiên liệu. Lực hấp dẫn tương tác kéo các vật thể về phía vật thể có khối lượng lớn nhất, giúp hệ thống 40 tổ máy nhanh chóng hội tụ về điểm chi phí thấp nhất 121.412 $/h sau 500 vòng lặp.

Hiệu ứng điểm van công suất ảnh hưởng như thế nào đến bài toán điều độ? Khi mở các van hơi nước để tăng công suất, hiện tượng nghẽn dòng tạo ra các đường cong chi phí lượn sóng phi tuyến với nhiều điểm cực tiểu cục bộ. Điều này khiến các phương pháp giải tích truyền thống như Lagrange hay Gradient dễ bị sai lệch nghiệm từ 3% đến 8%, đòi hỏi phải sử dụng các thuật toán tiến hóa như GSA để giải quyết.

Thuật toán GSA có ưu thế gì vượt trội so với GA và PSO? GSA làm việc trực tiếp trên không gian biến thực liên tục, không đòi hỏi mã hóa nhị phân phức tạp như GA và không bị phụ thuộc vào đạo hàm hàm mục tiêu. So với PSO, GSA có khả năng duy trì tính đa dạng quần thể tốt hơn, giảm 30% dung lượng bộ nhớ máy tính và cho kết quả chi phí thấp hơn khoảng 1,5% trên hệ thống chuẩn 13 tổ máy.

Ràng buộc cân bằng công suất và tổn thất lưới điện được xử lý ra sao trong nghiên cứu? Nghiên cứu sử dụng công thức tổn thất ma trận B để xác định chính xác tổn thất công suất truyền tải. Ràng buộc cân bằng công suất được giải quyết bằng cách thiết lập một máy phát làm nút cân bằng kết hợp với hàm phạt thích nghi, đảm bảo sai số cân bằng công suất luôn được kiểm soát chặt chẽ dưới ngưỡng 0,001 MW.

Thuật toán GSA có thể mở rộng cho các hệ thống tích hợp năng lượng tái tạo không? Hoàn toàn khả thi. Bằng việc bổ sung các hàm chi phí không chắc chắn cho nguồn điện gió và điện mặt trời, GSA có thể mở rộng giải bài toán điều độ cho các hệ thống điện lai quy mô trên 50 nút nguồn với tỷ lệ hội tụ thành công duy trì trên mức 95%.

Kết luận

• Luận văn giải quyết xuất sắc bài toán điều độ kinh tế phi tuyến phức tạp có xét hiệu ứng điểm van, đa nhiên liệu và vùng cấm vận hành. • Ứng dụng thành công thuật toán GSA trên 3 hệ thống điện chuẩn (10 máy phát, IEEE 13 máy và IEEE 40 tổ máy phát nhiệt điện). • Đạt hiệu quả tiết kiệm chi phí nhiên liệu vượt trội từ 0,85% đến 1,5% so với các phương pháp Gradient, GA và PSO. • Đảm bảo độ tin cậy và tính ổn định thuật toán rất cao với tỷ lệ hội tụ đạt 98% qua 100 lần thử nghiệm độc lập. • Xây dựng hoàn chỉnh bộ chương trình mô phỏng trên MATLAB có khả năng ứng dụng thực tiễn cao cho hệ thống điện Việt Nam.

Đóng góp quan trọng nhất của luận văn là chứng minh tính khả thi và ưu thế vượt trội của thuật toán GSA trong việc giải bài toán tối ưu hóa hệ thống điện phi tuyến quy mô lớn. Kế hoạch tiếp theo là phát triển thuật toán lai PSOGSA và ứng dụng thử nghiệm tại các trung tâm điều độ trong giai đoạn 2024 - 2025. Các đơn vị vận hành và nhà nghiên cứu quan tâm có thể khai thác mô hình thuật toán này để nâng cao hiệu quả kinh tế và an toàn cho hệ thống điện ngay hôm nay.