Tổng quan nghiên cứu
Khu vực duyên hải Bạc Liêu – Cà Mau sở hữu diện tích tự nhiên 53.206 ha, trong đó hơn 42.000 ha dành cho sản xuất nông nghiệp và nuôi trồng thủy sản với dân số trên 136.740 người. Tuy nhiên, hệ thống hạ tầng thủy lợi tại đây thường xuyên đối mặt với hiểm họa thiên tai nghiêm trọng. Thực tế sau cơn bão Linda, toàn bộ 55 km đê biển cùng hơn 800 km bờ bao địa phương đã bị hư hỏng, tràn bờ và xé rách do nền địa chất công trình đặc biệt phức tạp.
Vấn đề cốt lõi bắt nguồn từ nền đất yếu đồng bằng sông Cửu Long với tầng trầm tích Holoxen dày từ 10 m đến hơn 50 m. Lớp đất sét bùn bề mặt có dung trọng tự nhiên rất thấp khoảng 14,5 kN/m3, độ ẩm tự nhiên cao dao động từ 65% đến 75%, hệ số rỗng $e$ từ 1,5 đến 2,0 và góc ma sát trong $\varphi$ chỉ đạt 3° đến 5°. Khi đắp tuyến đê bằng chính vật liệu sét địa phương có độ bão hòa nước cao, công trình liên tục xảy ra các biến dạng nguy hiểm như lún sụp, phình trồi và trượt trồi cung tròn.
Mục tiêu của luận văn là nghiên cứu toàn diện điều kiện ổn định của tuyến đê bao đắp bằng vật liệu địa phương trên nền đất yếu Bạc Liêu – Cà Mau. Công trình tập trung giải quyết nguyên nhân gây mất ổn định ở 3 mặt cắt đại diện điển hình (MC1, MC1' phía biển và MC2 phía đồng), từ đó đề xuất giải pháp kỹ thuật nâng cao hệ số an toàn ổn định tổng thể cho hệ thống đê điều ven biển. Ý nghĩa thực tiễn của đề tài là cung cấp căn cứ tính toán chuẩn xác, giúp tối ưu hóa kết cấu công trình, giảm thiểu tối đa rủi ro vỡ đê và bảo vệ sinh kế cho hàng vạn hộ dân vùng dự án.
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 nền tảng cơ học đất và công trình trên nền đất yếu, kết hợp 2 lý thuyết nền tảng: lý thuyết trạng thái cân bằng giới hạn bền (Prandtl, Fellenius, Bishop) và lý thuyết sức chống cắt trượt từ biến của N. Maslov. Sức chống cắt của đất được thiết lập qua quan hệ ứng suất pháp tuyến và áp lực nước lỗ rỗng:
$$\tau_t = (\sigma - u) \tan\varphi + c_w$$
Trong đó bao gồm lực dính kết mềm keo nước $\Sigma_w$ và lực dính cấu trúc cứng $c_c$.
Nghiên cứu làm rõ 4 khái niệm trọng tâm:
- Khoáng vật sét Montmorillonite: Chiếm tỷ lệ vượt trội trong đất yếu đồng bằng sông Cửu Long ($Al_2O_3 \cdot 4SiO_2 \cdot H_2O$), có năng lượng hấp phụ tĩnh điện lên tới hàng trăm kN/m2, làm tăng độ trương nở và biến dạng trượt khi ngậm nước.
- Áp lực nước lỗ rỗng và ứng suất hữu hiệu: Hiện tượng áp lực nước lỗ rỗng $u$ phát sinh đột ngột trong quá trình gia tải nhanh khiến ứng suất hữu hiệu $\sigma' = \sigma - u$ suy giảm, dẫn đến phá hoại cắt tức thời.
- Các dạng mất ổn định đê đất: Phân định rõ 3 cơ chế phá hỏng gồm biến dạng lún sụp thân đê, phình trồi chân taluy và trượt trồi cắt qua cung trượt liên tục xé rách nền đất.
- Tải trọng an toàn ($q_{at}$) và tải trọng giới hạn ($q_{gh}$): Ngưỡng tải trọng ban đầu xuất hiện điểm biến dạng dẻo tại độ sâu $z \approx 0,5b$ và sức chịu tải tới hạn của khối đắp.
Phương pháp nghiên cứu
Nghiên cứu sử dụng nguồn dữ liệu khảo sát địa chất thực tế từ 30 hố khoan địa chất (từ HK1 đến HK30) dọc tuyến đê Bạc Liêu – Cà Mau. Cỡ mẫu nghiên cứu gồm hơn 60 mẫu đất nguyên dạng được phân tầng thí nghiệm: lớp 1 sét dẻo mềm (0 – 0,5 m), lớp 2 bùn sét hữu cơ yếu (0,5 – 10 m) và lớp 3 sét nửa cứng (> 10 m).
Phương pháp thí nghiệm bao gồm thí nghiệm nén ba trục cố kết không thoát nước có đo áp lực nước lỗ rỗng (CU), thí nghiệm cắt nhanh không thoát nước (UU) và thí nghiệm cắt cánh quay hiện trường ($C_u \le 3,0$ kN/m2). Dữ liệu sau đó được xử lý bằng phương pháp phần tử hữu hạn (FEM) kết hợp thuật toán phân tích ổn định mái dốc đa lớp STABL6. Lý do lựa chọn STABL6 là phần mềm cho phép mô hình hóa chính xác cung trượt trụ tròn Bishop, tính toán độc lập ứng suất hữu hiệu và áp lực nước lỗ rỗng trên từng phân thỏi đứng ($N_c$), khắc phục hoàn toàn hạn chế của các bảng tính đơn giản hóa trước đây.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Kết quả phân tích kiểm tra độ ổn định của tuyến đê Bạc Liêu – Cà Mau trên các mặt cắt thiết kế cũ đã chỉ ra nhiều bất cập nghiêm trọng về khả năng chịu lực:
- Mất ổn định theo tải trọng an toàn phân bố hình thang: Khi kiểm tra theo công thức giải tích đàn hồi - dẻo, hệ số an toàn tại mặt cắt mái phía đồng MC2 chỉ đạt $F = 0,969 < 1,0$; tại mặt cắt bờ biển MC1 đạt $F = 0,731 < 1,0$; đặc biệt tại mặt cắt MC1' hệ số $F$ tụt xuống mức nguy hiểm $0,373 < 1,0$. Toàn bộ các hệ số này đều không thỏa mãn điều kiện an toàn $F \ge 1,0$.
- Mất ổn định theo tải trọng an toàn phân bố hình chữ nhật: Giá trị hệ số an toàn tương ứng tại MC2 và MC1 chỉ đạt $F = 0,684$, còn tại MC1' chỉ đạt $F = 0,348$, thấp hơn quy chuẩn cho phép từ 31,6% đến 65,2%.
- Không đảm bảo theo tải trọng giới hạn hình băng Prandtl: Hệ số an toàn tính toán cho MC1 và MC2 đạt $F = 1,8$ (không đạt yêu cầu quy phạm $F = 2,0 - 4,0$ đối với công trình bảo vệ dân sinh quan trọng trên đất yếu); riêng MC1' hệ số an toàn chỉ đạt $F = 0,5$, nền đất rơi vào trạng thái phá hoại cắt hoàn toàn.
- Phát hiện lỗi cốt lõi trong quy trình tính toán cũ: Hồ sơ thiết kế ban đầu chỉ sử dụng chương trình macro V.Qbasic/Excel đơn giản để tính trượt cung tròn theo phương pháp áp lực trọng lực Terzaghi, hoàn toàn bỏ qua áp lực nước lỗ rỗng, sự suy giảm góc ma sát trong ($\varphi \approx 1^\circ - 3^\circ$) và lực dính kết cực nhỏ ($C \le 0,7$ kN/cm2) của lớp bùn sét dày 9,5 m dưới đáy đê.
Thảo luận kết quả
Khi thể hiện dữ liệu qua biểu đồ cung trượt quét từ phần mềm STABL6, tâm trượt nguy hiểm nhất luôn xuất hiện sâu bên trong lớp bùn sét hữu cơ (lớp 2), kéo dài từ đỉnh mái đê xuyên qua đáy móng và trồi lên ngoài chân taluy. Kết hợp với biểu đồ sóng Hindcast và đường bão hòa thấm qua thân đê đất đồng chất, áp lực sóng leo cùng gradien thủy lực dòng thấm $j = i \cdot \gamma_w$ tạo ra ứng suất trượt vượt xa giới hạn bền cấu trúc của đất sét Montmorillonite.
So với các công trình nghiên cứu đê điều tại miền Bắc, nền đất duyên hải Bạc Liêu – Cà Mau có độ sệt $B > 1,0$ (trạng thái chảy). Do đó, việc đắp đê dốc với hệ số mái dốc $m = 2$ mà không có biện pháp gia cường đã khiến vùng biến dạng dẻo mở rộng vượt quá 60% bề rộng đáy đê, giải thích hiện tượng đê liên tục bị trượt trồi và sập vỡ trong thực tế vận hành.
Đề xuất và khuyến nghị
Dựa trên kết quả mô phỏng và tính toán cơ học đất, 4 giải pháp công trình và phi công trình trọng tâm được đề xuất nhằm đảm bảo ổn định lâu dài:
- Áp dụng kết cấu bệ phản áp hai bên chân đê: Xây dựng bệ phản áp bằng đất đắp địa phương tại các mặt cắt xung yếu MC1, MC1' và MC2 nhằm gia tăng mômen giữ, triệt tiêu ứng suất trượt gây phình trồi. Mục tiêu kỹ thuật là nâng hệ số an toàn ổn định trượt $K_{min}$ lên $\ge 1,35$ trong giai đoạn vận hành thấm ổn định.
- Gia cường thân và đáy đê bằng vải địa kỹ thuật (Geotextile): Trải từ 1 đến 2 lớp vải địa kỹ thuật có cường độ chịu kéo cao tại mặt tiếp xúc giữa thân đê và nền bùn sét. Lớp vải địa đảm nhận toàn bộ ứng suất kéo tiếp xúc, phân phối đều tải trọng khối đắp, ngăn ngừa hiện tượng xé rách thân đê và tăng tốc độ cố kết của nền đất.
- Tối ưu hóa hình học mặt cắt và đắp thoải mái dốc: Giảm độ dốc mái đê phía biển từ $m = 2$ xuống $m = 3$ hoặc $m = 3,5$, kết hợp thiết kế thềm giảm sóng (TFS) để tiêu tán năng lượng sóng biển tác dụng lên mái thượng lưu, hạ thấp cao trình đường bão hòa trong thân đê.
- Quy chuẩn hóa quy trình thi công phân kỳ gia tải: Ban hành hướng dẫn thi công đắp đê nhiều đợt với chiều dày mỗi lớp đắp đầm nén không vượt quá 0,5 m. Khống chế tốc độ đắp nghiêm ngặt để áp lực nước lỗ rỗng kịp tiêu tán, đảm bảo độ tăng sức chống cắt không thoát nước $\Delta C_u$ của đất nền bắt kịp tải trọng gia tăng.
Timeline triển khai khuyến nghị nên chia làm 2 giai đoạn: Giai đoạn 1 (năm 2002 - 2004) hoàn thành xử lý bệ phản áp và trải vải địa kỹ thuật cho 55 km đê biển xung yếu; Giai đoạn 2 (sau năm 2004) đồng bộ hóa hệ thống bờ bao ngăn mặn và kênh trục dẫn nước ngọt do Sở Nông nghiệp & Phát triển Nông thôn cùng Ban Quản lý Dự án Thủy lợi chủ trì thực hiện.
Đối tượng nên tham khảo luận văn
- Kỹ sư thiết kế công trình thủy lợi, giao thông và đê điều ven biển: Cung cấp cẩm nang tra cứu chỉ tiêu cơ lý đất yếu đồng bằng sông Cửu Long, thuật toán phân tích ổn định mái đê bằng STABL6 và phương pháp tính toán bệ phản áp chuẩn xác.
- Ban quản lý dự án và các nhà thầu thi công hạ tầng đất nền: Nắm bắt các chỉ dẫn kỹ thuật thi công trải vải địa kỹ thuật, biện pháp lu lèn đất sét ngậm nước tại chỗ và kỹ thuật kiểm soát tốc độ đắp nền đê tránh sự cố trượt sạt.
- Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Địa kỹ thuật - Công trình trên đất yếu: Tài liệu tham khảo học thuật giá trị về lý thuyết cân bằng giới hạn, cơ chế phá hoại của khoáng vật Montmorillonite và mô hình tương tác đất - vải địa kỹ thuật.
- Cơ quan quản lý nhà nước về phòng chống thiên tai và biến đổi khí hậu: Nguồn tư liệu thực chứng giúp quy hoạch tuyến đê biển, đê bao ngăn mặn thích ứng với chế độ thủy triều biên độ cao (2,5 – 3,0 m) và nước biển dâng tại bán đảo Cà Mau.
Câu hỏi thường gặp
Nguyên nhân cơ bản khiến đê bao đắp bằng đất địa phương tại Bạc Liêu – Cà Mau dễ bị trượt vỡ là gì?
Nguyên nhân cốt lõi là nền đất chứa tầng bùn sét dày 0,5 – 10 m với hàm lượng khoáng Montmorillonite cao, góc ma sát trong rất thấp ($\varphi \approx 1^\circ - 3^\circ$), độ sệt chảy $B > 1,0$. Khi đắp đất đê lên trên, áp lực nước lỗ rỗng tăng nhanh làm triệt tiêu ứng suất hữu hiệu, dẫn đến mất ổn định dạng trượt trồi.
Vì sao các phương pháp tính toán ổn định đê trước đây tại địa phương không đạt độ tin cậy?
Các hồ sơ cũ sử dụng macro V.Qbasic/Excel đơn giản hóa giả thiết cung trượt tròn Fellenius thông thường, bỏ qua hoàn toàn áp lực nước lỗ rỗng, sự phân bố lại ứng suất dẻo và tính thấm không ổn định của đất sét dẻo chảy, dẫn đến đánh giá sai lệch hệ số an toàn thực tế ($F < 1,0$).
Bệ phản áp đóng vai trò kỹ thuật gì trong việc ổn định tuyến đê?
Bệ phản áp được đắp mở rộng ở chân mái taluy đóng vai trò đối trọng, tạo ra mômen giữ chống lại mômen gây trượt của khối thân đê chính. Giải pháp này giúp giữ phẳng chân đê, ngăn chặn hiện tượng đẩy phình trồi đất yếu và nâng hệ số an toàn $K_{min}$ lên trên 1,35.
Vải địa kỹ thuật hỗ trợ ổn định nền đê trên đất yếu theo nguyên lý nào?
Vải địa kỹ thuật hoạt động như một tầng phân cách và gia cường chịu kéo cường độ cao. Khi nền đất chịu tải trọng đắp, vải địa hấp thu ứng suất tiếp tuyến, ngăn ngừa sự xé rách cục bộ của thân đê, đồng thời lọc và hỗ trợ tiêu thoát nước lỗ rỗng nhanh chóng.
Tại sao tốc độ đắp đê lại quyết định trực tiếp đến sức chịu tải của nền móng?
Đất sét yếu có hệ số thấm cực nhỏ ($k \approx 10^{-6} - 10^{-8}$ cm/s), nước lỗ rỗng thoát ra rất chậm. Nếu tốc độ đắp quá nhanh, áp lực nước lỗ rỗng tích tụ lớn làm giảm sức kháng cắt; đắp theo từng phân đoạn thời gian giúp đất nền kịp cố kết và gia tăng độ bền cắt theo thời gian.
Kết luận
- Luận văn đã tổng hợp và xác lập bức tranh địa chất công trình chi tiết cho vùng đất yếu Bạc Liêu – Cà Mau với 5 phân vùng địa kỹ thuật đặc thù và tầng sét mềm dày đến 10 m.
- Chứng minh nguyên nhân mất ổn định của các tuyến đê hiện hữu qua tính toán 3 mặt cắt điển hình MC1, MC1', MC2 với hệ số an toàn thực tế đều dưới ngưỡng cho phép ($F = 0,348 - 0,969$).
- Khẳng định sự cần thiết phải tích hợp áp lực nước lỗ rỗng và ứng suất hữu hiệu thông qua phần mềm phần tử hữu hạn chuyên dụng như STABL6 thay thế các bảng tính đơn giản.
- Đề xuất thành công 3 giải pháp kỹ thuật công trình đồng bộ: mở rộng bệ phản áp, gia cố vải địa kỹ thuật thân đê và điều chỉnh độ dốc mái $m = 3$ kết hợp thềm giảm sóng.
- Đưa ra quy trình công nghệ thi công đắp phân đợt kiểm soát độ lún cố kết, đóng góp cơ sở khoa học quan trọng cho chiến lược củng cố đê biển quốc gia vùng đồng bằng sông Cửu Long. Các đơn vị tư vấn và ban quản lý dự án cần áp dụng ngay khung tính toán và giải pháp này vào các dự án nâng cấp hạ tầng ven biển để bảo đảm tính bền vững công trình.