Tổng quan về luận án

Khu vực dải ven biển đồng bằng sông Hồng (ĐBSH) trải dài qua các tỉnh Hải Phòng, Thái Bình, Nam Định và Ninh Bình là trung tâm kinh tế năng động nhưng sở hữu điều kiện địa chất công trình phức tạp bậc nhất miền Bắc Việt Nam. Địa tầng bề mặt chủ yếu là các trầm tích Đệ tứ trẻ có nguồn gốc sông - biển hỗn hợp ($am$), biển ($m$), biển - đầm lầy ($mb$), đặc biệt là tầng sét yếu thuộc phụ thống Holocen thượng ($amQ_2^{3}tb$). Các lớp đất này có đặc tính cơ lý bất lợi: độ ẩm tự nhiên ($W$) cao vượt giới hạn chảy ($W_L$), hệ số rỗng ban đầu ($e_0$) lớn, sức kháng cắt không thoát nước ($c_u$) thấp và tính nén lún rất mạnh. Dưới tác động kép của tải trọng công trình đất đắp quy mô lớn (đường cao tốc, đê biển, bãi cảng, khu công nghiệp) và kịch bản biến đổi khí hậu với mực nước biển dâng dự báo tới $0{,}8\text{ m}$ trong 100 năm tại Hải Hậu - Nam Định (theo Bộ Tài nguyên và Môi trường), hiện tượng lún lớn và lún kéo dài trở thành thách thức kỹ thuật cốt lõi.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự thiếu hụt các nghiên cứu thực nghiệm định lượng về "hiệu ứng quy mô kích thước mẫu" (scale effect: đường kính $\Phi$ và chiều cao $H$) tác động trực tiếp đến các thông số cố kết cốt lõi ($C_c, c_v, p_c$) của đất sét yếu ven biển ĐBSH. Trong thực tế thiết kế, các kỹ sư thường áp dụng trực tiếp kết quả từ mẫu nén Oedometer tiêu chuẩn ($\Phi = 61{,}8\text{ mm}, H = 20\text{ mm}$) hoặc giả định tỷ số hệ số cố kết ngang trên cố kết đứng $c_h/c_v = 2 \div 5$ theo khuyến cáo thô của tiêu chuẩn TCCS 41:2022 mà thiếu cơ sở thực nghiệm chuyên sâu cho từng phân tầng thạch học khu vực, dẫn đến sai số dự báo lún và lãng phí trong thiết kế xử lý bấc thấm (PVD).

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Thị Bích Hạnh (Viện Khoa học Công nghệ Xây dựng - IBST, 2024, Chuyên ngành Kỹ thuật Địa chất, Mã số 9 52 05 01, dưới sự hướng dẫn khoa học của PGS. Đoàn Thế Tường và PGS. Đỗ Minh Toàn) được xây dựng nhằm giải quyết triệt để vấn đề trên thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Các đặc tính cơ lý, thành phần khoáng vật và cấu trúc nền (CTN) đất loại sét yếu ven biển ĐBSH có quy luật phân bố không gian và phân hóa địa chất công trình như thế nào?
  • Câu hỏi nghiên cứu 2 (RQ2): Sự thay đổi kích thước hình học của mẫu thí nghiệm (đường kính $\Phi = 61{,}8\text{ mm}, 79{,}8\text{ mm}$; chiều cao $H = 20\text{ mm}, 30\text{ mm}, 40\text{ mm}$) làm biến đổi các đặc trưng nén cố kết ($C_c, c_v, p_c$) theo phương thẳng đứng và phương ngang ($c_h$) ra sao?
  • Câu hỏi nghiên cứu 3 (RQ3): Việc tích hợp các thông số cố kết hiệu chỉnh theo kích thước mẫu và tỷ số dị hướng $c_h/c_v$ vào bài toán tính lún cho các công trình đất đắp trên các dạng cấu trúc nền điển hình mang lại độ tin cậy và hiệu quả kinh tế - kỹ thuật định lượng như thế nào?
  • Giả thuyết khoa học 1 (H1): Tồn tại 3 kiểu cấu trúc nền đất yếu đặc trưng ven biển ĐBSH, trong đó kiểu CTN II chiếm ưu thế chi phối toàn bộ bức tranh biến dạng nền đắp.
  • Giả thuyết khoa học 2 (H2): Chiều cao mẫu $H$ có mối tương quan nghịch phi tuyến với hệ số cố kết $c_v$ ở mọi cấp áp lực nén (đặc biệt trong dải áp lực hữu hiệu $\sigma' = 100 \div 200\text{ kPa}$), trong khi tăng đường kính mẫu $\Phi$ làm tăng nhẹ $c_v$ và giảm nhẹ chỉ số nén $C_c$.
  • Giả thuyết khoa học 3 (H3): Xác định chính xác tỷ số $c_h/c_v$ bằng thí nghiệm cố kết hướng tâm (Rowe cell) giúp tối ưu hóa chiều sâu, bước cắm bấc thấm PVD và thời gian chờ lún cố kết so với phương pháp tính toán truyền thống.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa lý thuyết cố kết một chiều Terzaghi (1925, 1943), lý thuyết cố kết ba chiều Biot (1941), lý thuyết cố kết hướng tâm xung quanh giếng thoát nước thẳng đứng của Barron (1948) và Hansbo (1981). Đóng góp mang tính đột phá của nghiên cứu là xác lập hệ tương quan thực nghiệm định lượng giữa kích thước mẫu thí nghiệm với các chỉ tiêu cố kết, đồng thời mô hình hóa ứng xử lún của các công trình đất đắp thực tế trên nền đất sét yếu $amQ_2^{3}tb$ tại các trọng điểm Hải Phòng, Thái Bình, Nam Định và Ninh Bình.

Literature Review và Positioning

Lịch sử nghiên cứu cơ học đất và kỹ thuật nền móng ghi nhận nền tảng lý thuyết cố kết thấm kinh điển do Terzaghi (1925, 1943) đề xuất dựa trên giả định dòng thấm một chiều, môi trường bão hòa đồng nhất và hệ số nén thể tích $m_v$, hệ số thấm $k$ là hằng số. Để khắc phục các giới hạn của mô hình một chiều, Biot (1941) và Rendulic (1936) đã tổng quát hóa lý thuyết cố kết sang không gian ba chiều xét đến tính nén được của cốt đất và chất lưu lỗ rỗng. Trong lĩnh vực gia cố nền đất yếu bằng hệ thống tiêu thoát nước thẳng đứng, Barron (1948) đã xây dựng giải pháp giải tích toàn diện cho bài toán cố kết hướng tâm theo hai giả thiết biến dạng tự do (free strain) và biến dạng đều (equal strain). Tiếp đó, Hansbo (1981) hoàn thiện mô hình cố kết hướng tâm có xét đến sức cản thủy lực của bấc thấm ($F_r$) và vùng xáo động xung quanh thiết bị cắm ($F_s$).

Trên bình diện quốc tế, các tranh luận học thuật về hiệu ứng kích thước mẫu (sample size/scale effect) trong thí nghiệm Oedometer diễn ra sôi nổi với nhiều quan điểm đối lập:

  • Trường phái 1 (Nghịch biến $c_v - H$): Wang Yao et al. (2023) khi nghiên cứu đất sét hoàng thổ tại Tây An (Trung Quốc) với các tỷ lệ chiều cao/đường kính $H/\Phi = 0{,}32 \div 1{,}62$ ($H = 20 \div 100\text{ mm}$) khẳng định chỉ số nén $C_c$ giảm theo cấp số nhân khi $H/\Phi$ tăng, trong khi áp lực tiền cố kết $p_c$ tăng tuyến tính. Tương tự, Ronghua Hu, Ming Zhang & Jiaqi Wang (2023) thí nghiệm trên bùn nạo vét tại Thâm Quyến (Trung Quốc) với $H = 20, 30, 40\text{ mm}$ ($\Phi = 61{,}8\text{ mm}$) đã chỉ ra rằng $C_c$ giảm và $c_v$ giảm rõ rệt khi chiều cao mẫu tăng, đồng thời cảnh báo việc dùng mẫu tiêu chuẩn $H = 20\text{ mm}$ dẫn đến dự báo sai lệch lớn về tốc độ cố kết hiện trường.
  • Trường phái 2 (Đồng biến hoặc Không đổi): Al-Khafaji (2016) tại Đại học Salford (Anh) tiến hành thí nghiệm trên sét Kaolin bão hòa với đường kính $\Phi = 100, 150, 250\text{ mm}$ và chiều cao $H = 23 \div 200\text{ mm}$, lại đưa ra kết luận ngược lại rằng khi chiều cao mẫu tăng thì giá trị $c_v$ đo được có xu hướng tăng ở mỗi cấp gia tải. Trong khi đó, Kaothon et al. (2014) khảo sát đất sét mềm Băng Cốc (Thái Lan) với $\Phi = 35\text{ mm}$ và $60\text{ mm}$ ($H = 20\text{ mm}$) cho rằng các thông số $C_c, c_v, p_c$ biến đổi không đáng kể và quy sai khác nhỏ cho ma sát thành và xáo động cơ học.

Tại Việt Nam, các nghiên cứu địa kỹ thuật đất yếu đã được khởi xướng bởi Lê Trọng Thắng (1995) với công trình phân chia cấu trúc nền đất yếu khu vực Hà Nội; Nguyễn Viết Tình (2001) nghiên cứu trầm tích hồ - đầm lầy phụ hệ tầng Hải Hưng dưới ($Q_2^{1-2}hh_1$), chỉ ra tỷ số $c_h/c_v = 1{,}06 \div 4{,}01$ và sự tồn tại của gradient cố kết ban đầu; Bùi Đức Hải (2003) làm rõ tính từ biến của đất sét yếu; Phạm Thị Nghĩa & Phạm Văn Tỵ (2005) sử dụng thiết bị xuyên đo áp lực nước lỗ rỗng (CPTu) xác định $c_h/c_v = 1{,}3 \div 7{,}6$; Nguyễn Thị Nụ et al. (2009, 2014, 2020) nghiên cứu đất yếu ven biển đồng bằng sông Cửu Long và Bắc Bộ bằng hộp nén Rowe và CRS-R, xác định khoảng biến thiên $c_h/c_v = 1{,}35 \div 10{,}59$; Bùi Trường Sơn et al. (2020) phân tích đặc tính cố kết đất yếu Bắc Trung Bộ; Phú Truyền Nhật & Lê Vĩnh Ba (2020) và Võ Nhật Luân et al. (2021) đánh giá tính nén lún đất yếu tại TP. Hồ Chí Minh.

Vị thế học thuật của luận án: Luận án định vị chính xác vào khoảng trống khoa học chưa được giải quyết đầy đủ: thiết lập mối quan hệ định lượng giữa kích thước mẫu đa dạng ($\Phi = 61{,}8\text{ mm}, 79{,}8\text{ mm}$; $H = 20\text{ mm}, 30\text{ mm}, 40\text{ mm}$) với các thông số cố kết ($C_c, c_v, p_c$) chuyên biệt cho đất loại sét yếu Holocen $amQ_2^{3}tb$ ven biển ĐBSH. Công trình cung cấp minh chứng thực nghiệm ủng hộ luận điểm của Ronghua Hu et al. (2023) và Wang Yao et al. (2023), giải thích thấu đáo nguyên nhân ma sát thành và chiều dài đường thấm, đồng thời kết nối trực tiếp kết quả trong phòng với bài toán xử lý nền bấc thấm PVD cho công trình đắp thực tế.

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

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

Luận án kế thừa và mở rộng lý thuyết cố kết một chiều Terzaghi (1925) và lý thuyết thoát nước hướng tâm Barron (1948) - Hansbo (1981) qua các luận điểm khoa học vững chắc:

  1. Làm sáng tỏ bản chất cơ học của hiệu ứng quy mô mẫu: Chứng minh rằng trong thí nghiệm Oedometer, khi chiều cao mẫu $H$ tăng từ $20\text{ mm}$ lên $40\text{ mm}$, chiều dài đường thấm thoát nước tăng gấp đôi, làm gia tăng lực cản thấm nội tại và độ trễ tiêu tán áp lực nước lỗ rỗng $u$, dẫn đến hệ số cố kết thẳng đứng $c_v$ suy giảm ở mọi cấp áp lực nén.
  2. Quy luật biến thiên của áp lực tiền cố kết ($p_c$) và chỉ số nén ($C_c$): Khi tăng đường kính mẫu $\Phi$ từ $61{,}8\text{ mm}$ lên $79{,}8\text{ mm}$, diện tích tiếp xúc biên tăng nhưng tỷ lệ diện tích biên/thể tích mẫu giảm, làm giảm ảnh hưởng ma sát thành dao vòng, khiến giá trị $C_c$ giảm nhẹ và $p_c$ tăng nhẹ.
  3. Mô hình hóa mối tương quan thực nghiệm: Xây dựng phương trình hồi quy tương quan thực nghiệm giữa $c_v$ và chiều cao $H$ trong dải áp lực hữu hiệu $\sigma' = 100 \div 200\text{ kPa}$ ($1 \div 2\text{ kG/cm}^2$), tạo cơ sở lý thuyết để quy đổi chỉ tiêu phòng thí nghiệm về trạng thái nén của phân tầng thực tế ngoài hiện trường.

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

Khung phân tích của luận án là sự tích hợp đa ngành giữa: Địa tầng học trầm tích Đệ tứ $\rightarrow$ Cơ học đất thực nghiệm nâng cao $\rightarrow$ Mô hình hóa toán học địa kỹ thuật công trình đắp.

+-------------------------------------------------------------------------------+
|                      KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN                     |
+-------------------------------------------------------------------------------+
|  1. ĐẶC ĐIỂM TRẦM TÍCH ĐỆ TỨ & CẤU TRÚC NỀN (CTN I, CTN II, CTN III)          |
|     - Trầm tích Holocen amQ2^3tb, mbQ2^1-2hh1, mQ2^1-2hh3...                  |
|     - Phân tích thành phần hạt, khoáng vật sét (Illite, Kaolinite, Smectite) |
+---------------------------------------+---------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
|  2. THỰC NGHIỆM ĐẶC TRƯNG CỐ KẾT THEO QUY MÔ KÍCH THƯỚC MẪU                  |
|     - Biến số hình học: Đường kính Phi (61.8, 79.8 mm) x Chiều cao H (20, 30, 40 mm)|
|     - Sơ đồ thoát nước: Cố kết đứng (Oedometer) & Cố kết ngang (Rowe cell)   |
|     - Trích xuất thông số: Cc, Cr, pc, OCR, cv, ch, tỷ số ch/cv              |
+---------------------------------------+---------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
|  3. ỨNG DỤNG TÍNH TOÁN XỬ LÝ NỀN ĐẤT ĐẮP THỰC TẾ                              |
|     - Dự báo lún tổng St (lún tức thời, lún cố kết Sc, lún từ biến)           |
|     - Thiết kế hệ thống bấc thấm PVD: Khoảng cách D_PVD, thời gian U(t)       |
|     - Kiểm toán ổn định trượt cung tròn (Bishop/Fellenius) đa giai đoạn đắp   |
+-------------------------------------------------------------------------------+

Điều kiện biên và phạm vi áp dụng: Áp dụng cho đất loại sét yếu bão hòa nước thuộc phụ thống Holocen thượng ($amQ_2^{3}tb$) dải ven biển ĐBSH (Hải Phòng, Thái Bình, Nam Định, Ninh Bình); dải tải trọng nén $\sigma' \le 400\text{ kPa}$; sơ đồ thoát nước hai chiều và hướng tâm.

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

Thiết kế nghiên cứu

Luận án áp dụng triết lý thực chứng (Positivism) kết hợp phương pháp nghiên cứu đa cấp độ (Multi-level Experimental Design):

  • Cấp độ 1 (Khảo sát hiện trường): Thu thập dữ liệu địa chất công trình từ các lỗ khoan sâu ($15 \div 75\text{ m}$) tại Diêm Điền (Thái Bình), Hải Thịnh (Nam Định), Kim Sơn (Ninh Bình); khoan lấy mẫu nguyên trạng bằng ống mẫu thành mỏng chuyên dụng (Shelby tube) để bảo toàn tối đa kết cấu tự nhiên.
  • Cấp độ 2 (Thực nghiệm phòng thí nghiệm có kiểm soát biến số): Thiết kế ma trận thí nghiệm đối chứng thay đổi kích thước mẫu hình học đa thông số trên thiết bị nén một trục không nở hông Oedometer và hộp nén hướng tâm cải tiến kiểu Rowe theo tiêu chuẩn TCVN 4200:2012 và BS 1377: Part 6: 1990.
  • Cấp độ 3 (Mô phỏng và kiểm chứng tính toán): Thiết lập mô hình giải tích dự báo lún và ổn định mái dốc đắp cho 3 dạng cấu trúc nền tiêu biểu có xử lý bấc thấm PVD.
+-------------------------------------------------------------------------------+
|                MA TRẬN THỰC NGHIỆM BIẾN DẠNG CỐ KẾT TRONG PHÒNG               |
+-------------------+-----------------------------+-----------------------------+
| Kích thước mẫu    | Đường kính Phi = 61,8 mm    | Đường kính Phi = 79,8 mm    |
+-------------------+-----------------------------+-----------------------------+
| Chiều cao H = 20mm| Nhóm mẫu A1 (H20 - Phi 61.8)| Nhóm mẫu B1 (H20 - Phi 79.8)|
| Chiều cao H = 30mm| Nhóm mẫu A2 (H30 - Phi 61.8)| Nhóm mẫu B2 (H30 - Phi 79.8)|
| Chiều cao H = 40mm| Nhóm mẫu A3 (H40 - Phi 61.8)| Nhóm mẫu B3 (H40 - Phi 79.8)|
+-------------------+-----------------------------+-----------------------------+
| Cố kết hướng tâm  | Hộp nén Rowe Phi = 62 mm, H = 30 mm (Lõi cát tiêu chuẩn)  |
+-------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

  • Chuẩn bị và chế tạo dụng cụ bổ trợ: Gia công chế tạo bộ dao vòng thép không gỉ đồng bộ có kích thước chính xác tuyệt đối: $\Phi = 61{,}8\text{ mm}$ và $\Phi = 79{,}8\text{ mm}$ với các chiều cao cắt mẫu tương ứng $H = 20\text{ mm}, 30\text{ mm}, 40\text{ mm}$. Gia công đá thấm, giấy thấm và piston nén đồng bộ.
  • Quy trình gia tải: Cấp áp lực gia tải thiết kế chuẩn hóa: $\sigma' = 12{,}5; 25; 50; 100; 200; 400\text{ kPa}$ (tương ứng $0{,}125 \div 4{,}0\text{ kG/cm}^2$). Bước dỡ tải được tiến hành sau khi kết thúc cấp nén cuối cùng để xác định chỉ số nở $C_r$.
  • Giao thức đo đạc: Ghi nhận chuyển vị lún tự động và thủ công tại các mốc thời gian chuẩn: $t = 6'', 15'', 30'', 45'', 1', 2', 4', 8', 15', 30', 1\text{h}, 2\text{h}, 3\text{h}, 5\text{h}, 8\text{h}, 24\text{h}$ cho tới khi biến dạng đạt trạng thái ổn định (tốc độ biến dạng lún không vượt quá $0{,}01\text{ mm}/24\text{h}$).
  • Xác định thông số: Áp dụng song song hai phương pháp đồ giải Casagrande ($\log t$) xác định $t_{50}, c_v = \frac{0{,}197 S_{50}^2}{t_{50}}$ và Taylor ($\sqrt{t}$) xác định $t_{90}, c_v = \frac{0{,}848 H_{90}^2}{t_{90}}$; xác định áp lực tiền cố kết $p_c$ theo phương pháp dựng đường phân giác Casagrande trên đường cong $e - \log \sigma'$.
  • Thí nghiệm cố kết hướng tâm: Sử dụng hộp nén Rowe ($\Phi = 62\text{ mm}, H = 30\text{ mm}$), tạo giếng cát trung tâm đường kính $d_w = 10\text{ mm}$ (cát thạch anh tiêu chuẩn $0{,}09 \div 0{,}3\text{ mm}$, hệ số thấm $k \approx 2\times 10^{-2}\text{ cm/s}$), đo đồng thời biến dạng đứng và áp lực nước lỗ rỗng $u_b$ tại biên ngoài để trích xuất $c_h = \frac{T_{ri} \cdot \Phi^2}{t_{50}}$.

Data và phân tích

  • Đặc trưng địa chất khoáng vật: Phân tích nhiễu xạ tia X (XRD) và soi hiển vi điện tử xác định thành phần khoáng vật của sét $amQ_2^{3}tb$: Illite chiếm chủ đạo ($40 \div 55%$), Kaolinite ($20 \div 30%$), Montmorillonite/Smectite ($10 \div 20%$) và khoáng phi sét thạch anh (Quartz).
  • Chỉ tiêu cơ lý tự nhiên: Đất ở trạng thái dẻo chảy đến chảy, độ ẩm tự nhiên $W = 48{,}5 \div 65{,}2%$, giới hạn chảy $W_L = 42 \div 58%$, giới hạn dẻo $W_P = 22 \div 29%$, chỉ số dẻo $I_P = 20 \div 31$, hệ số rỗng tự nhiên $e_0 = 1{,}32 \div 1{,}85$, dung trọng tự nhiên $\gamma = 15{,}8 \div 17{,}2\text{ kN/m}^3$.
  • Xử lý số liệu và kiểm tra độ tin cậy: Toàn bộ số liệu biến dạng - thời gian được số hóa, lọc nhiễu, kiểm tra tính vững (robustness checks) qua việc loại bỏ các mẫu bị phá hoại kết cấu cục bộ, phân tích hồi quy tương quan phi tuyến bằng phần mềm chuyên dụng.

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

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

  1. Quy luật nghịch biến giữa hệ số cố kết thẳng đứng $c_v$ và chiều cao mẫu $H$: Kết quả thực nghiệm khẳng định ở tất cả các cấp áp lực nén ($\sigma' = 12{,}5 \div 400\text{ kPa}$), khi chiều cao mẫu tăng từ $H = 20\text{ mm} \rightarrow 30\text{ mm} \rightarrow 40\text{ mm}$ (trên cả hai đường kính $\Phi = 61{,}8\text{ mm}$ và $\Phi = 79{,}8\text{ mm}$), hệ số cố kết $c_v$ đều suy giảm rõ rệt. Trong dải áp lực hữu hiệu quan trọng $100 \div 200\text{ kPa}$ ($1 \div 2\text{ kG/cm}^2$), luận án đã thiết lập được hàm tương quan thực nghiệm chuẩn xác giữa $c_v$ và $H$.
  2. Ảnh hưởng của đường kính mẫu $\Phi$ đến các chỉ tiêu nén: Khi tăng đường kính mẫu từ $\Phi = 61{,}8\text{ mm}$ lên $\Phi = 79{,}8\text{ mm}$ (ở cùng chiều cao $H$):
    • Chỉ số nén $C_c$ có xu hướng giảm nhẹ (do giảm ma sát thành tương đối).
    • Áp lực tiền cố kết $p_c$ tăng nhẹ.
    • Hệ số cố kết $c_v$ tăng nhẹ ở các cấp áp lực.
  3. Xác lập tỷ số dị hướng thấm cố kết $c_h/c_v$ thực nghiệm: Thí nghiệm cố kết hướng tâm trên hộp nén Rowe xác định hệ số cố kết ngang $c_h$ của đất sét yếu $amQ_2^{3}tb$ ven biển ĐBSH luôn lớn hơn hệ số cố kết đứng $c_v$. Tỷ số $c_h/c_v$ thực nghiệm dao động trong khoảng $1{,}45 \div 2{,}85$ tùy thuộc cấp áp lực nén, phản ánh chính xác cấu trúc phân lớp vi mô của trầm tích sông - biển Đệ tứ.
  4. Phân chia 3 kiểu cấu trúc nền đất yếu ven biển ĐBSH:
    • Kiểu CTN I (Đất yếu mỏng phân bố nông): Chiều dày lớp sét yếu $< 5\text{ m}$, bên dưới là cát/sét chặt. Xử lý thuận lợi bằng đệm cát hoặc bấc thấm nông.
    • Kiểu CTN II (Đất yếu dày trung bình đến lớn, dạng đơn tầng hoặc đa tầng kẹp cát): Chiều dày sét yếu $5 \div 25\text{ m}$, là kiểu cấu trúc phổ biến nhất tại ven biển Thái Bình, Nam Định, Ninh Bình, gây lún cố kết lớn và nguy cơ mất ổn định trượt cao nhất dưới tải trọng đắp.
    • Kiểu CTN III (Đất yếu dạng thấu kính xen kẹp phức tạp): Phân bố tại các vùng cửa sông, bãi triều lầy ven biển Hải Phòng - Ninh Bình.
+----------------------------------------------------------------------------------------------------+
|               TỔNG HỢP SO SÁNH CÁC PHÁT HIỆN THỰC NGHIỆM VỚI NGHIÊN CỨU QUỐC TẾ                    |
+------------------------------+---------------------------+-----------------------------------------+
| Nghiên cứu & Địa bàn         | Đối tượng đất             | Kết luận về hiệu ứng kích thước mẫu     |
+------------------------------+---------------------------+-----------------------------------------+
| Luận án (2024, ĐBSH)         | Sét yếu Holocen amQ2^3tb  | H tăng -> cv giảm rõ rệt; Phi tăng ->   |
|                              | Phi = 61.8, 79.8; H=20..40| Cc giảm nhẹ, pc & cv tăng nhẹ.          |
+------------------------------+---------------------------+-----------------------------------------+
| Ronghua Hu et al. (2023, TQ) | Bùn nạo vét Thâm Quyến    | H tăng (20->40mm) -> cv giảm, Cc giảm.  |
|                              | Phi = 61.8; H = 20, 30, 40| Cảnh báo sai số khi dùng H = 20mm.      |
+------------------------------+---------------------------+-----------------------------------------+
| Wang Yao et al. (2023, TQ)   | Sét hoàng thổ Tây An      | H/Phi tăng -> Cc giảm hàm mũ, pc tăng.  |
|                              | H/Phi = 0.32 -> 1.62      |                                         |
+------------------------------+---------------------------+-----------------------------------------+
| Al-Khafaji (2016, Anh)       | Sét Kaolin bão hòa        | H tăng (23->200mm) -> cv tăng.          |
|                              | Phi = 100, 150, 250 mm    | (Trái ngược với đất ĐBSH do cấu trúc đất)|
+------------------------------+---------------------------+-----------------------------------------+
| Kaothon et al. (2014, TL)    | Sét mềm Băng Cốc          | Thay đổi Phi = 35 vs 60mm -> Cc, cv, pc |
|                              | Phi = 35, 60 mm; H = 20 mm| thay đổi không đáng kể.                 |
+------------------------------+---------------------------+-----------------------------------------+

Implications đa chiều

  • Ý nghĩa khoa học: Cung cấp bằng chứng thực nghiệm khẳng định tác động của hiệu ứng hình học mẫu nén đến quá trình tiêu tán áp lực nước lỗ rỗng trong môi trường sét mềm bão hòa, bổ sung dữ liệu cơ bản cho cơ học đất trầm tích nhiệt đới gió mùa.
  • Đổi mới phương pháp luận: Khuyến nghị các phòng thí nghiệm LAS-XD khi thí nghiệm đất sét mềm ven biển cần chuẩn hóa việc lựa chọn chiều cao mẫu $H = 30\text{ mm}$ hoặc $40\text{ mm}$ kết hợp hiệu chỉnh đường cong cố kết nhằm phản ánh sát thực tế hiện trường thay vì rập khuôn mẫu $H = 20\text{ mm}$.
  • Ứng dụng thực tiễn công trình đất đắp: Khi áp dụng tính toán xử lý bấc thấm PVD cho nền đắp cao tốc và đê biển trên kiểu CTN II:
    • Việc sử dụng giá trị $c_v$ hiệu chỉnh theo kích thước mẫu thực tế giúp tránh hiện tượng đánh giá quá cao tốc độ cố kết (overestimation), loại bỏ nguy cơ lún dư vượt ngưỡng sau khi thông xe.
    • Việc áp dụng tỷ số thực nghiệm $c_h/c_v = 1{,}45 \div 2{,}85$ (thay vì giả định $c_h/c_v = 4 \div 5$) giúp thiết kế khoảng cách cắm bấc thấm $D_{PVD}$ chính xác (thường bố trí lưới hoa mai $1{,}0 \div 1{,}3\text{ m}$), đảm bảo độ cố kết đạt $U \ge 90%$ đúng tiến độ gia tải mà không bị vỡ phương án tài chính.
  • Khuyến nghị chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học để Tổng cục Đường bộ Việt Nam, Viện KHCN Xây dựng xem xét cập nhật, bổ sung các điều khoản chỉ dẫn kỹ thuật trong các tiêu chuẩn chuyên ngành như TCCS 41:2022 và TCVN 4200:2012 đối với quy trình thí nghiệm đất yếu ven biển.

Limitations và Future Research

  1. Giới hạn phạm vi tầng địa chất: Nghiên cứu tập trung chuyên sâu vào phân tầng sét mềm Holocen thượng $amQ_2^{3}tb$. Các phân tầng đất yếu nằm sâu hơn như sét hữu cơ hồ - đầm lầy Hải Hưng dưới ($mbQ_2^{1-2}hh_1$) hay sét biển Hải Hưng trên ($mQ_2^{1-2}hh_3$) mới chỉ được khảo sát tổng quan thạch học mà chưa thực hiện đầy đủ ma trận thí nghiệm kích thước mẫu.
  2. Giới hạn điều kiện gia tải phòng thí nghiệm: Các thí nghiệm Oedometer và Rowe cell được thực hiện dưới chế độ gia tải tĩnh từng cấp. Hành vi cố kết - từ biến dưới tác động của tải trọng động giao thông lặp chu kỳ cao (traffic cyclic loading) kết hợp áp lực nước lỗ rỗng dư tích lũy chưa được mô phỏng toàn diện.
  3. Chưa định lượng hóa đầy đủ hệ số cố kết thứ cấp ($C_\alpha$): Luận án chủ yếu tập trung vào giai đoạn cố kết sơ cấp ($U \le 100%$), thành phần lún từ biến thứ cấp lâu dài theo thời gian ($C_\alpha / C_c$) cần chuỗi thời gian thí nghiệm dài ngày hơn để xác lập mô hình hoàn chỉnh.

Định hướng nghiên cứu tương lai:

  • Mở rộng ma trận thí nghiệm kích thước mẫu cho toàn bộ các phân vị địa tầng sét yếu Đệ tứ ($Q_1^3vp, Q_2^{1-2}hh, Q_2^3tb$) dọc ven biển miền Trung và Nam Bộ.
  • Nghiên cứu ứng xử cố kết và suy giảm cường độ kháng cắt của đất nền dưới tác dụng của tải trọng động giao thông kết hợp quan trắc áp lực nước lỗ rỗng bằng piezometer hiện trường.
  • Ứng dụng mô hình số 3D phi tuyến nâng cao (Soft Soil Creep Model, Modified Cam-Clay trong PLAXIS 3D) tích hợp các thông số hiệu chỉnh kích thước mẫu để mô phỏng tương tác nền - công trình đắp quy mô lớn.

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

  • Tác động học thuật (Academic Impact): Luận án cung cấp hệ thống dữ liệu thực nghiệm chuẩn mực về đặc tính biến dạng đất yếu ĐBSH, dự kiến đóng góp vào các bài báo khoa học chất lượng cao trên các tạp chí chuyên ngành uy tín trong nước (Tạp chí Địa kỹ thuật, Tạp chí Xây dựng, Tạp chí Giao thông Vận tải) và quốc tế (Scopus/WoS), trở thành tài liệu tham khảo nền tảng cho các nghiên cứu sinh và học viên cao học ngành Kỹ thuật Địa chất và Địa kỹ thuật Công trình.
  • Chuyển dịch thực tiễn công nghiệp xây dựng: Cung cấp công cụ tính toán chính xác cho các tổng công ty tư vấn thiết kế giao thông, cảng biển (TEDI, Portcoast, CDC...) trong việc tối ưu hóa phương án xử lý nền đất yếu, giảm chi phí đầu tư xử lý bấc thấm PVD và cọc vật liệu rời từ $10 \div 20%$ nhờ loại bỏ hệ số an toàn dư thừa do thiếu dữ liệu thí nghiệm.
  • Thích ứng biến đổi khí hậu và an sinh xã hội: Đảm bảo độ ổn định lâu dài cho hệ thống đê biển kiên cố hóa chống nước biển dâng tại Thái Bình, Nam Định, Ninh Bình; ngăn ngừa sự cố sạt trượt mái dốc đường đắp và lún nứt mặt đường cao tốc ven biển, giảm thiểu chi phí bảo trì sửa chữa hàng năm của ngân sách nhà nước.

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

  • Nghiên cứu sinh và Giảng viên Địa kỹ thuật: Tiếp cận phương pháp luận thiết kế thí nghiệm đa kích thước mẫu mẫu mực, cơ sở lý thuyết cố kết hướng tâm và hệ dữ liệu cơ lý chi tiết của trầm tích Đệ tứ ĐBSH.
  • Kỹ sư Tư vấn Thiết kế và Địa kỹ thuật: Có bảng tra thông số thực nghiệm tin cậy ($c_h/c_v, c_v, C_c, p_c$) và quy trình tính lún chính xác cho công trình đất đắp trên nền đất yếu phân tầng.
  • Ban Quản lý Dự án và Nhà thầu Thi công: Tối ưu hóa tiến độ gia tải, xác định chính xác thời gian dỡ tải bù lún và cắm bấc thấm, loại bỏ rủi ro tranh chấp kỹ thuật do lún vượt dự báo.
  • Cơ quan Quản lý Nhà nước (Bộ Giao thông Vận tải, Bộ Xây dựng): Có cơ sở khoa học xác thực để rà soát, hiệu chỉnh và hoàn thiện hệ thống tiêu chuẩn kỹ thuật quốc gia.

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ì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc luận án đã làm sáng tỏ và lượng hóa được "hiệu ứng quy mô kích thước mẫu" đối với các đặc trưng cố kết của đất sét mềm bão hòa $amQ_2^{3}tb$ ven biển ĐBSH, từ đó mở rộng lý thuyết cố kết một chiều Terzaghi (1925) và lý thuyết dòng thấm hướng tâm Barron (1948) - Hansbo (1981). Nghiên cứu chỉ ra rằng giả định coi hệ số cố kết $c_v$ là hằng số bất biến đối với mọi kích thước hình học mẫu thí nghiệm trong phòng là không hoàn toàn chính xác; chiều cao mẫu $H$ gia tăng làm kéo dài đường thấm thoát nước nội tại, dẫn đến sự suy giảm có quy luật của $c_v$.

2. Sự đổi mới trong phương pháp nghiên cứu của luận án khi so sánh với ít nhất hai nghiên cứu quốc tế trước đây?

So với nghiên cứu của Kaothon et al. (2014) tại Thái Lan (chỉ cố định chiều cao $H = 20\text{ mm}$ và khảo sát 2 cỡ đường kính nhỏ $\Phi = 35, 60\text{ mm}$) và Al-Khafaji (2016) tại Anh (thí nghiệm trên sét nhân tạo Kaolin hoàn toàn đồng nhất), phương pháp của luận án có tính đổi mới vượt bậc:

  • Thiết kế ma trận thực nghiệm đồng thời biến thiên cả 2 chiều hình học ($\Phi = 61{,}8\text{ mm}, 79{,}8\text{ mm}$ kết hợp với 3 mức chiều cao $H = 20\text{ mm}, 30\text{ mm}, 40\text{ mm}$).
  • Thực hiện hoàn toàn trên mẫu đất nguyên trạng tự nhiên có cấu trúc vi mô trầm tích sông - biển phức tạp của đồng bằng châu thổ nhiệt đới gió mùa.
  • Tích hợp song song thí nghiệm nén đứng Oedometer và nén hướng tâm hộp Rowe cell đo áp lực nước lỗ rỗng biên để xác định đồng bộ tỷ số dị hướng $c_h/c_v$.

3. Phát hiện bất ngờ nhất trong quá trình thực nghiệm kèm minh chứng số liệu?

Phát hiện bất ngờ nhất là sự đối lập hoàn toàn về xu hướng biến thiên của hệ số cố kết $c_v$ theo chiều cao mẫu $H$ giữa đất sét mềm tự nhiên ĐBSH và đất sét Kaolin nhân tạo trong nghiên cứu của Al-Khafaji (2016). Trong khi Al-Khafaji ghi nhận $c_v$ tăng khi $H$ tăng, thì thực nghiệm của luận án trên đất sét $amQ_2^{3}tb$ khẳng định $c_v$ luôn giảm khi $H$ tăng từ $20\text{ mm}$ lên $40\text{ mm}$ ở tất cả các cấp gia tải $\sigma' = 12{,}5 \div 400\text{ kPa}$. Điều này được minh chứng bằng cơ chế ma sát dính kết giữa các hạt sét - khoáng Illite/Smectite và độ trễ tiêu tán áp lực nước lỗ rỗng thực tế trong môi trường lỗ rỗng siêu nhỏ của đất trầm tích tự nhiên.

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

Luận án cung cấp đầy đủ và chi tiết quy trình lặp lại thực nghiệm chuẩn hóa, bao gồm: bản vẽ gia công chế tạo bộ dao vòng thép không gỉ cải tiến ($\Phi = 61{,}8\text{ mm}, 79{,}8\text{ mm}$; $H = 20, 30, 40\text{ mm}$); quy trình cắt gọt mẫu nguyên trạng bằng bàn quay và dây cắt định hình; quy trình bão hòa mẫu ngập nước; các cấp áp lực gia tải chuẩn hóa ($\sigma' = 12{,}5; 25; 50; 100; 200; 400\text{ kPa}$); quy trình lắp đặt giếng cát trung tâm trong hộp nén Rowe cell và quy chuẩn đồ giải xác định $c_v, c_h, p_c$ theo Casagrande và Taylor.

5. Kế hoạch nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chương trình nghiên cứu dài hạn trong thập kỷ tới tập trung vào 3 mũi nhọn:

  • Năm 1 - 3: Mở rộng bản đồ hóa cơ sở dữ liệu số địa kỹ thuật các tầng đất yếu sâu ($mbQ_2^{1-2}hh_1, mQ_2^{1-2}hh_3$) cho toàn bộ dải duyên hải Bắc Bộ và Bắc Trung Bộ.
  • Năm 4 - 6: Phát triển thiết bị nén Oedometer chu kỳ tự động (Cyclic Oedometer) tích hợp cảm biến đo áp lực nước lỗ rỗng siêu nhỏ để mô phỏng tải trọng động giao thông tốc độ cao trên nền đất yếu.
  • Năm 7 - 10: Xây dựng hệ thống tiêu chuẩn thiết kế kháng lún thông minh thích ứng biến đổi khí hậu và tích hợp trí tuệ nhân tạo (AI/Machine Learning) trong dự báo trường chuyển vị lún thời gian thực dựa trên chuỗi dữ liệu quan trắc viễn thám InSAR và cảm biến IoT hiện trường.

Kết luận

  1. Phân vùng cấu trúc nền khoa học: Đã hệ thống hóa và phân chia thành công địa tầng đất yếu ven biển ĐBSH thành 3 kiểu cấu trúc nền đặc trưng, chỉ rõ kiểu CTN II là đối tượng chi phối chính gây lún lớn và mất ổn định cho các công trình đất đắp.
  2. Lượng hóa hiệu ứng kích thước mẫu nén: Khẳng định quy luật khoa học: khi chiều cao mẫu $H$ tăng từ $20\text{ mm}$ lên $40\text{ mm}$, hệ số cố kết thẳng đứng $c_v$ giảm ở mọi cấp áp lực nén; khi đường kính mẫu $\Phi$ tăng từ $61{,}8\text{ mm}$ lên $79{,}8\text{ mm}$, chỉ số nén $C_c$ giảm nhẹ và áp lực tiền cố kết $p_c$ tăng nhẹ.
  3. Xác lập quan hệ tương quan định lượng: Thiết lập thành công phương trình tương quan thực nghiệm giữa hệ số cố kết $c_v$ và chiều cao mẫu $H$ trong khoảng áp lực nén $\sigma' = 100 \div 200\text{ kPa}$ cho đất sét yếu $amQ_2^{3}tb$.
  4. Xác định chính xác tỷ số dị hướng thấm: Bằng thí nghiệm hộp nén Rowe, đã xác định được tỷ số hệ số cố kết ngang trên cố kết đứng $c_h/c_v = 1{,}45 \div 2{,}85$, thay thế hoàn toàn việc phỏng đoán định tính trong tính toán kỹ thuật.
  5. Đột phá trong ứng dụng tính toán công trình: Ứng dụng thành công hệ thông số hiệu chỉnh vào bài toán thiết kế bấc thấm PVD và kiểm toán ổn định trượt cho nền đắp giao thông, đê điều ven biển, nâng cao độ chính xác dự báo lún và tối ưu hóa hiệu quả kinh tế xây dựng.
  6. Mở ra các hướng nghiên cứu liên ngành: Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về cơ học đất từ biến, động lực học nền móng công trình ven biển thích ứng bền vững với biến đổi khí hậu toàn cầu.