Tổng quan nghiên cứu

Khu vực Đồng bằng sông Cửu Long đang đối mặt với tình trạng biến hình lòng dẫn và xói lở bờ sông diễn biến hết sức phức tạp. Thống kê trên toàn tuyến sông Tiền ghi nhận 37 điểm xói lở nghiêm trọng, trong đó có 7 điểm xói lở mạnh với tốc độ vượt quá 10 m/năm, 18 điểm xói lở vừa từ 5 đến 10 m/năm và 12 điểm xói lở yếu dưới 5 m/năm. Điển hình tại huyện Hồng Ngự, tỉnh Đồng Tháp, đợt sạt lở lịch sử năm 1992 đã cướp đi sinh mạng của 10 người, làm bị thương 8 người, phá hủy 13 công sở và gây thiệt hại kinh tế trên 20 tỷ đồng. Gần đây hơn, vào tháng 2 năm 2010, đoạn bờ sông Tiền qua cù lao Long Khánh tiếp tục sạt lở nghiêm trọng với chiều dài gần 300 m và ăn sâu vào đất liền khoảng 30 m, đe dọa trực tiếp đến tính mạng cùng sinh kế của hàng trăm hộ dân.

Vấn đề cốt lõi của hiện tượng xói lở tại khu vực phân lưu cù lao Long Khánh - Hồng Ngự bắt nguồn từ sự tương tác phức tạp giữa dòng chảy lũ thượng nguồn và hình thái lòng sông cong gấp. Vào mùa lũ, lưu tốc dòng chảy đạt tới 1,6 m/s, duy trì liên tục từ 3 đến 4 tháng, vượt xa ngưỡng vận tốc xói cho phép của nền đất bờ là 0,6 m/s. Đề tài tập trung giải quyết mục tiêu đánh giá thực trạng thủy lực, mô phỏng quá trình vận chuyển bùn cát, dự báo diễn biến xói bồi lòng dẫn trong giai đoạn 2009 - 2010 và đề xuất giải pháp công trình kè mỏ hàn nhằm chuyển hướng dòng chủ lưu.

Phạm vi nghiên cứu trải dài khoảng 30 km dọc sông Tiền, từ biên giới Việt Nam - Campuchia qua thị trấn Tân Châu đến hạ lưu cù lao Long Khánh. Nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc bảo vệ hơn 2.000 m đường bờ xung yếu, giảm thiểu rủi ro phải di dời khẩn cấp cho hơn 50 hộ dân ven sông và cung cấp cơ sở khoa học định lượng cho công tác quy hoạch chỉnh trị sông bền vững tại tỉnh Đồng Tháp.

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 nền tảng thủy lực sông ngòi hiện đại kết hợp với cơ sở lý thuyết chuyên sâu của bộ mô hình số hai chiều MIKE 21C do Viện Thủy lực Đan Mạch phát triển. Khung lý thuyết cốt lõi bao gồm:

  • Lý thuyết thủy động lực học 2 chiều trung bình theo chiều sâu: Dựa trên hệ phương trình Saint-Venant mở rộng trên hệ tọa độ cong trực giao, tích hợp các thành phần gia tốc đối lưu, gradient áp suất, lực Coriolis, ma sát đáy Chézy và ứng suất rối Reynolds.
  • Mô hình dòng chảy vòng thứ cấp (Helical Flow): Giải thích sự xuất hiện của cấu trúc xoắn ốc 3 chiều tại các khúc sông uốn cong. Cường độ dòng xoắn tỷ lệ với vận tốc dòng chính và tỷ số độ sâu trên bán kính cong, tạo ra góc lệch ứng suất trượt đáy tác động trực tiếp vào chân mái dốc bờ sông.
  • Lý thuyết chuyển động bùn cát và biến đổi hình thái: Áp dụng tích hợp các công thức bùn cát đáy kinh nghiệm của Meyer-Peter & Müller (1947), Engelund-Hansen (1967) và Van Rijn (1984), kết hợp phương trình khuếch tán bùn cát lơ lửng trung bình theo độ sâu theo phân bố Rouse.

Các khái niệm chính được chuẩn hóa trong mô hình bao gồm: ứng suất cắt đáy tới hạn Shields, độ thô thủy lực hạt trầm tích, hệ số nhám Manning-Chézy và phương trình cân bằng thể tích hình thái lòng dẫn Exner.

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

Nghiên cứu kết hợp chặt chẽ giữa phương pháp đo đạc thực địa và phương pháp mô phỏng toán số tiên tiến:

  • Nguồn dữ liệu đầu vào: Thu thập dữ liệu địa hình lòng sông chi tiết từ 3 đợt khảo sát thực tế do Viện Khoa học Thủy lợi miền Nam thực hiện năm 2009 và Viện Kỹ thuật Biển thực hiện các năm 2008, 2012. Chuỗi số liệu thủy văn biên lỏng (lưu lượng thượng lưu và mực nước hạ lưu) cùng nồng độ bùn cát lơ lửng giai đoạn 2009 - 2011 được trích xuất đồng bộ từ mô hình thủy lực tổng thể MIKE 11 đoạn từ Tân Châu đến Mỹ Thuận.
  • Quy mô mẫu và thiết lập lưới tính: Miền tính toán được rời rạc hóa thành 23.453 ô lưới phi cấu trúc cong trực giao với 499 điểm nút theo phương dọc sông và 47 điểm nút theo phương ngang sông. Kích thước ô lưới dao động linh hoạt từ 40 đến 60 m theo chiều dọc và từ 15 đến 50 m theo chiều ngang, tập trung mật độ dày tại các vị trí uốn khúc và ngã ba phân lưu.
  • Phương pháp phân tích và hiệu chỉnh: Sử dụng phương pháp thể tích hữu hạn giải số các phương trình đạo hàm riêng. Quá trình hiệu chỉnh và kiểm định mô hình thủy lực (HD) và hình thái (SD) được thực hiện độc lập tại 3 mặt cắt ngang trắc diện chính (MC1, MC2, MC3) tương ứng với chuỗi thời gian thực đo các năm 2009 và 2010, đảm bảo sai số lưu lượng và cao trình đáy nằm trong giới hạn kiểm soát cho phép dưới 10%. Mô hình lưới cong trực giao được lựa chọn vì khả năng bám sát đường bờ tự nhiên và giảm thiểu đáng kể sai số số học so với lưới chữ nhật truyền thống.

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

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

Mô phỏng toán số trên mô hình MIKE 21C trong khoảng thời gian 24 tháng (2009 - 2010) đã làm sáng tỏ quy luật thủy động lực và biến hình lòng dẫn đoạn sông Tiền qua cù lao Long Khánh:

  • Bất đối xứng lưu lượng phân lưu: Nhánh Hồng Ngự (bờ tả) tiếp nhận tỷ lệ lưu lượng lớn hơn từ 55% đến 62% tổng lưu lượng sông Tiền trong mùa lũ, trong khi nhánh Long Khánh chỉ tiếp nhận khoảng 38% đến 45%. Điều này làm gia tăng mạnh mẽ áp lực thủy động học lên tuyến bờ phía thị xã Hồng Ngự.
  • Vận tốc dòng chảy vượt ngưỡng xói: Vận tốc dòng chính mùa lũ đạt cực đại từ 1,55 m/s đến 1,60 m/s tại các mặt cắt đỉnh cong, cao gấp 2,67 lần so với vận tốc tới hạn gây xói của hạt bùn cát đáy (0,60 m/s). Trong mùa kiệt, vận tốc dòng chảy duy trì trong khoảng 0,40 m/s đến 0,75 m/s.
  • Phát triển hố xói sâu dị thường: Xuất hiện hệ thống các hố xói sâu từ -30 m đến -40 m (cá biệt có điểm đạt -35 m tại vùng ngã ba phân lưu cù lao Long Khánh). Sự chênh lệch độ sâu giữa hố xói và bãi cạn lên tới 25 m tạo nên mái dốc bờ sông cực kỳ mất ổn định.
  • Biến đổi hình thái lòng dẫn sau 24 tháng: Quá trình xói sâu đáy sông diễn ra mạnh mẽ trong 12 tháng đầu tiên với chiều sâu xói cục bộ từ 1,2 m đến 2,5 m dọc theo lạch sâu. Ngược lại, khu vực hạ lưu bến phà Tân Châu và đuôi bãi cạn có xu thế bồi tụ với độ nâng cao trình đáy từ 0,5 m đến 1,8 m.

Thảo luận kết quả

Cơ chế gây sạt lở bờ sông nghiêm trọng tại khu vực cù lao Long Khánh được xác định là quá trình tổ hợp hai giai đoạn: xói sâu lòng dẫn do thủy lực học là tiền đề, và sạt trượt mái bờ do mất cân bằng cơ học đất là kết quả tất yếu. Dòng chảy uốn khúc sinh ra dòng chảy vòng thứ cấp (helical flow), làm lệch hướng ứng suất tiếp đáy và hướng trực tiếp tia dòng xoáy mang động năng lớn vào chân mái bờ tả. Khi chân bờ bị xói rỗng tạo hàm ếch, khối đất bên trên mất lực chống trượt và sụp đổ hàng loạt.

Khi đối chiếu với các nghiên cứu tương tự tại đoạn sông Tiền qua xã An Hiệp (nơi từng sạt lở ăn sâu 70 m vào năm 2011) và kè bảo vệ bờ Tân Châu, kết quả mô hình MIKE 21C phản ánh hoàn toàn tương thích về mặt động học. Tác động kết hợp của biên độ triều biển Đông (truyền sâu qua 6 cửa sông Tiền với tổng chiều rộng cửa khoảng 9 km, lớn hơn nhiều so với 5 km của 2 cửa sông Hậu) gây ra hiện tượng triều rút nhanh, tạo áp lực nước lỗ rỗng lớn trong thân bờ dốc, thúc đẩy nhanh quá trình trượt lở sau các đợt lũ rút.

Các dữ liệu nghiên cứu được biểu diễn trực quan thông qua chuỗi biểu đồ so sánh phân bố lưu tốc mặt cắt ngang giữa tính toán và thực tế tại các thời điểm triều lên - triều rút các năm 2009 - 2010, cùng các bản đồ bình đồ 2D thể hiện sự dịch chuyển đường bờ (banklines) và độ dịch chuyển lạch sâu qua các giai đoạn 3, 6, 12, 18 và 24 tháng.

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

Dựa trên kết quả mô phỏng thủy lực và kiểm chứng diễn biến lòng dẫn, nghiên cứu đề xuất hệ thống 4 giải pháp công trình và phi công trình đồng bộ:

  • Xây dựng hệ thống 2 kè mỏ hàn chỉnh trị dòng chảy: Bố trí 2 tuyến kè mỏ hàn trọng lực tại bờ phía thượng lưu cù lao Long Khánh. Kết quả mô phỏng lại trên MIKE 21C khẳng định giải pháp này giúp đẩy trục dòng chủ lưu dịch chuyển ra xa bờ từ 50 m đến 80 m, giảm lưu tốc ven bờ từ 1,60 m/s xuống dưới 0,55 m/s (thấp hơn ngưỡng xói an toàn 0,60 m/s) và tạo vùng nước quẩn gây bồi lắng bùn cát bảo vệ chân bờ trong vòng 12 đến 24 tháng. Chủ thể thực hiện: Ban Quản lý Dự án Nông nghiệp phối hợp Sở Nông nghiệp và Phát triển Nông thôn tỉnh Đồng Tháp.
  • Gia cố bảo vệ trực tiếp chân mái bờ: Triển khai thả thảm đá hộc bọc thép kết hợp thảm bê tông tự chèn linh hoạt tại các hố xói sâu có cao trình dưới -20 m nhằm giữ ổn định độ dốc bờ tự nhiên, ngăn chặn hiện tượng tạo hàm ếch đáy sông trước mùa lũ hàng năm. Chủ thể thực hiện: Chi cục Thủy lợi tỉnh Đồng Tháp.
  • Thiết lập mạng lưới trắc đạc và quan trắc cảnh báo sớm: Ứng dụng công nghệ đo sâu quét đa chùm (Multi-beam Echo Sounder) và thiết bị đo dòng ADCP định kỳ 6 tháng một lần (trước và sau mùa lũ) nhằm theo dõi sự dịch chuyển của lạch sâu và phát hiện sớm các hố xói mới phát sinh. Chủ thể thực hiện: Viện Kỹ thuật Biển và Trung tâm Dự báo Khí tượng Thủy văn khu vực Nam Bộ.
  • Hoàn thiện quy hoạch hành lang thoát lũ và tái định cư: Thiết lập ranh giới an toàn cách mép bờ sạt lở tối thiểu 30 m đến 50 m, tuyệt đối cấm xây dựng công trình kiên cố trên bãi bồi ven sông và hoàn thành di dời 100% các hộ dân trong vùng nguy cơ cao trong kế hoạch giai đoạn 2015 - 2020. Chủ thể thực hiện: Ủy ban Nhân dân huyện Hồng Ngự và Ủy ban Nhân dân tỉnh Đồng Tháp.

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

Nội dung và kết quả khoa học của luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng trọng tâm:

  • Cơ quan quản lý nhà nước và hoạch định chính sách: Cung cấp cơ sở dữ liệu khoa học định lượng cho Ủy ban Nhân dân tỉnh Đồng Tháp, Sở Nông nghiệp và Phát triển Nông thôn, Ban Chỉ huy Phòng chống Thiên tai trong việc lập quy hoạch chỉnh trị sông Tiền và phê duyệt kế hoạch vốn đầu tư các dự án bảo vệ bờ cấp bách.
  • Kỹ sư công trình thủy và đơn vị tư vấn thiết kế: Cung cấp tài liệu tham khảo thực chứng về quy trình số hóa công trình kè mỏ hàn, thông số tính toán trường lưu tốc, độ nhám đáy Chézy và cách bố trí hình học mỏ hàn tối ưu trong điều kiện sông vùng đồng bằng chịu ảnh hưởng bán nhật triều.
  • Các viện nghiên cứu và giảng viên chuyên ngành kỹ thuật tài nguyên nước: Là nguồn tư liệu học thuật giá trị phục vụ công tác giảng dạy, nghiên cứu chuyên sâu về động lực học sông ngòi, cơ chế tương tác bùn cát và phương pháp thiết lập lưới cong trực giao phi cấu trúc trong bộ phần mềm MIKE 21C.
  • Học viên cao học và sinh viên ngành Xây dựng Công trình Thủy: Hướng dẫn chi tiết phương pháp tiếp cận từ thu thập số liệu thủy văn, xử lý điều kiện biên MIKE 11 - MIKE 21C đến kỹ năng phân tích, kiểm định mô hình hình thái hai chiều đáp ứng chuẩn mực học thuật của Đại học Bách Khoa.

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

Mô hình MIKE 21C có ưu điểm gì vượt trội so với mô hình lưới vuông MIKE 21 trong nghiên cứu này? Mô hình MIKE 21C sử dụng hệ tọa độ cong trực giao với 23.453 ô lưới linh hoạt, cho phép các đường lưới uốn lượn bám sát chính xác đường bờ tự nhiên và hướng dòng chảy cong của sông Tiền. Điều này giúp loại bỏ hoàn toàn hiện tượng răng cưa biên và mô phỏng chính xác dòng chảy xoắn ốc 3 chiều mà mô hình lưới vuông thông thường không thể xử lý tốt.

Tại sao dòng chảy mùa lũ lại là nguyên nhân trực tiếp gây xói lở nghiêm trọng tại cù lao Long Khánh? Vào mùa lũ, lưu tốc dòng chảy đạt tới 1,60 m/s và kéo dài liên tục từ 3 đến 4 tháng, vượt xa ngưỡng xói cho phép của đất bờ là 0,60 m/s. Động năng dòng nước cực lớn ép sát bờ uốn cong, tạo ra các dòng xoáy đáy làm xói sâu chân bờ tới độ sâu -35 m, dẫn đến sụp đổ bờ hàng loạt.

Giải pháp kè mỏ hàn giúp bảo vệ bờ sông hoạt động theo nguyên lý nào? Kè mỏ hàn được bố trí nhô ra lòng dẫn để hướng dòng chủ lưu có vận tốc cao trên 1,50 m/s ra xa bờ từ 50 m đến 80 m. Tại khu vực giữa hai mỏ hàn và sau thân kè, vận tốc dòng chảy giảm mạnh xuống dưới 0,55 m/s, tạo điều kiện cho bùn cát lắng đọng tự nhiên, bồi đắp và gia cố chân bờ sông.

Độ tin cậy của kết quả mô phỏng trong luận văn được bảo đảm bằng phương pháp nào? Mô hình được hiệu chỉnh và kiểm định độc lập thông qua chuỗi số liệu thực tế đo đạc liên tục 2 năm (2009 - 2010). Kết quả tính toán lưu lượng và phân bố vận tốc được so sánh chi tiết tại 3 mặt cắt kiểm định chính (MC1, MC2, MC3) và trắc diện dọc lòng dẫn, cho thấy mức độ tương đồng cao với số liệu khảo sát thực tế của Viện Kỹ thuật Biển.

Nghiên cứu có xét đến tác động của chế độ thủy triều lên diễn biến xói bồi không? Nghiên cứu đã tích hợp toàn diện ảnh hưởng của chế độ bán nhật triều biển Đông thông qua chuỗi số liệu mực nước biên hạ lưu. Tác động của triều rút nhanh kết hợp với lũ thượng nguồn tạo ra gradient áp suất thấm cao trong bờ đất, là một trong những nguyên nhân chính kích hoạt các đợt sạt lở quy mô lớn.

Kết luận

  • Xác định chính xác cơ chế thủy động lực học: Luận văn đã làm rõ nguyên nhân gây sạt lở bờ sông Tiền đoạn qua cù lao Long Khánh do sự kết hợp giữa dòng chảy lũ vận tốc cao 1,60 m/s, cấu trúc dòng chảy xoắn đỉnh cong và sự xuất hiện của các hố xói sâu tới -35 m.
  • Ứng dụng thành công công nghệ mô hình hóa tiên tiến: Thiết lập thành công bộ công cụ mô phỏng 2 chiều MIKE 21C trên mạng lưới 23.453 ô cong trực giao, kiểm định chuẩn xác diễn biến thủy lực và biến hình lòng sông trong giai đoạn 2009 - 2010.
  • Đề xuất giải pháp công trình kè mỏ hàn hiệu quả: Chứng minh bằng số liệu tính toán định lượng rằng bố trí 2 tuyến kè mỏ hàn sẽ giảm lưu tốc ven bờ xuống dưới 0,55 m/s, đẩy dòng chủ lưu ra xa bờ 50 đến 80 m và thúc đẩy quá trình bồi tụ bảo vệ bờ tự nhiên.
  • Lộ trình triển khai khuyến nghị: Đề xuất kế hoạch hoàn thiện thiết kế kỹ thuật công trình bảo vệ bờ trong giai đoạn 12 đến 24 tháng, kết hợp thiết lập hành lang an toàn 30 m và quan trắc định kỳ 6 tháng một lần trước mùa mưa lũ.
  • Đóng góp học thuật và giá trị thực tiễn: Công trình là tài liệu khoa học mẫu mực về phương pháp chỉnh trị sông ngòi vùng đồng bằng phân lưu, kêu gọi các cấp quản lý và nhà nghiên cứu tiếp tục ứng dụng mô hình số để chủ động ứng phó biến đổi lòng dẫn trên toàn lưu vực sông Cửu Long.