Tổng quan nghiên cứu

Việt Nam hiện có khoảng 6.755 hồ chứa thủy lợi với tổng dung tích trữ nước đạt xấp xỉ 11 tỷ m³, trong đó có hơn 67% số hồ được xây dựng từ sau năm 1976. Phần lớn các công trình đầu mối là đập đất đá được thi công trong điều kiện công nghệ trước đây còn hạn chế. Thống kê của các tổ chức an toàn đập quốc tế chỉ ra rằng nguyên nhân sự cố liên quan đến dòng thấm chiếm khoảng 20%, trong khi tại Việt Nam tỷ lệ này chiếm khoảng 15,65% tổng số sự cố đập. Đặc biệt, theo nghiên cứu của Foster và cộng sự, có hơn 50% các sự cố vỡ đập do thấm xuất hiện ngay trong năm đầu tiên khi hồ bắt đầu tích nước. Thực tế thiết kế hiện nay thường chỉ chú trọng bài toán thấm ổn định khi mực nước đã đạt cao trình thiết kế mà chưa đánh giá đầy đủ giai đoạn tích nước ban đầu.

Đề tài nghiên cứu tập trung phân tích diễn biến dòng thấm và đánh giá mức độ ổn định của thân đập trong suốt quá trình tích nước tại công trình hồ chứa nước Ngàn Trươi, tỉnh Hà Tĩnh. Đây là công trình đầu mối thủy lợi cấp II trọng điểm với đập chính bằng đất đá nhiều khối có chiều cao lớn nhất 53,5 m, chiều dài đỉnh đập 345 m và dung tích toàn bộ đạt 775,70 triệu m³.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi tích hợp mô hình số phần tử hữu hạn với dữ liệu quan trắc thực địa tại 4 tuyến đo áp chuyên dụng. Kết quả tính toán giúp kiểm soát áp lực nước kẽ rỗng thặng dư, nâng cao độ chính xác dự báo an toàn đập lên trên 90%, đồng thời thiết lập cơ sở khoa học để tối ưu hóa quy trình tích nước an toàn cho hệ thống hồ chứa quy mô lớn.

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 nền tảng thủy lực ngầm và cơ học đất hiện đại để giải quyết bài toán thấm qua môi trường rỗng không đồng tính. Khung lý thuyết bao gồm hai mô hình toán học cốt lõi:

  1. Lý thuyết thấm mở rộng theo định luật Darcy cho môi trường bão hòa và không bão hòa hai chiều, giải phương trình vi phân đạo hàm riêng mô tả trường áp lực cột nước thấm thay đổi theo thời gian thực.
  2. Phương pháp cân bằng giới hạn theo lý thuyết Bishop mở rộng (phương pháp ứng suất hiệu quả) và phương pháp tổng ứng suất, phục vụ đánh giá hệ số ổn định trượt mái dốc thượng và hạ lưu đập.

Các khái niệm chuyên ngành nền tảng được chuẩn hóa gồm:

  • Áp lực nước kẽ rỗng: Đại lượng áp lực nước trong các lỗ rỗng của cốt đất, quyết định trực tiếp tới cường độ chống cắt hiệu quả của khối đất đắp.
  • Đường bão hòa: Biên trên của miền dòng thấm không áp trong thân đập, nơi áp suất chất lỏng bằng áp suất khí quyển.
  • Gradient thủy lực: Tỷ số tiêu hao cột nước trên một đơn vị chiều dài dòng thấm, chỉ số quyết định nguy cơ xói ngầm và biến hình thấm.
  • Hệ số trữ nước riêng: Đại lượng phản ánh thể tích nước được giải phóng hoặc hấp thụ bởi một đơn vị thể tích đất khi cột áp thủy lực thay đổi một đơn vị.
  • Hệ số an toàn chống trượt: Tỷ số giữa tổng mômen chống trượt và tổng mômen gây trượt dọc theo mặt trượt trụ tròn giả định.

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

Nguồn dữ liệu phục vụ nghiên cứu được trích xuất trực tiếp từ hệ thống 4 tuyến quan trắc thấm tại đập chính hồ Ngàn Trươi, gồm các mặt cắt ngang 3-3, 6-6, 9-9 và 12-12. Hệ thống thu thập số liệu tự động theo tiêu chuẩn quốc tế với hàng chục cảm biến đo áp lực nước kẽ rỗng dạng màng rung và ống đo mực nước bão hòa. Cơ sở dữ liệu bao gồm chuỗi quan trắc liên tục theo các chu kỳ dâng mực nước hồ từ mực nước chết cao độ +25,00 m lên mực nước dâng bình thường +52,00 m và mực nước lớn nhất kiểm tra +55,16 m.

Nghiên cứu sử dụng phương pháp phần tử hữu hạn thông qua module SEEP/W và SLOPE/W thuộc bộ phần mềm chuyên dụng GeoStudio. Cỡ mẫu nghiên cứu bao gồm 5 kịch bản tích nước điển hình với các tốc độ dâng nước thay đổi từ 0,2 m/ngày đến 1,0 m/ngày. Phương pháp chọn mẫu tập trung phân tích chuyên sâu tại mặt cắt ngang lớn nhất (mặt cắt 6-6 với chiều cao đập 53,5 m) do đây là vị trí xung yếu chịu gradient thấm và áp lực kẽ rỗng cao nhất.

Lý do lựa chọn phương pháp mô hình số phần tử hữu hạn là do khả năng mô phỏng chính xác miền thấm phức tạp gồm nhiều khối vật liệu đất đắp, lăng trụ tiêu nước, lõi chống thấm và màn chắn khoan phụt xi măng sâu trong nền đá. Thời gian thực hiện tính toán mô phỏng và hiệu chỉnh thông số mô hình kéo dài trong chu kỳ tích nước thực tế của công trình từ giai đoạn bắt đầu dâng nước đến khi đạt trạng thái vận hành ổn định năm 2022.

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 số kết hợp hiệu chỉnh với dữ liệu quan trắc thực địa đã đem lại những phát hiện quan trọng:

  1. Tồn tại hiện tượng trễ pha rõ rệt của áp lực nước kẽ rỗng so với tốc độ dâng của mực nước thượng lưu. Tại vùng lõi chống thấm, áp lực nước kẽ rỗng biến thiên chậm hơn khoảng 18% đến 22% so với tốc độ biến thiên mực nước ngoài hồ, chứng minh khả năng cản thấm vượt trội của vật liệu đất sét đầm nén.
  2. Đường bão hòa thực tế thu được từ chuỗi cảm biến quan trắc tại các giếng đo áp thấp hơn khoảng 10% đến 15% so với kết quả tính toán thuần lý thuyết ở điều kiện thấm ổn định truyền thống. Sự chênh lệch này giúp giảm áp lực đẩy ngược lên khối lăng trụ hạ lưu khoảng 12,5 kPa.
  3. Hệ số an toàn ổn định mái thượng lưu giảm mạnh nhất trong kịch bản tích nước nhanh rồi rút nước đột ngột (từ cao trình +52,00 m xuống +25,00 m), với hệ số an toàn giảm từ 1,85 xuống còn 1,38. Tuy nhiên, giá trị này vẫn cao hơn ngưỡng an toàn cho phép theo tiêu chuẩn quốc gia đối với công trình cấp II là 1,30.
  4. Tốc độ dâng mực nước tối ưu được xác lập là dưới 0,5 m/ngày. Ở tốc độ này, hơn 85% áp lực nước kẽ rỗng thặng dư được tiêu tán kịp thời qua hệ thống thoát nước hạ lưu, giữ gradient thấm tại chân khay chống thấm luôn nhỏ hơn 0,20, triệt tiêu nguy cơ xói ngầm.

Thảo luận kết quả

Sự sai khác giữa kết quả mô hình lý thuyết và dữ liệu quan trắc hiện trường bắt nguồn từ đặc tính bão hòa dần của đất đắp trong giai đoạn đầu tích nước. Khi nước hồ dâng, không khí bị giữ lại trong các lỗ rỗng của đất làm giảm hệ số thấm thực tế xuống khoảng 1,5 đến 2 lần so với trạng thái bão hòa hoàn toàn trong phòng thí nghiệm. Dữ liệu này được trực quan hóa rõ nét qua các biểu đồ phân bố áp lực cột nước tại các giếng quan trắc ALT5, ALT6, ALT7 và bảng đối chiếu hệ số an toàn trượt theo từng pha dâng nước.

So sánh với các nghiên cứu của Reinus về tiêu tán áp lực nước kẽ rỗng theo chỉ số vận tốc rút nước, kết quả tại hồ Ngàn Trươi hoàn toàn tương thích với quy luật dòng thấm không ổn định. Việc kết hợp số liệu quan trắc để hiệu chỉnh hệ số thấm $k_x, k_y$ trong module SEEP/W đã giúp nâng độ chính xác của bài toán ổn định dốc mái đập thêm 14,2% so với phương pháp tra bảng chỉ tiêu cơ lý quy ước. Điều này khẳng định vai trò sống còn của công tác hiệu chỉnh mô hình bằng số liệu thực tế trong quản lý an toàn đập.

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

Dựa trên kết quả phân tích động lực học dòng thấm, 4 nhóm giải pháp kỹ thuật và quản lý vận hành được đề xuất:

  1. Thiết lập quy trình kiểm soát tốc độ dâng nước hồ chứa: Khống chế tốc độ tích nước tối đa không vượt quá 0,45 m/ngày trong suốt chu kỳ tích nước ban đầu. Ban Quản lý Đầu tư và Xây dựng Thủy lợi 4 cùng đơn vị điều hành hồ cần phối hợp xả lũ và tích nước nhịp nhàng để duy trì hệ số an toàn mái đập luôn lớn hơn 1,35 trong thời gian 12 tháng đầu vận hành.
  2. Tự động hóa toàn diện hệ thống quan trắc thấm: Nâng cấp 100% các điểm đo thủ công sang hệ thống cảm biến thông minh kết nối SCADA truyền tín hiệu liên tục 1 giờ/lần trong mùa mưa lũ và 12 giờ/lần trong mùa khô. Kế hoạch hoàn thành lắp đặt và đồng bộ dữ liệu trong vòng 6 tháng do Công ty Quản lý Khai thác Công trình Thủy lợi chủ trì.
  3. Giám sát nghiêm ngặt gradient thấm tại các vị trí tiếp giáp: Thiết lập ngưỡng cảnh báo sớm tự động khi gradient thủy lực vượt quá 0,22 tại vùng tiếp giáp giữa khối đất đắp và màn khoan phụt vữa xi măng nền đập. Định kỳ 3 tháng một lần thực hiện đo kiểm tra lưu lượng thấm tại các rãnh tập trung hạ lưu.
  4. Cập nhật và hiệu chỉnh định kỳ mô hình số công trình: Ứng dụng phần mềm GeoStudio để chạy lại bài toán thấm và ổn định 6 tháng một lần dựa trên chuỗi số liệu quan trắc thực tế mới nhất. Chủ thể thực hiện là đơn vị tư vấn chuyên ngành phối hợp với Trường Đại học Thủy lợi nhằm đưa ra các kịch bản ứng phó trước mỗi mùa mưa bão.

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

  1. Kỹ sư thiết kế công trình thủy công: Luận văn cung cấp phương pháp luận chi tiết và quy trình tính toán dòng thấm không ổn định qua đập đất đá nhiều khối quy mô lớn, giúp tối ưu hóa kích thước hình học và bố trí lăng trụ thoát nước.
  2. Cán bộ quản lý và vận hành hồ chứa thủy lợi, thủy điện: Tài liệu là cẩm nang thực tiễn hướng dẫn quy trình theo dõi số liệu áp lực kẽ rỗng, giải đoán đường bão hòa thực tế và xây dựng quy chế tích nước an toàn cho hồ chứa.
  3. Cơ quan quản lý nhà nước về an toàn đập: Cung cấp bằng chứng khoa học và bộ tiêu chí định lượng về thấm để phục vụ công tác thẩm định, đánh giá rủi ro và ban hành các quy chuẩn kỹ thuật quốc gia về an toàn hồ đập.
  4. Giảng viên, học viên cao học và nghiên cứu sinh ngành Xây dựng Thủy lợi: Tài liệu tham khảo học thuật chuyên sâu về ứng dụng phương pháp phần tử hữu hạn FEM kết hợp kiểm chứng số liệu quan trắc hiện trường trong phân tích cơ học đất và thủy lực ngầm.

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

Tại sao phân tích thấm trong giai đoạn tích nước đầu tiên lại quan trọng hơn giai đoạn vận hành ổn định? Giai đoạn tích nước ban đầu tạo ra gradient thủy lực đột biến và áp lực nước kẽ rỗng thặng dư chưa kịp tiêu tán trong thân đập. Thống kê quốc tế chỉ ra hơn 50% sự cố vỡ đập do thấm xảy ra trong năm đầu tích nước, do đó việc phân tích giai đoạn này giúp ngăn ngừa triệt để nguy cơ sạt trượt mái và xói ngầm.

Đập chính hồ Ngàn Trươi có thông số kỹ thuật đặc thù như thế nào? Đập chính Ngàn Trươi là công trình cấp II bằng đất đá nhiều khối, chiều cao lớn nhất 53,5 m, chiều dài đỉnh 345 m, cao trình đỉnh +57,8 m. Hồ chứa có dung tích toàn bộ 775,70 triệu m³ với mực nước dâng bình thường ở cao độ +52,00 m, được xử lý chống thấm nền bằng màn khoan phụt vữa xi măng sâu.

Phần mềm GeoStudio SEEP/W giải quyết bài toán thấm không ổn định ra sao? 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 bão hòa và không bão hòa. Phần mềm cho phép tích hợp các hàm độ ẩm thể tích, quan hệ hệ số thấm biến thiên theo áp lực nước kẽ rỗng và điều kiện biên mực nước thay đổi linh hoạt theo thời gian thực.

Hệ thống quan trắc thấm tại đập Ngàn Trươi được bố trí theo nguyên tắc nào? Công trình được bố trí 4 tuyến đo áp đại diện tại các mặt cắt 3-3, 6-6, 9-9 và 12-12. Mỗi tuyến gồm các ống đo mực nước bão hòa và áp lực kế màng rung đặt so le từ thượng lưu, tim đập đến cơ hạ lưu theo tiêu chuẩn G7, đảm bảo theo dõi sát sao sự hình thành đường bão hòa.

Tốc độ dâng mực nước hồ chứa ảnh hưởng trực tiếp đến an toàn mái đập như thế nào? Tốc độ dâng nước càng nhanh thì áp lực nước kẽ rỗng bên trong khối đắp tích tụ càng lớn, làm giảm mạnh ứng suất pháp hiệu quả và lực ma sát trong của đất. Kiểm soát tốc độ dâng dưới 0,5 m/ngày giúp đất kịp cố kết, giữ hệ số an toàn ổn định mái dốc luôn đạt trên ngưỡng quy chuẩn 1,30.

Kết luận

  • Luận văn đã giải quyết toàn diện bài toán dòng thấm không ổn định trong thân đập đất đá nhiều khối giai đoạn tích nước ban đầu.
  • Xác lập thành công quy luật phân bố áp lực nước kẽ rỗng và chứng minh đường bão hòa thực tế thấp hơn từ 10% đến 15% so với tính toán lý thuyết cổ điển.
  • Làm chủ công nghệ mô hình hóa phần tử hữu hạn với SEEP/W và SLOPE/W, kết hợp hiệu chỉnh chuẩn xác bằng chuỗi dữ liệu quan trắc hiện trường 4 tuyến đo.
  • Đề xuất giới hạn tốc độ tích nước tối ưu dưới 0,45 m/ngày, đảm bảo hệ số an toàn chống trượt mái đập luôn lớn hơn 1,35 trong mọi kịch bản vận hành.
  • Định hướng lộ trình đến năm 2025: Nâng cấp hoàn chỉnh hệ thống giám sát SCADA tự động và mở rộng mô hình hóa trường thấm ba chiều cho toàn bộ cụm đầu mối.

Các cơ quan quản lý và đơn vị tư vấn thủy lợi cần áp dụng ngay phương pháp tích hợp mô hình số và số liệu quan trắc thực tế của luận văn để nâng cao năng lực quản trị rủi ro và đảm bảo an toàn tuyệt đối cho hệ thống hồ đập trên toàn quốc.