Tổng quan nghiên cứu

Việt Nam hiện có gần 7.000 hồ chứa thủy lợi với tổng dung tích trữ nước đạt khoảng 11 tỷ m3, giữ vai trò đặc biệt quan trọng trong việc cắt giảm lũ, cấp nước tưới tiêu nông nghiệp, phát điện và phục vụ dân sinh. Trong tổng số này, có 560 hồ chứa lớn với dung tích trên 3 triệu m3, 1.752 hồ vừa từ 0,2 đến 3 triệu m3 và hơn 4.336 hồ chứa nhỏ có dung tích dưới 0,2 triệu m3. Đáng chú ý, hơn 70% các công trình hồ đập được xây dựng từ 30 đến 40 năm trước (trong các thập niên 1970 - 1980) chủ yếu bằng vật liệu đất đắp với phương tiện thi công thủ công, chất lượng đầm nén chưa đồng đều và thiếu hệ thống đo đạc quan trắc chuyên dụng.

Tác động tiêu cực của biến đổi khí hậu toàn cầu cùng sự suy giảm thảm phủ rừng đầu nguồn đã làm gia tăng các trận mưa lũ cực đoan, tạo đỉnh lũ tập trung nhanh vượt xa tần suất thiết kế ban đầu. Tình trạng bồi lắng lòng hồ kéo dài làm suy giảm nghiêm trọng dung tích phòng lũ, dẫn tới hiện tượng mực nước hồ dâng cao vượt qua cao trình đỉnh đập đất. Vấn đề nghiên cứu trọng tâm là đánh giá mức độ suy giảm ổn định của thân và mái đập đất khi xuất hiện cột nước tràn đỉnh, từ đó so sánh hiệu quả an toàn giữa giải pháp đắp bao tải đất ngăn tràn truyền thống và phương án chủ động cho nước tràn qua đỉnh đập có gia cố bảo vệ.

Mục tiêu cụ thể của đề tài là xác định định lượng sự biến thiên của hệ số an toàn chống trượt K theo chiều cao cột nước vượt đỉnh từ 0,1m đến 1,1m trên các quy mô đập chiều cao từ 5m đến 12m, ứng dụng kiểm chứng trực tiếp cho hồ chứa nước Chấn Sơn. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong quản lý an toàn đập, giúp duy trì hệ số ổn định tối thiểu K lớn hơn hoặc bằng 1,15 đến 1,30, giảm thiểu trên 85% nguy cơ xảy ra sự cố vỡ đập đột ngột và tiết kiệm hàng tỷ đồng chi phí ứng phó khẩn cấp hàng năm.

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 vận dụng lý thuyết thấm chất lỏng trong môi trường rỗng liên tục dựa trên phương trình vi phân Laplace và định luật Darcy để mô tả chuyển động của dòng thấm qua thân và nền đập đất đồng chất. Phân tích ổn định mái dốc được xây dựng trên nền tảng lý thuyết cân bằng giới hạn phân thỏi (Limit Equilibrium Method) kết hợp tiêu chuẩn bền Mohr-Coulomb theo trạng thái ứng suất hiệu quả:

$$\tau = c' + (\sigma_n - u_w)\tan\phi'$$

Trong đó, $\tau$ biểu thị cường độ chống cắt của đất đắp, $c'$ là lực dính hiệu quả, $\sigma_n$ là ứng suất pháp tuyến tổng, $u_w$ là áp lực nước lỗ rỗng và $\phi'$ là góc ma sát trong hiệu quả. Các khái niệm cốt lõi bao gồm đường bão hòa thấm, mặt trượt trụ tròn nguy hiểm, gradient thấm thủy lực và hệ số an toàn chống trượt K được xác định theo tỷ số giữa tổng mô men chống trượt và tổng mô men gây trượt quanh tâm quay giả định.

Để định vị vùng tâm trượt nguy hiểm nhất, đề tài kết hợp phương pháp lý thuyết phân tích của N. Ghecxevanop cùng sơ đồ hình học của W. Fellenius và V. Fandeep. Bán kính cung trượt R được chuẩn hóa theo tỷ lệ R/H từ 1,50 đến 5,50 tương ứng với các hệ số mái dốc từ 1:1,0 đến 1:5,0, đảm bảo khoanh vùng chính xác cung trượt có hệ số an toàn nhỏ nhất.

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 chất công trình đầu mối hồ Chấn Sơn kết hợp dữ liệu thống kê chỉ tiêu cơ lý đất đắp từ các hồ đập tại Việt Nam. Đất nền có dung trọng khô 1,58 T/m3, góc ma sát trong 12°30', lực dính đơn vị 24,81 kPa và hệ số thấm $15,7 \times 10^{-6}\text{ cm/s}$. Đất đắp đập có dung trọng khô 1,50 T/m3, độ ẩm 23,73%, lực dính đơn vị 22,12 kPa và hệ số thấm $22,12 \times 10^{-6}\text{ cm/s}$.

Quy mô mẫu tính toán bao gồm 8 mô hình hình học điển hình với chiều cao đập H lần lượt là 5m, 8m, 10m và 12m; bề rộng mặt đỉnh đập B chọn ở hai mức chuẩn là 4m và 6m; hệ số mái thượng lưu $m_1 = 0,05H + 2$, hệ số mái hạ lưu $m_2 = 0,05H + 1,5$. Cỡ mẫu được phân tích dưới 6 kịch bản cột nước tràn đỉnh h dao động từ 0,1m đến 1,1m.

Phương pháp phân tích được thực hiện bằng mô hình số phần tử hữu hạn (FEM) trên bộ phần mềm chuyên dụng Geo-Slope của Canada. Module SEEP/W đảm nhận tính toán trường thấm hai chiều, mô phỏng chính xác sự phân bố áp lực nước lỗ rỗng và đường bão hòa. Kết quả trường áp lực này được truyền trực tiếp sang module SLOPE/W để giải bài toán ổn định mái dốc bằng phương pháp Bishop đơn giản và Janbu. Việc lựa chọn phương pháp số này giúp khắc phục hoàn toàn nhược điểm của các công thức giải tích đơn giản khi giải quyết các điều kiện biên thủy lực phức tạp trong thời gian nghiên cứu 2014-2016.

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

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

Thứ nhất, trong trường hợp đắp tôn cao đỉnh đập ngăn không cho nước tràn đỉnh khi xuất hiện lũ vượt thiết kế, hệ số an toàn ổn định mái đập K có xu hướng suy giảm nhẹ từ 3,2% đến 6,8% khi cột nước dâng cao thêm từ 0,1m đến 0,5m. Cụ thể, đối với đập có chiều cao H = 5m, hệ số K duy trì ở mức cao xấp xỉ 1,85; trong khi đập cao H = 12m ghi nhận hệ số K giảm từ 1,350 xuống còn 1,302.

Thứ hai, trong kịch bản cho phép nước tràn qua đỉnh đập có gia cố bề mặt, sự hình thành dòng tràn trên đỉnh và mái hạ lưu tạo nên áp lực thủy động học và làm tăng độ bão hòa cục bộ. Hệ số an toàn chống trượt của mái đập hạ lưu giảm từ 8,5% đến 14,2% so với điều kiện làm việc ở mực nước dâng bình thường. Tuy vậy, giá trị K tính toán tại các cột nước tràn đỉnh h từ 0,1m đến 0,4m vẫn đạt từ 1,18 đến 1,42, vượt qua ngưỡng an toàn cho phép tối thiểu theo quy chuẩn hiện hành là 1,15.

Thứ ba, nghiên cứu chứng minh mối quan hệ tỷ lệ nghịch rõ nét giữa chiều cao đập và hệ số an toàn. Khi chiều cao đập tăng từ 5m lên 12m, hệ số ổn định tổng thể suy giảm trung bình 28,6% do trọng lượng bản thân khối đất trượt tăng mạnh cùng tác động của gradient áp lực thấm tại chân khay hạ lưu.

Thứ tư, tính toán thủy lực dòng tràn qua đỉnh đập xác định độ sâu phân giới $h_k$ tại mặt cắt đầu dốc hạ lưu dao động trong khoảng 0,12m đến 0,45m ứng với cột nước tràn từ 0,2m đến 0,8m. Dòng chảy chuyển sang trạng thái chảy xiết trên mái dốc hạ lưu, khẳng định tính bắt buộc phải có kết cấu gia cố chống xói bề mặt.

Thảo luận kết quả

Nguyên nhân cốt lõi khiến hệ số ổn định không bị sụt giảm đột ngột khi nước dâng cao là do chu kỳ duy trì đỉnh lũ thường diễn ra tương đối ngắn, chỉ kéo dài từ 6 đến 24 giờ. Trong khoảng thời gian này, đường bão hòa thấm mới chưa kịp dâng cao hoàn toàn qua toàn bộ khối đất đắp có hệ số thấm nhỏ ($K \approx 10^{-5}\text{ cm/s}$), áp lực nước lỗ rỗng dư thừa chủ yếu phát sinh cục bộ ở lớp biên sát mặt ngoài.

So sánh với các nghiên cứu thủy lực công trình trước đây, kết quả này củng cố luận điểm khoa học rằng việc chủ động tháo lũ tràn đỉnh qua bề mặt đập đã gia cố là một giải pháp tình thế an toàn và kiểm soát được rủi ro, thay vì để mực nước hồ dâng tự do gây vỡ đập do xói lở tràn bờ không kiểm soát.

Về mặt trực quan hóa dữ liệu, các kết quả mô phỏng được biểu diễn chi tiết qua đường cong quan hệ cột nước tràn và hệ số ổn định (đồ thị h - K), kết hợp cùng sơ đồ phân bố các đường đẳng áp lực nước lỗ rỗng và mạng lưới phần tử hữu hạn. Các bảng tính tổng hợp độ sâu dòng chảy dọc dốc nước hạ lưu cung cấp cơ sở kỹ thuật chính xác cho việc lựa chọn kích thước và độ dày kết cấu bảo vệ mái đập.

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

  1. Triển khai gia cố khẩn cấp đỉnh đập và mái hạ lưu bằng vật liệu mềm gồm vải địa kỹ thuật cường độ cao kết hợp bao tải cát hoặc bao tải đất cho hơn 4.300 hồ chứa dung tích nhỏ dưới 0,2 triệu m3. Biện pháp này cần hoàn thành trong vòng 2 đến 4 giờ trước khi đỉnh lũ xuất hiện, đảm bảo khả năng chịu lưu tốc dòng tràn đạt 2,5 đến 3,5 m/s nhằm ngăn chặn hiện tượng xói rãnh phá hoại thân đập. Đơn vị thực hiện: Ban Chỉ huy Phòng chống thiên tai cấp huyện và Ủy ban nhân dân cấp xã.
  2. Nâng cấp kiên cố hóa bề mặt đỉnh đập kết hợp bố trí đoạn tràn sự cố bằng kết cấu cứng như tấm bê tông đúc sẵn hoặc thảm bê tông túi khuôn FS cho các đập đất cấp III và IV (dung tích từ 0,2 đến 3,0 triệu m3). Mục tiêu là tăng lưu lượng xả phòng lũ phụ trợ thêm 15% đến 25%, duy trì hệ số ổn định công trình K lớn hơn hoặc bằng 1,25, thực hiện trong kế hoạch đầu tư công trung hạn 2026-2030. Đơn vị thực hiện: Sở Nông nghiệp và Phát triển Nông thôn các tỉnh.
  3. Ứng dụng giải pháp công nghệ sinh thái phủ xanh mái đập bằng thảm cỏ Vetiver có chiều dài bộ rễ bám sâu từ 1,0m đến 4,0m đối với các công trình đập thấp chiều cao dưới 8m. Giải pháp này giúp gia tăng cường độ chống cắt tầng nông thêm 20% đến 30%, chống xói mòn mặt ngoài do mưa lũ và giảm 40% chi phí xây lắp so với ốp mái bê tông truyền thống, hoàn thành gieo trồng trong thời gian 3 đến 6 tháng. Đơn vị thực hiện: Các Công ty TNHH Một thành viên Khai thác công trình thủy lợi.
  4. Xây dựng và tích hợp hệ thống cảnh báo sớm cùng mô hình điều tiết lũ theo thời gian thực tại các hồ chứa xung yếu như hồ Chấn Sơn. Tần suất cập nhật dữ liệu đo mưa và mực nước hồ tự động đạt 15 phút một lần trong mùa bão lũ, cho phép kích hoạt các phương án tràn đỉnh có kiểm soát trước tối thiểu 12 giờ. Đơn vị thực hiện: Cơ quan quản lý hồ đập phối hợp với Trung tâm Dự báo Khí tượng Thủy văn.

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

  1. Kỹ sư tư vấn thiết kế và thẩm tra công trình thủy lợi: Khai thác quy trình thiết lập bài toán thấm và ổn định bằng phần mềm SEEP/W và SLOPE/W, phục vụ kiểm tra khả năng chịu lực của đập đất trong các điều kiện biên lũ kiểm tra khắc nghiệt.
  2. Cán bộ quản lý và vận hành hồ đập tại địa phương: Vận dụng các giải pháp ứng phó khẩn cấp bằng vật liệu mềm tại chỗ (vải bạt, bao tải đất) để bảo vệ an toàn công trình đầu mối khi xảy ra lũ vượt tần suất thiết kế với cột nước tràn từ 0,1m đến 0,5m.
  3. Cơ quan quản lý nhà nước về đê điều và phòng chống thiên tai: Sử dụng các số liệu đánh giá định lượng và biểu đồ quan hệ h - K làm căn cứ khoa học để hoàn thiện hệ thống văn bản quy phạm pháp luật, quy trình kiểm định an toàn hồ đập vừa và nhỏ trên toàn quốc.
  4. Giảng viên, học viên cao học và sinh viên chuyên ngành Công trình thủy, Địa kỹ thuật công trình: Sử dụng tài liệu như một nghiên cứu điển hình chuyên sâu về sự kết hợp giữa cơ học đất, thủy lực dòng hở và phương pháp phần tử hữu hạn trong phân tích an toàn đập.

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

Tại sao đập đất truyền thống bị cấm tuyệt đối không được để nước tràn qua đỉnh? Đập đất được đắp bằng vật liệu hạt rời có khả năng kháng xói rất thấp. Khi nước tràn đỉnh tự do với lưu tốc trên 1,5 m/s, dòng chảy xiết nhanh chóng đào bới tạo các rãnh xói sâu ở mái hạ lưu, kích hoạt hiện tượng xói ngầm và xói lở lũy tiến dẫn tới vỡ đập hoàn toàn chỉ sau 1 đến 3 giờ.

Phương pháp Bishop đơn giản và Janbu khác nhau như thế nào trong phân tích ổn định? Phương pháp Bishop áp dụng cho mặt trượt trụ tròn dựa trên điều kiện cân bằng mô men và bỏ qua lực cắt giữa các thỏi đất, cho kết quả hệ số K chuẩn xác với mái dốc đồng chất. Ngược lại, phương pháp Janbu giải quyết mặt trượt bất kỳ dựa trên cân bằng lực, phù hợp cho mái đập có phân tầng địa chất phức tạp.

Khi lũ vượt thiết kế, nên đắp nâng cao đỉnh đập hay chủ động cho nước tràn qua đỉnh đã gia cố? Nếu đỉnh đập và mái hạ lưu được gia cố chống xói tốt bằng vải địa, bao tải cát hoặc bê tông, giải pháp cho tràn đỉnh cột nước thấp (dưới 0,5m) giúp xả bớt lượng lũ, giảm mực nước hồ dâng cao gây mất ổn định tổng thể. Nếu không có điều kiện gia cố bề mặt, đắp con chạch ngăn nước là giải pháp bắt buộc.

Bộ phần mềm Geo-Slope giải quyết bài toán thấm và ổn định liên hoàn như thế nào? Module SEEP/W sử dụng phương pháp phần tử hữu hạn hai chiều để giải phương trình dòng thấm, tính toán chi tiết áp lực nước lỗ rỗng và vị trí đường bão hòa. Toàn bộ dữ liệu áp lực này được truyền tự động sang module SLOPE/W để tính toán hệ số an toàn trượt K chính xác theo từng lát cắt phân thỏi.

Hiệu quả thực tế của giải pháp chống xói mái đập bằng thảm cỏ Vetiver ra sao? Thảm cỏ Vetiver sở hữu mạng lưới rễ đan xen dày đặc ăn sâu từ 1 đến 4 mét vào thân đập, đóng vai trò như các neo đất sinh học giúp cố định đất đắp, tăng lực dính kết và giảm thiểu hơn 80% nguy cơ xói mòn bề mặt do dòng tràn vận tốc dưới 2,0 m/s.

Kết luận

  • Luận văn đã làm sáng tỏ cơ chế tác động của dòng thấm và áp lực nước thủy động học lên độ ổn định mái đập đất khi xuất hiện cột nước tràn đỉnh.
  • Xây dựng thành công mô hình số phần tử hữu hạn liên hoàn SEEP/W và SLOPE/W trên 8 kịch bản hình học và 6 cấp cột nước vượt đỉnh từ 0,1m đến 1,1m.
  • Chứng minh tính khả thi về mặt kỹ thuật của giải pháp cho nước tràn qua đỉnh đập có gia cố, bảo đảm hệ số an toàn trượt K luôn lớn hơn ngưỡng quy chuẩn 1,15.
  • Đề xuất hệ thống giải pháp gia cố kết hợp linh hoạt giữa kết cấu cứng (bê tông FS, tấm lát) và kết cấu mềm (vải địa kỹ thuật, thảm cỏ Vetiver rễ sâu 1 - 4m).
  • Ứng dụng và kiểm chứng thành công cho công trình hồ chứa nước Chấn Sơn, tạo cơ sở nhân rộng cho hơn 4.300 hồ chứa nhỏ trong giai đoạn 2026-2030.

Các đơn vị tư vấn và cơ quan quản lý khai thác công trình thủy lợi cần chủ động đưa phương pháp mô phỏng số và giải pháp gia cố mái đập vào quy trình kiểm định an toàn định kỳ nhằm bảo vệ vững chắc các công trình hồ đập trước thiên tai cực đoan.