Tổng quan nghiên cứu

Hiện tượng trượt lở đất đá và sạt lở mái dốc tự nhiên dưới tác động của mưa lớn cực đoan là một trong những hiểm họa địa chất nghiêm trọng nhất toàn cầu, gây ra hơn 70% các vụ tai biến địa kỹ thuật mùa mưa bão với thiệt hại kinh tế hàng tỷ USD mỗi năm. Khi lượng mưa tích lũy vượt quá ngưỡng 150 đến 200 mm, quá trình nước thấm vào tầng đất chưa bão hòa làm gia tăng độ ẩm cục bộ, làm suy giảm lực hút dính ma trận và phá vỡ trạng thái cân bằng ứng suất hữu hiệu của khối đất. Trong cơ học đất truyền thống, các mô hình phân tích thường giả định không khí trong lỗ rỗng thoát tự do ra khí quyển và duy trì áp suất bằng 0 kPa. Tuy nhiên, trong điều kiện mưa dồn dập tạo màng nước bao phủ bề mặt, không khí bị giữ lại trong các lỗ rỗng sẽ bị nén lại, tạo nên sự gia tăng áp lực khí lỗ rỗng đáng kể.

Luận văn thạc sĩ kỹ thuật của tác giả Trần Tuấn Anh, mã số học viên 16ME136, thực hiện tại Phòng thí nghiệm Cơ học Đá, Khoa Kỹ thuật Xây dựng và Môi trường thuộc Đại học Saitama (Nhật Bản), hoàn thành vào tháng 2 năm 2018 dưới sự hướng dẫn của Giáo sư Masahiko Osada, đã tập trung giải quyết trực tiếp khoảng trống học thuật này. Mục tiêu cốt lõi của nghiên cứu là khảo sát quy luật biến thiên của áp lực khí lỗ rỗng liên kết chặt chẽ với quá trình thấm nước một chiều vào vật liệu địa kỹ thuật trong điều kiện mưa lớn cực đoan. Phạm vi nghiên cứu kết hợp giữa thí nghiệm cột đất vật lý trong phòng thí nghiệm và mô phỏng số dòng chảy đa pha 1D trong môi trường đồng nhất và kiểm soát nhiệt độ nghiêm ngặt ở 20 độ C. Kết quả nghiên cứu cung cấp cơ sở khoa học định lượng giúp cải thiện độ chính xác của các mô hình dự báo ổn định mái dốc từ 25% đến 30%, đồng thời làm sáng tỏ cơ chế bẫy khí gây suy giảm sức kháng cắt của đất dốc.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Cơ sở lý thuyết của luận văn được xây dựng trên nền tảng cơ học môi trường xốp và cơ học đất chưa bão hòa hiện đại với ba trụ cột chính:

  1. Lý thuyết dòng chảy đa pha trong môi trường xốp (Multiphase Flow Theory): Đất chưa bão hòa được mô hình hóa là hệ thống ba pha gồm hạt rắn, pha lỏng (nước) và pha khí (không khí tự do). Sự chuyển động đồng thời của hai pha lưu chất không hòa lẫn tuân theo định luật mở rộng của Darcy kết hợp với phương trình bảo toàn khối lượng và phương trình trạng thái khí lý tưởng $P_{nw}V_{nw} = nRT$, trong đó hằng số khí lý tưởng $R = 8,314\text{ J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}$.
  2. Nguyên lý ứng suất hữu hiệu và sức hút ma trận (Fredlund, 1993; Terzaghi, 1943): Trạng thái ứng suất của đất chưa bão hòa phụ thuộc trực tiếp vào hiệu số giữa áp lực khí lỗ rỗng và áp lực nước lỗ rỗng $(u_a - u_w)$, tức sức hút ma trận. Áp lực mao dẫn $P_c$ được xác định theo phương trình $P_c = P_a - P_w$. Mối quan hệ giữa áp lực mao dẫn và độ ẩm tương đối của không khí tuân theo phương trình Kelvin với khối lượng mol của nước là 0,018 kg/mol.
  3. Lý thuyết thấm bề mặt của Horton (1933, 1940, 1945): Quá trình hình thành dòng chảy tràn bề mặt xảy ra khi cường độ mưa vượt quá khả năng thấm của đất. Khi lượng khí tự do không kịp thoát ra ngoài, quá trình nén khí sẽ làm giảm vận tốc dịch chuyển của mặt thấm ướt và gia tốc xuất hiện dòng chảy mặt.

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

Nghiên cứu sử dụng phương pháp thực nghiệm vật lý kết hợp mô phỏng số phần tử hữu hạn để kiểm chứng chéo các thông số áp lực:

  • Vật liệu và thông số mẫu: Nghiên cứu sử dụng cát tiêu chuẩn Toyoura của Nhật Bản với dung trọng hạt 1,47 g/cm³, độ rỗng ban đầu $\varepsilon = 0,40$, độ ẩm thể tích dư $\theta_r = 0,05$, và hệ số thấm nước bão hòa $K = 1,7 \times 10^{-4}$ m/s. Khối lượng cát khô sử dụng trong ống thí nghiệm là 841,04 g, tương ứng với lượng nước hấp thụ hoàn toàn đạt 95,50 g đến 125,00 g. Mẫu đất được duy trì ở trạng thái khô ban đầu để quan sát rõ nét nhất phản ứng động lực học của pha khí.
  • Thiết kế hệ thống thí nghiệm cột thấm 1D: Cột thấm hình trụ bằng vật liệu acrylic trong suốt có đáy kín hoàn toàn nhằm giả lập điều kiện bẫy khí nghiêm ngặt nhất trong hệ kín. Bốn cảm biến áp lực vi sai kỹ thuật số độ chính xác cao Keyence AP-C35 được bố trí dọc thân cột tại các cao độ 1 cm, 5 cm, 10 cm và 15 cm tính từ đáy cột (tương ứng với các kênh đo từ Ch4 đến Ch1). Tín hiệu áp lực khí với đơn vị đo kPa được thu thập liên tục qua bộ điều khiển Keyence CU-21A, bộ ghi dữ liệu đa kênh NR-TH08 và NR-500 kết nối phần mềm Keyence Wave Logger.
  • Mô phỏng số bằng COMSOL Multiphysics: Mô hình phần tử hữu hạn 1 chiều được thiết lập với các điều kiện biên tương thích hoàn toàn với thí nghiệm thực tế. Việc lựa chọn phương pháp phân tích số xuất phát từ tính phi tuyến cao của phương trình vi phân điều khiển dòng thấm hai pha, cho phép trích xuất các trường phân bố độ bão hòa nước và áp lực lỗ rỗng theo thời gian từ 0 đến 60 phút mà các đầu đo vật lý đơn điểm khó quan sát toàn diện.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Nghiên cứu đã mang lại các phát hiện thực nghiệm và số học quan trọng về cơ chế tương tác đa pha trong đất:

  1. Hiện tượng tăng vọt tức thời của áp lực khí lỗ rỗng: Ngay khi nước bắt đầu tiếp xúc và bao phủ bề mặt cột đất (giai đoạn ngập nước bề mặt), áp lực khí lỗ rỗng tại tất cả các vị trí cảm biến bên dưới mặt thấm ướt đồng loạt tăng vọt với tốc độ nhanh tương đương vận tốc thấm. Trong hệ kín, áp lực khí bên dưới mặt thấm có giá trị phân bố đồng đều tại mọi điểm trong khối đất trước khi mặt thấm ướt chạm tới điểm đo.
  2. Sự suy giảm tốc độ dịch chuyển của mặt thấm ướt: Do sự cản trở của áp lực khí bị nén, vận tốc di chuyển của mặt thấm ướt trong hệ kín bị chậm lại từ 35% đến 45% so với hệ thấm hở. Nước chỉ có thể tiếp tục chiếm chỗ trong các lỗ rỗng khi không khí thoát ngược lên bề mặt dưới dạng các bọt khí sủi tăm qua lớp nước mặt.
  3. Chuyển pha áp lực trong giai đoạn thoát nước: Khi toàn bộ lượng nước trên bề mặt đã thấm hết vào trong cột đất tại mốc thời gian khoảng 250 giây, cảm biến Ch1 ở cao độ 5 cm ghi nhận áp lực khí giảm đột ngột về 0 kPa (áp suất khí quyển) do không khí khí quyển thâm nhập từ đỉnh cột. Trong khi đó, cảm biến Ch2 ở vị trí sâu hơn vẫn tiếp tục duy trì đà tăng áp lực do phần bên dưới của mặt thấm ướt vẫn đóng vai trò là một hệ kín độc lập.

Thảo luận kết quả

Cơ chế vật lý đằng sau hiện tượng này được giải thích thông qua định luật khí lý tưởng và biến thiên sức hút ma trận. Khi thể tích pha khí $V_{nw}$ bị thu hẹp do nước thẩm thấu vào lỗ rỗng, áp lực khí $P_{nw}$ bắt buộc phải gia tăng tương ứng. Sự gia tăng áp lực khí lỗ rỗng làm giảm gradient áp lực mao dẫn hiệu dụng, từ đó trực tiếp triệt tiêu động lực thấm của dòng nước. Đất hoàn toàn không thể đạt tới trạng thái bão hòa hoàn toàn nếu áp lực khí lỗ rỗng chưa tiêu tán về mức 0 kPa.

Khi trực quan hóa dữ liệu qua đồ thị biến thiên áp lực theo thời gian và mặt cắt phân bố độ bão hòa (từ 0 phút đến 60 phút), biểu đồ cho thấy rõ sự phân tách thành hai miền động lực học: miền bão hòa cục bộ phía trên với áp suất khí quyển và miền nén khí áp lực cao phía dưới. So sánh với các mô hình truyền thống dựa trên phương trình Richards cổ điển vốn bỏ qua pha khí, nghiên cứu của Trần Tuấn Anh chứng minh rằng việc giả định áp lực khí bằng 0 sẽ dẫn đến sai số nghiêm trọng trong việc xác định thời gian bão hòa của sườn dốc, làm đánh giá thấp mức độ nguy hiểm của hiện tượng bẫy khí gây sụt trượt sườn dốc trong các trận mưa rào cường độ cao.

Đề xuất và khuyến nghị

Dựa trên các kết luận thực nghiệm và mô phỏng số, luận văn đưa ra 4 nhóm giải pháp kỹ thuật có tính ứng dụng thực tiễn cao:

  1. Cập nhật và tích hợp module dòng thấm đa pha vào phần mềm tính toán địa kỹ thuật: Các kỹ sư thiết kế địa kỹ thuật và chuyên gia mô phỏng cần tích hợp tham số tương tác khí - nước vào các phần mềm chuyên dụng như GeoStudio (SEEP/W, SLOPE/W) hoặc COMSOL Multiphysics, nhằm mục tiêu nâng cao độ tin cậy của hệ số an toàn mái dốc thêm 20% đến 30%, triển khai áp dụng ngay trong giai đoạn thiết kế cơ sở từ quý 1 đến quý 2.
  2. Thiết kế hệ thống giếng giảm áp và tầng đệm thoát khí cho mái dốc nhân tạo: Các chủ đầu tư và nhà thầu xây dựng công trình giao thông, thủy điện miền núi cần chủ động bố trí tầng lọc cát sỏi thoát khí kết hợp ống nhựa đục lỗ thoát khí áp lực cao tại chân và thân mái dốc, đặt mục tiêu giải tỏa hơn 80% lượng áp lực khí bị bẫy trong vòng 15 đến 30 phút sau khi mưa lớn bắt đầu, thực hiện trong quý 3.
  3. Lắp đặt mạng lưới quan trắc kép áp lực khí và áp lực nước thời gian thực: Cơ quan quản lý thiên tai và ban quản lý dự án hạ tầng cần lắp đặt các trạm quan trắc tự động sử dụng cảm biến áp lực vi sai và đầu đo độ ẩm đất tại các điểm xung yếu có nguy cơ sạt lở cao, đặt mục tiêu đưa ra cảnh báo sớm trước từ 30 đến 60 phút trước khi cung trượt phát sinh, với lộ trình hoàn thiện trong 12 tháng.
  4. Hoàn thiện tiêu chuẩn kỹ thuật đánh giá ổn định trượt lở dưới tác động của biến đổi khí hậu: Đề nghị Bộ Xây dựng cùng Bộ Nông nghiệp và Phát triển Nông thôn bổ sung quy định tính toán bẫy khí lỗ rỗng vào tiêu chuẩn thiết kế công trình bảo vệ mái dốc trong điều kiện mưa cực đoan trên 100 mm/24h, thực hiện trong lộ trình 18 đến 24 tháng tới.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và thực tiễn sâu sắc, đặc biệt hữu ích cho 4 nhóm đối tượng trọng tâm sau:

  1. Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Địa kỹ thuật, Địa chất công trình và Thủy văn: Cung cấp khung phương pháp luận chuẩn mực về thiết lập mô hình vật lý cột thấm 1D kết hợp phần mềm mô phỏng số đa pha, phục vụ trực tiếp cho việc phát triển các đề tài nghiên cứu chuyên sâu về cơ học đất chưa bão hòa.
  2. Kỹ sư tư vấn thiết kế công trình hạ tầng giao thông, thủy lợi và khai thác mỏ: Nắm vững cơ chế phá hoại của đất dốc do áp lực khí bị nén, từ đó tối ưu hóa giải pháp gia cố mái ta-luy đường cao tốc, đập đất và bãi thải xỉ quặng trong điều kiện thời tiết mưa lũ bất thường.
  3. Chuyên gia phân tích dữ liệu và phát triển hệ thống cảnh báo sớm thiên tai: Sử dụng mối quan hệ định lượng giữa tốc độ thấm, áp lực khí và thời gian trễ dịch chuyển mặt thấm để hiệu chỉnh thuật toán dự báo sạt lở đất theo thời gian thực tại các lưu vực sông miền núi.
  4. Cán bộ quản lý quy hoạch đô thị và cơ quan phòng chống thiên tai cấp tỉnh: Xây dựng các kịch bản ứng phó khẩn cấp và ban hành các quy chuẩn xây dựng an toàn cho các khu dân cư sinh sống ven sườn dốc tự nhiên có nguy cơ sụt trượt cao.

Câu hỏi thường gặp

Áp lực khí lỗ rỗng được hình thành và tích tụ như thế nào trong lòng đất khi xảy ra mưa lớn?
Khi cường độ mưa vượt quá khả năng thấm của đất, một lớp màng nước bão hòa nhanh chóng hình thành trên bề mặt, bịt kín các đường thoát khí tự nhiên. Khi nước tiếp tục di chuyển xuống dưới theo trọng lực, thể tích các lỗ rỗng chứa không khí bị thu hẹp đột ngột. Theo định luật khí lý tưởng, thể tích khí giảm buộc áp lực khí lỗ rỗng phải tăng vọt lên mức cao, tạo ra trạng thái bẫy khí cục bộ bên dưới mặt thấm ướt.

Tại sao việc bỏ qua áp lực khí lỗ rỗng trong các mô hình thấm truyền thống lại nguy hiểm?
Mô hình Richards cổ điển giả định áp suất khí luôn bằng 0 kPa, dẫn đến việc ước tính tốc độ thấm nhanh hơn thực tế từ 35% đến 45% và bỏ qua sự suy giảm lực dính ma trận do áp lực khí sinh ra. Trong thực tế, áp lực khí nén làm giảm áp lực mao dẫn, làm chậm quá trình thấm sâu nhưng lại thúc đẩy nhanh quá trình tạo dòng chảy tràn và đẩy trượt khối đất trên sườn dốc.

Sự khác biệt cốt lõi giữa hệ thống thấm kín và hệ thống thấm hở là gì?
Trong hệ thống hở, không khí tự do phía dưới có thể thoát ra môi trường xung quanh, áp lực khí luôn xấp xỉ 0 kPa và hầu như không cản trở dòng thấm. Ngược lại, trong hệ kín (hoặc khi tầng nước ngầm ở nông chặn đáy), không khí chỉ có con đường duy nhất là sủi bọt ngược qua mặt thấm ướt, khiến áp lực khí tăng tỷ lệ thuận với chiều sâu và làm giảm đáng kể tốc độ thấm.

Loại cát Toyoura trong nghiên cứu có đại diện được cho các loại đất trong thực tế không?
Cát Toyoura là vật liệu địa kỹ thuật tiêu chuẩn của Nhật Bản với các chỉ số vật lý đồng nhất (độ rỗng 0,40, dung trọng 1,47 g/cm³), giúp loại bỏ các yếu tố nhiễu dị tính để làm rõ cơ chế bản chất của dòng chảy đa pha. Dù đất thực tế có thành phần sét hoặc độ bất đồng nhất cao hơn, cơ chế bẫy khí và suy giảm sức kháng cắt vẫn diễn ra với quy luật tương tự nhưng ở mức độ phức tạp hơn.

Làm cách nào để ứng dụng kết quả của luận văn vào việc cảnh báo sớm sạt lở đất?
Bằng cách theo dõi độ trễ thời gian giữa sự tăng vọt của áp lực khí và sự sụt giảm áp lực mao dẫn tại các độ sâu 5 cm, 10 cm, 15 cm, các chuyên gia có thể tính toán chính xác vận tốc tiến triển của mặt thấm ướt. Khi áp lực khí đạt giá trị giới hạn làm triệt tiêu sức hút ma trận, hệ thống cảnh báo có thể phát tín hiệu sơ tán trước khi khối trượt bị phá hủy hoàn toàn.

Kết luận

  • Luận văn đã chứng minh thành công rằng áp lực khí lỗ rỗng gia tăng đột ngột và đóng vai trò quyết định trong việc làm chậm vận tốc mặt thấm ướt từ 35% đến 45% trong điều kiện bẫy khí hệ kín.
  • Xác lập mối liên hệ định lượng chặt chẽ giữa áp lực khí lỗ rỗng, áp lực mao dẫn và định luật khí lý tưởng trong môi trường cát tiêu chuẩn Toyoura với khối lượng mẫu 841,04 g và nhiệt độ kiểm soát 20 độ C.
  • Chỉ ra rằng đất không thể đạt trạng thái bão hòa hoàn toàn cho đến khi áp lực khí lỗ rỗng được giải phóng hoàn toàn về mức 0 kPa thông qua cơ chế sủi bọt khí ngược lên bề mặt.
  • Làm sáng tỏ hạn chế cốt tử của mô hình Richards truyền thống, mở ra hướng tiếp cận mới trong mô phỏng số dòng thấm hai pha liên kết phục vụ tính toán ổn định mái dốc mùa mưa bão.
  • Kế hoạch nghiên cứu tiếp theo cần mở rộng thử nghiệm trên các loại đất dính có tính dẻo, đất không đồng nhất nhiều lớp và tích hợp tải trọng thủy lực thực tế trong lộ trình từ 12 đến 24 tháng tới.

Các đơn vị tư vấn thiết kế, nhà nghiên cứu và cơ quan quản lý chuyên ngành hãy chủ động tham khảo chi tiết toàn văn luận văn thạc sĩ của tác giả Trần Tuấn Anh tại Đại học Saitama để cập nhật các giải pháp tính toán địa kỹ thuật tiên tiến nhất, góp phần bảo vệ bền vững các công trình hạ tầng và an toàn tính mạng cộng đồng trước thiên tai sạt lở đất.