Tổng quan nghiên cứu

Trong hệ thống truyền tải điện năng siêu cao áp từ 220kV đến 500kV, các sự cố quá điện áp do sét đánh và thao tác đóng cắt chiếm khoảng 50% đến 60% tổng số sự cố lưới điện, đe dọa trực tiếp đến an toàn vận hành. Do đó, công tác đo lường và ghi nhận chính xác các xung điện áp cao không chu kỳ là bài toán cấp thiết hàng đầu. Trong thực tế, việc đo trực tiếp các thông số điện áp cực lớn gặp nhiều trở ngại kỹ thuật và chi phí đắt đỏ, đòi hỏi phải sử dụng các bộ phân áp để đo gián tiếp. Tuy nhiên, sự xuất hiện của các phần tử ký sinh điện trở, điện cảm và điện dung trong thiết bị đo luôn gây ra hiện tượng méo dạng sóng xung.

Theo tiêu chuẩn quốc tế IEC 60-3, sai số cho phép khi đo giá trị biên độ của xung toàn sóng và xung cắt không được vượt quá 3,0%, trong khi các thí nghiệm chuyên sâu đòi hỏi mức sai số khắt khe hơn từ 1,5% đến 2,0%. Đề tài tập trung nghiên cứu mô hình toán học và khảo sát phổ Wavelet của các dạng xung quá điện áp không chu kỳ tiêu biểu, nhằm xác định chính xác dải tần số đáp ứng của thiết bị đo.

Nghiên cứu được triển khai dựa trên cơ sở dữ liệu thực nghiệm và chuẩn hóa xung điện áp tại khu vực TP. Hồ Chí Minh trong giai đoạn 2010 - 2012. Mục tiêu trọng tâm của luận văn là phát triển giải pháp phân tích phổ tần số hiệu dụng, tối ưu hóa băng thông đo lường và đưa sai số phép đo từ mức 3,0% theo chuẩn quy định giảm sâu xuống chỉ còn 0,1%, cung cấp cơ sở dữ liệu tin cậy cho các nhà sản xuất chế tạo thiết bị đo lường cao áp đạt độ chính xác tối ưu.

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 hệ thống lý thuyết xử lý tín hiệu hiện đại kết hợp mô hình mạch điện quá độ trong không gian hàm bình phương khả tích $L^2(\mathbb{R})$. Hai lý thuyết trụ cột được áp dụng gồm: Lý thuyết biến đổi Fourier (bao gồm biến đổi Fourier liên tục CTFT và rời rạc DTFT) và Lý thuyết biến đổi Wavelet (bao gồm biến đổi liên tục CWT, rời rạc DWT và cấu trúc Wavelet gói đa kênh). Mô hình nghiên cứu thiết lập hàm truyền đạt phức $H(j\omega)$ của bộ phân áp cao áp để liên kết giữa điện áp đầu vào $u_1(t)$ và điện áp đầu ra $u_2(t)$.

Các khái niệm cốt lõi được xây dựng gồm: phổ biên độ - tần số, hàm Wavelet mẹ (tập trung vào họ hàm Mexican Hat, Haar và Daubechies bậc 2 đến bậc 8), phân tích đa phân giải (MRA), lưới nhị tố dyadic với hệ số tỷ lệ $a = 2^m$ và hệ số dịch $b = 2^m n$, cùng hệ thống bộ lọc gương cầu phương 2 kênh theo thuật toán Mallat.

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

Nguồn dữ liệu nghiên cứu sử dụng tập mẫu gồm 5 dạng sóng xung điện áp chuẩn theo tiêu chuẩn IEC 60-3: xung sét chuẩn 1,2/50 µs, xung thao tác chuẩn không chu kỳ có độ dài đầu sóng từ vài chục đến hàng nghìn micro giây, xung thao tác dao động, xung cắt ở đầu sóng và xung cắt ở đuôi sóng. Cỡ mẫu bao gồm hàng nghìn điểm dữ liệu được số hóa với bước lấy mẫu thời gian ở mức micro giây, trích xuất từ các chu trình phóng điện thực nghiệm.

Phương pháp phân tích kết hợp giữa giải tích số và mô phỏng trên nền tảng MATLAB. Nghiên cứu sử dụng phương pháp dây cung để tìm nghiệm phi tuyến và công thức tích phân Simpson bậc 2 để tính toán gần đúng các tích phân phổ. Luận văn lựa chọn biến đổi Wavelet thay vì Fourier cổ điển vì biến đổi Fourier bị hạn chế bởi cửa sổ thời gian cố định, không thể biểu diễn cục bộ các đột biến tức thời có thời gian đầu sóng cực ngắn dưới 1,2 µs. Phương pháp Wavelet cho phép co giãn cửa sổ phân tích linh hoạt theo từng dải tần số. Toàn bộ quy trình tính toán, giải thuật và kiểm chứng mô hình được thực hiện liên tục trong thời gian 24 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, nghiên cứu đã xây dựng thành công điều kiện tồn tại toán học của dạng sóng xung sét thông qua hàm quan hệ $K = T_S / T_{dS}$. Bằng giải thuật dây cung với sai số giới hạn dưới 0,001%, mô hình đã tính toán chính xác các tham số thời gian đầu sóng và đuôi sóng cho các dạng xung cao áp không chu kỳ tiêu biểu.

Thứ hai, kết quả biến đổi Wavelet liên tục (CWT) với hàm cơ sở Mexican Hat trên các hệ số tỷ lệ $a$ biến thiên từ 10 đến $10^5$ (với $b = 1$) đã định vị rõ nét sự phân bố năng lượng tín hiệu. Ở các giá trị $a$ nhỏ dưới 50, phổ Wavelet phản ánh chi tiết các thành phần tần số cao tại sườn trước của xung, trong khi các giá trị $a$ lớn từ $10^3$ đến $10^5$ mô tả chính xác năng lượng ở vùng tần số thấp tại sườn sau của sóng.

Thứ ba, thông qua biến đổi Wavelet nghịch (IDWT) qua mô hình bộ phân áp, nghiên cứu đã xác định được dải tần số cắt thấp $f_t$ và tần số cắt cao $f_c$, thiết lập độ rộng băng thông $\Delta f$ tối ưu. Kết quả này giúp giảm triệt để sai số đo lường xung điện áp từ mức quy chuẩn 3,0% xuống ngưỡng 0,1%, tức nâng cao độ chính xác lên gấp 30 lần so với yêu cầu tiêu chuẩn IEC 60-3.

Thứ tư, tác giả đã thiết lập biểu thức toán học dự đoán nhanh phổ tần số cho xung sét chuẩn, giúp rút ngắn khoảng 60% thời gian tính toán mô phỏng so với các thuật toán phân tích số truyền thống.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự méo dạng tín hiệu khi đo xung cao áp là do các tham số ký sinh trong mạch phân áp tạo ra hiện tượng dịch pha ở tần số trên $10^6$ Hz. Biến đổi Wavelet với đặc tính định xứ không gian - thời gian đã bóc tách hoàn hảo các thành phần dao động quá độ này. So với các công trình xử lý tín hiệu trước đây sử dụng biến đổi Fourier nhanh (FFT) vốn có sai số dao động từ 1,5% đến 2,0% do hiện tượng rò rỉ phổ, phương pháp Wavelet của luận văn đạt độ ổn định vượt trội với sai số 0,1%.

Dữ liệu nghiên cứu được biểu diễn trực quan qua biểu đồ 3D thể hiện sự biến thiên của hệ số Wavelet theo thang tỷ lệ và trục thời gian, kết hợp cùng các bảng thống kê chi tiết giá trị điện áp $u_2^*(t)$ tại các dải tần số cụ thể của 5 dạng sóng chuẩn, minh chứng tính chuẩn xác của mô hình lý thuyết.

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

Thứ nhất, chuẩn hóa quy trình hiệu chuẩn thiết bị đo xung cao áp bằng giải thuật biến đổi Wavelet tại các trung tâm kiểm định và phòng thí nghiệm trọng điểm, hướng tới mục tiêu duy trì sai số đo lường dưới 0,15% trong vòng 6 tháng tới do Trung tâm Kỹ thuật Tiêu chuẩn Đo lường Chất lượng chủ trì.

Thứ hai, tích hợp biểu thức dự đoán nhanh phổ tần và module giải thuật Wavelet vào phần mềm điều khiển của các máy hiện sóng số chuyên dụng, giúp giảm thời gian phân tích dạng sóng từ 15 phút xuống dưới 2 phút trong lộ trình 12 tháng, do các doanh nghiệp chế tạo thiết bị đo lường phối hợp triển khai.

Thứ ba, cải tiến cấu trúc thiết kế mạch phân áp điện dung - điện trở nhằm mở rộng băng thông đáp ứng tần số từ 50 Hz đến trên 10 MHz, giảm thiểu 40% hiện tượng méo xung ở sườn đầu sóng trong vòng 18 tháng, thuộc trách nhiệm của các viện nghiên cứu kỹ thuật điện lực.

Thứ tư, tổ chức định kỳ hàng quý các chương trình đào tạo kỹ thuật chuyên sâu về phân tích dạng sóng quá điện áp cho đội ngũ kỹ sư điều độ và vận hành lưới điện truyền tải 110kV - 500kV, do Tập đoàn Điện lực điều phối thực hiện.

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

Nhóm kỹ sư thiết kế và chế tạo thiết bị đo lường điện cao áp có thể khai thác các thông số băng thông $\Delta f$ và hàm truyền đạt để phát triển các bộ phân áp có độ méo thấp, đạt độ chính xác dưới 0,1%.

Nhóm chuyên viên thí nghiệm và kiểm định tại các phòng thí nghiệm điện cao thế ứng dụng quy trình phân tích phổ Wavelet để đánh giá chính xác các tham số của dạng xung sét 1,2/50 µs và xung thao tác theo tiêu chuẩn IEC 60-3.

Nhóm kỹ sư quản lý vận hành và điều độ hệ thống điện truyền tải 220kV - 500kV sử dụng mô hình toán học để phân tích nguyên nhân sự cố quá điện áp không chu kỳ, phục vụ công tác điều tra sự cố lưới điện thực tế.

Nhóm giảng viên, học viên cao học và nghiên cứu sinh ngành Kỹ thuật Điện sử dụng tài liệu làm cơ sở tham khảo chuyên sâu về kết hợp phương pháp giải tích số (Simpson, dây cung) với lý thuyết Wavelet trong phân tích tín hiệu điện động.

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

Biến đổi Wavelet vượt trội hơn biến đổi Fourier như thế nào trong phân tích xung cao áp? Biến đổi Fourier trải đều thông tin trên toàn trục thời gian nên không thể định vị được các thời điểm biến thiên tức thời. Ngược lại, biến đổi Wavelet sử dụng các khung cửa sổ linh hoạt theo tần số, cho phép phân tích rõ nét các đỉnh nhọn sườn đầu sóng dưới 1,2 µs ở tần số cao và đuôi sóng 50 µs ở tần số thấp.

Luận văn đã cải thiện độ chính xác đo lường xung cao áp so với tiêu chuẩn IEC 60-3 ra sao? Tiêu chuẩn IEC 60-3 quy định sai số đo biên độ xung không được vượt quá 3,0%. Bằng cách ứng dụng biến đổi Wavelet nghịch qua bộ phân áp để tìm chính xác dải tần số đáp ứng, luận văn đã giảm mức sai số đo lường thực tế xuống chỉ còn 0,1%, đáp ứng các phép thử nghiệm đo lường ngặt nghèo nhất.

Tại sao hàm Wavelet Mexican Hat được ưu tiên sử dụng trong nghiên cứu này? Hàm Wavelet Mexican Hat có dạng sóng đối xứng và đạo hàm cấp hai của hàm Gaussian, rất tương đồng với đặc tính biến thiên đỉnh nhọn của xung quá điện áp không chu kỳ. Sự tương thích này giúp phép biến đổi đạt hệ số tương quan cao nhất khi quét các hệ số tỷ lệ $a$ từ 10 đến $10^5$.

Vai trò của phương pháp dây cung và công thức tích phân Simpson trong đề tài là gì? Phương pháp dây cung giải quyết bài toán tìm nghiệm điều kiện tồn tại hàm sóng phi tuyến với sai số dưới 0,001%. Trong khi đó, công thức tích phân Simpson bậc 2 chia đoạn nhỏ được sử dụng để tính toán số trị các biểu thức tích phân phổ phức tạp của tín hiệu, đảm bảo độ chuẩn xác cao cho mô hình.

Biểu thức dự đoán nhanh phổ tần số mang lại giá trị thực tế nào? Biểu thức này giúp kỹ sư xác định ngay các giới hạn tần số cắt $f_t$ và $f_c$ của xung sét chuẩn mà không cần chạy toàn bộ chu trình tính toán tích phân số phức tạp. Giải pháp này giúp rút ngắn khoảng 60% thời gian xử lý, hỗ trợ thiết kế nhanh các bộ lọc chống méo tín hiệu cho thiết bị đo.

Kết luận

  • Xây dựng thành công mô hình toán học và xác lập điều kiện tồn tại cho các dạng xung quá điện áp không chu kỳ trên lưới điện.
  • Ứng dụng xuất sắc biến đổi Wavelet thuận và nghịch với hàm Mexican Hat để khảo sát dải phổ tín hiệu, đưa sai số đo lường từ 3,0% xuống mức 0,1%.
  • Xác định chi tiết các thông số tần số cắt và băng thông đáp ứng cho 5 dạng sóng xung điện áp chuẩn theo tiêu chuẩn IEC 60-3.
  • Phát triển biểu thức dự đoán nhanh phổ tần số, tiết kiệm 60% thời gian tính toán và tối ưu hóa quy trình thiết kế bộ phân áp cao thế.
  • Xây dựng hệ thống mã nguồn thuật toán hoàn chỉnh trên MATLAB làm công cụ tính toán và mô phỏng thực nghiệm.

Hướng phát triển tiếp theo của đề tài là tích hợp giải thuật Wavelet vào các chip xử lý số DSP nhằm tự động hóa quá trình hiệu chuẩn thiết bị đo trên lưới điện thông minh trong giai đoạn 2026 - 2028. Các đơn vị chế tạo và thử nghiệm thiết bị điện cao áp nên ứng dụng ngay giải pháp phân tích phổ Wavelet để nâng cao chất lượng và độ tin cậy của phép đo.