Tổng quan nghiên cứu

Đập đất là loại hình công trình dâng nước phổ biến nhất trên toàn cầu, chiếm hơn 75% tổng số các đập vật liệu địa phương được xây dựng từ năm 1963 đến nay. Tại Việt Nam, hệ thống thủy lợi ghi nhận hàng trăm hồ chứa quy mô lớn sử dụng kết cấu đập đất như đập Tả Trạch cao 55m, đập Cấm Sơn cao 40m, đập Phú Ninh cao 36m và đập hồ Núi Cốc cao 25m. Mặc dù mang lại hiệu quả kinh tế vượt trội nhờ tận dụng vật liệu địa phương, công trình đập đất luôn tiềm ẩn nguy cơ mất an toàn nghiêm trọng do hiện tượng thấm qua thân và nền đập. Điển hình trong lịch sử thủy lợi Việt Nam, sự cố vỡ đập Suối Hành năm 1986 tại tỉnh Khánh Hòa xảy ra khi mực nước hồ dâng đột ngột 4,5m (từ cao trình 26,90m lên 31,40m chỉ trong 12 giờ), làm xuất hiện các dòng thấm ngầm tại cao trình 26m và hình thành 3 hang xói ngầm đường kính 20 đến 30cm. Tương tự, sự cố vỡ đập Suối Trầu năm 1977 bắt nguồn từ lỗ rò rỉ tại cao trình 16,7m xung quanh mang cống lấy nước, tạo thành 2 hành lang xói ngầm đường kính hơn 2m cuốn trôi thân đập.

Vấn đề cốt lõi đặt ra là xử lý triệt để tầng thấm mạnh ở nền có chiều dày không quá lớn, thường dưới 6m đến 7m. Mục tiêu nghiên cứu tập trung xác định vị trí và kích thước hình học tối ưu của chân răng chống thấm (tường răng) đắp bằng chính vật liệu á sét của thân đập với hệ số thấm dao động từ 0,0001 cm/s đến 0,00001 cm/s. Nghiên cứu giới hạn phạm vi cho loại đập đất đồng chất cấp III, đồng thời ứng dụng thực nghiệm tính toán trên công trình thực tế là đập Thượng Trí vào năm 2014. Việc xác định chính xác bề rộng đáy chân khay giúp giảm từ 15% đến 25% khối lượng đào đắp đất sét so với cách tính toán truyền thống, đồng thời đảm bảo Gradient thấm tiếp xúc luôn nằm trong giới hạn an toàn từ 4 đến 6, triệt tiêu nguy cơ xói ngầm và trượt mái hạ lưu.

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 được xây dựng dựa trên 2 nền tảng lý thuyết thủy lực và cơ học đất dòng thấm kinh điển:

Định luật Darcy về dòng chảy tầng trong môi trường xốp hạt mịn: Vận tốc thấm v tỉ lệ bậc nhất với Gradient thủy lực J theo công thức v = k * J, trong đó k là hệ số thấm của vật liệu (đối với đất thân đập kd = 10^-4 đến 10^-5 cm/s; đối với đất nền cát kn lớn gấp 10 đến 100 lần kd).

Phương trình vi phân dòng thấm phẳng hai chiều ổn định Laplace biểu diễn quan hệ thế năng cột nước trong môi trường thấm dị hướng: Đạo hàm riêng d/dx (kx * dH/dx) + d/dy (ky * dH/dy) = 0, kết hợp với nguyên lý biến phân cột nước khả dĩ nhằm cực tiểu hóa phiếm hàm năng lượng dòng thấm.

Mô hình giải tích biến đổi tương đương: Biến đổi hệ đập và nền không đồng chất thành sơ đồ đồng chất giả định bằng cách quy đổi chiều dày tầng thấm nền T thành Tbd = T * (kn / kd), từ đó xác định bề rộng trung bình của lõi tượng trưng lbd. Ba khái niệm chuyên ngành then chốt bao gồm: Tỷ số thấm kn/kd, Gradient thấm cho phép tại mặt tiếp xúc đáy tường răng Jcd_cp (dao động từ 4 đến 6 đối với đất á sét), và tỷ lệ vị trí chân khay x/B tương ứng với khoảng cách từ tim đập đến chân mái thượng lưu.

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

Nguồn dữ liệu nghiên cứu được thiết lập qua hệ thống 17 kịch bản mô phỏng toàn diện, bao gồm các tổ hợp thông số biến thiên: Chiều cao đập Hd thay đổi ở các mức 18m, 20m và 25m; chiều dày tầng thấm nền T khảo sát tại 3,5m và 4,5m; tỷ số hệ số thấm giữa nền và thân đập kn/kd thay đổi từ 10 đến 100. Cỡ mẫu 17 trường hợp đại diện đầy đủ cho điều kiện địa chất công trình hồ chứa quy mô vừa và nhỏ tại Việt Nam. Phương pháp chọn mẫu phân tầng theo cấp công trình thủy lợi cấp III đảm bảo tính đại diện cao.

Lý do lựa chọn phương pháp phân tích số là vì phương pháp phần tử hữu hạn (FEM) chia miền thấm thành mạng lưới các phần tử tam giác phẳng 3 nút, cho phép giải quyết chính xác sự tập trung ứng suất dòng thấm tại các góc gãy khúc của chân khay – nơi phương pháp giải tích 1 chiều không thể mô tả được. Tác giả sử dụng module Seep/W thuộc bộ phần mềm Geo-Slope với thuật toán Green & Corey để giải phương trình dòng thấm phi tuyến. Toàn bộ chuỗi nghiên cứu, thu thập số liệu và mô phỏng số được thực hiện có hệ thống trong suốt 12 tháng tại Trường Đại học Thủy Lợi.

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

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

Phát hiện 1 về vị trí chân răng: Vị trí bố trí chân răng tối ưu nhất trong thân đập là tại tỷ lệ x/B = 0,5 (khu vực trục tim hoặc dưới mép nước thượng lưu). Khi bố trí tại vị trí này, lưu lượng thấm qua thân đập giảm hơn 30% và đường bão hòa được hạ thấp nhanh chóng. Ngược lại, nếu bố trí chân khay lệch về phía hạ lưu tại khoảng x/B = 0,7 đến 1,0, Gradient thấm Jcd tăng vọt trên 40%, gây ướt mái và làm mất ổn định nghiêm trọng mái dốc hạ lưu.

Phát hiện 2 về sai lệch giữa các phương pháp: Phương pháp biến đổi tương đương thủy lực có xu hướng ước tính bề rộng đáy chân răng l quá lớn so với yêu cầu thực tế. Khi tỷ số thấm kn/kd tăng từ 10 lên 100 lần, phương pháp thủy lực cho ra kích thước đáy chân răng rộng hơn từ 20% đến 35% so với phương pháp phần tử hữu hạn Seep/W.

Phát hiện 3 về hệ số hiệu chỉnh: Luận văn đã xác lập được hàm quan hệ hệ số hiệu chỉnh eta = f(kn/kd, T). Kỹ sư có thể xác định bề rộng đáy chân răng thực tế theo công thức l' = eta * l, trong đó eta dao động từ 0,65 đến 0,85 tùy thuộc vào chiều dày T từ 3,5m đến 4,5m và tỷ số thấm kn/kd. Điều này giúp độ chính xác của phương pháp giải tích đạt mức tương đồng với mô hình số Seep/W với sai số dưới 3%.

Phát hiện 4 về ứng dụng thực tế: Tại công trình đập Thượng Trí, kích thước bề rộng đáy chân răng được tối ưu hóa ở mức 4,5m. Gradient thấm tiếp xúc tính toán đạt Jcd = 2,85, nhỏ hơn rất nhiều so với Gradient cho phép Jcd_cp = 4,0, đảm bảo an toàn tuyệt đối chống xói ngầm.

Thảo luận kết quả

Nguyên nhân dẫn đến sự chênh lệch kích thước giữa hai phương pháp là do phương pháp giải tích giả định các đường dòng thấm chuyển động song song dạng 1 chiều, trong khi thực tế dòng thấm bị uốn cong mạnh mẽ tại khu vực tiếp giáp giữa chân răng và nền thấm mạnh. Mô hình số Seep/W mô phỏng chính xác trường vận tốc hai chiều và sự tiêu tán áp lực nước lỗ rỗng cục bộ.

Toàn bộ dữ liệu nghiên cứu có thể được trực quan hóa xuất sắc qua 23 biểu đồ tương quan giữa Gradient thấm và hệ số hiệu chỉnh eta, cùng các bảng tổng hợp phân tích độ nhạy của 17 trường hợp mô phỏng. So với các nghiên cứu kinh điển trước đây vốn bắt buộc phải mở rộng chân khay rất lớn gây tốn kém, kết quả này chứng minh rằng việc áp dụng hệ số hiệu chỉnh eta vừa thỏa mãn điều kiện bền thấm của đất theo Tiêu chuẩn Quốc gia, vừa giảm thiểu đáng kể chi phí đào móng và khối lượng vật liệu sét đắp đập.

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

Thứ nhất, áp dụng ngay biểu đồ tra cứu hệ số hiệu chỉnh eta trong công tác tính toán thủy lực thiết kế chân răng cho các dự án đập đất cấp III và cấp IV. Giải pháp này giúp các đơn vị tư vấn thiết kế rút ngắn 40% thời gian lập dự toán và giảm từ 15% đến 20% khối lượng đào đắp đất nền trong vòng 6 tháng tới.

Thứ hai, chuẩn hóa quy trình mô phỏng thấm 2 chiều bằng phần mềm Seep/W đối với tất cả các đập có cột nước cao trên 15m xây dựng trên nền bồi tích có tỷ số thấm kn/kd vượt quá 20 lần. Cục Quản lý Đê điều và Phòng chống Thiên tai phối hợp với Viện Thủy điện và Năng lượng Tái tạo ban hành sổ tay hướng dẫn kỹ thuật trong năm 2026.

Thứ ba, kiểm soát chặt chẽ chỉ tiêu cơ lý của vật liệu đắp chân răng trong quá trình thi công. Đất á sét đắp chân khay phải có chỉ số dẻo Wn lớn hơn 7 và nhỏ hơn 20, độ ẩm đầm nén tối ưu dao động quanh mức 16% đến 18%, bảo đảm dung trọng khô và hệ số đầm nén K đạt từ 0,95 đến 0,98 nhằm triệt tiêu hoàn toàn các khe nứt tiếp xúc dọc theo đáy chân răng.

Thứ tư, thiết lập hệ thống quan trắc thấm tự động với tối thiểu 5 đến 7 ống đo áp lực nước lỗ rỗng đặt tại mặt cắt đáy chân khay và vai đập. Đơn vị quản lý khai thác công trình thủy lợi cần duy trì tần suất đo 1 lần mỗi tuần trong 12 tháng đầu tích nước hồ chứa để theo dõi sát sao sự phát triển của đường bão hòa.

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

Nhóm 1 - Kỹ sư tư vấn thiết kế công trình thủy lợi: Ứng dụng ngay công thức tính toán thủy lực kết hợp hệ số hiệu chỉnh eta để xác định nhanh mặt cắt hình học và kích thước đáy chân răng, tối ưu hóa hồ sơ thiết kế kỹ thuật thi công cho các dự án hồ chứa nước.

Nhóm 2 - Nhà thầu thi công và giám sát địa kỹ thuật: Nắm bắt các yêu cầu kỹ thuật nghiêm ngặt về đào mở mái móng chân khay (hệ số mái m > 1,0 đến 1,5), công nghệ đầm nén đất á sét liên kết với tầng không thấm, ngăn ngừa hiện tượng lún nứt cục bộ.

Nhóm 3 - Cán bộ quản lý, vận hành hồ đập: Hiểu rõ cơ chế xuất hiện các dòng thấm tập trung và hang ngầm nguy hiểm (qua bài học sự cố đập Suối Hành và Suối Trầu), từ đó xây dựng phương án phòng chống lụt bão và xử lý sự cố rò rỉ nước kịp thời.

Nhóm 4 - 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 công trình tham khảo chuẩn mực về sự kết hợp giữa lý thuyết cơ học môi trường liên tục, giải pháp giải tích tương đương và mô hình hóa phần tử hữu hạn hiện đại trong kỹ thuật công trình thủy.

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

Chân răng chống thấm đắp bằng chính đất thân đập có ưu điểm gì nổi bật so với các giải pháp khác? Giải pháp này tạo nên sự đồng nhất hoàn toàn về biến dạng lún giữa thân đập và chân răng, loại bỏ nguy cơ nứt gãy tiếp xúc. Việc tận dụng trực tiếp nguồn đất á sét tại mỏ vật liệu địa phương (hệ số thấm kd từ 10^-4 đến 10^-5 cm/s) giúp tiết kiệm từ 25% đến 35% chi phí so với phương án khoan phụt vữa xi măng - bentonite hay đóng cọc cừ thép.

Vị trí bố trí chân răng ở đâu trong thân đập là an toàn và hợp lý nhất? Vị trí tối ưu được xác định tại đoạn giữa thân đập tương ứng tỷ lệ x/B = 0,5 (dưới mép nước thượng lưu). Vị trí này đảm bảo triệt tiêu lưu lượng thấm và hạ thấp đường bão hòa hiệu quả nhất. Không được bố trí chân răng tại x/B từ 0,7 đến 1,0 vì Gradient thấm hạ lưu sẽ vượt ngưỡng an toàn gây sạt lở mái đập.

Tại sao phương pháp phần tử hữu hạn Seep/W lại cho kích thước chân răng nhỏ hơn phương pháp thủy lực? Mô hình Seep/W phân tích chính xác dòng thấm 2 chiều với mạng lưới phần tử tam giác liên tục, phản ánh đúng quy luật tiêu hao áp lực nước khi đi qua khúc uốn chân răng. Phương pháp thủy lực cổ điển dựa trên các giả thiết dòng chảy 1 chiều nên thường tính toán dư thừa từ 20% đến 35% bề rộng đáy để bù trừ sai số hình học.

Khi nào nên áp dụng biện pháp chân răng cắt qua tầng thấm mạnh ở nền đập? Giải pháp chân răng đất đắp áp dụng hiệu quả nhất khi chiều dày tầng thấm mạnh ở nền T nhỏ hơn 6m đến 7m và tầng không thấm nằm ở độ sâu máy đào cơ giới có thể tiếp cận được. Nếu tầng thấm dày hơn 10m đến 12m, giải pháp màng khoan phụt bentonite hoặc cừ chống thấm mới là phương án khả thi về mặt kỹ thuật.

Hệ số hiệu chỉnh eta được áp dụng vào quy trình tính toán thực tế như thế nào? Người thiết kế chỉ cần tính toán bề rộng đáy chân răng theo công thức giải tích cổ điển, sau đó nhân với hệ số hiệu chỉnh eta được tra cứu trực tiếp từ biểu đồ quan hệ f(kn/kd, T). Kích thước thu được vừa đảm bảo Gradient thấm tiếp xúc nhỏ hơn giới hạn cho phép từ 4 đến 6, vừa đạt độ chính xác tương đương chạy mô hình số phức tạp.

Kết luận

  • Luận văn làm sáng tỏ bản chất thủy lực của dòng thấm qua đập đất đồng chất có chân răng xử lý tầng thấm mạnh ở nền.
  • Xác lập vị trí chân răng tối ưu tại trục tim đập x/B = 0,5, giúp giảm lưu lượng thấm và kiểm soát Gradient đáy luôn dưới ngưỡng tới hạn.
  • Đề xuất thành công hệ số hiệu chỉnh eta, giải quyết triệt để nhược điểm tính toán kích thước chân răng quá lớn của phương pháp giải tích truyền thống.
  • Tối ưu hóa kích thước đáy chân khay công trình đập Thượng Trí ở mức 4,5m, đảm bảo an toàn vận hành và tiết kiệm hàng trăm triệu đồng chi phí đào đắp.
  • Mở ra hướng chuẩn hóa quy chuẩn tính toán thiết kế nền đập đất cho hệ thống hồ chứa thủy lợi Việt Nam giai đoạn 2026 - 2030.

Công trình đã cung cấp một công cụ tính toán thực hành chuẩn xác, kết hợp hài hòa giữa phương pháp số hiện đại và công thức giải tích kinh điển, mang lại đóng góp to lớn cho ngành kỹ thuật xây dựng công trình thủy lợi. Trong giai đoạn 12 đến 24 tháng tới, việc nhân rộng phương pháp tra cứu hệ số eta tại các ban quản lý dự án nông nghiệp và phát triển nông thôn sẽ nâng cao đáng kể chất lượng thẩm định an toàn đập. Hãy tham khảo và ứng dụng ngay các kết quả nghiên cứu chuyên sâu của luận văn để tối ưu hóa giải pháp chống thấm cho các công trình đập đất bền vững trong tương lai.