Tổng quan nghiên cứu

Khu vực đồng bằng châu thổ sông Cửu Long hàng năm vận chuyển khoảng 140 triệu tấn bùn cát, trong đó 80% tập trung vào mùa lũ. Tuy nhiên, sự gia tăng nhanh chóng của nhu cầu xây dựng đã biến cát lòng sông thành nguồn tài nguyên bị khai thác quá mức. Đề tài nghiên cứu tập trung giải quyết bài toán mâu thuẫn giữa nhu cầu khai thác khoáng sản phục vụ phát triển hạ tầng và hiểm họa sạt lở bờ sông đe dọa sinh kế dân cư tại tỉnh Vĩnh Long.

Mục tiêu cốt lõi của nghiên cứu là đánh giá định lượng tác động của hoạt động khai thác cát đến sự biến đổi hình thái đáy sông và diễn biến sạt lở bờ, từ đó đề xuất hệ thống giải pháp phòng chống thích hợp. Phạm vi không gian nghiên cứu trải dài trên tổng chiều dài 24 km thuộc đoạn sông Tiền (nhánh trái 12 km) và sông Cổ Chiên (nhánh phải 12 km), đi qua địa bàn thành phố Vĩnh Long và huyện Long Hồ, giới hạn trong tọa độ VN-2000 từ vĩ độ Bắc 1132000 đến 1144000 và kinh độ Đông từ 542000 đến 558000. Phạm vi thời gian đánh giá dựa trên chuỗi số liệu đo đạc thực tế giai đoạn 2009 - 2012 và xây dựng sơ đồ dự báo sạt lở đến năm 2020.

Nghiên cứu mang ý nghĩa thực tiễn lớn khi cung cấp cơ sở khoa học để quy hoạch 4,5 triệu m³ trữ lượng cát được cấp phép trên địa bàn tỉnh, đồng thời bảo vệ hệ thống hạ tầng đô thị cho khu vực có mật độ dân số cao lên tới 698 người/km². Kết quả tính toán giúp kiểm soát tốc độ biến dạng lòng dẫn, hạn chế nguy cơ sạt lở đường bờ vốn đã ghi nhận mức thoái lui nghiêm trọng từ 10m đến 30m tại các điểm xung yếu.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn xây dựng trên sự kết hợp liên ngành giữa cơ học đất, thủy lực học sông ngòi và động lực học trầm tích thông qua 3 khung lý thuyết nền tảng:

  • Lý thuyết dòng chảy hai chiều (2D) vùng nước nông: Ứng dụng hệ phương trình Navier-Stokes tích phân theo phương đứng nhằm mô phỏng chính xác trường lưu tốc và dao động mực nước dưới tác động của chế độ bán nhật triều không đều.
  • Lý thuyết vận chuyển bùn cát và biến đổi đáy sông: Dựa trên phương trình cân bằng bùn cát liên tục, xác định lưu tốc khởi động và lưu lượng bùn cát đáy bở rời dưới tác động của dòng chảy kết hợp sóng.
  • Lý thuyết cân bằng giới hạn mái dốc Mohr-Coulomb: Phân tích độ ổn định của khối trượt bờ sông khi chân mái dốc bị xói sâu hoặc bị nạo vét mất chân.

Năm khái niệm chuyên ngành then chốt được lượng hóa gồm:

  • Xâm thực ngang: Quá trình dòng nước bào mòn làm dịch chuyển đường bờ vào đất liền.
  • Xâm thực sâu: Hiện tượng đào khoét đáy sông tạo thành các hố xói sâu cục bộ.
  • Góc nghỉ tự nhiên của cát: Giá trị giới hạn ổn định của hạt cát, xác định ở trạng thái khô trung bình là 32,01° và trạng thái ngập nước tĩnh trung bình là 27,81°.
  • Vận tốc khởi động bùn cát: Vận tốc dòng chảy tối thiểu làm hạt trầm tích đáy (kích thước 0,1 - 0,25 mm) bắt đầu dịch chuyển.
  • Chế độ bán nhật triều: Chu kỳ dao động triều 2 lần nước lớn và 2 lần nước ròng mỗi ngày đêm, chi phối mạnh mẽ lưu tốc mùa kiệt từ 0,3 m/s đến 1,2 m/s.

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

Nguồn dữ liệu đầu vào bao gồm chuỗi quan trắc khí tượng thủy văn tại trạm Mỹ Thuận (lưu tốc cực đại mùa lũ đạt 2,45 m/s), bản đồ địa hình đáy sông tỷ lệ 1:50.000 lập năm 2009, ảnh viễn thám độ phân giải cao và hồ sơ địa kỹ thuật các lỗ khoan sâu đến 100m.

Phương pháp phân tích chủ đạo là mô hình hóa số trị bằng bộ phần mềm chuyên dụng MIKE 21 (gồm module MIKE 21 Flow Model FM tính thủy động lực và MIKE 21 ST mô phỏng dòng vận chuyển bùn cát), kết hợp công nghệ GIS MapInfo để chồng xếp bản đồ biến động lòng dẫn. Cỡ mẫu nghiên cứu bao trùm toàn bộ 8 mỏ cát được cấp phép với 19 phương tiện khai thác hoạt động trong năm 2011, phân tích cơ lý trên 4 tầng địa chất chủ yếu và thiết lập mô phỏng chi tiết trên 5 mặt cắt ngang thủy lực đại diện (từ I-I đến V-V).

Phương pháp chọn mẫu mặt cắt dựa trên nguyên tắc đại diện hình thái: mặt cắt qua đoạn uốn cong, đoạn thắt hẹp dòng chảy và khu vực tập trung mỏ khai thác lớn. Mô hình số lưới tam giác phi cấu trúc được lựa chọn vì khả năng mô phỏng biên bờ uốn lượn phức tạp và thích ứng với sự thay đổi độ sâu đáy sông từ 7m đến hơn 42m, mang lại độ chính xác vượt trội so với các mô hình 1D truyền thống. Quá trình nghiên cứu và xử lý mô phỏng được tác giả thực hiện trong 5 tháng (từ tháng 07/2012 đến tháng 11/2012).

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

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

Nghiên cứu chỉ ra 4 phát hiện quan trọng có tính định lượng về hình thái và động lực dòng sông:

  • Hạ thấp cao trình và biến dạng đáy sông cục bộ: Hoạt động nạo vét tại 8 khu vực mỏ làm biến đổi nghiêm trọng địa hình lòng dẫn. Thân cát lòng sông tuổi aQII33 có chiều dày trung bình 3,67m - 5,75m bị khai thác sâu, hình thành các hố trũng nhân tạo. Đáy sông tự nhiên tại các vùng thắt hẹp bị đào sâu tới -30m ở thượng lưu cầu Mỹ Thuận và sâu tới 42m tại hạ lưu rạch Cái Cối.

  • Mất ổn định mái dốc do vượt góc dốc giới hạn: Lớp đất bờ sông thuộc hệ tầng aQII32 có thành phần bùn sét dẻo chảy dày 10m - 15m với chỉ tiêu cơ lý rất yếu (lực dính c = 0,125 kg/cm², góc ma sát trong φ = 7°47'). Góc dốc ổn định cho phép theo tính toán lý thuyết chỉ đạt 20°, nhưng mái dốc bờ sông thực tế phổ biến từ 30° đến 50°, thậm chí vách đứng 70° - 80°, tạo nguy cơ trượt lở dạng cung tròn rất cao.

  • Tác động cộng hưởng của vận tốc dòng chảy và sóng tàu: Mùa lũ ghi nhận lưu tốc dòng chảy lên đến 2,45 m/s, vượt xa vận tốc khởi động xói của cát hạt mịn (chiếm 80,8% - 87,6% thành phần hạt đáy). Khi các mỏ hoạt động hết công suất, lưu tốc dòng chảy cục bộ tăng từ 12% đến 18%, kết hợp sóng do tàu thuyền tải trọng 2.000 - 3.500 tấn tạo áp lực ngang liên tục thúc đẩy xói chân bờ.

  • Tốc độ sạt lở đường bờ gia tăng nhanh: Tại các điểm nóng như đầu cồn An Bình, đường bờ đã bị xói lở và lùi sâu từ 10m đến 30m trong giai đoạn 2002 - 2009 trên chiều dài hơn 1.000m. Năm 2011, tổng khối lượng cấp phép là 1.400.000 m³ nhưng lượng khai thác thực tế báo cáo là 586.611 m³ (chiếm 41,9%), sự chênh lệch này cho thấy tình trạng khai thác không đều tạo ra sự mất cân bằng bùn cát nghiêm trọng.

Thảo luận kết quả

Nguyên nhân cốt lõi của sạt lở bờ sông bắt nguồn từ việc khai thác cát làm gia tăng độ dốc mái bờ dưới nước. Khi cát đáy bị lấy đi, gradient thủy lực tại sát bờ tăng lên, làm dòng chảy xoáy khoét sâu chân vách bờ, tạo hàm ếch khiến khối bùn sét dẻo chảy bên trên sụp đổ hoàn toàn.

So với các nghiên cứu động lực học sông ngòi trước đây chỉ tập trung vào dòng chảy mùa lũ, nghiên cứu này chứng minh rằng việc khai thác cát kết hợp dao động bán nhật triều mùa kiệt (biên độ triều chênh lệch mực nước từ -1m đến +1m) làm tăng tần suất xuất hiện dòng chảy ngược lưu tốc 1,2 m/s, liên tục làm suy yếu độ bền cắt của đất bờ dốc.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua chuỗi biểu đồ mặt cắt ngang địa hình đáy sông I-I đến V-V, cho phép so sánh chi tiết hình học lòng dẫn trước và sau khi khai thác. Bản đồ phân bố vector lưu tốc 2D và biểu đồ so sánh vận tốc khởi động hạt trầm tích minh chứng rõ ràng sự thay đổi dòng chảy tầng đáy, giúp các nhà quy hoạch dễ dàng nhận diện ranh giới xói lở.

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

Nhằm ổn định lòng dẫn và giảm thiểu thiệt hại sạt lở đến năm 2020, luận văn đề xuất 4 nhóm giải pháp đồng bộ:

  • Tái cơ cấu và siết chặt quy hoạch cấp phép khai thác: Cắt giảm từ 30% đến 50% công suất khai thác tại các khu vực mỏ nằm gần bờ lõm uốn khúc hoặc gần vùng sạt lở xung yếu như An Bình, Bình Hòa Phước; dừng cấp phép mới tại các nhánh sông có chiều rộng dưới 500m. Cơ quan thực hiện: Sở Tài nguyên và Môi trường tỉnh Vĩnh Long, thực hiện định kỳ hàng năm từ 2013.

  • Xây dựng công trình chỉnh trị và hướng dòng chủ động: Đầu tư xây dựng hệ thống kè mỏ hàn và kè hướng dòng tại các vị trí phân lưu sông Tiền - Cổ Chiên nhằm giảm 25% - 35% lưu tốc dòng áp bờ, đẩy dòng chủ lưu ra giữa lòng dẫn. Cơ quan thực hiện: Sở Nông nghiệp và Phát triển Nông thôn phối hợp Ban Quản lý Dự án Giao thông Thủy lợi, hoàn thành trong giai đoạn 2013 - 2018.

  • Gia cố bờ sông bằng giải pháp công trình thụ động: Thi công kè lát mái đá hộc kết hợp vải địa kỹ thuật hoặc tường cừ bê tông dự ứng lực tại các vách dốc đứng trên 35° thuộc Phường 1, Phường 5 và khu vực phà Đình Khao, nâng hệ số ổn định mái dốc K > 1,30. Cơ quan thực hiện: UBND thành phố Vĩnh Long và UBND huyện Long Hồ, triển khai giai đoạn 2014 - 2020.

  • Thiết lập cơ chế kiểm soát công nghệ cao và giới hạn tải trọng: Bắt buộc 100% phương tiện khai thác và vận tải cát phải gắn thiết bị giám sát hành trình GPS, đo sâu tự động; quy định tốc độ chạy tàu dưới 10 km/h tại các kênh hẹp để giảm 40% áp lực sóng va bờ. Cơ quan thực hiện: Công an Đường thủy phối hợp Sở Giao thông Vận tải tỉnh Vĩnh Long, áp dụng ngay từ quý I/2013.

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

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường (Sở TN&MT các tỉnh Đồng bằng sông Cửu Long): Ứng dụng số liệu tính toán trữ lượng thân cát (tổng tài nguyên cấp 333 đạt trên 28 triệu m³) và bản đồ dự báo xói lở để xây dựng quy hoạch khai thác khoáng sản an toàn, cấp phép khai thác bền vững.

  • Cơ quan quản lý phòng chống thiên tai và quy hoạch đô thị (Sở Nông nghiệp & PTNT, Viện Quy hoạch Thủy lợi): Sử dụng sơ đồ các điểm sạt lở nguy hiểm (đoạn 1.500m tại Bình Hòa Phước, 1.000m tại đầu cồn An Bình) để bố trí tái định cư dân cư ven sông và phân bổ vốn đầu tư bờ kè kiên cố.

  • Các kỹ sư tư vấn thiết kế công trình thủy lợi, giao thông thủy: Tham khảo các thông số địa kỹ thuật đất yếu (góc ma sát trong 7°47', lực dính 0,125 kg/cm²) và lưu tốc xói lũ 2,45 m/s để thiết kế kết cấu kè bảo vệ bờ sông tối ưu về mặt chi phí và kỹ thuật.

  • Giảng viên, nhà nghiên cứu và học viên cao học chuyên ngành Địa kỹ thuật, Kỹ thuật Địa chất: Tiếp cận phương pháp luận ứng dụng mô hình số trị MIKE 21 FM và ST kết hợp phân tích tương tác đất - nước phục vụ các đề tài nghiên cứu động lực học lòng dẫn châu thổ.

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

Hoạt động khai thác cát lòng sông gây sạt lở bờ theo cơ chế nào?

Khai thác cát làm hạ thấp cao trình đáy sông cục bộ từ 2m đến 4m, tạo hố trũng thu hút dòng chủ lưu. Quá trình này làm gia tăng vận tốc dòng chảy sát đáy, làm chân vách bờ dốc đứng vượt góc dốc ổn định 20°, gây hiện tượng xói ngầm hàm ếch và trượt đổ khối bùn sét bên trên.

Vì sao bờ sông tại khu vực nghiên cứu có khả năng tự chống xói lở rất kém?

Cấu trúc bờ sông được tạo bởi lớp sét bột pha cát tuổi aQII32 dày tới 15m ở trạng thái chảy và dẻo chảy. Đất có lực dính yếu 0,125 kg/cm² và góc ma sát trong chỉ 7°47', khi chịu tác động của lưu tốc dòng chảy mùa lũ lên tới 2,45 m/s sẽ nhanh chóng bị tan rã.

Mô hình MIKE 21 FM đóng góp gì nổi bật cho kết quả của luận văn?

Mô hình MIKE 21 FM thiết lập lưới tam giác phi cấu trúc phản ánh chính xác tương tác thủy lực 2 chiều và biến đổi bùn cát trên chiều dài 24 km sông. Mô hình cho phép mô phỏng định lượng sự thay đổi lưu tốc và xói bồi đáy sông giữa trường hợp có và không có hoạt động khai thác mỏ.

Áp lực sóng do tàu bè vận chuyển cát tác động ra sao đến bờ sông?

Tàu chở cát có tải trọng từ 2.000 đến 3.500 tấn khi lưu thông trên luồng hẹp tạo ra các đợt sóng va đập liên tục vào mái bờ. Áp lực động học này làm xói mòn lớp mùn thực vật bảo vệ mặt bãi bồi, đẩy nhanh tốc độ phá hủy liên kết đất tại mép nước ròng và nước lớn.

Đâu là giải pháp cấp bách nhất cần triển khai để bảo vệ bờ sông Vĩnh Long?

Giải pháp cấp bách nhất là khoanh vùng cấm khai thác tuyệt đối tại các hố xói sâu trên 30m và đầu cồn An Bình, đồng thời xây dựng công trình kè mỏ hàn kết hợp lát mái đá hộc tại các đoạn xung yếu nhằm giảm lưu tốc dòng áp bờ xuống dưới 1,0 m/s.

Kết luận

  • Đề tài lượng hóa thành công tác động của 8 khu vực mỏ khai thác cát đến quá trình biến đổi địa hình đáy và sạt lở bờ sông trên đoạn sông dài 24 km thuộc tỉnh Vĩnh Long.
  • Xác định nguyên nhân địa kỹ thuật chính: Mái dốc tự nhiên thực tế 30° - 80° vượt quá góc ổn định cho phép 20° của nền bùn sét dẻo chảy aQII32.
  • Ứng dụng hiệu quả bộ công cụ mô hình số MIKE 21 FM và ST để mô phỏng chính xác trường dòng chảy mùa lũ đạt lưu tốc 2,45 m/s và diễn biến vận chuyển bùn cát.
  • Xây dựng bản đồ phân vùng dự báo sạt lở bờ sông chi tiết đến năm 2020, làm căn cứ khoanh định các khu vực cấm và hạn chế khai thác khoáng sản.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, gắn liền với kế hoạch giám sát công nghệ GPS định kỳ giai đoạn 2013 - 2020.

Đóng góp khoa học lớn nhất của luận văn là việc kết hợp chặt chẽ giữa khảo sát địa chất công trình và mô hình hóa thủy động lực 2D, cung cấp phương pháp luận chuẩn xác cho công tác quản lý tài nguyên cát ven sông. Trong các giai đoạn tiếp theo, địa phương cần nhanh chóng thể chế hóa các đề xuất vào quy hoạch phát triển kinh tế - xã hội, cắm mốc cảnh báo sạt lở thực địa và lắp đặt hệ thống quan trắc tự động để bảo vệ an toàn tính mạng, tài sản của cộng đồng dân cư ven sông.