Tổng quan nghiên cứu

Sự cố vỡ đập đất do tác động của thiên tai và mưa bão cực đoan luôn để lại những thảm họa khốc liệt về người và tài sản. Trong lịch sử kỹ thuật thủy lợi, thảm kịch vỡ đập Bản Kiều tại Trung Quốc năm 1975 sau chuỗi ngày mưa lớn liên tục đã cướp đi sinh mạng của hàng chục nghìn người. Tại Việt Nam, sự cố vỡ đập Đầm Hà Động tại tỉnh Quảng Ninh vào ngày 29/10/2014 đã gây thiệt hại kinh tế ước tính khoảng 380 tỷ đồng, làm ngập úng hơn 3.500 ha đất canh tác và khiến 115 hộ dân phải di dời khẩn cấp. Tiếp đó, sự cố trượt lở mái đập Trường Lâm tại Thanh Hóa vào tháng 10/2013 cũng xuất phát từ việc nước mưa thấm sâu vào mái hạ lưu kết hợp lưu lượng nước hồ dâng cao.

Khu vực Nam Trung Bộ và Đông Nam Bộ hiện quản lý hơn 200 hồ chứa thủy lợi lớn nhỏ được đắp bằng vật liệu đất địa phương trong hơn 25 năm qua. Tuy nhiên, phần lớn các công trình này đang đối mặt với tình trạng xuống cấp, hệ thống tiêu thoát nước bị bồi lắng hoặc tắc nghẽn, dẫn đến nguy cơ mất ổn định nghiêm trọng khi xảy ra các đợt mưa kéo dài từ 2 đến 3 ngày.

Mục tiêu trọng tâm của nghiên cứu là xác lập cơ sở khoa học về động thái dòng thấm không ổn định, đánh giá mức độ suy giảm cường độ kháng cắt của đất đắp dưới tác động của mưa lớn kéo dài, và ứng dụng tính toán định lượng cho công trình hồ chứa nước Diên Trường, tỉnh Quảng Ngãi. Luận văn sử dụng công cụ mô phỏng số hiện đại nhằm xây dựng các kịch bản thấm bất lợi, từ đó đề xuất giải pháp kết cấu gia cố tối ưu nhằm nâng cao hệ số an toàn ổn định công trình từ 15% đến 30%, đảm bảo an toàn tuyệt đối cho vùng hạ du.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng cơ học đất và thủy lực ngầm hiện đại, tích hợp hai lý thuyết nền tảng cốt lõi:

  • Lý thuyết thấm bão hòa và không bão hòa: Vận tốc dòng thấm trong môi trường đất được thiết lập theo định luật thực nghiệm Darcy với công thức $v = k \cdot J$, trong đó $k$ là hệ số thấm và $J$ là gradient thủy lực. Quá trình dịch chuyển của nước mưa qua thân đập được mô tả tổng quát bằng phương pháp phần tử hữu hạn Galerkin giải phương trình vi phân dòng thấm không ổn định 2D:

$$\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 = S_s \frac{\partial h}{\partial t}$$

Trong đó $K_x, K_y$ lần lượt là hệ số thấm theo các phương tọa độ, $h$ là cột áp thủy lực tổng, $q$ là lưu lượng nguồn bổ cập từ mưa và $S_s$ là hệ số trữ nước riêng của tầng đất.

  • Lý thuyết cân bằng giới hạn và tiêu chuẩn phá hoại Mohr-Coulomb: Cường độ chống cắt huy động dọc theo mặt trượt tiềm năng được xác định qua hàm quan hệ: $\tau = c' + (\sigma_n - u)\tan\phi'$, với $c'$ là lực dính hiệu quả, $\phi'$ là góc ma sát trong hiệu quả, $\sigma_n$ là ứng suất pháp tổng và $u$ là áp lực nước lỗ rỗng. Khi mưa thấm sâu vào mái dốc, đại lượng $u$ chuyển từ giá trị âm sang giá trị dương, làm suy giảm trực tiếp ứng suất pháp hữu hiệu.

Nghiên cứu vận dụng đồng thời 4 khái niệm chuyên ngành then chốt: đường bão hòa thấm, áp lực nước lỗ rỗng dư, hệ số an toàn ổn định trượt mái dốc ($F_s$), và cơ chế thoát nước bằng tầng lọc ngược trọng lực.

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

Nghiên cứu sử dụng nguồn dữ liệu khảo sát địa kỹ thuật chi tiết tại đập Diên Trường bao gồm chỉ tiêu cơ lý mẫu đất đắp (dung trọng khô $\gamma = 16,8\text{ kN/m}^3$, lực dính $c = 18,5\text{ kPa}$, góc ma sát $\phi = 19^\circ$, hệ số thấm $k = 1,2 \times 10^{-5}\text{ cm/s}$), kết hợp chuỗi dữ liệu thủy văn lượng mưa đo đạc thực tế tại Quảng Ngãi giai đoạn 2010 - 2015.

Quy trình phân tích được thực hiện trên mặt cắt ngang thân đập điển hình với chiều cao lớn nhất qua 5 bước nghiêm ngặt:

  1. Rời rạc hóa miền tính toán thành lưới phần tử hữu hạn gồm hơn 2.400 phần tử tam giác và tứ giác trong module SEEP/W.
  2. Mô phỏng trường thấm chuyển tiếp dưới các trận mưa liên tục 24 giờ, 48 giờ và 72 giờ với cường độ mưa giả định từ 50 mm/ngày đến trên 200 mm/ngày.
  3. Truyền trường áp lực nước lỗ rỗng theo từng bước thời gian sang module SLOPE/W.
  4. Tính toán tìm kiếm cung trượt tròn nguy hiểm nhất theo các phương pháp phân thỏi Bishop rút gọn, Janbu và GLE (General Limit Equilibrium).
  5. Phương pháp chọn mẫu mặt cắt bất lợi kết hợp phương pháp số được lựa chọn thay thế cho giải tích cổ điển vì phản ánh chính xác tính phi đồng chất của các lớp vật liệu đắp và trạng thái thấm biến đổi theo thời gian thực.

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

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

Kết quả mô phỏng số trên mô hình đập Diên Trường phản ánh rõ rệt các quy luật cơ học khi xảy ra mưa kéo dài:

  • Sự biến động dâng cao của đường bão hòa: Sau 48 giờ đến 72 giờ mưa liên tục với lưu lượng thấm bổ cập $q = 1,0 \times 10^{-5}\text{ m}^3\text{/s/m}$, mặt thoáng dòng thấm bên trong thân đập dâng cao cục bộ từ 18% đến 35% so với trạng thái vận hành bình thường. Vùng bão hòa mở rộng bao trùm phần lớn khối lăng trụ mái hạ lưu, biến vùng đất chưa bão hòa thành đất no nước hoàn toàn.
  • Sự sụt giảm đột biến của hệ số an toàn ổn định ($F_s$): Tại mái hạ lưu, hệ số ổn định trượt $F_s$ giảm từ 1,52 (trạng thái khô ráo bình thường) xuống còn 1,08 sau 72 giờ mưa lớn, tương đương mức suy giảm 28,9%. Trị số này thấp hơn giới hạn an toàn cho phép $[F] = 1,15$ đến $1,25$ theo quy chuẩn công trình thủy lợi cấp III, đưa thân đập vào trạng thái cảnh báo nguy hiểm.
  • Xác định tâm trượt và mặt trượt nguy hiểm nhất: Mặt trượt giới hạn có xu hướng chuyển đổi từ dạng trượt nông bề mặt sang cung trượt sâu ăn sâu vào khối trung tâm thân đập, với diện tích khối trượt tăng 42% khi áp lực nước lỗ rỗng đạt đỉnh.
  • Hiệu quả vượt trội của kết cấu cải tạo: Khi bổ sung thiết bị thoát nước dạng ống khói kết hợp gối phẳng cát sỏi và thay thế khối đắp hạ lưu bằng vật liệu có hệ số thấm cao ($k > 1,0 \times 10^{-3}\text{ cm/s}$), đường bão hòa bị kéo tụt xuống sát đáy nền, giúp hệ số ổn định $F_s$ tăng trở lại mức 1,46 (tăng 35,1% so với hiện trạng hư hỏng).

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến mất ổn định mái dốc là do nước mưa ngấm vào thân đập đã triệt tiêu hoàn toàn lực hút dính chân không trong vùng đất không bão hòa, đồng thời làm tăng nhanh trọng lượng thể tích của khối đất từ trạng thái dung trọng tự nhiên lên dung trọng bão hòa ($\gamma_{bh} \approx 20,2\text{ kN/m}^3$). Khi áp lực nước lỗ rỗng $u$ tăng cao, sức kháng cắt của đất giảm mạnh theo định luật Mohr-Coulomb, làm mất cân bằng tĩnh học giữa mô-men gây trượt và mô-men chống trượt.

Hiện tượng này hoàn toàn tương đồng với cơ chế phá hoại được ghi nhận tại sự cố vỡ đập Đầm Hà Động năm 2014 và đập Trường Lâm năm 2013, nơi hệ thống thoát nước mái hạ lưu bị tê liệt khiến dòng thấm thoát lộ thiên trên mặt mái.

Toàn bộ diễn biến suy giảm hệ số $F_s$ theo thời gian mưa ($t = 0\text{h}, 24\text{h}, 48\text{h}, 72\text{h}$) có thể được trực quan hóa qua biểu đồ đường suy giảm phi tuyến, kết hợp bảng ma trận so sánh các phương án thoát nước giúp các kỹ sư dễ dàng lựa chọn giải pháp thiết kế tối ưu nhất.

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

Dựa trên các kết quả phân tích định lượng, luận văn đề xuất 4 nhóm giải pháp kỹ thuật cụ thể nhằm gia cố an toàn hồ đập:

  • Cải tạo và lắp đặt bổ sung hệ thống thoát nước kiểu ống khói kết hợp gối phẳng lọc ngược: Tiến hành đào bóc cục bộ và thay thế khối đắp số 3 ở lăng trụ hạ lưu bằng vật liệu đá dăm, sỏi sạn có hệ số thấm $k \ge 1,0 \times 10^{-2}\text{ cm/s}$. Giải pháp này giúp hạ thấp đường bão hòa ngay tại chân đập từ 40% đến 50%, do Ban Quản lý Đầu tư & Xây dựng Thủy lợi chủ trì thực hiện trong thời gian 6 đến 12 tháng.
  • Lắp đặt hệ thống quan trắc áp lực nước lỗ rỗng và cảnh báo tự động: Bố trí tối thiểu 6 cụm ống đo áp (piezometer) tại các vị trí mặt cắt xung yếu của đập Diên Trường, thiết lập tần suất đo tự động 15 phút/lần khi lượng mưa vượt ngưỡng 100 mm/24h. Do Công ty Khai thác Công trình Thủy lợi Quảng Ngãi vận hành từ quý II/2026.
  • Rà soát và kiểm toán lại hệ số an toàn ổn định đập theo kịch bản mưa cực đoan: Ứng dụng bộ phần mềm GeoStudio (SEEP/W và SLOPE/W) để đánh giá lại 100% các đập đất xung yếu tại miền Trung với chu kỳ mưa lũ 50 năm và 100 năm một lần. Do các Viện Nghiên cứu và Công ty Tư vấn Thiết kế Thủy lợi hoàn thành trong giai đoạn 2026 - 2028.
  • Quy chuẩn hóa quy trình duy tu, phát quang và tiêu thoát nước mặt định kỳ: Xử lý triệt để tình trạng bồi lấp rãnh thoát nước cơ mái và gia cố đá lát mái thượng lưu trước ngày 15 tháng 8 hàng năm, nhằm giảm thiểu trên 85% nguy cơ xói lở rãnh và tích tụ nước mặt trên thân đập do Chi cục Thủy lợi địa phương đôn đốc kiểm tra.

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

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng chính:

  • Kỹ sư tư vấn thiết kế và thẩm tra công trình thủy lợi - địa kỹ thuật: Nắm bắt quy trình thiết lập mô hình tính thấm không ổn định và kiểm toán ổn định mái dốc bằng phương pháp số, hỗ trợ tính toán chính xác kích thước tầng lọc và thiết bị tiêu nước cho các dự án đắp đập mới.
  • Ban Quản lý dự án và Đơn vị khai thác công trình thủy lợi: Ứng dụng quy trình đánh giá các kịch bản mưa nguy hiểm, phục vụ công tác lập phương án phòng chống thiên tai, bảo đảm an toàn hồ chứa trong mùa mưa bão tại địa phương.
  • Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Xây dựng Thủy lợi: Sử dụng làm tài liệu giảng dạy, học tập chuyên khảo về cơ chế tương tác giữa trường thấm và trường biến dạng - ứng suất trong môi trường đất đá rời rạc.
  • Cơ quan quản lý nhà nước về phòng chống thiên tai (Cục Thủy lợi, Sở Nông nghiệp và PTNT): Cung cấp cơ sở khoa học xác thực để rà soát, bổ sung và hoàn thiện các tiêu chuẩn kỹ thuật, quy chuẩn quốc gia về an toàn đập đất chịu tác động của biến đổi khí hậu.

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

Tại sao mưa kéo dài lại gây nguy hiểm cho đập đất hơn các trận mưa rào ngắn hạn?
Mưa ngắn hạn chủ yếu tạo dòng chảy tràn trên bề mặt. Ngược lại, mưa kéo dài từ 48 đến 72 giờ làm nước ngấm sâu vào thân đập, làm bão hòa lớp đất mặt dốc hạ lưu, triệt tiêu lực hút dính ma sát và gia tăng áp lực nước lỗ rỗng, khiến hệ số ổn định $F_s$ sụt giảm từ 20% đến 30%.

Module SEEP/W và SLOPE/W phối hợp xử lý bài toán ổn định thấm như thế nào?
Module SEEP/W sử dụng phương pháp phần tử hữu hạn để giải phương trình dòng thấm không ổn định theo thời gian thực, xác định trường áp lực nước lỗ rỗng $u$. Dữ liệu $u$ sau đó được tự động liên kết sang SLOPE/W để tính toán sức kháng cắt Mohr-Coulomb và tìm kiếm cung trượt nguy hiểm nhất.

Hệ thống thoát nước dạng ống khói (chimney drain) có ưu thế gì so với lăng trụ truyền thống?
Ống khói kết hợp gối phẳng cắt đứt hoàn toàn dòng thấm bên trong thân đập trước khi thoát ra mái dốc, giữ cho toàn bộ khối đắp hạ lưu luôn ở trạng thái khô ráo, giúp nâng cao hệ số ổn định công trình từ 25% đến 35% so với chỉ dùng lăng trụ chân đập.

Hệ số an toàn ổn định cho phép $[F]$ của đập đất cấp III được quy định bao nhiêu?
Theo tiêu chuẩn kỹ thuật quốc gia, đập đất cấp III với tổ hợp tải trọng cơ bản (mực nước dâng bình thường, có thấm ổn định) yêu cầu hệ số $[F] \ge 1,25$. Đối với tổ hợp tải trọng đặc biệt (mực nước lũ gia cường kết hợp mưa lớn), hệ số cho phép yêu cầu tối thiểu $[F] \ge 1,15$.

Kết quả tính toán tại đập Diên Trường có thể áp dụng cho các hồ chứa khác không?
Phương pháp luận mô hình hóa số và giải pháp kết cấu thoát nước hoàn toàn có thể nhân rộng cho hơn 200 hồ chứa thủy lợi tại vùng duyên hải miền Trung có điều kiện địa chất, vật liệu đắp đập và chế độ khí hậu mưa bão tương đồng.

Kết luận

  • Luận văn đã làm sáng tỏ cơ chế vật lý và cơ học của quá trình suy giảm ổn định mái đập đất dưới tác động của các trận mưa kéo dài từ 2 đến 3 ngày.
  • Chứng minh định lượng rằng áp lực nước lỗ rỗng tăng cao trong quá trình mưa làm suy giảm tới 28,9% hệ số an toàn ổn định $F_s$ của mái hạ lưu đập Diên Trường.
  • Khẳng định tính ưu việt của phương pháp phần tử hữu hạn qua bộ phần mềm GeoStudio (SEEP/W và SLOPE/W) trong việc giải quyết bài toán thấm không ổn định.
  • Đề xuất thành công giải pháp kỹ thuật bổ sung thiết bị thoát nước kiểu ống khói kết hợp gối phẳng lọc ngược, nâng hệ số an toàn công trình lên 1,46.
  • Cung cấp cơ sở khoa học thực tiễn vững chắc phục vụ công tác thiết kế, nâng cấp và quản lý an toàn hồ đập tại khu vực miền Trung trong giai đoạn 2026 - 2030.

Nghiên cứu mở ra hướng đi tiếp theo trong việc phát triển mô hình ghép 3D tương tác nhiệt - thủy - cơ và tích hợp dữ liệu viễn thám phục vụ cảnh báo thiên tai tự động. Các đơn vị quản lý hồ đập cần chủ động áp dụng các giải pháp mô phỏng số này để kiểm toán an toàn toàn diện cho các công trình thủy lợi hiện hữu ngay trong mùa mưa bão tới.