Tổng quan nghiên cứu

Hồ chứa nước là loại hình công trình thủy lợi phổ biến nhất tại Việt Nam, với hơn 550 hồ chứa loại vừa và lớn (dung tích từ 1 triệu m³, chiều cao đập trên 10m) tính đến năm 2000, cùng hàng ngàn hồ nhỏ khác. Trong số đó, đập vật liệu địa phương chiếm tới 84% tổng số đập có chiều cao dưới 100m, vượt xa đập bê tông trọng lực (12%) và đập bê tông vòm (4%). Tuy nhiên, thống kê của Middle Brooks cho thấy trên 60% sự cố công trình đất trên thế giới do thấm gây ra, cộng thêm khoảng 10% sự cố có tác nhân kích thích từ thấm. Riêng nguyên nhân liên quan đến xử lý nền chiếm trên 40% các sự cố đập trong và ngoài nước.

Luận văn thạc sĩ "Nghiên cứu giải pháp xử lý nền đập đất đầm nén hồ chứa nước Mỹ Lâm - Phú Yên" của tác giả Bùi Thị Lương, thực hiện dưới sự hướng dẫn của PGS.TS Trịnh Minh Thụ tại Trường Đại học Thủy lợi và hoàn thành tháng 5/2013, giải quyết trực tiếp bài toán này. Mục tiêu nghiên cứu gồm bốn nội dung cụ thể: tổng quan các giải pháp xử lý chống thấm cho nền; làm rõ cơ sở khoa học và thực tiễn của từng phương pháp; tính toán ổn định thấm cho mỗi phương án xử lý; và kiểm tra ổn định trượt, ổn định biến dạng lún cùng ổn định trượt mái khi nước rút nhanh với giải pháp được lựa chọn.

Phạm vi nghiên cứu tập trung vào đập đất đầm nén thuộc dự án hồ chứa nước Mỹ Lâm, tỉnh Phú Yên, với mực nước dâng bình thường +33,4m và mực nước chết +15m. Ý nghĩa thực tiễn của đề tài nằm ở chỗ đưa ra được cơ sở định lượng để so sánh ba nhóm giải pháp xử lý nền, giúp giảm lưu lượng thấm mất nước hồ chứa, hạ thấp đường bão hòa và khống chế gradient thấm trong giới hạn cho phép, qua đó nâng cao hệ số an toàn ổn định mái đập.

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

Khung lý thuyết áp dụng

Nền tảng lý thuyết thứ nhất là định luật thấm đường thẳng Darcy (1856): lưu lượng thấm Q qua mặt cắt vuông góc dòng chảy tỷ lệ thuận với diện tích mặt cắt F và độ chênh cột nước H, tỷ lệ nghịch với chiều dài đường thấm L, biểu diễn qua vận tốc thấm v = k·i. Đây là định luật áp dụng cho chế độ chảy tầng trong môi trường lỗ rỗng nhỏ.

Nền tảng thứ hai là lý thuyết thấm phi tuyến, áp dụng khi môi trường có hạt lớn (đường kính hạt từ 0,57cm đến 5,63cm) và dòng thấm chuyển sang chế độ chảy rối. Luận văn dẫn các công thức của G. Cơrôbe, P. và Pavlovxki (v = ε·A·√J), trong đó hệ số A xác định theo công thức kinh nghiệm của X. Izbas cho đá đổ.

Bên cạnh đó là lý thuyết biến hình thấm — khung phân loại các dạng phá hoại do dòng thấm gây ra. Năm khái niệm chính được sử dụng xuyên suốt gồm: xói ngầm cơ học (hạt nhỏ bị đẩy lọt qua khe hở giữa hạt lớn, làm tăng độ rỗng và lưu lượng thấm); xói tiếp xúc (khống chế bằng hệ số chuyển tiếp cỡ hạt ξ ≤ 3÷10); đẩy trồi đất, xảy ra khi gradient thấm cửa ra vượt gradient giới hạn J_gh; đùn đất tiếp xúc; và độ bền thấm chung J_K < J_KCP theo phương pháp VNIIG (Liên Xô cũ).

Mô hình nghiên cứu tổng thể là mô hình phần tử hữu hạn hai chiều, trong đó miền thấm được chia thành các phần tử tam giác — vùng cường độ thấm mạnh chia lưới mịn hơn. Nghiệm bài toán được tìm bằng cực tiểu hóa phiếm hàm năng lượng E(h) = ΣE_e trên toàn miền, ẩn số là cột nước h tại các nút lưới.

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

Nguồn dữ liệu gồm hai nhóm. Nhóm thứ nhất là tài liệu khảo sát địa hình, địa chất công trình hồ Mỹ Lâm do Công ty tư vấn và chuyển giao công nghệ - Trường Đại học Thủy lợi cung cấp, được hệ thống thành các chỉ tiêu cơ lý tính toán của đất nền và vật liệu đất đắp. Nhóm thứ hai là dữ liệu thứ cấp về sự cố đập, gồm thống kê 48 đập đất, đá lớn hoàn thành từ 1963 đến 2007 và các hồ sơ sự cố điển hình (Suối Hành, Suối Trầu, Am Chúa, Cà Giây, Phú Ninh, Vực Tròn).

Về cỡ mẫu tính toán: nghiên cứu xây dựng bốn trường hợp mô phỏng trên cùng một mặt cắt đại diện của đập chính — trường hợp 1 là nền thiên nhiên chưa xử lý, ba trường hợp còn lại tương ứng ba giải pháp xử lý. Cách chọn mẫu là chọn mẫu có chủ đích: mặt cắt được lựa chọn tại vị trí bất lợi nhất về chiều cao đập và điều kiện địa chất nền, vì đây là mặt cắt khống chế an toàn cho toàn tuyến.

Công cụ phân tích là bộ phần mềm địa kỹ thuật GEO-SLOPE với ba mô-đun: SEEP/W giải bài toán thấm (đường bão hòa, lưu lượng, đường đẳng gradient, lưới thấm), SIGMA/W phân tích ứng suất – biến dạng để xác định độ lún, và SLOPE/W tính hệ số ổn định mái theo phương pháp mặt trượt. Lý do lựa chọn phương pháp phần tử hữu hạn thay cho phương pháp cơ học chất lỏng hay thủy lực học cổ điển: hai phương pháp cổ điển chỉ giải được sơ đồ đơn giản hoặc cho kết quả trung bình, trong khi bài toán Mỹ Lâm có biên phức tạp, nhiều lớp vật liệu với hệ số thấm khác nhau và cần xác định gradient cục bộ tại cửa ra dòng thấm.

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

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

Thứ nhất, nền thiên nhiên chưa xử lý không đáp ứng yêu cầu ổn định thấm. Kết quả tính trường hợp 1 cho thấy đường bão hòa nằm cao hơn dự kiến thiết kế, lưu lượng thấm qua nền vượt ngưỡng cho phép, và vùng gradient tập trung xuất hiện tại khu vực dòng thấm thoát ra hạ lưu — đúng vùng nguy hiểm về xói ngầm cơ học mà lý thuyết cảnh báo. Đây là kết quả nhất quán với quan trắc thực tế: phần lớn lưu lượng thấm đo được tại các đập miền Trung thường lớn hơn nhiều lần so với lượng thấm tính toán thiết kế.

Thứ hai, cả ba giải pháp xử lý đều làm giảm lưu lượng thấm nhưng với mức độ chênh lệch rõ rệt. So sánh bốn trường hợp cho thấy tường hào xi măng – bentonite hạ thấp đường bão hòa mạnh nhất và cắt gần như triệt để dòng thấm qua tầng nền thấm nước, do tường cắm xuyên qua lớp thấm và tựa vào tầng cách nước. Tường nghiêng kết hợp sân phủ chỉ kéo dài đường viền thấm, làm giảm gradient trung bình chứ không cắt dòng, nên hiệu quả giảm lưu lượng thấp hơn. Khoan phụt tạo màng chống thấm cho hiệu quả trung gian, phụ thuộc mạnh vào khả năng tiếp nhận vữa của nền.

Thứ ba, tường hào xi măng – bentonite được chọn trên cả hai tiêu chí kinh tế và kỹ thuật. Dữ liệu so sánh này có thể trình bày hiệu quả qua một bảng đối chiếu bốn cột (lưu lượng thấm, gradient cửa ra, cao độ đường bão hòa, khối lượng – chi phí), kèm biểu đồ cột nhóm cho ba chỉ tiêu kỹ thuật để thấy rõ khoảng cách giữa các phương án.

Thứ tư, giải pháp được chọn thỏa mãn kiểm tra ổn định trượt và biến dạng. Với mực nước thượng lưu ở MNDBT +33,4m và mực nước hạ lưu +12m, hệ số ổn định mái hạ lưu tính bằng SLOPE/W đạt trị số lớn hơn hệ số an toàn cho phép; giá trị lún lớn nhất tính bằng SIGMA/W nằm trong giới hạn cho phép của đập đất đầm nén. Trường hợp mực nước rút nhanh từ MNDBT +33,4m xuống mực nước chết +15m — kịch bản bất lợi nhất cho mái thượng lưu — hệ số ổn định vẫn đạt yêu cầu.

Thảo luận kết quả

Nguyên nhân khiến tường hào bentonite vượt trội nằm ở cơ chế làm việc: đây là giải pháp cắt dòng theo phương đứng, thay đổi trực tiếp sơ đồ thấm chứ không chỉ kéo dài đường viền. Trong khi đó, tường nghiêng – sân phủ làm việc theo nguyên lý tăng tổng hệ số sức cản Σξᵢ của miền thấm, hiệu quả giảm dần khi tầng thấm nền dày.

Kết quả này phù hợp với xu hướng đã ghi nhận trong thực tiễn ngành: luận văn dẫn danh mục các công trình đã áp dụng thành công tường hào bentonite, khoan phụt (điển hình là Tân Giang – Ninh Thuận) và cọc đất xi măng. Điểm tương đồng với các nghiên cứu trước là khẳng định không có giải pháp vạn năng — mỗi phương pháp chỉ phù hợp với một kiểu cấu trúc đất nền nhất định.

Về biến dạng, kết quả tính lún của nghiên cứu củng cố nhận định lý thuyết rằng độ lún đỉnh đập tăng khi chiều cao đập và chiều dày nền tăng, nhưng không tỷ lệ tuyến tính — tỷ lệ giữa độ lún và chiều cao đập tăng nhanh hơn khi đập cao lên. Đây là cảnh báo quan trọng: tính khối lượng đắp bù lún theo quan hệ tuyến tính sẽ thiếu hụt vật liệu.

Ý nghĩa lớn nhất của kết quả là cung cấp quy trình kiểm chứng khép kín — thấm, trượt, lún, rút nước nhanh — thay vì chỉ dừng ở bài toán thấm như nhiều thiết kế trước đây. Bài học từ Suối Hành, Suối Trầu và Am Chúa cho thấy sự cố thường bắt nguồn từ việc bỏ sót một khâu kiểm tra.

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

Thứ nhất, áp dụng tường hào xi măng – bentonite làm giải pháp chống thấm nền chủ đạo cho đập Mỹ Lâm. Chủ thể thực hiện là đơn vị tư vấn thiết kế phối hợp ban quản lý dự án, với mục tiêu đưa lưu lượng thấm qua nền về dưới ngưỡng cho phép và giữ gradient cửa ra thấp hơn gradient giới hạn. Thi công theo trình tự panel sơ cấp – panel thứ cấp, kèm kiểm tra chất lượng vữa tại trạm ủ trộn trước mỗi ca.

Thứ hai, bổ sung tầng lọc ngược tại toàn bộ khu vực dòng thấm thoát ra hạ lưu. Tầng lọc cấu tạo nhiều lớp vật liệu không dính với đường kính hạt tăng dần theo chiều dòng thấm, hệ số chuyển tiếp giữa hai lớp kề nhau khống chế trong khoảng ξ = 3÷10. Đơn vị thiết kế cần đưa yêu cầu này thành chỉ dẫn kỹ thuật bắt buộc, hoàn thành trong giai đoạn thiết kế bản vẽ thi công.

Thứ ba, siết chặt công tác khảo sát vật liệu đất đắp với ba chỉ tiêu bắt buộc: độ tan rã, độ lún ướt và độ trương nở. Chính ba chỉ tiêu này đã bị bỏ sót trong khảo sát đập Suối Hành, dẫn tới đánh giá sai tính hoàng thổ của bãi vật liệu và góp phần gây vỡ đập năm 1986. Đơn vị khảo sát địa chất phải thí nghiệm đủ, đồng thời số lượng mẫu kiểm tra dung trọng phải đạt 100% quy định của tiêu chuẩn thay vì mức khoảng 10% từng ghi nhận tại Suối Trầu.

Thứ tư, thiết lập quy trình vận hành khống chế tốc độ hạ mực nước hồ. Đơn vị quản lý khai thác cần ban hành quy trình giới hạn tốc độ rút nước từ MNDBT về MNC, tránh lặp lại kịch bản đã gây 200 vị trí sạt trượt trong phạm vi 17km lòng hồ khi đập Teton vỡ ngày 5/6/1976.

Thứ năm, lắp đặt hệ thống quan trắc thấm và lún dài hạn. Bố trí ống đo áp theo mặt cắt đại diện và mốc quan trắc lún tại đỉnh, mái đập; đọc số liệu định kỳ trong ba năm đầu tích nước — giai đoạn nhiều đập phát sinh thấm mạnh, có công trình chỉ sau 1-2 năm khai thác.

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

Kỹ sư thiết kế công trình thủy lợi và địa kỹ thuật là nhóm hưởng lợi trực tiếp nhất. Luận văn cung cấp bộ thông số đầu vào đã kiểm chứng cho bốn trường hợp mô hình GEO-SLOPE, từ chỉ tiêu cơ lý đất nền đến điều kiện biên bài toán. Use case cụ thể: khi thiết kế đập đất trên nền thấm nước ở khu vực Nam Trung Bộ, có thể lấy sơ đồ tính và cách khai báo điều kiện biên của luận văn làm mẫu đối chiếu, rút ngắn đáng kể thời gian thiết lập mô hình.

Học viên cao học và nghiên cứu sinh ngành Xây dựng công trình thủy tìm thấy ở đây một mẫu luận văn có cấu trúc bốn chương chặt chẽ, đi từ tổng quan sự cố, cơ sở lý thuyết, phân tích giải pháp đến ứng dụng tính toán cho công trình cụ thể. Phần cơ sở lý thuyết về biến hình thấm — với đầy đủ công thức gradient giới hạn, chiều dày tầng gia trọng và độ bền thấm chung — là tài liệu tham khảo trực tiếp cho các đề tài tương tự.

Cán bộ quản lý dự án và đơn vị khai thác hồ chứa khai thác giá trị ở Chương 1: bảng thống kê 48 đập lớn và các bài học sự cố được phân tích tới tận nguyên nhân gốc. Use case: xây dựng danh mục kiểm tra rủi ro cho các hồ chứa đang vận hành, ưu tiên rà soát nhóm hồ vừa và nhỏ — nơi tập trung đa số sự cố dù dung tích không lớn, như trường hợp hồ 250.000m³ ở Hà Tĩnh vỡ làm trôi hơn 200m đường sắt Bắc Nam.

Nhà thầu thi công xử lý nền tham khảo phần mô tả công nghệ tường hào bentonite, khoan phụt và cọc đất xi măng, kèm hình ảnh thi công thực tế và danh mục công trình đã áp dụng. Điều này giúp lập biện pháp thi công và dự toán sát thực tế hơn khi đấu thầu các gói xử lý nền.

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

1. Vì sao thấm là nguyên nhân sự cố hàng đầu của đập đất? Vì đập đất và nền đập về bản chất đều thấm nước. Khi mực nước thượng lưu dâng, dòng thấm từ thượng về hạ lưu hình thành và tạo lực thủy động đẩy nổi làm mái dễ mất ổn định. Thống kê Middle Brooks cho thấy trên 60% sự cố công trình đất do thấm trực tiếp gây ra, cộng thêm 10% có tác nhân kích thích từ thấm — tổng cộng khoảng 70%.

2. Tường hào xi măng – bentonite khác gì tường nghiêng kết hợp sân phủ? Tường hào là giải pháp cắt dòng theo phương đứng, cắm xuyên tầng thấm để chặn dòng thấm nền. Tường nghiêng – sân phủ hoạt động theo nguyên lý kéo dài đường viền thấm, làm giảm gradient trung bình mà không cắt dòng. Kết quả tính toán tại Mỹ Lâm cho thấy tường hào hạ đường bão hòa mạnh hơn rõ rệt, đặc biệt khi tầng thấm nền dày.

3. Tại sao phải kiểm tra riêng trường hợp mực nước rút nhanh? Khi mực nước trên mái rút nhanh và hệ số thấm của đất nhỏ, áp lực kẽ rỗng trong khối đất hầu như không kịp tiêu tán, trong khi tác dụng phản áp giữ ổn định của khối nước bên ngoài mất đi. Chênh lệch này làm mái thượng lưu mất ổn định. Luận văn kiểm tra kịch bản rút từ +33,4m xuống +15m và xác nhận hệ số ổn định vẫn đạt yêu cầu.

4. Xói ngầm cơ học được phòng chống bằng cách nào? Vùng nguy hiểm nhất là nơi dòng thấm thoát ra hạ lưu. Biện pháp cơ bản là bố trí thiết bị thoát nước dạng tầng lọc ngược, gồm nhiều lớp cát, sỏi, cuội, đá dăm có đường kính hạt tăng dần theo chiều dòng thấm. Hệ số chuyển tiếp cỡ hạt giữa hai lớp kề nhau nên giữ trong khoảng 3 đến 10 để chống cả xói tiếp xúc.

5. Vì sao dùng phần tử hữu hạn thay cho phương pháp tính thấm cổ điển? Phương pháp cơ học chất lỏng chỉ giải được sơ đồ đơn giản và bế tắc khi điều kiện biên phức tạp. Phương pháp thủy lực chỉ cho đặc trưng trung bình của dòng thấm. Phần tử hữu hạn giải được miền có hình dạng và tính chất vật liệu phức tạp, đồng thời cho đầy đủ đường bão hòa, lưu tốc, gradient, cột nước và lưu lượng thấm tại mọi điểm.

Kết luận

  • Đập vật liệu địa phương chiếm 84% số đập cao dưới 100m tại Việt Nam, trong khi sự cố do nguyên nhân xử lý nền chiếm trên 40% — khoảng trống này chính là lý do tồn tại của nghiên cứu.
  • Luận văn xây dựng và giải bốn trường hợp tính toán trên GEO-SLOPE, chứng minh nền thiên nhiên hồ Mỹ Lâm chưa đáp ứng yêu cầu ổn định thấm và định lượng được hiệu quả của từng giải pháp xử lý.
  • Tường hào xi măng – bentonite được kiến nghị là giải pháp tối ưu cho đập Mỹ Lâm trên cả hai tiêu chí kinh tế và kỹ thuật.
  • Giải pháp chọn đã vượt qua trọn bộ ba kiểm tra: ổn định trượt mái hạ lưu tại MNDBT +33,4m, biến dạng lún, và ổn định trượt mái khi nước rút nhanh về MNC +15m.
  • Đóng góp chính là quy trình phân tích khép kín thấm – trượt – lún, có thể chuyển giao cho các đập đất đầm nén tương tự ở Nam Trung Bộ.

Hướng nghiên cứu tiếp theo nên mở rộng sang bài toán thấm không ổn định theo thời gian và đất không bão hòa, đồng thời đối chiếu kết quả mô hình với số liệu quan trắc thực tế sau tích nước. Kỹ sư và đơn vị quản lý đang xử lý nền đập đất trên nền thấm nước nên tham khảo trực tiếp bộ thông số và sơ đồ tính của luận văn để rút ngắn quá trình thiết lập mô hình.