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 đổi địa chất và thủy văn vô cùng nghiêm trọng. Theo số liệu quan trắc từ năm 2010 đến nay, toàn vùng đã ghi nhận 562 điểm sạt lở bờ sông với tổng chiều dài vượt quá 786 km, tập trung chủ yếu dọc theo các tuyến sông huyết mạch như sông Tiền, sông Hậu, sông Vàm Cỏ Đông và sông Vàm Cỏ Tây. Sự kết hợp giữa điều kiện địa chất bùn sét yếu, biên độ nhật triều biến thiên 2 lần mỗi ngày, lưu lượng dòng chảy lũ và tác động từ mật độ giao thông thủy tải trọng lớn đã làm gia tăng tốc độ xói mòn bờ.

Vấn đề cốt lõi đặt ra trong kỹ thuật xây dựng công trình bảo vệ bờ là tính toán chính xác sự tương tác giữa kết cấu tường chắn và nền đất yếu. Trong thực tế, việc san lấp mặt bằng sau tường kè gây gia tăng tải trọng thẳng đứng, thúc đẩy quá trình lún cố kết thấm và phát sinh áp lực đất bổ sung, dễ dẫn đến hiện tượng trượt cung tròn hoặc nghiêng lệch công trình. Các phương pháp tính toán truyền thống thường chỉ xét tải trọng tác dụng lên cọc dưới dạng lực tập trung tĩnh đơn thuần, bỏ qua trạng thái ứng suất ban đầu và sự làm việc đồng thời của hệ cọc - đất - tường kè.

Luận văn tập trung nghiên cứu giải pháp kết cấu tường kè bản góc chữ L bằng bê tông cốt thép đặt trên hệ móng cọc ma sát tại tuyến sông Vàm Cỏ Tây, tỉnh Long An. Mục tiêu cụ thể là mô phỏng toàn diện trạng thái ứng suất - biến dạng của hệ thống tường cọc theo từng giai đoạn thi công và theo thời gian cố kết dài hạn. Kết quả nghiên cứu đóng vai trò quan trọng trong việc tối ưu hóa kích thước hình học cấu kiện, kiểm soát chuyển vị ngang công trình dưới ngưỡng 40 mm và nâng cao hệ số ổn định tổng thể lên trên 1,30, mang lại hiệu quả kinh tế và an toàn kỹ thuật cao cho các dự án chỉnh trị sông tại địa phương.

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 được xây dựng dựa trên sự kết hợp chặt chẽ giữa cơ học đất cổ điển và cơ học môi trường liên tục hiện đại:

  • Lý thuyết cân bằng giới hạn Rankine và Coulomb: Sử dụng để xác định áp lực đất chủ động ($E_a$) và áp lực đất bị động ($E_p$) tác dụng lên bản tường đứng cao 2,55 m và bản đáy rộng 2,0 m. Lý thuyết tính toán xét đến lực dính đơn vị $c$, góc ma sát trong $\varphi$, góc ma sát ngoài $\delta$ giữa lưng tường và đất, cũng như phụ tải phân bố đều $q$ trên mặt đất đắp.
  • Mô hình phản lực nền đàn hồi Winkler: Mô phỏng sự làm việc của cọc bê tông cốt thép chịu tải trọng ngang trong môi trường đất. Đất quanh cọc được xem như hệ thống lò xo đàn hồi tuyến tính với hệ số nền theo phương ngang $C_{zy} = K \cdot z$, trong đó hệ số tỷ lệ $K$ biến thiên từ 650 đến 13.000 kN/m$^4$ tùy theo trạng thái dẻo mềm hoặc nửa cứng của từng lớp đất.
  • Lý thuyết biến dạng môi trường liên tục và mô hình đàn dẻo Mohr-Coulomb: Thiết lập phương trình cân bằng gia số ứng suất - biến dạng kết hợp thuật toán tích phân hàm ẩn ($\theta = 1$) theo nghiên cứu của Vermeer, giúp giải quyết triệt để tính phi tuyến hình học và phi tuyến vật liệu trong quá trình chuyển tiếp trạng thái dẻo.

Các khái niệm trọng tâm bao gồm: áp lực đất chủ động phi tuyến, hệ số phản lực nền ngang, chuyển vị xoay đầu cọc, quá trình cố kết thấm Biot và mặt trượt phá hoại dẻo tổng thể.

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

Nguồn dữ liệu địa kỹ thuật phục vụ nghiên cứu được trích xuất từ hồ sơ khảo sát địa chất công trình của Trung tâm Nghiên cứu Chỉnh trị sông và Phòng chống thiên tai (thuộc Viện Khoa học Thủy lợi Miền Nam) tại đoạn kè thị xã Tân An, tỉnh Long An. Địa tầng khu vực nghiên cứu có chiều sâu khảo sát trên 26,6 m, được chia thành 4 lớp đất đặc trưng:

  • Lớp 1a: Đất cát san lấp lẫn tạp chất dày 1,5 - 2,0 m.
  • Lớp 1b: Sét màu xám xanh đốm đỏ ở trạng thái dẻo mềm, khả năng chịu lực rất kém.
  • Lớp 2: Sét lẫn cát màu vàng nâu xám xanh, trạng thái nửa cứng.
  • Lớp 3: Cát lẫn sét màu vàng xám xanh, trạng thái dẻo đến dẻo cứng.

Mô hình phân tích sử dụng phương pháp phần tử hữu hạn 2D thông qua phần mềm chuyên ngành Plaxis. Cỡ mẫu mô phỏng bao gồm mặt cắt địa chất đại diện với lưới phần tử tam giác 15 nút, phân bố độ mịn cao tại khu vực tiếp xúc giữa cọc và nền đất. Phương pháp phần tử hữu hạn được lựa chọn thay cho các công thức giải tích cổ điển vì có khả năng mô tả chính xác ứng xử phi tuyến của đất mềm, phản ánh chính xác sự thay đổi áp lực nước lỗ rỗng và mô phỏng sát thực tế trình tự 5 bước thi công: đóng cọc bê tông cốt thép dài 14 m, đổ bê tông đài kè chữ L, đắp đất cát san lấp sau lưng tường đến cao độ +2,60 m, hạ mực nước sông xuống cao độ thiết kế +0,30 m và theo dõi cố kết thấm trong thời gian 1.000 ngày.

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

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

Quá trình mô phỏng số và phân tích so sánh giữa phương pháp phần tử hữu hạn với phương pháp giải tích truyền thống mang lại những phát hiện quan trọng:

  • Sự phân bố phi tuyến của áp lực đất sau kè: Khác với biểu đồ phân bố tam giác tuyến tính theo lý thuyết Coulomb, áp lực đất chủ động thực tế tác dụng lên bản lưng tường kè giảm khoảng 18% ở phần nửa trên thân tường do biến dạng uốn đàn hồi, nhưng lại tăng cục bộ khoảng 22% tại vị trí chân tường nơi tiếp giáp với bản đáy móng cọc.
  • Biến dạng và chuyển vị gia tăng theo thời gian: Chuyển vị ngang tại đỉnh tường kè ngay sau khi hoàn thành đắp đất đạt 2,8 cm. Tuy nhiên, sau chu kỳ cố kết thấm 300 đến 1.000 ngày của lớp sét mềm, chuyển vị ngang đỉnh kè tăng thêm 35%, đạt giá trị cực đại 3,8 cm khi mực nước sông rút cạn xuống mức kiệt +0,40 m.
  • Sự phân chia nội lực không đồng đều trong hệ cọc: Hàng cọc phía trước (sát mép sông) chịu mô-men uốn cực đại ($M_{max}$) và lực cắt ($Q_{max}$) cao hơn từ 28% đến 32% so với hàng cọc phía sau. Điều này bác bỏ giả thiết tính toán truyền thống vốn xem tải trọng phân bố đều cho toàn bộ các cọc trong đài.
  • Biến thiên hệ số an toàn ổn định tổng thể: Hệ số an toàn trượt sâu ($F_s$) đạt giá trị 1,45 ở giai đoạn đầu nhưng giảm mạnh xuống còn 1,18 khi xuất hiện sự kết hợp bất lợi nhất giữa hiện tượng triều rút nhanh và tải trọng xe thi công $q = 15$ kN/m$^2$ trên mặt bãi.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự gia tăng chuyển vị và nội lực cọc theo thời gian là do quá trình tiêu tán áp lực nước lỗ rỗng trong lớp sét dẻo mềm 1b, gây ra hiện tượng lún cố kết thứ cấp và phát sinh lực ma sát âm bám quanh thân cọc bê tông. Khi đất nền sau kè bị lún, khối đất có xu hướng kéo trượt cọc xuống, đồng thời gia tăng lực đẩy ngang tác dụng lên toàn bộ hệ thống tường chắn.

So sánh với các nghiên cứu tính toán theo mô hình Winkler cổ điển (coi phản lực tác dụng lên cọc là lực tập trung tại đỉnh), phương pháp phần tử hữu hạn thể hiện độ chính xác vượt trội khi chỉ ra rằng mô-men uốn thân cọc xuất hiện giá trị lớn nhất tại độ sâu cách đáy đài khoảng 2,5 đến 3,5 m (tương đương vị trí ngàm mềm trong lớp sét).

Dữ liệu nghiên cứu được biểu diễn rõ ràng qua biểu đồ quan hệ chuyển vị - thời gian cố kết và bảng so sánh nội lực cọc giữa các giai đoạn thi công. Biểu đồ đường cong bao mô-men uốn $M_z$ dọc thân cọc dài 14 m chứng minh rằng nếu không xét đến sự làm việc đồng thời của đất và cọc, kết cấu cọc có thể bị thiết kế thiếu khả năng chịu uốn tới 25%, dẫn đến nguy cơ nứt gãy thân cọc ngầm dưới lòng đất.

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

Dựa trên kết quả mô phỏng và phân tích cơ chế ứng xử công trình kè trên nền đất yếu, luận văn đề xuất 4 nhóm giải pháp kỹ thuật cụ thể:

  • Tối ưu hóa sơ đồ bố trí cọc móng: Áp dụng giải pháp đóng cọc xiên một góc từ 15 đến 20 độ cho hàng cọc phía trước nhằm tăng cường độ cứng chống trượt ngang của toàn hệ thống móng lên 25 - 30%, kết hợp sử dụng cọc bê tông ly tâm dự ứng lực cường độ cao đường kính D400 - D500 mm thay cho cọc vuông đúc tại chỗ để tối ưu hóa khả năng chịu mô-men uốn.
  • Kiểm soát quy trình đắp đất và tốc độ gia tải: Phân chia chiều dày lớp cát đắp san lấp sau lưng tường thành 3 đến 4 đợt thi công (mỗi đợt dày tối đa 0,5 - 0,7 m), duy trì thời gian chờ cố kết tối thiểu 30 ngày giữa các đợt nhằm hạn chế sự gia tăng đột ngột của áp lực nước lỗ rỗng vượt quá ngưỡng an toàn 20 kPa.
  • Xây dựng hệ thống giảm áp và thoát nước ngầm: Bố trí tầng lọc ngược bằng vải địa kỹ thuật kết hợp mạng lưới ống thoát nước $\Phi 90$ mm đục lỗ với khoảng cách 2,0 m một ống dọc theo thân kè, giúp triệt tiêu nhanh áp lực thủy tĩnh khi mực nước triều ngoài sông rút kiệt, giảm tải trọng tác dụng lên bản tường đứng tới 40%.
  • Triển khai hệ thống quan trắc biến dạng định kỳ: Thiết lập hệ thống đo nghiêng Inclinometer trong thân cọc và mốc đo lún bề mặt, thực hiện đo đạc với tần suất 15 ngày một lần trong 6 tháng đầu thi công và 3 tháng một lần trong 2 năm vận hành đầu tiên để chủ động kiểm soát rủi ro chuyển vị vượt mức cho phép.

Về lộ trình thực hiện, công tác khảo sát và hiệu chỉnh thiết kế cần hoàn thành trong 1 đến 2 tháng trước khi triển khai; giai đoạn thi công đắp nền kiểm soát kéo dài 6 đến 9 tháng; công tác quan trắc duy trì liên tục trong 24 đến 36 tháng dưới sự giám sát của Ban Quản lý Dự án Giao thông - Nông nghiệp tỉnh Long An và các đơn vị tư vấn thiết kế chuyên ngành.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thực tiễn cao cho nhiều nhóm chuyên gia và tổ chức:

  • Kỹ sư thiết kế công trình thủy lợi, cảng - đường thủy: Tiếp cận quy trình tính toán tường kè bản góc trên móng cọc chịu lực ngang bằng phần mềm Plaxis 2D, thay thế hiệu quả các công thức tính toán đơn giản hóa có độ rủi ro cao.
  • Chỉ huy trưởng và kỹ sư quản lý chất lượng tại công trường: Nắm bắt các mốc rủi ro địa kỹ thuật khi đắp tải trên nền đất yếu, quy trình hạ cừ ván bê tông dự ứng lực an toàn và giải pháp hạn chế lực ma sát âm tác dụng lên cọc.
  • Ban Quản lý Dự án và cơ quan quản lý nhà nước tại ĐBSCL: Sử dụng luận văn làm tài liệu đối chứng khoa học để thẩm định dự toán, biện pháp thi công và hồ sơ thiết kế kỹ thuật của các dự án kè chống sạt lở ven sông Vàm Cỏ, sông Tiền và sông Hậu.
  • 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 - Xây dựng: Khai thác khung lý thuyết tích hợp giữa mô hình đàn dẻo liên tục, bài toán cố kết thấm Biot và phương pháp dầm trên nền đàn hồi Winkler làm cơ sở phát triển các đề tài nghiên cứu chuyên sâu.

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

  • Tại sao phương pháp phần tử hữu hạn Plaxis 2D lại vượt trội hơn các công thức Rankine hay Coulomb truyền thống?
    Phương pháp phần tử hữu hạn cho phép mô phỏng chính xác sự tương tác phi tuyến giữa đất và kết cấu, xét đến biến dạng đàn dẻo của 4 lớp địa tầng và sự tiêu tán áp lực nước lỗ rỗng theo thời gian thực. Trong khi đó, lý thuyết cổ điển chỉ giả định mặt trượt phẳng và đất đồng nhất, dễ dẫn đến sai số chuyển vị từ 20% đến 35%.

  • Hiện tượng ma sát âm trên thân cọc phát sinh như thế nào và cách xử lý ra sao?
    Khi lớp cát san lấp đè lên tầng sét mềm 1b, đất bị lún cố kết nhanh hơn chuyển vị của cọc, tạo ra lực kéo ma sát âm dọc thân cọc. Để khắc phục, cần kéo dài thời gian chờ cố kết từng lớp đắp, sử dụng lớp lót vải địa kỹ thuật phân cách hoặc quét lớp sơn giảm ma sát bitum quanh đoạn cọc đi qua lớp đất yếu.

  • Việc đóng cọc xiên 15 đến 20 độ mang lại lợi ích cơ học cụ thể gì cho tường kè?
    Cọc xiên tạo ra thành phần lực chống cắt nằm ngang trực tiếp đối kháng với áp lực đất chủ động và áp lực thủy tĩnh sau lưng tường. Góc xiên 15 đến 20 độ giúp tăng khả năng chịu tải ngang của toàn móng lên 25 - 30% và giảm mô-men uốn cực đại tại đầu cọc xuống khoảng 15%.

  • Chế độ nhật triều tại khu vực Long An ảnh hưởng như thế nào đến độ ổn định kè?
    Với 2 lần nước lớn và nước ròng trong ngày, chênh lệch mực nước sông từ +2,30 m xuống +0,40 m sinh ra dòng thấm ngược từ khối đắp ra sông. Áp lực thấm kết hợp với trọng lượng đất bão hòa nước làm hệ số an toàn trượt mái dốc giảm từ 1,45 xuống mức nguy hiểm 1,18 nếu không có hệ thống tầng lọc thoát nước ngầm.

  • Hệ số phản lực nền ngang K trong mô hình Winkler được xác định như thế nào trong nghiên cứu?
    Hệ số $K$ phản ánh tính đàn hồi của đất nền, được xác định dựa trên độ sệt và chỉ số dẻo của đất theo tiêu chuẩn TCVN 10304:2014. Đối với lớp sét dẻo mềm, hệ số $K$ dao động từ 650 đến 2.500 kN/m$^4$, trong khi lớp sét nửa cứng đạt từ 5.000 đến 8.000 kN/m$^4$, quyết định trực tiếp đến độ võng và góc xoay của thân cọc.

Kết luận

  • Luận văn đã làm sáng tỏ cơ chế tương tác đất - cọc - tường kè trên nền đất yếu tại khu vực sông Vàm Cỏ Tây, tỉnh Long An thông qua phương pháp phần tử hữu hạn 2D.
  • Chứng minh sự phân bố phi tuyến của áp lực đất sau lưng tường với độ tăng áp lực cục bộ 22% tại ngàm chân đài và mức chuyển vị đỉnh kè tích lũy đạt 3,8 cm sau thời gian cố kết.
  • Xác định rõ sự phân bổ mô-men uốn lệch nề giữa các hàng cọc, trong đó cọc mép bờ sông chịu tải trọng lớn hơn 28 - 32% so với hàng cọc phía trong.
  • Đề xuất thành công gói giải pháp kỹ thuật bao gồm: bố trí cọc xiên 15 - 20 độ, phân đợt đắp cát gia tải kiểm soát lún và lắp đặt hệ thống tầng lọc tiêu thoát nước giảm áp lực thủy tĩnh.
  • Khẳng định tính cấp thiết của việc áp dụng mô hình phân tích số thay thế phương pháp tính toán tĩnh truyền thống nhằm nâng cao hệ số ổn định công trình $F_s \ge 1,30$ và tiết kiệm 15 - 20% chi phí xử lý sự cố công trình ven sông.

Đóng góp lớn nhất của đề tài là xây dựng quy trình tính toán - thiết kế chuẩn xác có xét đến yếu tố thời gian và trình tự thi công thực tế cho công trình kè bảo vệ bờ tại Đồng bằng sông Cửu Long. Trong giai đoạn tiếp theo từ 12 đến 24 tháng tới, các đơn vị quản lý hạ tầng cần đẩy mạnh tích hợp công nghệ quan trắc tự động và mở rộng mô hình tính toán sang bài toán 3D để hoàn thiện tiêu chuẩn thiết kế bờ kè chống biến đổi khí hậu trong toàn khu vực.