Tổng quan nghiên cứu
Khu vực Đồng bằng sông Cửu Long có cấu tầng địa chất đặc trưng với lớp bùn sét yếu bão hòa nước dày từ 15 m đến hơn 35 m, dẫn đến các công trình đắp nền đường thường ghi nhận độ lún tích lũy thực tế từ 50 cm đến 120 cm sau 3 đến 5 năm đưa vào khai thác. Vấn đề nghiên cứu trọng tâm xuất phát từ việc tính toán độ lún truyền thống dựa trên lý thuyết cố kết một chiều của Terzaghi thường giả định tải trọng phân bố rộng vô hạn, bỏ qua hiện tượng thoát nước theo phương ngang và biến dạng hông của khối đất. Sự đơn giản hóa này dẫn đến sai số lớn khi dự báo tốc độ lún theo thời gian và hoàn toàn không đánh giá được hiện tượng lún lệch giữa tim đường và mái taluy.
Mục tiêu cụ thể của công trình là thiết lập khung tính toán dự báo độ lún và độ lún lệch theo thời gian trên cơ sở bài toán cố kết thấm hai chiều kết hợp chia lớp phân tố, có xét đến tính nén ép của nước lỗ rỗng và độ lún ngắn hạn tức thời. Phạm vi nghiên cứu ứng dụng trực tiếp cho dự án nâng cấp, mở rộng Quốc lộ 1A đoạn Mỹ Thuận – Cần Thơ với chiều dài tuyến 31 km từ Km 2030 đến Km 2061, đi qua địa phận tỉnh Vĩnh Long và thành phố Cần Thơ, được hoàn thành đánh giá vào tháng 12 năm 2015.
Nghiên cứu mang ý nghĩa thực tiễn to lớn khi nâng cao độ chính xác trong dự báo độ lún ban đầu lên khoảng 15% đến 25%, giúp các đơn vị tư vấn thiết kế kiểm soát độ dốc ngang mặt đường và giảm thiểu hơn 30% nguy cơ nứt gãy mặt đường hoặc trượt sạt mái taluy do lún không đều trong quá trình khai thác lâu dài.
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ự tích hợp giữa lý thuyết cố kết thấm kinh điển Terzaghi, lý thuyết đàn hồi tuyến tính áp dụng cho môi trường đất và lý thuyết cố kết phẳng hai chiều Florin. Phương pháp tiếp cận mở rộng mô hình biến dạng liên tục nhiều pha gồm hạt rắn, chất lỏng và bọt khí, có kể đến tính hòa tan của khí theo định luật Henry với hệ số hòa tan H = 0,02.
Bốn khái niệm cốt lõi tạo nên nền tảng mô hình bao gồm:
- Áp lực nước lỗ rỗng thặng dư và nguyên lý ứng suất hữu hiệu tổng quát: Biểu diễn mối quan hệ giữa tổng ứng suất, ứng suất hữu hiệu khung cốt đất và áp lực nước lỗ rỗng.
- Phân tách chuyển vị tổng cộng thành hai thành phần độc lập: Chuyển vị do biến dạng hình dạng biến thiên theo mô đun trượt và chuyển vị do biến dạng thể tích biến thiên theo mô đun nén khối.
- Mô đun biến dạng thể tích tổng thể và mô đun biến dạng cắt ban đầu: Đánh giá trạng thái ứng suất biến dạng tức thời của nền đất bão hòa ngay sau khi gia tải khi nước lỗ rỗng chưa kịp thoát ra ngoài.
- Hệ số thấm dị hướng: Phản ánh tỷ số hệ số thấm phương ngang so với phương đứng, thường dao động trong khoảng từ 2 đến 5 lần trong trầm tích sét ĐBSCL.
Phương pháp nghiên cứu
Nguồn dữ liệu tính toán được thu thập từ 12 hố khoan khảo sát địa chất công trình có chiều sâu từ 25 m đến 40 m dọc tuyến Quốc lộ 1A đoạn Vĩnh Long – Cần Thơ, kết hợp các kết quả thí nghiệm nén cố kết oedometer tiêu chuẩn và nén ba trục không thoát nước CU.
Phương pháp phân tích chủ đạo là kết hợp giải tích toán học bằng chuỗi Fourier kép hai chiều với phương pháp chia nền đất thành các lớp phân tố mỏng có chiều dày từ 0,5 m đến 1,0 m. Thuật toán tích phân số được lập trình tự động hóa trên môi trường phần mềm tính toán Mathcad với số lượng phần tử chuỗi mở rộng trên 10 phần tử nhằm đảm bảo độ hội tụ và độ chính xác cao.
Lý do lựa chọn phương pháp giải tích lớp phân tố thay vì lấy giá trị trung bình điểm là vì việc tích phân giải tích thể tích vùng áp lực nước lỗ rỗng hai chiều dưới tải trọng hình thang của nền đắp rất phức tạp; chia lớp phân tố cho phép xác định chính xác ứng suất đẳng hướng và gradient thủy lực tại từng tọa độ không gian riêng biệt. Timeline nghiên cứu được thực hiện đồng bộ trong 5 tháng, từ ngày 06/07/2015 đến ngày 04/12/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à tính toán kiểm chứng trên nền đường đắp Quốc lộ 1A đoạn Vĩnh Long – Cần Thơ mang lại các phát hiện kỹ thuật nổi bật:
- Thành phần độ lún ban đầu tức thời do biến dạng hình dạng chiếm tỷ trọng từ 18% đến 25% trong tổng độ lún toàn bộ, chứng minh rằng biến dạng hông đóng vai trò then chốt ngay sau khi đắp nền đường mà bài toán một chiều Terzaghi đã bỏ sót.
- Tốc độ tiêu tán áp lực nước lỗ rỗng thặng dư theo mô hình hai chiều diễn ra nhanh hơn từ 30% đến 45% trong giai đoạn 100 đến 300 ngày đầu so với mô hình một chiều, nguyên nhân do sự xuất hiện của dòng thấm ngang hướng ra hai bên mái dốc với hệ số thấm phương ngang lớn gấp 2 đến 3 lần phương đứng.
- Giá trị độ lún ổn định cuối cùng tại tâm diện gia tải giữa phương pháp hai chiều theo lớp phân tố và phương pháp Terzaghi có sự tiệm cận rất cao, mức chênh lệch tuyệt đối ghi nhận dưới 5%.
- Hiện tượng lún lệch giữa tim đường và biên gia tải tại mép mái taluy ở tọa độ x = 9 m phát triển mạnh nhất trong 180 ngày đầu thi công với độ chênh lún lên đến 12 cm đến 18 cm, sau đó duy trì độ lệch ổn định theo thời gian dài hạn.
Thảo luận kết quả
Cơ chế của sự khác biệt giữa các phương pháp xuất phát từ sự phân bố lại ứng suất cắt và đường thấm nước lỗ rỗng. Khi tải trọng đắp dạng hình thang tác dụng, ứng suất nén không phân bố đều theo độ sâu mà suy giảm theo quy luật không gian hai chiều, đồng thời nước lỗ rỗng có xu hướng thoát nhanh qua biên mái dốc thoáng.
Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị quan hệ hệ số rỗng và logarit áp lực nén, biểu đồ đẳng mức áp lực nước lỗ rỗng thặng dư qua các mốc thời gian 30 ngày, 90 ngày, 360 ngày và bảng so sánh độ lún tổng hợp. Kết quả chỉ ra rằng mô hình hai chiều tính theo điểm trung bình thường đánh giá thấp độ lún ở các lớp đất mặt trên cùng khoảng 10% đến 15%, trong khi mô hình giải tích hai chiều theo lớp phân tố phản ánh chính xác từng đường cong lún thực tế. So với các công bố của những tác giả nghiên cứu về bài toán cố kết phẳng trước đây, phương pháp lớp phân tố kết hợp tính nén ép của pha lỏng đem lại lời giải liên tục và ổn định hơn.
Đề xuất và khuyến nghị
Dựa trên kết quả nghiên cứu, bốn giải pháp chiến lược được đề xuất cho công tác thiết kế và thi công nền đường trên đất yếu:
- Ứng dụng giải thuật cố kết hai chiều theo lớp phân tố vào phần mềm tính toán thiết kế kỹ thuật nhằm kiểm soát sai số dự báo lún dưới 8%, mục tiêu hoàn thành tích hợp trong 6 tháng do các đơn vị tư vấn thiết kế giao thông chủ trì thực hiện.
- Thiết lập hệ thống quan trắc áp lực nước lỗ rỗng đa tầng kết hợp mốc đo lún bề mặt định kỳ 7 ngày một lần trong 6 tháng đầu đắp nền, đảm bảo tốc độ lún không vượt quá ngưỡng an toàn 2 mm/ngày do Ban Quản lý Dự án và tư vấn giám sát vận hành.
- Áp dụng giải pháp bấc thấm kết hợp gia tải trước theo từng giai đoạn đắp nhằm rút ngắn thời gian đạt độ cố kết 80% từ 24 tháng xuống còn khoảng 8 đến 10 tháng, do các nhà thầu thi công hạ tầng triển khai tại hiện trường.
- Ban hành hướng dẫn kỹ thuật chuyên ngành về đánh giá biến dạng hông và lún lệch cho công trình giao thông khu vực Đồng bằng sông Cửu Long, đặt mục tiêu hoàn thiện khung tiêu chuẩn trong vòng 12 tháng do Bộ Giao thông Vận tải và các viện nghiên cứu chuyên ngành phụ trách.
Đố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 giao thông: Ứng dụng công thức giải tích hai chiều và các hàm thuật toán Mathcad để lập bảng tính dự báo lún lệch chính xác cho các dự án đắp nền đường, mố cầu và đê bao thủy lợi.
- Cán bộ quản lý dự án và kỹ sư tư vấn giám sát: Nắm vững quy luật tiêu tán áp lực nước lỗ rỗng để xây dựng biểu đồ tiến độ đắp gia tải hợp lý, xác định chính xác thời điểm dỡ tải gia tải trước và thảm bê tông nhựa mặt đường.
- Giảng viên và nghiên cứu sinh chuyên ngành Địa kỹ thuật Xây dựng: Sử dụng làm tài liệu tham khảo chuyên sâu về phương pháp giải phương trình vi phân đạo hàm riêng cố kết thấm phẳng có xét điều kiện biên phức tạp.
- Học viên cao học và sinh viên ngành Xây dựng Cầu đường: Tiếp cận phương pháp luận khoa học trong việc kết hợp số liệu khảo sát thực địa với mô hình toán học giải tích để giải quyết các bài toán cơ học đất thực tế.
Câu hỏi thường gặp
-
Tại sao phương pháp cố kết một chiều Terzaghi lại đánh giá chưa chuẩn xác thời gian lún của nền đường đắp? Phương pháp một chiều giả định chiều rộng tải trọng đắp là vô hạn và nước chỉ thoát theo phương đứng. Thực tế nền đường có bề rộng hữu hạn, nước lỗ rỗng thoát theo cả phương ngang ra hai bên taluy với hệ số thấm phương ngang cao gấp 2 đến 5 lần, làm tốc độ lún ban đầu diễn ra nhanh hơn đáng kể.
-
Việc tính toán theo lớp phân tố mang lại ưu điểm vượt trội gì so với phương pháp tính điểm trung bình? Phương pháp tính theo điểm trung bình giả định ứng suất nén đều trên toàn lớp, dễ gây sai lệch từ 10% đến 20% tại các vị trí biên. Phương pháp lớp phân tố chia nhỏ địa tầng thành từng dải từ 0,5 m đến 1,0 m, cho phép tính đúng gradient áp lực nước lỗ rỗng và lún lệch tại từng tọa độ không gian.
-
Tính nén ép của nước lỗ rỗng và khí hòa tan ảnh hưởng ra sao đến biến dạng nền đất? Trong đất yếu bão hòa chưa hoàn toàn, bọt khí và nước chịu nén ép tức thời khi chịu tải, tạo ra mô đun biến dạng thể tích ban đầu. Việc xét tính nén ép giúp tính toán chính xác thành phần lún tức thời ngay sau khi gia tải mà các mô hình truyền thống bỏ qua.
-
Hệ số thấm dị hướng tác động thế nào đến tiến độ cố kết của nền đất yếu? Trầm tích sét Đồng bằng sông Cửu Long có hệ số thấm ngang lớn hơn thấm đứng từ 2 đến 5 lần. Tỷ số dị hướng này thúc đẩy quá trình tiêu tán áp lực nước lỗ rỗng nhanh hơn từ 30% đến 45% trong 300 ngày đầu, giúp công trình sớm đạt độ ổn định mong muốn.
-
Thuật toán trong luận văn có thể áp dụng cho các loại đất nền và tải trọng khác không? Thuật toán giải tích chuỗi Fourier hai chiều hoàn toàn có thể tùy biến cho các dạng tải trọng đắp phức tạp, tải trọng móng băng hoặc nền san lấp công nghiệp, chỉ cần hiệu chỉnh thông số mô đun biến dạng, bề dày lớp phân tố và điều kiện biên thoát nước phù hợp.
Kết luận
- Luận văn đã giải quyết hoàn chỉnh bài toán ước lượng độ lún và lún lệch theo thời gian thông qua mô hình cố kết thấm hai chiều kết hợp chia lớp phân tố.
- Xác định chính xác thành phần lún ban đầu tức thời do biến dạng hình dạng chiếm từ 18% đến 25% tổng độ lún công trình đắp trên nền đất yếu.
- Chứng minh tốc độ cố kết thực tế diễn ra nhanh hơn mô hình một chiều từ 30% đến 45% trong giai đoạn đầu nhờ xét đến tính dị hướng của hệ số thấm.
- Cung cấp công cụ tính toán tự động hóa trên Mathcad với độ chính xác cao, kiểm chứng thành công trên tuyến Quốc lộ 1A đoạn Vĩnh Long – Cần Thơ dài 31 km.
- Đề ra định hướng nghiên cứu tiếp theo trong 12 tháng tới nhằm mở rộng thuật toán cho bài toán cố kết phi tuyến tính và tích hợp quan trắc tự động thời gian thực.
Để nâng cao chất lượng thiết kế và phòng ngừa rủi ro lún nứt công trình, các kỹ sư và nhà quản lý dự án hãy áp dụng ngay phương pháp tính toán cố kết hai chiều theo lớp phân tố vào quy trình thẩm tra và giám sát địa kỹ thuật hiện trường.