Tổng quan nghiên cứu

Trong bối cảnh chi phí nhiên liệu hóa thạch không ngừng biến động và áp lực cắt giảm phát thải khí nhà kính ngày càng gia tăng, việc tối ưu hóa hiệu suất tại các nhà máy điện đóng vai trò then chốt đối với an ninh năng lượng. Các nhà máy nhiệt điện truyền thống thường chỉ đạt hiệu suất chuyển đổi nhiệt năng sang điện năng từ 50% đến 60%, dẫn đến một lượng lớn năng lượng bị thất thoát dưới dạng nhiệt thải ra môi trường. Ngược lại, công nghệ đồng phát nhiệt - điện (Combined Heat and Power - CHP) có thể tận dụng lượng nhiệt tổn thất này để đồng thời cung cấp điện năng và nhiệt lượng, giúp nâng tổng hiệu suất sử dụng năng lượng lên đến 80% và tiết kiệm từ 20% đến 30% lượng nhiên liệu tiêu thụ.

Tuy nhiên, việc tích hợp các tổ máy đồng phát vào hệ thống điện đặt ra bài toán Điều độ kinh tế hỗn hợp nhiệt - điện (Combined Heat and Power Economic Dispatch - CHPED) vô cùng phức tạp. Thách thức lớn nhất nằm ở tính chất phụ thuộc lẫn nhau phi tuyến tính giữa công suất phát điện (MW) và nhiệt lượng sản xuất (MWth) trong vùng vận hành khả thi (Feasible Operation Region).

Luận văn thạc sĩ chuyên ngành Thiết bị, Mạng và Nhà máy điện của tác giả Phạm Ngọc Hiếu, được thực hiện tại Trường Đại học Bách Khoa - Đại học Quốc gia TP.HCM dưới sự hướng dẫn khoa học của TS. Huỳnh Châu Duy (hoàn thành tháng 6/2014), đã tập trung giải quyết bài toán CHPED. Mục tiêu trọng tâm của nghiên cứu là xây dựng và ứng dụng giải thuật Tối ưu bầy đàn cải tiến (Improved Particle Swarm Optimization - IPSO) kết hợp kỹ thuật pseudo-gradient và hệ số co nhằm cực tiểu hóa tổng chi phí nhiên liệu vận hành ($/h), đồng thời đáp ứng nghiêm ngặt các ràng buộc về cân bằng phụ tải điện (từ 160 MW đến 500 MW), phụ tải nhiệt (từ 110 MWth đến 400 MWth) và giới hạn kỹ thuật của các tổ máy. Nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc giảm hàng chục nghìn USD chi phí nhiên liệu hàng năm và nâng cao độ ổn định cho lưới điện hiện đại.

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 được xây dựng dựa trên sự kết hợp chặt chẽ giữa mô hình toán học điều độ hệ thống năng lượng và lý thuyết trí tuệ bầy đàn hiện đại:

  • Mô hình toán học bài toán CHPED: Hệ thống bao gồm ba nhóm nguồn phát chính là tổ máy phát điện thuần túy, tổ máy đồng phát nhiệt - điện và lò hơi sản xuất nhiệt thuần túy. Hàm mục tiêu là cực tiểu hóa tổng chi phí nhiên liệu toàn hệ thống, trong đó chi phí phát điện và sản xuất nhiệt thuần túy được mô hình hóa bằng các hàm bậc hai, còn chi phí của máy đồng phát là hàm đa thức bậc hai hai biến kết hợp giữa công suất điện và sản lượng nhiệt. Bài toán chịu các ràng buộc đẳng thức nghiêm ngặt về cân bằng công suất điện, cân bằng nhiệt lượng và các ràng buộc bất đẳng thức về vùng vận hành khả thi nhiệt - điện bị chặn bởi các đường biên kỹ thuật phức tạp.
  • Lý thuyết tối ưu bầy đàn (Particle Swarm Optimization - PSO): Dựa trên thuật toán tiến hóa do Kennedy và Eberhart phát triển năm 1995, mô phỏng hành vi xã hội của đàn chim tìm kiếm thức ăn. Mỗi cá thể đại diện cho một phương án phân bổ công suất, liên tục điều chỉnh vận tốc và vị trí dựa theo kinh nghiệm tốt nhất của bản thân (pbest) và kinh nghiệm tốt nhất của toàn bầy (gbest).
  • Kỹ thuật Pseudo-gradient và Hệ số co: Đối với các bài toán tối ưu phi lồi và không khả vi như CHPED, gradient giải tích thông thường không áp dụng được. Thuật toán IPSO sử dụng pseudo-gradient để xác định hướng dốc dương nhằm dẫn dắt các phần tử di chuyển nhanh chóng về phía nghiệm tối ưu. Đồng thời, hệ số co Clerc với thông số tổng hợp lớn hơn 4 (tiêu chuẩn chọn giá trị 4,1) được tích hợp vào phương trình vận tốc để kiểm soát dao động và đẩy nhanh tốc độ hội tụ toàn cục.

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

  • Nguồn dữ liệu và cấu hình hệ thống: Nghiên cứu sử dụng bộ dữ liệu kiểm chuẩn quốc tế gồm 4 hệ thống thử nghiệm: Hệ thống 4 máy phát (1 máy phát điện, 2 máy đồng phát, 1 máy phát nhiệt); Hệ thống 5 máy phát tích hợp thêm tuabin hơi cố định công suất 80 MW điện và tiêu thụ 120 MWth nhiệt; Hệ thống 5 máy phát mở rộng (1 máy phát điện, 3 máy đồng phát, 1 máy phát nhiệt); và Hệ thống 6 máy phát quy mô lớn (1 máy phát điện, 4 máy đồng phát, 1 máy phát nhiệt).
  • Cỡ mẫu và phương pháp thực nghiệm: Để đảm bảo tính khách quan và loại trừ tính bất định của thuật toán ngẫu nhiên, mỗi kịch bản phụ tải được thực hiện mô phỏng độc lập với cỡ mẫu 50 lần chạy (50 Monte Carlo runs), số vòng lặp tối đa thiết lập là 1.000 thế hệ trên phần mềm tính toán Matlab. Nghiên cứu được triển khai xuyên suốt trong 5 tháng thực nghiệm (từ tháng 1/2014 đến tháng 6/2014).
  • Lý do lựa chọn phương pháp phân tích: Việc phân tích thống kê đa chỉ số (chi phí nhỏ nhất, chi phí lớn nhất, chi phí trung bình, độ lệch tiêu chuẩn và thời gian xử lý) cho phép đánh giá toàn diện độ chính xác và tính ổn định của IPSO. Kết quả được so sánh đối chuẩn trực tiếp với nhiều thuật toán kinh điển như Nhân tử Lagrange (LR), Giải thuật Di truyền (GA, IGA-MU), Lập trình tiến hóa (EP), Mạng Neural Hopfield - Lagrange tăng cường (ALHN), Thuật toán Tìm kiếm hòa âm (HSA), Thuật toán Tìm kiếm bầy kiến cải tiến (LACSA) và Thuật toán Tìm kiếm chim cúc cu (CS).

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

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

  • Tiết kiệm chi phí nhiên liệu vượt trội: Tại hệ thống 4 máy phát với phụ tải 200 MW điện và 115 MWth nhiệt, IPSO đạt mức chi phí nhiên liệu tối ưu là 9.257,075 $/h. Mức chi phí này thấp hơn đáng kể so với phương pháp Nhân tử Lagrange (9.257,357 $/h), phương pháp HSA (9.257,102 $/h) và phương pháp EP (9.257,080 $/h).
  • Tốc độ hội tụ siêu nhanh: Nhờ bổ sung vector dẫn đường pseudo-gradient, IPSO rút ngắn từ 60% đến 90% số vòng lặp cần thiết để đạt điểm tối ưu. Cụ thể, ở kịch bản phụ tải PD = 160 MW và HD = 220 MWth của hệ thống 5 máy phát, IPSO đạt trạng thái hội tụ chỉ sau 8 vòng lặp. Tại hệ thống 6 máy phát với phụ tải lớn PD = 450 MW và HD = 400 MWth, thuật toán tìm ra điểm cực tiểu toàn cục chỉ sau 6 thế hệ lặp.
  • Độ ổn định cao qua nhiều kịch bản phụ tải: Qua 50 lần chạy độc lập ở hệ thống 4 máy phát phụ tải 175 MW điện và 110 MWth nhiệt, IPSO đạt chi phí tối ưu 8.555,274 $/h sau 80 vòng lặp. Tương tự, ở kịch bản PD = 225 MW và HD = 125 MWth, IPSO hội tụ vững chắc sau 40 vòng lặp.
  • Hiệu năng vượt trội trên hệ thống quy mô lớn: Đối với hệ thống 6 máy phát với phụ tải 500 MW điện và 300 MWth nhiệt, IPSO đạt mức chi phí tối ưu 22.009,369 $/h sau 120 vòng lặp, vượt trội hoàn toàn so với thuật toán PSO nguyên bản đạt mức 22.021,374 $/h (tiết kiệm 12,005 $/h, tương đương giảm hơn 105.000 USD chi phí vận hành mỗi năm nếu vận hành liên tục 8.760 giờ).

Thảo luận kết quả

Nguyên nhân chính giúp IPSO vượt qua các thuật toán heuristic khác là nhờ cơ chế pseudo-gradient đã khắc phục triệt để hiện tượng hội tụ sớm tại các cực trị địa phương – nhược điểm cố hữu của PSO chuẩn và GA. Trong các hệ thống năng lượng phức tạp, các đường biên của vùng vận hành khả thi nhiệt - điện tạo ra các miền không gian tìm kiếm không lồi. Pseudo-gradient đóng vai trò như một lực đẩy định hướng, kết hợp cùng hệ số co điều tiết vận tốc hạt, giúp bầy đàn tập trung khai thác sâu vào vùng nghiệm khả thi có chi phí thấp nhất mà không bị phân tán lãng phí.

Dữ liệu thực nghiệm được trực quan hóa rõ nét qua các biểu đồ quá trình hội tụ theo thang đo logarithmic và bảng so sánh 5 thông số thống kê. Các biểu đồ đặc tính cho thấy đường cong chi phí của IPSO giảm dốc đứng ngay từ 10 đến 50 vòng lặp đầu tiên, phản ánh khả năng xử lý nhanh chóng các hàm phạt vi phạm ràng buộc công suất và nhiệt lượng. Khi quy mô hệ thống tăng từ 4 máy lên 6 máy phát, độ lệch tiêu chuẩn của IPSO vẫn duy trì ở mức thấp, khẳng định tính khả thi cao khi ứng dụng vào các hệ thống điều độ thời gian thực.

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

  • Tích hợp giải thuật IPSO vào hệ thống phần mềm SCADA/EMS: Các trung tâm điều độ hệ thống điện và doanh nghiệp quản lý khu công nghiệp cần sớm triển khai thuật toán IPSO vào module tự động điều độ kinh tế thời gian thực. Mục tiêu là rút ngắn thời gian tính toán phân bổ công suất xuống dưới 1 giây và cắt giảm từ 0,1% đến 0,5% chi phí nhiên liệu mỗi ngày. Lộ trình thực hiện khuyến nghị trong vòng 6 đến 12 tháng do đội ngũ kỹ sư điều độ và chuyên gia CNTT đảm trách.
  • Mở rộng bài toán tính toán tổn thất lưới điện và mạng nhiệt: Các viện nghiên cứu và phòng kỹ thuật thuộc các tổng công ty điện lực cần nâng cấp mô hình toán học bằng cách tích hợp thêm ma trận hệ số tổn thất công suất B-matrix trên đường dây truyền tải và tổn thất nhiệt trên đường ống dẫn hơi. Mục tiêu nâng cao độ chính xác điều độ lên trên 98,5% với thời gian nghiên cứu hoàn thiện từ 12 đến 18 tháng.
  • Chuẩn hóa quy trình xác định vùng vận hành khả thi (FOR): Đội ngũ kỹ sư tại các nhà máy đồng phát cần tiến hành đo đạc thực nghiệm định kỳ 3 tháng/lần để cập nhật chính xác đường đặc tính nhiệt - điện của từng tổ máy. Mục tiêu loại bỏ hoàn toàn nguy cơ vận hành ngoài vùng an toàn kỹ thuật, giúp tăng tuổi thọ thiết bị lò hơi và tuabin thêm 5% đến 10%.
  • Phát triển mô hình điều độ đa mục tiêu kết hợp giảm phát thải: Cơ quan quản lý năng lượng và Bộ Công Thương cần ban hành các hướng dẫn kỹ thuật thúc đẩy điều độ hỗn hợp nhiệt - điện gắn liền với chỉ tiêu giảm phát thải CO2, SOx, NOx. Mục tiêu hướng đến việc cắt giảm từ 5% đến 12% lượng khí phát thải nhà kính toàn ngành điện vào năm 2026 thông qua các cơ chế định giá carbon và tối ưu hóa đa mục tiêu.

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

  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Điện, Tự động hóa: Luận văn cung cấp nền tảng toán học mẫu mực về tối ưu hóa phi tuyến, phương pháp thiết lập hàm phạt và kỹ thuật lập trình giải thuật metaheuristic hiện đại trên nền tảng Matlab.
  • Kỹ sư điều độ tại các Trung tâm Điều độ Hệ thống điện: Cung cấp phương pháp luận và công cụ tính toán phân bổ tối ưu công suất tức thời giữa các nguồn điện truyền thống và nguồn đồng phát công nghiệp trên lưới điện phân phối.
  • Ban quản lý kỹ thuật tại các nhà máy nhiệt điện đồng phát, khu công nghiệp: Cung cấp giải pháp kỹ thuật cụ thể để vận hành tối ưu các tổ máy phát điện và lò hơi phụ trợ, giúp tiết kiệm hàng trăm nghìn USD chi phí than, dầu và khí đốt hàng năm.
  • Chuyên gia phân tích chính sách và nhà phát triển giải pháp Lưới điện thông minh (Smart Grid): Nguồn tài liệu tham khảo giá trị để xây dựng các thuật toán điều khiển tự động cho các hệ thống vi lưới (Microgrid) tích hợp đa nguồn năng lượng điện - nhiệt.

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

  • Bài toán điều độ kinh tế hỗn hợp nhiệt - điện (CHPED) khác biệt thế nào so với điều độ kinh tế truyền thống (ED)? Bài toán ED truyền thống chỉ tập trung tối ưu hóa chi phí phát điện năng thuần túy. Ngược lại, bài toán CHPED phải đồng thời thỏa mãn cả hai nhu cầu phụ tải điện và nhiệt. Sự phụ thuộc ràng buộc phi tuyến giữa công suất điện và sản lượng nhiệt trong vùng vận hành khả thi khiến không gian tìm kiếm của CHPED trở nên phức tạp và phi lồi hơn rất nhiều.

  • Cải tiến then chốt của thuật toán IPSO so với PSO nguyên bản là gì? IPSO bổ sung kỹ thuật pseudo-gradient và hệ số co Clerc vào phương trình cập nhật vận tốc hạt. Pseudo-gradient tạo vector dẫn hướng có chủ đích giúp các hạt di chuyển thẳng tới vùng nghiệm tốt mà không đòi hỏi hàm mục tiêu phải khả vi, giúp giảm số vòng lặp từ hàng trăm xuống chỉ còn 6 đến 40 vòng lặp.

  • Thuật toán IPSO đã được kiểm chứng trên những hệ thống thử nghiệm nào? Luận văn đã kiểm chứng thuật toán trên 4 mô hình hệ thống chuẩn gồm hệ thống 4 máy phát, 5 máy phát có tuabin hơi công suất cố định 80 MW, 5 máy phát có 3 tổ máy CHP, và hệ thống 6 máy phát quy mô lớn. Mỗi hệ thống đều được chạy thử nghiệm 50 lần độc lập dưới nhiều mức phụ tải điện từ 160 MW đến 500 MW và phụ tải nhiệt từ 110 MWth đến 400 MWth.

  • Hiệu quả kinh tế mà IPSO mang lại so với các phương pháp khác cụ thể ra sao? IPSO luôn đem lại chi phí nhiên liệu thấp hơn các giải thuật đối chứng. Điển hình tại hệ thống 6 máy phát tải 500 MW điện và 300 MWth nhiệt, IPSO đạt chi phí 22.009,369 $/h, tiết kiệm hơn 12 $/h so với PSO chuẩn, giúp tiết kiệm hàng trăm nghìn USD chi phí nhiên liệu mỗi năm khi vận hành liên tục.

  • Thuật toán IPSO trong luận văn được lập trình trên công cụ nào? Toàn bộ thuật toán IPSO, cấu trúc dữ liệu và các hàm ràng buộc được lập trình và mô phỏng hoàn chỉnh trên môi trường Matlab. Kết quả thực nghiệm cho thấy thời gian xử lý trung bình rất nhanh, đáp ứng tốt yêu cầu tính toán của các bài toán điều độ thực tế.

Kết luận

  • Giải quyết thành công bài toán điều độ kinh tế hỗn hợp nhiệt - điện (CHPED) với cấu trúc hàm mục tiêu phi tuyến và các ràng buộc biên nhiệt - điện phức tạp.
  • Cải tiến và ứng dụng xuất sắc thuật toán IPSO tích hợp pseudo-gradient và hệ số co, nâng cao vượt bậc tốc độ hội tụ và khả năng thoát bẫy cực trị cục bộ.
  • Kiểm chứng toàn diện trên 4 hệ thống thử nghiệm với cỡ mẫu 50 lần chạy độc lập, chứng minh chi phí nhiên liệu luôn đạt mức tối ưu thấp nhất so với các thuật toán LR, GA, HSA, EP và PSO.
  • Rút ngắn thời gian tính toán vượt trội khi thuật toán đạt độ hội tụ chỉ sau 6 đến 120 vòng lặp, đáp ứng hoàn hảo tiêu chuẩn vận hành điều độ thời gian thực.
  • Đóng góp giải pháp khoa học và thực tiễn giá trị cao, mở ra hướng đi bền vững giúp các doanh nghiệp năng lượng tiết kiệm hàng trăm nghìn USD chi phí vận hành hàng năm.

Luận văn thạc sĩ của tác giả Phạm Ngọc Hiếu là công trình nghiên cứu công phu, tạo tiền đề quan trọng cho việc mở rộng mô hình tính toán tổn thất lưới điện và tối ưu hóa phát thải môi trường trong giai đoạn 2026-2027. Hãy tham khảo và áp dụng ngay các giải pháp từ luận văn này để nâng cao hiệu quả kinh tế và độ tin cậy cho hệ thống năng lượng của bạn.