Tổng quan nghiên cứu

Trong bối cảnh tái cơ cấu ngành năng lượng toàn cầu, chi phí nhiên liệu sơ cấp thường chiếm tới khoảng 60% đến 70% tổng chi phí sản xuất của các nhà máy điện. Việc tối ưu hóa vận hành thông qua công tác điều độ có thể giúp tiết kiệm hàng chục triệu USD mỗi năm cho toàn hệ thống. Bài toán Điều độ Kinh tế (Economic Dispatch - ED) đóng vai trò then chốt trong việc xác định chế độ vận hành tin cậy, tối ưu phân bổ công suất phát giữa các tổ máy nhằm đáp ứng biểu đồ phụ tải. Khi chuyển dịch sang mô hình thị trường điện cạnh tranh, mục tiêu của bài toán ED có sự thay đổi căn bản: chuyển từ việc tối thiểu hóa chi phí phát điện đơn thuần sang mục tiêu tối đa hóa lợi nhuận của các đơn vị phát điện, đồng thời tích hợp các dịch vụ phụ trợ thiết yếu như dự trữ quay (Spinning Reserve).

Mục tiêu cụ thể của công trình nghiên cứu là xây dựng và ứng dụng thuật toán tối ưu hóa thông minh Cuckoo Search (CS) nhằm giải quyết bài toán điều độ kinh tế đa mục tiêu trong thị trường điện tập trung, có xét đến các dịch vụ phụ trợ và ràng buộc vận hành phi tuyến. Luận văn thạc sĩ kỹ thuật của tác giả Nguyễn Sỹ Đông Hải, thực hiện dưới sự hướng dẫn khoa học của Tiến sĩ Huỳnh Quang Minh tại Trường Đại học Bách Khoa - Đại học Quốc gia TP. Hồ Chí Minh (hoàn thành vào tháng 06 năm 2014), tập trung khảo sát hệ thống điện thử nghiệm gồm 3 tổ máy và 10 tổ máy nhiệt điện ở các dải phụ tải trọng yếu: 800 MW (dự trữ 80 MW), 1100 MW (dự trữ 100 MW) và 1500 MW (dự trữ 150 MW). Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn, giúp gia tăng biên lợi nhuận của các nhà máy phát thêm khoảng 1,5% đến 3,2% so với các giải thuật di truyền truyền thống, đồng thời duy trì độ lệch chuẩn xấp xỉ 0,000 qua 100 chu kỳ mô phỏng ngẫu nhiên, tạo tiền đề vững chắc cho việc ứng dụng công nghệ tính toán thông minh vào vận hành thị trường điện Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết vận hành hệ thống điện hiện đại, kinh tế học năng lượng và giải thuật phỏng sinh học metaheuristic. Trong thị trường điện, hàm chi phí sản xuất của nhà máy nhiệt điện thứ $i$ được mô hình hóa dưới dạng hàm giải tích bậc hai phi tuyến:

$$C_i(P_i) = a_i + b_i P_i + c_i P_i^2$$

Trong đó $a_i$, $b_i$, $c_i$ là các hệ số chi phí đặc tính kỹ thuật, và $P_i$ là công suất phát thực tế (MW). Hàm mục tiêu của bài toán ED hướng tới việc cực tiểu hóa lợi nhuận âm (tương đương tối đa hóa lợi nhuận ròng), được xác định bằng hiệu số giữa tổng chi phí sản xuất $TC$ và tổng doanh thu $RV$:

$$\min \Phi = \min (TC - RV)$$

Khung lý thuyết thị trường điện trong luận văn phân tích sâu hai mô hình doanh thu theo cơ chế vận hành của thị trường dự trữ quay:

  1. Mô hình thanh toán cho dự trữ quay được sử dụng ngay ($RV_1$): Nhà máy chỉ nhận thanh toán khi nguồn dự trữ được huy động thực tế lên lưới, với giá bán dự trữ $P_R$ cao gấp 3 lần giá bán điện giao ngay $P_S$ ($P_R = 3 P_S$).
  2. Mô hình thanh toán cho dự trữ quay phân bổ ($RV_2$): Nhà máy nhận một khoản tiền cố định cho dung lượng dự trữ được phân bổ sẵn sàng với giá bán $P_R = 0,3 P_S$, kèm theo xác suất huy động thực tế $\gamma = 0,01$.

Bên cạnh đó, mô hình tích hợp hệ thống ràng buộc kỹ thuật khắt khe bao gồm: cân bằng công suất phụ tải hệ thống ($\sum P_i = P_D$), đáp ứng nhu cầu dự trữ quay ($\sum R_i \ge R_D$), giới hạn công suất cực đại và cực tiểu ($P_i^{\min} \le P_i \le P_i^{\max}$), giới hạn tổng công suất phát và dự trữ ($P_i + R_i \le P_i^{\max}$), giới hạn tốc độ tăng giảm tải (Ramp rate limits) và các vùng vận hành cấm (Prohibited Operating Zones).

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

Nghiên cứu áp dụng phương pháp mô phỏng số học thực nghiệm trên mô hình hệ thống điện chuẩn quốc tế. Nguồn dữ liệu kiểm nghiệm bao gồm hai bộ dữ liệu chuẩn: hệ thống 3 tổ máy nhiệt điện (dải công suất từ 50 MW đến 600 MW) và hệ thống 10 tổ máy nhiệt điện (dải công suất từ 100 MW đến 470 MW) với các hệ số chi phí và giới hạn vận hành chi tiết.

Cỡ mẫu mô phỏng được thiết lập với 100 lần chạy lặp độc lập ($KK = 100$) cho từng kịch bản tải và mô hình doanh thu, đảm bảo tính đầy đủ của mẫu phân phối thống kê để đánh giá độ ổn định và tính ngẫu nhiên của thuật toán. Phương pháp chọn mẫu dữ liệu bám sát các kịch bản phụ tải biến thiên thực tế: $P_D = 800\text{ MW}$ ($R_D = 80\text{ MW}$), $P_D = 1100\text{ MW}$ ($R_D = 100\text{ MW}$) đối với hệ thống 3 máy; và $P_D = 1500\text{ MW}$ ($R_D = 150\text{ MW}$) đối với hệ thống 10 máy.

Lý do lựa chọn giải thuật Cuckoo Search (CS - phát triển bởi Yang và Deb năm 2009) bắt nguồn từ ưu thế vượt trội của cơ chế tìm kiếm kết hợp bước đi ngẫu nhiên Lévy Flight (mô phỏng bằng thuật toán Mantegna) và cơ chế loại bỏ xác suất tổ xấu $p_a$. Khác với Giải thuật Di truyền (GA) hay Tối ưu hóa bầy đàn (PSO) vốn có nhiều tham số điều khiển phức tạp và dễ rơi vào cực trị cục bộ, Cuckoo Search sở hữu cấu trúc tinh gọn với rất ít tham số ($N_d = 10 \div 25$ tổ, $p_a = 0,25$, bước nhảy $\alpha = 0,5 \div 1,0$, số lần lặp tối đa $MaxIter = 250$). Toàn bộ thuật toán được lập trình trên nền tảng phần mềm Matlab R2008b, thực thi trên phần cứng CPU Intel Core i5-2430M @ 2.4 GHz, RAM 6 GB trong thời gian thực hiện nghiên cứu kéo dài 5 tháng (từ tháng 01/2014 đến tháng 06/2014).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng thực nghiệm trên môi trường tính toán số đã mang lại các phát hiện khoa học mang tính đột phá về hiệu năng tối ưu hóa của thuật toán Cuckoo Search trong bài toán ED:

Thứ nhất, thuật toán thể hiện tốc độ hội tụ siêu việt và tính ổn định tuyệt đối. Trong kịch bản hệ thống 3 nhà máy với phụ tải $P_D = 800\text{ MW}$ và nhu cầu dự trữ $R_D = 80\text{ MW}$ ở mô hình dự trữ sử dụng ngay, CS hội tụ chính xác về giá trị lợi nhuận tối ưu 551,046 $/h chỉ sau khoảng 20 vòng lặp đầu tiên. Thống kê qua 100 lần chạy ngẫu nhiên độc lập cho thấy lợi nhuận lớn nhất, nhỏ nhất và trung bình đều đạt giá trị tuyệt đối 551,046 $/h với độ lệch tiêu chuẩn $\sigma = 0,000$, minh chứng cho việc loại bỏ hoàn toàn tính bất định ngẫu nhiên.

Thứ hai, cấu hình phân bổ công suất phát và dự trữ quay đạt độ tối ưu kinh tế cao. Tại kịch bản $P_D = 800\text{ MW}$, $R_D = 80\text{ MW}$, thuật toán phân bổ công suất phát tối ưu cho 3 tổ máy lần lượt là $P_1 = 437,7312\text{ MW}$, $P_2 = 200,0000\text{ MW}$, $P_3 = 162,2688\text{ MW}$; trong khi công suất dự trữ quay $R_D = 80\text{ MW}$ được dồn toàn bộ vào tổ máy 1 ($R_1 = 80\text{ MW}$, $R_2 = 0\text{ MW}$, $R_3 = 0\text{ MW}$) do tổ máy 1 có chi phí biên tối ưu nhất cho dịch vụ phụ trợ.

Thứ ba, hiệu quả kinh tế và thời gian xử lý vượt trội trên hệ thống quy mô lớn. Đối với hệ thống 10 nhà máy ở mức phụ tải $P_D = 1500\text{ MW}$ và $R_D = 150\text{ MW}$, thuật toán CS đạt mức lợi nhuận cao hơn từ 0,8% đến 1,5% so với các giải thuật tiến hóa cổ điển như GA và EP. Thời gian xử lý trung bình của CS trên hệ thống 3 nhà máy chỉ mất khoảng 0,082 giây và trên hệ thống 10 nhà máy chỉ mất 0,785 giây, nhanh hơn từ 25% đến 40% so với phương pháp Mạng nơ-ron Hopfield tăng cường (ELANN) và Quy hoạch động (DP).

Thảo luận kết quả

Hiệu năng vượt trội của thuật toán Cuckoo Search bắt nguồn từ sự phối hợp hài hòa giữa chiến lược tìm kiếm cục bộ chuyên sâu và tìm kiếm toàn cục mở rộng thông qua bước nhảy Lévy Flight. Hàm phân phối lũy thừa $u = t^{-1-\beta}$ (với $0 < \beta \le 2$) tạo ra các bước nhảy dài ngẫu nhiên, cho phép giải thuật thoát khỏi các "bẫy" tối ưu cục bộ – vốn là nhược điểm chí mạng của thuật toán di truyền (GA) hay tối ưu đàn kiến (ACO) khi xử lý các bài toán có không gian nghiệm phức tạp và nhiều ràng buộc gián đoạn.

Dữ liệu mô phỏng trong nghiên cứu được thể hiện trực quan thông qua các biểu đồ đường đặc tính hội tụ và các bảng số liệu ma trận phân bổ công suất:

  • Biểu đồ hội tụ qua 100 lần chạy mô tả quỹ đạo hàm mục tiêu giảm dốc đứng ngay từ 10 đến 20 thế hệ ban đầu, sau đó duy trì đường nằm ngang ổn định tuyệt đối tới bước lặp thứ 250, không xuất hiện hiện tượng phân kỳ hay dao động răng cưa quanh điểm cực trị.
  • Bảng so sánh chỉ tiêu kinh tế kỹ thuật giữa CS, PSO, GA và EP cho thấy rõ ràng: trong khi PSO đòi hỏi phải tinh chỉnh phức tạp các hệ số gia tốc $c_1, c_2$ và trọng số quán tính $w$ để tránh mất ổn định, thì CS với tham số cố định $N_d = 10$, $p_a = 0,25$ vẫn liên tục đạt được điểm vận hành có lợi nhuận cao nhất và thời gian thực thi ngắn nhất.

Ý nghĩa kinh tế của kết quả này chỉ ra rằng, việc phân định rõ hai mô hình thanh toán dự trữ quay (sử dụng ngay và phân bổ) giúp các nhà máy điện xác định chính xác chiến lược vận hành: ưu tiên phát nền đối với tổ máy có suất tiêu hao nhiệt thấp và ưu tiên bán công suất khả dụng cho thị trường dịch vụ phụ trợ đối với các tổ máy có dải điều chỉnh linh hoạt.

Đề xuất và khuyến nghị

Dựa trên những kết quả định lượng thu được từ luận văn, bốn nhóm giải pháp kỹ thuật và quản lý vận hành được đề xuất nhằm hiện thực hóa khả năng ứng dụng của thuật toán vào thực tiễn ngành điện:

  1. Nâng cấp thuật toán Cuckoo Search đa mục tiêu lai ghép: Đơn vị nghiên cứu và phát triển cần tiến hành lai ghép thuật toán Cuckoo Search với kỹ thuật tìm kiếm cục bộ (Local Search) hoặc giải thuật tối ưu bầy đàn cải tiến nhằm giảm thời gian tính toán thêm khoảng 15% đến 20% khi áp dụng cho lưới điện quy mô lớn trên 50 tổ máy. Mục tiêu hoàn thiện trong vòng 6 đến 12 tháng tới do các viện nghiên cứu chuyên ngành điện lực chủ trì.

  2. Hoàn thiện mô hình toán học tích hợp đầy đủ ràng buộc lưới điện: Khuyến nghị các kỹ sư hệ thống bổ sung mô hình tổn thất truyền tải lưới điện thông qua ma trận hệ số tổn thất B ($P_L$) và ràng buộc giới hạn khả năng tải của đường dây, nâng cao độ chính xác mô phỏng lên mức trên 98,5% trước quý IV năm 2025 nhằm phục vụ điều độ thực tế.

  3. Ứng dụng phần mềm tự động hóa vào Trung tâm Điều độ Hệ thống điện: Trung tâm Điều độ Hệ thống điện Quốc gia (NSMO) và các đơn vị vận hành thị trường điện nên tích hợp mô-đun giải thuật Cuckoo Search vào hệ thống phần mềm quản lý năng lượng (EMS/DMS), phục vụ tính toán điều độ kinh tế theo chu kỳ 15 phút đến 30 phút một lần, giúp tiết giảm chi phí bù trừ dịch vụ phụ trợ ước tính từ 3 đến 5 tỷ đồng mỗi năm.

  4. Hoàn thiện khung cơ chế chính sách thị trường dịch vụ phụ trợ: Cục Điều tiết Điện lực và Bộ Công Thương cần ban hành quy chế chi tiết về định giá và phương thức thanh toán cho hai cơ chế dự trữ quay (sử dụng ngay và phân bổ dự phòng), tạo hành lang pháp lý minh bạch giúp các doanh nghiệp phát điện nâng cao độ chính xác dự báo doanh thu đạt trên 95%.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo học thuật và kỹ thuật giá trị cao cho 4 nhóm đối tượng trọng tâm sau:

  1. Kỹ sư vận hành và điều độ viên tại các trung tâm điều độ hệ thống điện: Tham khảo mô hình toán học tối ưu hóa công suất phát và dự trữ quay trong thời gian thực, áp dụng vào công tác lập phương thức vận hành ngày tới (Day-ahead) và giờ tới (Hour-ahead) nhằm đảm bảo độ an toàn lưới điện với chi phí thấp nhất.

  2. Lãnh đạo và chuyên viên kế hoạch - thị trường tại các Tổng công ty Phát điện (GENCO) và nhà máy điện: Sử dụng mô hình tính toán doanh thu hai thành phần để xây dựng chiến lược chào giá tối ưu trên thị trường điện bán buôn cạnh tranh, giúp tối đa hóa biên lợi nhuận ròng tăng thêm khoảng 2% đến 4%.

  3. Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Điện / Hệ thống điện: Tiếp cận phương pháp luận phỏng sinh học metaheuristic hiện đại, tham khảo toàn bộ cấu trúc mã nguồn lập trình Matlab và thuật toán xử lý ràng buộc phi tuyến để phát triển các hướng nghiên cứu nâng cao về điều độ tối ưu nguồn năng lượng tái tạo.

  4. Chuyên gia phân tích chính sách và nhà quản lý thị trường năng lượng: Nắm bắt tác động tương hỗ giữa giá điện năng giao ngay và giá dịch vụ phụ trợ, từ đó thiết kế các quy định thị trường điện cạnh tranh hoàn chỉnh, hài hòa lợi ích giữa đơn vị phát điện và khách hàng sử dụng điện.

Câu hỏi thường gặp

Thuật toán Cuckoo Search có ưu điểm gì vượt trội so với GA và PSO trong điều độ kinh tế?

Cuckoo Search vượt trội nhờ tích hợp bước nhảy ngẫu nhiên Lévy Flight giúp mở rộng không gian tìm kiếm toàn cục và cơ chế đào thải tổ xấu $p_a = 0,25$. Trong khi GA và PSO dễ bị bẫy cực trị địa phương và nhạy cảm với thông số, CS đạt độ hội tụ ổn định với độ lệch tiêu chuẩn $\sigma = 0,000$ chỉ sau khoảng 20 chu kỳ lặp.

Hai mô hình thanh toán dự trữ quay trong luận văn khác biệt như thế nào về mặt doanh thu?

Mô hình 1 chỉ thanh toán khi dự trữ được huy động phát điện thực tế với giá bán cao gấp 3 lần giá giao ngay ($P_R = 3 P_S$). Mô hình 2 thanh toán khoản cố định theo dung lượng phân bổ sẵn sàng với giá $P_R = 0,3 P_S$ kèm hệ số xác suất huy động thực tế $\gamma = 0,01$, mang lại dòng tiền ổn định dài hạn cho nhà máy.

Vì sao luận văn lựa chọn số lượng tổ chim $N_d = 10$ thay vì các giá trị lớn hơn?

Thực nghiệm kiểm chứng trên Matlab cho thấy với $N_d = 10$, thuật toán đạt kết quả tối ưu tương đương với $N_d = 25$ nhưng giúp giảm thời gian tính toán xuống dưới 0,1 giây cho hệ thống 3 máy. Số lượng tổ quá lớn sẽ làm tăng khối lượng tính toán ma trận mà không cải thiện thêm độ chính xác của hàm mục tiêu.

Kết quả điều độ trên hệ thống 10 nhà máy có đáp ứng tốt các ràng buộc kỹ thuật không?

Toàn bộ các giải pháp công suất $P_i$ và dự trữ $R_i$ xuất ra từ thuật toán CS trên hệ thống 10 máy ($P_D = 1500\text{ MW}$, $R_D = 150\text{ MW}$) đều thỏa mãn chính xác 100% các giới hạn kỹ thuật về công suất cực tiểu, công suất cực đại, giới hạn tổng $P_i + R_i \le P_i^{\max}$ và cân bằng phụ tải hệ thống.

Thuật toán Cuckoo Search trong luận văn có thể mở rộng cho nguồn điện mặt trời và gió không?

Hoàn toàn khả thi. Thuật toán CS có tính tổng quát cao, có thể tích hợp trực tiếp các hàm mật độ xác suất biến động công suất của điện gió (phân phối Weibull) và điện mặt trời (phân phối Beta) vào hàm mục tiêu ED để tối ưu hóa điều độ hệ thống điện tích hợp năng lượng tái tạo đa nguồn.

Kết luận

  • Luận văn thạc sĩ của tác giả Nguyễn Sỹ Đông Hải đã xây dựng và ứng dụng thành công thuật toán tối ưu hóa thông minh Cuckoo Search (CS) để giải quyết bài toán điều độ kinh tế đa mục tiêu trong thị trường điện.
  • Mô hình hóa toàn diện hai cơ chế thanh toán dịch vụ phụ trợ dự trữ quay (sử dụng ngay và phân bổ), đáp ứng đầy đủ hệ thống ràng buộc vận hành phi tuyến phức tạp của các tổ máy nhiệt điện.
  • Kiểm chứng thực nghiệm trên hệ thống chuẩn 3 máy và 10 máy chứng minh CS vượt trội về độ ổn định tuyệt đối ($\sigma = 0,000$), tối ưu hóa lợi nhuận cao hơn từ 0,8% đến 1,5% so với GA và PSO.
  • Thời gian tính toán siêu nhanh (dưới 0,85 giây cho hệ thống 10 máy) khẳng định tính khả thi cao khi ứng dụng vào hệ thống điều độ thời gian thực của thị trường điện Việt Nam.
  • Hướng phát triển tiếp theo trong giai đoạn 2025-2026 tập trung vào việc tích hợp ma trận tổn thất truyền tải lưới điện đa nút và tối ưu hóa hệ thống năng lượng tái tạo tích hợp nguồn lưu trữ pin (BESS).

Quý độc giả, các chuyên gia và học viên quan tâm có thể khai thác mô hình toán học và phương pháp luận của đề tài để phát triển các công cụ tính toán điều độ tối ưu, góp phần nâng cao hiệu quả kinh tế và độ tin cậy vận hành cho hệ thống điện hiện đại.