Luận án tiến sĩ: Tính chất cơ lý nền đất yếu Quảng Nam - Đà Nẵng

Tổng quan về luận án

Vùng đồng bằng ven biển Quảng Nam - Đà Nẵng giữ vị thế hạt nhân kinh tế, đô thị và du lịch trọng điểm của khu vực miền Trung và hành lang kinh tế Đông - Tây. Tuy nhiên, sự phát triển hạ tầng giao thông kết nối tại đây đối mặt với thách thức địa kỹ thuật phức tạp do sự hiện diện của các tầng trầm tích đất yếu Đệ Tứ đa nguồn gốc (trầm tích sông - biển - đầm lầy, vũng vịnh). Các tầng đất này có đặc tính nén lún lớn, khả năng thoát nước kém, cường độ chịu tải cực thấp, dẫn đến nguy cơ lún kéo dài và mất ổn định mái taluy nền đường đắp.

+-------------------------------------------------------------------------------+
|             BỐI CẢNH ĐỊA KỸ THUẬT VEN BIỂN QUẢNG NAM - ĐÀ NẴNG                |
|  - Trầm tích trẻ Đệ Tứ đa nguồn gốc: ambQ2(3), mbQ2(2), ambQ2(2), abmQ2(1)... |
|  - Mực nước ngầm nông (< 2m chiếm ưu thế), ăn mòn carbonic (79.67%)           |
+-------------------------------------------------------------------------------+
                                      │
                                      ▼
+-------------------------------------------------------------------------------+
|                    KHOẢNG TRỐNG NGHIÊN CỨU & THỰC TRẠNG                       |
|  - Thiếu hệ thống cơ sở dữ liệu cơ lý đồng bộ đặc thù miền Trung              |
|  - Thiết kế thường giả định thô sơ: Ch = (2 ÷ 5)Cv, bỏ qua cấu trúc dị hướng  |
|  - Hệ quả: Lún dư vượt hạn mức, sạt trượt mái taluy, lãng phí kinh phí gia cố |
+-------------------------------------------------------------------------------+
                                      │
                                      ▼
+-------------------------------------------------------------------------------+
|                      MỤC TIÊU & GIẢI PHÁP ĐỘT PHÁ                             |
|  - Xác lập bản chất cơ lý: Thành phần khoáng vật, hóa học, hữu cơ             |
|  - Định lượng chính xác Cv, Ch, tỷ số m = Ch/Cv (Thí nghiệm + Asaoka)         |
|  - Phân vùng cấu trúc nền 1/50.000 & tối ưu hóa giải pháp PVD, Giếng cát, Ráo |
+-------------------------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng các nghiên cứu địa kỹ thuật chuyên sâu về đất yếu ở Việt Nam trước đây hầu như chỉ tập trung tại Đồng bằng Bắc Bộ (ĐBBB) và Đồng bằng sông Cửu Long (ĐBSCL). Do điều kiện địa hình có độ dốc và mức độ chia cắt lớn hơn, chế độ thủy triều - hải văn đặc thù, đất yếu vùng ven biển Quảng Nam - Đà Nẵng có chiều dày biến động mạnh, diện phân bố đan xen phức tạp và hiếm khi lộ hẳn trên bề mặt. Trong thực tiễn thiết kế xử lý nền đường, các đơn vị tư vấn thường áp dụng quy trình 22TCN 262-2000 với các giả định kinh nghiệm thiếu kiểm chứng địa phương, đặc biệt là việc tạm ước tính hệ số cố kết thấm ngang qua hệ số cố kết thấm đứng $C_h = (2 \div 5)C_v$, hoặc thiếu vắng các thông số đầu vào quan trọng như $S_u, E_0, P_c, k_v, k_h, C_v, C_h, m=C_h/C_v, C_c, C_s$. Hậu quả là nhiều công trình huyết mạch (đường tránh Đà Nẵng, đường ven sông Tuyên Sơn - Túy Loan, đường gom cao tốc Đà Nẵng - Quảng Ngãi) xuất hiện hiện tượng trượt lở mái taluy và độ lún dư sau thi công vượt giới hạn cho phép.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

• RQ1: Nguồn gốc trầm tích, thành phần khoáng vật và hàm lượng vật chất hữu cơ chi phối tính chất cơ lý và khả năng tự nén chặt của đất yếu ven biển Quảng Nam - Đà Nẵng như thế nào?
  • H1: Đất yếu khu vực chịu sự kiểm soát của khoáng vật Kaolinit và Illit kết hợp độ ẩm cao, có hệ số tự nén chặt $k_d < 0.8$, phản ánh giai đoạn đang nén chặt dở dang.
• RQ2: Tỷ số cố kết thấm dị hướng $m = C_h/C_v$ biến thiên ra sao theo nguồn gốc địa tầng và cấp áp lực nén, và phương pháp phân tích ngược quan trắc lún có phản ánh chính xác hơn thí nghiệm trong phòng?
  • H2: Tỷ số $m = C_h/C_v$ không phải là hằng số cố định $(2 \div 5)$ mà phụ thuộc nghiêm ngặt vào cấu trúc trầm tích (biến thiên từ 1.3 đến 7.6), trong đó phân tích ngược Asaoka từ dữ liệu hiện trường tiệm cận ứng xử thực tế của nền xử lý tiêu thoát nước thẳng đứng.
• RQ3: Làm thế nào để phân chia cấu trúc nền (CTN) đất yếu một cách khoa học phục vụ tối ưu hóa lựa chọn giải pháp xử lý kỹ thuật?
  • H3: Tiếp cận theo lý thuyết Địa hệ Tự nhiên - Kỹ thuật cho phép tích hợp chiều dày, thế nằm và tính chất cơ lý thành các đơn vị cấu trúc nền chuyên biệt, chuẩn hóa bài toán thiết kế bấc thấm (PVD) và giếng cát (SD).

Khung lý thuyết của luận án được xây dựng dựa trên lý thuyết cố kết thấm một chiều của Karl Terzaghi (1925), lý thuyết cố kết bán kính của Barron (1948), định lý cố kết ba chiều của N. Carrillo (1942), và trường phái phân tích Địa hệ Tự nhiên - Kỹ thuật (Natural-Technical Geo-systems) của E.M. Sergeev và Bondarik. Phạm vi nghiên cứu khống chế dải đồng bằng ven biển kéo dài từ quận Liên Chiểu (Đà Nẵng) đến huyện Núi Thành (Quảng Nam), với chiều sâu khảo sát đến 30 m. Đóng góp đột phá của luận án là xác lập bộ dữ liệu cơ lý chuẩn tắc, xây dựng bản đồ địa kỹ thuật và cấu trúc nền tỷ lệ 1/50.000, giúp giảm sai số dự báo lún cố kết từ 20-30% so với phương pháp tính toán truyền thống.

Literature Review và Positioning

Nghiên cứu cơ học đất và xử lý nền đất yếu trên thế giới đã trải qua gần một thế kỷ phát triển với các mốc tiến hóa quan trọng. Khởi đầu từ mô hình cố kết một trục kinh điển của Karl Terzaghi (1925), lý thuyết này được mở rộng bởi Florin (1951) với mô hình "lực thể tích" có tính tổng quát cao, xét đến độ bền kiến trúc và tính nén của nước lỗ rỗng. Về bài toán tiêu thoát nước thẳng đứng kết hợp gia tải trước, N. Carrillo (1942) đã đặt nền móng bằng việc phân rã phương trình vi phân ba chiều thành tổ hợp nghiệm cố kết phẳng và hướng kính, tạo tiền đề cho Barron (1948) và sau đó là Hansbo (1981) hoàn thiện lời giải giải tích có xét đến hiệu ứng vùng xáo động quanh bấc thấm (smear zone) và sức cản dòng chảy qua lõi thoát nước (well resistance).

+─────────────────────────────────────────────────────────────────────────────+
|               TIẾN TRÌNH LÝ THUYẾT ĐỊA KỸ THUẬT & CỐ KẾT ĐẤT YẾU           |
+─────────────────────────────────────────────────────────────────────────────+
  1925: Terzaghi ────► Lý thuyết cố kết thấm 1D (Oedometer, cv)
  1938: Casagrande ──► Xác định áp lực tiền cố kết Pc, đường cong log(t)
  1942: Carrillo ────► Tách bài toán cố kết 3D = 1D đứng + bán kính ngang
  1948: Barron ──────► Lời giải cố kết bán kính cho giếng cát / PVD
  1957: Skempton ────► Tương quan Su/σ'v0 với chỉ số dẻo PI
  1974: Ladd & Foott ─► Phương pháp SHANSEP xác định sức kháng cắt theo OCR
  1981: Hansbo ──────► Bổ sung hiệu ứng xáo động (smear zone: ds/dw, kh/ks)
  2002: Bergado/Seah ─► Đo đạc Ch, m=Ch/Cv thực tế trên sét mềm Bangkok
  2018: Luận án Yến ─► Chuẩn hóa Ch, m=Ch/Cv cho trầm tích ĐBVB Quảng Nam - ĐN
+─────────────────────────────────────────────────────────────────────────────+

Trong đánh giá đặc tính nén lún và sức kháng cắt, Casagrande (1938) đã thiết lập phương pháp đồ thị xác định áp lực tiền cố kết $P_c$ trên hệ tọa độ $e - \log\sigma'$, xác định hệ số nén lún $C_c$ và hệ số cố kết $C_v$. Skempton (1957) đề xuất mối quan hệ thực nghiệm kinh điển giữa tỷ số sức kháng cắt không thoát nước trên ứng suất hữu hiệu tự nhiên với chỉ số dẻo: $$\frac{c_u}{\sigma'_{v0}} = 0.11 + 0.0037 \cdot PI$$ Ladd và Foott (1974) phát triển phương pháp SHANSEP nhằm chuẩn hóa việc đánh giá sức kháng cắt có xét đến hệ số quá cố kết ($OCR$). Hiroyuki Tanaka (2002) chỉ ra rằng chỉ số nén $C_c$ tương quan chặt chẽ với giới hạn chảy ($LL$), trong khi tính bất đẳng hướng của sức kháng cắt phụ thuộc vào cấu trúc vi mô và chỉ số dẻo $PI$.

Trong y văn địa kỹ thuật tồn tại hai luồng quan điểm tranh luận sâu sắc:

• Luồng quan điểm thứ nhất (Kinh nghiệm - Đơn giản hóa): Giả định tỷ số cố kết ngang trên cố kết đứng $m = C_h/C_v$ là một đại lượng ít biến động trong phạm vi hẹp ($2 \div 5$), cho phép nội suy trực tiếp trong thiết kế bấc thấm nhằm rút ngắn thời gian khảo sát.
• Luồng quan điểm thứ hai (Cấu trúc - Dị hướng thực nghiệm): Đại diện bởi Tavenas (1983), Seah và nnk (2004), Bergado và nnk (2002), khẳng định $C_h$ chịu sự kiểm soát khắt khe của vi cấu trúc trầm tích và mức độ xáo động khi cắm thiết bị thoát nước.

Tại Đông Nam Á, nghiên cứu của Bergado (2002) trên sét mềm Bangkok xác định $C_h = 2.8\text{ m}^2/\text{năm}$ khi giả thiết tỷ số đường kính vùng xáo động $d_s/d_w = 2$ và tỷ số thấm $k_h/k_s = 5$. Ngược lại, Seah và nnk (2004) tại cùng khu vực lại xác định $C_h = 0.75\text{ m}^2/\text{năm}$ với thông số $d_s/d_w = 2$ và $k_h/k_s = 1.4$. Tương tự, các nghiên cứu trên sét mềm Ariake (Nhật Bản) của Ohtsubo (1995) cho thấy khoáng vật Smectit chiếm ưu thế tạo nên độ nhạy rất cao, trong khi sét Singapore chiếm ưu thế bởi Kaolinit lại có độ nhạy trung bình.

Tại Việt Nam, các công trình của Đỗ Minh Toàn (1993, 2008), Nguyễn Thanh (1984), Nguyễn Văn Thơ (1999), Nguyễn Đình Thứ (2005, 2013), Nguyễn Thị Nụ (2015, 2016) đã làm sáng tỏ tính chất cơ lý của đất trầm tích ĐBBB ($mbQ_{2}^{1-2}hh$) và ĐBSCL ($amQ_{2}^{2-3}$). Phạm Thị Nghĩa (2005) chỉ ra rằng với bùn sét nguồn gốc biển - đầm lầy, tỷ số $C_h/C_v$ dao động mạnh từ $1.3 \div 7.6$. Tuy nhiên, tại dải ven biển Quảng Nam - Đà Nẵng, việc nghiên cứu còn manh mún, chưa có công trình nào đánh giá đồng bộ từ khoáng vật, hóa học, môi trường ăn mòn nước ngầm đến cơ chế cố kết dị hướng và phân vùng cấu trúc nền. Luận án đã định vị chính xác khoảng trống này, mở rộng tri thức địa kỹ thuật khu vực miền Trung và thiết lập cầu nối lý thuyết - thực nghiệm tương thích với các tiêu chuẩn quốc tế.

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

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

Luận án đã đóng góp vào nền tảng cơ học đất và địa kỹ thuật công trình qua việc mở rộng các lý thuyết cổ điển vào điều kiện trầm tích nhiệt đới ven biển miền Trung:

• Phát triển lý thuyết cố kết dị hướng địa phương: Kiểm chứng và mở rộng lời giải cố kết bán kính Barron - Hansbo bằng cách cung cấp hàm tương quan thực nghiệm giữa chỉ số nén $C_c$ với các chỉ tiêu vật lý cơ bản ($e_0, W, LL$) đặc thù cho các phân vị địa tầng $ambQ_{2}^{3}, mbQ_{2}^{2}, ambQ_{2}^{2}, abmQ_{2}^{1}$.
• Vận dụng chỉ số tự nén chặt Priklonsky ($k_d$): Luận giải bản chất cơ học của đất loại sét khu vực nghiên cứu đang trong giai đoạn tự nén chặt dở dang ($k_d = 0.13 \div 0.79$), giải thích căn nguyên tính nén lún cao và sức kháng cắt ban đầu thấp.
• Hoàn thiện mô hình trạng thái ứng suất hữu hiệu: Chứng minh sự suy giảm của tỷ số $S_u/\sigma'_{v0}$ theo chiều sâu đối với các thành tạo bùn sét hữu cơ và xác lập mối liên hệ giữa lực dính hữu hiệu $c'$, góc ma sát trong $\varphi'$ với hàm lượng hạt sét và chỉ số dẻo $PI$.

+─────────────────────────────────────────────────────────────────────────────+
|                    KHUNG LÝ THUYẾT & PHÂN TÍCH TỔNG HỢP                      |
+─────────────────────────────────────────────────────────────────────────────+
                                       ▲
  ┌────────────────────────────────────┼────────────────────────────────────┐
  │                                    │                                    │
┌─┴────────────────────────┐ ┌─────────┴────────────────┐ ┌─────────────────┴─┐
│     ĐỊA HỆ TN - KT       │ │    CƠ HỌC ĐẤT HIỆN ĐẠI   │ │  ĐỊA CHẤT TRẦM TÍCH │
│  (E.M. Sergeev/Bondarik) │ │ (Terzaghi/Barron/Casagr.)│ │   & KHOÁNG VẬT HỌC  │
│  Tương tác Nền - Công    │ │ Cố kết dị hướng (Ch, Cv),│ │  Nguồn gốc đa thành │
│  trình đường giao thông  │ │ sức kháng cắt (UU, CU)   │ │ tạo, Kaolinit-Illit │
└─┬────────────────────────┘ └─────────┬────────────────┘ └─────────────────┬─┘
  │                                    │                                    │
  └────────────────────────────────────┼────────────────────────────────────┘
                                       ▼
+─────────────────────────────────────────────────────────────────────────────+
|      BẢN ĐỒ CẤU TRÚC NỀN & MÔ HÌNH PHÂN TÍCH SỐ PLAXIS 8.5 / ARCMAP 10.2    |
+─────────────────────────────────────────────────────────────────────────────+

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

Khung phân tích của nghiên cứu tích hợp đa ngành dựa trên ba trụ cột lý thuyết:

1. Trụ cột Khoáng vật học - Thạch học trầm tích: Khảo sát thành phần khoáng vật sét (Kaolinit, Illit, Chlorit, Smectit), thành phần hóa học oxit ($SiO_2, Al_2O_3, Fe_2O_3$) và hàm lượng vật chất hữu cơ (MKN) để xác định nguồn gốc thành tạo.
2. Trụ cột Địa kỹ thuật thực nghiệm & Phân tích ngược: Kết hợp song song thí nghiệm nén không nở hông theo phương đứng ($C_v$) và phương ngang ($C_h$), thí nghiệm nén ba trục ($UU, CU$), thí nghiệm cắt cánh hiện trường ($VST$), thí nghiệm xuyên tiêu chuẩn ($SPT$) với phương pháp phân tích ngược quan trắc lún hiện trường của Asaoka (1978).
3. Trụ cột Địa tin học & Cơ học tính toán: Sử dụng hệ thông tin địa lý GIS (ArcMap 10.2) để tích hợp dữ liệu không gian, xây dựng bản đồ địa kỹ thuật và cấu trúc nền 1/50.000; sử dụng phần mềm phần tử hữu hạn Plaxis 8.5 mô phỏng ứng xử cố kết - biến dạng phi tuyến của nền đường đắp gia tải trước kết hợp PVD/giếng cát.

Điều kiện biên (Boundary conditions): Khung phân tích áp dụng cho đất dính mềm yếu bão hòa nước có chiều sâu phân bố $\le 30\text{ m}$, mực nước ngầm nông ($< 2\text{ m}$), dưới tác dụng của tải trọng đắp tĩnh và gia tải theo cấp giai đoạn trong thi công nền đường giao thông.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học Thực chứng (Positivism) kết hợp Hiện thực phê phán (Critical Realism), đi từ bản chất cấu trúc vi mô vật chất đến ứng xử vĩ mô công trình. Thiết kế nghiên cứu đa cấp độ (multi-level design) bao gồm:

• Cấp độ vĩ mô (Macro): Khảo sát điều kiện địa hình - địa mạo (6 kiểu địa hình chính), địa tầng kiến tạo và địa chất thủy văn toàn dải đồng bằng ven biển.
• Cấp độ trung mô (Meso): Khảo sát địa chất công trình qua hệ thống lỗ khoan sâu đến 30 m, thí nghiệm xuyên tiêu chuẩn SPT, cắt cánh VST, đo vẽ bản đồ chuyên đề tỷ lệ 1/50.000 trên nền tảng ArcMap 10.2.
• Cấp độ vi mô (Micro): Phân tích thành phần khoáng vật bằng nhiễu xạ tia X (XRD), phân tích hóa silicat, xác định hàm lượng hữu cơ bằng phương pháp nung (MKN), thí nghiệm cơ lý trong phòng chuyên sâu.

+─────────────────────────────────────────────────────────────────────────────+
|                     QUY TRÌNH NGHIÊN CỨU ĐỊA KỸ THUẬT RIGOROUS              |
+─────────────────────────────────────────────────────────────────────────────+
  [1] KHẢO SÁT & LẤY MẪU
      │── Ống mẫu thành mỏng Osterberg / Piston, bảo quản nghiêm ngặt
      └── Tiêu chí chất lượng Andresen (NGI): Độ biến dạng thể tích εv < 4%
  [2] THÍ NGHIỆM ĐA HỆ THỐNG
      │── Vi mô: XRD (Khoáng vật), Phân tích hóa học, Mất khi nung MKN (Hữu cơ)
      │── Cơ học trong phòng: Oedometer nén đứng (Cv) & ngang (Ch), Ba trục UU/CU
      └── Hiện trường: Vane Shear Test (Su-VST), Xuyên tiêu chuẩn (SPT N30)
  [3] TÍCH HỢP & PHÂN TÍCH NGƯỢC
      │── Tương quan thống kê: Cc = f(e0, W, LL); Cu = f(σ'v0, PI)
      │── Phân tích ngược Asaoka: Si = f(Si-1) từ mốc lún SP-01, SP-07
      └── Mô phỏng số FEM: Plaxis 8.5 (Soft Soil Model) & ArcMap 10.2 (GIS Map)
+─────────────────────────────────────────────────────────────────────────────+

Quy trình nghiên cứu rigorous

Nhằm triệt tiêu sai số do xáo động mẫu – vấn đề sống còn trong cơ học đất yếu, quy trình lấy mẫu và thí nghiệm tuân thủ chặt chẽ các hướng dẫn quốc tế:

• Kiểm soát tính nguyên trạng của mẫu: Mẫu được lấy bằng ống mẫu thành mỏng có van một chiều hoặc ống lấy mẫu piston. Độ nguyên trạng của mẫu được đánh giá theo tiêu chuẩn Viện Địa kỹ thuật Na Uy (NGI - Andresen và nnk), chỉ các mẫu có độ biến dạng thể tích $\varepsilon_v < 4%$ khi cố kết đẳng hướng trên máy ba trục ở áp lực hữu hiệu tự nhiên $\sigma'_{v0}$ mới được đưa vào tổ hợp số liệu tính toán.

"Dựa vào độ biến dạng khi nén cố kết hoặc độ biến dạng thể tích $\varepsilon_v$ khi cố kết đẳng hướng trên máy nén ba trục ở trạng thái ứng suất hữu hiệu tự nhiên, Andresen và nnk chỉ ra rằng các mẫu có độ biến dạng thể tích $\varepsilon_v < 4%$ được xem là đảm bảo tính nguyên trạng... (chất lượng mẫu: rất tốt $\varepsilon_v < 1%$; tốt $\varepsilon_v = 1 \div 2%$; khá tốt $\varepsilon_v = 2 \div 4%$; mẫu bị phá hoại $\varepsilon_v = 4 \div 8%$ và mẫu xấu $\varepsilon_v > 8%$)."

• Gia công và thí nghiệm cố kết ngang: Chế tạo mẫu theo phương nằm ngang trong phòng thí nghiệm để xác định trực tiếp hệ số thấm ngang $k_h$ và hệ số cố kết ngang $C_h(tp)$, so sánh trực tiếp với phương pháp Casagrande và Taylor trên mẫu đứng ($C_v$).
• Tam giác hóa phương pháp (Methodological Triangulation): Kiểm chứng chéo thông số $C_h$ thông qua 3 nguồn độc lập: (1) Thí nghiệm nén cố kết ngang trong phòng; (2) Phân tích ngược chuỗi số liệu quan trắc lún hiện trường theo phương pháp Asaoka tại các mốc quan trắc thực tế; (3) Mô phỏng số trên mô hình phần tử hữu hạn Plaxis 8.5.

Data và phân tích

Mẫu thí nghiệm và cơ sở dữ liệu địa kỹ thuật được tổng hợp từ hàng chục công trình thực tế, tiêu biểu là Dự án Đường cao tốc Đà Nẵng - Quảng Ngãi và Dự án nâng cấp mở rộng đường Nguyễn Tất Thành. Tổ hợp số liệu bao gồm hàng trăm mẫu đất nguyên trạng và 428 mẫu nước phân tích ăn mòn.

Phương pháp thống kê xác thực mô hình tương quan tuyến tính và phi tuyến được biểu diễn qua hệ số xác định $R^2$, cho phép đánh giá độ tin cậy của các công thức dự báo chỉ số nén $C_c$. Tại hiện trường, chuỗi số liệu quan trắc lún theo thời gian $S_i = f(S_{i-1})$ tại mốc SP-07 (Km12+480 cao tốc ĐN - QN) và mốc SP-01 (Km1+270 đường Nguyễn Tất Thành) được xử lý bằng thuật toán hồi quy sai phân bậc nhất để ngoại suy độ lún cố kết cuối cùng $S_\infty$ và hệ số cố kết ngang $C_h$.

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

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

1. Bản chất khoáng vật và hàm lượng hữu cơ đặc thù: Phân tích thạch học chỉ ra rằng bùn sét pha và bùn sét ven biển Quảng Nam - Đà Nẵng có thành phần khoáng vật sét chủ yếu là Kaolinit ($35 \div 50%$) và Illit ($20 \div 35%$), khoáng vật trương nở Smectit/Montmorillonit chỉ chiếm tỷ lệ nhỏ ($< 10%$). Điều này giải thích vì sao đất có độ dẻo và độ nhạy ở mức trung bình, thuận lợi hơn so với sét mềm ĐBSCL nhưng lại dễ suy giảm nhanh cường độ khi bão hòa nước. Hàm lượng vật chất hữu cơ (MKN) trong các phân vị $mbQ_{2}^{2}$ và $ambQ_{2}^{2}$ dao động từ $4.5%$ đến $14.2%$, phân bố không đồng đều theo chiều sâu.

+─────────────────────────────────────────────────────────────────────────────+
|             ĐẶC TÍNH ĂN MÒN NƯỚC DƯỚI ĐẤT (428 MẪU PHÂN TÍCH)               |
+─────────────────────────────────────────────────────────────────────────────+
  [█████████████████████████████████████░░░░░░░] 82.47% Ăn mòn rửa lũa (353/428)
  [████████████████████████████████░░░░░░░░░░░░] 79.67% Ăn mòn Carbonic (341/428)
  [░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░]  0.00% Ăn mòn Axit & Sulphat
+─────────────────────────────────────────────────────────────────────────────+

2. Môi trường địa hóa nước ngầm có tính xâm thực cao: Tổng hợp phân tích 428 mẫu nước dưới đất xác lập một phát hiện mang tính cảnh báo:

"Tính chất ăn mòn của nước dưới đất được xác định trên cơ sở kết quả phân tích TPHH của 10 mẫu nước dưới đất lấy từ các giếng sinh hoạt và tổng hợp các kết quả thí nghiệm của 418 mẫu nước lấy từ các hố khoan - đào của Liên đoàn Quy hoạch tài nguyên nước miền Trung... ăn mòn axit có 0/428 chiếm 0%; ăn mòn rửa lũa có 353/428 chiếm 82,47%; ăn mòn carbonic có 341/428 chiếm 79,67%; ăn mòn sulphat có 0/428 chiếm 0%."

Đặc tính này chứng minh giải pháp gia cố nền bằng chất kết dính vô cơ (xi măng, vôi) sẽ bị suy giảm hiệu quả nghiêm trọng nếu không có phụ gia kháng xâm thực carbonic và rửa lũa.

+──────────────────────────────────────────────────────────────────────────────+
|          QUAN HỆ DỊ HƯỚNG CỐ KẾT THẤM m = Ch/Cv THEO NGUỒN GỐC ĐỊA TẦNG      |
+──────────────────────────────────────────────────────────────────────────────+
|  Phân vị địa tầng  | Nguồn gốc thành tạo      | Tỷ số m = Ch/Cv thực nghiệm   |
+────────────────────+──────────────────────────+───────────────────────────────+
|  lbQ2(1-2)hh       | Hồ - đầm lầy             | 2.4 ÷ 3.5                     |
|  ambQ2(3)          | Sông - biển              | 1.5 ÷ 6.1                     |
|  mbQ2(2)           | Biển - đầm lầy           | 1.3 ÷ 7.6                     |
|  amQ2(2-3)         | Sông - biển (ĐBSCL)      | 1.35 ÷ 6.49 (Cố kết th.thường)|
+──────────────────────────────────────────────────────────────────────────────+
|  Kinh nghiệm cũ    | Quy trình 22TCN 262-2000 | 2.0 ÷ 5.0 (Cố định, ước tính) |
+──────────────────────────────────────────────────────────────────────────────+

3. Tính bất đẳng hướng cố kết và sự phụ thuộc áp lực: Kết quả thí nghiệm khẳng định tỷ số $m = C_h/C_v$ của đất yếu khu vực dao động rất rộng ($1.3 \div 7.6$) và giảm dần khi cấp áp lực nén gia tăng vượt qua áp lực tiền cố kết $P_c$. Việc áp dụng cứng nhắc hệ số $m = 2 \div 5$ dẫn đến sai lệch lớn trong bố trí cự ly bấc thấm $L$.

4. Xác lập định nghĩa và tiêu chí nhận diện đất yếu cục bộ: Trên cơ sở lý luận và thực nghiệm tương tác nền - công trình, luận án chuẩn hóa định nghĩa:

"Đất yếu là đất có khả năng chịu tải thấp (<50 - 100 kPa), độ bền kháng cắt được đặc trưng bằng cường độ lực dính đơn vị và góc ma sát trong thấp (c<10 kPa, φ<80), khả năng biến dạng được đặc trưng bằng hệ số nén lún lớn (a≥ 10 kPa-1) và môđun tổng biến dạng nhỏ (E0≤5000 kPa). Nếu không có biện pháp xử lý thích hợp thì không thể xây dựng công trình trên đó."

5. Phân loại 5 nhóm cấu trúc nền (CTN) và hệ thống bản đồ 1/50.000: Thành lập bản đồ phân bố đất yếu và bản đồ cấu trúc nền tỷ lệ 1/50.000, phân loại địa tầng thành 5 nhóm đất đá: Đá cứng, đá nửa cứng, đất rời, đất dính và đất yếu, làm cơ sở phân vùng xử lý kỹ thuật trực quan cho các nhà quản lý và tư vấn thiết kế.

+─────────────────────────────────────────────────────────────────────────────+
|               BẢN ĐỒ ĐỊA KỸ THUẬT & 5 NHÓM ĐẤT ĐÁ KHU VỰC                   |
+─────────────────────────────────────────────────────────────────────────────+
  [Nhóm 1] ĐÁ CỨNG: Granit Hải Vân, Núi Vú, Amphibolit (Rn = 1000 - 2000 kPa)
  [Nhóm 2] ĐÁ NỬA CỨNG: Cát kết, bột kết Ái Nghĩa, Tân Lâm (Rn = 300 - 700 kPa)
  [Nhóm 3] ĐẤT RỜI: Cát thô, vừa, sỏi cuội Pleistocen/Holocen (N30 = 10 - 50)
  [Nhóm 4] ĐẤT DÍNH: Sét pha, sét dẻo cứng nửa cứng (IL < 0.75, c = 15 - 35 kPa)
  [Nhóm 5] ĐẤT YẾU: Bùn sét, bùn sét pha (ambQ2, mbQ2) (e0 > 1, c < 10, Su < 35)
+─────────────────────────────────────────────────────────────────────────────+

Implications đa chiều

• Về mặt lý thuyết: Cung cấp hệ cơ sở dữ liệu địa kỹ thuật hoàn chỉnh đầu tiên cho dải đồng bằng ven biển miền Trung, lấp đầy khoảng trống dữ liệu giữa hai đầu châu thổ Bắc Bộ và Nam Bộ.
• Về phương pháp luận: Chuẩn hóa quy trình tam giác hóa xác định $C_h$ thông qua phân tích ngược quan trắc lún hiện trường kết hợp thí nghiệm trong phòng, loại bỏ tính tùy tiện trong nội suy địa kỹ thuật.
• Về ứng dụng thực tiễn: Tối ưu hóa sơ đồ bố trí thiết bị tiêu thoát nước thẳng đứng (PVD, Giếng cát). Giảm cự ly bố trí bấc thấm ở các tầng trầm tích $mbQ_{2}^{2}$ có $m$ thấp và giãn khoảng cách ở các tầng $ambQ_{2}^{3}$ có $m$ cao, giúp tiết kiệm $12 \div 18%$ chi phí xử lý nền đắp và đảm bảo kiểm soát độ lún dư theo tiêu chuẩn đường cao tốc.
• Về chính sách và tiêu chuẩn: Cung cấp luận cứ khoa học thực chứng để các cơ quan quản lý giao thông (Bộ GTVT, Sở GTVT Quảng Nam và Đà Nẵng) điều chỉnh chỉ dẫn kỹ thuật khảo sát, hạn chế rủi ro phá hoại nền đường ven biển.

Limitations và Future Research

Mặc dù đạt được những kết quả toàn diện, nghiên cứu vẫn tồn tại một số hạn chế khách quan:

1. Phạm vi độ sâu khảo sát: Chiều sâu nghiên cứu khống chế ở mức $\le 30\text{ m}$. Đối với các công trình cầu vượt, hầm chui có tải trọng cục bộ lớn truyền xuống tầng đá gốc sâu hơn, dữ liệu của luận án chỉ đóng vai trò tham khảo tầng mặt.
2. Ảnh hưởng thủy triều và động lực ven biển: Chưa tích hợp đầy đủ mô hình tương tác thủy lực động (coupled hydro-mechanical modeling) xét đến sự dao động mực nước ngầm theo chu kỳ triều cường và hiện tượng dâng cao do bão.
3. Thí nghiệm từ biến lâu dài: Các thí nghiệm nén cố kết thứ cấp ($C_\alpha$) trên đất than bùn và bùn sét hữu cơ chưa kéo dài đủ thời gian để xây dựng mô hình dự báo biến dạng từ biến trong vòng đời 50 - 100 năm của công trình.

Chương trình nghiên cứu tương lai (Future Research Agenda):

• Phát triển mô hình cấu trúc 3D liên kết giữa biến dạng nén lún và dòng thấm không ổn định dưới tác động của biến đổi khí hậu và nước biển dâng.
• Nghiên cứu cơ chế suy thoái liên kết xi măng - đất dưới tác động ăn mòn rửa lũa và carbonic trong môi trường nước dưới đất ven biển.
• Đánh giá khả năng ứng dụng công nghệ hút chân không kết hợp bấc thấm không gia tải (Vacuum Preloading) cho các đô thị nén ven biển nhằm giảm thiểu phát thải carbon từ vật liệu đắp.

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

+─────────────────────────────────────────────────────────────────────────────+
|                          MA TRẬN TÁC ĐỘNG ĐA LĨNH VỰC                       |
+─────────────────────────────────────────────────────────────────────────────+
| Lĩnh vực      | Cơ chế tác động              | Kết quả định lượng kỳ vọng   |
+---------------+------------------------------+------------------------------+
| Học thuật     | Cung cấp baseline data cho   | Trích dẫn chuẩn tắc trong các|
|               | địa kỹ thuật miền Trung      | luận án, bài báo ISI/Scopus  |
+---------------+------------------------------+------------------------------+
| Công nghiệp   | Tối ưu hóa cự ly PVD,        | Tiết kiệm 12 - 18% chi phí xử|
| Xây dựng      | rút ngắn thời gian chờ lún   | lý nền; triệt tiêu sự cố trượt|
+---------------+------------------------------+------------------------------+
| Quản lý nhà   | Chuẩn hóa chỉ dẫn kỹ thuật   | Ban hành cẩm nang địa kỹ thuật|
| nước & QH     | và bản đồ cấu trúc nền       | xây dựng hạ tầng ven biển    |
+---------------+------------------------------+------------------------------+
| Xã hội &      | Đảm bảo an toàn giao thông,  | Giảm chi phí duy tu sửa chữa,|
| Môi trường    | giảm khai thác mỏ vật liệu   | bảo tồn cảnh quan cồn cát    |
+─────────────────────────────────────────────────────────────────────────────+

Công trình mang lại tác động học thuật sâu sắc khi thiết lập hệ quy chiếu cơ lý đất yếu cho miền Trung, mở ra hướng nghiên cứu so sánh liên vùng với ĐBBB và ĐBSCL. Trong thực tiễn công nghiệp xây dựng, kết quả nghiên cứu đã được ứng dụng trực tiếp để giải quyết bài toán lún phức tạp tại các dự án trọng điểm quốc gia (Cao tốc Đà Nẵng - Quảng Ngãi, tuyến đường gom ven biển Đà Nẵng - Hội An). Về mặt xã hội, việc nâng cao độ tin cậy của thiết kế nền đường giúp giảm thiểu chi phí bảo trì định kỳ, ngăn ngừa tắc nghẽn giao thông do sự cố mặt đường và bảo vệ môi trường địa chất ven bờ.

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

• Nghiên cứu sinh và Giảng viên ngành Địa kỹ thuật - Địa chất công trình: Thụ hưởng hệ thống dữ liệu thực nghiệm toàn diện, phương pháp phân tích ngược Asaoka và mô hình phân vùng cấu trúc nền để phát triển các đề tài chuyên sâu.
• Các kỹ sư và Công ty tư vấn thiết kế giao thông (TEDI, CIENCO...): Tiếp cận bộ thông số cơ lý đầu vào ($C_c, C_v, C_h, m, S_u$) đáng tin cậy, chấm dứt việc tra bảng kinh nghiệm thiếu chuẩn xác, nâng cao tính cạnh tranh của hồ sơ thiết kế kỹ thuật.
• Nhà thầu thi công hạ tầng: Có căn cứ khoa học để lập biện pháp thi công đắp gia tải theo giai đoạn, kiểm soát tốc độ đắp và thời gian chờ cố kết, tránh nguy cơ sạt lở nền đường đang thi công.
• Sở Giao thông Vận tải và Sở Xây dựng các tỉnh miền Trung: Sở hữu công cụ bản đồ GIS 1/50.000 phục vụ quy hoạch hướng tuyến giao thông, tối ưu hóa tổng mức đầu tư hạ tầng công cộng.

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

+─────────────────────────────────────────────────────────────────────────────+
|                     GIẢI ĐÁP CÁC CÂU HỎI KHOA HỌC THEN CHỐT                 |
+─────────────────────────────────────────────────────────────────────────────+

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?

Luận án đã mở rộng lý thuyết cố kết bán kính của Barron (1948) và Hansbo (1981) thông qua việc lượng hóa tường minh tính bất đẳng hướng thấm ngang - đứng ($m = C_h/C_v$) gắn liền với điều kiện nguồn gốc trầm tích Đệ Tứ ven biển miền Trung. Nghiên cứu đã bác bỏ việc áp dụng giá trị cố định $m = 2 \div 5$ của tiêu chuẩn ngành, đồng thời tích hợp chỉ số tự nén chặt Priklonsky ($k_d$) để giải thích trạng thái cơ học dở dang của trầm tích trẻ.

2. Đột phá phương pháp luận thể hiện qua việc so sánh với các nghiên cứu trước đây như thế nào?

So với các nghiên cứu trước đây vốn chỉ dựa vào thí nghiệm nén một trục tiêu chuẩn trong phòng (dễ mất tính nguyên trạng) hoặc giả định tỷ số $C_h/C_v$, luận án đã thiết lập quy trình kiểm chứng chéo 3 cấp độ: Thí nghiệm hộp nén ngang chuyên dụng $\leftrightarrow$ Thí nghiệm hiện trường CPTu/VST $\leftrightarrow$ Phân tích ngược biến dạng lún thực tế Asaoka ($S_i = f(S_{i-1})$) trên phần mềm Plaxis 8.5. Phương pháp này vượt trội về độ tin cậy so với nghiên cứu đơn lẻ của Bergado (2002) hay Seah (2004) tại Bangkok.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là môi trường địa hóa nước dưới đất vùng nghiên cứu có mức độ ăn mòn rửa lũa (chiếm $82.47%$) và ăn mòn carbonic (chiếm $79.67%$) cực kỳ phổ biến, trong khi ăn mòn axit và sulphat hoàn toàn bằng $0%$. Điều này cảnh báo nguy cơ suy giảm độ bền nghiêm trọng của các giải pháp gia cố nền bằng cọc đất - xi măng thông thường, đòi hỏi bắt buộc phải sử dụng chất kết dính bền sunfat/carbonic hoặc kết hợp bấc thấm tiêu thoát nước truyền thống.

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

Luận án cung cấp chi tiết toàn bộ quy trình khảo sát lấy mẫu nguyên trạng bằng ống thành mỏng đạt chuẩn NGI ($\varepsilon_v < 4%$), phương pháp gia công mẫu nén ngang, thuật toán xử lý chuỗi số liệu quan trắc lún theo bước thời gian đều nhau $\Delta t$ để lập đồ thị hồi quy Asaoka, và quy trình phân nhóm cấu trúc nền trên ArcMap 10.2, đảm bảo khả năng tái lập và áp dụng cho các vùng ven biển khác.

5. Định hướng 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 tới định hướng: (1) Ứng dụng trí tuệ nhân tạo (Machine Learning) trong dự báo trường lún không gian từ dữ liệu CPTu lớn; (2) Hoàn thiện mô phỏng tương tác đất - kết cấu dưới tải trọng động giao thông nặng; (3) Phát triển các vật liệu địa kỹ thuật sinh học (bio-geotechnics) thân thiện với môi trường để gia cố đất bùn hữu cơ ven biển.

Kết luận

+─────────────────────────────────────────────────────────────────────────────+
|               TỔNG KẾT CÁC ĐÓNG GÓP CỐT LÕI CỦA LUẬN ÁN                    |
+─────────────────────────────────────────────────────────────────────────────+
  (1) Hệ thống hóa điều kiện ĐKT & phân loại 5 nhóm đất đá ven biển QN-ĐN
  (2) Xác lập bản chất khoáng vật (Kaolinit-Illit) & hóa học (Ăn mòn Carbonic)
  (3) Chuẩn hóa thông số cố kết dị hướng Ch, Cv & hàm tương quan Cc = f(e0,W,LL)
  (4) Đột phá phương pháp phân tích ngược Asaoka xác định Ch từ mốc quan trắc
  (5) Xây dựng bản đồ cấu trúc nền 1/50.000 phục vụ tối ưu hóa bấc thấm/giếng cát
  (6) Cung cấp luận cứ khoa học thực chứng phục vụ tiêu chuẩn thiết kế quốc gia
+─────────────────────────────────────────────────────────────────────────────+

Luận án tiến sĩ kỹ thuật của NCS. Nguyễn Thị Ngọc Yến đã hoàn thành xuất sắc mục tiêu nghiên cứu với 6 đóng góp mang tính bước ngoặt:

1. Hệ thống hóa toàn diện điều kiện địa kỹ thuật: Đánh giá đồng bộ từ địa hình - địa mạo, địa tầng trầm tích đến địa chất thủy văn dải ven biển Quảng Nam - Đà Nẵng, phân chia đất đá thành 5 nhóm cơ lý rõ ràng.
2. Làm sáng tỏ bản chất thành phần vật chất: Xác định ưu thế khoáng vật sét Kaolinit ($35 \div 50%$) và Illit ($20 \div 35%$), hàm lượng hữu cơ ($4.5 \div 14.2%$) và cảnh báo nguy cơ ăn mòn carbonic ($79.67%$) của nước ngầm đối với kết cấu móng/gia cố.
3. Chuẩn hóa thông số cố kết và kháng cắt: Cung cấp bảng giá trị đặc trưng cơ lý chuẩn tắc ($C_v, C_h, m, S_u, c', \varphi'$) và các phương trình hồi quy dự báo $C_c$, chấm dứt tình trạng thiết kế dựa trên số liệu tra bảng thiếu căn cứ.
4. Đột phá trong phương pháp phân tích ngược: Chứng minh tính ưu việt của phương pháp Asaoka trong việc xác định chính xác hệ số cố kết ngang $C_h$ và độ lún cuối cùng từ dữ liệu quan trắc hiện trường thực tế.
5. Thành lập hệ thống bản đồ chuyên đề 1/50.000: Xây dựng bản đồ phân bố đất yếu và bản đồ cấu trúc nền bằng công nghệ GIS, hỗ trợ trực quan hóa công tác lựa chọn giải pháp xử lý kỹ thuật (đắp gia tải, PVD, giếng cát, đệm cát).
6. Định hình bước chuyển đổi mô hình (Paradigm shift): Chuyển dịch tư duy thiết kế nền đường trên đất yếu từ phương pháp kinh nghiệm đơn thuần sang phương pháp tiếp cận hệ thống - địa hệ tự nhiên kỹ thuật, tạo di sản khoa học thực chứng bền vững cho ngành xây dựng công trình giao thông Việt Nam.