Tổng quan nghiên cứu

Thủy điện đóng vai trò trụ cột trong việc đảm bảo an ninh năng lượng quốc gia và thúc đẩy tăng trưởng kinh tế bền vững. Nhà máy thủy điện Đại Ninh nằm trên địa bàn hai tỉnh Lâm Đồng và Bình Thuận, sở hữu tổng công suất lắp đặt 300 MW gồm hai tổ máy phát điện 150 MW, từng là công trình thủy điện có quy mô công suất đứng thứ 5 tại Việt Nam khi đi vào vận hành từ năm 2008. Công trình có mực nước chết 860 m ứng với dung tích chết 68 triệu m3 và mực nước dâng bình thường 880 m tương ứng dung tích hồ chứa 320 triệu m3.

Mặc dù chỉ tiêu sản lượng điện năng được giao hàng năm đạt mức 1,2 tỷ kWh, nhưng trong 7 năm vận hành từ năm 2008 đến 2014, nhà máy chỉ có 2 năm đạt hoặc vượt chỉ tiêu thiết kế. Điển hình như năm 2011 sản lượng chỉ đạt hơn 905 triệu kWh và năm 2014 đạt gần 1,08 tỷ kWh. Nguyên nhân chủ yếu xuất phát từ phương thức vận hành hồ chứa theo kinh nghiệm và biểu đồ điều phối truyền thống, thiên về tính an toàn cơ học mà chưa được tối ưu hóa toán học theo biến động thủy văn thực tế.

Mục tiêu trọng tâm của nghiên cứu là xây dựng mô hình toán học tối ưu hóa công suất phát điện và doanh thu của nhà máy thủy điện Đại Ninh trong cả ngắn hạn lẫn dài hạn, đồng thời ứng dụng lý thuyết trò chơi để thiết lập chiến lược chào giá cạnh tranh trên thị trường điện. Kết quả nghiên cứu chứng minh phương án điều tiết tối ưu giúp nâng cao doanh thu và sản lượng điện tăng thêm 8,6% trong giai đoạn ngắn hạn, đồng thời triệt tiêu hoàn toàn lượng nước xả thừa qua đập tràn, góp phần nâng cao hiệu quả sử dụng tài nguyên nước và tối ưu hóa chi phí toàn hệ thống điện quốc gia.

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 tích hợp hai trụ cột toán học và kinh tế học hiện đại gồm lý thuyết tối ưu hóa phi tuyến và lý thuyết trò chơi trong kinh tế năng lượng. Khung lý thuyết này vận dụng 5 khái niệm nền tảng: bài toán quy hoạch toàn phương, điều kiện tối ưu Karush-Kuhn-Tucker, nón tiếp xúc của tập ràng buộc, trạng thái cân bằng Nash và hàm chi phí biên phát điện.

Bài toán điều tiết thủy điện được mô hình hóa dưới dạng quy hoạch toàn phương có dạng cực tiểu hóa hàm mục tiêu f(x) = (1/2) x^T Q x + c^T x với hệ ràng buộc tuyến tính Ax <= b và Aeq x = beq. Trong đó, ma trận Q là ma trận đối xứng nửa xác định dương, phản ánh đặc tính động lực học của lưu lượng qua turbin và dung tích hữu ích của hồ chứa. Điều kiện tối ưu Karush-Kuhn-Tucker đóng vai trò xác lập nghiệm tối ưu toàn cục thông qua hệ phương trình nhân tử Lagrange và đạo hàm theo hướng.

Trong phân tích thị trường phát điện cạnh tranh, lý thuyết trò chơi không hợp tác dạng chuẩn tắc được ứng dụng để tìm điểm cân bằng Nash giữa các nhà máy điện độc lập. Chi phí phát điện của mỗi tổ máy được mô hình hóa bằng hàm bậc hai Ci(P) = ai P^2 + bi P + ci, với chi phí biên là đạo hàm bậc nhất lambda_i = 2 ai P + bi. Trạng thái cân bằng Nash đạt được khi chi phí biên của tất cả các nhà máy tham gia thị trường bằng nhau, đảm bảo tối đa hóa lợi nhuận đơn vị và giảm thiểu tổng chi phí phát điện toàn hệ thống.

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

Dữ liệu nghiên cứu được trích xuất trực tiếp từ nhật ký vận hành thực tế của Nhà máy thủy điện Đại Ninh và biểu đồ phụ tải hệ thống điện Việt Nam. Cỡ mẫu nghiên cứu bao gồm hai chu kỳ tính toán rõ rệt: mẫu ngắn hạn gồm 480 mốc thời gian theo từng giờ trong 20 ngày đầu tháng 12 năm 2014, và mẫu dài hạn gồm 8.760 mốc thời gian tương ứng 8.760 giờ trong 3 năm thủy văn đặc trưng 2011, 2012 và 2013. Phương pháp chọn mẫu theo chuỗi thời gian liên tục giúp phản ánh đầy đủ chu kỳ thủy văn mùa khô, mùa lũ và biến động phụ tải.

Lý do lựa chọn phương pháp phân tích quy hoạch toàn phương giải bằng thuật toán điểm trong Interior-Point-Convex trên nền tảng phần mềm MATLAB là nhằm giải quyết bài toán quy hoạch phi tuyến cỡ lớn có số chiều ma trận ràng buộc lên tới 2.400 dòng và 960 cột trong bài toán ngắn hạn. Thuật toán này cho phép xử lý đồng thời 5 nhóm ràng buộc khắt khe: lưu lượng xả qua turbin từ 0 đến 55 m3/s, dung tích hồ chứa giới hạn từ 68 triệu m3 đến 320 triệu m3, dòng chảy sinh thái tối thiểu hạ du đạt 0,7 m3/s, và phương trình cân bằng dung tích nước tại từng mốc thời gian với sai số tiệm cận 0.

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

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

Thứ nhất, việc áp dụng mô hình quy hoạch toàn phương cho chu kỳ ngắn hạn 20 ngày đầu tháng 12 năm 2014 giúp tăng trưởng vượt bậc về hiệu quả kinh tế kỹ thuật. Doanh thu thực tế của nhà máy đạt 16,3 tỷ đồng tương ứng 29,7 triệu kWh sản lượng điện với giá bán cố định 550 đồng/kWh. Khi áp dụng giải thuật tối ưu, doanh thu đạt 17,7 tỷ đồng và sản lượng điện đạt 32,2 triệu kWh. Mức tăng trưởng đạt chính xác 8,6% về cả sản lượng lẫn doanh thu, trong khi lượng nước xả lãng phí qua cửa tràn giảm hoàn toàn về mức 0 m3/s.

Thứ hai, nghiên cứu xác lập được quỹ đạo điều khiển lưu lượng nước qua turbin và dung tích hồ chứa mang tính tối ưu cao. Trong 72 giờ đầu tiên, lưu lượng nước qua turbin được điều chỉnh giảm dần từ mức cực đại 55 m3/s xuống 20 m3/s, làm dung tích hồ giảm nhẹ từ 304 triệu m3 xuống 301 triệu m3 nhằm tích lũy cột nước hiệu dụng, sau đó duy trì ổn định mức phát 20 m3/s cho đến hết chu kỳ 480 giờ để đưa dung tích hồ chứa quay trở lại đúng mức ban đầu.

Thứ ba, quy luật vận hành tối ưu dài hạn 8.760 giờ được phân định rõ nét theo ba kịch bản thủy văn. Năm ít nước 2011 với lưu lượng dòng chảy về bình quân 17,8 m3/s chỉ cần vận hành 1 tổ máy ở mức lưu lượng 13 m3/s trong suốt 6 tháng đầu năm để duy tu tổ máy còn lại. Trong khi đó, các năm có lượng nước trung bình như năm 2012 với lưu lượng 28,2 m3/s và năm 2013 với lưu lượng 28,7 m3/s, lưu lượng qua turbin duy trì 20 đến 24 m3/s từ tháng 1 đến tháng 6, tăng dần vào tháng 7 và phát cực đại 40 đến 55 m3/s vào tháng 9 khi lũ chính vụ về.

Thứ tư, khi giải bài toán phân bổ công suất 350 MWh giữa hai nhà máy Đa Nhim và Đại Ninh theo mô hình cân bằng Nash, điểm cân bằng tối ưu xác định mức huy động của Đa Nhim là 157,6 MWh và Đại Ninh là 192,4 MWh. Khi đó, tổng chi phí phát điện của hai nhà máy đạt mức thấp nhất là 251.224.616 đồng, tổng doanh thu đạt 366.024.432 đồng, tạo ra lợi nhuận 53.103.848 đồng cho Đa Nhim và 61.695.968 đồng cho Đại Ninh.

Thảo luận kết quả

Các phát hiện nghiên cứu cho thấy sự khác biệt bản chất giữa phương pháp điều tiết tối ưu toán học và phương thức vận hành thủ công truyền thống. Phương pháp truyền thống thường xả nước phát điện đều đặn theo kinh nghiệm dẫn đến việc hồ chứa đầy sớm và phải mở cửa xả tràn lãng phí trong mùa mưa, hoặc thiếu hụt cột nước đẩy công suất phát giảm sâu trong mùa khô. Mô hình toán học đã giải quyết triệt để sự đánh đổi này thông qua việc quản lý cột nước áp lực cao nhất có thể tại từng thời điểm.

Dữ liệu vận hành được mô tả trực quan qua hệ thống đồ thị quỹ đạo dung tích hồ và biểu đồ phân bổ lưu lượng theo thời gian. So sánh với các nghiên cứu trước đây tại các cụm thủy điện bậc thang Sông Bung sử dụng ngôn ngữ Delphi, giải pháp lập trình trên MATLAB với hàm quadprog xử lý hệ ma trận 5N nhân 2N cho tốc độ hội tụ nhanh hơn và kiểm soát chặt chẽ ràng buộc đẳng thức bảo toàn khối lượng nước.

Đặc biệt, kết quả mô phỏng thị trường điện chỉ ra rằng nếu một nhà máy tự ý nâng giá chào bằng cách tăng hệ số chi phí từ 1,8 lên 2,0, sản lượng được huy động sẽ sụt giảm từ 192,4 MWh xuống 183,5 MWh. Ngược lại, nếu tối ưu hóa kỹ thuật để giảm hệ số chi phí thực tế xuống 1,6, sản lượng huy động tăng lên 202,3 MWh, đưa lợi nhuận của nhà máy Đại Ninh tăng vọt lên 62.195.384 đồng, đồng thời kéo giảm tổng chi phí phát điện toàn hệ thống xuống còn 243.438.231 đồng.

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

Dựa trên kết quả nghiên cứu lý thuyết và kiểm nghiệm số liệu thực nghiệm, bốn nhóm giải pháp hành động cụ thể được đề xuất nhằm nâng cao hiệu quả vận hành và thương mại hóa nguồn điện thủy điện:

Thứ nhất, tích hợp thuật toán tối ưu hóa quy hoạch toàn phương vào hệ thống SCADA điều khiển tự động của nhà máy thủy điện Đại Ninh. Ban Giám đốc và Phòng Kỹ thuật nhà máy cần chủ trì triển khai giải pháp này trong lộ trình 6 tháng, nhằm tự động điều phối lưu lượng xả qua turbin từ 0 đến 55 m3/s theo thời gian thực, đặt mục tiêu triệt tiêu 100% lượng nước xả tràn lãng phí và nâng sản lượng điện phát thêm tối thiểu 8,6%.

Thứ hai, chuẩn hóa quy trình điều tiết hồ chứa phân kỳ theo đặc thù thủy văn từng mùa. Phòng Điều độ vận hành cần áp dụng quy tắc vận hành linh hoạt trong vòng 12 tháng: duy trì lưu lượng 20 đến 24 m3/s trong 6 tháng mùa khô kết hợp lịch bảo dưỡng luân phiên một tổ máy 150 MW, và huy động công suất tối đa 40 đến 55 m3/s trong tháng 9, hướng tới mục tiêu duy trì sản lượng điện thương phẩm ổn định trên 1,2 tỷ kWh mỗi năm.

Thứ ba, thiết lập chiến lược chào giá tự động theo chi phí biên và cân bằng Nash trên thị trường điện giao ngay. Phòng Kế hoạch Kinh doanh cần xây dựng bảng chào giá theo từng chu kỳ 1 giờ dựa trên đường đặc tính chi phí biên trong thời gian 3 đến 6 tháng, nhằm đảm bảo công suất huy động bình quân đạt trên 190 MWh trên mỗi 350 MWh nhu cầu phụ tải hệ thống, tối đa hóa biên lợi nhuận ròng hàng tháng.

Thứ tư, Cục Điều tiết Điện lực thuộc Bộ Công Thương và Tập đoàn Điện lực Việt Nam cần mở rộng ứng dụng mô hình cân bằng Nash để tối ưu hóa công suất phát cho toàn bộ các cụm nhà máy thủy điện bậc thang trên toàn quốc. Lộ trình triển khai thực hiện trong 18 tháng, hướng tới mục tiêu giảm tổng chi phí huy động nguồn điện toàn quốc từ 251,2 triệu đồng xuống dưới 243,4 triệu đồng cho mỗi cụm phụ tải 350 MWh.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn sâu rộng cho 4 nhóm đối tượng chính:

Nhóm thứ nhất là các nhà quản lý, kỹ sư trưởng và chuyên viên vận hành tại các nhà máy thủy điện đa mục tiêu có quy mô từ 100 MW đến 300 MW như Thủy điện Đại Ninh, Đa Nhim, Trị An. Nhóm này có thể ứng dụng trực tiếp mô hình toán học và bộ mã nguồn MATLAB để lập kế hoạch điều tiết lưu lượng nước từ 0 đến 55 m3/s, loại bỏ nước xả tràn và gia tăng 8,6% doanh thu phát điện thực tế.

Nhóm thứ hai là cán bộ quản lý thị trường điện tại Cục Điều tiết Điện lực và Trung tâm Điều độ Hệ thống điện Quốc gia. Nghiên cứu cung cấp công cụ định lượng chuẩn xác để giải bài toán phân bổ tối ưu công suất phụ tải 350 MWh, thiết lập cơ chế giám sát hành vi thao túng giá và tối ưu hóa chi phí mua điện toàn hệ thống.

Nhóm thứ ba là các nhà nghiên cứu, giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Toán ứng dụng, Kỹ thuật hệ thống điện và Năng lượng tái tạo. Luận văn cung cấp khung tham chiếu hoàn chỉnh về cách chuyển đổi một bài toán kỹ thuật phức tạp thành bài toán quy hoạch toàn phương cỡ lớn với cấu trúc ma trận 2.400 dòng và 960 cột.

Nhóm thứ tư là các đơn vị tư vấn kỹ thuật điện lực và công ty đầu tư phát triển năng lượng. Nhóm này có thể khai thác phương pháp phân tích chuỗi thủy văn dài hạn 8.760 giờ để lập báo cáo khả thi kinh tế kỹ thuật và nâng cao độ tự chủ công nghệ mà không cần phụ thuộc vào phần mềm ngoại nhập đắt đỏ.

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

Tại sao Nhà máy thủy điện Đại Ninh cần áp dụng mô hình toán tối ưu thay cho biểu đồ điều phối truyền thống?

Biểu đồ điều phối truyền thống chủ yếu dựa trên kinh nghiệm lịch sử mang tính an toàn tĩnh, khiến nhà máy chỉ có 2 trong 7 năm đạt chỉ tiêu sản lượng 1,2 tỷ kWh. Mô hình toán tối ưu hóa quy hoạch toàn phương giúp điều tiết dòng chảy linh hoạt theo từng giờ, nâng sản lượng thực tế thêm 8,6% và triệt tiêu hoàn toàn nước xả thừa.

Việc ứng dụng quy hoạch toàn phương đem lại lợi ích kinh tế cụ thể như thế nào trong ngắn hạn?

Trong chu kỳ kiểm nghiệm 20 ngày với 480 giờ đầu tháng 12 năm 2014, phương pháp tối ưu đã nâng sản lượng điện từ 29,7 triệu kWh lên 32,2 triệu kWh. Với đơn giá bán điện 550 đồng/kWh, tổng doanh thu nhà máy tăng từ 16,3 tỷ đồng lên 17,7 tỷ đồng, tạo thêm 1,4 tỷ đồng giá trị thặng dư.

Lý thuyết trò chơi giải quyết vấn đề gì trong thị trường phát điện cạnh tranh?

Lý thuyết trò chơi dạng chuẩn tắc xác định điểm cân bằng Nash thông qua nguyên tắc cân bằng chi phí biên giữa các nhà máy. Tại điểm cân bằng khi điều động 350 MWh, nhà máy Đa Nhim phát 157,6 MWh và Đại Ninh phát 192,4 MWh, giúp hệ thống đạt chi phí thấp nhất là 251.224.616 đồng mà không nhà máy nào bị thiệt hại.

Khi một nhà máy thủy điện tối ưu hóa chi phí thì người tiêu dùng và xã hội được hưởng lợi ra sao?

Khi nhà máy Đại Ninh tối ưu kỹ thuật giúp giảm hệ số chi phí từ 1,8 xuống 1,6, giá chào điện giảm kéo theo tổng chi phí phát điện của hệ thống giảm từ 251,8 triệu đồng xuống 243,4 triệu đồng cho 350 MWh. Điều này trực tiếp hạ giá thành sản xuất điện và giảm áp lực tăng giá điện lên người tiêu dùng.

Mã nguồn và mô hình toán trong luận văn có thể chuyển giao cho các hồ thủy điện khác không?

Mô hình toán quy hoạch toàn phương và đoạn mã MATLAB trong luận văn có tính tổng quát hóa cao. Khi áp dụng cho các nhà máy thủy điện khác, người dùng chỉ cần thay đổi các hằng số đặc tính kỹ thuật k1, k2, dung tích hồ từ 68 đến 320 triệu m3 và dữ liệu thủy văn đầu vào tương ứng.

Kết luận

Nghiên cứu đã giải quyết trọn vẹn bài toán tích hợp giữa kỹ thuật thủy điện và kinh tế năng lượng thông qua 5 đóng góp nổi bật:

  • Xây dựng thành công mô hình toán học quy hoạch toàn phương cỡ lớn với thuật toán điểm trong tối ưu trên MATLAB để điều tiết hồ chứa thủy điện.
  • Tối ưu hóa hiệu quả vận hành ngắn hạn 20 ngày tại Thủy điện Đại Ninh, nâng doanh thu từ 16,3 tỷ lên 17,7 tỷ đồng và sản lượng từ 29,7 triệu lên 32,2 triệu kWh, tương đương mức tăng trưởng 8,6%.
  • Thiết lập quy luật vận hành tối ưu dài hạn 8.760 giờ cho từng năm thủy văn, đảm bảo tận dụng 100% lượng nước về hồ và loại bỏ hoàn toàn nước xả thừa qua đập.
  • Ứng dụng lý thuyết trò chơi và chi phí biên để xác lập trạng thái cân bằng Nash trong chào giá điện, tối ưu hóa phân bổ công suất 350 MWh giữa các nhà máy.
  • Cung cấp cơ sở khoa học giúp giảm tổng chi phí vận hành hệ thống phát điện toàn quốc từ 251,2 triệu đồng xuống 243,4 triệu đồng cho mỗi chu kỳ phụ tải tiêu chuẩn.

Trong giai đoạn 6 đến 18 tháng tới, các đơn vị quản lý cần tiến hành số hóa mô hình toán học vào quy trình điều độ thực tế để khai thác tối đa tiềm năng nguồn năng lượng tái tạo. Hãy liên hệ các chuyên gia toán ứng dụng và điều độ năng lượng để chuyển giao giải pháp, tối ưu hóa hiệu quả sản xuất kinh doanh và đóng góp tích cực vào sự phát triển bền vững của thị trường điện Việt Nam.