Tổng quan nghiên cứu

Hệ thống công trình thủy lợi đóng vai trò huyết mạch trong phát triển nông nghiệp và an ninh nguồn nước quốc gia. Theo thống kê của Ủy ban Quốc tế về Đập lớn (ICOLD), đập đất chiếm tới 61,8% trong tổng số hơn 39.000 đập lớn trên thế giới nhờ ưu điểm tận dụng vật liệu địa phương và chi phí thi công tối ưu. Tại Việt Nam, mạng lưới hồ chứa có quy mô hơn 6.080 công trình với tổng dung tích thiết kế đạt 12.477 triệu m³. Tuy nhiên, qua quá trình khai thác lâu dài, cả nước hiện có khoảng 1.200 hồ chứa bị hư hỏng, xuống cấp nghiêm trọng, trong đó sự cố thấm qua thân và nền đập là nguyên nhân hàng đầu đe dọa an toàn công trình.

Hồ chứa nước An Long thuộc xã Quế Phong, huyện Quế Sơn, tỉnh Quảng Nam (tọa độ địa lý 15°36'00" vĩ độ Bắc và 108°09'00" kinh độ Đông), được xây dựng từ năm 1985 và đưa vào vận hành năm 1987. Công trình có nhiệm vụ cấp nước tưới trực tiếp cho 250 ha đất canh tác nông nghiệp. Sau hơn 30 năm vận hành, đập đất hồ An Long xuất hiện hiện tượng thấm nghiêm trọng: đường bão hòa dâng cao lộ ra mái hạ lưu, xuất hiện vùng sình lầy tại bờ tả và nguy cơ trượt sạt cục bộ mái dốc.

Mục tiêu cốt lõi của nghiên cứu là đánh giá hiện trạng mất ổn định thấm, xác định rõ nguyên nhân gây rò rỉ và lựa chọn giải pháp kỹ thuật tối ưu nhằm gia cố chống thấm cho đập. Ý nghĩa thực tiễn của công trình thể hiện qua việc triệt tiêu nguy cơ vỡ đập, khôi phục 100% dung tích trữ nước hữu ích, bảo vệ tính mạng vùng hạ du và đảm bảo nguồn nước tưới ổn định cho 250 ha lúa, hoa màu của địa phương.

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 cơ bản trong cơ học đất và thủy lực công trình ngầm, khởi đầu từ phương trình bảo toàn năng lượng Bernoulli với cột nước tổng $h = u/\gamma_w + v^2/2g + z$. Do lưu tốc dòng thấm $v$ qua môi trường lỗ rỗng của đất rất nhỏ, thành phần cột nước vận tốc $v^2/2g$ được triệt tiêu, đưa phương trình về dạng cột nước đo áp $u = h + z$.

Khung lý thuyết phân tích dòng thấm dựa trên hai trạng thái vận động chính:

  • Định luật thấm tuyến tính Darcy: Ứng dụng cho dòng thấm chảy tầng trong môi trường đất bão hòa với công thức $Q = k \cdot J \cdot F$, trong đó $k$ là hệ số thấm (m/ngày đêm), $J$ là gradient thủy lực và $F$ là diện tích mặt cắt ướt.
  • Định luật thấm phi tuyến: Áp dụng khi vận tốc dòng thấm lớn vượt qua giới hạn chảy tầng sang chảy rối, xác định qua chỉ số Reynolds $Re = 7 \div 9$ theo thực nghiệm của Nikolai Pavlovsky. Các mô hình phi tuyến của Chezy-Kraxnoponxki ($v = k_k \cdot J^{1/2}$), Smơrêke ($v = k_s \cdot J^{1/n}$ với $n = 1 \div 2$) và Proni ($J = av + bv^2$) được sử dụng để phân tích các đới thấm nứt nẻ tập trung.

Các khái niệm then chốt bao gồm: dòng thấm không áp với mặt biên tự do, đường bão hòa, gradient thấm tới hạn gây xói ngầm và hiện tượng đẩy trồi cơ học tại chân mái hạ lưu.

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

Nguồn dữ liệu của đề tài được tổng hợp từ hồ sơ hoàn công năm 1987, tài liệu quan trắc địa kỹ thuật và số liệu khoan khảo sát hiện trường tại tuyến đập hồ An Long giai đoạn 2016-2017. Cỡ mẫu nghiên cứu bao gồm hệ thống các hố khoan địa chất dọc tim đập và vai đập tả, kết hợp lấy mẫu đất nguyên dạng thí nghiệm trong phòng để xác định các chỉ tiêu cơ lý: dung trọng khô, góc ma sát trong, lực dính kết và hệ số thấm $k$ của từng lớp địa tầng.

Phương pháp chọn mẫu tuân thủ theo tiêu chuẩn khảo sát địa kỹ thuật công trình thủy lợi, bố trí các hố khoan cắt ngang qua các vị trí xung yếu có biểu hiện sình lầy và nứt nẻ mặt đập.

Về phân tích tính toán, luận văn tích hợp ba nhóm phương pháp:

  1. Phương pháp giải tích và bán thực nghiệm: Phương pháp cơ học chất lỏng Pavlovsky, phương pháp hệ số sức kháng Trugaev và phương pháp tỷ lệ đường thẳng của Bligh - Lane ($L_{tt} = L_d + L_n/m \ge C \cdot H$ với $m = 1 \div 3,5$) nhằm sơ bộ xác định chiều dài đường viền thấm an toàn.
  2. Phương pháp mô hình số phần tử hữu hạn (FEM): Ứng dụng bộ phần mềm chuyên ngành GeoStudio 2007 (module SEEP/W để mô phỏng trường thấm và module SLOPE/W để tính toán hệ số ổn định trượt mái dốc). Lý do lựa chọn mô hình số FEM vì phương pháp này mô tả chính xác tính phi đồng nhất, bất đẳng hướng của các lớp đất đắp và điều kiện biên phức tạp khi có màn chống thấm.

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

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

Nghiên cứu đã làm rõ thực trạng kỹ thuật và lượng hóa các chỉ số an toàn của đập đất hồ chứa nước An Long qua các phát hiện cụ thể:

Thứ nhất, kết quả khảo sát địa kỹ thuật xác định nguyên nhân thấm xuất phát từ chất lượng thi công đắp đất đồng chất giai đoạn 1985-1987 chưa đạt độ chặt thiết kế tại một số phân đoạn, kết hợp với hiện tượng thấm tiếp xúc mạnh tại vị trí khớp nối giữa vai đập bờ tả với sườn đồi tự nhiên. Hệ số thấm cục bộ của thân đập đo được lên tới $k = 4,5 \times 10^{-4}$ cm/s, vượt ngưỡng cho phép đối với đập đất đồng chất.

Thứ hai, mô phỏng SEEP/W ở trạng thái hiện trạng khi chưa xử lý cho thấy đường bão hòa dâng rất cao, thoát ra tại cao trình 1/3 mái hạ lưu, tạo ra vùng sình lầy kéo dài và làm gradient thấm lối ra đạt $J_{ra} = 0,42$, vượt quá gradient giới hạn cho phép ($J_{cp} \approx 0,20 \div 0,25$). Điều này gây xói ngầm cơ học cuốn trôi các hạt mịn trong thân đập.

Thứ ba, tính toán ổn định trên SLOPE/W chỉ ra hệ số an toàn trượt mái hạ lưu trong điều kiện mực nước dâng bình thường chỉ đạt $K_{min} = 1,12$, không đáp ứng tiêu chuẩn an toàn công trình thủy lợi (yêu cầu $K \ge 1,35 \div 1,50$).

Thứ tư, khi áp dụng giải pháp khoan phụt màng chống thấm xi măng - sét sâu vào tầng đá gốc không thấm nước, lưu lượng thấm tổng cộng qua thân và nền đập giảm từ $Q = 14,8 \times 10^{-3}$ m³/s xuống còn $Q = 2,6 \times 10^{-3}$ m³/s (giảm hơn 82% lượng nước thất thoát). Đồng thời, đường bão hòa bị kéo tụt xuống dưới lăng trụ tiêu nước, triệt tiêu hoàn toàn hiện tượng sũng ướt mái hạ lưu.

Thảo luận kết quả

Nguyên nhân sự cố tại hồ An Long mang tính điển hình cho các công trình đập vật liệu địa phương xây dựng trong giai đoạn thập niên 1980 tại miền Trung: hạn chế về thiết bị đầm nện cơ giới, công tác xử lý tiếp giáp vai đập bằng thủ công và thiếu màng ngăn thấm chuyên dụng. So sánh với các sự cố vỡ đập nổi tiếng như đập Teton (Hoa Kỳ, cao 93 m vỡ năm 1976 do thấm qua khe nứt đá nền) hay đập Suối Trầu (Khánh Hòa vỡ năm 1977 do thấm xói ngầm quanh cống lấy nước), sự cố tại hồ An Long nếu không được can thiệp kịp thời sẽ dẫn đến kịch bản tạo hang rỗng xói ngầm tập trung và gây vỡ đập cục bộ.

Về mặt giải pháp, việc lựa chọn công nghệ khoan phụt vữa sét - xi măng tạo màng chống thấm liên tục có ưu thế vượt trội so với tường hào bentonite hay màng chống thấm HDPE. Giải pháp này không đòi hỏi phải hạ thấp mực nước hồ, không phá vỡ kết cấu thân đập hiện hữu và xử lý triệt để được cả tầng đá nứt nẻ dưới đáy đập.

Toàn bộ kết quả tính toán thấm và ổn định có thể được trực quan hóa thông qua bảng so sánh hệ số an toàn $K$ và biểu đồ phân bố đường đẳng áp lực trước và sau xử lý:

Kịch bản tính toán Đường bão hòa lộ mái Gradient thấm lối ra ($J_{ra}$) Lưu lượng thấm $Q$ ($10^{-3}$ m³/s) Hệ số ổn định $K_{min}$ Đánh giá an toàn
Trước khi xử lý Có (thoát ra ở mái) 0,42 14,8 1,12 Không an toàn
Sau khi khoan phụt màng Không (nằm sâu trong đập) 0,11 2,6 1,48 Đạt tiêu chuẩn

Dữ liệu mô phỏng SEEP/W cho thấy việc tạo lập màng phụt có chiều dày tương đương và hệ số thấm $k \le 1 \times 10^{-6}$ cm/s đã triệt tiêu hoàn toàn áp lực thấm cục bộ, đưa trạng thái làm việc của thân đập về vùng ổn định lâu dài.

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

Dựa trên kết quả nghiên cứu lý thuyết và mô phỏng thực tế, 4 nhóm giải pháp kỹ thuật và quản lý được đề xuất triển khai:

  1. Thi công màng chống thấm bằng công nghệ khoan phụt vữa xi măng - sét: Đơn vị quản lý dự án cần chỉ đạo nhà thầu chuyên ngành tổ chức khoan phụt một hàng hố khoan dọc theo tim đập với khoảng cách tim hố $a = 1,5 \div 2,0$ m, ngàm sâu vào tầng đá ít thấm tối thiểu 1,0 m. Mục tiêu đưa hệ số thấm của thân đập đạt $k \le 10^{-5} \div 10^{-6}$ cm/s, hoàn thành trong thời gian 6 đến 9 tháng trước mùa mưa bão.
  2. Cải tạo và khôi phục lăng trụ tiêu nước hạ lưu: Công ty TNHH MTV Khai thác Thủy lợi Quảng Nam tiến hành nạo vét, thay thế tầng lọc ngược bị bồi tắc bằng các lớp cát sỏi chọn lọc, đảm bảo khả năng thu và thoát nước ngấm với gradient dòng thoát $J \le 0,15$, thực hiện trong vòng 90 ngày sau khi hoàn thiện khoan phụt.
  3. Lắp đặt hệ thống quan trắc thấm tự động: Bố trí 5 đến 7 cụm ống đo áp (piezometer) tại mặt cắt đại diện vai tả và lòng sông để theo dõi áp lực nước lỗ rỗng và cao trình đường bão hòa định kỳ 7 ngày/lần trong mùa mưa lũ.
  4. Chuẩn hóa quy trình vận hành và ứng phó thiên tai: Chi cục Thủy lợi và Ban Chỉ huy Phòng chống thiên tai địa phương xây dựng kịch bản vận hành xả lũ an toàn khi lượng mưa vượt quá 200 mm/24h, đồng thời tổ chức kiểm tra thực địa định kỳ 6 tháng/lần nhằm phát hiện sớm các dị thường trên thân đập.

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

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

  1. Kỹ sư tư vấn thiết kế và địa kỹ thuật công trình thủy lợi: Ứng dụng quy trình tính toán mô hình số SEEP/W và SLOPE/W để thiết kế màng khoan phụt chống thấm cho các hồ chứa vừa và nhỏ có dung tích từ 1 đến 10 triệu m³.
  2. Đơn vị quản lý và khai thác công trình thủy lợi: Vận dụng phương pháp khảo sát, xác định nguyên nhân thấm qua vai đập để lập phương án bảo trì định kỳ cho các hồ đập có tuổi thọ trên 30 năm.
  3. Học viên cao học, nghiên cứu sinh chuyên ngành Địa kỹ thuật - Thủy công: Sử dụng luận văn như tài liệu học thuật tham khảo về phương pháp kết hợp giữa lý thuyết thấm phi tuyến, giải tích Pavlovsky và phương pháp phần tử hữu hạn trong giải quyết bài toán thấm thực tiễn.
  4. Cơ quan quản lý nhà nước về an toàn đập: Tham khảo các tiêu chí kỹ thuật và quy trình kiểm định để xây dựng cẩm nang gia cố cho hơn 1.200 hồ chứa đang trong tình trạng xuống cấp trên phạm vi toàn quốc.

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

Hiện tượng đường bão hòa lộ ra mái hạ lưu đập đất nguy hiểm như thế nào? Khi đường bão hòa thoát ra mặt mái dốc hạ lưu, áp lực nước lỗ rỗng làm giảm ứng suất hữu hiệu và sức chống cắt của đất, gây sình lầy và trượt mái dốc. Nếu gradient dòng thấm tại điểm thoát vượt quá ngưỡng cho phép ($J \ge 0,25$), hiện tượng xói ngầm cơ học sẽ cuốn trôi hạt đất, tạo hàm ếch và gây sập đổ thân đập trong thời gian ngắn.

Tại sao phương pháp khoan phụt vữa xi măng - sét được chọn thay vì tường hào bentonite? Khoan phụt vữa xi măng - sét phù hợp tuyệt đối với điều kiện địa hình chật hẹp tại vai đập hồ An Long và xử lý triệt để được tầng đá gốc nứt nẻ mà không cần đào phá thân đập. Ngược lại, tường hào bentonite đòi hỏi mặt bằng thi công rộng, máy móc cồng kềnh và chi phí huy động thiết bị cao hơn khoảng 30% đến 40% đối với công trình quy mô nhỏ.

Module SEEP/W trong bộ phần mềm GeoStudio giải quyết bài toán thấm ra sao? SEEP/W sử dụng phương pháp phần tử hữu hạn để giải phương trình vi phân đạo hàm riêng của dòng thấm 2 chiều trong môi trường đất bão hòa và không bão hòa. Phần mềm cho phép mô phỏng chính xác trường thế năng, vẽ đường bão hòa, xác định vận tốc dòng thấm và tính toán lưu lượng thấm qua từng lớp vật liệu có hệ số thấm $k$ khác nhau.

Gradient thấm lối ra bao nhiêu thì đảm bảo an toàn không bị xói ngầm? Theo tiêu chuẩn thiết kế đập đất TCVN, gradient thấm lối ra an toàn ($J_{ra}$) phải nhỏ hơn gradient thấm cho phép ($J_{cp} = J_{th}/K_s$). Đối với đất đắp pha sét tại hồ An Long, gradient giới hạn phá hoại $J_{th} \approx 0,50 \div 0,60$ và hệ số an toàn $K_s = 2,0$, do đó giá trị $J_{ra}$ sau khi xử lý phải đạt dưới 0,20 để ngăn ngừa tuyệt đối hiện tượng xói ngầm.

Sau khi xử lý chống thấm, hồ chứa nước An Long có đảm bảo cấp nước tưới lâu dài không? Sau khi hoàn thiện màng phụt chống thấm, lưu lượng rò rỉ giảm trên 82%, giúp bảo toàn dung tích trữ nước hữu ích của hồ trong suốt mùa khô. Công trình đảm bảo cung cấp nước tưới ổn định, liên tục cho toàn bộ 250 ha diện tích canh tác nông nghiệp của xã Quế Phong theo đúng nhiệm vụ thiết kế ban đầu.

Kết luận

  • Luận văn đã phân tích toàn diện thực trạng 30 năm khai thác hồ An Long, chỉ rõ nguyên nhân thấm do chất lượng đầm nén cục bộ và thấm tiếp xúc vai đập bờ tả.
  • Ứng dụng thành công lý thuyết thấm và phương pháp phần tử hữu hạn (GeoStudio SEEP/W và SLOPE/W) để lượng hóa chính xác trường thấm và nguy cơ mất ổn định trượt mái dốc.
  • Đề xuất giải pháp khoan phụt màng chống thấm vữa xi măng - sét ngàm vào tầng đá gốc, giúp giảm 82% lưu lượng thấm và nâng hệ số an toàn ổn định mái đập từ $K = 1,12$ lên $K = 1,48$.
  • Xây dựng quy trình công nghệ thi công, mạng lưới bố trí hố khoan và tiêu chí kiểm tra chất lượng màng phụt đạt tiêu chuẩn ngành.
  • Cung cấp mô hình tính toán chuẩn xác, có tính ứng dụng cao cho việc sửa chữa, nâng cấp mạng lưới 1.200 hồ chứa thủy lợi đang xuống cấp tại Việt Nam.

Đóng góp lớn nhất của công trình là kết hợp chặt chẽ giữa cơ sở lý thuyết địa kỹ thuật hiện đại với bài toán thực tế của địa phương, đưa ra phương án kỹ thuật khả thi và tiết kiệm chi phí đầu tư. Các đơn vị quản lý cần tiếp tục thực hiện quan trắc thấm định kỳ trong 12 tháng đầu sau thi công để đánh giá hiệu quả vận hành thực tế. Hãy chủ động liên hệ các chuyên gia địa kỹ thuật thủy lợi để nhận tài liệu chi tiết và chuyển giao quy trình tính toán cho công trình của bạn.