Phân Tích Ổn Định Trượt Sâu Công Trình Đắp Trên Nền Đất Yếu Xử Lý Bấc Thấm Gia Tải Trước: Nghiên Cứu Thực Nghiệm và Đánh Giá Đa Phương Pháp

Tóm tắt nghiên cứu (Research Summary)

  • Câu hỏi nghiên cứu chính (Research Question): Làm thế nào để dự báo chính xác độ gia tăng sức kháng cắt không thoát nước ($\Delta S_u$) của nền đất sét rất yếu trong quá trình cố kết dưới tải trọng đắp, và phương pháp phân tích ổn định trượt sâu nào phản ánh sát thực tế, đảm bảo an toàn tối ưu cho thi công phân đợt?
  • Phương pháp tiếp cận (Methodology Snapshot): Nghiên cứu tích hợp dữ liệu quan trắc hiện trường (đo lún mặt $SP$, lún sâu $E$, áp lực nước lỗ rỗng $Piezometer$, chuyển vị ngang $Inclinometer$) kết hợp thí nghiệm hiện trường ($VST, CPTu$) và thí nghiệm trong phòng ($CU, UU, DST$). Phân tích ổn định trượt sâu được thực hiện song song bằng 4 cách tiếp cận: Phương pháp giải tích ứng suất tổng với $S_u$ dự báo (theo 22 TCN 244-98 và SHANSEP qua phần mềm Slope/W), mô hình phần tử hữu hạn FEM (Plaxis 2D với chức năng $\Phi/c$-reduction), phương pháp ứng suất có hiệu dựa trên áp lực nước lỗ rỗng thực đo, và đồ thị quan trắc thực nghiệm Wakita & Matsuo (1994).
  • Kết quả then chốt (Key Findings): Sức kháng cắt không thoát nước gia tăng tuyến tính với độ cố kết ($U$) và ứng suất có hiệu ($\sigma'_v$), với tương quan $S_u \approx (0,24 - 0,26)\sigma'_v$ đối với lớp sét trên và $0,31\sigma'_v$ cho lớp sét dưới. Hệ số an toàn ổn định ($Fs$) phân tích theo ứng suất có hiệu cao hơn phân tích theo ứng suất tổng khoảng 35%. Đồ thị Wakita & Matsuo phản ánh xu hướng nhạy bén nhất với biến dạng ngang thực tế, đóng vai trò cảnh báo sớm rủi ro mất ổn định.
  • Ý nghĩa thực tiễn (Implications): Cung cấp quy trình kiểm soát an toàn thi công nền đắp theo từng phân đoạn cho các dự án hạ tầng, cảng biển và khu công nghiệp trên nền đất rất yếu tại khu vực Đồng bằng sông Cửu Long và TP. Hồ Chí Minh.

Bối cảnh và tầm quan trọng (Context & Significance)

1. Hiện trạng kiến thức và thách thức kỹ thuật

Việc xây dựng các công trình trọng điểm như cảng biển, sân bay, bãi chứa hàng container và đường dẫn trên nền đất yếu chịu nén lún lớn là một thách thức kỹ thuật thường trực tại Việt Nam, đặc biệt là khu vực hạ lưu sông Sài Gòn – Nhà Bè (Hiệp Phước). Địa tầng đặc trưng tại đây thường có lớp bùn sét yếu bão hòa nước dày từ 15 đến hơn 30m, với sức kháng cắt không thoát nước ban đầu cực nhỏ ($S_{u0} = 8 \div 10\text{ kPa}$), chỉ số dẻo cao ($PI > 40%$) và chỉ số sệt $I_L > 1$.

Phương pháp xử lý nền bằng bấc thấm ($PVD$) kết hợp đắp gia tải trước theo phân đợt (staged construction) là giải pháp kinh tế và phổ biến nhất so với các giải pháp cọc xi măng đất ($CDM$) hay cọc cứng. Tuy nhiên, vấn đề sống còn trong công nghệ này là kiểm soát nguy cơ phá hoại trượt sâu của khối đắp trong suốt quá trình đắp tôn cao tải trọng.

       Khối đắp gia tải (Staged Embankment)
================================================ Mặt đất tự nhiên (+3.1m CDL)
  Lớp cát đệm thoát nước & Vải ĐKT gia cường
------------------------------------------------
  Lớp sét dẻo mềm đến chảy (1a, 1b: dày ~21.3m)
  - Cắm Bấc thấm (PVD) lưới tam giác s = 1.5m
  - Quá trình cố kết: Δu tiêu tán ➔ σ'v tăng ➔ Su tăng (ΔSu)
------------------------------------------------
  Lớp sét dẻo cứng (1c: dày ~11.7m)
------------------------------------------------
  Lớp cát chịu lực (Lớp 2)
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Cung trượt sâu nguy hiểm

2. Khoảng trống nghiên cứu (Research Gap)

Trong thực hành kỹ thuật, đa phần các đơn vị tư vấn thiết kế chỉ tính toán ổn định ban đầu hoặc sử dụng các công thức giải tích đơn giản mà thiếu sự kiểm chứng đối chiếu với quá trình gia tăng sức kháng cắt thực tế ($\Delta S_u$) theo thời gian cố kết. Sự chênh lệch lớn giữa các lý thuyết dự báo ($SHANSEP$ của Ladd 1991, quy trình 22 TCN 244-98) và ứng xử thực địa tạo ra khoảng trống kỹ thuật lớn, tiềm ẩn nguy cơ sạt trượt thảm khốc (tương tự sự cố tại Kè Cảng Năm Căn hay Cảng Interflour).

3. Tính thời sự và tác động lan tỏa

Đề tài nghiên cứu cấp Trường mã số T-KTXD-2013-55 do PGS. Trần Xuân Thọ làm chủ nhiệm (Bộ môn Địa cơ Nền móng – Đại học Bách khoa TP.HCM) thực hiện trên đối tượng thực tế là Cảng Container Trung tâm Sài Gòn (SPCT - Hiệp Phước) mang tính cấp thiết cao. Nghiên cứu giải quyết bài toán tối ưu hóa tiến độ gia tải, xác định chính xác chiều cao đắp tới hạn cho từng phân đoạn mà không cần hệ số an toàn cưỡng bức quá mức gây lãng phí chi phí xử lý.

Phương pháp nghiên cứu và cách tiếp cận (Methodology & Approach)

Nghiên cứu ứng dụng phương pháp tiếp cận đa chiều, tích hợp giữa lý thuyết cơ học đất hiện đại, khảo sát thí nghiệm thực địa và mô phỏng số học.

                  QUY TRÌNH PHÂN TÍCH TỔNG HỢP
[DỮ LIỆU THỰC ĐỊA & LAB]                     [DỰ BÁO CƠ LÝ THEO THỜI GIAN]
              [ĐÁNH GIÁ ỔN ĐỊNH TRƯỢT SÂU ĐA PHƯƠNG PHÁP]
(1) Ứng suất tổng   (2) Mô hình FEM    (3) Ứng suất có     (4) Đồ thị quan trắc
  Slope/W - Bishop    Plaxis 2D          hiệu (Piezometer)   Wakita & Matsuo
  dựa trên Su dự báo  Phi/c-reduction     Slope/W             (1994)

1. Thiết kế nghiên cứu và đối tượng thực nghiệm

Dự án được khảo sát tại khu vực Phase 1-1 thuộc bãi Container Cảng SPCT. Nền đất yếu bao gồm lớp bùn sét dày 33m chia làm 3 phân lớp ($1a$: 7,2m; $1b$: 14,1m; $1c$: 11,7m). Bấc thấm bản dẹt ($PVD$) được bố trí theo sơ đồ tam giác đều cự ly $s = 1,5\text{ m}$, kết hợp 3 lớp vải địa kỹ thuật gia cường cường độ $200\text{ kN/m}$ tại đáy khối đắp và đắp cát gia tải qua 4 giai đoạn đạt cao độ $+8,13\text{m CDL}$.

2. Hệ thống thu thập dữ liệu quan trắc và thí nghiệm

  • Quan trắc nội hiện trường: Lắp đặt ống đo chuyển vị ngang $Inclinometer$ (I1, I2, I15); bản đo lún bề mặt (SP1, SP2); đo lún sâu đa điểm (E1); và đầu đo áp lực nước lỗ rỗng dây rung Piezometer tự động (P01, P01A) bố trí theo 3 tầng độ sâu.
  • Thí nghiệm kiểm chứng tại các mốc thi công: Thí nghiệm cắt cánh hiện trường ($VST$) hiệu chỉnh theo Bjerrum (1972) và xuyên tĩnh điện có đo áp lực nước lỗ rỗng ($CPTu$) tại các giai đoạn đắp khác nhau (mốc khảo sát ban đầu và mốc kiểm chứng ngày 13/06/2007).

3. Phương pháp dự báo sức kháng cắt $\Delta S_u$

  1. Theo quy chuẩn 22 TCN 244-98: $$\Delta S_u = U \cdot \Delta\sigma_v \cdot \tan\phi_{cu}$$ Trong đó độ cố kết $U = S_t / S_{fc}$ với độ lún cuối cùng $S_{fc}$ được nội suy chính xác qua phương pháp đồ thị quan trắc Asaoka (1978).
  2. Theo mô hình SHANSEP (Ladd, 1991): $$\frac{S_u}{\sigma'_{v0}} = S \cdot (OCR)^m$$ Với các thông số thực nghiệm xác lập cho bùn sét Hiệp Phước là $S = 0,24$ và $m = 0,8$.

4. Kỹ thuật phân tích ổn định trượt sâu

  • Phương pháp giải tích Bishop đơn giản: Phân tích ứng suất tổng và ứng suất có hiệu thông qua phần mềm GeoStudio (Slope/W) để quét tìm mặt trượt tròn nguy hiểm nhất.
  • Phương pháp phần tử hữu hạn (FEM - Plaxis 2D): Ứng dụng mô hình Soft Soil Model cho nền sét yếu và Mohr-Coulomb cho cát đắp. Quá trình tiêu tán áp lực nước lỗ rỗng qua hệ PVD được mô phỏng tương đương theo lý thuyết thấm Chai & Miura (1977) kết hợp nghiệm giải tích Hansbo (1979). Hệ số ổn định được xác định bằng thuật toán suy giảm độ bền $\Phi/c\text{-reduction}$.
  • Phương pháp thực nghiệm Wakita & Matsuo (1994): Thiết lập tỷ số giữa chuyển vị ngang lớn nhất tại chân mái dốc ($\gamma_s$) và độ lún tại tâm khối đắp ($\delta$) để tra trực tiếp hệ số an toàn $Fs$ tức thời.

Phát hiện chính và kết quả thực nghiệm (Major Findings)

1. Quy luật gia tăng sức kháng cắt không thoát nước thực tế

Nghiên cứu chỉ ra mối tương quan chặt chẽ giữa sức kháng cắt không thoát nước ban đầu $S_{u0}$ và sự gia tăng của nó theo chiều sâu:

  • Lớp 1a và 1b ($z \le 20\text{ m}$): $S_{u0} = 6 + 1,05z\text{ (kPa)}$
  • Lớp 1c ($z > 20\text{ m}$): $S_{u0} = 27 + 2,70z\text{ (kPa)}$

Khi gia tải cố kết, tương quan giữa $S_u$ dự báo và ứng suất hiệu quả $\sigma'_v$ thể hiện rõ đặc tính cố kết thường:

  • Theo công thức SHANSEP: $S_u = 0,24\sigma'_v$
  • Theo quy chuẩn 22 TCN 244-98: $S_u = 0,26\sigma'_v$ (cho lớp 1a, 1b) và $S_u = 0,31\sigma'_v$ (cho lớp 1c).

Kết quả kiểm chứng tại hố khoan kiểm tra $VS02$ (ngày 13/06/2007) cho thấy phương pháp 22 TCN 244-98 cho giá trị dự báo sát với thực tế cắt cánh hiện trường nhất, trong khi SHANSEP có xu hướng dự báo thấp hơn khoảng 15%.

2. Tương quan hệ số mũi xuyên $N^*_k$ từ thí nghiệm CPTu

Giá trị hệ số nón xuyên hiệu chỉnh $N^*_k$ thiết lập từ kết quả đối chiếu giữa $CPTu$ và $VST$ đạt $10,86$ ở giai đoạn khảo sát ban đầu và duy trì ở mức $10,01$ sau quá trình gia tải nén cố kết (sai số giảm $8,5%$). Điều này khẳng định tính ổn định cao của tương quan $CPTu$ phục vụ đánh giá nhanh sức kháng cắt tại hiện trường.

3. Đánh giá đa chiều về hệ số an toàn trượt sâu ($Fs$)

Bảng tổng hợp giá trị hệ số an toàn tại 7 mốc thời gian thi công (từ điểm A đến điểm G tương ứng chiều cao đắp từ $3,2\text{m}$ đến $7,8\text{m}$):

Điểm nghiên cứu Chiều cao đắp $H_t$ (m) Thời gian dồn (ngày) $Fs_1$ (Slope/W - $S_u$ dự báo) $Fs_2$ (Plaxis FEM - $\Phi/c$) $Fs_3$ (Ứng suất có hiệu - $u$ đo) $Fs_4$ (Đồ thị Wakita & Matsuo)
A (13/02/07) 3,20 101 1,72 2,49 2,66 -
B (16/03/07) 3,20 132 1,85 2,57 2,87 1,90
C (04/04/07) 4,80 151 1,69 1,92 2,31 1,70
D (25/05/07) 4,80 202 1,78 2,10 2,52 1,45
E (11/06/07) 7,00 219 1,49 1,48 1,79 1,38
F (17/08/07) 7,00 286 1,65 1,74 2,00 1,35
G (29/08/07) 7,80 298 1,57 1,56 1,84 1,32 
So sánh biến thiên Hệ số an toàn (Fs) qua các giai đoạn đắp:
Fs
     Điểm A      B       C       D   E   F   G

4. Phát hiện bất thường và diễn giải cơ chế

  • Chênh lệch phân tích ứng suất: Phân tích ổn định theo ứng suất có hiệu ($Fs_3$) luôn cho hệ số an toàn lớn nhất, vượt trội hơn phân tích theo ứng suất tổng ($Fs_1$) trung bình khoảng 35%. Nguyên nhân do giả định phân tán áp lực nước lỗ rỗng cục bộ trong phòng thí nghiệm cao hơn thực tế thoát nước qua bấc thấm.
  • Tính bảo thủ của phương pháp chuyển vị: Hệ số an toàn từ đồ thị Wakita & Matsuo ($Fs_4$) luôn cho giá trị nhỏ nhất ($1,32 \div 1,90$) và có xu hướng suy giảm liên tục khi độ lún đứng $\delta$ tích lũy tăng nhanh, ngay cả khi đất nền đang cố kết tốt. Điều này chứng minh đồ thị Wakita & Matsuo rất nhạy bén để đóng vai trò "ngưỡng chặn dưới" nhằm cảnh báo sớm rủi ro thi công.
  • Vị trí mặt trượt nguy hiểm: Quan trắc $Inclinometer$ (I1) xác nhận vị trí có chuyển vị ngang cực đại luôn tập trung ở cao độ $-6,0\text{m CDL}$ (khoảng $0,25$ chiều dày tầng sét yếu), hoàn toàn trùng khớp với đáy cung trượt nguy hiểm tìm được trong mô hình số FEM Plaxis.

Đóng góp khoa học và giá trị ứng dụng (Scientific & Practical Contributions)

                       CÁC GIÁ TRỊ ĐÓNG GÓP
[ĐÓNG GÓP HỌC THUẬT & LÝ THUYẾT]             [ỨNG DỤNG THỰC TIỄN & QUY PHẠM]

1. Đóng góp về mặt lý thuyết và học thuật

  • Chuẩn hóa bộ thông số cố kết - cường độ cho khu vực Hiệp Phước: Cung cấp bộ tham số thực nghiệm đáng tin cậy ($S = 0,24; m = 0,8; N^*_k = 10 \div 10,86$) cho loại đất sét bão hòa nước có tính dẻo cao tại miền Nam Việt Nam.
  • Làm sáng tỏ cơ chế gia tăng cường độ: Chứng minh sự tương đồng và khác biệt giữa lý thuyết SHANSEP quốc tế và phương pháp tính theo góc ma sát cố kết không thoát nước $\phi_{cu}$ theo quy chuẩn Việt Nam.

2. Đổi mới phương pháp luận và kiểm soát rủi ro

  • Thiết lập khung phân tích ổn định 4 chiều (4-Pillar Stability Framework), kết hợp giữa mô phỏng phần tử hữu hạn chuyển vị lớn và kiểm tra cân bằng giới hạn dựa trên dữ liệu Piezometer/Inclinometer trực tiếp từ hiện trường.

3. Ứng dụng thực tiễn trong thiết kế và thi công

  • Cho phép kỹ sư tự tin thiết kế biểu đồ đắp phân đợt (tốc độ đắp trung bình $0,5\text{ m/tuần}$) với chiều cao đắp mỗi đợt được kiểm soát chính xác dựa trên độ lún hồi tiếp qua phương pháp Asaoka.
  • Đưa ra khuyến nghị thiết kế an toàn: Nên sử dụng phương pháp giải tích Bishop đơn giản hóa kết hợp thông số ứng suất tổng $S_u$ dự báo để tính toán kiểm tra ổn định, vì phương pháp này mang tính thiên về an toàn và loại bỏ được sai số đo đạc áp lực nước lỗ rỗng.

Đối tượng quan tâm và lợi ích thực tiễn (Target Audience & Benefits)

  • Kỹ sư Địa kỹ thuật & Tư vấn thiết kế công trình:
    • Tiếp cận phương pháp dự báo $\Delta S_u$ chuẩn xác, không bị phụ thuộc hoàn toàn vào các thí nghiệm đắt đỏ trong phòng.
    • Lựa chọn hệ số an toàn $Fs$ hợp lý cho từng giai đoạn thi công nền đắp cảng biển, đường cao tốc.
  • Nhà thầu thi công hạ tầng & Quản lý dự án:
    • Nắm bắt công cụ kiểm soát an toàn thi công theo thời gian thực bằng đồ thị Wakita & Matsuo từ dữ liệu $Inclinometer$ và $Plate\text{ }Settlement$.
    • Tối ưu hóa thời gian chờ lún giữa các đợt đắp, rút ngắn tiến độ bàn giao mặt bằng mà vẫn triệt tiêu nguy cơ sạt lở.
  • Cơ quan quản lý nhà nước & Hội đồng thẩm định kỹ thuật:
    • Có cơ sở khoa học để đánh giá, thẩm tra các giải pháp xử lý nền đất yếu phức tạp; làm tài liệu tham khảo bổ sung cho việc soát xét các tiêu chuẩn chuyên ngành giao thông thủy và đường bộ.
  • Các nhà nghiên cứu và học viên sau đại học:
    • Nguồn dữ liệu thực nghiệm toàn diện về địa kỹ thuật công trình Cảng SPCT phục vụ nghiên cứu mô hình hóa phần tử hữu hạn nâng cao ($Plaxis\text{ }Soft\text{ }Soil\text{ }Creep, Cam-Clay$).

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

1. Phát hiện mang tính đột phá nhất của nghiên cứu này là gì?

Nghiên cứu chứng minh rằng phân tích ổn định theo ứng suất có hiệu dựa trên số liệu quan trắc áp lực nước lỗ rỗng thực tế cho hệ số an toàn $Fs$ lớn hơn tới 35% so với phân tích theo ứng suất tổng từ sức kháng cắt dự báo ($S_u$). Do đó, để an toàn trong thi công, kỹ sư bắt buộc phải sử dụng phương pháp ứng suất tổng làm tiêu chí nghiệm thu ổn định.

2. Phương pháp dự báo sức kháng cắt nào chuẩn xác nhất cho bùn sét Hiệp Phước?

Phương pháp dự báo theo 22 TCN 244-98 ($\Delta S_u = U \cdot \Delta\sigma_v \cdot \tan\phi_{cu}$) kết hợp độ lún cuối cùng theo Asaoka cho kết quả sát thực tế nhất với thí nghiệm cắt cánh kiểm chứng ($VST$). Công thức SHANSEP ($S=0,24; m=0,8$) có xu hướng dự báo thấp hơn khoảng 15%, thích hợp cho phương án thiết kế thận trọng.

3. Tại sao hệ số an toàn tra từ đồ thị Wakita & Matsuo lại có giá trị thấp nhất?

Đồ thị Wakita & Matsuo phản ánh trực tiếp tỷ số biến dạng ngang/đứng ($\gamma_s/\delta$). Khi đất lún sâu tích lũy theo thời gian, tỷ số này bị khuếch đại khiến hệ số an toàn trên biểu đồ suy giảm. Phương pháp này không tính trực tiếp đến sự gia tăng cường độ do cố kết, nên đóng vai trò là công cụ cảnh báo rủi ro biến dạng sớm cực kỳ an toàn trên công trường.

4. Kết quả nghiên cứu có thể mở rộng (generalize) cho các khu vực khác không?

Hoàn toàn có thể áp dụng quy trình tính toán cho toàn bộ khu vực Đồng bằng sông Cửu Long và ven biển Đông Nam Bộ có cấu trúc địa tầng tương đồng (lớp sét hữu cơ mặt bão hòa nước, $I_L > 1$). Tuy nhiên, cần xác định lại các tham số đặc trưng địa phương như $\phi_{cu}$, hệ số $N^*_k$ và chỉ số dẻo $PI$.

5. Những thiết bị quan trắc nào là tối thiểu cần phải có trên công trường xử lý bấc thấm gia tải?

Cần tối thiểu 3 nhóm thiết bị: Bản đo lún mặt ($SP$) để dự báo mức độ cố kết qua phương pháp Asaoka; Ống đo chuyển vị ngang ($Inclinometer$) tại chân ta-luy để kiểm soát cung trượt sâu; và Đầu đo áp lực nước lỗ rỗng ($Piezometer$) để theo dõi tốc độ tiêu tán áp lực nước trong suốt quá trình đắp cát.

Kết luận (Conclusion)

Báo cáo nghiên cứu khoa học của nhóm tác giả Bộ môn Địa cơ Nền móng – ĐH Bách Khoa TP.HCM đã giải quyết trọn vẹn bài toán cơ học phức tạp về kiểm soát ổn định trượt sâu công trình đắp gia tải trên nền đất yếu xử lý bấc thấm. Nghiên cứu khẳng định giá trị to lớn của việc tích hợp dữ liệu quan trắc hiện trường vào mô hình giải tích số, giúp chuyển hóa công tác quản lý chất lượng thi công từ trạng thái thụ động sang chủ động dự báo rủi ro chuẩn xác.

Trong tương lai, việc nghiên cứu mở rộng cần tích hợp thêm các cảm biến đo áp lực đất đáy khối đắp ($Earth\text{ }Pressure\text{ }Cell$), phân tích biến thiên chỉ số dẻo trong quá trình cố kết, và phát triển các thuật toán học máy (Machine Learning) để tự động hóa cảnh báo trượt lở thời gian thực trên các đại công trường hạ tầng giao thông biển.