Tổng quan nghiên cứu

Trong bối cảnh nền kinh tế toàn cầu tăng trưởng mạnh mẽ, nhu cầu tiêu thụ điện năng đã tăng đột biến khoảng 4% đến 6% mỗi năm, đặt ra thách thức nghiêm trọng cho công tác vận hành hệ thống truyền tải điện. Mục tiêu cốt lõi của các đơn vị vận hành hệ thống điện là tối ưu hóa chi phí sản xuất và phân phối điện năng, đồng thời phải đảm bảo tuyệt đối độ tin cậy và an ninh cung cấp điện. Bài toán phân bố công suất tối ưu có ràng buộc an ninh (Security Constrained Optimal Power Flow - SCOPF) là một bài toán quy hoạch phi tuyến quy mô lớn, phức tạp và có tính thực tiễn cao, đòi hỏi phải tìm ra điểm vận hành kinh tế nhất trong cả trạng thái vận hành cơ sở lẫn các trạng thái xảy ra sự cố đơn lẻ N-1.

Nghiên cứu của học viên Ngô Quang Sơn Hải dưới sự hướng dẫn khoa học của Phó Giáo sư Tiến sĩ Võ Ngọc Điều tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh (hoàn thành vào tháng 5 năm 2017) tập trung giải quyết bài toán SCOPF thông qua việc ứng dụng thuật toán tối ưu hóa bầy cá voi (Whale Optimization Algorithm - WOA). Mục tiêu trọng tâm của đề tài là xây dựng mô hình toán học hoàn chỉnh, lập trình giải thuật trên môi trường MATLAB và kiểm chứng hiệu năng trên ba hệ thống lưới điện chuẩn IEEE gồm 30 nút, 57 nút và 118 nút.

Nghiên cứu có ý nghĩa thực tiễn to lớn khi giúp giảm thiểu chi phí nhiên liệu máy phát từ 1,5% đến 3,2% so với các phương pháp giải thuật di truyền truyền thống, kiểm soát tổn thất công suất tác dụng trong khoảng 9,617 MW trên lưới 30 nút và duy trì điện áp tất cả các nút tải nghiêm ngặt trong giới hạn an toàn từ 0,95 đến 1,05 đơn vị tương đối (pu).

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 tích hợp hai nền tảng lý thuyết trọng tâm: lý thuyết tối ưu hóa phân bố công suất có ràng buộc an ninh và lý thuyết tối ưu hóa bầy đàn sinh học thông qua thuật toán WOA do Seyedali Mirjalili và Andrew Lewis công bố năm 2016.

Bài toán SCOPF được mô hình hóa toán học với hàm mục tiêu cực tiểu hóa tổng chi phí nhiên liệu của toàn bộ các tổ máy nhiệt điện trong hệ thống. Đường cong chi phí máy phát được thiết lập dưới dạng hàm bậc hai kết hợp thành phần hàm sin để phản ánh chính xác hiệu ứng điểm van công suất trong thực tế vận hành. Hệ thống ràng buộc bao gồm hai nhóm:

  • Nhóm ràng buộc đẳng thức: Các phương trình cân bằng công suất tác dụng và phản kháng tại từng nút mạng dựa trên định luật Kirchhoff, giải bằng phương pháp lặp dòng công suất Newton-Raphson.
  • Nhóm ràng buộc bất đẳng thức: Giới hạn công suất phát tác dụng và phản kháng của máy phát, ngưỡng điện áp nút phát (0,95 đến 1,1 pu), giới hạn điện áp nút tải (0,95 đến 1,05 pu), khả năng truyền tải dòng công suất biểu kiến cực đại trên các đường dây và phạm vi điều chỉnh nấc biến áp từ 0,9 đến 1,1 pu.

Thuật toán WOA mô phỏng hành vi xã hội và kỹ thuật săn mồi bằng lưới bong bóng của loài cá voi lưng gù. Khung giải thuật bao gồm ba cơ chế toán học then chốt: cơ chế bao vây con mồi, cơ chế tấn công bằng mạng lưới bong bóng xoắn ốc với xác suất chuyển đổi 50%, và cơ chế tìm kiếm con mồi ngẫu nhiên để mở rộng không gian tìm kiếm toàn cục.

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

Dữ liệu nghiên cứu được thu thập từ các hệ thống thử nghiệm tiêu chuẩn quốc tế của Viện Kỹ sư Điện và Điện tử (IEEE). Luận văn sử dụng cỡ mẫu gồm 3 mô hình lưới điện đại diện cho các quy mô khác nhau:

  • Mạng điện IEEE 30 nút: Gồm 6 máy phát, 41 đường dây truyền tải, 4 máy biến áp điều chỉnh nấc và 9 tụ bù công suất phản kháng, tổng công suất tải tác dụng là 283,4 MW và tải phản kháng là 126,2 MVar.
  • Mạng điện IEEE 57 nút: Đại diện cho lưới điện quy mô trung bình với 7 máy phát, 80 đường dây, 17 máy biến áp điều chỉnh nấc và 3 bộ tụ bù.
  • Mạng điện IEEE 118 nút: Đại diện cho hệ thống điện quy mô lớn với 54 nút phát điện, 186 đường dây, 9 máy biến áp và 14 trạm tụ bù chuyển mạch.

Phương pháp chọn mẫu có chủ đích được lựa chọn nhằm kiểm tra tính tổng quát, khả năng xử lý biến rời rạc và độ ổn định của thuật toán trên các không gian tìm kiếm từ 17 biến đến 131 biến điều khiển.

Phương pháp phân tích được thực hiện bằng cách lập trình trực tiếp thuật toán WOA trên phần mềm MATLAB kết hợp hàm giải trào lưu công suất Newton-Raphson. Các tham số thuật toán được cấu hình tối ưu: số lượng cá thể trong quần thể dao động từ 10 đến 40 phần tử, số vòng lặp cực đại từ 100 đến 250 thế hệ, và hệ số hàm phạt vi phạm ràng buộc được thiết lập trong khoảng từ 10.000 đến 1.000.000. Timeline nghiên cứu được thực hiện xuyên suốt trong 10 tháng, từ ngày 15 tháng 8 năm 2016 đến ngày 30 tháng 5 năm 2017.

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

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

Nghiên cứu đã chứng minh tính hiệu quả vượt trội của thuật toán WOA trong việc giải quyết bài toán SCOPF qua các kết quả thực nghiệm định lượng:

Thứ nhất, ở trạng thái vận hành cơ sở của mạng IEEE 30 nút, WOA đạt được tổng chi phí nhiên liệu tối ưu là 802,28 USD trên một giờ, với công suất phát tại nút cân bằng đạt 176,0021 MW và tổng tổn thất công suất truyền tải chỉ ở mức 9,617 MW. Mức chi phí này thấp hơn khoảng 1,8% so với phương pháp lập trình tiến hóa (EP) và thấp hơn 2,4% so với giải thuật di truyền tiêu chuẩn (GA).

Thứ hai, trong kịch bản sự cố N-1 nghiêm trọng khi ngắt đường dây truyền tải số 2-5 và đường dây số 1-3, thuật toán WOA tự động điều độ lại công suất các tổ máy phát (tăng công suất phát nút 2 lên 49,85 MW và nút 5 lên 21,34 MW) để giải tỏa nghẽn mạch. Tất cả các thông số điện áp nút tải đều được duy trì trong dải an toàn từ 0,962 pu đến 1,048 pu, không xuất hiện bất kỳ hiện tượng quá tải đường dây nào.

Thứ ba, đối với hệ thống lớn IEEE 57 nút và 118 nút, WOA thể hiện khả năng hội tụ nhanh chóng chỉ sau khoảng 120 đến 180 vòng lặp. Thời gian thực thi tính toán giảm từ 15% đến 20% so với phương pháp tối ưu hóa bầy đàn cổ điển (PSO) và phương pháp điểm nội cấu trúc (SEIP), trong khi vẫn kiểm soát chính xác 100% các biến điều khiển liên tục và rời rạc.

Thảo luận kết quả

Hiệu năng ấn tượng của WOA bắt nguồn từ sự cân bằng nhịp nhàng giữa hai giai đoạn thăm dò không gian tìm kiếm và khai thác cục bộ. Bằng cách giảm dần véc-tơ tham số A từ giá trị 2 về 0 trong suốt quá trình lặp, thuật toán dễ dàng thoát khỏi các cực trị địa phương – điểm yếu cố hữu của các phương pháp giải tích truyền thống như phương pháp mô hình mạng và điểm nội (NMIP).

Khi xem xét dữ liệu qua các biểu đồ hội tụ chi phí theo thời gian và bảng so sánh đa thông số, đường cong đặc tính của WOA cho thấy tốc độ giảm chi phí rất dốc trong 50 vòng lặp đầu tiên và ổn định vững chắc ở các thế hệ sau. Điều này minh chứng rằng mô hình săn mồi xoắn ốc logarit giúp các cá thể thu hẹp khoảng cách đến nghiệm tối ưu toàn cục nhanh hơn đáng kể so với bước nhảy ngẫu nhiên của thuật toán mô phỏng tôi luyện kim loại (SA). Kết quả nghiên cứu khẳng định WOA hoàn toàn đáp ứng tốt các yêu cầu khắt khe về kỹ thuật và kinh tế trong điều độ hệ thống điện hiện đại.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp cụ thể nhằm ứng dụng kết quả vào thực tế:

Thứ nhất, tích hợp trực tiếp module giải thuật WOA vào phần mềm quản lý năng lượng (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 hành động là tự động hóa tính toán SCOPF trực tuyến theo chu kỳ 15 phút, hướng tới mục tiêu cắt giảm từ 2% đến 3% tổng chi phí nhiên liệu phát điện hàng năm, thực hiện thử nghiệm trong vòng 12 tháng.

Thứ hai, mở rộng mô hình SCOPF có xét đến các nguồn năng lượng tái tạo biến đổi như điện mặt trời và điện gió. Tổng công ty Truyền tải điện Quốc gia cần chỉ đạo thực hiện đánh giá khả năng hấp thụ công suất và kiểm soát điện áp nút trong dải 0,95 đến 1,05 pu đối với các nút có tỷ trọng nguồn phân tán cao, hoàn thành trong giai đoạn 2024 - 2026.

Thứ ba, nâng cấp hệ thống phần cứng tính toán song song sử dụng công nghệ xử lý đồ họa (GPU Computing) kết hợp ngôn ngữ lập trình hiệu năng cao. Đơn vị kỹ thuật cần tối ưu hóa thời gian tính toán trào lưu công suất của lưới điện quy mô trên 500 nút xuống dưới 30 giây, triển khai trong vòng 6 tháng tới.

Thứ tư, chuẩn hóa quy trình cài đặt tham số giải thuật cho kỹ sư vận hành hệ thống điện. Các đơn vị nghiên cứu và vận hành cần phối hợp xây dựng bảng hướng dẫn tham số mẫu với quy mô quần thể từ 20 đến 50 cá thể và hệ số phạt 100.000 để đảm bảo thuật toán luôn hội tụ chính xác trong mọi kịch bản sự cố bất ngờ.

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

Tài liệu luận văn mang lại giá trị học thuật và ứng dụng sâu rộng cho 4 nhóm đối tượng chính:

Thứ nhất, các kỹ sư điều độ và chuyên gia vận hành hệ thống điện tại các tập đoàn điện lực. Luận văn cung cấp phương pháp luận chi tiết để tái điều độ công suất tổ máy, hiệu chỉnh nấc máy biến áp và điều khiển tụ bù nhằm xóa bỏ tình trạng nghẽn mạch truyền tải khi xảy ra sự cố đơn lẻ.

Thứ hai, các nhà nghiên cứu và giảng viên chuyên ngành Kỹ thuật điện, Tự động hóa và Toán ứng dụng. Đây là tài liệu tham khảo giá trị về cách thức mô hình hóa các bài toán kỹ thuật phi tuyến phức tạp và ứng dụng các thuật toán metaheuristic tiên tiến để giải quyết bài toán tối ưu đa chiều.

Thứ ba, học viên cao học và sinh viên đại học khối ngành năng lượng. Công trình là tài liệu mẫu mực về cấu trúc nghiên cứu, phương pháp lập trình mô phỏng trên nền tảng MATLAB và kỹ thuật xử lý dữ liệu lưới điện chuẩn IEEE.

Thứ tư, các nhà hoạch định chính sách và phát triển dự án hạ tầng năng lượng. Luận văn cung cấp cái nhìn định lượng về bài toán cân đối giữa chi phí đầu tư vận hành và biên độ an toàn lưới điện trong thị trường điện cạnh tranh.

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

Bài toán phân bố công suất tối ưu có ràng buộc an ninh khác gì so với bài toán tối ưu công suất truyền thống? Bài toán tối ưu truyền thống chỉ tìm điểm vận hành có chi phí thấp nhất ở trạng thái làm việc bình thường. Trong khi đó, bài toán có ràng buộc an ninh bắt buộc hệ thống phải thỏa mãn toàn bộ các giới hạn kỹ thuật về điện áp và dòng công suất ngay cả khi xảy ra các sự cố ngắt đường dây hoặc mất tổ máy phát.

Thuật toán bầy cá voi có ưu điểm gì vượt trội so với các thuật toán di truyền hoặc bầy đàn khác? Thuật toán bầy cá voi sở hữu cơ chế săn mồi xoắn ốc kết hợp tìm kiếm ngẫu nhiên độc đáo, chỉ cần tinh chỉnh hai tham số điều khiển chính. Nhờ đó, thuật toán có tốc độ hội tụ nhanh hơn khoảng 15%, không phụ thuộc vào đạo hàm hàm mục tiêu và hạn chế tối đa nguy cơ rơi vào cực tiểu cục bộ.

Làm thế nào để thuật toán xử lý được các biến điều khiển rời rạc như nấc máy biến áp và dung lượng tụ bù? Thuật toán sử dụng cơ chế làm tròn số học và hàm phạt trong hàm mục tiêu. Khi một biến rời rạc như nấc phân áp vượt ra ngoài bước chỉnh định từ 0,9 đến 1,1 pu, hệ số phạt có độ lớn từ 10.000 đến 1.000.000 sẽ tự động điều chỉnh hướng tìm kiếm của các cá thể về vùng nghiệm khả thi.

Luận văn đã kiểm chứng thuật toán trên những quy mô hệ thống điện nào? Đề tài đã thực hiện mô phỏng và kiểm chứng toàn diện trên 3 hệ thống chuẩn quốc tế của IEEE bao gồm: mạng IEEE 30 nút với 6 máy phát, mạng IEEE 57 nút với 7 máy phát và mạng IEEE 118 nút với 54 máy phát, bao quát từ lưới điện nhỏ đến lưới điện quy mô lớn.

Thời gian tính toán của giải thuật có đáp ứng được yêu cầu điều độ thời gian thực không? Với quy mô lưới điện IEEE 30 nút và 57 nút, thuật toán chỉ mất từ 100 đến 200 vòng lặp trên phần mềm MATLAB để đạt nghiệm tối ưu toàn cục. Khi được tích hợp trên các hệ thống máy chủ chuyên dụng của trung tâm điều khiển, thời gian tính toán hoàn toàn đáp ứng tốt chu kỳ điều độ thời gian thực.

Kết luận

  • Luận văn đã ứng dụng thành công thuật toán tối ưu hóa bầy cá voi để giải quyết trọn vẹn bài toán phân bố tối ưu công suất có ràng buộc an ninh trên các hệ thống điện phức tạp.
  • Kết quả tính toán trên mạng IEEE 30 nút đã tối ưu hóa chi phí nhiên liệu về mức 802,28 USD trên một giờ và kiểm soát tổn thất công suất ở mức 9,617 MW, vượt trội hơn các phương pháp GA và PSO.
  • Giải thuật đảm bảo duy trì ổn định điện áp tất cả các nút trong phạm vi an toàn từ 0,95 đến 1,05 pu trong mọi kịch bản sự cố mất điện đường dây đơn lẻ N-1.
  • Đề tài đã thiết lập thành công quy trình chuẩn hóa tham số với số cá thể từ 10 đến 40 và số vòng lặp từ 100 đến 250, đảm bảo tính hội tụ bền vững trên cả hệ thống quy mô lớn 118 nút.
  • Công trình mở ra hướng tiếp cận hiện đại, tin cậy cho công tác tự động hóa tính toán điều độ kinh tế và nâng cao an ninh năng lượng trong thị trường điện lực.

Trong giai đoạn 6 đến 12 tháng tới, hướng phát triển trọng tâm là mở rộng thuật toán cho lưới điện thực tế Việt Nam có tích hợp năng lượng tái tạo. Quý độc giả, kỹ sư và các nhà nghiên cứu quan tâm có thể tiếp tục khai thác mô hình toán học và mã nguồn thuật toán để nâng cao hiệu quả vận hành hệ thống điện bền vững.