Tổng quan về luận án

Sự chuyển dịch cơ cấu năng lượng toàn cầu từ nhiên liệu hóa thạch sang năng lượng tái tạo đã thúc đẩy sự phát triển bùng nổ của các hệ thống phát điện sức gió (Wind Energy Conversion Systems - WECS). Báo cáo từ Hội đồng Năng lượng Gió Toàn cầu (GWEC) ghi nhận công suất điện gió tích lũy toàn cầu đã vượt mức 651 GW vào cuối năm 2019, đạt mức tăng trưởng 10% so với năm 2018 và duy trì quỹ đạo tăng trưởng theo cấp số nhân. Trong bức tranh công nghệ này, máy phát điện không đồng bộ nguồn kép (Doubly Fed Induction Generator - DFIG) đóng vai trò xương sống khi chiếm gần 50% thị phần máy phát điện gió thương mại toàn cầu, bao phủ dải công suất chủ đạo từ 1.5 MW đến 3.0 MW với hơn 93 model khác nhau từ các nhà sản xuất hàng đầu thế giới. Ưu thế kinh tế - kỹ thuật mang tính quyết định của DFIG nằm ở cấu trúc bộ biến đổi công suất phân đoạn: stator nối trực tiếp với lưới điện ba pha trong khi rotor được cấp nguồn qua hai bộ biến đổi nghịch lưu độc lập (RSC - Rotor Side Converter và GSC - Grid Side Converter) nối chung tuyến một chiều (DC-link). Thiết kế này cho phép bộ biến đổi công suất chỉ cần định mức bằng 20% đến 30% tổng công suất máy phát trong khi vẫn đáp ứng dải điều chỉnh tốc độ rộng $\pm 30%$ quanh tốc độ đồng bộ (hệ số trượt $s = \pm 0.3$).

Tuy nhiên, DFIG là một hệ thống cơ điện phi tuyến mạnh, đa biến, liên kết chéo phức tạp và vận hành liên tục trong môi trường khắc nghiệt ngoài khơi lẫn đất liền. Mặc dù các phương pháp điều khiển vector truyền thống và điều khiển hiện đại (Backstepping, Sliding Mode, Adaptive Control) đã đạt được nhiều bước tiến, thực tiễn vận hành cho thấy hệ thống WECS-DFIG vẫn đối mặt với các sự cố không thể lường trước. Điển hình, sự cố cách điện cuộn dây stator chiếm từ 30% đến 40% tổng số hư hỏng của máy điện cảm ứng. Quan trọng hơn, sự biến thiên tham số cơ - điện, suy giảm cơ học của hộp số, biến động gió cực đoan và nhiễu loạn lưới điện thường kích hoạt các dao động phi tuyến bất thường mang bản chất hỗn loạn (chaos). Hiện tượng hỗn loạn trong hệ truyền động điện và lưới điện đã được chứng minh là nguồn gốc gây ra các thảm họa rã lưới diện rộng, tiêu biểu như sự cố mất điện tại Hoa Kỳ (1966) và sự cố "9.26" tại Trung Quốc (2005). Đối với DFIG, các quỹ đạo hỗn loạn không tuần hoàn làm gia tăng mệt mỏi cơ học, tạo sóng hài dòng điện nguy hiểm, gây quá áp trên bus một chiều và làm suy giảm nghiêm trọng độ tin cậy vận hành.

Khoảng trống nghiên cứu (Research Gap) cốt lõi hiện nay là: phần lớn các nghiên cứu về hỗn loạn trước đây chỉ tập trung vào động cơ đồng bộ kích từ nam châm vĩnh cửu (PMSM) hoặc các mạch nghịch lưu độc lập, trong khi các công bố về phân tích hành vi hỗn loạn và giải pháp triệt tiêu hỗn loạn cho hệ thống DFIG còn rất khiêm tốn. Các nghiên cứu hiện hữu thường tiếp cận đơn lẻ từng chi tiết phần cứng mà chưa đưa ra được một khung toán học mô hình hóa xét gộp toàn phần (lumped modeling) cho toàn bộ hệ thống DFIG tích hợp cả RSC, GSC và DC-link khi có sự xuất hiện của các thành phần bất định sinh hỗn loạn.

Luận án tiến sĩ kỹ thuật của tác giả Cồ Như Văn, dưới sự hướng dẫn khoa học của GS. Nguyễn Phùng Quang và PGS. Nguyễn Thanh Hải tại Đại học Bách khoa Hà Nội (mã số chuyên ngành: 9520216), đã giải quyết triệt để khoảng trống này thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Những cơ chế và thông số nào kích hoạt trạng thái vận hành hỗn loạn trong hệ thống phát điện sức gió DFIG?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để xây dựng mô hình toán học tổng quát xét gộp toàn phần các hành vi bất định và hỗn loạn của DFIG trên hệ tọa độ tựa theo điện áp lưới?
  • Câu hỏi nghiên cứu 3 (RQ3): Có thể thiết kế cấu trúc quan sát nào để nhận dạng và ước lượng trực tuyến các biến trạng thái cùng các thành phần bất định sinh hỗn loạn trong điều kiện đo lường thực tế bị nhiễu?
  • Câu hỏi nghiên cứu 4 (RQ4): Giải pháp điều khiển nào có khả năng dập tắt nhanh chóng dao động hỗn loạn, khôi phục trạng thái làm việc ổn định và nâng cao độ tin cậy toàn hệ thống?

Hệ thống giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Sự gia tăng quá mức của hệ số khuếch đại mạch vòng tốc độ ($K_{P\omega}$) hoặc biến động đột ngột của mô-men cơ học đầu vào ($T_L$) sẽ đẩy các giá trị riêng của ma trận Jacobian qua trục ảo, dẫn đến phân nhánh Hopf và hình thành tập hút lạ (strange attractor).
  • Giả thuyết 2 (H2): Thuật toán lọc Kalman mở rộng (EKF) với cấu trúc vector trạng thái mở rộng có thể ước lượng đồng thời và chính xác các biến dòng điện rotor ($i_{rd}, i_{rq}$), tốc độ góc ($\omega$) và các thành phần bất định ($z_1, z_2$) đại diện cho nhiễu loạn hỗn loạn.
  • Giả thuyết 3 (H3): Bộ điều khiển phản hồi trạng thái và bộ điều khiển phản hồi trễ (Delay Feedback Control - DFC) theo nguyên lý Pyragas có thể triệt tiêu thành công số mũ Lyapunov dương, đưa hệ thống DFIG từ trạng thái hỗn loạn về quỹ đạo tuần hoàn ổn định mà chỉ cần tiêu tốn năng lượng điều khiển nhỏ.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Lý thuyết Động lực học Phi tuyến (Nonlinear Dynamics Theory), Lý thuyết Hỗn loạn Lorenz/Poincaré, Tiêu chuẩn Tính ổn định Lyapunov, Lý thuyết Ước lượng Tối ưu Kalman và Lý thuyết Điều khiển Vector tựa theo điện áp lưới. Đóng góp đột phá của công trình được lượng hóa thông qua việc chuyển hóa chỉ số Lyapunov cực đại (Maximum Lyapunov Exponent - MLE) từ giá trị dương ($\lambda_1 > 0$, biểu thị hỗn loạn) về giá trị âm ($\lambda_1 < 0$, biểu thị ổn định xác lập), duy trì độ lệch điện áp DC-link trong giới hạn an toàn dưới 5% khi xảy ra sự cố. Phạm vi nghiên cứu bao quát toàn bộ động học cơ - điện của máy phát DFIG công suất lớn trong dải hệ số trượt $s \in [-0.3, +0.3]$, được kiểm chứng toàn diện từ mô phỏng số trên Matlab/Simulink đến thực nghiệm thời gian thực (Real-time Hardware-in-the-Loop) trên thiết bị công nghiệp Typhoon HIL 402 phối hợp với card xử lý tín hiệu số Texas Instruments C2000 DSP.

Literature Review và Positioning

Nghiên cứu về điều khiển máy điện không đồng bộ nguồn kép và lý thuyết hỗn loạn hình thành từ ba dòng chảy học thuật chính trong hơn một thế kỷ qua:

┌────────────────────────────────────────────────────────────────────────┐
│                        CÁC DÒNG CHẢY HỌC THUẬT                         │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
    ┌───────────────────────────────┼───────────────────────────────┐
    ▼                               ▼                               ▼
┌───────────────────────┐   ┌───────────────────────┐   ┌───────────────────────┐
│     DÒNG CHẢY 1       │   │      DÒNG CHẢY 2      │   │      DÒNG CHẢY 3      │
│Điều khiển Vector DFIG │   │  Lý thuyết Hỗn loạn   │   │  Điều khiển Hỗn loạn  │
│  Blaschke & Leonhard  │   │  Poincaré, Lorenz,    │   │  OGY (1990), Pyragas  │
│ (1970s), Quang (2010) │   │     Lyapunov (1992)   │   │  (1992), Chen (1998)  │
└───────────────────────┘   └───────────────────────┘   └───────────────────────┘
  1. Dòng chảy điều khiển định hướng trường và điều khiển phi tuyến DFIG: Khởi nguồn từ các công trình kinh điển của F. Blaschke (1972) và W. Leonhard (1985) về điều khiển định hướng từ thông (Field Oriented Control - FOC), sau đó được phát triển thành phương pháp tựa theo điện áp lưới cho DFIG bởi các công trình của Nguyen Phung Quang (2008, 2010), Pena et al. (1996), và Tapia et al. (2003). Các nghiên cứu này thiết lập khả năng điều khiển phân lập độc lập công suất tác dụng ($P_s$) và công suất phản kháng ($Q_s$) thông qua hai thành phần dòng điện rotor $i_{rd}$ và $i_{rq}$. Về sau, các thuật toán phi tuyến như điều khiển trượt (Sliding Mode Control - SMC), điều khiển cuốn chiếu (Backstepping Control), và điều khiển thích nghi mờ đã được áp dụng để tăng tính kháng nhiễu đối với sự thay đổi điện trở và điện cảm stator/rotor.
  2. Dòng chảy phân tích động lực học phi tuyến và hỗn loạn: Bắt đầu từ khám phá nền tảng của Henri Poincaré (1890) về tính nhạy cảm với điều kiện ban đầu trong bài toán ba vật thể, tiếp nối bởi Edward Lorenz (1963) khi mô hình hóa đối lưu khí quyển dẫn đến tập hút lạ Lorenz. Trong kỹ thuật điện tử công suất và truyền động điện, các nghiên cứu của Hamill, Deane, Jefferies (1992) và sau đó là Guanrong Chen, Z.G. Li (1998, 2002) trên động cơ PMSM đã chứng minh rằng các phần tử đóng cắt bán dẫn và liên kết phi tuyến giữa từ thông - mô men là môi trường thuận lợi làm bùng phát hiện tượng hỗn loạn.
  3. Dòng chảy điều khiển và triệt tiêu hỗn loạn: Được khai mở bởi phương pháp đột phá của Ott, Grebogi, và Yorke (OGY, 1990) dựa trên việc ổn định hóa các quỹ đạo tuần hoàn không ổn định (Unstable Periodic Orbits - UPOs) bằng các tác động xung nhỏ, và phương pháp điều khiển phản hồi trễ liên tục của K. Pyragas (1992).

Trong bức tranh tổng quan đó, một cuộc tranh luận học thuật sâu sắc đã nổ ra giữa hai trường phái:

  • Trường phái kỹ thuật điện truyền thống: Cho rằng các hiện tượng dao động bất thường tần số thấp (Sub-synchronous Resonance - SSR) hay biến dạng dòng điện trong hệ thống WECS-DFIG thuần túy là hệ quả của nhiễu ngẫu nhiên ngoài lưới, cộng hưởng tuyến tính giữa đường dây mang tụ bù và cảm kháng máy phát, hoặc do các khâu hiệu chỉnh PID bị sai lệch thông số. Giải pháp truyền thống do đó chỉ tập trung vào việc thiết kế các bộ lọc thụ động/tích cực hoặc bổ sung khâu bù quán tính ảo.
  • Trường phái động lực học phi tuyến: Minh chứng rằng nhiều dao động biên độ lớn không tuần hoàn thực chất là các chuyển động tất định (deterministic chaos) phát sinh từ các điểm phân nhánh (Hopf bifurcation hoặc period-doubling bifurcation) khi hệ thống chịu tải nặng, suy giảm cách điện cuộn dây, hoặc tham số điều khiển vượt ngưỡng ổn định. Các giải pháp lọc tuyến tính truyền thống hoàn toàn bất lực trước tập hút lạ hỗn loạn và thậm chí có thể khiến dao động phân kỳ nhanh hơn.
┌─────────────────────────────────────────────────────────────────────────┐
│              SO SÁNH VỊ TRÍ HỌC THUẬT VÀ ĐÓNG GÓP TIÊN PHONG            │
├──────────────────────┬─────────────────────────┬────────────────────────┤
│ Nghiên cứu quốc tế   │ Phạm vi & Hạn chế       │ Vị trí định vị của     │
│ điển hình            │                         │ Luận án Cồ Như Văn     │
├──────────────────────┼─────────────────────────┼────────────────────────┤
│ 1. Chen & Dong (1998)│ Khảo sát hỗn loạn trong │ Mở rộng từ hệ PMSM sang│
│    "Controlling      │ động cơ PMSM bậc 3;     │ hệ DFIG bậc cao phức   │
│    Chaos in PMSM"    │ không xét dòng năng     │ tạp có trao đổi công   │
│                      │ lượng 2 chiều qua lưới. │ suất hai chiều stator/ │
│                      │                         │ rotor và bus DC.       │
├──────────────────────┼─────────────────────────┼────────────────────────┤
│ 2. Ott, Grebogi,     │ Phương pháp OGY đòi hỏi │ Đề xuất mô hình xét gộp│
│    Yorke (1990)      │ tính toán ma trận Po-   │ toàn phần kết hợp EKF- │
│    "Controlling      │ incaré gián đoạn thời   │ DFC, ước lượng trực    │
│    Chaos"            │ gian thực; khó khả thi  │ tuyến và dập tắt hỗn   │
│                      │ trên vi xử lý WECS.     │ loạn liên tục trên DSP.│
└──────────────────────┴─────────────────────────┴────────────────────────┘

Luận án của tác giả Cồ Như Văn đã định vị chính xác tại điểm giao thoa giữa điều khiển truyền động nâng cao và lý thuyết hỗn loạn. Bằng việc xây dựng mô hình không gian trạng thái tích hợp toàn diện và áp dụng thành công bộ quan sát EKF mở rộng để bù thành phần bất định, công trình đã lấp đầy khoảng trống tồn tại nhiều năm giữa lý thuyết hỗn loạn trừu tượng và thực tiễn kỹ thuật điều khiển hệ thống điện gió DFIG công nghiệp.

Đó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à làm phong phú thêm Lý thuyết Động lực học Phi tuyến trong Kỹ thuật Điện và Lý thuyết Điều khiển Hiện đại thông qua các đóng góp chuyên sâu:

  1. Mở rộng lý thuyết hệ thống động lực Lorenz cho máy phát không đồng bộ nguồn kép: Công trình chứng minh rằng mô hình toán học của DFIG trong hệ tọa độ tựa theo điện áp lưới ($dq$) khi chịu tương tác chéo giữa các thành phần dòng $i_{rd}, i_{rq}$ và sai lệch tốc độ góc $\omega$ có thể quy đổi tương đương về cấu trúc của hệ động lực hỗn loạn ba chiều thỏa mãn điều kiện cần của hệ Lorenz: $$\begin{cases} \dot{x}_1 = -\sigma x_1 + \sigma y_1 \ \dot{x}_2 = \rho x_1 - x_2 - x_1 x_3 \ \dot{x}_3 = -\beta x_3 + x_1 x_2 \end{cases}$$ Trong đó, sự kết hợp phi tuyến bậc cao giữa dòng điện và từ thông đóng vai trò như các số hạng tích $x_1 x_2$ và $x_1 x_3$, tạo ra cơ chế phân kỳ cục bộ nhưng bị giam hãm trong không gian pha giới nội.

  2. Xác lập mô hình không gian trạng thái tổng quát có xét đến thành phần hỗn loạn: Hệ thống được biểu diễn dưới dạng mô hình liên tục trong không gian trạng thái: $$\frac{dx}{dt} = A x + B_s u_s + B_r u_r + D z(t)$$ Trong đó $x = [i_{rd}, i_{rq}, \psi_{sd}', \psi_{sq}']^T$ là vector biến trạng thái; $u_s, u_r$ là vector điện áp điều khiển stator và rotor; $z(t) = [z_1(t), z_2(t)]^T$ là vector đầu vào chưa xác định (unknown input) đại diện cho toàn bộ các thành phần phi tuyến hỗn loạn, biến thiên tham số và nhiễu loạn mô-men ngoài tác động lên hệ thống.

┌────────────────────────────────────────────────────────────────────────┐
│             KHUNG KHÁI NIỆM MÔ HÌNH HÓA XÉT GỘP TOÀN PHẦN             │
├────────────────────────────────────────────────────────────────────────┤
│                                                                        │
│   ┌─────────────────────┐               ┌──────────────────────────┐   │
│   │ Biến đổi tham số:   │               │ Nhiễu loạn cơ & sự cố:   │   │
│   │ ΔRs, ΔRr, ΔLs, ΔLr  │               │ ΔTL, lệch hộp số, sự cố  │   │
│   │ do nhiệt độ/tuổi thọ│               │ cách điện stator 30%-40% │   │
│   └──────────┬──────────┘               └────────────┬─────────────┘   │
│              │                                       │                 │
│              └───────────────────┬───────────────────┘                 │
│                                  ▼                                     │
│                 ┌─────────────────────────────────┐                    │
│                 │ VECTOR THÀNH PHẦN BẤT ĐỊNH XÉT  │                    │
│                 │       GỘP: z(t) = [z1, z2]T     │                    │
│                 └────────────────┬────────────────┘                    │
│                                  │                                     │
│                                  ▼                                     │
│                 ┌─────────────────────────────────┐                    │
│                 │   MÔ HÌNH KHÔNG GIAN TRẠNG THÁI │                    │
│                 │      DFIG ĐẶC TRƯNG HỖN LOẠN    │                    │
│                 │    dx/dt = Ax + Bsus + Brur + Dz│                    │
│                 └─────────────────────────────────┘                    │
│                                                                        │
└────────────────────────────────────────────────────────────────────────┘
  1. Mệnh đề toán học về ngưỡng phân nhánh (Bifurcation Propositions):
    • Mệnh đề 1 (P1): Tồn tại một ngưỡng tới hạn của hệ số tỷ lệ mạch vòng vận tốc $K_{P\omega}^*$, tại đó các cặp giá trị riêng của hệ kín vượt qua trục ảo từ nửa trái sang nửa phải mặt phẳng phức, kích hoạt phân nhánh Hopf làm triệt tiêu tính ổn định tiệm cận của điểm cân bằng tĩnh.
    • Mệnh đề 2 (P2): Thành phần bất định gộp $z(t)$ được giới nội trong một compact set $\Omega_z \subset \mathbb{R}^2$, đảm bảo sự tồn tại của tập hút lạ có số chiều fractal xác định.
    • Mệnh đề 3 (P3): Thuật toán phản hồi trễ $u_c(t) = K [x(t) - x(t - \tau)]$ với ma trận khuếch đại $K$ và thời gian trễ $\tau$ thỏa mãn miền giới hạn $U_0(K, \tau)$ sẽ đảm bảo chuyển dịch số mũ Lyapunov cực đại $\lambda_1$ từ miền thực dương về miền thực âm.

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

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết hiện đại: (1) Lý thuyết Điều khiển Vector tựa theo điện áp lưới; (2) Lý thuyết Ước lượng Trạng thái Tối ưu Phi tuyến Kalman (EKF); và (3) Lý thuyết Ổn định hóa Phản hồi Trễ Pyragas.

Sự đột phá trong phương pháp phân tích nằm ở Tư duy Mô hình hóa Xét gộp Toàn phần (Lumped Parameter Aggregation): Thay vì cố gắng phân tách và mô hình hóa chi tiết từng dạng lỗi vật lý riêng rẽ (như nứt răng hộp số, ngắn mạch vài vòng dây stator, hay trôi dạt tham số điện cảm tản $\sigma L_r$), luận án đã khéo léo ánh xạ toàn bộ các hiệu ứng phi tuyến gây bất ổn định này thành các tín hiệu ngoại suy tác động trực tiếp vào phương trình trạng thái của dòng điện rotor $i_{rd}$ và $i_{rq}$.

Điều kiện biên xác lập của khung phân tích:

  • Máy phát DFIG hoạt động trong dải tốc độ từ dưới đồng bộ đến trên đồng bộ ($s \in [-0.3, +0.3]$).
  • Điện áp lưới $u_s$ được đo và khóa pha chính xác thông qua vòng khóa pha (PLL) chất lượng cao, với trục $d$ định hướng hoàn toàn theo vector điện áp lưới ($u_{sd} = U_s, u_{sq} = 0$).
  • Bỏ qua hiện tượng bão hòa mạch từ sâu trong các phân tích trạng thái xác lập cơ bản.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học Thực chứng (Positivism) kết hợp với chủ nghĩa Hiện thực Phê phán (Critical Realism), sử dụng phương pháp diễn dịch toán học từ các định luật cơ bản của kỹ thuật điện (Maxwell, Newton) và kiểm chứng thực nghiệm bằng phương pháp nghiên cứu hỗn hợp (Mixed-Method Quantitative Experimental Design).

Thiết kế nghiên cứu đa tầng (Multi-level Architecture) được triển khai qua 3 cấp độ:

  • Cấp độ 1 (Thiết bị đơn lẻ - Machine Level): Mô hình hóa toán học chính xác máy điện DFIG với các tham số điện trở stator $R_s$, điện trở rotor $R_r$, hỗ cảm $L_m$, điện cảm tản stator $L_{s\sigma}$, điện cảm tản rotor $L_{r\sigma}$, mô-men quán tính $J$, hệ số ma sát nhớt $B$ và số đôi cực $n_p$.
  • Cấp độ 2 (Hệ thống biến đổi công suất - Converter & DC-link Level): Mô hình hóa chi tiết hai bộ nghịch lưu RSC và GSC với bộ lọc lưới $R_D, L_D, R_F, C_F$ và tụ điện bus một chiều $C_{DC}$, thiết lập các phương trình cân bằng năng lượng: $$\frac{3}{2} u_{Nd} i_{Nd} = u_{DC} i_{DC} = u_{DC} C_{DC} \frac{du_{DC}}{dt}$$
  • Cấp độ 3 (Cấp độ hệ thống điều khiển & quan sát - Control & Observation Level): Tích hợp hai bộ quan sát EKF-M1, EKF-M2 cùng hai cấu trúc điều khiển triệt tiêu hỗn loạn (Phản hồi trạng thái và Phản hồi trễ).
┌────────────────────────────────────────────────────────────────────────┐
│             KIẾN TRÚC THIẾT KẾ NGHIÊN CỨU THỰC NGHIỆM ĐA TẦNG          │
├────────────────────────────────────────────────────────────────────────┤
│                                                                        │
│   ┌────────────────────────────────────────────────────────────────┐   │
│   │                      MÔ PHỎNG SỐ OFFLINE                       │   │
│   │             (Matlab/Simulink R2020b / SimPowerSystems)         │   │
│   └───────────────────────────────┬────────────────────────────────┘   │
│                                   │ Kiểm chứng đối sánh               │
│                                   ▼                                    │
│   ┌────────────────────────────────────────────────────────────────┐   │
│   │               MÔ PHỎNG THỜI GIAN THỰC HIL 402                  │   │
│   ├───────────────────────────────┬────────────────────────────────┤   │
│   │       PHẦN CỨNG HIL           │         CARD ĐIỀU KHIỂN        │   │
│   │      (Typhoon HIL 402)        │        (TI C2000 DSP)          │   │
│   │  - Mô phỏng thực tế DFIG      │   - Chạy thuật toán EKF-M1/M2  │   │
│   │  - Mô phỏng nghịch lưu RSC/GSC│   - Chạy bộ điều khiển DFC     │   │
│   │  - Mô phỏng lưới & gió thật   │   - Xuất xung PWM & nhận tín   │   │
│   │  - Bước trích mẫu: 1 MHz      │     hiệu phản hồi dòng, áp     │   │
│   └───────────────────────────────┴────────────────────────────────┘   │
│                                                                        │
└────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình nhận dạng và điều khiển hỗn loạn trải qua 4 bước thực nghiệm nghiêm ngặt:

  1. Bộ công cụ chẩn đoán và nhận dạng đa chiều (Triangulation Diagnostics):

    • Phân tích miền thời gian (Time-series analysis): Quan sát dạng sóng không tuần hoàn, dao động tự kích biên độ lớn của $i_{rd}, i_{rq}, \omega$.
    • Mặt phẳng pha và không gian pha 3D (Phase portraits): Khảo sát quỹ đạo pha giữa $i_{rd} - \omega$, $i_{rd} - i_{rq} - \omega$ để phát hiện cấu trúc cuộn xoắn của tập hút lạ.
    • Biểu đồ phân nhánh (Bifurcation diagram): Khảo sát bước nhảy trạng thái khi quét thông số mô-men $T_L$ và hệ số khuếch đại $K_{P\omega}$.
    • Lát cắt Poincaré (Poincaré map): Lấy giao điểm của quỹ đạo trạng thái với siêu phẳng chọn trước trong không gian 3 chiều để xác định cấu trúc hình học fractal rời rạc.
    • Phân tích mật độ phổ Fourier (FFT Continuous Spectrum): Kiểm tra tính liên tục của phổ tần số (phổ dải rộng đặc trưng của chaos) thay vì các gai phổ vạch rời rạc của dao động điều hòa/chu kỳ.
    • Thuật toán định lượng số mũ Lyapunov cực đại (MLE): $$\lambda_1 = \lim_{n \to \infty} \frac{1}{n} \sum_{k=0}^{n-1} \ln |f'(x_k)|$$ Khi $\lambda_1 > 0$, hệ thống được khẳng định chắc chắn 100% đang rơi vào chế độ hỗn loạn tất định.

  2. Xây dựng hai cấu trúc bộ quan sát Kalman mở rộng (EKF Observers):

    • Phương pháp 1 (EKF-M1): Bổ sung biến trạng thái mở rộng $x_a = [i_{rd}, i_{rq}, \omega, z_1, z_2]^T$, tính toán ma trận Jacobi liên tục và cập nhật ma trận hiệp phương sai sai số $P(k)$ trực tuyến nhằm bóc tách đồng thời các biến trạng thái và 2 thành phần bất định $z_1, z_2$.
    • Phương pháp 2 (EKF-M2): Tách biệt khâu lọc trạng thái cơ bản và khâu ước lượng nhiễu loạn bậc cao, tối ưu hóa tốc độ tính toán cho vi xử lý thời gian thực.

Data và phân tích

Bộ thông số chuẩn hóa của máy phát DFIG công nghiệp được thiết lập trong môi trường mô phỏng và thực nghiệm phần cứng như sau:

  • Công suất định mức: $P_n = 1.5\text{ MW} - 2.0\text{ MW}$
  • Điện áp stator định mức: $V_s = 690\text{ V}$, tần số lưới $f = 50\text{ Hz}$
  • Điện trở cuộn dây stator: $R_s = 0.00706\ \Omega$; Điện cảm stator: $L_s = 0.00307\text{ H}$
  • Điện trở cuộn dây rotor: $R_r = 0.005\ \Omega$; Điện cảm rotor: $L_r = 0.00306\text{ H}$
  • Hỗ cảm stator - rotor: $L_m = 0.0029\text{ H}$
  • Điện cảm lọc phía lưới: $L_D = 0.0003\text{ H}$, Điện trở lọc $R_D = 0.005\ \Omega$
  • Tụ điện liên kết DC: $C_{DC} = 0.01\text{ F}$, Điện áp định mức bus một chiều $U_{DC} = 1200\text{ V}$
  • Mô-men quán tính: $J = 0.5\text{ kg}\cdot\text{m}^2$; Hệ số giảm chấn $B = 0.005\text{ N}\cdot\text{m}\cdot\text{s/rad}$; Số đôi cực $n_p = 2$

Kiểm tra độ bền vững (Robustness Checks) được thực hiện thông qua việc kích hoạt đồng thời:

  • Bước nhảy điện trở stator $\Delta R_s = +50%$ (mô phỏng sự suy giảm cách điện và tăng nhiệt độ cuộn dây).
  • Sụt áp lưới stator $\Delta V_s = -20%$.
  • Tăng đột biến mô-men cơ học đầu vào $T_L$ từ giá trị định mức lên dải quá tải kích hoạt hỗn loạn.

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

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

  1. Xác định chính xác các điểm rẽ nhánh kích hoạt hỗn loạn: Luận án chỉ ra rằng khi tăng thông số khuếch đại vận tốc $K_{P\omega}$ của bộ điều khiển PI truyền thống hoặc khi mô-men xoắn cơ học $T_L$ vượt ngưỡng giới hạn cho phép, hệ thống DFIG sẽ mất tính ổn định tuần hoàn thông qua chuỗi phân nhánh nhân đôi chu kỳ (period-doubling cascade) và rơi vào trạng thái hỗn loạn hoàn toàn.
  2. Khám phá định lượng thông qua số mũ Lyapunov: Ở chế độ vận hành bình thường, số mũ Lyapunov lớn nhất đạt giá trị âm ($\lambda_1 = -0.041 < 0$). Khi xuất hiện các thành phần bất định và sự cố kích thích, $\lambda_1$ nhanh chóng chuyển sang giá trị dương ($\lambda_1 = +0.385 > 0$), chứng minh hệ thống đã hình thành một tập hút lạ (strange attractor) giam hãm không gian trạng thái.
  3. Hiệu năng bóc tách bất định vượt trội của EKF-M1 và EKF-M2: Cả hai thuật toán EKF đều ước lượng bám sát hoàn hảo các biến dòng điện rotor thực ($i_{rd}, i_{rq}$) và tốc độ góc $\omega$ với sai số dưới $1.2%$. Quan trọng nhất, các thành phần chưa xác định $z_1, z_2$ được tái tạo với độ chính xác cao ngay cả khi có sự biến thiên lớn của điện trở $R_s$ và sụt giảm điện áp lưới $V_s$.
  4. Hiệu quả triệt tiêu hỗn loạn tức thời của bộ điều khiển phản hồi trễ (DFC): Với việc lựa chọn tối ưu ma trận khuếch đại $K$ (ví dụ $K = 15$) và thời gian trễ $\tau = 0.02\text{ s}$ (tương ứng chu kỳ lưới $50\text{ Hz}$), bộ điều khiển DFC đã dập tắt hoàn toàn các dao động hỗn loạn của dòng $i_{rd}, i_{rq}$, ổn định điện áp bus $U_{DC}$ về giá trị đặt $1200\text{ V}$ trong thời gian quá độ chỉ dưới $0.15\text{ s}$ mà không làm biến dạng sóng cơ bản.
┌────────────────────────────────────────────────────────────────────────┐
│                   SO SÁNH ĐỘNG THÁI TRƯỚC VÀ SAU ĐIỀU KHIỂN             │
├──────────────────────┬────────────────────────┬────────────────────────┤
│ Đại lượng khảo sát   │ Trạng thái hỗn loạn    │ Sau khi kích hoạt DFC  │
│                      │ (Chưa điều khiển)      │ (Đã triệt tiêu chaos)  │
├──────────────────────┼────────────────────────┼────────────────────────┤
│ Số mũ Lyapunov (MLE) │ Dương: λ1 = +0.385 > 0 │ Âm: λ1 = -0.052 < 0    │
│ Quỹ đạo pha không    │ Cuộn xoắn phân kỳ,     │ Quỹ đạo khép kín trơn  │
│ gian (ird - irq - ω) │ tập hút lạ đa chiều    │ tru, ổn định tuần hoàn │
│ Điện áp bus DC (UDC) │ Dao động méo dạng vô   │ Ổn định xác lập tại    │
│                      │ định 950V - 1450V      │ 1200V (sai số < 1%)    │
│ Dòng điện phía lưới  │ Méo hài nặng nề, THD   │ Sin chuẩn 50 Hz,       │
│                      │ vượt ngưỡng an toàn    │ THD < 2.5%             │
└──────────────────────┴────────────────────────┴────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Thiết lập mô hình toán học giải tích chuẩn mực cho phân tích động lực học hỗn loạn trong các hệ truyền động xoay chiều nguồn kép phức tạp, tạo tiền đề mở rộng lý thuyết điều khiển hỗn loạn cho các cấu trúc năng lượng tái tạo mới.
  • Về mặt phương pháp luận: Cung cấp quy trình phối hợp hoàn chỉnh giữa bộ quan sát ước lượng tối ưu EKF với kỹ thuật điều khiển phản hồi trễ Pyragas, chứng minh khả năng kiểm soát hệ phi tuyến bậc cao thông qua can thiệp năng lượng tối thiểu.
  • Về mặt ứng dụng thực tiễn: Toàn bộ thuật toán được biên dịch và chạy trực tiếp trên chip xử lý tín hiệu số TI C2000 DSP kết nối với thiết bị Typhoon HIL 402. Điều này chứng minh thuật toán hoàn toàn sẵn sàng để đóng gói thành các module phần mềm/firmware tích hợp vào các tủ biến tần DFIG công nghiệp hiện hành mà không cần thay đổi phần cứng máy phát.
  • Về mặt chính sách và an ninh năng lượng: Cung cấp công cụ dự báo sớm nguy cơ mất ổn định hỗn loạn cho các trung tâm điều độ hệ thống điện, hỗ trợ duy trì tiêu chuẩn vận hành lưới điện (Grid Codes) khi tỷ trọng thâm nhập của điện gió vượt mức 30% - 50%.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu (Boundary Conditions):

  1. Phạm vi mô hình lưới điện: Nghiên cứu hiện tại tập trung chủ yếu vào mô hình một máy phát kết nối với thanh cái vô cùng (Single-Machine Infinite Bus - SMIB), chưa mô phỏng chi tiết tương tác hỗn loạn đa máy giữa hàng trăm tuabin gió trong một trang trại điện gió quy mô lớn (Multi-turbine Wind Farm Interaction).
  2. Bản chất của mô hình xét gộp: Phương pháp xét gộp toàn phần gộp chung tất cả các nguồn kích động phi tuyến vào vector bất định $z(t)$, do đó không thực hiện định lượng riêng rẽ phần trăm đóng góp của từng lỗi cơ học cụ thể (như mẻ răng hộp số hay rơ lỏng ổ bi).
  3. Môi trường thực nghiệm: Nghiên cứu mới dừng lại ở mức mô phỏng thời gian thực Hardware-in-the-Loop trên thiết bị công nghiệp Typhoon HIL 402 và DSP C2000, chưa tiến hành thử nghiệm trên tuabin gió thực địa 2.0 MW do rào cản chi phí và an toàn vận hành thử nghiệm thực địa.

Các hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Hướng 1: Mở rộng mô hình toán học để khảo sát hiện tượng đồng bộ hóa hỗn loạn (Chaos Synchronization) và lan truyền hỗn loạn giữa các cụm trang trại điện gió ngoài khơi (Offshore Wind Farms) kết nối qua lưới truyền tải một chiều cao áp (HVDC).
  • Hướng 2: Tích hợp các thuật toán học sâu (Deep Reinforcement Learning) và mạng nơ-ron thích nghi vào khâu cập nhật ma trận hiệp phương sai của EKF nhằm tự động tinh chỉnh thông số bộ lọc khi các đặc tính cơ điện trôi dạt phi tuyến mạnh theo thời gian.
  • Hướng 3: Nghiên cứu ứng dụng lý thuyết vi phân bậc phân số (Fractional-order Calculus) trong thiết kế bộ điều khiển phản hồi trễ phân số (Fractional DFC) để tối ưu hóa hơn nữa thời gian dập tắt hỗn loạn.
  • Hướng 4: Triển khai thử nghiệm thực tế (Field Pilot Testing) trên các hệ thống tuabin gió thương mại quy mô nhỏ và vừa tại Việt Nam.

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

┌────────────────────────────────────────────────────────────────────────┐
│                     HỆ SINH THÁI TÁC ĐỘNG ĐA CHIỀU                     │
├────────────────────────────────────────────────────────────────────────┤
│                                                                        │
│  ┌─────────────────────────┐             ┌──────────────────────────┐  │
│  │   TÁC ĐỘNG HỌC THUẬT    │             │    TÁC ĐỘNG CÔNG NGHIỆP  │  │
│  │ Mở ra hướng nghiên cứu  │             │ Nâng cao tuổi thọ WECS,  │  │
│  │ mới về chaos trong DFIG,│             │ giảm chi phí O&M từ      │  │
│  │ dự kiến thu hút trích   │             │ 15% - 20% thông qua bảo  │  │
│  │ dẫn quốc tế cao (Scopus)│             │ trì dự đoán sớm.         │  │
│  └────────────┬────────────┘             └─────────────┬────────────┘  │
│               │                                        │               │
│               └───────────────────┬────────────────────┘               │
│                                   ▼                                    │
│                 ┌───────────────────────────────────┐                  │
│                 │   AN NINH NĂNG LƯỢNG & XÃ HỘI     │                  │
│                 │ Triệt tiêu nguy cơ rã lưới diện   │                  │
│                 │ rộng, ổn định tần số quốc gia,    │                  │
│                 │ thúc đẩy chuyển dịch Net Zero 2050│                  │
│                 └───────────────────────────────────┘                  │
│                                                                        │
└────────────────────────────────────────────────────────────────────────┘
  • Tác động học thuật: Đặt nền móng vững chắc cho hướng nghiên cứu điều khiển hỗn loạn ứng dụng trong thiết bị năng lượng tái tạo tại Việt Nam, mở rộng các diễn đàn học thuật chuyên sâu về phân tích dao động phi tuyến trong các tạp chí ISI/Scopus đầu ngành.
  • Chuyển đổi công nghiệp: Giúp các đơn vị chế tạo và vận hành trang trại điện gió chuyển đổi phương thức bảo trì từ "bảo trì định kỳ theo kế hoạch" sang "bảo trì dự đoán trên cơ sở chẩn đoán hỗn loạn sớm", giảm thiểu thời gian dừng máy ngoài dự kiến (downtime), tiết kiệm ước tính 15% - 20% chi phí vận hành và bảo dưỡng (O&M).
  • An ninh lưới điện quốc gia: Giảm thiểu rủi ro kích động dao động liên vùng và hiện tượng ngắt kết nối dây chuyền của các nhà máy điện gió khi xảy ra sự cố, bảo vệ an toàn cho hệ thống điện truyền tải 500 kV quốc gia.

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

  • Nghiên cứu sinh và Giảng viên ngành Điều khiển & Tự động hóa: Tiếp cận được một khung phương pháp luận toàn diện từ xây dựng mô hình toán, lập trình bộ lọc EKF đến kỹ thuật mô phỏng phần cứng HIL bán thực.
  • Kỹ sư R&D tại các tập đoàn sản xuất thiết bị điện gió (OEMs): Sở hữu giải pháp thuật toán điều khiển đã được kiểm chứng trên nền tảng DSP TI C2000, có thể thương mại hóa trực tiếp thành các tính năng cao cấp (Advanced Chaos-suppression Firmware).
  • Các đơn vị vận hành trang trại gió và Trung tâm Điều độ Hệ thống điện (A0/NSMO): Nhận được các công cụ phân tích tin cậy để giám sát chất lượng điện năng, nhận biết sớm các bất ổn định động học phi tuyến tiềm ẩn.

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? Trả lời: Đóng góp độc đáo nhất là việc mở rộng thành công Lý thuyết Hệ thống Động lực Phi tuyến Lorenz vào cấu trúc máy phát không đồng bộ nguồn kép DFIG và đề xuất nguyên lý "Mô hình hóa xét gộp toàn phần". Luận án đã chứng minh rằng các tương tác chéo phi tuyến dòng - áp - tốc độ trong hệ quy chiếu tựa theo điện áp lưới khi mất ổn định tham số sẽ tạo ra tập hút lạ ba chiều thỏa mãn các đặc trưng hỗn loạn kinh điển, từ đó cho phép gom các bất định phức tạp thành vector đầu vào chưa xác định $z(t)$ để giải quyết bằng bộ điều khiển hiện đại.

2. Đột phá phương pháp luận của luận án so với các nghiên cứu trước đây là gì? Trả lời: So với phương pháp OGY truyền thống (vốn đòi hỏi tính toán gián đoạn phức tạp trên mặt phẳng Poincaré khó áp dụng thực tế) hay các nghiên cứu trên động cơ PMSM của Chen & Dong (không xử lý được bài toán nối lưới hai chiều), luận án đã kết hợp đột phá giữa Bộ quan sát EKF đa biến ước lượng đồng thời trạng thái và thành phần bất định với Bộ điều khiển phản hồi trễ liên tục DFC. Cấu trúc này cho phép nhận dạng và dập tắt hỗn loạn hoàn toàn trong thời gian thực mà không làm tăng bậc của hệ kín và không cần đo lường các tín hiệu nhiễu loạn bên ngoài.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là hiện tượng hỗn loạn có thể xuất hiện ngay cả khi các tham số điện cơ của máy phát hoàn toàn nằm trong dải danh định, chỉ do sự gia tăng quá mức của hệ số khuếch đại mạch vòng tốc độ $K_{P\omega}$ trong bộ điều khiển PI truyền thống. Dữ liệu thực nghiệm HIL 402 chứng minh rằng khi $K_{P\omega}$ vượt ngưỡng, số mũ Lyapunov $\lambda_1$ lập tức chuyển từ $-0.041$ sang $+0.385$, làm méo dạng nghiêm trọng điện áp bus DC lên tới $1450\text{ V}$ và dòng điện rotor bị biến dạng phi tuần hoàn, điều mà các công cụ phân tích tuyến tính Bode/Nyquist thông thường không thể dự đoán được.

4. Quy trình sao chép và tái lập thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Có. Toàn bộ thông số chi tiết của máy phát DFIG ($R_s, R_r, L_s, L_r, L_m$), thông số bộ lọc lưới ($R_D, L_D, C_F$), các ma trận hiệp phương sai của bộ lọc Kalman ($Q, R$), cấu hình các kênh vào/ra trên thiết bị Typhoon HIL 402 và cấu trúc lập trình trên card DSP TMS320F28335/C2000 đều được trình bày minh bạch trong luận án, đảm bảo khả năng tái lập và kiểm chứng 100% kết quả nghiên cứu.

5. Lộ trình phát triển nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào? Trả lời: Lộ trình 10 năm được định hình qua 3 giai đoạn: (1) 2024–2026: Tích hợp trí tuệ nhân tạo (AI/Deep Learning) để tự thích nghi thông số EKF và DFC; (2) 2027–2030: Mở rộng mô hình kiểm soát hỗn loạn cho cụm trang trại điện gió ngoài khơi kết nối lưới điện HVDC-VSC; (3) 2031–2034: Chuẩn hóa bộ tiêu chuẩn kỹ thuật về chẩn đoán hỗn loạn và thương mại hóa toàn diện giải pháp nhúng trên các dòng tuabin gió thông minh thế hệ mới.

Kết luận

  1. Luận án đã chứng minh một cách khoa học và hệ thống về sự tồn tại tất yếu của hiện tượng hỗn loạn trong hệ thống phát điện sức gió sử dụng máy phát không đồng bộ nguồn kép (DFIG), khẳng định hỗn loạn là nguyên nhân gốc rễ dẫn đến nhiều sự cố suy giảm độ tin cậy và mất ổn định lưới điện.
  2. Đề xuất thành công mô hình toán học không gian trạng thái của DFIG trên hệ tọa độ tựa theo điện áp lưới theo phương pháp xét gộp toàn phần, phản ánh chính xác tác động của các thành phần bất định sinh hỗn loạn.
  3. Thiết kế thành công 02 cấu trúc bộ quan sát trạng thái tối ưu dựa trên thuật toán lọc Kalman mở rộng (EKF-M1 và EKF-M2), cho phép bóc tách và ước lượng trực tuyến chính xác các biến trạng thái cơ - điện cùng vector thành phần bất định $z(t)$.
  4. Xây dựng và thực thi thành công 02 giải pháp điều khiển triệt tiêu hỗn loạn dựa trên nguyên lý Phản hồi trạng thái và Bộ điều khiển phản hồi trễ (DFC), chuyển dịch thành công số mũ Lyapunov cực đại từ dương sang âm và đưa hệ thống về chế độ làm việc ổn định xác lập.
  5. Kiểm chứng toàn diện tính đúng đắn và hiệu năng của các thuật toán thông qua mô phỏng số offline trên Matlab/Simulink và thực nghiệm thời gian thực Hardware-in-the-Loop trên thiết bị công nghiệp Typhoon HIL 402 kết nối vi xử lý Texas Instruments C2000 DSP.
  6. Mở ra một phân ngành nghiên cứu mới về điều khiển trên cơ sở hỗn loạn trong các hệ thống năng lượng tái tạo, tạo tiền đề vững chắc cho việc nâng cao độ an toàn, tính kinh tế và hiệu quả vận hành của ngành công nghiệp điện gió toàn cầu.