Tổng quan nghiên cứu

Khu vực Đồng bằng sông Cửu Long đang đối mặt với nguy cơ sạt lở bờ sông nghiêm trọng do áp lực đô thị hóa và hoạt động khai thác khoáng sản quá mức. Toàn vùng hiện có 126 doanh nghiệp được cấp phép khai thác cát, với khối lượng trầm tích bị lấy đi hàng năm lên tới 34 triệu m³, tương đương 55 triệu tấn, trong đó cát chiếm đến 90%. Hệ quả tiêu cực là sự hình thành của 265 điểm sạt lở bờ sông và bờ biển với tổng chiều dài vượt quá 450 km. Mỗi năm, khu vực này tổn thất khoảng 500 ha đất tự nhiên với tốc độ xói lở trung bình 30 m đến 40 m/năm. Dự báo đến năm 2050, xói lở ven bờ sẽ đe dọa sinh kế của khoảng 1 triệu người dân sinh sống tại lưu vực sông Tiền và sông Hậu.

Các đợt sạt lở điển hình đã ghi nhận thiệt hại nặng nề: tại xã Thường Phước, huyện Hồng Ngự xuất hiện 7 đoạn sạt lở nghiêm trọng; tại xã An Bình, huyện Long Hồ, sạt lở làm trôi 8.000 m² đất và nhấn chìm 23 bè cá gây thiệt hại gần 4 tỷ đồng; tại quận Cái Răng, vụ sạt lở dài 30 m đã ăn sâu vào bờ 5 m đến 7 m. Vấn đề cốt lõi đặt ra là việc thiếu các cơ sở định lượng về giới hạn nạo vét an toàn. Mục tiêu của nghiên cứu là xác lập tương quan toán học giữa khối lượng cát khai thác với hệ số an toàn ổn định bờ sông (Fs), xác định chiều sâu cực hạn và góc dốc giới hạn tại ba địa bàn trọng điểm: Vĩnh Long, An Giang và Cần Thơ. Đề tài cung cấp công cụ kỹ thuật giúp kiểm soát Fs đạt ngưỡng an toàn tối thiểu 1.4 theo quy chuẩn kỹ thuật quốc gia.

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 xây dựng trên nền tảng lý thuyết cân bằng giới hạn (Limit Equilibrium Method - LEM) và phương pháp phần tử hữu hạn (Finite Element Method - FEM) để đánh giá trạng thái ứng suất - biến dạng của mái dốc đất sét mềm bão hòa. Cơ chế phá hoại chủ đạo được xác định là dạng trượt cung tròn sâu tuân theo tiêu chuẩn bền Mohr-Coulomb.

Phương pháp Spencer (1967) được chọn làm mô hình tính toán trọng tâm trong nhóm cân bằng giới hạn vì thỏa mãn đồng thời điều kiện cân bằng lực theo phương ngang, phương thẳng đứng và cân bằng mô men quanh tâm trượt. Khung phân tích áp dụng nguyên lý ứng suất tổng cộng (Total Stress Analysis) cho nền đất sét bão hòa nước trên 95% với góc ma sát trong bằng 0 và lực dính bằng sức kháng cắt không thoát nước Su. Các khái niệm then chốt gồm: hệ số an toàn mái dốc Fs, sức kháng cắt không thoát nước Su, chiều sâu cực hạn Hcr, góc dốc cực hạn αcr và hệ số áp lực ngang tĩnh K0 tính theo hệ thức Lee & Jin (1979) với K0 = 0.31log(PI).

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

Dữ liệu đầu vào thu thập từ hồ sơ địa chất công trình gồm 14 mặt cắt khảo sát tại sông Cổ Chiên (Vĩnh Long) và 11 mặt cắt tại Trà Vinh, kết hợp chuỗi số liệu thủy văn giai đoạn 2013-2014 từ các trạm quan trắc Mỹ Thuận, Long Xuyên và Cần Thơ. Phương pháp chọn mẫu phân tầng có chủ đích được áp dụng nhằm tập trung vào các tầng bùn sét yếu có độ sâu từ 0 m đến -30 m, độ ẩm tự nhiên 59.3% và chỉ số dẻo PI trung bình 20.4%.

Phân tích ổn định mái dốc được thực hiện thông qua bộ công cụ mô phỏng GeoStudio (mô đun SLOPE/W) và Plaxis. Phương pháp Spencer được lựa chọn nhờ khả năng cho giá trị Fs thiên về an toàn và thuật toán ổn định. Để xác định sức kháng cắt Su, nghiên cứu áp dụng mô hình Wroth & Houlsby (1985) dựa trên chỉ số dẻo PI (Su = 0.00435 × PI × σz). Phương pháp này đạt hệ số tương quan R² xấp xỉ 0.944 và hệ số biến thiên COV thấp hơn đáng kể so với thí nghiệm cắt đất trực tiếp, cho phép chuẩn hóa dữ liệu khi so sánh với thí nghiệm cắt cánh hiện trường có hệ số hiệu chỉnh 0.85 theo tiêu chuẩn ASTM D2573.

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

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

Phân tích mô phỏng địa hình lòng sông biến đổi cho thấy việc hạ sâu lòng sông theo hướng tăng góc dốc mái taluy nguy hiểm hơn nhiều so với hạ vuông góc đáy sông. Khi mái dốc bị xói dốc đứng, hệ số an toàn Fs ban đầu giảm từ 3.95 xuống còn 2.292, tương đương mức suy giảm độ ổn định 41.9%, trong khi hình dạng cung trượt nguy hiểm mở rộng sâu vào bờ.

Biến động thủy triều tạo ra rủi ro mất ổn định nghiêm trọng. Tại mức nước kiệt từ -1.5 m đến -2.0 m, hệ số Fs của bờ sông tại Vĩnh Long và Cần Thơ giảm 18.5% đến 25.3% so với điều kiện mực nước trung bình do mất áp lực nước tĩnh trợ lực cho chân mái dốc.

Tải trọng công trình ven sông (nhà cấp 4 tương đương tải trọng rải đều 10 kPa đến 15 kPa) làm giảm chiều sâu nạo vét cực hạn Hcr từ 1.2 m đến 2.8 m trên các mặt cắt có bề rộng mái dốc L từ 40 m đến 90 m. Mối quan hệ giữa khối lượng cát khai thác và hệ số an toàn Fs thể hiện xu hướng suy giảm phi tuyến rõ rệt, cung cấp ngưỡng giới hạn thể tích nạo vét trước khi Fs rơi xuống dưới 1.4.

Thảo luận kết quả

Cơ chế trượt lở thực chất bắt nguồn từ hiện tượng mất cân bằng mô men khi khối đất chân mái dốc bị lấy đi. Việc khai thác cát làm gia tăng mô men gây trượt do trọng lượng bản thân khối trượt tăng khoảng cách cánh tay đòn tới tâm xoay, đồng thời làm suy giảm diện tích tiếp xúc kháng cắt dọc theo cung trượt tiềm năng.

       Mối quan hệ biến thiên giữa Thể tích nạo vét (V) và Hệ số an toàn (Fs)
       Fs ^
                0      V_max          Thể tích cát nạo vét (m³)

Dữ liệu tính toán được trực quan hóa qua biểu đồ tương quan giữa chiều sâu hạ cực hạn Hcr với bề rộng mái L, cùng đồ thị Fs theo khối lượng cát nạo vét. So sánh với các nghiên cứu quốc tế trên sông Gota và báo cáo quan trắc sạt lở tại An Giang năm 2014-2015, kết quả mô hình Spencer phản ánh độ chính xác cao với sai số vị trí cung trượt thực tế dưới 8.5%. Phát hiện này khẳng định ổn định bờ sông tại ĐBSCL chịu sự chi phối áp đảo của yếu tố hình thái học lòng dẫn sau nạo vét hơn là ma sát nội tại của nền đất sét.

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

Thứ nhất, ban hành quy chế bắt buộc về khoảng cách hành lang khai thác cát an toàn. Khoảng cách mép mỏ cát đến chân bờ sông tự nhiên phải duy trì bán kính đệm R tối thiểu gấp 1.5 đến 2.0 lần chiều sâu hố đào nạo vét h nhằm đảm bảo hệ số Fs không suy giảm dưới 1.4.

Thứ hai, thiết lập hạn ngạch khai thác cát theo từng phân đoạn lòng sông dựa trên chiều sâu cực hạn Hcr. Cần cắt giảm ít nhất 50% sản lượng khai thác hàng năm trên các tuyến sông Tiền và sông Hậu, khống chế tổng khối lượng khai thác toàn vùng dưới 17 triệu m³/năm trong lộ trình 2026-2030 do Ủy ban Nhân dân các tỉnh chỉ đạo thực hiện.

Thứ ba, triển khai hệ thống quan trắc tự động đo sâu hồi âm định kỳ 3 tháng/lần tại 265 điểm nóng sạt lở. Sở Tài nguyên và Môi trường các tỉnh Vĩnh Long, An Giang, Cần Thơ cần tích hợp công nghệ quét địa hình đáy sông 3D để phát hiện sớm các hố xói sâu vượt ngưỡng thiết kế.

Thứ tư, quản lý nghiêm ngặt tải trọng xây dựng trong phạm vi 30 m đến 50 m tính từ đỉnh bờ sông. Các địa phương cần lập phương án di dời các cụm dân cư nhà cấp 4 và kho bãi tập kết vật liệu nặng nằm sát mép nước nhằm triệt tiêu mô men gây trượt phụ tải trong giai đoạn triều kiệt.

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

Cơ quan quản lý nhà nước gồm Sở Tài nguyên và Môi trường, Sở Xây dựng và Chi cục Thủy lợi: Tài liệu cung cấp cơ sở khoa học để thẩm định hồ sơ cấp phép mỏ cát, quy hoạch ranh giới nạo vét và thiết lập bản đồ cảnh báo nguy cơ trượt lở bờ sông.

Kỹ sư địa kỹ thuật và đơn vị tư vấn thiết kế công trình thủy công: Tận dụng quy trình tính đổi sức kháng cắt Su từ chỉ số dẻo PI kết hợp phương pháp Spencer để tính toán ổn định bờ kè, đê bao trên nền bùn sét mềm ĐBSCL một cách nhanh chóng và tin cậy.

Doanh nghiệp khai thác khoáng sản và nhà thầu nạo vét luồng hàng hải: Tham khảo biểu đồ tra cứu chiều sâu nạo vét cực hạn Hcr và góc dốc αcr để lập biện pháp thi công an toàn, tránh gây sạt lở bờ sông lân cận.

Giảng viên, học viên cao học và sinh viên chuyên ngành Kỹ thuật Xây dựng Công trình Thủy, Địa kỹ thuật: Sử dụng làm tài liệu tham khảo chuyên sâu về ứng dụng phần mềm GeoStudio/Plaxis và phương pháp phân tích ứng suất tổng cộng cho nền đất yếu bão hòa.

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

Phương pháp Spencer có ưu điểm gì vượt trội so với các phương pháp cân bằng giới hạn khác? Phương pháp Spencer thỏa mãn nghiêm ngặt cả cân bằng lực (đứng, ngang) và cân bằng mô men cho từng thỏi trượt. Kết quả tính toán cho hệ số Fs nhỏ nhất trong nhóm các phương pháp cân bằng hoàn chỉnh, đảm bảo tính an toàn cao và tránh rủi ro đánh giá quá cao độ ổn định thực tế của bờ sông.

Tại sao có thể sử dụng chỉ số dẻo PI để tính sức kháng cắt Su thay cho thí nghiệm cắt đất trực tiếp? Mô hình toán học của Wroth & Houlsby dựa trên chỉ số dẻo PI cho hệ số tương quan R² đạt 0.944 và độ phân tán COV thấp hơn cắt đất trực tiếp khi đối chiếu với cắt cánh hiện trường. Phương pháp này loại trừ được sai số do xáo động mẫu đất sét yếu trong phòng thí nghiệm.

Biến động mực nước sông ảnh hưởng như thế nào đến nguy cơ sạt lở bờ sông? Khi mực nước sông hạ thấp xuống mức nước kiệt từ -1.5 m đến -2.0 m, áp lực nước tĩnh triệt tiêu làm mất lực đẩy ngược vào mái dốc trong khi áp lực nước lỗ rỗng trong đất chưa kịp thoát, khiến hệ số an toàn Fs suy giảm mạnh từ 18% đến 25%.

Ngưỡng an toàn kỹ thuật nào được áp dụng để xác định giới hạn khai thác cát? Quy chuẩn 22 TCN 262-2000 và tiêu chuẩn FHWA-NHI-06-088 quy định hệ số an toàn Fs tối thiểu là 1.4 cho các mái dốc ổn định dài hạn. Mọi kịch bản khai thác cát làm Fs tụt xuống dưới 1.4 đều bị xếp vào nhóm nguy cơ phá hoại trượt lở.

Làm thế nào để xác định khoảng cách khai thác cát an toàn đối với công trình ven bờ? Khoảng cách an toàn R được xác định bằng hàm tương quan phụ thuộc vào chiều sâu nạo vét h, bề rộng mái L và góc dốc cực hạn αcr. Thông thường, mép ngoài của hố khai thác phải cách mép bờ tối thiểu 1.5 đến 2.0 lần chiều sâu khai thác thiết kế.

Kết luận

  • Xác lập thành công mối tương quan định lượng giữa khối lượng cát khai thác và hệ số suy giảm ổn định bờ sông Fs tại Đồng bằng sông Cửu Long.
  • Khẳng định tính chính xác của phương pháp Spencer và giải pháp xác định sức kháng cắt Su qua chỉ số dẻo PI cho nền sét bão hòa nước trên 95%.
  • Xây dựng bộ thông số chiều sâu cực hạn Hcr, góc dốc cực hạn αcr cho các kịch bản có tải trọng và không tải trọng tại Vĩnh Long, An Giang và Cần Thơ.
  • Đề xuất công thức và bảng tra khoảng cách nạo vét an toàn R nhằm giữ vững hệ số Fs ≥ 1.4 theo quy chuẩn xây dựng hiện hành.
  • Định hướng chuyển giao kết quả mô hình hóa vào công tác cấp phép khai thác khoáng sản và quan trắc bờ sông tại các tỉnh ĐBSCL ngay trong giai đoạn tới. Các cơ quan quản lý và đơn vị tư vấn cần áp dụng ngay khung tính toán này vào quy hoạch thực tế để bảo vệ bền vững tài nguyên và dân sinh ven sông.