Tổng quan nghiên cứu

Trong hệ thống năng lượng hiện đại, điện năng là dạng năng lượng đặc thù gần như không thể lưu trữ ở quy mô công nghiệp trong thời gian thực. Bài toán vận hành hệ thống điện đòi hỏi phải cân bằng liên tục giữa việc tối thiểu hóa chi phí nhiên liệu phát điện và đảm bảo độ tin cậy cung cấp điện tuyệt đối. Phân bố công suất tối ưu có xét đến ràng buộc an ninh (Security Constrained Optimal Power Flow - SCOPF) là bài toán phi tuyến phức tạp bậc cao, kết hợp giữa mục tiêu kinh tế của bài toán trào lưu công suất tối ưu (OPF) và yêu cầu kỹ thuật nghiêm ngặt nhằm duy trì trạng thái ổn định khi xảy ra sự cố theo tiêu chuẩn an ninh N-1.

Vấn đề nghiên cứu cốt lõi được đặt ra là sự bùng nổ về khối lượng tính toán khi quy mô lưới điện mở rộng. Khi khảo sát toàn diện 177 trường hợp sự cố đứt dây đơn trên mạng chuẩn IEEE 118 nút, thời gian giải bài toán SCOPF truyền thống có thể kéo dài lên tới 261 giờ (khoảng 11 ngày) do số lượng ràng buộc tăng theo hàm bậc hai của số trạng thái sự cố. Hơn nữa, các giải thuật tối ưu cổ điển dựa trên đạo hàm giải tích thường gặp bế tắc trước các hàm mục tiêu không khả vi hoặc dễ mắc kẹt tại các điểm cực trị địa phương.

Nghiên cứu tập trung phát triển giải thuật bầy đàn cải tiến (IPSO) kết hợp kỹ thuật Gradient giả (Pseudo-Gradient), hệ số co (Constriction Factor) và trọng số quán tính biến thiên theo thời gian (PSO-TVIW). Mục tiêu cụ thể là tối ưu hóa biểu đồ phát công suất của 54 tổ máy phát trên mạng chuẩn IEEE 118 nút, đồng thời kiểm chứng tính thích nghi trên mạng IEEE 6 nút và New England 39 nút. Kết quả đạt được giúp giải quyết triệt để bài toán an ninh tĩnh toàn diện, rút ngắn thời gian xử lý xuống còn khoảng 4 giờ, tạo cơ sở khoa học tin cậy cho công tác điều độ hệ thống điện thời gian thực.

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 xây dựng trên nền tảng lý thuyết phân bố công suất tối ưu (OPF) và lý thuyết an ninh hệ thống điện theo tiêu chuẩn N-1 của Carpentier. Mô hình toán học tổng quát hướng tới cực tiểu hóa tổng chi phí nhiên liệu của toàn bộ các máy phát trong hệ thống theo hàm chi phí bậc hai:

$$F(x, u) = \sum_{n=1}^{NG} (a_n + b_n P_n + c_n P_n^2)$$

Trong đó, $NG$ là tổng số tổ máy phát; $a_n, b_n, c_n$ là các hệ số chi phí nhiên liệu đặc trưng.

Hệ thống biến số bao gồm:

  1. Vector biến điều khiển $u$: Điện áp tại các nút máy phát ($V_g$), vị trí nấc phân áp của máy biến áp ($T_k$) và dung lượng bù công suất phản kháng của tụ bù tĩnh ($Q_c$).
  2. Vector biến phụ thuộc $x$: Góc pha điện áp, độ lệch điện áp tại các nút tải ($V_l$), công suất phản kháng ngõ ra của máy phát ($Q_g$) và dòng công suất truyền tải trên các nhánh ($S_l$).
  3. Hàm phạt đa thành phần: Được tích hợp vào hàm thích nghi (Fitness) với các trọng số phạt tương ứng gồm $K_p$ (công suất nút cân bằng), $K_q$ (công suất phản kháng), $K_v$ (điện áp tải) và $K_s$ (công suất biểu kiến nhánh truyền tải) cho cả trạng thái cơ sở lẫn trạng thái sự cố.

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

Nghiên cứu sử dụng tập dữ liệu kiểm chuẩn quốc tế gồm 3 mô hình mạng điện tiêu chuẩn:

  • Mạng IEEE 6 nút: Gồm 3 máy phát, 11 nhánh truyền tải, phục vụ kiểm chứng thuật toán ở chế độ cơ sở và sự cố đứt đường dây 2-3.
  • Mạng New England 39 nút: Gồm 10 máy phát, 46 đường dây, đại diện cho hệ thống lưới điện truyền tải khu vực.
  • Mạng IEEE 118 nút: Cỡ mẫu lớn gồm 54 máy phát, 186 nhánh truyền tải và toàn bộ 177 kịch bản sự cố đứt dây đơn N-1.

Phương pháp chọn mẫu thực nghiệm bao phủ toàn bộ 100% các phần tử đường dây thay vì chỉ khảo sát một vài trường hợp nghi ngờ đơn lẻ. Nghiên cứu lựa chọn phương pháp phân tích tối ưu meta-heuristic cải tiến kết hợp giải thuật phân rã DSCOPF (Decomposed SCOPF).

Lý do lựa chọn giải thuật IPSO đề xuất xuất phát từ nhu cầu khắc phục nhược điểm hội tụ chậm của PSO nguyên bản. Việc bổ sung Gradient giả cung cấp vector định hướng tức thời cho các hạt mà không đòi hỏi hàm mục tiêu phải khả vi liên tục. Đồng thời, việc áp dụng hệ số gia tốc nhận thức và xã hội $c_1 = c_2 = 2.05$, kết hợp hệ số giới hạn vận tốc $R$ trong khoảng từ 0.1 đến 0.5 giúp kiểm soát chính xác quỹ đạo chuyển động của các hạt trong không gian tìm kiếm đa chiều.

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 và tính toán trên các mạng điện tiêu chuẩn đã mang lại các kết quả vượt trội:

  1. Đột phá về tốc độ hội tụ và thời gian xử lý: Khi áp dụng thuật toán IPSO kết hợp phân rã DSCOPF trên mạng IEEE 118 nút với 177 kịch bản sự cố, thời gian tính toán giảm mạnh từ 261 giờ xuống chỉ còn khoảng 4 giờ, tương đương mức giảm hơn 98% thời gian thực thi so với phương pháp tính toán vét cạn thông thường.
  2. Tối ưu hóa chi phí nhiên liệu máy phát: Giải thuật IPSO đạt giá trị hàm chi phí thấp hơn từ 1,5% đến 3,8% so với các giải thuật di truyền tiêu chuẩn (GA) và giải thuật tiến hóa nâng cao (EPGA II). Điển hình tại mạng IEEE 118 nút, thuật toán ưu tiên huy động công suất tối đa của máy phát số 26 (đạt 414 MW với chi phí vận hành kinh tế 133,5 USD/h) thay vì duy trì công suất cao ở máy phát số 1 (75,14 USD/h tại mức phát tối thiểu 10 MW).
  3. Duy trì tuyệt đối độ an toàn lưới điện: Toàn bộ điện áp tại các nút tải được kiểm soát chặt chẽ trong dải cho phép từ 0,95 p.u đến 1,05 p.u. Chỉ số an toàn quá tải đường dây (SI) tại tất cả 177 trường hợp sự cố đứt dây đơn đều không xuất hiện hiện tượng vượt ngưỡng nhiệt truyền tải danh định.
  4. Tối ưu hóa chỉ số ưu tiên nguồn phát (Generator Index - GI): Phân bố chỉ số $GI_i$ của 54 tổ máy phát tiệm cận về giá trị 1 đối với các tổ máy có chi phí rẻ, chứng minh tính định hướng chính xác của vector Gradient giả.

Thảo luận kết quả

Hiệu năng vượt trội của giải thuật IPSO bắt nguồn từ sự phối hợp chặt chẽ giữa Gradient giả và hệ số co Clerc. Trong thuật toán PSO truyền thống, các hạt dễ bị mất phương hướng khi di chuyển trong các vùng phẳng của hàm mục tiêu hoặc bị dao động quá mức quanh cực trị địa phương. Cơ chế Gradient giả tạo ra một lực đẩy định hướng có điều kiện: khi hàm fitness cải thiện, hướng di chuyển được xác lập dương; khi hàm fitness suy giảm, hướng tìm kiếm lập tức bị triệt tiêu để chuyển hướng lân cận.

So với phương pháp quy hoạch phi tuyến tuần tự (SQP) vốn đòi hỏi ma trận đạo hàm liên tục, IPSO xử lý linh hoạt các ràng buộc phi tuyến và điểm làm việc không liên tục của máy biến áp. So với giải thuật tìm kiếm bầy vi khuẩn cải tiến (IBF) và thuật toán tìm kiếm hấp dẫn (GSA), IPSO có cấu trúc tham số tinh gọn hơn, giúp giảm áp lực bộ nhớ khi thực hiện lặp đa tầng.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị hội tụ của hàm Fitness theo từng thế hệ lặp, bảng đối chiếu phân bố công suất thực $P_g$ giữa 3 phương pháp (IPSO, GA, EPGA II), và biểu đồ cột so sánh chỉ số $GI$ giữa 54 tổ máy phát. Sự tương thích cao giữa các thông số chứng minh giải thuật đã hội tụ về nghiệm tối ưu toàn cục mà không vi phạm bất kỳ ràng buộc kỹ thuật nào.

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

Dựa trên kết quả đạt được, luận văn đưa ra 4 khuyến nghị kỹ thuật và lộ trình ứng dụng thực tiễn:

  1. Tích hợp mô-đun IPSO vào hệ thống phần mềm quản lý năng lượng (EMS/SCADA): Trung tâm Điều độ Hệ thống điện Quốc gia và các cấp điều độ miền nên ứng dụng thuật toán IPSO vào quy trình lập phương thức vận hành ngày tới và điều độ thời gian thực chu kỳ 15 phút, nhằm cắt giảm từ 2% đến 5% chi phí mua điện toàn hệ thống.
  2. Chuẩn hóa quy trình tự động quét an ninh tĩnh N-1: Tập đoàn Điện lực cần triển khai quy trình đánh giá tự động 100% kịch bản sự cố trên lưới truyền tải 220kV và 500kV dựa trên nền tảng giải thuật phân rã DSCOPF trong lộ trình chuyển đổi số giai đoạn 2026 - 2030.
  3. Mở rộng mô hình tính toán tích hợp năng lượng tái tạo: Các viện nghiên cứu chuyên ngành cần chủ trì nâng cấp thuật toán, đưa thêm các ràng buộc ngẫu nhiên về độ bất định công suất của nguồn điện gió và điện mặt trời (biến động 15% - 30%) vào bài toán SCOPF trong giai đoạn 2027 - 2028.
  4. Triển khai bài toán điều độ công suất phản kháng có ràng buộc an ninh (SCORPD): Các công ty truyền tải điện cần ứng dụng thuật toán kết hợp tối ưu hóa các thiết bị bù tĩnh SVC và nấc phân áp máy biến áp nhằm duy trì điện áp nút trong phạm vi hẹp từ 0,98 p.u đến 1,02 p.u, thực hiện hoàn thành trong vòng 12 đến 24 tháng tới.

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

Nội dung luận văn mang giá trị học thuật và ứng dụng cao cho 4 nhóm đối tượng trọng tâm:

  1. Kỹ sư điều độ và vận hành hệ thống điện: Nắm vững nguyên lý phân rã bài toán an ninh phức tạp, sử dụng công cụ tính toán để ra quyết định điều chuyển tải và huy động nguồn phát nhanh chóng khi xuất hiện nguy cơ sự cố.
  2. Chuyên gia quy hoạch và phân tích lưới điện: Tham khảo phương pháp đánh giá an toàn qua chỉ số quá tải đường dây (SI) và chỉ số ưu tiên nguồn ($GI$) để tối ưu hóa kế hoạch nâng cấp công suất nhánh truyền tải trên lưới điện thực tế.
  3. Học viên cao học, nghiên cứu sinh và giảng viên ngành Kỹ thuật điện: Sử dụng công trình như một tài liệu tham khảo chuẩn mực về mô hình toán học SCOPF, kỹ thuật lai ghép meta-heuristic và phương pháp xử lý hàm không khả vi.
  4. Đơn vị phát triển phần mềm mô phỏng hệ thống điện: Khai thác thuật toán IPSO với các tham số $c_1=c_2=2.05$ và hệ số co $C$ để xây dựng các mô-đun tối ưu hóa thương mại có tốc độ hội tụ cao.

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

Bài toán SCOPF khác biệt gì so với bài toán OPF truyền thống trong hệ thống điện? Bài toán OPF thông thường chỉ tìm điểm vận hành có chi phí thấp nhất ở trạng thái cơ sở. Trong khi đó, SCOPF mở rộng thêm toàn bộ các ràng buộc cân bằng và bất đẳng thức khi xảy ra sự cố N-1. Điều này đảm bảo hệ thống không bị sụp đổ điện áp hoặc quá tải đường dây khi mất đi bất kỳ một phần tử nào.

Tại sao giải thuật PSO cơ bản cần phải cải tiến bằng Pseudo-Gradient và hệ số co? PSO cơ bản dễ bị kẹt tại các cực tiểu địa phương và có tốc độ hội tụ chậm khi số chiều bài toán lớn. Kỹ thuật Gradient giả giúp xác định chính xác vector dẫn hướng cho hàm không khả vi, kết hợp hệ số co với $c_1 = c_2 = 2.05$ giúp kiểm soát vận tốc hạt, đảm bảo tìm kiếm tối ưu toàn cục.

Phương pháp phân rã DSCOPF giúp cải thiện hiệu năng tính toán như thế nào? Thay vì giải đồng thời 177 bài toán sự cố gây quá tải bộ nhớ, DSCOPF tách hệ thống thành bài toán cơ sở và các bài toán con sự cố. Thuật toán chỉ tái cấu trúc các ràng buộc cho 2 đến 3 sự cố xung yếu nhất thông qua vết cắt Benders, giúp rút ngắn thời gian tính toán từ 261 giờ xuống còn 4 giờ.

Chỉ số ưu tiên nguồn phát ($GI$) mang ý nghĩa thực tiễn gì trong vận hành kinh tế? Chỉ số $GI_i$ nằm trong đoạn từ 0 đến 1, định lượng mức độ ưu tiên huy động công suất của từng tổ máy phát. Thuật toán tối ưu sẽ tự động đẩy $GI$ của các nhà máy có suất tiêu hao nhiệt thấp tiến sát về 1 (phát cực đại) và giữ $GI$ của nguồn phát giá cao ở mức 0 (phát cực tiểu).

Thuật toán IPSO đã được kiểm chứng trên những hệ thống lưới điện tiêu chuẩn nào? Thuật toán đã được mô phỏng kiểm chứng thành công trên 3 hệ thống chuẩn quốc tế: mạng IEEE 6 nút (3 máy phát), mạng New England 39 nút (10 máy phát) và mạng IEEE 118 nút (54 máy phát, 186 đường dây, 177 kịch bản sự cố N-1), cho kết quả vượt trội so với GA và EPGA II.

Kết luận

  • Luận văn đã thiết lập thành công mô hình toán học toàn diện cho bài toán phân bố công suất tối ưu có xét đến ràng buộc an ninh N-1 đầy đủ.
  • Đề xuất giải thuật cải tiến IPSO kết hợp Gradient giả, hệ số co và giới hạn vận tốc $R \in [0.1, 0.5]$, xử lý hiệu quả các hàm mục tiêu phi tuyến phức tạp.
  • Ứng dụng giải thuật phân rã DSCOPF giúp giảm hơn 98% thời gian tính toán trên mạng chuẩn IEEE 118 nút, rút ngắn tiến trình giải từ 261 giờ xuống 4 giờ.
  • Kiểm soát ổn định toàn diện 177 trường hợp sự cố đứt dây đơn, tối ưu hóa chỉ số ưu tiên nguồn phát $GI$ và tiết kiệm chi phí nhiên liệu vượt trội so với GA và EPGA II.
  • Xác lập cơ sở khoa học vững chắc để tích hợp công cụ tính toán an ninh vào hệ thống điều độ thực tế giai đoạn 2026 - 2030.

Để nâng cao hiệu quả vận hành kinh tế và bảo đảm an ninh hệ thống điện, các kỹ sư và nhà nghiên cứu hãy tham khảo chi tiết toàn văn luận văn và ứng dụng giải thuật IPSO vào các bài toán tối ưu hóa nguồn lưới điện thực tế ngay hôm nay.