Phân Tích Cố Kết Nền Đất Yếu Có Gia Cố Bấc Thấm Xét Đến Độ Cản Thấm Dựa Vào Lý Thuyết Biến Dạng Lớn: Đột Phá Mới Trong Địa Kỹ Thuật Xây Dựng

Tóm tắt nghiên cứu

Trong kỹ thuật xử lý nền đất yếu phục vụ xây dựng hạ tầng giao thông và công trình thủy công, phương pháp cắm bấc thấm thẳng đứng (Prefabricated Vertical Drain - PVD) kết hợp gia tải trước đóng vai trò then chốt nhằm thúc đẩy quá trình thoát nước và tăng nhanh tốc độ cố kết. Tuy nhiên, các phương pháp dự báo lún kinh điển thường bỏ qua hiện tượng biến dạng lớn (large-strain) và xem khả năng thoát nước của bấc thấm ($q_w$) là một đại lượng không đổi.

Đề tài nghiên cứu khoa học do TS. Nguyễn Bá Phú chủ nhiệm (Khoa Kỹ thuật Xây dựng, Trường Đại học Công nghiệp TP.HCM) đã giải quyết bài toán cốt lõi: Làm thế nào để mô hình hóa chính xác ứng xử cố kết và dự báo độ lún của nền đất yếu gia cố PVD khi xét đồng thời sự suy giảm phi tuyến của khả năng thoát nước (độ cản thấm) theo chiều sâu, theo thời gian và hiệu ứng biến dạng lớn?

Bằng cách thiết lập hệ phương trình vi phân cố kết đối xứng trục trong hệ tọa độ đối lưu (Convective coordinates) kết hợp tọa độ Lagrangian, nghiên cứu đã tích hợp sự biến thiên phi tuyến của hệ số rỗng ($e$), hệ số thấm ($k_h$) và chỉ số nén ($C_c$). Kết quả nghiên cứu chỉ ra rằng:

  1. Độ cản thấm tăng phi tuyến theo chiều sâu và thời gian làm chậm trễ đáng kể tiến trình cố kết, đặc biệt rõ nét ở các tầng đất yếu có bề dày lớn.
  2. Khả năng thoát nước ban đầu $q_{w0} = 100\text{ m}^3/\text{năm}$ được xác định là ngưỡng tới hạn mà trên mức này, ảnh hưởng của độ cản thấm trở nên không đáng kể.
  3. Phương trình dự báo lún mới được đề xuất cho phép dự báo chính xác độ lún theo thời gian ngay từ giai đoạn đầu thi công đắp tải, khắc phục hoàn toàn nhược điểm của các phương pháp truyền thống như Asaoka hay Hyperbolic.

Bối cảnh và tầm quan trọng

Thực trạng kỹ thuật và khoảng trống nghiên cứu

Tại Việt Nam, các công trình hạ tầng giao thông trọng điểm, khu công nghiệp và đô thị ven biển tập trung mật độ cao tại khu vực Đồng bằng sông Cửu Long và Đồng bằng sông Hồng. Đặc điểm địa chất công trình tại các vùng này đặc trưng bởi tầng đất bùn, sét mềm bão hòa nước có chiều dày lớn (thay đổi phổ biến từ 8 đến 40 m), hệ số rỗng cao, khả năng chịu tải cực kỳ thấp và tính nén lún phức tạp.

Hiện nay, công tác tính toán thiết kế PVD tại Việt Nam chủ yếu áp dụng tiêu chuẩn TCVN 9355:2012 và 22 TCN 262-2000, vốn dựa trên nền tảng lý thuyết giải tích của Hansbo (1981). Lời giải cổ điển này giả định:

  • Khả năng thoát nước của bấc thấm $q_w$ là hằng số dọc theo chiều sâu.
  • Nền đất chịu biến dạng nhỏ (small-strain), các chỉ tiêu cơ lý như hệ số thấm và tính nén lún không đổi trong suốt thời gian cố kết.

Tuy nhiên, trong thực tế thi công:

  • Áp lực hông từ đất nền: Tăng dần theo độ sâu, ép chặt vỏ bọc vào lõi dẫn, làm co hẹp tiết diện dòng chảy.
  • Tắc nghẽn cơ học: Các hạt sét siêu mịn thâm nhập qua lớp vải địa kỹ thuật, lắng đọng trong rãnh dẫn nước.
  • Hiệu ứng uốn gập (crimping): Khi nền đất lún lớn (chiều dày nén lún có thể đạt 15–20%), bấc thấm bị uốn cong, gãy khúc làm giảm đột ngột khả năng dẫn nước.

Sự suy giảm này khiến khả năng thoát nước thực tế của bấc giảm mạnh theo cả chiều sâu ($z$) lẫn thời gian cố kết ($t$). Việc bỏ qua các yếu tố này và bỏ qua biến dạng lớn dẫn đến sai số nghiêm trọng trong dự báo lún thực tế (thực tế tại nhiều dự án quốc tế ghi nhận độ lệch giữa tính toán và quan trắc lên tới 200% – 600%).

Tính thời sự và ý nghĩa ứng dụng

Khống chế độ lún cố kết còn lại (độ lún dư) sau khi đưa công trình vào khai thác là yêu cầu kỹ thuật sống còn để tránh hiện tượng lún lệch mố cầu, hư hỏng mặt đường cao tốc và lãng phí chi phí duy tu. Nghiên cứu này cung cấp một công cụ giải tích - số chuẩn xác, mang tính thời sự cao, hỗ trợ đắc lực cho các nhà tư vấn thiết kế và quản lý dự án trong việc tối ưu hóa sơ đồ cắm bấc, rút ngắn thời gian gia tải và kiểm soát an toàn công trình.

Phương pháp luận và cách tiếp cận nghiên cứu

Đề tài triển khai cách tiếp cận toàn diện kết hợp giữa mô hình hóa giải tích lý thuyết, mô phỏng số và kiểm chứng thực nghiệm hiện trường:

1. Thiết lập hình học mô hình Unit Cell và hệ tọa độ biến dạng lớn

Nghiên cứu khảo sát một khối đất lăng trụ tròn đại diện (unit cell) có bán kính $r_e$, ở tâm bố trí bấc thấm có bán kính tương đương $r_w$, xung quanh bấc là vùng đất bị xáo trộn do thi công (smear zone) có bán kính $r_s$ và hệ số thấm $k_s$.

Để giải quyết bài toán biến dạng lớn khi chiều cao của khối đất thay đổi liên tục theo thời gian $\xi_0(t) = H - S_0(t)$, nghiên cứu tích hợp hai hệ tọa độ:

  • Hệ tọa độ nguyên thủy Lagrangian ($a$).
  • Hệ tọa độ chuyển dịch Convective ($\xi$), tuân theo quan hệ biến dạng Gibson: $$\frac{\partial \xi}{\partial a} = \frac{1 + e}{1 + e_0}$$

2. Mô hình hóa sự suy giảm phi tuyến của khả năng thoát nước ($q_w$)

Khả năng thoát nước được mô tả qua hàm phân bố phi tuyến đồng thời theo chiều sâu và thời gian: $$q_w(t, \xi) = q_{w0} \left(1 - \alpha \frac{\xi}{\xi_0(t)}\right)^2 e^{-\beta t}$$

  • Trong đó: $q_{w0}$ là khả năng thoát nước ban đầu tại mặt đất ($z = 0$); $\alpha$ là hệ số suy giảm theo chiều sâu ($0 \le \alpha \le 1$); $\beta$ là hệ số suy giảm theo thời gian ($day^{-1}$).

3. Tích hợp phi tuyến cơ lý đất nền và lời giải phương trình vi phân

Nghiên cứu liên kết chặt chẽ sự suy giảm hệ số rỗng với ứng suất hữu hiệu và hệ số thấm:

  • Quan hệ nén lún: $e = e_0 - C_c \log\left(\frac{\sigma'v}{\sigma'{v0}}\right)$
  • Quan hệ biến thiên hệ số thấm: $\log(k_h) = \log(k_{h0}) - \frac{e_0 - e}{C_k}$

Từ điều kiện cân bằng dòng thấm và định luật Darcy, phương trình vi phân phi tuyến bậc cao mô tả sự biến thiên của ứng suất hữu hiệu $\sigma'_v$ theo thời gian được thiết lập và giải số bằng phương pháp sai phân hữu hạn (FDM) kết hợp giải thuật Euler.

Thông số mô hình Ký hiệu Ý nghĩa vật lý / Nguồn gốc
Bán kính tương đương PVD $r_w$ Quy đổi từ tiết diện dẹt $a \times b$ theo quy tắc Hansbo
Hệ số xáo trộn $s = r_s / r_w$ Mức độ mở rộng của vùng đất bị xáo trộn do cắm bấc
Tỷ số thấm xáo trộn $k_h / k_s$ Phản ánh mức độ suy giảm tính thấm sát vách bấc thấm
Hệ số suy giảm chiều sâu $\alpha$ Mức độ bẹp tiết diện bấc do áp lực hông tăng theo độ sâu
Hệ số suy giảm thời gian $\beta, A_3$ Mức độ tắc nghẽn lõi bấc do hạt mịn và uốn nếp theo thời gian

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

1. Mức độ suy giảm $q_w$ theo chiều sâu làm trì hoãn mạnh mẽ tiến trình cố kết

Phân tích cho thấy khi hệ số $\alpha$ tăng từ $0$ (lý tưởng - Hansbo) lên $1.0$ (thoát nước bị triệt tiêu hoàn toàn ở đáy bấc), tốc độ cố kết của nền đất suy giảm nghiêm trọng.

  • Tại các lớp đất sát mặt đất ($z = 5\text{ m}$), ảnh hưởng của $\alpha$ tương đối nhỏ.
  • Tại các tầng đất sâu ($z = 15\text{ m} \div 20\text{ m}$), độ cố kết tại thời điểm 200 ngày có thể chênh lệch từ 30% đến 45% giữa trường hợp bấc lý tưởng và bấc bị cản thấm thực tế.
Mức độ ảnh hưởng của cản thấm theo thời gian:
Độ chênh lệch ΔU = U(α=0) - U(α=1) (%)

2. Xác lập ngưỡng khả năng thoát nước tới hạn $q_{w0} = 100\text{ m}^3/\text{năm}$

Nghiên cứu đã định lượng hóa mức độ chênh lệch độ cố kết trung bình $\Delta U = U(\alpha=0) - U(\alpha=1)$ dưới các dải lưu lượng bấc khác nhau:

  • Khi $q_{w0} = 10\text{ m}^3/\text{năm}$, $\Delta U$ đạt đỉnh trên 13%, gây trễ hạn cố kết nghiêm trọng.
  • Khi $q_{w0} \ge 100\text{ m}^3/\text{năm}$, độ chênh lệch $\Delta U < 2%$. Điều này khẳng định rằng khi bấc thấm đạt khả năng thoát nước từ $100\text{ m}^3/\text{năm}$ trở lên, độ cản thấm trở nên không đáng kể đối với tốc độ cố kết tổng thể của nền.

3. Ưu thế vượt trội của lý thuyết biến dạng lớn tại các dự án thực tế

Kiểm chứng tại hai đại công trình thực tế có chiều dày đất yếu lớn đã chứng minh tính chuẩn xác của mô hình đề xuất:

  • Dự án mở rộng Cảng Busan (Hàn Quốc - Chiều sâu bấc 46 m): Nền đường đắp cao 9 m trên tầng sét dày 46 m. Phương pháp cổ điển của Hansbo (1981) và phương pháp biến dạng nhỏ của Deng et al. (2013) đều dự báo độ lún vượt quá xa số liệu thực tế (dự báo lún tức thời quá nhanh). Ngược lại, đường cong dự báo theo mô hình biến dạng lớn phi tuyến của đề tài bám sát hoàn hảo chuỗi số liệu quan trắc lún suốt 400 ngày.
  • Dự án Sân bay Saga (Nhật Bản): Mô hình biến dạng lớn kết hợp hệ số suy giảm thời gian $\beta = 0.18\text{ day}^{-1}$ cho kết quả trùng khớp với độ lún đo đạc, trong khi các mô hình bỏ qua suy giảm $q_w$ đều gây sai số lớn.

4. Thiết lập công thức dự báo lún thực dụng theo thời gian

Dựa trên quan hệ phi tuyến giữa độ cản thấm $L(t)$ và nhân tố thời gian $T_h$, tác giả đã thiết lập phương trình dự báo lún hữu hạn: $$S(t) = S_{final} \cdot \left[ 1 - \exp\left( - \frac{8 \ln\left(\frac{1 - D_0 e^{-a_3 T_h}}{1 - D_0}\right)}{a_3 F(n)} \right) \right]$$ Phương trình này cho phép sử dụng chuỗi số liệu quan trắc độ lún $S(t_i)$ ở các giai đoạn đắp ban đầu để xác định tham số $F$, từ đó ngoại suy chính xác tiến trình lún và độ lún cuối cùng mà không bị phụ thuộc chủ quan vào việc chọn khoảng thời gian $\Delta t$ như phương pháp Asaoka.

Đóng góp khoa học và giá trị thực tiễn

Đóng góp về mặt học thuật

  • Bổ sung lý thuyết cố kết phi tuyến: Đề tài là một trong số rất ít công trình tích hợp thành công cả 3 yếu tố phi tuyến phức tạp trong cơ học đất: hiệu ứng biến dạng lớn (Lagrangian/Convective), sự thay đổi động của tính thấm/nén lún, và quy luật suy giảm phi tuyến không gian - thời gian của giếng thoát nước.
  • Đổi mới phương pháp tính toán: Xây dựng khung thuật toán số ổn định, cho phép giải quyết bài toán cố kết đối xứng trục có xét đến vùng xáo trộn và cản thấm mà không cần đơn giản hóa các điều kiện biên.

Ứng dụng thực tiễn trong kỹ thuật xây dựng

  • Chỉ dẫn lựa chọn vật liệu: Cung cấp cơ sở khoa học định lượng để kỹ sư lựa chọn chỉ tiêu kỹ thuật của bấc thấm (yêu cầu $q_w > 100\text{ m}^3/\text{năm}$ đối với công trình xử lý đất yếu sâu) nhằm triệt tiêu rủi ro tắc nghẽn và cản thấm.
  • Tối ưu hóa thiết kế và thi công: Hỗ trợ tính toán khoảng cách cắm bấc và thời gian gia tải tối ưu, giúp chủ đầu tư rút ngắn tiến độ chờ cố kết và tiết kiệm hàng tỷ đồng chi phí xử lý nền.
  • Đề xuất hoàn thiện tiêu chuẩn: Đóng vai trò là tài liệu khoa học nền tảng để xem xét, bổ sung và cập nhật tiêu chuẩn quốc gia TCVN 9355:2012 đối với các bài toán xử lý nền đất yếu có chiều dày lớn trên 15 m.

Đối tượng quan tâm và lợi ích ứng dụng

Nội dung nghiên cứu mang lại giá trị thiết thực cho nhiều nhóm đối tượng trong ngành xây dựng công trình:

  • Kỹ sư tư vấn thiết kế Địa kỹ thuật: Có thêm công cụ tính toán giải tích tin cậy để dự báo lún chính xác cho các dự án cao tốc, sân bay, bến cảng trên nền đất yếu sâu; khắc phục tình trạng thiết kế thiếu an toàn hoặc dư thừa lãng phí do dùng công thức cũ.
  • Nhà thầu thi công & Tư vấn giám sát: Ứng dụng công thức dự báo lún đề xuất để phân tích ngược số liệu quan trắc lún bàn đo lún tầng sâu, chủ động điều chỉnh kế hoạch đắp tải theo từng giai đoạn và xác định chính xác thời điểm dỡ tải gia tải trước.
  • Các nhà nghiên cứu, Giảng viên & Học viên Sau đại học: Tài liệu tham khảo học thuật giá trị cao trong lĩnh vực Địa kỹ thuật công trình, cơ học đất nâng cao, mô hình hóa toán học và phương pháp số trong địa kỹ thuật.
  • Cơ quan quản lý nhà nước & Ban soạn thảo tiêu chuẩn: Cung cấp dữ liệu thực nghiệm và lý thuyết xác thực để nâng cấp hệ thống quy chuẩn, tiêu chuẩn kỹ thuật xây dựng quốc gia.

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

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Nghiên cứu đã chỉ ra rằng trong nền đất yếu dày, khả năng thoát nước của bấc thấm ($q_w$) suy giảm mạnh theo dạng phi tuyến theo chiều sâu và thời gian, làm chậm đáng kể tiến trình cố kết ở các lớp đất sâu. Đồng thời, nghiên cứu xác lập ngưỡng $q_{w0} = 100\text{ m}^3/\text{năm}$ – mức lưu lượng tối thiểu giúp loại bỏ hoàn toàn tác động tiêu cực của độ cản thấm.

2. Phương pháp đề xuất có điểm gì ưu việt hơn phương pháp Hansbo truyền thống?

Lời giải Hansbo (1981) coi $q_w$ là hằng số và áp dụng lý thuyết biến dạng nhỏ, dẫn đến việc đánh giá quá nhanh tốc độ cố kết và sai lệch lớn về độ lún ở đất yếu dày. Phương pháp mới xét đến biến dạng lớn, sự giảm chiều dày lớp đất thực tế, sự suy giảm hệ số thấm $k_h$, chỉ số nén $C_c$ và sự suy giảm phi tuyến của $q_w$.

3. Kết quả nghiên cứu có thể khái quát hóa cho mọi loại nền đất yếu không?

Phương pháp hoàn toàn có thể áp dụng phổ quát cho các loại đất dính bão hòa nước (sét mềm, bùn sét, sét pha bùn) có tính nén lún cao và chiều dày lớn. Khi áp dụng cho từng loại đất cụ thể, kỹ sư chỉ cần xác định chính xác các thông số cơ bản ($e_0, C_c, C_k, k_{h0}$) thông qua các thí nghiệm nén cố kết tiêu chuẩn.

4. Vì sao phương pháp dự báo lún đề xuất lại tốt hơn phương pháp Asaoka và Hyperbolic?

Phương pháp Asaoka phụ thuộc nhiều vào việc lựa chọn bước thời gian $\Delta t$ và thường cho độ lún cuối cùng nhỏ hơn thực tế. Phương pháp Hyperbolic lại phụ thuộc lớn vào điểm gốc thời gian ban đầu. Phương trình đề xuất của đề tài xây dựng dựa trên bản chất vật lý của độ cản thấm $L(t)$, giúp xây dựng chính xác toàn bộ đường cong lún theo thời gian từ số liệu quan trắc sớm.

5. Trong thực tế thiết kế, làm thế nào để xác định các hệ số suy giảm $\alpha$ và $\beta$?

Hệ số $\alpha$ (suy giảm theo chiều sâu) và $\beta$ (suy giảm theo thời gian) có thể xác định thông qua thí nghiệm đo khả năng thoát nước của bấc thấm dưới áp lực hông tăng dần trong buồng áp lực chuyên dụng, hoặc thông qua phương pháp phân tích ngược (back-analysis) từ kết quả quan trắc của các ô đắp thử nghiệm (test embankment).

Kết luận

Đề tài nghiên cứu khoa học của TS. Nguyễn Bá Phú đã tạo nên một bước tiến quan trọng trong lĩnh vực Địa kỹ thuật xây dựng tại Việt Nam. Bằng việc giải quyết trọn vẹn sự tương tác phức tạp giữa lý thuyết biến dạng lớn và quy luật cản thấm phi tuyến của bấc thấm, nghiên cứu không chỉ làm sáng tỏ cơ chế cố kết thực tế của nền đất yếu dày mà còn cung cấp phương pháp dự báo lún có độ chính xác vượt trội.

Khuyến nghị các đơn vị thiết kế, chủ đầu tư các dự án giao thông trọng điểm nhanh chóng ứng dụng các kết quả và công thức dự báo đề xuất vào thực tiễn, góp phần nâng cao chất lượng công trình và tối ưu hóa hiệu quả kinh tế hạ tầng quốc gia.