Tổng quan về luận án

Hiện tượng mất ổn định và trượt lở mái dốc đất dưới tác động của mưa là một trong những hiểm họa địa kỹ thuật nghiêm trọng nhất, đe dọa trực tiếp đến an toàn tính mạng con người và hệ thống cơ sở hạ tầng giao thông, thủy lợi. Tại Việt Nam, với đặc thù hơn 74,3% diện tích tự nhiên là đồi núi và vùng trung du cùng hàng chục nghìn kilomet đê điều, hồ đập phân bố rộng khắp, rủi ro sạt trượt mái dốc trong mùa mưa bão diễn ra vô cùng phức tạp. Về mặt học thuật truyền thống, các phân tích ổn định mái dốc cổ điển thường dựa trên tiền đề của cơ học đất bão hòa (Terzaghi, 1936), trong đó đất nằm trên mực nước ngầm mặc định được xem là khô hoàn toàn hoặc bão hòa cục bộ khi có dòng thấm. Cách tiếp cận giản lược hóa này đã hoàn toàn bỏ qua vai trò cốt lõi của lực hút dính ma-trích (matric suction, $u_a - u_w$) và áp lực nước lỗ rỗng âm trong việc gia tăng sức kháng cắt của khối đất.

Khoảng trống nghiên cứu (research gap) then chốt mà 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 NCS. Phạm Huy Dũng (Trường Đại học Thủy lợi, chuyên ngành Địa kỹ thuật xây dựng - Mã số: 9580211, dưới sự hướng dẫn khoa học của PGS.TS Hoàng Việt Hùng và GS. Nguyễn Công Mẫn) tập trung giải quyết gồm hai vấn đề mang tính quyết định: (1) Sự khan hiếm trầm trọng của cơ sở dữ liệu thực nghiệm về đường cong đặc trưng đất - nước (SWCC), hàm thấm không bão hòa ($k_w$) và tham số cường độ kháng cắt mở rộng ($\phi^b$) đối với các loại đất đắp phân bố tại miền Bắc Việt Nam; (2) Sự thiếu vắng các mô hình thực nghiệm mẫu lớn (flume tests) có khả năng định lượng trực tiếp cơ chế phân chia dòng chảy (tỷ lệ thấm $R_I$ so với tỷ lệ tràn $R_R$) và quy luật tiêu tán lực hút dính theo thời gian thực dưới các kịch bản mưa có cường độ và thời gian biến thiên.

Mục tiêu tổng quát của luận án là xác lập cơ sở khoa học và định lượng hóa quy luật suy giảm hệ số ổn định ($F_s$) của mái dốc đất không bão hòa dưới tác động của các hình thái mưa khác nhau, từ đó ứng dụng vào công tác thiết kế, đánh giá an toàn công trình thủy lợi thực tế. Để hiện thực hóa mục tiêu này, hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được thiết lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Mối quan hệ tương quan giữa lực hút dính, độ ẩm thể tích và cường độ kháng cắt của các loại đất đắp điển hình miền Bắc Việt Nam biểu hiện như thế nào trong không gian ứng suất đa chiều?
  • Câu hỏi nghiên cứu 2 (RQ2): Quá trình xâm nhập của nước mưa vào mái dốc đất đắp diễn ra theo cơ chế thủy lực nào dưới sự chi phối của độ dốc mái ($m$), hệ số đầm chặt ($K$) và các dạng hình thái mưa (mưa ngắn hạn cường độ lớn - HI và mưa kéo dài cường độ nhỏ - LD)?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ khác biệt định lượng về hệ số an toàn ổn định mái dốc giữa phương pháp cơ học đất không bão hòa so với phương pháp đất bão hòa truyền thống trong điều kiện mưa thực tế là bao nhiêu?
  • Giả thuyết khoa học 1 (H1): Góc ma sát biểu kiến $\phi^b$ suy giảm phi tuyến rõ rệt khi lực hút dính $(u_a - u_w)$ vượt qua giá trị áp lực khí vào (Air Entry Value - AEV), khiến sức kháng cắt tăng chậm lại ở dải lực hút dính cao.
  • Giả thuyết khoa học 2 (H2): Hệ số đầm chặt $K$ và độ dốc mái $m$ đóng vai trò kiểm soát trực tiếp cường độ tràn ổn định ($R_R$), trong đó việc tăng độ đầm chặt làm giảm đột ngột lưu lượng thấm xâm nhập vào thân mái dốc.
  • Giả thuyết khoa học 3 (H3): Hình thái mưa kéo dài cường độ nhỏ (LD) gây suy giảm hệ số ổn định tổng thể sâu và nguy hiểm hơn so với hình thái mưa ngắn hạn cường độ lớn (HI) do khả năng làm bão hòa hoàn toàn đới không bão hòa nông.

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 nhiều pha của Fredlund & Morgenstern (1977), phương trình cường độ kháng cắt mở rộng Mohr-Coulomb (Fredlund et al., 1978), mô hình đường cong đặc trưng đất - nước của Fredlund & Xing (1994), mô hình ước lượng Modified Kovacs (MK) và phương pháp cân bằng giới hạn tổng quát (General Limit Equilibrium - GLE). Phạm vi nghiên cứu bao quát các loại đất đắp công trình phân bố tại đê 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à vật liệu mỏ đất Đại Phong. Kết quả nghiên cứu đóng góp một khối dữ liệu thực nghiệm nguyên bản và phương pháp tính toán chuẩn xác, giúp tối ưu hóa chi phí gia cố công trình từ 15% đến 30% đồng thời đảm bảo an toàn tuyệt đối cho hệ thống đê đập trước các kịch bản thời tiết cực đoan.

Literature Review và Positioning

Nghiên cứu về ổn định mái dốc đã trải qua sự tiến hóa sâu sắc từ các trường phái cơ học đất cổ điển sang lý thuyết môi trường đất không bão hòa nhiều pha. Nền tảng phân tích ổn định trượt truyền thống được định hình bởi các phương pháp phân mảnh cân bằng giới hạn (Limit Equilibrium Methods) như phương pháp mặt trượt tròn đơn giản của Fellenius (1936), Bishop đơn giản (1955), phương pháp cân bằng lực Janbu (1954), và các phương pháp thỏa mãn đồng thời cân bằng lực lẫn mô-men của Spencer (1967) hay Morgenstern-Price (1965). Mặc dù các phương pháp này đã được chuẩn hóa trong các tiêu chuẩn thiết kế quốc gia, điểm nghẽn lớn nhất của chúng nằm ở giả định đồng nhất về ứng suất hiệu quả của Terzaghi (1936): $\sigma' = \sigma - u_w$, vốn chỉ áp dụng chính xác cho môi trường đất bão hòa hai pha (hạt đất và nước).

Bước ngoặt lý thuyết diễn ra khi Bishop & Blight (1963) cùng Matyas & Radhakrishna (1968) chỉ ra rằng một biến trạng thái ứng suất đơn lẻ không thể mô tả đầy đủ ứng xử phức tạp của đất không bão hòa. Fredlund & Morgenstern (1977) đã đặt nền móng vững chắc bằng cách chứng minh đất không bão hòa là một hệ nhiệt động học bốn pha gồm: hạt đất, nước lỗ rỗng, khí lỗ rỗng và màng căng mặt ngoài khí - nước (contractile skin). Trên cơ sở đó, hai biến trạng thái ứng suất độc lập được xác lập là ứng suất pháp thực $(\sigma - u_a)$ và lực hút dính ma-trích $(u_a - u_w)$.

┌────────────────────────────────────────────────────────────────────────────┐
│                    TIẾN HÓA LÝ THUYẾT ỔN ĐỊNH MÁI DỐC                      │
├─────────────────────────────────────┬──────────────────────────────────────┤
│ TRƯỜNG PHÁI BÃO HÒA CỔ ĐIỂN         │ TRƯỜNG PHÁI KHÔNG BÃO HÒA HIỆN ĐẠI   │
│ - Terzaghi (1936): σ' = σ - uw      │ - Fredlund & Morgenstern (1977):     │
│ - Bishop (1955), Spencer (1967)     │   Hệ 4 pha, 2 biến độc lập           │
│ - Bỏ qua lực hút dính ma-trích      │   (σ - ua) và (ua - uw)              │
│ - Đánh giá quá an toàn / sai lệch   │ - Fredlund et al. (1978):            │
│   khi mực nước ngầm hạ thấp         │   τf = c' + (σ - ua)tanϕ' +          │
│                                     │        (ua - uw)tanϕb                │
└─────────────────────────────────────┴──────────────────────────────────────┘

Trong y văn thế giới, tồn tại cuộc tranh luận học thuật sâu sắc về tính chất của góc ma sát biểu kiến $\phi^b$ (đại lượng biểu thị tốc độ gia tăng sức kháng cắt theo lực hút dính). Ban đầu, Fredlund et al. (1978) đề xuất mô hình tuyến tính với $\phi^b$ là hằng số. Tuy nhiên, các thực nghiệm chuyên sâu sau đó của Escario & Saez (1986) trên đất sét Madrid, Satija (1978) trên đất sét Dhanauri, và Thu et al. (2006) trên đất bụi và sét kaolin đầm chặt đã phản bác quan điểm này. Các tác giả quốc tế chứng minh rằng $\phi^b$ mang bản chất phi tuyến cao: $\phi^b \approx \phi'$ ở dải lực hút dính thấp dưới giá trị khí vào AEV, nhưng suy giảm nhanh chóng xuống chỉ còn khoảng $\frac{1}{2}\phi' \div \frac{1}{3}\phi'$ khi $(u_a - u_w)$ vượt ngưỡng 150 - 200 kPa.

Về ảnh hưởng của mưa thấm, nghiên cứu kinh điển của Gasmo et al. (2000) và Rahardjo et al. (2002) trên các mái dốc đất tàn tích tại Singapore đã chứng minh rằng một trận mưa 80 mm/ngày làm giảm hệ số an toàn $F_s$ khoảng 25% chỉ sau 12 giờ, và cần tới 2 tuần sau khi dứt mưa để hệ số $F_s$ hồi phục lại 30%. Tương tự, các vụ sạt lở quy mô lớn tại Calabria (Ý) được phân tích bởi Gariano et al. (2015) hay thảm họa trượt lở Yabakei (tỉnh Oita, Nhật Bản, 2018) đều minh chứng cơ chế phá hoại tức thời do sự tiêu tán lực hút dính dưới tác động của mưa bão cực đoan.

Tại Việt Nam, các nghiên cứu của Hương (2013), Toan (2016), Việt et al. (2015, 2016) và Thắng (2017) đã bước đầu tiếp cận lý thuyết đất không bão hòa nhưng chủ yếu dừng lại ở các bài toán mô phỏng số hoặc thí nghiệm đơn lẻ ở quy mô mẫu nhỏ. Luận án của Phạm Huy Dũng đã định vị xuất sắc vị thế học thuật bằng cách lấp đầy khoảng trống giữa lý thuyết mô phỏng và thực tiễn hiện trường thông qua việc thiết lập một chuỗi nghiên cứu toàn diện: từ đo đạc hiện trường bằng căng kế chuyên dụng, thí nghiệm phòng xác định đường cong SWCC - hàm thấm - sức kháng cắt mở rộng, đến chế tạo hệ thống mô hình vật lý máng nghiêng tạo mưa tỷ lệ lớn và ứng dụng kiểm chứng cho 3 cụm công trình thủy lợi trọng điểm.

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

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

Luận án đã đóng góp những bước tiến quan trọng trong việc làm sáng tỏ cơ chế cơ học - thủy lực ghép nối của môi trường đất đắp không bão hòa:

  1. Mở rộng lý thuyết cường độ kháng cắt mở rộng cho đất đắp nhiệt đới: Luận án đã kiểm chứng và xác lập quy luật biến thiên phi tuyến của góc ma sát biểu kiến $\phi^b$ đối với các loại đất dính đầm chặt tại miền Bắc Việt Nam. Kết quả thực nghiệm khẳng định ở dải lực hút dính nhỏ hơn giá trị khí vào ($u_a - u_w < \text{AEV}$), góc $\phi^b$ xấp xỉ bằng góc ma sát trong hiệu quả $\phi'$ ($\phi^b \approx \phi'$). Khi lực hút dính tiếp tục gia tăng lên dải 200 - 300 kPa, $\phi^b$ suy giảm phi tuyến và duy trì tiệm cận quanh mức $12^\circ - 18^\circ$ (tương đương $0,4\phi' - 0,6\phi'$).
  2. Xác lập mô hình tương quan biến trạng thái: Chứng minh một cách thuyết phục rằng độ bền kháng cắt của đất không bão hòa không chỉ phụ thuộc độc lập vào ứng suất pháp thực $(\sigma - u_a)$ mà còn tương tác mật thiết với độ bão hòa ($S_r$) và cấu trúc lỗ rỗng hình thành qua quá trình đầm nén.
  3. Cơ chế dịch chuyển mặt trượt phá hoại do mưa: Luận án chứng minh rằng nước mưa xâm nhập làm dịch chuyển mặt trượt nguy hiểm từ vị trí cung trượt sâu truyền thống lên sát bề mặt mái dốc (cung trượt nông), giải thích bản chất vật lý của hàng loạt sự cố trượt lở mái đê, đập thực tế khi mực nước hồ/sông chưa đạt đỉnh lũ.
┌────────────────────────────────────────────────────────────────────────────┐
│                    KHUNG PHÂN TÍCH TOÀN DIỆN CỦA LUẬN ÁN                   │
│                                                                            │
│   ┌────────────────────────────────────────────────────────────────────┐   │
│   │                      LÝ THUYẾT 4 PHA ĐỘC LẬP                       │   │
│   │               (Fredlund & Morgenstern, 1977)                       │   │
│   │           Biến trạng thái: (σ - ua) và (ua - uw)                   │   │
│   └─────────────────┬────────────────────────────────┬─────────────────┘   │
│                     │                                │                     │
│                     ▼                                ▼                     │
│   ┌──────────────────────────────────┐ ┌───────────────────────────────┐   │
│   │     ỨNG XỬ THỦY LỰC KHÔNG BÃO HÒA│ │  SỨC KHÁNG CẮT MỞ RỘNG        │   │
│   │  - Mô hình SWCC Fredlund & Xing  │ │  τf = c' + (σ - ua)tanϕ' +    │   │
│   │  - Mô hình ước lượng Mod. Kovacs │ │       (ua - uw)tanϕb          │   │
│   │  - Hàm thấm kw(ψ) tích phân      │ │  - Góc ϕb phi tuyến           │   │
│   └─────────────────┬────────────────┘ └─────────────┬─────────────────┘   │
│                     │                                │                     │
│                     └────────────────┬───────────────┘                     │
│                                      ▼                                     │
│   ┌────────────────────────────────────────────────────────────────────┐   │
│   │                  MÔ HÌNH GHÉP NỐI THẤM - ỔN ĐỊNH                   │   │
│   │   - Dòng thấm không dừng (Transient Seepage: SEEP/W)               │   │
│   │   - Cân bằng giới hạn tổng quát (GLE / Morgenstern-Price: SLOPE/W) │   │
│   │   - Đánh giá Fs(t) theo kịch bản mưa HI (ngắn) và LD (dài)         │   │
│   └────────────────────────────────────────────────────────────────────┘   │
└────────────────────────────────────────────────────────────────────────────┘

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 mô hình lý thuyết cốt lõi nhằm giải quyết bài toán biến dạng và ổn định mái dốc:

  • Mô hình đường cong đặc trưng đất - nước (SWCC) của Fredlund & Xing (1994): $$\theta(\psi) = C(\psi) \frac{\theta_s}{\left[\ln\left(e + \left(\frac{\psi}{a}\right)^n\right)\right]^m}$$ Trong đó $C(\psi) = 1 - \frac{\ln(1 + \psi/\psi_r)}{\ln(1 + 10^6/\psi_r)}$ là hàm hiệu chỉnh; $a, n, m$ là các tham số hình dạng đường cong; $\psi = (u_a - u_w)$ là lực hút dính ma-trích.
  • Mô hình ước lượng Modified Kovacs (MK): Cho phép xác lập SWCC dựa trên các chỉ tiêu vật lý cơ bản (thành phần hạt $D_{10}, C_u$, giới hạn dẻo $W_p$, giới hạn chảy $W_L$, độ rỗng $e$). Luận án đã tinh chỉnh tham số để mô hình MK phản ánh chính xác đất đắp miền Bắc Việt Nam khi thiếu điều kiện thí nghiệm chuyên sâu.
  • Phương trình thấm không dừng Darcy mở rộng: $$\frac{\partial}{\partial x}\left(k_x \frac{\partial h}{\partial x}\right) + \frac{\partial}{\partial y}\left(k_y \frac{\partial h}{\partial y}\right) + Q = m_w \gamma_w \frac{\partial h}{\partial t}$$ Với $k_x, k_y$ là hệ số thấm không bão hòa thay đổi theo lực hút dính; $m_w = \frac{\partial \theta}{\partial \psi}$ là độ dốc của SWCC.

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan khoa học thực chứng (Positivism) kết hợp với chủ nghĩa hiện thực phê phán (Critical Realism). Thiết kế nghiên cứu được thực hiện theo cấu trúc đa tầng, kết hợp đa phương pháp (triangulation):

  • Cấp độ vi mô & trung mô (Micro/Meso-level): Thí nghiệm trong phòng trên các mẫu đất nguyên dạng và đầm chặt nhằm xác định bản chất tương tác hạt - nước - khí.
  • Cấp độ mô hình vật lý tỷ lệ lớn (Physical Modeling): Chế tạo hệ thống máng nghiêng trọng lực có kích thước lớn cùng dàn mưa nhân tạo để mô phỏng chính xác ứng xử thủy lực dòng tràn - dòng ngấm.
  • Cấp độ vĩ mô (Macro-level / Field-scale): Quan trắc thực địa bằng hệ thống căng kế tại các công trình đê đập và mô phỏng số kiểm chứng trên phần mềm chuyên dụng GeoStudio.
┌────────────────────────────────────────────────────────────────────────────┐
│                    THIẾT KẾ NGHIÊN CỨU ĐA TẦNG (TRIANGULATION)             │
├────────────────────────────────────────────────────────────────────────────┤
│ 1. THỰC NGHIỆM HIỆN TRƯỜNG                                                 │
│    - Quan trắc lực hút dính (ua - uw) bằng Tensometer 2725ARL-Jetfill      │
│    - Khảo sát đập Khau Piều, đập Chúc Bài Sơn, đập Bầu Lầy, đê hữu Cầu    │
│    - Hiệu chỉnh cột nước đo: Δp = 15,2 kPa ứng với ống dài 1,5 m           │
├────────────────────────────────────────────────────────────────────────────┤
│ 2. THÍ NGHIỆM PHÒNG CHUYÊN SÂU                                             │
│    - Bình chiết áp lực cao (Pressure Plate Extractor) + Đĩa gốm AEV cao    │
│    - Máy nén 3 trục tự động TRIAX50 (sơ đồ CD, CU, CW)                     │
│    - Máy cắt trực tiếp không bão hòa EDJ-2 cải tiến                        │
├────────────────────────────────────────────────────────────────────────────┤
│ 3. MÔ HÌNH VẬT LÝ MÁNG NGHIÊNG TỶ LỆ LỚN (FLUME TEST)                      │
│    - Máng trượt biến đổi độ dốc (m = 1:1, 1:2, 1:4)                        │
│    - Dàn phun mưa nhân tạo kiểm soát cường độ mưa chính xác                │
│    - Căng kế điện tử 2100F-Remote Tensometer + ICT Data Logger             │
├────────────────────────────────────────────────────────────────────────────┤
│ 4. MÔ PHỎNG SỐ VÀ BACK-ANALYSIS                                            │
│    - SEEP/W (Phân tích trường áp lực lỗ rỗng không dừng)                   │
│    - SLOPE/W (Phân tích ổn định cân bằng giới hạn tổng quát GLE)          │
└────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được chuẩn hóa với độ tin cậy và giá trị khoa học cao:

  • Đo đạc hiện trường: Sử dụng căng kế cơ học kiểu 2725ARL-Jetfill (Soilmoisture Equipment Corp, Mỹ). Các điểm quan trắc được bố trí tại đỉnh, mái và chân dốc ở độ sâu $0,3 \div 1,5\text{ m}$. Luận án áp dụng công thức hiệu chỉnh áp lực thủy tĩnh trong ống dẫn: $$u_w = u_{w,\text{đo}} + \gamma_w \cdot z_{\text{cột nước}}$$ (Ví dụ: chiều cao cột nước trong căng kế là $1,5\text{ m}$ thì giá trị hiệu chỉnh áp lực tương ứng là $15,2\text{ kPa}$).
  • Thí nghiệm xác định SWCC trong phòng: Sử dụng bình chiết áp lực khí cao (Pressure Plate Extractor) với đĩa gốm tiếp nhận khí cao (High Air Entry Ceramic Disk). Mẫu đất được bão hòa hoàn toàn trong chân không, sau đó áp các cấp áp lực khí nén tăng dần ($0, 10, 30, 100, 300, 500, 1000, 1500\text{ kPa}$) theo kỹ thuật dịch chuyển trục (Axis-translation technique) để ép nước thoát ra cho đến khi đạt trạng thái cân bằng thủy lực tại mỗi cấp áp lực.
  • Thí nghiệm sức kháng cắt không bão hòa: Thực hiện đồng thời trên máy nén 3 trục TRIAX50 và máy cắt trực tiếp EDJ-2. Các sơ đồ thí nghiệm cố kết thoát nước (CD) và độ ẩm không đổi (CW) được tiến hành với các cấp áp lực buồng khí $\sigma_3$, áp lực khí lỗ rỗng $u_a$ và áp lực nước lỗ rỗng $u_w$ được kiểm soát độc lập.
  • Mô hình máng thí nghiệm mưa nhân tạo: Chế tạo máng thép chuyên dụng kích thước lớn có trục quay tự do cho phép thay đổi hệ số mái dốc ($m = 1; 2; 4$). Hệ thống dàn mưa được thiết kế với đầu phun áp lực đồng đều. Bên trong khối đất đắp (mỏ Đại Phong ở các độ đầm chặt $K = 0,70; 0,85; 0,95; 0,97$), gắn hệ thống căng kế điện tử 2100F-Remote Tensometer kết nối trực tiếp với bộ xử lý trung tâm và phần mềm ICT Data Logger để thu thập số liệu tự động với tần suất 1 phút/lần.

Data và phân tích

Dữ liệu thực nghiệm của luận án phản ánh đầy đủ các đặc tính cơ lý - thủy lực của các đối tượng nghiên cứu:

Đối tượng / Vị trí Loại đất Khối lượng thể tích khô $\rho_d$ ($\text{g/cm}^3$) Độ ẩm tối ưu $W_{\text{opt}}$ (%) Giới hạn dẻo $W_p$ / Chảy $W_L$ (%) Lực dính $c'$ (kPa) Góc ma sát $\phi'$ (độ) Góc ma sát $\phi^b$ (độ) Áp lực khí vào AEV (kPa)
Đê hữu Cầu (Bắc Ninh) Sét pha nhẹ $1,62 \div 1,68$ $16,8$ $18,2 / 32,5$ $18,5$ $21^\circ 30'$ $13^\circ 15'$ $18,0$
Đập Khau Piều (Lạng Sơn) Sét pha nặng $1,55 \div 1,60$ $19,5$ $21,4 / 38,6$ $24,2$ $18^\circ 45'$ $11^\circ 30'$ $28,5$
Đập Chúc Bài Sơn (Quảng Ninh) Sét nhiều dăm sạn $1,70 \div 1,75$ $14,2$ $16,5 / 29,8$ $21,0$ $24^\circ 10'$ $15^\circ 40'$ $14,5$
Mỏ Đại Phong (thí nghiệm máng) Á sét $1,58 \div 1,72$ $17,5$ $19,0 / 34,2$ $22,0$ $22^\circ 00'$ $14^\circ 00'$ $22,0$

Các kiểm định độ tin cậy được thực hiện thông qua việc đối chiếu giữa số liệu đo đạc thực tế của căng kế ngoài hiện trường với kết quả dự báo từ đường cong SWCC ước lượng theo mô hình MK, sai số quân phương (RMSE) luôn nằm trong giới hạn kiểm soát chặt chẽ ($< 5%$).

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

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

  1. Quy luật phân chia dòng chảy tràn và dòng thấm mặt: Kết quả thí nghiệm máng khẳng định độ đầm chặt ($K$) và hệ số mái ($m$) quyết định hoàn toàn tỷ lệ nước mưa xâm nhập ($R_I$). Khi hệ số đầm chặt tăng từ $K = 0,70$ lên $K = 0,95$, thời gian bắt đầu xuất hiện dòng chảy tràn rút ngắn từ 45 phút xuống chỉ còn 8 phút; cường độ tràn ổn định ($R_R$) tăng từ 35% lên đến 82% tổng lượng mưa. Đất đầm chặt cao có cấu trúc lỗ rỗng mao dẫn nhỏ làm giảm vận tốc thấm ban đầu, khiến phần lớn nước mưa chuyển thành dòng chảy mặt gây xói lở bề mặt thay vì ngấm sâu tức thì.
  2. Động thái trễ và sự suy giảm phân đới của lực hút dính: Lực hút dính $(u_a - u_w)$ tại tầng nông sát mặt dốc ($0 \div 30\text{ cm}$) suy giảm gần như triệt để (từ $60 - 80\text{ kPa}$ về xấp xỉ $0\text{ kPa}$) chỉ sau $12 \div 24\text{ giờ}$ mưa liên tục. Tuy nhiên, ở độ sâu $> 1,0\text{ m}$, tốc độ suy giảm lực hút dính diễn ra có độ trễ thủy lực từ $24 \div 48\text{ giờ}$. Sau khi dứt mưa, quá trình bay hơi và tái phân bố ẩm diễn ra rất chậm; để phục hồi 50% giá trị lực hút dính ban đầu cần thời gian không mưa kéo dài từ $7 \div 15\text{ ngày}$.
  3. Sự tương phản phá hoại giữa hai hình thái mưa (HI vs LD):
    • Kiểu mưa HI (Heavy Intensity - Cường độ lớn, thời gian ngắn): Do cường độ mưa vượt xa hệ số thấm bão hòa ($I \gg k_s$), phần lớn nước tạo thành dòng chảy tràn gây xói mòn rãnh và phá hoại nông cục bộ. Hệ số ổn định $F_s$ suy giảm nhanh trong thời gian ngắn nhưng mức độ suy giảm tổng thể chỉ dao động từ $12% \div 18%$.
    • Kiểu mưa LD (Long Duration - Kéo dài, cường độ nhỏ hoặc trung bình): Do $I \approx k_s$, toàn bộ nước mưa thấm trọn vẹn vào thân đất ($R_I \to 100%$), triệt tiêu vùng áp lực nước lỗ rỗng âm, dâng cao mực bão hòa tạm thời (perched water table). Kết quả phân tích tại đập Khau Piều và đê hữu Cầu cho thấy kiểu mưa LD làm hệ số an toàn $F_s$ suy giảm mạnh nhất lên tới $28% \div 35%$, đẩy mái dốc vào trạng thái mất ổn định tiệm cận phá hoại ($F_s < 1,05$).
  4. Hiệu ứng gia cường của sức kháng cắt không bão hòa: Trong điều kiện vận hành bình thường (mực nước dâng bình thường, không mưa), việc tính toán xét đến lực hút dính giúp hệ số an toàn $F_s$ của mái dốc tăng từ $18% \div 32%$ so với phương pháp bão hòa cổ điển. Điều này giải thích tại sao nhiều mái dốc tự nhiên và công trình thực tế có góc nghiêng lớn hơn góc ma sát trong $\phi'$ nhưng vẫn đứng vững ổn định qua nhiều năm khi chưa gặp đợt mưa kéo dài cực đoan.
┌────────────────────────────────────────────────────────────────────────────┐
│                  SO SÁNH TÁC ĐỘNG CỦA HAI HÌNH THÁI MƯA                    │
├─────────────────────────────────────┬──────────────────────────────────────┤
│ KIỂU MƯA HI (Cường độ lớn, ngắn hạn)│ KIỂU MƯA LD (Kéo dài, cường độ nhỏ)  │
├─────────────────────────────────────┼──────────────────────────────────────┤
│ - Cường độ mưa: I >> ks             │ - Cường độ mưa: I ≈ ks               │
│ - Dòng chảy tràn: Lớn (RR > 80%)    │ - Dòng thấm thâm nhập: Lớn (RI ≈ 100%)│
│ - Độ sâu bão hòa: Nông (< 0,5 m)    │ - Độ sâu bão hòa: Sâu (> 1,5 - 3,0 m) │
│ - Cơ chế phá hoại: Xói mòn bề mặt,  │ - Cơ chế phá hoại: Bão hòa toàn phần,│
│   trượt nông cục bộ                 │   triệt tiêu lực hút dính, trượt sâu │
│ - Mức giảm hệ số an toàn: 12% - 18% │ - Mức giảm hệ số an toàn: 28% - 35%  │
│ - Nguy cơ tiềm ẩn: Trung bình       │ - Nguy cơ tiềm ẩn: Cực kỳ nguy hiểm  │
└─────────────────────────────────────┴──────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết & học thuật: Đặt nền móng dữ liệu thực nghiệm chuẩn tắc cho phân ngành cơ học đất không bão hòa tại Việt Nam; cung cấp phương trình hiệu chỉnh SWCC và thông số $\phi^b$ tin cậy cho các nghiên cứu tiếp theo.
  • Về quy chuẩn và phương pháp thiết kế: Đề xuất sửa đổi, bổ sung các tiêu chuẩn thiết kế đê, đập hiện hành (TCVN). Khuyến nghị bắt buộc đưa bài toán thấm không dừng và kịch bản mưa kéo dài (LD) vào quy trình thẩm tra an toàn ổn định công trình thủy lợi cấp III trở lên.
  • Về giải pháp công trình thực tiễn:
    • Thiết kế hệ thống tiêu thoát nước mặt kết hợp rãnh thu nước ngắt dòng trên mái dốc để hạn chế tối đa thời gian lưu nước của mưa kiểu LD.
    • Áp dụng các giải pháp bảo vệ bề mặt (trồng cỏ vetiver, phủ bạt địa kỹ thuật kết hợp khung bê tông lưới neo) nhằm bảo tồn lực hút dính ma-trích trong thân đắp và triệt tiêu nguy cơ xói ngầ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 một số giới hạn nghiên cứu (boundary conditions):

  1. Giới hạn đo đạc lực hút dính: Thiết bị căng kế cơ học và điện tử sử dụng trong nghiên cứu bị giới hạn đo tối đa ở mức $-90\text{ kPa}$ do hiện tượng sinh bọt khí (cavitation) của nước trong ống dẫn. Đối với các loại đất sét có độ chặt cao ở trạng thái rất khô, lực hút dính có thể lên tới hàng nghìn kPa chưa thể đo trực tiếp liên tục tại hiện trường mà phải dựa vào đường cong ngoại suy.
  2. Chưa xét đến hiện tượng trễ chu kỳ khô - ướt (Hysteresis): Thí nghiệm xác định SWCC và mô hình máng chủ yếu khảo sát quá trình hút nước (wetting path) đơn lẻ trong một đợt mưa, chưa mô phỏng toàn diện chu kỳ lặp lại nhiều lần của quá trình khô hạn - ngập mưa (wetting - drying cycles), vốn làm thay đổi dần cấu trúc vi mô của đất.
  3. Mô hình ghép nối cơ - thủy văn một chiều: Mô hình mô phỏng GeoStudio giải quyết bài toán thấm không dừng độc lập rồi chuyển trường áp lực nước lỗ rỗng sang bài toán ổn định (uncoupled/semi-coupled analysis), chưa tính toán đồng thời biến dạng thể tích co ngót - trương nở tương tác ngược trở lại trường thấm (fully coupled THMC analysis).

Chương trình nghiên cứu phát triển 10 năm tiếp theo:

  • Phát triển hệ thống cảm biến đo lực hút dính dải rộng (High-capacity tensometers lên tới $1500\text{ kPa}$) ứng dụng công nghệ vi cơ điện tử (MEMS) phục vụ quan trắc tự động thời gian thực tại các hồ đập trọng điểm quốc gia.
  • Nghiên cứu cơ chế nứt nẻ do co ngót trong mùa khô và vai trò dẫn dòng siêu tốc của các khe nứt này khi bước vào mùa mưa lũ.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning / Deep Learning) kết hợp dữ liệu lớn viễn thám và radar thời tiết để xây dựng bản đồ cảnh báo sớm nguy cơ sạt lở thời gian thực dựa trên ngưỡng mưa tích lũy và suy giảm lực hút dính động.

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

  • Tác động học thuật: Luận án mở ra hướng nghiên cứu mũi nhọn về địa kỹ thuật môi trường và cơ học đất không bão hòa tại các cơ sở đào tạo chuyên ngành xây dựng - thủy lợi tại Việt Nam. Các kết quả công bố trên các tạp chí chuyên ngành uy tín đã tạo tiền đề cho nhiều đề tài cấp Bộ và Nhà nước (tiêu biểu như Chương trình KC).
  • Chuyển đổi thực tiễn ngành xây dựng: Cung cấp giải pháp định lượng giúp các đơn vị tư vấn thiết kế tính toán chính xác mái dốc đắp đê đập, tránh lãng phí hàng trăm tỷ đồng do việc thiết kế hệ số mái quá thoải một cách cảm tính hoặc giảm thiểu tối đa các sự cố sạt trượt trong giai đoạn thi công và vận hành.
  • Tác động chính sách và an sinh xã hội: Cung cấp cơ sở khoa học vững chắc cho Ban Chỉ đạo Quốc gia về Phòng, chống thiên tai và Cục Quản lý Đê điều & Phòng, chống Thiên tai (Bộ Nông nghiệp và PTNT) trong việc ban hành các chỉ dẫn kỹ thuật ứng phó sự cố đê điều trong mùa bão lũ.

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

┌────────────────────────────────────────────────────────────────────────────┐
│                         ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                     │
├────────────────────────────────────────────────────────────────────────────┤
│ • NGHIÊN CỨU SINH & HỌC THUẬT: Kế thừa cơ sở dữ liệu SWCC, hàm thấm        │
│   và tham số ϕb cho đất miền Bắc Việt Nam; chuẩn hóa quy trình flume test. │
│                                                                            │
│ • CÁC NHÀ KHOA HỌC ĐỊA KỸ THUẬT: Khung phân tích kết hợp mô hình MK và     │
│   Fredlund & Xing giúp mở rộng nghiên cứu sang đất tàn tích sườn dốc tự    │
│   nhiên miền Trung và Tây Nguyên.                                          │
│                                                                            │
│ • CÔNG TY TƯ VẤN & THIẾT KẾ THỦY LỢI/GIAO THÔNG: Ứng dụng phần mềm Seep/W, │
│   Slope/W tính toán chính xác Fs, tối ưu kích thước hình học thân đê đập.  │
│                                                                            │
│ • CƠ QUAN QUẢN LÝ NHÀ NƯỚC (Bộ NN&PTNT, Bộ GTVT): Bộ công cụ đánh giá      │
│   ngưỡng an toàn đê đập theo lượng mưa dự báo, nâng cao năng lực ứng phó. │
└────────────────────────────────────────────────────────────────────────────┘

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 và chuẩn hóa lý thuyết cường độ kháng cắt mở rộng Mohr-Coulomb của Fredlund et al. (1978) cho các loại đất dính đầm chặt đa nguồn gốc tại Việt Nam thông qua việc chứng minh tính chất phi tuyến của góc ma sát biểu kiến $\phi^b$. Luận án chứng minh rõ: khi lực hút dính $(u_a - u_w) \le \text{AEV}$, sức kháng cắt gia tăng tuyến tính với $\phi^b \approx \phi'$; nhưng khi vượt qua AEV, $\phi^b$ suy giảm phi tuyến mạnh và tiệm cận về giá trị ổn định $\frac{1}{2}\phi' \div \frac{1}{3}\phi'$. Điều này khắc phục triệt để sai số đánh giá quá cao sức kháng cắt nếu giả định $\phi^b$ là hằng số cố định trong dải lực hút dính lớn.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế là gì?

So với các nghiên cứu trước đây (như Gasmo et al., 2000; Rahardjo et al., 2002 vốn chủ yếu dựa trên mô hình số hoặc quan trắc hiện trường thuần túy), luận án đã thiết kế và chế tạo thành công hệ thống mô hình vật lý máng nghiêng trọng lực tỷ lệ lớn tích hợp dàn phun mưa nhân tạo chính xác và hệ thống căng kế điện tử 2100F-Remote Tensometer kết nối ICT Data Logger. Mô hình này cho phép đo đạc đồng thời, liên tục với độ chính xác cao cả 3 đại lượng thủy lực: cường độ chảy tràn ($R_R$), lượng nước thấm vào thân đất ($R_I$) và sự biến thiên tức thời của áp lực nước lỗ rỗng âm theo chiều sâu dưới các góc dốc và độ đầm chặt khác nhau.

3. Phát hiện bất ngờ và phản trực giác nhất từ số liệu thực nghiệm là gì?

Phát hiện phản trực giác nhất là: Mưa lớn cường độ cao trong thời gian ngắn (kiểu HI) ít gây nguy cơ trượt lở tổng thể thân đê/đập hơn so với mưa nhỏ kéo dài (kiểu LD). Trực giác thông thường cho rằng mưa như trút nước sẽ gây sạt lở nghiêm trọng nhất. Tuy nhiên, dữ liệu thực nghiệm chứng minh mưa kiểu HI vượt quá hệ số thấm bão hòa của đất đầm chặt ($I \gg k_s$), khiến $80% - 85%$ lượng nước nhanh chóng biến thành dòng chảy tràn bề mặt. Ngược lại, mưa kiểu LD có cường độ tương thích với hệ số thấm ($I \approx k_s$), nước mưa thấm trọn vẹn $100%$ vào thân đắp, ngấm sâu và triệt tiêu hoàn toàn lực hút dính ma-trích trong toàn bộ khối đất, làm hệ số an toàn $F_s$ sụt giảm nghiêm trọng nhất (giảm đến 35%).

4. Quy trình lặp lại thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp chi tiết toàn bộ quy trình chế bị mẫu đất, các bước vận hành thí nghiệm bình chiết áp lực khí cao, quy chuẩn cắt 3 trục không bão hòa, bản vẽ thiết kế máng thí nghiệm, thông số kỹ thuật dàn mưa và sơ đồ lắp đặt căng kế hiện trường. Mọi tham số đầu vào của đất (thành phần hạt, giới hạn Atterberg, tỷ trọng, độ đầm chặt, hệ số thấm) đều được lập bảng chi tiết, cho phép các phòng thí nghiệm địa kỹ thuật độc lập tái lập hoàn toàn các kết quả đo đạc.

5. Khung nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chiến lược nghiên cứu dài hạn tập trung vào 3 trọng tâm: (1) Nghiên cứu bài toán tương tác biến dạng nứt nẻ co ngót - thấm - trượt lở của đê đập dưới tác động của biến đổi khí hậu cực đoan; (2) Chế tạo cảm biến đo lực hút dính dải đo siêu cao ($> 1500\text{ kPa}$) không bị ảnh hưởng bởi bọt khí; (3) Tích hợp mô hình số địa kỹ thuật với mạng lưới IoT quan trắc tự động và hệ thống trí tuệ nhân tạo để phát triển nền tảng dự báo trượt lở thông minh thời gian thực cấp quốc gia.

Kết luận

Luận án tiến sĩ kỹ thuật của NCS. Phạm Huy Dũng là công trình khoa học tiên phong, chuẩn mực và toàn diện trong lĩnh vực cơ học đất không bão hòa tại Việt Nam. Những giá trị cốt lõi và di sản học thuật mà công trình xác lập bao gồm:

  1. Xây dựng bộ cơ sở dữ liệu thực nghiệm hoàn chỉnh về đường cong đặc trưng đất - nước (SWCC), hàm thấm và các tham số cường độ kháng cắt mở rộng ($\phi^b, c'$) cho các loại đất đắp công trình thủy lợi tiêu biểu miền Bắc Việt Nam.
  2. Định lượng hóa thành công quy luật phân lưu thủy lực giữa dòng chảy tràn bề mặt và dòng thấm vào thân đất dưới ảnh hưởng đồng thời của hệ số đầm chặt ($K$) và hệ số mái ($m$).
  3. Làm sáng tỏ cơ chế tiêu tán lực hút dính ma-trích và xác định hình thái mưa kéo dài cường độ nhỏ (LD) là nhân tố nguy hiểm nhất làm suy giảm hệ số ổn định tổng thể của mái dốc đất đắp.
  4. Chuẩn hóa công cụ ước lượng SWCC thông qua mô hình Modified Kovacs (MK) hiệu chỉnh, mở ra giải pháp kỹ thuật tin cậy và tiết kiệm chi phí cho các dự án thiếu thiết bị thí nghiệm chuyên sâu.
  5. Khẳng định tính ưu việt của phương pháp tính toán ổn định trên cơ sở khoa học đất không bão hòa, góp phần nâng cao độ chính xác, an toàn và tối ưu hóa kinh tế trong thiết kế, xây dựng và quản lý vận hành hệ thống đê điều, hồ đập quốc gia.