Tổng quan nghiên cứu
Lưu vực sông Ba là một trong những mạng lưới sông ngòi lớn nhất khu vực Nam Trung Bộ và Tây Nguyên với tổng diện tích tự nhiên 13.900 km2, bao phủ địa bàn 4 tỉnh gồm Kon Tum, Gia Lai, Đắk Lắk và Phú Yên. Với chiều dài dòng chính 374 km cùng mật độ xây dựng công trình dày đặc đạt 53 km2/hồ (tổng số 254 hồ chứa thủy lợi và thủy điện), nguồn nước mặt tại đây đang chịu sức ép khai thác vô cùng khốc liệt. Tình trạng chuyển nước liên lưu vực, tiêu biểu như công trình hồ An Khê chuyển dòng sang sông Kôn với công suất 160 MW, cùng sự gia tăng nhu cầu nước sinh hoạt và nông nghiệp đã làm nảy sinh các xung đột gay gắt giữa thượng lưu và hạ lưu, giữa mục tiêu phát điện và cấp nước tưới trong mùa kiệt.
Vấn đề biến động dòng chảy theo mùa diễn ra cực đoan: mùa lũ chỉ kéo dài 3 đến 4 tháng nhưng chiếm 60% đến 65% tổng lượng nước cả năm, từng ghi nhận lưu lượng lũ lịch sử tại trạm Củng Sơn lên tới 20.700 m3/s. Trong khi đó, mùa khô kéo dài 8 đến 9 tháng khiến lưu lượng kiệt tại hạ du suy giảm nghiêm trọng. Mục tiêu cốt lõi của nghiên cứu là xây dựng bộ công cụ khoa học tích hợp đa mô hình nhằm dự báo định lượng tài nguyên nước mặt theo các hạn dự báo: ngắn hạn 5 ngày, trung hạn 10 ngày, theo tháng và theo mùa. Nghiên cứu tập trung giải quyết bài toán cân bằng nước tại các nút thủy văn xung yếu gồm An Khê, Ayun Hạ, Củng Sơn và Sông Hinh. Kết quả nghiên cứu có ý nghĩa quyết định trong việc tối ưu hóa quy trình vận hành của 164 hồ chứa đang hoạt động, đảm bảo an ninh nguồn nước và hỗ trợ thực thi hiệu quả Quyết định 878/QĐ-TTg cùng Quyết định 1879/QĐ-TTg của Thủ tướng Chính phủ về vận hành liên hồ chứa trên lưu vực sông Ba.
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 sự kết hợp chặt chẽ giữa 3 nền tảng lý thuyết và mô hình toán thủy văn hiện đại:
- Lý thuyết bảo toàn vật chất và cân bằng nước lưu vực: Dựa trên chu trình thủy văn tổng thể, mô phỏng quá trình chuyển hóa từ lượng mưa thành dòng chảy mặt, dòng chảy sát mặt và dòng chảy ngầm thông qua cấu trúc 4 bể chứa tương tác liên tục của mô hình MIKE-NAM (bể chứa tầng mặt, bể chứa sát mặt, bể chứa tầng rễ cây và bể chứa ngầm).
- Lý thuyết mạng nơ-ron nhân tạo (Artificial Neural Network - ANN): Ứng dụng cấu trúc mạng truyền thẳng nhiều lớp gồm lớp đầu vào (input), lớp ẩn (hidden) và lớp đầu ra (output) để giải quyết các mối quan hệ phi tuyến phức tạp giữa các chỉ số hoàn lưu khí quyển, bức xạ, nhiệt độ, bốc hơi và lưu lượng dòng chảy đến các trạm thủy văn.
- Hệ phương trình thủy động lực học một chiều Saint-Venant: Mô hình MIKE 11 giải hệ phương trình vi phân đạo hàm riêng gồm phương trình liên tục và phương trình động lượng bằng phương pháp sai phân hữu hạn 6 điểm ẩn để mô phỏng diễn toán dòng chảy và mực nước hạ du.
Các khái niệm chuyên ngành then chốt được lượng hóa bao gồm: mô đun dòng chảy năm trung bình đạt 22,3 đến 24,7 l/s.km2, mô đun dòng chảy kiệt cực tiểu đạt 0,22 l/s.km2, lưu lượng dòng chảy trung bình năm tại An Khê (33,4 m3/s) và Củng Sơn (277 m3/s), cùng nguyên lý phân bổ dòng chảy tĩnh qua các nút mạng trên mô hình MIKE BASIN.
Phương pháp nghiên cứu
Nguồn số liệu nghiên cứu được thu thập đồng bộ từ chuỗi quan trắc khí tượng thủy văn giai đoạn 40 năm (1977 - 2016) trên toàn lưu vực, bao gồm 20 trạm khí tượng đo mưa, nhiệt độ, độ ẩm, bốc hơi và 4 trạm thủy văn trọng điểm (An Khê, Củng Sơn, Sông Hinh, Krông Hnăng).
Về cỡ mẫu và phương pháp chọn mẫu, toàn bộ chuỗi số liệu 40 năm được phân chia khoa học theo phương pháp mẫu đại diện chuỗi thời gian: giai đoạn 1980 - 2000 (21 năm) phục vụ huấn luyện và hiệu chỉnh thông số mô hình; giai đoạn 2001 - 2015 (15 năm) phục vụ kiểm định tính độc lập của bộ công cụ trong điều kiện lưu vực đã xuất hiện nhiều công trình hồ chứa lớn; giai đoạn 2016 - 2018 sử dụng để thử nghiệm dự báo thực tế cho các mùa cạn và mùa lũ.
Lý do lựa chọn phương pháp tích hợp đa mô hình xuất phát từ thực tiễn mạng lưới quan trắc lưu vực bị cắt giảm (hiện chỉ còn 2 trạm An Khê và Củng Sơn đo lưu lượng) và chế độ dòng chảy bị điều tiết phức tạp bởi 254 công trình hồ đập. Mô hình ANN được chọn cho bài toán dự báo dài hạn tháng và mùa nhờ khả năng tự học sâu sắc quy luật khí hậu; mô hình MIKE-NAM giải quyết bài toán dự báo dòng chảy đến hồ ngắn hạn 5 đến 10 ngày từ dự báo mưa số trị; mô hình MIKE BASIN kết nối toàn diện mạng lưới phân bổ nước; mô hình MIKE 11 kiểm soát lan truyền thủy lực vùng hạ lưu.
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 dữ liệu và chạy thử nghiệm bộ công cụ dự báo đã đem lại các phát hiện khoa học quan trọng:
- Sự phân hóa sâu sắc của dòng chảy mặt: Dòng chảy tập trung cao độ vào mùa lũ (tháng 9 đến tháng 12 hoặc tháng 10 đến tháng 12) chiếm 60% đến 65% tổng lượng cả năm, trong đó tháng 11 có lượng dòng chảy lớn nhất, chiếm 24% đến 26% tổng lượng năm. Ngược lại, mùa kiệt kéo dài từ 8 đến 9 tháng có lưu lượng suy kiệt nghiêm trọng, với mô đun kiệt thấp kỷ lục đo được tại An Khê là 0,22 l/s.km2 và tại Củng Sơn là 0,62 l/s.km2 vào các năm hạn hán khốc liệt.
- Độ chính xác vượt trội của mô hình ANN trong dự báo dài hạn: Kết quả dự báo thử nghiệm tổng lượng dòng chảy tháng và mùa cạn, mùa lũ tại trạm An Khê và Củng Sơn giai đoạn 2017 - 2018 cho thấy mô hình ANN bám sát xu thế dòng chảy thực đo. Sai số tương đối của tổng lượng dòng chảy tháng phần lớn nằm trong giới hạn cho phép từ 10% đến 18%, đáp ứng tốt yêu cầu lập kế hoạch cấp nước trước 1 đến 6 tháng.
- Hiệu quả mô phỏng thủy văn của MIKE-NAM và MIKE BASIN: Trong giai đoạn hiệu chỉnh (1980 - 2000) và kiểm định (2001 - 2015), mô hình MIKE-NAM đạt hệ số tương quan Nash-Sutcliffe từ 0,75 đến 0,86 tại các tiểu lưu vực thượng nguồn. MIKE BASIN đã tái hiện thành công dòng chảy kiệt và quá trình điều tiết nước của 164 hồ chứa đang vận hành, khôi phục dòng chảy tại các sông nhánh thiếu số liệu như Krông Hnăng và Sông Hinh.
Thảo luận kết quả
Sự sai khác về chế độ dòng chảy trên lưu vực bắt nguồn từ yếu tố địa hình chắn gió của dải Trường Sơn. Sườn Tây Trường Sơn chịu chi phối của gió mùa Tây Nam có mùa mưa từ tháng 5 đến tháng 10, trong khi sườn Đông Trường Sơn chịu tác động của bão và gió mùa Đông Bắc có mùa mưa muộn từ tháng 9 đến tháng 12. Sự lệch pha này kết hợp với các hồ chứa lớn như Ka Nak (dung tích 313,7 triệu m3) và Sông Ba Hạ đã làm thay đổi hoàn toàn quy luật thủy văn tự nhiên.
Về mặt trực quan hóa dữ liệu, các kết quả được thể hiện rõ nét qua biểu đồ quá trình dòng chảy (hydrograph) so sánh giữa đường thực đo và đường tính toán qua các pha kiểm định, kết hợp cùng bảng ma trận sai số tương đối (Nash-Sutcliffe, sai số tổng lượng W%). Đồng thời, sơ đồ phân bổ dòng chảy được mô phỏng sinh động qua giao diện GIS của MIKE BASIN, giúp định vị chính xác mức độ thiếu hụt nguồn nước tại từng nút mạng hạ du.
So với các nghiên cứu trước đây chỉ sử dụng phương pháp thống kê tương quan đơn biến hoặc mô hình chuỗi thời gian ARIMA (vốn chỉ dựa vào tự tương quan dòng chảy quá khứ nên độ chính xác hạn chế), việc kết hợp mạng nơ-ron ANN với bộ mô hình MIKE đã tạo bước đột phá. Bộ công cụ không chỉ khắc phục nhược điểm thiếu trạm đo lưu lượng trên sông nhánh mà còn mô phỏng chân thực bức tranh vận hành điều tiết nước liên hồ, cung cấp số liệu dự báo tin cậy cho cả mùa lũ lẫn mùa cạn.
Đề xuất và khuyến nghị
Dựa trên kết quả nghiên cứu, 4 nhóm giải pháp chiến lược kèm mục tiêu định lượng và lộ trình triển khai cụ thể được đề xuất:
- Tự động hóa và vận hành bộ công cụ dự báo trực tuyến: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia phối hợp với các cơ quan chuyên môn tích hợp dữ liệu mưa dự báo số trị (NWP) vào hệ thống mô hình ANN và MIKE-NAM. Mục tiêu nâng độ chính xác dự báo dòng chảy đến các hồ chứa hạn 5 đến 10 ngày lên trên 85%, hoàn thành trong giai đoạn 2024 - 2026.
- Cập nhật và tối ưu hóa Quy trình vận hành liên hồ chứa: Bộ Tài nguyên và Môi trường cùng Ủy ban nhân dân 4 tỉnh trên lưu vực tiến hành rà soát quy chế vận hành theo Quyết định 878/QĐ-TTg. Bắt buộc các hồ thủy điện lớn duy trì lưu lượng xả tối thiểu bảo vệ môi trường hạ du tại Củng Sơn không dưới 28,7 m3/s trong các tháng kiệt cao điểm (tháng 3 và tháng 4), thực hiện đánh giá định kỳ hàng năm từ năm 2025.
- Hiện đại hóa mạng lưới quan trắc khí tượng thủy văn chuyên dùng: Cục Quản lý Tài nguyên nước và các chủ hồ đầu tư lắp đặt bổ sung 15 trạm đo mưa tự động và khôi phục 5 trạm thủy văn đo lưu lượng tự động trên các sông nhánh chính (Ia Yun, Krông Hnăng, Sông Hinh). Giải pháp nhằm giảm sai số ước tính diện mưa lưu vực từ 20% xuống dưới 10% trước năm 2027.
- Thiết lập nền tảng chia sẻ dữ liệu số dùng chung: Sở Tài nguyên và Môi trường 4 tỉnh phối hợp cùng Tập đoàn Điện lực Việt Nam (EVN) xây dựng hệ thống cơ sở dữ liệu thời gian thực, tự động cập nhật mực nước và dung tích của 254 hồ chứa với tần suất 1 giờ/lần, hoàn thiện đưa vào sử dụng trong năm 2025.
Đối tượng nên tham khảo luận văn
Nội dung và sản phẩm của luận văn mang giá trị thực tiễn cao cho 4 nhóm đối tượng chính:
- Cơ quan quản lý nhà nước về tài nguyên nước: Bộ Tài nguyên và Môi trường, Cục Quản lý Tài nguyên nước, Sở Tài nguyên và Môi trường các tỉnh Gia Lai, Phú Yên, Đắk Lắk, Kon Tum. Ứng dụng kết quả tính toán để xây dựng quy hoạch phân bổ nguồn nước mặt, giám sát hạn hán, kiểm soát xâm nhập mặn và ban hành lệnh điều hành xả lũ liên hồ chứa.
- Đơn vị quản lý và vận hành công trình thủy điện, thủy lợi: Công ty Thủy điện Sông Ba Hạ, Công ty Cổ phần Thủy điện Sông Hinh, Ban Quản lý Thủy điện An Khê - Ka Nak, Công ty TNHH MTV Khai thác Công trình Thủy lợi Ayun Hạ. Áp dụng mô hình MIKE-NAM và ANN để dự báo lượng nước về hồ trước 5 đến 10 ngày, từ đó lập kế hoạch tích trữ nước và phát điện tối ưu doanh thu.
- Cộng đồng nghiên cứu khoa học và cơ sở đào tạo: Giảng viên, nhà nghiên cứu, học viên cao học chuyên ngành Thủy văn học, Kỹ thuật tài nguyên nước, Biến đổi khí hậu tại các viện nghiên cứu và trường đại học. Sử dụng làm tài liệu tham khảo chuyên sâu về phương pháp luận tích hợp trí tuệ nhân tạo (ANN) với mô hình thủy lực - thủy văn (MIKE).
- Các tổ chức phát triển nông nghiệp và phòng chống thiên tai: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các địa phương, cơ quan quy hoạch nông nghiệp vùng Nam Trung Bộ. Khai thác dữ liệu dự báo kiệt và lũ để chủ động xây dựng lịch thời vụ, điều chỉnh cơ cấu cây trồng và lập bản đồ cảnh báo ngập lụt vùng hạ lưu.
Câu hỏi thường gặp
Tại sao việc dự báo tài nguyên nước trên lưu vực sông Ba lại phức tạp hơn các lưu vực khác?
Lưu vực rộng 13.900 km2 có địa hình chia cắt mạnh giữa sườn Đông và Tây Trường Sơn tạo ra chế độ mưa lệch pha hoàn toàn. Thêm vào đó, mật độ công trình rất dày đặc với 254 hồ chứa và các tuyến chuyển nước liên lưu vực như An Khê làm biến dạng dòng chảy tự nhiên, gây khó khăn lớn cho việc tính toán thủy văn truyền thống.
Mô hình mạng nơ-ron nhân tạo ANN thể hiện ưu điểm gì nổi bật trong nghiên cứu này?
ANN có khả năng xử lý phi tuyến tính cực mạnh giữa chuỗi số liệu khí hậu nhiều năm và lưu lượng dòng chảy. Mô hình khắc phục hoàn toàn nhược điểm của các phương pháp thống kê cũ như ARIMA, cho phép dự báo dài hạn tổng lượng nước theo tháng và theo mùa với sai số thấp từ 10% đến 18%.
Mô hình MIKE-NAM giải quyết bài toán dự báo ngắn hạn như thế nào?
MIKE-NAM mô phỏng liên tục dòng chảy qua 4 bể chứa vật lý đại diện cho từng tiểu lưu vực. Khi tiếp nhận dữ liệu mưa dự báo số trị thời gian thực, mô hình tính toán chuẩn xác lượng dòng chảy gia nhập lòng sông và lượng nước đổ về hồ chứa trong thời hạn từ 5 đến 10 ngày tới.
Kết quả nghiên cứu hỗ trợ trực tiếp như thế nào cho Quy trình vận hành liên hồ chứa 878/QĐ-TTg?
Nghiên cứu cung cấp công cụ định lượng dòng chảy dự báo đến từng công trình trọng điểm như hồ Ka Nak (313,7 triệu m3), Sông Hinh, Sông Ba Hạ. Nhờ đó, các nhà quản lý có căn cứ khoa học để điều phối dung tích cắt lũ mùa mưa và duy trì dòng chảy môi trường tối thiểu mùa khô.
Làm cách nào luận văn khắc phục được tình trạng thiếu hụt trạm đo thủy văn trên các sông nhánh?
Nghiên cứu sử dụng mô hình MIKE BASIN kết hợp MIKE-NAM để mô phỏng dòng chảy mưa và khôi phục chuỗi dữ liệu 40 năm (1977 - 2016) trên các nhánh sông Ia Yun, Krông Hnăng, Sông Hinh, giúp đồng bộ hóa dữ liệu toàn lưu vực một cách chính xác.
Kết luận
- Xây dựng thành công bộ công cụ tích hợp toàn diện gồm mạng nơ-ron nhân tạo ANN, mô hình thủy văn MIKE-NAM, mô hình cân bằng nước MIKE BASIN và mô hình thủy lực MIKE 11 cho lưu vực sông Ba diện tích 13.900 km2.
- Giải quyết triệt để bài toán dự báo đa thời hạn: từ dự báo ngắn hạn 5 đến 10 ngày phục vụ điều tiết hồ chứa đến dự báo dài hạn theo tháng và mùa với độ chính xác cao.
- Đánh giá định lượng sâu sắc bức tranh phân bố dòng chảy lũ (chiếm 60% đến 65% lượng nước năm) và dòng chảy kiệt cực tiểu (0,22 l/s.km2) dưới tác động của 254 hồ chứa.
- Cung cấp luận cứ khoa học vững chắc hỗ trợ đắc lực cho công tác thực thi Quy trình vận hành liên hồ chứa theo Quyết định 878/QĐ-TTg và phòng chống thiên tai tại 4 tỉnh Nam Trung Bộ - Tây Nguyên.
- Mở ra hướng đi tiên phong trong việc ứng dụng công nghệ trí tuệ nhân tạo kết hợp mô hình số trị thủy văn tại Việt Nam.
Đóng góp lớn nhất của luận văn là đã hoàn thiện một khung công cụ dự báo tài nguyên nước mặt có độ tin cậy cao và tính ứng dụng thực tiễn vượt trội. Để phát huy tối đa giá trị nghiên cứu, các cơ quan quản lý và đơn vị vận hành hồ chứa cần sớm triển khai số hóa, tích hợp bộ công cụ vào hệ thống tác nghiệp thời gian thực trong giai đoạn 2024 - 2026. Hãy liên hệ các cơ quan chuyên môn về tài nguyên nước hoặc tham khảo toàn văn luận văn tại Thư viện Đại học Thủy lợi để tiếp nhận chuyển giao công nghệ và tối ưu hóa hiệu quả quản trị nguồn nước ngay hôm nay.