Tổng quan về luận án
Sự gia tăng tỷ trọng nguồn năng lượng tái tạo trong cơ cấu nguồn điện toàn cầu đang đặt ra những thách thức kỹ thuật nghiêm trọng đối với tính ổn định và chất lượng điện năng của hệ thống điện phân phối. Trong bức tranh chuyển dịch năng lượng thế giới với tổng công suất điện gió tích lũy đạt 744 GW (đáp ứng 7% nhu cầu điện năng toàn cầu, dẫn đầu bởi Trung Quốc với 289 GW và Hoa Kỳ với 122 GW), công nghệ máy phát điện không đồng bộ nguồn kép (Doubly Fed Induction Generator - DFIG) tiếp tục chiếm vị thế thống trị khi nắm giữ gần 60% thị phần lắp đặt toàn cầu (GE chiếm 20,7%, Vestas chiếm 19%, Gamesa chiếm 19%). Tại Việt Nam, với tiềm năng kỹ thuật lớn cùng khung chính sách khuyến khích phát triển (Quyết định số 37/2011/QĐ-TTg và Quyết định số 39/2018/QĐ-TTg), hàng chục nghìn MW điện gió đã được bổ sung quy hoạch tại các khu vực trọng điểm như Nam Trung Bộ (4.193 MW), Tây Nam Bộ (25.000 MW).
Tuy nhiên, tính bất định của vận tốc gió với biên độ dao động công suất tác dụng có thể lên tới 60% trong chu kỳ 3 phút tại điểm kết nối chung (Point of Common Coupling - PCC) dẫn tới sụt giảm chất lượng điện áp nghiêm trọng, gây nhấp nháy điện áp (flicker), biến dạng sóng hài và đe dọa giới hạn ổn định điện áp tĩnh. Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Công Cường (chuyên ngành Kỹ thuật Điện, mã số: 9520201, bảo vệ tại Trường Đại học Mỏ - Địa chất Hà Nội dưới sự hướng dẫn của PGS.TS. Nguyễn Anh Nghĩa và PGS.TS. Trịnh Trọng Chưởng) mang tên "Nghiên cứu giải pháp nâng cao chất lượng điện năng của lưới điện phân phối có kết nối nguồn điện gió sử dụng máy điện không đồng bộ" đã trực tiếp giải quyết khoảng trống nghiên cứu then chốt này.
+--------------------------------------------------------+
| MÔ HÌNH HỆ THỐNG DFIG NỐI LƯỚI PHÂN PHỐI |
+--------------------------------------------------------+
|
+-------------------------------------------+------------------------------------------+
| |
v v
+---------------------------+ +---------------------------+
| STATOR DFIG (3-PHASE) | | ROTOR DFIG (SLIP RING) |
+---------------------------+ +---------------------------+
| (Đấu nối trực tiếp) | (Đấu qua bộ biến đổi)
| v
| +---------------------------+
| | RSC (Rotor Side Converter)|
| +---------------------------+
| |
| +---------------------------+
| | TỤ ĐIỆN DC-LINK |
| +---------------------------+
| |
| +---------------------------+
| | GSC (Grid Side Converter) |
| +---------------------------+
| |
+-------------------------------------------+------------------------------------------+
|
v
+-------------------------------+
| MBA 3 CUỘN DÂY / TĂNG ÁP |
+-------------------------------+
|
v
+-------------------------------+
| PCC (ĐIỂM ĐẤU NỐI CHUNG 22kV) |
+-------------------------------+
|
v
+-------------------------------+
| LƯỚI ĐIỆN PHÂN PHỐI 110kV/EVN |
+-------------------------------+
Khoảng trống nghiên cứu (Research Gap): Phần lớn các công trình điều khiển DFIG truyền thống phụ thuộc vào bộ điều khiển tỉ lệ - tích phân (PI) được chỉnh định bằng phương pháp thử sai hoặc quy tắc kinh điển Ziegler-Nichols (Z-N), vốn giả định mô hình tuyến tính hóa quanh một điểm làm việc tĩnh. Khi tích hợp vào lưới phân phối có tổng trở $Z_k = R_k + jX_k$ biến động và chịu tác động của nhiễu loạn ngẫu nhiên từ gió, các bộ điều khiển này bộc lộ độ vọt lố lớn, thời gian quá độ kéo dài ($>0,25\text{s}$) và làm suy giảm khả năng phát công suất phản kháng khi điện áp sụt giảm. Mặt khác, các thuật toán trí tuệ nhân tạo như Mạng nơ-ron nhân tạo (ANN), Giải thuật Di truyền (GA), Tối ưu bầy đàn (PSO), Tối ưu sói xám (GWO) hay Tối ưu trọng trường (GSA) tuy đã được áp dụng nhưng thường đối mặt với nhược điểm: kẹt tại cực trị cục bộ, tốc độ hội tụ chậm, chi phí tính toán cao hoặc cấu trúc phức tạp khó triển khai thời gian thực.
Câu hỏi nghiên cứu và Giả thuyết khoa học:
- RQ1: Thiết lập mối quan hệ giải tích chính xác giữa các thông số vật lý của máy điện DFIG ($R_s, R_r, L_m, L_s, L_r$) với thông số của bộ điều khiển PI trên hệ trục tọa độ quay $dq$ như thế nào nhằm đạt được phân ly hoàn toàn giữa điều khiển công suất tác dụng ($P$) và công suất phản kháng ($Q$)?
- RQ2: Thuật toán Tối ưu phản ứng hóa học (Chemical Reaction Optimization - CRO) có khả năng tối ưu hóa đồng thời các cặp thông số ($K_P, K_I$) cho cả bộ biến đổi phía máy phát (RSC) và phía lưới (GSC) để cực tiểu hóa phiếm hàm sai số tích phân theo thời gian nhân sai số tuyệt đối (ITAE) hay không?
- RQ3: Cấu hình điều khiển tối ưu hóa bằng CRO duy trì độ ổn định điện áp, triệt tiêu dao động công suất và đáp ứng các tiêu chuẩn kỹ thuật điện năng khắt khe (TCVN 10687-21, IEC 61400-21, Thông tư 39/2017/TT-BCT) trong các kịch bản sự cố n-1 máy biến áp và suy giảm vận tốc gió đột ngột như thế nào?
Giả thuyết nghiên cứu: Việc tích hợp mô hình giải tích hàm truyền rút gọn bậc hai đa thức Butterworth với cơ chế biến đổi thế năng - động năng của thuật toán tối ưu phản ứng hóa học CRO sẽ tạo ra bộ tham số điều khiển PI tối ưu toàn cục, rút ngắn thời gian quá độ xuống dưới $0,05\text{s}$, duy trì ổn định điện áp DC-link, giảm độ méo sóng hài (THD $< 3%$) và cải thiện vượt bậc đặc tính ổn định tĩnh $P-V$, $Q-V$ tại nút PCC.
Khung lý thuyết: Nghiên cứu dựa trên nền tảng Lý thuyết máy điện quay tổng quát (Park-Gorev transformation), Lý thuyết điều khiển tựa từ thông stator (Stator Flux Oriented Control - SFOC), Lý thuyết điều khiển gán điểm cực (Pole Placement Control) và Lý thuyết tối ưu hóa siêu phỏng sinh học / hóa học phân tử (Molecular Reaction Dynamics).
Phạm vi và Ý nghĩa: Luận án khảo sát đối tượng máy phát DFIG công suất 4 kW (thực nghiệm kiểm chứng) mở rộng lên quy mô công nghiệp 1,5 MW kết nối lưới phân phối 22 kV/110 kV, phân tích toàn diện các chế độ dưới đồng bộ ($s > 0$) và trên đồng bộ ($s < 0$) trong dải trượt $s = \pm 30%$.
Literature Review và Positioning
Cơ sở lý thuyết về điều khiển DFIG nối lưới đã chứng kiến sự phát triển mạnh mẽ qua nhiều thập kỷ với các trường phái tiếp cận khác nhau:
+-------------------------------------------------------------+
| TIẾN TRÌNH PHÁT TRIỂN CÁC PHƯƠNG PHÁP ĐIỀU KHIỂN HỆ DFIG |
+-------------------------------------------------------------+
|
+------------------------------------+---------------+------------------------------------+
| | |
v v v
+-----------------------+ +-----------------------+ +-----------------------+
| ĐIỀU KHIỂN KINH ĐIỂN | | ĐIỀU KHIỂN PHI TUYẾN | | THUẬT TOÁN TIẾN HÓA / |
| (PI / Z-N) | | (BACKSTEPPING / PBC) | | METAHEURISTICS |
+-----------------------+ +-----------------------+ +-----------------------+
| • Nhược điểm: Phụ | | • Nhược điểm: Cấu | | • GA/PSO: Dễ kẹt cực |
| thuộc điểm tĩnh, | | trúc phức tạp, khó | trị cục bộ, hội tụ |
| thời gian xác lập | | thực thi phần cứng, | chậm. |
| lớn (>0,25s). | | chưa tối ưu hóa | | • GWO/GSA: Thời gian |
| • Bị mất ổn định khi | | toàn diện tham số. | xác lập 0,38s - 0,47s.|
| có nhiễu mạnh. | | | | • ANN: Hộp đen, khó |
| | | | giải thích toán học. |
+-----------------------+ +-----------------------+ +-----------------------+
| | |
+------------------------------------+---------------+------------------------------------+
|
v
+-----------------------------------------+
| ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN: |
| KẾT HỢP SFOC dq + GÁN CỰC BUTTERWORTH |
| VỚI TỐI ƯU HÓA PHẢN ỨNG HÓA HỌC CRO |
| (Thời gian xác lập < 0,05s, ITAE min) |
+-----------------------------------------+
Tổng quan các dòng nghiên cứu chính
- Trường phái Điều khiển Tuyến tính và Vector cổ điển: Dựa trên kỹ thuật biến đổi tọa độ Park đưa hệ phương trình trạng thái xoay chiều ba pha về hệ trục tọa độ trực giao $d-q$ tựa theo từ thông stator (Stator Voltage/Flux Orientation - SVO/SFO) được phát triển bởi Leonhard (1985) và Pena et al. (1996). Các bộ điều khiển PI truyền thống áp dụng kỹ thuật giải mã phân ly (cross-coupling decoupling) điều khiển độc lập dòng $i_{dr}$ (công suất phản kháng $Q_s$) và dòng $i_{qr}$ (mô men điện từ $T_e$, công suất tác dụng $P_s$). Tuy nhiên, phương pháp Ziegler-Nichols hiệu chỉnh các hệ số $K_P, K_I$ thường dẫn đến hiện tượng quá dòng khi hòa lưới hoặc quá điện áp mạch DC-link khi có sự cố sụt áp.
- Trường phái Điều khiển Phi tuyến Hiện đại: Nhằm khắc phục tính biến thiên tham số và phi tuyến của DFIG, các kỹ thuật điều khiển trượt (Sliding Mode Control - SMC), điều khiển Backstepping (nghiên cứu trong nước của nhóm tác giả số [11], [12], [16], [17]) và điều khiển tựa thụ động Euler-Lagrange / Hamiltonian (Passivity-Based Control - PBC) [13], [21] đã được áp dụng. Mặc dù chứng minh được tính ổn định tiệm cận Lyapunov, các kỹ thuật này yêu cầu đo lường hoặc ước lượng trạng thái phức tạp, gây hiện tượng rung chấn góc kích đóng IGBT (chattering), đồng thời nghiên cứu của [21] chỉ giải quyết cho trường hợp phụ tải đối xứng, chưa đề cập đến sự cố không đối xứng và sụt giảm điện áp sâu.
- Trường phái Ứng dụng Trí tuệ Nhân tạo và Metaheuristics:
- Nghiên cứu của nhóm tác giả nước ngoài [42], [43] ứng dụng Giải thuật Di truyền (GA) và Tối ưu bầy đàn (PSO) để tìm cực trị $C_{p_max}$ trong điều khiển MPPT. Kết quả từ nghiên cứu [43] cho thấy PSO giúp cải thiện độ méo sóng hài dòng stator (THD) từ 17,67% (PI truyền thống) xuống 15,70%, nhưng thời gian tính toán lặp lại còn lớn.
- Nhóm nghiên cứu [50] ứng dụng thuật toán Tối ưu Trọng trường (Gravitational Search Algorithm - GSA) để tối ưu hóa bộ điều khiển DFIG theo hàm mục tiêu ITAE, đạt thời gian xác lập công suất tác dụng là 0,38s và công suất phản kháng là 0,47s, vượt trội so với GA và PSO (mất tới 0,72s).
- Nghiên cứu [51] áp dụng thuật toán Sói xám (Grey Wolf Optimizer - GWO) và nghiên cứu [52] áp dụng Mạng nơ-ron nhân tạo (ANN) nhằm tự động thích nghi thông số điều khiển. Tuy nhiên, ANN yêu cầu tập dữ liệu huấn luyện lớn, tính chất "hộp đen" làm giảm độ tin cậy trong các phân tích giải tích hệ thống điện nghiêm ngặt.
Định vị học thuật của Luận án
Luận án định vị tại điểm giao thoa giữa Lý thuyết điều khiển tựa từ thông stator giải tích chính xác cao và Thuật toán tối ưu hóa chuyển hóa năng lượng phân tử CRO (Lam & Li, 2010). So với nghiên cứu quốc tế của nhóm tác giả áp dụng GSA [50] (thời gian đáp ứng 0,38s) hay nhóm nghiên cứu sử dụng Fuzzy-Sliding Mode [20] (thời gian quá độ 0,25s), giải pháp đề xuất của luận án vượt lên một bước đột phá: đạt thời gian đáp ứng quá độ dưới 0,05s, triệt tiêu hoàn toàn dao động công suất phản kháng khi vận tốc gió biến thiên và duy trì điện áp thanh cái một chiều DC-link ổn định tuyệt đối quanh 1200V.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đã đóng góp vào lý thuyết Kỹ thuật Điều khiển Tự động và Kỹ thuật Hệ thống Điện ba luận điểm cốt lõi:
+--------------------------------------------------------------------------------+
| CHU TRÌNH 4 BƯỚC PHẢN ỨNG HÓA HỌC (CRO) TỐI ƯU HÓA BỘ ĐIỀU KHIỂN PI |
+--------------------------------------------------------------------------------+
|
+-----------------------------------------------+-------------------------------+
| (Phản ứng đơn phân tử) | (Phản ứng đa phân tử) |
v v v
+-------------------------------+ +--------------------+ +--------------------+
| VA CHẠM KHÔNG HIỆU QUẢ VỚI | | PHÂN HỦY PHÂN TỬ | | VA CHẠM ĐA PHÂN TỬ |
| THÀNH BÌNH (On-wall Collision)| | (Decomposition) | | KHÔNG HIỆU QUẢ |
+-------------------------------+ +--------------------+ +--------------------+
| Biến đổi cấu trúc cục bộ lân | | Tách 1 phân tử | | 2 phân tử va chạm |
| cận giải pháp: | | thành 2 phân tử | | trao đổi năng |
| $\omega' = \omega + N(0, \sigma)$| | mới để mở rộng | | lượng và cập nhật |
| Thế năng biến thành động năng.| | không gian tìm kiếm| | vùng nghiệm lân cận|
+-------------------------------+ +--------------------+ +--------------------+
| | |
+-----------------------------------------------+-------------------------------+
|
v
+------------------------------------+
| TỔNG HỢP PHÂN TỬ (Synthesis) |
+------------------------------------+
| Gộp 2 phân tử thành 1 phân tử có |
| thế năng thấp nhất (ITAE cực tiểu) |
| Hội tụ về bộ thông số tối ưu toàn |
| cục ($K_P^*, K_I^*$). |
+------------------------------------+
-
Thiết lập công thức toán học giải tích tường minh liên kết tham số điều khiển với thông số động học máy điện DFIG: Thay vì xem hệ thống kín là một đối tượng chưa biết hoàn toàn, luận án sử dụng phương pháp gán điểm cực Butterworth chuẩn hóa bậc hai: $$s^2 + 2\xi\omega_n s + \omega_n^2 = 0$$ Từ đó suy biến các biểu thức giải tích xác định miền giới hạn vật lý của $K_P$ và $K_I$ cho mạch vòng dòng điện trục $d-q$ và mạch vòng điện áp/công suất: $$K_{P_idr} = 2\xi\omega_n \sigma L_r - R_r$$ $$K_{I_idr} = \omega_n^2 \sigma L_r$$ với $\sigma = 1 - \frac{L_m^2}{L_s L_r}$ là hệ số tản từ. Điều này biến bài toán tìm kiếm từ không gian vô hạn thành bài toán tối ưu có ràng buộc biên xác định trên không gian trạng thái thực.
-
Chuyển dịch mô hình tối ưu hóa điều khiển DFIG sang Không gian Năng lượng Thế năng (Potential Energy Surface - PES): Ứng dụng nguyên lý thứ hai nhiệt động lực học trong thuật toán CRO, xem mỗi nghiệm $[K_P, K_I]$ là một cấu trúc phân tử $\omega$, thế năng của phân tử $PE(\omega)$ chính là giá trị của hàm mục tiêu ITAE: $$J_{\text{ITAE}} = \int_0^T t \cdot |e(t)| , dt$$ Động năng $KE(\omega)$ đại diện cho khả năng thoát khỏi các điểm cực trị địa phương (local minima).
-
Mô hình hóa giải tích quan hệ phi tuyến bậc 4 giữa điện áp nút PCC ($U_T$) và công suất phát ($P, Q$): $$U_T^4 + \left[2(PR_k + QX_k) - U_N^2\right]U_T^2 + (P^2 + Q^2)(R_k^2 + X_k^2) = 0$$ Chứng minh một cách chặt chẽ rằng sự biến thiên công suất tác dụng $\Delta P$ do gió gây ra có thể được bù đắp hoàn toàn bằng việc điều khiển phát/hấp thụ công suất phản kháng $Q$ của máy phát DFIG thông qua mạch rotor để ổn định điện áp thanh cái $U_T$.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đa tầng ba lý thuyết nền tảng:
- Lý thuyết Chuyển đổi Năng lượng Khí động học Betz-Lanchester: Xác định hệ số công suất giới hạn cực đại $C_{p_max} = 59,3%$ thông qua hàm phi tuyến của hệ số tốc độ đầu cánh $\lambda = \frac{\omega_m R}{v}$ và góc nghiêng pitch $\beta$: $$C_p(\lambda, \beta) = c_1 \left(\frac{c_2}{\lambda_i} - c_3\beta - c_4\right) e^{-\frac{c_5}{\lambda_i}} + c_6 \lambda$$
- Lý thuyết Biến đổi Tọa độ Trường vector xoay dq0: Khử ghép từ thông rotor $\psi_r$ và stator $\psi_s$ trong điều kiện tựa theo từ thông stator ($\psi_{ds} = \psi_s, \psi_{qs} = 0$), đưa hệ phương trình vi phân phi tuyến bậc cao về hai kênh điều khiển đơn biến độc lập: $$T_e = -\frac{3}{2} p \frac{L_m}{L_s} \psi_{ds} i_{qr}$$ $$Q_s = \frac{3}{2} \frac{v_s \psi_{ds}}{L_s} - \frac{3}{2} \frac{v_s L_m}{L_s} i_{dr}$$
- Lý thuyết Tối ưu hóa Động lực học Phản ứng Hóa học (CRO): Xây dựng 4 toán tử biến đổi phân tử cốt lõi:
- Va chạm không hiệu quả với thành bình (On-wall Ineffective Collision): Tìm kiếm cục bộ lân cận cấu trúc phân tử hiện tại.
- Phân hủy phân tử (Decomposition): Khi phân tử bị bẫy trong vùng thế năng sâu mà không cải thiện được ITAE, nó được tách thành 2 phân tử mới để khám phá vùng nghiệm mới.
- Va chạm đa phân tử không hiệu quả (Inter-molecular Ineffective Collision): Hai phân tử tương tác chia sẻ động năng, mở rộng không gian tìm kiếm.
- Tổng hợp phân tử (Synthesis): Hai phân tử kết hợp thành một phân tử ổn định có thế năng tối thiểu khi tìm được cấu hình tham số $K_P, K_I$ tiệm cận tối ưu toàn cục.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu áp dụng thế giới quan thực chứng (Positivism) kết hợp phương pháp mô hình hóa toán học giải tích cấu trúc sâu và kỹ thuật mô phỏng động học thời gian thực (Time-domain Simulation) trên nền tảng phần mềm tiêu chuẩn công nghiệp MATLAB/Simulink và công cụ chuyên dụng phân tích hệ thống điện (Power System Analysis Toolbox).
+---------------------------------------------------------+
| SƠ ĐỒ KHỐI CẤU TRÚC ĐIỀU KHIỂN ĐỘC LẬP RSC VÀ GSC |
+---------------------------------------------------------+
|
+--------------------------------------------------+--------------------------------------------------+
| |
v (BỘ BIẾN ĐỔI PHÍA MÁY PHÁT - RSC) v (BỘ BIẾN ĐỔI PHÍA LƯỚI - GSC)
+-------------------------------------------------------+ +-------------------------------------------------------+
| MẠCH VÒNG NGOÀI (Outer Loops): | | MẠCH VÒNG NGOÀI (Outer Loops): |
| • Đo tốc độ rotor $\omega_r \rightarrow$ So sánh | | • Đo điện áp DC-link $V_{dc} \rightarrow$ Bộ PI CRO |
| $\omega_r^* \rightarrow$ Bộ PI CRO $\rightarrow i_{qr}^*$ | | $\rightarrow i_{dg}^*$ |
| • Đo công suất phản kháng $Q_s \rightarrow$ So sánh | | • Đo công suất phản kháng lưới $Q_g \rightarrow$ Bộ |
| $Q_s^* \rightarrow$ Bộ PI CRO $\rightarrow i_{dr}^*$ | | PI CRO $\rightarrow i_{qg}^*$ |
| | | |
| MẠCH VÒNG TRONG (Inner Loops): | | MẠCH VÒNG TRONG (Inner Loops): |
| • So sánh $(i_{dr}^*, i_{qr}^*)$ với $(i_{dr}, i_{qr}) \rightarrow$ | | • So sánh $(i_{dg}^*, i_{qg}^*)$ với $(i_{dg}, i_{qg}) \rightarrow$ |
| Bộ PI CRO + Khử ghép giải tích $\rightarrow (v_{dr}^*, v_{qr}^*)$ | | Bộ PI CRO + Khử ghép $\rightarrow (v_{dg}^*, v_{qg}^*)$ |
| | | |
| KỸ THUẬT PHÁT XUNG (PWM): | | KỸ THUẬT PHÁT XUNG (PWM): |
| • Biến đổi $dq \rightarrow abc \rightarrow$ Điều chế PWM IGBT| | • Biến đổi $dq \rightarrow abc \rightarrow$ Điều chế PWM IGBT|
+-------------------------------------------------------+ +-------------------------------------------------------+
| |
+--------------------------------------------------+--------------------------------------------------+
|
v
+-----------------------------------+
| MÁY PHÁT ĐIỆN DFIG & LƯỚI 22/110kV|
+-----------------------------------+
Quy trình nghiên cứu rigorous
Quy trình thực hiện được kiểm soát nghiêm ngặt qua 4 giai đoạn logic:
- Thiết lập mô hình trạng thái DFIG liên tục: Xây dựng hệ phương trình vi phân mô tả điện áp và từ thông trên hệ trục dq: $$v_{ds} = R_s i_{ds} + \frac{d\psi_{ds}}{dt} - \omega_s \psi_{qs}$$ $$v_{qs} = R_s i_{qs} + \frac{d\psi_{qs}}{dt} + \omega_s \psi_{ds}$$ $$v_{dr} = R_r i_{dr} + \frac{d\psi_{dr}}{dt} - (\omega_s - \omega_r) \psi_{qr}$$ $$v_{qr} = R_r i_{qr} + \frac{d\psi_{qr}}{dt} + (\omega_s - \omega_r) \psi_{dr}$$
- Khởi tạo và cấu hình không gian tìm kiếm CRO: Thiết lập quần thể phân tử $N = 20 - 50$, các tham số động học: động năng ban đầu $Initial_KE$, ngưỡng phản ứng phân hủy $\alpha = 500$, ngưỡng tổng hợp $\beta = 10$, tỉ lệ suy giảm động năng $LossRate = 0,99$.
- Mô phỏng đa kịch bản (Scenario Testing):
- Kịch bản vận hành xác lập tại các vận tốc gió khác nhau ($v = 7\text{ m/s}, 11\text{ m/s}, 14\text{ m/s}$).
- Kịch bản sự cố nghiêm trọng: Mất một máy biến áp 110 kV trong trạm biến áp kép (sự cố $n-1$).
- Kịch bản sụt giảm vận tốc gió đột ngột (giảm tức thời 30% công suất phát).
- Phân tích ổn định tín hiệu nhỏ (Small-signal Stability Analysis): Tính toán giá trị riêng (Eigenvalues $\lambda_i = \sigma_i \pm j\omega_i$) của ma trận trạng thái $A$ và phân tích hệ số tham gia (Participation Factors) để đánh giá độ cản dịu dao động cơ điện.
Data và phân tích
Thông số hệ thống kiểm chứng cụ thể trong luận án:
- Hệ thống DFIG 4 kW: $P_n = 4\text{ kW}$, $V_n = 380\text{ V}$, $f = 50\text{ Hz}$, $R_s = 1,47,\Omega$, $R_r = 1,32,\Omega$, $L_s = 0,268\text{ H}$, $L_r = 0,268\text{ H}$, $L_m = 0,257\text{ H}$, số cặp cực $p = 2$.
- Hệ thống DFIG quy mô công nghiệp 1,5 MW: $P_n = 1,5\text{ MW}$, $V_s = 690\text{ V}$, $V_{dc} = 1200\text{ V}$, $f = 50\text{ Hz}$, $R_s = 0,00706\text{ p.u.}$, $R_r = 0,005\text{ p.u.}$, $L_{ls} = 0,171\text{ p.u.}$, $L_{lr} = 0,156\text{ p.u.}$, $L_m = 2,9\text{ p.u.}$, quán tính $H = 5,04\text{ s}$.
- Lưới điện phân phối mô phỏng: Trạm biến áp 110/22 kV gồm 2 máy biến áp làm việc song song công suất $2 \times 25\text{ MVA}$, đường dây trung áp 22 kV cấp điện cho phụ tải địa phương $P_{load} = 12\text{ MW}, Q_{load} = 3\text{ MVAR}$.
Bảng so sánh tham số điều khiển và chỉ số hàm mục tiêu ITAE giữa các phương pháp:
| Vị trí bộ điều khiển | Phương pháp Ziegler-Nichols (Z-N) | Phương pháp Tối ưu CRO (Luận án đề xuất) | Mức độ cải thiện ITAE (%) |
|---|---|---|---|
| Mạch vòng dòng $i_{dr}$ | $K_P = 4,215;; K_I = 125,4$ | $K_P = 12,842;; K_I = 312,65$ | $\downarrow 68,4%$ |
| Mạch vòng dòng $i_{qr}$ | $K_P = 3,890;; K_I = 110,2$ | $K_P = 11,954;; K_I = 298,40$ | $\downarrow 71,2%$ |
| Mạch vòng áp DC-link | $K_P = 8,120;; K_I = 45,6$ | $K_P = 24,350;; K_I = 115,80$ | $\downarrow 64,8%$ |
| Mạch vòng công suất $Q$ | $K_P = 1,250;; K_I = 15,2$ | $K_P = 5,620;; K_I = 48,90$ | $\downarrow 75,6%$ |
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Khả năng triệt tiêu hoàn toàn quá độ và rút ngắn thời gian xác lập: Dưới sự tối ưu hóa của thuật toán CRO, thời gian đáp ứng quá độ của dòng điện rotor ($i_{dr}, i_{qr}$) và công suất tác dụng ($P$), phản kháng ($Q$) giảm mạnh từ $0,25\text{s}$ (đối với phương pháp mờ trượt hoặc PI-ZN) xuống chỉ còn dưới $0,045\text{s}$ (nhanh hơn $5,5$ lần), không xuất hiện hiện tượng vọt lố (overshoot $\approx 0%$).
-
Duy trì ổn định tuyệt đối điện áp DC-link khi có sự cố nghiêm trọng: Khi mô phỏng kịch bản sự cố đứt một máy biến áp 110 kV tại trạm chính, điện áp thanh cái một chiều $V_{dc}$ chỉ dao động trong dải cực hẹp $\Delta V_{dc} < \pm 1,2%$ và trở về giá trị xác lập 1200V trong vòng $0,08\text{s}$, ngăn chặn hoàn toàn nguy cơ kích hoạt mạch bảo vệ Crowbar gây mất kết nối máy phát.
-
Cải thiện độ méo sóng hài (THD) và ổn định điện áp nút 22 kV: Độ méo sóng hài tổng của dòng điện phát lên lưới giảm xuống mức $1,85%$, thấp hơn nhiều so với ngưỡng khống chế $3%$ của Tiêu chuẩn IEC 61400-21 và Thông tư 39/2017/TT-BCT của Bộ Công Thương.
-
Nâng cao biên độ ổn định tĩnh trên đường cong P-V và Q-V: Đường đặc tính $P-V$ và $Q-V$ tại nút 22 kV chứng minh rằng giải thuật điều khiển tối ưu CRO giúp tăng giới hạn truyền tải công suất cực đại ($P_{max}$) thêm $14,6%$ trước khi chạm điểm sụp đổ điện áp (nose point), đồng thời tăng khả năng bơm công suất phản kháng tự động hỗ trợ lưới khi xảy ra ngắn mạch lân cận.
-
Phát hiện về ma trận giá trị riêng (Eigenvalue spectrum): Phần thực của tất cả các giá trị riêng cơ điện đều dịch chuyển sâu về phía bên trái trục ảo ($\sigma_i < -12,5$), hệ số cản dịu (damping ratio) đạt $\zeta > 0,15$, loại trừ hoàn toàn các dao động tần số thấp ($0,2 - 2\text{ Hz}$) giữa máy phát gió và hệ thống điện phân phối.
Implications đa chiều
- Ý nghĩa Lý thuyết: Khẳng định tính tương thích vượt trội của mô hình tối ưu hóa phân tử CRO trong việc giải các bài toán phi tuyến đa biến trong kỹ thuật điều khiển máy điện hiện đại.
- Ý nghĩa Thực tiễn: Cung cấp bộ thông số cài đặt ($K_P, K_I$) chuẩn mực cho các kỹ sư vận hành nhà máy điện gió sử dụng công nghệ DFIG mà không cần phải trải qua quá trình chỉnh định thử sai tốn kém trên hiện trường.
- Ý nghĩa Chính sách & Tiêu chuẩn: Làm căn cứ kỹ thuật để các cơ quan quản lý (EVN, Cục Điều tiết Điện lực) xây dựng quy chuẩn kỹ thuật đấu nối (Grid Code) cho các nguồn điện phân tán quy mô lớn.
Limitations và Future Research
- Giới hạn về mô hình bão hòa từ và sóng hài bậc cao: Luận án xây dựng mô hình dựa trên giả thiết khe hở không khí đồng nhất, bỏ qua hiện tượng bão hòa mạch từ và dòng rò phi tuyến ở tần số cao của máy điện.
- Giới hạn về điều kiện mất cân bằng pha cực hạn: Nghiên cứu tập trung giải quyết triệt để sự cố đối xứng và suy giảm vận tốc gió; các kịch bản sự cố chạm đất một pha không đối xứng sâu kèm theo biến dạng điện áp phi sin nghiêm trọng chưa được khảo sát chi tiết.
- Môi trường thử nghiệm thực nghiệm: Luận án mới chỉ thực nghiệm trên mô hình DFIG 4 kW phòng thí nghiệm; kết quả trên mô hình 1,5 MW được kiểm chứng thông qua mô phỏng số chi tiết (numerical simulation) mà chưa chạy thử nghiệm trên tuabin gió thực tế ngoài khơi do hạn chế kinh phí.
Hướng nghiên cứu tương lai:
- Tích hợp cấu trúc điều khiển CRO thời gian thực (Online Adaptive Real-time CRO) sử dụng chip xử lý tín hiệu số DSP/FPGA.
- Nghiên cứu mở rộng thuật toán CRO cho các cấu hình tuabin gió dùng máy phát đồng bộ nam châm vĩnh cửu (PMSG) kết nối bộ biến đổi toàn phần (Full-scale converter).
- Đánh giá khả năng tích hợp hệ thống lưu trữ năng lượng (BESS - Battery Energy Storage Systems) phối hợp cùng DFIG để tạo khả năng khởi động đen (Black-start) và duy trì tần số hệ sinh thái lưới vi mô (Microgrid).
Tác động và ảnh hưởng
- Tác động Học thuật: Cung cấp tài liệu tham khảo giá trị cao cho các chương trình đào tạo tiến sĩ, thạc sĩ chuyên ngành Kỹ thuật Điện, Kỹ thuật Điều khiển và Tự động hóa; mở ra hướng nghiên cứu mới về việc ứng dụng vật lý/hóa học tính toán vào kỹ thuật năng lượng.
- Chuyển dịch Công nghiệp: Giúp các chủ đầu tư dự án điện gió và các đơn vị EPC tối ưu hóa cấu hình vận hành bộ biến đổi, giảm thiểu tổn thất dừng máy do quá áp DC-link, ước tính nâng cao sản lượng điện phát hàng năm từ 1,5% đến 2,8%.
- Lợi ích Xã hội & Môi trường: Tăng cường khả năng giải tỏa công suất từ các nguồn năng lượng sạch, giảm thiểu việc cắt giảm công suất điện gió (curtailment), trực tiếp đóng góp vào cam kết Net Zero vào năm 2050 của Việt Nam tại COP26.
Đối tượng hưởng lợi
+-------------------------------------------------------------+
| CÁC ĐỐI TƯỢNG HƯỞNG LỢI CHÍNH |
+-------------------------------------------------------------+
|
+------------------------------------+---------------+------------------------------------+
| | |
v v v
+-----------------------+ +-----------------------+ +-----------------------+
| NGHIÊN CỨU SINH & | | KỸ SƯ VẬN HÀNH & | | CƠ QUAN HOẠCH ĐỊNH |
| HỌC GIẢ | | CÁC CHỦ ĐẦU TƯ ĐIỆN | | CHÍNH SÁCH & QUẢN LÝ |
+-----------------------+ +-----------------------+ +-----------------------+
| • Tiếp cận phương | | • Sở hữu bộ tham số | | • EVN, Điều độ Quốc |
| pháp tối ưu CRO | chỉnh định tối ưu | gia (A0): Có cơ sở kỹ |
| trong điều khiển | chuẩn xác cho DFIG. | thuật hoàn thiện Grid |
| phi tuyến máy điện. | | • Giảm sự cố cắt máy, | Code cho điện gió |
| • Kế thừa khung giải | kéo dài tuổi thọ tụ | phân tán. |
| tích dq-Butterworth.| DC-link & IGBT. | | • Tối ưu hóa vận hành |
+-----------------------+ +-----------------------+ +-----------------------+
- Nghiên cứu sinh và Giảng viên ngành Điện - Năng lượng tái tạo: Tiếp cận một phương pháp luận kết hợp hoàn hảo giữa mô hình hóa toán học giải tích cổ điển và thuật toán metaheuristic tiên tiến.
- Kỹ sư vận hành tại các Trang trại Điện gió (Wind Farms): Nắm bắt được bản chất tương tác động học giữa RSC, GSC với lưới phân phối để thiết lập chính xác các thông số bảo vệ rơ-le và tham số điều khiển bộ nghịch lưu.
- Tập đoàn Điện lực Việt Nam (EVN) và Trung tâm Điều độ Hệ thống điện Quốc gia (A0): Có luận cứ khoa học để thẩm định các yêu cầu kỹ thuật hòa lưới của các nhà máy điện gió, nâng cao độ tin cậy cung cấp điện.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của nghiên cứu là gì?
Đóng góp độc đáo nhất là việc thiết lập được biểu thức giải tích tường minh ràng buộc miền tham số điều khiển PI thông qua ma trận trạng thái DFIG kết hợp chuẩn đa thức Butterworth, sau đó ánh xạ trực tiếp bài toán tối ưu sai số động học ITAE vào không gian thế năng của thuật toán Tối ưu Phản ứng Hóa học (CRO), giải quyết triệt để sự đánh đổi giữa tốc độ đáp ứng và độ vọt lố trong điều khiển máy điện.
2. Đổi mới về phương pháp luận so với các nghiên cứu trước đây?
So với phương pháp kinh điển Ziegler-Nichols [22] và phương pháp mờ trượt [20] (thời gian quá độ $> 0,25\text{s}$), hay các thuật toán metaheuristics như GSA [50] (đáp ứng $0,38\text{s}$) và PSO [43] (dòng stator THD còn $15,7%$), phương pháp CRO của luận án vượt trội nhờ cơ chế chuyển hóa động năng - thế năng linh hoạt với 4 toán tử phân tử, cho phép quét toàn cục không gian nghiệm mà không bị mắc kẹt tại cực trị địa phương, rút ngắn thời gian quá độ xuống $0,045\text{s}$ và giảm THD xuống $1,85%$.
3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?
Phát hiện bất ngờ nhất là khi xảy ra sự cố mất 1 máy biến áp 110 kV tại trạm chính, bộ điều khiển GSC được tối ưu bằng CRO có khả năng tự động điều chỉnh góc pha dòng điện nghịch lưu để ghìm điện áp DC-link hầu như không biến động ($\Delta V_{dc} < 1,2%$), loại bỏ hoàn toàn sự can thiệp của mạch xả năng lượng Crowbar, giúp tuabin gió tiếp tục duy trì phát điện liên tục qua sự cố (Fault Ride-Through - FRT).
4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) không?
Có. Luận án trình bày toàn bộ hệ phương trình toán học giải tích hệ tọa độ $dq$, bảng thông số vật lý chi tiết của máy phát 4 kW và 1,5 MW, thông số đường dây, máy biến áp 3 cuộn dây, cùng mã giả (pseudocode) và lưu đồ thuật toán CRO với đầy đủ các hệ số toán tử ($\alpha, \beta, Initial_KE, LossRate$).
5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Lộ trình 10 năm mở rộng nghiên cứu bao gồm: (1) Phát triển chip điều khiển nhúng chuyên dụng tích hợp giải thuật CRO thời gian thực; (2) Mở rộng điều khiển phối hợp cho các trang trại điện gió ngoài khơi (Offshore Wind Farms) kết nối lưới truyền tải siêu cao áp một chiều HVDC; (3) Tích hợp công nghệ lưới điện thông minh tạo lập lưới (Grid-forming Control) cho DFIG trong bối cảnh hệ thống điện có quán tính thấp (Low-inertia Power Systems).
Kết luận
Luận án tiến sĩ của tác giả Nguyễn Công Cường đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra với 6 đóng góp cụ thể:
- Xây dựng hoàn chỉnh mô hình toán học giải tích phi tuyến của hệ thống máy phát điện không đồng bộ nguồn kép DFIG trên hệ trục tọa độ $dq$, làm rõ cơ chế phân bố dòng công suất $P_r \approx -s P_s$ giữa stator và rotor.
- Ứng dụng thành công thuật toán Tối ưu Hóa học CRO để xác định bộ thông số tối ưu toàn cục ($K_P, K_I$) cho các bộ điều khiển RSC và GSC, cực tiểu hóa hàm mục tiêu ITAE.
- Rút ngắn thời gian xác lập quá độ của hệ thống xuống dưới $0,05\text{s}$, loại bỏ hoàn toàn hiện tượng vọt lố và dao động cơ điện khi tốc độ gió biến thiên.
- Đảm bảo chất lượng điện năng vượt tiêu chuẩn quốc tế (IEC 61400-21) và quy định của Việt Nam (Thông tư 39/2017/TT-BCT) với độ méo sóng hài $\text{THD} = 1,85%$ và độ dao động điện áp cực tiểu.
- Chứng minh tính ổn định vững chắc của hệ thống qua phân tích ma trận giá trị riêng và nâng cao giới hạn ổn định tĩnh trên đường cong $P-V, Q-V$.
- Mở ra phương pháp tiếp cận mới trong việc ứng dụng thuật toán phỏng sinh học và hóa học phân tử để giải quyết các bài toán tối ưu hóa điều khiển phức tạp trong các nguồn năng lượng tái tạo phân tán.