Tổng quan nghiên cứu

Lưu vực sông Mã là một trong những hệ thống sông lớn tại khu vực miền núi phía Bắc và Bắc Trung Bộ với tổng diện tích tự nhiên khoảng 28.640 km², chiều dài dòng chính đạt 512 km và độ dốc bình quân lưu vực lên tới 17,6%. Khu vực thượng lưu sông Mã thuộc địa bàn các tỉnh Điện Biên và Sơn La có địa hình núi cao hiểm trở, mức độ chia cắt sâu và mạng lưới sông suối dạng cành cây với mật độ 0,66 km/km². Hàng năm, khu vực này thường xuyên hứng chịu từ 2 đến 3 đợt lũ lớn trong mùa mưa (kéo dài từ tháng 5 đến tháng 10), chiếm tới 85% tổng lượng mưa cả năm. Do địa hình dốc, dòng chảy lũ tập trung rất nhanh với thời gian hình thành và truyền lũ chỉ trong khoảng 6 đến 18 giờ, gây ra ngập lụt nghiêm trọng, lũ quét và xói lở đất, đe dọa trực tiếp đến tính mạng cùng sinh kế của nhân dân ven sông.

Trước tình hình đó, công tác dự báo lũ tại trạm thủy văn Xã Là thuộc Đài Khí tượng Thủy văn khu vực Tây Bắc trước đây chủ yếu dựa vào các phương pháp thống kê tương quan truyền thống kết hợp kinh nghiệm dự báo viên. Các phương án này tồn tại hạn chế lớn là thời gian dự kiến ngắn (chỉ đạt từ 6 đến 12 giờ) và độ chính xác suy giảm mạnh khi gặp các trận mưa cực đoan có cường suất lớn. Mục tiêu cốt lõi của nghiên cứu là ứng dụng mô hình toán thủy văn MIKE - NAM nhằm thiết lập công cụ tính toán, mô phỏng và dự báo chi tiết toàn bộ quá trình dòng chảy lũ tại trạm thủy văn Xã Là. Nghiên cứu thực hiện xử lý chuỗi dữ liệu thực đo 23 năm liên tục (giai đoạn 1994 - 2016) và tiến hành dự báo thử nghiệm trong mùa lũ giai đoạn 2016 - 2017. Kết quả đạt được mang ý nghĩa thực tiễn to lớn khi giúp kéo dài thời gian dự kiến cảnh báo sớm lên 24 đến 48 giờ, góp phần giảm thiểu ước tính khoảng 30% đến 40% thiệt hại do thiên tai lũ lụt cho các địa phương vùng thượng nguồn sông Mã.

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 ứng dụng khung lý thuyết mô hình thủy văn quan niệm thông số tập trung MIKE - NAM, một mô đun mưa - dòng chảy tiêu chuẩn thuộc hệ thống mô hình MIKE 11 do Viện Thủy lực Đan Mạch (DHI) phát triển. Cấu trúc mô hình MIKE - NAM hoạt động dựa trên nguyên lý cân bằng nước liên tục thông qua hệ thống 4 bể chứa vật lý tương tác theo chiều thẳng đứng và các bể chứa tuyến tính nằm ngang:

  1. Bể chứa ẩm bề mặt: Đặc trưng bởi lượng trữ ẩm tối đa Umax, đại diện cho lượng nước do thảm phủ thực vật chặn lại, lượng nước điền trũng và lớp sát mặt đất. Khi dung tích bể vượt ngưỡng Umax, lượng nước thừa sẽ hình thành dòng chảy tràn bề mặt hoặc thấm xuống các tầng dưới.
  2. Bể chứa tầng rễ cây: Giới hạn bởi dung tích ẩm tối đa Lmax, phản ánh lớp đất mà thực vật có thể hút nước để thoát hơi. Tỷ số L/Lmax biểu thị trạng thái ẩm thực tế của lưu vực, chi phối trực tiếp lượng bốc thoát hơi nước thực tế Ea so với bốc thoát hơi tiềm năng Ep cũng như lượng nước thấm bổ sung cho tầng ngầm.
  3. Bể chứa nước ngầm: Bao gồm tầng ngầm trên và tầng ngầm dưới, tiếp nhận lượng nước thấm qua ngưỡng bổ sung TG để duy trì dòng chảy kiệt và dòng chảy ngầm cho lòng sông.
  4. Cơ chế diễn toán dòng chảy: Dòng chảy sát mặt QIF được kiểm soát bởi hằng số thời gian CKIF và ngưỡng dưới TIF; dòng chảy tràn QOF được xác định qua hệ số dòng chảy mặt CQOF (biến thiên từ 0 đến 1) và ngưỡng TOF. Các thành phần này được tổng hợp và diễn toán qua hai hồ chứa tuyến tính liên tiếp với hằng số thời gian CK1 và CK2 để tính ra lưu lượng dòng chảy tổng cộng tại cửa ra của lưu vực.

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

Nghiên cứu thu thập và xử lý toàn bộ chuỗi số liệu khí tượng thủy văn trong 23 năm liên tục từ năm 1994 đến 2016, bao gồm 56 trận lũ thực đo có biên độ mực nước từ 1,0 m đến trên 4,0 m tại trạm Xã Là. Phương pháp chọn mẫu được thực hiện theo nguyên tắc chọn mẫu có chủ đích, phân loại thành 3 nhóm cấp lũ rõ rệt: nhóm lũ lớn (biên độ lũ trên 3,0 m), nhóm lũ trung bình (biên độ 2,0 - 3,0 m) và nhóm lũ nhỏ (biên độ 1,0 - 2,0 m). Việc phân nhóm này đảm bảo bộ thông số mô hình được tối ưu hóa chính xác cho từng điều kiện khí quyển và độ ẩm mặt đệm khác nhau.

Nguồn dữ liệu đầu vào bao gồm: số liệu đo mưa từng giờ từ mạng lưới các trạm khí tượng thủy văn thượng lưu (Điện Biên, Tuần Giáo, Mường Ẳng, Sông Mã, Sốp Cộp); số liệu mực nước và lưu lượng thực đo tại trạm thủy văn Xã Là; mô hình số độ cao DEM và bản đồ phân bố thảm phủ thực vật. Lý do lựa chọn mô hình MIKE - NAM là khả năng mô phỏng liên tục chu trình thủy văn với yêu cầu dữ liệu đầu vào trung bình, đặc biệt phù hợp với các lưu vực sông miền núi có mật độ trạm đo thưa thớt như vùng Tây Bắc. Các chỉ tiêu đánh giá bao gồm hệ số Nash-Sutcliffe (NASH), sai số lưu lượng đỉnh lũ (dQmax), sai số thời gian xuất hiện đỉnh (dt) và sai số tổng lượng lũ (dW). Timeline nghiên cứu phân định rõ ràng: chuỗi số liệu 1994 - 2015 phục vụ phân tích tương quan và hiệu chỉnh, kiểm định bộ thông số; chuỗi số liệu mùa lũ 2016 - 2017 được dùng độc lập để thử nghiệm dự báo nghiệp vụ.

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

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

Phân tích hình thế khí tượng synop trên 56 trận lũ trong 23 năm đã chỉ ra 3 nhóm nguyên nhân gây mưa lũ chủ đạo trên thượng lưu sông Mã: hoàn lưu bão và rãnh áp thấp bị nén chiếm tỷ lệ lớn nhất với 45% (25 trận); dải hội tụ nhiệt đới kết hợp các hình thế thời tiết khác chiếm 30% (17 trận); xoáy thấp phát triển từ mặt đất lên độ cao 5000 m chiếm 25% (14 trận). Đặc biệt, khoảng 42% số trận lũ lớn nhất trong chuỗi lịch sử tập trung vào tháng 8 hàng năm với lượng mưa ngày cực đại đạt từ 200 mm đến 300 mm.

Quá trình hiệu chỉnh bộ thông số mô hình MIKE - NAM trên các trận lũ lịch sử (điển hình như trận lũ ngày 18-28/7/2011 và 31/7-06/8/2015) đạt độ chính xác rất cao: hệ số Nash-Sutcliffe dao động từ 0,82 đến 0,91; sai số lưu lượng đỉnh lũ dQmax luôn duy trì dưới mức 10%; sai số thời gian xuất hiện đỉnh dt chỉ chênh lệch từ 1 đến 2 giờ.

Kết quả kiểm định độc lập trên các trận lũ lớn và trung bình (như trận lũ ngày 24-30/9/2008 và 15-24/8/2012) khẳng định tính vững chắc của mô hình với hệ số Nash-Sutcliffe đạt từ 0,80 đến 0,88, sai số tổng lượng lũ dW dưới 8,5%. Thử nghiệm dự báo tác nghiệp trong mùa lũ giai đoạn 2016 - 2017 chứng minh mô hình MIKE - NAM đạt tỷ lệ dự báo đủ tiêu chuẩn trên 85%, vượt trội hoàn toàn so với tỷ lệ chỉ đạt 63% đến 67% của các phương pháp tương quan thống kê truyền thống, đồng thời giảm sai số biên độ đỉnh lũ trung bình khoảng 18% đến 22%.

Thảo luận kết quả

Độ chính xác vượt trội của mô hình MIKE - NAM bắt nguồn từ việc mô phỏng đúng bản chất vật lý của quá trình bão hòa ẩm đất trên nền địa hình có độ dốc cao 17,6% của thượng lưu sông Mã. Khi thảm phủ thực vật và tầng rễ cây đạt trạng thái bão hòa ẩm (tỷ số L/Lmax xấp xỉ 1), cơ chế dòng chảy tràn bề mặt và dòng chảy sát mặt được kích hoạt mạnh mẽ, tạo nên đỉnh lũ nhọn và sườn lũ lên dốc đứng.

Diễn biến thủy văn được thể hiện trực quan qua các biểu đồ đường quá trình lưu lượng theo thời gian (hydrograph) so sánh giữa đường thực đo và tính toán, kết hợp đường lũy tích tổng lượng lũ và bảng đối chiếu sai số định lượng cho từng trận lũ. Các biểu đồ mô phỏng cho thấy đường lưu lượng tính toán bám sát chặt chẽ biên độ đỉnh thực tế cũng như nhánh rút của dòng chảy lũ.

Khi so sánh với các nghiên cứu áp dụng mô hình TANK hoặc MARINE trên lưu vực sông Đà và sông Hồng, mô hình MIKE - NAM tại trạm Xã Là thể hiện khả năng thích ứng cao hơn trong điều kiện mạng lưới đo mưa thưa thớt, không bị tích lũy sai số theo thời gian và giải quyết triệt để độ trễ dự báo. Kết quả này chứng minh tính khả thi của việc chuyển đổi từ các công cụ dự báo điểm thủ công sang hệ thống mô phỏng số trị tự động hóa hoàn toàn trong nghiệp vụ khí tượng thủy văn.

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

  1. Hiện đại hóa mạng lưới trạm đo mưa tự động: Sở Tài nguyên và Môi trường tỉnh Sơn La và Điện Biên cần phối hợp cùng Đài Khí tượng Thủy văn khu vực Tây Bắc khảo sát, lắp đặt bổ sung từ 8 đến 10 trạm đo mưa tự động tại các khu vực thượng nguồn trọng điểm như Tuần Giáo, Điện Biên Đông và Sốp Cộp trước năm 2028, nhằm giảm sai số không gian của lượng mưa bình quân lưu vực xuống dưới mức 5%.
  2. Tích hợp dữ liệu dự báo mưa định lượng số trị: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia chủ trì triển khai kết nối tự động các trường mưa dự báo số trị độ phân giải cao (từ các mô hình WRF, IFS) làm đầu vào trực tiếp cho mô hình MIKE - NAM trong giai đoạn 2027 - 2028, giúp kéo dài thời gian dự kiến cảnh báo lũ từ 12 giờ lên 24 đến 48 giờ.
  3. Chuẩn hóa quy trình dự báo tác nghiệp nghiệp vụ: Ban lãnh đạo Đài Khí tượng Thủy văn khu vực Tây Bắc cần ban hành quy chuẩn kỹ thuật vận hành mô hình MIKE - NAM cho trạm Xã Là, tổ chức tập huấn chuyên sâu để 100% dự báo viên làm chủ kỹ năng hiệu chỉnh thông số theo từng cấp lũ trong vòng 6 tháng kể từ khi nghiệm thu hệ thống.
  4. Xây dựng bản đồ số phân vùng cảnh báo ngập lụt: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn các tỉnh Điện Biên và Sơn La chủ động ứng dụng công nghệ GIS kết hợp kết quả mô phỏng lưu lượng đỉnh lũ của mô hình MIKE - NAM để lập bản đồ ngập lụt chi tiết đến cấp xã ven sông Mã, phấn đấu giảm thiểu ít nhất 25% thiệt hại về tài sản của nhân dân trong giai đoạn 2028 - 2030.

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

  1. Dự báo viên và kỹ thuật viên khí tượng thủy văn: Khai thác trực tiếp bộ thông số tối ưu và quy trình vận hành mô hình MIKE - NAM để áp dụng vào công tác ra bản tin cảnh báo lũ hàng ngày tại Đài Khí tượng Thủy văn khu vực Tây Bắc và các đài khu vực lân cận.
  2. Cơ quan quản lý nhà nước và Ban Chỉ huy Phòng chống thiên tai: Sử dụng các phân tích quy luật synop và kịch bản ngập lụt trong luận văn làm căn cứ khoa học để xây dựng phương án sơ tán dân cư, gia cố đê kè và điều hành liên hồ chứa giảm lũ an toàn cho hạ du.
  3. Giảng viên, nhà nghiên cứu và học viên cao học ngành Thủy văn học: Sử dụng toàn bộ phương pháp luận xử lý chuỗi số liệu 23 năm và kỹ thuật hiệu chuẩn mô hình dòng chảy miền núi làm tài liệu tham khảo giảng dạy, nghiên cứu mở rộng trên các lưu vực sông xuyên biên giới.
  4. Kỹ sư tư vấn quy hoạch xây dựng và thủy điện: Tham khảo các đặc trưng lưu lượng lũ thiết kế, thời gian truyền lũ 6 đến 18 giờ và đường lũy tích tổng lượng lũ để tính toán khẩu độ thoát nước cho cầu đường, đập dâng và các công trình thủy lợi ven sông Mã.

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

Mô hình MIKE - NAM có ưu điểm gì vượt trội so với các phương pháp dự báo truyền thống tại thượng lưu sông Mã? Mô hình MIKE - NAM mô phỏng liên tục toàn bộ đường quá trình lũ và phản ánh chính xác trạng thái ẩm đất, thay vì chỉ ước tính đỉnh lũ đơn lẻ như phương pháp tương quan thống kê. Kết quả thử nghiệm mùa lũ 2016 - 2017 cho thấy mô hình đạt tỷ lệ dự báo chuẩn trên 85%, cao hơn phương pháp truyền thống khoảng 18% đến 22% và kéo dài thời gian dự kiến từ 12 giờ lên 24 đến 48 giờ.

Những hình thế thời tiết nào gây ra các đợt mưa lũ nguy hiểm nhất tại thượng nguồn sông Mã? Phân tích 56 trận lũ trong 23 năm (1994 - 2016) khẳng định hoàn lưu bão kết hợp rãnh áp thấp bị nén là nguyên nhân hàng đầu, chiếm 45% tổng số trận lũ. Kế tiếp là dải hội tụ nhiệt đới chiếm 30% và xoáy thấp chiếm 25%. Các hình thế này tạo ra lượng mưa 100 đến 300 mm, sinh ra lũ lớn có biên độ trên 3,0 m tại trạm Xã Là.

Tại sao mô hình MIKE - NAM lại đặc biệt phù hợp với lưu vực sông miền núi Tây Bắc? Lưu vực thượng nguồn sông Mã có độ dốc bình quân lớn 17,6% và mạng lưới trạm đo thưa thớt. MIKE - NAM là mô hình thông số tập trung với cấu trúc toán học tinh gọn, nhu cầu dữ liệu trung bình nhưng mô tả chính xác cơ chế dòng chảy mặt và dòng sát mặt, giúp hệ số Nash-Sutcliffe luôn đạt mức cao từ 0,82 đến 0,91.

Hệ số Nash-Sutcliffe đạt từ 0,80 đến 0,91 mang ý nghĩa khoa học như thế nào? Chỉ số Nash-Sutcliffe là thước đo chuẩn quốc tế về mức độ tương thích giữa số liệu tính toán và thực đo thủy văn. Mức giá trị từ 0,80 đến 0,91 cùng sai số tổng lượng lũ dưới 8,5% khẳng định mô hình có độ tin cậy rất cao, phản ánh chính xác cả lưu lượng đỉnh, thời gian xuất hiện lũ và nhánh rút của các trận lũ lịch sử.

Làm thế nào để duy trì và nâng cao độ chính xác của mô hình trong điều kiện biến đổi khí hậu? Giải pháp then chốt là tích hợp chuỗi dữ liệu mưa số trị độ phân giải cao từ các mô hình khí hậu toàn cầu, kết hợp bổ sung từ 8 đến 10 trạm đo mưa tự động tại vùng núi cao. Việc cập nhật liên tục thông số mặt đệm sẽ giúp kiểm soát sai số đỉnh lũ dưới 5% trong điều kiện thời tiết cực đoan.

Kết luận

  • Ứng dụng thành công mô hình MIKE - NAM: Luận văn đã xây dựng hoàn chỉnh công cụ toán học mô phỏng liên tục dòng chảy lũ từ mưa cho lưu vực thượng lưu sông Mã tại trạm thủy văn Xã Là.
  • Độ tin cậy và độ chính xác xuất sắc: Kết quả hiệu chỉnh và kiểm định đạt hệ số Nash-Sutcliffe từ 0,80 đến 0,91, sai số lưu lượng đỉnh lũ dưới 10% và tỷ lệ dự báo nghiệp vụ mùa lũ 2016 - 2017 đạt trên 85%.
  • Nâng cao năng lực cảnh báo sớm: Kéo dài thời gian dự kiến dự báo từ 6 - 12 giờ của phương án truyền thống lên 24 - 48 giờ khi kết hợp nguồn dữ liệu mưa định lượng.
  • Lộ trình phát triển giai đoạn 2027 - 2030: Đẩy mạnh tích hợp mô hình dự báo mưa số trị, mở rộng mạng lưới quan trắc tự động và xây dựng bản đồ số cảnh báo ngập lụt chi tiết cho các huyện ven sông.
  • Khuyến nghị ứng dụng thực tiễn: Các cơ quan quản lý khí tượng thủy văn và Ban Chỉ huy Phòng chống thiên tai khu vực Tây Bắc cần sớm đưa bộ thông số mô hình vào quy trình dự báo tác nghiệp chính thức nhằm bảo vệ an toàn tính mạng và tài sản cho cộng đồng.