Tổng quan nghiên cứu

Thống kê lịch sử an toàn hồ đập thế giới ghi nhận 159 sự cố vỡ đập đất nghiêm trọng chỉ trong 30 năm đầu thế kỷ 20. Đáng chú ý, theo các báo cáo an toàn đập tại Mỹ, có tới 27 trong số 72 đập bị phá hoại hoàn toàn (chiếm tỷ lệ gần 40%) bắt nguồn từ hiện tượng dòng thấm tập trung gây xói ngầm. Tại Việt Nam, nhiều công trình thủy lợi trọng điểm từng đối mặt với nguy cơ thấm nghiêm trọng như đập Cà Giây xuất hiện thấm chân mái hạ lưu với lưu lượng từ 5 đến 7 lít/phút ngay trước khi nghiệm thu, hay các sự cố sình lầy, mạch sủi tại đập Am Chúa, hồ Vực Tròn, hồ Suối Hành và đập Kim Sơn.

Vấn đề cốt lõi đặt ra là trong thực tế địa hình, phần lớn các công trình đập vật liệu địa phương được xây dựng trên nền địa chất có đáy nằm nghiêng dốc về phía hạ lưu. Khi áp dụng các phương pháp thủy lực truyền thống, việc thiết lập bài toán thấm và ổn định mái dốc cho đáy nghiêng thường rất phức tạp, đòi hỏi thời gian mô phỏng kéo dài và thiếu các thông số cơ lý chi tiết trong giai đoạn tiền khả thi.

Mục tiêu nghiên cứu của luận văn là phân tích định lượng ảnh hưởng của góc nghiêng đáy đập đến lưu lượng thấm, vị trí đường bão hòa, gradient thủy lực thoát ra và hệ số ổn định mái dốc. Từ đó, đề tài thiết lập giải pháp chuyển đổi hình học từ sơ đồ đáy nghiêng sang sơ đồ đáy nằm ngang thông qua việc xác định chiều cao đập quy đổi tối ưu.

Nghiên cứu được thực hiện tại Trường Đại học Thủy Lợi trong giai đoạn 2010 đến 2013, thuộc chuyên ngành Xây dựng công trình thủy (Mã số: 60-58-40). Luận văn mang lại ý nghĩa thực tiễn lớn khi cung cấp bộ biểu đồ tra nhanh, giúp các kỹ sư rút ngắn hơn 60% thời gian tính toán sơ bộ mà vẫn đảm bảo độ chính xác kỹ thuật với sai số dưới 5%.

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 sự kết hợp chặt chẽ giữa thủy lực dòng thấm và cơ học đất hiện đại:

  • Lý thuyết thấm trong môi trường lỗ rỗng: Áp dụng định luật thực nghiệm Darcy xác lập năm 1856 về dòng thấm tuyến tính, kết hợp phương trình vi phân chuyển động của dòng nước dưới đất do Jucovxky đề xuất năm 1898 và phương trình vi phân chuyển động không ổn định của Boussinesq năm 1904. Bài toán thấm phẳng hai chiều được mô tả tổng quát qua phương trình Laplace và phương trình thế năng dòng thấm Muskat thiết lập năm 1937 đối với môi trường dị hướng, giải quyết mối quan hệ giữa hệ số thấm theo phương ngang và phương đứng.
  • Lý thuyết ổn định mái dốc: Dựa trên định luật sức bền cắt Mohr-Coulomb, xác định cường độ chống cắt của đất dựa trên ứng suất pháp hiệu quả, góc nội ma sát và lực dính đơn vị. Phương pháp cân bằng giới hạn phân thỏi Bishop rút gọn được áp dụng để tính toán hệ số an toàn ổn định mái dốc thông qua quá trình lặp phi tuyến. Thuật toán tìm kiếm mặt trượt nguy hiểm ngẫu nhiên Monte Carlo và phương pháp chuyển dịch điểm của Duncan và Celestino được tích hợp để xác định cung trượt tới hạn không tròn.
  • Các khái niệm then chốt: Lưu lượng thấm đơn vị qua thân đập ($10^{-6}\text{ m}^3/\text{s}$), gradient thủy lực thoát ra tại mái hạ lưu, đường bão hòa trong thân đập, hệ số an toàn ổn định trượt và chiều cao đập quy đổi tương đương.

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

Nghiên cứu ứng dụng phương pháp mô hình số phần tử hữu hạn (FEM) thông qua bộ công cụ GeoStudio (GeoSlope), tích hợp module SEEP/W tính thấm và SLOPE/W phân tích ổn định:

  • Cỡ mẫu và ma trận tổ hợp tính toán: Đề tài xây dựng một ma trận gồm 25 tổ hợp hình học đại diện. Chiều cao đập khảo sát gồm 5 cấp: 10 mét, 15 mét, 20 mét, 25 mét và 30 mét. Hệ số mái dốc khảo sát gồm 5 mức: 2.0, 2.5, 3.0, 3.5 và 4.0. Các góc nghiêng của nền dốc về phía hạ lưu khảo sát từ 0 độ đến 8 độ. Ba vị trí mặt cắt chiều cao quy đổi được phân tích gồm vị trí cách chân thượng lưu 1/3 chiều dài đáy đập, 1/2 chiều dài đáy đập và 2/3 chiều dài đáy đập.
  • Phương pháp chọn mẫu: Lựa chọn theo phương pháp quy hoạch thông số toàn phần, bao phủ toàn bộ dải kích thước hình học điển hình của các công trình đập đất quy mô vừa và nhỏ tại Việt Nam.
  • Lý do lựa chọn phương pháp: Phương pháp phần tử hữu hạn cho phép rời rạc hóa miền tính toán thành lưới các tam giác và tứ giác, giải hệ phương trình cân bằng dòng chảy tại các nút với độ chính xác cao. Việc kết hợp phương pháp số với phương pháp thủy lực giúp chuẩn hóa các kết quả phức tạp thành dạng biểu đồ tra cứu trực quan, phục vụ hiệu quả cho công tác tư vấn thiết kế.

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

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

Quá trình mô phỏng và phân tích số liệu trên 25 tổ hợp tính toán đã mang lại những kết quả định lượng rõ ràng:

  • Sự hạ thấp của đường bão hòa và gia tăng lưu lượng thấm: Khi góc nghiêng đáy đập dốc về hạ lưu tăng từ 0 độ lên 6 độ rồi 8 độ, đường bão hòa trong thân đập có xu hướng dốc xuống nhanh hơn về phía hạ lưu. Lưu lượng thấm đơn vị tăng từ 15% đến 28% tùy thuộc vào chiều cao đập. Ví dụ, với đập cao 20 mét có hệ số mái bằng 2.0, lưu lượng thấm tăng liên tục từ mức cơ sở khi đáy nằm ngang lên giá trị lớn hơn đáng kể khi góc nghiêng nền đạt 6 độ.
  • Gia tăng gradient thấm tại chân mái hạ lưu: Gradient thoát nước tại điểm tiếp giáp giữa mái hạ lưu và nền đất tăng mạnh theo góc dốc của đáy đập và chiều cao đập. Đối với đập cao 25 mét và 30 mét, giá trị gradient thấm tại chân đập tăng khoảng 18% đến 25%, tiệm cận ngưỡng gradient cho phép của đất đắp sét pha, làm tăng nguy cơ phát sinh mạch đùn, mạch sủi nếu không có tầng lọc tiêu nước.
  • Suy giảm hệ số an toàn ổn định mái dốc: Hệ số an toàn ổn định mái hạ lưu giảm từ 8% đến 15% khi góc dốc đáy tăng từ 0 độ đến 8 độ. Nguyên nhân do hướng dòng thấm nghiêng làm gia tăng áp lực thấm tác dụng lên các thỏi đất trượt ở phần hạ lưu đập.
  • Xác định vị trí chiều cao quy đổi tối ưu: Kết quả so sánh giữa mô hình đáy nghiêng thực tế và sơ đồ đáy ngang biến đổi cho thấy, vị trí chiều cao tại 1/2 chiều dài đáy đập mang lại kết quả lưu lượng và gradient thấm trùng khớp nhất, với sai số tính toán dưới 5%.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng thấm qua đập đáy nghiêng thể hiện ở chỗ: khi nền dốc về hạ lưu, quãng đường thấm bị rút ngắn cục bộ và thế năng dòng thấm dồn về chân mái hạ lưu mạnh hơn so với đập trên nền bằng phẳng. So sánh với các nghiên cứu thủy lực kinh điển của Dupuit hay Casagrande vốn chỉ áp dụng cho đáy nằm ngang, nghiên cứu này đã khỏa lấp khoảng trống khi lượng hóa được sự phân bố lại của áp lực nước lỗ rỗng dưới tác động của hình học đáy.

Toàn bộ hệ thống dữ liệu nghiên cứu được biểu diễn trực quan qua 20 bảng số liệu chi tiết và hơn 35 biểu đồ quan hệ hàm số. Các biểu đồ này mô tả chính xác tương quan giữa lưu lượng thấm với góc nghiêng, gradient thấm với góc nghiêng, hệ số an toàn với góc nghiêng, và tỷ số chiều cao quy đổi trên chiều cao thực tế theo từng hệ số mái dốc từ 2.0 đến 4.0. Việc trình bày dưới dạng họ đường cong giúp các kỹ sư dễ dàng nội suy thông số mà không cần lặp lại các bước chia lưới phần tử hữu hạn phức tạp.

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

Dựa trên các kết luận khoa học, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể nhằm nâng cao hiệu quả thiết kế và an toàn công trình:

  • Ứng dụng mặt cắt quy đổi tại vị trí giữa đáy đập trong thiết kế sơ bộ: Các đơn vị tư vấn thiết kế nên sử dụng giá trị chiều cao đập quy đổi tại vị trí 1/2 chiều dài đáy đập để chuyển đổi nhanh bài toán đập đáy nghiêng về bài toán đập đáy ngang. Giải pháp này giúp rút ngắn từ 60% đến 70% thời gian lập dự án tiền khả thi nhưng vẫn bảo đảm độ tin cậy cơ học cao.
  • Tăng cường kết cấu tiêu nước chân mái hạ lưu cho đập có đáy dốc: Đối với các công trình có góc nghiêng đáy đập lớn hơn hoặc bằng 4 độ, bắt buộc phải bố trí lăng trụ tiêu nước hoặc tầng lọc ngược bằng vật liệu rời có bề dày từ 1.5 đến 2.5 mét. Mục tiêu là khống chế gradient dòng thoát luôn nhỏ hơn giá trị giới hạn 0.5 đến 0.7, triệt tiêu nguy cơ xói ngầm.
  • Áp dụng công nghệ chống thấm hiện đại cho các đập hiện hữu: Đối với các đập đất đang vận hành trên nền dốc bị thấm mạnh hoặc sủi đất ở hạ lưu, các chủ quản lý cần triển khai ngay giải pháp khoan phụt vữa xi măng cao áp Jet-Grouting tạo màng chống thấm hoặc phủ màng địa kỹ thuật Geomembrane ở mái thượng lưu. Lộ trình xử lý cần hoàn thành trong vòng 12 đến 18 tháng sau khi phát hiện sự cố.
  • Chuẩn hóa bộ biểu đồ tra cứu vào tiêu chuẩn ngành: Kiến nghị cơ quan chuyên trách thuộc Bộ Nông nghiệp và Phát triển nông thôn bổ sung họ biểu đồ quan hệ tỷ số chiều cao quy đổi vào cẩm nang tính toán thiết kế thủy lợi, hướng tới mục tiêu đồng bộ hóa quy trình thẩm định hồ sơ đập đất trên nền nghiêng trong giai đoạn 24 tháng tới.

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

Nội dung và sản phẩm nghiên cứu của luận văn mang giá trị ứng dụng cao cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế công trình thủy lợi và thủy điện: Tra cứu trực tiếp các biểu đồ tương quan để xác định nhanh chiều cao đập thiết kế, mái dốc, lưu lượng thấm và gradient tới hạn trong giai đoạn khảo sát, lập báo cáo kinh tế - kỹ thuật.
  • Đơn vị quản lý, khai thác và bảo trì hồ đập: Nắm vững cơ chế phát sinh dòng thấm trên nền nghiêng để xây dựng phương án quan trắc, kịp thời phát hiện các dấu hiệu mất ổn định mái dốc và bão hòa chân đập trong mùa mưa lũ.
  • Giảng viên, nghiên cứu sinh và học viên cao học: Sử dụng tài liệu như một nghiên cứu mẫu mực về việc ứng dụng phương pháp phần tử hữu hạn giải quyết bài toán cơ học đất phức tạp và phương pháp phân tích ổn định mái dốc.
  • Cơ quan quản lý nhà nước và ban quản lý dự án: Có thêm cơ sở khoa học định lượng để thẩm định, đánh giá an toàn đập và phê duyệt phương án gia cố các hồ chứa nước xung yếu tại các địa phương.

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

  • Tại sao đáy đập nghiêng về hạ lưu lại làm tăng lưu lượng thấm và giảm ổn định mái dốc? Khi đáy đập dốc về phía hạ lưu, độ dốc thủy lực cục bộ tăng lên, làm dòng thấm thoát nhanh hơn và kéo lưu lượng thấm tăng thêm từ 15% đến 28%. Đồng thời, áp lực nước lỗ rỗng tập trung mạnh tại vùng chân đập hạ lưu, làm suy giảm sức kháng cắt của đất và khiến hệ số ổn định giảm từ 8% đến 15%.
  • Vị trí nào được khuyến nghị để xác định chiều cao đập quy đổi về đáy nằm ngang? Nghiên cứu đã khảo sát 3 vị trí gồm 1/3, 1/2 và 2/3 chiều dài đáy đập. Kết quả phân tích số khẳng định mặt cắt tại vị trí 1/2 chiều dài đáy đập mang lại sự tương thích tốt nhất về lưu lượng và gradient thấm so với đập đáy nghiêng thực tế, với sai số dưới 5%.
  • Phương pháp Bishop đơn giản có ưu thế gì trong bài toán ổn định mái dốc này? Phương pháp Bishop đơn giản xét đến điều kiện cân bằng mô-men của toàn bộ khối trượt và lực tương tác ngang giữa các thỏi đất. Phương pháp này giải quyết được tính siêu tĩnh của bài toán phân thỏi, cho kết quả hệ số an toàn chính xác, sát với thực tế sạt trượt và rất thuận lợi khi lập trình tính toán tự động.
  • Khi nào kỹ sư nên dùng biểu đồ tra cứu thay vì chạy phần mềm phần tử hữu hạn? Trong giai đoạn thiết kế sơ bộ hoặc tiền khả thi, khi chưa có đầy đủ số liệu khảo sát chi tiết về địa chất và hệ số thấm, kỹ sư nên sử dụng biểu đồ tra cứu trong luận văn để chọn nhanh kích thước đập tối ưu, giúp tiết kiệm hơn 60% thời gian thực hiện.
  • Giải pháp xử lý thấm khẩn cấp nào hiệu quả nhất cho các đập đất hiện hữu có đáy dốc? Biện pháp xử lý triệt để là thi công tầng phủ chống thấm mái thượng lưu kết hợp khoan phụt màng chống thấm Jet-Grouting vào nền. Đồng thời, tại mái hạ lưu cần bổ sung ngay lăng trụ phản áp và tầng lọc ngược dày từ 1.5 đến 2.5 mét để hạ thấp đường bão hòa.

Kết luận

  • Lượng hóa chính xác tác động tiêu cực của góc nghiêng đáy đập dốc về hạ lưu làm tăng lưu lượng thấm lên tới 28% và làm giảm hệ số an toàn ổn định mái dốc tới 15%.
  • Thiết lập thành công ma trận 25 tổ hợp mô phỏng số, xây dựng hệ thống 20 bảng số liệu và hơn 35 đồ thị tương quan chuẩn hóa cho các cấp chiều cao đập từ 10 mét đến 30 mét.
  • Khẳng định vị trí 1/2 chiều dài đáy đập là mặt cắt biến đổi hình học tối ưu nhất, cho phép chuyển đổi bài toán đáy nghiêng sang đáy ngang với sai số kỹ thuật dưới 5%.
  • Đóng góp giải pháp khoa học kết hợp giữa phương pháp mô phỏng số phần tử hữu hạn hiện đại và phương pháp tra cứu thủy lực đơn giản trong kỹ thuật thủy công.
  • Cung cấp công cụ tra cứu thực chứng giá trị cao cho các kỹ sư công trình thủy trong tính toán nhanh thấm và ổn định đập vật liệu địa phương.

Về định hướng tiếp theo, nghiên cứu cần được mở rộng cho các dạng kết cấu đập đất không đồng chất, đập có tường lõi sét và xét đến tác động của tải trọng động đất trong lộ trình từ 24 đến 36 tháng tới. Các đơn vị tư vấn và ban quản lý hồ đập hãy ứng dụng ngay các biểu đồ tra cứu của luận văn để nâng cao năng suất thiết kế và bảo đảm an toàn bền vững cho các công trình hồ chứa thủy lợi.