Tổng quan về luận án

Dải đồng bằng ven biển Quảng Nam - Đà Nẵng giữ vị trí chiến lược trọng yếu trong Vùng kinh tế trọng điểm miền Trung. Sự phát triển vượt bậc của hạ tầng giao thông đô thị và liên vùng đòi hỏi phải xây dựng các tuyến huyết mạch quy mô lớn qua những khu vực có điều kiện địa chất công trình phức tạp. Khác với các bồn trũng châu thổ rộng lớn và tương đối đồng nhất tại Đồng bằng Bắc Bộ (ĐBBB) hay Đồng bằng sông Cửu Long (ĐBSCL), địa tầng Đệ Tứ vùng duyên hải miền Trung đặc trưng bởi sự phân dị mạnh mẽ, địa hình chia cắt, chịu sự chi phối của các chu kỳ biển tiến - biển lùi và quá trình bồi tích đa nguồn gốc. Sự tồn tại của các thành tạo đất yếu trẻ tuổi Holocen và Pleistocen muộn với tính nén lún cao, khả năng thoát nước kém, cường độ chịu cắt thấp đã dẫn đến hàng loạt sự cố lún sụt nghiêm trọng, phá hoại trượt mái taluy và biến dạng kéo dài trên các tuyến công trình trọng điểm như đường tránh Đà Nẵng, tuyến ven sông Tuyên Sơn - Túy Loan hay đường cao tốc Đà Nẵng - Quảng Ngãi.

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 công tác khảo sát và thiết kế xử lý nền đường trước đây hầu như áp dụng máy móc các thông số cơ lý từ ĐBBB hoặc ĐBSCL, hoặc dựa vào các bảng tra kinh nghiệm mà thiếu vắng các nghiên cứu chuyên sâu, đồng bộ về thành phần vật chất, cấu trúc nền và đặc tính cố kết thấm bất đẳng hướng. Đặc biệt, thông số hệ số cố kết thấm ngang ($C_h$) và tỷ số biến dạng thấm $m = C_h/C_v$ vốn giữ vai trò quyết định đến khoảng cách bố trí bấc thấm (PVD) hoặc giếng cát (SD) thường bị giả định sơ lược theo quy chuẩn 22TCN262-2000 với khoảng dao động quá rộng ($C_h \approx 2 - 5 C_v$), dẫn đến tính trạng sai lệch lớn giữa dự báo lún cố kết và độ lún quan trắc thực tế tại hiện trường. Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Thị Ngọc Yến với tiêu đề "Nghiên cứu tính chất cơ lý nền đất yếu đồng bằng ven biển Quảng Nam - Đà Nẵng phục vụ xây dựng đường giao thông" (chuyên ngành Địa kỹ thuật xây dựng, mã số 9580211, Trường Đại học Thủy lợi, 2018) đã giải quyết triệt để bài toán này.

Các câu hỏi và giả thuyết nghiên cứu được thiết lập cụ thể:

  • RQ1: Nguồn gốc trầm tích và thành phần vật chất (khoáng vật, hóa học, cấp phối hạt, vật chất hữu cơ) chi phối như thế nào đến tính chất cơ lý và khả năng biến dạng của các phức hệ đất yếu tại khu vực Quảng Nam - Đà Nẵng?
  • RQ2: Đặc tính cố kết thấm theo phương ngang ($C_h$) và tỷ số bất đẳng hướng $m = C_h/C_v$ xác định qua thí nghiệm nén mẫu nằm ngang trong phòng tương thích ra sao với phân tích ngược từ số liệu quan trắc hiện trường?
  • RQ3: Cấu trúc nền (CTN) đất yếu khu vực được phân chia thành những đơn vị địa kỹ thuật nào và giải pháp xử lý nền tương ứng tối ưu theo góc độ kinh tế - kỹ thuật là gì?
  • H1: Tỷ số $m = C_h/C_v$ của đất yếu duyên hải miền Trung biến thiên có quy luật theo nguồn gốc trầm tích và cấp áp lực nén, khác biệt rõ rệt so với đất bùn sét ĐBBB và ĐBSCL.
  • H2: Việc tích hợp mô hình địa hệ Tự nhiên - Kỹ thuật (TN-KT) với phân tích số phần tử hữu hạn (FEM) và bản đồ GIS chuyên đề $1/50.000$ cho phép tối ưu hóa chính xác chiều cao đắp bù lún và cự ly thiết bị thoát nước thẳng đứng.

Khung lý thuyết của công trình dựa trên lý thuyết cố kết thấm một chiều cổ điển của Karl Terzaghi (1925), định lý Carrillo (1942) về phân rã bài toán cố kết không gian ba chiều, lý thuyết dòng thấm hướng tâm của Barron (1948), cơ chế phá hoại cắt Mohr-Coulomb, cùng lý thuyết Địa hệ Tự nhiên - Kỹ thuật (Bondarik, 1986; Trofimov, 2002). Phạm vi nghiên cứu trải dài dọc dải ven biển từ huyện Núi Thành (Quảng Nam) đến quận Liên Chiểu (Đà Nẵng) với chiều sâu khảo sát khống chế tới $30\text{ m}$. Công trình đã lập bản đồ cấu trúc nền địa kỹ thuật tỷ lệ $1/50.000$, thiết lập hệ phương trình tương quan thực nghiệm có độ tin cậy cao, mang lại giá trị định lượng trực tiếp giúp rút ngắn thời gian xử lý nền và tiết kiệm đáng kể kinh phí đầu tư hạ tầng giao thông.

Literature Review và Positioning

Nghiên cứu cơ chế ứng xử cơ học của đất loại sét mềm yếu bão hòa nước là một trong những chủ đề trọng tâm của cơ học đất hiện đại. Nền tảng cố kết thấm được định hình từ mô hình nén một trục không nở hông của Terzaghi (1925), sau đó được mở rộng bởi Florin (1948) với mô hình lực thể tích xét đến tính từ biến cốt đất và độ nén ép của nước lỗ rỗng, cùng quan điểm của Verigin (1952) về sự phân phối ứng suất ban đầu đồng thời giữa hạt đất và pha lỏng. Về phương diện xác định thông số tiền cố kết ($P_c$) và hệ số cố kết ($C_v$), phương pháp "đường cong phù hợp" bán logarit của Casagrande (1938) và phương pháp căn bậc hai thời gian của Taylor (1942) đã trở thành các tiêu chuẩn quốc tế kinh điển (ASTM D2435, BS 1377). Đến đầu thế kỷ 21, Hiroyuki Tanaka (2002) và các nhà khoa học Nhật Bản đã hoàn thiện phương pháp $3t$ trong tiêu chuẩn JGS-2000, chứng minh chỉ số nén $C_c$ có tương quan hàm mũ 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 lớn vào chỉ số dẻo ($PI$).

Một cuộc tranh luận khoa học lớn trong địa kỹ thuật là sự bất đẳng hướng của tính thấm và cố kết ($C_h$ so với $C_v$). Thí nghiệm thấm ngang của Rowe, Tavenas (1983) và Seah (2004) chỉ ra rằng cấu trúc phân lớp vi mô (laminated structure) của trầm tích sông - biển làm tăng đáng kể tính thấm ngang. Tuy nhiên, việc xác định $C_h$ trong phòng thường bị chỉ trích do hiệu ứng nhiễu loạn mẫu và kích thước buồng nén hạn chế. Để giải quyết mâu thuẫn này, Bergado et al. (2002) tại Viện Công nghệ Châu Á (AIT) khi nghiên cứu bùn sét mềm Bangkok trên nền gia tải xử lý bằng bấc thấm PVD đã áp dụng bài toán phân tích ngược từ quan trắc hiện trường, xác định $C_h = 2,8\text{ m}^2/\text{năm}$ với giả thiết tỷ số xáo động $d_s/d_w = 2,0$ và $k_h/k_s = 5,0$. Ngược lại, Seah et al. (2004) trong điều kiện tương tự tại Thái Lan lại tính ra $C_h = 0,75\text{ m}^2/\text{năm}$ với $k_h/k_s = 1,4$. Sự khác biệt sâu sắc giữa các kết quả quốc tế khẳng định tính cấp thiết của việc hiệu chỉnh thông số thực nghiệm theo từng môi trường địa chất cụ thể.

                  ┌────────────────────────────────────────────────────────┐
                  │ TỔNG THỂ ĐỊA HỆ TỰ NHIÊN - KỸ THUẬT (TN-KT)            │
                  └───────────────────────────┬────────────────────────────┘
                                              │
                     ┌────────────────────────┴────────────────────────┐
                     ▼                                                 ▼
     ┌───────────────────────────────┐                 ┌───────────────────────────────┐
     │   MÔI TRƯỜNG ĐỊA KỸ THUẬT     │                 │       TÁC ĐỘNG NHÂN SINH      │
     │ - Trầm tích đa nguồn gốc      │                 │ - Tải trọng nền đắp (Htk)     │
     │ - Khoáng vật sét (Illite,     │                 │ - Gia tải nén trước           │
     │   Kaolinite, Chlorite)        │                 │ - Bấc thấm (PVD) / Giếng cát  │
     │ - Anisotropy: Ch/Cv = 1.3-3.5 │                 │ - Bơm hút chân không          │
     └───────────────┬───────────────┘                 └───────────────┬───────────────┘
                     │                                                 │
                     └────────────────────────┬────────────────────────┘
                                              │ Tương tác phi tuyến & Thoát nước lỗ rỗng
                                              ▼
                             ┌──────────────────────────────────┐
                             │ ỨNG XỬ BIẾN DẠNG & CỐ KẾT THẤM   │
                             │ - Phân tích ngược Asaoka (Ch)    │
                             │ - Mô phỏng FEM Plaxis 8.5        │
                             │ - Kiểm soát lún dư & Ổn định Kmin│
                             └──────────────────────────────────┘

Tại Việt Nam, các công trình của Đỗ Minh Toàn (1993), Nguyễn Văn Thơ (1999), Bùi Đức Hải (2006) và Nguyễn Thị Nụ (2015, 2016) đã làm sáng tỏ tính chất cơ lý của đất yếu tại ĐBBB và ĐBSCL. Đất sét mềm ĐBSCL thuộc hệ tầng Cần Giờ ($amQ_2^{2-3}cg$) có hàm lượng sét cao, chứa khoáng vật Smectit trương nở, độ nhạy cao và tỷ số $C_h/C_v$ biến động rất rộng từ $1,0$ đến $12,35$. Trong khi đó, tại duyên hải miền Trung, các công trình của Bùi Hồng Trung (2006) và Đỗ Minh Toàn et al. (2008) chỉ dừng lại ở mức độ khảo sát mô tả sơ bộ tại nội thành Đà Nẵng, chưa lượng hóa được động thái cố kết thấm và sức kháng cắt của các trầm tích biển - đầm lầy ($mbQ_2^2$), sông - biển - đầm lầy ($ambQ_2^3$) kéo dài từ Đà Nẵng tới Quảng Nam. Luận án của Nguyễn Thị Ngọc Yến định vị chính xác vào khoảng trống này, mở rộng tri thức địa kỹ thuật khu vực bằng dữ liệu thực nghiệm toàn diện và so sánh trực tiếp với các bồn trũng bùn sét nổi tiếng thế giới như sét Bangkok (Thái Lan), sét Ariake (Nhật Bản) và sét Pusan (Hàn Quốc).

Đó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 hệ thống phân loại địa kỹ thuật của Lomtadze (1984) và F.D. Maslov vào điều kiện trầm tích duyên hải miền Trung, đưa ra định nghĩa chuẩn xác mang tính định lượng cao:

"Đấ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, φ<8°), khả năng biến dạng được đặc trưng bằng hệ số nén lún lớn ($a \ge 10\text{ kPa}^{-1}$) và môđun tổng biến dạng nhỏ ($E_0 \le 5000\text{ 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 đó."

Luận án đóng góp ba mệnh đề lý thuyết cốt lõi:

  1. Mệnh đề P1 (Về tính bất đẳng hướng cố kết thấm): Tỷ số thấm và cố kết ngang $m = C_h/C_v$ trong các trầm tích bùn sét pha biển - đầm lầy ($mbQ_2^2$) và sông - biển ($ambQ_2^3$) không phải là hằng số cố định mà giảm dần phi tuyến theo cấp áp lực nén gia tăng ($\sigma'$) do sự định hướng lại của các phiến khoáng vật sét dưới tác dụng của ứng suất hữu hiệu.
  2. Mệnh đề P2 (Về tương quan sức kháng cắt không thoát nước): 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 ($S_u/\sigma'_{v0}$) tăng tuyến tính theo độ sâu, tuân theo quy luật nén chặt trọng lực trong điều kiện cố kết bình thường ($OCR \approx 1,0$), song chịu sự khống chế cục bộ bởi hàm lượng vật chất hữu cơ và mức độ xi măng hóa thứ cấp.
  3. Mệnh đề P3 (Về tương quan nén lún): Chỉ số nén $C_c$ có mối tương quan đa biến tuyến tính chặt chẽ với hệ số rỗng tự nhiên ($e_0$), độ ẩm ($W$) và giới hạn chảy ($LL$), phản ánh trực tiếp nguồn gốc bồi tích và thành phần khoáng vật Ilit - Kaolinit đặc trưng của khu vực.

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

Khung phân tích của luận án kết hợp đa tầng 3 trục tiếp cận lý thuyết:

$$\text{Địa hệ TN-KT} = \mathcal{F}\left(\text{Trầm tích học - Khoáng vật}, \text{Cơ học đất thấm bất đẳng hướng}, \text{Mô hình số FEM Biến dạng - Cố kết}\right)$$

                                    KHUNG PHÂN TÍCH ĐỊA KỸ THUẬT ĐA TẦNG
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 1. TẦNG VI MÔ (Micro-scale)                                                                           │
  │    - Phân tích khoáng vật (XRD): Illite, Kaolinite, Chlorite                                          │
  │    - Thành phần hóa học (XRF) & Vật chất hữu cơ (Mất khi nung MKN)                                    │
  │    - Cấp phối hạt (Tỷ trọng kế & Rây sàng)                                                            │
  └───────────────────────────────────────────────────┬───────────────────────────────────────────────────┘
                                                      │ Chi phối
                                                      ▼
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 2. TẦNG TRUNG MÔ (Meso-scale: Thí nghiệm cơ lý)                                                       │
  │    - Cố kết thấm đứng (Cv) & Cố kết thấm ngang (Ch) bằng Oedometer cắt mẫu nằm ngang                  │
  │    - Kháng cắt ba trục UU, CU đo áp lực nước lỗ rỗng (uw)                                             │
  │    - Sức kháng cắt hiện trường qua cắt cánh (VST) và xuyên tiêu chuẩn (SPT)                           │
  └───────────────────────────────────────────────────┬───────────────────────────────────────────────────┘
                                                      │ Tích hợp & Kiểm chứng
                                                      ▼
  ┌───────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 3. TẦNG VĨ MÔ (Macro-scale: Mô phỏng & Không gian)                                                    │
  │    - Phân tích ngược Asaoka xác định Ch thực địa từ quan trắc lún bàn đo                              │
  │    - Mô hình hóa phần tử hữu hạn Biot 2D/3D (Plaxis 8.5) với Soft Soil Model                          │
  │    - Phân vùng cấu trúc nền (CTN) và thành lập bản đồ GIS ArcMap 10.2 (1/50.000)                      │
  └───────────────────────────────────────────────────────────────────────────────────────────────────────┘

Công trình xác định rõ điều kiện biên (boundary conditions): Chiều sâu ảnh hưởng thủy lực và ứng suất $Z \le 30\text{ m}$, tải trọng nền đắp chiều cao quy đổi $H_{tk} = 2,0 - 6,0\text{ m}$, mực nước ngầm áp lực tự do hoặc bán áp lực thuộc tầng chứa nước Holocen ($qh$) và Pleistocen ($qp$).

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng lập trường phương pháp luận Thực chứng tích hợp Hiện thực phê phán (Critical Realism) thông qua phương pháp Phân tích hệ thống (PTHT). Địa hệ TN-KT xem xét tương tác biện chứng giữa ba phân hệ: Phân hệ địa chất tự nhiên (đặc tính trầm tích, địa tầng, địa chất thủy văn), Phân hệ kỹ thuật công trình (hình thái nền đường, tải trọng, tiến độ đắp) và Phân hệ can thiệp cải tạo (bấc thấm PVD, đệm cát, giếng cát, gia tải nén trước).

Thiết kế nghiên cứu đa cấp độ (Multi-level design) kết hợp đồng bộ giữa khảo sát thực địa, lấy mẫu nguyên trạng, thí nghiệm phòng chuẩn hóa, quan trắc hiện trường dài hạn và mô phỏng số học FEM. Kích thước mẫu khảo sát bao quát $428$ mẫu nước phục vụ đánh giá ăn mòn hóa học, hàng trăm mẫu đất nguyên trạng đại diện cho $5$ nhóm đất đá chính trên dải duyên hải Quảng Nam - Đà Nẵng.

                              THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ VÀ ĐỐI SOÁT DỮ LIỆU
  ┌─────────────────────────┐       ┌─────────────────────────┐       ┌─────────────────────────┐
  │  KHẢO SÁT & LẤY MẪU     │       │   THÍ NGHIỆM TRONG PHÒNG│       │   QUAN TRẮC HIỆN TRƯỜNG │
  │ - 428 mẫu nước phân tích│       │ - Nén 1 trục đứng (Cv)  │       │ - Bàn đo lún định kỳ    │
  │ - Lỗ khoan sâu tới 30m  │ ────► │ - Nén mẫu ngang (Ch(tp))│ ◄───► │ - Đo ALNLR piezometer   │
  │ - Mẫu nguyên trạng      │       │ - Cắt ba trục UU, CU    │       │ - Cắt cánh VST & SPT    │
  │   (độ biến dạng ev < 4%)│       │ - Khoáng vật, Cấp hạt   │       │                         │
  └─────────────────────────┘       └────────────┬────────────┘       └────────────┬────────────┘
                                                 │                                 │
                                                 └────────────────┬────────────────┘
                                                                  ▼
                                            ┌───────────────────────────────────────────┐
                                            │      ĐỐI SOÁT & PHÂN TÍCH SỐ HỌC          │
                                            │ - Phân tích ngược Asaoka xác định Ch(ht)  │
                                            │ - Tương quan hồi quy thực nghiệm R2 > 0.85│
                                            │ - Mô phỏng FEM Plaxis 8.5 (Biot 2D)       │
                                            │ - Bản đồ hóa cấu trúc nền GIS 1/50.000    │
                                            └───────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Nhằm đảm bảo tính đại diện và tính nguyên trạng cao nhất của mẫu đất dính mềm yếu, quy trình lấy mẫu tuân thủ nghiêm ngặt tiêu chuẩn với ống lấy mẫu thành mỏng nòng đôi piston. Độ nguyên trạng của mẫu đất loại sét được đánh giá thông qua chỉ số biến dạng thể tích $\varepsilon_v$ khi cố kết đẳng hướng trong buồng ba trục ở trạng thái ứng suất hữu hiệu tự nhiên $\sigma'_{v0}$. Tất cả các mẫu đưa vào phân tích chuyên sâu đều đạt mức $\varepsilon_v \le 4%$ (đạt cấp chất lượng "rất tốt" đến "tốt" theo tiêu chuẩn Andresen và NGI).

Quy trình tam giác hóa dữ liệu (Data Triangulation) được thực hiện chặt chẽ qua 3 nguồn kiểm chứng độc lập:

  1. Xác định hệ số cố kết ngang trong phòng $C_{h(tp)}$ bằng hộp nén không nở hông Oedometer cải tiến gia công mẫu theo phương nằm ngang tự nhiên.
  2. Xác định sức kháng cắt hiện trường bằng thí nghiệm cắt cánh (VST) và thí nghiệm xuyên tiêu chuẩn (SPT).
  3. Xác định hệ số cố kết ngang thực địa $C_{h(ht)}$ thông qua bài toán phân tích ngược (Back-analysis) theo phương pháp Asaoka (1978) dựa trên chuỗi chu kỳ lún thực tế $S_i = f(S_{i-1})$ quan trắc tại các mốc bàn đo lún sâu (điển hình tại mốc SP-07 đoạn Km12+480 đường 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).

Data và phân tích

Đặc điểm thống kê mẫu và các thông số phân tích nâng cao:

  • Phân tích thạch học - khoáng vật: Phổ nhiễu xạ tia X (XRD) xác định thành phần khoáng vật sét; phương pháp nung xác định hàm lượng vật chất hữu cơ ($MKN$).
  • Phân tích biến dạng và cố kết: Xác định 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; hệ số cố kết đứng ($C_v$) theo phương pháp căn bậc hai thời gian $\sqrt{t}$ (Taylor) và $\log t$ (Casagrande).
  • Phân tích sức kháng cắt: Hệ thống nén ba trục tự động xác định các thông số tổng $\phi, c$ (sơ đồ UU) và thông số hữu hiệu $\phi', c'$ có đo áp lực nước lỗ rỗng $u_w$ (sơ đồ CU).
  • Xử lý địa thống kê và mô phỏng số: Sử dụng phần mềm ArcMap 10.2 biên tập cơ sở dữ liệu GIS các tầng địa kỹ thuật tỷ lệ $1/50.000$; phần mềm Plaxis 8.5 ứng dụng mô hình Soft Soil (SS) và Mohr-Coulomb (MC) giải phương trình cố kết thấm Biot liên kết hoàn toàn (fully coupled consolidation).

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

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

Công trình đã mang lại 5 phát hiện đột phá với số liệu minh chứng thực nghiệm phong phú:

                            QUY LUẬT BIẾN THIÊN TỶ SỐ ANISOTROPY m = Ch/Cv
  ┌────────────────────────────────────────┬─────────────────────────────────────────────────────────────┐
  │ THÀNH TẠO TRẦM TÍCH                    │ GIÁ TRỊ TỶ SỐ m = Ch/Cv VÀ ĐẶC ĐIỂM CỐ KẾT THẤM             │
  ├────────────────────────────────────────┼─────────────────────────────────────────────────────────────┤
  │ 1. Sông - Biển - Đầm lầy (ambQ2^3)     │ m = 1,45 - 2,85 (Trung bình: 1,85 - 2,20)                    │
  │    (Bùn sét pha, độ sệt IL > 0.75)     │ Ch(tp) biến thiên (1,85 - 4,20) * 10^-3 cm2/s               │
  ├────────────────────────────────────────┼─────────────────────────────────────────────────────────────┤
  │ 2. Biển - Đầm lầy (mbQ2^2)             │ m = 1,30 - 3,20 (Phụ thuộc cấp áp lực nén p)                │
  │    (Bùn sét giàu hữu cơ MKN 2-10%)     │ Khớp nối với phân tích ngược Asaoka tại SP-07 (m = 1,8-2,4) │
  ├────────────────────────────────────────┼─────────────────────────────────────────────────────────────┤
  │ 3. Sông - Đầm lầy ven biển (abmQ2^1)   │ m = 1,60 - 3,50 (Tầng dày 5 - 25m)                          │
  │    (Bùn sét, than bùn cục bộ)          │ Cố kết chậm, hệ số nén lún a >= 10 MPa^-1                   │
  ├────────────────────────────────────────┼─────────────────────────────────────────────────────────────┤
  │ 4. So sánh với ĐBBB & ĐBSCL            │ ĐBSCL: m đạt tới 1,00 - 12,35 (Quá cố kết/Khoáng Smectite)  │
  │                                        │ Miền Trung: m tập trung dải hẹp 1,30 - 3,50 (Ổn định hơn)   │
  └────────────────────────────────────────┴─────────────────────────────────────────────────────────────┘
  1. Quy luật biến thiên tỷ số cố kết thấm ngang/đứng ($m = C_h/C_v$): Kết quả thí nghiệm phòng trên mẫu nằm ngang xác định được tỷ số $m$ của bùn sét pha $ambQ_2^3$ và $mbQ_2^2$ dao động trong khoảng $1,30 - 3,20$. Khi đối chiếu với phân tích ngược từ quan trắc thực địa tại dự án Đường cao tốc Đà Nẵng - Quảng Ngãi qua đồ thị Asaoka ($S_i - S_{i-1}$), giá trị $C_h$ thực địa đạt từ $1,85\cdot 10^{-3}$ đến $3,65\cdot 10^{-3}\text{ cm}^2/\text{s}$, cho tỷ số $m = 1,65 - 2,40$. Phát hiện này bác bỏ quan điểm tùy tiện áp dụng hệ số $m = 4 - 5$ theo cận trên của tiêu chuẩn cũ, loại bỏ nguy cơ tính toán sai lệch thời gian cố kết bấc thấm.
  2. Đặc trưng khoáng vật và tính chất hóa lý phân dị: Thành phần khoáng vật sét trong bùn sét Quảng Nam - Đà Nẵng chiếm ưu thế bởi Kaolinit ($35 - 55%$) và Ilit ($25 - 40%$), hàm lượng Chlorit thấp ($5 - 15%$) và hầu như vắng bóng Montmorillonit trương nở cao (khác biệt căn bản với bùn sét ĐBSCL). Hàm lượng vật chất hữu cơ ($MKN$) biến động từ $2,15%$ đến $9,85%$, giảm dần theo chiều sâu. Nước dưới đất thuộc tầng $qh$ và $qp$ có tính ăn mòn rửa lũa ($82,47%$) và ăn mòn carbonic ($79,67%$), trong khi ăn mòn axít và sulphat chỉ xuất hiện cục bộ.
  3. Mô hình hóa tương quan thống kê thực nghiệm: Thiết lập thành công các hàm tương quan tuyến tính đa biến có độ tin cậy $R^2 > 0,85$ dự báo chỉ số nén $C_c$ và hệ số rỗng $e_0$:

$$C_c = 0,009(LL - 10) \quad \text{và} \quad C_c = 0,423(e_0 - 0,35)$$

  1. Sự biến đổi sức kháng cắt theo độ sâu: Cường độ kháng cắt không thoát nước xác định qua thí nghiệm cắt cánh $S_u$ (VST) và nén ba trục $c_u$ (UU) tăng tuyến tính theo chiều sâu với tỷ lệ $S_u/\sigma'_{v0} = 0,22 - 0,31$. Góc ma sát trong hữu hiệu đạt $\phi' = 14^\circ30' - 22^\circ15'$, cường độ lực dính hữu hiệu $c' = 5,1 - 11,2\text{ kPa}$, phản ánh đất ở trạng thái cố kết bình thường hoặc tiền cố kết nhẹ ($OCR = 1,0 - 1,2$).
  2. Cấu trúc phân vùng địa kỹ thuật phục vụ xử lý nền: Phân chia lãnh thổ nghiên cứu thành 5 đơn vị cấu trúc nền (CTN I đến CTN V) dựa trên sự phân bố tầng đất yếu so với chiều sâu vùng hoạt áp lực công trình giao thông.
                              BẢNG TỔNG HỢP CÁC CHỈ TIÊU CƠ LÝ ĐẶC TRƯNG
  ┌───────────────────────────┬──────────────┬──────────────┬──────────────┬──────────────┬──────────────┐
  │ Chỉ tiêu cơ lý            │ Ký hiệu      │ Đơn vị       │ Bùn cát pha  │ Bùn sét pha  │ Bùn sét      │
  │                           │              │              │ (ambQ2^2)    │ (mbQ2^2)     │ (abmQ2^1)    │
  ├───────────────────────────┼──────────────┼──────────────┼──────────────┼──────────────┼──────────────┤
  │ Độ ẩm tự nhiên            │ W            │ %            │ 28,4 - 34,2  │ 39,5 - 58,4  │ 62,5 - 85,0  │
  │ Giới hạn chảy             │ LL           │ %            │ 27,8 - 32,3  │ 35,0 - 45,4  │ 48,5 - 68,0  │
  │ Chỉ số dẻo                │ PI           │ %            │ 5,8 - 6,9    │ 10,5 - 19,2  │ 21,5 - 35,0  │
  │ Hệ số rỗng tự nhiên       │ e0           │ -            │ 0,76 - 0,84  │ 0,98 - 1,45  │ 1,52 - 2,15  │
  │ Khối lượng thể tích       │ gamma_w      │ kN/m3        │ 18,3 - 18,9  │ 15,8 - 17,2  │ 14,2 - 15,5  │
  │ Hệ số nén lún             │ a (0.5-1.0)  │ MPa^-1       │ 0,35 - 0,75  │ 0,85 - 1,80  │ 1,95 - 3,80  │
  │ Mô đun biến dạng          │ E0           │ kPa          │ 4500 - 6000  │ 2200 - 3800  │ 1200 - 2000  │
  │ Sức kháng cắt không TN    │ cuu / Su     │ kPa          │ 18,0 - 25,0  │ 11,0 - 18,0  │ 6,5 - 12,0   │
  │ Hệ số cố kết đứng         │ Cv           │ 10^-3 cm2/s  │ 1,85 - 3,20  │ 0,95 - 1,65  │ 0,45 - 0,85  │
  │ Tỷ số cố kết ngang/đứng   │ m = Ch/Cv    │ -            │ 1,45 - 2,10  │ 1,65 - 2,80  │ 1,80 - 3,50  │
  └───────────────────────────┴──────────────┴──────────────┴──────────────┴──────────────┴──────────────┘

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp bức tranh toàn cảnh về động thái cố kết của trầm tích duyên hải miền Trung, bổ sung nguồn dữ liệu chuẩn mực cho cơ sở dữ liệu địa kỹ thuật Đông Nam Á; kiểm chứng tính tương thích của định lý Carrillo và mô hình Barron trong môi trường đất dính phân lớp hạt mịn.
  • Về đổi mới phương pháp luận: Chuẩn hóa quy trình tích hợp giữa phương pháp nén mẫu ngang trong phòng và phân tích ngược Asaoka ngoài hiện trường để xác định $C_h$, giúp loại bỏ sai số chủ quan của người thiết kế.
  • Về thực tiễn kỹ thuật công trình: Cung cấp các thông số đầu vào chính xác tuyệt đối cho việc mô phỏng tính toán ổn định và độ lún trên Plaxis; tối ưu hóa cự ly cắm bấc thấm (rút ngắn từ $1,3 - 1,5\text{ m}$ xuống $1,0 - 1,2\text{ m}$ hoặc ngược lại tùy thuộc đơn vị cấu trúc nền), đảm bảo độ lún còn lại sau khi đưa vào khai thác $< 10 - 20\text{ cm}$ theo tiêu chuẩn đường cao tốc cấp cao.
  • Về hoạch định chính sách và quản lý: Kiến nghị Bộ Giao thông Vận tải và Sở Giao thông các tỉnh miền Trung ban hành phụ lục hướng dẫn khảo sát địa kỹ thuật chuyên biệt cho khu vực duyên hải, bắt buộc thí nghiệm nén ngang và quan trắc lún đo áp lực nước lỗ rỗng đối với nền đường đắp cao trên đất yếu.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, công trình vẫn tồn tại một số giới hạn nghiên cứu (boundary conditions & limitations):

  1. Thiết bị xuyên CPTu: Việc đo tiêu tán áp lực nước lỗ rỗng liên tục tại hiện trường bằng thiết bị Piezocone (CPTu) chưa được triển khai đồng loạt trên tất cả các vị trí do hạn chế về trang thiết bị công nghệ tại thời điểm khảo sát, chủ yếu dựa trên VST và SPT kết hợp phân tích ngược.
  2. Hiệu ứng tải trọng động giao thông: Mô hình phần tử hữu hạn Plaxis 8.5 chủ yếu phân tích bài toán cố kết thấm tĩnh dưới tải trọng đắp và xe cộ quy đổi tương đương, chưa mô phỏng hết hiệu ứng tích lũy biến dạng từ biến động (dynamic cyclic loading degradation) dưới tác động của các đoàn xe siêu trường siêu trọng chạy liên tục.
  3. Mô hình biến dạng từ biến thứ cấp ($C_\alpha$): Tác động của từ biến dẻo nhớt cốt đất trong các lớp than bùn hữu cơ dày chưa được tách biệt hoàn toàn khỏi pha cố kết sơ cấp.

Các hướng nghiên cứu tiếp nối trong tương lai (Future Research Agenda):

  • FR1: Ứng dụng mô hình nền đất nâng cao (như Creep-SCLAY1S hoặc Hypoplasticity for clays) để mô phỏng chính xác ứng xử từ biến dị hướng lâu dài của trầm tích bùn sét ven biển.
  • FR2: Triển khai đại trà thiết bị xuyên đo áp CPTu kết hợp đo sóng địa chấn (SCPTu) nhằm xây dựng bản đồ số hóa 3D các trường thông số $C_h, S_u, G_0$ liên tục theo độ sâu.
  • FR3: Nghiên cứu công nghệ xử lý đất yếu kết hợp gia tải chân không (Vacuum Consolidation) với cọc đất - xi măng đường kính lớn (CDM) cho các nút giao và mố cầu chịu tải trọng động phức tạp.

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

Nghiên cứu tạo ra tác động học thuật và kinh tế - xã hội to lớn. Về mặt học thuật, công trình là tài liệu tham khảo nền tảng cho các nghiên cứu sinh, học viên cao học và chuyên gia địa kỹ thuật công trình biển và ven bờ. Về mặt kinh tế - kỹ thuật, việc chính xác hóa hệ số $C_h$ và phân vùng cấu trúc nền đã giúp các dự án hạ tầng lớn như Cao tốc Đà Nẵng - Quảng Ngãi và các trục vành đai đô thị Đà Nẵng - Hội An tiết kiệm hàng chục tỷ đồng chi phí xử lý nền, rút ngắn thời gian chờ lún từ 3 - 6 tháng, hạn chế triệt để các rủi ro sạt trượt mái đắp trong mùa mưa bão miền Trung. Về môi trường, việc tối ưu hóa giải pháp bấc thấm kết hợp gia tải đã giảm thiểu nhu cầu đào bỏ thay đất cục bộ, bảo vệ cảnh quan tự nhiên và các hệ sinh thái đất ngập nước ven bờ.

                   CHỈ SỐ TÁC ĐỘNG VÀ ĐỊNH LƯỢNG LỢI ÍCH CỦA CÔNG TRÌNH
  ┌─────────────────────────────────┬────────────────────────────────────────────────────────────┐
  │ LĨNH VỰC TÁC ĐỘNG               │ ĐỊNH LƯỢNG TÁC ĐỘNG & KẾT QUẢ CỤ THỂ                       │
  ├─────────────────────────────────┼────────────────────────────────────────────────────────────┤
  │ 1. Tiến độ & Thời gian chờ lún  │ Rút ngắn 3 - 6 tháng nhờ tính toán chính xác Ch và cự ly PVD│
  ├─────────────────────────────────┼────────────────────────────────────────────────────────────┤
  │ 2. Chi phí đầu tư xử lý nền     │ Giảm 12 - 18% chi phí nhờ tránh đào thay đất lãng phí     │
  ├─────────────────────────────────┼────────────────────────────────────────────────────────────┤
  │ 3. Chất lượng khai thác         │ Lún dư sau 15 năm kiểm soát < 10 - 20cm, không gãy trượt   │
  ├─────────────────────────────────┼────────────────────────────────────────────────────────────┤
  │ 4. Dữ liệu bản đồ hóa           │ Thiết lập bản đồ GIS 1/50.000 bao phủ dải Núi Thành-Liên C.│
  └─────────────────────────────────┴────────────────────────────────────────────────────────────┘

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

  • Nghiên cứu sinh và Giới học thuật Địa kỹ thuật: Tiếp cận phương pháp luận phân tích hệ thống (PTHT) chuẩn mực, khai thác bộ dữ liệu cơ lý thực nghiệm đầy đủ về trầm tích Holocen - Pleistocen duyên hải miền Trung.
  • Các Kỹ sư Tư vấn Thiết kế và Nhà thầu Thi công: Sở hữu bộ thông số đầu vào đáng tin cậy ($C_h, C_v, m, c', \phi', S_u$) và hướng dẫn tính toán xử lý nền bằng PVD/giếng cát trên phần mềm chuyên dụng (Plaxis, Geo-Slope), tránh các sự cố lún nứt công trình.
  • Chủ đầu tư và Cơ quan Quản lý Nhà nước (Bộ Giao thông Vận tải, Sở Xây dựng): Có cơ sở khoa học vững chắc để phê duyệt dự án, lập định mức kinh tế - kỹ thuật, tối ưu hóa tổng mức đầu tư và ban hành quy chuẩn kỹ thuật địa phương.

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 lý thuyết cố kết thấm bất đẳng hướng của Barron (1948) và định lý Carrillo (1942) vào địa tầng trầm tích đa nguồn gốc miền duyên hải. Luận án đã lượng hóa chính xác tỷ số $m = C_h/C_v$ không phải là một hằng số cố định mà tuân theo hàm suy giảm phụ thuộc cấp áp lực nén và nguồn gốc trầm tích ($mbQ_2^2, ambQ_2^3, abmQ_2^1$), đồng thời chứng minh sự đồng quy giữa thí nghiệm nén ngang Oedometer trong phòng và phân tích ngược Asaoka hiện trường.

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

So với các nghiên cứu của Bergado et al. (2002) tại Thái Lan hay Nguyễn Thị Nụ (2015) tại ĐBSCL vốn chỉ áp dụng đơn lẻ thí nghiệm hiện trường hoặc nén lún một chiều thông thường, luận án đã thực hiện phương pháp tam giác hóa: Gia công mẫu nén ngang chuyên dụng $\rightarrow$ Thí nghiệm cắt ba trục CU đo áp lực nước lỗ rỗng $\rightarrow$ Phân tích ngược Asaoka từ số liệu quan trắc bàn đo lún sâu $\rightarrow$ Mô phỏng phần tử hữu hạn Biot (Plaxis 8.5) và bản đồ hóa GIS $1/50.000$ theo quan điểm địa hệ Tự nhiên - Kỹ thuật.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ?

Phát hiện bất ngờ nhất là bùn sét ven biển Quảng Nam - Đà Nẵng dù có độ sệt cao ($I_L > 0,75$) và hệ số rỗng lớn ($e_0 = 1,0 - 1,5$), song tỷ số cố kết thấm ngang $m = C_h/C_v$ chỉ dao động trong dải hẹp $1,30 - 3,20$, hoàn toàn không xuất hiện các giá trị cực đoan ($m = 6 - 12$) như quan sát thấy ở bùn sét ĐBSCL. Bằng chứng được xác nhận qua phổ XRD cho thấy khoáng vật sét chủ yếu là Kaolinit và Ilit dạng tấm ổn định, hàm lượng Smectit trương nở không đáng kể.

4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?

Hoàn toàn đầy đủ. Luận án mô tả chi tiết: (a) Kỹ thuật cắt gọt mẫu thí nghiệm nén theo phương ngang từ ống mẫu nguyên trạng nòng mỏng; (b) Trình tự gia tải từng cấp ứng suất $\Delta \sigma'$; (c) Thuật toán vẽ và nội suy đường cong lún theo phương pháp Asaoka để trích xuất $C_h$; (d) Các bước thiết lập thông số địa kỹ thuật đầu vào cho mô hình Soft Soil trong Plaxis.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình tập trung vào 3 trục: (1) Số hóa 3D toàn diện tầng địa kỹ thuật duyên hải miền Trung tích hợp AI/Machine Learning trong dự báo lún cố kết; (2) Nghiên cứu ứng xử suy giảm độ bền cắt của đất mềm dưới tác động của biến đổi khí hậu, nước biển dâng và xâm nhập mặn; (3) Phát triển vật liệu kết dính vô cơ thân thiện môi trường (geopolymer từ tro bay, xỉ lò cao) trong công nghệ cọc đất - xi măng xử lý nền đường ngập nước.

Kết luận

  1. Làm sáng tỏ một cách có hệ thống điều kiện địa kỹ thuật dải đồng bằng ven biển Quảng Nam - Đà Nẵng, phân chia đất đá thành 5 nhóm cơ bản và xây dựng thành công bản đồ phân bố đất yếu cùng bản đồ cấu trúc nền địa kỹ thuật tỷ lệ $1/50.000$ trên nền tảng ArcMap 10.2.
  2. Xác định chi tiết thành phần vật chất (khoáng vật Kaolinit - Ilit chiếm ưu thế, vật chất hữu cơ $2 - 10%$, hóa học nước có tính ăn mòn rửa lũa và carbonic cao) chi phối trực tiếp đến bản chất biến dạng của các phức hệ trầm tích bùn sét Holocen - Pleistocen.
  3. Lượng hóa chính xác hệ số cố kết thấm ngang ($C_h$) và tỷ số biến dạng thấm bất đẳng hướng $m = C_h/C_v = 1,30 - 3,20$ thông qua sự kết hợp chuẩn mực giữa thí nghiệm mẫu nằm ngang trong phòng và bài toán phân tích ngược Asaoka từ quan trắc lún hiện trường.
  4. Thiết lập hệ phương trình tương quan thực nghiệm có độ tin cậy cao giữa chỉ số nén $C_c$, hệ số rỗng $e_0$, độ ẩm $W$, giới hạn chảy $LL$ và sức kháng cắt không thoát nước $S_u$, cung cấp công cụ dự báo nhanh cho công tác tư vấn thiết kế.
  5. Vận dụng xuất sắc mô hình địa hệ Tự nhiên - Kỹ thuật và mô phỏng số Plaxis 8.5 để đề xuất và tối ưu hóa các giải pháp xử lý nền đường đắp (bấc thấm PVD, giếng cát SD, gia tải nén trước), giải quyết triệt để các bài toán lún nứt, trượt lở mái dốc, mang lại hiệu quả kinh tế - kỹ thuật vượt trội cho các công trình giao thông trọng điểm miền Trung.