Tổng quan nghiên cứu

Việt Nam nằm trong khu vực nhiệt đới gió mùa nóng ẩm, ghi nhận mật độ dông sét thuộc nhóm cao trên thế giới với khoảng 1,5 đến 2 triệu cú sét đánh mỗi năm. Trên hệ thống lưới điện phân phối trung áp từ 22kV đến 35kV, sự cố do quá điện áp cảm ứng và sét đánh trực tiếp chiếm tỷ lệ từ 40% đến 60% tổng số sự cố thoáng qua và vĩnh cửu. Cùng với sự phát triển của hệ thống tự động hóa lưới điện, các thiết bị điện tử công suất có độ dự trữ cách điện nhạy cảm như SCADA, biến tần, tụ bù Thyristor và bộ lọc sóng hài ngày càng phổ biến. Điều này đặt ra bài toán cấp bách về việc nâng cao tính chính xác trong tính toán phối hợp cách điện, khi các phương pháp giải tích truyền thống không còn phản ánh đúng bản chất phi tuyến phức tạp của thiết bị chống sét.

Nghiên cứu được thực hiện nhằm giải quyết mục tiêu mô hình hóa chính xác đáp ứng động của chống sét van oxide kim loại (Metal Oxide Varistor - MOV) không khe hở và xây dựng mô hình nguồn phát xung sét chuẩn $8/20\mu s$ trên nền tảng Matlab/Simulink. Phạm vi nghiên cứu tập trung vào thiết bị chống sét van cấp phân phối trung áp, kiểm chứng qua thông số thực tế của 5 nhà sản xuất thiết bị điện hàng đầu và khảo sát sự cố trên tuyến đường dây dài 15km tại Việt Nam. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, giúp các kỹ sư điện đánh giá chính xác độ suy giảm điện áp quá độ, giảm thiểu nguy cơ quá điện áp trên 30%, từ đó nâng cao độ tin cậy cung cấp điện, cải thiện rõ rệt các chỉ số SAIDI và SAIFI cho lưới điện phân phối.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Cơ sở lý thuyết của đề tài xây dựng trên nền tảng vật lý chất bán dẫn gốm đa tinh thể ZnO và các mô hình quá độ điện từ cao tần. Thành phần cơ bản của biến trở MOV chứa trên 90% mol ZnO kết hợp với khoảng 10% các oxit kim loại như $\text{Bi}_2\text{O}_3$, $\text{Sb}_2\text{O}_3$, $\text{CoO}$, $\text{MnO}$ và $\text{Cr}_2\text{O}_3$. Các hạt ZnO có kích thước hạt từ $15\mu m$ đến $100\mu m$ với điện trở suất thấp ($0,1 - 10,\Omega\cdot\text{cm}$), được bao bọc bởi lớp biên tiếp giáp giàu Bismuth có điện trở suất rất cao ($10^{10} - 10^{12},\Omega\cdot\text{cm}$). Cấu trúc này tạo thành chuỗi vi biến trở đối xứng với điện áp rào thế khoảng 2,0V đến 3,5V tại mỗi biên hạt, quyết định đặc tính V-I phi tuyến biểu diễn qua phương trình $I = K \cdot V^\alpha$, trong đó hệ số phi tuyến $\alpha$ đạt từ 10 đến trên 50.

Khung lý thuyết của luận văn phân tích sâu 5 mô hình toán học chống sét van: mô hình điện trở phi tuyến quy ước, mô hình điện trở kết hợp điện cảm phi tuyến, mô hình Schmidt, mô hình Mardira và mô hình phụ thuộc tần số IEEE (IEEE Working Group 3.4.11). Nghiên cứu định lượng hóa các khái niệm cốt lõi theo tiêu chuẩn IEC 60099-4 và ANSI C62.11 bao gồm: Điện áp vận hành liên tục cực đại (MCOV/Uc), Điện áp dư ($U_{\text{res}}$), Quá điện áp tạm thời (TOV), Mức chịu xung cơ bản (BIL) và 3 biên hạn bảo vệ tối thiểu gồm đầu sóng ($\text{PM}_1 \ge 20%$), xung sét ($\text{PM}_2 \ge 20%$) và xung đóng cắt ($\text{PM}_3 \ge 15%$).

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

Nghiên cứu sử dụng nguồn dữ liệu kỹ thuật từ catalog của 5 nhà chế tạo uy tín quốc tế gồm Cooper, GE Tranquell, Siemens, ELPRO và Ohio-Brass. Cỡ mẫu nghiên cứu gồm 19 tập dữ liệu kiểm nghiệm phóng điện độc lập ở các mức dòng xung danh định 1kA, 3kA, 5kA, 10kA và 20kA. Phương pháp chọn mẫu là chọn mẫu có mục đích (purposive sampling), tập trung vào các chống sét van cấp phân phối có dải điện áp định mức phổ biến từ 1kV đến 35kV trên lưới điện thực tế.

Lý do lựa chọn phương pháp mô hình hóa và mô phỏng số trên Matlab/Simulink thay vì phương pháp giải tích truyền thống là do quá trình phóng sét tạo ra các dao động điện từ có thời gian đầu sóng siêu ngắn cỡ nano giây đến micro giây. Phương pháp giải tích thông thường coi điện trở MOV là hằng số hoặc hàm tĩnh, dẫn đến bỏ qua hiệu ứng phụ thuộc tần số. Bằng cách thiết lập bộ phát xung dòng không chu kỳ kết hợp mạch lọc hai nhánh phi tuyến $A_0 - A_1$ qua thông số điện cảm $L_1$ và điện trở $R_1$, phương pháp mô phỏng số cho phép tính toán chính xác hiện tượng trễ V-I và sự vọt lố điện áp dư khi gặp sóng sét dốc đứng. Toàn bộ chu trình tính toán, tinh chỉnh tham số và mô phỏng được thực hiện liên tục trong thời gian nghiên cứu 2 năm học thuậ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 và phân tích đáp ứng động của chống sét van MOV mang lại các phát hiện cốt lõi:

  • Độ chính xác vượt trội của mô hình phụ thuộc tần số: Khi kiểm chứng trên cả 5 hãng sản xuất (ELPRO, GE, Siemens, Cooper, Ohio-Brass) tại dòng xung chuẩn $8/20\mu s$ với biên độ từ 1kA đến 20kA, mô hình phụ thuộc tần số IEEE cải tiến cho kết quả điện áp dư ($U_{\text{res}}$) có sai số dưới 5% so với dữ liệu thực nghiệm. Đặc biệt, mô hình điện cảm phi tuyến đối với các xung đầu dốc nhanh ghi nhận mức sai số cực thấp, chỉ khoảng 1%.
  • Hiện tượng tăng điện áp dư theo độ dốc đầu sóng: Kết quả mô phỏng chỉ ra rằng khi thời gian đạt đỉnh của dòng phóng giảm nhanh từ $8\mu s$ xuống $1,3\mu s$, điện áp dư trên hai đầu cực chống sét van tăng xấp xỉ 6%. Điều này khẳng định điện cảm nội tại và điện cảm dây nối đóng vai trò triệt tiêu hoặc làm trầm trọng thêm mức quá điện áp đặt lên cách điện thiết bị.
  • Hiệu quả phân tán dòng sét theo cự ly: Khảo sát thực nghiệm trên tuyến đường dây phân phối 15km cho thấy khi sét đánh tại vị trí cách nguồn 15km, 10km và 1km, việc lắp đặt chống sét van giúp dập tắt hoàn toàn dòng xung quá độ, kẹp điện áp đỉnh từ mức trên 200kV xuống ngưỡng an toàn từ 45kV đến 65kV, bảo vệ toàn diện cho hệ thống máy biến áp và phụ tải phía hạ lưu.

Thảo luận kết quả

Nguyên nhân tạo nên sự khác biệt giữa mô hình phụ thuộc tần số và mô hình tĩnh truyền thống xuất phát từ cơ chế phân bố mật độ hạt dẫn trong vùng tiếp giáp $p-n$ của các hạt ZnO. Khi xung dòng có tốc độ tăng dốc cao ($di/dt$ lớn), trở kháng của bộ lọc $R_1 - L_1$ trong mô hình IEEE tăng mạnh, buộc phần lớn dòng điện đi qua phần tử phi tuyến $A_0$ thay vì chia sẻ đồng đều cho $A_1$. Cơ chế này phản ánh chính xác hiện tượng vật lý thực tế khi điện thế rào $\Phi_L$ suy giảm nhanh và các hạt dẫn trôi tự do bị giới hạn bởi thời hằng tái hợp.

Dữ liệu mô phỏng được trực quan hóa thông qua hệ thống 26 biểu đồ sóng điện áp - dòng điện theo miền thời gian trên khối hiển thị Scope của Simulink, đồng thời tổng hợp qua các bảng so sánh sai số điện áp kẹp ở từng mức dòng xung 5kA, 10kA và 20kA. So sánh với các công bố khoa học trước đây vốn ghi nhận sai số tĩnh từ 15% đến 20%, phương pháp đề xuất trong luận văn đã tối ưu hóa tham số điện cảm $L_1 = 15d/n,(\mu\text{H})$ (với $d$ là chiều cao thân chống sét và $n$ là số cột đĩa song song), giúp khớp hoàn toàn đặc tuyến phóng điện catalog của nhà sản xuất. Điều này chứng minh rằng việc bố trí chống sét van càng gần thiết bị cần bảo vệ (khoảng cách dưới 10m) sẽ hạn chế tối đa điện cảm ký sinh trên dây dẫn, giữ biên hạn bảo vệ $\text{PM}_2$ luôn vượt mức an toàn 20%.

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

Dựa trên kết quả nghiên cứu và tính toán thực nghiệm, bốn giải pháp công nghệ trọng tâm được đề xuất nhằm nâng cao hiệu quả bảo vệ chống sét trên lưới điện phân phối:

  1. Chuẩn hóa quy trình lựa chọn chống sét van MOV theo mô hình phụ thuộc tần số: Yêu cầu các Công ty Điện lực và đơn vị tư vấn thiết kế áp dụng mô hình toán học hai nhánh phi tuyến $A_0 - A_1$ để tính toán dung lượng hấp thụ năng lượng (đạt mức 2,5 đến 5,0 kJ/kV theo MCOV), hoàn thành áp dụng đồng bộ trong quý 2 năm 2027 nhằm thay thế hoàn toàn phương pháp tra bảng tĩnh cổ điển.
  2. Tối ưu hóa khoảng cách lắp đặt và giảm chiều dài dây nối tiếp địa: Đội ngũ kỹ sư quản lý vận hành cần điều chỉnh vị trí lắp đặt chống sét van sao cho khoảng cách dây dẫn nối từ đầu cực thiết bị đến máy biến áp hoặc thiết bị đóng cắt không vượt quá 10m đến 15m. Giải pháp này giúp giảm sụt áp cảm ứng đầu cực tối thiểu 12%, thực hiện kiểm tra định kỳ 6 tháng/lần tại tất cả các trạm biến áp phân phối.
  3. Triển khai chương trình đo kiểm tra định kỳ dòng rò điện trở ($I_{\text{ref}}$): Đơn vị thí nghiệm điện lực áp dụng quy trình đo dòng rò trong dải từ 0,05mA đến 1,0mA trên mỗi $\text{cm}^2$ tiết diện đĩa MOV tại mức điện áp $1,25 \times \text{MCOV}$, mục tiêu cắt giảm 40% sự cố nổ vỡ chống sét do già hóa nhiệt trước mùa mưa bão năm 2027.
  4. Số hóa và tích hợp thư viện mô phỏng chống sét vào hệ thống quản lý SCADA/DMS: Trung tâm điều độ hệ thống điện phối hợp với các viện nghiên cứu chuyên ngành xây dựng cơ sở dữ liệu số hóa 100% các loại chống sét van đang vận hành, thực hiện phân tích tự động tọa độ sét đánh và biên hạn bảo vệ trong giai đoạn 2027 - 2028.

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

Công trình nghiên cứu mang giá trị thực tiễn và học thuật cao, là tài liệu tham khảo thiết thực cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế và quy hoạch hệ thống điện: Ứng dụng mô hình toán và thuật toán mô phỏng trong luận văn để tính toán chính xác mức phối hợp cách điện, xác định dung lượng dòng xung định mức (5kA hoặc 10kA) và kiểm toán biên hạn bảo vệ ($\text{PM}_1, \text{PM}_2, \text{PM}_3$) cho các trạm biến áp phân phối 22/0,4kV xây mới.
  • Chuyên gia quản lý và vận hành lưới điện tại các Tổng công ty Điện lực: Sử dụng các kịch bản mô phỏng lan truyền sóng sét ở các cự ly 1km, 10km, 15km để xây dựng phương án lắp đặt chống sét van tối ưu trên đường dây trên không, ngăn ngừa sự cố nhảy máy cắt Recloser do sét đánh cảm ứng.
  • Giảng viên, học viên cao học và sinh viên ngành Kỹ thuật điện: Khai thác cấu trúc mô hình khối Mask Editor trong Matlab/Simulink và mô hình nguồn phát xung dòng $8/20\mu s$ phục vụ công tác giảng dạy, nghiên cứu chuyên sâu về quá độ điện từ và kỹ thuật điện cao áp.
  • Các đơn vị sản xuất và thử nghiệm thiết bị điện cao thế: Tham chiếu quy trình vi cấu trúc gốm ZnO và phương pháp đo kiểm tra đặc tính $U_{\text{res}}$, TOV theo tiêu chuẩn IEC 60099-4 và ANSI C62.11 nhằm kiểm soát chất lượng xuất xưởng của các lô sản phẩm đĩa MOV.

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

1. Vì sao chống sét van ZnO không khe hở (MOV) thay thế hoàn toàn loại SiC có khe hở trước đây?

Chống sét van MOV sở hữu hệ số phi tuyến $\alpha$ vượt trội (từ 10 đến trên 50 so với 3 đến 5 của SiC), cho phép thiết bị vận hành liên tục mà không cần khe hở phóng điện. Nhờ đó, MOV triệt tiêu hoàn toàn dòng điện duy trì sau phóng sét, giảm thời gian đáp ứng xuống mức nano giây và tăng khả năng hấp thu năng lượng quá độ lên trên 2,5 kJ/kV.

2. Mô hình phụ thuộc tần số của IEEE thể hiện đặc tính động của chống sét van như thế nào?

Mô hình IEEE sử dụng hai nhánh điện trở phi tuyến $A_0$ và $A_1$ được phân tách bởi bộ lọc $R_1 - L_1$. Đối với xung sét có sườn trước dốc đứng (tần số cao), cảm kháng của $L_1$ tăng cao, ép dòng điện dồn qua nhánh $A_0$, làm tăng điện áp kẹp ở đầu cực lên khoảng 6%, phản ánh chính xác hiệu ứng trễ $V-I$ trong thực tế.

3. Chiều dài dây dẫn nối chống sét van ảnh hưởng như thế nào đến mức bảo vệ cách điện của thiết bị?

Điện cảm ký sinh của dây dẫn nối chống sét van ước tính khoảng $1\mu\text{H}/\text{m}$. Khi dòng sét có độ dốc cao $di/dt = 10\text{kA}/\mu\text{s}$ đi qua, mỗi mét dây dẫn tạo ra điện áp rơi cảm ứng $L \cdot (di/dt) = 10\text{kV}$. Nếu dây nối dài thêm 2m, điện áp tổng đặt lên cách điện máy biến áp sẽ tăng thêm 20kV, làm suy giảm nghiêm trọng biên hạn bảo vệ.

4. Các tiêu chuẩn quốc tế quy định biên hạn bảo vệ tối thiểu cho máy biến áp phân phối là bao nhiêu?

Theo tiêu chuẩn ANSI C62.2, biên hạn bảo vệ tối thiểu được thiết lập nhằm đảm bảo an toàn vận hành cách điện bao gồm: Biên hạn bảo vệ đầu sóng $\text{PM}_1 \ge 20%$, Biên hạn bảo vệ mức xung cơ bản $\text{PM}_2 \ge 20%$ và Biên hạn bảo vệ xung đóng cắt $\text{PM}_3 \ge 15%$. Các giá trị này tính toán dựa trên tỷ lệ giữa độ bền cách điện thiết bị và điện áp dư của chống sét van.

5. Làm thế nào để xác định các thông số $R_1$ và $L_1$ của mô hình chống sét van trong Matlab?

Thông số được tính toán khởi tạo dựa trên kích thước vật lý qua công thức $L_1 = 15d/n,(\mu\text{H})$ và $R_1 = 65d/n,(\Omega)$, với $d$ là chiều cao thiết bị tính bằng mét và $n$ là số cột song song. Sau đó, kỹ sư tiến hành hiệu chỉnh giá trị $L_1$ thông qua mô phỏng xung dòng chuẩn $8/20\mu s$ để điện áp dư mô hình khớp với bảng catalog của nhà sản xuất với sai số dưới 5%.

Kết luận

  • Xây dựng thành công mô hình toán và khối mô phỏng máy phát xung sét $8/20\mu s$ kết hợp chống sét van MOV phụ thuộc tần số trên nền tảng Matlab/Simulink.
  • Kiểm chứng độ chính xác cao của mô hình trên 5 thương hiệu thiết bị điện hàng đầu thế giới (Cooper, GE, Siemens, ELPRO, Ohio-Brass) với sai số điện áp dư dưới 5% ở dải dòng từ 1kA đến 20kA.
  • Lượng hóa thành công ảnh hưởng của độ dốc xung sét đến sự gia tăng 6% điện áp dư và chứng minh hiệu quả dập tắt quá áp khi lắp đặt chống sét van trên tuyến dây 15km.
  • Cung cấp bộ công cụ trực quan hỗ trợ đắc lực cho công tác đào tạo kỹ sư và tính toán phối hợp cách điện trên lưới điện phân phối trung áp 22kV - 35kV.
  • Định hướng phát triển giai đoạn 2027 - 2028: Mở rộng nghiên cứu đáp ứng động của chống sét van đối với xung sét dốc đứng $1/2\mu s$, phân tích trên lưới cáp ngầm đô thị và tích hợp thuật toán trí tuệ nhân tạo để chẩn đoán già hóa đĩa gốm MOV.

Hãy ứng dụng ngay mô hình mô phỏng chống sét van phụ thuộc tần số vào quy trình thiết kế và vận hành lưới điện để tối ưu hóa chi phí đầu tư, bảo vệ an toàn tuyệt đối cho các thiết bị điện tử công suất và nâng cao độ tin cậy cung cấp điện.