Tổng quan nghiên cứu

Lưu vực sông Lam là hệ thống sông lớn thứ hai tại khu vực miền Trung Việt Nam với tổng diện tích tự nhiên đạt 27.200 km² và chiều dài dòng chính lên tới 531 km, trong đó phần diện tích thuộc lãnh thổ Việt Nam chiếm khoảng 17.730 km² trải rộng trên địa bàn hai tỉnh Nghệ An và Hà Tĩnh. Chế độ thủy văn tại khu vực này vô cùng khắc nghiệt khi lượng dòng chảy trong mùa lũ tập trung chiếm từ 61,1% đến 75,1% tổng lượng dòng chảy cả năm. Địa hình sườn dốc phía Tây đón gió trực tiếp từ biển kết hợp với dải đồng bằng hẹp ven biển làm cho lũ tập trung nhanh, biên độ lũ lớn và khả năng thoát nước kém. Những trận lũ lịch sử tiêu biểu như năm 1978 với mực nước đỉnh lũ tại trạm Nam Đàn lên tới 964 cm và trạm Linh Cảm đạt 775 cm đã gây ngập lụt nghiêm trọng, làm vỡ nhiều tuyến đê và phá hủy hạ tầng dân sinh.

Mục tiêu cốt lõi của nghiên cứu là đánh giá thực trạng diễn biến thủy văn dòng chảy lũ, ứng dụng công nghệ mô hình toán thủy văn - thủy lực hiện đại nhằm lượng hóa khả năng phòng chống lũ của hệ thống đê dọc sông Lam trước các kịch bản lũ lịch sử và lũ thiết kế. Nghiên cứu khai thác chuỗi số liệu đo đạc thực tế giai đoạn từ năm 1960 đến năm 2011 để phân tích không gian ngập lụt tại vùng hạ lưu. Kết quả tính toán cung cấp luận cứ khoa học thực tiễn giúp nâng cao chỉ số an toàn đê điều, hỗ trợ chính quyền địa phương giảm thiểu tới 40% nguy cơ ngập úng diện rộng và bảo vệ vững chắc tính mạng cùng tài sản của hơn 3 triệu người dân trong vùng dự án.

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 hệ thống lý thuyết thủy văn và thủy lực hiện đại, bao gồm lý thuyết hình thành dòng chảy mặt từ mưa rào và hệ phương trình vi phân Saint-Venant một chiều mô tả dòng chảy không ổn định trong mạng lưới sông kênh kết hợp hệ phương trình nước nông Navier-Stokes hai chiều mô tả quá trình lan truyền sóng lũ trên đồng bằng bãi tràn.

Mô hình nghiên cứu được xây dựng trên nền tảng bộ công cụ chuyên dụng của Viện Thủy động lực học Đan Mạch với cấu trúc tích hợp:

  • Mô hình mưa - dòng chảy thủy văn liên tục mô phỏng quá trình bốc thoát hơi nước, bổ cập nước ngầm và sinh dòng chảy mặt cho từng tiểu lưu vực.
  • Mô hình thủy lực một chiều tính toán chế độ dòng chảy không ổn định trong hệ thống mạng lưới sông chính và các phụ lưu cấp 1.
  • Mô hình thủy lực hai chiều phân tích chi tiết trường vận tốc và độ sâu ngập lụt trên bề mặt địa hình đồng bằng.
  • Mô hình tích hợp liên kết động giữa dòng chính một chiều và đồng bằng hai chiều qua các kết cấu bờ bao, đê ngăn lũ.

Các khái niệm then chốt được lượng hóa bao gồm: mực nước chân lũ và đỉnh lũ, mô đun dòng chảy đỉnh lũ lớn nhất, tần suất lũ thiết kế chu kỳ lặp lại 1% và 2%, hệ số nhám lòng dẫn Manning trong phạm vi từ 0,025 đến 0,035.

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

Nguồn dữ liệu đầu vào phục vụ mô hình bao gồm chuỗi quan trắc khí tượng thủy văn tại 33 trạm đo mưa với mật độ trung bình 537 km²/trạm và 10 trạm thủy văn cấp I, cấp II phân bố dọc sông Cả, sông La, sông Hiếu và sông Nậm Mộ. Dữ liệu hình học được xây dựng từ hơn 120 mặt cắt ngang lòng sông thực đo kết hợp mô hình số độ cao địa hình độ phân giải cao.

Phương pháp chọn mẫu chuỗi số liệu tập trung vào việc trích xuất các trận lũ lịch sử và lũ đặc biệt lớn điển hình gồm các năm 1978, 2002, 2005, 2007, 2010 và 2011. Các trận lũ này phản ánh đầy đủ các tổ hợp thời tiết cực đoan nguy hiểm nhất như bão đổ bộ kết hợp dải hội tụ nhiệt đới và không khí lạnh tăng cường.

Lý do lựa chọn phương pháp mô hình hóa tích hợp một chiều và hai chiều bắt nguồn từ tính chất phức tạp của địa hình vùng hạ lưu sông Lam. Mô hình cho phép mô phỏng chính xác sự tương tác dòng chảy tràn bờ khi mực nước vượt cao trình đê, khắc phục triệt để nhược điểm của các phương pháp thủy lực một chiều truyền thống. Quá trình hiệu chỉnh mô hình được thực hiện trên chuỗi dữ liệu các trận lũ năm 2002, 2007, 2010 và kiểm định độc lập với các trận lũ năm 2005, 2011, mang lại chỉ số tương quan Nash-Sutcliffe đạt trên 0,85 tại tất cả các trạm khống chế trọng điểm.

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

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

Quá trình mô phỏng và phân tích diễn biến lũ trên lưu vực sông Lam mang lại 4 phát hiện khoa học quan trọng:

Thứ nhất, lũ trên lưu vực sông Lam có tính phân hóa không gian và thời gian rất sâu sắc. Lượng dòng chảy lũ của 3 tháng liên tục lớn nhất chiếm từ 50% đến 60% tổng lượng năm, riêng tháng đỉnh lũ (tháng IX hoặc X) chiếm từ 20% đến 24% lượng nước cả năm. Phụ lưu sông La có mật độ lưới sông dày đặc từ 0,87 đến 0,91 km/km² cùng độ dốc lưu vực lên tới 28,2%, khiến thời gian tập trung lũ nhanh hơn 1,5 lần và có tính chất ác liệt hơn hẳn so với thượng lưu sông Cả có độ dốc từ 10,6% đến 20,1%.

Thứ hai, kết quả mô phỏng trận lũ lịch sử năm 1978 trên hệ thống đê hiện trạng xác định mực nước tại Nam Đàn vượt báo động III tới 2,64 m, gây ra hiện tượng tràn cục bộ tại hàng loạt vị trí xung yếu dọc đê Tả Lam và Hữu Lam với chiều dài ngập tràn vượt trên 15 km, dẫn đến diện tích ngập sâu trên 1,0 m ở vùng hạ du đạt hơn 35.000 ha.

Thứ ba, việc nâng cấp đê Tả Lam lên cấp II chống lũ tần suất 1% và kiên cố hóa 45,9 km mặt đê giúp bảo vệ an toàn cho thành phố Vinh và huyện Hưng Nguyên. Tuy nhiên, việc khép kín tuyến đê làm mực nước đỉnh lũ trong lòng dẫn tăng cục bộ từ 0,15 m đến 0,35 m tại các đoạn co hẹp ven sông.

Thứ tư, các tổ hợp mưa lớn có lượng mưa 7 ngày tích lũy từ 800 mm đến 1.400 mm ở đồng bằng là nguyên nhân trực tiếp kích hoạt các đỉnh lũ kép kéo dài từ 5 đến 9 ngày, gây quá tải toàn diện khả năng tiêu thoát tự nhiên ra biển.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến tình trạng ngập lụt nghiêm trọng vùng hạ du sông Lam là do địa hình đồi núi chiếm 80% diện tích lưu vực với tỷ lệ đất trống đồi núi trọc còn ở mức xấp xỉ 20%, làm suy giảm khả năng giữ nước tự nhiên của thảm phủ thực vật. Đồng thời, hệ thống hồ chứa thủy điện ở thượng nguồn chưa có đủ dung tích phòng lũ chuyên biệt dành riêng cho việc cắt giảm lũ hạ lưu.

So sánh với các nghiên cứu quy hoạch thủy lợi sông Cả năm 2004 và quản lý lũ năm 2012, nghiên cứu này đã tiến xa hơn khi lượng hóa cụ thể ranh giới và độ sâu ngập 2D thay vì chỉ tính toán đường mặt nước một chiều dọc sông. Dữ liệu nghiên cứu được biểu diễn trực quan qua các biểu đồ đường quá trình lưu lượng và mực nước thực đo so với tính toán tại các trạm Dừa, Linh Cảm, Nam Đàn, cùng hệ thống bảng tổng hợp diện tích phân vùng ngập úng theo các cấp độ sâu từ 0,5 m đến trên 2,0 m. Bản đồ ngập lụt số hóa thể hiện chi tiết vị trí các mặt cắt tràn đê xung yếu, cung cấp cơ sở tin cậy cho việc thiết kế cao trình đỉnh đê quy hoạch.

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

Dựa trên kết quả tính toán thủy lực và đánh giá rủi ro thiên tai, luận văn đưa ra 5 nhóm giải pháp đồng bộ:

  1. Nâng cấp và gia cố công trình đê điều: Ủy ban nhân dân hai tỉnh Nghệ An và Hà Tĩnh cùng Sở Nông nghiệp và Phát triển nông thôn tập trung triển khai nâng cấp tuyến đê Tả Lam thành đê cấp II (chống lũ thiết kế tần suất 1%), đắp cơ gia cố 12 km thân đê, xử lý triệt để 753 m nền đê xung yếu có nguy cơ sủi mạch đùn, và kiên cố hóa 81,3 km đê Hữu Lam và đê Bích Hào đạt chuẩn đê cấp III (chống lũ tần suất 2%) trong giai đoạn từ năm 2025 đến năm 2030.

  2. Phục hồi và phát triển rừng phòng hộ đầu nguồn: Chi cục Kiểm lâm phối hợp cùng chính quyền các huyện Kỳ Sơn, Tương Dương, Vũ Quang đẩy mạnh trồng mới tối thiểu 10.000 ha rừng mỗi năm, nâng tỷ lệ che phủ rừng toàn lưu vực sông Lam đạt trên 55% nhằm giảm thiểu 10% đến 15% lưu lượng đỉnh lũ đổ về hạ du.

  3. Hiện đại hóa hệ thống quan trắc và dự báo cảnh báo sớm: Đài Khí tượng Thủy văn khu vực Bắc Trung Bộ lắp đặt thêm 15 đến 20 trạm đo mưa tự động tại các khu vực vùng núi hiểm trở thượng nguồn sông Nậm Mộ và sông Ngàn Phố trước năm 2027, rút ngắn khoảng cách lưới trạm xuống dưới 300 km²/trạm và kéo dài thời gian dự báo cảnh báo lũ trước từ 24 đến 48 giờ.

  4. Vận hành tối ưu quy trình liên hồ chứa thủy điện: Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn tỉnh Nghệ An chỉ đạo các đơn vị quản lý hồ chứa Bản Vẽ, Khe Bố chủ động duy trì dung tích phòng lũ đón trước từ 15% đến 20% dung tích hữu ích khi có dự báo bão, nhằm cắt giảm đỉnh lũ hạ lưu an toàn.

  5. Thực hiện các biện pháp quản lý rủi ro phi công trình: Chính quyền các địa phương ven sông ban hành bản đồ phân vùng rủi ro ngập lụt tỷ lệ 1/25.000, cắm mốc cảnh báo ngập sâu tại 100% các khu dân cư trũng thấp và tổ chức diễn tập phương án sơ tán định kỳ hàng năm trước ngày 15 tháng 8.

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

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng chính:

  1. Cán bộ quản lý đê điều và quy hoạch thủy lợi: Khai thác các thông số tính toán mực nước lũ tần suất 1% và 2% cùng danh mục các đoạn đê xung yếu để lập dự án đầu tư nâng cấp công trình đê kè hạ lưu sông Lam đúng tiêu chuẩn kỹ thuật.

  2. Kỹ sư và chuyên gia mô hình hóa thủy văn - thủy lực: Sử dụng bộ thông số kiểm định của mô hình mưa - dòng chảy và thủy lực dòng chính làm tài liệu tham khảo chuẩn khi thiết lập các bài toán lan truyền sóng lũ phức tạp cho các lưu vực sông khu vực Bắc Trung Bộ.

  3. Ban Chỉ huy Phòng thủ dân sự các cấp tại Nghệ An và Hà Tĩnh: Vận dụng bản đồ ngập lụt mô phỏng ứng với các kịch bản lũ lịch sử và lũ thiết kế để xây dựng phương án cứu hộ cứu nạn, phân bố kho vật tư dự phòng và điều phối lực lượng sơ tán dân cư an toàn.

  4. Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Kỹ thuật Tài nguyên nước, Quản lý Thiên tai: Dùng làm tài liệu tham khảo học thuật chuyên sâu về phương pháp ứng dụng công nghệ mô hình toán tích hợp giải quyết bài toán quy hoạch kiểm soát lũ lưu vực sông.

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

  1. Mô hình toán tích hợp một chiều và hai chiều mang lại ưu thế gì so với mô hình một chiều truyền thống khi tính lũ sông Lam? Mô hình tích hợp giúp mô tả chính xác sự tương tác thủy lực giữa lòng dẫn chính 531 km và bãi tràn đồng bằng, xác định cụ thể tọa độ và lưu lượng tràn qua đê, trong khi mô hình một chiều đơn thuần không thể mô phỏng đúng phạm vi lan truyền ngập lụt hai chiều trên diện rộng.

  2. Tại sao lũ trên sông La lại dâng cao và nguy hiểm hơn so với thượng lưu sông Cả? Sông La có độ dốc bình quân lưu vực lên tới 22,6% - 28,2%, mật độ lưới sông đạt 0,87 - 0,91 km/km² cùng lòng sông ngắn và dốc, khiến thời gian tập trung dòng chảy sau mưa nhanh gấp gần 2 lần so với thượng nguồn sông Cả vốn có nhiều nếp uốn khúc.

  3. Hệ thống đê hạ lưu sông Lam chịu tải như thế nào khi gặp trận lũ lịch sử tương đương năm 1978? Khi gặp trận lũ năm 1978 với đỉnh lũ Nam Đàn đạt 964 cm, hệ thống đê hiện trạng bị tràn tại nhiều vị trí với tổng chiều dài tràn vượt 15 km, gây ngập sâu hơn 35.000 ha đất canh tác và khu dân cư hạ lưu nếu không được nâng cấp kịp thời.

  4. Việc nâng cấp tuyến đê Tả Lam lên cấp II mang lại hiệu quả bảo vệ ra sao? Tuyến đê Tả Lam sau khi được nâng chuẩn cấp II và kiên cố hóa 45,9 km mặt đê sẽ đảm bảo an toàn tuyệt đối trước các trận lũ thiết kế tần suất 1%, bảo vệ toàn diện trung tâm kinh tế thành phố Vinh và hàng trăm cơ sở hạ tầng trọng điểm ven sông.

  5. Giải pháp phi công trình nào cần được ưu tiên triển khai ngay trước mùa mưa bão? Ưu tiên hàng đầu là hiện đại hóa mạng lưới quan trắc tự động với việc bổ sung 15 đến 20 trạm đo mưa chuyên dùng và thiết lập quy chế phối hợp vận hành xả lũ liên hồ chứa theo thời gian thực, giúp giảm thiểu trên 70% các thiệt hại bất ngờ về người và tài sản.

Kết luận

  • Luận văn đã phân tích toàn diện các đặc trưng thủy văn và hình thái lũ trên lưu vực sông Lam diện tích 27.200 km², xác định dòng chảy 3 tháng mùa lũ chiếm tới 50% - 60% tổng lượng nước cả năm.
  • Thiết lập và kiểm định thành công bộ mô hình toán tích hợp thủy văn - thủy lực với chỉ số Nash-Sutcliffe đạt trên 0,85, tái hiện chính xác diễn biến trận lũ lịch sử năm 1978 và các đợt lũ lớn gần đây.
  • Đánh giá chi tiết năng lực phòng chống lũ của hệ thống đê điều hiện trạng, chỉ ra các vị trí co hẹp và nguy cơ tràn bờ khi xuất hiện lũ tần suất 1%.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, trọng tâm là kiên cố hóa 45,9 km đê Tả Lam, xử lý 753 m nền đê yếu và nâng độ che phủ rừng lên trên 55%.
  • Xác lập lộ trình hành động cụ thể cho giai đoạn 2025–2030 nhằm hoàn thiện hệ thống phòng thủ thiên tai bền vững thích ứng với biến đổi khí hậu.

Các cơ quan quản lý và đơn vị nghiên cứu có thể tiếp tục khai thác cơ sở dữ liệu mô hình toán của luận văn để tích hợp vào hệ thống dự báo ngập lụt thời gian thực, góp phần nâng cao hiệu quả phòng chống thiên tai cho khu vực Bắc Trung Bộ.