Tổng quan nghiên cứu

Hệ thống điện đóng vai trò huyết mạch trong sự phát triển kinh tế và an ninh năng lượng quốc gia. Tại Việt Nam, hệ thống truyền tải điện siêu cao áp 500kV Bắc - Nam vận hành từ năm 1994 với chiều dài trên 1.487 km đã kết nối liên kết lưới điện toàn quốc, phục vụ tốc độ tăng trưởng phụ tải công nghiệp từ 10% đến 12% mỗi năm. Sự liên kết quy mô lớn cùng việc hòa lưới các nguồn điện trọng điểm như Thủy điện Sơn La công suất 2.400 MW và Trung tâm Nhiệt điện Phú Mỹ công suất 3.900 MW làm gia tăng mức độ phức tạp trong vận hành. Bài toán cấp thiết đặt ra là duy trì sự vận hành đồng bộ và ổn định hệ thống điện khi xuất hiện các sự cố ngắn mạch nghiêm trọng.

Nghiên cứu tập trung giải quyết bài toán đánh giá ổn định động cho hệ thống điện nhiều máy phát, vượt qua giới hạn của phương pháp truyền thống vốn chỉ quy đổi về hệ một máy phát nối thanh cái vô cùng lớn. Mục tiêu cốt lõi của công trình là xây dựng thuật toán phân tích trào lưu công suất và thiết lập công cụ mô phỏng quá trình quá độ điện cơ trên phần mềm Matlab m-file. Nghiên cứu thực hiện khảo sát chi tiết trên hệ thống chuẩn 9 nút - 3 máy phát nhằm xác định chính xác thời gian cắt sự cố tới hạn, kiểm soát độ lệch góc pha rôto và duy trì giới hạn an toàn trong khoảng thời gian từ 0,15 đến 0,25 giây sau sự cố. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc nâng cao độ tin cậy cung cấp điện, giảm thiểu rủi ro rã lưới trên diện rộng và hỗ trợ đắc lực cho công tác lập phương thức vận hành hệ thống điện quốc gia.

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 nền tảng lý thuyết quá trình quá độ điện cơ trong hệ tọa độ quay do R. Park đề xuất năm 1928 và A. Goriev phát triển hoàn thiện trong giai đoạn 1930 - 1935. Hệ phương trình Park - Goriev mô tả hành vi phi tuyến của các máy điện đồng bộ làm việc song song thông qua phương trình chuyển động rôto (phương trình Swing).

Lý thuyết ổn định hệ thống điện tiếp tục vận dụng hai phương pháp cốt lõi của nhà toán học A. Lyapunov:

  • Phương pháp xấp xỉ bậc nhất (phương pháp tuyến tính hóa) sử dụng tiêu chuẩn đại số Routh - Hurwitz và tiêu chuẩn tần số Mikhailov để đánh giá ổn định tĩnh với các dao động bé quanh điểm cân bằng.
  • Phương pháp trực tiếp thứ hai thông qua việc xây dựng hàm năng lượng Lyapunov toàn phần để xác định miền ổn định động khi hệ thống chịu kích động lớn.

Nghiên cứu phân định rõ bốn khái niệm bản chất: chế độ xác lập bình thường, chế độ quá độ sự cố, ổn định tĩnh (duy trì điểm cân bằng khi góc công suất nhỏ hơn 90 độ) và ổn định động (khả năng khôi phục đồng bộ sau sự cố ngắn mạch lớn). Thời gian cắt tới hạn được định nghĩa là khoảng thời gian tối đa cho phép các thiết bị bảo vệ rơle và máy cắt loại trừ điểm ngắn mạch mà không làm hệ thống mất đồng bộ.

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

Nguồn dữ liệu nghiên cứu sử dụng hệ thống chuẩn 9 nút và 3 máy phát điện đồng bộ (mô hình WSCC/IEEE), bao gồm 3 máy phát, 3 máy biến áp nâng áp, 6 đường dây truyền tải và 3 nút phụ tải tập trung. Cỡ mẫu nghiên cứu bao quát 9 kịch bản tăng trưởng phụ tải liên tục từ 10% đến 90% tại các nút số 5, số 6 và số 8. Phương pháp chọn mẫu áp dụng kỹ thuật mô phỏng kịch bản có chủ đích đối với dạng sự cố ngắn mạch ba pha đối xứng tại nút số 7 kèm theo thao tác cô lập đường dây truyền tải 5-7.

Phương pháp phân tích trào lưu công suất sử dụng thuật toán lặp Newton - Raphson kết hợp thuật toán Gauss - Seidel với ngưỡng sai số hội tụ điện áp đặt ở mức 0,0001 pu. Lý do lựa chọn thuật toán Newton - Raphson là nhờ tốc độ hội tụ bậc hai vượt trội và khả năng xử lý ma trận Jacobi mạnh mẽ đối với lưới điện nhiều nút. Quá trình quá độ phi tuyến được giải tích số từng bước trên môi trường lập trình m-file của phần mềm Matlab trong giai đoạn nghiên cứu hoàn tất vào tháng 01/2013, giúp trích xuất trực quan các quỹ đạo góc pha và vận tốc góc rôto.

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à phân tích dữ liệu trên hệ thống điện 9 nút - 3 máy phát mang lại các kết quả thực nghiệm nổi bật:

  • Thuật toán trào lưu công suất Newton - Raphson đạt trạng thái hội tụ tuyệt đối sau 4 đến 5 bước lặp, xác lập chính xác điện áp tại đầu cực máy phát 1 ở mức 1,040 pu, máy phát 2 ở mức 1,025 pu và máy phát 3 ở mức 1,025 pu trước thời điểm xuất hiện ngắn mạch.
  • Xác định thời gian cắt sự cố tới hạn tại nút 7 là 0,18 giây. Khi thời gian cắt thực tế duy trì ở mức 0,15 giây, góc lệch pha giữa máy phát 2 và máy phát 3 so với máy phát chuẩn 1 nhanh chóng dao động tắt dần và thiết lập trạng thái cân bằng mới.
  • Khi thời gian cắt sự cố kéo dài vượt ngưỡng 0,22 giây (tăng 22,2% so với ngưỡng tới hạn), góc lệch pha giữa các máy phát vượt quá 180 độ, vận tốc góc gia tăng mất kiểm soát dẫn đến hiện tượng trượt bước và mất ổn định động toàn hệ thống.
  • Khảo sát độ nhạy phụ tải cho thấy khi công suất tải tại các nút 5, 6, 8 tăng 30%, thời gian cắt tới hạn giảm từ 0,18 giây xuống còn 0,14 giây (giảm 22,2%). Khi phụ tải tăng mức tối đa 90%, thời gian cắt tới hạn bị thu hẹp xuống dưới 0,09 giây, làm suy giảm biên độ dự trữ ổn định hơn 50% so với chế độ vận hành cơ sở.

Thảo luận kết quả

Cơ chế mất ổn định động phát sinh do trong thời gian tồn tại ngắn mạch ba pha tại nút 7, công suất điện từ phát ra từ máy phát 2 và máy phát 3 giảm sâu về gần bằng 0, trong khi công suất cơ của tuabin giữ nguyên giá trị định mức. Sự chênh lệch công suất cơ và công suất điện từ sinh ra mômen gia tốc cực lớn, tích lũy động năng dư thừa vào roto máy phát.

Dữ liệu mô phỏng được biểu diễn trực quan thông qua đồ thị quỹ đạo trạng thái vận tốc góc tương đối giữa các máy phát và biểu đồ biến thiên góc lệch pha theo thời gian thực từ 0 đến 5 giây. Kết quả này được so sánh và đối chiếu trực tiếp trên bảng ma trận tổng dẫn nút rút gọn trước sự cố, trong sự cố và sau khi tách đường dây 5-7. So với các nghiên cứu áp dụng tiêu chuẩn năng lượng đơn giản hay quy đổi về hệ một máy nối thanh cái vô cùng lớn, việc giải trực tiếp hệ phương trình vi phân đa máy bằng Matlab phản ánh chính xác 100% sự tương tác qua lại giữa các cụm máy phát, tăng độ tin cậy đánh giá lên hơn 25%.

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

Dựa trên các kết quả tính toán ổn định động, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật có tính ứng dụng cao:

  • Rút ngắn thời gian tác động tổng định thời của hệ thống bảo vệ rơle so lệch kỹ thuật số và máy cắt khí quyển xuống dưới 0,10 giây (100 mili giây) đối với toàn bộ các tuyến đường dây 220kV - 500kV trước năm 2025 do Tập đoàn Điện lực Việt Nam và Tổng công ty Truyền tải điện Quốc gia chủ trì thực hiện.
  • Nâng cấp và lắp đặt thiết bị tự động điều chỉnh kích từ tác động mạnh kết hợp bộ ổn định hệ thống điện (AVR/PSS) cho 100% các tổ máy phát điện có công suất trên 100 MW trong lộ trình 2024 - 2026, nhằm dập tắt nhanh dao động góc pha công suất trong vòng 1,5 đến 2 giây sau kích động.
  • Ứng dụng công cụ phần mềm mô phỏng Matlab m-file và tích hợp các thuật toán trí tuệ nhân tạo (mạng nơron, logic mờ) vào công tác tính toán chế độ vận hành 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 trước năm 2027.
  • Thiết lập quy trình tự động sa thải phụ tải theo tần số và điện áp với các ngưỡng ngắt tải linh hoạt từ 10% đến 30% tổng công suất khi xảy ra sự cố sụt giảm điện áp trên 15% tại các nút phụ tải tập trung.

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

Nội dung luận văn mang lại giá trị chuyên môn sâu rộng cho các nhóm đối tượng cụ thể:

  • Kỹ sư điều độ và vận hành hệ thống điện tại các Trung tâm Điều độ Hệ thống điện Quốc gia (A0, A1, A2, A3), áp dụng các kịch bản thời gian cắt tới hạn để xây dựng phương thức vận hành an toàn cho lưới điện truyền tải.
  • Kỹ sư thiết kế và tư vấn quy hoạch tại các Công ty Tư vấn Xây dựng Điện (PECC1, PECC2, PECC3, PECC4), sử dụng thuật toán phân tích trào lưu và mô hình 9 nút để tính toán dự phòng ổn định cho các dự án tích hợp nguồn điện mới.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Điện tại các trường đại học kỹ thuật, làm tài liệu tham khảo giảng dạy về lý thuyết quá trình quá độ điện cơ và lập trình mô phỏng Matlab.
  • Chuyên gia phát triển phần mềm mô phỏng năng lượng, kế thừa cấu trúc thuật toán m-file để xây dựng các module tính toán ổn định động tự động hóa.

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

Sự khác biệt bản chất giữa ổn định tĩnh và ổn định động trong hệ thống điện là gì?
Ổn định tĩnh là khả năng hệ thống tự phục hồi về trạng thái cân bằng ban đầu sau những kích động nhỏ và ngẫu nhiên như biến thiên phụ tải dưới 5%. Ngược lại, ổn định động liên quan đến phản ứng quá độ sau những kích động lớn như ngắn mạch ba pha, đòi hỏi hệ thống duy trì được tính đồng bộ trong khoảng thời gian từ 3 đến 5 giây đầu tiên.

Thời gian cắt tới hạn (CCT) có vai trò như thế nào trong thiết kế bảo vệ rơle?
Thời gian cắt tới hạn xác định giới hạn thời gian tối đa để hệ thống bảo vệ tác động ngắt điểm sự cố mà không làm rôto máy phát mất đồng bộ. Trong mô hình 9 nút, thời gian cắt tới hạn là 0,18 giây; nếu rơle tác động chậm hơn mức này, toàn bộ hệ thống sẽ rơi vào tình trạng sụp đổ ổn định điện cơ.

Vì sao phương pháp Newton - Raphson vượt trội hơn Gauss - Seidel trong bài toán phân bố công suất?
Phương pháp Newton - Raphson có tốc độ hội tụ bậc hai, chỉ cần từ 4 đến 5 bước lặp để đạt sai số 0,0001 pu bất kể quy mô mạng điện. Trong khi đó, thuật toán Gauss - Seidel có tốc độ hội tụ bậc một tuyến tính, đòi hỏi hàng chục đến hàng trăm bước lặp đối với các lưới điện lớn có nhiều nút máy phát.

Mức tăng phụ tải ảnh hưởng cụ thể như thế nào đến khả năng ổn định động?
Gia tăng phụ tải từ 10% đến 90% tại các nút tải làm suy giảm điện áp vận hành và đẩy góc công suất ban đầu của máy phát lên gần giới hạn 90 độ. Điều này làm thu hẹp biên độ diện tích hãm, khiến thời gian cắt tới hạn giảm mạnh từ 0,18 giây xuống dưới 0,09 giây khi tải tăng 90%.

Mô hình mô phỏng bằng Matlab m-file trong luận văn có thể mở rộng cho hệ thống thực tế không?
Mô hình hoàn toàn có khả năng mở rộng cho các hệ thống điện quy mô hàng trăm nút và hàng chục máy phát điện. Cấu trúc chương trình m-file được module hóa linh hoạt, cho phép nhập ma trận tổng dẫn kích thước lớn và tích hợp thêm các bộ điều khiển PSS, FACTS để phân tích chế độ vận hành thực tế của lưới điện Việt Nam.

Kết luận

  • Công trình đã xây dựng thành công bộ công cụ mô phỏng phân tích trào lưu công suất và ổn định động đa máy trên nền tảng Matlab m-file với độ chính xác cao.
  • Khảo sát toàn diện hệ thống chuẩn 9 nút - 3 máy phát, xác định chính xác ngưỡng thời gian cắt tới hạn 0,18 giây cho sự cố ngắn mạch trên đường dây 5-7.
  • Đánh giá định lượng tác động của việc gia tăng phụ tải từ 10% đến 90% đối với độ dự trữ ổn định và thời gian cắt tới hạn của các tổ máy phát đồng bộ.
  • Đề xuất lộ trình kỹ thuật giai đoạn 2024 - 2026 về trang bị hệ thống AVR/PSS và hoàn thiện các giải thuật bảo vệ rơle tốc độ cao dưới 0,10 giây.
  • Kế hoạch nghiên cứu giai đoạn 2027 - 2030 định hướng tích hợp mạng nơron nhân tạo và logic mờ vào bài toán nhận dạng trạng thái ổn định động trên lưới điện truyền tải siêu cao áp thực tế. Các kỹ sư, nhà nghiên cứu và học viên quan tâm có thể khai thác mô hình toán và mã nguồn mô phỏng của luận văn để phục vụ công tác đào tạo và ứng dụng thực tiễn trong ngành điện.