Tổng quan nghiên cứu

Sự chuyển dịch năng lượng từ các nguồn nhiên liệu hóa thạch sang năng lượng tái tạo đang diễn ra mạnh mẽ trên quy mô toàn cầu và tại Việt Nam. Theo định hướng tại Quyết định số 11/2017/QĐ-CP của Thủ tướng Chính phủ, các dự án điện mặt trời quy mô lớn liên tục được cấp phép và xây dựng. Điển hình là cụm trang trại quang điện tại tỉnh Ninh Thuận với diện tích 500 ha, công suất lắp đặt 304 MW, cung cấp sản lượng khoảng 650 MWh hàng năm hòa vào hệ thống điện quốc gia. Tuy nhiên, việc tích hợp nguồn điện mặt trời công suất lớn thay thế các tổ máy phát điện đồng bộ truyền thống đang đặt ra những thách thức kỹ thuật lớn về chế độ xác lập và chế độ quá độ.

Vấn đề cốt lõi của nghiên cứu là đánh giá định lượng mức độ ảnh hưởng của nguồn quang điện công suất lớn đến tính ổn định tĩnh và ổn định động của lưới điện truyền tải. Mục tiêu cụ thể là xây dựng mô hình toán học chi tiết, tính toán chế độ làm việc và mô phỏng kiểm chứng khả năng vận hành ổn định của hệ thống 3 máy 9 bus chuẩn IEEE khi tích hợp hệ thống pin năng lượng mặt trời công suất lớn thay thế máy phát điện truyền thống với tỷ lệ thâm nhập lên tới 51%.

Phạm vi nghiên cứu tập trung vào mô hình lưới điện truyền tải cấp điện áp 230 kV và 18 kV, thực hiện hoàn tất vào tháng 3 năm 2019. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng giúp kiểm soát thời gian giải trừ sự cố ngắn mạch dưới 0,15 giây và khống chế độ lệch điện áp nút trong ngưỡng an toàn ±5%, tạo tiền đề kỹ thuật vững chắc cho việc mở rộng nguồn năng lượng tái tạo trong tương lai.

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 trên nền tảng lý thuyết ổn định hệ thống điện và mô hình hóa thiết bị điện hiện đại:

  • Lý thuyết ổn định hệ thống điện: Khảo sát hai trạng thái cốt lõi là ổn định tĩnh và ổn định động. Ổn định tĩnh được phân tích thông qua phương pháp dao động bé dựa trên ma trận trạng thái tuyến tính hóa. Ổn định động được phân tích qua phương trình chuyển động rotor và tiêu chuẩn cân bằng diện tích nhằm xác định góc quay tới hạn khi xảy ra các kích động lớn.
  • Mô hình máy phát điện đồng bộ hai trục: Hệ thống máy phát được mô hình hóa chi tiết trên hệ trục trực giao d-q thông qua 4 phương trình vi phân phi tuyến kết hợp hệ thống kích từ loại I (Exciter Type 1) có tính đến hiện tượng bão hòa từ trường.
  • Mô hình toán học trạm pin năng lượng mặt trời: Mô hình hóa mạch tương đương của tế bào quang điện bao gồm dòng quang điện, dòng bão hòa, điện trở nối tiếp và điện trở sun. Hệ thống tích hợp thuật toán dò tìm điểm công suất cực đại (MPPT) cùng mô hình nghịch lưu biến tần DC-AC điều chế độ rộng xung (PWM) chuyển đổi sang khung tọa độ d-q.
  • Các khái niệm cốt lõi: Công suất tác dụng P, công suất phản kháng Q, góc lệch pha rotor δ, tần số định mức 50 Hz, điện áp nút U và hằng số quán tính quay H.

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

Nghiên cứu sử dụng hệ thống dữ liệu chuẩn IEEE 3 máy 9 bus với cấp công suất cơ bản 100 MVA, điện áp thanh cái 230 kV và 18 kV. Mẫu nghiên cứu tích hợp gồm 363 tổ hợp pin mặt trời công suất 500 kVA đấu nối song song tại Bus số 2. Mỗi khối 500 kVA được cấu thành từ 140 module ghép nối tiếp và 20 nhánh ghép song song, sử dụng module quang điện 175 Wp với điện áp hở mạch 43,8 V và dòng điện ngắn mạch 4,58 A được trích xuất từ thư viện DIgSILENT PowerFactory.

Phương pháp chọn mẫu dựa trên cấu trúc chuẩn IEEE 3 máy 9 bus bởi đây là mô hình kinh điển đại diện đầy đủ cho hệ thống điện truyền tải đa nguồn, có tính liên kết vùng và đa dạng phụ tải. Phương pháp phân tích bao gồm:

  • Phân tích miền tần số: Tuyến tính hóa hệ phương trình vi phân đại số (DAE) tại điểm làm việc cân bằng để thành lập ma trận trạng thái A kích thước 3x3 và tính toán các giá trị riêng (Eigenvalues).
  • Mô phỏng miền thời gian: Ứng dụng phần mềm Matlab để giải tích phân số phân đoạn các phương trình vi phân phi tuyến, khảo sát phản ứng quá độ của hệ thống trong thời gian 10 giây đối với sự cố ngắn mạch 3 pha và sự biến thiên bức xạ mặt trời.

Lý do lựa chọn phương pháp này là vì sự kết hợp giữa phân tích miền tần số và mô phỏng miền thời gian cho phép đánh giá toàn diện cả tính chất suy giảm dao động tín hiệu bé lẫn hành vi quá độ phi tuyến phức tạp của hệ thống. Timeline nghiên cứu được thực hiện đồng bộ trong thời gian 12 tháng, nghiệm thu vào tháng 3 năm 2019.

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 phần mềm chuyên dụng đã ghi nhận các phát hiện khoa học quan trọng:

  • Tính khả thi của tỷ lệ thâm nhập 51%: Việc thay thế hoàn toàn máy phát đồng bộ số 2 bằng 363 tổ hợp pin quang điện 500 kVA (tổng công suất định mức đạt xấp xỉ 163,35 MW) tại Bus 2 vẫn đảm bảo sự cân bằng công suất tác dụng và phản kháng toàn hệ thống ở chế độ xác lập 100 MVA.
  • Đảm bảo ổn định tĩnh tuyệt đối: Kết quả phân tích miền tần số cho thấy toàn bộ các giá trị riêng của ma trận trạng thái hệ thống đều có phần thực mang giá trị âm. Các cặp nghiệm phức liên hợp đặc trưng cho chế độ dao động cơ điện đều có tỷ số cản dịu vượt mức 0,1, minh chứng các dao động nhỏ bị triệt tiêu hoàn toàn sau thời gian dưới 2 giây.
  • Giới hạn ổn định động khi ngắn mạch 3 pha: Khi mô phỏng ngắn mạch 3 pha tại Bus 2 (dạng sự cố chiếm 5-10% trong thực tế nhưng nguy hiểm nhất), hệ thống vẫn giữ được tính ổn định nếu thời gian cắt sự cố được khống chế dưới 0,15 giây. Sau khi cô lập đường dây sự cố, góc lệch rotor của các máy phát còn lại nhanh chóng ổn định về điểm cân bằng mới sau 3,5 giây dao động.
  • Khả năng thích ứng với biến thiên bức xạ: Khi cường độ bức xạ mặt trời suy giảm đột ngột 40% (từ 1000 W/m² xuống 600 W/m²), điện áp tại các nút tải chỉ dao động trong biên độ 3,8% và bộ điều khiển nghịch lưu đã tái lập điện áp định mức sau 1,2 giây.

Thảo luận kết quả

Khi tích hợp nguồn điện mặt trời qua bộ biến tần tĩnh thay cho máy phát đồng bộ truyền thống, hằng số quán tính quay H của hệ thống bị suy giảm cục bộ, làm tăng tốc độ biến thiên góc rotor khi xảy ra sự cố. Tuy nhiên, hệ thống vẫn duy trì được độ ổn định động nhờ vào tốc độ phản ứng cực nhanh của bộ điều chế xung PWM trên biến tần DC-AC cùng sự hỗ trợ công suất phản kháng từ hệ thống kích từ loại I của các máy phát số 1 và số 3.

So với các công bố khoa học trước đây vốn chỉ tập trung vào tỷ lệ tích hợp phân tán nhỏ từ 10% đến 20%, nghiên cứu này đã chứng minh hệ thống điện truyền tải hoàn toàn có thể hấp thụ tỷ lệ điện mặt trời tập trung lên tới 51% mà không gây sụp đổ điện áp. Toàn bộ dữ liệu nghiên cứu được biểu diễn trực quan qua bảng ma trận trị riêng và các biểu đồ đặc tính góc công suất theo thời gian (P-δ), minh chứng rõ ràng sự cân bằng giữa diện tích tăng tốc A1 và diện tích giảm tốc A2 theo đúng tiêu chuẩn cân bằng diện tích.

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

Dựa trên kết quả nghiên cứu thực nghiệm, các giải pháp kỹ thuật và chính sách được đề xuất cụ thể như sau:

  • Nâng cấp hệ thống rơle bảo vệ đường dây: Trang bị rơle kỹ thuật số tác động nhanh và máy cắt hiện đại tại các trạm biến áp 110 kV và 230 kV đấu nối nguồn năng lượng tái tạo, đảm bảo thời gian loại trừ ngắn mạch 3 pha dưới 0,10 giây. Chủ thể thực hiện là Tập đoàn Điện lực Việt Nam trong giai đoạn 2020-2023.
  • Bắt buộc tích hợp chế độ điều khiển hỗ trợ lưới (Grid-forming): Áp dụng chuẩn điều khiển điện áp và tần số tự động cho các bộ biến tần nghịch lưu DC-AC tại các nhà máy điện mặt trời, nhằm duy trì độ lệch điện áp thanh cái không vượt quá ±3% khi bức xạ dao động trên 30%. Chủ thể thực hiện là các chủ đầu tư nhà máy điện mặt trời từ năm 2021-2024.
  • Lắp đặt hệ thống lưu trữ năng lượng pin (BESS): Bổ sung trạm tích trữ năng lượng quy mô 15-20% công suất đỉnh của nhà máy quang điện để cung cấp quán tính ảo và ổn định tần số hệ thống trong dải 49,8 - 50,2 Hz. Đơn vị chủ trì là Trung tâm Điều độ Hệ thống điện Quốc gia với lộ trình đến năm 2025.
  • Hoàn thiện khung quy chuẩn kỹ thuật đấu nối: Ban hành các tiêu chuẩn kỹ thuật bắt buộc về khả năng duy trì vận hành khi sụt giảm điện áp (Fault Ride-Through - FRT) cho các nguồn năng lượng tái tạo công suất trên 50 MW. Đơn vị chủ trì là Bộ Công Thương trong giai đoạn 2020-2022.

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

Công trình nghiên cứu mang giá trị tham khảo chuyên môn cao cho 4 nhóm đối tượng chính:

  • Kỹ sư vận hành và điều độ hệ thống điện: Nắm vững phương pháp tính toán giới hạn ổn định, phân tích chế độ quá độ và xác định thời gian cắt sự cố tối ưu khi lưới điện tích hợp nguồn quang điện tỷ lệ cao trên 50%.
  • Chuyên gia tư vấn thiết kế dự án năng lượng tái tạo: Ứng dụng mô hình toán học chi tiết của 363 tổ hợp PV 500 kVA và cấu hình biến tần để tối ưu hóa thiết kế kỹ thuật trạm điện mặt trời hòa lưới cấp 110 kV - 230 kV.
  • Giảng viên, học viên cao học và sinh viên ngành Kỹ thuật điện: Sử dụng tài liệu như một nguồn tham khảo chuẩn mực về phương pháp giải phương trình trạng thái phi tuyến, thuật toán tích phân số và kỹ thuật mô phỏng hệ thống điện trên phần mềm Matlab.
  • Cơ quan quản lý nhà nước và hoạch định chính sách: Khai thác các căn cứ khoa học và số liệu kiểm chứng thực tế để ban hành các quy định đấu nối, quy hoạch nguồn điện tái tạo phù hợp với định hướng phát triển năng lượng bền vững.

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

  • Tỷ lệ thâm nhập điện mặt trời 51% ảnh hưởng như thế nào đến quán tính hệ thống? Nguồn pin mặt trời kết nối qua biến tần không có các bộ phận quay cơ học, làm giảm tổng quán tính quay của toàn hệ thống. Tuy nhiên, tính toán cho thấy các máy phát đồng bộ còn lại kết hợp bộ điều khiển biến tần phản ứng nhanh vẫn duy trì tốt sự ổn định của hệ thống.

  • Tại sao nghiên cứu lại tập trung vào sự cố ngắn mạch 3 pha? Dù ngắn mạch 3 pha chỉ chiếm 5-10% tổng số sự cố thực tế so với mức 70-90% của ngắn mạch 1 pha chạm đất, đây là dạng sự cố nghiêm trọng nhất làm công suất điện phát sụt giảm về 0, gây nguy cơ mất đồng bộ cao nhất cho toàn lưới điện.

  • Phương pháp dao động bé đánh giá ổn định tĩnh như thế nào? Phương pháp tiến hành tuyến tính hóa hệ phương trình vi phân quanh điểm vận hành xác lập để trích xuất ma trận trạng thái A. Khi tất cả các giá trị riêng đều có phần thực âm, hệ thống được khẳng định đạt trạng thái ổn định tĩnh bền vững trước các nhiễu loạn nhỏ.

  • Biến thiên bức xạ mặt trời ảnh hưởng ra sao đến chế độ làm việc của lưới? Khi bức xạ giảm từ 1000 W/m² xuống 600 W/m², dòng ngắn mạch và công suất tác dụng sụt giảm tỷ lệ thuận. Hệ thống điều khiển nghịch lưu cần can thiệp trong thời gian dưới 1,5 giây để bù công suất phản kháng và ổn định điện áp thanh cái trong dải cho phép ±5%.

  • Tiêu chuẩn cân bằng diện tích đóng vai trò gì trong phân tích ổn định động? Tiêu chuẩn cân bằng diện tích cho phép so sánh động năng tích lũy trong quá trình tăng tốc rotor với năng lượng tiêu tán trong quá trình giảm tốc sau sự cố. Nếu diện tích giảm tốc lớn hơn hoặc bằng diện tích tăng tốc, hệ thống đảm bảo duy trì ổn định động mà không bị mất đồng bộ.

Kết luận

  • Nghiên cứu đã xây dựng thành công mô hình toán học tích hợp 363 tổ hợp pin năng lượng mặt trời 500 kVA (tổng công suất xấp xỉ 163,35 MW) vào hệ thống chuẩn IEEE 3 máy 9 bus.
  • Khẳng định tính khả thi kỹ thuật của việc thay thế hoàn toàn một nhà máy điện truyền thống bằng nguồn quang điện công suất lớn với tỷ lệ thâm nhập 51%.
  • Xác định chính xác giới hạn thời gian giải trừ ngắn mạch 3 pha dưới 0,15 giây nhằm bảo đảm ổn định động cho hệ thống truyền tải điện 230 kV.
  • Chứng minh khả năng giữ vững điện áp và tần số của lưới điện khi xảy ra biến thiên bức xạ mặt trời đột ngột 40%.
  • Đóng góp cơ sở lý thuyết và thực tiễn vững chắc cho công tác quy hoạch, điều độ và vận hành hệ thống điện tích hợp năng lượng tái tạo tại Việt Nam.

Hướng nghiên cứu tiếp theo trong giai đoạn 2025-2030 là mở rộng phân tích trên hệ thống điện thực tế quy mô lớn IEEE 39 bus và tích hợp thêm hệ thống pin lưu trữ năng lượng BESS. Độc giả, kỹ sư và các nhà nghiên cứu quan tâm có thể khai thác toàn văn tài liệu để phục vụ công tác nghiên cứu chuyên sâu và ứng dụng vào thực tiễn quản lý vận hành lưới điện.