Tổng quan nghiên cứu

Trong hệ thống năng lượng hiện đại, các sự cố điện như sét đánh, ngắn mạch hay rò rỉ điện áp chiếm tới hơn 70% nguyên nhân gây hư hại thiết bị và đe dọa an toàn tính mạng con người. Để kiểm soát rủi ro, hệ thống nối đất giữ vai trò cốt lõi trong việc tản nhanh dòng điện sự cố vào lòng đất. Tuy nhiên, việc tính toán phân bố điện thế truyền thống bằng các phương pháp số có lưới như sai phân hữu hạn (FDM) hay phần tử hữu hạn (FEM) thường gặp rào cản lớn về thời gian tạo lưới phức tạp và sai số suy biến hình học tại các vùng biên không gian phức tạp.

Luận văn thạc sĩ chuyên ngành Kỹ thuật điện của tác giả Trần Phước Chung, 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 PGS. Vũ Phan Tú và bảo vệ thành công ngày 22 tháng 08 năm 2020, tập trung giải quyết triệt để hạn chế trên. Mục tiêu trọng tâm của công trình là nghiên cứu, thiết lập và ứng dụng phương pháp không lưới hàm cơ sở bán kính (Meshfree-RBF) để mô phỏng chính xác sự phân bố điện thế trên hệ thống nối đất.

Phạm vi nghiên cứu bao gồm mô phỏng phân bố thế 2D và 3D từ các mô hình cơ bản như điện cực hình cầu bán kính 0.5m, cọc tiếp địa dài 5m trong không gian khảo sát 40m x 30m, cho đến lưới nối đất hình vuông 70m x 70m, lưới hình chữ L và hệ thống tiếp địa thực tế của Trạm biến áp 110kV khu công nghiệp Cổng Xanh, tỉnh Bình Dương. Ý nghĩa học thuật và thực tiễn của công trình thể hiện ở việc rút ngắn thời gian tạo lưới về mức 0 giây, đồng thời đạt độ chính xác vượt trội với sai số trung bình tiệm cận mức $4.95 \times 10^{-8}$, mở ra giải pháp mô phỏng số tối ưu cho kỹ thuật điện cao áp.

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 2 nền tảng lý thuyết cốt lõi: lý thuyết trường điện từ tĩnh trong môi trường dẫn điện đồng nhất và lý thuyết giải tích số không lưới dựa trên hàm cơ sở bán kính. Bản chất của bài toán phân bố điện thế khi có dòng điện ngắn mạch tản vào đất được mô tả bởi phương trình vi phân đạo hàm riêng Laplace trong không gian 3 chiều:

$$\nabla^2 u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} = 0$$

Luận văn phân loại hệ thống tiếp địa theo 3 chức năng độc lập: nối đất an toàn bảo vệ con người, nối đất làm việc ổn định điểm trung tính và nối đất chống sét tản dòng sét với quy chuẩn khoảng cách ly tối thiểu 5m. Để giải phương trình vi phân mà không cần rời rạc hóa miền tính toán thành các phần tử lưới, phương pháp Meshfree-RBF xấp xỉ trường thế $u(x)$ qua tổ hợp tuyến tính của các hàm bán kính đối xứng quanh các điểm tâm cố định. Ba kiểu hàm cơ sở bán kính tiêu biểu được áp dụng và so sánh chuyên sâu bao gồm:

  • Hàm Gaussian (GA): $\Phi(r) = e^{-(\varepsilon r)^2}$
  • Hàm Multiquadric (MQ): $\Phi(r) = \sqrt{1 + (\varepsilon r)^2}$
  • Hàm Inverse Multiquadric (IMQ): $\Phi(r) = \frac{1}{\sqrt{1 + (\varepsilon r)^2}}$

Trong đó, $r = |x - x_k|_2$ là khoảng cách Euclidean giữa điểm khảo sát và điểm nút, còn $\varepsilon$ là tham số hình dạng chi phối độ dốc và tính hội tụ của nghiệm xấp xỉ.

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

Nghiên cứu ứng dụng phương pháp mô phỏng số giải tích trên nền tảng phần mềm MATLAB. Nguồn dữ liệu kiểm chứng bao gồm bài toán chuẩn biên Dirichlet - Neumann với 289 điểm nút thực nghiệm (lưới $17 \times 17$) trên miền đơn vị $D = [0,1]^2$, kết hợp tập điểm đánh giá toàn cục có độ phân giải bước chia $1/40$.

Luận văn khảo sát 4 kỹ thuật bố trí điểm nút không gian: phân bố đồng nhất (Uniform), phân bố ngẫu nhiên (Random Uniform), phân bố Chebyshev và phân bố chuỗi giả ngẫu nhiên Halton. Việc lựa chọn phương pháp phân tích này xuất phát từ khả năng loại bỏ hoàn toàn ma trận liên kết phần tử vốn tiêu tốn tài nguyên tính toán của FEM, đồng thời giải quyết bài toán nội suy dữ liệu phân tán nhiều chiều một cách trực tiếp.

Nhằm xác định tham số hình dạng $\varepsilon$ tối ưu trong khoảng từ 1.0 đến 4.0 với bước quét mịn 0.01, nghiên cứu áp dụng thuật toán kiểm định chéo loại bỏ một phần tử (LOOCV). Toàn bộ quá trình nghiên cứu, xây dựng thuật toán và mô phỏng kiểm chứng được hoàn thành trong khung thời gian 6 tháng (từ tháng 02/2020 đến tháng 08/2020).

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 thực nghiệm và phân tích định lượng đã mang lại 4 phát hiện khoa học quan trọng:

  1. Hiệu năng vượt trội của kiểu phân bố Halton kết hợp hàm Gaussian: Trong bài toán thử nghiệm 2D, kiểu bố trí điểm nút Halton kết hợp hàm GA đạt sai số trung bình thấp nhất ở mức $4.955599 \times 10^{-8}$ và sai số cực đại là $1.14 \times 10^{-7}$. Kết quả này nâng cao độ chính xác lên hơn 85% so với kiểu phân bố Uniform (sai số trung bình $5.252432 \times 10^{-7}$). Ngược lại, hàm IMQ kết hợp phân bố Chebyshev cho sai số lớn nhất lên tới $6.244227 \times 10^{-6}$.
  2. Quy luật lan truyền thế của điện cực hình cầu: Với quả cầu bán kính 0.5m đặt trong không gian $40\text{m} \times 30\text{m}$, điện thế bề mặt đạt cực đại 1.0 pu tại tâm $x = 0$. Khi tăng độ sâu chôn điện cực từ 0m xuống 1.4m so với mặt đất, biên độ điện áp bề mặt giảm mạnh hơn 60%, giúp triệt tiêu đáng kể nguy cơ xuất hiện vùng gradient điện thế nguy hiểm.
  3. Giá trị tham số hình dạng $\varepsilon$ tối ưu cho cọc tiếp địa: Đối với cọc nối đất đơn dài 5m, việc khảo sát $\varepsilon$ từ 0.5 đến 9.0 cho thấy các giá trị $\varepsilon < 1.0$ hoặc $\varepsilon > 4.5$ làm điện thế biên lệch khỏi giá trị chuẩn 1.0 pu. Thuật toán LOOCV đã xác định chính xác điểm cực tiểu sai số tại $\varepsilon = 1.16$.
  4. Độ tương thích cao trên lưới tiếp địa phức tạp: Mô phỏng lưới nối đất vuông kích thước $70\text{m} \times 70\text{m}$ và lưới hình chữ L cho thấy kết quả của Meshfree-RBF hoàn toàn tương thích với phần mềm chuyên dụng ETAP, phương pháp FEM và FDM với độ lệch tương đối dưới 2.5%. Đặc biệt, khi bố trí thêm cọc tiếp địa tại các điểm biên của lưới 70m x 70m, biên độ điện thế tại các góc giảm rõ rệt, phân bố thế trên bề mặt đất trở nên bằng phẳng hơn.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp không lưới RBF đạt độ chuẩn xác cao là nhờ tính liên tục mượt mà cấp vô hạn của hàm cơ sở bán kính, cho phép nội suy trực tiếp trường thế mà không chịu ảnh hưởng bởi hiện tượng méo phần tử như trong FDM hay FEM. Dữ liệu mô phỏng được trực quan hóa sinh động thông qua các đồ thị phân bố đẳng thế 3D, biểu đồ mặt cắt điện áp 2D tại các vị trí tọa độ chiến lược ($y = 0\text{m}$, $y = 3.5\text{m}$, $y = 35\text{m}$) và các bảng so sánh sai số định lượng.

Tại Trạm biến áp 110kV KCN Cổng Xanh, mô phỏng đã bóc tách rõ ràng các vị trí có điện áp bước và điện áp tiếp xúc cao nhất tại các điểm rìa ngoài của trạm ($y = -27\text{m}$, $y = 9\text{m}$, $y = 21\text{m}$). Kết quả này khẳng định phương pháp không chỉ giải quyết tốt về mặt lý thuyết toán học mà còn phản ánh chính xác hiện tượng vật lý tản dòng trong các công trình công nghiệp thực tế.

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

Dựa trên kết quả đạt được, luận văn đưa ra 4 khuyến nghị kỹ thuật mang tính thực tiễn cao:

  • Chuẩn hóa quy trình tính toán Meshfree-RBF: Đề nghị Tập đoàn Điện lực Việt Nam (EVN) và các tổng công ty tư vấn thiết kế điện chuẩn hóa thuật toán Meshfree-RBF trong giai đoạn 2026 - 2027 để thay thế một phần các phép tính gần đúng truyền thống, hướng tới mục tiêu giảm 30% thời gian xử lý dữ liệu hình học trạm biến áp từ 110kV đến 500kV.
  • Tích hợp mô-đun tối ưu hóa tham số $\varepsilon$ tự động: Các nhóm nghiên cứu tại các trường đại học kỹ thuật cần phát triển công cụ phần mềm nội địa hóa, tích hợp thuật toán LOOCV tự động nhận diện tham số $\varepsilon \in [1.0, 1.5]$ cho từng cấu trúc điện cực, đảm bảo sai số mô phỏng luôn kiểm soát dưới $10^{-6}$ trong vòng 12 tháng tới.
  • Cải tạo lưới tiếp địa Trạm biến áp 110kV Cổng Xanh: Ban quản lý vận hành lưới điện Bình Dương cần triển khai bổ sung các cọc tiếp địa bổ trợ dài 3m - 5m tại các khu vực góc xung yếu ($y = -27\text{m}$ và $y = 21\text{m}$) nhằm giảm 15% đến 20% biên độ điện áp bước cục bộ trước quý 4 năm 2026.
  • Mở rộng mô hình tính toán cho môi trường đất đa lớp: Bộ Công Thương cùng các viện năng lượng cần tài trợ triển khai dự án nghiên cứu mở rộng thuật toán Meshfree-RBF cho môi trường đất không đồng nhất từ 2 đến 3 tầng điện trở suất, hoàn thiện bộ chỉ dẫn thiết kế tiếp địa quốc gia hoàn thành trong năm 2027.

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

Công trình luận văn là nguồn tài liệu học thuật và kỹ thuật giá trị cao cho 4 nhóm đối tượng:

  • Kỹ sư thiết kế trạm biến áp và đường dây cao thế: Nắm vững công cụ tính toán phân bố thế hiện đại, tối ưu hóa cấu hình lưới tiếp địa 110kV - 500kV, từ đó tiết giảm 10% - 15% chi phí vật liệu kim loại đồng và thép chôn ngầm mà vẫn đảm bảo tiêu chuẩn IEEE 80.
  • Giảng viên và nhà nghiên cứu chuyên ngành Kỹ thuật điện - Toán ứng dụng: Khai thác cơ sở lý thuyết toán học của phương pháp không lưới RBF, hàm cơ sở Gaussian, Multiquadric và giải thuật LOOCV để ứng dụng giải các phương trình đạo hàm riêng phi tuyến khác trong kỹ thuật điện từ.
  • Cán bộ quản lý an toàn và vận hành lưới điện: Nắm bắt quy luật phân bố điện thế tiếp xúc và điện thế bước thực tế tại các trạm biến áp công nghiệp, nâng cao hiệu quả đánh giá an toàn điện cho công nhân vận hành trong điều kiện sự cố ngắn mạch.
  • Học viên cao học và sinh viên ngành Hệ thống điện: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về phương pháp nghiên cứu khoa học, kỹ thuật mô phỏng trên nền tảng MATLAB với 7 chương nghiên cứu chi tiết và hệ thống dữ liệu kiểm chứng phong phú.

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

1. Phương pháp không lưới RBF có ưu thế gì vượt trội so với FEM và FDM?
Phương pháp Meshfree-RBF loại bỏ hoàn toàn giai đoạn phân chia lưới phần tử phức tạp, giúp thời gian tạo lưới giảm về 0 giây. Thuật toán chỉ sử dụng các điểm nút phân tán để nội suy, tránh triệt để tình trạng méo lưới khi mô phỏng cấu trúc phức tạp và đạt sai số cực nhỏ dưới $10^{-7}$.

2. Hàm cơ sở bán kính nào mang lại độ chính xác cao nhất cho bài toán thế điện trường?
Kết quả khảo sát khẳng định hàm Gaussian (GA) kết hợp phân bố điểm Halton cho độ chính xác cao nhất với sai số trung bình $4.955599 \times 10^{-8}$. Các hàm khác như Multiquadric (MQ) và Inverse Multiquadric (IMQ) vẫn đảm bảo tính hội tụ nhưng có sai số dao động trong khoảng từ $10^{-7}$ đến $10^{-6}$.

3. Vì sao việc lựa chọn tham số hình dạng $\varepsilon$ lại quyết định chất lượng mô phỏng?
Tham số $\varepsilon$ quyết định độ dốc và độ bao phủ của hàm cơ sở bán kính. Nếu $\varepsilon < 1.0$ hoặc $\varepsilon > 4.5$, kết quả tính toán sẽ sai lệch đáng kể so với điều kiện biên 1.0 pu ban đầu. Do đó, việc ứng dụng thuật toán LOOCV để tìm ra $\varepsilon = 1.16$ là bước bắt buộc để nghiệm hội tụ chính xác.

4. Luận văn đã ứng dụng phương pháp vào công trình thực tế cụ thể nào?
Luận văn đã tiến hành mô phỏng chi tiết hệ thống nối đất phức hợp của Trạm biến áp 110kV khu công nghiệp Cổng Xanh, tỉnh Bình Dương, xác định chính xác các điểm có điện áp bước và điện áp tiếp xúc cao nhất tại tọa độ $y = -27\text{m}$ và $y = 9\text{m}$.

5. Việc thay đổi độ sâu chôn điện cực ảnh hưởng thế nào đến an toàn điện trên mặt đất?
Khi tăng độ sâu chôn điện cực hình cầu bán kính 0.5m từ 0m xuống 1.4m dưới lòng đất, điện áp cực đại trên mặt đất giảm hơn 60% so với giá trị đỉnh 1.0 pu ban đầu. Điều này giúp giảm thiểu nguy hiểm từ điện áp tiếp xúc và điện áp bước cho người di chuyển trên mặt bằng trạm.

Kết luận

  • Luận văn đã xây dựng thành công bộ thuật toán không lưới Meshfree-RBF giải phương trình Laplace 3D cho bài toán phân bố thế trong hệ thống nối đất.
  • Chứng minh tính ưu việt của phân bố điểm Halton kết hợp hàm Gaussian với sai số trung bình đạt mức kỷ lục $4.955599 \times 10^{-8}$.
  • Tối ưu hóa thành công tham số hình dạng $\varepsilon = 1.16$ qua thuật toán LOOCV, giải quyết bài toán cọc tiếp địa 5m và lưới nối đất $70\text{m} \times 70\text{m}$.
  • Kiểm chứng thực tế thành công trên Trạm biến áp 110kV KCN Cổng Xanh (Bình Dương), cho kết quả tương thích với phần mềm ETAP và phương pháp FEM với độ lệch dưới 2.5%.
  • Đóng góp giải pháp mô phỏng số tiên tiến không cần chia lưới, tạo tiền đề vững chắc cho việc thiết kế an toàn hệ thống điện cao áp.

Lộ trình phát triển tiếp theo giai đoạn 2026 - 2028 sẽ tập trung mở rộng phương pháp cho môi trường đất đa tầng phức tạp và tích hợp mã nguồn thành phần mềm thương mại. Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm có thể khai thác chi tiết các thuật toán và bảng số liệu mô phỏng trong toàn văn công trình luận văn để ứng dụng vào công tác chuyên môn.