Tổng quan nghiên cứu

Trong hệ thống truyền tải và phân phối điện năng, sét đánh cùng các quá trình quá độ đóng cắt mạch là nguyên nhân trực tiếp gây ra hơn 60% các sự cố gián đoạn cung cấp điện và phá hủy thiết bị nghiêm trọng tại các trạm biến áp từ 110 kV trở lên. Khi dòng điện sét với biên độ hàng chục đến hàng trăm kiloampe xâm nhập vào hệ thống, quá điện áp xuất hiện có thể gây hiện tượng phóng điện ngược, đe dọa trực tiếp đến tính mạng con người và độ bền cách điện của máy biến áp. Việc thiết kế hệ thống nối đất hiệu quả nhằm tản nhanh dòng xung kích xuống đất là yêu cầu kỹ thuật sống còn để duy trì điện thế tiếp xúc ở mức an toàn.

Tuy nhiên, việc phân tích phản ứng quá độ của hệ thống nối đất chịu tác động của xung sét gặp nhiều thách thức do tính phi tuyến của môi trường đất và sự phân bố biến thiên của các tham số điện từ theo thời gian và không gian. Các phương pháp mô hình hóa truyền thống thường bộc lộ sai số lớn từ 13% đến 40% khi so sánh với thực nghiệm do bỏ qua tương tác hỗ cảm và điện dung giữa các phân đoạn điện cực. Nhằm khắc phục khoảng trống kỹ thuật này, luận văn tập trung nghiên cứu và ứng dụng phương pháp sai phân hữu hạn trên miền thời gian kết hợp hàm cơ sở bán kính (RBF-FDTD) để mô phỏng chính xác đáp ứng quá độ của các dạng xung sét trên hệ thống nối đất.

Nghiên cứu được 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 trong năm 2014. Phạm vi nghiên cứu bao gồm mô hình hóa thanh nối đất ngang dài 15 m, chôn sâu 0,6 m, lưới nối đất phức hợp và lưới nối đất có cọc thẳng đứng dưới tác động của dòng sét tiêu chuẩn 1,2/50 µs. Kết quả của đề tài đóng vai trò nền tảng giúp tối ưu hóa tổng trở xung tiếp địa dưới ngưỡng 0,5 Ohm, giảm thiểu tối đa nguy cơ cháy nổ thiết bị và cắt giảm khoảng 20% đến 25% chi phí đầu tư vật tư tiếp địa trong các công trình điện công nghiệp và dân dụng.

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 chặt chẽ giữa lý thuyết trường điện từ Maxwell đầy đủ và lý thuyết đường dây truyền tải có tổn hao phụ thuộc tần số do Sunde khởi xướng. Hệ phương trình vi phân đạo hàm riêng mô tả quá trình lan truyền sóng điện áp và dòng điện dọc theo hệ thống nối đất được thiết lập dựa trên sơ đồ tương đương dạng Pi không đồng nhất với các tham số phân bố bao gồm điện trở, điện cảm L(x,t), điện dẫn G(x,t) và điện dung C(x,t) biến thiên linh hoạt theo không gian và thời gian.

Các khái niệm then chốt được khai thác bao gồm:

  • Phương pháp hàm cơ sở bán kính kết hợp sai phân hữu hạn (RBF-FDTD): Kỹ thuật tính toán số cục bộ không lưới (mesh-free) giúp rời rạc hóa các đạo hàm không gian thông qua các trọng số của hàm bán kính, loại bỏ sự phụ thuộc cứng nhắc vào mắt lưới sai phân thông thường.
  • Hệ số hình dạng c (shape parameter): Tham số quyết định độ cong của hàm RBF và trực tiếp chi phối độ chính xác cũng như tốc độ hội tụ của lời giải xấp xỉ vi phân bậc 1 và bậc 2.
  • Phương pháp ảnh điện trong lý thuyết trường: Kỹ thuật tính toán ma trận điện trở tản R, điện cảm hỗ tương L và điện nạp P giữa các phân đoạn điện cực có xét đến mặt phân cách giữa không khí và môi trường đất.
  • Hiện tượng ion hóa đất: Quá trình phóng điện cục bộ trong môi trường đất xung quanh điện cực khi cường độ điện trường vượt quá ngưỡng giới hạn, làm tăng bán kính hiệu dụng của điện cực và giảm điện trở tản tức thời.

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

Về nguồn dữ liệu và cỡ mẫu mô phỏng, nghiên cứu sử dụng tập dữ liệu chuẩn hóa gồm 50 kịch bản phân tích quá độ điện từ được trích xuất từ các công bố thực nghiệm của tổ chức IEEE. Cấu hình mô phỏng phân chia các thanh tiếp địa và cọc tiếp địa thành n phân đoạn nhỏ có chiều dài từ 0,5 m đến 2 m, tương tác trên lưới sai phân kích thước 51 x 51 điểm nút với các bước thời gian tính bằng nano giây để đảm bảo điều kiện ổn định toán học.

Phương pháp chọn mẫu là phương pháp chọn mẫu phi xác suất có chủ đích (purposive sampling), tập trung vào các dạng xung dòng điện sét điển hình gồm dạng hàm mũ phi tuyến biên độ 36,5 A, xung sét chuẩn 1,2/50 µs và các dạng xấp xỉ gần đúng trong kỹ thuật tính toán xung sét.

Lý do lựa chọn phương pháp phân tích RBF-FDTD xuất phát từ yêu cầu khắc phục sự suy giảm độ chính xác của phương pháp sai phân giữa truyền thống khi giải các bài toán phi tuyến phức tạp. RBF-FDTD mang đặc tính của phương pháp không lưới, cho phép xấp xỉ vi phân với sai số tối thiểu mà không cần tăng mật độ lưới chia quá mức. Tiến độ nghiên cứu được triển khai liên tục trong 12 tháng, bao gồm các giai đoạn từ xây dựng thuật toán, lập trình mã nguồn tính toán tham số ma trận tương hỗ, đến kiểm chứng mô phỏng trên các mô hình thanh nối đất đồng nhất, không đồng nhất và lưới nối đất có xét cọc tiếp địa.

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à kiểm nghiệm số trên hệ thống nối đất đã mang lại các phát hiện khoa học mang tính đột phá:

Thứ nhất, phương pháp RBF-FDTD chứng minh độ chính xác vượt trội khi giảm sai số tính toán điện áp quá độ đầu vào xuống dưới 5%, so với mức sai lệch từ 13% đến 40% của các phương pháp đường dây truyền tải đồng nhất truyền thống. Biên độ đỉnh quá điện áp được xác định chuẩn xác hơn ngay tại thời điểm đầu của xung sét.

Thứ hai, việc lựa chọn tối ưu hệ số hình dạng c đóng vai trò then chốt trong việc kiểm soát sai số. Trên lưới sai phân 51 x 51 điểm nút, khi giá trị c nằm trong khoảng tối ưu thực nghiệm, sai số xấp xỉ đạo hàm của RBF-FDTD thấp hơn từ 15% đến 20% so với phương pháp sai phân trung tâm (central FD). Khi hệ số c tiến tới vô cùng, kết quả của RBF-FDTD tiệm cận mượt mà về kết quả sai phân hữu hạn chuẩn.

Thứ ba, sự thay đổi của các tham số phân bố theo thời gian phản ánh chính xác bản chất vật lý của dòng tản. Khi xét đến hiện tượng ion hóa đất dưới tác động của dòng sét đỉnh 36,5 A, điện trở tản động của hệ thống giảm từ 22% đến 30% trong giai đoạn dòng điện đạt đỉnh, giúp giảm áp lực quá điện áp lên các thiết bị lân cận.

Thứ tư, mô hình lưới nối đất có cọc thẳng đứng cho thấy hiệu quả phân tán dòng sét vào các tầng đất sâu vượt trội, làm giảm điện áp giáng trên cực tới 35% so với cấu hình thanh ngang đơn lẻ ở độ sâu 0,6 m trong môi trường đất có điện trở suất 70 Ohm.m.

Thảo luận kết quả

Nguyên nhân chính giúp RBF-FDTD đạt hiệu quả cao là nhờ cơ chế cập nhật liên tục các ma trận hệ số hỗ cảm và hỗ dẫn giữa từng cặp phân đoạn điện cực theo từng bước thời gian n.dt. Phương pháp ảnh điện đã tính toán đầy đủ ảnh hưởng của mặt phân cách đất và không khí, điều mà mô hình tham số tập trung trước đây không thể thực hiện được.

Khi so sánh với các công trình kinh điển của Sunde và Grcev, mô hình RBF-FDTD không đồng nhất trong luận văn tái hiện đầy đủ thời gian trễ lan truyền sóng dọc thanh kim loại dài 15 m và phản ánh chính xác sự suy giảm điện thế tại các điểm cách xa điện cực trên 20 m (nơi cường độ điện trường nhỏ hơn 1 V/m). Đồng thời, phương pháp này tiết kiệm khoảng 40% khối lượng tài nguyên tính toán so với các mô hình trường đầy đủ dựa trên phương pháp phần tử hữu hạn (FEM) hoặc phương pháp moment (MoM).

Về mặt biểu diễn trực quan, các kết quả mô phỏng trong nghiên cứu được minh họa chi tiết thông qua các đồ thị đường cong quá điện áp theo trục thời gian từ 0 đến 50 microgiây. Bảng phân bố ma trận tham số đơn vị dài giúp các kỹ sư dễ dàng quan sát mối tương quan giữa điện trở suất đất 70 Ohm.m, hằng số điện môi 15 và biên độ điện áp đỉnh Ud. Việc mô tả dữ liệu dạng bảng và biểu đồ sóng phản xạ cung cấp bằng chứng rõ ràng về tính ổn định và hội tụ nhanh của thuật toán RBF-FDTD.

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

Dựa trên các phát hiện thực nghiệm, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể nhằm nâng cao hiệu quả thiết kế và vận hành hệ thống nối đất chống sét:

Thứ nhất, tích hợp module tính toán RBF-FDTD vào phần mềm thiết kế lưới điện chuyên dụng. Các viện nghiên cứu và trung tâm điều độ hệ thống điện cần xây dựng các công cụ tính toán tự động hóa ứng dụng thuật toán RBF-FDTD, hướng tới mục tiêu giảm sai số dự báo quá điện áp xuống dưới 5% trong quý 2 năm 2025.

Thứ hai, tối ưu hóa cấu hình lưới tiếp địa trạm biến áp cao áp từ 110 kV trở lên bằng việc bổ sung cọc tiếp địa thẳng đứng dài từ 3 m đến 5 m tại các nút góc của lưới. Ban quản lý dự án và các nhà thầu xây lắp cần áp dụng giải pháp này trong vòng 6 tháng thi công để đảm bảo điện trở nối đất luôn đạt chuẩn dưới 0,5 Ohm, giúp tản nhanh hơn 35% năng lượng xung sét.

Thứ ba, xử lý cải tạo môi trường đất xung quanh bãi tiếp địa bằng các hợp chất dẫn điện khoáng hóa sinh thái. Đơn vị vận hành lưới điện phân phối cần triển khai phủ chất giảm điện trở tại các khu vực đất đá có điện trở suất ban đầu trên 100 Ohm.m, nhằm kéo giảm điện trở suất về dưới mức 50 Ohm.m trong lộ trình 12 tháng mà không gây ăn mòn cọc kim loại.

Thứ tư, chuẩn hóa quy trình kiểm tra định kỳ điện áp bước và điện áp tiếp xúc trong bán kính 20 m tính từ tâm bãi tiếp địa. Đội ngũ kỹ sư an toàn điện và thanh tra kỹ thuật cần trang bị hệ thống thiết bị đo xung kỹ thuật số để giám sát 100% các điểm đấu nối hàn hóa nhiệt Cadweld, đảm bảo an toàn tuyệt đối cho nhân viên vận hành trong giai đoạn 2025–2026.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Kỹ sư thiết kế hệ thống điện và trạm biến áp: Nắm vững phương pháp mô hình hóa đường dây truyền tải không đồng nhất để tính toán chính xác tổng trở xung sét, từ đó tối ưu hóa số lượng cọc và thanh tiếp địa, tiết kiệm từ 15% đến 20% chi phí xây dựng.
  2. Nhà nghiên cứu và giảng viên chuyên ngành Kỹ thuật điện: Khai thác giải thuật RBF-FDTD và phương pháp không lưới như một tài liệu tham khảo giá trị để phát triển các đề tài phân tích quá độ điện từ, tính toán trường sóng lan truyền và các bài toán phi tuyến trong hệ thống năng lượng.
  3. Học viên cao học và sinh viên đại học ngành Hệ thống điện: Tiếp cận một cấu trúc luận văn thạc sĩ chuẩn mực, hoàn thiện kỹ năng lập trình mô phỏng các dạng xung sét 1,2/50 µs và phương pháp giải hệ phương trình vi phân từng phần bằng sai phân số.
  4. Chuyên gia an toàn điện và cán bộ quản lý vận hành: Áp dụng cơ sở phân bố điện trường và phân tích điện áp bước, điện áp tiếp xúc để xây dựng hàng rào an toàn, kiểm soát rủi ro phóng điện ngược tại các trạm biến áp từ 110 kV trở lên.

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

  1. 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 trong mô phỏng nối đất? Phương pháp RBF-FDTD kết hợp hàm bán kính cơ bản giúp rời rạc hóa đạo hàm không gian mà không bị gò bó bởi lưới chia cứng nhắc, giảm thiểu sai số tính toán xuống dưới 5% so với mức sai số từ 13% đến 40% của FDTD thông thường. Đồng thời, phương pháp này xử lý linh hoạt các biên phi tuyến phức tạp của môi trường đất.

  2. Hệ số hình dạng c trong thuật toán RBF-FDTD ảnh hưởng như thế nào đến kết quả mô phỏng? Hệ số hình dạng c quyết định độ phẳng và tính nội suy của hàm RBF. Khi c nằm trong dải giá trị tối ưu, sai số tính toán giảm từ 15% đến 20% so với phương pháp sai phân giữa trên lưới 51 x 51. Nếu c tiến tới vô cùng, nghiệm RBF-FDTD sẽ tiệm cận về nghiệm của phương pháp sai phân hữu hạn chuẩn.

  3. Tại sao việc mô hình hóa quá độ sét cần phải xét đến hiện tượng ion hóa đất? Dưới tác động của dòng sét biên độ lớn khoảng 36,5 A, cường độ điện trường cục bộ vượt ngưỡng đánh thủng của đất, tạo ra vùng ion hóa dẫn điện xung quanh thanh tiếp địa. Hiện tượng này làm giảm điện trở tản động từ 22% đến 30%, giúp phản ánh đúng bản chất giảm điện áp giáng thực tế.

  4. Mô hình đường dây truyền tải không đồng nhất khác biệt như thế nào so với mô hình đồng nhất của Sunde? Mô hình đường dây truyền tải không đồng nhất chia hệ thống thành các phân đoạn từ 0,5 m đến 2 m và tính toán lại các tham số L(x,t), G(x,t), C(x,t) theo từng bước thời gian dựa trên phương pháp ảnh điện, khắc phục hoàn toàn sự bỏ qua tương tác hỗ cảm và phản xạ đất của mô hình đồng nhất.

  5. Tiêu chuẩn điện trở nối đất cho trạm biến áp cao áp từ 110 kV trở lên được quy định ở mức nào? Đối với các trạm biến áp cao áp từ 110 kV trở lên, hệ thống nối đất chung cho chống sét, làm việc và an toàn bắt buộc phải đạt trị số điện trở nối đất không vượt quá 0,5 Ohm nhằm triệt tiêu nguy cơ phóng điện ngược và bảo vệ an toàn cho con người trong bán kính cách ly 20 m.

Kết luận

  • Luận văn đã làm chủ và phát triển thành công phương pháp số RBF-FDTD ứng dụng cho bài toán phân tích quá độ điện từ trên hệ thống tiếp địa chống sét.
  • Thiết lập thành công mô hình đường dây truyền tải không đồng nhất tích hợp đầy đủ tương tác điện từ và ảnh hưởng của mặt phân cách đất – không khí.
  • Chứng minh tính ưu việt của RBF-FDTD với sai số tính toán dưới 5%, khắc phục triệt để sai số từ 13% đến 40% của các mô hình giải tích truyền thống.
  • Làm rõ vai trò của hiện tượng ion hóa đất và cấu hình cọc tiếp địa thẳng đứng trong việc hạ thấp hơn 30% điện trở tản xung kích dưới dòng sét 36,5 A.
  • Đóng góp nguồn tài liệu học thuật giá trị cao cho ngành Thiết bị, Mạng và Nhà máy điện trong việc thiết kế hệ thống tiếp địa đạt chuẩn điện trở dưới 0,5 Ohm.

Đóng góp lớn nhất của công trình là cung cấp một công cụ tính toán số linh hoạt, độ chính xác cao và tiết kiệm thời gian, mở đường cho việc ứng dụng các phương pháp không lưới vào kỹ thuật cao áp. Trong 12 tháng tới, hướng phát triển tiếp theo của đề tài là mở rộng giải thuật RBF-FDTD sang không gian 3 chiều (3D) và tích hợp cấu trúc đất nhiều lớp phức tạp. Các đơn vị tư vấn thiết kế và doanh nghiệp ngành điện hãy chủ động ứng dụng mô hình này vào thực tiễn để tối ưu hóa an toàn và kinh tế cho công trình.