Tổng quan nghiên cứu

Khu vực Đồng bằng sông Cửu Long nói chung và thành phố Cần Thơ nói riêng có cấu trúc địa chất đặc trưng bởi tầng đất sét yếu bồi tích trẻ (Holocen) dày từ 15m đến hơn 30m, với độ ẩm tự nhiên cao vượt mức 50% đến 85% và sức kháng cắt rất thấp. Trong thiết kế nền móng, ổn định mái dốc và công trình kè bảo vệ bờ sông, module biến dạng (E) là thông số cơ lý đầu vào cốt lõi quyết định độ chính xác của bài toán dự báo độ lún và chuyển vị công trình. Thí nghiệm nén ba trục cố kết không thoát nước (CU) được đánh giá là phương pháp mô phỏng gần sát nhất trạng thái làm việc và ứng xử chịu lực tự nhiên của đất nền. Tuy nhiên, do chi phí cao và quy trình thí nghiệm kéo dài từ 5 đến 15 ngày cho mỗi tổ mẫu, các dự án quy mô vừa và nhỏ tại địa phương thường thiếu hụt dữ liệu nén ba trục chi tiết.

Nghiên cứu của tác giả Đỗ Thuận Phát thực hiện năm 2014 tập trung giải quyết bài toán cấp thiết này với mục tiêu thiết lập mối tương quan thực nghiệm giữa module biến dạng cát tuyến không thoát nước ($E_{u50}$) và ứng suất hữu hiệu theo phương đứng do trọng lượng bản thân ($\sigma'_v$) cho 3 loại đất sét phổ biến tại Cần Thơ. Đề tài được giới hạn trong phạm vi xử lý số liệu địa kỹ thuật thực tế trên địa bàn thành phố và ứng dụng mô phỏng số cho công trình kè sông Cần Thơ. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn cao, giúp giảm thiểu sai số tính toán chuyển vị ngang tường kè từ 18% đến 24%, đồng thời tiết kiệm 30% chi phí khảo sát bổ sung cho các đơn vị tư vấn thiết kế tại khu vực Tây Nam Bộ.

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 kinh điển và lý thuyết phần tử hữu hạn hiện đại:

  • Thuyết phá hoại Mohr - Coulomb: Mô tả trạng thái cân bằng giới hạn của phân tố đất khi chịu tổ hợp ứng suất chính, xác định sức kháng cắt qua phương trình $\tau_{gh} = \sigma \tan\varphi + c$. Vòng tròn ứng suất Mohr biểu diễn sự gia tăng ứng suất lệch $(\sigma_1 - \sigma_3)$ đến điểm phá hoại trong buồng nén ba trục.
  • Lý thuyết cố kết thấm của Terzaghi và lý thuyết bán không gian đàn hồi: Phân tích quá trình tiêu tán áp lực nước lỗ rỗng thặng dư và biến dạng tức thời của nền đất dưới tải trọng công trình.
  • Các khái niệm then chốt:
    • Module biến dạng cát tuyến không thoát nước ($E_{u50}$): Đoạn dốc nối từ gốc tọa độ đến điểm ứng với 50% ứng suất lệch cực hạn trên biểu đồ quan hệ ứng suất - biến dạng.
    • Ứng suất hữu hiệu thẳng đứng ($\sigma'_v$): Thành phần ứng suất do trọng lượng bản thân các lớp đất bên trên truyền qua cốt hạt khoáng, tính bằng $\sigma_v - u$.
    • Tương quan sức chống cắt không thoát nước ($S_u$): Khung tham chiếu thực nghiệm của Duncan & Buchignani (1976) với tỷ số $E_u/S_u$ dao động từ 150 đến 1500 tùy thuộc chỉ số dẻo $I_p$ và hệ số quá cố kết OCR.

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

Nghiên cứu áp dụng quy trình thực nghiệm kết hợp phân tích thống kê và mô phỏng số với các tiêu chí chặt chẽ:

  • Cỡ mẫu và phương pháp chọn mẫu: Bộ dữ liệu bao gồm hơn 100 mẫu đất sét nguyên dạng được lấy từ các hố khoan khảo sát địa chất công trình tại Cần Thơ. Mẫu được chọn theo phương pháp phân tầng đại diện cho 3 trạng thái cơ lý: sét trạng thái nhão (độ sệt $I_L > 1,0$), sét dẻo nhão ($0,75 < I_L \le 1,0$) và sét dẻo cứng ($0,25 < I_L \le 0,5$). Mẫu thí nghiệm có kích thước tiêu chuẩn hình trụ tròn $38 \times 76$ mm, bọc màng cao su và cố kết đẳng hướng trước khi nén không thoát nước với tốc độ biến dạng kiểm soát $1% - 2%/\text{phút}$.
  • Lý do lựa chọn phương pháp phân tích: Phương pháp hồi quy tuyến tính được lựa chọn để chuẩn hóa quan hệ $E_{u50} = f(\sigma'_v)$ giúp kỹ sư dễ dàng tra cứu. Phương pháp so sánh đối chiếu với tiêu chuẩn TCVN 9362:2012 và nghiên cứu quốc tế của Mitchell & Gardner (1975), Ladd & Edgers (1972) nhằm xác thực độ tin cậy. Cuối cùng, phương pháp mô phỏng số bằng phần mềm Plaxis 2D được sử dụng để kiểm chứng định lượng ứng xử biến dạng của kè sông Cần Thơ trong 4 kịch bản gán thông số độ cứng.
  • Timeline nghiên cứu: Toàn bộ công tác thu thập mẫu, thí nghiệm phòng và mô phỏng số hoàn thành trong giai đoạn từ tháng 02/2014 đến tháng 06/2014.

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

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

Phân tích xử lý dữ liệu từ hơn 100 mẫu nén ba trục CU đã mang lại các phát hiện khoa học mang tính ứng dụng cao:

  • Xác lập hàm tương quan tuyến tính cho 3 phân tầng sét:
    • Lớp sét trạng thái nhão: Giá trị $E_{u50}$ dao động từ 1.250 kPa đến 3.400 kPa tương ứng dải ứng suất $\sigma'_v$ từ 20 kPa đến 75 kPa, tỷ lệ gia tăng theo chiều sâu đạt trung bình 28 kPa/m.
    • Lớp sét trạng thái dẻo nhão: Giá trị $E_{u50}$ phân bố trong khoảng 3.200 kPa đến 7.600 kPa khi $\sigma'_v$ tăng từ 40 kPa đến 115 kPa.
    • Lớp sét trạng thái dẻo cứng: Độ cứng tăng vọt rõ rệt, $E_{u50}$ biến thiên từ 8.500 kPa đến trên 21.000 kPa khi ứng suất $\sigma'_v$ vượt 120 kPa.
  • Quy luật tỷ số $E_{u50}/S_u$: Đối với đất sét yếu bồi tích Cần Thơ có hệ số quá cố kết $\text{OCR} \approx 1,0 - 1,2$, tỷ số $E_{u50}/S_u$ đạt giá trị trung bình từ 250 đến 450, hoàn toàn tương thích với các nghiên cứu của Bjerrum (250 đến 500) và Ladd & Edgers (280 đến 420).
  • Tác động rõ rệt trong bài toán mô phỏng chuyển vị: Kết quả phân tích phần tử hữu hạn 4 trường hợp mô hình kè sông Cần Thơ chỉ ra rằng: khi áp dụng module biến dạng $E_{u50}$ thay đổi theo chiều sâu (tăng theo $\sigma'_v$), chuyển vị ngang lớn nhất của tường kè giảm 18,5% đến 24,2% và độ lún mặt đất sau lưng kè giảm từ 14,3% đến 21,0% so với trường hợp gán module không đổi theo chiều sâu.

Thảo luận kết quả

Sự gia tăng tuyến tính của $E_{u50}$ theo $\sigma'_v$ phản ánh bản chất nén chặt tự nhiên của đất nền dưới áp lực địa tĩnh: chiều sâu càng lớn, liên kết giữa các hạt khoáng sét càng chặt chẽ, hệ số rỗng giảm dần khiến khả năng kháng biến dạng trượt tăng lên.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các biểu đồ phân tán (scatter plots) kết hợp đường xu hướng tuyến tính có hệ số tương quan $R^2 > 0,82$, cùng các bảng tổng hợp so sánh sai số nội lực và chuyển vị giữa 4 trường hợp mô phỏng. Khi so sánh với phương pháp quy đổi module $E_0$ từ thí nghiệm nén cố kết oedometer theo hệ số $m_k$ của TCVN 9362:2012 hoặc suy từ thí nghiệm xuyên tĩnh CPT ($E = 7 \cdot q_c$), việc áp dụng trực tiếp $E_{u50}$ từ thí nghiệm ba trục CU giúp mô tả chính xác độ cứng chống cắt không thoát nước thực tế, loại bỏ hiện tượng đánh giá quá non độ lún tức thời và chuyển vị ngang trong các bài toán hố đào và bờ kè.

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

Dựa trên kết quả thực nghiệm và phân tích mô phỏng, tác giả đưa ra 4 khuyến nghị hành động cụ thể:

  • Ứng dụng trực tiếp hàm tương quan $E_{u50} = f(\sigma'_v)$ vào thiết kế: Đơn vị tư vấn thiết kế cần áp dụng các phương trình tương quan đã xác lập để gán thông số đầu vào cho mô hình Mohr - Coulomb trong Plaxis và GeoStudio. Mục tiêu giảm sai số dự báo chuyển vị xuống dưới 10%, thực hiện ngay trong giai đoạn thiết kế cơ sở và kỹ thuật (thời gian từ 15 đến 30 ngày) do kỹ sư địa kỹ thuật phụ trách.
  • Kích hoạt thông số tăng độ cứng theo chiều sâu ($E_{\text{increment}}$): Trong mô phỏng phần tử hữu hạn, bắt buộc kích hoạt tham số biến thiên module theo chiều sâu thay vì gán giá trị cố định duy nhất cho toàn bộ tầng sét dày trên 10m. Giải pháp này giúp tối ưu hóa kết cấu tường kè, giảm từ 8% đến 15% khối lượng cốt thép cọc ván dự ứng lực, do đơn vị thẩm tra thiết kế kiểm soát.
  • Chuẩn hóa tỷ lệ thí nghiệm nén ba trục trong khảo sát: Chủ đầu tư các dự án công trình cấp II trở lên tại khu vực Cần Thơ cần quy định tỷ lệ mẫu thí nghiệm nén ba trục CU đạt tối thiểu 10% đến 15% trên tổng số mẫu cơ lý lấy từ hiện trường. Quy định áp dụng trong chu kỳ khảo sát 2024–2026 do các nhà thầu khảo sát địa chất thực thi.
  • Xây dựng hệ thống cơ sở dữ liệu địa kỹ thuật số: Sở Xây dựng thành phố Cần Thơ phối hợp với các trường đại học triển khai bản đồ số hóa tương quan địa kỹ thuật tích hợp dữ liệu ba trục CU, SPT và CPT. Mục tiêu rút ngắn 25% thời gian tra cứu và tính toán ban đầu, hoàn thiện trong lộ trình 12 tháng.

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

  • Kỹ sư thiết kế địa kỹ thuật và nền móng: Khai thác ngân hàng dữ liệu tương quan $E_{u50} - \sigma'_v$ đã được kiểm chứng để tính toán ổn định bờ kè, đê bao ngăn triều và móng cọc công trình ven sông mà không cần chờ đợi kết quả thí nghiệm kéo dài.
  • Chuyên viên khảo sát địa chất công trình: Nắm vững phương pháp chuẩn bị mẫu $38 \times 76$ mm, quy trình vận hành buồng nén ba trục CU đối xứng trục và kỹ thuật trích xuất module cát tuyến $E_{u50}$ chuẩn xác từ biểu đồ ứng suất lệch.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Địa kỹ thuật: Sử dụng làm tài liệu tham khảo tiêu biểu về phương pháp luận kết hợp giữa nghiên cứu thực nghiệm cơ học đất, phân tích thống kê toán học và kiểm chứng mô phỏng số bằng Plaxis 2D.
  • Đơn vị thẩm tra và chủ đầu tư dự án hạ tầng Tây Nam Bộ: Sử dụng các chỉ số đối chứng trong luận văn làm cơ sở kỹ thuật độc lập để đánh giá tính chuẩn xác của báo cáo khảo sát và phương án gia cố nền đất yếu.

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

Module biến dạng $E_{u50}$ từ thí nghiệm ba trục CU khác biệt thế nào so với module $E_0$ từ thí nghiệm nén cố kết?
Module $E_{u50}$ phản ánh độ cứng cắt không thoát nước ở trạng thái ứng suất đa trục có nở hông thực tế của đất nền. Trong khi đó, module $E_0$ từ thí nghiệm nén cố kết oedometer là biến dạng một chiều bị chặn hông hoàn toàn, thường phản ánh độ nén lún thể tích dài hạn thay vì ứng xử biến dạng trượt ngang của tường kè.

Tại sao cần phải xét sự thay đổi của module biến dạng theo chiều sâu khi tính toán kè sông?
Ứng suất bản thân $\sigma'_v$ tăng dần theo độ sâu làm tăng mật độ hạt và độ cứng của đất sét. Nếu giả định module không đổi, mô hình sẽ tính toán sai lệch lớn: chuyển vị đỉnh kè sông Cần Thơ có thể bị khuếch đại từ 18% đến 24%, dẫn đến việc thiết kế thừa kích thước kết cấu chịu lực gây lãng phí kinh phí.

Tỷ số $E_{u50}/S_u$ của đất sét khu vực Cần Thơ có phù hợp với các nghiên cứu quốc tế không?
Kết quả phân tích hơn 100 mẫu cho thấy tỷ số $E_{u50}/S_u$ của đất sét Cần Thơ tập trung từ 250 đến 450. Giá trị này hoàn toàn nằm trong dải tương quan thực nghiệm đã được công nhận toàn cầu của Bjerrum (250 đến 500) và Ladd & Edgers (280 đến 420) đối với đất sét cố kết bình thường có chỉ số dẻo $I_p > 30%$.

Các hàm tương quan trong đề tài có áp dụng được cho các tỉnh lân cận ở Đồng bằng sông Cửu Long không?
Công thức hoàn toàn có thể áp dụng cho các khu vực lân cận có cùng nguồn gốc trầm tích sông - biển Holocen như Hậu Giang, Vĩnh Long và Sóc Trăng với mức sai số dự báo dưới 12%. Tuy nhiên, đối với các tầng sét cổ Pleistocene có độ cố kết cao hơn, cần tiến hành thí nghiệm hiệu chỉnh cục bộ.

Khi nào kỹ sư bắt buộc phải sử dụng thí nghiệm nén ba trục CU thay cho thí nghiệm cắt phẳng trực tiếp?
Thí nghiệm ba trục CU là bắt buộc đối với các công trình ven sông quy mô lớn, hố đào sâu hoặc mái dốc phức tạp. Phương pháp này cho phép kiểm soát áp lực nước lỗ rỗng và hướng ứng suất chính, khắc phục hoàn toàn nhược điểm áp đặt mặt trượt cưỡng bức của thí nghiệm cắt phẳng.

Kết luận

  • Luận văn đã xác lập thành công hệ phương trình tương quan tuyến tính giữa module biến dạng $E_{u50}$ và ứng suất hữu hiệu $\sigma'_v$ cho 3 nhóm đất sét đặc trưng tại Cần Thơ dựa trên 100 mẫu thí nghiệm nén ba trục CU.
  • Khẳng định giá trị thực nghiệm tỷ số $E_{u50}/S_u$ dao động từ 250 đến 450, mang lại độ tin cậy cao theo các tiêu chuẩn quốc tế và khu vực.
  • Chứng minh rõ nét hiệu quả của việc xét biến thiên module theo chiều sâu trong phần mềm Plaxis, giúp giảm sai số chuyển vị đỉnh kè từ 18,5% đến 24,2% và độ lún từ 14,3% đến 21,0%.
  • Đóng góp nguồn dữ liệu tham chiếu quan trọng cho ngành xây dựng Cần Thơ, giải quyết bài toán thiếu hụt thông số nén ba trục trong các dự án vừa và nhỏ.
  • Định hình quy trình tích hợp giữa cơ học đất thực nghiệm và mô phỏng số trong phân tích ổn định công trình ven sông.

Đề tài là công trình nghiên cứu giá trị, đóng góp giải pháp kỹ thuật thiết thực cho việc xử lý nền đất yếu tại Đồng bằng sông Cửu Long. Trong vòng 6 đến 12 tháng tới, các cơ quan chuyên môn nên tiếp tục mở rộng ngân hàng dữ liệu thực nghiệm cho các vùng phụ cận. Hãy tham khảo và ứng dụng ngay các hàm tương quan này vào hồ sơ thiết kế công trình của bạn để tối ưu hóa chi phí và đảm bảo an toàn kết cấu!