Tổng quan nghiên cứu

Trong hệ thống năng lượng hiện đại, việc tối ưu hóa hiệu suất truyền tải trên các đường dây trên không là ưu tiên kỹ thuật hàng đầu. Khi dòng điện xoay chiều tần số 50 Hz chạy qua hệ thống dây dẫn, sự xuất hiện của dòng điện xoáy (dòng điện Foucault) cùng hiệu ứng bề mặt và hiệu ứng gần khiến mật độ dòng phân bố không đồng đều trên tiết diện dẫn điện. Hiện tượng này làm tăng điện trở xoay chiều hiệu dụng từ 10% đến 15% so với điện trở một chiều, trực tiếp gây phát nhiệt cục bộ và gia tăng tổn thất công suất truyền tải trên lưới điện cao thế.

Vấn đề cốt lõi của nghiên cứu bắt nguồn từ việc mô hình hóa toán học các tương tác điện từ phức tạp trong không gian đa dây dẫn. Các phương pháp giải tích truyền thống hoàn toàn bất lực trước hình học biên phức tạp, trong khi phương pháp phần tử hữu hạn (FEM) đơn thuần đòi hỏi mật độ chia lưới rất lớn, làm bùng nổ chi phí tính toán và yêu cầu tài nguyên phần cứng vượt trội. Ngược lại, phương pháp không lưới toàn cục lại dễ dẫn đến hiện tượng ma trận suy biến. Mục tiêu cụ thể của luận văn là phát triển một giải thuật tổ hợp lai giữa phương pháp không lưới cục bộ sử dụng hàm cơ sở xuyên tâm (Kernel RBF) và phương pháp phần tử hữu hạn truyền thống (FEM) để giải hệ phương trình vi - tích phân thế từ vector hai chiều miền thời gian điều hòa.

Phạm vi nghiên cứu tập trung vào hệ thống đường dây truyền tải 110 kV thực tế thuộc Công ty Điện lực Lâm Đồng, điển hình là xuất tuyến đấu nối trạm biến áp 110/22 kV Bảo Lộc 2 và trạm Lộc Châu, sử dụng dây dẫn nhôm lõi thép ACSR 240/32 mm² theo tiêu chuẩn TCVN 5046 trong giai đoạn thực hiện từ tháng 9 năm 2023 đến tháng 6 năm 2024. Ý nghĩa của nghiên cứu thể hiện qua việc nâng cao độ chính xác tính toán với sai số hội tụ được kiểm soát dưới 1.5%, giảm tải thời gian xử lý dữ liệu từ 20% đến 30%, đồng thời cung cấp cơ sở khoa học định lượng để các kỹ sư điện đánh giá và cải tiến phương thức vận hành hệ thống điện thực tế.

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 xây dựng trên nền tảng hệ phương trình Maxwell ở trạng thái biến thiên điều hòa theo thời gian cho môi trường dẫn điện tuyến tính đồng nhất, bỏ qua dòng điện dịch. Trọng tâm là phương trình liên hệ giữa thế từ vector A(x, y) và mật độ dòng điện tổng Jt qua định luật Ampere và định luật cảm ứng Faraday. Mật độ dòng điện tổng được phân tách thành hai thành phần cơ bản: mật độ dòng điện nguồn Js liên quan đến gradient của điện thế vô hướng phi, và mật độ dòng điện xoáy cảm ứng Je = -j * omega * sigma * A, trong đó omega = 2 * pi * f là tần số góc (với f = 50 Hz) và sigma là độ dẫn điện của vật liệu.

Lý thuyết hiệu ứng bề mặt xác định độ sâu thâm nhập tiêu chuẩn delta theo công thức delta = sqrt(2 / (omega * sigma * mu)), phản ánh mức độ suy giảm dòng điện xoáy xuống còn khoảng 37% so với bề mặt. Đối với dây nhôm có độ dẫn điện sigma = 35 * 10^6 S/m và độ từ thẩm tương đối mu_r = 1, độ sâu bề mặt ở tần số 50 Hz đạt khoảng 8.5 mm; trong khi ở kim loại đồng với sigma = 58 * 10^6 S/m, độ sâu này giảm xuống chỉ còn khoảng 6.6 mm. Dây dẫn ACSR được mô hình hóa với cấu trúc bện nhiều sợi xoắn ngược chiều quanh lõi thép chịu lực nhằm phân bố dòng đều hơn và hạn chế hiện tượng tập trung dòng bề mặt.

Lý thuyết xấp xỉ không lưới sử dụng các hàm cơ sở xuyên tâm cục bộ (Local RBF) kết hợp các dạng hàm kinh nghiệm như Multiquadric (MQ), Inverse Multiquadric (IMQ) và Gauss (GA). Thông số hình dạng tối ưu của hàm RBF được xác định thông qua giải thuật kiểm định chéo loại một điểm (LOOCV), đảm bảo tính ổn định và triệt tiêu sai số số trị khi giải hệ phương trình vi phân đạo hàm riêng cấp hai.

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

Nguồn dữ liệu nghiên cứu sử dụng bộ thông số kỹ thuật chuẩn hóa của dây nhôm lõi thép ACSR 240/32 mm² (đường kính ngoài danh định, điện trở suất, độ từ thẩm) kết hợp với sơ đồ mặt bằng và thiết kế hình học cột néo, cột đỡ thực tế từ lưới điện 110 kV tỉnh Lâm Đồng.

Về phương pháp chọn mẫu và rời rạc hóa, toàn bộ miền tính toán hai chiều bao gồm dây dẫn và không gian không khí xung quanh được tạo lưới sơ bộ bằng các phần tử tam giác tuyến tính chuẩn với tổng số 6.120 phần tử và hệ thống 25 nút toàn cục đến hàng nghìn điểm nút rời rạc. Phương pháp chọn mẫu nút cục bộ được thực hiện bằng cách trích xuất 9 điểm nút lân cận xung quanh mỗi điểm nguồn khảo sát, tạo thành các miền ảnh hưởng chồng lấn cục bộ.

Lý do lựa chọn phương pháp tổ hợp Kernel RBF-FE xuất phát từ bản chất của phương trình vi - tích phân điều hòa. Toán tử vi phân Laplace kết hợp số hạng dòng xoáy được giải quyết tối ưu bằng phương pháp không lưới RBF cục bộ mà không cần tích phân ma trận độ cứng phức tạp, trong khi toán tử tích phân dòng điện trên diện tích mặt cắt ngang của dây dẫn được xử lý chính xác tuyệt đối bằng các hàm nội suy bậc một của phương pháp phần tử hữu hạn. Toàn bộ giải thuật được lập trình và giải số trên nền tảng phần mềm MATLAB trong khung thời gian 9 tháng liên tục.

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

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

Quá trình kiểm chứng mô hình trên bài toán điểm chuẩn với lưới 6.120 phần tử tam giác đã khẳng định tính vượt trội của phương pháp Kernel RBF-FE khi đồ thị sai số toàn cục giảm nhanh hơn hẳn so với FEM truyền thống khi tăng số lượng nút rời rạc.

Đối với mô hình một dây dẫn tròn đơn lẻ đặt trong không khí mang mật độ dòng định mức RMS là J0 = 3.0 * 10^6 A/m² ở tần số 50 Hz, phân bố mật độ dòng điện tổng dọc theo đường kính thể hiện rõ hiệu ứng bề mặt: giá trị thấp nhất tại tâm dây dẫn chỉ đạt Jmin = 2.25 * 10^6 A/m², trong khi tại vùng biên ngoài cùng vọt lên mức Jmax = 5.70 * 10^6 A/m², tương đương mức chênh lệch biên độ lên tới 153.3%.

Khi khảo sát hệ thống đa dây dẫn trên đường dây 110 kV Bảo Lộc 2, kết quả mô phỏng trên 11 vị trí đo dọc đường kính từng pha dây dẫn ACSR 240/32 mm² chỉ ra các quy luật đặc thù:

  • Hệ thống ba dây dẫn bố trí thẳng đứng: Phân bố từ trường và dòng xoáy bị lệch trục nghiêm trọng do hiệu ứng gần, pha ở giữa chịu tác động cảm ứng tương hỗ mạnh nhất từ hai pha biên, dẫn đến mật độ dòng xoáy cục bộ cao hơn pha ngoài khoảng 8.4%.
  • Hệ thống sáu dây dẫn mạch kép bố trí dọc: Tương tác điện từ giữa hai mạch làm gia tăng biên độ mật độ dòng xoáy cực đại tại các bề mặt dây dẫn đối diện nhau thêm 11.2% so với trường hợp mạch đơn.
  • Hệ thống sáu dây dẫn mạch kép bố trí tam giác: Tính đối xứng hình học giúp phân tán đường sức từ trường đều hơn, làm giảm đỉnh mật độ dòng điện xoáy cục bộ từ 5.0% đến 7.2% so với cấu hình bố trí thẳng đứng cùng thông số vận hành.

Thảo luận kết quả

Hiện tượng phân bố dòng không đồng đều được giải thích rõ ràng bởi sự tương tác phi tuyến giữa từ thông xoay chiều tự thân và từ thông giao hỗ giữa các pha lân cận. Dòng điện xoáy Je sinh ra có xu hướng ngược chiều với dòng nguồn Js ở phần lõi bên trong và cùng chiều ở lớp biên ngoài, đẩy phần lớn dòng điện ra sát bề mặt ngoài của dây dẫn. Dữ liệu tính toán từ mô phỏng được biểu diễn trực quan thông qua các biểu đồ phân bố đường sức thế từ vector A(x, y) hai chiều dạng đẳng trị và bảng đối sánh mật độ dòng điện tại 11 điểm nút hình học trên tiết diện dây. Các bảng số liệu cho thấy sự tương quan chặt chẽ giữa khoảng cách hình học giữa các pha và biên độ dòng xoáy cảm ứng.

So với các nghiên cứu sử dụng FEM thuần túy trước đây, phương pháp tổ hợp Kernel RBF-FE giảm triệt để kích thước ma trận hệ thống toàn cục. Ma trận thưa nhận được có số điều kiện ổn định, loại bỏ hoàn toàn nguy cơ phân kỳ nghiệm thường gặp ở phương pháp RBF toàn cục có bậc ma trận vượt quá 10^15. Thời gian chạy thuật toán trên máy tính cá nhân tiết kiệm được khoảng 28% trong khi vẫn bảo đảm độ lệch chuẩn kết quả dưới 1.2% so với nghiệm giải tích chuẩn. Kết quả này mang ý nghĩa thực tiễn to lớn trong việc xác định chính xác điện trở xoay chiều RAC, từ đó định lượng chính xác tổn hao công suất trên toàn tuyến dây truyền tải.

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

Dựa trên các kết quả tính toán định lượng dòng điện xoáy, bốn nhóm giải pháp cụ thể được đề xuất nhằm nâng cao hiệu quả kỹ thuật cho hệ thống truyền tải điện:

Thứ nhất, tối ưu hóa cấu hình hình học treo dây trên cột điện. Ban Quản lý Dự án và Phòng Kỹ thuật thuộc Công ty Điện lực Lâm Đồng cần ưu tiên áp dụng thiết kế xà cột dạng tam giác đối xứng thay thế cho dạng thẳng đứng đối với các dự án xây mới hoặc cải tạo đường dây 110 kV có hệ số mang tải thường xuyên trên 80% công suất định mức. Mục tiêu là giảm từ 5% đến 7% độ lệch dòng xoáy cục bộ giữa các pha, thời gian triển khai khuyến nghị trong giai đoạn 2024 - 2026.

Thứ hai, chuẩn hóa quy cách lựa chọn dây dẫn phân đoạn và bện chéo. Các đơn vị tư vấn thiết kế và chế tạo cáp cần kiểm soát nghiêm ngặt bước xoắn bện của các lớp sợi nhôm ngoài theo tiêu chuẩn TCVN 5046, duy trì đường kính mỗi sợi con trong khoảng 4.0 mm đến 5.0 mm để đảm bảo dòng điện xoay chiều 50 Hz tận dụng trọn vẹn độ sâu thâm nhập 8.5 mm, qua đó giữ tỷ số trở kháng RAC/RDC luôn dưới ngưỡng 1.05.

Thứ ba, tích hợp module tính toán số RBF-FE vào hệ thống quản lý lưới điện thông minh. Trung tâm Điều độ Hệ thống điện cần phối hợp với các nhóm nghiên cứu tin học hóa giải thuật tổ hợp RBF-FEM thành công cụ phụ trợ trong phần mềm SCADA/DMS trong vòng 12 đến 18 tháng tới, phục vụ công tác giám sát trực tuyến nhiệt độ dây dẫn và dự báo chính xác tổn thất công suất tác dụng theo tải thực tế.

Thứ tư, điều chỉnh khoảng cách an toàn pha - pha tại các khoảng cột chịu rung động lớn. Đội Quản lý vận hành đường dây cao thế cần thực hiện rà soát định kỳ hàng năm các khoảng cột néo của tuyến 110 kV Bảo Lộc 2, cân chỉnh khoảng cách tương hỗ giữa các dây dẫn đạt ngưỡng tối ưu từ 3.5 m đến 4.5 m để triệt tiêu ảnh hưởng của hiệu ứng gần gây méo dạng phân bố dòng điện.

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

Công trình nghiên cứu mang giá trị thực tiễn và học thuật chuyên sâu, phù hợp với bốn nhóm đối tượng trọng tâm sau:

Kỹ sư thiết kế và quy hoạch hệ thống điện: Luận văn cung cấp công cụ tính toán trường điện từ chính xác cao để các kỹ sư lựa chọn tiết diện dây ACSR, thiết kế khoảng cách pha trên cột đỡ và tối ưu hóa kết cấu hình học tuyến dây 110 kV, 220 kV và 500 kV nhằm hạn chế phát nhiệt.

Cán bộ quản lý vận hành và điều độ lưới điện: Cung cấp phương pháp luận để đánh giá chính xác tổn thất kỹ thuật do hiệu ứng bề mặt và hiệu ứng gần gây ra, từ đó xây dựng biểu đồ mang tải tối ưu và dự báo các điểm quá nhiệt nguy hiểm trên đường dây truyền tải.

Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật điện và Toán ứng dụng: Là tài liệu tham khảo giá trị về cách thiết lập mô hình toán vi - tích phân thế từ vector, kỹ thuật kết hợp phương pháp không lưới cục bộ RBF với FEM, cùng giải thuật tối ưu hóa tham số hình dạng LOOCV để giải các bài toán trường vật lý phức tạp.

Các nhà sản xuất thiết bị điện và cáp dẫn: Cung cấp cơ sở lý thuyết định lượng về dòng điện Foucault và độ sâu thâm nhập để tối ưu hóa quy trình chế tạo lõi bện dẫn điện, thanh cái phân phối công suất lớn và máy biến áp công nghiệp.

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

Dòng điện xoáy gây ra những tác hại chính nào đối với đường dây truyền tải 110 kV trên không? Dòng điện xoáy làm dòng điện xoay chiều bị đẩy dồn ra lớp ngoài dây dẫn do hiệu ứng bề mặt và hiệu ứng gần, khiến diện tích dẫn điện hiệu dụng bị thu hẹp. Hệ quả là điện trở xoay chiều RAC tăng cao hơn từ 10% đến 15% so với điện trở một chiều RDC, làm gia tăng tổn thất công suất nhiệt Joule và gây suy giảm khả năng tải định mức của đường dây.

Tại sao luận văn lại kết hợp phương pháp không lưới RBF với phương pháp phần tử hữu hạn FEM? Phương pháp không lưới RBF cục bộ xuất sắc trong việc giải toán tử vi phân Laplace mà không chịu sự phụ thuộc vào chất lượng chia lưới hình học, trong khi FEM lại xử lý hoàn hảo toán tử tích phân dòng điện trên diện tích mặt cắt dây dẫn. Sự kết hợp Kernel RBF-FE giúp loại bỏ hiện tượng suy biến ma trận, nâng cao độ chính xác và giảm 25% đến 30% thời gian tính toán.

Độ sâu bề mặt của dây dẫn nhôm và dây dẫn đồng ở tần số công nghiệp 50 Hz khác nhau như thế nào? Theo công thức tính toán độ sâu thâm nhập, ở tần số 50 Hz, độ sâu bề mặt của nhôm đạt khoảng 8.5 mm do có độ dẫn điện 35 * 10^6 S/m. Trong khi đó, đồng có độ dẫn điện cao hơn đạt 58 * 10^6 S/m nên độ sâu bề mặt co lại chỉ còn khoảng 6.6 mm, đồng nghĩa với việc dòng xoay chiều trong dây đồng bị tập trung dày đặc hơn sát bề mặt ngoài.

Dây dẫn nhôm lõi thép ACSR khắc phục tác hại của dòng điện xoáy bằng cơ chế nào? Dây ACSR được cấu tạo từ nhiều sợi nhôm nhỏ có đường kính từ 4 mm đến 5 mm bện xoắn chéo nhiều lớp quanh một lõi thép chịu lực. Việc chia nhỏ tiết diện thành nhiều sợi độc lập giúp mỗi sợi đều nằm trọn vẹn trong vùng độ sâu thâm nhập 8.5 mm, triệt tiêu đáng kể dòng điện xoáy khép vòng bên trong và phân phối dòng điện đồng đều hơn trên toàn tiết diện.

Cấu hình mạch kép bố trí tam giác có ưu thế gì vượt trội so với bố trí thẳng đứng trên đường dây 110 kV? Bố trí tam giác tạo ra cấu trúc không gian có tính đối xứng cao về mặt tương hỗ từ trường giữa ba pha. Kết quả tính toán mô phỏng chứng minh cấu hình tam giác giúp làm phẳng mật độ từ thông, giảm đỉnh mật độ dòng điện xoáy cực đại từ 5% đến 7% so với cấu hình thẳng đứng, giảm hiện tượng lệch pha dòng điện và cân bằng nhiệt độ vận hành giữa các dây dẫ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 - tích phân thế từ vector hai chiều mô tả chính xác tương tác dòng điện xoáy trên hệ thống đa dây dẫn đường dây trên không 110 kV.
  • Phát triển và hoàn thiện giải thuật số tổ hợp lai Kernel RBF-FE với stencil 9 điểm cục bộ, kết hợp phương pháp LOOCV để tối ưu hóa tham số hình dạng, đạt tốc độ hội tụ nhanh vượt trội trên lưới 6.120 phần tử kiểm chuẩn.
  • Định lượng chi tiết mức độ gia tăng mật độ dòng điện bề mặt lên tới 153.3% do hiệu ứng bề mặt trên dây dẫn tròn mang tải RMS J0 = 3.0 * 10^6 A/m² ở tần số công nghiệp 50 Hz.
  • Đánh giá toàn diện 4 cấu hình dây dẫn thực tế của lưới điện 110 kV Lâm Đồng, chỉ rõ cấu hình mạch kép tam giác giúp hạ thấp mật độ dòng xoáy từ 5% đến 7% so với dạng bố trí thẳng đứng.
  • Mở đường cho các nghiên cứu tiếp theo mở rộng áp dụng thuật toán Kernel RBF-FE cho lưới siêu cao áp 220 kV - 500 kV, tích hợp đồng thời trường nhiệt - điện từ và cấu trúc sợi bện xoắn 3D thực tế.

Để tiếp cận toàn bộ dữ liệu mô phỏng, bảng ma trận hệ số và mã nguồn giải thuật trên phần mềm MATLAB, quý độc giả và các kỹ sư hệ thống điện có thể liên hệ trực tiếp đơn vị lưu trữ luận văn tại Thư viện Trường Đại học Bách khoa - ĐHQG TP. Hồ Chí Minh để khai thác chi tiết công trình nghiên cứu khoa học này.