Tổng quan nghiên cứu

Khu vực Thành phố Hồ Chí Minh có sự phân bố rộng khắp của các tầng trầm tích đất yếu tuổi Holocene với chiều dày lớp bùn sét lên tới 30m. Trong thực tế xây dựng hạ tầng tại khu vực Quận 2 (nay thuộc Thành phố Thủ Đức), hiện tượng biến dạng lớn, lún nứt kéo dài đã xảy ra ở nhiều dự án trọng điểm như Đại lộ Đông Tây hay dự án nâng cấp mở rộng đường Nguyễn Hữu Cảnh, gây phát sinh hàng trăm tỷ đồng chi phí xử lý và kéo dài tiến độ thi công. Thách thức lớn đối với các kỹ sư địa kỹ thuật là việc xác định chính xác sức chống cắt không thoát nước ($S_u$) – thông số quyết định trực tiếp đến sức chịu tải của nền móng và độ ổn định hố đào.

Các thí nghiệm chuyên sâu như nén ba trục không cố kết không thoát nước ($S_{u,UU}$) và cắt cánh hiện trường ($S_{u,VS}$) thường tốn kém chi phí và đòi hỏi thiết bị phức tạp, trong khi thí nghiệm cắt trực tiếp ($S_{u,DS}$) lại phổ biến và có chi phí thấp hơn từ 40% đến 60%. Luận văn tập trung nghiên cứu đặc điểm địa kỹ thuật đặc trưng, ứng dụng phần mềm GMS phiên bản 7.1 để mô hình hóa không gian chiều dày tầng đất bùn sét, đồng thời thiết lập các phương trình tương quan định lượng giữa $S_{u,DS}$ với $S_{u,UU}$ và $S_{u,VS}$ trong phạm vi độ sâu từ 0m đến 16m tại Quận 2. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn, giúp tối ưu hóa từ 20% đến 30% chi phí khảo sát địa chất, rút ngắn 50% thời gian đánh giá sơ bộ địa tầng và nâng cao độ an toàn cho các công trình xây dự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 được xây dựng trên nền tảng lý thuyết địa thống kê (Geostatistics) với phương pháp nội suy Kriging thông thường (Ordinary Kriging) và lý thuyết hồi quy thống kê (Regression Analysis). Phương pháp Ordinary Kriging dựa trên giả thuyết hàm ngẫu nhiên dừng bậc hai và cấu trúc hàm Variogram thực nghiệm $\gamma(h) = \frac{1}{2N(h)} \sum [Z(x) - Z(x+h)]^2$ cùng hàm hiệp phương sai Covariance để ước lượng tối ưu chiều dày lớp đất tại các vị trí chưa khoan mà không tạo ra sai số hệ thống, với mô hình Spherical có các tham số thực nghiệm: Nugget = 0,34371546, Contribution = 43,535 và Range = 0,05169.

Lý thuyết hồi quy tuyến tính đơn và phi tuyến tính áp dụng phương pháp bình phương cực tiểu (Least Squares Method) để xác lập mối liên hệ giữa các biến $S_u$. Các mô hình được kiểm định nghiêm ngặt với khoảng tin cậy $CI = 95%$, mức ý nghĩa thống kê $p < 0,05$ và hệ số xác định $R^2 \ge 50%$. Ba khái niệm cốt lõi bao gồm: sức chống cắt không thoát nước từ thí nghiệm cắt trực tiếp phẳng ($S_{u,DS}$), thí nghiệm nén ba trục ($S_{u,UU}$), thí nghiệm cắt cánh tại hiện trường ($S_{u,VS}$), kết hợp cùng hai chỉ tiêu phân loại đất quan trọng là hệ số rỗng tự nhiên ($e$) và chỉ số độ sệt ($LI$).

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

Nguồn dữ liệu của luận văn được tổng hợp từ 28 hố khoan khảo sát địa chất công trình chuẩn tọa độ WGS-84 cùng hệ thống 578 mẫu thí nghiệm cắt trực tiếp ở độ sâu dưới 28m, 305 điểm đo cắt cánh hiện trường ở độ sâu dưới 16m và 26 mẫu nén ba trục UU ở độ sâu dưới 14m từ các đơn vị tư vấn uy tín như Trung tâm Nghiên cứu Công nghệ và Thiết bị Công nghiệp - Trường Đại học Bách Khoa TP.HCM, Công ty Tư vấn SCE và Công ty Kỹ thuật Xây dựng Tung Feng.

Phương pháp chọn mẫu phân tầng có chủ đích (Stratified Purposive Sampling) được thực hiện bằng cách chia dữ liệu theo từng khoảng độ sâu 2m (từ 0-2m, 2-4m cho đến 14-16m) và theo 9 dải độ sệt $LI$ từ 1,00 đến 1,90. Để đảm bảo tính đại diện và công suất kiểm định thống kê đạt $Power = 95%$ ($\beta = 0,05$) với mức sai lầm loại I $\alpha = 0,05$ và hệ số tương quan kỳ vọng $r = 0,90$, cỡ mẫu tối thiểu được tính toán theo công thức chuyên ngành cho ra $N_{min} = 9$. Trên thực tế, cỡ mẫu tại mỗi dải phân tích $N_{DS}$ đều đạt từ 14 đến 100 mẫu, vượt xa ngưỡng yêu cầu. Việc lựa chọn công cụ phần mềm Minitab phiên bản 16 và GMS phiên bản 7.1 cho phép xử lý đồng bộ chuỗi dữ liệu trong giai đoạn nghiên cứu năm 2012, đảm bảo tính chuẩn xác và loại trừ hoàn toàn các sai số ngẫu nhiên.

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

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

Phân tích thống kê 28 hố khoan và hàng trăm mẫu thí nghiệm đã xác lập bộ chỉ tiêu cơ lý đặc trưng của bùn sét Quận 2: độ ẩm tự nhiên trung bình $\bar{W} = 79,90%$, khối lượng thể tích tự nhiên $\bar{\gamma} = 1,50\text{ g/cm}^3$, tỷ trọng $\bar{\Delta} = 2,63\text{ g/cm}^3$, hệ số rỗng $\bar{e} = 2,17 > 1,5$, chỉ số dẻo $\bar{I_p} = 31,76%$, độ sệt $\bar{LI} = 1,33 > 1$, góc ma sát trong $\bar{\varphi} = 4,01^\circ$ và lực dính $\bar{c} = 0,06\text{ kG/cm}^2$ (~6 kPa). Các chỉ số này khẳng định lớp bùn sét thuộc nhóm đất yếu có độ rỗng cao, trạng thái chảy dẻo và sức chịu tải rất kém.

Mô hình nội suy Ordinary Kriging cho thấy chiều dày lớp bùn sét tăng dần rõ rệt từ phía Đông Nam sang Tây Bắc. Tại khu vực Cát Lái và Thạnh Mỹ Lợi, chiều dày bùn sét xấp xỉ 0m, trong khi tại khu vực Thảo Điền và An Phú, chiều dày phát triển mạnh từ 14m lên hơn 28m với chiều dày trung bình toàn quận đạt 11,62m. Kiểm chứng thực tế tại Trường Quốc tế Anh (British International School, số 215 Nguyễn Văn Hưởng, Thảo Điền) ghi nhận chiều dày thực tế 19,0m so với giá trị nội suy 18,0m (độ sai lệch chỉ -5,26%), và tại Chùa Huệ Nghiêm 2 (An Phú) là 14,0m so với 16,0m (sai số 14,28%).

Thiết lập thành công phương trình tương quan tuyến tính giữa thí nghiệm cắt trực tiếp và nén ba trục theo độ sâu $z \le 14\text{m}$: $$S_{u,UU} = -0,148 + 1,293 S_{u,DS}$$ với hệ số xác định $R^2 = 98,2%$ và $p < 0,05$.

Thiết lập thành công phương trình tương quan phi tuyến đa thức bậc hai giữa thí nghiệm cắt trực tiếp và cắt cánh hiện trường theo độ sâu $z \le 16\text{m}$: $$S_{u,VS} = -29,96 + 8,505 S_{u,DS} - 0,2408 S_{u,DS}^2$$ với $R^2 = 94,5%$ và $p < 0,05$ (trong khi mô hình tuyến tính đơn chỉ đạt $R^2 = 15,8%$ nên bị loại bỏ).

Theo chỉ số độ sệt trong khoảng $LI \in [1,0; 1,8]$, phương trình tuyến tính đạt chuẩn là $S_{u,UU} = -39,92 + 4,044 S_{u,DS}$ với $R^2 = 76,3%$. Mô hình tương quan giữa $S_{u,DS}$ và $S_{u,VS}$ theo độ sệt bị loại bỏ do $R^2$ thấp (46%) và hàm logarit có $p > 0,05$.

Thảo luận kết quả

Kết quả phân tích chỉ ra rằng sức chống cắt từ thí nghiệm cắt cánh hiện trường ($S_{u,VS}$) luôn cao hơn đáng kể so với thí nghiệm cắt trực tiếp trong phòng ($S_{u,DS}$), với tỷ lệ chênh lệch từ 1,14 lần đến 3,50 lần (đạt cực đại tại độ sâu 4-6m với $S_{u,VS} = 36,32\text{ kPa}$ so với $S_{u,DS} = 10,37\text{ kPa}$). Nguyên nhân cốt lõi do quá trình lấy mẫu nguyên dạng từ độ sâu lòng đất lên mặt đất gây xáo động cấu trúc trầm tích sét, giải phóng ứng suất hữu hiệu và phá hủy một phần liên kết xi măng tự nhiên giữa các hạt sét. Thí nghiệm cắt cánh hiện trường đo đạc trực tiếp trong lòng đất nên phản ánh chính xác trạng thái ứng suất nguyên trạng và tính dị hướng của khối đất.

Đối với thí nghiệm ba trục, giá trị $S_{u,UU}$ cao hơn $S_{u,DS}$ từ 1,20 đến 1,37 lần nhờ buồng áp lực ba trục tái lập được một phần áp lực buồng vây, hạn chế sự phá hoại cục bộ dọc theo mặt cắt cưỡng bức phẳng như trong hộp cắt trực tiếp.

Dữ liệu nghiên cứu được trình bày trực quan thông qua đồ thị đường khớp hồi quy (Fitted Line Plot) với dải tin cậy 95% Confidence Interval (CI) và dải dự báo 95% Prediction Interval (PI) trên Minitab 16, kết hợp bảng ma trận thống kê đa biến và bản đồ bình đồ đẳng trị 3D trên nền tảng GMS 7.1. Sự tương đồng cao với các công bố địa kỹ thuật tại TP.HCM (như các nghiên cứu của trường Đại học Bách Khoa ghi nhận $R^2$ dao động 97%-99%) khẳng định tính quy luật và độ tin cậy vượt trội của các mô hình toán học đã thiết lập.

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

Áp dụng phương trình phi tuyến $S_{u,VS} = -29,96 + 8,505 S_{u,DS} - 0,2408 S_{u,DS}^2$ và phương trình tuyến tính $S_{u,UU} = -0,148 + 1,293 S_{u,DS}$ vào công tác tính toán thiết kế sơ bộ nền móng công trình và hố đào sâu tại khu vực Quận 2. Chủ thể thực hiện là các kỹ sư tư vấn địa kỹ thuật, với mục tiêu cắt giảm từ 20% đến 30% khối lượng thí nghiệm cắt cánh và nén ba trục trong giai đoạn thiết kế cơ sở, áp dụng ngay trong năm 2026.

Tích hợp bản đồ nội suy Kriging 3D về chiều dày tầng bùn sét vào quy hoạch không gian ngầm và hệ thống hạ tầng kỹ thuật đô thị Thành phố Thủ Đức. Sở Xây dựng TP.HCM và các Ban Quản lý Dự án cần khoanh vùng các khu vực có chiều dày bùn sét lớn hơn 20m tại Thảo Điền và An Phú nhằm đưa ra yêu cầu bắt buộc về giải pháp xử lý nền đất yếu (bấc thấm hút chân không, cọc đất xi măng), hoàn thành triển khai trong vòng 6 đến 12 tháng.

Chuẩn hóa hệ số suy giảm độ bền do xáo động mẫu trong quy trình thí nghiệm cắt trực tiếp tại các phòng kiểm định LAS-XD. Các đơn vị thí nghiệm cần đưa tỷ số chuyển đổi thực nghiệm từ 1,2 đến 1,3 lần vào báo cáo địa chất công trình để phản ánh đúng sức chịu tải thực tế, hoàn thành việc ban hành cẩm nang kỹ thuật nội bộ trong thời hạn 3 tháng tới.

Mở rộng cơ sở dữ liệu địa kỹ thuật công trình cho các tầng sâu trên 16m và các khu vực lân cận như Quận 7, Quận 9, huyện Nhà Bè. Các viện nghiên cứu và trường đại học chuyên ngành cần phối hợp triển khai thu thập bổ sung 500 hố khoan sâu và 1.000 mẫu thí nghiệm nén ba trục cố kết triệt để (CU, CD) trong kế hoạch nghiên cứu 24 tháng tới.

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

Kỹ sư thiết kế địa kỹ thuật và công trình ngầm: Sử dụng các phương trình tương quan để ngoại suy nhanh thông số $S_{u,UU}$ và $S_{u,VS}$ từ dữ liệu cắt trực tiếp thông thường, phục vụ tính toán sức chịu tải cọc ma sát, kiểm tra ổn định trượt cung tròn bờ kênh sông Sài Gòn và thiết kế hệ tường vây tầng hầm.

Chủ đầu tư và nhà thầu xây dựng: Tham khảo để lập kế hoạch khảo sát địa chất tối ưu, giảm thiểu 15% đến 25% ngân sách dành cho các thí nghiệm hiện trường đắt đỏ mà vẫn đảm bảo độ tin cậy khoa học phục vụ thẩm định dự án.

Các đơn vị tư vấn khảo sát địa chất công trình: Tiếp cận quy trình chuẩn về ứng dụng phương pháp địa thống kê Ordinary Kriging và phần mềm GMS 7.1 nhằm tự động hóa việc lập bình đồ phân bố địa tầng 3D phục vụ các dự án khu đô thị mới.

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 phương pháp luận toán thống kê, kỹ thuật kiểm định mô hình hồi quy phi tuyến và bộ cơ sở dữ liệu mẫu phong phú của vùng đất yếu Holocene TP.HCM cho các đề tài nghiên cứu chuyên sâu.

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

Sức chống cắt không thoát nước đóng vai trò gì trong xây dựng trên nền bùn sét?

Sức chống cắt không thoát nước ($S_u$) là chỉ tiêu then chốt phản ánh khả năng chịu tải của đất sét yếu bão hòa nước trong giai đoạn thi công nhanh khi nước lỗ rỗng chưa kịp thoát. Thông số này quyết định trực tiếp đến sức chịu tải của móng, độ ổn định thành hố đào và tính toán áp lực đất tác dụng lên kết cấu chắn giữ.

Vì sao sức chống cắt từ thí nghiệm cắt cánh hiện trường luôn lớn hơn thí nghiệm cắt trong phòng?

Thí nghiệm cắt cánh ($S_{u,VS}$) đo trực tiếp trong khối đất nguyên trạng, giữ nguyên ứng suất tự nhiên và cấu trúc gắn kết của các hạt khoáng vật sét. Ngược lại, mẫu đất thí nghiệm cắt trực tiếp trong phòng ($S_{u,DS}$) bị xáo động do quá trình khoan lấy mẫu, vận chuyển và giải phóng ứng suất, khiến độ bền chống cắt bị suy giảm từ 1,14 đến 3,50 lần.

Tại sao mô hình tương quan giữa thí nghiệm cắt trực tiếp và cắt cánh theo độ sâu lại là hàm phi tuyến bậc hai?

Khi phân tích tương quan giữa $S_{u,DS}$ và $S_{u,VS}$, mô hình tuyến tính đơn chỉ đạt hệ số xác định $R^2 = 15,8%$ (bị loại bỏ). Việc sử dụng hàm đa thức bậc hai $S_{u,VS} = -29,96 + 8,505 S_{u,DS} - 0,2408 S_{u,DS}^2$ phản ánh chính xác sự thay đổi phi tuyến của ứng suất dính bám theo độ sâu, nâng hệ số $R^2$ lên $94,5%$ với $p < 0,05$.

Mô hình nội suy Ordinary Kriging có độ chính xác như thế nào trong dự báo chiều dày bùn sét?

Phương pháp Ordinary Kriging trên phần mềm GMS 7.1 có độ tin cậy rất cao nhờ tối thiểu hóa phương sai sai số. Kiểm chứng đối chiếu tại 28 hố khoan thực tế ở Quận 2 cho thấy độ sai lệch chiều dày chỉ dao động từ -5,26% tại phường Thảo Điền đến 14,28% tại phường An Phú, hoàn toàn nằm trong giới hạn cho phép của tiêu chuẩn khảo sát xây dựng.

Vì sao tương quan giữa các loại thí nghiệm theo độ sệt không được khuyến nghị sử dụng rộng rãi?

Mặc dù tương quan $S_{u,UU} - S_{u,DS}$ theo độ sệt $LI \in [1,0; 1,8]$ đạt $R^2 = 76,3%$, nhưng mô hình tương quan $S_{u,VS} - S_{u,DS}$ theo $LI$ chỉ đạt $R^2 = 46,0%$ và mô hình logarit có $p > 0,05$. Độ sệt chịu ảnh hưởng lớn bởi độ ẩm cục bộ và giới hạn Atterberg trong phòng thí nghiệm, do đó không phản ánh ổn định bằng biến số độ sâu địa tầng.

Kết luận

  • Xác lập hệ thống chỉ tiêu cơ lý đặc trưng của bùn sét Quận 2 với độ ẩm tự nhiên $79,90%$, hệ số rỗng $2,17$, độ sệt $1,33$, khẳng định trạng thái đất rất yếu và độ nén lún cao.
  • Xây dựng thành công mô hình 3D Ordinary Kriging bằng phần mềm GMS 7.1, làm rõ quy luật biến thiên chiều dày bùn sét từ 0m (Thạnh Mỹ Lợi) lên hơn 28m (Thảo Điền) với sai số kiểm chứng chỉ từ 5,26%.
  • Thiết lập phương trình hồi quy tuyến tính $S_{u,UU} = -0,148 + 1,293 S_{u,DS}$ ($R^2 = 98,2%$) cho độ sâu $z \le 14\text{m}$, cho phép ngoại suy tin cậy kết quả nén ba trục từ thí nghiệm cắt trực tiếp.
  • Thiết lập phương trình hồi quy phi tuyến bậc hai $S_{u,VS} = -29,96 + 8,505 S_{u,DS} - 0,2408 S_{u,DS}^2$ ($R^2 = 94,5%$) cho độ sâu $z \le 16\text{m}$, giải quyết triệt để bài toán ước tính sức chống cắt hiện trường.
  • Đóng góp giải pháp toán thống kê thực tiễn giúp ngành xây dựng tối ưu hóa chi phí khảo sát địa chất và nâng cao độ an toàn nền móng công trình trên vùng đất yếu TP.HCM.

Trong giai đoạn 2026-2030, hướng nghiên cứu tiếp theo cần mở rộng thu thập dữ liệu tại các độ sâu trên 16m và liên kết mạng lưới địa tầng trên toàn địa bàn Thành phố Thủ Đức. Các tổ chức tư vấn thiết kế và chủ đầu tư hãy ứng dụng ngay bộ công thức tương quan thực nghiệm này vào quy trình thẩm tra, thiết kế móng để tối ưu hóa hiệu quả kinh tế và đảm bảo tính bền vững lâu dài cho công trình.