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 báo động khi mỗi năm có khoảng 34 triệu m3 trầm tích, tương đương 55 triệu tấn bị khai thác từ lòng sông, trong đó cát chiếm đến 90%. Toàn vùng hiện ghi nhận 265 điểm sạt lở bờ sông và bờ biển với tổng chiều dài hơn 450 km. Hiện tượng xói lở làm biến mất khoảng 500 ha đất màu mỡ mỗi năm với tốc độ lấn sâu 30-40 m/năm. Hoạt động khai thác cát không kiểm soát tại 126 mỏ cát được cấp phép dọc các nhánh sông chính đã làm thay đổi nghiêm trọng hình thái dòng dẫn, gây suy giảm độ ổn định của các mái dốc bờ sông trên nền địa chất bùn sét yếu. Hậu quả thực tế thể hiện rõ qua các sự cố sạt lở nghiêm trọng, điển hình như vụ sạt lở ngày 29/10/2012 tại xã An Bình, huyện Long Hồ, tỉnh Vĩnh Long làm trôi 8.000 m2 đất cùng 27 bè và hầm cá, gây thiệt hại gần 4 tỷ đồng; hay vụ sạt lở ngày 19/6/2015 tại quận Cái Răng, TP Cần Thơ làm đứt gãy 30 m đường giao thông và ăn sâu vào đất liền 5-7 m.

Mục tiêu cốt lõi của luận văn là nghiên cứu định lượng mối tương quan giữa khối lượng khai thác cát và mức độ ổn định của bờ sông Đồng bằng sông Cửu Long. Đề tài tập trung xác định hệ số an toàn (Fs), chiều sâu xói cực hạn (Hcr), góc dốc cực hạn và khoảng cách khai thác cát an toàn. Phạm vi nghiên cứu bao gồm 3 địa bàn trọng điểm: sông Cổ Chiên (Vĩnh Long), sông Hậu (An Giang) và sông Cần Thơ (Cần Thơ), kết hợp số liệu kiểm chứng tại Trà Vinh với chuỗi dữ liệu thủy văn thực đo giai đoạn 2013-2014. Nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp công cụ dự báo xói lở, phục vụ quy hoạch khai thác khoáng sản bền vững và bảo vệ an toàn sinh kế cho hơn 1 triệu người dân vùng hạ lưu châu thổ trước năm 2050.

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 đồng thời các lý thuyết cơ học đất và địa kỹ thuật công trình hiện đại:

  • Phương pháp Cân bằng Giới hạn (Limit Equilibrium Method) theo thuật toán Spencer (1967): Đây là phương pháp tối ưu khi thỏa mãn đầy đủ cả hai điều kiện cân bằng lực (phương ngang, phương đứng) và cân bằng mô men đối với từng thỏi phân tố trượt. Tiêu chuẩn FHWA-NHI-06-088 khuyến nghị áp dụng phương pháp Spencer để đảm bảo kết quả chính xác và chặt chẽ hơn các phương pháp đơn giản như Fellenius hay Bishop.
  • Phương pháp Phần tử Hữu hạn (Finite Element Method): Ứng dụng kỹ thuật giảm dần các thông số độ bền kháng cắt (c-phi reduction) để xác định trạng thái ứng suất - biến dạng và tìm hệ số an toàn ổn định khi chuyển vị đạt giá trị phá hoại cực hạn.
  • Lý thuyết phân tích ứng suất tổng cộng cho nền đất dính bão hòa nước: Áp dụng với điều kiện đất sét yếu có độ bão hòa nước lớn hơn 95%, giả định góc ma sát trong tổng cộng bằng 0 và sức kháng cắt không thoát nước Su đóng vai trò là lực dính thuần túy.
  • Mô hình xác định sức chống cắt không thoát nước Su: Kết hợp mô hình thực nghiệm Worth and Houlsby (1985) dựa trên chỉ số dẻo PI thông qua công thức toán học và mô hình trượt Mohr-Coulomb tiêu chuẩn.

Các khái niệm then chốt được chuẩn hóa gồm: Hệ số an toàn mái dốc Fs (với Fs tối thiểu bằng 1.4 theo quy chuẩn 22 TCN 262 - 2000 của Bộ Giao thông Vận tải), chiều sâu bờ dốc cực hạn Hcr (tương ứng Fs đạt ngưỡng tới hạn 1.1), hệ số hạ sâu cực hạn Rcr và góc dốc cực hạn của mái bờ sông.

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

  • Nguồn dữ liệu và quy mô khảo sát: Nghiên cứu thu thập dữ liệu địa chất công trình từ hệ thống các lỗ khoan khảo sát sâu đến 30 m dọc tuyến bờ sông Cổ Chiên, sông Hậu và sông Cần Thơ. Cỡ mẫu nghiên cứu gồm 25 bộ chỉ tiêu cơ lý đất nguyên dạng (14 mẫu tại Vĩnh Long và 11 mẫu tại Trà Vinh) nhằm phân tích, đối sánh các chỉ số sức kháng cắt Su. Dữ liệu thủy văn được trích xuất từ các trạm quan trắc Mỹ Thuận, Long Xuyên và Cần Thơ với chuỗi đo liên tục 730 ngày trong hai năm 2013-2014.
  • Phương pháp chọn mẫu: Lựa chọn mẫu đất sét yếu bão hòa đại diện (thuộc nhóm đất dẻo mềm đến chảy, phân loại sét bùn OH-OL và CL-CH) tại các khu vực bờ sông chịu tác động trực tiếp của hoạt động nạo vét cát.
  • Phương pháp phân tích và lý do lựa chọn: Sử dụng bộ phần mềm chuyên ngành GeoStudio (mô đun SLOPE/W) và Plaxis 2D. Lý do lựa chọn thuật toán Spencer là vì phương pháp này cho kết quả hệ số an toàn có độ tin cậy cao, thiên về an toàn và đạt hệ số tương quan R2 = 0.9887 so với mô hình phần tử hữu hạn. Bên cạnh đó, việc xác định Su thông qua chỉ số dẻo PI được lựa chọn vì mang lại tỷ số giá trị trung bình đạt 0.944 (tiệm cận chuẩn 1.0) và hệ số biến thiên COV chỉ 0.23, thấp hơn đáng kể so với phương pháp cắt đất trực tiếp trong phòng thí nghiệm (COV = 0.38).
  • Tiến trình nghiên cứu: Quá trình khảo sát hiện trường, thu thập dữ liệu thủy văn và thực hiện các mô phỏng cơ học đất được hoàn thành xuyên suốt giai đoạn 2015-2016.

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

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

  • Hoạt động nạo vét cát làm sụt giảm nghiêm trọng hệ số an toàn bờ sông: Mô phỏng hình thái thay đổi địa hình lòng sông cho thấy kịch bản gia tăng góc dốc bờ do xói mòn thủy động nguy hiểm hơn nhiều so với kịch bản hạ sâu đáy phẳng. Tại mặt cắt sông Cổ Chiên với bề rộng mái dốc L = 40 m và chiều sâu ban đầu H = 10 m, hệ số an toàn Fs giảm mạnh từ 3.292 xuống 2.292 (sụt giảm 30.38% độ ổn định) khi chiều sâu hạ thấp 4 m.
  • Xác lập giới hạn chiều sâu cực hạn và góc dốc cực hạn: Tại mặt cắt sông thuộc tỉnh Vĩnh Long với bề rộng mái dốc L = 90 m và chiều sâu H = 14.96 m ở điều kiện không tải, chiều sâu bờ dốc cực hạn đạt Hcr = 22.5 m (độ hạ sâu cực hạn Delta Hcr = 7.54 m) trước khi mái dốc chạm ngưỡng phá hoại Fs = 1.1. Khi có thêm tải trọng công trình xây dựng phân bố đều trên đỉnh bờ, chiều sâu cực hạn giảm xuống còn 18.2 m (mức suy giảm độ sâu an toàn lên tới 19.11%).
  • Mối quan hệ tỷ lệ nghịch chặt chẽ giữa thể tích khai thác cát và hệ số an toàn: Kết quả mô hình hóa chứng minh thể tích cát lấy đi tỷ lệ nghịch tuyến tính với hệ số Fs. Để duy trì hệ số Fs >= 1.4 theo quy chuẩn kỹ thuật, khối lượng cát tối đa được phép khai thác tại phân đoạn thân cát khảo sát không được vượt quá 120.000 m3 đối với bề rộng mái dốc 40 m. Đồng thời, khoảng cách khai thác cát an toàn R từ mép bờ sông phải duy trì tối thiểu 45-60 m tương ứng với chiều sâu khai thác h từ 5-10 m.
  • Tác động bất lợi từ biên độ dao động mực nước sông: Mực nước thấp nhất trong mùa kiệt tại trạm Mỹ Thuận (-1.85 m) và trạm Long Xuyên (-1.20 m) làm triệt tiêu lực đẩy nổi tĩnh học, khiến hệ số an toàn Fs của bờ sông giảm thêm từ 15% đến 22% so với điều kiện mực nước triều dâng thông thường.

Thảo luận kết quả

  • Cơ chế phá hoại: Quá trình nạo vét cát làm thay đổi cấu trúc chân mái dốc, tạo nên các hố sâu cục bộ và làm tăng góc dốc tự nhiên alpha. Khi áp lực thủy tĩnh suy giảm trong các kỳ triều rút kiệt, mô men gây trượt MD tăng đột biến. Trong khi đó, các tầng đất bùn sét yếu có sức kháng cắt Su thấp (chỉ từ 9.3 kPa đến 25 kPa) không đủ tạo ra mô men kháng trượt MR cần thiết, dẫn đến sự trượt lở mái bờ theo các cung trượt tròn sâu.
  • So sánh học thuật: Nghiên cứu giải quyết toàn diện khoảng trống của các công bố quốc tế trước đây như nghiên cứu của Luo và cộng sự (2007) hay báo cáo của Ủy ban Nước Trung ương Ấn Độ (2007) vốn chỉ dừng lại ở việc đánh giá biến đổi thủy văn chung chung mà chưa lượng hóa được chiều sâu cực hạn Hcr và góc dốc cực hạn alpha_cr theo từng cấp mực nước và tải trọng bờ.
  • Mô tả dữ liệu trực quan: Các kết quả phân tích trong luận văn được hệ thống hóa trực quan thông qua biểu đồ hồi quy tuyến tính R2 = 0.9887 giữa phương pháp Spencer và phần tử hữu hạn, kết hợp với các bảng tra cứu ma trận tham số (H, L, WL, Fs) và các lát cắt mô phỏng hình dạng cung trượt nguy hiểm trên phần mềm GeoStudio.

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

  • Ban hành quy chuẩn khoảng cách và độ sâu khai thác cát an toàn: Ủy ban Nhân dân và Sở Tài nguyên & Môi trường các tỉnh Vĩnh Long, An Giang, Cần Thơ cần quy định bắt buộc khoảng cách nạo vét mỏ cát an toàn tối thiểu từ 45 m đến 60 m tính từ mép bờ hoặc chân kè hiện hữu. Đồng thời, khống chế chiều sâu nạo vét h không vượt quá 5-8 m nhằm duy trì hệ số an toàn bờ dốc Fs luôn lớn hơn hoặc bằng 1.4. Thời gian hoàn thành việc rà soát và cắm mốc quy hoạch mỏ trong vòng 6 tháng.
  • Kiểm soát nghiêm ngặt tải trọng xây dựng trên đỉnh bờ sông: Sở Xây dựng các địa phương cần thiết lập hành lang an toàn bảo vệ bờ sông với bề rộng cách mép bờ từ 20 m đến 30 m. Tuyệt đối cấm xây dựng các công trình kiên cố hoặc chất tải vượt quá 15 kN/m2 trong phạm vi cung trượt nguy hiểm đã được xác định. Kế hoạch thanh tra và xử lý vi phạm cần được thực hiện định kỳ 3 tháng/lần.
  • Lắp đặt hệ thống trạm quan trắc tự động và cảnh báo sớm: Sở Khoa học & Công nghệ chủ trì phối hợp cùng các viện nghiên cứu chuyên ngành lắp đặt 10-15 trạm quan trắc biến dạng tự động và đo áp lực nước lỗ rỗng tại các điểm nóng sạt lở trọng yếu. Đồng thời, ứng dụng công nghệ đo sâu quét sườn đa tia để cập nhật bản đồ địa hình lòng sông theo chu kỳ 30 ngày/lần. Lộ trình triển khai thí điểm hoàn thành trong 12 tháng.
  • Cắt giảm sản lượng và phân bổ hạn ngạch khai thác theo sức chịu tải lòng dẫn: Bộ Tài nguyên & Môi trường cần rà soát, cắt giảm từ 30% đến 40% tổng sản lượng khai thác cát thương mại trên toàn vùng ĐBSCL (hiện ở mức 34 triệu m3/năm) để đưa về mức cân bằng với lượng trầm tích bồi đắp thực tế từ thượng nguồn. Kế hoạch cắt giảm được thực hiện theo lộ trình từng năm trong giai đoạn 2026-2030.

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

  • Cơ quan quản lý nhà nước (Sở Tài nguyên & Môi trường, Sở Nông nghiệp & Phát triển Nông thôn các tỉnh ĐBSCL): Nắm bắt cơ sở khoa học định lượng để thẩm định báo cáo đánh giá tác động môi trường của các dự án khai thác cát, xác lập hành lang bảo vệ bờ sông và quản lý hiệu quả 265 điểm sạt lở xung yếu.
  • Doanh nghiệp khai thác khoáng sản và nhà thầu nạo vét đường thủy: Ứng dụng quy trình tính toán khối lượng cát cho phép và góc dốc nạo vét giới hạn nhằm xây dựng phương án thi công an toàn, duy trì hệ số ổn định Fs >= 1.4, hạn chế tối đa nguy cơ sạt lở gây thiệt hại tài sản hàng tỷ đồng.
  • Kỹ sư tư vấn thiết kế địa kỹ thuật và công trình thủy: Áp dụng phương pháp tính nhanh sức kháng cắt không thoát nước Su từ chỉ số dẻo PI và quy trình mô phỏng bằng phương pháp Spencer trên GeoStudio để tính toán ổn định bờ kè, đê bao trên nền đất sét bão hòa có độ sâu 10-30 m.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng / Địa kỹ thuật: Sử dụng làm tài liệu tham khảo chuyên sâu về phương pháp cân bằng giới hạn, phân tích ứng suất tổng cộng và bộ dữ liệu thực nghiệm địa chất thủy văn chi tiết tại khu vực Tây Nam Bộ.

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

  • Tại sao nghiên cứu lại ưu tiên lựa chọn phương pháp Spencer trong phân tích ổn định bờ sông? Phương pháp Spencer (1967) thỏa mãn đồng thời tất cả các phương trình cân bằng lực thẳng đứng, lực nằm ngang và mô men quay của từng phân tố trượt. Kết quả kiểm chứng thực nghiệm cho thấy phương pháp Spencer cho hệ số an toàn Fs = 2.072, mang tính chặt chẽ và thiên về an toàn hơn hẳn các phương pháp Fellenius hay Bishop, đồng thời đạt độ tương quan tuyệt vời R2 = 0.9887 so với mô hình phần tử hữu hạn FEM.

  • Việc tính toán sức chống cắt Su qua chỉ số dẻo PI có đảm bảo độ chính xác so với thí nghiệm hiện trường không? Phương pháp tính Su qua chỉ số dẻo PI theo công thức của Worth and Houlsby cho kết quả rất tin cậy. Dữ liệu thực nghiệm trên 14 mẫu tại Vĩnh Long và 11 mẫu tại Trà Vinh chứng minh tỷ số giá trị trung bình của phương pháp này đạt 0.944 và hệ số biến thiên COV chỉ ở mức 0.23 (thấp hơn nhiều so với mức 0.38 của thí nghiệm cắt đất trực tiếp), giúp loại bỏ sai số do xáo động mẫu đất dính bão hòa.

  • Hoạt động nạo vét cát lòng sông gây mất ổn định bờ theo cơ chế nào là nguy hiểm nhất? Cơ chế hạ sâu lòng sông theo hướng gia tăng góc dốc mái bờ là nguy hiểm nhất. Dưới tác động kết hợp của máy hút cát và dòng chảy thủy động, góc dốc tăng lên làm gia tăng nhanh chóng mô men gây trượt, khiến hệ số an toàn Fs sụt giảm tới 30.38% (từ mức 3.292 xuống còn 2.292), gây sụt lún mái dốc nhanh hơn nhiều so với kịch bản hạ sâu đáy phẳng đồng đều.

  • Ngưỡng hệ số an toàn và khoảng cách tối thiểu khi khai thác cát gần bờ được xác định như thế nào? Căn cứ theo tiêu chuẩn 22 TCN 262 - 2000, hệ số an toàn ổn định nhỏ nhất cho phép của bờ sông là Fs = 1.4. Dựa trên điều kiện địa chất và thủy văn khu vực ĐBSCL, nghiên cứu chỉ ra rằng để đảm bảo Fs >= 1.4, tàu hút cát phải duy trì khoảng cách cách bờ tối thiểu từ 45 m đến 60 m tương ứng với chiều sâu nạo vét h từ 5 m đến 10 m.

  • Sự hiện diện của tải trọng công trình nhà ở trên đỉnh bờ ảnh hưởng ra sao đến chiều sâu khai thác cực hạn? Tải trọng nhà cấp 4 phân bố trên bờ làm gia tăng đáng kể tải trọng gây trượt. Tại mặt cắt sông Cổ Chiên với bề rộng mái dốc L = 90 m, sự xuất hiện của tải trọng xây dựng làm giảm chiều sâu cực hạn Hcr từ 22.5 m xuống còn 18.2 m (giảm khoảng 19.11%), khiến bờ sông dễ bị sạt lở ở độ sâu khai thác cạn hơn so với bờ tự nhiên không có tải.

Kết luận

  • Luận văn đã xây dựng thành công mô hình toán học định lượng chính xác mối tương quan giữa khối lượng khai thác cát và hệ số an toàn ổn định Fs trên nền đất sét yếu Đồng bằng sông Cửu Long.
  • Khẳng định tính chuẩn xác của phương pháp Spencer trong phân tích ổn định trượt bờ sông và chứng minh tính ưu việt của công thức xác định sức kháng cắt Su từ chỉ số dẻo PI với hệ số biến thiên COV chỉ 0.23.
  • Thiết lập bộ thông số kỹ thuật giới hạn bao gồm chiều sâu cực hạn Hcr từ 18.2 m đến 22.5 m, góc dốc cực hạn và khoảng cách khai thác cát an toàn từ 45 m đến 60 m.
  • Đóng góp luận cứ khoa học thực tiễn vững chắc phục vụ lộ trình cắt giảm 30-40% sản lượng khai thác cát quá mức, góp phần bảo vệ hơn 450 km bờ sông trước thảm họa xói lở.
  • Trong giai đoạn 2026-2030, các nghiên cứu tiếp theo cần mở rộng mô hình hóa tích hợp áp lực động dòng lũ và kịch bản nước biển dâng; các cơ quan hữu quan cần nhanh chóng đưa bộ công cụ tính toán này vào quy chuẩn cấp phép khai thác mỏ thực tế.