Tổng quan nghiên cứu

Quá điện áp cảm ứng do sét đánh lân cận là một trong những nguyên nhân hàng đầu gây ra sự cố phóng điện cách điện, ngắn mạch và gián đoạn cung cấp điện trên các hệ thống lưới điện truyền tải và phân phối. Theo các số liệu đo đạc thực nghiệm của Krider và Guo, năng lượng bức xạ điện từ sinh ra từ xung sét phóng điện đầu tiên có thể đạt tới mức 2x10^10 W, trong khi các xung phóng điện kế tiếp duy trì ở mức xấp xỉ 3x10^9 W. Khi kênh sét phóng điện xuống đất gần đường dây, trường điện từ bức xạ cực mạnh này sẽ cảm ứng trực tiếp lên các dây dẫn, tạo nên sóng quá điện áp lan truyền với biên độ có thể vượt quá mức chịu đựng cách điện cơ bản của đường dây phân phối từ 22 kV đến 110 kV.

Nghiên cứu tập trung giải quyết bài toán cốt lõi: nâng cao độ chính xác và hiệu quả tính toán biên độ cũng như dạng sóng điện áp cảm ứng sét trên đường dây cao áp và phân phối trên không. Luận văn đặt ra mục tiêu phát triển và hoàn thiện phương pháp sai phân hữu hạn trong miền thời gian kết hợp hàm bán kính cơ sở (RBF-FDTD), đồng thời tối ưu hóa thông số hình dạng c để vượt qua các hạn chế về sai số phân tán số của phương pháp sai phân truyền thống.

Phạm vi nghiên cứu được thực hiện toàn diện trong giai đoạn từ tháng 08 năm 2013 đến tháng 06 năm 2014 tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh. Các mô hình khảo sát bao gồm đường dây một pha và ba pha tiêu chuẩn đặt trong không gian đa chiều, cùng trường hợp kiểm chứng thực tế trên tuyến đường dây 173 Đa Nhim – 171 Hạ Sông Pha. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số tính toán quá điện áp từ 15% đến 25% so với phương pháp FDTD cổ điển, cung cấp công cụ tính toán tin cậy phục vụ công tác thiết kế cách điện và phối hợp bảo vệ chống sét cho hệ thống đ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 trên nền tảng hệ phương trình vi phân điện từ trường Maxwell kết hợp lý thuyết đường dây truyền tải đa vật dẫn của Agrawal để phân tích sự tương tác giữa trường điện từ ngoài và dây dẫn. Mô hình trường điện từ kênh sét được khảo sát toàn diện qua 4 mô hình dòng sét kinh điển: mô hình Đường dây truyền tải (TL), mô hình suy giảm tuyến tính (MTLL), mô hình suy giảm theo hàm mũ (MTLE) và mô hình Diendorfer – Uman (DU). Trường điện từ tổng hợp được tách thành hai thành phần chính: điện trường dọc sinh ra điện áp tới và điện trường ngang sinh ra điện áp nhiễu. Trường hợp mặt đất có tổn hao hữu hạn được mô hình hóa chính xác thông qua công thức Cooray – Rubinstein.

Khung lý thuyết tính toán số tích hợp 3 khái niệm nền tảng:

  • Phương pháp sai phân hữu hạn trong miền thời gian (FDTD) dựa trên lưới Yee để rời rạc hóa các phương trình sóng.
  • Hàm bán kính cơ sở (RBF) với ba dạng hàm chính: Multiquadrics (MQ), Inverse Multiquadrics (IMQ) và Gaussians (GA).
  • Hệ số hình dạng c (shape parameter), đại lượng quyết định độ phẳng và tính chính xác của phép xấp xỉ không gian trong phương pháp không lưới (Mesh-free).

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

Nguồn dữ liệu đầu vào sử dụng thông số dòng sét phóng ngược ghi nhận từ Trung tâm Không gian Kennedy của NASA với biên độ đỉnh dòng điện phóng thứ nhất đạt I01 = 9,9 kA và dòng điện phóng thứ hai đạt I02 = 7,5 kA, thời gian đầu sóng tăng trưởng cực nhanh trong khoảng 0,072 microgiây. Dữ liệu chuẩn đối sánh được trích xuất từ phần mềm chuyên dụng LIOV (Lightning Induced Overvoltage) và các công bố quốc tế uy tín trên IEEE của nhóm tác giả Rachidi và Nucci.

Phương pháp phân tích số sử dụng mô hình không gian 1 chiều và 2 chiều với cỡ mẫu tính toán gồm 500 đến 1000 phần tử nút không gian được chia đều với bước lưới Delta x = 5 m đến 10 m. Bước thời gian lấy mẫu Delta t được thiết lập nghiêm ngặt theo điều kiện ổn định Courant: Delta t <= Delta x / (c * sqrt(n)), với c = 3x10^8 m/s là vận tốc ánh sáng và n là số chiều không gian, đảm bảo bước thời gian luôn nhỏ hơn 1,67x10^-8 giây để triệt tiêu hoàn toàn sự mất ổn định số. Phương pháp RBF-FDTD được lựa chọn nhờ khả năng nội suy mượt mà các đạo hàm không gian bậc cao, giúp triệt tiêu hiện tượng méo dạng sóng quá điện áp khi lan truyền dọc theo chiều dài đường dây. Toàn bộ quy trình mô phỏng và thuật toán tối ưu hóa thông số c ngẫu nhiên được lập trình thực hiện xuyên suốt timeline 10 tháng nghiên cứu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Nghiên cứu đã chứng minh tính ưu việt vượt trội của thuật toán RBF-FDTD thông qua các kết quả định lượng cụ thể:

  • Độ chính xác xấp xỉ đạo hàm: Phương pháp RBF-FDTD kết hợp hàm Inverse Multiquadrics (IMQ) cho kết quả tiệm cận tốt nhất với phần mềm chuẩn LIOV, giúp giảm sai số biên độ đỉnh điện áp cảm ứng từ 8,4% của phương pháp FDTD truyền thống xuống còn dưới 2,1%, tương đương mức cải thiện độ chính xác hơn 75%.
  • Quy luật biến thiên của thông số c: Khi giá trị thông số hình dạng c tiến dần ra vô cùng, hệ số sai phân kx1 của RBF-FDTD hội tụ chính xác về giá trị sai phân trung tâm 1/(2*Delta x) của FDTD cổ điển. Tuy nhiên, tại khoảng giá trị c tối ưu được xác định bởi thuật toán ngẫu nhiên, sai số tính toán đạt cực tiểu.
  • Ảnh hưởng của độ dẫn điện mặt đất: Khi điện trở suất đất tăng từ trường hợp đất dẫn điện lý tưởng sang đất có tổn hao với độ dẫn điện 0,01 S/m đến 0,001 S/m, biên độ điện trường ngang tăng lên đáng kể, dẫn đến điện áp cảm ứng cực đại tại khoảng cách r = 100 m đến 500 m tăng từ 18% đến 32% và xuất hiện hiện tượng đảo cực tính ở đuôi sóng.
  • Hiệu quả bảo vệ của dây chống sét trên đường dây ba pha: Mô phỏng trên tuyến 173 Đa Nhim – 171 Hạ Sông Pha tại vị trí x = 500 m cho thấy, việc bố trí dây chống sét treo phía trên giúp cảm ứng ngược một điện áp trái dấu, làm giảm điện áp cảm ứng tổng trên các dây pha từ 42% đến 48% đối với cả cấu hình bố trí dây thẳng đứng và nằm ngang.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự vượt trội của phương pháp RBF-FDTD là do các hàm cơ sở như Multiquadrics và Inverse Multiquadrics có tính chất trơn cục bộ, cho phép xấp xỉ gradient điện trường biến thiên cực nhanh tại đầu sóng sét mà không đòi hỏi phải chia lưới không gian quá mịn như FDTD. Khi so sánh với các công bố của Rachidi trên tạp chí IEEE Transactions on Electromagnetic Compatibility, dạng sóng điện áp tính toán bởi RBF-FDTD hoàn toàn trùng khớp cả về thời gian đạt đỉnh (khoảng 1,5 đến 2,0 microgiây) và độ dốc sườn sóng.

Trong thực tế phân tích kỹ thuật, các dữ liệu này được thể hiện trực quan nhất thông qua đồ thị 2D biểu diễn dạng sóng điện áp theo thời gian tại các vị trí quan sát x = 0 m, 250 m, 500 m và 1000 m, kết hợp cùng bảng thống kê so sánh sai số tương đối giữa bốn biến thể thuật toán (MQ, IMQ, GA và FDTD) đối chiếu với LIOV. Kết quả khẳng định RBF-FDTD là công cụ toán học mạnh mẽ, giải quyết triệt để sự đánh đổi giữa khối lượng bộ nhớ tính toán và độ chính xác mô phỏng trong các bài toán quá độ điện từ trường.

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

  • Ứng dụng tích hợp thuật toán RBF-FDTD vào phần mềm tính toán quá độ lưới điện: Tập đoàn Điện lực Việt Nam (EVN) và các công ty tư vấn điện nên chủ trì triển khai tích hợp mô-đun RBF-FDTD vào các chương trình mô phỏng nội bộ trong vòng 6 đến 12 tháng tới, nhằm mục tiêu kiểm soát sai số tính toán quá điện áp dưới ngưỡng 3% cho toàn bộ lưới phân phối 22 kV và 35 kV.
  • Tối ưu hóa cấu hình treo dây chống sét và dây pha: Các Ban Quản lý dự án lưới điện cần áp dụng ngay việc lắp đặt bổ sung dây chống sét có nối đất định kỳ cho các khoảng cột đi qua khu vực có mật độ sét đánh cao (trên 4 lần/km2/năm). Mục tiêu giảm thiểu từ 40% đến 50% nguy cơ phóng điện sự cố, hoàn thành đánh giá lại thiết kế cho các tuyến xung yếu trong thời gian 12 tháng.
  • Khảo sát và cải tạo điện trở suất tiếp địa đất: Các đơn vị truyền tải và phân phối điện cần thực hiện đo đạc, lập bản đồ điện trở suất đất định kỳ 6 tháng một lần dọc theo các tuyến đường dây, tiến hành xử lý hóa chất giảm điện trở tiếp địa cột đạt trị số dưới 10 Ohm tại các vùng đất đồi núi đá có độ dẫn điện kém (dưới 0,001 S/m).
  • Tự động hóa quy trình xác định thông số hình dạng c: Các nhóm nghiên cứu tại các trường đại học kỹ thuật cần tiếp tục hoàn thiện thuật toán tối ưu hóa thông minh (như giải thuật di truyền GA hoặc tối ưu bầy đàn PSO) để tự động hóa hoàn toàn việc trích xuất hệ số c tối ưu trong vòng 3 đến 6 tháng, nâng cao tốc độ tính toán lên gấp 2 đến 3 lần.

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

  • Kỹ sư thiết kế và quy hoạch hệ thống điện: Nắm vững phương pháp tính toán điện áp cảm ứng chính xác để lựa chọn mức cách điện định mức, khoảng cách pha - đất và cấu hình dây chống sét tối ưu cho các dự án đường dây trên không từ 22 kV đến 220 kV.
  • Cán bộ quản lý vận hành và điều độ lưới điện: Ứng dụng các phân tích về sự suy giảm điện trường và quá điện áp cảm ứng để xây dựng quy trình bảo dưỡng tiếp địa, phân tích nguyên nhân sự cố rơle bảo vệ tác động do sét đánh ngoài đường dây trong bán kính 500 m.
  • Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật điện: Sử dụng luận văn như tài liệu học thuật chuyên sâu về sự kết hợp giữa phương pháp giải tích điện từ trường Maxwell và các phương pháp tính toán số hiện đại (Mesh-free RBF và FDTD).
  • Các nhà phát triển phần mềm mô phỏng kỹ thuật: Khai thác chi tiết các công thức toán học sai phân rời rạc và thuật toán tối ưu tham số c để xây dựng các công cụ phần mềm tính toán quá điện áp lan truyền mới thay thế cho các thuật toán sai phân truyền thống.

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

  • Phương pháp RBF-FDTD có ưu điểm gì vượt trội so với phương pháp FDTD truyền thống? Phương pháp RBF-FDTD kết hợp hàm bán kính cơ sở giúp nội suy không gian chính xác hơn, khắc phục hoàn toàn hiện tượng tán sắc số của lưới FDTD cổ điển. Kết quả thực nghiệm cho thấy phương pháp này giảm sai số tính toán biên độ đỉnh sóng quá điện áp từ 8,4% xuống dưới 2,1% so với phần mềm chuẩn LIOV mà không đòi hỏi phải chia lưới quá mịn làm tăng thời gian xử lý.

  • Thông số hình dạng c trong hàm RBF đóng vai trò như thế nào đối với độ chính xác của bài toán? Thông số hình dạng c quyết định độ dốc và phạm vi ảnh hưởng của hàm nội suy bán kính. Khi c tiến tới vô cùng, công thức RBF-FDTD sẽ tiệm cận về sai phân FDTD thông thường. Bằng cách sử dụng thuật toán tối ưu ngẫu nhiên để xác định giá trị c phù hợp cho từng loại hàm MQ, IMQ hay GA, sai số tính toán của bài toán điện từ trường được triệt tiêu về mức tối thiểu.

  • Tại sao điện trường ngang lại đóng vai trò quan trọng trong việc xác định điện áp cảm ứng sét? Mặc dù biên độ của điện trường ngang nhỏ hơn nhiều so với điện trường dọc, nhưng điện trường ngang xuất hiện do tổn hao mặt đất hữu hạn và tác động trực tiếp dọc theo chiều dài dây dẫn. Thành phần này quyết định trực tiếp đến điện áp nhiễu, làm thay đổi đáng kể dạng sóng, độ dốc đầu sóng và gây ra hiện tượng đổi dấu cực tính ở đuôi sóng điện áp cảm ứng.

  • Mô hình dòng sét nào được đánh giá là tối ưu nhất để tính toán điện áp cảm ứng gần đường dây? Mô hình đường dây truyền tải (TL) với giả thiết dòng sét lan truyền không suy giảm vận tốc (v khoảng 1,3x10^8 m/s) là mô hình tối ưu nhất để ước lượng biên độ đỉnh điện trường và điện áp cảm ứng ở khoảng cách gần dưới 5 km trong vài microgiây đầu tiên, nhờ công thức toán học đơn giản, tường minh và độ chính xác cao.

  • Dây chống sét giúp làm giảm điện áp cảm ứng trên dây pha như thế nào? Khi sóng điện từ kênh sét quét qua đường dây, dây chống sét được nối đất tại các vị trí cột sẽ sinh ra một điện áp cảm ứng ngược dấu với điện áp cảm ứng trên dây pha. Nhờ hiện tượng ghép nối điện từ hỗ cảm, thành phần điện áp này triệt tiêu một phần điện áp cảm ứng tổng, giúp giảm biên độ quá điện áp trên các dây pha từ 40% đến 50%.

Kết luận

  • Xây dựng thành công công thức sai phân trong miền thời gian kết hợp hàm bán kính cơ sở (RBF-FDTD) cho cả 3 dạng hàm Multiquadrics, Inverse Multiquadrics và Gaussians để giải phương trình điện từ trường kênh sét.
  • Chứng minh quy luật hội tụ toán học chặt chẽ: phương pháp FDTD truyền thống chính là một trường hợp đặc biệt của RBF-FDTD khi thông số hình dạng c tiến đến vô cùng.
  • Phát triển thuật toán ngẫu nhiên tối ưu hóa thông số hình dạng c, giúp nâng cao độ chính xác mô phỏng điện áp cảm ứng trên đường dây một pha và ba pha vượt trội so với FDTD khi đối chiếu với chuẩn LIOV.
  • Khảo sát toàn diện ảnh hưởng của các yếu tố vật lý thực tế: vận tốc dòng sét (1,3x10^8 m/s), biên độ dòng phóng (9,9 kA), điện trở suất đất có tổn hao và hiệu ứng vầng quang đến dạng sóng quá điện áp lan truyền.
  • Ứng dụng thành công mô hình tính toán trên tuyến đường dây thực tế 173 Đa Nhim – 171 Hạ Sông Pha tại các mốc chiều dài 0 m, 500 m và 1000 m, làm cơ sở khoa học tin cậy cho thiết kế bảo vệ chống sét lưới điện.

Kế hoạch phát triển tiếp theo trong vòng 12 tháng tới sẽ mở rộng thuật toán RBF-FDTD cho hệ thống lưới điện phân phối phức tạp có cấu trúc phân nhánh và xét đến ảnh hưởng phi tuyến của điện trở nối đất cột khi có dòng sét lớn đi qua. Quý độc giả và các kỹ sư quan tâm có thể ứng dụng trực tiếp khung thuật toán này để nâng cao độ an toàn tin cậy cho công trình đường dây truyền tải điện.