Tổng quan về luận án

Hiện tượng trượt lở đất đá tại các vùng đồi núi nhiệt đới gió mùa là một trong những hiểm họa địa chất phức tạp và tàn phá nặng nề nhất đối với sinh mạng, cơ sở hạ tầng giao thông và kinh tế - xã hội. Tại khu vực Nam Trung Bộ, đặc biệt là địa bàn miền núi huyện Khánh Vĩnh, tỉnh Khánh Hòa (diện tích tự nhiên $1.165\text{ km}^2$, độ cao tuyệt đối biến thiên từ $30\text{ m}$ đến $1.800\text{ m}$), các đợt mưa bão kéo dài cực đoan thường xuyên kích hoạt hàng loạt khối trượt trên các sườn dốc tự nhiên và mái dốc nhân tạo cắt xẻ taluy dọc các tuyến giao thông huyết mạch. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Địa chất (Mã số: 62.01) của tác giả Nguyễn Thanh Danh, dưới sự hướng dẫn khoa học của PGS. Đậu Văn Ngọ và TS. Tạ Quốc Dũng tại Trường Đại học Bách khoa – Đại học Quốc gia TP. Hồ Chí Minh (2019), mang tính tiên phong khi giải quyết đồng thời hai bài toán cốt lõi: cơ chế vật lý - địa kỹ thuật nội tại của quá trình mất ổn định mái dốc dưới tác động thấm của nước mưa trong môi trường đất không bão hòa, và phương pháp luận thành lập bản đồ dự báo tai biến trượt lở đất ở quy mô vùng thông qua tích hợp các mô hình xác suất thống kê không gian tiên tiến trên nền tảng GIS và R.

Khoảng trống nghiên cứu (research gap) trọng tâm được luận án định vị rõ ràng: Phần lớn các nghiên cứu dự báo trượt lở trước đây tại Việt Nam và trên thế giới chủ yếu chỉ xây dựng "Bản đồ nguy cơ trượt lở" (Landslide Susceptibility Index - LSI) dựa trên các yếu tố tĩnh nội tại (nền địa chất, độ dốc, hướng dốc, cao độ, lớp phủ thực vật) hoặc định lượng yếu tố mưa một cách thuần túy cảm tính qua phương pháp chuyên gia AHP (Masanat, 1998; Phạm Văn Hùng, 2011; Nguyễn Thám, 2012). Ngược lại, các phân tích cơ học sườn dốc thường chỉ khảo sát thấm một chiều thẳng đứng đơn giản mà thiếu vắng sự liên kết ghép đôi 2 chiều giữa sự tiêu tán áp lực nước lỗ rỗng âm (độ hút dính matric suction) với sự suy giảm sức chống cắt thực tế của đất vỏ phong hóa tàn tích (Hoàng Đại Toàn, 2012; Nguyễn Văn Hoàng, 2011). Luận án đã lấp đầy khoảng trống này bằng cách tích hợp trực tiếp động lực kích hoạt bên ngoài (lượng mưa lớn nhất năm - MP và chỉ số ẩm ướt địa hình mở rộng - TWI) vào quy trình xác suất thống kê để xây dựng "Bản đồ tai biến trượt lở" (Landslide Hazard Index - LHI) đích thực.

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

  • Câu hỏi nghiên cứu 1 (RQ1): Quá trình thấm nước mưa làm biến đổi trường áp lực nước lỗ rỗng âm và sức chống cắt của đất tàn tích không bão hòa trên mái dốc như thế nào theo thời gian và cường độ mưa?
  • Câu hỏi nghiên cứu 2 (RQ2): Ngưỡng cường độ mưa và thời gian mưa tới hạn kích hoạt trạng thái mất ổn định (hệ số an toàn $\text{FS} < 1,0$) tương ứng với các tham số hình học mái dốc (chiều cao $H$, góc dốc $\alpha$) và hệ số thấm bão hòa ($k_s$) là bao nhiêu?
  • Câu hỏi nghiên cứu 3 (RQ3): Mô hình xác suất thống kê không gian nào (FR, SI, WoE, LR) và tổ hợp yếu tố ảnh hưởng nào thông qua thuật toán trung bình mô hình Bayesian (BMA) đem lại độ tin cậy và độ chính xác dự báo tối ưu trên địa bàn nghiên cứu?
  • Giả thuyết 1 (H1): Sự gia tăng độ ẩm thể tích và triệt tiêu độ hút dính $(u_a - u_w)$ trong tầng đất tàn tích không bão hòa khi gặp mưa kéo dài làm suy giảm đột ngột cường độ chống cắt theo tiêu chuẩn Mohr-Coulomb mở rộng, khiến hệ số an toàn $\text{FS}$ giảm mạnh nhất ở nhóm đất có hệ số thấm nhỏ ($k_s \approx 10^{-7}\text{ m/s}$) sau 3 đến 5 ngày mưa liên tục.
  • Giả thuyết 2 (H2): Việc tích hợp yếu tố kích hoạt động lực lượng mưa lớn nhất năm (MP) và chỉ số ẩm ướt địa hình (TWI) cùng kỹ thuật chuẩn hóa phân bố chuẩn kết hợp phân vùng độ lệch chuẩn (Standard Deviation) sẽ tối ưu hóa diện tích dưới đường cong ROC ($\text{AUC} > 0,8$) của các mô hình dự báo không gian.

Khung lý thuyết của nghiên cứu được xây dựng vững chắc trên nền tảng Cơ học đất không bão hòa của Fredlund & Rahardjo (1993), lý thuyết thấm không bão hòa 2 chiều của Richards (1931) kết hợp đường cong đặc trưng đất - nước SWCC (van Genuchten, 1980; Fredlund & Xing, 1994), phương pháp cân bằng giới hạn lát cắt tổng quát Morgenstern - Price (1965), và lý thuyết xác suất không gian Bayesian (Bonham-Carter, 1994; Hoeting et al., 1999).

Về quy mô và ý nghĩa: Cơ sở dữ liệu của luận án bao gồm 231 điểm trượt lở thực tế được điều tra, khảo sát hiện trường, giải đoán ảnh vệ tinh độ phân giải cao và Google Earth trên toàn bộ diện tích $1.165\text{ km}^2$ huyện Khánh Vĩnh; kết hợp cùng hệ thống 11 lớp bản đồ chuyên đề đa nguồn. Đóng góp đột phá của luận án nằm ở việc lượng hóa chính xác ngưỡng mưa kích hoạt phá hủy trượt lở ($I \ge 10\text{ mm/h}$, thời gian mưa liên tục $t \ge 5-7\text{ ngày}$, chiều cao mái dốc $H \ge 10\text{ m}$, góc dốc $\alpha \ge 45^\circ$, hệ số thấm $k_s = 10^{-7} - 10^{-5}\text{ m/s}$), đồng thời thiết lập khung dự báo không gian đạt độ chính xác kiểm chứng định lượng vượt trội ($\text{AUC} = 0,8 - 0,9$).

Literature Review và Positioning

Các công trình nghiên cứu trượt lở đất do mưa trên thế giới đã phát triển qua hai nhánh tiếp cận chính: phương pháp thực nghiệm ngưỡng mưa kinh nghiệm và phương pháp cơ học địa kỹ thuật thấm - ổn định. Trong nhánh thực nghiệm, các công trình kinh điển của Caine (1980), Guzzetti et al. (2008), và Dahal & Hasegawa (2008) đã thiết lập các hàm tương quan dạng lũy thừa $I = \alpha \cdot D^{-\beta}$ (với $I$ là cường độ mưa, $D$ là thời gian mưa) để xác định ngưỡng xảy ra trượt lở trên quy mô toàn cầu hoặc khu vực. Tuy nhiên, tranh luận học thuật nảy sinh khi Shuin et al. (2012) chỉ ra rằng đường quan hệ $I - D$ hoàn toàn không phản ánh được trạng thái độ ẩm tích lũy tiền thời của đất và đặc tính cơ lý riêng biệt của từng loại thổ nhưỡng; do đó không thể phân biệt được các đợt mưa lớn có gây trượt và không gây trượt.

Để giải quyết mâu thuẫn này, nhánh nghiên cứu cơ học đất không bão hòa (Fredlund & Rahardjo, 1993; Gasmo et al., 2000; Tsaparas et al., 2002; Rahardjo et al., 2005, 2008, 2010; Collins & Znidarcic, 2004) đã chứng minh vai trò quyết định của gradient thủy lực ngấm Darcy và sự suy giảm áp lực nước lỗ rỗng âm. Gasmo et al. (2000) phát hiện tốc độ thấm ban đầu vào mái dốc có thể lớn gấp nhiều lần hệ số thấm bão hòa do gradient thủy lực ban đầu rất cao, làm giảm nhanh độ hút dính và sức chống cắt. Rahardjo et al. (2002, 2008) khẳng định lượng nước ngấm thay đổi từ 40% đến 100% tổng lượng mưa và lượng mưa tiền thời đóng vai trò sống còn trong việc tạo ra áp lực nước lỗ rỗng bất lợi, trong đó đất có hệ số thấm thấp đòi hỏi lượng mưa tích lũy trước đó dài ngày hơn so với đất có hệ số thấm cao. Lo & Lee (2010) cùng Tsai & Chen (2011) cũng chỉ ra rằng mô hình mưa có cường độ tăng dần hoặc cường độ không đổi là kịch bản nguy hiểm nhất đối với độ ổn định mái dốc.

Về phương diện thành lập bản đồ phân vùng nhạy cảm và tai biến trượt lở trên thế giới, sự đối lập thể hiện rõ giữa các phương pháp định tính/bán định tính chuyên gia như AHP (Analytic Hierarchy Process) vốn mang nặng tính chủ quan (Masanat, 1998; Nguyễn Thám et al., 2012) với các phương pháp định lượng xác suất thống kê tích hợp GIS như Tỷ số tần suất - FR (Lee & Talib, 2005; Yilmaz, 2009), Chỉ số thống kê - SI (Bùi Tiến Diệu et al., 2011), Trọng số chứng cứ - WoE (Bonham-Carter, 1994; Youssef et al., 2015), và Hồi quy Logistic - LR (Yilmaz, 2009; Jaafari et al., 2014).

                 TIẾP CẬN DỰ BÁO TAI BIẾN TRƯỢT LỞ
                                 │
         ┌───────────────────────┴───────────────────────┐
         ▼                                               ▼
Nhánh Cơ học Địa kỹ thuật                        Nhánh Không gian Thống kê
- Caine (1980), Guzzetti (2008):                 - Masanat (1998), Phạm Hùng (2011):
  Đường cong ngưỡng mưa I-D (thiếu ẩm tiền thời)   AHP bán định tính (chủ quan)
- Fredlund & Rahardjo (1993), Gasmo (2000):      - Lee & Talib (2005), Yilmaz (2009):
  Thấm không bão hòa 2D & tiêu tán Matric Suction  FR, WoE, LR thống kê không gian (thiếu động lực mưa)
         │                                               │
         └───────────────────────┬───────────────────────┘
                                 ▼
                    ĐỊNH VỊ CỦA LUẬN ÁN (2019)
         Ghép đôi mô hình thấm-ổn định Geostudio (SEEP/W + SLOPE/W)
                                 +
        Tích hợp yếu tố động lực Mưa (MP) & TWI mở rộng vào WoE/LR/BMA

Luận án của Nguyễn Thanh Danh đã định vị nghiên cứu tại điểm giao thoa giữa cơ học vi mô và xác suất thống kê vĩ mô. So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  1. Nghiên cứu của Isik Yilmaz (2009) tại Tokat, Thổ Nhĩ Kỳ: Yilmaz sử dụng FR, LR và ANN với 8 yếu tố thuần túy địa hình và địa chất, bỏ qua yếu tố lượng mưa cực đoan và chỉ số độ ẩm địa hình tích hợp thấm. Luận án của Nguyễn Thanh Danh đã mở rộng lên 11 yếu tố, bổ sung yếu tố động lực MP và chỉ số TWI thủy văn cải tiến, đồng thời áp dụng thuật toán BMA để loại bỏ hiện tượng overfitting trong LR mà Yilmaz chưa xử lý được.
  2. Nghiên cứu của Youssef et al. (2015) tại Lưu vực Wadi Itwad, Ả Rập Xê Út: Youssef áp dụng WoE và Dempster-Shafer trên vùng khí hậu khô hạn với lượng mưa trung bình năm. Trong khi đó, luận án tập trung vào vùng nhiệt đới gió mùa có độ ẩm cao với lớp vỏ phong hóa laterit hóa dày ($0,5 - 12\text{ m}$), đưa vào tham số lượng mưa lớn nhất năm (MP) và xác lập quy trình chuẩn hóa phân bố dữ liệu trước khi phân vùng bằng độ lệch chuẩn, giải quyết triệt để vấn đề bất đồng nhất thang đo giá trị chỉ số LHI giữa các mô hình.

Đó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à phát triển ứng dụng lý thuyết Cơ học đất không bão hòa của Fredlund, Morgenstern & Widger (1978) và Fredlund & Rahardjo (1993) vào điều kiện địa chất công trình - vỏ phong hóa nhiệt đới Nam Trung Bộ. Cụ thể, phương trình sức kháng cắt trên mặt phá hủy không bão hòa được xác lập:

$$\tau = c' + (\sigma_f - u_a)\tan\phi' + (u_a - u_w)_f \tan\phi^b$$

Trong đó, đóng góp lý thuyết then chốt là việc làm sáng tỏ quy luật suy giảm phi tuyến của thành phần lực dính mở rộng $c = c' + (u_a - u_w)\tan\phi^b$ khi độ hút dính $(u_a - u_w)$ tiệm cận về 0 dưới tác động của dòng thấm hai chiều không dừng. Luận án chứng minh bằng mô hình toán rằng sự mất ổn định mái dốc nhiệt đới không bắt nguồn từ sự gia tăng ứng suất cắt gây trượt do trọng lượng bản thân đất tăng khi ẩm hóa, mà nguyên nhân chủ yếu và quyết định là sự sụt giảm tức thời của sức kháng cắt hiệu dụng do áp lực nước lỗ rỗng âm bị triệt tiêu từ mặt thấm đi sâu vào khối trượt.

TIÊU CHUẨN PHÁ HỦY MOHR-COULOMB MỞ RỘNG TRONG ĐẤT KHÔNG BÃO HÒA:
  τ = c' + (σ_f - u_a)tan(φ') + (u_a - u_w)_f tan(φ^b)
       │         │                     │
       │         │                     └─ Thành phần tăng cường sức kháng cắt do Matric Suction
       │         └─ Ứng suất pháp thực hiệu dụng
       └─ Lực dính hiệu dụng của cốt đất (khi u_a - u_w = 0)

Luận án tạo ra sự chuyển dịch nhận thức khoa học (paradigm shift): Chuyển đổi căn bản từ việc đánh giá nguy cơ tĩnh (Landslide Susceptibility) sang đánh giá tai biến động có xét đến thời gian và biến cố kích hoạt (Landslide Hazard) bằng cách tích hợp trực tiếp tham số lượng mưa lớn nhất năm (MP) và trường thấm của đất vào các mô hình xác suất toán học.

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

Khung phân tích của luận án tích hợp đa tầng từ cấp độ phần tử vật lý đến cấp độ cảnh quan địa không gian:

  1. Mô hình ghép đôi phân tích thấm - ổn định (Coupled Seepage-Stability Analysis): Sử dụng phương pháp phần tử hữu hạn (FEM) giải phương trình thấm Richards 2 chiều không dừng trong mô đun SEEP/W, sau đó chuyển trường áp lực nước lỗ rỗng âm theo từng bước thời gian sang mô đun SLOPE/W để tính toán hệ số an toàn $\text{FS}$ theo phương pháp cân bằng giới hạn lát cắt nghiêm ngặt Morgenstern - Price, thỏa mãn đồng thời cả cân bằng lực ($F_f$) và cân bằng mô men ($F_m$).
  2. Mô hình tối ưu hóa thống kê đa biến Bayesian (Bayesian Model Averaging - BMA): Thay vì lựa chọn tổ hợp biến độc lập theo kinh nghiệm chuyên gia, luận án lần đầu tiên ứng dụng thuật toán BMA trong môi trường R để quét toàn bộ không gian mô hình của 11 yếu tố ảnh hưởng, định lượng xác suất hậu nghiệm (posterior probability) của từng biến số và xác lập các tổ hợp biến tối ưu nhất.
  3. Quy trình chuẩn hóa phân bố và phân vùng độ lệch chuẩn (Standard Deviation Normalization Framework): Giải quyết rào cản toán học khi so sánh các mô hình thống kê khác nhau bằng cách đưa các dải giá trị thô của chỉ số LSI/LHI về phân bố chuẩn trên đoạn $[0, 1]$, cho phép phân chia khách quan 5 cấp độ tai biến: Rất thấp, Thấp, Trung bình, Cao, và Rất cao.
  4. Điều kiện biên xác định (Boundary Conditions): Vỏ phong hóa nhiệt đới phát triển trên đá trầm tích lục nguyên Jura hệ tầng La Ngà ($J_{2}ln$) và đá phun trào Creta hệ tầng Đèo Bảo Lộc ($J_{3}đbl$), Nha Trang ($Knt$), Đơn Dương ($K_{2}dd$); chiều cao mái dốc $H = 10 - 30\text{ m}$; góc dốc $\alpha = 15^\circ - 60^\circ$; mực nước ngầm ban đầu ở mức sâu $H_w$; hệ số thấm bão hòa $k_s = 10^{-7} - 10^{-5}\text{ m/s}$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ lập trường nhận thức luận thực chứng sau hiện thực (post-positivism) kết hợp hiện thực phản biện (critical realism), khẳng định hiện tượng trượt lở là kết quả tương tác nhân quả khách quan giữa các điều kiện biên địa chất tĩnh và các xung lực thủy văn động, có thể đo lường và lượng hóa chính xác bằng thực nghiệm toán - lý.

Thiết kế nghiên cứu đa cấp độ (multi-level research design):

  • Cấp độ vi mô/cục bộ: Khảo sát cơ học đất không bão hòa, đo vẽ mặt cắt địa kỹ thuật taluy, thí nghiệm xác định đường cong SWCC, hàm thấm thủy lực $k(u_a - u_w)$ và các chỉ tiêu cơ lý mẫu đất tàn tích ($eQ$), sườn tích ($dQ$).
  • Cấp độ vĩ mô/vùng: Số hóa và phân tích không gian 11 lớp yếu tố trên toàn huyện Khánh Vĩnh với kích thước ô lưới raster chuẩn hóa trên nền tảng phần mềm GIS mã nguồn mở ILWIS.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được thực hiện nghiêm ngặt qua 4 bước:

  1. Xây dựng bản đồ kiểm kê trượt lở (Landslide Inventory Map): Thu thập và số hóa chính xác 231 điểm trượt lở thực tế thông qua định vị GPS khảo sát thực địa dọc các tuyến giao thông chính (tỉnh lộ 2, đoạn đèo Khánh Lê - Hòn Giao thuộc quốc lộ 27C Nha Trang - Đà Lạt tại Km 42+800, các điểm VS32, VS43, WS37, VS39, VS347, VS349), đối soát qua ảnh viễn thám và ảnh vệ tinh đa thời gian.
  2. Trích xuất và số hóa 11 biến số môi trường (Landslide Conditioning & Triggering Factors):
    • Mô hình số độ cao (DEM, m)
    • Góc dốc (SLOPE, độ)
    • Hướng dốc (ASPECT)
    • Độ uốn cong bề mặt mái dốc (CURVATURE)
    • Chỉ số ẩm ướt địa hình (TWI) tính theo Montgomery & Dietrich (1994) và Borga et al. (1998) có xét đến quan hệ thấm: $\text{TWI} = \ln(a / T\tan\beta)$
    • Thạch học (LITHOLOGY) với 5 phân hệ thạch học chính
    • Khoảng cách đến đường giao thông (DFR, m)
    • Khoảng cách đến sông suối (DFD, m)
    • Khoảng cách đến đứt gãy kiến tạo (DFF, m)
    • Chỉ số thực vật khác biệt chuẩn hóa (NDVI)
    • Lượng mưa lớn nhất năm (MP, mm) từ chuỗi số liệu quan trắc khí tượng trạm Nha Trang, Cam Ranh, Khánh Vĩnh và Khánh Sơn giai đoạn 1977 - 2016.
  3. Kiểm tra tính độc lập và hiện tượng đa cộng tuyến: Sử dụng ma trận hệ số tương quan Pearson, chỉ số dung sai (Tolerance - TOL) và hệ số phóng đại phương sai (Variance Inflation Factor - VIF) cho 11 yếu tố. Tất cả các biến được chọn đều thỏa mãn $\text{VIF} < 5$ và $\text{TOL} > 0,2$, chứng minh không tồn tại hiện tượng đa cộng tuyến gây sai lệch mô hình.
  4. Đối chuẩn tam giác (Triangulation) và Kiểm chứng mô hình (Validation): Phân chia ngẫu nhiên bộ dữ liệu 231 điểm trượt thành 2 tập độc lập: 70% dữ liệu huấn luyện (Training dataset) để xây dựng bản đồ trọng số và 30% dữ liệu kiểm chứng (Validation dataset). Đánh giá độ chính xác thông qua đường cong tỷ lệ thành công (Success Rate Curve) và đường cong tỷ lệ dự báo (Prediction Rate Curve) dựa trên chỉ số diện tích dưới đường cong (Area Under Curve - AUC).
                      QUY TRÌNH NGHIÊN CỨU TOÀN DIỆN
 ┌───────────────────────────────┐     ┌───────────────────────────────┐
 │   231 ĐIỂM TRƯỢT LỞ THỰC TẾ   │     │  11 LỚP BẢN ĐỒ YẾU TỐ ẢNH HƯỞNG│
 │ (Khảo sát GPS + Viễn thám)   │     │ (DEM, SLOPE, MP, TWI, DFR...) │
 └──────────────┬────────────────┘     └──────────────┬────────────────┘
                │                                     │
                ▼                                     ▼
        Chia ngẫu nhiên                       Kiểm tra đa cộng tuyến
       ┌────────┴────────┐                       (TOL > 0.2, VIF < 5)
       ▼                 ▼                            │
  70% Huấn luyện   30% Kiểm chứng                     │
       │                 │                            │
       └────────┬────────┘                            │
                │                                     │
                ▼                                     ▼
 ┌─────────────────────────────────────────────────────────────────────┐
 │       MÔ HÌNH HÓA XÁC SUẤT THỐNG KÊ & TỐI ƯU HÓA BMA (TRÊN R)       │
 │   - Tỷ số tần suất (FR)               - Chỉ số thống kê (SI)        │
 │   - Trọng số chứng cứ (WoE)           - Hồi quy Logistic (LR)       │
 └──────────────────────────────────┬──────────────────────────────────┘
                                    │
                                    ▼
 ┌─────────────────────────────────────────────────────────────────────┐
 │ Chuẩn hóa phân bố + Phân vùng 5 cấp độ bằng Độ lệch chuẩn (Std Dev) │
 └──────────────────────────────────┬──────────────────────────────────┘
                                    │
                                    ▼
 ┌─────────────────────────────────────────────────────────────────────┐
 │ KIỂM CHỨNG BẰNG ĐƯỜNG CONG ROC/AUC (Success Rate & Prediction Rate) │
 │                      Kết quả: AUC = 0.8 - 0.9                       │
 └─────────────────────────────────────────────────────────────────────┘

Data và phân tích

Đặc tính cơ lý của các nhóm đất vỏ phong hóa tàn tích - sườn tích tại vùng nghiên cứu được thí nghiệm và tổng hợp chính xác trong Bảng 2.2 của luận án:

Chỉ tiêu cơ lý Ký hiệu & Đơn vị Đất Sét ($eQ$) Đất Sét pha ($dQ$) Đất Sét pha lẫn dăm sạn ($pdQ$)
Thành phần hạt sạn $%$ $1$ $5$ $34$
Thành phần hạt cát $%$ $32$ $61$ $35$
Thành phần hạt bụi $%$ $29$ $14$ $13$
Thành phần hạt sét $%$ $38$ $20$ $18$
Giới hạn chảy ($LL$) $%$ $46$ $34$ $36,3$
Giới hạn dẻo ($PL$) $%$ $28$ $22$ $23,9$
Chỉ số dẻo ($PI$) $%$ $18$ $12$ $12,4$
Độ sệt ($B$) $-$ $0,06$ $-0,06$ $-0,03$
Độ ẩm tự nhiên ($W$) $%$ $29$ $20,68$ $23,53$
Dung trọng tự nhiên ($\gamma$) $\text{kN/m}^3$ $17,4$ $17,3$ $18,0$
Dung trọng bão hòa ($\gamma_{sat}$) $\text{kN/m}^3$ $18,5$ $18,9$ $19,3$
Dung trọng khô ($\gamma_d$) $\text{kN/m}^3$ $13,5$ $14,2$ $14,6$
Tỷ trọng ($G_s$) $-$ $2,68$ $2,69$ $2,73$
Độ lỗ rỗng ($n$) $%$ $0,496$ $0,472$ $0,465$
Hệ số rỗng ($e$) $-$ $0,985$ $0,894$ $0,870$
Độ bão hòa ($S_r$) $%$ $79$ $62$ $74$
Góc ma sát trong ($\phi'$) độ ($^\circ$) $16$ $17$ $18$
Lực dính hiệu dụng ($c'$) $\text{kN/m}^2$ (kPa) $21$ $15$ $18$
Hệ số thấm bão hòa ($k_s$) $\text{m/s}$ $2,09 \times 10^{-7}$ $1,55 \times 10^{-6}$ $2,03 \times 10^{-6}$

Phần mềm chuyên dụng được áp dụng bao gồm: Bộ công cụ GeoStudio 2007 (SEEP/W cho phân tích phần tử hữu hạn thấm không dừng; SLOPE/W cho phân tích hệ số an toàn cân bằng giới hạn); Hệ thống thông tin địa lý mã nguồn mở ILWIS; và Gói phần mềm thống kê R (thực hiện mô hình BMA với gói BMS / bicreg).

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

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

  1. Xác lập định lượng ngưỡng mưa kích hoạt phá hủy mái dốc theo thời gian: Kết quả mô phỏng số ghép đôi SEEP/W và SLOPE/W cho thấy hệ số an toàn $\text{FS}$ suy giảm mạnh theo hàm phi tuyến đối với thời gian mưa và cường độ mưa. Khi cường độ mưa đạt ngưỡng $I = 10\text{ mm/h}$, sau thời gian mưa liên tục kéo dài từ 5 đến 7 ngày, trường áp lực nước lỗ rỗng âm trong mái dốc bị triệt tiêu hoàn toàn, kéo tụt hệ số an toàn từ $\text{FS} > 1,35$ xuống dưới ngưỡng cân bằng giới hạn ($\text{FS} < 1,0$), dẫn đến phá hủy sạt trượt mái dốc.

  2. Nghịch lý độ nhạy hệ số thấm ($k_s$) đối với thời gian mưa: Phân tích cơ học chỉ ra vai trò then chốt có tính phản trực giác của hệ số thấm đất: Với các trận mưa kéo dài trên 3 ngày, hệ số an toàn của các mái dốc cấu tạo bởi đất có hệ số thấm nhỏ ($k_s \approx 2,09 \times 10^{-7}\text{ m/s}$, đất sét tàn tích) suy giảm mạnh hơn và duy trì trạng thái bão hòa bất lợi lâu hơn so với mái dốc cấu tạo bởi đất có hệ số thấm lớn ($k_s \approx 2,03 \times 10^{-6}\text{ m/s}$, đất sét pha lẫn dăm sạn). Đất thấm lớn cho phép nước ngấm nhanh nhưng cũng thoát nước nhanh, trong khi đất thấm bé giữ nước lại ở tầng sát mặt, tạo nên mặt trượt nông vô cùng nguy hiểm.

  3. Mối tương quan hình học mái dốc và nguy cơ trượt lở: Nguy cơ phá hủy sạt lở tăng vọt khi chiều cao mái dốc $H \ge 10\text{ m}$ và góc dốc $\alpha \ge 45^\circ$. Với góc dốc $\alpha < 30^\circ$, ngay cả dưới vũ lượng mưa cực đại kéo dài 7 ngày, hệ số an toàn $\text{FS}$ vẫn duy trì $> 1,15$.

  4. Xếp hạng độ chính xác vượt trội của các mô hình xác suất thống kê: Đánh giá định lượng trên 231 điểm trượt lở bằng đường cong ROC cho thấy cả 4 phương pháp đều đạt độ chính xác xuất sắc với $\text{AUC} = 0,8 - 0,9$. Trong đó, phương pháp Trọng số chứng cứ (WoE) đạt hiệu quả dự báo cao nhất, kế tiếp là Tỷ số tần suất (FR), Chỉ số thống kê (SI), và Hồi quy Logistic (LR). Mặc dù FR và SI có độ chính xác thấp hơn một bậc không đáng kể, nhưng có ưu thế thuật toán đơn giản, xử lý ma trận raster cực nhanh trên GIS.

  5. Xác định 4 yếu tố chi phối cốt lõi và 5 tổ hợp mô hình tối ưu bằng BMA: Ứng dụng thuật toán Bayesian Model Averaging đã chứng minh trong số 11 yếu tố khảo sát, có 8 yếu tố liên quan trực tiếp đến trượt lở, và xác lập được 4 yếu tố trọng yếu nhất gồm:

    • Cao độ (DEM)
    • Khoảng cách đến đường giao thông (DFR)
    • Lượng mưa lớn nhất năm (MP)
    • Góc dốc (SLOPE)

    Đồng thời, BMA xác lập 5 mô hình tổ hợp tối ưu:

    • Mô hình 1: DEM + DFR + MP + SLOPE (Tối ưu nhất về hiệu quả dự báo do sử dụng ít biến nhất - 4 yếu tố).
    • Mô hình 2: DEM + DFR + MP + SLOPE + LITHOLOGY.
    • Mô hình 3: DEM + DFR + MP + SLOPE + LITHOLOGY + DFF (Tối ưu nhất về độ chính xác dự báo tuyệt đối).
    • Mô hình 4: DEM + DFR + MP + SLOPE + DFF.
    • Mô hình 5: DEM + DFR + MP + SLOPE + DFD.
                    5 TỔ HỢP MÔ HÌNH TỐI ƯU TỪ BMA
┌────────────┬──────────────────────────────────────────┬────────────────────────┐
│  Mô hình   │      Tổ hợp các yếu tố ảnh hưởng         │  Đánh giá hiệu năng    │
├────────────┼──────────────────────────────────────────┼────────────────────────┤
│ Mô hình 1  │ DEM + DFR + MP + SLOPE                   │ TỐI ƯU HIỆU QUẢ        │
│            │ (4 yếu tố cốt lõi)                       │ (Ít biến, tốc độ cao)  │
├────────────┼──────────────────────────────────────────┼────────────────────────┤
│ Mô hình 2  │ DEM + DFR + MP + SLOPE + LITHOLOGY       │ Trung gian             │
├────────────┼──────────────────────────────────────────┼────────────────────────┤
│ Mô hình 3  │ DEM + DFR + MP + SLOPE + LITHOLOGY + DFF │ TỐI ƯU ĐỘ CHÍNH XÁC    │
│            │ (6 yếu tố đầy đủ)                        │ (AUC cao nhất)         │
├────────────┼──────────────────────────────────────────┼────────────────────────┤
│ Mô hình 4  │ DEM + DFR + MP + SLOPE + DFF             │ Mở rộng đứt gãy        │
├────────────┼──────────────────────────────────────────┼────────────────────────┤
│ Mô hình 5  │ DEM + DFR + MP + SLOPE + DFD             │ Mở rộng thủy văn       │
└────────────┴──────────────────────────────────────────┴────────────────────────┘

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Công trình chứng minh tính tất yếu của việc tích hợp cơ học đất không bão hòa vào các bài toán phân vùng tai biến địa chất diện rộng; cung cấp cơ sở chuyển đổi từ LSI sang LHI một cách toán học chặt chẽ.
  • Về mặt phương pháp luận: Đề xuất quy trình chuẩn mực: Chuẩn hóa phân bố dữ liệu $\rightarrow$ Phân vùng độ lệch chuẩn (Standard Deviation) $\rightarrow$ Tối ưu hóa biến bằng BMA $\rightarrow$ Kiểm chứng chéo bằng AUC. Quy trình này hoàn toàn có thể nhân rộng cho các vùng địa chất phức tạp khác tại Đông Nam Á và Nam Á.
  • Về ứng dụng thực tiễn và quy hoạch giao thông: Cung cấp tọa độ và khoanh vùng chính xác các điểm nóng tai biến dọc tuyến Quốc lộ 27C (đoạn đèo Khánh Lê - Hòn Giao) và Tỉnh lộ 2, giúp ngành giao thông vận tải thiết kế chính xác hệ thống thoát nước bề mặt, gia cố taluy bằng neo ứng suất trước hoặc tường chắn đất có neo tại các mái dốc $\alpha \ge 45^\circ$, $H \ge 10\text{ m}$.
  • Về chính sách phòng chống thiên tai: Cung cấp cho Ủy ban Nhân dân huyện Khánh Vĩnh và tỉnh Khánh Hòa bản đồ số tỷ lệ lớn làm công cụ ra quyết định di dời dân cư cục bộ trong mùa bão lũ khi lượng mưa dự báo vượt ngưỡng $10\text{ mm/h}$ liên tục trong 5 ngày.

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 các giới hạn khoa học:

  1. Độ phân giải của số liệu mưa không gian: Tham số lượng mưa lớn nhất năm (MP) được nội suy không gian từ 6 trạm đo mưa khí tượng khu vực (Nha Trang, Cam Ranh, Sơn Thái, Khánh Phú, Khánh Bình, Tô Hạp), chưa phản ánh trọn vẹn hiện tượng mưa tiểu vi khí hậu cục bộ tại các đỉnh núi cao trên $1.500\text{ m}$ (núi Hòn Giao, Chư Tông).
  2. Mô hình địa chất công trình đồng nhất hóa: Mô hình phân tích thấm số 2D giả định lớp đất tàn tích tương đối đồng nhất về mặt dị hướng, chưa mô phỏng chi tiết sự hiện diện của các khe nứt tế vi, hang rỗng sinh học (macropores) và đới tiếp xúc bất chỉnh hợp giữa lớp đất phong hóa với mặt đá gốc granit/ryolit chưa phong hóa.
  3. Đặc tính tiền ẩm biến đổi theo thời gian thực: Nghiên cứu sử dụng phân tích theo các kịch bản mưa tĩnh mà chưa kết nối trực tiếp với hệ thống quan trắc áp lực nước lỗ rỗng và độ ẩm đất tự động theo thời gian thực (Real-time IoT Piezometer Sensors).

Chương trình nghiên cứu tiếp nối trong 5-10 năm tới bao gồm:

  • Mở rộng mô hình thấm ghép đôi sang không gian 3 chiều (3D Transient Unsaturated Seepage) kết hợp mô hình phá hủy mặt trượt phức hợp phi tròn.
  • Tích hợp dữ liệu lượng mưa radar thời gian thực (Doppler Weather Radar) và ảnh viễn thám độ ẩm đất vệ tinh (SMAP) vào thuật toán BMA để xây dựng hệ thống cảnh báo sớm trượt lở trực tuyến theo thời gian thực (Dynamic Early Warning System).
  • Ứng dụng các thuật toán học máy và trí tuệ nhân tạo thế hệ mới (Deep Learning, XGBoost, CatBoost) lai ghép với thuật toán tối ưu hóa BMA.

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

  • Tác động học thuật: Luận án là công trình tiêu biểu tạo cầu nối giữa cơ học địa kỹ thuật công trình và khoa học thông tin địa lý môi trường tại Việt Nam, mở ra hướng nghiên cứu liên ngành chuyên sâu với tiềm năng trích dẫn cao trên các tạp chí chuyên ngành như Landslides, Engineering Geology, Geomorphology, và Catena.
  • Chuyển đổi ngành công nghiệp xây dựng và giao thông: Cung cấp tiêu chuẩn thiết kế mái dốc thực tế cho các dự án hạ tầng giao thông qua vùng đồi núi hiểm trở Nam Trung Bộ; tối ưu hóa chi phí xử lý nền móng và taluy bằng cách định lượng chính xác ngưỡng phá hủy thay vì áp dụng hệ số an toàn dư thừa quá mức.
  • Ảnh hưởng chính sách cấp tỉnh và quốc gia: Đóng góp trực tiếp vào Đề án tổng thể "Điều tra, đánh giá và phân vùng cảnh báo nguy cơ trượt lở đất đá các vùng miền núi Việt Nam" của Bộ Tài nguyên và Môi trường; giúp chính quyền địa phương tích hợp bản đồ tai biến vào quy hoạch sử dụng đất cấp huyện đến năm 2030.

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

  • Nghiên cứu sinh và Giới học thuật Địa kỹ thuật - Địa chất: Thụ hưởng cơ sở lý thuyết hoàn chỉnh về cơ học đất không bão hòa, phương trình mở rộng Mohr-Coulomb, quy trình tối ưu BMA và bộ số liệu kiểm chứng thực tế tại Khánh Vĩnh.
  • Kỹ sư Tư vấn thiết kế Giao thông và Thủy lợi: Ứng dụng trực tiếp thông số ngưỡng mưa $10\text{ mm/h}$, hệ số thấm $k_s$, góc dốc $\alpha \ge 45^\circ$, chiều cao $H \ge 10\text{ m}$ để tính toán giải pháp ổn định taluy đường và mái dốc hồ đập.
  • Các nhà Hoạch định Chính sách và Ban Chỉ huy Phòng chống Thiên tai: Sử dụng bản đồ phân vùng 5 cấp độ tai biến để xây dựng phương án sơ tán dân cư, phân bổ nguồn lực ứng phó khẩn cấp và cấp phép xây dựng công trình dân sinh an toàn.

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 thuyết Sức kháng cắt đất không bão hòa của Fredlund & Rahardjo (1993) kết hợp lý thuyết thấm không dừng Richards 2D vào việc giải thích tường minh cơ chế sụt giảm độ an toàn mái dốc tàn tích nhiệt đới do mưa. Luận án đã lượng hóa chính xác sự triệt tiêu của độ hút dính $(u_a - u_w)$ dẫn đến suy giảm sức kháng cắt thực tế dọc theo mặt trượt tiềm năng, chứng minh bản chất của trượt lở do mưa nhiệt đới là quá trình suy thoái sức bền kháng cắt do bão hòa ẩm chứ không thuần túy là quá trình tăng tải trọng bản thân.

2. Đổi mới phương pháp luận của luận án thể hiện ra sao khi so sánh với các nghiên cứu trước đây?

So sánh với nghiên cứu của Yilmaz (2009) tại Thổ Nhĩ Kỳ và Bùi Tiến Diệu et al. (2011) tại Hòa Bình, đổi mới phương pháp luận của luận án nằm ở 3 điểm:

  1. Chuyển đổi từ bản đồ nguy cơ tĩnh (LSI) sang bản đồ tai biến động (LHI) có tích hợp lượng mưa lớn nhất năm (MP) và chỉ số ẩm ướt địa hình TWI mở rộng.
  2. Lần đầu tiên ứng dụng thuật toán trung bình mô hình Bayesian (BMA) trên R để giải quyết bài toán lựa chọn biến tối ưu, loại bỏ yếu tố chủ quan của con người.
  3. Thiết lập quy trình chuẩn hóa giá trị chỉ số LHI về phân bố chuẩn và phân vùng bằng độ lệch chuẩn (Standard Deviation), giải quyết triệt để vấn đề so sánh tương thích giữa các mô hình FR, SI, WoE, LR.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là: Trong các đợt mưa kéo dài ($> 3\text{ ngày}$), đất có hệ số thấm nhỏ ($k_s \approx 2,09 \times 10^{-7}\text{ m/s}$, đất sét) lại chịu sự sụt giảm hệ số an toàn $\text{FS}$ nghiêm trọng và nguy cơ trượt lở cao hơn nhiều so với đất có hệ số thấm lớn ($k_s \approx 2,03 \times 10^{-6}\text{ m/s}$, đất sét pha lẫn dăm sạn). Cơ chế là do đất thấm bé cản trở sự thoát nước trọng lực, hình thành tầng ứ đọng nước bão hòa cục bộ sát bề mặt mái dốc, gây phá hủy trượt nông tức thì.

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

Hoàn toàn có. Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình: từ bảng chỉ tiêu cơ lý đất thí nghiệm (Bảng 2.2), tọa độ 231 điểm trượt lở, phương pháp nội suy 11 lớp chuyên đề trên GIS, hệ số phương trình hồi quy LR, ma trận trọng số FR/SI/WoE (Bảng 4.4, 4.5, 4.6), mã nguồn cấu hình BMA trên R, cho đến các thông số đầu vào trong mô đun SEEP/W và SLOPE/W (Geostudio 2007).

5. Kế hoạch nghiên cứu 10 năm được vạch ra từ nền tảng luận án này là gì?

Kế hoạch 10 năm định hướng: Xây dựng hệ thống tự động hóa tích hợp mô hình thấm 3D thời gian thực với mạng lưới quan trắc mưa vệ tinh/radar và trạm cảm biến IoT áp lực nước lỗ rỗng ngầm; phát triển nền tảng WebGIS trí tuệ nhân tạo cảnh báo sớm nguy cơ tai biến trượt lở theo thời gian thực (Real-time Landslide Early Warning WebGIS Platform) cho toàn bộ 14 tỉnh thành ven biển miền Trung và Tây Nguyên.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và khoa học mối quan hệ hữu cơ giữa quá trình thấm nước mưa không bão hòa với sự suy giảm sức chống cắt và độ ổn định mái dốc đất tàn tích phong hóa huyện Khánh Vĩnh, tỉnh Khánh Hòa.
  2. Xác lập chính xác luận điểm khoa học số 1: Tai biến trượt lở mái dốc đất tàn tích ($k_s = 10^{-7} - 10^{-5}\text{ m/s}$) bắt đầu bị kích hoạt phá hủy khi đồng thời hội đủ các điều kiện: Chiều cao $H \ge 10\text{ m}$, góc dốc $\alpha \ge 45^\circ$, cường độ mưa đạt ngưỡng $I \ge 10\text{ mm/h}$ và thời gian mưa kéo dài liên tục từ 5 đến 7 ngày.
  3. Xác lập luận điểm khoa học số 2: Các mô hình xác suất thống kê tích hợp GIS đều đạt độ tin cậy và độ chính xác xuất sắc ($\text{AUC} = 0,8 - 0,9$); trong đó phương pháp Trọng số chứng cứ (WoE) cho kết quả dự báo chính xác nhất; mô hình 4 yếu tố cốt lõi (DEM, DFR, MP, SLOPE) do BMA xác lập là mô hình tối ưu nhất về hiệu quả dự báo.
  4. Mở ra bước chuyển đổi nhận thức khoa học căn bản từ đánh giá bản đồ nguy cơ trượt lở tĩnh (LSI) sang bản đồ tai biến trượt lở động (LHI) thông qua việc kết hợp yếu tố kích hoạt lượng mưa lớn nhất năm (MP) và chỉ số ẩm ướt địa hình TWI mở rộng.
  5. Mở ra 3 nhánh nghiên cứu chuyên sâu mới: Mô hình hóa thấm 3D tương tác nứt nẻ; Ứng dụng BMA kết hợp Machine Learning trong địa kỹ thuật môi trường; Thiết lập hệ thống cảnh báo sớm trượt lở thông minh dựa trên dữ liệu khí tượng viễn thám thời gian thực.
  6. Cung cấp sản phẩm bản đồ phân vùng tai biến trượt lở 5 cấp độ tin cậy cao, đóng vai trò là tài liệu khoa học nền tảng phục vụ trực tiếp công tác quy hoạch phát triển kinh tế - xã hội, bảo vệ an toàn các tuyến giao thông trọng điểm và phòng chống giảm nhẹ thiên tai bền vững tại khu vực Nam Trung Bộ.