Tổng quan nghiên cứu

Trong hệ thống điện hiện đại, đảm bảo sự cân bằng giữa nguồn phát và phụ tải tiêu thụ là điều kiện tiên quyết để duy trì tính liên tục và chất lượng điện năng. Theo quy trình vận hành và điều độ hệ thống điện, tần số lưới chuẩn luôn được kiểm soát nghiêm ngặt quanh mốc 50 Hz, với cấp điều tần I dao động trong khoảng 50 Hz ± 0,2 Hz và cấp điều tần II trong khoảng 50 Hz ± 0,5 Hz. Khi xảy ra các sự cố nghiêm trọng như mất tổ máy phát công suất lớn hay sự cố đường dây truyền tải, tần số có thể sụt giảm nhanh chóng vượt ngưỡng an toàn. Để ngăn chặn thảm họa rã lưới diện rộng, hành động sa thải phụ tải khẩn cấp thông qua rơ le tần số thấp với thời gian tác động cực nhanh khoảng 0,04 giây (tương đương 2 chu kỳ) là giải pháp bắt buộc.

Tuy nhiên, các phương pháp sa thải phụ tải truyền thống (UFLS) hiện hành tại các trạm biến áp 110kV và 220kV chỉ tập trung vào yếu tố "lượng" công suất cần cắt mà chưa đánh giá đúng mức yếu tố "chất" – tức vị trí không gian và vai trò cấu trúc của từng nút phụ tải trên lưới. Nghiên cứu của luận văn tập trung giải quyết bài toán tối ưu hóa vị trí sa thải phụ tải trong hệ thống điện thông qua phương pháp phân tích ma trận Jacobian thu gọn và khảo sát động học góc quay rotor cực máy phát. Nghiên cứu được thử nghiệm trên mô hình lưới điện 37 nút với 9 máy phát điện, tổng công suất truyền tải cơ bản đạt 804,1 MW. Mục tiêu trọng tâm là xác định chính xác các nút nhạy cảm để sa thải lượng công suất tối thiểu nhưng mang lại tốc độ hồi phục ổn định nhanh nhất, giúp rút ngắn thời gian dao động từ 8,7% đến 15,9% và giảm thiểu tối đa tổn thất kinh tế cho người sử dụng điện.

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 ba trụ cột lý thuyết cốt lõi trong kỹ thuật hệ thống điện:

Thứ nhất, lý thuyết phân bố công suất phi tuyến và ma trận Jacobian thu gọn. Từ phương trình truyền tải công suất phức tại mỗi nút phụ tải $P_i - jQ_i = V_i^* I_i$, hệ phương trình vi phân đại số phi tuyến được tuyến tính hóa qua phép khai triển chuỗi Taylor, hình thành ma trận Jacobian gồm bốn ma trận con $J_1, J_2, J_3, J_4$. Bằng kỹ thuật phân lập công suất phản kháng xem như $\Delta Q = 0$, ma trận Jacobian thu gọn được xác lập theo công thức $J_R = J_1 - J_2 \cdot J_4^{-1} \cdot J_3$, thể hiện mối tương quan trực tiếp giữa độ lệch góc điện áp $\Delta \delta$ và độ biến thiên công suất tác dụng $\Delta P$.

Thứ hai, lý thuyết trị riêng (Eigenvalue) và không gian riêng (Eigenspace). Dựa trên phép phân tích ma trận nghịch đảo $J_R^{-1}$ qua hệ thống vector riêng bên trái $\Gamma_i$ và vector riêng bên phải $\Phi_i$, quan hệ vi phân góc tại nút $i$ được mô tả qua biểu thức $\Delta \delta_i = \sum (\Phi_{ki} \Gamma_{ik} / \lambda_i) \Delta P_k$. Trị riêng $\lambda_i$ phản ánh độ cứng vững tĩnh của nút: khi $\lambda_i > 0$, hệ thống đảm bảo ổn định; giá trị $\lambda_i$ càng lớn thì cùng một lượng biến thiên công suất $\Delta P$ chỉ gây ra độ lệch góc $\Delta \delta$ rất nhỏ.

Thứ ba, phương trình vi phân chuyển động rotor máy phát (Swing Equation): $(2H / \omega_0) \cdot (d\omega / dt) = P_m - P_e$. Luận văn sử dụng độ biến thiên góc rotor và tốc độ góc làm tiêu chuẩn then chốt để đánh giá mức độ ổn định động, khắc phục nhược điểm của việc chỉ giám sát tần số toàn cục.

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

Nghiên cứu áp dụng phương pháp chọn mẫu toàn diện (comprehensive census sampling) trên toàn bộ hệ thống điện thử nghiệm chuẩn gồm 37 nút và 9 máy phát điện. Cỡ mẫu phân tích bao gồm đầy đủ 36 giá trị trị riêng tương ứng với 36 nút độc lập (sau khi đã loại trừ nút cân bằng hệ thống) cùng toàn bộ các thông số điện áp nút dao động từ 0,9970 pu đến 1,0300 pu và góc lệch pha ban đầu từ -7,28 độ đến 2,83 độ.

Phương pháp phân tích được lựa chọn kết hợp hai công cụ chuyên dụng: phần mềm MATLAB dùng để giải thuật toán đại số ma trận cấp 36x36, tính toán vector riêng và phổ trị riêng với độ chính xác cao; phần mềm mô phỏng hệ thống điện PowerWorld Simulator dùng để mô phỏng động học miền thời gian (time-domain simulation) với bước tích phân siêu nhỏ từ 0,001 giây đến 0,01 giây. Lý do lựa chọn kết hợp hai phương pháp này là nhằm đảm bảo tính chính xác tuyệt đối trong khâu xếp hạng độ nhạy nút tĩnh, đồng thời kiểm chứng trực quan đáp ứng quá độ động lực học của máy phát điện trong các kịch bản sự cố sụt giảm công suất nguồn từ 11 MW đến 55 MW. Toàn bộ quá trình tính toán và mô phỏng được hoàn thiện vào tháng 10/2014.

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

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

Quá trình tính toán giải tích và mô phỏng số trên mô hình 37 nút đã mang lại bốn phát hiện khoa học mang tính đột phá:

Thứ nhất, phân tích phổ trị riêng đã định lượng chính xác mức độ mạnh yếu của toàn bộ 36 nút phụ tải. Giá trị trị riêng phân bố trải rộng từ 0,0011 (nút yếu nhất, xếp hạng 36) đến 1,3176 (nút mạnh nhất, xếp hạng 1). Điển hình, nút TIM69 đạt trị riêng 0,1859 (hạng 5), WEBER69 đạt 0,1541 (hạng 7), HANNAH69 đạt 0,0955 (hạng 11), WOLEN69 đạt 0,0720 (hạng 14), HALE69 đạt 0,0607 (hạng 15), và HISKY69 đạt 0,0573 (hạng 16).

Thứ hai, trong kịch bản sự cố sụt giảm nguồn phát tại nút BLT138 từ 140 MW xuống 100 MW (thiếu hụt 40 MW, công suất cần sa thải là 39 MW), cắt tải tại các nút có trị riêng khác nhau tạo ra thời gian hồi phục khác biệt rõ rệt. Khi cắt tải tại HISKY69 (trị riêng 0,0573), hệ thống mất 46 giây để ổn định; cắt tải tại WOLEN69 (trị riêng 0,0720) mất 44 giây; trong khi cắt tải tại TIM69 (trị riêng 0,1859) hệ thống phục hồi hoàn toàn chỉ sau 42 giây, nhanh hơn 4 giây (tương đương cải thiện 8,7% thời gian đáp ứng).

Thứ ba, trong kịch bản sự cố nặng tại nút BLT138 sụt giảm từ 140 MW xuống 85 MW (thiếu hụt 55 MW, sa thải 60 MW), sa thải tải tại nút WEBER69 (trị riêng 0,1541) giúp lưới ổn định ở thời điểm 37 giây, vượt trội hoàn toàn so với nút HANNAH69 mất 41 giây và nút HALE69 mất 44 giây, rút ngắn thời gian hồi phục tới 7 giây (cải thiện 15,91%).

Thứ tư, khi xảy ra sự cố dao động nhỏ tại nút ROGER69 từ 38 MW giảm xuống 27 MW (thiếu hụt 11 MW, sa thải 14,5 MW), thời gian hồi phục tại cả 6 vị trí nút khảo sát đều tiệm cận mức 29 giây đến 30 giây. Điều này chứng minh rằng ưu thế tối ưu hóa vị trí theo trị riêng phát huy tác dụng mạnh mẽ nhất trong các sự cố quy mô vừa và lớn.

Thảo luận kết quả

Các phát hiện thực nghiệm được minh chứng rõ nét khi dữ liệu được tổng hợp qua bảng ma trận độ nhạy và hệ thống biểu đồ đường cong chuyển động rotor máy phát theo thời gian. Trên đồ thị động học quá độ từ 0 đến 60 giây, biên độ dao động góc quay rotor tại các nút có trị riêng cao suy giảm theo hàm mũ với độ dốc lớn hơn hẳn, triệt tiêu hiện tượng vọt lố góc lệch và duy trì góc rotor ổn định quanh ngưỡng dưới 23,10 độ.

Nguyên nhân vật lý cốt lõi là do nút có trị riêng $\lambda_i$ lớn sở hữu độ liên kết điện thế mạnh mẽ với các máy phát lân cận, khiến phản hồi cơ học $P_m - P_e$ được cân bằng tức thời mà không gây sốc điện áp cục bộ. So với các giải pháp sa thải tải truyền thống chỉ dựa vào rơ le tần số đơn thuần hay các nghiên cứu ứng dụng mạng nơ-ron nhân tạo (ANN) đòi hỏi tập dữ liệu quá khứ khổng lồ, giải pháp phân tích ma trận Jacobian thu gọn mang tính giải tích tất định, tính toán thời gian thực với độ tin cậy tuyệt đối, giúp chuyển đổi căn bản chiến lược sa thải phụ tải từ thụ động sang tối ưu hóa chủ động cả về lượng và chất.

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

Dựa trên kết quả nghiên cứu, luận văn đề xuất bốn nhóm giải pháp hành động cụ thể nhằm hiện đại hóa công tác vận hành hệ thống điện:

Thứ nhất, tích hợp module thuật toán phân tích ma trận Jacobian thu gọn vào hệ thống phần mềm SCADA/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 (A1, A2, A3). Hệ thống cần tự động tính toán và cập nhật bảng xếp hạng trị riêng các nút phụ tải theo chu kỳ thời gian thực từ 5 đến 15 phút, thiết lập sẵn danh mục nút ưu tiên cắt tải tối ưu.

Thứ hai, nâng cấp và lập trình lại thông số chỉnh định cho 100% rơ le kỹ thuật số bảo vệ tần số thấp (tiêu biểu như rơ le SEL-351) tại các trạm biến áp 110kV và 220kV. Cài đặt các cấp thời gian trễ linh hoạt từ 0,04 giây (2 chu kỳ) kết hợp phân tầng thứ tự cắt tải theo trọng số trị riêng thay vì chỉ dựa vào phân loại tải chính trị hoặc kinh tế đơn thuần, hoàn thành triển khai trong giai đoạn 2025 - 2027.

Thứ ba, Tập đoàn Điện lực Việt Nam (EVN) cần ban hành bộ quy chuẩn kỹ thuật mới về sa thải phụ tải thông minh (ILS). Quy chuẩn cần kết hợp đồng thời hai điều kiện biên: xác định lượng phụ tải cần cắt theo tốc độ suy giảm tần số $df/dt$ và lựa chọn vị trí nút cắt theo hệ số nhạy trị riêng $\lambda_i$, hướng tới mục tiêu giảm ít nhất 10% tổng công suất phụ tải bị cắt không cần thiết trong mỗi đợt sự cố.

Thứ tư, đầu tư mở rộng hệ thống đo lường góc pha diện rộng WAMS/PMU tại tất cả các thanh cái trạm 220kV và 500kV với tốc độ lấy mẫu từ 20 đến 50 mẫu mỗi giây, cung cấp dữ liệu góc quay rotor trực tiếp hỗ trợ thuật toán điều khiển thích nghi trong giai đoạn 2026 - 2030.

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

Luận văn là tài liệu tham khảo chuyên sâu, mang lại giá trị thực tiễn cao cho bốn nhóm đối tượng:

Nhóm 1: Các kỹ sư điều độ viên và chuyên gia vận hành hệ thống điện tại các Trung tâm Điều độ A0, A1, A2, A3. Nghiên cứu cung cấp phương pháp luận định lượng giúp xây dựng kịch bản điều hành khẩn cấp, rút ngắn từ 8% đến 15% thời gian khôi phục hệ thống khi xảy ra sự cố lớn.

Nhóm 2: Kỹ sư thiết kế phương thức và cấu hình rơ le bảo vệ tại các Công ty Truyền tải điện và Công ty Điện lực. Luận văn cung cấp hướng dẫn chi tiết về cách phối hợp logic sa thải phụ tải trên các thiết bị rơ le kỹ thuật số đa chức năng như SEL-351 tại các trạm 110kV và 220kV.

Nhóm 3: 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, Hệ thống Điện tại các trường đại học kỹ thuật. Đây là tài liệu mẫu mực về ứng dụng đại số tuyến tính nâng cao, giải tích ma trận Jacobian trên MATLAB và mô phỏng ổn định động trên PowerWorld Simulator cho lưới 37 nút 9 máy phát.

Nhóm 4: Các nhà hoạch định chính sách thị trường điện và quản lý an ninh năng lượng. Tài nguyên này hỗ trợ thiết lập cơ chế điều tiết phụ tải kỹ thuật minh bạch, đảm bảo tính công bằng và hiệu quả kinh tế trong môi trường thị trường điện cạnh tranh.

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

Tại sao không thể chỉ dựa vào tần số lưới điện để ra quyết định sa thải phụ tải?

Tần số là đại lượng biến thiên toàn cục phản ánh sự mất cân bằng công suất nhưng không thể hiện được vị trí xung yếu. Kết quả mô phỏng cho thấy khi cắt cùng 39 MW tải, cắt tại nút TIM69 chỉ mất 42 giây để hồi phục, trong khi cắt tại HISKY69 mất 46 giây. Do đó, chỉ dựa vào tần số sẽ bỏ qua vai trò cấu trúc của vị trí nút.

Ma trận Jacobian thu gọn đóng vai trò gì trong việc xác định nút phụ tải mạnh hay yếu?

Ma trận Jacobian thu gọn $J_R$ thiết lập quan hệ vi phân tuyến tính giữa góc điện áp và công suất tác dụng. Phổ trị riêng $\lambda_i$ của $J_R$ phân loại 36 nút với giá trị từ 0,0011 đến 1,3176. Nút có trị riêng càng lớn thì biến thiên công suất $\Delta P$ càng ít gây biến động góc $\Delta \delta$, chứng minh nút đó chịu tải vững chắc nhất.

Góc quay rotor máy phát có ý nghĩa như thế nào trong đánh giá ổn định hệ thống?

Góc rotor phản ánh trực tiếp trạng thái cân bằng cơ - điện qua phương trình chuyển động $(2H/\omega_0) \cdot (d\omega/dt) = P_m - P_e$. Theo dõi góc rotor giúp đánh giá chính xác dao động quá độ; ví dụ tại nút LAUF69, chọn đúng nút cắt tải giúp góc dao động rotor được dập tắt ổn định dưới ngưỡng 23,10 độ.

Phương pháp tối ưu sa thải phụ tải đề xuất có khả thi với lưới điện Việt Nam không?

Phương pháp hoàn toàn tương thích với quy trình điều độ 3 cấp của Việt Nam ở các ngưỡng 50 Hz ± 0,2 Hz và 50 Hz ± 0,5 Hz. Thuật toán có thể nhúng trực tiếp vào các rơ le hiện hữu tại trạm 220kV với thời gian trễ chỉ 0,04 giây kết hợp SCADA trung tâm mà không cần thay thế hạ tầng đắt đỏ.

Khi nào phương pháp sa thải phụ tải theo trị riêng phát huy hiệu quả cao nhất?

Phương pháp phát huy tối đa ưu thế trong các sự cố quy mô vừa và lớn với lượng công suất sụt giảm từ 40 MW đến 55 MW, giúp rút ngắn thời gian ổn định từ 4 đến 7 giây. Với các sự cố nhỏ dưới 11 MW, thời gian hồi phục giữa các nút chỉ chênh lệch dưới 1 giây.

Kết luận

Nghiên cứu của luận văn thạc sĩ đã mang lại những đóng góp khoa học và thực tiễn mang tính then chốt:

  • Xác lập cơ sở khoa học định lượng yếu tố "chất" của phụ tải thông qua phân tích phổ trị riêng của ma trận Jacobian thu gọn $J_R$, giải quyết triệt để hạn chế của phương pháp cắt tải truyền thống.
  • Chứng minh thực nghiệm trên mô hình 37 nút 9 máy phát rằng việc tối ưu hóa vị trí sa thải phụ tải giúp rút ngắn thời gian phục hồi lưới điện lên tới 15,91% (từ 44 giây xuống 37 giây khi sa thải 60 MW).
  • Đưa chuyển động góc quay rotor cực máy phát làm tiêu chuẩn kép song hành cùng tần số để đánh giá chính xác độ ổn định động học sau sự cố.
  • Hoàn thiện bộ công cụ kết hợp giải tích ma trận cấp 36x36 trên MATLAB và mô phỏng động học quá độ miền thời gian trên PowerWorld Simulator.
  • Đề xuất lộ trình công nghệ chuyển đổi hệ thống sa thải phụ tải thông minh kết hợp thiết bị đo lường WAMS/PMU cho hệ thống điện quốc gia giai đoạn 2025 - 2030.

Các đơn vị vận hành điều độ lưới điện, viện nghiên cứu và doanh nghiệp truyền tải điện cần nhanh chóng tiếp cận, ứng dụng khung phương pháp luận của đề tài để nâng cao độ tin cậy cung cấp điện và bảo vệ an toàn hệ thống năng lượng quốc gia.