Tổng quan nghiên cứu

Tại các khu vực đồng bằng và ven biển Việt Nam, trầm tích đất yếu chiếm tới hơn 60% diện tích tự nhiên, đặt ra thách thức lớn cho các công trình hạ tầng giao thông, khu công nghiệp và cảng biển. Những lớp bùn sét này có độ rỗng lớn, tính nén lún cao và sức chịu tải cực kỳ thấp, thường dẫn đến nguy cơ lún kéo dài hoặc trượt trồi khi thi công tải trọng lớn. Biện pháp xử lý nền bằng phương pháp gia tải đất đắp truyền thống đòi hỏi khối lượng vật liệu san lấp rất lớn và dễ gây ra hiện tượng dịch chuyển ngang hướng ra ngoài chân khối đắp, làm mất ổn định các công trình lân cận. Ngược lại, phương pháp cố kết chân không duy trì áp lực âm liên tục khoảng 80 kPa lại tạo ra xu hướng chuyển dịch đất nền hướng vào trong tâm khu vực xử lý.

Nghiên cứu tập trung giải quyết bài toán kiểm soát và triệt tiêu chuyển vị ngang của nền đất yếu thông qua việc kết hợp đồng thời tải trọng đất đắp và áp lực hút chân không cùng hệ thống bấc thấm PVD. Mục tiêu trọng tâm là làm sáng tỏ cơ chế biến dạng, thiết lập thông số không thứ nguyên để dự báo định lượng hướng chuyển vị và đánh giá các yếu tố ảnh hưởng trực tiếp đến sự ổn định của nền. Phạm vi nghiên cứu được thực hiện tại dự án Nhà máy chế biến gỗ MDF VRG Kiên Giang thuộc Khu công nghiệp Thạnh Lộc, với tổng diện tích bãi chứa gỗ là 24.836 m2, tập trung phân tích chuyên sâu trên khu vực 6 có diện tích 4.362 m2 với chiều dày tầng bùn sét yếu từ 15 m đến 20 m. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp rút ngắn thời gian cố kết từ 40% đến 50%, tiết kiệm 30% đến 45% khối lượng đất đắp và bảo đảm an toàn tuyệt đối cho kết cấu hạ tầng liền kề.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nền tảng lý thuyết của đề tài dựa trên sự tích hợp giữa lý thuyết cố kết thấm và cơ học đất biến dạng dẻo. Quá trình thoát nước cố kết một chiều theo phương đứng được mô tả qua mô hình của Terzaghi kết hợp với lời giải cố kết hướng tâm qua giếng thấm bấc thấm PVD của Barron và hiệu chỉnh của Hansbo. Bấc thấm có kích thước mặt cắt 100 mm x 4 mm được quy đổi thành đường kính tương đương khoảng 66 mm, có xét đến ảnh hưởng của vùng xáo trộn xung quanh ống dẫn với đường kính vùng xáo trộn gấp 5 đến 6 lần kích thước tương đương của ống chèn. Độ cố kết tổng thể của khối đất được xác lập thông qua phương trình Carrillo, tích hợp đồng thời quá trình tiêu tán áp lực nước lỗ rỗng theo cả phương thẳng đứng và phương ngang.

Về mặt ứng xử cơ học, lộ trình ứng suất hữu hiệu trên mặt phẳng p'-q' chỉ ra sự khác biệt bản chất: gia tải đất đắp làm tăng ứng suất tổng và đưa lộ trình ứng suất tiến sát đường phá hoại cắt, trong khi áp lực chân không duy trì ứng suất tổng không đổi, làm giảm áp lực nước lỗ rỗng và đưa lộ trình ứng suất dịch chuyển cách xa đường phá hoại, giúp bảo đảm hệ số an toàn trượt. Biến dạng phi tuyến của các lớp đất bùn sét yếu được mô phỏng bằng mô hình Modified Cam-Clay với các thông số độ dốc đường trạng thái tới hạn M, chỉ số nén lún, chỉ số nở và tỷ số rỗng ban đầu. Lớp cát đệm và các lớp sét pha dẻo cứng bên dưới được mô hình hóa bằng mô hình Mohr-Coulomb với các tham số lực dính hiệu quả, góc ma sát trong và mô đun đàn hồi. Nghiên cứu áp dụng hệ số chuyển vị không thứ nguyên DR, được định nghĩa bằng tỷ số giữa chuyển vị ngang trung bình và độ lún bề mặt cuối cùng, kết hợp với tỷ số tải trọng RL và tốc độ đắp tải LR từ 2,5 kPa/ngày đến 10,0 kPa/ngày để khảo sát định lượng ứng xử của nền.

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

Nghiên cứu sử dụng phương pháp mô phỏng số dựa trên phương pháp phần tử hữu hạn thông qua phần mềm địa kỹ thuật Geostudio Sigma/W kết hợp với dữ liệu quan trắc thực địa. Nguồn dữ liệu địa chất được trích xuất từ báo cáo khảo sát gồm 10 hố khoan từ HK1 đến HK11 với các chỉ tiêu cơ lý chi tiết của tầng bùn sét trạng thái chảy có độ ẩm tự nhiên từ 65% đến 85% và chỉ số nén Cc từ 0,6 đến 1,2. Hệ thống quan trắc hiện trường bao gồm các bàn đo lún bề mặt, đầu đo áp lực nước lỗ rỗng piezometer và các ống đo độ nghiêng inclinometer được bố trí theo chiều sâu tới 20 m nhằm ghi nhận liên tục sự biến thiên áp lực và chuyển vị trong suốt chu kỳ gia tải 120 ngày.

Lý do lựa chọn phương pháp phân tích phần tử hữu hạn 2D bài toán biến dạng phẳng là nhằm mô phỏng chính xác sự tương tác phức tạp giữa áp lực chân không phân bố đẳng hướng và tải trọng đắp phân lớp theo thời gian. Để đưa hệ thống bấc thấm 3D đối xứng trục vào mô hình 2D, tác giả áp dụng phương pháp chuyển đổi hệ số thấm tương đương của Indraratna và Chai, hiệu chỉnh hệ số thấm ngang và thấm đứng của khối đất tương đương mà vẫn bảo đảm chính xác tốc độ tiêu tán áp lực lỗ rỗng. Cỡ mẫu nghiên cứu bao gồm bộ dữ liệu thực tế tại khu vực 6 và 24 kịch bản mô phỏng số mở rộng, kết hợp giữa 6 mức tỷ số áp lực RL từ 0,5 đến 2,5 và 4 dải tốc độ đắp LR. Timeline nghiên cứu được thực hiện đồng bộ trong khoảng thời gian từ tháng 01/2015 đến tháng 06/2015.

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 và kiểm chứng dữ liệu quan trắc hiện trường đã mang lại ba phát hiện kỹ thuật cốt lõi:

Thứ nhất, tồn tại cơ chế bù trừ chuyển vị ngang rõ nét giữa áp lực chân không và tải trọng đất đắp. Tải chân không duy trì ở mức 80 kPa tạo ra lực hút kéo đất nền co cụm về phía tâm khối xử lý với chuyển vị ngang hướng vào lớn nhất tại bề mặt đất. Khi chất tải đất đắp bên trên, ứng suất cắt làm đất có xu hướng phình ra ngoài chân khối đắp. Khi kết hợp đồng thời hai tác nhân này, chuyển vị ngang trung bình của nền đất yếu giảm mạnh từ mức biến dạng trượt dương từ 15 cm đến 25 cm của phương pháp đắp đơn thuần xuống mức ổn định dao động quanh giá trị âm từ -2 cm đến -5 cm, triệt tiêu hơn 75% chuyển vị đẩy ngang gây nguy hiểm cho khu vực xung quanh.

Thứ hai, mối quan hệ giữa hệ số chuyển vị DR với tỷ số áp lực RL và tốc độ gia tải LR có tính quy luật chặt chẽ. Khi tỷ số áp lực RL nhỏ hơn 1,0, tác dụng co ngót của áp lực chân không chiếm ưu thế, dẫn đến hệ số DR nhận giá trị âm. Khi RL vượt quá 1,5, tải trọng đất đắp chiếm ưu thế khiến DR chuyển nhanh sang giá trị dương với biên độ chuyển vị ngang tăng vọt. Bên cạnh đó, khi tốc độ gia tải LR tăng gấp 4 lần, từ 2,5 kPa/ngày lên 10,0 kPa/ngày, áp lực nước lỗ rỗng thặng dư tức thời tăng thêm khoảng 35% đến 45%, thúc đẩy biến dạng dẻo trượt ngang phát triển sớm hơn và làm tăng biên độ phình hông của khối đất nền.

Thứ ba, mô hình số xây dựng bằng GeoStudio phản ánh mức độ tương thích vượt trội so với thực tế quan trắc tại dự án MDF Kiên Giang. Độ lún bề mặt tính toán theo mô hình đạt 1,82 m, so khớp gần như hoàn hảo với số liệu đo lún thực tế 1,78 m với sai số dưới 3%. Áp lực nước lỗ rỗng trong các lớp đất sét yếu tiêu tán đạt trên 85% sau 120 ngày thi công liên tục qua 4 đợt gia tải đất đắp kết hợp hút chân không, chứng minh tính tin cậy tuyệt đối của giải pháp tính toán chuyển đổi hệ số thấm tương đương.

Thảo luận kết quả

Cơ chế chuyển vị ngang trong bài toán kết hợp mang tính chất dị thời và phức tạp: biến dạng trượt do tải đất đắp mang bản chất biến dạng dẻo cắt tức thời xảy ra ngay khi chất tải, trong khi chuyển vị co ngót hướng tâm do áp lực chân không lại phụ thuộc vào tốc độ thoát nước cố kết diễn ra chậm dần theo thời gian. Do đó, trong giai đoạn đầu của từng đợt gia tải, đường cong chuyển vị thực tế luôn có xu hướng dịch chuyển nhẹ về phía ngoài trước khi co dần về phía trong tâm khi áp lực nước lỗ rỗng tiêu tán hoàn toàn.

Để trực quan hóa các quy luật phức tạp này trong công tác thiết kế, toàn bộ tập hợp dữ liệu được tổng hợp thành bảng tra cứu các giá trị DR ứng với từng cặp thông số RL và LR, đồng thời được biểu diễn qua biểu đồ quan hệ 2D giữa hệ số chuyển vị DR và tỷ số áp lực RL theo các đường cong đẳng tốc độ gia tải LR. Biểu đồ phân bố chuyển vị ngang theo chiều sâu thể hiện rõ dạng đường cong chữ S mượt mà, trong đó chuyển vị đạt cực đại âm ở độ sâu từ 0 m đến 3 m và triệt tiêu dần ở độ sâu chân bấc thấm 18 m đến 20 m. Những kết quả này hoàn toàn tương đồng với các nghiên cứu thực nghiệm tỷ lệ lớn của Ong và Chai, đồng thời giải thích rõ ràng hơn các hiện tượng nứt trượt cục bộ từng được ghi nhận trong một số dự án hạ tầng trước đây do chất tải quá nhanh.

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

Dựa trên các kết luận rút ra từ mô hình phân tích và kiểm chứng hiện trường, nghiên cứu đề xuất bốn nhóm giải pháp kỹ thuật có tính ứng dụng cao:

  1. Khống chế tỷ số tải trọng tối ưu: Đơn vị tư vấn thiết kế cần khống chế tỷ số áp lực RL trong ngưỡng từ 0,9 đến 1,2 đối với các công trình xử lý nền sát khu dân cư hoặc công trình hiện hữu. Mức tỷ số này tương ứng với áp lực đất đắp duy trì từ 70 kPa đến 95 kPa trên nền áp lực chân không 80 kPa, bảo đảm hệ số chuyển vị DR dao động an toàn trong khoảng từ -0,05 đến +0,05 để triệt tiêu biến dạng phình ngang.
  2. Kiểm soát tốc độ gia tải đất đắp: Nhà thầu thi công bắt buộc duy trì tốc độ đắp tải LR không vượt quá 5,0 kPa/ngày, chia nhỏ khối đắp thành các lớp mỏng từ 20 cm đến 30 cm theo 3 đến 4 giai đoạn gia tải. Thời gian chờ cố kết giữa mỗi đợt đắp phải đạt từ 15 đến 20 ngày nhằm bảo đảm độ tiêu tán áp lực nước lỗ rỗng thặng dư đạt tối thiểu 70% trước khi đắp lớp tiếp theo.
  3. Nâng cao chất lượng thi công màng kín khí: Đơn vị thi công chuyên nghiệp cần thi công rãnh khóa biên bentonite ngập sâu tối thiểu 1,5 m vào tầng đất sét tự nhiên có tính thấm thấp, kết hợp sử dụng lớp vải địa kỹ thuật bảo vệ không dệt để ngăn ngừa triệt để hiện tượng rách màng hoặc thủng lỗ kim, duy trì áp lực chân không ổn định từ 75 kPa đến 80 kPa liên tục suốt chu kỳ xử lý từ 90 đến 120 ngày.
  4. Áp dụng quan trắc tự động thời gian thực: Chủ đầu tư và tư vấn giám sát cần thiết lập hệ thống quan trắc chuyển vị ngang bằng inclinometer tự động với tần suất đo 2 ngày/lần trong giai đoạn gia tải và 7 ngày/lần trong giai đoạn duy trì tải. Khi tốc độ chuyển vị ngang vượt ngưỡng cảnh báo 2,0 mm/ngày, phải dừng ngay việc đắp tải để kiểm tra độ kín của hệ thống hút chân không.

Đố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 cho bốn nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế địa kỹ thuật và nền móng: Cung cấp công thức tính toán chuyển đổi hệ số thấm tương đương và biểu đồ chuẩn DR - RL để tối ưu hóa chiều cao đắp tải, giảm thiểu rủi ro tính toán chuyển vị hông trong các dự án xử lý nền đất yếu ven biển.
  • Chỉ huy trưởng và kỹ sư giám sát thi công hạ tầng: Cung cấp quy trình phân chia đợt đắp, kỹ thuật khóa rãnh màng kín khí chân không và quy chuẩn kiểm soát tốc độ gia tải nhằm bảo đảm an toàn tuyệt đối cho các công trình lân cận trong phạm vi ảnh hưởng 30 m đến 50 m.
  • Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Địa kỹ thuật - Công trình: Làm tài liệu tham khảo chuyên sâu về mô hình hóa bài toán cố kết thấm kết hợp trong GeoStudio Sigma/W, kết nối lý thuyết cơ học đất dẻo Modified Cam-Clay với bài toán thực tế.
  • Ban quản lý dự án và chủ đầu tư hạ tầng khu công nghiệp, cảng biển: Cung cấp cơ sở khoa học để thẩm tra biện pháp thi công, tiết kiệm từ 20% đến 30% tổng chi phí xử lý nền và rút ngắn tiến độ tổng thể của các dự án quy mô từ 20.000 m2 đến 50.000 m2.

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

Phương pháp hút chân không kết hợp đất đắp mang lại ưu điểm vượt trội gì so với gia tải thông thường?

Phương pháp kết hợp giúp gia tăng áp lực hữu hiệu lên tới hơn 160 kPa gần như tức thời mà không gây nguy cơ mất ổn định trượt sâu. Áp lực chân không 80 kPa tạo hiệu ứng co ngót thể tích, giúp triệt tiêu từ 70% đến 80% chuyển vị đẩy ngang do tải trọng đất đắp gây ra, đồng thời giảm khoảng 40% khối lượng cát đắp và rút ngắn thời gian xử lý xuống còn 3 đến 4 tháng.

Hệ số chuyển vị không thứ nguyên DR có ý nghĩa như thế nào trong thiết kế?

Hệ số DR là tỷ số giữa chuyển vị ngang trung bình và độ lún bề mặt cuối cùng. Thông số này phản ánh xu hướng biến dạng tổng thể của khối đất nền: khi DR mang dấu âm thì nền đất co cụm vào trong tâm, khi DR mang dấu dương thì đất phình ra ngoài. Việc kiểm soát DR ở mức xấp xỉ 0 giúp kỹ sư loại bỏ hoàn toàn nguy cơ nứt gãy công trình liền kề.

Làm thế nào để xác định tốc độ gia tải đất đắp an toàn trong quá trình thi công?

Tốc độ gia tải đất đắp cần được xác định dựa trên tốc độ tiêu tán áp lực nước lỗ rỗng đo đạc bằng piezometer. Trong điều kiện bùn sét yếu tại Kiên Giang, tốc độ đắp an toàn được khuyến nghị từ 2,5 kPa/ngày đến 5,0 kPa/ngày. Tốc độ này bảo đảm ứng suất cắt không vượt quá sức chống cắt không thoát nước tăng thêm của khối đất trong quá trình cố kết.

Tại sao cần quy đổi bài toán bấc thấm từ mô hình 3D đối xứng trục sang mô hình 2D biến dạng phẳng?

Trong một dự án xử lý nền, hàng nghìn tim bấc thấm được cắm theo dạng lưới tam giác hoặc ô vuông 3D. Việc mô phỏng chi tiết từng tim bấc trong bài toán 3D đòi hỏi tài nguyên tính toán khổng lồ. Phương pháp quy đổi hệ số thấm tương đương của Indraratna đưa hệ thống về dạng dải thấm phẳng 2D nhưng vẫn bảo đảm sai số độ lún và chuyển vị dưới 5%.

Khi nền đất có tầng cát thấm nước nằm xen kẹt thì hiệu quả hút chân không bị ảnh hưởng ra sao?

Các thấu kính cát hoặc tầng cát xen kẹt có hệ số thấm lớn sẽ dẫn không khí và nước ngầm từ bên ngoài vào vùng xử lý, gây thất thoát áp lực chân không nghiêm trọng từ mức thiết kế 80 kPa xuống dưới 40 kPa. Khi đó, kỹ sư bắt buộc phải thi công tường hào sét bentonite hoặc cọc xi măng đất bao quanh chu vi để cô lập hoàn toàn tầng thấm.

Kết luận

  • Luận văn đã phân tích toàn diện cơ chế bù trừ chuyển vị ngang giữa áp lực chân không 80 kPa và tải trọng đất đắp trên nền bùn sét yếu, cung cấp giải pháp hữu hiệu để kiểm soát biến dạng công trình.
  • Thiết lập thành công mối quan hệ định lượng giữa hệ số chuyển vị không thứ nguyên DR với tỷ số tải trọng RL và tốc độ đắp LR thông qua 24 kịch bản mô phỏng số chi tiết.
  • Kiểm chứng mô hình phần tử hữu hạn bằng phần mềm GeoStudio trên số liệu thực tế dự án Nhà máy MDF Kiên Giang, đạt độ chính xác cao với sai số độ lún bề mặt dưới 3% so với thực tế 1,78 m.
  • Xác lập ngưỡng tốc độ đắp an toàn dưới 5,0 kPa/ngày và tỷ số áp lực tối ưu từ 0,9 đến 1,2 nhằm bảo đảm chuyển vị ngang trung bình của nền đất duy trì trong giới hạn an toàn.
  • Đóng góp bộ hướng dẫn thực hành thi công và bố trí quan trắc hiện trường, tạo tiền đề nâng cao chất lượng thi công các công trình hạ tầng trên nền đất yếu.

Trong giai đoạn 2026 - 2028, các nghiên cứu tiếp theo cần mở rộng đánh giá mô hình 3D phi tuyến xét đến ảnh hưởng của biến dạng từ biến thứ cấp và kiểm chuẩn trên các dự án đường cao tốc quy mô lớn tại Đồng bằng sông Cửu Long. Các kỹ sư và chuyên gia địa kỹ thuật quan tâm có thể ứng dụng ngay quy trình tính toán chuyển đổi hệ số thấm tương đương và biểu đồ chuẩn DR - RL vào các dự án xử lý nền thực tế để tối ưu hóa chi phí và bảo đảm an toàn công trình.