Tổng quan nghiên cứu

Công nghệ truyền tải điện cao áp một chiều (HVDC) giữ vai trò then chốt trong việc giải quyết bài toán truyền tải điện năng công suất lớn trên các khoảng cách xa vượt trên 800 km. Trên thế giới, các dự án tiêu biểu như tuyến Itaipu - Sao Paulo vận hành ở cấp điện áp ±600 kV với tổng công suất truyền tải 6300 MW, hay tuyến Rihand - Delhi tại Ấn Độ công suất 1500 MW ở cấp điện áp ±500 kV trên chiều dài 814 km đã chứng minh tính ưu việt vượt trội của công nghệ một chiều so với lưới điện xoay chiều truyền thống. Mặc dù sở hữu những lợi thế như giảm thiểu tổn thất đường dây, không chịu ảnh hưởng bởi hiệu ứng bề mặt và không đòi hỏi bù công suất phản kháng, các hệ thống HVDC lại đối mặt với thách thức kỹ thuật lớn liên quan đến hiện tượng phóng điện vầng quang (corona discharge).

Khi điện áp vận hành vượt qua ngưỡng khởi tạo vầng quang, không khí xung quanh bề mặt dây dẫn bị ion hóa mạnh mẽ, tạo ra các đám mây điện tích không gian mang năng lượng cao. Dưới tác dụng của điện trường tĩnh điện, các ion này di chuyển liên tục về phía mặt đất và các công trình lân cận, tạo nên vùng điện trường ion hóa phức tạp tập trung trong phạm vi bán kính lên tới 3.5 lần chiều cao treo dây (3.5H). Sự tích tụ điện tích này không chỉ gây thất thoát công suất trên đường dây mà còn tiềm ẩn nguy cơ mất cân bằng điện tích tự nhiên, gây nhiễu sóng vô tuyến và tác động đến sức khỏe cộng đồng.

Nghiên cứu được hoàn thành vào tháng 06 năm 2013 tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh, tập trung giải quyết bài toán phân tích và tính toán phân bố điện trường ion hóa xung quanh các cấu trúc đường dây truyền tải HVDC. Luận văn đặt mục tiêu thiết lập quy trình tính toán số chuẩn xác dựa trên phương pháp phần tử hữu hạn (FEM) và thuật toán lưới thích nghi, cung cấp cơ sở dữ liệu định lượng phục vụ công tác thiết kế cách điện, tối ưu hóa hành lang an toàn và giảm thiểu rủi ro môi trường cho các dự án lưới điện tương lai.

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 vận dụng nền tảng vật lý phóng điện chất khí dựa trên lý thuyết ion hóa do va chạm của Townsend và công thức thực nghiệm của Peek để xác định điện trường khởi tạo vầng quang bề mặt dây dẫn ($E_0$). Cụ thể, đối với cấu trúc dây dẫn hình trụ bán kính $R$ đặt song song với mặt phẳng đất, điện trường tới hạn được xác định qua quan hệ thực nghiệm với giá trị cơ sở xấp xỉ 30 kV/cm đến 31.53 kV/cm tùy thuộc vào độ hoàn thiện bề mặt và bán kính dây.

Mô hình toán học mô tả không gian điện trường ion hóa bao gồm hệ phương trình vi phân đạo hàm riêng phi tuyến (PDEs), kết hợp chặt chẽ giữa phương trình Poisson biểu diễn mối liên hệ giữa điện thế ($\varphi$) và mật độ điện tích khối ($\rho$), phương trình liên tục của dòng điện tích và định luật Ohm mở rộng cho môi trường ion hóa:

  • Phương trình Poisson: $\nabla^2 \varphi = -\frac{\rho}{\varepsilon_0}$ (trong đó $\varepsilon_0 \approx 8.854 \times 10^{-12}\text{ F/m}$ là hằng số điện môi chân không).
  • Phương trình mật độ dòng ion: $\vec{j} = k \rho \vec{E}$ (với $k$ là độ linh động của ion, $\vec{E} = -\nabla \varphi$ là cường độ điện trường).
  • Phương trình liên tục ở trạng thái xác lập: $\nabla \cdot \vec{j} = 0$.

Để đơn giản hóa bài toán mà vẫn đảm bảo độ tin cậy kỹ thuật, nghiên cứu áp dụng các giả thuyết kinh điển bao gồm: giả thuyết Kaptzov (cường độ điện trường trên bề mặt dây dẫn duy trì không đổi ở giá trị khởi tạo $E_0$ khi xảy ra vầng quang) và giả thuyết Deutsch (sự hiện diện của điện tích không gian chỉ làm thay đổi độ lớn mà không làm biến dạng hướng của đường sức điện trường tĩnh). Các cấu trúc đường dây được phân loại chi tiết thành ba mô hình chính: đường dây đơn cực (Monopolar), đường dây lưỡng cực (Bipolar) và đường dây đồng cực (Homopolar).

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

Nghiên cứu ứng dụng phương pháp phần tử hữu hạn (Finite Element Method - FEM) dạng 2D kết hợp giải thuật tạo lưới tam giác Delaunay thích nghi. Phương pháp này chia miền khảo sát không gian liên tục thành hàng chục nghìn phần tử tam giác nhỏ, cho phép tùy biến mật độ lưới dày đặc tại các khu vực có gradient điện thế biến thiên mạnh gần bề mặt dây dẫn và giãn thưa dần về phía biên xa.

Quy trình tính toán số được xây dựng theo thuật toán lặp tương hỗ gồm 4 bước chính:

  • Bước 1: Khởi tạo điều kiện biên Dirichlet ($\varphi = V_{\text{app}}$ trên bề mặt dây dẫn, $\varphi = 0$ tại mặt đất) và điều kiện biên Neumann trên biên mở.
  • Bước 2: Thiết lập lưới tam giác Delaunay thích nghi trên miền không gian hai chiều có kích thước mở rộng từ 50 m đến 100 m quanh tim tuyến, với cỡ mẫu kiểm thử mô phỏng gồm hơn 15.000 phần tử hữu hạn bậc một và bậc hai.
  • Bước 3: Giải phương trình Poisson tìm phân bố điện thế ban đầu, sau đó cập nhật mật độ điện tích không gian dọc theo các đường sức thông qua phương pháp giải tích vi phân.
  • Bước 4: Kiểm tra điều kiện hội tụ của cường độ điện trường bề mặt với sai số tương đối giới hạn dưới mức $10^{-4}$; nếu chưa đạt, lặp lại việc hiệu chỉnh giá trị mật độ điện tích tại nút.

Toàn bộ thuật toán được cài đặt và kiểm chứng thông qua phần mềm mô phỏng chuyên dụng COMSOL Multiphysics kết hợp các đoạn mã xử lý dữ liệu ma trận. Quá trình thu thập thông số hình học, thiết lập mô hình và phân tích dữ liệu được thực hiện liên tục trong khoảng thời gian 12 tháng, từ ngày 02/07/2012 đến ngày 21/06/2013.

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

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

Phân tích số liệu mô phỏng thu được từ các mô hình đường dây truyền tải HVDC đem lại 4 phát hiện quan trọng:

  • Sự gia tăng đột biến của điện trường mặt đất do điện tích không gian: Tại khu vực mặt đất ngay dưới tim đường dây đơn cực mang điện áp danh định, sự xuất hiện của dòng thác ion khiến cường độ điện trường tổng hợp tăng từ 30% đến 45% so với trạng thái tĩnh điện thuần túy khi không có vầng quang.
  • Ưu thế triệt tiêu điện trường của cấu trúc lưỡng cực (Bipolar): Nhờ hiện tượng tái kết hợp giữa các ion dương (độ linh động $k^+ \approx 1.4 \times 10^{-4}\text{ m}^2/\text{V}\cdot\text{s}$) và ion âm ($k^- \approx 1.8 \times 10^{-4}\text{ m}^2/\text{V}\cdot\text{s}$) trong khoảng không gian giữa hai cực trái dấu, mật độ điện tích chạm đất của cấu trúc Bipolar giảm khoảng 35% so với cấu trúc Homopolar cùng cấp điện áp.
  • Quy luật dịch chuyển vị trí cực đại của cường độ điện trường: Cường độ điện trường lớn nhất tại mặt đất không nằm chính xác ở tâm đường dây mà phân bố lệch sang hai bên tim tuyến một khoảng cách từ 0.5H đến 1.2H (với H là chiều cao treo dây từ 15 m đến 25 m), sau đó suy giảm nhanh chóng khi khoảng cách vượt quá 3.5 lần chiều cao H.
  • Hiệu quả bảo vệ của hệ thống dây chắn điện trường: Khi bố trí bổ sung các dây chắn nối đất đặt bên dưới các dây dẫn mang điện ở khoảng cách hợp lý, cường độ điện trường và mật độ dòng ion cực đại tại bề mặt đất suy giảm đáng kể, ghi nhận mức giảm từ 25% đến 40% tùy thuộc vào vị trí và số lượng dây chắn.

Thảo luận kết quả

Cơ chế vật lý của sự gia tăng điện trường mặt đất bắt nguồn từ việc các hạt mang điện cùng dấu bị đẩy ra khỏi lớp ion hóa mỏng quanh dây dẫn và trôi tự do dưới tác dụng của lực điện trường. Quá trình này hình thành nên một lớp điện tích thể tích trong không gian, đóng vai trò như một nguồn sinh trường thứ cấp làm biến dạng phân bố điện thế ban đầu.

Khi so sánh với các nghiên cứu công bố quốc tế của các tác giả thuộc Viện Kỹ thuật Điện và Điện tử (IEEE) và các phòng thí nghiệm điện cao áp uy tín, kết quả giải thuật FEM thích nghi trong luận văn cho thấy sự tương đồng cao với độ sai lệch tương đối dưới 5%. Các kết quả nghiên cứu được trực quan hóa rõ nét thông qua hệ thống đồ thị đường cong phân bố điện trường $E(x)$ theo phương ngang tại các cao trình $y = 0\text{ m}$ và $y = 1.5\text{ m}$, bảng so sánh các giá trị cực đại $E_{\text{max}}$, $J_{\text{max}}$ giữa các cấu hình dây, cùng biểu đồ không gian 2D thể hiện các đường đẳng thế bao quanh dây dẫn.

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

Nhằm ứng dụng hiệu quả kết quả nghiên cứu vào thực tiễn quy hoạch, thiết kế và vận hành lưới điện truyền tải một chiều, 4 giải pháp trọng tâm được đề xuất:

  • Tối ưu hóa chiều cao treo dây và khoảng cách giữa các cực: Áp dụng phương pháp tính toán phần tử hữu hạn để xác định chiều cao treo dây tối thiểu từ 18 m đến 22 m cho cấp điện áp ±500 kV, nhằm kiểm soát cường độ điện trường tại mặt đất luôn duy trì dưới ngưỡng an toàn 25 kV/m và mật độ dòng ion dưới 100 nA/m² theo tiêu chuẩn quốc tế; giải pháp cần được Tập đoàn Điện lực Việt Nam (EVN) và các công ty tư vấn điện triển khai áp dụng trong giai đoạn thiết kế tiền khả thi 2026 - 2028.
  • Lắp đặt hệ thống dây chắn phụ trợ tại các khu vực nhạy cảm: Bố trí các cụm dây chắn nối đất bên dưới hành lang tuyến khi đường dây đi qua khu vực gần khu dân cư hoặc công trình hạ tầng, hướng tới mục tiêu giảm thiểu từ 25% đến 40% cường độ điện trường bề mặt; kế hoạch thực hiện kéo dài trong chu kỳ 12 - 18 tháng của giai đoạn thi công xây lắp.
  • Hoàn thiện khung quy chuẩn kỹ thuật quốc gia về hành lang an toàn HVDC: Ban hành quy định bắt buộc về khoảng cách ly an toàn tối thiểu tương đương 3.5 lần chiều cao treo dây (khoảng 50 m đến 70 m) đối với các dự án truyền tải một chiều; Bộ Công Thương và Cục Điều tiết Điện lực chủ trì soạn thảo trong lộ trình 2026 - 2027.
  • Phát triển phần mềm mô phỏng 3D tích hợp điều kiện khí tượng: Nâng cấp thuật toán mô phỏng FEM từ 2D lên 3D có tính đến vận tốc gió từ 2 m/s đến 10 m/s và độ ẩm không khí, nhằm nâng cao độ chính xác dự báo mật độ điện tích không gian lên trên 90%; các viện nghiên cứu chuyên ngành và trường đại học thực hiện trong thời gian 24 tháng.

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

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

  • Kỹ sư thiết kế và quản lý vận hành đường dây truyền tải: Nắm vững phương pháp tính toán phân bố trường ion hóa cho các cấp điện áp từ ±350 kV đến ±600 kV, hỗ trợ lựa chọn cấu hình dây, khoảng cách pha - đất và bố trí dây chắn tối ưu nhằm giảm tổn thất vầng quang.
  • Chuyên gia an toàn môi trường và thẩm định dự án năng lượng: Tiếp cận các công cụ và số liệu định lượng về cường độ điện trường và mật độ dòng ion trong phạm vi bán kính 50 m đến 100 m, phục vụ công tác lập báo cáo đánh giá tác động môi trường (ĐTM).
  • Giảng viên và nhà nghiên cứu kỹ thuật điện cao áp: Sử dụng tài liệu như một cẩm nang chuyên sâu về giải thuật phần tử hữu hạn thích nghi 2D và phương pháp giải hệ phương trình đạo hàm riêng phi tuyến trong kỹ thuật điện.
  • Học viên cao học và sinh viên ngành Hệ thống điện: Tham khảo quy trình xây dựng mô hình mô phỏng số trên COMSOL Multiphysics, tạo nền tảng vững chắc để thực hiện các đề tài nghiên cứu liên quan trong chu kỳ đào tạo 1 đến 2 năm.

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

  • Hiện tượng vầng quang trên đường dây HVDC khác biệt thế nào so với đường dây xoay chiều AC? Trên đường dây xoay chiều 50 Hz, sự đổi chiều liên tục của điện áp khiến các hạt mang điện bị giữ lại trong một vùng hẹp lân cận dây dẫn. Ngược lại, trên đường dây một chiều HVDC, các ion duy trì chuyển động liên tục về một hướng cố định dưới tác dụng của điện trường tĩnh, lấp đầy không gian giữa dây dẫn và mặt đất và làm tăng cường độ điện trường mặt đất lên 30% đến 45%.

  • Tại sao phương pháp phần tử hữu hạn (FEM) lại vượt trội trong giải bài toán trường ion hóa? Phương pháp phần tử hữu hạn cho phép rời rạc hóa miền không gian phức tạp thành hàng chục nghìn phần tử tam giác thích nghi, giúp xử lý chính xác các điều kiện biên phi tuyến tại bề mặt dây dẫn và mặt đất. So với các phương pháp giải tích truyền thống, sai số tính toán của FEM giảm xuống dưới mức 5%.

  • Điện tích không gian ảnh hưởng thế nào đến cường độ điện trường dưới mặt đất? Sự xuất hiện của đám mây điện tích không gian làm biến dạng mạnh mẽ đường sức tĩnh điện, khiến cường độ điện trường tại mặt đất gia tăng từ 30% đến 45% và đạt giá trị cực đại tại vị trí cách tâm đường dây một khoảng từ 0.5H đến 1.2H (với H là độ cao treo dây từ 15 m đến 25 m).

  • Dây chắn điện trường mang lại hiệu quả cụ thể ra sao trong thực tế vận hành? Dây chắn nối đất đặt bên dưới các dây pha đóng vai trò hấp thụ và phân tán bớt các đường sức điện trường cũng như dòng ion di chuyển xuống đất. Các tính toán số chứng minh giải pháp này giúp suy giảm trực tiếp từ 25% đến 40% giá trị cường độ điện trường cục bộ tại bề mặt đất.

  • Cấu trúc đường dây lưỡng cực (Bipolar) có ưu thế gì hơn so với cấu trúc đơn cực (Monopolar)? Cấu trúc Bipolar vận hành với hai dây dẫn mang cực tính trái dấu (+ và -), tạo điều kiện cho các ion dương và âm tái kết hợp tự nhiên trong không gian. Quá trình này giúp triệt tiêu một lượng lớn hạt mang điện tự do, làm giảm mật độ dòng ion chạm đất khoảng 35% so với cấu trúc chỉ mang một cực tính.

Kết luận

  • Xây dựng thành công thuật toán số hóa hoàn chỉnh dựa trên phương pháp phần tử hữu hạn (FEM) kết hợp lưới tam giác Delaunay thích nghi để giải quyết hệ phương trình Poisson - ion hóa phi tuyến cho các loại hình đường dây HVDC.
  • Định lượng chính xác mức độ tác động của điện tích không gian, khẳng định điện trường mặt đất tăng từ 30% đến 45% dưới ảnh hưởng của dòng thác ion vầng quang.
  • Xác định rõ quy luật phân bố không gian của trường ion hóa với vùng ảnh hưởng tập trung trong phạm vi 3.5 lần chiều cao treo dây và vị trí cực đại xuất hiện lệch tâm từ 0.5H đến 1.2H.
  • Chứng minh tính khả thi và hiệu quả vượt trội của giải pháp dây chắn điện trường trong việc cắt giảm từ 25% đến 40% cường độ điện trường chạm đất, bảo vệ an toàn cho các công trình dân sinh.
  • Đóng góp nguồn tài liệu kỹ thuật giá trị cao làm cơ sở mở rộng nghiên cứu mô hình 3D trong giai đoạn 12 đến 24 tháng tới, thúc đẩy tiến trình làm chủ công nghệ truyền tải điện siêu cao áp ±500 kV tại Việt Nam; các kỹ sư và chuyên gia có thể khai thác ngay bộ thuật toán này để tối ưu hóa thiết kế các dự án năng lượng mới.