Tổng quan nghiên cứu

Hồ thủy điện Thác Bà là công trình thủy điện đầu tiên tại miền Bắc Việt Nam, được khởi công xây dựng vào ngày 19 tháng 8 năm 1964 với công suất thiết kế 108 MW và sản lượng điện năng bình quân hàng năm đạt 400 triệu kWh. Tọa lạc trên lưu vực sông Chảy thuộc địa bàn huyện Yên Bình, tỉnh Yên Bái, công trình giữ vai trò chiến lược đối với an ninh năng lượng quốc gia và an sinh xã hội khu vực trung du, miền núi phía Bắc. Tuy nhiên, các tác động nghiêm trọng của biến đổi khí hậu toàn cầu đang khiến chế độ thủy văn trở nên cực kỳ bất thường. Trên bình diện quốc tế, các hiện tượng thời tiết cực đoan đã từng gây ra những thiệt hại kinh tế nặng nề; điển hình như đợt hạn hán năm 2001 tại Brazil làm sụt giảm 1,5% GDP tương đương thiệt hại 10 tỷ USD do thiếu hụt công suất thủy điện, hay tình trạng cạn kiệt tại hồ Volta khiến mực nước hạ thấp hơn 1,5 mét dưới ngưỡng an toàn, đe dọa trực tiếp 60% sản lượng điện quốc gia Ghana.

Tại Việt Nam, sự biến động dòng chảy theo mùa tạo ra bài toán phi tuyến tính phức tạp trong công tác điều tiết hồ chứa. Trong cơ chế thị trường phát điện cạnh tranh, nhà máy thủy điện Thác Bà phải đồng thời giải quyết hai mục tiêu mâu thuẫn: tối đa hóa doanh thu phát điện thương mại và duy trì xả dòng chảy tối thiểu nhằm phục vụ sản xuất nông nghiệp cũng như nguồn nước sinh hoạt cho toàn bộ vùng hạ du. Đề tài luận văn tập trung nghiên cứu ứng dụng hệ suy luận mờ nơron thích nghi ANFIS (Adaptive Neuro-Fuzzy Inference System) để dự báo chính xác mực nước hồ Thác Bà theo chuỗi thời gian 12 tháng. Mục tiêu cụ thể là xây dựng mô hình toán học thông minh có khả năng xấp xỉ chính xác các mối quan hệ phi tuyến giữa lưu lượng dòng chảy đến hồ và cao trình mực nước, đạt sai số bình phương cực tiểu dưới 0,05 mét, nâng cao hiệu quả vận hành kinh tế và bảo đảm an toàn tuyệt đối cho công trình đầu mối.

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 trên nền tảng tích hợp giữa hai trường phái tính toán mềm cốt lõi: Mạng nơron nhân tạo (Artificial Neural Network - ANN) và Hệ logic mờ (Fuzzy Logic). Mạng nơron nhân tạo, khởi nguồn từ mô hình tế bào thần kinh McCulloch-Pitts và các nguyên lý học có giám sát dựa trên quy tắc Delta của Widrow (1962) cùng thuật toán Perceptron của Rosenblatt (1958), sở hữu năng lực xử lý song song mạnh mẽ và khả năng tự học dữ liệu thông qua hiệu chỉnh ma trận trọng số. Tuy nhiên, mạng nơron thuần túy mang đặc tính hộp đen, thiếu tính tường minh trong biểu diễn tri thức. Ngược lại, lý thuyết Logic mờ của Zadeh với các hệ suy diễn mờ nhiều đầu vào - một đầu ra (MISO) theo mô hình Mamdani hoặc Takagi-Sugeno cho phép mô tả tri thức chuyên gia dưới dạng các luật ngôn ngữ Nếu - Thì trực quan, nhưng lại thiếu khả năng tự động học hỏi và điều chỉnh các hàm thuộc khi điều kiện môi trường thay đổi.

Để khắc phục triệt để các nhược điểm này, mô hình mạng nơron mờ thích nghi ANFIS do Roger Jang đề xuất năm 1993 được lựa chọn làm khung lý thuyết chủ đạo. Kiến trúc ANFIS bao gồm cấu trúc 5 tầng chức năng riêng biệt:

  • Tầng 1 (Mờ hóa): Chuyển đổi các giá trị đầu vào rõ thành các mức độ thuộc mờ thông qua các hàm thuộc dạng Gauss hoặc dạng chuông Bell với tập tham số tiền đề.
  • Tầng 2 (Nhân luật): Thực hiện phép toán T-norm (thường là tích) để tính toán độ kích hoạt của từng luật mờ If-Then.
  • Tầng 3 (Chuẩn hóa): Tính tỷ lệ độ mạnh của từng luật so với tổng cường độ của tất cả các luật trong hệ thống.
  • Tầng 4 (Kết luận): Áp dụng phương trình hàm tuyến tính bậc nhất Takagi-Sugeno đối với các biến đầu vào dựa trên tập tham số hệ quả.
  • Tầng 5 (Giải mờ): Tính tổng toàn bộ tín hiệu đầu ra để cho ra giá trị dự báo mực nước dạng số rõ ràng.

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

Phương pháp nghiên cứu thực nghiệm được triển khai dựa trên chuỗi số liệu thủy văn lịch sử của hồ thủy điện Thác Bà trên sông Chảy. Bộ dữ liệu bao gồm lưu lượng dòng chảy tự nhiên bình quân tháng đổ về hồ, mực nước thượng lưu, mực nước hạ lưu, dung tích hữu ích và cột nước phát điện cực đại.

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu sử dụng tổng cộng 120 mẫu dữ liệu tháng được thu thập liên tục trong chu kỳ quan trắc 10 năm. Phương pháp chọn mẫu phân tầng theo mùa được áp dụng để chia dữ liệu thành hai giai đoạn thủy văn đặc trưng: mùa lũ (từ tháng 6 đến tháng 10) với lưu lượng biến động dữ dội và mùa cạn (từ tháng 11 đến tháng 5 năm sau) với dòng chảy duy trì ở mức thấp. Toàn bộ tập dữ liệu được phân chia theo tỷ lệ chuẩn 70% dành cho quá trình huấn luyện (training) và 30% dành cho quá trình kiểm tra, kiểm định độc lập (testing).

Lý do lựa chọn phương pháp phân tích ANFIS sử dụng thuật toán học lai (Hybrid Learning Algorithm) xuất phát từ yêu cầu khắt khe về tốc độ hội tụ và độ chính xác toán học. Thuật toán học lai kết hợp phương pháp ước lượng bình phương cực tiểu (Least Squares Estimation - LSE) ở bước lan truyền tiến để tối ưu hóa các tham số kết luận tại tầng 4, cùng với phương pháp giảm độ dốc sai số (Gradient Descent) ở bước lan truyền ngược để tinh chỉnh các tham số tiền đề tại tầng 1. Thuật toán khởi tạo ma trận hiệp phương sai với hệ số dương lớn gamma bằng 10^6, giúp hệ thống rút ngắn hơn 40% thời gian tính toán so với thuật toán lan truyền ngược truyền thống và hoàn toàn loại bỏ nguy cơ rơi vào điểm cực tiểu địa phương.

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

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

Quá trình huấn luyện và kiểm thử thực nghiệm trên 5 cấu hình mô hình ANFIS khác nhau (từ mô hình A đến mô hình E) với các tổ hợp biến đầu vào từ độ trễ thời gian t-1, t-2 đến t-12 đã mang lại những kết quả nổi bật sau:

Thứ nhất, việc kết hợp đồng thời biến lưu lượng dòng chảy tự nhiên quá khứ và cao trình mực nước hồ tại các thời điểm liền kề mang lại độ chính xác dự báo vượt trội. Trong 5 cấu hình thử nghiệm, mô hình E đạt hệ số tương quan R2 đạt mức trên 0,96, phản ánh sự trùng khớp gần như hoàn hảo giữa chuỗi số liệu dự báo và số liệu đo đạc thực tế của nhà máy thủy điện Thác Bà.

Thứ hai, sai số bình phương trung bình cực tiểu (RMSE) của mô hình ANFIS đạt mức dưới 0,042, giảm hơn 35% so với mô hình mạng nơron truyền thẳng 3 lớp thông thường và giảm 48% so với phương pháp hồi quy thống kê truyền thống. Sai số tuyệt đối trung bình về mực nước thực tế luôn duy trì ở mức biên độ hẹp dưới 0,28 mét.

Thứ ba, mô hình ANFIS thể hiện năng lực thích ứng xuất sắc trong các tháng mùa lũ cao điểm khi lưu lượng dòng chảy sông Chảy biến động đột ngột từ 300 m3/s lên trên 1200 m3/s. Sai lệch tương đối của dự báo đỉnh lũ chỉ dao động trong khoảng từ 2,1% đến 4,2%, cung cấp khoảng thời gian cảnh báo sớm an toàn từ 24 đến 48 giờ trước khi lũ về hồ chứa.

Thảo luận kết quả

Sự vượt trội của mô hình ANFIS bắt nguồn từ cơ chế tối ưu hóa tự động hình dạng và vị trí của các hàm thuộc dạng chuông Bell ở tầng 1, kết hợp với tính tuyến tính cục bộ của các luật mờ Takagi-Sugeno ở tầng 4. Điều này giúp mô hình hấp thụ trọn vẹn đặc tính phi tuyến và tính bất định cao của chuỗi thời gian thủy văn miền núi.

Trong thực tế phân tích học thuật, các kết quả này được trực quan hóa sinh động thông qua biểu đồ đường chuỗi thời gian so sánh giữa đường mực nước đo thực và các đường dự báo từ mô hình A đến mô hình E trong suốt 12 tháng của năm thủy văn. Biểu đồ chỉ ra rằng đường dự báo của mô hình ANFIS bám sát từng điểm uốn và chu kỳ dao động của mực nước thực tế. Đồng thời, bảng tổng hợp sai số LSE và hệ số xác định giữa các mô hình đã chứng minh rõ ràng việc tăng số lượng luật mờ hợp lý giúp giảm thiểu tổn thất dự báo một cách tối đa.

Khi so sánh với các nghiên cứu thủy văn công bố trước đây dựa trên phương pháp quy hoạch động ngẫu nhiên, giải pháp ANFIS trong luận văn cho phép nhà máy nâng cao hiệu quả kinh tế phát điện thêm khoảng 12% nhờ khả năng tích trữ nước tối ưu vào cuối mùa lũ, đồng thời cắt giảm tới 25% rủi ro thiếu hụt nguồn nước tưới tiêu cho hạ du tỉnh Yên Bái trong các đợt hạn hán kéo dài.

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

Dựa trên kết quả nghiên cứu toán học và thực nghiệm tại công trình thủy điện Thác Bà, các giải pháp và khuyến nghị hành động cụ thể được đề xuất như sau:

  1. Tích hợp phần mềm dự báo ANFIS vào hệ thống SCADA điều hành: Ban Giám đốc Nhà máy Thủy điện Thác Bà chủ trì triển khai cài đặt mô hình toán vào máy chủ điều độ trung tâm trong thời hạn 6 tháng, đặt mục tiêu rút ngắn thời gian xử lý dữ liệu dự báo thời gian thực xuống dưới 3 giây mỗi chu kỳ tính toán.
  2. Mở rộng hệ thống 15 trạm đo mưa và quan trắc dòng chảy tự động: Công ty Cổ phần Thủy điện Thác Bà phối hợp cùng Sở Tài nguyên và Môi trường tỉnh Yên Bái đầu tư lắp đặt bổ sung dọc lưu vực sông Chảy trong vòng 12 tháng, nhằm nâng cao độ tin cậy của dữ liệu đầu vào đạt tỷ lệ trên 98% và khống chế sai số quan trắc dưới 2%.
  3. Ban hành quy chế phối hợp xả nước đa mục tiêu giữa phát điện và thủy lợi: Ủy ban Nhân dân tỉnh Yên Bái và Tập đoàn Điện lực Việt Nam xây dựng khung quy chế vận hành hồ chứa tối ưu vào quý 2 năm 2027, bảo đảm cam kết duy trì lưu lượng xả hạ du tối thiểu 150 m3/s trong các tháng đổ ải vụ đông xuân, đáp ứng 100% nhu cầu nước canh tác nông nghiệp.
  4. Đào tạo nâng cao năng lực ứng dụng công nghệ tính toán mềm cho đội ngũ kỹ thuật: Phòng Kỹ thuật công nghệ của nhà máy tổ chức 4 khóa tập huấn chuyên sâu trong vòng 3 tháng cho toàn bộ 25 kỹ sư điều độ, hướng tới làm chủ 100% kỹ năng tái huấn luyện tham số mạng nơron mờ khi xuất hiện các biến số thủy văn mới.

Đố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ị ứng dụng thiết thực cho 4 nhóm đối tượng trọng tâm:

  1. Kỹ sư vận hành và điều độ tại các nhà máy thủy điện: Nắm vững phương pháp luận xây dựng kế hoạch tích nước và phát điện tối ưu dài hạn lẫn ngắn hạn, hỗ trợ công tác chào giá cạnh tranh trên thị trường điện lực nhằm gia tăng doanh thu hàng năm từ 8% đến 15%.
  2. Cơ quan quản lý tài nguyên nước và phòng chống thiên tai địa phương: Khai thác công cụ mô phỏng toán học để thiết lập các kịch bản điều tiết liên hồ chứa, nâng cao năng lực ứng phó ngập lụt vùng hạ du lưu vực sông Chảy với độ tin cậy kỹ thuật đạt trên 95%.
  3. Học viên cao học và nghiên cứu sinh ngành Công nghệ thông tin và Trí tuệ nhân tạo: Tiếp cận một công trình nghiên cứu hoàn chỉnh về kiến trúc mạng nơron mờ 5 tầng, thuật toán học lai LSE kết hợp lan truyền ngược, làm tài liệu tham khảo cho các đề tài xử lý chuỗi thời gian phi tuyến.
  4. Giảng viên và nhà nghiên cứu tại các trường đại học khối kỹ thuật: Sử dụng bộ khung lý thuyết logic mờ Mamdani/Sugeno và thuật toán thích nghi ANFIS làm bài giảng chuyên đề hoặc dữ liệu kiểm chứng cho các công trình nghiên cứu khoa học cấp bộ và cấp quốc gia.

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

  1. Mô hình ANFIS có ưu điểm gì vượt trội so với mạng nơron truyền thống trong dự báo mực nước? Mô hình ANFIS kết hợp cấu trúc 5 tầng của nơron nhân tạo với hệ luật mờ Sugeno tường minh, giúp loại bỏ hoàn toàn nhược điểm hộp đen khó giải thích của mạng nơron truyền thống. Trong thực nghiệm tại hồ Thác Bà, ANFIS đạt sai số dự báo thấp hơn 35% so với mạng truyền thẳng nhờ khả năng tự động tối ưu hóa hàm thuộc dạng chuông.

  2. Thuật toán học lai vận hành như thế nào để tối ưu hóa quá trình huấn luyện? Thuật toán học lai phân tách quá trình thành hai chiều: chiều tiến sử dụng phương pháp bình phương cực tiểu LSE để cố định tham số tiền đề và tìm tham số kết luận ở tầng 4; chiều lùi sử dụng Gradient Descent cập nhật tham số tiền đề ở tầng 1. Quy trình này giúp tốc độ hội tụ nhanh hơn 40% và tránh bẫy cực tiểu địa phương.

  3. Biến đầu vào nào có ảnh hưởng quyết định nhất đến độ chính xác của mô hình dự báo? Biến đầu vào quyết định nhất là lưu lượng dòng chảy bình quân tháng kết hợp cùng cao trình mực nước hồ ở các thời điểm quá khứ t-1 và t-2. Kết quả phân tích 5 mô hình thực nghiệm cho thấy việc tích hợp đầy đủ chuỗi dữ liệu 12 tháng giúp nâng hệ số xác định R2 từ 0,82 lên mức kỷ lục 0,96.

  4. Kết quả dự báo của ANFIS hỗ trợ như thế nào cho bài toán chào giá trên thị trường điện? Dự báo chính xác mực nước giúp nhà máy thủy điện công suất 108 MW chủ động phân bổ lượng nước phát điện vào các khung giờ cao điểm có giá mua điện cao nhất. Nhờ đó, doanh nghiệp có thể tối đa hóa doanh thu phát điện thêm 10% đến 15% mỗi năm mà vẫn bảo đảm các chỉ tiêu an toàn kỹ thuật.

  5. Mô hình ANFIS xây dựng cho hồ Thác Bà có thể mở rộng áp dụng cho các hồ chứa khác không? Mô hình hoàn toàn có thể nhân rộng cho hơn 80 hồ thủy điện lớn nhỏ trên toàn lãnh thổ Việt Nam như Hòa Bình, Sơn La, Tuyên Quang hay Trị An. Khi áp dụng cho lưu vực mới, chỉ cần thay đổi tập dữ liệu huấn luyện và hiệu chỉnh lại miền giá trị của các hàm thuộc mờ cho phù hợp với đặc thù địa phương.

Kết luận

  • Luận văn đã xây dựng thành công mô hình suy luận mờ nơron thích nghi ANFIS 5 tầng hoàn chỉnh, giải quyết xuất sắc bài toán dự báo mực nước chuỗi thời gian phi tuyến tính tại hồ thủy điện Thác Bà.
  • Chứng minh tính hiệu quả vượt bậc của thuật toán học lai LSE kết hợp Gradient Descent, giúp giảm hơn 35% sai số bình phương trung bình và nâng hệ số tương quan thực nghiệm R2 lên trên 0,96.
  • Cung cấp cơ sở khoa học vững chắc giúp nhà máy thủy điện Thác Bà với công suất 108 MW tối ưu hóa kế hoạch phát điện thương mại, đồng thời bảo đảm 100% nhu cầu cung cấp nước an sinh cho hạ du.
  • Đề xuất lộ trình công nghệ 12 tháng rõ ràng nhằm tích hợp giải pháp toán học thông minh vào hệ thống điều độ và giám sát vận hành thực tế tại doanh nghiệp.
  • Mở ra tiềm năng ứng dụng rộng lớn của các công nghệ tính toán mềm trong quản lý tài nguyên nước, thủy văn công trình và thích ứng với biến đổi khí hậu trên phạm vi cả nước.

Công trình luận văn là một đóng góp học thuật nghiêm túc và có giá trị thực tiễn to lớn đối với ngành công nghệ thông tin và kỹ thuật năng lượng. Các cơ quan quản lý, doanh nghiệp thủy điện và cộng đồng nghiên cứu khoa học quan tâm có thể khai thác, ứng dụng giải pháp ANFIS này nhằm nâng cao chất lượng quản trị và phát triển nguồn năng lượng xanh bền vững.