Tổng quan nghiên cứu

Trong bối cảnh thị trường năng lượng hiện đại đang chuyển dịch mạnh mẽ theo hướng phát thải thấp, bài toán điều tiết sản lượng điện nhằm cân bằng lợi nhuận doanh nghiệp và giảm thiểu ô nhiễm môi trường trở thành thách thức cấp thiết. Luận văn thạc sĩ chuyên ngành Toán ứng dụng với mã số 8 46 01 12, thực hiện tại Trường Đại học Khoa học thuộc Đại học Thái Nguyên năm 2019 dưới sự hướng dẫn khoa học của GS. Lê Dũng Mưu, đã giải quyết thấu đáo bài toán tối ưu hóa phức hợp này. Nghiên cứu tập trung vào vấn đề cạnh tranh sản lượng bất hợp tác giữa các doanh nghiệp phát điện theo mô hình Nash - Cournot, đồng thời tích hợp thêm ràng buộc chung về giới hạn phí xả thải môi trường.

Mục tiêu cụ thể của luận văn là phát triển một tiếp cận mới dựa trên lý thuyết bài toán cân bằng tách (Split Equilibrium Problem - SEP), từ đó xây dựng một thuật toán lặp xấp xỉ không gian hiệu quả kèm chứng minh toán học chặt chẽ về sự hội tụ nghiệm. Phạm vi nghiên cứu được xác lập trong không gian Hilbert thực hữu hạn chiều, mô phỏng thị trường điện lực gồm 2 công ty với 4 tổ máy phát điện vận hành trong dải công suất chiến lược từ 0 đến 100 đơn vị sản lượng trong giai đoạn 2018 - 2019. Về mặt ý nghĩa thực tiễn và định lượng, phương pháp đề xuất giúp các nhà quản lý thị trường kiểm soát 100% các tiêu chuẩn phát thải carbon và tro xỉ than, đồng thời tối ưu hóa lợi nhuận toàn ngành với mức giảm thiểu chi phí xử lý môi trường ước tính đạt khoảng 18,5% so với các mô hình vận hành tự do truyền thống.

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

Khung lý thuyết áp dụng

Nền tảng học thuật của đề tài được xây dựng trên sự giao thoa giữa Lý thuyết trò chơi bất hợp tác và Giải tích lồi hiện đại. Thứ nhất, mô hình cân bằng Nash - Cournot nguyên bản từ nghiên cứu của Nikaido và Isoda năm 1955 được sử dụng để mô tả trạng thái cân bằng cạnh tranh, nơi không một doanh nghiệp nào có thể đơn phương thay đổi sản lượng phát điện để gia tăng lợi nhuận. Thứ hai, lý thuyết bài toán chấp nhận tách (Split Feasibility Problem) do Censor và Elving khởi xướng năm 1994 được mở rộng thành bài toán cân bằng tách kết hợp quy hoạch lồi.

Khung nghiên cứu vận hành dựa trên 4 khái niệm cốt lõi:

  1. Không gian tập lồi đa diện trong không gian Euclid n chiều biểu diễn miền chiến lược của các doanh nghiệp.
  2. Toán tử gần kề (Proximal Operator) của hàm chi phí môi trường với tham số điều khiển dương.
  3. Ánh xạ đơn điệu và giả đơn điệu (Pseudomonotone) đảm bảo tính ổn định của tập nghiệm cân bằng.
  4. Khái niệm dưới vi phân xấp xỉ (Epsilon-subdifferential) cho phép tính toán hướng dịch chuyển tối ưu mà không đòi hỏi hàm mục tiêu phải khả vi trơn.

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

Nghiên cứu sử dụng phương pháp mô phỏng số học kết hợp giải tích toán học thuần túy trên tập dữ liệu thực nghiệm gồm 2 tổng công ty phát điện quản lý 4 nhà máy nhiệt điện than. Cỡ mẫu nghiên cứu được thiết lập bao gồm 4 biến quyết định sản lượng tương ứng với 4 tổ máy, vận hành trong không gian Euclid 4 chiều với miền ràng buộc hộp từ 0 đến 100 đơn vị công suất.

Phương pháp chọn mẫu có chủ đích đại diện cho cấu trúc thị trường bán buôn điện cạnh tranh dạng thiểu quyền (Oligopoly), trong đó công ty thứ nhất sở hữu 3 nhà máy và công ty thứ hai sở hữu 1 nhà máy. Lý do lựa chọn phương pháp phân tích giải tích biến phân và thuật toán phép chiếu lặp là vì mô hình thực tế chứa đựng các ràng buộc phi tuyến liên kết giữa hai không gian tín hiệu và mục tiêu, khiến các giải thuật quy hoạch tuyến tính cổ điển hoàn toàn bất khả thi. Toàn bộ quy trình giải tích, thiết lập dãy lặp và kiểm chứng điều kiện hội tụ được thực hiện xuyên suốt trong khoảng thời gian 6 tháng từ tháng 10 năm 2018 đến tháng 04 năm 2019.

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

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

Quá trình phân tích giải tích và thực nghiệm thuật toán trên mô hình thị trường điện lực đã mang lại 4 phát hiện khoa học quan trọng:

  1. Thuật toán cân bằng tách đề xuất đạt tốc độ hội tụ cao, dãy nghiệm xấp xỉ tiến về nghiệm tối ưu sau khoảng 25 đến 35 bước lặp với sai số chuẩn Euclid tuyệt đối nhỏ hơn 0,0001, khẳng định tính khả thi vượt trội khi tính toán thực tế.
  2. Việc áp dụng ràng buộc chi phí môi trường thông qua toán tử tuyến tính chuyển vị đã giúp cắt giảm khoảng 18,5% tổng lượng phát thải than tiêu thụ tại 4 nhà máy phát điện so với kịch bản cân bằng Nash thông thường không có ràng buộc sinh thái.
  3. Với thông số giá trần thị trường bằng 10 và hệ số nhạy cảm giá bằng 0,1, sản lượng phát điện cân bằng tối ưu của công ty thứ nhất (sở hữu 3 nhà máy) chiếm tỷ trọng 67,4% tổng cung, trong khi công ty thứ hai (1 nhà máy) cung ứng 32,6% còn lại của thị trường.
  4. Nghiên cứu đã chứng minh thành công 4 bổ đề kỹ thuật và 1 định lý hội tụ toàn cục, chỉ ra rằng dãy nghiệm luôn bị chặn và hội tụ mạnh về nghiệm tối ưu của bài toán cân bằng lồi tách dưới các giả thiết đơn điệu và liên tục Lipschitz của toán tử gradient.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp thuật toán cân bằng tách vận hành ổn định là nhờ sự kết hợp nhịp nhàng giữa phép chiếu trực giao lên tập ràng buộc chiến lược và toán tử gần kề của hàm phạt ô nhiễm. Khi giải quyết bài toán tối ưu phi tuyến nhiều biến, cách tiếp cận này đã tách rời bài toán tìm điểm cân bằng thị trường phức tạp thành các bài toán con đơn giản hơn tại từng bước lặp.

So sánh với các nghiên cứu tiền tiêu của Santos và Scheimberg năm 2011 cũng như công trình của Lê Thị Hải Yến và GS. Lê Dũng Mưu năm 2016, thuật toán trong luận văn này có ưu thế nổi bật khi sử dụng kỹ thuật dưới vi phân xấp xỉ. Cải tiến này giúp loại bỏ hoàn toàn yêu cầu tính toán đạo hàm cấp một chính xác tại mỗi vòng lặp, từ đó giảm khoảng 28% thời gian xử lý của bộ vi xử lý máy tính.

Dữ liệu thực nghiệm của mô hình có thể được trực quan hóa tối ưu thông qua một bảng ma trận công suất biểu diễn phân bổ tải chi tiết giữa 4 tổ máy và một biểu đồ đường suy giảm logarit thể hiện tốc độ triệt tiêu sai số theo từng bước lặp. Về mặt ý nghĩa kinh tế, mô hình cung cấp một công cụ định lượng chuẩn xác giúp cơ quan quản lý năng lượng điều tiết thị trường điện cạnh tranh một cách minh bạch, hài hòa lợi ích kinh tế của doanh nghiệp với mục tiêu bảo vệ môi trường sống.

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

Dựa trên các kết quả lý thuyết và thực nghiệm vững chắc, luận văn đưa ra 4 nhóm giải pháp mang tính ứng dụng cao cho ngành năng lượng:

  1. Tích hợp giải thuật cân bằng tách vào hệ thống phần mềm điều độ thị trường điện lực: Giao Cục Điều tiết Điện lực chủ trì triển khai thử nghiệm thuật toán vào phần mềm quản lý thị trường phát điện trong vòng 12 đến 18 tháng tới, nhằm mục tiêu tối ưu hóa công suất phát và giảm thiểu ít nhất 15% tổng phát thải tro xỉ than trên toàn hệ thống lưới điện quốc gia.
  2. Xây dựng khung chính sách thuế môi trường lũy tiến theo hàm toán tử tuyến tính: Bộ Tài chính phối hợp cùng Bộ Tài nguyên và Môi trường áp dụng mô hình hàm phí ô nhiễm trong giai đoạn 2020 - 2022, thiết lập mức phí linh hoạt gắn liền với sản lượng tiêu thụ nhiên liệu của từng nhà máy để thu hồi khoảng 20% chi phí ngoại ứng tiêu cực do ô nhiễm nhiệt điện gây ra.
  3. Hiện đại hóa quy trình lập kế hoạch sản xuất tại các doanh nghiệp phát điện: Các tập đoàn và tổng công ty phát điện như EVN Genco 1, Genco 2 và Genco 3 cần ứng dụng công thức tính toán dưới vi phân xấp xỉ trong vòng 6 đến 9 tháng để tái cơ cấu định mức tiêu hao nhiên liệu tại 100% các tổ máy, giúp cải thiện biên lợi nhuận ròng từ 8% đến 12%.
  4. Mở rộng đề tài nghiên cứu sang các nguồn năng lượng tái tạo bất định: Các viện nghiên cứu toán học và năng lượng tiếp tục phát triển mô hình cân bằng tách cho hệ thống tích hợp điện gió và điện mặt trời trong lộ trình 24 tháng, hướng tới kiểm soát độ biến thiên công suất và đảm bảo mục tiêu 30% tỷ trọng năng lượng sạch vào năm 2030.

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

Công trình luận văn thạc sĩ này là nguồn tư liệu học thuật và ứng dụng chuyên sâu dành cho 4 nhóm đối tượng trọng tâm:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Toán ứng dụng, Tối ưu hóa và Giải tích biến phân: Nghiên cứu cung cấp tài liệu tham khảo mẫu mực về kỹ thuật chứng minh sự hội tụ của thuật toán chiếu tách và khai thác toán tử gần kề trong không gian đa chiều.
  2. Kỹ sư vận hành hệ thống điện và chuyên viên trung tâm điều độ hệ thống điện quốc gia: Tài liệu hỗ trợ trực tiếp việc thiết lập thuật toán lập lịch huy động công suất các nhà máy điện tối ưu và phân bổ dòng công suất trên lưới điện truyền tải.
  3. Chuyên viên hoạch định chính sách tại các cơ quan quản lý năng lượng và môi trường: Hỗ trợ các nhà quản lý xây dựng công cụ mô phỏng định lượng để dự báo chính xác hành vi phát điện của các tập đoàn khi áp dụng các mức thuế bảo vệ môi trường mới.
  4. Giảng viên đại học giảng dạy môn Lý thuyết trò chơi, Quy hoạch phi tuyến và Kinh tế lượng: Sử dụng các ví dụ thực tế về thị trường điện 4 nhà máy trong luận văn làm bài giảng tình huống sinh động cho sinh viên và học viên sau đại học.

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

1. Bài toán cân bằng tách khác biệt như thế nào so với bài toán cân bằng Nash truyền thống? Mô hình cân bằng Nash truyền thống chỉ tìm kiếm điểm chiến lược tối ưu nội tại giữa các doanh nghiệp độc lập. Bài toán cân bằng tách mở rộng thêm một không gian mục tiêu mới, liên kết thông qua một toán tử tuyến tính để bắt buộc điểm cân bằng vừa phải tối ưu hóa cạnh tranh thị trường, vừa phải thỏa mãn ràng buộc tối thiểu hóa một hàm chi phí phụ như phí môi trường.

2. Tại sao luận văn lại sử dụng kỹ thuật dưới vi phân xấp xỉ thay vì đạo hàm giải tích thông thường? Trong thực tế sản xuất điện, hàm lợi nhuận và hàm chi phí môi trường thường có tính chất không khả vi trơn hoặc chứa điểm gãy góc tại các ngưỡng công suất giới hạn. Kỹ thuật dưới vi phân xấp xỉ cho phép thuật toán tính toán được các vector hướng giảm mà không cần hàm số có đạo hàm, giúp tăng tính ổn định và tránh hiện tượng phân kỳ toán học.

3. Thuật toán lặp của luận văn xử lý bài toán thị trường điện 4 tổ máy như thế nào? Thuật toán phân tách 4 tổ máy thành các biến độc lập thuộc miền ràng buộc hộp từ 0 đến 100 đơn vị công suất. Tại mỗi bước lặp, thuật toán tính toán giá điện theo hàm tuyến tính, cập nhật vector sản lượng qua toán tử chiếu trực giao và điều chỉnh độ lệch bằng toán tử gần kề của phí ô nhiễm cho đến khi sai số đạt mức dưới 0,0001.

4. Những điều kiện toán học nào bảo đảm thuật toán luôn luôn hội tụ về nghiệm tối ưu? Thuật toán được chứng minh hội tụ tuyệt đối khi thỏa mãn 4 giả thiết: hàm cân bằng là lồi nửa liên tục dưới, tập dưới vi phân xấp xỉ khác rỗng và bị chặn, hàm thỏa mãn tính chất giả đơn điệu trên tập nghiệm và toán tử gradient của hàm khoảng cách liên tục Lipschitz với hằng số chuẩn bị chặn.

5. Doanh nghiệp nhiệt điện than có thể áp dụng ngay mô hình này để giảm chi phí hay không? Doanh nghiệp hoàn toàn có thể ứng dụng ngay mô hình bằng cách nạp các thông số kỹ thuật thực tế gồm định mức tiêu hao than cho một đơn vị điện năng và biểu phí xả thải môi trường vào ma trận bài toán. Việc này giúp doanh nghiệp tự động tính toán được sản lượng phát điện tối ưu hàng ngày để đạt lợi nhuận tối đa sau thuế.

Kết luận

Bản luận văn thạc sĩ của tác giả Nguyễn Thành Huế đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn nổi bật:

  • Hệ thống hóa toàn diện các kiến thức nền tảng về giải tích lồi, phép chiếu và toán tử đơn điệu trong không gian Euclid hữu hạn chiều.
  • Thiết lập thành công mô hình toán học tích hợp bài toán cân bằng Nash - Cournot với bài toán chấp nhận lồi tách có ràng buộc môi trường.
  • Xây dựng thuật toán lặp xấp xỉ hiệu quả sử dụng kỹ thuật dưới vi phân và toán tử gần kề không yêu cầu đạo hàm trơn.
  • Cung cấp các chứng minh giải tích toán học chuẩn xác về tính bị chặn và định lý hội tụ toàn cục của dãy nghiệm lặp.
  • Mô phỏng thành công bài toán thực tế trên thị trường điện năng 4 nhà máy, mở ra giải pháp định lượng mẫu mực cho nền kinh tế năng lượng xanh.

Trong giai đoạn 2020 - 2025, hướng phát triển tiếp theo của đề tài là mở rộng thuật toán lên không gian Hilbert vô hạn chiều và tích hợp các yếu tố ngẫu nhiên của thị trường năng lượng tái tạo. Đây là công trình học thuật có giá trị khoa học cao, đóng góp thiết thực vào kho tàng nghiên cứu Toán ứng dụng và tiến trình chuyển đổi năng lượng bền vững tại Việt Nam.