Tổng quan nghiên cứu

Trong hệ thống điện hiện đại, sự gia tăng nhanh chóng của các trung tâm phụ tải và nguồn điện phân tán khiến mạng lưới vận hành ngày càng tiệm cận giới hạn kỹ thuật cho phép. Thống kê vận hành lưới điện chỉ ra rằng sự cố ngắn mạch một pha chiếm từ 70% đến 90% tổng số sự cố, trong khi ngắn mạch ba pha tuy chỉ chiếm khoảng 5% đến 10% nhưng lại tiềm ẩn nguy cơ nghiêm trọng nhất dẫn đến mất đồng bộ. Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết bài toán đánh giá ổn định quá độ và ổn định tín hiệu nhỏ trong hệ thống điện nhiều máy phát liên kết, khắc phục triệt để nhược điểm của các phương pháp quy đổi tương đương truyền thống vốn dễ làm sai lệch các chế độ dao động cục bộ.

Mục tiêu cụ thể của công trình tập trung vào việc xây dựng hệ phương trình vi phân phi tuyến hoàn chỉnh mô tả máy phát đồng bộ có xét đến ảnh hưởng của cuộn kích từ, cuộn cản và bộ ổn định hệ thống điện PSS; từ đó đánh giá định lượng tính ổn định tĩnh và ổn định động trên hệ thống một máy nối thanh cái vô cùng lớn cũng như hệ thống kiểm chứng chuẩn hai vùng bốn máy của Kundur với tổng công suất định mức 2220 MVA. Phạm vi nghiên cứu được thực hiện trong không gian mô phỏng chuyên sâu, bám sát các điều kiện vận hành của lưới điện Việt Nam với trục truyền tải siêu cao áp 500 kV Bắc - Nam. Ý nghĩa khoa học và thực tiễn của đề tài thể hiện qua việc xác định chính xác thời gian loại trừ sự cố tới hạn trong khoảng 0,15s đến 0,20s và cải thiện hệ số suy giảm dao động lên trên 10%, tạo nền tảng vững chắc giúp các kỹ sư điều độ hoạch định phương thức vận hành an toàn, ngăn ngừa thảm họa rã lưới diện rộng gây thiệt hại hàng triệu USD cho nền kinh tế.

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 ứng dụng lý thuyết ổn định Lyapunov và phương pháp phân tích dao động bé tuyến tính hóa quanh điểm làm việc xác lập. Mô hình máy phát đồng bộ được xây dựng dựa trên phép biến đổi Park đưa hệ phương trình stator và rotor về hệ tọa độ quay đồng bộ dq0. Mô hình chi tiết bậc cao khảo sát đầy đủ phản ứng phần ứng, từ thông cuộn kích từ dọc trục d ($L_{fd}$), một cuộn cản dọc trục d ($L_{1d}$) và hai cuộn cản ngang trục q ($L_{1q}, L_{2q}$).

Khung lý thuyết tích hợp ba khái niệm nền tảng:

  • Ổn định tĩnh (Small-signal stability): Khả năng hệ thống duy trì trạng thái đồng bộ khi chịu các kích động nhỏ xuất hiện liên tục từ sự biến thiên phụ tải.
  • Ổn định động (Transient stability): Khả năng hệ thống tự phục hồi về trạng thái cân bằng xác lập sau những nhiễu loạn lớn như ngắn mạch 3 pha hoặc sa thải đột ngột đường dây truyền tải.
  • Tiêu chuẩn cân bằng diện tích (Equal Area Criterion - EAC): Nguyên lý cân bằng giữa năng lượng gia tốc $A_1$ tích lũy trong quá trình ngắn mạch và năng lượng giảm tốc $A_2$ giải phóng sau khi cắt sự cố dựa trên phương trình góc quay rotor: $2H \frac{d^2\delta}{dt^2} = P_m - P_e$.

Hệ thống kích từ được mô hình hóa cùng bộ điều khiển PSS với hệ số khuếch đại $K_A = 200$ và hằng số quán tính $H = 3,5$, thiết lập nên vector trạng thái 10 biến mô tả đầy đủ đáp ứng cơ - điện.

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

Nghiên cứu sử dụng nguồn dữ liệu kỹ thuật chuẩn từ mô hình hệ thống điện hai vùng bốn máy của Kundur, bao gồm 4 tổ máy phát đồng bộ công suất 555 MVA mỗi máy, cấp điện áp danh định 24 kV và tần số làm việc 60 Hz (hoặc 50 Hz quy đổi). Quy mô mẫu thử nghiệm này được lựa chọn có chủ đích vì cấu trúc hai vùng liên kết qua đường dây truyền tải dài tái hiện hoàn hảo hiện tượng dao động liên vùng tần số thấp (từ 0,1 Hz đến 2,0 Hz), phản ánh chân thực đặc tính của lưới điện thực tế.

Lý do lựa chọn phương pháp phân tích miền thời gian và phân tích ma trận không gian trạng thái kích thước 10x10 thay vì phương pháp tương đương hóa một máy là nhằm bảo toàn toàn bộ các biến trạng thái nội tại của từng máy phát riêng lẻ. Các phương pháp tích phân số từng bước được áp dụng gồm thuật toán Euler, Runge-Kutta bậc 2, Runge-Kutta bậc 4 và phương pháp phân đoạn liên tiếp với bước thời gian tích phân $\Delta t = 0,02s$ đến $0,1s$. Toàn bộ thuật toán được lập trình bằng ngôn ngữ M-file trong môi trường phần mềm MATLAB, hoàn thành chuỗi đánh giá toàn diện trong giai đoạn nghiên cứu năm 2015.

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

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

Quá trình phân tích không gian trạng thái và mô phỏng số đã mang lại các phát hiện cốt lõi sau:

Thứ nhất, phân tích trị riêng của ma trận trạng thái hệ thống $E = [A] + [B][C]^{-1}[D]$ cho thấy khi hệ số mômen cản $K_D \ge 0$, toàn bộ 10 trị riêng $\lambda_i$ đều có phần thực âm, khẳng định hệ thống đạt trạng thái ổn định tín hiệu nhỏ tuyệt đối. Ngược lại, khi $K_D = -10$, xuất hiện các trị riêng có phần thực dương, dẫn đến mất ổn định không chu kỳ hoặc mất ổn định dao động.

Thứ hai, các mode dao động được phân loại rõ rệt thành 3 nhóm: mode điện, mode cơ và mode kích từ. Trong khi các mode liên quan đến biến $V_2, V_3, V_s$ là mode tắt dần không dao động, thì mode góc pha rotor $\delta$ và điện áp kích từ $V_1$ biểu hiện tần số dao động tự nhiên cao với hệ số suy giảm dao động ban đầu rất thấp (dưới 0,03).

Thứ ba, việc tích hợp bộ ổn định hệ thống điện PSS với độ lợi $K_A = 200$ đã tăng cường hệ số giảm chấn của các mode cơ điện nguy hiểm từ mức dưới 0,05 lên mức tối ưu trên 0,15, giúp rút ngắn thời gian dập tắt dao động quá độ từ hơn 4,5 giây xuống còn dưới 1,8 giây (giảm khoảng 60% thời gian quá độ).

Thứ tư, khảo sát ổn định động với sự cố ngắn mạch ba pha cho thấy góc cắt tới hạn $\delta_c$ và thời gian cắt giới hạn $t_{cr}$ có mối tương quan phi tuyến chặt chẽ. Khi thời gian tồn tại ngắn mạch vượt quá 0,18 giây, diện tích gia tốc $A_1$ vượt qua diện tích giảm tốc cực đại $A_2$, khiến góc rotor phân kỳ nhanh chóng và hệ thống rơi vào trạng thái mất đồng bộ.

Thảo luận kết quả

Cơ chế vật lý dẫn đến sự dao động kéo dài khi thiếu PSS là do tác động bất lợi của bộ tự động điều chỉnh điện áp (AVR) độ lợi cao, vô tình tạo ra thành phần mômen cản âm trên trục máy phát. Dữ liệu đánh giá ổn định có thể được trực quan hóa tối ưu thông qua biểu đồ quỹ đạo góc pha rotor $\delta(t)$ theo thời gian thực (từ 0 đến 5 giây) kết hợp cùng bảng thống kê so sánh 10 trị riêng ứng với các giá trị $K_D = -10, 0, 10$. Biểu đồ quỹ tích nghiệm sẽ minh họa rõ sự dịch chuyển của các cặp nghiệm phức từ nửa mặt phẳng bên phải sang nửa mặt phẳng bên trái khi có sự can thiệp của bộ PSS.

So sánh với các nghiên cứu tương đương hóa một máy nối thanh cái vô cùng lớn, phương pháp đánh giá từng máy riêng lẻ của luận văn phản ánh độ chính xác cao hơn khoảng 15% đến 20% đối với các tương tác dao động qua lại giữa các máy phát thuộc hai vùng khác nhau. Kết quả này cung cấp cái nhìn chi tiết về sự phân bố công suất và điện áp tại từng nút máy phát mà các phương pháp giản lược truyền thống không thể phát hiện.

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

Dựa trên kết quả nghiên cứu, bốn giải pháp kỹ thuật cụ thể được đề xuất nhằm nâng cao giới hạn ổn định cho hệ thống điện:

  • Thiết lập cấu hình bảo vệ rơle tác động siêu tốc: Đơn vị quản lý vận hành cần cài đặt rơle bảo vệ so lệch đường dây và khoảng cách với thời gian tác động loại trừ sự cố ngắn mạch ba pha dưới 0,12 giây (120 ms) trên các cung đoạn 220 kV - 500 kV trọng yếu, do Trung tâm Điều độ Hệ thống điện Quốc gia (A0) và Công ty Truyền tải điện chỉ đạo hoàn thành trong lộ trình 6 đến 12 tháng.
  • Chuẩn hóa và tối ưu thông số bộ ổn định PSS: Tiến hành hiệu chỉnh tham số bộ kích từ với hệ số khuếch đại $K_A = 200$, kết hợp các khâu trễ - vượt pha $T_1 - T_4$ nhằm đảm bảo hệ số giảm chấn $\xi \ge 0,10$ cho toàn bộ các tổ máy nhiệt điện và thủy điện có công suất từ 100 MW trở lên, thực hiện bởi các nhà máy điện trong vòng 18 tháng.
  • Ứng dụng hệ thống đo lường góc pha diện rộng WAMS/PMU: Triển khai lắp đặt các thiết bị đo lường đồng bộ pha (PMU) với tốc độ lấy mẫu 50 mẫu/giây tại tất cả các trạm biến áp nút liên kết vùng, phục vụ giám sát trực tiếp trị riêng và nhận diện sớm các dao động tần số thấp từ 0,2 Hz đến 1,0 Hz, do Tập đoàn Điện lực Việt Nam (EVN) chủ trì trong giai đoạn 24 tháng.
  • Chu kỳ hóa công tác tính toán ổn định động trực tuyến: Ban hành quy chuẩn tính toán mô phỏng miền thời gian định kỳ với bước chia tích phân $\Delta t \le 0,02s$ trước mỗi đợt chuyển đổi phương thức vận hành mùa khô và mùa mưa, hướng tới mục tiêu giảm thiểu trên 95% nguy cơ quá độ dẫn đến mất liên lạc giữa các miền.

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

Luận văn là tài liệu chuyên khảo mang giá trị học thuật và ứng dụng cao cho 4 nhóm đối tượng:

  • Kỹ sư điều độ và phương thức hệ thống điện: Ứng dụng trực tiếp các phương trình trạng thái và giới hạn thời gian cắt sự cố để tính toán giới hạn truyền tải ổn định, lập phương thức vận hành an toàn cho các đường dây liên kết 500 kV liên miền.
  • Kỹ sư thiết kế và quy hoạch nguồn lưới: Sử dụng mô hình chuẩn hai vùng bốn máy (2220 MVA) để thẩm định tính khả thi kỹ thuật, đánh giá tác động ổn định khi đấu nối các cụm nhà máy điện mới quy mô từ 500 MW đến 1200 MW vào lưới điện quốc gia.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Kỹ thuật điện: Tham khảo cấu trúc toán học chi tiết của hệ thống 10 phương trình vi phân tuyến tính hóa, ma trận không gian trạng thái và phương pháp lập trình mô phỏng trên MATLAB phục vụ công tác giảng dạy và nghiên cứu chuyên sâu.
  • Chuyên gia phát triển phần mềm mô phỏng năng lượng: Khai thác các thuật toán tích phân số phân đoạn liên tiếp và Runge-Kutta bậc 4 làm nền tảng xây dựng các module thuật toán đánh giá ổn định động trong các phần mềm tính toán lưới điện chuyên dụng.

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

Phương pháp phân tích ổn định từng máy riêng lẻ có ưu điểm gì vượt trội so với phương pháp biến đổi tương đương một máy?

Phương pháp đánh giá từng máy riêng lẻ giữ nguyên toàn bộ hệ phương trình vi phân nội tại của từng tổ máy phát. Điều này giúp nhận diện chính xác các mode dao động cục bộ giữa các máy phát trong cùng một nhà máy và mode dao động liên vùng tần số thấp từ 0,1 Hz đến 2,0 Hz, loại bỏ sai số ước lượng từ 15% đến 20% so với phương pháp tương đương hóa thanh cái vô cùng lớn.

Tại sao sự cố ngắn mạch ba pha chiếm tỷ lệ chỉ 5% đến 10% nhưng luôn là đối tượng khảo sát trọng tâm trong ổn định động?

Mặc dù ngắn mạch một pha chiếm 70% đến 90% tần suất thực tế, ngắn mạch ba pha làm suy giảm điện áp nút về mức gần bằng 0, triệt tiêu hoàn toàn khả năng truyền tải công suất điện $P_e$ trên đường dây. Điều này tạo ra độ chênh lệch công suất cơ - điện cực đại, khiến rotor máy phát tăng tốc dữ dội và đe dọa trực tiếp đến tính ổn định của toàn hệ thống.

Hệ số mômen cản $K_D$ ảnh hưởng trực tiếp thế nào đến các trị riêng của ma trận không gian trạng thái?

Hệ số mômen cản $K_D$ quyết định dấu của phần thực trong các nghiệm trị riêng cơ điện. Khi $K_D > 0$, phần thực của trị riêng luôn âm, đảm bảo biên độ dao động tắt dần theo hàm mũ. Khi $K_D < 0$, hệ thống xuất hiện trị riêng có phần thực dương, biến các kích động nhỏ thành dao động khuếch đại không kiểm soát dẫn đến mất ổn định.

Bộ ổn định hệ thống điện PSS triệt tiêu dao động cơ điện bằng cơ chế nào?

Bộ PSS trích xuất tín hiệu độ lệch tốc độ góc $\Delta \omega_r$ hoặc độ lệch công suất gia tốc $\Delta P_a$, đưa qua các khâu lọc và sớm pha để tác động trực tiếp vào mạch điều khiển kích từ $E_{fd}$. Quá trình này tạo ra một thành phần mômen điện đồng pha với độ lệch tốc độ, bù đắp phần mômen cản âm do AVR gây ra và tăng hệ số suy giảm dao động lên trên 10%.

Vì sao mô hình toán học cần tích hợp cả cuộn kích từ và 3 cuộn cản trên các trục dq?

Việc mô hình hóa chi tiết cuộn kích từ cùng cuộn cản trục d ($L_{1d}$) và hai cuộn cản trục q ($L_{1q}, L_{2q}$) cho phép mô tả chính xác hiện tượng bão hòa từ thông và quá trình tắt dần của dòng điện xoáy trong thân rotor. Nhờ đó, mô hình tái hiện chuẩn xác đáp ứng điện từ siêu quá độ trong khoảng thời gian từ 0,03s đến 0,05s ngay sau sự cố.

Kết luận

  • Xây dựng thành công hệ mô hình toán học giải tích bậc cao cho máy phát đồng bộ tích hợp đầy đủ cuộn kích từ, cuộn cản và bộ điều khiển PSS với hệ thống 10 biến trạng thái.
  • Đánh giá toàn diện ổn định tín hiệu nhỏ thông qua phân tích ma trận trị riêng, làm sáng tỏ vai trò quyết định của hệ số mômen cản $K_D$ đối với sự suy giảm dao động.
  • Định lượng chính xác giới hạn ổn định động trên hệ thống chuẩn hai vùng bốn máy 2220 MVA bằng các giải thuật tích phân số Runge-Kutta và phương pháp cân bằng diện tích.
  • Xác lập ngưỡng thời gian cắt ngắn mạch giới hạn trong khoảng 0,15s đến 0,20s nhằm đảm bảo an toàn vận hành cho các đường dây liên kết.
  • Cung cấp công cụ tính toán M-file tin cậy trên nền tảng MATLAB, có khả năng mở rộng trực tiếp cho việc phân tích các chế độ vận hành phức tạp của hệ thống điện Việt Nam.

Định hướng tiếp theo tập trung vào việc tối ưu hóa tốc độ thuật toán để ứng dụng đánh giá ổn định trực quan theo thời gian thực trong 12 tháng tới, đồng thời tích hợp các nguồn năng lượng tái tạo gió và mặt trời vào mô hình ổn định nhiều máy trong vòng 24 tháng. Các đơn vị quản lý vận hành, viện nghiên cứu và trường đại học quan tâm có thể khai thác toàn văn công trình luận văn để ứng dụng vào công tác đào tạo và nâng cao độ tin cậy cung cấp điện quốc gia.