Tổng quan về luận án
Sự chuyển dịch cơ cấu năng lượng toàn cầu hướng tới mục tiêu trung hòa carbon đã thúc đẩy sự bùng nổ của các nguồn năng lượng tái tạo phân tán, đặc biệt là điện mặt trời quang điện (PV) và điện gió (WT). Tại Việt Nam, quy mô công suất lắp đặt nguồn năng lượng tái tạo đã có bước phát triển nhảy vọt: tính đến cuối năm 2020, điện mặt trời áp mái đạt 7.760 MW (chiếm 10,1% tổng công suất hệ thống), các trang trại điện mặt trời đạt 8.732 MW (chiếm 11,39%), và đến tháng 11/2021 điện gió đạt 3.980 MW (chiếm 5,19%). Tổng công suất lắp đặt điện mặt trời năm 2021 của Việt Nam cán mốc 10.909 MW, đứng thứ 3 toàn cầu. Tuy nhiên, tính chất bất định, biến thiên ngẫu nhiên và suy giảm quán tính quay cơ điện từ các nguồn phân tán này đặt ra thách thức nghiêm trọng: hiện tượng quá điện áp cục bộ tại các nút tải, quá tải đường dây truyền tải/phân phối, trào lưu công suất ngược (reverse power flow) về máy biến áp hạ áp 22/0,4 kV không có bộ điều áp dưới tải (OLTC), cùng rủi ro sa thải công suất nguồn sạch gây lãng phí kinh tế.
Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: phần lớn các nghiên cứu trong và ngoài nước tập trung giải bài toán trào lưu công suất tối ưu (Optimal Power Flow - OPF) hoặc tối ưu hóa vị trí - dung lượng nguồn phân tán trên lưới truyền tải trung/cao áp 3 pha đối xứng (như các mô hình chuẩn IEEE-30, IEEE-33 nút), mà bỏ qua đặc thù của lưới điện hạ áp 0,4 kV phân phối 1 pha cấp cho phụ tải dân dụng. Quy định pháp lý hiện hành tại Việt Nam (Thông tư số 39/2015/TT-BCT của Bộ Công Thương) áp dụng mức trần cứng nhắc, khống chế công suất nguồn phân tán 3 pha không vượt quá 30% công suất định mức của trạm biến áp hạ áp (ví dụ tối đa 96 kVA cho trạm 320 kVA), gây rào cản lớn cho việc giải phóng tiềm năng năng lượng sạch. Luận án tiến sĩ chuyên ngành Kỹ thuật Điều khiển và Tự động hóa (Mã số: 952 02 16) của nghiên cứu sinh Đặng Thành Trung tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam đã giải quyết triệt để bài toán này.
Các câu hỏi nghiên cứu và giả thuyết trọng tâm bao gồm:
- Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xác định chính xác tổng công suất phát cực đại của nguồn PV và WT trên từng pha của lưới điện phân phối hạ áp 0,4 kV mà vẫn đảm bảo điện áp nút $U_i \in [209\text{V}, 231\text{V}]$ (sai lệch $\pm 5%$) và dòng điện nhánh không vượt quá $I_{\max}$ khi phụ tải biến thiên từ cực tiểu (10%) đến cực đại (100%)?
- Câu hỏi nghiên cứu 2 (RQ2): Mô hình điều phối tích hợp giữa lưới phân phối hạ áp và lưới truyền tải trung/cao áp cần được thiết lập ra sao khi tính đến khả năng giải tỏa công suất ngược, chi phí kinh tế phát điện của các tổ máy truyền thống và hệ số có lợi vị trí ($\alpha_i$)?
- Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để phối hợp mạng nơ-ron nhân tạo (ANN) dự báo đa thông số chuỗi thời gian ngắn hạn (công suất PV, WT, phụ tải và giá điện) với thuật toán tối ưu hóa điều khiển phụ tải linh hoạt (Demand-Side Management) và trạm sạc xe điện (EV/BESS) nhằm nâng cao độ tin cậy cung cấp điện (SAIDI, LOLE)?
Giả thuyết nghiên cứu (Hypotheses):
- $H_1$: Tối ưu hóa phân bố không gian và công suất nguồn phân tán theo từng nút 1 pha cho phép nâng công suất tiếp nhận của lưới hạ áp vượt xa ngưỡng giới hạn 30% định mức máy biến áp mà không vi phạm bất kỳ ràng buộc kỹ thuật nào.
- $H_2$: Tích hợp bộ lưu trữ năng lượng từ pin xe điện (EV) kết hợp giàn pin PV nối lưới dưới sự điều khiển tối ưu sẽ giảm thiểu đáng kể thời gian mất điện trung bình hàng năm (SAIDI) của lưới điện khi vận hành ở chế độ tách đảo sự cố (islanding mode).
Phạm vi nghiên cứu bao quát từ mô hình toán giải tích 1 pha lưới hạ áp 0,4 kV (trạm biến áp 320 kVA, bán kính cấp điện 400 m, dây dẫn A-120 đến A-35) đến các cấu hình lưới điện 3 pha thử nghiệm 5 nút, 23 nút, IEEE-24 nút và IEEE-30 nút, kết hợp dữ liệu khí tượng và thị trường điện cạnh tranh với độ phân giải 30 phút/mẫu.
Literature Review và Positioning
Tổng quan y văn quốc tế và trong nước cho thấy ba luồng nghiên cứu (research streams) chính:
Luồng 1 - Tối ưu hóa vị trí và công suất nguồn phân tán (Distributed Generation Siting & Sizing): Các nghiên cứu của T. T. Nguyen (2018) [1] và tác giả trong nước [2, 3] sử dụng thuật toán Di truyền (GA) để tối ưu tổn thất công suất tác dụng $\Delta P$ và độ lệch điện áp trên lưới 22 kV và lưới IEEE-33 nút. Tuy nhiên, các tác giả đều thiết lập giả định lý tưởng hóa: phụ tải không biến thiên đột biến và lưới điện 3 pha hoàn toàn đối xứng. Trên trường quốc tế, các công trình của Moradi & Abedini (2012) [20], Yammani et al. (2016) [21], và Viral & Khatod (2015) [22] tích hợp giải thuật tối ưu hóa bầy đàn (PSO) hoặc mô phỏng Monte Carlo (MCS) để tính toán chi phí vòng đời và độ tin cậy. Dẫu vậy, những mô hình này chỉ áp dụng cho cấp trung áp, thiếu vắng mô hình trào lưu công suất giải tích cho lưới phân phối hạ áp 1 pha có dây trung tính - nơi trực tiếp gánh chịu dao động điện áp do nguồn PV áp mái gây ra.
Luồng 2 - Dự báo công suất phát tái tạo và phụ tải chuỗi thời gian: Nhóm nghiên cứu Mellit & Kalogirou (2008) [24], Pedro & Coimbra (2012) [29], và Bacher et al. (2009) [30] đã thực hiện các đánh giá toàn diện giữa các phương pháp: mạng nơ-ron nhân tạo (ANN), mô hình chuỗi thời gian (AR, ARMA, ARIMA) và mô hình dự báo thời tiết số trị (NWP). Kết luận thực nghiệm đồng thuận rằng ANN chiếm ưu thế tuyệt đối về độ chính xác phi tuyến. Tiếp đó, Shi et al. (2012) [26] và Chen et al. (2011) [27] đã áp dụng mạng lan truyền ngược (BPNN) và phân cụm theo ngày thời tiết (quang, mây, mưa) với chu kỳ lấy mẫu 30 phút. Đối với điện gió, Catalao et al. (2011) [36] và Foley et al. (2012) [37] chứng minh mạng nơ-ron kết hợp tối ưu tham số cho sai số MAPE nhỏ nhất. Dù vậy, các nghiên cứu dự báo này chủ yếu đứng độc lập như một công cụ giám sát, chưa được cấu trúc thành khối đầu vào trực tiếp phục vụ bộ điều khiển tối ưu phụ tải thời gian thực.
Luồng 3 - Điều khiển tối ưu Microgrid, phụ tải linh hoạt và lưu trữ năng lượng: Các nghiên cứu của Ha & Saha (2014) [56], Marra et al. (2014) [58], và Al-Saadi et al. (2018) [61] phân tích hệ thống lưu trữ BESS, sạc xe điện (EV) và máy phát diesel trong lưới điện nhỏ một chiều (DC microgrid) hoặc xoay chiều tần số cao. Nhóm tác giả Zhang et al. (2017) [63] và Sedghi et al. (2015) [65] sử dụng GA và PSO để giải bài toán đa mục tiêu kinh tế - kỹ thuật. Tuy nhiên, các công trình quốc tế thường dựa trên hạ tầng tiên tiến có sẵn bộ biến tần thông minh (Smart Inverters) hỗ trợ điều khiển công suất phản kháng $Q$ và trạm biến áp tự động điều áp dưới tải (OLTC) - điều kiện hoàn toàn không tương thích với bối cảnh lưới điện hạ áp tại các nước đang phát triển như Việt Nam.
So sánh đối sánh quốc tế:
- Trường hợp nước Đức [94]: Tính đến năm 2015, Đức lắp đặt gần 40 GW điện mặt trời, trong đó 22 GW đấu nối trực tiếp vào lưới hạ áp (90% là hệ thống áp mái $< 30\text{ kW}$). Đức giải quyết bài toán dung lượng tiếp nhận nhờ tiêu chuẩn kỹ thuật nghiêm ngặt về dây dẫn tiết diện lớn, trạm biến áp phân phối thông minh có OLTC và biến tần tham gia ổn định điện áp theo đường cong $Q(U)$.
- Trường hợp đảo Jeju (Hàn Quốc) [41-42]: Hệ thống điện cô lập với tỷ trọng điện gió cao áp dụng mô hình mạng lan truyền ngược BPNN dự báo công suất để điều độ các tổ máy nhiệt điện và pin tích năng BESS quy mô lớn.
Định vị học thuật của luận án: Luận án đã lấp đầy khoảng trống y văn bằng cách xây dựng một khung giải pháp tổng thể, khép kín: thiết lập mô hình tính toán giải tích khả năng tiếp nhận công suất cực đại trên từng pha lưới hạ áp 0,4 kV kết hợp với bài toán tối ưu trào lưu công suất OPF đa cấp (lưới cao - trung - hạ áp); đồng thời tích hợp mạng ANN huấn luyện bằng thuật toán Levenberg-Marquardt dự báo 30 phút/lần với cơ chế điều khiển tối ưu đóng/cắt phụ tải và phân bổ xả điện trạm sạc EV hỗ trợ ổn định điện áp và nâng cao độ tin cậy lưới điện.
Đó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ú hệ thống lý thuyết điều khiển tối ưu trong kỹ thuật hệ thống điện:
- Mở rộng lý thuyết trào lưu công suất tối ưu (Optimal Power Flow - OPF): Phát triển khung toán học OPF truyền thống của Carpentier (1962) và lý thuyết điểm trong (Interior Point / Barrier Method) bằng cách thiết lập cấu trúc ràng buộc phân tầng kép. Khung toán học mới không chỉ tối thiểu hóa hàm chi phí nhiên liệu của các nhà máy truyền thống $\min \sum (a_i P_i^2 + b_i P_i + c_i)$ mà còn lồng ghép ràng buộc biên điện áp 1 pha $U_{i,\min} \le U_{i,\text{phase}} \le U_{i,\max}$ và hệ số có lợi vị trí $\alpha_i$ phản ánh đặc tính không gian địa lý của nguồn tái tạo.
- Phát triển lý thuyết thích nghi bầy đàn (PSO Dynamic Boundary Correction): Mở rộng nền tảng thuật toán bầy đàn của Kennedy & Eberhart (1995) [29] bằng cơ chế hiệu chỉnh vector vận tốc triệt tiêu ($v_d = 0$) kết hợp phản xạ biên giới hạn khi hạt vi phạm không gian tìm kiếm $[x_{\min}, x_{\max}]^D$. Cải tiến này khắc phục triệt để hiện tượng "bị bẫy tại biên" (boundary lock-in) và hiện tượng hội tụ sớm (premature convergence) trong không gian biến số tối ưu rời rạc của bài toán phân bổ nguồn phân tán.
- Xác lập mô hình toán học giải tích cho trào lưu công suất ngược 1 pha: Xây dựng mô hình quan hệ phi tuyến giữa phụ tải cực tiểu $S_{L,\min}$, tổng trở đường dây $Z = R + jX$, điện áp nút nguồn $U_0$ và tổng công suất phát ngược cực đại $S_{PV} + S_{WT}$, chứng minh bằng giải tích điều kiện giới hạn ổn định điện áp cho các phụ tải dân dụng đơn pha.
Khung phân tích độc đáo
Khung phân tích của đề tài tích hợp ba trụ cột lý thuyết: (1) Lý thuyết mạch điện phi tuyến và trào lưu công suất nâng cao, (2) Lý thuyết trí tuệ tính toán (Computational Intelligence) bao gồm giải thuật di truyền thích nghi và mạng nơ-ron nơ-ron học sâu, và (3) Lý thuyết xác suất độ tin cậy hệ thống điện (Reliability Engineering).
Định nghĩa các khái niệm cốt lõi:
- Khả năng tiếp nhận công suất phân tán (Hosting Capacity): Ngưỡng công suất phát cực đại $\max \sum P_{PV_WT,i}$ mà lưới điện hạ áp 1 pha có thể hấp thụ tại chế độ tải 10% công suất định mức mà không làm điện áp bất kỳ nút nào vượt ngưỡng 231 V hoặc dòng điện nhánh vượt dòng định mức dây dẫn.
- Hệ số có lợi vị trí ($\alpha_i$): Trọng số không gian định lượng tiềm năng tài nguyên bức xạ/tốc độ gió và chi phí giải phóng mặt bằng tại nút $i$, cho phép thuật toán tối ưu ưu tiên phân bổ dung lượng tại các nút có hiệu suất kinh tế - kỹ thuật cao nhất.
- Phụ tải dịch chuyển được thời gian (Shiftable Demand Load - LC): Các phụ tải có thể thay đổi chu kỳ đóng/cắt thông qua tín hiệu điều khiển tự động mà không làm thay đổi tổng điện năng tiêu thụ trong ngày (ví dụ: máy bơm nước nông nghiệp, bình nước nóng, chu trình sạc xe điện).
Điều kiện biên (Boundary Conditions): Khung phân tích được xác lập trên điều kiện vận hành chuẩn của lưới phân phối Việt Nam: tần số định mức 50 Hz, điện áp danh định 220 V pha - trung tính, giới hạn sai lệch điện áp tĩnh $\pm 5%$, hệ số công suất phụ tải $\cos\varphi \ge 0,95$, và giới hạn nhiệt của dây dẫn trên không trần theo tiêu chuẩn TCVN.
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ủ thế giới quan thực chứng (positivism paradigm) với phương pháp định lượng diễn dịch nghiêm ngặt. Thiết kế nghiên cứu phối hợp đa tầng (multi-level hierarchical framework), kết hợp mô hình hóa giải tích toán học thuần túy với mô phỏng số học (numerical simulation) trên nền tảng MATLAB R2018a.
+---------------------------------------------------------------------------------------+
| QUY TRÌNH NGHIÊN CỨU |
+---------------------------------------------------------------------------------------+
| [TẦNG 1: THU THẬP & MÔ HÌNH HÓA THÔNG SỐ] |
+---------------------------------------------------------------------------------------+
| [TẦNG 2: TỐI ƯU HÓA CÔNG SUẤT TIẾP NHẬN BẰNG GA & PSO] |
+---------------------------------------------------------------------------------------+
| [TẦNG 3: DỰ BÁO NGẮN HẠN BẰNG MẠNG NƠ-RON LEVENBERG-MARQUARDT] |
+---------------------------------------------------------------------------------------+
| [TẦNG 4: ĐIỀU KHIỂN ĐÓNG CẮT PHỤ TẢI & TÍCH HỢP TRẠM SẠC EV] |
+---------------------------------------------------------------------------------------+
Quy trình nghiên cứu rigorous
Quy trình giải quyết bài toán tối ưu công suất nguồn phân tán 1 pha hạ áp được thực thi qua các bước giải tích chặt chẽ:
- Mô hình hóa mạch thay thế 1 pha: Tại chế độ nguồn phân tán chưa phát, trào lưu công suất trên nhánh $k-m$ và điện áp nút được tính toán giải tích:
$$\Delta S_{k-m} = \frac{P_{m}^2 + Q_{m}^2}{U_m^2} (R_{k-m} + j X_{k-m})$$
$$U_k = \sqrt{\left(U_m + \frac{P_m R_{k-m} + Q_m X_{k-m}}{U_m}\right)^2 + \left(\frac{P_m X_{k-m} - Q_m R_{k-m}}{U_m}\right)^2}$$
- Khi tích hợp nguồn PV và WT tại các nút $1 \dots 10$: Công suất nút được tái cấu trúc: $S'm = (P{L,m} - P_{PV,m} - P_{WT,m}) + j (Q_{L,m} - Q_{PV,m} - Q_{WT,m})$.
- Thiết lập bài toán tối ưu phi tuyến có ràng buộc:
$$\max f = \sum_{i=1}^n P_{PV_WT,i}$$
Ràng buộc kỹ thuật bắt buộc:
$$\begin{cases} U_{\min} \le U_i \le U_{\max} \quad (209\text{ V} \le U_i \le 231\text{ V}, \forall i \in [1, 10]) \ I_{k-m} \le I_{\max, k-m} \quad (\forall k-m \in \text{nhánh}) \end{cases}$$
- Cấu trúc giải thuật Di truyền (GA):
- Độ dài chuỗi nhị phân nút $i$: $x_i = (P_{\max,i} - P_{\min,i}) / P_{\text{step}}$, với $P_{\text{step}} = 1\text{ kW}$, thỏa mãn $x_i \le 2^{m_i} - 1$. Chiều dài nhiễm sắc thể: $m = \sum_{i=1}^{10} m_i$.
- Quần thể ban đầu: $\text{pop_size} = 50 - 100$.
- Toán tử chọn lọc: Bánh xe Roulette với xác suất chọn $p_i = f_i / \sum f$ và vị trí tích lũy $q_i = \sum_{j=1}^i p_j$.
- Toán tử lai ghép: Điểm lai ngẫu nhiên $\text{pos} \in [1, m-1]$ với xác suất lai $p_c \in [0,7, 0,9]$.
- Toán tử đột biến: Đột biến đảo bit với xác suất $p_m \in [0,01, 0,05]$.
- Xử lý ràng buộc: Áp dụng hàm phạt tĩnh (penalty function). Bất kỳ nhiễm sắc thể nào vi phạm ngưỡng điện áp hoặc dòng điện tại kịch bản tải cực tiểu 10% sẽ bị gán $f_i = 0$.
Data và phân tích
Thông số kỹ thuật mô hình lưới thử nghiệm 1 pha: Trạm biến áp 320 kVA, điện áp nút nguồn $U_0 = 220\text{ V}$ (hoặc đặt $230\text{ V}$). Lưới phân phối gồm 10 nút phụ tải với chiều dài mỗi phân đoạn 40 m:
- Nhánh 0-1: Dây nhôm A-120 ($R = 0,011,\Omega, X = 0,012,\Omega, I_{\max} = 375\text{ A}$)
- Nhánh 1-5: Dây nhôm A-95 ($R = 0,014,\Omega, X = 0,012,\Omega, I_{\max} = 320 - 375\text{ A}$)
- Nhánh 5-7: Dây nhôm A-70 ($R = 0,018,\Omega, X = 0,012,\Omega, I_{\max} = 265\text{ A}$)
- Nhánh 7-9: Dây nhôm A-50 ($R = 0,026,\Omega, X = 0,013,\Omega, I_{\max} = 215\text{ A}$)
- Nhánh 9-10: Dây nhôm A-35 ($R = 0,037,\Omega, X = 0,013,\Omega, I_{\max} = 170\text{ A}$)
Dữ liệu phụ tải nút: Công suất phụ tải tác dụng tại 10 nút lần lượt là: Nút 1 (6 kW), Nút 2 (6 kW), Nút 3 (6 kW), Nút 4 (8 kW), Nút 5 (6 kW), Nút 6 (10 kW), Nút 7 (4 kW), Nút 8 (12 kW), Nút 9 (6 kW), Nút 10 (12 kW); hệ số công suất đồng nhất $\cos\varphi = 0,95$.
Cấu trúc mạng nơ-ron Levenberg-Marquardt (LM-ANN): Mạng truyền thẳng đa lớp (Multi-Layer Perceptron) với thuật toán cập nhật trọng số:
$$w_{k+1} = w_k - \left[ J^T J + \mu I \right]^{-1} J^T e$$
Trong đó $J$ là ma trận Jacobian của vector sai số đạo hàm theo trọng số, $\mu$ là tham số kết hợp tự động điều chỉnh giữa phương pháp hạ gradient và phương pháp Gauss-Newton, đảm bảo tốc độ huấn luyện hội tụ siêu nhanh cho dữ liệu chuỗi thời gian 48 chu kỳ/ngày (30 phút/chu kỳ).
Phát hiện đột phá và implications
Những phát hiện then chốt
Kết quả mô phỏng và thực nghiệm số mang lại các bằng chứng định lượng mang tính đột phá:
+---------------------------------------------------------------------------------------+
| SO SÁNH CÔNG SUẤT TIẾP NHẬN NGUỒN TÁI TẠO HẠ ÁP |
+---------------------------------------------------------------------------------------+
| QUY ĐỊNH HIỆN HÀNH (Thông tư 39/2015/TT-BCT) |
| ████████████ 96 kVA (Trần cứng 30% công suất trạm 320 kVA) |
+---------------------------------------------------------------------------------------+
| KẾT QUẢ TỐI ƯU HÓA CỦA LUẬN ÁN (Thuật toán GA / PSO) |
| ██████████████████████████████████ 267 kW (Tăng 178,1% công suất tiếp nhận) |
+---------------------------------------------------------------------------------------+
- Phá vỡ giới hạn tiếp nhận quy chuẩn: Luận án chứng minh trạm biến áp 320 kVA có thể tiếp nhận an toàn tới 89 kW trên mỗi pha, tương đương tổng công suất 3 pha đạt 267 kW (chiếm 83,4% công suất định mức trạm). Con số này vượt trội hoàn toàn so với mức trần 96 kVA (30%) quy định tại Thông tư 39/2015/TT-BCT. Toàn bộ điện áp nút tại kịch bản khắc nghiệt nhất (tải 10%, nguồn phát 100%) đều nằm nghiêm ngặt trong dải an toàn $209\text{ V} - 231\text{ V}$.
- Phân bố không gian tối ưu tập trung tại đầu tuyến: Thuật toán GA chỉ ra cấu hình phân bổ tối ưu của pha A đạt cực đại tại 89 kW với sự phân bổ: Nút 1 (14 kW), Nút 2 (65 kW), Nút 3 (6 kW), Nút 4 (2 kW), Nút 5 (2 kW), Nút 6-10 (0 kW). Kết quả chứng minh về mặt vật lý: công suất nguồn tái tạo tập trung lớn ở các nút gần trạm nguồn (nơi tiết diện dây lớn A-120, A-95 và trở kháng thấp) cho phép giải tỏa công suất ngược qua trạm biến áp mà không gây sụt áp hoặc quá áp tích lũy về cuối tuyến.
- Vượt trội của PSO so với GA trong tối ưu hóa hệ thống điện: Khi thử nghiệm trên các hệ thống IEEE-24 và IEEE-30 nút, giải thuật PSO với cơ chế điều chỉnh vận tốc triệt tiêu đạt trạng thái hội tụ tại bước lặp thứ 35-42, nhanh hơn 38% so với GA (hội tụ tại bước lặp 60-75) và loại bỏ hoàn toàn hiện tượng dao động quanh nghiệm cực trị cục bộ.
- Hiệu quả giảm thiểu SAIDI bằng trạm sạc EV: Mô phỏng Monte Carlo trên lưới phân phối tích hợp chứng minh rằng việc khai thác dung lượng xả từ trạm sạc xe điện kết hợp điện mặt trời áp mái trong chế độ vận hành tách đảo sự cố giúp giảm chỉ số thời gian mất điện trung bình (SAIDI) từ 24,5 giờ/khách hàng/năm xuống còn 8,2 giờ/khách hàng/năm (giảm 66,5% thời gian gián đoạn cấp điện).
Implications đa chiều
- Về mặt học thuật và phương pháp luận: Cung cấp một phương pháp luận tích hợp hoàn chỉnh từ cấp độ giải tích vi mô (từng pha hạ áp) đến điều độ vĩ mô (lưới truyền tải trung/cao áp có điều độ OPF nguồn truyền thống), mở ra hướng nghiên cứu mới cho bài toán tối ưu hóa hệ thống điện đa tầng bất đối xứng.
- Về mặt thực tiễn vận hành công nghiệp: Cung cấp công cụ phần mềm tính toán tự động cho các Công ty Điện lực (PC) trực thuộc EVN để cấp phép đấu nối nguồn điện mặt trời mái nhà và điện gió phân tán một cách khoa học, thay vì áp dụng máy móc quy tắc tỷ lệ phần trăm cố định.
- Về mặt hoạch định chính sách: Cung cấp cơ sở khoa học định lượng vững chắc cho Bộ Công Thương và Cục Điều tiết Điện lực trong việc sửa đổi, cập nhật Thông tư số 39/2015/TT-BCT và xây dựng cơ chế mua bán điện trực tiếp (DPPA), thúc đẩy phát triển xe điện gắn liền với biểu giá điện theo thời gian sử dụng (TOU).
Limitations và Future Research
Mặc dù đạt được những đóng góp học thuật và ứng dụng to lớn, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:
- Giới hạn mô hình hóa pha: Luận án phân tích chi tiết trên từng pha độc lập với giả định dòng trên dây trung tính được triệt tiêu tương đối hoặc đã được cân bằng pha sơ bộ; chưa mô hình hóa liên kết hỗ cảm điện từ phi tuyến đồng thời giữa 3 pha và dây trung tính trong điều kiện mất cân bằng pha nghiêm trọng.
- Giới hạn dữ liệu chuỗi thời gian khí tượng: Bộ dữ liệu huấn luyện nơ-ron LM-ANN chủ yếu dựa trên các kịch bản thời tiết điển hình (quang, mây, mưa) tại một số vùng địa lý nhất định, chưa bao quát đầy đủ các hiện tượng thời tiết cực đoan (dông lốc đột ngột, sương mù kéo dài).
- Giới hạn kinh tế pin xe điện: Bài toán điều khiển xả trạm sạc EV chưa tính đến hàm suy giảm tuổi thọ pin (battery degradation cost) theo số chu kỳ sạc/xả sâu.
Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):
- Mở rộng mô hình trào lưu công suất giải tích 3 pha 4 dây không đối xứng hoàn toàn kết hợp điều khiển cân bằng pha tự động bằng thiết bị chuyển pha tĩnh (Static Phase Balancer).
- Phát triển các thuật toán học sâu thời gian thực (Deep Reinforcement Learning - DRL như DDPG, SAC) trong điều khiển phân tán thích nghi cho hệ thống hàng triệu thiết bị IoT và trạm sạc xe điện thông minh (V2G).
- Nghiên cứu bài toán tối ưu hóa thị trường điện ngang hàng (Peer-to-Peer - P2P Energy Trading) ứng dụng công nghệ chuỗi khối (Blockchain) trên nền tảng lưới phân phối hạ áp.
- Tích hợp công nghệ pin nhiên liệu Hydro xanh (Green Hydrogen Fuel Cell) làm nguồn tích trữ năng lượng dài hạn kết hợp với các nguồn điện gió ngoài khơi quy mô lớn.
Tác động và ảnh hưởng
- Tác động học thuật: Các bài báo khoa học xuất bản từ luận án trên các tạp chí chuyên ngành uy tín trong nước và quốc tế (thuộc danh mục Scopus/SCIE) tạo tiền đề cho việc trích dẫn trong các nghiên cứu về tích hợp năng lượng tái tạo và lưới điện thông minh (Smart Grid).
- Chuyển dịch ngành công nghiệp năng lượng: Cung cấp giải pháp kỹ thuật trực tiếp giúp giải tỏa công suất nghẽn mạch tại các "điểm nóng" năng lượng tái tạo (như khu vực miền Trung và miền Nam Việt Nam), tối ưu hóa vốn đầu tư hàng tỷ USD của các nhà đầu tư điện áp mái và điện gió.
- Ảnh hưởng chính sách quốc gia: Đóng góp luận cứ khoa học vào Quy hoạch phát triển điện lực quốc gia (Quy hoạch điện VIII), định hình lộ trình phát triển hệ sinh thái xe điện quốc gia gắn liền với hạ tầng lưới phân phối tự động hóa.
- Lợi ích môi trường và xã hội: Việc nâng cao khả năng tiếp nhận năng lượng sạch thêm 178,1% trên mỗi trạm phân phối góp phần cắt giảm hàng triệu tấn phát thải khí nhà kính $CO_2$, thúc đẩy 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â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ền tảng nào?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Trào lưu Công suất Tối ưu (OPF) truyền thống của Carpentier bằng cách thiết lập mô hình toán tối ưu hóa phân tầng liên kết: giải mã không gian trào lưu công suất phi tuyến 1 pha hạ áp kết nối lồng ghép với bài toán tối ưu kinh tế nguồn truyền thống trung/cao áp 3 pha có xét đến hệ số có lợi vị trí ($\alpha_i$). Mô hình giải quyết triệt để bài toán trào lưu công suất ngược (reverse power flow) từ phụ tải phân tán lên lưới truyền tải trong điều kiện lưới điện không có thiết bị OLTC.
2. Điểm cải tiến phương pháp luận cốt lõi khi so sánh với các nghiên cứu tiền nhiệm quốc tế?
Trả lời: So với nghiên cứu của Moradi & Abedini (2012) [20] và Yammani et al. (2016) [21] (chỉ tối ưu trên lưới trung áp chuẩn IEEE 3 pha đối xứng), luận án đã thiết lập quy trình giải tích chi tiết trên từng pha hạ áp 0,4 kV với tiết diện dây dẫn suy giảm theo khoảng cách cột. Đồng thời, so với thuật toán PSO nguyên bản của Kennedy & Eberhart (1995), luận án bổ sung toán tử cưỡng bức vận tốc $v_d = 0$ khi vi phạm biên $[x_{\min}, x_{\max}]$, loại bỏ hoàn toàn hiện tượng dao động kẹt biên trong bài toán tối ưu dung lượng nguyên.
3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng định lượng đi kèm?
Trả lời: Kết quả bất ngờ nhất là dung lượng tiếp nhận cực đại của lưới hạ áp có thể đạt 89 kW/pha, tương đương 267 kW/3 pha cho trạm biến áp 320 kVA (chiếm 83,4% công suất định mức), cao gấp gần 2,8 lần giới hạn 30% (96 kVA) theo Thông tư 39/2015/TT-BCT mà vẫn bảo toàn biên độ điện áp trong dải $209\text{ V} - 231\text{ V}$. Điều kiện tiên quyết là công suất phải được phân bổ lệch, tập trung chủ yếu tại Nút 1 (14 kW) và Nút 2 (65 kW) ở đầu xuất tuyến thay vì chia đều cho các nút.
4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?
Trả lời: Hoàn toàn đầy đủ. Luận án cung cấp chi tiết: bảng thông số trở kháng đơn vị ($r_0, x_0$), chiều dài khoảng cột (40 m), ma trận phụ tải từng nút, mã nguồn giải thuật GA/PSO trên nền MATLAB R2018a với các tham số: kích thước quần thể $\text{pop_size} = 50$, xác suất lai $p_c = 0,8$, đột biến $p_m = 0,02$, tham số gia tốc PSO ($c_1 = 0,8, c_{\max} = 1,62$), cấu trúc mạng nơ-ron truyền thẳng 2 lớp và thuật toán lan truyền ngược Levenberg-Marquardt.
5. Khung chương trình nghiên cứu 10 năm được định hình như thế nào?
Trả lời: Lộ trình 10 năm tập trung vào 3 trụ cột: (1) Nâng cấp mô hình giải tích sang hệ thống phân phối 3 pha 4 dây mất cân bằng động thời gian thực; (2) Ứng dụng trí tuệ nhân tạo tăng cường sâu (Multi-Agent Deep Reinforcement Learning) điều khiển phân tán hàng triệu xe điện (V2G); (3) Xây dựng thị trường giao dịch năng lượng phân tán ngang hàng (Peer-to-Peer) dựa trên công nghệ sổ cái phân tán (DLT).
Kết luận
- Xác lập phương pháp luận giải tích tiên phong: Luận án đã xây dựng thành công mô hình toán học và thuật toán tối ưu hóa xác định công suất cực đại của nguồn điện gió, điện mặt trời phát lên 1 pha lưới điện phân phối hạ áp 0,4 kV, kiểm soát nghiêm ngặt điện áp nút $U_i \in [209\text{ V}, 231\text{ V}]$ và dòng điện nhánh $I_{k-m} \le I_{\max}$.
- Đột phá về khả năng tiếp nhận công suất sạch: Bằng chứng thực nghiệm khẳng định lưới điện hạ áp trạm biến áp 320 kVA có thể hấp thụ an toàn 267 kW công suất nguồn phân tán 3 pha (đạt 83,4% công suất đặt của trạm), phá vỡ giới hạn trần 30% (96 kVA) của quy định hiện hành tại Việt Nam.
- Tối ưu hóa đa tầng lồng ghép kinh tế - kỹ thuật: Thiết lập thành công mô hình điều độ tối ưu trào lưu công suất OPF cho lưới truyền tải trung/cao áp kết hợp giải tỏa công suất ngược từ hạ áp, tự động điều chỉnh nguồn truyền thống nhằm tối thiểu hóa chi phí phát điện toàn hệ thống và tích hợp hệ số có lợi vị trí ($\alpha_i$).
- Khung dự báo và điều khiển phụ tải linh hoạt: Kết hợp mạng nơ-ron nhân tạo huấn luyện bằng thuật toán Levenberg-Marquardt dự báo chuỗi thời gian 30 phút/mẫu với thuật toán di truyền điều khiển đóng/cắt phụ tải dịch chuyển, nâng cao tính tự cân bằng của hệ thống điện phân phối.
- Mô hình tích hợp xe điện nâng cao độ tin cậy: Đề xuất cấu trúc trạm sạc xe điện kết hợp điện mặt trời áp mái vận hành ở chế độ tách đảo, chứng minh định lượng qua mô phỏng Monte Carlo khả năng cắt giảm chỉ số thời gian mất điện SAIDI từ 24,5 giờ xuống 8,2 giờ/khách hàng/năm.
- Đóng góp nền tảng cho tiến trình chuyển dịch năng lượng: Toàn bộ kết quả nghiên cứu là tài liệu tham khảo học thuật giá trị cao, đồng thời là cẩm nang kỹ thuật thực tiễn thúc đẩy tự động hóa lưới điện thông minh, phục vụ đắc lực cho chiến lược phát triển bền vững và an ninh năng lượng quốc gia.