Tổng quan về luận án

Nghiên cứu phụ tải và dự báo phụ tải ngắn hạn (Short-Term Load Forecasting - STLF) đóng vai trò sống còn trong việc đảm bảo an ninh năng lượng, tối ưu hóa chi phí vận hành hệ thống điện và quyết định trực tiếp đến hiệu quả kinh tế của các tác nhân tham gia thị trường điện cạnh tranh. Trong bối cảnh Việt Nam chuyển đổi mạnh mẽ từ mô hình độc quyền sang Thị trường bán buôn điện cạnh tranh (Vietnam Wholesale Electricity Market - VWEM) theo định hướng của Bộ Công Thương, bài toán STLF không còn dừng lại ở cấp độ kỹ thuật điều độ thuần túy mà đã trở thành công cụ chiến lược sống còn để tối ưu hóa giá chào, phân bổ danh mục hợp đồng sai khác (Contract for Differences - CFD) và giảm thiểu rủi ro biến động giá cận biên thị trường (System Marginal Price - SMP).

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn vận hành: các mô hình kinh tế lượng truyền thống và thuật toán chuỗi thời gian cổ điển thường bộc lộ sai số lớn khi đối mặt với dữ liệu chuỗi thời gian phụ tải mang tính phi tuyến cao, chịu sự tác động đồng thời của biến đổi khí hậu nhiệt đới gió mùa và đặc thù lệch pha chu kỳ văn hóa giữa lịch dương và lịch âm (điển hình là biến động phụ tải dị biệt trong dịp Tết Nguyên Đán tại các tháng 01 và 02 hàng năm). Hơn nữa, vấn đề dữ liệu đo đếm bị gián đoạn, mất gói tin, xuất hiện điểm ngoại lai (outliers) hoặc giá trị bằng 0 giả tạo (false zeros) do sự cố thiết bị đo đếm và đường truyền viễn thông chưa được xử lý bằng các công cụ toán học tự động hóa mà vẫn phụ thuộc nặng nề vào kinh nghiệm thủ công.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học mang tính hệ thống:

  • RQ1: Làm thế nào để tự động hóa hoàn toàn quy trình nhận dạng, làm sạch và khôi phục các chuỗi dữ liệu phụ tải bị thiếu, lỗi hoặc ngoại lai với độ chính xác cao?
  • RQ2: Cơ chế tích hợp giữa Biểu đồ phụ tải chuẩn hóa đơn vị (Standardized Load Profiles - SLP) và Hồi quy Vector Hỗ trợ (Support Vector Regression - SVR) có thể triệt tiêu được các sai lệch do yếu tố phi tuyến của thời tiết và sự chênh lệch chu kỳ âm - dương lịch hay không?
  • RQ3: Mô hình STLF đề xuất tác động như thế nào đến chiến lược phân bổ sản lượng hợp đồng Qc và tối ưu hóa lợi nhuận của Tổng công ty Điện lực trong thị trường VWEM?
  • H1: Mô hình học máy dựa trên nguyên lý cực tiểu hóa rủi ro cấu trúc (SVR) kết hợp hàm nhân RBF vượt trội hơn các phương pháp thống kê cổ điển trong việc khôi phục dữ liệu phi tuyến ngẫu nhiên.
  • H2: Việc sử dụng biểu đồ phụ tải chuẩn hóa đơn vị (SLP) 24 chu kỳ/ngày làm đầu vào cho mô hình hồi quy SVR sẽ làm giảm đáng kể sai số dự báo trong các giai đoạn chuyển tiếp chu kỳ phụ tải phức tạp (Tết Âm lịch và các kỳ nghỉ lễ).
  • H3: Sai số dự báo STLF tiệm cận mức tối ưu sẽ giảm thiểu trực tiếp chi phí bù trừ sai lệch công suất và nâng cao biên lợi nhuận của đơn vị phân phối điện trên thị trường giao ngay.

Khung lý thuyết của công trình được xây dựng trên sự giao thoa giữa Lý thuyết Học Thống kê (Statistical Learning Theory của Vapnik), Phương pháp luận Chuỗi thời gian Box-Jenkins và Lý thuyết Cân bằng Thị trường Điện lực. Đóng góp đột phá của luận án là thiết kế thành công kiến trúc lai ghép thông minh SLP-SVR, giải quyết triệt để bài toán khôi phục dữ liệu lỗi và nâng cao độ chính xác dự báo 24/48 chu kỳ phụ tải ngày tới và tuần tới. Phạm vi nghiên cứu bao trùm chuỗi số liệu vận hành thực tế tại Tổng công ty Điện lực TP.HCM (EVNHCMC) từ giai đoạn 1997–2018 (trong đó tập trung xử lý chi tiết bộ dữ liệu đo đếm chu kỳ 30–60 phút giai đoạn 2014–2016) và dữ liệu vận hành VWEM từ Trung tâm Điều độ Hệ thống điện Quốc gia (A0), đối chiếu trên quy mô 05 Tổng công ty Điện lực trực thuộc Tập đoàn Điện lực Việt Nam (EVN).

Literature Review và Positioning

Nghiên cứu dự báo phụ tải ngắn hạn trên thế giới đã trải qua ba giai đoạn tiến hóa phương pháp luận rõ rệt: trường phái kinh tế lượng cổ điển, trường phái chuỗi thời gian thống kê học và trường phái trí tuệ nhân tạo/học máy tiên tiến.

Giai đoạn đầu tập trung vào các phương pháp hồi quy tuyến tính đa biến và san bằng số mũ (Winter, 1960). Tuy nhiên, các phân tích của Barakat et al. (1990) đã chỉ ra rằng mô hình nhu cầu năng lượng tại các khu vực đang phát triển với tốc độ đô thị hóa nhanh rất khó nắm bắt nếu chỉ áp dụng đơn thuần phương pháp làm mịn hàm mũ của Winter. Để khắc phục điều này, El-Keib et al. (1995) đã đề xuất phương pháp lai ghép kết hợp làm mịn hàm mũ với phân tích phổ điện và mô hình tự điều chỉnh tự hồi quy.

Trường phái chuỗi thời gian Box-Jenkins với các biến thể AR, MA, ARMA, ARIMA và SARIMA đánh dấu bước chuyển mang tính triết lý sâu sắc: "hãy để dữ liệu tự nói" thay vì áp đặt các giả định hành vi kinh tế định trước. Tuy vậy, các tranh luận học thuật lớn đã nổ ra xoay quanh việc xử lý tính mùa vụ và cấu trúc chuỗi thời gian trong STLF. Weron (2006) và Taylor (2008) cho rằng trong STLF, chu kỳ dự báo ngắn hơn rất nhiều so với độ dài chu kỳ năm, do đó không nhất thiết phải mô hình hóa tính mùa vụ năm mà cần tập trung nắm bắt tính chu kỳ kép (intra-day và intra-week) thông qua mô hình SARIMA mùa kép hoặc làm mịn hàm mũ mở rộng. Ngược lại, Soares và Medeiros (2008) cùng Dordonnat et al. (2008) lập luận rằng cần phải tích hợp tường minh xu hướng dài hạn (như hàm phụ thuộc GDP) và mô hình hóa chu kỳ năm bằng các tổ hợp hàm điều hòa sin-cosin tương tự phân tích chuỗi Fourier.

Một cuộc tranh luận phương pháp luận khác diễn ra giữa cách tiếp cận đơn phương trình toàn cục và cách tiếp cận đa phương trình cục bộ cho 24 giờ trong ngày. Ramanathan et al. (1997), Cottet và Smith (2003) chứng minh rằng việc chia tách đồ thị phụ tải thành 24 chuỗi thời gian độc lập cho 24 giờ riêng biệt mang lại độ chính xác cao hơn hẳn mô hình đơn phương trình do triệt tiêu được tương quan chéo phức tạp giữa các giờ cao điểm và thấp điểm.

Về biến ngoại sinh khí tượng, khảo sát của Hippert et al. (2001) và Weron (2006) cho thấy nhiệt độ là biến số tác động phi tuyến mạnh mẽ nhất. Engle et al. (1986), Cancelo et al. (2008) và Moral-Carcedo và Vicéns-Otero (2005) chỉ ra rằng tương quan giữa nhiệt độ và phụ tải mang đặc tính phi tuyến chữ U hoặc chữ V: phụ tải tăng vọt ở cả hai cực nhiệt độ thấp (sưởi ấm) và nhiệt độ cao (điều hòa không khí). Để giải quyết phi tuyến này, Sánchez-Úbeda và Wehenkel (1998) đề xuất mô hình phân đoạn nhiệt độ (vùng lạnh $<15^\circ\text{C}$, vùng trung gian $15-17^\circ\text{C}$, và vùng nóng $>17^\circ\text{C}$), trong khi Taylor và Buizza (2003) phát triển mô hình học máy tổ hợp dự báo phụ tải dựa trên đa kịch bản khí tượng.

Tại Việt Nam, các nghiên cứu trước đây như đề tài cấp Bộ của Trần Kỳ Phúc (2008) sử dụng mạng Nơ-ron đa lớp lan truyền ngược (MLP) và mạng tự tổ chức Kohonen; Nguyễn Hoàng Việt et al. (2007) áp dụng mạng Wavelet cho các ngày đặc biệt; Trần Thị Hoàng Oanh et al. (2007) sử dụng mạng Neural song tuyến; hay Lê Đình Dương ứng dụng thuật toán xử lý dữ liệu nhóm (GMDH) xây dựng mạng nơ-ron đa thức. Mặc dù đạt được những kết quả học thuật nhất định, các nghiên cứu này còn nặng tính phương pháp luận mô phỏng, chưa giải quyết được bài toán làm sạch dữ liệu tự động, chưa xử lý được độ trễ âm - dương lịch đặc thù và chưa gắn kết với cơ chế thanh toán thực tế của thị trường điện bán buôn. Luận án đã định vị chính xác khoảng trống này để phát triển giải pháp tích hợp toàn diệ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 Lý thuyết Học Thống kê (Statistical Learning Theory) của Vladimir Vapnik thông qua việc ứng dụng hồi quy vector hỗ trợ (SVR) trên không gian đặc trưng nhiều chiều vào chuỗi thời gian năng lượng có nhiễu lớn. Khác với nguyên lý Cực tiểu hóa Rủi ro Thực nghiệm (Empirical Risk Minimization - ERM) của mạng nơ-ron truyền thống dễ dẫn đến hiện tượng quá khớp (overfitting) và bẫy cực trị địa phương, SVR triển khai nguyên lý Cực tiểu hóa Rủi ro Cấu trúc (Structural Risk Minimization - SRM) nhằm tìm kiếm một siêu phẳng tối ưu có khả năng tổng quát hóa cao nhất.

       +-------------------------------------------------------------+
       |   KHUNG PHÂN TÍCH TÍCH HỢP ĐỘC ĐÁO: KIẾN TRÚC HYBRID SLP-SVR |
       +-------------------------------------------------------------+
                                      |
     +--------------------------------+--------------------------------+
     |                                                                 |
     v                                                                 v
+-----------------------------+                   +----------------------------------+
| PHÂN HỆ XỬ LÝ VÀ ĐIỀN ĐẦY   |                   | PHÂN HỆ TIỀN XỬ LÝ KHÔNG GIAN TẢI|
| DỮ LIỆU THIẾU (DATA ENGINE) |                   | (SLP PROFILING ENGINE)           |
+-----------------------------+                   +----------------------------------+
| - Thu thập chuỗi Pmax, Atong|                   | - Phân rã chuỗi thời gian:       |
| - Phân loại cơ chế mất dữ   |                   |   Trend, Cycle, Season, Irregular|
|   liệu (MCAR / MAR / NMAR)  |                   | - Chuẩn hóa đồ thị 24 chu kỳ/ngày|
| - Nhận dạng dị thường       |                   | - Khử nhiễu biến động dương lịch |
|   (Distance / Stat / Dev)   |                   |   và âm lịch (Tết Nguyên Đán)    |
| - Tái tạo chuỗi qua mô hình |                   | - Xây dựng cơ sở dữ liệu vector  |
|   SVR với hàm mất mát eps   |                   |   tải chuẩn hóa SLP              |
+-----------------------------+                   +----------------------------------+
     |                                                                 |
     +--------------------------------+--------------------------------+
                                      |
                                      v
          +-------------------------------------------------------+
          | PHÂN HỆ DỰ BÁO PHỤ TẢI NGẮN HẠN HỌC MÁY TỐI TẬP (SVR) |
          +-------------------------------------------------------+
          | - Ánh xạ không gian đặc trưng phi tuyến qua Kernel    |
          |   Gaussian RBF thỏa mãn Định lý Mercer                |
          | - Tối ưu hóa điều kiện KKT và biến lỏng (Slack vars)  |
          | - Tối ưu bộ tham số siêu cấp (C, epsilon, gamma)      |
          | - Dự báo chính xác phụ tải 24/48 chu kỳ ngày tới      |
          +-------------------------------------------------------+
                                      |
                                      v
          +-------------------------------------------------------+
          | PHÂN HỆ ỨNG DỤNG VẬN HÀNH THỊ TRƯỜNG ĐIỆN VWEM        |
          +-------------------------------------------------------+
          | - Tối ưu bảng chào sản lượng Qc theo Thông tư 45/BCT   |
          | - Cân đối danh mục sai khác CFD và giá thị trường SMP |
          | - Cực đại hóa doanh thu, kiểm soát rủi ro tài chính   |
          +-------------------------------------------------------+

Mô hình lý thuyết được xây dựng thông qua ba mệnh đề (propositions):

  • Mệnh đề 1 (Đặc tính thưa của nghiệm): Nghiệm của bài toán hồi quy phụ tải phi tuyến có thể biểu diễn tối ưu thông qua một tập con thưa gồm các Vector Hỗ trợ (Support Vectors) nằm trên hoặc ngoài biên của ống sai số $\varepsilon$-insensitive, cho phép nén thông tin và tăng tốc độ xử lý dữ liệu lớn.
  • Mệnh đề 2 (Tính bất biến của biểu đồ chuẩn hóa): Phụ tải điện khi được chuẩn hóa thành dạng đồ thị đơn vị (SLP) sẽ triệt tiêu được độ lệch quy mô tuyệt đối, bảo toàn hình thái tiêu thụ đặc trưng theo loại ngày, từ đó cho phép ánh xạ chính xác chu kỳ phụ tải âm lịch sang hệ quy chiếu dương lịch.
  • Mệnh đề 3 (Cơ chế truyền dẫn rủi ro giá điện): Độ lệch tuyệt đối giữa phụ tải dự báo và phụ tải thực tế có mối liên hệ phi tuyến bậc hai với chi phí bù trừ sai lệch công suất trên thị trường giao ngay SMP khi đơn vị phân phối không tối ưu hóa lượng hợp đồng Qc.

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: (1) Phương pháp phân tích chuỗi thời gian kinh tế lượng để trích xuất các thành phần xu hướng (trend), chu kỳ (cyclical), mùa vụ (seasonal) và bất quy tắc (irregular); (2) Thuật toán học máy SVR với hàm nhân Gauss; và (3) Kỹ thuật chuẩn hóa đồ thị phụ tải đơn vị (SLP).

Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ: áp dụng cho hệ thống điện phân phối cấp vùng và đô thị với chu kỳ đo đếm 30–60 phút, độ trễ dự báo từ 24 giờ đến 07 ngày (ngắn hạn), trong điều kiện khí hậu nhiệt đới gió mùa có phân hóa hai mùa mưa - nắng rõ rệt và thị trường điện vận hành theo cơ chế chào giá tập trung có thanh toán bù trừ qua hợp đồng sai khác CFD.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị theo hệ hình thực chứng (positivism) với lập trường nhận thức luận khách quan, tiếp cận hoàn toàn dựa trên dữ liệu định lượng thực tế (data-driven empirical research). Thiết kế nghiên cứu theo mô hình đa cấp độ (multi-level design):

  • Cấp độ vĩ mô: Hệ thống điện toàn quốc và thị trường điện bán buôn VWEM.
  • Cấp độ trung mô: 05 Tổng công ty Điện lực phân phối (EVNHCMC, EVNHANOI, EVNCPC, EVNNPC, EVNSPC).
  • Cấp độ vi mô: Các đồ thị phụ tải 24/48 chu kỳ ngày và từng giờ riêng biệt.

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu được chuẩn hóa theo chuẩn quốc tế nhằm giải quyết triệt để ba cơ chế khuyết thiếu dữ liệu: Thiếu hoàn toàn ngẫu nhiên (Missing Completely at Random - MCAR), Thiếu ngẫu nhiên (Missing at Random - MAR) và Thiếu không ngẫu nhiên (Not Missing at Random - NMAR).

Quy trình phát hiện phần tử ngoại lai kết hợp đồng thời ba kỹ thuật:

  1. Xác định theo khoảng cách (Distance-Based metrics).
  2. Xác định theo phân phối thống kê chuẩn và phân phối Chi-bình phương (Statistical-Based).
  3. Xác định theo độ sai lệch đặc trưng cơ bản của tập phần tử (Deviation-Based).

Luận án áp dụng phương pháp tam giác đạc dữ liệu (data triangulation) bằng cách đối chiếu chéo số liệu công suất cực đại ($P_{\max}$), công suất cực tiểu ($P_{\min}$), điện năng tổng tiêu thụ ($A_{\text{tổng}}$), biến thiên nhiệt độ thực tế ($t^\circ\text{C}$), độ ẩm, bức xạ mặt trời và các mốc chỉ số phụ tải hệ thống từ A0.

Data và phân tích

Bộ dữ liệu thực nghiệm bao gồm chuỗi dữ liệu phụ tải theo giờ và theo tháng của EVNHCMC từ năm 1997 đến 2018. Để kiểm chứng thuật toán khôi phục dữ liệu lỗi, tác giả thiết lập kịch bản thử nghiệm trên chuỗi đo đếm hoàn chỉnh từ ngày 01/01/2014 đến 28/02/2016 và chủ động mô phỏng các tình huống lỗi điển hình ghi nhận trong thực tế:

"04 số liệu điện năng (ngày 10 – 13/01/2015) bị lỗi ghi nhận bằng 0 – không đo đếm được; 03 số liệu điện năng (ngày 12 – 14/02/2015) bị lỗi ghi nhận bằng rỗng – không lưu trữ được và 03 số liệu điện năng (ngày 05 – 07/03/2015) bị lỗi ghi cao bất thường."

          +-----------------------------------------------------------------------+
          |     THIẾT KẾ THỰC NGHIỆM ĐIỀN ĐẦY VÀ KHÔI PHỤC DỮ LIỆU ĐO ĐẾM PHỤ TẢI |
          +-----------------------------------------------------------------------+
                                              |
       +--------------------------------------+--------------------------------------+
       |                                      |                                      |
       v                                      v                                      v
+------------------------------+ +------------------------------+ +------------------------------+
| KỊCH BẢN 1: GIÁ TRỊ BẰNG 0   | | KỊCH BẢN 2: GIÁ TRỊ RỖNG     | | KỊCH BẢN 3: NGOẠI LAI CAO    |
| (FALSE ZEROES)               | | (MISSING / NULL DATA)        | | (TOO HIGH / ANOMALIES)       |
+------------------------------+ +------------------------------+ +------------------------------+
| 04 ngày liên tiếp            | | 03 ngày liên tiếp            | | 03 ngày liên tiếp            |
| (10/01/2015 - 13/01/2015)    | | (12/02/2015 - 14/02/2015)    | | (05/03/2015 - 07/03/2015)    |
| Nguyên nhân: Sự cố truyền tín| | Nguyên nhân: Sự cố bộ lưu trữ| | Nguyên nhân: Nhiễu cảm biến, |
| hiệu công tơ về Kho dữ liệu  | | viễn thông, đứt đoạn gói tin | | lỗi đo đếm xung công tơ      |
+------------------------------+ +------------------------------+ +------------------------------+
       |                                      |                                      |
       +--------------------------------------+--------------------------------------+
                                              |
                                              v
          +-----------------------------------------------------------------------+
          | GIẢI THUẬT HỒI QUY VECTOR HỖ TRỢ TỐI ƯU (SVR DATA RECONSTRUCTION)     |
          | - Hàm lỗi: eps-insensitive loss function kết hợp biến lỏng xi_n       |
          | - Không gian nhân: Gaussian RBF Kernel k(x, x') = exp(-gamma||x-x'||^2)|
          | - Kết quả: Tái lập chuỗi hoàn chỉnh, triệt tiêu sai số cục bộ         |
          +-----------------------------------------------------------------------+

Về mặt toán học, thuật toán SVR tối thiểu hóa hàm mục tiêu có xét biến lỏng: $$\min_{w, b, \xi, \hat{\xi}} \frac{1}{2} |w|^2 + C \sum_{n=1}^{N} (\xi_n + \hat{\xi}_n)$$ với các ràng buộc điều kiện Karush-Kuhn-Tucker (KKT): $$\begin{cases} t_n - y(x_n) \le \varepsilon + \xi_n \ y(x_n) - t_n \le \varepsilon + \hat{\xi}_n \ \xi_n, \hat{\xi}_n \ge 0 \end{cases}$$

Chuyển đổi sang bài toán đối ngẫu Lagrange: $$\max_{a, \hat{a}} \sum_{n=1}^{N} t_n(a_n - \hat{a}n) - \varepsilon \sum{n=1}^{N} (a_n + \hat{a}n) - \frac{1}{2} \sum{n=1}^{N} \sum_{m=1}^{N} (a_n - \hat{a}_n)(a_m - \hat{a}_m) k(x_n, x_m)$$ với hàm nhân Gauss thỏa mãn điều kiện Mercer: $$k(x_i, x_j) = \exp\left(-\gamma |x_i - x_j|^2\right)$$ với ràng buộc: $0 \le a_n, \hat{a}n \le C$ và $\sum{n=1}^{N} (a_n - \hat{a}_n) = 0$.

Toàn bộ thuật toán xử lý dữ liệu, trích xuất chuỗi thời gian ACF/PACF, thuật toán Rừng ngẫu nhiên (Random Forest Regression dựa trên Bagging và CART của Breiman, 2001), mô hình SARIMA và mô hình SVR được lập trình và mô phỏng trên nền tảng MATLAB chuyên sâu.

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

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

  1. Khả năng tự động hóa khôi phục dữ liệu lỗi đạt độ chính xác vượt trội: Mô hình hồi quy SVR với tham số tối ưu hóa ($C, \varepsilon, \gamma$) kết hợp đồ thị SLP đã khôi phục thành công các chuỗi dữ liệu bị mất (04 ngày mất dữ liệu tháng 01/2015, 03 ngày rỗng tháng 02/2015 và 03 ngày ngoại lai tháng 03/2015). Sai số tuyệt đối trung bình (MAD) và sai số phần trăm tuyệt đối trung bình (MAPE) của dữ liệu khôi phục so với dữ liệu gốc đạt dưới $1.5%$, vượt trội hoàn toàn so với các phương pháp truyền thống như nội suy tuyến tính hay gán nhãn trung bình 4 tuần gần nhất ($p < 0.001$).
  2. Hóa giải triệt để bài toán biến động phụ tải dị biệt Tết Âm lịch: Bằng việc sử dụng biểu đồ phụ tải chuẩn hóa đơn vị (SLP) làm dữ liệu đầu vào, mô hình loại bỏ được hiện tượng méo dạng biểu đồ trong tháng 01 và tháng 02 qua các năm mà không cần phải thực hiện các bước quy đổi thủ công phức tạp như các phần mềm thương mại ngoại nhập.
  3. Tính ưu việt của giải pháp thưa (Sparse Solution) trong SVR so với Random Forest và SARIMA: Thử nghiệm đối sánh cho thấy mô hình SVR chỉ sử dụng một tập hợp con nhỏ các Support Vectors thay vì toàn bộ không gian mẫu, giúp giảm thời gian tính toán $45%$ so với mạng Neural lan truyền thẳng nhiều lớp và giảm hiện tượng phân mảnh dữ liệu so với cây quyết định CART trong Random Forest.
  4. Xác lập ngưỡng tương quan nhiệt độ phi tuyến tại đô thị nhiệt đới: Kết quả phân tích chuỗi dữ liệu EVNHCMC chỉ ra mối tương quan phi tuyến chặt chẽ giữa nhiệt độ môi trường và đỉnh phụ tải $P_{\max}$. Khác với các nước hàn đới (có đồ thị phụ tải tăng ở cả 2 đầu nhiệt độ), tại TP.HCM phụ tải biến thiên đồng pha mạnh mẽ khi nhiệt độ vượt ngưỡng $28^\circ\text{C}$ do hiệu ứng sử dụng thiết bị làm mát đồng loạt, với hệ số tương quan đạt đỉnh tại các giờ trưa và đầu giờ chiều.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp khung phương pháp luận mới trong việc xử lý chuỗi thời gian năng lượng có tính mùa vụ đa tầng và chu kỳ lịch hỗn hợp (Solar-Lunar Calendar), mở rộng ứng dụng của lý thuyết học máy thống kê trong kinh tế năng lượng.
  • Về mặt phương pháp: Tạo ra một pipeline chuẩn tắc từ khâu tiền xử lý, tự động làm sạch dữ liệu ngoại lai, chuẩn hóa hình thái phụ tải đến mô hình hóa dự báo có thể chuyển giao cho nhiều bài toán chuỗi thời gian hạ tầng khác (giao thông, cấp thoát nước, viễn thông).
  • Về mặt thực tiễn và chính sách: Giúp các Tổng công ty Điện lực đáp ứng nghiêm ngặt quy định của Cục Điều tiết Điện lực (ERAV) theo Thông tư 45/2018/TT-BCT:

"Thông tư 45/2018/TT-BCT ngày 15 tháng 11 năm 2018 của Bộ Công Thương về Quy định vận hành Thị trường bán buôn điện cạnh tranh" quy định "các Tổng công ty Điện lực phải công bố bảng chào Qc tháng N+1 24 chu kỳ trước ngày 20 của tháng N".

Mô hình dự báo có độ chính xác cao cho phép đơn vị phân phối lập kế hoạch chào mua sản lượng hợp đồng sai khác Qc tối ưu, hạn chế tối đa rủi ro phải mua điện bù trừ trên thị trường giao ngay với giá SMP tăng vọt trong các giờ cao điểm, từ đó bảo toàn dòng tiền và tối ưu hóa lợi nhuận kinh doanh.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật và thực tiễn vững chắc, công trình vẫn tồn tại một số giới hạn nghiên cứu:

  1. Độ phụ thuộc vào chất lượng dự báo khí tượng: Mô hình STLF phụ thuộc vào nguồn dữ liệu dự báo nhiệt độ đầu vào. Trong bối cảnh biến đổi khí hậu cực đoan, độ tin cậy của các bản tin dự báo thời tiết thường suy giảm mạnh sau khoảng thời gian 05 đến 07 ngày, gây khó khăn cho công tác lập kế hoạch trung hạn.
  2. Hiện tượng thời tiết dị thường quy mô lớn: Các pha khí hậu cực đoan như El Niño kéo dài hoặc La Niña gây mưa bất thường chưa được tích hợp như một biến trạng thái ngẫu nhiên độc lập trong mô hình hồi quy.
  3. Chi phí tính toán của ma trận Kernel: Khi kích thước tập dữ liệu huấn luyện mở rộng lên quy mô hàng chục năm với chu kỳ lấy mẫu 15–30 phút, việc tính toán ma trận Gram trong không gian đối ngẫu SVR đòi hỏi tài nguyên bộ nhớ lớn.
  4. Phạm vi kiểm chứng: Mặc dù đã đối chiếu trên 05 Tổng công ty Điện lực, thực nghiệm chuyên sâu tập trung chủ yếu vào lưới điện phân phối đô thị (EVNHCMC) với tỷ trọng phụ tải sinh hoạt và thương mại dịch vụ cao, chưa bao quát đầy đủ đặc thù phụ tải công nghiệp nặng chuyên biệt.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra các hướng phát triển:

  • Ứng dụng các kiến trúc Học Sâu tiên tiến (Deep Learning) như mạng Bộ nhớ Ngắn-Dài (LSTM), Gated Recurrent Units (GRU) và Temporal Fusion Transformers (TFT) để khai thác các phụ thuộc thời gian siêu dài.
  • Tích hợp động thái phụ tải phân tán (Distributed Energy Resources - DER) và nguồn điện mặt trời mái nhà tự sản tự tiêu vào biểu đồ tải chuẩn hóa SLP.
  • Phát triển thuật toán tối ưu hóa siêu tham số thích nghi trực tuyến (Online Adaptive Hyperparameter Tuning) dựa trên thuật toán Di truyền (Genetic Algorithm) hoặc Tối ưu hóa Bầy đàn (PSO) để cập nhật trọng số mô hình theo thời gian thực.
  • Mở rộng mô hình dự báo xác suất (Probabilistic Load Forecasting) cung cấp các khoảng tin cậy thay vì dự báo điểm đơn lẻ.

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

Công trình tạo ra những tác động lan tỏa mạnh mẽ trên nhiều phương diện:

+-----------------------------------------------------------------------------------+
|                        MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG ĐA TẦNG                      |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [HỌC THUẬT & NGHIÊN CỨU]                    [DOANH NGHIỆP & NGÀNH ĐIỆN]          |
|  - Cung cấp mô hình tham chiếu cho           - Ứng dụng trực tiếp tại EVNHCMC và  |
|    các nghiên cứu STLF trong nước              05 Tổng công ty Điện lực.          |
|  - Tiềm năng trích dẫn cao trong các         - Giảm chi phí sai lệch công suất,   |
|    tạp chí ISI/Scopus chuyên ngành             tối ưu hóa danh mục hợp đồng Qc.   |
|    năng lượng và hệ thống thông minh.        - Tự động hóa giám sát Kho dữ liệu.  |
|                                                                                   |
|                                                                                   |
|  [CHÍNH SÁCH & CƠ QUAN QUẢN LÝ]             [XÃ HỘI & MÔI TRƯỜNG]                |
|  - Hỗ trợ ERAV hoàn thiện khung kỹ thuật     - Giảm thiểu rủi ro thiếu điện sự cố |
|    giám sát thị trường VWEM.                   và sa thải phụ tải khẩn cấp.       |
|  - Thúc đẩy lộ trình vận hành thị trường     - Tối ưu hóa vận hành nhà máy phát,  |
|    bán lẻ điện cạnh tranh (VREM).              cắt giảm phát thải CO2 biên.       |
|                                                                                   |
+-----------------------------------------------------------------------------------+
  • Tác động học thuật: Đặt nền móng lý thuyết vững chắc cho bài toán xử lý dữ liệu đo đếm bất đối xứng và dự báo chuỗi thời gian năng lượng tại các quốc gia đang phát triển; ước tính tạo ra tiềm năng trích dẫn khoa học cao trong cộng đồng nghiên cứu kỹ thuật hệ thống điện và khoa học dữ liệu ứng dụng.
  • Tác động ngành và doanh nghiệp: Triển khai trực tiếp thành công cụ phần mềm tại EVNHCMC, giúp tự động hóa quá trình giám sát dữ liệu đo đếm, giảm thiểu chi phí nhân công xử lý dữ liệu thủ công, nâng cao độ chính xác dự báo phụ tải, mang lại lợi ích tài chính hàng tỷ đồng mỗi năm thông qua tối ưu hóa giá mua điện bình quân.
  • Tác động chính sách: Cung cấp cơ sở khoa học định lượng cho Cục Điều tiết Điện lực (Bộ Công Thương) trong việc đánh giá năng lực dự báo phụ tải của các đơn vị thành viên thị trường, làm tiền đề hoàn thiện các thông tư hướng dẫn vận hành thị trường bán buôn và chuẩn bị cho thị trường bán lẻ điện cạnh tranh (VREM).
  • Lợi ích kinh tế - xã hội: Việc dự báo phụ tải chuẩn xác giúp hệ thống điện vận hành an toàn, giảm thiểu sự cố quá tải đường dây, giảm lượng công suất dự phòng quay không cần thiết của các nhà máy nhiệt điện, qua đó giảm tiêu hao nhiên liệu hóa thạch và cắt giảm phát thải khí nhà kính.

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

  • Nghiên cứu sinh và học giả: Tiếp cận một khung lý thuyết và phương pháp luận chặt chẽ về xử lý dữ liệu khuyết thiếu và ứng dụng học máy phi tuyến trong bài toán năng lượng; khai thác các hướng nghiên cứu mở về mô hình hóa chu kỳ lịch hỗn hợp.
  • Giảng viên và chuyên gia kinh tế năng lượng: Sử dụng các kết quả nghiên cứu và dẫn chứng thực nghiệm từ thị trường điện Việt Nam làm học liệu giảng dạy cao học và tài liệu tham khảo chuyên sâu.
  • Kỹ sư R&D và chuyên viên điều độ hệ thống điện: Khai thác các thuật toán và mã nguồn mô phỏng SVR-SLP để tích hợp vào các hệ thống quản lý năng lượng (EMS), hệ thống quản lý phân phối (DMS) và các phần mềm dự báo thương mại tại các Tổng công ty Điện lực.
  • Nhà hoạch định chính sách và cơ quan điều tiết (ERAV/EVN): Sở hữu công cụ đánh giá độc lập, minh bạch để kiểm tra bảng chào sản lượng của các đơn vị phân phối, hỗ trợ công tác quy hoạch phát triển nguồn và lưới điện đồng bộ.

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?

Đóng góp độc đáo nhất là việc thiết lập nguyên lý kết hợp giữa Biểu đồ phụ tải chuẩn hóa đơn vị (SLP) với Hồi quy Vector Hỗ trợ (SVR) trên cơ sở mở rộng Lý thuyết Học Thống kê của Vapnik. Luận án đã giải quyết bài toán biểu diễn không gian đặc trưng chuỗi thời gian bị nhiễu và đứt đoạn bằng cách ánh xạ chuỗi phụ tải thực tế sang không gian Hilbert nhiều chiều thông qua hàm nhân Gaussian RBF, đồng thời triệt tiêu độ lệch quy mô bằng kỹ thuật chuẩn hóa 24 chu kỳ tải đơn vị.

2. Tính đổi mới phương pháp luận thể hiện như thế nào khi so sánh với các nghiên cứu tiền nhiệm?

So với các nghiên cứu cổ điển sử dụng mô hình Box-Jenkins (SARIMA) hay mạng nơ-ron đa lớp (MLP) của Trần Kỳ Phúc (2008) và Taylor (2008), phương pháp của luận án mang tính đột phá ở 2 điểm:

  • Tự động hóa hoàn toàn khâu phát hiện ngoại lai và điền đầy dữ liệu thiếu thông qua tối ưu hóa lồi có ràng buộc (KKT conditions), chấm dứt việc xử lý thủ công theo cảm tính chuyên gia.
  • Không cần bẻ gãy hay phân tách dữ liệu phức tạp để xử lý giai đoạn Tết Âm lịch như các phần mềm thương mại ngoại nhập mà vẫn duy trì sai số dự báo ở mức tối ưu nhờ đặc tính thích nghi của không gian vector SLP.

3. Phát hiện thực nghiệm nào bất ngờ và mang lại giá trị thực tế cao nhất?

Phát hiện bất ngờ nhất là tính ổn định vượt trội của hình thái đồ thị phụ tải chuẩn hóa đơn vị (SLP) ngay cả trong những giai đoạn biến động mạnh về sản lượng điện năng do thay đổi thời tiết hoặc nghỉ lễ kéo dài. Dữ liệu thực nghiệm chứng minh rằng biến động của tổng nhu cầu ($A_{\text{tổng}}$) và công suất đỉnh ($P_{\max}$) có thể rất lớn, nhưng tỷ trọng phân bổ công suất tương đối giữa 24 chu kỳ trong ngày tuân theo các quy luật phân phối xác định, cho phép mô hình SVR nắm bắt được đặc trưng phụ tải cốt lõi mà không bị quá khớp bởi các xung nhiễu ngẫu nhiên.

4. Luận án có cung cấp quy trình và giao thức tái lập nghiên cứu (Replication Protocol) hay không?

Luận án cung cấp đầy đủ các phương trình toán học thiết lập bài toán tối ưu hóa đối ngẫu, quy trình tiền xử lý dữ liệu khuyết thiếu (MCAR/MAR/NMAR), các tiêu chí lựa chọn hàm nhân Gaussian, phương pháp xác định tham số ($C, \varepsilon, \gamma$) và cấu trúc các module mô phỏng trên nền tảng MATLAB. Bộ dữ liệu thử nghiệm từ EVNHCMC và các kịch bản lỗi giả lập được mô tả chi tiết, cho phép các nhà nghiên cứu độc lập hoàn toàn có thể lập trình và tái tạo lại toàn bộ kết quả thực nghiệm.

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 tới tập trung vào 3 trọng tâm:

  1. Mở rộng kiến trúc lai sang các mạng nơ-ron học sâu thế hệ mới (Deep Spatio-Temporal Graph Neural Networks) để mô hình hóa đồng thời tương quan không gian giữa các nút lưới điện và chuỗi thời gian.
  2. Tích hợp dữ liệu lớn từ hệ thống đo đếm thông minh (Advanced Metering Infrastructure - AMI) với chu kỳ thu thập 5–15 phút nhằm phục vụ bài toán điều khiển phụ tải linh hoạt (Demand Response) theo thời gian thực.
  3. Thích ứng mô hình với thị trường bán lẻ điện cạnh tranh hoàn chỉnh (VREM), nơi người tiêu dùng đầu cuối tham gia trực tiếp vào thị trường điện thông qua các biểu giá biến đổi theo thời gian thực (Real-time Pricing).

Kết luận

  1. Xây dựng thành công công cụ tự động hóa xử lý dữ liệu đo đếm lỗi: Luận án đã giải quyết căn bản bài toán điền đầy dữ liệu thiếu, loại bỏ điểm ngoại lai và khôi phục chuỗi thời gian đo đếm gián đoạn với độ chính xác cao, thay thế hoàn toàn các thao tác thủ công phụ thuộc kinh nghiệm.
  2. Đề xuất và hiện thực hóa giải thuật lai ghép mới SLP-SVR: Ứng dụng thành công biểu đồ phụ tải chuẩn hóa đơn vị làm đầu vào cho mô hình hồi quy vector hỗ trợ, tạo ra bước nhảy vọt về độ chính xác dự báo STLF trong các điều kiện vận hành phức tạp và các dịp lễ Tết cổ truyền.
  3. Thực hiện khảo sát đối sánh phương pháp luận đa diện: Đánh giá toàn diện và định lượng ưu - nhược điểm của các dòng thuật toán từ thống kê cổ điển (ARIMA, SARIMA, San bằng số mũ Winter), trí tuệ nhân tạo (Mạng nơ-ron truyền thẳng, GMDH) đến máy học tiên tiến (Random Forest, SVR), xác lập cơ sở khoa học vững chắc cho việc lựa chọn hàm hồi quy tối ưu.
  4. Giải quyết xuất sắc bài toán tối ưu hóa kinh tế trong thị trường điện VWEM: Kết nối trực tiếp kết quả dự báo phụ tải ngắn hạn với chiến lược chào sản lượng hợp đồng Qc theo quy định của Thông tư 45/2018/TT-BCT, giúp đơn vị phân phối điện kiểm soát hiệu quả rủi ro biến động giá thị trường SMP và cực đại hóa lợi nhuận kinh tế.
  5. Mở ra các hướng nghiên cứu liên ngành giàu tiềm năng: Đặt nền tảng cho việc tích hợp mô hình học máy vào hệ thống lưới điện thông minh (Smart Grid), quản lý năng lượng tái tạo phân tán và vận hành thị trường bán lẻ điện cạnh tranh tại Việt Nam trong các giai đoạn phát triển tiếp theo.