Tổng quan về luận án

Nghiên cứu địa kỹ thuật công trình ven biển và các vùng đồng bằng ngập lũ tại Việt Nam luôn đối mặt với thách thức cơ học phức tạp từ các tầng trầm tích sét yếu đệ tứ Holocene ($Q_{IV}^{1,2}$ và $Q_{IV}^{2,3}$) có niên đại từ 2.000 đến 10.000 năm, chiều dày biến động từ 8 m đến 40 m. Các tầng sét này có độ ẩm tự nhiên ($w_n$) xấp xỉ hoặc vượt quá giới hạn chảy ($w_L$), hệ số rỗng ban đầu ($e_0$) rất lớn và sức kháng cắt không thoát nước ($s_u$) cực kỳ thấp. Để xử lý nền đất yếu phục vụ xây dựng hạ tầng cảng biển và khu công nghiệp quy mô lớn, giải pháp sử dụng bấc thấm đứng (Prefabricated Vertical Drains - PVD) kết hợp gia tải trước bằng đất đắp (Surcharge) hoặc gia tải hút chân không (Vacuum consolidation) là công nghệ chủ đạo. Tuy nhiên, một khoảng trống nghiên cứu (research gap) nghiêm trọng tồn tại trong suốt nhiều thập kỷ: các tiêu chuẩn thiết kế hiện hành tại Việt Nam như TCVN 4200:2012, TCVN 9355:2012, TCVN 9842:2013 và 22 TCN 262-2000 chủ yếu quy định thí nghiệm cố kết nén một chiều gia tải từng cấp truyền thống (Incremental Loading - IL). Phương pháp IL tiêu tốn từ 7 đến 10 ngày cho một mẫu thử, cung cấp các điểm dữ liệu rời rạc trên đường cong $e-\log \sigma'_v$, làm lu mờ điểm uốn xác định áp lực tiền cố kết ($\sigma'_c$), đồng thời các công thức tính toán hiện hành thường đơn giản hóa địa tầng nhiều lớp phức tạp về một lớp tương đương với các thông số trung bình cố định.

Luận án tiến sĩ kỹ thuật của NCS. Nguyễn Công Oanh với đề tài "Nghiên cứu ứng dụng một số thông số trong thí nghiệm cố kết tốc độ biến dạng không đổi (CRS) vào phân tích bài toán cố kết có sử dụng bấc thấm trong điều kiện Việt Nam" dưới sự hướng dẫn khoa học của GS. Trần Thị Thanh (2019) tại Viện Khoa học Thủy lợi Miền Nam đã giải quyết triệt để khoảng trống học thuật và thực tiễn này. Luận án đặt ra 3 câu hỏi nghiên cứu cốt lõi:

  1. Cơ chế biến dạng và ứng xử thoát nước lỗ rỗng của đất sét yếu Holocene Việt Nam dưới thí nghiệm nén cố kết tốc độ biến dạng không đổi (Constant Rate of Strain - CRS) diễn ra như thế nào ở tốc độ chuẩn hóa $\dot{\varepsilon} = 0,02%/\text{phút}$ so với thí nghiệm IL truyền thống?
  2. Mối tương quan định lượng giữa các thông số biến dạng phi tuyến ($C_{c1}, C_{c2}, C_r, \sigma'_c, c_v, c_h$) từ thí nghiệm CRS với các thí nghiệm hiện trường (CPTu, FVT, DST) được xác lập ra sao?
  3. Mô hình toán học sai phân hữu hạn (FDM) đa lớp cải tiến có khả năng phản ánh chính xác lịch sử gia tải theo thời gian, hiệu ứng vùng xáo động (smear zone), sức cản lõi thấm (well resistance) và chiều dài bấc thấm thay đổi ($L_{\text{PVD}} < H_{\text{clay}}$) hay không?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • H1: Đất sét yếu Holocene tại các đồng bằng ven biển Việt Nam tồn tại trạng thái quá cố kết tự nhiên nhẹ ($\text{OCR} > 1,0$) do quá trình già hóa (aging/thixotropy) và từ biến lịch sử với áp lực quá cố kết $\text{POP} = \sigma'c - \sigma'{v0}$ dao động từ 20 kPa đến 60 kPa;
  • H2: Áp lực tiền cố kết $\sigma'_c(\text{CRS})$ từ thí nghiệm CRS lớn hơn có quy luật từ 5% đến 25% (trung bình 16%, $R^2 = 0,85$) so với $\sigma'_c(\text{IL})$, và tốc độ biến dạng $\dot{\varepsilon} = 0,02%/\text{phút}$ đảm bảo tỷ số áp lực nước lỗ rỗng đáy mẫu $u_b/\sigma_v < 30%$, phản ánh chính xác trạng thái ứng suất thực tế tại hiện trường;
  • H3: Việc tích hợp các thông số CRS phi tuyến đa pha vào phần mềm sai phân hữu hạn CONSOPRO do tác giả tự lập trình sẽ triệt tiêu sai số dự báo độ lún và áp lực nước lỗ rỗng dư so với phương pháp đơn lớp truyền thống.

Khung lý thuyết của luận án tích hợp lý thuyết cố kết thấm một chiều phi tuyến của Terzaghi (1925), lý thuyết dòng thấm hướng tâm qua bấc thấm của Barron (1948) và Hansbo (1979, 1981, 1997, 2011), cơ sở lý thuyết thí nghiệm CRS của Wissa et al. (1971) cùng mô hình chuyển tiếp nén lún đa đoạn của Suzuki (2004). Nghiên cứu được kiểm chứng trên chuỗi dữ liệu thực địa quy mô lớn tại 8 dự án hạ tầng trọng điểm trải dài từ Bắc vào Nam: Cảng Hải Phòng Giai đoạn 2 (Chùa Vẽ), Cảng Quốc tế Cái Mép - ODA, Cảng SPCT Hiệp Phước, Nhà máy Khí Cà Mau, Cảng Quốc tế Thị Vải, Cần Giuộc (Long An), Bình Chánh và Đa Phước (TP.HCM). Kết quả chứng minh thí nghiệm CRS rút ngắn thời gian thử nghiệm trong phòng từ 7 ngày xuống còn 24–48 giờ (tiết kiệm hơn 75% thời gian), đồng thời nâng cao độ chính xác dự báo độ lún cố kết dài hạn đạt mức tương thích trên 90% so với hệ thống quan trắc hiện trường đa tầng.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp nén cố kết trong phòng bắt đầu từ các nghiên cứu nền tảng của Crawford (1964), Byrne (1969) và Smith (1969) nhằm khắc phục nhược điểm kéo dài thời gian của hộp nén Oedometer cổ điển. Wissa et al. (1971) đã chính thức đặt nền móng giải tích cho thí nghiệm tốc độ biến dạng không đổi (CRS), được chuẩn hóa quốc tế qua ASTM D4186-12 và JIS A1227-2009. Umehara (1975, 1979, 1983) mở rộng phương pháp xác định thông số cố kết có xét đến trọng lượng bản thân mẫu đất và ứng dụng thành công tính toán độ lún đê chắn sóng trên nền đất sét giàu bụi tại Nhật Bản. Ngược lại, Leroueil (1983, 1985, 1996) khi nghiên cứu ảnh hưởng của tốc độ biến dạng lên áp lực tiền cố kết ($\sigma'_y$) trên các loại sét Champlain (Canada) đã chỉ ra rằng ứng xử ứng suất - biến dạng chịu sự chi phối chặt chẽ của tốc độ gia tải, từ đó Leroueil (1983) và Korhonen (1997) lập luận hoài nghi rằng giá trị $\sigma'_c$ từ thí nghiệm IL truyền thống phản ánh sát hơn điều kiện hiện trường so với CRS.

Một cuộc tranh biện học thuật lớn thứ hai tập trung vào việc lựa chọn tốc độ biến dạng thích hợp. Gorman et al. (1976, 1981) đề xuất duy trì tỷ số $u_b/\sigma_v$ trong khoảng 3% đến 30%. Jia et al. (2010) khi thực hiện 114 thí nghiệm CRS trên đất sét Ariake (Nhật Bản) chứng minh tốc độ $\dot{\varepsilon} = 0,02%/\text{phút}$ mang lại sự tương thích cao nhất về $\sigma'_c$ và $c_v$ giữa CRS và IL. Chung et al. (2012) nghiên cứu tính dị hướng của đất sét Ariake chỉ ra tỷ số thấm ngang/đứng $k_h/k_v$ dao động quanh 1,65 ở tốc độ $0,02%/\text{phút}$. Ngược lại, Kassim et al. (2016) cảnh báo rằng tốc độ biến dạng quá thấp sẽ kích hoạt từ biến làm phóng đại hệ số cố kết $c_v(\text{CRS})$ bất hợp lý, khuyến nghị tốc độ chuẩn hóa $0,0125 - 0,100\text{ mm/phút}$ tùy thuộc hàm lượng hạt sét (35% - 68%). Về mặt lý thuyết cố kết hướng tâm qua bấc thấm, sau các lời giải kinh điển của Rendulic (1936), Biot (1941), Barron (1948), Yoshikuni & Nakanode (1974), Hansbo (1979, 1981, 1997, 2011) đã hoàn thiện công thức giải tích tính đến sức cản lõi thấm và vùng xáo động dựa trên giả thiết biến dạng bằng nhau (equal strain). Tuy nhiên, nghiệm của Hansbo bị giới hạn nghiêm ngặt trong điều kiện nền một lớp đồng nhất với tải trọng tức thời hoặc tăng tuyến tính.

Vị thế học thuật của luận án được định vị vững chắc khi giải quyết đồng thời hai nút thắt:

  1. Cung cấp luận cứ thực nghiệm bác bỏ quan điểm cho rằng $\sigma'_c(\text{CRS})$ không thể áp dụng cho tính toán thiết kế tại Việt Nam;
  2. Vượt qua giới hạn của lời giải Hansbo cổ điển bằng cách số hóa mô hình sai phân hữu hạn đa lớp nén lún phi tuyến, trực tiếp liên kết các thông số CRS với hệ thống thoát nước thẳng đứng có xét đến chiều dài bấc thấm hữu hạn ($L_{\text{PVD}} < H_{\text{clay}}$). So sánh với các nghiên cứu quốc tế điển hình như dự án Sân bay Quốc tế Kansai và cải tạo đất Isogo Yokohama của Suzuki (2004), đê biển Banjarmasin Indonesia của Suzuki (2008), và nghiên cứu đất sét Ariake của Jia et al. (2010), công trình của NCS. Nguyễn Công Oanh mở rộng phạm vi ứng dụng sang môi trường trầm tích châu thổ nhiệt đới gió mùa có tính phân lớp sâu sắc và độ nhạy cảm cao tại Việt Nam.

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

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

Luận án đã mở rộng và thách thức các giả định cơ bản trong cơ học đất cố kết cổ điển thông qua ba đóng góp lý thuyết nền tảng:

Hệ số rỗng (e)
 ^
         σ'v0  σ'c    (3-5)σ'c
  1. Mô hình hóa phi tuyến đường cong nén lún ba giai đoạn ($C_r - C_{c1} - C_{c2}$): Thách thức giả định truyền thống của Terzaghi và Casagrande vốn coi chỉ số nén $C_c$ là một hằng số duy nhất trong giai đoạn cố kết thường (NC). Dữ liệu CRS liên tục chỉ ra rằng đường cong $e-\log \sigma'_v$ có dạng hình chữ S đặc trưng: khi ứng suất vượt qua áp lực tiền cố kết $\sigma'c$, chỉ số nén tức thời $C{c1}$ tăng vọt và biến thiên rất mạnh trong khoảng ứng suất $\sigma'_c \le \sigma'_v \le (3 \div 5)\sigma'c$, sau đó độ dốc mới suy giảm và tiệm cận về một hằng số ổn định $C{c2}$. Việc áp dụng một giá trị $C_c$ duy nhất sẽ gây sai số lớn trong giai đoạn đầu của quá trình gia tải công trình.
  2. Xác lập bản chất quá cố kết thực tế của trầm tích Holocene: Trái ngược với các quan niệm trước đây cho rằng sét yếu bờ biển Việt Nam ở trạng thái chưa cố kết hoặc cố kết bình thường ($\text{OCR} \le 1,0$), luận án chứng minh mẫu nguyên dạng lấy bằng ống Piston luôn có $\text{OCR} > 1,0$ và $\text{POP} = \sigma'c - \sigma'{v0} = 20 \div 60\text{ kPa}$. Hiện tượng $\text{OCR} < 1,0$ trong các báo cáo khảo sát cũ hoàn toàn bắt nguồn từ sự xáo động mẫu khi sử dụng ống lấy mẫu thành mỏng Shelby thông thường (khẳng định nhận định của Takemura et al., 1997 và Bùi Tấn Mẫn, 2003).
  3. Cơ chế kiểm soát tốc độ biến dạng đối với thông số cố kết: Luận án mở rộng lý thuyết của Leroueil (1996) khi chứng minh rằng ở tốc độ $\dot{\varepsilon} = 0,02%/\text{phút}$, hệ số cố kết $c_v(\text{CRS})$ tiệm cận giá trị $c_v(\text{NC})$ của thí nghiệm IL nhưng cho phép đo đạc liên tục sự suy giảm của hệ số thấm $k_v$ và hệ số rỗng $e$ theo ứng suất hữu hiệu thực tế.

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

Khung phân tích của luận án kết hợp hài hòa ba lý thuyết lớn: Lý thuyết biến dạng một chiều phi tuyến của Wissa et al. (1971), Lý thuyết cố kết thấm hướng tâm có xét đến vùng xáo động của Hansbo (1981), và Phương pháp sai phân hữu hạn (FDM) đa lớp mở rộng từ cấu trúc của Suzuki (2004).

Phương trình vi phân chi phối quá trình tiêu tán áp lực nước lỗ rỗng dư $u(r,z,t)$ trong mô hình đơn nguyên (unit cell) có dạng:

$$\frac{\partial \bar{u}}{\partial t} = c_h \left( \frac{\partial^2 u}{\partial r^2} + \frac{1}{r}\frac{\partial u}{\partial r} \right) + c_v \frac{\partial^2 \bar{u}}{\partial z^2}$$

Với các điều kiện biên và biến số thực nghiệm được cập nhật động theo từng bước thời gian $\Delta t$ và từng phân tầng $\Delta z$:

  • Bán kính ảnh hưởng tương đương của bấc thấm: $D_e = 1,05 \times D$ (bố trí tam giác) hoặc $D_e = 1,13 \times D$ (bố trí hình vuông), với $D$ là khoảng cách bấc thấm.
  • Đường kính tương đương của bấc thấm dẹt: $d_w = \frac{2(a+b)}{\pi}$, với $a, b$ là chiều rộng và chiều dày của bản bấc.
  • Hệ số cản trở dòng chảy do xáo động và chiều dài bấc thấm: $\mu = \ln\left(\frac{n}{s}\right) + \frac{k_h}{k'_h}\ln(s) - \frac{3}{4} + \pi z (2l - z) \frac{k_h}{q_w}$, trong đó $n = \frac{D_e}{d_w}$, $s = \frac{d_s}{d_w}$, $k'_h$ là hệ số thấm ngang trong vùng xáo động, $q_w$ là khả năng thoát nước của bấc thấm.
  • Cập nhật thông số phi tuyến: Tại mỗi nút lưới sai phân, ứng suất hữu hiệu $\sigma'v(z,t) = \sigma{v0}(z) + \Delta\sigma_v(z,t) - \bar{u}(z,t)$. Nếu $\sigma'_v \le \sigma'_c$, mô đun biến dạng tính theo $C_r$; nếu $\sigma'_c < \sigma'_v \le (3\div 5)\sigma'c$, tính theo $C{c1}$; nếu $\sigma'_v > (3\div 5)\sigma'c$, tính theo $C{c2}$. Đồng thời hệ số cố kết $c_v, c_h$ tự động chuyển đổi từ trạng thái quá cố kết $c_v(\text{OC})$ sang trạng thái cố kết thường $c_v(\text{NC})$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng (Positivism) và khoa học kỹ thuật suy diễn - thực nghiệm (Empirical-Deductive Engineering Approach). Luận án triển khai thiết kế đa mức (multi-level design) kết hợp chặt chẽ giữa 3 cấp độ: Thử nghiệm trong phòng thí nghiệm kiểm soát cao, Khảo sát hiện trường chuyên sâu, và Quan trắc ngược (Back-analysis) trên các công trình thực tế.

Ma trận mẫu nghiên cứu bao gồm hàng trăm mẫu đất sét nguyên dạng được lấy tại các độ sâu từ 2 m đến 38 m tại 8 phân vùng địa kỹ thuật trọng điểm. Tiêu chuẩn chọn mẫu tuyệt đối: chỉ sử dụng mẫu lấy bằng ống lấy mẫu Piston thủy lực có màng ngăn chân không nhằm khống chế biến dạng thể tích ban đầu $\varepsilon_0$ khi đưa về ứng suất hiện trường $\sigma'_{v0}$ dưới mức 6% (mức phân loại A theo thang chất lượng mẫu Lunne et al., 1997).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu tuân thủ các tiêu chuẩn kỹ thuật nghiêm ngặt nhất:

  • Thí nghiệm CRS: Tiến hành theo ASTM D4186-12 / JIS A1227-2009 trên mẫu đường kính $D = 60\text{ mm}$, chiều cao $H = 25\text{ mm}$. Bôi trơn thành dao vòng bằng mỡ silicon khử chân không; bão hòa mẫu với áp lực ngược (back-pressure) duy trì ổn định ở mức 100 kPa; vận hành máy nén điều khiển biến dạng tự động với tốc độ cố định $\dot{\varepsilon} = 0,02%/\text{phút}$ ($0,005\text{ mm/phút}$). Đo liên tục tải trọng nén trục $\sigma_v$ và áp lực nước lỗ rỗng đáy mẫu $u_b$ thông qua đầu đo áp lực điện tử kết nối máy tính ghi nhận dữ liệu theo thời gian thực.
  • Thí nghiệm gia tải từng cấp IL: Thực hiện theo TCVN 4200:2012 với tỷ số gia tải $\Delta\sigma'/\sigma' = 1,0$, mỗi cấp tải duy trì 24 giờ, đo lún ở các thời điểm tiêu chuẩn để so sánh đối chứng trực tiếp trên các mẫu cùng độ sâu.
  • Thí nghiệm CPTu và Cắt cánh hiện trường (FVT): Thiết bị CPTu hình nón diện tích mũi $A_c = 10\text{ cm}^2$, góc đỉnh $60^\circ$, vận tốc xuyên $1 \div 2\text{ cm/s}$, đo đồng thời sức kháng mũi $q_c$, ma sát hông $f_s$ và áp lực nước lỗ rỗng tại vai côn $u_2$. Thí nghiệm cắt cánh FVT tự khoan lỗ quay với tốc độ chuẩn $6^\circ/\text{phút}$ xác định sức kháng cắt không thoát nước hiện trường $s_{uf}$.
  • Hệ thống quan trắc hiện trường: Bố trí đồng bộ bao gồm bàn đo lún mặt (Settlement Plate - SP), vòng từ đo lún sâu đa tầng (Magnetic Extensometer - EXT), đầu đo áp lực nước lỗ rỗng kiểu dây rung (Vibrating Wire Piezometer - PIEZ) và hộp đo áp lực đất (Earth Pressure Cell - EPC).

Data và phân tích

Đặc trưng trầm tích tại các khu vực nghiên cứu thể hiện qua biểu đồ dẻo Casagrande với các điểm dữ liệu phân bố bám sát đường A-line ($PI = 0,73(LL - 20)$):

  • Cái Mép: Hàm lượng sét từ 60% đến 80%, độ ẩm tự nhiên $w_n = 70 \div 110%$, giới hạn chảy $w_L = 75 \div 115%$, chỉ số dẻo $PI = 45 \div 70%$.
  • Hải Phòng: Hàm lượng sét thấp hơn ($35 \div 45%$), hàm lượng bụi cao ($45 \div 55%$), $w_n = 45 \div 65%$.
  • Hiệp Phước & Cà Mau: Tổng hàm lượng bụi và sét ($d < 0,075\text{ mm}$) chiếm từ 90% đến 100%, tạo nên môi trường nén lún cực kỳ nhạy cảm.

Phần mềm CONSOPRO (đã được đăng ký bản quyền tác giả) được lập trình trên nền tảng ngôn ngữ kỹ thuật Fortran/Visual Basic, cấu trúc thành các module xử lý chuyên biệt:

  • Module STRAIN01: Tính toán độ biến dạng nén lún đơn vị $\Delta\varepsilon_i$ cho từng phân tầng đất dựa trên 3 độ dốc $C_r, C_{c1}, C_{c2}$;
  • Module STRAIN01A và STRAIN01B: Xử lý riêng biệt các phân lớp nằm trong vùng cắm bấc thấm ($z \le L_{\text{PVD}}$) và các phân lớp nằm dưới chân bấc thấm ($z > L_{\text{PVD}}$);
  • Module STRAIN02: Tính toán chuyển vị lún tổng cộng và điều chỉnh gradient áp lực cục bộ;
  • Module UNGSUATTRUNGBINH: Cập nhật phân bố ứng suất hữu hiệu bình quân $\bar{\sigma}'_v(z,t)$ theo từng bước thời gian;
  • Module HESORONG: Tự động tính toán sự suy giảm của hệ số rỗng $e(z,t)$ và cập nhật hệ số thấm $k_v, k_h$ tương ứng.

Các phương trình hồi quy thực nghiệm được thiết lập cho đất sét yếu Việt Nam với độ tin cậy thống kê cao:

  1. Tương quan chỉ số nén ổn định $C_{c2}$: $$C_{c2} = 0,01(w_n - 5) \quad \text{hoặc} \quad C_{c2} = 0,0125 w_n \quad (R^2 = 0,88)$$ $$C_{c2} = 0,0115(w_L - 14) \quad \text{và} \quad C_{c2} = 0,016 PI$$
  2. Tương quan áp lực tiền cố kết giữa hai phương pháp: $$\sigma'_c(\text{CRS}) = 1,05 \div 1,25 , \sigma'_c(\text{IL}) \quad (\text{Giá trị bình quân: } \sigma'_c(\text{CRS}) = 1,16 , \sigma'_c(\text{IL}), , R^2 = 0,85)$$
  3. Tương quan sức kháng cắt không thoát nước chuẩn hóa: $$\frac{s_{uf}}{\sigma'{v0}} = \left(\frac{s{un}}{\sigma'{v0}}\right){\text{NC}} \times \text{OCR}^m = 0,22 \times \text{OCR}^{0,85}$$
  4. Tương quan sức kháng xuyên mũi côn CPTu hiệu chỉnh: $$\sigma'c = \frac{q_T - \sigma{v0}}{N_{kt}} \quad \text{với hệ số nón } N_{kt} = 10 \div 20 \text{ (trung bình } 13,5\text{)}$$

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

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

      CRS (Tốc độ ε_dot = 0.02%/phút)              IL (Gia tải từng cấp 24h)
  1. Minh chứng thực nghiệm về sự vượt trội của $\sigma'_c(\text{CRS})$: Luận án chứng minh một cách thuyết phục rằng giá trị $\sigma'_c(\text{CRS})$ cao hơn $\sigma'_c(\text{IL})$ trung bình 16% không phải là một sai số kỹ thuật hay "ảo giác thí nghiệm" do tốc độ nhanh gây ra, mà là giá trị đại diện trung thực cho ứng suất tiền cố kết tại hiện trường. Trong thí nghiệm IL, khoảng thời gian chờ 24 giờ cho mỗi cấp tải đã làm biến dạng từ biến thứ cấp (secondary creep) tích tụ sớm ngay trong giai đoạn nén đàn hồi, khiến điểm gãy $\sigma'_c$ bị dịch chuyển sang trái và suy giảm cường độ giả tạo.
  2. Khám phá quy luật biến thiên của chỉ số nén dẻo $C_{c1}$: Luận án phát hiện hiện tượng phá hủy cấu trúc đất (structural breakdown) nguyên sinh diễn ra dữ dội ngay sau khi ứng suất vượt quá $\sigma'c$, đẩy hệ số nén $C{c1}$ lên các giá trị cực đại trước khi nén chặt dần để ổn định tại $C_{c2}$. Nếu mô hình hóa bằng một giá trị $C_c$ không đổi sẽ làm sai lệch tốc độ lún giai đoạn đầu của công trình từ 20% đến 35%.
  3. Hiện tượng tiêu tán áp lực nước lỗ rỗng hai pha trong bài toán bấc thấm lửng ($L_{\text{PVD}} < H_{\text{clay}}$): Tại các phân lớp có bấc thấm ($z \le L_{\text{PVD}}$), nước lỗ rỗng thoát theo phương ngang hướng tâm vào bấc thấm với tốc độ rất nhanh, áp lực nước lỗ rỗng dư tiêu tán trên 80% chỉ sau 6 đến 9 tháng. Tuy nhiên, tại tầng sét nằm dưới chân bấc thấm ($z > L_{\text{PVD}}$), nước lỗ rỗng bị buộc phải thoát theo phương thẳng đứng với đường thấm kéo dài, tạo ra hiện tượng ứ đọng áp lực nước lỗ rỗng dư cục bộ kéo dài hàng năm, là nguyên nhân gốc rễ gây lún trồi và lún dư sau thi công.
  4. Sự gia tăng sức kháng cắt không thoát nước sau gia tải: Cường độ kháng cắt $s_u$ đo lại bằng FVT và DST2 sau quá trình hút chân không và đắp gia tải tại Cảng Cái Mép và SPCT Hiệp Phước tăng từ $12 \div 18\text{ kPa}$ lên tới $45 \div 68\text{ kPa}$, hoàn toàn phù hợp với đường dự báo tăng trưởng sức kháng cắt theo mô hình $\Delta s_u = 0,22 \times U% \times \Delta\sigma'_v$ được cài đặt trong CONSOPRO.

Implications đa chiều

  • Về mặt học thuật và phương pháp luận: Thiết lập một chuẩn mực mới trong việc xử lý số liệu thí nghiệm cơ học đất tại Việt Nam; cung cấp một chuỗi phương trình hồi quy địa phương có độ tin cậy cao kết nối giữa thí nghiệm CRS trong phòng với các thí nghiệm hiện trường tiên tiến CPTu, FVT.
  • Về mặt kỹ thuật thực hành: Cho phép các đơn vị tư vấn thiết kế rút ngắn thời gian khảo sát - thí nghiệm phục vụ lập báo cáo khả thi (FS) và thiết kế kỹ thuật từ hàng tháng xuống vài tuần; loại bỏ hiện tượng tính thừa chiều cao đắp gia tải do đánh giá thấp áp lực tiền cố kết $\sigma'_c$.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học vững chắc để Bộ Xây dựng và Bộ Giao thông Vận tải xem xét bổ sung tiêu chuẩn quốc gia về "Phương pháp thí nghiệm nén cố kết tốc độ biến dạng không đổi CRS" và cập nhật quy trình tính toán xử lý đất yếu nhiều lớp trong các tiêu chuẩn TCVN liên quan.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Chưa tích hợp biến dạng từ biến thứ cấp (Creep): Luận án tập trung chuyên sâu vào bài toán cố kết sơ cấp thoát nước lỗ rỗng, chưa cài đặt phương trình từ biến phi tuyến ($C_\alpha$) theo thời gian dài hạn (hàng chục năm sau khi áp lực nước lỗ rỗng đã tiêu tán hoàn toàn).
  2. Giới hạn trong không gian biến dạng một chiều (1D Strain): Mô hình toán FDM trong CONSOPRO giả định biến dạng thuần túy theo phương đứng, chưa xét đến biến dạng trượt ngang (lateral displacement) tại khu vực chân taluy đắp của các công trình đắp cao.
  3. Mô hình hóa khả năng thoát nước của bấc thấm: Giả thiết sức cản lõi thấm $q_w$ là hằng số hoặc tuyến tính, chưa mô phỏng chi tiết hiện tượng gập bấc thấm (kinking/clogging) dưới biến dạng lún cực lớn (> 2 m).
  4. Cơ chế suy giảm hệ số thấm trong vùng xáo động: Tỷ số $k_h/k'_h$ và đường kính vùng xáo động $d_s$ được lấy cố định theo kinh nghiệm ($d_s = (2 \div 3)d_w$, $k_h/k'_h = 2 \div 5$), chưa xét đến sự biến đổi động của vi cấu trúc hạt sét quanh mỏ neo bấc thấm trong quá trình thi công cắm bấc.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) được định hình gồm 5 hướng:

  • Mở rộng mô hình CONSOPRO lên hệ không gian 2D/3D tích hợp phần tử hữu hạn phi đàn dẻo từ biến (Visco-elastoplastic FEM);
  • Nghiên cứu cơ chế suy giảm khả năng thoát nước $q_w$ của bấc thấm dẹt chịu biến dạng uốn nếp lớn trong điều kiện gia tải chân không sâu;
  • Xây dựng phần mềm trí tuệ nhân tạo (Machine Learning) kết hợp mạng nơ-ron sâu để tự động ánh xạ dữ liệu CPTu liên tục sang biểu đồ thông số CRS cho toàn bộ vùng châu thổ;
  • Đánh giá ảnh hưởng của biến đổi khí hậu và xâm nhập mặn đến tính chất hóa lý và cấu trúc cố kết của sét yếu Holocene;
  • Thiết lập hệ thống trạm quan trắc IoT truyền dữ liệu độ lún và áp lực nước lỗ rỗng thời gian thực tại các siêu dự án cảng biển để hiệu chuẩn liên tục mô hình số.

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

  • Tác động học thuật: Đóng góp các công trình khoa học chất lượng cao công bố tại các hội nghị địa kỹ thuật quốc tế uy tín (có chỉ số ISBN) và các tạp chí chuyên ngành (có chỉ số ISSN); cung cấp một công cụ tính toán mã nguồn độc lập CONSOPRO phục vụ công tác đào tạo sau đại học và nghiên cứu chuyên sâu tại Viện Khoa học Thủy lợi Việt Nam.
  • Chuyển đổi công nghiệp hạ tầng: Ứng dụng thành công tại các dự án trọng điểm quốc gia như Cảng Container Quốc tế Cái Mép - ODA (đón tàu trọng tải 200.000 DWT), Cảng SPCT Hiệp Phước, Nhà máy Khí Cà Mau và Cảng Hải Phòng Giai đoạn 2. Việc dự báo chính xác diễn biến độ lún và áp lực nước lỗ rỗng giúp nhà thầu (như TOA Corporation) và các Ban Quản lý Dự án (PMU 85, PMU Hàng Hải 2) kiểm soát an toàn tốc độ đắp gia tải, dỡ tải đúng thời điểm, tiết kiệm hàng chục tỷ đồng chi phí đắp bù lún và rút ngắn tiến độ thi công toàn dự án từ 2 đến 4 tháng.
  • Lợi ích xã hội và môi trường: Việc tối ưu hóa chiều sâu cắm bấc thấm và khối lượng đất đắp gia tải giúp giảm thiểu hàng triệu mét khối cát san lấp khai thác từ lòng sông, góp phần bảo vệ tài nguyên thiên nhiên và hệ sinh thái ven sông ngòi đồng bằng sông Cửu Long.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một hệ số liệu thực nghiệm chuẩn mực, phương pháp luận nghiên cứu kết hợp FDM và quy trình thí nghiệm CRS hiện đại; kế thừa mã nguồn và thuật toán CONSOPRO để phát triển các mô hình cố kết bậc cao.
  • Kỹ sư Địa kỹ thuật và Đơn vị Tư vấn Thiết kế: Sở hữu cẩm nang tra cứu tương quan địa kỹ thuật tin cậy ($C_{c1}, C_{c2}, C_r, \sigma'_c, s_u$) cho các vùng đất sét yếu trọng điểm Việt Nam; có phương pháp luận vững chắc để áp dụng thí nghiệm CRS vào hồ sơ thiết kế xử lý nền đất yếu bằng bấc thấm kết hợp hút chân không.
  • Chủ đầu tư và Ban Quản lý Dự án Giao thông/Thủy lợi: Giảm thiểu rủi ro pháp lý và tài chính phát sinh từ sự cố lún nứt công trình sau khi đưa vào vận hành; tối ưu hóa chi phí đầu tư và chi phí duy tu bảo dưỡng định kỳ nền bãi, cầu cảng.
  • Cơ quan Quản lý Nhà nước (Bộ Xây dựng, Bộ NN&PTNT, Bộ GTVT): Có căn cứ khoa học xác thực để thẩm định các giải pháp xử lý nền đất yếu phức tạp; làm tiền đề biên soạn và ban hành Tiêu chuẩn Quốc gia mới về thí nghiệm CRS.

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 lý thuyết độc đáo nhất là việc mở rộng nghiệm giải tích cố kết hướng tâm của Hansbo (1981) thành mô hình sai phân hữu hạn đa tầng nén lún phi tuyến ba giai đoạn ($C_r - C_{c1} - C_{c2}$). Luận án đã phá vỡ giả định đơn giản hóa nền đồng nhất một lớp của Hansbo, tích hợp sự biến thiên phi tuyến liên tục của hệ số cố kết $c_v(z,t), c_h(z,t)$, hệ số rỗng $e(z,t)$ và áp lực quá cố kết $\text{POP}(z)$ theo chiều sâu và thời gian, giải quyết trọn vẹn bài toán bấc thấm lửng ($L_{\text{PVD}} < H_{\text{clay}}$) mà lý thuyết Hansbo cổ điển không thể mô phỏng chính xác.

2. Đổi mới phương pháp luận của luận án thể hiện qua điểm nào khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó?

So với nghiên cứu của Leroueil (1983, 1985) trên sét Champlain (vốn cho rằng CRS làm sai lệch áp lực tiền cố kết do hiệu ứng tốc độ và khuyến nghị dùng IL) và nghiên cứu của Jia et al. (2010) trên sét Ariake (chỉ dừng lại ở phạm vi đối sánh thống kê trong phòng thí nghiệm):

  • Luận án của NCS. Nguyễn Công Oanh đã thực hiện một bước tiến phương pháp luận vượt bậc khi xây dựng một quy trình khép kín từ khâu lấy mẫu Piston chất lượng cao ($\varepsilon_0 < 6%$), chuẩn hóa tốc độ biến dạng $\dot{\varepsilon} = 0,02%/\text{phút}$ ($u_b/\sigma_v < 30%$), trích xuất trực tiếp các thông số CRS đưa vào phần mềm FDM tự phát triển (CONSOPRO), và kiểm chứng chéo với chuỗi số liệu quan trắc hiện trường đa thông số (bàn đo lún mặt, nhện từ đo lún sâu, piezometer đo áp lực rỗng, hộp đo áp lực đất) tại 4 đại dự án thực tế ở Việt Nam.

3. Phát hiện bất ngờ nhất trong quá trình nghiên cứu là gì và có số liệu thực nghiệm nào chứng minh?

Phát hiện bất ngờ nhất là toàn bộ các lớp sét yếu trầm tích Holocene tại miền Nam và miền Bắc Việt Nam đều ở trạng thái quá cố kết tự nhiên ($\text{OCR} > 1,0$) với áp lực quá cố kết $\text{POP} = 20 \div 60\text{ kPa}$, hoàn toàn trái ngược với niềm tin phổ biến trước đây của nhiều kỹ sư thực hành cho rằng đất sét châu thổ Việt Nam chưa cố kết xong ($\text{OCR} < 1,0$). Dữ liệu thực nghiệm chứng minh:

  • 100% mẫu lấy bằng ống Piston đạt cấp chất lượng cao đều cho $\text{OCR}$ dao động từ $1,15$ đến $2,50$ (Hình 2-3).
  • Áp lực tiền cố kết $\sigma'c(\text{CRS})$ luôn nằm cao hơn đường ứng suất hữu hiệu bản thân $\sigma'{v0}$ một khoảng cách biệt rõ ràng: $\Delta\sigma' = 45\text{ kPa}$ tại Hải Phòng, $30\text{ kPa}$ tại Cái Mép, Hiệp Phước, Bình Chánh, Đa Phước, và $20\text{ kPa}$ tại Cà Mau (Hình 2-2).

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

Có, luận án cung cấp quy trình tái lập cực kỳ chi tiết tại Bảng 1-1, Bảng 1-3 và Mục 1.1:

  • Quy định đường kính dao vòng $60\text{ mm}$, chiều cao mẫu $25\text{ mm}$, bôi trơn thành bằng dầu chân không chuyên dụng, lắp giấy thấm và đá thấm khử bọt khí hai mặt;
  • Áp lực bão hòa ngược (back-pressure) chuẩn xác 100 kPa;
  • Tốc độ biến dạng nén trục cố định $\dot{\varepsilon} = 0,02%/\text{phút}$;
  • Quy chuẩn đầu đo áp lực đáy mẫu $u_b$ đảm bảo tỷ số $u_b/\sigma_v < 30%$;
  • Cấu trúc các module tính toán và sơ đồ khối thuật toán của phần mềm CONSOPRO (STRAIN01, STRAIN01A, STRAIN01B, STRAIN02, UNGSUATTRUNGBINH, HESORONG) được trình bày minh bạch tại Chương 3 (Hình 3-2 đến Hình 3-6), cho phép các nhà nghiên cứu độc lập lập trình lại toàn bộ hệ thống giải thuật.

5. Chương trình 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 10 năm được vạch ra với 4 trụ cột chiến lược:

  1. Năm 1-3: Phát triển hoàn thiện phần mềm CONSOPRO-3D tích hợp phần tử hữu hạn phi đàn dẻo và từ biến thứ cấp ($C_\alpha$);
  2. Năm 4-5: Nghiên cứu thực nghiệm ảnh hưởng của nhiệt độ và xử lý chân không sâu kết hợp gia nhiệt đến động học tiêu tán áp lực nước lỗ rỗng trong bấc thấm;
  3. Năm 6-8: Số hóa bản đồ địa kỹ thuật thông minh cho toàn vùng đồng bằng sông Cửu Long và đồng bằng sông Hồng, tích hợp AI để dự báo tức thời các thông số CRS từ dữ liệu CPTu số;
  4. Năm 9-10: Xây dựng và đề xuất ban hành Tiêu chuẩn Quốc gia hoàn chỉnh về "Thiết kế xử lý nền đất yếu bằng bấc thấm đa tầng kết hợp gia tải chân không dựa trên thông số thí nghiệm CRS".

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Công Oanh đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra, đánh dấu một bước chuyển dịch mô thức (paradigm shift) quan trọng trong lĩnh vực Địa kỹ thuật Xây dựng tại Việt Nam qua 6 đóng góp cốt lõi:

  1. Xác lập cơ sở khoa học và thực tiễn vững chắc chứng minh tính khả thi, tính chính xác và hiệu quả vượt trội của thí nghiệm cố kết tốc độ biến dạng không đổi (CRS) theo chuẩn ASTM D4186-12 trong việc xác định các chỉ tiêu cơ lý đất sét yếu Holocene tại Việt Nam.
  2. Chứng minh thực nghiệm bản chất quá cố kết tự nhiên ($\text{OCR} > 1,0$, $\text{POP} = 20 \div 60\text{ kPa}$) của sét yếu bờ biển Việt Nam, đập tan nhận định sai lầm tồn tại lâu nay về trạng thái chưa cố kết do ảnh hưởng của xáo động mẫu.
  3. Thiết lập hệ thống công thức tương quan thực nghiệm có độ tin cậy cao ($R^2 = 0,85 \div 0,88$) giữa các thông số nén phi tuyến CRS ($C_{c1}, C_{c2}, C_r, \sigma'_c, c_v$) với các chỉ tiêu cơ bản ($w_n, w_L, PI$) và các thí nghiệm hiện trường CPTu, FVT.
  4. Cải tiến thành công lời giải giải tích cổ điển của Hansbo (1981), xây dựng mô hình toán sai phân hữu hạn đa lớp nén lún phi tuyến ba giai đoạn phản ánh sát thực tế địa tầng và lịch sử gia tải theo thời gian.
  5. Phát triển và đăng ký bản quyền phần mềm chuyên dụng CONSOPRO, giải quyết xuất sắc bài toán cố kết phức hợp có bấc thấm dẹt kết hợp gia tải đất đắp/hút chân không và bài toán bấc thấm lửng ($L_{\text{PVD}} < H_{\text{clay}}$).
  6. Kiểm chứng thành công mô hình trên 4 siêu dự án hạ tầng cảng biển và năng lượng trọng điểm quốc gia, mang lại giá trị kinh tế - kỹ thuật to lớn, rút ngắn tiến độ thí nghiệm 75% và nâng cao độ an toàn ổn định công trình.

Công trình của NCS. Nguyễn Công Oanh đặt nền móng học thuật vững chắc, mở ra các hướng nghiên cứu hiện đại về cơ học đất phi tuyến và để lại di sản thực tiễn lâu dài cho ngành xây dựng công trình trên nền đất yếu tại Việt Nam và khu vực Đông Nam Á.