Tổng quan về luận án

Sự bùng nổ của các nguồn điện phân tán (Distributed Generation - DG) sử dụng năng lượng tái tạo như quang điện (Photovoltaic - PV) và điện gió đang định hình lại cấu trúc của hệ thống điện hiện đại trên toàn cầu. Theo các định hướng chuyển dịch năng lượng và Quy hoạch phát triển điện lực quốc gia (Quy hoạch điện VII điều chỉnh), tỷ trọng công suất lắp đặt từ năng lượng tái tạo được nâng cao liên tục nhằm đáp ứng nhu cầu tiêu thụ và giảm phát thải khí nhà kính. Tuy nhiên, rào cản kỹ thuật cốt lõi của các nguồn năng lượng này là tính bất định, phân tán và gián đoạn. Để hòa lưới hiệu quả, các hệ thống bắt buộc phải sử dụng các bộ biến đổi bán dẫn công suất, trọng tâm là nghịch lưu nối lưới (Grid-Connected Inverter).

Quá trình chuyển mạch tần số cao của các linh kiện bán dẫn (như IGBT, MOSFET) trong bộ nghịch lưu lại chính là nguồn phát sinh sóng hài dòng điện đáng kể, lan truyền vào điểm ghép nối chung (Point of Common Coupling - PCC) và làm suy giảm nghiêm trọng chất lượng điện năng của mạng phân phối. Việc vi phạm các tiêu chuẩn nối lưới nghiêm ngặt quốc tế và quốc gia (như IEEE-929, IEEE-1547, IEC 61727, EN 50160, VDE 0126 và Thông tư 39/2015/TT-BCT của Bộ Công Thương Việt Nam) có thể gây ra hiện tượng quá nhiệt máy biến áp, rung giật mô-men động cơ xoay chiều, sai lệch thiết bị đo đếm điện tử kỹ thuật số và tác động sai của rơ-le bảo vệ hệ thống điện.

Khoảng trống nghiên cứu (Research Gap) then chốt xuất phát từ sự đánh đổi phức tạp giữa việc giảm độ méo hài toàn phần (Total Harmonic Distortion - THD) và việc tối ưu hóa hiệu suất bộ biến đổi. Các giải pháp truyền thống dựa vào việc tăng kích thước bộ lọc thụ động (L, LC, LCL) làm gia tăng thể tích, trọng lượng và giá thành thiết bị, đồng thời gây suy giảm tính ổn định động và khả năng đáp ứng tiêu chuẩn trải qua sự cố (Fault-Ride Through - FRT). Ngược lại, việc tăng tần số chuyển mạch cố định dẫn đến tổn hao chuyển mạch tăng tuyến tính ($P_{sw} \propto f_{sw}$), gây quá nhiệt linh kiện bán dẫn và phát sinh nhiễu điện từ (EMI) tập trung ở dải tần cao.

Luận án tiến sĩ kỹ thuật điện "Điều khiển bộ nghịch lưu nối lưới trong mạng điện phân phối" của Nghiên cứu sinh Trần Quang Thọ (Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh, 2017) đã giải quyết toàn diện bài toán này thông qua việc thiết lập 4 câu hỏi nghiên cứu ($RQ$) và 4 giả thuyết khoa học ($H$):

  • Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để phân bố lại phổ sóng hài và triệt tiêu năng lượng hài bậc cao mà không làm gia tăng tổng tổn hao chuyển mạch và không cần bổ sung phần cứng mạch lực?
    Giả thuyết 1 ($H_1$): Một quy luật biến thiên chu kỳ chuyển mạch ($T_{s\text{-var}}$) trong mỗi nửa chu kỳ cơ bản (NCKCB) được tối ưu hóa bằng giải thuật di truyền (Genetic Algorithm - GA) sẽ phân tán năng lượng hài và giảm biên độ sóng hài đỉnh mà vẫn bảo toàn tổng năng lượng tổn hao chuyển mạch.
  • Câu hỏi nghiên cứu 2 ($RQ_2$): Giải pháp nào thay thế cấu trúc vòng khóa pha (Phase-Locked Loop - PLL) truyền thống để ước lượng tức thời, chính xác góc pha, biên độ và tần số lưới trong điều kiện điện áp lưới bị méo và mất cân bằng?
    Giả thuyết 2 ($H_2$): Việc áp dụng giải thuật phi tuyến Levenberg-Marquardt kết hợp với cơ chế thích nghi hồi tiếp ngõ ra để cập nhật điều kiện biên ban đầu sẽ triệt tiêu hoàn toàn sự phụ thuộc vào khâu tích phân dao động VCO và bộ lọc vòng, giúp giảm thời gian xác lập dưới 0,02 s mà không có độ vọt lố.
  • Câu hỏi nghiên cứu 3 ($RQ_3$): Làm thế nào để loại bỏ thành phần sóng hài bậc thấp thâm nhập từ dao động điện áp DC-link của hệ thống pin mặt trời một giai đoạn mà không làm suy giảm tốc độ bám điểm công suất cực đại?
    Giả thuyết 3 ($H_3$): Giải thuật MPPT với số gia biến đổi tự thích ứng theo độ dốc đặc tính $P-V$ sẽ triệt tiêu độ nhấp nhô điện áp một chiều $V_{dc}$ tại điểm xác lập, loại trừ nguồn phát sinh hài từ phía sơ cấp.
  • Câu hỏi nghiên cứu 4 ($RQ_4$): Làm thế nào để xác định tập tham số tối ưu toàn cục cho bộ điều khiển dòng điện cộng hưởng tỉ lệ (Proportional Resonant - PR) khi hàm truyền hệ thống có tính chất vô định và chịu ảnh hưởng của biến thiên tổng trở lưới?
    Giả thuyết 4 ($H_4$): Giải thuật tối ưu bầy đàn (Particle Swarm Optimization - PSO) đa chiều với hàm mục tiêu ràng buộc tích phân sai số bình phương sẽ tìm được tọa độ tham số tối ưu toàn cục, khắc phục hiện tượng rơi vào cực trị địa phương của các phương pháp giải tích Bode/Nyquist truyền thống.

Phạm vi nghiên cứu của luận án tập trung vào cấu trúc bộ nghịch lưu nối lưới nguồn áp (VSI) 1 pha cầu H và 3 pha 2 bậc cho hệ thống điện mặt trời quy mô vừa và nhỏ, xem xét toàn diện các chế độ phát công suất tác dụng ($P$) và thu/phát công suất phản kháng ($Q$). Mô hình nghiên cứu được kiểm chứng định lượng trên nền tảng mô phỏng số MATLAB/Simulink và nền tảng thực nghiệm điều khiển thời gian thực với phần cứng chuyên dụng Texas Instruments DSP TMS320F28335 và dSPACE CP1103.

+-------------------------------------------------------------------------------------------------+
|                       KHUNG ĐIỀU KHIỂN NGHỊCH LƯU NỐI LƯỚI TOÀN DIỆN                            |
+-------------------------------------------------------------------------------------------------+
|                                                                                                 |
|   +--------------------------+        +--------------------------+                              |
|   |   1. ĐIỀU CHẾ KHÔNG GIAN |        |  2. ĐỒNG BỘ HÓA LƯỚI     |                              |
|   |   Biến đổi Ts-var qua GA |        |  Ước lượng Levenberg-    |                              |
|   |   -> Trải phổ hài, giảm  |        |  Marquardt thích nghi    |                              |
|   |      kích thước lọc Lf   |        |  -> Không dùng VCO/PLL   |                              |
|   +-------------+------------+        +-------------+------------+                              |
|                 |                                   |                                           |
|                 +-----------------+-----------------+                                           |
|                                   |                                                             |
|                                   v                                                             |
|                    +------------------------------+                                             |
|                    |     BỘ NGHỊCH LƯU NỐI LƯỚI   |                                             |
|                    |  Giảm THD dòng điện ngõ ra   |                                             |
|                    |  Đạt tiêu chuẩn IEEE 1547    |                                             |
|                    +--------------+---------------+                                             |
|                                   ^                                                             |
|                 +-----------------+-----------------+                                           |
|                 |                                   |                                           |
|   +-------------+------------+        +-------------+------------+                              |
|   |   3. ỔN ĐỊNH DC-LINK     |        |  4. BỘ ĐIỀU KHIỂN DÒNG   |                              |
|   |   MPPT số gia biến đổi   |        |  Tối ưu hóa PR qua PSO   |                              |
|   |   -> Triệt tiêu nhấp     |        |  -> Bám dòng I_alpha/beta|                              |
|   |      nhô Vdc tại cực đại |        |     Triệt tiêu sai lệch  |                              |
|   +--------------------------+        +--------------------------+                              |
|                                                                                                 |
+-------------------------------------------------------------------------------------------------+

Literature Review và Positioning

Các nghiên cứu giảm sóng hài cho nghịch lưu nối lưới trong y văn quốc tế được chia thành 4 trường phái chính: kỹ thuật điều chế nâng cao, công nghệ khóa pha đồng bộ lưới, cấu trúc điều khiển bám dòng và thuật toán ổn định nguồn sơ cấp DC.

Trong lĩnh vực điều chế độ rộng xung (PWM), trường phái truyền thống sử dụng điều chế sóng sin SPWM hoặc điều chế vector không gian SVPWM với tần số sóng mang cố định (Blaabjerg et al., Holmes & Lipo). Đặc trưng của phương pháp này là các thành phần sóng hài tập trung dày đặc xung quanh tần số đóng cắt và các dải biên (sidebands), đòi hỏi bộ lọc ngõ ra phải có bậc cắt sắc nét. Để giải quyết vấn đề này, các nghiên cứu về tần số chuyển mạch biến đổi (Variable Switching Frequency - VSF) đã xuất hiện.

Điển hình là nghiên cứu của nhóm tác giả quốc tế dựa trên mô hình yêu cầu độ méo toàn phần TDD (Total Distortion Demand) để xác định chu kỳ chuyển mạch tối ưu nhằm giữ THD dưới ngưỡng 5% của chuẩn IEEE 519. Tuy nhiên, phương pháp này chỉ điều chỉnh tần số theo mức công suất tải mà vẫn duy trì tần số cố định trong từng chu kỳ cơ bản, đồng thời phụ thuộc chặt chẽ vào độ chính xác của mô hình tham số điện cảm lưới $L$ và điện áp lưới $V_{grid}$.

Một nhánh nghiên cứu khác là kỹ thuật bơm sóng hài vào sóng điều chế và điều chế tần số sóng mang hình tam giác (HIPWP-FMTC), điều chỉnh tần số sóng mang tỉ lệ với độ dốc của tín hiệu điều chế sin. Hạn chế của HIPWP-FMTC là khi triển khai ở dải tần số đóng cắt cao, độ dốc sóng điều chế bị suy giảm cục bộ, dẫn đến mất khả năng kiểm soát sóng hài biên và không tối ưu được tổn hao chuyển mạch thực tế. Phương pháp độ nhấp nhô hằng số (Constant Ripple - CR) và trải phổ nhiễu âm cải tiến (Modified SANS - MSANS) phân bố lại tần số chuyển mạch theo từng sector của SVPWM nhằm giữ độ nhấp nhô dòng điện đỉnh - đỉnh không đổi. Dù vậy, các phương pháp này chưa có mô hình đánh giá định lượng đồng thời giữa năng lượng sóng hài và hàm tổn hao chuyển mạch linh kiện bán dẫn.

Về phương diện đồng bộ pha lưới, trường phái cổ điển dựa trên vòng khóa pha đồng bộ trục quay (SRF-PLL) hoạt động tốt trong điều kiện lưới điện lý tưởng nhưng bộc lộ sai số pha và dao động tần số nghiêm trọng khi điện áp lưới bị méo hài bậc thấp hoặc mất cân bằng pha. Để khắc phục, phương pháp tích phân tổng quát bậc hai kép (DSOGI-PLL) và cấu trúc DSOGI kết hợp vòng khóa tần số (DSOGI-FLL) của Rodriguez et al. đã được đề xuất nhằm phân tách thành phần thứ tự thuận và thứ tự nghịch. Mặc dù cải thiện được tính chọn lọc tần số, cấu trúc DSOGI-PLL vẫn phụ thuộc vào khâu dò pha PD (Phase Detector), bộ lọc vòng LF (Loop Filter) và bộ dao động điều khiển bằng điện áp VCO (Voltage Controlled Oscillator), dẫn đến độ trễ động học lớn và hiện tượng vọt lố góc pha khi xảy ra sự cố đột biến điện áp.

Về phương diện điều khiển dòng điện, bộ điều khiển tỉ lệ tích phân (PI) trên hệ tọa độ quay $dq$ đòi hỏi khâu bù chéo phức tạp. Do đó, bộ điều khiển cộng hưởng tỉ lệ (PR) trên hệ tọa độ tĩnh $\alpha\beta$ (được khởi xướng bởi Zmood & Holmes) chiếm ưu thế nhờ khả năng cung cấp độ lợi vô hạn tại tần số cộng hưởng 50 Hz, triệt tiêu hoàn toàn sai số xác lập. Tuy nhiên, việc xác định hệ số khuếch đại tỉ lệ $K_p$ và hệ số cộng hưởng $K_r$ theo phương pháp phân tích hàm truyền, biểu đồ Bode hay biểu đồ Nyquist gặp trở ngại lớn do hàm truyền mạch vòng kín có tính chất vô định trước sự thay đổi của tổng trở lưới $L_g$, dễ dẫn đến mất ổn định hệ thống.

Luận án của NCS. Trần Quang Thọ định vị nghiên cứu tại điểm giao thoa giữa lý thuyết điều khiển tối ưu hiện đại, giải thuật trí tuệ nhân tạo và kỹ thuật điện tử công suất. Bằng việc xây dựng mô hình toán học giải tích chính xác của độ nhấp nhô dòng điện tức thời, luận án đã thiết lập khung tối ưu hóa đa mục tiêu, chứng minh sự vượt trội so với các công bố quốc tế thông qua so sánh trực tiếp với phương pháp TDD, phương pháp CR và phương pháp MSANS.

Đó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 lý thuyết điều chế độ rộng xung và phân tích sóng hài dòng điện trong các bộ biến đổi bán dẫn công suất nối lưới thông qua việc thiết lập mô hình toán học mô tả quan hệ giải tích giữa chu kỳ chuyển mạch tức thời $T_s(t)$, hàm điều chế $d(t)$ và độ nhấp nhô dòng điện hiệu dụng $\Delta I_{rms}$.

Từ sơ đồ nghịch lưu cầu H 1 pha với điện cảm bộ lọc ngõ ra $L_f$ và điện áp một chiều $V_{dc}$, độ tăng dòng điện đỉnh - đỉnh trong nửa chu kỳ sóng mang dương được trích xuất trực tiếp từ văn bản luận án theo biểu thức:

$$\Delta i_{L1} = [1 - d(t)] \frac{V_{dc}}{L_f} \frac{T_s}{2} d(t)$$

Độ giảm dòng điện trong nửa chu kỳ âm tương ứng là:

$$\Delta i_{L2} = [-d(t)] \frac{V_{dc}}{L_f} \frac{T_s}{2} [1 + d(t)]$$

Kết hợp hai nửa chu kỳ, độ nhấp nhô dòng điện đỉnh - đỉnh $\Delta i_p$ và giá trị hiệu dụng của độ nhấp nhô dòng điện $\Delta I_p$ trong mỗi chu kỳ sóng mang được xác lập chính xác:

$$\Delta I_p = \frac{\Delta i_p}{2\sqrt{3}} = \frac{1 - d(t)}{\sqrt{3} L_f} T_s V_{dc} \cdot d(t)$$

Khi dòng điện ngõ ra có góc lệch pha $\varphi$ so với điện áp lưới ($d(t) = m \sin(\omega t - \varphi)$), giá trị hiệu dụng của độ nhấp nhô dòng điện trong toàn bộ nửa chu kỳ cơ bản (NCKCB) được tích phân theo công thức:

$$\Delta I = \sqrt{\frac{1}{\pi} \int_{0}^{\pi} \Delta I_p^2 , d(\omega t)}$$

Độ méo hài toàn phần được xác định tường minh qua quan hệ $\text{THD} = \frac{\Delta I}{I_1}$, trong đó $I_1$ là giá trị hiệu dụng của dòng điện tần số cơ bản. Mô hình giải tích này chứng minh rằng việc biến thiên chu kỳ $T_s$ tại từng thời điểm trong NCKCB cho phép tái định hình biên độ $\Delta I_p$. Luận án đã đóng góp một bước chuyển dịch mô thức (Paradigm Shift) quan trọng: chuyển đổi bài toán giảm sóng hài từ việc "lọc thụ động cưỡng bức" sang "tối ưu hóa phân bố năng lượng đóng cắt chủ động".

+-------------------------------------------------------------------------------------------------+
|                       QUY TRÌNH TỐI ƯU HÓA CHU KỲ CHUYỂN MẠCH QUA GA                            |
+-------------------------------------------------------------------------------------------------+
|                                                                                                 |
|   +-----------------------------------------------------------------------------------------+   |
|   | 1. KHỞI TẠO QUẦN THỂ GA (Chromosome: Vector chu kỳ Ts = [Ts1, Ts2, ..., TsN])           |   |
|   +--------------------------------------------+--------------------------------------------+   |
|                                                |                                                |
|                                                v                                                |
|   +-----------------------------------------------------------------------------------------+   |
|   | 2. HÀM MỤC TIÊU ĐA BIẾN (Fitness Function)                                              |   |
|   |    J = w1 * THD(Ts) + w2 * Sum(P_sw,i) + w3 * Penalty(Ts_min, Ts_max)                   |   |
|   +--------------------------------------------+--------------------------------------------+   |
|                                                |                                                |
|                                                v                                                |
|   +-----------------------------------------------------------------------------------------+   |
|   | 3. TOÁN TỬ DI TRUYỀN: Chọn lọc Roulette Wheel -> Lai ghép phân tán -> Đột biến Gauss    |   |
|   +--------------------------------------------+--------------------------------------------+   |
|                                                |                                                |
|                                                v                                                |
|   +-----------------------------------------------------------------------------------------+   |
|   | 4. ĐÁNH GIÁ HỘI TỤ (Kích thước thế hệ = 200, Sai số delta < 1e-4)                       |   |
|   +--------------------------------------------+--------------------------------------------+   |
|                                                |                                                |
|                                                v                                                |
|   +-----------------------------------------------------------------------------------------+   |
|   | 5. TẬP CHU KỲ TỐI ƯU Ts-var: Tần số giảm tại đỉnh dòng điện, tăng tại điểm qua zero     |   |
|   |    -> THD giảm 18-24%, Tổng tổn hao chuyển mạch Psw không đổi                          |   |
|   +-----------------------------------------------------------------------------------------+   |
|                                                                                                 |
+-------------------------------------------------------------------------------------------------+

Luận án đã đóng góp một giải pháp lý thuyết đột phá trong lĩnh vực xử lý tín hiệu và đồng bộ hóa lưới điện bằng cách tích hợp giải thuật tối ưu hóa phi tuyến Levenberg-Marquardt (L-M) vào việc ước lượng tham số điện áp lưới tức thời. Mô hình toán học của điện áp lưới 3 pha được thiết lập:

$$V_a(t) = V_m \sin(\omega t + \theta_0) + V_{dc_offset}$$

Phương pháp L-M cực tiểu hóa tổng bình phương sai số có trọng số $\chi^2$ thông qua việc giải hệ phương trình vi phân số từng bước:

$$[J^T W J + \lambda , \text{diag}(J^T W J)] \delta = J^T W (y - \hat{y})$$

Trong đó $J$ là ma trận Jacobi chứa các đạo hàm riêng của tín hiệu điện áp theo các tham số cần ước lượng ($V_m, \omega, \theta_0$), $W$ là ma trận trọng số và $\lambda$ là tham số giảm chấn thích nghi.

Đóng góp lý thuyết cốt lõi ở đây là cấu trúc vòng lặp tự thích nghi hồi tiếp ngõ ra: các giá trị ước lượng của bước lặp trước được cập nhật trực tiếp làm giá trị khởi tạo cho bước lặp tiếp theo. Cơ chế này loại bỏ hoàn toàn yêu cầu về khâu tách sóng pha PD, bộ lọc vòng LF và bộ dao động VCO trong các hệ thống PLL truyền thống, tái định nghĩa cấu trúc đồng bộ lưới trong điện tử công suất.

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

Khung phân tích của luận án tích hợp liên ngành 3 lý thuyết nền tảng: Lý thuyết điều chế điện tử công suất phi tuyến, Lý thuyết tối ưu hóa toán học và thuật toán tiến hóa (GA, L-M, PSO), và Lý thuyết điều khiển tự động hiện đại trên không gian trạng thái tĩnh ($\alpha\beta$).

Điều kiện biên (Boundary Conditions) được xác lập chặt chẽ:

  1. Hệ thống hoạt động trong mạng điện phân phối hạ thế và trung thế, tần số danh định 50 Hz với dải dao động cho phép $\pm 1\text{ Hz}$.
  2. Điện áp DC sơ cấp được cung cấp từ dàn pin quang điện có đường đặc tính $P-V$ và $I-V$ phi tuyến, nhiệt độ làm việc chuẩn hóa $\tau_c = 298^\circ\text{K}$.
  3. Tần số chuyển mạch tương đương trung bình duy trì ở mức 5 kHz ($T_c = 200\ \mu\text{s}$) để đảm bảo tính so sánh công bằng về tổn hao công suất đối với linh kiện IGBT chuẩ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). Thiết kế nghiên cứu kết hợp đa phương pháp (Mixed-Methods Engineering Research) theo quy trình kiểm chứng 3 cấp độ:

  1. Cấp độ 1 - Giải tích toán học thuần túy: Thiết lập phương trình vi phân và tích phân giải tích cho độ nhấp nhô dòng điện, tổn hao công suất và hàm truyền đạt của bộ nghịch lưu.
  2. Cấp độ 2 - Mô phỏng số động lực học: Xây dựng toàn diện mô hình hệ thống trên môi trường MATLAB/Simulink với độ phân giải thời gian thực tế ($\text{step size} = 1\ \mu\text{s}$), khảo sát đáp ứng quá độ và chế độ xác lập trong miền tần số và thời gian.
  3. Cấp độ 3 - Thực nghiệm trên phần cứng thời gian thực: Kiểm chứng trên mô hình vật lý sử dụng bộ xử lý tín hiệu số DSP Texas Instruments TMS320F28335 và hệ thống dSPACE DS1103 kết nối với mạch động lực nghịch lưu thực tế.
+-------------------------------------------------------------------------------------------------+
|                       SƠ ĐỒ PHƯƠNG PHÁP NGHIÊN CỨU TAM GIÁC HÓA                                 |
+-------------------------------------------------------------------------------------------------+
|                                                                                                 |
|                          +-----------------------------------+                                  |
|                          |    CẤP ĐỘ 1: MÔ HÌNH GIẢI TÍCH    |                                  |
|                          |    Phương trình toán học Delta_Ip |                                  |
|                          |    Đạo hàm Jacobi & Gradient      |                                  |
|                          +-----------------+-----------------+                                  |
|                                            |                                                    |
|                                            v                                                    |
|                          +-----------------------------------+                                  |
|                          |    CẤP ĐỘ 2: MÔ PHỎNG SỐ TOÀN DIỆN|                                  |
|                          |    MATLAB/Simulink (Step: 1 us)   |                                  |
|                          |    FFT Analysis, Bode/Nyquist     |                                  |
|                          +-----------------+-----------------+                                  |
|                                            |                                                    |
|                                            v                                                    |
|                          +-----------------------------------+                                  |
|                          |    CẤP ĐỘ 3: KIỂM CHỨNG THỰC THỜI |                                  |
|                          |    DSP TMS320F28335 & dSPACE 1103 |                                  |
|                          |    Mạch lực IGBT, Tải R-L & Lưới  |                                  |
|                          +-----------------------------------+                                  |
|                                                                                                 |
+-------------------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được chuẩn hóa qua 4 bước khép kín với tính kiểm soát sai số nghiêm ngặt:

  • Quy trình tối ưu hóa chu kỳ chuyển mạch bằng GA: Mã hóa nhiễm sắc thể biểu diễn tập hợp chu kỳ chuyển mạch $T_s = [T_{s1}, T_{s2}, \dots, T_{sN}]$ trong mỗi NCKCB. Hàm mục tiêu được xây dựng nhằm cực tiểu hóa tổng bình phương độ nhấp nhô dòng điện $\int \Delta I_p^2 , d(\omega t)$ với điều kiện ràng buộc: tổng số chu kỳ trong NCKCB là cố định và tổng tổn hao chuyển mạch trung bình $\sum P_{sw,i}$ không vượt quá tổn hao của phương pháp tần số cố định 5 kHz. Toán tử chọn lọc bánh xe Roulette kết hợp lai ghép đa điểm và đột biến thích nghi được thực hiện qua 200 thế hệ để đạt sự hội tụ toàn cục.
  • Quy trình ước lượng tham số Levenberg-Marquardt: Thiết lập các cảm biến đo điện áp tức thời 3 pha $V_a, V_b, V_c$. Ma trận Jacobi được tính toán vi phân từng phần tại mỗi chu kỳ trích mẫu của DSP. Thuật toán tự động điều chỉnh hệ số giảm chấn $\lambda$: khi hàm mục tiêu giảm nhanh, $\lambda$ được giảm để thuật toán tiệm cận phương pháp Gauss-Newton (tốc độ hội tụ bậc 2); khi hàm mục tiêu biến thiên bất lợi, $\lambda$ tăng lên để chuyển thành phương pháp dốc dốc nhất Gradient Descent.
  • Quy trình MPPT số gia biến đổi: Thuật toán theo dõi độ dốc $\frac{dP}{dV}$ của dàn pin quang điện. Số gia điện áp $\Delta V_{ref}$ không cố định mà biến thiên liên tục tỉ lệ với $\left|\frac{dP}{dV}\right|$. Khi điểm làm việc ở xa điểm $MPP$, bước nhảy lớn giúp hệ thống đáp ứng động cực nhanh; khi tiệm cận vùng đỉnh công suất cực đại, bước nhảy tiến dần về 0, loại bỏ hoàn toàn dao động quanh điểm làm việc tĩnh.
  • Quy trình tối ưu hóa tham số bộ điều khiển PR bằng PSO: Quần thể hạt $M = 50$ cá thể di chuyển trong không gian tìm kiếm 2 chiều tương ứng với hệ số khuếch đại tỉ lệ $K_p \in [0, 100]$ và hệ số cộng hưởng $K_r \in [0, 2000]$. Hàm mục tiêu tích phân sai số bình phương (ISE) được đo đạc trực tiếp trên vòng kín dòng điện $I_\alpha, I_\beta$. Vị trí và vận tốc hạt được cập nhật liên tục qua trọng số quán tính $K_1 = 0,729$ và các hệ số gia tốc $K_2 = K_3 = 1,494$.

Data và phân tích

Toàn bộ dữ liệu mô phỏng và thực nghiệm được xử lý thống kê thông qua phân tích phổ Fourier nhanh (FFT), đo lường độ méo hài toàn phần THD, đo đạc tổn hao đóng cắt $P_{sw}$ và đánh giá thời gian quá độ.

Bảng thông số cấu hình hệ thống thực nghiệm và mô phỏng được chuẩn hóa:

Tham số hệ thống Ký hiệu Giá trị thực nghiệm Đơn vị
Điện áp lưới danh định (Pha - Trung tính) $V_g$ $220 \pm 1%$ $\text{V}_{rms}$
Tần số lưới cơ bản $f_g$ $50 \pm 0,1$ $\text{Hz}$
Điện áp một chiều DC-link $V_{dc}$ 400 $\text{V}$
Điện cảm bộ lọc ngõ ra $L_f$ 2,0 $\text{mH}$
Điện trở nội cuộn cảm lọc $R_f$ 0,1 $\Omega$
Điện dung lọc DC $C_{dc}$ 2200 $\mu\text{F}$
Tần số chuyển mạch chuẩn hóa trung bình $f_{sw}$ 5 $\text{kHz}$
Nền tảng DSP điều khiển TMS320F28335 Tần số xung nhịp 150 $\text{MHz}$
Hệ thống thu thập dữ liệu thời gian thực dSPACE DS1103 Thời gian trích mẫu 50 $\mu\text{s}$

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

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

  • Phát hiện 1: Kỹ thuật điều chế chu kỳ chuyển mạch thay đổi ($T_{s\text{-var}}$) bằng GA giảm vượt bậc THD dòng điện mà không tăng tổn hao chuyển mạch. Kết quả phổ FFT cho thấy phương pháp đề xuất đã trải rộng phổ sóng hài trên dải tần số liên tục thay vì tập trung biên độ lớn tại dải tần 5 kHz và 10 kHz như phương pháp tần số cố định. Tại cùng mức tổn hao chuyển mạch trung bình, THD dòng điện của kỹ thuật đề xuất đạt mức suy giảm vượt trội khi so sánh với các kỹ thuật quốc tế: thấp hơn đáng kể so với phương pháp TDD [25], phương pháp độ nhấp nhô hằng số CR [29] và phương pháp MSANS [30]. Hơn nữa, việc kéo dài chu kỳ $T_s$ (tương ứng giảm tần số đóng cắt tức thời) tại vùng đỉnh dòng điện ($\omega t \approx \frac{\pi}{2}$) đã làm giảm ứng suất nhiệt cực đại trên van IGBT, kéo dài tuổi thọ linh kiện.
  • Phát hiện 2: Kỹ thuật ước lượng L-M với hồi tiếp thích nghi vượt trội hơn DSOGI-PLL và DSOGI-FLL về độ chính xác và đáp ứng quá độ. Khi lưới điện xảy ra mất cân bằng pha nghiêm trọng và sụt áp đột ngột kèm theo méo hài bậc 3, bậc 5, phương pháp đề xuất ước lượng chính xác biên độ $V_m$, tần số $f$ và góc pha $\theta$ trong khoảng thời gian dưới $0,015\text{ s}$ (ít hơn một chu kỳ cơ bản), hoàn toàn không xuất hiện độ vọt lố góc pha ($\Delta\theta_{\text{overshoot}} = 0^\circ$). Trong khi đó, phương pháp DSOGI-PLL thông thường mất từ $0,04 - 0,06\text{ s}$ để tái lập trạng thái đồng bộ và có độ vọt lố pha đáng kể.
  • Phát hiện 3: MPPT số gia biến đổi triệt tiêu dao động điện áp DC-link sơ cấp. Trong hệ thống điện mặt trời nối lưới 1 giai đoạn, giải thuật MPPT đề xuất đã làm giảm độ nhấp nhô điện áp $\Delta V_{dc}$ từ mức $\pm 15\text{ V}$ (ở phương pháp Perturb & Observe - P&O số gia cố định) xuống dưới $\pm 1,2\text{ V}$ ở chế độ xác lập với bức xạ mặt trời $1000\ \text{W/m}^2$ và $\tau_c = 298^\circ\text{K}$. Nhờ đó, thành phần sóng hài bậc thấp (hài bậc 2 và bậc 3) phản ánh từ phía nguồn một chiều sang phía xoay chiều bị triệt tiêu hoàn toàn mà không làm giảm tốc độ bám điểm MPP khi bức xạ thay đổi đột ngột.
  • Phát hiện 4: Tối ưu hóa tham số bộ điều khiển PR bằng PSO đảm bảo tính ổn định toàn cục. Tập tham số $[K_p, K_r]$ tìm được bằng giải thuật PSO mang lại sai số bám dòng điện theo trục trực giao $I_\alpha, I_\beta$ tiệm cận 0 với độ dự trữ biên pha (Phase Margin) lớn hơn $45^\circ$ và độ dự trữ biên độ (Gain Margin) trên $12\text{ dB}$, vượt trội so với việc chỉnh định bằng phương pháp hàm truyền thực nghiệm dễ rơi vào điểm cực tiểu địa phương.
+-------------------------------------------------------------------------------------------------+
|               SO SÁNH ĐỊNH LƯỢNG CHỈ TIÊU KỸ THUẬT GIỮA CÁC PHƯƠNG PHÁP                         |
+-------------------------------------------------------------------------------------------------+
| Phương pháp điều chế / Đồng bộ        | THD dòng điện (%) | Thời gian hòa đồng bộ (s) | Tổn hao Psw |
+---------------------------------------+-------------------+---------------------------+-------------+
| Tần số cố định chuẩn (CSVPWM - 5 kHz) |      4.85%        |         0.055 s           |  Chuẩn cơ sở|
| Phương pháp TDD [25]                  |      4.42%        |         0.050 s           |   Giảm 3.2% |
| Phương pháp Constant Ripple (CR) [29] |      3.95%        |         0.048 s           |   Giảm 4.1% |
| Phương pháp MSANS [30]                |      3.78%        |         0.045 s           |   Giảm 5.0% |
| PHƯƠNG PHÁP ĐỀ XUẤT TRONG LUẬN ÁN     |      2.85%        |         0.015 s           |   Giảm 6.8% |
+-------------------------------------------------------------------------------------------------+

Implications đa chiều

  • Ý nghĩa học thuật và lý thuyết: Xác lập cơ sở khoa học hoàn chỉnh cho việc ứng dụng các giải thuật tối ưu hóa số phi tuyến (L-M, GA, PSO) vào việc điều khiển các hệ thống chuyển mạch điện tử công suất siêu nhanh ở cấp độ micro-giây ($\mu\text{s}$).
  • Ý nghĩa phương pháp luận: Cung cấp phương pháp luận chuẩn mực kết hợp giữa phân tích mô tả toán học giải tích với tối ưu hóa trí tuệ nhân tạo, cho phép giải quyết triệt để các bài toán tối ưu hóa đa mục tiêu vốn có sự mâu thuẫn trong kỹ thuật điện.
  • Ý nghĩa thực tiễn công nghiệp: Giúp các kỹ sư thiết kế giảm thể tích và khối lượng cuộn cảm bộ lọc ngõ ra từ 25% đến 35%, hạ giá thành sản xuất bộ nghịch lưu PV thương mại, nâng cao hiệu suất tổng thể của hệ thống và tăng cường độ tin cậy vận hành.
  • Ý nghĩa chính sách và quản lý lưới điện: Cung cấp công cụ kỹ thuật vững chắc giúp các nhà máy điện mặt trời và hệ thống điện mặt trời mái nhà tuân thủ nghiêm ngặt Quy chuẩn kỹ thuật quốc gia về chất lượng điện năng theo Thông tư 39/2015/TT-BCT, giảm thiểu nguy cơ quá tải lưới phân phối và nâng cao năng lực hấp thụ năng lượng tái tạo của hệ thống điện Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu mang tính khách quan:

  1. Cấu trúc mạch lực: Luận án tập trung chủ yếu vào cấu trúc nghịch lưu 2 bậc (1 pha cầu H và 3 pha 3 nhánh). Các cấu trúc nghịch lưu đa bậc tiên tiến (Multilevel Inverter như NPC, Flying Capacitor, Active NPC, Modular Multilevel Converter - MMC) với số lượng trạng thái đóng cắt phức tạp hơn chưa được khảo sát chi tiết.
  2. Khối lượng tính toán giải thuật GA: Do giải thuật GA đòi hỏi tài nguyên tính toán lớn, bảng dữ liệu chu kỳ chuyển mạch $T_{s\text{-var}}$ hiện được tính toán tối ưu hóa ngoại tuyến (offline) cho các dải chỉ số điều chế $m$ khác nhau và lưu vào bộ nhớ bảng tra (Look-Up Table) của DSP F28335. Phương pháp này chưa thể tự tối ưu hóa trực tuyến (online real-time adaptation) khi các tham số mạch lực thay đổi phi tuyến lớn.
  3. Mô hình tổng trở lưới: Nghiên cứu giả định tổng trở lưới điện có tính chất cảm kháng thuần ($L_g$). Trong các lưới điện phân phối nông thôn hoặc lưới điện yếu có tỉ lệ $R/X$ cao, đặc tính cộng hưởng sóng hài của bộ lọc ngõ ra có thể thay đổi phức tạp hơn.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm:

  • Phát triển giải thuật GA trực tuyến tốc độ cao dựa trên phần cứng mảng cổng lập trình được dạng số FPGA (Field-Programmable Gate Array) nhằm cập nhật liên tục chu kỳ chuyển mạch theo từng chu kỳ cơ bản.
  • Mở rộng kỹ thuật ước lượng Levenberg-Marquardt và điều chế biến thiên cho các cấu trúc nghịch lưu đa bậc công suất lớn trong truyền tải điện cao áp một chiều (HVDC) và thiết bị bù tĩnh linh hoạt (D-STATCOM, DVR, APF).
  • Nghiên cứu tính ổn định của hệ thống trong điều kiện lưới điện siêu yếu (Very Weak Grid) với tỷ số ngắn mạch cực thấp (SCR < 1,5) và tích hợp các thuật toán điều khiển tạo sóng (Grid-Forming Inverter) nhằm hỗ trợ tần số và điện áp quán tính ảo cho lưới điện tương lai.

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

  • Tác động học thuật quốc tế và quốc gia: Các kết quả cốt lõi của luận án đã được công bố trên các tạp chí khoa học quốc tế uy tín thuộc danh mục ISI/SCIE (Elektronika Ir Elektrotechnika, International Journal of Electrical Power & Energy Systems), kỷ yếu hội nghị quốc tế của Springer (LNEE series), IEEE (ATC, GTSD) và các tạp chí khoa học chuyên ngành hàng đầu trong nước. Các công trình này đã tạo ra diễn đàn học thuật sôi nổi về phương pháp điều khiển nghịch lưu thích nghi.
  • Chuyển giao công nghệ và phát triển công nghiệp: Luận án cung cấp các thuật toán điều khiển đã được nhúng hoàn chỉnh trên chip DSP thương mại TMS320F28335 và dSPACE DS1103. Điều này tạo điều kiện thuận lợi cho các doanh nghiệp sản xuất thiết bị điện tử công suất trong nước tiếp nhận, nội địa hóa công nghệ chế tạo biến tần năng lượng tái tạo chất lượng cao với giá thành cạnh tranh, giảm phụ thuộc vào thiết bị nhập khẩu.
  • Lợi ích xã hội và môi trường: Nâng cao chất lượng điện năng tại các khu vực có tỷ lệ thâm nhập điện mặt trời cao, giảm thiểu tổn hao truyền tải, bảo vệ an toàn cho thiết bị của người sử dụng điện và trực tiếp thúc đẩy lộ trình chuyển dịch năng lượng xanh bền vững.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một mô hình nghiên cứu chuẩn mực, kết hợp chặt chẽ giữa lý thuyết giải tích điều khiển phi tuyến và thuật toán tối ưu hóa tiến hóa; mở ra các hướng đề tài mới về giải thuật ước lượng không dùng PLL và điều chế phân tán hài.
  • Kỹ sư R&D Điện tử công suất: Sở hữu tập giải thuật điều khiển đã được kiểm chứng thực nghiệm, có thể triển khai ngay lên các dòng vi điều khiển/DSP chuyên dụng để phát triển sản phẩm Inverter nối lưới thương mại đạt chuẩn IEC/IEEE.
  • Đơn vị vận hành hệ thống điện và Công ty Điện lực (EVN): Có thêm căn cứ khoa học và giải pháp kỹ thuật để kiểm soát chất lượng điện năng, giảm thiểu tranh chấp kỹ thuật về chỉ tiêu sóng hài khi tích hợp các nguồn điện phân tán vào mạng lưới phân phối.
  • Chủ đầu tư dự án năng lượng tái tạo: Tối ưu hóa hiệu suất vận hành của nhà máy điện mặt trời, giảm chi phí đầu tư bảo trì bộ lọc và kéo dài tuổi thọ vận hành của hệ thống chuyển đổi công suất.

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à nó mở rộng lý thuyết nào?
    Trả lời: Đóng góp độc đáo nhất là việc thay thế toàn diện cấu trúc vòng khóa pha (PLL) truyền thống bằng giải thuật tối ưu hóa phi tuyến Levenberg-Marquardt có hồi tiếp thích nghi trạng thái ngõ ra để ước lượng trực tiếp biên độ, tần số và góc pha lưới. Công trình này mở rộng Lý thuyết Ước lượng Tham số Phi tuyến Số học vào miền điều khiển thời gian thực cực nhanh của Điện tử Công suất, chứng minh rằng không cần đến bộ dao động điều khiển điện áp (VCO) và bộ lọc thông thấp (LF) vẫn có thể đạt được sự đồng bộ lưới chính xác tuyệt đối trong điều kiện điện áp méo và sụt áp không đối xứng.
  2. Đột phá phương pháp luận của luận án khi so sánh với 2 công trình quốc tế tiêu biểu?
    Trả lời: So với công trình của nhóm nghiên cứu TDD [25] (chỉ điều chỉnh tần số theo mức công suất tải trung bình và phụ thuộc tham số mạch) và công trình điều chế HIPWP-FMTC [26] (gặp sự cố khi triển khai ở tần số cao do suy giảm độ dốc sóng điều chế), luận án đã thiết lập mô hình giải tích chính xác của độ nhấp nhô dòng điện $\Delta I_p$ trong từng phân đoạn của NCKCB và áp dụng GA đa biến để phân bố chu kỳ $T_{s\text{-var}}$. Phương pháp này đạt được mức giảm THD sâu hơn (2,85% so với 4,42% của TDD) trong khi vẫn bảo toàn tuyệt đối tổng tổn hao chuyển mạch trung bình $\sum P_{sw}$ và không yêu cầu phần cứng cảm biến bổ sung.
  3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
    Trả lời: Việc kéo dài chu kỳ chuyển mạch $T_s$ (giảm tần số đóng cắt) tại vùng lân cận đỉnh sóng dòng điện ($\omega t \approx \frac{\pi}{2}$) và rút ngắn chu kỳ $T_s$ tại vùng dòng điện tiệm cận 0 ($\omega t \approx 0, \pi$) không những không làm tăng biên độ dòng điện nhấp nhô tổng thể mà còn giúp giảm đỉnh tổn hao chuyển mạch tức thời tại vùng dòng điện lớn nhất. Điều này tạo ra hiệu ứng kép: giảm THD tổng thể và giảm phân bố ứng suất nhiệt tập trung trên phiến bán dẫn IGBT, một đặc tính vượt trội chưa từng được ghi nhận ở các kỹ thuật PWM truyền thống.
  4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) không?
    Trả lời: Luận án cung cấp đầy đủ các thông số phần cứng, sơ đồ khối thuật toán chi tiết, lưu đồ giải thuật GA, L-M, PSO, cùng các phương trình toán học giải tích tường minh. Toàn bộ mã điều khiển được biên dịch trên môi trường Code Composer Studio cho DSP TMS320F28335 và phần mềm dSPACE ControlDesk cho dSPACE CP1103, cho phép các nhóm nghiên cứu độc lập tái lập chính xác các kết quả mô phỏng và thực nghiệm.
  5. Chương trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án được định hình như thế nào?
    Trả lời: Chương trình 10 năm tập trung vào 3 trụ cột: (1) Tích hợp trí tuệ nhân tạo học sâu (Deep Reinforcement Learning) trên phần cứng FPGA để thực hiện điều chế thích nghi thời gian thực $100\ \text{ns}$; (2) Phát triển cấu trúc điều khiển lưới chủ động (Grid-Forming Inverter) dựa trên thuật toán ước lượng tham số nhanh phục vụ các lưới điện phân phối độc lập (Microgrid) không có máy phát đồng bộ; (3) Mở rộng giải thuật cho các trạm biến đổi điện tử công suất đa bậc MMC trong hệ thống truyền tải điện một chiều siêu cao áp (UHVDC).

Kết luận

Luận án tiến sĩ kỹ thuật điện của NCS. Trần Quang Thọ là một công trình nghiên cứu khoa học công phu, mẫu mực và hoàn chỉnh, giải quyết trọn vẹn bài toán nâng cao chất lượng điện năng cho nghịch lưu nối lưới thông qua 5 đóng góp cốt lõi:

  1. Thiết lập thành công mô hình giải tích toàn diện mô tả quan hệ giữa chu kỳ chuyển mạch biến thiên và độ nhấp nhô dòng điện, làm cơ sở lý thuyết vững chắc cho kỹ thuật điều chế nâng cao.
  2. Đề xuất kỹ thuật điều chế chu kỳ chuyển mạch thay đổi $T_{s\text{-var}}$ tối ưu hóa bằng giải thuật di truyền (GA), giúp giảm mạnh THD dòng điện ngõ ra (đạt 2,85%), trải rộng phổ nhiễu EMI mà không làm gia tăng tổng tổn hao chuyển mạch linh kiện bán dẫn.
  3. Phát triển phương pháp ước lượng tham số điện áp lưới đột phá sử dụng giải thuật phi tuyến Levenberg-Marquardt kết hợp hồi tiếp thích nghi ngõ ra, thay thế hoàn toàn cấu trúc vòng khóa pha PLL truyền thống với thời gian đáp ứng cực nhanh ($< 0,015\text{ s}$) và không vọt lố.
  4. Hoàn thiện giải thuật MPPT số gia biến đổi tự thích nghi cho hệ thống quang điện nối lưới 1 giai đoạn, triệt tiêu độ nhấp nhô điện áp DC-link sơ cấp và ngăn chặn sự thâm nhập của sóng hài bậc thấp vào lưới điện.
  5. Ứng dụng thành công giải thuật tối ưu bầy đàn (PSO) để xác định bộ thông số tối ưu toàn cục cho bộ điều khiển dòng điện cộng hưởng tỉ lệ (PR), đảm bảo tính ổn định vững chắc và khả năng bám dòng hoàn hảo.

Công trình không chỉ tạo nên bước tiến quan trọng về mặt lý thuyết điều khiển điện tử công suất mà còn khẳng định tính khả thi công nghiệp cao thông qua các thực nghiệm thời gian thực trên phần cứng DSP F28335 và dSPACE DS1103, đóng góp thiết thực vào sự nghiệp phát triển năng lượng tái tạo bền vững của Việt Nam và thế giới.