Tổng quan nghiên cứu

Hệ thống đê điều tại Việt Nam đóng vai trò then chốt trong công tác phòng chống thiên tai với tổng chiều dài khoảng 3.400 km đê biển và 1.300 km đê cửa sông. Tại khu vực miền Trung, lưu vực sông Cả có diện tích rộng 27.200 km2 với tổng lượng dòng chảy trung bình năm đạt từ 21 đến 24 tỷ m3, đổ ra biển Đông qua cửa thoát duy nhất là Cửa Hội. Vùng cửa sông ven biển này là đầu mối kinh tế trọng điểm nhưng thường xuyên phải đối mặt với nguy cơ thiên tai kép từ lũ thượng nguồn và tác động của biển cả như bão, triều cường với biên độ dao động lên đến 3,72 m và nước dâng do bão từng đạt mốc kỷ lục 2,5 m.

Thực tiễn xây dựng và quản lý đê điều cho thấy việc phân định ranh giới giữa đê sông và đê cửa sông vẫn còn nhiều bất cập, thiếu cơ sở khoa học đồng bộ. Đê sông chịu tải trọng chủ yếu từ mực nước lũ thượng nguồn, trong khi đê cửa sông phải chịu tổ hợp tải trọng thủy động lực học phức tạp gồm lũ, sóng biển, dòng triều và nước dâng. Sự thiếu rõ ràng trong phân định dẫn tới nguy cơ thiết kế sai tiêu chuẩn, gây lãng phí nguồn vốn đầu tư hoặc mất an toàn công trình.

Mục tiêu của nghiên cứu là xây dựng bộ tiêu chí khoa học và thiết lập phương pháp mô hình hóa thủy động lực nhằm xác định chính xác ranh giới đê sông và đê cửa sông tại khu vực Cửa Hội, sông Cả, tỉnh Nghệ An. Phạm vi không gian nghiên cứu bao quát mạng lưới hạ lưu sông Cả, sông La, sông Lam và mở rộng ra vùng biển ngoài Cửa Hội với cự ly 70 km ra khơi và 40 km dọc bờ biển. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc tối ưu hóa chi phí đầu tư cho các dự án nâng cấp đê trọng điểm, bảo vệ an toàn cho hơn 1.000.000 người dân và trên 45.600 ha đất canh tác ven sông Lam.

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 động lực học sông biển và cơ học chất lỏng chuyên sâu để phân tích sự tương tác giữa dòng chảy lũ và sóng triều. Lý thuyết truyền sóng triều vào sông được xây dựng trên hệ phương trình liên tục và phương trình chuyển động Saint-Venant kết hợp định luật Chézy. Thông qua phép tuyến tính hóa lực ma sát đáy, nghiệm dao động của sóng triều tắt dần khi truyền ngược về phía thượng lưu được mô tả qua hàm số mũ suy giảm $e^{\lambda x}$, trong đó hệ số tắt dần $\lambda$ đạt giá trị khoảng $3 \times 10^{-5},\text{m}^{-1}$, chứng minh rằng sau quãng đường 100 km biên độ sóng triều giảm tới 95% và chỉ còn lại 5%.

Khung lý thuyết kết hợp ba khái niệm nền tảng:

  • Cửa sông dạng phễu (estuary): Khu vực có chế độ thủy triều mạnh khống chế hình thái dòng chảy với vận tốc dòng triều đạt từ 20 đến 180 cm/s.
  • Ngưỡng độ lớn triều cho phép $[a]$: Tiêu chí giới hạn thủy lực với giá trị $[a] \le 0,5,\text{m}$ để xác định điểm kết thúc ảnh hưởng rõ rệt của động lực biển vào mùa lũ.
  • Tổ hợp tải trọng thủy động lực học tương tác: Trạng thái kết hợp bất lợi nhất giữa đường quá trình lưu lượng lũ thiết kế tần suất 1% và mực nước triều dâng cực trị.

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

Nguồn dữ liệu phục vụ nghiên cứu được thu thập toàn diện từ 9 trạm khí tượng và 27 trạm thủy văn trong lưu vực sông Cả với chuỗi số liệu quan trắc liên tục trên 40 năm (giai đoạn 1956–2004), kết hợp chuỗi dữ liệu hoa sóng và gió bão 12 năm (1997–2009) tại Cửa Hội và đảo Hòn Ngư. Cỡ mẫu nghiên cứu bao gồm toàn bộ các trận bão đổ bộ khu vực Thanh - Nghệ - Tĩnh với tần suất trung bình 1,7 cơn bão/năm và các trận lũ lịch sử điển hình.

Phương pháp chọn mẫu tập trung vào hai trận lũ lớn thực đo để phục vụ hiệu chỉnh và kiểm định mô hình số: trận lũ từ ngày 28/10 đến 12/11/2008 và trận lũ từ ngày 27/09 đến 07/10/2009 tại trạm thủy văn Chợ Tràng và Cửa Hội. Lý do lựa chọn bộ công cụ mô hình toán MIKE của Viện Thủy lực Đan Mạch (DHI) là khả năng tích hợp vượt trội giữa mô hình thủy lực một chiều MIKE 11 HD (sử dụng sơ đồ sai phân hữu hạn 6 điểm ẩn Abbott-Ionescu cho mạng sông) và mô hình thủy động lực hai chiều MIKE 21FM HD (sử dụng lưới tam giác linh hoạt cho vùng cửa biển 70 km x 40 km). Phương pháp này cho phép mô phỏng chính xác sự biến dạng phi tuyến của sóng triều và nước dâng khi gặp dòng lũ thượng nguồn có lưu lượng đỉnh vượt 5.200 m3/s.

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 thủy động lực và phân tích số liệu thực đo đã chỉ ra các phát hiện khoa học mang tính bước ngoặt:

  • Đặc trưng thủy triều và nước dâng cực trị: Tại Cửa Hội, chế độ triều là dạng triều pha tạp với biên độ lớn nhất đạt 3,72 m, trong khi tại Chợ Tràng biên độ triều mùa cạn đạt 2,84 m. Mực nước dâng do bão lớn nhất ghi nhận được là 2,5 m trong cơn bão số 7 năm 1989; khi kết hợp với đỉnh triều 1,0 m đã tạo ra mực nước tổng cộng 3,50 m, trực tiếp phá hủy tuyến đê Hội Thống.
  • Quy luật diễn biến hình thái cửa sông dạng phễu: Cửa Hội tồn tại lưỡi cát bồi tụ phát triển mạnh theo hướng Đông Nam dưới tác động của sóng hướng Đông Nam và Đông Bắc. Lưỡi cát kéo dài tạo dạng móc câu làm thu hẹp mặt cắt thoát lũ, gây xói lở nghiêm trọng 12 km bờ phía Bắc thuộc xã Nghi Hải (bờ biển lấn sâu hàng trăm mét trong 15 năm) và bồi lắng mạnh ở bờ Nam thuộc xã Xuân Hội.
  • Dịch chuyển ranh giới thủy động lực: Phân tích đường bao mực nước lớn nhất giữa hai tổ hợp kịch bản lũ và triều cho thấy biên độ dao động mực nước do triều giảm mạnh từ 3,22 m tại Cửa Hội xuống dưới 0,5 m khi tiến sâu vào nội địa. Ranh giới đê sông và đê cửa sông trên dòng sông Lam được xác định chính xác tại vị trí cách Cửa Hội khoảng 15 km về phía thượng lưu (khu vực nằm giữa Bến Thủy và Chợ Tràng), nơi mức độ chi phối của sóng và triều biển giảm hơn 80% so với ngoài cửa sông.

Thảo luận kết quả

Cơ chế thủy động lực giải thích hiện tượng trên là do hình thái lòng sông mở rộng dạng phễu ở hạ lưu làm khuếch đại sóng triều ban đầu, nhưng khi dòng triều đi sâu qua Bến Thủy, lực cản ma sát đáy sông kết hợp cùng áp lực thủy tĩnh của dòng lũ thượng nguồn đổ về với lưu lượng trên 5.000 m3/s đã triệt tiêu hoàn toàn năng lượng sóng biển.

Trong báo cáo khoa học, toàn bộ kết quả mô phỏng được thể hiện chi tiết thông qua các biểu đồ trắc dọc mực nước dọc tuyến sông Lam, bảng thống kê sai số mực nước đỉnh lũ hiệu chỉnh mô hình tại trạm Chợ Tràng với độ lệch dưới 0,15 m, và bản đồ vector trường vận tốc dòng chảy hai chiều tại Cửa Hội trong cả hai pha triều lên và triều rút. So với các hướng tiếp cận thống kê truyền thống vốn đòi hỏi chuỗi số liệu đo liên tục tốn kém, phương pháp mô hình hóa kết hợp 1D-2D đã khắc phục hoàn toàn sự thiếu hụt số liệu quan trắc, cung cấp cơ sở phân tích ứng suất dòng chảy và tải trọng sóng chân thực hơn 30% so với phương pháp tính toán giải tích đơn thuần.

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

Dựa trên kết quả nghiên cứu khoa học, 4 nhóm giải pháp chiến lược được đề xuất nhằm nâng cao hiệu quả quy hoạch và bảo vệ công trình đê điều:

  • Chuẩn hóa bộ tiêu chuẩn phân định đê điều cấp quốc gia: Bộ Nông nghiệp và Phát triển nông thôn cùng Tổng cục Thủy lợi cần ban hành văn bản quy phạm kỹ thuật, áp dụng tiêu chí độ lớn triều cho phép $[a] \le 0,5,\text{m}$ kết hợp mô hình số thủy động lực để rà soát, phân loại lại toàn bộ 1.300 km đê cửa sông trên toàn quốc trước năm 2026.
  • Đầu tư kiên cố hóa các tuyến đê cửa sông xung yếu: Ủy ban nhân dân hai tỉnh Nghệ An và Hà Tĩnh cần ưu tiên bố trí nguồn vốn nâng cấp 10 km đê Hội Thống và hơn 7 km đê Hữu Lam theo tiêu chuẩn thiết kế đê cửa sông chịu tải trọng sóng biển và nước dâng 2,5 m, mục tiêu hoàn thành gia cố mặt đê và kè mỏ hàn chữ T trong giai đoạn 2026–2028.
  • Xây dựng mạng lưới quan trắc thủy văn tự động vùng cửa sông: Cục Khí tượng Thủy văn cần bổ sung 4 trạm đo mực nước và độ mặn tự động theo thời gian thực tại Cửa Hội, Bến Thủy, Chợ Tràng và Linh Cảm trước năm 2027 nhằm nâng cao độ chính xác dự báo ngập lụt lên trên 90%.
  • Định kỳ nạo vét luồng tiêu thoát lũ Cửa Hội: Ban Quản lý Dự án Đầu tư Xây dựng Công trình Thủy lợi tỉnh Nghệ An chủ trì thực hiện nạo vét định kỳ lưỡi cát chắn cửa sông theo chu kỳ 2 năm một lần, đảm bảo duy trì mặt cắt thoát lũ sông Cả với lưu lượng thiết kế tối thiểu 5.200 m3/s.

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

Công trình nghiên cứu cung cấp giá trị học thuật và thực tiễn chuyên sâu cho 4 nhóm đối tượng:

  • Kỹ sư tư vấn thiết kế công trình thủy lợi, thủy điện và công trình biển: Nắm bắt phương pháp tính toán tổ hợp tải trọng thủy động lực học, áp dụng bộ phần mềm MIKE 11 HD và MIKE 21FM HD để tối ưu hóa kết cấu mặt cắt đê, giảm thiểu 15–20% khối lượng vật liệu xây dựng.
  • Cán bộ quản lý quy hoạch đê điều tại các Chi cục Thủy lợi địa phương: Sử dụng cơ sở khoa học để phân định chính xác ranh giới đê sông và đê cửa sông, từ đó xây dựng kế hoạch phân bổ ngân sách duy tu cho 104,2 km tuyến đê Tả Lam và 19,2 km đê La Giang.
  • Giảng viên, nhà nghiên cứu và học viên cao học chuyên ngành Kỹ thuật công trình biển và Thủy văn: Sử dụng luận văn như tài liệu tham khảo mẫu mực về quy trình thiết lập, hiệu chỉnh và kiểm định mô hình liên kết thủy lực một chiều và hai chiều cho các vùng cửa sông ven biển phức tạp.
  • Ban Chỉ huy Phòng chống thiên tai và Tìm kiếm cứu nạn cấp tỉnh: Khai thác kịch bản mô phỏng ngập lụt để xây dựng phương án sơ tán dân cư khẩn cấp khi xảy ra tổ hợp bão cấp 12 kết hợp triều cường cao 3,5 m tại vùng hạ du sông Cả.

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

  • Tiêu chí cốt lõi nào được sử dụng để phân định ranh giới giữa đê sông và đê cửa sông? Tiêu chí cốt lõi dựa trên tương tác thủy động lực học giữa sông và biển, cụ thể là độ lớn triều cho phép $[a] \le 0,5,\text{m}$ trên đường bao mực nước lũ kết hợp triều. Điểm ranh giới là vị trí mà tại đó tác động của sóng triều và nước dâng giảm xuống dưới ngưỡng quy định, nhường chỗ cho dòng chảy lũ một chiều từ sông chi phối hoàn toàn.

  • Tại sao không thể áp dụng chung một tiêu chuẩn thiết kế cho đê sông và đê cửa sông? Đê sông chủ yếu chịu tải trọng áp lực nước tĩnh và vận tốc dòng chảy lũ xuôi dòng. Ngược lại, đê cửa sông phải chịu đồng thời tải trọng sóng vỗ, sóng tràn, triều cường và nước dâng do bão lên đến 2,5 m. Việc áp dụng tiêu chuẩn đê sông cho vùng cửa sông sẽ gây mất ổn định mái đê, dẫn đến vỡ đê như sự cố đê Hội Thống năm 1989.

  • Mô hình MIKE 11 HD và MIKE 21FM HD đóng vai trò gì trong đề tài này? MIKE 11 HD giải phương trình thủy lực một chiều cho toàn bộ mạng lưới nhánh sông thượng nguồn sông Cả, sông La, sông Hiếu. MIKE 21FM HD mô phỏng trường dòng chảy và sóng hai chiều ngoài vịnh và cửa biển Cửa Hội. Việc liên kết hai mô hình giúp phản ánh trung thực bài toán tương tác thủy động lực học phi tuyến tại vùng giáp ranh.

  • Lưỡi cát ngầm tại Cửa Hội ảnh hưởng như thế nào đến chế độ thoát lũ? Lưỡi cát phát triển bồi tụ dài hàng trăm mét theo hướng Đông Nam chắn ngang cửa sông làm co hẹp lòng thoát lũ. Hiện tượng này làm tăng mức ngập lụt cục bộ vùng hạ du từ 0,3 đến 0,5 m trong các trận lũ lớn, đồng thời làm lệch dòng chảy gây xói lở 12 km đê kè ven bờ biển Nghi Hải.

  • Phương pháp nghiên cứu trong luận văn có thể mở rộng cho các cửa sông khác không? Hoàn toàn có thể mở rộng áp dụng cho tất cả các hệ thống cửa sông dạng phễu hoặc dạng doi cát tại miền Trung và đồng bằng sông Hồng. Quy trình thiết lập mô hình kết hợp 1D-2D và tiêu chí phân định $[a] \le 0,5,\text{m}$ cung cấp một khung phương pháp chuẩn mực, có thể chuyển giao trực tiếp cho các lưu vực sông Mã, sông Gianh hay sông Thu Bồn.

Kết luận

  • Đề tài đã xây dựng thành công cơ sở khoa học và hệ thống tiêu chí định lượng chuẩn xác nhằm phân định ranh giới giữa đê sông và đê cửa sông dựa trên quan điểm tương tác thủy động lực học sông biển.
  • Ứng dụng thành công chuỗi mô hình liên kết MIKE 11 HD và MIKE 21FM HD được hiệu chuẩn từ chuỗi số liệu thực đo trên 40 năm để mô phỏng chính xác các kịch bản lũ và triều dâng cực trị.
  • Xác định cụ thể vị trí ranh giới đê sông - đê cửa sông trên sông Lam cách Cửa Hội 15 km về phía thượng nguồn, làm căn cứ điều chỉnh cấp kỹ thuật cho hơn 170 km đê điều tại Nghệ An và Hà Tĩnh.
  • Đề xuất lộ trình 4 nhóm giải pháp kỹ thuật, hoàn thiện quy chuẩn quốc gia và nâng cấp các tuyến đê xung yếu với mục tiêu hoàn thành đồng bộ trong giai đoạn 2026–2028.
  • Cung cấp nguồn tài liệu học thuật và thực tiễn giá trị cao cho các cơ quan quản lý đê điều, đơn vị tư vấn thiết kế và các nhà khoa học ngành kỹ thuật công trình biển.

Hãy tham khảo và ứng dụng ngay các phương pháp mô hình hóa thủy động lực học tiên tiến này vào quy hoạch đê điều tại địa phương để nâng cao năng lực chống chịu thiên tai và bảo vệ an toàn bền vững cho các vùng đồng bằng ven biển.