Tổng quan nghiên cứu

Quá trình đô thị hóa nhanh chóng và nhu cầu phát triển cơ sở hạ tầng tại các dải đồng bằng ven biển trên thế giới đang đối mặt với thách thức địa kỹ thuật nghiêm trọng do sự hiện diện của các tầng đất sét yếu có chiều dày lớn từ 30 mét đến hơn 70 mét. Điển hình như tại khu vực cửa sông Nakdong và vùng phát triển Cảng mới Busan, trầm tích sét mềm có sức chịu tải cực thấp, hệ số rỗng cao và khả năng thấm nước rất kém. Giải pháp kỹ thuật phổ biến và tiết kiệm chi phí nhất hiện nay là sử dụng bấc thấm thẳng đứng kết hợp tải trọng đắp gia tải trước nhằm rút ngắn chiều dài thoát nước theo phương ngang từ hàng chục mét xuống còn khoảng 0,8 mét đến 1,5 mét, qua đó thúc đẩy nhanh quá trình cố kết và tăng cường cường độ chống cắt của đất nền.

Tuy nhiên, các phương pháp tính toán và thiết kế truyền thống thường giả định khả năng thoát nước của bấc thấm là hằng số không đổi và bỏ qua sự xáo động của đất xung quanh bấc do quá trình đóng quả ép tạo ra. Thực tế thi công cho thấy sự suy giảm khả năng thoát nước của bấc và sự hình thành vùng xáo động làm chậm đáng kể tốc độ tiêu tán áp lực nước lỗ rỗng, dẫn đến sai số dự báo độ lún và thời gian gia tải có thể lên tới 25% đến 40%. Mục tiêu trọng tâm của nghiên cứu này là xây dựng hệ thống lời giải giải tích và mô hình số trị tiên tiến dựa trên lý thuyết biến dạng lớn, tích hợp đồng thời sự suy giảm khả năng thoát nước phi tuyến theo độ sâu và thời gian cùng với sự biến thiên phi tuyến của hệ số thấm và tính nén lún trong vùng xáo động. Kết quả nghiên cứu được kiểm chứng qua các dữ liệu quan trắc thực tế tại Sân bay Saga, Cảng mới Busan, tuyến cao tốc Pacific Highway và các bãi đắp thử nghiệm tại Malaysia, Trung Quốc. Nghiên cứu mang lại ý nghĩa to lớn trong việc nâng cao độ tin cậy của dự báo độ lún lên trên 90%, giảm thiểu nguy cơ lún lệch kéo dài và tối ưu hóa thời gian chờ gia tải từ 3 tháng đến 6 tháng cho các siêu dự án hạ tầng.

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

Khung lý thuyết áp dụng

Khung lý thuyết của nghiên cứu được phát triển dựa trên nền tảng lý thuyết cố kết thấm hướng tâm kinh điển của Barron và Hansbo, kết hợp với lý thuyết biến dạng lớn cho dòng thấm hướng tâm. Mô hình giải tích xem xét một ô đơn vị đối xứng trục bao gồm bấc thấm ở trung tâm, bao quanh bởi vùng xáo động và vùng đất nguyên dạng, phía trên có lớp đệm cát thoát nước dày khoảng 0,5 mét. Để mô tả chính xác ứng xử cơ học của đất nền chịu tải trọng đắp lớn, mô hình đất yếu Soft Soil và tiêu chuẩn phá hoại Mohr-Coulomb được tích hợp vào các phân tích số nhằm phản ánh mối quan hệ phi tuyến giữa chỉ số nén lún, chỉ số trương nở và trạng thái ứng suất hiệu quả.

Ba khái niệm học thuật cốt lõi chi phối toàn bộ mô hình nghiên cứu bao gồm:

  • Vùng xáo động: Khu vực đất xung quanh bấc bị phá vỡ cấu trúc trong quá trình ép quả ép, có bán kính biến thiên từ 1,5 đến 5,0 lần bán kính tương đương của quả ép. Trong vùng này, hệ số thấm ngang suy giảm mạnh xuống chỉ còn khoảng 0,2 đến 0,5 lần hệ số thấm tự nhiên, trong khi hệ số nén lún thể tích tăng từ 1,2 đến 2,0 lần.
  • Sức cản giếng và sự suy giảm khả năng thoát nước: Lưu lượng dòng chảy dọc theo thân bấc suy giảm phi tuyến theo thời gian do hiện tượng uốn nếp khi đất lún lớn và giảm theo chiều sâu dưới tác dụng của ứng suất nén hông tăng dần từ 50 kPa đến hơn 400 kPa.
  • Phương pháp biến dạng phẳng tương đương: Chuyển đổi bài toán 3D hoặc đối xứng trục của mạng lưới bấc thấm về bài toán 2D biến dạng phẳng bằng cách hiệu chỉnh hệ số thấm đứng tương đương và bề rộng tương đương của phần tử đặc.

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

Nghiên cứu sử dụng nguồn dữ liệu đa dạng bao gồm chuỗi số liệu quan trắc hiện trường dài hạn từ 5 dự án đắp thử nghiệm quy mô lớn tại các quốc gia gồm Nhật Bản, Hàn Quốc, Úc, Trung Quốc và Malaysia. Cụ thể tại Cảng mới Busan, dữ liệu được thu thập từ 12 mặt cắt địa kỹ thuật tại Phân kỳ 2-5 với tầng đất yếu dày từ 30 mét đến 60 mét. Cỡ mẫu thực nghiệm trong phòng bao gồm hơn 40 mẫu đất sét nguyên dạng được lấy bằng ống mẫu thành mỏng piston tại các độ sâu từ 2,0 mét đến 35,0 mét, phục vụ cho các thí nghiệm nén cố kết Oedometer tiêu chuẩn và thí nghiệm cố kết trên mẫu khối lớn có đường kính 300 mm. Phương pháp chọn mẫu phi xác suất có chủ đích được áp dụng nhằm tập trung vào các tầng sét bùn có độ dẻo cao với giới hạn chảy từ 45% đến 85% và hệ số rỗng ban đầu dao động từ 1,4 đến 2,8.

Phương pháp phân tích chủ đạo là sự kết hợp giữa giải tích toán học bằng phương pháp sai phân hữu hạn và mô phỏng số thông qua phương pháp phần tử hữu hạn. Lý do lựa chọn phương pháp phân tích này là vì các giải pháp giải tích đóng truyền thống không thể giải quyết được hệ phương trình vi phân phi tuyến bậc hai khi các thông số độ thấm, tính nén lún và khả năng thoát nước thay đổi liên tục theo không gian và thời gian trong điều kiện biến dạng lớn. Toàn bộ quy trình thu thập dữ liệu, lập trình thuật toán giải tích và kiểm chứng mô hình mô phỏng số được thực hiện nghiêm ngặt theo tiến trình nghiên cứu kéo dài 36 thá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 đã mang lại 4 phát hiện đột phá về mặt lý thuyết và ứng dụng thực tiễn trong công nghệ gia cố nền đất yếu:

Thứ nhất, nghiên cứu đã thiết lập thành công mô hình phân bố phi tuyến dạng parabol cho hệ số thấm ngang và hệ số nén lún thể tích trong vùng xáo động. Kết quả phân tích chỉ ra rằng khi xét đến sự gia tăng hệ số nén lún lên 1,5 lần tại sát bề mặt bấc thấm, tốc độ cố kết thực tế của nền đất bị chậm lại từ 18% đến 32% so với mô hình Hansbo kinh điển vốn chỉ giả định hệ số nén lún là một hằng số cố định.

Thứ hai, việc áp dụng lý thuyết biến dạng lớn đã làm sáng tỏ cơ chế suy giảm khả năng thoát nước của bấc thấm theo thời gian dưới tác động của biến dạng đứng lớn. Tại bãi đắp thử nghiệm Sân bay Saga, khả năng thoát nước của bấc bị suy giảm hơn 65% sau 180 ngày thi công khi độ biến dạng đứng vượt quá ngưỡng 20%. Mô hình đề xuất đã đạt hệ số tương quan rất cao với giá trị R bình phương trên 0,94 khi so sánh với độ lún đo được ngoài hiện trường lên đến 1,85 mét.

Thứ ba, phương pháp hệ số thấm đứng tương đương cải tiến có xét đến sự suy giảm khả năng thoát nước phi tuyến theo chiều sâu đã giải quyết xuất sắc bài toán cố kết cho tầng đất dày hơn 45 mét tại Cảng mới Busan. Áp lực nước lỗ rỗng thặng dư tại độ sâu 20 mét tiêu tán từ 120 kPa ban đầu xuống dưới 15 kPa sau 400 ngày gia tải, đạt độ trùng khớp ấn tượng với sai số dự báo độ lún dưới 5% so với số liệu đo đạc thực tế.

Thứ tư, mô hình biến dạng phẳng 2D sử dụng phần tử đặc đã chứng minh tính hiệu quả vượt trội trong phân tích đa bấc thấm. Bằng cách cân bằng thể tích nước thoát ra và tỷ lệ diện tích tương đương, mô hình 2D đã giảm hơn 85% thời gian tính toán máy tính so với mô hình 3D hoàn chỉnh, trong khi vẫn duy trì độ chính xác cao với sai số độ lún bề mặt chỉ dao động từ 3% đến 6% tại các dự án đắp thử nghiệm ở miền Đông Trung Quốc và Malaysia.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự chậm trễ cố kết trong các mô hình thực tế xuất phát từ hai yếu tố vật lý tương hỗ. Quá trình đóng quả ép bằng búa rung hoặc kích tĩnh đã phá vỡ hoàn toàn cấu trúc hạt sét tự nhiên, làm giảm kích thước các vi lỗ rỗng và làm tăng tính nén lún thể tích của lớp đất sát bấc. Đồng thời, ứng suất nén hông gia tăng theo độ sâu từ 20 kPa ở mặt đất lên hơn 300 kPa tại các tầng sâu đã bóp nghẹt diện tích mặt cắt rãnh thoát nước của lõi nhựa, gây ra hiện tượng nghẽn dòng cục bộ.

So sánh với các nghiên cứu trước đây của Hansbo, Walker và Indraratna hay Basu, các mô hình trước thường đánh giá quá cao tốc độ thoát nước trong 60 đến 90 ngày đầu do giả định bấc thấm có khả năng xả lý tưởng. Cách tiếp cận của luận văn đã lấp đầy khoảng trống này bằng các hàm phân bố phi tuyến liên tục. Về mặt trực quan hóa, toàn bộ dữ liệu quan trắc và tính toán có thể được biểu diễn một cách rõ nét thông qua biểu đồ đường cong quan hệ giữa độ lún theo thời gian trên hệ tọa độ bán logarit, biểu đồ phân bố áp lực nước lỗ rỗng thặng dư theo chiều sâu từ 0 đến 50 mét, cùng bảng tổng hợp so sánh sai số tương đối giữa phương pháp cổ điển với sai số 22% đến 38% và phương pháp đề xuất với sai số chỉ từ 3% đến 6%.

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

Nhằm chuyển hóa các kết quả nghiên cứu thành quy chuẩn thực hành trong kỹ thuật xây dựng công trình trên nền đất yếu, 4 nhóm giải pháp cụ thể được đề xuất như sau:

Tích hợp công thức phân bố phi tuyến của hệ số thấm và tính nén lún vùng xáo động vào quy trình tính toán thiết kế chi tiết để xác định khoảng cách cắm bấc thấm tối ưu từ 1,0 mét đến 1,3 mét theo mạng tam giác, hướng tới mục tiêu đạt độ cố kết 90% trong thời gian gia tải dưới 180 ngày. Giải pháp này do Kỹ sư trưởng thiết kế địa kỹ thuật chủ trì hoàn thành trong giai đoạn lập hồ sơ thiết kế bản vẽ thi công trong thời hạn 60 ngày trước ngày khởi công.

Kiểm tra và đánh giá nghiêm ngặt khả năng thoát nước thực tế của bấc thấm dưới các cấp áp lực nén hông từ 100 kPa đến 400 kPa trong phòng thí nghiệm, đảm bảo khả năng thoát nước duy trì của vật liệu bấc không thấp hơn 15 mét khối mỗi năm ở điều kiện làm việc khắc nghiệt nhất. Đơn vị tư vấn thí nghiệm độc lập chịu trách nhiệm thực hiện và hoàn thành chứng nhận chất lượng trong vòng 30 ngày trước khi nhập vật tư đại trà.

Áp dụng phương pháp phần tử đặc trong mô hình biến dạng phẳng 2D để mô phỏng sự tương tác giữa nền đất và mạng lưới bấc thấm, cho phép tối ưu hóa chiều cao đắp gia tải từ 3,0 mét đến 5,5 mét và kiểm soát tốc độ đắp không vượt quá 10 kPa mỗi tuần nhằm ngăn chặn sự cố trượt trồi thân đắp. Ban Quản lý dự án cùng Đơn vị tư vấn giám sát có trách nhiệm kiểm soát chỉ tiêu này liên tục trong suốt tiến độ thi công 12 tháng.

Thiết lập hệ thống quan trắc địa kỹ thuật tự động hóa bao gồm các bàn đo lún đa tầng và đầu đo áp lực nước lỗ rỗng dạng màng rung tại các độ sâu cách nhau mỗi 5 mét để thu thập dữ liệu với tần suất 24 giờ một lần, phục vụ công tác phân tích ngược và điều chỉnh tiến độ dỡ tải kịp thời. Nhà thầu thi công chuyên trách nền móng phối hợp cùng đơn vị quan trắc duy trì vận hành hệ thống này trong suốt 18 tháng gia tải và chờ lún.

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

Nội dung chuyên sâu của luận văn mang lại giá trị thực tiễn và học thuật cao cho 4 nhóm đối tượng chính:

Kỹ sư thiết kế địa kỹ thuật và công trình cảng biển: Tài liệu cung cấp công thức tính toán tường minh để xác định chính xác thông số vùng xáo động và lựa chọn loại bấc thấm phù hợp, phục vụ thiết kế các dự án lấn biển, san nền khu công nghiệp và nền đường ô tô đắp trên đất bùn dày đến 40 mét.

Chuyên viên mô phỏng số và phân tích phần tử hữu hạn: Nghiên cứu cung cấp phương pháp quy đổi tương đương từ bài toán 3D phức tạp sang bài toán biến dạng phẳng 2D bằng phần tử đặc, giúp tiết kiệm hơn 80% thời gian phân tích số và tối ưu hóa tài nguyên phần cứng máy tính trong các bài toán quy mô lớn.

Chủ đầu tư, Giám đốc ban quản lý dự án và Tư vấn giám sát hạ tầng: Giúp nắm bắt sâu sắc cơ chế tiêu tán áp lực nước lỗ rỗng và độ lún theo thời gian để hoạch định tiến độ thi công đắp nền cao trên 4 mét, kiểm soát rủi ro mất ổn định công trình và quản lý dòng tiền dự án hiệu quả.

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 - Địa kỹ thuật môi trường: Đây là nguồn tài liệu tham khảo chuẩn mực về lý thuyết biến dạng lớn, cung cấp hệ thống dữ liệu thực nghiệm phong phú từ 5 dự án quốc tế để phát triển các hướng nghiên cứu mở rộng trong tương lai.

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

Vùng xáo động do đóng bấc thấm ảnh hưởng như thế nào đến tốc độ cố kết của nền đất yếu? Quá trình ép quả ép chuyên dụng vào lòng đất làm xáo trộn cấu trúc đất sét trong phạm vi bán kính gấp 2,0 đến 4,5 lần bán kính quả ép. Sự xáo động này làm suy giảm hệ số thấm ngang từ 50% đến 75% và làm tăng độ nén lún thể tích, khiến thời gian đạt độ cố kết 90% bị kéo dài thêm từ 20% đến 35% so với lý thuyết thông thường.

Tại sao bắt buộc phải tính toán sự suy giảm khả năng thoát nước của bấc thấm theo độ sâu? Khi cắm bấc ở độ sâu lớn từ 25 mét đến 55 mét, áp lực nén hông từ đất xung quanh tăng cao vượt mức 250 kPa sẽ ép chặt vỏ bọc lọc vào các rãnh lõi nhựa. Thực tế tại Cảng mới Busan cho thấy khả năng thoát nước ở chân bấc giảm hơn 45% so với đầu bấc, khiến áp lực nước lỗ rỗng tiêu tán chậm hơn đáng kể nếu không xét đến yếu tố này.

Sự khác biệt cốt lõi giữa lý thuyết biến dạng lớn và biến dạng nhỏ trong bài toán bấc thấm là gì? Lý thuyết biến dạng nhỏ giả định chiều dày lớp đất không đổi nên chỉ phù hợp khi độ biến dạng đứng dưới 10%. Trong các dự án đất sét siêu mềm có độ lún từ 1,5 mét đến 3,2 mét, lý thuyết biến dạng lớn cập nhật liên tục sự thay đổi hình học và hiện tượng gập bấc, giúp giảm sai số dự báo lún từ 25% xuống dưới 5%.

Mô hình biến dạng phẳng 2D sử dụng phần tử đặc có ưu điểm gì vượt trội so với các phần tử thoát nước 1D? Mô hình phần tử đặc 2D trực tiếp gán các đặc tính hình học tương đương và chỉ tiêu cơ lý của bấc cũng như vùng xáo động vào lưới phần tử. Phương pháp này giúp mô phỏng chính xác đồng thời độ lún và sự tiêu tán áp lực nước lỗ rỗng với độ tin cậy tương đương mô hình 3D nhưng giảm hơn 80% thời gian tính toán.

Làm thế nào để xác định chính xác các thông số đầu vào cho mô hình phi tuyến trong thực tế công trường? Kỹ sư có thể kết hợp thí nghiệm xuyên tĩnh đo áp lực nước lỗ rỗng CPTU với thí nghiệm nén lún Oedometer để xác định chỉ số nén và chỉ số thấm. Trong thực tế thiết kế, tỷ lệ bán kính vùng xáo động trên bán kính quả ép thường được lựa chọn trong khoảng 2,0 đến 3,5 dựa trên kết quả phân tích ngược từ các ô đắp thử nghiệm.

Kết luận

  • Xây dựng thành công mô hình giải tích tích hợp phân bố phi tuyến của hệ số thấm và tính nén lún trong vùng xáo động.
  • Ứng dụng xuất sắc lý thuyết biến dạng lớn để mô phỏng hiện tượng suy giảm khả năng thoát nước theo thời gian và độ sâu.
  • Đề xuất phương pháp hệ số thấm đứng tương đương và mô hình biến dạng phẳng 2D phần tử đặc với sai số kiểm chứng dưới 6%.
  • Xác thực thành công độ tin cậy trên 5 dự án đắp thử nghiệm thực tế tại Nhật Bản, Hàn Quốc, Úc, Trung Quốc và Malaysia.
  • Tối ưu hóa các thông số thiết kế giúp rút ngắn từ 15% đến 25% thời gian xử lý nền và tiết kiệm đáng kể chi phí đầu tư.

Đóng góp lớn nhất của luận văn là hoàn thiện cơ sở khoa học và cung cấp bộ công cụ tính toán số tin cậy cho kỹ thuật xử lý nền đất yếu bằng bấc thấm PVD dưới tải trọng đắp lớn. Trong lộ trình 12 đến 24 tháng tới, hướng nghiên cứu cần tiếp tục mở rộng cho các bài toán kết hợp gia tải hút chân không và tải trọng động giao thông. Các đơn vị tư vấn thiết kế và nhà thầu hạ tầng hãy chủ động ứng dụng mô hình giải tích phi tuyến và phương pháp 2D phần tử đặc này vào các dự án cảng biển và đường cao tốc ngay từ hôm nay để đảm bảo chất lượng công trình bền vững và an toàn tuyệt đối.