Tổng quan nghiên cứu

Lưu vực sông Đà sở hữu diện tích tự nhiên rộng lớn lên tới 52.900 km², trong đó phần diện tích nằm trên lãnh thổ Việt Nam chiếm 26.800 km², tương đương 50,7% toàn bộ lưu vực. Địa hình khu vực mang đặc trưng vùng núi cao với độ cao trung bình 965 m và độ dốc bình quân bề mặt lên tới 36,8%, tạo nên mạng lưới thủy văn có dòng chảy lũ tập trung nhanh, cường suất lũ lớn và diễn biến thủy lực hết sức phức tạp. Sông Đà đóng vai trò là chi lưu lớn nhất của hệ thống sông Hồng, cung cấp phần lớn lượng nước đổ về hạ lưu đồng bằng Bắc Bộ. Do đó, việc kiểm soát chế độ dòng chảy và vận hành an toàn công trình thủy điện Hòa Bình giữ vị trí then chốt trong chiến lược phòng chống lũ lụt cho thủ đô Hà Nội và bảo đảm an ninh năng lượng quốc gia.

Sự biến đổi khí hậu toàn cầu kết hợp với những thay đổi mạnh mẽ trong hiện trạng thảm phủ bề mặt lưu vực đã khiến quy luật lũ sông Đà ngày càng biến động bất thường. Công tác vận hành đa mục tiêu hồ chứa Hòa Bình đòi hỏi phải có thông tin dự báo lưu lượng dòng chảy đến hồ trước 48 giờ với độ tin cậy cao. Tuy nhiên, các phương pháp dự báo truyền thống dạng mô hình toán kinh nghiệm hoặc mô hình hộp đen bộc lộ nhiều hạn chế khi áp dụng cho lưu vực có độ dốc lớn và thiếu hụt chuỗi dữ liệu thực đo đồng bộ.

Nghiên cứu tập trung giải quyết bài toán trên bằng cách xây dựng và ứng dụng công cụ tính toán tích hợp giữa mô hình thủy văn tham số phân bố MARINE và mô hình thủy lực một chiều IMech1D, kết hợp thuật toán lọc Kalman rời rạc. Phạm vi nghiên cứu bao phủ 26.800 km² lưu vực sông Đà thuộc Việt Nam, tích hợp chuỗi số liệu từ 36 trạm đo mưa khí tượng thủy văn và 121 mặt cắt ngang lòng dẫn sông chính. Hệ thống mang lại ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp công cụ dự báo thủy văn thủy lực hiện đại phục vụ đắc lực cho công tác điều hành hồ chứa và phòng chống thiên tai.

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 ba trụ cột lý thuyết toán - cơ học thủy khí và thống kê ngẫu nhiên hiện đại:

Mô hình thủy văn tham số phân bố MARINE: Khung lý thuyết của MARINE dựa trên hệ phương trình bảo toàn khối lượng và phương trình sóng động học hai chiều của Saint-Venant để mô phỏng sự hình thành dòng chảy tràn trên sườn dốc lưu vực từ mưa rào. Quá trình tổn thất do thấm được tính toán chi tiết bằng mô hình thấm phi tuyến Green-Ampt dựa trên định luật Darcy, phản ánh trực tiếp các thông số áp suất mao dẫn màng ẩm, độ dẫn thủy lực bão hòa và độ ẩm ban đầu của các lớp đất.

Mô hình thủy lực một chiều IMech1D: Dòng chảy không ổn định trong mạng lưới sông hở và quá trình ngập lụt tại các ô chứa (ô ruộng) được mô phỏng bằng hệ phương trình vi phân phi tuyến Saint-Venant một chiều đầy đủ, gồm phương trình liên tục và phương trình động lượng có xét đến lực ma sát cản đáy Manning và trao đổi lưu lượng phụ. Quá trình trao đổi nước qua đê, đập tràn được mô hình hóa theo các quy luật dòng chảy qua đập tràn tự do và đập tràn ngập.

Lý thuyết lọc Kalman rời rạc: Thuật toán lọc Kalman tuyến tính được ứng dụng như một kỹ thuật tối ưu hóa ngẫu nhiên trong không gian trạng thái nhằm kết hợp giữa giá trị dự báo tiên nghiệm từ mô hình với chuỗi giá trị đo đạc thực tế tại các trạm kiểm tra, từ đó đưa ra ước lượng hậu nghiệm chính xác và liên tục hiệu chỉnh sai số dự báo.

Các khái niệm chính bao gồm: dòng chảy mặt sườn dốc, quá trình truyền lũ trong lòng dẫn, lưu lượng gia nhập khu giữa, và ma trận hiệp phương sai sai số.

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

Nghiên cứu khai thác nguồn dữ liệu không gian và thủy văn đồ sộ: Bản đồ số độ cao DEM kích thước ô lưới 50 m x 50 m trích xuất từ bản đồ địa hình tỷ lệ 1:50.000; bản đồ phân loại đất tỷ lệ 1:1.000.000 phân chia thành 29 loại đất và nhóm gọn lại thành 4 nhóm đất cơ bản theo chuẩn quốc tế; bản đồ hiện trạng sử dụng đất tỷ lệ 1:100.000 với 13 loại thảm phủ; chuỗi số liệu mưa quan trắc đồng bộ tại 36 trạm khí tượng thủy văn; cùng hệ thống dữ liệu hình học 121 mặt cắt ngang đo đạc lòng dẫn sông Đà từ Mường Tè về đến đập Hòa Bình.

Phương pháp chọn mẫu không gian được triển khai toàn diện trên toàn bộ diện tích 26.800 km² bằng kỹ thuật GIS, phân chia lưu vực thành 10 tiểu lưu vực chức năng (6 tiểu lưu vực chứa sông chính và 4 tiểu lưu vực chứa sông nhánh trọng điểm). Vùng ảnh hưởng lượng mưa của 36 trạm quan trắc được thành lập theo phương pháp đa giác Thiessen.

Lý do lựa chọn phương pháp phân tích: Mô hình phân bố MARINE cho phép tái hiện chi tiết tính không đồng nhất về mặt không gian của địa hình và thổ nhưỡng; mô hình IMech1D sử dụng lược đồ sai phân ẩn 4 điểm Preissmann kết hợp thuật toán tối ưu hóa khử nút theo bậc liên kết giúp giải hệ phương trình đại số tuyến tính quy mô lớn trên mạng lưới sông phức tạp với độ ổn định cao và thời gian xử lý nhanh chóng; thuật toán lọc Kalman giúp tự động hiệu chỉnh độ lệch tức thời mà không cần can thiệp thủ công vào các hệ số vật lý.

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

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

Quá trình số hóa và tích hợp không gian lưu vực sông Đà trên nền DEM 50 m x 50 m đã hỗ trợ xác định chính xác mạng lưới sông suối và hướng dòng chảy, phân chia thành công 10 tiểu lưu vực bộ phận. Sự phân bố không gian lượng mưa qua 36 trạm quan trắc cho thấy tâm mưa tập trung lớn nhất tại vùng sườn tây dãy Hoàng Liên Sơn (lưu vực Nậm Na và Nậm Mu) với lượng mưa trung bình năm đạt từ 2500 mm đến 2997 mm, đóng góp lượng dòng chảy lũ khu giữa đột biến với môđun đỉnh lũ lên đến 2000 - 3000 l/s/km².

Phát hiện đột phá về mặt tính toán số trị trong IMech1D là thuật toán tự động khử đuổi các nút có bậc liên kết thấp. Đối với mạng sông Hồng - Thái Bình mở rộng gồm 2093 đoạn sông, 1103 nút và 278 ô ruộng, phương án khử đuổi các nút hợp lưu bậc nhỏ hơn hoặc bằng 10 đã làm giảm thời gian tính toán cho một cơn lũ lớn từ 2 giờ 50 phút xuống còn 32 phút 50 giây, tương đương mức giảm thời gian tính toán lên tới 80,7%, tạo điều kiện tiên quyết cho việc ứng dụng mô hình trong tác nghiệp dự báo thời gian thực.

Kết quả kiểm định mô hình IMech1D qua 12 bài toán kiểm định mẫu tiêu chuẩn của các phòng thí nghiệm thủy lực châu Âu đều cho kết quả trùng khớp với nghiệm giải tích. Việc kết nối mô hình MARINE (tính toán lưu lượng gia nhập sườn dốc khu giữa) và IMech1D (tính toán lan truyền sóng lũ trong lòng dẫn 121 mặt cắt) đã tái hiện chính xác diễn biến trận lũ lịch sử năm 2006 và mô phỏng tác nghiệp mùa lũ năm 2009. Thời gian truyền lũ thực tế từ Lai Châu về Hòa Bình được mô phỏng chuẩn xác trong khoảng từ 6 đến 12 giờ khi có hồ chứa, rút ngắn đáng kể so với 12 đến 24 giờ trong điều kiện tự nhiên.

Ứng dụng kỹ thuật lọc Kalman rời rạc đã tối ưu hóa quá trình cập nhật trạng thái dòng chảy, giúp giảm độ lệch sai số giữa lưu lượng mô phỏng và lưu lượng thực đo tại các trạm kiểm tra thủy văn từ khoảng 15% đến 25% trong các giai đoạn lũ lên nhanh có cường suất đạt tới 77,4 cm/h.

Thảo luận kết quả

Khác với các mô hình thủy lực thương mại đóng gói như MIKE 11 hay HEC-RAS vốn khó can thiệp mã nguồn và hạn chế bản quyền, IMech1D cho phép tùy biến toàn diện cấu trúc liên kết mạng sông vòng và tích hợp linh hoạt 16 đường cong điều chỉnh hệ số nhám Manning theo mực nước thực đo tại các trạm kiểm tra chuẩn. Việc kết nối trực tiếp với mô hình phân bố MARINE đã giải quyết triệt để bài toán cung cấp lưu lượng nhập dòng bên dọc theo hai bờ sông chính từ các tiểu lưu vực sườn dốc cao.

Dữ liệu mô phỏng trong nghiên cứu có thể được trực quan hóa tối ưu thông qua các đồ thị phân bố không gian lưu lượng gia nhập (bản đồ Raster 2D) kết hợp với các bảng tra cứu thủy lực từng mặt cắt và biểu đồ đường quá trình lưu lượng theo thời gian (Q-t) tại tuyến đập Hòa Bình. Đồ thị so sánh giữa đường quá trình thực đo và tính toán trước - sau khi áp dụng lọc Kalman thể hiện rõ khả năng hấp thụ nhiễu quan trắc và bù đắp sai số pha của đỉnh lũ, khẳng định tính ưu việt của phương pháp luận tích hợp.

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

Nhằm nâng cao hiệu quả phòng chống thiên tai và tối ưu hóa điều hành hệ thống hồ chứa thủy điện trên lưu vực sông Đà, nghiên cứu đưa ra 4 khuyến nghị then chốt:

Tự động hóa chuỗi truyền nhận dữ liệu khí tượng thủy văn thời gian thực: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần chủ trì đầu tư nâng cấp mạng lưới 36 trạm đo mưa tự động và truyền dữ liệu trực tuyến về máy chủ mô hình MARINE-IMech1D với tần suất cập nhật 15 phút một lần, hoàn thành trong giai đoạn 2026 - 2027 nhằm đáp ứng yêu cầu dự báo lũ cực ngắn.

Nâng cấp độ phân giải mô hình không gian và bản đồ thổ nhưỡng: Viện Cơ học phối hợp cùng Cục Quản lý Đê điều và Phòng chống Thiên tai thực hiện cập nhật bản đồ DEM từ độ phân giải 50 m xuống 10 m bằng công nghệ viễn thám LiDAR, đồng thời chi tiết hóa 29 phân lớp đất trên toàn bộ 26.800 km² diện tích lưu vực trong vòng 24 tháng tới để giảm sai số tính toán thấm Green-Ampt xuống dưới 8%.

Mở rộng áp dụng thuật toán lọc Kalman mở rộng hoặc lọc Kalman tổ hợp: Tập đoàn Điện lực Việt Nam và Viện Khoa học Công nghệ Việt Nam cần phối hợp nghiên cứu triển khai thuật toán lọc ngẫu nhiên nâng cao cho toàn bộ 121 mặt cắt sông chính và các hồ chứa thượng nguồn, hoàn thiện trước mùa lũ năm 2028 nhằm tự động hóa quá trình hiệu chỉnh hệ số nhám lòng dẫn.

Xây dựng quy trình vận hành liên hồ chứa thời gian thực dựa trên dự báo sớm 48 giờ: Ban Chỉ đạo Quốc gia về Phòng chống Thiên tai cần ban hành khung hướng dẫn kỹ thuật áp dụng công cụ mô hình kết nối để điều hành liên hồ chứa Sơn La - Hòa Bình - Tuyên Quang, hướng tới mục tiêu cắt giảm từ 10% đến 15% đỉnh lũ cho hạ du đồng bằng sông Hồng và gia tăng hiệu suất phát điện từ 3% đến 5% hàng năm.

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

Nghiên cứu mang lại giá trị thực tiễn và học thuật cao cho 4 nhóm đối tượng chính:

Kỹ sư điều hành và quản lý hồ chứa thủy điện: Cán bộ kỹ thuật tại Công ty Thủy điện Hòa Bình và Thủy điện Sơn La có thể khai thác mô hình để tính toán chính xác lưu lượng nước về hồ trước 48 giờ, hỗ trợ xây dựng kịch bản xả lũ, tích nước và phát điện an toàn, hiệu quả.

Chuyên gia dự báo khí tượng thủy văn: Đội ngũ dự báo viên tại các Trung tâm Dự báo Khí tượng Thủy văn có thể sử dụng kết quả và quy trình ghép nối MARINE-IMech1D như một công cụ nghiệp vụ đáng tin cậy để nâng cao chất lượng bản tin cảnh báo lũ trên hệ thống sông Đà và sông Hồng.

Cơ quan quản lý nhà nước về phòng chống thiên tai và đê điều: Các nhà hoạch định chính sách tại Cục Quản lý Đê điều và Ban Chỉ đạo Phòng chống Thiên tai có thêm cơ sở khoa học để lập phương án phân lũ, chậm lũ và gia cố các trọng điểm đê điều xung yếu tại vùng đồng bằng Bắc Bộ.

Học viên cao học, nghiên cứu sinh ngành Cơ học Thủy khí và Kỹ thuật Tài nguyên Nước: Người làm nghiên cứu có thể sử dụng luận văn như tài liệu tham khảo chuyên sâu về phương pháp số giải phương trình Saint-Venant, kỹ thuật khử ma trận thưa tối ưu và phương pháp lọc Kalman trong mô hình hóa thủy văn.

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

Tại sao cần kết nối mô hình MARINE với mô hình IMech1D?

Mô hình MARINE đóng vai trò thu gom nước mưa và tính toán dòng chảy tràn trên 26.800 km² sườn dốc lưu vực, sau đó cung cấp giá trị lưu lượng gia nhập khu giữa dọc theo các triền sông. IMech1D tiếp nhận nguồn lưu lượng phụ này cùng lưu lượng biên thượng lưu để mô phỏng chính xác sự lan truyền sóng lũ qua 121 mặt cắt thủy lực phức tạp về đến hồ Hòa Bình.

Mô hình thấm Green-Ampt trong phần mềm MARINE hoạt động như thế nào?

Dựa trên định luật Darcy và thế áp mao dẫn, mô hình Green-Ampt chia quá trình thấm thành 3 giai đoạn rõ rệt qua 4 nhóm đất chính. Mô hình giúp xác định chính xác thời điểm đất đạt trạng thái bão hòa và bắt đầu sinh dòng chảy mặt khi cường độ mưa vượt qua tốc độ thấm tiềm năng, đảm bảo tính chân thực vật lý.

Thuật toán khử nút trong IMech1D mang lại lợi thế gì vượt trội?

Thuật toán tự động tìm và khử đuổi các nút có bậc liên kết thấp trước khi khử Gauss, giúp giải nhanh hệ phương trình thủy lực đại số tuyến tính quy mô lớn. Đối với mạng lưới 2093 đoạn sông và 1103 nút, thuật toán giúp giảm thời gian tính toán gần 7 lần, từ 170 phút xuống còn 32 phút 50 giây.

Kỹ thuật lọc Kalman nâng cao chất lượng dự báo lưu lượng ra sao?

Lọc Kalman thực hiện ước lượng tiên nghiệm trạng thái dòng chảy và kết hợp hiệu chỉnh hậu nghiệm dựa trên độ lệch giữa giá trị mô phỏng với dữ liệu thực đo tại các trạm kiểm tra. Nhờ đó, mô hình tự động triệt tiêu các sai số pha và biên độ ngay cả khi cường suất lũ tăng vọt tới 77,4 cm/h.

Luận văn xử lý bài toán thiếu hụt số liệu mưa không gian như thế nào?

Nghiên cứu tích hợp chuỗi dữ liệu từ 36 trạm quan trắc khí tượng thủy văn, sử dụng phương pháp đa giác Thiessen để phân chia vùng ảnh hưởng không gian dạng lưới Raster 50 m x 50 m. Cách tiếp cận này phản ánh chân thực các tâm mưa lớn tại vùng núi cao Hoàng Liên Sơn với lượng mưa trên 2500 mm/năm.

Kết luận

  • Xây dựng thành công bộ công cụ tích hợp liên hoàn giữa mô hình thủy văn phân bố MARINE và mô hình thủy lực một chiều IMech1D cho lưu vực sông Đà 26.800 km².
  • Chuẩn hóa toàn diện chuỗi dữ liệu không gian GIS độ phân giải 50 m x 50 m, tích hợp 10 tiểu lưu vực, 4 nhóm đất, 13 loại thảm phủ và 36 vùng ảnh hưởng mưa Thiessen.
  • Đột phá về mặt thuật toán số trị với kỹ thuật khử nút mạng sông, rút ngắn thời gian giải hệ phương trình 2093 đoạn sông xuống 32 phút 50 giây, giảm gần 7 lần.
  • Ứng dụng thành công kỹ thuật lọc Kalman rời rạc giúp hiệu chỉnh liên tục sai số dòng chảy, nâng cao độ chính xác dự báo lưu lượng về hồ Hòa Bình trước 48 giờ.
  • Kiểm định xuất sắc qua 12 bài toán mẫu quốc tế và chứng minh hiệu quả thực tế qua mô phỏng trận lũ năm 2006 cũng như tác nghiệp thực tế mùa lũ năm 2009.

Trong giai đoạn tiếp theo từ năm 2026 đến năm 2027, hệ thống cần được tiếp tục mở rộng tích hợp nguồn dữ liệu mưa dự báo số trị phân giải cao và phát triển giao diện dự báo tự động. Hãy ứng dụng ngay giải pháp kết nối mô hình thủy văn - thủy lực tiên tiến này để nâng tầm hiệu quả quản lý lưu vực và điều hành hồ chứa an toàn trong tương lai!