Tổng quan nghiên cứu

Việt Nam nằm trong khu vực khí hậu nhiệt đới gió mùa nóng ẩm, là môi trường đặc biệt thuận lợi cho sự hình thành và phát triển của các hiện tượng dông sét cường độ mạnh. Theo thống kê thực tế từ ngành điện lực, hàng năm hệ thống điện ghi nhận hàng ngàn vụ sự cố, trong đó hơn 60% trường hợp bắt nguồn từ tác động trực tiếp hoặc lan truyền của dông sét. Sự cố nghiêm trọng xảy ra vào lúc 19h41 ngày 30 tháng 05 năm 2016 là một minh chứng điển hình, khi sét đánh làm bật đồng thời hai đường dây 500kV Hà Tĩnh – Đà Nẵng và Vũng Áng – Đà Nẵng, gây mất liên kết lưới điện quốc gia và làm gián đoạn cung cấp điện diện rộng tại khu vực miền Trung và miền Nam trong khoảng 50 phút.

Trên lưới điện phân phối, quá điện áp do xung sét lan truyền là mối đe dọa nguy hiểm hàng đầu đối với cách điện của máy biến áp và các thiết bị đóng cắt đắt tiền. Việc đo đạc, thử nghiệm thực địa phản ứng của thiết bị chống sét van trung áp dưới tác động của xung sét thực tế gặp nhiều rào cản lớn về kỹ thuật an toàn và chi phí vận hành. Nhằm giải quyết triệt để thách thức này, nghiên cứu tập trung xây dựng mô hình chống sét van oxit kim loại trung áp (MVLA) dựa trên cấu trúc P-K (Pinceti - Giannettoni) trong môi trường Matlab/Simulink và đề xuất phương án bố trí bảo vệ tối ưu cho trạm biến áp phân phối 22kV tại Huyện Phú Tân, Tỉnh Cà Mau vào tháng 03 năm 2018. Kết quả nghiên cứu mang lại ý nghĩa thực tiễn to lớn khi cung cấp công cụ mô phỏng có độ chính xác cao với sai số điện áp dư dưới 5,1%, giúp các công ty điện lực tối ưu hóa chi phí đầu tư bảo vệ và nâng cao chỉ số độ tin cậy cung cấp đ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 được xây dựng dựa trên sự kết hợp của lý thuyết vi cấu trúc gốm bán dẫn oxit kim loại (ZnO MOV) và lý thuyết sóng quá độ lan truyền trong hệ thống điện. Bản chất của biến trở ZnO gồm ma trận các hạt oxit kẽm dẫn điện tốt với đường kính từ 15 đến 100 micromet, bao quanh bởi lớp biên oxit kim loại có điện trở cao. Mỗi điểm tiếp giáp hạt tạo thành một vi biến trở với điện áp rào khoảng 2,0V đến 3,5V. Khi điện áp lưới ở mức bình thường, thiết bị có trở kháng rất cao gần như hở mạch với dòng rò chỉ từ 0 đến 2A; khi xuất hiện sóng sét đột biến, trở kháng giảm tức thì về mức rất thấp để xả dòng xung xuống đất.

Khung lý thuyết của đề tài tập trung vào 4 khái niệm và mô hình kỹ thuật trọng tâm:

  1. Mức chịu đựng xung cơ bản (BIL) của thiết bị phân phối nhằm xác định ngưỡng phối hợp cách điện an toàn.
  2. Mô hình sóng xung sét tiêu chuẩn dạng 8/20 micro giây với phương trình phân rã hàm mũ kép.
  3. Mô hình chống sét van phụ thuộc tần số của IEEE và mô hình giản lược P-K.
  4. Điện áp dư danh định ứng với các mức biên độ dòng phóng 5kA, 10kA và 20kA.

Mô hình P-K đơn giản hóa cấu trúc IEEE bằng cách sử dụng hai điện trở phi tuyến A0, A1 ngăn cách bởi một cuộn cảm L1. Giá trị của L1 được tính toán trực tiếp từ điện áp định mức và tỷ số điện áp dư ở các dạng sóng xung khác nhau mà không cần bóc tách các thông số vật lý phức tạp bên trong thiết bị.

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

Nguồn dữ liệu đầu vào của luận văn được thu thập từ catalogue kỹ thuật chính hãng của hai dòng chống sét van trung áp 18kV phổ biến trên thị trường: loại AZG2 do Hãng Cooper Power Systems sản xuất và loại EVP do Hãng Ohio Brass sản xuất. Đồng thời, đề tài thu thập số liệu vận hành thực tế về sơ đồ mặt bằng và cấu hình đấu nối của các trạm biến áp phân phối đại diện tại Huyện Phú Tân, Tỉnh Cà Mau.

Phương pháp phân tích chủ đạo là mô hình hóa toán học kết hợp mô phỏng số trên nền tảng Matlab/Simulink. Tác giả sử dụng công cụ Curve Fitting Toolbox để tìm hàm xấp xỉ các thông số phi tuyến của bộ phát xung dòng điện và lập trình giao diện đóng gói Mask Subsystem để tự động hóa tính toán thông số mô hình.

Nghiên cứu lựa chọn cỡ mẫu gồm 2 thiết bị chống sét van từ 2 nhà sản xuất uy tín toàn cầu và tiến hành kiểm tra trên 3 cấp biên độ dòng sét thử nghiệm tiêu chuẩn là 5kA, 10kA và 20kA, cùng với 2 mô hình trạm biến áp (trạm 1 máy biến áp và trạm 2 máy biến áp). Phương pháp mô phỏng số được lựa chọn vì đây là giải pháp khả thi duy nhất cho phép tái hiện chính xác các quá trình quá độ diễn ra trong thời gian micro giây mà không làm gián đoạn hoạt động của lưới điện thực. Toàn bộ quá trình thu thập số liệu, xây dựng thuật toán và kiểm chứng mô hình được hoàn thành trong giai đoạn từ cuối năm 2017 đến tháng 03 năm 2018.

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 định lượng đã mang lại 4 phát hiện quan trọng:

  1. Bộ nguồn phát xung dòng sét dạng 8/20 micro giây được thiết lập thành công trên Matlab, đáp ứng chính xác dạng sóng hàm mũ kép tiêu chuẩn ở mọi cấp dòng điện thử nghiệm từ 5kA, 10kA đến 20kA với thời gian đạt đỉnh 8 micro giây và thời gian suy giảm một nửa giá trị đỉnh là 20 micro giây.
  2. Mô hình P-K đạt độ chính xác cao vượt trội so với mô hình mặc định trong thư viện SimPowerSystems của Matlab. Đối với chống sét van AZG2 của Hãng Cooper ở điện áp định mức 18kV, mô hình P-K ghi nhận sai số điện áp dư chỉ từ 2,3% đến 5,1%, trong khi mô hình mặc định của Matlab có sai số dao động từ 4,5% đến 8,3%.
  3. Đối với chống sét van EVP của Hãng Ohio Brass, mô hình P-K thể hiện khả năng thích ứng xuất sắc khi sai số điện áp dư ở xung dòng 10kA đạt mức ấn tượng chỉ 0,2%, so với mức sai số 6,8% của mô hình Matlab. Tính chung trên toàn dải đo 5kA đến 20kA, sai số của mô hình P-K luôn duy trì trong khoảng 0,2% đến 7,3%, hoàn toàn đáp ứng yêu cầu kỹ thuật chuẩn với sai số dưới 10%.
  4. Khoảng cách lắp đặt chống sét van có ảnh hưởng mang tính quyết định đến mức quá điện áp tác động lên đầu cực máy biến áp. Khi chống sét van được lắp trực tiếp tại đầu cực máy biến áp (khoảng cách 0 mét), điện áp xung được ghìm ở mức an toàn nhất. Khi tăng chiều dài dây nối lên 4 mét và 12 mét trong cấu hình trạm cấp điện cho nhiều máy biến áp, biên độ điện áp tại đầu cực máy biến áp tăng vọt đáng kể do ảnh hưởng của điện cảm trên đường dây đấu nối.

Thảo luận kết quả

Nguyên nhân chính giúp mô hình P-K đạt độ chính xác vượt bậc nằm ở cấu trúc mạch tương đương phụ thuộc tần số. Bằng việc phân tách thành hai phần tử phi tuyến A0 (đặc trưng cho quá trình quá độ xung dốc) và A1 (đặc trưng cho xung sét thông thường) liên kết qua độ tự cảm L1, mô hình P-K đã tái lập chuẩn xác hiện tượng trễ pha và cảm kháng ký sinh của gốm oxit kim loại khi dòng điện có độ dốc biến thiên lớn. Ngược lại, mô hình mặc định trong SimPowerSystems của Matlab chỉ sử dụng phương trình hàm mũ 3 đoạn thuần túy không xét đến đáp ứng tần số cao, dẫn đến độ lệch lớn khi dòng xung tăng nhanh.

So với các mô hình kinh điển từ Đại học Kỹ thuật Quốc gia Athens hay Đại học Queensland, mô hình P-K trong luận văn thể hiện ưu thế vượt trội về tính ứng dụng kỹ thuật thực tế. Các mô hình quốc tế trước đây đòi hỏi nhiều thông số vật lý hình học bên trong đĩa MOV như chiều cao cột, tiết diện tiếp giáp và điện dung nội tại – những thông số vốn thuộc bí mật công nghệ của nhà sản xuất và không bao giờ xuất hiện trong tài liệu thương mại. Mô hình P-K chỉ yêu cầu 3 thông số cơ bản có sẵn trong catalogue: điện áp định mức, điện áp dư ở xung 8/20 micro giây 10kA và điện áp dư ở xung dốc 10kA.

Trong thực tế tính toán kỹ thuật, toàn bộ kết quả này có thể được biểu diễn trực quan thông qua các biểu đồ đáp ứng thời gian – điện áp dạng sóng kép và bảng tổng hợp sai số phần trăm định lượng. Nhờ đó, kỹ sư vận hành có thể nhanh chóng so sánh biên độ điện áp dư cực đại với mức chịu đựng xung sét của máy biến áp để ra quyết định lắp đặt bảo vệ chính xác.

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

Nhằm nâng cao hiệu quả chống sét cho hệ thống phân phối điện thực tế, luận văn đưa ra 4 giải pháp trọng tâm:

  1. Ứng dụng ngay công cụ mô phỏng chống sét van P-K trên Matlab/Simulink vào công tác tính toán quá điện áp định kỳ, đặt mục tiêu duy trì sai số mô phỏng thiết kế dưới 5,0% trong toàn bộ các dự án cải tạo lưới điện trung áp bắt đầu từ quý 3 năm 2018. Chủ thể thực hiện là Phòng Kỹ thuật và Ban Quản lý Dự án thuộc Công ty Điện lực Cà Mau.
  2. Thi công rút ngắn chiều dài dây dẫn kết nối từ chống sét van đến đầu cực máy biến áp phân phối về khoảng cách dưới 1,0 mét (lý tưởng nhất là gá lắp trực tiếp 0 mét tại cọc sứ máy biến áp), nhằm giảm điện áp quá độ tác động lên cách điện máy biến áp từ 15% đến 25% trong vòng 6 tháng triển khai. Chủ thể thực hiện là Đội Quản lý Vận hành Lưới điện Huyện Phú Tân.
  3. Chuẩn hóa quy trình nghiệm thu vật tư thiết bị chống sét van trung áp 18kV và 24kV, bắt buộc nhà thầu cung cấp đầy đủ thông số điện áp dư ở cả hai chế độ xung 8/20 micro giây và xung dốc 10kA trong hồ sơ kỹ thuật trước quý 4 năm 2018. Chủ thể thực hiện là Phòng Vật tư và Ban Đấu thầu Tổng công ty Điện lực miền Nam.
  4. Tổ chức các khóa đào tạo chuyển giao công nghệ mô phỏng lưới điện trên phần mềm Matlab cho 100% kỹ sư điều độ và kỹ thuật viên lưới điện phân phối trong kế hoạch đào tạo chuyên môn năm 2019. Chủ thể thực hiện là Trung tâm Đào tạo ngành điện phối hợp cùng các chuyên gia thuộc Trường Đại học Sư phạm Kỹ thuật TP.HCM.

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

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

  1. Kỹ sư thiết kế và quy hoạch lưới điện tại các công ty điện lực: Nắm vững công cụ mô phỏng để tính toán tọa độ lắp đặt và khoảng cách bảo vệ an toàn của chống sét van, áp dụng trực tiếp cho các dự án xây dựng mới hoặc nâng cấp trạm biến áp phân phối 22kV.
  2. 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 làm tài liệu học tập và nghiên cứu chuyên khảo về mô hình hóa quá điện áp khí quyển, kỹ thuật phân tích phần tử phi tuyến và phương pháp ứng dụng Matlab/Simulink trong hệ thống điện.
  3. Đơn vị sản xuất và cung ứng thiết bị chống sét oxit kim loại: Hiểu rõ yêu cầu về dữ liệu đầu vào của các mô hình tính toán hiện đại để chuẩn hóa và công bố các thông số catalogue phù hợp với nhu cầu thiết kế thực tế của khách hàng.
  4. Đơn vị tư vấn và kiểm định an toàn điện: Vận dụng mô hình P-K để đánh giá mức độ suy giảm chất lượng và dự báo rủi ro phóng điện chọc thủng cách điện của các trạm biến áp đang vận hành tại các khu vực có mật độ sét cao.

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

Tại sao mô hình P-K lại có độ chính xác cao hơn mô hình mặc định trong Matlab SimPowerSystems?
Mô hình P-K tích hợp cấu trúc mạch tương đương phụ thuộc tần số gồm hai nhánh phi tuyến A0 và A1 liên kết qua điện cảm L1. Cấu trúc này phản ánh chính xác đặc tính cảm kháng ký sinh của gốm oxit kẽm khi chịu xung dòng điện có tốc độ tăng cực nhanh, điều mà mô hình hàm mũ tĩnh của Matlab không thể mô phỏng được.

Dạng sóng xung sét 8/20 micro giây mang ý nghĩa kỹ thuật gì trong nghiên cứu?
Dạng sóng 8/20 micro giây là chuẩn quy ước quốc tế mô phỏng dòng sét cảm ứng lan truyền trên đường dây trung áp, với thời gian tăng trưởng điện áp từ 0 lên đỉnh trong 8 micro giây và thời gian suy giảm xuống 50% biên độ đỉnh trong 20 micro giây, đại diện cho tác động sét phổ biến nhất trên lưới điện phân phối.

Khoảng cách lắp đặt giữa chống sét van và máy biến áp có tác động thực tế như thế nào?
Khi khoảng cách dây dẫn từ chống sét van đến máy biến áp tăng từ 0 mét lên 4 mét hoặc 12 mét, độ sụt áp trên cảm kháng đường dây phụ thuộc vào độ dốc dòng điện ($L \cdot di/dt$) sẽ cộng thêm vào điện áp dư của van, làm tăng đáng kể nguy cơ đánh thủng cách điện máy biến áp.

Kỹ sư vận hành cần thu thập những thông số nào để thiết lập mô hình chống sét van P-K?
Người dùng chỉ cần trích xuất 3 thông số định mức từ catalogue của nhà sản xuất: điện áp danh định của chống sét van (kV), điện áp dư tại dòng sét 10kA sóng 8/20 micro giây (kV), và điện áp dư tại dòng sét xung dốc 10kA (kV) mà không cần can thiệp vào cấu trúc bên trong.

Mô hình nghiên cứu tại Huyện Phú Tân có thể áp dụng cho các khu vực khác không?
Mô hình có tính tổng quát hóa cao và hoàn toàn có thể áp dụng cho mọi lưới điện phân phối trung áp trên cả nước, chỉ cần thay đổi các thông số đường dây và thông số thiết bị chống sét van tương ứng với điều kiện vận hành tại địa phương đó.

Kết luận

  • Xây dựng thành công công cụ phát xung sét tiêu chuẩn 8/20 micro giây và mô hình chống sét van trung áp P-K trên nền tảng Matlab/Simulink với độ tin cậy vượt bậc.
  • Khẳng định mô hình P-K đạt độ chính xác cao với sai số điện áp dư chỉ từ 0,2% đến 5,1%, thấp hơn đáng kể so với mô hình mặc định của Matlab (sai số lên đến 8,5%).
  • Đơn giản hóa toàn diện quy trình mô phỏng khi chỉ sử dụng 3 tham số kỹ thuật sẵn có trong catalogue nhà sản xuất, giải quyết triệt để rào cản thông số vật lý ẩn.
  • Chứng minh bằng số liệu định lượng về việc lắp đặt chống sét van trực tiếp tại đầu cực máy biến áp là giải pháp kỹ thuật tối ưu nhất để triệt tiêu quá áp lan truyền.
  • Định hướng giai đoạn tiếp theo cần mở rộng mô phỏng cho các dạng xung sét phức tạp và các trạm biến áp thông minh; các đơn vị quản lý vận hành điện lực nên chủ động tích hợp công cụ này vào quy trình kiểm tra an toàn lưới điện ngay hôm nay.