Tổng quan nghiên cứu

Lưu vực sông Hồng – Thái Bình là lưu vực sông liên quốc gia lớn thứ hai tại Việt Nam với tổng diện tích tự nhiên đạt khoảng 169.000 km2, trong đó phần diện tích thuộc lãnh thổ Việt Nam chiếm 87.840 km2 (khoảng 51,3%), phần lãnh thổ Trung Quốc chiếm 81.200 km2 (khoảng 48%) và Lào chiếm 1.100 km2 (khoảng 0,7%). Đây là nơi sinh sống của hơn 29,26 triệu người trải dài trên 26 tỉnh, thành phố, tập trung mật độ dân cư đồng bằng Bắc Bộ lên tới 1.000 người/km2. Trong lịch sử 100 năm qua, đồng bằng sông Hồng – Thái Bình đã phải gánh chịu 26 trận lũ lớn, điển hình là thảm họa vỡ đê vào tháng 8 năm 1971 khiến khoảng 100.000 người thiệt mạng. Những đợt mưa lũ cực đoan năm 1986 với lượng mưa 300 mm đến 400 mm, năm 1996 với đỉnh lũ Hà Nội đạt 12,43 m vượt báo động 3 là 0,93 m trong 6 ngày, cùng các trận ngập úng lịch sử năm 2008 và 2010 (khiến 46 người chết và 21 người mất tích) đã cho thấy tính cấp thiết của công tác phòng chống thiên tai.

Trước yêu cầu thực tiễn đó, nghiên cứu được triển khai nhằm xây dựng một công nghệ dự báo lũ tác nghiệp hoàn chỉnh, nâng cao độ tin cậy và kéo dài thời gian dự kiến cảnh báo từ 24 giờ đến 48 giờ. Đề tài tập trung giải quyết bài toán thủy lực phức tạp khi hệ thống xuất hiện thêm các công trình hồ chứa lớn như Tuyên Quang (năm 2007) và Sơn La (năm 2011), đồng thời xử lý các điểm nghẽn về thiếu hụt dữ liệu quan trắc từ thượng nguồn Trung Quốc. Phạm vi nghiên cứu bao quát mạng lưới sông Hồng – Thái Bình với 25 sông chính, tập trung vào chuỗi số liệu mùa lũ liên tục từ ngày 15 tháng 6 đến ngày 15 tháng 9 giai đoạn 2006–2012. Kết quả nghiên cứu mang lại ý nghĩa kinh tế – xã hội to lớn, hỗ trợ đắc lực cho các cơ quan quản lý đê điều và vận hành liên hồ chứa giảm thiểu rủi ro ngập lụt cho vùng kinh tế trọng điểm Bắc Bộ.

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 tổ hợp ba mô hình toán thủy văn – thủy lực hiện đại nhằm mô phỏng toàn diện chu trình dòng chảy trên lưu vực:

  • Mô hình thủy văn MIKE NAM: Là mô hình khái niệm thông số tập trung mô phỏng quá trình mưa – dòng chảy liên tục dựa trên phương trình cân bằng nước qua 4 bể chứa thẳng đứng: bể chứa mặt (dung tích tối đa $U_{max}$), bể chứa tầng rễ cây (dung tích $L_{max}$), bể ngầm tầng trên ($BF_U$) và bể ngầm tầng dưới ($BF_L$). Các thông số động lực học như hệ số dòng chảy tràn ($CQOF$), hệ số dòng chảy sát mặt ($CQIF$), hằng số thời gian tập trung nước ($CK_{1,2}$) và ngưỡng dòng chảy mặt ($TOF$) được sử dụng để tính toán lưu lượng lũ gia nhập từ các tiểu lưu vực với ngưỡng dòng chảy tối thiểu đạt 0,4 mm/h.
  • Mô hình thủy lực 1 chiều MIKE 11 HD: Dựa trên hệ phương trình vi phân phi tuyến Saint-Venant mô tả dòng chảy một chiều không ổn định trong lòng dẫn (gồm phương trình liên tục bảo toàn khối lượng và phương trình động lượng bảo toàn động lượng). Mô hình sử dụng sơ đồ sai phân hữu hạn ẩn 6 điểm Abbott-Ionescu trên lưới so le với điều kiện ổn định Courant dao động từ 10 đến 15, cho phép xử lý ma trận phân nhánh phức tạp tại các hợp lưu và phân lưu.
  • Mô hình thủy lực 2 chiều MIKE 21: Ứng dụng giải hệ phương trình nước nông 2 chiều bằng thuật toán ADI (Alternating Direction Implicit) nhằm mô phỏng chính xác trường sóng triều truyền từ Vịnh Bắc Bộ vào các cửa sông ven biển, cung cấp điều kiện biên mực nước hạ lưu.

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

Nghiên cứu sử dụng nguồn dữ liệu thực đo phong phú với cỡ mẫu bao gồm chuỗi quan trắc khí tượng thủy văn liên tục qua 6 mùa lũ (2006–2011) để hiệu chỉnh, kiểm định và mùa lũ năm 2012 để thử nghiệm tác nghiệp. Cỡ mẫu hình học bao gồm mạng sơ đồ tính với 52 nhánh sông, 792 mặt cắt trắc ngang địa hình được chuẩn hóa, kết hợp thông số dung tích phòng lũ của các hồ chứa thủy điện Sơn La (9,26 tỷ m3), Hòa Bình (9,45 tỷ m3) và Tuyên Quang (2,26 tỷ m3).

Phương pháp chọn mẫu phân tầng theo không gian và thời gian được lựa chọn nhằm tối ưu hóa việc phân chia các lưu vực bộ phận theo mạng lưới trạm đo mưa thực tế. Phương pháp phân tích số học kết hợp kiểm tra chéo (cross-validation) qua các năm độc lập và phân cấp mực nước (thấp, trung bình, cao) được lựa chọn vì tính chất thủy lực lòng dẫn sông Hồng thay đổi mạnh theo mùa và chịu tác động của hiện tượng nước vật. Mức độ phù hợp của mô hình được định lượng thông qua chỉ số hiệu quả Nash-Sutcliffe ($NASH$) và sai số đỉnh lũ ($ΔH_{max}$, $ΔQ_{max}$).

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

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

Quá trình nghiên cứu và thử nghiệm tính toán đã ghi nhận nhiều kết quả vượt trội:

  • Tối ưu hóa mô hình mưa – dòng chảy NAM: Việc phân chia lại mạng lưới lưu vực bộ phận giúp mô phỏng chính xác lưu lượng nhập lưu về các hồ chứa lớn. Chỉ số $NASH$ trong giai đoạn hiệu chỉnh và kiểm định tại các trạm thượng nguồn hồ Sơn La, Hòa Bình và Tuyên Quang năm 2010–2011 đều đạt mức cao từ 0,82 đến 0,91, với sai số lưu lượng đỉnh lũ duy trì dưới 8,5%.
  • Nâng cao độ chính xác mô phỏng thủy lực MIKE 11: Kết quả kiểm tra chéo trong giai đoạn 2006–2011 tại 3 trạm khống chế trọng điểm (Sơn Tây, Hà Nội, Phả Lại) đạt chỉ số $NASH$ dao động từ 0,85 đến 0,94. Khi áp dụng bộ hệ số nhám Manning phân cấp theo từng cấp mực nước thực tế, sai số mực nước đỉnh lũ giảm từ 0,35 m xuống dưới 0,15 m, cải thiện 57% độ lệch so với phương pháp dùng hệ số nhám cố định.
  • Hiệu quả cập nhật sai số tác nghiệp thời gian thực: Thuật toán hiệu chỉnh sai số biên và sai số pha giúp triệt tiêu độ trễ truyền lũ. Dự báo thử nghiệm mùa lũ năm 2011 và 2012 với thời gian dự kiến 24 giờ cho thấy sai số mực nước tại trạm Hà Nội giảm từ 0,22 m xuống còn 0,11 m (giảm 50% sai số), trong khi dự báo 48 giờ đạt độ tin cậy trên 84%.
  • Kết nối biên truyền triều MIKE 21: Mô hình 2 chiều khu vực Vịnh Bắc Bộ đã tái hiện chính xác dao động triều tại các cửa sông Đáy, Ba Lạt, Lạch Giang, giúp giảm 18% sai số lan truyền mực nước ngược dòng trong các đợt triều cường kết hợp lũ lớn năm 1996 và 2002.

Thảo luận kết quả

Sự biến đổi lòng dẫn và tác động điều tiết của hệ thống bậc thang thủy điện sông Đà (đặc biệt khi hồ Sơn La phát điện năm 2011) đã làm thay đổi hoàn toàn đường quan hệ lưu lượng – mực nước $Q = f(H)$ tại trạm thủy văn Lai Châu và Tạ Bú. Hiện tượng nước vật từ hồ chứa đòi hỏi mô hình thủy lực phải cập nhật liên tục đường quá trình xả lũ thực tế thay vì dựa vào các hàm ngoại suy truyền thống.

Trong các báo cáo chuyên ngành, dữ liệu thủy văn được trực quan hóa tối ưu qua biểu đồ kết hợp đường quá trình mực nước thực đo và tính toán ($H - t$) theo thời gian thực tại các trạm khống chế, đi kèm bảng đối chiếu chỉ số $NASH$, sai số pha đỉnh lũ ($Δt$) và sai số biên độ ($ΔH$). So với các công trình nghiên cứu trước đây chỉ dừng lại ở mô phỏng tĩnh hoặc nghiên cứu quy hoạch phân lũ sông Đáy, mô hình đề xuất đã tích hợp thành công mô-đun cập nhật sai số trực tiếp, giải quyết triệt để bài toán thiếu hụt thông tin quan trắc từ 81.200 km2 diện tích lưu vực phía Trung Quốc.

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

Nhằm hoàn thiện công nghệ dự báo lũ và đưa vào ứng dụng thực tiễn bền vững, nghiên cứu đề xuất 4 giải pháp trọng tâm:

  1. Hiện đại hóa mạng lưới đo đạc tự động: Lắp đặt bổ sung 45 trạm đo mưa tự động và 15 trạm quan trắc mực nước số hóa trên các nhánh sông Đà, sông Thao và sông Lô trước năm 2028. Cơ quan thực hiện: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia, nhằm mục tiêu giảm sai số dự báo mưa đầu vào từ 20% xuống dưới 8%.
  2. Xây dựng cổng dữ liệu liên hồ chứa theo thời gian thực: Tự động hóa việc truyền nhận thông tin lưu lượng xả và mực nước hồ Sơn La, Hòa Bình, Tuyên Quang với tần suất cập nhật 1 giờ/lần trong giai đoạn 2027–2029. Cơ quan chủ trì: Tập đoàn Điện lực Việt Nam phối hợp với Ban Chỉ đạo Quốc gia về Phòng, chống thiên tai.
  3. Cập nhật định kỳ trắc ngang lòng dẫn và địa hình ngập lụt: Tổ chức đo đạc lại 792 mặt cắt sông định kỳ 3 năm/lần bằng công nghệ đo sâu Sonar đa tia nhằm cập nhật kịp thời hiện tượng xói lở, bồi lắng và khai thác cát. Đơn vị thực hiện: Cục Đường thủy nội địa Việt Nam và Viện Quy hoạch Thủy lợi.
  4. Mở rộng chia sẻ thông tin thủy văn xuyên biên giới: Thúc đẩy đàm phán song phương nhằm tiếp nhận dữ liệu xả lũ thời gian thực từ các hồ chứa trên lãnh thổ Trung Quốc, hướng tới kéo dài thời gian dự kiến dự báo lũ lên 72 giờ vào năm 2030 dưới sự điều phối của Bộ Tài nguyên và Môi trường.

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

Công trình nghiên cứu mang giá trị học thuật và thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng:

  1. Dự báo viên khí tượng thủy văn: Tham khảo quy trình thiết lập mô hình tích hợp MIKE NAM – MIKE 11 HD, phương pháp hiệu chỉnh hệ số nhám theo cấp mực nước và thuật toán cập nhật sai số tác nghiệp thời gian thực phục vụ bản tin dự báo hàng ngày.
  2. Kỹ sư vận hành công trình thủy điện và hồ chứa: Ứng dụng các kịch bản diễn toán lũ về hồ để xây dựng phương án xả lũ hợp lý, tối ưu hóa dung tích đón lũ và đảm bảo an toàn tuyệt đối cho công trình đầu mối.
  3. Cán bộ quản lý đê điều và phòng chống thiên tai: Sử dụng kết quả dự báo mực nước tại các trạm trọng điểm (Sơn Tây, Hà Nội, Phả Lại) để chủ động lập kế hoạch hộ đê, vận hành hệ thống phân lũ sông Đáy và phương án sơ tán dân cư vùng trũng thấp.
  4. Giảng viên, học viên cao học và sinh viên ngành Tài nguyên nước: Khai thác làm tài liệu giảng dạy, nghiên cứu chuyên sâu về mô hình hóa thủy lực sông ngòi phức tạp, phân tích bất định và phương pháp giải số phương trình thủy động lực học.

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

1. Điểm mới nổi bật của luận văn so với các nghiên cứu dự báo lũ sông Hồng trước đây là gì?
Luận văn đã tích hợp chuỗi công nghệ liên hoàn từ mô hình mưa – dòng chảy (NAM), thủy lực 1D (MIKE 11) đến thủy lực 2D (MIKE 21), đồng thời cập nhật sự xuất hiện của 2 đại công trình hồ chứa Tuyên Quang và Sơn La cùng thuật toán hiệu chỉnh sai số tác nghiệp thời gian thực.

2. Tại sao cần phân chia hệ số nhám Manning theo từng cấp mực nước thay vì dùng một giá trị cố định?
Do đặc thù hình học lòng sông Hồng có bãi bồi rộng và độ gồ ghề thay đổi lớn giữa lòng chính và bãi sông. Việc phân cấp hệ số nhám theo 3 cấp mực nước giúp tăng độ chính xác tính toán, nâng chỉ số $NASH$ lên trên 0,90 và giảm sai số đỉnh lũ xuống dưới 0,15 m.

3. Mô hình xử lý thách thức thiếu hụt số liệu từ phần lưu vực phía Trung Quốc như thế nào?
Nghiên cứu kết hợp phương pháp hồi quy tương quan chuỗi thủy văn trạm biên giới với chức năng cập nhật sai số pha và sai số biên trong mô-đun dự báo của MIKE 11, giúp tự động điều chỉnh đường quá trình mực nước tính toán tiệm cận sát với giá trị thực đo.

4. Chỉ số Nash-Sutcliffe ($NASH$) đạt từ 0,85 đến 0,94 có ý nghĩa như thế nào trong dự báo thủy văn tác nghiệp?
Trong khoa học thủy văn, chỉ số $NASH > 0,75$ được xếp loại rất tốt. Kết quả đạt từ 0,85 đến 0,94 khẳng định mô hình mô phỏng cực kỳ chuẩn xác cả về biên độ đỉnh lũ, thời gian xuất hiện đỉnh và tổng lượng dòng chảy mùa lũ qua các năm kiểm định độc lập.

5. Thời gian dự kiến dự báo của công nghệ này đạt hiệu quả cao nhất trong bao lâu?
Công nghệ đạt độ tin cậy tối ưu trong khoảng thời gian dự báo 24 giờ đến 48 giờ. Sai số mực nước dự báo 24 giờ tại Hà Nội chỉ dao động trong khoảng 0,11 m đến 0,15 m, đáp ứng hoàn toàn tiêu chuẩn nghiệp vụ của cơ quan khí tượng thủy văn quốc gia.

Kết luận

  • Luận văn đã xây dựng thành công chuỗi công nghệ dự báo lũ tích hợp trên hệ thống sông Hồng – Thái Bình bằng bộ mô hình MIKE (NAM, 11 HD, 21), đáp ứng yêu cầu dự báo tác nghiệp thực tế.
  • Cập nhật hoàn chỉnh cơ sở dữ liệu hình học gồm 52 nhánh sông, 792 mặt cắt ngang và quy trình vận hành liên hồ chứa Sơn La, Hòa Bình, Tuyên Quang.
  • Đạt độ chính xác cao trong mô phỏng thủy văn và thủy lực với chỉ số hiệu quả $NASH$ từ 0,85 đến 0,94 tại các trạm khống chế trọng điểm vùng đồng bằng Bắc Bộ.
  • Thiết lập quy trình tự động cập nhật sai số biên và sai số pha, giúp giảm hơn 40% sai số dự báo mực nước thời hạn 24 giờ và 48 giờ.
  • Các cơ quan quản lý và đơn vị nghiên cứu cần tiếp tục ứng dụng, mở rộng thử nghiệm công nghệ này trong các mùa lũ tiếp theo giai đoạn 2026–2030 nhằm nâng cao năng lực ứng phó thiên tai toàn diện cho lưu vực.