Tổng quan về luận án

Sự phát triển vượt bậc của các hệ thống lưới điện thông minh (Smart Grid) và nhu cầu tiêu thụ năng lượng gia tăng trong các thập kỷ qua đã đặt ra yêu cầu nghiêm ngặt về độ tin cậy và tính liên tục trong cung cấp điện. Trong cấu trúc hạ tầng năng lượng, máy biến áp (MBA) là phần tử trọng yếu thực hiện nhiệm vụ truyền tải và phân phối điện năng. Khi một máy biến áp gặp sự cố, hệ thống bảo vệ rơle sẽ tác động cô lập thiết bị, nhưng quá trình gián đoạn cung cấp điện thường kéo theo thiệt hại kinh tế nặng nề, gia tăng chi phí đại tu và gây áp lực lên toàn bộ hệ thống điện.

Luận án tiến sĩ kỹ thuật với đề tài "Xây dựng giải pháp chẩn đoán sự cố trong máy biến áp 3 pha sử dụng các phương pháp xử lý tín hiệu thông minh" do nghiên cứu sinh Đào Duy Yên thực hiện tại Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên, dưới sự hướng dẫn khoa học của PGS.TS. Trần Hoài Linh và PGS.TS. Trần Xuân Minh (chuyên ngành Kỹ thuật điều khiển và tự động hóa, mã số: 9520216), đã giải quyết trực diện bài toán chẩn đoán sớm và phân loại chính xác các khiếm khuyết cơ - điện tiềm ẩn trong máy biến áp phân phối 3 pha.

flowchart TD
    A["Mô hình hóa PTHH 3D trong ANSYS (Maxwell + Mechanical Workbench)"] --> B["Mô phỏng 6 kịch bản vận hành x 3 cấp tải (50%, 80%, 100%)"]
    B --> C["Trích chọn 15 vector đặc trưng đa chiều (Điện từ & Rung động 3 trục x-y-z)"]
    C --> D["Mạng Nơ-ron MLP huấn luyện bằng thuật toán Levenberg-Marquardt"]
    D --> E["Nhận dạng & Phân loại sự cố (Độ chính xác 100%)"]
    F["Thiết bị đo thực nghiệm PSoC CY8C29566 + Cảm biến MPU-6050"] --> G["Kiểm chứng thực địa tại Trạm ĐH Công nghiệp 3 (15/9/2020)"]
    G -.-> E

Bối cảnh khoa học và khoảng trống nghiên cứu (Research Gap)

Trong thực tiễn vận hành và các y văn chuyên ngành quốc tế, các giải pháp chẩn đoán sự cố máy biến áp kinh điển chủ yếu tập trung vào các máy biến áp lực truyền tải dung lượng rất lớn thông qua ba nhóm phương pháp chính:

  1. Phân tích tín hiệu điện (dòng điện, điện áp, sóng hài, đặc tính elip RLC): Thường chỉ phát hiện được khi sự cố chạm chập cuộn dây đã tiến triển ở mức nghiêm trọng (trên 15% đến 30% số vòng dây), kém nhạy đối với các hư hỏng cơ khí vi mô ban đầu.
  2. Phân tích khí hòa tan trong dầu (Dissolved Gas Analysis - DGA): Dựa trên các tiêu chuẩn IEC, tỷ số Dornenburg hay tỷ số Rogers; tuy hiệu quả với sự cố nhiệt phân hay phóng điện hồ quang nhưng chỉ áp dụng được cho máy biến áp dầu, cảm biến khí đắt tiền, độ trôi cao và không thể phản ánh tức thời các biến dạng cơ khí cục bộ.
  3. Phân tích đáp ứng tần số rung (Vibration Frequency Response Analysis - VFRA): Mới dừng lại ở các thử nghiệm ngoại tuyến (off-line), máy biến áp làm việc ở chế độ không tải hoặc quy mô phòng thí nghiệm, thiếu vắng mô hình toán học giải tích kết hợp trường điện từ - cơ học đa vật lý hoàn chỉnh.

Đặc biệt, phần lớn các công trình nghiên cứu trong và ngoài nước ít tập trung vào phân hệ máy biến áp phân phối cấp điện áp 22/0.4kV – đối tượng có số lượng vận hành lớn nhất trên lưới điện nhưng lại tiềm ẩn nguy cơ sự cố cao do biến động phụ tải phức tạp. Khoảng trống nghiên cứu then chốt chính là sự thiếu hụt một phương pháp luận tích hợp, cho phép kết hợp đồng thời tín hiệu điện trường và phản ứng rung động cơ học phi tuyến ba chiều ($x, y, z$) trong điều kiện máy mang tải biến thiên, nhằm nhận dạng các dạng sự cố từ giai đoạn khởi phát vi mô (chỉ chạm 2 vòng dây hoặc nới lỏng cơ cấu cơ khí).

Câu hỏi nghiên cứu và giả thuyết khoa học

Luận án được cấu trúc nhằm giải quyết hệ thống câu hỏi nghiên cứu và kiểm chứng các giả thuyết khoa học tường minh sau:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một mô hình toán học và mô phỏng số 3D chính xác phản ánh đầy đủ sự tương tác phi tuyến giữa trường điện từ và ứng suất - chuyển vị cơ học vỏ máy biến áp phân phối dưới tác động của các khiếm khuyết cơ - điện khác nhau?
  • Câu hỏi nghiên cứu 2 (RQ2): Những thông số đặc trưng đa miền nào (kết hợp dòng điện, điện áp, lực điện từ và phổ chuyển vị rung động ba chiều) có khả năng phân tách tối ưu các trạng thái sự cố tiềm ẩn?
  • Câu hỏi nghiên cứu 3 (RQ3): Cấu trúc mạng nơ-ron nhân tạo và giải thuật tối ưu hóa nào mang lại độ chính xác, tốc độ hội tụ và tính khái quát hóa cao nhất để phân loại các dạng sự cố từ tập dữ liệu trích chọn?
  • Giả thuyết khoa học 1 (H1): Sự xuất hiện của các khiếm khuyết cơ học (lỏng bu lông gá, nới lỏng cuộn dây) và sự cố điện (chập từ 2 vòng dây đến 10% số vòng dây) sẽ làm biến dạng phân bố từ thông tản, gây ra các sóng hài lực điện từ và dẫn tới sự thay đổi dị thường có tính quy luật của vector chuyển vị rung động vỏ máy ở dải tần số 100 Hz và các bội số bậc cao.
  • Giả thuyết khoa học 2 (H2): Việc tích hợp đồng thời 15 thông số đặc trưng điện - cơ vào mô hình mạng nơ-ron Perceptron đa lớp (MLP) tối ưu hóa bằng giải thuật Levenberg - Marquardt sẽ phân loại chính xác tuyệt đối các kịch bản sự cố mà các phương pháp truyền thống đơn biến không thể phân biệt.

Khung lý thuyết và phạm vi nghiên cứu

Nghiên cứu được xây dựng trên nền tảng:

  • Lý thuyết trường điện từ Maxwell tổng quát kết hợp Lý thuyết cơ học môi trường liên tục và Hiện tượng từ giảo (Magnetostriction).
  • Phương pháp phần tử hữu hạn (Finite Element Method - FEM) được hiện thực hóa qua bộ công cụ ANSYS Maxwell 3D và ANSYS Mechanical Workbench.
  • Lý thuyết trí tuệ nhân tạo và nhận dạng mẫu dựa trên mạng nơ-ron Perceptron đa tầng (MultiLayer Perceptron - MLP) với giải thuật huấn luyện Levenberg – Marquardt.

Phạm vi và quy mô thực nghiệm:

  • Đối tượng nghiên cứu: Máy biến áp phân phối 3 pha ngâm dầu dung lượng 400kVA, tỷ số biến áp 22/0.4kV, tổ đấu dây Y/Y0, sử dụng lõi thép kỹ thuật JGH100.
  • Tập kịch bản mô phỏng: 06 trạng thái làm việc (01 trạng thái bình thường và 05 trạng thái sự cố: chập 2 vòng dây cao áp pha B, chập 5% số vòng dây cao áp pha B, chập 10% số vòng dây cao áp pha B, nới lỏng cuộn dây theo phương trục, nới lỏng bu lông gá cuộn dây) được khảo sát ở 3 mức mang tải: 50%, 80% và 100% định mức.
  • Tập dữ liệu: Tổng cộng 234 mẫu dữ liệu số hóa, phân chia thành 180 mẫu huấn luyện (Training Set) và 54 mẫu kiểm tra (Testing Set).
  • Thực nghiệm thực tế: Xây dựng phần cứng thiết bị đo chuyên dụng sử dụng vi điều khiển PSoC CY8C29566 kết hợp cảm biến gia tốc MPU-6050, kiểm chứng rung động thực địa tại Trạm biến áp Đại học Công nghiệp 3 thuộc Công ty Điện lực Thái Nguyên.

Literature Review và Positioning

Tổng quan các nhánh nghiên cứu chính

Tổng quan y văn quốc tế và trong nước về chẩn đoán sự cố máy biến áp cho thấy sự phân hóa thành ba trường phái kỹ thuật chủ đạo:

                                  CÁC TRƯỜNG PHÁI CHẨN ĐOÁN SỰ CỐ MBA

1. Nhóm phương pháp dựa trên phân tích tham số điện và sóng hài

Cách tiếp cận cổ điển sử dụng các đại lượng dòng điện, điện áp và tổng trở để phát hiện hư hỏng. Soares và cộng sự (2008) [11] đã mô hình hóa máy biến áp 345/16kV 250MVA bằng mạch tương đương RLC trên phần mềm PSIM, sử dụng quỹ đạo elip quan hệ giữa hiệu điện áp sơ - thứ cấp và dòng điện pha đầu vào. Kết quả chỉ ra khi cuộn cao áp chập 30% và 60% số vòng, góc lệch trục lớn của elip biến thiên lần lượt 11.71% và 16.36%, độ lệch tâm dịch chuyển 0.08% và 0.21%.

Gomes và các cộng sự (2012) [19] áp dụng Support Vector Machines (SVM) để nhận dạng sự cố chập vòng dây trên lưới 100MVA với 144 kịch bản sự cố (tải từ 60% đến 90%, điện trở ngắn mạch 0.1Ω đến 5Ω, tỷ lệ chập từ 1% đến 25%), đạt độ chính xác phân loại dao động từ 92.2% đến 96.96%.

Khaled và cộng sự (2013) [25] áp dụng biến đổi Wavelet rời rạc (DWT) kết hợp mạng nơ-ron xác suất (PNN) trên máy biến áp 115/23kV 50MVA, phân biệt sự cố ngắn mạch ngoài (độ chính xác 83.3%) và chập vòng trong (độ chính xác 96.3% khi số vòng chập từ 10% đến 90%).

Moradzadeh và cộng sự (2014) [60] sử dụng Wavelet Daubechies bậc 4 (db4) để bóc tách phổ năng lượng dòng từ hóa quá độ trên máy biến áp 125V/24V qua dao động ký TDS2024B, cảnh báo sự cố khi sai lệch phổ vượt quá 1%.

2. Nhóm phương pháp phân tích nồng độ khí hòa tan trong dầu (DGA)

DGA được đánh giá là tiêu chuẩn công nghiệp hiệu quả cho chẩn đoán suy thoái cách điện dầu. Các tác giả Duval (2002) [40], Abu-Elanien (2010) [44] và các nghiên cứu trong nước của Nguyễn Hữu Công (2009) [3] đã hệ thống hóa các phương pháp tỷ số khí:

  • Phương pháp tỷ số Dornenburg [24, 28]: Sử dụng 4 tỷ số của 5 loại khí ($R_1 = \text{CH}_4/\text{H}_2$, $R_2 = \text{C}_2\text{H}_2/\text{C}_2\text{H}_4$, $R_3 = \text{C}_2\text{H}_2/\text{CH}_4$, $R_4 = \text{C}_2\text{H}_6/\text{C}_2\text{H}_2$) để định danh các dạng quá nhiệt, phóng điện vầng quang và phóng hồ quang.
  • Phương pháp tỷ số Rogers/IEC: Chuẩn hóa ngưỡng nồng độ tới hạn của các hydrocarbon phân hủy ($\text{H}_2: 100\text{ ppm}$, $\text{CH}_4: 120\text{ ppm}$, $\text{CO}: 350\text{ ppm}$, $\text{C}_2\text{H}_2: 35\text{ ppm}$, $\text{C}_2\text{H}_4: 50\text{ ppm}$).

3. Nhóm phương pháp phân tích đáp ứng tần số rung (VFRA)

Nhánh nghiên cứu cơ học dựa trên hiện tượng biến dạng cơ khí gây dịch chuyển phổ rung vỏ máy. Tang và cộng sự (2017) [38, 41] khảo sát máy biến áp 220kVA trong dải tần số 230Hz – 626Hz, thiết lập chỉ số độ lệch tuyệt đối trung bình (Average Absolute Deviation - AAD): $$\text{AAD} = \frac{1}{N} \sum_{f=f_s}^{f_e} \frac{|A_1(f) - A_0(f)|}{\max(A_0(f), A_1(f))}$$ Ji và cộng sự (2018) [31] sử dụng ANSYS mô phỏng máy biến áp 220kV khi cuộn hạ áp pha C bị nới lỏng hướng trục, chứng minh rằng sự suy giảm lực ép cuộn dây làm xuất hiện các đỉnh cộng hưởng rung động đặc trưng tại các điểm đo trên bề mặt thùng dầu.

Nghiên cứu Phương pháp Đối tượng Độ chính xác / Kết quả Hạn chế chính
Soares et al. (2008) Đặc tính elip RLC (PSIM) MBA 345/16kV 250MVA Lệch trục 11.71% - 16.36% Chỉ phát hiện sự cố lớn (chập 30%-60% vòng), không đo rung động
Gomes et al. (2012) Tín hiệu dòng pha + SVM MBA 100MVA 92.2% – 96.96% Kém nhạy với sự cố cơ khí thuần túy (lỏng bulông, nới lỏng dây)
Khaled et al. (2013) DWT + Mạng nơ-ron PNN MBA 115/23kV 50MVA 96.3% (sự cố trong) Yêu cầu số vòng chập tối thiểu >10%, không phản ánh trạng thái mang tải thấp
Tang et al. (2017) VFRA miền tần số (AAD) MBA 220kVA Nhận diện biến dạng cuộn Thử nghiệm không tải, không kết hợp với mô hình đa vật lý điện từ
Luận án (Đào Duy Yên, 2021) FEM 3D (ANSYS) + MLP (Levenberg-Marquardt) MBA 22/0.4kV 400kVA 100% (180 mẫu học, 54 mẫu kiểm tra) Phát hiện chập siêu vi mô (2 vòng) và lỏng cơ khí dưới 3 cấp tải định mức

Định vị học thuật và những tiến bộ vượt bậc

So sánh với các nghiên cứu tiền nhiệm, luận án của Đào Duy Yên đã tạo ra một bước tiến vượt bậc mang tính hệ thống:

  • Khắc phục sự phân tách rời rạc giữa chẩn đoán điện và chẩn đoán cơ bằng cách xây dựng một khung mô phỏng ghép nối điện từ - cấu trúc hoàn chỉnh trên ANSYS Maxwell và Mechanical.
  • Vượt qua ngưỡng giới hạn phát hiện sự cố của Gomes (chập >1%) và Khaled (chập >10%) khi mô hình hóa thành công và nhận diện chính xác sự cố chập chỉ 02 vòng dây cuộn cao áp pha B.
  • Bổ sung trường hợp khảo sát thực nghiệm khi máy biến áp phân phối làm việc đầy tải và mang tải biến thiên ngoài hiện trường, thay vì chỉ giới hạn ở trạng thái không tải như các công trình của Tang hay Ji.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và hợp nhất các lý thuyết vật lý nền tảng trong kỹ thuật điện và cơ học kết cấu:

  1. Lý thuyết Điện từ học Maxwell: Khảo sát bài toán trường điện từ biến thiên phi tuyến theo thời gian trong không gian 3 chiều: $$\nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t}$$ $$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$$ $$\nabla \cdot \mathbf{B} = 0, \quad \nabla \cdot \mathbf{D} = \rho$$ Thiết lập mối quan hệ giữa thế vector từ $\mathbf{A}$ ($\mathbf{B} = \nabla \times \mathbf{A}$) và mật độ dòng điện $\mathbf{J}$ qua ma trận tenxơ từ trở $[v]$: $$\mathbf{J} = \nabla \times ([v] \nabla \times \mathbf{A})$$
  2. Lý thuyết Lực từ và Tenxơ ứng suất Maxwell: Lực điện động phân bố tác dụng lên cuộn dây và lõi từ được xác định thông qua tích có hướng của mật độ dòng điện và cảm ứng từ $\mathbf{F} = \int_V (\mathbf{J} \times \mathbf{B}) dV$, kết hợp tenxơ ứng suất Maxwell $[T]$ trên bề mặt phần tử: $$T_{ij} = \mu_0 \left( H_i B_j - \frac{1}{2} \delta_{ij} (\mathbf{H} \cdot \mathbf{B}) \right)$$
  3. Lý thuyết Từ giảo và Cơ học động lực học kết cấu: Hiện tượng từ giảo trong các lá thép kỹ thuật JGH100 của lõi từ tạo ra biến dạng dài tỷ lệ bậc hai với cảm ứng từ $(\Delta l / l \approx k B^2 \approx k' U^2)$, dẫn tới tần số rung động cưỡng bức cơ bản luôn gấp đôi tần số lưới điện ($f_{\text{rung}} = 2 f_{\text{lưới}} = 100\text{ Hz}$ khi $f = 50\text{ Hz}$).
[Phương trình Maxwell 3D: J = ∇ × (v ∇ × A)]
[Tenxơ ứng suất Maxwell T_ij & Lực Lorentz: F = J × B]
[Phương trình Động lực học Cơ học: [M]{Ü} + [C]{Ů} + [K]{U} = {Fa}]
[Trích chọn 15 vector đặc trưng biến dạng 3 trục x, y, z]
[Mạng nơ-ron MLP + Cập nhật trọng số Levenberg-Marquardt]

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp phương trình liên kết đa vật lý giữa kích thích điện từ và phản ứng động lực học vỏ máy biến áp: $$[M]{\ddot{U}} + [C]{\dot{U}} + [K]{U} = {F_a}$$ Trong đó:

  • $[M]$: Ma trận khối lượng kết cấu hệ thống vỏ, lõi và cuộn dây.
  • $[C]$: Ma trận cản động lực học kết cấu.
  • $[K]$: Ma trận độ cứng của hệ đàn hồi tuyến tính.
  • ${U}$: Vector chuyển vị nút phần tử trong không gian 3 chiều ($u_x, u_y, u_z$).
  • ${F_a}$: Vector ngoại lực điện từ nút được chuyển tiếp trực tiếp từ kết quả tính toán tenxơ lực Maxwell trong ANSYS Maxwell.

Điều kiện biên và tính toàn vẹn hình học:

  • Biên từ trường: Thiết lập thế vector từ $\mathbf{A} = 0$ tại bề mặt bao quanh không gian mô phỏng (Balloon boundary condition).
  • Biên cơ học: Ràng buộc cố định (Fixed Support) tại 4 vị trí chân đế gá đặt máy biến áp trên sàn trạm; liên kết tiếp xúc đàn hồi giữa cuộn dây, lõi thép và vỏ máy qua môi trường truyền dẫn dầu cách điện.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) và chủ nghĩa hiện thực phản biện (Critical Realism), kết hợp giữa mô phỏng số học giải tích bậc cao và kiểm chứng thực nghiệm tại hiện trường.

Quy trình nghiên cứu chi tiết

  1. Thiết lập mô hình phần tử hữu hạn:
    • Sử dụng các phần tử tứ diện bậc hai 8 nút và phần tử lập phương phi tuyến 20 nút để mô tả chính xác dạng hình học phức tạp của 3 trụ lõi từ và các bối dây đồng hình xuyến.
    • Gán đường đặc tính từ hóa vật liệu phi tuyến $B-H$ của thép lá kỹ thuật điện JGH100.
    • Thiết lập các trường hợp sự cố hình học: mô hình nới lỏng cuộn dây bằng cách mở rộng khe hở không gian giữa các lớp dây; mô hình ngắn mạch bằng cách cô lập điện áp và ngắn mạch trực tiếp 2 vòng dây, 5% và 10% số vòng dây cuộn cao áp pha B.
  2. Quy trình trích chọn đặc trưng: Từ kết quả mô phỏng quá độ theo thời gian (Time-transient analysis) với bước tích phân thời gian siêu nhỏ ($\Delta t = 0.0001\text{s}$), 15 thông số đặc trưng được trích xuất tự động bao gồm:
    • Các đại lượng điện: Dòng điện hiệu dụng 3 pha sơ cấp và thứ cấp ($I_{1A}, I_{1B}, I_{1C}, I_{2a}, I_{2b}, I_{2c}$).
    • Các thành phần lực điện từ: Lực hướng trục và lực xuyên tâm tác động lên cuộn dây và gông từ lõi thép theo 3 trục tọa độ ($F_x, F_y, F_z$).
    • Các đại lượng chuyển vị cơ học: Biên độ chuyển vị dao động cực đại và độ lệch chuẩn trên bề mặt vỏ máy biến áp theo 3 trục ($U_x, U_y, U_z$).

Dữ liệu và giải thuật huấn luyện mạng nơ-ron

Tập dữ liệu gồm 234 mẫu được chuẩn hóa phi tuyến về dải $[0, 1]$. Cấu hình mạng nơ-ron MLP được khảo sát với cấu trúc:

  • Lớp đầu vào: 15 nơ-ron tương ứng với 15 đặc trưng trích chọn.
  • Lớp ẩn: Khảo sát số lượng nơ-ron ẩn tăng dần từ 1 đến 5 để đánh giá dung lượng học và hiện tượng quá khớp (Overfitting).
  • Lớp đầu ra: 01 nơ-ron đại diện cho trạng thái làm việc (gán nhãn mã hóa từ trạng thái bình thường đến các dạng sự cố cụ thể).
  • Giải thuật huấn luyện: Thuật toán Levenberg – Marquardt (LM) với công thức cập nhật véc-tơ trọng số $\mathbf{w}$: $$\mathbf{w}_{k+1} = \mathbf{w}_k - \left[ \mathbf{J}^T \mathbf{J} + \mu \mathbf{I} \right]^{-1} \mathbf{J}^T \mathbf{e}$$ Trong đó: $\mathbf{J}$ là ma trận Jacobian chứa các đạo hàm riêng bậc nhất của sai số mạng theo trọng số; $\mathbf{e}$ là vector sai số đầu ra; $\mu$ là tham số giảm chấn tự thích nghi (khi $\mu$ nhỏ thuật toán tương đương Gauss-Newton, khi $\mu$ lớn thuật toán tương đương bước giảm cực đại Gradient Descent).

Phát hiện đột phá và implications

Những phát hiện then chốt

Dữ liệu mô phỏng và đo lường thực nghiệm mang lại 5 phát hiện mang tính đột phá khoa học:

  1. Hiệu năng nhận dạng đạt độ chính xác tuyệt đối 100%: Trong toàn bộ quá trình thử nghiệm cấu hình mạng MLP với số nơ-ron ẩn từ 1 đến 5, mạng nơ-ron với cấu trúc 15 đầu vào - 3 nơ-ron ẩn - 1 đầu ra (15-3-1) cho kết quả tối ưu vượt trội. Toàn bộ 180 mẫu huấn luyện và 54 mẫu kiểm nghiệm độc lập đều đạt sai số nhận dạng nhỏ hơn ngưỡng dung sai $0.5$, đạt tỷ lệ phân loại chính xác $100%$ mà không xảy ra hiện tượng quá khớp.
              HIỆU NĂNG NHẬN DẠNG CỦA CÁC KIẾN TRÚC MẠNG MLP
               1 nơ-ron  2 nơ-ron  3 nơ-ron  4 nơ-ron  5 nơ-ron
                             (Cấu hình 15-3-1 tối ưu)

  1. Quy luật biến dạng từ trường và lực cơ học phi đối xứng khi chập vòng dây: Khi xảy ra sự cố chập chỉ 2 vòng dây trên cuộn cao áp pha B, dòng điện quẩn nội tại trong các vòng chập tăng vọt, gây méo cục bộ mật độ từ thông tản. Kết quả mô phỏng trên ANSYS cho thấy lực điện động hướng kính tác dụng lên cuộn pha B mất tính đối xứng, làm biên độ chuyển vị rung động theo phương $y$ trên vỏ máy biến áp tăng vọt gấp 2.8 lần so với chế độ vận hành bình thường.

  2. Sự dịch chuyển phổ rung động của sự cố cơ học thuần túy: Đối với sự cố nới lỏng bu lông gá và nới lỏng cuộn dây theo phương trục, các đại lượng điện áp và dòng điện đầu cực gần như không đổi (dòng tải không biến thiên vượt ngưỡng kích hoạt rơle quá dòng). Tuy nhiên, chuyển vị rung động theo phương $z$ (phương thẳng đứng dọc trục cuộn dây) thể hiện sự suy giảm tần số dao động tự nhiên, làm gia tăng mạnh mẽ biên độ hài bậc 2 ($200\text{ Hz}$) và hài bậc 3 ($300\text{ Hz}$).

  3. Tương quan tuyến tính giữa công suất mang tải và biên độ rung động thực địa: Thực nghiệm tại Trạm ĐH Công nghiệp 3 vào ngày 15/9/2020 ghi nhận dữ liệu đo rung bằng cảm biến MPU-6050 gắn trên vỏ máy biến áp và dữ liệu đo xa công suất tác dụng tổng ($P_{\text{tổng}}$) từ hệ thống SCADA của Điện lực Thái Nguyên. Sau khi chuẩn hóa về dải $[0, 1]$, đồ thị biến thiên độ rung bám sát hoàn hảo theo đồ thị phụ tải $P_{\text{tổng}}$ trong dải từ 9h00 đến 17h00, minh chứng rằng phản ứng rung động tỷ lệ thuận chặt chẽ với bình phương dòng điện phụ tải định mức $I^2$.

Biên độ
chuẩn hóa
           09:00      11:00      13:00      15:00      17:00 (Thời gian)
  1. Tính nhạy vượt trội so với các công trình quốc tế: So với phương pháp elip RLC của Soares (2008) chỉ nhận biết khi chập $\ge 30%$ số vòng, và giải thuật SVM của Gomes (2012) có độ chính xác $92.2% - 96.96%$, việc kết hợp 15 đặc trưng điện - cơ đa vật lý của luận án cho phép phát hiện sự cố chập vòng ngay từ mức $\approx 0.1%$ số vòng (2 vòng dây trên tổng số hàng nghìn vòng dây của cấp điện áp 22kV) với độ tin cậy $100%$.

Implications đa chiều

  • Về mặt lý thuyết: Luận án đã thiết lập thành công mô hình tích hợp trường điện từ - cấu trúc giải tích phần tử hữu hạn cho MBA phân phối, chứng minh rằng sự kết hợp giữa biến dạng cơ học cục bộ và phân bố lực Maxwell tạo ra một không gian vector đặc trưng có khả năng phân tách trạng thái sự cố vượt trội so với các mô hình mạch đẳng trị RLC tập trung.
  • Về mặt phương pháp luận: Cung cấp một quy trình khép kín từ mô phỏng số học FEM đa vật lý để tạo mẫu huấn luyện (giải quyết triệt để bài toán khan hiếm dữ liệu sự cố thực tế trong công nghiệp) đến tối ưu hóa mạng nơ-ron nhận dạng và tích hợp hệ thống phần cứng đo lường nhúng.
  • Về mặt thực tiễn vận hành: Cho phép các đơn vị truyền tải và phân phối điện chuyển dịch từ chiến lược bảo dưỡng định kỳ thụ động (Time-Based Maintenance - TBM) sang bảo dưỡng dự đoán theo tình trạng thực tế (Condition-Based Maintenance - CBM), giúp ngăn ngừa các sự cố nổ vỡ máy biến áp, giảm thiểu thời gian mất điện SAIDI/SAIFI.

Limitations và Future Research

Nghiên cứu duy trì tính khách quan khoa học khi thẳng thắn thừa nhận các giới hạn kỹ thuật nội tại:

  1. Phụ thuộc vào dữ liệu mô phỏng số cho các kịch bản sự cố nguy hiểm: Do yêu cầu an toàn nghiêm ngặt của hệ thống điện, các kịch bản sự cố chạm chập vòng dây phá hủy không thể thực hiện thực nghiệm cưỡng bức trên máy biến áp đang vận hành tại hiện trường, mà chủ yếu dựa trên bộ tham số mô phỏng số học chính xác cao trong ANSYS.
  2. Ảnh hưởng của môi trường dầu cách điện và sóng âm: Mô hình phần tử hữu hạn hiện tại ghép nối trường điện từ và biến dạng cơ học rắn của lõi/cuộn dây/vỏ máy, nhưng chưa tích hợp đầy đủ tương tác chất lưu - kết cấu (Fluid-Structure Interaction - FSI) phản ánh sự truyền dẫn sóng âm và dòng chảy đối lưu nhiệt của dầu máy biến áp.
  3. Giới hạn phần cứng cảm biến gia tốc: Cảm biến MPU-6050 sử dụng công nghệ vi cơ điện tử MEMS có dải tần số đo hữu hạn và nhạy cảm với nhiệt độ môi trường cao trên bề mặt thùng dầu nếu không được bù nhiệt độ chuyên dụng.

Định hướng nghiên cứu phát triển

  • Tích hợp mô hình ghép nối 3 trường (Điện từ - Nhiệt thủy lực - Cơ học rung động): Xây dựng mô hình FSI mô phỏng đầy đủ ảnh hưởng của độ nhớt và đối lưu dầu biến áp đến quá trình truyền sóng rung động từ lõi/cuộn dây ra thành vỏ.
  • Ứng dụng học sâu (Deep Learning) và Mạng nơ-ron tích chập (CNN): Ứng dụng biến đổi biểu đồ thời gian - tần số (như CWT, STFT) của tín hiệu rung để tạo ảnh phổ (Spectrogram), sau đó sử dụng mạng CNN 2D nhận dạng tự động không cần qua bước trích chọn đặc trưng thủ công.
  • Phát triển thiết bị giám sát thông minh Edge AI trên nền tảng IoT: Tích hợp các thuật toán nén mạng nơ-ron trực tiếp lên vi xử lý nhúng (như STM32, ESP32 hoặc FPGA) kết hợp cảm biến rung áp điện công nghiệp chuẩn cách điện cao, truyền dữ liệu không dây thời gian thực qua chuẩn LoRaWAN / 4G-5G về trung tâm điều độ SCADA.
  • Mở rộng dải đối tượng nghiên cứu: Mở rộng mô hình toán và cơ sở dữ liệu nhận dạng cho các dòng máy biến áp lực 110kV, 220kV và 500kV có cấu trúc lõi từ và bệ đỡ phức tạp.

Tác động và ảnh hưởng

                               CÁC TRỤ CỘT TÁC ĐỘNG CỦA NGHIÊN CỨU
  • Tác động học thuật: Đóng góp một phương pháp luận mẫu mực trong việc ứng dụng công cụ mô phỏng phần tử hữu hạn đa vật lý kết hợp kỹ thuật trí tuệ nhân tạo để giải quyết bài toán chẩn đoán kỹ thuật trong ngành Kỹ thuật Điều khiển và Tự động hóa.
  • Chuyển giao công nghiệp: Mở ra giải pháp chế tạo các module giám sát độ rung trực tuyến gắn ngoài (Non-invasive Online Monitoring) với chi phí sản xuất thấp hơn từ 5 đến 10 lần so với các hệ thống phân tích DGA nhập ngoại đắt đỏ, dễ dàng triển khai đại trà cho hàng vạn trạm biến áp phân phối thuộc Tập đoàn Điện lực Việt Nam (EVN).
  • Lợi ích kinh tế - xã hội: Việc phát hiện sớm nguy cơ chập vòng và nới lỏng cơ khí giúp kéo dài tuổi thọ máy biến áp thêm từ 5 đến 10 năm, ngăn chặn kịp thời các thảm họa cháy nổ trạm điện, bảo vệ an toàn cho con người và giảm thiểu tối đa tổn thất kinh tế do gián đoạn điện năng trong sản xuất công nghiệp.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Điện - Tự động hóa: Thụ hưởng một khung lý thuyết hoàn chỉnh, các công thức giải tích trường - lực tương hỗ và bộ chương trình huấn luyện mạng nơ-ron MLP trên MATLAB để tiếp tục mở rộng các hướng nghiên cứu chuyên sâu.
  • Kỹ sư vận hành và điều độ hệ thống điện (EVN / Trung tâm Điều độ HTĐ Quốc gia): Tiếp cận một công cụ chẩn đoán trực quan, tin cậy, hỗ trợ ra quyết định cô lập và sửa chữa thiết bị kịp thời dựa trên dữ liệu đo lường định lượng thay vì cảm tính.
  • Các nhà sản xuất máy biến áp trong nước (EEMC, ABB, THIBIDI, Đông Anh): Sử dụng mô hình mô phỏng ANSYS Maxwell và Mechanical của luận án để tối ưu hóa thiết kế kết cấu cơ khí, bulông gá và lực ép gông từ, hạn chế tối đa độ ồn và rung động ngay từ khâu thiết kế chế tạo ban đầu.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập mô hình liên kết giải tích phi tuyến giữa bài toán trường điện từ Maxwell 3D và bài toán động lực học kết cấu đàn hồi Navier-Cauchy thông qua tenxơ lực Maxwell để giải thích bản chất biến dạng phổ rung động $100\text{ Hz}$ và các sóng hài bậc cao trên vỏ máy biến áp phân phối. Luận án đã mở rộng lý thuyết điện từ học cổ điển từ chỗ chỉ phân tích tham số điện sang lý thuyết tương tác đa vật lý (Electro-Magneto-Mechanical Interaction), chứng minh rằng các biến dạng vi mô của cuộn dây ($0.1%$ số vòng) tạo ra sự tái phân bố lực từ trường cục bộ đủ lớn để phát hiện thông qua chuyển vị không gian 3 chiều.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu trước đây được thể hiện như thế nào?

So với các nghiên cứu của Gomes (2012) [19] (chỉ dùng tín hiệu dòng điện) và Tang (2017) [38] (chỉ đo đáp ứng tần số rung VFRA đơn biến), luận án đã:

  1. Đột phá trong việc tổng hợp vector đặc trưng đa miền 15 chiều kết hợp đồng thời cả biến điện từ và biến động lực học cơ học 3 trục ($x, y, z$).
  2. Ứng dụng công cụ chia lưới phần tử hữu hạn 3D với các phần tử bậc cao (8 nút tứ diện cong và 20 nút lập phương) trên ANSYS để giải quyết trọn vẹn bài toán thiếu hụt mẫu dữ liệu sự cố thực nghiệm.
  3. Tích hợp giải thuật Levenberg – Marquardt giúp mạng MLP hội tụ nhanh chóng và đạt độ chính xác phân loại tuyệt đối $100%$ so với mức $92.2% - 96.9%$ của SVM truyền thống.

3. Phát hiện bất ngờ và thú vị nhất thu được từ dữ liệu là gì?

Phát hiện thú vị nhất là sự suy giảm mạnh mẽ của tính đối xứng lực điện từ theo phương hướng kính khi xảy ra sự cố chập chỉ 02 vòng dây. Dù số lượng vòng dây bị chạm chiếm tỷ lệ cực kỳ nhỏ trên tổng số vòng dây cuộn cao áp 22kV, lực điện từ cục bộ sinh ra do dòng ngắn mạch nội tại đã làm thay đổi hoàn toàn vector chuyển vị trên vỏ máy biến áp theo phương ngang ($y$), tạo nên một dấu vân tay rung động (Vibration Signature) rõ rệt mà các cảm biến dòng điện và rơle bảo vệ truyền thống hoàn toàn bỏ qua.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp một quy trình tái lặp vô cùng tường minh tại phần phụ lục, bao gồm:

  • Toàn bộ thông số kích thước hình học, thông số vật liệu thép silic JGH100 (đường cong $B-H$, hệ số tổn hao lõi), thông số dây đồng cách điện của máy biến áp phân phối 400kVA 22/0.4kV Y-Y0.
  • Mã nguồn chương trình trích chọn đặc trưng từ kết quả ANSYS Maxwell/Mechanical.
  • Mã nguồn kịch bản huấn luyện và kiểm thử mạng nơ-ron MLP trên môi trường MATLAB Neural Network Toolbox.
  • Sơ đồ nguyên lý mạch điện tử phần cứng chi tiết của hệ thống đo rung (khối vi điều khiển PSoC CY8C29566, cảm biến MPU-6050, thời gian thực DS1307, mạch nguồn, màn hình LCD và giao tiếp truyền thông PL2303).

5. Kế hoạch nghiên cứu 10 năm tiếp theo từ nền tảng luận án được định hình ra sao?

Kế hoạch nghiên cứu 10 năm được vạch ra theo lộ trình 3 giai đoạn:

  • Giai đoạn 1 (1-3 năm): Hoàn thiện thuật toán bù nhiệt độ cho cảm biến MEMS, tích hợp mô hình tương tác dầu - kết cấu (FSI) 3 pha trên ANSYS Fluent/Mechanical.
  • Giai đoạn 2 (3-6 năm): Thương mại hóa thiết bị đo thông minh chuẩn công nghiệp (Industrial IoT Edge Node), triển khai thử nghiệm trên quy mô 100 trạm biến áp phân phối tại các công ty điện lực trực thuộc EVN.
  • Giai đoạn 3 (6-10 năm): Xây dựng hệ thống cơ sở dữ liệu lớn (Big Data) về chữ ký rung động máy biến áp toàn quốc, ứng dụng các mô hình học sâu Transformer và Generative AI để tự động dự báo tuổi thọ còn lại (Remaining Useful Life - RUL) của thiết bị trong hệ thống lưới điện thông minh quốc gia.

Kết luận

Công trình nghiên cứu luận án tiến sĩ của tác giả Đào Duy Yên đã hoàn thành xuất sắc các mục tiêu khoa học đề ra với các dấu ấn cốt lõi:

  1. Xây dựng thành công mô hình phần tử hữu hạn 3D trên phần mềm ANSYS mô tả chính xác máy biến áp phân phối 400kVA 22/0.4kV Y-Y0 ở chế độ bình thường và 05 kịch bản sự cố cơ - điện điển hình dưới 3 cấp mang tải khác nhau.
  2. Đề xuất bộ vector 15 đặc trưng đa chiều kết hợp giữa tín hiệu điện trường và rung động động lực học 3 trục ($x, y, z$), tạo bước đột phá trong kỹ thuật trích chọn dấu hiệu sự cố.
  3. Thiết kế mạng nơ-ron Perceptron đa lớp (MLP) cấu hình tối ưu 15-3-1 sử dụng thuật toán tối ưu Levenberg – Marquardt, đạt độ chính xác nhận dạng tuyệt đối 100% trên toàn bộ 234 tập mẫu học và kiểm tra.
  4. Chế tạo thành công thiết bị phần cứng đo độ rung thực nghiệm sử dụng vi xử lý PSoC CY8C29566 và cảm biến gia tốc MPU-6050, kiểm chứng thành công mối tương quan chặt chẽ giữa mức độ mang tải $P_{\text{tổng}}$ và biên độ rung động thực tế tại Trạm ĐH Công nghiệp 3 thuộc Điện lực Thái Nguyên.
  5. Mở ra hướng tiếp cận mới trong chẩn đoán không xâm lấn (Non-invasive Diagnosis), đóng góp nền tảng kỹ thuật quan trọng phục vụ quá trình hiện đại hóa và chuyển đổi số hạ tầng lưới điện phân phối thông minh tại Việt Nam.