Tổng quan về luận án

Mất ổn định mái dốc và trượt lở đất dưới tác động của mưa là hiểm họa địa kỹ thuật nghiêm trọng hàng đầu tại Việt Nam cũng như trên toàn cầu. Với 74,3% diện tích tự nhiên là đồi núi và trung du, cùng mạng lưới hạ tầng thủy lợi đồ sộ gồm hơn 6.000 hồ chứa đập đất và hàng chục nghìn kilômét đê điều, sự an toàn của các công trình đất đắp chịu áp lực thủy triều, lũ bão và mưa cực đoan là bài toán sống còn. 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" (Chuyên ngành: Địa kỹ thuật xây dựng; Mã số: 9580211) do nghiên cứu sinh Phạm Huy Dũng thực hiện 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ại Trường Đại học Thủy lợi (Hà Nội, 2020) đã tiếp cận trực diện vấn đề cốt lõi này bằng nền tảng cơ học đất không bão hòa (Unsaturated Soil Mechanics).

Trong thực hành kỹ thuật truyền thống, các phân tích ổn định mái dốc dựa trên cơ học đất bão hòa cổ điển của Karl Terzaghi thường giả thiết đất bão hòa hoàn toàn dưới mực nước ngầm và khô hoàn toàn trên mực nước ngầm, hoàn toàn bỏ qua vai trò của áp lực nước lỗ rỗng âm (lực hút dính ma trận - matric suction). Khoảng trống nghiên cứu (research gap) lớn nhất hiện nay là sự thiếu hụt nghiêm trọng cơ sở dữ liệu thực nghiệm về đường cong đặc trưng đất - nước (SWCC), hàm hệ số thấm không bão hòa ($k_w$) và các tham số kháng cắt phi tuyến của đất đắp công trình miền Bắc Việt Nam. Đồng thời, cơ chế thấm xâm nhập của nước mưa vào mái dốc có độ dốc ($m$) và độ đầm chặt ($K$) khác nhau, cùng quy luật suy giảm tức thời và biến thiên trễ của lực hút dính trong và sau khi dứt mưa chưa được thực nghiệm hóa một cách định lượng bằng mô hình vật lý quy mô lớn.

Luận án giải quyết hai câu hỏi nghiên cứu trung tâm:

  1. Quy luật biến đổi của lực hút dính ma trận $(u_a - u_w)$ và cường độ kháng cắt $\tau_f$ của các loại đất đắp đê, đập miền Bắc Việt Nam diễn ra như thế nào trong điều kiện độ ẩm biến thiên?
  2. Cơ chế phân bổ giữa lượng nước mưa chảy tràn bề mặt ($RR$) và thấm xâm nhập ($RI$) tương tác ra sao với độ dốc mái, hệ số đầm chặt, từ đó kiểm soát sự suy giảm hệ số ổn định an toàn ($FS$) dưới các kịch bản mưa cường độ cao thời gian ngắn (HI) và mưa kéo dài cường độ nhỏ (LD)?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • $H_1$: Áp lực nước lỗ rỗng âm đóng góp trực tiếp vào thành phần lực dính biểu kiến của đất không bão hòa thông qua góc ma sát biểu kiến $\phi^b$, làm gia tăng đáng kể hệ số an toàn $FS$ so với trạng thái bão hòa quy ước.
  • $H_2$: Dưới tác động của mưa, sự triệt tiêu lực hút dính ma trận diễn ra bất đối xứng theo chiều sâu và hình thái trận mưa, trong đó các trận mưa kéo dài (LD) gây ra nguy cơ mất ổn định nghiêm trọng hơn các trận mưa rào cục bộ (HI) do làm ẩm sâu tầng trượt tiềm năng.

Nghiên cứu tích hợp khung lý thuyết cơ học liên tục nhiều pha của Fredlund & Morgenstern (1977), mô hình đường cong đặc trưng đất - nước Fredlund & Xing (1994), mô hình Modified Kovacs (MK), phương trình cường độ kháng cắt mở rộng Fredlund et al. (1978) và phương pháp cân bằng giới hạn tổng quát (General Limit Equilibrium - GLE). Đóng góp mang tính đột phá của công trình là thiết lập hệ thống máng mô hình trọng lực biến thiên góc nghiêng tích hợp dàn mưa nhân tạo và mạng lưới cảm biến đo áp lực nước lỗ rỗng từ xa (Remote Tensometer 2100F), từ đó định lượng hóa tỉ lệ nước mưa xâm nhập ($RI$) đạt từ 14% đến 34% và chứng minh quy luật trễ thủy lực sau mưa. 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 và mỏ đất Đại Phong (Quảng Ninh).

Literature Review và Positioning

Lịch sử phát triển của các phương pháp phân tích ổn định mái dốc bắt đầu từ các kỹ thuật phân mảnh cân bằng giới hạn tĩnh học 2D như Fellenius (1936), Bishop (1955), Janbu (1954), Spencer (1967), Morgenstern-Price (1965) và Sarma (1973). Mặc dù các phương pháp phần tử hữu hạn (FEM) và phân tích giới hạn dẻo (Limit Analysis) phát triển mạnh mẽ (Griffiths & Lane, 1999), phương pháp cân bằng giới hạn tổng quát (GLE) vẫn giữ vị trí chủ đạo trong tiêu chuẩn thiết kế quốc tế nhờ khả năng thỏa mãn đồng thời cân bằng lực và cân bằng mô-men. Tuy nhiên, rào cản kỹ thuật của các phương pháp này trong thế kỷ XX là việc chỉ sử dụng cặp tham số cường độ kháng cắt hiệu quả $(c', \phi')$ của đất bão hòa.

Trường phái cơ học đất không bão hòa hiện đại hình thành bước ngoặt khi Fredlund & Rahardjo (1993) chuẩn hóa việc 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 ma trận $(u_a - u_w)$. Tranh luận học thuật quốc tế tập trung vào hai quan điểm đối nghịch về xử lý lực hút dính trong tính toán ổn định:

  1. Quan điểm "Lực dính toàn phần" (Total Cohesion Approach) của Ching et al. (1984), xem lực hút dính làm tăng trực tiếp lực dính của đất theo biểu thức $c = c' + (u_a - u_w)\tan\phi^b$.
  2. Quan điểm "Cường độ kháng cắt mở rộng" (Extended Shear Strength Theory) của Fredlund et al. (1978), phân tách không gian ứng suất thành mặt phá hoại 3D phi tuyến.

Trên bình diện quốc tế, các công trình kinh điển đã chứng minh sự suy giảm ổn định mái dốc do mưa. Gasmo et al. (2000) phân tích mái dốc đất tàn tích tại Singapore và chỉ ra rằng trận mưa cường độ 80 mm/ngày làm suy giảm 25% hệ số an toàn $FS$ chỉ sau 12 giờ, và cần tới 2 tuần sau khi tạnh mưa để $FS$ hồi phục 30%. Rahardjo et al. (2002) nghiên cứu mái dốc tàn tích granite tại Singapore chịu mưa 80 mm/h trong 4 giờ, khẳng định chiều sâu mặt trượt nguy hiểm bị kéo cạn lên gần mặt dốc do sự triệt tiêu lực hút dính tầng nông. Tại châu Âu, Gariano et al. (2015) ghi nhận 1.466 vụ trượt lở đất tại Calabria (Ý) từ năm 1921 đến 2010, nhấn mạnh tần suất sạt lở tăng đột biến do biến đổi khí hậu làm gia tăng các trận mưa tích lũy. Tại Nhật Bản, thảm họa sạt lở Yabakei (tỉnh Oita, 2018) với thể tích khối trượt 220 m $\times$ 160 m $\times$ 35 m là minh chứng điển hình cho sự phá hủy liên kết mao dẫn của đất không bão hòa sau các đợt mưa kéo dài.

Tại Việt Nam, các nghiên cứu của Hoàng & Khang (2011) tại Cốc Pài (Hà Giang), Hương (2013) tại Quảng Ninh - Yên Bái - Ninh Thuận, Việt et al. (2015, 2016) tại Yên Bái, và Thắng (2017) tại đê tả Đuống đã đặt những viên gạch đầu tiên. Tuy vậy, các tác giả chủ yếu áp dụng mô hình số ngoại suy hoặc sử dụng thiết bị thí nghiệm gửi ra nước ngoài (như Đại học Ibaraki, Nhật Bản). Luận án của Phạm Huy Dũng (2020) định vị chính xác khoảng trống này bằng cách tự chủ công nghệ thí nghiệm trong nước, kết hợp thí nghiệm ba trục TRIAX50, cắt trực tiếp EDJ-2, bình chiết áp lực cao với đĩa sứ tiếp nhận khí cao (High Air Entry Ceramic Disc), hệ thống Jetfill Tensometer 2725ARL và mô hình hóa vật lý tương tác mưa - mái dốc đất đắp.

Đó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 trực tiếp lý thuyết ứng suất hiệu quả của Karl Terzaghi (1936) sang miền không bão hòa. Terzaghi định nghĩa biến trạng thái ứng suất kiểm soát đất bão hòa:

"Ứng suất tại bất kỳ điểm nào của một mặt cắt qua khối đất có thể được tính toán từ ứng suất chính tổng $\sigma_1, \sigma_2, \sigma_3$ tác dụng tại điểm đó... Tất cả các ảnh hưởng đo được do ứng suất thay đổi như nén, xoắn và thay đổi sức kháng cắt đều là do các thay đổi về ứng suất hiệu quả $\sigma'_1, \sigma'_2, \sigma'_3$."

Luận án kế thừa và phát triển phương trình trạng thái cường độ kháng cắt mở rộng theo Fredlund et al. (1978):

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

Trong đó, đóng góp lý thuyết đặc thù là làm sáng tỏ tính phi tuyến mạnh mẽ của góc ma sát biểu kiến $\phi^b$ đối với đất đắp công trình thủy lợi miền Bắc Việt Nam:

  • Khi lực hút dính $(u_a - u_w) \le AEV$ (Air-Entry Value - Giá trị khí vào, dao động từ 11,8 kPa đến 40,0 kPa), góc ma sát biểu kiến xấp xỉ góc ma sát trong hiệu quả: $\phi^b \approx \phi'$.
  • Khi lực hút dính tăng cao trong vùng chuyển tiếp và tàng dư ($(u_a - u_w) > 50 \div 70\text{ kPa}$), $\phi^b$ suy giảm phi tuyến và tiệm cận về khoảng $\frac{1}{2}\phi' \div \frac{1}{3}\phi'$.

Sự thay đổi này được chứng minh bằng thực nghiệm trên đất sét Dhanauri (Satija, 1978: $\phi' = 29^\circ, \phi^b = 16,2^\circ$), đất sét Madrid (Escario & Saez, 1986: $\phi' = 22,5^\circ, \phi^b = 16,1^\circ$), và các mẫu đất thực tế tại đê hữu Cầu, đập Khau Piều và đập Chúc Bài Sơn.

Mô hình lý thuyết đề xuất cơ chế phá hoại trượt mái dốc đất không bão hòa gồm 3 mệnh đề (propositions):

  • Mệnh đề 1: Áp lực nước lỗ rỗng âm duy trì một màng sức căng mặt ngoài $T_s$ với bán kính cong $R_s$, tạo nên lực liên kết hạt bổ sung tỉ lệ thuận với $(u_a - u_w)\tan\phi^b$.
  • Mệnh đề 2: Khi mưa xâm nhập, mặt ẩm di chuyển xuống dưới làm bão hòa cục bộ tầng đất mặt, triệt tiêu $(u_a - u_w) \to 0$, khiến thành phần kháng cắt do lực hút dính biến mất, gây phá hủy trượt nông cục bộ trước khi xuất hiện mực nước ngầm dâng cao.
  • Mệnh đề 3: Tốc độ suy giảm hệ số ổn định $FS$ bị chi phối bởi độ dốc đường cong đặc trưng đất - nước ($m_w$) và đạo hàm của hàm thấm $k_w(\psi)$.
   [Trạng thái ban đầu: Đất KBH]
   Lực hút dính (ua - uw) cao ---> Sức căng mặt ngoài Ts ---> Tăng cường độ kháng cắt tau_f ---> FS cao (an toàn)
                 |
                 | (Nước mưa xâm nhập: Lượng thấm RI = 14% - 34%)
                 v
   [Quá trình mưa: Di chuyển mặt thấm]
   Độ ẩm thể tích theta tăng ---> Triệt tiêu lực hút dính (ua - uw -> 0) ---> Giảm góc phi_b
                 |
                 v
   [Mất ổn định trượt lở]
   Gia tăng trọng lượng thể tích gamma ---> Giảm tau_f về c' + (sigma - ua)tan(phi') ---> FS < 1.0 (Phá hoạ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 lý thuyết địa kỹ thuật nền tảng:

  1. Lý thuyết cơ học liên tục 4 pha của Fredlund & Morgenstern (1977): Đất gồm hạt khoáng, nước, khí và màng phân cách khí - nước (contractile skin).
  2. Mô hình cấu trúc lỗ rỗng và đặc trưng giữ nước Modified Kovacs (MK): Ước lượng SWCC thông qua đường kính hạt $D_{10}$, hệ số đồng đều $C_u$, giới hạn chảy $W_L$, giới hạn dẻo $W_P$ và tỷ trọng $G_s$.
  3. Định luật thấm Darcy mở rộng cho môi trường không bão hòa: Vận tốc dòng thấm $v_w = -k_w(\psi)\nabla h$, trong đó hàm thấm $k_w(\psi)$ là một hàm phi tuyến suy giảm hàng bậc phân bội theo lực hút dính $\psi$.

Điều kiện biên giới hạn (boundary conditions) được xác lập nghiêm ngặt: áp lực khí lỗ rỗng được giả thiết đẳng áp với khí quyển ($u_a = 0$ áp kế), giới hạn đo trực tiếp của căng kế chân không đạt $-90\text{ kPa}$, và dòng thấm bề mặt chuyển đổi linh hoạt giữa điều kiện biên lưu lượng ($q \le k_s$) và điều kiện biên cột nước áp lực ($h = 0$) khi xuất hiện dòng tràn ($RR$).

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng thực nghiệm (Positivism / Post-positivism) kết hợp chặt chẽ giữa:

  • Thí nghiệm trong phòng quy chuẩn cao (Advanced Laboratory Testing).
  • Thí nghiệm hiện trường trên công trình thực tế (In-situ Field Monitoring).
  • Mô hình vật lý kích thước lớn có kiểm soát biến số (Large-scale Physical Flume Modeling).
  • Mô phỏng số phần tử hữu hạn kết hợp cân bằng giới hạn tổng quát (Coupled Seepage-Slope Stability Numerical Modeling qua GeoStudio: SEEP/W và SLOPE/W).

Thiết kế đa tầng (multi-level design) được tiến hành qua 3 cấp độ:

  • Cấp độ vi mô / phần tử: Xác định thành phần hạt, giới hạn Atterberg, tỷ trọng $G_s$, hệ số rỗng $e$, SWCC, hàm thấm và thông số kháng cắt $(c', \phi', \phi^b)$.
  • Cấp độ mô hình vật lý: Máng thí nghiệm thép kích thước lớn ($L \times W \times H$) có khớp quay điều chỉnh hệ số mái $m = 1; 2; 4$, tích hợp dàn tạo mưa kiểm soát lưu lượng và độ đầm chặt chế bị đất $K = 0,70 \div 0,95$.
  • Cấp độ công trình thực địa: Quan trắc áp lực nước lỗ rỗng âm và phân tích mặt cắt thực tế tại 3 cụm công trình trọng điểm.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu tuân thủ các giao thức chuẩn hóa quốc tế:

  • Đo lực hút dính hiện trường: Sử dụng căng kế Jetfill 2725ARL (Soilmoisture Equipment Corp, Hoa Kỳ) với đầu gốm cao áp, đo đạc định kỳ tại đỉnh, giữa và chân mái dốc công trình đập Khau Piều, Bầu Lầy, Chúc Bài Sơn trong các đợt quan trắc mùa khô và mùa mưa năm 2017. Độ sâu đặt đầu đo từ $0,3\text{ m}$ đến $1,5\text{ m}$. Áp dụng công thức hiệu chỉnh cột nước trong ống dẫn căng kế: với chiều sâu $1,5\text{ m}$, trị số hiệu chỉnh áp lực là $\Delta u_w = 15,2\text{ kPa}$.
  • Xác định SWCC trong phòng: Sử dụng bình chiết áp lực cao (Pressure Plate Extractor) kết hợp đĩa gốm có giá trị bọt khí cao (High Air Entry Ceramic Disc - 15 bar). Mẫu đất được bão hòa hoàn toàn bằng buồng hút chân không, sau đó tăng từng cấp áp lực khí nén ($10, 20, 50, 100, 200, 400, 800, 1500\text{ kPa}$) để é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.
  • Xác định cường độ kháng cắt không bão hòa: Vận hành song song máy nén ba trục TRIAX50 (thực hiện các sơ đồ CW - Constant Water Content, CD - Consolidated Drained, CU - Consolidated Undrained) và máy cắt trực tiếp biến tốc EDJ-2 để xác lập đường bao phá hoại Mohr-Coulomb mở rộng trong không gian 3D.
  • Mô hình máng thí nghiệm mưa: Dàn phun mưa tự động hiệu chuẩn bằng hệ thống van áp kế, đo lưu lượng nước tràn bề mặt $Q_T$ (lít/phút), cường độ tràn ổn định $RR$ (mm/h) và tỷ lệ xâm nhập $RI$ (%). Sự biến thiên của lực hút dính trong thân mô hình được ghi nhận liên tục thời gian thực thông qua hệ thống Remote Tensometer 2100F và bộ chuyển đổi dữ liệu kết nối phần mềm chuyên dụng ICT.

Tam giác đạc (triangulation) được thực hiện đồng thời giữa: dữ liệu thí nghiệm bình chiết áp lực cao - mô hình ước lượng giải tích MK - kết quả đo căng kế hiện trường, triệt tiêu sai số hệ thống và bảo đảm tính vững (validity) của các hàm SWCC xây dựng được.

Data và phân tích

Tham số cơ lý / Thí nghiệ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) Đất mỏ Đại Phong (Quảng Ninh)
Thành phần hạt (%): Cát / Bột / Sét 18,2 / 46,5 / 35,3 12,4 / 51,2 / 36,4 22,1 / 44,8 / 33,1 15,6 / 48,2 / 36,2
Giới hạn Atterberg: $W_L$ / $W_P$ / $I_P$ (%) 38,5 / 21,2 / 17,3 41,2 / 22,6 / 18,6 36,8 / 20,4 / 16,4 37,4 / 21,0 / 16,4
Khối lượng thể tích khô max $\rho_{d\max}$ (g/cm³) 1,72 1,68 1,75 1,74
Độ ẩm tối ưu $w_{opt}$ (%) 16,8 18,2 15,9 16,2
Hệ số thấm bão hòa $k_s$ (cm/s) $3,2 \times 10^{-5}$ $1,8 \times 10^{-5}$ $4,5 \times 10^{-5}$ $2,5 \times 10^{-5}$
Giá trị khí vào (AEV) $\psi_a$ (kPa) 28,0 35,0 22,0 25,0
Lực dính hiệu quả $c'$ (kPa) 18,5 22,0 16,0 19,0
Góc ma sát trong hiệu quả $\phi'$ (độ) 23,5° 21,0° 25,0° 24,0°
Góc ma sát biểu kiến trung bình $\phi^b$ (độ) 13,5° 12,0° 14,5° 13,8°

Phân tích số liệu sử dụng phần mềm GeoStudio để giải bài toán ngẫu hợp thấm - biến dạng - ổn định:

  1. Mô-đun SEEP/W: Giải phương trình vi phân dòng thấm phi tuyến không dừng 2D của Richard: $$\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^2 \gamma_w \frac{\partial h}{\partial t}$$
  2. Mô-đun SLOPE/W: Truyền trường áp lực nước lỗ rỗng $u_w(x, y, t)$ từ SEEP/W sang để tính toán hệ số an toàn $FS(t)$ theo phương pháp Spencer và Morgenstern-Price, quét tìm tâm trượt nguy hiểm nhất qua từng bước thời gian $\Delta t = 0,5\text{ giờ}$.

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

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

  1. Quy luật phân phối dòng tràn và dòng thấm trên mái dốc đất đắp: Kết quả thí nghiệm mô hình máng với đất mỏ Đại Phong chứng minh rõ ràng: Tỷ lệ nước mưa thấm vào mái dốc ($RI$) chỉ chiếm từ 14% đến 34% tổng lượng mưa, phần còn lại từ 66% đến 86% chuyển thành dòng chảy tràn bề mặt ($RR$). Tỷ lệ này phụ thuộc phi tuyến vào độ đầm chặt ($K$) và độ dốc mái ($m$):

    • Khi hệ số đầm chặt tăng từ $K = 0,70$ lên $K = 0,95$, cường độ tràn ổn định $RR$ tăng từ 28,4 mm/h lên 46,2 mm/h, làm giảm tỷ lệ thấm xâm nhập $RI$.
    • Khi mái dốc càng thoải ($m$ tăng từ 1 lên 4), vận tốc chảy tràn giảm, kéo dài thời gian lưu nước trên bề mặt, làm tăng lượng nước mưa thấm sâu vào thân dốc.

  2. Cơ chế suy giảm và hiện tượng "trễ thủy lực" của lực hút dính ma trận: Quan trắc bằng mạng lưới cảm biến 2100F ghi nhận: Trong suốt quá trình mưa, lực hút dính tại tầng sát mặt ($0 \div 30\text{ cm}$) giảm đột ngột từ các giá trị ban đầu ($60 \div 80\text{ kPa}$) về tiệm cận $0\text{ kPa}$. Tuy nhiên, tại các tầng sâu ($> 80\text{ cm}$), lực hút dính duy trì độ trễ đáng kể. Sau khi dứt mưa $1 \div 3\text{ ngày}$, quá trình bốc hơi bề mặt tái lập lại lực hút dính ở lớp ngoài cùng, nhưng nước mưa tích tụ bên trong tiếp tục thấm sâu xuống dưới do trọng lực, khiến lực hút dính tại các tầng sâu tiếp tục suy giảm và đạt giá trị cực tiểu sau khi mưa đã kết thúc nhiều ngày.

  3. Sự khác biệt mang tính quyết định giữa hai kiểu mưa HI và LD đối với hệ số an toàn:

    • Kiểu mưa HI (Cường độ lớn, thời gian ngắn): Nước mưa vượt quá tốc độ thấm $k_s$, phần lớn chảy tràn bề mặt. Hệ số an toàn $FS$ suy giảm cục bộ ở bề mặt nhưng đáy cung trượt sâu hầu như không bị ảnh hưởng. $FS$ giảm nhanh trong vài giờ đầu nhưng phục hồi nhanh chóng sau khi dứt mưa.
    • Kiểu mưa LD (Cường độ nhỏ, kéo dài nhiều ngày): Toàn bộ lượng mưa thấm hết vào thân đất ($RI \to 100%$). Mặt ẩm thấm sâu xuống cung trượt nguy hiểm, triệt tiêu lực hút dính trên diện rộng. Hệ số an toàn $FS$ của đê hữu Cầu, đập Khau Piều và đập Chúc Bài Sơn suy giảm tuyến tính liên tục, đạt giá trị nguy hiểm nhất sau $48 \div 72\text{ giờ}$ và duy trì ở mức thấp kéo dài, tạo nguy cơ vỡ đập, vỡ đê tiềm tàng.

  4. Khẳng định sự sai lệch khi bỏ qua cơ học đất không bão hòa: Khi tính toán ổn định ở trạng thái ban đầu chưa mưa, việc xét đến thành phần lực hút dính $(u_a - u_w)\tan\phi^b$ giúp hệ số an toàn $FS$ thực tế cao hơn từ 15% đến 38% so với tính toán theo cơ học đất bão hòa thuần túy. Điều này giải thích thỏa đáng vì sao nhiều mái dốc tự nhiên và công trình đất đắp dốc đứng trong mùa khô vẫn đứng vững an toàn dù tính theo lý thuyết cổ điển cho kết quả $FS < 1,0$.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án bổ sung hoàn thiện hệ thống cơ sở dữ liệu đường đặc trưng giữ nước SWCC và tham số kháng cắt 3D cho các phân lớp đất dính nguồn gốc trầm tích và phong hóa tại miền Bắc Việt Nam, tạo tiền đề để chuẩn hóa các mô hình toán địa kỹ thuật trong điều kiện khí hậu nhiệt đới gió mùa.
  • Về phương pháp luận: Xác lập quy trình thí nghiệm đồng bộ từ hiện trường (căng kế Jetfill) đến phòng thí nghiệm (bình chiết áp lực cao, máng tạo mưa) có tính lặp lại cao, chuyển giao được cho các cơ sở nghiên cứu chuyên ngành.
  • Về ứng dụng thực tiễn công trình: Cung cấp công cụ ước lượng SWCC nhanh dựa trên các chỉ tiêu cơ lý phổ thông qua mô hình MK, giúp kỹ sư thiết kế không cần thực hiện các thí nghiệm đo SWCC đắt đỏ mà vẫn có đủ tham số đầu vào tin cậy cho bài toán phân tích ổn định.
  • Về chính sách và tiêu chuẩn: Khuyến nghị sửa đổi, bổ sung các Tiêu chuẩn Quốc gia (TCVN) về thiết kế đê, đập thủy lợi: bắt buộc phải đưa kịch bản mưa kéo dài (LD) và sự suy giảm lực hút dính không bão hòa vào các trường hợp tính toán kiểm tra ổn định công trình cấp đặc biệt và cấp I.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các hạn chế học thuật và rào cản kỹ thuật:

  • Giới hạn đo của thiết bị căng kế hiện trường Jetfill và Tensometer dạng cơ học/thủy lực bị chặn ở ngưỡng bọt khí $-90\text{ kPa}$. Trong điều kiện khô hạn gay gắt, lực hút dính của đất sét có thể lên tới hàng nghìn kPa nhưng chưa thể đo trực tiếp liên tục tại hiện trường.
  • Mô hình máng thí nghiệm vật lý tuy có kích thước lớn nhưng vẫn chịu hiệu ứng thành biên (side-wall boundary effect) và chưa tích hợp được ảnh hưởng của thảm phủ thực vật (rễ cây hút nước và gia cường cơ học) cũng như sự xuất hiện của các khe nứt nẻ co ngót (desiccation cracks) trong đất đắp.
  • Mới chỉ khảo sát 4 loại đất dính đại diện tại miền Bắc Việt Nam, chưa bao phủ hết các nhóm đất đỏ bazan Tây Nguyên hay đất cát ven biển miền Trung.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi trọng tâm:

  1. Nghiên cứu cơ chế nứt nẻ bề mặt do nhiệt độ cao kết hợp với mưa cực đoan (Desiccation cracking - Infiltration coupling) và sự hình thành các dòng thấm ưu tiên (preferential flow).
  2. Phát triển thiết bị đo lực hút dính dải rộng (High-capacity tensometers lên tới $1.500\text{ kPa}$) và cảm biến độ ẩm sợi quang (fiber-optic moisture sensors) chôn vĩnh cửu trong thân đập.
  3. Mở rộng nghiên cứu ảnh hưởng của hệ rễ thực vật (Bio-slope engineering) đến việc duy trì lực hút dính ma trận trong mùa mưa lũ.
  4. Xây dựng thuật toán học máy (Machine Learning) kết hợp mạng thần kinh vật lý (PINNs) để dự báo thời gian thực hệ số an toàn mái dốc từ dữ liệu mưa vệ tinh và cảm biến áp lực nước lỗ rỗng.

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

  • Giá trị học thuật: Các kết quả của luận án mở đường cho chuỗi công bố trên các tạp chí khoa học chuyên ngành uy tín trong nước (Tạp chí Khoa học & Công nghệ Thủy lợi, Tạp chí Địa kỹ thuật) và quốc tế, đóng góp một khối dữ liệu thực nghiệm quý giá cho cộng đồng cơ học đất không bão hòa khu vực Đông Nam Á.
  • Chuyển đổi ngành công nghiệp xây dựng thủy lợi: Giúp các tổng công ty tư vấn thiết kế (như HEC1, HEC2, các Viện Quy hoạch Thủy lợi) tối ưu hóa kích thước hình học thân đê, đập và mái đào kênh dẫn. Việc tính toán chính xác lực hút dính giúp tránh thiết kế quá dư thừa trong mùa khô nhưng đồng thời cảnh báo chính xác độ suy giảm an toàn trong mùa mưa, tiết kiệm hàng chục tỷ đồng chi phí xử lý gia cố mái dốc không cần thiết.
  • Ảnh hưởng chính sách quản lý thiên tai: Cung cấp cơ sở khoa học để Tổng cục Phòng chống Thiên tai, Cục Quản lý Đê điều và Phòng chống Thiên tai xây dựng các bản đồ phân vùng cảnh báo sạt lở đê điều theo ngưỡng mưa tích lũy ($P_{10}$), nâng cao năng lực ứng phó sự cố trong mùa bão lũ.

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

  • Nghiên cứu sinh và học viên cao học Địa kỹ thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh, quy chuẩn từ lý thuyết cơ học môi trường liên tục đến kỹ thuật vận hành thiết bị đo áp lực nước lỗ rỗng âm hiện đại.
  • Các nhà khoa học và giảng viên chuyên ngành: Sở hữu bộ dữ liệu thực nghiệm đối sánh (benchmark dataset) về SWCC, hàm thấm và góc ma sát biểu kiến $\phi^b$ cho các loại đất nhiệt đới Việt Nam.
  • Kỹ sư thiết kế và tư vấn Địa kỹ thuật - Thủy lợi: Ứng dụng quy trình tính toán ngẫu hợp SEEP/W và SLOPE/W, sử dụng mô hình Modified Kovacs (MK) để ước lượng nhanh các tham số không bão hòa phục vụ thẩm tra, thẩm định công trình.
  • Cơ quan quản lý nhà nước và Ban quản lý hồ đập: Có được công cụ đánh giá rủi ro trượt lở dựa trên quan trắc thời gian thực áp lực nước lỗ rỗng âm và kịch bản mưa dự báo.

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 độc đáo nhất là việc thực nghiệm hóa và chuẩn hóa định luật cường độ kháng cắt mở rộng Fredlund et al. (1978) cho đất dính đầm chặt công trình thủy lợi miền Bắc Việt Nam, chứng minh sự biến thiên phi tuyến của góc ma sát biểu kiến $\phi^b$ theo lực hút dính ma trận. Luận án đã mở rộng lý thuyết ứng suất hiệu quả đơn pha kinh điển của Terzaghi (1936) sang lý thuyết trạng thái ứng suất đa pha độc lập, định lượng hóa giá trị tăng cường lực dính $c_{add} = (u_a - u_w)\tan\phi^b$ trong việc duy trì ổn định mái dốc.

2. Điểm cải tiến mang tính đột phá về mặt phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu trước đây vốn phụ thuộc vào mô hình số đơn thuần hoặc phải gửi mẫu ra nước ngoài (như Hương, 2013 gửi mẫu sang Nhật Bản; Gasmo et al., 2000 chỉ mô phỏng số), luận án đã tự chủ thiết kế, chế tạo thành công hệ thống máng mô hình dốc trọng lực kích thước lớn tích hợp dàn mưa nhân tạo kiểm soát lưu lượng và hệ thống thu thập dữ liệu tự động Remote Tensometer 2100F. Điều này cho phép đo đạc trực tiếp sự phân tách giữa dòng tràn bề mặt ($RR$) và dòng thấm ($RI$), ghi nhận hiện tượng trễ thủy lực thời gian thực của áp lực nước lỗ rỗng âm.

3. Phát hiện bất ngờ nhất thu được từ số liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là hiện tượng trượt lở mái dốc đất đắp nguy hiểm nhất không xảy ra trong các trận mưa rào xối xả cường độ cực lớn (HI), mà diễn ra trong và sau các trận mưa dầm dề kéo dài cường độ nhỏ (LD). Số liệu chứng minh với mưa HI, lượng nước chủ yếu chảy tràn bề mặt ($RR$ đạt tới 86%), màng nước mặt đóng vai trò cản trở thấm sâu. Ngược lại, mưa LD cho phép $100%$ lượng nước mưa thẩm thấu vào cấu trúc rỗng ($RI \to 100%$), triệt tiêu hoàn toàn lực hút dính ở tầng sâu của mặt trượt tiềm năng, khiến hệ số an toàn $FS$ suy giảm sâu và không thể tự phục hồi trong nhiều ngày sau khi tạnh mưa.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp chi tiết toàn bộ quy trình chế tạo máng thí nghiệm (vật liệu, khớp xoay tự do, đồng hồ đo lưu lượng, kích thước hình học), thông số vận hành dàn mưa nhân tạo, quy trình bão hòa mẫu và nén khí từng cấp trong bình chiết áp lực cao 15 bar, cùng thuật toán hiệu chỉnh cột nước căng kế $\Delta u_w = \gamma_w h_{tube}$. Các phòng thí nghiệm Địa kỹ thuật hoàn toàn có thể nhân bản chuẩn xác các bước thực nghiệm này.

5. Lộ trình phát triển học thuật 10 năm được vạch ra như thế nào?

Lộ trình 10 năm hướng tới việc xây dựng bộ Bản đồ số quốc gia về các tham số đất không bão hòa (SWCC, $\phi^b$, $k_w$) cho toàn bộ các vùng địa chất công trình Việt Nam; tích hợp cảm biến đo lực hút dính không dây chuẩn IoT vào thân đê, đập trọng điểm; và cập nhật toàn diện Tiêu chuẩn Quốc gia về tính toán ổn định mái dốc theo hệ tiêu chuẩn tương thích Eurocode 7 mở rộng cho đất không bão hòa.

Kết luận

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

  1. Xây dựng bộ cơ sở dữ liệu thực nghiệm toàn diện về đường cong đặc trưng đất - nước (SWCC) và tham số cường độ kháng cắt không bão hòa $(c', \phi', \phi^b)$ cho các loại đất đắp tiêu biểu tại miền Bắc Việt Nam (đê hữu Cầu, đập Khau Piều, đập Chúc Bài Sơn, mỏ đất Đại Phong).
  2. Chế tạo thành công hệ thống máng mô hình nghiêng và dàn mưa nhân tạo tự động, cho phép thực nghiệm hóa quá trình tương tác thủy - cơ giữa mưa và mái dốc đất đắp có kiểm soát độ dốc và độ đầm chặt.
  3. Định lượng hóa tỷ lệ nước mưa xâm nhập ($RI$ từ 14% đến 34%) và dòng chảy tràn bề mặt ($RR$ từ 66% đến 86%), làm sáng tỏ quy luật trễ của áp lực nước lỗ rỗng âm sau khi dừng mưa.
  4. Đánh giá bản chất cơ chế phá hoại của các hình thái mưa khác nhau, chứng minh mức độ nguy hiểm vượt trội của kiểu mưa kéo dài cường độ nhỏ (LD) so với kiểu mưa thời gian ngắn cường độ lớn (HI).
  5. Chuẩn hóa và đề xuất công cụ giải tích ước lượng SWCC dựa trên mô hình Modified Kovacs (MK), tạo giải pháp kỹ thuật ứng dụng thực tế nhanh chóng, tin cậy cho ngành xây dựng thủy lợi.
  6. Ứng dụng thành công khung phân tích ngẫu hợp SEEP/W và SLOPE/W để đánh giá định lượng hệ số an toàn cho 3 cụm công trình đê, đập thực tế, chứng minh sự gia tăng an toàn tự nhiên từ 15% đến 38% nhờ lực hút dính ma trận trong mùa khô và cảnh báo nguy cơ sạt trượt chính xác trong mùa mưa.

Công trình tạo bước chuyển dịch mô thức (paradigm shift) từ cơ học đất bão hòa truyền thống sang cơ học đất không bão hòa hiện đại trong thực tế kỹ thuật công trình tại Việt Nam, mở ra các hướng nghiên cứu liên ngành bền vững và để lại giá trị ứng dụng lâu dài cho công tác phòng chống giảm nhẹ thiên tai trượt lở đất trên phạm vi cả nước.