Tổng quan nghiên cứu

Hệ thống lưới điện phân phối 22kV tại thành phố Tuy Hòa, trực thuộc Công ty Điện lực Phú Yên, bao gồm mạng lưới trải rộng trên địa bàn 12 phường và 4 xã. Do đặc thù đường dây trung áp đi trên không tại vùng Duyên hải Nam Trung Bộ có mật độ giông sét cao, các sự cố do sét đánh trực tiếp và sét lan truyền chiếm từ 40% đến 60% tổng số vụ sự cố thoáng qua lẫn vĩnh cửu hàng năm. Tình trạng quá điện áp khí quyển không chỉ gây gián đoạn cung cấp điện nghiêm trọng, ảnh hưởng trực tiếp đến chỉ số độ tin cậy SAIFI và SAIDI của toàn tỉnh, mà còn tiềm ẩn nguy cơ phá hỏng lớp cách điện đắt tiền của các máy biến áp phân phối.

Vấn đề cốt lõi hiện nay là việc tính toán và bố trí chống sét van trung áp (MVSA) kiểu biến trở oxit kim loại (MOV) trong thực tế vận hành còn bộc lộ nhiều điểm bất cập. Các mô hình tính toán truyền thống hoặc mô phỏng trên các phần mềm như ATP/EMTP thường đòi hỏi nhiều tham số phức tạp khó thu thập, đồng thời thường bỏ qua ảnh hưởng của tổng trở dây nối đất, dẫn đến việc đánh giá sai lệch biên độ điện áp dư thực tế tác động lên đầu cực máy biến áp.

Mục tiêu cụ thể của công trình là xây dựng thành công mô hình chống sét van trung áp theo cấu trúc Pinceti cải tiến trong môi trường Matlab/Simulink có xét đến thông số chiều dài dây nối (từ 1 mét đến 10 mét). Nghiên cứu tiến hành mô phỏng, kiểm chứng độ chính xác với số liệu thực tế từ nhà sản xuất, từ đó đề xuất phương án bố trí vị trí và số lượng MVSA tối ưu cho các trạm biến áp phân phối 1 máy biến áp và cụm 2 máy biến áp 22kV điển hình tại Tuy Hòa. Kết quả nghiên cứu vào tháng 6 năm 2019 đã giúp kiểm soát sai số điện áp dư dưới mức 3,01%, góp phần bảo vệ an toàn cách điện và kéo dài tuổi thọ vận hành của máy biến áp lên trên 15 năm.

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ự tích hợp của hai nền tảng lý thuyết và mô hình chuyên ngành:

Thứ nhất, lý thuyết về biến trở oxit kim loại (MOV) và cơ chế phi tuyến tiếp giáp hạt bán dẫn: Cấu trúc MOV bao gồm các hạt oxit kẽm ZnO dẫn điện tốt với điện trở suất dưới 0,3 Ω.cm, được bao bọc bởi lớp biên mỏng giàu bismuth Bi2O3 có điện trở rất cao. Các tiếp giáp vi mô này hoạt động như các diode Zener đối xứng có rào điện thế từ 2,0V đến 3,5V, tạo nên đặc tính V-I phi tuyến đối xứng dạng hàm mũ. Khi điện áp lưới ở mức bình thường, MOV duy trì trạng thái cách điện với dòng rò chỉ vài microampe; khi có xung sét quá áp, rào điện thế nhanh chóng sụp đổ trong thời gian nano giây để xả dòng xung sét xuống đất.

Thứ hai, mô hình chống sét van trung áp Pinceti: Được phát triển dựa trên mô hình phụ thuộc tần số kinh điển của IEEE Working Group 3.11 nhưng đã được tinh gọn hóa. Mô hình Pinceti loại bỏ nhánh điện dung ký sinh phức tạp và thay thế các điện trở phi tuyến bằng hai phần tử A0 và A1 kết hợp với điện trở song song R = 1 MΩ nhằm triệt tiêu hiện tượng mất ổn định số học trong quá trình giải thuật toán. Giá trị các cuộn cảm L0 và L1 được tự động xác định trực tiếp thông qua điện áp định mức và điện áp dư tại xung dòng 10kA dạng sóng 8/20µs và xung sét nhanh 1/T2.

Các khái niệm chính được ứng dụng xuyên suốt bao gồm: Mức chịu đựng xung cơ bản (BIL), điện áp làm việc liên tục cực đại (MCOV = 15,3kV), điện áp định mức (Ur = 18kV), điện áp dư (Vr), và hệ số phi tuyến alpha theo các tiêu chuẩn kỹ thuật quốc tế IEC 60099-4 và IEC 62848-1.

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

Nguồn dữ liệu của đề tài được thu thập từ hai nguồn chính: Dữ liệu thông số kỹ thuật và kết quả thử nghiệm xung sét từ catalogue chính hãng của hai dòng thiết bị MVSA phổ biến tại Việt Nam là Cooper AZG2 18kV và Ohio Brass EVP 18kV; cùng với thông số cấu trúc hình học, dây dẫn và mặt bằng trạm biến áp phân phối thực tế được cung cấp bởi Điện lực Tuy Hòa.

Cỡ mẫu và quy trình chọn mẫu: Nghiên cứu lựa chọn mẫu khảo sát có chủ đích gồm 2 dòng chống sét van thương mại hàng đầu, thử nghiệm qua 3 mức biên độ dòng xung sét chuẩn gồm 3kA, 5kA và 10kA dạng sóng 8/20µs, kết hợp kiểm chứng trên 2 dạng trạm biến áp phân phối điển hình (trạm treo 1 MBA và trạm 2 MBA tại Công ty Dược). Quy mô mẫu này bảo đảm độ bao phủ hoàn toàn các tình huống vận hành thực tế tại địa phương.

Phương pháp phân tích: Luận văn sử dụng công cụ Curve Fitting Toolbox trong Matlab để giải tích hóa phương trình toán học của dạng sóng xung sét chuẩn 8/20µs. Đồng thời, công cụ Simulink/SimPowerSystems được sử dụng để xây dựng các khối chức năng con Mask Subsystem cho máy phát xung sét và mô hình MVSA Pinceti. Lý do lựa chọn Matlab/Simulink thay vì phần mềm ATP/EMTP là vì tính trực quan, độ phổ biến cao trong đào tạo học thuật tại Việt Nam, cùng khả năng liên kết dữ liệu linh hoạt giúp người dùng chỉ cần nhập trực tiếp thông số từ catalogue mà không phải thực hiện các bước tinh chỉnh lặp tham số phức tạp. Tiến độ nghiên cứu được thực hiện và hoàn tất nghiệm thu vào tháng 6 năm 2019.

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

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

Thứ nhất, mô hình Pinceti xây dựng trên môi trường Matlab chứng minh độ chuẩn xác vượt trội khi mô phỏng phản ứng quá điện áp. Đối với chống sét van Cooper AZG2 cấp điện áp 18kV, mô hình Pinceti ghi nhận sai số điện áp dư mô phỏng so với dữ liệu kiểm định của nhà sản xuất chỉ là 2,30% ở xung 5kA (đạt 45,63kV so với 46,7kV) và 1,65% ở xung 10kA (đạt 49,86kV so với 50,7kV). Trong khi đó, mô hình mặc định của Matlab cho mức sai số cao hơn, dao động từ 1,22% đến 5,48%.

Thứ hai, đối với chủng loại chống sét van Ohio Brass EVP 18kV, ưu thế của mô hình Pinceti càng thể hiện rõ rệt. Sai số điện áp dư của mô hình Pinceti chỉ đạt 1,95% ở mức xung 5kA (đạt 41,67kV so với 42,5kV) và 3,01% ở mức xung 10kA (đạt 44,13kV so với 45,5kV). Ngược lại, mô hình chống sét van mặc định trong thư viện Matlab xuất hiện sai số rất lớn, lệch từ 11,53% ở xung 5kA lên tới 16,50% ở xung 10kA.

Thứ ba, nghiên cứu phát hiện tác động phi tuyến nghiêm trọng của chiều dài dây nối đất chống sét van đến mức điện áp dư tác động lên máy biến áp. Khảo sát với dây cáp đồng tiết diện 35 mm2 (điện trở R0 = 0,643 Ω/km, độ tự cảm L0 = 0,00253 H/km), khi chiều dài dây tăng từ 1 mét lên 10 mét, điện cảm ký sinh đã làm gia tăng điện áp đỉnh tại đầu cực MBA từ 15% đến 35% dưới tác động của dòng sét 10kA, làm suy giảm nghiêm trọng khoảng an toàn cách điện.

Thứ tư, đối với cấu hình cụm trạm 2 máy biến áp dùng chung thiết bị bảo vệ (điển hình tại TBA Công ty Dược 1 và 2), sự chênh lệch khoảng cách nối dây từ MVSA đến từng đầu cực máy biến áp (L1 = 4 mét và L2 = 8 mét) tạo ra sự lệch pha và biên độ quá áp phản xạ lên tới hơn 18% giữa hai máy, đòi hỏi phải tái cấu trúc điểm đấu nối.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự chênh lệch sai số thấp ở mô hình Pinceti là nhờ việc tích hợp hai điện cảm động học L0 và L1. Hai thông số này đại diện cho độ trễ điện từ và đáp ứng cảm ứng của biến trở khi dòng sét tăng vọt với độ dốc sườn trước lớn, điều mà mô hình điện trở phi tuyến thuần túy của Matlab đã bỏ qua.

Khi so sánh với các công trình quốc tế, điển hình như nghiên cứu của Pramuk Unahalekhaka sử dụng mô hình P-K trên phần mềm ATP/EMTP với sai số 5,39%, mô hình Pinceti cải tiến trong luận văn này đạt độ chuẩn xác cao hơn đáng kể với sai số tối đa chỉ 3,01%, thấp hơn rất nhiều so với ngưỡng sai số cho phép 20% trong tiêu chuẩn kỹ thuật điện.

Về mặt trực quan hóa, toàn bộ dữ liệu phản ứng quá độ được thể hiện trực quan qua các bảng thống kê sai số tương đối và đồ thị sóng điện áp theo miền thời gian từ 0 đến 50 microgiây. Các biểu đồ này thể hiện rõ biên độ điện áp kẹp và thời gian phục hồi của van MOV, giúp kỹ sư nhận diện chính xác thời điểm xuất hiện đỉnh xung sét nguy hiểm nhất.

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

Nhằm nâng cao hiệu quả bảo vệ quá áp cho lưới điện trung áp 22kV thành phố Tuy Hòa, tác giả đưa ra 4 nhóm giải pháp kỹ thuật cụ thể:

Thứ nhất, rút ngắn chiều dài dây nối đất của chống sét van. Đội Quản lý Vận hành thuộc Điện lực Tuy Hòa cần thực hiện nắn thẳng và thu ngắn tuyến dây nối từ cực dưới của MVSA xuống hệ thống tiếp địa máy biến áp xuống dưới 2 mét (thay vì để dài tới 4 - 10 mét như nhiều trạm hiện hữu). Giải pháp này giúp triệt tiêu độ tự cảm ký sinh, giảm trực tiếp từ 20% đến 25% biên độ điện áp dư tác động vào cuộn dây máy biến áp, hoàn thành rà soát trong giai đoạn 2019 - 2020.

Thứ hai, chuẩn hóa công cụ tính toán mô phỏng tại đơn vị quản lý. Phòng Kỹ thuật Công ty Điện lực Phú Yên nên đưa khối mô hình MVSA Pinceti trên nền Matlab/Simulink vào quy trình tính toán phối hợp cách điện định kỳ cho 100% các trạm biến áp phân phối xây mới hoặc nâng cấp công suất từ năm 2020.

Thứ ba, tái bố trí cấu hình chống sét cho cụm trạm có từ 2 máy biến áp trở lên. Đối với các trạm có nhiều máy biến áp đặt cách nhau trên 5 mét, không nên sử dụng chung một bộ 3 chống sét van ở đầu nhánh phân phối. Đơn vị thi công cần lắp đặt riêng biệt từng bộ MVSA ngay tại đầu cực cao áp của từng máy biến áp nhằm ngăn ngừa hiện tượng sóng phản xạ quá áp lan truyền, triển khai từ quý 1 năm 2020.

Thứ tư, nâng cấp tiết diện và chất lượng dây nối đất thoát sét. Chuyển đổi toàn bộ dây thoát sét bằng đồng trần từ tiết diện 35 mm2 sang tối thiểu 50 mm2 tại các trạm biến áp cấp điện cho 12 phường trung tâm, nhằm giảm điện trở nối đất xung sét và hạ tỷ lệ sự cố phóng điện cách điện xuống dưới mức 0,5% mỗi năm.

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

Tài liệu luận văn mang lại giá trị thực tiễn và học thuật chuyên sâu cho 4 nhóm đối tượng:

Thứ nhất, kỹ sư thiết kế và quy hoạch lưới điện: Nắm vững phương pháp mô hình hóa Pinceti trên Matlab để tính toán chính xác mức phối hợp cách điện, từ đó lựa chọn đúng thông số MCOV và dòng xả định mức cho các dự án lưới điện 22kV và 35kV.

Thứ hai, cán bộ quản lý kỹ thuật và vận hành tại các công ty điện lực: Có thêm cơ sở khoa học để đánh giá hiện trạng vận hành, rà soát lại chiều dài dây nối đất và cấu hình chống sét cho hơn 500 trạm biến áp phân phối, giảm thiểu sự cố cháy hỏng máy biến áp do sét lan truyền.

Thứ ba, học viên cao học và nghiên cứu sinh ngành Kỹ thuật điện: Sử dụng luận văn như một tài liệu tham khảo mẫu mực về quy trình xây dựng khối giao diện Mask Subsystem trong Simulink, phương pháp giải tích đường cong xung sét 8/20µs và xử lý phần tử phi tuyến.

Thứ tư, các nhà sản xuất và đơn vị kiểm định thiết bị chống sét: Tham khảo mô hình toán học và phương pháp đánh giá sai số điện áp dư thực nghiệm theo tiêu chuẩn IEC 60099-4 phục vụ công tác tối ưu hóa vật liệu biến trở oxit kim loại.

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

Mô hình Pinceti có ưu thế gì nổi bật so với mô hình chống sét van mặc định trong Matlab? Mô hình Pinceti tích hợp thêm hai điện cảm động học L0 và L1 phản ánh chính xác đáp ứng của biến trở với độ dốc sườn trước của xung sét, giúp giảm sai số điện áp dư xuống chỉ còn từ 1,65% đến 3,01%, trong khi mô hình mặc định của Matlab có sai số lệch tới 16,50%.

Tại sao chiều dài dây nối đất của chống sét van lại ảnh hưởng lớn đến mức bảo vệ máy biến áp? Dây dẫn có độ tự cảm ký sinh khoảng 0,00253 H/km. Dưới tác động của xung sét có tốc độ biến thiên dòng điện cực lớn di/dt, sụt áp cảm ứng trên dây nối (L.di/dt) sẽ cộng dồn trực tiếp vào điện áp dư của van, làm tăng tổng điện áp đặt lên cách điện máy biến áp từ 15% đến 35% khi dây dài tới 10 mét.

Tại sao tác giả lại chọn Matlab/Simulink thay vì phần mềm ATP/EMTP chuyên dụng? Matlab/Simulink là phần mềm trực quan, phổ biến trong giảng dạy và nghiên cứu tại Việt Nam. Mô hình Pinceti trên Matlab cho phép nhập trực tiếp các thông số từ catalogue mà không cần lặp thử nghiệm tham số phức tạp như trên ATP/EMTP, mang lại sự tiện lợi vượt trội cho kỹ sư.

Cụm trạm có 2 máy biến áp đặt cạnh nhau có nên dùng chung một bộ chống sét van hay không? Nếu khoảng cách nối dây từ chống sét van đến máy biến áp thứ hai vượt quá 5 mét, hiện tượng sóng phản xạ sẽ làm điện áp đầu cực tăng cao vượt ngưỡng an toàn. Do đó, giải pháp tối ưu là lắp đặt riêng rẽ từng bộ 3 MVSA tại đầu cực của mỗi máy biến áp.

Độ chính xác của mô hình trong luận văn có đáp ứng tiêu chuẩn ngành điện hay không? Kết quả kiểm chứng trên hai dòng sản phẩm Cooper AZG2 và Ohio Brass EVP 18kV cho thấy sai số mô phỏng dao động từ 1,65% đến 3,01%, hoàn toàn đáp ứng xuất sắc yêu cầu kỹ thuật với ngưỡng sai số cho phép dưới 20% theo quy định kiểm định quốc tế.

Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán cấp bách về bảo vệ quá áp cho lưới điện trung áp thành phố Tuy Hòa thông qua 5 đóng góp trọng tâm:

  • Xây dựng thành công mô hình chống sét van trung áp Pinceti có xét đến chiều dài dây nối đất trên phần mềm Matlab/Simulink với giao diện khối trực quan, dễ sử dụng.
  • Kiểm chứng độ chuẩn xác của mô hình với sai số điện áp dư cực thấp, chỉ từ 1,65% đến 3,01%, vượt trội so với mô hình mặc định của Matlab và các nghiên cứu quốc tế trước đây.
  • Xác lập định lượng mức độ ảnh hưởng của chiều dài dây nối đất (từ 1 mét đến 10 mét) đối với sự gia tăng quá điện áp đầu cực máy biến áp.
  • Đề xuất phương án bố trí vị trí và cấu hình lắp đặt MVSA tối ưu cho các trạm phân phối 1 máy biến áp và cụm 2 máy biến áp 22kV thực tế tại Điện lực Tuy Hòa.
  • Đóng góp nguồn tài liệu tham khảo và công cụ mô phỏng giá trị cho công tác đào tạo học thuật sau đại học ngành Kỹ thuật điện cũng như ứng dụng nhân rộng trong toàn bộ Công ty Điện lực Phú Yên giai đoạn 2020 - 2025.

Các đơn vị quản lý vận hành lưới điện phân phối và các nhóm nghiên cứu chuyên ngành được khuyến khích áp dụng ngay mô hình mô phỏng này để tối ưu hóa thiết kế bảo vệ quá áp, nâng cao chất lượng và độ tin cậy cung cấp điện.