Tổng quan nghiên cứu

Trong hệ thống lưới điện phân phối, các sự cố do sét đánh gián tiếp chiếm tỷ lệ áp đảo với hơn 80% tổng số vụ quá điện áp lan truyền, đặc biệt là khi điểm phóng điện nằm cách đường dây trên 30 mét. Sét là hiện tượng phóng điện tự nhiên mang năng lượng khổng lồ với dòng điện đỉnh có thể đạt từ 30kA đến hơn 100kA, tạo ra trường điện từ biến thiên dữ dội trong không gian. Khi trường điện từ này cảm ứng lên các dây dẫn phân phối trung áp và 110kV ở độ cao từ 6 đến 15 mét, nó sinh ra xung quá điện áp cảm ứng tức thời lên đến hàng trăm kilovolt. Hệ quả thực tế là gây chọc thủng cách điện chuỗi sứ, nhảy máy cắt, kích nổ chống sét van và phá hỏng các thiết bị điều khiển điện - điện tử nhạy cảm, gây thiệt hại kinh tế hàng tỷ đồng mỗi năm.

Vấn đề cốt lõi của bài toán là làm sao tính toán chính xác biên độ và dạng sóng của điện áp cảm ứng mà không làm tiêu tốn quá nhiều tài nguyên máy tính. Các phương pháp mô phỏng số truyền thống như sai phân hữu hạn trong miền thời gian (FDTD) hay phương pháp kết hợp hàm bán kính cơ sở hai chiều (RBF-FDTD) trước đây đòi hỏi phải tối ưu hóa đồng thời cả 2 hệ số hình dạng cho không gian và thời gian, đồng thời bắt buộc phải có nghiệm giải tích đối chứng.

Để giải quyết triệt để hạn chế này, luận văn thạc sĩ chuyên ngành Kỹ thuật điện của học viên Nguyễn Hoàng Bảo, thực hiện tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh dưới sự hướng dẫn của Phó Giáo sư Vũ Phan Tú từ tháng 09/2023 đến tháng 01/2024, đã đề xuất phương pháp đường RBF (1D RBF-Line) cùng thuật toán tự động xác định hệ số hình dạng tối ưu c độc lập. Nghiên cứu tập trung khảo sát trên đường dây 1 pha, 3 pha và tuyến đường dây thực tế 110kV Đa Nhim - Hạ Sông Pha dài 1000 mét. Kết quả mang lại ý nghĩa đột phá với độ chính xác cao khi sai số RMS duy trì dưới 1% và thời gian tính toán giảm hơn 50% so với phương pháp RBF-FDTD 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 dựa trên sự liên kết chặt chẽ giữa 3 trụ cột lý thuyết kinh điển và hiện đại trong kỹ thuật điện từ:

Thứ nhất, mô hình kỹ thuật kênh sét đường dây truyền tải (Transmission Line Model - TL) do Uman và McLain phát triển, kết hợp hàm giải tích Heidler để mô tả dạng sóng dòng điện đáy kênh sét cho cả lần sét đánh đầu (First Stroke - FS với đỉnh 30kA) và các lần đánh kế tiếp (Subsequent Stroke - SS với đỉnh 12kA, thời gian đầu sóng 0.25 microgiây). Mô hình này giả định xung dòng điện di chuyển lên đám mây với vận tốc không đổi từ 1.0x10^8 m/s đến 1.5x10^8 m/s mà không bị suy giảm biên độ.

Thứ hai, lý thuyết trường điện từ kênh sét trên mặt đất dẫn có tổn hao. Luận văn ứng dụng công thức tích phân giải tích miền thời gian của Barbosa - Paulino (năm 2007) và công thức Cooray - Rubinstein (CR) miền thời gian mới phát triển để tính toán thành phần điện trường ngang, khắc phục hoàn toàn sự phức tạp của tích phân Sommerfeld kinh điển.

Thứ ba, mô hình ghép nối điện từ Agrawal, biểu diễn sự tương tác giữa trường điện từ phân tán và dây dẫn bằng hệ phương trình vi phân đạo hàm riêng kết hợp các nguồn áp dọc phân bố đại diện cho điện trường tiếp tuyến.

Các khái niệm then chốt trong nghiên cứu bao gồm: điện áp cảm ứng sét gián tiếp, hàm bán kính cơ sở (với 4 dạng: Gaussians - GA, Inverse Multiquadrics - IMQ, Inverse Quadrics - IQ, Multiquadrics - MQ), hệ số hình dạng tối ưu c, mặt đất dẫn tổn hao (lossy ground) và hiệu ứng phóng điện vầng quang (corona).

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

Nguồn dữ liệu đầu vào của đề tài sử dụng bộ thông số chuẩn hóa quốc tế với dòng điện sét SS biên độ 12kA, điện trở suất của đất biến thiên từ 0.01 S/m đến 0.001 S/m, hằng số điện môi tương đối bằng 10, cùng thông số hình học thực tế của tuyến đường dây 110kV mạch kép 173 Đa Nhim - 171 Hạ Sông Pha.

Cỡ mẫu tính toán được thiết lập thông qua việc rời rạc hóa không gian đường dây chiều dài 1000 mét thành 200 đến 500 đoạn vi phân với bước lưới không gian Δx từ 2 đến 5 mét, và bước thời gian mô phỏng Δt từ 5 đến 10 nano giây trong khoảng thời gian khảo sát tổng thể là 10 microgiây. Phương pháp chọn mẫu lưới chia được thực hiện đồng đều, tuân thủ nghiêm ngặt điều kiện ổn định Courant-Friedrichs-Lewy (CFL).

Phương pháp phân tích cốt lõi là phương pháp đường RBF (Mesh-free 1D), thay thế phép xấp xỉ đạo hàm không gian bậc 1 và bậc 2 bằng tổ hợp tuyến tính các hàm cơ sở bán kính trên sơ đồ vi phân 3 điểm cách đều. Điểm đột phá là thuật toán tìm hệ số hình dạng tối ưu cx dựa trên việc tối thiểu hóa ma trận thưa A(c) kết hợp khai triển chuỗi Taylor số dư cấp 4 và cấp 5 của điện áp và dòng điện. Toàn bộ thuật toán được lập trình mô phỏng trên nền tảng phần mềm kỹ thuật số từ ngày 04/09/2023 đến 18/12/2023.

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

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

Phát hiện thứ nhất về độ chính xác và tính tương thích của 4 dạng hàm RBF: Cả 4 hàm GA, IMQ, IQ và MQ đều cho kết quả điện áp cảm ứng bám sát tuyệt đối với phần mềm chuẩn quốc tế LiOV và nghiệm sai phân FDTD. Trong đó, hàm Multiquadrics (MQ) thể hiện hiệu năng xuất sắc nhất với sai số hiệu dụng RMS tại điểm giữa đường dây (x = 500m) chỉ dừng lại ở mức 0.35%, thấp hơn khoảng 40% so với hàm GA truyền thống.

Phát hiện thứ hai về hiệu quả tiết kiệm tài nguyên tính toán: Bằng cách chuyển đổi từ bài toán tối ưu 2 tham số (cx và ct) trong RBF-FDTD xuống chỉ còn 1 tham số không gian cx trong phương pháp đường RBF, thuật toán giúp giảm hơn 55% khối lượng bộ nhớ RAM tiêu thụ và rút ngắn thời gian xử lý ma trận trên máy tính hơn 50%, đồng thời loại bỏ 100% sự phụ thuộc vào nghiệm giải tích mẫu.

Phát hiện thứ ba về độ nhạy của các thông số vật lý kênh sét: Khi vận tốc lan truyền dòng điện kênh sét tăng từ 1.0x10^8 m/s lên 1.5x10^8 m/s, biên độ đỉnh của điện áp cảm ứng tại điểm x = 0m tăng vọt thêm 28.5%. Ngược lại, khi khoảng cách từ điểm sét đánh đến tâm đường dây tăng từ 100 mét lên 500 mét, giá trị điện áp đỉnh cảm ứng giảm mạnh hơn 65%.

Phát hiện thứ tư trên mô hình đường dây 3 pha và đường dây thực tế 110kV: Sự hiện diện của dây chống sét treo phía trên (Ground Wire - GW) đóng vai trò lá chắn cảm ứng vô cùng hiệu quả, giúp triệt tiêu và làm giảm biên độ điện áp cảm ứng trên các dây pha từ 25% đến 35% tùy thuộc vào cấu hình bố trí dây pha dạng ngang hay dọc.

Thảo luận kết quả

Khi xem xét tính chất đất dẫn có tổn hao thông qua công thức Barbosa - Paulino và Cooray - Rubinstein, thành phần điện trường ngang sinh ra có chiều ngược với điện trường trên đất dẫn lý tưởng. Điều này giải thích tại sao biên độ điện áp cảm ứng trên đất có điện trở suất cao bị suy giảm từ 15% đến 22% so với đất dẫn hoàn hảo, đồng thời thời gian tăng trưởng đầu sóng bị kéo dài ra. Kết quả mô phỏng cho thấy sự tương đồng trên 98% khi đối chiếu với các công bố khoa học của Barbosa (năm 2007) và Rachidi (năm 1993).

Về mặt trực quan hóa dữ liệu, các kết quả tính toán có thể được biểu diễn trực quan thông qua biểu đồ dạng sóng điện áp theo thời gian v(t) tại các tọa độ then chốt x = 0m (đầu đường dây), x = 250m, x = 500m (tâm đối xứng) và x = 1000m (cuối đường dây). Đi kèm với đó là bảng tổng hợp đối chuẩn sai số RMS và thời gian chạy mô phỏng giữa phương pháp FDTD cổ điển, RBF-FDTD và phương pháp đường RBF đề xuất, giúp làm nổi bật ưu thế vượt bậc của cấu trúc thuật toán mới.

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

Thứ nhất, tích hợp module tính toán phương pháp đường RBF vào hệ thống phần mềm quản lý lưới điện phân phối (DMS/SCADA): Đơn vị tư vấn thiết kế và các công ty điện lực cần triển khai tích hợp mã nguồn RBF-Line nhằm tự động hóa việc tính toán quá áp cảm ứng thời gian thực với mục tiêu xử lý dưới 2 giây cho một kịch bản sét đánh. Thời gian thực hiện dự kiến trong giai đoạn 2024 - 2025 do Trung tâm Nghiên cứu Phát triển PECC2 và các đơn vị thành viên EVN chủ trì.

Thứ hai, tối ưu hóa vị trí lắp đặt chống sét van (Surge Arrester - SA) trên các xuất tuyến 110kV và trung áp: Các kỹ sư vận hành cần ưu tiên bổ sung chống sét van không khe hở (ZnO) tại các cung đoạn đường dây đi qua khu vực đồi núi có điện trở suất đất lớn hơn 1000 Ω.m. Mục tiêu giảm ít nhất 30% số vụ sự cố phóng điện bề mặt chuỗi sứ do sét đánh gián tiếp trong vòng 12 tháng sau khi lắp đặt.

Thứ ba, chuẩn hóa khoảng cách treo dây chống sét và cấu hình hình học cột điện: Ban quản lý dự án lưới điện khi lập báo cáo kinh tế kỹ thuật các tuyến đường dây 110kV mới cần áp dụng mô hình ghép nối Agrawal kết hợp RBF để tính toán khoảng cách tối ưu giữa dây chống sét và dây pha. Target cần đạt là nâng cao hệ số che chắn điện từ thêm 20%, hoàn thành áp dụng trong quý IV/2024.

Thứ tư, ban hành quy trình tính toán hệ số hình dạng tối ưu c không cần nghiệm chuẩn: Nhóm nghiên cứu tại các trường đại học khối kỹ thuật cần phối hợp với Viện Năng lượng chuẩn hóa quy trình xác định hệ số cx thông qua ma trận thưa A(c) thành tiêu chuẩn kỹ thuật nội bộ, hướng tới kiểm soát sai số mô phỏng dưới 1% cho toàn bộ các đề án quy hoạch lưới điện truyền tải và phân phối đến năm 2026.

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

Nhóm thứ nhất là các kỹ sư thiết kế và chuyên gia phân tích hệ thống điện tại các tổng công ty điện lực và công ty tư vấn xây dựng điện: Tài liệu này cung cấp công cụ toán học trực tiếp giúp tính toán quá điện áp cảm ứng gián tiếp, từ đó lựa chọn đúng chủng loại cách điện, định mức chống sét van và tối ưu hóa chi phí đầu tư xây dựng trạm biến áp, đường dây.

Nhóm thứ hai là các nhà nghiên cứu, giảng viên và học viên sau đại học thuộc ngành Kỹ thuật điện, Năng lượng và Toán ứng dụng: Luận văn là tài liệu học thuật toàn diện về phương pháp không lưới (Mesh-free RBF), giải pháp ghép nối phương trình vi phân riêng phần Maxwell với đường dây truyền tải, và ứng dụng biến đổi trường điện từ trong miền thời gian.

Nhóm thứ ba là các kỹ sư phát triển phần mềm mô phỏng kỹ thuật (CAE/CAD): Đội ngũ lập trình có thể khai thác thuật toán giải ma trận thưa A(c) và phương pháp rời rạc hóa 1D RBF để nhúng vào các thư viện tính toán điện từ trường nhanh, thay thế cho các thuật toán FDTD truyền thống cồng kềnh.

Nhóm thứ tư là các cán bộ điều độ hệ thống điện và quản lý an toàn lưới: Sử dụng các số liệu đánh giá tác động của độ cao dây dẫn, vận tốc kênh sét và điện trở suất đất để xây dựng kịch bản phòng ngừa sự cố mất điện diện rộng vào mùa mưa bão có mật độ dông sét cao.

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

Sét đánh gián tiếp nguy hiểm như thế nào so với sét đánh trực tiếp trên lưới điện phân phối?
Sét đánh gián tiếp ở khoảng cách từ 30 đến 500 mét xảy ra với tần suất cao gấp 4 đến 5 lần so với sét đánh trực tiếp. Quá trình này cảm ứng ra các xung điện áp từ 50kV đến trên 200kV lên dây dẫn, đủ sức làm đánh thủng cách điện không khí và hư hỏng thiết bị trung - hạ áp trong thực tế.

Tại sao phương pháp đường RBF lại ưu việt hơn phương pháp RBF-FDTD truyền thống?
Phương pháp RBF-FDTD cũ đòi hỏi tối ưu 2 hệ số cx và ct và phải có sẵn nghiệm giải tích đối chứng. Phương pháp đường RBF mới chỉ tối ưu 1 hệ số cx qua ma trận thưa và chuỗi Taylor, giúp giảm hơn 50% thời gian tính toán và áp dụng được cho mọi tuyến dây thực tế không có nghiệm mẫu.

Điện trở suất của mặt đất ảnh hưởng ra sao đến biên độ điện áp cảm ứng?
Mặt đất có tổn hao làm suy giảm thành phần điện trường ngang do sinh ra phản xạ ngược pha. Đất có điện trở suất cao (như vùng đồi núi 1000 Ω.m) khiến biên độ đỉnh điện áp cảm ứng giảm khoảng 15% đến 22% nhưng làm biến dạng và kéo dài đuôi sóng quá áp.

Hàm bán kính cơ sở nào cho kết quả mô phỏng tối ưu nhất trong bài toán này?
Trong 4 hàm khảo sát gồm GA, IMQ, IQ và MQ, hàm Multiquadrics (MQ) thể hiện độ chính xác vượt trội nhất. Sai số RMS của hàm MQ chỉ đạt khoảng 0.35% tại điểm giữa đường dây, mang lại đường cong điện áp mượt mà và bám sát dạng sóng chuẩn của chương trình LiOV.

Luận văn đã kiểm chứng thuật toán trên tuyến đường dây thực tế nào?
Nghiên cứu đã áp dụng tính toán thành công trên tuyến đường dây 110kV mạch kép 173 Đa Nhim - 171 Hạ Sông Pha dài 1000 mét. Kết quả chứng minh dây chống sét giúp giảm từ 25% đến 35% biên độ quá áp cảm ứng trên 3 dây pha, khẳng định độ tin cậy tuyệt đối của phương pháp.

Kết luận

Luận văn thạc sĩ của tác giả Nguyễn Hoàng Bảo đã giải quyết xuất sắc bài toán tính toán quá điện áp cảm ứng do sét đánh gián tiếp thông qua các đóng góp trọng tâm:

  • Đề xuất thành công phương pháp đường RBF (1D RBF-Line) giúp giảm chiều tính toán và tiết kiệm hơn 50% thời gian mô phỏng số.
  • Phát triển thuật toán tự động xác định hệ số hình dạng tối ưu c độc lập dựa trên ma trận thưa A(c) và sai số Taylor cấp cao mà không cần nghiệm mẫu.
  • Ứng dụng chuẩn xác các công thức Barbosa - Paulino và Cooray - Rubinstein miền thời gian để mô phỏng điện trường ngang trên nền đất dẫn có tổn hao.
  • Đánh giá định lượng toàn diện ảnh hưởng của 4 yếu tố then chốt: vận tốc kênh sét, biên độ dòng điện đỉnh, độ cao treo dây và hiệu ứng vầng quang corona.
  • Kiểm chứng thực nghiệm thành công trên đường dây 110kV Đa Nhim - Hạ Sông Pha với sai số RMS cực thấp dưới 1%.

Trong giai đoạn 2024 - 2026, hướng phát triển tiếp theo của nghiên cứu sẽ tập trung mở rộng mô hình cho lưới điện đa nhánh phức tạp và tích hợp phần tử chống sét van phi tuyến. Các đơn vị tư vấn và vận hành lưới điện được khuyến nghị sớm ứng dụng phương pháp này để nâng cao độ tin cậy cung cấp điện quốc gia.