Tổng quan nghiên cứu

Sự phát triển mạnh mẽ của nền kinh tế thúc đẩy nhu cầu tiêu thụ điện năng tăng trưởng trung bình từ 10% đến 12% mỗi năm tại Việt Nam, kéo theo việc mở rộng mạng lưới truyền tải điện siêu cao áp 500kV đi sâu vào các trung tâm phụ tải và vùng dân cư đông đúc. Tuy nhiên, sự xuất hiện của các đường dây truyền tải điện áp cao phát sinh trường điện từ tần số cực thấp 50Hz trong không gian xung quanh, làm dấy lên những lo ngại sâu sắc về an toàn sức khỏe con người, hiện tượng nhiễu loạn thông tin viễn thông và suy hao vầng quang điện.

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học và ứng dụng phương pháp phần tử hữu hạn kết hợp thuật toán lưới thích nghi để mô phỏng, tính toán chính xác sự phân bố cường độ điện trường và cảm ứng từ quanh đường dây truyền tải điện. Phạm vi nghiên cứu tập trung khảo sát tuyến đường dây siêu cao áp 500kV Cai Lậy - Ô Môn tại khu vực Đồng bằng Sông Cửu Long, phân tích không gian từ mặt đất lên đến độ cao 13 mét.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc cung cấp công cụ tính toán số tin cậy trên nền tảng phần mềm kỹ thuật, giúp giảm thiểu sai số thiết kế cơ học truyền thống xuống dưới 3%. Đồng thời, kết quả mô phỏng hỗ trợ đắc lực cho các cơ quan quản lý năng lượng kiểm soát nghiêm ngặt ngưỡng cường độ điện trường khu vực công cộng dưới 5 kV/m và cảm ứng từ dưới 100 microTesla theo quy định an toàn hiện hành.

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 nền tảng hệ phương trình đạo hàm riêng Maxwell nhằm mô tả mối quan hệ tương hỗ giữa bốn vector đặc trưng của trường điện từ gồm cường độ điện trường, cảm ứng từ, mật độ dòng điện dịch và cường độ từ trường. Từ hai định luật cơ bản của Faraday và Ampere, phương trình truyền sóng Helmholtz cho môi trường không đồng nhất được thiết lập dưới dạng phức điều hòa ở tần số 50Hz, cho phép tách biệt độc lập bài toán phân bố điện trường và từ trường trong không gian hai chiều.

Khung lý thuyết vận dụng ba khái niệm trụ cột: trường tĩnh điện trạng thái xác lập trong môi trường điện môi không đồng nhất, bức xạ không ion hóa tần số cực thấp, và nguyên lý bảo toàn năng lượng trường tĩnh điện thông qua phiếm hàm năng lượng cực tiểu. Nghiên cứu xác lập các điều kiện biên Dirichlet và Neumann chặt chẽ tại bề mặt tiếp giáp giữa các môi trường chất khác nhau như dây dẫn kim loại, không khí và mặt đất phẳng, bảo đảm tính liên tục của thành phần tiếp tuyến và bước nhảy của thành phần pháp tuyến.

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

Nguồn dữ liệu tính toán sử dụng thông số thực tế của tuyến đường dây 500kV Cai Lậy - Ô Môn với điện áp danh định 500kV, cấu hình dây dẫn phân pha, độ võng dây và khoảng cách pha theo hai dạng hình học: mạch đơn bố trí hình tam giác và mạch kép bố trí đối xứng thẳng đứng.

Nghiên cứu ứng dụng phương pháp phần tử hữu hạn kết hợp kỹ thuật phần dư có trọng số Galerkin để giải hệ phương trình vi phân đạo hàm riêng. Miền khảo sát không gian hai chiều có kích thước bề rộng 70 mét và chiều cao 35 mét được rời rạc hóa thành các phần tử tam giác ba nút bậc nhất. Lý do lựa chọn phương pháp phần tử hữu hạn là nhờ khả năng thích ứng linh hoạt với các miền hình học phức tạp, cho phép gán thuộc tính vật liệu điện môi cục bộ và đạt độ hội tụ chính xác cao hơn hẳn phương pháp sai phân hay giải tích cổ điển.

Để tối ưu hóa tài nguyên tính toán, giải thuật tạo lưới Delaunay thích nghi được phát triển dựa trên hàm điều khiển lực cân bằng nút. Trong mô hình thử nghiệm, lưới thích nghi chỉ cần 758 nút và 1356 phần tử, giảm hơn 52% số lượng nút tính toán so với lưới Delaunay đều cổ điển (yêu cầu 1594 nút và 2929 phần tử). Toàn bộ thuật toán được lập trình và giải ma trận đại số tuyến tính trên môi trường MATLAB trong thời gian thực hiện 5 tháng.

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

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

Thứ nhất, cường độ điện trường đạt giá trị cực đại ngay tại vị trí tim tuyến và suy giảm nhanh chóng khi tiến ra mép hành lang an toàn. Tại độ cao 1 mét so với mặt đất (vị trí tiếp xúc cơ thể người), cường độ điện trường cực đại của đường dây mạch kép 500kV đạt khoảng 4,2 kV/m, thấp hơn khoảng 18% so với đường dây mạch đơn (đạt xấp xỉ 5,1 kV/m), bảo đảm nằm trong ngưỡng an toàn cho phép 5 kV/m đối với khu vực công cộng.

Thứ hai, cường độ từ trường phụ thuộc tuyến tính vào dòng điện phụ tải; ở chế độ vận hành định mức, mức cảm ứng từ tại độ cao 1 mét dao động từ 15 đến 22 microTesla, thấp hơn rất nhiều so với ngưỡng giới hạn phơi nhiễm công cộng 100 microTesla do Ủy ban Quốc tế về Bảo vệ Bức xạ Không Ion hóa ban hành.

Thứ ba, giải thuật lưới Delaunay thích nghi giúp tối ưu hóa thời gian xử lý ma trận trên phần mềm MATLAB xuống còn khoảng 0,68 giây, rút ngắn hơn 53% thời gian so với mức 1,45 giây của lưới phân bố đều thông thường. Sai số tính toán giữa phương pháp phần tử hữu hạn và phương pháp mật độ điện tích truyền thống chỉ dao động trong khoảng 1,2% đến 2,5%.

Thứ tư, khi tăng chiều cao treo dây dẫn từ mốc 10 mét lên 13 mét so với mặt đất, cường độ điện trường bề mặt tiếp xúc giảm trung bình khoảng 32%, tạo ra vùng đệm an toàn rộng hơn cho các hoạt động dân sinh bên dưới đường dây.

Thảo luận kết quả

Sự suy giảm điện trường và từ trường rõ rệt ở cấu hình mạch kép được giải thích bởi hiện tượng triệt tiêu tương hỗ giữa các vector trường của các pha đối xứng khi được hoán vị thứ tự pha hợp lý. Điều này khẳng định việc sắp xếp không gian của các dây dẫn mang dòng có ý nghĩa quyết định đối với mức độ bức xạ trường điện từ.

Dữ liệu mô phỏng được biểu diễn trực quan thông qua đồ thị phân bố không gian hai chiều theo phương ngang từ -35 mét đến +35 mét, kết hợp bản đồ đường đồng mức điện thế và biểu đồ ma trận màu ba chiều thể hiện mật độ năng lượng trường. Bảng đối chiếu số liệu tại 7 mức độ cao khác nhau gồm 1m, 2m, 3m, 4m, 5m, 10m và 13m cung cấp cái nhìn chi tiết về sự biến thiên của điện từ trường theo phương thẳng đứng.

So với các công bố quốc tế của Viện Kỹ sư Điện và Điện tử cũng như Hội đồng Lưới điện Lớn Quốc tế, kết quả đạt được có sự tương đồng cao về mặt quy luật phân bố. Phương pháp phần tử hữu hạn chứng minh tính ưu việt vượt trội trong việc xử lý các bài toán biên không gian mở và môi trường điện môi nhiều lớp so với các phương pháp giải tích giản lược.

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

Thứ nhất, tối ưu hóa thứ tự hoán vị pha cho tất cả các tuyến truyền tải mạch kép 500kV chuẩn bị xây dựng mới nhằm triệt tiêu tối đa trường điện từ tổng hợp, hướng tới mục tiêu giảm từ 15% đến 25% cường độ điện trường tại tâm hành lang tuyến, do các Ban Quản lý dự án điện triển khai ngay trong giai đoạn lập báo cáo nghiên cứu khả thi trong vòng 6 tháng.

Thứ hai, nâng cao độ tĩnh không dây dẫn tại các điểm giao cắt khu dân cư tập trung từ mức tối thiểu 10 mét lên trên 14 mét, giúp giảm hơn 30% cường độ điện trường tiếp xúc mặt đất, do các đơn vị tư vấn thiết kế cập nhật vào hồ sơ kỹ thuật trong khung thời gian từ 6 đến 12 tháng.

Thứ ba, lắp đặt bổ sung hệ thống dây chống sét kết hợp dây che chắn tiếp địa hoặc trồng thảm cây xanh tán thấp có độ ẩm cao dọc hành lang bảo vệ nhằm phân tán và hấp thụ điện trường, bảo đảm duy trì cường độ điện trường tại khu vực sinh hoạt dưới mức 3,5 kV/m, do Tổng công ty Truyền tải điện Quốc gia phối hợp chính quyền địa phương thực hiện định kỳ hàng năm.

Thứ tư, tích hợp công cụ mô phỏng số phần tử hữu hạn trên nền tảng MATLAB vào quy trình tự động hóa đánh giá tác động môi trường của Tập đoàn Điện lực Việt Nam, hoàn thành việc xây dựng quy chuẩn và ban hành sổ tay hướng dẫn tính toán nội bộ trong thời hạn 18 tháng.

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

Kỹ sư thiết kế đường dây và trạm biến áp: Tiếp cận công cụ tính toán số chính xác để tối ưu hóa hình học tháp điện, khoảng cách pha và chiều cao treo dây, giúp tiết kiệm từ 5% đến 10% chi phí vật liệu xây lắp trong khi vẫn bảo đảm an toàn bức xạ điện từ.

Cán bộ quản lý vận hành và điều độ lưới điện: Sử dụng hệ thống bảng tra số liệu cường độ điện trường và cảm ứng từ tại các cấp điện áp 110kV, 220kV, 500kV để xây dựng quy trình bảo hộ lao động và khoảng cách an toàn khi công nhân làm việc gần đường dây đang mang điện.

Chuyên gia đánh giá tác động môi trường và nhà hoạch định chính sách: Khai thác các dữ liệu đối chiếu với tiêu chuẩn quốc tế như Tổ chức Y tế Thế giới và Ủy ban Quốc tế về Bảo vệ Bức xạ Không Ion hóa để xây dựng báo cáo an toàn môi trường và hoàn thiện quy chuẩn kỹ thuật quốc gia về hành lang bảo vệ lưới điện cao áp.

Giảng viên, học viên cao học và sinh viên ngành Kỹ thuật điện: Sử dụng mô hình giải thuật lưới Delaunay thích nghi và mã nguồn giải phương trình đạo hàm riêng trên MATLAB làm tài liệu giảng dạy, học tập và phát triển các hướng nghiên cứu chuyên sâu về trường nhiệt cáp ngầm hoặc trường điện từ ba chiều.

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

Phương pháp phần tử hữu hạn có ưu thế gì vượt trội so với các phương pháp giải tích truyền thống khi tính toán điện từ trường? Phương pháp phần tử hữu hạn cho phép rời rạc hóa các miền không gian hình học phức tạp và giải quyết chính xác các điều kiện biên không đồng nhất của hệ thống dây dẫn thực tế. Điều này giúp loại bỏ các giả định đơn giản hóa thái quá, giảm sai số tính toán xuống dưới 3% so với mức sai số trên 10% của các công thức giải tích cổ điển.

Thuật toán lưới Delaunay thích nghi mang lại lợi ích cụ thể nào cho bài toán mô phỏng? Thuật toán tự động tái bố trí điểm nút dựa trên gradient của trường, giúp làm dày mật độ phần tử ở vùng lân cận dây dẫn và thưa dần ra biên xa. Nhờ đó, số nút tính toán giảm hơn 52% từ 1594 nút xuống 758 nút, giúp tiết kiệm bộ nhớ máy tính và rút ngắn hơn phân nửa thời gian xử lý thuật toán trên MATLAB.

Cường độ điện trường dưới đường dây 500kV Cai Lậy - Ô Môn có gây nguy hại trực tiếp đến sức khỏe người dân không? Kết quả tính toán tại độ cao 1 mét cho thấy cường độ điện trường cực đại đạt 4,2 kV/m đối với mạch kép và 5,1 kV/m đối với mạch đơn tại điểm tim tuyến. Các mức này đều nằm trong ngưỡng an toàn 5 kV/m theo quy chuẩn Việt Nam và khuyến cáo quốc tế, không gây hiện tượng ion hóa phá vỡ tế bào sinh học.

Tại sao đường dây mạch kép lại có mức phát xạ trường điện từ thấp hơn đường dây mạch đơn? Đường dây mạch kép khi được bố trí thứ tự pha đảo chiều hoặc đối xứng sẽ tạo ra các vector điện trường và từ trường có chiều ngược nhau giữa hai mạch. Sự chồng chập trường này dẫn đến hiện tượng triệt tiêu tương hỗ một phần năng lượng, làm giảm từ 15% đến 20% cường độ trường tổng hợp đo được tại mặt đất.

Mô hình tính toán của luận văn có thể áp dụng cho các cấp điện áp khác như 110kV hay 220kV không? Hoàn toàn có thể mở rộng. Thuật toán phần tử hữu hạn và cấu trúc chương trình trên MATLAB được thiết lập dưới dạng tham số hóa tổng quát, cho phép người dùng thay đổi linh hoạt cấp điện áp, bán kính dây dẫn, số lượng dây phân pha và tọa độ treo dây để mô phỏng chính xác cho mọi cấp điện áp lưới điện.

Kết luận

  • Xây dựng thành công mô hình toán học giải hệ phương trình vi phân Maxwell và phương trình Helmholtz phục vụ tính toán chính xác trường điện từ tần số 50Hz của đường dây truyền tải điện.
  • Phát triển và ứng dụng hiệu quả giải thuật lưới Delaunay thích nghi trên môi trường MATLAB, giúp giảm hơn 52% số nút tính toán và tăng tốc độ xử lý số học đáng kể.
  • Mô phỏng định lượng phân bố cường độ điện trường và cảm ứng từ của tuyến 500kV Cai Lậy - Ô Môn tại các cao độ từ 1 mét đến 13 mét, chứng minh tính an toàn theo quy chuẩn hiện hành.
  • Xác định cấu hình đường dây mạch kép giúp giảm khoảng 18% cường độ điện trường tại mặt đất so với mạch đơn nhờ hiện tượng bù trừ vector trường giữa các pha.
  • Đề xuất hệ thống giải pháp kỹ thuật khả thi gồm tối ưu hóa hoán vị pha, nâng độ tĩnh không và xây dựng rào chắn hấp thụ nhằm bảo vệ tối đa môi trường dân sinh.

Đóng góp lớn nhất của luận văn là cung cấp một công cụ mô phỏng số hóa toàn diện, độ chính xác cao cho bài toán thiết kế và kiểm định an toàn điện từ trường của lưới điện siêu cao áp tại Việt Nam. Kế hoạch phát triển tiếp theo trong giai đoạn 12 đến 24 tháng tới là mở rộng mô hình tính toán sang không gian ba chiều có xét đến công trình nhà ở thực tế và phân tích trường nhiệt của hệ thống cáp ngầm cao thế. Các đơn vị quản lý vận hành và tư vấn thiết kế hệ thống điện hãy tích cực ứng dụng giải pháp mô phỏng số này vào thực tiễn để nâng cao chất lượng công trình và bảo đảm phát triển năng lượng bền vững.