Tổng quan về luận án

Sự phát triển nhanh chóng của các đô thị ven biển và mạng lưới hạ tầng giao thông trọng điểm tại dải đồng bằng ven biển (ĐBVB) Quảng Nam - Đà Nẵng đặt ra những thách thức kỹ thuật phức tạp đối với công tác nền móng công trình. Khu vực này có cấu trúc địa chất Đệ Tứ đa nguồn gốc, nơi các thành tạo đất yếu trầm tích trẻ tuổi Holocen và Pleistocen muộn phân bố không đồng đều, bị che phủ cục bộ và có sự biến thiên lớn về bề dày.

Khoảng trống nghiên cứu (research gap) lớn nhất trong giai đoạn qua là tình trạng thiếu vắng các công trình nghiên cứu cơ bản, đồng bộ và có hệ thống về mối quan hệ hữu cơ giữa thành phần vật chất (TPVC), điều kiện thành tạo trầm tích và tính chất cơ lý (TCCL) của đất yếu duyên hải miền Trung. Các đơn vị tư vấn thiết kế công trình giao thông phần lớn phải ngoại suy dữ liệu từ Đồng bằng Bắc Bộ (ĐBBB) hoặc Đồng bằng sông Cửu Long (ĐBSCL), hoặc áp dụng cơ học các chỉ tiêu tra bảng định tính theo quy trình 22TCN 262-2000 mà thiếu cơ sở thực nghiệm địa phương. Điều này dẫn đến hàng loạt sự cố biến dạng lún không đều, trượt trôi mái taluy và lún dư vượt giới hạn cho phép sau khi đưa vào khai thác tại các tuyến giao thông trọng điểm như đường vành đai Đà Nẵng, đường ven sông Tuyên Sơn - Túy Loan.

Luận án tiến sĩ kỹ thuật chuyên ngành Địa kỹ thuật Xây dựng (Mã số: 9580211) của tác giả Nguyễn Thị Ngọc Yến, dưới sự hướng dẫn của PGS.TS. Bùi Văn Trường và NGND.GS.TSKH. Nguyễn Thanh tại Trường Đại học Thủy lợi (2018), được thực hiện nhằm giải quyết trọn vẹn bài toán khoa học và thực tiễn này. Nghiên cứu tập trung vào 4 câu hỏi nghiên cứu (Research Questions) và 4 giả thuyết khoa học (Hypotheses) cốt lõi:

  • RQ1: Điều kiện địa chất và quy luật phân bố không gian của các thành tạo trầm tích đất yếu tại dải ĐBVB Quảng Nam - Đà Nẵng có những đặc trưng gì khác biệt so với các đồng bằng châu thổ lớn?
    • H1: Đất yếu khu vực nghiên cứu mang tính đa nguồn gốc (sông, biển, đầm lầy, vũng vịnh), phân bố dạng thấu kính không liên tục, bị chi phối mạnh mẽ bởi cấu trúc móng sụt lún bậc thang và mạng lưới đứt gãy kiến tạo trẻ.
  • RQ2: Mối quan hệ tương hỗ giữa thành phần khoáng vật, hàm lượng vật chất hữu cơ, thành phần hóa học với các chỉ tiêu cơ lý biến dạng và chống cắt của đất yếu diễn ra theo quy luật nào?
    • H2: Thành phần khoáng vật sét ưu thế Kaolinit và Illit kết hợp với hàm lượng hữu cơ biến động ($2 - 10%$) tạo nên trạng thái tự nén chặt vừa ($K_g = 0{,}13 - 0{,}79$), quyết định tính nén lún cao và tính bất đẳng hướng thấm rõ rệt.
  • RQ3: Tỷ số cố kết thấm ngang trên đứng ($m = C_h / C_v$) trong phòng và ngoài hiện trường sai khác như thế nào so với các khung quy phạm thiết kế hiện hành?
    • H3: Tỷ số $m = C_h / C_v$ thực tế nằm trong khoảng hẹp hơn ($1{,}5 - 3{,}5$) so với giả định định tính ($2 - 5$) của quy phạm 22TCN 262-2000, đòi hỏi phương pháp phân tích ngược từ quan trắc hiện trường để kiểm soát chính xác tốc độ cố kết.
  • RQ4: Mô hình cấu trúc nền nào phản ánh đúng nhất cơ chế làm việc của nền đường đắp và tối ưu hóa giải pháp xử lý kỹ thuật?
    • H4: Việc tích hợp mô hình phần tử hữu hạn với cơ sở dữ liệu cấu trúc nền theo chiều sâu cho phép dự báo chính xác độ lún theo thời gian, chiều cao đắp bù lún và nâng cao hệ số ổn định tổng thể mái taluy.

Về khung lý thuyết, công trình tích hợp Lý thuyết Hệ thống Địa chất Công trình - Địa kỹ thuật (Địa hệ Tự nhiên - Kỹ thuật) của E.M. Bondarik (1981) và Schellnhuber & Wenzel (1998) với Nguyên lý ứng suất hữu hiệu và Lý thuyết cố kết thấm của Karl Terzaghi (1925), bổ sung lý thuyết cố kết không gian của N. Carrillo (1942) và Barron (1948).

Đóng góp đột phá của luận án được định lượng cụ thể: Xây dựng bản đồ phân bố đất yếu và bản đồ cấu trúc nền tỷ lệ 1/50.000 cho toàn dải ven biển từ quận Liên Chiểu (Đà Nẵng) đến huyện Núi Thành (Quảng Nam); thiết lập hệ thống phương trình tương quan thực nghiệm chuẩn xác giữa các chỉ tiêu biến dạng - sức kháng cắt; giảm sai số dự báo lún cố kết từ $25 - 40%$ xuống dưới $10%$; tối ưu hóa khoảng cách bấc thấm (PVD) và giếng cát, giúp tiết kiệm chi phí xử lý nền móng và ngăn ngừa sự cố công trình. Phạm vi nghiên cứu bao quát chiều sâu khống chế địa tầng đến 30 m, xử lý dữ liệu từ 428 mẫu phân tích hóa học nước ngầm, 116 điểm quan trắc thủy văn và hàng trăm tổ hợp mẫu cơ lý thí nghiệm chuyên sâu.

Literature Review và Positioning

Lịch sử 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 trường phái lý thuyết kinh điển. Nền móng lý thuyết cố kết một chiều được Karl Terzaghi thiết lập từ năm 1925, dựa trên nguyên lý áp lực hữu hiệu $\sigma' = \sigma - u_w$ (trong đó $\sigma$ là ứng suất tổng, $\sigma'$ là ứng suất hữu hiệu và $u_w$ là áp lực nước lỗ rỗng). Lý thuyết này sau đó được Arthur Casagrande (1936, 1938) hoàn thiện qua việc định nghĩa áp lực tiền cố kết ($P_c$) và phương pháp logarit thời gian ($\log t$) để xác định hệ số cố kết thẳng đứng ($C_v$), cùng với phương pháp căn bậc hai thời gian ($\sqrt{t}$) của D.W. Taylor (1942, 1948).

Trong lĩnh vực cố kết đa chiều và thoát nước đứng, N. Carrillo (1942) đã đưa ra định lý phân tách phương trình vi phân cố kết ba chiều thành tổ hợp nghiệm cố kết thẳng đứng và hướng tâm, mở đường cho mô hình tính toán giếng cát và bấc thấm của Barron (1948). Vấn đề xác định hệ số cố kết ngang ($C_h$) và độ bất đẳng hướng thấm ($k_h / k_v$) tiếp tục được Rowe, Tavenas (1983) và Bergado et al. (2002) phát triển thông qua các thiết bị thấm ngang chuyên dụng và phương pháp phân tích ngược quan trắc hiện trường.

Về sức kháng cắt, A.W. Skempton (1957) thiết lập tương quan 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 thẳng đứng ($c_u / \sigma'_v$) với chỉ số dẻo ($PI$) theo biểu thức thực nghiệm $c_u / \sigma'_v = 0{,}11 + 0{,}0037 \cdot PI$. Tiếp đó, L. Bjerrum (1972, 1973), Azzouz (1983), Kulhawy & Mayne (1990) và phương pháp SHANSEP của Ladd & Foott (1974) đã làm sáng tỏ tính bất đẳng hướng của sức kháng cắt và ảnh hưởng của lịch sử ứng suất ($OCR$).

Trong bức tranh học thuật quốc tế, hai cuộc tranh luận lý thuyết lớn vẫn tồn tại:

  1. Bản chất của ứng suất tiền cố kết ($P_c$) và tỷ số quá cố kết ($OCR$): Trường phái cơ học cổ điển thuần túy (Casagrande, 1936) xem $P_c$ là ứng suất lớn nhất mà tầng đất từng chịu trong quá khứ do tải trọng địa tầng. Ngược lại, trường phái địa kỹ thuật hiện đại (Bjerrum, 1967; Janbu, 1977; Hanzawa, 1989) chứng minh rằng đất yếu có $OCR > 1$ có thể hình thành do quá trình nén thứ cấp lâu dài, hiện tượng hóa già (aging), biến đổi địa hóa và phản ứng xi măng hóa tự nhiên giữa các khoáng vật sét mà không nhất thiết phải trải qua quá trình bóc mòn dỡ tải.
  2. Xác định hệ số cố kết ngang ($C_h$) và vùng xáo động (smear zone): Các mô hình lý thuyết của Barron hay Hansbo giả định tỷ số $C_h / C_v$ là hằng số vật liệu. Tuy nhiên, các nhà nghiên cứu thực nghiệm (Bergado et al., Seah et al.) chỉ ra rằng tỷ số này phụ thuộc nghiêm trọng vào mức độ phá hủy kết cấu tự nhiên khi cắm thiết bị thoát nước, tỷ lệ đường kính vùng xáo động ($d_s / d_w$) và tỷ số suy giảm hệ số thấm ($k_h / k_s$), khiến các giá trị đo trong phòng thường sai lệch lớn so với vận hành thực tế.

Tại Việt Nam, các nghiên cứu của Đỗ Minh Toàn (1993), Nguyễn Văn Thơ (1999), Nguyễn Viết Tình (2001), Bùi Đức Hải (2006) và Nguyễn Thị Nụ (2015) đã làm sáng tỏ đặc tính đất yếu tại ĐBBB và ĐBSCL. Tuy nhiên, khu vực duyên hải miền Trung hoàn toàn bị bỏ ngỏ về mặt hệ thống dữ liệu.

So sánh với các nghiên cứu quốc tế điển hình:

  • Trường hợp đất sét mềm Bangkok (Thái Lan): Nghiên cứu của Bergado et al. (2002) tại độ sâu $4 - 8\text{ m}$ xác định $C_h = 2{,}8\text{ m}^2/\text{năm}$ với giả thiết $d_s / d_w = 2$ và $k_h / k_s = 5$; đất có hàm lượng khoáng Montmorillonite/Smectite cao, độ nhạy lớn, hàm lượng vôi tối ưu để gia cố là $5 - 10%$ (Balasubramaniam, 1988).
  • Trường hợp đất sét biển Singapore: Nghiên cứu của Seah et al. (2004) xác định $C_h = 0{,}75\text{ m}^2/\text{năm}$ với $d_s / d_w = 2$ và $k_h / k_s = 1{,}4$; đất có thành phần khoáng vật sét ưu thế là Kaolinite với độ nhạy từ thấp đến trung bình (Ohtsubo et al., Tanaka, 2002).

Công trình của Nguyễn Thị Ngọc Yến định vị chính xác khoảng trống khoa học này, chứng minh đất yếu ĐBVB Quảng Nam - Đà Nẵng có thành phần khoáng vật tương đồng với đất sét Singapore (giàu Kaolinit và Illit) nhưng mang tính phân tầng phức tạp của trầm tích sông - biển hỗn hợp miền Trung, thiết lập hệ thông số $C_h, C_v, m$ chuẩn xác cho vùng nghiên cứu.

Đó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à làm sâu sắc thêm Lý thuyết Hệ thống Địa chất Công trình - Địa kỹ thuật của E.M. Bondarik (1981) và GS.TSKH. Nguyễn Thanh (1984) thông qua việc xây dựng khái niệm "Địa hệ Tự nhiên - Kỹ thuật" (TN - KT) áp dụng cho công trình giao thông ven biển. Nghiên cứu đã chứng minh rằng tính chất cơ lý xây dựng của đất yếu không phải là các đại lượng tĩnh cô lập, mà là kết quả của sự tương tác động giữa 4 phân hệ: Thạch quyển (thành phần khoáng vật, thạch học, cấu trúc móng), Thủy quyển (chế độ thủy văn, mực nước ngầm, tính ăn mòn hóa học), Sinh - Khí quyển (quá trình tích tụ vật chất hữu cơ, biến đổi khí hậu) và Hệ thống Kỹ thuật (tải trọng nền đắp, biện pháp tiêu thoát nước đứng).

Đồng thời, luận án đóng góp vào Lý thuyết Cố kết và Cơ học Đất không bão hòa / bão hòa qua việc bổ sung định nghĩa chuẩn xác về đất yếu phục vụ công trình giao thông:

"Đất yếu là đất có khả năng chịu tải thấp ($< 50\text{ kPa}$), độ bền kháng cắt thấp ($c < 10\text{ kPa}, \varphi < 8^\circ$), khả năng biến dạng được đặc trưng bằng hệ số nén lún lớn ($a > 10^{-3}\text{ kPa}^{-1}$) và môđun biến dạng nhỏ ($E_0 < 5000\text{ kPa}$), không phù hợp làm nền trực tiếp cho công trình."

Luận án cũng xác nhận luận điểm của Hanzawa (1989) và Bjerrum (1973) trên đối tượng bùn sét duyên hải miền Trung: Hiện tượng tiền cố kết biểu kiến ($OCR > 1$) tại các lớp bùn sét pha $mbQ_2^3$ và $ambQ_2^2$ chủ yếu bắt nguồn từ quá trình nén thứ cấp lâu dài và liên kết xi măng hóa tự nhiên chứ không phải do sự thay đổi tải trọng địa tầng đơn thuần.

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

Khung phân tích của luận án tích hợp liên ngành ba trụ cột lý thuyết:

  1. Lý thuyết Cấu trúc nền và Môi trường Địa chất: Phân tích nguồn gốc trầm tích học kết hợp tướng đá - cổ địa lý.
  2. Lý thuyết Cố kết Thấm Đa chiều và Trạng thái Tới hạn: Ứng dụng các sơ đồ thí nghiệm nén một trục, ba trục ($UU, CU$) và mô hình thoát nước hướng tâm của Barron - Hansbo.
  3. Phương pháp Phân tích Ngược (Back-Analysis): Kết hợp giải tích lý thuyết, thuật toán đồ thị Asaoka và mô hình hóa số phần tử hữu hạn (FEM).
                      ┌────────────────────────────────────────┐
                      │    ĐỊA HỆ TỰ NHIÊN - KỸ THUẬT (TN-KT)  │
                      │  (Môi trường Địa chất - Nền đường đắp) │
                      └──────────────────┬─────────────────────┘
                                         │
        ┌────────────────────────────────┼────────────────────────────────┐
        ▼                                ▼                                ▼
┌──────────────────────┐      ┌──────────────────────┐      ┌──────────────────────┐
│  THÀNH PHẦN VẬT CHẤT │      │ CƠ HỌC BIẾN DẠNG     │      │ CẤU TRÚC NỀN & PHÂN  │
│  & ĐIỀU KIỆN ĐỊA CHẤT│      │ & SỨC KHÁNG CẮT      │      │ TÍCH NGƯỢC HIỆN TRƯỜNG│
├──────────────────────┤      ├──────────────────────┤      ├──────────────────────┤
│• TPKV (Kaolinit,     │      │• Thí nghiệm Oedometer│      │• Bản đồ ĐKT 1/50.000 │
│  Illit, Chlorit)     │      │  (Casagrande, Taylor,│      │  (ArcMap 10.2)       │
│• TPHH & Hữu cơ       │      │  JGS-2000 3t)        │      │• Quan trắc lún SP-01,│
│• Thủy văn & Ăn mòn   │      │• Nén 3 trục (UU, CU) │      │  SP-07 (Asaoka)      │
│  (428 mẫu nước)      │      │• Cắt cánh VST, SPT   │      │• Mô phỏng số FEM     │
│• Tự nén chặt Kg      │      │• Tương quan Ch, Cv, m│      │  (Plaxis 8.5)        │
└──────────┬───────────┘      └──────────┬───────────┘      └──────────┬───────────┘
           │                             │                             │
           └─────────────────────────────┼─────────────────────────────┘
                                         ▼
                      ┌────────────────────────────────────────┐
                      │    TỐI ƯU HÓA THIẾT KẾ XỬ LÝ NỀN       │
                      │    (PVD, Giếng cát, Đắp bù lún)        │
                      └────────────────────────────────────────┘

Điều kiện biên (Boundary conditions) của khung phân tích: Áp dụng cho các thành tạo trầm tích Đệ Tứ ven biển miền Trung có chiều sâu phân bố $\le 30\text{ m}$, mực nước ngầm nằm nông ($0 - 2\text{ m}$), chịu tác động của tải trọng nền đường đắp có chiều cao từ $2 - 6\text{ m}$.

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 lập trường triết học thực chứng (Positivism) kết hợp hiện thực phản tư (Critical Realism), đi từ quan sát thực nghiệm vi mô (khoáng vật, cấu trúc lỗ rỗng) đến quy luật cơ học vĩ mô và mô hình hóa toán học. Thiết kế đa cấp độ (Multi-level design) được triển khai nhất quán:

  • Cấp độ Vĩ mô: Phân vùng địa kỹ thuật và cấu trúc nền toàn dải ĐBVB Quảng Nam - Đà Nẵng (dài hơn 100 km từ Núi Thành đến Liên Chiểu) trên nền tảng GIS.
  • Cấp độ Trung mô: Thiết lập các mặt cắt địa kỹ thuật đặc trưng dọc theo các công trình tuyến thực tế (đường Nguyễn Tất Thành nối dài, đường cao tốc Đà Nẵng - Quảng Ngãi).
  • Cấp độ Vi mô: Phân tích thành phần khoáng vật bằng nhiễu xạ tia X (XRD), phân tích nhiệt vi sai (DTA), xác định hàm lượng vật chất hữu cơ bằng phương pháp sấy nung (MKN) và hóa học nước dưới đất.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý mẫu tuân thủ nghiêm ngặt các tiêu chuẩn kiểm soát chất lượng quốc tế. Độ nguyên trạng của mẫu đất dính mềm yếu trong phòng thí nghiệm được kiểm soát theo tiêu chí biến dạng thể tích của Andresen et al.:

"Mẫu có độ biến dạng thể tích $\varepsilon_v < 4%$ 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 $\sigma'_0$ đượ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 - 2%$; khá tốt $\varepsilon_v = 2 - 4%$; mẫu bị xáo động phá hoại khi $\varepsilon_v > 4%$)."

Quy trình thí nghiệm và công cụ phân tích bao gồm:

  • Thí nghiệm nén cố kết không nở hông (Oedometer): Xác định hệ số nén lún $a$, chỉ số nén $C_c$, chỉ số nở $C_s$, áp lực tiền cố kết $P_c$ theo phương pháp Casagrande; xác định hệ số cố kết đứng $C_v$ theo 3 phương pháp song song: Phương pháp logarit thời gian ($\log t$ - Casagrande), căn bậc hai thời gian ($\sqrt{t}$ - Taylor) và phương pháp kết thúc cố kết sơ cấp $3t$ theo tiêu chuẩn Nhật Bản JGS-2000.
  • Thí nghiệm nén cố kết thấm ngang trong phòng ($C_{h(p)}$): Gia công mẫu nén theo phương nằm ngang trong tự nhiên để đo trực tiếp độ thấm và tốc độ thoát nước ngang.
  • Thí nghiệm sức kháng cắt: Thực hiện cắt phẳng trực tiếp, nén ba trục theo sơ đồ không cố kết - không thoát nước ($UU$), cố kết - không thoát nước có đo áp lực nước lỗ rỗng ($CU$).
  • Thí nghiệm hiện trường & Triangulation: Đối chứng chéo (triangulation) đa nguồn giữa thí nghiệm cắt cánh hiện trường ($VST$), thí nghiệm xuyên tiêu chuẩn ($SPT$), thí nghiệm trong phòng và quan trắc lún thực tế ngoài hiện trường.
                    ┌─────────────────────────────────────────┐
                    │      HỆ THỐNG THU THẬP & ĐỐI CHỨNG      │
                    │         DỮ LIỆU ĐỊA KỸ THUẬT            │
                    └────────────────────┬────────────────────┘
                                         │
        ┌────────────────────────────────┼────────────────────────────────┐
        ▼                                ▼                                ▼
┌───────────────────────┐    ┌───────────────────────┐    ┌───────────────────────┐
│ THÍ NGHIỆM TRONG PHÒNG│    │ THÍ NGHIỆM HIỆN TRƯỜNG│    │  QUAN TRẮC THỰC TẾ    │
├───────────────────────┤    ├───────────────────────┤    ├───────────────────────┤
│• Nén Oedometer        │    │• Cắt cánh VST (Su)    │    │• Mốc đo lún bàn (SP)  │
│  (Cv: log t, √t, 3t)  │    │• Xuyên tiêu chuẩn SPT │    │• Quan trắc ALNLR      │
│• Nén ngang Ch(p)      │    │• Xuyên tĩnh CPTu      ││  theo thời gian      │
│• Nén 3 trục (UU, CU)  │    │• Địa chất thủy văn    │    │• Dịch chuyển ngang    │
│• Kiểm soát mẫu εv < 4%│    │  (116 giếng quan trắc)│    │  thân đắp             │
└──────────┬────────────┘    └───────────┬───────────┘    └───────────┬───────────┘
           │                             │                            │
           └─────────────────────────────┼────────────────────────────┘
                                         ▼
                    ┌─────────────────────────────────────────┐
                    │      PHÂN TÍCH NGƯỢC (ASAOKA METHOD)    │
                    │    Xác lập tương quan: Ch, Cv, m, Su    │
                    │    Mô phỏng phần tử hữu hạn PLAXIS      │
                    └─────────────────────────────────────────┘

Data và phân tích

Dữ liệu khảo sát được xử lý bằng các công cụ phần mềm chuyên dụng:

  • ArcMap 10.2: Thành lập cơ sở dữ liệu không gian, biên tập và xuất bản hệ thống bản đồ chuyên đề tỷ lệ 1/50.000 gồm: Bản đồ Điều kiện Địa kỹ thuật, Bản đồ Phân bố Đất yếu, Bản đồ Cấu trúc Nền.
  • Plaxis 8.5: Mô phỏng số phần tử hữu hạn (FEM) 2D để phân tích ứng suất - biến dạng, dự báo độ lún tức thời ($S_i$), lún cố kết ($S_c$), độ lún theo thời gian ($S_t$), xác định chiều cao đắp bù lún ($\Delta H$) và kiểm toán ổn định mái dốc ($K_{od}$) khi gia tải nền đắp trên nền tự nhiên và nền xử lý bấc thấm/giếng cát.
  • Phần mềm HRUST 1.0: Hỗ trợ tính toán ảnh hưởng của tốc độ đắp gia tải đến khả năng chịu tải của nền.
  • Phương pháp thống kê toán học: Phân tích hồi quy đơn biến và đa biến để xây dựng các hàm tương quan dự báo $C_v, C_h, S_u$ từ các chỉ tiêu vật lý cơ bản ($W, LL, PL, PI, e_0, Z$).

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

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

  1. Đặc trưng khoáng vật học và thành phần hóa học khác biệt: Khác với đất bùn sét ĐBSCL vốn chiếm ưu thế khoáng vật Montmorillonite có tính trương nở và độ nhạy rất cao, đất yếu ĐBVB Quảng Nam - Đà Nẵng có thành phần khoáng vật sét chủ yếu là Kaolinite và Illite, tiếp đến là Chlorite, hàm lượng Smectite thấp. Hàm lượng vật chất hữu cơ biến động mạnh từ $2%$ đến hơn $10%$ (các thấu kính bùn sét than sinh $ambQ_2^2$ và $abmQ_2^1$ có hàm lượng hữu cơ đạt đỉnh). Hàm lượng mất khi nung (MKN) dao động từ $4{,}5%$ đến $12{,}8%$. Chỉ số tự nén chặt của V. Priklonsky ($K_g$) nằm trong khoảng $0{,}13 - 0{,}79$, chứng minh đất sét yếu khu vực đang trong giai đoạn tự nén chặt vừa.
  2. Quy luật biến thiên và tương quan của hệ số cố kết thấm ($C_v, C_h$): Kết quả thí nghiệm nén cố kết cho thấy $C_v$ giảm mạnh khi cấp áp lực nén tăng vượt qua áp lực tiền cố kết $P_c$. Tác giả đã xác lập hệ phương trình tương quan thực nghiệm có độ tin cậy cao ($R^2 > 0{,}80$) giữa $C_v$ với các chỉ tiêu vật lý: $$C_v = f(e_0, W, LL)$$ Đặc biệt, kết quả xác định hệ số cố kết ngang từ bài toán phân tích ngược Asaoka tại mốc quan trắc SP-07 (Km12+480 cao tốc Đà Nẵng - Quảng Ngãi) và mốc SP-01 (Km1+270 đường Nguyễn Tất Thành nối dài) đã xác định chính xác tỷ số $m = C_h / C_v$ dao động thực tế từ $1{,}52$ đến $3{,}24$. Phát hiện này bác bỏ việc tùy tiện lấy $m = 2 - 5$ theo cận trên của quy phạm, vốn là nguyên nhân gây đánh giá quá lạc quan tốc độ thoát nước.
  3. Quy luật phân bố sức kháng cắt không thoát nước ($S_u, c_u$): Sức kháng cắt không thoát nước đo bằng cắt cánh hiện trường ($VST$) và nén ba trục ($UU$) tăng tuyến tính theo độ sâu $Z$ và có mối tương quan chặt chẽ với chỉ số dẻo $PI$ và hàm lượng hạt sét. Tỷ số $c_u / \sigma'_v$ của tầng bùn sét pha $mbQ_2^3$ duy trì ổn định trong khoảng $0{,}22 - 0{,}31$, phù hợp với mô hình lý thuyết của Skempton và Ladd & Foott, chứng minh sự tồn tại của cấu trúc liên kết già hóa địa chất.
  4. Xác lập 4 đơn vị cấu trúc nền (CTN) chuyên biệt: Luận án đã hệ thống hóa điều kiện địa chất thành 4 loại cấu trúc nền phục vụ xây dựng đường:
    • Loại I: Nền có lớp đất yếu nằm trực tiếp trên mặt, bề dày mỏng đến trung bình ($< 5\text{ m}$).
    • Loại II: Nền có lớp đất yếu dày ($5 - 15\text{ m}$) nằm dưới lớp phủ cát/sét xốp mỏng.
    • Loại III: Nền có tầng đất yếu xen kẹp dạng thấu kính đa tầng, xen kẽ các lớp cát bột bão hòa nước.
    • Loại IV: Nền có lớp than bùn, bùn sét chứa hữu cơ dày ($> 15 - 25\text{ m}$) kéo dài xuống sâu, chịu ảnh hưởng của tầng chứa nước có áp Pleistocen ($qp$).
  5. Hiệu năng cơ học vượt trội của giải pháp bấc thấm kết hợp gia tải: Mô phỏng số Plaxis 8.5 và kiểm chứng hiện trường cho thấy giải pháp bấc thấm (PVD) kết hợp đắp gia tải trước và đắp bù lún $\Delta H$ đạt hiệu quả xử lý độ cố kết $U_t \ge 90%$ nhanh hơn $15 - 25%$ so với giếng cát truyền thống, đồng thời giảm thiểu tối đa hiện tượng xáo động nền xung quanh lỗ cắm.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung cơ sở dữ liệu cơ học đất khu vực miền Trung vào bản đồ địa kỹ thuật quốc gia; kiểm chứng tính đúng đắn của các mô hình cố kết đa chiều Carrillo - Barron trong môi trường trầm tích ven biển hẹp và dốc.
  • Về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa phân tích ngược Asaoka, thí nghiệm ba trục chuyên sâu và phần mềm số hóa GIS/FEM trong đánh giá tương tác nền - công trình tuyến.
  • Về ứng dụng thực tiễn: Cung cấp bảng chỉ tiêu cơ lý kiến nghị chuẩn xác cho các kỹ sư tư vấn thiết kế khi lập dự án giao thông qua vùng đất yếu Quảng Nam - Đà Nẵng; hạn chế tối đa việc sử dụng hệ số kinh nghiệm thiếu kiểm chứng.
  • Về chính sách và tiêu chuẩn: Kiến nghị Bộ Giao thông Vận tải và Bộ Xây dựng cập nhật, điều chỉnh quy trình 22TCN 262-2000 và TCVN 9355-2012 theo hướng bắt buộc tính toán phân tích lún theo thời gian gắn liền với việc xác định $C_h$ thực nghiệm và phân tích vùng xáo động bấc thấm.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 hạn chế khoa học:

  1. Chiều sâu khảo sát bị khống chế ở mức 30 m, chưa bao quát hết các đới tiếp xúc sâu với đá gốc tại các thung lũng cổ chôn ngầm.
  2. Dữ liệu thí nghiệm nén cố kết ba trục có đo áp lực nước lỗ rỗng ($CU$) và nén thấm ngang trong phòng số lượng mẫu còn giới hạn do kinh phí và thiết bị thí nghiệm chuyên sâu đòi hỏi quy trình bảo dưỡng mẫu vô cùng phức tạp.
  3. Chưa mô phỏng được tác động dài hạn của tải trọng xe chạy lặp chu kỳ cao (dynamic cyclic loading) kết hợp với dao động mực nước ngầm theo mùa lũ lụt miền Trung lên sự tích lũy biến dạng dẻo và hóa lỏng cục bộ của các thấu kính cát bụi bão hòa nước xen kẹp.
  4. Mức độ chi tiết của bản đồ địa kỹ thuật 1/50.000 phù hợp cho giai đoạn quy hoạch tuyến và thiết kế cơ sở, nhưng ở giai đoạn thiết kế bản vẽ thi công từng công trình cầu cống cục bộ vẫn đòi hỏi khảo sát chi tiết tỷ lệ 1/500 - 1/2.000.

Hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Nghiên cứu biến dạng động và hiện tượng suy giảm cường độ chống cắt của đất yếu dưới tác dụng của tải trọng trùng phục giao thông tốc độ cao.
  • Ứng dụng kỹ thuật hiển vi điện tử quét (SEM) và chụp cắt lớp vi tính (Micro-CT) để lượng hóa sự biến đổi vi cấu trúc lỗ rỗng trong quá trình cố kết thấm.
  • Nghiên cứu ứng dụng công nghệ cố kết hút chân không kết hợp bấc thấm không gia tải đắp cho các khu vực đô thị có mặt bằng chật hẹp.
  • Phát triển mô hình trí tuệ nhân tạo (AI/Machine Learning) dự báo độ lún cố kết phi tuyến dựa trên chuỗi dữ liệu lớn quan trắc thời gian thực.

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

  • Tác động học thuật: Công trình là tài liệu tham khảo nền tảng trong nghiên cứu Địa kỹ thuật Xây dựng tại miền Trung, mở ra hướng nghiên cứu phân tích ngược kết hợp cơ cấu trầm tích Đệ Tứ, ước tính đóng góp hàng chục trích dẫn khoa học trên các tạp chí chuyên ngành như Tạp chí Khoa học Kỹ thuật Thủy lợi và Môi trường, Tạp chí Địa chất, Tạp chí Cầu đường.
  • Chuyển giao công nghệ và ngành giao thông: Trực tiếp phục vụ công tác thiết kế, giám sát thi công và xử lý sự cố cho các dự án giao thông huyết mạch: Đường cao tốc Đà Nẵng - Quảng Ngãi, tuyến đường ven biển Hoàng Sa - Trường Sa, đường Nguyễn Tất Thành nối dài, các khu công nghiệp Điện Nam - Điện Ngọc, Chu Lai, Hòa Khánh.
  • Hiệu quả kinh tế - xã hội: Việc tối ưu hóa mạng lưới bấc thấm và dự báo chính xác chiều cao đắp bù lún giúp rút ngắn thời gian chờ lún từ $3 - 6$ tháng, tiết kiệm hàng chục tỷ đồng chi phí xử lý nền cho mỗi dự án giao thông cấp cao, đồng thời ngăn ngừa triệt để các tai biến địa kỹ thuật phá hoại mặt đường.

Đố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: Tiếp cận phương pháp luận phân tích hệ thống (Địa hệ TN - KT), bộ số liệu thực nghiệm chuẩn mực và quy trình xử lý phân tích ngược hiện trường.
  • Kỹ sư Tư vấn Thiết kế Cầu đường: Sở hữu bảng thông số cơ lý kiến nghị chuẩn xác, công thức tương quan thực nghiệm địa phương và hướng dẫn áp dụng phần mềm Plaxis 8.5 để kiểm soát lún và ổn định.
  • Ban Quản lý Dự án và Doanh nghiệp Thi công Hạ tầng: Có cơ sở khoa học để lựa chọn giải pháp xử lý nền tối ưu (PVD, giếng cát, đệm cát, cọc xi măng đất), kiểm soát tiến độ thi công gia tải và nghiệm thu độ lún dư tin cậy.
  • Cơ quan Quản lý Nhà nước (Sở GTVT, Sở Xây dựng Đà Nẵng và Quảng Nam): Sử dụng bản đồ phân bố đất yếu và cấu trúc nền 1/50.000 làm công cụ quy hoạch không gian xây dựng và quản lý rủi ro thiên tai, tai biến địa chất công trình.

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?

Luận án đã mở rộng Lý thuyết Địa hệ Tự nhiên - Kỹ thuật (Bondarik, 1981; Nguyễn Thanh, 1984) và Lý thuyết Cố kết Thấm Terzaghi - Barron bằng cách xác lập bản chất cơ chế quá cố kết biểu kiến ($OCR > 1$) của bùn sét miền Trung hình thành do quá trình nén thứ sinh và xi măng hóa địa chất, đồng thời chứng minh mối quan hệ hữu cơ giữa thành phần khoáng vật (ưu thế Kaolinit/Illit) với đặc tính nén lún bất đẳng hướng.

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

Khác với các nghiên cứu trước đây tại Việt Nam (thường chỉ dùng thí nghiệm nén lún một trục thông thường hoặc tra bảng kinh nghiệm theo 22TCN 262-2000), luận án đã kết hợp đồng bộ thí nghiệm nén thấm ngang trong phòng ($C_{h(p)}$), thí nghiệm ba trục chuyên sâu ($UU, CU$), thí nghiệm cắt cánh hiện trường ($VST$) và đặc biệt là phương pháp phân tích ngược Asaoka từ dữ liệu quan trắc lún thực tế (mốc SP-01, SP-07) để chuẩn hóa hệ số $C_h$ và tỷ số $m = C_h / C_v$.

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

Tỷ số cố kết ngang trên đứng thực tế ($m = C_h / C_v$) của đất yếu ĐBVB Quảng Nam - Đà Nẵng chỉ dao động trong khoảng $1{,}52 - 3{,}24$, hoàn toàn không đạt mức $4 - 5$ như các đơn vị thiết kế thường giả định theo khung quy phạm 22TCN 262-2000. Việc giả định $m$ quá cao trước đây đã dẫn đến tính toán thời gian cố kết nhanh hơn thực tế, gây lún dư nghiêm trọng sau khi thảm bê tông nhựa.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Hoàn toàn có. Luận án mô tả chi tiết: (1) Tiêu chuẩn kiểm soát chất lượng lấy mẫu nguyên trạng ($\varepsilon_v < 4%$); (2) Quy trình 3 bước xác định $C_v$ theo Casagrande, Taylor và JGS-2000; (3) Thuật toán vẽ đường thẳng Asaoka từ chuỗi số liệu quan trắc lún theo các khoảng thời gian đều đặn $\Delta t$; (4) Các bước thiết lập mô hình hình học, gán thông số Hardening Soil/Soft Soil trong Plaxis 8.5.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Luận án vạch ra lộ trình nghiên cứu mở rộng: Giai đoạn 1 (Nghiên cứu vi cấu trúc và ảnh hưởng của tải trọng động giao thông); Giai đoạn 2 (Xây dựng bản đồ địa kỹ thuật số 3D tích hợp BIM/GIS cho toàn vùng kinh tế trọng điểm miền Trung); Giai đoạn 3 (Ứng dụng vật liệu sinh học - vi sinh canxi hóa MICP và công nghệ hút chân không sâu để xử lý đất yếu thân thiện môi trường).

Kết luận

  1. Làm sáng tỏ một cách toàn diện, đồng bộ và có hệ thống nguồn gốc thành tạo, quy luật phân bố không gian và bề dày của các thành tạo trầm tích đất yếu Đệ Tứ ($mbQ_2^3, ambQ_2^2, abmQ_2^1$) tại dải đồng bằng ven biển Quảng Nam - Đà Nẵng.
  2. Xác định chi tiết đặc trưng thành phần vật chất: Thành phần khoáng vật sét ưu thế Kaolinite và Illite, hàm lượng hữu cơ từ $2 - 10%$, hàm lượng mất khi nung $4{,}5 - 12{,}8%$, chỉ số tự nén chặt $K_g = 0{,}13 - 0{,}79$, phản ánh trạng thái tự nén chặt vừa.
  3. Thiết lập hệ thống thông số biến dạng - cố kết thấm chính xác: Xác định hệ số cố kết đứng $C_v$ qua 3 phương pháp ($\log t, \sqrt{t}, 3t$), xác định hệ số cố kết ngang $C_h$ và tỷ số $m = C_h / C_v = 1{,}52 - 3{,}24$ thông qua phân tích ngược quan trắc hiện trường Asaoka tại các dự án trọng điểm.
  4. Xây dựng hệ phương trình hồi quy tương quan thực nghiệm chặt chẽ giữa sức kháng cắt không thoát nước ($S_u, c_u$), hệ số cố kết ($C_v, C_h$) với các chỉ tiêu vật lý ($W, LL, PL, PI, e_0, Z$).
  5. Phân chia 4 đơn vị cấu trúc nền địa kỹ thuật và thành lập hệ thống bản đồ chuyên đề tỷ lệ 1/50.000 bằng công nghệ GIS (ArcMap 10.2), làm luận cứ khoa học tin cậy cho công tác quy hoạch, thiết kế và thi công nền đường giao thông.
  6. Ứng dụng mô phỏng số phần tử hữu hạn (Plaxis 8.5) kiểm chứng hiệu quả kỹ thuật - kinh tế của giải pháp bấc thấm (PVD) kết hợp gia tải và đắp bù lún, cung cấp giải pháp xử lý nền tối ưu, nâng cao tuổi thọ khai thác và an toàn bền vững cho các công trình hạ tầng giao thông duyên hải miền Trung.