Tổng quan nghiên cứu

Khu vực phía Nam Thành phố Hồ Chí Minh, đặc biệt là Quận 7, Nhà Bè và Bình Chánh, đang trải qua tốc độ đô thị hóa nhanh chóng với hơn 60% diện tích xây dựng nằm trên các tầng trầm tích bùn sét bão hòa nước có chiều dày biến thiên từ 15m đến hơn 35m. Vấn đề trọng tâm trong công tác thiết kế địa kỹ thuật tại đây là tính toán chính xác độ lún và tốc độ cố kết của nền đất yếu khi áp dụng giải pháp gia tải trước. Hiện nay, các phòng thí nghiệm cơ học đất sử dụng nhiều loại dao vòng nén cố kết có kích thước khác nhau theo tiêu chuẩn thương mại như đường kính 50.8mm, 61.8mm và 79.8mm. Sự sai khác về đường kính và chiều cao mẫu nén dẫn đến hiện tượng ma sát biên và thay đổi mức độ đại diện vi cấu trúc đất, gây ra sai số dự báo lún từ 15% đến 30% so với diễn biến thực tế tại công trình.

Nghiên cứu tập trung giải quyết bài toán cốt lõi: phân tích định lượng ảnh hưởng của kích thước mẫu nén trong thí nghiệm nén cố kết một chiều đến các thông số biến dạng của đất dính yếu bão hòa nước, từ đó ứng dụng vào bài toán dự báo độ lún công trình gia tải trước tại khu vực dự án Riviera Point, Phường Tân Phú, Quận 7, Thành phố Hồ Chí Minh trong giai đoạn khảo sát và thi công năm 2015. Mục tiêu cụ thể là xác định tổ hợp kích thước mẫu nén tối ưu giúp giảm thiểu sai số giữa kết quả giải tích, mô phỏng số trên phần mềm Plaxis 2D và số liệu quan trắc lún thực tế ngoài hiện trường.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc chuẩn hóa quy trình chế bị mẫu nén trong phòng, giúp các kỹ sư tư vấn kiểm soát độ lún dư sau xử lý dưới ngưỡng cho phép 10cm, đồng thời rút ngắn thời gian chờ cố kết từ 12 tháng xuống còn khoảng 8 đến 9 tháng, mang lại hiệu quả kinh tế vượt trội cho các dự án hạ tầng đô thị.

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 dựa trên nền tảng lý thuyết cố kết thấm một chiều kinh điển của Karl Terzaghi (1925), mô tả quá trình tiêu tán áp lực nước lỗ rỗng thặng dư và gia tăng ứng suất hữu hiệu theo thời gian qua phương trình vi phân đạo hàm riêng bậc hai:

$$\frac{\partial u}{\partial t} = C_v \frac{\partial^2 u}{\partial z^2}$$

Trong đó, hệ số cố kết thẳng đứng $C_v$ đóng vai trò quyết định đến tốc độ lún của nền đất dưới các cấp tải trọng từ 12.5 kPa đến 400 kPa. Song song với đó, nghiên cứu tích hợp phương pháp phân tích ổn định trượt sâu Bishop đơn giản hóa và lý thuyết cân bằng giới hạn Prandtl - Terzaghi thông qua hệ số sức chịu tải $N_c$ nhằm kiểm soát khả năng trượt trồi của khối đắp gia tải trước khi chiều cao đắp vượt quá 3.5m.

Khung phân tích tập trung vào bốn khái niệm cơ học đất then chốt:

  1. Áp lực tiền cố kết ($P_c$): Ngưỡng ứng suất lịch sử lớn nhất mà đất từng chịu đựng, phân định trạng thái cố kết bình thường ($OCR = 1$) và quá cố kết ($OCR > 1$).
  2. Chỉ số nén ($C_c$) và chỉ số nở ($C_s$): Độ dốc của đường cong nén lún trên hệ tọa độ $e - \log(p)$, phản ánh mức độ giảm thể tích lỗ rỗng khi gia tải.
  3. Hệ số nén thể tích ($m_v$): Biến dạng thể tích trên một đơn vị gia tăng ứng suất hữu hiệu.
  4. Hệ số nén thứ cấp ($C_\alpha$): Biến dạng từ biến do sự tái sắp xếp cấu trúc hạt đất sau khi áp lực nước lỗ rỗng thặng dư đã tiêu tán hoàn toàn.

Quy trình thí nghiệm và đánh giá tuân thủ nghiêm ngặt các tiêu chuẩn kỹ thuật hiện hành, đặc biệt là tiêu chuẩn quốc gia TCVN 4200-2012 về phương pháp xác định tính nén lún trong phòng thí nghiệm.

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

Nguồn dữ liệu của nghiên cứu được thu thập từ 25 lỗ khoan khảo sát địa chất công trình tại dự án Riviera Point (Quận 7, TP.HCM), với chiều sâu khoan khảo sát tối đa đạt 50m. Mẫu đất nghiên cứu là đất dính yếu thuộc lớp sét pha bùn và bùn sét màu xám đen, trạng thái chảy đến dẻo chảy, phân bố liên tục từ độ sâu 1.5m đến 22.0m. Phương pháp lấy mẫu nguyên dạng sử dụng ống mẫu thành mỏng chuyên dụng có đường kính ngoài 100mm, đường kính trong 95mm và chiều dài 2000mm, được bịt kín bằng màng cao su và sáp parafin nhằm duy trì độ ẩm tự nhiên $W = 68.5%$ và dung trọng tự nhiên $\gamma = 15.8\text{ kN/m}^3$.

Phương pháp phân tích thực nghiệm bao gồm:

  • Thí nghiệm nén 3 trục cố kết không thoát nước (CU) có đo áp lực nước lỗ rỗng ở 3 cấp áp lực buồng: 50 kPa, 100 kPa và 150 kPa để xác định chính xác góc ma sát trong $\varphi = 15^\circ$ và lực dính hiệu quả $c' = 12.5\text{ kPa}$.
  • Thí nghiệm nén cố kết oedometer đa cấp với 5 tổ hợp kích thước dao vòng: diện tích tiết diện $20\text{ cm}^2$ (đường kính $5.08\text{ cm}$), $30\text{ cm}^2$ (đường kính $6.18\text{ cm}$), $50\text{ cm}^2$ (đường kính $7.98\text{ cm}$) kết hợp các chiều cao mẫu biến thiên $2.0\text{ cm}$, $3.0\text{ cm}$ và $4.0\text{ cm}$.
  • Mô phỏng phần tử hữu hạn 2D trên phần mềm chuyên ngành Plaxis 2D kết hợp mô-đun Slope/W của GeoStudio để kiểm tra hệ số an toàn ổn định mái dốc $K_{at} \ge 1.35$.
  • Đối soát hồi quy với chuỗi số liệu quan trắc lún thực tế bằng bàn đo lún bề mặt và mốc định vị GPS trong thời gian theo dõi liên tục 180 ngày (từ tháng 1/2015 đến tháng 6/2015).

Lý do lựa chọn phương pháp tích hợp này là nhằm loại bỏ các giả định lý thuyết đơn thuần, tạo cầu nối trực tiếp giữa thí nghiệm vi mô trong phòng và ứng xử biến dạng vĩ mô ngoài hiện trường.

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

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

Nghiên cứu đã chỉ ra những phát hiện thực nghiệm quan trọng về mối tương quan giữa kích thước hình học mẫu nén và các đặc trưng biến dạng cơ học của đất bùn sét:

  1. Ảnh hưởng rõ rệt của đường kính mẫu nén đến chỉ số nén ($C_c$): Khi giữ nguyên chiều cao mẫu $2.0\text{ cm}$ và gia tăng đường kính từ $5.08\text{ cm}$ ($20\text{ cm}^2$) lên $6.18\text{ cm}$ ($30\text{ cm}^2$) và $7.98\text{ cm}$ ($50\text{ cm}^2$), tỷ lệ biến dạng tương đối tăng từ 12.4% lên 18.2% dưới cùng cấp áp lực nén 200 kPa. Giá trị chỉ số nén $C_c$ tăng từ 0.76 ở mẫu nhỏ lên 0.92 ở mẫu đường kính lớn (tăng xấp xỉ 21.0%), trong khi áp lực tiền cố kết $P_c$ dao động trong khoảng hẹp từ 0.35 đến 0.50 kg/cm2.

  2. Quy luật biến thiên của hệ số cố kết ($C_v$) theo kích thước dao vòng: Mẫu nén đường kính nhỏ ($5.08\text{ cm}$) cho giá trị hệ số cố kết $C_v$ cao hơn từ 15% đến 25% so với mẫu đường kính lớn ($7.98\text{ cm}$) ở các cấp tải trọng ban đầu dưới 100 kPa. Nguyên nhân là do ở mẫu nhỏ, nước lỗ rỗng có xu hướng thoát nhanh dọc theo khe hở tiếp xúc giữa đất và thành dao vòng kim loại.

  3. Tác động của chiều cao mẫu nén ($h = 2\text{ cm}, 3\text{ cm}, 4\text{ cm}$): Khi tăng chiều cao mẫu nén lên $4.0\text{ cm}$ với đường kính $7.98\text{ cm}$, hệ số nén thể tích $m_v$ thể hiện tính liên tục và ổn định cao hơn, giảm hiện tượng nhiễu loạn cục bộ do ma sát hai đầu mặt đá thấm. Mẫu cao $4.0\text{ cm}$ cho giá trị hệ số nén thứ cấp $C_\alpha$ ổn định quanh mức 0.035, phản ánh chính xác bản chất từ biến dẻo của đất yếu hữu cơ.

  4. Độ chính xác vượt trội của mẫu $d \times h = 7.98\text{ cm} \times 4.0\text{ cm}$ trong mô hình dự báo: Kết quả tính toán giải tích và mô phỏng số Plaxis 2D sử dụng thông số từ mẫu $7.98\text{ cm} \times 4.0\text{ cm}$ cho độ lún ổn định cuối cùng đạt 132.5cm, sai số chỉ 3.8% so với độ lún quan trắc hiện trường (137.8cm). Ngược lại, việc sử dụng thông số từ mẫu tiêu chuẩn kích thước nhỏ ($5.08\text{ cm} \times 2.0\text{ cm}$) dẫn đến độ lún dự báo chỉ đạt 108.2cm, gây thiếu hụt tới 21.5% độ lún thực tế.

Thảo luận kết quả

Sự khác biệt về kết quả nén lún giữa các kích thước mẫu bắt nguồn từ hai cơ chế vật lý cơ bản: hiệu ứng ma sát thành và tính đại diện cấu trúc trầm tích.

Về mặt ma sát thành, lực dính và ma sát tiếp xúc giữa đất sét dẻo và thành trong dao vòng làm tiêu hao từ 10% đến 18% ứng suất thẳng đứng truyền vào mẫu đất nhỏ ($5.08\text{ cm}$), khiến áp lực nén thực tế tác dụng lên khung hạt đất bị suy giảm đáng kể. Khi tăng đường kính lên $7.98\text{ cm}$, tỷ số diện tích xung quanh trên diện tích đáy ($4/d$) giảm 36.3%, giúp giảm thiểu tổn thất ứng suất do ma sát biên.

Về mặt cấu trúc, đất yếu tại khu vực Quận 7 chứa nhiều vi thấu kính cát mịn và tàn tích hữu cơ không đồng nhất. Mẫu nén thể tích lớn ($V \approx 200\text{ cm}^3$ ở mẫu $7.98\text{ cm} \times 4\text{ cm}$) bảo toàn được đầy đủ các đặc trưng cấu trúc vĩ mô và vi mô tự nhiên hơn so với mẫu nhỏ ($V \approx 40.5\text{ cm}^3$). Phát hiện này hoàn toàn tương đồng với các công bố quốc tế của Brendan C. O'Kelly (2009) trên đất than bùn Ireland và S. Bharati (2011) trong khảo sát địa kỹ thuật công trình biển.

Trong báo cáo kỹ thuật, dữ liệu biến dạng nén lún có thể được trực quan hóa tối ưu thông qua hệ thống biểu đồ đường cong nén lún $e - \log(p)$, đồ thị quan hệ $C_v - \log(p)$ qua các cấp áp lực từ 12.5 kPa đến 400 kPa, và biểu đồ chữ S biểu diễn tốc độ lún theo thời gian so sánh giữa đường quan trắc mốc lún GPS với đường mô phỏng Plaxis 2D.

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

Dựa trên các kết luận thực nghiệm vững chắc, nghiên cứu đưa ra 4 giải pháp và kiến nghị hành động cụ thể nhằm nâng cao chất lượng thiết kế xử lý nền đất yếu:

  1. Chuẩn hóa kích thước dao vòng trong phòng thí nghiệm địa kỹ thuật: Khuyến nghị các phòng thí nghiệm LAS-XD trên toàn quốc ưu tiên áp dụng dao vòng nén cố kết có kích thước đường kính tối thiểu $7.98\text{ cm}$ và chiều cao $4.0\text{ cm}$ khi thực hiện thí nghiệm cho các dự án xây dựng trên nền bùn sét yếu ven sông, ven biển; đặt mục tiêu giảm sai số dự báo lún xuống dưới 5% trong vòng 12 tháng áp dụng.
  2. Hiệu chỉnh hệ số ma sát thành cho thiết bị nén thương mại hiện hữu: Đề xuất các đơn vị tư vấn thiết kế áp dụng hệ số hiệu chỉnh tăng tải trọng từ 1.08 đến 1.15 khi phải sử dụng thông số nén từ dao vòng nhỏ ($5.08\text{ cm}$ hoặc $6.18\text{ cm}$) nhằm bù đắp lượng ứng suất tiêu hao do ma sát thành, hoàn thành việc xây dựng bảng tra hệ số trong quý 3.
  3. Ứng dụng quy trình mô phỏng số Plaxis 2D kết hợp phương pháp Asaoka: Yêu cầu các kỹ sư thiết kế tích hợp mô hình Plaxis 2D với phương pháp ngoại suy Asaoka dựa trên chuỗi quan trắc 3 ngày/lần trong 60 ngày đầu đắp tải, đảm bảo kiểm soát chặt chẽ hệ số ổn định trượt sâu $K_{at} \ge 1.40$ và áp lực nước lỗ rỗng thặng dư dưới 25 kPa.
  4. Cập nhật tiêu chuẩn thiết kế và thi công nền móng: Kiến nghị Viện Khoa học Công nghệ Xây dựng (IBST) và Bộ Xây dựng xem xét bổ sung các quy định cụ thể về hiệu ứng kích thước mẫu nén cho đất dính bão hòa nước vào dự thảo sửa đổi tiêu chuẩn TCVN 4200 trong giai đoạn 2026-2027.

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

Luận văn thạc sĩ này là tài liệu tham khảo có giá trị học thuật và thực tiễn cao cho 4 nhóm đối tượng trọng tâm sau:

  1. Kỹ sư tư vấn thiết kế địa kỹ thuật và công trình ngầm: Nắm bắt cơ sở lựa chọn thông số đầu vào ($C_c, C_s, C_v, E_{oed}$) chính xác cho các bài toán phân tích phần tử hữu hạn; ứng dụng trực tiếp vào việc thiết kế đắp gia tải trước kết hợp bấc thấm (PVD) hoặc hút chân không cho các khu đô thị và đường cao tốc.
  2. Cán bộ quản lý và kỹ thuật viên phòng thí nghiệm LAS-XD: Nắm vững kỹ thuật gọt tỉa mẫu bùn sét kích thước lớn ($7.98\text{ cm} \times 4\text{ cm}$) mà không gây xáo động kết cấu; tối ưu hóa quy trình bảo dưỡng bão hòa mẫu và kiểm soát tốc độ tăng tải từ 12.5 kPa đến 400 kPa.
  3. Học viên cao học, nghiên cứu sinh và giảng viên ngành Kỹ thuật Xây dựng Công trình Ngầm - Địa kỹ thuật: Sử dụng làm tài liệu giảng dạy chuyên sâu về môn học Cơ học đất nâng cao, Mô hình hóa địa kỹ thuật và Phân tích tương tác đất - kết cấu; làm tiền đề mở rộng các đề tài nghiên cứu về từ biến và phá hoại trượt dẻo.
  4. Ban quản lý dự án và chủ đầu tư hạ tầng giao thông, đô thị: Có công cụ khoa học để thẩm tra, phản biện hồ sơ thiết kế xử lý nền đất yếu; dự báo chính xác khối lượng đất đắp bù lún (thường chiếm 15% đến 25% tổng khối lượng đắp) và tối ưu hóa tiến độ hoàn thiện mặt đường.

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

Kích thước mẫu nén ảnh hưởng như thế nào đến chỉ số nén $C_c$ của bùn sét?

Khi tăng đường kính mẫu từ $5.08\text{ cm}$ lên $7.98\text{ cm}$, diện tích chịu nén tăng gấp 2.5 lần giúp giảm thiểu tỷ lệ lực cản ma sát tiếp xúc ở thành dao vòng từ khoảng 18% xuống dưới 7%. Do đó, năng lượng gia tải truyền trực tiếp vào khung hạt đất lớn hơn, làm tăng chỉ số nén $C_c$ từ 0.76 lên 0.92, phản ánh đúng khả năng nén lún thực tế của đất.

Tại sao mẫu nén có kích thước $d \times h = 7.98\text{ cm} \times 4.0\text{ cm}$ lại cho kết quả dự báo lún chính xác nhất?

Kích thước $7.98\text{ cm} \times 4.0\text{ cm}$ tạo ra thể tích mẫu đủ lớn ($V \approx 200\text{ cm}^3$) để chứa đựng đầy đủ các vi thấu kính cát và cấu trúc hữu cơ tự nhiên, đồng thời duy trì tỷ lệ $h/d \approx 0.5$ tối ưu. Khi đưa thông số này vào Plaxis 2D, sai số độ lún so với quan trắc thực tế tại Quận 7 chỉ là 3.8% (132.5cm so với 137.8cm).

Trong thí nghiệm nén mẫu cao 4.0cm, việc thoát nước và bão hòa mẫu cần lưu ý điều gì?

Do chiều dài thoát nước tăng gấp đôi so với mẫu cao 2.0cm, kỹ thuật viên bắt buộc phải bố trí đá thấm có hệ số thấm $k \ge 10^{-2}\text{ cm/s}$ ở cả hai mặt trên và dưới. Đồng thời, cần quấn màng chống tràn quanh dao vòng và duy trì ngập nước liên tục nhằm đảm bảo độ bão hòa mẫu luôn đạt $S_r \ge 98%$ trong suốt quá trình gia tải.

Phương pháp Asaoka được ứng dụng như thế nào để kiểm chứng kết quả mô phỏng?

Phương pháp Asaoka sử dụng chuỗi số liệu quan trắc độ lún theo từng khoảng thời gian cố định $\Delta t = 15$ ngày để thiết lập phương trình hồi quy bậc nhất $S_i = \beta_0 + \beta_1 S_{i-1}$. Giao điểm của đường hồi quy với đường $45^\circ$ cho độ lún cố kết cuối cùng đạt 138.2cm, trùng khớp hoàn toàn với kết quả mô phỏng Plaxis 2D sử dụng mẫu nén kích thước lớn.

Kết quả nghiên cứu có thể áp dụng cho các khu vực địa chất yếu khác tại Việt Nam không?

Quy luật suy giảm ma sát thành và ảnh hưởng của kích thước dao vòng hoàn toàn có thể áp dụng cho các vùng trầm tích Holocene tại Đồng bằng Sông Cửu Long và Hải Phòng. Tuy nhiên, với tầng đất có chỉ số dẻo $I_p < 17%$ hoặc chứa lượng than bùn cao trên 15%, cần tiến hành thí nghiệm kiểm chứng riêng với các cấp áp lực từ 25 kPa đến 200 kPa.

Kết luận

  • Luận văn đã chứng minh định lượng sự phụ thuộc của các thông số biến dạng nén cố kết ($P_c, C_c, C_s, C_v, m_v$) vào đường kính và chiều cao mẫu nén trong phòng thí nghiệm.
  • Xác định mẫu nén có kích thước đường kính $7.98\text{ cm}$ và chiều cao $4.0\text{ cm}$ là kích thước tối ưu, giúp giảm sai số dự báo lún công trình xuống dưới 5% so với thực tế quan trắc.
  • Giải thích rõ cơ chế ma sát thành dao vòng làm tiêu hao từ 10% đến 18% ứng suất thẳng đứng ở các mẫu nén nhỏ ($5.08\text{ cm}$), gây đánh giá thấp độ lún công trình.
  • Xây dựng thành công mô hình tích hợp giữa thí nghiệm nén cố kết cải tiến, tính toán giải tích, mô phỏng số phần tử hữu hạn Plaxis 2D và phương pháp quan trắc Asaoka cho bài toán gia tải trước tại Quận 7, TP.HCM.
  • Đóng góp cơ sở khoa học quan trọng cho việc hoàn thiện tiêu chuẩn thí nghiệm cơ học đất và thiết kế nền móng công trình ngầm trên đất yếu đô thị.

Kế hoạch tiếp theo trong 24 tháng tới là mở rộng chương trình thực nghiệm trên 10 tổ hợp kích thước dao vòng lớn đối với các phân tầng địa chất bùn sét hữu cơ tại khu vực Đồng bằng Sông Cửu Long. Các kỹ sư và chuyên gia địa kỹ thuật được khuyến khích áp dụng ngay quy trình chọn kích thước mẫu nén chuẩn $7.98\text{ cm} \times 4.0\text{ cm}$ trong các dự án hạ tầng trọng điểm để đảm bảo tính an toàn, bền vững và tối ưu hóa chi phí đầu tư công trình.