Tổng quan nghiên cứu

Việt Nam nằm trọn trong khu vực nhiệt đới gió mùa ẩm với cường độ dông sét hoạt động vô cùng phức tạp và khắc nghiệt. Thống kê khí tượng ghi nhận số ngày dông cực đại tại một số khu vực lên đến 113,7 ngày tại Đồng Phú và số giờ dông đạt đỉnh 433,18 giờ tại Mộc Hóa. Trong thực tế vận hành hệ thống năng lượng, khoảng 50% tổng số sự cố lưới điện hàng năm bắt nguồn trực tiếp hoặc gián tiếp từ sét. Điển hình như sự cố ngày 04 tháng 06 năm 2001, sét đánh phá hủy máy cắt 220kV tại Nhà máy Thủy điện Hòa Bình đã làm tan rã mạch liên kết lưới điện miền Bắc, gây gián đoạn cung cấp điện nghiêm trọng.

Vấn đề cốt lõi đặt ra là sự bùng nổ của các thiết bị điện tử công suất nhạy cảm như bộ bù tĩnh SVC, hệ thống SCADA, bộ lưu điện UPS và các bộ lọc sóng hài trên lưới phân phối trung áp. Các thiết bị này sở hữu độ dự trữ cách điện rất thấp, dễ dàng bị đánh thủng khi xuất hiện xung quá điện áp lan truyền. Thiết bị chống sét van oxit kim loại không khe hở (MOV) đóng vai trò lá chắn then chốt bảo vệ phối hợp cách điện. Tuy nhiên, việc thử nghiệm thực nghiệm điện áp dư tại các phòng thí nghiệm chuyên ngành đòi hỏi chi phí đắt đỏ và khó triển khai diện rộng trên lưới thực tế.

Nghiên cứu tập trung vào mục tiêu xây dựng, chuẩn hóa và đánh giá độ chính xác của các mô hình toán học mô phỏng chống sét van trung áp dạng MOV phụ thuộc tần số trong môi trường Matlab/Simulink. Luận văn thực hiện kiểm chứng trên các dải điện áp định mức 18kV và 21kV với các biên độ dòng sét tiêu chuẩn 3kA, 5kA, 10kA và 20kA theo dạng sóng 8/20 µs. Công trình mang lại giá trị thực tiễn cao khi cung cấp công cụ mô phỏng chính xác với sai số điện áp dư duy trì dưới ngưỡng 10%, giúp các kỹ sư điện tối ưu hóa bài toán bảo vệ quá điện áp mà chỉ cần dựa vào các thông số cơ bản trong catalogue của nhà sản xuất.

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 lý thuyết vật lý chất bán dẫn oxit kim loại đa tinh thể gốm ZnO và lý thuyết phối hợp cách điện trong hệ thống điện (tiêu chuẩn Basic Impulse Level - BIL). Thành phần chính của đĩa biến trở MOV gồm hơn 90% ZnO kết hợp với các oxit kim loại phụ trợ như Bi2O3, Sb2O3, CoO, MnO và Cr2O3 qua quá trình thiêu kết nhiệt độ cao từ 1000°C đến 1400°C. Cấu trúc vi mô của MOV tạo nên hàng triệu tiếp giáp bán dẫn vi mô với điện thế rào khoảng 2,0V đến 3,5V tại biên hạt, mật độ hạt dẫn đạt khoảng 2x10^17 hạt/cm3 và điện trở suất vùng hạt đạt 0,1 đến 10 Ohm.cm trong khi điện trở suất vùng biên đạt 10^10 đến 10^12 Ohm.cm.

Đặc tính quan hệ điện áp - dòng điện (V-I) phi tuyến của MOV được xác định thông qua phương trình hàm mũ kinh điển biểu diễn mối quan hệ giữa dòng điện, điện áp và hệ số phi tuyến alpha. Để mô tả hành vi động phụ thuộc tần số khi có dòng sét xung kích dốc đứng đi qua, khung lý thuyết phân tích sâu 3 mô hình cấu trúc mạch tương đương:

  1. Mô hình tiêu chuẩn IEEE: Bao gồm hai nhánh điện trở phi tuyến A0 và A1 được ngăn cách bởi bộ lọc thông thấp R-L, phản ánh sự thay đổi trở kháng theo độ dốc sườn trước sóng sét.
  2. Mô hình Pinceti: Cải tiến đơn giản hóa từ mô hình IEEE bằng cách loại bỏ tụ điện ký sinh và gộp các điện trở chống mất ổn định số học thành một điện trở R giá trị 1 MegaOhm, tính toán thông số trực tiếp từ điện áp dư Vr8/20 và Vr1/T2.
  3. Mô hình P-K (Pinceti - Khalifa): Tối ưu hóa mạch tương đương bằng cách chỉ giữ lại độ tự cảm L1 giữa hai nhánh phi tuyến A0 và A1, giải quyết triệt để bài toán đáp ứng động với xung dòng có thời gian đầu sóng từ 0,5 µs đến 8 µs.

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

Luận văn kết hợp phương pháp nghiên cứu tài liệu chuyên sâu, phương pháp chuyên gia và phương pháp mô hình hóa mô phỏng số trên công cụ Matlab/Simulink. Nguồn dữ liệu thực nghiệm được trích xuất trực tiếp từ tài liệu kỹ thuật chính hãng của hai nhà sản xuất thiết bị điện hàng đầu thế giới: dòng chống sét van EVP của hãng Ohio Brass và dòng AZG2 của hãng Cooper.

Nghiên cứu áp dụng phương pháp chọn mẫu chủ đích với quy mô gồm 16 kịch bản thử nghiệm toàn diện: khảo sát 2 cấp điện áp trung áp thực tế (18kV và 21kV) tương ứng 2 chủng loại thiết bị, kích thích dưới 4 mức biên độ dòng xung sét tiêu chuẩn (3kA, 5kA, 10kA và 20kA) với dạng sóng chuẩn 8/20 µs và dạng sóng xung nhanh 1/T2 µs. Lý do lựa chọn phần mềm Matlab/Simulink cùng công cụ Mask Editor là vì khả năng xử lý các bài toán quá độ phi tuyến phức tạp, tính trực quan cao, không phụ thuộc vào các phần mềm đóng gói bản quyền đắt tiền và dễ dàng tích hợp vào các bài toán phân tích lưới điện phân phối tổng thể. Toàn bộ quá trình nghiên cứu, xây dựng thuật toán và chạy mô phỏng kiểm chứng được hoàn thành trong giai đoạn 6 tháng từ tháng 08 năm 2015 đến tháng 02 năm 2016.

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 trên 16 bộ số liệu kiểm thử đã mang lại những kết quả kỹ thuật nổi bật:

Thứ nhất, mô hình P-K chứng minh độ chính xác vượt trội nhất trong tất cả các mô hình được khảo sát. Đối với chống sét van hãng Cooper (cấp điện áp 18kV và 21kV), sai số biên độ điện áp dư của mô hình P-K dao động cực thấp, chỉ từ 0,04% đến 7,00%. Đối với thiết bị của hãng Ohio Brass, sai số điện áp dư của mô hình P-K duy trì ổn định ở mức xấp xỉ 2,20%. Toàn bộ các kết quả này đều nằm sâu dưới ngưỡng sai số cho phép của ngành điện là 10%.

Thứ hai, nghiên cứu xác nhận rõ rệt hiện tượng phụ thuộc tần số của chống sét van MOV. Khi thời gian đạt đỉnh của xung dòng sét giảm từ 8 µs xuống 1,3 µs (xung dốc đứng), biên độ điện áp dư sinh ra trên đầu cực van chống sét tự động tăng thêm xấp xỉ 6% đối với cùng một giá trị dòng điện đỉnh. Điều này chứng minh các mô hình điện trở phi tuyến tĩnh thông thường trong thư viện cơ bản của Matlab hoàn toàn không thể phản ánh chính xác trạng thái quá độ thực tế.

Thứ ba, việc ứng dụng các thuật toán giải tích trực tiếp để tính toán độ tự cảm L1 và L0 theo công thức của Pinceti và P-K đã loại bỏ hoàn toàn quá trình dò tìm thông số lặp phức tạp của mô hình IEEE truyền thống. Các giá trị điện cảm tính toán giúp mạch lọc phân chia dòng điện chính xác: ở sườn sóng dốc cao, phần lớn dòng điện đi qua nhánh điện trở A0 tạo điện áp kẹp tức thời, trong khi ở sườn sóng thoải, dòng điện phân bố đều qua nhánh A1.

Thảo luận kết quả

Cơ chế vật lý dẫn đến sự khác biệt về điện áp dư được giải thích bởi sự tương tác giữa điện cảm ký sinh bên trong khối gốm ZnO, cảm kháng của đầu dây nối và điện thế rào tại lớp biên hạt. Khi dòng sét có tốc độ biến thiên dòng điện (di/dt) cực lớn phóng qua, điện kháng của cuộn cảm L1 trong mạch tương đương tăng cao, ngăn cản dòng điện tức thời đi qua nhánh A1 và ép dòng sét tập trung chạy qua nhánh A0. Do đặc tuyến của A0 có mức điện áp cao hơn A1 ứng với cùng trị số dòng điện, điện áp dư tổng trên hai đầu MOV sẽ tự động vọt lố để triệt tiêu năng lượng xung sét.

Dữ liệu mô phỏng được hệ thống hóa trực quan thông qua các bảng so sánh sai số chi tiết và biểu đồ dạng sóng điện áp dư theo thời gian. Khi đối chiếu đồ thị xung điện áp dư giữa giá trị catalogue nhà sản xuất công bố và kết quả chạy trên Matlab, đường cong mô phỏng của mô hình P-K gần như trùng khít với đường đặc tính thực tế của cả hai hãng Ohio Brass và Cooper ở mọi mức xung 3kA, 5kA, 10kA và 20kA. So với mô hình ATP truyền thống vốn đòi hỏi các phương trình số mũ phức tạp dễ gây bất ổn định giải thuật số, mô hình khối chức năng Simulink tự xây dựng mang lại độ tin cậy và tính ổn định thuật toán vượt bậc.

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

Dựa trên các phát hiện thực nghiệm và kết quả mô phỏng, nghiên cứu đưa ra 4 khuyến nghị hành động cụ thể nhằm nâng cao hiệu quả bảo vệ quá điện áp trên hệ thống phân phối:

  1. Chuẩn hóa quy trình mô phỏng quá độ (Chủ thể: Tập đoàn Điện lực Việt Nam và các Tổng công ty Điện lực): Ban hành hướng dẫn kỹ thuật áp dụng mô hình chống sét van P-K trên nền tảng Matlab/Simulink để tính toán phối hợp cách điện cho các đường dây trung áp 22kV và 35kV. Mục tiêu giảm sai số tính toán quá áp xuống dưới 5% trong các dự án cải tạo lưới điện, hoàn thành triển khai trước quý 4 năm 2026.
  2. Thiết lập cơ sở dữ liệu số hóa thông số MOV (Chủ thể: Các Công ty Điện lực địa phương): Tiến hành thu thập và số hóa đầy đủ các thông số điện áp định mức Vr, điện áp dư Vr8/20 và Vr1/T2 của hơn 15.000 bộ chống sét van đang vận hành trên lưới. Thời gian hoàn thành trong vòng 12 tháng nhằm phục vụ công tác giám sát chất lượng thiết bị định kỳ.
  3. Nâng cấp cấu hình bảo vệ cho thiết bị nhạy cảm (Chủ thể: Các đơn vị tư vấn thiết kế điện): Bắt buộc tích hợp mô hình mô phỏng quá áp xung sét sườn dốc (thời gian đầu sóng từ 0,5 µs đến 2 µs) khi thiết kế trạm phân phối có lắp đặt biến tần, hệ thống SCADA và bộ bù SVC từ năm 2027. Giải pháp này giúp cắt giảm ít nhất 30% nguy cơ suy giảm cách điện và sự cố nổ thiết bị thứ cấp.
  4. Tích hợp bộ công cụ mô phỏng vào đào tạo chuyên sâu (Chủ thể: Các trường Đại học và Viện nghiên cứu khối kỹ thuật): Đưa các khối mô hình chống sét van phụ thuộc tần số (Masked Simulink Subsystems) vào chương trình thực hành môn học Kỹ thuật điện cao áp và Mô phỏng hệ thống điện ngay trong năm học 2026 - 2027.

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

Công trình nghiên cứu mang tính ứng dụng kỹ thuật cao, đặc biệt hữu ích cho 4 nhóm đối tượng trọng tâm sau:

  1. Kỹ sư vận hành và điều độ hệ thống điện: Giúp nắm vững cơ chế phản ứng động của thiết bị chống sét van trước các đợt dông sét thực tế, từ đó phân tích chính xác nguyên nhân sự cố chạm chập, phóng điện và xây dựng phương án bảo vệ tối ưu cho lưới điện phân phối trung áp.
  2. Kỹ sư tư vấn thiết kế trạm và đường dây: Cung cấp phương pháp luận chuẩn xác để lựa chọn trị số điện áp định mức và mức cách điện xung cơ bản (BIL) cho máy biến áp và thiết bị đóng cắt, đảm bảo tính kinh tế và an toàn kỹ thuật cao mà không cần đầu tư phần mềm chuyên dụng đắt tiền.
  3. Giảng viên và nhà nghiên cứu chuyên ngành Kỹ thuật Điện: Sử dụng tài liệu như một công trình tham khảo chuẩn mực về phương pháp mô hình hóa phần tử phi tuyến phụ thuộc tần số, phục vụ công tác giảng dạy chuyên đề cao áp và phát triển các đề tài nghiên cứu chuyên sâu.
  4. Học viên cao học và sinh viên ngành Điện - Điện tử: Cung cấp tài liệu học tập chi tiết, hướng dẫn từng bước từ khâu thiết lập phương trình toán, đóng gói giao diện Mask Editor trong Simulink đến kỹ năng phân tích dạng sóng quá độ trong hệ thống điện.

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

Mô hình P-K có ưu điểm gì vượt trội so với mô hình IEEE chuẩn?

Mô hình P-K loại bỏ các nhánh tụ điện phức tạp và chỉ giữ lại cuộn cảm L1 giữa hai phần tử phi tuyến A0 và A1. Nhờ đó, kỹ sư có thể xác định ngay các thông số thông qua công thức giải tích từ dữ liệu catalogue nhà sản xuất mà không cần thực hiện các thuật toán lặp phức tạp, giảm thiểu sai số tính toán xuống mức 0,04% đến 7,00%.

Cần những thông số nào trong datasheet để thiết lập mô hình chống sét van?

Người sử dụng chỉ cần trích xuất 3 thông số cơ bản có sẵn trong catalogue của mọi nhà sản xuất: Điện áp định mức của van chống sét (Vr), điện áp dư tương ứng dòng xung sét tiêu chuẩn 10kA dạng sóng 8/20 µs (Vr8/20), và điện áp dư tương ứng dòng sét xung nhanh 10kA dạng sóng 1/T2 µs (Vr1/T2).

Tại sao khối điện trở phi tuyến mặc định trong Matlab không mô phỏng chính xác sóng sét dốc đứng?

Khối điện trở phi tuyến mặc định của Matlab chỉ mô tả quan hệ V-I tĩnh theo hàm mũ cố định. Khi gặp xung sét có sườn trước tăng nhanh (dưới 1,3 µs), điện áp dư thực tế bị đội lên khoảng 6% do hiệu ứng động và điện kháng ký sinh, điều mà khối tĩnh không thể tái hiện được nếu thiếu cấu trúc lọc R-L phân nhánh.

Giới hạn sai số tiêu chuẩn cho phép trong mô phỏng chống sét van là bao nhiêu?

Theo các tiêu chuẩn kỹ thuật điện quốc tế và thực hành tính toán quá điện áp, sai số điện áp dư giữa kết quả mô phỏng so với giá trị đo đạc thực nghiệm của nhà sản xuất được chấp nhận khi nằm dưới ngưỡng 10%. Mô hình trong luận văn đạt độ chính xác cao với sai số thực tế chỉ từ 0,04% đến 7,00%.

Mô hình này có thể áp dụng cho các cấp điện áp khác ngoài 18kV và 21kV không?

Hoàn toàn có thể áp dụng rộng rãi. Thuật toán xác định thông số của mô hình P-K và Pinceti được chuẩn hóa theo hệ đơn vị tương đối và các biểu thức toán học tổng quát, cho phép kỹ sư nhập thông số datasheet của bất kỳ cấp điện áp nào từ hạ áp, trung áp (15kV, 22kV, 35kV) đến cao áp (110kV, 220kV).

Kết luận

  • Luận văn đã làm sáng tỏ bản chất vật lý và cơ chế hoạt động phi tuyến phức tạp của lớp biên hạt gốm ZnO trong thiết bị chống sét van MOV dưới tác động của dòng xung sét.
  • Xây dựng thành công bộ công cụ mô phỏng chống sét van phụ thuộc tần số hoàn chỉnh trên nền tảng phần mềm Matlab/Simulink với giao diện đóng gói Mask Editor thân thiện, trực quan.
  • Chứng minh mô hình P-K đạt độ chuẩn xác cao nhất với sai số điện áp dư chỉ từ 0,04% đến 7,00% đối với hãng Cooper và khoảng 2,20% đối với hãng Ohio Brass, vượt trội hơn mô hình IEEE truyền thống.
  • Đề xuất thuật toán xác định thông số mạch tương đương trực tiếp từ dữ liệu catalogue nhà sản xuất, giúp giải quyết triệt để rào cản phụ thuộc vào các phần mềm đóng gói đắt tiền.
  • Mở ra hướng nghiên cứu mở rộng mô hình cho các bài toán phân tích quá điện áp đóng cắt, quá điện áp duy trì và lan truyền sét trên toàn bộ sơ đồ lưới điện phân phối phức tạp.

Đóng góp lớn nhất của luận văn là đã lấp đầy khoảng trống nghiên cứu mô hình hóa chống sét van trung áp tại Việt Nam, mang lại một công cụ tính toán số tin cậy và kinh tế cho ngành điện. Trong giai đoạn 2026 - 2030, hướng phát triển tiếp theo sẽ tập trung tích hợp mô hình vào các bài toán thời gian thực (Real-time Digital Simulator) và đánh giá mức độ lão hóa của đĩa van ZnO sau nhiều lần phóng điện sét. Các đơn vị quản lý vận hành lưới điện và các viện nghiên cứu hãy chủ động ứng dụng ngay mô hình này vào công tác tính toán phối hợp cách điện thực tế để nâng cao độ tin cậy cung cấp điện và bảo vệ an toàn cho hệ thống năng lượng quốc gia.