Tổng quan về luận án

Nghiên cứu ứng xử cơ học và biến dạng cố kết của đất sét mềm yếu được xử lý bằng giải pháp gia tải trước kết hợp thoát nước thẳng đứng (Prefabricated Vertical Drains - PVD) là bài toán cốt lõi trong kỹ thuật địa kỹ thuật hiện đại. Tại Việt Nam, các dải trầm tích trẻ Đệ Tứ phân bố rộng khắp vùng đồng bằng Bắc Bộ (khoảng 15.000 $\text{km}^2$) và đồng bằng Nam Bộ với chiều dày tầng đất yếu biến thiên từ vài mét đến hơn 50 m. Đặc trưng cơ lý điển hình của nhóm đất này là trạng thái bão hòa nước cao ($S_r > 88 - 100%$), độ ẩm tự nhiên lớn ($W = 50 - 75%$, cá biệt than bùn lên tới $80 - 140%$), hệ số rỗng ban đầu cao ($e_0 = 1,2 - 2,07$), sức chống cắt không thoát nước rất thấp ($c_u < 35\text{ kPa}$, $N_{SPT} < 5$, $q_c < 0,1\text{ MPa}$) và hệ số thấm nhỏ ($k = 10^{-7} - 10^{-9}\text{ cm/s}$). Việc xây dựng công trình giao thông, thủy công và công nghiệp trên nền đất này thường đối mặt với nguy cơ lún lớn, lún lệch kéo dài và trượt sụt mất ổn định.

Khoảng trống nghiên cứu (Research Gap) then chốt xuất phát từ sự mất cân xứng giữa lý thuyết tính toán cố kết cổ điển với ứng xử thực tế tại hiện trường. Mặc dù công nghệ PVD kết hợp gia tải trước đã được ứng dụng phổ biến, hàng loạt sự cố nghiêm trọng vẫn xảy ra:

  • Đoạn Km82+500 – Km83+500 cao tốc Nội Bài – Lào Cai (2014) xuất hiện vết nứt vòng cung kéo dài gần 1 km do túi bùn xen kẹp và tốc độ lún nhanh vượt dự báo;
  • Đường dẫn cầu Văn Thánh (2003) lún tiếp diễn $2 - 3\text{ cm/tháng}$, vượt quá 50 cm so với tính toán thiết kế sau 2 năm khai thác;
  • Đoạn đắp cao $5,4\text{ m}$ Nam cầu Trìa Km732+100 Quốc lộ 1A (2001) trượt trồi trên cung trượt rộng 26 m do không kiểm soát tốc độ đắp theo diễn biến gia tăng sức kháng cắt;
  • Cầu Hoàng Long Km0+620 (1999) đắp đến chiều cao $6,8\text{ m}$ thì lún sụt $1,6 - 1,8\text{ m}$ dài 140 m do tốc độ tăng tải vượt quá tốc độ cố kết;
  • Đường dẫn cầu Đồng Niên trên Quốc lộ 5 (1995) lún thực tế $2,64 - 3,3\text{ m}$ sau 3 năm (gần gấp đôi mức dự báo $1,7\text{ m}$);
  • Kho cảng LPG Thị Vải (2002–2005) lún với tốc độ $1,1 - 1,2\text{ m/tháng}$ khi bắt đầu chịu tải.
+-----------------------------------------------------------------------------------+
|               KHỦNG HOẢNG TÍNH TOÁN & THI CÔNG NỀN ĐẤT YẾU TẠI VIỆT NAM            |
+-----------------------------------------------------------------------------------+
|  1. Đánh giá sai lệch Độ Cố Kết (ĐCK) giữa Áp Lực Nước Lỗ Rỗng và Đo Lún Hiện Trường |
|  2. Dự báo thiếu tin cậy Sức Chống Cắt Không Thoát Nước (Su) theo tiến trình tải   |
|  3. Lựa chọn sai mô hình cấu trúc đất (Constitutive Models) trong mô phỏng PTHH    |
+-----------------------------------------------------------------------------------+

Luận án của NCS. Nguyễn Hồng Trường (2017) với đề tài "Nghiên cứu đặc tính cơ học của đất yếu được xử lý bằng gia tải trước kết hợp thoát nước thẳng đứng" (Chuyên ngành Địa kỹ thuật Xây dựng, Mã số: 62-58-60-01; Người hướng dẫn: PGS.TS. Nguyễn Hữu Thái, TS. Nguyễn Tiếp Tân; Cơ sở đào tạo: Trường Đại học Thủy Lợi) đã giải quyết triệt để các câu hỏi và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Nguyên nhân dẫn đến sự sai lệch giữa phương pháp đánh giá độ cố kết theo áp lực nước lỗ rỗng (ALNLR) và theo số liệu quan trắc lún (Asaoka, Hyperbolic)? Đâu là vị trí quan trắc tối ưu phản ánh đúng độ cố kết trung bình của nền?
  • Câu hỏi nghiên cứu 2 (RQ2): Các hàm thực nghiệm quốc tế dự báo chỉ số nén ($C_c$) và chỉ số nở ($C_s$) có áp dụng chính xác cho đất sét yếu đồng bằng Bắc Bộ và Nam Bộ không, và mô hình nào có độ tin cậy cao nhất?
  • Câu hỏi nghiên cứu 3 (RQ3): Phạm vi ứng dụng của công thức (V.6) trong Quy trình 22TCN 262-2000 về dự tính sức chống cắt không thoát nước ($S_u$) là gì, và làm thế nào để xác định $S_u$ chính xác tại thời điểm cố kết bất kỳ?
  • Câu hỏi nghiên cứu 4 (RQ4): Trong các mô hình đàn - dẻo hiện đại (Mohr-Coulomb, Soft Soil, Cam-Clay cải tiến), mô hình nào phản ánh đúng nhất quan hệ ứng suất – biến dạng và cố kết thấm của sét yếu khi dùng phần tử hữu hạn (PTHH)?

Hệ thống giả thuyết tương ứng:

  • Giả thuyết 1 (H1): Độ cố kết xác định theo độ lún thực tế phản ánh toàn bộ biến dạng của khối đất, trong khi ALNLR chịu ảnh hưởng cục bộ bởi vùng xáo trộn ($d_s$) và khoảng cách đến bấc thấm; tồn tại điểm đặc trưng trong lưới bấc thấm mà tại đó ĐCK đo được tiệm cận ĐCK trung bình.
  • Giả thuyết 2 (H2): Hàm dự báo $C_c, C_s$ dựa trên độ ẩm tự nhiên ($W$) và hệ số rỗng ban đầu ($e_0$) bản địa hóa sẽ giảm sai số trung bình xuống dưới $15%$ so với sai số trên $35%$ của các công thức quốc tế truyền thống.
  • Giả thuyết 3 (H3): Sự gia tăng sức chống cắt không thoát nước $\Delta S_u(t)$ phụ thuộc phi tuyến vào độ cố kết hiện thời $U(t)$ và hệ số cố kết thực nghiệm $(m, \alpha)$, vượt qua giới hạn của công thức tuyến tính giản lược trong 22TCN 262-2000.
  • Giả thuyết 4 (H4): Mô hình Cam-Clay cải tiến (Modified Cam-Clay - MCC) với đường phá hoại dẻo biến cứng Ellipsoid là mô hình số thích hợp nhất để mô phỏng sự suy giảm thể tích lỗ rỗng và phát triển ứng suất hữu hiệu theo thời gian.

Phạm vi thực nghiệm bao gồm: 02 mô hình vật lý (MHVL) tỷ lệ lớn đối với đất sét yếu Yên Nghĩa - Hà Nội (MHVL 1 đạt cố kết $U=50%$, MHVL 2 đạt $U=90%$; đo lún và đo ALNLR bằng Piezometer dây rung tại các độ sâu $z = 25, 50, 75\text{ cm}$ và khoảng cách tim bấc $r = 15, 50\text{ cm}$ trong suốt 165 ngày); tổ hợp thí nghiệm nén ba trục cố kết không thoát nước (CU), nén cố kết Oedometer, cắt cánh trong phòng và hiện trường trên 13 dự án trọng điểm (Cao tốc Cầu Giẽ – Ninh Bình, Hà Nội – Hải Phòng, Nội Bài – Lào Cai, Tân Vũ – Lạch Huyện, Nhiệt điện Long Phú 1, Kho LPG Thị Vải, Đồng Tháp, Cần Thơ).

Literature Review và Positioning

Khung lý thuyết về cố kết thấm xung quanh hệ thống thoát nước thẳng đứng khởi nguồn từ lý thuyết cố kết một chiều của Karl Terzaghi (1925), giả định dòng thấm đơn hướng và hệ số thấm $k$ không đổi. Để mở rộng cho bài toán không gian đối xứng trục quanh giếng cát và bấc thấm:

  • Carrilo (1942) đề xuất nguyên lý tích phân phân tách độ cố kết tổng hợp 3 chiều ($U_{vh}$): $$U_{vh} = 1 - (1 - U_h)(1 - U_v)$$ trong đó $U_h$ và $U_v$ lần lượt là độ cố kết theo phương ngang và phương đứng.
  • Barron (1948) xây dựng nền tảng giải tích cố kết đối xứng trục với hai giả thiết biên: "biến dạng đều" (Equal Vertical Strain) và "biến dạng tự do" (Free Strain). Barron đưa ra hàm khoảng cách hình học: $$F(n) = \frac{n^2}{n^2 - 1}\ln(n) - \frac{3n^2 - 1}{4n^2} \approx \ln(n) - \frac{3}{4}$$ với $n = d_e / d_w$, trong đó $d_e$ là đường kính ảnh hưởng ($d_e = 1,05S$ cho lưới tam giác, $d_e = 1,13S$ cho lưới vuông), $d_w$ là đường kính tương đương của bấc thấm ($d_w = 2(a+b)/\pi$).
  • Hansbo (1981, 1997) hoàn thiện lời giải giải tích khi đưa vào xét đồng thời vùng xáo trộn (Smear Zone) do cần xuyên cắm bấc gây ra và sức cản dòng chảy dọc trục bấc thấm (Well Resistance): $$U_h = 1 - \exp\left(-\frac{8T_h}{F}\right)$$ $$F = F(n) + F_s + F_r$$ trong đó hệ số xáo trộn $F_s = \left(\frac{k_h}{k_s} - 1\right)\ln\left(\frac{d_s}{d_w}\right)$, và hệ số cản giếng $F_r = \pi z (2l - z) \frac{k_h}{q_w}$.
                 TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT CỐ KẾT THẤM
Terzaghi (1925)       Carrilo (1942)       Barron (1948)           Hansbo (1981, 1997)
[Cố kết 1D] -----> [Tổ hợp Uh & Uv] ---> [Đối xứng trục] ---> [Xét Smear Zone Fs & Cản Fr]

Các cuộc tranh luận học thuật cốt lõi

Tranh luận 1: Sự biến thiên của hệ số thấm ($k$) và hệ số cố kết ngang ($C_h$) trong quá trình ép nén. Trường phái cổ điển (Hansbo 1981; Jamiolkowski et al. 1991) giả định $k_h$ và $C_h$ là các hằng số không đổi trong toàn bộ quá trình gia tải. Ngược lại, các nghiên cứu thực nghiệm của Lê Bá Vinh và cs. (2010), Nguyễn Viết Tình (2010) chứng minh rằng khi hệ số rỗng $e$ giảm mạnh dưới áp lực gia tải, hệ số thấm $k$ giảm theo hàm logarit phi tuyến, khiến đường cong lún thực tế bị kéo dài đáng kể ở pha sau ($U > 80%$).

Tranh luận 2: Tính dị hướng thấm và tỷ số $C_h/C_v$. Trong Quy trình 22TCN 262-2000, tỷ số $C_h/C_v$ được cho phép dao động trong khoảng rất rộng ($2 - 5$). Các nghiên cứu chuyên sâu của Nguyễn Thị Nụ (2012) tại ven biển đồng bằng sông Cửu Long chỉ ra rằng với áp lực nén $< 0,5\text{ kPa}$, đất bùn sét hữu cơ có $C_h/C_v = 6 - 7$, nhưng khi nén từ $0,5 - 4\text{ kPa}$, tỷ số này giảm về $3 - 3,7$. Với bùn sét pha, tỷ số biến thiên từ $3,8 - 5$ xuống $2,6 - 3$. Việc chọn sai tỷ số này dẫn đến sai số dự báo thời gian dỡ tải hàng tháng trời.

Định vị so sánh với các công trình quốc tế

  1. Nghiên cứu đất sét Bangkok (Bergado et al., 1996; Indraratna et al., 2005): Phân tích ảnh hưởng của vùng xáo trộn $d_s = (2 - 3)d_m$ và $k_h/k_s = 1,5 - 5$ trên đất sét mềm Bangkok. Luận án của NCS. Nguyễn Hồng Trường kế thừa mô hình xáo trộn của Bergado nhưng định lượng hóa riêng cho các trầm tích sét vịnh biển - đầm lầy Hà Nội ($bmQ_2IV$) và Nam Bộ ($ambQIV$).
  2. Dự án lấn biển Sân bay Changi, Singapore (Choa et al., 1981; Bo et al., 2003): Ứng dụng PVD kết hợp gia tải chân không quy mô lớn, khẳng định phương pháp Asaoka cho độ tin cậy cao ở giai đoạn cố kết muộn ($U > 60%$) nhưng đánh giá thấp lún tức thời ban đầu.
  3. Công trình Cảng Busan và Sân bay Incheon, Hàn Quốc (Kim et al., 2005): So sánh hiệu quả giữa giếng cát, cọc cát đầm chặt và PVD trên nền sét biển dày $20 - 70\text{ m}$, chứng minh PVD vượt trội về kinh tế và độ ổn định khi được kiểm soát tốc độ đắp theo lộ trình cố kết.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã tạo ra bước tiến quan trọng trong cơ học đất và kỹ thuật xử lý nền khi giải mã và định lượng hóa cơ chế tương quan giữa quá trình tiêu tan áp lực nước lỗ rỗng, sự dịch chuyển trạng thái nén lún và sự tăng trưởng sức kháng cắt không thoát nước ($S_u$).

                      KHUNG MÔ HÌNH HÓA ỨNG XỬ ĐẤT SÉT YẾU
                      
+-----------------------------------------------------------------------------------+
|               TRẠNG THÁI BAN ĐẦU: $e_0, W, LL, PL, OCR \approx 1,0, S_{u0}$        |
+-----------------------------------------------------------------------------------+
                 [Cơ chế Cố kết Thấm Hansbo: $d_s, k_h/k_s, F(n), q_w$]
+-----------------------------------------------------------------------------------+
|              TIÊU TAN ÁP LỰC NƯỚC LỖ RỖNG DƯ: $\Delta u(r, z, t) \to 0$           |
|              TĂNG ỨNG SUẤT CÓ HIỆU: $\sigma'_v(t) = \sigma_{v0} + \Delta\sigma - u|
+-----------------------------------------------------------------------------------+
[Biến dạng Nén thể tích]                                         [Tăng Cường độ Chống cắt]
$S_c(t) = \frac{H_c}{1+e_0} C_c \log\left(\frac{\sigma'_0+\Delta\sigma'}{\sigma'_0}\right)$    $S_u(t) = S_{u0} + m \cdot U(t) \cdot \Delta\sigma'_v \cdot \tan\varphi_{cu}$
(Cam-Clay Yield Surface M)                                       (Xác lập qua Thí nghiệm Cắt cánh)

Luận án đã mở rộng và hoàn thiện lý thuyết trạng thái giới hạn (Critical State Soil Mechanics) của Roscoe & Burland (1968) và Casagrande (1936) khi áp dụng vào môi trường sét bùn trầm tích Holocen Việt Nam:

  1. Làm rõ và mở rộng công thức (V.6) trong Quy chuẩn 22TCN 262-2000: Quy chuẩn cũ áp dụng công thức tính sức chống cắt gia tăng: $$S_u = S_{u0} + \Delta\sigma_v \cdot \tan\varphi_{cu}$$ vốn chỉ đại diện cho trạng thái lý tưởng khi độ cố kết đạt $U = 100%$ và chỉ xét tải trọng đẳng hướng. Luận án đã phát triển công thức tổng quát xác định $S_u(t)$ tại thời điểm cố kết $U(t)$ bất kỳ: $$S_u(t) = S_{u0} + m \cdot U(t) \cdot \Delta\sigma_v \cdot \tan\varphi_{cu}$$ hoặc biểu diễn theo hàm thực nghiệm số mũ thích hợp: $$S_u(t) = S_{u0} + \alpha \cdot (\sigma'_v)^m$$ trong đó các tham số $m$ và $\alpha$ được chuẩn hóa từ ma trận thực nghiệm cắt cánh hiện trường và mô hình vật lý.
  2. Xác lập tương quan chuyển dịch bề mặt dẻo (Yield Surface): Chứng minh thực nghiệm rằng dưới tải trọng gia tải nhiều cấp, quỹ đạo ứng suất có hiệu dịch chuyển vượt qua áp lực tiền cố kết $p_c$, đưa đất từ trạng thái quá cố kết nhẹ ($OCR \approx 1,0 - 1,2$) vào miền cố kết bình thường ($OCR = 1,0$), làm thay đổi đột ngột mô đun đàn dẻo từ $C_s$ sang $C_c$.

Khung phân tích độc đáo

Khung phân tích tích hợp đa chiều của luận án kết hợp 3 nhánh lý thuyết chính:

  • Lý thuyết dòng thấm hướng tâm và đối xứng trục có xét xáo trộn (Barron-Hansbo);
  • Lý thuyết cơ học đất trạng thái tới hạn và biến dạng dẻo tăng bền (Modified Cam-Clay);
  • Lý thuyết hồi quy phi tuyến đa biến địa thống kê.

Đặc biệt, luận án thiết lập điều kiện biên tương đương (Equivalent Plane-Strain Permeability) để chuyển đổi bài toán 3D đối xứng trục thực tế của mạng bấc thấm sang mô hình biến dạng phẳng 2D trong phần mềm Plaxis, thông qua hiệu chỉnh hệ số thấm tương đương $k_{ve}$: $$k_{ve} = \frac{2 B^2}{3 R^2 \left[\ln(R/r_w) - 3/4\right]} k_h$$ cho phép bảo toàn chính xác tốc độ tiêu tan áp lực nước lỗ rỗng và tốc độ phát triển lún mà không làm suy giảm tính ổn định số của lưới phần tử hữu hạn.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ nghiêm ngặt triết học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism) trong địa kỹ thuật công trình. Thiết kế nghiên cứu đa tầng bậc tích hợp giữa mô hình vật lý quy mô lớn trong phòng, thí nghiệm cơ học đất chuyên sâu và mô phỏng số phần tử hữu hạn kiểm chứng bằng dữ liệu quan trắc hiện trường thực tế.

+-----------------------------------------------------------------------------------+
|                            QUY TRÌNH NGHIÊN CỨU TOÀN DIỆN                         |
+-----------------------------------------------------------------------------------+
[MÔ HÌNH VẬT LÝ (MHVL)]       [THÍ NGHIỆM ĐỊA KỸ THUẬT]         [MÔ PHỎNG SỐ PLAXIS 2D]
- Kích thước khối đất chuẩn    - Oedometer (e-log p, Cc, Cs)     - Mohr-Coulomb (MC)
- Cắm bấc thấm PVD quy chuẩn   - Nén 3 trục CU (c', phi', u)     - Soft Soil (SS)
- Piezometer dây rung & đo lún  - Cắt cánh (Field/Lab Vane Su)    - Modified Cam-Clay (MCC)
- 2 cấp tải: U=50% và U=90%    - Thấm biến đổi (kh, kv)          - Chuyển đổi thấm kve
                 [KIỂM NGHIỆM ĐỐI CHIẾU 13 CÔNG TRÌNH THỰC TẾ]
                 - Cao tốc Cầu Giẽ - Ninh Bình, Hà Nội - Hải Phòng
                 - Cao tốc Nội Bài - Lào Cai, Đường Tân Vũ - Lạch Huyện
                 - Nhiệt điện Long Phú 1, Kho cảng LPG Thị Vải...

Quy trình nghiên cứu thực nghiệm

Thiết kế và vận hành 02 Mô hình Vật lý (MHVL) chịu tải trọng gia tải tĩnh 2 cấp:

  • Chuẩn bị mẫu đất: Đất sét dẻo nhão lấy tại khu vực Yên Nghĩa (Hà Nội) được nhào trộn đồng nhất ở độ ẩm tự nhiên $W = 62,5%$, khối lượng thể tích $\gamma = 1,58\text{ T/m}^3$, hệ số rỗng ban đầu $e_0 = 1,71$, độ bão hòa $G = 94,8%$. Đất được đưa vào thùng mô hình kích thước lớn, đầm nén kiểm soát dung trọng để đảm bảo tính đồng nhất tuyệt đối.
  • Bố trí bấc thấm và thiết bị quan trắc: Lắp đặt bấc thấm tiết diện $100 \times 4\text{ mm}$ ở trung tâm mô hình. Hệ thống đo lường bao gồm:
    • Đầu đo áp lực nước lỗ rỗng dây rung (Vibrating Wire Piezometer) đặt tại các độ sâu $z = 25\text{ cm}, 50\text{ cm}, 75\text{ cm}$ và phân bố tại hai cự ly bán kính $r = 15\text{ cm}$ và $r = 50\text{ cm}$ tính từ tim bấc thấm;
    • Thiết bị đọc tự động Portable Vibrating Wire Readout Box kết nối phần mềm Logview;
    • Bản đo lún bề mặt và đo lún phân tầng;
    • Thiết bị cắt cánh cầm tay (Model T0174/B1) đo trực tiếp $S_u$ theo các chu kỳ cố kết.
  • Chế độ chất tải: Gia tải tĩnh phân bố đều theo 2 cấp tải trọng nhằm ngăn chặn phá hoại cắt tức thời. MHVL 1 vận hành đến khi độ cố kết đạt $U = 50%$; MHVL 2 vận hành kéo dài 165 ngày đêm để đạt độ cố kết $U \ge 90%$.

Data và kỹ thuật phân tích

Hệ thống thí nghiệm trong phòng xác lập ma trận thông số địa kỹ thuật:

  • Thí nghiệm nén cố kết Oedometer: Xác định đường cong nén $e - \log p'$, áp lực tiền cố kết $p_c$ theo phương pháp Casagrande, chỉ số nén $C_c$, chỉ số nở $C_s$, và hệ số cố kết $C_v, C_h$ theo phương pháp Logarit thời gian (Casagrande) và Căn bậc hai thời gian (Taylor).
  • Thí nghiệm nén ba trục cố kết không thoát nước (CU) có đo áp lực nước lỗ rỗng: Xác định quỹ đạo ứng suất có hiệu ($p' - q$), các thông số độ bền chống cắt hữu hiệu $c' = 14\text{ kPa}, \varphi' = 11^\circ$, và hệ số áp lực nước lỗ rỗng Skempton $A_f$.
  • Thí nghiệm thấm cột nước giảm dần: Xác định hệ số thấm thẳng đứng $k_v = 1,2 \times 10^{-8}\text{ cm/s}$ và thấm ngang $k_h = 1,8 \times 10^{-8}\text{ cm/s}$, suy ra tỷ số dị hướng $k_h/k_v = 1,5$.
  • Phân tích số phần tử hữu hạn: Sử dụng phần mềm Plaxis 2D với lưới phần tử tam giác 15 nút (15-node triangular elements), thực hiện phân tích ghép nối cố kết thấm - biến dạng (Coupled Consolidation Analysis) so sánh 3 mô hình cấu trúc đất: Mohr-Coulomb (MC), Soft Soil (SS) và Cam-Clay cải tiến (Modified Cam-Clay - MCC).

Phát hiện đột phá và implications

Những phát hiện then chốt

+-----------------------------------------------------------------------------------+
|               TỔNG HỢP 4 PHÁT HIỆN ĐỘT PHÁ TỪ NGHIÊN CỨU THỰC NGHIỆM              |
+-----------------------------------------------------------------------------------+
| 1. Điểm quan trắc đại diện ĐCK trung bình: r = 0,5S và z = 0,5H                   |
| 2. Sai số hàm dự báo Cc, Cs bản địa hóa giảm từ >35% xuống <10%                   |
| 3. Tăng trưởng Su(t) phi tuyến: Xác lập chuẩn xác bộ tham số (m, alpha)           |
| 4. Mô hình Cam-Clay cải tiến (MCC) vượt trội trong mô phỏng cố kết thấm PVD       |
+-----------------------------------------------------------------------------------+

Phát hiện 1: Giải mã bản chất sai khác giữa phương pháp ALNLR và phương pháp Dự báo lún (Asaoka/Hyperbolic) – Đề xuất vị trí quan trắc tối ưu

Khi phân tích song song dữ liệu quan trắc từ MHVL và các đại dự án giao thông (Tân Vũ – Lạch Huyện, Nội Bài – Lào Cai, Long Phú 1), tác giả phát hiện:

  • Ở giai đoạn đầu gia tải ($U < 60%$), phương pháp quan trắc ALNLR phản ánh tốc độ cố kết nhanh hơn thực tế do nước lỗ rỗng thoát cục bộ xung quanh bấc thấm, trong khi phương pháp Asaoka và Hyperbolic chưa hội tụ đủ chuỗi số liệu lún để cho kết quả chính xác.
  • Ở giai đoạn sau ($U > 80%$), phương pháp dự báo lún Asaoka đạt độ tin cậy rất cao (sai số $< 3%$), trong khi đầu đo ALNLR thường bị trơ hoặc chịu ảnh hưởng của hiện tượng tiêu tán áp lực cục bộ không đại diện cho toàn khối đất.
  • Giải pháp đột phá: Luận án chứng minh điểm quan trắc lún và ALNLR phản ánh chuẩn xác nhất Độ Cố Kết Trung Bình ($U_{avg}$) của toàn khối nền là điểm nằm ở khoảng cách $r = 0,5S$ (trung điểm giữa hai bấc thấm) và độ sâu $z = 0,5H$ (giữa chiều dày tầng sét yếu). Khuyến nghị này chấm dứt việc bố trí cảm biến tùy tiện sát bấc hoặc tại biên làm sai lệch nghiệm thu công trình.

Phát hiện 2: Thiết lập hệ hàm dự báo chỉ số nén $C_c$ và chỉ số nở $C_s$ chuẩn xác cho các vùng đồng bằng Việt Nam

Phân tích đối chiếu 12 công thức quốc tế kinh nghiệm (Skempton 1944, Terzaghi & Peck 1967, Nishida 1956, Azzouz et al. 1976...) trên bộ dữ liệu địa chất hàng trăm mẫu đất tại 13 dự án cho thấy sai số của các hàm quốc tế lên tới $35 - 55%$. Luận án đã lựa chọn và hiệu chỉnh các hàm dự báo tối ưu:

  • Đồng bằng Bắc Bộ: Hàm tương quan theo độ ẩm tự nhiên $W$ và hệ số rỗng $e_0$: $$C_c = 0,009(W - 10) \quad \text{hoặc} \quad C_c = 0,35(e_0 - 0,5)$$ cho sai số trung bình so với thí nghiệm trực tiếp giảm xuống chỉ còn $8,6 - 11,2%$.
  • Đồng bằng Nam Bộ: Tỷ số $C_c/C_s$ được lượng hóa chính xác trong dải $7,5 - 9,2$ (trung bình $C_c = 8,35 C_s$), khắc phục tình trạng ước tính tùy tiện trong tính toán lún đàn hồi dỡ tải.

Phát hiện 3: Định lượng hóa chính xác tốc độ gia tăng sức chống cắt không thoát nước $S_u(t)$

Kết quả thí nghiệm cắt cánh tại MHVL và kiểm nghiệm tại Dự án Tân Vũ – Lạch Huyện, Nhà máy Nhiệt điện Long Phú 1 chỉ ra:

  • Tại thời điểm cố kết $U = 50%$, công thức (V.6) trong 22TCN 262-2000 đánh giá quá cao cường độ đất thực tế (dự tính $S_u$ cao hơn thực tế $28 - 34%$), tiềm ẩn nguy cơ gây trượt sụt công trình khi đắp tăng tải sớm.
  • Áp dụng công thức đề xuất của luận án: $$S_u(t) = S_{u0} + m \cdot U(t) \cdot \Delta\sigma'v \cdot \tan\varphi{cu}$$ với $m = 0,75 - 0,85$ cho đất sét Hà Nội và $m = 0,68 - 0,78$ cho sét yếu Nam Bộ, kết quả dự báo khớp hoàn toàn với số liệu cắt cánh hiện trường (độ lệch $< 5%$).
        SO SÁNH DỰ BÁO SỨC KHÁNG CẮT Su TẠI CÁC THỜI ĐIỂM CỐ KẾT
+-----------------------+-------------------+--------------------+--------------------+
| Trạng thái Cố kết     | 22TCN 262-2000    | Công thức Luận án  | Cắt cánh Thực tế   |
+-----------------------+-------------------+--------------------+--------------------+
| Ban đầu (U = 0%)      | 12,5 kPa          | 12,5 kPa           | 12,5 kPa           |
| Giai đoạn 1 (U = 50%) | 26,8 kPa (+114%)  | 19,8 kPa (+58%)    | 19,2 kPa (+54%)    |
| Giai đoạn 2 (U = 94%) | 36,4 kPa (+191%)  | 35,1 kPa (+181%)   | 34,8 kPa (+178%)   |
+-----------------------+-------------------+--------------------+--------------------+

Phát hiện 4: Xác lập mô hình cấu trúc đất tối ưu trong tính toán số phần tử hữu hạn (PLAXIS)

So sánh 3 mô hình đàn dẻo trên cùng bài toán cố kết bấc thấm:

  1. Mô hình Mohr-Coulomb (MC): Đánh giá thấp độ lún cố kết cuối cùng từ $25 - 40%$, không mô phỏng được sự suy giảm thể tích lỗ rỗng phi tuyến;
  2. Mô hình Soft Soil (SS): Mô phỏng lún tương đối tốt nhưng đánh giá tốc độ tiêu tan áp lực nước lỗ rỗng chậm hơn thực tế;
  3. Mô hình Cam-Clay cải tiến (Modified Cam-Clay - MCC): Thể hiện sự trùng khớp vượt trội với đường cong quan trắc lún và ALNLR theo thời gian trên cả MHVL và các mặt cắt thực tế tại cao tốc Nội Bài – Lào Cai và đường Tân Vũ – Lạch Huyện.

Limitations và Future Research

Nghiên cứu duy trì tính khách quan khoa học khi thẳng thắn thừa nhận các giới hạn biên:

  1. Hiệu ứng kích thước mô hình vật lý (Scale Effect): Dù đã kiểm tra bằng PTHH để giảm thiểu ảnh hưởng của thành biên mô hình, ma sát thành hộp và giới hạn biên cứng vẫn có thể gây phân tán ứng suất nhỏ so với nền đất vô hạn tại hiện trường.
  2. Bỏ qua biến dạng từ biến thứ cấp ($C_\alpha$): Luận án tập trung chuyên sâu vào biến dạng cố kết sơ cấp ($S_c$) do tiêu tan áp lực nước lỗ rỗng, chưa tích hợp mô hình lưu biến dẻo nhớt (Viscoplasticity) để dự báo lún từ biến sau 30–50 năm khai thác.
  3. Hiện tượng suy giảm khả năng thoát nước của bấc thấm theo thời gian (Clogging & Bending): Trong mô hình số, khả năng thoát nước $q_w$ của bấc thấm được giả định không đổi, chưa xét đến hiện tượng gập gãy bấc do lún lớn ($S > 2\text{ m}$) và hạt sét mịn bọc nghẹt vải lọc sau nhiều năm.

Chương trình nghiên cứu 10 năm tiếp theo

  • Mở rộng 1: Phát triển mô hình cấu trúc đàn - nhớt - dẻo (Creep-SCLAY1S) tích hợp hiệu ứng phá hủy cấu trúc tự nhiên (Structure Degradation) và tính dị hướng cảm ứng của sét mềm ven biển.
  • Mở rộng 2: Nghiên cứu bài toán cố kết thấm 3 chiều phi tuyến xét đến hiện tượng nghẹt bấc thấm (Well Clogging) và độ uốn gãy phi tuyến của lõi nhựa dưới biến dạng lún cực hạn.
  • Mở rộng 3: Tích hợp công nghệ trí tuệ nhân tạo (Machine Learning/Physics-Informed Neural Networks - PINNs) để cập nhật thời gian thực các thông số $C_h, k_h, S_u$ từ dữ liệu quan trắc tự động IoT tại công trường.

Tác động và ảnh hưởng

+-----------------------------------------------------------------------------------+
|                        TÁC ĐỘNG ĐA CHIỀU CỦA CÔNG TRÌNH                           |
+-----------------------------------------------------------------------------------+
|  HỌC THUẬT: Định chuẩn phương pháp luận kiểm tra và mô phỏng cố kết sét yếu       |
|  KINH TẾ:   Tiết kiệm 15-25% chi phí xử lý nền; chống lãng phí bấc thấm & gia tải  |
|  KỸ THUẬT:  Xóa bỏ sự cố trượt lún do tăng tải sớm (như cầu Trìa, Văn Thánh)     |
|  CHÍNH SÁCH: Cung cấp cơ sở khoa học sửa đổi Quy trình 22TCN 262-2000             |
+-----------------------------------------------------------------------------------+
  • Tác động học thuật: Đóng góp hệ thống dữ liệu thực nghiệm chuẩn xác về đất sét yếu Việt Nam vào cơ sở dữ liệu địa kỹ thuật khu vực Đông Nam Á; cung cấp phương pháp luận mẫu mực về kết hợp mô hình vật lý và phần tử hữu hạn.
  • Tác động ngành và kinh tế: Giúp các tổng công ty tư vấn thiết kế (TEDI, VCE, Thủy Lợi) xác định chính xác thời điểm dỡ tải đắp gia tải, rút ngắn tiến độ thi công từ 2–4 tháng cho các dự án đường cao tốc, tiết kiệm hàng chục tỷ đồng chi phí bảo trì xử lý lún bù sau lún.
  • An toàn công trình: Ngăn ngừa nguy cơ phá hoại trượt sụt nền đắp do đắp tăng tải vượt quá tốc độ tăng sức kháng cắt của nền đất, đảm bảo an toàn tuyệt đối cho người và phương tiện trong giai đoạn thi công.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Nhà khoa học Địa kỹ thuật: Kế thừa bộ dữ liệu thực nghiệm nén ba trục, cắt cánh và mô hình số đã được kiểm chứng để phát triển các mô hình cơ học đất nâng cao.
  • Kỹ sư Thiết kế Địa kỹ thuật: Áp dụng trực tiếp hệ hàm dự báo $C_c, C_s$ và công thức tính $S_u(t)$ chính xác trong giai đoạn lập báo cáo nghiên cứu khả thi khi thiếu kinh phí thí nghiệm chuyên sâu.
  • Tư vấn Giám sát & Quản lý Dự án: Nắm vững vị trí đặt thiết bị quan trắc ($r=0,5S; z=0,5H$) và tiêu chuẩn phối hợp giữa biểu đồ Asaoka và đo áp lực nước lỗ rỗng để phê duyệt nghiệm thu dỡ tải an toàn.
  • Cơ quan Quản lý Nhà nước (Bộ Giao thông Vận tải, Bộ Nông nghiệp & PTNT, Bộ Xây dựng): Sử dụng các kết luận khoa học làm căn cứ cập nhật, sửa đổi các tiêu chuẩn, quy chuẩn kỹ thuật quốc gia về xử lý nền đất yếu.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Đóng góp lý thuyết độc đáo nhất là việc lượng hóa bản chất phi tuyến của sự gia tăng sức chống cắt không thoát nước $\Delta S_u(t)$ theo thời gian và mức độ cố kết $U(t)$, đồng thời hiệu chỉnh công thức (V.6) trong quy chuẩn 22TCN 262-2000 bằng cách đưa vào hệ số điều chỉnh thực nghiệm $m$ và quan hệ hàm số mũ $\alpha(\sigma'_v)^m$. Nghiên cứu đã chứng minh rằng cường độ chống cắt không tăng tuyến tính thuần túy với tải trọng đắp mà bị chi phối bởi mức độ tiêu tan áp lực nước lỗ rỗng thực tế trong cấu trúc đất dính.

2. Đột phá về phương pháp luận so với các nghiên cứu trước đây là gì?

So với các nghiên cứu trước đây (như Chang & Duncan, Hans-Georg Kempfert chỉ tập trung vào bài toán hố đào, hoặc các nghiên cứu cố kết đơn lẻ của Lê Bá Vinh, Nguyễn Thị Nụ):

  • Luận án đã thiết kế đồng bộ mô hình vật lý tỷ lệ lớn đa tầng quan trắc kết hợp giải tích và mô hình số PLAXIS với mô hình Cam-Clay cải tiến.
  • Xây dựng thành công thuật toán chuyển đổi hệ số thấm tương đương $k_{ve}$ bảo toàn chính xác động lực học tiêu tan áp lực nước lỗ rỗng khi chuyển từ bài toán 3D sang 2D Plane-Strain.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất?

Phát hiện bất ngờ nhất là sự suy giảm độ tin cậy nghiêm trọng của phương pháp đo áp lực nước lỗ rỗng (ALNLR) ở giai đoạn cố kết muộn ($U > 80%$), trái ngược với quan niệm truyền thống cho rằng đo ALNLR luôn chính xác hơn đo lún. Ở giai đoạn này, tốc độ biến thiên của ALNLR rất nhỏ, dễ bị nhiễu bởi độ nhạy cảm biến và hiệu ứng bấc thấm, trong khi phương pháp ngoại suy độ lún Asaoka lại đạt độ chính xác gần như tuyệt đối (độ lệch $< 3%$).

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) không?

Có. Luận án mô tả chi tiết quy trình chế bị mẫu đất bùn sét nhân tạo đồng nhất, kích thước thùng mô hình vật lý, quy cách cắm PVD, vị trí gắn cảm biến Piezometer dây rung bằng tọa độ $(r, z)$, các bước gia tải từng cấp, tốc độ tăng tải khống chế và phương pháp đo cắt cánh kiểm chứng $S_u$ theo từng mốc thời gian cố kết.

5. Khung nghiên cứu 10 năm tới được định hình như thế nào?

Khung nghiên cứu định hình việc tích hợp lý thuyết biến dạng từ biến phi tuyến (Non-linear Creep), mô hình hóa 3D động học nghẹt bấc thấm (Bio-chemical and Physical Clogging) và phát triển công nghệ quan trắc thông minh ứng dụng học máy để tự động tối ưu hóa lộ trình chất tải nền đắp trên đất yếu thời gian thực.

Kết luận

Công trình nghiên cứu của NCS. Nguyễn Hồng Trường đã đạt được 5 đóng góp đột phá có giá trị khoa học và thực tiễn bền vững:

  1. Xác lập quy chuẩn vị trí quan trắc hiện trường: Định vị chính xác tọa độ đại diện cho độ cố kết trung bình của khối nền đất yếu xử lý bằng bấc thấm tại trung điểm khoảng cách bấc ($r = 0,5S$) và trung điểm chiều sâu tầng đất sét ($z = 0,5H$).
  2. Bản địa hóa hệ hàm dự báo chỉ số nén $C_c$ và $C_s$: Xây dựng bộ công thức dự báo $C_c, C_s$ với độ chính xác cao (sai số $< 10%$) từ các chỉ tiêu vật lý đơn giản ($W, e_0$) cho hai vùng đồng bằng kinh tế trọng điểm Bắc Bộ và Nam Bộ.
  3. Hoàn thiện công thức dự tính sức chống cắt không thoát nước $S_u(t)$: Làm rõ giới hạn nguy hiểm của công thức (V.6) trong 22TCN 262-2000 và đề xuất công thức mới cho phép xác định chính xác $S_u$ tại bất kỳ thời điểm cố kết $U(t)$ nào, loại trừ triệt để nguy cơ trượt sụt do tăng tải quá nhanh.
  4. Chuẩn hóa mô hình số trong tính toán cố kết PTHH: Chứng minh mô hình Cam-Clay cải tiến (Modified Cam-Clay) là mô hình cấu trúc đất tối ưu nhất để phân tích bài toán cố kết thoát nước thẳng đứng, vượt trội hoàn toàn so với mô hình Mohr-Coulomb truyền thống.
  5. Mở ra hướng tiếp cận toàn diện cho cơ học đất ứng dụng: Kết hợp nhuần nhuyễn giữa thực nghiệm mô hình vật lý, lý thuyết giải tích đối xứng trục và mô phỏng số 2D/3D, tạo tiền đề vững chắc cho công tác thiết kế và thi công nền đắp hạ tầng giao thông trên nền đất yếu tại Việt Nam.