Tổng quan nghiên cứu

Sự phát triển mạnh mẽ của các nguồn năng lượng tái tạo đang định hình lại cấu trúc vận hành của hệ thống điện hiện đại. Tại Việt Nam, nguồn năng lượng gió sở hữu tiềm năng khai thác to lớn, đóng vai trò quan trọng trong chiến lược giảm phát thải carbon và bổ sung nguồn cung cho các khu vực xa lưới điện quốc gia. Tuy nhiên, việc tích hợp các trang trại điện gió quy mô lớn vào lưới truyền tải 110kV đặt ra nhiều thách thức nghiêm trọng về mặt kỹ thuật, đặc biệt là nguy cơ mất ổn định điện áp và sụt giảm độ dự trữ công suất phản kháng. Trong điều kiện vận hành bình thường, tổn thất điện năng trên các đường dây quá tải có thể vượt mức 7% đến 10%, đồng thời hiện tượng dao động công suất do tính bất định của tốc độ gió làm gia tăng áp lực lên hệ thống điều độ.

Mục tiêu cốt lõi của nghiên cứu là xây dựng công cụ tính toán và đánh giá toàn diện mức độ ổn định điện áp của lưới truyền tải khi có sự tham gia của các máy phát điện gió không đồng bộ nguồn kép (DFIG). Đề tài tập trung giải quyết bài toán xác định vị trí đấu nối tối ưu, định lượng mức độ xâm nhập công suất tác dụng và công suất phản kháng, đồng thời so sánh tác động ổn định điện áp giữa nguồn điện gió với các nhà máy thủy điện và nhiệt điện truyền thống. Phạm vi nghiên cứu được thực nghiệm trên mô hình lưới điện truyền tải 110kV điển hình trong khung thời gian từ tháng 2 năm 2014 đến tháng 6 năm 2014 tại Trường Đại học Bách Khoa – ĐHQG TP.HCM. Kết quả nghiên cứu mang lại giá trị thực tiễn cao, hỗ trợ nâng biên độ truyền tải công suất thêm 15% đến 25% mà vẫn đảm bảo điện áp các thanh cái duy trì trong ngưỡng cho phép dao động trong dải ±10% điện áp định mứ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 được xây dựng trên nền tảng lý thuyết ổn định điện áp tĩnh và lý thuyết phân tích giá trị riêng (modal analysis) của hệ thống điện. Hai mô hình nền tảng được áp dụng xuyên suốt bao gồm: mô hình toán học máy phát điện không đồng bộ nguồn kép (DFIG) trong hệ tọa độ đồng bộ dq tựa theo vectơ điện áp lưới và mô hình phân bố công suất cải biên dựa trên thuật toán Newton-Raphson. Trong mô hình này, nút máy phát gió không được đối xử như nút PV thông thường do giới hạn công suất của bộ biến đổi nguồn back-to-back PWM không đủ khả năng duy trì cố định điện áp cực máy. Thay vào đó, nút gió được mô hình hóa dưới dạng nút tải PQ đặc thù, trong đó công suất tác dụng $P_e$ và công suất phản kháng $Q_e$ được cập nhật liên tục sau mỗi bước lặp thông qua kiểm tra giới hạn dòng điện rotor $I_{2d}$ và $I_{2q}$.

Các khái niệm then chốt cấu thành khung phân tích gồm: ma trận Jacobian rút gọn $J_R = J_{QV} - J_Q J_P^{-1} J_{PV}$, trị riêng nhỏ nhất $\lambda_{\min}$ đặc trưng cho khoảng cách đến điểm sụp đổ điện áp, vectơ riêng phải và trái biểu diễn dạng thức dao động (mode), và hệ số tham gia của các nút phụ tải. Tiêu chuẩn ổn định xác định rằng hệ thống giữ được ổn định điện áp khi và chỉ khi tất cả các trị riêng của ma trận $J_R$ đều mang giá trị dương ($\lambda_i > 0$). Khi $\lambda_{\min}$ tiến dần về 0, hệ thống rơi vào trạng thái nguy hiểm và có thể dẫn tới sụp đổ điện áp cục bộ hoặc toàn phần.

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

Nguồn dữ liệu phục vụ nghiên cứu bao gồm các thông số kỹ thuật thực tế của cụm tuabin gió DFIG công suất lớn, kết hợp với bộ thông số đường dây, máy biến áp và phụ tải của mạng điện truyền tải 110kV. Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu đã thiết lập và khảo sát toàn diện 11 kịch bản vận hành chi tiết (từ Trường hợp 0 đến Trường hợp 11) trên cấu trúc lưới đa nút. Phương pháp chọn mẫu có chủ đích được áp dụng để khảo sát các nút đại diện gồm nút 12, nút 10 và nút 3 nhằm phản ánh đầy đủ các vị trí đấu nối có khoảng cách điện và đặc tính ngắn mạch khác nhau, cùng 2 trường hợp thay thế nguồn phát tại nút 5 và nút 2.

Lý do lựa chọn phương pháp phân tích modal Q-V thay vì phương pháp mô phỏng động theo miền thời gian xuất phát từ yêu cầu tính toán nhanh, trực quan và độ chính xác cao trong bài toán đánh giá tĩnh. Mô phỏng động đòi hỏi thời gian giải hệ phương trình vi phân phi tuyến kéo dài nhiều phút và tiêu tốn tài nguyên tính toán lớn, trong khi phân tích modal cho phép trích xuất trực tiếp giá trị $\lambda_{\min}$ và hệ số tham gia của từng nút chỉ qua một chu kỳ tính toán ma trận. Toàn bộ thuật toán được lập trình tự động hóa trên nền tảng phần mềm MATLAB, hoàn thành trong giai đoạn 5 tháng nghiên cứu nghiêm ngặt.

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 11 kịch bản vận hành đã mang lại 4 phát hiện kỹ thuật mang tính đột phá:

Thứ nhất, vị trí kết nối trang trại gió quyết định trực tiếp đến giới hạn truyền tải công suất cực đại của hệ thống. Khi kết nối nguồn gió tại nút 12 (vị trí có liên kết lưới mạnh), mức độ tăng tải tối đa của hệ thống đạt mức cao hơn khoảng 22.4% so với kịch bản cơ sở chưa kết nối gió (Trường hợp 0), với trị riêng nhỏ nhất $\lambda_{\min}$ duy trì ở mức an toàn 0.185. Ngược lại, khi kết nối tại nút 3 (vị trí xa nguồn), mức tăng tải chỉ đạt khoảng 8.5% trước khi xuất hiện hiện tượng sụt áp nghiêm trọng.

Thứ hai, mức độ xâm nhập công suất tác dụng của gió ($P_{wind}$) có mối quan hệ phi tuyến với độ dự trữ công suất phản kháng. Khi nâng công suất phát của gió lên mức 30% đến 40% công suất tải toàn hệ thống, lượng công suất phản kháng $Q$ mà máy phát DFIG tiêu thụ từ lưới tăng nhanh từ 15% đến 28%, kéo theo sự suy giảm của $\lambda_{\min}$ với tốc độ 12% mỗi chu kỳ tăng tải nếu không có biện pháp bù cưỡng bức.

Thứ ba, việc thay thế nguồn phát truyền thống (thủy điện hoặc nhiệt điện tại nút 5 và nút 2) bằng nguồn điện gió làm suy giảm biên độ ổn định điện áp tổng thể từ 9.2% đến 14.6%. Nguyên nhân xuất phát từ việc máy phát DFIG thiếu khả năng hỗ trợ quán tính cơ học và dự trữ công suất phản kháng nội tại thấp hơn các máy phát đồng bộ cực lồi truyền thống.

Thứ tư, hiệu chỉnh hệ số công suất $\cos\varphi$ của máy phát gió từ 0.90 trễ sang 0.98 sớm giúp cải thiện trị riêng $\lambda_{\min}$ lên thêm 25.3%, ngăn ngừa hiện tượng mất ổn định điện áp khi phụ tải đỉnh gia tăng đột biến.

Thảo luận kết quả

Hiện tượng sụp đổ điện áp khi tích hợp điện gió bắt nguồn từ sự mất cân đối công suất phản kháng cục bộ. Khi các tuabin gió tăng công suất phát tác dụng, dòng điện trên các xuất tuyến 110kV tăng cao, làm tổn thất công suất phản kháng dòng điện trên đường dây ($I^2X$) tăng vọt. Bộ điều áp dưới tải của máy biến áp (OLTC) cần từ 2 đến 4 phút để chuyển nấc phân áp, nhưng mỗi lần chuyển nấc trong điều kiện điện áp thấp lại làm dòng điện sơ cấp tăng thêm, đẩy nhanh quá trình bão hòa từ và dẫn đến sụp đổ điện áp dây chuyền.

Các kết quả nghiên cứu có thể được biểu diễn trực quan thông qua bảng so sánh mức độ tăng tải cực đại giữa 11 trường hợp và đồ thị phổ trị riêng modal $\lambda_i$. Việc đối chiếu đường cong phân bố độ nhạy V-Q tại các thanh cái trên biểu đồ cột minh họa rõ nét sự dịch chuyển của điểm làm việc tới giới hạn sụp đổ. So sánh với các công bố quốc tế về máy phát phân tán, nghiên cứu này đã làm rõ cơ chế tương tác đa biến giữa công suất tác dụng và phản kháng của DFIG, cung cấp bằng chứng thực nghiệm giá trị cho bài toán vận hành lưới điện truyền tải khu vực.

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

Dựa trên các kết luận tính toán khoa học, 4 giải pháp trọng tâm được đề xuất nhằm đảm bảo an toàn vận hành hệ thống điện:

  1. Quy hoạch và tối ưu hóa vị trí đấu nối nguồn gió: Tập đoàn Điện lực Việt Nam và các đơn vị tư vấn thiết kế cần áp dụng thuật toán phân tích modal để xác định chỉ số nhạy cảm điện áp trước khi phê duyệt điểm đấu nối, ưu tiên các nút có liên kết mạch vòng mạnh nhằm nâng mức tăng tải khả dụng lên trên 20%, áp dụng ngay trong giai đoạn thẩm định dự án từ 6 đến 12 tháng.
  2. Kiểm soát nghiêm ngặt hệ số công suất phát: Các chủ đầu tư nhà máy điện gió bắt buộc phải cài đặt cấu hình bộ biến đổi DFIG vận hành trong dải $\cos\varphi$ từ 0.95 trễ đến 0.98 sớm, hoặc trang bị bổ sung hệ thống bù tĩnh FACTS/STATCOM có dung lượng bù tối thiểu bằng 20% công suất đặt của nhà máy, duy trì liên tục trong suốt vòng đời dự án 20 năm.
  3. Hoàn thiện quy chuẩn kỹ thuật hòa lưới (Grid Code): Cục Điều tiết Điện lực cần ban hành quy định bắt buộc về khả năng vượt qua sự cố điện áp thấp (LVRT) và đặt ngưỡng giới hạn dự trữ ổn định điện áp $\lambda_{\min} \ge 0.10$ đối với tất cả các trang trại gió đấu nối cấp điện áp 110kV trở lên, hoàn thành lộ trình pháp lý trong vòng 18 tháng.
  4. Tích hợp module cảnh báo ổn định tĩnh thời gian thực: Trung tâm Điều độ Hệ thống điện Quốc gia cần triển khai thuật toán ma trận Jacobian rút gọn vào hệ thống SCADA/EMS, quét chu kỳ 15 đến 30 giây để phát hiện sớm các nút có độ nhạy V-Q âm, hoàn thành lắp đặt thử nghiệm trong vòng 2 năm.

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

Nội dung luận văn mang tính ứng dụng chuyên sâu, là tài liệu tham khảo thiết thực cho 4 nhóm đối tượng:

  • Kỹ sư vận hành và điều độ hệ thống điện: Nắm vững phương pháp tính ma trận modal để đánh giá nhanh trạng thái an toàn điện áp, đưa ra quyết định sa thải phụ tải hoặc huy động nguồn bù chính xác trong các tình huống khẩn cấp trên lưới 110kV.
  • Kỹ sư tư vấn thiết kế và quy hoạch nguồn điện: Khai thác mô hình toán học DFIG cải biên và lưu đồ giải thuật Newton-Raphson để tính toán trào lưu công suất, lựa chọn vị trí đấu nối trang trại gió tối ưu về mặt kinh tế - kỹ thuật.
  • Chủ đầu tư và đơn vị quản lý dự án năng lượng tái tạo: Hiểu rõ các ràng buộc kỹ thuật về công suất phản kháng và tiêu chuẩn điện áp tại điểm đấu nối, từ đó lựa chọn công nghệ tuabin và thiết bị bù phù hợp, giảm thiểu rủi ro bị từ chối phát công suất lên lưới.
  • 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 làm tài liệu nghiên cứu chuyên đề về ổn định hệ thống điện, phát triển các thuật toán điều khiển tối ưu cho nguồn phân tán trên nền tảng phần mềm MATLAB.

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

Tại sao máy phát điện gió DFIG lại được mô hình hóa thành nút PQ thay vì nút PV trong bài toán trào lưu công suất? Do công suất của bộ biến đổi điện tử công suất PWM bị giới hạn dòng điện cực đại, DFIG không thể cung cấp lượng công suất phản kháng vô hạn để giữ cố định điện áp đầu cực như máy phát đồng bộ. Nút DFIG được mô hình hóa thành nút PQ cải biên với $P_e$ và $Q_e$ thay đổi linh hoạt sau mỗi vòng lặp dựa trên đặc tính trượt của rotor và điện áp nút thực tế.

Phương pháp phân tích modal Q-V có ưu điểm gì vượt trội so với phân tích ổn định động? Phân tích modal Q-V dựa trên việc tuyến tính hóa ma trận Jacobian rút gọn, giúp tính toán tức thời các trị riêng $\lambda_i$ mà không cần giải hệ phương trình vi phân phức tạp theo thời gian. Phương pháp này chỉ ra trực quan khoảng cách đến điểm mất ổn định thông qua biên độ $\lambda_{\min}$ và định vị chính xác các nút nhạy cảm nhất trên lưới 110kV.

Hiện tượng sụp đổ điện áp trên lưới truyền tải khi tích hợp điện gió xảy ra do nguyên nhân chủ yếu nào? Nguyên nhân cốt lõi là sự thiếu hụt công suất phản kháng nghiêm trọng khi hệ thống truyền tải công suất tác dụng lớn đi xa. Quá trình này làm gia tăng tổn thất phản kháng $I^2X$ trên đường dây, kết hợp với đặc tính hút công suất phản kháng của máy phát không đồng bộ và thời gian đáp ứng trễ từ 2 đến 4 phút của bộ điều áp OLTC.

Mức độ thâm nhập của điện gió ảnh hưởng thế nào đến khả năng chịu tải của mạng điện? Khi điện gió kết nối tại nút tối ưu có dung lượng ngắn mạch lớn, mức tăng tải cực đại có thể nâng thêm từ 15% đến 25%. Tuy nhiên, nếu tỷ lệ xâm nhập vượt ngưỡng cho phép mà không có nguồn bù hỗ trợ, trị riêng $\lambda_{\min}$ sẽ suy giảm nhanh về 0, khiến hệ thống mất ổn định điện áp ngay ở mức phụ tải thấp.

Việc thay thế nhà máy nhiệt điện hoặc thủy điện bằng trang trại gió gây ra thách thức gì? Các nhà máy truyền thống sở hữu máy phát đồng bộ có khả năng duy trì điện áp mạnh mẽ và quán tính cơ học cao. Khi thay thế bằng máy phát gió DFIG, hệ thống bị suy giảm mức dự trữ ổn định điện áp từ 9% đến 14%, đòi hỏi phải bổ sung các thiết bị điều khiển điện áp nhanh để bù đắp sự thiếu hụt này.

Kết luận

  • Xây dựng thành công thuật toán và chương trình tính toán phân bố công suất cải biên cho lưới truyền tải 110kV tích hợp máy phát gió DFIG trên nền tảng MATLAB.
  • Chứng minh tính hiệu quả vượt trội của phương pháp phân tích modal Q-V trong việc xác định định lượng ranh giới ổn định điện áp thông qua trị riêng nhỏ nhất $\lambda_{\min}$.
  • Khảo sát chi tiết 11 kịch bản vận hành, xác định chính xác vị trí đấu nối tối ưu giúp nâng mức tăng tải của lưới điện lên thêm 22.4% mà vẫn duy trì an toàn điện áp.
  • Làm rõ tác động tiêu cực của việc gia tăng công suất phát gió và suy giảm hệ số $\cos\varphi$ đến mức độ sụt giảm dự trữ công suất phản kháng toàn hệ thống.
  • Đề xuất các giải pháp kỹ thuật khả thi về quy hoạch, điều khiển bộ biến đổi và bổ sung thiết bị bù tĩnh phục vụ công tác thiết kế và vận hành thực tế.

Trong giai đoạn tiếp theo từ năm 2024 đến 2026, hướng phát triển trọng tâm là mở rộng thuật toán sang phân tích ổn định quá độ thời gian thực kết hợp hệ thống lưu trữ năng lượng pin (BESS). Các kỹ sư và nhà nghiên cứu quan tâm có thể áp dụng ngay khung phương pháp luận này để nâng cao độ tin cậy và an toàn cho các dự án năng lượng tái tạo đang triển khai.