Tổng quan về luận án

Phân tích ổn định mái dốc là một trong những bài toán cốt lõi của kỹ thuật địa kỹ thuật công trình thủy lợi, giao thông và xây dựng. Tại Việt Nam, với hơn 74,3% diện tích tự nhiên là địa hình đồi núi và trung du cùng mạng lưới hơn 6.000 hồ đập và hàng chục nghìn kilômét đê điều, hiện tượng sạt trượt mái dốc dưới tác động của lượng mưa nhiệt đới gió mùa diễn ra đặc biệt phức tạp. Các thảm họa trượt lở nghiêm trọng tại đập Đầm Hà Động (Quảng Ninh), đập Đầm Thìn (Phú Thọ), đê tả sông Chu (Thanh Hóa) hay khu vực miền Trung vào tháng 10/2020 đã chứng minh tính cấp bách của việc giải quyết triệt để bài toán ổn định công trình đất đắp. Luận án tiến sĩ kỹ thuật "Nghiên cứu ảnh hưởng của mưa đến ổn định mái dốc đất không bão hòa" của nghiên cứu sinh Phạm Huy Dũng (ngành Địa kỹ thuật xây dựng, mã số 9580211, Trường Đại học Thủy lợi, năm 2021) đã tiên phong chuyển đổi mô hình tính toán từ cơ học đất bão hòa truyền thống sang cơ học đất không bão hòa tiên tiến.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án nhắm tới là sự thiếu hụt dữ liệu thực nghiệm về đường cong đặc trưng đất - nước (Soil-Water Characteristic Curve - SWCC), cơ chế suy giảm lực hút dính (matric suction) và mối tương quan phi tuyến giữa lượng mưa xâm nhập với sự suy thoái cường độ kháng cắt của đất đắp tại miền Bắc Việt Nam. Các mô hình thiết kế đê đập hiện hành thường mặc định đất ở trạng thái bão hòa hoàn toàn hoặc khô hoàn toàn, bỏ qua đóng góp của áp lực nước lỗ rỗng âm. Luận án đặt ra hệ thống câu hỏi và giả thuyết nghiên cứu có cấu trúc chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố và suy giảm lực hút dính $(u_a - u_w)$ trong thân đê, đập diễn ra như thế nào dưới các chế độ mưa có cường độ và thời gian khác nhau?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ ảnh hưởng của hệ số đầm chặt $(K)$ và độ dốc mái $(m)$ đến cường độ nước chảy tràn bề mặt và lượng ẩm thấm vào khối đất không bão hòa là bao nhiêu?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ chính xác hóa hệ số an toàn ổn định mái dốc $(F_s)$ khi tích hợp lý thuyết cường độ kháng cắt mở rộng so với phương pháp cổ điển?
  • Giả thuyết nghiên cứu 1 (H1): Lực hút dính đóng góp một thành phần ma sát biểu kiến $\phi^b$ làm tăng đáng kể cường độ kháng cắt của đất, nhưng giá trị này suy giảm phi tuyến theo thời gian mưa kéo dài.
  • Giả thuyết nghiên cứu 2 (H2): Trận mưa kéo dài cường độ nhỏ (LD) gây suy giảm hệ số ổn định nghiêm trọng và tạo mặt trượt sâu hơn so với trận mưa thời gian ngắn cường độ lớn (HI).
  • Giả thuyết nghiên cứu 3 (H3): Mô hình ước lượng Modified Kovacs (MK) hiệu chỉnh cho phép dự báo chính xác SWCC của đất dính miền Bắc mà không cần thí nghiệm áp lực cao phức tạp.

Khung lý thuyết của công trình dựa trên lý thuyết cơ học môi trường liên tục đa pha của Fredlund & Morgenstern (1977) và mô hình cường độ kháng cắt không bão hòa của Fredlund, Morgenstern & Widger (1978). Phạm vi nghiên cứu bao quát các công trình trọng điểm: đê hữu Cầu (Bắc Ninh), đập Khau Piều (Lạng Sơn), đập Chúc Bài Sơn (Quảng Ninh) và đất đắp mỏ Đại Phong, kết hợp mô hình vật lý máng nghiêng kích thước lớn và mô phỏng số 2D Geostudio.

Literature Review và Positioning

Lịch sử phát triển của bài toán ổn định mái dốc trải qua ba giai đoạn chính: phương pháp cân bằng giới hạn (Limit Equilibrium Methods - LEM), phương pháp phân tích giới hạn (Limit Analysis) và phương pháp số (Finite Element Method - FEM). Trường phái cân bằng giới hạn phân mảnh tĩnh học phát triển từ các mô hình giả định đơn giản không xét lực tương tác giữa các thỏi như Ordinary/Fellenius (1936), mô hình chỉ thỏa mãn cân bằng mô-men như Bishop đơn giản (1955), mô hình thỏa mãn cân bằng lực như Janbu đơn giản (1954), đến các phương pháp chặt chẽ thỏa mãn đồng thời cân bằng lực và mô-men như Spencer (1967), Morgenstern-Price (1965) và Sarma (1973). Tuy nhiên, các phương pháp này ban đầu chỉ giới hạn trong khuôn khổ cơ học đất bão hòa với nguyên lý ứng suất hiệu quả của Terzaghi (1936):

$$\sigma' = \sigma - u_w$$

Trong cơ học đất không bão hòa, cuộc tranh luận học thuật kéo dài nhiều thập kỷ tập trung vào biến trạng thái ứng suất. Bishop (1959) từng đề xuất phương trình ứng suất hiệu quả duy nhất chứa tham số $\chi$:

$$\sigma' = (\sigma - u_a) + \chi(u_a - u_w)$$

Tuy nhiên, các bằng chứng thực nghiệm của Bishop & Blight (1963), Burland (1965) và Fredlund & Morgenstern (1977) đã bác bỏ tính duy nhất của $\chi$, chứng minh sự cần thiết phải sử dụng hai biến trạng thái ứng suất độc lập: ứng suất pháp thực $(\sigma - u_a)$ và lực hút dính $(u_a - u_w)$.

Về mặt đường cong đặc trưng đất - nước (SWCC), các nghiên cứu kinh điển của Brooks & Corey (1964), van Genuchten (1980) và Fredlund & Xing (1994) đã mô tả mối quan hệ giữa hàm lượng nước và lực hút dính thông qua các thông số hình học: giá trị khí vào (Air Entry Value - AEV), độ dốc vùng chuyển tiếp và độ ẩm dư $(\theta_r)$. Đối với cường độ kháng cắt, trong khi Fredlund et al. (1978) giả thiết góc ma sát biểu kiến $\phi^b$ là hằng số tuyến tính, thì các nghiên cứu của Escario & Saez (1986) trên đất sét Madrid và Satija (1978) trên đất sét Dhanauri đã chỉ ra bản chất phi tuyến rõ rệt của $\phi^b$ khi lực hút dính vượt qua ngưỡng AEV.

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

  • Nghiên cứu của Gasmo, Rahardjo & Fredlund (2000): Khi phân tích mái dốc đất tàn tích tại Singapore, tác giả nhận thấy "Trận mưa với cường độ 50mm/ngày sẽ làm giảm hệ số an toàn mái dốc khoảng 25% sau thời gian 12 giờ, trong khi đó phải mất hai tuần lễ sau khi mưa thì hệ số an toàn mới tăng ngược trở lại 30%".
  • Nghiên cứu của Rahardjo et al. (2002): Khảo sát mái dốc tàn tích chịu mưa liên tục 24 giờ với cường độ 30 mm/h chỉ ra rằng sự dâng cao của áp lực nước lỗ rỗng làm chiều sâu cung trượt nguy hiểm bị nông hóa rõ rệt, bất kể chiều cao hình học của mái dốc.

Luận án của Phạm Huy Dũng đã định vị xuất sắc trong dòng chảy học thuật khi giải quyết bài toán thiếu hụt dữ liệu SWCC tại Việt Nam (vốn trước đó chỉ có các công bố sơ khai của Hằng et al., 2010; Hương, 2013; Toàn, 2016), đồng thời cung cấp bằng chứng thực nghiệm từ mô hình máng tạo mưa nhân tạo quy mô lớn mà các nghiên cứu trước đây tại Việt Nam chưa thực hiện được.

Đó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à kiểm chứng lý thuyết cường độ kháng cắt Mohr-Coulomb mở rộng cho các loại đất dính đầm chặt vùng nhiệt đới:

$$\tau_f = c' + (\sigma - u_a)\tan\phi' + (u_a - u_w)\tan\phi^b$$

Trong đó:

  • $\tau_f$: Cường độ kháng cắt của đất không bão hòa (kPa).
  • $c'$: Lực dính hiệu quả (kPa).
  • $(\sigma - u_a)$: Ứng suất pháp thực (kPa).
  • $(u_a - u_w)$: Lực hút dính hay áp lực nước lỗ rỗng âm (kPa).
  • $\phi'$: Góc ma sát trong hiệu quả (độ).
  • $\phi^b$: Góc ma sát biểu kiến, phản ánh tốc độ gia tăng sức chống cắt theo lực hút dính (độ).

Luận án đã chứng minh sự chuyển dịch mô hình (paradigm shift) từ hàm tuyến tính sang hàm phi tuyến của $\phi^b$:

  • Khi lực hút dính $(u_a - u_w) \le \text{AEV}$, góc ma sát biểu kiến đạt cực đại xấp xỉ bằng góc ma sát trong: $\phi^b \approx \phi'$.
  • Khi $(u_a - u_w) > \text{AEV}$ và tăng lên vùng $200 - 300\text{ kPa}$, $\phi^b$ suy giảm nhanh chóng và tiệm cận về khoảng $(1/3 - 1/2)\phi'$.
                [Độ bão hòa hoàn toàn S = 100%]
                               │
               (Thấm & Thoát nước mặt)
                               ▼
        ┌──────────────────────────────────────────────┐
        │     Hệ thống Đất Không Bão Hòa (4 pha)       │
        │ - Hạt khoáng  - Nước  - Khí  - Màng co mỏng  │
        └──────────────────────┬───────────────────────┘
                               │
            ┌──────────────────┴──────────────────┐
            ▼                                     ▼
┌──────────────────────────────┐     ┌──────────────────────────────┐
│       Mô hình Ứng suất       │     │     Thủy lực & Thấm dột      │
│  - Ứng suất pháp: (σ - ua)   │     │  - Định luật Darcy biến tính │
│  - Lực hút dính: (ua - uw)   │     │  - Đường cong SWCC (Hàm MK)  │
└──────────────┬───────────────┘     └──────────────┬───────────────┘
               │                                    │
               └──────────────────┬─────────────────┘
                                  ▼
             ┌─────────────────────────────────────────┐
             │       Phương trình Sức kháng cắt        │
             │   τf = c' + (σ-ua)tanϕ' + (ua-uw)tanϕb  │
             └────────────────────┬────────────────────┘
                                  ▼
             ┌─────────────────────────────────────────┐
             │     Phân tích Cân bằng giới hạn GLE     │
             │        (Cân bằng Lực & Mô-men)          │
             │   => Đánh giá Hệ số An toàn Fs(t)       │
             └─────────────────────────────────────────┘

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

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

  1. Lý thuyết thấm không bão hòa: Mở rộng định luật Darcy với hệ số thấm là một hàm phi tuyến biến thiên theo lực hút dính $k_w = k_w(u_a - u_w)$ suy dẫn từ SWCC qua mô hình Fredlund, Xing & Huang (1994).
  2. Lý thuyết sức kháng cắt đa biến: Tách biệt ảnh hưởng của khung hạt khoáng và sức căng bề mặt của màng co khí - nước (contractile skin).
  3. Phương pháp cân bằng giới hạn tổng quát (General Limit Equilibrium - GLE): Thiết lập phương trình vi phân cân bằng đồng thời mô-men ($F_m$) và cân bằng lực ($F_f$) với hàm lực tương tác giữa các thỏi $f(x)$ biến thiên.

Điều kiện biên và phạm vi áp dụng (boundary conditions) được xác định rõ: Đất dính đầm chặt có độ đầm chặt $K = 0,90 - 0,97$, độ dốc mái $m = 1:1 - 1:4$, chịu các hình thái mưa phân bố theo thời gian thực tế tại miền Bắc Việt Nam.

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ủ thế giới quan thực chứng (positivism) kết hợp chủ nghĩa hiện thực phản biện (critical realism), áp dụng phương pháp hỗn hợp định lượng từ thực nghiệm vi mô trong phòng, mô hình vật lý tỷ lệ lớn đến mô phỏng số hiện trường.

Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm:

  • Tầng 1 (Vật liệu & Cơ lý vi mô): Phân tích thành phần cấp phối hạt, giới hạn Atterberg (giới hạn chảy $w_L$, giới hạn dẻo $w_P$, chỉ số dẻo $I_P$), khối lượng thể tích khô $\rho_d$, độ ẩm tối ưu $w_{opt}$.
  • Tầng 2 (Đặc trưng không bão hòa): Đo lực hút dính hiện trường, xác lập SWCC và hàm thấm, kiểm tra sức kháng cắt ba trục và cắt trực tiếp không bão hòa.
  • Tầng 3 (Mô hình vật lý tương đồng): Máng thí nghiệm biến đổi góc nghiêng có dàn tạo mưa kiểm soát lưu lượng.
  • Tầng 4 (Mô hình số hóa công trình thực tế): Phân tích trường thấm phi dừng (SEEP/W) kết hợp ổn định mái dốc (SLOPE/W).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được thiết kế đạt độ tin cậy và giá trị khoa học cao:

  • Đo lực hút dính hiện trường: Sử dụng thiết bị căng kế chuyên dụng kiểu Jet-fill 2725ARL (Soilmoisture Equipment Corp) với ống đo áp có khả năng làm đầy bọt khí nhanh, đo chính xác lực hút dính lên tới 80 - 100 kPa. Kết nối hệ thống cảm biến Remote Tensiometer 2100F và bộ chuyển đổi dữ liệu thông minh ICT Data Logger ghi nhận dữ liệu liên tục theo thời gian thực.
  • Xác định SWCC trong phòng: Sử dụng bình chiết áp lực cao (Pressure Plate Extractor) trang bị đĩa tiếp nhận khí cao (High Air Entry Ceramic Disc) chịu áp lực 15 bar (1.500 kPa). Mẫu đất được bão hòa hoàn toàn trong chân không trước khi gia tải từng cấp áp lực khí nén $u_a$ để đuổi nước lỗ rỗng, xác định chính xác độ ẩm thể tích $\theta$ tại từng cấp lực hút dính $(u_a - u_w)$.
  • Thí nghiệm sức kháng cắt: Thực hiện trên hệ thống máy nén ba trục TRIAX50 và máy cắt phẳng trực tiếp tự động EDI-2 với các sơ đồ cố kết thoát nước (CD) và độ ẩm không đổi (CW).
┌────────────────────────────────────────────────────────────────────────┐
│                        QUY TRÌNH THỰC NGHIỆM RIGOROUS                  │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
    ┌───────────────────────────────┼───────────────────────────────┐
    ▼                               ▼                               ▼
┌───────────────────────┐ ┌───────────────────────┐ ┌───────────────────────┐
│  Thí nghiệm Căng kế   │ │ Bình chiết Áp lực   │ │  Máy nén Ba trục    │
│  Jet-fill 2725ARL /   │ │ High-Air Entry      │ │  TRIAX50 & EDI-2    │
│  Remote Tens. 2100F   │ │ Ceramic Disc 15 bar │ │  (Sơ đồ CD & CW)    │
└───────────┬───────────┘ └───────────┬───────────┘ └───────────┬───────────┘
            │                         │                         │
            ▼                         ▼                         ▼
┌───────────────────────┐ ┌───────────────────────┐ ┌───────────────────────┐
│ Đo lực hút dính hiện  │ │ Xây dựng SWCC &       │ │ Xác định c', ϕ' và    │
│ trường (0 - 100 kPa)  │ │ Ước lượng mô hình MK  │ │ hàm góc ma sát ϕb    │
└───────────┬───────────┘ └───────────┬───────────┘ └───────────┬───────────┘
            │                         │                         │
            └─────────────────────────┼─────────────────────────┘
                                      ▼
                        ┌───────────────────────────┐
                        │   Máng mô hình Vật lý     │
                        │ (Dàn mưa + Cảm biến ICT)  │
                        └─────────────┬─────────────┘
                                      ▼
                        ┌───────────────────────────┐
                        │    Mô phỏng Số hóa        │
                        │ (SEEP/W & SLOPE/W - 2D)   │
                        └───────────────────────────┘

Data và phân tích

Đặc tính mẫu đất nghiên cứu tại các công trình thực tế được tổng hợp nghiêm ngặt:

  • Đất đắp đê hữu Cầu: Đất á sét, dung trọng tự nhiên $\gamma = 18,2\text{ kN/m}^3$, lực dính hiệu quả $c' = 16,5\text{ kPa}$, góc ma sát trong $\phi' = 19,5^\circ$. Lực hút dính đo tại đỉnh đê đạt $18,4 - 42,0\text{ kPa}$ (đợt 1) và $25,1 - 56,2\text{ kPa}$ (đợt 2).
  • Đất đắp đập Khau Piều: Đất sét pha dẻo cao, $c' = 24,0\text{ kPa}$, $\phi' = 16,0^\circ$, độ ẩm tối ưu $w_{opt} = 17,2%$, dung trọng khô lớn nhất $\rho_{d,max} = 1,72\text{ g/cm}^3$.
  • Đất đắp đập Chúc Bài Sơn: Đất á sét nặng, $c' = 28,5\text{ kPa}$, $\phi' = 18,2^\circ$, AEV xác định từ SWCC đạt xấp xỉ $32\text{ kPa}$.
  • Vật liệu đắp mỏ Đại Phong (thí nghiệm máng tạo mưa): Khảo sát tại 3 cấp đầm chặt $K = 0,90; 0,95; 0,97$ tương ứng dung trọng khô $\rho_d = 1,51; 1,60; 1,63\text{ g/cm}^3$.

Công cụ phân tích: Phần mềm GeoStudio (Canada) với mô-đun SEEP/W (giải phương trình vi phân dòng thấm 2D không bão hòa bằng phần tử hữu hạn) và SLOPE/W (phân tích ổn định bằng phương pháp Morgenstern-Price và Spencer).

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

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

  1. Quy luật phi tuyến của góc ma sát biểu kiến $\phi^b$: Thí nghiệm cắt trực tiếp và ba trục khẳng định khi lực hút dính $(u_a - u_w) < 40\text{ kPa}$ (vùng trước AEV), $\phi^b \approx \phi'$. Khi $(u_a - u_w)$ tăng từ $50\text{ kPa}$ lên $200 - 300\text{ kPa}$, $\phi^b$ giảm mạnh xuống còn $7^\circ - 10^\circ$, tương đương $(1/3 - 1/2)\phi'$. Điều này bác bỏ giả định góc $\phi^b$ không đổi của các tiêu chuẩn cũ.
  2. Vai trò khống chế của hệ số đầm chặt ($K$) đối với tỷ lệ thấm và dòng chảy tràn: Trong mô hình máng tạo mưa, với đất có độ đầm chặt cao $K = 0,97$, cường độ nước chảy tràn bề mặt đạt mức ổn định $i_t = 0,148\text{ mm/giây}$, tỷ lệ nước mưa xâm nhập chỉ chiếm $26,2%$. Ngược lại, khi độ đầm chặt thấp $K = 0,90$, cường độ tràn giảm xuống $i_t = 0,085\text{ mm/giây}$, lượng mưa thâm nhập vào thân mái dốc tăng vọt lên tới $57,3%$, làm suy giảm nhanh chóng lực hút dính nội tại.
  3. Ảnh hưởng dị biệt của các hình thái mưa (Rainfall Patterns): Kiểu mưa kéo dài cường độ nhỏ (LD - Long Duration) nguy hiểm hơn nhiều so với kiểu mưa ngắn cường độ lớn (HI - High Intensity). Kiểu mưa LD cho phép toàn bộ lượng nước mưa thâm nhập sâu vào khối đắp mà không tạo dòng tràn, làm triệt tiêu hoàn toàn lực hút dính bề mặt, đẩy mặt trượt nguy hiểm ăn sâu vào thân đập và làm giảm hệ số an toàn $F_s$ từ $25%$ đến $38%$.
  4. Hiệu ứng trễ phục hồi lực hút dính sau khi dứt mưa: Dữ liệu cảm biến ICT ghi nhận sau khi kết thúc mưa 3 ngày, lực hút dính tại đới sát mặt dốc chỉ phục hồi chưa đầy $15 - 20%$ giá trị ban đầu. Sự duy trì áp lực nước lỗ rỗng cao trong thời gian dài lý giải hiện tượng trượt lở đất thường xảy ra sau khi mưa đã tạnh nhiều ngày.
Trường hợp công trình Trạng thái bão hòa truyền thống ($F_s$) Trạng thái không bão hòa trước mưa ($F_s$) Sau mưa lớn HI ($F_s$) Sau mưa kéo dài LD ($F_s$) Mức suy giảm $F_s$ do mưa LD (%)
Đê hữu Cầu 1,185 1,462 1,280 1,195 -18,26%
Đập Khau Piều 1,220 1,584 1,342 1,210 -23,61%
Đập Chúc Bài Sơn 1,310 1,695 1,415 1,265 -25,37% 
Hệ số Ổn định (Fs)
  1.8 ┼─────────────────────────────────────────────────────────────
      │ [1.695] Đập Chúc Bài Sơn (KBH ban đầu)
  1.6 ┼───▲─────────────────[1.584] Đập Khau Piều (KBH ban đầu)
      │   │                      ▲
  1.4 ┼───┼──────────────────────┼──────────────────[1.415] Sau mưa HI (Chúc Bài Sơn)
      │   │                      │                     ▲
  1.2 ┼───┼──────────────────────┼─────────────────────┼────────[1.210] Sau mưa LD (Khau Piều)
      │   ▼                      ▼                     ▼            ▼
  1.0 ┼─[1.185] Đê Cầu (BH)  [1.220] Khau Piều (BH) [1.310] CBS(BH) [1.195] Sau mưa LD (Đê Cầu)
      └─────────────────────────────────────────────────────────────

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ tham số thực nghiệm hoàn chỉnh về SWCC và góc ma sát $\phi^b$ cho đất dính miền Bắc Việt Nam, bổ sung vào cơ sở dữ liệu Địa kỹ thuật quốc tế cho vùng nhiệt đới gió mùa.
  • Về phương pháp luận: Chuẩn hóa quy trình kết hợp đo đạc lực hút dính tự động bằng căng kế điện tử với thí nghiệm mô hình máng tạo mưa có thể điều chỉnh độ dốc.
  • Về ứng dụng thực tiễn: Kiến nghị các cơ quan tư vấn thiết kế công trình thủy lợi (như HEC, VEC) chuyển đổi thuật toán kiểm tra ổn định, bổ sung trường hợp tính toán ứng với tổ hợp mưa bất lợi kéo dài nhiều ngày (LD).
  • Về chính sách phòng chống thiên tai: Cung cấp cơ sở khoa học để Tổng cục Phòng chống Thiên tai và Bộ Nông nghiệp & PTNT xây dựng ngưỡng cảnh báo trượt lở đê đập dựa trên lượng mưa tích lũy 10 ngày ($P_{10}$) kết hợp kiểm soát độ đầm chặt đắp thân đê ($K \ge 0,95$).

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 giới hạn nghiên cứu:

  1. Hiệu ứng tỷ lệ mô hình vật lý: Thí nghiệm máng tạo mưa trong phòng dù có kích thước lớn nhưng chưa phản ánh đầy đủ hiệu ứng ứng suất địa tầng thực tế của các đập đất cao trên 30 m.
  2. Hiện tượng tạo bọt khí trong căng kế: Thiết bị căng kế Jet-fill chỉ đo chính xác lực hút dính dưới 90 kPa; khi đất quá khô, hiện tượng xâm thực khí (cavitation) xảy ra làm gián đoạn chuỗi đo liên tục.
  3. Mô hình dòng thấm 2D: Phân tích ổn định trong GeoStudio mới chỉ dừng lại ở bài toán biến dạng phẳng 2D, chưa xét đến hiệu ứng không gian 3D tại các khúc quanh đê hoặc vai đập tiếp giáp địa hình phức tạp.
  4. Bỏ qua ảnh hưởng của thảm phủ thực vật: Chưa định lượng tác động của hệ rễ cây trong việc gia cường cơ học và hút nước làm tăng lực hút dính sinh học (bio-suction).

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 nghiên cứu mô phỏng liên kết 3D Thủy - Cơ học (Coupled Hydro-Mechanical 3D Modeling).
  • Ứng dụng công nghệ mạng cảm biến không dây IoT đo áp lực nước lỗ rỗng âm thời gian thực tích hợp hệ thống AI cảnh báo sớm rủi ro vỡ đê đập.
  • Nghiên cứu cơ chế bio-suction của các loài cỏ bảo vệ mái dốc (như cỏ Vetiver) trong việc hạn chế dòng thấm nhiệt đới.

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

Công trình của Phạm Huy Dũng tạo ra tác động lan tỏa mạnh mẽ trong cả giới hàn lâm lẫn ngành công nghiệp xây dựng thủy lợi. Về mặt học thuật, các công bố trích xuất từ luận án trên các tạp chí chuyên ngành uy tín đã trở thành tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh và học viên cao học ngành Địa kỹ thuật tại Việt Nam.

Về mặt công nghiệp và chuyển giao kỹ thuật, kết quả nghiên cứu đã được trực tiếp ứng dụng trong các dự án nâng cấp đê điều và gia cố hồ chứa trọng điểm tại các tỉnh Bắc Ninh, Lạng Sơn và Quảng Ninh. Việc tối ưu hóa hệ số đầm chặt $K$ và cấu tạo tầng lọc tiêu thoát nước mặt đã giúp tiết kiệm hàng tỷ đồng chi phí xử lý sạt trượt cục bộ, đồng thời bảo vệ an toàn tính mạng cho hàng trăm nghìn người dân sinh sống tại vùng hạ lưu các đập thủy lợi trong mùa mưa bão.

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

  • Nghiên cứu sinh & Nhà khoa học Địa kỹ thuật: Tiếp cận hệ phương pháp thực nghiệm SWCC tiên tiến, khung lý thuyết mở rộng Mohr-Coulomb và cơ sở dữ liệu đất không bão hòa chuẩn hóa.
  • Kỹ sư tư vấn thiết kế đê đập: Nắm vững công cụ tính toán chính xác hóa hệ số an toàn mái dốc thông qua phần mềm SEEP/W và SLOPE/W, tránh tình trạng thiết kế thừa lãng phí hoặc thiếu an toàn cục bộ.
  • Doanh nghiệp thi công xây dựng: Nhận thức rõ tầm quan trọng của việc kiểm soát hệ số đầm chặt ($K \ge 0,95$) trong việc ngăn ngừa nước mưa thâm nhập thân công trình.
  • Nhà quản lý & Cơ quan hoạch định chính sách: Sử dụng các phát hiện về kiểu mưa nguy hiểm (LD) và thời gian trễ phục hồi lực hút dính để xây dựng kịch bản ứng phó, vận hành xả lũ và cảnh báo thiên tai chuẩn xác.

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ì? Đó là việc chứng minh và định lượng hóa bản chất phi tuyến của góc ma sát biểu kiến $\phi^b$ trong phương trình sức chống cắt của Fredlund et al. (1978) đối với đất dính đầm chặt miền Bắc. Luận án chỉ rõ $\phi^b$ không phải là hằng số mà suy giảm mạnh từ $\phi^b \approx \phi'$ xuống $(1/3 - 1/2)\phi'$ khi lực hút dính vượt qua giá trị khí vào AEV ($30 - 50\text{ kPa}$).

2. Đổi mới phương pháp luận nổi bật so với các nghiên cứu trước đây? So sánh với các nghiên cứu của Gasmo et al. (2000) và Rahardjo et al. (2002) chỉ dựa thuần túy vào mô phỏng số, luận án đã tự thiết kế, chế tạo thành công hệ thống máng thí nghiệm góc nghiêng tự do kết hợp dàn mưa nhân tạo kiểm soát lưu lượng và tích hợp cảm biến ICT thu thập dữ liệu lực hút dính tự động theo thời gian thực.

3. Phát hiện nào gây bất ngờ nhất về mặt dữ liệu thực nghiệm? Phát hiện về sự vượt trội của kiểu mưa kéo dài cường độ nhỏ (LD) trong việc phá hủy mái dốc so với mưa lớn thời gian ngắn (HI). Kiểu mưa LD có tỷ lệ nước mưa xâm nhập đạt tới trên $57%$, triệt tiêu hoàn toàn đới lực hút dính sâu trong thân đập, gây sụt giảm hệ số an toàn $F_s$ lên tới $38%$.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) không? Luận án cung cấp chi tiết toàn bộ thông số hình học máng thí nghiệm, sơ đồ mạch cảm biến căng kế 2100F, quy trình bão hòa mẫu đất và các bước hiệu chuẩn áp lực khí trong bình chiết 15 bar, cho phép tái lập quy trình thực nghiệm với độ chính xác cao.

5. Định hướng lộ trình nghiên cứu 10 năm tiếp theo? Lộ trình 10 năm tập trung vào bài toán tương tác 3D Thủy - Cơ - Nhiệt (THM), tích hợp cảm biến đo hút dính nano không dây trong đê biển thông minh và tích hợp dữ liệu viễn thám vệ tinh InSAR để dự báo sạt trượt thời gian thực.

Kết luận

Luận án tiến sĩ của Phạm Huy Dũng đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Xác lập hệ cơ sở dữ liệu thực nghiệm toàn diện đầu tiên về SWCC, lực hút dính và sức kháng cắt không bão hòa cho các loại đất dính công trình thủy lợi miền Bắc Việt Nam.
  2. Kiểm chứng và hiệu chỉnh thành công mô hình ước lượng Modified Kovacs (MK) và phương trình Fredlund & Xing (1994), tạo công cụ dự báo SWCC tin cậy khi không có điều kiện thí nghiệm áp lực cao.
  3. Chế tạo thành công hệ thống thiết bị thực nghiệm máng nghiêng tạo mưa hiện đại, giải mã cơ chế thâm nhập của nước mưa và quy luật chảy tràn phụ thuộc vào hệ số đầm chặt ($K$) và độ dốc mái ($m$).
  4. Làm sáng tỏ quy luật biến thiên áp lực nước lỗ rỗng âm và sự suy giảm hệ số ổn định dưới các kịch bản mưa khác nhau, khẳng định tính chất nguy hiểm cực hạn của dạng mưa kéo dài LD.
  5. Ứng dụng thành công lý thuyết vào tính toán thực tế cho đê hữu Cầu, đập Khau Piều và đập Chúc Bài Sơn, chính xác hóa hệ số an toàn và nâng cao độ tin cậy trong thiết kế công trình đất đắp.

Công trình tạo ra bước tiến quan trọng trong việc chuyển đổi tư duy kỹ thuật địa kỹ thuật tại Việt Nam, mở ra các hướng nghiên cứu mới về công trình thủy lợi thông minh thích ứng với biến đổi khí hậu toàn cầu.